一种在炮检距域中求取折射剩余静校正量的方法

摘要

本发明是地球物理勘探中利用初至时间在基础静校正量的基础上求取折射剩余静校正量的方法。用基础静校正量计算所有炮检距的激发点或接收点在多个炮检距中的不同的剩余静校正量，计算分解炮检对剩余静校正量，计算所有炮检对剩余校正量的均方根误差；计算一个激发点和接收点在共炮检距域中的最终剩余静校正量，用激发点和接收点的剩余静校正量到已知的基础静校正量上，重复迭代直到均方根误差值达到给定范围完成共炮检距域中折射剩余静校正。本发明避免了常规在共激发点域、共接收点域和共中心点域计算剩余静校正量的不足，提高了短波长静校正的精度，改善长波长静校正的精度，适应二维勘探和三维勘探。

主权项

一种在炮检距域中求取折射剩余静校正量的方法，特点是通过以下步骤实现：1)采集地震数据，抽取共炮检距道集G：所述的抽取共炮检距道集G是在炮检距D和浮动范围d以内，抽取所有初至和相应炮检对的中点坐标生成；2)应用基础静校正量：3)计算炮检对剩余静校正量：所述的计算炮检对剩余静校正量是在炮检距道集G中，以一个炮检对的中点坐标为圆心，以平滑半径R形成一个圆，计算道集G中所有炮检对中点坐标在该圆中的时间t1的平均值，形成时间t2，时间t1与时间t2时差就是该炮检对的剩余静校正量dt；所述的时间t2利用以下公式计算：$t2=\frac{1}{N}\underset{i}{\overset{N}{\Sigma }}t{1}_i---\left(1\right)$式中：t1_i是指在道集G中，其它炮检对中点坐标到当前炮检对中点坐标的距离小于平滑半径R的时间t1；N是相应的炮检对的总个数；i是序号；4)分解炮检对剩余校正量：5)采用另一个炮检距重复迭代步骤1)到4)，直到计算完所有炮检距；得到激发点或接收点在多个炮检距中的不同的剩余静校正量；所述的另一个炮检距是炮检距D加两倍的浮动范围d；6)利用以下公式计算所有炮检对剩余校正量dt的均方根误差rms：$rms=\sqrt{\frac{1}{N}\underset{i}{\overset{N}{\Sigma }}{\mathrm{dt}}_i^2}---\left(2\right)$式中：N是步骤1)到5)中所有道集G中所有炮检对的总个数；i是序号；7)利用以下公式计算一个激发点在共炮检距域中的最终剩余静校正量ts：$ts=\frac{1}{Ns}\underset{i}{\overset{Ns}{\Sigma }}{\mathrm{ds}}_i---\left(3\right)$式中：Ns是同一个激发点在步骤1)到4)中产生剩余静校正量的总个数；ds_i是激发点剩余静校正量；i是序号；利用以下公式计算一个接收点在共炮检距域中的最终剩余静校正量tr：$tr=\frac{1}{Nr}\underset{i}{\overset{Nr}{\Sigma }}{\mathrm{dr}}_i---\left(4\right)$式中：Nr是同一个接收点在步骤1)到4)中产生剩余静校正量的总个数；dr_i是接收点剩余静校正量；i是序号；8)应用激发点和接收点的剩余静校正量叠加到已知的基础静校正量上，重复迭代步骤1)到7),直到前后两次迭代的rms差值达到给定范围，完成在共炮检距域中的折射剩余静校正量。