一种核磁测井数据时间域分析的计算方法与流程

本发明涉及地球物理测井技术领域，特别是涉及到核磁测井资料处理解释应用领域。

背景技术：

流体分析与计算是测井的基本任务之一，核磁测井是一种新型的裸眼井测井方法，它能为地层的油气评价提供与岩性无关的孔隙度、毛管束缚水、泥质束缚水、可动流体饱和度、渗透率等。跟常规测井相比，它提高了对地层油气评价的准确性，具有独特的技术优势。

核磁共振分为极化和弛豫两个过程。地层中的氢原子核可以看作一个个小磁针，自然状态下，它的方向随机分布，对外不显极性，宏观极化矢量M为零。施加磁场进行极化时，小磁针逐渐定向排列，一段时间后，宏观极化矢量M达到最大值，该过程的时间常数称为纵向弛豫时间T₁。磁场的作用时间称为极化时间TW。根据宏观极化矢量M建立的数学模型，等待时间TW大于3倍的纵向弛豫时间T₁，就能保证极化率在95％以上。极化完成后，施加一个90°偏转脉冲，宏观极化矢量开始衰减。此时，再施加一个180°反向脉冲，宏观极化矢量反向汇聚，形成一个可观测的信号。重复这个过程得到回波串，相邻两次观测信号的时间称为回波间隔TE。回波串按照指数规律衰减，该过程的时间常数为横向弛豫时间T₂。通常，地层岩石的T₂不是一个值，而是一个分布谱，称为T₂谱。

双极化时间TW测井资料的处理解释的主要根据轻烃(天然气和轻质油)与水的纵向弛豫时间T₁的差异。通常，轻烃有比较长的T₁，而水T₁较短。对于水而言，较短的极化时间就足以使其完全磁化；而轻质油与天然气则需要较长的极化时间。所以，如果有轻烃存在，长、短极化时间得到的T₂谱就会有明显差异。两个T₂谱相减，水的信号可以基本相互抵消，而油与气的信号则余留在差谱之中，由此识别和计算油气。

国外核磁处理软件(哈里伯顿公司DPP软件)时间域分析模块包括两个步骤：T2T1_Event、TDA_comp，但是，受软件平台的限制，影响了该项技术的应用。为此，发明了一种核磁测井数据时间域分析的数学方法，解决了以上技术问题。

技术实现要素：

本发明的目的是提供一种针对双极化时间(TW)核磁共振测井资料，进行定性识别和定量计算的核磁测井数据时间域分析的数学方法。

本发明的目的可通过如下技术措施来实现：核磁测井数据时间域分析的数学方法，该核磁测井数据时间域分析的数学方法包括：步骤1，选择核磁测井长短极化时间测量的两组回波串数据进行相减，得到回波串差；步骤2，确定核磁测井回波串差符合的指数拟合模型，利用非常快速模拟退火算法对回波串差进行指数拟合，得到流体相T₂特征值及视孔隙度；步骤3，选取目的层段流体相T₂特征值的平均值作为储层流体相的T₂特征值，设成常量，重新计算流体相的视孔隙度；以及步骤4，对流体相的视孔隙度进行极化校正得到流体相的真孔隙度。

本发明的目的还可通过如下技术措施来实现：

在步骤2中，对于单相流体，即油或气，利用单指数拟合，两相流体，即油和气，利用双指数拟合。

在步骤2中，假设存在单相流体，设为X相，则按照(1)式进行单指数拟合：

其中，EDIF(t)为回波串差，t为时间，φ_x1为X相流体孔隙度，T_2x1为X相流体T₂特征值；

假设存在两相流体，设为X、Y相，则按照(2)式进行双指数拟合：

其中，EDIF(t)为回波串差，t为时间，φ_x为X相流体孔隙度，φ_y为Y相流体孔隙度，T_2x为X相流体T₂特征值，T_2y为Y相流体T₂特征值。

步骤2包括：

(a)给定模型参数一个变化范围，在这个范围内产生一个初始解，作为当前最优解，并计算目标函数值；

(b)采用依赖于温度Cauchy分布产生新模型，即对当前最优解做一随机扰动，产生一个新解，并计算新解的目标函数值；

(c)计算目标函数值增量ΔE，若ΔE<0，则接受新解为目前最优解；若ΔE>0，则以概率P＝[1-(1-h)ΔE/t]^1/(1-h)进行接收，h为接收概率公式的参数；

(d)在温度t下，重新一定次数的扰动和接收过程，即重复步骤b到步骤c；

(e)若内循环次数达到设定值N，则依据温度更新函数降低温度t，其中，t₀是初始温度，M为回波串个数，α为降温系数；

(f)重复步骤b到步骤e，直到满足所有的收敛条件为止，若执行结果超过若干步未发生变化，则退出循环，输出结果。

在步骤a中，单相T_2x1特征值范围为32ms-2000ms，孔隙度范围为1％-20％，两相T_2x、T_2y特征值范围分别为32ms-200ms，200-2000ms，孔隙度范围均为1％-20％。单相计算的目标函数为：两相计算的目标函数为：采用a参数代表φ_x1，b参数代表T_2x1，其初始解为a_I、b_I，初始解的目标函数值为：F_I＝|a_I×e^-t/b_I-EDIF(t)|。

在步骤b中，新解各参数按照公式a_pre＝a_I+(a_E-a_S)*myrand产生；其中，a_pre为a参数新解，a_I为a参数初始解，a_S、a_E为a参数范围的最小、最大值，myrand为随机数，ΔE＝F_pre-F_I，其中，F_pre为新解的目标函数值，F_I为初始解目标函数值。

在步骤c中，h为接收概率公式的参数，设定为-5。

在步骤e中，N为内循环循环截止次数，设定为3，t₀是初始温度，设定为100，M为回波串个数，根据实际回波串个数设定，α为降温系数，设定为0.98。

步骤4包括：

根据核磁测井响应机理可知，视孔隙度和真孔隙度的关系为：

其中，分别为油和气的视孔隙度，φ_o、φ_g分别为油和气的真孔隙度，HI_o、HI_g分别为油和气的含氢指数，T_1o、T_1g分别为油气的纵向弛豫时间；Twl、Tws分别是核磁测井观测模式中的极化时间；

在得到油和气的视孔隙度后，根据式(3)、(4)，计算出油气真孔隙度φ_o、φ_g。

本发明中的核磁测井数据时间域分析的数学方法，研究分析了时间域分析的数学模型，建立了通过指数拟合来进行油气识别和计算的方法。在没有国外软件的情况下，采用本发明可以进行时间域分析，计算得到轻质油气的真孔隙度，实现了核磁测井数据时间域分析的快速处理解释。

附图说明

图1为本发明的核磁测井数据时间域分析的数学方法一具体实施例的流程图；

图2为本发明的一具体实施例中两相T₂搜索结果与计算示意图；

图3为本发明的一具体实施例中两相孔隙度计算结果对比示意图。

具体实施方式

为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合所附图式，作详细说明如下。

如图1所示，图1为本发明的核磁测井数据时间域分析的数学方法一具体实施例的流程图。在步骤101，选择核磁测井长短极化时间测量的两组回波串数据进行相减，得到回波串差。流程进入到步骤102。

在步骤102，利用非常快速模拟退火算法对回波串差进行指数拟合，得到流体相T₂特征值及视孔隙度，其中，单相流体(油或气)利用单指数拟合，两相流体(油和气)利用双指数拟合。假设存在单相流体，设为X相，则按照(1)式进行单指数拟合：

其中，EDIF(t)为回波串差，φ_x1为X相流体孔隙度，T_2x1为X相流体T₂特征值。

假设存在两相流体，设为X、Y相，则按照(2)式进行双指数拟合：

其中，EDIF(t)为回波串差，t为时间，φ_x为X相流体孔隙度，φ_y为Y相流体孔隙度，T_2x为X相流体T₂特征值，T_2y为Y相流体T₂特征值。流程进入到步骤103。

在步骤103，选取目的层段流体相T₂特征值的平均值作为储层流体相的T₂特征值，设成常量，重新计算流体相的视孔隙度。具体包括以下步骤：

(a)给定模型参数一个变化范围，在这个范围内产生一个初始解，作为当前最优解，并计算目标函数值。

(b)采用依赖于温度Cauchy分布产生新模型，即对当前最优解做一随机扰动，产生一个新解，并计算新解的目标函数值；

(c)计算目标函数值增量ΔE，若ΔE<0，则接受新解为目前最优解；若ΔE>0，则以概率P＝[1-(1-h)ΔE/t]^1/(1-h)进行接收，h为接收概率公式的参数；

(d)在温度t下，重新一定次数的扰动和接收过程，即重复步骤b到步骤c；

(e)若内循环次数达到设定值N，则依据温度更新函数降低温度t，其中，t₀是初始温度，M为回波串个数，α为降温系数；

(f)重复步骤b-步骤e，直到满足所有的收敛条件为止，若执行结果超过若干步未发生变化，则退出循环，输出结果。

所述步骤a中，单相T_2x1特征值范围为32ms-2000ms，孔隙度范围为1％-20％，两相T_2x、T_2y特征值范围分别为32ms-200ms，200-2000ms，孔隙度范围均为1％-20％。单相计算的目标函数为：两相计算的目标函数为：例如，a参数代表φ_x1，b参数代表T_2x1，其初始解为a_I、b_I，初始解的目标函数值为：

F_I＝|a_I×e^-t/b_I-EDIF(t)|。

所述步骤b中，新解各参数按照公式a_pre＝a_I+(a_E-a_S)*myrand产生。其中，a_pre为a参数新解，a_I为a参数初始解，a_S、a_E为a参数范围的最小、最大值，myrand为随机数，ΔE＝F_pre-F_I，其中，F_pre为新解的目标函数值，F_I为初始解目标函数值。

所述步骤c中，h为接收概率公式的参数，一般设定为-5。

所述步骤e中，N为内循环循环截止次数，一般设定为3，t₀是初始温度，一般设定为100，M为回波串个数，根据实际回波串个数设定，α为降温系数，一般设定为0.98。流程进入到步骤104。

在步骤104，对流体相的视孔隙度进行极化校正得到流体相的真孔隙度。根据核磁测井响应机理可知，视孔隙度和真孔隙度的关系为：

其中，分别为油和气的视孔隙度，φ_o、φ_g分别为油和气的真孔隙度，HI_o、HI_g分别为油和气的含氢指数，T_1o、T_1g分别为油气的纵向弛豫时间。

因此，步骤104得到油气的相孔隙度后，根据式(3)、(4)，就可以计算出油气真孔隙度φ_o、φ_g。

图2是两相搜索与计算的结果。第一道为深度道，第二、三道为长短极化时间下测量的T₂谱。第四道为回波串差，第五道为差谱道，点线为X相T₂搜索结果，实线为Y相T₂搜索结果，第六道点线曲线为X相孔隙度结果，实线的为Y相孔隙度搜索结果。第五道粗点线为设定的X相T₂值，粗实线为设定的Y相T₂值，分别为100ms和800ms。最右一道为重新计算的两相计算的视孔隙度。

图3为两相孔隙度结算的结果。第一道为深度道，第二、三道为长短TW下的T₂谱。第四道为未进行时间域分析前的结果，第五、六道为本发明与国外软件DPP时间域分析后孔隙度结果。通过对比可以看出，本发明与国外软件孔隙度计算结果一致，符合地质实际情况。