第３８卷第２期 华 南 地 震 ＳＯＵＴＨ ＣＨＩＮＡ １ｏＵ１ＬＮＡＬ ＯＦ ＳＥＩＳＭｏＬｏＧＹ Ｖｏ１．３８．Ｎ０．２ ２０１８年６月 Ｊｕｎ．，２０１８ 张旭东，尹成，陈玺．垂直观测系统海洋地震资料信号处理研究【Ｊ１．华南地震，２０１８，３８（２）：１－７．［ＺＨＡＮＧ Ｘｕｄｏｎｇ，ＹＩＮ Ｃｈｅｎｇ，ＣＨＥＮ Ｘｉ Ｓｔｕｄｙ ｏｎ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｉｎｇ ｏｆ Ｍａｒｉｎｅ Ｓｅｉｓｍｉｃ Ｄａｔａ ｉｎ Ｖｅｒｔｉｃａｌ Ｏｂｓｅｒｖａｔｉｏｎ Ｓｙｓｔｅｍ［Ｊ１．Ｓｏｕｔｈ Ｃｈｉｎａ ｊｏｕＦｎａｌ ｏｆ ｓｅｉｓｍｏｌｏｇｙ，２０１８，３８（２）：１－７】 垂直观测系统海洋地震资料信号处理研究 张旭东 ，尹成 ，陈玺 （１．油气藏地质及开发工程国家重点实验室西南石油大学．成都６１０５００： ２．国土资源部海底矿产资源重点实验室广州海洋地质调查局．广州５１００７５） 摘要：海洋地震资料处理中，通常采用的提取远场子波方法是从地震资料中直接提取和根据自由气泡震荡理 论中不相干子波叠加算法由近场子波推算远场子波，但这两种方法都不精确 首先分析了由ＯＢＳ模拟垂直观 测地震系统实测得到远场子波，进而获取零相位算子（包括零相位算子、压制气泡算子和压制虚反射算子）的 原理和实际效果。利用理论模型和实际数据对所述方法进行了检验，得到的结果不仅消除虚反射据处理效果来看此方法是切实可靠的 关键词：垂直观测系统；零相位算子：虚反射：气泡效应 中图分类号：Ｐ６３１．４６ 文献标志码：Ａ 文章编号：１００１—８６６２（２０１８）０２—０００１—０７ ＤＯＩ：１０．１３５１２／ｊ．ｈｎｄｚ．２０１８．０２．００１ 气泡效应 影响，而且整形为零相位。通过上述处理消除了非零相位数据对高分辨率处理和解释标定的影响．从实际数 、Ｓｔｕｄｙ ｏｎ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｉｎｇ ｏｆ Ｍａｒｉｎｅ Ｓｅｉｓｍｉｃ Ｄａｔａ ｉｎ Ｖｅｒｔｉｃａｌ Ｏｂｓｅｒｖａｔｉｏｎ Ｓｙｓｔｅｍ ＺＨＡＮＧ Ｘｕｄｏｎｇ ，ＹＩＮ Ｃｈｅｎｇ ，ＣＨＥＮ Ｘｉ （１．Ｓｔａｔｅ Ｋｅｙ Ｌａｂｏｒａｔｏｒｙ ｏｆ Ｏｉｌ ａｎｄ Ｇａｓ Ｒｅｓｅｒｖｏｉｒ Ｇｅｏｌｏｇｙ ａｎｄ Ｅｘｐｌｏｉｔｔｉａｏｎ，Ｓｏｕｔｈｗｅｓｔ Ｐｅｔｒｏｌｅｕｍ Ｕｎｉｖｅｒｓｉｔｙ， Ｃｈｅｎｇｄｕ ６１０５００，Ｃｈｉｎａ；２．Ｋｅｙ Ｌａｂｏｒａｔｏｒｙ ｆＭａｒｉｏｎｅ Ｍｉｎｅｒａｌ Ｒｅｓｏｕｒｃｅｓ，Ｍｉｎｉｓｔｒｙ ｆＬａｎｄ ｏａｎｄ Ｒｅｓｏｕｒｃｅｓ， Ｇｕａｎｇｚｈｏｕ Ｍａｒｉｎｅ Ｇｅｏｌｏｇｉｃａｌ Ｓｕｒｖｅｙ，Ｇｕａｎｇｚｈｏｕ ５１００７５，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｉｎ ｏｃｅａｎ ｓｅｉｓｍｉｃ ｄａｔａ ｐｒｏｃｅｓｓｉｎｇ，ｔｈｅ ｃｏｍｍｏｎｌｙ ｕｓｅｄ ｍｅｔｈｏｄ ｏｆ ｅｘｔｒａｃｔｉｎｇ ｆａｒ ｆｉｅｌｄ ｗａｖｅｌｅｔ ｉｓ ｔｏ ｄｉｒｅｃｔｌｙ ｅｘｔｒａｃｔ ｔｈｅ ｆａｒ－ｆｉｅｌｄ ｗａｖｅｌｅｔ ｆｒｏｍ ｔｈｅ ｓｅｉｓｍｉｃ ｄａｔａ ａｎｄ ｃａｌｃｕｌａｔｅ ｔｈｅ ｆａｒ——ｉｆｅｌｄ ｗａｖｅｌｅｔ ｆｒｏｍ ｔｈｅ ｎｅａｒ——ｆｉｅｌｄ ｗａｖｅｌｅｔ ａｃｃｏｒｄｉｎｇ ｔｏ ｔｈｅ ｉｒｒｅｌｅｖａｎｔ ｗａｖｅｌｅｔ ｓｕｐｅｒｉｍｐｏｓｉｔｉ０ｎ ａｌｇｏｒｉｔｈｍ ｉｎ ｔｈｅ ｆｒｅｅ－ｂｕｂｂｌｅ ｏｓｃｉｌｌａｔｉｏｎ ｔｈｅｏｒｙ． Ｈｏｗｅｖｅｒ，ｎｅｉｔｈｅｒ ｍｅｔｈｏｄ ｉｓ ａｃｃｕｒａｔｅ．Ｉｎ ｔｈｉｓ ｐａｐｅｒ，ｔｈｅ ｐｒｉｎｃｉｐｌｅ ａｎｄ ｔｈｅ ｐｒａｃｔｉｃａｌ ｅｆｆｅｃｔ ｏｆ ｔｈｅ ｚｅｒｏ—ｐｈａｓｅ ｏｐｅｒａｔｏｒ（ｉｎｃｌｕｄｉｎｇ ｚｅｒｏ—ｐｈａｓｅ ｏｐｅｒａｔｏｒ，ｄｅ—ｂｕｂｂｌｅ ｏｐｅｒａｔｏｒ ａｎｄ ｄｅ—ｇｈｏｓｔ ｏｐｅｒａｔｏｒ）ｏｂｔａｉｎｅｄ ｂｙ ｍｅａｓｕｒｉｎｇ ｔｈｅ ｆａｒ—ｆｉｅｌｄ ｗａｖｅｌｅｔ ｏｂｔａｉｎｅｄ ｆｒｏｍ ｔｈｅ ｓｉｍｕｌａｔｅｄ ｖｅｒｔｉｃａｌ ｃａｂｌｅ ｓｅｉｓｍｉｃ ｓｙｓｔｅｍ ｏｆ ＯＢＳ ａｒｅ ａｎａｌｙｚｅｄ．Ｕｓｉｎｇ ｔｈｅ ｔｈｅｏｒｅｔｉｃａｌ ｍｏｄｅｌ ａｎｄ ｔｈｅ ａｃｔｕａｌ ｄａｔａ．ｔｈｅ ｍｅｔｈｏｄ ｄｅｓｃｒｉｂｅｄ ｉｎ ｔｈｉｓ ｐａｐｅｒ ｉｓ ｔｅｓｔｅｄ．Ｔｈｅ ｒｅｓｕｌｔｓ ｏｂｔａｉｎｅｄ ｎｏｔ ｏｎｌｙ 收稿日期：２０１６—１０—２６ 基金项目：中国地质调查局项目（ＤＤ２０１８９６４２） 作者简介：张旭东（１９８０一），男，硕士，高级工程师，主要从事地震资料成像及反演研究。 Ｅ－ｍａｉｌ：ｚｈａｎｇｘｄ９１ １＠１２６．ＣＯＩＴＩ． ＿２ 华 南 地 震 ３８卷 ｅｌｉｍｉｎａｔｅ ｔｈｅ ｅｆｆｅｃｔｓ ｏｆ ｇｈｏｓｔ ａｎｄ ｂｕｂｂｌｅ，ｂｕｔ ａｌｓｏ ｓｈａｐｅ ｔｈｅ ｓｅｉｓｍｉｃ ｄａｔａ ｔｏ ｚｅｒｏ ｐｈａｓｅ．Ｔｈｅ ａｂｏｖｅ ｐｒｏｃｅｓｓ ｅｌｉｍｉｎａｔｅｓ ｔｈｅ ｉｎｆｌｕｅｎｃｅ ｏｆ ｎｏｎ—ｚｅｒｏ ｐｈａｓｅ ｄａｔａ ｏｎ ｈｉｇｈ—ｒｅｓｏｌｕｔｉｏｎ ｐｒｏｃｅｓｓｉｎｇ ａｎｄ ｉｎｔｅｒｐｒｅｔａｔｉｏｎ ｃａｌｉｂｒａｔｉｏｎ，ａｎｄ ｔｈｅ ｍｅｔｈｏｄ ｉｓ ｒｅｌｉａｂｌｅ ｆｒｏｍ ｔｈｅ ａｃｔｕａｌ ｄａｔａ ｐｒｏｃｅｓｓｉｎｇ ｅｆｆｅｃｔ． Ｋｅｙｗｏｒｄｓ：Ｖｅｒｔｉｃａｌ ｃａｂｌｅ；Ｚｅｒｏ ｐｈａｓｅ ｏｐｅｒａｔｏｒ；Ｇｈｏｓｔ；Ｂｕｂｂｌｅ ｅｆｆｅｃｔ ０引言 地震子波提取方法总体上可以分为确定性子 波提取方法和统计性子波提取方法 确定性子波 空间最大尺寸为３０ ｎ１．声波在水中速度取１ ５００ ｍ／ｓ．对于最大频率２００ Ｈｚ的子波．求得的远场距 离为１ ２０ ｍ．而实际垂直观测系统距离震源最近距 离也有１ ８００ ｍ，完全满足远场要求。 提取方法不需要对反射系数分布作任何假设，能 得到较为准确的子波 统计性子波需要对地下地 震资料和反射系数分布进行假设．得到子波精度 与假设条件满足程度有关 传统方法利用海洋地震资料直接提取地震子波， 受不确定因素影响较大，效果不太理想 近年来。 一些海洋地震采集通过在阵上放置近场检波器获 取近场子波．再根据自由气泡震荡理论中不相干子 波叠加算法求解方程组．获得远场子波 但是上述 两种方法均基于一定假设．得到的远场子波精度有 限 ２０１５年广州海洋地质调查局在南海进行了模拟 垂直观测系统试验．利用垂直观测系统在深水海底 实测获得精确远场子波 利用精确远场子波可以进 行确定性子波提取．而且利用实测得到的远场子波 求取零相位化算子（包括零相位算子、压制气泡算子 和压制虚反射算子）并用于实际地震资料１８－９１。取得了 不错的效果 为后续地震资料处理和地质解释提供 高质量的资料 ｝ １垂直观测系统数据及提取地震子波 方法 海洋地震勘探的气阵列中各之间激发互 相不影响．各自单独产生单震源子波．在近场 检波器中记录到的信号．是经过传播距离ｒ后直接 到达检波器的单震源子波和经过传播距离尺后 接收到的虚反射叠加而成．如图１所示 当检波 器和阵距离较近时，也就是近场区时．ｒ和Ｒ相 差较大，由于球面扩散使得单震源子波与虚反 射的振幅差异较大．随着检波器和震源之问距离 的增加，虚反射的传播路径逐渐接近震源子波路 径，这时两者受到的球面扩散引起的衰减近似相 ｎ ２ 同，当这个距离ｄ满足经验公式：ｄ＞ ： ， Ａ… （Ｊ 其中厂为激发子波频率：Ｄ为激发子波最小波长： 为水中声波速度。对于本次试验来说．气阵列 １震源子波和虚反射传播路径 ２０１５年广州海洋地质调查局利用多台ＯＢＳ （Ｏｃｅａｎ ｂｏｔｔｏｍ ｓｅｉｓｍｏｇｒａｐｈ 模拟实现垂直观测系 统 ＯＢＳ模拟垂直观测系统总长度２００ ｍ．接收节 点数量３个．节点间距２５ ｍ．其整体结构如图２ 所示．ＯＢＳ采样频率最高达１ ｋＨｚ．动态范同大于 １２０ ｄＢ，最长＿丁作时间可达２４ ｄ（４ ｍｓ） ＯＢＳ最 大工作水深５ ０００ ｎ／．可以记录水听器分量（Ｐ分 量）、垂直分量（ｚ分量）、水平分量（ 分量）和水 平分量（ｙ分量） ＯＢＳ模拟垂直观测系统采集参 数见表１ ２５ｍ ２５ｍ 冈２垂直观测系统结构图 Ｆｉｇ．２ Ｓｔｒｕｃｔｕｒｅ ｏｆ ｖｅｒｔｉｃａｌ ｃａｂｌｅ 第２ 张 东等：垂血脱测系统 海洋地震资料处 ｌｆ１的心川 表１垂直缆地震系统采集参数 Ｔａｂｌｅ １ Ａｃｑｕｉｓｉｔｉｏｎ ｐａｒａｍｅｔｅｒ ｏｆ ｖｅｒｔｉｃａｌ ｃａｂｌｅ ｓｅｉｓｍｉｃ ｓｙｓｔｅｍ 点数 ３ 源 放深度／ｌｌ ５ 节 ，，￣ＩＮｈｑ ／ｍ ｉ 录长度／ｍｓ 炮间距／ｎ 气怆类， 采样率／ｍ 伧阵容 本次试验…ｔ，０ＢＳ模拟 直观测系统与常规拖 缆联合采集．扶得ｒ天然气水合物立体探测地震数 试验 采集测线ｌ７条．每条长ｌ７ ｋｍ．采集面 １５４ ｋｍ ．垂直观测系统投放位置位于采集Ｔ区 ＩｆＩ心 ３为９４９ ６００测线Ｌ｛ｌ１ 个接收节点水Ｉ听器 分量（Ｐ分 ）的共接收点道集数据 从 ３中可以 ７亍ｆｆ；，Ｊ：、ＬＩＪ、下 个ＯＢＳ记录的有效波、多次 也是垂直观测系统采集地震资料优于常舰地震资料 之处。ＯＢＳ接收节点的直达波到达Ｉ１，￣１ｉｉ１与其节点深 度吻合很好，直达波和海底反射波同州轴卡｝Ｉ百尢 涉。在远偏移 处，直达波与虚反射乖Ｉｌ初至波分离 明显。经过反复比较．从１号ＯＢＳ ｊ 选取离（）Ｂｓ 模拟乖直观测地震系统最近的９４９ ６００测线数抛． 波、折射波郁非常清晰．而几．时１于接收节点靠近海 底且位置＿牛ＩＩ对 定．采集数据具 较高信噪比．这 将其直达波拉平，选取波形稳定的一段叠加后俄取 ５００ Ｉｌｌｓ．即为提取的远场子波 【】 ４） 篁 兰 百 未 ＿ ｒ＾ ＿０ 主 １ ０ ３（】 ５．０ ７（） ３ ＯＢＳ接收（测线９４９６００）水听器分黾 接收点道集数 （ａ．１ ＯＢＳ，ｈ．２号ＯＢＳ，ｒ．３ （） ） Ｆｉｇ．３（ ‘）ｌｎｌｌｌＯＩ１ ｒｅ‘’ｅｉｖｉｎｇ ｐｏｉｎｌ ｇａｔｈｅｒ ｄａｔａ ｒｅｃｅｉｒｅｄ ｂｙ ０ＢＳ（１ａｎｅ ９４９６００）ｌ１ｙ‘ｌｒｌ１Ｉ）Ｉｌｏｌｌｔ￣ｔ‘‘）ｉｉｉｐｏｎ￣ｌｌ［ （ａ．Ｎｏ．１ ＯＢＳ，｝）＿ＮＯ．２ ＯＢＳ， Ｎｏ．３ ＯＢＳ） ２远场子波的应用 远场子波最主要Ｖ，Ｊ￣ｉ！Ｌ，ｑＪ足求收零卡Ｈ值化 ｒ． 再把求取的零卡Ｈ位化箅了应川到地震资料处 『『｝． 实现地震资料的零卡日位化处 、地震资料零卡Ｉ｛位 化算予Ｆ｝１爪制气泡箅子＿不１Ｉ整形箅－，两部分绀成．憋 形箅子又包括零相位整形箅予和 制虚反射箅 ２．１压制气泡算子求取 ＭＩｈ］／『ｌ１ｓ 气阵列ｍ具有小同容量、不同型 ‘，能够 输Ⅲ 同能量、不同频率的 波形经过优化组 合而成．ｈｉ以实现卡Ｈ十』ＪＩ】强 最大限度爪制气泡， ４ ＩｒＥＩ ［脱测系统数据提取的子波 Ｆｉｇ．４ Ｗａｗｄｅｔ ｅｘｔｒａｃｔｅｄ ｆｒｏｎｌ ｖｅｉ＇ｌｉｃａｌ ｃａｂｌｅ（１ａｌａ ４ 华 南 地 震 ３８卷 从而获得具有所需频带特征的子波．最终实现对 目标层位的有效勘探１２５］ 尽管阵组合在设计时就 充分考虑了对气泡的压制．但如图５ａ所示的气 子波还是会受到气泡效应的影响．在频谱上气泡 相加特性可知．输入的带电缆鬼波的远场子波的 相位谱（图７ａ）与零相位化算子的相位谱（图７ｂ） 是相反的。二者相加的结果应趋于零相位（图７ｃ） 这一特性可作为检验零相位算子（图８）正确性的一 效应表现在低频端，见图６ａ．这是地震勘探不需 要的．而且气泡效应具有明显的周期性．可以通 过预测反褶积来消除 另外．从褶积过程中相位 种方法．也可以将求取的零相位化算子应用于输 入子波（图５ａ），得到期望的零相位信号（图９）来 验证求取算子的正确性 ３０Ｉ ，ｌ １Ｉ１ ０ ｝ 寸ｌｆ ４／ Ｏ． ３ ０Ｏ ａ ＪｌｏＣ ３１ ）０ ｌｏ【 １厂 －３ ＩＩ １Ｉ１ ０ 寸ｆ 习， ２ｔ ｌｏ ３｛： １Ｏ０ ▲ ｆ ｌ －ｂ ３ ） 图５带电缆鬼波的远场子波（ａ）及去气泡后的带电缆鬼波的远场子波（ｂ） Ｆｉｇ．５ Ｆａｒ－ｆｉｅｌｄ ｗａｖｅｌｅｔ ｗｉｔｈ ｇｈｏｓｔ ｗａｖｅ ｏｆ ｃａｂｌｅ（ａ）ａｎｄ ｆａｒ—ｆｉｅｌｄ ｗａｖｅｌｅｔ ｗｉｔｈ ｃａｂｌｅ ｇｈｏｓｔ ｗａｖｅ ａｆｔｅｒ ｂｕｂｂｌｅ（ｂ） 名 罟 号 蠢 ｌ 图６ ａ、ｂ分别为图５中ａ、ｂ所示子波频谱 图７ ａ．去气泡后相位谱ｂ．零相位化算子相位谱ｃ．零相位子波相位谱 Ｆｉｇ．７（ａ）Ｐｈａｓｅ ｓｐｅｃｔｒｕｍ ｏｆ ｄｅ ｂｕｂｂｌｅ（ｂ）Ｐｈａｓｅ ｓｐｅｃｔｒｕｍ ｏｆ ｚｅｒｏ ｐｈａｓｅ ｏｐｅｒａｔｏｒ（ｃ）Ｐｈａｓｅ ｓｐｅｃｔｒｕｍ ｏｆ ｚｅｒｏ ｐｈａｓｅ ｗａｖｅｌｅｔ 第２期 张旭东等：垂直观测系统在海洋地震资料处理中的应用 时间／ｍｓ ６０ 一ｌ５Ｏ 一１００ －５０ ｝ ５０ １００ ｌ！ Ｏ－ ５ ｌｑｒ ．■▲Ｊ …Ｖ ， － ６０ 冈８零相位化算子 Ｆｉｇ．８ Ｚｅｒｏ ｐｈａｓｅ ｏｐｅｒａｔｏｒ 时问／ｍｓ ２（ ）Ｏ ４５ １Ｃ １０ （ １Ｃ ｌ０ ２（） ０ １００ Ｊ Ｉ ＪＵ －＿Ｉ，．一 ４５ 图９零相位子波 Ｆｉｇ．９ Ｚｅｒｏ ｐｈａｓｅ ｗａｖｅｌｅｔ 首先地震记录的褶积模型为： （ ）＝（￡Ｊ（ ） ｅ（ ）＋ｎ（ ） （１） 这样最终的零相位化算子可以南上面两个算 子得到： 其中ｘ（ｔ）为地震记录，ｗ（ｔ）为地震子波，ｅ（ｔ） 为地层脉冲响应，ｎ（ ）为随机噪音。当我们假设为 ｎ（ｔ）Ｏ，震源波形已知时，公式（１）可以表示为： ｘ（ｔ）＝ｔｏ（ｔ） ｅ（ｔ） 这样公式（２）可以改写为： （ ）＝ｓ（ｔ） （ ）＋ｅ（ ） （３） （２） 后的地震记录 ＿厂（ ）＝０（ ） Ｃ（ｔ） （６） 将／（丁）与地震记录 （ｆ）褶积，得到零相位化 ２．２整形算子求取 海洋地震资料中．由于震源激发和电缆检波 器均在海面以下．容易产生虚反射。由于海水与 海上地震勘探可以记录震源组合的远场子波． 其中ｓ（ｔ）为记录的震源组合远场子波，采用 震源子波反褶积方法既公式（４）求取去气泡因子 ａ（ｔ）。 １ ｒ ｘ空气形成了反射系数近似为一１的自由表面，因此 检波点除了接收由震源激发．并经过地层反射， 直达检波点的一次反射波（图ｉＯａ）以外，还会接收 到伴随在一次反射波后面的三种干扰波（图ｌＯｂ、 图１０ｃ、图ｌＯｄ） 图１０ｂ为先经过海面再经过地 Ｉ ｒ （） （０） 一 １… （ｍ一）Ｉ １ ｎ（ ）Ｉ一 （ｒｘ ｎ十１） ｒ㈩ 层接收到的虚反射．图ｌＯｃ为先经过地层再经过 １ ｒｘｘ（ｍ）ｒｘｘ（ｍ一１）… （０）Ｊｌ。（ｍ）Ｊ【ｒｘ （ｎ＋ｍ）ｊ 经过反褶积压制气泡效应后，子波的频谱得 到明显改善，见图６。这样得到的子波可以采用公 式（５）来求取整形算子ｃ（ｔ）。 ｆ（Ｉ＋Ａ） ，（０） 【 ｎ （ｍ） ｍ，（１） … ｒｂｌ，（ｍ） ］ｆ ｃ（Ｏ）］ 海面接收到的虚反射．图ｌＯｄ为先经过海面再经 过地层然后又经过海面接收到的虚反射。本文与 垂直观测系统同时采集的拖缆数据震源深度和电 缆深度均为５ ｍ．这就使得震源激发的虚反射与 电缆检波器接收的虚反射叠加在一起。理论模型 试算效果见图１１。图１１ａ为反射系数序列，图 ｒｂｂ（１） （Ｉ＋Ａ）ｒｂｂ（０）… ｒｔ￣（ｍ－１） ｃ（１）ｌ一 ｒ｝ （” 一１） ···（１＋ｈ）ｒ （０）Ｊ【ｃ（ｍ）Ｊ ｆ ｂ（Ｏ）］ １ ０ … １ ｌｂ为受虚反射影响反射系数序列，图１ｌｃ为雷 克子波与受虚反射影响反射系数序列褶积得到合 成地震记录．图１ｌｄ和图１１ｅ分别为合成地震记 录零相位化和压制虚反射后结果，图１ ｌｇ为压制 虚反射前后频谱 ｌ １（５） 其中ｒ 为输入子波的自相关函数，ｍ为自相 关函数的长度，ｃ（ ）为待求的整形算子。 ６ 华 南 地 震 ｆａ ’ｌ ｌｂ ｆ ｃ 『ｄ ｌ● ｌｅ ｌ· ‘ ＝ ： ■－－—＿－ ■＿‘■‘·＿－－－ 《＝＝ ｒ ＜ ■■■－——＿－ 亡＝＝ ． ／ ／ ａｃ 、Ｖ ／ 图ｌＯ检波器接收虚反射物理过程 Ｆｉｇ．１０ Ｐｈｙｓｉｃａｌ ｐｒｏｃｅｄｕｒｅｓ ｏｆ ｇｈｏｓｔ ｒｅｃｅｉｖｅｄ ｂｙ ｄｅｔｅｃｔｏｒ 图ｌ１合成地震记录零相位化及压制虚反射结果 Ｆｉｇ．１ １ Ｄｅ—ｇｈｏｓｔ ａｎｄ ｚｅｒｏ ｐｈａｓｅ ｒｅｓｕｌｔ ｂｙ ｓｙｎｔｈｅｔｉｃ ｓｅｉｍｏｇｒａｍ ３实际应用效果分析 本文所采用的实际资料位于南海北部陆坡区 ４结语 在海洋地震资料处理中．实际获得的气子 的琼东南盆地，琼东南盆地存在多个生烃凹陷． 该盆地已证实是我国近海大型含气区之一 图１２ 是远场子波零相位化（压制气泡算子、零相位算 子、压制虚反射算子）前后效果，通过震源子波反 褶积处理，远场子波的气泡和虚反射得到很好压 制，子波整形为零相位叠加剖面，波组特征变得 更加清晰，同相轴连续性得到明显改善 波并不是人们期望的形态，而是被虚反射和气泡 效应改造过的子波。利用本文的ＯＢＳ模拟垂直观 测系统提取远场子波处理技术可以基本消除虚反 射和气泡效应影响，同时子波整形为零相位 该 方法获得的零相位资料相位稳定、可靠．能为地 震资料拼接处理、高分辨率处理等打好基础．同 时为准确标定油气藏顶底界面、精确落实圈闭范 围和尖灭追踪等精细解释工作提供保障 但该方 法仍然残留部分虚反射．值得今后进一步研究 第２ ｆｊｃ Ｊ张 尔等： 观测系统 海ｆ 地旋资料处 ｌｆｌ的Ｊ、 川 ７ １２ Ｊ 制席反射种ｒ冬ｆｌ１倚化 （ 、Ｊ （下）稳加削 ｌ：ｉｇ．Ｉ ２．ｒ （）、 Ｐ ｒＪａｙ Ｉｌｒｏｌ＇ｉｌｅｓ（）ｆ－ｔｈ ｒ’一ｍ·ｈｉ ｗａ、，ｅｌｅｌ ｂｅｆｂｌ ｅ ａ Ｉ ａｔｌ＇ｅｔ’ｔｌ （１ｅ Ｉｍｈ１）ｌｃ ｔｉｌｌ（］ｇｈｏｓｔ 参考文献： ＩＩ１个水此．Ｉ钧 小 ． Ｉｌ４发， ．地 数 处Ｊ：ｌｔＴ力’法ＩＭＩ． ＩＬ ｊ ：￣／ｉｌｌｔ １ ，１１：：１１版社，２００７：５７～１００ Ｉ１６１ Ｋｅｌｈ ｒ Ｊ Ｂ， ｔ）ｄｌｌｔ－ｒ ｌ １． Ｉｍｌ，Ｉ－１ —ｔ‘．Ｉ】ａ１ｉ（Ｊ ｌｓ ＩＪｌ ，Ｉ ｌｌｌｌｄｔｑ ｗ ＩＩｔ·ｒ ＰＸ１）ｈ，ｓｉｏｎ Ｉｌｆ ｋｐｐｌ１ｉ｛·ＩＩ ｌ’Ｉ１、ｓｉ｛’ ． Ｉ２ｌ ｆＪ｝尔马滋． ．地 资料分析一地震资料处 、反演和解 １ Ｉ ７ ｌ Ｚｉｏｌｋｏｗｓｋｉ Ａ．Ｐａｒｋｓ（；，Ｈａｔｌ＜）ｎ ｌ Ｊ’ｔ ｔ ａ１． ｒｌ ｓｉｇｎａｈ． ０１ ｉｌｌｌ Ｖ：ｃｏｍ１）ｕｌｉｏｔ１ ｆｒｏｍ Ｉｌｌ！ｉｌｌ’一ｆｉｅｈｌ Ｉｎｔ＇ｉＬ＇；－ＩＩ１ ￣ｌｌｌ（ ｈｉｓ ａｉ ｒ’一ｇｕｎ ａｒｌ’ｔ１．释（Ｉ：升＂‘）Ｉ ．刘怀ＩｌＩ， 兜破，童思友．等｝币．北京： Ｈ Ｉ－．￣Ｌ 版 ，２００６：ｌ２２一ｌ７８． ｉｎ｛·ｌｍｌｉ￣ｉｇ ｉｎｔｅ］’ａｃｔｉｏｎｓ…．（ ｏ１）ｌＬｖｓｉ（·ｓ，１ ９８２ ４７（１（））：ｌ ４ｌ３一ｌ ４２１． Ｉ３１孙抓川， Ｉ＊１．－ 长晓渝，等．海洋拖揽地震资料处理关 键技术…．人然气１ ＇ｌｋ，２００７，２７（Ｓ１）：ｌ９２一Ｊ９４． ｉ亮，张义波，张淑协，等．深海拖揽 维地震资料处 技术ｌＪ１．ｆｎ巾地球物 勘探．２００６，４１（５）：５ｌ９－５２４． 抓 ， 浩，拿乃 Ｊ＇ｔ’．海洋地震数 处 技术探泔…． ， ＩＩ地球物理勘探．２００７，４２（Ｓ）．８—１３． １ｌ ８１ Ｚ．Ｉ）ｌｋｏｗｓｋｉ Ａ，Ｊｏｈｎｓｌｏｎ Ｊ．（］ｏｍｐｍａｔｉｏｎ ｏｆ ｆａ ｙ’一 ｈＩ ｄ＿ｒ ｇｎｉ１ ｓｉｇｎａｔｕｒｅｓ ｌｉｉ ｇｕｎ—ＩＩＩＯＵ］１ＩＰｔｌ ＩｎＰＳＳＩＩ ｒ（ Ｉ｜ｌｅｔｌｓＩｉｒ￣ ＩＩ１￣＇１ｉｔｓ…． Ｅｘｌｍｎｄｅ｛Ｉ ＡＩ）ｓｔｒａ（’ｔｓ ｏｆ ６６ｔｈ Ａｉｍ Ｌｍ］ｈｉｔ（ ｉ｜ｉｌｌ １９９６：１３—１６． （ Ｍｔｇ． 【１ ９】Ｊ，｝ｌｕＩＳＯＤ Ｄ １．ＵＩＩ（ｔＣＩ’ｓｉａ１］（１ｉｎｇ ａ．ｒ－ｇｕｎ Ｉ）ｕｌ｝ｌ｝ｌｅ ｌｌＩ ｌｌａｖｉｏｌ’…． ｛．ｅｏｐｈｙｓｉ（ ｓ，１９９４，５９｛ｌ １）：１ ７２９一ｌ ７３４． 『６１ ：炎仕．沈 非，ｒ 垠．地震数 处 巾父＿丁子波 位 特 的探泔ＩＪ］．７ １地球物瑚勘探，２００８，４３（Ｓ２）：Ｉ２１一 ｌ２４． 『２０】Ｚｉｏｌｋ，【）ｗｓｋｉ Ａ．Ｍｅｔ１＇．．；ｕｒｅｉｆｌｅＩＩＩ ｏｌ’ｄ－ｒ－ｇｔＩＩＩ Ｉ）ｕｂ１）ｈ．Ｉｌｓ‘ ｉｌｌａｌｉｏｎｓ …．Ｇｅｏｐｈｙｓｉ｛·ｓ，ｌ９９８，６３（６）：２ ００９一１ ０２４． 『７ｌ陈擘格．海 气怆 源地震予波研究进爬ｌＪＩ．中『 ｌ海Ｉ２ 汕 ，１９９４，８（１）：６０—６５． ｆ２１ｌ赵波，俞寿ＪＪＪＪ，聂勋孵，等．ｆ｝｝模拟反 ＿１＝！！方法埂』 川 …． ｛｝｝｛地球物理勘探．１９９６，３ｌ（１）：ｌ（）２—１ １６． Ｉ２２ｌ郭向宁。 必元，萤敏 ． 介卡Ｈ化ｊ 波的十ｌｊ化 钎：硬 校正ｌ Ｊ１．－￣ｉ￣１ｔｔ地球物理助探，１９９８．３３（２）：２１４－２２１． Ｉ２３１梁）匕河．地 子波捉取ｈ’法／ｉＪｌ：３：［１ｊｌ＿ｍｍ物探，ｌ９９８．３７ （１）：３１－３９． ｆ８ｌ ｆ｛景英，蒋进勇．地震子波零相位化改造的廊垌ｌＪＪ．勘探 地球物 进腱，２００３．２６（２）：ｌ２ｌ＿ｌ２３． １９俞外｝９１Ｊ』｛．离分辨牢地震勘探ｌ…．北 ： 油ｒ＿＞Ｉｋｔｌ｛版｝ｆ：， ｌ９９３：１—１８９． ｉ１０ｌ刘金俊， 修…，刘怀ｔｌ１．—】’．波零棚位化、反 积‘ｊ地 录分封辛 的关系…．海洋地质前ｆｆ｝，２０００，１６（４）：５－８． １ｌ １ ｌ畅怀存，Ｉ ，Ｅ ．海洋地震；ｔＪＪ探【１ｌ 气震源激发特性 ｆ究ＩＪ１． 物探．２００４．４３（４）：３２３—３２６． Ｉ２４ｌ杨培杰，印兴耀．地震子波提取方法综述…．ｒ ｊ１地球 物娌勘探，２００８．４３（Ｉ）：ｌ２３—１２７． Ｉ２５１李大１＿Ｉ ．坫于高阶统计的濉合＿卞Ｈ化ｒ波 汁搜 处 ｌＪ１．勘探地球物 逃展，２００５，２８（２）：１ ｌ７一ｌ２１． 『１２ｌ陈浩１水，个海燕，刘军，等．丛于近场测量的气阵列模 拟远场予波【Ｊ１． 油地球物理勘探，２００５，４０（６）：７０３— ７０７． Ｉ２６１仑阳发．牟水比， 濮．交ｆ 也震子波提取技＇￣ｌＪ１． ｌ Ｉ大学学报（『Ｉ然科 版），２００５，２９（５）：３３—３６． １２７１：ｔ－ｑ：君，１ ｆＩｆ！．海上地震资料 警车ｌｌｆ 化处州技术研究…． 石油物探，２０Ｉ２．５ｌ（４）：４０２—４０７． １ １３ｌ 成洲，陈浩＿林，牛宏年ｒ． １ 近场测量气阼列远场 ｒ。皮模拟软件研发…．物探 备，２００８，ｌ８（１）：ｌ１一ｌ７． Ｉ２８１黄缱守，仉忠良，ｌ ｊ巍．等．人然 水合物眼 缆地 系统关键技术技术’ 试验…．新能源进爬，２０１６，４（３）： ２ｌ９－２２４． Ｉ １ ４ ＩＩ川１（１Ｈ＇Ｍ．Ｍｏｄｅｌｉｎｇ ｏｆ ＧＩ ｇｕｎ ｓｉｇｎａｌｌｉｒｅｓ ｆＪ１．（；ｅｏＩ，ｌＤｓｉｃｓ Ｉ）ｒｏｓｌ｝ｅｃｌｉｎｇ．１９９２，４０（７）：７２ｌ一７４７． 『ｌ５Ｉ Ｃｈｅｒｔ Ｈ １．，Ｎｉ Ｃ ．Ｓｉｍｕｌａｔｉｏｎ ａ［１｛Ｉ ａｐｐｔｉｃａｔｉ（）ｉｉ（）ｒ ｆａｒ—ｆｉｅｌｄ Ｉ２９１何Ｉ见，张建 ｆ１．海洋地震垂ｆ（缆蚬洲系统的蚶线照ＩＪｌｊ分 析 ＩＩ．ＣＴ理论与应ｊｆｊ研充，２（）ｌ５，２４（５）：６８９—７（）０． Ｉ３０１陈金海．胤旧良．徐金祥．胀反射和ｌ海水呜 Ｉ 扰的消 除ｌＪｆ．海洋 油，２０００，３（２）：２９—３４． ｗａｖｅｌｅｔ ｆｏｒ ａｉｒ ｇｕｎ ｉｌｌ‘ｒａｙ ＩＪ １．０订（；ｅｏｐｈｙｓｉ｛‘ａｌ Ｐｒｏｓｔ）ｅ（’ ｉｇ （ｉｎ ＣＩ１ＩＩ１ＰｓＰ），２００８，４３（６）：６２３—６２５．
