具有弱夹层的封闭油藏中的瞬态渗流

多层油藏分为层间无窜流的油藏和层间有窜流的油藏。本文研究了层间有弱渗透夹层的两层油藏中两种情形的瞬态渗流,包括单层打开的瞬态渗流和一层注入、一层采出的两层油藏的瞬态渗流。用面积平均方法求出了弱渗透夹层封闭油藏中瞬态渗流的平均压降。分析表明:在渗流前期,层间窜流量逐渐增大,在后期层间窜流量趋于稳定;在单层开采的情况下:初期两层压降不同,非打开层的压降有滞后现象,并且隔层渗透率越小,油井工作强度越弱,滞后时间越长;后期两层压降趋于一致,可以作为一层油藏处理。在一注一采的油藏中,当注采量相同时,两层压降在后期趋于不同的稳定值;当注采量不同时,两层压降在后期趋于一致,但不能达到稳定。     

孔隙介质中稠油流体非线性渗流方程

为揭示稠油流体在油藏孔隙中渗流特性,基于力学平衡方程,建立了描述稠油流体渗流特征的非线性渗流方程,对油藏孔隙中稠油渗流过程及启动机理进行了深入分析,着重分析了边界层、流体屈服应力以及表面力对渗流过程的影响.结果表明,Hagen-Poiseuille定律需经修正方能描述稠油流体流动,边界层、流体屈服应力以及表面力对稠油渗流影响非常显著.孔隙中,稠油启动压力梯度来源于其屈服应力、表面力,边界层加剧了渗流非线性程度,实际稠油油藏开发中,要充分掌握稠油渗流非线性特征.      

基于渗透率连续变化的低渗透多孔介质非线性渗流模型

基于低渗透多孔介质渗透率的渐变理论,确定了能精确描述低渗透多孔介质渗流特征的非线性运动方程,并通过实验数据拟合．验证了非线性运动方程的有效性。非线性渗流速度关于压力梯度具有连续-阶导数,方便于工程计算;由此建立了低渗透多孔介质的单相非线性径向渗流数学模型,并巧妙采用高效的Douglas-Jones预估一校正有限差分方法求得了其数值解。数值结果分析表明：非线性渗流模型为介于拟线性渗流模型和达西渗流模型之间的一种中间模型或理想模型,非线性渗流模型和拟线性渗流模型均存在动边界;拟线性渗流高估了启动压力梯度的影响,使得动边界的移动速度比实际情况慢得多;非线性越强,地层压力下降的范围越小,地层压力梯度越陡峭,影响地层压力的敏感性减弱,而影响地层压力梯度的敏感性增强。     

破碎岩体非等温渗流的非线性动力学研究

分别从固体及流体导热的能量方程出发,导出破碎岩体非等温渗流的能量本构方程, 结合渗流的连续性方程、运动方程、状态方程等建立了破碎岩体非等温渗流的一维非线性动力学方程组;结合Mathcad软件计算得到了系统的无量纲化平衡态, 利用逐次亚松弛迭代法分析了对应于不同参数时平衡态的稳定性;指出非等温渗流系统存在鞍结分岔及折叠突变, 与等温渗流相比, 考虑温度场的破碎岩体渗流动力系统更容易发生渗流突变.      

弹塑性变形油藏中多相渗流的数值模拟

基于流固耦合力学理论,建立了弹性变形油藏中多相渗流的数学模型,该模型考虑了渗流与变形的耦合作用,以及注采交变载荷作下油藏多孔介质的弹性变形特征,给出了耦合数值模拟方法和算例。     

变形三重介质低渗透油藏垂直裂缝井的非线性流动分析

不仅考虑低渗透油藏具有启动压力梯度的渗流特征,还考虑应力敏感地层中介质的变形;发展了Cinco-Ley H.提出的有限导流垂直裂缝井双线性流理论,将流体在垂直裂缝与地层中形成的流动划分为两个区域—垂直裂缝中的线性流区域和变形三重介质低渗透油藏中的非线性流区域;由此建立了变形三重介质低渗透油藏有限导流垂直裂缝井的非线性流...     

复杂边界非均质渗流场流线分布研究

建立了考虑源(汇)影响的含有不渗透区域复杂边界条件下非均质油藏稳定渗流的数学模型。利用扰动边界元方法求解数学模型,获得了地层中任意一点的压力公式.在此基础上,提出了流线场的生成方法。绘制了考虑非均质性、复杂边界和不渗透区域影响的流线分布图,并分析了流线分布的特征。通过分析表明,渗流场的非均质性和不渗透区域的存在都对流线分布存在较大的影响。利用本文方法产生的流线分布图能够较为直观地反映出油藏流体在注采井间的运动轨迹,为优化井网和注入方案提供了重要依据。     

轴对称二维定常和不定常渗流

本文研究碳酸盐岩地层条件下单相液体轴对称二维定常和不定常渗流。给出了注、采平衡时底水块状油藏单井极限产量的计算公式,借助该式可预估最大无水产油量。找到了双重介质中弹性渗流的线源解,它是不完整井试井资料分析、计算的理论基础。     

峰后岩石非Darcy渗流的分岔行为研究

煤矿采动围岩大多处于峰后应力状态或破碎状态，其渗流一般不符合Darcy定律，为非Darcy渗流系统．峰后岩石非Darcy渗流系统的失稳和分岔是煤矿突水和煤与瓦斯突出动力灾害发生的根源．文中用谱截断方法建立了Ahmed-Sunada型非Darcy渗流系统的降阶动力学方程，再由变量代换得到以无量纲变量表示的平衡态附近的演化方程，分析了系统的分岔条件，给出了系统的各种吸引子图案，并结合采矿工程实际，用非线性数学的观点揭示了煤矿突水和煤与瓦斯突出的机理．研究表明：当非Darcy渗流系统渗流特性和边界压力的初始值满足一定条件时，系统由平衡转向不稳定，即存在跨临界Hopf分岔和切分岔，并且，系统的动力学响应不随渗透特性连续变化，即该系统存在突变性．     

渗流力学的近况、动向和展望(续)

六、非均质地层中的渗流 地层的宏观非均质性及其对渗流的影响越来越为人们所认识和重视。不可压缩均质液体在非均质地层中的渗流归结为求解变系数椭圆型方程;当液体为可压缩时,则归结为求解变系数抛物型方程;一般情况下无法求得解析解,多年来只获得了地层参数按特定规律分布时的一些解答。在非均质液体渗流情形下则是更加难以求解的非线性数学物理问题。迄今对非均质地层的渗流问题只能求其近似解。大型高速电子数字计算机使这类问题的数值求解变得越来越有效和方便。     

考虑二次梯度项及动边界的双重介质低渗透油藏流动分析

在传统试井模型的非线性偏微分方程中根据弱可压缩流体的假设,忽略了二次梯度项,对于低渗透油藏这种方法是有疑问的.低渗透问题一个显著的特点就是流体的流动边界随着时间不断向外扩展.为了更好地研究双重介质低渗透油藏中流体的流动问题,考虑了二次梯度项及活动边界的影响,同时考虑了低渗透油藏的非达西渗流特征,建立了双重介质低渗透油藏流动模型.采用Douglas-Jones预估-校正差分方法获得了无限大地层定产量生产时模型的数值解,分别讨论了不同参数变化时压力的变化规律及活动边界随时间的传播规律,还分析了考虑和忽略二次梯度项影响时模型数值解之间的差异随时间的变化规律,做出了典型压力曲线图版,这些结果可用于实际试井分析.     

Oil Recovery by Low Salinity Water Injection into a Reservoir: A New Study of Tertiary Oil Recovery Mechanism

Low salinity water injections for oil recovery have shown seemingly promising results in the case of clay-bearing sandstones saturated with asphaltic crude oil. Reported data showed that low salinity water injection could provide up to 20% pore volume (PV) of additional oil recovery for core samples and up to 25% PV for reservoirs in near wellbore regions, compared with brine injection at the same Darcy velocity. The question remains as to whether this additional recovery is also attainable in reservoirs. The answer requires a thorough understanding of oil recovery mechanism of low salinity water injections. Numerous hypotheses have been proposed to explain the increased oil recovery using low salinity water, including migration of detached mixed-wet clay particles with absorbed residual oil drops, wettability alteration toward increased water-wetness, and emulsion formation. However, many later reports showed that a higher oil recovery associated with low salinity water injection at the common laboratory flow velocity was neither necessarily accompanied by migration of clay particles, nor necessarily accompanied by emulsion. Moreover, increased water-wetness has been shown to cause the reduction of oil recovery. The present study is based on both experimental and theoretical analyses. Our study reveals that the increased oil recovery is only related to the reduction of water permeability due to physical plugging of the porous network by swelling clay aggregates or migrating clay particles and crystals. At a fixed apparent flow velocity, the value of negative pressure gradient along the flow path increases as the water permeability decreases. Some oil drops and blobs can be mobilized under the increased negative pressure gradient and contribute to the additional oil recovery. Based on the revealed mechanism, we conclude that low salinity water injection cannot be superior to brine injection in any clay-bearing sandstone reservoir at the maximum permitted injection pressure. Through our study of low salinity water injection, the theory of tertiary oil recovery has been notably improved.     

Experimental displacement patterns in a 2 × 2 quadrant block with permeability and wettability heterogeneities—problems for numerical modelling

The effect of heterogeneities on miscible and immiscible flood displacements in 2D bead packs in quadrant form, 2 × 2 block heterogeneity, with either a permeability or a wettability contrast is the subject of this paper. The physical processes occurring during miscible and immiscible flow and displacement within permeability and wettability quadrant bead pack models have been studied experimentally. This geometry occurs in a number of situations relevant to hydrocarbon production: particularly faults where adjacent rocks have large permeability contrasts with rapid changes, in the laboratory with core butting, in reservoir simulation where grid blocks have different permeability and in reservoirs having near-wellbore damage problems. The model quadrants 1–4, had 1 and 4 and 2 and 3 with identical properties, either in permeability or wettability. Reported are complete unit mobility miscible displacements, then the effects of viscosity differences (mobility modifiers) and finally immiscible displacements on displacement patterns for initial linear injection. The experiments demonstrate that nodal flow occurs for both miscible and immiscible flow, but for immiscible flow there are boundary effects due to capillary pressure differences created by water saturation changes or wettability contrasts which can leave patches of isolated fluid within a quadrant. The displacement patterns for the different models and fluids change significantly with the viscosity and wettability changes, particularly for the immiscible displacements. This is due to the changing capillary pressure between the quadrant blocks as the water saturation change. These are difficult to address in numerical modelling but should be accounted for. Other effects include coupling of all physical processes governing the flow through the node and creations of microzones of trapped residual oil. Our displacement patterns can therefore be a valuable verification benchmark tool for numerical modelling and a calibration data source for those wishing to simulate the effects of capillary pressure under differing wettability conditions and for those investigating upscaling modelling procedures. However, the possible loss of physical reality when averaging must always be considered.     

变形双重介质广义流动分析

对于碳酸盐油藏和低渗油藏的渗流问题，传统的研究方法都是假设地层渗透率是常数，这假设，对于地层渗透率是压力敏感的情况，对压力的空间变化和瞬时变化将导致较大的误差。本文研究了应力敏感地层中双重介质渗流问题的压力不稳定响应，不仅考虑了储层的双重介质特征，而且考虑了应力敏感地层中介质的变形，建立了应力敏感地层双重介质的数学模型，渗透率依赖于孔隙压力变化的流动方程是强非线性的，采用Douglas-Jones预估－校正法获得了只有裂缝发生形变定产量生产时无限大地层的数值解及定产量生产岩块与裂隙同时发生形变时无限大地层的数值解，并探讨了变形参数和双重介质参数变化时压力的变化规律，给出几种情况下典型压力曲线图版，这些结果可用于实际试井分析。     

A Numerical Model of Tracer Transport in a Non-isothermal Two-Phase Flow System for {\text{ CO}}_2 Geological Storage Characterization

For the purpose of characterizing geologically stored $\text{ CO}_{2}Air sparging is an in situ soil/groundwater remediation technology, which involves the injection of pressurized air through air sparging well below the zone of contamination. To investigate the rate-dependent flow properties during multistep air sparging, a rule-based dynamic two-phase flow model was developed and applied to a 3D pore network which is employed to characterize the void structure of porous media. The simulated dynamic two-phase flow at the pore scale or microscale was translated into functional relationships at the continuum-scale of capillary pressure?Csaturation (P ^c?CS) and relative permeability??saturation (K _r?CS) relationships. A significant contribution from the air injection pressure step and duration time of each air injection pressure on both of the above relationships was observed during the multistep air sparging tests. It is observed from the simulation that at a given matric potential, larger amount of water is retained during transient flow than that during steady flow. Shorter the duration of each air injection pressure step, there is higher fraction of retained water. The relative air/water permeability values are also greatly affected by the pressure step. With large air injection pressure step, the air/water relative permeability is much higher than that with a smaller air injection pressure step at the same water saturation level. However, the impact of pressure step on relative permeability is not consistent for flows with different capillary numbers (N _ca). When compared with relative air permeability, relative water permeability has a higher scatter. It was further observed that the dynamic effects on the relative permeability curve are more apparent for networks with larger pore sizes than that with smaller pore sizes. In addition, the effect of pore size on relative water permeability is higher than that on relative air permeability.     

Two-phase flow for a horizontal well penetrating a naturally fractured reservoir with edge water injection

This paper examines the two-phase flow for a horizontal well penetrating a naturally fractured reservoir with edge water injection by means of a fixed streamline model. The mathematical model of the vertical two-dimensional flow or oil-water for a horizontal well in a medium with double-porosity is established, and whose accurate solutions are obtained by using the characteristic method. The saturation distributions in the fractured system and the matrix system as well as the formula of the time of water free production are presented. All these results provide a theoretical basis and a computing method for oil displacement by edge water from naturally fractured reservoirs.     

低渗透多孔介质中的非线性渗流理论简

文中论述了低渗透性多孔介质中非线性渗流理论的几个问题,阐明了渗流流体的性质,指出了多孔介质对流体通过的选择性,提出了新的非线性渗流方程,用实验资料对其进行了验证,分析了该方程演变功能,表明它可以描述各种渗流规律.该方程的各项参数都可从实验中直接得到,应用方便,并且参数的物理意义明确.     

Numerical modeling of the process of oil displacement from reservoirs with zonal nonuniformity

In extracting oil from nonuniform reservoirs a considerable fraction remains unrecovered from the zones of lesser permeability. The mechanism of displacement of oil from reservoirs with zonal nonuniformity is investigated within the framework of the two-dimensional Muskat-Meres model of combined oil, water and gas flow [1]. A wholly conservative difference scheme implicit in the saturations and pressure is used for the calculations. Various reservoir exploitation regimes are considered with the object of seeking means of improving the characteristics of the process.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 177–180, September–October, 1987.     

Modeling the effects of capillary pressure with the presence of full tensor permeability and discrete fracture models using the mimetic finite difference method

Capillary dominated flow or imbibition—whether spontaneous or forced—is an important physical phenomena in understanding the behavior of naturally fractured water-driven reservoirs (NFR’s). When the water flows through the fractures, it imbibes into the matrix and pushes the oil out of the pores due to the difference in the capillary pressure. In this paper, we focus on modeling and quantifying the oil recovered from NFR’s through the imbibition processes using a novel fully implicit mimetic finite difference (MFD) approach coupled with discrete fracture/discrete matrix (DFDM) technique. The investigation is carried out in the light of different wetting states of the porous media (i.e., varying capillary pressure curves) and a full tensor representation of the permeability. The produced results proved the MFD to be robust in preserving the physics of the problem, and accurately mapping the flow path in the investigated domains. The wetting state of the rock affects greatly the oil recovery factors along with the orientation of the fractures and the principal direction of the permeability tensor. We can conclude that our novel MFD method can handle the fluid flow problems in discrete-fractured reservoirs. Future works will be focused on the extension of MFD method to more complex multi-physics simulations.