海上地球物理浅地层剖面测量环境噪声分析.docxVIP

海上地球物理浅地层剖面测量环境噪声分析 平层勘探技术起源于19世纪60年代，是海上地球物理调查的重要手段之一。该探测技术基于水声学原理, 与反射地震学相似, 通过发射换能器将控制信号转换为不同频率声波信号向水下发射, 声波在传播过程中遇到地层分层界面以及地层沉积面等声阻抗界面后, 经反射返回接收换能器, 由于反射回来的信号带有地层信息, 所以利用该原理可以获得反映浅地层声学特征的剖面。目前, 该技术广泛应用于天然气水合物探测、近岸海洋地质环境调查、工程地质调查、区域地质调查、大洋科考等方面, 可获得海底浅部地层结构和构造, 为基础地质调查和矿产资源开发提供科学依据。 浅地层剖面仪是在测深仪的基础上发展起来的, 其设备经历了从固定频率和线性调频(chirp技术) 到声学参量阵式的发展。浅地层剖面测量在实际作业中容易受到周围环境的干扰, 如涌浪、机械振动、螺旋桨转动以及电力干扰等次生干扰源的影响, 很大程度上降低了采集资料的信噪比和分辨率, 进而影响对地质结构的解释, 因此, 必须对原始资料进行有效处理才能获得高质量的数据剖面。 对于浅地层剖面资料的处理主要是参照常规地震资料的处理方法, 而单独针对浅地层剖面数据处理方法的研究不多, Quinn等在已知震源特征前提下采用相关、反褶积以及预测滤波等方法对chirp型浅剖数据处理方法进行了系统研究, 吕国涛等利用图像处理方法对浅剖图像进行了滤波处理。本文结合前人的研究成果综合总结阐述了浅地层剖面测量从资料采集到后处理等阶段用来消除噪音和提高分辨率的数据处理方法。 1 劳动原则和设备介绍 1.1 反射波的强度 通常把海底地质结构看作一个层状的模型 (图1) , 不同介质层中密度、压强等不同导致声波在各层介质中传播速度就会有差异。近似声波向下垂直入射, 海水层作为第1种介质, 它的密度为ρ1, 声波在其中传播的速度为v1, 假设海底下各层介质密度和速度分别是ρ1、v1、ρ2、v2、…ρn、vn。当声波传播到两种不同介质分界面时, 一部分声波在分界面处发生反射, 另一部分沿法线方向继续传播, 在下一分界面处继续发生反射和透射。衡量各界面的反射的强度的参数称为反射系数, 其数学表达式如下: 式中:R为反射系数; ρ1、v1为第1种介质的密度和声波在其中的传播速度; ρ2、v2为第2种介质的密度和声波在其中的传播速度。 可以看出, 获得反射声波的强弱取决于界面两侧物质的差异。当ρ2v2-ρ1v1的值较大, 即界面两侧物质差异较大时, 接收到的反射信号就比较强, 反之则比较弱。由于接收到的反射信号携带了水下地层分界面的大量有用的地质信息, 所以通过观测记录并分析海底沉积物对于声波的反射特性, 就可以了解沉积物的地质属性, 直观地描述地质构造。 1.2 海底浅部区域波场观测系统 本文中信号采集设备为ATLAS PARA-SOUND 70 (以下简称为P70) 型深水参量浅地层剖面仪, 其原理是利用声波在介质中传播的非线性特性, 即2个沿同一方向传播的高频初始声波在远场介质中可以获得差频的声波信号。具体如下: 参量阵换能器在高压下同时向水底发射2个频率接近的高频声波主频信号 (f1, f2) , 在换能器以下会产生一系列二次频率声波信号, 如f1, f2, (f1+f2) , (f1-f2) , 2f1, 2f2, 声波信号改变2个主频 (f1, f2) 的频率就可以控制差频波 (f1-f2) 的频率, 当f1, f2的频率非常接近时, 差频 (f1-f2) 的频率相对低, 但是其依然保持高频时的波束角不变, 因此, 差频信号具有很强的沉积层穿透能力以及横向分辨能力, 可以用来精细探测海底浅部地层结构, 同时, 反射的主频声波 (f1, f2) 信号还可以用于精确的水深测量。 P70浅地层剖面仪 (见图2) 是目前国际上最先进的、最灵活的、功能强大的浅地层剖面仪之一, 换能器固定安装于船底, 不仅可以发射连续波 (continuous wave) 脉冲, 也可以发射脉冲调制的波形 (chirped pulse和barker coded pulse) , 该系统主工作频率在18~40kHz之间, 二次频率在0.5~6kHz之间, 工作水深10~11 000m, 其具有良好的地层穿透能力和地层分辨能力, 最大穿透深度可达200m, 垂向地层分辨率达到15cm, 其换能器共有256个发射通道, 每个通道使用98Khz的初始采样率对测量数据进行采样并且都具有独立的24位A/D转换;该设备还具有发射角自动控制功能, 因此, 可以通过计算海底坡度, 使窄波束很精确的入射到海底达到最佳剖面效果, 增加回波强度, 特别当海底地形变化复杂时, P70对海底跟踪能力很强。 2 噪比的基础 信噪比和分辨率是衡量野外采集数