[发明专利]一种在炮检距域中求取折射剩余静校正量的方法

说明书

技术领域

本发明涉及地球物理勘探静校正方法，是一种利用初至时间，在基础静校正量的基础上求取折射剩余静校正量的方法。

背景技术

地震勘探中静校正处理主要是解决地震波在近地表介质旅行时间变化影响反射波叠加成像问题。初至波反映近地表地层变化，利用初至波的折射静校正方法或层析静校正方法可以反演出相对精确的近地表模型，解决大部分长波长和中长波长的静校正问题，在解决短波长静校正方面，还存在明显不足。短波长静校正直接影响地震剖面成像的质量，是非常重要的。为了进一步提高短波长静校正量的精度，大都采取了求取剩余静校正量的方法。从使用数据的类型和原理上来分，常规求取剩余静校正量的方法有两个：一个是利用反射波信息和最优叠加反射波同相轴的原理求取地震道的剩余静校正量，再分解该地震道的剩余静校正量，得到激发点和接收点的剩余静校正量。该方法需要精度较高的叠加速度场，而得到精度较高的叠加速度场本身就是一个比较困难的问题，尤其在复杂的近地表区域，静校正量严重影响叠加速度的分析结果。叠加速度分析又是一个交互迭代过程，所以用该方法计算剩余静校正量的速度较慢；另一个方法是利用初至信息或折射波信息，在共激发点域、共接收点域和共中点域中，利用拟合初至波成一条光滑的曲线原理来求取激发点和接收点的剩余静校正量。由于使用的数据仅仅是初至和静校正量，所以该方法计算速度较快。但是该方法要求初至是连续的，在应用静校正量后，相邻地震道的初至不能出现较大跳动。为了避免相同炮检距的初至来自不同的层或折射界面倾斜的情况，在计算中要求对初至进行区分左右支，对每支的初至进行多项式拟求取剩余静校正量。该方法是在共激发点域、共接收点域或共中点域中进行的，比较适合二维勘探，在三维勘探中效果不明显。

发明内容

本发明目的在于提供一种利用初至信息，适合二维勘探和三维勘探，进一步提高短波长静校正量的精度，快速求取在炮检距域中求取折射剩余静校正量的方法。

本发明通过以下步骤实现：

1）采集地震数据，抽取共炮检距道集G：

所述的抽取共炮检距道集是在炮检距D和浮动范围d以内，抽取所有初至和相应炮检对的中点坐标，形成一个共炮检距道集G。

2）应用基础静校正量：

所述的应用基础静校正量是在共炮检距道集G中，把已知激发点和接收点基础静校正量的和，应用到初至上，形成时间t1。

3）计算炮检对剩余静校正量：

所述的计算炮检对剩余静校正量是在炮检距道集G中，以一个炮检对的中点坐标为圆心，以平滑半径R形成一个圆，计算道集G中所有炮检对中点坐标在该圆中的时间t1的平均值，形成时间t2，时间t1与时间t2时差就是该炮检对的剩余静校正量dt；

所述的时间t2利用以下公式计算：

t2=1NΣiNt1i---(1)]]>

式中：t1_i是指在道集G中，其它炮检对中点坐标到当前炮检对中点坐标的距离小于平滑半径R的时间t1；

N是相应的炮检对的总个数；i是序号；

所述的平滑半径R取最大炮检距的一半；

所述的平滑半径R取1000-2000m；

4）分解炮检对剩余校正量：

所述的分解炮检对剩余校正量是：按照激发点的比例系数a和接收点的比例系数，分解炮检对的剩余校正量dt到相应的激发点ds和接收点dr中。

步骤4）中所述比例系数a取值范围0至1；相应的接收点的比例系数是1减a，取值范围是1至0。

步骤4）中所述激发点ds=dt*a；接收点dr=dt*(1-a)。

5）采用另一个炮检距D+2d，重复迭代步骤1）到4），直到计算完所有炮检距；得到激发点或接收点在多个炮检距中的不同的剩余静校正量；

6）利用以下公式计算所有炮检对剩余校正量dt的均方根误差rms：