当前位置:首页 >> 能源/化工 >>

P核磁测井数据处理与参数文件命名规范(暂定稿)


QQ:270073230

川庆钻探工程公司测井公司(重庆)

测井技术发展中心方法研究所

1

MRIL-Prime型核磁共振测井数据处理 与参数文件命名规范(暂定稿)
一、MRIL数据处理解释流程
原始回波串 抽取、 分组 重新排序

ECHOSTRIP

(数据反演、拟合)

T2T1EVENT (T2、 值搜索) T1

PROCESS-T2D(时间域-深度域转换)

T2 TOOLKIT(岩石物理参数计算) 常 规 测 井 数 据

PRE_MRIN

T2_DIFSEL

SWB_MRIN

TDA_COMP

MRIAN

TDA/MRIAN

综合解释结果

二、SPLIT MCLS
作用:
把含5组数据的文件(时间域文件)分成几个含3组数据的文件。

运行方法:
Model Launcher→Numar→SPLIT mcls→输入m.cls文件 →输出文件名(abc、adc、dec)

说明:
DPP主窗底部目录必须是MRIL文件所在目录(见图1) 。 用CLS Editor查看该新文件中的CACT和GRP发生了变化。

2

图 1 DPP 工作站主界面 在 DPP 主界面双击 Model Launcher,出现以下界面,点击 Numar,然后根据需要,选取所要运 行的处理程序(见图 2) 。

图 2 Numar 处理主界面

3

三、RESEQ
作用:
把SPLIT mcls处理过的文件中的CACT和GRP进行排序。界面见图3。

图3 RESEQ排序程序模块运行界面(D9TWE2) 运行方法:
Model Launcher→Numar→RESEO 输入SPLIT处理过的CLS文件; 点击EDIT→CURVE MAP,确定输入、输出曲线名(重定向); 点击OPTION→create parameter database,生成新的参数文件,或File→Open parameter database,打开已存在的参数文件; 点击OPTION→RUN PARAMETER EDITOR,对参数文件进行编辑,从cact reseq db.txt中找到文件 编号,放在参数ACT-SEQ中; 点击SAVE SESSION FILE,生成各种文件集合,以备再处理之用,最后点击RUN,运行此程序

说明:
1、该程序只能运行一次,若出错,必须从SPLIT mcls开始重做。
4

2、参数文件是/desktop/dpp/etc/cact reseq db.txt 3、若有SESSION文件,当需要重新运行时,可调出再运行。 用CLS Editor会发现该新文件中CACT和GRP发生变化,CACT已重新排序。

四、ECHO STRIP
作用:
对核磁测量信号进行回波串拟合,获得各bin的孔隙度及T2拟合线;对TDA方法还输出差谱EDIF 。程序运行界面见图4。

图4 ECHO STRIP程序运行界面

运行方法:
Model Launcher→Numar→ECHO_STRIP 输入上一步处理过的CLS文件; 点击EDIT→CURVE MAP,确定输入、输出曲线名(重定向); 点击OPTION→create parameter database,生成新的参数文件,或File→Open parameter database,打开已存在的参数文件; 点击OPTION→RUN PARAMETER EDITOR,对参数文件进行编辑; 点击SAVE SESSION FILE,生成各种文件集合,以备再处理之用,最后点击RUN,运行此程序

参数设置:
ACT-FLG——模式编号(从desktop/dpp/etc/dpp act db.txt中查找) 1、做平均处理的参数 RUNAVA——对A组平均处理(隐含16,范围1-80)
5

RUNAVB——对B组平均处理(隐含16,范围1-80) RUNAVPR——对PRO6组平均处理(隐含4,范围1-80) RUNAVED——对差谱EDIP平均处理(隐含4,范围1-8) 2、对A组做相位校正的参数 PHSCORA——对A组做相对校正(=1做校正,=0不做) PTHRESA——A组孔隙度门槛(隐含0,范围1-10) NLVLPHA——相位校正的累计水平(隐含5,范围1-12) FECHOPHA——A组第一个做组相校正的回波序号(隐号2,范围1-5) LECHOPHA——A组最后一个做相位校正的回波序号(隐含9,范围3-1000) 3、对B组做相位校正的参数 PHSCORB——对B组做相位校正(=1做校正,=0不做) PTHRESB——B组孔隙度门槛(隐含0,范围1-10) NLVLPHB——B组相位校正的累计水平(隐含5,范围1-12) FECHOPHB——B组第一个做组相校正的回波序号(隐号2,范围1-5) LECHOPHB——B组最后一个做相位校正的回波序号(隐含9,范围3-1000) 4、BAND鉴别参数(只在有坏BAND时使用) DBANDO——鉴别BAND0(隐含0,0=不做,1=做) DBAND1——鉴别BAND1(隐含0,0=不做,1=做) DBAND2——鉴别BAND2(隐含0,0=不做,1=做) DBAND3——鉴别BAND3(隐含0,0=不做,1=做) DBAND4——鉴别BAND4(隐含0,0=不做,1=做) 5、MAP处理参数 FECHOA——A组第一个做处理的回波序号(隐含2,范围1-51) LECHOA——A组最后一个做处理的回波序号(隐含400,范围10-1000) FECHOB——B组第一个做处理的回波序号(隐含2,范围1-51) LECHOB——B组最后一个做处理的回波序号(隐含400,范围10-1000) NECHOA——A组回波个数 NECHOB——B组回波个数 ERFLG——提高分辨率标志(=1进入提高分辨率处理,=0不做)

TE 0.9ms 1.2ms 2.4ms 3.6ms 4.8ms 6.0ms 输入曲线
AMPN:前3个回波的幅度 B1:井下测量的原始B1。

MRIL—P相位校正的回波序号 2—12 2—9 1—6 1—5 1—4 1—3

6

B1MO:温度校正后的B1。 CACT:观测模式中的一个周期内回波的序号。 CECH:回波个数。 DCOF:由回波串计算的直流偏置(未经刻度) 。 DLOD:井下张力。 E1CO:第一个回波校正因子。 E2CO:第二个回波校正因子。 GAIN:增益。 GAMM:自然伽马(未经刻度) 。 GRP:回波串识别符。 OFFS:直流偏置均值(未经刻度) 。 NOIS:直流偏置的标准偏差(未经刻度) 。 IEDC:回波间直流偏置的均值(未经刻度) 。 PWCO:功率校正因子。 RAMP:原始回波幅度(经过所有校正,未经累加) 。 SACO:矿化度校正因子。 SEQN:观测模式序号。 STAT:数据传输状态。 TMCO:温度校正因子。 TMP1:发射器温度。 TMP3:探头温度(单位:℃) 。 TTIM:地面系统记录的回波串传输时间。

输出曲线
1、回波串A处理后的输出曲线 DEPTHA:回波串A的参考深度(累加窗口的中心) AMPA:回波串A的幅度(经过刻度、增益、功率、矿化度、温度校正、没加累加) PHASA:回波串A的相位(经过刻度、增益、功率、矿化度、温度校正、没加累加) REALA:回波信号实部(直角坐标系) ,经过刻度,没做校正,没做累加 IMAGA:回波信号虚部(直角坐标系) ,经过刻度,没做校正,没做累加 AVIA:累加后的实部 AVIA:累加后的虚部 REALCA:经过相位校正以后的ACRA,并做增益、功率、矿化度、温度校正 IMAGCA:经过相位校正以后的AVIA,并做增益、功率、矿化度、温度校正 ECHOA:重复REALCA(仅用于显示) 2、回波串A拟合处理后的质量控制曲线 PHANA:A组的平均相位角 PHNOA:噪声(虚部)的标准偏差,经累加和相位校正 PHERA:噪声(虚部)的均值,经累加和相位校正 SNRA:回波串的信噪比 CHIA:拟合误差的标准偏差(拟合质量) 3、回波串A反演拟合后的输出曲线
7

P01A:T2=1.0ms的孔隙度幅度 P02A:T2=2.0ms的孔隙度幅度 P1A:T2=4.0ms的孔隙度幅度 P2A:T2=8.0ms的孔隙度幅度 P3A:T2=16.0ms的孔隙度幅度 P4A:T2=32.0ms的孔隙度幅度 P5A:T2=64.0ms的孔隙度幅度 P6A:T2=128.0ms的孔隙度幅度 P7A:T2=256.0ms的孔隙度幅度 P8A:T2=512.0ms的孔隙度幅度 P9A:T2=1024.0ms的孔隙度幅度 P10A:T2=2048.0ms的孔隙度幅度 PBINA:孔隙度幅度向量 4、回波串B处理后的输出曲线 DEPTHB:回波串B的参考深度(累加窗口的中心) AMPB:回波串B的幅度(经过刻度、增益、功率、矿化度、温度校正,没做累加) PHASB:回波串B的相位(经过刻度、增益、功率、矿化度、温度校正,没做累加) REALB:回波信号实部(直角坐标系) ,经过刻度,没做校正,没做累加 IMAGB:回波信号虚部(直角坐标系) ,经过刻度,没做校正,没做累加 AVRB:累加后的实部 AVIB:累加后的虚部 REALCB:经过相位校正以后的AVRB,并做增益、功率、矿化度、温度校正 IMAGCB:经过相位校正以后的AVIB,并做增益、功率、矿化度、温度校正 ECHOB:重复REALCB(仅用于显示) 5、回波串B拟合处理后的质量控制曲线 PHANB:B组的平均相位角 PHNOB:噪声(虚部)的标准偏差,经累加和相位校正 PHERB:噪声(虚部)的均值,经累加和相位校正 SNRB:回波串的信噪比 CHIB:拟合误差的标准偏差(拟合质量) 6、回波串B反演拟合后的输出曲线 P01B:T2=1.0ms的孔隙度幅度 P02B:T2=2.0ms的孔隙度幅度 P1B:T2=4.0ms的孔隙度幅度 P2B:T2=8.0ms的孔隙度幅度 P3B:T2=16.0ms的孔隙度幅度 P4B:T2=32.0ms的孔隙度幅度 P5B:T2=64.0ms的孔隙度幅度 P6B:T2=128.0ms的孔隙度幅度 P7B:T2=256.0ms的孔隙度幅度 P8B:T2=512.0ms的孔隙度幅度
8

P9B:T2=1024.0ms的孔隙度幅度 P10B:T2=2048.0ms的孔隙度幅度 PBINB:孔隙度幅度向量 7、回波串C(RP06)处理后的输出曲线 DEPTHC:回波串C的参考深度(累加窗口的中心) AMPC:回波串C的幅度(经过刻度、增益、功率、矿化度、温度校正,没做累加) PHASC:回波串C的相位(经过刻度、增益、功率、矿化度、温度校正,没做累加) REALC:回波信号实部(直角坐标系) ,经过刻度,没做校正,没做累加 IMAGC:回波信号虚部(直角坐标系) ,经过刻度,没做校正,没做累加 AVRC:累加后的实部 AVIC:累加后的虚部 REALCC:经过功率、矿化度、温度校正后的AVRC 8、回波串C拟合处理后的质量控制曲线 SNRC:回波串的信噪比 CHIC:拟合误差的标准偏差(拟合质量)

9、回波串C反演拟合后的输出曲线
P1PR:T2=0.5ms的孔隙幅度 P2PR:T2=1.0ms的孔隙幅度 P3PR:T2=2.0ms的孔隙幅度 P4PR:T2=4.0ms的孔隙幅度 P5PR:T2=8.0ms的孔隙幅度 P6PR:T2=16.0ms的孔隙幅度 P7PR:T2=256.0ms的孔隙幅度 PBINPR:孔隙度幅度向量 10、双TW模式回波串的差(用于TDA分析) EDIF:回波幅度的差(A组-B组) 11、回波串的差处理之后的质量控制曲线 CHIED:拟合误差的标准偏差(拟合质量) 12、回波串差反演拟合后的输出曲线 P01ED:T2=1.0ms的孔隙度幅度 P02ED:T2=2.0ms的孔隙度幅度 P1ED:T2=4.0ms的孔隙度幅度 P2ED:T2=8.0ms的孔隙度幅度 P3ED:T2=16.0ms的孔隙度幅度 P4ED:T2=32.0ms的孔隙度幅度 P5ED:T2=64.0ms的孔隙度幅度 P6ED:T2=128.0ms的孔隙度幅度 P7ED:T2=256.0ms的孔隙度幅度 P8ED:T2=512.0ms的孔隙度幅度 P9ED:T2=512.0ms的孔隙度幅度 P10ED:T2=512.0ms的孔隙度幅度
9

PBINED:孔隙度幅度向量 13、增益(用于质量控制) GAINA:A组的增益(经过累加) GAINB:B组的增益(经过累加) GAINC:C组的增益(经过累加) 14、反演拟合后的质量控制曲线 YFITA:回波串A的理论值 YFITB:回波串B的理论值 YFITPR:回波串C的理论值 YFITATR:回波串ATR的理论值 YFITBTR:回波串BTR的理论值 15、频带相位曲线 PHANA0:A组0频带相位角 PHANA1:A组1频带相位角 PHANA2:A组2频带相位角 PHANA3:A组3频带相位角 PHANB0:B组0频带相位角 PHANB1:B组1频带相位角 PHANB2:B组2频带相位角 PHANB3:B组3频带相位角 16、回波串的差增强因子 EDBST:用于增强回波串的差 17、深度对比曲线 MPHIX:核磁共振总孔隙度(仅用于深度对比)

说明
1、质量控制 ERFLG为提高分辨率标志,当其等于0时,各个做平均处理的参数RUNAVA、RUNAVB、RUNAVPR对不 同模式有最小值,不能下调,且不计算纵向分辨率RAAPER。当ERFLG为1时,启动SETOR方法,参 数RUNAVA、RUNAVB、RUNAVPR可以下调,与NLVLPH一起变化,NLVLPH按2倍增加,目的是:CHI<2 ,PHNO≈1.0,PHER≈0,SNR>5,RAAPER较合理。 2、ECHO STRIP运行完后,有几个临时文件形成,即“MWS……” ,占很大空间,应删掉。 3、经常使用的处理参数:RUNAVA、RUNAVB、RUNAVPR、NLVLPHA、FECHOPHA、LECHOPHA、NLVLPHB 、FECHOPHB、LECHOPHB、FECHOA、LECHOA、FECHOB、LECHOB、NECHOA、NECHOB、ERFLG。其余参 数一般不用动。

五、PROCESS T2D
作用:
将核磁时间域文件转为深度域文件,以备下步处理使用。

运行方法:
Model Launcher→Numar→PROCESS T2D
10

打开输入文件; 给定输出文件; 输入采样间距以及相应的输入配置文件; 判断是否选择延迟(Y或N) ,若选择N,则根据提示输入延迟深度;

注意:
1、先查原文件中的TDEP是公制还是英制,再根据公英制以及文件中有几组数据选择配置文件, 该程序不能进行公英制转制。具体如下: “2组数据且英制” “2组数据且公制” “3组数据且英制” “3组数据且公制” m.m2d /…/…/config 2 /…/…/config 2 /…/…/config 3 /…/…/config 3 grp ft.m2d grp m.m2d grp ft.m2d grp m.m2d 3 grp

2、用CLS Editor可看出时间文件转为了深度文件。我们一般选择:/…/…/config

六、T2 TOOLKIT
作用:
将各区间孔隙度滤波,计算MSIGTA、MPHITA、MBVITA和渗透率,输出T2分布,计算理论烃性质 。

运行:
Model Launcher→Numar→T2 TOOLKIT 输入上一步处理过的CLS文件; 点击EDIT→CURVE MAP,确定输入、输出曲线名(重定向); 点击OPTION→create parameter database,生成新的参数文件,或File→Open parameter database,打开已存在的参数文件; 点击OPTION→RUN PARAMETER EDITOR,对参数文件进行编辑; 点击SAVE SESSION FILE,生成各种文件集合,以备再处理之用,最后点击RUN,运行此程序

参数设置:
ACT-FLG——模式编号(从desktop/dpp/etc/dpp act db.txt中查找) 。 OSFLG——是否老类型的谱(=0不是,=1是) 。 1、关于滤波的参数 FILTFLG——是否滤波(=0不是,=1是) 。 FLTTYP——滤波类型(=0为Hanning,=1为Triangular,=2为Block) 。 FLTLEN——滤波长度(隐度5,范围3-21) 当采样间距为0.5ft,则令FLTLEN=5; 当采样间距为0.25ft,则令FLTLEN=9; 当采样间距为0.1m,则令FLTLEN=7; T2SPLICE——用于把PR06与A或B组合并的的T2时间。 2、关于计算BVI的参数 T2CUTOFF——束缚流体T2截止值(隐含33) 。 若无岩心数据,则砂岩用33,碳酸岩用90。
11

SBVICOEF——SBVI系数(隐含0.062) 。 若无岩心数据,则砂岩用0.062,碳酸岩用0.009。 3、关于计算渗透率的参数 CCOEF——Coates模型的渗透率系数(隐含10) 。 MUDTYP——泥浆类型(=0水基,=1油基) CTHRES——渗透率的最小孔隙度门槛(隐含0.05) 。 SBCOEF——Ben Swanson模型的渗透率系数(×10 ) 。 SBEXP——Ben Swanson模型的渗透率指数。 FFICUT——可动流体截止值。 T2STRTSB——PERMSB的低T2 Bin时间。 T2ENDSB——PERMSB的高T2 Bin时间。 CUMSUMTH——PERMSB的累计门槛。 SBFITTYP——PERMSB的拟合模型。 0=为Cubic Spline 1=为2nd order poly. 4、关于理论烃性质的参数 HCFLG——计算理论烃性质的标志。 (=0不计算,=1计算) 。 TOOLTYP——核磁仪器类型(=0为C型6",=1为C型4.5",=2为P型) 。 FRQFLG——仪器频率(=0低频,=1高频) 。 TE——回波间距(单位:msec) 。 FPRESS——地层压力(单位:PSI) 。 FTEMPF——地层温度(单位 5、自然伽马刻度参数 RGRGAIN——自然伽马增益系数。 RGROFFS——自然伽马偏移量。 ) 。 OILVISC——油的粘度(单位:cp) 。
-7

输出曲线
1、C组(PRO6)的输出曲线 P1PRF:滤波后的T2=0.5ms的孔隙度幅度 P2PRF:滤波后的T2=1.0ms的孔隙度幅度 P3PRF:滤波后的T2=2.0ms的孔隙度幅度 P4PRF:滤波后的T2=4.0ms的孔隙度幅度 P5PRF:滤波后的T2=8.0ms的孔隙度幅度 P6PRF:滤波后的T2=16ms的孔隙度幅度 P7PRF:滤波后的T2=256ms的孔隙度幅度 MCBW:泥质束缚水孔隙度(PRO6:P0.5ms+P1ms+P2ms) PRSPEC:PRO6的T2孔隙度分布 2、B组(短TW或长TE)的输出曲线 PO1BF:滤波后的T2=1.0ms的孔隙度幅度 PO2BF:滤波后的T2=2.0ms的孔隙度幅度 P1BF:滤波后的T2=4.0ms的孔隙度幅度
12

P2BF:滤波后的T2=8.0ms的孔隙度幅度 P3BF:滤波后的T2=16ms的孔隙度幅度 P4BF:滤波后的T2=32ms的孔隙度幅度 P5BF:滤波后的T2=64ms的孔隙度幅度 P6BF:滤波后的T2=128ms的孔隙度幅度 P7BF:滤波后的T2=256ms的孔隙度幅度 P8BF:滤波后的T2=512ms的孔隙度幅度 P9BF:滤波后的T2=1024ms的孔隙度幅度 P10BF:滤波后的T2=2048ms的孔隙度幅度 MSIGTB:总孔隙度,即C组(P0.5+P1+P2+P4)+B组(P8+……+P2048) MPHITB:有效孔隙度(P4PR+P2B+……+P1OB,即4—2048ms) TBSPEC:总孔隙度T2分布,即C组与B组合并(0.5-2048ms) 3、A组(长TW或短TE)的输出曲线 PO1AF:滤波后的T2=1.0ms的孔隙度幅度 PO2AF:滤波后的T2=2.0ms的孔隙度幅度 P1AF:滤波后的T2=4.0ms的孔隙度幅度 P2AF:滤波后的T2=8.0ms的孔隙度幅度 P3AF:滤波后的T2=16ms的孔隙度幅度 P4AF:滤波后的T2=32ms的孔隙度幅度 P5AF:滤波后的T2=64ms的孔隙度幅度 P6AF:滤波后的T2=128ms的孔隙度幅度 P7AF:滤波后的T2=256ms的孔隙度幅度 P8AF:滤波后的T2=512ms的孔隙度幅度 P9AF:滤波后的T2=1024ms的孔隙度幅度 P10AF:滤波后的T2=2048ms的孔隙度幅度 MSIGTA:总孔隙度,即C组(P0.5+P1+P2+P4)+A组(P8+……+P2048) MPHITA:有效孔隙度(P4PR+P2A+……+P1OA,即4—2048ms) TASPEC:总孔隙度T2分布,即C组与A组合并(0.5-2048ms) CBVITA:截止值BVI,用A组合并的谱计算 SBVITA:谱加权BVI,用A组合并的谱计算 MBVITA:最大BVI(CBVITA与SBVITA中的最大者) MPERM:渗透率估计值(Coates模型,用MPHITA) PERMSB:渗透率估计值(Swanson方法) 4、回波串差的输出曲线 PO1EDF:滤波后的T2=1.0ms的孔隙度幅度 PO2EDF:滤波后的T2=2.0ms的孔隙度幅度 P1EDF:滤波后的T2=4.0ms的孔隙度幅度 P2EDF:滤波后的T2=8.0ms的孔隙度幅度 P3EDF:滤波后的T2=16ms的孔隙度幅度 P4EDF:滤波后的T2=32ms的孔隙度幅度 P5EDF:滤波后的T2=64ms的孔隙度幅度
13

P6EDF:滤波后的T2=128ms的孔隙度幅度 P7EDF:滤波后的T2=256ms的孔隙度幅度 P8EDF:滤波后的T2=512ms的孔隙度幅度 P9EDF:滤波后的T2=1024ms的孔隙度幅度 P10EDF:滤波后的T2=2048ms的孔隙度幅度 DSPHIED:回波串差EDIF的总孔隙度 EDSPEC:回波串差EDIF的T2分布 5、油气NMR特性理论值曲线 DOI:MRIL探测深度 MGRAD:MRIL磁场梯度 RHOGAS:甲烷密度估算 T1GAS:甲烷的T1 HIGAS:甲烷的含氢指数 T2GAS:甲烷的T2 DIGAS:甲烷扩散系数(X10-5) T1OIL:油的T1 T2OIL:油的T2 DOOIL:油的扩散系数(X10-5) 6、EDIF时间域分析(TDF)后的曲线 PHIGAF:气的视孔隙体积(长TW) PHIOAF:油的视孔隙体积(长TW) PHIWAF:水的视孔隙体积(长TW) PHIGF:气相的孔隙体积 PHIOF:油相的孔隙体积 PHIWF:水相的孔隙体积 PHIFPL:TDA完全极化的液相孔隙度 BWTR:TDA水的孔隙体积 BOIL:TDA油的孔隙体积 BGAS:TDA气的孔隙体积 TDAMSIG:极化与含氢指数校正后的总孔隙度 TDAMPHI:极化与含氢指数校正后的有效孔隙度 TDAPERM:TDAA渗透率(Coates模型用TDAMPHI) MGR:刻度与滤波后的伽马射线

七、核磁流体分析MRIAN
(一) 、PRE MRIAN(优化总孔隙度)
作用:
结合常规测井中子、密度、电阻率、自然伽马等曲线,计算出几种孔隙度,同时利用核磁资料 计算束缚水饱和度SWBMRI。
14

运行方法:
Model Launcher→Numar→PRE-MRIAN 输入上一步处理过的CLS文件; 点击EDIT→CURVE MAP,确定输入、输出曲线名(重定向); 点击OPTION→create parameter database,生成新的参数文件,或File→Open parameter database,打开已存在的参数文件; 点击OPTION→RUN PARAMETER EDITOR,对参数文件进行编辑; 点击SAVE SESSION FILE,生成各种文件集合,以备再处理之用,最后点击RUN,运行此程序 注意:TMP3(磁铁温度)单位应为DEGC,即℃; NPHI必须为百分数,以及是使用岩性刻度,什么仪器类型。

参数设置:
1、SET UP参数 0=HES CNL 1=HES SNL 2=HES DSN 3=HES DSEN 4=SWS NPHI(CNT-H) 5=SWS NPHI(CNT-G) 6=SWS NPHI(CNT-A) 7=SWS TNPH/NPOR CURVE 8=SWS APS APXX CURCES 9=SWS APS FPXX CURVES 10=SPERRY 6.75in CNP 11=SPERRY 8in CNP 12=SPERRY 4.75in CTN 13=SPERRY 6.75in CTN 14=SPERRY 8in CTN 15=ANADRILL6.5in CDN 16=ANADRILL 8in CDN 17=BPB CNT 18=BAKER-ATLAS 2420 CN 19=BAKER-ATLAS 2435 CN INNEUTFL——中子输入骨架值选项(0=灰岩,1=砂岩,2=白云岩) OUTNEUTFL——中子输出骨架值选项(0=灰岩,1=砂岩,2=白云岩) TPORFL——总孔隙度选项 0=中子密度交会 1=密度 2=中子 3=声波 NTOOLFL——中子仪器类型

15

4=声波-密度平均 5=声波中子交会 6=MSIGTA 7=TDAMSIG 8=DMRP 9=外来孔隙度 DPORFL——密度输入骨架(0=灰岩2.71,1=砂岩2.65,2=白云岩2.87) 2、流体/骨架 参数 SURFT——平均地表温度(单位:DEGF) BHT——井底温度(单位:DEGF) TD——井深(单位:FT) PRESSG——地层压力梯度(单位:PSI/ft) PRESSF——视地层压力(单位:PSI) PRESSD——视地层压力的深度(单位:ft) RWREF——地层水电阻率 TWREF——地层水电阻率的参考温度(单位:DEGF) RMFREF——泥浆滤液电阻率 TMFREF——泥浆滤液电阻率的参考温度(单位:DEGF) ARCHIE——阿尔奇常数(a) EXPHI——胶结指数 RHOFL——流体密度(单位:g/cm3) DMA——视骨架密度(单位:g/cm3) RHOSAN——砂岩骨架密度(单位:g/cm3) RHOLIM——灰岩骨架密度(单位:g/cm3) RHODOL——白云岩骨架密度(单位:g/cm3) DTMA——视声波骨架值(单位:μs/ft) SONICFL——声波方程选择(0=RAIGA-CLEMENCEAU,1=CHAPMAN MODIFIED)

输入曲线:
MSIGTA:MRIL的总孔隙度(单位:%) MPHITA:MRIL的有效孔隙度(单位:%) MBVITA:MRIL的束缚体积(单位:%) TDAMSIG:TDA的总孔隙度(单位:%) TDAMPHI:TDA的有效孔隙度(单位:%) DMRP:DMR的总孔隙度(单位:%) DMREPHI:DMR的有效孔隙度(单位:%) NPHI:中子孔隙度(单位:%) RHOB:体积密度 PE:光电吸收截面 DPOR:密度孔隙度(单位:%) DT:声波时差 RT:真电阻率
16

TPOREXT:外来总孔隙度(单位:%)

输出曲线:
FTEMP:地层温度 RWAT:地层温度下的地层水电阻率 RMFAT:地层温度下的泥浆滤液电阻率 SALW:地层水矿化度(单位:Kppm) SALMF:泥浆滤液矿化度(单位:Kppm) RHOW:地层水密度 RHOMF:泥浆滤液密度 RHOBR:从DPOR计算的骨架密度 PHID:计算的密度孔隙度 PHIN:计算的中子孔隙度 PHIS:计算的声波孔隙度 PHIX:计算的中子密度交会孔隙度 RHOMA:视骨架密度 UMA:视骨架Pe值 TPOR:总孔隙度 RWA:视地层水电阻率 CWA:视地层水电导率 SWBMRI:MPHITA确定的束缚水饱和度 BN:束缚水中子值 DPND:中子与密度孔隙度的差 DPTD:声波与密度孔隙度的差 RPND:中子与密度孔隙度的比值 DPDM:密度孔隙度与MRIL孔隙度MPHITA的差 MFFI:总的自由流体指数 VQTZ:视骨架石英含量 VCAL:视骨架方解石含量 VDOL:视骨架白云质含量

(三) 、SWB-MRIAN(优化束缚水饱和度)
作用:
利用除了核磁资料之外的现有曲线计算所有可能的束缚水饱和度,并且选择最优的束缚水饱和 度作为最终的束缚水饱和度SWB。

运行:
Model Launcher→Numar→SWB-MRIAN 输入上一步处理过的CLS文件; 点击EDIT→CURVE MAP,确定输入、输出曲线名(重定向); 点击OPTION→create parameter database,生成新的参数文件,或File→Open parameter database,打开已存在的参数文件;
17

点击OPTION→RUN PARAMETER EDITOR,对参数文件进行编辑; 点击SAVE SESSION FILE,生成各种文件集合,以备再处理之用,最后点击RUN,运行此程序

参数设置
1、流体参数 SURFT——平均地表温度(单位:DEGF) BHT——井底温度(单位:DEGF) TD——井深(单位:ft) PRESSG——地层压力梯度(单位:PSI/ft) PRESSF——视地层压力(单位:PSI) PRESSD——视地层压力的深度(单位:ft) 2、粘土参数 GRB——泥岩的GR GRF——纯岩性的GR BNB——泥岩的BN BNF——纯泥岩的BN STEIBER——STEIBER的相关系数 DPNDB——泥岩的DPND DPNDF——纯岩性的DPND RHOMAB——泥岩的RHOMA RHOMAF——纯岩性的RHOMA DPTDB——泥岩的DPTD DPTDF——纯岩性的DPTD RTB——泥岩的RT RTF——纯岩性的RT EPHLSLP——中子密度有效孔隙度的斜率(=-1表示不用EPHI) RPDDB——泥岩的RPND RPNDF——纯岩性的RPND XB——泥岩的X XF——纯岩性的X MRIFLG——SWBMRI的计算标志(=0不计算,=1计算) MSIGFLG——SWBMSIG的计算标志(=0不计算,=1计算) SPFLG——SWBSP的计算标志(=0不计算,=1计算) DMRFLG——DMRSWB的计算标志(=0不计算,=1计算) HLWGT——HODGES-LEHMAN加权方法 0=平均 1=最小的SWB 2=最大的SWB SWBFLG——束缚水饱和度选择标志 0=HODGES-LEHMAN 1=最小值 2=最大值
18

3=SWBMRI 4=SWBMSIG 5=SWBGR 6=SWBNEUT 7=SWBDPND 8=SWBRHOMA 9=SWBDPTD 10=SWBRT 11=SWBSP 12=SWBEPHI 13=SWBRPND 14=SWBX 15=DMRSWB TSWBFL——TPOR与SWB标志 0=TPOR 1=SWB 2=MPHI WCLAY——RPHICLAY的W指数 PHIMAX——预计的总孔隙度最大值 输入曲线 PHID:计算的密度孔隙度 PHIN:计算的中子孔隙度 PHIS:计算的声波孔隙度 PHIX:计算的中子密度交会孔隙度 RHOMA:视骨架密度 TPOR:总孔隙度 RWA:视地层水电阻率 CWA:视地层水电导率 SWBMRI:MPHITA确定的束缚水饱和度 BN:束缚水中子值 DPND:中子与密度孔隙度的差 DPTD:声波与密度孔隙度的差 RPND:中子与密度孔隙度的比值 输出曲线 SWBMRI:MPHITA计算的泥质束缚水饱和度 SWBMSIG:MSIGA计算的泥质束缚水饱和度 SWBGR:GR计算的泥质束缚水饱和度 SWBNEUT:BN计算的泥质束缚水饱和度 SWBDPND:DPND计算的泥质束缚水饱和度 SWBRHOMA:RHOMA计算的泥质束缚水饱和度 SWBDPTD:DPTD计算的泥质束缚水饱和度
19

SWBRT:RT计算的泥质束缚水饱和度 SWBSP:SP计算的泥质束缚水饱和度 SWBMIN:最小计算的泥质束缚水饱和度 SWBMAX:最大计算的泥质束缚水饱和度 SWBMED:计算的泥质束缚水饱和度均值 SWB:用户选择的泥质束缚水饱和度 CCW:等效泥质束缚水电导率 QV:阳离子交换容量 WW:含水饱和度100%时的视W变量(WQ在SWT=1的边界) WI:束缚水条件下视W变量(WQ在SWT=SWTIRR的边界) WQ:双水模型下的W系数 PHICLAY:视泥质水孔隙度 RPHICLAY:视泥质水孔隙度与总孔隙度的比值 RWE:等效地层水电阻率(由SP计算) CWE:等效地层水电导率(由SP计算)

(三) 、MRIAN(核磁流体分析)
作用:
使用双水模型,综合核磁资料和常规测井资料,提供总孔隙度和有效孔隙度、总含水体积和有 效含水体积、不可动水体积、自由流体体积、总含水饱和度和有效含水饱和度、粘土水饱和度、渗 透率。

运行方法:
Model Launcher→Numar→MRIAN 输入上一步处理过的CLS文件; 点击EDIT→CURVE MAP,确定输入、输出曲线名(重定向); 点击OPTION→create parameter database,生成新的参数文件,或File→Open parameter database,打开已存在的参数文件; 点击OPTION→RUN PARAMETER EDITOR,对参数文件进行编辑; 点击SAVE SESSION FILE,生成各种文件集合,以备再处理之用,最后点击RUN,运行此程序

参数设置:
1、粘土/渗透率参数 RPINDX——RPHICALY-SWB交会图指数 TSWBFL——TPOR-SWB校正标志 0=TPOR 1=SWB 2=MPHI CSS——SW与SXO关系系数 FIXW——固定W的逻辑标志(=0不固定,=1固定) PCOEF——渗透率系数(Coates公式系数) PTHRES——渗透率门槛
20

POREXP——孔隙度指数( Coates公式系数) BDFEXP——BVI/FFI指数( Coates公式系数) PHIMAX——预计的最大孔隙度 2、流体系数 SURFT——地表温度(单位:DEGF) BHT——井底温度(单位:DEGF) TD——井深(单位:ft) PRESSG——地层压力梯度(单位:PSI/ft) PRESSF——视地层压力(单位PSI) PRESSD——视地层压力的深度(单位:ft) RWREF——地层水电阻率 TWREF——地层水电阻率的参考温度(单位:DEFG) RMFREF——泥浆滤液电阻率 TMFREF——泥浆滤液电阻率的参考温度(单位:DEFG) WFLAG——参数W的逻辑标志 0=原始MRIAN 1=线性SWT 2=非引性SWT EXPHI——胶结指数 EXSAT——饱和度指数 WOFF——使用SWT与W线性关系的偏移 WGAIN——使用SWT与W线性关系的增益 3、稠油参数 HOFLAG——稠油逻辑处理标志(=0不能,=1能) TEB——B组回波间距 CFRQ——仪器中心频率 TOOLTYP——核磁仪器类型 0=MRLL—6″探头 1=MRLL—4.5″探头 2=MRLL Prime 6″探头 3=MRLL Prime 4.875″探头 GFACT——磁场梯度的倍数

输入曲线
SWB:束缚水饱和度 SWBMRI:来自MPHITA的束缚水饱和度 TPOR:总孔隙度 MSIGTA:核磁总孔隙度 MPHITA:核磁有效孔隙度 MBVITA:核磁束缚水体积 TDAMSIG:TDA分析得到的总孔隙度 TDAMPHI:TDA分析得到的有效孔隙度 DMRP:DMR分析得到的总孔隙度 DMREPHI:DMR分析得到的有效孔隙度 RT:真电阻率

21

TMP3:核磁探头的温度 P1PRF:PR06 BIN 1#(T2=0.5ms)的孔隙度 P2PRF:PR06 BIN 1#(T2=1ms)的孔隙度 P3PRF:PR06 BIN 1#(T2=2ms)的孔隙度 P4PRF:PR06 BIN 1#(T2=4ms)的孔隙度 P2AF:A组的 BIN 2#(T2=8ms)的孔隙度 P3AF:A组的 BIN 3#(T2=16ms)的孔隙度 P4AF:A组的 BIN 4#(T2=32ms)的孔隙度 P5AF:A组的 BIN 5#(T2=64ms)的孔隙度 P6AF:A组的 BIN 6#(T2=128ms)的孔隙度 P7AF:A组的 BIN 7#(T2=256ms)的孔隙度 P8AF:A组的 BIN 8#(T2=512ms)的孔隙度 P9AF:A组的 BIN 9#(T2=1024ms)的孔隙度 P10AF:A组的 BIN 10#(T2=2048ms)的孔隙度 P2AF:B组的 BIN 2#(T2=8ms)的孔隙度 P2AF:B组的 BIN 3#(T2=16ms)的孔隙度 P2AF:B组的 BIN 4#(T2=32ms)的孔隙度 P2AF:B组的 BIN 5#(T2=64ms)的孔隙度 P2AF:B组的 BIN 6#(T2=1288ms)的孔隙度 P2AF:B组的 BIN 7#(T2=256ms)的孔隙度 P2AF:B组的 BIN 8#(T2=512ms)的孔隙度 P2AF:B组的 BIN 9#(T2=1024ms)的孔隙度 P2AF:B组的 BIN 10#(T2=2048ms)的孔隙度

输出曲线
FTEMP:地层温度 RWAT:地层温度下的地层水电阻率 RMFAT:地层温度下的泥浆滤液电阻率 SALW:地层水矿化度(单位:kppm) SALMF:泥浆滤液矿化度(单位:kppm) RHOM:地层水密度 RHOMF:泥浆滤液密度 SWBCOR:校正过的束缚水饱和度 TPORCOR:校正过的总孔隙度 PHIT:来自MRIAN的总孔隙度 EPOR:有效孔隙度 MSIG:来自MRIAN的核磁总孔隙度 MPHI:来自MRIAN的核磁有效孔隙度 MBVI:来自MRIAN的核磁束缚水体积 WW:含水饱和度100%时的视W变量(WQ在SWT=1的边界) WI:束缚水条件下的视W变量(WQ在SWT=Swirr的边界) WQ:双水模型中估算的W系数 RO:地层电阻率(100%含水) RI:地层电阻率(束缚水状态) CBVWT:总孔隙度中的含水孔隙体积(TPOR×SWT)

22

CBVWE:有效孔隙度中的含水孔隙体积(CBVWT-CLPOR) MBVIT:总束缚水(MBVITA+泥质束缚水) MSWT:总含水饱和度 MSWE:有效含水饱和度(CBVWE/EPOR) MSWIT:总束缚水饱和度(MBVIT/TPOR) MFFID:总的自由流体指数 CCW:等效泥质束缚水导电率(单位:ms) RTR:双水模型重建的RT SPR:双水模型重建的SP(单位:mv) -3 2 SWBNEW:渗透率(用MPHI) (单位:×10 μm ) -3 2 MPERMW:水的渗透率(用MPHI) (单位:×10 μm ) -3 2 MPERMH:烃的渗透率(用MPHI) (单位:×10 μm ) -3 2 PMRI:渗透率(用TPOR) (单位:×10 μm ) -3 2 PMRIW:水的渗透率(用TPOR) (单位:×10 μm ) -3 2 PMRIH:烃的渗透率(用TPOR) (单位:×10 μm )

八、时间域TDA分析
通常轻烃有较长的T1,水的T1较短,因此对孔隙水而言,较短的极化时间足以磁化,而轻烃则 需较长的极化时间才能完全磁化。理论上讲,两个T2谱相减,水信号抵消,油与气的信号则余 留在差谱中,由此识别油气,这就是差谱法。但实际上由于噪声影响,这种差谱的定性方法不 可靠。在应用中需要做时间域分析TDA(TIME DOMAIN ANALYSIS) ,来利用双TW信息完成对地层 流体的识别和油气定量评价。

(一) 、T2T1-EVENT
作用:
用于从差异较大的区间中,对可能存在的流体进行T1和T2的搜索。

运行方法:
Model Launcher→Numar→T2T1EVENT 输入上一步处理过的CLS文件; 点击EDIT→CURVE MAP,确定输入、输出曲线名(重定向); 点击OPTION→create parameter database,生成新的参数文件,或File→Open parameter database,打开已存在的参数文件; 点击OPTION→RUN PARAMETER EDITOR,对参数文件进行编辑; 点击SAVE SESSION FILE,生成各种文件集合,以备再处理之用,最后点击RUN,运行此程序

参数设置
1、SET UP参数 T1SEARCH——T1搜索的标志(=0不搜索,=1搜索) T2SEARCH——T2搜索的标志(=0不搜索,=1搜索) NUMECHO——差谱中回波个数 TE——回波间距(单位:ms) TWB——短等待时间(单位:ms) TWA——长等待时间(单位:ms)

23

2、T1搜索参数 T1STKW——T1累加的CPMG个数(隐含40,范围1-200) T1FBLN——最大孔隙度对比的第一个区间(隐含6,范围4-12) T1LBLN——最大孔隙度对比的最后一个区间(隐含8,范围4-12) MINT1——预计的最小T1(隐含1500,范围10-6000,单位:ms) MAXT1——预计的最大T1(隐含6600,范围10-18000,单位:ms) 3、T2搜索参数 T2STKW——T2累加的CPMG个数(隐含40,范围1-200) FECHO——T2搜索的第一个回波(隐含2,范围1-100) LECHO——T2搜索的最后一个回波(隐含400,范围100-1000) T2XMIN——第一相的T2最小值(隐含10,范围3-2000,单位:ms) T2XMAX——第一相的T2最大值(隐含100,范围30-2000,单位:ms) T2YMIN——第二相的T2最小值(隐含200,范围30-2000,单位:ms) T2YMAX——第二相的T2最大值(隐含900,范围30-2000,单位:ms) MNPHIT2——最小的做搜索的孔隙度(单位:%) T2DELTA——T2搜索的步长(单位:ms) T2EPS——T2搜索的误差容限(单位:ms) T2MXITR——T2预测的最大迭代次数(隐含99,范围3-99999) T2COMPS——进行搜索的T2的个数(隐含2,范围1-2)

输入曲线
DEPTHA:A组与B组的平均参考深度 P1A:第一区间的孔隙度幅度(长TW) P2A:第二区间的孔隙度幅度(长TW) P3A:第三区间的孔隙度幅度(长TW) P4A:第四区间的孔隙度幅度(长TW) P5A:第五区间的孔隙度幅度(长TW) P6A:第六区间的孔隙度幅度(长TW) P7A:第七区间的孔隙度幅度(长TW) P8A:第八区间的孔隙度幅度(长TW) P9A:第九区间的孔隙度幅度(长TW) P10A:第十区间的孔隙度幅度(长TW) P11A:第十一区间的孔隙度幅度(长TW) P12A:第十二区间的孔隙度幅度(长TW) P1B:第一区间的孔隙度幅度(短TW) P2B:第二区间的孔隙度幅度(短TW) P3B:第三区间的孔隙度幅度(短TW) P4B:第四区间的孔隙度幅度(短TW) P5B:第五区间的孔隙度幅度(短TW) P6B:第六区间的孔隙度幅度(短TW) P7B:第七区间的孔隙度幅度(短TW) P8B:第八区间的孔隙度幅度(短TW) P9B:第九区间的孔隙度幅度(短TW) P10B:第十区间的孔隙度幅度(短TW) P11B:第十一区间的孔隙度幅度(短TW)

24

P12B:第十二区间的孔隙度幅度(短TW) EDBST:用于增强回波串的差 EDIF:长短TW的回波串差

输出曲线
T1APP:T1的估算值(单位:ms) T1ERR:T1搜索的误差最小值 DEPTHT1:T1搜索输出的深度 T2X1:单相搜索的T2 PORX1:单相搜索的孔隙度 T2ERRX:单相搜索的误差最小值 T2X:双相搜索的T2(短分量) PORX:短分量的孔隙度 T2Y:双相搜索的T2(长分量) PORY:长分量的孔隙度 T2ERRXY:双相搜索的误差最小值

说明
1、T2T1EVENT的处理文件为时间域文件*.m.cls。 2、在处理之前,用02-echostrip-ab-event.spc和09pte-t2t1event-qcpor.spc显示有关曲线,找出 哪几个区间差异较大,用于T1搜索。 3、对于孔隙度很小的井段,T1APP、T2X、T2Y意义不大。对于孔隙度较大的井段,若T2X1(单相) 与T2X、T2Y中的某一条基本重叠,说明那一部分占主要地位;若T2X1(单相)界于T2X、T2Y中间, 说明存在两相流体。 4、若T2ERR等误差曲线数值高,可适当增大MNPHIT2,做重新处理。

(二) 、TDA—COMP
作用: 运行方法:
Model Launcher→Numar→TDA-COMP 输入上一步处理过的CLS文件; 点击EDIT→CURVE MAP,确定输入、输出曲线名(重定向); 点击OPTION→create parameter database,生成新的参数文件,或File→Open parameter database,打开已存在的参数文件; 点击OPTION→RUN PARAMETER EDITOR,对参数文件进行编辑; 点击SAVE SESSION FILE,生成各种文件集合,以备再处理之用,最后点击RUN,运行此程序 3.参数设置 SET UP参数 FECHOTDA——用于TDA分析的第一个回波序号 LECHOTDA——用于TDA分析的最后一个回波序号 TE——回波间距 TWB——短等待时间 TWA——长等待时间 FLUIDFLG——流体存在标志 1=气与油
25

2=气与水 3=油与水 4=气 5=油 6=水 流体的核磁参数 T1GAS——气的T1 T2GAS——气的T2 HIGAS——气的含氢指数 T1OLL——油的T1 T2OLL——油的T2 HIOLL——油的含氢指数 T1WTR——水的T1 T2WTR——水的T2 HIWTR——水的含氢指数

输入曲线
EDIF:长TW组与短TW组的差谱

输出曲线
PHIGU:气的孔隙度估算值(未经T1和HI校正) PHIOU:油的孔隙度估算值(未经T1和HI校正) PHIWU:水的孔隙度估算值(未经T1和HI校正) PHIGA:气的视孔隙度(长TW) PHIOA:油的视孔隙度(长TW) PHIWA:水的视孔隙度(长TW) PHIG:气的孔隙度估算值(经过T1和HI校正) PHIO:油的孔隙度估算值(经过T1和HI校正) PHIW:水的孔隙度估算值(经过T1和HI校正) YFITTDA:TDA回波串EDIF的理论值 SNRT:累加后EDIF的信噪比

3、做TDA的流程
ECHO-STRI P
处 理 之 后 用 02-echostrip-ab-event.spc 09pre-t2t1event-qcpor.spc 显示有关曲线,找出哪几个区间差异较大,用于T1搜索。 和

T2-TOOLKI T

令HCFLG=1 计算理论烃性质,从其输出中主要使用T1GAS、HIGAS。 因T2GAS与TE有关, 故此处意义不大。 因OILVISC不易得到准确数值, 故此处T1OIL 也意义不大。

九、扩散分析DIFAN
通常,水的扩散系数比较大,而高粘度原油的扩散系数比水小。观测的横向弛豫时间T2是 流体的扩散系数D、回波间距TE、以及磁场梯度G的函数。对于固定的G,改变TE,高粘度油与自
26

由水的T2将发生不同程度的变化,即自由水的T2比高粘度油以更快的速度减小。通过合理地选 择TE,甚至可以在T2分布上把高粘度原油与自由水完全分开。比较长、短TE的T2分布,找出油 、水的特征信号,从而识别流体的方法叫做移谱法。移谱法只能是定性的,而且也不可靠,因 为油、水的扩散系数等参数以及谱位移的大小都不能直接获得。在应用中,基于上述原理,通 过扩散分析或扩散增强分析来实现对高粘度油的识别和定量评价。

(一)T2-DIFSEL
作用:
通过对T2分布谱进行分析,计算出长、短TE的几何平均值、峰值、半峰值,为运行DIFAN提供数 据。

运行方法:
Model Launcher→Numar→T2-DIFSEL 输入上一步处理过的CLS文件; 点击EDIT→CURVE MAP,确定输入、输出曲线名(重定向); 点击OPTION→create parameter database,生成新的参数文件,或File→Open parameter database,打开已存在的参数文件; 点击OPTION→RUN PARAMETER EDITOR,对参数文件进行编辑; 点击SAVE SESSION FILE,生成各种文件集合,以备再处理之用,最后点击RUN,运行此程序

参数设置
1、确定峰值的限制参数 AMINA——用于短TE组确定峰值的最小孔隙度(单位:%,隐含1.2,范围0.5—5) AMINB——用于长TE组确定峰值的最小孔隙度(单位:%,隐含1.2,范围0.5—5) 2、确定几何平均值的参数 T2GMLA——短TE组计算几何平均T2GMS的下限(单位:ms,隐含32,范围4-2048) T2GMUA——短TE组计算几何平均T2GMS的上限(单位:ms,隐含1024,范围4-2048) T2GMLB——长TE组计算几何平均T2GML的下限(单位:ms,隐含32,范围4-2048) T2GMUB——长TE组计算几何平均T2GMS的上限(单位:ms,隐含32,范围4-2048) PMINA——短TE组计算几何平均T2GMS的最小孔隙度界限 PMINB——长TE组计算几何平均T2GML的最小孔隙度界限 T2SELFLG——确定DIFAN中T2用哪组数的选择标志,同时控制T2S和T2L 1=用T2GM 2=用T2DIF1 3=用T2DIF

输入曲线
MPHITA:核磁的有效孔隙度 TASPEC:短TE的T2分布 TBSPEC:长TE的T2分布

输出曲线
T2GMS:短TE组的T2几何平均值(单位:ms) T2DIF1S:短TE组的T2m峰值(单位:ms) T2DIFS:短TE组的T2半峰值(单位:ms) T2S:用于扩散分析的短TE的T2值(单位:ms) T2GML:长TE组的T2几何平均值(单位:ms)
27

T2DIF1L:长TE组的T2峰值(单位:ms) T2DIFL:长TE组的T2半峰值(单位:ms) T2L:用于扩散分析的长TE的T2值(单位:ms)

说明
T2-DIFSEL运行之后,可用07-t2difsel·spc在屏幕出图,观察参数的处理效果。T2GM、T2DIF1 、T2DIF应主要在有可动流体的部分计算。

(二)、DIFAN
作用:
使用双TE(A组为短TES,B组为长TEL),做差谱分析。

运行方法:
Model Launcher→Numar→DIFAN→ 输入上一步处理过的CLS文件; 点击EDIT→CURVE MAP,确定输入、输出曲线名(重定向); 点击OPTION→create parameter database,生成新的参数文件,或File→Open parameter database,打开已存在的参数文件; 点击OPTION→RUN PARAMETER EDITOR,对参数文件进行编辑; 点击SAVE SESSION FILE,生成各种文件集合,以备再处理之用,最后点击RUN,运行此程序

参数设置
1、SET UP参数 TES:短回波间距 TEL:长回波间距 CFRQ:仪器中心频率(B组长TE) SURFT——地表温度(单位:DEGF) BHT——井底温度(单位:DEGF) TD——井深(单位:ft) PRESSG--地层压力梯度(单位:PSI/ft) PRESSF--视地层压力(单位:PSI) PRESSD--视地层压力的深度(单位:ft) TOOLTYP——核磁仪器类型 0=MRIL-C 6"探头 1=MRIL-C 4.5"探头 2=MRIL Prime 6" 探头 3=MRIL Prime 4.875" 探头 TFLAG——温度单位标志(隐含0,范围0-1) 0=华氏度DEGF 1=摄氏度DEGC T2LIM——为固有T2的上限(隐含0,表示无限制) 2、流体参数 T2MIN——自由流体的最小固有T2

28

T2HY——烃的固有T2 RDDWHY——烃/水的扩散系数比值 ROSF——残余油饱和度 3、调节参数 GFACT——磁场梯度调节因子 T2SFACT——短T2调节因子 4、相对渗透率系数 RELPFL——相对渗透率模型(=0为Minas,=1为General) WVIS——地层水粘度(单位:cp) OVIS——地层中油的粘度(单位:cp)

输入曲线
FTEMP:地温数据 MPHITA:短TE的核磁有效孔隙度 MBVITA:最大束缚体积 MPERM:Coates模型渗透率 T2S:用于扩散分析的短TE的T2值 T2L:用于扩散分析的长TE的T2值 TMPS:核磁探头温度

输出曲线
T2INT:固有T2驰豫时间 T2DW:T2扩散的水界限 D:地层流体的扩散系数 DW:地层水的扩散系数 RDDW:D/DW的比值 IT2:固有T2的倒数,等于1/T2INT DIFSW:出自DIFAN的含水饱和度 DIFBVW:出自DIFAN的束缚水体积 DIFKRO:出自DIFAN的油相对渗透率 DIFKRW:出自DIFAN的水相对渗透率 DIFWCUT:出自DIFAN的预测产水率

29

核磁共振DPP处理参数文件命名规范
现在以X井为例,说明如下:

一 原始文件
原始深度域文件:X.1.cls:X→井名,1→第几趟测井,cls→文件格式。 原始时域文件:X.1.m.cls: X→井名,1→第几趟测井,m→时域标识,cls→文件格式。 要求提供原始数据基本信息清单,包括:深度范围和层位。



资料处理
(一) Split_MCLS 若为双 TW,双 TE 测井模式,在处理过程中要进行文件分割: X.1.abc.m.cls→1→第几趟测井,abc→分割出的回波,m→时域标识,cls→文件格式。 X.1.adc.m.cls→1→第几趟测井,abc→分割出的回波,m→时域标识,cls→文件格式。 (二)RESEQ 若为双 TW,双 TE 测井模式,在处理过程要支 GRP 和 CACT 进行重排: abc 组,参数文件命名 R1ABC.jbv,R→RESEQ,1→第几趟测井数据,ABC→回波。 adc 组,参数文件命名 R1ADC.jbv,R→RESEQ,1→第几趟测井数据,ADC→回波。 原则上要求保存*.ses 文件,命名方式为: abc 组,参数文件命名 R1ABC.ses,R→RESEQ,1→第几趟测井数据,ABC→回波。 adc 组,参数文件命名 R1ADC.ses,R→RESEQ,1→第几趟测井数据,ADC→回波。 (三)ECHO_STRIP abc 组,参数文件命名 ES1ABC.jbv,ES→ECHO_STRIP,1→第几趟测井数据,ABC→回波。 adc 组,参数文件命名 ES1ADC.jbv,ES→ECHO_STRIP,1→第几趟测井数据,ADC→回波。 原则上要求保存*.ses 文件,命名方式为: abc 组,参数文件命名 ES1ABC.ses,ES→ECHO_STRIP,1→第几趟测井数据,ABC→回波。 adc 组,参数文件命名 ES1ADC.ses,ES→ECHO_STRIP,1→第几趟测井数据,ADC→回波。 (四)T2_TOOLKIT abc 组,参数文件命名 T2T1ABC.jbv,T2T→T2_TOOLKIT,1→第几趟测井数据,ABC→回波。 adc 组,参数文件命名 T2T1ADC.jbv,T2T→T2_TOOLKIT,1→第几趟测井数据,ADC→回波。 原则上要求保存*.ses 文件,命名方式为:

30

abc 组,参数文件命名 T2T1ABC.ses,T2T→T2_TOOLKIT,1→第几趟测井数据,ABC→回波。 adc 组,参数文件命名 T2T1ADC.ses,T2T→T2_TOOLKIT,1→第几趟测井数据,ADC→回波。 (五)其它处理模块在*.jbv 和*.ses 参数文件命名中的简写 表 1 DPP 核磁处理模块参数文件命名规范

处理模块
PRE_MRIN SWB_MRIN MRIAN TDA_COMP TDA/MRIAN T2_DIFSEL

简写 PM SM MR TC TM T2D

ABC 组回波 PM1ABC.jbv PM1ABC.ses SM1ABC.jbv SM1ABC.ses MR1ABC.jbv MR1ABC.ses TC1ABC.jbv TC1ABC.ses TM1ABC.jbv TM1ABC.ses T2D1ABC.jbv T2D1ABC.ses

ADC 组回波 PM1ADC.jbv PM1ADC.ses SM1ADC.jbv SM1ADC.ses MR1ADC.jbv MR1ADC.ses TC1ADC.jbv TC1ADC.ses TM1ADC.jbv TM1ADC.ses T2D1ADC.jbv T2D1ADC.ses

注意;都是假设第 1 趟的测井数据。

声明: 本文档隶属《核磁共振(MRIL-P)测井资料处理方法与应用研究》项目

31


相关文章:
核磁共振测井资料处理及解释规范
4页 免费 P核磁测井数据处理与参数文... 31页 8财富值如要投诉违规内容,请到百度文库投诉中心;如要提出功能问题或意见建议,请点击此处进行反馈。 ...
核磁测井原理与解释
暂无评价|0人阅读|0次下载|举报文档 核磁测井原理与...4.3多维 NMR 流体特征描述 NMR 数据采集和处理方面...NMR 信号幅度是流体参数 T1、 T2、T1/T2 D 中...
测井数据标准化
暂无评价|0人阅读|0次下载|举报文档测井数据标准化...标准化工作是为测井解释储集层参数,进行油藏描述的...比如在测井曲 线名称的命名方式上均有所不同, ...
核磁共振测井简介
投入商业应用的核磁共振测井仪器的世界知 名测井服务...Schlumberger--CMR、Halliburton--MRIL-P、Baker hugh...[10]张宫. 核磁共振测井数据处理方法与软件开发研究...
核磁测井
核磁测井原理与解释 12页 7下载券 MRIL-P核磁测井...这两个 参数都根据测井目标预测的地层流体性质...3.3 多维 NMR 流体特征描述 NMR 数据采集和处理...
P型核磁处理软件
原始的 P核磁测井数据用 CLS 解编模块进行解编...保留 WIS 文件名中的 abc、adc、cde、bce 将有...本处理不需要参数文件,为与其他处理框架一致,单击...
页岩气单井测井资料处理与解释规范0625
与国际标准一致性程度的标识(征求意见稿) 2012.6....测井资料处理解释规范 SY/T 6617-2005 核磁共振测井...式中: Gs——吸附气含量,m /t; 3 vl p ......
测井采集、解释规程
测井资料 采集与处理解释操作规程 (送审稿) 中国...种类多,处理、解释的方法多,相应的标准、作业文件多...(七)MRIL-P核磁共振成像测井仪??? (八)水泥...
陆上地震勘探数据处理技术规程
陆上地震勘探数据处理技术规程 1 范围 本标准规定了...型号、线(束)号、文件号及地 震记录采集有关参数...c) 工区内的平均速度曲线或地震测井曲线。 d) 以...
常用测井曲线符号及单位(最规范版)
暂无评价|0人阅读|0次下载|举报文档常用测井曲线符号...核磁共振渗透率 标准回波数据 T2 分布数据 总孔隙度...PORG 气指数 CLASS CHR 阳离子交换能力与含氢量的...
更多相关标签:
核磁共振测井 | 核磁共振测井原理 | 核磁测井 | 核磁共振测井问题分析 | 核磁共振测井仪 | 核磁共振测井发展过程 | 核磁共振测井技术分析 | 核磁共振测井仪器 |