首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 87 毫秒
1.
气藏压裂水平井非稳态流动分析   总被引:1,自引:0,他引:1  
应用Green函数和Newman积原理,建立了压裂水平井与气藏耦合的非稳态流动模型,并给出求解方法。所建立的模型可应用多种约束条件,考虑了井筒压降和加速度的影响,并适用于各向异性气藏。实例计算表明:水平井在早期非稳态阶段的流量比拟稳态阶段大得多。早期非稳态流动阶段,各条裂缝之间未发生干扰,流量的大小与裂缝条数成正比例关...  相似文献   

2.
高岳  王涛  严子铭  柳占立  庄茁 《力学学报》2022,54(8):2248-2268
本文研究了在页岩气高效开采中钻井完井和水力压裂缝网改造的关键力学问题. 提出了页岩多孔弹性介质的本构、强度和断裂韧性的各向异性模型, 指出了微观均匀假设与微观各向同性假设在页岩多孔弹性本构中的简化应用, 给出了横观各向同性多孔弹性岩石材料常数的简化测量方案, 讨论了基于修正的能量最大释放率下裂纹扩展的弱面模型, 提出了裂纹扩展禁止区现象. 阐述了钻井完井过程中的多孔弹性介质井壁稳定性和剪切破坏的时间效应, 提供了多种破坏模式下井壁许可压力范围的显式表达式, 并与传统广义胡克定律得到结果对比. 给出了水平井水力压裂缝网改造技术, 包含水力压裂的大物模实验技术、水力压裂过程中基于扩展有限元与有限体积法的耦合流体/固体/裂缝扩展的数值模拟方法, 并与黏性主导的水力压裂解析解结果对比, 针对性分析了川渝地区的水力压裂施工现场实践应用问题. 建立了基于数据驱动的页岩气采收率预测方法, 指出了机器学习中极限梯度爬升法在小数据集情况下的有效性.   相似文献   

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

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

5.
水平井和水力压裂是页岩气开发中的关键技术。对水力压裂中多簇裂缝同时扩展的物理过程进行了数值模拟。采用扩展有限元法(XFEM)模拟岩石中裂缝沿着任意路径扩展,采用有限体积法(FVM)模拟裂缝中流体的流动,并且考虑井筒中流体流动以及在各簇裂缝间的流量动态分配。通过牛顿迭代对全耦合物理过程进行数值求解,重点研究了初始长度不同的两条裂缝的扩展过程,证明较大的射孔摩阻能促进两条裂缝的同时扩展。并通过算例证明了本方法的精度和有效性。  相似文献   

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

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

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

9.
不仅考虑低渗透油藏具有启动压力梯度的渗流特征,还考虑应力敏感地层中介质的变形;发展了Cinco-Ley H.提出的有限导流垂直裂缝井双线性流理论,将流体在垂直裂缝与地层中形成的流动划分为两个区域—垂直裂缝中的线性流区域和变形三重介质低渗透油藏中的非线性流区域;由此建立了变形三重介质低渗透油藏有限导流垂直裂缝井的非线性流...  相似文献   

10.
两套节点格林元嵌入式离散裂缝模型数值模拟方法   总被引:1,自引:1,他引:0  
程林松  杜旭林  饶翔  曹仁义  贾品 《力学学报》2022,54(10):2892-2903
对于原始嵌入式离散裂缝模型(EDFM), 在计算包含裂缝单元的基质网格内的压力分布时采用了线性分布假设, 这导致了油藏开发早期对非稳态窜流量的计算精度不足. 因此, 本文提出了一种两套节点格林元法的EDFM数值模拟方法. 两套节点格林元法的核心思想是将压力节点与流量节点区分开, 一套压力节点设置在单元顶点, 另一套流量节点设置在网格边的中点, 满足局部物质守恒、具有二阶精度的同时, 可适用于任意网格类型. 本文将两套节点格林元法与EDFM耦合, 采用了非稳态渗流控制方程的边界积分形式推导了基质网格与裂缝网格之间传质量的新格式, 代替了线性分布假设以提高模拟精度; 此外, 修正后的EDFM能适应任意形态的基质网格剖分, 拓展了原始EDFM仅适用于矩形基质网格、难以考虑复杂油藏边界的局限性. 研究表明: 通过对比商业模拟软件tNavigator? LGR模块与原始EDFM, 验证了本文模型具有较高的早期计算精度; 以复杂油藏边界?缝网?SRV分区模型为例, 通过对比SFEM-COMSOL商业模拟软件, 验证了本文模型处理复杂问题的适应性. 本文研究可用于裂缝性油藏开发动态的精确模拟.   相似文献   

11.
Double-diffusive natural convection within a multilayer anisotropic porous medium is studied numerically and analytically. The domain composed of two horizontal porous layers is subjected to a uniform horizontal heat flux and a vertical mass flux, where only the lower one is thermally anisotropic. Darcy model with classical Boussinesq approximation is used in formulating the mathematical model. The effect of thermal anisotropy and the relative width of the two layers on the flow and transfers is illustrated with characterising the transitions from the diffusive to the convective solution. Results were well compared with respect to a developed analytical approach, based on a parallel flow approximation for thermally anisotropic multilayer media.  相似文献   

12.
The elastic fracture behavior of a plate subjected to uniform stress surrounding two equal cracks inclined at an angle is investigated. The orientation of the crack plane with applied stress can be varied. Among the cases are: (1) two cracks inclined symmetrically with respect to the vertical and horizontal applied stress, and (2) one crack is horizontal while the other is inclined to the vertical applied stress. The strain energy density criterion is used for determining the combined crack and load arrangement that correspond to the lowest critical load at global instability. The direction of crack initiation is also determined. Quantitative results pertaining to the fracture characteristics are given in graphical forms.  相似文献   

13.
The elastic fracture behavior of a plate subjected to uniform stress surrounding two equal cracks inclined at an angle is investigated. The orientation of the crack plane with applied stress can be varied. Among the cases are: (1) two cracks inclined symmetrically with respect to the vertical and horizontal applied stress, and (2) one crack is horizontal while the other is inclined to the vertical applied stress. The strain energy density criterion is used for determining the combined crack and load arrangement that correspond to the lowest critical load at global instability. The direction of crack initiation is also determined. Quantitative results pertaining to the fracture characteristics are given in graphical forms.  相似文献   

14.
The probability distributions of sand particles' lift-off and incident velocities in a wind-blown sand flux play very important roles in the simulation of the wind-blown sand movement. In this paper, the vertical and the horizontal speeds of sand particles located at 1.0 mm above a sand-bed in a wind-blown sand flux are observed with the aid of Phase Doppler Anemometry (PDA) in a wind tunnel. Based on the experimental data, the probability distributions of not only the vertical lift-off speed but also the lift-off velocity as well as its horizontal component and the incident velocity as well as its vertical and horizontal components can be obtained by the equal distance histogram method. It is found, according to the results of the χ 2-test for these probability distributions, that the probability density functions (pdf's) of the sand particles' lift-off and incident velocities as well as their vertical components are described by the Gamma density function with different peak values and shapes and the downwind incident and lift-off horizontal speeds, respectively, can be described by the lognormal and the Gamma density functions. These pdf's depend on not only the sand particle diameter but also the wind speed.  相似文献   

15.
This paper presents comprehensive correlation equations of the local Nusselt numbers and surface shear stresses for laminar forced convection, natural convection, and mixed convection on vertical and horizontal flat plates which are maintained with uniform wall temperature or uniform surface heat flux. The correlation for pure forced convection and pure natural convection are very accurate for any Prandtl number between 0.001 and infinity. The correlations for mixed convection coincide very well with the numerical results over the entire regimes of mixed convection intensity and Prandtl number for the eight cases of the two plates with distinct thermal boundary conditions and buoyancy-assisting and-opposing flow conditions.  相似文献   

16.
An analytical model for studying double-diffusive natural convection within a multilayer anisotropic porous medium is developed and validated with respect to a direct numerical silmulation. The studied domain is composed of two horizontal porous layers where the lower one is thermally anisotropic and is submitted to a uniform horizontal heat flux and a vertical mass flux. The assumption of parallel flow is validated and the effect of anisotropy on dynamic transitions is investigated. To cite this article: R. Bennacer et al., C. R. Mecanique 332 (2004).  相似文献   

17.
A new approach is presented for the pressure calculation around a horizontal well or a partially penetrating vertical well in a finite height reservoir. The new method matches the asymptotic solutions close to- and far from the well. To this end, the reservoir is divided into two distinct volumes with ellipsoidal and elliptically cylindrical symmetry, respectively. Comparison with other techniques shows that the new expressions yield excellent results. The method is also applicable to reservoirs containing one elliptical, vertical or horizontal fracture, and can be used with account of a zone of altered permeability around the well or fracture due to formation damage or stimulation.  相似文献   

18.
The probability distributions of sand particles' lift-off and incident velocities in a wind-blown sand flux play very important roles in the simulation of the wind-blown sand movement. In this paper, the vertical and the horizontal speeds of sand particles located at 1.0 mm above a sand-bed in a wind-blown sand flux are observed with the aid of Phase Doppler Anemometry (PDA) in a wind tunnel. Based on the experimental data, the probability distributions of not only the vertical lift-off speed but also the lift-off velocity as well as its horizontal component and the incident velocity as well as its vertical and horizontal components can be obtained by the equal distance histogram method. It is found, according to the results of the χ2-test for these probability distributions, that the probability density functions (pdf's) of the sand particles' lift-off and incident velocities as well as their vertical components are described by the Gamma density function with different peak values and shapes and the downwind incident and lift-off horizontal speeds, respectively, can be described by the lognormal and the Gamma density functions. These pdf's depend on not only the sand particle diameter but also the wind speed. The project supported by the National Natural Science Foundation of China (10532040) and the Hundred Talents Project, the Knowledge Innovation Project of Chinese Academy of Sciences (KZCX2-304). The English text was polished by Keren Wang.  相似文献   

19.
In this paper we study the pressure drop in a hydraulic fracture after shut-in of a water injection well. The pressure transient behavior depends on fracture closure, lateral stress, rock elasticity and fracture fluid leak-off. Under the assumption that horizontal cross-sections of a vertical fracture do not depend on the vertical variable, we formulate a mathematical model which allows for determination of both pore pressure and elastic rock displacements jointly with the fracture aperture and fracture fluid pressure. An analytical consideration is performed for the case of an ideal very long fracture with the same aperture along its full length. In the general case, fracture closure is analyzed numerically.  相似文献   

20.
Data recorded during a well test is interpreted for formation parameters, such as permeability, by comparing the measured pressure transient with that predicted by a mathematical model of the system. In single-phase homogeneous situations, this model is based on the linear diffusion equation. Despite the simplicity of this equation, it is often difficult or laborious to construct exact solutions, in some cases because the relevant problem is of mixed boundary-value character. Here we describe an approximate calculational method which, with little effort, gives good results in a variety of well test problems of this type.We first prove a rather general perturbation theorem, and then apply it in three illustrative cases. The first example provides a test of our basic result for the case of a fully penetrating vertical fracture, for which an exact (though nonelementary) solution is known. The second example is that of partially penetrating or horizontal wells; it is shown that our general result can provide a mathematical basis for the so-called pressure-averaging technique in which the pressure-flux relationship is approximated by assuming that the flux is uniform along the well and then computing the spatial average of the wellbore pressure. Our third example is a new result for the steady-state pressure drop due to flow into a small circular hole on the surface of an impermeable cylinder. This example has relevance for the Schlumberger RFT tool, which withdraws fluid from the formation via a circular probe which penetrates the mudcake surrounding the borehole wall. Numerical results for the shape factor which represents the effect of the borehole curvature are provided.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号