CN111794738A - 一种基于阵列感应与各向异性联合反演地层倾角的方法 - Google Patents

一种基于阵列感应与各向异性联合反演地层倾角的方法 Download PDF

Info

Publication number
CN111794738A
CN111794738A CN202010647456.2A CN202010647456A CN111794738A CN 111794738 A CN111794738 A CN 111794738A CN 202010647456 A CN202010647456 A CN 202010647456A CN 111794738 A CN111794738 A CN 111794738A
Authority
CN
China
Prior art keywords
anisotropy
conductivity
stratum
array induction
logging
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
CN202010647456.2A
Other languages
English (en)
Other versions
CN111794738B (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 Petroleum Corp
China Petroleum Logging Co Ltd
Original Assignee
China National Petroleum Corp
China Petroleum Logging Co Ltd
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 Petroleum Corp, China Petroleum Logging Co Ltd filed Critical China National Petroleum Corp
Priority to CN202010647456.2A priority Critical patent/CN111794738B/zh
Publication of CN111794738A publication Critical patent/CN111794738A/zh
Application granted granted Critical
Publication of CN111794738B publication Critical patent/CN111794738B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Geology (AREA)
  • Geochemistry & Mineralogy (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Fluid Mechanics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于阵列感应与各向异性联合反演地层倾角的方法,首先基于阵列感应电阻率测井资料,推导出在倾斜非均质地层中双线圈系的视电导率和水平电导率计算公式;引入无限厚的横向各向异性地层中的理想电位电极系视电导率与水平电导率的关系公式;建立目的层段的各向异性参数计算模型;利用阵列感应电阻率测井来处理得到目的层段的视电导率、水平电导率与各向异性建立函数关系式,形成阵列感应与各向异性联合反演地层倾角的方法。在新井中应用本发明的方法,得到目的层段的地层倾角与电成像测井处理得到地层倾角数据相关性好,所以通过阵列感应电阻率测井资料来准确计算致密砂岩地层倾角精度高,为今后寻找致密气藏圈闭提供了可靠资料。

Description

一种基于阵列感应与各向异性联合反演地层倾角的方法
技术领域
本发明属于石油勘探中致密砂岩储层的测井评价技术领域,涉及一种在致密砂岩中利用常规测井和阵列感应测井联合计算出的各向异性来反演出地层倾角的方法。
背景技术
非常规致密砂岩储层的地层倾角参数主要是依靠地层倾角测井仪器直接测量得到和电成像测井精细解释处理得到,很少利用常规基础测井资料联合反演算出。目前,基于电磁感应原理,主要是研究阵列感应受倾角影响特性分析和基于阵列感应测井受倾角影响的校正方法研究,很少利用电导率测井曲线与地层各向异性对地层倾角进行连续定量的反演。然而,地层倾角在致密砂岩气藏中是十分重要的参数,因为通过测量致密砂岩气藏的地层倾角来研究各种地质问题,对指导油气田的勘探开发具有重要意义。为了深入研究致密砂岩的地层倾角特征,来创新了一套基于阵列感应和各向异性反演出有效的地层倾角计算方法,国内外学者未尝试对这套方法进行研究。但事实上,大多数致密砂岩气藏绝对属于非均匀介质地层环境。尽管如此,目前在非均质地层中阵列感应测井因受地层倾角影响导致电流线不再关于仪器轴旋转对称,而是沿空间复杂轨迹流动,使阵列感应合成处理结果异常。经典Maxwell微分方程为:
Figure BDA0002573672690000011
Figure BDA0002573672690000012
式中,E为电场强度,V/m;B为磁感应强度,Wb/m2;H为磁场强度,A/m;D为电位移矢量,C/m2;J为电流密度,A/m2;ω为交变电流的角频率,ω=2πf,f为仪器工作频率,Hz。然而,我们知道现实中的阵列感应测井要从根本上适用于非均质地层储层条件。因此,详细研究阵列感应测井仪器的非均质地层条件下的影响因素,才能解决感应测井在非均质地层中存在的问题。
本发明基于阵列感应电阻率测井资料,推导出在倾斜非均质地层中双线圈系的视电导率和水平电导率计算公式;其次引入无限厚的横向各向异性地层中的理想电位电极系视电导率与水平电导率的关系公式;最后利用阵列感应电阻率测井来处理得到目的层段的视电导率和水平电导率与各向异性建立函数关系式,形成阵列感应与各向异性联合反演地层倾角的方法。致密砂岩地层倾角计算模型是利用研究区阵列感应与各向异性联合反演得到,与电成像测井处理得到地层倾角数据进行对比,验证这套创新方法有效可行,为油田勘探开发中寻找有经济价值的致密气藏圈闭提供了可靠资料。
发明内容
本发明的目的在于提供一种基于阵列感应与各向异性联合反演地层倾角的方法,在倾斜非均质地层中视电导率、水平电导率和各向异性推导的基础上,通过测井资料建立基于阵列感应与各向异性联合反演地层倾角的方法,来精细评价致密砂岩的地层倾角,更为有效地寻找出有经济价值的致密气藏圈闭,在勘探开发前期中起指导作用。
本发明是通过下述技术方案来实现的。
一种基于阵列感应与各向异性联合反演地层倾角的方法,包括以下步骤:
步骤101,首先基于阵列感应电阻率测井资料,推导出在倾斜非均质地层中双线圈系的视电导率和水平电导率计算公式;
步骤102,考虑到地层倾斜对测井电阻率的影响,引入无限厚的横向各向异性地层中的理想电位电极系视电导率与水平电导率的关系公式;
步骤103,由于受构造地层倾角的影响,引用受地层倾角影响较大的测井数据,并建立目的层段的各向异性参数计算公式;
步骤104,利用阵列感应电阻率测井来处理得到目的层段的视电导率和水平电导率与各向异性建立函数关系式,形成阵列感应与各向异性联合反演地层倾角的方法。
所述步骤101,在倾斜非均质地层中,通过推导Maxwell方程组,可以利用阵列感应电阻率测井得到视电导率计算方法。
双线圈系的视电导率模型:
Figure BDA0002573672690000021
水平电导率模型:
Figure BDA0002573672690000022
式中,L为双线圈间距,m;μ为地层的磁导率,H/m;ω为交变电流的角频率,ω=2πf,f是仪器工作频率,HZ;kh为电磁波沿地层水平方向传播的波数,1/m;δ为理想电位电极系视电导率与水平电导率的比值,无因次;i为线圈电流强度的虚数单位,
Figure BDA0002573672690000023
引入无限厚的横向各向异性地层中的理想电位电极系视电导率与水平电导率的关系公式。
视电导率与水平电导率的比值模型:
Figure BDA0002573672690000031
式中:δ为目的层段视电导率与水平电导率的比值,无因次;λ为研究区目的层段的各向异性,无因次;θ为研究区目的层段地层倾角,°。
所述步骤103,利用自然伽玛能谱测井处理得到的泥质含量与研究区的各向异性建立函数关系式如下:
Figure BDA0002573672690000032
式中:λ表示目的层段的各向异性,无因次;Vsh表示研究区致密砂岩中泥质含量,%;a、b为模型的经验系数,无因次。
所述模型公式中的系数a、b利用岩心实验测量得到样品的各向异性,通过自然伽玛能谱测井处理得到的泥质含量来拟合得到。
所述步骤104,建立阵列感应与各向异性联合反演地层倾角模型如下:
Figure BDA0002573672690000033
式中:σa基于阵列感应电阻率测井得到视电导率,s/m;σh沿利用沿地层水平方向传播的电磁波波数反演出水平电导率,s/m。
本发明的有益效果在于,实现了在倾角非均质地层中,利用阵列感应电阻率测井来处理得到目的层段的视电导率、水平电导率与各向异性建立函数关系式,形成阵列感应与各向异性联合反演地层倾角的方法,有效地寻找有经济价值的致密气藏圈闭,指导油田前期勘探开发。
附图说明
图1为本发明提供的一种基于阵列感应与各向异性联合反演地层倾角的方法流程图;
图2为本发明实施例中提供K井区致密砂岩中泥质含量与电阻率各向异性关系图;
图3为本发明实施例中提供K-27井基于阵列感应与各向异性联合反演地层倾角参数的结果验证图;
图4为本发明实施例中提供K-27井基于阵列感应与各向异性联合反演地层倾角与成像测井处理得到地层倾角数据之间对比图。
具体实施方式
下面通过对本实施例中的具体实施情况做进一步详细说明,以支持本发明所要解决的技术问题,但并不作为对发明做任何限制的依据。
如图1所示,本发明实施例提供的一种基于阵列感应与各向异性联合反演地层倾角的方法,包括如下步骤:
步骤101,首先基于阵列感应电阻率测井资料,推导出在倾斜非均质地层中双线圈系的视电导率和水平电导率计算公式;
步骤102,考虑到地层倾斜对测井电阻率的影响,引入无限厚的横向各向异性地层中的理想电位电极系视电导率与水平电导率的关系公式;
步骤103,由于受构造地层倾角的影响,引用受地层倾角影响较大的测井数据,并建立目的层段的各向异性参数计算公式;
步骤104,利用阵列感应电阻率测井来处理得到目的层段的视电导率和水平电导率与各向异性建立函数关系式,形成阵列感应与各向异性联合反演地层倾角的方法。
下面,通过对本实施例的具体实施情况做进一步详细说明,以支持本发明所要解决的技术问题,按照以下步骤进行操作:
步骤一,选取区域内从上到下连续进行过电成像测井测量的井为K-27井,全井段开展电成像测井资料处理并对地层进行勾画,获取不同深度变化下致密砂岩的地层倾角、倾向及方位,为后续的地层倾角参数精细处理准备数据。
步骤二,利用研究区块的K-27井测有阵列感应电阻率测井资料,在倾斜非均质地层中,通过推导Maxwell方程组,可以利用阵列感应电阻率测井得到视电导率计算方法如下:
双线圈系的视电导率模型:
Figure BDA0002573672690000041
水平电导率模型:
Figure BDA0002573672690000042
式中,L为双线圈间距,m;μ为地层的磁导率,H/m;ω为交变电流的角频率,ω=2πf,f是仪器工作频率,HZ;kh为电磁波沿地层水平方向传播的波数,1/m;δ为理想电位电极系视电导率与水平电导率的比值,无因次;i为线圈电流强度的虚数单位,
Figure BDA0002573672690000043
步骤三,利用步骤二中分析得到的视电导率计算公式,引入无限厚的横向各向异性地层中的理想电位电极系视电导率与水平电导率的关系公式。
视电导率与水平电导率的比值模型:
Figure BDA0002573672690000044
式中:δ为目的层段视电导率与水平电导率的比值,无因次;λ为研究区目的层段的各向异性,无因次;θ为研究区目的层段地层倾角,°。
步骤四,再利用研究区块多口井的自然伽玛能谱测井处理得到的泥质含量与研究区的各向异性建立函数关系式如下:
Figure BDA0002573672690000051
式中:λ表示目的层段的各向异性,无因次;Vsh表示研究区致密砂岩中泥质含量,%;a、b为模型的经验系数,无因次。
步骤五,结合阵列感应、自然能谱及电成像测井,利用步骤二中式(1)、(2)确定目的层段的视电导率、水平电导率;在此基础上,结合步骤三中式(3)、(4)来确定视电导率与水平电导率的比值、各向异性与泥质含量关系式,利用步骤五中式(5)来反演出致密砂岩的地层倾角模型:
Figure BDA0002573672690000052
式中:σa基于阵列感应电阻率测井得到视电导率,s/m;σh沿利用沿地层水平方向传播的电磁波波数反演出水平电导率,s/m。
图3为本发明提供的K-27井基于阵列感应与各向异性联合反演地层倾角参数的计算结果,图4为本发明提供的K-27井基于阵列感应与各向异性联合反演地层倾角与成像测井处理得到地层倾角的结果基本一致,验证了这套方法的准确可靠,可以有效识别出致密砂岩储层的地层特征,指导测井资料解释精细评价。
最后应说明的是,以上具体实施方式仅用以说明本发明的技术方案而非限制,尽管参照实例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (6)

1.一种基于阵列感应与各向异性联合反演地层倾角的方法,其特征在于,包括以下步骤:
步骤101,首先基于阵列感应电阻率测井资料,推导出在倾斜非均质地层中双线圈系的视电导率和水平电导率计算公式;
步骤102,考虑到地层倾斜对测井电阻率的影响,引入无限厚的横向各向异性地层中的理想电位电极系视电导率与水平电导率的关系公式;
步骤103,结合研究区的实际情况分析,由于受构造地层倾角的影响,对测井电阻率值造成偏差,因此引用受地层倾角影响较大的测井数据,并建立目的层段的各向异性参数计算模型;
步骤104,利用阵列感应电阻率测井来处理得到目的层段的视电导率、水平电导率与各向异性建立函数关系式,形成阵列感应与各向异性联合反演地层倾角的方法。
2.如权利要求1所述的一种基于阵列感应与各向异性联合反演地层倾角的方法,其特征在于,所述步骤101,在倾斜各向异性地层中,通过推导Maxwell方程组,利用阵列感应测井得到视电导率计算方法;
双线圈系的视电导率模型:
Figure FDA0002573672680000011
水平电导率模型:
Figure FDA0002573672680000012
式中,L为双线圈间距,m;μ为地层的磁导率,H/m;ω为交变电流的角频率,ω=2πf,f是仪器工作频率,Hz;kh为电磁波沿地层水平方向传播的波数,1/m;δ为理想电位电极系视电导率与水平电导率的比值,无因次;i为线圈电流强度的虚数单位,
Figure FDA0002573672680000013
3.如权利要求1所述的一种基于阵列感应与各向异性联合反演地层倾角的方法,其特征在于,引入无限厚的横向各向异性地层中的理想电位电极系视电导率与水平电导率的关系公式;
视电导率与水平电导率的比值模型:
Figure FDA0002573672680000014
式中:δ为目的层段视电导率与水平电导率的比值,无因次;λ为研究区目的层段的各向异性,无因次;θ为研究区目的层段地层倾角,°。
4.如权利要求1所述的一种基于阵列感应与各向异性联合反演地层倾角的方法,其特征在于,所述步骤103,利用自然伽玛能谱测井处理得到的泥质含量与研究区的各向异性建立函数关系式如下:
Figure FDA0002573672680000021
式中:λ表示目的层段的各向异性,无因次;Vsh表示研究区致密砂岩中泥质含量,%;a、b为模型的经验系数,无因次。
5.如权利要求1所述的一种基于阵列感应与各向异性联合反演地层倾角的方法,其特征在于,所述模型公式中的系数a、b利用岩心实验测量得到样品的各向异性,通过自然伽玛能谱测井处理得到的泥质含量来拟合得到。
6.如权利要求1所述的一种基于阵列感应与各向异性联合反演地层倾角的方法,其特征在于,所述步骤104,建立阵列感应与各向异性联合反演地层倾角模型如下:
Figure FDA0002573672680000022
式中:σa基于阵列感应电阻率测井得到视电导率,s/m;σh沿利用沿地层水平方向传播的电磁波波数反演出水平电导率,s/m。
CN202010647456.2A 2020-07-07 2020-07-07 一种基于阵列感应与各向异性联合反演地层倾角的方法 Active CN111794738B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010647456.2A CN111794738B (zh) 2020-07-07 2020-07-07 一种基于阵列感应与各向异性联合反演地层倾角的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010647456.2A CN111794738B (zh) 2020-07-07 2020-07-07 一种基于阵列感应与各向异性联合反演地层倾角的方法

Publications (2)

Publication Number Publication Date
CN111794738A true CN111794738A (zh) 2020-10-20
CN111794738B CN111794738B (zh) 2023-07-25

Family

ID=72811247

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010647456.2A Active CN111794738B (zh) 2020-07-07 2020-07-07 一种基于阵列感应与各向异性联合反演地层倾角的方法

Country Status (1)

Country Link
CN (1) CN111794738B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112610204A (zh) * 2020-12-21 2021-04-06 中国地质大学(北京) 一种地层视电导率测量方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2373565Y (zh) * 1999-03-25 2000-04-12 辽河石油勘探局测井公司 一种不带辅助增益的地层倾角测井仪
CN103291278A (zh) * 2012-03-02 2013-09-11 中国石油集团长城钻探工程有限公司 一种用于径向分量刻度的三维阵列感应刻度装置
CN104122603A (zh) * 2013-04-24 2014-10-29 中国石油集团长城钻探工程有限公司 对三维阵列感应测井仪器中的线圈系进行校准的方法
CN104808248A (zh) * 2015-04-22 2015-07-29 长江大学 基于电成像测井全井眼图像的地层产状检测方法
CN106842313A (zh) * 2015-12-04 2017-06-13 中国石油化工股份有限公司 基于方位叠前地震数据的各向异性参数反演方法
CN107784159A (zh) * 2017-09-19 2018-03-09 中国石油天然气集团公司 一种储层电阻率各向异性系数的确定方法
CN109416411A (zh) * 2016-05-17 2019-03-01 沙特阿拉伯石油公司 使用来自倾斜天线的电磁(em)脉冲确定各向异性和地层倾角
CN111058834A (zh) * 2019-12-06 2020-04-24 中国石油大学(华东) 基于瞬变多分量感应测井的各向异性地层倾角确定方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2373565Y (zh) * 1999-03-25 2000-04-12 辽河石油勘探局测井公司 一种不带辅助增益的地层倾角测井仪
CN103291278A (zh) * 2012-03-02 2013-09-11 中国石油集团长城钻探工程有限公司 一种用于径向分量刻度的三维阵列感应刻度装置
CN104122603A (zh) * 2013-04-24 2014-10-29 中国石油集团长城钻探工程有限公司 对三维阵列感应测井仪器中的线圈系进行校准的方法
CN104808248A (zh) * 2015-04-22 2015-07-29 长江大学 基于电成像测井全井眼图像的地层产状检测方法
CN106842313A (zh) * 2015-12-04 2017-06-13 中国石油化工股份有限公司 基于方位叠前地震数据的各向异性参数反演方法
CN109416411A (zh) * 2016-05-17 2019-03-01 沙特阿拉伯石油公司 使用来自倾斜天线的电磁(em)脉冲确定各向异性和地层倾角
CN107784159A (zh) * 2017-09-19 2018-03-09 中国石油天然气集团公司 一种储层电阻率各向异性系数的确定方法
CN111058834A (zh) * 2019-12-06 2020-04-24 中国石油大学(华东) 基于瞬变多分量感应测井的各向异性地层倾角确定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
徐波: "基于常规电测井资料的各向异性储层水平井测井解释", 《测井技术》 *
徐波: "水平井各向异性储层测井解释方法研究", 《中国优秀博士学位论文全文数据库 工程科技1辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112610204A (zh) * 2020-12-21 2021-04-06 中国地质大学(北京) 一种地层视电导率测量方法
CN112610204B (zh) * 2020-12-21 2023-02-10 中国地质大学(北京) 一种地层视电导率测量方法

Also Published As

Publication number Publication date
CN111794738B (zh) 2023-07-25

Similar Documents

Publication Publication Date Title
US10816494B2 (en) Device for determining petrophysical parameters of an underground formation
RU2382193C2 (ru) Многомерная инверсия данных каротажа с использованием различных масштабов и способ построения изображения глубокозалегающего пласта
CN110685600B (zh) 一种用于地质导向的钻头调整预测方法
US6950749B2 (en) Method for resistivity anisotropy determination in near vertical wells
BR112018072718B1 (pt) Método e aparelho para estimativa de um parâmetro invertido de uma formação subsuperficial
US10732315B2 (en) Real-time inversion of array dielectric downhole measurements with advanced search for initial values to eliminate non-uniqueness
US20170276832A1 (en) Estimating Parameters of Archie's Law and Formation Texture Information
US8005619B2 (en) Method of determining reservoir parameters
CN110927794A (zh) 一种致密储层裂缝识别与孔隙度定量计算的方法
CN112578474A (zh) 一种圈定覆盖区砂岩型铀矿找矿远景区的地物化探组合方法
RU2722861C1 (ru) Способ расчета статических поправок
CN112363226A (zh) 一种非常规油气有利区地球物理预测方法
CN111794738A (zh) 一种基于阵列感应与各向异性联合反演地层倾角的方法
CN106842299A (zh) 一种基于地震属性的裂缝定量化预测的方法
US10114141B2 (en) Subsurface resistivity modeling with electromagnetic fields using a conductive casing
CN116168224A (zh) 基于成像砾石含量的机器学习岩相自动识别方法
CN107939385B (zh) 定量计算极化值及应用的方法
Deng et al. A new index used to characterize the near-wellbore fracture network in naturally fractured gas reservoirs
RU2253886C1 (ru) Способ геофизической разведки для определения нефтепродуктивности трещинных карбонатных коллекторов в трехмерном межскважинном пространстве
US9377554B2 (en) Pore parameters and hydraulic parameters from electric impedance spectra
CN116430464A (zh) 提高时频电磁法频率域资料处理质量的方法
CN108019207A (zh) 一种对称电磁波电阻率的测量方法
RU2314554C1 (ru) Способ размещения наклонных и горизонтальных нефтегазовых скважин на основе спектральной декомпозиции геофизических данных
CN110927819B (zh) 一种裂缝发育程度表征方法
RU2253885C1 (ru) Способ определения нефтепродуктивности пористых коллекторов в трехмерном межскважинном пространстве

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