地球物理学报  2020, Vol. 63 Issue (2): 766-777   PDF    
三维OBS探测实验中炮点和OBS位置校正新方法
杨富东1,2,4, 张佳政1,2, 杜峰3, 王强1,2,4, 庞新明1,2,4, 赵明辉1,2,4, 丘学林1,2,4     
1. 中国科学院边缘海与大洋地质重点实验室, 南海海洋研究所, 广州 510301;
2. 中国科学院南海生态环境工程创新研究院, 广州 510301;
3. 湖北省地震局, 地震预警湖北省重点实验室, 武汉 430071;
4. 中国科学院大学, 地球与行星科学学院, 北京 100049
摘要:炮点和海底地震仪(ocean bottom seismometer,OBS)位置校正是三维地震数据处理的基本环节,也是获取高精度三维速度结构的关键所在.本文基于南海洋陆转换带(Continental-Oceanic-Transition zone,COT)IODP367/368钻探区开展的三维OBS深地震探测数据,开展了炮点及OBS位置校正研究,新的校正方法主要体现在三个方面:(1)利用连续三个炮点的平均航向对中间炮点进行位置校正,更真实地反映气枪枪阵与船体之间的软连接状态;(2)根据"滑动窗口"思想将海水声学速度阈值划分成N等份,通过循环测试获得全局最优的OBS位置校正结果,改进了前人方法只能获得局部最优解的问题;(3)针对单条测线穿过的OBS,通过加入其临近测线的直达水波走时,构成视双测线OBS位置校正法,提高了校正精度.49台OBS位置校正结果表明,除3台单测线法校正的OBS在垂直测线方向存在较大不确定性外,其余误差范围均为35 m左右.本文改进的OBS位置校正方法,不仅提高了单条测线穿过的OBS位置校正精度,保证后续三维地震结构研究的可靠性,而且为今后类似的OBS位置校正提供了经验和借鉴.
关键词: OBS位置校正      炮点位置校正      视双测线校正      三维速度结构      误差分析     
A new method for shots and OBSs' relocation applying in three-dimensional seismic survey
YANG FuDong1,2,4, ZHANG JiaZheng1,2, DU Feng3, WANG Qiang1,2,4, PANG XinMing1,2,4, ZHAO MingHui1,2,4, QIU XueLin1,2,4     
1. Key Laboratory of Ocean and Marginal Sea Geology, South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou 510301, China;
2. Innovation Academy of South China Sea Ecology and Environmental Engineering, Chinese Academy of Sciences, Guangzhou 510301, China;
3. Hubei Key Laboratory of Earthquake Early Warning, Hubei Earthquake Agency, Wuhan 430071, China;
4. College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
Abstract: The position correction of shots and OBSs is the basic procedure of three-dimensional (3D) seismic data processing, and also a key step in acquiring high-precision 3D velocity structure. Based on the deep seismic data of the ocean bottom seismometer (OBS) experiment carried out in the IODP 367/368 drilling area on the Continental-Oceanic-Transition zone (COT) of South China Sea (SCS), this paper improved the previous shots and OBSs' relocation program mainly with three aspects. (1) The position of the intermediate shot is corrected by the average heading angle of three consecutive shots, which more realistically reflects the soft connection between the air gun array and the hull. (2) According to the Sliding Window idea, the water wave velocity threshold was divided into N equal parts, and the global optimal OBS position correction result could be obtained through the loop test, which solved the problem that the previous method could only obtain the local optimal solution; (3)To improve the accuracy of single line OBS relocation, travel-times were picked from in-line shots and nearest off-line shots, where a quasi-crossline model was built to improve the correction accuracy. 49 OBSs' positions were finally relocated. Except that the three OBSs corrected by single line method had large uncertainty in the vertical line direction, the remaining semi-major axis of error ellipses of all the OBSs were about 35 m. This approach not only makes better use of the high-quality data of the 3D OBS experiment in the IODP drilling area, but also provides accurate data base for the subsequent acquisition of the 3D velocity structure, and can provide solution to shots and OBS relocation of air-gun OBS survey in the future.
Keywords: OBS relocation    Shot relocation    Quasi-crossline relocation method    3D Seismic velocity structure    Error analysis    
0 引言

南海作为西太平洋最大的边缘海之一,是我国走向深海研究的重要突破口(汪品先,2009),而南海西北次海盆东部发育的洋陆转换带(Continent-Ocean Transition zone,COT),记录着南海如何从陆地演变成洋盆的重要构造信息,并且发育着深水-超深水含油气盆地,是研究南海科学问题和资源问题的重要区域之一.前人研究发现,南海北部COT的岩石圈结构有着典型的三维空间变化特征,相邻平行剖面的地质构造变化剧烈(Ding et al., 2011Cameselle et al., 2015).国际大洋发现计划(IODP)367/368航次的钻探结果显示,该区基底具有同裂谷或前裂谷的砾岩,没有发现蛇蚊石化的上地幔出露,从裂谷阶段到海底扩张阶段是一个快速的演变过程(Larsen et al., 2018).为了精细解读COT张裂-破裂机制,2018年在南海COT区域第一次开展了三维OBS(Ocean Bottom Seismometer)深地震探测实验,旨在揭示南海由伸展、张裂、破裂到海底扩张的三维变化特征;其中OBS的准确海底位置是正确拾取震相的重要保障,也是最终获取高精度三维速度结构的关键所在.

在海上开展OBS地震探测实验时,OBS按设计点投放,以自由落体形式沉入海底,由于受惯性和海流等因素的影响会偏离投放点,因此,OBS位置校正是地震资料处理的首要环节.前人在求取OBS真实落点方面做了许多尝试,从最初的借助声学收发系统进行OBS定位(Creager and Dorman, 1982),逐渐演变为目前较为流行的基于直达水波走时数据的OBS位置校正方法(Nakamura et al., 1987王彦林等,2007薛彬等,2008).前人利用西南印度洋中脊、南海中央次海盆、南海东北部马尼拉俯冲前缘和马里亚纳海沟的主动源OBS数据,对OBS位置校正方法进行了不断地改进和完善.研究发现,对于交叉测线控制的OBS,其位置校正结果具有较高的精度,而单条测线控制的OBS则在垂直测线方向具有较大的不确定性(敖威等,2010张莉等,2013Du et al., 2018陈瀚等,2019).

本文针对2018年南海COT三维地震探测实验中大量OBS站位只有单条测线穿过的特殊情况(图 1),为了提高OBS位置校正结果的精度,提出了一种新的位置校正方案,即借助临近视测线(定义其为与OBS垂直距离大于5 km的测线)的直达水波走时信息,采用视双测线的校正手段替代先前精度较低的单测线校正方法,达到接近于精度较高的交叉测线校正的效果.首先以交叉测线上的OBS28为例,对新的位置校正方案进行测试,验证新方案的可行性后,将其应用于所有单测线穿过的OBS.在航时有限、多数OBS只有单条放炮测线的情况下,该视双测线OBS位置校正方法,不仅提高了OBS位置校正精度,同时为今后类似的三维OBS地震探测实验的位置校正提供了经验和借鉴.

图 1 南海北部洋陆转换带三维OBS深地震探测实验 红色圆圈为OBS设计位置,黑色直线为设计放炮测线,五角星为IODP钻孔位置. Fig. 1 3D OBS seismic experiment at COT zone in northern SCS Red circles represent designed OBS locations; black lines are designed shooting lines, yellow stars are IODP drilling sites.
1 数据采集与处理

2018年5—7月,国家基金委南海北部地球物理共享航次搭载中国科学院南海海洋研究所的“实验2”号科考船,在南海北部COT区域(114.5°E—117.5°E,17°N—20°N)成功开展三维OBS深地震探测实验(图 1).此次实验使用的是中国科学院地质与地球物理研究所研制的便携式OBS(郝天珧和游庆瑜,2011),采样率分别为100 Hz和250 Hz,共投放52台,最终成功回收49台,回收率达到94.23%;震源为中国科学院南海海洋研究所的4支大容量Bolt枪组成的枪阵,总容量达6000 in3(98.32 L).此次实验放炮测线共计11条,有效激发炮点数量为8751次,放炮间隔为90 s;放炮期间船速约为5节,炮间距约200 m.在成功回收的49台OBS中,交叉测线通过的OBS为19台(图 1),而仅有一条测线穿过的OBS达到30台,占回收总数一半以上.

室内地震数据预处理主要包括导航数据的整理和OBS记录数据的格式转换.首先,通过RAW2UKOOA程序提取原始导航数据中记录的炮号、精确的格林尼治标准时间(GMT时间)和经纬度坐标等信息,并依据UKOOA标准格式写入到炮点信息文件中,形成此次实验的炮点导航文件;其次,利用RAW2SAC程序对OBS的原始数据进行解编并转化为SAC格式,对故障数据进行修复(张佳政等,2018);然后使用地震分析软件SAC(Seismic Analysis Code,William and Joseph, 1991)进行查看,通过均衡和滤波等处理,在确认可见信噪比较高的气枪信号后,结合炮点导航文件,利用SAC2Y程序截取SAC格式数据中的有效信号段,使之转换成国际通用的SEGY格式,并利用地震处理软件SU(Seismic Unix, Stockwell and Cohen, 2002)进行可视化处理,得到综合地震记录剖面(赵明辉等,2004).

炮点位置校正和直达水波走时拾取是OBS位置校正前的关键步骤,这里先对二者进行简单介绍.

1.1 炮点位置校正

图 2a所示,气枪枪阵与差分全球定位系统(DGPS)接收点的距离是炮点位置校正的首要原因,“实验2”号船DGPS到船尾距离为17 m,连接船尾与枪阵的拖缆长35 m,真实的炮点位置应位于DGPS记录位置的后方;其次,气枪枪阵与船体之间为拖缆为软连接,船只受海流、海风等因素发生摇摆,尤其是在拐弯的时候,枪阵位置并非落在前后两炮DGPS记录位置连线的向后延伸直线上(图 2b虚线).因此,本研究对原有的炮点航向校正方法(敖威等,2010)进行改进,利用连续三炮的平均航向来校正炮点位置,从而更真实地反映气枪枪阵与船体之间的软连接状态.

图 2 炮点位置校正方法示意图 (a) “实验2”号综合地球物理科考船主要工作单元示意图;(b)连续3炮的位置示意图. Fig. 2 Schematic diagram of shot position correction method (a) Schematic diagram of major devices installed on R/V Shiyan 2; (b) Schematic diagram of continuous three shots′ positions.

具体炮点校正流程如下:以图 2b的炮点x位置校正为例,首先分别求取其与x-1和x+1的连线航向,并将其平均值作为气枪枪阵在船体后方的航向;然后依据气枪枪阵到DGPS的间距,将DGPS记录的炮点位置校正到真实的气枪枪阵中心.

1.2 直达水波走时拾取

本文利用Upicker软件(Wilcock, 2011)进行OBS直达水波走时拾取.通常情况下,近偏移距的直达水波(Pw)表现为初至震相,容易识别和拾取;但当OBS位置与放炮测线不重合时,直达水波有时会晚于其它震相,成为后至震相,此时其信噪比较低,不易识别和拾取.因此,对Upicker软件进行了改进,增加了交叉点走时和理论直达水波走时显示功能(图 3),用于辅助直达水波的识别和拾取.

图 3 OBS28台站直达水波走时拾取 (a)主测线L2直达水波走时;(b)主测线H3直达水波走时;(c)视测线H2直达水波走时;(d)主测线L2直达水波走时拾取;(e)主测线H3直达水波走时拾取;(f)视测线H2直达水波走时拾取.横轴为炮号, 纵轴为折合时间, 折合速度为3 km·s-1.红色点为直达水波到时;白色点为测线交点;绿色曲线为直达水波理论走时;蓝色点为测线交点最近两炮间距.滤波参数flow-pass=3 Hz, fhigh-pass=15 Hz. Fig. 3 The travel-time picks of direct water waves in the seismic record section of OBS28 (a) Direct water wave travel-times in main line L2; (b) Direct water wave travel-times in main line H3; (c) Direct water wave travel-times in off line H2; (d) Direct water wave travel-times picked in main line L2; (e) Direct water wave travel-times picked in main line H3; (f) Direct water travel-times time picked in off line H2. The horizontal axis is shot number, the vertical axis is the reduced time, reduced velocity is 3 km·s-1. Red dots represent shots in shooting lines; white dots are the intersection of two lines; green line is theoretical travel-times of direct water wave; blue dots are the distance between two shots that are nearest to the intersection of two shooting lines.Filters, flow-pass=3 Hz, fhigh-pass=15 Hz.

这里以OBS28为例,简单介绍其主测线L2、H3和视测线H2的直达水波走时拾取.主测线L2和H3近偏移距的直达水波为初至震相,信噪比高,极易识别和拾取;远偏移距时则埋藏于地壳折射波震相(Pg)以下,可依据初至震相外延和理论直达水波(假设OBS投放点即为其落点时计算所得)进行识别和拾取(图 3d3e).对于视测线H2,其直达水波完全埋藏于Pg震相之下,信噪比较差,此时可根据与其相交的主测线L2上识别的直达水波以及该测线的理论直达水波,快速准确地识别和拾取其直达水波走时(图 3f).

2 视双测线OBS位置校正 2.1 校正原理及程序优化

本文采用的OBS位置校正方法是在前人所用方法(敖威等,2010张莉等,2013Du et al., 2018)的基础之上,经过改进形成的.相同之处:均是结合OBS记录的直达水波走时信息,采用最小二乘法和蒙脱卡洛法经过多次正反演迭代后,取直达水波走时残差均方根(Root Mean Square, RMS)最小值对应的随机点作为OBS的最佳落点,并同时求得最佳的海水声学速度和时间校正量.不同之处在于:(1)原有方法在反演海水声速时,将其初始阈值(郑红波等,2009)设为一个固定范围(通常根据物理海洋学所测得海水声学速度值设置为1480~1520 m·s-1),此时程序在手动(敖威等,2010张莉等,2013)或自动(Du et al., 2018)调整时间校正量时,就有可能造成反演声速值刚落入该阈值范围便停止反演,从而导致所获得的OBS落点位置、声速和时间校正量为局部最优的结果.为解决这一问题,我们参照地震数据处理中常用的“滑动窗口”思想,将海水声学速度初始阈值依据给定的步长划分为N等份,然后经过循环测试直至N个速度阈值均参与反演,最终选取其中RMS最小值所对应的结果作为最优结果,从而获得全局最优的OBS落点位置、声速和时间校正量(图 4图 5d).(2)原有方法主要应用于交叉测线或单条测线穿过的OBS,本次研究对单条测线通过的OBS位置校正加入视测线的直达水波走时,构成视双测线OBS位置校正,从而提高单测线穿过的OBS的位置精度.

图 4 OBS位置校正流程(据Du et al., 2018修改) Fig. 4 Flow chart of updated OBS relocation method (after Du et al., 2018)
图 5 OBS28位置校正结果精度对比 (a) OBS28台站相关测线位置.其中黑色圆圈为OBS28设计位置;(b) OBS28台站RMS分布.其中蓝色、红色、黑色曲线分别为单测线、交叉测线和视双测线校正结果;(c) OBS28台站位置误差椭圆分布.其中黑色圆圈为OBS28设计位置,蓝色、红色、黑色椭圆分别为单测线、交叉测线和视双测线校正结果;(d) OBS28台站声学速度反演结果.其中蓝色、红色、黑色圆圈分别为单测线、交叉测线和视双测线反演结果;实心点为最小RMS值对应的声速. Fig. 5 Precision contrasts of OBS28 relocation results (a) Locations of shooting lines and OBS28 station. Black circle represents designed position of OBS28; (b) Distributions of RMSs, OBS28. Blue, red and black lines are results of single line relocation, crossline relocation and quasi-crossline relocation, respectively; (c) Error ellipses of relocation results, OBS28. Black circle represents designed position of OBS28; blue, red and black ellipses are results of single line relocation, crossline relocation and quasi-crossline relocation, respectively; (d) Inversion results of water wave velocity, OBS28. Blue, red and black dots represent single line inverting method, crossline inverting method and quasi-crossline inverting method, respectively; solid dots are water wave velocities with minimum RMS.
2.2 可行性分析

以交叉测线穿过的OBS28台站为例,通过对比分析单测线、交叉测线和视双测线的OBS位置校正结果,探讨了视双测线OBS位置校正方法的可行性与精度.在本次对比试验中,单测线方法以L2测线为例,交叉测线方法以L2、H2测线为例,视双测线方法以L2、H3测线为例(图 5a).结果表明:使用三种方法校正的OBS28台站直达水波走时RMS均呈现以0 ms为中心的正态分布特征,视双测线校正RMS分布情况与交叉测线更加一致(图 5b);视双测线和单测线最终校正结果RMS值均与交叉测线相近(表 1),但视双测线校正所得的OBS位置更接近于交叉测线的校正结果,且误差椭圆长轴相较于单测线法有大幅缩小,在校正精度上也有明显的提高(图 5c);全局搜索中三种校正方法的RMS值均经历先变小后变大趋势,其中视双测线法与交叉测线法的变化趋势更接近,均表现为以最小值为中心的对称形态,反演声速结果均为1500.5 m·s-1,单测线法的变化趋势则表现为以最小值为中心的非对称形态,反演声速为1501.5 m·s-1,与前两者存在一定差异(图 5d).此外,对比校正前后OBS28台站主测线剖面的综合地震记录剖面发现,相对于位置校正前不对称的直达水波震相,经视双测线法校正后,直达水波震相呈现完全对称的“M”型(图 6),进一步表明该方法在位置校正中的准确性(敖威等,2010).因此,对于本次三维地震探测实验中仅有单测线穿过的OBS,采用视双测线位置校正方法显著地提高了校正精度.

表 1 OBS28台站位置校正结果对比 Table 1 Comparison of relocation results, OBS28
图 6 OBS 28台站沿L2测线的综合地震剖面 (a)位置校正前;(b)位置校正后;横轴为偏移距,纵轴为折合时间,折合速度8 km·s-1.带通滤波器f (Hz) = 1, 3, 8, 15. Fig. 6 Seismic record sections for OBS28 along the line L2 (a) Before OBS relocation; (b) After OBS relocation; the horizontal axis is shot number, the vertical axis is the reduced time, reduced velocity is 8 km·s-1. Band-pass filters, f (Hz) = 1, 3, 8, 15.
3 结果与讨论 3.1 OBS位置校正结果

本文对此次三维地震探测实验的49台OBS进行了位置校正(表 2图 7).除18台OBS应用交叉测线校正方法外,在分析其余30台单条测线穿过的OBS的地震数据时发现,其中28台OBS在临近测线上能够识别清晰可靠的直达水波震相,满足视双测线校正,另外3台OBS(OBS01、OBS02和OBS19)由于临近测线的距离过远,直达水波震相难以识别,故采用单测线法校正.

表 2 OBS位置校正结果 Table 2 OBS relocation result
图 7 南海洋陆转换带三维OBS探测实验的炮点与OBS位置校正结果 红色圆点为OBS投放位置;白色圆点为OBS校正后位置;黑色点为校正后炮点位置;黄色五角星为IODP 367/368航次钻孔位置. Fig. 7 Shots and OBS relocation results in the 3D OBS seismic experiment at COT zone in northern SCS The red and white circles indicate deployed and corrected OBSs′ locations, respectively; black dots are the relocated shots′ positions; yellow stars are sites of IODP Leg 367/368.

从校正结果可知:(1)所有OBS台站的RMS均呈现以0 ms为中心的正态分布特征,且RMS值均为4 ms左右,与设定误差扰动值相当,说明反演结果可靠,不存在系统误差.(2)视双测线和交叉测线的误差椭圆轴长约为35 m,而且长短轴比值接近1,说明OBS位置校正结果在各个方向的不确定性相当;单测线校正结果的误差椭圆长短轴比值约为20,表明其垂直测线方向的不确定性较大,这些结果与敖威等(2010)所取得认识一致.(3)反演所得研究区内水体最小声速为1494.0 m·s-1,最大声速为1503.1 m·s-1,平均声速约1499.3 m·s-1;与该区域表层(0 km)及深层(2.5 km)环流对比可见(Lan et al., 2013, 2015),局部水体结构在三维空间中变化较大,既有垂向分层特征,也存在横向的剧烈变化.(4)视双测线校正的大部分OBS台站在不同测线上二次时钟校正量(tadjust)基本一致,但也有部分台站存在几十毫秒的差距,可能与视测线上的直达水波信噪比低,且沿主测线方向产生震相叠加,最终造成走时拾取整体超前或滞后有关.(5)各OBS台站漂移距离(投放点与校正位置的距离)在75~541.7 m之间,平均值为217.8 m,整体偏移量远小于西南印度洋中脊三维OBS探测实验结果,与南海中央次海盆、巴士海峡等地三维OBS探测实验结果较为接近.

3.2 新位置校正方法的优势

与前人研究中的OBS位置校正方法相比,本文提出的视双测线OBS位置方法具有两方面的优势.(1)可以自动反演并获得全局最优解.在基于直达水波走时的OBS位置校正过程中,最初的方法是通过手动调整时间校正量来保证反演的海水声学速度落入其经验范围内(敖威等,2010张莉等,2013),耗时长、人为因素大,多次运行获得的位置校正结果可能有所差异.因此,Du等(2018)对程序进行了改进,使其可以自动调整时间校正量,减少了人为因素并节省了时间.然而,两者均未对海水声速经验区间进行全局搜索和分析,只是简单地将反演声速首次落入经验区间时的反演结果作为最优结果,从而造成两者得到的均为局部最优值.本次研究采用“滑动窗口”思想设计的全局校正方法,使得经验区间内所有取值均参与了反演,最后取所有窗口中RMS最小值所对应的结果作为全局最优值,从而解决了此前只能获得局部最优的问题.(2)提高了OBS位置校正结果的可靠性.以OBS28台站测试结果(表 1)为例,可以清楚地发现,虽然单测线法和视双测线法得到的RMS值均与交叉测线相差不大,但是视双测线法和交叉测线法的OBS定位结果更为接近,单测线法的OBS定位结果在垂直测线方向与前两者存在较大误差,而且其校正结果不确定度相比前两者存在显著差距(表 1).因此,该方法对于单测线OBS位置的校正精度有明显的改进,有利于提高后续三维速度结构模拟的精度.

3.3 L5测线震相初步识别

此处以NW-SE向横跨陆壳及洋壳的主测线L5为例,选取该测线上位于交叉测线位置的4台OBS(05,09,16和18),在OBS位置校正工作完成后,重新生成其SEGY格式数据,然后以8 km·s-1的折合速度生成各台站沿L5测线垂直分量的综合地震记录剖面(图 8),并对其震相进行了初步识别.由图 8可见,各台站所记录质量较好,信噪比高,震相清晰,震相延伸距离可达40~80 km,可识别出较为清晰的Pg、PmP和Pn震相(分别来自地壳的折射波、莫霍面的反射波和地幔的折射波);Pg震相起伏受多波束水深地形影响,OBS05台站中在Offset为-20 km处,对应着OBS05台站NW方向的地形隆起(图 8a图 7),在OBS09台站中也有显示(图 8b).Pg、PmP与Pn三个震相的交叉点位置,可以初步地指示地壳厚度,区分陆壳与洋壳.由4个台站的综合地震记录剖面来看,OBS05和OBS09两个台站(图 8a8b),交叉点位于Offset为40~50 km处,其下方的地壳厚度较厚,属于减薄的陆壳范围;而OBS16和OBS18两个台站(图 8c8d),三个震相的交叉点位于Offset为25~30 km处,其下方的地壳厚度较薄,属于正常的洋壳范围;因此,可初步推断南海北部洋陆转换带可能位于OBS09台站附近(18.4°N)(图 7),具体位置需要通过下一步深部速度结构的模拟结果进行确定.

图 8 4个OBS台站沿L5测线的垂直分量综合地震记录剖面 (a),(b),(c)和(d)分别为OBS05,OBS09,OBS16和OBS18台站沿L5测线垂直分量综合地震记录剖面图;横轴为偏移距,纵轴为折合走时,折合速度8 km·s-1.带通滤波器f (Hz) = 1, 3, 8, 15. Fig. 8 Vertical-component seismic record sections of 4 OBSs along the line L5 (a), (b), (c) and (d) represent seismic record section for OBS05, OBS09, OBS15 and OBS18 along the line L5, respectively; the horizontal axis is shot number, the vertical axis is the reduced time, the reduced velocity is 8 km·s-1. Band-pass filters, f (Hz)=1, 3, 8, 15.
4 结论

南海IODP367/368航次钻探区三维OBS深地震探测实验共投放52台OBS,最终成功回收49台,是一次非常成功的实验.本文在前人研究基础上,改进与应用新的方法,完成了此次实验炮点和OBS位置校正,获得以下几点认识和结论:

(1) 本研究利用连续三个炮点的平均航向将DGPS记录的炮点位置校正到了枪阵中心,从而更加真实地反映气枪枪阵与船体之间的软连接状态;根据“滑动窗口”思路改进了OBS位置校正程序,将原来的海水声学速度阈值划分为足够小间距的N等份,然后通过循环测试来获取全局最优的位置校正结果.

(2) 创新性地提出了视双测线OBS位置校正方法.位置校正结果表明,走时残差均方根(RMS)控制在5.8 ms以内,海水中声学速度反演结果稳定于1500 m·s-1;交叉测线法、视双测线法和单测线法的OBS位置校正误差椭圆长轴区间分别为19.6~91.5 m、24.4~35.7 m、441.5~489.4 m;采用视双测线新方法校正的28台OBS位置精度相较于先前的单测线法提升了一个数量级.

(3) 沿L5测线综合地震记录剖面数据质量好,信噪比高,震相清晰,记录了多组P波震相(Pg、PmP、Pn),震相延伸距离可达40~80 km,为下一步三维速度结构模拟奠定了坚实的数据基础.

致谢  感谢数据处理以及成图中使用的GMT绘图软件(Wessel and Smith, 1998).本研究的数据采集得到国家自然科学基金委员会共享航次计划(航次编号NORC2018-08)的资助,该航次由中国科学院南海海洋研究所“实验2”号科考船实施,在此一并致谢.
References
Ao W, Zhao M H, Qiu X L, et al. 2010. The correction of shot and OBS position in the 3D seismic experiment of the SW Indian Ocean Ridge. Chinese Journal of Geophysics (in Chinese), 53(12): 2982-2991. DOI:10.3969/j.issn.0001-5733.2010.12.022
Cameselle A L, Ranero C R, Franke D, et al. 2015. The continent-ocean transition on the northwestern South China Sea. Basin Research, 29(S1): 73-95. DOI:10.1111/bre.12137
Chen H, Qiu X L, He E Y, et al. 2019. Accurate measurement and inversion for the seafloor positions of Hadal landers. Chinese Journal of Geophysics (in Chinese), 62(5): 1744-1754. DOI:10.6038/cjg2019M0009
Creager K C, Dorman L M. 1982. Location of instruments on the seafloor by joint adjustment of instrument and ship positions. Journal of Geophysical Research, 87(B10): 8379-8388. DOI:10.1029/JB087iB10p08379
Ding W W, Li M B, Zhao L H, et al. 2011. Cenozoic tectono-sedimentary characteristics and extension model of the Northwest Sub-basin, South China Sea. Geoscience Frontiers, 2(4): 509-517. DOI:10.1016/j.gsf.2011.05.010
Du F, Zhang J Z, Yang F D, et al. 2018. Combination of least square and Monte Carlo methods for OBs relocation in 3D seismic survey near Bashi channel. Marine Geodesy, 41(5): 494-515. DOI:10.1080/01490419.2018.1479993
Hao T Y, You Q Y. 2011. Progress of homemade OBS and its application on ocean bottom structure survey. Chinese Journal of Geophysics (in Chinese), 54(12): 3352-3361. DOI:10.3969/j.issn.0001-5733.2011.12.033
Lan J, Zhang N N, Wang Y. 2013. On the dynamics of the South China Sea deep circulation. Journal of Geophysical Research:Oceans, 118(3): 1206-1210. DOI:10.1002/jgrc.20104
Lan J, Wang Y, Cui F J, et al. 2015. Seasonal variation in the South China Sea deep circulation. Journal of Geophysical Research:Oceans, 120(3): 1682-1690. DOI:10.1002/2014JC010413
Larsen H C, Mohn G, Nirrengarten M, et al. 2018. Rapid transition from continental breakup to igneous oceanic crust in the South China Sea. Nature Geoscience, 11(10): 782-789. DOI:10.1038/s41561-018-0198-1
Nakamura Y, Donoho P L, Roper P H, et al. 1987. Large-offset seismic surveying using ocean-bottom seismographs and air guns; instrumentation and field technique. Geophysics, 52(12): 1601-1611. DOI:10.1190/1.1442277
Stockwell J W, Cohen J K. 2002. The New SU User's Manual. Golden, CO: Center for Wave Phenomena.
Wang P X. 2009. Toward scientific breakthrough in the South China Sea. Journal of Tropical Oceanography (in Chinese), 28(3): 1-4.
Wang Y L, Yan P, Zheng H B, et al. 2007. Timing and positioning corrections of ocean bottom seismograph data. Journal of Tropical Oceanography (in Chinese), 26(5): 40-46.
Wessel P, Smith W H F. 1998. New, improved version of generic mapping tools released. Eos, Transactions American Geophysical Union, 79(47): 579. DOI:10.1029/98EO00426
Wilcock W. 2011. TlPicker Users' Manual. University of Washington. http://faculty.washington.edu/wilcock/?page_id=39.
William C T, Joseph E T. 1991. SAC-Seismic Analysis Code Users Manual. Livermore, CA: Lawrence Livermore National Laboratory, 1-153.
Xue B, Ruan A G, Li X Y, et al. 2008. The seismic data corrections of short period auto-floating ocean bottom seismometer. Journal of Marine Sciences (in Chinese), 26(2): 98-102.
Zhang J Z, Qiu X L, Zhao M H, et al. 2018. Abnormal data retrieval of three-dimensional OBS survey at the Bashi Channel area of the South China Sea. Chinese Journal of Geophysics (in Chinese), 61(4): 1529-1538. DOI:10.6038/cjg2018L0268
Zhang L, Zhao M H, Wang J, et al. 2013. Correction of OBS position and recent advances of 3D seismic exploration in the central sub-basin of South China Sea. Earth Science-Journal of China University of Geosciences (in Chinese), 38(1): 33-42. DOI:10.3799/dqkx.2013.004
Zhao M H, Qiu X L, Xia K Y, et al. 2004. Onshore-offshore seismic data processing and preliminary results in NE South China Sea. Journal of Tropical Oceanography (in Chinese), 23(1): 58-63.
Zheng H B, Yan P, Wang Y L, et al. 2009. Ocean thermohaline structure from processing multi-channel seismic reflection data of ocean water mass. Journal of Tropical Oceanography (in Chinese), 28(1): 48-55.
敖威, 赵明辉, 丘学林, 等. 2010. 西南印度洋中脊三维地震探测中炮点与海底地震仪的位置校正. 地球物理学报, 53(12): 2982-2991. DOI:10.3969/j.issn.0001-5733.2010.12.022
陈瀚, 丘学林, 贺恩远, 等. 2019. 深渊着陆器坐底位置的精确测量和反演计算. 地球物理学报, 62(5): 1744-1754. DOI:10.6038/cjg2019M0009
郝天珧, 游庆瑜. 2011. 国产海底地震仪研制现状及其在海底结构探测中的应用. 地球物理学报, 54(12): 3352-3361. DOI:10.3969/j.issn.0001-5733.2011.12.033
汪品先. 2009. 南海——我国深海研究的突破口. 热带海洋学报, 28(3): 1-4. DOI:10.3969/j.issn.1009-5470.2009.03.001
王彦林, 阎贫, 郑红波, 等. 2007. OBS记录的时间和定位误差校正. 热带海洋学报, 26(5): 40-46. DOI:10.3969/j.issn.1009-5470.2007.05.007
薛彬, 阮爱国, 李湘云, 等. 2008. SEDIS IV型短周期自浮式海底地震仪数据校正方法. 海洋学研究, 26(2): 98-102. DOI:10.3969/j.issn.1001-909X.2008.02.014
张佳政, 丘学林, 赵明辉, 等. 2018. 南海巴士海峡三维OBS探测的异常数据恢复. 地球物理学报, 61(4): 1529-1538. DOI:10.6038/cjg2018L0268
张莉, 赵明辉, 王建, 等. 2013. 南海中央次海盆OBS位置校正及三维地震探测新进展. 地球科学-中国地质大学学报, 38(1): 33-42.
赵明辉, 丘学林, 夏戡原, 等. 2004. 南海东北部海陆联测地震数据处理及初步结果. 热带海洋学报, 23(1): 58-63. DOI:10.3969/j.issn.1009-5470.2004.01.008
郑红波, 阎贫, 王彦林, 等. 2009. 研究海水温盐结构的反射地震方法. 热带海洋学报, 28(1): 48-55. DOI:10.3969/j.issn.1009-5470.2009.01.008