论文部分内容阅读
随着当今地球物理方法的不断发展,重力勘探作为矿产勘探中的重要方法有其不可替代的地位和价值。而物性反演则是重力资料解释中的重点和难点。反演技术的关键在于充分利用约束条件和先验地质信息尽量避免多解性。目前,遗传算法在地球物理反演中的应用主要见于地震波场和电磁场反演方面,而在重力位场反问题中的应用较少。注意到与传统的最优化方法相比较,具有全局收敛、搜索效率高、不计算目标函数偏导数以及隐含并行计算等优点,这些将会在重力反演中发挥作用。所以本文重点研究了重力异常三维约束反演,在反演技术中引入了遗传算法和多种约束手段,对比了不同模型下约束反演的理论效果,并通过处理实际资料验证了约束反演的良好效果。首先在研究了重力场基本理论的基础上,重点推导了立方体场源的重力异常计算方法,得到了立方体重力异常正演公式,编写C语言程序实现正演计算。并且深入研究了如何优化多立方体正演计算速度,为反演打下基础。其次研究了遗传算法和Levenberg-Marquardt最小二乘法。编写了两种算法的MATLAB程序,构建了目标函数并对地下空间进行了合理的网格划分。并将深度加权引入到目标函数中,克服了核函数的趋肤效应。最后加入了正定约束和参数上下限约束,使得反演结果更加符合地球物理场定律。再次分别设计了单一立方体模型、双立方体模型等多个理论模型,运用了遗传算法对各模型进行反演,并加入光滑约束、深度加权、正定约束、参数约束。对比分析了有无约束时的反演效果,肯定了约束效果。最后将该反演方法用于实际重力探勘资料的处理,验证该方法的有效性。