NET平台P型核磁操作手册操作手册.docx

1、NET平台P型核磁操作手册操作手册附件： Forward.NET平台P型核磁共振成像分析软件操作手册渤海钻探工程有限公司测井分公司资料评价中心2008年12月一、P型核磁共振测前设计操作手册（一）界面登录1、点击桌面上的图标，进入Forward.NET主界面；2、找到精细解释，点击成象，选择 P核磁共振测前设计模块；3、如果是首次设计，在测前设计图标处单击鼠标右键，在下拉菜单出现后选择新建进入新设计；如果是修改设计，双击测前设计图标，在文件选择对话框出现后选择设计文件，进入修改。（二）操作说明通过输入油气粘度、温度、压力即可计算各相流体在不同TE条件下的T1、T2，从而选择合适的观测模式。设计

2、中可以有两种操作方式：直接选择仪器观测模式，观察极化周期和T2分布；在极化曲线图可编辑的情况下，移动等待时间，在不同的等待时间设置下选择仪器观测模式，观察极化周期和T2分布。界面的左方显示极化曲线图、T2谱和结论报告；右方是测井地质参数设置；界面上方有四个功能键。1、极化曲线图 根据选定的模式、型号和参数自动绘出极化曲线，极化曲线图有两种操作方式 1）、利用功能键“选择观测模式”选定模式和型号，确认后，极化曲线图显示该仪器的极化曲线和等待时间（粉红色虚线）。 2）、将“能够编辑图形”键关闭，将光标移至等待时间线（粉红色虚线），当光标变为左右箭头时可以左右拖动光标修改等待时间，此时在合理的等待时

3、间处会弹出“自动搜索观测模式”对话框，列出该等待时间附近可能对应的模式型号供选择。如果等待时间（图中粉红色虚线）不够，可单击功能键“增加TW”，此时光标变成+TW模式，拖动光标到适当位置按下鼠标左键，新的等待时间线将出现在图上。图面上最多可出现三条等待时间线。2、T2谱分析图 极化曲线确定后，各组的理论T2谱分析被绘出。油、气、水和泥浆滤液的选择在参数设置的“油藏流体类型”栏下。3、测前设计结论 设计中选择的观测模式及参数设置在本结论中集总为报告。4、功能键说明 选择观测模式：点击选择观测模式，弹出观测模式对话框，用户可选择测井仪器，观测模式和具体型号进行设计。 增加TW：单击功能键“增加TW

4、”，此时光标变成+TW模式，拖动光标到适当位置按下鼠标左键，新的等待时间线将出现在图上。图面上最多可出现三条等待时间线。 显隐设置参数对话栏：这是一个开关，点击它，参数设置栏在显示和隐藏之间变换以保证对图形观察的空间。新增或修改参数：单击新增或修改参数允许用户添加新的泥浆参数或仪器参数。当前只允许添加油基泥浆，新泥浆类型名一定以OBM三个字符开头。单击增加，在固定行标题下出现新空行，给出油基泥浆类型名称和参数，键入保存，确定，新的泥浆类型将被纳入参数设置的泥浆类型表，供设计中选择。参数对极化周期和扩散系数的影响由选定的公式1和公式2确定，这两个公式仅用于油基泥浆。在增加仪器参数前，要先增加仪器

5、类型，打开核磁模版目录中的tooltype.txt文件,在文件末尾字符串”New”的下边输入新的仪器类型及该仪器所能测的观测模式，然后单击仪器参数按钮，可增加新仪器的参数，增加后保存。 二、P型核磁共振解释软件操作手册Forward.NET平台P型核磁共振解释软件可以分别处理时间域（*.m.wis）和深度域(*.d.wis)的核磁共振测井数据。整个运行过程都是由“P型核磁共振成像分析”模块来完成。启动模块，选定WIS文件后，模块首先对WIS文件作检查，根据检查结果作对处理内容的自动定位 如果WIS文件是双TW双TE的时间域原始数据首先要进行文件拆分，将其拆分为单TW双TE或单TE双TW的新WI

6、S文件，运行的处理步骤为“文件分离”；如果WIS文件是单TW双TE或单TE双TW的时域数据，运行的处理步骤为“回波拆分”、“时域分析计算”和“时深转换”等；如果WIS文件是深度域数据则可运行的处理步骤为“岩石物理计算”、“与常规测井结合确定总孔隙度”、“与常规测井结合确定泥质束缚水饱和度”、“与常规测井结合计算含油饱和度”等。（一）、文件分离1、功能 将DTWE数据文件分解成DTW、DTE数据文件，每组数据文件包括三组回波串，以便于后续测井资料评价。当MRIL-Prime型核磁共振测井仪以DTWE模式工作时，采集的回波串包括A、B、PRO6（或C）、D、E五组，其中A、B两组是由短TE（1.2

7、 ms或0.9 ms）采集的双TW模式；D、E则是由长TE采集的双TW模式，所以A与D组合可以得到一个双TE观测，B与E是另一个双TE观测（只是由于它们都是短TW采集到的，通常不使用）。PR06组（或C组）用于粘土束缚水的观测，其采集参数是固定的，即TE0.6 ms，TW=0.2 s，NE=10，RA50。因此，五组回波串的分组方式有abc、adc、cde、bce四种。2、输入文件 原始的*.m.wis(时间域DTWE数据文件)3、输出文件 自动形成*.abc.m.wis文件或*.acd.m.wis文件或*.cde.m.wis文件或*.bce.m.wis文件，分别保存在选择的目录下的同名文件夹

8、内。其中 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

9、)点击开始处理，在弹出的“文件提取”界面选择需要得到的观测方式和存储目录；5)选择“确定”运行程序，自动生成与选择观测方式对应的文件夹。（二）、回波拆分1、功能 核磁测井原始资料是时域数据且原始回波按等时记录于RAMP, 回波拆分模块将记录回波按组拆分，同时完成对回波串累加和作相位校正；对回波串多指数拟合，得到离散的T2分布及各区间孔隙度；为TDA_COMP的计算提供长、短TW回波串的差谱EDIF；输出孔隙度MPHIX用作深度校正。注意,该程序不能对bce、cde组数据进行处理。2、输入曲线 GRP MRIL-P各采集组的编码CACT MRIL-P各采集回波串编码（0-39）CECH 采集设置

10、中的回波串个数RAMP ECHO幅度（%）RPHA ECHO相位角（度）TDEP 回波串采集对应的电缆深度码PWCO 功率校正因子SACO 矿化度校正因子TMCO 温度校正因子GAIN 输入增益3、输出曲线 A组校正后的回波串的幅度、相位、实部、虚部B组校正后的回波串的幅度、相位、实部、虚部各组各部分的离散T2值4、运行说明 运行方式与文件分离模块相似。本处理不需要参数文件，但为与其他处理框架一致，单击处理按钮弹出模式选择对话框。模块根据读出数据的回波个数自动搜索可能的模式型号，等待时间和回波间隔显示为空。用户选定本次测井具体的模式型号后等待时间和回波间隔会自动填入。模式型号选择后，确定相位校

11、正参数和解谱参数，如果用户需要改变这个值，可钩选该方框活化下面的两个输入框，然后键入新值，然后单击确定进入运行。相位角校正：将通道1和通道2的数据旋转一个角度，旋转之后，一道将主要包含NMR信号（其幅度是回波串的实数部分），另外一道将主要包含噪音（其幅度是回波串的虚数部分）, 回波串的实数部分被转换为T2分布。用下式计算角度这里i是回波串中的第i个回波，k是用做相位角计算的回波数目，模块隐含用第2到底9个回波计算。如果用户需要改变这个值，可钩选该方框活化下面的两个输入框，然后键入新值。平滑系数：滤波平滑系数。解谱首波：解谱中，时间轴上小于解谱首波的回波将被忽略。解谱末波：解谱中，时间轴上大于解

12、谱末波的回波将被忽略。（三）时深转换1、功能把时间域的*.m.wis文件转换成以深度域的*.d.wis文件。2、运行说明点击按钮弹出选择转换位置对话框 “将转换的深度域曲线写入新的WIS文件”，输入相应的目录及文件名，选择确定后开始转换；“将转换的深度域曲线写入存在的WIS文件”，保存文件对话框弹出后用户可选择硬盘上存在的WIS文件。此时模块不创建新文件，但所选择的WIS文件必须是深度域的文件，否则转换失败，你收到信息 “核磁仪器偏置深度”：该深度由测量过程中测量的延迟引起。“GR仪器偏置深度”：该深度由测量过程中测量的延迟引起。 在根据原始记录的时间、深度记录转换之后，将根据两者之差作深度移

13、动。(四)、岩石物理计算1、功能 计算核磁各种孔隙度；计算核磁渗透率；T2谱的合并，T2的0.5，1，2，4ms bin数据由TE=0.6ms的回波串拟合结果提供，8，2048ms数据由TE0.6ms的回波串拟合结果提供；由离散的T2分布得到连续的T2谱，（200个元素的向量）；计算油气的NMR特性（根据温度压力及理论公式）；对所有的孔隙度曲线进行滤波处理；对TDA-COMP的计算结果进行滤波。2、输入曲线ECHOSTRIP和PROCESS-T2D处理后的输出曲线。3、输出曲线各组的T2谱和各种孔隙度、渗透率流体核磁特性等曲线。4、输入参数 FILTFLG 滤波选择 0=NO,1=YES FL

14、TTYP 滤波方式0=汉宁窗,1=三角形,2=矩形窗 FLTLEN 孔隙度滤波长度（采样点数） T2CUTOFF T2 BVI 截止值 SBVICOEF 台阶状模型系数（SBVI） T2SPLICE T2谱拟合时前几个BIN从PR谱拼合 CCOEF 渗透率系数（COATES模型） CTHRES 渗透率计算时最小孔隙截止因子 TE 回波间隔 一般取0.9 5、运行说明 在做完回波拆分后应用本模块对回波串做时间域分析。在作岩石物理计算处理以前，处理合并带常规GR、孔隙度、电阻率资料，并将核磁的GAMM曲线与常规GR曲线对比，进行核磁测井资料校深；亦可直接作岩石物理计算处理，这样软件可以自动拟合得到

15、一条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）输出曲线T1APP T1的估算值T1ERR T1搜索误差最小值DEPTHT1 T1搜索的输出深度T2X1 单相搜索的T2PORX1 单相搜索的孔隙度T2ERRX 单相搜索误差最小值T2X 双相

16、搜索的短分量T2PORX 双相搜索的短分量孔隙度T2Y 双相搜索的长分量T2PORY 双相搜索的长分量孔隙度T2ERRXY 双相搜索误差最小值（4）输入参数TISEARCH 是否进行T1搜索，0不搜索，1搜索T2SEARCH 是否进行T2搜索，0不搜索，1搜索NUMECHO 回波差的回波个数Te 回波间隔TWB 记录一个完整的短Tw回波串记录的时间TWA 记录一个完整的长Tw回波串记录的时间TISTKW T1累加次数（40的倍数，次数越多信噪比越高）TIFBIN T1搜索的第一个binTIFBIN T1搜索的最后一个binMINTI T1的最小时间范围（范围越大搜索精度越高）MAXTI T1的

17、最大时间范围T2STKW T2累加次数（40的倍数）FECHO 搜索T2的第一个回波LECHO 搜索T2的最后一个回波T2XMIN 1相的T2最小值T2XMAX 1相的T2最大值T2YMIN 2相的T2最小值，T2XMAXT2YMINT2YMAX 2相的T2最大值MNPHIT2 计算长的最小孔隙度门槛值T2DELIA T2的搜索步长，一般取1.0T2EPS T2的搜索的误差容限T2MXITR T2的搜索的最多迭代次数T2COMPS 搜索几组T2值（单相或两相）2、时域分析计算（1）功能 时域分析（TDA）依赖于这样一个实事，即不同的流体有不同的极化率，或不同的T1驰豫时间。气和轻质油（粘度小于

18、5cp）通常比水的T1要长得多。时域分析提供 冲洗带流体类型；经含烃指数校正后的核磁孔隙度；各组分流体饱和度的分析（冲洗带）；进入TDA计算的数据已经过对双TW双TE的分离，经过了回波拆分。（2）输入曲线EDIF A组、B组回波串之差（3）输出曲线PHIGU 未作任何校正的气相孔隙体积PHIOU 未作任何校正的油相孔隙体积PHIWU 未作任何校正的水相孔隙体积PHIGA 未作HI校正的气相孔隙体积PHIOA 未作HI校正的油相孔隙体积PHIWA 未作HI校正的水相孔隙体积PHIG 气相孔隙体积PHIO 油相孔隙体积PHIW 水相孔隙体积（4）输入参数T1GAS 气的T1值 T2GAS 气的T2

19、值 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）点

20、击开始处理，完成T1、T2搜索；4）选择“时域分析计算”模块，编辑参数并处理，完成时域分析。 （6）流体T1、T2计算器说明 如果用户没有参数的经验值，可采用理论参数作初始处理。单击功能按钮 “计算流体T1、T2值”，弹出“计算核磁特性参数”对话框，这是一个输入参数的辅助工具，请注意对话框右侧的计算结果除气体密度外均为输入参数值，如果用户有经验参数可直接输入。如果需要本计算器来计算，可选定油藏类型、测井仪器、地层温度、地层压力、油的粘度、扩散系数、气体密度。输入参数栏说明 油藏类型 砂岩层（按渗透率分类）、碳酸盐岩。测井仪器 MRIL6” MRIL4，7-8”选择计算气体密度方法 气体密度由两

21、种方法输入， “用户输入”将使用用户在输入参数栏气体密度输入值计算T1、T2。“由温度压力计算”将根据用用户在输入参数栏输入的温度压力计算气体密度。 按下后，计算出右侧的所有参数。在精细解释深度棒上选定一个层，此按钮被活化，按下后将计算结果写入该层段参数。按下后将计算结果写入所有层段参数。按下后，所有输入状态被保留，退出计算核磁特性参数对话框，下一次打开计算核磁特性参数时保留参数被显示。如果用对话框的关闭X，则输入状态不保留。(六) 扩散分析在DTE测井模式下，可根据TEL下最大水线T2DW右侧谱峰的有无、大小识别评价油气层。扩散分析计算分为确定可动流体T2、扩散分析预处理和计算含水饱和度三部

22、完成的。1、确定可动流体T2（1）功能 对T2谱进行分析，计算出长、段TE的几何平均值、峰值、半峰值、为扩散分析计算提供数据。（2）输入曲线 TASPEC 短TE的T2谱TBSPEC 长TE的T2谱（3）输出曲线 T2GMS 短TE组的T2几何平均值T2DIF1S 短TE组孔隙度峰值对应的T2T2DIFS 短TE组孔隙度半幅点的T2T2S 用于扩散分析短TE组的T2T2GML 长TE组的T2几何平均值T2DIF1L 长TE组孔隙度峰值对应的T2T2DIFL 长TE组孔隙度半幅点的T2T2L 用于扩散分析长TE组的T2（4）输入参数 AMINA 用于短TE组确定峰值的最小孔隙度 AMINB 用于

23、长TE组确定峰值的最小孔隙度 T2GMLA 短TE组计算几何平均T2GMS的下限T2GMUA 短TE组计算几何平均T2GMS的上限T2GMLB 长TE组计算几何平均T2GML的下限T2GMUB 长TE组计算几何平均T2GML的上限PMINA 短TE组计算几何平均T2GMS的最小孔隙度界限PMINB 长TE组计算几何平均T2GML的最小孔隙度界限T2SELFLG 扩散分析计算时使用哪组T2,1-T2GM 2-T2DIF1 3-T2DIF2、扩散分析预处理（1）功能 计算流体扩散系数，为扩散分析计算提供数据（2）输入曲线 T2S 用于扩散分析的短TE的T2值T2L 用于扩散分析的长TE的T2值TB

24、SPEC 长TE的T2谱（3）输出曲线 T2INT 固有T2值T2DW 水的T2界限(长TE)D 地层流体的扩散系数DW 水的扩散系数RDDW D与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.

25、875(现场提供)T2LIM 固有T2的上限(0表示无限制)GFACT 磁场剃度调节因子(现场提供)T2SFACT 短T2调节因子3、扩散分析计算含水饱和度（1）功能 对双TE（A组为短TES，B组为长TEL）做移谱分析。（2）输入曲线 MPHITA 有效孔隙度(短TE)MSIGTA 有效孔隙度(短TE)MBVITA 毛管束缚水孔隙度(短TE)MPERM Coates模型渗透率T2S 用于扩散分析的短TE的T2值T2L 用于扩散分析的长TE的T2值T2INT 固有T2值IT2 固有T2值的倒数RDDW D与DW比值（3）输出曲线 DIFSW 水饱和度DIFBVW 水的孔隙体积DIFKRO 油的

26、相对渗透率DIFKRW 水的相对渗透率DIFWCUT 持水率估计值（4）输入参数T2MIN 自由流体的最小固有T2T2HY 烃的固有T2RDDWHY 烃与水的扩散系数的比值ROSF 残余油饱和度(%)T2SFACT 短T2调节因子RELPFL 相对渗透率模型,1-米那斯方法 2-普通方法 3-琼斯方程 4-彼尔逊方程 WVIS 地层水粘度(cp)OVIS 地层油的粘度(cp)（5）运行说明 为确定各种流体参数，可使用扩散分析图版进行分析。单击功能按钮 “显示确定扩散分析图版”。图版显示到文档中，双击图版，将X轴设置为RDDW（流体扩散度与地层水扩散度的比值），将Y轴设置为IT2(流体固有T2值

27、的倒数)，将“能够编辑图形”键关闭，图中兰色虚线为视水线标尺，上边是100%视水线标尺，下边是0%视水线标尺，移动光标到视水线标尺处，当光标变为上下箭头时，上下移动光标，可拖动视水线标尺，移动时程序自动控制两条标尺保持平行，参数框中的参数也随之改变，单击流体参数框中的写入选中层字段可将确定的参数写入到选中的层段中，单击写入所有层字段可将确定的参数写入到所有层段中。（七）标准T2测井处理(MRIAN)主要是利用核磁共振测井资料与常规的孔隙度、电阻率资料结合进行综合分析、提供准确的储层物性参数、原状地层含油饱和度计算。包括三个解释模块，分别为与常规测井结合求总孔隙度模块、与常规测井结合求泥质束缚水

28、饱和度模块和与常规测井结合求含油饱和度模块。1、与常规测井结合求总孔隙度（1）功能 准确求取总孔隙度t，计算从核磁得到的Swb，提供后续处理所需的各种交会图数据。（2）输入曲线 MSIGTA 总孔隙度(由深度域曲线得到) MPHITA 有效孔隙度(由深度域曲线得到) MBVITA BVI值 TDAMSIG 总孔隙度(由时间域曲线得到) TDAMPHI 有效孔隙度(由时间域曲线得到) MCBW 粘土束缚水孔隙度 NPHI 补偿中子 RHOB 体积密度 DT 声波时差 RT 深探测电阻率（3）输出曲线 FTEMP 地层温度RWAT 地层温度下的地层水电阻率RMFAT 地层温度下的泥浆滤液电阻率SA

29、LW 地层水矿化度(来自电阻率-矿化度图版)SALMF 泥浆滤液矿化度(来自电阻率-矿化度图版)RHOW 地层水密度PHID 计算的密度孔隙度PHIN 计算的中子孔隙度PHIS 计算的声波孔隙度PHIX 中子-密度交会孔隙度TPOR 地层总孔隙度RWA 视地层水电阻率CWA 视地层水导电率SWBMRI MPHITA确定的束缚水饱和度DPND 中子与密度孔隙度的差DPTD 声波与密度孔隙度的差RPND 中子与密度孔隙度的比值DPDM 密度与MRIL孔隙度的差（4）输入参数 TFLAF 温度单位 1=华氏度 2=摄氏度SURFT 地表温度 BHT 井底温度 TD 井深(米)PRESSG 地层压力梯

30、度(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 阿尔奇AEXPHI 阿尔奇MRHOFL 流体密度DMA 骨架视密度DTMA 骨架时差SONICFL 计算声波孔隙度的公式, 1= Raiga-Clemenceau 2=Chapman-Modified2、与常规测井结合求泥质束缚水饱和度（1）功能 计算除核磁测井以外所有其它方法计算的Swb，并将Swb与核磁的Swb进行比较，选择一个最佳的Swb。在进行Swb-MRIAN解释前要对前一步骤解释结果作一些交会图