CN110454153A - 一种核磁共振测井弛豫反演方法 - Google Patents
一种核磁共振测井弛豫反演方法 Download PDFInfo
- Publication number
- CN110454153A CN110454153A CN201910768327.6A CN201910768327A CN110454153A CN 110454153 A CN110454153 A CN 110454153A CN 201910768327 A CN201910768327 A CN 201910768327A CN 110454153 A CN110454153 A CN 110454153A
- Authority
- CN
- China
- Prior art keywords
- echo
- magnetic resonance
- nuclear magnetic
- inverting
- relaxation
- 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 49
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 35
- 238000001228 spectrum Methods 0.000 claims abstract description 15
- 230000002902 bimodal effect Effects 0.000 claims abstract description 6
- 238000012360 testing method Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims 1
- 238000009499 grossing Methods 0.000 abstract description 5
- 238000000354 decomposition reaction Methods 0.000 abstract description 4
- 239000011435 rock Substances 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002969 morbid Effects 0.000 description 1
- 230000005311 nuclear magnetism Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000000053 physical method Methods 0.000 description 1
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/32—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electron or nuclear magnetic resonance
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- 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
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- High Energy & Nuclear Physics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种核磁共振测井弛豫反演方法,包括以下步骤:S1,利用高斯概率分布建立大孔占优双峰模型,设置核磁共振(NMR)采集回波串参数,正演得到衰减的回波串;S2,采用吉洪诺夫(Tikhonov)正则化方法,结合广义交叉校验(GCV)准则,利用最小二乘拟合使得病态方程求解最小化;S3,设置核磁共振回波间隔、横向弛豫时间(T2)谱布点方式,加入不同白噪声,采用基于GCV的Tikhonov方法对回波串进行反演得到T2谱,解决了核磁共振测井采集回波串后反演T2谱的问题,相对传统奇异值分解法或模平滑方法难以选取最佳折衷参数的方法,本方法采取自动选取最优化参数,即使在信噪比100、30、20甚至低至10的情况下,反演的T2谱也能与模型吻合,提高了精度。
Description
技术领域
本发明涉及地球物理测井技术领域,具体涉及一种核磁共振测井弛豫反演方法。
背景技术
核磁共振测井作为测井的一种特殊手段,与电阻率、放射性、声波测井一样也是对地层岩石的一种物理测量,它所采集到的衰减回波串是岩石孔隙中流体弛豫信号的反应。地质油藏需要的孔隙流体属性以及物性属性等都是从采集的回波串中得来的。岩石的核磁共振信息能提供许多流体、物性、岩性信息,获得储层信息的其中一个难点是对岩石回波信号进行反演,反演手段归属的方程是病态方程,也就是说在反拉普拉斯变换中,输入端的细小变化如核磁共振测井采集中的噪音等对反演结果会有很大的影响。所以,必须采用未知解的某种先验信息如数据的非负限制、光滑性、信噪比等来约束,以最佳形式得到岩石中的信息。
在核磁共振反演T2谱领域,国内有很多外学者提出了不同的反演手段,Prammer提出了基于奇异值分解的反演方法,Butler等提出了基于非负约束的模平滑方法,Boqin等提出了基于蒙特卡罗的反演方法。翁爱华提出了高分辨率反演方法,王才志,姜瑞忠等发展了奇异值分解方法,王忠东等提出了基于整体迭代修正的反演方法,谭茂金等提出了基于遗传算法的反演方法。以上反演方法均受到一些方面的制约如最佳平滑因子选择、如何自动选择最佳折衷参数、算法收敛性和速度以及对解估计的评价,并且最难的一点也是大部分算法欠缺的一点就是算法会受低信噪比的严重影响。
发明内容
本发明的目的是在于提供一种核磁共振测井弛豫反演方法,采用吉洪诺夫正则化方法,并且通过引入广义交叉校验准则对解进行评价,在原问题解得基础上加入平滑因子,从而阻尼矩阵高频分量中的噪声,并力求在其中找出一个最佳折衷解。
为解决上述发明的目的,本发明提供技术方案如下:
一种核磁共振测井弛豫反演方法,包括如下步骤:
S1,利用高斯概率分布建立大孔占优双峰模型,设置核磁共振(NMR)T2谱模型及模拟采集回波串参数,正演得到衰减的回波串;
S2,采用吉洪诺夫(Tikhonov)正则化方法,结合广义交叉校验(GCV)准则,利用最小二乘拟合使得病态方程求解最小化;
S3,设置核磁共振回波间隔、横向弛豫时间(T2)谱布点方式,加入不同白噪声,采用基于GCV的Tikhonov方法对回波串进行反演得到T2谱。
进一步地,所述步骤S1中,设置核磁共振T2谱模型及模拟采集回波串参数,具体步骤如下:
通过高斯概率分布建立双孔模型,设置预先要采集的回波串个数,简化的回波串采集模式为:
式中,M(ti)为第i个回波的幅度值(i=1,...,n);
ti为n个同样回波间隔的递减时间;
εi为第n个回波的噪声;
m为事先的横向弛豫时间;
T2为横坐标上的对数分布;
a(Tj)为病态方程解出来的弛豫时间Tj对应的T2幅度。
进一步地,所述步骤S2中,采用吉洪诺夫(Tikhonov)正则化方法,结合广义交叉校验(GCV)准则,具体步骤如下:
将解的极小范数作为附加条件代入到病态方程中,此时相当于原来的问题中加入约束条件下的最小问题,原问题则可转化为:
正则化方法的实质是在解的范数最小和数据的拟合程度之间达到最优折衷解;关于正则参数的选择,为了兼顾算法的收敛与速度,以及收敛问题的真实解,采用基于数据误差后验取法中的广义交叉较验准则,是对算子方程Ax=y的离散形式给出的,
其中A=(aij)m×n;y=(y1,y2,…,ym)T;
这样取a*满足V(α*)=minV(α)此法源于估计理论中选择最佳模型的PRESS准则,但比它更为收敛。
进一步地,所述步骤S3中,对回波串进行反演,具体步骤如下:
采用对数均匀布点,个数为32,分布0.5-2048(ms),回波间隔:TE=0.6ms,之前采集的回波个数为1024个,设置4种不同的信噪比噪声,分别为100、30、20、10,反演的方程离散化后为:
利用吉洪诺夫正则化方法和交叉校验准则参数,结合反演的参数,对采集的回波串进行反演,得到反应储层孔隙结构的T2谱。
本发明相对于现有技术的有益效果是:
本发明提出了一种核磁共振测井弛豫反演方法,采用吉洪诺夫正则化方法,并且通过引入广义交叉校验准则对解进行评价,在原问题解得基础上加入平滑因子,从而阻尼矩阵高频分量中的噪声,并力求在其中找出一个最佳折衷解。解决核磁共振测井采集回波串后反演T2谱的问题,相对传统奇异值分解法或模平滑方法难以选取最佳折衷参数的方法,本方法采取自动选取最优化参数,最终提供的T2谱相比模型,更吻合,精度更高。
附图说明
图1是本发明实施例中建立的核磁共振大孔占优双峰孔隙的模型图;
图2是本发明实施例中采集的回波串图;
图3是本发明实施例中利用吉洪诺夫正则化与广义交叉校验准则的参数选择图。
图4是本发明实施例中信噪比为100、30、20、10时本反演方法在三种不同正则化参数情况下反演结果的对比图。
图5是本发明实施例中一种核磁共振测井弛豫反演方法的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
结合图1、图2、图3、图4、图5所示,本发明一种核磁共振测井弛豫反演方法包括以下步骤:
步骤S1,利用高斯概率分布建立大孔占优双峰模型,设置核磁共振(NMR)采集回波串参数,正演得到衰减的回波串;
具体地,建立大孔占优双峰模型,正演核磁共振(NMR)采集回波串:
通过高斯公式建立双孔模型,如图1中左图所示,然后设置预先要采集的回波串个数,简化的回波串采集模式为:
式中,M(ti)是第i个回波的幅度值(i=1,...,n),ti是n个同样回波间隔的递减时间,εi是第n个回波的噪声,m是事先的横向弛豫时间T2在横坐标上的对数分布,a(Tj)是病态方程解出来的弛豫时间Tj对应的T2幅度,采集的回波串如图2所示。
步骤S2,采用吉洪诺夫(Tikhonov)正则化方法,结合广义交叉校验(GCV)准则,利用最小二乘拟合使得病态方程求解最小化;
具体地,采用吉洪诺夫(Tikhonov)正则化结合广义交叉校验(GCV)准则,利用最小二乘拟合使得病态方程求解最小化:
吉洪诺夫(Tikhonov)正则化方法,最核心的是将解的极小范数作为附加条件代入到病态方程中,此时相当于原来的问题中加入约束条件下的最小问题,原问题则可转化为:
正则化方法的实质是在解的范数最小和数据的拟合程度之间达到最优折衷解;其次关于正则参数的选择,为了兼顾算法的收敛与速度,以及是否收敛问题的真实解,采用基于数据误差后验取法中的广义交叉较验准则,这一准则通常是对算子方程Ax=y的离散形式给出的,其中A=(aij)m×n;y=(y1,y2,…,ym)T
这样可以取a*满足V(α*)=minV(α)此法源于估计理论中选择最佳模型的PRESS准则,但比它更为收敛。通过此方法选择的最优正则化因子如图3所示。
步骤S3,设置核磁共振回波间隔、横向弛豫时间(T2)谱布点方式,加入不同白噪声,采用基于GCV的Tikhonov方法对回波串进行反演得到T2谱。
具体地,设置回波间隔、T2谱布点方式、加入不同噪声,对回波串进行反演得到T2谱:
采用对数均匀布点,个数为32,分布0.5-2048(ms),回波间隔:TE=0.6ms,之前采集的回波个数为1024个,设置4种不同的信噪比噪声,分别为100、30、20、10,反演的方程离散化后为:
利用吉洪诺夫正则化方法和交叉校验准则参数,结合反演的参数,对采集的回波串进行反演,得到反应储层孔隙结构的T2谱,最终得到的反演T2谱如图4所示。
以上所述仅为本发明的实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (4)
1.一种核磁共振测井弛豫反演方法,其特征在于,包括如下步骤:
S1,利用高斯概率分布建立大孔占优双峰模型,设置核磁共振(NMR)T2谱模型及模拟采集回波串参数,正演得到衰减的回波串;
S2,采用吉洪诺夫(Tikhonov)正则化方法,结合广义交叉校验(GCV)准则,利用最小二乘拟合使得病态方程求解最小化;
S3,设置核磁共振回波间隔、横向弛豫时间(T2)谱布点方式,加入不同白噪声,采用基于GCV的Tikhonov方法对回波串进行反演得到T2谱。
2.如权利要求1所述的一种核磁共振测井弛豫反演方法,所述步骤S1中,设置核磁共振T2谱模型及模拟采集回波串参数,具体步骤如下:
通过高斯概率分布建立双孔模型,设置预先要采集的回波串个数,简化的回波串采集模式为:
式中,M(ti)为第i个回波的幅度值(i=1,...,n);
ti为n个同样回波间隔的递减时间;
εi为第n个回波的噪声;
m为事先的横向弛豫时间;
T2为横坐标上的对数分布;
a(Tj)为病态方程解出来的弛豫时间Tj对应的T2幅度。
3.如权利要求2所述的一种核磁共振测井弛豫反演方法,其特征在于,所述步骤S2中,采用吉洪诺夫(Tikhonov)正则化方法,结合广义交叉校验(GCV)准则,具体步骤如下:
将解的极小范数作为附加条件代入到病态方程中,此时相当于原来的问题中加入约束条件下的最小问题,原问题则可转化为:
正则化方法的实质是在解的范数最小和数据的拟合程度之间达到最优折衷解;关于正则参数的选择,为了兼顾算法的收敛与速度,以及收敛问题的真实解,采用基于数据误差后验取法中的广义交叉较验准则,是对算子方程Ax=y的离散形式给出的,
其中A=(aij)m×n;y=(y1,y2,…,ym)T;
4.根据权利要求3所述的一种核磁共振测井弛豫反演方法,其特征在于,所述步骤S3中,对回波串进行反演,具体步骤如下:
采用对数均匀布点,个数为32,分布0.5-2048(ms),回波间隔:TE=0.6ms,之前采集的回波个数为1024个,设置4种不同的信噪比噪声,分别为100、30、20、10,反演的方程离散化后为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910768327.6A CN110454153B (zh) | 2019-08-20 | 2019-08-20 | 一种核磁共振测井弛豫反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910768327.6A CN110454153B (zh) | 2019-08-20 | 2019-08-20 | 一种核磁共振测井弛豫反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110454153A true CN110454153A (zh) | 2019-11-15 |
CN110454153B CN110454153B (zh) | 2022-11-04 |
Family
ID=68487897
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910768327.6A Active CN110454153B (zh) | 2019-08-20 | 2019-08-20 | 一种核磁共振测井弛豫反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110454153B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111898734A (zh) * | 2020-07-10 | 2020-11-06 | 中国科学院精密测量科学与技术创新研究院 | 一种基于mlp的nmr弛豫时间反演方法 |
CN112799141A (zh) * | 2021-01-05 | 2021-05-14 | 中国石油天然气集团有限公司 | 一种快速二维核磁共振测井信号处理及t1t2谱反演方法 |
CN112903567A (zh) * | 2021-04-15 | 2021-06-04 | 中国石油大学(北京) | 孔径测量方法、装置、设备及存储介质 |
CN113743682A (zh) * | 2021-11-03 | 2021-12-03 | 中国科学院精密测量科学与技术创新研究院 | 一种基于有监督深度神经网络的nmr弛豫时间反演方法 |
CN114137623A (zh) * | 2020-09-04 | 2022-03-04 | 中国石油化工股份有限公司 | 核磁测井仪观测模式确定方法、存储介质以及电子设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5363041A (en) * | 1992-12-31 | 1994-11-08 | Schlumberger Technology Corporation | Determining bound and unbound fluid volumes using nuclear magnetic resonance pulse sequences |
US5517115A (en) * | 1993-12-16 | 1996-05-14 | Numar Corporation | Efficient processing of NMR echo trains |
US20080036457A1 (en) * | 2005-03-18 | 2008-02-14 | Baker Hughes Incorporated | NMR Echo Train Compression |
CN106156505A (zh) * | 2016-07-05 | 2016-11-23 | 中国科学技术大学 | 一种基于正交匹配追踪算法的核磁共振t2谱反演方法 |
CN106291723A (zh) * | 2016-07-25 | 2017-01-04 | 中国石油大学(北京) | 基于双参数正则化的核磁共振回波数据反演方法及装置 |
CN109633762A (zh) * | 2019-01-07 | 2019-04-16 | 吉林大学 | 基于正弦函数的相关性约束条件联合反演重磁数据的方法 |
-
2019
- 2019-08-20 CN CN201910768327.6A patent/CN110454153B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5363041A (en) * | 1992-12-31 | 1994-11-08 | Schlumberger Technology Corporation | Determining bound and unbound fluid volumes using nuclear magnetic resonance pulse sequences |
US5517115A (en) * | 1993-12-16 | 1996-05-14 | Numar Corporation | Efficient processing of NMR echo trains |
US20080036457A1 (en) * | 2005-03-18 | 2008-02-14 | Baker Hughes Incorporated | NMR Echo Train Compression |
CN106156505A (zh) * | 2016-07-05 | 2016-11-23 | 中国科学技术大学 | 一种基于正交匹配追踪算法的核磁共振t2谱反演方法 |
CN106291723A (zh) * | 2016-07-25 | 2017-01-04 | 中国石油大学(北京) | 基于双参数正则化的核磁共振回波数据反演方法及装置 |
CN109633762A (zh) * | 2019-01-07 | 2019-04-16 | 吉林大学 | 基于正弦函数的相关性约束条件联合反演重磁数据的方法 |
Non-Patent Citations (2)
Title |
---|
林君等: "基于不等式约束的磁共振信号T_2谱多指数分解法及算法改进", 《吉林大学学报(地球科学版)》 * |
葛成等: "核磁共振T_2谱反演影响因素", 《物探与化探》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111898734A (zh) * | 2020-07-10 | 2020-11-06 | 中国科学院精密测量科学与技术创新研究院 | 一种基于mlp的nmr弛豫时间反演方法 |
CN111898734B (zh) * | 2020-07-10 | 2023-06-23 | 中国科学院精密测量科学与技术创新研究院 | 一种基于mlp的nmr弛豫时间反演方法 |
CN114137623A (zh) * | 2020-09-04 | 2022-03-04 | 中国石油化工股份有限公司 | 核磁测井仪观测模式确定方法、存储介质以及电子设备 |
CN114137623B (zh) * | 2020-09-04 | 2024-04-30 | 中国石油化工股份有限公司 | 核磁测井仪观测模式确定方法、存储介质以及电子设备 |
CN112799141A (zh) * | 2021-01-05 | 2021-05-14 | 中国石油天然气集团有限公司 | 一种快速二维核磁共振测井信号处理及t1t2谱反演方法 |
CN112903567A (zh) * | 2021-04-15 | 2021-06-04 | 中国石油大学(北京) | 孔径测量方法、装置、设备及存储介质 |
CN113743682A (zh) * | 2021-11-03 | 2021-12-03 | 中国科学院精密测量科学与技术创新研究院 | 一种基于有监督深度神经网络的nmr弛豫时间反演方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110454153B (zh) | 2022-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110454153A (zh) | 一种核磁共振测井弛豫反演方法 | |
Ge et al. | Determination of nuclear magnetic resonance T 2 cutoff value based on multifractal theory—An application in sandstone with complex pore structure | |
Li et al. | Coupled time‐lapse full‐waveform inversion for subsurface flow problems using intrusive automatic differentiation | |
AU2014216339B9 (en) | Estimating adsorbed gas volume from NMR and dielectric logs | |
CN104516018B (zh) | 一种地球物理勘探中岩性约束下的孔隙度反演方法 | |
Yu et al. | Porosity estimation in kerogen-bearing shale gas reservoirs | |
US10107930B2 (en) | Hybrid saturation recovery-inversion recovery pulse sequence for improved NMR logging of boreholes | |
CN106855636A (zh) | 基于碳酸盐岩储层露头的原型地质模型地震正演方法 | |
CN107219553B (zh) | 基于gr分频反演的暗河充填预测方法 | |
CN106154323A (zh) | 基于地震拓频处理的相控随机反演薄储层预测方法 | |
CN1458536A (zh) | 由短等待时间测量值增强的cpmg测量值的变换方法 | |
CN106066492A (zh) | 一种核磁共振测井孔隙度回波间隔影响的校正方法 | |
Killingbeck et al. | Multimodal layered transdimensional inversion of seismic dispersion curves with depth constraints | |
CN104727813A (zh) | 一种源储共生型地层的孔隙度测量方法 | |
US10705241B2 (en) | Determining sea water resistivity | |
Grana et al. | Petroelastic and geomechanical classification of lithologic facies in the Marcellus Shale | |
Ding et al. | Reliability analysis of seismic attribute in the detection of fault-karst | |
CN107765307B (zh) | 礁滩型储集体建模方法及系统 | |
Deng et al. | Interpretation of well logs and core data via Bayesian inversion | |
Yang et al. | Controlled-source electromagnetic modeling using a high-order finite-difference time-domain method on a nonuniform grid | |
CN114254505A (zh) | 渗透率模型的构建方法、装置、电子设备及存储介质 | |
CN110320227A (zh) | 一种二维核磁共振d-t2谱反演方法及装置 | |
TAN et al. | A Joint Inversion Method for NMR Dual‐TW Logging Data and Fluid Typing | |
Anand et al. | Toward Accurate Reservoir Characterization from New-Generation NMR Logging | |
Shetty et al. | New multiphysics, multiscale inversion for imaging petrophysical properties in anisotropic, laminated formations |
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 |