SEDIS IV型短周期自浮式海底地震仪数据校正方法

利用15台SEDISIV型短周期自浮式海底地震仪在南海中、北部地壳深部结构调查中所获得的资料,探讨了海底地震仪数据校正的方法和校正后的效果,结果表明:使用该地震仪所获得的原始资料经过放炮时间、炮点坐标数据局部化、海底地震仪位置误差以及记录时间漂移4方面的校正后,数据更趋合理,误差显著降低。放炮时间的校正消除了时钟漂移和时间延迟的误差;炮点坐标数据局部化处理消除了炮点位置整体趋势性偏移的现象;试错法进行位置误差和记录时间的精细校正时,时间漂移的校正量值约为几个到十几个毫秒,位置校正的量值仅在几米到数百米之间,实测数据所绘曲线的形态和位置都与理论曲线十分吻合,可见校正后误差显著降低。     

海底地震仪的时间同步技术

海底地震仪(OBS)是记录海底地震数据的主要仪器.在我国,OBS仍处于研究与实验阶段.时间同步的精准度是海底地震仪的重要指标,直接影响了地震数据分析和地震数据反演的准确性.以法国MicroOBS _Plus为例,对海底地震仪时间同步技术进行深入分析,详细介绍其时间同步原理及实现方法,为OBS在海底地震地壳深部结构调查研究等方面的应用提供参考依据.     

海底地震仪重新定位:以东沙群岛海域为例

利用海底地震仪(OBS)进行探测时,一般采用自由下落的方式将OBS布设在海底。由于海流和海底地形的影响,OBS的位置一般都偏离设定的位置。OBS重新定位是OBS数据处理的基础,不正确的位置信息将导致错误的观测系统信息,从而影响后续的处理效果。直达水波包含了OBS的位置信息,一般利用直达水波的走时,采用最小二乘的方式来确定OBS的实际位置。炮点位置精度、放炮时间延迟、炮点分布方式、OBS时钟漂移、走时拾取误差以及海水速度变化等因素都会影响重新定位的精度。文章采用数值方式研究这些因素对重新定位精度的影响,并对东沙群岛海域实际的OBS站位进行精细重新定位。数值结果表明,直线排列的炮点不能很好地反演出OBS的位置,应该采用十字或者井字分布的炮点;在偏移距小于10km时,海水速度结构的精度对重新定位的影响可以忽略;对OBS重新定位影响最大的是由放炮时间延迟、时钟漂移以及走时误差共同构成的右端项误差,5ms的右端项误差可引起40m左右的偏差。实际数据重新定位的结果较好,符合数值研究的结果。     

海底地震仪在天然气水合物勘探中的应用综述

海底地震仪(OBS)直接布设在海底,在地球动力学研究、油气勘探以及地震学研究等众多领域中应用广泛。随着采集技术的提高,海底地震仪也越来越广泛地应用到水合物的勘探中。文中从OBS数据采集、数据处理以及速度反演三个方面对OBS在水合物勘探中的应用现状进行综述,重点在于综述目前的应用现状和OBS使用的难点。数据采集方面主要综述了在OBS水合物勘探中几种常见的观测系统;OBS数据处理方面主要综述OBS数据预处理和OBS数据镜像成像;OBS速度反演方面主要综述常用的反演方法及其优缺点;最后,在以下几个方面对OBS在水合物勘探中的应用进行了展望:(1)优化观测系统的设计;(2)发展更先进的波场分离算法;(3)基于全波形反演以及基于面波的速度反演。     

深水环境下OBS的二次定位技术

在海洋环境下,OBS（海底地震仪）投放后,由于潮汐、涌浪以及海流等因素的影响而漂移,OBS资料处理的首要任务就是进行OBS的二次定位。深水环境下一般使用直达波二次定位技术,直达波定位的传统方法是采用网格节点法计算直达波线性动校正,计算工作量非常大。列举了OBS落点漂移的多种情形,模拟了不同情况下直达波的线性动校正曲线,分析了影响线性动校正曲线变化的各个因素,建立了判断落点方向的标准,提出了十字交叉法反演OBS落点位置,对实测资料进行落点位置反演,获得了准确的定位结果。此方法明显减少了计算工作量,可以快速获得准确的OBS落点位置。     

利用海底地震仪数据分析台风对海底环境噪音的影响

在海底布设的海底地震仪（OBS）能比较清晰地记录到海底的环境噪音，而台风可以直接或间接的产生在海底传播的弹性波，从而影响海底的环境噪音，并在较大程度上影响OBS的数据记录。本文通过分析台风对工作区的整个影响过程中OBS记录数据的振幅变化，再选择合适的滤波方式，首次发现台风产生的风浪及涌浪在短周期海底地震仪的记录数据上有良好的表现特征，指出了台风对海底环境噪音的另一种可能的影响方式，并由此得出：1）台风产生的风浪和涌浪对海底环境噪音的影响模式不同；2）风浪和涌浪所加强的海底环境噪音的范围和程度不同；3）短周期OBS可以比较清晰的记录涌浪信息，其周期主要是6—8 s，且能量稳定（简称“8秒现象”）。这三点结论为后期的海洋地震研究和海洋学其他研究提供经验与借鉴。     

海底地形校正在正确拾取OBS震相中的重要作用

海底地震仪(ocean bottom seismometer,OBS)探测是获得海底深部地壳结构的首选方法。正确拾取OBS记录中折射/反射震相,对获得准确的深部速度结构非常重要。当OBS布设在崎岖海底时,起伏的海底地形会影响OBS地震记录剖面中Pg、Pm P、Pn等岩石圈内部震相的展布特征,如南海东部次海盆中横穿珍贝-黄岩海山链的A4M4地震测线和西南印度洋中脊横穿扩张脊的Y3Y4测线,强烈的地形高差变化增加了这些测线上OBS台站地震剖面的震相识别难度。在震相拾取之前,通过地形校正方法消除海底地形对震相的影响;地形校正后,根据震相视速度及其展布趋势可以准确地识别震相,有效地提高了震相识别的可靠性。地形校正方法在上述地区的震相识别与速度结构研究方面发挥了重要作用,同时为今后在其他地形变化复杂地区的OBS震相识别提供了经验与借鉴。     

海底地震仪实测信号特性分析

置于海底数百米至数千米的海底地震仪(OBS)的实测信号相比陆地地震仪具有不同的特性;由于水的作用或记录信号源频率的不同,短周期OBS和宽频带OBS记录的信号又有明显的差异。文章对南海西南次海盆地震探测期间记录的人工气枪震源和天然地震实测信号进行了时频分析,结果如下。1)气枪作业后在海底激发两种噪声:一是水的波动不断叠加形成的长波,周期50s左右,以水平分量为主;二是高频噪声,主要是OBS底座细微晃动引起的。2)宽频OBS对于水下移动目标激发海底波动具有很好的探测能力,特别是水平分量可以获得大振幅且周期特征清晰的记录,并能够指示方向。3)宽频OBS能记录到清晰的天然地震信号,为研究调查区岩石圈结构增添了更多的信息,短周期OBS对远震直达P波有很好的记录。国产宽频I-4C型OBS碰巧记录了日本M9.0级大震。     

国产OBS的时间矫正方法——以2019年台湾海峡地壳结构海陆探测实验为例

时间服务系统对利用走时层析成像方法进行地下介质速度结构反演至关重要。海底地震仪(ocean bottom seismometer, OBS)工作期间由于没有GPS时间接入, 其时间误差(包括守时误差和授时误差)主要来源于内部石英晶振的准确程度, 受到外部环境变化以及开关机等因素影响。长期实践发现, 部分国产OBS在记录气枪信号以及天然地震信号时存在较大的时间偏差。本文对2019年福建及台湾海峡地壳结构海陆探测实验所获得的53台次国产OBS记录进行了时间服务系统矫正。其中, 针对OBS授时误差, 利用出海前不中断采集的一致性试验和运输船运输过程中产生的晃动互相关进行时间矫正; 针对守时误差, 采用计算实际采样频率与理论采样频率偏差进行矫正; 通过对比矫正前后OBS记录到的天然地震信号, 进行秒级别的检测。结果表明, 经过以上步骤矫正的OBS数据, 其时间记录的准确性得到了显著提高, 从而降低了震相识别、走时拾取的时间误差, 为标准化国产OBS数据采集作业流程提供了重要参考。     

南沙地块海底地震仪转换横波震相识别最新进展

通过挖掘海底地震仪记录的转换横波信息,可以促进对地壳岩性、均质性的认识。以南海南沙地块OBS973-1剖面上18个台站的数据为实例,阐明了海底地震仪横波震相识别的方法。首先对三分量OBS(Ocean Bottom Seismometer)地震数据进行带通滤波、维纳滤波、极化滤波等去噪处理,然后利用能量扫描法求得极化角进行水平分量坐标旋转,求取最佳径向分量。数据处理结果表明,相对于OBS记录的磁罗盘方位角,能量扫描法求取的极化角更为准确可靠。最后通过OBS973-1地震剖面垂直分量与径向分量上的纵横波走时对比、质点运动轨迹、速度模型试算等手段,进一步确定了转换横波震相的类别。在南沙地块OBS探测中成功地在10个台站中识别出了PgSs、PnSc、Pms等震相,不仅可以为下一步横波速度结构模拟提供坚实的数据基础,而且可以为今后OBS转换横波在其他地区的有效应用与推广提供借鉴。     

基于背景噪声互相关技术的海底地震仪的时钟校正方法在西南印度洋脊的应用

Seismic monitoring using ocean bottom seismometers(OBS) is an efficient method for investigating earthquakes in mid-ocean ridge far away from land. Clock synchronization among the OBSs is difficult without direct communication because electromagnetic signals cannot propagate efficiently in water. Time correction can be estimated through global positioning system(GPS) synchronization if clock drift is linear before and after the deployment. However, some OBSs in the experiments at the southwest Indian ridge(SWIR) on the Chinese DY125-34 cruise had not been re-synchronized from GPS after recovery. So we attempted to estimate clock drift between each station pairs using time symmetry analysis(TSA) based on ambient noise cross-correlation. We tested the feasibility of the TSA method by analyzing daily noise cross-correlation functions(NCFs) that extract from the data of another OBS experiment on the Chinese DY125-40 cruise with known clock drift and the same deployment site. The results suggest that the NCFs' travel time of surface wave between any two stations are symmetrical and have an opposite growing direction with the date. The influence of different band-pass filters,different components and different normalized methods was discussed. The TSA method appeared to be optimal for the hydrophone data within the period band of 2–5 s in dozens of km-scale interstation distances. A significant clock drift of ~2 s was estimated between OBSs sets through linear regression during a 108-d deployment on the Chinese cruise DY125-34. Time correction of the OBS by the ambient noise cross-correlation was demonstrated as a practical approach with the appropriate parameters in case of no GPS re-synchronization.     

Combination of Least Square and Monte Carlo Methods for OBS Relocation in 3D Seismic Survey Near Bashi Channel

Abstract

The relocation of ocean bottom seismometers (OBSs) is a key step in analyzing the three-dimensional seismic tomographic structure of crust and mantle. In order to get the accurate location of OBSs on the seafloor, we analyze the travel times of direct water waves emitted by air-guns. The Monte Carlo and least square methods have been adopted to calculate the true OBS location. The secondary time correction is necessary if the arrivals of direct water waves show overall time drift during relocation which maybe originates from remnant of linear clock drift correction and average errors of travel time picking, mean water velocity assumption, and experiment geometry. We have improved the original OBS relocation procedure which we used previously for other experiments by deliberateness of a secondary time correction and automatically approaching the really mean water velocity. A series of synthetic tests are carried out firstly to document the feasibility of our procedure and then it is applied on a real experiment. In here, we relocate 28 OBSs in total were relocated in 3D seismic survey near Bashi Channel. Relocation results show that the drifting distances for the 28 OBSs range from 65 to 1136 m between the deployed and relocated locations deduced by relocation results. The Pearson correlation coefficient between OBS drifting direction and sea current direction is 0.79, indicating that the two sets of data are highly linearly related and further manifest the sea current as the most possible driving force for OBS drifting during landing on the seafloor but its detailed influence mechanism is unclear by now. This research is necessary and critical for velocity structure modeling, and the optimal relocation program provides valuable experiences for 3D seismic survey in other area.     

台湾海峡OBS地震记录中二次Ps震相的特征及应用

在海底地震仪(ocean bottom seismometer, OBS)广角地震记录剖面上, 经常可以见到震相清晰且连续的多次波信号, 多次波和初至波是由相同的震源信号产生的, 也是地壳真实结构的反映。但是在通常的OBS数据处理过程中, 经常将多次波作为无效信号剔除掉, 对其属性及应用的研究比较少。文章通过对台湾海峡南部OBS探测测线HXN01数据的处理, 对多个台站记录到的二次震相进行了识别与拾取, 并以OBS0106台站为例, 对识别出的二次Ps震相进行了系统的研究分析, 发现二次Ps震相的波形特征和质点运动轨迹与初至震相相似, 但波形最大振幅值明显大于初至震相。通过Rayinvr射线追踪方法模拟, 确定了二次Ps震相的主要反射层, 并发现加入二次震相后, 台站下方浅部沉积层射线覆盖密度有显著提升, 射线覆盖的区域也明显增加, 为沉积层精细结构的反演提供了更为丰富的数据基础。另外, 对理论模型的地壳结构进行加入二次Ps震相前后的反演测试, 结果显示加入二次Ps震相数据后, 沉积层的界面深度误差得到明显的改善。     

OBS海洋环境信号分析与应用

为研究海洋环境信号在OBS(Ocean Bottom Seismograph)原始数据中的规律及应用,根据OBS原始数据的波形及频谱特征,将研究区划分为5个时间段,依次为旧涌浪阶段、风浪渐强阶段、风浪全盛阶段、风浪消退阶段和新涌浪阶段。结合海洋天气预报,认为上述现象是由偏南风风浪对海流的影响造成的。参考野外地震数据采集记录班报,得到各阶段的时长和距离,计算风浪渐强、全盛和消退阶段OBS附近海流的平均速度。结果表明:OBS原始资料中浅海海洋环境噪音增强的主要因素是风浪,且风浪引起的噪音信号的波形变化特征是渐进式的;OBS可用于接收某种特殊阶段(如台风、海啸等)的噪音信号,并根据噪音信号的波形特征、频谱变化规律和持续时间估算该阶段的海流速度变化。     

High resolution velocity analysis of ocean bottom seismometer data by theτ-p method

A method of high resolution seismic velocity analysis for ocean bottom seismometer (OBS) records is applied to the study of the shallow oceanic crust, especially sedimentary and basement layers. This method is based on the direct-p mapping and the-sum inversion. We use data obtained from a 1989 airgun-OBS experiment in the northern Yamato Basin, Japan Sea and derive P- and S-wave velocity functions that can be compared with the seismic reflection profiles. Using split-spread profile records, we obtain interface dips and true interval velocities from the OBS data. These results show good agreement with the reflection profile records, the acoustic velocities of core samples, and sonic log profiles. We also present a method for estimating errors in the derived velocity functions by calculating covariance of the derived layers' thicknesses. The estimated depth errors are about 150 m at shallow depths, which is close to the seismic wavelength used. The high resolution of this method relies on accurate determination of shot positions by GPS, spatially dense seismic observations, and the use of unsaturated reflected waves arriving after the direct water wave that are observed on low-gain component records.     

Velocity structure and gas hydrate saturation estimation on active margin off SW Taiwan inferred from seismic tomography

This paper presents results of a seismic tomography experiment carried out on the accretionary margin off southwest Taiwan. In the experiment, a seismic air gun survey was recorded on an array of 30 ocean bottom seismometers (OBS) deployed in the study area. The locations of the OBSs were determined to high accuracy by an inversion based on the shot traveltimes. A three-dimensional tomographic inversion was then carried out to determine the velocity structure for the survey area. The inversion indicates a relatively high P wave velocity (Vp) beneath topographic ridges which represent a series of thrust-cored anticlines develop in the accretionary wedge. The bottom-simulating reflectors (BSR) closely follow the seafloor and lies at 325 ± 25 m within the well-constrained region. Mean velocities range from ~1.55 km/s at the seabed to ~1.95 km/s at the BSR. We model Vp using an equation based on a modification of Wood’s equation to estimate the gas hydrate saturation. The hydrate saturation varies from 5% at the top ~200 m below the seafloor to 25% of pore space close to the BSR in the survey area.     

Velocity-Interface Structure of the Southwestern Ryukyu Subduction Zone from EW9509-1 OBS/MCS Data

A wide-angle seismic survey, combining ocean-bottom seismometers (OBS) and multi-channel seismic (MCS) profiling, was implemented in the southwestern Ryukyu subduction zone during August and September 1995. In this paper, we present the data analysis of eight OBSs and the corresponding MCS line along profile EW9509-1 from this experiment. Seismic data modeling includes identification of refracted and reflected arrivals, initial model building from velocity analysis of the MCS data, and simultaneous and layer-stripping inversions of the OBS and MCS arrivals. The velocity-interface structure constructed along profile EW9509-1 shows that the northward subduction of the Philippine Sea Plate has resulted in a northward thickening of the sediments of the Ryukyu Trench and the Yaeyama accretionary wedge north of the trench. The boundary between the subducting oceanic crust and the overriding continental crust (represented by a velocity contour of 6.75 km/s) and a sudden increase of the subducting angle (from 5 degrees to 25 degrees) are well imaged below the Nanao Basin. Furthermore, velocity undulation and interface variation are found within the upper crust of the Ryukyu Arc. Therefore, the strongest compression due to subduction and a break-off of the slab may have occurred and induced the high seismicity in the forearc region. This revised version was published online in November 2006 with corrections to the Cover Date.