一种核磁共振测井T₂谱100%纯水谱的确定方法

一种核磁共振测井T₂谱100%纯水谱的确定方法
【专利摘要】本发明公开了一种核磁共振测井T2谱100％纯水谱的确定方法，包括：首先根据已知的岩石物理实验特征建立勘探地区储层束缚水饱和度Swi与T2几何均值T2g的关系，然后利用球管薄膜模型逐点计算测井T2谱束缚水饱和度Swi及T2几何均值T2g，并建立核磁共振离心谱，在核磁共振离心谱基础上构建100％纯水核磁共振谱。当储层含油气时核磁信号与100％纯水核磁信号存在显著的差异，利用这种差异对储层流体进行准确的识别。本发明利用核磁共振测井资料构造100％纯水谱能够反映储层真实的孔隙结构，且通过对比实测核磁共振谱与构造纯水谱的差异有效识别储层流体性质，回避了由于原油性质、观测模式导致的核磁共振测井流体识别困难的难题，提高了核磁共振测井资料的应用，增强了解决地质问题的能力。
【专利说明】
-种核磁共振测井T2谱100%纯水谱的确定方法
技术领域
[0001] 本发明设及岩屯、核磁共振、核磁共振测井及储层评价领域，特别设及一种核磁共 振测井T2谱100 %纯水谱的确定方法。
【背景技术】
[0002] 核磁共振是磁场中原子核对电磁波的一种响应。核磁共振测井的原始数据是幅度 随时间衰减的自旋回波串，自旋回波串包含了储层物性、孔隙类型、孔径大小、流体类型及 其分布等十分丰富的信息，可W表示为多个衰减指数的叠加：
[0003]
[0004] 其中:M(t)-一核磁共振测井测量得到的随时间t变化的回波幅度；
[0005] T2,i-第i个组分的横向弛豫时间，i = l，. . .，N;
[0006] Mi(0)-对应于T2,i组分的核磁共振孔隙度，i = l，...，N。
[0007] 通过适当的刻度，根据自旋回波串信号强度可W得到核磁孔隙度，并通过多指数 拟合自旋回波串，可W得到每一深度点的7'2,广--机分布曲线，即核磁共振测井T2谱。该测井 T2谱横坐标为T2时间值，纵坐标为孔隙度分量。
[000引核磁共振测井在评价储层品质方面具有独特的优势，通过解谱得到测井T2谱，可 W直观的反映孔径大小及其分布。当储层饱含水时，大孔喉对应较大T2值，小孔喉对应较小 T2值，从而达到评价孔隙结构的目的。但是当储层含油气时，计算的核磁孔隙度偏低，测井T2 谱的幅度和峰值位置不再反映储层孔喉分布，若直接使用测井T2谱进行储层品质的评价必 然是不恰当的，所W很有必要对测井T2谱进行含控校正。
[0009] 目前国内研究人员从核磁共振实验及测井T2谱出发对储层含控校正做了相关研 究工作。
[0010] ①优化观测模式。在原始数据采集上，可W采用特殊的观测模式W便提供孔隙度 校正信息。
[0011] ②数值模拟方法。参见2016年1月第2期《石油勘探与开发》杂志中，胡法龙、周灿灿 等著作的《核磁共振测井构造水谱流体识别方法及应用》的文章，其中记载了数值模拟方法 的应用，利用Iw和Te测量下的T2分布，构造 Iw长等待时间和Te下完全含水状态下T2分布，从而 达到储层评价的目的。
[0012] ③其他数学统计方法。参见2000年6月第3期《测井技术》杂志中，高敏等著作的《用 核磁共振测井资料评价储层的孔隙结构》的文章，其中记载了区域性T2g与各孔隙结构参数 关系，进一步消除储层含控的影响;参见2012年2月02期《石油天然气学报》杂志中，肖飞、何 宗斌等著作的《核磁共振测井连续表征储层孔隙结构方法研究》的文章，其中记载了通过测 井T2谱计算孔隙结构等相关参数，然后与实验值对比，最后利用含油饱和度作为校正系数 消除二者之间的误差;参见2010年石油技术交流会上，刘忠华、周灿灿等著作的《一种利用 核磁共振测井资料评价储层孔隙结构的新方法》的文章，其中记载了核磁=孔隙组分百分 含量与孔隙结构类型的关系，用w准确识别高产油气层，运种方法适用于多相流体存在的 地层。
[0013] W上提出的各种含控校正方法具有各自的优势和局限性。第1种方法的优点是可 W从原始数据采集上获得高质量的油气水分离信号，但其应用上的局限性是需要慎重进行 测前设计，但目前来说国内测前设计方面研究较少。方法2的优点是可W从根本上得到完全 含水的T2分布，但是算法等复杂，不易推广。最后一种方法区域应用性较强，不具有普遍规 律，没有从根本上解决含控对谱分布的影响。
[0014] 为了解决核磁测井应用中T2谱受含控影响不能很好反映储层品质的实际难题，应 该采用全新的研究思路与方法，避开岩屯、实验或其他信息标定的传统思路与方法，提出确 定核磁共振测井含控校正的一种可靠方法和技术，克服现行方法的技术局限性。

【发明内容】

[0015] 本发明的目的是利用测井T2谱构造100%纯水谱而实现含控校正，利于后续流体 识别及储层品质评价等。
[0016] 为实现上述目的，本发明提供了一种核磁共振测井T2谱100%纯水谱的确定方法， 包括W下步骤：
[0017] 步骤(1)通过岩屯、核磁共振实验获取实验核磁共振T2饱和谱和离屯、谱数据；
[001引步骤（2)根据步骤（1)得到的岩屯、核磁共振T2饱和谱和离屯、谱实验数据验数据计 算岩屯、束缚水饱和度Swi及核磁共振T2谱几何均值T2g并建立二者关系，采用最小二乘法确定 待定系数a、b;
[0019] 步骤(3)利用球管薄膜模型计算核磁共振T2谱的岩屯、束缚水饱和度Swi，并构建离 屯、谱；
[0020] 步骤(4)根据步骤(2)确定的待定系数a、b，计算含水饱和度Sw=l时的核磁共振T2 谱几何均值T2g;
[0021] 步骤(5)借助正态分布函数在核磁共振T2谱上构建一个自由水谱；
[0022] 步骤(6)，将步骤(3)构建的离屯、谱与步骤(5)构建的自由水谱相加得到一个总谱， 利用规划求解的方法使得构建的总谱孔隙度批和几何均值T2g/与步骤(3)中计算岩屯、束缚 水饱和度Swi产生的中间参数测井孔隙度與t及步骤(4)中确定的核磁共振T2谱几何均值T2g 值误差最小，此时得到的总谱即为100%纯水谱。
[0023] 进一步的，所述步骤（1)中，核磁共测井振T2实验按照《岩样核磁共振参数实验室 测量规范SY/T 6490-2007》标准规定的流程进行。
[0024] 进一步的，所述步骤(2)中岩屯、束缚水饱和度Swi及核磁共振实验T2谱几何均值T2g， 具体按照如下公式进行计算：
[0025]
[0026]
[0027]其中：與细相为岩屯、离屯、谱孔隙分量;取幼岩屯、饱和谱孔隙分量;T2,功实验T2谱的 横向弛豫时间。
[002引进一步的，所述步骤(2)中，建立岩屯、束缚水饱和度Swi及核磁共振T2谱几何均值T2g 的关系如下：
[0029]
[0030] 进一步的，在进行步骤(3)操作时，需对核磁测井资料处理得到标准T2谱，处理方 法对于P型核磁参考《DPP核磁处理技术手册》。
[0031] 在进行步骤(3)操作之前，还需利用中子-密度曲线进行测井孔隙度祭t的计算，通 过下式实现：
[0032]
[003；3]其中：與DEJV为测井密度孔隙度;机:亂为测井中子孔隙度。
[0034] 所述测井孔隙度取的计算按照石油天然气行业标准《生产测井解释规程SY/ T5793-1993》规定的流程进行。
[0035] 进一步的，所述步骤(3)利用球管薄膜模型计算测井T2谱的束缚水饱和度Swi并构 建离屯、谱，是根据球管薄膜模型确定的权系数计算得到：
[0036]
[0037]
[00；3 引
[0039]
[0040] 其中：取swU为离屯、谱分量；辦如为束缚水孔隙度;W功权系数分量;m为系数;Swi 为束缚水饱和度;取t为测井孔隙度;沪i为测井T2谱孔隙分量;T2,i为测井T2谱的横向弛豫时 间。
[0041] 进一步的，所述步骤(5)借助正态分布函数在测井T2谱上构建一个自由水谱：
[0042]
[0043] 其中：A为构建自由水谱的幅度值;U为构建自由水谱的期望值;0为构建自由水谱 的标准差;X为构建自由水谱各组分的时间值，与测井T2谱时间一致;f(x)为构建自由水谱 的幅值。
[0044] 本发明将自由水谱与离屯、谱相加得到一个总谱，利用规划求解方法使得构建的总 谱孔隙度取f和测井T2谱几何均值T2g/与步骤(3)中的测井孔隙度餐t及步骤(4)中确定的T2g 值误差最小，此时得到的总谱即为100%纯水谱，
[0045] 本发明实现了多相流体存在情况下利用核磁共振测井资料构造100%纯水谱，纯 水谱能够反映储层真实的孔隙结构，并且通过对比实测核磁共振谱与构造纯水谱的差异能 够有效识别储层流体性质，在构建水谱基础上可W进行储层品质的评价，回避了由于原油 性质、观测模式等导致的核磁共振测井流体识别困难的难题，并且可W通过对比构造水谱 与实测谱之间的差异达到流体识别的目的。提高了核磁共振测井资料的应用，增强了解决 地质问题的能力。
【附图说明】
[0046] 图1为本发明一种核磁共振测井T2谱100%纯水谱的确定方法流程示意图。
[0047] 图2为本发明一种核磁共振测井T2谱100 %纯水谱的确定方法实施例T2g确定方法 示意图。
[004引图3为本发明一种核磁共振测井T2谱100%纯水谱的确定方法实施例m值确定方法 示意图。
[0049] 图4为本发明一种核磁共振测井T2谱100%纯水谱的确定方法实施例单峰型示意 图。
[0050] 图5本发明一种核磁共振测井T2谱100%纯水谱的确定方法实施例双峰型示意图。
[0051] 图6本发明一种核磁共振测井T2谱100%纯水谱的确定方法实施例双峰型示意图。
[0052] 图7为本发明一种核磁共振测井T2谱100%纯水谱的确定方法计算机流程图。
[0053] 图8为本发明一种核磁共振测井T2谱100%纯水谱的确定方法某地区某井实际处 理的效果图。
[0054] 图9为本发明一种核磁共振测井T2谱100%纯水谱的确定方法某地区计算的孔隙 度与岩屯、分析孔隙度结果对比图。
[0055] 图10为本发明一种核磁共振测井T2谱100%纯水谱的确定方法某地区计算的T2g与 岩屯、分析T2g结果对比图。
【具体实施方式】
[0056] 下面结合附图和实施例对发明作进一步的详细说明，但并不作为对发明做任何限 制的依据。
[0057] 通过对核磁共振T2谱的形态研究，在大量实验数据和实测资料的基础上，提出了 一种核磁共振测井T2谱100%纯水谱的确定方法，该方法的流程示意图如图1所示。为使本 发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详 细描述。
[005引第一步，按照《岩样核磁共振参数实验室测量规范SY/T 6490-2007》标准规定的流 程进行岩屯、核磁T2实验，获取实验核磁共振&谱饱和谱和离屯、谱数据。
[0059] 第二步，根据第一步得到的岩屯、核磁共振T2饱和谱和离屯、谱实验数据计算束缚水 饱和度Swi及核磁共振T2谱几何均值T2g，并建立二者关系，如图2所示，具体确定方法如下：
[0060]
[0061]
[0062]
[0063] 其中：取Swa为岩屯、离屯、谱孔隙分量;斬为岩屯、饱和谱孔隙分量;T2,i为实验T2谱的 横向弛豫时间；
[0064] 具体实施过程中，采用最小二乘法确定待定系数a、b，a值取0.878，b值取1.149。
[0065] 第=步，利用球管薄膜模型计算测井T2谱的束缚水饱和度Swi并构建离屯、谱，具体 确定方了-
[0066]
[0067]
[006引
[0069]
[0070]
[007。 其中：從碱i为离屯、谱分量；從wi为束缚水孔隙度;Wi为权系数分量;m为系数；Swi 为束缚水饱和度;辑￡为测井孔隙度;取呪酣为测井密度孔隙度;與为测井中子孔隙度;特左 为测井T2谱孔隙分量;T2,功测井T2谱的横向弛豫时间。
[0072] 我们发现，利用球管模型时权系数曲线更加平滑，取值更加合理：当T2小于某一个 值时，孔隙空间充满了束缚流体，此时权系数为1，随着孔喉半径的增大，不可动流体所占空 间比例不断减小，权系数趋近于0,权系数是T2值的函数。对于要准确构造离屯、谱，则上述参 数m值取值很关键，我们通过大量实验数据分析，发现m值与T2截止值T2cutDff有很好的线性关 系，即：
[0073] m = c*T2cut〇ff+d
[0074] 如图3所示，具体实施过程中，采用最小二乘法确定待定系数系数c、d，c值取 0.196，d值取-0.364;当缺少地区经验T2截止值时，默认m取2.46。
[0075] 第四步，根据第二步中确定的T2g~Swi关系W及待定系数a、b，计算細=1时的核磁 共振T2谱几何均值T2g。
[0076] 第五步，借助正态分布函数在测井T2谱上构建一个自由水谱：
[0077]
[0078] 其中：A为构建自由水谱的幅度值;U为构建自由水谱的期望值;0为构建自由水谱 的标准差;X为构建自由水谱各组分的时间值，与测井T2谱时间一致;f(x)为构建自由水谱 的幅值。
[0079] 第六步，将第五步构建的自由水谱与第=步构建的离屯、谱相加得到一个总谱，具 体实施过程中，我们采用规划求解的方法，使得构建的总谱孔隙度取寸日核磁共振几何均值 T2g/与第=步中的测井孔隙度巧t及第四步中确定的T2g值误差最小，此时得到的总谱即为 100 %纯水谱，即：
[0080]
[0081]
[0082] 对于要准确构造纯水谱，则上述A、u、〇 =个参数的确定也是至关重要，在具体实施 过程中，采用规划求解方法确定，但是通过大量数据分析发现，=个参数都有一个明确的取 值范围，A值取值在0.35附近，U值取值在1.70附近，0取值在0.44附近。
[0083] 按照W上思路，利用构造水谱法我们对某地区的部分岩屯、做了计算，如图4-6,岩 屯、数据如表1。
[0084] 同时我们在geolog中编写了该方法的程序，实现了自动处理。流程图如图7所示。
[0085] 表 1 「OORAl
[0087]
[0088] 利用构造水谱法自动构造核磁共振测井T2谱100%纯水谱方法W及自主开发的软 件包，我们对某油田某井进行了处理，如图8所示为某油田某井实际处理的效果图，处理的 效果图如下：
[0089] 图中第一道为CAL(井径曲线）、SP(自然电位曲线)和GR(自然伽马曲线），第二道为 深度，第=道为=孔隙度曲线，第四道为电阻率曲线，第五道为实测T2谱孔隙分布谱，第六 道和第屯道分别为构造离屯、谱与纯水谱。从图上可W看出，利用构造水谱法得到的100%纯 水谱在含水地层与实测谱形态基本一致，说明该方法的可信性。
[0090] 和现有的通过含控校正得到纯水谱的方法相比，构造水谱法可W快速便捷的逐点 计算每个深度点的100%纯水谱，受流体类型等的影响较小。利用自主开发的软件包，我们 就可W得到一系列连续的和测井井段相匹配的纯水谱，进而利用该孔隙分布更加准确的对 储层品质及流体类型进行评价。利用构造法，我们对某地区的岩屯、做了计算，图9中的两条 边界线分别是±5%为误差的分界线，图10中的两条边界线是W±5%为误差分界线，从图9 和图10上可W看出构造法的结果可W很好的替代岩屯、分析的结果。
[0091] 最后所应说明的是，W上实施例仅用W说明本发明的技术方案而非限制，尽管参 照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可W对本发明 的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均涵盖在 本发明的权利要求范围当中。
【主权项】
1. 一种核磁共振测井1'2谱100 %纯水谱的确定方法，其特征在于，所述方法包括以下步 骤： 步骤(1 )，对岩心进行核磁共振T2实验，获取核磁共振!^饱和谱和离心谱实验数据； 步骤(2)，根据步骤（1)得到的岩心核磁共振T2饱和谱和离心谱实验数据计算岩心束缚 水饱和度SW1及核磁共振T2谱几何均值T2g，并建立二者关系，采用最小二乘法确定待定系数 a、b; 步骤（3)，利用球管薄膜模型计算核磁共振^谱的岩心束缚水饱和度SW1，并构建离心 谱； 步骤(4)，根据步骤(2)确定的待定系数a、b，计算含水饱和度Sw=l时的核磁共振T2谱几 何均值T2g; 步骤(5)，借助正态分布函数在核磁共振!^谱上构建一个自由水谱； 步骤(6)，将步骤(3)构建的离心谱与步骤(5)构建的自由水谱相加得到一个总谱，利用 规划求解的方法使得构建的总谱孔隙度#和几何均值T2f/与步骤(3)中计算岩心束缚水饱 和度SW1产生的中间参数测井孔隙度炉t及步骤(4)中确定的核磁共振T 2谱几何均值T2g值误 差最小，此时得到的总谱即为100 %纯水谱。2. 如权利要求1所述的方法，其特征在于，所述步骤(1)中，岩心核磁共振^实验按照《岩 样核磁共振参数实验室测量规范SY/T 6490-2007》标准规定的流程进行。3. 如权利要求1所述的方法，其特征在于，所述步骤(2)中，岩心束缚水饱和度SW1及核磁 共振T2谱几何均值T 2g，具体按照如下公式进行计算：其中：识Swi,i为离心谱孔隙分量；_为饱和谱孔隙分量;T2>1为实验^谱的横向弛豫时 间。4. 如权利要求1所述的方法，其特征在于，所述步骤(2)中，建立岩心束缚水饱和度SW1及 核磁共振!^谱几何均值T 2g的关系如下：5. 如权利要求1所述的方法，其特征在于，在进行步骤(3)操作之前，还需利用中子-密 度曲线进行测井孔隙度街的计算，通过下式实现：其中：识0SW为测井密度孔隙度;识C継为测井中子孔隙度。6. 如权利要求5所述的方法，其特征在于，所述测井孔隙度识t的计算按照石油天然气行 业标准《生产测井解释规程SY/T 5793-1993》规定的流程进行。7. 如权利要求1所述的方法，其特征在于，所述步骤(3)中计算的束缚水饱和度并构建 离心谱，是根据球管薄膜模型确定的权系数计算得到的：其中，炉swi,i力离心谱分量;为束缚水孔隙度;Wi为权系数分量;m为系数;Swi为束缚 水饱和度;炉t为测井孔隙度;#i为测井T2谱孔隙分量;T2, i为测井T2谱的横向弛豫时间。8. 如权利要求1所述的方法，其特征在于，所述步骤(5)借助正态分布函数在核磁共振!^ 谱上构建一个自由水谱：其中：A为构建自由水谱的幅度值;u为构建自由水谱的期望值;〇为构建自由水谱的标 准差;x为构建自由水谱各组分的时间值，与测井!^谱时间一致;f(x)为构建自由水谱的幅 值。9. 如权利要求1所述的方法，其特征在于，通过所述步骤(6)中建立的100%纯水谱进行 储层品质评价，并通过与测井!^谱对比进行储层流体识别。
【文档编号】E21B49/08GK106050225SQ201610394292
【公开日】2016年10月26日
【申请日】2016年6月6日
【发明人】王明方, 万金彬, 李慧莹, 程道解, 白松涛, 黄科, 许思勇, 金小慧, 王琪, 张旭, 黄娅, 田扬, 赵建斌, 任龙, 曾静波, 何羽飞, 袁野
【申请人】中国石油天然气集团公司, 中国石油集团测井有限公司