油藏复杂驱动体系物理模拟相似准则研究进展

综述了水驱油、化学驱油、蒸汽驱、混相驱及非混相驱等复杂油藏驱替体系物理模拟相似准则的研究现状.系统介绍了相似参数敏感性分析方法,定义了表征目标函数对相似参数依赖程度的敏感因子,结合数值方法便可定量地确定复杂体系主要相似参数,并实际应用于工程问题.讨论了水驱和聚合物驱物理模拟应优先满足的相似准则及其随参数范围的变化.      

Theoretical analysis of microscopic oil displacement mechanism by viscoelastic polymer solution

The microscopic oil displacement mechanism in viscoelastic polymer flooding is theoretically analyzed with mechanical method. The effects of viscoelasticity of polymer solution on such three kinds of residual oil as in pore throat, in sudden expansion pore path, and in dead end are analyzed. Results show that the critical radius of mobile residual oil for viscoelastic polymer solution is larger than that for viscous polymer solution, which makes the oil that is immobile in viscous polymer flooding displaced under the condition of viscoelastic polymer solution. The viscous polymer solution hardly displaces the oil in dead ends. However, when the effect of viscoelasticity is considered, the residual oil in sudden expansion pore paths and dead ends can be partly displaced. A dimensionless parameter is suggested to denote the relative dominance of gravity and capillary pressure. The larger the dimensionless parameter, the more accurate the increment expressions.     

Effects of Physical Parameter Range on Dimensionless Variable Sensitivity in Water Flooding Reservoirs

The similarity criterion for water flooding reservoir flows is concerned with in the present paper. When finding out all the dimensionless variables governing this kind of flow, their physical meanings are subsequently elucidated. Then, a numerical approach of sensitivity analysis is adopted to quantify their corresponding dominance degree among the similarity parameters. In this way, we may finally identify major scaling law in different parameter range and demonstrate the respective effects of viscosity, permeability and injection rate.The project supported by the Innovative Project of CAS (KJCX-SW-L08) and the National Basic Research Program of China(973)     

Designing of Partial Similarity Models and Evaluation Method in Polymer Flooding Experiment

Based on the scaling criteria of polymer flooding reservoir obtained in our previous work in which the gravity and capillary forces, compressibility, non-Newtonian behavior, absorption, dispersion, and diffusion are considered, eight partial similarity models are designed. A new numerical approach of sensitivity analysis is suggested to quantify the dominance degree of relaxed dimensionless parameters for partial similarity model. The sensitivity factor quantifying the dominance degree of relaxed dimensionless parameter is defined. By solving the dimensionless governing equations including all dimensionless parameters, the sensitivity factor of each relaxed dimensionless parameter is calculated for each partial similarity model; thus, the dominance degree of the relaxed one is quantitatively determined. Based on the sensitivity analysis, the effect coefficient of partial similarity model is defined as the summation of product of sensitivity factor of relaxed dimensionless parameter and its relative relaxation quantity. The effect coefficient is used as a criterion to evaluate each partial similarity model. Then the partial similarity model with the smallest effect coefficient can be singled out to approximate to the prototype. Results show that the precision of partial similarity model is not only determined by the number of satisfied dimensionless parameters but also the relative relaxation quantity of the relaxed ones.     

??????????????????????????????

采用相似准数敏感性分析方法,对各局部相似 模型中所放松相似准数的敏感性进行分析,得到应尽量减少放松的相似准数. 提出了定量评 价局部相似模型优劣的影响系数,通过比较该影响系数的大小即可定量评价各局部相似模型. 结果表明：水驱局部相似模型的精度不仅取决于所满足的相似准数的个数,还与相似准数的 相对放松量有关.     

Analytical and Experimental Study to Predict the Residual Resistance Factor on Polymer Flooding Process in Fractured Medium

The major objectives of this study are to analytically and experimentally determine the residual resistance factor in the fractured medium based on the polymer solution properties and operational conditions. The parameters considered in this study are the polymer concentration, power law constitutive equation parameter, and salt concentration, sulfonation content of polymer, temperature, and molecular weight of the water soluble polymers which are used in polymer flooding for enhanced oil recovery. The results indicated that residual resistance factor in fractured medium is dependent on the coil overlap parameter and power law equation parameter of polymer. The coil overlap parameter is a dimensionless number consists of intrinsic viscosity and polymer concentration. Since intrinsic viscosity is a function of polymer diameter in medium conditions, to predict the residual resistance factor in fracture medium, an experimental correlation is generated for determination of the molecular diameter of polymer based on polymer molecular weight, temperature, salt concentration, and sulfonation content.     

Experimental Investigation of Synergy of Components in Surfactant/Polymer Flooding Using Three-Dimensional Core Model

Surfactant/polymer (SP) floods have significant potentials to recover remaining oil after water flooding. Their efficiency can be maximized by fully utilizing synergistic effect of polymer and surfactant. Various components adsorbed on the rock matrix due to chromatographic separation can significantly weaken the synergistic effect. Due to scale and dimensional problems, it is hard to investigate chromatographic separation among various components using one-dimensional natural cores. This study compared the adsorption difference between artificial and natural cores and developed a three-dimensional artificial core model of a 1/4 5-spot configuration to simulate oil recovery in multilayered reservoirs with high, middle and low permeability for each layer. Sampling wells were established to monitor pressures, and effluent fluids were acquired to measure interfacial tension (IFT) and viscosity. Then, distances of synergy of polymer and surfactant in three layers were evaluated. Meanwhile, electrodes were set in the model to measure oil saturation variation with resistance changes at different locations. Through comparison with IFT values, the contribution of improved swept volume and oil displacement efficiency to oil recovery during SP flooding could be known. Results showed that injected 0.65 PV of SP could improve oil recovery by 21.56% when water cut reached 95% after water flooding. The retention ratio of polymer viscosity was kept 55.3% at the outlet, but IFT was only 2 mN/m within the 3/10 injector–producer spacing during SP injection. Although subsequent water flooding could result in surfactant desorption and the IFT became 10^?2?mN/m within the 3/10 injector–producer spacing, the IFT turned to 2?mN/m at the half of the model. The enhanced displacement efficiency by reducing IFT only worked within three-tenth location of the model in the high permeability layer, while the enlarged swept volume contributed much in the other areas.

     

Investigating the Temperature Dependency of Oil and Water Relative Permeabilities for Heavy Oil Systems

A look into the literature on the temperature dependency of oil and water relative permeabilities reveals contradictory reports. There are some publications reporting shifts in the water saturation range as well as variations in the relative permeability curves by temperature. On the other hand, some authors have blamed the experimental artifacts, viscous instabilities and fingering issues for these variations. We have performed core flooding experiments to further investigate this issue. Glass bead packs and sand packs were used as the porous media, and Athabasca bitumen with varying viscosities was displaced by hot water at differing temperatures. The unsteady-state method of relative permeability measurement was applied and the experimental data were history matched by a simulator that is tailor made to predict the relative permeabilities. The matches were obtained by varying the relative permeability correlation parameters. The results indicated that the initial water saturation has a direct relation with temperature, while residual oil saturation generally drops at higher temperatures. Although the water saturation range shifts, no direct and unique trend for either oil or water relative permeability is justified. The spread in relative permeabilities especially in the case of higher permeable cores suggests that viscous instabilities are present. As the same saturation shift happens by only changing the oil viscosity, the relative permeability variations with temperature can be attributed to oil to water viscosity ratio changes with temperature. Temperature dependency of relative permeabilities is more related to experimental artifacts, viscous fingering and viscosity changes than fundamental flow properties.     

Unstable Immiscible Displacements in Oil-Wet Rocks

Displacement of a viscous fluid by a lower viscosity immiscible fluid (such as waterflood of a viscous oil) in a porous medium is unstable. The displacement front generates viscous fingers which lead to low oil recovery efficiency. These fingers are much smaller in width than typical reservoir simulation grid blocks, and capturing their effect in reservoir simulation is important. A dimensionless scaling group (viscous finger number) had been suggested in the past, which has a power-law relationship with the breakthrough recovery and cumulative recovery in unstable core floods. The relative permeability used in large grid block simulations had been modified to so-called pseudo-relative permeability on the basis of the dimensionless group, thus incorporating the effect of fingers in waterflood predictions. However, the previous proposed models were constructed from experiments in only water-wet rocks. This paper extends the recent viscous fingering models to oil-wet systems. Sandstone cores were treated to alter the wettability to oil-wet. Adverse viscosity water floods were performed in oil-wet cores. Viscosity ratio, velocity and diameter were varied. It is shown that the previously developed viscous finger number does not work for the oil-wet experiments. The correlating dimensionless number is modified for oil-wet systems; it is also different from the dimensionless group identified by Peters and Flock (Soc Petroleum Eng, 1981. doi: 10.2118/8371-PA) for oil-wet cores. A pseudo-relative permeability model has been developed for oil-wet cores. Corefloods have been matched by the new pseudo-relative permeability model to determine the model parameters. This pseudo-relative permeability model can be used in reservoir simulations of water and polymer floods in viscous oil-wet reservoirs.     

Pore-Scale Monitoring of Wettability Alteration by Silica Nanoparticles During Polymer Flooding to Heavy Oil in a Five-Spot Glass Micromodel

It is well known that the oil recovery is affected by wettability of porous medium; however, the role of nanoparticles on wettability alteration of medium surfaces has remained a topic of debate in the literature. Furthermore, there is a little information of the way dispersed silica nanoparticles affect the oil recovery efficiency during polymer flooding, especially, when heavy oil is used. In this study, a series of injection experiments were performed in a five-spot glass micromodel after saturation with the heavy oil. Polyacrylamide solution and dispersed silica nanoparticles in polyacrylamide (DSNP) solution were used as injected fluids. The oil recovery as well as fluid distribution in the pores and throats was measured with analysis of continuously provided pictures during the experiments. Sessile drop method was used for measuring the contact angles of the glass surface at different states of wettability after coating by heavy oil, distilled water, dispersed silica nanoparticles in water (DSNW), polyacrylamide solution, and DSNP solution. The results showed that the silica nanoparticles caused enhanced oil recovery during polymer flooding by a factor of 10%. The distribution of DSNP solution during flooding tests in pores and throats showed strong water-wetting of the medium after flooding with this solution. The results of sessile drop experiments showed that coating with heavy oil, could make an oil-wet surface. Coating with distilled water and polymer solution could partially alter the wettability of surface to water-wet and coating with DSNW and DSNP could make a strongly water-wet surface.     

非牛顿流体在非均质油藏渗流压力场实验

在非均质油藏模型上进行非牛顿流体流动物理模拟实验，对比研究水驱、聚合物驱和交联聚合物对提高石油采收率的影响．通过布置高精度的压差传感器测量不同驱替过程模型中的渗流压力场的动态变化，成胶后的交联聚合物封堵了高渗条区，改变了油藏内流体流动方向，驱替出低渗区内油，提高了采收率．     

Mathematical model of two-phase fluid nonlinear flow in low-permeability porous media with applications

IntroductionItisasuccessfulexampleinadevelopmentstoryofscienceandtechnologyformechanicsoffluidsinporousmediatocombinewithengineeringtechnology .Fieldsinfluencedbythemechanicsinvolveddevelopmentofoil_gasandgroundwaterresources,controlonseawaterintrusionandsubsidenceandgeologichazards,geotechnicalengineeringandbioengineering ,andairlineindustry[1～ 7].Aproblemonnonlinearflowinlow_permeabilityporousmediaisbutonlyabasiconeindifferentkindsofengineeringfields,butalsooneoffrontlineresearchfieldsofmod…     

A Mathematical Model for Hysteretic Two-Phase Flow in Porous Media

We develop a mathematical model for hysteretic two-phase flow (of oil and water) in waterwet porous media. To account for relative permeability hysteresis, an irreversible trapping-coalescence process is described. According to this process, oil ganglia are created (during imbibition) and released (during drainage) at different rates, leading to history-dependent saturations of trapped and connected oil. As a result, the relative permeability to oil, modelled as a unique function of the connected oil saturation, is subject to saturation history. A saturation history is reflected by history parameters, that is by both the saturation state (of connected and trapped oil) at the most recent flow reversal and the most recent water saturation at which the flow was a primary drainage. Disregarding capillary diffusion, the flow is described by a hyperbolic equation with the connected oil saturation as unknown. This equation contains functional relationships which depend on the flow mode (drainage or imbibition) and the history parameters. The solution consists of continuous waves (expansion waves and constant states), shock waves (possibly connecting different modes) and stationary discontinuities (connecting different saturation histories). The entropy condition for travelling waves is generalized to include admissible shock waves which coincide with flow reversals. It turns out that saturation history generally has a strong influence on both the type and the speed of the waves from which the solution is constructed.     

High Resolution Imaging of Unstable,Forced Imbibition in Berea Sandstone

An experimental investigation is presented of immiscible, high-mobility ratio forced imbibition in a representative linear homogeneous sandstone. Water floods with mobility ratios (from 1 to 155) at various water injection rates were conducted. Fine-scale (order mm³) in situ water saturation history was collected via X-ray computed tomography (CT). Three-dimensional images were constructed of stable displacement and the initiation and growth of unstable water fingers. Interestingly, viscous fingers do not lead the displacement front by significant distances, counter to experience in miscible systems. In this homogeneous porous medium, both water (displacing phase) injection rate and oil (displaced phase) viscosity have an obvious effect on the stability of the water front. As the oil viscosity and displacement rate increase, the water front becomes less stable. In addition, the so-called shock mobility ratio, as computed from steady-state relative permeability, is found to be predictive regarding displacement front stability. When the shock mobility ratio is greater than 1, the displacement is always unstable. Steady-state relative permeability, however, is found to be a function of viscosity ratio for unstable displacements.     

Study of Microscopic and Macroscopic Displacement Behaviors of Polymer Solution in Water-Wet and Oil-Wet Media

Performance of a polymer flood process requires the knowledge of rheological behavior of the polymer solution and reservoir properties such as rock wettability. To provide a better understanding of effects of polymer chemistry and wettability on the performance of a polymer flood process, a comprehensive experimental study was conducted using a two-dimensional glass micromodel. A series of water and polymer flood processes were carried out at different polymer molecular weights, degrees of polymer hydrolysis, and polymer concentrations in both water-wet and oil-wet systems. Image processing technique was applied to analyze and compare microscopic and macroscopic displacement behaviors of polymer solution in each experiment. From micro-scale observations, the configuration of connate water film, polymer solution trapping, flow of continuous and discontinuous strings of polymer solution, piston-type displacement of oil, snap-off of polymer solution, distorted flow of polymer solution, emulsion formation, and microscopic pore-to-pore sweep of oil phase were observed and analyzed in the strongly oil-wet and water-wet media. Rheological experiments showed that a higher polymer molecular weight, degree of hydrolysis, and concentration result in a higher apparent viscosity for polymer solution and lower oil–polymer viscosity ratio. It is also shown that these parameters have different impacts on the oil recovery in different wettabilities. Moreover, a water-wet medium generally had higher recovery in contrast with an oil-wet medium. This experimental study illustrates the successful application of glass micromodel techniques for studying enhanced oil recovery (EOR) processes in five-spot pattern and provides a useful reference for understanding the displacement behaviors in a typical polymer flood process.     

Mathematical simulation of micellar-polymer displacement of oil from flooded strata

The micellar-polymer method of increasing the oil recovery from strata [1] is currently regarded as promising. The method consists of injecting into an oil stratum, which has previously undergone ordinary flooding, a relatively small amount, a slug, of micellar solution (5–10% of of the pore volume), which is propelled through the stratum by slugs of a highly viscous buffer fluid (aqueous solution of a polymer). In turn, the system of slugs is propelled from the injection points to the extraction wells by the water used for ordinary flooding. The displacement of the oil that remains after flooding in the stratum is achieved by a decrease in the coefficient of surface tension at the boundaries of the micellar solution with the oil and the water to the value 10^–2-10^–3 dyn/cm, which leads to a decrease in the amount of fixed oil and also to a control of the mobility of the fluids, which is achieved by varying the concentrations of the components of the micellar solution and the buffer fluid. The main components of micellar solutions are: a hydrocarbon fluid (oil or its fractions), water, surface-active substances. The relationships between the main components, and also the addition of salts and alcohol to the water component have a strong influence on the interaction between the solution and the stratal oil and water [2]. The micellar solution considered in the present paper dissolves oil but does not mix with water; the relationships between the components in it are characteristic of the solutions used to increase oil recovery from strata.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 84–93, November–December, 1982.     

SPONTANEOUS CAPILLARY IMBIBITION MODEL IN TIGHT OIL RESERVOIRS 1)

部分致密油井压后关井一段时间，压裂液返排率普遍低于30%，但是致密油气井产量反而越高，这与压裂液毛细管力渗吸排驱原油有关。然而，致密油储层致密，物性差，渗流机理复杂，尚没有形成统一的自发渗吸模型。本文基于油水两相非活塞式渗流理论，建立了压后闷井期间压裂液在毛细管力作用下自发渗吸进入致密油储层的数学模型，采用数值差分方法进行求解，并分析了相关影响因素。结果显示渗吸体积、渗吸前缘移动距离与渗吸时间的平方根呈线性正相关关系，与经典Handy渗吸理论模型预测结果一致，说明毛细管力自发渗吸模型可靠性较高。数值计算结果表明毛细管水相扩散系数是致密储层自发渗吸速率的主控参数，毛细管水相扩散系数越高，自发渗吸速率越大。毛细管水相扩散系数随着含水饱和度先增加后减小；随着束缚水饱和度、油相和水相端点相对渗透率增加而增加；随着相渗特征指数、油水黏度比和残余油饱和度增加而减小。该研究有助于深入认识致密油储层压裂液渗吸机理，对优化返排制度、提高致密油井产量具有重要意义。     

Pore-scale Simulation of Water Alternate Gas Injection

We use a three-dimensional mixed-wet random network model representing Berea sandstone to compute displacement paths and relative permeabilities for water alternating gas (WAG) flooding. First we reproduce cycles of water and gas injection observed in previously published experimental studies. We predict the measured oil, water and gas relative permeabilities accurately. We discuss the hysteresis trends in the water and gas relative permeabilities and compare the behavior of water-wet and oil-wet media. We interpret the results in terms of pore-scale displacements. In water-wet media the water relative permeability is lower during water injection in the presence of gas due to an increase in oil/water capillary pressure that causes a decrease in wetting layer conductance. The gas relative permeability is higher for displacement cycles after first gas injection at high gas saturation due to cooperative pore filling, but lower at low saturation due to trapping. In oil-wet media, the water relative permeability remains low until water-filled elements span the system at which point the relative permeability increases rapidly. The gas relative permeability is lower in the presence of water than oil because it is no longer the most non-wetting phase.