用传播矩阵法研究层状交错地层中的多分量感应测井

为了实现交错沉积等复杂环境中的电磁场数值模拟,本文在常规横向同性模型的基础上引入了电导率主轴坐标系相对地层坐标系的层理方位角和倾角,建立了交错地层模型.并利用传播矩阵法建立了一维层状交错地层模型中的多分量感应测井仪器响应的正演模拟算法.首先将频率-波数域中的电磁场分解为上行和下行模式波,给出了任意朝向的磁偶极子在无限大地层中模式波的解析解.进一步通过引入地层界面上的透射、局部反射以及广义反射系数矩阵,推导了一维层状地层中的模式波表达式.在此基础上,利用二维Gauss-Legendre积分实现了Fourier逆变换,得到了可用于多分量感应测井模拟的频率-空间域磁场并矢格林函数.最后,通过多个数值模拟结果考察了井眼倾角、层理方位角和倾角变化对多分量感应测井响应的影响.     

用积分方程法模拟各向异性地层中三维电性异常体的电磁响应

本文研究并建立了一种模拟各向异性地层中三维电性异常体电磁响应的积分方程算法.首先讨论了并矢Green函数及其相关积分的计算,将水平层状各向异性地层中的电场并矢Green函数分解成含有奇异项的直达波与非奇异的来自各个层界面的反射和透射波两个部分,再应用等效体积单元和表面积分技术对积分方程的奇异核进行离散化处理以便提高离散方程的精度.然后为了节省计算机内存以及计算时间,引入基于Krylov子空间的迭代算法求解积分方程的离散化矩阵方程.最后通过与现有文献中的结果作对比从而检验了所述算法的有效性,并结合具体算例考察分析了地层的各向异性对三维电性异常体电磁响应的影响特征和规律.     

水平层状非均质TI地层中仪器偏心情况下三维感应测井响应高效数值模拟与响应特征分析

本文针对中国石油集团测井有限公司(CPL)三维阵列感应测井仪器(Three Dimensional Induction Tool,TDIT)的真实结构与参数,利用2.5维数值模式匹配技术研究建立水平层状非均质TI地层中仪器偏心情况下三维感应测井响应的高效数值模拟算法.首先,利用2.5维数值模式匹配技术给出水平层状非均质TI模型中偏心磁流源并矢Green函数的半解析解.在此基础上,结合三维阵列感应测井仪器参数以及测量过程中出现的仪器偏心与旋转等实际情况,通过叠加原理推导出该仪器响应的有效计算公式.最后,通过数值模拟结果分析考察仪器偏心、仪器自旋角等参数对测井响应的影响以及变化特征,为多分量感应资料处理与解释提供理论基础.     

环状刻槽钻铤中随钻方位电磁波响应的混合算法

为了提高随钻电磁波测井仪器的抗震能力,仪器传感器中所有发射与接收线圈均安装在柱状金属钻铤上的环状刻槽中,同时槽中充填高磁导率的铁氧体,以提高发射和接收效率.本文针对多环状刻槽钻铤中发射与接收线圈分布,研究建立一套随钻电磁响应的混合模拟算法.首先,根据环状刻槽几何尺寸以及轴向边界,将整个钻铤与外部的井眼和地层划分成多个不同电导率和磁导率分布的轴对称非均质各向同性层状模型,并应用混合法建立层状非均质地层中电磁场的半解析解与各个线圈上感应电动势的计算方法.其次,利用摄动原理给出随钻电磁测井中振幅比与相位差的空间灵敏度计算公式.然后,通过数值结果考察铁氧体磁导率和地层电导率变化对发射线圈的发射功率和接收线圈上感应电动势的影响,并确定各个接收线圈上感应电动势的动态变化范围,为仪器设计过程中选择合适的铁氧体提供理论依据.最后,通过空间灵敏度分布特征考察振幅比与相位差在探测特征上的差异.     

基于耦合标势与矢势的有限体积法模拟非均匀各向异性地层中多分量感应测井三维响应

为快速有效地研究、考察各向异性地层条件下多分量感应测井的响应特征,本文利用电场标势与矢势的有限体积法研制出三维频率域电磁场响应的数值模拟算法,克服由低频发射或高阻地层产生的低感应数问题,有效提高了三维电磁数值模拟算法的应用范围和计算效率.首先利用电场的标势与矢势将Maxwell方程转化为满足库仑规范条件的耦合势Helmholtz方程,以Yee氏交错非均匀网格中不同位置上的节点为中心建立四种控制体积单元,通过对控制体积单元中电磁场与电导率的积分平均实现耦合势方程和磁偶极子旋度的离散,并得到一个对角占优的大型稀疏复线性代数方程组,然后,通过不完全LU分解预处理和稳定双共轭梯度法快速求解离散方程.数值结果证明了该算法的有效性,并进一步考查了仪器偏心、倾斜井、垂直裂缝等复杂条件下多分量感应的响应特征.     

倾斜各向异性地层中多分量电磁波测井响应三维时域有限差分（FDTD）算法

利用三维时域有限差分技术（FDTD）研究建立倾斜各向异性地层中多分量电磁测井响应的数值模拟算法.首先,利用Yee氏交错差分格式和线性内插公式对各向异性介质中时域Maxwell方程进行离散化,得到各向异性介质中各个交错节点上电磁场的时域递推公式,同时将面积加权平均与旋转矩阵技术相结合确定各个网格节点的等效电导率张量,此外,单轴各向异性吸收边界条件（UPML）被用于降低求解区域外边界反射.在此基础上,利用单频正弦磁偶源计算时域电磁场空间分布以及特殊的双方程算法（2E2U）从数值结果中快速提取电磁场幅度和相位,确定多分量电磁测井仪器的三维响应.最后,通过数值结果证明该算法的有效性,并系统考察不同层厚、地层倾角以及侵入带半径等变化对多分量电磁测井响应的影响.     

水平层状介质中双侧向测井资料的迭代Tikhonov正则化反演

根据非线性反演理论与Morozov偏差原理研究建立从双侧向测井(DLL)资料中同时重构地层原状电阻率、侵入带电阻率、侵入半径、层界面位置以及井眼泥浆电阻率的迭代正则化算法.首先利用Tikhonov正则化反演理论将双侧向测井资料的反演问题转化为含有稳定泛函的非线性目标函数的极小化问题,并利用Gauss-Newton算法确定极小化解.为得到稳定的反演结果并有效实现测井资料的最佳拟合,在迭代过程中将Morozov偏差原理和Cholesky分解技术相结合,建立了一套后验选择正则化因子的方法.最后通过理论模型和大庆油田实际测井资料的处理结果,验证了该算法能够取得更为满意的反演效果.     

水平层状非均质横向同性地层中阵列多分量感应测井响应的快速计算

应用模式匹配算法研究建立水平层状非均质横向同性地层中多分量感应测井响应的快速算法.首先,利用Fourier级数展开法将多分量感应响应的数值模拟转化为三个轴对称问题,并利用电阻率径向导数的奇异表达式,引入两个附加奇异微分算子,用于描述柱状分界面上的积累面电荷对共面线圈系电磁响应的影响.然后通过模式匹配算法求解轴对称问题,得到水平层状非均质横向同性地层中多分量感应磁场的半解析解以及测井响应计算方法,最后通过数值模拟结果对该算法进行检验并进一步考察阵列多分量感应仪器的响应特征.     

利用2.5维数值模式匹配算法高效高精度建立垂直井眼中多分量阵列感应井眼校正库

本文将应用2.5维数值模式匹配算法研究建立柱状横向同性地层中偏心条件下磁流源并矢Green函数的高效算法,并利用该算法高效计算多分量阵列感应井眼校正库.首先,通过Fourier展开技术将偏心条件下磁流源并矢Green函数的数值模拟转化为一系列的轴对称问题;然后利用模式匹配算法求解轴对称问题,得到柱状介质中偏心条件下磁流源并矢Green函数的半解析解;最后,结合中国石油集团测井有限公司开发的三维感应测井仪器(TDIT)的结构参数,推导出复合线圈系视电导率张量的计算公式.在此基础上,根据井场实际地层资料归纳出理论地层模型相应参数的变化范围,高效高精度建立一套直接应用于工业生产的多分量阵列感应井眼校正库,并利用井眼校正库详细地分析考察各模型参数对TDIT响应的影响.     

垂直井眼中三维阵列感应测井资料自适应井眼校正

三维感应测井受井眼环境影响很大,如何消除这些不利影响是有效获取各向异性地层电导率真值并提高资料处理和解释质量的一项重要工作.本文以一维柱状各向异性地层模型为基础,通过最小平方反演技术从实际测量资料中快速反演井径、偏心距、泥浆电导率、地层水平电导率和各向异性系数等模型参数,设法减小井眼环境对测量结果的影响.首先,利用仪器偏心情况下不同旋转角电导率张量的关系方程,建立旋转角提取与校正方法,获取无旋转角三维感应测井数据.在此基础上,根据事先计算好的井眼校正库并结合多维非线性有限元逼近技术,建立快速计算三维感应测井响应以及Jacobi导数矩阵的插值公式,并利用自适应正则化迭代反演技术不断修改模型参数,实现理论合成资料与输入资料的最佳拟合.最后,通过反演出的模型参数计算三维感应测井资料所有分量的校正量,实现三维感应资料的井眼校正.理论模型和实际资料的处理结果均证实,在低阻井眼泥浆情况下,一维柱状模型中的反演方法仍然可以提取出较可靠的地层水平和垂直电导率,且井眼环境校正对于受井眼环境影响相对较小的测井曲线具有非常好的校正效果.