基于Matlab的伪随机扫描信号的仿真与研究.pdfVIP

石油天然气学报 (江汉石油学院学报) 2010年12月 第32卷 第6期 JournalofOilandGasTechnology (J。JPI) Dec．2010 Vo1．32 No．6 基于 Matlab的伪随机扫描信号的仿真与研究 蹇清平 ，余新平，刘益成 (长江大学电子信息学院，湖北荆州434023) [摘要]针对可控震源扫描信号的要求和伪随机序列的特点，提 出了一种基于Matlab的伪随机扫描信号的 设计方法，该方法是将生成的长时间 M 序 列通过 频谱整形 的方 法，形 成频带 限制、 自相关性 良好 的伪随 机扫描信号。在 Matlab平 台上完成 了扫描信号的仿真实验 。仿真结果表 明，伪随机扫描信号 的 白相关性 良好 ，信号的相关噪声、扫描 时间和频带宽度均满足可控震源 的要求 。 [关键词]可控震源；频谱整形 ；伪随机扫描 ；Matlab [中图分类号]P631．44 [文献标识码]A [文章编号]1000—9752 (2010)O6—0387—04 可控震源扫描信号是振幅均衡，长时间有限带宽的信号 ]，它主要依靠长时问的振动激发来得到相 对弱的带宽有限的地震信号 ，信号的持续时间很长 (8～64s)，远大于地下各层面问的反射时差，并 且由于高频信号在地层中的传播衰减很大 ，因而要求信号的频带宽度 限制在 3～250Hz。可控震源的核 心技术是脉冲 缩 ，对于可控震源的原始记录 ，必须对扫描信号进行压缩才能类 比于脉冲源记录，并用 于后续的处理。自可控震源问世 以来，相关算法一直是信号压缩的主要方式 2l】，是 由可控震源原始记录 与扫描信号互相关形成相关子波来实现的，相关子波的分辨率决定了可控震源地震勘探 的分辨率 ，相关 子波应该具有非常尖锐 的主瓣和尽可能小 的旁瓣和相关噪声 ，才能获得地震勘探的高分辨率。可控震源 有线性扫描、非线性扫描、伪随机扫描 3种扫描方式 ，其中伪随机扫描是非线性扫描方式中的一种 。伪 随机信号具有 良好的 自相关性 ，它又是周期的，其长度可控制 ，但是其频谱较宽 ，针对其特点和缺陷， 笔者研究了一种基于频谱整形 的伪随机扫描信号的设计方法 。 1 伪随机序列及其 自相关函数 伪随机序列 由0和 1两个元素组成二元 (或二进制)序列 ，它可以用波形进行模拟。通常把二元序 列 (或波形)中的元素称为比特或码元，而把它的持续 时间称为 比特长度 。伪 随机序列具有 良好 的 “随机性”，它的自相关 函数接近 白噪声的 自相关函数 ，伪随机序列的伪随机性表现在它有一定 的规律， 即预先的可确定性和可重复性，使它易于从其他信号和干扰 中分离出来，实现相关接收或匹配接收，故 有好的抗干扰性能。 以相关函数作为划分标准 ，将伪随机序列分为 3类 ：①最佳 PACF (周期 自相关 函数)，如 M 序 列 、差分级序列等 ；② 良好 的PACF和 PCCF (周期互相关 函数)，如 Gold序列 、Kasami序列；③大 线性复杂度 ，如 GWM 序列 。 为了产生满足可控震源要求的伪随机扫描信号 ，选择了具有最佳 自相关函数 、序列长度可控制 的， 且易于实现的 M 序列来研究 。 M 序列是最长线性反馈移位寄存器序列的简称 ，它是 由带线性反馈的移存器产生的周期最长 的一 种序列 。M序列生成器的结构如图 1所示。它 由 个移位寄存器级联构成 ，当离散时间增加一个时间 单位 ，移存器的内容向右移一位 ，输出一个码元 ，同时反馈路径产生的反馈信号移进最左边的寄存器 。 在移位寄存器的初始状态不全为零，且反馈路径由本原多项式决定的前提下 ，移位寄存器回到初始状态 [收稿日期]201o～07—02 [基金项 目]中国石油科技创新基金项目 (2OO8D一5006—03—07)。 [作者简介]蹇清平 (1986一)，女 ，2008年大学毕业，硕士生，现主要从事信号与信息处理方面的研究工作。 ● ·388 · 石油天然气学报 (江汉石油学院学报) 2010年 12月 前会 连 续 出现 (2 一 1)个 不 同的非 零 状 态 ，所以M 序列的周期也是 (2一1)；但 在一个周期内，它又是纯随机的。 M序列的自相关函数为 ： fI1 J一 0