三维密度分布遗传算法反演.pdfVIP

第２９卷第１期 大地 测 量 与地 球 动 力 学 Ｖｏｌ．２９Ｎｏ．１ ２００９年２月 ＪＯＵＲＮＡＬＯＦＧＥＯＤＥＳＹＡＮＤＧＥＯＤＹＮＡＭＩＣＳ 　 Ｆｅｂ．，２００９ 　　文章编号：１６７１５９４２（２００９）０１００４１０５ 三维密度分布的遗传算法反演 柯小平　王　勇　许厚泽 （中国科学院测量与地球物理研究所，武汉　４３００７７） 摘　要　提出把棱柱体单元的重力正演计算方法与遗传算法相结合，进行三维密度结构反演。介绍了目标函数 的构建及遗传算法过程，通过模型试算证明了该方法的有效性和稳定性，并利用该方法反演了青藏高原东缘三维 密度分布。 关键词　重力；反演；遗传算法；密度；青藏高原 中图分类号：Ｐ３１３．３　　　　文献标识码：Ａ ３ＤＤＥＮＳＩＴＹＩＮＶＥＲＳＩＯＮＷＩＴＨＧＥＮＥＴＩＣＡＬＧＯＲＩＴＨＭ ＫｅＸｉａｏｐｉｎｇ，ＷａｎｇＹｏｎｇａｎｄＸｕＨｏｕｚｅ （ＩｎｓｔｉｔｕｔｅｏｆＧｅｏｄｅｓｙａｎｄＧｅｏｐｈｙｓｉｃｓ，ＣＡＳ，Ｗｕｈａｎ　４３００７７） Ａｂｓｔｒａｃｔ　Ａｎｅｗｉｎｖｅｒｓｉｏｎｍｅｔｈｏｄｆｏｒ３Ｄｄｅｎｓｉｔｙｓｔｒｕｃｔｕｒｅｉｓｐｒｏｐｏｓｅｄ，ｗｈｉｃｈｃｏｍｂｉｎｅｓｔｈｅｍｅｔｈｏｄｏｆｇｒａｖｉｔｙ ｆｏｒｗａｒｄｏｆ３Ｄｒｅｃｔａｎｇｕｌａｒｐｒｉｓｍｗｉｔｈｇｅｎｅｔｉｃａｌｇｏｒｉｔｈｍ，ｔｈｅａｐｐｒｏａｃｈｏｆｃｏｎｓｔｒｕｃｔｉｎｇｔｈｅｏｂｊｅｃｔｉｖｅｆｕｎｃｔｉｏｎａｎｄｔｈｅ ｐｒｏｃｅｄｕｒｅｏｆｇｅｎｅｔｉｃａｌｇｏｒｉｔｈｍａｒｅｉｎｔｒｏｄｕｃｅｄ．Ｔｈｅｖａｌｉｄｉｔｙａｎｄｓｔａｂｉｌｉｔｙｏｆｔｈｉｓｍｅｔｈｏｄａｒｅｐｒｏｖｅｄｂｙｓｙｎｔｈｅｔｉｃｄａ ｔａ．ＷｉｔｈｔｈｉｓｍｅｔｈｏｄｔｈｅｔｈｒｅｅｄｉｍｅｎｓｉｏｎａｌｄｉｓｔｒｉｂｕｔｉｏｎｏｆｄｅｎｓｉｔｙｉｎｅａｓｔｅｒｎＱｉｎｇｈａｉＴｉｂｅｔｐｌａｔｅａｕｉｓｉｎｖｅｒｔｅｄｆｒｏｍ Ｂｏｕｇｕｅｒｇｒａｖｉｔｙａｎｏｍａｌｙ． Ｋｅｙｗｏｒｄｓ：ｇｒａｖｉｔｙ；ｉｎｖｅｒｓｉｏｎ；ｇｅｎｅｔｉｃａｌｇｏｒｉｔｈｍ；ｄｅｎｓｉｔｙ；ＱｉｎｇｈａｉＴｉｂｅｔｐｌａｔｅａｕ 遗传算法（ＧＡ）作为全局最优化算法，其主要优 １　引言 点是不用求目标函数的偏导数，不用解大型矩阵方 利用重力进行三维密度反演最常用的方法是阻 程组，而且易于加入约束条件，方法易于移植。目 ［１］ 尼最小二乘法 ，但是，如果重力观测数据存在观 前，遗传算法在地球物理学中的应用主要集中在地 测噪声，则该方法在反演过程中会遇到方程的病态 ［２，３］ 震方面 ，在重力方面主要是利用遗传算法进行 性问题，并且该方法强烈地依赖于初始模型，如果初 地壳结构反演［４－８］，在密度反演方面的研究较 始模型好，收敛速度就快，反演结果就好；如果初始 ［９－１１］ 少 。本文把棱柱体单元的重力正演计算方法 模型差，就有可能导致发散，即使收敛，也有可能趋 与遗传算法相结合，进行三维密度结构反演，介绍了 于局部极值，通常是利用地震波速转换得到的密度 目标函数的构建及遗传算法过程，通过模型试算及 或其他地球物理资料作为约束。在先验信息缺乏的