第３９卷第５期 煤田地质与勘探 Ｖ＿０１．３９ＮＯ．５ ２０１１年１Ｏ月 ＣＯＡＬＧＢ０】Ｌ（１ＧＹ＆Ｅ日＇【ＤＲ 删 Ｏｃｔ．２Ｏｌ１ 文章编号：１００１—１９８６（２０１ １）０５－００６９—０３ ＶＳＰ波场分离的径向中值滤波方法 聂爱兰 伸煤科工集团西安研究院，陕西西安７１００７７） 摘要：在ＶＳＰ资料中，中值滤波法是一种简单有效的波场分离方法，主要适用于简单线性波场分 离。提出并实现的径向中值滤波方法，除具有常规中值滤波方法的优势外，对具有发散状多方向 线性同相轴轨迹的干扰波场，其波场分解可一次完成。通过井中地震（ＶＳＰ）实例测试表明，该方法 能够更有效地分离ＶＳＰ上下行波场，以及Ｐ波和Ｐ—ＳＶ波分解，且在压制面波或提取面波处理方 面具有潜在的应用前景。 关键词：中值滤波；径向扫描；波场分离 中图分类号：Ｐ６３１ 文献标识码：Ａ ＤＯＩ：１０．３９６９￣．ｉｓｓｎ．１００１．１９８６．２０１１．０５．０１７ ＶＳＰ ｗａｖｅ ｆｉｅｌｄｓ ｓｅｐａｒａｔｉｏｎ ｕｓｉｎｇ ｒａｄｉａｌ ｍｅｄｉａｎ ｆｉｌｔｅｒｉｎｇ ｍｅｔｈｏｄ ＮＩＥＡｉｌａｎ （ｘ／＇ａｎＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅ，ＣｈｉｎａＣｏａｌＴｅｃｈｎｏｌｏｇｙａｎｄＥｎｇｉｎｅｅｒｉｎｇＧｒｏｕｐ ｃ０ ，Ｘｉ＇ａｎ ７１００７７，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｓｅｐａｒａｔｉｎｇ ｕｐ ｇｏｉｎｇ ａｎｄ ｄｏｗｎ ｇｏｉｎｇ ｗａｖｅ ｆｉｅｌｄ ｅｆｉｆｃｉｅｎｔｌｙ ｉｓ ｔｈｅ ｍｏｓｔ ｆｕｎｄａｍｅｎｔａｌ ａｎｄ ｃｒｕｃｉａｌ ｐａｒｔ ｉｎ ＶＳＰ ｄａｔａ ｐｒｏｃｅｓｓｉｎｇ．Ｔｈｅ ｍｅｄｉａｎ ｆｉｌｔｅｒｉｎｇ ｉｓ ａ ｓｉｍｐｌｅ ｂｕｔ ｅｆｆｅｃｔｉｖｅ ｍｅｔｈｏｄ ｆｏｒ ｗａｖｅ ｆｉｅｌｄ ｓｅｐａｒａｔｉｏｎ，ｂｕｔ ｃｏｎｖｅｎ－ ｔｉｏｎａｌ ｍｅｄｉａｎ ｆｉｌｔｅｒｉｎｇ ｍｅｔｈｏｄ ｉｓ ｏｎｌｙ ｓｕｉｔａｂｌｅ ｆｏｒ ｓｉｍｐｌｅ ｌｉｎｅａｒ ｗａｖｅ ｆｉｅｌｄ ｓｅｐａｒａｔｉｏｎ．Ｉｎ ｔｈｉｓ ｐａｐｅｒ，ａ ｒａｄｉａｌ ｍｅｄｉａｎ ｉｆｌｔｅｒｉｎｇ ｍｅｔｈｏｄ ｉｓ ｐｒｏｐｏｓｅｄ ａｎｄ ｒｅａｌｉｚｅｄ ａｎｄ ｔｈｅ ｎｅｗ ｍｅｔｈｏｄ ａｍｅｎｄ ｔｈｅ ｏｌｄ ｏｎｅ ｗｈｅｎ ｒａｄｉａｌ ｌｉｎｅａｒ ｉｎｔｅｒｆｅｒｅｎｃｅ ｗａｖｅ ｆｉｅｌｄ ｅｘｉｓｔｓ ｉｎ ｄａｔａ，ｉｎ ｔｈｉｓ ｓｉｔｕａｔｉｏｎ，ｔｈｅ ｎｅｗ ｍｅｔｈｏｄ ｃａｎ ｓｅｐａｒａｔｅ ｗａｖｅ ｆｉｅｌｄｓ ｉｎ ｊｕｓｔ ｏｎｅ ｔｉｍｅ．ＶＳＰ ｄａｔａ ｐｒｏｃｅｓｓｉｎｇ ｔｅｓｔ ｓｈｏｗｓ ｔｈａｔ ｔｈｉｓ ｎｅｗ ｍｅｔｈｏｄ ｉｓ ｍｏｒｅ ｅｆｉｆｃｉｅｎｔ ｉｎ ｓｅｐａｒａｔｉｎｇ ｕｐ ｇｏｉｎｇ ａｎｄ ｄｏｗｎ ｇｏｉｎｇ ｗａｖｅ ｆｉｅｌｄ ｏｆ ＶＳＰ ｄａｔａ，ａｎｄ ａｌｓｏ ｇｏｏｄ ｉｎ ｓｅｐａｒａｔｉｎｇ Ｐ ａｎｄ Ｐ—ＳＶ ｗａｖｅ ｆｉｅｌｄ ｏｆ ＶＳＰ ｄａｔａ，ａｎｄ ｗｉｌｌ ｈａｖｅ ａ ｇｏｏｄ ａｐｐｌｉｃａｔｉｏｎ ｉｎ ｓｕｒｆａｃｅ ｗａｖｅ ｓｕｐｐｒｅｓｓｉｏｎ ａｎｄ ｓｕｒｆａｃｅ ｗａｖｅ ｅｘｔｒａｃｔｉｏｎ． Ｋｅｙ ｗｏｒｄｓ：ＶＳＰ；ｍｅｄｉｎａ ｉｆｌｔｅｒｉｎｇ；ｒａｄｉａｌ ｓｃａｎ；ｗａｖｅ ｆｉｅｌｄ ｓｅｐａｒａｔｉｏｎ 井中ＶＳＰ方法既能够高精度测定地层速度和获 移，得到上、下行波剖面。常规中值滤波方法主要 取精细的“时间．深度”标尺，亦能够利用ＶＳＰ反子波 适用于简单线性波场分离，本文提出的径向中值滤 算子及地层吸收等参量关系，高可信度地提高地面 波方法，除具有常规中值滤波方法的性能外，对具 地震剖面的有效分辨率（包括消除多次波）。ＶＳＰ资 有发散状多方向线性同相轴轨迹的干扰波场，其波 料与地面地震资料相比，具有反射振幅强、信噪比 场分解处理可一次性完成，并且可进一步对复杂的 高、初至波清晰可辨等优势。因此，在复杂构造地 Ｐ波和Ｐ—ＳＶ波场进行精细分解，为多波参数分析提 震精细勘探中，ＶＳＰ技术再度得到重视［１ 】。 供可靠的波场剖面信息。 在ＶＳＰ资料中高质量地分离上行波和下行波是 Ｉ方法原理 基础和关键内容。中值滤波法是一种简单有效的波 Ｉ．Ｉ 中值滤波基本算法 场分离方法，常规中值滤波法可分为５个步骤【 】： 中值滤波是基于排序统计理论的一种能有效抑 ａ．拾取各地震道下行直达波初至时间；ｂ．依据各道 制噪声的非线性信号处理技术。它把数字图像或数 初至时间做静态时移，将下行波排齐；ｃ．沿深度轴 字序列中一点的值用该点的一个邻域中各点值的中 按所选跨度（即道数）实施中值滤波运算，得到下行 值代替［引。基本算法及步骤为：ａ．设有一组数 １， 波；ｄ．原始剖面与下行波剖面相减得到上行波剖 ２，…，ｘ ）；ｂ．对这Ⅳ个数，按其数值大小排序； 面；ｅ．对上行波和下行波剖面分别做反向时 ｃ．取重排序后的中间数值为输出： 收稿日期：２０１０．１２．２４ 作者简介：聂爱兰（１９７８一），女，江西宜春人，硕士，工程师，从事煤田地震资料处理和反演工作 ・７０・ 煤田地质与勘探 第３９卷 Ｙ＝ｍｅｄｉａｎ｛ｘ１，ｘ２，…， ｝ ｆｘｉ ２ ， 为奇数 ＝ ， ｌ １［ ／２）＋ ／２＋１）］， 为偶数 中值滤波具有算法简单，去噪明显等优点： ａ．中值滤波绝对阻止噪声峰值； ｂ．中值滤波不改变阶跃函数在空间、时间上的 位置； ｃ．中值滤波对三角函数的平滑能力与滤波跨 度有关； ｄ．中值滤波对方波平滑能力与滤波跨度有关， 当跨度足够长时，方波被完全平滑。 １．２径向中值滤波法 径向中值滤波法，是基于地震资料中各波场的 视速度差异，沿强线性同相轴方向运用中值滤波法 分离出线性干扰，进而利用“减去法”衰减强线性干 扰，从而增强有效信号，得到期望的波场记录。其 实现步骤为（图１）： ， 、 ＶＳＰ道集记录读人，径向扫描斜率 参数下限 Ｐ ，上限 ，跨度参数 ２埘＋１，径向扫描斜率间隔 【 初始化道参数ｊ 。：１ 】 ～ｌ Ｉ 【初始化采样点参数＾ ＝１】 ｎ 径向扫描：记录每一斜率的振幅值序列 ａｍｐ｛／｝，及振幅和的绝对值：Ａ，，／ｔｐ） 』 求取，４ 的最大值，确定对应 的斜率≯为该采样点的线性斜率 、 Ｕ 求取斜率 对应的序歹 Ｕ＞的中值，、 记为ｒｅｅｄ：‘ ）；并得到衰减斜率 该点 的振幅值：ｔｒａ（ｉ，ｊ） ｘ（ｉｄ）－ｍｅｄ（ｉ，ｊ） 一一， Ｙ ： Ｎ 图１ ＶＳＰ波场分离的径向中值滤波法流程图 Ｆｉｇ．１ Ｔｈｅ ｆｌｏｗｃｈａｒｔ ｏｆ ｒａｄｉａｌ ｍｅｄｉａｎ ｆｉｌｔｅｒｉｎｇ ｍｅｔｈｏｄ ｆｏｒ ＶＳＰ ｗａｖｅ ｆｉｅｌｄｓ ｓｅｐａｒａｔｉｏｎ ａ对于地震道上的点（ｆ， （ｆ记录道号；Ｊ为采样 点号），在记录上选取２ｍ＋ｌ道（ｆ为２ｍ＋ｌ道的中心 记录道号），形成该采样点的径向扫描道集； ｂ．给定斜率扫描范围，在扫描道集里面逐点扫 描，当某个扫描斜率的振幅值的和的绝对值最大时， 则将该斜率定义为需分离的线性干扰的斜率； ｃ．将该斜率记录下来的振幅数列进行中值滤 波，其结果为该点的新记录； ｄ．对所有输入道集记录逐点逐道重复ｂ和ｃ过 程，分离出线性干扰记录； ｅ．用输入道集记录减去被分离出的线性干扰 记录，从而得到对线性干扰的衰减，完成径向中值 滤波。 其中径向扫描斜率范围有正负，定义沿着顺序 道时间增大时为正，反之为负。 ２模型测试 对于ＶＳＰ波场上、下行波分离，径向中值滤波 最主要的处理参数一是时空域内径向扫描斜率取值 范围；二是空间跨度（道数）。时空域内径向扫描斜 率取值范围，应视其下行波和上行波视速度值变化 情况而定；径向空间跨度对处理效果有重要影响。 通过多组模型进行空间跨度测试分析，当跨度较小 时，容易出现假频；当跨度较大时，对波场的信息 都有损失。一般道集径向空间跨度建议采用９～１３道。 图２ａ为模拟的非零偏ＶＳＰ垂直分量记录，观 测井段为２００～１ ７００ ｍ，道距１０ ｍ，偏移距３００ ｍ。 对模型运用径向中值滤波法是以分离出下行波场得 到上行波场为目标。依据原始记录中下行波场的视速 度和测试分析，径向扫描斜率取值范围为１～１０ ｎｌｓ／道， 空间跨度取值为１１道效果最佳。图２ｂ为处理结果， 下行波得到很好的分离，和图２ａ原始记录中相位和 能量保持了一致；图２ｃ的上行波场可为后续ＶＳＰ 处理提供可靠信息。 套 ： ２ （ａ） （ｂ） （ｃ） 图２ ＶＳＰ理论垂直分量记录（ａ）和径向中值滤波分离出 的下行波剖面（ｂ）及上行波场剖面（ｃ） Ｆｉｇ．２ Ｔｈｅ ＶＳＰ ｔｈｅｏｒｅｔｉｃ ｖｅｒｔｉｃａｌ ｃｏｍｐｏｎｅｎｔ ｒｅｃｏｒｄ（ａ）ａｎｄ ｉｔｓ ｒｅｓｕｌｔ ｓｅｃｔｉｏｎｓ ｂｙ ｔｈｅ ｒａｄｉａｌ ｍｅｄｉａｎ ｆｉｌｔｅｒｉｎｇ ｔｈｅ ｄｏｗｎ ｇｏｉｎｇ ｗａｖｅ ｉｆｅｌｄ ｓｅｃｔｉｏｎ（ｂ）ｎａｄ ｈｔｅ ｕｐ ｇｏｉｎｇ ｗａｖｅ ｉｆｅｌｄ ｓｅｃｔｉｏｎ（ｃ） ３实例资料处理效果测试 图３ａ为非零偏ＶＳＰ原始ｚ分量记录，观测井段 第５期 聂爱兰：ＶＳＰ波场分离的径向中值滤波方法 １ １ ｌ ２ ２ ３ ３ ３ ・７１ ・ ∞ＩＩＩ＼厦苔 Ｉ ５００～４ ５６０ ｎｌ，道距２０ ｒｆｌ，偏移距２ ０００ ｍ。由该 Ｏ Ｏ Ｏ Ｏ Ｏ Ｏ Ｏ Ｏ 枷咖瑚枷咖伽渤 为１ｌ道，上行Ｐ波和上行Ｐ—ＳＶ波的分离结果。 记录剖面可见，除常规目标信号Ｐ波（下行直达Ｐ波、 上行反射Ｐ波）之外，与地层岩性密切相关的转换横 波Ｐ—ＳＶ在此记录剖面中也具有可识别的品质特征 （见图中标注）。下面以该实测ＶＳＰ记录为例，就径 向中值滤波方法用于ＶＳＰ记录上、下行波分离，以 及Ｐ波和Ｐ—ｓｖ波分解处理步骤和测试结果作如下 简述。 （ａ） （ｂ） （ｃ） 图３ ＶＳＰ实测垂直分量记录（ａ）和径向中值滤波方法分 离出的下行波剖面（ｂ）及结果上行波场剖面（ｃ） Ｆｉｇ．３ Ｔｈｅ ＶＳＰ ａｃｔｕａｌ ｖｅｒｔｉｃａｌ ｃｏｍｐｏｎｅｎｔ ｒｅｃｏｒｄ（ａ）ａｎｄ ｉｔｓ ｒｅｓｕｌｔ ｓｅｃｔｉｏｎｓ ｂｙ ｔｈｅ ｒａｄｉａｌ ｍｅｄｉａｎ ｆｉｌｔｅｒｉｎｇ ｔｈｅ ｄｏｗｎ ｇｏｉｎｇ ｗａｖｅ ｉｆｅｌｄ ｓｅｃｔｉｏｎ（ｂ）ａｎｄ ｈｔｅ ｕｐ ｇｏｉｎｇ ｗａｖｅ ｉｆｅｌｄ ｓｅｃｔｉｏｎ（ｃ） Ａ— 行Ｐ波；Ｂ＿一匕行Ｐ波；Ｄ—下行Ｐ—ＳＶ波；Ｉ）＿—士行Ｐ．ＳＶ波 ３．１上、下行波场分离 图３ａ中视下行波场为干扰，需被分离，径向扫 描斜率范围取２～１５ ｍｓ／Ｊ￣，空间跨度为１１道。图３ｂ 为对应的下行波分离结果。图３ｃ是图３ａ记录“减去” ３ｂ记录所得的上行波记录剖面。 由图３处理结果可见，图３ｂ中几乎没有上行波 成份，下行Ｐ波和Ｐ．ＳＶ波能量在记录剖面上得到 了很好的增强，并且与图３ａ相对比，两者的下行Ｐ 波初至和下行Ｐ—ＳＶ波的相位特征具有高保真度及 较好的一致性，这一特点相对其他二维滤波方法有 一定优势（如Ｆ—Ｋ方法）。 ３．２ Ｐ波和Ｐ—ＳＶ波分离 图３ｂ中下行Ｐ波和下行Ｐ—ＳＶ波，特别是图 ３ｃ中的上行Ｐ波和上行Ｐ—ＳＶ波相互重叠，即便是 不提取转换波信息（Ｐ—ｓｖ）而只考虑纵波（Ｐ）问题时， 对下行波子波提取、上行波成像处理及解释都将造 成一定的影响。因此，将Ｐ波和Ｐ．ＳＶ波进一步分 解是精细波场分离的目标之一。对于上、下行波场 中Ｐ波和Ｐ—ＳＶ波分离，可依据图３ａ所示ｚ分量记 录中纵波和转换横波的上、下行波视速度（或时差） 的取值区间（差异）进行分离。图４是图３ｂ对应的径 向扫描斜率范围取８～１５ ｍｓ／道、空间跨度为１１道， 下行Ｐ波和下行Ｐ．ＳＶ波的分离结果。图５是图３ｃ 对应的径向扫描斜率范围取８—１５ ｍｓ／道、空间跨度 由图４可见，下行Ｐ波初至和下行Ｐ—ＳＶ波能 量得到了增强，并且两者相对比，下行Ｐ波初至和 下行Ｐ—ＳＶ波的相位特征具有高度一致性。图４ｂ中 下行Ｐ波初至被有效分离，下行Ｐ—ＳＶ波与下行Ｐ 波初至的交点更清晰。图５ａ中，浅层三组上行Ｐ波 同相轴波组和图３ｂ中的相位一致，并且能量得到了 增强；图５ｂ中上行Ｐ波被分离出去，上行Ｐ—ｓｖ波 能量得到增强，极易识别。分离出来的Ｐ波和Ｐ—ＳＶ 波场记录都可以提取各自的波场信息，这为后续的 ＶＳＰ处理提供可靠资料。 １ ０００ ｌ ４００ １ ８００ 兽２ ２００ 要２ ６００ １ １ １ 兽＼厘鲁 ２ Ｏ ４ ８ ２ ３ ０００ ∞ ∞ ∞ ∞ ∞ ∞ ３ ４００ ３ ８００ 图４对图３ｂ记录用径向中值滤波分离出的下行Ｐ 波（ａ）和下行Ｐ．ＳＶ波（ｂ）剖面 Ｆｉｇ．４ Ｔｈｅ ｄｏｗｎ ｇｏｉｎｇ Ｐ ｗａｖｅ ｆｉｅｌｄ ｓｅｃｔｉｏｎ（ａ）ａｎｄ ｔｈｅ ｄｏｗｎ ｇｏｉｎｇ Ｐ—ＳＶ ｗａｖｅ ｉｆｅｌｄ ｓｅｃｔｉｏｎ（ｂ）ｗｉｔｈ ｔｈｅ ｒａｄｉａｌ ｍｅｄｉａｎ ｉｆｌ— ｔｅｒｉｎｇ ｍｅｔｈｏｄ ｆｏｒ ｔｈｅ ｄｏｗｎ ｇｏｉｎｇ ｗａｖｅ ｆｉｅｌｄ ｓｅｃｔｉｏｎ ３４００ ３ ８ｏｏ （ａ） （ｂ） 图５对图３ｃ记录用径向中值滤波分离出的上行Ｐ 波（ａ）和上行Ｐ－ＳＶ波（ｂ）剖面 Ｆｉｇ．５ Ｔｈｅ ｕｐ ｇｏｉｎｇ Ｐ ｗａｖｅ ｉｆｅｌｄ ｓｅｃｔｉｏｎ（ａ）ａｎｄ ｔｈｅ ｕｐ ｇｏｉｎｇ Ｐ－ＳＶ ｗａｖｅ ｉｆｅｌｄ ｓｅｃｔｉｏｎ（ｂ）ｗｉｔｈ ｔｈｅ ｒａｄｉａｌ ｍｅｄｉａｎ ｉｆｌｔｅｒｉｎｇ ｍｅｔｈｏｄｆｏｒｔｈｅｕｐｇｏｉｎｇｗａｖｅｆｉｅｌｄ ｓｅｃｔｉｏｎ ４结语 理论和实际资料处理结果表明，径向中值滤波 法除了能保持中值滤波特性之外，还具有如下特点： （下转第７５页） ２ ３ ６ ０ 第５期 王喜迁等：高密度电法在岩溶探测中的应用 ・７５－ 器Ｆ Ｄ２ 。篚 图９成果解释图 Ｆｉｇ．９ Ｉｎｔｅｒｐｒｅｔａｔｉｏｎ ｒｅｓｕｌｔｓ 粘土充填区溶蚀区 ９５９ ９６９ ９７９ ９８９ ９９９ ｌ ００９ １ ０ｌ９ １ ０２９ ｌ ０３９ｌ ０４９ ｌ ０５９ ｌ ０６９ ｌ ０７９ １ ０８９ ｌ ０９９ Ｋ１２１＋１０７ ｂ．Ｋ１２１＋０７６一＋０８１段浅部（Ⅲ）存在一富水并 有粘土充填的松散区，实地观察发现此区以前是一 旧水沟。 法断面结果反映的是灰岩与上覆第四系地层的平均 电性反映，当灰岩表面溶蚀发育而呈锯齿状时，断 面测量资料能反映灰岩面总体起伏状况。地形起伏 地段的断面实测结果，应进行地形校正，以消除地 形影响，使测量断面与实际断面更接近。－ 参考文献 ［１】傅良魁．电法勘探教程【Ｍ】．北京：地质出版社，１９８３． ｃ．Ｋ１２１＋０８９一十０９３段（Ⅳ）存在一个充水的溶蚀 区，顶部深约５ ｍ。 ｄ．Ｋ１２０＋９９４段公路中线向右９ ＩＴＩ处浅部有一小 的粘土充填区；Ｋｌ２１＋ｏ０７—０１１段公路中线向右９ ｍ 处浅部有一小的粘土充填区；Ｋ１２Ｏ＋９７１一Ｋｌ２Ｏ＋９７２ 公路中线向右３ ｍ处浅部有一小的粘土充填区。 ４结论 ［２】何伟，于鹏，张罗磊，等．高密度电法在探测地下金属与非金 属组合管线中的应用［Ｊ】．工程地球物理学报，２００８，５（ｉ）： ９５—９８． 高密度电法兼具剖面法和测深的功能，具有点 距小、数据采集密度大的特点，能较直观、形象地 反映断面电性异常体的形态、产状等。在灰岩分布 区应用高密度电法探测第四系土洞发育、岩溶、断 裂发育情况等，能取得较好的地质效果。高密度电 （上揍第７１页） ［３顾中华．高密度电法探测溶洞成果精度分析【３］Ｊ】．勘察科学科 技，２００８（４）：２７－２９． 【４】李成香，王建军，唐诗群．应用综合物探方法探测隐伏岩溶漏 斗［Ｊ】．工程勘察，２００８（３）：９９－１０１． ［５王建军，强建科，李成香，等．高密度电法在地面塌陷勘察中 ５］的应用【Ｊ】．资源环境与工程，２００８，２２（特刊）：２３２－２３４． ａ．径向中值滤波不需拾取地震道下行直达波初至 解，为多波参数分析提供可靠的波场剖面信息；对 叠前地面地震反射资料，可分离任意方向的线性干 扰，信噪分离能力强、振幅保真度高，是反射资料高 分辨处理的理想方法。 时间，也不需依据各道初至做静态时移将下行波排齐。 ｂ．常规中值滤波法是沿深度轴实施中值滤波， 是一维滤波运算；径向中值滤波是在道方向以及时 间方向上进行中值滤波，是二维滤波运算。 ｃ．常规中值滤波法最基础的工作是在拾取地 参考文献 ［１】朱光明．垂直地震剖面方法【Ｍ】．北京：石油工业出版社， １９８８． 震道下行直达波初至时间之后进行的，对资料的信 噪比要求高；而径向中值滤波法只要地震资料线性 干扰的视速度可识别就可以实施，相较而言，对地 【２】郭建．ＶＳＰ技术应用现状及发展趋势［Ｊ】．勘探地球物理进展， ２００４，２７（１）：１＿Ｉ８． 震资料的信噪比要求就低些，比常规中值滤波法有 更广的适用性。 【３】涂齐崔．ＶＳＰ射线追踪模拟与处理方法研究及应用【Ｄ】．青岛： 中国海洋大学，２００６． ｄ．常规中值滤波法在需要消除某组下行波时， 先将该组下行波排齐后再进行中值滤波，要消除多 组下行波就需进行多次中值滤波运算；而径向中值滤 波只需一次即可完成不同视速度的多组下行波分离。 ｅ．反射地震资料中，面波一般都视为干扰，在 处理时都对其进行压制；而径向中值滤波法可以将 地震资料中的面波分离出来，分离出来的波场有利 于研究地表层低降速带结构［９－１０ｌ。 ［４】于海超．非线性滤波技术在地震信号处理中的应用【Ｄ】．大 庆：大庆石油学院，２００７． ［５】唐金良，曹辉．中值滤波在井间地震资料处理中的应用【Ｊ】．石 油物探，２０ｏ５，４４（１）：４７－５０． 【６】夏洪瑞，董江伟．环形非对称中值滤波技术在可控震源资料 去噪处理中的应用【Ｊ】．石油物探，２００９，４８（５）：５１０－５１４． ［７】李承楚．关于中值滤波的理论基础【Ｊ】．石油地球物理勘探， １９８６，２１（４）：３７２－３７９． 【８】胡光瑞，郑志航，戚飞虎．二维数字信号处理Ⅱ【Ｍ】．北京： 科学出版社。１９８５． 径向中值滤波法是一种非线性滤波技术，不仅 对ＶＳＰ资料，对上下行波场分解处理可一次性完成， 亦可进一步对复杂的Ｐ波和Ｐ—ＳＶ波场进行精细分 ［９］周熙蓑，王利国．利用瑞利面波调查表层结构的可行性探讨 ［Ｊ］．石油地球物理勘探，２００４，３９（２）：１８１—１８６． ［１Ｏ】成锁．近地表调查中的面波方法与利用井资料提高分辨率的 技术［Ｄ】．长春：吉林大学，２００８．