CN108829994B - 多层页岩拟三维压裂模型缝高求解方法 - Google Patents

多层页岩拟三维压裂模型缝高求解方法 Download PDF

Info

Publication number
CN108829994B
CN108829994B CN201810655968.6A CN201810655968A CN108829994B CN 108829994 B CN108829994 B CN 108829994B CN 201810655968 A CN201810655968 A CN 201810655968A CN 108829994 B CN108829994 B CN 108829994B
Authority
CN
China
Prior art keywords
fracture
net
perforation
fracturing
height
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.)
Expired - Fee Related
Application number
CN201810655968.6A
Other languages
English (en)
Other versions
CN108829994A (zh
Inventor
崔晓娜
杨二龙
宋考平
王梅
董驰
张晓彤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeast Petroleum University
Original Assignee
Northeast Petroleum University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northeast Petroleum University filed Critical Northeast Petroleum University
Priority to CN201810655968.6A priority Critical patent/CN108829994B/zh
Publication of CN108829994A publication Critical patent/CN108829994A/zh
Application granted granted Critical
Publication of CN108829994B publication Critical patent/CN108829994B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/25Methods for stimulating production
    • E21B43/26Methods for stimulating production by forming crevices or fractures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mining & Mineral Resources (AREA)
  • Theoretical Computer Science (AREA)
  • Geology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Fluid Mechanics (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明涉及的是多层页岩拟三维压裂模型缝高求解方法,具体为:一、建立多层页岩地层压裂缝高物理模型,将模型中各层页岩地层所对应的层面顶底深度数据等输入计算机;二、对多层页岩地层压裂缝高物理模型进行网格划分;三、分别计算从射孔位置开始起裂,裂缝向上或向下生长到每个网格节点所对应的裂缝尖端应力强度因子为尖端网格节点所在页岩地层断裂韧性时的裂缝内流体净压力数值;四、给定任意一个射孔位置处进入裂缝的流体压力,求解净压力;五、用净压力到三得到的净压力矩阵中去对比寻找;步骤六、变化射孔位置处进入裂缝的流体压力,重复四五,实现不同压裂压力下裂缝缝高的计算。本发明可以准确计算和预测页岩气井压裂改造形成的裂缝缝高。

Description

多层页岩拟三维压裂模型缝高求解方法
技术领域
本发明涉及页岩储层压裂改造技术,具体涉及多层页岩拟三维压裂模型缝高求解方法。
背景技术
页岩气已经被视为替代煤和常规石油天然气的重要可持续能源。页岩气的全球总资源量约为4.56×1014m3,约占全球非常规油气资源的50%,相当于常规天然气的1.4倍,目前已被公认为是改变世界能源格局,使全球油气勘探展现良好发展前景的重要能源。
我国能源局于2016年9月14日发布的《页岩气发展规划(2016-2020)》表明,“十三五”将对页岩气开发及其配套开采技术发展给予高度重视和政策倾斜,2020年力争实现页岩气产量300×108m3。据统计2015年我国天然气供应存在大约600×108m3左右的缺口需要依靠进口,预计到2020年,国内天然气缺口将达1000×108m3。若能高效开发页岩气资源对我国国民经济发展和国家能源战略供给意义重大。
页岩气的有效开发需要通过水平井钻井结合大规模水力压裂技术,其中水力压裂裂缝的缝高是描述裂缝形态的一个重要参数。压裂设计时对裂缝高度的求解准确与否直接影响到压裂后的效果,甚至决定了压裂作业能否成功。
如求解的裂缝高度小于实际能够产生的裂缝高度,在压裂设计时,为了达到一定的设计缝高,需要加大泵注排量,由此可能产生压裂裂缝穿透隔层或遮挡地层延伸,造成许多不必要的生产事故,另外,由于额外增加了不必要的泵注排量,还会造成施工成本增加。
如求解的裂缝高度大于实际能够产生的裂缝高度,在压裂设计时,为了减小裂缝缝高,需要减小泵注排量,由此可能导致压裂后的裂缝缝高无法达到预期设计要求,使压裂改造体积减小,压裂效果较差甚至造成压裂失败。
所以,压裂设计时需要准确的计算压裂能够形成的裂缝缝高,才能保证压裂作业顺利开展。
以往研究人员对压裂裂缝高度的求解主要采用拟三维模型和全三维模型两种方法。由于全三维模型计算效率较拟三维模型差,并且在工业化应用,特别是大规模压裂大尺度裂缝模拟时,拟三维模型表现出更好的适用性,所以拟三维压裂模型是石油工业首选的高效裂缝求解模型。
学者们对拟三维模型中裂缝高度的求解,大致可以归为两类:
一类是简单的裂缝高度求解,假定裂缝一直在储层内部延伸,即只考虑地层含有一层储层的情况。这种高度求解模型偏于简单,但与实际地层状况相差甚远,也导致该模型求解的裂缝高度存在较大误差。
另一类是复杂的裂缝高度求解,假定裂缝可以穿透储层以外的上、下地层隔层或遮挡层延伸。这种高度求解模型较单层储层的假设更接近实际地层情况,但目前该类模型对于地层层数的假定都是有限的,而且一般只有三层,即包含储层和山、下遮挡层。
对于页岩地层,大量的研究和发表的论文已经表明,页岩的层理大量发育,这样页岩地层实际上是纵向多层地层的集合,压裂时,裂缝必然需要在多层地层中延伸形成一定的裂缝缝高。这也说明现有的压裂模型,无论是单层缝高模型,还是多层有限缝高模型都无法适用于实际页岩地层压裂的裂缝高度求解。
发明内容
本发明的目的是提供多层页岩拟三维压裂模型缝高求解方法,这种多层页岩拟三维压裂模型缝高求解方法用于解决现有技术中对页岩地层压裂裂缝高度的求解存在误差或不适应的问题。
本发明解决其技术问题所采用的技术方案是:这种多层页岩拟三维压裂模型缝高求解方法:
步骤一、建立多层页岩地层压裂缝高物理模型,将多层页岩地层压裂缝高物理模型中各层页岩地层所对应的层面顶底深度数据、地应力数据、岩石断裂韧性以及射孔深度数据输入计算机,做为待用数据;
步骤二、对多层页岩地层压裂缝高物理模型进行网格划分,以注入压裂液的射孔深度为起始位置,向下和向上分别划分500个网格,并对各个网格节点赋值,分别填充步骤一中输入对应的属性数据,向下是向井深增加方向,向上是向井深减小方向;
步骤三、根据断裂力学理论,分别计算从射孔位置开始起裂,裂缝向上或向下生长到每个网格节点所对应的裂缝尖端应力强度因子为尖端网格节点所在页岩地层断裂韧性时的裂缝内流体净压力数值P net(x);计算时采用积分方式,积分路径包括裂缝延伸所穿越的所有地层,并对各段地层进行分段求解,实现多层地层的叠加;
断裂力学中确定裂缝尖端应力强度因子的公式为:
Figure 907086DEST_PATH_IMAGE001
式中:C为裂缝半缝高,m;P net(x)为裂缝内任意位置的流体净压力,其数值等于流体压力与地层水平最小主应力的差值,MPa;
求解得到的流体净压力数值是一个m×n=500×500的矩阵,m代表从射孔位置向下的所有节点位置数据,n代表从射孔位置向上的所有节点位置数据;
步骤四、给定任意一个射孔位置处进入裂缝的流体压力,再将射孔位置处进入裂缝的流体压力与所在深度地层的水平最小主应力做差求解出净压力P net(x);
步骤五、用步骤四中求解得到的净压力P net(x)去步骤三的m×n=500×500的净压力矩阵中去对比寻找,首先增加m对应的网格值令其为1,然后不断增加n,找到第一个大于步骤四中求解得到的净压力P net(x)的位置即为此时的裂缝高度,裂缝下端点为射孔下侧1位置对应的深度,裂缝上端点为射孔上侧搜寻到位置所对应的深度;如m=1找不到满足的网格节点,增加m数值,直到其值等于500,结束寻找;寻找过程中,裂缝上端或下端有一侧到达m=500或n=500的位置即结束寻找,说明裂缝高度已经延伸到顶界面或底界面;
步骤六、变化射孔位置处进入裂缝的流体压力,重复步骤四-步骤五,实现不同压裂压力下裂缝缝高的计算,最后绘制裂缝缝高生长随射孔位置处进入裂缝的流体压力变化的关系曲线;
如给定射孔位置处进入裂缝的流体压力计算出的净压力P net(x)小于步骤三的 n=500×500的净压力矩阵中所有数据,说明射孔位置处进入裂缝的流体压力数值太小,裂缝不会发生起裂延伸。
本发明具有以下有益效果:
1、通过本发明提供的多层页岩拟三维压裂模型缝高求解方法,可以准确计算和预测页岩气井压裂改造形成的裂缝缝高,为压裂设计提供依据;
2、通过本发明提供的多层页岩拟三维压裂模型缝高求解方法,可以克服以往缝高计算模型无法考虑裂缝缝高在多层属性不同的页岩地层中延伸的缺点,使缝高计算结果更真实可靠;
3、通过本发明提供的多层页岩拟三维压裂模型缝高求解方法,可以避免由于预测裂缝高度小于实际裂缝高度而盲目加大泵注排量所产生的压裂裂缝穿透隔层或遮挡地层延伸现象,减少不必要的事故和降低施工成本。
4、通过本发明提供的多层页岩拟三维压裂模型缝高求解方法,可以避免由于预测裂缝高度大于实际裂缝高度而盲目减小泵注排量所产生的裂缝缝高无法达到预期设计要求现象,避免压裂参数设计不合理引起的压裂效果较差甚至是压裂失败。
附图说明
图1为多层页岩地层压裂缝高物理模型。
示意图只绘制了地层有6层情况,模型不局限于6层。C为裂缝半缝高,m;P net(x)为裂缝内任一点流体净压力,MPa;h i 为第i层地层的厚度,i=1,2,…, n,m;s i 为第i层地层的水平最小主应力,i=1,2, …, n,MPa;K ICi 为第i层地层岩石的断裂韧性,i=1,2, …, n,MPa·m1/2K I+ K I+ 分别表示裂缝下端和上端的应力强度因子,MPa·m1/2
具体实施方式
下面对本发明作进一步的说明:
这种多层页岩拟三维压裂模型缝高求解方法如下:
步骤一、建立多层页岩地层压裂缝高物理模型,将多层页岩地层压裂缝高物理模型中各层页岩地层所对应的层面顶底深度数据、地应力数据、岩石断裂韧性以及射孔深度数据输入计算机,做为待用数据;
步骤二、对多层页岩地层压裂缝高物理模型进行网格划分,以注入压裂液的射孔深度为起始位置,向下和向上分别划分500个网格,并对各个网格节点赋值,分别填充步骤一中输入对应的属性数据,向下是向井深增加方向,向上是向井深减小方向;
步骤三、根据断裂力学理论,分别计算从射孔位置开始起裂,裂缝向上或向下生长到每个网格节点所对应的裂缝尖端应力强度因子为尖端网格节点所在页岩地层断裂韧性时的裂缝内流体净压力数值P net(x);计算时采用积分方式,积分路径包括裂缝延伸所穿越的所有地层,并对各段地层进行分段求解,实现多层地层的叠加;
断裂力学中确定裂缝尖端应力强度因子的公式为:
Figure 721458DEST_PATH_IMAGE001
式中:C为裂缝半缝高,m;P net(x)为裂缝内任意位置的流体净压力,其数值等于流体压力与地层水平最小主应力的差值,MPa;
求解得到的流体净压力数值是一个m×n=500×500的矩阵,m代表从射孔位置向下的所有节点位置数据,n代表从射孔位置向上的所有节点位置数据;
步骤四、给定任意一个射孔位置处进入裂缝的流体压力,再将射孔位置处进入裂缝的流体压力与所在深度地层的水平最小主应力做差求解出净压力P net(x);
步骤五、用步骤四中求解得到的净压力P net(x)去步骤三的m×n=500×500的净压力矩阵中去对比寻找,首先增加m对应的网格值令其为1,然后不断增加n,找到第一个大于步骤四中求解得到的净压力P net(x)的位置即为此时的裂缝高度,裂缝下端点为射孔下侧1位置对应的深度,裂缝上端点为射孔上侧搜寻到位置所对应的深度;如m=1找不到满足的网格节点,增加m数值,直到其值等于500,结束寻找;寻找过程中,裂缝上端或下端有一侧到达m=500或n=500的位置即结束寻找,说明裂缝高度已经延伸到顶界面或底界面;
步骤六、变化射孔位置处进入裂缝的流体压力,重复步骤四-步骤五,实现不同压裂压力下裂缝缝高的计算,最后绘制裂缝缝高生长随射孔位置处进入裂缝的流体压力变化的关系曲线;
如给定射孔位置处进入裂缝的流体压力计算出的净压力P net(x)小于步骤三的 n=500×500的净压力矩阵中所有数据,说明射孔位置处进入裂缝的流体压力数值太小,裂缝不会发生起裂延伸。

Claims (1)

1.一种多层页岩拟三维压裂模型缝高求解方法,其特征在于包括如下步骤:
步骤一、建立多层页岩地层压裂缝高物理模型,将多层页岩地层压裂缝高物理模型中各层页岩地层所对应的层面顶底深度数据、地应力数据、岩石断裂韧性以及射孔深度数据输入计算机,做为待用数据;
步骤二、对多层页岩地层压裂缝高物理模型进行网格划分,以注入压裂液的射孔深度为起始位置,向下和向上分别划分500个网格,并对各个网格节点赋值,分别填充步骤一中输入对应的属性数据,向下是向井深增加方向,向上是向井深减小方向;
步骤三、根据断裂力学理论,分别计算从射孔位置开始起裂,裂缝向上或向下生长到每个网格节点所对应的裂缝尖端应力强度因子为尖端网格节点所在页岩地层断裂韧性时的裂缝内流体净压力数值P net(x);计算时采用积分方式,积分路径包括裂缝延伸所穿越的所有地层,并对各段地层进行分段求解,实现多层地层的叠加;
断裂力学中确定裂缝尖端应力强度因子的公式为:
Figure DEST_PATH_IMAGE001
式中:C为裂缝半缝高,m;P net(x)为裂缝内任意位置的流体净压力,其数值等于流体压力与地层水平最小主应力的差值,MPa;
求解得到的流体净压力数值是一个m×n=500×500的矩阵,m代表从射孔位置向下的所有节点位置数据,n代表从射孔位置向上的所有节点位置数据;
步骤四、给定任意一个射孔位置处进入裂缝的流体压力,再将射孔位置处进入裂缝的流体压力与所在深度地层的水平最小主应力做差求解出净压力P net(x);
步骤五、用步骤四中求解得到的净压力P net(x)去步骤三的m×n=500×500的净压力矩阵中去对比寻找,首先增加m对应的网格值令其为1,然后不断增加n,找到第一个大于步骤四中求解得到的净压力P net(x)的位置即为此时的裂缝高度,裂缝下端点为射孔下侧1位置对应的深度,裂缝上端点为射孔上侧搜寻到位置所对应的深度;如m=1找不到满足的网格节点,增加m数值,直到其值等于500,结束寻找;寻找过程中,裂缝上端或下端有一侧到达m=500或n=500的位置即结束寻找,说明裂缝高度已经延伸到顶界面或底界面;
步骤六、变化射孔位置处进入裂缝的流体压力,重复步骤四-步骤五,实现不同压裂压力下裂缝缝高的计算,最后绘制裂缝缝高生长随射孔位置处进入裂缝的流体压力变化的关系曲线;
如给定射孔位置处进入裂缝的流体压力计算出的净压力P net(x)小于步骤三的m×n=500×500的净压力矩阵中所有数据,说明射孔位置处进入裂缝的流体压力数值太小,裂缝不会发生起裂延伸。
CN201810655968.6A 2018-06-23 2018-06-23 多层页岩拟三维压裂模型缝高求解方法 Expired - Fee Related CN108829994B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810655968.6A CN108829994B (zh) 2018-06-23 2018-06-23 多层页岩拟三维压裂模型缝高求解方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810655968.6A CN108829994B (zh) 2018-06-23 2018-06-23 多层页岩拟三维压裂模型缝高求解方法

Publications (2)

Publication Number Publication Date
CN108829994A CN108829994A (zh) 2018-11-16
CN108829994B true CN108829994B (zh) 2022-08-23

Family

ID=64137504

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810655968.6A Expired - Fee Related CN108829994B (zh) 2018-06-23 2018-06-23 多层页岩拟三维压裂模型缝高求解方法

Country Status (1)

Country Link
CN (1) CN108829994B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114086946B (zh) * 2020-08-24 2023-08-22 中国石油天然气股份有限公司 压裂缝的缝高确定方法
CN113849999B (zh) * 2021-09-26 2022-08-23 西南石油大学 一种多层地层多缝压裂缝高预测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102518417A (zh) * 2011-12-16 2012-06-27 中国石油大学(北京) 一种页岩气藏水力压裂排量的确定方法
CN102852516A (zh) * 2012-04-19 2013-01-02 北京大学 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置
CN103670358A (zh) * 2013-11-25 2014-03-26 北京科技大学 水力压裂裂缝在砂泥岩薄互层地质界面断裂扩展判识方法
CN106844909A (zh) * 2017-01-05 2017-06-13 西南石油大学 一种致密储层压裂改造体积区的计算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102518417A (zh) * 2011-12-16 2012-06-27 中国石油大学(北京) 一种页岩气藏水力压裂排量的确定方法
CN102852516A (zh) * 2012-04-19 2013-01-02 北京大学 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置
CN103670358A (zh) * 2013-11-25 2014-03-26 北京科技大学 水力压裂裂缝在砂泥岩薄互层地质界面断裂扩展判识方法
CN106844909A (zh) * 2017-01-05 2017-06-13 西南石油大学 一种致密储层压裂改造体积区的计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"A Rigorous Hydraulic-Fracture Equilibrium-Height Model for Multilayer Formations";Songxia Liu等;《SPE Production & Operations》;20180531;第214-234页 *
"页岩储层压裂裂缝力学行为数值模拟研究";杨海;《中国优秀博硕士论文全文数据库 工程科技II辑》;20180515(第05期);第B019-10页 *

Also Published As

Publication number Publication date
CN108829994A (zh) 2018-11-16

Similar Documents

Publication Publication Date Title
CN111322050B (zh) 一种页岩水平井段内密切割暂堵压裂施工优化方法
US20160376885A1 (en) Method and Apparatus for Performance Prediction of Multi-Layered Oil Reservoirs
CN102930345B (zh) 一种基于梯度算法的自适应井网优化方法
CN100999992B (zh) 稠油蒸汽吞吐后期转蒸汽驱方法
CN102777157B (zh) 一种co2驱油气水异井注入油藏混驱开发方法
CN103899285B (zh) 多层砂岩油藏近极限含水期轮替水驱方法
CN104295292A (zh) 多层叠置煤层气***开采井设计方法
CN109359332A (zh) 一种浅薄层油藏数值模拟模型建立方法及转蒸汽驱的方法
CN109751037A (zh) 一种常压页岩气藏高频变排量体积压裂方法
CN108829994B (zh) 多层页岩拟三维压裂模型缝高求解方法
ZHANG et al. Main controlling factors and exploration direction of Permian to Triassic reservoir in the central sag of Junggar Basin
CN111827936B (zh) 一种批钻滚动式井群开采天然气水合物的***及方法
CN106545319A (zh) 断块油藏水驱后转人工气顶-边水双向驱提高采收率方法
Mlacnik et al. Using well windows in full-field reservoir simulation
US11649711B2 (en) Method for improving recovery ratio of braided well pattern of hugely thick or multi-layer oil and gas reservoir
CN111594132B (zh) 一种流体注入开发巨厚或多层油气藏的编织驱替井网
CN114357766A (zh) 一种长直井段井网整体体积压裂优化设计方法
CN110359892B (zh) 一种基于油藏条件的蒸汽驱方案优化设计方法
CN106600693A (zh) 含油饱和度的四维建模方法
CN114086946A (zh) 压裂缝的缝高确定方法
Mohammed et al. Combining Statistical Analysis with Simulation to Optimize Unconventional Completions-Upper and Lower Montney Formations, Canada
CN111594133B (zh) 基于多水平缝弓形井开发多层低渗油气藏的编织状井网
CN118148563A (zh) 一种超稠油采收方法
CN108979611A (zh) 一种油气藏储层水平缝水平井钻完井压裂改造方法
CN115455632B (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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220823