[发明专利]基于单台、双台和双事件数据二维Lg波Q值层析成像方法

权利要求书

1.基于单台、双台和双事件数据二维Lg波Q值层析成像方法，其特征在于，包括如下步骤：

(a)地震资料预处理，获得Lg波真振幅和信噪比，设定信噪比的门限值，选择用于LgQ层析成像的Lg谱数据体；

(b)从单台数据中提取双台和双事件数据；

(c)基于振幅数据的Lg波Q值层析成像。

2.根据权利要求1所述的基于单台、双台和双事件数据二维Lg波Q值层析成像方法，其特征在于，在步骤(a)中：数据预处理分为4个步骤：(a1)分离Lg波列与P波前噪声；(a2)计算Lg波和噪声的Fourier振幅谱；(a3)拾取谱振幅值；(a4)去噪处理。

3.根据权利要求2所述的基于单台、双台和双事件数据二维Lg波Q值层析成像方法，其特征在于，在步骤(a1)中：为了识别Lg波列并选择有效的Lg波时窗，使用递归的零相位四阶Butterworth滤波器对地震记录进行滤波处理，其拐角频率设定为1Hz，从而得到增强的Lg震相以便于确定其群速度窗口；或者在宽频带地震图上卷积一个全球范围标准地震台网的短周期仪器响应，并通过地震图可视化来选择Lg波群速度窗；因为全球范围标准地震台网仪器响应的高通属性，可识别的离散的Lg波在全球范围标准地震台网地震图上清晰可见；通过对全部记录的可视化检测，可以确定Lg波列的群速度窗口为3.6-3.0或3.5-2.9km/s；噪声序列在P波初至前与Lg波列等时窗内提取。

4.根据权利要求3所述的基于单台、双台和双事件数据二维Lg波Q值层析成像方法，其特征在于，在步骤(a2)中：确定Lg波列和噪声的群速度窗和/或时间窗，在宽频带地震图上轻松获得有效波列的持续时间；在波列两端各增加时窗长度的10％用于余弦镶边消减；对镶边的信号进行快速Fourier变换，得到Lg波和噪声的振幅谱。

5.根据权利要求4所述的基于单台、双台和双事件数据二维Lg波Q值层析成像方法，其特征在于，在步骤(a3)中：在0.05至10.0Hz的频带范围内选取参考频率点，并使之在对数域中为等间隔，即，Δlog₁₀f_ref＝0.04；对于每个参考频率(f_ref)，与其对应的Lg波振幅谱(A(f_ref))采用均方根的方法来求取，即，其中f_i∈[f_ref-Δf,f_ref+Δf]，这里的Δf可由Δlog₁₀(Δf)＝0.02来确定。

6.根据权利要求5所述的基于单台、双台和双事件数据二维Lg波Q值层析成像方法，其特征在于，在步骤(a4)中：假定Lg时窗内的波列为Lg波与噪声的能量叠加，并认为Lg信号与背景噪声的谱振幅是不相关的，利用下式估计每个频率点上的Lg信号：

Asig2(fref)=Aobs2(fref)-·Anoi2(fref),---(3-1)]]>

其中A表示谱振幅和窗口长度，下标sig、obs和noi分别指真振幅、观测值以及噪声。