延拓；
[0048] 延拓步长的确定方法如下：
[0049] 设高速层顶界面高程为H。，各个检波器高程为H1，首先对各检波器高程进行排序， 高程排序后为H1, H2.........Hni，则各延拓步长分别为Hni-H0, Hni I-Hni.........H2-H1。本步骤 能够自动分析所需延拓步长，以最高效的延拓方式完成数据变换。
[0050] 5)对延拓后的数据按照检波点位置进行分选，形成共检波点数据集。
[0051] 6)对共检波点数据集按照3)，4)步骤进行处理，处理完成后，重新选排为共炮道 集数据。
[0052] 该方法为适用于纵横向剧烈变速的非一致性波动方程延拓静校正方法，适用于解 决复杂构造下非一致性静校正问题。本方法使用层析反演获得的速度场信息，利用原始单 炮地震数据进行波动方程延拓，将地表采集的地震数据延拓到某一个基准面上。从共炮集 记录开始使用有限差分法进行延拓，从地表延拓到高速层顶界面，然后采用相移法从高速 层顶界面延拓到固定基准面上。将延拓后的数据进行道集抽取，抽取成共检波点道集数据， 然后再采用有限差分法和相移法对共检波点道集记录进行波场延拓。完成共检波点数据延 拓后将数据重新抽取为共炮道集数据。在进行波场延拓的过程中，加入随速度场复杂度自 动调整的速度场平滑参数，使用经过平滑处理后的速度场进行延拓，可以解决有限差分法 在大倾角下准确度低的问题。该方法经过模型和实际资料验证，计算效果稳定可靠。
[0053] 实施例2。一种自适应波动方程波场延拓静校正方法，采用以下具体的实施步骤：
[0054] 1)建立起伏地表模型资料，模型近地表高程差超过500m，并伴有严重横向变速， 横向速度从1500m/s到3000m/s，如图1所示。
[0055] 2)针对图1的模型进行正演模拟，观测系统参数为：中间放炮，每炮300道，道距 30米，最小偏移距为0,炮距50米，获得的单炮记录如图2所示。
[0056] 3)对图2的单炮记录进行波动方程的波场延拓静校正处理，得到静校正后的单炮 记录如图3所示。
[0057] 4)采用层析方法获得某实际资料的近地表速度模型，如图4所示。
[0058] 5)选取实际地震资料，地震资料时间长度为5000ms，时间采样间隔为2ms，共1000 道，如图5所示。
[0059] 6)应用本发明的波场延拓静校正技术对图5的原始单炮进行处理，得到处理后的 单炮，如图3所示。
[0060] 7)选取实际地震资料叠加剖面，结合图4的近地表速度模型对实际资料进行层析 静校正，静校正后叠加如图7所示。
[0061] 8)选取实际地震资料叠加剖面，结合图4的近地表速度模型对实际资料进行本发 明的波场延拓静校正，静校正后叠加如图8所示。
[0062] 对比图2与图3,自适应波动方程波场延拓静校正方法处理后，单炮各层的反射同 相轴都很好的恢复了双曲线形态，与理论双曲线吻合度高，与此同时，通过波场延拓消除了 一些噪首，资料保幅性提商。
[0063] 对比图5与图6,经过自适应波动方程波场延拓静校正方法处理后，原始单炮初至 及同向轴的扭曲现象得到了改善，资料信噪比略有提高。对比图7与图8,相比于常规的层 析静校正技术，在相同近地表速度模型的前提下，波场延拓静校正叠加后同相轴更加连续， 构造形态更加自然，而在相同位置，层析静校正同相轴凌乱，表现出静校正问题比较严重。
【主权项】
1. 一种自适应波动方程波场延拓静校正方法，其特征在于：包括以下步骤： 步骤1:采用原始的单炮地震数据； 步骤2 :使用常规的初至波层析方法获得近地表的速度结构； 步骤3 :采用有限差分法将地表获得的波场延拓到高速层顶界面 常规有限差分法延拓公式为：式中：P为波场值，v为速度，a_j、bj为系数，j= 1,2,…，n; 改进后的有限差分法延拓方程为：式中^r为模型自适应平滑速度，在进行;7计算前，需用对当前选取速度点三维空间 一定范围内速度进行分析，获得其横向变化系数r= | (Vo-vD/vol，当r〈0. 2时，不修改原速 度场，当0. 2彡r时，按照w= 5* (r-0. 2)计算平滑范围； 采用改进后的延拓方程，在延拓的过程中判断周围速度场的变化程度，自动进行平滑 处理，使有限差分法在遇到速度突变时造成的空间假频现象得到压制； 步骤4 :采用相移法将波场从高速层顶界面延拓到地表，对数据进行延拓时采用变步 长延拓法进行延拓； 步骤5 :对延拓后的数据进行分选，按照检波点选排为共检波点数据集； 步骤6 :再采用有限差分法和相移法对共检波点道集记录进行波场延拓；处理完成后， 重新选排为共炮道集数据。2. 根据权利要求1所述的自适应波动方程波场延拓静校正方法，其特征在于：所述变 步长延拓法的具体步骤是： 首先分析高速层顶界面到地表的高差，根据各高差差值设计延拓步长；然后分步延 拓； 延拓步长的确定方法如下： 设高速层顶界面高程为H。，各个检波器高程为私，首先对各检波器高程进行排序，高程 排序后为氏，H2.........Hm，则各延拓步长分别为Hm_H。，HmfHm.........氏-氏。3. 根据权利要求1所述的自适应波动方程波场延拓静校正方法，其特征在于：步骤5 中，所述对延拓后的数据进行分选，按照检波点选排为共检波点数据集的具体步骤为：将延 拓后的数据进行道集抽取，抽取成共检波点道集数据。4. 根据权利要求1所述的自适应波动方程波场延拓静校正方法，其特征在于：在步骤6 中，按照步骤3所述的常规的有限差分法延拓公式和步骤4所述的相移法对共检波点道集 记录进行波场延拓。5.根据权利要求1所述的自适应波动方程波场延拓静校正方法，其特征在于：在步骤6 中，按照步骤3所述的改进后的有限差分法延拓公式和步骤4所述的相移法对共检波点道 集记录进行波场延拓。
【专利摘要】本发明公开了一种自适应波动方程波场延拓静校正方法。使用层析反演获得的速度场信息，利用原始单炮地震数据进行波动方程延拓，将地表采集的地震数据延拓到某一个基准面上。从共炮集记录开始使用有限差分法进行延拓，从地表延拓到高速层顶界面，然后采用相移法从高速层顶界面延拓到固定基准面上。将延拓后的数据进行道集抽取，抽取成共检波点道集数据，然后再采用有限差分法和相移法对共检波点道集记录进行波场延拓。完成共检波点数据延拓后将数据重新抽取为共炮道集数据。在进行波场延拓的过程中，加入随速度场复杂度自动调整的速度场平滑参数，使用经过平滑处理后的速度场进行延拓，可以解决有限差分法在大倾角下准确度低的问题。
【IPC分类】G01V1/36
【公开号】CN105093318
【申请号】CN201410203113
【发明人】冮明川, 韩站一, 马晓义, 王延光, 单联瑜, 王兴谋, 赵胜天, 胡晓婷, 孙天真, 潘树林, 孙兴刚
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司胜利油田分公司物探研究院
【公开日】2015年11月25日
【申请日】2014年5月14日