共查询到20条相似文献，搜索用时 125 毫秒
1.
本文提出了在τ-P域实现三维叠前深度偏移的计算方法，该方法把地震道集由时-空域变换到τ-P域，然后对共P数据体进行相移或相移加插值偏移，将偏移后的各共P数据体叠加得到叠前深度偏移结果．该算法能够实现陡倾地层的正确归位，降低偏移处理维数，减少计算机内存需求，易于并行处理，但是计算量和所花费的计算机时间仍然很大．为此我们又提出了针对目的层的分时（层）偏移方法，由于只考虑目的层，可节省大量的计算时间．文中以二维实际地震资料及三维数值模型资料的处理说明了上述方法的可行性与有效性． 相似文献
2.
3.
本文介绍一种基于τ-P变换的新的叠后深度偏移方法。这一方法是在均匀介质条件下根据波动方程推导出的。其基本思想是:在τ-p域中,当波场分解为它的平面波成分时,改变其τ、p值,再通过τ-p反变换回到深度域中,从而完成了地下界面的偏移归位。 方法所需的计算时间基本是与完成τ-p正、反变换所需时间相同,是一种快速的偏移方法。它可用于任意形状的地下界面,甚至在强干扰背景的条件下,也能提高资料的信噪比和改善信号的连续性。这一方法是在深度域中进行τ-p变换,故它能适应速度存在纵、横变化的地区。为了论证本方法的效果,本文给出了不同的正演理论模型和偏移结果的例子。 相似文献
4.
本文对用爆炸地震资料研究地壳与上地幔结构及速度分布的τ法反演进行了探讨,并且在Bessonova提出的平均值方法的基础上,就减少解的不确定性及简化反演过程提出了新的作法:1.τ（p）的上下界各自独立地参与反演而后再将所得到的结果合并起来;2.在τ法反演中,抽去地壳模型中的低速层。根据以上两种计算方案,改进方法,取得较好的效果。 本文还用Bessonova的平均值方法计算了我国东部地区地壳-上地幔爆炸地震探测剖面的实例。该地区地壳厚度约为31 km,在18 km深处存在着低速层,厚度约6 km,速度为6 km/s。 相似文献
5.
《地球物理学进展》2017,(2)
使用常规的Wiener反褶积必须假设震源子波在地层旅行过程中是平稳的即一成不变的,这个前提条件与实际野外地震资料采集差别较大,而基于Gabor变换反褶积技术考虑到地震能量的衰减、子波的形变等非平稳性特征.地震道在Gabor域可因式分解成三项即震源子波、衰减函数和反射系数,该技术设计POU窗函数,并利用此函数在Gabor域对地震信号进行局部时频分解.Gabor域反褶积算法在Gabor域通过除以衰减函数和震源子波的乘积来估算地层反射系数,然后再做Gabor反变换可求得时间域的地层反射系数.理论模型的测试和实际地震资料的应用均表明,与Wiener反褶积相比较,基于Gabor变换反褶积可补偿中深层的能量衰减并因此拓宽有效频带和提高时间分辨率. 相似文献
6.
《应用地球物理》2019,(4)
叠前深度偏移的广泛应用产生了大量的深度域地震数据,提供了深度域地震处理的数据基础,并为探索直接基于深度域地震数据的储层弹性参数反演提供了可能性。由于深度域地震子波存在与速度密切相关的特征,直接提取深度域地震子波存在困难,本文将深度域速度变换与脊回归方法结合,实现了深度域地震子波的确定性提取,将提取的深度域常速地震子波构建成依赖速度的空变地震子波,基于加权叠加思想,提出了深度域地震记录合成方法。再结合非线性共轭梯度算法,建立了基于深度域地震数据和模型驱动的深度域叠前波阻抗反演流程,不再进行额外的时间域和深度域之间的数据转换。深度域合成地震数据与实际地震资料的叠前波阻抗反演取得了较好的结果,表明该方法是可行的。该方法可作为深度域正演分析和叠前弹性参数反演的基础方法,为深度域反演方法的稳定性、精度和效率的进一步研究创造了条件。 相似文献
7.
从导出的地震定标律和地震破裂过程的断裂力学模式出发,得到了震源平均位移、速度和加速度谱的表达式,进而又推导出震源的峰值位移d_m,速度v_m和加速度a_m的表达式:d_m=k_dM_0~2/~3τ_0~(2/3),v=k_vM_0~(1/3)τ_0~(4/3),a_m=k_aτ_0~2式中M_0是地震矩,τ_0是环境剪应力值,k_d,k_v,k_a为适当的常数.我们选用了66个地震的观测资料,这些地震的矩震级范围包括了从1级左右的极微震,3-5级的小震,直到6-7级的大震;地震矩从10~9-10~(20)Nm,跨越了10个数量级,并用这些地震检验了上述公式. 令所有地震的平均应力为5MPa,定出常数k_d,k_v,k_a,进而由速度和加速度观测资料求得66个地震的环境剪应力τ_0值,这些数值相当稳定.多数极微震的τ_0值在2-4MPa之间;小震的τ_0值多数在4-8MPa左右;大震的τ_0值为10MPa左右。τ_0值对震源深度和断层类型有明显的依赖性.一般深度很浅的小震和极微震,τ_0值很低;正断层地震的τ_0值相对较低;逆断层地震的τ_0值较高;走滑断层地震的τ_0值则居中. 相似文献
8.
《地球物理学进展》2015,(6)
反褶积是叠前地震数据处理中的重要环节,反褶积效果的好坏很大程度上依赖于地震子波的准确性.早期的反褶积处理大多数都是基于Robinson提出的平稳褶积模型,即地震子波是时不变的,但实际上由于地下介质的吸收衰减作用,地震子波是随时间不断变化的,这说明要进一步改善反褶积,使用时变的地震子波是必要的,因此本文提出了一种从地震资料中直接提取时变子波的方法.具体地,首先对单道地震数据做S变换求出其时频谱,进而得到其时变振幅谱,然后利用谱模拟技术从求得的地震记录振幅谱中拟合出每一时刻的子波振幅谱,在子波是零相位假设的前提下,拟合出的时变子波振幅谱即是所求频率域的时变子波,本文最后利用正演的单道地震记录和实际资料分别验证了所求频率域时变子波的准确性. 相似文献
9.
地震预警作为一种能够有效减轻地震灾害的手段已经被世界上越来越多的国家和地区所研究,并在实际应用中取得了显著的减灾实效.利用汶川主震及其余震P波和S波初期部分的信息,研究了最终地震震级与4个预警参数的相关性.考察的4个参数是位移幅值(Pd)、速度平方积分(IV2)、P波卓越周期和特征周期(τp和τc).使用的数据集时间跨度从2008年5月12日至2008年10月4日,共计218次震级大于等于4级的地震事件,包括主震8.0、最大余震6.5和7次大于等于6.0级的事件.Pd、IV2和τc在震级小于等于6.5级时与最终震级具有较好的相关性,没有出现震级饱和现象,验证了前人的统计结果.但是在估算主震震级时,都存在不同程度的低估现象,IV2参数尤其明显.另外,Pd和IV2统计曲线在M6级附近存在着一个斜率变化,并且Pd结果与前人统计结果相差较大.τp的统计结果表明该参数不适合应用到汶川地区的地震预警系统中. 相似文献
10.
《地球物理学进展》2020,(3)
全波形反演具有高精度成像能力,然而由理论走向实际应用还存在很多问题.全波形反演通过匹配波形来更新模型,数据的波形与地震子波有直接的关系.本文介绍了一种针对时间域波形反演的子波估计方法,并将其应用到全波形反演中.由于在频率域子波反演可表示为一个线性优化问题,因此本文先给定一个试探子波,在时间域通过有限差分法正演得到地震波场,并在频率域与观测数据比较分析,反演得到预测子波.此外,本文还简要地介绍了如何使用伴随状态法计算全波形反演的梯度.数值实验证明,本文方法反演得到地震子波与真实子波具有很好的吻合度,并在全波形反演中取得了不错的效果.相比不依赖于子波的方法,提前进行子波估计可提高反演效率.此外,本文方法适用于全波形反演中常使用的多尺度反演策略. 相似文献
11.
Qingju Wu Xiufen Zheng Jiatie Pan Fengxue Zhang Guangcheng Zhang Institute of Geophysics China Earthquake Administration Beijing China 《地震学报(英文版)》2009,(4):425-429
In this paper,we present wavelet transformation method to measure interstation phase velocity. We use Morlet wavelet function as mother wavelet to filter two seismograms at various period of interest,and correlate the wavelet filtered seismograms to form cross-correlogram. If both wavelet filtered signals are in phase at that period,the phase of the cross-correlogram is a minimum. Using 3-spline interpolation to transform cross-correlation matrix to a phase velocity verse period image,it is convenient for u... 相似文献
12.
13.
MONIKA BREITZKE 《Geophysical Prospecting》1992,40(1):31-70
In spite of a geometrical rotation into radial and transverse parts, two- or three-component in-seam seismic data used for underground fault detection often suffer from the problem of overmoding ‘noise’. Special recompression filters are required to remove this multimode dispersion so that conventional reflection seismic data processing methods, e.g. CMP stacking techniques, can be applied afterwards. A normal-mode superposition approach is used to design such multimode recompression filters. Based on the determination of the Green's function in the far-field, the normal-mode superposition approach is usually used for the computation of synthetic single- and multi-mode (transmission) seismograms for vertically layered media. From the filter theory's point of view these Green's functions can be considered as dispersion filters which are convolved with a source wavelet to produce the synthetic seismograms. Thus, the design of multimode recompression filters can be reduced to a determination of the inverse of the Green's function. Two methods are introduced to derive these inverse filters. The first operates in the frequency domain and is based on the amplitude and phase spectrum of the Green's function. The second starts with the Green's function in the time domain and calculates two-sided recursive filters. To test the performance of the normal-mode superposition approach for in-seam seismic problems, it is first compared and applied to synthetic finite-difference seismograms of the Love-type which include a complete solution of the wave equation. It becomes obvious that in the case of one and two superposing normal modes, the synthetic Love seam-wave seismograms based on the normal-mode superposition approach agree exactly with the finite-difference data if the travel distance exceeds two dominant wavelengths. Similarly, the application of the one- and two-mode recompression filters to the finite-difference data results in an almost perfect reconstruction of the source wavelet already two dominant wavelengths away from the source. Subsequently, based on the dispersion analysis of an in-seam seismic transmission survey, the normal-mode superposition approach is used both to compute one- and multi-mode synthetic seismograms and to apply one- and multimode recompression filters to the field data. The comparison of the one- and two-mode synthetic seismograms with the in-seam seismic transmission data reveals that arrival times, duration and shape of the wavegroups and their relative excitation strengths could well be modelled by the normal-mode superposition approach. The one-mode recompressions of the transmission seismograms result in non-dispersive wavelets whose temporal resolution and signal-to-noise ratio could clearly be improved. The simultaneous two-mode recompressions of the underground transmission data show that, probably due to band-limitation, the dispersion characteristics of the single modes could not be evaluated sufficiently accurately from the field data in the high-frequency range. Additional techniques which overcome the problem of band-limitation by modelling all of the enclosed single-mode dispersion characteristics up to the Nyquist frequency will be mandatory for future multimode applications. 相似文献
14.
15.
P. H. NIELSEN 《Geophysical Prospecting》1978,26(2):399-406
Calculation of reflection seismograms in the frequency domain makes it possible to treat absorption correctly, i.e. makes it possible to include dispersion (= causality) and frequency dependent reflection coefficients. In the paper it is shown how such a calculation can be done in a reasonably efficient manner. 相似文献
16.
The authors proposed a method for obtaining high-quality acceleration seismograms from velocity type seismograms of digital Seismographic network, and took as an example the analysis and processing of the seismograms of a same earthquake that was simultaneously recorded by velocity seismograph CTS1-EDAS24 and strong motion seismograph EST-Q4128 installed in Jixian Station, Tianjin. The calculation steps and the processing method have been discussed in detail. From the analysis and the comparison of the obtained results, it is concluded that the proposed method is simple and effective, and it broadens the application of digital seismographic network. 相似文献
17.
Shallow SH-wave reflections are far from routine, although their study can provide insights into important properties of near-surface materials that cannot be inferred from P-wave data alone. Difficulties in separating SH-wave reflections from Love waves are generally considered the major obstacle to progress in shallow SH-wave seismic reflection. This may be the case in surveys undertaken at great depths, but it is not necessarily true for reflection data gathered at shallow and ultra-shallow depths. This paper shows that when SH-wave data possess wavelengths greater than the thickness of the superficial layer, Love waves are not greatly dispersed. In this case, misinterpretation between parts of reflection hyperbolae and waveguide arrivals is sufficiently limited. In a one-layer model earth, which well approximates typical situations of the near-surface underground, the most energetic modes (the lowermost modes) of the dispersed surface waves have a dominant frequency band that falls below the wavelet spectrum of the shallow reflections; therefore, they can be filtered out in the frequency domain. Higher modes, although their spectral content overlaps that of the reflections, exhibit small amplitudes on seismograms and leave strong reflections unaffected.We present field examples from three different sites where we were able to obtain ultra-shallow reflections (< 3 m) in unconsolidated sediments. The high level of resolution (vertical resolution up to 15 cm) suggests that SH-wave reflection imaging has the potential to complement other high-resolution techniques, such as P-wave reflection and ground-penetrating radar (GPR) imaging, allowing a better and more complete characterization of the near-surface environments. 相似文献
18.
For years, reflection coefficients have been the main aim of traditional deconvolution methods for their significant informational content. A method to estimate seismic reflection coefficients has been derived by searching for their amplitude and their time positions without any other limitating assumption. The input data have to satisfy certain quality constraints like amplitude and almost zero phase noise—ghosts, reverberations, long period multiples, and diffracted waves should be rejected by traditional processing. The proposed algorithm minimizes a functional of the difference between the spectra of trace and reflectivity in the frequency domain. The estimation of reflection coefficients together with the consistent “wavelet’ is reached iteratively with a multidimensional Newton-Raphson technique. The residual error trace shows the behavior of the process. Several advantages are then obtainable from these reflection coefficients, like conversion to interval velocities with an optimum calibration either to the well logs or to the velocity analysis curves. The procedure can be applied for detailed stratigraphic interpretations or to improve the resolution of a conventional velocity analysis. 相似文献
19.
In this paper, a novel data denoising method is proposed for seismic exploration with a vibrator which produces a chirp-like signal. The method is based on fractional wavelet transform (FRWT), which is similar to the fractional Fourier transform (FRFT). It can represent signals in the fractional domain, and has the advantages of multi-resolution analysis as the wavelet transform (WT). The fractional wavelet transform can process the reflective chirp signal as pulse seismic signal and decompose it into multi-resolution domain to denoise. Compared with other methods, FRWT can offer wavelet transform for signal analysis in the timefractional-frequency plane which is suitable for processing vibratory seismic data. It can not only achieve better denoising performance, but also improve the quality and continuity of the reflection syncphase axis. 相似文献
20.
R.-G. FERBER 《Geophysical Prospecting》1985,33(2):212-223
A main problem in computing reflection coefficients from seismograms is the instability of the inversion procedure due to noise. This problem is attacked for two well-known inversion schemes for normal-incidence reflection seismograms. The crustal model consists of a stack of elastic, laterally homogeneous layers between two elastic half-spaces. The first method, which directly computes the reflection coefficients from the seismogram is called “Dynamic Deconvolution”. The second method, here called “Inversion Filtering”, is a two-stage procedure. The first stage is the construction of a causal filter by factorization of the spectral function via Levinson-recursion. Filtering the seismogram is the second stage. The filtered seismogram is a good approximation for the reflection coefficients sequence (unless the coefficients are too large). In the non-linear terms of dynamic deconvolution and Levinson-recursion the noise could play havoc with the computation. In order to stabilize the algorithms, the bias of these terms is estimated and removed. Additionally incorporated is a statistical test for the reflection coefficients in dynamic deconvolution and the partial correlation coefficients in Levinson-recursion, which are set to zero if they are not significantly different from noise. The result of stabilization is demonstrated on synthetic seismograms. For unit spike source pulse and white noise, dynamic deconvolution outperforms inversion filtering due to its exact nature and lesser computational burden. On the other hand, especially in the more realistic bandlimited case, inversion filtering has the great advantage that the second stage acts linearly on the seismogram, which allows the calculation of the effect of the inversion procedure on the wavelet shape and the noise spectrum. 相似文献