一种核磁测井数据时间域分析的计算方法技术

本发明专利技术公开了一种核磁测井数据时间域分析的计算方法，该方法包括：步骤1，选择核磁测井长短极化时间测量的两组回波串数据进行相减，得到回波串差；步骤2，确定核磁测井回波串差符合的指数拟合模型，利用非常快速模拟退火算法对回波串差进行指数拟合，得到流体相T

A computational method for time domain analysis of nuclear magnetic logging data

The invention discloses a method for analyzing and computing NMR logging data in time domain, the method includes: Step 1, choose the length of time of two groups of NMR polarization measurements on echo data were subtracted, get echo; step 2, determine the NMR echo difference index fitting model with the use of very fast simulated annealing the algorithm of echo on exponential fitting by fluid phase T

【技术实现步骤摘要】
一种核磁测井数据时间域分析的计算方法
本专利技术涉及地球物理测井
，特别是涉及到核磁测井资料处理解释应用领域。
技术介绍
流体分析与计算是测井的基本任务之一，核磁测井是一种新型的裸眼井测井方法，它能为地层的油气评价提供与岩性无关的孔隙度、毛管束缚水、泥质束缚水、可动流体饱和度、渗透率等。跟常规测井相比，它提高了对地层油气评价的准确性，具有独特的技术优势。核磁共振分为极化和弛豫两个过程。地层中的氢原子核可以看作一个个小磁针，自然状态下，它的方向随机分布，对外不显极性，宏观极化矢量M为零。施加磁场进行极化时，小磁针逐渐定向排列，一段时间后，宏观极化矢量M达到最大值，该过程的时间常数称为纵向弛豫时间T1。磁场的作用时间称为极化时间TW。根据宏观极化矢量M建立的数学模型，等待时间TW大于3倍的纵向弛豫时间T1，就能保证极化率在95％以上。极化完成后，施加一个90°偏转脉冲，宏观极化矢量开始衰减。此时，再施加一个180°反向脉冲，宏观极化矢量反向汇聚，形成一个可观测的信号。重复这个过程得到回波串，相邻两次观测信号的时间称为回波间隔TE。回波串按照指数规律衰减，该过程的时间常数为横向弛豫时间T2。通常，地层岩石的T2不是一个值，而是一个分布谱，称为T2谱。双极化时间TW测井资料的处理解释的主要根据轻烃(天然气和轻质油)与水的纵向弛豫时间T1的差异。通常，轻烃有比较长的T1，而水T1较短。对于水而言，较短的极化时间就足以使其完全磁化；而轻质油与天然气则需要较长的极化时间。所以，如果有轻烃存在，长、短极化时间得到的T2谱就会有明显差异。两个T2谱相减，水的信号可以基本相互抵消，而油与气的信号则余留在差谱之中，由此识别和计算油气。国外核磁处理软件(哈里伯顿公司DPP软件)时间域分析模块包括两个步骤：T2T1_Event、TDA_comp，但是，受软件平台的限制，影响了该项技术的应用。为此，专利技术了一种核磁测井数据时间域分析的数学方法，解决了以上技术问题。
技术实现思路
本专利技术的目的是提供一种针对双极化时间(TW)核磁共振测井资料，进行定性识别和定量计算的核磁测井数据时间域分析的数学方法。本专利技术的目的可通过如下技术措施来实现：核磁测井数据时间域分析的数学方法，该核磁测井数据时间域分析的数学方法包括：步骤1，选择核磁测井长短极化时间测量的两组回波串数据进行相减，得到回波串差；步骤2，确定核磁测井回波串差符合的指数拟合模型，利用非常快速模拟退火算法对回波串差进行指数拟合，得到流体相T2特征值及视孔隙度；步骤3，选取目的层段流体相T2特征值的平均值作为储层流体相的T2特征值，设成常量，重新计算流体相的视孔隙度；以及步骤4，对流体相的视孔隙度进行极化校正得到流体相的真孔隙度。本专利技术的目的还可通过如下技术措施来实现：在步骤2中，对于单相流体，即油或气，利用单指数拟合，两相流体，即油和气，利用双指数拟合。在步骤2中，假设存在单相流体，设为X相，则按照(1)式进行单指数拟合：其中，EDIF(t)为回波串差，t为时间，φx1为X相流体孔隙度，T2x1为X相流体T2特征值；假设存在两相流体，设为X、Y相，则按照(2)式进行双指数拟合：其中，EDIF(t)为回波串差，t为时间，φx为X相流体孔隙度，φy为Y相流体孔隙度，T2x为X相流体T2特征值，T2y为Y相流体T2特征值。步骤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，其中，t0是初始温度，M为回波串个数，α为降温系数；(f)重复步骤b到步骤e，直到满足所有的收敛条件为止，若执行结果超过若干步未发生变化，则退出循环，输出结果。在步骤a中，单相T2x1特征值范围为32ms-2000ms，孔隙度范围为1％-20％，两相T2x、T2y特征值范围分别为32ms-200ms，200-2000ms，孔隙度范围均为1％-20％。单相计算的目标函数为：两相计算的目标函数为：采用a参数代表φx1，b参数代表T2x1，其初始解为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，t0是初始温度，设定为100，M为回波串个数，根据实际回波串个数设定，α为降温系数，设定为0.98。步骤4包括：根据核磁测井响应机理可知，视孔隙度和真孔隙度的关系为：其中，分别为油和气的视孔隙度，φo、φg分别为油和气的真孔隙度，HIo、HIg分别为油和气的含氢指数，T1o、T1g分别为油气的纵向弛豫时间；Twl、Tws分别是核磁测井观测模式中的极化时间；在得到油和气的视孔隙度后，根据式(3)、(4)，计算出油气真孔隙度φo、φg。本专利技术中的核磁测井数据时间域分析的数学方法，研究分析了时间域分析的数学模型，建立了通过指数拟合来进行油气识别和计算的方法。在没有国外软件的情况下，采用本专利技术可以进行时间域分析，计算得到轻质油气的真孔隙度，实现了核磁测井数据时间域分析的快速处理解释。附图说明图1为本专利技术的核磁测井数据时间域分析的数学方法一具体实施例的流程图；图2为本专利技术的一具体实施例中两相T2搜索结果与计算示意图；图3为本专利技术的一具体实施例中两相孔隙度计算结果对比示意图。具体实施方式为使本专利技术的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合所附图式，作详细说明如下。如图1所示，图1为本专利技术的核磁测井数据时间域分析的数学方法一具体实施例的流程图。在步骤101，选择核磁测井长短极化时间测量的两组回波串数据进行相减，得到回波串差。流程进入到步骤102。在步骤102，利用非常快速模拟退火算法对回波串差进行指数拟合，得到流体相T2特征值及视孔隙度，其中，单相流体(油或气)利用单指数拟合，两相流体(油和气)利用双指数拟合。假设存在单相流体，设为X相，则按照(1)式进行单指数拟合：其中，EDIF(t)为回波串差，φx1为X相流体孔隙度，T2x1为X相流体T2特征值。假设存在两相流体，设为X、Y相，则按照(2)式进行双指数拟合：其中，EDIF(t)为回波串差，t为时间，φx为X相流体孔隙度，φy为Y相流体孔隙度，T2x为X相流体T2特征值，T2y为Y相流体T2特征值。流程进入到步骤103。在步骤103，选取目的层段流体相T2特征值的平均值本文档来自技高网...

【技术保护点】
核磁测井数据时间域分析的数学方法，其特征在于，该核磁测井数据时间域分析的数学方法包括：步骤1，选择核磁测井长短极化时间测量的两组回波串数据进行相减，得到回波串差；步骤2，确定核磁测井回波串差符合的指数拟合模型，利用非常快速模拟退火算法对回波串差进行指数拟合，得到流体相T

【技术特征摘要】
1.核磁测井数据时间域分析的数学方法，其特征在于，该核磁测井数据时间域分析的数学方法包括：步骤1，选择核磁测井长短极化时间测量的两组回波串数据进行相减，得到回波串差；步骤2，确定核磁测井回波串差符合的指数拟合模型，利用非常快速模拟退火算法对回波串差进行指数拟合，得到流体相T2特征值及视孔隙度；步骤3，选取目的层段流体相T2特征值的平均值作为储层流体相的T2特征值，设成常量，重新计算流体相的视孔隙度；步骤4，对流体相的视孔隙度进行极化校正，得到流体相的真孔隙度。2.根据权利要求1所述的核磁测井数据时间域分析的数学方法，其特征在于，在步骤2中，对于单相流体，即油或气，利用单指数拟合，两相流体，即油和气，利用双指数拟合。3.根据权利要求2所述的核磁测井数据时间域分析的数学方法，其特征在于，在步骤2中，假设存在单相流体，设为X相，则按照(1)式进行单指数拟合：其中，EDIF(t)为回波串差，t为时间，φx1为X相流体孔隙度，T2x1为X相流体T2特征值；假设存在两相流体，设为X、Y相，则按照(2)式进行双指数拟合：其中，EDIF(t)为回波串差，t为时间，φx为X相流体孔隙度，φy为Y相流体孔隙度，T2x为X相流体T2特征值，T2y为Y相流体T2特征值。4.根据权利要求3所述的核磁测井数据时间域分析的数学方法，其特征在于，步骤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，其中，t0是初始温度，M为回波串个数，α为降温系数；(f)重复...

【专利技术属性】
技术研发人员：李忠新，郭迎春，陈利，刘磊，李志鹏，王玮，王瑞，刘峰，高劲松，张书凡，孙玉花，何伟，段悦，赵林，单宝忠，杨明林，李坚，程培涛，班丽，姜官波，刘阳，李靓，陈云佳，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司胜利油田分公司勘探开发研究院，
类型：发明
国别省市：山东,37