[发明专利]一种地震散射波成像方法

权利要求书

1.一种地震散射波成像方法，包括下述步骤：

1)读入地震数据；

2)计算每一输入道每个采样点映射到散射点道集中所在的等效偏移距位置，抽取散射点道集；所述的共散射点道集形成采用下述方法处理：

ti=2xhvmig(x2+h2-hei2)12]]>

其中h：半偏移距；x：共散射点道集CSP到共中心点CMP的距离；he：等效偏移距；v：散射点速度；t：波从震源到散射点再到接收点的旅行时；i＝1，2，3…，共散射点道集中地震道号；

按照上述公式将输入道数据映射到共散射点道集上；

3)对步骤(2)的结果在共散射点道集上进行叠加速度谱分析；

4)对步骤(3)速度分析数据进行横向和纵向插值，为变速共散射点道集形成提供速度参数；

5)对原始输入数据，根据步骤(4)得到的速度信息，重新抽取共散射点道集，然后进行步骤(3)，直到获得合理的偏移速度信息；

6)对步骤(5)获得的共散射点道集与速度分析结果，进行能量归位，得到偏移后的时间剖面，即得成像；

其中对于散射波，其成像方法依据下述：

Pout(x,τ,t=0)=Δhe2πΣhe[cosθvrmsrρ(t)*Pin(x,he,t-rvrms)]]]>

其中：P_out(x，τ，t＝0)为成像结果，为CSP道集的输入波场，ρ(t)表示对输入波场的时间导数，符号*为褶积算子，x为散射点道集的位置，h_e为散射点道集中并置等效震源和等效接收点到散射点在地表投影的距离，Δh_e为等效偏移距间隔，∑表示在最佳偏移孔径范围内求和；对于纵波，其成像方法依据下述：

Pout(x,τ,t=0)=Δhe2πΣhe[cosθvsemrρ(t)*Pin(x,he,t-rvsem)].]]>