裂缝性低渗透油藏流固耦合渗流分析

在低渗透油田的开发过程中，油藏流体渗流和储层岩土之间存在明显的耦合作用。本文首先研究给出了低渗裂缝性储层孔渗参数的等效方法，然后将渗流力学和岩土力学相结合，给出了低渗透裂缝性储层流固耦合渗流的数学模型，该模型不仅可以反映基质孔渗参数在开发中的变化，而且更能反映裂缝开度变化所引起的渗透率变化，而这对于低渗透裂缝性油田而言十分重要。最后对一实际井网进行了流固耦合油藏数值模拟，给出了开发过程中孔渗参数的变化及其耦合效应对油田开发的影响．     

Numerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured Reservoirs

A numerical method as well as a theoretical study of non-Darcy fluid flow through porous and fractured reservoirs is described. The non-Darcy behavior is handled in a three-dimensional, multiphase flow reservoir simulator, while the model formulation incorporates the Forchheimer equation for describing single-phase or multiphase non-Darcy flow and displacement. The non-Darcy flow through a fractured reservoir is handled using a general dual-continuum approach. The numerical scheme has been verified by comparing its results against those of analytical methods. Numerical solutions are used to obtain some insight into the physics of non-Darcy flow and displacement in reservoirs. In addition, several type curves are provided for well-test analyses of non-Darcy flow to demonstrate a methodology for modeling this type of flow in porous and fractured rocks, including flow in petroleum and geothermal reservoirs.     

The Experimental and Numerical Studies on Gas Production from Hydrate Reservoir by Depressurization

A set of experimental system to study hydrate dissociation in porous media is built and some experiments on hydrate dissociation by depressurization are carried out. A mathematical model is developed to simulate the hydrate dissociation by depressurization in hydrate-bearing porous media. The model can be used to analyze the effects of the flow of multiphase fluids, the kinetic process and endothermic process of hydrate dissociation, ice-water phase equilibrium, the variation of permeability, convection and conduction on the hydrate dissociation, and gas and water productions. The numerical results agree well with the experimental results, which validate our mathematical model. For a 3-D hydrate reservoir of Class 3, the evolutions of pressure, temperature, and saturations are elucidated and the effects of some main parameters on gas and water rates are analyzed. Numerical results show that gas can be produced effectively from hydrate reservoir in the first stage of depressurization. Then, methods such as thermal stimulation or inhibitor injection should be considered due to the energy deficiency of formation energy. The numerical results for 3-D hydrate reservoir of Class 1 show that the overlying gas hydrate zone can apparently enhance gas rate and prolong life span of gas reservoir.     

Permeability Changes During the Evolution of a Geothermal Field Due to the Dissolution and Precipitation of Quartz

The changes in permeability and porosity associated with quartz deposition in an evolving geothermal reservoir are investigated. We review the processes associated with permeability changes in a geothermal reservoir and also review recent work on quartz solubility and deposition rates. Porosity and permeability changes are calculated for two reservoir models. The first is a generic model of a Taupo Volcanic Zone geothermal reservoir and the second is based on the model of the deep circulation system at Kakkonda published by Hanano. We find that when a reservoir experiences BPD conditions during part of its lifetime the lower reservoir becomes impermeable and a hydrothermal circulation system is only present in the upper part of the reservoir.     

3D Coupled Displacement Discontinuity and Finite Element Analysis of Reservoir Behavior during Production in Semi-infinite Domain

In this paper, we apply a DDM-FEM coupling method to reservoir simulation to evaluate the reservoir behavior over a compacting oil reservoir in half space. We use displacement discontinuity methods to account for the reservoir surroundings, and finite element methods in the fully coupled simulation of the reservoir itself.     

利用改进的ADI方法实现单机千万节点的油水两相渗流数值模拟

针对实际油藏的非均质分布特征及其复杂的边界条件,本文通过引入迭代参数的压缩因子和放大因子,对现有的交替方向迭代法(ADI)进行改进,提出一种适用于大规模油藏数值模拟的新算法．改进的ADI 方法计算精度可靠,且与现有的算法相比计算效率有所提高．更为关键的是,ADI 算法将求解三维压力方程的七对角矩阵分解为三个方向的三对角矩阵依次迭代求解,所需的运算存储量大幅减少,最大的计算规模有了大幅的提升．使用改进的ADI 方法,在单机上成功实现了千万节点的油水两相渗流数值模拟．计算实例表明,在同等单机硬件条件下,改进ADI 算法的最大计算规模是现有算法的1.7 倍以上．     

Effect of Gas Injection from the Surface of a Body on Its Interaction with a Temperature Inhomogeneity in a Supersonic Flow

The flow pattern in the shock layer and the aerodynamic characteristics of a hemisphere in unsteady axisymmetric interaction with a closed spherical hot-gas region embedded in the oncoming supersonic flow in the presence of intense injection of gas from the body surface into the shock layer are studied on the basis of the inviscid perfect gas model. Two cases are considered, namely, (1) when the radius of the permeable surface is greater than that of the temperature inhomogeneity and (2) when the injection is localized in the vicinity of the forward stagnation point and the permeable region is smaller the inhomogeneity.     

Pore-Level Modeling of Gas and Condensate Flow in Two- and Three-Dimensional Pore Networks: Pore Size Distribution Effects on the Relative Permeability of Gas and Condensate

We present a mechanistic model of retrograde condensation processes in two- and three-dimensional capillary tube networks under gravitational forces. Condensate filling-emptying cycles in pore segments and gas connection–isolation cycles are included. With the pore-level distribution of gas and condensate in hand, we determine their corresponding relative permeabilities. Details of pore space and displacement are subsumed in pore conductances. Solving for the pressure field in each phase, we find a single effective conductance for each phase as a function of condensate saturation. Along with the effective conductance for the saturated network, the relative permeability for each phase is calculated. Our model porous media are two- and three-dimensional regular networks of pore segments with distributed size and square cross-section. With a Monte Carlo sampling we find the optimum network size to avoid size effects and then we investigate the effect of network dimensionality and pore size distribution on the relative permeabilities of gas and condensate.     

Positive Effect of Flow Velocity on Gas–Condensate Relative Permeability: Network Modelling and Comparison with Experimental Results

Positive velocity dependency of relative permeability of gas–condensate systems, which has been observed in many different core experiments, is now well acknowledged. The above behaviour, which is due to two-phase flow coupling in condensing systems at low interfacial tension (IFT) conditions, was simulated using a 3D pore network model. The steady-dynamic bond network model developed for this purpose was also equipped with a novel anchoring technique, which was based on the equivalent hydraulic length concept adopted from fluid flow through pipes. The available rock data on the co-ordination number, capillary pressure, absolute permeability, porosity and one set of measured relative permeability curves were utilised to anchor the capillary, volumetric and flow characteristics of the constructed network model to those properties of the real core sample. Then the model was used to predict the effective permeability values at other IFT and velocity levels. There is a reasonable quantitative agreement between the predicted and measured relative permeability values affected by the coupling rate effect.     

可聚合添加剂和极压添加剂对矿物油极压抗磨和抗疲劳性能的影响

选用重极压齿轮油和二聚酸 /司本 - 80添加剂制备了多种不同配比的油样 ,进行了疲劳和胶合试验 ,同时分析了磨损表面形貌和抗磨机理 .结果表明 :聚合添加剂的耐疲劳性能较好 ,其与极压添加剂经合理复配可以更好地提高油样的耐疲劳及抗胶合综合性能 .这是两类添加剂协同作用的结果     

Numerical Studies on the Mixing of CH4 and Kerosene Injected into a Supersonic Flow with H2 Pilot Injection

Two-fluid model and divisional computation techniques were used. The multi-species gas fully N-S equations were solved by upwind TVD scheme. Liquid phase equations were solved by NND scheme. The phases-interaction ODE equations were solved by 2nd Runge-Kutta approach. The favorable agreement is obtained between computational results and PLIF experimental results of iodized air injected into a supersonic flow. Then, the numerical studies were carried out on the mixing of CH ₄ and kerosene injected into a supersonic flow with H ₂ pilot injection. The results indicate that the penetration of kerosene approaches maximum when it is injected from the second injector. But the kerosene is less diffused compared with the gas fuels. The free droplet region appears in the flow field. The mixing mechanism of CH ₄ with H ₂ pilot injection is different from that of kerosene. In the staged duct, H ₂ can be entrained into both recirculation zones produced by the step and injectors. But CH ₄ can only be carried into the recirculation between the injectors. Therefore, initiations of H ₂ and CH ₄ can occur in those regions. The staged duct is better in enhancing mixing and initiation with H ₂ pilot flame.     

The Effect of Partial Premixing and Heat Loss on the Reacting Flow Field Prediction of a Swirl Stabilized Gas Turbine Model Combustor

This work addresses the prediction of the reacting flow field in a swirl stabilized gas turbine model combustor using large-eddy simulation. The modeling of the combustion chemistry is based on laminar premixed flamelets and the effect of turbulence-chemistry interaction is considered by a presumed shape probability density function. The prediction capabilities of the presented combustion model for perfectly premixed and partially premixed conditions are demonstrated. The effect of partial premixing for the prediction of the reacting flow field is assessed by comparison of a perfectly premixed and partially premixed simulation. Even though significant mixture fraction fluctuations are observed, only small impact of the non-perfect premixing is found on the flow field and flame dynamics. Subsequently, the effect of heat loss to the walls is assessed assuming perfectly premixing. The adiabatic baseline case is compared to heat loss simulations with adiabatic and non-adiabatic chemistry tabulation. The results highlight the importance of considering the effect of heat loss on the chemical kinetics for an accurate prediction of the flow features. Both heat loss simulations significantly improve the temperature prediction, but the non-adiabatic chemistry tabulation is required to accurately capture the chemical composition in the reacting layers.