随钻超深前视电阻率测井二次场耦合势三维有限体积法高效数值仿真  

3D efficient finite-volume simulation of the LWD ultra-deep look aheadresistivity measurement using the secondary field coupled potentials

作　　者：王亚洲 汪宏年[1] 莫修文[2] 康庄庄[1] 殷长春[2] WANG YaZhou;WANG HongNian;MO XiuWen;KANG ZhuangZhuang;YIN ChangChun(International Center for Computational Method and Software,College of Physics,Jilin University,Changchun 130012,China;College of Geo-exploration Science and Technology,Jilin University,Changchun 130026,China)

机构地区：[1]吉林大学物理学院计算方法与软件国际中心,长春130012 [2]吉林大学地球探测科学与技术学院,长春130026

出　　处：《地球物理学报》2022年第12期4981-4996,共16页Chinese Journal of Geophysics

基　　金：中国科学院先导专项A(XDA14020102);国家自然科学基金项目(42174150)联合资助

摘　　要：随钻超深前视电阻率测井是一种新型井中电磁测量技术,能够实现复杂储层探测和开发中井眼轨迹优化与地质导向.本文将基于二次场耦合势三维有限体积法建立一套非均质各向异性地层中随钻超深前视多分量电阻率测井高效数值模拟算法.首先,利用Yee氏交错网格、非均质网格上的等效电导率标准化技术以及三维有限体积法实现二次场耦合势Helmholtz方程离散,并采用传输线法建立背景电磁场并矢Green函数库,同时结合二维Newton插值技术确定背景电磁场,提高异常体内散射电流场计算效率.在此基础上,针对仪器移动测量方式的特点,通过并行直接稀疏求解器PARDISO求解离散方程,结合并行计算同时确定不同位置发射源产生的电磁场.最后,数值结果验证了算法的有效性,显示了二次场方法具有更高的计算精度;建库插值方法显著提高了背景场计算效率,拓展了二次场方法的应用价值.The LWD ultra-deep look ahead resistivity measurement is a new type of electromagnetic measurement technology in borehole.It is mainly used to solve the problem of wellbore trajectory optimization and geo-steering in complex reservoir exploration and development.In this paper,an efficient numerical algorithm is proposed to simulate the tool response in a 3D inhomogeneous anisotropic formation,which uses the 3D finite volume method(3D FV)based on the secondary field coupled potentials.The Yee′s staggered grid,equivalent conductivity standardization and 3D FV are applied to discretize the Helmholtz equations of the coupled potentials.Meanwhile,the background EM field tensor Green′s function library is established through transmission line method(TLM),and the scattered electrical currents distributed in abnormal bodies are calculated by the 2D Newton interpolation technique to enhance efficiency.Furthermore,the Parallel Direct Sparse Solver(PARDISO)is applied to invert the discrete matrix,and the EM fields generated by the transmitters at different locations are simultaneously determined by parallel computing.Finally,the numerical results verify the effectiveness of the algorithm and show that the secondary field method has higher computational accuracy;the library interpolation method greatly improves the efficiency of background field calculation and expands the application value of the secondary field method.

关 键 词：二次场 三维有限体积法 二维Newton插值技术 超深前视电阻率测井 

分 类 号：P631[天文地球—地质矿产勘探]