共查询到20条相似文献，搜索用时 15 毫秒
1.
We present a new method for the extraction and removal of the source wavelet from the reflection seismogram. In contrast to all other methods currently in use, this one does not demand that there be any mathematically convenient relationship between the phase spectrum of the source wavelet and the phase spectrum of the earth impulse response. Instead, it requires a fundamental change in the field technique such that two different seismograms are now generated from each source-receiver pair: the source and receiver locations stay the same, but the source used to generate one seismogram is a scaled version of the source used to generate the other. A scaling law provides the relationship between the two source signatures and permits the earth impulse response to be extracted from the seismograms without any of the usual assumptions about phase. We derive the scaling law for point sources in an homogeneous isotropic medium. Next, we describe a method for the solution of the set of three simultaneous equations and test it rigorously using a variety of synthetic data and two types of synthetic source waveform: damped sine waves and non-minimum-phase air gun waveforms. Finally we demonstrate that this method is stable in the presence of noise. 相似文献
2.
A.M. ZIOLKOWSKI 《Geophysical Prospecting》1986,34(3):383-408
The article provides a theoretical basis for the extension of the method of scaling law deconvolution to three dimensions using airgun arrays as a sound source. Earlier papers by the author required the dimensions of the scaled sources to be different while the depths and firing pressures were maintained the same in order to preserve the same dynamics of the scaled sources at scaled time. However, this forces the source ghost to be considered as part of the impulse response of the earth rather than as part of the downgoing source wave. And, in fact, the dynamics of the scaled sources are not the same at the same depth because the ghost reflection modulates the behaviour of the oscillating bubbles generated by the airguns, and this modulation does not scale. To force the sources to scale properly, including the ghost interaction, the larger source must be put at greater depth, where hydrostatic pressure is greater, and the initial firing pressure must be adjusted accordingly. Thus, the depth, initial firing pressure and gun volume are all variables. The interaction among guns in scaled airgun arrays also scales exactly if the geometry of an array and the depth of its deployment are scaled by the same factor. 相似文献
3.
本文首先分析了地震波在黏弹介质的传播规律,基于黏弹介质地震波动方程总结了时变子波振幅谱和相位谱的关系,从而得出结论,准确估计子波相位谱初值和不同时刻的子波振幅谱是实现时变子波准确提取的必要条件.在此基础上,针对传统方法限制子波振幅谱形态且受限于分段平稳假设的问题,提出了一种利用EMD(Empirical Mode Decomposition)和子波振幅谱与相位谱关系的时变子波提取方法,根据子波对数振幅谱光滑连续而反射系数对数振幅谱振荡剧烈的特点,采用EMD方法将不同时刻地震记录的对数振幅谱分解为一组具有不同振荡尺度的模态分量,通过滤除振荡剧烈分量、重构光滑连续分量提取时变子波振幅谱;再应用子波振幅谱和相位谱的关系提取时变子波相位谱,将分别提取的振幅谱和相位谱逐点进行合成,最终实现时变子波的准确提取.本文方法不需要求取Q值,适用于变Q值的情况,具有良好的抗噪性能.数值仿真和叠后实际资料处理结果表明,相比传统的分段提取方法,利用本文方法提取的时变子波准确度更高,研究成果对提高地震资料分辨率具有重要意义. 相似文献
4.
Air-gun arrays are used in marine-seismic exploration. Far-field wavelets in subsurface media represent the stacking of single air-gun ideal wavelets. We derived single air-gun ideal wavelets using near-field wavelets recorded from near-field geophones and then synthesized them into far-field wavelets. This is critical for processing wavelets in marine- seismic exploration. For this purpose, several algorithms are currently used to decompose and synthesize wavelets in the time domain. If the traveltime of single air-gun wavelets is not an integral multiple of the sampling interval, the complex and error-prone resampling of the seismic signals using the time-domain method is necessary. Based on the relation between the frequency-domain phase and the time-domain time delay, we propose a method that first transforms the real near-field wavelet to the frequency domain via Fourier transforms; then, it decomposes it and composes the wavelet spectrum in the frequency domain, and then back transforms it to the time domain. Thus, the resampling problem is avoided and single air-gun wavelets and far-field wavelets can be reliably derived. The effect of ghost reflections is also considered, while decomposing the wavelet and removing the ghost reflections. Modeling and real data processing were used to demonstrate the feasibility of the proposed method. 相似文献
5.
M.M. ROKSANDIý 《Geophysical Prospecting》1986,34(8):1154-1166
Two factors are responsible for the fact that an extended marine source array performs better than a point source: 1. a higher degree of transmission of the radiated seismic energy through the water-sediment interface providing a better penetration; 2. filtering effects. The higher degree of transmission is due to: (a) the directivity of extended sources, (b) the lower reflection coefficient at the water-sediment interface for seismic waves radiated from an extended source array than for spherical seismic waves radiated from a point source, (c) the lower amplitude decay of the pulses from an extended source than from a point source. In addition, signature characteristic of an extended source array and Fresnel zone of waves generated by such a source differ from those corresponding to a point source. The propagating wavelet radiated from a point source array may not be, in a sedimentary sequence below the sea-floor, the linear combination of wavelets emitted from point sources. In such cases, there is a noticeable difference between the performance of a field-implemented source array and that of the corresponding simulated source array. The performances of the field-implemented and simulated extended receiver arrays can be identical if the recording system is adequate and the processing technique appropriate. 相似文献
6.
为研究地震子波相位对反射系数序列反演的影响,在自回归滑动平均(ARMA)模型描述子波的基础上,提出采用z域对称映射ARMA模型零极点的方法构造了一系列相同振幅谱、不同相位谱的地震子波,并结合谱除法对人工合成地震记录进行反射系数序列反演.理论分析表明,子波相位估计不准时反射系数序列反演结果中残留一个纯相位滤波器,该纯相位滤波器的相位谱为真实子波和构造子波的相位谱之差.采用丰度和变分作为评价方法,在反演结果中确定出真实的或准确的反射系数序列.仿真实验和实际数据处理结果也验证了子波相位对反射系数序列反演的影响规律和评价方法的有效性,为进一步提高反射系数序列反演结果精度指明了研究方向. 相似文献
7.
小波变换与地震信号特征分析 总被引：6，自引：0，他引：6
地震波的瞬时信号,如瞬时振幅、瞬时频率及瞬时相位等,是研究地球介质的重要参数。根据小波定义,对Morlet 小波进行了修正,并对修正后的小波形态进行了深入讨论。理论分析表明,小波变换效果受到整形参数、小波长度、中心频率、频带宽度及小波个数等参数的制约,特别是整形参数与小波中心频率及频带之间关系对小波变换起到决定性作用。在地震波信号的实际处理中,可选取恰当的整形参数,同时采用合适的小波中心频率以避免小波变换对信号产生的遗漏和冗余。文中给出了实际地震记录处理的示例。 相似文献
8.
利用中、小容量气枪组成的立体气枪阵列延迟激发震源和海底地震仪(OBS)在我国北部浅海海域开展了人工地震深部地球物理探测试验.基于水深条件和压制水体虚反射、提升低频能量的需要,使气枪震源有足够的输出能量和高品质子波特性,研究了立体气枪阵列延迟激发震源工作机理,经远场子波理论模拟优选了组合参数并进行了海上试验工作.结果表明,中、小容量气枪组成的立体气枪阵列延迟激发震源,适应了浅水海域的激发环境,降低了由虚反射造成的局部陷波和干扰作用,有效地改善了OBS信号的品质,获得了Ps,Pg,PmP,Pn等多种震相.创新了由中、小容量气枪组成的立体气枪阵列延迟激发震源在浅海区OBS探测中的应用,也填补了南黄海海域深地震探测数据的空白,为南黄海、渤海深部地壳结构研究及含油气盆地形成演化研究提供了重要的基础资料. 相似文献
9.
Deconvolution is an essential step for high-resolution imaging in seismic data processing. The frequency and phase of the seismic wavelet change through time during wave propagation as a consequence of seismic absorption. Therefore, wavelet estimation is the most vital step of deconvolution, which plays the main role in seismic processing and inversion. Gabor deconvolution is an effective method to eliminate attenuation effects. Since Gabor transform does not prepare the information about the phase, minimum-phase assumption is usually supposed to estimate the phase of the wavelet. This manner does not return the optimum response where the source wavelet would be dominantly a mixed phase. We used the kurtosis maximization algorithm to estimate the phase of the wavelet. First, we removed the attenuation effect in the Gabor domain and computed the amplitude spectrum of the source wavelet; then, we rotated the seismic trace with a constant phase to reach the maximum kurtosis. This procedure was repeated in moving windows to obtain the time-varying phase changes. After that, the propagating wavelet was generated to solve the inversion problem of the convolutional model. We showed that the assumption of minimum phase does not reflect a suitable response in the case of mixed-phase wavelets. Application of this algorithm on synthetic and real data shows that subtle reflectivity information could be recovered and vertical seismic resolution is significantly improved. 相似文献
10.
A new approach to deconvolution has been developed to improve the attenuation of multiple energy. This approach to deconvolution is unique in that it not only eliminates the usual assumptions of a minimum phase lag wavelet and a random distribution of impulses, but also overcomes the noise limitation of the homomorphic deconvolution and its inherent instability to phase computation. We attempt to analyse the continuous alteration of the acoustic waveform during the propagation through a linear system. Based on the results of this analysis, the surface-related measurements are described as a convolution of the impulse response of the system with the non-stationary forward wavelet which includes all multiple terms generated within the system. The amplitude spectrum of the forward wavelet is recovered from the amplitude spectrum of the recorded signal, using the difference between the rate of decay of the source wavelet and the duration of the measurement. The phase spectrum of the forward wavelet is estimated using the Hilbert transform and the fact that the mixed phase lag wavelet can be presented as a convolution of the minimum and maximum phase lag wavelets. The multiples are discriminated from primaries by comparison of the phase spectrum of the seismic signal and the inverse of the forward wavelet. Therefore, the technique is called phase inversion deconvolution (PID). This approach requires no velocity information in order to recognize and attenuate multiple energy. Therefore, primary energy is recovered in the near-offset region where the velocity differential between primary and multiple energies is very small. 相似文献
11.
A. N. BOWEN 《Geophysical Prospecting》1986,34(3):366-382
Wiener ‘spiking’ deconvolution of seismic traces in the absence of a known source wavelet relies upon the use of digital filters, which are optimum in a least-squares error sense only if the wavelet to be deconvolved is minimum phase. In the marine environment in particular this condition is frequently violated, since bubble pulse oscillations result in source signatures which deviate significantly from minimum phase. The degree to which the deconvolution is impaired by such violation is generally difficult to assess, since without a measured source signature there is no optimally deconvolved trace with which the spiked trace may be compared. A recently developed near-bottom seismic profiler used in conjunction with a surface air gun source produces traces which contain the far-field source signature as the first arrival. Knowledge of this characteristic wavelet permits the design of two-sided Wiener spiking and shaping filters which can be used to accurately deconvolve the remainder of the trace. In this paper the performance of such optimum-lag filters is compared with that of the zero-lag (one-sided) operators which can be evaluated from the reflected arrival sequence alone by assuming a minimum phase source wavelet. Results indicate that the use of zero-lag operators on traces containing non-minimum phase wavelets introduces significant quantities of noise energy into the seismic record. Signal to noise ratios may however be preserved or even increased during deconvolution by the use of optimum-lag spiking or shaping filters. A debubbling technique involving matched filtering of the trace with the source wavelet followed by optimum-lag Wiener deconvolution did not give a higher quality result than can be obtained simply by the application of a suitably chosen Wiener shaping filter. However, cross correlation of an optimum-lag spike filtered trace with the known ‘actual output’ of the filter when presented with the source signature is found to enhance signal-to-noise ratio whilst maintaining improved resolution. 相似文献
12.
Linear and especially non-linear analyses of spatially extended structures, such as pipelines and bridges, often requires specification of time histories of ground motion at an array of closely spaced points. As the number of dense accelerograph arrays worldwide is small, and the number of earthquake observations is limited, synthetic motions with desired characteristics become necessary. This paper presents a method for synthesizing such motions, which is an extension of the SYNACC method, developed first in the early 1970s for synthetic accelerations, velocities and displacements at a point, and later extended to synthetic near surface strains, rotations and curvatures of ground motion at a point. It consists of unfolding in time a site specific Fourier amplitude spectrum of ground acceleration, obtained by an empirical scaling model, by representing the ground motion as a superposition of traveling wavelets of Love and Rayleigh waves and body waves, which propagate with phase and group velocities consistent with the dispersion characteristic of the site geology, approximated by parallel layers. Uniform hazard Fourier spectra or any specified target Fourier spectrum can also be used. Derivations of the point strains, rotations and curvatures are also presented. The method is illustrated for scenario M6.5 and M7.5 earthquakes and three dispersion models. 相似文献
13.
The amplitude spectrum of ground penetrating radar (GPR) reflection data acquired with a particular antenna set is normally concentrated over a spectral bandwidth of a single octave, limiting the resolving power of the GPR wavelet. Where variously-sized GPR targets are located at numerous depths in the ground, it is often necessary to acquire several profiles of GPR data using antennas of different nominal frequencies. The most complete understanding of the subsurface is obtained when those frequency-limited radargrams are jointly interpreted, since each frequency yields a particular response to subsurface reflectivity. The application of deconvolution to GPR data could improve image quality, but is often hindered by limited spectral bandwidth.We present multiple-frequency compositing as a means of combining data from several frequency-limited datasets and improving the spectral bandwidth of the GPR profile. A multiple-frequency composite is built by summing together a number of spatially-coincident radargrams, each acquired with antennae of different centre frequency. The goal of the compositing process is therefore to produce a composite radargram with balanced contributions from frequency-limited radargrams and obtain a composite wavelet that has properties approximating a delta function (i.e. short in duration and having a broad, uniform spectral bandwidth).A synthetic investigation of the compositing process was performed using Berlage wavelets as proxies for GPR source pulses. This investigation suggests that a balanced, broad bandwidth, effective source pulse is obtained by a compositing process that equalises the spectral maxima of frequency-limited wavelets prior to summation into the composite. The compositing of real GPR data was examined using a set of 225, 450 and 900 MHz GPR common offset profiles acquired at a site on the Waterloo Moraine in Ontario, Canada. The most successful compositing strategy involved derivation of scaling factors from a time-variant least squares analysis of the amplitude spectra of each frequency-limited dataset. Contributions to the composite from each nominal acquisition frequency are clear, and the trace averaged amplitude spectrum of the corresponding composite is broadened uniformly over a bandwidth approaching two-octaves. Improvements to wavelet resolution are clear when a composite radargram is treated with a spiking deconvolution algorithm. Such improvement suggests that multiple-frequency compositing is a useful imaging tool, and a promising foundation for improving deconvolution of GPR data. 相似文献
14.
Time‐lapse seismic surveying has become an accepted tool for reservoir monitoring applications, thus placing a high premium on data repeatability. One factor affecting data repeatability is the influence of the rough sea‐surface on the ghost reflection and the resulting seismic wavelets of the sources and receivers. During data analysis, the sea‐surface is normally assumed to be stationary and, indeed, to be flat. The non‐flatness of the sea‐surface introduces amplitude and phase perturbations to the source and receiver responses and these can affect the time‐lapse image. We simulated the influence of rough sea‐surfaces on seismic data acquisition. For a typical seismic line with a 48‐fold stack, a 2‐m significant‐wave‐height sea introduces RMS errors of about 5–10% into the stacked data. This level of error is probably not important for structural imaging but could be significant for time‐lapse surveying when the expected difference anomaly is small. The errors are distributed differently for sources and receivers because of the different ways they are towed. Furthermore, the source wavelet is determined by the sea shape at the moment the shot is fired, whereas the receiver wavelet is time‐varying because the sea moves significantly during the seismic record. 相似文献
15.
Klauder wavelet removal before vibroseis deconvolution 总被引：1，自引：0，他引：1
The spiking deconvolution of a field seismic trace requires that the seismic wavelet on the trace be minimum phase. On a dynamite trace, the component wavelets due to the effects of recording instruments, coupling, attenuation, ghosts, reverberations and other types of multiple reflection are minimum phase. The seismic wavelet is the convolution of the component wavelets. As a result, the seismic wavelet on a dynamite trace is minimum phase and thus can be removed by spiking deconvolution. However, on a correlated vibroseis trace, the seismic wavelet is the convolution of the zero-phase Klauder wavelet with the component minimum-phase wavelets. Thus the seismic wavelet occurring on a correlated vibroseis trace does not meet the minimum-phase requirement necessary for spiking deconvolution, and the final result of deconvolution is less than optimal. Over the years, this problem has been investigated and various methods of correction have been introduced. In essence, the existing methods of vibroseis deconvolution make use of a correction that converts (on the correlated trace) the Klauder wavelet into its minimum-phase counterpart. The seismic wavelet, which is the convolution of the minimum-phase counterpart with the component minimum-phase wavelets, is then removed by spiking deconvolution. This means that spiking deconvolution removes both the constructed minimum-phase Klauder counterpart and the component minimum-phase wavelets. Here, a new method is proposed: instead of being converted to minimum phase, the Klauder wavelet is removed directly. The spiking deconvolution can then proceed unimpeded as in the case of a dynamite record. These results also hold for gap predictive deconvolution because gap deconvolution is a special case of spiking deconvolution in which the deconvolved trace is smoothed by the front part of the minimum-phase wavelet that was removed. 相似文献
16.
Mihailo D. Trifunac 《地震工程与结构动力学》1973,2(1):59-69
Prediction of response spectra for earthquake engineering purposes is considered from a new point of view based on the dislocation theory of earthquakes. It is shown that the traditional scaling of response spectra by the predicted peak acceleration should be limited to the high-frequency end of the spectrum, and that the peak acceleration in the near field is not strongly correlated with earthquake magnitude. The amplitude of the long-period end of response spectrum at source to station distances greater than about 10 source dimensions should be scaled with seismic moment, while for distances less than about one source dimension this amplitude should be proportional to the permanent ground displacement. To reconcile the existing extensive data on seismicity of active regions based on magnitude scale, it is shown that magnitude can be used to determine approximately the seismic moment. 相似文献
17.
为了对比天然源面波勘探不同台阵布局的探测效果, 筛选出探测成果可靠、 效率高和便于野外施工的天然源面波勘探台阵阵形, 在天水市黄土覆盖区的同一场地分别用4种常见的阵形进行数据采集试验, 并对各种阵形数据使用空间自相关法或扩展空间自相关法提取相应的频散曲线, 通过反演得到了试验点地下的浅层速度结构模型. 分析对比试验结果表明: 4种台阵提取的频散曲线数值很相近; 频散谱能量集中度较高的是嵌套式等边三角形和圆形台阵, L形和直线形台阵相对分散; L形台阵低频段(4—8 Hz)比直线形台阵差, 其高频段(8—40 Hz)比直线形台阵好. 针对直线形台阵在高频段信噪比较低的情况, 在确保探测成果可靠性的前提下, 为了提高探测效率, 提出了在同一直线形台阵开展天然源与人工源面波联合勘探的数据采集方法. 实验结果证实, 这种联合勘探方法不仅可弥补直线形台阵高频段的不足, 确保探测精度和结果的可靠性, 而且还能实现“高低”频兼顾, 即“深浅”兼顾. 相似文献
18.
D. A. Yuen G. Erlebacher O. V. Vasilyev D. E. Goldstein M. Fuentes 《Pure and Applied Geophysics》2004,161(11-12):2231-2244
— Today, wavelets are recognized to have a wide range of useful properties that allow them to treat effectively multifacet problems, such as data compression, scale-localization analysis, feature extraction, statistics, numerical simulation, visualization, and communication. Second-generation wavelets represent a recent improvement of the wavelet algorithm, allowing for greater flexibility in the spatial domain and other computational advantages. We will show how these wavelets can be employed to extract large-scale coherent structures from (1) three-dimensional turbulent flows and (2) high Rayleigh number thermal convection. We will discuss the concept of modelling via decomposition into coherent and incoherent fields, taking into account the effect of the incoherent field via statistical modelling. We will discuss wavelet properties and how they can be utilized and integrated in handling large-scale problems in earthquake physics and other nonlinear phenomena in the geosciences. 相似文献
19.
Sándor Szalai László Szarka 《Journal of Applied Geophysics》2011,75(1):1-8
Recent two-dimensional multielectrode measurements are restricted to only a few geoelectric arrays. Realizing that specific features of nearly 90 other arrays are totally ignored, all original intentions as published about the development of new geoelectric arrays were reviewed. Apart from arrays, either already applied in two-dimensional geoelectric arrays or impossible to be applied in such systems, 61 forgotten once-developed arrays were found. These provide altogether 102 various solutions, which would be able to increase the efficiency of two-dimensional multielectrode measurements in some respect. 46 array solutions are able to enhance the depth of investigation; 9/11 array solutions give better vertical/horizontal resolution; 17 array solutions provide better planview images; 8 array solutions are worth applying in areas with limited access; 11 array solutions may reduce the effect of near-surface inhomogeneities. By reviving these forgotten arrays, it will be possible to develop versatile multielectrode systems, which are more adaptive to the diverse field needs. 相似文献
20.
Selecting, scaling and matching accelerograms are critically important to engineering design and assessment, enabling structural response to be determined with greater confidence and through fewer analyses than if unscaled accelerograms are employed. This paper considers the response of an 8‐storey multiple‐degree‐of‐freedom reinforced concrete structure to accelerograms selected, linearly scaled or spectrally matched using five different techniques. The first method consists of selecting real records on the basis of seismological characteristics, while the remaining methods make an initial selection on the basis of magnitude and spectral shape before (1) scaling to the target spectral acceleration at the initial period; (2) scaling to the target spectrum over a range of periods; (3) using wavelet adjustments to match the target spectrum and (4) using wavelet adjustments to match multiple target spectra for multiple damping ratios. The analyses indicate that the number of records required to obtain a stable estimate of the response decreases drastically as one moves through these methods. The exact number varies among damage measures and is related to the predictability of the damage measure. For measures such as peak roof and inter‐storey drift, member end rotation and the Park and Ang damage index, as few as one or two records are required to estimate the response to within ±5% (for a 64% confidence level) if matching to multiple damping ratios is conducted. Bias checks are made using predictive equations of the expected response derived from the results of 1656 nonlinear time‐domain analyses of the structure under the action of unscaled accelerograms. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献