核磁共振测井资料处理及解释规范.docx
- 文档编号:25491159
- 上传时间:2023-06-09
- 格式:DOCX
- 页数:26
- 大小:129.11KB
核磁共振测井资料处理及解释规范.docx
《核磁共振测井资料处理及解释规范.docx》由会员分享,可在线阅读,更多相关《核磁共振测井资料处理及解释规范.docx(26页珍藏版)》请在冰豆网上搜索。
核磁共振测井资料处理及解释规范
核磁共振测井资料处理及解释规范
I范围
本标准规定了MRIL-C型、MRIL-C/TP型和MRIL-Prime型核磁共振测井数据处理和解释的技术要求。
本标准适用于MRIL-C型、MRIL-C/TP型和MRIL-Prime型核磁共振测井数据的处理和解释。
2规范性引用文件
下列文件中的条款通过本标准的引用而成为本标准的条款。
凡是注日期的引用文件,其随后所有的修改单(不包括勘误的内容)或修订版均不适用于本标准,然而,鼓励根据本标准达成协议的各方研究是否可使用这些文件的最新版本。
凡是不注日期的引用文件,其最新版本适用于本标准。
SY/T5132测井原始资料质量要求
SY/T5360裸眼井单井测井数据处理流程
3解释软件
解释软件包括:
——express解释软件;
——DPP解释软件。
4测井资料质量检验
4.1依据SY/T5132规定对测井原始资料进行质量检验。
4.2检查对比原始测井资料与编辑回放测井资料的一致性。
5数据合并及深度校正
5.1数据合并
测井资料处理前,应将程序中所用到的测井数据转换成统一的数据格式,并合并为一个文件。
5.2深度校正
用核磁共振测井并测的自然伽马曲线进行深度校正。
6 MRIL-C型、MRIL-C/TP型核磁共振测井资料处理
6.1处理流程
MRIL-C型、MRIL-C/TP型资料处理流程如图1。
图1MRIL-C型、MRIL-C/TP型资料处理流程图
6.2回波处理(MRILPOST)
6.2.1回波处理流程如图2.
图2回波处理流程图
6.2.2对回波串进行反演拟合,得到T2分布、核磁共振有效孔隙度、地层束缚水孔隙度和可动流体孔隙度等。
6.2.3输入曲线主要包括:
——ECHO:
长等待时间原始回波串,单位为毫秒(ms);
——ECHOB:
短等待时间原始回波串,单位为毫秒(ms)。
6.2.4输入参数主要包括:
-STEP:
开关控制选择,表示暂停或继续;
-DEPTH:
深度信息;
-BIN:
用拟合回波串所用Bin的个数;
-ECHO:
计算T2分布的原始回波申序号、回波个数和回波间隔;
-MODE:
显示操作模式(浏览或记录);
-SCALE:
设置比例;
-FILTER:
设置低通滤波和平均值参数。
6.2.5输出曲线主要包括:
——MPHI:
核磁共振有效孔隙度,以百分数表示;
——MBVM:
可动流体孔隙度,以百分数表示;
——MBVI:
束缚水孔隙度,以百分数表示;
——MPERM:
渗透率,单位为毫达西(mD或10-3μm2);
——PPOR:
各弛豫组分的孔隙度,以百分数表示。
6.3差谱处理(MRDEA)
6.3.1差谱处理流程如图3.
图3差谱处理流程图
6.3.2利用长、短等待时间测量的T2谱的差值,定量评价孔隙中流体的类型及含量。
6.3.3油、气的弛豫时间参数(TlOIL,T20IL,T2GAS)应符合地区规律。
6.3.4与观测模式有关的参数(TE,TWL,TWS)应与实际观测模式相一致。
6.3.5输入曲线主要包括:
-ECHODIFF:
长、短等待时间T2谱的差值;
-MPHE:
用长等待时间计算的核磁共振有效孔隙度,以百分数表示;
-TEMP:
地层温度,单位为华氏度(oF)。
6.3.6输人参数主要包括:
-NECHO:
回波个数;
-T10IL:
油的纵向弛豫时间,单位为毫秒(ms);
-T20IL:
油的横向弛豫时间,单位为毫秒(ms);
-T2GAS:
气的横向弛豫时间,单位为毫秒(ms);
-TE:
回波间隔,单位为毫秒(ms);
-TWL:
长等待时间,单位为毫秒(ms);
-TWS:
短等待时间,单位为毫秒(ms)。
6.3.7输出曲线主要包括:
-BVGXO:
冲洗带含气孔隙度,以百分数表示;
-BVOXO:
冲洗带含油孔隙度,以百分数表示;
-BVWXO:
冲洗带含水孔隙度,以百分数表示;
-MPHEDEA:
烃校正后的核磁共振有效孔隙度,以百分数表示。
6.4综合流体分析(MRAX)
6.4.1综合流体分析处理流程如图4.
图4综合流体分析处理流程图
6.4.2综合流体分析是用核磁共振测井资料结合电阻率、中子、密度等测井资料计算地层的总孔隙度、有效孔隙度、束缚水孔隙度、自由流体孔隙度、含烃孔隙度和渗透率等参数.
6.4.3用于计算流体性质饱和度的参数(RW,RMT,M,A)选取按SY/T5360执行。
6.4.4输入曲线主要包括:
——CNC:
中子测井曲线,以百分数表示;
——FLAG:
求Swb标识符;
——MPHE:
有效孔隙度,以百分数表示;
——MBVI:
束缚水孔隙度,以百分数表示;
——POR:
外部输入孔隙度,以百分数表示;
——RT:
深电阻率,单位为欧姆米(Ω·m);
——SWBCORR:
所选的动态Swb,以百分数表示。
6.4.5输人参数主要包括:
——EDEP:
处理层段的底界深度,单位为米(m);
——SDEP:
处理层段的顶界深度,单位为米(m);
——FD:
流体密度,单位为克每立方厘米(g/cm3);
——RW:
地层水电阻率,单位为欧姆米(Ω·m);
——MPHECL:
粘土水孔隙度,以百分数表示;
——COEF,计算渗透率的系数;
——FLOOR:
束缚水饱和度的最小值,以小数表示;
——POROP:
计算总孔隙度选择符;
——SWBOPT:
计算束缚水饱和度选择符。
6.4.6输出曲线主要包括:
——CBVWE:
有效含水孔隙度,以百分数表示;
——CBVWT:
总含水孔隙度,以百分数表示;
——CLPOR:
粘土水孔隙度,以百分数表示;
——EPOR:
有效孔隙度,以百分数表示;
——MPHECO:
烃校正后的有效孔隙度,以百分数表示;
——PERM:
渗透率,单位为毫达西(mD或10-3μm2);
——SW:
含水饱和度,以百分数表示;
——SWB:
粘土水饱和度,以百分数表示;
——SWIA:
束缚水饱和度,以百分数表示。
7 MRIL-Prime型核磁共振测井资料处理
7.1处理流程
MRIL-Prime型资料处理流程如图5.
图5 MRIL-Prime型资料处理流程图
7.2岩石物理参数计算
7.2.1岩石物理参数计算处理解释流程如图6.
图6岩石物理参数计算处理解释流程图
7.2.2分组SPLIT_MCLS:
把含五组(a,b,c.d,e)回波串数据的文件分成若干个含三组数据
的文件,分组方式有abc,dec,adc和bec四种。
7.2.3排序RESEQ;将与测井观测模式相一致的ID码输入到参数中,对分组后文件中的回波序号(CACT)和回波串识别符(GRP)进行重新排序。
7.2.4回波拟合ECHE_STRIP。
7.2.4.1进行回波串多指数拟合.获得离散的T2分布及各区间孔隙度;计算长、短等待时间回波串的幅度差。
7.2.4.2观测模式编号(ACT-FLG)参数应与实际测井观测模式编号相一致。
7.2.4.3累加平均处理参数(RUNAVA,RUNAVB,RUNAVPR,RUNAVED)、回波拟合处理的第一个和最后一个序号参数(FECHOA,LECHOA,FECHOB,LECHOB,FECHOPR.LE-CHOPR)、平滑规则化参数(SMOOTHA,SMOOTHB,SMOOTHPR,SMOOTHED)和相位校正的回波序号的第一个和最后一个序号参数(FECHOPHA)应与观测模式提供的参数相一致。
7.2.4.4输入曲线主要包括:
-CACT:
观测模式中的一个周期内回波的序号;
-CECH:
回波个数;
-GRP:
回波串识别符;
-PWCO:
功率校正因子;
-RAMP:
原始回波幅度,以百分数表示;
-RPHA:
原始回波相位,以百分数表示。
7.2.4.5输人参数主要包括:
-ACT-FLG:
观测模式编号;
-RUNAVA,RUNAVB,RUNAVPR和RUNAVED:
A,B,C组和EDIF组的累加平均处理参数;-FECHOA,LECHOA,FECHOB,LECHOB,FECHOPR和LECHOPR:
A,B.C组第一个和最后一个做回波拟合处理的序号;
-SMOOTHA.SMOOTHB.SMOOTHPR和SMOOTHED:
A,B,C组和EDIF组平滑规则化参数;
-FECHOPHA,LECHOPHA,FECHOPHB和LECHOPHB:
A和B组第一个和最后一个做相位校正的回波序号。
7.2.4.6输出曲线主要包括:
-AMPA,AMPB和AMPC:
回波A,B和c组的幅度,以百分数表示;
-PHASA,PHASB和PHASC:
回波串A,B和C组的相位,单位为度(o);
-AVRA,AVRB和AVRC:
A,B和C组累加后的实部,以百分数表示;
-AVIA,AVIB和AVIC:
A,B和C组累加后的虚部,以百分数表示;
-REALCA:
A组经过相位、增益、功率、矿化度和温度校正以后的实部,以百分数表示;
-IMAGCA:
A组经过相位、增益、功率、矿化度和温度校正以后的虚部,以百分数表示;
-PO1A,P02A.P1A—P10A:
A组拟合的T2=lms,2ms,4ms,8ms,16ms,32ms,64ms,128ms,256ms,512ms,1024ms,2048ms的孔隙度,以百分数表示;
-PO1B,P02B,P1B--P10B:
B组拟合的T2=4ms,8ms,16ms,32ms,64ms,128ms,256ms,
512ms,1024ms,2048ms的孔隙度,以百分数表示;
-PIPR~P7PR:
C组拟合的T2=0.5ms,1ms,2ms,4ms,8ms,16ms,256ms的孔隙度,以百分数表示;
-EDIF:
回波幅度的差,以百分数表示。
7.2.5时一深转换Process_t2d。
7.2.5.1时一深转换是将时间域文件转为深度域文件。
7.2.5.2时一深转换采用英制时采样间距为0.5ft或0.25ft;采用公制时采样间距为0.1m。
7.2.5.3时一深转换应选择与其组数和深度单位相一致的配置文件。
7.2.6岩石物理参数计算T2_toolkit。
7.2.6.1 T2_toolkit可计算总孔隙度、有效孔隙度、束缚水含量和渗透率。
7.2.6.2滤波长度参数(FLTLEN)的取值应与采样间距相一致。
7.2.6.3计算束缚水的参数(T2CUTOFF,SBVIOEF)应根据岩心实验数据确定.如无实验数据,
砂岩层T2CUTOFF=33ms,SBVIOEF=0.062;碳酸盐岩层T2CUTOFF=90ms,SBVIOEF=0.009。
7.2.6.4 Coates渗透率模型的系数(CCOEF)应根据岩心实验数据确定,如无实验数据采用区域经验参数。
7.2.6.5计算理论烃的参数(FPRESS,FTEMP,OILVISC)应根据井深和地区规律确定。
7.2.6.6输入曲线应用回波拟合的输出曲线作为输入曲线。
7.2.6.7输人参数主要包括:
-FLTLEN:
滤波长度;
-T2SPLICE:
C组与A或B组合并的T2时间,单位为毫秒(ms);
--T2CUTOFF:
束缚流体T2截止值,单位为毫秒(ms);
-SBVICOEF:
计算束缚水含量(SBVI)的系数;
-TE:
回波间隔,单位为毫秒(ms);
-CCOEF:
Coates渗透率模型的系数;
-FFICUT:
可动流体截止值,以百分数表示。
7.2.6.8输出曲线主要包括。
-MCBW:
粘土水含量,以百分数表示;
-TASPEC:
A组T2分布,单位为毫秒(ms);
-MSIGTA:
用A组计算的总孔隙度,以百分数表示;
-MPHITA:
用A组计算的有效孔隙度,以百分数表示;
-CBVITA:
用合并后的A组计算的束缚水含量,以百分数表示;
-MPERM:
渗透率,单位为毫达西(mD或10-3μm2);
-EDSPEC:
回波串差EDIF的T2分布,以百分数表示;
-BOIL:
滤波后的含油孔隙度,以百分数表示;
-BGAS:
滤波后的含气孔隙度,以百分数表示;
-BWTR:
滤波后的含水孔隙度,以百分数表示;
-TDAMSIG:
极化与含氢指数校正后的总孔隙度,以百分数表示;
-TDAMPHI:
极化与含氢指数校正后的有效孔隙度,以百分数表示。
7.3差谱分析
7.3.1利用两组不同等待时间的双TW数据,做地层流体识别和油气定量评价。
7.3.2差谱分析处理流程如图7.
图7差谱分析处理流程图
7.3.3 弛豫时间搜索T2Tl_ EVENT。
7.3.3.1从双等待时间(TW)差异较大的区间中,对可能存在的流体进行T1和T2搜索。
7.3.3.2回波间隔(TE)、等待时间(TW)等参数应与测井同步生成的测井模式参数一致。
7.3.3.3流体T1,T2弛豫时间的第一个和最后一个区间序号应是双等待时间差异较大的区间。
7.3.3.4输入曲线主要包括:
-P1A~P12A:
第1~12区间的孔隙度(长TW),以百分数表示;
-P1B~P12B:
第1~12区间的孔隙度(短TW),以百分数表示;
-EDBST:
回波串差增强因子;
-EDIF:
长、短TW的回波串差,以百分数表示。
7.3.3.5输入参数主要包括:
-NUMECHO:
回波串差的回波个数;
_TE:
回波串差的回波间隔,单位为毫秒(ms);
-TWA:
长等待时间,单位为毫秒(ms);
-TWB:
短等待时间,单位为毫秒(ms);
-T1FBIN:
孔隙度最大差异对比的第一个区间序号;
-T1LBIN:
孔隙度最大差异对比的最后一个区间序号;
-FECHO:
T2搜索的第一个回波;
-LECHO:
T2搜索的最后一个回波。
7.3.3.6输出曲线主要包括:
-T1APP:
T1的估算值,单位为毫秒(ms);
-T1ERR:
T1搜索的误差最小值;
-T2X1:
单相搜索的T2.单位为毫秒(ms);
-PORX1:
单相搜索的孔隙度,以百分数表示;
-T2ERRX:
单相搜索的误差最小值;
-T2X:
双相搜索的短分量T2,单位为毫秒(ms);
-PORX:
短分量的孔隙度,以百分数表示;
-T2Y:
双相搜索的长分量T2,单位为毫秒(ms);
-PORY:
长分量的孔隙度,以百分数表示;
-T2ERRXY:
双相搜索的误差最小值。
7.3.4时间域分析TDA_COMP。
7.3.4.1通过对时间域的回波串幅度差进行处理,计算油、气、水的孔隙度。
7.3.4.2流体类型参数(FLUIDFLG)选择应与地层流体类型实际情况基本相符。
7.3.4.3有关流体T1.T2的参数应与T2T1_EVENT搜索的结果相一致。
7.3.4.4输入曲线(EDIF):
长、短等待时间的回波串幅度差,以百分数表示。
7·3.4.5输人参数主要包括:
-TWA:
长等待时间,单位为毫秒(ms);
-TWB:
短等待时间,单位为毫秒(ms);
-FLUIDFLG:
流体类型选择符;
-TIGAS:
气的纵向弛豫时间,单位为毫秒(ms);
-T2GAS:
气的横向弛豫时间,单位为毫秒(ms);
-TlOIL:
油的纵向弛豫时间,单位为毫秒(ms);
-T20IL:
油的横向弛豫时间,单位为毫秒(ms);
-T1WTR:
水的纵向弛豫时间,单位为毫秒(ms);
-T2WTR:
水的横向弛豫时间,单位为毫秒(ms).
7.3.4.6输出曲线主要包括:
-PHIG:
气的孔隙度,以百分数表示;
-PHIO:
油的孔隙度,以百分数表示;
-PHTW:
水的孔隙度,以百分数表示;
7.4移谱分析
7.4.l对两组不同回波间隔的T2分布进行扩散分析或扩散增强分析,来定量计算含水饱和度和油水相对渗透率。
7.4.2移谱分析处理流程如图8.
图8移谱分析处理流程图
7.4.3T2特征值选择T2_DIFSEL。
7.4.3.1计算出长、短TE分布谱的几何平均值、峰值、半峰值及用于扩散分析的横向弛豫时间。
7.4.3.2选择的长TE和短TE几何平均值的上限、下限值参数,应有利于反映不同性质流体的特征。
7.4.3.3输入曲线主要包括:
-MPHITA:
核磁共振有效孔隙度,以百分数表示;
-TBSPEC:
长TE回波串的T2分布,单位为毫秒(ms);
-TASPEC:
短TE回波串的T2分布,单位为毫秒(ms)。
7.4.3.4输入参数主要包括;
-T2GMUA:
短TE组计算几何平均值的上限,单位为毫秒(ms);
-T2GMLA:
短TE组计算几何平均值的下限,单位为毫秒(ms);
-T2GMUB:
长TE组计算几何平均值的上限,单位为毫秒(ms);
-T2GMLB:
长TE组计算几何平均值的下限,单位为毫秒(ms);
-T2SELFLG:
确定DIFAN中T2用哪组数据的选择标志,同时控制T2S和T2L。
7.4.3.5输出曲线主要包括:
-T2GMS:
短TE组的T2几何平均值,单位为毫秒(ms);
-T2DIFIS,短TE组孔隙度峰值对应的T2,单位为毫秒(ms);
-T2DIFS,短TE组孔隙度半峰值对应的T2.单位为毫秒(ms);
-T2S1用于扩散分析的T2值(短TE),单位为毫秒(ms);
-T2GML:
长TE组的T2几何平均值,单位为毫秒(ms);
-T2DIFIL:
长TE组孔隙度峰值对应的T2,单位为毫秒(ms);
-T2DIFL:
长TE组孔隙度半峰值对应的T2,单位为毫秒(ms);
-T2L:
用于扩散分析的T2值(长TE),单位为毫秒(ms)。
7.4.4扩散分析DIFAN。
7.4.4.1扩散分析是使用两组不同回波间隔的测井数据做扩散分析。
7.4.4.2计算流体饱和度的参数(T2MIN,T2HY,ROSF,RDDWHY)和计算相对渗透率参数(RELPFL,WVIS,OVIS)应符合地区流体性质的规律。
7.4.4.3输入曲线主要包括:
-MPHITA:
有效孔隙度(短TE),以百分数表示;
-MBVITA:
地层束缚水含量,以百分数表示;
-MPERM:
Coates模型渗透率,单位为毫达西(mD或10-3μm2);
-T2L:
用于扩散分析的长TE的T2值,单位为毫秒(ms);
-T2S:
用于扩散分析的短TE的T2值,单位为毫秒(ms)。
7.4.4.4输人参数主要包括:
-TEL:
长回波间隔,单位为毫秒(ms);
-TES:
短回波间隔,单位为毫秒(ms);
-T2HY:
烃的固有T2.单位为毫秒(ms);
-RDDWHY:
烃与水的扩散系数比值;
-T2SFACT:
水的最大T2值调节因子;
-RELPFL:
相对渗透率模型选择符;
-OVIS:
地层中油的粘度,单位为厘泊(cP)。
7.4.4.5输出曲线主要包括:
-DIFSW:
含水饱和度,以小数表示;
-DIFBVW:
束缚水孔隙度,以百分数表示;
-DIFKRO:
油相对渗透率;
-DIFKRW:
水相对渗透率。
7.5综合流体分析
7.5.1综合流体分析处理流程如图9.
图9综合流体分析处理流程图
7.5.2预处理PRE—MRIAN。
7.5.2.1利用核磁共振测井资料和其他测井资料计算视粘土水饱和度、总孔隙度、视地层水电阻率等参数,识别岩性。
7.5.2.2岩石骨架参数(RHOSAN,RHOLIM,RHOOL,DTMA)和计算流体的参RWREF,RMFREF,EXPHI,ARCHIE)的选取按SY/T5360执行。
7.5.2.3输人曲线主要包括:
-MSIGTA:
核磁总孔隙度,以百分数表示;
-MPHITA:
核磁有效孔隙度,以百分数表示;
-MBVITA:
束缚水孔隙度,以百分数表示;
-TDAMSIG:
TDA总孔隙,以百分数表示;
-TDAMPHI:
TDA有效孔隙度,以百分数表示;
-NPHI:
中子孔隙度,以百分数表示;
-RHOB:
体积密度,单位为克每立方厘米(g/cm3);
-PE:
光电吸收截面,单位为靶恩每电子(b/e);
-DT:
声波时差,单位为微秒每英尺(μs/ft);
-RT:
地层深探测电阻率,单位为欧姆米(Ω·m)。
7.5.2.4输人参数主要包括:
-TPORFL:
总孔隙度选择标志,应根据流体性质和井眼状况选取;
-RWREF:
地层水电阻率,单位为欧姆米(Ω·m);
-TWREF:
地层水电阻率的参考温度,单位为华氏度(oF);
-RMFREF:
钻井液滤液电阻率,单位为欧姆米(Ω·m);
-ARCHIE:
阿尔奇公式中的系数(a);
-EXPHI:
阿尔奇公式中的胶结指数(m);
-RHOFL:
流体密度,单位为克每立方厘米(g/cm3);
-DMA:
视骨架密度,单位为克每立方厘米(g/cm3);
-RHOSAN:
砂岩骨架密度,单位为克每立方厘米(g/cm3);
-RHOLIM:
灰岩骨架密度,单位为克每立方厘米(g/cm3);
-RHODOL:
白云岩骨架密度,单位为克每立方厘米(g/cm3);
-DTMA:
视声波骨架值,单位为微秒每英尺(μs/ft)。
7.5.2.5输出曲线主要包括:
-PHID:
密度孔隙度,以百分数表示;
-PHIN:
中子孔隙度,以百分数表示;
-PHIS,声波孔隙度,以百分数表示;
-PHIX:
中子密度交会孔隙度,以百分数表示;
-BN:
束缚水中子孔隙度,以百分数表示;
-DPND:
中子孔隙度与密度孔隙度的差,以百分数表示;
-DPTD:
声波孔隙度与密度孔隙度的差,以百分数表示;
-RPND:
中子孔隙度与密度孔隙度的比值;
-DPDM:
密度孔隙度与MRIL孔隙度的差,以百分数表示;
-TPOR:
总孔隙度,以百分数表示;
-RWA:
视地层水电阻率,单位为欧姆米(Ω·m);
-CWA:
视地层水电导率,单位为毫西门子每米(
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 核磁共振 测井 资料 处理 解释 规范