-
在钻井过程中,液相侵入会导致页岩地层压力上升并引起岩石的压力穿透效应[1],进而造成井壁失稳及套管挤毁等一系列井下事故的发生[2],严重制约了页岩储层的有效开发。国内外学者针对页岩地层的液相侵入进行了大量的研究,指出在渗透压力、化学作用力以及毛细管作用力的驱动下,液相迅速侵入并充满井眼周围的裂缝网络,在压差作用下致使储层基质吸水[3]。随着侵入时间的延长,裂缝表面渗透、吸水带范围不断扩大,直至达到含水平衡并饱和[4]。液相通过微裂缝进入页岩内部后,水化黏土矿物胶结物,岩石力学性质发生明显变化,流固耦合现象显著[5]。
页岩地层是典型的复杂孔隙介质,针对复杂孔隙介质中的渗流问题,Blaskovich等提出了双孔双渗模型,该模型假设基质和裂缝均可渗透,但裂缝是输送流体的主要通道[6]。Noorishad提出了一种同时考虑流体在基岩和裂隙内流动的离散裂隙基岩模型[7]。Meyer等将离散裂缝模型引入页岩储层来描述流体在裂缝内的流动规律[8]。针对渗流引起的流固耦合现象,卢义玉等分析了热流固耦合作用对页岩渗透特性的影响[9]。康永尚等通过对页岩岩心的流-固耦合物理模拟实验分析了岩心渗透率随压力的变化关系[10]。渗流将会引起流固耦合作用,而流固耦合作用同时会影响渗流过程。本文在将人工裂缝及天然裂缝视为两个不同压力系统的基础上,考虑了流固耦合作用对液相渗流的影响,建立了页岩地层复杂孔隙介质液相侵入数学模型,并对计算结果及流固耦合作用对渗流的影响进行了分析,研究结果对于科学认识流固耦合作用下的页岩地层液相侵入机理具有一定的指导意义。
-
为建立页岩地层多重介质液相侵入渗流模型,在建立模型之前做如下假设:
(1)液相在侵入过程中为单相流;
(2)液体在基质及裂缝系统中的渗流均满足达西渗流规律;
(3)只考虑岩石基质及裂缝系统的压缩系数;
(4)液相侵入过程中温度保持恒定;
(5)渗流过程中忽略重力作用;
(6)基质及人工裂缝与井筒之间存在流体交换。
考虑基质与裂缝之间的窜流,并结合运动方程、状态方程[15]代入连续性方程中,可得页岩地层基质及裂缝系统液相渗流的基本方程。
基质渗流方程
$$ \rho \frac{{{k_{\rm{m}}}}}{\mu }{\rm{div}}(\nabla {p_{\rm{m}}}) = - {\phi _{\rm{m}}}{C_{\rm{m}}}\frac{{\partial {p_{\rm{m}}}}}{{\partial t}} + \rho {\alpha _{\rm{m}}}\frac{{\partial {\varepsilon _{\rm{m}}}}}{{\partial t}} + \left[ {\frac{{{\alpha _1}{k_{{\rm{nf}}}}}}{\mu }({p_{{\rm{nf}}}} - {p_{\rm{m}}}) + \frac{{{\alpha _3}{k_{{\rm{hf}}}}}}{\mu }({p_{{\rm{hf}}}} - {p_{\rm{m}}})} \right]\rho + \frac{{{k_{\rm{m}}}}}{{\mu r}}\left( {{p_{\rm{w}}} - {p_{\rm{m}}}} \right)\rho $$ (1) 天然裂缝渗流方程
$$ \rho \frac{{{k_{{\rm{nf}}}}}}{\mu }{\rm{div}}(\nabla {p_{{\rm{nf}}}}) = - {\phi _{{\rm{nf}}}}{C_{{\rm{nf}}}}\frac{{\partial {p_{{\rm{nf}}}}}}{{\partial t}} + \left[ {\frac{{{\alpha _2}{k_{{\rm{hf}}}}}}{\mu }({p_{{\rm{hf}}}} - {p_{{\rm{nf}}}}) - \frac{{{\alpha _1}{k_{{\rm{nf}}}}}}{\mu }({p_{{\rm{nf}}}} - {p_{\rm{m}}})} \right]\rho $$ (2) 人工裂缝渗流方程
$$ \rho \frac{{{k_{{\rm{hf}}}}}}{\mu }{\rm{div}}(\nabla {p_{{\rm{hf}}}}) = - {\phi _{{\rm{hf}}}}{C_{{\rm{hf}}}}\frac{{\partial {p_{nf}}}}{{\partial t}} - \left[ {\frac{{{\alpha _2}{k_{hf}}}}{\mu }({p_{{\rm{hf}}}} - {p_{{\rm{nf}}}}) + \frac{{{\alpha _3}{k_{{\rm{hf}}}}}}{\mu }({p_{{\rm{hf}}}} - {p_{\rm{m}}})} \right]\rho + \frac{{{k_{{\rm{hf}}}}}}{{\mu r}}({p_{\rm{w}}} - {p_{{\rm{hf}}}})\rho $$ (3) 人工裂缝向基质的窜流因子α1 [11 ]
$$ {\alpha _1} = 8\left( {1 + \frac{2}{{s_{\rm{n}}^2}}} \right) $$ (4) 天然裂缝向基质的窜流因子α2 [12]
$$ {\alpha _2} = \frac{{12{k_{{\rm{hf}}}}}}{{s_{\rm{n}}^2{k_{{\rm{nf}}}}}} $$ (5) 人工裂缝向基质的窜流因子α3
$$ {\alpha _3} = 8\left( {1 + \frac{2}{{s_{\rm{h}}^2}}} \right) $$ (6) 式中,km为基质渗透率,D;knf为天然裂缝渗透率,D;khf为人工裂缝渗透率,D;pm为基质内压力,MPa;pnf为天然裂缝渗流压力,MPa;phf为人工裂缝渗流压力,MPa;αm为基质有效应力系数,无量纲;εm为基质体积应变,无量纲;α1为天然裂缝向基质的窜流因子,m-2;α2为人工裂缝向天然裂缝的窜流因子,m-2;α3为人工裂缝向基质的窜流因子,m-2;ϕm、ϕnf、ϕhf分别为岩石基质、天然裂缝及人工裂缝的孔隙度,%;pw为静水柱压力,MPa;r为渗流距离,此处亦可表征井眼至渗流节点的距离,m;Cm、Cnf、Chf分别为基质、天然裂缝及人工裂缝的压缩系数,MPa-1;Sn、Sh分别为天然裂缝、人工裂缝间距,m。
-
液相侵入过程中页岩变形控制方程为[13]
$$ G{u_i}_{,jj} + \frac{G}{{1 - 2v}}{\rm{ }}{u_{j,ji}} - {\alpha _{\rm{b}}}{p_{{\rm{m}},i}} - {\beta _{{\rm{bh}}}}{p_{{\rm{hf}},i}} - {\beta _{{\rm{bh}}}}{p_{{\rm{nf}},i}} + {f_i} = 0 $$ (7) 式中,ui、uj为位移分量,m;f为体积力,MPa;αb、βbh、βbn为Biot系数,无量纲;G为剪切模量,MPa;ν为泊松比,无量纲。
在正压差作用下,液相进入岩石基质以及裂缝系统,改变地层的有效应力场,进而影响基质及裂缝系统的渗透率、孔隙度物性参数,进一步加剧液相侵入程度[14]。因此,分析钻开页岩地层后液相侵入过程中的有效地应力分布规律必须要考虑岩石渗透率、孔隙度的应力敏感性[15]。
基质孔隙度与应力的变化关系满足[16]
$$ {\phi _{\rm{m}}} + \frac{{{\phi _{{\rm{m0}}}} + {\varepsilon _v}}}{{1 + {\varepsilon _v}}} $$ (8) εv反映岩石基质的变形,在液相侵入的过程中,岩石基质的变形由应力变化引起。由应力变化引起的岩石骨架体积应变为
$$ {\varepsilon _{\rm{v}}} = \frac{{\sigma '}}{{{K_{\rm{v}}}}} $$ (9) 基质渗透率动态变化关系满足[17]
$$ \frac{{{k_{\rm{m}}}}}{{{k_{{\rm{m0}}}}}} = \frac{1}{{1 + {\varepsilon _{\rm{v}}}}}{\left( {1 + \frac{{{\varepsilon _v}}}{{{\phi _{{\rm{m}}0}}}}} \right)^3} $$ (10) 裂缝中的渗流符合立方定律,其渗透率满足
$$ {k_{\rm{f}}} = \frac{{e_{\rm{h}}^2}}{{12}} $$ (11) 裂缝渗透率受应力环境的影响极大,裂缝的动态开度满足[18]
$$ {e_{\rm{h}}} = {u_{{\rm{f0}}}} + \Delta {u_{{\rm{fn}}}} $$ (12) $$ \Delta \sigma {'_{{\rm{fn}}}} = {K_{{\rm{fn}}}}\Delta {u_{{\rm{fn}}}} $$ (13) 式中,ϕm0为基质初始状态孔隙度,%;εv为岩石骨架体积应变,无量纲;σ'为基质有效应力,MPa;Kv为岩石体积模量,MPa;km0为基质初始渗透率,D;eh为裂缝开度,m;uf0为裂缝初始开度,m;Δufn为裂缝的法向位移变化,m;Δσ'fn为裂缝有效应力变化值,MPa;Kfn为裂缝法向刚度,Pa/m。
式(12)、(13)代入式(11)后可分别算出天然裂缝系统渗透率及人工裂缝系统渗透率随应力的变化。
在上述模型中,总应力σ均被假设为常数,即作用在基质或裂缝系统中的有效应力变化等于基质或裂缝系统中的压力变化[19]
$$ d\sigma = d\sigma ' - dp = 0 $$ (14) 式中,dσ'表示有效应力,MPa。
-
Navier型方程初始条件及边界条件为[16]
ui(r, 0)=u0
σij(r, 0)=σ0
ui(+∞, t)=uib
σij(+∞, t)=Fib
渗流方程的初始条件及边界条件为
pm(r, 0)=phf(r, 0)=pnf(r, 0)=p0
pm(+∞, t)=phf(+∞, t)=pnf(+∞, t)=p0
pm(rw, t)=phf(rw, t)=pnf(rw, t)=pw
表 1 渗流模型计算参数
Table 1. Calculation parameters of seepage model
参数 数值 地层压力/MPa 20 井底压力/MPa 25 岩石体积
模量/MPa2.5×104 泊松比 0.3 基质Biot系数 0.1 天然裂缝
Biot系数0.3 人工裂缝
Biot系数0.6 基质压缩系数/
MPa-13×10–4 天然裂缝
压缩系数/MPa-12.5×10–4 人工裂缝压缩
系数/MPa-14×10–4 基质初始
渗透率/mD0.001 液相密度/(kg·m-3) 1000 液相黏度/(Pa·s) 0.001 天然裂缝间距/m 20 人工裂缝间距/m 100 天然裂缝
初始开度/m2×10–5 人工裂缝
初始开度/m2×10–3 天然裂缝法向
刚度/(Pa·m-1)1.1×109 人工裂缝法向
刚度/(Pa·m-1)1.1×109 基质初始
孔隙度/%2 天然裂缝
初始孔隙度/%0.4 人工裂缝
初始孔隙度/%0.4 运用有限元软件COMSOL对建立的流固耦合作用下页岩多重孔隙介质液相侵入数学模型进行数值求解。其求解步骤如图 1所示。计算结果见图 2~图 4。
在液相侵入初期,液相沿着人工裂缝快速推进,导致人工裂缝内的液相压力迅速上升;天然裂缝内的液相压力在初期上升速度小于人工裂缝,这是由于人工裂缝的初始开度大于天然裂缝,液相在人工裂缝内的渗流起主导作用。当液相前沿推进到天然裂缝处并在压力作用下往裂缝深部推进后,天然裂缝内的液相压力也将在短时间内迅速上升,但始终低于人工裂缝内的液相压力。由于基质渗透率远低于裂缝系统的渗透率,因此液相侵入基质的速度非常缓慢,裂缝系统对基质内的液相渗流具有非常显著的影响。随着液相的进一步侵入,裂缝系统和基质系统将会存在压力差,在压力差的作用下液相通过裂缝面向基质渗透[16],进而带动基质渗流压力升高。
-
选取文献数据对本文建立的理论模型进行验证。计算参数[21]见表 2。
表 2 二维瞬态渗流模型计算参数
Table 2. Calculation parameters of two-dimensional transient seepage model
参数 数值 初始孔隙度/% 30 井筒半径/m 0.1 裂缝渗透率/mD 1 000 黏度/(Pa·s) 1.84×10−5 基质渗透率/mD 1 地层厚度/m 1 地层压力/Pa 1×10−5 计算时间/d 10 图 5为应用本文所建立的理论模型计算所得的基质渗流压力与文献中数据对比,二者吻合程度较高。但文献没有将人工裂缝与天然裂缝视为2个不同的压力系统。在考虑人工裂缝及天然裂缝均向基质窜流的情况下,本文理论模型的数值解与文献中模型计算所得解析解有所差异。
根据图 2所示的计算结果可见基质渗流压力随径向距离的演化在近井壁地带出现拐点,分析认为拐点出现的原因是由于基质渗透率远低于裂缝,导致在初始时刻液相沿裂缝渗流,当裂缝内渗流压力远大于基质后,液相开始从裂缝及井筒内窜流至基质,此时在裂缝渗流压力以及井底液柱压力的共同作用下,基质渗流压力出现突然升高的现象。
-
为进一步探寻液相压力在基质及裂缝系统内的传递,扩大计算时长至1 h,并比较三重基质内100 s、10 min及1 h时的压力分布,如图 6所示。基质中的液相压力在初期上升缓慢,在液相侵入一段时间后,液相锋面逐渐前移并占据整个裂缝系统,整个储层的压力趋于稳定,此时,基质受裂缝的影响不如初期明显,液相压力将迅速上升,基质内的渗流压力在第10 min时已逐渐逼近同时间内的裂缝系统。由于裂缝渗透率远高于基质渗透率,因此液相在裂缝系统内的传播非常迅速,液相侵入10 min后,人工裂缝及天然裂缝内的渗流压力即趋于稳定。在液相侵入1 h后,基质系统和裂缝系统的压力差减小[16],页岩裂缝-基质内的渗流压力趋于一致,基质及裂缝系统间的液相交换速度变缓,最终保持一个稳定的量。在压力稳定前,压力波在基质内的传播速度始终低于其在裂缝系统内的传播速度。
图 6 页岩裂缝-基质渗流压力不同时刻对比
Figure 6. Comparison between the seepage pressure in shale fracture and that in matrix at different moments
分析页岩裂缝-基质内两种工况下的渗流压力,如图 7所示。
图 7 2种工况下页岩裂缝-基质渗流压力对比
Figure 7. Comparison between the seepage pressure in shale fracture and that in matrix under two working conditions
流固耦合情况下的渗流压力始终大于同时段内非耦合工况下的渗流压力。这是由于流固耦合情况下裂缝和基质的渗透率均呈现增加趋势[16]。液相侵入过程中,裂缝中的渗流压力大于孔隙压力,高渗流压力将诱发裂缝变形,裂缝系统开度增加,且由于裂缝扩张、新微细裂纹的产生以及缝间交会,增加了岩心裂缝的复杂程度、裂缝体积及其连通性[22],渗流能力增强。基质受到孔隙压力和围压的双重作用,液相侵入导致超孔隙压力,进而引起基质渗透率的增加[16]。这将导致流固耦合作用下的液相前缘在相同时间内推进到更远的位置,并向储层深部渗透。
-
(1)本文建立了页岩地层复杂孔隙介质液相侵入数学模型。该模型在将人工裂缝及天然裂缝视为两个不同的压力系统的基础上考虑了液相侵入过程中的流固耦合作用,并考虑了人工裂缝与天然裂缝以及裂缝系统及基质间的窜流,更符合页岩地层真实渗流情况。
(2)人工裂缝的初始渗透率较大,将导致基质及天然裂缝内初期渗流压力增加缓慢,随着压力差的减小,复杂孔隙介质的渗流压力将趋于一致,基质及裂缝系统间的液相窜流速率放缓,并最终趋于稳定;流固耦合作用对渗流影响显著,裂缝系统及基质的渗透率在流固耦合作用下增大,同时间段内耦合工况下的渗流压力大于非耦合工况下的渗流压力。
Study on the invasion mechanism of liquid phase into shale formation under the effect of fluid-solid coupling
-
摘要: 页岩地层液相侵入将会导致井壁失稳等井下事故的发生,但目前关于页岩地层液相侵入及流固耦合现象研究尚不深入。建立了流固耦合作用下的页岩地层液相侵入数学模型,并对模型求解结果进行了验证和分析。结果表明:侵入初期,液相在人工裂缝内的渗流占据主导作用,但随着侵入时间的延长,复杂孔隙介质内的压力差将逐渐减弱,基质及裂缝系统内的渗流压力将趋于一致;流固耦合作用在液相侵入过程中影响显著,裂缝系统及基质的渗透率在流固耦合作用下增大,同时间段内耦合工况下的渗流压力大于非耦合工况下的渗流压力。将基质、天然裂缝、人工裂缝视为3个不同压力系统基础上得到的研究结果对科学认识流固耦合作用下的页岩地层液相侵入机理具有一定的指导意义。Abstract: The invasion of liquid phase into shale formation can lead to downhole accidents, e.g. borehole instability. So far, however, the invasion of liquid phase into shale formation and the fluid-solid coupling phenomenon have not been studied deeply. The mathematical model for the invasion of liquid phase into shale formation under the effect of fluid-solid coupling was established. In addition, its solution results were verified and analyzed. It is shown that in the early stage of liquid invasion, its seepage in artificial fractures is dominant. With the extension of invasion time, however, the pressure difference in complex porous medium decreases gradually, and the seepage pressure in the matrix tends to be consistent with that in the fracture system. The effect of fluid-solid coupling in the process of liquid invasion is significant, and the permeability of fracture system and matrix increases under the effect of fluid-solid coupling. And in the same time interval, the seepage pressure under the working condition of coupling is higher than that of non-coupling. In conclusion, the research results which are obtained by taking matrix, natural fracture and artificial fracture as three different pressure systems play a guiding role for understanding the invasion mechanisms of liquid phase into shale formation under the effect of fluid-solid coupling scientifically.
-
Key words:
- shale /
- complex porous medium /
- fluid-solid coupling /
- liquid invasion /
- mechanism
-
表 1 渗流模型计算参数
Table 1. Calculation parameters of seepage model
参数 数值 地层压力/MPa 20 井底压力/MPa 25 岩石体积
模量/MPa2.5×104 泊松比 0.3 基质Biot系数 0.1 天然裂缝
Biot系数0.3 人工裂缝
Biot系数0.6 基质压缩系数/
MPa-13×10–4 天然裂缝
压缩系数/MPa-12.5×10–4 人工裂缝压缩
系数/MPa-14×10–4 基质初始
渗透率/mD0.001 液相密度/(kg·m-3) 1000 液相黏度/(Pa·s) 0.001 天然裂缝间距/m 20 人工裂缝间距/m 100 天然裂缝
初始开度/m2×10–5 人工裂缝
初始开度/m2×10–3 天然裂缝法向
刚度/(Pa·m-1)1.1×109 人工裂缝法向
刚度/(Pa·m-1)1.1×109 基质初始
孔隙度/%2 天然裂缝
初始孔隙度/%0.4 人工裂缝
初始孔隙度/%0.4 表 2 二维瞬态渗流模型计算参数
Table 2. Calculation parameters of two-dimensional transient seepage model
参数 数值 初始孔隙度/% 30 井筒半径/m 0.1 裂缝渗透率/mD 1 000 黏度/(Pa·s) 1.84×10−5 基质渗透率/mD 1 地层厚度/m 1 地层压力/Pa 1×10−5 计算时间/d 10 -
[1] 刘厚彬, 孟英峰, 李皋, 李平, 邓元洲.泥页岩水化作用对岩石强度的影响[J].钻采工艺, 2010, 33(6):18-20. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zcgy201006006 LIU Houbin, MENG Yingfeng, LI Gao, LI Ping, DENG Yuanzhou. Theoretical simulation and experimental evaluation of the effect of hydration on the shale rock strength[J]. Drilling & Production Technology, 2010, 33(6): 18-20. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zcgy201006006 [2] 梁利喜, 丁乙, 刘向君, 许丽.硬脆性泥页岩井壁稳定渗流-力化耦合研究[J].特种油气藏, 2016, 23(2):140-143. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=tzyqc201602034 LIANG Lixi, DING Yi, LIU Xiangjun, XU Li. Seepage- mechanochemistry coupling of wellbore stability in hard- brittle shale[J]. Special Oil & Gas Reservoirs, 2016, 23(2): 140-143. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=tzyqc201602034 [3] 曹园, 邓金根, 蔚宝华, 赵凯.深部泥页岩水化特性研究[J].科学技术与工程, 2014, 14(6):118-120. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=kxjsygc201406026 CAO Yuan, DENG Jingen, WEI Baohua, ZHAO Kai. Hydration properties research of deep formation shale[J]. Science Technology and Engineering, 2014, 14(6): 118-120. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=kxjsygc201406026 [4] 胡进科, 李皋, 孟英峰.页岩气钻井过程中的储层保护[J].天然气工业, 2012, 32(12):66-70. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=trqgy201212014 HU Jinke, LI Gao, MENG Yingfeng. Reservoir protection in the process of shale gas drilling[J]. Nature Gas Industry, 2012, 32(12): 66-70. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=trqgy201212014 [5] 康毅力, 皇凡生, 游利军, 俞杨烽, 练章华.钻井液浸泡页岩裂缝宽度的模拟及应用[J].石油钻采工艺, 2014, 36(5):41-46. http://www.syzcgy.cn/oa/DArticle.aspx?type=view&id=201405011 KANG Yili, HUANG Fansheng, YOU Lijun, YU Yangfeng, LIAN Zhanghua. Simulation and application of shale fracture width immersed in drilling fluid[J]. Oil Drilling & Production Technology, 2014, 36(5): 41-46. http://www.syzcgy.cn/oa/DArticle.aspx?type=view&id=201405011 [6] BLASKOVICH F T, CAIN G M, SONIER F, WALDREN D. A multicomponent isothermal system for efficient reservoir simulation[C]. Middle East Oil Technical Conference and Exhibition, Bahrain, SPE 11480, 1983. [7] NOORISHAD J, MEHRAN M. An upstream finite element method for solution of transient transport equation in fractured porous media[J]. Water Resources Research, 1982, 18(3): 588-596. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1029/WR018i003p00588 [8] MEYER B R, BAZAN L W. A discrete fracture network model for hydraulically induced fractures - theory, parametric and case studies[C]. SPE Hydraulic Fracturing Technology Conference, Texas, 2011. [9] 卢义玉, 刘小川, 汤积仁, 黄飞.热流固耦合作用下页岩渗透特性实验[J].重庆大学学报, 2016, 39(1):65-71. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=cqdxxb201601009 LU Yiyu, LIU Xiaochuan, TANG Jiren, HUANG Fei. Effects of heat flow fluid-solid coupling on the characteristics of shale permeability[J]. Journal of Chongqing University, 2016, 39(1): 65-71. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=cqdxxb201601009 [10] 康永尚, 邓泽, 王红岩, 刘洪林, 袁春林, 赵群.流-固耦合物理模拟实验及其对页岩压裂改造的启示[J].地球科学, 2016, 41(8):1376-1383. http://qikan.cqvip.com/Qikan/Article/Detail?id=669890903 KANG Yongshang, DENG Ze, WANG Hongyan, LIU Honglin, YUAN Chunlin, ZHAO Qun. Fluid-solid coupling physical experiments and their implications for fracturing stimulations of shale gas reservoirs[J]. Earth Science, 2016, 41(8): 1376-1383. http://qikan.cqvip.com/Qikan/Article/Detail?id=669890903 [11] WU Yu, LIU Jishan, ELSWORTH D, CHEN Zhongwei. Dual poroelastic response of a coal seam to CO2 injection[J]. International Journal of Greenhouse Gas Control, 2010, 4(4): 668-678. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=59312a0b4aae9f482e88208fc4f9039f [12] 顾岱鸿, 丁道权, 刘军, 丁志文, 刘锦华, 朱智.三重介质页岩气藏分段压裂水平井产能预测模型[J].大庆石油地质与开发, 2016, 35(1):158-165. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqsydzykf201601032 GU Daihong, DING Daoquan, LIU Jun, DING Zhiwen, LIU Jinhua, ZHU Zhi. Productivity predicting model of the staged fractured horizontal well in triple-media shale gas reservoirs[J]. Petroleum Geology and Oilfield Development in Daqing, 2016, 35(1): 158-165. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqsydzykf201601032 [13] 赵强.裂缝流固耦合渗流机理研究[D].成都: 西南石油大学, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10615-1015599183.htm ZHAO Qiang. Study on fluid-solid coupling flow mechanism of fracture[D]. Chengdu: Southwest Petroleum University, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10615-1015599183.htm [14] 曾辉.致密砂岩气藏压裂缝耦合损害机理研究[D].成都: 西南石油大学, 2017. http://cdmd.cnki.com.cn/Article/CDMD-10615-1017108614.htm ZENG Hui. Coupled damage mechanisms of hydraulic fracture in tight gas reservoirs[D]. Chengdu: Southwest Petroleum University, 2017. http://cdmd.cnki.com.cn/Article/CDMD-10615-1017108614.htm [15] 罗兵.致密气层气体钻开井周应力动态分布规律研究[D].成都: 西南石油大学, 2014. http://cdmd.cnki.com.cn/Article/CDMD-10615-1014409728.htm LUO Bing. Study on the dynamic distribution law of persimal stress in tight gas layer drilling[D]. Chengdu: Southwest Petroleum University, 2014. http://cdmd.cnki.com.cn/Article/CDMD-10615-1014409728.htm [16] 任文希.孔隙-裂缝性致密砂岩液相侵入过程流固耦合作用研究[D].成都: 西南石油大学, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10615-1015599431.htm REN Wenxi. Fluid-solid interaction in liquid-phase intrusion of pore-fractured tight sandstone[D]. Chengdu: Southwest Petroleum University, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10615-1015599431.htm [17] 冉启全, 李士伦.流固耦合油藏数值模拟中物性参数动态模型研究[J].石油勘探与开发, 1997, 24(3):61-65. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199700493745 RAN Qiquan, LI Shilun. Study on dynamic models of reservoir parameters in the coupled simulation of multiphase flow and reservoir deformation[J]. Petroleum Exploration and Development, 1997, 24(3): 61-65. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199700493745 [18] ZHAO Zhihong, JING Lanru, NERETNIEKS I, MORENO L. Analytical solution of coupled stress-flow-transport processes in a single rock fracture[J]. Computers & Geosciences, 2011, 37(9): 1437-1449. http://www.sciencedirect.com/science/article/pii/S0098300411000902 [19] 吴宇.煤层中封存二氧化碳的双重孔隙力学效应研究[D].北京: 中国矿业大学, 2010. WU Yu. Dual poroelastic response of coal to CO2 sequestration[D]. Beijing: China University of Mining and Technology, 2010. [20] 程万.三维空间下裂缝性页岩储层水力裂缝扩展机理研究[D].北京: 中国石油大学(北京), 2016. CHENG Wan. Mechanism of hydraulic fracture propagation in fractured shale reservoir in three dimensional space[D]. Beijing: China University of Petroleum, 2016. [21] 范鑫.页岩基质流固耦合数值模拟与人工裂缝参数优化研究[D].北京: 中国石油大学(北京), 2016. FAN Xin. Research on hydro-mechanical coupled simulation in shale matrix and parameter optimization of hydraulic fracture[D]. Beijing: China University of Petroleum, 2016. [22] 宋付权, 胡箫, 纪凯, 黄小荷.考虑流固耦合影响的页岩力学性质和渗流规律[J].天然气工业, 2017, 37(7):69-75. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=trqgy201707010 SONG Fuquan, HU Xiao, JI Kai, HUANG Xiaohe. Effect of fluid–solid coupling on shale mechanics and seepage laws[J]. Natural Gas Industry, 2017, 37(7): 69-75. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=trqgy201707010 -