CN115828683A - 一种含天然裂缝地层裂缝缝宽的计算方法 - Google Patents
一种含天然裂缝地层裂缝缝宽的计算方法 Download PDFInfo
- Publication number
- CN115828683A CN115828683A CN202211502769.4A CN202211502769A CN115828683A CN 115828683 A CN115828683 A CN 115828683A CN 202211502769 A CN202211502769 A CN 202211502769A CN 115828683 A CN115828683 A CN 115828683A
- Authority
- CN
- China
- Prior art keywords
- fracture
- unit
- crack
- calculating
- width
- 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
Images
Abstract
本发明涉及一种含天然裂缝地层裂缝缝宽的计算方法,属于水力压裂数值模拟技术领域。为了克服现有技术中的问题,本发明提供一种含天然裂缝地层裂缝缝宽的计算方法,通过获取储层地质力学参数和储层岩石力学参数建立含天然裂缝储层有限元物理模型,建立含天然裂缝储层的裂缝起裂和扩展数学模型,开展水力压裂数值模拟,获取裂缝扩展形态,基于裂缝扩展变形后有限元物理模型,获取水力裂缝单元的坐标信息并划分裂缝边,求取裂缝宽度。本发明克服传统基于裂缝图像识别由于像素点不连续导致误差的缺点。
Description
技术领域
本发明涉及一种含天然裂缝地层裂缝缝宽的计算方法,属于水力压裂数值模拟技术领域。
背景技术
利用有限元方法模拟计算水力裂缝参数是一种常见的数值模拟手段,对于含天然裂缝地层,采用有限元方法模拟裂缝扩展后裂缝形态复杂,对于形成的裂缝可能存在主裂缝、一级裂缝甚至二级裂缝多级裂缝体系,导致裂缝宽度统计和计算困难,目前对于裂缝宽度的提取是通过在两侧壁面上标记两点,通过获取两点间的距离,来反映裂缝宽度的大小,但是该方法存在由于像素点不连续导致误差的缺点。
发明内容
为了克服现有技术中的问题,本发明提供一种含天然裂缝地层裂缝缝宽的计算方法。
本发明解决上述技术问题所提供的技术方案是:一种含天然裂缝地层裂缝缝宽的计算方法,包括以下步骤:
步骤一、:获取目标储层的地质力学参数以及目标储层岩石力学的基本资料,所述基本资料包括储层地应力大小、地应力方向、天然裂缝特征、储层岩石杨氏模量、泊松比、抗拉强度;
步骤二、根据步骤一获取的参数,建立含天然裂缝的有限元物理模型;
步骤三、步骤一和步骤二建立的模型,建立含天然裂缝储层的裂缝起裂和扩展数学模型;
步骤四、根据步骤二和步骤三建立的模型,开展裂缝扩展数值模拟分析,获取裂缝扩展形态;
步骤五、基于裂缝扩展变形后有限元物理模型,获取水力裂缝单元的坐标信息并划分裂缝边;
步骤六、基于裂缝边求取裂缝宽度。
进一步的技术方案是,所述步骤1目标储层的地质力学参数包括地应力大小、地应力方位、天然裂缝长度和方向。
进一步的技术方案是,所述步骤1目标储层的岩石力学参数包括杨氏模量、泊松比、抗拉强度,通过岩石力学实验获取。
进一步的技术方案是,所述步骤3含天然裂缝储层的裂缝起裂和扩展数学模型包括:
(1)裂缝单元应力场模型:
式中:σ为应力矢量,Pa;σn为裂缝单元法向应力,Pa;σs、σt为裂缝单元切向应力,Pa;ε为应变,无因次;K为单元刚度矩阵,无因次;dn为裂缝单元法向位移,m;ds、dt裂缝单元切向位移,m;To为裂缝单元厚度,m;
(2)裂缝单元流体流动模型:
式中:q为通过裂缝单元单位面积的体积流率密度矢量,m2/s;d为裂缝单元的开度,m;μ为压裂液粘度,Pa·s;▽p为单元切向流动流体压力梯度,Pa/m;qt为流体流入裂缝单元上表面的流量,m3/s;qb为流体流入裂缝单元下表面的流量,m3/s;ct为裂缝单元上表面的流量滤失系数,m3/(Pa·s);cb为裂缝单元下表面的流量滤失系数,m3/(Pa·s);pt为裂缝单元上表面的流体压力,MPa;pb为裂缝单元下表面的流体压力,MPa;pi为裂缝单元上表面的流体压力,MPa。
(3)裂缝起裂准则:
(4)裂缝扩展准则:
GT=Gn+Gs+Gt
GT=GC时,裂缝扩展
式中:Gn为裂缝单元法向能量释放率,Pa·m;Gs、Gt为裂缝单元切向能量释放率,Pa·m;为裂缝单元法向临界能量释放率,Pa·m;为裂缝单元切向临界释放率,Pa·m;η为材料本身有关的常数,取2.284,无因次;GC为裂缝单元混合模式下的总临界能量释放率,Pa·m。
进一步的技术方案是,所述步骤5中获取水力裂缝单元的坐标信息并划分裂缝边的具体过程为:提取裂缝单元法向位移大于零的裂缝单元,确定该裂缝单元的节点分别为Node1、Node2、Node3和Node4。
进一步的技术方案是,所述步骤6中求取裂缝宽度的具体过程为:分别计算Node1、Node2之间宽度,Node3、Node4之间有效宽度,根据流体流动规律,裂缝单元入口端流体压力大于出口端流体压力,选取节点间最大值作为裂缝宽度
本发明具有以下有益效果:本发明克服传统基于裂缝图像识别由于像素点不连续导致误差的缺点。
附图说明
图1是本发明的流程框图;
图2是节点示意图;
图3是裂缝识别图;
图4是裂缝宽度分布图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明的一种含天然裂缝地层裂缝缝宽的计算方法,包括以下步骤:
步骤一、获取目标储层的地质力学参数以及目标储层岩石力学的基本资料,所述基本资料包括储层地应力大小、地应力方向、天然裂缝特征、储层岩石杨氏模量、泊松比、抗拉强度;
步骤二、根据步骤一获取的参数,建立含天然裂缝的有限元物理模型;
步骤三、步骤一和步骤二建立的模型,建立含天然裂缝储层的裂缝起裂和扩展数学模型;
含天然裂缝储层的裂缝起裂和扩展数学模型包括:
(1)裂缝单元应力场模型:
式中:σ为应力矢量,Pa;σn为裂缝单元法向应力,Pa;σs、σt为裂缝单元切向应力,Pa;ε为应变,无因次;K为单元刚度矩阵,无因次;dn为裂缝单元法向位移,m;ds、dt裂缝单元切向位移,m;To为裂缝单元厚度,m;
(2)裂缝单元流体流动模型:
式中:q为通过裂缝单元单位面积的体积流率密度矢量,m2/s;d为裂缝单元的开度,m;μ为压裂液粘度,Pa·s;▽p为单元切向流动流体压力梯度,Pa/m;qt为流体流入裂缝单元上表面的流量,m3/s;qb为流体流入裂缝单元下表面的流量,m3/s;ct为裂缝单元上表面的流量滤失系数,m3/(Pa·s);cb为裂缝单元下表面的流量滤失系数,m3/(Pa·s);pt为裂缝单元上表面的流体压力,MPa;pb为裂缝单元下表面的流体压力,MPa;pi为裂缝单元上表面的流体压力,MPa。
(3)裂缝起裂准则:
(4)裂缝扩展准则:
GT=Gn+Gs+Gt
GT=GC时,裂缝扩展
式中:Gn为裂缝单元法向能量释放率,Pa·m;Gs、Gt为裂缝单元切向能量释放率,Pa·m;为裂缝单元法向临界能量释放率,Pa·m;为裂缝单元切向临界释放率,Pa·m;η为材料本身有关的常数,取2.284,无因次;GC为裂缝单元混合模式下的总临界能量释放率,Pa·m;
步骤四、根据步骤二和步骤三建立的模型,开展裂缝扩展数值模拟分析,获取裂缝扩展形态;
步骤五、基于裂缝扩展变形后有限元物理模型,获取水力裂缝单元的坐标信息并划分裂缝边;
步骤六、基于裂缝边求取裂缝宽度。
以上所述,并非对本发明作任何形式上的限制,虽然本发明已通过上述实施例揭示,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (6)
1.一种含天然裂缝地层裂缝缝宽的计算方法,其特征在于,包括以下步骤;
步骤一、:获取目标储层的地质力学参数以及目标储层岩石力学的基本资料,所述基本资料包括储层地应力大小、地应力方向、天然裂缝特征、储层岩石杨氏模量、泊松比、抗拉强度;
步骤二、根据步骤一获取的参数,建立含天然裂缝的有限元物理模型;
步骤三、步骤一和步骤二建立的模型,建立含天然裂缝储层的裂缝起裂和扩展数学模型;
步骤四、根据步骤二和步骤三建立的模型,开展裂缝扩展数值模拟分析,获取裂缝扩展形态;
步骤五、基于裂缝扩展变形后有限元物理模型,获取水力裂缝单元的坐标信息并划分裂缝边;
步骤六、基于裂缝边求取裂缝宽度。
2.根据权利要求1所述的一种含天然裂缝地层裂缝缝宽的计算方法,其特征在于,所述步骤1目标储层的地质力学参数包括地应力大小、地应力方位、天然裂缝长度和方向。
3.根据权利要求1所述的一种含天然裂缝地层裂缝缝宽的计算方法,其特征在于,所述步骤1目标储层的岩石力学参数包括杨氏模量、泊松比、抗拉强度,通过岩石力学实验获取。
4.根据权利要求1所述的一种含天然裂缝地层裂缝缝宽的计算方法,其特征在于,所述步骤3含天然裂缝储层的裂缝起裂和扩展数学模型包括:
(1)裂缝单元应力场模型:
式中:σ为应力矢量,Pa;σn为裂缝单元法向应力,Pa;σs、σt为裂缝单元切向应力,Pa;ε为应变,无因次;K为单元刚度矩阵,无因次;dn为裂缝单元法向位移,m;ds、dt裂缝单元切向位移,m;To为裂缝单元厚度,m;
(2)裂缝单元流体流动模型:
式中:q为通过裂缝单元单位面积的体积流率密度矢量,m2/s;d为裂缝单元的开度,m;μ为压裂液粘度,Pa·s;为单元切向流动流体压力梯度,Pa/m;qt为流体流入裂缝单元上表面的流量,m3/s;qb为流体流入裂缝单元下表面的流量,m3/s;ct为裂缝单元上表面的流量滤失系数,m3/(Pa·s);cb为裂缝单元下表面的流量滤失系数,m3/(Pa·s);pt为裂缝单元上表面的流体压力,MPa;pb为裂缝单元下表面的流体压力,MPa;pi为裂缝单元上表面的流体压力,MPa;
(3)裂缝起裂准则:
(4)裂缝扩展准则:
GT=Gn+Gs+Gt
GT=GC时,裂缝扩展
式中:Gn为裂缝单元法向能量释放率,Pa·m;Gs、Gt为裂缝单元切向能量释放率,Pa·m;Gn C为裂缝单元法向临界能量释放率,Pa·m;Gs C为裂缝单元切向临界释放率,Pa·m;η为材料本身有关的常数,取2.284,无因次;GC为裂缝单元混合模式下的总临界能量释放率,Pa·m。
5.根据权利要求1所述的一种含天然裂缝地层裂缝缝宽的计算方法,其特征在于,所述步骤5中获取水力裂缝单元的坐标信息并划分裂缝边的具体过程为:提取裂缝单元法向位移大于零的裂缝单元,确定该裂缝单元的节点分别为Node1、Node2、Node3和Node4。
6.根据权利要求1所述的一种含天然裂缝地层裂缝缝宽的计算方法,其特征在于,所述步骤6中求取裂缝宽度的具体过程为:分别计算Node1、Node2之间宽度,Node3、Node4之间有效宽度,根据流体流动规律,裂缝单元入口端流体压力大于出口端流体压力,选取节点间最大值作为裂缝宽度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211502769.4A CN115828683A (zh) | 2022-11-28 | 2022-11-28 | 一种含天然裂缝地层裂缝缝宽的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211502769.4A CN115828683A (zh) | 2022-11-28 | 2022-11-28 | 一种含天然裂缝地层裂缝缝宽的计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115828683A true CN115828683A (zh) | 2023-03-21 |
Family
ID=85532246
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211502769.4A Pending CN115828683A (zh) | 2022-11-28 | 2022-11-28 | 一种含天然裂缝地层裂缝缝宽的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115828683A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116956769A (zh) * | 2023-06-30 | 2023-10-27 | 中国地质大学(北京) | 随机裂隙地层中水力压裂裂纹扩展的数值模拟方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109359376A (zh) * | 2018-10-10 | 2019-02-19 | 北京科技大学 | 水力压裂裂缝在页岩储层天然裂缝界面扩展判识方法 |
CN112012710A (zh) * | 2019-05-31 | 2020-12-01 | 中国石油化工股份有限公司 | 一种砂泥岩交互地层水平井压裂三维裂缝扩展模拟方法 |
CN112036098A (zh) * | 2020-09-15 | 2020-12-04 | 中国石油大学(华东) | 一种深层油气藏水力裂缝扩展数值模拟的方法 |
CN114722682A (zh) * | 2022-05-05 | 2022-07-08 | 西昌学院 | 一种页岩储层水平井暂堵压裂多裂缝竞争起裂预测方法 |
-
2022
- 2022-11-28 CN CN202211502769.4A patent/CN115828683A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109359376A (zh) * | 2018-10-10 | 2019-02-19 | 北京科技大学 | 水力压裂裂缝在页岩储层天然裂缝界面扩展判识方法 |
CN112012710A (zh) * | 2019-05-31 | 2020-12-01 | 中国石油化工股份有限公司 | 一种砂泥岩交互地层水平井压裂三维裂缝扩展模拟方法 |
CN112036098A (zh) * | 2020-09-15 | 2020-12-04 | 中国石油大学(华东) | 一种深层油气藏水力裂缝扩展数值模拟的方法 |
CN114722682A (zh) * | 2022-05-05 | 2022-07-08 | 西昌学院 | 一种页岩储层水平井暂堵压裂多裂缝竞争起裂预测方法 |
Non-Patent Citations (1)
Title |
---|
罗波: "川西须五段砂泥岩互层水力裂缝扩展规律研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》, no. 01, pages 1 - 93 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116956769A (zh) * | 2023-06-30 | 2023-10-27 | 中国地质大学(北京) | 随机裂隙地层中水力压裂裂纹扩展的数值模拟方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109117540B (zh) | 一种解决大坝混凝土力学参数反演不唯一性的概率统计分析方法 | |
US11542801B2 (en) | Optimized design method for temporary blocking agent to promote uniform expansion of fractures produced by fracturing in horizontal wells | |
Sun et al. | Investigating the effect of improved fracture conductivity on production performance of hydraulically fractured wells: Field-case studies and numerical simulations | |
CN108681635B (zh) | 一种致密储层体积压裂可压性评价方法 | |
GB2563241A (en) | Hydraulic fracturing simulation | |
CN113076676B (zh) | 非常规油气藏水平井压裂缝网扩展与生产动态耦合方法 | |
CN109064561B (zh) | 基于三维拟连续介质水力压裂模型的支撑剂运移模拟方法 | |
CN110348032B (zh) | 层理发育页岩地层水力裂缝高度数值模拟方法 | |
CN108051035A (zh) | 基于门控循环单元的神经网络模型的管网漏损识别方法 | |
CN106894814A (zh) | 复杂断块油藏高含水后期剩余油二次富集的快速识别方法 | |
CN115828683A (zh) | 一种含天然裂缝地层裂缝缝宽的计算方法 | |
Sagmo et al. | Numerical simulations of the NREL S826 airfoil | |
CN103913772A (zh) | 基于储层地质力学参数的微地震事件正演模拟方法 | |
CN108647449A (zh) | 一种基于絮凝动力学的粘性泥沙运动数值模拟方法 | |
CN109033504A (zh) | 一种油水井套管损坏预测方法 | |
CN106202649A (zh) | 考虑软基固结和混凝土徐变的闸首施工仿真方法 | |
Hewitt et al. | Obstructed and channelized viscoplastic flow in a Hele-Shaw cell | |
CN110863810B (zh) | 一种耦合页岩气藏水力压裂返排生产过程一体化模拟方法 | |
CN115510778A (zh) | 一种陆相页岩储层无限级压裂工艺优化方法及系统 | |
CN110516407B (zh) | 一种裂缝性储层水平井段内多簇压裂裂缝复杂度计算方法 | |
CN117291059A (zh) | 一种高温高压有水气藏可采储量及水侵参数计算方法 | |
CN109871507A (zh) | 正交各向异性煤层裂隙绝对渗透率计算方法 | |
CN112765854B (zh) | 一种路面内部裂缝数量预测方法 | |
CN111159794A (zh) | 多裂隙类岩石试样力学性质的几何损伤流变分析方法 | |
CN110909477B (zh) | 一种酸蚀裂缝粗糙程度定量表征方法 |
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 |