CN109723435B - 确定页岩气藏水平井重复压裂前现地应力场的方法及设备 - Google Patents
确定页岩气藏水平井重复压裂前现地应力场的方法及设备 Download PDFInfo
- Publication number
- CN109723435B CN109723435B CN201811517655.0A CN201811517655A CN109723435B CN 109723435 B CN109723435 B CN 109723435B CN 201811517655 A CN201811517655 A CN 201811517655A CN 109723435 B CN109723435 B CN 109723435B
- Authority
- CN
- China
- Prior art keywords
- stress field
- crack
- interference
- fracture
- horizontal well
- 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 47
- 238000011065 in-situ storage Methods 0.000 title claims abstract description 44
- 239000011148 porous material Substances 0.000 claims abstract description 49
- 230000003068 static effect Effects 0.000 claims abstract description 40
- 239000011435 rock Substances 0.000 claims abstract description 39
- 238000009826 distribution Methods 0.000 claims abstract description 28
- 238000004519 manufacturing process Methods 0.000 claims abstract description 25
- 239000007789 gas Substances 0.000 claims description 51
- 238000011160 research Methods 0.000 claims description 17
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 238000004891 communication Methods 0.000 claims description 8
- 230000009467 reduction Effects 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 7
- 239000003345 natural gas Substances 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- -1 hydrogen Chemical class 0.000 claims description 4
- 239000001257 hydrogen Substances 0.000 claims description 4
- 230000011218 segmentation Effects 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 2
- 230000008859 change Effects 0.000 abstract description 8
- 230000006698 induction Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000010276 construction Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000003079 shale oil Substances 0.000 description 2
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明实施例提供了一种确定页岩气藏水平井重复压裂前现地应力场的方法及设备。所述方法包括：获取水平井多段簇主裂缝干扰应力场，根据岩石力学动态参数及静态参数，确定原始地应力场，根据页岩气藏水平井生产规律，确定孔隙压力降干扰应力场；将上述三个应力场进行叠加，获取现地应力场分布。本发明实施例提供的确定页岩气藏水平井重复压裂前现地应力场的方法及设备，通过考虑缝内净压力分布影响下的多段簇主裂缝干扰应力场，明确单井产量随生产时间变化关系并确定孔隙压力降干扰应力场，可以获得不同段簇主裂缝分布特征和不同生产时间孔隙压力降对现地应力场的影响。
Description
技术领域
本发明实施例涉及非常规油气储层增产改造技术领域，尤其涉及一种确定页岩气藏水平井重复压裂前现地应力场的方法及设备。
背景技术
水平井分段压裂改造技术的快速发展使页岩气得到有效开发，其单井产量常表现为初期产量高、递减快、稳产期短的特征。而重复压裂技术是增加页岩气单井产量稳产期、减少递减速率、提高最终可采资源量的有效手段。以美国为例，自2000年以来约600口页岩油气水平井的重复压裂作业可显著提高页岩油气井的产量，有一半的井在重复压裂后的初始产量能达到初次压裂的80％以上，递减率更低。目前重复压裂的主要作用是张开老缝和起裂新缝，而重复压裂前地应力场分布规律是明确裂缝张开或起裂的先决条件，对暂堵重复压裂工艺参数的选取起到决定性作用，进一步达到增能挖潜效果。目前国内外研究学者获取页岩气藏水平井重复压裂前地应力未考虑页岩初次压裂多段簇裂缝干扰及缝内净压力分布对裂缝干扰的影响。因此，建立一种准确可靠的方法认识页岩气藏水平井重复压裂前地应力场分布情况，用来指导页岩气藏水平井重复压裂施工优化设计，就成为业界亟待解决的技术问题。
发明内容
针对现有技术存在的上述问题，本发明实施例提供了一种确定页岩气藏水平井重复压裂前现地应力场的方法及设备。
第一方面，本发明的实施例提供了一种确定页岩气藏水平井重复压裂前现地应力场的方法，包括：获取水平井多段簇主裂缝干扰应力场，根据岩石力学动态参数及静态参数，确定原始地应力场，根据页岩气藏水平井生产规律，确定孔隙压力降干扰应力场；将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加，获取现地应力场分布。
进一步地，所述获取水平井多段簇主裂缝干扰应力场，包括：采用每段裂缝的缝内净压力的加权平均值作为缝内净压力值，结合二维裂缝干扰应力场计算方法，获取水平井多段簇主裂缝干扰应力场。
进一步地，所述二维裂缝干扰应力场计算方法，包括：
σy诱导＝μs(σx诱导+σz诱导)
其中，σx诱导、σy诱导及σz诱导分别为x，y，z方向上的干扰应力；pn为裂缝壁面上净压力；r为裂缝中心距离二维空间中干扰应力场研究目标点的距离；r1为裂缝底部距离二维空间中干扰应力场研究目标点的距离；r2为裂缝顶部距离二维空间中干扰应力场研究目标点的距离；θ为二维空间中干扰应力场研究目标点偏离裂缝中心的角度；θ1为二维空间中干扰应力场研究目标点偏离裂缝底部的角度；θ2为二维空间中干扰应力场研究目标点偏离裂缝顶部的角度；c为缝高的一半；μs为静态泊松比。
进一步地，所述采用每段裂缝的缝内净压力的加权平均值作为缝内净压力值，包括：
Li＝i*Lf/n
其中，pn0为最大缝宽处裂缝壁面上净压力；Li为距离井筒的裂缝长度；γ为裂缝内净压力分布指数；G为岩石剪切模量；w为缝长对应的缝宽；Lf为等效缝长；n为缝长分割段数；μs为静态泊松比；pni为第i段裂缝的缝内净压力；pn为裂缝壁面上净压力；Es为静态杨氏模量。
进一步地，所述根据岩石力学动态参数及静态参数，确定原始地应力场，包括：
其中，σz为深度H处的垂向应力；ρr(h)为随深度变化的上覆岩体密度；H为压裂层位深度；g为重力加速度；σH和σh分别为最大、最小水平主应力；μs为岩石静态泊松比；KH和Kh分别为最大、最小水平地应力方向的构造应力系数；α为Biot系数；pp为孔隙压力。
进一步地，所述根据页岩气藏水平井生产规律，确定孔隙压力降干扰应力场，包括：
其中，Pe为原始地层压力；Pwf为井底流压；A和B为二项式产能方程系数；q为产量；qn为归一化拟产量；Gp为累积产气量；t为物质平衡时间；Zi为原始地层孔隙压力下的天然气偏差系数，采用Dranchuk-Abu-Kassem方法获取；Z为当前地层孔隙压力下的天然气偏差系数；G为原始地质储量；Δpp为孔隙压力降低值；P为当前底层压力。
进一步地，所述将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加，获取现地应力场分布，包括：
其中，σz为深度H处的垂向应力；g为重力加速度；σH和σh分别为最大、最小水平主应力；σ′H和σh′分别为现地应力场的最大、最小水平主应力；σx诱导(i)和σy诱导(i)分别为第i段裂缝在x，y方向上的干扰应力；μs为静态泊松比；α为Biot系数；Δpp为孔隙压力降低值；ρs为静态岩体密度。
第二方面，本发明的实施例提供了一种确定页岩气藏水平井重复压裂前现地应力场的装置，包括：
初始应力场获取模块，用于获取水平井多段簇主裂缝干扰应力场，根据岩石力学动态参数及静态参数，确定原始地应力场，根据页岩气藏水平井生产规律，确定孔隙压力降干扰应力场；
现地应力场获取模块，用于将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加，获取现地应力场分布。
第三方面，本发明的实施例提供了一种电子设备，包括：
至少一个处理器；以及
与处理器通信连接的至少一个存储器，其中：
存储器存储有可被处理器执行的程序指令，处理器调用程序指令能够执行第一方面的各种可能的实现方式中任一种可能的实现方式所提供的确定页岩气藏水平井重复压裂前现地应力场的方法。
第四方面，本发明的实施例提供了一种非暂态计算机可读存储介质，非暂态计算机可读存储介质存储计算机指令，计算机指令使计算机执行第一方面的各种可能的实现方式中任一种可能的实现方式所提供的确定页岩气藏水平井重复压裂前现地应力场的方法。
本发明实施例提供的确定页岩气藏水平井重复压裂前现地应力场的方法及设备，通过考虑缝内净压力分布影响下的多段簇主裂缝干扰应力场，明确单井产量随生产时间变化关系并确定孔隙压力降干扰应力场，可以获得不同段簇主裂缝分布特征和不同生产时间孔隙压力降对现地应力场的影响。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做一简单的介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的确定页岩气藏水平井重复压裂前现地应力场的方法流程图；
图2为现有技术提供的典型压裂施工曲线示意图；
图3为本发明实施例提供的页岩长水平段多段簇裂缝干扰应力场变化示意图；
图4为本发明实施例提供的页岩气藏水平井重复压裂前原地应力场与现地应力场对比示意图；
图5为本发明实施例提供的确定页岩气藏水平井重复压裂前现地应力场的装置结构示意图；
图6为本发明实施例提供的电子设备的实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。另外，本发明提供的各个实施例或单个实施例中的技术特征可以相互任意结合，以形成可行的技术方案，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时，应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围之内。
本发明实施例提供了一种确定页岩气藏水平井重复压裂前现地应力场的方法，参见图1，该方法包括：
101、获取水平井多段簇主裂缝干扰应力场，根据岩石力学动态参数及静态参数，确定原始地应力场，根据页岩气藏水平井生产规律，确定孔隙压力降干扰应力场；具体地，获取岩石力学动态参数及静态参数包括：
岩石力学动态参数通过测井数据(纵波时差、横波时差、泥质含量、密度、孔隙度)解释获取；进一步地，页岩气藏储层中发育的天然裂缝与层理导致测井资料响应产生异常(密度测井、声波测井、电阻率测井等等)，可间接反映对动态岩石力学参数的影响。其中动态岩石力学参数计算公式如下：
其中，μd为动态泊松比；Δts和Δtp分别为岩石的横波和纵波时差；Ed为动态杨氏模量；ρ为岩石密度。然后，建立岩石力学参数动静态转换关系：静态岩石力学参数通过室内实验同步测试获取，对同一研究目标井层(岩石)，其动静态参数(杨氏模量E、泊松比μ)之间存在满足μs＝A1+K1μd与Es＝A2+K2Ed的转换关系。利用线性回归方法对同一研究目标井层(岩石)的岩石力学参数进行回归，可得到动静态参数转换关系；其中，μs为静态泊松比；μd为动态泊松比；A1和K1分别为动静态泊松比转换关系常数；Es为静态杨氏模量；Ed为动态杨氏模量；A2和K2分别为动静态弹性模量转换关系常数。
102、将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加，获取现地应力场分布。
在上述实施例的基础上，本发明实施例中提供的确定页岩气藏水平井重复压裂前现地应力场的方法，所述获取水平井多段簇主裂缝干扰应力场，包括：采用每段裂缝的缝内净压力的加权平均值作为缝内净压力值，结合二维裂缝干扰应力场计算方法，获取水平井多段簇主裂缝干扰应力场。
在上述实施例的基础上，本发明实施例中提供的确定页岩气藏水平井重复压裂前现地应力场的方法，所述二维裂缝干扰应力场计算方法，包括：
σy诱导＝μs(σx诱导+σz诱导) (3)
其中，σx诱导、σy诱导及σz诱导分别为x，y，z方向上的干扰应力；pn为裂缝壁面上净压力；r为裂缝中心距离二维空间中干扰应力场研究目标点的距离；r1为裂缝底部距离二维空间中干扰应力场研究目标点的距离；r2为裂缝顶部距离二维空间中干扰应力场研究目标点的距离；θ为二维空间中干扰应力场研究目标点偏离裂缝中心的角度；θ1为二维空间中干扰应力场研究目标点偏离裂缝底部的角度；θ2为二维空间中干扰应力场研究目标点偏离裂缝顶部的角度；c为缝高的一半；μs为静态泊松比。页岩长水平段多段簇裂缝干扰应力场的变化请参见图3，图3中包括：干扰应力差介于0至2MPa阶段301、干扰应力差介于-2至-4MPa阶段302和干扰应力差介于-4至-6MPa阶段303。在图3中，横纵轴的单位为米，在裂缝缝长为0米，水平井筒方向3300米处，干扰应力差最大(超过6MPa)，其余位置压力差逐渐减小。
在上述实施例的基础上，本发明实施例中提供的确定页岩气藏水平井重复压裂前现地应力场的方法，所述采用每段裂缝的缝内净压力的加权平均值作为缝内净压力值，包括：
Li＝i*Lf/n (7)
其中，pn0为最大缝宽处裂缝壁面上净压力；Li为距离井筒的裂缝长度；γ为裂缝内净压力分布指数；G为岩石剪切模量；w为缝长对应的缝宽；Lf为等效缝长；n为缝长分割段数；μs为静态泊松比；pni为第i段裂缝的缝内净压力；pn为裂缝壁面上净压力；Es为静态杨氏模量。该实施例中的具体步骤如下：
将给定的等效缝长Lf分割成n等份，当i＝0,1,2,3,…,n时，缝长Li＝i*Lf/n；给定等效主裂缝最大缝宽值，由式(5-2)和(5-3)计算i＝0时的缝内净压力，此时缝长L0＝0；由式(5-1)计算i＝1,2,3,…,n时的缝内净压力pn,i；将计算的每一等份(i＝0,1,2,3,…,n时)所对应的缝内净压力pn,i进行加权平均得到缝内净压力pn，计算公式为式(6)。然后，将计算的缝内净压力pn代入(2)式，联立(3)式和(4)式求得多段簇裂缝应力干扰应力场。
在上述实施例的基础上，本发明实施例中提供的确定页岩气藏水平井重复压裂前现地应力场的方法，所述根据岩石力学动态参数及静态参数，确定原始地应力场，包括：
其中，σz为深度H处的垂向应力；ρr(h)为随深度变化的上覆岩体密度；H为压裂层位深度；g为重力加速度；σH和σh分别为最大、最小水平主应力；μs为岩石静态泊松比；KH和Kh分别为最大、最小水平地应力方向的构造应力系数；α为Biot系数；pp为孔隙压力。在此基础上请参见图2，在图2中，a曲线代表致密岩石，b曲线代表微缝高渗岩石，压力增大到F点时进入破裂阶段，此时加入前置液，之后进入加砂阶段，在该阶段加入携砂液，之后进入顶替阶段。在E点顶替阶段结束后就进入停泵阶段，在停泵阶段内裂缝闭合，裂缝闭合在C点处结束，在E点和C点之间包括管内摩阻和净裂缝延伸压力两个阶段。初次压裂施工曲线上裂缝闭合压力pc近似为最小水平主应力σh，进而由式(8-2)确定最小水平地应力方向的构造应力系数Kh。对于页岩水平井，破裂压力计算式表述为pf＝3σH-σz-αpp+σt，可根据初次压裂施工曲线破裂压力数值确定最大水平地应力方向的构造应力系数KH。
在上述实施例的基础上，本发明实施例中提供的确定页岩气藏水平井重复压裂前现地应力场的方法，所述根据页岩气藏水平井生产规律，确定孔隙压力降干扰应力场，包括：
其中，Pe为原始地层压力；Pwf为井底流压；A和B为二项式产能方程系数；q为产量；qn为归一化拟产量；Gp为累积产气量；t为物质平衡时间；Zi为原始地层孔隙压力下的天然气偏差系数，采用Dranchuk-Abu-Kassem方法获取；Z为当前地层孔隙压力下的天然气偏差系数；G为原始地质储量；Δpp为孔隙压力降低值；P为当前底层压力。在(9)式和(10)式之间还存在孔隙压力变化引起水平主应力的改变的情况。具体地，水平主应力的改变如下：
其中，ΔσH和Δσh分别为最大和最小水平主应力方向上的应力改变值；α为Biot系数；Δpp为地层孔隙压力降低值。
此外，天然气偏差系数Z采用Dranchuk-Abu-Kassem方法计算，具体包括：
采用牛顿迭代法对Z进行迭代计算，其中A1＝0.3265，A2＝-1.0700，A3＝-0.5339，A4＝0.01569，A5＝-0.05165，A6＝0.5457，A7＝-0.7361，A8＝0.1844，A9＝0.1056，A10＝0.6134，A11＝0.7210。其中，Tpr为对比温度；ρpr为对比密度。
在上述实施例的基础上，本发明实施例中提供的确定页岩气藏水平井重复压裂前现地应力场的方法，所述将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加，获取现地应力场分布，包括：
其中，σz为深度H处的垂向应力；g为重力加速度；σH和σh分别为最大、最小水平主应力；σ′H和σ′h分别为现地应力场的最大、最小水平主应力；σx诱导(i)和σy诱导(i)分别为第i段裂缝在x，y方向上的干扰应力；μs为静态泊松比；α为Biot系数；Δpp为孔隙压力降低值；ρs为静态岩体密度。
页岩气藏水平井重复压裂前原地应力场与现地应力场的对比效果可以参见图4，图4中包括：垂向应力曲线401、原最小水平主应力曲线402、原最大水平主应力曲线403、现最小水平主应力曲线404和现最大水平主应力曲线405。由图4中可见，深度(单位m)的范围是从2800m到4000m，垂向应力从60MPa变化到63MPa；原最小水平主应力曲线402从52MPa变化到60MPa；原最大水平主应力曲线403从52MPa变化到60MPa；现最小水平主应力曲线404从51MPa变化到60MPa；现最大水平主应力曲线405从51MPa变化到60MPa。由此可见，最大及最小水平主应力的范围稍有扩大，当前的最小水平主应力有较大波动。
本发明实施例提供的确定页岩气藏水平井重复压裂前现地应力场的方法，通过考虑缝内净压力分布影响下的多段簇主裂缝干扰应力场，明确单井产量随生产时间变化关系并确定孔隙压力降干扰应力场，可以获得不同段簇主裂缝分布特征和不同生产时间孔隙压力降对现地应力场的影响。
本发明各个实施例的实现基础是通过具有处理器功能的设备进行程序化的处理实现的。因此在工程实际中，可以将本发明各个实施例的技术方案及其功能封装成各种模块。基于这种现实情况，在上述各实施例的基础上，本发明的实施例提供了一种确定页岩气藏水平井重复压裂前现地应力场的装置，该装置用于执行上述方法实施例中的确定页岩气藏水平井重复压裂前现地应力场的方法。参见图5，该装置包括：
初始应力场获取模块501，用于获取水平井多段簇主裂缝干扰应力场，根据岩石力学动态参数及静态参数，确定原始地应力场，根据页岩气藏水平井生产规律，确定孔隙压力降干扰应力场；
现地应力场获取模块502，用于将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加，获取现地应力场分布。
本发明实施例提供的确定页岩气藏水平井重复压裂前现地应力场的装置，采用初始应力场获取模块和现地应力场获取模块，通过考虑缝内净压力分布影响下的多段簇主裂缝干扰应力场，明确单井产量随生产时间变化关系并确定孔隙压力降干扰应力场，可以获得不同段簇主裂缝分布特征和不同生产时间孔隙压力降对现地应力场的影响。
本发明实施例的方法是依托电子设备实现的，因此对相关的电子设备有必要做一下介绍。基于此目的，本发明的实施例提供了一种电子设备，如图6所示，该电子设备包括：至少一个处理器(processor)601、通信接口(Communications Interface)604、至少一个存储器(memory)602和通信总线603，其中，至少一个处理器601，通信接口604，至少一个存储器602通过通信总线603完成相互间的通信。至少一个处理器601可以调用至少一个存储器602中的逻辑指令，以执行如下方法：获取水平井多段簇主裂缝干扰应力场，根据岩石力学动态参数及静态参数，确定原始地应力场，根据页岩气藏水平井生产规律，确定孔隙压力降干扰应力场；将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加，获取现地应力场分布。
此外，上述的至少一个存储器602中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。例如包括：获取水平井多段簇主裂缝干扰应力场，根据岩石力学动态参数及静态参数，确定原始地应力场，根据页岩气藏水平井生产规律，确定孔隙压力降干扰应力场；将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加，获取现地应力场分布。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (6)
1.一种确定页岩气藏水平井重复压裂前现地应力场的方法，其特征在于，包括：
获取水平井多段簇主裂缝干扰应力场，根据岩石力学动态参数及静态参数，确定原始地应力场，根据页岩气藏水平井生产规律，确定孔隙压力降干扰应力场；
将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加，获取现地应力场分布；
所述获取水平井多段簇主裂缝干扰应力场，包括：
采用每段裂缝的缝内净压力的加权平均值作为缝内净压力值，结合二维裂缝干扰应力场计算方法，获取水平井多段簇主裂缝干扰应力场；
所述将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加，获取现地应力场分布，包括：
其中，
为深度H处的垂向应力；g为重力加速度；
和
分别为最大、最小水平主应力；
和
分别为现地应力场的最大、最小水平主应力；
和
分别为第i段裂缝在x，y方向上的干扰应力；
为静态泊松比；
为Biot系数；
为孔隙压力降低值；
为静态岩体密度；
所述二维裂缝干扰应力场计算方法，包括：
（2）
其中，
、
及
分别为x，y，z方向上的干扰应力；pn为裂缝壁面上净压力；r为裂缝中心距离二维空间中干扰应力场研究目标点的距离；r1为裂缝底部距离二维空间中干扰应力场研究目标点的距离；r2为裂缝顶部距离二维空间中干扰应力场研究目标点的距离；
为二维空间中干扰应力场研究目标点偏离裂缝中心的角度；
为二维空间中干扰应力场研究目标点偏离裂缝底部的角度；
为二维空间中干扰应力场研究目标点偏离裂缝顶部的角度；c为缝高的一半；
为静态泊松比；
所述采用每段裂缝的缝内净压力的加权平均值作为缝内净压力值，包括：
Li=i*Lf/n
其中，pn0为最大缝宽处裂缝壁面上净压力；Li为距离井筒的裂缝长度；
为裂缝内净压力分布指数；G为岩石剪切模量；
为缝长对应的缝宽；Lf为等效缝长；n为缝长分割段数；
为静态泊松比；pni为第i段裂缝的缝内净压力；pn为裂缝壁面上净压力；Es为静态杨氏模量；
其中，将给定的等效缝长Lf分割成n等份，当i=0,1,2,3,…,n时，缝长Li=i*Lf/n；给定等效主裂缝最大缝宽值，由式（5-2）和（5-3）计算i=0时的缝内净压力，此时缝长L0=0；由式（5-1）计算i=1,2,3,…,n时的缝内净压力pn,i；将计算的每一等份i=0,1,2,3,…,n时所对应的缝内净压力pn,i进行加权平均得到缝内净压力pn，计算公式为式（6），将计算的缝内净压力pn代入（2）式，联立（3）式和（4）式求得多段簇裂缝应力干扰应力场。
2.根据权利要求1所述的确定页岩气藏水平井重复压裂前现地应力场的方法，其特征在于，所述根据岩石力学动态参数及静态参数，确定原始地应力场，包括：
其中，
为深度H处的垂向应力；
为随深度变化的上覆岩体密度；H为压裂层位深度；g为重力加速度；
和
分别为最大、最小水平主应力；
为岩石静态泊松比；
和
分别为最大、最小水平地应力方向的构造应力系数；
为Biot系数；
为孔隙压力。
3.根据权利要求2所述的确定页岩气藏水平井重复压裂前现地应力场的方法，其特征在于，所述根据页岩气藏水平井生产规律，确定孔隙压力降干扰应力场，包括：
其中，Pe为原始地层压力；Pwf为井底流压；A和B为二项式产能方程系数；q为产量；qn为归一化拟产量；Gp为累积产气量；t为物质平衡时间；Zi为原始地层孔隙压力下的天然气偏差系数，采用Dranchuk-Abu-Kassem方法获取；Z为当前地层孔隙压力下的天然气偏差系数；G为原始地质储量；
为孔隙压力降低值；P为当前底层压力。
4.一种确定页岩气藏水平井重复压裂前现地应力场的装置，其特征在于，包括：
初始应力场获取模块，用于获取水平井多段簇主裂缝干扰应力场，根据岩石力学动态参数及静态参数，确定原始地应力场，根据页岩气藏水平井生产规律，确定孔隙压力降干扰应力场；
现地应力场获取模块，用于将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加，获取现地应力场分布；
所述获取水平井多段簇主裂缝干扰应力场，包括：
采用每段裂缝的缝内净压力的加权平均值作为缝内净压力值，结合二维裂缝干扰应力场计算方法，获取水平井多段簇主裂缝干扰应力场；
所述将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加，获取现地应力场分布，包括：
其中，
为深度H处的垂向应力；g为重力加速度；
和
分别为最大、最小水平主应力；
和
分别为现地应力场的最大、最小水平主应力；
和
分别为第i段裂缝在x，y方向上的干扰应力；
为静态泊松比；
为Biot系数；
为孔隙压力降低值；
为静态岩体密度；
所述二维裂缝干扰应力场计算方法，包括：
其中，
、
及
分别为x，y，z方向上的干扰应力；pn为裂缝壁面上净压力；r为裂缝中心距离二维空间中干扰应力场研究目标点的距离；r1为裂缝底部距离二维空间中干扰应力场研究目标点的距离；r2为裂缝顶部距离二维空间中干扰应力场研究目标点的距离；
为二维空间中干扰应力场研究目标点偏离裂缝中心的角度；
为二维空间中干扰应力场研究目标点偏离裂缝底部的角度；
为二维空间中干扰应力场研究目标点偏离裂缝顶部的角度；c为缝高的一半；
为静态泊松比；
所述采用每段裂缝的缝内净压力的加权平均值作为缝内净压力值，包括：
Li=i*Lf/n
其中，pn0为最大缝宽处裂缝壁面上净压力；Li为距离井筒的裂缝长度；
为裂缝内净压力分布指数；G为岩石剪切模量；
为缝长对应的缝宽；Lf为等效缝长；n为缝长分割段数；
为静态泊松比；pni为第i段裂缝的缝内净压力；pn为裂缝壁面上净压力；Es为静态杨氏模量；
其中，将给定的等效缝长Lf分割成n等份，当i=0,1,2,3,…,n时，缝长Li=i*Lf/n；给定等效主裂缝最大缝宽值，由式（5-2）和（5-3）计算i=0时的缝内净压力，此时缝长L0=0；由式（5-1）计算i=1,2,3,…,n时的缝内净压力pn,i；将计算的每一等份i=0,1,2,3,…,n时所对应的缝内净压力pn,i进行加权平均得到缝内净压力pn，计算公式为式（6），将计算的缝内净压力pn代入（2）式，联立（3）式和（4）式求得多段簇裂缝应力干扰应力场。
5.一种电子设备，其特征在于，包括：
至少一个处理器、至少一个存储器、通信接口和总线；其中，
所述处理器、存储器、通信接口通过所述总线完成相互间的通信；
所述存储器存储有能够被所述处理器执行的程序指令，所述处理器调用所述程序指令，以执行如权利要求1至4任一项所述的方法。
6.一种非暂态计算机可读存储介质，其特征在于，所述非暂态计算机可读存储介质存储计算机指令，所述计算机指令使所述计算机执行如权利要求1至4中任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811517655.0A CN109723435B (zh) | 2018-12-12 | 2018-12-12 | 确定页岩气藏水平井重复压裂前现地应力场的方法及设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811517655.0A CN109723435B (zh) | 2018-12-12 | 2018-12-12 | 确定页岩气藏水平井重复压裂前现地应力场的方法及设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109723435A CN109723435A (zh) | 2019-05-07 |
| CN109723435B true CN109723435B (zh) | 2022-07-22 |
Family
ID=66295234
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811517655.0A Active CN109723435B (zh) | 2018-12-12 | 2018-12-12 | 确定页岩气藏水平井重复压裂前现地应力场的方法及设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109723435B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113536706A (zh) * | 2021-07-15 | 2021-10-22 | 中国石油天然气股份有限公司 | 一种页岩油储层水平井分簇设计方法 |
| CN113280951B (zh) * | 2021-07-22 | 2021-10-08 | 中国科学院地质与地球物理研究所 | 一种建立峡谷区斜坡地应力场分布的方法 |
| CN114462272B (zh) * | 2022-01-20 | 2024-08-23 | 西南石油大学 | 一种深层复杂构造下页岩气水平井井眼轨迹优化方法 |
| CN116838308B (zh) * | 2023-08-11 | 2024-08-20 | 同济大学 | 一种重复压裂工艺优化方法及系统 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105201479A (zh) * | 2015-10-09 | 2015-12-30 | 西南石油大学 | 一种页岩储层水平井分段压裂射孔簇参数优化设计方法 |
| WO2018195909A1 (zh) * | 2017-04-26 | 2018-11-01 | 西南石油大学 | 一种低渗透油气藏直井体积压裂储层改造体积预测方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9618652B2 (en) * | 2011-11-04 | 2017-04-11 | Schlumberger Technology Corporation | Method of calibrating fracture geometry to microseismic events |
| CN105735960B (zh) * | 2016-03-22 | 2017-05-17 | 西南石油大学 | 一种低渗透油气藏水平井分段多簇压裂簇间距优化方法 |
| CN106593390B (zh) * | 2016-12-23 | 2019-09-20 | 西南石油大学 | 一种模拟水平井分段压裂裂缝扩展轨迹的方法 |
| CN106979000B (zh) * | 2017-05-06 | 2019-03-01 | 东北石油大学 | 水平井分段多簇压裂各射孔簇破裂压力计算方法 |
| CN108663414A (zh) * | 2018-03-22 | 2018-10-16 | 长江大学 | 一种确定压裂液饱和度范围的方法及系统 |
-
2018
- 2018-12-12 CN CN201811517655.0A patent/CN109723435B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105201479A (zh) * | 2015-10-09 | 2015-12-30 | 西南石油大学 | 一种页岩储层水平井分段压裂射孔簇参数优化设计方法 |
| WO2018195909A1 (zh) * | 2017-04-26 | 2018-11-01 | 西南石油大学 | 一种低渗透油气藏直井体积压裂储层改造体积预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109723435A (zh) | 2019-05-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109723435B (zh) | 确定页岩气藏水平井重复压裂前现地应力场的方法及设备 | |
| CN108868748B (zh) | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 | |
| CN106991236B (zh) | 一种基于四维地应力动态变化的重复压裂选井选层方法 | |
| CN105913155B (zh) | 考虑应力干扰及压裂液滤失的致密油产能预测方法及系统 | |
| CN113076676B (zh) | 非常规油气藏水平井压裂缝网扩展与生产动态耦合方法 | |
| CN113011048B (zh) | 一种致密砾岩油藏水平井重复压裂模拟方法 | |
| CN109522634A (zh) | 一种致密气多段体积压裂水平井数值分析方法 | |
| CN114427432B (zh) | 一种气藏剩余气开发潜力确定方法 | |
| CN109684685B (zh) | 一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法 | |
| CN106126866A (zh) | 基于地质力学模型的油气井防砂方法及装置 | |
| CN107145671B (zh) | 一种油藏数值模拟方法及系统 | |
| CN110929462B (zh) | 一种使用生产数据计算低渗气藏真实压力及气藏储量的方法 | |
| CN107965315A (zh) | 一种低渗透油藏产能预测方法 | |
| CN115324557A (zh) | 基于多因素分析预测压裂诱发套管变形风险程度的方法 | |
| CN109726450A (zh) | 确定页岩气藏水平井重复压裂时机的方法及设备 | |
| Sesetty et al. | Modeling and analysis of stimulation for fracture network generation | |
| CN115522918A (zh) | 深层砂岩油藏射孔井出砂压差剖面预测方法 | |
| CN115270533A (zh) | 重复压裂设计方法、装置、存储介质及电子设备 | |
| US10538995B2 (en) | Hydrocarbon production method comprising a productivity index of the wells subject to thermal effect | |
| CN110850470A (zh) | 一种利用地震资料计算地应力差异系数的方法 | |
| CN107916917B (zh) | 提高水力压裂双翼裂缝缝长的方法 | |
| CN106285657B (zh) | 用于确定缝洞型油藏溶洞发生坍塌事件的方法及系统 | |
| Semenov et al. | History matching of integrated reservoir simulation model for green field offshore Abu Dhabi | |
| CN111188613B (zh) | 一种致密气藏气井井控半径确定方法及系统 | |
| CN115324556A (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 |