CN106869896B - 一种低渗透油气藏直井体积压裂储层改造体积预测方法 - Google Patents
一种低渗透油气藏直井体积压裂储层改造体积预测方法 Download PDFInfo
- Publication number
- CN106869896B CN106869896B CN201710280797.9A CN201710280797A CN106869896B CN 106869896 B CN106869896 B CN 106869896B CN 201710280797 A CN201710280797 A CN 201710280797A CN 106869896 B CN106869896 B CN 106869896B
- Authority
- CN
- China
- Prior art keywords
- stress
- reservoir
- volume
- mpa
- fracture
- 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 41
- 239000011148 porous material Substances 0.000 claims abstract description 33
- 239000012530 fluid Substances 0.000 claims abstract description 26
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims description 22
- 230000035699 permeability Effects 0.000 claims description 14
- 230000004048 modification Effects 0.000 claims description 13
- 238000012986 modification Methods 0.000 claims description 13
- 230000009466 transformation Effects 0.000 claims description 13
- 238000006073 displacement reaction Methods 0.000 claims description 10
- 238000011065 in-situ storage Methods 0.000 claims description 10
- 238000010008 shearing Methods 0.000 claims description 7
- 239000004215 Carbon black (E152) Substances 0.000 claims description 6
- 238000007906 compression Methods 0.000 claims description 6
- 229930195733 hydrocarbon Natural products 0.000 claims description 6
- 150000002430 hydrocarbons Chemical class 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 239000011435 rock Substances 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 abstract description 3
- 230000002349 favourable effect Effects 0.000 abstract description 3
- 230000007812 deficiency Effects 0.000 abstract 1
- 230000007115 recruitment Effects 0.000 abstract 1
- 206010017076 Fracture Diseases 0.000 description 67
- 208000010392 Bone Fractures Diseases 0.000 description 60
- 238000012544 monitoring process Methods 0.000 description 14
- 230000000694 effects Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 239000009096 changqing Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- ZZUFCTLCJUWOSV-UHFFFAOYSA-N furosemide Chemical compound C1=C(Cl)C(S(=O)(=O)N)=CC(C(O)=O)=C1NCC1=CC=CO1 ZZUFCTLCJUWOSV-UHFFFAOYSA-N 0.000 description 1
- 230000001965 increasing effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK 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 OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种低渗透油气藏直井体积压裂储层改造体积预测方法,依次包括以下步骤:(1)计算水力裂缝在三维空间中产生的诱导应力;(2)计算压裂液滤失后的地层孔隙压力;(3)计算压裂液滤失后的地层孔隙弹性应力;(4)将上述三种应力场与原地应力场叠加获得新的地应力场,计算叠加后储层空间中的三向有效主应力的大小与方向;(5)计算储层空间中天然裂缝张开破裂判定系数M、天然裂缝剪切破裂区域判定系数S,从而预测体积压裂储层改造体积。本发明用于预测直井体积压裂后的储层改造体积,更具有操作性和准确性,为低渗透油气藏直井体积压裂压后效果评估与产量预测提供了有利的理论依据,克服了现有技术的不足。
Description
技术领域
本发明属于油气田开发领域,具体涉及一种低渗透油气藏储层直井体积压裂储层改造体积的预测方法。
背景技术
低渗透储层物性差,孔隙结构复杂、面孔率低、吼道细小,常规压裂技术很难达到预期的增产效果。体积压裂技术是改造低渗透油气田,实现增产的一项重要工艺措施。水平井体积压裂是低渗透油气藏增产改造的主要技术手段,但我国的低渗透油气储层以陆相沉积为主,部分油田区块中(如长庆油田、吉林油田的部分区块)发育的低渗透储层具有成藏面积小,纵向上小层多,厚度薄的特征,水平井体积改造技术表现出一定的不适应性,使得直井体积压裂技术被大规模的应用于此类储层的开发。准确的压后储层改造体积解释结果有利于合理评价直井体积压裂效果,优化压裂施工设计和精确预测压后产量。目前国内外常用的体积压裂裂缝现场监测手段主要有微地震监测、测斜仪监测和分布式声传感裂缝监测。其中,微地震监测能够检测到复杂裂缝网络体的方位、缝长、缝宽、缝高和倾角,是目前应用较多的裂缝监测手段,但该方法技术成本高,不适合大规模多井次应用。若采用储层改造体积数学模型预测体积压裂储层改造体积,则可大大降低成本。现有的储层改造体积预测方法主要包括半解析法、产量拟合法和离散裂缝网络模拟法。半解析法和产量拟合法实施时对基础数据要求较高,前者需要采用目标区块部分井的微地震监测数据进行校准,才能得出可靠的预测结果,而后者则必须根据压后产量数据预测储层改造体积大小。离散裂缝网络模拟法虽不需要事先获取微地震监测数据和压后产量数据,但该方法的基础理论较为复杂,计算量大,模型收敛性差,不易实施。由此可知,现有的储层改造体积预测方法均存在一定程度的不足,不适合大量应用于现场实践。
综上所述,目前需要的低渗透油气藏直井体积压裂储层改造体积预测方法应具有以下两个特点:(1)实施时对基础数据要求较低,不需要事先获得大量的基础数据如微地震监测数据和压后产量数据等,就可直接预测体积压裂储层改造体积;(2)基础理论完善,思路简洁直观,可操作性强,且能获得准确的预测结果。
发明内容
本发明的目的在于提供一种低渗透油气藏直井体积压裂储层改造体积预测方法,用于预测直井体积压裂后的储层改造体积,更具有操作性和准确性,为低渗透油气藏直井体积压裂压后效果评估与产量预测提供了有利的理论依据,克服了现有技术的不足。
为达到以上技术目的,本发明提供了以下技术方案。
首先,建立三维水力裂缝干扰模式下的地应力场计算模型,先分别计算出储层空间中任意一点处的水力裂缝诱导应力、压裂液滤失后的地层孔隙压力和孔隙弹性应力,再基于弹性力学基础理论将上述三种应力场与原始地应力场进行叠加获得三维水力裂缝干扰模式下的地应力场;其次,计算三维水力裂缝干扰模式下储层空间中的三向有效主应力的大小与方向;最后,在存在三维水力裂缝干扰的复杂地应力场条件下,计算分析储层空间中的天然裂缝是否发生张开或剪切破裂,将这两种破裂方式所波及到的总体积等效为储层改造体积。
一种低渗透油气藏直井体积压裂储层改造体积预测方法,依次包括以下步骤:
(1)计算水力裂缝在三维空间中产生的诱导应力。基于三维位移不连续理论,建立三维水力裂缝诱导应力场计算模型(Wu K.Numerical modeling of complex hydraulicfracture development in unconventional reservoirs[D].The University of Texasat Austin,2014)。计算三维空间中任意点i处的水力裂缝诱导应力大小,计算公式如下:
式中:
σxx i、σyy i、σzz i、σxy i、σxz i、σyz i分别为三维空间中任意点i处水力裂缝产生的正向诱导应力与剪切诱导应力,MPa;
均为边界影响系数;
N为水力裂缝被划分的总单元个数;
DsL j、DsH j、Dn j分别为水力裂缝单元j上的切向与法向位移不连续量,可由式(2)中给出的条件,结合式(1)反求得到;
σzz j为水力裂缝单元j上受到的沿原始最小地应力方向(平行于z轴)的正应力边界值,一般为水力裂缝壁面受到的净压力,MPa;
σxz j、σyz j均为水力裂缝单元j上受到的剪切应力边界值,一般为0,MPa;
pnet j为水力裂缝单元j上受到的净压力,MPa。
(2)计算压裂液滤失后的地层孔隙压力。低渗透储层基质渗透率极低,可忽略压裂液向基质中的滤失。一般情况下,储层中发育的天然裂缝是压裂液滤失的主要通道,在低渗透储层中仅考虑压裂液沿天然裂缝的滤失行为。由于地层原始流体性质和储层岩石结构特征存在差异,在低渗透油藏与低渗透气藏中分别采用不同的方法计算压裂液滤失导致的地层孔隙压力变化值。具体计算公式如下(Warpinski N R,Wolhart S L,Wright CA.Analysis and Prediction of Microseismicity Induced by Hydraulic Fracturing[J].SPE Journal,2004,9(01):24-33):
低渗透油藏:
低渗透气藏:
式中:
Pi为压裂液滤失后地层中任意点i处的孔隙压力,MPa;
P0为地层原始孔隙压力,MPa;
Pf为水力裂缝缝内液体压力,MPa;
λ、μ为任意滤失点i在三维水力裂缝椭球体坐标系中的坐标值,m2;
b为水力裂缝在井筒处的半缝高,m;
φ为天然裂缝孔隙度,无因次;
μl为压裂液粘度,mPa·s;
c为天然裂缝压缩系数,MPa-1;
k为储层渗透率,D;
t为滤失时间,s;
L为水力裂缝长度,m;
xi为任意滤失点i在三维直角坐标系中沿缝长方向(x轴方向)的坐标值,m;
yi为任意滤失点i在三维直角坐标系中沿缝高方向(y轴方向)的坐标值,m。
(3)计算压裂液滤失后的地层孔隙弹性应力。压裂液的滤失会增加局部储层的孔隙压力,这种作用将会扰动水力裂缝周围的应力。在低渗透气藏储层中,可忽略孔隙弹性应力影响,但在油藏储层中,此孔隙弹性应力不能被忽略,其计算公式如下(Smith MB.Stimulation design for short,precise hydraulic fractures[J].Society ofPetroleum Engineers Journal,1985,25(03):371-379):
式中:
Δσi为任意点i处地层孔隙压力变化产生的孔隙弹性应力,MPa;
A为孔弹性常数,无因次;
ν地层岩石泊松比,无因次;
为描述流体侵入的尺寸和形状的参数,无因次;
α为Biot系数,无因次;
h为天然裂缝半缝高,m。
(4)将上述三种应力场与原地应力场叠加获得新的地应力场,计算叠加后储层空间中的三向有效主应力的大小与方向。首先,基于弹性力学基础理论将上述三种诱导应力场与原地应力场进行叠加,其计算公式如下:
其次,将方程(6)代入方程(7)并求解方程(7)得到各应力场叠加后储层空间中的三向有效主应力的大小。
然后,计算各应力场叠加后储层空间中的三向有效主应力的方向,其计算公式如下:
式中:
σxx为水力裂缝产生的沿原始最大水平地应力方向(平行于x轴方向)的诱导应力,具体各点处的应力值为σxx i,MPa;
σyy为水力裂缝产生的沿原始垂向地应力方向(平行于y轴方向)的诱导应力,具体各点处的应力值为σyy i,MPa;
σzz为水力裂缝产生的沿原始最小水平地应力方向(平行于z轴方向)的诱导应力,具体各点处的应力值为σzz i,MPa;
σxy、σxz、σyz均为水力裂缝产生的剪切诱导应力,具体各点处的应力值为σxy i、σxz i、σyz i,MPa;
P为压裂液滤失后的地层中的孔隙压力,具体各点处的应力值为Pi,MPa;
Δσ为压裂液滤失后的地层中的孔隙弹性应力,具体各点处的应力值为Δσi,MPa;
σh为原始最小水平地应力,MPa;
σH为原始最大水平地应力,MPa;
σV为原始垂向地应力,MPa;
σi(i=1,2,3)为以上三种应力场与原地应力场叠加以后储层空间中的三向有效主应力的大小,其中σ1>σ2>σ3,MPa;
βxi、βyi、βzi(i=1,2,3)分别为以上三种应力场与原地应力场叠加以后储层空间中的各向有效主应力σi(i=1,2,3)分别与x、y、z轴方向的夹角,0°~90°。
本发明中的所有计算以压应力为正,张应力为负。
(5)计算储层空间中天然裂缝张开破裂判定系数M、天然裂缝剪切破裂区域判定系数S,从而预测体积压裂储层改造体积,确定其长度、宽度和高度,获得储层改造体积大小。天然裂缝发育是低渗透油气藏获得理想体积压裂改造效果的必要条件。天然裂缝是储层力学上的薄弱环节,储层压裂改造过程中天然裂缝更易先于基岩发生张开或剪切破裂,从而在三维储层空间中形成复杂裂缝网络体,因此可将这两种破裂方式的总波及体积等效为储层改造体积。由于储层改造体积形态较为复杂,为方便计算储层改造体积大小,按照现阶段的常规计算思路将其等效成长方体。按照长方体计算体积大小的方式(“长×宽×高”)计算储层改造体积大小。
M为天然裂缝张开破裂判定系数,M>0则代表在水力压裂过程中该处的天然裂缝会发生张开破裂,其计算公式为:
S为天然裂缝剪切破裂区域判定系数,S>0则代表在水力压裂过程中该处的天然裂缝会发生剪切破裂,其计算公式为:
式中:
σn为天然裂缝壁面受到的正应力,MPa;
为天然裂缝面法向矢量与应力场叠加后的最大有效主应力σ1的夹角,0°~90°;
θ为天然裂缝面法向矢量与应力场叠加后的最小有效主应力σ3的夹角,0°~90°;
α为天然裂缝倾角即天然裂缝面法向矢量与y轴的夹角,0°~90°;
为天然裂缝走向与原始最大水平地应力(x轴方向)的夹角,0°~90°;
τ为天然裂缝面受到的剪切应力,MPa;
kf为天然裂缝壁面摩擦系数,无因次;
σo为天然裂缝内聚力,MPa。
储层空间中S>0以及M>0的总波及体积即为直井体积压裂储层改造体积,通过观察该储层改造体积的俯视图和侧视图确定其长度、宽度和高度,按照“长×宽×高”方式计算获得储层改造体积大小。
本发明中涉及的计算公式和参数较多,若一一列出会显得过于累赘,为保证本发明的简洁直观,因此仅列出了主要的计算公式和参数,对于未给出的计算公式和参数则列出了相应的参考文献。
与现有技术相比,本发明的有益效果在于:(1)基于弹性力学基础理论和三维位移不连续法建立了考虑三维水力裂缝干扰模式下的复杂地应力场数学模型,计算分析了天然裂缝在复杂地应力场条件下的张开和剪切破裂行为,并基于体积压裂储层改造体积力学形成机理,将以上两种破裂方式的总波及体积等效为储层改造体积,给出了一种更加准确直观的直井体积压裂储层改造体积预测方法。(2)本方法简洁直观,可实施性强,在预测储层改造体积时不需要事先获得微地震监测数据或压后产量数据,弥补了现有的计算方法对基础数据要求较高的不足。
附图说明
图1为水力裂缝三维位移不连续法单元划分示意图。
图2为Q30井压裂后本方法预测结果俯视图与微地震事件监测结果俯视图。
图3为Q30井压裂后本方法预测结果侧视图与微地震事件监测结果侧视图。
具体实施方式
以下结合附图及现场运用实例,对本发明进一步详细说明。
以东部某油田的一口致密砂岩油藏直井(Q30)为例,该井完钻井深2550m,储层深度为2442.4~2446.2m,厚度为3.8m,油层平均孔隙度为9.2%,平均渗透率为0.62mD,属于低孔、低渗储层。本井采用大排量、低砂比、大液量低黏液体体系(以滑溜水为主)的技术做法进行体积压裂,最大限度使储层改造体积最大化,提高储层动用程度,从而提高单井产量。储层内天然裂缝发育,裂缝倾角较大约为80°,与最大水平地应力方向的夹角约为15°。其它基本参数见下表1所示。
表1 Q30井基本参数表
步骤1,采用表1中的数据,运用公式(1)、(2)计算出水力裂缝在三维地层空间中产生的诱导应力。图1展示了水力裂缝按照三维位移不连续法进行单元划分的具体方式以及地层中任意点i处水力裂缝产生的正向诱导应力和切向诱导应力。
步骤2,采用表1中的数据,运用公式(3)计算压裂液滤失后的地层孔隙压力(若是低渗透气藏则使用公式(4)进行计算)。
步骤3,采用表1中的数据,运用公式(5)计算压裂液滤失后的地层孔隙弹性应力(若是低渗透气藏,压裂液滤失产生的地层孔隙弹性应力可忽略)
步骤4,采用表1中的数据,结合步骤1~3的计算结果,运用公式(6)~(8)计算以上三种应力场与原地应力叠加后的新地应力场的空间三向有效主应力的大小和方向。
步骤5,采用表1中的数据,结合步骤4中的计算结果,运用公式(9)~(10)预测油井Q30体积压裂后的储层改造体积,通过观察该储层改造体积预测结果的俯视图和侧视图确定其长度、宽度和高度,按照“长×宽×高”方式计算获得储层改造体积大小。
观察图2和图3可知,采用以上步骤进行计算,预测Q30井体积压裂后的储层改造体积大小约为154×104m3(长300m,宽78m,高66m)。Q30井压裂施工结束后,微地震监测解释结果显示储层中形成了一个改造体积约为141×104m3(长300m,宽76m,高62m)的复杂裂缝网络体,与本发明预测结果十分接近。同时,Q30井压裂后微地震事件监测结果俯视图和侧视图与本方法预测结果的俯视图和侧视图基本重合。说明本发明提出的低渗透油气藏直井体积压裂储层改造体积预测方法较为合理,可为储层压裂施工参数优化设计,压后效果评价和压后产量精确预测提供有利的参考。
Claims (6)
1.一种低渗透油气藏直井体积压裂储层改造体积预测方法,依次包括以下步骤:
(1)计算水力裂缝在三维空间中产生的诱导应力;
(2)计算压裂液滤失后的地层孔隙压力;
(3)计算压裂液滤失后的地层孔隙弹性应力;
(4)将上述三种应力场与原地应力场叠加获得新的地应力场,计算叠加后储层空间中的三向有效主应力的大小与方向;
(5)计算储层空间中天然裂缝张开破裂判定系数M、天然裂缝剪切破裂区域判定系数S,从而预测体积压裂储层改造体积,确定其长度、宽度和高度,获得储层改造体积大小。
2.如权利要求1所述的一种低渗透油气藏直井体积压裂储层改造体积预测方法,其特征在于,所述步骤(1)计算水力裂缝在三维空间中产生的诱导应力,是指基于三维位移不连续理论,建立三维水力裂缝诱导应力场计算模型,计算三维空间中任意点i处的水力裂缝诱导应力大小,计算公式如下:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<mi>i</mi>
</msup>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
<mo>,</mo>
<mi>s</mi>
<mi>L</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mrow>
<mi>s</mi>
<mi>L</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
<mo>,</mo>
<mi>s</mi>
<mi>H</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mrow>
<mi>s</mi>
<mi>H</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
<mo>,</mo>
<mi>n</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mi>n</mi>
</msub>
<mi>j</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
<mi>i</mi>
</msup>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
<mo>,</mo>
<mi>s</mi>
<mi>L</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mrow>
<mi>s</mi>
<mi>L</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
<mo>,</mo>
<mi>s</mi>
<mi>H</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mrow>
<mi>s</mi>
<mi>H</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
<mo>,</mo>
<mi>n</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mi>n</mi>
</msub>
<mi>j</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>z</mi>
<mi>z</mi>
</mrow>
</msub>
<mi>i</mi>
</msup>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mi>z</mi>
<mo>,</mo>
<mi>s</mi>
<mi>L</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mrow>
<mi>s</mi>
<mi>L</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mi>z</mi>
<mo>,</mo>
<mi>s</mi>
<mi>H</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mrow>
<mi>s</mi>
<mi>H</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mi>z</mi>
<mo>,</mo>
<mi>n</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mi>n</mi>
</msub>
<mi>j</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
<mi>i</mi>
</msup>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
<mo>,</mo>
<mi>s</mi>
<mi>L</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mrow>
<mi>s</mi>
<mi>L</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
<mo>,</mo>
<mi>s</mi>
<mi>H</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mrow>
<mi>s</mi>
<mi>H</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
<mo>,</mo>
<mi>n</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mi>n</mi>
</msub>
<mi>j</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>z</mi>
</mrow>
</msub>
<mi>i</mi>
</msup>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>x</mi>
<mi>z</mi>
<mo>,</mo>
<mi>s</mi>
<mi>L</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mrow>
<mi>s</mi>
<mi>L</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>x</mi>
<mi>z</mi>
<mo>,</mo>
<mi>s</mi>
<mi>H</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mrow>
<mi>s</mi>
<mi>H</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>x</mi>
<mi>z</mi>
<mo>,</mo>
<mi>n</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mi>n</mi>
</msub>
<mi>j</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>y</mi>
<mi>z</mi>
</mrow>
</msub>
<mi>i</mi>
</msup>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>y</mi>
<mi>z</mi>
<mo>,</mo>
<mi>s</mi>
<mi>L</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mrow>
<mi>s</mi>
<mi>L</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>y</mi>
<mi>z</mi>
<mo>,</mo>
<mi>s</mi>
<mi>H</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mrow>
<mi>s</mi>
<mi>H</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>F</mi>
<mrow>
<mi>y</mi>
<mi>z</mi>
<mo>,</mo>
<mi>n</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msubsup>
<msup>
<msub>
<mi>D</mi>
<mi>n</mi>
</msub>
<mi>j</mi>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>z</mi>
<mi>z</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>=</mo>
<msup>
<msub>
<mi>p</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>t</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>z</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>y</mi>
<mi>z</mi>
</mrow>
</msub>
<mi>j</mi>
</msup>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
式中:
σxx i、σyy i、σzz i、σxy i、σxz i、σyz i分别为三维空间中任意点i处水力裂缝产生的正向诱导应力与剪切诱导应力,MPa;
均为边界影响系数;
N为水力裂缝被划分的总单元个数;
DsL j、DsH j、Dn j分别为水力裂缝单元j上的切向与法向位移不连续量;
σzz j为水力裂缝单元j上受到的沿原始最小地应力方向(平行于z轴)的正应力边界值,MPa;
σxz j、σyz j均为水力裂缝单元j上受到的剪切应力边界值,一般为0,MPa;
pnet j为水力裂缝单元j上受到的净压力,MPa。
3.如权利要求1所述的一种低渗透油气藏直井体积压裂储层改造体积预测方法,其特征在于,所述步骤(2)计算压裂液滤失后的地层孔隙压力,具体计算公式如下:
低渗透油藏:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>P</mi>
<mi>i</mi>
</msup>
<mo>=</mo>
<msub>
<mi>P</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>&xi;</mi>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msqrt>
<mrow>
<mi>&mu;</mi>
<mo>+</mo>
<msup>
<mi>b</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&xi;</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>g</mi>
<mi>&pi;</mi>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>g</mi>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<msub>
<mi>&phi;&mu;</mi>
<mi>l</mi>
</msub>
<mi>c</mi>
</mrow>
<mrow>
<mi>k</mi>
<mi>t</mi>
</mrow>
</mfrac>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
低渗透气藏:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>P</mi>
<mi>i</mi>
</msup>
<mo>=</mo>
<msub>
<mi>P</mi>
<mi>o</mi>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>o</mi>
</msub>
<mo>)</mo>
</mrow>
<mfrac>
<mi>&eta;</mi>
<msub>
<mi>&eta;</mi>
<mi>L</mi>
</msub>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&eta;</mi>
<mi>L</mi>
</msub>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<mn>2</mn>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>o</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>t</mi>
</mrow>
<mrow>
<msub>
<mi>&mu;</mi>
<mi>l</mi>
</msub>
<msup>
<mi>&phi;L</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mi>a</mi>
<mi>r</mi>
<mi> </mi>
<mi>sinh</mi>
<mrow>
<mo>(</mo>
<msqrt>
<mfrac>
<mrow>
<msup>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>4</mn>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</msqrt>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中:
Pi为压裂液滤失后地层中任意点i处的孔隙压力,MPa;
P0为地层原始孔隙压力,MPa;
Pf为水力裂缝缝内液体压力,MPa;
λ、μ为任意滤失点i在三维水力裂缝椭球体坐标系中的坐标值,m2;
b为水力裂缝在井筒处的半缝高,m;
φ为天然裂缝孔隙度,无因次;
μl为压裂液粘度,mPa·s;
c为天然裂缝压缩系数,MPa-1;
k为储层渗透率,D;
t为滤失时间,s;
L为水力裂缝长度,m;
xi为任意点i在三维直角坐标系中沿缝长方向(x轴方向)的坐标值,m;
yi为任意点i在三维直角坐标系中沿缝高方向(y轴方向)的坐标值,m。
4.如权利要求1所述的一种低渗透油气藏直井体积压裂储层改造体积预测方法,其特征在于,所述步骤(3)计算压裂液滤失后的地层孔隙弹性应力,其计算公式如下:
式中:
Δσi为任意点i处地层孔隙压力变化产生的孔隙弹性应力,MPa;
A为孔弹性常数,无因次;
ν地层岩石泊松比,无因次;
为描述流体侵入的尺寸和形状的参数,无因次;
α为Biot系数,无因次;
h为天然裂缝半缝高,m。
5.如权利要求1所述的一种低渗透油气藏直井体积压裂储层改造体积预测方法,其特征在于,所述步骤(4)计算各应力场叠加后储层空间中的三向有效主应力的大小与方向,过程如下:
首先,将上述三种诱导应力场与原地应力场进行叠加:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mi>H</mi>
</msub>
<mo>+</mo>
<mi>&alpha;</mi>
<mi>P</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>&sigma;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mi>V</mi>
</msub>
<mo>+</mo>
<mi>&alpha;</mi>
<mi>P</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>&sigma;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>z</mi>
<mi>z</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>z</mi>
<mi>z</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mi>h</mi>
</msub>
<mo>+</mo>
<mi>&alpha;</mi>
<mi>P</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>&sigma;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>z</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>z</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>y</mi>
<mi>z</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>y</mi>
<mi>z</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
将上式代入下式并求解:
σi 3-(σxxT+σyyT+σzzT)σi 2+(σxxTσyyT+σxxTσzzT+σyyTσzzT-σxyT 2-σxzT 2-σyzT 2)σi
-(σxxTσyyTσzzT-σxxTσyzT 2-σyyTσzxT 2-σzzTσxyT 2+2σxyTσxzTσyzT)=0 (i=1,2,3)
从而得到各应力场叠加后储层空间中的三向有效主应力的大小;
然后,通过下式计算各应力场叠加后储层空间中的三向有效主应力的方向:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>arc</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<msqrt>
<mrow>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>B</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
</mrow>
</msqrt>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>arc</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>B</mi>
<mi>i</mi>
</msub>
<msqrt>
<mrow>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>B</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
</mrow>
</msqrt>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>z</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>arc</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<msqrt>
<mrow>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>B</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
</mrow>
</msqrt>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
<mi>T</mi>
</mrow>
</msub>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>y</mi>
<mi>z</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>z</mi>
<mi>T</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>z</mi>
<mi>T</mi>
</mrow>
</msub>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>y</mi>
<mi>z</mi>
<mi>T</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
<mi>T</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:
σxx为水力裂缝产生的沿原始最大水平地应力方向(平行于x轴方向)的诱导应力,具体各点处的应力值为σxx i,MPa;
σyy为水力裂缝产生的沿原始垂向地应力方向(平行于y轴方向)的诱导应力,具体各点处的应力值为σyy i,MPa;
σzz为水力裂缝产生的沿原始最小水平地应力方向(平行于z轴方向)的诱导应力,具体各点处的应力值为σzz i,MPa;
σxy、σxz、σyz均为水力裂缝产生的剪切诱导应力,具体各点处的应力值为σxy i、σxz i、σyz i,MPa;
P为压裂液滤失后的地层中的孔隙压力,具体各点处的应力值为Pi,MPa;
Δσ为压裂液滤失后的地层中的孔隙弹性应力,具体各点处的应力值为Δσi,MPa;
σh为原始最小水平地应力,MPa;
σH为原始最大水平地应力,MPa;
σV为原始垂向地应力,MPa;
σi(i=1,2,3)为以上三种应力场与原地应力场叠加以后储层空间中的三向有效主应力的大小,其中σ1>σ2>σ3,MPa;
βxi、βyi、βzi(i=1,2,3)分别为以上三种应力场与原地应力场叠加以后储层空间中的各向有效主应力σi(i=1,2,3)分别与x、y、z轴方向的夹角,0°~90°。
6.如权利要求1所述的一种低渗透油气藏直井体积压裂储层改造体积预测方法,其特征在于,所述步骤(5)计算储层空间中天然裂缝张开破裂判定系数M、天然裂缝剪切破裂区域判定系数S,从而预测体积压裂储层改造体积,过程如下:
式中:
σn为天然裂缝壁面受到的正应力,MPa;
为天然裂缝面法向矢量与应力场叠加后的最大有效主应力σ1的夹角,0°~90°;
θ为天然裂缝面法向矢量与应力场叠加后的最小有效主应力σ3的夹角,0°~90°;
α为天然裂缝倾角即天然裂缝面法向矢量与y轴的夹角,0°~90°;
为天然裂缝走向与原始最大水平地应力(x轴方向)的夹角,0°~90°;
τ为天然裂缝面受到的剪切应力,MPa;
kf为天然裂缝壁面摩擦系数,无因次;
σo为天然裂缝内聚力,MPa;
储层空间中S>0和M>0的总波及体积即为直井体积压裂储层改造体积。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710280797.9A CN106869896B (zh) | 2017-04-26 | 2017-04-26 | 一种低渗透油气藏直井体积压裂储层改造体积预测方法 |
PCT/CN2017/082411 WO2018195909A1 (zh) | 2017-04-26 | 2017-04-28 | 一种低渗透油气藏直井体积压裂储层改造体积预测方法 |
US16/069,173 US10627543B2 (en) | 2017-04-26 | 2017-04-28 | Method for predicting reservoir reform volume after vertical well volume fracturing of low-permeability oil/gas reservoir |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710280797.9A CN106869896B (zh) | 2017-04-26 | 2017-04-26 | 一种低渗透油气藏直井体积压裂储层改造体积预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106869896A CN106869896A (zh) | 2017-06-20 |
CN106869896B true CN106869896B (zh) | 2018-05-08 |
Family
ID=59161734
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710280797.9A Active CN106869896B (zh) | 2017-04-26 | 2017-04-26 | 一种低渗透油气藏直井体积压裂储层改造体积预测方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US10627543B2 (zh) |
CN (1) | CN106869896B (zh) |
WO (1) | WO2018195909A1 (zh) |
Families Citing this family (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109507723A (zh) * | 2017-09-15 | 2019-03-22 | 中国石油化工股份有限公司 | 微地震压裂裂缝模型的压裂体积计算方法及*** |
CN107832524B (zh) * | 2017-11-07 | 2021-01-01 | 中国石油天然气股份有限公司 | 一种确定储层改造体积的方法、装置及*** |
CN108280275B (zh) * | 2018-01-09 | 2021-07-06 | 中国石油大学(华东) | 一种致密砂岩水力压裂缝高预测方法 |
CN108561106B (zh) * | 2018-03-12 | 2020-09-08 | 中国石油天然气股份有限公司 | 一种人工缝控储量提高采收率的油气开采方法 |
CN108984877A (zh) * | 2018-07-03 | 2018-12-11 | 中国石油化工股份有限公司胜利油田分公司石油工程技术研究院 | 一种致密储层复杂压裂裂缝形成条件的评估方法 |
CN108805365B (zh) * | 2018-07-05 | 2022-04-22 | 西南石油大学 | 一种非均质页岩储层裂缝诱导应力预测方法 |
CN109162701B (zh) * | 2018-09-12 | 2021-08-10 | 西南石油大学 | 一种煤层裸眼井破裂压力预测方法 |
CN109723435B (zh) * | 2018-12-12 | 2022-07-22 | 长江大学 | 确定页岩气藏水平井重复压裂前现地应力场的方法及设备 |
CN110298057B (zh) * | 2019-04-04 | 2022-04-05 | 西南石油大学 | 一种超临界二氧化碳压裂裂缝延伸计算方法 |
CN110175723B (zh) * | 2019-05-31 | 2020-04-21 | 西南石油大学 | 一种页岩气井重复压裂裂缝起裂方向与延伸路径预测方法 |
CN110490241B (zh) * | 2019-08-12 | 2022-05-06 | 中国石油大学(北京) | 一种水平井参数优化方法及装置 |
CN110929413B (zh) * | 2019-12-04 | 2020-11-10 | 西南石油大学 | 一种模拟裂缝性储层中水力裂缝壁面不连续性行为的方法 |
CN112505759B (zh) * | 2020-03-10 | 2022-09-09 | 中海油能源发展股份有限公司 | 一种低渗砂岩储层含气饱和度分布的预测方法 |
CN111396038B (zh) * | 2020-04-22 | 2023-10-03 | 中国石油天然气集团有限公司 | 一种水平单井产量的预测方法 |
CN111472752B (zh) * | 2020-05-27 | 2020-11-10 | 西南石油大学 | 一种利用录井数据的井漏裂缝宽度随钻估算方法 |
CN113803041B (zh) * | 2020-06-12 | 2023-06-20 | 中国石油化工股份有限公司 | 一种体积压裂设计方法及*** |
CN111502593B (zh) * | 2020-06-17 | 2020-11-06 | 西南石油大学 | 一种缝内暂堵转向压裂相变暂堵剂用量确定方法 |
CN111859260A (zh) * | 2020-07-03 | 2020-10-30 | 中国科学技术大学 | 一种利用压裂施工数据获得srv区域的方法 |
CN111927421B (zh) * | 2020-08-18 | 2022-07-15 | 中国石油化工股份有限公司 | 一种基于页岩气不稳定线性流理论的重复压裂选井方法 |
CN112127882B (zh) * | 2020-11-02 | 2021-05-25 | 西南石油大学 | 一种裂缝性地层钻井液漏失动态裂缝宽度计算方法 |
CN114508334B (zh) * | 2020-11-17 | 2024-05-31 | 中国石油化工股份有限公司 | 一种基于三维地应力场分布的溶洞循缝沟通技术确定方法 |
CN114592840B (zh) * | 2020-12-04 | 2023-10-27 | 中国石油天然气股份有限公司 | 暂堵压裂方法及其应用 |
US11960046B2 (en) * | 2021-01-22 | 2024-04-16 | Saudi Arabian Oil Company | Method for determining in-situ maximum horizontal stress |
CN112989721B (zh) * | 2021-03-04 | 2021-09-07 | 东北石油大学 | 一种致密油藏体积压裂水平井改造体积的快速计算方法 |
CN113420264B (zh) * | 2021-07-20 | 2022-07-05 | 西安石油大学 | 一种垂直井井壁破裂压力计算方法、装置和计算机设备 |
CN113417632B (zh) * | 2021-08-02 | 2023-03-24 | 中国石油大学(北京) | 基于压裂后压力确定煤储层原始地层压力的方法及装置 |
CN114441337B (zh) * | 2022-01-14 | 2022-10-18 | 中国石油大学(北京) | 一种金属断裂韧性裂纹尖端张开位移的检测方法 |
CN115577645B (zh) * | 2022-12-08 | 2023-04-18 | 中国石油大学(华东) | 一种燃爆压裂破裂范围预测模型的构建方法及预测方法 |
CN117077572B (zh) * | 2023-10-16 | 2024-01-26 | 西安石油大学 | 一种页岩油储层多簇裂缝扩展均匀程度的定量表征方法 |
CN117494446B (zh) * | 2023-11-09 | 2024-06-04 | 西南石油大学 | 一种判别二维水力裂缝与天然裂相遇后延伸行为的方法 |
CN117648886B (zh) * | 2024-01-25 | 2024-04-16 | 新疆石油管理局有限公司 | 高压脉冲射流储层改造压力预测方法、装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009063189A2 (en) * | 2007-11-13 | 2009-05-22 | Halliburton Energy Services, Inc. | Methods for geomechanical fracture modeling |
CN104278980A (zh) * | 2013-07-09 | 2015-01-14 | 中国石油天然气股份有限公司 | 一种采用改造体积优化致密油水平井缝网参数的方法 |
CN104712301A (zh) * | 2015-02-14 | 2015-06-17 | 中国石油集团西部钻探工程有限公司 | 一种确定火山岩储层压裂规模的方法 |
EP2966602A1 (en) * | 2014-07-11 | 2016-01-13 | Services Petroliers Schlumberger | Regional stress inversion using frictional faults |
CN106321051A (zh) * | 2015-07-01 | 2017-01-11 | 中国石油化工股份有限公司 | 一种用于优化多段压裂水平井网络裂缝参数的方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150205006A1 (en) * | 2010-03-25 | 2015-07-23 | Schlumberger Technology Corporation | Downhole modeling using inverted pressure and regional stress |
US20140358510A1 (en) * | 2013-05-29 | 2014-12-04 | Chevron U.S.A. Inc. | System and method for characterizing uncertainty in subterranean reservoir fracture networks |
CN105735960B (zh) * | 2016-03-22 | 2017-05-17 | 西南石油大学 | 一种低渗透油气藏水平井分段多簇压裂簇间距优化方法 |
-
2017
- 2017-04-26 CN CN201710280797.9A patent/CN106869896B/zh active Active
- 2017-04-28 US US16/069,173 patent/US10627543B2/en not_active Expired - Fee Related
- 2017-04-28 WO PCT/CN2017/082411 patent/WO2018195909A1/zh active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009063189A2 (en) * | 2007-11-13 | 2009-05-22 | Halliburton Energy Services, Inc. | Methods for geomechanical fracture modeling |
CN104278980A (zh) * | 2013-07-09 | 2015-01-14 | 中国石油天然气股份有限公司 | 一种采用改造体积优化致密油水平井缝网参数的方法 |
EP2966602A1 (en) * | 2014-07-11 | 2016-01-13 | Services Petroliers Schlumberger | Regional stress inversion using frictional faults |
CN104712301A (zh) * | 2015-02-14 | 2015-06-17 | 中国石油集团西部钻探工程有限公司 | 一种确定火山岩储层压裂规模的方法 |
CN106321051A (zh) * | 2015-07-01 | 2017-01-11 | 中国石油化工股份有限公司 | 一种用于优化多段压裂水平井网络裂缝参数的方法 |
Also Published As
Publication number | Publication date |
---|---|
US20190377101A1 (en) | 2019-12-12 |
US10627543B2 (en) | 2020-04-21 |
WO2018195909A1 (zh) | 2018-11-01 |
CN106869896A (zh) | 2017-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106869896B (zh) | 一种低渗透油气藏直井体积压裂储层改造体积预测方法 | |
CN112561144B (zh) | 一种致密油压裂水平井产能主控因素评判与产能预测方法 | |
CN109441422B (zh) | 一种页岩气井间距优化开采方法 | |
CN109522634B (zh) | 一种致密气多段体积压裂水平井数值分析方法 | |
CN105735960B (zh) | 一种低渗透油气藏水平井分段多簇压裂簇间距优化方法 | |
WO2016192077A1 (zh) | 一种致密气压裂水平井数值试井模型建立求解方法 | |
Du et al. | Modeling hydraulic fracturing induced fracture networks in shale gas reservoirs as a dual porosity system | |
CN113901681B (zh) | 一种全寿命周期页岩气储层双甜点三维可压性评估方法 | |
Cui et al. | Multidomain two-phase flow model to study the impacts of hydraulic fracturing on shale gas production | |
CN109324345B (zh) | 叠合盆地油气成藏期岩石孔隙度恢复方法 | |
CN104865610B (zh) | 一种深部储层岩心的数值表征方法 | |
CN114755310B (zh) | 一种裂缝性储层岩石力学层演变规律预测方法 | |
Zhang et al. | Pressure–dependent fracture permeability of marine shales in the Northeast Yunnan area, Southern China | |
CN112727534B (zh) | 基于真三向应力和渗透率动态变化的瓦斯钻孔布置方法 | |
CN105822298A (zh) | 基于产气指数的页岩气层绝对无阻流量的获取方法 | |
CN110656915B (zh) | 一种页岩气多段压裂水平井多工作制度产能预测方法 | |
CN115267905B (zh) | 一种复杂构造区钻井工程中裂缝性井漏预测方法 | |
Xi et al. | Research on the influence of strike-slip fault slippage on production casing and control methods and engineering application during multistage fracturing in deep shale gas wells | |
CN110568160A (zh) | 一种油气储层岩石的综合评价方法及装置 | |
CN107965316B (zh) | 一种提高高瓦斯低渗透单一煤层抽采效果的方法 | |
CN109211745A (zh) | 一种富含有机质泥页岩孔隙度的演化过程的恢复方法 | |
CN110863814B (zh) | 巨厚型生屑灰岩油藏单井分段比采液指数确定方法及装置 | |
CN112560246A (zh) | 一种目标井散点地层压力系数的预测方法 | |
Zhou et al. | Fluid flow through porous sandstone with overprinting and intersecting geological structures of various types | |
Zhang et al. | Staged fracturing of horizontal wells in continental tight sandstone oil reservoirs: a case study of Yanchang Formation in Western Ordos Basin, China |
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 |