共查询到20条相似文献，搜索用时 437 毫秒
1.
Numerical Green’s function method:Application to quantifying ground motion variations of M7 earthquakes 
下载免费PDF全文
下载免费PDF全文 The concept of "numerical Green’s functions" (NGF or Green’s function database) is developed. The basic idea is: a large seismic fault is divided into subfaults of appropriate size, for which synthetic Green’s functions at the surface (NGF) are calculated and stored. Consequently, ground motions from arbitrary kinematic sources can be simulated, rapidly, for the whole fault or parts of it by superposition. The target fault is a simplified, vertical model of the Newport-Inglewood fault in the Los Angeles basin. This approach and its functionality are illustrated by investigating the variations of ground motions (e.g. peak ground velocity and synthetic seismograms) due to the source complexity. The source complexities are considered with two respects: hypocenter location and slip history. The results show a complex behavior, with dependence of absolute peak ground velocity and their variation on source process directionality, hypocenter location, local structure, and static slip asperity location. We concluded that combining effect due to 3-D structure and finite-source is necessary to quan- tify ground motion characteristics and their variations. Our results will facilitate the earthquake hazard assessment projects. 相似文献
2.
Three-component seismic and geoelectrical in-mine surveys were carried out in Lyukobanya colliery near Miskolc, Hungary to determine the in situ petrophysical parameter distributions and to detect inhomogeneities in the coal seam. The seismic measurements comprise an underground vertical seismic profile, using body waves, and an in-seam seismic amplitude-depth distribution and transmission survey, using channel waves. The geoelectrical measurements are based on the drift- and seam-sounding method. Interval traveltime-, amplitude-, multiple-filter- and polarization analysis methods are applied to the seismic data. They lead to a five-layer model for the strata including the coal seam. The coal seam and two underlying beds act as a seismic waveguide. The layer sequence supports the propagation of both normal and leaky mode channel waves of the Love- and Rayleigh type. A calculation of the total reflected energy for each interface using Knott's energy coefficients shows that the velocity ranges of high reflection energy and of normal and leaky mode wavegroups coincide. The excitation of wavegroups strongly depends on the seismic source. A simultaneous inversion of a geoelectrical drift- and seam-sounding survey prevents misinterpretations of the seismic data by clearly identifying the low-velocity coal seam as a high-resistivity bed. Calculations of dispersion and sounding curves improve the resolution of the slowness and resistivity in each layer. Both diminished amplitudes and distortions in the polarization of transmission seismo-grams and decreasing resistivities in a geoelectrical pseudosection of the coal seam are related to an inhomogeneity. A calculation of synthetic seismograms for Love and Rayleigh channel waves with the finite-difference and the Alekseev-Mikhailenko method agrees well with the field data for the main features, i.e., particular arrivals in the wave train, wavegroups, velocities and symmetries or asymmetries. This in-mine experiment demonstrates that the simultaneous acquisition, processing and interpretation of seismic and geoelectrical data improve the lithological interpretation of petrophysical parameter distributions. Coal seam inhomogeneities can also be detected more reliably by the two independent surveys than by one alone. 相似文献
3.
Field data from two-component in-seam seismic measurements are used to study roadway modes and their interaction with reflected seam waves. Using the multiple-filter technique to investigate the dispersion behaviour of the different waves, it can be shown that the roadway modes disperse very similarly to the related transmitted seam waves. However, because of the free surface of the coal face, the dispersion curves of the roadway modes show a velocity reduction and a slight shift to lower frequencies compared to those of the related transmitted seam waves. Polarization analysis using hodograms, rectilinearity and polarization angle confirms these results. The parameters found by polarization analysis can be used to design polarization filters which help to separate roadway modes and reflected events in the traveltime range of superposition in the presented field case. 相似文献
4.
In-seam seismic surveys with channel waves have been widely used in the United Kingdom and elsewhere to map coal-seams and to detect anomalous features such as dirt bands, seam thinning and thickening, and particularly in-seam faulting. Although the presence of cleat-induced anisotropy has been recognized in the past, almost all previous analyses have assumed homogeneous isotropic or transversely isotropic coal-seams. Channel waves, however, exhibit properties which cannot be fully explained without introducing anisotropy into the coal-seam. In particular, Love-type channel waves are observed for recording geometries where, in a homogeneous isotropic or transversely isotropic structure, the source would not be expected to excite transverse motion. Similarly, modes of channel-wave propagation display the coupled three-component motion of generalized modes in anisotropic substrates, which would not be expected for Rayleigh and Love wave motion in isotropy or in transversely isotropic media with azimuthal isotropy. We model the observed in-seam seismic channel waves with synthetic seismograms to gain an understanding of the effects of cleat-induced anisotropy on the behaviour of channel waves. The results show a reasonable good match with the observations in traveltime, relative amplitudes, dispersion characteristics and particle motions. We demonstrate that anisotropy in the surrounding country rocks contributes significantly to the coupling of channel wave particle motion, although its effect is not as strong as the anisotropy in the coal-seam. We conclude that the effects of cleat- and stress-induced anisotropy are observed and can be modelled with synthetic seismograms, and that anisotropy must be taken into account for the detailed interpretation of channel waves. 相似文献
5.
Synthetic seismograms for finite sources in spherically symmetric Earth using normal-mode summation 下载免费PDF全文
下载免费PDF全文 Normal-mode summation is the most rapidly used method in calculating synthetic seismograms. However, normal-mode summation is mostly applied to point sources. For earthquakes triggered by faults extending for as long as several 100 km, the seismic waves are usually simulated by point source summation. In this paper, we attempt to follow a different route, i.e., directly calculate the excitation of each mode, and use normal-mode summation to obtain the seismogram. Furthermore, we assume the finite source to be a ‘‘line source' and numerically calculate the transverse component of synthetic seismograms for vertical strike-slip faults. Finally, we analyze the features in the Love waves excited by finite faults. 相似文献
6.
Determining the focal mechanism of earthquakes helps us to better define faults and understand the stress regime. This technique can be helpful in the oil and gas industry where it can be applied to microseismic events. The objective of this paper is to find double couple focal mechanisms, excluding scalar seismic moments, and the depths of small earthquakes using data from relatively few local stations. This objective is met by generating three‐component synthetic seismograms to match the observed normalized velocity seismograms. We first calculate Green's functions given an initial estimate of the earthquake's hypocentre, the locations of the seismic recording stations and a 1D velocity model of the region for a series of depths. Then, we calculate the moment tensor for different combinations of strikes, dips and rakes for each depth. These moment tensors are combined with the Green's functions and then convolved with a source time function to produce synthetic seismograms. We use a grid search to find the synthetic seismogram with the largest objective function that best fits all three components of the observed velocity seismogram. These parameters define the focal mechanism solution of an earthquake. We tested the method using three earthquakes in Southern California with moment magnitudes of 5.0, 5.1 and 4.4 using the frequency range 0.1–2.0 Hz. The source mechanisms of the events were determined independently using data from a multitude of stations. Our results obtained, from as few as three stations, generally match those obtained by the Southern California Earthquake Data Center. The main advantage of this method is that we use relatively high‐frequency full‐waveforms, including those from short‐period instruments, which makes it possible to find the focal mechanism and depth of earthquakes using as few as three stations when the velocity structure is known. 相似文献
7.
High-speed train seismology has come into being recently. This new kind of seismology uses a high-speed train as a repeatable moving seismic source. Therefore, Green's function for a moving source is needed to make theoretical studies of the high-speed train seismology. Green's function for three-dimensional elastic wave equation with a moving point source on the free surface is derived. It involves a line integral of the Green's function for a fixed point source with different positions and corresponding time delays. We give a rigorous mathematical proof of this Green's function. According to the principle of linear superposition, we have also obtained the Green's function for a group of moving sources which can be regarded as a model of a traveling high-speed train. Based on a temporal convolution, an analytical formula for other moving sources is also given. In terms of a moving Gaussian source, we deal with the issue of numerical calculations of the analytical formula. Applications to modelling of a traveling high-speed train are presented. We have considered both the land case and the bridge case for a traveling high-speed train. The theoretical seismograms show different waveform features for these two cases. 相似文献
8.
9.
Xiao-Ping Li 《Geophysical Prospecting》1997,45(6):945-961
A deterministic pure phase shift filter (PPSF) is developed to extract the fundamental mode from multimode surface Love waves. Because of different phase velocities of modes and hence different phase traveltimes for a fixed travel distance between source and receiver, the deterministic PPSF can be computed, provided that the dispersion relation of the medium is estimated from the existing transmission data. The process consists of (a) applying the deterministic PPSF to the multimode wave (this step of the filtering process results in a time series in which the amplitudes of the fundamental mode appear at acausal times and the amplitudes of higher modes appear at causal times); (b) setting amplitude values equal to zero for positive times; (c) applying the inverse PPSF to the filtered signal. By using such a deterministic PPSF process, the higher modes almost disappear. The method is applied to synthetic multimode data computed by the normal-mode summation method. 相似文献
10.
Dispersion analysis is an important part of in-seam seismic data processing, and the calculation accuracy of the dispersion curve directly influences pickup errors of channel wave travel time. To extract an accurate channel wave dispersion curve from in-seam seismic two-component signals, we proposed a time–frequency analysis method based on single-trace signal processing; in addition, we formulated a dispersion calculation equation, based on S-transform, with a freely adjusted filter window width. To unify the azimuth of seismic wave propagation received by a two-component geophone, the original in-seam seismic data undergoes coordinate rotation. The rotation angle can be calculated based on P-wave characteristics, with high energy in the wave propagation direction and weak energy in the vertical direction. With this angle acquisition, a two-component signal can be converted to horizontal and vertical directions. Because Love channel waves have a particle vibration track perpendicular to the wave propagation direction, the signal in the horizontal and vertical directions is mainly Love channel waves. More accurate dispersion characters of Love channel waves can be extracted after the coordinate rotation of two-component signals. 相似文献
11.
A local plane-wave approach of generalized diffraction tomography in heterogeneous backgrounds, equivalent to Kirchhoff summation techniques when applied in seismic reflection, is re-programmed to act as repeated synthetic aperture radar (SAR) imaging for seismic prestack depth migration. Spotlight-mode SAR imaging quickly provides good images of the electromagnetic reflectivity of the ground via fast Fourier transform (FFT)-based signal processing. By calculating only the Green's functions connecting the aircraft to the centre of the illuminated patch, scattering structures around that centre are also recovered. SAR technology requires us to examine seismic imaging from the local point of view, where the quantity and quality of the available information at each image point are what are important, regardless of the survey geometry. When adapted to seismics, a local image of arbitrary size and sampling is obtained by FFT of seismic energy maps in the scattering wavenumber domain around each node of a pre-calculated grid of Green's functions. These local images can be used to generate a classic prestack depth-migrated section by collecting only their centres. However, the local images also provide valuable information around the centre, as in SAR. They can therefore help to pre-analyse prestack depth migration efficiently, and to perform velocity analysis at a very low cost. The FFT-based signal-processing approach allows local, efficient and automatic control of anti-aliasing, noise and resolution, including optimized Jacobian weights. Repeated local imaging could also be used to speed up migration, with interpolation between local images associated with a coarse grid of Green's functions, as an alternative to interpolation of Green's functions. The local images may, however, show distortions due to the local plane-wave approximation, and the velocity variations across their frame. Such effects, which are not necessarily a problem in SAR, should be controlled and corrected to further enhance seismic imaging. Applications to realistic models and to real data show that, despite the distortion effects, the local images can yield similar information to prestack depth migration, including common-image-point gathers for velocity analyses and AVO/AVA effects, at a much lower cost when a small target is considered. 相似文献
12.
Synthetic seismograms generated by solving the seismic wave equation using numerical methods are being widely used in seismology. For fully three-dimensional seismic structure models, the generation of these synthetic seismograms may require large amount of computing resources. Conventional high-performance computer clusters may not provide a cost-effective solution to this type of applications. The newly emerging cloud-computing platform provides an alternative solution. In this paper, we describe our implementation of a synthetic seismogram generator based on the reciprocity principle using the Windows Azure cloud application framework. Our preliminary experiment shows that our cloud-based synthetic seismogram generator provides a cost-effective and numerically efficient approach for computing synthetic seismograms based on the reciprocity principle. 相似文献
13.
We present the theory and numerical results for interferometrically interpolating 2D and 3D marine surface seismic profiles data. For the interpolation of seismic data we use the combination of a recorded Green's function and a model‐based Green's function for a water‐layer model. Synthetic (2D and 3D) and field (2D) results show that the seismic data with sparse receiver intervals can be accurately interpolated to smaller intervals using multiples in the data. An up‐ and downgoing separation of both recorded and model‐based Green's functions can help in minimizing artefacts in a virtual shot gather. If the up‐ and downgoing separation is not possible, noticeable artefacts will be generated in the virtual shot gather. As a partial remedy we iteratively use a non‐stationary 1D multi‐channel matching filter with the interpolated data. Results suggest that a sparse marine seismic survey can yield more information about reflectors if traces are interpolated by interferometry. Comparing our results to those of f‐k interpolation shows that the synthetic example gives comparable results while the field example shows better interpolation quality for the interferometric method. 相似文献
14.
XIAOHONG XU 《Geophysical Prospecting》1990,38(2):139-168
A multichannel lattice filter structure is utilized to represent seismic waves propagating in adjacent layers in an elastic medium. Using this model, an explicit time-domain solution for arbitrary source and receiver locations is obtained as an ARMA (AutoRegressive and Moving-Average) process. The lattice and ARMA structures have given rise to an effective algorithm for the calculation of offset/downhole synthetic seismograms. A large range of recently developed offset/downhole seismic survey geometries, such as the ‘Yo-Yo’ arrangement, can thus be simulated. In addition, the explicit solutions for upgoing and downgoing waves provide new insight into the properties of general downhole seismic signals, including wave-mode conversion effects and multiple reflections. Furthermore, offset/downhole seismograms generated by a line source (i.e. 2D point source) can also be constructed by superposition of plane waves with different incidence angles. Synthetic seismograms generated using a different source-receiver arrangement indicate that the properties especially associated with offset/downhole seismic signals can be predicted by this modelling method. These properties include arrival times, amplitude attenuation and wave-mode conversion effects. Finally, utilizing this numerical modelling method to a real downhole survey with Yo-Yo geometry may lead to a proper data acquisition and processing procedure, and improves the interpretation confidence of the field section. 相似文献
15.
高分辨率的面波频散谱成像是浅层地震勘探领域基于频散性质反演横波速度结构中的一个关键步骤.在天然地震探测领域,仅利用两个台站记录的线性信号比较法(LSC),被广泛用来计算面波的频散谱,并用于大尺度的面波层析成像.然而互相关的成像方式会造成频散谱在低频端较低的分辨率.非线性信号比较法(NLSC)利用指数函数克服了这个问题,同时极大地提高了频散谱的成像分辨率.然而,在研究中我们发现,仅利用两个台站的地震记录,并不能将面波的频散特性完整地考虑在内,导致LSC和NLSC方法对高阶模式的成像存在较大的误差.由于主动源面波勘探多道采集的方式,基于信号比较理论的多道信号比较法(MSC)充分利用多道地震信号,可以获得准确的多模式成像,然而该方法需要计算任意两道的频散谱并叠加,存在冗余的计算,导致计算效率较低.因此,本研究对MSC方法进行了相应的改进,通过追踪炮集记录上的面波波组提高了原方法的计算效率,同时,利用理论频散曲线进行叠加分析,验证了改进的MSC方法的正确性和有效性.通过与相移法、LSC和NLSC方法的对比分析,展示了MSC方法是一种准确的、高分辨率的面波多模式频散谱成像方法.实际地震资料的应用,揭示了MSC方法在浅层地震勘探中用于提取面波多模式频散信息的巨大潜力. 相似文献
16.
１９８０年和１９９２年分别在昆明和北京白家疃台架设了ＩＤＡ观测仪器，１９８５年开始又陆续在北京，兰州，海拉尔等地布设和完善了ＣＤＳＮ数字化地震观测台网，为开展长周期面波研究提供了方便。本文结合国内外研究现状，论述了超长周期地震观测在地球自由振荡，震源物理面波理论地震图合成及非均匀地球介质结构研究中的应用。同时结合１９８９年５月６日发生在秘鲁的Ｍｓ７．０地震，介绍了利用简正振型资料开发地震矩张量反演 相似文献
17.
True amplitude migration is one of the most important procedures of seismic data processing. As a rule it is based on the decomposition of the velocity model of the medium into a known macrovelocity component and its sharp local perturbations to be determined. Under this decomposition the wavefield can be considered as the superposition of an incident and reflected/scattered waves. The single scattering approximation introduces the linear integral operator that connects the sharp local perturbations of the macrovelocity model with the multishot/multioffset data formed from reflected/scattered waves. We develop the pseudoinverse of this operator using the Gaussian beam based decomposition of acoustic Green's functions. The computation of this pseudoinverse operator is done pointwise by shooting Gaussian beams from the target area towards the acquisition system. The numerical implementation of the pseudoinverse operator was applied to the synthetic data Sigsbee2A. The results obtained demonstrate the high quality of the true amplitude images computed both in the smooth part of the model and under the salt body. 相似文献
18.
Consistency of the spatial autocorrelation method with seismic interferometry and its consequence 总被引：1，自引：0，他引：1
We have cross‐checked the conventional theory of the spatial autocorrelation method and the consequence of seismic interferometry: the retrieval of the elastodynamic Green's function. Their mutual consistency is almost complete. The basic formulas of the conventional spatial autocorrelation theory can be derived by an alternative approach based on the retrieval of the elastodynamic Green's function. The only discrepancy is found with the average of the complex coherence function over azimuth in a wavefield dependent on azimuth. It is hypothesized, in discussion, that this discrepancy is due to the way of representing the wavefield in the background theory of seismic interferometry that can produce only wavefields moderately dependent on azimuth and that the mentioned consequence of seismic interferometry can also only make sense in a wavefield moderately dependent on azimuth. Our field experiment with a wavefield dependent on azimuth showed that the consequence of seismic interferometry in the logical framework of the conventional spatial autocorrelation theory is appropriate under such degrees of approximation as the measure proposed in this study, i.e., the deviation of the total dispersion curves is between about 10 and 16 per cent at the maximum from those averaged over azimuth. The acceptance of the retrieval of Green's function gives a proper physical meaning to the complex coherence function: the real part of the elastodynamic Green's function normalized by its zero‐offset version. This makes it possible to take a deterministic approach rather than the statistical one on which the conventional spatial autocorrelation method is based and gives fruitful new aspects and perspectives. For example, the formula for the multi‐mode case is given and the possibility of exploration of two or three dimensional velocity structures is suggested. 相似文献
19.
Bao-Li Wang 《应用地球物理》2018,15(1):118-124
Accurately detecting the arrival time of a channel wave in a coal seam is very important for in-seam seismic data processing. The arrival time greatly affects the accuracy of the channel wave inversion and the computed tomography (CT) result. However, because the signal-to-noise ratio of in-seam seismic data is reduced by the long wavelength and strong frequency dispersion, accurately timing the arrival of channel waves is extremely difficult. For this purpose, we propose a method that automatically picks up the arrival time of channel waves based on multi-channel constraints. We first estimate the Jaccard similarity coefficient of two ray paths, then apply it as a weight coefficient for stacking the multichannel dispersion spectra. The reasonableness and effectiveness of the proposed method is verified in an actual data application. Most importantly, the method increases the degree of automation and the pickup precision of the channel-wave arrival time. 相似文献
20.
An in-seam fan shoot was conducted in 1981 over a 300 m × 500 m panel of the Schwalbach seam at Ensdorf by a team from Prakla-Seismos AG of Hannover under contract to Saarbergwerke AG, Saarbrücken. The object was to study SH-mode propagation in the coal seam waveguide. The high quality dataset retrieved provides a general and valuable test bed with which to compare the performance of in-seam seismic velocity analysers. Five different dispersion analysers are demonstrated using the Schwalbach data. They are all based on the a priori assumption of coal seam homogeneity and isotropy. Space or time windows limit the resolution of the Fourier moving-window analyser, the migration based phase-velocity analyser, and the double Fourier transformer. The other two analysers, the maximum entropy moving-window analyser and the phase-moveout analyser, achieve noise-limited super-resolution by predicting the probable behaviour of the wavefield outside the window. The coal seam's characteristics conform to those predicted for a simplified model based on proposals by Elsen, Rüter and Schott of Westfälische Berggewerkschaftskasse, Bochum. The slight discrepancy between theoretical and actual dispersion characteristics could be reduced by increasing the model's complexity. However, there would be no material gain without testing the validity of the signal processing assumptions of seam isotropy and homogeneity. 相似文献