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

摘要

本发明公开了一种地震散射波成像方法，克服了传统的反射波成像方法适用范围窄、无法处理实际面临的多种地质构造和地形成像的缺陷，可以对复杂断块、岩脉、透镜状矿体等地质体成像，增加成像效果；并且所得的散射道集更平滑、能量分布合理、信噪比高、容易获得偏移速度、消除了转换波的阶梯状噪音。进一步的，在上述基础上本发明公开了一种地震散射波成像系统，用于地震散射波成像，具有适用范围广、成像效果好、散射道集平滑、信噪比高的优点，并且可以针对转换波共散射点散射道集消除其阶梯状噪音。

主权项

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)].]]>