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

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

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

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

低渗透煤层气有限导流垂直裂缝井压力动态分析

作为典型的致密多孔介质,煤岩储层已被证实存在启动压力梯度.根据煤层气垂直裂缝井的双线性流动机制,综合考虑启动压力梯度和井筒储存效应的影响,建立了一个新的低渗透煤层气有限导流垂直裂缝井双线性流动数学模型,采用Laplace变换和Stehfest数值反演方法对数学模型进行了求解,并分析了无因次启动压力梯度等参数对无因次井底压力及其导数曲线的影响规律.分析结果表明:典型的低渗透煤层气垂直裂缝井双线性流动曲线可划分为早期续流段、双线性流段、煤层线性流段、过渡流段和煤层边界线性流段5个特征阶段,其中由于启动压力梯度存在的影响,无因次井底压力及其导数曲线自煤层线性流段开始出现明显上翘,且启动压力梯度值越大,曲线上翘趋势越明显;此外,煤层边界线性流段呈现为单位斜率的直线,而非1/2斜率的线性流段直线.这些结果表现出启动压力梯度对低渗透煤层气垂直裂缝井双线性流动的影响,可用于指导现场煤层气井试井分析.     

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

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

井筒气液两相流对致密气压裂水平井试井的影响

多段压裂水平井技术是目前开采致密气最常用的方法之一,在致密气压裂水平井试井测试中常常伴随着一定的产水量,井筒气液两相流会增加井筒流体的流动阻力,加大井筒流体流动对试井解释的影响.为了明确井筒气液两相流对致密气藏压裂水平井试井的影响,提高产水致密气压裂水平井的试井解释精度,建立了一种井筒气液两相流与地层渗流耦合的试井模型,采用数值方法对模型进行求解,获得了考虑井筒气液两相流的压裂水平井试井理论曲线、压力场分布及裂缝产量分布.研究结果表明:井筒气液两相流会增加试井理论曲线中压力和压力导数值,造成靠近入窗点的压力扩散要快于远离入窗点的压力扩散,引起靠近入窗点的裂缝产量要高于远离入窗点的裂缝产量.现场实例分析进一步说明,不考虑井筒两相流可能会对产水压裂水平井的试井解释结果产生很大误差,主要表现为水平井筒假设为无限大导流能力会使得拟合得到的表皮系数偏大,将测试点视为入窗点会使得拟合得到的原始地层压力偏小.所建立的考虑井筒两相流的压裂水平井试井模型为产水致密气井试井资料的正确解释提供了重要技术保障.      

低渗油层压裂水平井两相流研究

依据压裂水平井不同流动区域的流动规律, 将压裂水平井的渗流分为裂缝中的高速非达西流 动区、裂缝控制影响的椭圆渗流区、远离裂缝的基质非达西渗流区, 考虑启动压力梯度的影 响, 对压裂水平井两相渗流进行了分析, 得到了低渗透油层压裂水平井的产量公式. 研究结 果表明, 裂缝的导流能力越大, 压裂水平井的产量越高. 但随着开采时间的增加, 其产量递减幅度越 大; 压裂裂缝长度越小, 压裂水平井的初始产量越高. 但随着生产时间的推移, 压裂裂缝的 长度越大, 产量的递减幅度越小; 中间裂缝长, 两翼裂缝短的情况下, 压裂水平井的产量最 高.     

考虑多重运移机制耦合页岩气藏压裂水平井数值模拟

页岩作为典型的微纳尺度多孔介质,游离气与吸附气共存,传统的达西定律已无法准确描述气体在页岩微纳尺度的运移规律.基于双重介质模型和离散裂缝模型构建页岩气藏分段压裂水平井模型,其中基岩中考虑气体的黏性流、Knudsen 扩散以及气体在基岩孔隙表面的吸附解吸,吸附采用Langmuir等温吸附方程;裂缝中考虑黏性流和Knudsen扩散,在此基础上建立基岩-裂缝双重介质压裂水平井数学模型并采用有限元方法对模型进行求解.结果表明,基岩固有渗透率越小,表面扩散和Knudsen扩散的影响越大,反之则越小;人工裂缝的性质包括条数、开度、半长以及间距,主要影响压裂水平井生产早期,随着人工裂缝参数值的增加,压裂水平井产能增加,累产气量也越大.其次,页岩气藏压裂诱导缝和天然裂缝的发育程度对页岩气藏的产能有很大的影响,水平井周围只有人工裂缝,周围天然裂缝不开启或不发育时,页岩气藏的水平井的产能较低.      

气藏压裂水平井非稳态流动分析

应用Green函数和Newman积原理,建立了压裂水平井与气藏耦合的非稳态流动模型,并给出求解方法。所建立的模型可应用多种约束条件,考虑了井筒压降和加速度的影响,并适用于各向异性气藏。实例计算表明：水平井在早期非稳态阶段的流量比拟稳态阶段大得多。早期非稳态流动阶段,各条裂缝之间未发生干扰,流量的大小与裂缝条数成正比例关...     

致密砂岩垂直裂缝井CO2混相驱试井模型研究

向地层注入CO2可以有效地提高致密砂岩原油采收率,常规的试井解释数学模型不能满足致密砂岩无限导流垂直裂缝井CO2混相驱试井解释的需求．因此,基于渗流力学基本原理建立考虑应力敏感影响的无限导流垂直裂缝井CO2混相驱试井解释数学模型,利用Laplace变换、摄动变换和Stehfest数值反演的方法进行求解,编程绘制典型特征曲线并进行敏感性分析．研究表明：该模型典型特征曲线共划分为八个流动阶段．由于应力敏感效应的影响,径向流阶段内、外区压力导数曲线不再是0.5和0.5M13水平线,而是呈“上翘”的曲线,并且应力敏感系数越大,曲线“上翘”越明显;混相区压力导数曲线符合幂律型变化指数规律且高于(1-n)/(3-n)斜率直线;内区、混相区半径和M12的变化都会使得外区压力曲线升高;通过该模型可以有效地对致密砂岩压裂井CO2混相驱试井资料进行解释．     

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

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

缝洞型油藏三维离散缝洞数值试井模型

缝洞型碳酸盐岩油藏发育着大尺度的溶洞和裂缝,非均质性极强,缝洞型碳酸盐岩油藏问题的研究成为了世界级难题之一.由于大尺度溶洞和裂缝对储层的流体流动起主导作用,因此,基于连续介质理论的双重介质或三重介质模型已不适合其中流体流动的描述.根据大型缝洞分布地质特征,探索性地提出了一种板块组合的复合架构离散缝洞模型描述该类油藏中的流体流动,将三维空间大裂缝用板块描述,溶洞用高渗透率和高孔隙度不规则多面体团块描述.将裂缝面用二维三角形单元离散, 溶洞和基质用三维四面体离散, 利用三维混合单元有限元法对建立的不定常渗流模型进行求解,得到了三维渗流条件下的试井理论曲线及压力场分布.通过对试井理论曲线特征的分析, 获得了各敏感参数对试井曲线的影响规律.通过对1口井的实际测试资料解释结果的分析,并与实际地震反射资料及生产实际资料的对比,发现本文所建立的模型可以较好地反映裂缝和溶洞的地质动态状况,并与实际生产状况具有较好的一致性.这一结果说明了所建模型的正确性以及测试资料分析结果的可靠性.      

A new computation model of fractured horizontal well coupling with reservoir

Based on Green's functions and Newman's product principle, pressure drop formula was derived for considering simultaneous production of fractures and horizontal wellbore in unsteady state. A reservoir/fractured horizontal well coupling model is developed for finite conductivity condition that can be solved by the combination of quasi‐Newton method and PSO (Particle Swarm Optimization) algorithm. The solution of a practical example shows that lots of factors can affect the productivity of the fractured horizontal well. The number of the fractures has an optimizing range. Different fractures have different flow rates and the production in the central fracture is the lowest. The distribution of production rate in the wellbore has a wave‐like shape due to the influence of fractures. The closer fractures are the lower flow rate in well segments. Considering the horizontal wellbore production, the production rates in the fractures are asymmetric and the corresponding pressure drop curve is smoother. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical investigation of dual-porosity model with transient transfer function based on discrete-fracture model

Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed.Bottom hole pressure response curves and the pressure field are obtained by solving the model equations with the finite-element method. By analyzing bottom hole pressure curves and the fluid flow in the pressure field, seven flow stages can be recognized on the curves. An upscaling method is developed to compare with the dual-porosity model(DPM). The comparisons results show that the DPM overestimates the inter-porosity coefficient λ and the storage factor ω. The analysis results show that fracture conductivity plays a leading role in the fluid flow. Matrix permeability influences the beginning time of flow from the matrix to fractures. Fractures density is another important parameter controlling the flow. The fracture linear flow is hidden under the large fracture density.The pressure propagation is slower in the direction of larger fracture density.     

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.     

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

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

On Vaporizing Water Flow in Hot Sub-Vertical Rock Fractures

Water injection into unsaturated fractured rock at above-boiling temperatures gives rise to complex fluid flow and heat transfer processes. Examples include water injection into depleted vapor-dominated geothermal reservoirs, and emplacement of heat-generating nuclear wastes in unsaturated fractured rock. We conceptualize fractures as two-dimensional heterogeneous porous media, and use geostatistical techniques to generate synthetic permeability distributions in the fracture plane. Water flow in hot high-angle fractures is simulated numerically, taking into account the combined action of gravity, capillary, and pressure forces, and conductive heat transfer from the wall rocks which gives rise to strong vaporization. In heterogeneous fractures boiling plumes are found to have dendritic shapes, and to be subject to strong lateral flow effects. Fractures with spatially-averaged homogeneous permeabilities tend to give poor approximations for vaporization behavior and liquid migration patterns. Depending on water flow rates, rock temperature, and fracture permeability, liquid water can migrate considerable distances through fractured rock that is at above-boiling temperatures and be only partially vaporized.     

Transient Pressure Behavior of Reservoirs with Discrete Conductive Faults and Fractures

Fractures and faults are common features of many well-known reservoirs. They create traps, serve as conduits to oil and gas migration, and can behave as barriers or baffles to fluid flow. Naturally fractured reservoirs consist of fractures in igneous, metamorphic, sedimentary rocks (matrix), and formations. In most sedimentary formations both fractures and matrix contribute to flow and storage, but in igneous and metamorphic rocks only fractures contribute to flow and storage, and the matrix has almost zero permeability and porosity. In this study, we present a mesh-free semianalytical solution for pressure transient behavior in a 2D infinite reservoir containing a network of discrete and/or connected finite- and infinite-conductivity fractures. The proposed solution methodology is based on an analytical-element method and thus can be easily extended to incorporate other reservoir features such as sealing or leaky faults, domains with altered petrophysical properties (for example, fluid permeability or reservoir porosity), and complicated reservoir boundaries. It is shown that the pressure behavior of discretely fractured reservoirs is considerably different from the well-known Warren and Root dual-porosity reservoir model behavior. The pressure behavior of discretely fractured reservoirs shows many different flow regimes depending on fracture distribution, its intensity and conductivity. In some cases, they also exhibit a dual-porosity reservoir model behavior.