首页期刊简介编委会编辑部过刊浏览征稿简则期刊信息下载中心欢迎订阅 English

引 en

分形离散裂缝页岩气藏多级压裂水平井数值试井模型

吴明录
，丁明才

1. 中国石油大学(华东)石油工程学院,山东青岛 266580；

A numerical well testing interpretation model for multiple fractured horizontal well in fractured shale gas reservoir based on fractal discrete fracture network

WU Minglu
， DING Mingcai

1. School of Petroleum Engineering in China University of Petroleum(East China),Qingdao 266580 , China；

作者简介:

吴明录(1978-),男,副教授,博士,硕士生导师,研究方向为油气田开发工程。wuminglu@upc.edu.cn。

中图分类号:TE377

文献标识码:A

文章编号:1673-5005(2020)03-0098-07

DOI:10.3969/j.issn.1673-5005.2020.03.011

全文
评论
参考文献
出版信息

参考文献 1

WARREN J E,ROOT P J.The behavior of naturally fractured reservoirs[J].Society of Petroleum Engineers Journal,1963,3(3):245-255.

查找原文

参考文献 2

LAUBACH S E.Fracture patterns in low-permeability-sandstone gas reservoir rocks in the rocky mountain re-gion[R].SPE 21853-MS,1991.

查找原文

参考文献 3

PRIEST S D.Discontinuity analysis for rock engineering [M].Dordrecht:Springer Science & Business Media,1993:395-410.

查找原文

参考文献 4

VILLAESCUSA E.Statistical modelling of rock jointing:Conference on probabilistic methods in geotechnical engi-neering[C].Rotterdam:A A Balkema Publishers,1993.

查找原文

参考文献 5

TAMAGAWA T,MATSUURA T,ANRAKU T,et al.Construction of fracture network model using static and dynamic data[R].SPE 77741,2002.

查找原文

参考文献 6

DARCEL C,BOUR O,DAVY P,et al.Connectivity properties of two-dimensional fracture networks with sto-chastic fractal correlation [J].Water Resources Re-search,2003,39(10):1-13.

查找原文

参考文献 7

DAVY P,SORNETTE A,SORNETTE D.Some conse-quences of a proposed fractal nature of continental faulting [J].Nature,1990,348(6296):56-58.

查找原文

参考文献 8

KIM T H,SCHECHTER D S.Estimation of fracture po-rosity of naturally fractured reservoirs with no matrix po-rosity using fractal discrete fracture networks [J].SPE Reservoir Evaluation & Engineering,2009,12(2):232-242.

查找原文

参考文献 9

吴明录,丁明才,姚军,等.分形离散裂缝数值试井解释模型[J].东北石油大学学报,2016,40(6):114-120.WU Minglu,DING Mingcai,YAO Jun,et al.Numerical well-testing interpretation model based on fractal discrete fracture[J].Journal of Northeast Petroleum University,2016,40(6):114-120.

查找原文

参考文献 10

OZKAN E,RAGHAVAN R.New solutions for well-test analysis problems:part 1—analytical considerations [J].SPE Formation Evaluation,1991,6(3):359-368.

查找原文

参考文献 11

ZERZAR A,BETTAM Y.Interpretation of multiple hy-draulically fractured horizontal wells in closed systems [R].SPE 84888-MS,2004.

查找原文

参考文献 12

BETTAM Y,ZERZAR A,TIAB D.Interpretation of multi-hydraulically fractured horizontal wells in natural-ly fractured reservoirs[R].SPE 95535-MS,2005.

查找原文

参考文献 13

ESCOBAR F H,MARTÍNEZ J A,MONTEALEGRE-MADERO M.Pressure transient analysis for a reservoir with a finite-conductivity fault[J].Ciencia,TecnologÍay Futuro,2013,5(2):5-17.

查找原文

参考文献 14

FAN D,YAO J,SUN H,et al.A composite model of hydraulic fractured horizontal well with stimulated reser-voir volume in tight oil & gas reservoir[J].Journal of Natural Gas Science and Engineering,2015,24:115-123.

查找原文

参考文献 15

朱光普,姚军,樊冬艳,等.页岩气藏压裂水平井试井分析[J].力学学报,2015,47(6):945-954.ZHU Guangpu,YAO Jun,FAN Dongyan,et al.Pres-sure transient analysis of fractured horizontal well in shale Gas reservoir[J].Chinese Journal of Theoretical and Applied Mechanics,2015,47(6):945-954.

查找原文

参考文献 16

樊冬艳,姚军,孙海,等.页岩气藏分段压裂水平井不稳定渗流模型[J].中国石油大学学报(自然科学版),2014,38(5):116-123.FAN Dongyan,YAO Jun,SUN Hai,et al.Transient flow model of stage-fractured horizontal wells in shale gas reservoirs[J].Journal of China University of Petro-leum(Edition of Natural Science),2014,38(5):116-123.

查找原文

参考文献 17

姚军,孙海,樊冬艳,等.页岩气藏运移机制及数值模拟[J].中国石油大学学报(自然科学版),2013,37(1):91-98.YAO Jun,SUN Hai,FAN Dongyan,et al.Transport mechanisms and numerical simulation of shale gas reser-voirs[J].Journal of China University of Petroleum(Edi-tion of Natural Science),2013,37(1):91-98.

查找原文

参考文献 18

MOINFAR A.Development of an efficient embedded discrete fracture model for 3D compositional reservoir simulation in fractured reservoirs[D].Texas:The Uni-versity of Texas at Austin,2013.

查找原文

参考文献 19

SHAKIBA M.Modeling and simulation of fluid flow in naturally and hydraulically fractured reservoirs using em-bedded discrete fracture model [ D].Texas:The Uni-versity of Texas at Austin,2014.

查找原文

目录contents

摘要Abstract
关键词Keywords
1 分形离散裂缝随机建模方法
2 数学模型的建立
2.1 基本假设
2.2 基质系统渗流数学模型
2.3 人工裂缝渗流数学模型
3 数学模型的求解与验证
3.1 模型的求解
3.2 一致性验证
4 试井曲线特征及参数敏感性
4.1 试井曲线特征
4.2 参数敏感性
4.2.1 裂缝密度 α
4.2.2 裂缝分形维数D_c 和D_l
4.2.3 裂缝主发育方向 θ₀
4.2.4 Fisher
4.2.5 Langmuir吸附体积V_L 和吸附压力p_L
5 应用实例
6 结论
参考文献

摘要

基于分形理论,考虑天然裂缝的长度、方向、中心位置、开度等分布规律建立分形离散裂缝网络模型,并在此基础上建立多级压裂水平井试井模型。采用有限差分和嵌入式离散裂缝相结合的方法求解该模型,得到水平井的压力响应及其影响规律。结果表明:裂缝密度、分形维数、裂缝主发育方向和 Fisher 常数等参数能够较好地控制试井曲线形态; 裂缝密度、中心维数和长度维数决定裂缝的数量,裂缝密度和中心维数越大(或长度维数越小),裂缝总数越多,水平井的双重介质特征越显著,即中期拟压力导数下凹越明显;裂缝主发育方向和 Fisher 常数决定裂缝的走向及其一致性程度,Fisher 常数越大,裂缝的方向越一致,渗流时有效渗流面积越小,则晚期拟压力导数曲线越高;Langmuir 吸附体积和吸附压力对试井曲线也具有一定的影响,二者越大,中期和晚期拟压力和拟压力导数曲线整体越低。

Abstract

A discrete fracture network model was established based on the fractal theory, considering the length, direction, center position and width of natured fractures, and then a numerical well testing model for multiple fractured horizontal well in fractured shale gas reservoir was obtained. The model was solved by a finite difference method combined with embedded discrete fractures. The pressure response and its influence characteristics of horizontal wells were studied. The analysis re- sults show that the pressure response curve can be well controlled by parameters such as the fracture density, the fractal di- mension, the main propagation direction of fractures and the Fisher constant. The total number of fractures is dependent on the fracture density and the fracture’s center and length dimensions. The larger the fracture density and the fracture center di- mension (or the smaller the fracture length dimension), the greater the total number of fractures. As a consequence, the du- al-medium characteristics of the horizontal wells can be more visible, or in other words, the pressure derivative curve can de- press more obviously during the intermediate period. The analysis also indicates that the orientation and its consistency of the fractures are dependent on the propagation direction of the main fracture and the Fisher constant. The pressure derivative curve can be higher with larger Fisher constant during the late flow period due to smaller effective flow area caused by more consistent fracture orientation. Furthermore, The Langmuir adsorption volume and adsorption pressure also have certain influ- ence on the pressure response. The pressure and its derivative curves can be lower together with larger Langmuir adsorption volume and adsorption pressureduring the intermediate and late flow periods.

关键词

水平井；数值试井；多级压裂；页岩气；离散裂缝；分形

Keywords

horizontal well ； numerical well testing ； multiple fractured ； shale gas ； discrete fracture ； fractal

页岩气储层多有裂缝发育,早期的研究通常基于连续性介质假设^[1],Laubach等^[2] 认为储层中的裂缝应该是离散的、不规则的、经常成簇出现的,并提出了离散裂缝网络(DFN) ^[3-5] 的概念,但他们所提出的采用随机建模生成裂缝的方法具有极大的不确定性。 Darcel等^[6]、Davy等^[7]、Kim等^[8] 提出了一种应用分形理论随机生成离散裂缝网络的方法, 采用这种方法生成的裂缝称为分形离散裂缝模型(FDFN),研究表明分形离散裂缝模型更加接近裂缝性油藏的实际情况,而且分形参数对裂缝的整体分布具有很好的控制作用,但并未对分形裂离散裂缝模型的动态特征进行研究。吴明录等^[9] 将分形离散裂缝模型与试井解释相结合,提出了分形离散裂缝数值试井模型,研究了其不稳定压力响应特征, 研究发现分形参数能够很好地控制试井曲线的形态,研究结果为分形离散裂缝模型应用于裂缝性油藏建模和生产动态分析具有重要意义,但该项研究只针对普通砂岩油藏垂直井情形。页岩气藏渗流机制复杂,而且通常采用水平井多级压裂的方式进行开采。目前的页岩气藏压裂水平井试井理论,以采用点源解、线源解和格林函数等方法进行解析求解的方法为主^[10-13],该方法在建立模型过程常常忽略页岩气藏储层的各向异性、非均质性以及复杂流动机制。以离散裂缝模型为基础的数值试井方法^[14-18]由于求解压裂水平井压力动态时需要根据裂缝分布划分大量非结构化网格,计算量大,大大影响了现场实际应用。笔者考虑页岩气吸附解吸、黏性流、基质与裂缝间的窜流等流动机制,将分形离散裂缝网络应用于页岩气藏多级压裂水平井,建立数值试井模型,并采用嵌入式离散裂缝方法进行求解, 研究水平井的拟压力响应特征及分形参数对其影响规律。
1 分形离散裂缝随机建模方法
目前,分形离散裂缝网络建模的方法主要有两种:第一种采用随机建模方法以裂缝中心分布为基础;第二种采用随机建模方法以裂缝长度分布为基础^[6]。前者具有诸多优点:首先裂缝中心分布易于随机模拟生成;其次在得到裂缝中心分布的情况下, 依据裂缝长度分布、裂缝方向特征便于得到裂缝的总体分布。因此本文中采用基于裂缝中心分布的建模方法。基于裂缝中心分布的离散裂缝网络模型的建模过程可以参照参考文献^[6]和^[9]。
2 数学模型的建立
2.1 基本假设
假设:①储层均质、水平、等厚,且基质渗透率各向同性;②储层中的气体流动符合达西定律;③储层中的裂缝是分形离散的,且所有裂缝垂直贯穿整个储层;④自由气储存在裂缝和基质孔隙中,吸附气吸附满足Langmuir等温吸附定律; ⑤忽略重力和毛管力,开井生产前储层各处压力相等,均为原始地层压力; ⑥水平井为射孔完井,天然气只由人工裂缝进入井筒,定产量生产。
2.2 基质系统渗流数学模型
考虑页岩基质颗粒表面吸附解吸,基质系统渗流控制方程为

\frac{\partial}{\partial t} (φ_{m} C_{p} ρ_{m_{B}} + (1 - φ_{m}) q_{m, a h s}) = \nabla \cdot (ρ_{m g} \frac{k_{m}}{μ_{r}} \nabla p)

(1)

其中
$q_{m, als} = \frac{ρ_{s} M_{g} Z_{g}}{Z_{B} V_{std} {(p_{L} + p)}^{2}}, ρ_{m g} = \frac{M_{e p}}{Z_{B} R T}$
式中,ρ_s为岩石密度,g/cm ³;V_L 为Langmuir吸附体积,cm ³/g;p为压力,10 ^-1MPa;p_L为Langmuir吸附压力,10 ^-1MPa;Vstd为标准状况(0℃,0.1 MPa)下天然气的摩尔体积,cm³/mol;Z_g为气体压缩因子;μ_g为气体黏度,mPa·s;φ_m为基质孔隙度;C_g为气体压缩系数,cm ³/(10^-1 MPa);ρ_mg为基质中页岩气的密度,g/cm³;M_g为天然气摩尔质量,g/mol;R为热力学常数,R=8.314 J/(mol·K);T为储层温度,K。
定义拟压力为

ψ = \int_{p_{0}}^{p} \frac{2 p}{Z_{g} μ_{g}} d p

(2)

将拟压力带入式(1),可得到基质系统渗流控制方程的拟压力形式为

γ_{m} \frac{\partial ψ_{m}}{\partial t} = \nabla \cdot (β k_{m} \nabla ψ_{m})

(3)

其中
$β = \frac{M_{g}}{R T}, γ_{m} = β φ_{m} C_{g} μ_{g} + \frac{ρ_{s} M_{g} Z_{g} μ_{g} (1 - φ_{m}) V_{L} p_{L}}{V_{s t d}}$
初始条件为

ψ_{f i} = ψ_{h i} = ψ_{m i}

(4)

封闭外边界为

{\frac{\partial ψ_{m}}{\partial n} |}_{Γ} = 0

(5)

式中, $Γ$ 为外边界; n为边界 $Γ$ ₁ 的外法向方向的单位向量。
2.3 人工裂缝渗流数学模型
裂缝中渗流控制方程为

\frac{\partial}{\partial t} (ρ_{h g} φ_{m} C_{g}) = \nabla \cdot (ρ_{h g} \frac{k_{h}}{μ_{g}} \nabla p_{h})

(6)

其中
$ρ_{h g} = \frac{M_{g} p_{h}}{Z_{g} R T}$
将拟压力代入式(6),得到裂缝系统渗流控制方程的拟压力形式为

γ_{h} \frac{\partial ψ_{h}}{\partial t} = \nabla \cdot (β k_{h} \nabla ψ_{h})

(7)

其中
$γ_{h} = β φ_{m} C_{g} μ_{g}$
将裂缝嵌入到基质系统中,裂缝与基质网格单元间窜流量的表达式为

q_{m - h} = T_{m - h} β (ψ_{m} - ψ_{h})

(8)

式中,T_m-h为基质与人工裂缝间的窜流系数^[9]。
根据假设,天然气只通过人工裂缝进入井筒,每条人工裂缝流入井筒的天然气量为

q_{h - w} = I_{W} β μ_{g} (ψ_{h} - ψ_{w})

(9)

式中,I_W 为人工裂缝到井筒的井指数^[20]。
当水平井定产量生产时,总产量为井筒中压缩流体膨胀量和由人工裂缝流向井筒的流量的总和,有

Q_{g} = \sum_{i = 1}^{N_{h}} I_{W} β μ_{g} (ψ_{h i} - ψ_{w}) - C \frac{β μ_{g}}{2} \frac{\partial ψ_{w}}{\partial t}

(10)

式中,C为井筒储存系数,cm ³/(10 ^-1MPa)。
3 数学模型的求解与验证
3.1 模型的求解
应用嵌入式离散裂缝方法对以上模型进行求解。根据Shakiba ^[19]的研究,在将离散裂缝嵌入到基质中时,网格m与网格n之间的流动系数计算式为

T_{m - n} = \frac{T_{m - f} T_{n - f}}{T_{m - f} + T_{n - f}}

(11)

式中,T_m-f和T_n-f为基质网格m、n与所嵌入的裂缝之间的流动系数,可由基质与裂缝的渗透率计算得到。
3.2 一致性验证
为验证裂缝参数对解的一致性,取相同的分形参数:裂缝中心维数D_c=1.8,裂缝主发育方向 θ₀=π/4,Fisher常数K=25,裂缝密度 α=2.5,裂缝长度维数D_l=1.3。通过随机模拟生成4 个分形离散裂缝模型,如图1 所示。模拟区域为1600 m×1600 m, 水平井长度为1000 m,等间隔分6 段压裂,每段裂缝半长为200 m。求解该模型,得到水平井拟压力响应曲线如图2 所示。
图1 分形离散裂缝模型
Fig.1 Fractal discrete fracture models
图2 不同裂缝网络模型的拟压力响应曲线
Fig.2 Pseudo-pressure response curves for different fractal discrete fracture models
由图2 可见,在相同的分形参数条件下,水平井的拟压力和拟压力导数曲线基本重合,但由于这4 个模型由随机建模生成,具有一定的差异性(图1), 水平井在生产过程中形成的压降漏斗不尽相同,因此其试井曲线也存在细微差别。正如弹性储容比、窜流系数等参数能够控制连续双重孔隙介质模型中试井曲线的形态一样,通过对以上4 个模型的计算表明,虽然局部裂缝分布各不相同,但只要总体分形参数相同,则由不同的离散裂缝模型得到的压力响应曲线基本相同,也就是说分形参数对压力响应曲线的形态具有较好的控制作用。
4 试井曲线特征及参数敏感性
基本参数为 θ₀=π/4、K=25、α=2.5、D_l=1.3、 D_c=1.8。
4.1 试井曲线特征
假设6 条人工裂缝,且沿水平井筒呈对称分布(图1),得到该模型的拟压力响应曲线如图3 所示。
由图3 可见,水平井的压力响应包括6 个阶段: ①为井筒储存阶段,压力及其导数为斜率为1 的直线;②为早期(第一)过渡流动阶段;③为线性流阶段,天然裂缝间的页岩气通过基质线性流流向人工裂缝,压力和压力导数均为斜率为0.5 的直线;④为中期(第二)过渡流阶段;在流动阶段⑤,人工裂缝周围基质中的流体开始沿着天然裂缝向人工裂缝流动,并通过人工裂缝流向水平井。由于非均匀分布的天然裂缝对渗流场造成的影响,流线和等势线受裂缝分布影响较大,该阶段不是严格的径向流(流线指向半径方向,等势线为近似同心圆),因此压力导数曲线也不是水平线;⑥为边界反映阶段,反映外边界对压力响应的影响,当外边界封闭时,压力和压力导数曲线为斜率为1 的直线。
图3 分形离散裂缝页岩气藏多段压裂水平井的压力响应曲线
Fig.3 Pseudo-pressure response curves for multiple fractured horizontal well in fractal discrete fracture shale gas reservoir
4.2 参数敏感性
4.2.1 裂缝密度 α
依次取 α=0.5、2.5、4.5,分别建立离散裂缝网络多段压裂水平井模型,如图4 所示。采用本文中方法对数学模型进行求解,得到水平井的试井曲线如图5 所示。
由图4 可见,裂缝密度 α 主要影响裂缝的数目, α 越大,裂缝数目越多。由图5 可见,裂缝密度 α 主要影响流动阶段⑤, α 越大,压力导数曲线的“下凹” 越深。在本文中,基质除了作为储容空间外,还具有一定的导流能力(在经典的Warren-Root双重孔隙介质模型中,基质主要视为流体的储容空间,通常不考虑其内部的渗流)。而且在裂缝较少时,基质仍然是主要的渗流介质,而裂缝作为连通基质的高导流通道;当裂缝较多时,裂缝才成为主要的渗流介质。因此 α 越大,裂缝数目越多时,裂缝所占的总系统的弹性孔隙体积越大,双重介质的特征也就越明显,表现在试井曲线上,即是压力导数曲线中期段 “下凹”越深。
图4 不同裂缝密度的分形离散裂缝网络模型
Fig.4 Fractal discrete fracture network models with different fracture dimension density
图5 裂缝密度对拟压力响应曲线的影响
Fig.5 Effect of fracture distribution density on pseudo-pressure response curves
4.2.2 裂缝分形维数D_c 和D_l
由式(1)可见,一定长度范围内的裂缝数目与 α、D_c、D_l有关: α 和D_c 越大,D_l越小,裂缝数目越多。因此D_c 与 α 对试井曲线的影响规律大致相同; D_l 与 α 对试井曲线的影响规律大致相反(图5、6)。
4.2.3 裂缝主发育方向 θ₀
依次取θ₀=0、π/4、π/2,分别建立离散裂缝网络多段压裂水平井模型,如图7 所示。采用本文中方法对数学模型进行求解,得到水平井的试井曲线如图8 所示。
由图8 可见,裂缝主发育方向 θ₀ 对流动阶段⑤ 和⑥ 均有一定程度的影响。在阶段⑥ 之前,流体以垂直于人工裂缝方向(沿X方向)的流动为主,当 θ₀=0 时,由于裂缝的沟通作用储层在水平方向上的渗透性最强,其对应的压力导数曲线最低;而在边界流阶段(⑥),流体以垂直于水平井的流动(沿Y方向)为主,当 θ₀=π/2 时,储层在垂直方向上的渗透性最强,其对应的压力导数曲线最低。
图6 裂缝中心维数和长度维数对拟压力响应曲线的影响
Fig.6 Effect of fracture center dimension and length dimension on pseudo-pressure response curves
图7 不同地应力主方向的分形离散裂缝网络
Fig.7 Fractal discrete fracture network with different formation principal stress direction
图8 裂缝主发育方向对拟压力响应曲线的影响Fig.8 Effect of formation principal stress direction on pseudo-pressure response curves
4.2.4 Fisher
常数K Fisher常数主要用于描述裂缝方向与裂缝发育主方向的偏离程度,K越大,裂缝总体发育方向的一致性越强。依次取K=5、15、45,分别建立离散裂缝网络多段压裂水平井模型,如图9 所示。采用本文中方法对数学模型进行求解,得到水平井的试井曲线如图10 所示。
由图10 可见,Fisher常数K主要影响边界流阶段(⑥),且K越大,压力导数曲线的晚期段越高。这是因为K越大,裂缝总体发育方向一致性越强(图9(c)),这是裂缝的沟通作用使有效渗流面积变小的缘故。
图9 不同Fisher常数的分形离散裂缝网络模型
Fig.9 Fractal discrete fracture network models with different Fisher constant
图10 Fisher常数对拟压力响应曲线的影响
Fig.10 Effect of Fisher constant on pseudo-pressure response curves
4.2.5 Langmuir吸附体积V_L 和吸附压力p_L
依次取V_L 为0.0064、0.0224、0.0384 m ³/kg(p_L=1 MPa)及p_L 为1、3、10 MPa(V_L=0.0064 m ³/kg)对数学模型进行求解,得到水平井的试井曲线, 结果如图11 所示。
V_L 和p_L 反映基质颗粒表面对页岩气的吸附解吸能力。 V_L 越大,给定温度条件下单位质量基质饱和时,吸附的页岩气体积越大;p_L 越大,水平井生产过程中地层压力降低时,页岩气从基质颗粒表面解吸扩散到基质孔隙中的时间越早。由图11 可见,V_L 和p_L 主要影响阶段⑤之后的拟压力和拟压力导数曲线。 V_L 和p_L 越大,单位质量基质颗粒表面吸附的页岩气体积越多,当压力降低到吸附压力以下时, 就会有越多的吸附气从基质颗粒表面解吸扩散到基质孔隙中,缓解水平井的压力降落,因此拟压力和拟压力导数曲线越低。
图11 Langmuir吸附体积和吸附压力对拟压力响应曲线的影响
Fig.11 Effect of Langmuir adsorption volume and adsorption pressure on pseudo-pressure response curves
5 应用实例
为了验证本文中模型的可行性,用Ecrin软件模拟了双重孔隙介质气藏中一口多级压裂水平井的拟压力数据,主要参数:裂缝系统的平均渗透率为15×10 ^-3 μm²,气藏的弹性储容比为8×10 ^-4,窜流系数为2×10 ^-6;水平井长度为1200 m;人工裂缝5 条等间距压裂,裂缝半长为75 m。应用本文中方法计算得到水平井的拟压力动态,将其与以上模拟曲线进行拟合解释(图12),得到气藏的分形参数为:α=2.65,D_c=1.64,D_l=1.12,θ₀=36.5°,K=14.25。由图12 可见,由本文中方法得到的理论曲线除早期受网格效应影响导致拟压力导数拟合较差外,总体与模拟曲线变化规律一致,整体拟合效果较好,证明本文模型应用于现场具有可行性。
图12 本文中方法与Ecrin模拟曲线拟合对比
Fig.12 Comparison of proposed method with Ecrin software
6 结论
(1)分形参数可较好地控制水平井的试井曲线形态,证明本文中模型对离散裂缝性页岩气藏试井解释具有一定的适用性。
(2)裂缝密度和裂缝分形维数主要通过影响裂缝总数目对试井曲线产生影响,裂缝密度越大、中心维数越大(或长度分布维数越小),裂缝总数目越大,裂缝的沟通作用越明显,压力导数曲线中期段“下凹”越深,双重介质的特征也就越明显。
(3)Fisher常数K主要影响边界流阶段,其值越大,压力导数曲线的晚期段越高。裂缝主发育方向 θ₀对边界流和之前的流动阶段均有一定程度的影响:在边界流之前,流体以垂直于人工裂缝方向的流动为主,当离散裂缝方向与人工裂缝垂直( θ₀=0) 时,压力导数曲线最低;在边界流阶段,流体以垂直于水平井的流动为主,当离散裂缝方向与水平井垂直(θ₀=π/2)时,压力导数曲线最低。
(4)Langmuir吸附体积和吸附压力越大,基质颗粒表面对页岩气的吸附能力越强,当压力降低到吸附压力以下时,就会有越多的吸附气解吸扩散到基质孔隙中,压力和压力导数曲线越低。
参考文献
- [1] WARREN J E,ROOT P J.The behavior of naturally fractured reservoirs[J].Society of Petroleum Engineers Journal,1963,3(3):245-255.
- [2] LAUBACH S E.Fracture patterns in low-permeability-sandstone gas reservoir rocks in the rocky mountain re-gion[R].SPE 21853-MS,1991.
- [3] PRIEST S D.Discontinuity analysis for rock engineering [M].Dordrecht:Springer Science & Business Media,1993:395-410.
- [4] VILLAESCUSA E.Statistical modelling of rock jointing:Conference on probabilistic methods in geotechnical engi-neering[C].Rotterdam:A A Balkema Publishers,1993.
- [5] TAMAGAWA T,MATSUURA T,ANRAKU T,et al.Construction of fracture network model using static and dynamic data[R].SPE 77741,2002.
- [6] DARCEL C,BOUR O,DAVY P,et al.Connectivity properties of two-dimensional fracture networks with sto-chastic fractal correlation [J].Water Resources Re-search,2003,39(10):1-13.
- [7] DAVY P,SORNETTE A,SORNETTE D.Some conse-quences of a proposed fractal nature of continental faulting [J].Nature,1990,348(6296):56-58.
- [8] KIM T H,SCHECHTER D S.Estimation of fracture po-rosity of naturally fractured reservoirs with no matrix po-rosity using fractal discrete fracture networks [J].SPE Reservoir Evaluation & Engineering,2009,12(2):232-242.
- [9] 吴明录,丁明才,姚军,等.分形离散裂缝数值试井解释模型[J].东北石油大学学报,2016,40(6):114-120.WU Minglu,DING Mingcai,YAO Jun,et al.Numerical well-testing interpretation model based on fractal discrete fracture[J].Journal of Northeast Petroleum University,2016,40(6):114-120.
- [10] OZKAN E,RAGHAVAN R.New solutions for well-test analysis problems:part 1—analytical considerations [J].SPE Formation Evaluation,1991,6(3):359-368.
- [11] ZERZAR A,BETTAM Y.Interpretation of multiple hy-draulically fractured horizontal wells in closed systems [R].SPE 84888-MS,2004.
- [12] BETTAM Y,ZERZAR A,TIAB D.Interpretation of multi-hydraulically fractured horizontal wells in natural-ly fractured reservoirs[R].SPE 95535-MS,2005.
- [13] ESCOBAR F H,MARTÍNEZ J A,MONTEALEGRE-MADERO M.Pressure transient analysis for a reservoir with a finite-conductivity fault[J].Ciencia,TecnologÍay Futuro,2013,5(2):5-17.
- [14] FAN D,YAO J,SUN H,et al.A composite model of hydraulic fractured horizontal well with stimulated reser-voir volume in tight oil & gas reservoir[J].Journal of Natural Gas Science and Engineering,2015,24:115-123.
- [15] 朱光普,姚军,樊冬艳,等.页岩气藏压裂水平井试井分析[J].力学学报,2015,47(6):945-954.ZHU Guangpu,YAO Jun,FAN Dongyan,et al.Pres-sure transient analysis of fractured horizontal well in shale Gas reservoir[J].Chinese Journal of Theoretical and Applied Mechanics,2015,47(6):945-954.
- [16] 樊冬艳,姚军,孙海,等.页岩气藏分段压裂水平井不稳定渗流模型[J].中国石油大学学报(自然科学版),2014,38(5):116-123.FAN Dongyan,YAO Jun,SUN Hai,et al.Transient flow model of stage-fractured horizontal wells in shale gas reservoirs[J].Journal of China University of Petro-leum(Edition of Natural Science),2014,38(5):116-123.
- [17] 姚军,孙海,樊冬艳,等.页岩气藏运移机制及数值模拟[J].中国石油大学学报(自然科学版),2013,37(1):91-98.YAO Jun,SUN Hai,FAN Dongyan,et al.Transport mechanisms and numerical simulation of shale gas reser-voirs[J].Journal of China University of Petroleum(Edi-tion of Natural Science),2013,37(1):91-98.
- [18] MOINFAR A.Development of an efficient embedded discrete fracture model for 3D compositional reservoir simulation in fractured reservoirs[D].Texas:The Uni-versity of Texas at Austin,2013.
- [19] SHAKIBA M.Modeling and simulation of fluid flow in naturally and hydraulically fractured reservoirs using em-bedded discrete fracture model [ D].Texas:The Uni-versity of Texas at Austin,2014.

基本信息

中图分类号: TE377
文献标识码: A
DOI: 10.3969/j.issn.1673-5005.2020.03.011
文章编号: 1673-5005(2020)03-0098-07

基金信息

国家自然科学基金重大项目(51490654)；

引用信息

吴明录,丁明才. 分形离散裂缝页岩气藏多级压裂水平井数值试井模型[J]. 中国石油大学学报(自然科学版),2020,44(3):98-104.

WU Minglu, DING Mingcai. A numerical well testing interpretation model for multiple fractured horizontal well in fractured shale gas reservoir based on fractal discrete fracture network[J]. Journal of China University of Petroleum(Edition of Natural Science),2020,44(3):98-104.

稿件历史

收稿日期: 2019-06-20

参考文献

[1] WARREN J E,ROOT P J.The behavior of naturally fractured reservoirs[J].Society of Petroleum Engineers Journal,1963,3(3):245-255.
[2] LAUBACH S E.Fracture patterns in low-permeability-sandstone gas reservoir rocks in the rocky mountain re-gion[R].SPE 21853-MS,1991.
[3] PRIEST S D.Discontinuity analysis for rock engineering [M].Dordrecht:Springer Science & Business Media,1993:395-410.
[4] VILLAESCUSA E.Statistical modelling of rock jointing:Conference on probabilistic methods in geotechnical engi-neering[C].Rotterdam:A A Balkema Publishers,1993.
[5] TAMAGAWA T,MATSUURA T,ANRAKU T,et al.Construction of fracture network model using static and dynamic data[R].SPE 77741,2002.
[6] DARCEL C,BOUR O,DAVY P,et al.Connectivity properties of two-dimensional fracture networks with sto-chastic fractal correlation [J].Water Resources Re-search,2003,39(10):1-13.
[7] DAVY P,SORNETTE A,SORNETTE D.Some conse-quences of a proposed fractal nature of continental faulting [J].Nature,1990,348(6296):56-58.
[8] KIM T H,SCHECHTER D S.Estimation of fracture po-rosity of naturally fractured reservoirs with no matrix po-rosity using fractal discrete fracture networks [J].SPE Reservoir Evaluation & Engineering,2009,12(2):232-242.
[9] 吴明录,丁明才,姚军,等.分形离散裂缝数值试井解释模型[J].东北石油大学学报,2016,40(6):114-120.WU Minglu,DING Mingcai,YAO Jun,et al.Numerical well-testing interpretation model based on fractal discrete fracture[J].Journal of Northeast Petroleum University,2016,40(6):114-120.
[10] OZKAN E,RAGHAVAN R.New solutions for well-test analysis problems:part 1—analytical considerations [J].SPE Formation Evaluation,1991,6(3):359-368.
[11] ZERZAR A,BETTAM Y.Interpretation of multiple hy-draulically fractured horizontal wells in closed systems [R].SPE 84888-MS,2004.
[12] BETTAM Y,ZERZAR A,TIAB D.Interpretation of multi-hydraulically fractured horizontal wells in natural-ly fractured reservoirs[R].SPE 95535-MS,2005.
[13] ESCOBAR F H,MARTÍNEZ J A,MONTEALEGRE-MADERO M.Pressure transient analysis for a reservoir with a finite-conductivity fault[J].Ciencia,TecnologÍay Futuro,2013,5(2):5-17.
[14] FAN D,YAO J,SUN H,et al.A composite model of hydraulic fractured horizontal well with stimulated reser-voir volume in tight oil & gas reservoir[J].Journal of Natural Gas Science and Engineering,2015,24:115-123.
[15] 朱光普,姚军,樊冬艳,等.页岩气藏压裂水平井试井分析[J].力学学报,2015,47(6):945-954.ZHU Guangpu,YAO Jun,FAN Dongyan,et al.Pres-sure transient analysis of fractured horizontal well in shale Gas reservoir[J].Chinese Journal of Theoretical and Applied Mechanics,2015,47(6):945-954.
[16] 樊冬艳,姚军,孙海,等.页岩气藏分段压裂水平井不稳定渗流模型[J].中国石油大学学报(自然科学版),2014,38(5):116-123.FAN Dongyan,YAO Jun,SUN Hai,et al.Transient flow model of stage-fractured horizontal wells in shale gas reservoirs[J].Journal of China University of Petro-leum(Edition of Natural Science),2014,38(5):116-123.
[17] 姚军,孙海,樊冬艳,等.页岩气藏运移机制及数值模拟[J].中国石油大学学报(自然科学版),2013,37(1):91-98.YAO Jun,SUN Hai,FAN Dongyan,et al.Transport mechanisms and numerical simulation of shale gas reser-voirs[J].Journal of China University of Petroleum(Edi-tion of Natural Science),2013,37(1):91-98.
[18] MOINFAR A.Development of an efficient embedded discrete fracture model for 3D compositional reservoir simulation in fractured reservoirs[D].Texas:The Uni-versity of Texas at Austin,2013.
[19] SHAKIBA M.Modeling and simulation of fluid flow in naturally and hydraulically fractured reservoirs using em-bedded discrete fracture model [ D].Texas:The Uni-versity of Texas at Austin,2014.

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。

分形离散裂缝页岩气藏多级压裂水平井数值试井模型

A numerical well testing interpretation model for multiple fractured horizontal well in fractured shale gas reservoir based on fractal discrete fracture network

摘要

Abstract

关键词

Keywords

1 分形离散裂缝随机建模方法

2 数学模型的建立

2.1 基本假设

2.2 基质系统渗流数学模型

(1)

(2)

(3)

(4)

(5)

2.3 人工裂缝渗流数学模型

(6)

(7)

(8)

(9)

(10)

3 数学模型的求解与验证

3.1 模型的求解

(11)

3.2 一致性验证

4 试井曲线特征及参数敏感性

4.1 试井曲线特征

4.2 参数敏感性

4.2.1 裂缝密度 α

4.2.2 裂缝分形维数Dc 和Dl

4.2.3 裂缝主发育方向 θ0

4.2.4 Fisher

4.2.5 Langmuir吸附体积VL 和吸附压力pL

5 应用实例

6 结论

参考文献

基本信息

基金信息

引用信息

稿件历史

参考文献

4.2.2 裂缝分形维数D_c 和D_l

4.2.3 裂缝主发育方向 θ₀

4.2.5 Langmuir吸附体积V_L 和吸附压力p_L