[发明专利]地表一致性和非一致性联合静校正方法

说明书

本发明涉及一种地表一致性和非一致性联合静校正方法，属于地球物理勘探的地震资料处理领域。本发明主要包括如下步骤：1)地震资料预处理；2)建立近地表速度模型；3)计算地表平均高程；4)对地表以上空气层以地表速度进行填充；5)计算炮点端和检波点端地表一致性静校正量；6)进行地表一致性静校正；7)将地震资料从时间‑空间域变换到τ‑p域；8)计算截距时间校正量；9)进行截距时间校正；10)将地震资料从τ‑p域变换到时间‑空间域。应用本发明，可以结合地表一致性和地表非一致性静校正方法的优点，克服单一方法的不足，得到准确的静校正结果。

技术领域

本发明涉及一种地表一致性和非一致性联合静校正方法，属于地球物理勘探的地震资料处理领域。

背景技术

黄土塬地区地表高程起伏剧烈，近地表低速层厚度在横向上变化极大。在黄土塬地区进行地震勘探时，受巨厚的近地表低速层的影响，地表记录的地震反射波同相轴极度扭曲，无法直接进行动校正和叠加成像。

为了降低近地表低速层对深部反射波同相轴形态的影响，基于地表一致性假设的静校正方法被广泛采用。地表一致性假设认为近地表低速层的速度远小于基岩速度，反射波在近地表低速层中沿垂直方向传播。基于地表一致性假设的静校正方法对地表同一个位置记录的地震资料应用相同的时间校正量。然而，不同深度界面的反射波在巨厚的近地表低速层中沿不同的倾斜路径传播，具有不同的旅行时间。

在τ-p域（截距时间-射线参数域），同一p值（射线参数值）对应的地震反射波在近地表低速层中具有相同的传播路径。对同一p值的地震反射波应用相同的校正量，对不同p值的地震反射波应用不同的校正量，就可以根据反射波在近地表低速层中的传播路径进行校正，而不是根据地表位置进行校正。然而，在黄土塬地区采集的地震资料中的反射波同相轴不仅受反射波传播路径的影响还受地形起伏的影响。在τ-p域直接对地震资料进行校正，会受地形的影响产生极大的误差，严重影响静校正的效果。

发明内容

针对现有技术存在的上述缺陷，本发明提出了一种地表一致性和非一致性联合静校正方法，假定反射波在地表与地表平均高程之间沿垂直路径传播，而在地表平均高程与近地表低速层底界面间沿倾斜路径传播，对垂直路径应用地表一致性静校正，对倾斜路径应用地表非一致性静校正，解决了反射波在黄土塬地区近地表低速层中沿垂直方向校正的问题和射线参数受地形起伏影响的问题。

本发明是采用以下的技术方案实现的：本发明所述一种地表一致性和非一致性联合静校正方法，包括如下步骤：

步骤一：对陆上地震资料进行预处理，预处理包括振幅均衡、去除面波、压制随机噪声；

步骤二：利用初至走时数据建立近地表速度模型，近地表速度模型采用相同大小的正方形网格离散；

步骤三：计算地表平均高程；

步骤四：对近地表速度模型中地表以上的空气层以地表速度进行填充；

步骤五：计算炮点到地表平均高程的地表一致性静校正量；炮点到地表平均高程的地表一致性静校正量计算公式为：

，（1）

式中，t_s1是炮点到地表平均高程的地表一致性静校正量，i是炮点与地表平均高程间离散网格的序号，n₁是炮点与地表平均高程间离散网格的最大序号，dz是离散网格的厚度，v_i是第i个离散网格的速度；

步骤六：计算检波点到地表平均高程的地表一致性静校正量；检波点到地表平均高程的地表一致性静校正量计算公式为：

，（2）