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

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

页岩有机质纳米孔隙气体吸附与流动规律研究

页岩储层孔隙结构复杂, 气体赋存方式多样. 有机质孔隙形状对受限空间气体吸附和流动规律的影响尚不明确, 导致难以准确认识页岩气藏气体渗流机理. 为解决该问题, 本文首先采用巨正则蒙特卡洛方法模拟气体在不同形状有机质孔隙(圆形孔隙、狭长孔隙、三角形孔隙、方形孔隙)内吸附过程, 发现不同形状孔隙内吸附规律符合朗格缪尔单层吸附规律, 分析了绝对吸附量、过剩吸附浓量、气体吸附参数随孔隙尺寸、压力的变化, 研究了孔隙形状对气体吸附的影响. 在明确不同形状有机质孔隙内气体热力学吸附规律基础上, 建立不同形状有机质孔隙内吸附气表面扩散数学模型和考虑滑脱效应的自由气流动数学模型, 结合分子吸附模拟结果研究了不同孔隙形状、孔隙尺寸有机质孔隙内吸附气流动与自由气流动对气体渗透率的贡献. 结果表明, 狭长孔隙内最大吸附浓度和朗格缪尔压力最高, 吸附气表面扩散能力最弱. 孔隙半径5 nm以上时, 吸附气表面扩散对气体渗透率影响可忽略. 本文研究揭示了页岩气藏实际生产过程中有机质孔隙形状对页岩气吸附和流动能力的影响机制.      

页岩气藏压裂水平井试井分析

页岩气藏资源丰富,开发潜力巨大,已成为目前研究的热点.与常规气藏相比,页岩气藏运移机制复杂,流动模式呈非线性,有必要考虑页岩气的吸附解吸,天然微裂缝的应力敏感性,人工裂缝内的非达西流等非线性因素对压裂水平井压力响应的影响. 基于双重介质和离散裂缝混合模型,分别采用Langmuir等温吸附方程描述吸附解吸,渗透率指数模型描述应力敏感,Forchheimer方程描述非达西效应,建立页岩气藏压裂水平井数值试井模型. 运用伽辽金有限元法对模型进行求解.根据试井特征曲线,划分流动阶段,着重分析非线性因素对压力响应的影响.结果表明:页岩气藏压裂水平井存在压裂裂缝线性流、压裂裂缝径向流、地层线性流、系统径向流及封闭边界影响5 种流动阶段.吸附解吸的影响发生窜流之后,Langmuir吸附体积增大,拟压力导数曲线凹槽更加明显,系统径向流出现时间与压力波传播到边界时间均延迟;天然裂缝系统的应力敏感性主要影响试井曲线的晚期段,拟压力和拟压力导数曲线均表现为上翘,应力敏感效应越强,上翘幅度越大;高速非达西效应对早期段影响较大,非达西效应越强,拟压力降幅度越大,试井曲线上翘.与解析解的对比以及矿场实例验证了模型的正确性与适用性.      

页岩气开采中的井底压力并行计算

根据页岩气流动特点建立了考虑混合气体高压物性参数、渗透率与孔隙度随压力变化的页岩气流动方程,通过定义拟压力函数将页岩气流动的偏微分方程线性化。针对页岩气开发采用水平井多段压裂技术,采用Newman乘积原理得到地层拟压力流动方程的解析解表达式。依据解析解的特征将解析解分解成适合并行计算的无限求和及积分形式,提出了一套基于CUDA的页岩气地层压力算法,将地层拟压力函数解析解划分为多个并行度较高的步骤,利用GPU的并行计算能力,设计每个步骤的CUDA核函数,在英特尔i3 540CPU(3.07GHz主频,4GB内存)和NVIDIA的GTX 550显卡上,计算了页岩气的井底压力,分析了井底压力特征。实验结果表明,页岩吸附影响曲线变化剧烈程度,而扩散主要影响曲线发生变化的时间,在GPU上的页岩气压力计算可达近40倍的加速比。     

矩形页岩气藏压裂直井解析解产能评价

为了准确评价矩形渗流区域内页岩气藏压裂直井的产能,加强气井生产过程中的科学性管理。通过运用Fick扩散定律和达西渗流规律,建立页岩气藏压裂直井基质扩散方程和裂缝渗流方程。结合时间叠加原理,通过Green函数给出矩形封闭边界和矩形定压边界下无限导流垂直裂缝井不稳定渗流的解析解。研究结果明确了页岩气藏矩形边界压裂直井的渗流特征及不同扩散特征下的压力响应,提出了矩形页岩气藏压裂直井产能评价方法,达到了快速评价及准确预测气井产能的目的。矩形页岩气藏的流动能分为7个阶段。通过将本文提出的矩形封闭边界模型与实际生产数据拟合,得到单井产量平均拟合误差为11.6%,压力平均拟合误差为8.3%,累计产气量平均拟合误差为7.12%。     

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

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

页岩气在超低渗介质中的渗流行为

基于页岩气藏复杂孔隙结构和页岩气在纳米孔隙表面的滑脱和吸附-解吸附等现象,通过引入表观渗透率,修正Darcy渗流模型,得到了页岩气渗流本构方程. 将计算结果与Darcy模型计算结果进行了比照,结果表明,在产量定常情形下,基于非Darcy渗流模型得到的井底压力高于Darcy流结果; Darcy流模型得到的压力衰减速度较快,改进后的渗流模型更能准确描述与表征页岩气渗流过程;研究结果可为页岩气藏的经营与管理提供基础参数.      

页岩膨胀性能评价与拟合模型建立

近年来,页岩气作为一种非常规天然气备受关注。在页岩气藏钻井与生产过程中,页岩膨胀是非常关键的问题之一。在页岩矿物组成分析的基础上,采用实际页岩样品和模拟页岩样品,利用高温高压页岩膨胀仪测量了页岩样品在蒸馏水中的膨胀性能,分析了温度、压力与粘土含量等因素对页岩膨胀性能的影响。结果表明,温度、压力对页岩的膨胀性能均有影响。在实验温度和压力区间内,温度越高、压力越低,页岩在蒸馏水中的膨胀越严重。页岩中的粘土矿物是影响页岩膨胀的最主要因素,粘土含量越高,页岩的膨胀率越大。页岩膨胀性能可以用非线性Log-Normal方程来定量描述。     

页岩气开采中的若干力学前沿问题

页岩气的开采涉及破裂和收集输运两个关键过程.如何实现2000,m以下、复杂地应力作用下、多相复杂介质组分的页岩层内网状裂纹的形成,同时将孔洞、缝隙中的游离、吸附气体进行高效收集,涉及到诸多的核心力学问题.这一工程过程涵盖了力学前沿研究的诸多领域:介质和裂纹从纳米尺度到千米尺度的空间跨越,游离、吸附气体输运过程中微秒以下的时间尺度事件到历经数年开采的时间尺度跨越,不同尺度上流体固体的相互作用,以及压裂过程中通过监测信息反演内部破坏状态等.针对近年来我们国家页岩气勘探开发工作所取得的成就及后续发展中面临的前沿力学问题,在综合介绍页岩气藏的基本特征和开发技术的基础上,以页岩气开采中的若干力学前沿问题为主线,从页岩力学性质及其表征方法、页岩气藏实验模拟技术、页岩气微观流动机制及流固耦合特征、水力压裂过程数值模拟方法、水力压裂过程微地震监测技术、高效环保的无水压裂技术等6个方面的最新研究进展进行了总结和展望,结合页岩气藏开发的工程实践, 深入探究了其中力学关键问题,以期对从事页岩气领域的开发和研究的从业人员提供理论基础, 同时,该方面的内容对力学学科、尤其是岩土力学领域的科研工作也具有重要指导价值.      

Mechanical frontiers in shale-gas development

页岩气的开采涉及破裂和收集输运两个关键过程.如何实现2000,m以下、复杂地应力作用下、多相复杂介质组分的页岩层内网状裂纹的形成,同时将孔洞、缝隙中的游离、吸附气体进行高效收集,涉及到诸多的核心力学问题.这一工程过程涵盖了力学前沿研究的诸多领域：介质和裂纹从纳米尺度到千米尺度的空间跨越,游离、吸附气体输运过程中微秒以下的时间尺度事件到历经数年开采的时间尺度跨越,不同尺度上流体固体的相互作用,以及压裂过程中通过监测信息反演内部破坏状态等.针对近年来我们国家页岩气勘探开发工作所取得的成就及后续发展中面临的前沿力学问题,在综合介绍页岩气藏的基本特征和开发技术的基础上,以页岩气开采中的若干力学前沿问题为主线,从页岩力学性质及其表征方法、页岩气藏实验模拟技术、页岩气微观流动机制及流固耦合特征、水力压裂过程数值模拟方法、水力压裂过程微地震监测技术、高效环保的无水压裂技术等6个方面的最新研究进展进行了总结和展望,结合页岩气藏开发的工程实践, 深入探究了其中力学关键问题,以期对从事页岩气领域的开发和研究的从业人员提供理论基础, 同时,该方面的内容对力学学科、尤其是岩土力学领域的科研工作也具有重要指导价值.     

Apparent Permeability and Gas Flow Behavior in Carboniferous Shale from the Qaidam Basin,China: An Experimental Study

Shale can act as an unconventional gas reservoir with low permeability and complex seepage characteristics. Study of the apparent permeability and percolation behavior of shale gas is important in understanding the permeability of shale reservoirs, to evaluate formation damage, to develop gas reservoirs, and to design wells. This study simulated methane percolation at 298.15 K under inlet pressures ranging from 0.2 to 4 MPa and a constant outlet pressure of 0.1 MPa to investigate shale gas percolation behavior and apparent permeability. Five representative shale cores from the Carboniferous Hurleg and Huitoutala formations in the eastern Qaidam Basin, China, were analyzed. Each experiment measured the volume flow rate of methane and the inlet pressure. Pseudopressure approach was used to analyze high-velocity flow in shale samples, and apparent permeability at different pressures was calculated using the traditional method. A nonlinear apparent permeability model that considers diffusion and slippage is established from theory and experimental data fitting, and the shale gas flow characteristics affected by slippage and diffusion are analyzed. The results indicate that the pseudopressure formulation that considers the effect of gas properties on high-velocity flow produces a more accurate linear representation of the experimental data. The apparent gas permeability of shale consists of contributions from Darcy permeability, slippage, and diffusion. The apparent permeability and gas flow behavior in the studied shales strongly depended on pressure. The diffusion contribution increased greatly as pressure decreased from 2 to 0.2 MPa, and the smaller the shale permeability, the greater the relative contribution of diffusion flow. At pressures greater than 2 MPa, slip flow contributes \(\sim \)20% of the total flux, Darcy flow contributes up to 70%, and diffusion makes only a minor contribution. This study provides useful information for future studies of the mechanism of shale gas percolation and the exploration and development of Qaidam Basin shale gas specifically.     

页岩气藏的双重介质-离散裂缝模型

考虑页岩气藏开发中渗流的多尺度效应,提出了一个基于裂缝-孔隙双重介质的离散裂缝模型．在该模型中,基质、天然裂缝和人工压裂裂缝采用各自控制方程独立计算,不同介质之间通过流量交换相互关联．为分析模型可靠性,分别和基于渗透率粗化及压裂裂缝导流能力无穷大的模型对比．数值算例显示,伴随着网格细分,该模型与精确渗透率粗化模型具有相同计算精度,两者收敛速度均较快,但该模型易推广到多相流动问题,而等压模型对产量将有所高估．研究了地质参数和工艺参数对气井产量的影响规律．计算结果表明天然裂缝渗透率及基质孔隙扩散系数对产气速率有着重要影响,产气速率伴随着人工压裂裂缝导流能力、长度以及数目的增加而增加,但是增加幅度会逐步趋缓．     

Comprehensive Seepage Simulation of Fluid Flow in Multi-scaled Shale Gas Reservoirs

Shale gas seepage behaviour is a multi-field/-scale problem and makes transient pressure analysis a very challenging task. Non-Darcy flow in nanopores is prominent due to the broken of continuity hypothesis. Slippage effect and Knudsen diffusion are two important seepage mechanisms in nanopores, while recent studies show surface diffusion is another important transporting mechanism on surface of nanopores. Porous kerogen system contains large amounts of dissolved gas, which should not be overlooked. In this study, a comprehensive mathematical model was established by pseudo-quadruple porosity medium conception, coupling the effects of slippage flow, Knudsen diffusion, surface diffusion, ad-/desorption and gas transferring from kerogen to nanopore system, while fluid flow in fractures/macropores is described by Darcy’s law. Transient pressure behaviours of a multiple fractured horizontal well in box-shaped shale gas reservoir were studied, with nine possible flow regimes divided and parameters sensitivity analysed. Adsorbed constant and dissolved constant were defined to reflect the amount of adsorbed gas and dissolved gas, respectively. Research shows that adsorbed gas and dissolved gas are two important gas storage forms, neither of which should be neglected. The study can not only help us understand fluid flow mechanisms in nanopores from microscopic perspective, but enable us to analyse production performance and determine key operational parameters from macroscopic perspective.     

Microscopic mechanism of periodical electroosmosis in reservoir rocks

Based on the electric double layer (EDL) theory and the momentum equation governing the electroosmosis flow, this paper presents an analytical solution to the periodical electroosmosis with a parallel straight capillary bundle model of reservoir rocks to reveal the microscopic mechanism of the electroosmotic flows in rocks. The theory shows that both the frequency dispersion characteristics of the macroscopic electroosmotic Darcy velocity in unsealed rocks and the electroosmotic pressure coefficient in sealed rocks depend on the porosity and electrochemical properties of reservoir rocks. The mathematical simulation indicates that the distribution of the periodical electroosmotic velocity is wavelike in the rock pore. The greater the porosity is, the greater electroosmotic the Darcy velocity and the smaller electroosmotic pressure coefficient are generated. The module values of the electroosmotic Darcy velocity and the electroosmotic pressure coefficient increase with the decreasing solution concentration or the increasing cation exchange capacity without affecting the phase of the electroosmotic Darcy velocity.     

Investigation on nonlinear multi-scale effects of unsteady flow in hydraulic fractured horizontal shale gas wells

A unified mathematical model is established to simulate the nonlinear unsteady percolation of shale gas with the consideration of the nonlinear multi-scale effects such as slippage, diffusion, and desorption. The continuous inhomogeneous models of equivalent porosity and permeability are proposed for the whole shale gas reservoir including the hydraulic fracture, the micro-fracture, and the matrix regions. The corresponding semi-analytical method is developed by transforming the nonlinear partial differential governing equation into the integral equation and the numerical discretization. The non-linear multi-scale effects of slippage and diffusion and the pressure dependent effect of desorption on the shale gas production are investigated.     

A Fractal Model for Gas–Water Relative Permeability in Inorganic Shale with Nanoscale Pores

A reliable gas–water relative permeability model in shale is extremely important for the accurate numerical simulation of gas–water two-phase flow (e.g., fracturing fluid flowback) in gas-shale reservoirs, which has important implication for the economic development of gas-shale reservoir. A gas–water relative permeability model in inorganic shale with nanoscale pores at laboratory condition and reservoir condition was proposed based on the fractal scaling theory and modified non-slip boundary of continuity equation in the nanotube. The model not only considers the gas slippage in the entire Knudsen regime, multilayer sticking (near-wall high-viscosity water) and the quantified thickness of water film, but also combines the real gas effect and stress dependence effect. The presented model has been validated by various experiments data of sandstone with microscale pores and bulk shale with nanoscale pores. The results show that: (1) The Knudsen diffusion and slippage effects enhance the gas relative permeability dramatically; however, it is not obviously affected at high pressure. (2) The multilayer sticking effect and water film should not be neglected: the multilayer sticking would reduce the water relative permeability as well as slightly decrease gas relative permeability, and the film flow has a negative impact on both of the gas and water relative permeability. (3) The increased fractal dimension for pore size distribution or tortuosity would increase gas relative permeability but decrease the water relative permeability for a given saturation; however, the effect on relative permeability is not that notable. (4) The real gas effect is beneficial for the gas relative permeability, and the influence is considerable when the pressure is high enough and when the nanopores of bulk shale are mostly with smaller size. For the stress dependence, not like the intrinsic permeability, none of the gas or water relative permeability is sensitive to the net pressure and it can be ignored completely.     

Can Moderately Rarefied Gas Transport Through Round and Flat Tight Channels of Fractured Porous Media be Described Accurately?

We study the generation and flow of foam through rough-walled, fractured marble rocks that mimic natural fracture systems in carbonate reservoirs. Flow was isolated to the fracture network because of the very low rock permeability of the marble samples and foam generated in situ during co-injection of surfactant solution and gas. The foam apparent viscosities were calculated at steady pressure gradients for a range of gas fractions, and similar to foam flow in porous media, we identified two flow regimes for foam flow in fractures: a high-quality flow regime only dependent on liquid velocity and a low-quality flow regime determined by the gas and liquid velocities. Variations in local fluid saturation during co-injection were visualized and quantified using positron emission tomography combined with computed tomography.