The invention provides a high-precision method and system for calculating reflection coefficient, including acquiring seismic records and attenuated wavelets of target objects, converting them into frequency domain to obtain seismic records, attenuated wavelets and reflection coefficients in frequency domain, and Gabor transformation of seismic records, attenuated wavelets and reflection coefficients in frequency domain to obtain Gabor transform. Seismic records, attenuated wavelet and reflection coefficients in the region; parity decomposition of attenuated wavelet to obtain the parity component of attenuated wavelet; high-order travel time calculation of the parity component of attenuated wavelet to obtain the parity component of high-order travel time attenuated wavelet in the Gabor domain; according to the reflection coefficient and seismic records in the Gabor domain. High-order travel-time attenuation wavelet odd and even components can be obtained, and high-order travel-time reflection coefficient odd and even components can be obtained; then high-precision reflection coefficient sequence can be obtained; and seismic imaging can be carried out to obtain the seismic profile of the target object, which can improve the resolution of seismic exploration and seismic profile imaging accuracy, thus better serve seismic exploration.
【技术实现步骤摘要】
高精度反射系数求取方法和系统
本专利技术涉及地震勘探
，尤其是涉及高精度反射系数求取方法和系统。
技术介绍
在地震波传播过程中，地层的粘弹特性会对震源激发的地震波造成能量衰减。对煤炭勘探而言，比较普遍的一种现象是，上覆强反射地层阻碍地震波向下传播，影响目标煤层的勘探和开发。通常的地震波处理方法是基于弹性框架研究的，因此不能很好解决此类地质问题。实际地下介质并不是理想的弹性体，具有粘滞性，对地震波有很大的吸收衰减作用。地震成像的常用方法一般是通过求解地下介质的反射系数序列来进行成像，在计算地下介质反射系数序列的过程中必然要考虑来自地层对地震波的吸收衰减作用。目前，国内外通过计算衰减函数与地震子波获得反射系数的方法分为两大类，一种是反Q滤波，主要是沿时间方向做地震波的能量补偿，但不适用于沿波场偏移方向进行外推；另一种是反褶积方法，通过反褶积方法求解反射系数序列的过程主要是转换在频率域中进行的，往往在计算的过程中缺少了某些时间域中的信息，不利于反射系数序列的精确求解。因此，需要一种方法结合反射系数在时间域与频率域中的变化特点来计算反射系数序列，伽伯反褶积便应运而生...
【技术保护点】
1.一种高精度反射系数求取方法，其特征在于，所述方法包括：获取目标对象的地震记录和衰减子波；将所述地震记录和所述衰减子波转换到频域中，得到频域中的地震记录和频域中的衰减子波；根据所述频域中的地震记录和所述频域中的衰减子波，得到频域中的反射系数；对所述频域中的地震记录、所述频域中的衰减子波和所述频域中的反射系数做伽伯变换，得到伽伯域中的地震记录、伽伯域中的衰减子波和伽伯域中的反射系数；对所述衰减子波进行奇偶分解，得到衰减子波奇分量和衰减子波偶分量；结合高阶走时计算，对所述衰减子波奇分量和所述衰减子波偶分量中地震走时进行高阶走时计算，得到伽伯域中高阶走时衰减子波奇分量和伽伯域中...
【技术特征摘要】
1.一种高精度反射系数求取方法，其特征在于，所述方法包括：获取目标对象的地震记录和衰减子波；将所述地震记录和所述衰减子波转换到频域中，得到频域中的地震记录和频域中的衰减子波；根据所述频域中的地震记录和所述频域中的衰减子波，得到频域中的反射系数；对所述频域中的地震记录、所述频域中的衰减子波和所述频域中的反射系数做伽伯变换，得到伽伯域中的地震记录、伽伯域中的衰减子波和伽伯域中的反射系数；对所述衰减子波进行奇偶分解，得到衰减子波奇分量和衰减子波偶分量；结合高阶走时计算，对所述衰减子波奇分量和所述衰减子波偶分量中地震走时进行高阶走时计算，得到伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量；根据所述伽伯域中的反射系数、所述伽伯域中的地震记录、所述伽伯域中高阶走时衰减子波奇分量和所述伽伯域中高阶走时衰减子波偶分量，得到高阶走时反射系数奇分量和高阶走时反射系数偶分量；根据所述高阶走时反射系数奇分量和所述高阶走时反射系数偶分量，得到高精度反射系数序列；根据所述高精度反射系数序列进行地震成像，得到所述目标对象的地震剖面。2.根据权利要求1所述的高精度反射系数求取方法，其特征在于，所述获取目标对象的地震记录和衰减子波，包括：根据地震子波、衰减函数和反射系数的褶积结果得到所述地震记录，根据所述地震子波和所述衰减函数的合并结果得到所述衰减子波；根据下式计算所述地震记录：x(t)＝w(t)*a(t)*r(t)其中，x(t)为所述地震记录，w(t)为所述地震子波，a(t)为所述衰减函数，r(t)为所述反射系数；根据下式计算所述衰减子波：wa(t)＝w(t)*a(t)其中，wa(t)为所述衰减子波，w(t)为所述地震子波，a(t)为所述衰减函数。3.根据权利要求2所述的高精度反射系数求取方法，其特征在于，所述将所述地震记录和所述衰减子波转换到频域中，得到频域中的地震记录、频域中的衰减子波，包括：根据下式计算所述频域中的地震记录：x(f)＝w(f)a(f)r(f)其中，x(f)为所述频域中的地震记录，w(f)为频域中的地震子波，a(f)为频域中的衰减函数，r(f)为频域中的反射系数；根据下式计算所述频域中的衰减子波：wa(f)＝w(f)a(f)其中，wa(f)为所述频域中的衰减子波，w(f)为频域中的地震子波，a(f)为频域中的衰减函数。4.根据权利要求3所述的高精度反射系数求取方法，其特征在于，所述根据所述频域中的地震记录和所述频域中的衰减子波，得到频域中的反射系数，包括：根据下式计算所述频域中的反射系数：其中，r(f)为所述频域中的反射系数，x(f)为所述频域中的地震记录，wa(f)为所述频域中的衰减子波。5.根据权利要求1所述的高精度反射系数求取方法，其特征在于，所述对所述频域中的地震记录、所述频域中的衰减子波和所述频域中的反射系数做伽伯变换，得到伽伯域中的地震记录、伽伯域中的衰减子波和伽伯域中的反射系数，包括：根据下式计算所述伽伯域中的地震记录：其中，G[x](t,f)为所述伽伯域中的地震记录，为对所述地震记录做小时窗傅里叶变换；根据下式计算所述伽伯域中的衰减子波：G[wa](t,f)＝Gn[w(f)a(f)]其中，G[wa](t,f)为所述伽伯域中的衰减子波，w(f)为频域中的地震子波，a(f)为频域中的衰减函数；根据下式计算所述伽伯域中的反射系数：其中，G[r](t,f)为所述伽伯域中的反射系数，G[x](t,f)为所述...
【专利技术属性】
技术研发人员：柳倩男，彭苏萍，赵惊涛，崔晓芹，
申请(专利权)人：中国矿业大学北京，
类型：发明
国别省市：北京,11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。