地震波全波及散射非线性全局最优化反演方法研究.pdfVIP

地震波全波及散射非线性全局最优化反演方法研究 李信富李小凡张美根 (中国科学院地质与地球物理研究所北京100029) 地震物理反演理论与方法研究是地球物理学中—个非常重要的研究领域，其进展对人们正确认识和了解地球深部的细 节及物性具有重大意义，对地球动力学的研究也具有不可忽视的推动作用。地球物理反演理论中线性反演方法已经成熟， 而且在理论上也是最完整、应用最广泛的反演方法。然而，大多数的地球物理问题都是非线性的，虽然已有许多经典的非 线性反演方法，但在精度、收敛性和收敛速度方面已不能满足我们的要求，因此必须发展新的非线性反演方法。 地震波散射层析成像是—种高分辨地震成像技术，可提供地球内部各种尺度非均匀体的细致的三维图像。地震波散射 成像将可提供有关岩石层、上下地幔及核幔边界的非均匀性和详细构造的宝贵的，有时甚至是唯一的信息。这些信息对地 球内部研究通常是最为关键的。尽管理论上散射层析成像方法的分辨能力远高于常规地震方法(这种差别甚至可达数量 级)；然而，实现地震波总场与散射场的精确分离绝非易事。因此，必须发展既可处理散射信息，又不必进行波场分离的 方法(如全波方法等)来解决这一难题。 一般而言，与常规地震反演方法相比，全波反演可最大限度地利用所有类型波动的运动学和动力学信息，全波场反演 无需进行波场的识别与分离，5I而极大地减轻处理中的困难，克服分离各类波场的困难，减小引入人为误差的可能性，可 以提高成像的分辨率。为此，我们采用了全波反演。—个好的反演方法应该具有的特点是：反演方法不依赖于初始模型， 即无论初始模型精度如何，只要观测数据精确可靠，其反演结果都将趋近于真实结构。另外，地震波全波及散射反演所涉 及的问题非常复杂，其高维数、多参数及高度非线性等问题，现有的反演方法已经无能为力。因此，我们必须发展适用于 上述问题的非线性反演方法。 本文拟以数学中的共轭方向理论、不动点理论、混沌控制及分形理论为基础，构造一种全新的地震波散射非线性直 接反演方法。由于我们所研究问题的高度复杂性，目标函数的光滑性和连续性得不到保障，甚至有时连日标函数的解析表 达式也难以得到，使得我们必须发展适用于散射及全波反演的非线性直接反演方法。Powell方法是一种有效的局部直接搜 索法，它具有较快的收敛速度。混沌与分形本身就是—种非线性理论，利用混沌变量的随机性、遍历性和规律性特点可以 进行优化搜索，目混沌优化方法容易跳出局部最优点。但是，混沌方法搜索速度幔，效率低，耗时长，耗费资源大，不能 够满足我们的要求。为此，本文利用Powell方法增强解附近的收敛速度和搜索能力，同时引进混沌搜索，以克服Powell 方法所无法避免的局部极值俘获问题，从而使本方法具有高速目全局收敛的性能。我们在谱元法正演模拟的基础匕，利用 谱元法正演计算和分形插值方法合成全波理论地震图，通过细合Powell与混沌优化方法最小化—个拟合差函数得到模型 参数。 I 这—进展将为研究地球物理反演理论及方法研究领域中长期困扰人们的一系列非唯—性问题开辟一条新的途径。由 于本研究所涉及的讨算方法在处理非均匀介质问题时需要大量计算资源，因此与本研究有关的数值计算，都在浪潮天梭 10000大型并行计算集群服务器系统上进行。 参考文献 l efficientmethodfor theminimumofafunctionofseveralvariableswithout l lPowell，MJ．D．An finding calculating derivatives·ComputerJournal,1964(7)：155-162· 12l连长云等．混沌理论及其在地学中的应用前景．吉林地质．第14卷．第2期．1995(6)：31～36．