P型核磁共振测井原理及应用

1、0 p 型核磁共振测井推广应用型核磁共振测井推广应用 目 录前前 言言.1一、基本测井原理及仪器介绍一、基本测井原理及仪器介绍.11、核磁的测井原理及仪器介绍.12、质量控制.4二、二、p 型核磁的处理解释方法型核磁的处理解释方法.61、处理流程.62.时域分析 （tda）.73.扩散分析 （difan）.9三、对三、对 p 型核磁测井处理系统的改进型核磁测井处理系统的改进.111、builder 开发工具包.112、解释系统的二次开发及应用.12四、实际应用四、实际应用.141前前 言言1945 年斯坦福大学的 bloch 教授和哈佛大学的 purcell 教授领导的两个小组相继发现了核磁共

2、振现象。从此开创了一门新的学科核磁共振波谱学，经过了近年的探索和发展，核磁共振已广泛应用于化学、石油、地学、食品和医学等领域。年，公司利用概念，利用梯度磁场和自旋回波方法，设计开发了全新的磁共振成像测井（） ，并于年月正式投入油田商业服务，斯仑贝谢公司于年把以贴井壁磁体为核心的核工业磁测井（）推向商业服务。哈里伯顿公司收购公司以后，逐步推出了一系列核磁共振测井仪器，如：挂接在excell2000 下的 p 型核磁共振仪（mril-prime） ，地层测试器 rdt 发展出称为核磁实验室的模块（mri lab）以及随钻核磁共振测井仪（mril-wd tool） 。核磁共振测井可以提供十分丰富的地

3、层信息，能够定量确定有效孔隙度，自由流体孔隙度、束缚水孔隙度、孔径分布以及渗透率等参数。核磁测井的应用范围在不断扩大，在油田注水开发过程中，可用于确定油层水淹程度、驱替效率、剩余油饱和度、产层性质、可采储量及采收率等。在复杂岩性碳酸盐岩、火成岩储层，在低孔低渗低电阻储层中，寻找气层，区分油、气界面。在裂缝性油藏的综合评价，提供强有力的信息。一、基本测井原理一、基本测井原理及仪器介绍及仪器介绍1、核磁的测井原理及仪器介绍、核磁的测井原理及仪器介绍2000 年我公司引进了哈里伯顿的 mril-prime 核磁测井仪器，与阿特拉斯公司的 c 型核磁共振仪器同属于 numar 公司的产品，采用与 c

4、型仪器相同的测井原理，都属于居中测量，采用的也是 cpmg 脉冲序列，测量方式上也采用标准 t2、差谱、移谱测井，因为可以用 0.6ms 回波间隔进行测井，能测量地层粘土孔隙，所以可以得到地层的总孔隙度。mril-prime 型仪器与 mril-c 型结构是一样，仪器探头由三个基本部分组成：一个永久磁铁，一个射频脉冲（rf）发射器及一个射频接收器，但 mril-prime 仪器的永久磁铁较 c 型加长，增强了极化能力。与 c 型仪器相比具有更多的优势。p 型核磁共振测井仪采用 9 个观测频率，它在地层中形成 9 个厚度为 1mm 间隔 1mm 高度 24in2的探测圆环，使得仪器对地层的探测体

5、积增大，得到的地下信息更准确，更丰富。9 个频率分为 5个频带 0、1、2、3、4，中心频率分别为 590、620、650、680、760khz，最后一个频带（760khz）只用来进行粘土束缚体积测量。频率的增加，意味着在相同的时间内，它比单频或双频仪器采集更多的信息，换个方式说，如果要采集相同的信息，所用时间要少 3 到 4 倍，这是因为其中一个频率在等待的时候，其它几个频带可以进行数据的采集，大大提高了测井的速度。见图 6。p 型核磁共振测井分为单 tw/单 te、双 tw/单 te、单 tw/双 te、双 tw/双 te 四种采集次序，每张图都显示了一个完整的采集周期，分别如图 14 所

6、示：图 2、双 tw/单 te 采集次序3 单 tw/单 te：a-a-a-a-pr-a-a-a-a-pr双 tw/单 te：a-a-b-b-pr-a-a-b-b-pr（a：长等待时间，b：短等待时间）单 tw/双 te：a-a-pr-a-a-pr-a-a-pr-a-a-pr双 tw/双 te：a-a-b-b-pr-d-d-e-e-pr- a-a-b-b-pr-d-d-e-e-pr，其中 a 与 d 为相同的长等待时间，b 与 e 有相同的短等待时间，但回波间隔不同，b 为短回波间隔，e 为长回波间隔，所以可以将之分解为两个双 tw/单 te 和一个单 tw/双 te。如图 5。 图 1：单

7、tw/单 te 采集次序图 3：单 tw/双 te 采集次序图 4：双 tw/双 te 采集次序4abc：双 tw/短 te+prdec：双 tw/长 te+pradc：双 te+pr由于可以进行多个频率采集，p 型核磁共振测井设计了很多观测模式，在测井施工以前，根据地区和油气藏特征，通过测前设计，选用所需要的观测模式。其中 tw/单 te 模式分为 6 种；双tw/单 te 共有 14 种；单 tw/双 te 模式有 24 种；双 tw/双 te 模式有 8 种。模式的划分是根据等待时间和回波间隔的不同，短等待时间最小为 1.0s，长等待时间最长为 15s，回波间隔最小为0.9ms，最长为

8、6.0ms，每种模式都有回波间隔为 0.6ms 的 10 个回波，对粘土束缚水进行测量。根据 c 型核磁共振测井应用的经验及习惯，推荐选用以下模式：单 tw/单 te 模式：d9tp12（tw：12s，te：0.9ms，echo：500）双 tw/单 te 模式：d9tw（长 tw：12.98s，短 tw：1.0s，te：0.9ms，echo：500）单 tw/双 te 模式：d9te310（短 te：0.9ms ，长 te：3.6ms ，tw：10s，echo：500、125） d9te412（短 te：0.9ms ，长 te：4.5ms ，tw：10s，echo：500、100）双 tw/

9、双 te 模式：d9twe3（te：0.9、3.6ms ，tw：12.98、1.0s，echo：500、125） d9twe4（te：0.9、4.5ms ，tw：12.98、1.0s，echo：500、100）mril-prime 仪器的主要技术参数：仪器探头外径：47/8和 6仪器探测范围：6-8.5(47/8探头)图 5：双 tw/双 te4abc、dec、adc 图 6：mril-prime 仪器示意图57-16（6探头）泥浆电阻率：大于 0.02m最大工作温度：350f（约等于 176 摄氏度）磁共振频率：500khz-800khz2、质量控制、质量控制 图 7 ：p 型核磁共振测井控

10、制面板（a 组）如图 7 所示，p 型核磁采用了较为完善的质量监控手段，可以分为以下几个方面增益：增益是仪器对周围环境的感应，gain 值应比较稳定，a、b、c、d、e 五组增益值应比较相近，增益值一般受泥浆电阻率和井眼温度影响比较大，一般随泥浆电阻率的减小和温度的增高而减小，如果增益值减小到 100 以下，仪器将无法工作。 b1 和 b1mod：b1 是 cpmg 脉冲的强度，它使 90 度的质子产生 180 度偏转，这个值是已被刻度好的，这个值必须做温度校正，下井后操作员调整 b1a、b1b 使 b1moda、b1modb 等于 b1，误差应在5%以内。 chi 值：多指数拟合程度，这个值

11、应小于 2，如果信噪比低，chi 值可以达到 3。 噪声指示：（offset、noise、ringing、ienoise）noise、ienoise 应小于10，offset、ringing 应大于-30 而小于 30。6探针电压：相位校正信息：pher（虚部均值）：应等于零；phno（回波串虚部标准偏差）：约等于1；phco（相位角）：在测量井段内应比较稳定。 图 8：p 型核磁共振测井质量控制图 8 包括了不同组和不同频率的 gain、b1、b1mog、探针电压、回波串相位特性、噪声指示、chi 以及测量得到的有效孔隙度和束缚孔隙度。如果任何一个质量控制指示值超出它所允许的范围，它所对应的

12、颜色就会变暗。7二、二、p 型核磁的处理解释方法型核磁的处理解释方法1、处理流程、处理流程 图 9：mril-p 核磁共振测井资料处理流程p 型核磁处理流程如上图所示。下面对处理项目进行简单的介绍split_mcls：将双 tw 双 te（max2twte）测井文件分解，一般拆成三个文件，分别为 2 个回波间隔不同的 dtw 和 1 个 dte 文件。reseq：将前面拆开的文件分别重写代码。echostrip：运用标准模块处理所有回波。经过叠加平均和相位校正把回波转为孔隙度（bin） ，如果资料信噪比低，可以增加回波串的叠加个数（8 的倍数） ，但这就要牺牲纵向分辨率，叠加次数越多，纵向分辨

13、率就越差。在时间域做 t2 拼接：pr 组为 0.5、1、2、4ms，其它组为 8 到2048ms，得到总孔隙度，这个孔隙度没有进行余弦滤波。t2t1_event:估算自由流体的 t1 以及最多两相流体的 t2。tda_comp：这个模块的基本思想是认为差谱后的信号都是油气的信号，从而得出油、气、水三相流体的孔隙度。双 tw 双 te（max2twte）：split_mcls、reseq，然后按不同测量方式进行处理。标准标准t2t2：echostrip、processt2d、t2_toolkit、pre_mrian、swb_mrian、mrian。dtwdtw：echostrip、t2t1_e

14、vent、tda_comp、processt2d、t2_toolkit、pre_mrian、swb_mrian、tda/mrian。dtedte：echostrip、processt2d、t2_toolkit、pre_mrian、swb_mrian、t2_difsel、difan、difan/mrian。可不做 tda 和 difan 分析。8processt2d：时间域转换为深度域，以下的步骤全部在深度域进行。t2_toolkit：为后续处理做各种相关的计算。包括孔隙度（bin）的滤波、粘土束缚水、渗透率、总孔隙度的计算以及解谱。完成 t2_toolkit 以后，对资料进行综合分析，需要输入

15、常规资料，电阻率、中子、密度曲线是必须的。综合分析模块总称为 mrian，分为三个步骤：pre_mrian、swb_mrian、mrian。pre_mrian：完善总孔隙度，计算核磁共振的束缚水含量 swbmri，为交会图准备数据。swb_mrian：更准确的计算束缚水含量。mrian：计算输出最终处理解释结果。t2_difsel：通过 te 分布范围，计算长短回波间隔的几何平均值，为 difan 输出参数。difan：扩散分析，也可称为稠油分析。p 型核磁数据可以单独处理分析或与常规数据结合分析，当单独从核磁数据解释，它可以提供侵入带孔隙、渗透率、流体类型和流体饱和度，p 型核磁提供了两种单

16、独分析数据的计算模块，在时间域里运行的时域分析（tda）模块和在深度域里运行的扩散分析（difan）模块。 2.时域分析时域分析 （tda）时域分析使用的是由双 tw 采集的数据，进行油气识别，包括两个步骤：t2t1_event、tda_comp，它能有效的从测量数据中搜索出自由流体的 t1 以及最多两相流体的 t2，来得到油、气、水的孔隙体积。它的基础是不同的流体有不同的极化速率，或不同的 t1 弛豫时间，气和轻质油的 t1（粘度小于 5cp）通常比水的 t1 长，时域分析可以提供以下成果（1）冲洗带的流体类型（2）含气储层校正的 mril 孔隙度（如果不进行校正，由于气比较长的 t1 和低

17、氢指数，使测得的孔隙度偏小）（3）对轻质油校正的 mril 孔隙度（4）完成冲洗带的流体饱和度分析tda 是从 dsm(differential spectrum method)发展起来的，差谱主要是定性检测气层的存在，基本原理如图 11 所示，tda 分析是在时间域把长等待回波串与短等待回波串相减，形成差谱（必须大于 1.5cp） ，而不是在 t2 域相减，因此它比 dsm 有两个关键优势在时间域计算的两个回波串的差更可靠，这个差谱再转换成 t2 分布tda 对未完全极化和氢指示影响能进行较好的校正9图 10 是 tda 分析的 t1 搜索基本原理t1fbin、t1lbin：t1 搜索的起始

18、、结束 bin。twlu、twsu：在给定的搜索范围内，长短等待时间的孔隙度。由核磁共振孔隙度公式我们知道：twlu=*hi*（1-e*（-tw/t1） ）所以可以用下面公式计算出自由流体的 t1：twlu/twsu=（*hi*（1-e*（-twl/t1） ） ）/（*hi*（1-e*（-tws/t1） ） ）t2 搜索：tda 分析认为差谱后的信号都是油气的信号，在时间域把长等待回波串与短等待回波串相减，形成差谱 edif，见图 12：如果认为差谱由两相流体组成（油、气） ，采用的是双指数拟合：edif=x *e*（-t/t2x）+y *e*（-t/t2y）如果认为只有一相流体（油或气） ，

19、则采用单指数拟合：edif=x1 *e*（-t/t2x1）通过给出的流体相数目、起止回波、迭代限制等参数搜索得到两相流体的 t2x 和 t2y 或单相流体的 t2x1。图 11： dsm 差谱示意图 图 12： tda 差谱示意图 以两相流体为例，前面得到的 t2x 和 t2y 分别设定为 t2gas、t2oil，t1 设为 t1oil。从而拟合得到油、气的孔隙度。经过含氢指数和 t1 校正，最后得到 tda 分析的油、气、水孔隙体积，这 图 10：t1 搜索示意图10些孔隙体积反映的应该是冲洗带的流体。如果已知施工地区油的粘度、温度，可以计算一个理论的 t2，与搜索得到的 t2 进行对比，如

20、果两者吻合，就可以认定为油。图 13 是 tda 分析的成果图，第 6 道是 tda 处理的流体孔隙度结构，第 7 道是核磁共振结合常规资料处理得到的流体孔隙结构。3.扩散分析扩散分析 （difan） 扩散分析采用移谱测井数据，用于中等粘度的油品的识别与计算。不同流体的扩散性质不同，是由于梯度磁场的不均匀性，导致 t2 随回波间隔 te 改变，其关系可以用方程表示： 12/)(21212intgrtecdtta g 磁场强度 r 氢核的旋磁比da 分子的视扩散系数 te 回波间隔t2int 当磁场强度为 0 时固有的弛豫时间c 经验常数，对核磁为 1.08方程中除了回波间隔 te，所有的都是常

21、量。移谱方法（ssm）是利用不同的回波间隔，不同的流体在 t2 上的变化不同的一种定性识别流体的方法。difan 是定量分析流体扩散性质的一个经验模式，它可以搜索出流体的长、短回波间隔的 t2 值 t2l、t2s。它主要解决 t1 差别不太明显，不易使用 tda，或扩散性差别很小，不易使用双 te（ssm）的情况，它通过搜索得到不同流体在长短回波间隔的 t2 分布上位置的不同，定量计算孔隙中含水体积和含油体积。对于低粘度的碳氢化合物如轻质油和凝析油，不适用 difan，因为油水的扩散系数相差很小，对于粘度大于 100cp 的高粘度的油（重油） ，也不适用 difan，因为残余油的 t2int

22、值和不可动水的 t2int 值相差很小。图 14 是扩散分析处理结果，第 6 道是经过difan 处理得到的不同流体的孔隙体积，第 7 道是核磁共振结合常规资料处理得到的流体孔隙结构。11 图 13：tda、mrian 处理成果图12 图 14：difan、mrian 处理成果图三、对三、对 p 型核磁测井处理系统的改进型核磁测井处理系统的改进国内多家测井公司在引进哈里伯顿公司的 mril-p 型核磁采集系统的同时，都相应地引进了其解释系统 dpp。其核磁测井资料的谱差分测井和谱位移定量计算为我们提供了较准确的冲洗带饱和度参数。通常，在不同的地区，采用与地质情况相符合的含油饱和度模型才能获得准

23、确的结果。但是，dpp 软件在含水饱和度模型中，仅提供双水模型，且其粘土水的附加导电性是一温度的函数。事实上，由于同期地层在不同地区埋深不同（温度不同）或沉积环境的差异均会导致粘土水的电导率变化。而且，在成果图上无法标注符合国内解释标准的解释结论，因此，它不能满足生产中的需要，鉴于这种情况，利用其中的开发工具包 builder 进行二次开发，完善了双水饱和度模型，并填加了适合于不同地区的饱和度模型；增加了标注符合国内解释标准的各种结论等，解决了引进软件在后处理上不适应国内生产需要。1、builder 开发工具包开发工具包哈里伯顿公司的石油物理解释平台为用户提供了功能较强大的用户开发工具包，它为

24、用户提供了两种开发语言（fortran 和 c）,并为用户建立模型提供了良好的输入、输出接口，它要求用户指定输入、输出曲线、参数，通过用户程序的计算输出各种处理结果，并为用户提供一标准的dpp 运行界面。完全能满足在实际生产中二次开发的需要。其主要的开发流程如图 15 所示，其中的解释模型参数申明文件（mds）为指定测井解释程序的输入曲线、输出曲线和各种参数，以及申明在解释程序中需要的各种变量。而模型使用的语言文件(mdl 文件)可以采用 c 语言，也可以采用 fortran 语言进行编写，编译连接就可以得到用户需要的各种解释程序。mds 文件的包含较多的对资料处理方法和处理参数进行予定义模式

25、，但应用较多的包含以下部分：（1）、测井数据输入输出，表示为logs|曲线名|指定输入输出标记|曲线单位|指定曲线数据类型|曲线描述|指示特殊处理标记|其中在特殊处理的标记中可以指定处理矢量数据的大小(sz)和处理数据的采样率(step)。（2）、解释模型参数输入，表示为：13cnst|参数名|缺省值|输入或输出|最小值|最大值|单位|参数描述|参数组序号|特殊处理标记|其中的特殊处理标记为是否定义最大和最小值，或采用处理的数据文件中对应的最大、最小值。（3）、参数组序号grpn|组名|序号|它的重要作用是在运行解释系统时根据不同的组序号需要选择不同的解释参数。（4）、运行程序时需要的变量申明

26、darr|变量名|数据类型|变量维数大小|描述|内部名标记|在指明各种输入、输出以及参数后，需要写一个数据处理模块，其处理模块mdl文件为一标准的fortran或c语言的处理程序，利用开发工具包中的builder创建的fortran解释程序结构如下：if (mdlint) then 处理初始read (topwl, (f12.4) welltp 读入数据文件顶深read (botwl, (f12.4) wellbt 读入数据文件底深 加入用户代码 end ifif (zonint) then 是否是需要处理的深度段pass = 0 在多步处理中使用. 加入用户代码endif 由 mdl 编写的

27、解释模型代码if (lstzon) then 本段程序为判断深度段处理是否完成代码mypass = mypass + 1pass = mypassxindep = xznstrsetcur = .true.newpas = .true. (只对多步处理有效)xotdep = absentendifif (lstmdl) then 本段为判断是否为最后的模型处理 加入代码endif1000 continueend在解释程序的编写中，必须要注意的是，由于 dpp 系统中可以在 mds 文件中选择处理顺序可以从井底井口，也可以从井口井底。因此，在处理程序中要采用深度控制的处理程序一定要注意这一点，否

28、则在一些情况下将得不到正确的结果。图 15、builder 开发示意图142、解释系统的二次开发及应用、解释系统的二次开发及应用利用开发工具包，在 dpp 系统中，改进了 dpp 系统中的 mrain 的双水饱和度模型和增加了ariche、waxman-smith、patchett-herrick 饱和度解释模型。下面主要对改进的双水模型的进行介绍。在哈里伯顿核磁综合解释模块 mrain 中，采用的双水模型计算饱和度方法如下： （1）)cc(ssc)s(cwcwwtwbwwwttt其中：ct-电导率，a-曲折度系数，t-总孔隙度，swt-总含水饱和度，ccw-粘土水电导率， cw-地层水电导率

29、，swb-粘土水饱和度，w-是系数其中的系数 w 是coates 和 howard 在 1992 提出的与有效束缚水饱和度有关的系数，表示为：w=0.4swire+1.65 其中:swire= （2）mphimbviswire而粘土水的电导率采用的是一与温度有关的函数，函数为：)7 .16t()4 .504t(000216. 0ccw其中：t-为温度，单位为华氏度。可以看出，在粘土水的计算中，粘土水的导电能量只与温度有关；事实上，相同温度下的地层，由于沉积条件和粘土矿物的不同，它们的粘土水的电阻率有较大的差异，它并不只与地层的温度有关。通常情况下，利用 p 型核磁测井的总孔隙度测井模式可以得到地层的粘土束缚水，而对于区域内压实好的无毛