1、第 卷 第期 年月 ，中文引用格式：张体强近岸海域浅地层剖面数据处理技术及应用工程地球物理学报，（）：英文引用格式：，（）：近岸海域浅地层剖面数据处理技术及应用张体强（上海东海海洋工程勘察设计研究院有限公司，上海 ）收稿日期：基金项目：国家海洋局东海分局青年科技基金（编号：）作者简介：张体强（），男，工程师，主要从事海洋测绘和海洋工程勘察工作。：摘要：浅地层剖面探测技术是应用声学原理高效获取海底浅部地质结构的探测技术，在海洋科学研究、海洋工程调查中发挥着重要作用。本文针对近岸海域浅地层剖面随机噪声和多次波发育的问题开展研究，首先分析了发育成因和影响因素，在此基础上，采用空间域滤波和频率域滤波相
2、结合的方式对随机噪声衰减，通过中值滤波法对浅地层剖面影像的预处理，消除孤立、随机的噪声点；再采用二维矩形预测算子在频率域滤波，滤除噪声频带，衰减背景噪声，提高信噪比；然后，利用一次波与多次波之间的周期相关性，采用改进型的预测反褶积对多次波处理，并给出预测算子长度和预测步长的选取方式，提高浅地层剖面数据成像质量。在近岸海域实际工程应用中，通过本方法和常规方法的效果对比发现，本方法可有效衰减随机噪声，压制多次波效果明显，展示更多隐藏的海底管线和浅部地质结构，说明了本方法的有效性、可靠性和实用性。关键词：近岸海域；浅地层剖面；随机噪声；多次波中图分类号：文献标识码：文章编号：（）：（，）：，第期张体
3、强：近岸海域浅地层剖面数据处理技术及应用 ，：；引言随着海洋强国战略持续推进，近岸涉海工程建设越来越多。为确保工程顺利安全施工，建设前期需探明海底管线状况和海底地质结构。浅地层剖面仪作为一种重要的海底声学探测设备，以其高分辨率、可穿透地层、可连续走航式探测、操作方便等优点，在海洋调查中的应用越来越广泛。近年来，浅地层剖面探测仪器已有了较大发展，探测数据仍主要依赖于仪器配套的软件处理，主要包括数据回放、粗差修正、参数修改等功能，少数软件具有辅助层界划分功能。浅地层剖面数据采集时，易受仪器设备与外部环境噪声的干扰影响，产生随机噪声，数据信噪比低；再者，部分海域海底地层复杂，声呐信号在自由界面与地层
4、之间来回反射，形成多次波，导致浅地层剖面数据成像模糊，有效信息被多次波掩盖，影响海底管线识别和海底地层划分。施凤应用滑动均值滤波、改进的中值滤波等方法对浅地层剖面图噪声处理，以此提高图像信噪比，提出了改进预测反褶积法多次波压制方法，突出有效波；王国强研究了浅地层剖面多次波自动搜索与压制方法，设计了多次波压制流程；刘宏扬等、赵铁虎等详细介绍了几种常见的浅地层剖面干扰因素及声学特征，给出了合理化建议；赵波等提出了基于同相轴追踪的多次波剔除与多道维纳滤波相结合的衰减方法，验证了该方法在多次波预测与匹配衰减方法的有效性。由于目前涉海工程大多位于近岸海域，相比于远海，外界环境更为复杂，海流更为紊乱，随机
5、噪声也相对更为发育；再者，近岸海域水深较浅，地层更为复杂，多次波也相对更为发育。因此，有必要在前人研究的基础上，对近岸海域复杂条件下的浅地层剖面数据处理技术进行有针对性的研究。本文针对近岸海域浅地层剖面随机噪声和多次波发育的问题开展深入研究，对其产生的机理进行分析，采用了空间域和频率域滤波相结合的方式衰减随机噪声，提高信噪比；再利用预测反褶积滤波法压制多次波，进而提高浅地层剖面数据成像质量。随机噪声和多次波发育成因浅地层剖面仪主要包括震源装置、收发基阵和记录显示系统，其中震源装置、收发基阵位于固定或拖拽于船底，采用连续走航式探测，以声学剖面图形记录。在近岸海域探测过程中，浅地层剖面影像容易受到
6、船体自身、船速、安装位置、海况、海流、过往船舶等多种因素的影响，在浅地层剖面影像上出现一些随机的、离散的、孤立的像素点，即随机噪声（图）。船体不稳，船速不固定，安装螺旋桨越近，海况越差，顺流逆流掌控不好，过往船舶越多，随机噪声越明显。浅地层剖面探测是单道自激自收，遇到海底面直接反射叫做一次“真”波，但由于声波在海水面、海底面，海底地层面等介质中传播形成波阻抗的强反射界面，经再次反射在剖面图上呈现非垂 工 程 地 球 物 理 学 报（）第 卷图随机噪声发育的浅地层剖面影像 图多次波发育的浅地层剖面影像 直反射引起的波叫做二次声波，多次反射形成多次波，。多次波与一次波走时相近时，会干扰甚至掩盖一次
7、波有效信号，影响成像质量，尤其在近岸海域，水深越浅，多次波越发育（图），若不能有效压制，将会导致深埋管线和地层信息被覆盖而无法识别。浅地层剖面数据处理技术随机噪声处理是浅地层剖面数据首先需要解决的重要问题。为了在不损失有效波的情况下尽可能消除随机噪声，本文先采用空间域滤波对浅地层剖面影像进行预处理，消除孤立、随机的噪声点，再进行频率域滤波噪声衰减。中值滤波是基于排序统计理论的非线性平滑技术。基本原理是：将数字图像中一点的像素值用该点邻域中各点像素值的中值代替，从而消除孤立噪声点 。随机噪音在频率域的频带比有效波宽，可滤除其对应频带去噪。本文引用李国以二维矩形预测算子代替常规随机噪声衰减技术（，
8、）使用的一维预测算子，并利用浅地层剖面数据的相关性增强先验信息，以提高数据处理质量。首先，对浅地层剖面数第期张体强：近岸海域浅地层剖面数据处理技术及应用据进行频率域转换，再在每个频段数据选择合适时窗和算子长度，沿方向做预测滤波去噪，然后将处理后的数据转换回空间域，完成滤波去噪。数据时窗选择和预测算子长度决定了随机噪音的效果。数据时窗是每次处理的道总数，当有效波同相轴为直线时，预测误差较小；当有效波同相轴为曲线时，预测误差较大。预测算子越小，随机噪声衰减越大，但有效波损失越多；预测算子越大，有效波损失越少，但随机噪声衰减也越小。因此，在数据处理时，应选取不同预测算子长度尝试，确定最佳长度，提高信
9、噪比。多次波压制技术多次波在近岸海域浅地层剖面探测中是干扰海底管线和地层划分影响最大的波。根据多次波具有周期的特性，采用预测反褶积对浅地层剖面数据中多次波进行处理，其基本原理是利用一次波与多次波之间的相关性，预测出多次波，再减去预测 出 的 多 次 波，从 而 得 到 一 次 反 射 波 信号，。浅地层剖面的褶积模型可近似表示为：（）（）其中，（）是地震记录；是反射系数序列；是地震子波。反褶积是这个褶积过程的逆过程，其中较常用是预测反褶积。预测反褶积公式如下：（）（）（）（）其中，（）为时间为时的预测值；为预测步长（）；（）为预测算子；为预测算子长度（）。由式（）可以得出，多次波压制效果跟预测
10、算子长度和预测步长有关。在实际数据处 理 中，一 般 预 测 滤 波 算 子 长 度 取 ；采用自相关函数来选择预测步长，如当预测滤波算子长度为，预测步长为时，自相关长度不得小于。由于预测反褶积对浅地层剖面数据中一阶多次波压制效果较好，但当浅地层剖面数据中存在二阶、三阶多次波等高阶多次波时，效果减弱。为增强压制效果，可在预测得到一阶多次波后，再和预测算子做褶积得到二阶多次波，以此类推，得到高阶多次波；然后将多次波叠加移除，得到一次波。该方法对压制复杂介质情况下的多次波效果较好。但预测的高阶多次波数越髙，删除的信息也越多，有可能损失有效一次波。因此，在实际数据处 理 时，应 尝 试 不 同 阶
11、数 效 果，选 择 最 佳阶数，。应用案例以某沿岸建设工程为例，为保障项目前期方案制定，及后期施工顺利实施，需对整个工区范围的海底管线和地质条件进行调查。工区为沿海段，水深，海底表面以淤泥质滩涂为主。本项目采用 浅地层剖面系统，重量轻，方便携带和运输，适合在小船浅水作业，采集频率 。由于工区位于近海岸海域，探测时海况较差，时常为避让船只偏离测线，改变航速，以至于随机噪声和多次波较为发育，掩盖了海底面下的管线识别和地层划分。图（）为浅地层剖面原始记录影像图；图（）为单独采用中值滤波的浅地层剖面影像图；图（）为单独采用频率域滤波后的影像图，其中，数据时窗和预测算子长度分别为 和；图（）为本文空间域
12、滤波和频率域滤波相结合处理后的影像图。通过对比图各图可知，图（）的背景噪声发育，地层信息模糊不清，界面不连续；图（）中虽然背景噪声仍然存在，但影像得到明显改善，地层界面逐渐显现，孤立、随机的噪声点基本消除；图（）的背景噪声在很大程度上得到了压制，地层界面更为清晰，由于孤立、随机的噪声点的存在导致部分地层不连续；图（）的背景噪声基本消除，地层界面清晰可辨，地层连续性较好，展示了更多的细节信息，海底界面辨识别度高。图为经过多次波压制技术处理后的影像图，其 中，预 测 算 子 长 度 和 预 测 步 长 ，虽然海底面和地层信息变得相对模糊，但在影像图中的和两个位置的异常反射弧较为明显，走时分别为 和
13、 ，通过不同测线多次验证，在相同位置均出现了类似的异常反射弧，并与已有海底管线资料对比，验证了此两处为海底管线。本实例的应用说明了本方法的有效性和可靠性。工 程 地 球 物 理 学 报（）第 卷图地层剖面原始记录及采用不同滤波处理后的浅地层剖面影像 第期张体强：近岸海域浅地层剖面数据处理技术及应用图多次波压制处理后的影像 结论本文通过对近岸海域浅地层随机噪声和多次波发育的问题开展研究，首先分析了产生随机噪声的主要干扰因素和多次波发育的原因；在此基础上，采用空间域中值滤波和频率域滤波相结合的方式衰减随机噪声，提高信噪比；然后，采用改进型的预测反褶积对浅地层剖面数据多次波处理，增强有效信息的可识别
14、性。在实际资料处理中采用本方法，浅地层剖面影像质量得到明显改善，辨识度高，可展示更多海底面以下管线和地层的信息，为方案设计和后续施工提供技术支撑。参考文献：郭军，陈宏文，李邵波，等探测海底内部地质构造的浅地 层 剖 面 系 统 测 绘 地 理 信 息，（）：张彬彬，赵国兴，吴永亭，等浅地层剖面仪噪声特征分析及采集参数 优选 海洋科 学 进 展，（）：杨国明，朱俊江，赵冬冬，等浅地层剖面探测技术及应用海洋科学，（）：施凤浅地层剖面数据精处理关键技术研究武汉：武汉大学，王 国 强浅 地层剖面数 据 多次波 自 动 搜 索 与 压 制北京：中国地质大学，刘宏扬，张异彪，李斌，等 浅地层剖面测量常见干
15、扰因素分析海岸工程，（）：赵铁虎，张志珣，许枫 浅水区浅地层剖面测量典型问题分 析 物 探 化 探 计 算 技 术，（）：赵波，谭军，李金山，等 浅地层剖面的自由界面多次波预测与衰减中国海洋大学学报（自然科学版），（）：，：，（）：李平，杜军浅地层剖面探测综述海洋通报，（）：李振春，张军华地震数据处理方法东营：石油大学出版社，吴玉莲图像处理的中值滤波方法及其应用西安：西安电子科技大学，郝继杰保中心值滤波及其效果分析青岛：山东科技大学，：，（）：李国发全三维随机噪声衰减技术石油地球物理勘探，（）：王兆旗 预测反褶积及其相关问题分析 西安：长安大学，（）：张军华，缪彦舒，郑旭刚，等预测反褶积去多次波几个理论问题探讨物探化探计算技术，（）：唐海敏，苗庆库，王银 预测反褶积在压制周期性多次反射波中的应用研究 铁道勘察，（）：赵夏，金文昱，李振勇 应用反褶积技术压制多次波的实例分析石油地球物理勘探，（）：
