CN112668219A - 一种水合物储层水力压裂裂缝转向扩展模拟方法 - Google Patents
一种水合物储层水力压裂裂缝转向扩展模拟方法 Download PDFInfo
- Publication number
- CN112668219A CN112668219A CN202011545155.5A CN202011545155A CN112668219A CN 112668219 A CN112668219 A CN 112668219A CN 202011545155 A CN202011545155 A CN 202011545155A CN 112668219 A CN112668219 A CN 112668219A
- Authority
- CN
- China
- Prior art keywords
- setting
- reservoir
- hydraulic fracturing
- fracture
- hydrate reservoir
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000004088 simulation Methods 0.000 title claims abstract description 36
- 239000012530 fluid Substances 0.000 claims abstract description 23
- 238000009826 distribution Methods 0.000 claims abstract description 10
- 230000008569 process Effects 0.000 claims description 21
- 238000006073 displacement reaction Methods 0.000 claims description 14
- 239000011148 porous material Substances 0.000 claims description 12
- 230000035699 permeability Effects 0.000 claims description 10
- 238000002347 injection Methods 0.000 claims description 7
- 239000007924 injection Substances 0.000 claims description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000012821 model calculation Methods 0.000 claims description 3
- 238000009738 saturating Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 239000000463 material Substances 0.000 abstract description 2
- 238000012805 post-processing Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 5
- 230000005484 gravity Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 150000004677 hydrates Chemical class 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- NMJORVOYSJLJGU-UHFFFAOYSA-N methane clathrate Chemical compound C.C.C.C.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O NMJORVOYSJLJGU-UHFFFAOYSA-N 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种水合物储层水力压裂裂缝转向扩展模拟方法，包括如下步骤：步骤100：建立水合物储层几何模型；步骤200：对几何模型进行网格划分；步骤300：对模型赋值材料参数；步骤400：对物理模型中边界条件和初始条件进行设置，开展模拟；步骤500：对模拟结果后处理。本发明利用扩展有限元方法，对水合物储层中水力压裂裂缝转向扩展进行数值模拟，分析射孔角度、地应力分布、压裂液排量等参数对裂缝扩展路径的影响。
Description
技术领域
本发明涉及水合物储层水力压裂数值模拟领域，是一种利用扩展有限元发模拟水力压裂的方法。
背景技术
海域天然气水合物储层以泥质粉砂型为主，渗透率低，开采难度大。目前常用的注热、降压等开采方法，生产成本过高，难以实现水合物的商业化开采。水力压裂是一种有效的增产改造方法，掌握水合物储层中水力压裂裂缝转向的扩展规律，可为海域水合物的高效开发提供理论支撑。
发明内容
本发明的目的在于提供一种水合物储层水力压裂裂缝转向扩展模拟方法，以模拟不同地应力分布、射孔方位角和压裂液排量等参数条件下水力裂缝转向扩展行为，明确各参数对裂缝转向扩展程度的影响，掌握裂缝转向扩展路径，有利于深入理解水力裂缝在水合物储层中的扩展规律。
为达到上述目的，本发明具体提供下述技术方案：
一种水合物储层水力压裂裂缝转向扩展模拟方法，包括如下步骤：
步骤100：建立几何模型，包括初始裂缝；
步骤200：对几何模型进行网格划分；
步骤300：对物理模型参数进行设置；
步骤400：对物理模型中边界条件和初始条件进行设置，开展模拟；
步骤500：对模拟结果进行后处理，后处理过程包括查看不同射孔角度、不同压裂液排量和不同地应力分布情况下的裂缝转向扩展情况。
作为本发明的一种优选方案，所述几何模型包括储层模型和裂缝模型，其中，步骤100的具体设置过程为：
建立储层物理模型为正方形，边长为50m，在中心位置设置初始对称裂缝，缝长为1m；
对裂缝进行初始角度设定，角度为与最大主应力方向的夹角，分别设置为30°、45°和60°。
作为本发明的一种优选方案，步骤200的具体设置过程为：
对整个储层进行四边形网格划分，网格边长尺寸为0.1m。
作为本发明的一种优选方案，步骤300的具体设置过程为：
对储层设置基本物性参数，包括弹性模量、泊松比、孔隙度、渗透率、饱和度、抗拉强度、滤失率等；
对压裂液设置基本物性参数，包括粘度、注入储层的排量。
作为本发明的一种优选方案，步骤400中具体设置过程为：
(a)采用的模型计算方法为：
(1)平衡方程：
多孔介质的应力平衡可用虚功原理表示，储层介质的虚功与作用在该介质上作用力所产生的虚功相等，即；
其中，t—面力，N；
f—体力，N；
δε—虚位移；
δu—虚应变；
u—虚应变；
Dep—弹塑性矩阵；
(2)连续方程：
其中，k′—初始渗透率张量和液体密度的乘积；
kr—比例渗透系数，dless；
Kw—水的体积模量，MPa；
φ一孔隙度，dless；
g一重力加速度，m/s2；
(3)有限元方程：
定义形函数
其中，
—单元节点的位移，m；
将形函数带入平衡方程，得：
其中
k＝∫VBTDepBdV (5)
df＝∫VNu TdfdV+∫sNu TdtdS (7)
流量边界条件：
其中，n—流量边界的单位法向量；
qw—单位时间流过边界的水流量，m3/d；
(b)设定边界条件和初始条件
对储层四条边界设置垂向位移为0，边界孔压为定压0MPa，整个储层的孔隙被流体饱和，初始孔压为0；
(c)设定压裂液注入时长
设定压裂液的注入时间，定排量注入；
(d)提交模型，开展模拟
作为本发明的一种优选方案，步骤500的具体设置过程为：
对模拟结果进行后处理，后处理过程包括提取水合物储层的水力压裂裂缝扩展路径，查看不同射孔角度、不同压裂液排量和不同地应力分布情况下的裂缝转向扩展情况。
本发明与现有技术相比较具有如下有益效果：
本发明利用扩展有限元方法模拟了水合物储层的水力压裂裂缝扩展模拟，可以真实反映裂缝的扩展路径，对不同参数对水合物储层水力裂缝的扩展情况进行了模拟，对各影响因素进行分析，确定了各参数对裂缝扩展的影响程度。
附图说明
为了更清楚地说明本发明的实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地，下面描述中的附图仅仅是示例性的，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图引伸获得其它的实施附图。
图1为本发明实施例中的流程示意图；
图2本发明实施例中的水合物储层水力压裂裂缝转向扩展模拟结果图；
图3本发明实施例中提取模拟结果的裂缝扩展路径；
图4本发明实施例中不同排量不同射孔角度不同地应力分布条件下的裂缝扩展模拟结果；
图5本发明实施例中的裂缝扩展对比图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，一种水合物储层水力压裂裂缝转向扩展模拟方法，包括如下步骤：
步骤100：建立几何模型，包括初始裂缝；
步骤200：对几何模型进行网格划分；
步骤300：对物理模型参数进行设置；
步骤400：对物理模型中边界条件和初始条件进行设置，开展模拟；
步骤500：对模拟结果进行后处理，后处理过程包括查看不同射孔角度、不同压裂液排量和不同地应力分布情况下的裂缝转向扩展情况。
所述几何模型包括储层模型和裂缝模型，其中，步骤100的具体设置过程为：
建立储层物理模型为正方形，边长为50m，在中心位置设置初始对称裂缝，缝长为1m；
对裂缝进行初始角度设定，角度为与最大主应力方向的夹角，分别设置为30°、45°和60°。
步骤200的具体设置过程为：
对整个储层进行四边形网格划分，网格边长尺寸为0.1m。
步骤300的具体设置过程为：
对储层设置基本物性参数，包括弹性模量、泊松比、孔隙度、渗透率、饱和度、抗拉强度、滤失率等；
对压裂液设置基本物性参数，包括粘度、注入储层的排量。
步骤400中具体设置过程为：
(a)采用的模型计算方法为：
(1)平衡方程：
多孔介质的应力平衡可用虚功原理表示，储层介质的虚功与作用在该介质上作用力所产生的虚功相等，即；
其中，t—面力，N；
f—体力，N；
δε—虚位移；
δu—虚应变；
u—虚应变；
Dep—弹塑性矩阵；
(2)连续方程：
其中，k′—初始渗透率张量和液体密度的乘积；
kr—比例渗透系数，dless；
Kw—水的体积模量，MPa；
φ一孔隙度，dless；
g一重力加速度，m/s2；
(3)有限元方程：
定义形函数
其中，
—单元节点的位移，m；
将形函数带入平衡方程，得：
其中
k＝∫VBTDepBdV (5)
df＝∫VNu TdfdV+∫sNu TdtdS (7)
流量边界条件：
其中，n—流量边界的单位法向量；
qe—单位时间流过边界的水流量，m3/d；
(b)设定边界条件和初始条件
对储层四条边界设置垂向位移为0，边界孔压为定压0MPa，整个储层的孔隙被流体饱和，初始孔压为0；
(c)设定压裂液注入时长
设定压裂液的注入时间，定排量注入；
(d)提交模型，开展模拟。
步骤500的具体设置过程为：
对模拟结果进行后处理，后处理过程包括提取水合物储层的水力压裂裂缝扩展路径，查看不同射孔角度、不同压裂液排量和不同地应力分布情况下的裂缝转向扩展情况。
图2为水合物储层水力压裂裂缝转向扩展模拟结果图，图3为提取模拟结果的裂缝扩展路径，图4为不同排量不同射孔角度不同地应力分布条件下的裂缝扩展模拟结果，图5为裂缝扩展对比图，通过对其中的裂缝进行提取，并将多条扩展裂缝进行对比而形成。
如下表所示，下表为本实施例中得到个附图结构的参数设定。
表1模拟实例中所使用的参数
以上实施例仅为本申请的示例性实施例，不用于限制本申请，本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内，对本申请做出各种修改或等同替换，这种修改或等同替换也应视为落在本申请的保护范围内。
Claims (6)
1.一种水合物储层水力压裂裂缝转向扩展模拟方法，其特征在于，包括如下步骤：
步骤100：建立几何模型，包括初始裂缝；
步骤200：对几何模型进行网格划分；
步骤300：对物理模型参数进行设置；
步骤400：对物理模型中边界条件和初始条件进行设置，开展模拟；
步骤500：对模拟结果进行后处理，后处理过程包括查看不同射孔角度、不同压裂液排量和不同地应力分布情况下的裂缝转向扩展情况。
2.根据权利要求1所述的一种水合物储层水力压裂裂缝转向扩展模拟方法，其特征在于：所述几何模型，其中，步骤100的具体设置过程为：
建立储层物理模型为正方形，边长为50m，在中心位置设置初始对称裂缝，缝长为1m；
对裂缝进行初始角度设定，角度为与最大主应力方向的夹角，分别设置为30°、45°和60°。
3.根据权利要求1所述的一种水合物储层水力压裂裂缝转向扩展模拟方法，其特征在于：对几何模型进行网格划分，步骤200的具体设置过程为：
对整个储层进行四边形网格划分，网格边长尺寸为0.1m。
4.根据权利要求1所述的一种水合物储层水力压裂裂缝转向扩展模拟方法，其特征在于：步骤300的具体设置过程为：
对储层设置基本物性参数，包括弹性模量、泊松比、孔隙度、渗透率、饱和度、抗拉强度、滤失率等；
对压裂液设置基本物性参数，包括粘度、注入储层的排量。
5.根据权利要求1所述的一种水合物储层水力压裂裂缝转向扩展模拟方法，其特征在于：步骤400中具体设置过程为：
(a)采用的模型计算方法为：
(1)平衡方程：
多孔介质的应力平衡可用虚功原理表示，储层介质的虚功与作用在该介质上作用力所产生的虚功相等，即；
其中，t-面力，N；
f-体力，N；
δε-虚位移；
δu-虚应变；
u-虚应变；
Dep-弹塑性矩阵
(2)连续方程：
其中，k′-初始渗透率张量和液体密度的乘积；
kr-比例渗透系数，dless；
Kw-水的体积模量，MPa；
φ一孔隙度，dless；
g一重力加速度，m/s2；
(3)有限元方程：
定义形函数
其中，
-单元节点的位移，m；
将形函数带入平衡方程，得：
其中
k＝∫VBTDepBdV (5)
df＝∫VNu TdfdV+∫sNu TdtdS (7)
流量边界条件：
其中，n-流量边界的单位法向量；
qw-单位时间流过边界的水流量，m3/d；
(b)设定边界条件和初始条件
对储层四条边界设置垂向位移为0，边界孔压为定压0MPa，整个储层的孔隙被流体饱和，初始孔压为0；
(c)设定压裂液注入时长
设定压裂液的注入时间，定排量注入；
(d)提交模型，开展模拟。
6.根据权利要求5所述的一种水合物储层水力压裂裂缝转向扩展模拟方法，其特征在于：步骤500中具体设置过程为：
对模拟结果进行后处理，后处理过程包括提取水合物储层的水力压裂裂缝扩展路径，查看不同射孔角度、不同压裂液排量和不同地应力分布情况下的裂缝转向扩展情况。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011545155.5A CN112668219A (zh) | 2020-12-24 | 2020-12-24 | 一种水合物储层水力压裂裂缝转向扩展模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011545155.5A CN112668219A (zh) | 2020-12-24 | 2020-12-24 | 一种水合物储层水力压裂裂缝转向扩展模拟方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112668219A true CN112668219A (zh) | 2021-04-16 |
Family
ID=75408193
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011545155.5A Pending CN112668219A (zh) | 2020-12-24 | 2020-12-24 | 一种水合物储层水力压裂裂缝转向扩展模拟方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112668219A (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6076046A (en) * | 1998-07-24 | 2000-06-13 | Schlumberger Technology Corporation | Post-closure analysis in hydraulic fracturing |
| US20150160371A1 (en) * | 2013-12-06 | 2015-06-11 | Schlumberger Technology Corporation | Gpu accelerated deflation in geomechanics simulator |
| CN105114065A (zh) * | 2015-08-06 | 2015-12-02 | 中国石油天然气股份有限公司 | 水力压裂裂缝延伸过程的模拟方法 |
| US20180355701A1 (en) * | 2017-06-07 | 2018-12-13 | Rockfield Software Limited | Hydraulic fracturing simulation |
-
2020
- 2020-12-24 CN CN202011545155.5A patent/CN112668219A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6076046A (en) * | 1998-07-24 | 2000-06-13 | Schlumberger Technology Corporation | Post-closure analysis in hydraulic fracturing |
| US20150160371A1 (en) * | 2013-12-06 | 2015-06-11 | Schlumberger Technology Corporation | Gpu accelerated deflation in geomechanics simulator |
| CN105114065A (zh) * | 2015-08-06 | 2015-12-02 | 中国石油天然气股份有限公司 | 水力压裂裂缝延伸过程的模拟方法 |
| US20180355701A1 (en) * | 2017-06-07 | 2018-12-13 | Rockfield Software Limited | Hydraulic fracturing simulation |
Non-Patent Citations (2)
| Title |
|---|
| 刘祖岩等: "材料变形力学", vol. 1, 哈尔滨工业大学出版社, pages: 95 - 96 * |
| 郭亚兵: "致密砂岩气藏暂堵转向压裂技术研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》, no. 03, 15 March 2017 (2017-03-15), pages 019 - 226 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109902918B (zh) | 一种页岩压裂自支撑裂缝缝网渗透率计算方法 | |
| CN107622165A (zh) | 一种页岩气水平井重复压裂产能计算方法 | |
| CN105404735B (zh) | 裂缝与基质对超低渗透油藏单井产量贡献率的定量评价方法 | |
| CN106382111B (zh) | 增加页岩气压裂裂缝复杂性的方法 | |
| CN110348032B (zh) | 层理发育页岩地层水力裂缝高度数值模拟方法 | |
| CN109751025B (zh) | 一种提高深层页岩气全尺度裂缝支撑体积的压裂方法 | |
| CN112541287A (zh) | 疏松砂岩压裂充填防砂增产调剖一体化设计方法 | |
| CN106372377B (zh) | 一种细粉砂油层充填防砂方法 | |
| CN107589020A (zh) | 一种基于应力路径的水力压裂试验方法 | |
| CN111577269A (zh) | 一种基于离散元流固耦合的多簇压裂裂缝形态预测方法 | |
| Fang et al. | Research of phenolic crosslinker gel for profile control and oil displacement in high temperature and high salinity reservoirs | |
| Li et al. | Investigation of sand production mechanisms using DEM with fluid flow | |
| CN113536706A (zh) | 一种页岩油储层水平井分簇设计方法 | |
| CN111101913A (zh) | 基于离散元的砂砾岩水力压裂裂缝穿砾过程描述方法 | |
| Bauer et al. | Impact of acidizing-induced wormholes in chalk on rock strength | |
| CN112100707A (zh) | 一种穿层压裂图版构建方法 | |
| CN112668219A (zh) | 一种水合物储层水力压裂裂缝转向扩展模拟方法 | |
| CN110374574A (zh) | 直井定向井大规模压裂控制裂缝复杂程度的方法 | |
| CN110991016B (zh) | 不规则边界油藏两口体积压裂水平井渗流模型的建立方法 | |
| CN112012710A (zh) | 一种砂泥岩交互地层水平井压裂三维裂缝扩展模拟方法 | |
| CN115935843A (zh) | 一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法 | |
| CN113738351B (zh) | 一种断缝体油藏物理模型的制作方法和实验方法 | |
| CN115099062A (zh) | 一种致密油藏蓄能压裂工艺参数设计方法 | |
| CN110850470A (zh) | 一种利用地震资料计算地应力差异系数的方法 | |
| Edwards et al. | Field normalization of formation mechanical properties for use in sand control management |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |