CN106874544B - 一种页岩储层改造体积的地质表征方法 - Google Patents
一种页岩储层改造体积的地质表征方法 Download PDFInfo
- Publication number
- CN106874544B CN106874544B CN201710006624.8A CN201710006624A CN106874544B CN 106874544 B CN106874544 B CN 106874544B CN 201710006624 A CN201710006624 A CN 201710006624A CN 106874544 B CN106874544 B CN 106874544B
- Authority
- CN
- China
- Prior art keywords
- mpa
- fracture
- stress
- formula
- reservoir
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (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
本发明公开一种页岩储层改造体积的地质表征方法,包括以下步骤:获取各项储层物性参数、各项岩石力学参数;计算各个测井点对应的脆性因子Brit、新型脆性因子Bn、储层断裂韧性指数Kn、天然裂缝张开难易指数Pn、裂缝穿过评价指数Cn;建立各测井数据点对应的新型评价因子En,并判断;最后通化裂缝扩展波及长度Lo、波及高度Ho和波及宽度Wo获取储层改造体积。本发明的有益效果在于:依靠静态生产数据,采用测井等动态地质数据进行储层改造体积的预估,实现了压裂施工前对井网布署和压裂有利段的制定与选择,具有较强的实用性。
Description
技术领域
本发明涉及一种页岩储层改造体积的地质表征方法,属于油田开发技术领域。
背景技术
页岩储层改造体积(SRV)是单井获得高产的评价指标,对储层改造体积(SRV)的预测与表征具有重要意义。储层改造体积的表征和计算方法已从定性向定量发展,目前主要包括产量拟合法、半经验公式法、离散网络模型法、微地震监测法和半解析法。产量拟合法采用均值模型未考虑页岩储层天然裂缝发育引起的各向异性,且因页岩储层分段压裂方式不能得到单段产能数据,进而不能有效确定单段压裂SRV的展布范围。半经验公式法采用拉伸断裂Tx、Ty、剪切断裂Sx、Sy等经验常数,同样不能反映各水平井乃至水平井各段差异。离散网络模型法通过考虑压裂液、支撑剂性能和施工泵注程序等工程影响因素,模拟单次水力压裂过程中缝网的延伸过程进而确定SRV三维形态,但无法考虑多级或同步压裂过程中的相互干扰效应,并且不能在储层认识和压裂施工前对SRV进行预测和表征。微地震监测法结果较为可靠,但成本较高,不能做到对工区各单井进行监测。半解析法同样基于各向渗透率k、有效孔隙度φ等经验常数,缺乏变动性,且扩散系数需微地震数据校准。
缝网压裂技术作为页岩储层改造的重要手段,主要通过利用储层中天然裂缝等岩石弱面开启并连通形成复杂网络达到有效的页岩储层改造效果。对于储层可压裂性与压裂规模的预测等方面的工作应着重在页岩储层地质和压裂过程力学特征的分析上,而目前已有的研究则更多侧向压裂施工条件和工艺等方面,这可能因为忽略储层地质分析使得施工过程耗费大量财力和物力并可能破坏页岩储层。而仅仅通过地质和力学特征对页岩储层进行可压性评价却不能定量描述复杂缝网大致的延伸规模。因此,为使得在压裂施工时均取得较好的改造效果,必须建立在储层可压性评价基础上对复杂缝网延伸规模进行预测,即通过地质参数对储层改造体积进行定量表征。
页岩可压性评价主要通过岩心分析和测井解释结果,通过选取多个关键参数进行储层特征、力学性质的全面评价,适合现场操作并简单有效,亦可结合室内实验分析结果动静态转换;再通过连续可靠的地质数据嵌入岩石断裂力学理论中,消除该方法对高成本的微地震监测方法的依赖,最终对储层改造体积(SRV)进行地质表征。
发明内容
本发明所要解决的技术问题是提供一种克服半解析法经验常数的非确定性和对微地震的依赖,使得在压裂施工前对研究区各单井的改造程度有定量的认识,并有效预测改造体积的空间分布,为后续施工作为前期指导或施工印证提供可靠依据的页岩储层改造体积的地质表征方法
本发明解决上述技术问题所采用的技术方案是一种页岩储层改造体积的地质表征方法,该方法包括以下步骤:
S100、利用统计得到的岩心室内实验分析和单井测井解释结果获取各项储层物性参数、各项岩石力学参数;
S200、根据步骤S100得到的各项参数,计算各个测井数据点对应的脆性因子Brit、新型脆性因子Bn;
S300、再计算出各测井数据点对应的储层断裂韧性指数Kn、天然裂缝张开难易指数Pn、裂缝穿过评价指数Cn;
S400、通过下式建立各测井数据点对应的新型评价因子En;当En>0.5时,该测井数据点判定为有效数据点,进行下一步骤,否则判定为无效数据点;
式中:En为新型评价因子;
S500、计算步骤S400中有效数据点的裂缝扩展的波及长度Lo、波及高度Ho和波及宽度Wo,最后通过下式获取储层改造体积;
式中:VSRV为储层改造体积,其单位为m3;Lo为波及长度,m,Ho为波及高度,其单位为m,Wo为波及宽度,其单位为m,n为压裂段数;
所述步骤S100的具体过程为:
S101、对岩心室内实验分析结果处理,获取主力产气层单井垂直深度h、孔隙度φ、密度、硅质矿物组分含量VSi、碳酸盐岩矿物组分含量VCarb、岩石泊松比ν、杨氏模量E、体积模量K、剪切模量G、抗张强度St、孔隙压力Pp、垂向应力σz、水平最大主应力σH和水平最小主应力σh,单井测井解释得到无铀伽马γ、声波时差、中子CNL和密度ρ,所述声波时差为横波Δts、纵波Δtp;
S102、通过下式多元拟合获取孔隙度φ、泥质含量Vsh;
ρ=a1+b1ρlog
φlog=a2+b2γ+c2Δtp+d2CNL+e2ρlog
φ=a3+b3φlog
式中:ρ为静态密度值,其单位为g/cm3,a1、b1为密度拟合参数,ρlog为动态密度值,其单位为g/cm3,φlog为动态孔隙度值,其单位为小数,a2、b2、c2、d2、e2、a3、b3为孔隙度拟合参数,γ为无铀伽马值,其单位为API,CNL为中子值,其单位为%,Δts、Δtp为横波、纵波时差,其单位为s/m,φ为静态孔隙度值,其单位为小数;ISh为归一化无铀伽马值,γ为无铀伽马值,其单位为API,γmin、γmax为全区或单井测井无铀伽马最小值和最大值,其单位为API,GCUR为地层年代经验系数,新地层取3.7,老地层取2;
S103、通过下式多元拟合获取岩石力学参数;
式中:Cma为骨架体积压缩系数,其单位为MPa-1,Δtsma、Δtpma为岩石骨架横波时差、纵波时差,可由岩心实验得到,其单位为s,ρma为岩石骨架密度,可由岩心实验得到,其单位为kg/m3,Cb为体积压缩系数,即体积模量K的倒数,其单位为MPa-1,α为有效应力系数,ν为泊松比,Δts、Δtp为横波、纵波时差,其单位为s/m,ρ为静态密度值,其单位为g/cm3,E为杨氏模量,其单位为MPa,K为体积模量,其单位为MPa,G为第剪切模量,其单位为MPa;
Pp=αp×gh
St为抗张强度,其单位为MPa,Vsh为泥质含量,其单位为%,K为常数,取值12.26,Pp为孔隙压力,其单位为MPa,αp为区域压力系数,g为重力加速度,m/s2,h为垂直深度,其单位为m,σz为垂向主应力,其单位为MPa,σH、σh为水平最大和最小主应力,其单位为MPa,α为有效应力系数,β1、β2为水平最大、最小主应力方向的构造应力系数,可由岩心力学实验获取。
进一步的是,所述步骤S200的具体过程为:基于常规测井解释成果,提取脆性矿物组分的数据,通过下式获取各个测井数据点对应的脆性因子Brit、新型脆性因子Bn;
Brit=E/ν
式中:Vrit为矿物成分脆性指数,其单位为%,Vj为第j类脆性矿物含量,其单位为%,Brit为脆性因子,无量纲,Vmin、Vmax为全区或单井最小、最大矿物成分脆性指数,其单位为%,VRIT为归一化脆性矿物含量,无量纲,Bmin、Bmax为全区或单井最小、最大脆性因子,无量纲,BRIT为归一化脆性因子,无量纲,Bn为新型脆性因子,无量纲。
进一步的是,所述步骤S300的具体过程为:
S301、通过下式获取储层断裂韧性指数Kn;
KIC=0.217σn+0.0059St 3+0.0923St 2+0.517St-0.3322
KIIC=0.0956σh+0.1383St-0.082
KRIT=KICKIIC
式中:St为抗张强度,其单位为MPa,KIC、KIIC为I型、II型断裂韧性,KRIT为断裂韧性因子,无量纲,σh为最小主应力,MPa,KRITmax、KRITmin为全区或单井I型最大、最小断裂韧性因子,无量纲,Kn为断裂韧性指数;
S302、通过下式获取各测井数据点对应的天然裂缝张开难易指数Pn;
式中:σn为天然裂缝壁面受到的正应力,其单位为MPa,σz为垂向主应力,其单位为MPa,σH、σh为水平最大和最小主应力,其单位为MPa,ψ为天然裂缝与水平最大主应力的夹角,其单位为度,P为天然裂缝临界张开缝内压力,其单位为MPa,为天然裂缝壁面在三维空间的单位法向矢量,Pn为归一化的天然裂缝临界张开难易指数,其单位为MPa,Pmax、Pmin为区内或单井天然裂缝张开缝内最大、最小压力,其单位为MPa;
S303、通过下式获取各测井数据点对应的裂缝穿过评价指数Cn;
式中:ω为天然裂缝倾角,其单位为°,σH、σh为水平最大和最小主应力,其单位为MPa,σn为天然裂缝壁面受到的正应力,其单位为MPa,K0为天然裂缝穿过条件因子,大小为τ,方向向量为的剪应力矢量,大小为στ,方向向量为的剪应力矢量,So为天然裂缝内聚力,其单位为MPa,u为天然裂缝壁面摩擦系数,其系数为0~1,δ为水力裂缝与天然裂缝间的逼近角θ,τ为天然裂缝壁面剪应力大小,其单位为MPa,στ为天然裂缝壁面受水力裂缝尖端应力作用下的剪应力,其单位为MPa,Cn为天然裂缝穿过评价指数。
进一步的是,所述步骤S500的具体过程为:
S501、根据各单井储层地质特征设计单段泵注施工总液量Q和排量q、井底压力Pi、设计压裂液粘度μ、设计施工时间t和滤失系数;
S502、通过下式计算垂直于缝长的断面内,椭圆形断面的最大缝宽w;
式中:w为的椭圆断面最大缝宽,其单位为mm,p(x)为缝内压力,其单位为MPa,f为储层有效厚度与迭代波及高度之比,Δσ为水平主应力差,其单位为MPa;
S503、建立缝内压力流动方程和裂缝的波及高度Ho表征模型如下:
式中:q(x)为缝内流变量,其单位为m3/min,μ为压裂液粘度,其单位为mPa·s;Hp为储层有效厚度,其单位为m;
S504、建立关于波及长度Lo的流体连续性方程如下:
式中:λ(x,t)为压裂液滤失量,其单位为m3,tp(x)为从点x开始滤失的时间,其单位为s,Cx为滤失系数,其单位为m/min0.5,A(x,t)为裂缝断面面积,其单位为m2;
S505、建立y-z方向的二维裂缝诱导应力场如下:
式中:σ1、σ2、σ3为最大、中间和最小应力方向的诱导应力,其单位为MPa,Δσo为水平诱导应力差,其单位为MPa,r为裂缝中心距离诱导点A的距离,其单位为m,r1为裂缝底部距离诱导点A的距离,其单位为m,r2为裂缝顶部距离诱导点A的距离,其单位为m,θo为诱导点A偏离裂缝中心的角度,其单位为°,θ1为诱导点A偏离裂缝底部的角度,其单位为°,θ2为诱导点A偏离裂缝顶部的角度,其单位为°;
S506、通过诱导应力转向公式计算向非主应力方向的转向半径,从而得到相同生产模式下的各测井数据点对应的波及长度Lo、波及高度Ho、波及宽度Wo;
S507、最后通过下式获取储层改造体积;
式中:VSRV为储层改造体积,其单位为m3;Lo为波及长度,其单位为m,Ho为波及高度,其单位为m,Wo为波及宽度,其单位为m,n为压裂段数。
与现有技术相比,本发明的有益效果在于:首先,充分考虑了地层纵横向连续性变化的非均质性的影响,对不同地质特征的数据点逐一进行验证;其次,尽可能消除了对实验或高经济消耗的电测方法的应用,而是利用常规测井数据对物性和岩石力学基本特征进行计算,从而具有较广的实用度;进而,充分考虑各质点在压裂过程中的可行性,弥补了现有方法未能考虑单段不同区间可压性的不足;最后,依靠静态生产数据,采用测井等动态地质数据进行储层改造体积的预估,实现了压裂施工前对井网布署和压裂有利段的制定与选择,具有较强的实用性。
附图说明
图1为储层压裂改造规律示意图;
图2为实施例压裂后微地震裂缝监测平面图;
图3为实施例计算得到的裂缝三维分布图。
具体实施方式
下面结合实施例对本发明的具体实施方式做进一步的描述,并不因此将本发明限制在所述的实施例范围之中。
以四川盆地上奥陶统五峰组至下志留统龙马溪组某一区块的某单井页岩储层为例;
(1)获取其常规测井解释结果如表1所示;
表1水平段分段测井原始数据
分段 | 纵波时差 | 密度 | 无铀伽马 | 垂深 |
段 | tp(s) | ρ(g/cm3) | KTH | h(m) |
17 | 214.539 | 2.606 | 100.063 | 2648.597 |
16 | 212.987 | 2.587 | 97.299 | 2650.384 |
15 | 206.21 | 2.614 | 93.593 | 2649.778 |
14 | 208.822 | 2.586 | 87.315 | 2648.785 |
13 | 211.713 | 2.572 | 83.181 | 2648.656 |
12 | 213.806 | 2.547 | 80.07 | 2648.17 |
11 | 213.211 | 2.533 | 70.757 | 2647.431 |
10 | 215.992 | 2.515 | 69.452 | 2645.769 |
9 | 216.328 | 2.507 | 66.172 | 2645.216 |
8 | 218.421 | 2.511 | 54.8 | 2645.519 |
7 | 215.883 | 2.491 | 54.811 | 2643.677 |
6 | 212.19 | 2.521 | 57.676 | 2640.056 |
5 | 216.595 | 2.492 | 58.488 | 2636.624 |
4 | 216.487 | 2.5 | 61.063 | 2632.626 |
3 | 208.635 | 2.527 | 53.622 | 2628.203 |
2 | 216.358 | 2.5 | 70.585 | 2623.694 |
1 | 212.302 | 2.525 | 68.624 | 2619.715 |
(2)计算储层基本物性参数,采用表1数据,运用下列公式,孔隙度、脆性矿物含量均可通过多元拟合得到。其中对岩心室内实验分析结果处理,获取主力产气层单井垂直深度h,孔隙度密度,硅质矿物、碳酸盐岩矿物组分含量VSi,VCarb;再通过与测井数据对比建立多元拟合关系,将动态数据静态化,结果如表2所示。
公式为:ρ=a1+b1ρlog
φlog=a2+b2γ+c2Δtp+d2CNL+e2ρlog
φ=a3+b3φlog
式中:ρ为静态密度值,其单位为g/cm3,a1、b1为密度拟合参数,ρlog为动态密度值,其单位为g/cm3,φlog为动态孔隙度值,其单位为小数,a2、b2、c2、d2、e2、a3、b3为孔隙度拟合参数,γ为无铀伽马值,其单位为API,CNL为中子值,其单位为%,Δts、Δtp为横波、纵波时差,其单位为s/m,φ为静态孔隙度值,其单位为小数;ISh为归一化无铀伽马值,γ为无铀伽马值,其单位为API,γmin、γmax为全区或单井测井无铀伽马最小值和最大值,其单位为API,GCUR为地层年代经验系数,新地层取3.7,老地层取2,Cma为骨架体积压缩系数,其单位为MPa-1,Δtsma、Δtpma为岩石骨架横波时差、纵波时差,可由岩心实验得到,其单位为s,ρma为岩石骨架密度,可由岩心实验得到,其单位为kg/m3,Cb为第i个测井数据点对应的体积压缩系数,即体积模量K的倒数,其单位为MPa-1。
表2水平段分段测井物性数据
(3)计算岩石力学参数,同样通过表1数据,运用下列公式,计算各项岩石力学参数,如表3所示;其中,由于横波数据获取成本较高,可通过纵波时差与体积密度建立经验公式得到,再将其与进行偶极声波测井解释的结果进行拟合对比。对岩心室内实验分析结果处理,获取主力产气层岩石泊松比ν、杨氏模量E、体积模量K、剪切模量G、抗张强度St、孔隙压力Pp、垂向应力σz、水平最大主应力σH和水平最小主应力σh;利用单井测井解释声波时差和密度资料,计算岩石力学参数;
公式为:
式中:Cma为骨架体积压缩系数,MPa-1,Δtsma、Δtpma为岩石骨架横波时差、纵波时差,可由岩心实验得到,其单位为s,ρma为岩石骨架密度,可由岩心实验得到,其单位为kg/m3,Cb为体积压缩系数,即体积模量K的倒数,其单位为MPa-1,α为有效应力系数,ν为泊松比,Δts、Δtp为横波、纵波时差,其单位为s/m,ρ为静态密度值,其单位为g/cm3,E为杨氏模量,其单位为MPa,K为体积模量,其单位为MPa,G为剪切模量,其单位为MPa;
Pp=αp×gh
式中:St为抗张强度,其单位为MPa,Vsh为泥质含量,其单位为%,K为常数,取值12.26,Pp为孔隙压力,其单位为MPa,αp为区域压力系数,g为重力加速度,其单位为m/s2,h为垂直深度,其单位为m,σz为垂向主应力,其单位为MPa,σH、σh为水平最大和最小主应力,其单位为MPa,α为有效应力系数,β1、β2为水平最大、最小主应力方向的构造应力系数,可由岩心力学实验获取。
表3水平段分段测井岩石力学解释数据
(4)基于常规测井解释成果,提取脆性矿物组分的数据,通过下式获取各个测井数据点对应的脆性因子Brit、新型脆性因子Bn;
Brit=E/ν
式中:Vrit为矿物成分脆性指数,其单位为%,Vj为第j类脆性矿物含量,其单位为%,Brit为脆性因子,无量纲,Vmin、Vmax为全区或单井最小、最大矿物成分脆性指数,其单位为%,VRIT为归一化脆性矿物含量,无量纲,Bmin、Bmax为全区或单井最小、最大脆性因子,无量纲,BRIT为归一化脆性因子,无量纲,Bn为新型脆性因子,无量纲。
(5)通过下式获取储层断裂韧性指数Kn;
KIC=0.217σn+0.0059St 3+0.0923St 2+0.517St-0.3322
KIIC=0.0956σh+0.1383St-0.082
KRIT=KICKIIC
式中:St为抗张强度,其单位为MPa,KIC、KIIC为I型、II型断裂韧性,KRIT为第i个测井数据点对应的断裂韧性因子,无量纲,σh为最小主应力,MPa,KRITmax、KRITmin为全区或单井I型最大、最小断裂韧性因子,无量纲,Kn为断裂韧性指数;
(6)通过下式获取各测井数据点对应的天然裂缝张开难易指数Pn;
式中:σn为天然裂缝壁面受到的正应力,其单位为MPa,σz为垂向主应力,其单位为MPa,σH、σh为水平最大和最小主应力,其单位为MPa,ψ为天然裂缝与水平最大主应力的夹角,其单位为度,P为天然裂缝临界张开缝内压力,其单位为MPa,为天然裂缝壁面在三维空间的单位法向矢量,Pn为归一化的天然裂缝临界张开难易指数,其单位为MPa,Pmax、Pmin为区内或单井天然裂缝张开缝内最大、最小压力,其单位为MPa;
(7)通过下式获取各测井数据点对应的裂缝穿过评价指数Cn;
式中:ω为天然裂缝倾角,其单位为°,σH、σh为水平最大和最小主应力,其单位为MPa,σn为天然裂缝壁面受到的正应力,其单位为MPa,K0为天然裂缝穿过条件因子,大小为τ,方向向量为的剪应力矢量,大小为στ,方向向量为的剪应力矢量,So为天然裂缝内聚力,其单位为MPa,u为天然裂缝壁面摩擦系数,其系数为0~1,δ为水力裂缝与天然裂缝间的逼近角θ,τ为天然裂缝壁面剪应力大小,其单位为MPa,στ为天然裂缝壁面受水力裂缝尖端应力作用下的剪应力,其单位为MPa,Cn为天然裂缝穿过评价指数。
(8)通过下式建立各测井数据点对应的新型评价因子En;当En>0.5时,该测井数据点判定为有效数据点,即理存在可改造性,进行下一步骤,否则判定为无效数据点;其结果如表4:
式中:En为第i个测井数据点对应的新型评价因子;
表4水平段分段测井评价因子解释数据
结果表明从压裂段尺度而言均具备较好的可压性。
(9)根据各单井储层地质特征设计单段泵注施工总液量Q和排量q、井底压力Pi、设定施工参数:流体粘度μ为8mPa·s,施工时间t设定为3h,滤失系数取0.00071m/min0.5;
(10)通过下式计算垂直于缝长的断面内,椭圆形断面的最大缝宽w;
式中:w为椭圆断面最大缝宽,其单位为mm,p(x)为缝内压力,其单位为MPa,f为储层有效厚度与迭代波及高度之比,Δσ为水平主应力差,其单位为MPa;
(11)建立建立缝内压力流动方程和裂缝波及高度Ho表征模型如下:
式中:q(x)为缝内流变量,其单位为m3/min,μ为压裂液粘度,其单位为mPa·s;Hp为储层有效厚度,其单位为m;
(12)、建立关于波及长度Lo(即x)的流体连续性方程如下:
式中:λ(x,t)为压裂液滤失量,其单位为m3,tp(x)为从点x开始滤失的时间,其单位为s,Cx为滤失系数,其单位为m/min0.5,A(x,t)为裂缝断面面积,其单位为m2;
(13)建立y-z方向的二维裂缝诱导应力场如下:
式中:σ1、σ2、σ3为最大、中间和最小应力方向的诱导应力,其单位为MPa,Δσo为水平诱导应力差,其单位为MPa,r为裂缝中心距离诱导点A的距离,其单位为m,r1为裂缝底部距离诱导点A的距离,其单位为m,r2为裂缝顶部距离诱导点A的距离,其单位为m,θo为诱导点A偏离裂缝中心的角度,其单位为°,θ1为诱导点A偏离裂缝底部的角度,其单位为°,θ2为诱导点A偏离裂缝顶部的角度,其单位为°。
(14)通过诱导应力转向公式计算向非主应力方向的转向半径,从而得到相同生产模式下的各测井数据点对应的波及长度Lo、波及高度Ho、波及宽度Wo;
(15)最后通过下式获取储层改造体积;其结果如表5:
式中:VSRV为储层改造体积,其单位为m3;Lo为波及长度,其单位为m,Ho为波及高度,其单位为m,Wo为波及宽度,其单位为m,n为压裂段数。
表5地质表征SRV计算表
与现有技术相比,本发明的有益效果在于:首先,充分考虑了地层纵横向连续性变化的非均质性的影响,对不同地质特征的数据点逐一进行验证;其次,尽可能消除了对实验或高经济消耗的电测方法的应用,而是利用常规测井数据对物性和岩石力学基本特征进行计算,从而具有较广的实用度;进而,充分考虑各质点在压裂过程中的可行性,弥补了现有方法未能考虑单段不同区间可压性的不足;最后,依靠静态生产数据,采用测井等动态地质数据进行储层改造体积的预估,实现了压裂施工前对井网布署和压裂有利段的制定与选择,具有较强的实用性。
Claims (4)
1.一种页岩储层改造体积的地质表征方法,其特征在于,该方法包括以下步骤:
S100、利用统计得到的岩心室内实验分析和单井测井解释结果获取各项储层物性参数、各项岩石力学参数;
S200、根据步骤S100得到的各项参数,计算各个测井数据点对应的脆性因子Brit、新型脆性因子Bn;
S300、再计算出各测井数据点对应的储层断裂韧性指数Kn、天然裂缝张开难易指数Pn、裂缝穿过评价指数Cn;
S400、通过下式建立各测井数据点对应的新型评价因子En;当En>0.5时,该测井数据点判定为有效数据点,进行下一步骤,否则判定为无效数据点;
<mrow>
<msub>
<mi>E</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mroot>
<mrow>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msub>
<mi>K</mi>
<mi>n</mi>
</msub>
<msub>
<mi>P</mi>
<mi>n</mi>
</msub>
<msub>
<mi>C</mi>
<mi>n</mi>
</msub>
</mrow>
<mn>4</mn>
</mroot>
</mrow>
式中:En为新型评价因子;
S500、计算步骤S400中有效数据点的裂缝扩展的波及长度Lo、波及高度Ho和波及宽度Wo,最后通过下式获取储层改造体积;
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>S</mi>
<mi>R</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mi>&pi;</mi>
<mrow>
<mn>6</mn>
<msup>
<mi>n</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>L</mi>
<mi>o</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>W</mi>
<mi>o</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>H</mi>
<mi>o</mi>
</msub>
</mrow>
式中:VSRV为储层改造体积,其单位为m3;Lo为波及长度,其单位为m,Ho为波及高度,其单位为m,Wo为波及宽度,其单位为m,n为压裂段数;
所述步骤S100的具体过程为:
S101、对岩心室内实验分析结果处理,获取主力产气层单井垂直深度h、孔隙度φ、密度ρ、硅质矿物组分含量VSi、碳酸盐岩矿物组分含量VCarb、岩石泊松比ν、杨氏模量E、体积模量K、剪切模量G、抗张强度St、孔隙压力Pp、垂向应力σz、水平最大主应力σH和水平最小主应力σh,单井测井解释得到无铀伽马γ、声波时差、中子CNL和密度ρ,所述声波时差为横波Δts、纵波Δtp;
S102、通过下式多元拟合获取孔隙度φ、泥质含量Vsh;
ρ=a1+b1ρlog
φlog=a2+b2γ+c2Δtp+d2CNL+e2ρlog
φ=a3+b3φlog
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>S</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&gamma;</mi>
<mo>-</mo>
<msub>
<mi>&gamma;</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&gamma;</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msup>
<mn>2</mn>
<mrow>
<mi>G</mi>
<mi>C</mi>
<mi>U</mi>
<mi>R</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>S</mi>
<mi>h</mi>
</mrow>
</msub>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<msup>
<mn>2</mn>
<mrow>
<mi>G</mi>
<mi>C</mi>
<mi>U</mi>
<mi>R</mi>
</mrow>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
式中:ρ为静态密度值,其单位为g/cm3,a1、b1为密度拟合参数,ρlog为动态密度值,其单位为g/cm3,φlog为动态孔隙度值,其单位为小数,a2、b2、c2、d2、e2、a3、b3为孔隙度拟合参数,γ为无铀伽马值,其单位为API,CNL为中子值,其单位为%,Δts、Δtp为横波、纵波时差,其单位为s/m,φ为静态孔隙度值,其单位为小数;ISh为归一化无铀伽马值,γmin、γmax为全区或单井测井无铀伽马最小值和最大值,其单位为API,GCUR为地层年代经验系数,新地层取3.7,老地层取2;
S103、通过下式多元拟合获取岩石力学参数;
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<msup>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>s</mi>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>4</mn>
<msup>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>p</mi>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<mn>3</mn>
<msup>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>s</mi>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>&CenterDot;</mo>
<msup>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>p</mi>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>&alpha;</mi>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mi>b</mi>
</msub>
</mfrac>
</mrow>
<mrow>
<mi>v</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>0.5</mn>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;t</mi>
<mi>s</mi>
</msub>
<mo>/</mo>
<msub>
<mi>&Delta;t</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;t</mi>
<mi>s</mi>
</msub>
<mo>/</mo>
<msub>
<mi>&Delta;t</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>E</mi>
<mo>=</mo>
<mi>&rho;</mi>
<mo>&times;</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mo>&times;</mo>
<msup>
<msub>
<mi>&Delta;t</mi>
<mi>s</mi>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>4</mn>
<mo>&times;</mo>
<msup>
<msub>
<mi>&Delta;t</mi>
<mi>p</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<msub>
<mi>&Delta;t</mi>
<mi>s</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>&Delta;t</mi>
<mi>s</mi>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>&Delta;t</mi>
<mi>p</mi>
</msub>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>K</mi>
<mo>=</mo>
<mi>&rho;</mi>
<mfrac>
<mrow>
<mn>3</mn>
<msup>
<msub>
<mi>&Delta;t</mi>
<mi>s</mi>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>4</mn>
<msup>
<msub>
<mi>&Delta;t</mi>
<mi>p</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>3</mn>
<msup>
<msub>
<mi>&Delta;t</mi>
<mi>s</mi>
</msub>
<mn>2</mn>
</msup>
<mo>&CenterDot;</mo>
<msup>
<msub>
<mi>&Delta;t</mi>
<mi>p</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>G</mi>
<mo>=</mo>
<mfrac>
<mi>E</mi>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
式中:Cma为骨架体积压缩系数,其单位为MPa-1,Δtsma、Δtpma为岩石骨架横波时差、纵波时差,可由岩心实验得到,其单位为s,ρma为岩石骨架密度,可由岩心实验得到,其单位为kg/m3,Cb为体积压缩系数,即体积模量K的倒数,其单位为MPa-1,α为有效应力系数,ν为泊松比,Δts、Δtp为横波、纵波时差,其单位为s/m,ρ为静态密度值,其单位为g/cm3,E为杨氏模量,其单位为MPa,K为体积模量,其单位为MPa,G为剪切模量,其单位为MPa;
<mrow>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>0.0045</mn>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>0.008</mn>
<mi>E</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
</mrow>
<mi>K</mi>
</mfrac>
</mrow>
Pp=αp×gh
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>H</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>v</mi>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>v</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mi>z</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;P</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&alpha;P</mi>
<mi>p</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mi>z</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;P</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>h</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>v</mi>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>v</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mi>z</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;P</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&alpha;P</mi>
<mi>p</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mi>z</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;P</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
式中:St为抗张强度,其单位为MPa,Vsh为泥质含量,其单位为%,K为常数,取值12.26,Pp为孔隙压力,其单位为MPa,αp为区域压力系数,g为重力加速度,其单位为m/s2,h为垂直深度,其单位为m,σz为垂向主应力,其单位为MPa,σH、σh为水平最大和最小主应力,其单位为MPa,α为有效应力系数,β1、β2为水平最大、最小主应力方向的构造应力系数,可由岩心力学实验获取。
2.根据权利要求1所述的一种页岩储层改造体积的地质表征方法,其特征在于,所述步骤S200的具体过程为:基于常规测井解释成果,提取脆性矿物组分的数据,通过下式获取各个测井数据点对应的脆性因子Brit、新型脆性因子Bn;
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>r</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mi>j</mi>
<mi>n</mi>
</munderover>
<msub>
<mi>V</mi>
<mi>j</mi>
</msub>
</mrow>
Brit=E/ν
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>R</mi>
<mi>I</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>r</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>V</mi>
<mi>max</mi>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>B</mi>
<mrow>
<mi>R</mi>
<mi>I</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>B</mi>
<mrow>
<mi>r</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>B</mi>
<mi>max</mi>
</msub>
<mo>-</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>R</mi>
<mi>I</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>R</mi>
<mi>I</mi>
<mi>T</mi>
</mrow>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
式中:Vrit为矿物成分脆性指数,其单位为%,Vj为第j类脆性矿物含量,其单位为%,Brit为脆性因子,无量纲,Vmin、Vmax为全区或单井最小、最大矿物成分脆性指数,其单位为%,VRIT为归一化脆性矿物含量,无量纲,Bmin、Bmax为全区或单井最小、最大脆性因子,无量纲,BRIT为归一化脆性因子,无量纲,Bn为新型脆性因子,无量纲。
3.根据权利要求2所述的一种页岩储层改造体积的地质表征方法,其特征在于,所述步骤S300的具体过程为:
S301、通过下式获取储层断裂韧性指数Kn;
KIC=0.217σn+0.0059St 3+0.0923St 2+0.517St-0.3322
KIIC=0.0956σh+0.1383St-0.082
KRIT=KICKIIC
<mrow>
<msub>
<mi>K</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>R</mi>
<mi>I</mi>
<mi>T</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>R</mi>
<mi>I</mi>
<mi>T</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>R</mi>
<mi>I</mi>
<mi>T</mi>
<mi>max</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>R</mi>
<mi>I</mi>
<mi>T</mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
式中:St为抗张强度,其单位为MPa,KIC、KIIC为I型、II型断裂韧性,KRIT为断裂韧性因子,无量纲,σh为最小主应力,MPa,KRITmax、KRITmin为全区或单井I型最大、最小断裂韧性因子,无量纲,Kn为断裂韧性指数;
S302、通过下式获取各测井数据点对应的天然裂缝张开难易指数Pn;
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mover>
<mi>n</mi>
<mo>&RightArrow;</mo>
</mover>
<msup>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&sigma;</mi>
<mi>z</mi>
</msub>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<msub>
<mi>&sigma;</mi>
<mi>H</mi>
</msub>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<msub>
<mi>&sigma;</mi>
<mi>h</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mn>2</mn>
</msup>
<msup>
<mover>
<mi>n</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>P</mi>
<mo>=</mo>
<mrow>
<mo>|</mo>
<msub>
<mi>&sigma;</mi>
<mi>n</mi>
</msub>
<mo>|</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mover>
<mi>n</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&omega;</mi>
<mo>,</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&omega;</mi>
<mi>sin</mi>
<mi>&psi;</mi>
<mo>,</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&omega;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&psi;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>P</mi>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mi>max</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
式中:σn为天然裂缝壁面受到的正应力,其单位为MPa,σz为垂向主应力,其单位为MPa,σH、σh为水平最大和最小主应力,其单位为MPa,ψ为天然裂缝与水平最大主应力的夹角,其单位为度,P为天然裂缝临界张开缝内压力,其单位为MPa,为天然裂缝壁面在三维空间的单位法向矢量,Pn为天然裂缝临界张开难易指数,其单位为MPa,Pmax、Pmin为区内或单井天然裂缝张开缝内最大、最小压力,其单位为MPa;
S303、通过下式获取各测井数据点对应的裂缝穿过评价指数Cn;
<mrow>
<mi>&tau;</mi>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<msub>
<mi>&sigma;</mi>
<mi>v</mi>
</msub>
<mn>2</mn>
</msup>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&omega;</mi>
<mo>+</mo>
<msup>
<msub>
<mi>&sigma;</mi>
<mi>H</mi>
</msub>
<mn>2</mn>
</msup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&omega;sin</mi>
<mn>2</mn>
</msup>
<mi>&psi;</mi>
<mo>+</mo>
<msup>
<msub>
<mi>&sigma;</mi>
<mi>h</mi>
</msub>
<mn>2</mn>
</msup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&omega;cos</mi>
<mn>2</mn>
</msup>
<mi>&psi;</mi>
<mo>-</mo>
<msup>
<msub>
<mi>&sigma;</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>&tau;</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mi>h</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mfrac>
<mi>&delta;</mi>
<mn>2</mn>
</mfrac>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mfrac>
<mi>&delta;</mi>
<mn>2</mn>
</mfrac>
<mo>,</mo>
<mi>&delta;</mi>
<mo>=</mo>
<mi>&theta;</mi>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mo>|</mo>
<mo>|</mo>
<mover>
<mi>&tau;</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>+</mo>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>&tau;</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>o</mi>
</msub>
<mo>-</mo>
<mi>u</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mi>h</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>cos</mi>
<mn>3</mn>
</msup>
<mfrac>
<mi>&delta;</mi>
<mn>2</mn>
</mfrac>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mi>n</mi>
</msub>
<mo>&rsqb;</mo>
<mo><</mo>
<mn>0</mn>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0.5</mn>
<mo>,</mo>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
<mo>></mo>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<mn>2</mn>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<mo>,</mo>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
<mo><</mo>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中:ω为天然裂缝倾角,其单位为°,σH、σh为水平最大和最小主应力,其单位为MPa,σn为天然裂缝壁面受到的正应力,其单位为MPa,K0为天然裂缝穿过条件因子,大小为τ,方向向量为的剪应力矢量,大小为στ,方向向量为的剪应力矢量,So为天然裂缝内聚力,其单位为MPa,u为天然裂缝壁面摩擦系数,取值范围0~1,δ为水力裂缝与天然裂缝间的逼近角θ,τ为天然裂缝壁面剪应力大小,其单位为MPa,στ为天然裂缝壁面受水力裂缝尖端应力作用下的剪应力,其单位为MPa,Cn为天然裂缝穿过评价指数。
4.根据权利要求3所述的一种页岩储层改造体积的地质表征方法,其特征在于,所述步骤S500的具体过程为:
S501、根据各单井储层地质特征设计单段泵注施工总液量Q和排量q、井底压力Pi、设计压裂液粘度μ、设计施工时间t和滤失系数;
S502、通过下式计算垂直于缝长的断面内,椭圆形断面的最大缝宽w;
<mrow>
<mi>w</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>E</mi>
</mfrac>
<mo>{</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mn>2</mn>
<mi>&pi;</mi>
</mfrac>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&sigma;</mi>
</mrow>
<mrow>
<msub>
<mi>p</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mi>h</mi>
</msub>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msup>
<mi>cos</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>f</mi>
<mi> </mi>
<mi>l</mi>
<mi>n</mi>
<mfrac>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>f</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mi>f</mi>
</mfrac>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
式中:w为椭圆断面最大缝宽,其单位为mm,p(x)为缝内压力,其单位为MPa,f为储层有效厚度与迭代波及高度之比,Δσ为水平主应力差,其单位为MPa;
S503、建立缝内压力流动方程和裂缝的波及高度Ho表征模型如下:
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mn>64</mn>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>&mu;</mi>
</mrow>
<mrow>
<mi>&pi;</mi>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>w</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mn>64</mn>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>&mu;</mi>
</mrow>
<mrow>
<mi>&pi;</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>w</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msub>
<mi>K</mi>
<mrow>
<mi>I</mi>
<mi>C</mi>
</mrow>
</msub>
<msqrt>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>&Delta;&sigma;H</mi>
<mi>p</mi>
</msub>
</mrow>
<mrow>
<mi>&pi;</mi>
<msqrt>
<mrow>
<msup>
<mi>H</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<msub>
<mi>H</mi>
<mi>p</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中:q(x)为缝内流变量,其单位为m3/min,μ为压裂液粘度,其单位为mPa·s;Hp为储层有效厚度,其单位为m;
S504、建立关于波及长度Lo的流体连续性方程如下:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mo>-</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>&lambda;</mi>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&lambda;</mi>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
<mo>=</mo>
<mn>2</mn>
<mi>H</mi>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
<mfrac>
<msub>
<mi>C</mi>
<mi>x</mi>
</msub>
<msqrt>
<mrow>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>t</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
式中:λ(x,t)为压裂液滤失量,其单位为m3,tp(x)为从点x开始滤失的时间,其单位为s,Cx为滤失系数,其单位为m/min0.5,A(x,t)为裂缝断面面积,其单位为m2;
S505、建立y-z方向的二维裂缝诱导应力场如下:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>p</mi>
<mfrac>
<mi>r</mi>
<mi>H</mi>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msup>
<mi>H</mi>
<mn>2</mn>
</msup>
<mrow>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mn>3</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<msub>
<mi>sin&theta;</mi>
<mn>0</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mfrac>
<mn>3</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>p</mi>
<mo>&lsqb;</mo>
<mfrac>
<mi>r</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mfrac>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&sigma;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mi>p</mi>
<mfrac>
<mi>r</mi>
<mi>H</mi>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msup>
<mi>H</mi>
<mn>2</mn>
</msup>
<mrow>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mn>3</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<msub>
<mi>sin&theta;</mi>
<mn>0</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mfrac>
<mn>3</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>p</mi>
<mo>&lsqb;</mo>
<mfrac>
<mi>r</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mfrac>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&sigma;</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;&sigma;</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mi>p</mi>
<mfrac>
<mi>r</mi>
<mi>H</mi>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msup>
<mi>H</mi>
<mn>2</mn>
</msup>
<mrow>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mn>3</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<msub>
<mi>sin&theta;</mi>
<mn>0</mn>
</msub>
<mi>sin</mi>
<mfrac>
<mn>3</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mi>p</mi>
<mo>&lsqb;</mo>
<mfrac>
<mi>r</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mfrac>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mi>arctan</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mi>x</mi>
<mi>y</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>arctan</mi>
<mo>&lsqb;</mo>
<mfrac>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>y</mi>
<mo>-</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mo>,</mo>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mi>arctan</mi>
<mo>&lsqb;</mo>
<mfrac>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>H</mi>
<mo>-</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>r</mi>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>,</mo>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>+</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>,</mo>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
式中:σ1、σ2、σ3为最大、中间和最小应力方向的诱导应力,其单位为MPa,Δσo为水平诱导应力差,其单位为MPa,r为裂缝中心距离诱导点A的距离,其单位为m,r1为裂缝底部距离诱导点A的距离,其单位为m,r2为裂缝顶部距离诱导点A的距离,其单位为m,θo为诱导点A偏离裂缝中心的角度,其单位为°,θ1为诱导点A偏离裂缝底部的角度,其单位为°,θ2为诱导点A偏离裂缝顶部的角度,其单位为°;
S506、通过诱导应力转向公式计算向非主应力方向的转向半径,从而得到相同生产模式下的各测井数据点对应的波及长度Lo、波及高度Ho、波及宽度Wo;
S507、最后通过下式获取储层改造体积;
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>S</mi>
<mi>R</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mi>&pi;</mi>
<mrow>
<mn>6</mn>
<msup>
<mi>n</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>L</mi>
<mi>o</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>W</mi>
<mi>o</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>H</mi>
<mi>o</mi>
</msub>
</mrow>
式中:VSRV为储层改造体积,其单位为m3;Lo为波及长度,其单位为m,Ho为波及高度,其单位为m,Wo为波及宽度,其单位为m,n为压裂段数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710006624.8A CN106874544B (zh) | 2017-01-05 | 2017-01-05 | 一种页岩储层改造体积的地质表征方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710006624.8A CN106874544B (zh) | 2017-01-05 | 2017-01-05 | 一种页岩储层改造体积的地质表征方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106874544A CN106874544A (zh) | 2017-06-20 |
CN106874544B true CN106874544B (zh) | 2018-01-30 |
Family
ID=59165629
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710006624.8A Active CN106874544B (zh) | 2017-01-05 | 2017-01-05 | 一种页岩储层改造体积的地质表征方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106874544B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107705215A (zh) * | 2017-09-25 | 2018-02-16 | 西南石油大学 | 一种页岩储层重复压裂选井选段方法 |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109209356B (zh) * | 2017-07-06 | 2021-08-31 | 中国石油化工股份有限公司 | 一种基于张破裂和剪破裂确定地层可压裂性的方法 |
CN108019205B (zh) * | 2017-09-14 | 2021-04-30 | 中国石油天然气股份有限公司 | 一种确定储层可压裂性指数的方法及装置 |
CN107977480B (zh) * | 2017-10-18 | 2021-04-30 | 中国石油化工集团有限公司 | 一种页岩气储层产气性快速评价方法 |
CN107977505A (zh) * | 2017-11-28 | 2018-05-01 | 兰州大学 | 一种前期降水衰退系数k确定的新方法 |
CN108343433B (zh) * | 2018-02-28 | 2019-11-05 | 西南石油大学 | 用于计算页岩微裂缝缝宽变化下气体质量传输参数的方法 |
CN108645596B (zh) * | 2018-05-11 | 2019-12-27 | 北京石油化工学院 | 一种用于评价流体对岩石造缝能力的实验方法 |
CN110569512B (zh) * | 2018-06-05 | 2021-08-20 | 中国石油化工股份有限公司 | 一种泥页岩中渗透性夹层发育程度定量表征方法 |
CN111859603B (zh) * | 2019-04-30 | 2022-11-01 | 中国石油天然气股份有限公司 | 砂岩裂缝改造方式的评定方法和评定装置 |
CN110288233B (zh) * | 2019-06-26 | 2023-08-08 | 重庆科技学院 | 一种基于模糊灰色关联法的深层页岩气可压性评价方法 |
CN110864966B (zh) * | 2019-10-12 | 2022-08-05 | 中国石油天然气股份有限公司 | 一种适用于裂缝型岩石的可压性评价方法 |
CN113756793B (zh) * | 2020-06-05 | 2023-09-26 | 中国石油天然气股份有限公司 | 页岩油开采方式的确定方法、装置、设备及可读存储介质 |
CN111859260B (zh) * | 2020-07-03 | 2024-07-16 | 中国科学技术大学 | 一种利用压裂施工数据获得srv区域的方法 |
CN114109374B (zh) * | 2020-08-31 | 2023-07-25 | 中国石油天然气股份有限公司 | 页岩气储层靶窗位置的确定方法及装置 |
CN113033021A (zh) * | 2021-04-21 | 2021-06-25 | 北京艾迪博科油气技术有限公司 | 致密储层的改造评估方法和装置 |
CN115263272A (zh) * | 2022-07-19 | 2022-11-01 | 华美孚泰油气增产技术服务有限责任公司 | 通过测井确定水平最大主应力的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103437746A (zh) * | 2013-06-28 | 2013-12-11 | 东方宝麟科技发展(北京)有限公司 | 一种水平井多段段内多缝体积压裂方法 |
CN104775810A (zh) * | 2015-03-03 | 2015-07-15 | 西南石油大学 | 一种页岩气储层可压性评价方法 |
CN105114047A (zh) * | 2015-07-22 | 2015-12-02 | 西南石油大学 | 一种利用测井资料优选体积压裂射孔层位的方法 |
CN105550410A (zh) * | 2015-12-07 | 2016-05-04 | 西南石油大学 | 计算页岩储层水力压裂倾斜裂缝诱导应力的方法 |
-
2017
- 2017-01-05 CN CN201710006624.8A patent/CN106874544B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103437746A (zh) * | 2013-06-28 | 2013-12-11 | 东方宝麟科技发展(北京)有限公司 | 一种水平井多段段内多缝体积压裂方法 |
CN104775810A (zh) * | 2015-03-03 | 2015-07-15 | 西南石油大学 | 一种页岩气储层可压性评价方法 |
CN105114047A (zh) * | 2015-07-22 | 2015-12-02 | 西南石油大学 | 一种利用测井资料优选体积压裂射孔层位的方法 |
CN105550410A (zh) * | 2015-12-07 | 2016-05-04 | 西南石油大学 | 计算页岩储层水力压裂倾斜裂缝诱导应力的方法 |
Non-Patent Citations (4)
Title |
---|
《Simultaneous hydraulic fracturing of ultra-low permeability sandstone reservoirs in China: Mechanism and its field test》;REN Lan等;《中南大学学报(英文版)》;20150415;第22卷(第4期);第1427-1436页 * |
《低渗透油藏体积压裂数值模拟研究》;蔡田田;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20131215(第12期);第18-20页第4.2节 * |
《储层改造体积预测模型的研究与应用》;翁定为等;《石油钻探技术》;20160413;第44卷(第1期);第95-100页 * |
《页岩气储层可压性评价新方法》;赵金洲等;《天然气地球科学》;20150610;第26卷(第6期);第1165-1172页 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107705215A (zh) * | 2017-09-25 | 2018-02-16 | 西南石油大学 | 一种页岩储层重复压裂选井选段方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106874544A (zh) | 2017-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106874544B (zh) | 一种页岩储层改造体积的地质表征方法 | |
US10571605B2 (en) | Method and system for modeling and simulating a fractured reservoir | |
US20180355701A1 (en) | Hydraulic fracturing simulation | |
Wang et al. | Poroelastic and poroplastic modeling of hydraulic fracturing in brittle and ductile formations | |
CN108681635B (zh) | 一种致密储层体积压裂可压性评价方法 | |
Abdollahipour et al. | DEM simulation of confining pressure effects on crack opening displacement in hydraulic fracturing | |
Zou et al. | 3-D numerical simulation of hydraulic fracturing in a CBM reservoir | |
Gao et al. | Exploration of non-planar hydraulic fracture propagation behaviors influenced by pre-existing fractured and unfractured wells | |
CN112502701B (zh) | 一种低渗透储层综合地质-工程的分类评价方法 | |
Liu et al. | Geomechanical simulation of partially propped fracture closure and its implication for water flowback and gas production | |
CN108665105B (zh) | 一种导水裂缝带发育高度预计方法 | |
CN103278614A (zh) | 一种动静态岩石力学参数校正方法及装置 | |
EA014144B1 (ru) | Способы и устройства для прогнозирования добычи углеводородов с места заложения скважины | |
Feng et al. | Modeling of curving hydraulic fracture propagation from a wellbore in a poroelastic medium | |
Ji et al. | An assessment method for shale fracability based on fractal theory and fracture toughness | |
Zhang et al. | An evaluation method of volume fracturing effects for vertical wells in low permeability reservoirs | |
CN108319738A (zh) | 一种页岩气井产量预测方法 | |
CN105931125B (zh) | 一种致密油分段多簇体积压裂水平井产量预测方法 | |
Zhao et al. | Numerical study of the effect of natural fractures on shale hydraulic fracturing based on the continuum approach | |
Ju et al. | Adaptive finite element-discrete element method for numerical analysis of the multistage hydrofracturing of horizontal wells in tight reservoirs considering pre-existing fractures, hydromechanical coupling, and leak-off effects | |
He et al. | Parametric study and dimensional analysis on prescribed hydraulic fractures in cave mining | |
CN107290799B (zh) | 一种岩石可压性的确定方法 | |
CN116127675A (zh) | 一种页岩油水平井体积压裂最大可采储量的预测方法 | |
CN116128083A (zh) | 一种页岩油水平井体积压裂裂缝体积的定量表征方法 | |
Liu et al. | Rock crack propagation mechanism of oriented perforation hydraulic fracture under different perforation parameters |
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 |