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

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

Numerical Simulation of Hydraulic Fracturing Water Effects on Shale Gas Permeability Alteration

Hydraulic fracturing has been recognized as the necessary well completion technique to achieve economic production from shale gas formation. However, following the fracturing, fluid–wall interactions can form a damaged zone nearby the fracture characterized by strong capillarity and osmosis effects. Here, we present a new reservoir multi-phase flow model which includes these mechanisms to predict formation damage in the aftermath of the fracturing during shut-in and production periods. In the model, the shale matrix is treated as a multi-scale porosity medium including interconnected organic, inorganic slit-shaped, and clay porosity fields. Prior to the fracturing, the matrix holds gas in the organic and the inorganic slit-shaped pores, water with dissolved salt in the inorganic slit-shaped pores and the clay pores. During and after fracturing, imbibition causes water invasion into the matrix, and then, the injected water–clay interaction may lead to clay-swelling pressure development due to osmosis. The swelling pressure gives additional stress to slit-shaped pores and cause permeability reduction in the inorganic matrix. We develop a simulator describing a system of three pores, two phases (aqueous and gaseous phases), and three components (\(\hbox {H}_{2}\hbox {O}, \hbox {CH}_{4}\), and salt), including osmosis and clay-swelling effect on the permeability. The simulation of aqueous-phase transport through clay shows that high swelling pressure can occur in clays as function of salt type, salt concentration difference, and clay-membrane efficiency. The new model is used to demonstrate the damage zone characteristics. The simulation of two-phase flow through the shale formation shows that, although fracturing is a rapid process, fluid–wall interactions continue to occur after the fracturing due to imbibition mechanism, which allows water to penetrate into the inorganic pore network and displace the gas in-place near the fracture. This water invasion leads to osmosis effect in the formation, which cause clay swelling and the subsequent permeability reduction. Continuing shale–water interactions during the production period can expand the damage zone further.     

水力压裂解析模型裂缝扩展参数敏感性分析

水力压裂形成复杂裂缝网络是致密储层油气开采的重要技术,掌握水压裂缝扩展机理是控制压裂行为和优化压裂效果的关键．水压裂缝动态扩展行为涉及储层岩体、注入压裂液、压裂实施工艺等方面,其中水力压裂扩展时间、压裂液流体动力粘度系数、压裂液流体注入流速、储层岩石剪切模量成为决定裂缝扩展长度和裂缝开度的重要因素．本研究采用KGD、PKN两类等高解析模型对主控因素的参数敏感性进行分析,直观、快速、可靠地获得水压裂缝扩展长度、张开度动态演化行为的量化数值．研究发现,压裂持续开展过程中水压裂缝扩展长度呈线性增长、开度逐渐趋于稳定,高流体动力粘度导致裂缝难扩展、形成较大裂缝开度,通过增加压裂液流体注入流速可同时增加裂缝扩展长度和开度,较高的岩石剪切模量将降低水压裂缝的开度．通过对比两类解析模型在不同参数下的水压裂缝扩展结果,分析压裂参数与裂缝扩展的相关性和敏感系数,讨论水力压裂解析模型的裂缝扩展参数敏感性．     

MODEL AND METHOD OF WELL TEST ANALYSIS FOR WELLS WITH VERTICAL FRACTURE

NomenclatureB———volumecoefficientofformationfluid,dimensionless;C———wellborestoragecoefficient,m3 MPa;Cft ———comprehensivecompressioncoefficientoffluidinfracture ,1 MPa;Ct ———comprehensivecompressioncoefficientofformationfluid,1 MPa ;h———heightofformation ,m ;K———formationpermeability ,μm2 ;Kf ———fracturepermeability ,μm2 ;KfWf ———fractureconductivity ,μm2 ·cm ;Lf ———halfoffracturelength,m ;N———thetotalnumberofselectedtimepoint,dimensionless;P———for…     

基于离散裂缝的多段压裂水平井数值试井模型及应用

水平井压裂技术已经成为开发低渗透油气藏、页岩气藏和致密气场等非常规油气藏的关键技术。基于离散裂缝模型,对裂缝进行简化,建立了二维多段压裂水平井有限导流数值试井模型,利用有限元方法求解模型,获得多段压裂水平井试井理论曲线和压力场特征。分析表明:多段压裂水平井的试井理论曲线一共分为七个阶段:井筒储存段、裂缝线性流段、裂缝-地层双线性流段、裂缝干扰段、地层线性流段、系统径向流段和边界作用段,其中裂缝-地层双线性流段和裂缝干扰是其典型特征。分析了裂缝数量、裂缝间距、裂缝不对称、裂缝不等长和裂缝部分缺失等因素对试井理论曲线的影响,结果表明:裂缝数量和裂缝间距对试井理论曲线的影响最大。较多的裂缝、较大裂缝间距、对称的裂缝和等长的裂缝有利于降低压裂水平井井底的流动阻力,提高产能。将建立的数值试井模型应用于四川盆地一口多段压裂水平井的压力恢复测试的数值试井解释,结果表明:本文建立的模型可以较好的拟合压力恢复测试数据,可以获得裂缝的导流能力和裂缝长度,为压裂效果评价和压裂设计提供指导。      

3-D Fracture Propagation Simulation and Production Prediction in Coalbed

NｏｍｅｎｃｌａｔｕｒｅＡ(ｘ,ｔ)—Ｃｒｏｓｓ＿ｓｅｃｔｉｏｎａｒｅａｏｆｆｒａｃｔｕｒｅａｔｔｉｍｅｔａｎｄｓｉｔｅｘ,ｍ2 ;Ｂｇ,Ｂｗ —Ｃｏａｌｂｅｄｇａｓａｎｄｗａｔｅｒｖｏｌｕｍｅｆａｃｔｏｒ,ｄｉ ｍｅｎｓｉｏｎｌｅｓｓ;ｂ—Ｌａｎｇｍｕｉｒｆａｃｔｏｒ,(ＭＰａ) - 1 ;Ｃｇ,ＣＲ,Ｃｗ —Ｃｏａｌｂｅｄｇａｓ,ｒｏｃｋａｎｄｗａｔｅｒｃｏｍ ｐｒｅｓｓｉｂｉｌｉｔｙｃｏｅｆｆｉｃｉｅｎｔ,(ＭＰａ) - 1 ;Ｃｐ —Ｓａｎｄｒａｔｉｏ ,ｄｉｍｅｎｓｉｏｎｌｅｓｓ;Ｃｔ(ｘ ,ｔ)—Ｔｏｔａｌｆｒａｃｔｕｒｉｎ…     

页岩水力压裂的关键力学问题

页岩气是指以吸附和游离时而还有流体相的形态状态赋存于泥页岩中的非常规天然气。页岩气开采成为我国绿色能源开发的新领域。尽管北美页岩气革命取得了成功,但目前仍仅有预期产量5~15%的油气采收率,问题出在什么地方？水力压裂被认为是提高采收率的关键一环,但水力压裂过程中复杂缝网的形成和力学控制机理尚不清楚,这给力学家提出了巨大的挑战和机遇。结合本课题组近期研究成果,本文从理论、计算和实验三个方面对页岩水力压裂中的关键基础力学问题进行介绍和总结。主要内容包括发展页岩人工裂缝扩展的大型物理模型实验,建立页岩本构模型和断裂力学理论,发展耦合断裂力学和流体力学的裂缝网扩展数值模拟方法。页岩水力压裂研究对发展断裂力学的实验技术、理论模型和数值模拟方法起到推动作用,对提高我国页岩气高效开采技术具有重要意义。     

页岩储层体积压裂的地应力变化研究

页岩储层属于致密超低渗透储层,需改造形成复杂缝网才有经济产能.体积压裂是页岩储层增产改造的主要措施,而地应力场特别是水平主应力差值是体积压裂的关键控制因素. 理论研究表明:(1)当初始两向水平主应力差较小时,容易形成缝网,反之不易产生缝网;(2)人工裂缝的形成能够改变地层初始应力场. 因此应在前人研究的基础上优化设计压裂方式,以克服和翻转初始水平主应力差值,产生体积缝网.基于此,建立了页岩气藏水平井体积压裂数值模型,模型中采用多孔介质流固耦合单元模拟页岩基质的行为,采用带有孔压的"cohesive"单元描述水力裂缝的性质,模型对"Texas Two-Step" 压裂方法进行了数值模拟,模拟结果得到了压裂过程中地层应力场的分布及其变化,模拟结果和解析公式计算结果吻合良好.模拟结果表明:(1)裂缝的产生减弱了地层应力场的各向异性;(2 对于低水平应力差页岩储层,采用"Texas Two-Step"压裂方法可以产生缝网. 对于采用"Texas Two-Step"压裂方法无法产生缝网的高应力差页岩储层,提出了三次应力"共振" 和四次应力"共振" 压裂方法并进行了数值模拟,模拟结果得到了压裂过程中页岩储层应力场的分布及其变化,得到了缝网形成的区域,模拟结果表明:(1)对于高应力差页岩储层,采用"Texas Two-Step" 压裂方法无法产生缝网;(2)对于高应力差页岩储层,三次应力"共振" 和四次应力"共振"压裂方法是有效的体积压裂缝网形成的方法.      

Experimental Method to Simulate Coal Fines Migration and Coal Fines Aggregation Prevention in the Hydraulic Fracture

A simple and effective experimental method is proposed to simulate coal fines migration through the proppant pack; such migration inevitably occurs during the process of fracturing fluid flowback or dewatering and gas production in coalbed methane (CBM) reservoirs. The damage to conductivity caused by coal fines migration in the pack and the factors affecting such migration are analyzed. A dispersion agent of coal fines applicable to hydraulic fracturing in CBM is optimized, consequently solving the problem of coal fines aggregation and retention in the proppant pack. Discharging coal fines with water or water-based fracturing fluid from the proppant pack can be difficult because of the adsorption and hydrophobicity of coal fines. Thus, coal fines are likely to aggregate and be retained in the proppant pack, thereby resulting in pore throat plugging, which causes serious damage to fracture conductivity. Two percent coal fines can reduce propped fracture conductivity by 24.4 %. The mobility and retention of coal fines in the proppant pack are affected by proppant size, proppant type, flowback rate, and coal fines property. When flowback rate exceeds the critical value, coal fines can be discharged from the pack, consequently reducing damage to propped fracture conductivity. More importantly, the steady discharging of coal fines requires steady dewatering and gas production to avoid flow shock, which causes pressure disturbance to drive coal fines in a remote formation. The optimized dispersant FSJ-02 employed in this paper can effectively change the wettability and surface potential of coal fines to improve their suspension and dispersion in water-based fracturing fluid. The recovery rate of coal fines increased by 31.5 %, whereas conductivity increased by 13.3 %.     

On the problem of fluid leakoff during hydraulic fracturing

While a hydraulic fracture is propagating, fluid flow and associated pressure drops must be accounted for both along the fracture path and perpendicularly, into the formation that is fractured, because of fluid leakoff. The accounting for the leakoff shows that it is the main factor that determines the crack length. The solved problem is useful for the technology of hydraulic fracturing and a good example of mass transport in a porous medium. To find an effective approach for the solution, the thin crack is represented here as the boundary condition for pore pressure spreading in the formation. Earlier such model was used for heat conduction into a rock massif from a seam under injection of hot water. Of course, the equations have other physical sense and mathematically they are somewhat different. The new plane solution is developed for a linearized form that permits the application of the integral transform. The linearization itself is analogous to the linearization of the natural gas equation using the real gas pseudo-pressure function and where the flux rates are held constant and approximations are introduced only into the time derivatives. The resulting analytical solution includes some integrals that can be calculated numerically. This provides rigorous tracking of the created fracture volume, leakoff volume and increasing fracture width. The solutions are an advance over existing discreet formulations and allow ready calculations of the resulting fracture dimensions during the injection of the fracturing fluid.     

低渗透油藏爆炸压裂产能分析及压裂改造区域优化设计

基于爆炸压裂裂缝分布规律,提出爆炸压裂缝网双重介质复合流动产能模型,应用Laplace变换Stehfest数值反演,得到了定压条件下封闭外边界低渗透油藏爆炸压裂生产井产能表达式。在模型正确性验证的基础上结合某低渗透油藏储层特征参数研究了爆炸压裂改造区域参数对封闭边界油藏产量的影响,同时对爆炸压裂改造改造体积优化设计进行了研究。研究结果表明,爆炸压裂改造区域半径主要影响生产中期产能,改造区域渗透率对生产早期和中期影响比较大,且对于实例油藏爆炸压裂改造比为0.1时效果最好。     

Experiment and simulation of slurry flow in irregular channels to understand proppant transport in complex fractures

Slurry flow and proppant placement in irregular fractures are crucial to evaluate hydraulic fracturing stimulation but need to be better understood. This study aims to investigate how irregular fracture affects proppant transport and distribution using laboratory experiments and micro-scale numerical models. The unresolved method of the computational fluid dynamics (CFD) and the discrete element method (DEM) considers Saffman lift force, Magnus force, and virtual mass force to accurately capture the frequent interaction between proppant and slickwater. Experimental results validated the reliability of the optimized CFD-DEM model and calibrated primary parameters. The effects of crack height and width, bending angle, and distance between the crack and inlet on particle distribution were studied. The results indicated that the improved numerical method could rationally simulate proppant transport in fractures at a scale factor. The small crack height causes downward and upward flows, which wash proppant to the fracture rear and form isolated proppant dunes. A wider region in the fracture is beneficial to build up a large dune, and the high dune can hinder particle transport into the fracture rear. When the crack is close to the inlet, the primary fracture without proppants will close to hinder gas production. The smaller the bending angle, the smaller the proppant dune. A regression model can precisely predict the dune coverage ratio. The results fundamentally understand how complex fractures and natural cracks affect slurry flow and proppant distribution.     

COHESIVE ZONE FINITE ELEMENT-BASED MODELING OF HYDRAULIC FRACTURES

Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.     

页岩气高效开采的可压裂度和射孔簇间距预测

页岩储层的可压裂性是影响页岩气产量的关键因素. 本文基于断裂力学理论, 以高围压下岩石层理弱面的剪切破坏为主要研究对象, 依据岩石抗拉强度和层理弱面抗剪强度的比值关系, 首先提出了可压裂度的概念, 给出了无量纲的定性曲线图, 涵盖了岩石脆性矿物质含量, 粘性主导和韧性主导裂缝尖端流体压强、射孔簇分布间距的综合地质与工程因素. 接着提出了一种新的表征高围压下页岩可压裂度的无量纲参数, 在保证达到充分解吸附的最小压裂间距前提下, 依据该参数可计算水平井压裂中的射孔簇间距, 可作为工程参考指标. 本文将断裂力学理论结合水力压裂高效开采页岩气工程, 具有力学理论意义和工程应用前景.      

Problem of propagation of a gas fracture in a porous medium

The problem of gas fracture formation in a porous medium is investigated. An inertialess viscous polytropic gas flow along the fracture is considered. The assumption of small fracture width with respect to the height and length makes it possible to adopt the vertical plane cross-section hypothesis on the basis of which the dependence of the gas pressure inside the fracture on its width can be reduced to a linear law. Initially, the soil surrounding the fracture is soaked with oil-bearing fluid. During fracturing the reservoir gas penetrates into the soil mass and displaces the fluid. A closed system of equations, which describes the evolution of the fracture opening, the depth of gas penetration into the reservoir, and the gas velocities inside the fracture, is constructed. The limiting regimes of gas seepage into the surrounding reservoir are considered and a one-parameter family of self-similar solutions of the system is given for each. The asymptotics of the solution in the neighborhood of the fracture nose is investigated and analytic expressions for the fracture length are obtained. The solution of the problem of gas fracture is compared with the hydraulic fracturing problem in an analogous formulation within the framework of the plane cross-section hypothesis.     

Simulation of Pressure Transient Behavior for Asymmetrically Finite-Conductivity Fractured Wells in Coal Reservoirs

Based on Fick’s law in matrix and Darcy flow in cleats and hydraulic fractures, a new semi-analytical model considering the effects of boundary conditions was presented to investigate pressure transient behavior for asymmetrically fractured wells in coal reservoirs. The new model is more accurate than previous model proposed by Anbarci and Ertekin, SPE annual technical conference and exhibition, New Orleans, 27–30 Sept 1998 because new model is expressed in the form of integral expressions and is validated well through numerical simulation. (1) In this paper, the effects of parameters including fracture conductivity, coal reservoir porosity and permeability, fracture asymmetry factor, sorption time constant, fracture half-length, and coalbed methane (CBM) viscosity on bottomhole pressure behavior were discussed in detail. (2) Type curves were established to analyze both transient pressure behavior and flow characteristics in CBM reservoir. According to the characteristics of dimensionless pseudo pressure derivative curves, the process of the flow for fractured CBM wells was divided into six sub-stages. (3) This paper showed the comparison of transient steady state and pseudo steady state models. (4) The effects of parameters including transfer coefficient, wellbore storage coefficient, storage coefficient of cleat, fracture conductivity, fracture asymmetry factor, and rate coefficient on the shape of type curves were also discussed in detail, indicating that it is necessary to keep a bigger fracture conductivity and fracture symmetry for enhancing well production and reducing pressure depletion during the hydraulic fracturing design.     

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

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

基于混合方法的二维水力压裂数值模拟

水力压裂在页岩气开采中被广泛使用,采用数值方法研究压裂机理具有重要意义.基于连续-非连续单元法(CDEM) 和中心型有限体积法(FVM),提出解决水力压裂流固耦合问题的二维混合数值计算模型.该混合模型中,使用CDEM 求解应力场和裂缝扩展过程,使用FVM 求解裂隙渗流场.应力场裂缝扩展和渗流场均使用显式迭代求解, 并通过相互之间数据交换实现流固耦合.通过与KGD 理论模型进行对比, 验证数值模型的正确性.通过与颗粒离散元数值结果进行对比,验证数值模型的有效性.通过计算复杂缝网压裂模型,研究水力压裂机理,并说明该数值模型在水力压裂模拟中具有很好的前景.      

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

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

Modeling of Coal Fine Migration During CBM Production in High-Rank Coal

During CBM (coalbed methane) production, the interaction of coal fracture surface with water flow commonly generates and starts coal fine flow. Part of flowing coal fines deposit in coal fracture system due to water production reduction and methane production increase. The fine sedimentation results in the reduction of coal permeability and well productivity. Despite the increasing awareness of the importance of fine migration, limited research has been carried out on the flow model of coal fine coupled with water and gas. In this paper, a flow model of coal fine is established coupled with water and gas flow, taking coal fine generation, migration and sedimentation process into consideration. Then, case simulations are conducted to illustrate effects of water production schedule, permeability performance and gas content on production performance in flow model. The simulation results indicate that methane rate with the lowest initial water rate is observed to have the highest production in late production period. This is mainly due to the reason that the low water flow cannot generate and start the flow of coal fine. Further, the case with high initial water production has faster gas and water flow rate, thus higher coal fine generation rates, which can improve well productivity at earlier production period. As water production declines quickly, both permeability and production performance decrease, which leads to the loss of well productivity. Meanwhile, higher gas content will lead to a faster water production decline at late production period. This indicates that a portion of coal fines plugged in the fracture as water production deceases and the CBM reservoir with high gas content should not adopt a high initial water production schedule.