共查询到20条相似文献，搜索用时 257 毫秒
1.
常规的地震定位方法通常需要拾取地震记录的初至,当初至不明显或被较高水平的噪声淹没时精度较低.本文采用基于三维高斯射线束的偏移成像方法对震源进行定位,较好地解决了该问题.通过三维高斯射线束对台站记录进行偏移归位,并将各台站成像结果的交点作为地震能量释放的中心位置;当各台站成像结果不能交于一点时,采用三维空间高斯滤波方法可实现震源位置的自动获取.提出的变网格计算方案极大地减少了计算量,显著地提高了成像精度和计算效率.利用首都圈地震台网数据,对涿鹿、滦县以及房山三个地震事件进行试算,结果表明:基于变网格三维高斯束偏移成像的地震定位方法自动化程度很高,而且具有较好的抗噪能力,特别适合处理低信噪比资料的地震定位问题. 相似文献
2.
3.
R. L. Nowack 《Pure and Applied Geophysics》2003,160(3-4):487-507
— In this paper, an overview of the calculation of synthetic seismograms using the Gaussian beam method is presented accompanied by some representative applications and new extensions of the method. Since caustics are a frequent occurrence in seismic wave propagation, modifications to ray theory are often necessary. In the Gaussian beam method, a summation of paraxial Gaussian beams is used to describe the propagation of high-frequency wave fields in smoothly varying inhomogeneous media. Since the beam components are always nonsingular, the method provides stable results over a range of beam parameters. The method has been shown, however, to perform better for some problems when different combinations of beam parameters are used. Nonetheless, with a better understanding of the method as well as new extensions, the summation of Gaussian beams will continue to be a useful tool for the modeling of high-frequency seismic waves in heterogeneous media. 相似文献
4.
Robert L. Nowack Wang-Ping Chen Ulrich Kruse Saptarshi Dasgupta 《Pure and Applied Geophysics》2007,164(10):1921-1936
We carry out a sequence of numerical tests to understand conditions under which rapid changes in crustal thickness can be
reliably imaged by teleseismic body waves. Using the finite-difference method over a 2-D grid, we compute synthetic seismograms
resulting from a planar P-wavefield incident below the grid. We then image the Moho using a migration scheme based on the Gaussian beam representation
of the wavefield. The use of Gaussian beams for the downward propagation of the wavefield is particularly advantageous in
certain geologically critical cases such as overthrusting of continental lithosphere, resulting in the juxtaposition of high-velocity
mantle material over crustal rocks. In contrast to ray-based methods, Gaussian beam migration requires no special treatment
to handle such heterogeneities. Our results suggest that with adequate station spacing and signal-to-noise ratios, offsets
of the Moho, on the order of 10 km in height, can be reliably imaged beneath thickened crust at depths of about 50 km. Furthermore,
even sharp corners and edges are faithfully imaged when precise values of seismic wave speeds are available. Our tests also
demonstrate that flexibility in choices of different types of seismic phases is important, because any single phase has trade-offs
in issues such as spatial resolution, array aperture, and amplitude of signals. 相似文献
5.
— The effects of interfaces and velocity gradients on wide-angle seismic attributes are investigated using synthetic seismograms. The seismic attributes considered include envelope amplitude, pulse instantaneous frequency, and arrival time of selected phases. For models with interfaces and homogeneous layers, head waves can propagate which have lower amplitudes, as well as frequency content, compared to the direct arrivals. For media with interfaces and velocity gradients, higher amplitude diving waves and interference waves can also occur. The Gaussian beam and reflectivity methods are used to compute synthetic seismograms for simple models with interfaces and gradients. From the results of these methods, seismic attributes are obtained and compared. It was found that both methods were able to simulate wide-angle seismic attributes for the simple models considered. The advantage of using the Gaussian beam method for seismic modeling and inversion is that it is fast and also asymptotically valid for laterally varying media. 相似文献
6.
本文通过格林函数反褶积方法，由台湾峡谷附近的记录资料预测峡谷区的强地面运动的时程曲线．峡谷区的理论格林函数应用２．５维ＳＨ混合方法求解．通过格林函数反褶积方法得到的峡谷区费丛１及费丛２台的位移、速度、加速度和实际资料对比，取得了满意结果．计算结果表明，峡谷底部的峰值加速度相对峡谷边缘为最小；在靠近震中的一侧，峡谷的加速度的最大振幅相对比另一侧大．还给出了其余３个台的预测结果，研究了它们的加速度傅里叶谱和反应谱． 相似文献
7.
本文提出了在计算机上实现分块三维地壳模型及利用加权最小二乘拟合生成平缓光滑的三维速度函数的方法,给出了适用于分块、块内速度连续变化的三维模型中Cauchy射线追踪的新算法,简介了基于上述方法反射线的基本理论所编制的合成三维理论地震图的程序包RSSGTD.给出的两个盆地状模型的算例表明,所使用的模型生成方法具有模拟复杂地壳结构的能力;与三维样条函数方法比较,最小二乘拟合方法能给出更加适合射线方法合成地震图计算的速度函数,并且内存小、计算速度快;所给出的Cauchy射线追踪算法能够适合块状模型中任何体波射线的追踪. 相似文献
8.
9.
The application of Maslov asymptotic theory in a general 3-D mixed subspace of 6-D complex phase space is proposed to obtain the integral superpositions of Gaussian packets and beams. The ray method and the superposition of plane waves (Maslov method of Chapman and Drumond [7]) are special limiting cases of the above mentioned approach. The same high-frequency asymptotic expansion formulae for seismic body waves were derived previously in [8] using the Gaussian beam method. 相似文献
10.
2.5-D modeling and inversion techniques are much closer to reality than the simple and traditional 2-D seismic wave modeling
and inversion. The sensitivity kernels required in full waveform seismic tomographic inversion are the Fréchet derivatives
of the displacement vector with respect to the independent anisotropic model parameters of the subsurface. They give the sensitivity
of the seismograms to changes in the model parameters. This paper applies two methods, called ‘the perturbation method’ and
‘the matrix method’, to derive the sensitivity kernels for 2.5-D seismic waveform inversion. We show that the two methods
yield the same explicit expressions for the Fréchet derivatives using a constant-block model parameterization, and are available
for both the line-source (2-D) and the point-source (2.5-D) cases. The method involves two Green’s function vectors and their
gradients, as well as the derivatives of the elastic modulus tensor with respect to the independent model parameters. The
two Green’s function vectors are the responses of the displacement vector to the two directed unit vectors located at the
source and geophone positions, respectively; they can be generally obtained by numerical methods. The gradients of the Green’s
function vectors may be approximated in the same manner as the differential computations in the forward modeling. The derivatives
of the elastic modulus tensor with respect to the independent model parameters can be obtained analytically, dependent on
the class of medium anisotropy. Explicit expressions are given for two special cases—isotropic and tilted transversely isotropic
(TTI) media. Numerical examples are given for the latter case, which involves five independent elastic moduli (or Thomsen
parameters) plus one angle defining the symmetry axis. 相似文献
11.
12.
-- Fault zones are thought to consist of regions with reduced seismic velocity. When sources are located in or close to these low-velocity zones, guided seismic head and trapped waves are generated which may be indicative of the structure of fault zones at depth. Observations above several fault zones suggest that they are common features of near fault radiation, yet their interpretation may be highly ambiguous. Analytical methods have been developed to calculate synthetic seismograms for sources in fault zones as well as at the material discontinuities. These solutions can be used for accurate modeling of wave propagation in plane-parallel layered fault zone structures. However, at present it is not clear how modest deviations from such simplified geometries affect the generation efficiency and observations of trapped wave motion. As more complicated models cannot be solved by analytical means, numerical methods must be employed. In this paper we discuss 3-D finite-difference calculations of waves in modestly irregular fault zone structures. We investigate the accuracy of the numerical solutions for sources at material interfaces and discuss some dominant effects of 3-D structures. We also show that simple mathematical operations on 2-D solutions generated with line sources allow accurate modeling of 3-D wave propagation produced by point sources. The discussed simulations indicate that structural discontinuities of the fault zone (e.g., fault offsets) larger than the fault zone width affect significantly the trapping efficiency, while vertical properly gradients, fault zone narrowing with depth, small-scale structures, and moderate geometrical variations do not. The results also show that sources located with appropriate orientations outside and below a shallow fault zone layer can produce considerable guided wave energy in the overlying fault zone layer. 相似文献
13.
2-D Versus 3-D Magnetotelluric Data Interpretation 总被引：6，自引：0，他引：6
In recent years, the number of publications dealing with the mathematical and physical 3-D aspects of the magnetotelluric method has increased drastically. However, field experiments on a grid are often impractical and surveys are frequently restricted to single or widely separated profiles. So, in many cases we find ourselves with the following question: is the applicability of the 2-D hypothesis valid to extract geoelectric and geological information from real 3-D environments? The aim of this paper is to explore a few instructive but general situations to understand the basics of a 2-D interpretation of 3-D magnetotelluric data and to determine which data subset (TE-mode or TM-mode) is best for obtaining the electrical conductivity distribution of the subsurface using 2-D techniques. A review of the mathematical and physical fundamentals of the electromagnetic fields generated by a simple 3-D structure allows us to prioritise the choice of modes in a 2-D interpretation of responses influenced by 3-D structures. This analysis is corroborated by numerical results from synthetic models and by real data acquired by other authors. One important result of this analysis is that the mode most unaffected by 3-D effects depends on the position of the 3-D structure with respect to the regional 2-D strike direction. When the 3-D body is normal to the regional strike, the TE-mode is affected mainly by galvanic effects, while the TM-mode is affected by galvanic and inductive effects. In this case, a 2-D interpretation of the TM-mode is prone to error. When the 3-D body is parallel to the regional 2-D strike the TE-mode is affected by galvanic and inductive effects and the TM-mode is affected mainly by galvanic effects, making it more suitable for 2-D interpretation. In general, a wise 2-D interpretation of 3-D magnetotelluric data can be a guide to a reasonable geological interpretation. 相似文献
14.
v--vAn estimation of local site effects in the Volvi basin as derived from observation and modeling is presented in this paper. The Volvi basin is located in the Mygdonian graben in northern Greece near the city of Thessaloniki. This test site has been studied and instrumented in the framework of the "EURO-SEISTEST" and "EURO-SEISMOD" projects funded by the European Union, aimed at improving knowledge of the influence of the local geology on the seismic response of a target area. In this context we calculate synthetic seismograms along a 2-D profile intersecting the graben, instrumented and accurately investigated with a geophysical survey and geotechnical tests. The seismic wavefield from the source to the target area has been computed with the modal summation method, while inside two of the investigated models representing the 2-D section, the wavefield has been numerically propagated with the finite-difference method. We compare the results of the two simulations, both in the time and frequency domain. We also compare the results with experimental data related to an event recorded by the Reftek network installed in the target area. This permits a better understanding of how the structural features of the 2-D models affect the seismic wavefield, especially in the frequency range between 2 Hz and 4 Hz, where one can observe differences between the simulations and the observations. This means that the general features of the models are able to reproduce observed amplification effects, apart from some discrepancies due to still unresolved structural features of the site. 相似文献
15.
丽江-攀枝花-者海地带二维地壳结构及其构造意义 总被引：26，自引：8，他引：26
本文利用攀西地区通过攀枝花的东西向剖面爆炸地震资料,进行了震相的重新识别和二维射线追踪与理论图计算.结果表明,沿剖面地表附近有4个低速区和若干高速带,它们与地质和构造有很好的对应关系.渡口附近的高速岩体一直延伸到了上地壳的底部,形成一个统一的地垒状构造,该高速体与攀枝花成矿岩体相关,并推断华坪及其以西地带也是找矿的有利地区.中地壳下部有一厚度约9km的低速层,它可能是壳内的韧性剪切带.低速层顶部深度为27.0-29.5km,与研究地区的居里面深度及天然地震震源深度的分布基本符合.剖面东段中地壳顶部还有一层很薄的低速层,反映了构造带两侧运动的不对称性.地壳厚度为53-56km,构造带中部的Moho界面没有明显的上隆. 相似文献
16.
文章提出了考虑剪切变形弹塑性刚度影响的多弹簧模型的空间梁柱单元,用于反复加载下钢构件的数值模拟。应用多轴应力状态下的塑性应力-应变关系理论,在单元模型中考虑了弹塑性区域剪切变形对单元的弹塑性刚度的影响,针对单元模型的塑性区长度和弹簧布置两个参数,文中给出了合理建议取值。数值模拟分析表明,所提出的单元模型能够很好地模拟钢构件的弹塑性性能。在此基础上,以多高层钢结构商业设计软件MTS为平台,进行三维钢框架结构弹塑性动力时程分析模块的开发。最后,文章对一纯钢框架结构足尺振动台试验进行数值模拟,模拟分析结果表明,本文所提出的多弹簧单元模型及开发的动力分析模块能够较好地模拟钢结构在地震作用下的弹塑性性能。 相似文献
17.
在给出真实模型和相应光滑背景模型的情况下,如何计算扰动模型(散射体)产生的散射波场是一个有实际意义的正演问题.在Gabor变换域描述散射体,且入射波场为短时宽带信号时,散射波场可以在频率域用高斯束或时间域用高斯波包描述.相对于波动方程方法,高斯束和高斯波包的计算效率更高;背景模型光滑时,高斯束和高斯波包方法的精度也接近波动方程方法.文中导出了声波假设下应用高斯束和高斯波包计算散射波的方法.测试分析了高斯波包的计算精度.给出了一般散射体的散射波模拟策略.同时针对一个理论模型完成了本文方法计算散射波的实验,实验结果表明高斯波包散射波计算方法是有效可行的. 相似文献
18.
19.
A global cross-section of the Earth parallel to the tectonic equator (TE) path, the great circle representing the equator of net lithosphere rotation, shows a difference in shear wave velocities between the western and eastern flanks of the three major oceanic rift basins. The low-velocity layer in the upper asthenosphere, at a depth range of 120 to 200 km, is assumed to represent the decoupling between the lithosphere and the underlying mantle. Along the TE-perturbed (TE-pert) path, a ubiquitous LVZ, about 1,000-km-wide and 100-km-thick, occurs in the asthenosphere. The existence of the TE-pert is a necessary prerequisite for the existence of a continuous global flow within the Earth. Ground-shaking scenarios were constructed using a scenario-based method for seismic hazard analysis (NDSHA), using realistic and duly validated synthetic time series, and generating a data bank of several thousands of seismograms that account for source, propagation, and site effects. Accordingly, with basic self-organized criticality concepts, NDSHA permits the integration of available information provided by the most updated seismological, geological, geophysical, and geotechnical databases for the site of interest, as well as advanced physical modeling techniques, to provide a reliable and robust background for the development of a design basis for cultural heritage and civil infrastructures. Estimates of seismic hazard obtained using the NDSHA and standard probabilistic approaches are compared for the Italian territory, and a case-study is discussed. In order to enable a reliable estimation of the ground motion response to an earthquake, three-dimensional velocity models have to be considered, resulting in a new, very efficient, analytical procedure for computing the broadband seismic wave-field in a 3-D anelastic Earth model. 相似文献
20.
L.A. Tabarovsky M.M. Goldman M.B. Rabinovich K.-M. Strack 《Journal of Applied Geophysics》1996,35(4):261-284
Understanding, using, or eliminating three-dimensional (3-D) effects in electromagnetic methods of geophysics are critical requirements. Numerous achievements in 3-D modeling sometimes give the impression that they are widely available today in geophysical practice. This is not necessarily true. Existing 3-D modeling packages prove that we know how to perform 3-D modeling. However, the computer resources and costs involved make the practical application of 3-D EM modeling in geophysical applications very limited.A practical compromise, or even alternative, is represented by 2.5-D modeling characterized by the use of a 3-D source in a 2-D medium. This combination allows one to mathematically describe the related boundary value problem as a sequence of independent two-dimensional problems. The typical technique leading to such a split formulation is Fourier analysis. That is why the individual terms of a split solution are often referred to as harmonics.Although each independent problem is two-dimensional, the algorithmic implementation of finite differences or integral equations for the higher harmonics has some specific features not present in the classical 2-D cases. In this paper, a hybrid scheme consisting of a combination of the finite difference technique with the integral equation approach for transient fields is described. Evaluation of algorithm accuracy is presented and a transient logging technique application is considered. The algorithm is fast and easily implemented on a personal computer 相似文献