obs技术在南海的应用

obs技术在南海的应用

0陆相中脊设计范畴海底地震仪（obs）直接接触海底面。与多道地震的接收相比有很大优势。首先，可以观察到交叉波，通常使用四分量地震检测器。其次，它用于研究和获取深度地壳结构信息。第三，频率范围obs也可以进行海床的自然地震观测（浮秋宇等，2003a，2003b；陈刚等，2004）。obs广泛应用于石油工程调查（陈刚等，2009）、洋中脊柱结构（牛伟等，2009；陈刚等，2010；朱亚洲等，2012年；周志勇等，2009；陆学林等，2011年；赵明辉等，2011；吴振丽等，2011年；魏晓东等，2011年；黄海波等，2011年；陈刚等，2012年；陆川川等，2011）等。例如，在南海、黄海、西南部印度洋和南北极的地震波观测和勘探中，有许多实际的研究和研究工作，例如进口设备和中国、欧洲和美国、印度、南北极的研究。在这些研究中，不仅有进口设备，还有中国、福建省的、江西省的、江西省的、广东省的、，以及福建省的、江西省的和浙江省的。因此，在海洋操作、数据处理和研究方法上有许多复杂的问题。为了规范作业，提高回收率和数据处理效率，确保结果的可靠性和精度，有必要总结相关技术，为未来的更好地开展工作提供参考。本文以实际工作为例（西南印度洋中脊人工起源2d段），系统地总结了obs地震观测的海上操作方法、数据处理技术和流程、建模正反演和模型评价方法。1obs投放和开放系统正常设置OBS的海上作业主要分5步:(1)测线设计和震源设计;(2)OBS参数设置;(3)OBS投放;(4)炸测作业和(5)OBS回收.OBS记录一般要求同时获得海底反射、沉积层折射、地壳折射、莫霍面反射和上地幔折射等地震信号,因而要求炸测时间间隔较大,以防止前炮的后续震相被后炮的浅层震相所覆盖.实践中常采用炮间距150~300m,航速4~5kn,相应的放炮时间间隔为80~120s.OBS站位间距尽可能小些,一般不超过30km.为了得到深部信息,要求枪阵容量大于5000立方英寸,工作压力120kg/cmOBS投放前要进行充电、参数设置、密封和上浮系统测试.实践中设数据采样频率为100~500Hz,其它参数包括开始记录时间、终止记录时间和自动上浮时间.自动上浮设置是在正常的声学释放系统由于各种原因失效时,采用的补救手段.另外,需要配置罗经用于确定OBS水平分量在海底的方位,配置闪光灯和旗子,以便在夜间和阳光强烈的条件下寻找在海面上漂浮的OBS.OBS投放是根据事先设定的点位,利用GPS定位来进行,要求提前降低船速,投放时船速不得高于2kn,使OBS以较小的水平速度平稳入水.同时要及时记录实际投放点的坐标(由于船在运动,投放点与设计点之间总是有偏差),其重要性在于实际投放点对后面的回收意义重大,而且常作为数据处理的初始坐标.OBS回收作业,正常情况下船开到投放点上风向1km处,关闭船舶的螺旋桨和主机(防止螺旋桨搅住释放器探头和噪声对声学释放的干扰),把声学释放器探头放入水中(保证深过船底),发送声学释放命令.可以灵活设计回收路线,尽可能使每个站位有2次以上的回收机会.发现目标和打捞上船时要分别记录坐标,这有利于判断海流方向为后续工作中的目标搜索提供参考.2其它常用数据格式OBS数据处理主要包括炮点位置和时间校正、OBS位置校正、OBS数据格式转换、OBS时钟漂移校正、增益恢复、滤波及预测反褶积(李湘云等,2007;王彦林等,2007;Oshidaetal.,2008;薛彬等,2008;夏常亮,2009;敖威等,2010;夏少红等,2011),为了更准确的识别震相,有时还需要进行水深静校正(Xianetal.,1999).常用数据格式采用SEGY及SU格式文件,因为可处理这类格式的商用软件较多,方便使用.下面具体阐述(以西南印度洋中脊50°E附近OBS地震试验(LiandChen,2010;Zhaoetal.,2013)的某测线为例).2.1船炮点和炮点位置校正用GPS对放炮时间进行控制,使用精确到毫秒的精密计时器,与震源直接相联,如果有困难,可以将多余的OBS或陆地地震仪(带GPS接收器)安置在船上,记录放炮产生的甲板振动脉冲.由于各种原因,航迹(炮点位置)相对设计测线总会产生一定的偏差(一般为几米到十几米),直接影响震相走时的拾取.校正方法是将船上GPS所记录的位置校正到实际的枪阵震源中心位置(熬威等,2010):2.2位置校正前的地震记录自由落体投放的OBS下落速度较小(约1m/s),受海流的影响,其在海底的实际落点位置会偏离设计点,这种偏离现象对OBS记录剖面后中偏移距较小的震相(如直达水波)的拾取有非常大的影响(如图1a为未做任何校正的地震记录),当水深较大时偏移可能更大.因此对OBS的位置进行校正是必要的(如图1c,位置校正前的地震记录;图1d,位置校正后的地震记录).3D地震试验与2D地震试验采用不同的校正方法.2.2.1实验1:青少年暴力型当有交叉放炮测线经过OBS上方时,通常可以得到较准确的OBS位置(夏常亮,2009;敖威等;2010;夏少红等,2011).方法是用读取所有站位OBS数据垂直分量中与该站位位置最近的5炮(或更多)的初至波走时,精确到毫秒,并进行时间漂移和滤波延迟校正(后述),然后减去校正后的放炮时间,差值即为直达水波的单程走时.迭代方法方程为2.2.2obs位置的确定与3D不同,对于2D试验,首先用最小二乘法将所有炮点归一到同一直线上形成剖面.OBS初始平面坐标由投放点的坐标确定,深度由水深测量确定.然后将其投影到上述炮点测线,得到局部化坐标以及偏离实际点的误差(薛彬等,2008).在此基础上计算小偏移距附近的直达水波的理论走时曲线并与实测记录相拟合,对OBS位置进行人工微调.2.3obs时间漂移校正OBS投放前其记录器参数设置使用GPS授时,回收后立即再次用GPS对时,确定总的时钟漂移量(薛彬等,2008).一般认为时间漂移是线性的,因而可以由工作起始时间、终止时间和总漂移量,计算出OBS每个记录道(对应一个炮点)的时间漂移.实际工作表明时间漂移的校正量约为几个到十几个毫秒.然后对地震记录剖面进行校正,这个过程中可再次对OBS位置进行微调.经过上述OBS位置和时间的校正,小偏移距附近海底面反射震相的形态得到了较好改善(图2d).2.4增益恢复、滤波及预测反褶积对所有站位的OBS地震记录,可使用软件包SeismicUnix(SU)(CohenandStockwell,1995)进行必要的处理.如最小相位带通滤波(图1b),增益恢复处理(图1c)等.滤波要根据气枪主频进行(赵明辉等,2008),可以使用F-K滤波消除紧随直达波的干扰波,噪声干扰严重时可采用一致性滤波(MilkereitandSpencer,1989),还可采用预测反褶积压制多次波(图1e)便于识别沉积层和上下地壳分界面的反射震相.2.5深水静校正处理OBS地震记录是以共接收点记录的折合时间剖面的形式来表达,常用的折合速度为6.0km/s或8.0km/s,前者用于突出地壳内折射波Pg,后者用于突出上地幔顶层内折射波Pn.但当地形变化较大时,对震相识别造成一定影响.为此,可以采用水深静校正处理,从每道记录的时间轴上减去水层内垂直走时.图1f给出了水深静校正之后的地震记录,该OBS位于地形隆起处,海底面起伏高差达2000m,可以看出校正后Pg的视速度接近真实,利于震相识别.通常水深静校正用于震相识别,震相拾取仍使用先前的地震记录.3初始模型与速度参考研究区历史地震剖面、水深数据等资料,用网格法建立初始模型,用试错法正演手动修改各层的速度和分界面变化等模型参数.固定界面,以全局走时误差为目标函数,采用阻尼最小二乘法,由浅到深,逐层反演速度,得到最优化速度模型.最终模型可分为最佳网格参数和最简洁结构模型、最少参数或结构优先模型、最佳网格参数和结构优先模型以及最少参数或最简洁结构模型(ZeltandSmith,1992;ZeltandForsyth,1994;Zelt,1999)等4种.举例如下.测线初始模型分为4层(图2a),分别为海水层(层内速度为1.5km/s),2km厚的上地壳(洋壳层2,顶面速度2.4km/s,底面速度6.4km/s),4km厚的下地壳(洋壳层3,顶面速度6.4km/s,底面速度7.0km/s)和上地幔(顶面速度7.8km/s,底面速度8.2km/s).上述分层和速度参考了多波束水深数据、标准洋壳模型(Kennett,1982;Whiteetal.,1992)和SWIR66°E及57°E处的P波速度模型(MinshullandWhite,1996;Mulleretal.,1997,1999,2000;Minshulletal.,2006).网格划分方法是,各层网格节点水平间距为5km.垂直方向上,除海水层外,设定层内速度为线性变化,由顶面与底面的速度差和厚度,自动确定网格划分.层内不设速度间断面.采用走时模拟和反演方法来构建速度模型,使用软件为RayInvr(ZeltandSmith,1992).在初始模型基础上,通过上述反演方法得到最终速度模型(图2b).4模型评价由于模型参数存在非唯一性和不确定性,模型评价主要分为两种:模型间接评价和模型直接评价(Zelt,1999).4.1速度模型的计算模型间接评价主要评价模型的分辨率和模型参数的不确定性,是最简单的模型评价方法,用统计图表来表达.(1)用表格给出模型统计分析结果(表1),包括震相拾取数、模型可以约束的震相数及其百分比、每种震相的RMS(rootmeansquare)走时拟合误差和chi-square、所有走时的RMS与chi-square(ZeltandSmith,1992).(2)画出每一条射线和每一震相的计算走时图,同时将其与拾取的震相对比(图3).(3)画出最终速度模型图的同时计算并给出射线密度分布图(图4).(4)用扰动方法进行模型分辨率测试(图5).这个方法比射线密度图更好,因为同一角度或同一方向的射线密度大并不能代表较高的分辨率,不同方向和偏移距的较少的射线也能有较高的分辨率.另外常用的技术还有检测板分析(Xiaetal.,2010)等.4.2算法2:射线追踪拟合,多参数不确定性测试,通用模型评价直接模型评价主要是通过改变模型参数来检验模型的唯一性(Zelt,1999),其操作较为复杂和耗时.主要有三种方法:(1)单一参数不确定性测试:手动修改某一参数,然后进行射线追踪拟合,使用F-test来检验模型是否相同,如果模型不同,当前改变量则为该模型参数的不确定性.(2)多参数不确定性测试:手动修改多个参数,然后进行射线追踪拟合,原理同单参数不确定性测试,本文使用测线的多参数不确定性测试见表2.(3)通用模型评价:在合理的范围内系统的改变任一模型参数,再做反演,判断对模型的影响.如改变震相的不确定性,速度和深度节点,目标函数或线性系统的射线出射偏移距参数,或模型网格节点或层数等,可以判断模型的唯一性和检验初始模型对最终模型的影响.5深水静校正处理本文根据多年实际经验,以西南印度洋中脊一条实测OBS广角地震剖面为例,系统地介绍和总结了OBS数据采集处理建模反演和模型评价的流程和方法,为以后更好地开展相关工作提供参考.主要结论如下:(1)海上作业前要精心设计测线和工作程序,重点是震源参数和记录参数.OBS回收成功率是关键问题,工作要到位细致,特别要在投放前进行必要的测试.(2)OBS数据处理主要包括炮点的时间和位置校正、OBS的位置校正、OBS数据格式转换、OBS时钟漂移校正、增益恢复、滤波及预测反褶积处理,如果地形起伏大,还需进行水深静校正.(3)震相识别是建模和反演的又一关键.除了从折合速度剖面上直接人工识别外,要结合正演模拟反过来验证与拾取.(4)参考研究区历史资料,特别是地震和水深资料建立初始模型,用网格化方法对初始模型参数化,用射线追踪正演和试错法确定速度模型和地壳结构的基本轮廓,然后固定界面,对速度分布逐层进行自动反演.如是3D资料,最好先进行2D建模,然后再进行3D层析成像自动反演.(5)用多种手段进行模型评价是非常必要的,得到模型的不确定性和分辨率,保证可靠性和唯一性.牛雄伟,阮爱国,吴振利,等.2014.海底地震仪实用技术探讨.地球物理学进展,29(3):1418-1425,doi:10.6038/pg20140358.NIUXiong-wei,RUANAi-guo,WUZhen-li,etal.2014.ProgressonpracticalskillsofOceanBott