小波变换多尺度地震走时层析成像方法.pdfVIP

地 球 学 报 第20卷增刊 ACTAGEOfrIF．NTIASINICA Vol20‰ 1 999茫 型!!!生!生壁堕垒坐竺兰!兰!!!型兰!竺i !!!! 小波变换多尺度地震走时层析成像方法+ 裴正林余钦范 (中国地质大学，北京) 摘要复杂模型的地震走时层析成像是一难度大且急需解决问题。由于复杂模型存在多个 局部极小，常规的迭代反演方法不能得到好的解。本文提出了小波变换多尺度地震走时反演方 法，数值试验和实际资料的反演结果表明．小波变换多尺度方法不仅有效地改善了常规线性化 反演方法易陷入局部极小的不足，而且由于小渡多尺度反演方法实质上是一种阶梯式逐次逼近 的反演方法，也是大扰动速度异常体反演成像的有效方法。另外，小渡多尺度方法还能够提供 不同分辨率的反演图像，有益于图像的解释。 关键词地震走时层析成像小波变换 多尺度分析 井间地震层析成像(CT)是非线性反问题。常规的线性化CT反演方法只有当初始模 型在目标函数的全局极小的邻域内(在真实模型附近)时，近似线性化的关系才成立。所以， 只有当引入很好的背景速度场先验信息约束时，线性化CT反演结果才能接近全局最优解， 否则易陷人局部极小。地震走时非线性反演的困难主要是目标函数中存在许多局部极小， 因而影响了迭代反演达到全局极小。这些局部极小点产生的主要原因包括3个方面：①局 部线性化近似；②观测数据中含有噪声；⑧实际介质的非均匀性。Bunks等L1o提出了多重网 格方法，该方法使反演结果稳定地收敛到全局最小，且加速反演计算。小波变换也是一种多 分辨变换，它能够将信号的低频成分映射到由大尺度小渡基张成的粗采样子空间，把高频成 分映射到由小尺度小波基张成的较细采样子空间。因此，小波变换为多重网格法提供了一 种自然框架。本文应用小波变换对观测走时数据进行多尺度分解，这样～来，大尺度的走时 目标函数局部极小少，而且互相分得很开，迭代反演方法可以得到比较好的模型低频成分即 背景速度分布，即大尺度反演结果给H{了一个接近全局极小的初始模型，这样就保证了在小 尺度时解稳定、快速的收敛到全局极小。实验结果表明，小波变换多尺度方法不仅有效地改 善了常规线性化反演方法易陷入局部极小的不足，而且由于小波多尺度反演方法实质上是 一种阶梯式逐次逼近的反演方法，它也是大扰动速度异常体反演成像的有效方法。另外，小 波多尺度方法还能够提供不同分辨率的反演图像，更有益于图像的解释。 1小波多尺度分析 则尺度函数^，。(z)和小波函数分别为以。。(z)分别为 ·“九五”国家科技攻关项甘资助(95—913—060203) 第一作者：裴正林，男，1962年生，高级工程师，从事地球物理研究，邮编：100083 增刊 裴正林等：小波变换多尺度地震走时层析成像方法831 ≠…(_72)=2一j／2≠(2一k—n) (1) 心．。(z)=2一，4驴(2一Jx一”) (2) 其中，J，n∈Z表示为尺度函数和小波函数的伸缩和平移。对于正交小波变换，函数 (z)满足二尺度方程 ≠(o)=2h。4．(2x一”) (3) 庐(z)=殛。≯(2z一”) (4) 其中，{h。}和{岛}分别为尺度函数和小波函数对应的滤波器。对于二进正交小波变 换，它们满足下面的关系 g。=(一1)“^l-。 (5) 对于任意的f(z)∈L2(R)，有如下的小波展开式