CN112412445B - 一种确定致密气藏混合井型的数字化表征方法 - Google Patents

一种确定致密气藏混合井型的数字化表征方法 Download PDF

Info

Publication number
CN112412445B
CN112412445B CN202011500091.7A CN202011500091A CN112412445B CN 112412445 B CN112412445 B CN 112412445B CN 202011500091 A CN202011500091 A CN 202011500091A CN 112412445 B CN112412445 B CN 112412445B
Authority
CN
China
Prior art keywords
well
coefficient
type
directional
horizontal
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
Application number
CN202011500091.7A
Other languages
English (en)
Other versions
CN112412445A (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.)
Shaanxi Yanchang Petroleum Group Co Ltd
Original Assignee
Shaanxi Yanchang Petroleum Group 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 Shaanxi Yanchang Petroleum Group Co Ltd filed Critical Shaanxi Yanchang Petroleum Group Co Ltd
Priority to CN202011500091.7A priority Critical patent/CN112412445B/zh
Publication of CN112412445A publication Critical patent/CN112412445A/zh
Application granted granted Critical
Publication of CN112412445B publication Critical patent/CN112412445B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP 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
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • 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
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Geometry (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Mathematical Analysis (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明具体的涉及一种确定致密气藏混合井型的数字化表征方法,该方法以渗透率比值、地层系数和储能系数这3个物性指标作为约束条件,计算计划部署井在各层位的井型部署条件;若各层位只满足定向井部署条件,则第k口井的井型为定向井;若存在任意层位满足水平井部署条件,则第k口井的井型为水平井;计算得到所有井的井型后,若计算结果只包含定向井,则混合井型的类型为直井+定向井;若计算结果只包含水平井,则混合井型的类型为直井+水平井;若计算结果同时包含定向井和水平井,则混合井型的类型为直井+定向井+水平井。本发明提升了井型确定过程的科学性和智能化程度,对致密气藏有效开发具有重要价值。

Description

一种确定致密气藏混合井型的数字化表征方法
技术领域
本发明涉及致密气藏开采技术领域,具体的涉及一种确定致密气藏混合井型的数字化表征方法。
背景技术
鄂尔多斯盆地广泛发育致密气藏,东西南北不同区域发育不同的主力产层,主力产层呈现出砂体纵横向空间展布复杂、多层叠置、跨度大、非均质强等特征。采用单一井型难以形成有效开发,例如:若采用直井分层开发,由于单层产量很低,导致生产效益差;若采用水平井进行整体开采,为了实现对主力产层的充分控制,需要大量的水平井,导致开发成本高。考虑到物性的区域性与层系间的差异性,为实现各主力产层的最大限度动用,通常采用丛式井组进行开发。
丛式井组由直井+定向井/水平井等井型混合构成,直井作为骨架井无需进行井型确定,致密气藏开发井数多,导致井型确定任务繁重;定向井/水平井的选择主要依靠天然气工作者的个人认识,设计结果存在一定差异。如何提升混合井型确定过程的科学性和智能化程度,对致密气藏开发意义重大。
发明内容
本发明旨在针对上述问题,提出一种确定致密气藏混合井型的数字化表征方法。
本发明的技术方案在于:
一种确定致密气藏混合井型的数字化表征方法,以渗透率比值、地层系数和储能系数这3个物性指标作为约束条件,计算计划部署井在各层位的井型部署条件;若各层位只满足定向井部署条件,则第k口井的井型为定向井;若存在任意层位满足水平井部署条件,则第k口井的井型为水平井;计算得到所有井的井型后,若计算结果只包含定向井,则混合井型的类型为直井+定向井;若计算结果只包含水平井,则混合井型的类型为直井+水平井;若计算结果同时包含定向井和水平井,则混合井型的类型为直井+定向井+水平井。
一种确定致密气藏混合井型的数字化表征方法,该方法具体如下:
sgn(x-x0)函数的表达式为:
式中:sgn(x-x0)为符号函数;x为数值变量;x0为任意有效的固定数值;
步骤1:利用公式(1),构建非负函数、负值函数和等于函数:
f(x)=sgn(1+sgnx) (2)
g(x)=1-sgn(1-sgnx)-sgnx (3)
h(x)=sgn(1+sgnx)+sgn(1-sgnx)-1 (4)
式中:f(x)为非负函数,当x≥0时f(x)=1,当x<0时f(x)=0;g(x)为负值函数,当x<0时g(x)=1,当x≥0时g(x)=0;h(x)为等于函数,当x=0时g(x)=1,当x<0或者x>0时g(x)=0。
步骤2:构建物性指标判断系数
式中:cα为渗透率比值判断系数,无量纲;Kik为第k口井在第i层处的平均渗透率,mD;Kjk为第k口井在第j层处的平均渗透率,j∈[1,m],j≠i,mD;Kik/Kjk为渗透率比值,无量纲;α为平均渗透率比值下限,无量纲;cβ为地层系数判断系数,mD·m;hik为第k口井在第i层处的平均有效储层厚度,m;Kikhik为地层系数,mD·m;β为平均地层系数下限,mD·m;cγ为储能系数判断系数,m;为第k口井在第i层处的平均孔隙度,无量纲;Sgik为第k口井在第i层处的平均含气饱和度,无量纲;/>为储能系数,m;γ为储能系数下限,m。
步骤3:根据公式(2)、公式(3)和公式(5),建立第k口井在第i层满足的井型部署条件表达式:
式中:(wt)ik为第k口井在第i层满足的井型部署条件;H为水平井部署条件,其为不等于0的常数;D为定向井部署条件,其为不等于0的常数,且满足D≠H;f(cα)为渗透率比值判断系数的非负函数;f(cβ)为地层系数判断系数的非负函数;f(cγ)为储能系数判断系数的非负函数;g(cα)为渗透率比值判断系数的负值函数;g(cβ)为地层系数判断系数的负值函数;g(cγ)为储能系数判断系数的负值函数。
步骤4:构建单井井型判断系数:
式中:cD为定向井判断系数;cH为水平井判断系数;为各层的井型部署条件之和;根据公式(6)计算得到第k口井在所有层位满足的井型部署条件后,根据公式(4)和公式(7),建立第k口井的井型确定表达式:
(wt)k=wD×h(cD)+wH×h(cH) (8)
式中:(wt)k为第k口井的井型;wD为定向井;wH为水平井;h(cD)为定向井判断系数的等于函数;h(cH)为水平井判断系数的等于函数;
根据公式(8)的计算结果,当(wt)k=wD时代表第k口井的井型为定向井,当(wt)k=wH时代表第k口井的井型为水平井。
一种确定致密气藏混合井型的数字化表征方法,还包括步骤5:定义混合井型判断系数:
式中:cV+D为“直井+定向井”判断系数;cV+H为“直井+水平井”判断系数;cV+D+H为“直井+定向井+水平井”判断系数;为井型之和;nH为水平井的数量;
根据公式(8)计算得到所有井的井型后,根据公式(4)和公式(9),建立确定混合井型的表达式:
wt=(wV+wD)×h(cV+D)+(wV+wH)×h(cV+H)+(wV+wD+wH)×h(cV+D+H) (10)
式中:wt为混合井型,wV为直井;wD为定向井;wH为水平井;wV+wD为“直井+定向井”;wV+wH为“直井+水平井”;wV+wD+wH为“直井+定向井+水平井”;h(cV+D)为“直井+定向井”判断系数的等于函数;h(cV+H)为“直井+水平井”判断系数的等于函数;h(cV+D+H)为“直井+定向井+水平井”判断系数的等于函数;
根据公式(10)的计算结果,当wt=wV+wD时代表混合井型为“直井+定向井”,当wt=wV+wH时代表混合井型为“直井+水平井”,当wt=wV+wD+wH时代表混合井型为“直井+定向井+水平井”。
本发明的技术效果在于:
本发明从渗透率比值、地层系数和储能系数这3个物性指标出发,通过构建非负函数、负值函数、等于函数和判断系数,实现了单井井型和混合井型的数学计算,避免了现有井型确定主要依靠个人经验的缺陷,提升了井型确定过程的科学性和智能化程度,对致密气藏有效开发具有重要价值。
具体实施方式
一种确定致密气藏混合井型的数字化表征方法,具体过程如下:
步骤1:统计目标区域的主力层位数量m,确定计划部署定向井/水平井数量n
目标区域位于鄂尔多斯盆地东南部,为致密气藏发育区,自上而下发育4个主力层位,即m=4,分别是盒8、山1、山2和本溪,目前已在该区域部署了直井作为骨架井,计划采用丛式井组进行开发,其中定向井/水平井数量n为4,总井数为5。
步骤2:统计目标区域的物性参数
表1目标区域物性参数统计表
步骤3:确定目标区域的物性指标下限α、β、γ值、D和H值、wV、wD和wH值为保证水平井的开发效果,适合部署水平井的产层需具备较好的渗流能力和储集能力,物性指标下限过低则水平井开发效果不佳,物性指标下限过高则产层不容易满足,根据致密气藏的水平井开发经验,α通常介于30~300,β通常介于2~100,γ通常介于0.2~1。在本实施例中,取
平均渗透率比值下限α=50;
平均地层系数下限β=3;
储能系数下限γ=0.4;
关于D和H值,根据目标区域定向井和水平井的平均投资进行赋值,公式为:
式中:D为定向井部署条件,其为不等于0的常数,且满足D≠H;H为水平井部署条件,其为不等于0的常数;CD为定向井的平均投资,万元;CH为水平井的平均投资,万元。
本实施例目标区域的CD=900,CH=2300,计算得到D=0.28,H=0.72;
关于wV、wD和wH值,根据目标区域直井、定向井和水平井的平均投资进行赋值,公式为:
式中:wV为直井;wD为定向井;wH为水平井;CV为直井的平均投资,万元;
本实施例目标区域的CV=800,计算得到wV=0.2,wD=0.225,wH=0.575。
步骤4:令k=1。
步骤5:令i=1。
步骤6:计算物性指标判断系数
根据表1,计算得到Kik/Kjk分别等于0.69、0.50和0.91,Kikhik=2.32,第i层渗透率比值的最小值min(Kik/Kjk)=0.50;
利用公式(5)计算得到物性指标判断系数cα、cβ和cγ,得到:cα=-49.5,cβ=-0.6821,cγ=0.0392。
步骤7:利用公式(2)和公式(3)计算当x=cα、x=cβ和x=cγ时的非负函数和负值函数,得到f(cα)=0,f(cβ)=0,f(cγ)=1,g(cα)=1,g(cβ)=1,g(cγ)=0。
步骤8:利用公式(6)计算第k口井在第i层满足的井型部署条件。
步骤9:令i=2到m,重复步骤6到步骤8,计算第k口井在第1至第m层满足的井型部署条件,得到(wt)1k=0.28,(wt)2k=0.28,(wt)3k=0.28,(wt)4k=0.28。
步骤10:利用公式(7)计算单井井型判断系数,得到cD=0,cH=-0.4375。
步骤11:利用公式(4)计算x=cD和x=cH时的等于函数,得到h(cD)=1,h(cH)=0。
步骤12:利用公式(8)计算第k口井的井型(wt)k,得到(wt)k=0.225,此结果表明第k口井的井型为定向井。
步骤13:令k=2到n,重复步骤5到步骤12,计算剩余计划部署井的井型当k=2到n时,计算第k口井的井型,得到:(wt)2=wD=0.225,(wt)3=wD=0.225,(wt)4=wH=0.575。
统计水平井数量,得到nH=1。
步骤14:根据公式(9)计算混合井型判断系数,计算得到cV+D=0.35,cV+H=-1.05,cV+D+H=0。
步骤15:通过公式(4)计算x=cV+D、x=cV+H和x=cV+D+H时的等于函数,得到:h(cV+D)=0,h(cV+H)=0,h(cV+D+H)=1。
步骤16:利用公式(10)计算混合井型的类型,得到wt=1。由于wV+wD+wH=1,故计算结果表明混合井型的类型为“直井+定向井+水平井”。

Claims (2)

1.一种确定致密气藏混合井型的数字化表征方法,其特征在于:以渗透率比值、地层系数和储能系数这3个物性指标作为约束条件,计算计划部署井在各层位的井型部署条件;若各层位只满足定向井部署条件,则第k口井的井型为定向井;若存在任意层位满足水平井部署条件,则第k口井的井型为水平井;计算得到所有井的井型后,若计算结果只包含定向井,则混合井型的类型为直井+定向井;若计算结果只包含水平井,则混合井型的类型为直井+水平井;若计算结果同时包含定向井和水平井,则混合井型的类型为直井+定向井+水平井;
具体为:
步骤1:构建物性指标判断系数
式中:cα为渗透率比值判断系数,无量纲;Kik为第k口井在第i层处的平均渗透率,mD;Kjk为第k口井在第j层处的平均渗透率,j∈[1,m],j≠i,mD;Kik/Kjk为渗透率比值,无量纲;α为平均渗透率比值下限,无量纲;cβ为地层系数判断系数,mD·m;hik为第k口井在第i层处的平均有效储层厚度,m;Kikhik为地层系数,mD·m;β为平均地层系数下限,mD·m;cγ为储能系数判断系数,m;为第k口井在第i层处的平均孔隙度,无量纲;Sgik为第k口井在第i层处的平均含气饱和度,无量纲;/>为储能系数,m;γ为储能系数下限,m;
步骤2:建立第k口井在第i层满足的井型部署条件表达式
式中:(wt)ik为第k口井在第i层满足的井型部署条件;H为水平井部署条件,其为不等于0的常数;D为定向井部署条件,其为不等于0的常数,且满足D≠H;f(cα)为渗透率比值判断系数的非负函数;f(cβ)为地层系数判断系数的非负函数;f(cγ)为储能系数判断系数的非负函数;g(cα)为渗透率比值判断系数的负值函数;g(cβ)为地层系数判断系数的负值函数;g(cγ)为储能系数判断系数的负值函数;
步骤3:构建单井井型判断系数:
式中:cD为定向井判断系数;cH为水平井判断系数;为各层的井型部署条件之和;
根据公式(6)计算得到第k口井在所有层位满足的井型部署条件后,建立第k口井的井型确定表达式:
(wt)k=wD×h(cD)+wH×h(cH) (8)
式中:(wt)k为第k口井的井型;wD为定向井;wH为水平井;h(cD)为定向井判断系数的等于函数;h(cH)为水平井判断系数的等于函数;
根据公式(8)的计算结果,当(wt)k=wD时代表第k口井的井型为定向井,当(wt)k=wH时代表第k口井的井型为水平井。
2.根据权利要求1所述确定致密气藏混合井型的数字化表征方法,其特征在于:还包括定义混合井型判断系数:
式中:cV+D为“直井+定向井”判断系数;cV+H为“直井+水平井”判断系数;cV+D+H为“直井+定向井+水平井”判断系数;为井型之和;nH为水平井的数量;
根据公式(8)计算得到所有井的井型后,根据公式(4)和公式(9),建立确定混合井型的表达式:
wt=(wV+wD)×h(cV+D)+(wV+wH)×h(cV+H)+(wV+wD+wH)×h(cV+D+H) (10)
式中:wt为混合井型,wV为直井;wD为定向井;wH为水平井;wV+wD为“直井+定向井”;wV+wH为“直井+水平井”;wV+wD+wH为“直井+定向井+水平井”;h(cV+D)为“直井+定向井”判断系数的等于函数;h(cV+H)为“直井+水平井”判断系数的等于函数;h(cV+D+H)为“直井+定向井+水平井”判断系数的等于函数;
根据公式(10)的计算结果,当wt=wV+wD时代表混合井型为“直井+定向井”,当wt=wV+wH时代表混合井型为“直井+水平井”,当wt=wV+wD+wH时代表混合井型为“直井+定向井+水平井”。
CN202011500091.7A 2020-12-18 2020-12-18 一种确定致密气藏混合井型的数字化表征方法 Active CN112412445B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011500091.7A CN112412445B (zh) 2020-12-18 2020-12-18 一种确定致密气藏混合井型的数字化表征方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011500091.7A CN112412445B (zh) 2020-12-18 2020-12-18 一种确定致密气藏混合井型的数字化表征方法

Publications (2)

Publication Number Publication Date
CN112412445A CN112412445A (zh) 2021-02-26
CN112412445B true CN112412445B (zh) 2023-11-28

Family

ID=74782394

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011500091.7A Active CN112412445B (zh) 2020-12-18 2020-12-18 一种确定致密气藏混合井型的数字化表征方法

Country Status (1)

Country Link
CN (1) CN112412445B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2243400A (en) * 1995-11-02 2000-06-01 Mcgarry, John E. Producing well artificial lift system control
CN103413030A (zh) * 2013-07-24 2013-11-27 中国石油天然气股份有限公司 一种缝洞型碳酸盐岩气藏动态分析方法及系统
CN104747144A (zh) * 2015-02-06 2015-07-01 中国石油天然气股份有限公司 一种基于天然气多层系分压合求的地质选层方法
CN108343420A (zh) * 2017-12-20 2018-07-31 中国石油天然气股份有限公司 一种多因素协同分析的工厂化作业大井组布井方法
CN110043254A (zh) * 2019-05-11 2019-07-23 中国海洋石油集团有限公司 一种基于电缆地层测试资料地层有效渗透率的获取方法
CN110334431A (zh) * 2019-07-02 2019-10-15 西南石油大学 一种低渗透致密气藏单井控制储量计算及剩余气分析方法
CN111485865A (zh) * 2020-06-04 2020-08-04 陕西延长石油(集团)有限责任公司 一种耦合储能系数和地层系数的气井井距计算方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2013224747B1 (en) * 2013-04-12 2014-07-03 Future Energy Innovations Pty Ltd Detection and Collection System for Fugitive Gases and Effluent Liquids Leaking from Around Drilled Wellheads

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2243400A (en) * 1995-11-02 2000-06-01 Mcgarry, John E. Producing well artificial lift system control
CN103413030A (zh) * 2013-07-24 2013-11-27 中国石油天然气股份有限公司 一种缝洞型碳酸盐岩气藏动态分析方法及系统
CN104747144A (zh) * 2015-02-06 2015-07-01 中国石油天然气股份有限公司 一种基于天然气多层系分压合求的地质选层方法
CN108343420A (zh) * 2017-12-20 2018-07-31 中国石油天然气股份有限公司 一种多因素协同分析的工厂化作业大井组布井方法
CN110043254A (zh) * 2019-05-11 2019-07-23 中国海洋石油集团有限公司 一种基于电缆地层测试资料地层有效渗透率的获取方法
CN110334431A (zh) * 2019-07-02 2019-10-15 西南石油大学 一种低渗透致密气藏单井控制储量计算及剩余气分析方法
CN111485865A (zh) * 2020-06-04 2020-08-04 陕西延长石油(集团)有限责任公司 一种耦合储能系数和地层系数的气井井距计算方法

Also Published As

Publication number Publication date
CN112412445A (zh) 2021-02-26

Similar Documents

Publication Publication Date Title
CN105740574B (zh) 一种古地貌剖面恢复方法
CN103236086B (zh) 一种顾及地表水文上下文的多尺度dem建模方法
CN107622328A (zh) 一种页岩气藏分段多簇压裂水平井产量预测方法
CN110984950B (zh) 一种注采井网井位优化部署的方法
CN105545275B (zh) 气田气井部署方法和装置
CN106285621A (zh) 基于多因素决策海上油田老区新井产能的方法
CN112412445B (zh) 一种确定致密气藏混合井型的数字化表征方法
Shi et al. Urban three-dimensional expansion and its driving forces—A case study of Shanghai, China
CN105092820A (zh) 一种生物碎屑灰岩储层非均质性研究方法
CN114091287B (zh) 基于复杂网络理论评价裂缝连通性及优化裂缝参数的方法
CN107038516A (zh) 一种中渗复杂断块油藏水驱开发效果定量评价方法
CN106503407A (zh) 存在部分连通断层的线性水侵油藏的试井分析方法及装置
CN104574513A (zh) 基于几何形态的准确刻画夹层三维分布的表征方法
CN110991933A (zh) 一种山地页岩气资源评价方法及系统
CN109325251B (zh) 一种gis数据的室内三维自动化建模方法
CN112258603A (zh) 三因素复合影响规律分析的三轴版图绘制方法及其用途
CN108932397A (zh) 一种双重介质油藏数值模拟调参方法
Andreychouk et al. Geomorphology and hydrogeology of an exposed evaporite dome: the Dumre karst area, Central Albania
CN116702261A (zh) 一种光伏建筑群最优化布局方法及系统
CN109086528A (zh) 在civil3d中按排序规则自动命名网状道路的方法
Ahmad Soil classification and micromorphology: A case study of Cholistan Desert
CN107704646A (zh) 一种致密储层体积改造后的建模方法
CN106682327A (zh) 基于bim的大坝实时理论渗压场分布展示方法
CN112561223A (zh) 油气井增产潜力的定量评价方法
CN110968963A (zh) 一种基于电磁理论计算m值的改进三孔隙度模型及方法

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