Analytical Modeling for Gravity Segregation in Gas Improved Oil Recovery of Tilted Reservoirs

Stone’s model for gravity segregation in gas improved oil recovery (IOR) indicates the distance that injected gas and water travel together before the segregation being completed (length of complete segregation). This model is very useful for co-injection of water and gas into horizontal depleted reservoirs. A proof by Rossen and van Duijn showed that Stone’s model applies to steady-state gas–liquid flow, and also foam flow, in horizontal reservoirs as long as the standard assumptions of fractional flow theory (incompressible flow, Newtonian mobilities, local equilibrium) are applied. However, until now, there has been no analytical study on the length of segregation when co-injection of water and gas occurs in tilted reservoirs. In this article, in order to extent the validity of Stone’s model to tilted reservoirs, governing equations of fluids displacement based on fractional flow theory are solved by the method of characteristics, MOC. The results are then compared to Stone’s model and to the results of a three-dimensional finite-difference compositional reservoir simulator. This study shows that Stone’s model should be corrected for tilted reservoirs and that the presented math proof can model gravity segregation in gas IOR of tilted reservoirs, appropriately. The effect of co-injecting of water and gas into tilted reservoirs on recovery efficiency is also examined.     

A Numerical Study of Microscale Flow Behavior in Tight Gas and Shale Gas Reservoir Systems

Various attempts have been made to model flow in shale gas systems. However, there is currently little consensus regarding the impact of molecular and Knudsen diffusion on flow behavior over time in such systems. Direct measurement or model-based estimation of matrix permeability for these “ultra-tight” reservoirs has proven unreliable. The composition of gas produced from tight gas and shale gas reservoirs varies with time for a variety of reasons. The cause of flowing gas compositional change typically cited is selective desorption of gases from the surface of the kerogen in the case of shale. However, other drivers for gas fractionation are important when pore throat dimensions are small enough. Pore throat diameters on the order of molecular mean free path lengths will create non-Darcy flow conditions, where permeability becomes a strong function of pressure. At the low permeabilities found in shale gas systems, the dusty-gas model for flow should be used, which couples diffusion to advective flow. In this study we implement the dusty-gas model into a fluid flow modeling tool based on the TOUGH+ family of codes. We examine the effects of Knudsen diffusion on gas composition in ultra-tight rock. We show that for very small average pore throat diameters, lighter gases are preferentially produced at concentrations significantly higher than in situ conditions. Furthermore, we illustrate a methodology which uses measurements of gas composition to more uniquely determine the permeability of tight reservoirs. We also describe how gas composition measurement could be used to identify flow boundaries in these reservoir systems. We discuss how new measurement techniques and data collection practices should be implemented in order to take advantage of this method. Our contributions include a new, fit-for-purpose numerical model based on the TOUGH+ code capable of characterizing transport effects including permeability adjustment and diffusion in micro- and nano-scale porous media.     

压裂页岩气藏多尺度耦合流动数值模拟研究

在碳达峰的国策背景之下,页岩气成为传统能源向绿色清洁低碳能源转型的重要过渡和能源支点.压后页岩气藏流体流动力学成为高效开发页岩气的关键力学问题.文章将小尺度低导流天然裂缝等效升级为连续介质,建立有机质-无机质-天然裂缝三重连续介质模型,同时对大尺度高导流裂缝采用离散裂缝模型刻画,嵌入天然裂缝连续介质中,构建多重连续/离散裂缝模型.综合考虑吸附气的非平衡非线性解吸附和表面扩散,自由气的黏性流和克努森扩散,给出页岩气在多尺度复杂介质中的非线性耦合流动数学模型.提出多尺度扩展有限单元法对离散裂缝进行显式求解,创新性构建三类加强形函数捕捉离散裂缝的局部流场特征,解决了压后页岩海量裂缝及多尺度流动通道的流动模拟难题.文章提出的模型和方法既能准确刻画高导流裂缝对渗流的影响,又克服了海量多尺度离散裂缝导致计算量增大的问题.通过算例展示了压后页岩各连续介质的压力衰减规律,发现裂缝中自由气、有机质中自由气、无机质中吸附气依次滞后的压力(浓度)扩散现象,重点分析了吸附气表面扩散系数、自由气克努森扩散系数、天然裂缝连续介质渗透率和吸附气解吸附速率对页岩气产量的影响.文章重点解决压后页岩多尺度流动通道的表征和...     

Reservoir storage and containment of greenhouse gases

This paper deals with the problem of disposal of industrial waste greenhouse gases (CO₂) into deep reservoirs. The simulator TOUGH2 was used to model the injection of 100 kg/s of CO₂ for 10 years into an aquifer 3 km deep with the object of evaluating the long-term storage prospects for this gas. Depending on the permeability structure above the injection point, some gas may escape to the surface. In the most favourable case, all of the gas dissolves into the water, and the resulting dense fluid settles in the aquifer over several thousand years. Consequently, underground storage provides a promising sink for reducing CO₂ emissions to the atmosphere.     

The Dual-Reciprocity Boundary Element Analysis for Hydraulically Fractured Shale Gas Reservoirs Considering Diffusion and Sorption Kinetics

This work presents a new application of boundary element method (BEM) to model fluid transport in unconventional shale gas reservoirs with discrete hydraulic fractures considering diffusion, sorption kinetics and sorbed-phase surface diffusion. The fluid transport model consists of two governing partial differential equations (PDEs) written in terms of effective diffusivities for free and sorbed gases, respectively. Boundary integral formulations are analytically derived using the fundamental solution of the Laplace equation for the governing PDEs and Green’s second identity. The domain integrals arising due to the time-dependent function and nonlinear terms are transformed into boundary integrals employing the dual-reciprocity method. This transformation retains the domain-integral-free, boundary-integral-only character of standard BEM approaches. In the proposed solution, the free- and sorbed-gas flow in the shale matrix is solved simultaneously after coupling the fracture flow equation of free gas. Well production performance under the effect of relaxation phenomenon due to delayed responses of sorbed gas under nonequilibrium sorption condition is rigorously captured by imposing the zero-flux condition at fracture–matrix interface for the sorbed-gas transport equation. The validity of proposed solution is verified using several case studies through comparison against a commercial finite-element numerical simulator.

     

页岩气和致密砂岩气藏微裂缝气体传输特性

页岩气和致密砂岩气藏发育微裂缝,其开度多在纳米级和微米级尺度且变化大,因此微裂缝气体传输机理异常复杂.本文基于滑脱流动和努森扩散模型,分别以分子之间碰撞频率和分子与壁面碰撞频率占总碰撞频率的比值作为滑脱流动和努森扩散的权重系数,耦合这两种传输机理,建立了微裂缝气体传输模型. 该模型考虑微裂缝形状和尺度对气体传输的影响. 模型可靠性用分子模拟数据验证.结果表明:(1)模型能够合理描述微裂缝中所有气体传输机理,包括连续流动,滑脱流动和过渡流动;(2)模型能够描述不同开发阶段,微裂缝中各气体传输机理对传输贡献的逐渐变化过程;(3)微裂缝形状和尺度影响气体传输,相同开度且宽度越大的微裂缝,气体传输能力越强,且在高压和微裂缝大开度的情况下表现更明显.      

A Mathematical Model for Determining Oil Migration Characteristics in Low-Permeability Porous Media Based on Fractal Theory

Shale reservoirs are characterized by very low permeability in the scale of nano-Darcy. This is due to the nanometer scale of pores and throats in shale reservoirs, which causes a difference in flow behavior from conventional reservoirs. Slip flow is considered to be one of the main flow regimes affecting the flow behavior in shale gas reservoirs and has been widely studied in the literature. However, the important mechanism of gas desorption or adsorption that happens in shale reservoirs has not been investigated thoroughly in the literature. This paper aims to study slip flow together with gas desorption in shale gas reservoirs using pore network modeling. To do so, the compressible Stokes equation with proper boundary conditions was applied to model gas flow in a pore network that properly represents the pore size distribution of typical shale reservoirs. A pore network model was created using the digitized image of a thin section of a Berea sandstone and scaled down to represent the pore size range of shale reservoirs. Based on the size of pores in the network and the pore pressure applied, the Knudsen number which controls the flow regimes was within the slip flow regime range. Compressible Stokes equation with proper boundary conditions at the pore’s walls was applied to model the gas flow. The desorption mechanism was also included through a boundary condition by deriving a velocity term using Langmuir-type isotherm. It was observed that when the slip flow was activated together with desorption in the model, their contributions were not summative. That, is the slippage effect limited the desorption mechanism through a reduction of pressure drop. Eagle Ford and Barnett shale samples were investigated in this study when the measured adsorption isotherm data from the literature were used. Barnett sample showed larger contribution of gas desorption toward gas recovery as compared to Eagle Ford sample. This paper has produced a pore network model to further understand the gas desorption and the slip flow effects in recovery of shale gas reservoirs.

     

Understanding Oil and Gas Flow Mechanisms in Shale Reservoirs Using SLD–PR Transport Model

More and more attention has been paid to the oil and gas flow mechanisms in shale reservoirs. The solid–fluid interaction becomes significant when the pores are in the nanoscale. The interaction changes the fluid’s physical properties and leads to different flow mechanisms in shale reservoirs from those in conventional reservoirs. By using a Simplified Local Density–Peng Robinson transport model, we consider the density and viscosity profiles, which result from solid–fluid interaction. Gas rarefaction effect is negligible at high pressure, so we assume it is viscous flow. Considering the density- and viscosity-changing effects, we proposed a slit permeability model. The velocity profiles are obtained by this newly established model. This proposed model is validated by matching the density profile and velocity profile from molecular dynamic simulation. Then, the effects of pressure and pore size on gas and oil flow mechanisms are also studied in this work. The results show that both gas and oil exhibit enhanced flow rates in nanopores. Gas-phase flow in nanopores is dominated by the density-changing effect (adsorption), while the oil-phase flow is mainly controlled by the viscosity-changing effect. Both gas and oil permeability quickly decrease to the Darcy permeability when the slit aperture becomes large. The results reported in this work are representative and should significantly help us understand the mechanisms of oil and gas flow in shale reservoirs.     

天然气驱长岩心室内实验研究

低渗透油藏注水开发效果差、采收率低,而采用气驱技术是动用此类难采储量的有效方法之一。本文利用长岩心实验模型,进行了物理模拟研究,得到了该油藏在纯气驱、纯水驱、完全水驱后气水交替驱、原始状态下气水交替驱和油藏目前注水倍数下气水交替驱等方式下的采收率和压力等变化情况,为油藏选择合理的开采方式提供了依据,并且为进一步的数值模拟工作提供了基础数据。     

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

利用ADINA软件针对某工程实例隧道建立了有限元模型,然后对其开挖过程 中地表沉降进行了预测并分析了隧道及冻结壁的应力与位移场,其结果与工程实测值相符; 在此模型基础上,对隧道水平局部冻结法施工引起地表变形的相关因素进行了模拟分析,结 果表明:冻结法施工时,增加冻结壁厚度对地表沉降影响并不大,而采用扩大冻结壁的范围或 减小隧道半径及埋深的方法可有效地控制地表沉降,这对今后的冻结法施工设计有借鉴意义.     

Pressure Transient Analysis for Multi-stage Fractured Horizontal Wells in Shale Gas Reservoirs

This article presents the PTA on the multi-stage fractured horizontal well in shale gas reservoirs incorporating desorption and diffusive flow in the matrix. Currently, most PTA models are simply based on Darcy flow both in natural fractures and matrix without considering the mechanisms of desorption and diffusion in shale matrix. Source function and Laplace transform with the numerical discrete method are employed to solve the mathematical model. The solution is presented in the Laplace domain so that the wellbore storage effect and skin factor can be easily incorporated by convolution. Type curves are plotted with Stehfest algorithm and different flow regimes are identified. The presented model could be used to interpret pressure signals more accurately for shale gas reservoirs.     

A SURFACE SUBSIDENCE MODEL OF COAL SEAM MINING IN MOUNTAIN AREA BASED ON PLATE BENDING DEFORMATION1)

分析了长壁工作面采煤上覆岩层破坏、移动的 “三带” 理论,认为地表的下沉变形受控于弯曲下沉带,且地表的最大下沉量远小于弯曲带岩层的厚度,弯曲带岩层的变形符合板的小挠度弯曲问题。在此基础上,根据山区地形地貌的特征及山体坡向与煤层倾向的关系,将地表山体简化为顺向坡、反向坡、沟槽、山岭四种情况,用薄板理论中的半逆解法分别建立了四种地貌情形下的地表下沉盆地模型。模型虽然作了一些简化假设,但体现了煤层倾角、埋深、岩层物理力学性质、地貌形状等影响开采沉陷的主要因素。     

Combining TOUGH2 and FLAC3D to Solve Problems in Underground Gas Storage

This paper describes modeling studies assessing the feasibility of increasing the maximum storage pressure in several underground natural gas storage reservoirs. This required an assessment of the potential for gas transport in the caprock and the geomechanical response to pressure change in the storage reservoir. To solve this problem in an efficient manner, two-phase flow (TOUGH2) and geomechanical (FLAC3D) models were combined in series. The TOUGH2 model was calibrated to fit pressure data collected on-site, from both the reservoir and caprock. The mechanical response of the caprock to increased storage pressure was modeled using FLAC3D, allowing assessment of the induced stresses in formations surrounding the reservoirs. We focused on two sites. In the first, field data were obtained from a deep borehole above the gas reservoir, which provided indirect observations of the geomechanical response of the caprock to pressure changes in the reservoir. In the second, open boreholes intersecting two thin caprock units immediately above the reservoir allowed gas flow to a shallower unit, significantly impacting the modeled fracture gradient.     

页岩气井组分比例变化规律研究

页岩气开发过程中,生产井产出气的组分比例会随时间发生变化.本文基于组分模型数值模拟研究了生产井中甲烷组分比例变化的规律.研究表明,吸附气、渗透率与孔隙度影响页岩气组分比例的瞬态响应特征. 吸附气显著影响组分比例的变化规律,吸附量的大小决定组分比例的变化值及组分比例导数曲线的上下位置. 渗透率影响组分比例初期变化规律,但在后期,不同渗透率对瞬态组分比例规律的影响基本一致.孔隙度对组分比例变化及其导数曲线的影响与吸附气的影响类似,但在生产初期,孔隙度对组分比例的影响要小于吸附气对组分比例的影响. 本文的研究提供了一种进行页岩地层参数评价的新方法.      

Laboratory Investigation of LNAPL Migration in Double-Porosity Soil Under Fractured Condition Using Digital Image Analysis

Gas production from shale gas reservoirs plays a significant role in satisfying increasing energy demands. Compared with conventional sandstone and carbonate reservoirs, shale gas reservoirs are characterized by extremely low porosity, ultra-low permeability and high clay content. Slip flow, diffusion, adsorption and desorption are the primary gas transport processes in shale matrix, while Darcy flow is restricted to fractures. Understanding methane diffusion and adsorption, and gas flow and equilibrium in the low-permeability matrix of shale is crucial for shale formation evaluation and for predicting gas production. Modeling of diffusion in low-permeability shale rocks requires use of the Dusty gas model (DGM) rather than Fick’s law. The DGM is incorporated in the TOUGH2 module EOS7C-ECBM, a modified version of EOS7C that simulates multicomponent gas mixture transport in porous media. Also included in EOS7C-ECBM is the extended Langmuir model for adsorption and desorption of gases. In this study, a column shale model was constructed to simulate methane diffusion and adsorption through shale rocks. The process of binary \(\hbox {CH}_{4}{-}\hbox {N}_{2}\) diffusion and adsorption was analyzed. A sensitivity study was performed to investigate the effects of pressure, temperature and permeability on diffusion and adsorption in shale rocks. The results show that methane gas diffusion and adsorption in shale is a slow process of dynamic equilibrium, which can be illustrated by the slope of a curve in \(\hbox {CH}_{4}\) mass variation. The amount of adsorption increases with the pressure increase at the low pressure, and the mass change by gas diffusion will decrease due to the decrease in the compressibility factor of the gas. With the elevated temperature, the gas molecules move faster and then the greater gas diffusion rates make the process duration shorter. The gas diffusion rate decreases with the permeability decrease, and there is a limit of gas diffusion if the permeability is less than \(1.0\,\times \,10^{-15}\, \hbox { m}^{2}\). The results can provide insights for a better understanding of methane diffusion and adsorption in the shale rocks so as to optimize gas production performance of shale gas reservoirs.     

Relative Permeability of Gas for Unconventional Reservoirs

Relative permeability of gas gains great significance in exploring unconventional gas. This paper developed a universal relative permeability model of gas, which is applicable for unconventional gas reservoirs such as coal, tight sandstone and shale. The model consists of the absolute relative permeability of gas and the gas slippage permeability. In the proposed model, the effects of water saturation and mean pore pressure on gas slippage permeability are taken into account. Subsequently, the evaluation of the model with existing model is done and then the validation of the model is made with data of tight sandstones, coals and shales from published literatures. The modeling results illustrate that a strong power-law relationship between relative permeability of gas and water saturation and the contribution of gas slippage permeability to relative permeability is determined by water saturation and mean pore pressure simultaneously. Furthermore, a sensitivity analysis of the impact of the parameters in the model is conducted and their effects are discussed.     

海陆过渡相页岩气藏不稳定渗流数学模型

海陆过渡相页岩常与煤层和砂岩呈互层状产出, 储层连续性较差、横向变化快、非均质性强, 水力压裂技术是其获得经济产量的关键手段. 然而, 目前缺乏有效的海陆过渡相页岩气藏不稳定渗流数学模型, 对其渗流特征分析及储层参数评价不利. 针对这一问题, 首先建立海陆过渡相页岩气藏压裂直井渗流数学模型, 其次采用径向复合模型来反映强非均质性, 采用Langmuir等温吸附方程来描述气体的解吸和吸附, 分别采用双重孔隙模型和边界元模型模拟天然裂缝和水力裂缝, 建立并求解径向非均质的页岩气藏压裂直井不稳定渗流数学模型, 分析海陆过渡相页岩气藏不稳定渗流特征, 并进行数值模拟验证和模型分析应用. 分析结果表明, 海陆过渡相页岩气藏不稳定渗流特征包括流动早期阶段、双线性流、线性流、内区径向流、页岩气解吸、内外过渡段、外区径向流及边界控制阶段. 将本模型应用在海陆过渡相页岩气试井过程中, 实际资料拟合效果较好, 其研究成果可为同类页岩气藏的压裂评价提供一些理论支撑, 具有较好应用前景.      

含水合物沉积物多相渗流特性研究进展

天然气水合物作为一种储量大、无污染的清洁能源近些年受到了广泛关注. 近20年来,中国进行了较大范围的陆海域天然气水合物储层勘探与储量预测.2017年,中国地质调查局牵头对南海神狐海域的天然气水合物进行了基于降压渗流原理的试验性开采.国内外已进行的水合物试采工程面临着气体产量低、出砂较多等问题,其最主要的原因之一是开发过程中沉积物内复杂多相渗流机理尚不明晰.本文综述了平行毛细管模型、Kozeny模型等广泛应用于天然气水合物开发渗流分析的理论模型,对比分析了水合物开发多尺度渗流过程模拟方法,简述了国内外含水合物沉积物渗透率测试、渗流过程中沉积物物性演变以及水合物开采室内模拟等方面的渗流实验进展,总结了矿场尺度的天然气水合物储层开采过程中产气数值模拟手段,展望了多相渗流模型、储层原位含水合物样品室内测试及结构与物性演化、矿场尺度数值模拟与水平井压裂技术等应用研究的未来方向与挑战.      

油页岩原位注热开采储层有效热解区变化规律分析

针对油页岩原位注热开采过程中储层有效热解区变化规律不清,实际热解效果无法准确判断难题,采用数值模拟方法,以抚顺油页岩储层为研究对象,建立了油页岩原位注热开采热流固耦合力学模型,与前人结果对比,验证了模型可靠性。重点考察水力压裂裂缝通道短路问题,分析得到了油页岩原位注热开采过程中储层有效热解区、储层有效热解区中地应力、注汽压力及沉降量随注热时间变化规律。结果表明,过热蒸汽沿水力压裂裂缝流动不会出现裂缝通道短路现象,过热蒸汽可通过水力压裂裂缝加速油页岩储层热解;采用过热蒸汽对流加热油页岩储层效率高,只需1年能使96%的油页岩储层达到热解所需温度;油页岩储层有效热解区中部形成应力集中区,最大地应力为21.6 MPa;热解后靠近注热井处岩层发生沉降,热解2年后最大沉降量达0.85 m。所得结论对现场油页岩原位注热开采有参考意义。