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

权利要求书

1.地表一致性和非一致性联合静校正方法，其特征在于，包括如下步骤：

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

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

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

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

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

，（1）

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

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

，（2）

式中，t_r1是检波点到地表平均高程的地表一致性静校正量，j是检波点与地表平均高程间离散网格的序号，m₁是检波点与地表平均高程间离散网格的最大序号，dz是离散网格的厚度，v_j是第j个离散网格的速度；

步骤七：对预处理后的地震记录进行地表一致性静校正；地表一致性静校正的计算公式为：

，（3）

式中，d₁是地表一致性静校正后的地震记录，x_s是炮点x轴坐标，x_r是检波点x轴坐标，t_sr¹是地表一致性静校正后的反射波走时，d是预处理后的地震记录，t_sr是地表一致性静校正前的反射波走时，t_s1是炮点到地表平均高程的地表一致性静校正量，t_r1是检波点到地表平均高程的地表一致性静校正量；

步骤八：将地表一致性静校正后的地震资料从时间-空间域变换到τ-p域；采用下式进行变换：

，（4）

式中，m是变换后的地震资料，τ是截距时间，p是射线参数，d₁是地表一致性静校正后的地震记录，x_s是炮点x轴坐标，x_r是检波点x轴坐标，t_sr¹是地表一致性静校正后的反射波走时；

步骤九：计算射线参数对应的截距时间校正量；采用下式计算截距时间校正量：

，（5）

式中，Δτ是截距时间校正量，k是射线经过的近地表速度模型离散网格序号，n_z是射线经过的近地表速度模型离散网格序号的最大值，dz是离散网格的厚度，q_k是射线经过第k个网格时的垂直慢度；

步骤十：对变换后的地震资料进行截距时间校正；采用下式进行截距时间校正：

，（6）

式中，m’是截距时间校正后的地震资料，τ’是校正后的截距时间，p是射线参数，m是变换后的地震资料，τ是截距时间，Δτ是截距时间校正量；

步骤十一：将截距时间校正后的地震资料从τ-p域变换到时间-空间域；采用下式进行变换：

，（7）

式中，d’是变换后时间-空间域的地震资料，x_s是炮点x轴坐标，x_r是检波点x轴坐标，t_sr²是地表非一致性静校正后的反射波走时，m’是截距时间校正后的地震资料，τ’是校正后的截距时间，p是射线参数。