CN103487837A - 拟饱含水核磁共振自旋回波信号的分解与合成方法 - Google Patents
拟饱含水核磁共振自旋回波信号的分解与合成方法 Download PDFInfo
- Publication number
- CN103487837A CN103487837A CN201310418890.3A CN201310418890A CN103487837A CN 103487837 A CN103487837 A CN 103487837A CN 201310418890 A CN201310418890 A CN 201310418890A CN 103487837 A CN103487837 A CN 103487837A
- Authority
- CN
- China
- Prior art keywords
- relaxation
- saturation
- signals
- relaxation signals
- full
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
拟饱含水核磁共振自旋回波信号的分解与合成方法。根据已知的饱和度数据和对应的弛豫信号,优化拟合确定不同油样饱和度对弛豫过程的影响系数。对于不同含油饱和度的弛豫信号,用同体积的盐水取代其中的油样,通过引入所述的影响系数,表征不同含油饱和度油样对弛豫信号的影响,据此写出不同含油饱和度下拟饱含水核磁共振弛豫信号。反演拟饱含水岩心的自旋回波弛豫信号,可以在一定程度上消除不同流体类型和不同饱和度的影响,提高核磁共振信号对岩石孔隙结构和流体的识别能力。本发明实现拟饱含水岩心核磁共振自旋回波信号的合成与分解,达到提高核磁共振数据评价岩石孔隙结构与流体识别的精度与可靠性。
Description
技术领域
本发明属于石油勘探技术领域,应用于地层岩石孔隙流体的识别和孔隙结构分析。
技术背景
核磁共振自旋回波信号是核磁共振弛豫谱反演的基础数据,核磁共振弛豫谱可以用于地层岩石孔隙流体的识别和孔隙结构分析。在实际应用中,地层孔隙流体与孔隙结构对自旋回波信号的影响相互叠加,所以,目前直接反演自旋回波信号得到的弛豫谱分布,实际上包含了岩石孔隙流体和孔隙结构相互叠加的信息,导致弛豫谱数据在孔隙流体识别和孔隙结构分析上都出现困难。
现有的自旋回波信号反演方法,都是不加区别的反演包含孔隙流体信息和孔隙结构信息的总回波信号。然后,根据反演弛豫谱的分布形式和数值特征分析孔隙流体类型和孔隙结构特征。尽管通过改进反演算法,可以提高弛豫谱分布的精度和稳定性,但是,由于没有实现孔隙流体信息和孔隙结构信息的分离反演,所以,无法从根本上提高识别孔隙流体和分析孔隙结构的精度。
发明内容
本发明目的在于公开一种核磁共振弛豫信号的分解与合成方法,这种方法可以提高核磁共振数据对岩石孔隙结构和流体类型的识别能力。
概括的说,首先在实验流程中分别测量研究流体样品(油样和盐水)的横向自由弛豫时间,配置不同饱和度的油样,并测量在不同油样饱和度下岩样的横向弛豫时间;实验获得数据后,在数据处理中采用“挖补法”,将不同含油饱和度下岩样的弛豫信号与完全饱含水岩样的弛豫信号进行优化拟合,得到不同饱和度油样对岩石总弛豫信号的影响系数(即流体分离系数),如此,通过分解合成得到不同饱和度下岩石拟完全含水弛豫信号。本发明获得的岩石拟完全含水弛豫信号,在一定程度上消除了不同流体类型和饱和度对总弛豫信号的影响,提高岩石弛豫信号对岩石孔隙结构和流体的识别能力。
本发明通过以下技术方案实现:
一种拟饱含水核磁共振自旋回波信号的分解与合成方法,其特征在于,该方法包括实验流程和数据分析两个部分,其中,
(一)所述实验流程,即进行的核磁实验包括步骤有:
1.1、准备好若干实验岩芯、一定矿化度的盐水和实验油样,并对岩芯进行必要的预处理,完成实验前的准备。(该部分属本领域普通技术人员具备的公知常识以及掌握的实验常规操作,属于现有技术)
1.2、测量或者查询出实验配置盐水和实验油样(比如变压器油,或煤油)的自由横向弛豫时间。
1.3、分不同饱和度阶段,向岩心中注入油样,驱替原先占据孔隙的部分盐水。在各个不同的饱和度阶段,记录下它们饱和度数据,同时测量对应的不同盐水饱和度岩芯的弛豫回波信号。
继而,(二)所述数据分析包括步骤有:
2.1、“挖补法”用以消除油相驰豫信号
使用挖补法,去掉含油的体积,用相同的体积的水替代,将饱含油样和盐水的弛豫信号恢复成拟完全饱含盐水的弛豫信号,即方程(2),
式中,M(t),挖补法恢复的饱含水总信号,也即拟完全含水弛豫信号,无单位;T2w,水的体积弛豫时间,ms;Mo用第一个回波信号代替;在水相和油相信号的前面加上修正系数Vc,方程2变成,
2.2、确定不同流体弛豫信号的分离系数Vc
方程3中的M(t)相当于挖补法确定的岩心饱含水的弛豫信号,设M(t)与实测饱含水岩心弛豫信号Mw(t)之间的均方差为Qa(Vc),
方程中分子部分是Vc的函数。在分子中对Vc求导,并令导数等于零,得到确定Vc的方程(方程5),
2.3、根据不同饱和度So的核磁共振实验数据,计算Vc值,得到含油饱和度(So)与Vc的关系:
2.4、计算岩心拟完全含水的弛豫信号。
2.5、反演岩样拟饱含水弛豫信号,计算岩样的拟饱含水T2谱分布,根据这个T2谱分布计算岩样的孔隙结构参数。
上述技术方案实现了拟饱含水岩心核磁共振自旋回波信号的合成与分解,达到提高核磁共振数据评价岩石孔隙结构与流体识别的精度与可靠性。其中公开的实用化的算法,分离核磁共振自旋回波信号中的包含的孔隙流体信息与孔隙结构信息,分别得到孔隙流体的自旋回波信号与孔隙结构的自旋回波信号,并分别反演孔隙流体自旋回波信号和孔隙结构自旋回波信号,得到孔隙流体的弛豫谱分布和孔隙结构的弛豫谱。在此基础上,分别使用孔隙流体的弛豫谱和孔隙结构的弛豫谱,研究岩石的孔隙流体和孔隙结构,提高对岩石孔隙识别和孔隙结构分析的精度。
附图说明
图1本发明整个方法实施流程示意图
图2为“挖补法”示意图
图3为研究流体的自由横向弛豫时间测量图
图4为Zhen53-240样品饱含盐水的弛豫回波及T2谱反演结果(A:弛豫回波,B:反弹T2谱分布)
图5为Zhen53240盐岩不同含水饱和度的T2谱分布
图6为不同饱和度弛豫回波信号分离后的反演结果对比(Zhen-53-241)
图7为信号分离算法的软件化编程框图
具体实施方式
实施例1(理论依据)
氢核在岩石孔隙中发生的弛豫过程,就是在外加磁场的作用下,氢核与孔隙壁之间相互作用而发生的能量传递和转化的过程。氢核在孔隙中的弛豫过程,与氢核弛豫的边界条件有关,而这个边界条件由孔隙内壁构成,所以,氢核的弛豫过程与岩石的孔隙结构有关,这是核磁共振弛豫信号评价岩石孔隙结构的基础。另一方面,岩石弛豫信号显著受到岩石中所有流体的影响,包括盐水和油气。而能够反映岩石孔隙结构的弛豫信号,主要来源于孔隙壁附近的流体。实际的计算结果表明,当岩石孔隙中存在多种流体时,各种流体的影响叠加在岩石样品的总弛豫信号上,会引起岩石弛豫信号的严重畸变,甚至完全掩盖孔隙结构对弛豫信号的影响。如果不剔除多种流体叠加在弛豫信号上的影响,将会严重降低弛豫信号及相应的反演结果对岩石孔隙结构的评价精度。本次发明对不同含油饱和度下的弛豫信号进行对比分析,得到不同含油饱和度下油样对总弛豫信号的影响程度,确定油样弛豫信号的分离系数,采用“挖补法”,恢复岩样拟完全含水弛豫信号,在一定程度上消除不同流体类型和不同含油饱和度对总弛豫信号的影响,提高弛豫信号对岩样孔隙结构和流体的识别精度。
图1为本发明的实施流程图。
步骤1,分别测量和查询研究流体的自由横向弛豫时间(图2)。对于不同矿化度的盐水,在必要时,需要采用杯状容器,测量盐水在各种矿化度下的自由横向弛豫时间。核磁共振实验仪器采用低场核磁共振分析仪器,例如,牛津公司的MARAN2S核磁分析仪器,采用CPMG(Carr-Purcel1-Meiboom-Gil1)脉冲序列,测量流体的自旋回波信号,并用SVD(Singularvalue Decomposition,奇异值分解)方法反演回波信号,获得研究自由流体样品的T2谱分布。除非明确说明,以下的核磁共振测量方法和数据反演方法与此类似。
本次研究作了煤油和盐水的体积弛豫测量(图2),用几何平均值计算煤油和盐水的横向弛豫时间分别为T2o=875ms,T2w=773ms。
步骤2,对所有的岩芯,完全饱和盐水,测量饱含盐水岩样的弛豫信号,并将弛豫信号反演成T2谱分布。
步骤3,对所有的岩芯,依次分不同饱和度阶段注入油样,记录饱和度数据。在每个饱和度下,测量饱含盐水及油样岩芯的弛豫信号,并将弛豫信号反演成T2谱分布。
步骤4。现有的理论研究已经证明,在饱含油气和盐水的岩石孔隙中,岩石样品的总自旋回波信号由油样信号与盐水信号叠加而成。据此,测量的总自旋回波信号可以表示成(方程(1)),
式中,Mt(t)是含油岩心的测量弛豫总信号相对幅度,无单位;so,含油饱和度,小数;Mo,岩心的初始弛豫强度,无单位;T2o,油的体积横向弛豫时间,ms;Mw(t),水相的弛豫信号,无单位;t,测量时间,ms。
本发明对于饱含油样和盐水岩心的弛豫信号,采用“挖补法”恢复成拟完全含水岩心的弛豫信号(图3)。就是说,挖去油样占据的流体体积,同时,补入同等体积的盐水,同体积油样和盐水对总弛豫信号的影响差异,用分离系数表示。“挖补法”的数学表达式写成方程2。
其中,M(t),“挖补法”恢复的饱含水总信号,也即拟完全含水弛豫信号,无单位;Mt(t)是饱含油样和盐水岩芯实测总弛豫信号,标准化以后无单位;So,含油饱和度,小数;Mo,岩心的初始弛豫强度,无单位;T2o,油的体积横向弛豫时间,ms;T2w,盐水的自由弛豫时间,ms;t,测量时间,ms。Vc,油样弛豫信号分离系数,无单位。
由于Mo在实验中无法测量的,实用中可以用第一个回波信号代替。同时,考虑到原油和盐水在含氢指数上的差异,水相和油相信号的前面加上修正系数Vc。方程2变成,
本次研究作了煤油和盐水的体积弛豫测量,用几何平均值计算煤油和盐水的横向弛豫时间分别为T2o=875ms,T2w=773ms。
步骤5优化对比确定不同饱和度油样的弛豫信号分离系数。同一块岩样完全饱含盐水的弛豫信号,是恢复该块岩样在不同含油饱和度下,拟完全含水弛豫信号的转换依据和标准。换句话说,就是设法最大限度把饱含盐水和油样的芯弛豫信号,使用一种实用化的方法,尽可能恢复成完全含水的弛豫信号。本发明使用分离系数的方法,使方程2恢复的拟完全含水弛豫信号,最大程度接近该块岩芯的实测完全含水弛豫信号,目标函数为方程4。
其中,Qa(Vc)是弛豫信号的优化拟合方差,无单位;方程4中的M(t)就是“挖补法”确定的岩芯拟饱含盐水的弛豫信号;Mw(t),实测饱含盐水岩芯的弛豫信号,无单位;Vc,油样弛豫信号分离系数,无单位。
方程4的分子是Vc的函数,在方程4中,对Vc求导,并令导数等于零,可以得到确定Vc的方程(方程5)。
步骤6,确定分离系数(Vc)实用化计算方法。方程4,方程5解决了在确定饱和度下,计算Vc的方法,但是,并不能实现Vc的实用化计算。既然Vc反映油样对弛豫信号的影响程度,可以推测,Vc应该与含油饱和度So有关。据此,建立饱和度So与Vc之间的拟合公式。
步骤7,在确定Vc以后,根据方程2得到岩样的拟完全含水弛豫信号,并反演这个信号,得到岩样拟完全含水的T2谱分布。拟完全含水岩芯的弛豫信号在一定程度上消除了流体类型和不同油样饱和度对弛豫信号的影响,所以,提高了弛豫信号及相应的T2谱分布对对岩石孔隙结构和流体的刻画能力。
实施例2(具体操作)
首先,式样获得实验数据:
第1步,取3块实验岩芯,分别进行预处理,诸如包括洗油,洗盐,烘干,抽真空等,皆属于本领域常规技术。
第2步,配置具有一定矿化度的盐水,本次使用盐水的矿化度为80000ppm。
第3步,针对配置的盐水和实验油样,分别测量盐水和油样的自由横向弛豫时间。本次实验实施盐水和油样(变压器油)的自由横向弛豫时间分别为773ms(T2w)和875ms(T20)
第4步,用盐水饱和经过预处理的岩心。经过抽真空处理后的岩芯,排除了气泡的影响,岩芯完全被盐水饱和。
第5步,对饱含盐水的岩心测量横向弛豫回波数据,测量参数:Te=0.2ms,0.4ms,Tw=1,3,6,8ms。计算结果表明,Te=0.2ms,Tw=3ms的采集参数组合,可以测量到比较真实有效的弛豫回波数据(图4)。此时,岩芯的含水饱和度为100%。
第6步,分不同饱和度阶段向岩心中注入油样,驱替原先占据孔隙的部分盐水。在每个饱和度阶段,记录饱和度数据(参见表1),同时测量具有不同盐水饱和度岩芯的弛豫回波信号。反演回波信号,得到不同饱和度岩芯的T2谱分布(图5)。
接着,对实验数据处理:
第7步,根据方程2,写出不同饱和度下,实测岩芯弛豫信号恢复到拟饱含水弛豫信号的表达式,即“挖补法”恢复岩芯饱含水弛豫信号表达式。
第8步,将岩心拟饱含水弛豫信号与实测饱含水弛豫信号进行优化拟合,确定不同饱和度油样的信号分离系数(Vc)。
第9步,根据不同油样饱和度弛豫信号确定的分离系数Vc,建立分离系数Vc与油样饱和度S0之间的拟合关系,见方程式6;
第10步,反演岩样拟饱含水弛豫信号,计算岩样的拟饱含水T2谱分布,根据这个T2谱分布计算岩样的孔隙结构参数。
本实施例采用核磁共振实验测量技术,实验仪器采用牛津仪器公司的MARAN Ultra核磁实验仪,数据采集及数据反演均采用该仪器的配套软件。本次发明实施案例,取得了拟完全含水岩心的弛豫信号,对拟完全含水岩芯弛豫信号的反演,并得到T2谱分布。将这个T2谱分布与压汞数据进行方差计算的结果证明,拟完全含水岩芯的弛豫信号改善了对岩石孔隙结构与孔隙流体的识别能力。
表1岩心编号及饱和度数据
本发明的实施案例选用某油田的3块真实岩心(表1)。本发明的技术实施计算结果表明,其中2块岩芯的拟饱含盐水弛豫信号数据改善了对岩石孔隙结构的评价和流体的识别效果(表2,图4)。
表2油水信号分离前后的反演精度分析
拟饱含水弛豫信号的反演分析为:
在Zhen53-240,Zhen53-241岩芯核磁数据中,以压汞数据为标准,对比研究弛豫回波信号在油水信号分离前后的反演精度(表2)。表2的计算结果表明,油水信号分离前后,反演结果对岩石结构评价精度改善幅度在5%-27%之间,平均改善幅度为10.7%(图6)。图6中,纵轴是汞饱和度和T2谱分布的相对值,横轴是毛管力和T2弛豫时间转换的伪毛管力,单位MPa。T2弛豫时间在不同饱和度下转换成伪毛管力时,转换系数不同,这里,都是采用在反演结果与毛管力最大吻合下的转换系数(核磁弛豫时间与伪毛管力之间的转换系数)。
本发明实现了孔隙流体信号与孔隙结构信号的分离。可望在核磁共振岩心分析和核磁共振测井中提高核磁共振数据在识别孔隙流体和评价孔隙结构方面的应用效果。在算法软件化以后,可以提高算法的运用效率。
本发明实施的有效和高效方法是信号分离算法的软件化,图7为信号分离算法的软件化编程框图。
Claims (1)
1.一种拟饱含水核磁共振自旋回波信号的分解与合成方法,其特征在于,该方法包括实验流程和数据分析两个部分,其中,
(一)所述实验流程,即进行的核磁实验包括步骤有:
1.1、准备好若干实验岩芯、一定矿化度的盐水和实验油样,并对岩芯进行预处理;
1.2、测量或者查询出实验配置盐水和实验油样的自由横向弛豫时间;
1.3、分不同饱和度阶段,向岩心中注入油样,驱替原先占据孔隙的部分盐水,在各个不同的饱和度阶段,记录下它们饱和度数据,同时测量对应的不同盐水饱和度岩芯的弛豫回波信号;
继而,(二)所述数据分析包括步骤有:
2.1、“挖补法”用以消除油相驰豫信号
使用挖补法,去掉含油的体积,用相同的体积的水替代,将饱含油样和盐水的弛豫信号恢复成拟完全饱含盐水的弛豫信号,即方程(2),
式中,M(t),挖补法恢复的饱含水总信号,也即拟完全含水弛豫信号,无单位;T2w,水的体积弛豫时间,ms;Mo用第一个回波信号代替;在水相和油相信号的前面加上修正系数Vc,方程2变成,
2.2、确定不同流体弛豫信号的分离系数Vc
方程3中的M(t)相当于挖补法确定的岩心饱含水的弛豫信号,设M(t)与实测饱含水岩心弛豫信号Mw(t)之间的均方差为Qa(Vc),
方程中分子部分是Vc的函数,在分子中对Vc求导,并令导数等于零,得到确定Vc的方程(5)
2.3、根据不同饱和度So的核磁共振实验数据,计算Vc值,得到含油饱和度(So)与Vc的关系:
2.4、计算岩心拟完全含水的弛豫信号;
2.5、反演岩样拟饱含水弛豫信号,计算岩样的拟饱含水T2谱分布,根据这个T2谱分布计算岩样的孔隙结构参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310418890.3A CN103487837B (zh) | 2013-09-13 | 2013-09-13 | 拟饱含水核磁共振自旋回波信号的分解与合成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310418890.3A CN103487837B (zh) | 2013-09-13 | 2013-09-13 | 拟饱含水核磁共振自旋回波信号的分解与合成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103487837A true CN103487837A (zh) | 2014-01-01 |
CN103487837B CN103487837B (zh) | 2016-11-23 |
Family
ID=49828206
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310418890.3A Expired - Fee Related CN103487837B (zh) | 2013-09-13 | 2013-09-13 | 拟饱含水核磁共振自旋回波信号的分解与合成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103487837B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105785457A (zh) * | 2016-05-13 | 2016-07-20 | 吉林大学 | 基于双极性脉冲的地面核磁共振横向弛豫时间的测量方法 |
CN110109037A (zh) * | 2019-06-11 | 2019-08-09 | 东南大学 | 一种面向多弛豫组分的核磁共振信号的反演方法 |
CN110306977A (zh) * | 2018-03-27 | 2019-10-08 | 中国石油化工股份有限公司 | 核磁共振t2谱含烃校正方法及*** |
CN112526622A (zh) * | 2020-12-24 | 2021-03-19 | 中国石油天然气集团有限公司 | 基于成像测井孔隙度谱的拟核磁回波数据计算方法 |
CN112731541A (zh) * | 2020-12-21 | 2021-04-30 | 同济大学 | 核磁回波数据优化反演方法及孔隙结构参数确定方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1193387A (zh) * | 1995-06-26 | 1998-09-16 | 纽马公司 | 采用扩散和弛豫测井测量技术评估地岩的核磁共振***和方法 |
WO2006132861A1 (en) * | 2005-06-03 | 2006-12-14 | Baker Hughes Incorporated | Pore-scale geometric models for interpetation of downhole formation evaluation data |
US20090248311A1 (en) * | 2008-03-26 | 2009-10-01 | Baker Hughes Incorporated | Determination of irreducible water cut-off using two dimensional nuclear magnetic resonance data |
CN102434152A (zh) * | 2011-12-05 | 2012-05-02 | 中国石油天然气股份有限公司 | 一种储层含油饱和度的计算方法 |
-
2013
- 2013-09-13 CN CN201310418890.3A patent/CN103487837B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1193387A (zh) * | 1995-06-26 | 1998-09-16 | 纽马公司 | 采用扩散和弛豫测井测量技术评估地岩的核磁共振***和方法 |
WO2006132861A1 (en) * | 2005-06-03 | 2006-12-14 | Baker Hughes Incorporated | Pore-scale geometric models for interpetation of downhole formation evaluation data |
US20090248311A1 (en) * | 2008-03-26 | 2009-10-01 | Baker Hughes Incorporated | Determination of irreducible water cut-off using two dimensional nuclear magnetic resonance data |
CN102434152A (zh) * | 2011-12-05 | 2012-05-02 | 中国石油天然气股份有限公司 | 一种储层含油饱和度的计算方法 |
Non-Patent Citations (4)
Title |
---|
ELLIOT GRUNEWALD,ET AL.: "A laboratory study of NMR ralzxation times and pore coupling in heterogeneous media", 《GEOPHYSICS》 * |
刘堂晏等: "岩心弛豫信号系数分离方法", 《中国地球物理2009》 * |
徐红军等: "利用核磁共振技术研究岩心的润湿性", 《内蒙古石油化工》 * |
翁爱华等: "油水双相饱和孔隙模型核磁特性理论研究", 《地球物理学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105785457A (zh) * | 2016-05-13 | 2016-07-20 | 吉林大学 | 基于双极性脉冲的地面核磁共振横向弛豫时间的测量方法 |
CN110306977A (zh) * | 2018-03-27 | 2019-10-08 | 中国石油化工股份有限公司 | 核磁共振t2谱含烃校正方法及*** |
CN110109037A (zh) * | 2019-06-11 | 2019-08-09 | 东南大学 | 一种面向多弛豫组分的核磁共振信号的反演方法 |
CN112731541A (zh) * | 2020-12-21 | 2021-04-30 | 同济大学 | 核磁回波数据优化反演方法及孔隙结构参数确定方法 |
CN112731541B (zh) * | 2020-12-21 | 2022-05-03 | 同济大学 | 核磁回波数据优化反演方法及孔隙结构参数确定方法 |
CN112526622A (zh) * | 2020-12-24 | 2021-03-19 | 中国石油天然气集团有限公司 | 基于成像测井孔隙度谱的拟核磁回波数据计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103487837B (zh) | 2016-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102096107B (zh) | 一种根据声波时差和密度反演孔隙扁度进行储层渗透性评价的方法 | |
CN103884633B (zh) | 一种确定岩石渗透率的方法及装置 | |
CN103487837A (zh) | 拟饱含水核磁共振自旋回波信号的分解与合成方法 | |
CN101892837B (zh) | 地层因数确定方法及含油饱和度确定方法 | |
CN109613612A (zh) | 一种碳酸盐岩颗粒滩精细刻画与预测方法 | |
CN104278989B (zh) | 一种获取低孔低渗储层饱和度指数的方法 | |
CN104564042B (zh) | 一种页岩储层脆性的评价方法 | |
CN104819923A (zh) | 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 | |
CN104075974A (zh) | 一种利用低场核磁共振精确测定页岩孔隙度的方法 | |
CN104215652B (zh) | 确定油气饱和度的方法和装置 | |
CN104453874A (zh) | 一种基于核磁共振的砂砾岩储层含油饱和度的计算方法 | |
CN104343445B (zh) | 识别气层与水层、高气油比储层及含气饱和度的方法 | |
CN104863574A (zh) | 一种适用于致密砂岩储层的流体识别方法 | |
CN111897011B (zh) | 储层孔隙特征确定方法、装置及设备 | |
CN103513270A (zh) | 一种基于岩石声学特性的气层识别评价方法及装置 | |
CN103064112A (zh) | 一种检测地层含油气性的吸收方法 | |
CN106154343A (zh) | 计算致密油储层的含油饱和度的方法 | |
CN102913240B (zh) | 一种储层流体识别方法 | |
CN106644879B (zh) | 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置 | |
CN102156297A (zh) | 基于砂岩油藏叠后地震数据的流体替换方法 | |
CN104714249A (zh) | 直接提取流体因子的新方法 | |
CN103604836B (zh) | 一种测定天然气水合物储层饱和度的方法及设备 | |
Meadows | Enhancements to Landro's method for separating time-lapse pressure and saturation changes | |
CN106290103A (zh) | 一种页岩气储层中粘土微孔孔隙度的测定方法 | |
CN111381292B (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161123 Termination date: 20190913 |