首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Synthetic seismograms are usually computed for reflections from vertical incidence of P waves for a horizontally layered medium. In actual practice the angle of incidence departs from the vertical, as receivers are usually located at some distance from the source. At angles other than the vertical, the conversion of P- to S-wave energy and changes in the reflection coefficient affect the shape of the synthetic seismograms. The effect of non-vertical incidence on synthetic seismograms is examined in this paper. Seismograms at non-vertical incidence have been computed using the plane-wave approach of Haskell (1953) for a layered medium. The use of plane waves is an approximation to the actual case of spherical wavefronts from a surface source. Using plane-wave theory, the expected wave forms as a function of angle of incidence were computed numerically for several simple models. The results indicate that the synthetic seismograms do not change significantly for angles of incidence between o and 25 degrees. For larger angles the changes in the wave forms may be severe. The effect is more pronounced for high-velocity layers than for low-velocity layers.  相似文献   

2.
The Kunetz-Claerbout equation for the acoustic transmission problem in a layered medium in its original form establishes the relation between the transmission and the reflec tion response for P-waves in an horizontally layered medium and with vertical incidence. It states that the reflection seismogram due to an impulsive source at the surface, is one side of the autocorrelation of the seismogram due to an impulsive source at depth and a surface receiver. By adapting Claerbout's formulation to the transmission of SH-waves, the Kunetz-Claerbout equation also holds for reflection and transmission coefficients dependent on the incident angle. Thus, earthquake geometry SH-transmission seismograms can be used to caculate corresponding pseudoreflection seismograms which, in turn, can be inverted for the impedance structure using the Levinson algorithm. If the average incidence angle is known, a geometrical correction on the resulting impedance model can improve the resolution of layer thicknesses. In contrast to the inversion of reflection seismograms, the Levinson algorithm is shown to yield stable results for the inversion of transmission seismograms even in the presence of additive noise. This noise stabilization is inherent to the Kunetz-Claerbout equation. Results of inverted SH-wave microearthquake seismograms from the Swabian Jura, SW Germany, seismic zone obtained at recording site Hausen im Tal have been compared with sonic-log data from nearby exploration drilling at Trochtelfingen. The agreement of the main structural elements is fair to a depth of several hundred metres.  相似文献   

3.
We consider the problem of acoustic pulse propagation through a layered medium with a reflector at one end. The fluctuations in the medium properties are assumed to be strong, i.e. of finite amplitude, rapid in comparison to the typical wavelength and to have statistical structure. The depth of the reflector is assumed to be large in comparison to the wavelength. In this regime, simple formulae for the reflected pulse and its arrival time at the surface are obtained. The amplitude of the pulse is broadened and attenuated as a result of multiple scattering: the fine-layered structure of the medium can be characterized by a single constant which appears in the formula for the limiting waveform and which measures the size of the fluctuations in the medium. Within the theory, the commonly observed discrepancy between the integrated sonic traveltime and the seismic traveltime can be studied and understood. The theory is a natural extension of the long-wavelength effective medium theory of Backus. The analysis is rigorous and based on the invariant embedding technique.  相似文献   

4.
A hybrid seismic modelling technique has been developed to investigate complex geological phenomena. Those parts of a geological structure which are too complicated to be treated theoretically are studied by two-dimensional physical models; other sections of the structure which can be treated theoretically, i.e., inhomogeneities in the vertical direction, are modelled by computer methods. A feedback process is used to combine the results of both physical and computer modelling. Horizontally layered coal-seam models are presented to test the hybrid modelling technique for normal incidence. A comparison of the hybrid seismograms with pure synthetic seismograms shows an acceptable conformity for normal incidence. A hybrid zero-offset section is shown to investigate a complex geological structure in the Ruhr coalfield in Germany.  相似文献   

5.
The propagation of acoustic waves through a periodic layered medium is analyzed by an eigenvalue decomposition of the propagator matrix. This reveals how the velocity and attenuation of the layered medium vary as function of the periodic structure, material parameters and frequency. There are two important parameters which control the wave propagation in the periodic medium: the reflection coefficient and the ratio between one‐way traveltimes of the two parts of the cyclic layered medium. For low frequencies (large values of wavelength to layer thickness), the layered structure behaves as an effective medium, then there is a transition zone, and for higher frequencies (small values of wavelength to layer thickness) the medium is described by the time‐average velocity. In this paper we mostly concentrate on the transition zone between an effective medium and time‐average medium regimes. The width of the transition zone increases with larger values of the reflection coefficient. The transition zone corresponds to a blocking regime for which the transmission response of the layered structure is close to zero. For even higher frequencies, the time‐average medium is replaced by a new transition zone, and then again a time‐average medium. This pattern is periodically repeated with higher frequencies. For small values of the reflection coefficient, the transition between effective medium and time‐average medium occurs around a value of wavelength to layer thickness equal to 4.  相似文献   

6.
Griddispersioningeneratingfinite-differencessyntheticseismogramsAbdolrahimJavaherian(InstituteofGeophysics,theUniversityofTeh...  相似文献   

7.
An indirect boundary integral method to obtain the three-dimensional response of an infinitely long, layered, viscoelastic valley of arbitrary cross-section embedded in a layered viscoelastic half-space is presented. The valley is excited by homogeneous plane waves impinging at an oblique angle with respect to the axis of the valley. The method and associated computer programs are tested by comparison with available results in the limiting two-dimensional case of incidence normal to the axis of the valley. Additional comparisons with previous three-dimensional results obtained by a hybrid finite element-boundary integral method for cylindrical valleys subjected to obliquely incident waves show large differences. However, the results obtained here for an infinitely long valley appear to be in some agreement with earlier results for an elongated prolate semi-ellipsoidal valley and with results obtained by a discrete wavenumber boundary element approach. An extensive bibliography on the dynamic response of valleys is also presented.  相似文献   

8.
川南、滇北地震台网区域壳幔介质分层结构的初步研究   总被引:1,自引:1,他引:1  
本文利用川南、滇北21个地震台短周期地震仪的深源远地震波形作为观测资料,通过调整地壳和上地幔介质模型并计算相应的综合地震图的方式,直至获得与观测波形可相比拟的综合地震图。据此,推断该地震台网区域的地壳中下部由一组高速和低速相间的薄层所组成,且大部分台站下的地壳厚度达到或超过54km。此外,文中还研究了该区域壳幔介质的吸收特点。  相似文献   

9.
利用合成地震图研究东北地区壳幔介质结构   总被引:1,自引:0,他引:1       下载免费PDF全文
林美  龚绍京 《地震学报》1991,13(4):480-488
本文采用计算合成地震图与分布在东北地区的11个地震台记录的深源远震波形对比分析的方法,研究该区壳幔结构.结果表明,东北地区地下100km 以上壳幔介质为分层结构,并存在有低速层.莫霍界面较清晰.地壳厚度在31.8——35.8km 之间.   相似文献   

10.
本文将计算综合地震图所用的反射法,推广于层状介质中的剪切位错内源情况。  相似文献   

11.
地震泼由多种频率成份组成,在其传播过程中。由于地球介质的滤波作用,随着传播距离的增大。高频成份逐渐被衰减。周期越大的成份传播越远。因此,在地震观测上,不同频带的地震仪用于观测不同震中距范围的地震,使仪器较好地响应,就我国台网而言,短周期仪用于地方震、近震的监测,中长周期仪和长周期仪用于远震,极远震的监测,其中长周期仪偏重于记录极远震。但作者在实际分析工作中注意到,短周期仪的记录应用于远震,极远震的分析,显示出独特的作用,充分挖掘和利用短仪资料,将会更加丰富远震、极远震的震相资料。本文以高台地震台的资料为依据,从震中距和震级角度,与中长仪和长仪作比较,对短仪记录远震、极远震的震相及特征作了统计分析与初步探讨,并对新疆地区“影区”地震S波的记录、日本地区地震_PPcP等震相的出现提出了现象的存在,给出了相应的观测结果,以待步一步研究。  相似文献   

12.
The indirect boundary element method is used to study the 3D dynamic response of an infinitely long alluvial valley embedded in a saturated layered half-space for obliquely incident SV waves. A wave-number transform is first applied along the valley’s axis to reduce a 3D problem to a 2D plane strain problem. The problem is then solved in the section perpendicular to the axis of the valley. Finally, the 3D dynamic responses of the valley are obtained by an inverse wave-number transform. The validity of the method is confirmed by comparison with relevant results. The differences between the responses around the valley embedded in dry and in saturated poroelastic medium are studied, and the effects of drainage conditions, porosity, soil layer stiffness, and soil layer thickness on the dynamic response are discussed in detail resulting in some conclusions.  相似文献   

13.
Numerical examples of high-frequency synthetic seismograms of body waves in a 2-D layered medium with complex interfaces (faults, wedges, curvilinear, corrugated) are presented. The wave field modeling algorithm combines the possibilities of the ray method and the edge wave superposition method. This approach preserves all advantages of the ray method and eliminates restrictions related to diffraction by boundary edges and to caustic effects in singular regions. The method does not require two-point ray tracing (source-to-receiver), and the position of the source, as well as the type of source, and the position of receivers can be chosen arbitrarily. The memory and the time required for synthetic seismogram computation are similar to ray synthetic seismograms. The computation of the volume of the medium (the Fresnel volume or Fresnel zones), which gives the essential contribution to the wave field, is included in the modeling program package. In the case of complicated irregular interface (or a layered medium with a regular ray field at the last interface), the method displays a high accuracy of wave field computation. Otherwise, the method can be considered a modification of the ray method with regularization by the superposition of edge waves.  相似文献   

14.
We consider nine earthquakes in the Koyna-Warna reservoir region on the western side of the Peninsular India. The deviatoric moment tensors of these earthquakes have been evaluated by minimizing the least-squares misfit between observed and synthetic seismograms. We use broadband seismograms of observatories at KARD and PUNE which are at distances of nearly 50 and 150 km, respectively, from the epicenters. Both surface wave inversion and the difference between the arrival times of SH and SV show the presence of an anisotropic crust. However, we have obtained an equivalent isotropic structure by improving the published crustal structures of this area through inversion of surface wave group velocity data. The deviatoric moment tensors of the earthquakes are decomposed into two components: double-couple and compensated linear vector dipoles (CLVD). The double-couple components of all the nine earthquakes show normal faulting with minor strike slip; the T axis is consistently subhorizontal with an average azimuth of 260.6° and the P axis is nearly vertical. The fault planes of six events give average strike direction and dip, respectively as 194.0° and 51.8° and are associated with the main fault of the area. The other three events lie in the southern part of this area and have strike direction between SSE and SE which is parallel to the tectonic features in this part. The CLVD component is generally within 20 percent of the total moment tensor. Recent studies show that anisotropy can produce source mechanism with CLVD up to 30 percent and can also cause high pore fluid pressure leading to fault instability more rapidly compared to conventional mechanism in an isotropic medium. It appears that the anisotropic crust, noted in the present work, is generating the CLVD component and also gives the proper environment to trigger earthquakes by reservoirs through pore fluid pressure.  相似文献   

15.
地震动入射角度对地下结构地震响应的影响   总被引:1,自引:0,他引:1  
考虑土-结构接触面效应和场地初始静应力影响,基于大型商用有限元软件ANSYS和粘弹性人工边界条件,采用动力松弛法的分析思路,建立了一种地震动斜入射条件下地下结构的接触非线性动力反应分析模型和方法,并讨论了地震动入射角度对地下结构动力反应的影响。结果表明:地震波斜入射使得结构的整体反应发生明显变化;随着入射角度的增加,节点的水平向应力反应明显增大,竖向应力峰值较小,增大程度也相对较小;节点的位移峰值随输入加速度峰值的增大也有一定的变化。因此,在分析近源地震作用下的地下结构动力响应时,需要考虑地震动的非一致输入问题。  相似文献   

16.
A new time-domain method is introduced for the calculation of theoretical seismograms which include frequency dependent effects like absorption. To incorporate these effects the reflection and transmission coefficients become convolutionary operators. The method is based on the communication theory approach and is applicable to non-normal incidence plane waves in flat layered elastic media. Wave propagation is simulated by tracking the wave amplitudes through a storage vector inside the computer memory representing a Goupillaud earth model discretized by equal vertical transit times. Arbitrary numbers of sources and receivers can be placed at arbitrary depth positions, while the computational effort is independent of that number. Therefore, the computation of a whole plane-wave vertical seismic profile is possible with no extra effort compared to the computation of the surface seismogram. The new method can be used as an aid to the interpretation of plane-wave decomposed reflection data where the whole synthetic vertical seismic profile readily gives the interpreter the correct depth position of reflection events.  相似文献   

17.
Prediction of acoustic impedances below bottom-of-hole developed in the paper is based on VSP data processing. The pulse form and reflections are deduced from records on the vertical array by the method of subtraction. A deterministic prediction error and a wave-shaping Wiener filter are used to transform the pulse train and reflections into a short symmetrical pulse. Additional broadening of the pulse spectrum is achieved by stacking of seismograms from shots with various dominant frequencies after zero-phase transformation of downgoing waves. The inversion of seismograms is made after ideal spiking by means of a subtraction procedure. Acoustic impedances deduced are closely related to lithology of rocks below bottom-of-hole.  相似文献   

18.
In this paper, structure models of the crust and upper mantle beneath each station have been obtained by way of fitting synthetic seismograms with P waveforms of deep focus teleseismic records from the 11 stations in the Northeastern Region of China. We have studied the structure in the region based on those models. Our results show that the medium of the crust and upper mantle is a layered structure with alternate high and low velocity layers within about 100 km under the region. The crustal thickness is about 31.8–35.8 km. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 471–479, 1991. This work is one part of the project funded by the State Seismological Bureau of China.  相似文献   

19.
While velocity contrasts are responsible for most of the events recorded in our data, the long wavelength behavior of the velocity model is responsible for the geometrical shape of these events. For isotropic acoustic materials, the wave dependency on the long (wave propagation) and short (scattering) wavelength velocity components is stationary with the propagation angle. On the other hand, in representing a transversely isotropic with a vertical symmetry axis medium with the normal moveout velocity, the anellepticity parameter η, the vertical scaling parameter δ, and the sensitivity of waves vary with the polar angle for both the long and short wavelength features of the anisotropic dimensionless medium parameters (δ and η). For horizontal reflectors at reasonable depths, the long wavelength features of the η model is reasonably constrained by the long offsets, whereas the short wavelength features produce very week reflections at even reasonable offsets. Thus, for surface acquired seismic data, we could mainly invert for smooth η responsible for the geometrical shape of reflections. On the other hand, while the δ long wavelength components mildly affects the recorded data, its short wavelength variations can produce reflections at even zero offset, with a behavior pattern synonymous to density. The lack of the long wavelength δ information will mildly effect focusing but will cause misplacement of events in depth. With low enough frequencies (very low), we may be able to recover the long wavelength δ using full waveform inversion. However, unlike velocity, the frequencies needed for that should be ultra‐low to produce long‐wavelength scattering‐based model information as δ perturbations do not exert scattering at large offsets. For a combination given by the horizontal velocity, η, and ε, the diving wave influence of η is absorbed by the horizontal velocity, severely limiting the η influence on the data and full waveform inversion. As a result, with a good smooth η estimation, for example, from tomography, we can focus the full waveform inversion to invert for only the horizontal velocity and maybe ε as a parameter to fit the amplitude. This is possibly the most practical parametrization for inversion of surface seismic data in transversely isotropic with vertical symmetry axis media.  相似文献   

20.
2011年1月19日发生的安庆MS4.8地震,在较远范围引起了明显震感.通过仔细辨认此次地震不同频段的波形记录,发现在一些台站记录到了清晰的SmS震相.采用CAP方法反演了其震源机制解,并利用理论地震图与实际观测数据对比进一步确认了SmS震相的存在.结合已有实验结果,推断在人类有感频率(1Hz至几Hz)范围内,震中距70——200 km内SmS震相往往是振幅最大的震相,是引起远距离有感的主要原因.利用理论地震图研究了震源深度对SmS震相的影响. 结果表明,震源深度可改变SmS震相发育的临界震中距,震源越深对应的临界震中距越小.进一步通过对比SmS与S的振幅比,定性讨论了震源机制解中倾角、滑动角和地壳浅层衰减等因素对SmS震相的影响.结果表明,倾角和滑动角对其振幅比影响呈现比较复杂的关系,而SmS与S振幅比随地壳浅层衰减的增大而减小,说明浅层衰减对SmS震相影响较大.综合分析认为,除当地地壳速度结构的影响外,SmS震相的发育受到震源机制解的倾角、滑动角,震源深度以及地壳浅层衰减等多种因素的影响.   相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号