速度总是在相差 约0. 3m/s的两个速度之间来回跳变，这显然与实际地震资料采集时采集船基本保持匀速 航行的情况不符。这个问题可以通过从P2/94导航文件中提取相应的放炮时间信息得到改 善，P2/94导航文件是P1/90导航文件的基础文件，但它记录的放炮时间精度达到了 0. 1秒。 P2/94文件是同样是英国勘探委员会UKOOA制定的海上地震勘探导航格式文件，用ASCII码 编写，每行80个字符。按照P2/94导航文件格式规定：每行以E1000开头的数据记录了放 炮时间?目息，其中第7 _22列是炮线名，第24-31列是放炮编号，第33-48列是地震记录标识 符，第50-57列是放炮日期，第59-66列是放炮时间，第68-70列是震源识别号。图5是按图 3所示的同一炮线的Ρ2/94导航文件计算的放炮时间间隔，数值主要分布在7. 2~8. 0秒之 间的以〇. 1秒间隔的条带上，分选更好。图6是用图5中的放炮时间间隔求得的采集船航 速，数值大致从2. 38m/s逐渐变化到2. 58m/s，相邻两炮船速的变化减小，更符合实际情况。
[0038] 2)选择该炮线上的某一个炮点，对采集船所采集到的该炮点激发所获得的地震数 据按照采集缆进行分组。
[0039] 3)选择一个采集缆，计算炮点到采集缆上各接收点的偏移距，以最小偏移距接收 点为原点，计算该采集缆各接收点在炮点激发时刻即t = 0时的弧长坐标馬，i?f，···， ：共n+1个接收点。
[0040] 其中，"炮点到接收点的偏移距"指的就是炮点到接收点的距离，这是行业内的惯 用叫法。
[0041] 如7所示，图中曲线代表海面上的一条采集缆，曲线上的圆圈代表采集缆上的接 收器，当采集缆受采集船牵引移动时，采集缆上的每个接收器相互牵引，沿图中箭头的方向 移动，因此，可以定义炮点激发时（t = 0)接收点的弧长坐标：以采集缆上靠近采集船一端 的首个接收点为原点，原点处的接收点的弧长坐标为〇,记作笔，而采集缆上除原点外的任 一接收点k的弧长坐标起为：
[0043] 式中，JTj表示第j个接收点到第j-Ι个接收点的弧长；X 分别表示第j个接收 点的横坐标和纵坐标；Xj i分别表示第j-Ι个接收点的横坐标和纵坐标。由于放炮时海 面上采集缆的形态难以精确描述，实际应用中以两点间的直线距离代替弧长。因此，在t = O时刻采集缆上接收点k的位置可以描述为：采集缆上距离原点的折线距离为 < 的点。
[0044] 4)计算t时刻采集缆上各接收点相对于炮点激发时弧长坐标原点的弧长坐标^ :
[0046] 式中，V表示采集船航行速度，其取值由公式（2)获得，即先确定该炮点为所选炮 线上的第几个炮，然后由公式（2)计算出该炮数据采集时采集船的平均航行速度，此时的V 即等于该平均航行速度；t表示当前地震数据的记录时间，其中k = 0, 1，2,…，η。
[0047] 5)以弧长坐标耗和t时刻接收点上的地震数据为原始样本，样条插值if上的地震 波场值（k = 0, 1，2,…，η)，此插值结果即为t时刻位移校正后的地震数据；
[0048] 6)按照记录时间间隔逐渐改变t值，重复实施步骤4)和步骤5)，直到记录时间结 束，由此获得在步骤2)所选的炮点下，对步骤3)所选的采集缆上各接收点的位移校正。
[0049] 类似地，参照上述过程，可以依次完成：在步骤2)所选的炮点下对全部采集缆上 各接收点的位移校正，在步骤1)所选的炮线上的各炮点下对全部采集缆上各接收点的位 移校正，以及在全部炮线上的各炮点下对全部采集缆上各接收点的位移校正。
[0050] 本申请还对上述位移校正方法的效果进行了验证：海上时移地震数据经过接收点 位移校正以后，消除了因两次采集时船速不同带来的成像误差影响，使得资料的一致性更 好。为了检测两次成像结果的可重复性，需要对标志层位置的时窗内的成像结果做可重复 性测量。这些检测手段包括计算频谱、相关系数、时移误差、振幅比和均方根振幅差等。目 前最常用的重复性测量方法是计算NRMS (Normalized Root-Mean-Square difference，归 一化均方根差异分析）。NRMS值是两次成像的振幅差的平均均方根振幅除以两次成像的平 均均方根振幅和，BP
[0052] 其中，B是前一次采集资料所成的像，M是后一次采集资料所成的像，rms算子定义 如下
[0054] 其中，X1是成像剖面时窗里的振幅，N是时窗里的样点个数。
[0055] 时移地震用NRMS值来检测两次处理结果的可重复性。如果NRMS = 0,表示数据 完全一致，如果NRMS = 2,说明数据反相关，如果B和M只包含随机噪音，NRMS=n/I。NRMS 值受相位和振幅差、时移误差和噪音的影响，值越小说明可重复性越好。
[0056] 图8是海上某实际工区2003年和2013年两次时移地震采集分布范围，选取两条 相邻的主测线分别记作2003-79线和2013-08线，对其叠加剖面做可重复性测量，检验地 震资料接收点位移校正的应用效果。图9和图10分别是2003-79线和2013-08线的叠加 剖面，其中左边是没有做接收点位移校正的原始数据的叠加剖面，中间是校正后数据的叠 加剖面，右边是校正前后的差异剖面。设计3个不同时窗检测校正前后的NRMS值，其中 1. 2~I. 5s时窗的NRMS值校正前为0. 582531，校正后降为0. 582055 ;1. 5~I. 7s时窗的 NRMS值校正前为0. 566818,校正后降为0. 462104,1. 9~2. 2s的时窗则从0. 603000降为 0. 600378。可见在保持相同处理参数的情况下，时移地震数据经接收点位移校正后NRMS值 减小，表明本发明提高了时移地震两次采集资料的一致性。
[0057] 上述各实施例仅用于说明本发明，其中方法的实施步骤等都是可以有所变化的， 凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围 之外。
【主权项】
1. 一种海上时移地震接收点位移校正方法，包括以下步骤： 1) 选择一条炮线，计算采集船在采集该炮线上各炮的地震数据时其自身的平均航行速 度； 2) 选择该炮线上的某一个炮点，对采集船所采集到的该炮点激发所获得的地震数据按 照采集缆进行分组； 3) 选择一个采集缆，计算炮点到采集缆上各接收点的偏移距，以最小偏移距接收点为 原点，计算该采集缆各接收点在炮点激发时刻即t= 0时的弧长坐标<，<，…，<，···， C，共n+1个接收点； 4) 计算t时刻采集缆上各接收点相对于炮点激发时弧长坐标原点的弧长坐标巧： 5) 以弧长坐标/?丨和t时刻接收点上的地震数据为原始样本，样条插值及；1上的地震波场 值，此插值结果即为t时刻位移校正后的地震数据，其中k= 0, 1，2,…，η; 6) 按照记录时间间隔逐渐改变t值，重复实施步骤4)和步骤5)，直到记录时间结束， 由此获得在步骤2)所选的炮点下，对步骤3)所选的采集缆上各接收点的位移校正。2. 如权利要求1所述的一种海上时移地震接收点位移校正方法，其特征在于：所述步 骤1)中，对于包括N个炮点的炮线，由下式计算出前N-1个炮数据采集时采集船的平均航 行逋麽.式中，t表示第i炮数据采集时采集船的平均航行速度；?\表示第i炮的激发时间；T1+1 表示第i+1炮的激发时间分别表示第i炮激发时采集船的横坐标和纵坐标；x1+1、y1+1 分别表示第i+1炮激发时采集船的横坐标和纵坐标；其中，i= 1，2,...，N-1 ; 第N炮采集时采集船的平均航行速度VN取值为VNρ3. 如权利要求1所述的一种海上时移地震接收点位移校正方法，其特征在于：所述步 骤3)中，所选的采集缆上的各接收点在炮点激发时刻即t= 0时的弧长坐标if的计算公式 为-式中，h表示第j个接收点到第j-Ι个接收点的弧长；X 分别表示第j个接收点的 横坐标和纵坐标；Xji、yji分别表示第j-Ι个接收点的横坐标和纵坐标。4. 如权利要求1所述的一种海上时移地震接收点位移校正方法，其特征在于：所述步 骤4)中，t时刻采集缆上各接收点相对于炮点激发时弧长坐标原点的弧长坐标每的计算公 式为： /?；=/<：'-Vi 式中，V表示采集船航行速度，其值取自步骤1)，即先确定该炮点为所选炮线上的第几 个炮，然后根据步骤1)计算出该炮数据采集时采集船的平均航行速度赋予V值；t表示当 前地震数据的记录时间，其中k= 0, 1，2,…，η。
【专利摘要】本发明涉及一种海上时移地震接收点位移校正方法，包括以下步骤：1)选择一条炮线，计算采集船采集各炮数据时自身的平均航行速度；2)选择该炮线上的某炮点，对采集船所采集到的该炮点激发所获得的地震数据按照采集缆分组；3)选择一个采集缆，计算炮点到采集缆上各接收点的偏移距，以最小偏移距接收点为原点，计算该采集缆各接收点在炮点激发时刻的弧长坐标；4)计算t时刻采集缆上各接收点相对于炮点激发时弧长坐标原点的弧长坐标；5)以弧长坐标和t时刻接收点上的地震数据为原始样本，样条插值上的地震波场值，此插值结果即为t时刻位移校正后的地震数据；6)按照记录时间间隔逐渐改变t值，重复实施步骤4)和步骤5)，直到记录时间结束，由此获得在所选的炮点下，对所选的采集缆上各接收点的位移校正。
【IPC分类】G01V1/24, G01V1/22
【公开号】CN105259570
【申请号】CN201510770047
【发明人】朱振宇, 薛东川, 张金淼, 陈宝书, 王小六, 何洋洋, 孙文博, 陈剑军, 刘志鹏, 欧阳炀, 王清振, 桑淑云, 李丽霞, 杨俊 , 王艳冬, 汪小将
【申请人】中国海洋石油总公司, 中海油研究总院
【公开日】2016年1月20日
【申请日】2015年11月12日