基于保角变换三谱域混合相位子波相位谱估计.pdfVIP

062 SPG/SEG 深圳 2011 国际地球物理会议论文集 基于保角变换的三谱域混合 相位子波相位谱估计 ① ① ② 于永才 王尚旭 戚鹏飞 ( ①中国石油大学 ( 北京 ) C N P C 物探重点实验室； ②中国地质大学 ( 北京 ) 地球物理与信息技术学院 ) 摘要 高阶统计量是近年来研究地震子波估计的比较热门的理论，由于保留了相位信 息，因此可以用来估计混合相位子波。基于三谱的混合相位地震子波估计往往因为三谱 相位卷绕的原因，导致子波相位谱求解产生偏差，进而影响了混合相位子波估计的效 果。针对这一问题，本文在 M U 算法和 P a n - N i k i a s 算法的基础上，提出了一种基于保 角变换的三谱域子波相位谱求解方法。通过缩小 F o u r i e r 相位谱的取值范围，有效避免 了三谱相位发生卷绕的情况，从而消除了原相位谱估计中三谱相位卷绕的影响。 关键词 混合相位 子波 三谱 保角变换 谱估计中，从而避免了三谱相位卷绕对混合相 1 前言 位子波相位谱估计的影响。 高阶谱域混合相位子波相位谱估计主要分 2 三谱域地震子波相位谱估计原 为两类，一类是在双谱域实现的，另一类是在 理 三谱域实现的。 在 双 谱 域 中 ， Brillinger[1] ， Lii 等 [2], 根据 Robinson 褶积模型，有 Matsuoka等[3] [4] , Tekalp等 分别提出了基于双谱的 s (t) w(t) r(t) n(t) (1) 信号相位谱估计算法。然而这些算法，除了Lii 其中：s(t)、w(t)、r(t)和n(t)分别表示地震记录、 等的算法以外，都会受到双谱相位卷绕的影 地震子波、反射系数序列和噪声，“* ”为褶 响。为此，Marron 等[5]通过引入相位解卷矩 积算子。由累积量的Bartlett-Brillinger-Rosenblatt 阵，提出了一种相位谱估计最小二乘算法； 公式[12]，可以得到式(1)的四阶累积量表达式 Kang 等[6]提出了三种非相位解卷递推算法； c ( , , ) c ( , , ) m ( , , ) 4s 1 2 3 4r 1 2 3 4w 1 2 3 (2) Zhang[7, 8]应用Lp范数最小化解卷方法来解决双 c ( , , ) 4n 1 2 3 谱域中的相位卷绕问题；Yu等[9]提出了一种基 其中：c (τ ,τ ,τ)、c (τ ,τ ,τ )和 c (τ ,τ ,τ )分别 4s 1 2 3 4r 1 2 3 4n 1 2 3 于保角变换的双谱域最小二乘子波相位谱估计 表示地震记录、反射系数序列和噪声的四阶累 方法。 积量，m (τ