摘要
为了进一步提高广域电磁法正演求解精度,采用基于电流密度法向连续为依据的后验误差估计策略的自适应有限元技术,实现了三维广域电磁法正演模拟。首先,从麦克斯韦方程出发,推导了可控场源满足的电场方程,采用非结构四面体网格的矢量有限元技术对电场方程进行离散,从而形成大型稀疏复线性方程组;其次,设计电流密度法向连续的后验误差估计策略,并将其融入大型稀疏复线性有限元方程组,实现了高精度三维广域电磁法正演计算;最后,采用一种优化Brent-Dekker求根算法,提高广域视电阻率计算的稳定性,并利用层状介质模型验证了自适应算法的正确性,同时设计了典型的地形和球状地电模型,分析广域电磁响应特征。研究结果表明,自适应有限元法能够避免人为设计网格带来的低效率问题,通过对重点区域进行局部网格加密提升正演模拟精度;同时,在同等收发距情况下,广域电磁法相比可控源音频大地电磁(CSAMT)具有更大勘探深度与勘探范围,受非平面波的影响更小。
- 单位