共查询到20条相似文献,搜索用时 31 毫秒
1.
This paper is directed at modeling layered media. We extend the plane-wave normal-incidence state-space model developed by Mendel, Nahi and Chan in 1979, to the non-normal incidence case. To do this we introduce a shifting principle, a zero-offset wavefront, and zero-offset travel times for different layers. We also develop an algorithm for obtaining a synthetic line source reflection seismogram. In this algorithm non-normal incidence plane-wave seismograms are summed over a range of incident angles. The algorithm is based on a modified version of Sommerfield's (1896) theorem. Simulations of acoustic and elastic media are included which illustrate the applicability of our plane-wave and line source seismograms for both elastic and acoustic cases. 相似文献
2.
F. AMINZADEH 《Geophysical Prospecting》1989,37(8):893-906
Data from offshore Norway is used to study applications of elastic VSP modelling in detecting shear waves and observing the effects of successive mode conversion in field-recorded VSP data. The shear-wave velocities and densities from log data are used in conjunction with compressional wave velocities determined from surface seismic and log data in the VSP modelling. The time domain non-normal incidence elastic VSP modelling technique of Aminzadeh and Mendel is used as the modelling algorithm. Two surface seismograms are computed first. One is the vertical component and the other is the horizontal component for plane waves that have specified incident angles. A downward continuation method is then applied to generate seismograms at different depth points. The collection of these seismograms constitutes non-normal incidence VSPs. Both vertical and horizontal components of VSP data can be obtained by this procedure. In this paper non-normal incidence VSPs are generated for a 12.5° incident plane wave. The modelling results of layered earth systems of thin layers and thick layers are both compared with field data, and the effect of mode conversions in thin layers is observed. Several events in the field data can be explained by this elastic VSP modelling. Comparison of the model data and field data enabled a probable tube wave or out-of-plane event to be identified, the removal of which significantly improved the final VSP section. This study also shows how the VSP data helped the interpretation of the surface 3D data. 相似文献
3.
叠前同时反演是油气探测的一种有效工具.其理论基础是平面P波Zoeppritz方程计算的反射系数的近似,是入射角的函数.叠前同时反演可以利用三项或两项Fatti方程进行反演分析.本文针对实际油田的测井数据,利用反射率法模拟了仅包含P波一次反射记录,包含P波一次反射和P波层间多次波记录以及全波场地震记录,再利用叠前同时反演对合成地震记录进行反演研究.研究结果表明,在大偏移距处P波主要反射受到其它模式波的污染,从而影响了叠前同时反演结果的精度.对于薄互层介质当转换波影响严重时,使用小角度数据的两项AVO反演比使用大角度数据的三项AVO反演更合理可靠. 相似文献
4.
Margarita Luneva 《Pure and Applied Geophysics》1996,148(1-2):113-136
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. 相似文献
5.
A new approximate method to calculate the space-time acoustic wave motion generated by an impulsive point source in a horizontally layered configuration is presented. The configuration consists of a stack of fluid layers between two acoustic half-spaces where the source and the receiver are located in the upper half-space. A distorted-wave Born approximation is introduced; the important feature of the method is the assumption of a background medium with vertical varying root-mean-square acoustic wave speed. A closed-form expression for the scattered field in space and time as a function of the contrast parameters is deduced. The result agrees closely with rigorously calculated synthetic seismograms. In the inverse scheme the wave speed and mass density can be reconstructed within a single trace. Results of the inversion scheme applied to synthetic data are shown. 相似文献
6.
B.O. Ruud 《Studia Geophysica et Geodaetica》2006,50(3):479-498
The calculation of reflection and transmission coefficients of plane waves at a plane interface between two homogeneous anelastic
media may become ambiguous because it is not always obvious how to determine the sign of the vertical component of the slowness
vector of the scattered waves. For elastic media, the sign is determined by applying so-called radiation condition when the
slowness vector is complex-valued, but it has long been known that this approach does not work satisfactorily for anelastic
media. Other approaches have been suggested, e.g., by requiring that the reflection and transmission coefficients should vary
continuously with increasing incident angles, or by relating the sign to the direction of the energy flux. In the present
paper, it is shown that these approaches may give different results, and that the results can be inconsistent with the elastic
case even for weak attenuation. Instead, it is demonstrated that the ambiguity in the reflection coefficient can be resolved
by expressing the seismic response of a point source over an interface as a superposition of plane waves and their reflection
coefficients, and solving the resulting integral by the saddle point approximation. Although the saddle point itself (point
of stationary phase) does not provide new insight, the ambiguity is removed by considering the steepest descent path through
the point. Ray synthetic seismograms computed by this method compare well with synthetics computed by the reflectivity method,
which does not suffer from the above-mentioned ambiguity since the integration path is taken along the real axis.
This paper concentrates on the isotropic case, but it is discussed how the result may be extended to layered transversely
isotropic media. The suggested approach, derived for a point source and plane layers, does not directly apply to 2-D or 3-D
laterally inhomogeneous media, or to media of general anisotropy. A generalization of the result found is that the sign of
the vertical slowness components should be chosen according to the energy flux direction for subcritical incidence and according
to the radiation condition for supercritical incidence, even if this creates a discontinuity in the coefficients at the critical
incidence angle. Such a discontinuity is sometimes necessary to get results which are consistent with the elastic case. It
is discussed how the generalized result can be obtained by applying certain continuity criteria for the sub-and supercritical
angle intervals, but the validity of this approach for general models remains to be proved. 相似文献
7.
R.-G. FERBER 《Geophysical Prospecting》1988,36(8):857-877
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. 相似文献
8.
R. -G. Ferber 《Surveys in Geophysics》1989,10(2-4):133-154
The goal of seismic reflection surveys is the derivation of petrophysical subsurface parameters from surface measurements. Today's well established technique in data acquisition, as well as processing terms, is based on the acoustic approximation to the real world's wave propagation. In recent years a lot of work has been done to extend the technique to the elastic approximation. There was especially an important trend towards elastic inversion techniques operating on plane-wave seismograms, called simultaneous P-SV inversion (or short P-SV inversion) within this paper. Being still under investigation, some important aspects of P-SV inversion concerning data acquisition as well as pre-processing, should be pointed out. To fit the assumptions of P-SV inversion schemes, at least a two-dimensional picture of the reflected wavefield with vertical and in-line horizontal receivers has to be recorded. Moreover, the theoretical work done suggests that in addition to a survey with a compressional wave source, a second survey should be done using sources radiating vertically polarized shear waves, is needed. Finally, proper slant stacking must be performed to get plane-wave seismograms. The P/S separated plane-wave seismograms are then well prepared for feeding into the inversion algorithms. P/S separated planewave seismograms are then well prepared for feeding into the inversion algorithm.s In this paper, a tutorial overview of the data acquisition and pre-processing in accordance with the P-SV inversion philosophy is given and illustrated using synthetic seismograms. A judgement on the feasibility of the P-SV inversion philosophy must be left to ongoing research. 相似文献
9.
This paper presents results of testing an efficient ray generation scheme needed whenever ray synthetic seismograms are to be computed for layered models with more than 10‘ thick’layers. Our ray generation algorithm is based on the concept of kinematically equivalent waves (the kinematic analogs) having identical traveltimes along different ray-paths between the source and the receiver, both located on the surface of the model. These waves, existing in any medium composed of laterally homogeneous parallel layers, interfere at any location along the recording surface, thereby producing a composite wavelet whose amplitude and shape depend directly on the number of kinematic analogs (the multiplicity factor). Hence, explicit knowledge of the multiplicity factor is crucial for any analysis based on the amplitude and shape of individual wavelets, such as wavelet shaping, Q estimation, or linearized wavelet inversion. For unconverted waves, such as those discussed in this paper, the multiplicity factor can be computed analytically using formulae given in the Appendix; for converted waves, the multiplicity factor should be computed numerically, using the algorithm employed for the computation of the seismograms presented in a previous paper by one of the authors. 相似文献
10.
F. SCHERBAUM 《Geophysical Prospecting》1987,35(7):787-802
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. 相似文献
11.
A first-order one-way wave system has been created based on characteristic analysis of the acoustic wave system and optimization of the dispersion relation. We demonstrate that this system is equivalent to a third-order scalar partial-differential equation which, for a homogeneous medium, reduces to a form similar to the 45° paraxial wave equation. This system describes accurately waves propagating in a 2D heterogeneous medium at angles up to 75°. The one-way wave system representing downgoing waves is used for a modified reverse time migration method. As a wavefield extrapolator in migration, the downgoing wave system propagates the reflection events backwards to their reflectors without scattering at the discontinuities in the velocity model. Hence, images with amplitudes proportional to reflectivity can be obtained from this migration technique. We present examples of the application of the new migration method to synthetic seismic data where P-P reflections P-SV converted waves are present. Absorbing boundaries, useful in the generation of synthetic seismograms, have been constructed by using the one-way wave system. These boundaries absorb effectively waves impinging over a wide range of angles of incidence. 相似文献
12.
The Transition Between the Scale Domains of Ray and Effective Medium Theory and Anisotropy: Numerical Models 总被引:1,自引:0,他引:1
The anisotropy of a periodically layered isotropic medium is numerically modeled in order to study the effect of the scale
of heterogeneity on seismic observations. An important motivation is to delineate the wavelength ranges over which a pulse
propagating obliquely through the structure will be described by either ray (short wavelength) or effective medium (long wavelength)
theory. The same band-limited pulse is propagated obliquely at a variety of incidence angles through a compositionally uniform
layered structure as a function of the layer thicknesses. The resulting seismograms display similar behavior to that encountered
for normal incidence including the effects of stop- and pass-bands. Velocities determined from time picks on these seismograms
show a large difference in velocities between the long and short wavelength limits as has been previously demonstrated for
normal incidence propagation. The bulk of the transition between these two limits is independent of incidence angle and occurs
when the ratio between the wavelength and the layering thickness is near a value of 10. Two more geologically reasonable models
show that these effects are diminished with smaller contrasts between the layers. 相似文献
13.
In order to trace a ray between known source and receiver locations in a perfectly elastic medium, the take-off angle must be determined, or equialently, the ray parameter. In a viscoelastic medium, the initial value of a second angle, the attenuation angle (the angle between the normal to the plane wavefront and the direction of maximum attenuation), must also be determined. There seems to be no agreement in the literature as to how this should be done. In computing anelastic synthetic seismograms, some authors have simply chosen arbitrary numerical values for the initial attenuation angle, resulting in different raypaths for different choices. There exists, however, a procedure in which the arbitrariness is not present, i.e., in which the raypath is uniquely determined. It consists of computing the value of the anelastic ray parameter for which the phase function is stationary (Fermat's principle). This unique value of the ray parameter gives unique values for the take-off and attenuation angles. The coordinates of points on these stationary raypaths are complex numbers. Such rays are known as complex rays. They have been used to study electromagnetic wave propagation in lossy media. However, ray-synthetic seismograms can be computed by this procedure without concern for the details of complex raypath coordinates. To clarify the nature of complex rays, we study two examples involving a ray passing through a vertically inhomogeneous medium. In the first example, the medium consists of a sequence of discrete homogeneous layers. We find that the coordinates of points on the ray are generally complex (other than the source and receiver points which are usually assumed to lie in real space), except for a ray which is symmetric about an axis down its center, in which case the center point of the ray lies in real space. In the second example, the velocity varies continuously and linearly with depth. We show that, in geneneral, the turning point of the ray lies in complex space (unlike the symmetric ray in the discrete layer case), except if the ratio of the velocity gradient to the complex frequency-dependent velocity at the surface is a real number. We also present a numerical example which demonstrates that the differences between parameters, such as arrival time and raypath angles, for the stationary ray and for rays computed by the above-mentioned arbitrary approaches can be substantial. 相似文献
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.
Xiaoshan Wang Xiangdong Feng Guiling Diao Libing Wang Lingling Cai Earthquake Administration of Hebei Province Shijiazhuang China Hongshan Benchmark Seismic Station Longyao China 《地震学报(英文版)》2009,(6):667-673
Using the maximum amplitude ratios of vertical component of P and S waves recorded by a regional network, 921 focal mechanisms of Dayao earthquake doublet sequence are determined by means of synthetic seismograms of a point source of dislocation in a plane layered medium.Among them,389 focal mechanisms are in the aftershock sequence of M6.2 earthquake occurred on 21 July,2003 and the other 532 focal mechanisms are in the aftershock sequence of M6.1 earthquake occurred on 16 October,2003 in Dayao,Yunnan.The ... 相似文献
16.
Summary The hybrid ray-reflectivity method is applied to the problem of the transmission of the reflected wave field through a thin high-velocity layer (or through a thin stack of high velocity layers), situated in the overburden of the reflector. In the hybrid ray-reflectivity method, the standard ray method is applied in the smooth parts of the model, and the reflectivity method is used locally at the thin high-velocity layer. With the exception of small epicentral distances, the standard ray method itself fails in such computations. The reason is that a considerable part of the energy for overcritical angles of incidence may be tunneled through the thin high-velocity layer along complex ray-paths, corresponding to inhomogeneous waves. The reflectivity method, applied locally at the thin high-velocity layer, automatically includes all inhomogeneous wave contributions. Thus, the hybrid ray-reflectivity method removes fully the limitations of the standard ray method, but still retains its main advantages, such as its applicability to 2-D and 3-D complex layered structures, flexibility, and low-cost computations. In the numerical examples, the hybrid ray-reflectivity synthetic seismograms are compared with standard ray synthetic seismograms and with full reflectivity computations. The numerical examples show that the hybrid ray-reflectivity method describes the tunneling of seismic energy through a thin high-velocity layer with sufficient accuracy. 相似文献
17.
18.
19.
20.
Synthetic seismograms for finite sources in spherically symmetric Earth using normal-mode summation 下载免费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. 相似文献