多层自适应交叉近似算法在复杂地层测井响应积分方程数值模拟应用.pdfVIP

第28卷增刊 电波科学学报 V01．28，Sup． 2013年8月 CHINESEJOURNALOFRADIOSCIENCE August，2013 多层自适应交叉近似算法在复杂地层测井响应 积分方程数值模拟中的应用 兰天聂在平 (电子科技大学电子工程学院，四川成都611731) 摘 要 在复杂地层随钻测并(LWD)响应的数值仿真中，利用多层自适应交叉近似 算法(MLACA)对积分方程的数值求解进行压缩和加速．通过这个方法，对应于长 井段众多测量点的测井响应只需在同一套模型和网格中计算．也就是说，在此长井 段中无论测井仪器如何移动，非均匀地层、井眼、金属钻铤等都只需要剖分一次．由 于这个方法是纯代数的，它的应用与格林函数无关．通过对秩亏子矩阵的LU分解， 大大节省了内存，加速了积分方程求解．数值实例表明，在复杂地层LwD测井响应 数值仿真中，该方法的应用将显著提升求解效率． 关键词 积分方程；随钻测井；多层自适应交叉近似；矩量法 法(MLACA)来提高数值求解效率．一般来说， 引 言 MLACA与LU分解类似，通过压缩子矩阵来减少 在复杂地层环境中随钻测井响应的数值仿真广 内存和复杂度．MLACA的应用不依赖于格林函数 泛应用于测并仪器设计和测井数据解释．在测井响 和基函数的选择，所以它能够很好地集成到现有矩 应的数值仿真中，经常使用基于微分方程的有限元 量法或多层快速多极子算法的程序中．在低频应用 方法(FEM)u1和时域有限差分方法(FDTD)[2]，基 中，它仅具有O(Nlog(N))的复杂度．本文还利用 于积分方程的矩量法(MoM)[3]，数值模式匹配算法块对角化预条件来提高迭代求解的收敛性． (NMM)[41等半数值半解析混合方法，以及共轭梯 1测井理论及模型 度快速傅里叶变换(CG—FFHT)[5]．半数值半解析 主要应用于对称的二维结构；当维数提高，计算复杂 随钻测井模型如图1，在充满泥浆的井眼中有 度将迅速提升．FEM和FDTD在传播空间直接离一个金属钻铤，发射线圈发射电磁波透过井眼以及 散麦克斯韦方程组，所以可以用来求解任意的三维 周围的地层，两个接收线圈接收电磁波信号．最后， 模型．但是，由于基于对传播空间的剖分，未知量数 我们通过获取两个接收线圈电压的幅度比和相位差 目极大；而且，只要激励源的位置改变，网格必须要 来决定钻井过程中的地质导向． 重新剖分并重新生成系数矩阵。在实际应用中，这 种处理浪费了大量的时间。同时，FEM和FDTD需 要设置截断边界，导致截断误差．FDTD还将面对 曲面拟合的困难以及网格色散等问题．相比之下， 积分方程的矩量法不存在这些问题．对于一次很长 的井段中众多测量点上的测井响应，可以使用同样 的一套模型和网格来求解，阻抗矩阵也只需要求解 一次，使得计算时间大大节省．为了进一步减少内 存和计算时间，本文还利用多层自适应交叉近似算 图1随钻测井模型 资助项目：国家自然科学基金重点项目(编号61231001) 联系人：兰天E—mail：lantianee@qq．corn 12 电 波科学 学报 第28卷 我们把钻铤当做理想电导体(PEC)，同时，非均得 匀地层作为损耗介质，它的电导率和介电常数分别 I}彤“lI—l|汐“一z“”I}≤elI汐“|I，(8)