基于互联网地震短信发布平台的构建

基于互联网地震短信发布平台的构建 地震层位自动拾取方法 3冯智慧 刘财冯王典王兆国 ( )吉林大学地球探测科学与技术学院 ,吉林长春 130026 冯智慧 , 刘财 , 冯 , 王典 , 王兆国 . 基于互四阶累积量一维切片的地震层位自动拾取方法 . 石油地球物理勘 探 , 2011 , 46( 1) :58,63 摘要 为了解决强相干高斯色噪声干扰下地震层位无法正确拾取的问题 ,文中提出了一种基于互四阶累积量 一维切片的地震层位自动拾取方法 ,阐述了该方法的基本原理和具体算法 。该方法利用地震层位上横向波形 的相似性以及互四阶累积量一维切片对高斯噪声不敏感的特点 ,实现了同相轴的时间延迟估计 ; 同时对直接影 响地震层位拾取精度的相关参数进行了设定 。理论模型计算明 ,该方法对相干高斯色噪声有很强的抑制能 力 ,信噪比达到 - 10dB 。实际资料应用结果与原始数据的地质层位吻合较好 ,断层 、断点位置准确 ,其拾取精度 不亚于人工方法的精度 ,而效率远远高于人工方法 。 信噪比 相干高斯色噪声关键词 互四阶累积量 一维对角切片 层位自动拾取 中图分类号 : P631 文献标识码 : A 切片值 ,并根据一维切片值的变化情况来检测地震 信号 。当两道间有效地震信号重合时 ,由柯西 —施 1 引言 瓦茨公式可知 ,它们的互四阶累积量一维切片值达 如今 ,地震勘探采集的数据量日趋增大 ,传统的 到最大值 ,从而确定地震层位位置 。本文提出的方 人工层位拾取方法不仅效率较低 ,而且容易受到人 法不受两道间高斯噪声是否是独立关系的限制 ,克 为因素及噪声的影响 。为了克服传统人工拾取的弱 服了传统互相关方法由相干高斯色噪声带来的诸多 点 ,许多地震层位的自动提取方法应运而生 ,如神经 问题 。理论模型计算及实际地震数据的处理结果表 [ 1 ,2 ] [ 3 ] [ 4 ] 明 ,该方法能有效抑制相干高斯色噪声 ,是一种有效 网络方法、倾角扇形方法、模式识别方法、链[ 5 ] [ 6,8 ] 匹配算法、互相关方法等 。在上述方法中应 的地震层位拾取方法 。 用最广泛的是互相关方法 。但是此法存在一个致命 性缺点 ,即当干扰噪声为强相干噪声时 ,对两道地震 2 层位拾取的高阶累积量理论基础 信号进行互相关运算 ,噪声间的互相关值要大于有 效地震信号间的互相关值 ,致使拾取的有效信号间 2 . 1 地震信号高阶累积量理论基础 的延迟时出现误差 。 通常假设地震为一零均值平稳随机过程本文提出的基于互四阶累积量一维切片的地震 [ 9 ] ( ) x n, 可用如下褶积模型描述一道地震记录层位自动拾取方法 ,主要是利用地震层位上横向波 q 形的相似性以及互四阶累积量一维切片对高斯噪声 ( )( )( )( ) ( ) ( )x n = s n + v n = h ir n - i+ v n? i = 0 不敏感的特点进 行地 震 层位 拾取 。首先 通过 精 确 ( ) 1 “层位标定”确定参考道和层位追踪的起点 ,再从参( ) ( ) ( ) 式中 : s n为有效地震信号 ; h n为地震子波 ; r n考道中截取一段地震信号 ,计算这部分地震信号与 ( ) 为反射系数 ; v n为加性不噪声 。 其他道上滑动时窗中地震信号的互四阶累积量一维 当应用高阶累积量理论对地震记录进行处理 3 吉林省长春市民主大街 938 号吉林大学地球探测科学与技术学院 ,130026 。Email :f englezhi1981 @163 . co m 本文于 2010 年 1 月 26 日收到 ,最终修改稿于同年 11 月 10 日收到 。 [ 10 ] 时 , 一般应有如下假设: ( ) ( ) 1有效信号 s n一般可假设为均值为零 、独 3 互四阶累积量一维切片方法原理 立同分布的非高斯过程 , 其相应的 k 阶高阶累积量 值 c 存在 ;3 . 1 互四阶累积量一维切片方法原理ks ( ) ( ) ( ) 2v n为与 s n统计独立的零均值 、因果的 设地震剖面中地震记录 A 和 B 的波形采样信 ( ) ( ) 加性高斯色噪声 。 号为 x n和 y n, 其中 A 为参考道 , B 为要追踪的 ( ) ( ) ( ) ( ) 对于 k 道地震记录 x ni = 1 , 2 ,, k, x n地震道 。根 据地 震 信号 的传 播性 质 , 假 设n和x i i [ 11 ]( ) y n满足如下的信号模型 的 k 阶联合累积量可以表示为 (τττ) ( ) C,, ,= c um[ xn,( )( ) ( ) x, x, , x 1 2 k - 1 k 1 x n= s n+ vn( )1 6 1 2 k ( τ) ( τ) ( ) 2 x2 n +1 ,, x k n +k - 1 ]( )( ) ( ) ( )7 y n= s n - D+ vn2 ( ) ( ) τττ: s n为 A 道的有效地震信号 ;vn为该道的其中 式中 :,,,表 示延 迟 时 间 ; 等 式 右 侧 符 号 1 1 2 k - 1 ( ) ( ) 观测随机高斯色噪声 , 独立于 s n; s n - D为 B“c um ”表示做“ k 阶累积量”运算 , k ?3 时称为互高 k ( ) ( ) 道的有效地震信号 ; D 为 y n相对于 x n真实的 阶累积量 。 零均值地震记录的互二阶 、互三阶和互 ( ) 时间延迟 ; vn为该道的观测随机高斯色噪声 , 独四阶累 2 [ 12 , 13 ]( ) 积量分别为 立于 s n - D。根据互高阶累积量一维切片的性质 (τ) ( ) ( τ) ( ) E[ x nx n +]3 1 Ck1 k2 1 = x , x 可知 :含有高斯噪声的地震信号的互高阶累积量一 k1 k2 (ττ)( ) ( τ) C 1 ,2 = E[ x k1 nx k2 n +1 ×x , x , x维切片 与 有 效 信 号 的 互 高 阶 累 积 量 一 维 切 片 相 k1 k2 k3 [ 16 , 17 ] ( τ) x n +]( )k3 2 4 等。 [ 18 ] (τττ) C,,x , x , x x 1 2 3 Tugnait 提出了基于互四阶累积量一维切片 k1 k2 k3 k4 ( ) ( τ) ( τ) ( τ) E[ x k1 nx k2 n +1 x k3 n +2 x k4 n +3 ]= 的时延估计函数表达式 。式中采用了可变长度 ( ) ( τ) ( τ) ( τ) - E[ x nx n +] E[ x n +x n +]k1 k2 1 k3 2 k4 3 时窗 。此法有一个明显缺点 ,即在实际应用中存在 ( ) ( τ) ( τ) ( τ) - E[ x nx n +] E[ x n +x n +]k1 k3 2 k2 1 k4 3 受坏道影响大 、计算结果失稳 、时移量拾取不足的问 ( ) ( τ) ( τ) ( τ) E[ x nx n +] E[ x n +x n +]k1 k4 3 k2 1 k3 2 - 题 。为了避免上述问题 ,本文采用稳健性更好的固 ( )5 定长度时窗作为估计函数的分析时窗 ; 并且根据实 τττ 式中 :,,表示 3 道不同地震记录的延迟时间 ;1 2 3 际需要对延迟时参数的位置作了相应的改变 。 (τ) k, k, k, k?{ i = 1 , 2 ,C1 ,, k } ; 互四阶累积量一维切片时延估计函数的离散化 1 2 3 4 x , xk1 k2 (ττ) (τττ) C,和 C,,分别表示x , x , x 1 2 x , x , x , x 1 2 3 形式为 k1 k2 k3 k1 k2 k3 k4 ( ) τ) τ) ( τ) (( ) ( ) ( J = | cu m[ x n, x n, y n +, y n +]|多个地震记录 x i n的互二阶 、互三阶和互四阶累 4 积量 。( ) 8 2 . 2 高阶累积量一维切片τ其中变量表示时间延迟的估计值 。由柯西 —施瓦 [ 20 ] ( ) 由式 5可知 , 互四阶累积量是一个有三个参数 (τ) τ茨公式知 , 当= D 时 , J 有最大值 。 变量的四维函数 , 在有限数据情况下相对于互相关 ( ) ( ) ( τ) ( τ) cum[ x n, x n, y n +, y n +]的离散化 4 计 运算其估计方差较大 。当互四阶累积量只有一个参算公式为 ( τττ) τ数变量时 如= 0 ,==, 称为互四阶累积量 1 2 3 一维τ) ( τ) ( ) ( ) ( c um[ x n, x n, y n +, y n +]4 [ 14 ] N N 切片。互四阶累积量的一维切片降低了估 计方2 1 1 2 2 ()τ( ) ( ) ( ) τ x ny n + x nyn +- 2 = - ?? 差 , 而且利于实际应用中的图形显示 。互高阶 累积 N N n = 1 n = 1 [ 15 ] N N 量一维切片具有以下性质: 1 2 2 1 ((τ) x n) y n +( )9 ?? ( ) 1互高阶累积量一维切片仍然包含了地震信 N N n = 1 n = 1 号的互高阶累积量信息 ; 式中 : N 为 做 累 积 量 计 算 的 分 析 时 窗 长 度 ; 变 量 τ( ) 2零均值的加性高斯噪声的一维切片仍然为 表示时间延迟的估计值 。 零 , 即含有零均值加性高斯噪声的地震记录的一维 3 . 2 地震层位拾取过程的相关参数设定 切片值等于有效地震信号的一维切片值 。 在进行地震层位自动拾取前需要指定相关的参 () 片方法则能抑制相干噪声的影响 。始位置 种子点、参考道 、倾角时窗范围 。 () 1在实际应用中可以通过层位标定 ,确定参考 为了验证互四阶累积量一维切片方法比互相关 ( ) 道上的地震层位拾取起始位置 种子点,而且层位 ,首先建 方法对相干高斯色噪声有更强的抑制能力 标定的准确性决定了层位拾取的准确性 。 立理论模型 。本文采用双层倾斜地质模型进行定量 () 计算 ,模型的参数由表 1 给出 ,选取主频为 55 Hz 的 2参考道的选择方式有三种 : ?设定一个固 () 定的参考道 ,其余各道与该道进行互四阶累积量运 Ricker 子波制作合成地震记录 图 2a,其中不考虑 能量损失及 A V O 效应 。 算 ,估计出时间延迟值 ; ?以前一道作为参考道 ,进 行邻近道间的时间延迟值估计运算 ; ?在层位拾取 表 1 倾斜地质模型参数 之前 , 首 先 对 地 震 记 录 进 行 平 均 分 组 , 每 组 含 有 埋深速度密度地层倾角 层位 3 3,5道地震记录 ,每组最后一道作为下一组的参考 m m/ s g/ cm 采样点/ 道 道 。实际应用中 ,可根据具体情况选择上述一种参 1 300 2800 1 5 考道选择方式 。 2 3500 — [ 20 ] () 3倾角扫描的范围一般根据实际地层的倾 在理论合成地震记录基础上加入相干高斯色噪角来估计 。 ( ) ( ) 声 ,信 噪 比 分 别 为 0dB 图 2 b 、- 5 dB 图 2c 、 () - 10 dB 图 2 d。由于噪声间的相关性 ,在地震剖面 中出现了许多“假同相轴”,随着信噪比的降低 ,它们 对地震层位拾取的干扰增强 。 ( ) 我们 从 参 考 道 中 截 取 的 一 段 地 震 信 号 x n { n = 1 , 2 , , N} , N = 141 , 分别计算这段信号与其他 道上相应滑动时窗内信号的互四阶累积量一维切片 图 1 地震层位自动拾取示意图 值和互相关值 。图 3 、图 4 、图 5 分别为三种信噪比情 3 . 3 互四阶累积量一维切片方法实现过程 况下两种方法的地震层位拾取结果与误差情况 。在 互四阶累积量一维切片方法的实现过程如下 : 信噪比为 0dB 时 ,两种方法的拾取结果均未出现误 () 1首先通过层位标定确定种子点 ,从参考道 A ( ) 差 ,都能准确地拾取层位 图 3。但 信噪 比 降低 到 ( ) 上截取长度为 N 的地震信号 x n{ 1 , 2 , , N } , N - 5dB时 ,由 图 4 可 见 : 互 四 阶 累 积 量 法 在 第 46 个 的取值应该覆盖 A 道上有效信号的波形 , 参考道的 CDP 点出现 6 个采样点的拾取误差 ;互相关方法分别 选择可以根据具体情况选取三种方法中的一种 ;在第 40 、43 、46 个 CDP 点出现 30 个点 、6 个点 、6 个点 ( ) 2对于待测地震道 B , 用时窗长度为 N 的数的误差 。当信噪比降低到 - 10dB 时 ,互四阶累积量 ( τ) 在实际地层的倾角范, N } , 据 y n +{ n = 1 , 2 , 法除了少数几个 CDP 点出现很小的拾取误差外 ,其 (τ) 围 [ 0 , D]内移动 , 计算 J ;max 他各个 CDP 点均拾取正确 ,而互相关方法的拾取结 ( ) (τ) τ(τ) 3对于 J ,= 0 , 1 ,, D, 将 J 中出现 max 果误差很大 ,已经不能有效抑制相干噪声 。 的最大值作为检测道待测地震信号的标志 , 并把出 下面对上述结果做进一步分析 。以 拾取 图 2c 现最大值的点作为同相轴上的点 。中第 40 个 CD P 的地震层位为例 ,图 6a 为互相关方 法对应的互相关曲线图 ,图 6 b 为新方法对应的互四 阶累积量一维切片曲线图 。从图 6a 可知 ,由于相干4 理论模型计算及结果分析 高斯色噪声的干扰 ,互相关值随检测时窗的滑动具 由累积量性质可知 ,互相关方法可以很好地抑有波浪式的变化 , 在 440 m s 处出现极值 , 但正确层 位是在 410 ms 处 。这说明由于互相关方法不能有效 制独立噪声的干扰 ,在实际应用中也证明了这一点 。 抑制强相干高斯色噪声的干扰 , 噪声间的延迟时作 由文献 [ 21 ]可知 ,当地震记录中存在大量相关噪声 图 2 理论地震记录 ( ) ( ) () ( ) a未含噪声 ; b含 0dB 相干噪声 ; c含 - 5dB 相干噪声 ; d含 - 10dB 相干噪声 () ( )图 5 SN R = - 10dB 时层位拾取结果 a及其误差 b 为有效信号间的延迟时被拾取 , 致使出现 30 ms 的 层位拾取误差 。同样是对图 2c 中第 40 个 CD P 的 地震层位进行拾取 ,从图 6 b 可见 ,互四阶累积量一 维切片随检测时窗滑动具有近似单调的变化 ,且在 ( ) 极大值处 410 m s正确拾取了同相轴上的点 。这说明新方法 对 相 干 高 斯 色 噪 声 有 更 强 的 抑 制 能 力 。 此外 ,与一般统计量存 在估 计方 差 一样 , 在 有限 数 据的情况下 ,相干高斯色噪声的互四阶累积量一维 切片并非完全为零 ,这是在信噪比很低时该方法出 现较小拾取误差的原 因 。在模 型计 算中 我 们还 对 信噪比低于 - 1 0 dB 的 情 况 进 行 了 计 算 , 但 两 种 方 () ( ) 图 6 互相关值 a与互四阶累积量一维切片值 b随时窗滑动的变化曲线 R 和 J 分别为互相关值和互四阶累积量一维切片值 , T无用数据对处理结果的干扰 。从图中可以看到 法的拾取结果 均 出 现 了 较 大 误 差 。因 此 , 为 了 有1 效地利用互四阶累积 量 一维 切片 方法 , 应保 证 观 分别在第 10 、第 26 、第 47 个 CD P 处出现断层 ,而且 测信号能够达到一定 的 信噪 比 , 所 以当 信噪 比 很 整个层位横向上波形存在一定的变化 ; T在第 472 低时 , 可提前采用滤波 等预 处 理方 法提 高待 测 信个 CD P 出现断距较大的断层 , 横向波形 变 化相 对 号的信噪比 。T较小 。两 个 目 的 层 位 附 近 均 出 现 干 扰 层 , 尤 其1 T附近的干扰层对层位追踪的影响较大 。图 7 b 为 1 应用互四阶累积量一维切片方法的地震层位自动拾 5 实际数据层位参数提取 取结果和原始剖面的叠合图 ,如图所示 ,拾取的同相 轴光滑连续 ,与实际的地质层位一致 ,而且对断层断 图 7a 是某一实际地震剖面 ,图中同相轴 T1 、T2 为所要追踪的目 的层 位 。在 进 行地 震数 据处 理 之 点的拾取也比较准确 ,并没有出现拾取结果远离目 前 ,可以通过对目的层 T、T的自动拾取 , 把感兴 1 2 的层位的现象 ,其精度不亚于人工拾取方法精度 ,而 趣的地震数据提取出来 ,既减少了计算量 ,又避免了 且在效率上要远远高于后者 。 图 7 地震层位拾取剖面对比图 ( ) ( ) a实际地震剖面图 ; b层位拾取结果与原始剖面的叠合图 际地震信号分析实例说明该拾取方法有效 ,利用互 四阶累积量一维切片函数描述不同地震道间信号的 6 结束语 相关性 ,可以克服传统互相关方法由相干高斯色噪 声带来的诸多问题 。相对于已广泛应用的互相关方 本文利用地震层位横向波形的相似性以及高阶 法 ,本文提出的互四阶累积量一维切片方法是一种 累积量对高斯噪声不敏感的特点 ,给出了基于互四 阶累积量一维切片自动拾取地震层位的方法 。由数 初步的尝试 , 还需要更加深入研究 , 从而进一步提 ( 值模拟结果可知 , 当信噪比很低时 即使信噪比达 高高阶累积量理论 在地 震 信号 分析 和处 理 中的 实 ) 到 - 10dB,该方 法仍 能 有效 地拾 取 地震 层 位 。实 用性 。 [ 10 ] 戴永寿 ,王俊岭 ,王伟伟等 . 基于高阶累积量 A RMA 参 考 文 献 模型线性非线性结合的地震子波提取方法研 究 . 地 () 球物理学报 ,2008 ,51 6:1851,1859 [ 1 ] Fra nklin Kemp ,赵蓓蓓译 . 应用神经网络分支和边界 Dai Yo ng2sho u , Wa ng J un2ling , Wa ng Wei2wei et al . 值搜索 法 的 地 震 层 位 追 踪 系 统 . 地 质 科 学 译 丛 , Sei smic wavelet ext ractio n via cumulant2ba sed A RMA () 1993 ,10 4:93,98 mo del app roach wit h linea r and no nlinear co mbina2 Fra nklin Kemp , Zhao Pei2p ei t ra nslated. The system ( ) tio n. Chi nese J ou rn al o f Geo p h y s ics , 2008 , 51 6 : of sei smic ho rizo n t racki ng by neural net wo r k s branch 1851,1859 and bo undar y val ue sea rching. J ou rn al o f Geoscie nce [ 11 ] 邱天爽 ,张旭秀 ,李小兵等 . 统计信号处理 . 北京 : 电 ( ) T ra ns l at i ons , 1993 , 10 4:93,98 子工业出版社 ,2004 [ 2 ] H ua ng Ko u2 Yua n , 陈 炳 文 译 . 用 于 地 震 层 位 拾 取 的 [ 12 ] 尹成 ,熊晓军 ,张白林等 . 四阶累积量在剩余静校正中 Hopfield 神经网络 . 美国 勘探地 球 物 理 学 家 学 会 第 () 的应用研究 . 天然气工业 ,2004 ,24 12:48,50 67 届年会集 ,1997 Yin Cheng , Xio ng Xiao2j un , Zhang Bai2li n et al . Appli2 Kou2 Yuan Huang ,Chen Bing2wen t ranslated. Hop f iel d catio n of fo urt h2o rder cumula nt in re sidual static co r2 () rectio n. N at u ral Gas I n d ust r y , 2004 , 24 12:48,50 熊neu ral net w orks to be use d i n e x t ract i n g seismic hori z on. 晓军 ,尹成 ,张白林等 . 基于四阶互累积量的小断层 自Cor p us of S EG 67 t h annual meeti n g , 1997 [ 13 ] ) ( 动识别方法 . 石油地球物理勘探 ,2004 ,39 增刊: [ 3 ] 孙奇政 . 基于人工神经网络方法的断层两侧地层层位 158,164 追踪的研究 [ 博士学位论文 ] . 中南大学 ,2007 ,22,25 Xio ng Xiao2j un , Yin Cheng , Zhang Bai2lin et al . A Sun Qi2zheng. T he l a y e r t rack i n g o f f a ul t ’s s i des met ho d of a uto matic small f ault s identificatio n ba sed base d on t he a rt i f ici al ne u ral net w o rks [ D ] . Cent ral ( ) on cro ss2fourt h2o rder cumulant s. O G P , 2004 , 39 Sup: So ut h U niver sit y ,2007 ,22,25 158,164 [ 4 ] 董恩清 ,朱光明 ,吴文奎. 应用模式识别自动追踪地震 赵清明 ,陈西宏 ,张敏. 基于四阶累积量一维切片谱线 () 剖面同相轴. 西安石油学院学报 ,1998 ,13 2:12,16 [ 14 ] () 增强方法. 弹箭与制导学报 ,2004 ,24 3:236,238 Do ng En2qing , Zhu Guang2ming , Wu Wen2kui . Appli2 Zhao Qing2ming , Chen Xi2ho ng , Zha ng Min. Ba sed catio n of pat ter n reco gnitio n to a uto matically t racking o ne2dimensio nal slice of t he ho s va ria ble step size a2 sei smic p rofile event s. J ou rnal o f X i ’a n S hi y ou U ni 2 ( ) () ve rs i t y Nat ural Science Edition,1998 ,13 2:12,16 dap tive algo rit hm fo r line enhancement . J ou rn al o f 许建华 ,王世库 ,赵廷寿 . 用链匹配算法自动拾取同 P roj ect i les R ockets M iss i les an d Gui d ance , 2004 , [ 5 ] () 相轴 . 石油物探 ,1990 ,29 2:13,23 ( ) 24 3:236,238 Xu J ia n2hua , Wa ng Shi2ku , Zhao Ti ng2sho u. Picking 柯宏发 ,陈永光 ,张哲峰 . 基于四阶累积量对角切片 [ 15 ] up t he event s a uto matically by chain matching algo2 () 的短波自适应通信信号检测 . 电子学报 ,2006 ,34 3: rit hm. Geo p h ys ic al P ros pect i n g f or Pet role u m , 419,423 ( ) 1990 , 29 2:13,23 Ke Ho ng2f a , Chen Yo ng2gua ng , Zhang Zhe2f eng. H F 王一重 ,冯 ,王俊祥等 . 相干技术在构造形态自动识 adap tive radio co mmunicatio n signal detectio n ba sed [ 6 ] ( ) 别中的应 用 . 吉 林 大 学 学 报 地 球 科 学 版 , 2008 , 38 o n fo urt h2o r der cumulant s diago nal slice . A ct a Elec2 () 增刊:85,89 ( ) t ronic a S i nic a , 2006 , 34 3:419,423 Wa ng Yi2zho ng , Feng Xua n , Wa ng J un2xiang et al . Mendel J M . Tuto rial o n higher2o r der stati stic s in sig2 [ 16 ] A n applicatio n of co herence technique in st r uct ure nal p rocessing and system t heo r y. T heo ret ic al Res ul ts ( ) an d S ome A p p l ic at i ons P roc I E E E , 79 3:278,305 mo dalit y a uto matic reco gnitio n. J ou rn al o f J i l i n U ni2 Chr yso sto mo s L Nikia s , J er r y M Mendel . Signal ( ) ( ) ve rs i t y Eart h Science Editio n , 2008 , 38 Sup: [ 17 ] p rocessing wit h higher2o r der sp ect ra . I E E E T ra ns 85,89 ( ) S i g n al P rocessi n g , 1993 , 10 3:10,37 [ 7 ] 陈洁 . 超宽带雷达信号处理及成像方法研究 [ 博士学 Tugnait J K. Time delay esti matio n wit h unkno w n 位论文 ] . 中国科学院研究生院 ,2007 ,26,39 [ 18 ] spatially co r relat ed Ga ussian noi se. I E E E T ra ns S i g2 Chen J ie . R ese a rc h on s i g n al p rocess i n g an d i m a g i n g ( ) nal P rocess i n g , 1993 , 41 2:549,558 o f ul t ra2w i deba n d ra d a r [ D ] . Beijing : The Graduate 张贤达 . 时间序列分析 ———高阶统计量方法 . 北京 : U niver sit y of Chine se Academy of Sciences , 2007 , [ 19 ] 清华大学出版社 ,1996 26,39 覃天 . 多属性相干方法识别微小断层 . 工程地球物理 [ 8 ] 程建远 . 三维地震资料微机解释性处理技术 .北京 : [ 20 ] () 学报 ,2008 ,5 5:538,542 石油工业出版社 ,2002 Qin Tian. Identification of mino r fault based on multi2at2 [ 9 ] 戴永寿 ,郑德玲 ,魏磊 ,霍志勇 . 高阶统计量地震子波 t ribute seismic co herence technology. Chi nese J ou rnal of ( ) 估计建模 . 石 油 地 球 物 理 勘 探 , 2006 , 41 5 : 514 , () En gi nee ri n g Geop hysics , 2008 , 5 5:538,542 518 ,540 陆基孟 . 地震勘探原理 . 北京 :石油工业出版社 ,2001 Dai Yo ng2sho u , Zheng De2ling , Wei L ei , H uo Zhi2 [ 21 ] yo ng . Estimatio n of high2o r der stati stic sei smic wave2 ( ) let . O G P , 2006 , 41 5:514,518 , 540 ()本文编辑 :金文昱 ( ) 李晓东 高级工程师 ,1963 年生 ; 1991 年毕业于石油大学 华东计算 ,发表论文 10 余篇 。 项目一项 机专业 ;曾就职于地震资料处理员 、地球物理师 、采集项目技术总监 潘树林 讲师 ,1979 年生 ;2008 年毕业于成都理工大学获得地球探 等岗位 ,参加和主持过多项部级 ,局级科研项目 ,在专业期刊及相关 测与信息技术专业博士学位 ; 主要从事多波静校正和去噪方法的 学术会议上发表论文多篇 ,主要从事物探采集技术支持工作 。 研究 ,已公开发表论文 20 余篇 ; 目前在西南石油大学资源与环境 徐维秀 高级工程师 , 1966 年生 ; 1989 年本科毕业于重庆大学应用 学院从事地震勘探相关教学和科研工作 。 数学专业 ,获学士学位 ,1998 年获南京大学计算机应用专业硕士学 于承业 高级工程师 , 1964 年生 ; 1986 年毕业于大庆石油学院勘查 位 ,2007 年获同济大学控制理论与控制工程专业博士学位 ; 长期从 地球物理专业 ,获学士学位 ; 一直在大庆勘探开发研究院从事地震 事数字信号处理及其方法研究 、相关软件开发等工作 ; 现在中国石 资料处理和方法研究工作 ,研究方向为 Q 值估算 , 反 Q 滤波 , A VO 化胜利石油管理局地球物理勘探开发公司从事地震数据采集技术 和油气预测方法 。 ( 研究及相关软件开发工作 。 杜启振 教授 , 博士生导师 ,1969 年生 ; 1991 年毕业于石油大学 华 ( ) ) 张 虎 硕士 ,1983 年生 ; 2004 年毕业于石油大学 华东石油工程 东勘查地球物理专业 , 1997 年获该校应用地球物理专业硕士学 专业 位 ,获双学士学位 ,2007 年毕业于中国科学技术大学固体地球物 ,2003 年获清华大学固体力学专业博士学位 ;现在中国石油大学 ( ) 理专业 ,获硕士学位 ; 现在东方地球物理公司物探技术研究中心从 华东地球资源与信息学院从事地震全弹性波传播与联合成像方 事地震数据采集方法研究和计算机软件开发工作 ; 已发表论文 5 面的教学与研究工作 。 篇 ,申请专利一项 。 陈建阳 工程师 ,1978 年生 ; 2000 年毕业于曲阜师范大学物理工程 ( ) 易昌华 高级工程师 , 1968 年生 ; 1991 年毕业于华东地质学院工程 学院 ,2008 年毕业于中国地质大学 北京能源学院 , 获博士学位 ; 测量专业 ,获学士学位 ; 主要从事石油地震勘探测量导航工作 , 现 2010 年于中国石油勘探开发研究院博士后流动站出站 , 现在中国 任中国石油东方地球物理公司装备事业部测量总工程师 。 石油勘探开发研究院从事油藏描述与储层建模工作 。 赵 虎 博士研究生 代双河 高级工程师 , 1983 年生 ; 2004 年本科毕业于中国矿业大学 , 1966 年生 ; 1988 年毕业于西南石油学院地球 ( ) 地质工程专业 ,2007 年获中国矿业大学地球探测与信息技术专业 物理勘探专业 ,获学士学位 ;2008 年毕业于中国地质大学 北京能 硕士学位 源学院矿产普查与勘探专业 ,2010 年于成都理工大学地球探测与信息技术专业攻读 , 获博士学位 。长期在中国石油集团 博士学位 ; 现在西南石油大学从事地震资料采集处理与解释方面 东方地球物理公司从事石油勘探研究工作 , 现在北京大学地质学 的教学和科研工作 。 博士后流动站的地球与空间科学学院从事科研工作 。 郭向宇 高级工程师 , 1961 年生 ; 1983 年毕业于华东石油学院物探 文晓涛 副教授 ,1976 年生 ;1997 年毕业于成都理工学院物探系 , 获 专业 ,获学士学位 ,2003 年获中国石油大学物探专业博士学位 ; 长 学士学位 ,2003 年获硕士学位 , 并留校任教 , 2006 年获博士学位 , 期从事地震数据处理技术的开发和应用研究 , 近几年参与了三分 研究方向为储层地震预测及流体识别 。 量地震数据处理等油藏地球物理技术的研究和应用工作 。 李国发 副教授 ,博士 ,1966 年生 ; 1987 年毕业于长春地质学院应用 ( ) ( ) 熊章强 教授 ,1963 年生 ;2006 年毕业于中国地质大学 武汉地球探 地球物理系 ,2003,2004 年 ,在中国矿业大学 北京从事博士后阶 ( ) 测与信息技术专业 ,获博士学位 ;现在中南大学地球科学与信息物理 段研究工作 ,2005 年起在中国石油大学 北京资源与信息学院从 学院从事教学和科研工作 ,研究方向为地震波正反演研究及应用 。 事地震勘探教学与科研工作 ; 研究方向为高分辨率地震资料处理 ( 高建军 博士研究生 , 1983 年生 ; 2006 年毕业 于中 国石油 大学 华 和复杂构造地震成像 。 ) ( ( ) 东信息与计算科学专业 ,获学士学位 ,2006 年在中国石油大学 北 孙夕平 博士后 ,1976 年生 ; 1998 毕业于石油大学 华东勘查地球 ) 京攻读地球探测与信息技术专业硕士课程 ,2008 年在资源与信息 物理专业 ,2003 年获该校地球探测与信息技术专业博士学位 ,现从 事复杂地质条件下地震勘探机理 、处理和解释方法研究工作 。 学院获准提前攻读地质资源与地质工程专业博士学位 ; 近期主要 从事多维不规则地震数据规则化重建方面的研究工作 。 许辉群 博士研究生 ,1981 年生 ;2007 年毕业于长江大学 ,获地球探 测与信息技术专业硕士学位 , 同年留校任教 ; 现 为中 国石油 大学 高 秦 工程师 , 1983 年生 ;2005 年毕业于成都理工大学地球物理 ( ) 学专业 ,获理学学士学位 ;现在东方地球物理公司研究院长庆分院 北京物探专业博士研究生 , 主要从事地震信号处理及储层流体 识别方法研究 。 从事地震资料处理分析工作 。 杨 宁 博士研究生 ,1981 年生 ; 2008 年毕业于成都理工大学 , 获硕 黄福喜 博士研究生 , 1975 年生 ; 1998 年毕业于西南石油大学石油 与天然气地质勘查专业 , 现为成都理工大学沉积学专业博士研究 士学位 ;现在成都理工大学油气藏地质与开发工程国家重点实验 生 ,主要从事沉积学 、成藏分析与油气地质条件综合研究工作 , 发 室从事地震正 、反演方法以及各向异性研究 。 冯智慧 博士研究生 , 1981 年生 ; 2006 年毕业于南昌航空工业学院 表论文 10 余篇 。 李建慧 博士研究生 ,1982 年生 ; 2005 年毕业于昆明理工大学资源 信息管理与信息系统专业 ,获学士学位 ;2008 年毕业于吉林大学地 球探测科学与技术学院固体地球物理专业 ,获硕士学位 ; 现为吉林 环境与城乡规划管理专业 ,获学士学位 ,2008 年获中国矿业大学地 球探测与信息技术专业硕士学位 ; 现为中南大学地球探测与信息 大学地球探测科学与技术学院固体地球物理专业博士研究生 , 主 要从事地震信号处理及地球物理勘探方法研究 。 技术专业博士研究生 ,主要从事瞬变电磁法正演研究与资料处理 工作 。 高少武 高级工程师 ,1965 年生 ; 1986 年毕业于西安地质学院 , 获学 士学位 ,1995 年毕业于成都理工学院 , 获硕士学位 , 现在成都理工 索孝东 高级工程师 ,1965 年生 ; 长期从事重磁电及化探资料的综合 () 研究工作 ,现为中国地质大学 武汉石油勘探构造方向博士研究生 。 大学攻读博士学位 ; 在各种国际国内会议和杂志上发表论文 20 余 篇 ,国家专利局授权发明专利 2 个 , 受理发明专利 8 个 ; 现在东方 孙 娅 1984 年生 ; 2007 年毕业于中南大学 , 获地球信息科学与技 地球物理公司物探技术研究中心从事物探处理方法研究与软件开 术工学学士学位 ,同年被保送到中南大学地球探测与信息技术专 发工作 。 业 ,2010 年获工学硕士学位 ,主要从事地球物理电磁法勘探以及数 ( 张 华 博士研究生 , 1979 年生 ; 2004 年毕业 于中 国矿业 大学 徐 值模拟方面的研究 ;现为中南大学地质资源与工程学院博士研究 ) ( 州勘查技术与工程专业 ,获学士学位 ,2007 年获中国矿业大学 徐 生 ,目前是美国 Rice 大学地球科学系为期两年的访问学者 ,主要研 ) ( )州地球探测与信息技术专业硕士学位 ; 现在中国石油大学 北京 究地震波成像方面的技术 。 攻读博士学位 ,研究领域为地震数据处理和成像方法 。 刘振武 教授级高工 ,享受政府津贴专家 ; 毕业于西南石油学院石油 李桂花 1977 年生 ; 2001 年毕业于长安大学应用地 球物 理专业 , 工程专业 ,曾就读英国 Herio t - Wet t 大学 ,获石油工程专业博士学 2008 毕业于西安科技大学地质工程专业资源探测与信息技术方 位 ;主编出版了《二十一世纪初中国油气关键技术》、《油气应用基 向 ,获工学博士学位 ; 现在山东科技大学从事地球物理专业和勘察 础研究展望》《、企业技术创新管理》和《中国石油" 十五" 科技进展 技术与工程专业的教学与科研工作 , 目前主持国家青年科学基金 丛书》等专著 ;现任中国石油天然气集团公司咨询中心副主任 。 U nive r sit y of Tech nolo gy , Che ngdu Cit y , Sic h ua n gle f reque ncy i nt erf ere nce f reque ncy wa s p ropo sed Pro vi nce , 610059 , Chi na i n t hi s p ap e r . The ti me dela y a nd t he a mp lit ude of t he si ngle f reque ncy were calc ulat ed di rect l y u si ng An a utomat ic seismic hor izon picking2up method ti me delay a dap tive fo r mula a nd a mp lit ude fo r mula ba sed on one dimensional sl ice of cross f ourth order re sp ectivel y . 1 1 1 2hui, L iu Ca i, cumulant. Feng Zhi, Feng XuanThe met ho d i s ha r mle ss fo r t he nea r by eff ec2 1 1 Wang D ianand Wang Zhao2guo. O G P , 2011 , tive si gnal s a nd ca n rai se t he si gnal to noi se ratio of ( ) 46 1: 58,63 t he dat a a ro und t he si ngle f reque ncy i nt erf ere nce . In o r der to sol ve t he p ro ble m i n w hich t he Wit h t hi s met ho d t he f reque ncy of t he si ngle i nt er2 sei smic ho rizo n ca n no t co r rect l y picke d up unde r f ere nce ca n be f a st e sti mat ed a nd t he t a r get of sup2 t he st ro ng co here nt Ga u ssia n colo red noi se , a n a u2 p re ssio n of such i nt e rf e re nce wa s reac he d. The t ri2 to matic sei smic ho rizo n pic ki ng2up met ho d ba sed al calc ulatio n re sult s fo r t he synt hetic dat a a nd o n o ne di me n sio nal slice of cro ss fo ur t h o r der cu2 fiel d dat a sho we d t he eff ective ne ss a nd f ea si bilit y mula nt wa s p ropo sed i n t hi s p ap er . The ba sic p ri n2 of t he met ho d . cip le a nd e xact al go rit h m of t he met ho d we re al so Key words : sei smic acqui sitio n , si ngle f reque ncy ela bo rat ed i n t hi s p ap e r . The si mila rit y of t he sei s2 i nt e rf ere nce , si gnal to noi se ratio n , co si ne f unctio n mic ho rizo n o ver t he ho rizo nt al di rectio n a nd t he app ro xi matio n , a dap tive met ho d i n se n sitivit y of t he cro ss fo ur t h o r de r c umula nt to 1 . Che ngdu U nive r sit y of Tec h nolo gy , Che ngdu Ga u ssia n colo red noi se were utilized i n realizi ng t he Cit y , Sich ua n Pro vi nce , 610059 , Chi na ti me delay e sti matio n of t he sei smic eve nt s , mea n2 2 . Geop hy sical Tec h nique Re sea rc h Ce nt re , B GP w hile t he releva nt p a ra met er s w hich di rect l y aff ect Inc . of CN PC , Zh uozho u Cit y , He bei Pro vi nce , t he pic ki ng2up acc uracy of t he sei smic ho rizo n were 072751 , Chi na set . The calculatio n fo r t he t heo retical mo del sho w s t hat t he met ho d ha s ver y st ro ng re st rai ni ng Opt imist ic wavelet ba sis select ion in se ismic signal 1 , 2 1a bilit y fo r t he co here nt Ga u ssia n colo red noi se a nd noise el iminat ion . Zhang Hua, Chen Xiao2hong 2 it s si gnal to noi se ratio reache s 210 dB . The app lica2 ( ) and Yang Ha i2yan. O G P , 2011 , 46 1: 70,75 tio n eff ect s of t he fiel d dat a mat c h well wit h t he The eff ect of utilizi ng wa velet t ec h nolo gy to geolo gical ho rizo n s a nd t he locatio n s of t he f a ult s co nduct sei smic si gnal noi se eli mi natio n i s relat ed a nd t he f a ult poi nt s we re accurat e . It i s no ticed to selectio n of op ti mi stic wa velet ba si s. The eff ect t hat t he picki ng2up accuracy of t he met ho d co ul d of u si ng wa velet t ec h nolo gy to eli mi nat e t he noi se rep lace s t hat of t he ma n ual met ho d , ho weve r it s i n sei smic si gnal i s relat ed to t he selectio n of op ti2 efficie ncy i s muc h hi ghe r t ha n t he ma nual met ho d. mal wavelet ba si s. Ba sed o n t he a nal y si s of cha rac2 Key words : cro ss fo ur t h o r de r c umula nt , o ne di2 t eri stic s of sei smic si gnal , t he p ap er di sc u sse d t he me n sio nal dia go nal slice , a uto matic ho rizo n pick2 selectio n crit erio n of wa velet ba si s t heo reticall y , i ng2up , si gnal to noi se ratio , co here nt Ga u ssia n co nduct e d t he co r relatio n a nd selectio n of wavelet colo re d noi se ba si s f ro m t he t e st view a nd fi nall y p re se nt e d t he 1 . Colle ge of Geo e xp lo ratio n Scie nce a nd Tech nolo2 sugge stio n t hat diff e re nt wavelet ba si s f unctio n gy , J ili n U nive r sit y , Cha ngch un Cit y , J ili n Pro v2 sho ul d be select e d i n deco mpo sitio n a nd reco n2 i nce , 130026 ,Chi na st r uctio n of wa velet t ra n sfo r m to a ssure t he accu2 racy of t he reco n st r uct e d si gnal a nd e nha nce t he Ident if icat ion and suppression of single f requency in2 p roce ssi ng cap a bilit y fo r sei smic si gnal . The new terf erence by using cosine f unct ion a da pt ive method. 1 , 2 2 2 2app roach c ha nge d t he p a st ro uti ne wa y t hat t he Gao Shao2wu, Zhao Bo, Zhu Shu2yun,Luo Guo2an 1 sa me wavelet ba si s f u nctio n wa s u sed i n deco mpo2 ( ) an d He Zhen2hua. O G P , 2011 , 46 1: 64,69 sitio n a nd reco n st r uctio n of wavelet t ra n sfo r m. The st ro ng si ngle f reque ncy i nt erf ere nce a2 The bio r2 . 4 wavelet ba si s w hich i s suit a ble fo r ro und 50 Hz ha s bee n alway s t he st udy o bject i n sei smic si gnal noi se eli mi natio n i s fi nall y select e d i n sei smic dat a p roce ssi ng. The cur re nt sof t wa re mo dule s w hic h a re u se d to supp re ss 50 Hz st ro ng t hi s p ap er . The co r relatio n re sult s of fiel d sei smic si ngle f reque ncy i nt erf ere nce have so me def ect s , dat a p roce ssi ng sho w t hat app l yi ng t he bio r2 . 4 such a s lo w efficie ncy , la r ge ti me2co n sumi ng dur2 wavelet ba si s i n sei smic dat a p roce ssi ng no t o nl y a2 i ng t he calculatio n a nd bei ng ha r mf ul to eff ective c hieve s t he bet t er eff ect of noi se eli mi natio n , but si gnal s a nd so o n . The co si ne f u nctio n a dap tive al2 ca n e xt ract t he eff ective co mpo ne nt of lo w a nd hi gh go rit h m w hich ca n be u se d to f a st e sti mat e t he si n2 f reque ncy sei smic si gnal f ro m re sult dat a , broa de n