论文部分内容阅读
本文主要对大地电磁(Magnetotelluric)数值模拟中现有正反演方法进行了改进,在有限元法(FEM)的基础上引入了迭代有限元,通过将地下计算空间划分为目标算区(内区)及迭代算区(外区)两部分,在外区中采用迭代消去的处理技术,将外区中的能量叠加到内区的边界上,通过此种处理方式既减少了有限元最终求解方程组的阶数,节省了计算内存,提高了程序运行效率,又保证了大地电磁正演计算时将边界条件视为无穷远的条件,提高了正演计算的精度。该算法相当于将有限元与边界元相结合,在针对相同的内区时,外区的迭代矩阵只需计算一次,这样就给后续的反演工作提供便利,提高了大地电磁的反演速度。通过研究迭代有限元算法的基本原理与特点,建立了基于迭代有限元算法的正演模型,对模型进行理论试算,检验了算法的正确性,并对算法的几个重要参数进行了分析,研究了不同参数值的选取对迭代有限元算法精度的影响。同时将智能优化领域的全局优化算法-差分进化算法引入到大地电磁反演中,充分利用该算法在处理优化问题时进行全空间并行搜索,算法执行效率高,鲁棒性强的优点,避开了以往在进行大地电磁反演时需要计算偏导数矩阵,易忽略细节信息等缺点。通过对标准差分进化算法的研究分析,针对其在搜索过程中易陷入局部最优及收敛慢等缺点,从操作算子的改进和算法策略集成等层面对其进行了改进,采用双算子策略动态集成的方式提升算法的性能,并将改进后的算法利用标准函数进行了测试,此外,还将改进差分进化算法(TDE)与标准差分进化算法,模拟退火及遗传算法的寻优性能进行了对比,结果表明,TDE算法具有寻优速度快,精度高,鲁棒性好,易逃离局部最优的优点。在上述基础上,结合当前大地电磁反演的正则化理论,通过选择合理的反演目标函数,建立了性能良好的智能优化非线性反演模型。基于Matlab平台编写了迭代有限元正演及改进差分进化算法的大地电磁正反演程序,并采取室内数值模拟与野外实测数据验证的方式对模型进行了评价,结果表明,所建立的智能优化模型具较好的反演效果,具备野外推广的能力。