[发明专利]一种基于浮动基准面的山区地震资料静校正量计算方法

权利要求书

1.一种基于浮动基准面的山区地震资料静校正量计算方法，其特征在于，包括如下步骤：

步骤一：获取初至走时数据和近地表速度模型；

步骤二：计算风化层底界面和浮动基准面；风化层底界面的计算公式如下：

其中，z_down表示风化层底界面的高程，x表示水平方向的坐标值，w_x表示积分半径，z_m表示近地表速度模型中替换速度所在的层界面的高程，x’表示在积分区间内的坐标值；

浮动基准面的计算公式如下：

其中，z_up表示浮动基准面的高程，x表示水平方向的坐标值，w_x表示积分半径，z_surface表示地表的高程，x’表示在积分区间内的坐标值，Δz是一个修正量，其作用是使浮动基准面的高程小于地表的高程；

步骤三：计算地震波从浮动基准面传播到风化层底界面的垂向旅行时间，计算公式如下：

其中，t₁表示地震波从浮动基准面传播到风化层底界面的垂向旅行时间，x表示水平方向的坐标值，z表示高程，z_up表示浮动基准面的高程，z_down表示风化层底界面的高程，v(x,z)表示风化层在空间(x,z)处的速度；

步骤四：计算地震波的剩余走时；其计算公式可以表示为：

t_rest(s,r)＝t_raw(s,r)-t₁(x_s)-t₁(x_r)，(4)

其中，s表示炮点序号，r表示检波点序号，x_s表示炮点在水平方向的坐标值，x_r表示检波点在水平方向的坐标值，t_rest表示地震波的剩余走时，t_raw表示步骤一中获取的初至走时，t₁表示地震波从浮动基准面传播到风化层底界面的垂向旅行时间；

步骤五：选择最小偏移距和基准点；从初至走时数据中识别刚产生折射波的偏移距,将该偏移距作为最小偏移距，并将偏移距小于该最小偏移距的初至走时数据剔除，使之不参与后续步骤的计算；将高程最小的炮点或检波点作为基准点；

步骤六：计算相邻检波点间的延迟时间；

步骤七：计算相邻炮点间的延迟时间；

步骤八：计算检波点到基准点的延迟时间；其计算公式如下：

其中，r表示检波点序号，t₂(r)表示第r个检波点到基准点的延迟时间，r₀表示距离基准点最近的检波点的序号，i表示r₀和r-1之间的检波点序号，Δt_i表示第i+1和第i个检波点间的延迟时间；

步骤九：计算炮点到基准点的延迟时间；其计算公式如下：

其中，s表示炮点序号，t₂(s)表示第s个炮点到基准点的延迟时间，s₀表示距离基准点最近的炮点的序号，j表示s₀和s-1之间的炮点序号，Δt_j表示第j+1和第j个炮点间的延迟时间；

步骤十：计算炮点端和检波点端的静校正量，炮点端的静校正量的计算公式如下：

其中，s表示炮点序号，x_s表示炮点在水平方向的坐标值，z_down表示风化层底界面的高程，z_base表示基准面的高程，v_m表示替换速度，z表示高程，t_static表示静校正量，t₁表示地震波从浮动基准面传播到风化层底界面的垂向旅行时间，t₂(s)表示炮点到基准点的延迟时间；

检波点端的静校正量的计算公式如下：

其中，r表示检波点序号，x_r表示检波点在水平方向的坐标值，z_down表示风化层底界面的高程，z_base表示基准面的高程，v_m表示替换速度，z表示高程，t_static表示静校正量，t₁表示地震波从浮动基准面传播到风化层底界面的垂向旅行时间,t₂(r)表示检波点到基准点的延迟时间；

步骤十一：输出静校正量；

进一步地，所述步骤一中，近地表速度模型由初至走时层析方法获得。