MATLAB软件在鹤岗地震台地震观测记录中的应用

随着数字化地震计的应用和发展,数字观测资料中不仅记录了宽频带地震信息,也包含了大量的干扰信号。在实际工作中,对地震信息的分析处理影响很大,有时还会造成地震波形记录不完整或地震信号淹没在干扰信号中,很难再进行深入的分析和应用。因此,利用MATLAB软件设计IIR数字滤波器对不同记录中的干扰波进行滤波器设计,达到消除干扰波的目的。实例表明,此方法在台站应用效果良好。     

Filter formulation and wavefield separation of cross-well seismic data

Multichannel filtering to obtain wavefield separation has been used in seismic processing for decades and has become an essential component in VSP and cross-well reflection imaging. The need for good multichannel wavefield separation filters is acute in borehole seismic imaging techniques such as VSP and cross-well reflection imaging, where strong interfering arrivals such as tube waves, shear conversions, multiples, direct arrivals and guided waves can overlap temporally with desired arrivals. We investigate the effects of preprocessing (alignment and equalization) on the quality of cross-well reflection imaging wavefield separation and we show that the choice of the multichannel filter and filter parameters is critical to the wavefield separation of cross-well data (median filters, f–k pie-slice filters, eigenvector filters). We show that spatial aliasing creates situations where the application of purely spatial filters (median filters) will create notches in the frequency spectrum of the desired reflection arrival. Eigenvector filters allow us to work past the limits of aliasing, but these kinds of filter are strongly dependent on the ratio of undesired to desired signal amplitude. On the basis of these observations, we developed a new type of multichannel filter that combined the best characteristics of spatial filters and eigenvector filters. We call this filter a ‘constrained eigenvector filter’. We use two real data sets of cross-well seismic experiments with small and large well spacing to evaluate the effects of these factors on the quality of cross-well wavefield separation. We apply median filters, f–k pie-slice filters and constrained eigenvector filters in multiple domains available for these data sets (common-source, common-receiver, common-offset and common-midpoint gathers). We show that the results of applying the constrained eigenvector filter to the entire cross-well data set are superior to both the spatial and standard eigenvector filter results.     

均衡多道1范数匹配多次波衰减的方法与应用研究

反馈迭代法多次波衰减分为预测和相减两个步骤.在匹配相减过程中,当多次波和一次波同相轴不满足正交性,以及去除多次波后的地震记录不满足能量最小的情况下,最小二乘自适应匹配滤波方法不能获得正确的匹配,因此多次波的衰减产生误差.针对这两个问题,本文提出了基于1范数最小的均衡多道自适应匹配滤波法.该方法通过在空间方向上对地震记录进行均衡,避免了多次波和一次波同相轴正交的假设条件.同时利用了1范数对于大异常值保持稳健的特点,因此可以有效地解决去除多次波后的地震记录能量最小的问题.通过对模型数据和野外实际数据的多次波压制结果显示,该方法可以有效、准确地衰减多次波.     

基于多道卷积信号盲分离的多次波自适应相减方法

本文将多次波自适应相减问题表示为一个多道卷积信号的盲分离问题.利用2D卷积核来表示预测多次波和实际多次波之间的差异,并采用分离出的一次波信号的非高斯性最大化作为优化目标,我们提出一种基于多道卷积信号盲分离的多次波自适应相减算法.为了求解上述非线性优化问题,所提方法将其转化为一个迭代线性优化问题,采用迭代最小二乘方法加以实现.由于采用了多道卷积信号盲分离模型,所提方法能够适应预测和真实多次波之间在时间及空间上的变化.通过对简单模型数据、Pluto数据和实际数据进行处理,验证了所提算法的有效性.     

用小波变换能量线性度方法识别天然地震与爆破或塌方

为了提取天然地震和爆破或塌方记录波形在震源深度、震源尺度、震源破裂机制、地震波传播途径、地震波衰减等方面的差异特征信息,本对山东数字化台网记录的天然地震和爆破或塌方波形进行了小波多分辨率的能量线性度特征分析,提出了用小波变换能量线性度方法识别天然地震与爆破或塌方事件．结果表明：在精细结构小波分解信号“能量”线性度方面,天然地震主要集中在-2．0～1．0之间,爆破或塌方主要集中在2．0～3．4之间;在精细结构小波分解信号“能量”最大值对应的小波分解尺度方面,爆破或塌方主要集中在4～5,频段集中在0．7～3．1Hz之间,而天然地震主要集中在1～2,频段集中在6．25～25Hz之间．     

Perfect partial reconstructions for multiple simultaneous sources

A major focus of research in the seismic industry of the past two decades has been the acquisition and subsequent separation of seismic data using multiple sources fired simultaneously. The recently introduced method of signal apparition provides a new take on the problem by replacing the random time-shifts usually employed to encode the different sources by fully deterministic periodic time-shifts. In this paper, we give a mathematical proof showing that the signal apparition method results in optimally large regions in the frequency–wavenumber space where exact separation of sources is achieved. These regions are diamond shaped and we prove that using any other method of source encoding results in strictly smaller regions of exact separation. The results are valid for arbitrary number of sources. Numerical examples for different number of sources (three, respectively, four sources) demonstrate the exact recovery of these diamond-shaped regions. The implementation of the theoretical proofs in the field is illustrated by the results of a conducted field test.     

基于余弦函数的自适应单频干扰消除

在野外地震数据采集过程中,如果地震测线上空有高压输电线通过时,在地震记录中就存在一个50HZ左右强单频干扰.这种干扰在地震记录的整个时间段上具有很强的能量,并且从浅层到深层频率、相位和振幅可能不相同,但是它们的变化不太大,可以认为基本保持不变.采用频率域压制方法很难完全消除这种干扰,仅仅是噪声和信号同时压制一定的倍数,因此并没有改善单频分量频率附近的信噪比.我们提出了基于余弦函数的自适应单频干扰消除方法,以有效消除地震数据单频干扰.在时间域内采用频率、时延和振幅的余弦函数来逼近这种单频干扰并从地震记录中减去,以实现消除这种单频干扰的目的.单频干扰余弦函数的振幅使用公式直接计算,单频干扰频率和时延采用自适应方法估算.本方法的最大优点就是能够有效消除干扰,而并不伤害干扰附近有效信号,因此在单频干扰频率分量附近数据信噪比提高了.通过实际数据的试算,说明了本方法的有效性和可行性.     

基于地震信号波形形态差异的面波噪声稀疏优化分离方法

实际地震信号通常可表示为具有波形特征差异的多种基本波形信号的线性组合,如叠前道集中的工频干扰噪声与有效波信号、面波噪声与体波信号等.选择单一数学变换方法,往往不易实现地震信号的稀疏表示.近年来发展的形态成分分析理论,通过联合多种数学变换,可实现对复杂信号的稀疏表示.本文根据单道地震记录中面波与体波信号波形结构特征的差异性,提出一种基于形态成分分析的面波噪声衰减方法.针对面波的低频、窄带以及频散特性选择一维平稳小波变换作为其稀疏表示字典,而针对体波波形的局部相关特性选择局部离散余弦变换作为其稀疏表示字典,建立基于双波形字典的形态成分分析模型,通过求解该稀疏优化问题获得最终的信噪分离结果.理论模型和实际地震资料处理证实该方法不仅能够衰减单炮地震记录中的强面波干扰噪声,同时能够更好地保护有效信号的波形特征与频谱带宽,为地震资料的后续处理和分析提供良好的数据基础.     

Improved risk‐targeted performance‐based seismic design of reinforced concrete frame structures

This paper presents a procedure for seismic design of reinforced concrete structures, in which performance objectives are formulated in terms of maximum accepted mean annual frequency (MAF) of exceedance, for multiple limit states. The procedure is explicitly probabilistic and uses Cornell's like closed‐form equations for the MAFs. A gradient‐based constrained optimization technique is used for obtaining values of structural design variables (members' section size and reinforcement) satisfying multiple objectives in terms of risk levels. The method is practically feasible even for real‐sized structures thanks to the adoption of adaptive equivalent linear models where element‐by‐element stiffness reduction is performed (2 linear analyses per intensity level). General geometric and capacity design constraints are duly accounted for. The procedure is applied to a 15‐storey plane frame building, and validation is conducted against results in terms of drift profiles and MAF of exceedance, obtained by multiple‐stripe analysis with records selected to match conditional spectra. Results show that the method is suitable for performance‐based seismic design of RC structures with explicit targets in terms of desired risk levels.     

基于能量比的天然地震与人工爆破自动识别算法研究

研究了从天然地震和人工爆破事件的波形记录中提取出来的能量比特征在天然地震和人工爆破事件的自动识别中的有效性及适用性。对波形记录进行了4层小波变换,然后对变换得到的小波系数提取能量比特征,最后利用支持向量分类机ν-SVC进行识别效果检验。实验证明,由bior2.2小波包分解后提取出来的能量比特征对天然地震和人工爆破事件的识别效果很好,可用于实际的自动识别系统作为识别判据之一。     

首都圈数字地震台网对微弱爆破信号的检测能力

利用首都圈数字地震台网接收人工地震信号,进行地下结构研究具有重要意义.但人工震源释放的能量小,激发的地震波以短周期为主,因此本文较全面地研究了地震台网对短周期微弱信号（1～20 Hz）的检测能力：(1) 分析了台网的背景噪声,结果表明基岩台址的地震台噪声比沉积盖层台址的地震台噪声低约13 dB,这相当于近1个震级的检测阈值;夜间的噪声比白天低约5 dB;噪声有逐年增高的趋势,2006年比2001年噪声提高约4 dB.(2 )分析了在台网内进行的药量为25 kg的陆地井下爆破实验,一次爆破相当于0.69级（M_L）的天然地震,有18个地震台可辨认爆炸产生的Pg、Pm或Pc波;离爆破点218 km的基岩台,仍可以接收到振幅只有1.6 nm 的Pm波,这个结果可为地震勘探实际工作提供参考.(3) 研究了台网外核爆试验的信号特征,2006年发生在朝鲜的地下核试验是一次检验台网检测微弱信号能力的好机会.波形记录经1～5Hz滤波后,台网中噪声小的18个基岩台可以清晰辨认核爆破产生的P波或Lg波,P波平均振幅为16 nm,计算的平均震级为m_b4.3,和NEIC给出的震级相同;分析还表明背景噪声是影响台站信号检测能力的主要因素之一.     

带有多道相关的抛物线Radon变换法分离P-P、P-SV波

因为多波多分量地震勘探中P-P波和P-SV波通常混杂在一起,所以较好地分离P-P波和P-SV波能够提高数据处理和解释的质量.抛物线Radon变换法在分离P-P波和P-SV波时取得了一定效果,但是在离散叠加的计算过程中会带来假频,这些假频会干扰波场分离.本文针对这一问题,将多道相关算法引入抛物线Radon变换,发展了带有多道相关的抛物线Radon变换法.该方法利用叠加信号具有相似性的特点,依据多道相关中衡量多道信号相似性的能量比标准,对叠加过程加以控制,压制变换中出现的假频.本文用该方法对合成地震记录进行了波场分离,取得了较好的效果.     

利用小波变换分析朝鲜核爆波形特征

基于小波变换方法,对2009-2017年黑龙江省牡丹江地震台（MDJ）记录的5次朝鲜核爆信号进行Daubechies小波8层多尺度分解,通过对比分析近似和细节部分信息,发现核爆信号进行小波分解后,在不同层存在相同特征,且有别于地震信号分解后在相应各层的特征。可见,利用该方法可有效识别朝鲜核爆。     

基于压缩感知和宽带俞式低通整形滤波器的地震低频信息特征分析与补偿

碳酸盐岩储集层已成为世界石油新发现储量的重要组成部分,识别该类储层对地震数据的信噪比、分辨率以及成像精度提出了更高的要求.本文从地震低频信号缺失的问题出发,首先研究了低频信号缺失对子波、合成地震记录和波阻抗反演的影响,其次分析了深层碳酸盐岩裂缝储层中弱信号低频缺失的特征.针对低频信号缺失问题,本文利用压缩感知理论,并结合反射系数的稀疏特性,提出了自适应计算L1范数权重因子的方法,同时构建了改进的宽带俞式低通整形滤波器,在不影响地震高频信号的同时对地震弱信号进行低频补偿.结果表明,缺失低频信号,会使子波旁瓣变大,合成记录出现假同相轴,厚层波阻抗反演畸变,深层碳酸盐岩裂缝储层弱信号难以识别;而本文方法有效地补偿了深层碳酸盐岩裂缝储层弱信号10Hz以下的频率成分,使得波组反射特征更加清晰,深层弱信号成像质量得到改善,为进一步有效识别深层碳酸盐岩裂缝储层建立了基础.     

压制地震勘探随机噪声的分段时频峰值滤波方法

提高勘探资料信噪比是地震勘探的主要内容之一.本项研究利用短时能量将记录分成能量均衡的若干段,获得减少尺度变换误差的分段时频峰值滤波方法,并运用端点拓展和与地震信号特征匹配的多级时窗参数改进时频峰值滤波精度.理论模型和共炮点资料处理结果表明,分段时频峰值滤波很好地消除了尺度变化误差引起的信号波形阶梯状畸变,能够在压制强随...     

编码震源提高地震探测能力的野外实验研究

地球深部结构研究是地震预报和地球动力学研究的基础，传统人工地震受到探测深度、分辨率或震源能量的环保限制，探索利用小能量震源进行深部高分辨率探测技术具有重要意义.本文借鉴雷达探测中的编码原理，开展了震源编码探测实验，利用人工震源向地下发送具有独特特征的随机编码序列脉冲，作为一次等效激发，利用准确记录的编码震源的源函数与地震连续观测的互相关进行“解码”分析，获取地下信息，形成编码地震探测技术.实验表明编码地震能够大幅度提高地震探测能力，利用小能量震源获得远距离的探测和高分辨率，是一种具有重要发展前景的人工探测技术.     

Orthogonal vibroseis sweeps

Vibroseis is a method that imparts coded seismic energy into the ground. The energy is recorded with geophones and then processed using the known (coded) input signal. The resulting time‐domain representation of vibroseis data is an impulsive wavetrain with wavelet properties consistent with the coded input signal convolved with the earth's reflectivity series. Historically, vibratory seismic surveys collect data from one source location at a time, summing one or more sources at each location. We present a method of designing orthogonal sweeps using the concept of combisweeps. The orthogonal sweeps allow simultaneous recording and later separation of two or more unique source locations. Orthogonality of sweeps permits separation of the data into unique source‐location field records by a conventional correlation procedure. The separation power of the orthogonal sweeps is demonstrated by a comparison between separated data and data acquired with one vibrator. Separation noise was at a negligible level for our demonstration data sets when two vibrators were located 50 m to 200 m apart. Coincident generation and recording of two vibroseis sweeps at different locations would allow almost double the amount of data to be recorded for a given occupation time and requires only half the storage medium.     

Salt body detection from seismic data via sparse representation

In seismic interpretation and seismic data analysis, it is of critical importance to effectively identify certain geologic formations from very large seismic data sets. In particular, the problem of salt characterization from seismic data can lead to important savings in time during the interpretation process if solved efficiently and in an automatic manner. In this work, we present a novel numerical approach that is able to automatically segmenting or identifying salt structures from a post‐stack seismic data set with a minimum intervention from the interpreter. The proposed methodology is based on the recent theory of sparse representation and consists in three major steps: first, a supervised learning assisted by the user which is performed only once, second a segmentation process via unconstrained ?₁ optimization, and finally a post‐processing step based on signal separation. Furthermore, since the second step only depends upon local information at each time, the whole process greatly benefits from parallel computing platforms. We conduct numerical experiments in a synthetic 3D seismic data set demonstrating the viability of our method. More specifically, we found that the proposed approach matches up to 98.53% with respect to the corresponding 3D velocity model available in advance. Finally, in appendixes A and B, we present a convergence analysis providing theoretical guarantees for the proposed method.     

空间局部加权回归自适应TFPF

时频峰值滤波(TFPF)算法是一种非常有效的去噪方法.但是传统的TFPF采用的单一窗长,并且仅沿时间方向进行滤波,忽略了信号的空间信息,并且TFPF近似等效成一个时不变的低通滤波器,不能追踪快速变化的信号.针对这些问题,引入空间局部加权回归自适应TFPF(SLWR-ATFPF).鉴于随机噪声在各个位置的方向随机性,以及有效信号在各个位置的方向确定性,首先利用空间局部加权回归(SLWR),对含噪信号进行空间加权,从而使加权之后的信号包含空间信息.然后,再引入凸集和Viterbi的思想,对空间加权之后的信号进行自适应滤波.从而,完成时空域二维自适应滤波.将SLWR-ATFPF应用于合成记录和实际的共炮点记录,实验结果表明,改进的方法与原算法相比,能够在压制低信噪比(SNR)随机噪声的同时更好地保留有效信号.     

基于VMD的随掘地震超前探测信号谱白化方法研究

随掘地震超前探测震源具有连续、非可控的特点,需要使用互相关技术将连续震源记录转换为等效脉冲震源记录以获取有效反射信号。但随掘地震超前探测信号往往包含一些能量较强且频带较窄的优势频率成分,使得互相关处理会引入较严重的干扰。为此,本文提出将基于变分模态分解（VMD）的希尔伯特谱白化方法应用于随掘地震超前探测信号。该方法首先利用VMD将地震信号分解为若干个本征模态函数（IMF）,再对各IMF应用希尔伯特变换进行时频分解,最后使用白化滤波器对其希尔伯特谱进行谱白化。数值模拟结果表明基于VMD的希尔伯特谱白化方法能够在保持各道信号一致性的同时均衡信号的不同频率成分,有效压制互相关结果中的虚假同相轴,加快主峰旁瓣衰减。将本文方法应用于安徽淮北杨庄煤矿某巷道实际数据,成功探测了掌子面前方存在的断层构造,表明该方法具有较好的实用性。