王浩森2,
汪宏年1,,,
杨守文1,
殷长春3
1. 吉林大学物理学院计算方法与软件国际中心, 长春 130012
2. 河北建筑工程学院, 河北张家口 075000
3. 吉林大学地球探测科学与技术学院, 长春 130026
基金项目: 中国科学院先导专项A(XDA14020102)，国家重点研发项目(2017YFC0601805，2018YFF01013304)和自然科学基金项目(41574110)联合资助
详细信息
作者简介: 于蕾, 女, 1996年生, 博士生, 2018年毕业于吉林大学物理学院, 目前在吉林大学物理学院攻读博士学位, 主要研究方向为电磁场数值模拟与反演成像.E-mail: 378621327@qq.com
通讯作者: 汪宏年, 男, 1962年生, 吉林大学物理学院教授, 博士生导师, 主要从事石油测井、可控源电磁与海洋电磁等方面的电磁场数值模拟与反演等方面的研究.E-mail: wanghn@jlu.edu.cn
中图分类号: P631收稿日期:2020-08-25
修回日期:2021-03-28
上线日期:2021-09-10
Numerical simulation of the LWD orthogonal azimuth electromagnetic tool with annular antenna recesses by the 3D finite volume method
YU Lei1,,WANG HaoSen2,
WANG HongNian1,,,
YANG ShouWen1,
YIN ChangChun3
1. International Center for Computational Method and Software, College of Physics, Jilin University, Changchun 130012, China
2. Hebei Institute of Architecture and Civil Engineering, Zhangjiakou Hebei 075000, China
3. College of Geo-Exploration Sciences and Technology, Jilin University, Changchun 130026, China
More Information
Corresponding author: WANG HongNian,E-mail:wanghn@jlu.edu.cn
MSC: P631--> Received Date: 25 August 2020
Revised Date: 28 March 2021
Available Online: 10 September 2021
摘要
摘要:针对由多环状天线槽构成的方位随钻电磁波测井仪器真实结构，基于柱坐标系下耦合势Helmholtz方程和三维有限体积法研究建立一套相应的三维电磁响应数值模拟算法.首先，采用柱坐标系下三维不规则Yee氏交错网格对计算区域进行剖分，保证了钻铤表面附近的剖分网格与非光滑刻槽钻铤表面同形；并利用电场延拓边界条件逼近钻铤表面的理想电导体边界条件，保证非光滑刻槽钻铤表面耦合势方程离散精度.在此基础上，应用标准均质化技术确定非均质网格的等效电导率和磁导率，并借助三维有限体积法对柱坐标系中耦合势方程以及环状天线槽内的发射天线进行离散，得到一个大型非对称稀疏代数方程，通过PARDISO直接求解器求解离散方程.最后，通过与柱状均匀地层中模式匹配算法所得的数值结果的对比验证该三维算法的有效性，并进一步通过数值结果研究考察金属钻铤和天线槽对方位随钻电磁响应的影响，以及层状TI地层模型中该仪器的响应特征.
关键词: 方位随钻电磁波测井/
刻槽钻铤/
铁氧体材料/
环状天线槽/
三维有限体积法/
电场延拓边界条件
Abstract:All antennas of the sensor of the LWD orthogonal azimuth electromagnetic tool (LWD-OAEMT) are installed in the annular recesses on the grooved metal drill collar. The recesses are often filled with ferrite or other magnetic materials to seal the coils. In this paper, we apply the Helmholtz equations about coupling potentials and three-dimensional finite volume method (3D FVM) with cylindrical grids to simulate the responses of the LWD tool consisting of multiple antenna recesses on the non-smooth grooved metal mandrel surface. Firstly, irregular Yee's staggered grids in cylindrical-coordinate system are used to guarantee the grid near the mandrel surface is conformal with the non-smooth grooved tool geometry. The electrical field continuation technique is employed to approximate the perfectly electric conductor (PEC) boundary condition in order to efficiently improve the discretization precision of the coupling potential on the grooved drill collar surface. On this basis, the equivalent conductivity and magnetic permeability in a heterogeneous element are determined by the standardization technique, and the 3D FVM is employed to discretize the coupling potential equations in the cylindrical coordinates and the transmitting antennas in the annular recesses. Then, the Parallel Direct Sparse Solver (PARDISO) is used to solve the discrete equations to ensure the accuracy of numerical simulation and increase the stability of solutions. Finally, the 3D FVM is validated with comparison of results by numerical model matching (NMM) algorithm in homogeneous formation. In terms of numerical results, we further examine the influence of metal drill collar and antenna recesses on the LWD tool and its response characteristics in the layered TI formation model with different borehole dipping angles.
Key words:LWD orthogonal azimuth electromagnetic tool/
Grooved metal drill collar/
Ferrite materials/
Annular recesses/
Three-dimensional finite volume method/
Electric continuation boundary condition (ECBC)
PDF全文下载地址:
http://www.geophy.cn/data/article/export-pdf?id=981d7f24-f42c-4948-a213-1719dffbc9b6
