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

说明书

本发明涉及一种基于浮动基准面的山区地震资料静校正量计算方法，属于地震资料静校正处理技术领域。本发明包括如下步骤：获取初至走时数据和近地表速度模型；计算风化层底界面和浮动基准面；计算地震波从浮动基准面传播到风化层底界面的垂向旅行时间；计算地震波的剩余走时；选择最小偏移距和基准点；计算相邻检波点间的延迟时间；计算相邻炮点间的延迟时间；计算检波点和炮点到基准点的延迟时间；计算炮点端和检波点端的静校正量；输出计算的静校正量。本发明解决了浮动基准面以上介质速度估计不准对静校正量计算的影响，且避免了使用近偏移距初至走时数据计算静校正量。

技术领域

本发明涉及一种基于浮动基准面的山区地震资料静校正量计算方法，属于地震资料静校正处理技术领域。

背景技术

在近地表条件复杂的地区，静校正问题十分突出。静校正的效果直接影响到最终的叠加剖面的质量。现有的静校正量计算方法将反演的近地表速度模型中风化层底界面以上介质的速度变化和地表起伏对地震波走时的影响作为长波长静校正量，然后利用反射波或初至波剩余走时来计算短波长静校正量。

在地表高程剧烈变化的山区，近地表浅部地层的速度往往由于观测系统中检波点间距过大和相邻检波点间高程变化剧烈而估计不准。利用近地表浅部地层速度不准的速度模型计算得到的长波长静校正量存在较大的误差，会进一步地影响短波长静校正量的准确性，最终使得静校正后的地震剖面的成像质量较低。

发明内容

针对现有技术存在的上述缺陷，本发明提出了一种基于浮动基准面的山区地震资料静校正量计算方法，其解决了在近地表条件复杂地区，现有方法静校正量计算不准确的问题。

本发明是采用以下的技术方案实现的：本发明所述一种基于浮动基准面的山区地震资料静校正量计算方法，包括如下步骤：

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

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

，（1）

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

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

，（2）

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

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

，（3）

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

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

，（4）

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