核磁共振测井井下模拟器设计＊

张嘉伟 宋公仆 赛 芳 薛志波

(中海油田服务股份有限公司油技事业部 北京)

核磁共振测井井下模拟器设计＊

张嘉伟 宋公仆 赛 芳 薛志波

(中海油田服务股份有限公司油技事业部 北京)

核磁共振目前作为一种高端技术已经广泛应用于实际地层测试中。文章主要介绍了一种核磁共振测井仪井下模拟器的相关设计。通过该模拟器可以生成核磁共振测井仪在实际地层中所测得的数据,并按照一定的通讯协议格式传输到地面处理系统进行相应的分析处理。通过该井下模拟器,地面控制处理系统可以在没有实际核磁共振仪器的情况下很方便地进行地面软件的调试与处理算法的测试,从而极大地提高了地面系统软件编写与前期测试的效率。

核磁共振;回波信号;T2谱反演;孔隙度;渗透率;扩散系数

0 引 言

核磁共振测井方法可直接测量地层孔隙中可动流体的信息,可定量确定自由流体、束缚水、渗透率及孔径分布,其孔隙测量不受岩石骨架矿物成分的影响,因此目前颇受测井行业的广泛应用[1、2]。本文主要介绍一种核磁共振测井仪井下模拟器的相关设计。主要描述了模拟器的整体设计构架、各个关键模块的功能简介、工作流程设计程序结构,并详细说明了仿真数据生成的模型。最后对该模拟器的实际应用进行简要介绍。目前该模拟器已经应用于实际核磁共振地面系统的调试中,并取得了很好的运用效果。

1 核磁共振井下模拟器设计

1.1 井下模拟器设计方案

核磁井下仪器模拟器主要由三部分构成,包括数据传输短节简化部分、通讯电路和测井数据容量增强模块(以下简称数据增强模块)。通讯电路与遥传短节之间进行交互,通讯内容包括命令接收和数据上传;数据增强模块则用于产生与仪器采集相类似的数据,提供给通讯与遥传交互。如图1所示。

图1 核磁井下仪器模拟器功能框图

模拟器主要完成以下工作：通讯板接收传输短节发送来的曼切斯特码,进行解码;通讯板接收数据增强模块生成的数据,并将数据编成曼切斯特码,再通过不同的模式发送给传输短节;通讯板对传输短节发送来的命令解码后,下传给数据增强模块;通讯板产生基本模拟测井数据;数据增强模块通过软件建立孔隙度模型、渗透率模型等数据的算法,并在其中加入噪声及振铃等干扰,产生扩展模拟测井数据;数据增强模块在接收到数据请求命令后,将符合要求的模拟数据上传给通讯板;为操作人员提供方便的人机接口。

1.1.1 核磁井下仪器通讯电路的结构

通讯电路通过测井仪器总线实现与地面系统的数据和命令的交互功能,其典型工作流程是将模式时序表信息下发到事件控制电路,将采集到的辅助参数信息和采集处理电路获得的回波信息上传给地面系统。该电路主要由辅助参数测量、板间通讯模块(与事件控制采集电路通讯)、与地面系统通讯模块等组成,如图2所示。

图2 模拟器通讯电路功能框图

图3 仪器模拟器通讯板硬件构架

按照通讯板的功能要求,在图3中描述了通讯板的硬件构架。由图3中可以看出,通讯板主要由DSP+FPGA构成,其中的DSP负责辅助参量的数据采集、与数据增强模块的CAN通讯控制、与测井仪器的通讯接口的控制等;FPGA则主要是实现模式2、模式5和模式7传输方式(EDIB总线协议定义的模式)的数据编码和模式2传输方式的命令解码,同时还完成对扩展CAN控制器SJA1000的逻辑控制[3]。1.1.2 模拟器数据增强模块说明

如图1所示,模拟器中数据增强模块是其核心部分也是与仪器主控通讯电路的最大区别所在。数据增强模块存在的主要原因是由上传回波数据的类型和容量要求所决定的。测井回波数据的来源有三类,分别是：①对岩性分析仪静态数据进行动态化处理获得的数据;②根据设定的测井模型参数和NMR原理计算出的回波数据,这类数据需要预设模型参数和噪声参数,并进行动态化处理;③实际测井数据。

仪器主控通讯电路受到功耗、高温环境和电路板尺寸的限制,其计算速度和存储资源有限,而这三类数据无论从计算速度和存储资源看,都超出了通讯电路的限度,因此不可能仅由仪器主控电路来完全提供如上三类测井回波数据,数据容量增强模块即担负着提供回波数据这一功能[4]。

1.2 模拟器工作流程

核磁井下仪器模拟器是对井下仪器的工作的模拟,其基本功能是产生模拟测井数据传送到地面系统,从这一点看来井下仪器的工作流程基本上就是模拟器的工作流程。图4描述了模拟器工作的基本流程。

首先模拟器在上电伊始需进行相应的初始化操作,包括复位缓存,初始化完成后,等待接收地面系统的命令。若是收到参数下载或参数修改命令,模拟器则将参数表存储或更新;若是收到扫频或刻度命令,模拟器则产生数据后直接上传(根据现有协议不需要数据传输命令);若是收到采集命令,则由通讯板产生基本测井模拟数据或由增强板产生扩展测井模拟数据,等待发送;若是收到数据传输命令,则将准备好的数据上传。

图4 核磁井下仪器模拟器工作的基本流程

1.3 通讯电路与数据增强模块通讯

由图4的工作流程可知,通讯板在收到地面系统下发的命令及参数后,通过CAN总线将命令及参数发送给工控机,待工控机将数据生成并准备好上传,也通过CAN总线将准备好的数据上传给通讯板,通讯板则将数据编码为曼彻斯特码后上传给地面系统。因此,通讯板与数据增强模块之间的通讯是通过CAN总线来实现的,TMS320F2812片内自带一个CAN控制器,只需扩展隔离和驱动器即可。CAN总线的最大传输速率为1 M,通过计算可得实际传输速率为530 kbps,而M5和M7传输通道(EDIB总线协议定义的通道类别M5、M7,都以标准差分不归零曼码进行传输)的传输速率仅为93.75 kbps,数据从工控机传输到通讯板的速率远大于通讯板通过EDIB总线上传的速率,而通讯板上扩充了512K×16位的RAM作为存储缓冲区,能够存储工控机上传的数据并且不产生上溢,可以满足要求[5]。仪器要求的最小 Te(回波间隔时间)为0.4 ms,这里按0.1 ms的 Te来计算,需要传输的一个有用回波信号数据量为两个字,一个字幅值和一个字相位,由此计算出工控机与通讯板之间要求的传输速率最小为：2×16bit/0.1 ms=320 kbit/s=20 K words/s。

1.4 模拟器硬件结构设计

图5 模拟器的硬件构架

依据前面所述模拟器的功能和工作流程设计了图5所示的模拟器的硬件构架。地面系统与模拟器间的通讯包括命令下发和数据上传,均通过EDIB总线(一种有线电缆测井总线协议)来实现。通讯板中,EDIB总线协议在FPGA中实现,完成编码、解码工作,并通过变压器和收发器实现与 EDIB接口的匹配连接。DSP协同FPGA一起完成对DC电源和开关的控制工作,并提供CAN接口与数据增强模块通讯。数据增强模块的硬件平台是基于PC机结构,液晶触摸屏和键盘为用户提供了良好的人机界面,CAN接口使增强模块与通讯板可通过CAN协议传输数据和命令。

2 模拟器仿真数据生成模型设计

模拟器有效仿真数据生成主要是在数据增强模块里实现的,它相当于井下模拟器的仿真信号数据库。根据地面系统下发给模拟器的相关指令参数,生成相应的仿真回波信号数据并通过CAN总线传输给通讯电路,由通讯电路将其转换成 EDIB总线对应的M5、M7通道数据类型后传输给地面系统进行处理和分析。

2.1 仿真数据生成模型上位机设计

模拟器仿真数据生成模型上位机软件主要实现两个方面的功能：(1)与地面系统之间的通信,接收地面系统的指令并根据不同的指令实现不同的功能,生成相应的数据回传至地面系统;(2)上位机软件还需要实现扫频、刻度以及回波数据的生成,修改,导入,导出等功能。其主要程序模块介绍如下：①CAN驱动程序：模拟器通讯电路与上位机通过CAN接口相连,采用PCA82C250芯片作为CAN接口控制芯片实现模拟器通讯电路与数据增强模块之间的通讯交互。②上位机软件的线程设计：地面系统通过发送指令,模拟器仿真数据生成模型上位机软件接收到指令后开始仿真,并将生产的数据回传至地面系统。地面系统与上位机软件之间的通信协议简图如图6所示。③回波数据的生成：回波数据的生成是重点,也是难点。回波数据的生成需要根据核磁共振的原理来建立一个核磁模型,输入不同的对象和不同的参数(例如回波的个数值Ne、回波的间隔时间 Te,等待时间TW等)来进行仿真生成不同对象的回波数据,并回传至地面系统。重点是需要建立核磁共振的模型,对模型的仿真存在着不同的观测模式,产生的回波数据根据不同的观测模式有不同的数据格式。

图6 上位机软件通信协议简图

2.2 仿真数据生成模型

仿真数据的生成模型主要分为：孔隙度模型和渗透率模型两种。上位机软件通过对这两种模型的生成公式计算出相应孔隙度与渗透率对应的回波数据,并回传给地面系统进行处理与分析。

2.2.1 孔隙度模型设计

核磁共振测井确定地层孔隙度的依据来自观测信号强度与孔隙流体中氢核含量的对应关系。如果观测信号能够正确地反映宏观磁化强度 M,那么,它在零时刻的数值大小将与地层孔隙中的含氢总量成正比。由此经过恰当的标定,即可把零时刻的信号强度(FID或回波串)标定为岩层孔隙度[5]。通过刻度处理后,由T2分布可直接得到孔隙度,即

这意味着观测的孔隙度可以被分解成不同弛豫时间区间的孔隙度,即得到孔隙度分布 P1,P2,…,P8(或至 P10)。它们是与 T2i(i=1,…,n)对应的各孔隙系统在观测到的总孔隙系统中所占的比重。由孔隙度分布信息,可以进一步了解储层质量,例如当孔隙度分布集中在比较小的弛豫时间上,即 P1,P2等占优时,说明储层以微孔为主,如果是碎屑岩,则意味着骨架颗粒很细;当孔隙度分布集中在比较大的弛豫时间上,即P7,P8等占优时,说明该储层以大孔为主,对碎屑岩,意味着骨架颗粒很粗。孔隙度模型下原始数据的计算,就是根据预先定义好的孔隙度(P1,P2,…,Pn),及对应的弛豫时间(T2i)区间分布,再叠加一些实际测井中的噪声,模拟生成原始测井数据。

回波原始数据的模拟计算和回波拟合是一个相反的过程。

2.2.2 渗透率模型设计

核磁共振渗透率的测量是基于试验和理论模型及其相互关系的结合。当这些模型中或关系式中其他因素保持常量时,渗透率随联通孔隙度的增加而增加。渗透率的单位达西(D)是一个面积单位。从岩石理论应用的实际考虑,可以认为渗透率与某些几何尺寸的平方成正比。毛管压力曲线和渗透率之间的相关性表明这一尺寸就是孔隙喉道尺寸。核磁共振虽然测量孔隙体的尺寸,但是在几乎所有砂岩地层和部分碳酸盐岩地层中,孔隙尺寸和孔喉尺寸有很强的相关性。渗透率是反映孔隙介质(岩石)允许通过流体的能力。迄今,估计渗透率的方法都是间接的,核磁共振也一样。通过对渗透率与核磁共振特性之间的相关分析,来建立相应的渗透率模型。利用岩石核磁共振的弛豫特性及扩散测试结果,已经建立了多种有关渗透率的经验公式。渗透率与孔隙度及岩石比表面积有关。基本表达式是所谓的K ozeny公式,即：

式中,K是渗透率,mD,φ为孔隙度,%;S/V为岩石的比表面积,(cm2/cm3);Γ为“结构因子”,或“弯曲因子”,无量纲,其量子决定与孔隙的形状以及单位长度内多孔固体中流体流过的路径。利用K ozeny方程,通过岩石核磁共振弛豫时间与岩石孔隙比表面积的相关性,可以建立岩石渗透率的方法。下面主要介绍Coates束缚水-渗透率模型,其模型主要公式：

式中,FFI为自由流体的孔隙度;BVI为束缚水孔隙度;C为系数。Coates模型中,孔隙尺寸参数是通过T2截止值隐含输入的,它确定 FFI和BVI的比值,同时FFI/BVI值得变化描述了S/V的变化。工作在该模式下时,用户需输入总孔隙度(φ),T2截止值,自由流体孔隙度(FFI),束缚水孔隙度(BVI),模拟器便可以在满足要求的情况下随机计算孔隙度与对应的弛豫分布,模拟出仿真数据。

模拟出的数据与孔隙度模型是一样的,不同的是需满足一定的条件。

3 结束语

通过该模拟器,地面系统可以很方便地进行前期控制与采集处理程序的开发和测试。在开发地面系统前期,成功摆脱了传统必须挂接实际井下仪器进行程序开发与试验的束缚,大大提高了地面系统的开发效率。在后期应用阶段当地面系统与井下仪器挂接联调出现故障时,通过该模拟器可以快速查找故障,节省了大量维修时间。目前该方案已经成功应用于核磁共振仪器的测试体系中,并得到了很好的应用效果。

[1] George Coates,肖立志,Manfred Prammer著,孟繁莹 译.核磁共振测井原理与应用[M].北京：石油工业出版社,2007

[2] 邵维志,庄 升,丁娱娇.一种新型核磁共振测井仪——MREx[J].石油仪器,2004,18(2)

[3] 肖立志.核磁共振测井资料解释与应用导论[M].北京：石油工业出版社,2001

[4] 肖立志.核磁共振成像测井与岩石核磁共振及其应用[M].北京：科学出版社,1998

[5] 肖立志,谢然红.核磁共振测井仪器的最新进展与未来发展方向[J].测井技术,2003,27(4)

Design of simulator based on NMR downhole tools.

Zhang Jiawei,Song Gongpu,Sai Fang,Xue Zhibo.

As an advanced technology,NMR is widely used in actual well logging.In this paper,a design of simulator based on NMR downhole tools is introduced.Through this simulator,we can generate data which is similar with actual data logged from well logging.And through a certain communication protocol,we can transmit simulated data to surface system for analysis and process.Through this simulator,surface system is very convenient to debug software and test processing arithmetic.This simulator improves the efficiency of software-writing and previous test obviously,so it is valuable in application.

NMR(Nuclear magnetic resonance);echoes;T2 spectrum inversion;porosity;permeability;diffusion coefficient

TN912

B

1004-9134(2011)06-0007-04

中海油企业发展基金(H04010701W070552)

张嘉伟,男,1982年生,电子技术工程师,现在中海油服油技事业部从事核磁共振测井仪器研制工作。邮编：101149

2011-05-04

高红霞)

PI,2011,25(6)：7～10

·开发设计·