[发明专利]一种提高地下岩溶探测精度的方法

权利要求书

1.一种提高地下岩溶探测精度的方法，其特征在于：所述方法包括以下步骤：

①分频处理：输入叠后纯波保幅地震数据，然后对叠后纯波保幅地震数据进行分频处理，得到一系列频率从低到高的分频数据体；

②相干数据处理：分别对所述不同的分频数据体使用相干算法进行运算，得到一系列的分频相干数据体；

③加权数据处理：对所述一系列的分频相干数据体利用最优加权因子进行加权数据处理，得到重构数据体。

2.根据权利要求1所述的提高地下岩溶探测精度的方法，其特征在于：所述步骤①包括以下步骤：

(11)开展研究区的测井资料和叠后纯波保幅地震数据的精细标定，即将岩溶发育层段标定在叠后纯波保幅地震数据剖面上，并将该岩溶发育层段作为目标段，对该目标段进行频谱分析，得到研究区的叠后纯波保幅地震数据的有效频带范围，再根据所述有效频带范围设计分频的起、止频带范围；

(12)如果井资料满足建模条件，则利用已知不同规模的岩溶地质体进行正演计算，以确定出研究区不同规模岩溶地质体的调谐频率，该调谐频率应包含在有效频带范围之内；

正演计算的公式如下：

f＝V/4H

其中，H为岩溶地质体模型的厚度，f为调谐频率，V为层速度；

如果井资料不满足建模条件，则利用反演方法求取调谐频率。通过对所述叠后纯波保幅地震数据经过相干算法计算后得到的相干数据体沿岩溶发育层位提取相干切片，选取切片平面上不同发育规模的岩溶所造成的相干值异常点范围，并对该范围的叠后纯波保幅地震数据体进行时频分析，以此统计得到不同发育规模的岩溶地质体的调谐频率；

(13)对步骤(12)中得到的所述不同发育规模的岩溶地质体的调谐频率进行统计分析得到频谱分析成果，划分出不同发育规模的岩溶地质体的调谐频率集中分布的频段，再根据所述频谱分析成果设计分频数据的参数，所述参数包括分频个数和滤波频段；

(14)应用离散傅里叶变换把所述叠后纯波保幅地震数据变换到频率域，然后按步骤(13)设计的所述参数利用频率域滤波方法生成一系列不同频段的分频数据体。

3.根据权利要求2所述的提高地下岩溶探测精度的方法，其特征在于：所述步骤②包括以下步骤：

(21)选取相干算法，然后对相干算法中的x、y方向上的视倾角α和β进行参数测试，找到能准确识别出测试区域内的岩溶形态时的α和β的值，所述测试区域为研究区的一部分；

(22)应用该α和β的值对研究区内所述的一系列不同频段的分频数据体进行相干计算，得到一系列的分频相干数据体。

4.根据权利要求3所述的提高地下岩溶探测精度的方法，其特征在于：所述步骤③包括以下步骤：

(31)求取所述分频相干数据体的加权因子；所述加权因子是通过等值法或者非等值法进行计算得到的；然后从所有加权因子中选出较好的加权因子，所述较好的加权因子是指能综合突出要预测的中、小规模岩溶地质体边界的加权因子；

(32)针对所述较好的加权因子进行覆盖目标地质体在内的小范围数据重构后求取其相对应的评价因子，然后找到最小的评价因子，所述最小的评价因子对应的加权因子即为最优加权因子；

所述评价因子的计算公式如下：

B‾=Σi=1nMi/M′in≥2]]>

式中，为评价因子，M_i为第i点重构数据体的目标地质体的相干值，M′_i为第i点叠后纯波的分频相干数据体的目标地质体的相干值，n为测试计算用到的部分数据体中的所有样点个数；

(33)对所述分频相干数据体采用步骤(32)得到的最优加权因子进行研究区数据重构得到重构数据体；

进行数据重构的计算公式如下：

Y=Σi=1nkiCi]]>

其中：Y为加权数据处理后的重构数据体；C_i为第i个加权频率相干数据体；其中k_i为第i个加权频率相干数据体的最优加权因子，且0≤k_i。