论文部分内容阅读
大地电磁法作为一种重要的地球物理勘探手段,现已被广泛应用于研究地壳和上地幔地质构造以及深部矿产勘探。作为连接地球物理观测与解释的桥梁,反演方法研究一直都是学者关注的热点。本文首先将智能优化算法——果蝇优化算法引入到大地电磁数据反演,避免线性化迭代方法需要计算偏导数矩阵、对初始模型依赖等缺点。果蝇优化算法具有原理简单、控制参数少,容易程序实现等优点,通过对标准果蝇优化算法的分析,发现其在处理高维、多峰的目标函数时存在收敛缓慢、易陷入局部极值的情况,为此对其进行了改进,加入了差分进化算法的交叉操作和变异操作,增加了果蝇的种群多样性,提高全局优化能力,同时利用变异尺度因子,将果蝇的固定搜索步长方式改为逐步递减的搜索步长,以达到平衡算法全局优化和局部优化的目的。利用多个测试函数对改进的果蝇优化算法进行测试,并与标准果蝇算法以及差分进化算法结果进行了比较,结果表明改进的果蝇优化算法具有寻优快,优化精度高,不易早熟收敛的优点。在此基础上,结合大地电磁反演理论,利用改进果蝇优化算法对大地电磁一维模型进行反演,并利用不同噪声水平的模型对算法进行了测试,结果表明改进的果蝇优化算法能有效的处理大地电磁数据,反演结果精度高,算法鲁棒性好。本文还对基于贝叶斯理论的统计反演方法进行了研究,对非线性贝叶斯反演的基本原理、目前常见的非线性数值采样方法进行了归纳和总结。贝叶斯反演理论将反演模型参数看成是随机变量,反演的结果是统计意义上的后验概率分布,能直观的对结果进行评价。基于变维反演的思想,利用可逆跳跃马尔科夫链蒙特卡洛方法对一维大地电磁数据进行反演。贝叶斯反演结果基于大量样本,因此采样速度的快慢对算法具有重要影响,为了加快算法收敛速度,利用改进的并行回火技术,将副本之间的相邻交换方式替换为随机交换方式,使得算法能够快速对整个空间进行采样,获得解的大量样本。结果表明变维反演能够有效的对层状介质进行自动分层,有效减少人为因素的干扰,并行回火技术能够加速采样过程收敛。