CN100349013C - 核磁共振测井t2谱t2截止值的确定方法 - Google Patents
核磁共振测井t2谱t2截止值的确定方法 Download PDFInfo
- Publication number
- CN100349013C CN100349013C CNB2005100722329A CN200510072232A CN100349013C CN 100349013 C CN100349013 C CN 100349013C CN B2005100722329 A CNB2005100722329 A CN B2005100722329A CN 200510072232 A CN200510072232 A CN 200510072232A CN 100349013 C CN100349013 C CN 100349013C
- Authority
- CN
- China
- Prior art keywords
- value
- spectrum
- logging
- area
- match
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 87
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 70
- 230000003595 spectral effect Effects 0.000 title claims description 42
- 238000001228 spectrum Methods 0.000 claims abstract description 101
- 238000009826 distribution Methods 0.000 claims description 39
- 230000002902 bimodal effect Effects 0.000 claims description 36
- 239000012530 fluid Substances 0.000 claims description 27
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 14
- 230000001186 cumulative effect Effects 0.000 claims description 11
- 238000007620 mathematical function Methods 0.000 claims description 3
- 230000000877 morphologic effect Effects 0.000 abstract description 6
- 239000011435 rock Substances 0.000 description 33
- 238000010586 diagram Methods 0.000 description 25
- 230000005311 nuclear magnetism Effects 0.000 description 23
- 238000005259 measurement Methods 0.000 description 21
- 238000002474 experimental method Methods 0.000 description 14
- 230000035699 permeability Effects 0.000 description 10
- 239000011148 porous material Substances 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 8
- 238000013461 design Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 229920006395 saturated elastomer Polymers 0.000 description 4
- 230000010354 integration Effects 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000007635 classification algorithm Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种核磁共振测井T2谱T2截止值的确定方法,包括以下步骤:步骤1、对核磁共振测井T2谱按形态特征进行分类;步骤2、根据不同类型的测井T2谱拟合得出其离心谱;步骤3、将离心谱进行面积积分,得出离心谱的面积;步骤4、计算测井T2谱的面积;步骤5、将测井T2谱的面积进行累加,当累加面积等于或与离心谱面积最相近时,此时对应的T2值即为T2截止值。由此,本发明实现了直接利用核磁共振测井得到的T2谱而确定每个深度点的T2截止值。
Description
技术领域
本发明涉及一种核磁共振测井T2谱T2截止值的确定方法,尤其是一种对核磁测井T2谱采用形态特征分类并逐点确定出每个采样点上测井T2谱的T2截止值的方法。
技术背景
核磁共振测井的原始数据是随时间衰减的自旋-回波串,自旋-回波串包含了储层物性、孔隙类型、孔径大小、流体类型及其分布等十分丰富的信息。自旋-回波串可以表示为多个衰减指数的迭加:
式中:Echo(t)----核磁测井测量得到的随时间t变化的回波幅度;
T2i---第i个组分的横向弛豫时间,i=1,…,N;
Φi---对应于T2i组分的核磁共振孔隙度,i=1,…,N;
通过对核磁测井仪器测量得到的自旋-回波串(即式(1))进行多指数拟合,可以求出每一测量点的T2i--ΦI的分布曲线,即所谓的核磁共振测井T2分布谱。该T2分布谱的横坐标为T2值,纵坐标为信号相对幅度,总的信号幅度与岩石孔隙度有关。在单相流体条件下,T2分布定性反映孔隙大小分布。即大孔径孔隙对应较大的T2值,小孔径孔隙对应较小的T2值。
核磁共振测井是目前唯一一种可以区分孔隙中流体赋存状态的测井方法,因此具有独特的技术优势。通过对核磁测井T2谱进行处理,可以得到岩石孔隙大小及其分布,从而更加直接的确定储层的束缚水饱和度和渗透率等。然而与其它测井新技术(MDT、FMI等)相比,核磁测井技术无论从应用的深度和发挥的地质效果上都还存在很大的差距,其中的一个重要的因素就是现行的T2截止值的确定方法还存在应用上的困难,核磁测井T2截止值的确定技术存在缺陷和局限性。
T2截止值是核磁测井资料分析与资料应用的一个重要的基础参数,它是计算储层可动流体或束缚流体含量的关键参数,也是定量评价储层孔隙结构、渗透率等的基础参数。所谓T2截止值,是指T2谱上存在一个确定界限值,大于这个数值所对应的孔隙中的流体称之为自由流体,在储层条件下是可以流动的。小于这个数值对应的孔隙流体称之为(毛管)束缚流体,是不能流动的。即
式中,
T2cutoff——T2截止值;
BVI——毛管束缚水孔隙体积;
FFI——自由流体孔隙体积(或称自由流体指数)。
T2截止值作为表征储层毛管孔隙中流体特性的一个基础参数,可以通过实验室测量确定。国外的研究人员通过大量岩心(主要来自北美地区)实验测量,认为碎屑岩储层岩石的T2截止值为33ms,碳酸盐岩储层岩石为92ms(Morriss等,1993),即T2截止值在已知岩性的情况下基本上是一个常数数值。事实上,目前国内外基本上也是采用这两个截止值对核磁测井资料进行处理和应用的。
自从上世纪90年代中期核磁测井技术引入到中国以来,关于岩心T2截止值的核磁共振实验室测量结果却不支持国外同行普遍认同的T2截止值是某个常数数值的结论。特别是在碎屑岩储层中,T2截止值的变化范围相当大,根本不是一个定值。T2截止值随地区或储层不同而差异很大。图1为我国部分油田砂岩储层根据岩心核磁共振实验得到的T2截止值数据汇总,从图1中可以明显看出T2截止值并不是一个固定的值,也看不出数据有趋向于33ms(碎屑岩地层)的迹象,对国外适用的33ms(碎屑岩地层)在我国的实际地层中变得不适用了。
目前在国内,大家基本上认可了砂岩碎屑岩储层的T2截止值不是33ms这样一个定值的客观事实,但是关于砂岩储层的T2截止值为什么会发生如此大的变化、是什么因素控制了砂岩储层的T2截止值的变化规律仍没有明确的认识。因此导致了核磁测井资料处理和应用的一个难题和尴尬局面:尽管已经认识到砂岩储层的T2截止值不是一个定值,但由于无法找到一个确定T2截止值的适用方法,在核磁测井的资料处理时,仍然无奈地沿用国外的做法,即仍然使用33ms作为砂岩储层统一的T2截止值。这样做的一个严重后果就是,技术人员对核磁测井与T2截止值有关的处理结果持怀疑态度,因此制约了核磁测井技术的应用。
目前已发表的文献中有关求取T2截止值的方法有以下几种:
①.岩心实验室核磁共振实验法(谱面积法)。如图2所示,先对饱和水的岩心进行核磁共振测量,得到横向弛豫时间的饱和T2分布谱;然后在一定的压力条件下将饱和水的岩心在离心机上甩出其中的自由水,只剩下束缚水后,再对岩心进行核磁共振测量得到离心后的T2分布谱,两者相比较。在测井T2谱中从低端(小的T2值)开始进行面积积分S1,当积分面积S1=S2(离心谱总面积)时,测井T2谱中对应T2值即为T2截止值(记为T2cutoff)。这种方法是根据T2截止值的定义所确定的,也是准确的确定T2截止值的方法。离心谱面积与测井T2谱面积的比值为束缚水饱和度(用符号Swi表示)。
②.毛管压力曲线刻度法。首先测量岩心的孔隙度和毛管压力曲线,由毛管压力曲线确定近似于油藏的束缚水饱和度Swir,而岩心的束缚流体体积为φSwir,然后按图3的方法刻度T2截止值。按图4所示的迭代图解法优化T2截止值。图4中,横坐标为油层测井深度范围岩心刻度T2截止值,纵坐标为同一油层深度范围的T2谱和用常规岩心NMR分析确定的T2截止值所求的BVI之间的累计误差。最小累计误差所对应的T2截止值为地层最佳T2截止值。进行岩样刻度时,应该收集地区有代表性取心资料,并将取心资料进行深度归位处理。因取心资料在深度上可能存在误差,所以应与经过校正的测井资料进行仔细对比。
③.渗透率交会图法。在遇到取心困难,取心和岩心分析时间长以及不易获得有代表性的岩心样品时,可利用渗透率交会图技术确定地层的T2截止值。由NMR测井深度范围内的测井曲线获得一系列T2截止值。通过改变T2截止值并计算其相应的Timur-Coates渗透率值(KNMR),然后与同一深度范围MDT(模块式地层测试仪)压降测量得出的渗透率值(KMDT)绘制交会图,选择交会图上线性相关性最好的NMR测井渗透率相应的T2截止值为最佳T2截止值。
④.其他数学统计方法。参见2004年2月《石油地球物理勘探》杂志中,高楚桥,何宗斌等著作的《核磁共振T2截止值与毛管压力的关系》的文章,其中记载了T2截止值与毛管压力Pc之间存在指数关系,得出T2c=A/PcB。参见2004年9月《地球物理学进展》杂志中,汪中浩,章成广,肖承文等著作的《低渗透储层T2截止值实验研究》的文章,其中记载了T2截止值与岩石物性参数(
)之间存在一定的关系,得出
上一部分所论述的确定T2截止值各种方法具有各自的优势和局限性。前两种方法的优点是可以获得准确的T2截止值。但其应用上得局限性是需要岩心样品,都需要岩心实验分析或要用到其它的辅助测井方法,成本高、周期长。由于层段的限制,当有限的岩心实验数据不存在良好的规律(如一批实验数据变化规律很复杂或根本没有什么变化规律可寻),即使得到了这些实验数据也无法应用到测井资料的连续处理。方法3和方法4都是基于地区特征得出的经验公式,难以应用于其他地区,而且公式中的参数也不是唯一的,视地层情况不同而变化。总之,现行确定T2截止值各种方法的一个明显的局限性就是需要岩心实验或其他来源的信息进行标定。况且不论很多时候根本找不出某井或某储层T2截止值实验或标定结果的统计规律,即使能做到这点,也是地区规律并且往往存在很大误差,无法普遍推广与使用,而且其可靠程度与精度也得不到保证。
为了解决核磁测井应用中T2截止值准确确定的实际难题,应该采用全新的研究思路与方法,即避开岩心实验或其他信息标定的传统思路与方法,研究直接根据核磁共振测井T2谱(饱和T2谱)的形态特征确定每个测量点的T2截止值。提供确定核磁共振测井T2截止值的一种可靠方法和技术,克服现行方法的技术局限性。推动核磁测井技术的应用,更大程度地发挥其技术优势,使这项技术在油田勘探开发中发挥重要的作用。通过研究,我们认为在只有核磁测井T2谱的情况下,是能够根据测井T2谱的形态来确定每个深度点的T2截止值的。
发明内容
本发明的目的是直接利用核磁共振测井得到的T2谱而确定每个深度点的T2截止值。
为实现上述目的,本发明提供了一种核磁共振测井T2谱T2截止值的确定方法,包括以下步骤:
步骤1、对核磁共振测井T2谱按形态特征进行分类;
步骤2、根据不同类型的测井T2谱拟合得出其离心谱;
步骤3、将离心谱进行面积积分,得出离心谱的面积;
步骤4、计算测井T2谱的面积;
步骤5、将测井T2谱的面积进行累加,当累加面积等于或与离心谱面积最相近时,此时对应的T2值即为T2截止值。
所述步骤5之后还具有步骤6、将离心谱的面积除以测井T2谱的面积,即为束缚水饱和度的值。所述核磁共振井T2谱分为双峰型、单峰型和三峰型。所述双峰型分为双峰小孔型、双峰大孔型,单峰型分为单峰小孔型、单峰大孔型。所述步骤2中单峰大孔型测井T2谱T2截止值的计算方法为利用数学函数对大孔段进行拟合,其数学形式为:
其中:
A-------拟合T2谱中的幅值极大值,A值取T2分布谱中的幅值极大值;
B-------某个T2时间值,B值取幅度值减小到某个极小时的T2值;
C-------某个待定的参数;
x-----各组份的T2时间值与测井T2谱的T2时间值相同;
f(x)-----拟合得到的幅度值;
所述步骤2中单峰小孔型测井T2谱T2截止值的计算方法为,假设孔隙中的流体都是束缚流体,将孔隙赋予一个固定的T2截止值。所述固定的T2截止值为33ms。
所述步骤2中双峰型测井T2谱T2截止值的计算方法为,利用正态函数在测井T2谱上拟合出一个离心谱,其数学形式为:
其中:
A为拟合T2谱幅度值,A值取T2分布谱中第一个峰的幅值,即使拟合的函数图形和T2分布谱有相同的幅度;
μ为拟合T2谱的期望值,μ值取第一个峰值所对应的T2时间值;
σ为拟合T2谱的方差;
x为各组份的T2时间值与测井谱的T2时间值相同;
f(x)为拟合得到的幅度值;
所述步骤2中三峰型测井T2谱T2截止值的计算方法为利用正态函数在测井T2谱上拟合出一个离心谱,其数学形式为:
其中:
A为拟合T2谱幅度值,A值取T2分布谱中第一个峰的幅值,即使拟合的函数图形和T2分布谱有相同的幅度;
μ为拟合T2谱的期望值,μ值取第一个峰值所对应的T2时间值;
σ为拟合T2谱的方差;
x为各组份的T2时间值与测井谱的T2时间值相同;
f(x)为拟合得到的幅度值;
所述面积的计算方法为谱龙贝格求积法进行面积积分。
所述步骤5中T2截止值的计算方法为将测井T2谱的面积进行累加,当累加面积等于或与离心谱面积最相近时,此时对应的T2值即为T2截止值。上述各种类型均采用此种方法计算T2截止值。
因此,本发明实现了直接利用核磁共振测井得到的T2谱而确定每个深度点的T2截止值。
下面结合附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为新疆、华北以及西北部分地区的岩心实验的T2截止值统计图。
图2为离心法确定T2截止值示意图。
图3为毛管压力曲线刻度法确定T2截止值示意图。
图4为BVI累计误差与T2截止值关系曲线图。
图5为本发明核磁共振测井T2谱T2截止值的确定方法双峰小孔型岩心核磁共振测量结果示意图。
图6为本发明核磁共振测井T2谱T2截止值的确定方法双峰小孔型实际核磁测井仪器测量结果示意图。
图7为本发明核磁共振测井T2谱T2截止值的确定方法双峰大孔型岩心核磁共振测量结果示意图。
图8为本发明核磁共振测井T2谱T2截止值的确定方法双峰大孔型实际核磁测井仪器测量结果示意图。
图9为本发明核磁共振测井T2谱T2截止值的确定方法单峰大孔型岩心核磁共振测量结果示意图。
图10为本发明核磁共振测井T2谱T2截止值的确定方法单峰大孔型实际核磁测井仪器测量结果示意图。
图11为本发明核磁共振测井T2谱T2截止值的确定方法单峰小孔型岩心核磁共振测量结果示意图。
图12为本发明核磁共振测井T2谱T2截止值的确定方法单峰小孔型实际核磁测井仪器测量结果示意图。
图13为本发明核磁共振测井T2谱T2截止值的确定方法三峰型岩心核磁共振测量结果示意图。
图14为本发明核磁共振测井T2谱T2截止值的确定方法三峰型实际核磁测井仪器测量结果示意图。
图15为本发明核磁共振测井T2谱T2截止值的确定方法实施例双峰型示意图1。
图16为本发明核磁共振测井T2谱T2截止值的确定方法实施例双峰型示意图2。
图17为本发明核磁共振测井T2谱T2截止值的确定方法实施例单峰型示意图。
图18为本发明核磁共振测井T2谱T2截止值的确定方法实施例三峰型示意图。
图19为本发明核磁共振测井T2谱T2截止值的确定方法新疆油田某井实际处理的效果图。
图20为本发明核磁共振测井T2谱T2截止值的确定方法西北地区某井实际处理的效果图。
图21为本发明核磁共振测井T2谱T2截止值的确定方法新疆等四个地区计算的T2截止值与岩心分析的T2截止值结果对比图。
图22为本发明核磁共振测井T2谱T2截止值的确定方法新疆等四个地区计算的束缚水饱和度与岩心分析的束缚水饱和度结果的对比图。
图23为本发明核磁共振测井T2谱T2截止值的确定方法计算机计算流程示意图。
图24为本发明核磁共振测井T2谱T2截止值的确定方法流程示意图。
具体实施方式
通过我们多年的研究,在大量实验数据和实测资料的基础上,提出了一种基于核磁测井T2谱形态特征来确定T2截止值的方法,简称形态法。形态法确定T2截止值的流程示意图如图24所示:
第一步,对核磁共振测井的T2谱按形态特征进行分类。
由于不同孔隙结构的岩心饱和谱与离心谱特征差异较大,经过对比分析,将T2饱和谱分成以下三大类五小类,分别为:a.双峰型(又分为双峰大孔型和双峰小孔型);b.单峰型(又分为单峰大孔型和单峰小孔型);c.三峰型。对应不同的孔隙结构类型,计算T2截止值的方法也有差别。对于低孔低渗储层岩心,T2分布谱形态主要为单峰小孔型、双峰小孔型、双峰大孔型,对于高孔、中一高渗储层岩石T2分布谱则通常以单峰大孔型为主。
各类谱的特征和形态描述如下:
a.双峰型(亚类:双峰小孔型、双峰大孔型)
双峰小孔型:如图5和图6所示,为双峰小孔型岩心核磁共振测量示意图和实际核磁测井仪器测量示意图。T2分布谱呈双峰分布,在整个T2分布谱上有两个幅值极大值,且第一个峰的幅值要大于第二个峰的幅值,由于小孔隙所对应的T2时间(第一个峰值点的T2值)要小于大孔隙的T2时间(第二个峰值点的T2值),所以该类型被称为双峰小孔型。
双峰大孔型:如图7和图8所示,为双峰大孔型岩心核磁共振测量示意图和实际核磁测井仪器测量示意图。T2分布谱呈双峰分布,在整个T2分布谱上有两个幅值极大值,且第一个峰的幅值要小于第二个峰的幅值,由于小孔隙所对应的T2时间(第一个峰值点的T2值)要小于大孔隙的T2时间(第二个峰值点的T2值),所以该类型被称为双峰大孔型。
b.单峰型(亚类:单峰大孔型、单峰小孔型)
单峰大孔型:如图9和图10所示,为单峰大孔型岩心核磁共振测量示意图和实际核磁测井仪器测量示意图。T2分布谱呈单峰分布,在整个T2分布谱上只有一个幅值极大值,且该极大值所对应的T2时间比较靠后(一般大于100ms),也就是在大孔隙段的某个T2值达到极大,所以称为单峰大孔型。
单峰小孔型:如图11和图12所示,为单峰小孔型岩心核磁共振测量示意图和实际核磁测井仪器测量示意图。T2分布谱呈单峰分布,在整个T2分布谱上只有一个幅值极大值,且该极大值所对应的T2时间比较靠前(一般大于100ms),也就是在小孔隙段的某个T2值达到极大,所以称为单峰小孔型。
C.三峰型如图13和图14所示,为三峰型岩心核磁共振测量示意图和实际核磁测井仪器测量示意图。
T2分布谱呈三峰分布,在整个T2分布谱上有三个幅值极大值,且最大的峰值不定,可能是第一个峰,可能是第二个峰,也可能是第三个峰。
第二步,针对各类T2分布谱,采用不同的方法确定其T2截止值,具体确定方法如下:
a.双峰型的T2截止值的确定方法:
利用正态函数在测井T2谱上拟合出一个离心谱,其数学形式为:
式中:A-----拟合T2谱幅度值;μ-----拟合T2谱的期望值;σ-----拟合T2谱的方差;x-----各组份的T2时间值(和测井T2谱的T2时间值相同);f(x)-----拟合得到的幅度值。
对于要准确确定T2截止值,则上述A、μ、σ三个参数的确定至关重要,我们的确定方法如下:
A值取T2分布谱中第一个峰的幅值,即使拟合的函数图形和T2分布谱有相同的幅度;μ值取第一个峰值所对应的T2时间值;σ值和以第一个峰值的一半所对应的T2分布谱上两个T2值的宽度有一定的相关性。
b.单峰大孔型的T2截止值的确定方法:
利用数学函数对大孔段进行拟合,其数学形式为:
式中:A-------拟合T2谱中的幅值极大值;B-------某个T2时间值;C-------某个待定的参数;x-----各组份的T2时间值(和测井T2谱的T2时间值相同);f(x)-----拟合得到的幅度值。
如前所述,上述A、B、C三个参数的确定对确定T2截止值也是至关重要。研究发现,A值取T2分布谱中的幅值极大值,B值取幅度值减小到某个极小时的T2值,C值也和以第一个峰值的一半所对应的T2分布谱上两个T2值的宽度有一定的相关性。
c.单峰小孔型的T2截止值的确定方法:
单峰小孔型所对应的孔隙结构多为微毛细管孔隙或互不连通的死孔隙,流体被完全封闭在其中,不能流动,这样的孔隙结构中的流体是不能被采出的或采出量非常有限,在实际生产中也没有工业价值。因此,我们可以给这类孔隙赋予一个固定的T2截止值例如33ms,认为这种孔隙中的流体都是束缚流体。
d.三峰型的T2截止值的确定方法
三峰型的T2截止值的确定采用和双峰型相同的正态函数来求取,其方法如双峰型的求取方法。
有了以上各分类的算法后,用面积积分法来计算T2截止值以及束缚水饱和度:
①.按上述方法拟合得出离心谱;
②.将离心谱按龙贝格(Romberg)求积法进行面积积分,得出的值相当于离心谱的面积;
③.用同样的算法计算测井T2谱的面积;
④.将测井T2谱的面积进行累加,当累加面积等于或与离心谱面积最相近时,此时对应的T2值即为T2截止值;
⑤.将离心谱的面积除以测井T2谱的面积,得到的值即为束缚水饱和度。
按照以上的思路,我们用Visual Fortran编写了该方法的程序,实现了自动处理。流程图如图23所示。
利用形态法我们对新疆等地区的部分岩心做了计算,各分类的拟合效果图如下:
双峰型:如图15和16所示,对小孔段直接拟合,积分求出其面积即为束缚流体的面积,利用该面积除以总谱面积得到的是束缚水饱和度;对总谱从小T2值开始做面积积分,当累计面积和束缚流体面积相等时,此时的T2值即为T2截止值。
单峰型:如图17所示,对大孔段拟合后,积分求出其面积,然后用总谱面积减去大孔段面积,得到的即为束缚流体的面积,利用该面积除以总谱面积得到的是束缚水饱和度;对总谱从小T2值开始做面积积分,当累计面积和束缚流体面积相等时,此时的T2值即为T2截止值。
三峰型:如图18所示,对小孔段直接拟合,积分求出其面积即为束缚流体的面积,利用该面积除以总谱面积得到的是束缚水饱和度;对总谱从小T2值开始做面积积分,当累计面积和束缚流体面积相等时,此时的T2值即为T2截止值。
利用形态法自动确定核磁共振测井T2截止值方法以及自主开发的软件包,我们对新疆油田某井和西北地区某井进行了处理,如图19和20所示为新疆油田某井实际处理的效果图和西北地区某井实际处理的效果图,处理的效果图如下:
图中第一道为CALI(井径曲线)和GR(自然伽马曲线);第二道为深度;第三道中SWICAL为利用形态法计算得到的束缚水饱和度曲线,SWI 33为以33ms作为T2截止值时计算得到的束缚水饱和度曲线,CSWI为岩心分析得到的束缚水饱和度;第四道中T2_AMP为核磁共振测井得到的T2谱,T2CUTCAL为利用形态法计算得到的T2截止值曲线,T2CUT33为以33ms作为T2截止值时计算得到的T2截止值曲线,T2CUT_CORE为岩心分析得到的T2截止值。从图上可以看出,利用形态法计算得到的T2截止值以及利用截止值计算的束缚水饱和度要好于以固定值33ms作为T2截止值的结果,形态法计算的结果跟接近岩心分析的结果,这说明,可以在没有取心的情况下,直接利用形态法来确定T2截止值而不需依赖岩心实验。
和现有的确定T2截止值的方法相比,形态法可以逐点计算出测井井段每个深度点上的T2截止值,而传统方法由于是对岩心做实验,所以只能是确定出取心点处的T2截止值。利用上述形态法以及自主开发的软件包,我们就可以对核磁共振测井数据的每个深度点上的测井T2谱确定出其T2截止值,这样就得到了一条连续的和测井井段相匹配的T2截止值曲线,进而利用该曲线可以计算出每个深度点上的束缚水饱和度以及渗透率等。利用形态法,我们对新疆等四个地区的岩心做了计算,图21中的两条边界线分别是以±10ms为误差的分界线,图22中的两条边界线分别是以±5%为误差的分界线,从图21和图22上可以看出利用形态法的结果可以很好的代替岩心分析的结果。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (4)
1、一种核磁共振测井T2谱T2截止值的确定方法,其特征在于包括以下步骤:
步骤1、对核磁共振测井T2谱按形态为双峰型、单峰型和三峰型;
所述双峰型分为双峰小孔型、双峰大孔型,单峰型分为单峰小孔型、单峰大孔型;
步骤2、根据不同类型的测井T2谱拟合得出其离心谱;
单峰大孔型测井T2谱T2截止值的计算方法为利用数学函数对大孔段进行拟合,其数学形式为:
其中:
A-------拟合T2谱中的幅值极大值,A值取T2分布谱中的幅值极大值;
B-------某个T2时间值,B值取幅度值减小到某个极小时的T2值;
C-------某个待定的参数;
x-----各组份的T2时间值与测井T2谱的T2时间值相同;
f(x)-----拟合得到的幅度值;
单峰大孔型测井T2谱T2截止值的计算方法为,将测井T2谱的面积减去拟合谱面积即为离心谱面积;
单峰小孔型测井T2谱T2截止值的计算方法为,假设孔隙中的流体都是束缚流体,将孔隙赋予一个固定的T2截止值;
双峰型测井T2谱T2截止值的计算方法为,利用正态函数在测井T2谱上拟合出一个离心谱,其数学形式为:
其中:
A为拟合T2谱幅度值,A值取T2分布谱中第一个峰的幅值,即使拟合的函数图形和T2分布谱有相同的幅度;
μ为拟合T2谱的期望值,μ值取第一个峰值所对应的T2时间值;
σ为拟合T2谱的方差;
x为各组份的T2时间值与测井T2谱的T2时间值相同;
f(x)为拟合得到的幅度值;
三峰型测井T2谱T2截止值的计算方法为利用正态函数在测井T2谱上拟合出一个离心谱,其数学形式为:
其中:
A为拟合T2谱幅度值,A值取T2分布谱中第一个峰的幅值,即使拟合的函数图形和T2分布谱有相同的幅度;
μ为拟合T2谱的期望值,μ值取第一个峰值所对应的T2时间值;
σ为拟合T2谱的方差;
x为各组份的T2时间值与测井T2谱的T2时间值相同;
f(x)为拟合得到的幅度值;
步骤3、将离心谱进行面积积分,得出离心谱的面积,所述面积的计算方法为龙贝格求积法进行面积积分;
步骤4、计算测井T2谱的面积;
步骤5、将测井T2谱的面积进行累加,当累加面积等于或与离心谱面积最相近时,此时对应的T2值即为T2截止值。
2、根据权利要求1所述的核磁共振测井T2谱T2截止值的确定方法,其特征在于:所述步骤5之后还具有步骤6、将离心谱的面积除以测井T2谱的面积,即为束缚水饱和度的值。
3、根据权利要求1所述的核磁共振测井T2谱T2截止值的确定方法,其特征在于:所述固定的T2截止值为33ms。
4、根据权利要求1所述的核磁共振测井T2谱T2截止值的确定方法,其特征在于:所述步骤5中的T2截止值的计算方法为,将测井T2谱的面积进行累加,当累加面积等于或与离心谱面积最相近时,此时对应的T2值即为T2截止值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100722329A CN100349013C (zh) | 2005-05-27 | 2005-05-27 | 核磁共振测井t2谱t2截止值的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100722329A CN100349013C (zh) | 2005-05-27 | 2005-05-27 | 核磁共振测井t2谱t2截止值的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1869733A CN1869733A (zh) | 2006-11-29 |
CN100349013C true CN100349013C (zh) | 2007-11-14 |
Family
ID=37443444
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005100722329A Expired - Fee Related CN100349013C (zh) | 2005-05-27 | 2005-05-27 | 核磁共振测井t2谱t2截止值的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100349013C (zh) |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102253069B (zh) * | 2011-04-20 | 2014-01-15 | 中国石油天然气股份有限公司 | 一种根据核磁共振t2谱确定渗透率的方法和装置 |
CN104278989B (zh) * | 2013-07-02 | 2017-02-15 | 中国石油天然气股份有限公司 | 一种获取低孔低渗储层饱和度指数的方法 |
CN103437759B (zh) * | 2013-08-09 | 2015-11-25 | 中国石油集团川庆钻探工程有限公司 | 非实验测量天然气层t2截止值的方法 |
CN103675722B (zh) * | 2013-11-27 | 2016-05-25 | 中国石油大学(华东) | 岩石t2-g实验采集参数自动匹配方法 |
CN104237957B (zh) * | 2014-10-20 | 2017-01-25 | 中国石油大学(华东) | 一种核磁共振测井t2截止值的确定方法 |
CN104316554B (zh) * | 2014-10-27 | 2016-07-06 | 中国石油天然气股份有限公司 | 储层孔隙水可动性测试方法 |
CN104932027A (zh) * | 2015-05-06 | 2015-09-23 | 中国石油大学(北京) | 基于核磁共振测井的储层分类方法 |
CN104990854B (zh) * | 2015-07-06 | 2017-07-07 | 中国石油天然气股份有限公司 | 确定束缚水饱和度的方法及装置 |
CN105629316B (zh) * | 2016-01-13 | 2018-10-12 | 中国石油大学(北京) | 获得原状地层的流体径向油水饱和度变化的方法及装置 |
CN105804734B (zh) * | 2016-04-12 | 2018-10-16 | 中国石油天然气集团公司 | 一种利用核磁共振测井识别稠油储层的方法 |
CN105866160B (zh) * | 2016-06-16 | 2017-09-01 | 中国石油大学(华东) | 一种基于纵波约束的核磁共振t2截止值计算方法 |
CN106324688B (zh) * | 2016-09-23 | 2018-04-17 | 中国石油大学(北京) | 一种储层束缚水饱和度确定方法及装置 |
CN106351652B (zh) * | 2016-11-25 | 2019-03-29 | 中国地质大学(北京) | 一种含烃储集层核磁共振测井t2谱形态校正方法 |
US10545105B2 (en) * | 2017-07-27 | 2020-01-28 | Saudi Arabian Oil Company | Estimating formation properties using saturation profiles |
CN110306977B (zh) * | 2018-03-27 | 2021-10-29 | 中国石油化工股份有限公司 | 核磁共振t2谱含烃校正方法及*** |
CN110850491B (zh) * | 2018-08-20 | 2021-11-02 | 中国石油天然气股份有限公司 | T2谱反演方法、装置和存储介质 |
CN110160927B (zh) * | 2019-05-16 | 2021-12-10 | 中海油田服务股份有限公司 | 一种确定束缚水饱和度的方法及装置 |
CN110410058B (zh) * | 2019-06-20 | 2023-01-31 | 中国石油化工集团有限公司 | 一种校正岩心实验结果刻度二维核磁测井的方法 |
CN110244369B (zh) * | 2019-06-28 | 2020-10-13 | 中国石油大学(北京) | 储层束缚和可动流体分布确定方法、装置及*** |
CN111289552A (zh) * | 2020-01-19 | 2020-06-16 | 中国地质大学(北京) | 煤的核磁共振截止值确定方法以及装置 |
CN113253354B (zh) * | 2021-05-12 | 2024-06-14 | 中国石油天然气集团有限公司 | 一种核磁t2谱截止值的确定方法、装置、介质及电子设备 |
CN114112856A (zh) * | 2021-12-07 | 2022-03-01 | 核工业北京化工冶金研究院 | 一种含铀砂岩全尺寸孔径结构表征方法及*** |
CN116859473A (zh) * | 2023-05-23 | 2023-10-10 | 贵州省油气勘查开发工程研究院 | 面向非常规储层含水饱和度的形态学和分形学预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1193387A (zh) * | 1995-06-26 | 1998-09-16 | 纽马公司 | 采用扩散和弛豫测井测量技术评估地岩的核磁共振***和方法 |
WO1999047939A1 (en) * | 1998-03-20 | 1999-09-23 | Numar Corporation | System and method for identification of hydrocarbons using enhanced diffusion |
CN1253294A (zh) * | 1998-11-05 | 2000-05-17 | 施卢默格控股有限公司 | 用于计算自旋-自旋驰豫时间分布的设备和方法 |
US6225803B1 (en) * | 1998-10-29 | 2001-05-01 | Baker Hughes Incorporated | NMR log processing using wavelet filter and iterative inversion |
US6316940B1 (en) * | 1999-03-17 | 2001-11-13 | Numar Corporation | System and method for identification of hydrocarbons using enhanced diffusion |
-
2005
- 2005-05-27 CN CNB2005100722329A patent/CN100349013C/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1193387A (zh) * | 1995-06-26 | 1998-09-16 | 纽马公司 | 采用扩散和弛豫测井测量技术评估地岩的核磁共振***和方法 |
WO1999047939A1 (en) * | 1998-03-20 | 1999-09-23 | Numar Corporation | System and method for identification of hydrocarbons using enhanced diffusion |
US6225803B1 (en) * | 1998-10-29 | 2001-05-01 | Baker Hughes Incorporated | NMR log processing using wavelet filter and iterative inversion |
CN1253294A (zh) * | 1998-11-05 | 2000-05-17 | 施卢默格控股有限公司 | 用于计算自旋-自旋驰豫时间分布的设备和方法 |
US6316940B1 (en) * | 1999-03-17 | 2001-11-13 | Numar Corporation | System and method for identification of hydrocarbons using enhanced diffusion |
Also Published As
Publication number | Publication date |
---|---|
CN1869733A (zh) | 2006-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100349013C (zh) | 核磁共振测井t2谱t2截止值的确定方法 | |
Ge et al. | Determination of nuclear magnetic resonance T 2 cutoff value based on multifractal theory—An application in sandstone with complex pore structure | |
CN103884633B (zh) | 一种确定岩石渗透率的方法及装置 | |
Chang et al. | Effective porosity, producible fluid and permeability in carbonates from NMR logging | |
CN102434152B (zh) | 一种储层含油饱和度的计算方法 | |
CN100504446C (zh) | 通过比较不同勘测深度处的核磁共振响应来探测烃的方法 | |
CN106066492B (zh) | 一种核磁共振测井孔隙度回波间隔影响的校正方法 | |
CN104075974A (zh) | 一种利用低场核磁共振精确测定页岩孔隙度的方法 | |
MX2008008999A (es) | Pelicula de capas multiples resistente a la intemperie. | |
US7603237B2 (en) | Method for analyzing data having shared and distinct properties | |
CN102042011B (zh) | 利用常规测井资料构建伪核磁t2谱的方法 | |
US10545105B2 (en) | Estimating formation properties using saturation profiles | |
CN105114064A (zh) | 确定致密砂岩储层饱和度的方法 | |
CN104777181A (zh) | 致密油核磁共振t2截止值及流体饱和度确定方法、装置 | |
CN106443800B (zh) | 一种中基性火山岩的核磁孔隙度校正方法 | |
CN107688037A (zh) | 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 | |
CN111827968A (zh) | 一种基于核磁共振测井的储层非均质性评价方法及装置 | |
CN205297569U (zh) | 用于确定致密砂岩饱和度指数的装置 | |
CN111077574A (zh) | 一种确定地层弹性参数的方法、装置及*** | |
CN108548765B (zh) | 一种变粘土骨架的孔隙度计算方法 | |
CN113267809A (zh) | I类页岩储层预测方法及装置 | |
CN113027421B (zh) | 油基钻井液环境下核磁共振测井t2谱形态校正方法 | |
US11821862B2 (en) | Method for measuring the spatial water permeability profile of porous media by using non-destructive nuclear magnetic resonance technique | |
Lowden | A New Method For Separating Lithologies And Estimating Thickness-Weighted Permeability Using Nmr Logs | |
CN114428049B (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 |
Granted publication date: 20071114 |