CN113090258A - 基于测井数据的深层页岩气水平井压裂差异化设计方法 - Google Patents
基于测井数据的深层页岩气水平井压裂差异化设计方法 Download PDFInfo
- Publication number
- CN113090258A CN113090258A CN202110570418.6A CN202110570418A CN113090258A CN 113090258 A CN113090258 A CN 113090258A CN 202110570418 A CN202110570418 A CN 202110570418A CN 113090258 A CN113090258 A CN 113090258A
- Authority
- CN
- China
- Prior art keywords
- calculating
- evaluation index
- horizontal
- index
- shale gas
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000013461 design Methods 0.000 title claims abstract description 27
- 230000004069 differentiation Effects 0.000 title claims abstract description 13
- 238000011156 evaluation Methods 0.000 claims abstract description 84
- 239000011435 rock Substances 0.000 claims abstract description 63
- 238000011161 development Methods 0.000 claims abstract description 62
- 238000010276 construction Methods 0.000 claims abstract description 24
- 230000008569 process Effects 0.000 claims abstract description 13
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 12
- 230000018109 developmental process Effects 0.000 claims description 61
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 44
- 239000011707 mineral Substances 0.000 claims description 44
- 206010017076 Fracture Diseases 0.000 claims description 36
- 208000010392 Bone Fractures Diseases 0.000 claims description 34
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 21
- 229910052799 carbon Inorganic materials 0.000 claims description 21
- 239000011575 calcium Substances 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 12
- 230000000704 physical effect Effects 0.000 claims description 11
- 239000012530 fluid Substances 0.000 claims description 9
- 230000035699 permeability Effects 0.000 claims description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 239000004927 clay Substances 0.000 claims description 8
- 235000013305 food Nutrition 0.000 claims description 7
- 239000011148 porous material Substances 0.000 claims description 6
- 229920006395 saturated elastomer Polymers 0.000 claims description 6
- 239000002734 clay mineral Substances 0.000 claims description 4
- 229910001576 calcium mineral Inorganic materials 0.000 claims description 3
- 229910001748 carbonate mineral Inorganic materials 0.000 claims description 3
- 229910052570 clay Inorganic materials 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 230000027648 face development Effects 0.000 claims description 3
- 229910052683 pyrite Inorganic materials 0.000 claims description 3
- 239000011028 pyrite Substances 0.000 claims description 3
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 claims description 3
- 239000000758 substrate Substances 0.000 claims description 3
- 230000001965 increasing effect Effects 0.000 abstract description 5
- 239000007789 gas Substances 0.000 description 43
- 230000000694 effects Effects 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 238000012351 Integrated analysis Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000035558 fertility Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 235000013310 margarine Nutrition 0.000 description 1
- 239000003264 margarine Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- 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
-
- 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
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Theoretical Computer Science (AREA)
- Geophysics (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)
Abstract
本发明公开了基于测井数据的深层页岩气水平井压裂差异化设计方法,包括步骤:S1,计算水平井段的岩石力学参数、储层物性参数;S2,对步骤S1中仅用纵横波计算的岩石力学参数进行修正,计算得到修正后的岩石力学参数;S3,计算地层弱面发育参数;S4,计算目标深层页岩气井水平段的综合脆性指数、缝网建造能力指数和水平段资源丰度评价指数;S5,计算目标深层页岩气井水平段的“地质‑工程”综合评价指数,并进行曲线绘制;S6,对目标深层页岩气井水平段进行综合可压性分级,并根据综合可压性分级后的等级为不同水平井段优选最优的压裂工艺方法等;本发明可显著降低深层页岩气水平井缝网压裂成本、提高施工效率,改善增产效果。
Description
技术领域
本发明涉及非常规油气增产改造领域,更为具体的,涉及基于测井数据的深层页岩气水平井压裂差异化设计方法。
背景技术
我国深层页岩气资源丰富,是天然气资源增产的重要潜力之一,主要分布在四川盆地南部和东南部、涪陵平桥以及荆门等地区。但由于深层页岩处于高温、高压以及高应力差环境下,三向应力复杂、水平应力差大、温敏塑性强、以及储层非均质,压裂形成复杂缝网难度大等一系列问题,给深层页岩气的开发带来巨大的困难与挑战。因此,需要对储层做出准确认识,对不同水平段进行合理压裂工艺设计,实现压裂的高效性、合理性。
目前,国内外学者关于深层页岩气缝网压裂开展了研究,其中包括页岩的岩性、脆性指数、杨氏模量、泊松比、地应力差异、天然裂缝、层理、孔隙度、断裂韧性、成岩作用以及应力应变曲线等进行分析储层改造效果。陶祖文(题目:四川盆地永川区块五峰组—龙马溪组深层页岩可压性评价,期刊:天然气技术与经济,发表时间:2020)分析岩石力学参数以及相应的成分分析,建立考虑脆性指数、水平差应力系数、成岩作用、断裂韧性对深层页岩可压性的影响。曾波(题目:天然裂缝对页岩气水平井压裂的影响及工艺调整,期刊:中国石油勘探,发表时间:2020)综合分析已压裂井压裂效果的基础上,结合深层页岩储层地质工程特点,采用室内试验与数值模拟相结合的方式,优化了压裂工艺和关键参数,以提高缝网复杂程度、增大裂缝改造体积、维持裂缝长期导流能力。侯冰(题目:页岩暂堵转向压裂水力裂缝扩展物模试验研究,期刊:辽宁石油化工大学学报,发表时间:2020)发现页岩纹层面与水力主裂缝面距离较近时,可实现对纹层面的大面积沟通,纹层面张开或发生滑移,形成裂缝网络。
近年来,页岩气开发规模迅速扩大。然而,现有在现场应用的压裂工艺方法未能很好的针对不同区域、不同井段使用不同的压裂方法,仅是大规模使用比较单一的压裂方法进行压裂施工,成本较高但效果不理想。而且目前对影响缝网建造能力因素考虑不足,对岩石脆性评价指标较为单一,未能反映出矿物含量变化引起的岩石致密程度变化及内在机理。除此以外,忽略了岩石综合脆性与缝网建造能力之间的本质关系,更未考虑综合可压性与资源丰度的相关性,故地质-工程综合评价指数是深层页岩气差异化压裂设计的基础。
因此,亟需建立一种基于测井数据的深层页岩气水平井差异化压裂设计方法,同过对深层页岩储层复杂缝网形成评价与表征,根据不同井段综合可压性指数不同,优选出最优的压裂工艺方法。
发明内容
本发明的目的在于克服现有技术的不足,提供基于测井数据的深层页岩气水平井压裂差异化设计方法,可显著降低深层页岩气水平井缝网压裂成本、提高施工效率,改善增产效果。
本发明的目的是通过以下方案实现的:
基于测井数据的深层页岩气水平井压裂差异化设计方法,包括步骤:
S1,基于目标深层页岩气井的测井数据,计算水平井段的岩石力学参数、储层物性参数;
S2,基于测井数据中的物性参数,对步骤S1中仅用纵横波计算的岩石力学参数进行修正,建立考虑测井数据中的物性参数的影响计算模型,并通过该影响计算模型计算得到修正后的岩石力学参数;
S3,基于岩性与天然弱面参数之间的关系,基于矿物组分计算地层弱面发育参数;
S4,根据步骤S2和步骤S3计算得到的修正后的岩石力学参数、地层弱面发育参数,计算目标深层页岩气井水平段的综合脆性指数、缝网建造能力指数,根据步骤S1计算得到的储层物性参数,计算水平段资源丰度评价指数;
S5,基于步骤S4中计算得到的综合脆性指数、缝网建造能力指数、资源丰度评价指数,计算目标深层页岩气井水平段的“地质-工程”综合评价指数,并进行曲线绘制;
S6,基于步骤S5中绘制的“地质-工程”综合评价指数曲线,对目标深层页岩气井水平段进行综合可压性分级,并根据综合可压性分级后的等级为不同水平井段优选最优的压裂工艺方法。
进一步地,在步骤S1中,计算水平井段的岩石力学参数包括计算泊松比、计算弹性模量的步骤;在计算泊松比的步骤中,利用声波测井获取的横波速度、纵波速度按照如下公式计算泊松比:
在计算弹性模量的步骤中,利用声波测井获取的横波速度、纵波速度、岩石密度按照如下公式计算弹性模量:
式中:ν为岩石泊松比,无量纲;vp为声波测井纵波速度,单位m/s;vs为声波测井横波速度,单位m/s;E为岩石弹性模量,单位Pa;ρ为岩石密度,单位kg/m3。
进一步地,在步骤S1中,计算水平井段的储层物性参数包括计算储层孔隙度、计算总有机碳含量TOC的步骤;在计算储层孔隙度的步骤中,利用测井数据中的中子孔隙度、密度孔隙度按照如下公式计算得到储层孔隙度:
φ=φCNL·Pφ+φρ·(1-Pφ)
在计算总有机碳含量TOC的步骤中,利用测井数据中的密度孔隙度按照如下公式计算得到总有机碳含量:
TOC=α1φρ+α2fClay+α3fCa+α4fSi
式中:φ为储层孔隙度,%;φCNL为测井中子孔隙度,%;Pφ为孔隙度计算比例因子,无量纲;φρ分别为测井密度孔隙度,%;TOC为总有机碳含量,%;φρ分别为测井密度孔隙度,%;fClay、fCa、fSi分别为黏土矿物、钙质矿物、硅质矿物组分含量,%;α1、α2、α3、α4为总有机碳含量计算系数,无量纲。
进一步地,在步骤S2中,包括如下子步骤:
S21,按照如下公式计算页岩基质等效弹性力学参数:
Ke=(KV+KR)/2
Ge=(GV+GR)/2
式中:Ke、Ge分别为V-R-H模型计算得到的骨架岩石等效体积模量和切变模量,MPa;i为不同矿物种类,包括硅质矿物、碳酸盐矿物、黄铁矿及黏土等;Ki、Gi分别为第i类岩矿成分的基质岩石弹性模量,MPa;KR、GR分别为Reuss平均弹性模量,MPa;KV、GV分别为Voight平均弹性模量,MPa;fi为第i类脆性矿物含量,%;
S22,按照如下公式计算饱和气体岩石体积模量Ksat、剪切模量Gsat:
Ksat=Ke(1-α)+α2p,Gsat=Ge(1-α)
式中:p为随体积应变和流量增大的压力,MPa;α为Biot系数;φe为有效孔隙度,%;Km、Gm、Kf分别为岩石基质体积模量、切变模量与孔隙气体体积模量,MPa;
然后将储层的气水两相均视为流体,利用含水饱和度Sw、含气饱和度Sg,Sw=100%-Sg,按照如下公式计算得到孔隙流体体积模量:
式中:Kf为裂缝渗透率,KW为水渗透率,Kg为气体渗透率;
S23,修正后纵横波波速计算,建立饱和流体岩石体积模量Ksat、剪切模量Gsat、密度ρsat对其纵横波vp、vs速度表达式:
式中:ρ为岩石密度,g/cm3;
S24,根据修正后计算得到的纵横波,计算出目标储层岩石力学参数杨氏模量和泊松比。
进一步地,在步骤S3中,计算地层弱面发育参数包括计算天然裂缝发育强度和计算水平层理发育强度的步骤;
在计算天然裂缝发育强度的步骤中,通过对天然裂缝发育强度与硅质矿物组分含量、钙质矿物组分含量关系进行回归,按照如下公式计算得到天然裂缝发育强度:
BF=ξSifSi+ξCafCa+ξ
式中:BF为天然裂缝发育强度,m-1;ξSi为硅质矿物对天然裂缝发育贡献系数,m-1;ξCa为钙质矿物对天然裂缝发育贡献系数,m-1;ξ为天然裂缝发育本底系数,m-1;
在计算水平层理发育强度的步骤中,对水平层理发育强度与硅质矿物、钙质矿物含量关系进行回归,按照如下公式计算得到水平层理发育强度:
式中:BD为水平层理发育强度,m-1;ζSi为硅质矿物对水平层理发育贡献系数,m-1;ζCa为钙质矿物对水平层理发育贡献系数,m-1;ζ为水平层理发育本底系数,m-1。
进一步地,在步骤S4中,在计算综合脆性指数的步骤中,按照如下公式计算得到综合脆性指数FM:
式中:FM为综合脆性指数,无量纲;AM为基质脆性,无量纲;AMmax、AMmin、Emax、Emin、vmax、vmin分别为基质脆性、杨氏模量、泊松比的最大值、最小值;
在计算缝网建造能力指数的步骤中,按照如下公式计算得到缝网建造能力指数FNF:
式中:FNF为缝网建造能力指数,无量纲;BD、BF分别为水平层理发育指数和天然裂缝发育指数;BDmax、BDmin、BFmax、BFmin分别为水平层理和天然裂缝发育强度的最大值、最小值;
在计算资源丰度评价指数的步骤中,按照如下公式计算得到资源丰度评价指数FG:
FG=a1φ'+b1T'
式中:a1,b1为物性权重系数;Tmax,Tmin最大、最小TOC含量,%;To为等效总有机碳含量,%;T′为无量纲化的总有机碳/饱和度含量;φmax,φmin为工区最大、最小孔隙度,%;φe为有效孔隙度,%;φ′为无量纲化的孔隙度。
进一步地,在步骤S5中,计算目标深层页岩气井水平段的“地质-工程”综合评价指数包括计算地质评价指数CG、计算工程评价指数CE和计算工程评价指数门限极值或地质评价指数门限极值Υi;其中,在计算地质评价指数CG的步骤中,将资源丰度评价指数FG作为地质评价指数CG,即CG=FG;
工程评价指数CE则按照如下公式计算得到:
CE=λFFNF+(1-λF)FM
式中:λF为缝网因子权重系数,无量纲;
按照如下公式计算工程评价指数门限极值或地质评价指数门限极值Υi:
式中:i代表工程或者地质;Ci max、Ci min为工程或地质评价指数最大值与最小值,无量纲;ψi为工程评价指数或地质评价指数门限极值权重系数,无量纲。
进一步地,在步骤S6中,通过将计算出的地质评价指数和工程评价指数分别与工程评价指数门限极值和地质评价指数门限极值对比,对目标深层页岩气井水平段按照如下表格中关系进行差异化压裂设计:
分级依据 | 分级指标 | 差异化设计方法 |
(C<sub>G</sub>≥Υ<sub>G</sub>)且(C<sub>E</sub>≥Υ<sub>E</sub>) | 优 | 缝网压裂 |
(C<sub>G</sub>≥Υ<sub>G</sub>)且(C<sub>E</sub><Υ<sub>E</sub>) | 良 | 密切割缝内转向压裂 |
(C<sub>G</sub><Υ<sub>G</sub>)且(C<sub>E</sub>≥Υ<sub>E</sub>) | 中 | 长段段内转向压裂 |
(C<sub>G</sub><Υ<sub>G</sub>)且(C<sub>E</sub><Υ<sub>E</sub>) | 差 | 不压裂 |
。
本发明的有益效果包括:
本发明可显著降低深层页岩气水平井缝网压裂成本、提高施工效率,改善增产效果。具体的,本发明充分考虑深层页岩储层的“地质-工程”特征对压裂缝网形成的共同影响,基于目标水平井的测井数据,并通过孔隙度、矿物组分等物性参数,对传统仅用纵横波计算的岩石力学参数进行修正,最终获得修正后的岩石力学参数、地层弱面发育参数和储层物性参数,建立综合脆性指数、缝网建造能力指数以及资源丰度评价指数,进而计算出深层页岩储层水平段“地质-工程”综合评价指数,并对水平井段进行分级,优选出井段所需压裂工艺。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的方法流程图;
图2为本发明实施案例井水平段地质评价指数曲线图;
图3为本发明实施案例井水平段工程评价指数曲线图;
图4为本发明实施案例井水平段差异化设计曲线图。
具体实施方式
本说明书中所有实施例公开的所有特征,或隐含公开的所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以以任何方式组合和/或扩展、替换。
如图1~4所示,基于测井数据的深层页岩气水平井压裂差异化设计方法,已知深层页岩气井常规测井数据,包括声波横波速度、声波纵波速度、岩石密度、黏土矿物组分含量、硅质矿物组分含量、钙质矿物组分含量、中子孔隙度、密度孔隙度等数据。
基于测井数据的深层页岩气水平井压裂差异化设计方法,包括步骤:
S1,基于目标深层页岩气井的测井数据,计算水平井段的岩石力学参数、储层物性参数;
在步骤S1中,计算水平井段的岩石力学参数包括计算泊松比、计算弹性模量的步骤;在计算泊松比的步骤中,利用声波测井获取的横波速度、纵波速度按照如下公式计算泊松比:
在计算弹性模量的步骤中,利用声波测井获取的横波速度、纵波速度、岩石密度按照如下公式计算弹性模量:
式中:ν为岩石泊松比,无量纲;vp为声波测井纵波速度,单位m/s;vs为声波测井横波速度,单位m/s;E为岩石弹性模量,单位Pa;ρ为岩石密度,单位kg/m3。
在步骤S1中,计算水平井段的储层物性参数包括计算储层孔隙度、计算总有机碳含量TOC的步骤;在计算储层孔隙度的步骤中,利用测井数据中的中子孔隙度、密度孔隙度按照如下公式计算得到储层孔隙度:
φ=φCNL·Pφ+φρ·(1-Pφ)
在计算总有机碳含量TOC的步骤中,利用测井数据中的密度孔隙度按照如下公式计算得到总有机碳含量:
TOC=α1φρ+α2fClay+α3fCa+α4fSi
式中:φ为储层孔隙度,%;φCNL为测井中子孔隙度,%;Pφ为孔隙度计算比例因子,无量纲;φρ分别为测井密度孔隙度,%;TOC为总有机碳含量,%;φρ分别为测井密度孔隙度,%;fClay、fCa、fSi分别为黏土矿物、钙质矿物、硅质矿物组分含量,%;α1、α2、α3、α4为总有机碳含量计算系数,无量纲。
S2,基于测井数据中的物性参数,对步骤S1中仅用纵横波计算的岩石力学参数进行修正,建立考虑测井数据中的物性参数的影响计算模型,并通过该影响计算模型计算得到修正后的岩石力学参数;
在步骤S2中,包括如下子步骤:
S21,按照如下公式计算页岩基质等效弹性力学参数:
Ke=(KV+KR)/2
Ge=(GV+GR)/2
式中:Ke、Ge分别为V-R-H模型计算得到的骨架岩石等效体积模量和切变模量,MPa;i为不同矿物种类,包括硅质矿物、碳酸盐矿物、黄铁矿及黏土等;Ki、Gi分别为第i类岩矿成分的基质岩石弹性模量,MPa;KR、GR分别为Reuss平均弹性模量,MPa;KV、GV分别为Voight平均弹性模量,MPa;fi为第i类脆性矿物含量,%;
S22,按照如下公式计算饱和气体岩石体积模量Ksat、剪切模量Gsat:
Ksat=Ke(1-α)+α2p,Gsat=Ge(1-α)
式中:p为随体积应变和流量增大的压力,MPa;α为Biot系数;φe为有效孔隙度,%;Km、Gm、Kf分别为岩石基质体积模量、切变模量与孔隙气体体积模量,MPa;
然后将储层的气水两相均视为流体,利用含水饱和度Sw、含气饱和度Sg,Sw=100%-Sg,按照如下公式计算得到孔隙流体体积模量:
式中:Kf为裂缝渗透率,KW为水渗透率,Kg为气体渗透率;
S23,修正后纵横波波速计算,建立饱和流体岩石体积模量Ksat、剪切模量Gsat、密度ρsat对其纵横波vp、vs速度表达式:
式中:ρ为岩石密度,g/cm3;
S24,根据修正后计算得到的纵横波,计算出目标储层岩石力学参数杨氏模量和泊松比。
S3,基于岩性与天然弱面参数之间的关系,基于矿物组分计算地层弱面发育参数;
在步骤S3中,计算地层弱面发育参数包括计算天然裂缝发育强度和计算水平层理发育强度的步骤;
在计算天然裂缝发育强度的步骤中,通过对天然裂缝发育强度与硅质矿物组分含量、钙质矿物组分含量关系进行回归,按照如下公式计算得到天然裂缝发育强度:
BF=ξSifSi+ξCafCa+ξ
式中:BF为天然裂缝发育强度,m-1;ξSi为硅质矿物对天然裂缝发育贡献系数,m-1;ξCa为钙质矿物对天然裂缝发育贡献系数,m-1;ξ为天然裂缝发育本底系数,m-1;
在计算水平层理发育强度的步骤中,对水平层理发育强度与硅质矿物、钙质矿物含量关系进行回归,按照如下公式计算得到水平层理发育强度:
式中:BD为水平层理发育强度,m-1;ζSi为硅质矿物对水平层理发育贡献系数,m-1;ζCa为钙质矿物对水平层理发育贡献系数,m-1;ζ为水平层理发育本底系数,m-1。
S4,根据步骤S2和步骤S3计算得到的修正后的岩石力学参数、地层弱面发育参数,计算目标深层页岩气井水平段的综合脆性指数、缝网建造能力指数,根据步骤S1计算得到的储层物性参数,计算水平段资源丰度评价指数;
在步骤S4中,在计算综合脆性指数的步骤中,按照如下公式计算得到综合脆性指数FM:
式中:FM为综合脆性指数,无量纲;AM为基质脆性,无量纲;AMmax、AMmin、Emax、Emin、vmax、vmin分别为基质脆性、杨氏模量、泊松比的最大值、最小值;
在计算缝网建造能力指数的步骤中,按照如下公式计算得到缝网建造能力指数FNF:
式中:FNF为缝网建造能力指数,无量纲;BD、BF分别为水平层理发育指数和天然裂缝发育指数;BDmax、BDmin、BFmax、BFmin分别为水平层理和天然裂缝发育强度的最大值、最小值;
在计算资源丰度评价指数的步骤中,按照如下公式计算得到资源丰度评价指数FG:
FG=a1φ'+b1T'
式中:a1,b1为物性权重系数;Tmax,Tmin最大、最小TOC含量,%;To为等效总有机碳含量,%;T′为无量纲化的总有机碳/饱和度含量;φmax,φmin为工区最大、最小孔隙度,%;φe为有效孔隙度,%;φ′为无量纲化的孔隙度。
S5,基于步骤S4中计算得到的综合脆性指数、缝网建造能力指数、资源丰度评价指数,计算目标深层页岩气井水平段的“地质-工程”综合评价指数,并进行曲线绘制;
在步骤S5中,计算目标深层页岩气井水平段的“地质-工程”综合评价指数包括计算地质评价指数CG、计算工程评价指数CE和计算工程评价指数门限极值或地质评价指数门限极值Υi;其中,如图2所示,在计算地质评价指数CG的步骤中,将资源丰度评价指数FG作为地质评价指数CG,即CG=FG;
如图3所示,工程评价指数CE则按照如下公式计算得到:
CE=λFFNF+(1-λF)FM
式中:λF为缝网因子权重系数,无量纲;
按照如下公式计算工程评价指数门限极值或地质评价指数门限极值Υi:
式中:i代表工程或者地质;Ci max、Ci min为工程或地质评价指数最大值与最小值,无量纲;ψi为工程评价指数或地质评价指数门限极值权重系数,无量纲。
S6,基于步骤S5中绘制的“地质-工程”综合评价指数曲线,对目标深层页岩气井水平段进行综合可压性分级,并根据综合可压性分级后的等级为不同水平井段优选最优的压裂工艺方法。
在步骤S6中,通过将计算出的地质评价指数和工程评价指数分别与工程评价指数门限极值和地质评价指数门限极值对比,对目标深层页岩气井水平段按照如下表格中关系进行差异化压裂设计,如图4所示,如下表1、2:
表1深层页岩气水平井差异化压裂设计对应表
分级依据 | 分级指标 | 差异化设计方法 |
(C<sub>G</sub>≥Υ<sub>G</sub>)且(C<sub>E</sub>≥Υ<sub>E</sub>) | 优 | 缝网压裂 |
(C<sub>G</sub>≥Υ<sub>G</sub>)且(C<sub>E</sub><Υ<sub>E</sub>) | 良 | 密切割缝内转向压裂 |
(C<sub>G</sub><Υ<sub>G</sub>)且(C<sub>E</sub>≥Υ<sub>E</sub>) | 中 | 长段段内转向压裂 |
(C<sub>G</sub><Υ<sub>G</sub>)且(C<sub>E</sub><Υ<sub>E</sub>) | 差 | 不压裂 |
表2目标深层页岩气井水平井段压裂工艺表
深度/m | 分级指标 | 结论与压裂方案建议 |
3600-3769 | 差 | 不压裂 |
3769-3956 | 中 | 长段段内转向压裂 |
3956-4124 | 优 | 缝网压裂 |
4124-4288 | 良 | 密切割缝内转向压裂 |
4288-4430 | 优 | 缝网压裂 |
4430-4464 | 良 | 密切割缝内转向压裂 |
4464-4500 | 优 | 缝网压裂 |
4500-4624 | 良 | 密切割缝内转向压裂 |
4624-4867 | 优 | 缝网压裂 |
根据深层页岩气井水平段分级,为不同水平井段优选出针对性的压裂工艺方法:优——地质工程情况皆佳,可压性最好,优选缝网压裂;良——地质情况较好,工程情况稍差,优选密切割缝内转向压裂;中——地质情况稍差,工程情况较好差,可压性较II类次之,优选长段段内转向压裂;差——地质工程情况皆差,可压性最差,可不进行压裂。
本发明充分考虑深层页岩储层中“地质-工程”特征对压裂缝网形成的共同影响,基于目标水平井的测井数据,并通过孔隙度、矿物组分等物性参数,对传统仅用纵横波计算的岩石力学参数进行修正,最终获得修正后的岩石力学参数、地层弱面发育参数和储层物性参数,建立综合脆性指数、缝网建造能力指数以及资源丰度评价指数,进而计算出深层页岩储层水平段“地质-工程”综合评价指数,并对水平井段进行分级,优选出井段所需压裂工艺。
本发明未涉及部分均与现有技术相同或可采用现有技术加以实现。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
除以上实例以外,本领域技术人员根据上述公开内容获得启示或利用相关领域的知识或技术进行改动获得其他实施例,各个实施例的特征可以互换或替换,本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
本发明功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,在一台计算机设备(可以是个人计算机,服务器,或者网络设备等)以及相应的软件中执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、或者光盘等各种可以存储程序代码的介质,进行测试或者实际的数据在程序实现中存在于只读存储器(Random Access Memory,RAM)、随机存取存储器(Random Access Memory,RAM)等。
Claims (8)
1.基于测井数据的深层页岩气水平井压裂差异化设计方法,其特征在于,包括步骤:
S1,基于目标深层页岩气井的测井数据,计算水平井段的岩石力学参数、储层物性参数;
S2,基于测井数据中的物性参数,对步骤S1中仅用纵横波计算的岩石力学参数进行修正,建立考虑测井数据中的物性参数的影响计算模型,并通过该影响计算模型计算得到修正后的岩石力学参数;
S3,基于岩性与天然弱面参数之间的关系,基于矿物组分计算地层弱面发育参数;
S4,根据步骤S2和步骤S3计算得到的修正后的岩石力学参数、地层弱面发育参数,计算目标深层页岩气井水平段的综合脆性指数、缝网建造能力指数,根据步骤S1计算得到的储层物性参数,计算水平段资源丰度评价指数;
S5,基于步骤S4中计算得到的综合脆性指数、缝网建造能力指数、资源丰度评价指数,计算目标深层页岩气井水平段的“地质-工程”综合评价指数,并进行曲线绘制;
S6,基于步骤S5中绘制的“地质-工程”综合评价指数曲线,对目标深层页岩气井水平段进行综合可压性分级,并根据综合可压性分级后的等级为不同水平井段优选最优的压裂工艺方法。
3.根据权利要求1所述的基于测井数据的深层页岩气水平井压裂差异化设计方法,其特征在于,在步骤S1中,计算水平井段的储层物性参数包括计算储层孔隙度、计算总有机碳含量TOC的步骤;在计算储层孔隙度的步骤中,利用测井数据中的中子孔隙度、密度孔隙度按照如下公式计算得到储层孔隙度:
φ=φCNL·Pφ+φρ·(1-Pφ)
在计算总有机碳含量TOC的步骤中,利用测井数据中的密度孔隙度按照如下公式计算得到总有机碳含量:
TOC=α1φρ+α2fclay+α3fCa+α4fSi
式中:φ为储层孔隙度,%;φCNL为测井中子孔隙度,%;Pφ为孔隙度计算比例因子,无量纲;φρ分别为测井密度孔隙度,%;TOC为总有机碳含量,%;φρ分别为测井密度孔隙度,%;fC1ay、fCa、fSi分别为黏土矿物、钙质矿物、硅质矿物组分含量,%;α1、α2、α3、α4为总有机碳含量计算系数,无量纲。
4.根据权利要求1所述的基于测井数据的深层页岩气水平井压裂差异化设计方法,其特征在于,在步骤S2中,包括如下子步骤:
S21,按照如下公式计算页岩基质等效弹性力学参数:
Ke=(KV+KR)/2
Ge=(GV+GR)/2
式中:Ke、Ge分别为V-R-H模型计算得到的骨架岩石等效体积模量和切变模量,MPa;i为不同矿物种类,包括硅质矿物、碳酸盐矿物、黄铁矿及黏土等;Ki、Gi分别为第i类岩矿成分的基质岩石弹性模量,MPa;KR、GR分别为Reuss平均弹性模量,MPa;KV、GV分别为Voight平均弹性模量,MPa;fi为第i类脆性矿物含量,%;
S22,按照如下公式计算饱和气体岩石体积模量Ksat、剪切模量Gsat:
Ksat=Ke(1-α)+α2p,Gsat=Ge(1-α)
式中:p为随体积应变和流量增大的压力,MPa;α为Biot系数;φe为有效孔隙度,%;Km、Gm、Kf分别为岩石基质体积模量、切变模量与孔隙气体体积模量,MPa;
然后将储层的气水两相均视为流体,利用含水饱和度Sw、含气饱和度Sg,Sw=100%-Sg,按照如下公式计算得到孔隙流体体积模量:
式中:Kf为裂缝渗透率,KW为水渗透率,Kg为气体渗透率;
S23,修正后纵横波波速计算,建立饱和流体岩石体积模量Ksat、剪切模量Gsat、密度ρsat对其纵横波vp、vs速度表达式:
式中:ρ为岩石密度,g/cm3;
S24,根据修正后计算得到的纵横波,计算出目标储层岩石力学参数杨氏模量和泊松比。
5.根据权利要求1所述的基于测井数据的深层页岩气水平井压裂差异化设计方法,其特征在于,在步骤S3中,计算地层弱面发育参数包括计算天然裂缝发育强度和计算水平层理发育强度的步骤;
在计算天然裂缝发育强度的步骤中,通过对天然裂缝发育强度与硅质矿物组分含量、钙质矿物组分含量关系进行回归,按照如下公式计算得到天然裂缝发育强度:
BF=ξSifSi+ξCafCa+ξ
式中:BF为天然裂缝发育强度,m-1;ξSi为硅质矿物对天然裂缝发育贡献系数,m-1;ξCa为钙质矿物对天然裂缝发育贡献系数,m-1;ξ为天然裂缝发育本底系数,m-1;
在计算水平层理发育强度的步骤中,对水平层理发育强度与硅质矿物、钙质矿物含量关系进行回归,按照如下公式计算得到水平层理发育强度:
式中:BD为水平层理发育强度,m-1;ζSi为硅质矿物对水平层理发育贡献系数,m-1;ζCa为钙质矿物对水平层理发育贡献系数,m-1;ζ为水平层理发育本底系数,m-1。
6.根据权利要求1所述的基于测井数据的深层页岩气水平井压裂差异化设计方法,其特征在于,在步骤S4中,在计算综合脆性指数的步骤中,按照如下公式计算得到综合脆性指数FM:
式中:FM为综合脆性指数,无量纲;AM为基质脆性,无量纲;AMmax、AMmin、Emax、Emin、vmax、vmin分别为基质脆性、杨氏模量、泊松比的最大值、最小值;
在计算缝网建造能力指数的步骤中,按照如下公式计算得到缝网建造能力指数FNF:
式中:FNF为缝网建造能力指数,无量纲;BD、BF分别为水平层理发育指数和天然裂缝发育指数;BDmax、BDmin、BFmax、BFmin分别为水平层理和天然裂缝发育强度的最大值、最小值;
在计算资源丰度评价指数的步骤中,按照如下公式计算得到资源丰度评价指数FG:
FG=a1φ'+b1T'
式中:a1,b1为物性权重系数;Tmax,Tmin最大、最小TOC含量,%;To为等效总有机碳含量,%;T′为无量纲化的总有机碳/饱和度含量;φmax,φmin为工区最大、最小孔隙度,%;φe为有效孔隙度,%;φ′为无量纲化的孔隙度。
7.根据权利要求1所述的基于测井数据的深层页岩气水平井压裂差异化设计方法,其特征在于,在步骤S5中,计算目标深层页岩气井水平段的“地质-工程”综合评价指数包括计算地质评价指数CG、计算工程评价指数CE和计算工程评价指数门限极值或地质评价指数门限极值Υi;其中,在计算地质评价指数CG的步骤中,将资源丰度评价指数FG作为地质评价指数CG,即CG=FG;
工程评价指数CE则按照如下公式计算得到:
CE=λFFNF+(1-λF)FM
式中:λF为缝网因子权重系数,无量纲;
按照如下公式计算工程评价指数门限极值或地质评价指数门限极值Υi:
8.根据权利要求7所述的基于测井数据的深层页岩气水平井压裂差异化设计方法,其特征在于,在步骤S6中,通过将计算出的地质评价指数和工程评价指数分别与工程评价指数门限极值和地质评价指数门限极值对比,对目标深层页岩气井水平段按照如下表格中关系进行差异化压裂设计:
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110570418.6A CN113090258B (zh) | 2021-05-25 | 2021-05-25 | 基于测井数据的深层页岩气水平井压裂差异化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110570418.6A CN113090258B (zh) | 2021-05-25 | 2021-05-25 | 基于测井数据的深层页岩气水平井压裂差异化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113090258A true CN113090258A (zh) | 2021-07-09 |
CN113090258B CN113090258B (zh) | 2023-06-16 |
Family
ID=76664690
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110570418.6A Active CN113090258B (zh) | 2021-05-25 | 2021-05-25 | 基于测井数据的深层页岩气水平井压裂差异化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113090258B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113504569A (zh) * | 2021-09-09 | 2021-10-15 | 中国石油集团川庆钻探工程有限公司 | 阵列声波测井识别和评价岩体弱面的方法、系统、设备及介质 |
CN113622905A (zh) * | 2021-07-20 | 2021-11-09 | 中国地质大学(武汉) | 一种基于多因素综合分析的页岩储层脆性评价方法 |
CN114021466A (zh) * | 2021-11-10 | 2022-02-08 | 西南石油大学 | 基于返排数据和机器学习的页岩气有效缝网体积预测方法 |
CN114186440A (zh) * | 2022-02-16 | 2022-03-15 | 西南石油大学 | 一种地质-工程“双轨制”页岩可压性综合评价方法 |
CN115422718A (zh) * | 2022-08-10 | 2022-12-02 | 中海石油(中国)有限公司海南分公司 | 一种快速优选压裂层段的方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080159075A1 (en) * | 2006-12-28 | 2008-07-03 | Underhill William B | Technique and system for performing a cross well survey |
US20120203525A1 (en) * | 2011-02-08 | 2012-08-09 | Adrian Rodriguez Herrera | Three-dimensional modeling of parameters for oilfield drilling |
CN104832169A (zh) * | 2015-05-30 | 2015-08-12 | 重庆地质矿产研究院 | 水平井两井同步或异步多段分簇压裂室内实验井筒装置及方法 |
CN105719340A (zh) * | 2016-01-15 | 2016-06-29 | 西南石油大学 | 一种油气储层中人工压裂裂缝的三维建模方法 |
US20160222283A1 (en) * | 2013-01-04 | 2016-08-04 | Carbo Ceramics, Inc. | Proppant having non-uniform electrically conductive coatings and methods for making and using same |
CN108678772A (zh) * | 2018-07-17 | 2018-10-19 | 湖南科技大学 | 一种用于富水砂砾层巷道的多级层次支护装置及其施工方法 |
CN109931048A (zh) * | 2019-03-27 | 2019-06-25 | 南智(重庆)能源技术有限公司 | 油气井油套管完整性检测方法及评价系统 |
CN109964004A (zh) * | 2017-08-08 | 2019-07-02 | 中国科学院地质与地球物理研究所 | 一种油气田现场随钻岩屑扫描系统及方法 |
CN111425180A (zh) * | 2020-04-01 | 2020-07-17 | 中国矿业大学 | 一种ucg-ccs覆岩裂隙发育高度预测方法 |
CN112282751A (zh) * | 2020-12-01 | 2021-01-29 | 西南石油大学 | 一种地质工程三维耦合的致密油气水平井开采检测方法 |
-
2021
- 2021-05-25 CN CN202110570418.6A patent/CN113090258B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080159075A1 (en) * | 2006-12-28 | 2008-07-03 | Underhill William B | Technique and system for performing a cross well survey |
US20120203525A1 (en) * | 2011-02-08 | 2012-08-09 | Adrian Rodriguez Herrera | Three-dimensional modeling of parameters for oilfield drilling |
US20160222283A1 (en) * | 2013-01-04 | 2016-08-04 | Carbo Ceramics, Inc. | Proppant having non-uniform electrically conductive coatings and methods for making and using same |
CN104832169A (zh) * | 2015-05-30 | 2015-08-12 | 重庆地质矿产研究院 | 水平井两井同步或异步多段分簇压裂室内实验井筒装置及方法 |
CN105719340A (zh) * | 2016-01-15 | 2016-06-29 | 西南石油大学 | 一种油气储层中人工压裂裂缝的三维建模方法 |
CN109964004A (zh) * | 2017-08-08 | 2019-07-02 | 中国科学院地质与地球物理研究所 | 一种油气田现场随钻岩屑扫描系统及方法 |
CN108678772A (zh) * | 2018-07-17 | 2018-10-19 | 湖南科技大学 | 一种用于富水砂砾层巷道的多级层次支护装置及其施工方法 |
CN109931048A (zh) * | 2019-03-27 | 2019-06-25 | 南智(重庆)能源技术有限公司 | 油气井油套管完整性检测方法及评价系统 |
CN111425180A (zh) * | 2020-04-01 | 2020-07-17 | 中国矿业大学 | 一种ucg-ccs覆岩裂隙发育高度预测方法 |
CN112282751A (zh) * | 2020-12-01 | 2021-01-29 | 西南石油大学 | 一种地质工程三维耦合的致密油气水平井开采检测方法 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113622905A (zh) * | 2021-07-20 | 2021-11-09 | 中国地质大学(武汉) | 一种基于多因素综合分析的页岩储层脆性评价方法 |
CN113504569A (zh) * | 2021-09-09 | 2021-10-15 | 中国石油集团川庆钻探工程有限公司 | 阵列声波测井识别和评价岩体弱面的方法、系统、设备及介质 |
CN113504569B (zh) * | 2021-09-09 | 2021-12-03 | 中国石油集团川庆钻探工程有限公司 | 阵列声波测井识别和评价岩体弱面的方法、系统、设备及介质 |
CN114021466A (zh) * | 2021-11-10 | 2022-02-08 | 西南石油大学 | 基于返排数据和机器学习的页岩气有效缝网体积预测方法 |
CN114186440A (zh) * | 2022-02-16 | 2022-03-15 | 西南石油大学 | 一种地质-工程“双轨制”页岩可压性综合评价方法 |
CN114186440B (zh) * | 2022-02-16 | 2022-04-29 | 西南石油大学 | 一种地质-工程“双轨制”页岩可压性综合评价方法 |
US11598205B1 (en) | 2022-02-16 | 2023-03-07 | Southwest Petroleum University | Method for comprehensive evaluation of shale fracability under the geology-engineering “double-track” system |
CN115422718A (zh) * | 2022-08-10 | 2022-12-02 | 中海石油(中国)有限公司海南分公司 | 一种快速优选压裂层段的方法 |
CN115422718B (zh) * | 2022-08-10 | 2023-07-07 | 中海石油(中国)有限公司海南分公司 | 一种快速优选压裂层段的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113090258B (zh) | 2023-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113090258A (zh) | 基于测井数据的深层页岩气水平井压裂差异化设计方法 | |
CN110163533B (zh) | 一种页岩气缝网综合可压性评价方法 | |
Wang et al. | A novel experimental approach for fracability evaluation in tight-gas reservoirs | |
CN104007482B (zh) | 一种基于各向异性有效场的泥页岩岩石物理模型方法 | |
CN110714753A (zh) | 一种页岩储层可压性评价方法 | |
CN108009705A (zh) | 一种基于支持向量机技术的页岩储层可压性评价方法 | |
CN104268377B (zh) | 一种基于煤岩工业组分的脆性指数确定方法 | |
CN109632459B (zh) | 一种页岩可压裂性评价方法 | |
CN113625364B (zh) | 一种基于双重校正的泥页岩地层孔隙压力的计算方法 | |
CN110288233B (zh) | 一种基于模糊灰色关联法的深层页岩气可压性评价方法 | |
KR101902779B1 (ko) | 물리검층 자료를 이용한 s파 속도 산출방법 | |
Mokhtari et al. | Impacts of stress, natural and induced fractures on mechanical properties of organic-rich shales | |
Guo et al. | Classification and evaluation on shale gas reservoir for Wufeng-Longmaxi Formation in Chuannan Area, Sichuan Basin | |
CN111353218B (zh) | 一种煤层气-致密气储层合压性的测井定量评价方法 | |
Zhonghua et al. | Evaluation method of the least horizontal principal stress by logging data in anisotropic fast formations | |
CN112862259B (zh) | 一种海陆过渡相页岩储层立体动用开发评价方法 | |
CN108875115B (zh) | 一种确定岩石强度的方法 | |
CN113221232B (zh) | 用于裂缝性砂岩厚层水力压裂施工参数优化设计的方法 | |
CN115030714A (zh) | 一种致密储层工程甜点的识别方法 | |
CN116146176A (zh) | 基于地质-工程因素的深部煤层气储层测井定量分类方法 | |
CN110515124B (zh) | 页岩气储层水平井测井横波时差预测方法 | |
CN115387784B (zh) | 一种表征页岩油储层水平井改造工程甜点的评价方法 | |
CN115600354A (zh) | 一种非均质储层甜点识别评价方法 | |
Li et al. | Segmentation Study of Deep Shale Gas Horizontal Wells of the South Sichuan Shale Gas | |
Lv et al. | Determination of the red bed soft rock’s apparent preconsolidation pressure |
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 | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 100007 No. 9 North Main Street, Dongcheng District, Beijing, Dongzhimen Applicant after: PETROCHINA Co.,Ltd. Address before: 100011 Beijing city Dongcheng District No. 16 Andrew Applicant before: PETROCHINA Co.,Ltd. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |