王鹏飞，田 豆，王 轲

（1.西安石油大学，陕西 西安 710065；2.陕西省油气井测控技术重点实验室，陕西 西安 710065）

在电阻率测井中，传统的一维建模缺乏对三维地层建模的能力，而三维建模数值计算方法形式复杂，计算时间较长，对计算机性能要求较高。为了快速、准确地计算复杂地层的电导率，建立一种模拟三维电阻率测井的二维地层电阻率测井模型。

1 构建有限差分模型

假设电磁波沿y轴方向均匀传播，则一般的导波结构可以表示为

为了抵消时间变化因子 ejωt，式 （1） 与 （2）中沿y方向的导数使用来代替，本文以下关于y方向的导数都用-jk来表示。

已知麦克斯韦方程

对电磁场分量进行归一化则可以得到

通过沿传播方向方向进行傅里叶变换，将场矢量和源矢量从时域转换为频域，将三维问题转化为二维问题，有限差分法使用Ex，Ev，Ez和Ex，Ev，Ez分量对2D空间中的磁（或电） 偶极子的场进行建模。

对函数进行傅里叶变换，得

其中f（y） 表示电场强度E、磁场强度H、电流密度J、磁化强度M。

在交错的Yee网格结构中，各电场分量位于网格单元棱边，各磁场分量位于网格单元面中心。

在应用有限差分法进行建模计算时，其仿真时间不但取决于2D建模计算方法的复杂度，而且与逆傅里叶变换的正交点数相关，点数越多，计算时间越长。我们通过截断积分来提高积分效率，设置积分区间为 [-kym，kym]，kym为最大场振幅的0.1%。

求解（10）式可得时域电磁场分量。

2 数值模拟计算

为了对交错网格有限差分法及计算效果进行验证，我们在MATLAB平台上构建数值计算地层模型，分别构建了地层断层模型以及地层裂缝模型。首先对本文提出的算法进行验证。

图1 2.5D与1D数值解对比图

图1为均质三层地层模型与1D模型视电导率计算结果对比图。图1中两曲线基本重合，有限差分法数值计算结果正确。

图2为不同的傅里叶逆变换正交点数数值解与解析解视电导率对比图，通过综合对比计算精度与计算时间，本文中所有模型的Ky数值定为81。

设计一种地层断层模型进行计算，设计的模型二维剖面如图3所示。其中目的层电阻率R2=R5=10Ω·m，围岩电阻率 R1=R3=R4=R6=1Ω·m，仪器工作频率为f=20kHz，收发间距为L=1m，采样间隔为0.4m，计算域为40m。目的层厚H=4m，断距h=1m，a为断层倾角。仪器位于上目的层中沿x方向水平移动，穿过断层界面到达下目的层，模型中仪器初始位置位于上目的层z方向中间位置；仪器到达下目的层后，距上边界1m，距下边界3m。

图2 数值解与解析解对比图

图3 断层剖面示意图

图4为不同断层倾角a的断层视电导率曲线。图4中曲线表明，当断层倾角a越小，仪器穿过断层面时的视电导率越小；随着a的不断增大，断层面处的视电导率随之增大，“犄角”不断变小，直至消失。曲线中“犄角”现象的出现是因为断层倾角小，断层面处下目的层对上目的层产生影响，导致视电导率下降；当倾角变大时，这种影响逐渐减弱，当倾角为80°时，下目的层的影响完全消失。当仪器进入下目的层时，视电导率增大，出现这种结果的原因是因为仪器在进入下目的层后，仪器与目的层和围岩边界的距离变小，受到围岩影响的效应增强，导致下目的层的视电导率增大。

图4 断层视电导率

利用三层地层水平井数值计算结果。其中模型计算域为40m，目的层厚为4m，目的层电导率为σt=0.1s/m，围岩电导率为σs=1s/m，仪器工作频率为f=20kHz，收发间距为L=1m，采样间隔为0.4m。图4中曲线在目的层与围岩交界面处因角度以及采样间隔过大的原因出现“犄角”效应。此结果表明应用有限差分法进行水平井探测，层边界明显，计算效果好。

3 结语

在对复杂地层进行电磁测井时，一维模型难以构建复杂三维地层，三维模型又存在结构复杂，计算矩阵规模巨大，计算消耗时间长，对计算机性能要求高。本文在MATLAB软件平台上应用提出交错网格有限差分法能够很好的以二维的形式模拟三维地层建模，对地层中的断层有着很好的探测效果，能够真实、准确的反映地层电导率。计算时间与三维计算软件COMSOL相比大幅减少，表明交错网格有限差分法形式简单、计算效率高。