各向异性地层随钻方位电磁波测井响应数值模拟.pdfVIP

中图分类号：TE122+P631 单位代码：10425 学 号： 各向异性地层随钻方位电磁波 测井响应数值模拟 Numerical simulation of azimuthal electromagnetic LWD response in anisotropic formation 学科专业： 地球探测与信息技术 研究方向： 测井理论、方法与技术 作者姓名： 胡云云 指导教师： 范宜仁 教授 二〇一三年六月 Numerical simulation of azimuthal electromagnetic LWD response in anisotropic formation A Thesis Submitted for the Degree of Master Candidate：Hu Yunyun Supervisor ：Prof. Fan Yiren School of Geosciences China University of Petroleum (East China) 摘 要 随着油气勘探复杂程度的不断增大，大斜度井、水平井的广泛应用以及国外随钻方 位测井技术的迅猛发展，国内急需开展各向异性地层随钻方位电磁波测井响应的模拟研 究以及其方位特性的研究，为实时地质导向以及储层评价提供依据；本文分别模拟了一 维地层模型，以及考虑了复杂井眼环境的三维地层模型，分析其响应特征，并就随钻方 位电磁波的方位指向性进行了模拟分析。 首先，针对一维层状各向异性地层模型，将发射线圈等效为磁偶极子源，并将磁偶 极子源分解为水平磁偶极子源与垂直磁偶极子源，分别推导两个方向磁偶极子源产生场 的解析解，进而通过叠加获得任意方向磁偶极子源电磁场的分布；并针对 Bessel 函数 无穷积分被积函数的震荡性，分别选择高斯积分与 Hankel 变换，对 Bessel 函数无穷积 分准确求解，从而完成各向异性地层随钻电磁波测井响应的快速模拟。 针对实际钻井过程中复杂的测井环境，建立合理的地层模型，并考虑仪器以及井眼 形态特征采用柱坐标，推导三维各向异性地层时域电磁响应的差分计算格式；通过坐标 转换方法将倾斜单轴各向异性电阻率转换到柱坐标下；在径向以及纵向上均采用非均匀 网格剖分，对线圈以及井眼精细剖分，在界面处采用“阶梯近似”，以及采用介电常数 放大等策略实现各向异性地层条件下随钻电磁波测井响应的精确模拟。 基于各向异性地层随钻电磁波测井响应数值模拟技术，模拟研究了各向异性地层随 钻电磁波测井响应特征，以及层厚、围岩、泥浆电阻率、井眼、井斜角、电阻率各向异 性系数、工作频率以及线圈距等参数对仪器响应的影响；并通过与传统随钻仪器响应特 征比较，对随钻方位电磁波测井仪器的方位指向性进行了模拟分析。 为提高三维随钻电磁波测井响应的模拟效率，本文采用了无条件稳定的交替隐式时 域差分方法（ADI-FDTD ），克服了传统时域有限差分（FDTD ）中时间步长受空间步长 的约束；并将两种方法的模拟结果进行对比，验证了算法的正确性。 关键词：随钻电磁波，各向异性，方位指向性，时域有限差分，ADI-FDTD i Numerical simulation of azimuthal electromagnetic LWD response in anisotropic formation Hu Yunyun(Ge