CN109752773B - 用于深部岩体古应力场模拟的非均质力学参数场确定方法 - Google Patents

用于深部岩体古应力场模拟的非均质力学参数场确定方法 Download PDF

Info

Publication number
CN109752773B
CN109752773B CN201811448918.7A CN201811448918A CN109752773B CN 109752773 B CN109752773 B CN 109752773B CN 201811448918 A CN201811448918 A CN 201811448918A CN 109752773 B CN109752773 B CN 109752773B
Authority
CN
China
Prior art keywords
density
rock
fracture
numerical simulation
fracture density
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
CN201811448918.7A
Other languages
English (en)
Other versions
CN109752773A (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 University of Petroleum East China
Original Assignee
China University of Petroleum East China
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 University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201811448918.7A priority Critical patent/CN109752773B/zh
Publication of CN109752773A publication Critical patent/CN109752773A/zh
Application granted granted Critical
Publication of CN109752773B publication Critical patent/CN109752773B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了用于深部岩体古应力场模拟的非均质力学参数场确定方法,包括对矿场区域内的一口钻井取的全直径岩心进行参数测量;选取岩石力学参数,并首次赋值;根据岩石力学参数的值建立全直径岩心的数值模拟模型;计算数值模拟模型上的裂缝密度和全直径岩心上的裂缝密度,当两者接近时,则数值模拟模型成功标定全直径岩心所对应单井的力学参数,当两者不接近时,调整岩石力学参数赋值直至两者接近为止;采用相同的方法对矿场区域内的所有钻井也建立数值模拟模型,最后将所有数值模拟模型的岩石力学参数值作为输入,通过序贯高斯插值的方法,确定矿场区域的非均质力学参数场。该方法适用范围广、简单易行、模拟准确度高更适合于现场运用。

Description

用于深部岩体古应力场模拟的非均质力学参数场确定方法
技术领域
本发明涉及的深部岩体古应力场模拟技术领域,具体涉及一种用于深部岩体古应力场模拟的非均质力学参数场确定方法。
背景技术
研究古构造运动和古应力场对评价裂缝性储层有重要的意义,裂缝和溶洞的充填特征与期次与古应力场有直接的关系,另外,利用古应力场可以有效的预测矿场区块裂缝分布,确定裂缝的延伸方位等,对寻找裂缝性储层以及油田的勘探开发有重要意义。
目前研究古应力场的方法有岩石声发射,声发射技术已经普遍应用于古应力的测试,但是该技术对于浅层岩石是适用的,对于深层或者经受古应力值较高的岩石适用性有限。以及岩石力学实验推断法,就是通过给定与天然条件相似的人工模拟条件,来测定岩石的破裂强度及其物理性质,并记录其破裂时的压力差,最后把该值近似地当作被模拟天然条件下岩石断裂时古构造应力的大小。该方法关键在于人工模拟条件的设定,但时间因素没有模拟成功且进展缓慢。还有数学解析法,根据共轭剪切角大小和断裂带中动力变质岩的剪切角来估算剪应力的大小。该方法简单易行而且可以达到数量级的精度,但是只适用于共轭角小于90°的情况,即在脆性破裂的条件下才成立。此外,还有破裂准则反演、测井资料解释和理论推导等方法研究古应力场。
另外,有限单元法是地应力场数值模拟经常采用的手段,主要针对连续介质进行近似计算的模拟方法,其一般思路:首先离散单个完整的地质体,形成连续网格模型,以节点相连接,并附上相应的岩石力学参数;进一步通过求解位移、应力和应变来研究目的工区的连续场函数;在充分了解边界受力条件和节点平衡条件的基础上,建立相应的方程组(节点位移-未知量,总体刚度矩阵-系数)。采用插值法计算各个节点的位移,从而求解每个单元的应力、应变值。
目前古应力场数值模拟一般仅利用一套平均的地层参数(弹性模量、泊松比、岩石密度等),而且模拟过程中假设目的层岩石是均质的、各向同性的。但实际地层各向异性,具有不同岩性,不同储层序列,所处的古应力场也是有差异的。所以,为了精确的模拟深部岩体古应力场,确定非均质地层力学参数场很有必要。
发明内容
针对上述技术问题,本发明的目的是提供一种适用范围广、简单易行、更适合于现场运用的一种用于深部岩体古应力场模拟的非均质力学参数场确定方法。
为了实现上述目的,本发明采用的技术方案如下:一种用于深部岩体古应力场模拟的非均质力学参数场确定方法,包括如下步骤:
S01: 设在矿场区域内钻井,钻井数量为N口;
S02:令i=1;
S03:对第i口钻井取全直径岩心,并对第j个全直径岩心的基本参数进行测量;
S04:选取建立第k个数值模拟模型所需的岩石力学参数,并对岩石力学参数首次赋值;
S05:根据岩石力学参数的赋值建立第k个数值模拟模型;
S06:统计计算第j个全直径岩心的单井裂缝密度,并将该单井裂缝密度定义为第g个实际裂缝密度;
调整第k个数值模拟模型使其产生裂缝,统计计算第k个数值模拟模型上的裂缝密度,并将该裂缝密度定义为第t个模拟裂缝密度;
S07:将第g个实际裂缝密度与第t个模拟裂缝密度进行拟合,
如果第g个实际裂缝密度与第t个模拟裂缝密度相近,则说明建立第k个数值模拟模型所赋予的岩石力学参数赋值正确,第k个数值模拟模型成功标定第i个单井不同岩性段的力学参数,并执行下一步;
如果第g个实际裂缝密度与第t个模拟裂缝密度不相近,则调整岩石力学参数的赋值,返回S05;
S08:令i=i+1;
S09:如果i>N,则执行下一步,否则返回S03;
S10:确定了N个钻井的岩石力学参数的值之后,将N个钻井的岩石力学参数值作为输入,通过序贯高斯插值的方法,确定矿场区域的非均质力学参数场。
作为优选:所述S03中的基本参数包括岩性、裂缝发育条数、裂缝密度、裂缝产状和裂缝的倾角方位。通过地质描述获得这些参数后,可以对岩心进行定性甚至定量的描述,用于区分不同岩心段。特别是,裂缝发育条数、裂缝密度等裂缝表征参数,可作为标定模拟所需力学参数的约束条件。
作为优选:所述S03中的岩石力学参数包括岩石密度、弹性模量和泊松比。对不同岩性和不同储层序列赋予不同的岩石力学参数,从而得到更加准确的模型。对该模型进行求解得到位移、应力和应变来研究响应的连续场函数,也就是本发明所针对的深部岩体古应力场模拟。
所述S07中所述的第g个实际裂缝密度与第t个模拟裂缝密度相近是指两者的误差不大于0.1.
两者的误差不大于0.1,本发明是确定用于深部岩体古应力场模拟的非均质力学参数场,而且这种确定方法是建立在古应力状态下的裂缝密度基本不变的情况下,以单井裂缝密度为约束,通过多次反演、拟合确定古应力的大小。然而,模拟的裂缝密度与岩心实际的裂缝密度不可能完全一致,因此需要给定一个误差范围,在该误差范围内,认为模拟裂缝密度与实际裂缝密度一致。
相比现有技术,本发明至少具有如下有益效果:
1、通常在应力场模拟过程中,假设目的层岩石是均质的、各向同性的,所以模拟仅仅只使用一套力学参数,但实际地层各向异性,具有不同岩性,不同储层序列,所处的古应力场也是有差异的。因此导致模拟不准确,现实参考性也非常有限。本发明中为了精确的模拟深部岩体古应力场,确定非均质地层力学参数场很有必要,对不同岩性和不同储层序列赋予不同的岩石力学参数,在已标定多口井模拟所需的力学参数的前提下,通过序贯高斯插值等地质统计学方法,插值拟合出整个矿场的力学参数分布,这种非均质力学参数场,模拟深部岩体古应力更加精确。
2、本发明提供的这种用于深部岩体古应力场模拟的非均质力学场确定方法,可以用于矿场尺度裂缝分布预测。
3、本发明提供的用于深部岩体古应力场模拟的非均质力学参数场确定方法,适用范围广、简单易行、模拟准确度高更适合于现场运用。
附图说明
图1为本发明方法的流程图。
图2为钻井数值模拟模型示意图。
图3为实施例1中BC1-12井数值模拟模型与BC1-12井对应的全直径岩心对比图。
图4为实施例1中矿场尺度弹性模量力学参数分布图。
图5为实施例1中不同储层裂缝体密度预测结果,图5a为风化淋滤带裂缝体密度预测结果,图5b为缝洞发育带裂缝体密度预测结果,图5c半充填发育带为裂缝体密度预测结果,图5d为致密发育带裂缝体密度预测结果。
具体实施方式
为了使本领域的技术人员可以更好地理解本发明。
本发明先通过对矿场区域内一个钻井取的全直径岩心进行全面的地质描述和参数测量,设定数值模型所需的力学参数,再建立单井岩心段的数值模拟模型,即根据不同岩性、不同储层序列(风化淋滤带、缝洞发育带、半充填裂缝发育带以及致密带)的差异,分段建立数值模型。统计计算所取的岩心段单井裂缝密度,以单井裂缝密度为约束,反复调整数值模拟模型过程中的力学参数场与实际岩心观察结果对比拟合,定量统计数值模拟模型的裂缝密度,直至数值模拟模型在设定力学参数场下产生的裂缝密度与取心段的裂缝相近时,则数值模拟模型的力学参数场可以看似标定古应力场模拟所需的力学参数。随后,依相同的方法,对矿场区域内的他钻井建立数值模拟模型。最后,根据序贯高斯插值的方法拟合矿场尺度非均质力学参数场分布。
实施例1:为了进一步说明该技术方法的有效性,以非洲乍得H区块Bongor盆地Baobab潜山地区为例,对本发明方法做进一步详细说明,由图2技术路线实施流程图可知,本发明的具体步骤如下:
S01:设在该潜山地区共钻井总数N口;
S02:令i=1;
S03:对第i个钻井即BC1-12井取全直径岩心,并对该全直径岩心的基本参数进行测量;对BC1-12井取的全直径岩心进行全面的地质描述和参数测量;
此处的参数包括岩心取心段、岩性、裂缝发育条数、裂缝密度、裂缝产状、裂缝充填程度、裂缝终止规律、裂缝的类型(成岩缝、构造缝等)和裂缝的倾角方位。
S04:选取建立第k个数值模拟模型所需的岩石力学参数,并对岩石力学参数首次赋值;
对第j个全直径岩心进行地质观察,测量岩心的尺寸等。同时划分不同岩性、不同储层序列并首次假设数值模型所需的力学参数。具体设定的第j个全直径岩心的岩石力学参数如下表1:
表1 岩石力学参数设定
岩性 密度(kg/cm<sup>3</sup>) 弹性模量(Gpa) 泊松比
二长花岗质碎裂岩 2700 36.8 0.21
花岗质碎裂岩 2710 35.2 0.22
二长花岗岩 2720 34 0.23
二长岩 2660 33.8 0.24
混合花岗岩 2840 25.2 0.31
混合花岗岩 2860 21.7 0.33
S05:根据岩石力学参数的赋值建立第k个数值模拟模型;
参见图1,从上到下不同岩性不同储层序列(风化淋滤带、缝洞发育带、半充填裂缝发育带以及致密带)分别分段建立不同的数值模型。
S06:统计计算第j个全直径岩心的单井裂缝密度,并将该单井裂缝密度定义为第g个实际裂缝密度;统计计算全直径岩心的单井裂缝密度方法属于现有技术,此处不做详述;
调整第k个数值模拟模型使其产生裂缝,统计计算第k个数值模拟模型上的裂缝密度,并将该裂缝密度定义为第t个模拟裂缝密度;
S07:将第g个实际裂缝密度与第t个模拟裂缝密度进行拟合,参见图3,第k个数值模拟模型在该力学参数条件下,第k个数值模拟模型的破裂带、裂缝密集带、致密带跟第j个全直径岩心相吻合,证明了古应力场模拟的力学参数的正确性;也即说明建立第k个数值模拟模型所赋予的岩石力学参数赋值正确,第k个数值模拟模型成功标定第i个单井不同岩性段的力学参数,并执行下一步;
如果出现第g个实际裂缝密度与第t个模拟裂缝密度不相近的情况,则需要调整岩石力学参数的赋值,返回S05;如何调整岩石力学参数的赋值,是现有技术,不是本发明的发明点,在此不做详述;
S08:令i=i+1;
S09:如果i>N,则执行下一步,否则返回S03;
S10:根据确定的区域内所有井的岩石力学参数的值,将五个钻井的岩石力学参数值作为输入,通过序贯高斯插值的方法,确定矿场区域的非均质力学参数场。,如图四矿场尺度弹性模量分布图,通过地质统计学插值拟合后,可以得到整个矿场尺区域的弹性模量值,从图中可以看出,不同区域的弹性模量不同,成功的考虑了岩石实际地层的非均质性,确定了该区域的非均质力学参数场,通过确定的非均质力学参数场,可以得出该矿场区域内任何一个位置的钻井时,从该井中取得的全直径岩心的岩石力学参数,具体到实际应用中,可以在非均质力学参数场的基础上,进行矿场尺度的裂缝分布预测。
在该方法反演标定的非均质力学参数场的基础上,利用模拟的古应力场,进行矿场尺度的裂缝分布预测,参见图5a~5d 。风化淋滤带、缝洞发育带、半充填发育带以及致密发育带这四个储层裂缝发育规律一致,表现为:构造高点为裂缝体密度高值区,主要分布在潜山中部以及西北部,潜山两翼为裂缝体密度低值区,主要分布在潜山东北部以及西南部,与应力分布规律一致;风化淋滤带的裂缝体密度远远大于致密带的裂缝体密度。
需要说明的是,以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解;其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (4)

1.一种用于深部岩体古应力场模拟的非均质力学参数场确定方法,其特征在于,包括如下步骤:
S01: 设在矿场区域内钻井,钻井数量为N个;
S02:令i=1;
S03:对第i口钻井取全直径岩心,并对第j个全直径岩心的基本参数进行测量;
S04:选取建立第k个数值模拟模型所需的岩石力学参数,并对岩石力学参数首次赋值;
S05:根据岩石力学参数的赋值建立第k个数值模拟模型;
S06:统计计算第j个全直径岩心的单井裂缝密度,并将该单井裂缝密度定义为第g个实际裂缝密度;
调整第k个数值模拟模型使其产生裂缝,统计计算第k个数值模拟模型上的裂缝密度,并将该裂缝密度定义为第t个模拟裂缝密度;
S07:将第g个实际裂缝密度与第t个模拟裂缝密度进行拟合,
如果第g个实际裂缝密度与第t个模拟裂缝密度相近,则说明建立第k个数值模拟模型所赋予的岩石力学参数赋值正确,第k个数值模拟模型成功标定第i个单井不同岩性段的力学参数,并执行下一步;
如果第g个实际裂缝密度与第t个模拟裂缝密度不相近,则调整岩石力学参数的赋值,返回S05;
S08:令i=i+1;
S09:如果i>N,则执行下一步,否则返回S03;
S10:确定了N个钻井的岩石力学参数的值之后,将N个钻井的岩石力学参数值作为输入,通过序贯高斯插值的方法,确定矿场区域的非均质力学参数场。
2.如权利要求1所述的非均质力学参数场确定方法,其特征在于:所述S03中的基本参数包括岩性、裂缝发育条数、裂缝密度、裂缝产状和裂缝的倾角方位。
3.如权利要求1所述的非均质力学参数场确定方法,其特征在于:所述S03中的岩石力学参数包括岩石密度、弹性模量和泊松比。
4.如权利要求1-3任一项所述的非均质力学参数场确定方法,其特征在于:所述S07中所述的第g个实际裂缝密度与第t个模拟裂缝密度相近是指两者的误差小于等于0.1。
CN201811448918.7A 2018-11-30 2018-11-30 用于深部岩体古应力场模拟的非均质力学参数场确定方法 Active CN109752773B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811448918.7A CN109752773B (zh) 2018-11-30 2018-11-30 用于深部岩体古应力场模拟的非均质力学参数场确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811448918.7A CN109752773B (zh) 2018-11-30 2018-11-30 用于深部岩体古应力场模拟的非均质力学参数场确定方法

Publications (2)

Publication Number Publication Date
CN109752773A CN109752773A (zh) 2019-05-14
CN109752773B true CN109752773B (zh) 2020-11-13

Family

ID=66403391

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811448918.7A Active CN109752773B (zh) 2018-11-30 2018-11-30 用于深部岩体古应力场模拟的非均质力学参数场确定方法

Country Status (1)

Country Link
CN (1) CN109752773B (zh)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080288172A1 (en) * 2000-07-14 2008-11-20 Schlumberger Technology Corporation Simulation method and apparatus for determining subsidence in a reservoir
CN102852516A (zh) * 2012-04-19 2013-01-02 北京大学 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置
CN103278614A (zh) * 2013-04-28 2013-09-04 北京大学 一种动静态岩石力学参数校正方法及装置
CN103913774A (zh) * 2014-04-02 2014-07-09 西南石油大学 基于微地震事件的储层地质力学参数反演方法
US20150161303A1 (en) * 2013-12-11 2015-06-11 Jerome Lewandowski Method for Automated Rock Physics Modeling
CN104730596A (zh) * 2015-01-25 2015-06-24 中国石油大学(华东) 一种基于多尺度因素约束的离散裂缝建模方法
CN104865610A (zh) * 2015-05-26 2015-08-26 中国石油化工股份有限公司胜利油田分公司石油工程技术研究院 一种深部储层岩心的数值表征方法
CN105334536A (zh) * 2015-12-01 2016-02-17 中国石油大学(华东) 致密砂岩储层网状裂缝系统有效性评价方法
CN105629308A (zh) * 2014-11-07 2016-06-01 中国石油化工股份有限公司 相控非均质力学参数地应力方法
CN107784150A (zh) * 2017-01-19 2018-03-09 中国地质大学(北京) 一种裂缝封闭性预测方法
CN107957596A (zh) * 2017-12-20 2018-04-24 刘敬寿 一种基于应力释放率的储层裂缝预测方法
CN108442922A (zh) * 2018-03-26 2018-08-24 刘敬寿 一种水平井最优钻井轨迹预测方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080288172A1 (en) * 2000-07-14 2008-11-20 Schlumberger Technology Corporation Simulation method and apparatus for determining subsidence in a reservoir
CN102852516A (zh) * 2012-04-19 2013-01-02 北京大学 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置
CN103278614A (zh) * 2013-04-28 2013-09-04 北京大学 一种动静态岩石力学参数校正方法及装置
US20150161303A1 (en) * 2013-12-11 2015-06-11 Jerome Lewandowski Method for Automated Rock Physics Modeling
CN103913774A (zh) * 2014-04-02 2014-07-09 西南石油大学 基于微地震事件的储层地质力学参数反演方法
CN105629308A (zh) * 2014-11-07 2016-06-01 中国石油化工股份有限公司 相控非均质力学参数地应力方法
CN104730596A (zh) * 2015-01-25 2015-06-24 中国石油大学(华东) 一种基于多尺度因素约束的离散裂缝建模方法
CN104865610A (zh) * 2015-05-26 2015-08-26 中国石油化工股份有限公司胜利油田分公司石油工程技术研究院 一种深部储层岩心的数值表征方法
CN105334536A (zh) * 2015-12-01 2016-02-17 中国石油大学(华东) 致密砂岩储层网状裂缝系统有效性评价方法
CN107784150A (zh) * 2017-01-19 2018-03-09 中国地质大学(北京) 一种裂缝封闭性预测方法
CN107957596A (zh) * 2017-12-20 2018-04-24 刘敬寿 一种基于应力释放率的储层裂缝预测方法
CN108442922A (zh) * 2018-03-26 2018-08-24 刘敬寿 一种水平井最优钻井轨迹预测方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
Numerical Investigation on the Heat Extraction Capacity of Dual Horizontal Wells in Enhanced Geothermal Systems Based on the 3-D THM Model;Zhixue Sun et al.;《energies》;20180124;第1-19页 *
Research on Stress Sensitivity of Fractured Carbonate Reservoirs Based on CT Technology;Yongfei Yang et al.;《energies》;20171110;第1-15页 *
三维岩石力学参数场随机模拟预测方法及应用;曹正林等;《现代地质》;20091231;第23卷(第6期);第1126-1130页 *
多因素约束下的致密砂岩气藏离散裂缝;冯建伟等;《中国石油大学学报( 自然科学版)》;20160229;第40卷(第1期);第18-26页 *
鄂尔多斯盆地大牛地气田致密砂岩储层三维岩石力学参数场研究;陆诗阔等;《天然气地球科学》;20151031;第26卷(第10期);第1844-1850页 *
非均质储层三维构造应力场模拟方法;朱传华等;《吉林大学学报(地球科学版)》;20160930;第46卷(第5期);第1580-1588页 *
页岩储层构造应力场模拟与裂缝分布预测方法及应用;丁文龙等;《地学前缘(中国地质大学(北京);北京大学)》;20160331;第23卷(第2期);第63-74页 *

Also Published As

Publication number Publication date
CN109752773A (zh) 2019-05-14

Similar Documents

Publication Publication Date Title
CA2820942C (en) A method to improve reservoir simulation and recovery from fractured reservoirs
CN109655903B (zh) 页岩层横波速度预测方法及系统
CN105445791A (zh) 一种基于多种地震属性的地层孔隙压力预测方法
CN113534291A (zh) 岩石力学层约束下的低渗透储层不同尺度裂缝定量预测方法
CN105653815A (zh) 一种基于岩石物理模型理论的定量解释油藏流体分布方法
CN110795878A (zh) 一种隧道涌水量预测方法
CN105386756B (zh) 一种应用应变量计算脆性地层孔隙度的方法
CN105629308A (zh) 相控非均质力学参数地应力方法
CN109236255B (zh) 一种水平井压裂潜力评价方法和装置
CN105093331B (zh) 获取岩石基质体积模量的方法
CN113821956B (zh) 一种深层页岩储层现今地应力结构扰动量的评价方法
CN102288996B (zh) 叠前afi反演定量预测含气概率的方法
CN113158315A (zh) 一种基于静力触探数据的岩土体参数三维非平稳条件随机场建模方法
CN110727031B (zh) 一种基于三维叠前地震反演结果的地应力获取方法
CN107526117B (zh) 基于自动编码和超限学习联合网络的声波速度预测方法
CN104834007A (zh) 地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法
Altun et al. Geostatistical interpolation for modelling SPT data in northern Izmir
CN109752773B (zh) 用于深部岩体古应力场模拟的非均质力学参数场确定方法
CN116306118A (zh) 一种裂缝建模方法
EP0864882A2 (en) Method for estimating or simulating parameters of a stratum structure
CN110954949A (zh) 一种致密砂岩软孔隙度分布反演方法
CN111830593A (zh) 一种储层物性预测方法及装置
CN115144929A (zh) 油气藏储量确定方法、装置、存储介质及电子设备
CN112780253B (zh) 裂缝型储层的预测及评价方法
CN115542399A (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