首页|应用非一致性网格和高阶有限元的井地电磁三维各向异性正演算法

应用非一致性网格和高阶有限元的井地电磁三维各向异性正演算法

扫码查看
地下介质通常表现为电各向异性,增加了井地电磁法响应的复杂性.因此,针对各向异性介质模型进行正演模拟并总结规律,对于正确解释电磁勘探的观测数据至关重要.目前,对复杂三维地质模型进行精确建模是正演中的难点之一.本文使用一种基于八叉树网格的高精度建模方法,通过将六面体单元划分为八个较小的单元来局部细化网格,能够使用较少数量单元的同时保证建模的精度,从而降低计算成本.此外,为提高数值解的精度,在计算时使用高阶形函数.结合八叉树网格和高阶形函数实现了井地电磁法三维各向异性正演方法.通过对比各向异性层状模型的数值解与解析解,验证了算法的正确性.使用八叉树网格实现复杂各向异性油气模型精确建模,将异常体设置不同的水平和垂直电阻率,计算结果表明垂直电阻率对响应的影响大于水平电阻率.研究成果对实际的野外施工具有一定的指导作用,为进一步进行各向异性介质的三维反演打下基础.
Three-dimensional anisotropic forward modeling algorithm for borehole-to-surface electromagnetic method using non-conforming meshes and high-order finite element
Electrical anisotropy typically exhibited by subsurface media increases the complexity of borehole-to-surface electromagnetic(BSEM)responses.Therefore,forward modeling for anisotropic media models and summarizing their characteristics are crucial for accurately interpreting EM exploration data.Currently,precise modeling for complex three-dimensional geological structures is one of the challenges in forward modeling.We employ a high-precision modeling method based on octree grids.By subdividing hexahedral elements into eight smaller elements to locally refine the grids,this method ensures the accuracy of modeling with fewer cells,thereby reducing computational costs.Additionally,to improve the accuracy of numerical solutions,we use higher-order shape functions during calculations.The three-dimensional anisotropic forward modeling method for borehole-to-surface EM surveys is achieved by a combination of octree grids and higher-order shape func-tions.The comparison of numerical solutions with analytical solutions of anisotropic layered models validates the correctness of the algorithm.We utilize octree grids to achieve precise modeling of complex anisotropic oil and gas models and set different horizontal and vertical resistivities for anomalous bodies.The results indicate that vertical resistivity has a greater influence on the response than horizontal resistivity.These findings provide practical guidance for field applications and lay the groundwork for further three-dimensional inversion of anisotropic media.

anisotropyoctreefinite element methodhigh-order shape function

秦策、付文亮、王绪本、米晓利、杨云见、赵宁

展开 >

河南理工大学物理与电子信息学院,河南焦作 454003

成都理工大学地球物理学院,四川成都 610059

东方地球物理公司综合物化探经理部,河北涿州 072751

各向异性 八叉树 有限单元法 高阶形函数

2024

石油地球物理勘探
东方地球物理勘探有限责任公司

石油地球物理勘探

CSTPCD北大核心
影响因子:1.766
ISSN:1000-7210
年,卷(期):2024.59(6)