CN110185427B - 一种缝内暂堵条件下天然裂缝开启时机的获取方法 - Google Patents
一种缝内暂堵条件下天然裂缝开启时机的获取方法 Download PDFInfo
- Publication number
- CN110185427B CN110185427B CN201910387470.0A CN201910387470A CN110185427B CN 110185427 B CN110185427 B CN 110185427B CN 201910387470 A CN201910387470 A CN 201910387470A CN 110185427 B CN110185427 B CN 110185427B
- Authority
- CN
- China
- Prior art keywords
- fracture
- crack
- hydraulic fracture
- natural
- hydraulic
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 239000012530 fluid Substances 0.000 claims abstract description 53
- 238000004364 calculation method Methods 0.000 claims abstract description 24
- 230000008569 process Effects 0.000 claims abstract description 19
- 230000000704 physical effect Effects 0.000 claims abstract description 4
- 239000011435 rock Substances 0.000 claims description 22
- 230000015572 biosynthetic process Effects 0.000 claims description 11
- 238000005086 pumping Methods 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 1
- 238000013461 design Methods 0.000 abstract description 4
- 239000003795 chemical substances by application Substances 0.000 description 6
- 230000000694 effects Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000002981 blocking agent Substances 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B33/00—Sealing or packing boreholes or wells
- E21B33/10—Sealing or packing boreholes or wells in the borehole
- E21B33/13—Methods or devices for cementing, for plugging holes, crevices, or the like
- E21B33/138—Plastering the borehole wall; Injecting into the formation
Abstract
本发明公开一种缝内暂堵条件下天然裂缝开启时机的获取方法，包括以下步骤：根据现场地质资料得到地层物性参数；再将水力裂缝划分为长度相等的N个单元体并依次编号，即每个单元体的长度为L/N；同时以缝内暂堵时刻作为初始时刻t0，将总计算时间t划分为间隔相同的m个时间节点，相邻时间节点的间隔时间为t/m；计算初始时刻水力裂缝内各个单元体的宽度；计算第k个时间节点水力裂缝内流体压力；再计算第k个时间节点水力裂缝上、下两侧天然裂缝入口处所受闭合应力；判断准则进行判断天然裂缝是否开启。本发明原理可靠，计算精度高，能够精确计算出缝内暂堵转向压裂过程中天然裂缝的开启时机，进而为压裂方案设计提供有效指导。
Description
技术领域
本发明涉及一种缝内暂堵条件下天然裂缝开启时机的获取方法，属于石油天然气勘探开发领域。
背景技术
水力压裂技术是低渗透油气藏增产改造的重要措施。水力压裂是利用地面高压泵组，以超过地层吸收能力的排量将压裂液泵入地层来产生水力裂缝，然后继续注入带有支撑剂(砂粒)的压裂液，使裂缝继续延伸并在其中充填支撑剂，当压裂液返排后，在地层压力作用下，支撑剂在裂缝中起到支撑裂缝的作用，阻止裂缝闭合，从而在地层中形成具有一定长度和流动能力的填砂裂缝。
缝内暂堵转向压裂是水力压裂的一种形式，具体是指在压裂过程中泵入暂堵剂对水力裂缝尖端进行暂时封堵，人为限制水力裂缝尖端向前延伸，迫使水力裂缝内部流体压力大幅上升，进而开启水力裂缝周围的天然裂缝，达到增加压裂改造范围的目的。因此，准确获取缝内暂堵条件下天然裂缝开启时机对于天然裂缝延伸过程预测和暂堵转向压裂工艺设计具有十分重要的意义。
暂堵失效是指压裂过程中水力裂缝内部压裂液突破尖端暂堵区域，导致暂堵段失去封堵作用，水力裂缝继续沿原路径向前延伸的现象；一般暂堵失效发生的条件为暂堵段两侧压力之差达到临界值，该临界值也被称为暂堵强度，由暂堵剂本身性质所决定。
天然裂缝是相对于人工裂缝而言的，指的是地层中由于地壳运动或其它自然因素而天然存在的一类裂缝。在水力压裂过程中，当水力裂缝向前延伸时通常会与天然裂缝相遇，此时存在两种可能情况：水力裂缝直接穿过天然裂缝沿原路径向前延伸，或水力裂缝沿天然裂缝所在路径向前延伸；缝内暂堵转向压裂主要适用于第一种情况，即当水力裂缝穿过天然裂缝后天然裂缝仍然保持闭合，随后通过泵入暂堵剂提升水力裂缝内部流体压力，进而迫使天然裂缝开启。另外，根据水力裂缝与天然裂缝相交时的相对位置，可将相交过程划分为两种类型：正交(垂直相交)和非正交。
诱导应力是指当材料某一位置受外力作用时，材料其它位置因抵抗该外力而诱发产生的作用力。针对水力压裂而言，在压裂过程中水力裂缝长度和宽度不断增加，其对周围岩石的挤压作用也不断增强，导致岩石内部所产生的诱导应力不断增大，该诱导应力也将间接影响天然裂缝的开启过程。
针对水力压裂过程中天然裂缝开启时机的获取方法，国内外学者已经做了许多相关研究，但大多数只针对于水力裂缝尖端与天然裂缝相遇时刻进行开启判断，没有具体分析当水力裂缝穿过天然裂缝之后，在尖端实施暂堵情况下天然裂缝的开启时机；同时，针对水力裂缝与天然裂缝非正交情况，学者们仅直接假设水力裂缝某一侧的天然裂缝将发生开启，未将两侧天然裂缝的受力情况进行对比，因此这些方法都不能很好地反映缝内暂堵条件下天然裂缝的实际开启过程。
发明内容
本发明主要是克服现有技术中的不足之处，提出一种缝内暂堵条件下天然裂缝开启时机的获取方法，该方法原理可靠，计算精度高，能够精确计算出缝内暂堵转向压裂过程中天然裂缝的开启时机，进而为压裂方案设计提供有效指导。
本发明解决上述技术问题所提供的技术方案是：一种缝内暂堵条件下天然裂缝开启时机的获取方法，包括以下步骤：
步骤S10、根据现场地质资料得到地层物性参数，并测得水力裂缝的缝长L；
步骤S20、再将水力裂缝划分为长度相等的N个单元体并依次编号，即每个单元体的长度为L/N；同时以缝内暂堵时刻作为初始时刻t0，将总计算时间t划分为间隔相同的m个时间节点，相邻时间节点的间隔时间为t/m；
步骤S30、计算初始时刻水力裂缝内各个单元体的宽度；
步骤S40、计算第k个时间节点水力裂缝内流体压力；
步骤S50、再计算第k个时间节点水力裂缝上、下两侧天然裂缝入口处所受闭合应力；
步骤S60、基于上述步骤S40和S50的计算结果，通过以下判断准则进行判断天然裂缝是否开启；
若开启，则时间节点k所对应的时刻t0+kt/m就是天然裂缝的开启时刻；
若未开启，则令k＝k+1，重复步骤S40-S50，直至天然裂缝开启或暂堵段失效；
所述判断准则为：
式中：Pk为第k个时间节点水力裂缝内流体压力，MPa；
为上侧天然裂缝入口处所受闭合应力，MPa；
为下侧天然裂缝入口处所受闭合应力，MPa；Pc为暂堵段封堵强度，MPa；Pr为地层流体压力，MPa。
进一步的技术方案是，所述步骤S30的计算公式为：
式中：p0为初始时刻水力裂缝内流体压力，MPa；σh为地层最小水平主应力，MPa；G为地层岩石剪切模量，MPa；υ为地层岩石泊松比，无因次；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；dij为裂缝单元i与裂缝单元j中点之间的距离，m；H为水力裂缝高度，m；α，β分别为经验系数，取α＝1，β＝2.3；i，j为水力裂缝单元体编号；Wi 0为第i个水力裂缝单元体在初始时刻的宽度，m。
进一步的技术方案是，所述步骤40的具体计算过程为：
步骤S401、根据下式计算第k个时间节点水力裂缝内预估流体压力；
式中：p0为初始时刻水力裂缝内流体压力，MPa；pk-1为第(k-1)时刻水力裂缝内的实际流体压力；
为第k个时间节点水力裂缝内预估流体压力，MPa；
步骤S402、根据上述计算得到的预估流体压力和下式计算第k个时间节点水力裂缝各单元体的预估宽度；
式中：
为第k个时间节点水力裂缝内预估流体压力，MPa；σh为地层最小水平主应力，MPa；G为地层岩石剪切模量，MPa；υ为地层岩石泊松比，无因次；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；dij为裂缝单元i与裂缝单元j中点之间的距离，m；H为水力裂缝高度，m；α，β分别为经验系数，取α＝1，β＝2.3；i，j为水力裂缝单元体编号；
为第k个时间节点水力裂缝各单元体的预估宽度，m；
步骤S403、通过下式计算预估宽度的误差α；
式中：
为第k个时间节点水力裂缝各单元的预估宽度，m；
为第k-1个时间节点水力裂缝各单元的宽度，m；H为水力裂缝高度，m；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；Q为缝内暂堵后压裂液的泵注排量，m3/s；Δt为相邻时间节点的间隔时间，s；i为水力裂缝单元体编号；α为误差；
步骤S404、设定求解精度ε，将上述得到的误差α与求解精度ε进行比较；
若α≤ε，则步骤S402和步骤S403中计算得到的
分别为第k个时间节点水力裂缝内流体压力和各个单元体的宽度；若α＞ε，则采用下式重新预估流体压力并重复步骤402-步骤404，直至满足α≤ε；
式中：ε为求解精度；
为第k个时间节点水力裂缝内预估流体压力，MPa；α为误差。
进一步的技术方案是，所述步骤S50的计算公式为：
式中：
为第k个时间节点水力裂缝上侧天然裂缝入口处所受闭合应力，MPa；
为第k个时间节点水力裂缝下侧天然裂缝入口处所受闭合应力，MPa；σH为地层最大水平主应力，MPa；σh为地层最小水平主应力，MPa；ω为水力裂缝与天然裂缝之间的夹角，度；dui为上侧天然裂缝入口单元中点与水力裂缝单元i中点之间的距离，m；dli为下侧天然裂缝入口单元中点与水力裂缝单元i中点之间的距离，m；H为水力裂缝高度，m；α，β为经验系数，取α＝1，β＝2.3；
为第k个时间节点水力裂缝单元体i的宽度，m；Cui，Cli分别为上侧和下侧天然裂缝入口单元体相对于水力裂缝单元体i的形状系数。
进一步的技术方案是，所述上侧、下侧天然裂缝入口单元体相对于水力裂缝单元体i的形状系数通过以下步骤得到；
(4)以第1个水力裂缝单元体中心点为原点，以水力裂缝长度方向为X轴，以过原点且垂直于水力裂缝壁面方向为Y轴，建立全局坐标系；
(5)上侧和下侧天然裂缝入口单元体中点在全局坐标系中的坐标可表示为：
式中：
为在全局坐标系中，上侧天然裂缝入口单元体中点坐标；
为在全局坐标系中，下侧天然裂缝入口单元体中点坐标；
为在全局坐标系中，水力裂缝与天然裂缝相交点的横坐标；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；ω为水力裂缝与天然裂缝之间的夹角，度；
(6)上侧和下侧天然裂缝入口单元体中点在以水力裂缝单元体i中点为基准的局部坐标系中的坐标可表示为：
式中：xui,yui为在局部坐标系中，上侧天然裂缝入口单元体中点坐标；xli,yli为在局部坐标系中，下侧天然裂缝入口单元体中点坐标；
为水力裂缝单元体i在全局坐标系中的坐标；
为在全局坐标系中，水力裂缝与天然裂缝相交点的横坐标；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；ω为水力裂缝与天然裂缝之间的夹角，度；
(4)再将步骤(3)中的公式带入到下式中求解得到上侧、下侧天然裂缝入口单元体相对于水力裂缝单元体i的形状系数；
Cij＝2G[-f1+yij(f2sin2γij-f3cos2γij)]
式中：
为第k个时间节点裂缝单元j所受正应力，MPa；G为地层岩石剪切模量，MPa；υ为地层岩石泊松比，无因次；dij为裂缝单元i与裂缝单元j中点之间的距离，m；H为水力裂缝高度，m；α，β为经验系数，取α＝1，β＝2.3；i，j为水力裂缝单元体编号；
为第i个水力裂缝单元体在初始时刻的宽度，m；Cij为裂缝单元j相对于裂缝单元i的形状系数；γij为裂缝单元i相对于裂缝单元j的偏转角；a为裂缝单元体半长，m；xij,yij为裂缝单元j中点在以裂缝单元i中点为基准的局部坐标系中的坐标值。
本发明具有以下优点：本发明原理可靠，计算精度高，能够精确计算出缝内暂堵转向压裂过程中天然裂缝的开启时机，进而为压裂方案设计提供有效指导。
附图说明
图1为本发明的流程框图。
具体实施方式
下面结合实施例和附图对本发明做更进一步的说明。
本发明的一种缝内暂堵条件下天然裂缝开启时机的获取方法，包括以下步骤：
步骤S10、根据现场地质资料得到地层物性参数，并测得水力裂缝的缝长L；
步骤S20、再将水力裂缝划分为长度相等的N个单元体并依次编号，即每个单元体的长度为L/N；同时以缝内暂堵时刻作为初始时刻t0，将总计算时间t划分为间隔相同的m个时间节点，相邻时间节点的间隔时间为t/m；
步骤S30、根据下式计算初始时刻水力裂缝内各个单元体的宽度；
式中：p0为初始时刻水力裂缝内流体压力，MPa；σh为地层最小水平主应力，MPa；G为地层岩石剪切模量，MPa；υ为地层岩石泊松比，无因次；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；dij为裂缝单元i与裂缝单元j中点之间的距离，m；H为水力裂缝高度，m；α，β分别为经验系数，取α＝1，β＝2.3；i，j为水力裂缝单元体编号；
为第i个水力裂缝单元体在初始时刻的宽度，m；
步骤S40、计算第k个时间节点水力裂缝内流体压力，其具体为：
步骤S401、根据下式计算第k个时间节点水力裂缝内预估流体压力；
式中：p0为初始时刻水力裂缝内流体压力，MPa；pk-1为第(k-1)时刻水力裂缝内的实际流体压力；
为第k个时间节点水力裂缝内预估流体压力，MPa；
步骤S402、根据上述计算得到的预估流体压力和下式计算第k个时间节点水力裂缝各单元体的预估宽度；
式中：
为第k个时间节点水力裂缝内预估流体压力，MPa；σh为地层最小水平主应力，MPa；G为地层岩石剪切模量，MPa；υ为地层岩石泊松比，无因次；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；dij为裂缝单元i与裂缝单元j中点之间的距离，m；H为水力裂缝高度，m；α，β分别为经验系数，取α＝1，β＝2.3；i，j为水力裂缝单元体编号；
为第k个时间节点水力裂缝各单元体的预估宽度，m；
步骤S403、通过下式计算预估宽度的误差α；
式中：
为第k个时间节点水力裂缝各单元的预估宽度，m；
为第k-1个时间节点水力裂缝各单元的宽度，m；H为水力裂缝高度，m；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；Q为缝内暂堵后压裂液的泵注排量，m3/s；Δt为相邻时间节点的间隔时间，即t/m，s；i为水力裂缝单元体编号；α为误差；
步骤S404、设定求解精度ε，将上述得到的误差α与求解精度ε进行比较；
求解精度一般取5％，这个求解精度主要是看求解过程中对于结果精度的要求，求解出的裂缝宽度越接近真实值，误差α就越小，如果没满足求解精度值就需要继续迭代；
若α≤ε，则步骤S402和步骤S403中计算得到的
分别为第k个时间节点水力裂缝内流体压力和各个单元体的宽度；若α＞ε，则采用下式重新预估流体压力并重复步骤402-步骤404，直至满足α≤ε；
式中：ε为求解精度；
为第k个时间节点水力裂缝内预估流体压力，MPa；α为误差；
步骤S50、再通过以下公式计算第k个时间节点水力裂缝上、下两侧天然裂缝入口处所受闭合应力；
式中：
为第k个时间节点水力裂缝上侧天然裂缝入口处所受闭合应力，MPa；
为第k个时间节点水力裂缝下侧天然裂缝入口处所受闭合应力，MPa；
σH为地层最大水平主应力，MPa；σh为地层最小水平主应力，MPa；ω为水力裂缝与天然裂缝之间的夹角，度；dui为上侧天然裂缝入口单元中点与水力裂缝单元i中点之间的距离，m；dli为下侧天然裂缝入口单元中点与水力裂缝单元i中点之间的距离，m；H为水力裂缝高度，m；α，β为经验系数，取α＝1，β＝2.3；
为第k个时间节点水力裂缝单元体i的宽度，m；Cui，Cli分别为上侧和下侧天然裂缝入口单元体相对于水力裂缝单元体i的形状系数；
上述上侧、下侧天然裂缝入口单元体相对于水力裂缝单元体i的形状系数通过以下步骤得到；
步骤S501、以第1个水力裂缝单元体中心点为原点，以水力裂缝长度方向为X轴，以过原点且垂直于水力裂缝壁面方向为Y轴，建立全局坐标系；
步骤S502、上侧和下侧天然裂缝入口单元体中点在全局坐标系中的坐标可表示为：
式中：
为在全局坐标系中，上侧天然裂缝入口单元体中点坐标；
为在全局坐标系中，下侧天然裂缝入口单元体中点坐标；
为在全局坐标系中，水力裂缝与天然裂缝相交点的横坐标；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；ω为水力裂缝与天然裂缝之间的夹角，度；
步骤S503、上侧和下侧天然裂缝入口单元体中点在以水力裂缝单元体i中点为基准的局部坐标系中的坐标可表示为：
式中：xui,yui为在局部坐标系中，上侧天然裂缝入口单元体中点坐标；xli,yli为在局部坐标系中，下侧天然裂缝入口单元体中点坐标；
为水力裂缝单元体i在全局坐标系中的坐标；
为在全局坐标系中，水力裂缝与天然裂缝相交点的横坐标；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；ω为水力裂缝与天然裂缝之间的夹角，度；
步骤S504、再将步骤(3)中的公式带入到下式中求解得到上侧、下侧天然裂缝入口单元体相对于水力裂缝单元体i的形状系数；
Cij＝2G[-f1+yij(f2sin2γij-f3cos2γij)]
式中：
为第k个时间节点裂缝单元j所受正应力，MPa；G为地层岩石剪切模量，MPa；υ为地层岩石泊松比，无因次；dij为裂缝单元i与裂缝单元j中点之间的距离，m；H为水力裂缝高度，m；α，β为经验系数，取α＝1，β＝2.3；i，j为水力裂缝单元体编号；
为第i个水力裂缝单元体在初始时刻的宽度，m；Cij为裂缝单元j相对于裂缝单元i的形状系数；γij为裂缝单元i相对于裂缝单元j的偏转角；a为裂缝单元体半长，即L/2N，m；xij,yij为裂缝单元j中点在以裂缝单元i中点为基准的局部坐标系中的坐标值；
步骤S60、基于上述步骤S40和S50的计算结果，通过以下判断准则进行判断天然裂缝是否开启；
若开启，则时间节点k所对应的时刻t0+kt/m就是天然裂缝的开启时刻；
若未开启，则令k＝k+1，重复步骤S40-S50，直至天然裂缝开启或暂堵段失效；
所述判断准则为：
式中：Pk为第k个时间节点水力裂缝内流体压力，MPa；
为上侧天然裂缝入口处所受闭合应力，MPa；
为下侧天然裂缝入口处所受闭合应力，MPa；Pc为暂堵段封堵强度，MPa；Pr为地层流体压力，MPa。
本实施例中初始时刻水力裂缝内各个单元体宽度计算公式的根据以下步骤得到：
1)初始时刻水力裂缝各单元体的宽度与其内部流体压力之间的关系式可表示为：
Cij＝2G[-f1+yij(f2sin2γij-f3cos2γij)] (2)
式中：
为第k个时间节点裂缝单元j所受正应力，MPa；G为地层岩石剪切模量，MPa；υ为地层岩石泊松比，无因次；dij为裂缝单元i与裂缝单元j中点之间的距离，m；H为水力裂缝高度，m；α，β为经验系数，通常取α＝1，β＝2.3；i，j为水力裂缝单元体编号；
为第i个水力裂缝单元体在初始时刻的宽度，m；Cij为裂缝单元j相对于裂缝单元i的形状系数；γij为裂缝单元i相对于裂缝单元j的偏转角；a为裂缝单元体半长(即L/2N)，m；xij,yij为裂缝单元j中点在以裂缝单元i中点为基准的局部坐标系中的坐标值；该局部坐标系以裂缝单元i中点为原点，以裂缝长度方向为X轴，以垂直于裂缝壁面方向为Y轴。
2)由于水力裂缝总是趋向于沿垂直于最小水平主应力方向延伸，裂缝单元j在外部将受最小水平主应力作用，在内部将受流体压力作用，因此其所受正应力可表示为：
式中：p0为初始时刻水力裂缝内流体压力，由暂堵剂实际泵注过程所决定，MPa；σh为地层最小水平主应力，MPa；
2)以第1个水力裂缝单元中心点为原点，以水力裂缝长度方向为X轴，以过原点且垂直于水力裂缝壁面方向为Y轴，建立全局二维直角坐标系；基于此坐标系，第i个水力裂缝单元中点的坐标可表示为：
式中：
为裂缝单元i在全局坐标系中的坐标值；i为裂缝单元的编号；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个。
3)根据局部坐标与全局坐标之间的变换规律，裂缝单元j中点在以裂缝单元i中点为基准的局部坐标系中的坐标可表示为：
式中：xij,yij为裂缝单元j中点在以裂缝单元i中点为基准的局部坐标系中的坐标值；i，j为水力裂缝单元体编号；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个。
4)将式(4)代入式(1)，式(6)代入式(2)、式(3)，即可得到初始时刻水力裂缝内各单元体的宽度
与其内部流体压力之间关系式：
本实施例中第k个时间节点水力裂缝上、下两侧天然裂缝入口处所受闭合应力的计算公式由以下步骤得到：
首先，水力裂缝上侧和下侧天然裂缝是指同一条天然裂缝的两翼，水力裂缝一般从天然裂缝中间部位穿过，从而将原本连续的天然裂缝一分为二，天然裂缝的两翼分别位于水力裂缝两侧(此处用上侧和下侧加以区分)。当水力裂缝与天然裂缝正交时，根据对称性，两侧天然裂缝将同时开启；当水力裂缝与天然裂缝非正交时，两侧天然裂缝具有开启先后顺序，因此在确定暂堵后天然裂缝开启时机时，还需要同步判断是哪一侧天然裂缝优先发生开启，此过程对于天然裂缝开启时机的确定非常重要。
类似地，为保证计算统一性和适应数值求解的需要，两侧天然裂缝也将被视为由长度为L/N的多个单元体组成，但计算过程只针对两侧天然裂缝入口处的第一个单元体进行。另外，水力裂缝周围一般存在着大量天然裂缝，此处以仅存在一条天然裂缝情况为例来说明求解过程，当存在多条天然裂缝时整体计算方法类似。
其次，天然裂缝所受闭合应力是指迫使天然裂缝保持闭合状态的作用力，可分为地层正应力和水力裂缝诱导应力两个部分，其中地层正应力部分可表示为：
式中：
——地层作用在天然裂缝壁面上的正应力，MPa；地层对于天然裂缝的作用力可分为正应力和切应力，仅正应力会迫使天然裂缝闭合。
σH——地层最大水平主应力，MPa；
σh——地层最小水平主应力，MPa；
ω——水力裂缝与天然裂缝之间的夹角，度。
水力裂缝诱导应力部分仍然采用式(7)表示，将式(7)与式(8)相叠加，即可得到下式：
式中：
为第k个时间节点水力裂缝上侧天然裂缝入口处所受闭合应力，MPa；
为第k个时间节点水力裂缝下侧天然裂缝入口处所受闭合应力，MPa；σH为地层最大水平主应力，MPa；σh为地层最小水平主应力，MPa；ω为水力裂缝与天然裂缝之间的夹角，度；dui为上侧天然裂缝入口单元中点与水力裂缝单元i中点之间的距离，m；dli为下侧天然裂缝入口单元中点与水力裂缝单元i中点之间的距离，m；H为水力裂缝高度，m；α，β为经验系数，取α＝1，β＝2.3；
为第k个时间节点水力裂缝单元体i的宽度，m；Cui，Cli分别为上侧和下侧天然裂缝入口单元体相对于水力裂缝单元体i的形状系数。
以上所述，并非对本发明作任何形式上的限制，虽然本发明已通过上述实施例揭示，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，当可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。
Claims (5)
1.一种缝内暂堵条件下天然裂缝开启时机的获取方法，其特征在于，包括以下步骤：
步骤S10、根据现场地质资料得到地层物性参数，并测得水力裂缝的缝长L；
步骤S20、再将水力裂缝划分为长度相等的N个单元体并依次编号，即每个单元体的长度为L/N；同时以缝内暂堵时刻作为初始时刻t0，将总计算时间t划分为间隔相同的m个时间节点，相邻时间节点的间隔时间为t/m；
步骤S30、计算初始时刻水力裂缝内各个单元体的宽度；
步骤S40、计算第k个时间节点水力裂缝内流体压力；
步骤S50、再计算第k个时间节点水力裂缝上、下两侧天然裂缝入口处所受闭合应力；
步骤S60、基于上述步骤S40和S50的计算结果，通过以下判断准则进行判断天然裂缝是否开启；
若开启，则时间节点k所对应的时刻t0+kt/m就是天然裂缝的开启时刻；
若未开启，则令k＝k+1，重复步骤S40-S50，直至天然裂缝开启或暂堵段失效；
所述判断准则为：
式中：Pk为第k个时间节点水力裂缝内流体压力，MPa；
为上侧天然裂缝入口处所受闭合应力，MPa；
为下侧天然裂缝入口处所受闭合应力，MPa；Pc为暂堵段封堵强度，MPa；Pr为地层流体压力，MPa。
2.根据权利要求1所述的一种缝内暂堵条件下天然裂缝开启时机的获取方法，其特征在于，所述步骤S30的计算公式为：
式中：p0为初始时刻水力裂缝内流体压力，MPa；σh为地层最小水平主应力，MPa；G为地层岩石剪切模量，MPa；υ为地层岩石泊松比，无因次；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；dij为裂缝单元i与裂缝单元j中点之间的距离，m；H为水力裂缝高度，m；α，β分别为经验系数，取α＝1，β＝2.3；i，j为水力裂缝单元体编号；Wi 0为第i个水力裂缝单元体在初始时刻的宽度，m。
3.根据权利要求1所述的一种缝内暂堵条件下天然裂缝开启时机的获取方法，其特征在于，所述步骤40的具体计算过程为：
步骤S401、根据下式计算第k个时间节点水力裂缝内预估流体压力；
式中：p0为初始时刻水力裂缝内流体压力，MPa；pk-1为第(k-1)时刻水力裂缝内的实际流体压力；
为第k个时间节点水力裂缝内预估流体压力，MPa；
步骤S402、根据上述计算得到的预估流体压力和下式计算第k个时间节点水力裂缝各单元体的预估宽度；
式中：
为第k个时间节点水力裂缝内预估流体压力，MPa；σh为地层最小水平主应力，MPa；G为地层岩石剪切模量，MPa；υ为地层岩石泊松比，无因次；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；dij为裂缝单元i与裂缝单元j中点之间的距离，m；H为水力裂缝高度，m；α，β分别为经验系数，取α＝1，β＝2.3；i，j为水力裂缝单元体编号；
为第k个时间节点水力裂缝各单元体的预估宽度，m；
步骤S403、通过下式计算预估宽度的误差α；
式中：
为第k个时间节点水力裂缝各单元的预估宽度，m；Wi k-1为第k-1个时间节点水力裂缝各单元的宽度，m；H为水力裂缝高度，m；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；Q为缝内暂堵后压裂液的泵注排量，m3/s；Δt为相邻时间节点的间隔时间，s；i为水力裂缝单元体编号；α为误差；
步骤S404、设定求解精度ε，将上述得到的误差α与求解精度ε进行比较；
若α≤ε，则步骤S402和步骤S403中计算得到的
分别为第k个时间节点水力裂缝内流体压力和各个单元体的宽度；若α＞ε，则采用下式重新预估流体压力并重复步骤402-步骤404，直至满足α≤ε；
式中：ε为求解精度；
为第k个时间节点水力裂缝内预估流体压力，MPa；α为误差。
4.根据权利要求1或2所述的一种缝内暂堵条件下天然裂缝开启时机的获取方法，其特征在于，所述步骤S50的计算公式为：
式中：
为第k个时间节点水力裂缝上侧天然裂缝入口处所受闭合应力，MPa；
为第k个时间节点水力裂缝下侧天然裂缝入口处所受闭合应力，MPa；σH为地层最大水平主应力，MPa；σh为地层最小水平主应力，MPa；ω为水力裂缝与天然裂缝之间的夹角，度；dui为上侧天然裂缝入口单元中点与水力裂缝单元i中点之间的距离，m；dli为下侧天然裂缝入口单元中点与水力裂缝单元i中点之间的距离，m；H为水力裂缝高度，m；α，β为经验系数，取α＝1，β＝2.3；Wi k为第k个时间节点水力裂缝单元体i的宽度，m；Cui，Cli分别为上侧和下侧天然裂缝入口单元体相对于水力裂缝单元体i的形状系数。
5.根据权利要求4所述的一种缝内暂堵条件下天然裂缝开启时机的获取方法，其特征在于，所述上侧、下侧天然裂缝入口单元体相对于水力裂缝单元体i的形状系数通过以下步骤得到；
(1)以第1个水力裂缝单元体中心点为原点，以水力裂缝长度方向为X轴，以过原点且垂直于水力裂缝壁面方向为Y轴，建立全局坐标系；
(2)上侧和下侧天然裂缝入口单元体中点在全局坐标系中的坐标表示为：
式中：
为在全局坐标系中，上侧天然裂缝入口单元体中点坐标；
为在全局坐标系中，下侧天然裂缝入口单元体中点坐标；
为在全局坐标系中，水力裂缝与天然裂缝相交点的横坐标；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；ω为水力裂缝与天然裂缝之间的夹角，度；
(3)上侧和下侧天然裂缝入口单元体中点在以水力裂缝单元体i中点为基准的局部坐标系中的坐标表示为：
式中：xui,yui为在局部坐标系中，上侧天然裂缝入口单元体中点坐标；xli,yli为在局部坐标系中，下侧天然裂缝入口单元体中点坐标；
为水力裂缝单元体i在全局坐标系中的坐标；
为在全局坐标系中，水力裂缝与天然裂缝相交点的横坐标；L为水力裂缝总长度，m；N为水力裂缝所划分的单元体数量，个；ω为水力裂缝与天然裂缝之间的夹角，度；
(4)再将步骤(3)中的公式带入到下式中求解得到上侧、下侧天然裂缝入口单元体相对于水力裂缝单元体i的形状系数；
Cij＝2G[-f1+yij(f2sin2γij-f3cos2γij)]
式中：
为第k个时间节点裂缝单元j所受正应力，MPa；G为地层岩石剪切模量，MPa；υ为地层岩石泊松比，无因次；dij为裂缝单元i与裂缝单元j中点之间的距离，m；H为水力裂缝高度，m；α，β为经验系数，取α＝1，β＝2.3；i，j为水力裂缝单元体编号；Wi 0为第i个水力裂缝单元体在初始时刻的宽度，m；Cij为裂缝单元j相对于裂缝单元i的形状系数；γij为裂缝单元i相对于裂缝单元j的偏转角；a为裂缝单元体半长，m；xij,yij为裂缝单元j中点在以裂缝单元i中点为基准的局部坐标系中的坐标值；该局部坐标系以裂缝单元i中点为原点，以裂缝长度方向为X轴，以垂直于裂缝壁面方向为Y轴。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910387470.0A CN110185427B (zh) | 2019-05-10 | 2019-05-10 | 一种缝内暂堵条件下天然裂缝开启时机的获取方法 |
| US16/550,336 US10844710B1 (en) | 2019-05-10 | 2019-08-26 | Method for acquiring opening timing of natural fracture under in-slit temporary plugging condition |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910387470.0A CN110185427B (zh) | 2019-05-10 | 2019-05-10 | 一种缝内暂堵条件下天然裂缝开启时机的获取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110185427A CN110185427A (zh) | 2019-08-30 |
| CN110185427B true CN110185427B (zh) | 2020-06-30 |
Family
ID=67714328
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910387470.0A Active CN110185427B (zh) | 2019-05-10 | 2019-05-10 | 一种缝内暂堵条件下天然裂缝开启时机的获取方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US10844710B1 (zh) |
| CN (1) | CN110185427B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111379537B (zh) * | 2020-03-20 | 2021-04-13 | 西南石油大学 | 一种定量评价暂堵剂封堵强度的方法 |
| CN112033812B (zh) * | 2020-08-12 | 2022-11-08 | 成都北方石油勘探开发技术有限公司 | 一种水力剪切压裂剪胀导流能力的测试方法及系统 |
| CN112065351B (zh) * | 2020-08-25 | 2021-09-17 | 中国石油大学(北京) | 水力裂缝内暂堵体信息的一体化确定方法、装置、设备 |
| CN114592823A (zh) * | 2020-12-04 | 2022-06-07 | 中国石油天然气股份有限公司 | 暂堵转向材料用量的确定方法及应用 |
| CN114592840B (zh) * | 2020-12-04 | 2023-10-27 | 中国石油天然气股份有限公司 | 暂堵压裂方法及其应用 |
| CN115387755B (zh) * | 2022-08-09 | 2023-06-30 | 中国石油大学(华东) | 一种co2地质封存期间沿断层泄漏的暂堵方法 |
| CN115680559B (zh) * | 2022-11-08 | 2023-08-08 | 西南石油大学 | 一种基于承压分散的裂缝性储层堵漏配方优选和设计方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106777663A (zh) * | 2016-12-12 | 2017-05-31 | 西南石油大学 | 一种考虑天然裂缝的压裂液滤失速度计算方法 |
| CN107313773A (zh) * | 2017-08-31 | 2017-11-03 | 中国石油化工股份有限公司 | 一种确定缝内暂堵转向所需封堵强度的方法 |
| CN107609258A (zh) * | 2017-09-07 | 2018-01-19 | 西南石油大学 | 一种页岩重复压裂转向裂缝起裂压力的计算方法 |
| CN108805365A (zh) * | 2018-07-05 | 2018-11-13 | 西南石油大学 | 一种非均质页岩储层裂缝诱导应力预测方法 |
| CN109505576A (zh) * | 2017-09-13 | 2019-03-22 | 中国石油化工股份有限公司 | 页岩水力压裂三维全耦合离散裂缝网络模拟方法及系统 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4850431A (en) * | 1988-05-06 | 1989-07-25 | Halliburton Company | Method of forming a plurality of spaced substantially parallel fractures from a deviated well bore |
| US7066284B2 (en) * | 2001-11-14 | 2006-06-27 | Halliburton Energy Services, Inc. | Method and apparatus for a monodiameter wellbore, monodiameter casing, monobore, and/or monowell |
| US8839875B2 (en) * | 2009-12-28 | 2014-09-23 | Ben M. Enis | Method and apparatus for sequestering CO2 gas and releasing natural gas from coal and gas shale formations |
| US10927604B2 (en) * | 2016-04-01 | 2021-02-23 | Board of Regents of the University of the Nevada System of Higher Education, on behalf of the University Nevada, Reno | Systems and methods for enhancing energy extraction from geothermal wells |
-
2019
- 2019-05-10 CN CN201910387470.0A patent/CN110185427B/zh active Active
- 2019-08-26 US US16/550,336 patent/US10844710B1/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106777663A (zh) * | 2016-12-12 | 2017-05-31 | 西南石油大学 | 一种考虑天然裂缝的压裂液滤失速度计算方法 |
| CN107313773A (zh) * | 2017-08-31 | 2017-11-03 | 中国石油化工股份有限公司 | 一种确定缝内暂堵转向所需封堵强度的方法 |
| CN107609258A (zh) * | 2017-09-07 | 2018-01-19 | 西南石油大学 | 一种页岩重复压裂转向裂缝起裂压力的计算方法 |
| CN109505576A (zh) * | 2017-09-13 | 2019-03-22 | 中国石油化工股份有限公司 | 页岩水力压裂三维全耦合离散裂缝网络模拟方法及系统 |
| CN108805365A (zh) * | 2018-07-05 | 2018-11-13 | 西南石油大学 | 一种非均质页岩储层裂缝诱导应力预测方法 |
Non-Patent Citations (2)
| Title |
|---|
| How Natural Fractures Could Affect Hydraulic-Fracture Geometry;Arash Dahi-Taleghanl等;《SPE Journal》;20140228;第161-171页 * |
| Numerical Modeling of Multistranded-Hydraulic-Fracture Propagation:Accounting for the Interaction Between Induced and Natural Fractures;Arash Dahi-Taleghani等;《SPE Journal》;20110930;第575-581页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| US20200355069A1 (en) | 2020-11-12 |
| US10844710B1 (en) | 2020-11-24 |
| CN110185427A (zh) | 2019-08-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110185427B (zh) | 一种缝内暂堵条件下天然裂缝开启时机的获取方法 | |
| Lian et al. | A study on casing deformation failure during multi-stage hydraulic fracturing for the stimulated reservoir volume of horizontal shale wells | |
| EP3271547B1 (en) | Stacked height growth fracture modeling | |
| Warpinski | Measurement of width and pressure in a propagating hydraulic fracture | |
| CN111236908A (zh) | 一种适用在低渗透致密气藏中的多段压裂水平井产能预测模型及产能敏感性分析的方法 | |
| CN113821953B (zh) | 基于动态地应力和剩余油特征的暂堵重复压裂设计方法 | |
| CN110219630A (zh) | 一种裂缝性砂岩气藏压裂水平井的压裂液返排量计算方法 | |
| CN109033504B (zh) | 一种油水井套管损坏预测方法 | |
| CN108090291B (zh) | 一种预测承压水悬挂式帷幕基坑涌水量的计算方法 | |
| CN104847379B (zh) | 一种奥灰顶部注浆加固工作面底板的方法 | |
| CN104653156B (zh) | 一种碎屑岩油气藏开发单元划分方法 | |
| CN107366530A (zh) | 一种深层页岩气藏增产方法及其应用 | |
| CN113177322B (zh) | 一种压裂单井控制储量计算方法 | |
| RU2666573C1 (ru) | Способ разработки нефтяной залежи с проведением повторного гидроразрыва пласта с изменением направления трещины | |
| CN109800523B (zh) | 一种水力压裂暂堵条件下非正交天然裂缝开启判定方法 | |
| CN113389534A (zh) | 一种水平井密切割压裂裂缝扩展预测与设计参数优化方法 | |
| CN106777663B (zh) | 一种考虑天然裂缝的压裂液滤失速度计算方法 | |
| CN111125905A (zh) | 耦合油藏流体流动的二维裂缝网络扩展模型及其模拟方法 | |
| CN115659736A (zh) | 一种深层页岩气水平井缝内转向压裂缝网扩展计算方法 | |
| CN114647957A (zh) | 一种基于酸压“改造区”理论的酸压效果评价方法 | |
| CN114592840B (zh) | 暂堵压裂方法及其应用 | |
| CN110714742B (zh) | 一种提高底水凝析气藏采收率的方法 | |
| Bazan et al. | Innovative analysis of treatment well data determines cluster efficiency, assesses fracture complexity, and provides fracture geometry in unconventional reservoirs | |
| US20230034507A1 (en) | Formation stimulation with acid etching model | |
| CN109339776B (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 | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20190830 Assignee: CHENGDU LEARN-PRACTICES TECHNOLOGY Co.,Ltd. Assignor: SOUTHWEST PETROLEUM University Contract record no.: X2023980040276 Denomination of invention: A Method for Obtaining the Opening Time of Natural Cracks under Temporary Sealing Conditions Granted publication date: 20200630 License type: Common License Record date: 20230825 |
|
| OL01 | Intention to license declared | ||
| OL01 | Intention to license declared |