[发明专利]一种地表主动源反射波干涉成像方法

权利要求书

1.一种地表主动源反射波干涉成像方法，其特征在于，所述方法包括以下步骤：

步骤(1)、获取共成像点道集地震数据；

步骤(2)、获取基于一致性谱相关干涉重构的共成像点虚反射道集地震数据；

步骤(3)、获取基于互相关干涉重构的共成像点道集地震数据；

步骤(4)、获取基于水平叠加的目标层反射波成像结果；

所述步骤(1)包括：

步骤(1.1)、检查地震野外班报、SPS文件和地震原始数据；

步骤(1.2)、将地震野外班报和SPS文件中的观测系统信息加载到地震原始数据中，计算生成共中心点号、偏移距和方位角的观测系统信息，获取具有观测系统信息的地震数据；

步骤(1.3)、利用加载完观测系统信息的地震数据，采用自动加人工修正的方法拾取初至波走时信息，并利用初至波走时信息开展折射静校正处理，计算折射静校正后的地震数据；

步骤(1.4)、利用折射静校正后的地震数据，采用自适应噪声压制方法开展随机噪声、异常极值、侧面反射波、绕射波和地震面波的噪声压制处理，生成噪声压制处理后的地震数据；

步骤(1.5)、利用噪声压制处理后的地震数据，开展地表一致性振幅补偿处理，弥补震源点和接收点产生的振幅非均匀性，获取振幅补偿后的地震数据；

步骤(1.6)、利用振幅补偿后的地震数据，开展地表一致性反褶积和脉冲反褶积处理，生成反褶积处理后的地震数据；

步骤(1.7)、利用反褶积处理后的地震数据，开展分频去噪处理，压制步骤(1.6)中产生的高频噪声干扰，获取高频噪声压制后的地震数据。

2.根据权利要求1所述的一种地表主动源反射波干涉成像方法，其特征在于，所述步骤(2)包括：

步骤(2.1)、利用高频噪声压制后的地震数据，针对目标层开展高精度的速度分析，获取目标层速度场数据；

步骤(2.2)、利用步骤(1.7)生成的高频噪声压制后的地震数据和步骤(2.1)生成的目标层速度场数据，开展局部拉伸清零的动校正处理，获取局部拉伸清零动校正处理后的地震数据；

步骤(2.3)、利用局部拉伸清零动校正处理后的地震数据，通过水平叠加，获取目标层的反射波场数据；

步骤(2.4)、利用步骤(1.7)生成的地震数据与步骤(2.3)获取的反射波场数据，开展一致性谱相关地震干涉处理，以求取代表动校正量的虚反射道集。

3.根据权利要求2所述的一种地表主动源反射波干涉成像方法，其特征在于，所述步骤(2.4)中求取代表动校正量的虚反射道集的公式为：

其中，*代表了复共轭，ω代表了角频率，α代表了稳定因子，R_ZERO(ω)代表了步骤(2.3)生成的目标层位反射波场数据，G(ω)代表了步骤(1.7)生成的高频噪声压制后的地震数据，V(ω)代表了步骤(2.4)重构的虚反射道集。

4.根据权利要求3所述的一种地表主动源反射波干涉成像方法，其特征在于，所述步骤(3)具体为：利用步骤(2.4)生成的虚反射道集与步骤(1.7)获取的高频噪声压制后的地震数据，开展基于互相关的地震干涉处理，以重构相位校正后的共成像点道集地震数据。

5.根据权利要求4所述的一种地表主动源反射波干涉成像方法，其特征在于，所述步骤(3)中求取相位校正后的共成像点道集地震数据的公式为：

C(ω)＝G(ω)V(ω)^*

其中，C(ω)代表了相位校正后的共成像点道集地震数据。

6.根据权利要求4或5所述的一种地表主动源反射波干涉成像方法，其特征在于，所述步骤(4)具体为：针对步骤(3)生成的相位校正后的共成像点道集地震数据，开展水平叠加处理，计算生成针对目标层的反射波成像结果。

7.根据权利要求6所述的一种地表主动源反射波干涉成像方法，其特征在于，所述步骤(4)中求取针对目标层的反射波成像结果的公式为：

I(ω)＝∑C(ω)

其中，∑代表了叠加求和，I(ω)代表了通过水平叠加生成的成像结果。