脉冲编码震源的匹配冲击技术

常规Vibroseis可控震源Chirp扫描的地震响应存在明显的旁瓣效应,二元m-序列伪随机扫描的解码地震剖面中存在严重的相关噪声,而基于伪随机编码的匹配扫描方法可有效压制解码地震剖面中的相关噪声,提高了地震记录的信噪比.为了将匹配伪随机编码方法应用到冲击式震源中,提出了基于脉冲编码震源的匹配冲击技术.以匹配扫描方法为理论基础给出了匹配冲击信号的编码与解码过程,以有限元差分法对匹配冲击技术和常规线性扫描冲击技术进行了地震数值模拟.结果表明,匹配冲击的综合解码地震剖面中,不存在线性扫描冲击解码地震剖面中能量较强的相关噪声干扰,剖面具有很高的信噪比和分辨率.地震波激发采集对比实验验证了匹配冲击技术在压制相关噪声及远偏移距地震道随机噪声方面相比线性扫描冲击技术所具有的优势.     

伪随机编码震源信号的地震响应

常规的可控震源Chirp扫描信号的地震响应存在着相关信号旁瓣大、分辨率低的缺点.本文利用数值模拟技术研究可控震源伪随机编码信号地震响应的规律和特点.伪随机编码信号的地震响应剖面的旁瓣电平与常规的Chirp扫描信号情形相比有明显地降低,但是伪随机编码信号地震响应剖面的相关噪声仍然没有得到有效地压制.利用伪随机编码信号的周期性特点把一个长地震记录以信号周期为间隔进行分时叠加可以提高地震记录的信噪比和地震数据的利用率.     

可控震源匹配扫描方法研究

常规的可控震源Chirp扫描信号的地震响应存在着相关信号旁瓣大、分辨率低的缺点;用二元伪随机编码信号作为可控震源的驱动信号,地震响应剖面中存在严重的相关噪声,降低了地震记录的信噪比.本文提出的可控震源匹配扫描方法,利用匹配伪随机序列偶调制正弦载波信号作为可控震源驱动信号进行扫描激发,得到两个匹配的原始地震记录,对两个匹配的原始地震记录分别解码后再进行叠加形成最终的综合解码地震剖面.地震数值模拟的结果表明,该方法基本不存在常规Chirp扫描情形的旁瓣效应,可有效避免m-序列编码扫描情形的相关噪声,其最后形成的综合解码地震剖面可以与脉冲震源的炮集记录相媲美.     

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

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

基于三元伪随机编码的冲击夯可控震源控制系统设计

冲击夯可控震源具有移动方便、质量轻、输出力大等优点,在城市工程勘察以及中浅层地震勘探等领域具有一定的应用前景.其激励信号主要基于二元伪随机编码技术,由于旁瓣效应导致相关干扰严重、地震记录信噪比低等问题,制约了冲击夯震源的推广应用.本文针对传统控制方法的不足,开展了基于三元伪随机的冲击夯可控震源编码技术和编码匹配控制方法研究,并在此基础上设计和实现了相应的控制系统.野外实验结果表明:相比于二元伪随机编码控制方法,三元伪随机编码控制方法得到的地震子波的主旁瓣比提高了4.04 dB;震源信号穿透力强,在测线远端处能量约为30磅锤击震源的100倍,约为二元伪随机方法的4倍;单炮地震记录中干扰抑制能力强,信噪比显著提升,子波频带宽,提高了分辨率,在一定程度上提高了地震勘探效果.     

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

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

地震勘探中相控阵震源的方向特性研究

电磁驱动式可控震源在城市浅层地震勘探中所面临的最突出的困难是微弱的反射信号常常淹没在很强的背景噪声之中.为了提高地震记录的信噪比,可以利用多台可控震源阵列实施相位控制形成定向地震波束以增强地震波的能量.本文讨论这种相控阵震源的波束形成机制.引入了地震波场的边际能量密度的概念,利用地震波场的时间切片技术,对模型空间各个方向上的能量强度进行了定量分析.用有限差分法对相控阵震源Chirp信号扫描的地震响应进行了数值模拟.当定向地震波束的汇聚带与观测排列的空间范围相一致时,相控阵震源合成地震记录的能量强度要显著高于单个可控震源情形的能量强度,波形振幅的均匀性要明显优于常规组合激发震源情形波形振幅的均匀性.     

一种基于同相轴追踪的炮混叠记录分离方法的研究与应用

海上地震勘探中为提高采集效率会采取连续采集方案,导致地震记录尾端与下一炮起始记录重合,这种现象对浅层地震处理影响不大,但是会完全掩盖深层地震信号.由于混叠噪声与真实信号属于同方向及同震源采集,其频谱完全重叠,常规的噪声压制方法无法完全压制混叠噪声.本文提出了一种基于同相轴追踪的混叠记录分离方法,利用时域上混叠噪声与真实地震信号时距曲线的曲率差异得到混叠噪声的模型,再利用最小平方约束反演方法进行自适应相减,使有效信号误差最小,最终完成混叠炮集分离.通过理论分析与正演测试,混叠衰减的关键在于求取可靠的混叠模型,而自适应衰减则控制细微振幅与相位差异.通过对实际地震资料进行处理,并与F-K滤波方法进行了对比,证实了该方法的有效性.     

井间地震方向性震源模拟方法

为了能够在水平方向上传播更远的距离,井间地震震源辐射能量具有水平方向强的特性,这是井间地震资料大入射角反射波能量强的主要原因.为了研究井间地震震源方向特性对地震信号的影响,本文提出了圆形震源阵列模拟方法.基于声波方程的一阶应力一速度方程的交错网格有限差分算法,进行了均匀各向同性无限介质的直达波和反射波正演数值模拟.数值模拟结果表明,井间方向性震源辐射能量的方向特性影响了地震记录的能量分布,为进一步研究方向性震源的能量补偿提供了理论依据.     

基于Seislet变换的稀疏约束多震源最小二乘逆时偏移

多震源地震采集技术允许一次性激发不同位置处的震源,得到来自多个震源的混合地震数据,该技术采集效率高,能有效降低采集成本.多震源地震数据成像效率高,但在偏移剖面中会引入串扰噪声,影响成像精度.最小二乘偏移常被用于压制多震源地震数据成像中的串扰噪声,但常规的最小二乘偏移并不能很好的消除串扰噪声对成像结果的影响,难以满足成像精度的要求.因此,为了保证反演的稳定性并改善反演结果,根据反射系数在Seislet域的稀疏性,本文引入了Seislet变换作为变换域稀疏约束的变换算子,实现了基于Seislet变换的稀疏约束多震源最小二乘逆时偏移,数值实验表明该方法能有效压制串扰噪声.     

MAXIMUM LIKELIHOOD ESTIMATION OF SEISMIC REFLECTION COEFFICIENTS*

A seismic trace is modeled as a moving average (MA) process both in signal and noise: a signal wavelet convolved with a reflection coefficient series plus colored random noise. Seismic reflection coefficients can be estimated from seismic traces using suitable estimation algorithms if the input wavelet is known and vice versa. The maximum likelihood (ML) algorithm is used to estimate the system order and the reflection coefficients. The system order is related to the arrival time of the latest signal in a complex seismic reflection event. The least-squares (LS) method does not provide such information. The ML algorithm makes assumptions only about the Gaussian nature of the noise. It is better suited for seismic applications since the LS method inherits the white noise assumption. The Gauss-Newton (G-N) and Newton-Raphson (N-R) optimization algorithms are used to obtain the ML and the LS estimates. Reflection coefficient estimations are affected by the choice of sampling rate of seismic data. Theoretically, the optimum choice in system identification is the Nyquist rate. Experience with synthetic data confirms the theory. In practice, good estimates of reflection coefficients are possible only up to certain pulse separations (or, equivalently, orders). This is mostly due to numerical problems with the optimization algorithms used and partly due to the limited bandwidth of seismic signals. Good estimates from data simulated using three airgun array pulses recorded with 6–128 Hz filter setting are possible up to about 40.0 ms pulse separations. Successful estimations from pinchout and thin layer simulations and well controlled offshore “bright-spots” are given.     

IDENTIFICATION OF SEISMIC REFLECTIONS USING SINGULAR VALUE DECOMPOSITION*

Singular value decomposition (SVD) is applied to the identification of seismic reflections by using two different models: the impulse response model, where a seismic trace is assumed to consist of a known signal pulse convolved with a reflection coefficient series plus noise, and the delayed pulse model, where the seismic signal is assumed to consist of a small number of delayed pulses of known shape and with unknown amplitudes and arrival times. SVD clearly shows how least-squares estimation of the reflection coefficients may become unstable, since a division by the singular values is required. Two methods for stabilizing this procedure are investigated. The inverse of the singular values may be replaced by zeros when they are less than a given threshold. This is called the SVD cut-off method. Alternatively, we may use ridge regression which in filter design corresponds to assuming white noise. Statistical methods are used to compute an optimal SVD cut-off level and also to compute an optimal weighting parameter in ridge regression. Numerical studies indicate that the use of SVD cut-off or ridge regression stabilizes the least-squares procedure, but that the results are inferior to maximum-likelihood estimation where the noise is assumed to be filtered white noise. For the delayed pulse model, we use a linearization procedure to iteratively update the estimates of both the reflection amplitudes and the arrival times. In each step, the optimal SVD cut-off method is used. Confidence regions for the estimated reflection amplitudes and arrival times are also computed. Synthetic data examples demonstrate the effectiveness of this method. In a real data example, the maximum-likelihood method assuming an impulse response model is first used to obtain initial estimates of the number of reflections and their amplitudes and traveltimes. Then the iterative procedure is used to obtain improved estimates of the reflection amplitudes and traveltimes.     

PULSE CODING IN SEISMIC PROSPECTING SOSIE AND SEISCODE *

The conventional seismic technique is subject to a recording time following each transmission of energy, in which it is forbidden to release any new pulse. The recording time depends on the deepest reflection of interest, and is often 10 seconds or more in actual practice. To each transmission corresponds one record, i.e. a fixed amount of data which cannot be increased in a given time. Pulse coding allows us to go beyond this limit, by transmitting several times during the normal recording time. The procedure gives as many records as there are pulses, but they overlap, each event being repeated every time there is a pulse. It is possible to process the composite record back to its usual appearance with all events in their proper place if the time breaks are accurately known and make up a code such that the unavoidable noise generated by the process be kept, on the final section, below the ambient noise. The processing is quite similar to that of records made from vibrating sources, though faster in practice. The additional information can be devoted to a saving of time and money as the same profile may be recorded in a shorter time; or to an improvement of quality of the section due to a higher order of coverage, a multiplication of the ray paths and a closer sampling of the reflectors. It is also possible to record information in several planes at the same time, and to work out a 3-dimensional restitution, without loss of production. The process applies to all kinds of sources provided they can be triggered according to the code with sufficient accuracy. Depending on the source and conditions of implementation, the method benefits from other advantages such as better resolution, increased flexibility, and better coupling. Two different names have been given to the process, Sosie and Seiscode , which apply to slightly different parameters for the sequence of pulses. Sosie is more useful at sea, while the normal scope for Seiscode is onshore. Both names are trademarks for SNPA.     

Coherence between the fields of low-frequency seismic noise in Japan and California

The relationship is considered between the statistics of the field of low-frequency seismic noise which was synchronously recorded by two broadband seismic networks in Japan (78 stations) and California (81 stations). The analysis is based on the data for seven years of observations (2008–2014). For each network, the daily time series of the median values are constructed for five parameters of seismic noise: kurtosis (excess), minimal normalized entropy of the distribution of the squared wavelet coefficients, generalized Hurst exponent, support width of the singularity spectrum, and index of linear predictability. The median values for each parameter were calculated on a daily basis over all the stations of the networks and resulted in a time series containing 2557 data points of the integral characteristics of the noise with a daily time step. The use of the median values of the noise parameters avoids considering the effects of the gaps in recording by individual stations and provides the continuous time series as the integral characteristic of the whole network. Next, for each network, the aggregate signals were calculated for the obtained five-variate time series. By construction, the aggregate signal is a scalar signal which maximally accumulates the most general variations that are simultaneously present in all the analyzed signals and, at the same time, rejects the components that are only characteristic of a single process. The final step of the analysis consists in estimating the evolution of the quadratic spectrum in the moving time window with a length of one year. It is shown that during the considered interval of the observations, the coherence is characterized by the increasing linear trend, which independently supports the previous conclusion about the enhancement of the synchronization between the parameters of the global seismic noise.     

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.     

一种神经网络改进小波的地震数据随机噪声去除方法

地震资料的有效信号反射弱,且易受多次波的影响,不可避免地存在随机噪声干扰。提出一种基于神经网络改进小波的地震数据随机噪声去除方法,采用神经网络模型,识别出随机噪声信号,对该信号进行小波包分解,获取多类别随机噪声信号,采用级联BP神经网络模型提取出多类别随机噪声信号,实现地震数据的随机信号压制。实验结果显示,这种改进小波方法对地震数据随机噪声信号的去噪效果较好,在复杂沉积地质结构被探测介质的地震数据随机噪声压制方面具有较强的适用性。     

REDUCTION OF HARMONIC DISTORTION IN VIBRATORY SOURCE RECORDS*

Due to non-linear effects, the swept frequency signals (sweeps) transmitted into the subsurface by vibrators are contaminated by harmonics. Upon correlation of the recorded seismograms, these harmonics lead to noise trains which are particularly disturbing in the case of down-sweeps. The method described in this paper—which can be regarded as a generalization of Sorkin's approach to the suppression of even order harmonics—allows elimination, from the final vibratory source seismogram, of harmonics of the sweep up to any desired order. It requires that not one single signal but rather a series of M signals is employed where each signal has an initial phase differing from that of the previous one of the series by the phase angle 2πM. Prior to stacking, the seismograms generated with the different signals have to be brought into the form they would have if they had been generated with the same signal. The method seems also to be capable of reducing the correlation noise if sign-bit recording techniques are used.