[发明专利]高精度反射系数求取方法和系统

权利要求书

1.一种高精度反射系数求取方法，其特征在于，所述方法包括：

获取目标对象的地震记录和衰减子波；

将所述地震记录和所述衰减子波转换到频域中，得到频域中的地震记录和频域中的衰减子波；

根据所述频域中的地震记录和所述频域中的衰减子波，得到频域中的反射系数；

对所述频域中的地震记录、所述频域中的衰减子波和所述频域中的反射系数做伽伯变换，得到伽伯域中的地震记录、伽伯域中的衰减子波和伽伯域中的反射系数；

对获取的所述目标对象的衰减子波进行奇偶分解，得到衰减子波奇分量和衰减子波偶分量；

结合高阶走时计算，对所述衰减子波奇分量和所述衰减子波偶分量中地震走时进行高阶走时计算，得到伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量；

根据所述伽伯域中的反射系数、所述伽伯域中的地震记录、所述伽伯域中高阶走时衰减子波奇分量和所述伽伯域中高阶走时衰减子波偶分量，得到高阶走时反射系数奇分量和高阶走时反射系数偶分量；

根据所述高阶走时反射系数奇分量和所述高阶走时反射系数偶分量，得到高精度反射系数序列；

根据所述高精度反射系数序列进行地震成像，得到所述目标对象的地震剖面；

其中，所述对获取的所述目标对象的衰减子波进行奇偶分解，得到衰减子波奇分量和衰减子波偶分量，包括：

根据下式对获取的所述目标对象的衰减子波进行奇偶分解：

w_a(t)＝2w_o(t)+2w_e(t)

其中，w_o(t)为所述衰减子波奇分量，w_e(t)为所述衰减子波偶分量，w_a(t)为获取的所述目标对象的衰减子波；

所述结合高阶走时计算，对所述衰减子波奇分量和所述衰减子波偶分量中地震走时进行高阶走时计算，得到伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量，包括：

根据下式得到所述伽伯域中高阶走时衰减子波奇分量和所述伽伯域中高阶走时衰减子波偶分量：

G[w_a](t,f)＝G[w_o](t,f)+G[w_e](t,f)

其中，G[w_o](t,f)为所述伽伯域中高阶走时衰减子波奇分量，G[w_e](t,f)为所述伽伯域中高阶走时衰减子波偶分量，G[w_a](t,f)为所述伽伯域中的衰减子波。

2.根据权利要求1所述的高精度反射系数求取方法，其特征在于，所述获取目标对象的地震记录和衰减子波，包括：

根据地震子波、衰减函数和反射系数的褶积结果得到所述地震记录，根据所述地震子波和所述衰减函数的合并结果得到所述衰减子波；

根据下式计算所述地震记录：

x(t)＝w(t)*a(t)*r(t)

其中，x(t)为所述地震记录，w(t)为所述地震子波，a(t)为所述衰减函数，r(t)为所述反射系数；

根据下式计算所述衰减子波：

w_a(t)＝w(t)*a(t)

其中，w_a(t)为所述衰减子波，w(t)为所述地震子波，a(t)为所述衰减函数。

3.根据权利要求2所述的高精度反射系数求取方法，其特征在于，所述将所述地震记录和所述衰减子波转换到频域中，得到频域中的地震记录和频域中的衰减子波，包括：

根据下式计算所述频域中的地震记录：

x(f)＝w(f)a(f)r(f)

其中，x(f)为所述频域中的地震记录，w(f)为频域中的地震子波，a(f)为频域中的衰减函数，r(f)为频域中的反射系数；

根据下式计算所述频域中的衰减子波：

w_a(f)＝w(f)a(f)

其中，w_a(f)为所述频域中的衰减子波，w(f)为频域中的地震子波，a(f)为频域中的衰减函数。