NET平台P型核磁操作手册操作手册.docx
《NET平台P型核磁操作手册操作手册.docx》由会员分享，可在线阅读，更多相关《NET平台P型核磁操作手册操作手册.docx（33页珍藏版）》请在冰豆网上搜索。
NET平台P型核磁操作手册操作手册
附件：
Forward.NET平台P型核磁共振成像分析软件
操作手册
渤海钻探工程有限公司测井分公司资料评价中心
2008年12月
一、P型核磁共振测前设计操作手册
（一）界面登录
1、点击桌面上的
图标，进入Forward.NET主界面；
2、找到精细解释，点击成象，选择P核磁共振测前设计模块；
3、如果是首次设计，在测前设计图标处单击鼠标右键，在下拉菜单出现后选择新建进入新设计；如果是修改设计，双击测前设计图标，在文件选择对话框出现后选择设计文件，进入修改。
（二）操作说明
通过输入油气粘度、温度、压力即可计算各相流体在不同TE条件下的T1、T2，从而选择合适的观测模式。
设计中可以有两种操作方式：
直接选择仪器观测模式，观察极化周期和T2分布；在极化曲线图可编辑的情况下，移动等待时间，在不同的等待时间设置下选择仪器观测模式，观察极化周期和T2分布。
界面的左方显示极化曲线图、T2谱和结论报告；右方是测井地质参数设置；界面上方有四个功能键。
1、极化曲线图
根据选定的模式、型号和参数自动绘出极化曲线，极化曲线图有两种操作方式
1）、利用功能键“选择观测模式”选定模式和型号，确认后，极化曲线图显示该仪器的极化曲线和等待时间（粉红色虚线）。
2）、将“能够编辑图形”键关闭，将光标移至等待时间线（粉红色虚线），当光标变为左右箭头时可以左右拖动光标修改等待时间，此时在合理的等待时间处会弹出“自动搜索观测模式”对话框，列出该等待时间附近可能对应的模式型号供选择。
如果等待时间（图中粉红色虚线）不够，可单击功能键“增加TW”，此时光标变成+TW模式，拖动光标到适当位置按下鼠标左键，新的等待时间线将出现在图上。
图面上最多可出现三条等待时间线。
2、T2谱分析图
极化曲线确定后，各组的理论T2谱分析被绘出。
油、气、水和泥浆滤液的选择在参数设置的“油藏流体类型”栏下。
3、测前设计结论
设计中选择的观测模式及参数设置在本结论中集总为报告。
4、功能键说明
选择观测模式：
点击选择观测模式，弹出观测模式对话框，用户可选择测井仪器，观测模式和具体型号进行设计。
增加TW：
单击功能键“增加TW”，此时光标变成+TW模式，拖动光标到适当位置按下鼠标左键，新的等待时间线将出现在图上。
图面上最多可出现三条等待时间线。
显隐设置参数对话栏：
这是一个开关，点击它，参数设置栏在显示和隐藏之间变换以保证对图形观察的空间。
新增或修改参数：
单击新增或修改参数允许用户添加新的泥浆参数或仪器参数。
当前只允许添加油基泥浆，新泥浆类型名一定以OBM三个字符开头。
单击增加，在固定行标题下出现新空行，给出油基泥浆类型名称和参数，键入保存，确定，新的泥浆类型将被纳入参数设置的泥浆类型表，供设计中选择。
参数对极化周期和扩散系数的影响由选定的公式1和公式2确定，这两个公式仅用于油基泥浆。
在增加仪器参数前，要先增加仪器类型，打开核磁模版目录中的tooltype.txt文件,在文件末尾字符串”New”的下边输入新的仪器类型及该仪器所能测的观测模式，然后单击仪器参数按钮，可增加新仪器的参数，增加后保存。
二、P型核磁共振解释软件操作手册
Forward.NET平台P型核磁共振解释软件可以分别处理时间域（*.m.wis）和深度域（*.d.wis）的核磁共振测井数据。
整个运行过程都是由“P型核磁共振成像分析”模块来完成。
启动模块，选定WIS文件后，模块首先对WIS文件作检查，根据检查结果作对处理内容的自动定位
如果WIS文件是双TW双TE的时间域原始数据首先要进行文件拆分，将其拆分为单TW双TE或单TE双TW的新WIS文件，运行的处理步骤为“文件分离”；
如果WIS文件是单TW双TE或单TE双TW的时域数据，运行的处理步骤为“回波拆分”、“时域分析计算”和“时深转换”等；
如果WIS文件是深度域数据则可运行的处理步骤为“岩石物理计算”、“与常规测井结合确定总孔隙度”、“与常规测井结合确定泥质束缚水饱和度”、“与常规测井结合计算含油饱和度”等。
（一）、文件分离
1、功能
将DTWE数据文件分解成DTW、DTE数据文件，每组数据文件包括三组回波串，以便于后续测井资料评价。
当MRIL-Prime型核磁共振测井仪以DTWE模式工作时，采集的回波串包括A、B、PRO6（或C）、D、E五组，其中A、B两组是由短TE（1.2ms或0.9ms）采集的双TW模式；D、E则是由长TE采集的双TW模式，所以A与D组合可以得到一个双TE观测，B与E是另一个双TE观测（只是由于它们都是短TW采集到的，通常不使用）。
PR06组（或C组）用于粘土束缚水的观测，其采集参数是固定的，即TE＝0.6ms，TW=0.2s，NE=10，RA＝50。
因此，五组回波串的分组方式有abc、adc、cde、bce四种。
2、输入文件
原始的*.m.wis（时间域DTWE数据文件）
3、输出文件
自动形成*.abc.m.wis文件或*.acd.m.wis文件或*.cde.m.wis文件或*.bce.m.wis文件，分别保存在选择的目录下的同名文件夹内。
其中abc代表短回波间隔的双TW观测方式；
cde代表长回波间隔的双TW观测方式；
acd代表长等待时间的双TE观测方式；
bec代表短等待时间的双TE观测方式。
拆分以后的文件均以ABC三组方式出现
在双TW观测模式下A为长等待时间组，B为短等待时间组。
在双TE观测模式下A为短回波间隔组，B为长回波间隔组。
保留WIS文件名中的abc、adc、cde、bce将有利于区分观测模式与文件的对应关系。
4、运行方式
1）单击主界面的精细解释，选择成象；
2）双击P型核磁共振成像分析模块；
3）选中要分解的时间域DTWE数据文件（*.m.wis），模块检查后将自动定位到“文件分离”；
4）点击开始处理
，在弹出的“文件提取”界面选择需要得到的观测方式和存储目录；
5）选择“确定”运行程序，自动生成与选择观测方式对应的文件夹。
（二）、回波拆分
1、功能
核磁测井原始资料是时域数据且原始回波按等时记录于RAMP,回波拆分模块将记录回波按组拆分，同时完成对回波串累加和作相位校正；对回波串多指数拟合，得到离散的T2分布及各区间孔隙度；为TDA_COMP的计算提供长、短TW回波串的差谱EDIF；输出孔隙度MPHIX用作深度校正。
注意,该程序不能对bce、cde组数据进行处理。
2、输入曲线
GRPMRIL-P各采集组的编码
CACTMRIL-P各采集回波串编码（0-39）
CECH采集设置中的回波串个数
RAMPECHO幅度（%）
RPHAECHO相位角（度）
TDEP回波串采集对应的电缆深度码
PWCO功率校正因子
SACO矿化度校正因子
TMCO温度校正因子
GAIN输入增益
3、输出曲线
A组校正后的回波串的幅度、相位、实部、虚部
B组校正后的回波串的幅度、相位、实部、虚部
各组各部分的离散T2值
4、运行说明
运行方式与文件分离模块相似。
本处理不需要参数文件，但为与其他处理框架一致，单击处理按钮弹出模式选择对话框。
模块根据读出数据的回波个数自动搜索可能的模式型号，等待时间和回波间隔显示为空。
用户选定本次测井具体的模式型号后等待时间和回波间隔会自动填入。
模式型号选择后，确定相位校正参数和解谱参数，如果用户需要改变这个值，可钩选该方框活化下面的两个输入框，然后键入新值，然后单击确定进入运行。
相位角校正
：
将通道1和通道2的数据旋转一个角度φ，旋转之后，一道将主要包含NMR信号（其幅度是回波串的实数部分），另外一道将主要包含噪音（其幅度是回波串的虚数部分）,回波串的实数部分被转换为T2分布。
用下式计算角度φ
这里i是回波串中的第i个回波，k是用做相位角计算的回波数目，模块隐含用第2到底9个回波计算。
如果用户需要改变这个值，可钩选该方框活化下面的两个输入框，然后键入新值。
平滑系数：
滤波平滑系数。
解谱首波：
解谱中，时间轴上小于解谱首波的回波将被忽略。
解谱末波：
解谱中，时间轴上大于解谱末波的回波将被忽略。
（三）时深转换
1、功能
把时间域的*.m.wis文件转换成以深度域的*.d.wis文件。
2、运行说明
点击
按钮弹出选择转换位置对话框
“将转换的深度域曲线写入新的WIS文件”，输入相应的目录及文件名，选择确定后开始转换；
“将转换的深度域曲线写入存在的WIS文件”，保存文件对话框弹出后用户可选择硬盘上存在的WIS文件。
此时模块不创建新文件，但所选择的WIS文件必须是深度域的文件，否则转换失败，你收到信息
“核磁仪器偏置深度”：
该深度由测量过程中测量的延迟引起。
“GR仪器偏置深度”：
该深度由测量过程中测量的延迟引起。
在根据原始记录的时间、深度记录转换之后，将根据两者之差作深度移动。
（四）、岩石物理计算
1、功能
计算核磁各种孔隙度；计算核磁渗透率；T2谱的合并，T2的0.5，1，2，4msbin数据由TE=0.6ms的回波串拟合结果提供，8，…，2048ms数据由TE≠0.6ms的回波串拟合结果提供；由离散的T2分布得到连续的T2谱，（200个元素的向量）；计算油气的NMR特性（根据温度压力及理论公式）；对所有的孔隙度曲线进行滤波处理；对TDA-COMP的计算结果进行滤波。
2、输入曲线
ECHOSTRIP和PROCESS-T2D处理后的输出曲线。
3、输出曲线
各组的T2谱和各种孔隙度、渗透率流体核磁特性等曲线。
4、输入参数
FILTFLG滤波选择0=NO,1=YES
FLTTYP滤波方式0=汉宁窗,1=三角形,2=矩形窗
FLTLEN孔隙度滤波长度（采样点数）
T2CUTOFFT2BVI截止值
SBVICOEF台阶状模型系数（SBVI）
T2SPLICET2谱拟合时前几个BIN从PR谱拼合
CCOEF渗透率系数（COATES模型）
CTHRES渗透率计算时最小孔隙截止因子
TE回波间隔一般取0.9
5、运行说明
在做完回波拆分后应用本模块对回波串做时间域分析。
在作岩石物理计算处理以前，处理合并带常规GR、孔隙度、电阻率资料，并将核磁的GAMM曲线与常规GR曲线对比，进行核磁测井资料校深；亦可直接作岩石物理计算处理，这样软件可以自动拟合得到一条GR曲线，然后再合并曲线校深，但是此时T2谱将发生变形，需要再运行一次该模块。
（五）、时域分析
1、时域分析（搜索T1T2）
（1）功能
确定油的T1、T2，气的T2或水的T1、T2值。
（2）输入曲线
DEPTHA回波串A的参考深度
P01A、P02A、P1A、…、P10A回波串A拟合后的各时间带的孔隙度
P01B、P02B、P1B、…、P10B回波串A拟合后的各时间带的孔隙度
EDBST回波串差的增强因子
（3）输出曲线
T1APPT1的估算值
T1ERRT1搜索误差最小值
DEPTHT1T1搜索的输出深度
T2X1单相搜索的T2
PORX1单相搜索的孔隙度
T2ERRX单相搜索误差最小值
T2X双相搜索的短分量T2
PORX双相搜索的短分量孔隙度
T2Y双相搜索的长分量T2
PORY双相搜索的长分量孔隙度
T2ERRXY双相搜索误差最小值
（4）输入参数
TISEARCH是否进行T1搜索，0不搜索，1搜索
T2SEARCH是否进行T2搜索，0不搜索，1搜索
NUMECHO回波差的回波个数
Te回波间隔
TWB记录一个完整的短Tw回波串记录的时间
TWA记录一个完整的长Tw回波串记录的时间
TISTKWT1累加次数（40的倍数，次数越多信噪比越高）
TIFBINT1搜索的第一个bin
TIFBINT1搜索的最后一个bin
MINTIT1的最小时间范围（范围越大搜索精度越高）
MAXTIT1的最大时间范围
T2STKWT2累加次数（40的倍数）
FECHO搜索T2的第一个回波
LECHO搜索T2的最后一个回波
T2XMIN1相的T2最小值
T2XMAX1相的T2最大值
T2YMIN2相的T2最小值，T2XMAXT2YMAX2相的T2最大值
MNPHIT2计算长的最小孔隙度门槛值
T2DELIAT2的搜索步长，一般取1.0
T2EPST2的搜索的误差容限
T2MXITRT2的搜索的最多迭代次数
T2COMPS搜索几组T2值（单相或两相）
2、时域分析计算
（1）功能
时域分析（TDA）依赖于这样一个实事，即不同的流体有不同的极化率，或不同的T1驰豫时间。
气和轻质油（粘度小于5cp）通常比水的T1要长得多。
时域分析提供
冲洗带流体类型；经含烃指数校正后的核磁孔隙度；各组分流体饱和度的分析（冲洗带）；进入TDA计算的数据已经过对双TW双TE的分离，经过了回波拆分。
（2）输入曲线
EDIFA组、B组回波串之差
（3）输出曲线
PHIGU未作任何校正的气相孔隙体积
PHIOU未作任何校正的油相孔隙体积
PHIWU未作任何校正的水相孔隙体积
PHIGA未作HI校正的气相孔隙体积
PHIOA未作HI校正的油相孔隙体积
PHIWA未作HI校正的水相孔隙体积
PHIG气相孔隙体积
PHIO油相孔隙体积
PHIW水相孔隙体积
（4）输入参数
T1GAS气的T1值
T2GAS气的T2值
HIGAS气的含烃指数
T1OIL油的T1值
T2OIL油的T2值
HIOIL油的含烃指数
T1WTR水的T1值
T2WTR水的T2值
HIWTR水的含烃指数
FECHOTDA做TDA计算首回波位置
LECHOTDA做TDA计算末回波位置
FLUIDFLG流体类型。
1油气混合相，2气相，3气水混合相，4油相，
5油水混合相，6水相（选择5）
（5）运行方式
1）时间域数据（*.m.wis）在完成“回波拆分”后，选择“时域分析（搜索T1T2）”模块；
2）点击参数编辑
，弹出“可视参数编辑”，在编辑器内点击
，可以显示参数名称，根据需要编辑参数并保存；
3）点击开始处理
，完成T1、T2搜索；
4）选择“时域分析计算”模块，编辑参数并处理，完成时域分析。
（6）流体T1、T2计算器说明
如果用户没有参数的经验值，可采用理论参数作初始处理。
单击功能按钮
“计算流体T1、T2值”，弹出“计算核磁特性参数”对话框，这是一个输入参数的辅助工具，请注意对话框右侧的计算结果除气体密度外均为输入参数值，如果用户有经验参数可直接输入。
如果需要本计算器来计算，可选定油藏类型、测井仪器、地层温度、地层压力、油的粘度、扩散系数、气体密度。
输入参数栏说明
油藏类型砂岩层（按渗透率分类）、碳酸盐岩。
测井仪器MRIL6”MRIL4，7-8”
选择计算气体密度方法气体密度由两种方法输入，“用户输入”将使用用户在输入参数栏气体密度输入值计算T1、T2。
“由温度压力计算”将根据用用户在输入参数栏输入的温度压力计算气体密度。
按下后，计算出右侧的所有参数。
在精细解释深度棒上选定一个层，此按钮被活化，按下后将计算结果写入该层段参数。
按下后将计算结果写入所有层段参数。
按下后，所有输入状态被保留，退出计算核磁特性参数对话框，下一次打开计算核磁特性参数时保留参数被显示。
如果用对话框的关闭X，则输入状态不保留。
（六）扩散分析
在DTE测井模式下，可根据TEL下最大水线T2DW右侧谱峰的有无、大小识别评价油气层。
扩散分析计算分为确定可动流体T2、扩散分析预处理和计算含水饱和度三部完成的。
1、确定可动流体T2
（1）功能
对T2谱进行分析，计算出长、段TE的几何平均值、峰值、半峰值、为扩散分析计算提供数据。
（2）输入曲线
TASPEC短TE的T2谱
TBSPEC长TE的T2谱
（3）输出曲线
T2GMS短TE组的T2几何平均值
T2DIF1S短TE组孔隙度峰值对应的T2
T2DIFS短TE组孔隙度半幅点的T2
T2S用于扩散分析短TE组的T2
T2GML长TE组的T2几何平均值
T2DIF1L长TE组孔隙度峰值对应的T2
T2DIFL长TE组孔隙度半幅点的T2
T2L用于扩散分析长TE组的T2
（4）输入参数
AMINA用于短TE组确定峰值的最小孔隙度
AMINB用于长TE组确定峰值的最小孔隙度
T2GMLA短TE组计算几何平均T2GMS的下限
T2GMUA短TE组计算几何平均T2GMS的上限
T2GMLB长TE组计算几何平均T2GML的下限
T2GMUB长TE组计算几何平均T2GML的上限
PMINA短TE组计算几何平均T2GMS的最小孔隙度界限
PMINB长TE组计算几何平均T2GML的最小孔隙度界限
T2SELFLG扩散分析计算时使用哪组T2,1-T2GM2-T2DIF13-T2DIF
2、扩散分析预处理
（1）功能
计算流体扩散系数，为扩散分析计算提供数据
（2）输入曲线
T2S用于扩散分析的短TE的T2值
T2L用于扩散分析的长TE的T2值
TBSPEC长TE的T2谱
（3）输出曲线
T2INT固有T2值
T2DW水的T2界限（长TE）
D地层流体的扩散系数
DW水的扩散系数
RDDWD与DW比值
IT2固有T2值的倒数
（4）输入参数
TES短回波间隔
TEL长回波间隔
FRQ仪器中心频率（现场提供）
SURFT地表温度
BHT井底温度
TFLAF温度单位,1-华氏度2-摄氏度
TD井深（米）
PRESSG地层压力梯度（PSI/米）
PRESSF视地层压力（PSI）
PRESSD视地层压力的深度（米）
TOOLTYPE核磁仪器类型,1-MRIL-C6"MRIL-C4.5"2-MRIL-PRIME6"3-MRIL-PRIME4.875"（现场提供）
T2LIM固有T2的上限（0表示无限制）
GFACT磁场剃度调节因子（现场提供）
T2SFACT短T2调节因子
3、扩散分析计算含水饱和度
（1）功能
对双TE（A组为短TES，B组为长TEL）做移谱分析。
（2）输入曲线
MPHITA有效孔隙度（短TE）
MSIGTA有效孔隙度（短TE）
MBVITA毛管束缚水孔隙度（短TE）
MPERMCoates模型渗透率
T2S用于扩散分析的短TE的T2值
T2L用于扩散分析的长TE的T2值
T2INT固有T2值
IT2固有T2值的倒数
RDDWD与DW比值
（3）输出曲线
DIFSW水饱和度
DIFBVW水的孔隙体积
DIFKRO油的相对渗透率
DIFKRW水的相对渗透率
DIFWCUT持水率估计值
（4）输入参数
T2MIN自由流体的最小固有T2
T2HY烃的固有T2
RDDWHY烃与水的扩散系数的比值
ROSF残余油饱和度（%）
T2SFACT短T2调节因子
RELPFL相对渗透率模型,1-米那斯方法2-普通方法3-琼斯方程4-彼尔逊方程
WVIS地层水粘度（cp）
OVIS地层油的粘度（cp）
（5）运行说明
为确定各种流体参数，可使用扩散分析图版进行分析。
单击功能按钮
“显示确定扩散分析图版”。
图版显示到文档中，双击图版，将X轴设置为RDDW（流体扩散度与地层水扩散度的比值），将Y轴设置为IT2（流体固有T2值的倒数），将“能够编辑图形”键关闭，图中兰色虚线为视水线标尺，上边是100%视水线标尺，下边是0%视水线标尺，移动光标到视水线标尺处，当光标变为上下箭头时，上下移动光标，可拖动视水线标尺，移动时程序自动控制两条标尺保持平行，参数框中的参数也随之改变，单击流体参数框中的写入选中层字段可将确定的参数写入到选中的层段中，单击写入所有层字段可将确定的参数写入到所有层段中。
（七）标准T2测井处理（MRIAN）
主要是利用核磁共振测井资料与常规的孔隙度、电阻率资料结合进行综合分析、提供准确的储层物性参数、原状地层含油饱和度计算。
包括三个解释模块，分别为与常规测井结合求总孔隙度模块、与常规测井结合求泥质束缚水饱和度模块和与常规测井结合求含油饱和度模块。
1、与常规测井结合求总孔隙度
（1）功能
准确求取总孔隙度фt，计算从核磁得到的Swb，提供后续处理所需的各种交会图数据。
（2）输入曲线
MSIGTA总孔隙度（由深度域曲线得到）
MPHITA有效孔隙度（由深度域曲线得到）
MBVITABVI值
TDAMSIG总孔隙度（由时间域曲线得到）
TDAMPHI有效孔隙度（由时间域曲线得到）
MCBW粘土束缚水孔隙度
NPHI补偿中子
RHOB体积密度
DT声波时差
RT深探测电阻率
（3）输出曲线
FTEMP地层温度
RWAT地层温度下的地层水电阻率
RMFAT地层温度下的泥浆滤液电阻率
SALW地层水矿化度（来自电阻率-矿化度图版）
SALMF泥浆滤液矿化度（来自电阻率-矿化度图版）
RHOW地层水密度
PHID计算的密度孔隙度
PHIN计算的中子孔隙度
PHIS计算的声波孔隙度
PHIX中子-密度交会孔隙度
TPOR地层总孔隙度
RWA视地层水电阻率
CWA视地层水导电率
SWBMRIMPHITA确定的束缚水饱和度
DPND中子与密度孔隙度的差
DPTD声波与密度孔隙度的差
RPND中子与密度孔隙度的比值
DPDM密度与MRIL孔隙度的差
（4）输入参数
TFLAF温度单位1=华氏度2=摄氏度
SURFT地表温度
BHT井底温度
TD井深（米）
PRESSG地层压力梯度（PSI/米）
PRESSF地表压力（PSI）
PRESSD地表深度（米）
NEMOVE中子测井刻度校正量
INNEUTFL中子骨架1=砂岩2=灰岩3=白云岩
TPORFL用什么方式确定总孔隙度
1=中子2=-密度交会3=密度4=中子5=声波6=声波-7=密度
8=声波-中子9=核磁10=校正后的核磁孔隙度
RWREF参考点处地层水电阻率
TWREF参考点温度
RMFREF参考点处泥浆电阻率
TMFREF参考点泥浆温度
ARCHIE阿尔奇A
EXPHI阿尔奇M
RHOFL流体密度
DMA骨架视密度
DTMA骨架时差
SONICFL计算声波孔隙度的公式,1=Raiga-Clemenceau2=Chapman-Modified
2、与常规测井结合求泥质束缚水饱和度
（1）功能
计算除核磁测井以外所有其它方法计算的Swb，并将Swb与核磁的Swb进行比较，选择一个最佳的Swb。
在进行Swb-MRIAN解释前要对前一步骤解释结果作一些交会图
升级会员 