基于Python的直流电阻率法三维有限元数值模拟  

Three-Dimensional Finite Element Numerical Simulation of Direct Current Resistivity Method Based on Python

作　　者：石靖 庞满莹 梁卓 崔杨洋 陈汉波[1] 喻佳敏 李耀昕 

机构地区：[1]桂林理工大学地球科学学院,广西 桂林 [2]江西省地质局物化探大队,江西 南昌 [3]中国冶金地质总局第三地质勘查院,山西 太原

出　　处：《地球科学前沿(汉斯)》2024年第12期1604-1615,共12页Advances in Geosciences

基　　金：广西自然科学基金(2024GXNSFBA010257);广西人才基地专项基金(桂科AD23026231);江西省地质局青年科学技术带头人培养计划项目(2024JXDZKJRC10)资助。

摘　　要：利用Python平台强大的科学计算库,拟实现电导率分块均匀的直流电阻率法三维有限元数值模拟。基于Galerkin法给出稳定电流场三维双电源总场电位边值问题的控制方程,采用六面体单元剖分和三线性插值进行有限单元分析,将控制方程化为线性方程组,分别采用压缩稀疏列矩阵(CSC)与压缩稀疏行矩阵(CSR)存储总体系数矩阵,对比科学计算库SciPy及PyPardiso的求解速度与准确度,试算层状介质模型及单一、组合异常体模型,结果表明:H型层状介质的视电阻率测深解析解与有限元数值解最大相对误差低于3.00%,平均相对误差1.54%,本文算法有效;在Python平台求解地球物理正演中大型稀疏线性方程组时,PyPardiso与SciPy精度一致,当方程未知个数达37万时,PyPardiso比SciPy求解速度快约66倍;以CSC与CSR形式存储的大型稀疏线性系统,在同一求解器中求解效率相当。Utilizing the powerful scientific computing libraries on the Python platform, this study aims to implement a three-dimensional finite element numerical simulation for direct current (DC) resistivity with uniform conductivity blocks. By applying the Galerkin method, we derived control equations for the boundary value problem of a steady-state electric field with dual-source potential. Hexahedral meshing and trilinear interpolation were used for finite element analysis, converting the control equations into linear algebraic systems. The global coefficient matrix was stored in Compressed Sparse Column (CSC) and Compressed Sparse Row (CSR) formats to evaluate the speed and accuracy of solvers in SciPy and PyPardiso. Test cases included layered media models and single and combined anomaly models. Results show that, for H-type layered media, the maximum relative error between the apparent resistivity depth profiling analytical solution and the finite element numerical solution was below 3.00%, with an average relative error of 1.54%, validating the algorithm. When solving l

关 键 词：直流电阻率法 有限单元法 数值模拟 压缩存储 PYTHON 

分 类 号：TP3[自动化与计算机技术—计算机科学与技术]