CN109374757A - 应用声波幅度定量化处理评价火成岩储层有效性的方法 - Google Patents

应用声波幅度定量化处理评价火成岩储层有效性的方法 Download PDF

Info

Publication number
CN109374757A
CN109374757A CN201811280231.7A CN201811280231A CN109374757A CN 109374757 A CN109374757 A CN 109374757A CN 201811280231 A CN201811280231 A CN 201811280231A CN 109374757 A CN109374757 A CN 109374757A
Authority
CN
China
Prior art keywords
wave amplitude
ratio
section
stoneley
effective 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.)
Granted
Application number
CN201811280231.7A
Other languages
English (en)
Other versions
CN109374757B (zh
Inventor
张国强
谭忠健
郝仲田
张璋
李东
李欣
任宏
岳明亮
伍鹏
祁晓
刘如明
涂春赵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Tianjin Branch
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Tianjin Branch
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Offshore Oil Corp CNOOC, CNOOC China Ltd Tianjin Branch filed Critical China National Offshore Oil Corp CNOOC
Priority to CN201811280231.7A priority Critical patent/CN109374757B/zh
Publication of CN109374757A publication Critical patent/CN109374757A/zh
Application granted granted Critical
Publication of CN109374757B publication Critical patent/CN109374757B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/48Processing the detected response signal, e.g. electronic circuits specially adapted therefor by amplitude comparison
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04Analysing solids
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/024Mixtures
    • G01N2291/02441Liquids in porous solids
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/024Mixtures
    • G01N2291/0245Gases in porous solids

Landscapes

  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Signal Processing (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

一种应用声波幅度定量化处理评价火成岩储层有效性的方法:在深度区间内进行阵列声波测井;对阵列声波测井资料采用均方根法进行波幅计算;获取同口井同岩性条件下,有效储层段与非有效储层段的波幅计算值的比值;利用斯通利波、横波幅度比值指示的有效厚度计算波幅产能指数;利用斯通利波幅度比和横波幅度比,建立基于波幅比值的储层有效性评价图版;依据波幅产能指数、斯通利波及横波幅度比值将储层有效性划分为三个区间,分别为I类储层区、II类储层区及III类储层区;将单井声波测井资料处理后,根据处理结果在图版中所处区域,对储层有效性及产能进行快速预判。本发明能够解决火成岩储层有效性评价问题,应用效果好,适用度高,能够显著降低勘探开发成本。

Description

应用声波幅度定量化处理评价火成岩储层有效性的方法
技术领域
本发明涉及一种火成岩储层有效性的评价方法。特别是涉及一种应用阵列声波幅度定量化处理评价火成岩储层有效性的方法。
背景技术
随着勘探技术的进步,火成岩油气藏作为陆相断陷盆地隐蔽油气藏的重要组成,近年来已成为我国东部盆地深化勘探和增加储量的重要领域。渤海湾盆地多个区块均在火成岩地层中获得良好油气显示,表明渤海湾盆地火成岩储层潜力巨大。
目前,火成岩储层有效性评价主要依靠电成像及阵列声波资料。当地层中发育裂缝、孔洞时,由于泥浆侵入的影响,裂缝、孔洞往往具有暗色正弦线及暗色斑点的电成像测量表征。但同时,由于地层层理、层界面等地质体在电成像资料上与裂缝具有相同的表征,因此应用电成像资料识别裂缝存在较高程度的多解性。而当裂缝被泥质、黄铁矿等暗色矿物充填时,应用电成像资料也难以对其有效性进行评价;阵列声波资料虽然能够通过波幅度衰减对裂缝及孔洞有效性进行评价,但由于仪器、井况等的不同,往往只能进行单井分析,无法对多井资料进行联立对比应用,为井间储层有效性评价研究带来了难题。
目前针对上述解释评价需求,火成岩储层有效性评价方法存在较多的不确定性及复杂性,且无法实现井间对比,从而需要研究应用声波幅度定量化处理评价火成岩储层有效性的方法。
发明内容
本发明所要解决的技术问题是,提供一种能够有效提升火成岩储层有效性评价准确度的应用阵列声波幅度定量化处理评价火成岩储层有效性的方法。
本发明所采用的技术方案是:一种应用声波幅度定量化处理评价火成岩储层有效性的方法,包括如下步骤:
1)在给定深度区间内进行阵列声波测井,得到阵列声波全波列数据;
2)对阵列声波全波列数据分别采用均方根法进行全井段的横波幅度和斯通利波幅度计算;
3)在同口井同岩性条件下,将有效储层段的横波幅度和斯通利波幅度的计算结果与非有效储层段的横波幅度和斯通利波幅度的计算结果进行比值计算,获得横波幅度比值与斯通利波幅度比值;
4)以横波幅度比值为横坐标,以斯通利波幅度比值为纵坐标,建立基于波幅比值的储层有效性评价图版;
5)依据横波幅度比值、斯通利波幅度比值将有效性评价图版划分为三个区间,定义横波幅度比值小于0.7、斯通利波幅度比值小于0.3的区间为I类储层区;定义横波幅度比值小于0.8且大于或等于0.7、斯通利波幅度比值小于0.8且大于或等于0.3的区间为II类储层区;定义横波幅度比值大于或等于0.8、斯通利波幅度比值大于或等于0.8的区间为III类储层区;
6)重复步骤1)~步骤3)得到单井阵列声波全波列数据处理结果,根据处理结果在储层有效性评价图版中所处的区间,快速预判单井储层有效性,并进行井间储层有效性对比。
步骤2)是采用均方根法对阵列声波全波列数据进行波幅度计算:
其中,是任一接收器接收到的阵列声波全波列数据的均方根波幅度,Ten是计算终止窗位置,Tst是计算起始窗位置,x是模式波幅度,t是时间变量。
步骤3)是采用如下公式获得横波幅度比值:
其中,RatioAmpx是有效储层段横波幅度计算值与非有效储层段横波幅度计算值的比值,AmpxiE是N个接收器中的第i个接收器在有效储层段的横波幅度计算值,AmpxiN是N个接收器中的第i个接收器在非有效储层段的横波幅度计算值,N是参与横波幅度计算的接收器总数;
步骤3)是采用如下公式获得斯通利波幅度比值:
其中,RatioAmpst是有效储层段斯通利波幅度计算值与非有效储层段斯通利波幅度计算值的比值,AmpxiE是N个接收器中的第i个接收器在有效储层段的斯通利波幅度计算值,AmpxiN是N个接收器中的第i个接收器在非有效储层段的斯通利波幅度计算值,N是参与斯通利波幅度计算的接收器总数。
本发明的应用阵列声波幅度定量化处理评价火成岩储层有效性的方法,能够解决火成岩储层有效性评价存在多解性及复杂性的问题,通过对阵列声波幅度进行定量化处理进行储层有效性评价,能够将多井阵列声波资料进行对比应用,有效提升火成岩储层有效性评价准确度。本发明的方法应用效果好,适用度高,能够显著降低勘探开发成本,同时为后续类似地层条件下的储层评价工作提供了思路方法。
附图说明
图1是本发明应用阵列声波幅度定量化处理评价火成岩储层有效性的方法的流程图;
图2是利用横波幅度比值和斯通利波幅度比值建立的基于波幅比值的储层有效性评价图版;
图3a是渤海BZ区块某井火成岩段阵列声波幅度处理成果图;
图3b是采用本发明中方法处理该井阵列声波全波列资料所得成果图;
图4a是渤海PL区块某井火成岩段阵列声波幅度处理成果图;
图4b是采用本发明中方法处理该井阵列声波全波列资料所得成果图。
具体实施方式
下面结合实例和附图对本发明的应用阵列声波幅度定量化处理评价火成岩储层有效性的方法做出详细说明。
本发明的应用阵列声波幅度定量化处理评价火成岩储层有效性的方法,通过对阵列声波测井获取的阵列声波全波列资料进行声波幅度定量化处理分析,得到有效储层段与非有效储层段的声波幅度计算值的比值。应用横波幅度比值和斯通利波幅度比值建立基于波幅比值的储层有效性评价图版,根据横波波幅比值、斯通利波幅度比值在图版中所处的区域,对储层有效性进行快速预判,并进行井间储层有效性对比。
如图1所示,本发明的应用声波幅度定量化处理评价火成岩储层有效性的方法,包括如下步骤:
1)在给定深度区间内进行阵列声波测井,得到阵列声波全波列数据;
2)对阵列声波全波列数据分别采用均方根法进行全井段的横波幅度和斯通利波幅度计算;
采用均方根法对阵列声波全波列数据进行波幅度计算的公式如下:
其中,是任一接收器接收到的阵列声波全波列数据的均方根波幅度,Ten是计算终止窗位置,Tst是计算起始窗位置,x是模式波幅度,t是时间变量。
3)在同口井同岩性条件下,将有效储层段的横波幅度和斯通利波幅度的计算结果与非有效储层段的横波幅度和斯通利波幅度的计算结果进行比值计算,获得横波幅度比值与斯通利波幅度比值;其中,
采用公式(2)获得横波幅度比值:
其中,RatioAmpx是有效储层段横波幅度计算值与非有效储层段横波幅度计算值的比值,AmpxiE是N个接收器中的第i个接收器在有效储层段的横波幅度计算值,AmpxiN是N个接收器中的第i个接收器在非有效储层段的波横幅度计算值,N是参与横波幅度计算的接收器总数。
采用公式(3)获得斯通利波幅度比值:
其中,RatioAmpst是有效储层段斯通利波幅度计算值与非有效储层段斯通利波幅度计算值的比值,AmpxiE是N个接收器中的第i个接收器在有效储层段的斯通利波幅度计算值,AmpxiN是N个接收器中的第i个接收器在非有效储层段的斯通利波幅度计算值,N是参与斯通利波幅度计算的接收器总数。
4)以横波幅度比值为横坐标、以斯通利波幅度比值为纵坐标,建立基于波幅比值的储层有效性评价图版;
图2为渤海区块利用斯通利波幅度比值和横波幅度比值,建立的基于波幅比值的储层有效性评价图版;
5)依据横波幅度比值、斯通利波幅度比值将有效性评价图版划分为三个区间。定义横波幅度比值小于0.7、斯通利波幅度比值小于0.3的区间为I类储层区;定义横波幅度比值小于0.8且大于或等于0.7、斯通利波幅度比值小于0.8且大于或等于0.3的区间为II类储层区;定义横波幅度比值大于0.8、斯通利波幅度比值大于0.8的区间为III类储层区;
6)重复步骤1)~步骤3)得到单井阵列声波全波列数据处理结果,根据处理结果在储层有效性评价图版中所处的区间,快速预判单井储层有效性,并进行井间储层有效性对比。
下面,通过实际例子对本发明的应用声波幅度定量化处理评价火成岩储层有效性的方法作进一步说明。
实例1:图3a为渤海BZ区块某井火成岩段声波幅度处理成果图,图3b为采用本发明的方法处理该井资料所得成果图。该井图示段岩性主要为花岗岩,电成像观察该段孔洞、裂缝均有发育,但依据常规资料难以识别有效储层位置。波幅处理结果表明本井有效储层段主要集中于3932-3980m处(图3a实线框内),该段内横波幅度、斯通利波幅度降低明显,指示该段储层有效性较好。将该段进行波幅比值处理并投射至储层有效性评价图版,投射结果位于I类储层区内,如图3b所示,表明储层具备较好的产出能力。本井针对3879.00-3998.66m进行测试,采用6.35mm油嘴放喷,测试压差4.475Mpa,日产油100.48方,日产气107272方,与图版判断结果一致。
实例2:图4a为渤海PL区块某井火成岩段声波幅度处理成果图,图4b为采用本发明的方法处理该井资料所得成果图。该井图示段岩性主要为花岗岩,电成像资料表明本井储集空间以裂缝为主,未见孔洞,由于裂缝所占空间体积较小,因此依据常规资料难以识别有效储层位置。波幅处理结果在1500-1515m、1520-1540m处出现两段较明显衰减,常规资料表明1500-1515m段井眼扩径严重(图4a实线框内),声波幅度降低与井径扩径程度呈正比,因此推断声波幅度降低主要受扩径影响,不能反映储层有效性信息。1520-1540m段波幅度变化与井径扩径成反相关(图4a虚线框内),表明本段声波幅度受井眼扩径影响较小,对该段进行波幅比值计算并投影至储层有效性评价图版中,投影结果表明该有效储层段落在II类储层区内,如图4b所示,产出能力有限。本井针对1500.00-1565.00m进行测试,测试压差7.484Mpa,日产水12.3方,与图版判断结果一致。

Claims (3)

1.一种应用声波幅度定量化处理评价火成岩储层有效性的方法,其特征在于,包括如下步骤:
1)在给定深度区间内进行阵列声波测井,得到阵列声波全波列数据;
2)对阵列声波全波列数据分别采用均方根法进行全井段的横波幅度和斯通利波幅度计算;
3)在同口井同岩性条件下,将有效储层段的横波幅度和斯通利波幅度的计算结果与非有效储层段的横波幅度和斯通利波幅度的计算结果进行比值计算,获得横波幅度比值与斯通利波幅度比值;
4)以横波幅度比值为横坐标,以斯通利波幅度比值为纵坐标,建立基于波幅比值的储层有效性评价图版;
5)依据横波幅度比值、斯通利波幅度比值将有效性评价图版划分为三个区间,定义横波幅度比值小于0.7、斯通利波幅度比值小于0.3的区间为I类储层区;定义横波幅度比值小于0.8且大于或等于0.7、斯通利波幅度比值小于0.8且大于或等于0.3的区间为II类储层区;定义横波幅度比值大于或等于0.8、斯通利波幅度比值大于或等于0.8的区间为III类储层区;
6)重复步骤1)~步骤3)得到单井阵列声波全波列数据处理结果,根据处理结果在储层有效性评价图版中所处的区间,快速预判单井储层有效性,并进行井间储层有效性对比。
2.根据权利要求1所述的应用声波幅度定量化处理评价火成岩储层有效性的方法,其特征在于,步骤2)是采用均方根法对阵列声波全波列数据进行波幅度计算:
其中,是任一接收器接收到的阵列声波全波列数据的均方根波幅度,Ten是计算终止窗位置,Tst是计算起始窗位置,x是模式波幅度,t是时间变量。
3.根据权利要求1所述的应用声波幅度定量化处理评价火成岩储层有效性的方法,其特征在于,步骤3)是采用如下公式获得横波幅度比值:
其中,RatioAmpx是有效储层段横波幅度计算值与非有效储层段横波幅度计算值的比值,AmpxiE是N个接收器中的第i个接收器在有效储层段的横波幅度计算值,AmpxiN是N个接收器中的第i个接收器在非有效储层段的横波幅度计算值,N是参与横波幅度计算的接收器总数;
步骤3)是采用如下公式获得斯通利波幅度比值:
其中,RatioAmpst是有效储层段斯通利波幅度计算值与非有效储层段斯通利波幅度计算值的比值,AmpxiE是N个接收器中的第i个接收器在有效储层段的斯通利波幅度计算值,AmpxiN是N个接收器中的第i个接收器在非有效储层段的斯通利波幅度计算值,N是参与斯通利波幅度计算的接收器总数。
CN201811280231.7A 2018-10-30 2018-10-30 应用声波幅度定量化处理评价火成岩储层有效性的方法 Active CN109374757B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811280231.7A CN109374757B (zh) 2018-10-30 2018-10-30 应用声波幅度定量化处理评价火成岩储层有效性的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811280231.7A CN109374757B (zh) 2018-10-30 2018-10-30 应用声波幅度定量化处理评价火成岩储层有效性的方法

Publications (2)

Publication Number Publication Date
CN109374757A true CN109374757A (zh) 2019-02-22
CN109374757B CN109374757B (zh) 2021-06-29

Family

ID=65390692

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811280231.7A Active CN109374757B (zh) 2018-10-30 2018-10-30 应用声波幅度定量化处理评价火成岩储层有效性的方法

Country Status (1)

Country Link
CN (1) CN109374757B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110173259A (zh) * 2019-05-14 2019-08-27 中石化石油工程技术服务有限公司 一种斯通利波能量补偿校正的方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080173077A1 (en) * 2001-12-04 2008-07-24 Labcyte Inc. Acoustic assessment of characteristics of a fluid relevant to acoustic ejection
CN102200008A (zh) * 2010-03-26 2011-09-28 中国石油天然气股份有限公司 一种基于电成像测井的储层有效性识别方法
RU2433957C1 (ru) * 2010-02-17 2011-11-20 Федеральное государственное образовательное учреждение высшего профессионального образования "Дальневосточный государственный технический рыбохозяйственный университет" Способ получения льдосодержащей суспензии из морской воды и установка для его осуществления
CN102373923A (zh) * 2010-08-20 2012-03-14 中国石油天然气股份有限公司 一种储层识别方法
CN104100262A (zh) * 2014-07-22 2014-10-15 中国海洋石油总公司 一种利用三维定量荧光参数图版识别储层流体性质的方法
CN104314558A (zh) * 2014-08-21 2015-01-28 中国石油集团川庆钻探工程有限公司 用斯通利波能量损失度判别储层有效性的方法
CN104345346A (zh) * 2014-10-30 2015-02-11 中国石油天然气股份有限公司 一种获取裂缝宽度的方法
CN106568846A (zh) * 2016-10-25 2017-04-19 长江大学 岩心裂缝宽度下快慢横波衰减系数测量方法
CN107153096A (zh) * 2017-06-23 2017-09-12 中国石油大学(北京) 页岩储层成缝能力测试方法及装置

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080173077A1 (en) * 2001-12-04 2008-07-24 Labcyte Inc. Acoustic assessment of characteristics of a fluid relevant to acoustic ejection
RU2433957C1 (ru) * 2010-02-17 2011-11-20 Федеральное государственное образовательное учреждение высшего профессионального образования "Дальневосточный государственный технический рыбохозяйственный университет" Способ получения льдосодержащей суспензии из морской воды и установка для его осуществления
CN102200008A (zh) * 2010-03-26 2011-09-28 中国石油天然气股份有限公司 一种基于电成像测井的储层有效性识别方法
CN102373923A (zh) * 2010-08-20 2012-03-14 中国石油天然气股份有限公司 一种储层识别方法
CN104100262A (zh) * 2014-07-22 2014-10-15 中国海洋石油总公司 一种利用三维定量荧光参数图版识别储层流体性质的方法
CN104314558A (zh) * 2014-08-21 2015-01-28 中国石油集团川庆钻探工程有限公司 用斯通利波能量损失度判别储层有效性的方法
CN104345346A (zh) * 2014-10-30 2015-02-11 中国石油天然气股份有限公司 一种获取裂缝宽度的方法
CN106568846A (zh) * 2016-10-25 2017-04-19 长江大学 岩心裂缝宽度下快慢横波衰减系数测量方法
CN107153096A (zh) * 2017-06-23 2017-09-12 中国石油大学(北京) 页岩储层成缝能力测试方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘之的 等: ""火成岩裂缝有效性测井评价"", 《西南石油大学学报(自然科学版)》 *
张承森 等: ""数字声波幅度比值方法在评价碳酸盐岩水平井储层中的应用"", 《长江大学学报(自然版)》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110173259A (zh) * 2019-05-14 2019-08-27 中石化石油工程技术服务有限公司 一种斯通利波能量补偿校正的方法

Also Published As

Publication number Publication date
CN109374757B (zh) 2021-06-29

Similar Documents

Publication Publication Date Title
CN103256046B (zh) 非常规油气藏水平井全缝长压裂参数模拟的方法及装置
CN1217203C (zh) 通过野外测量标定确定应力-速度耦合系数的方法
CN103258091B (zh) 非常规油气藏水平井段三维岩体力学模型建立的方法及装置
CN106556866B (zh) 一种利用岩石物理模板对储层定量解释的方法
CN105134189A (zh) 一种LogGMIR测井地应力识别优质储层的方法
CN106154334A (zh) 基于网格搜索的井下微地震事件实时反演定位方法
CN101231346A (zh) 应用地震波速度预测煤、岩体物理力学参数方法
CN104747183A (zh) 一种碳酸盐岩储层综合分类方法
Gupta et al. Petrophysical characterization of the Woodford shale
CN105629308A (zh) 相控非均质力学参数地应力方法
CN110954944A (zh) 一种断层圈闭含油高度地震预测方法
CN105134185A (zh) 储层流体性质识别方法
CN112177605B (zh) 一种海域天然气水合物有利聚集区主控因素的确定方法
CN112946751B (zh) 一种预测不同储层平面分布的方法及系统
CN112946743B (zh) 区分储层类型的方法
CN112068198B (zh) 基于地震波全波形特征的裂缝破裂尺度的描述方法
CN111381292A (zh) 一种预测砂岩含烃储层的测井解释方法与装置
CN112185469B (zh) 一种预测海域天然气水合物有利聚集区的方法
CN109577969A (zh) 一种基于岩石压缩系数计算碳酸盐岩地层孔隙压力的方法
CN109374757A (zh) 应用声波幅度定量化处理评价火成岩储层有效性的方法
Yu et al. Rock physics modeling of heterogeneous carbonate reservoirs: porosity estimation and hydrocarbon detection
CN115857047B (zh) 一种地震储层综合预测方法
CN112068197B (zh) 基于地震波动力学参数的裂缝破裂尺度的描述方法
CN113376692B (zh) 致密砂岩气水平井压裂改造方案优化方法及装置
CN109339771A (zh) 一种页岩油气层孔隙压力预测方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant