基于相场法的深层干热岩储层水力压裂模拟研究

戴一凡; 侯冰; 廖志豪

doi:10.11911/syztjs.2024047

基于相场法的深层干热岩储层水力压裂模拟研究

戴一凡^{1, 2,},
侯冰^{1, 2, 3, ,},
廖志豪^{1, 2, 3}

1.
油气资源与工程全国重点实验室（中国石油大学（北京））, 北京 102249
2.
石油工程教育部重点实验室（中国石油大学（北京））, 北京 102249
3.
中国石油大学（北京）克拉玛依校区石油学院, 新疆克拉玛依 834000

基金项目: 国家重点研发计划项目“深部含煤岩系超临界CO₂穿层压裂–驱替–封存评价技术研究”（编号：2022YFE0129800）、国家自然科学基金项目“超深碳酸盐岩断缝体储层酸压裂缝与天然裂缝–溶洞相互作用机理”（编号：52074311）、国家自然科学基金重点项目“提高超深大科度井压裂效率的关键力学问题研究”（编号：52334001）、中国石油大学（北京）青年拔尖人才项目“超深层近井裂缝扩展转向流固耦合力学机理研究”（编号：ZX20230042）联合资助。

详细信息

作者简介:
戴一凡（1996—），男，江苏泰州人，2018年毕业于中国石油大学（华东）石油工程专业，中国石油大学（北京）石油与天然气工程专业在读博士研究生，主要从事岩石力学方面的研究。E-mail: 819673313@qq.com

通讯作者:
侯冰，binghou@vip.163.com

中图分类号: P314; TE357.1
计量
- 文章访问数: 190
- HTML全文浏览量: 61
- PDF下载量: 89
出版历程
- 收稿日期: 2023-11-30
- 修回日期: 2024-02-29
- 录用日期: 2024-04-20
- 网络出版日期: 2024-04-24
- 刊出日期: 2024-04-02

Simulation of Hydraulic Fracturing in Deep Hot Dry Rock Reservoir Based on Phase-Field Method

DAI Yifan^{1, 2,},
HOU Bing^{1, 2, 3, ,},
LIAO Zhihao^{1, 2, 3}

1.
State Key Laboratory of Petroleum Resources and Engineering(China University of Petroleum(Beijing)), Beijing, 102249, China
2.
MOE Key Laboratory of Petroleum Engineering(China University of Petroleum(Beijing)), Beijing, 102249, China
3.
Department of Petroleum Engineering, China University of Petroleum-Beijing at Karamay, Karamay, Xinjiang, 834000, China

摘要

摘要:
为认识深层干热岩储层水力裂缝起裂与扩展规律，基于相场法建立了流−固−热耦合裂缝扩展数值模型，分析了干热岩储层水力压裂裂缝扩展力学行为，研究了温度、排量及天然裂缝等因素对水力裂缝扩展的影响。结果表明，基于相场法建立的流−固−热耦合裂缝扩展数值模型判别准则简单，计算准确；干热岩储层水力压裂过程中，热应力呈现为拉应力，有助于增加裂缝开度，对裂缝扩展起辅助作用；压裂液与地层温度差越大，排量越大，热应力影响效果越显著。水力压裂裂缝遇到天然裂缝后会沟通并开启天然裂缝，再次起裂时受地应力与天然裂缝控制。研究结果可为深层干热岩储层水力压裂提供参考。
- 干热岩 /
- 相场法 /
- 裂缝扩展 /
- 多物理场耦合 /
- 数值模拟 /
- 自适应网格 /
- 天然裂缝
Abstract:
In order to study the hydraulic fracture initiation and propagation mechanisms in deep hot dry rock (HDR) reservoirs, a fluid-solid-heat coupling numerical model for fracture propagation simulation was established based on the phase-field method. This model was used to analyze the mechanical behavior of hydraulic fracture propagation in HDR reservoirs, as well as the effects of factors such as temperature, pump rate, and natural fractures on the hydraulic fracture propagation. The results show that this fluid-solid-heat coupling numerical model has a relatively simple criterion with high calculation accuracy. During the process of hydraulic fracturing in HDR reservoirs, thermal stress acts as tensile stress, which contributes to increasing fracture aperture and facilitating fracture propagation. The larger temperature difference between fracturing fluid and formation, and the higher pumping rates, the more significant effect of the thermal stress. After encountering natural fractures, hydraulic fractures will communicate and initiate those natural fractures, and the re-initiation of fractures is controlled by in-situ stress and natural fractures. The findings of this study can provide a good reference for guiding hydraulic fracturing in deep HDR reservoirs.
- hot dry rock /
- phase-field method /
- fracture propagation /
- multi-physical field coupling /
- numerical simulation /
- adaptive mesh /
- natural fractures

HTML全文

深层干热岩储层具有温度高、硬度大和基质渗透率低等特点，需采用水力压裂等方法改造致密干热岩储层，构建复杂的裂缝网络体系，提高循环工质换热体积，实现高效循环取热^[1–6]。干热岩水力压裂是一个多物理场耦合问题，需要考虑位移、流体压力、温度等多个物理量之间的相互影响^[7–17]。H. Hofmann等人^[18]通过研究发现，复杂裂缝网络有助于提高干热岩储层开采时的经济效益。Zhang Yanjun等人^[19]提出，水平井分段压裂技术可提高循环工质的换热效率，从而提高裂缝型干热岩储层的开采效益。谢紫霄等人^[20]通过干热岩压裂物理模拟试验，发现水力裂缝沟通天然裂缝后能显著提高裂缝渗流系数，实现地热能的高效开发。干热岩储层裂缝扩展过程中受多因素综合影响，且温度影响不容忽略。

相场法是一种基于系统能量最小化原理判断裂缝扩展的有限元方法，在处理多物理场耦合问题时具有网格处理简单、结果准确性高和计算量小等优点，近年来被广泛应用于水力压裂数值模拟领域^[21–26]。侯冰等人^[23]基于相场法建立了流−固−化耦合酸压裂缝扩展模型，分析了孔洞结构对酸压裂缝的影响。N. Noii等人^[25]提出适用于正交各向异性的相场裂缝描述方法，对相场裂缝建模进行了初步分析。

相场法在处理复杂地层多物理场耦合裂缝扩展模拟方面具有显著优势。基于相场法建立了流−固−热耦合裂缝扩展数值模型，采用自适应网格方法平衡计算成本与计算精度，通过解析解验证模拟结果的准确性；在此基础上，研究了干热岩储层水力裂缝的起裂与扩展规律，分析了地质与工程因素对干热岩储层水力裂缝扩展的影响，以期为深层干热岩储层水力压裂提供参考。

1. 深层干热岩储层水力压裂数值模型

基于相场法构建流−固−热耦合裂缝扩展数值模型时，其核心思想是引入无量纲相场参数ϕ描述空间点的物理状态。ϕ的取值在 [0,1]区间，ϕ= 0 表示材料断裂，ϕ= 1 表示材料未断裂， 0 < ϕ < 1 表示材料处于裂缝过渡区域。该参数可以描述裂缝形态，同时可作为指示参数耦合多物理场方程。

1.1 控制方程

1.1.1 相场−位移控制方程

在热弹性多孔介质的相场裂缝模型能量泛函基础上^[24–25]，将裂缝扩展问题转化为系统能量最小化问题并求解，可得相场−位移控制方程强形式：

$\begin{split} - [(1 - k)\phi {\bf{\varepsilon }}({{u}}) - {G_{\mathrm{c}}}\varepsilon \nabla \phi - \frac{{{G_{\mathrm{c}}}}}{\varepsilon }1 - \phi - 2( {{\alpha _{\rm{B}}} } & - 1 )\phi \left( {p - {p_0}} \right) \nabla \cdot {{u}}+ 2\phi \nabla \left( {p - {p_0}} \right){{u}} - 2\left( {{\alpha _{\mathrm{T}}}{K_{\mathrm{d}}} + C} \right)\phi \left( {\theta - {\theta_0}} \right)\nabla \cdot u + \\ & 2C\phi \nabla \left( {\theta - {\theta_0}} \right){{u}}]{\partial _t}\phi = 0\end{split}$

(1)

$\begin{split} - \nabla \cdot [((1 - k){\phi ^2} + k){\sigma ^ + }({{u}}) + {\sigma ^ - }({{u}})] + (\alpha - 1)\nabla ({\phi ^2}p) + {\phi ^2}\nabla p + ({\alpha _{\mathrm{T}}}{K_{\mathrm{d}}} + {C_{\mathrm{T}}})\nabla ({\phi ^2}\theta) + {\phi ^2}\nabla \theta = 0 \end{split}$

(2)

式中：u为位移，m；k为正则化参数；ε(u)为应变；ε为正则化参数；α_B为Biot系数；α_T为热膨胀系数，℃^－1 ；p为流体压力，Pa；p₀为初始压力，Pa； $\theta$ 为温度，℃； $\theta_0$ 为初始温度，℃；C为界面热应力常数，Pa/℃；σ⁺(u)为裂缝驱动力，Pa；σ^－(u)为裂缝非驱动力，Pa；在二维问题中，K_d=μ+λ，μ，λ为Lame系数；G_c为临界能量释放率，Pa·m。

1.1.2 温度场控制方程

基于相场参数的定义，可构建指示函数：

${X}\mathrm{_f}=1-\phi\text{，}{X}_{\mathrm{r}}=\phi$

(3)

基于局部热平衡假设，同一空间点处温度值具有唯一性。引入指示函数后，得到温度场控制方程的有限元强形式：

$\begin{split} {{(}}{{{X}}_{\mathrm{f}}}{C_{\mathrm{f}}} + {{{X}}_{\text{r}}} & {C_{\mathrm{r}}}) \frac{{\partial \theta}}{{\partial t}} = ({{{X}}_{\mathrm{f}}}{K_{\mathrm{f}}} + {{{X}}_{\text{r}}}{K_{\mathrm{r}}})\left[ \frac{\partial }{{\partial x}}\left( {\frac{{\partial \theta}}{{\partial x}}} \right)\right. +\\ & \left. \frac{\partial }{{\partial y}}\left( {\frac{{\partial \theta}}{{\partial y}}} \right) \right] +v\nabla \theta + {Q_{\mathrm{i}}}\end{split}$

(4)

式中：C_f为压裂液热容，kJ·m³/℃ ；C_r为储层热容，kJ·m³/℃ ；K_f为压裂液热传导系数，W/(m·℃)；K_r为地层岩石热传导系数，W/(m·℃)； v为流体流速，m/s；Q_i为注入热源，J。

1.1.3 流体压力控制方程

对于流体场，其物质守恒方程为：

$\partial_t\left(\rho\mathrm{_f}\varphi^*\right)+\nabla\cdot\left(\rho_{\mathrm{f}}v\right)=q$

(5)

式中：ρ_f为液体密度，g/cm³；φ^*为孔隙度；q为单位时间内注入的压裂液质量，g。

计算流体流速时忽略重力影响，表达式为：

$v = - \frac{K}{\eta }\Delta p$

(6)

式中：K为渗透率，D；η为液体黏度，mPa·s。

$\varphi^*=\varphi_0^*+\alpha\mathrm{_B}\nabla\cdot u+\frac{1}{M}(p-p_0)+K\mathrm{_T}(\theta -\theta_0)$

(7)

式中： $\varphi_0^*$ 为初始时刻孔隙度；M为Biot系数；K_T为孔隙热膨胀系数，℃⁻¹。

流体压力控制方程的强形式为：

$\rho\mathrm{_f}\partial_t\left(\alpha_{\mathrm{B}}\nabla\cdot u+\frac{1}{M}p+K_{\mathrm{T}}\theta \right)-\rho_{\mathrm{f}}\frac{K}{\eta}\nabla\cdot\Delta p=q$

(8)

1.2 自适应网格

模型中以相场参数为指示参数，引入自适应网格技术，在保证模拟精度的前提下，减少模型的计算量。自适应网格算法流程如图1所示。

图 1 自适应网格处理流程

Figure 1. Adaptive mesh processing procedure

下载: 全尺寸图片幻灯片

自适应网格可对同一时间步内的裂缝扩展状态进行多次求解，但与全局细分网格相比节点数量少、计算简单，计算精度仍保持良好^[26]。

1.3 多物理场耦合求解

模型共包括固体场模块、温度场模块和流体场模块等3个模块，模型中包括4个解变量，分别为相场参数（ϕ）、位移（u）、温度（T）和流体压力（p）。在固体场模块中，采用全耦合方法求解相场参数与位移；利用温度场模块、流体场模块分别求解温度和压力值，并采用固定应力迭代的方式与固体场耦合。具体流程如下：

1）在任意时间步tⁿ⁺¹，给定上一个时间步的解θⁿ，uⁿ，ϕⁿ和pⁿ。

2）基于固定应力法，构建迭代序列l=0，1，2…，定义初始解变量u_l₌₀=uⁿ，ϕ_l₌₀=ϕⁿ，p_l₌₀=pⁿ，θ_l₌₀=θⁿ。

3）给定θⁿ，uⁿ，ϕⁿ，pⁿ，u_l，ϕ_l和p_l，求解式（4），得到θ_l₊₁。

4）给定θⁿ，uⁿ，ϕⁿ，pⁿ，u_l，ϕ_l和θ_l，求解式（8）获得p_l₊₁。

5）给定θⁿ，uⁿ，ϕⁿ，pⁿ，p_l和θ_l，求解式（1）和式（2），可得{u_l₊₁，ϕ _l₊₁}。

6）若max{||p_l+1−p_l||，||θ_l+1−|θ_l||，||u_l+1−u_l||，||ϕ_l+1−ϕ_l||}≤TOL_fs，则迭代收敛（TOL_fs为固定应力法容差，取10^－8），令pⁿ⁺¹=p_l+1，θⁿ⁺¹=θ_l+1，uⁿ⁺¹=u_l+1，ϕⁿ⁺¹=ϕ_l+1；若迭代不收敛，则令u_l+1=u_l，ϕ_l+1=ϕ_l，p_l+1=p_l，θ_l+1=θ_l，重复步骤3）—5），直至满足收敛条件。

7）检查网格系统是否满足计算精度要求，若不满足，则更新自适应网格系统。

8）输出θⁿ⁺¹，uⁿ⁺¹， ϕⁿ⁺¹和pⁿ⁺¹作为时间步tⁿ⁺¹的解变量，并进行下一个时间步的计算。

1.4 初始条件与边界条件

模拟二维条件下的裂缝扩展过程，忽略重力影响与垂向应变。模型初始状态及边界条件如图2所示。

图 2 相场裂缝模型初始状态

Figure 2. Initial condition of phase-field fracture model

下载: 全尺寸图片幻灯片

模型长30 cm，宽30 cm。设置天然裂缝等结构时，将对应位置相场参数固定为0，其余岩石基质区域处相场参数设置为1，进行模拟分析，根据相场参数判断裂缝扩展情况。设定压裂液温度为20 ℃，岩石弹性模量为40 GPa，泊松比为0.20，岩石密度为2.50 g/cm³，Biot系数为0.80，渗透率为1 mD，地层岩石热容为2.0 MJ·m³/℃ ，压裂液热容4.2 MJ·m³/℃ ，地层岩石热传导系数为0.63 W/(m·℃)。

1.5 解析解验证

相场裂缝模型中，最小网格尺寸受网格细分次数控制。采用自适应网格方法，全局细分5次、局部细分3次时，最小网格尺寸为1.17 mm；采用全局网格细分方法，全局细分8次时，最小网格尺寸同样为1.17 mm。固定边界位移，分别采用全局细分网格与自适应网格对裂缝扩展进行模拟，模拟结果如图3所示。

图 3 固定边界位移条件下干热岩水力压裂数值模拟

Figure 3. Numerical simulation of hydraulic fracturing in HDR reservoir under fixed boundary displacement

下载: 全尺寸图片幻灯片

无限大线弹性物体内有一条缝内压力均匀分布的裂缝，其裂缝宽度的计算公式为^[27]：

$w = \frac{{4\left( {1 - {\nu ^2}} \right){p_{\text{f}}}}}{E}\sqrt {{c^2} - {y^2}}$

(9)

式中：w为裂缝宽度，mm；p_f 为裂缝内流体压力，Pa；c为裂缝半缝长，mm；y为裂缝内某点到裂缝中心点的距离，mm； $\nu$ 为岩石泊松比；E为岩石弹性模量，GPa。

将流体压力代入式（9），计算裂缝宽度解析解，并与不同网格方案下裂缝宽度对比，结果如图4所示。

图 4 沿y轴方向缝宽分布情况

Figure 4. Fracture width distribution along y-axis direction

下载: 全尺寸图片幻灯片

模拟结果表明，裂缝宽度的数值解与解析解吻合程度较好，自适应网格与全局细分网格的模拟结果相接近。自适应网格方法节点少，算力消耗低，计算相场裂缝模型时有显著优势。

2. 地质与工程因素对裂缝扩展的影响

2.1 温度

温度场的变化会造成岩石体积变化，在系统内部互相约束，从而产生热应力^[28–29]。为研究温度对水力裂缝起裂与扩展的影响，设置初始地层温度分别为20，100和200 ℃，压裂液注入温度均为20 ℃，排量均为20 mL/min，边界应力分别为10和7 MPa，裂缝长度变化趋势如图5所示。

图 5 不同地层温度下干热岩储层水力压裂缝长变化趋势

Figure 5. Variation trend of hydraulic fracture length in HDR reservoir with different formation temperature

下载: 全尺寸图片幻灯片

模拟过程中，由于压裂液温度低于地层温度，热应力呈现为拉应力，且地层温度越高，热应力越大；由于水力裂缝以拉伸缝为主，热应力在裂缝扩展过程中起辅助作用。压裂液排量相同时，地层温度越高，裂缝扩展速度越快，起裂压力越低。真实干热岩储层中，岩石不同矿物热膨胀时存在各向异性，热应力方向与大小更为复杂，更容易产生局部应力集中；若地层中微裂缝发育，热应力易激活并诱导微裂缝起裂，提高水力裂缝复杂程度，且热应力越大，微裂缝激活程度越高。

2.2 压裂液排量

为探究不同压裂液排量对水力裂缝扩展的影响，设置压裂液排量分别为10，20和30 mL/min，地层初始温度均为200 ℃，压裂液温度均为20 ℃，边界应力分别为10和7 MPa，注液30 s后，模拟结果如图6所示。

图 6 不同压裂液排量下干热岩储层水力压裂模拟结果

Figure 6. Simulation results of hydraulic fracturing in HDR reservoir under different pumping rates of fracturing fluids

下载: 全尺寸图片幻灯片

从图6可以看出，压裂液排量为10 mL/min时，缝长为11.63 cm；压裂液排量为20 mL/min时，缝长为20.26 cm；压裂液排量为30 mL/min时，缝长为27.48 cm。研究表明，随着压裂液排量增加，液体憋压速度更快，裂缝周边岩石温度降低更明显。干热岩储层水力裂缝扩展受流体压力和热应力的共同驱动，裂缝扩展速度更快，裂缝宽度也随之增大，更有利于后续循环介质的注入与热交换。

2.3 天然裂缝夹角

为探究天然裂缝对水力裂缝扩展的影响，设置天然裂缝与水平最大主应力方向的夹角分别为90°，45°和20°，地层初始温度均为200 ℃，压裂液温度均为20 ℃，压裂液排量为20 mL/min，边界应力分别为10和7 MPa，模拟结果如图7所示。

图 7 不同天然裂缝夹角下干热岩储层水力压裂模拟结果

Figure 7. Simulation results of hydraulic fracturing in HDR reservoir with different angles of natural fractures

下载: 全尺寸图片幻灯片

水力压裂裂缝起裂后沿水平最大主应力方向扩展，遇到天然裂缝后发生转向并沟通天然裂缝。天然裂缝与水平最大主应力方向夹角为90°时，水力裂缝转向不明显，沟通天然裂缝后在天然裂缝两侧同时起裂并扩展，扩展角度与最大水平地应力夹角约为45°，2条裂缝呈分离趋势。岩石温度在水力压裂裂缝与天然裂缝沟通处降幅较大，主要原因为压裂液在该区域流动情况更复杂，传热效率高。

天然裂缝与水平最大主应力方向夹角为45°时，水力裂缝转向效果较明显，沟通天然裂缝后在天然裂缝两侧起裂并扩展，左右两侧扩展速度差异明显，左侧扩展速度更快，这是因为该侧与水平最大主应力方向夹角较小，扩展速度更快。另外，岩石温度在水力裂缝与天然裂缝右侧夹角处降低幅度较大，主要是该侧天然裂缝与水力压裂裂缝距离较近，传热效率更高。

天然裂缝与水平最大主应力方向夹角为20°时，天然裂缝诱导水力压裂裂缝产生转向现象。水力压裂裂缝沟通天然裂缝后在天然裂缝两侧起裂并扩展，左右两侧扩展速度差异明显，左侧扩展速度更快，裂缝沿水平最大主应力方向继续扩展。岩石温度在水力压裂裂缝与天然裂缝右侧夹角处降低幅度较大。

基于上述模拟结果可知，水力压裂裂缝遇到天然裂缝后易沟通天然裂缝，天然裂缝与水平最大主应力夹角越小，转向现象越明显。水力压裂裂缝与天然裂缝沟通位置流动情况复杂，传热效率高。水力压裂裂缝再次起裂方向受到地应力与天然裂缝方向的共同作用，在水平最大主应力方向夹角较小的方向扩展迅速，在水平最大主应力方向夹角较大的方向扩展较慢，且该侧岩石降温幅度更大。

2.4 多天然裂缝干扰下裂缝扩展

为探究多天然裂缝干扰下干热岩储层水力裂缝扩展力学行为，设置不同地应力与排量进行模拟。单侧天然裂缝条数为2条，天然裂缝与水平最大主应力方向夹角为45°，地层初始温度为200 ℃，压裂液温度为20 ℃，压裂液排量分别为20 和50 mL/min，边界应力差分别为3和8 MPa，模拟结果如图8所示。

图 8 不同天然裂缝密度下干热岩储层水力压裂模拟结果

Figure 8. Simulation results of hydraulic fracturing in HDR reservoir with different densities of natural fractures

下载: 全尺寸图片幻灯片

地层中存在多条天然裂缝时，水力压裂裂缝优先与距离最近的天然裂缝发生相互作用。地应力差为3 MPa、压裂液排量为20 mL/min时，水力压裂裂缝在靠近天然裂缝时轻微转向并开启初条天然裂缝，随后在天然裂缝尖端二次起裂；地应力差为8 MPa、压裂液排量为20 mL/min时，水力压裂裂缝靠近天然裂缝无明显转向现象，开启初条天然裂缝后在裂缝尖端二次起裂；地应力差为8 MPa、压裂液排量增加至50 mL/min时，水力压裂裂缝贯穿初条天然裂缝并继续扩展，由于水力压裂裂缝贯穿天然裂缝后会消耗水力能量，因此仅能贯穿单条天然裂缝。

3. 结　论

1）基于相场法建立了流−固−热多场耦合裂缝扩展数值模型，根据能量最小化原理计算裂缝形态，无需额外添加裂缝扩展准则，模型搭建简单，可用于干热岩储层水力裂缝扩展模拟研究。

2）干热岩储层水力压裂过程中，热应力体现为拉应力，有助于增加裂缝开度，在裂缝扩展过程中起辅助作用。压裂液与地层温度差越高，排量越大，热应力影响效果越显著；若地层中微裂缝发育，热应力易激活并诱导微裂缝起裂，提高水力压裂裂缝的复杂程度。

3）地层中存在多条天然裂缝时，水力压裂裂缝优先与距离最近的天然裂缝发生相互作用。水力压裂裂缝遇到天然裂缝后易沟通天然裂缝，天然裂缝与水平最大主应力方向夹角越小，地应力差越小，转向现象越明显。水力压裂裂缝再次起裂时沿天然裂缝方向起裂，且向水平最大主应力方向偏转。

图 1 自适应网格处理流程

Figure 1. Adaptive mesh processing procedure

下载: 全尺寸图片幻灯片

图 2 相场裂缝模型初始状态

Figure 2. Initial condition of phase-field fracture model

下载: 全尺寸图片幻灯片

图 3 固定边界位移条件下干热岩水力压裂数值模拟

Figure 3. Numerical simulation of hydraulic fracturing in HDR reservoir under fixed boundary displacement

下载: 全尺寸图片幻灯片

图 4 沿y轴方向缝宽分布情况

Figure 4. Fracture width distribution along y-axis direction

下载: 全尺寸图片幻灯片

图 5 不同地层温度下干热岩储层水力压裂缝长变化趋势

Figure 5. Variation trend of hydraulic fracture length in HDR reservoir with different formation temperature

下载: 全尺寸图片幻灯片

图 6 不同压裂液排量下干热岩储层水力压裂模拟结果

Figure 6. Simulation results of hydraulic fracturing in HDR reservoir under different pumping rates of fracturing fluids

下载: 全尺寸图片幻灯片

图 7 不同天然裂缝夹角下干热岩储层水力压裂模拟结果

Figure 7. Simulation results of hydraulic fracturing in HDR reservoir with different angles of natural fractures

下载: 全尺寸图片幻灯片

图 8 不同天然裂缝密度下干热岩储层水力压裂模拟结果

Figure 8. Simulation results of hydraulic fracturing in HDR reservoir with different densities of natural fractures

下载: 全尺寸图片幻灯片

参考文献(29)

[1]	陈作,许国庆,蒋漫旗. 国内外干热岩压裂技术现状及发展建议[J]. 石油钻探技术,2019,47(6):1–8. CHEN Zuo, XU Guoqing, JIANG Manqi. The current status and development recommendations for dry hot rock fracturing technologies at home and abroad[J]. Petroleum Drilling Techniques, 2019, 47(6): 1–8.
[2]	谢文苹,路睿,张盛生,等. 青海共和盆地干热岩勘查进展及开发技术探讨[J]. 石油钻探技术,2020,48(3):77–84. XIE Wenping, LU Rui, ZHANG Shengsheng, et al. Progress in hot dry rock exploration and a discussion on development technology in the Gonghe Basin of Qinghai[J]. Petroleum Drilling Techniques, 2020, 48(3): 77–84.
[3]	BUJAKOWSKI W, BARBACKI A, MIECZNIK M, et al. Modelling geothermal and operating parameters of EGS installations in the Lower Triassic sedimentary formations of the central Poland area[J]. Renewable Energy, 2015, 80: 441–453. doi: 10.1016/j.renene.2015.02.018
[4]	思娜,叶海超,牛新明,等. 油气钻井技术在干热岩开发中的适应性分析[J]. 石油钻探技术,2019,47(4):35–40. doi: 10.11911/syztjs.2019042 SI Na, YE Haichao, NIU Xinming, et al. Analysis on the adaptability of oil and gas drilling technologies in development for hot dry rocks[J]. Petroleum Drilling Techniques, 2019, 47(4): 35–40. doi: 10.11911/syztjs.2019042
[5]	曾义金. 干热岩热能开发技术进展与思考[J]. 石油钻探技术,2015,43(2):1–7. ZENG Yijin. Technical progress and thinking for development of hot dry rock (HDR) geothermal resources[J]. Petroleum Drilling Techniques, 2015, 43(2): 1–7.
[6]	陈作,张保平,周健,等. 干热岩热储体积改造技术研究与试验[J]. 石油钻探技术,2020,48(6):82–87. CHEN Zuo, ZHANG Baoping, ZHOU Jian, et al. Research and test on the stimulated reservoir volume technology of hot dry rock[J]. Petroleum Drilling Techniques, 2020, 48(6): 82–87.
[7]	杨英涛,温庆志,段晓飞,等. 通道压裂裂缝导流能力数值模拟研究[J]. 石油钻探技术,2016,44(6):104–110. YANG Yingtao, WEN Qingzhi, DUAN Xiaofei, et al. Numerical simulation for flow conductivity in channeling fractures[J]. Petroleum Drilling Techniques, 2016, 44(6): 104–110.
[8]	KOLAWOLE O, ISPAS I. Interaction between hydraulic fractures and natural fractures: current status and prospective directions[J]. Journal of Petroleum Exploration and Production Technology, 2020, 10(4): 1613–1634. doi: 10.1007/s13202-019-00778-3
[9]	WARPINSKI N R, TEUFEL L W. Influence of geologic discontinuities on hydraulic fracture propagation (includes associated papers 17011 and 17074 )[J]. Journal of Petroleum Technology, 1987, 39(2): 209–220. doi: 10.2118/13224-PA
[10]	朱万成,魏晨慧,田军,等. 岩石损伤过程中的热–流–力耦合模型及其应用初探[J]. 岩土力学,2009,30(12):3851–3857. ZHU Wancheng, WEI Chenhui, TIAN Jun, et al. Coupled thermal-hydraulic-mechanical model during rock damage and its preliminary application[J]. Rock and Soil Mechanics, 2009, 30(12): 3851–3857.
[11]	GHASSEMI A, ZHANG Q. A transient fictitious stress boundary element method for porothermoelastic media[J]. Engineering Analysis with Boundary Elements, 2004, 28(11): 1363–1373. doi: 10.1016/j.enganabound.2004.05.003
[12]	金泰宇. 三维粗糙裂缝网络钻井液漏失流固耦合模型研究[J]. 石油钻探技术,2024,52(1):69–77. JIN Taiyu. Study on three-dimensional fluid-solid coupling model of drilling fluid leakage in rough fracture network[J]. Petroleum Drilling Techniques, 2024, 52(1): 69–77.
[13]	张文平,许争鸣,吕泽昊,等. 深层页岩欠平衡钻井气液固三相瞬态流动传热模型研究[J]. 石油钻探技术,2023,51(5):96–105. ZHANG Wenping, XU Zhengming, LYU Zehao, et al. Research on a transient flow heat transfer model of gas-liquid-solid three-phase flow for unbalanced drilling in deep shale wells[J]. Petroleum Drilling Techniques, 2023, 51(5): 96–105.
[14]	李浩,徐怀民,王千军,等. 准东地区平地泉组微量元素地球化学特征及油气地质意义[J]. 断块油气田,2023,30(2):277–285. LI Hao, XU Huaimin, WANG Qianjun, et al. Geochemical characteristics and petroleum geological significance of trace elements of Pingdiquan Formation in eastern Junggar Basin[J]. Fault-Block Oil & Gas Field, 2023, 30(2): 277–285.
[15]	HUANG Liuke, HE Rui, YANG Zhaozhong, et al. Exploring hydraulic fracture behavior in glutenite formation with strong heterogeneity and variable lithology based on DEM simulation[J]. Engineering Fracture Mechanics, 2023, 278: 109020. doi: 10.1016/j.engfracmech.2022.109020
[16]	彭井宏,周军,胡承强,等. 热力耦合作用下地下盐岩储气库注采运行稳定性研究[J]. 断块油气田,2023,30(5):858–867. PENG Jinghong, ZHOU Jun, HU Chengqiang, et al. Study on the injection-production stability of underground salt rock gas storage under thermo-mechanical coupling[J]. Fault-Block Oil & Gas Field, 2023, 30(5): 858–867.
[17]	李文拓,罗鸣,黄洪林,等. 高温高压小井眼水平井环空ECD综合计算模型[J]. 石油钻采工艺,2023,45(3):259–268. LI Wentuo, LUO Ming, HUANG Honglin, et al. Comprehensive calculation model of annular ECD for high-temperature high-pressure slim-hole horizontal wells[J]. Oil Drilling & Production Technology, 2023, 45(3): 259–268.
[18]	HOFMANN H, BABADAGLI T, ZIMMERMANN G. Hot water generation for oil sands processing from enhanced geothermal systems: process simulation for different hydraulic fracturing scenarios[J]. Applied Energy, 2014, 113: 524–547. doi: 10.1016/j.apenergy.2013.07.060
[19]	ZHANG Yanjun, GUO Liangliang, LI Zhengwei, et al. Electricity generation and heating potential from enhanced geothermal system in Songliao Basin, China: different reservoir stimulation strategies for tight rock and naturally fractured formations[J]. Energy, 2015, 93(Part 2): 1860-1885.
[20]	谢紫霄,黄中伟,熊建华,等. 天然裂缝对干热岩水力压裂裂缝扩展的影响规律[J]. 天然气工业,2022,42(4):63–72. XIE Zixiao, HUANG Zhongwei, XIONG Jianhua, et al. Influence of natural fractures on the propagation of hydraulic fractures in hot dry rock[J]. Natural Gas Industry, 2022, 42(4): 63–72.
[21]	易良平,胡滨,李小刚,等. 基于相场法的煤砂互层水力裂缝纵向延伸计算模型[J]. 煤炭学报,2020,45(增刊2):706–716. YI Liangping, HU Bin, LI Xiaogang, et al. Calculation model of hydraulic crack vertical propagation in coal-sand interbedded formation based on the phase field method[J]. Journal of China Coal Society, 2020, 45(supplement 2): 706–716.
[22]	路千里. 基于相场法的页岩水力压裂裂缝扩展规律研究[D]. 成都:西南石油大学,2018. LU Qianli. Fracture propagation modelling of shale hydraulic fracturing based on phase field method[D]. Chengdu: Southwest Petroleum University, 2018.
[23]	侯冰,戴一凡,范濛,等. 基于相场法的酸压裂缝连通孔洞数值模拟[J]. 石油学报,2022,43(6):849–859. HOU Bing, DAI Yifan, FAN Meng, et al. Numerical simulation of pores connection by acid fracturing based on phase field method[J]. Acta Petrolei Sinica, 2022, 43(6): 849–859.
[24]	NOII N, FAN Meng, WICK T, et al. A quasi-monolithic phase-field description for orthotropic anisotropic fracture with adaptive mesh refinement and primal–dual active set method[J]. Engineering Fracture Mechanics, 2021, 258: 108060. doi: 10.1016/j.engfracmech.2021.108060
[25]	NOII N, WICK T. A phase-field description for pressurized and non-isothermal propagating fractures[J]. Computer Methods in Applied Mechanics and Engineering, 2019, 351: 860–890. doi: 10.1016/j.cma.2019.03.058
[26]	HEISTER T, WICK T. Pfm-cracks: A parallel-adaptive framework for phase-field fracture propagation[J]. Software Impacts, 2020, 6: 100045. doi: 10.1016/j.simpa.2020.100045
[27]	TRANTER C J. The opening of a pair of coplanar Griffith cracks under internal pressure[J]. The Quarterly Journal of Mechanics and Applied Mathematics, 1961, 14(3): 283–292. doi: 10.1093/qjmam/14.3.283
[28]	周舟,金衍,卢运虎,等. 干热岩地热储层钻井和水力压裂工程技术难题和攻关建议[J]. 中国科学:物理学力学天文学,2018,48(12):124621. ZHOU Zhou, JIN Yan, LU Yunhu, et al. Present challenge and prospects of drilling and hydraulic fracturing technology for hot dry rock geothermal reservoir[J]. Scientia Sinica: Physica, Mechanica & Astronomica, 2018, 48(12): 124621.
[29]	周舟,金衍,曾义金,等. 青海共和盆地干热岩地热储层水力压裂物理模拟和裂缝起裂与扩展形态研究[J]. 吉林大学学报(地球科学版),2019,49(5):1425–1430. ZHOU Zhou, JIN Yan, ZENG Yijin, et al. Experimental study on hydraulic fracturing physics simulation, crack initiation and propagation in hot dry rock geothermal reservoir in Gonghe Basin, Qinghai[J]. Journal of Jilin University(Earth Science Edition), 2019, 49(5): 1425–1430.

施引文献(1)

期刊类型引用(1)

侯冰，廖志豪，张庄，罗加伦，琚宜文，王文. 水力压裂物理模拟方法的数字化和智能化发展综述. 辽宁石油化工大学学报. 2025(02): 1-12 .

百度学术

其他类型引用(0)

资源附件(0)

图(8)

计量

文章访问数: 190
HTML全文浏览量: 61
PDF下载量: 89
被引次数: 1

1. 深层干热岩储层水力压裂数值模型
1.1 控制方程
1.1.1 相场−位移控制方程
1.1.2 温度场控制方程
1.1.3 流体压力控制方程
1.2 自适应网格
1.3 多物理场耦合求解
1.4 初始条件与边界条件
1.5 解析解验证
2. 地质与工程因素对裂缝扩展的影响
2.1 温度
2.2 压裂液排量
2.3 天然裂缝夹角
2.4 多天然裂缝干扰下裂缝扩展
3. 结　论

1. 深层干热岩储层水力压裂数值模型
1.1 控制方程
1.1.1 相场−位移控制方程
1.1.2 温度场控制方程
1.1.3 流体压力控制方程
1.2 自适应网格
1.3 多物理场耦合求解
1.4 初始条件与边界条件
1.5 解析解验证
2. 地质与工程因素对裂缝扩展的影响
2.1 温度
2.2 压裂液排量
2.3 天然裂缝夹角
2.4 多天然裂缝干扰下裂缝扩展
3. 结　论

参考文献(29)

施引文献

资源附件(0)

基于相场法的深层干热岩储层水力压裂模拟研究

作者简介: 戴一凡（1996—），男，江苏泰州人，2018年毕业于中国石油大学（华东）石油工程专业，中国石油大学（北京）石油与天然气工程专业在读博士研究生，主要从事岩石力学方面的研究。E-mail: 819673313@qq.com

通讯作者: 侯冰，binghou@vip.163.com

计量

出版历程