共查询到20条相似文献，搜索用时 31 毫秒
1.
Traveltime computation by wavefront-orientated ray tracing 总被引：1，自引：0，他引：1
For multivalued traveltime computation on dense grids, we propose a wavefront‐orientated ray‐tracing (WRT) technique. At the source, we start with a few rays which are propagated stepwise through a smooth two‐dimensional (2D) velocity model. The ray field is examined at wavefronts and a new ray might be inserted between two adjacent rays if one of the following criteria is satisfied: (1) the distance between the two rays is larger than a predefined threshold; (2) the difference in wavefront curvature between the rays is larger than a predefined threshold; (3) the adjacent rays intersect. The last two criteria may lead to oversampling by rays in caustic regions. To avoid this oversampling, we do not insert a ray if the distance between adjacent rays is smaller than a predefined threshold. We insert the new ray by tracing it from the source. This approach leads to an improved accuracy compared with the insertion of a new ray by interpolation, which is the method usually applied in wavefront construction. The traveltimes computed along the rays are used for the estimation of traveltimes on a rectangular grid. This estimation is carried out within a region bounded by adjacent wavefronts and rays. As for the insertion criterion, we consider the wavefront curvature and extrapolate the traveltimes, up to the second order, from the intersection points between rays and wavefronts to a gridpoint. The extrapolated values are weighted with respect to the distances to wavefronts and rays. Because dynamic ray tracing is not applied, we approximate the wavefront curvature at a given point using the slowness vector at this point and an adjacent point on the same wavefront. The efficiency of the WRT technique is strongly dependent on the input parameters which control the wavefront and ray densities. On the basis of traveltimes computed in a smoothed Marmousi model, we analyse these dependences and suggest some rules for a correct choice of input parameters. With suitable input parameters, the WRT technique allows an accurate traveltime computation using a small number of rays and wavefronts. 相似文献
2.
为了解决三维复杂介质的射线追踪,本文改变了传统的三维层状地层的建模描述方式,提出了块状结构的建模描述方法,结合三角形面片来描述地质界面,可以构造非常复杂的三维地质模型．为了满足射线追踪的需要,本文对模型界面内的法向量进行光滑处理,光滑后的法向量在界面内是连续变化的．在块状模型的基础上,本文运用三角形的面积坐标,提出了几种试射角度的修正方法：细分三角形法、分割三角形法和子三角形法,计算表明子三角形法最好．文中给出了三维块状模型和射线追踪实例． 相似文献
3.
A new 3D wavefield modelling approach based on dynamic ray tracing is presented. This approach is called wavefront construction, and it can be used in 3D models with constant or smoothly varying material properties (S- and P-velocity and density) separated by smooth interfaces. Wavefronts consisting of rays arranged in a triangular network are propagated stepwise through the model. At each time step, the differences in a number of parameters are checked between each pair of rays on the wavefront. New rays are interpolated whenever this difference between pairs of rays exceeds some predefined maximum value. A controlled sampling of the wavefront at all time steps is thus obtained. Receivers are given multiple-event values by interpolation when the wavefronts pass them. The strength of the wavefront construction method is that it is robust and efficient. 相似文献
4.
The 4 × 4 T -propagator matrix of a 3D central ray determines, among other important seismic quantities, second-order (parabolic or hyperbolic) two-point traveltime approximations of certain paraxial rays in the vicinity of the known central ray through a 3D medium consisting of inhomogeneous isotropic velocity layers. These rays result from perturbing the start and endpoints of the central ray on smoothly curved anterior and posterior surfaces. The perturbation of each ray endpoint is described only by a two-component vector. Here, we provide parabolic and hyperbolic paraxial two-point traveltime approximations using the T -propagator to feature a number of useful 3D seismic models, putting particular emphasis on expressing the traveltimes for paraxial primary reflected rays in terms of hyperbolic approximations. These are of use in solving several forward and inverse seismic problems. Our results simplify those in which the perturbation of the ray endpoints upon a curved interface is described by a three-component vector. In order to emphasize the importance of the hyperbolic expression, we show that the hyperbolic paraxial-ray traveltime (in terms of four independent variables) is exact for the case of a primary ray reflected from a planar dipping interface below a homogeneous velocity medium. 相似文献
5.
The behaviour of the actual polarization of an electromagnetic wave or elastic S–wave is described by the coupling ray theory, which represents the generalization of both the zero–order isotropic and anisotropic ray theories and provides continuous transition between them. The coupling ray theory is usually applied to anisotropic common reference rays, but it is more accurate if it is applied to reference rays which are closer to the actual wave paths. In a generally anisotropic or bianisotropic medium, the actual wave paths may be approximated by the anisotropic–ray–theory rays if these rays behave reasonably. In an approximately uniaxial (approximately transversely isotropic) anisotropic medium, we can define and trace the SH (ordinary) and SV (extraordinary) reference rays, and use them as reference rays for the prevailing–frequency approximation of the coupling ray theory. In both cases, i.e. for the anisotropic–ray–theory rays or the SH and SV reference rays, we have two sets of reference rays. We thus obtain two arrivals along each reference ray of the first set and have to select the correct one. Analogously, we obtain two arrivals along each reference ray of the second set and have to select the correct one. In this paper, we suggest the way of selecting the correct arrivals. We then demonstrate the accuracy of the resulting prevailing–frequency approximation of the coupling ray theory using elastic S waves along the SH and SV reference rays in four different approximately uniaxial (approximately transversely isotropic) velocity models. 相似文献
6.
We develop the true‐amplitude prestack migration of multicomponent data based on the use of elastic Gaussian beams for walkaway vertical seismic profile (VSP) acquisition systems. It consists in a weighted summation of multishot data with specific weights, computed by tracing elastic Gaussian beams from each imaging point of the target area towards the sources and receivers. Each pair of beams may be connected with either a pair of P‐rays (PP‐image) or the P‐ray towards sources and the S‐ray to receivers (PS‐image) and is uniquely determined by dip (the angle of the bisector between the rays and the vertical direction) and opening (the angle between the rays) angles. Shooting from the bottom towards the acquisition system helps to avoid well‐known troubles, in particular multipathing for the imaging conditions in complex velocity models. The ability to fix the dip angle and implement summation over opening angles leads to the so‐called selective images that contain mostly interfaces with desired slopes. On the other hand, a set of images computed for a range of opening angles by summation over all available dip angles is used as input of an AVO‐like inversion procedure for the recovery of elastic parameters. The feasibility of this imaging procedure is verified by synthetic data for 2D realistic elastic models. 相似文献
7.
A first-order perturbation theory for seismic isochrons is presented in a model independent form. Two ray concepts are fundamental in this theory, the isochron ray and the velocity ray, for which I obtain first-order approximations to position vectors and slowness vectors. Furthermore, isochron points are connected to a shot and receiver by conventional ray fields. Based on independent perturbation of the shot and receiver ray I obtain first-order approximations to velocity rays. The theory is applicable for 3D inhomogeneous anisotropic media, given that the shot and receiver rays, as well as their perturbations, can be generated with such model generality.
The theory has applications in sensitivity analysis of prestack depth migration and in velocity model updating. Numerical examples of isochron and velocity rays are shown for a 2D homogeneous VTI model. The general impression is that the first-order approximation is, with some exceptions, sufficiently accurate for practical applications using an anisotropic velocity model. 相似文献
8.
3D multivalued travel time and amplitude maps 总被引：2，自引：0，他引：2
An algorithm for computing multivalued maps for travel time, amplitude and any other ray related variable in 3D smooth velocity models is presented. It is based on the construction of successive isochrons by tracing a uniformly dense discrete set of rays by fixed travel-time steps. Ray tracing is based on Hamiltonian formulation and includes computation of paraxial matrices. A ray density criterion ensures uniform ray density along isochrons over the entire ray field including caustics. Applications to complex models are shown. 相似文献
9.
The coupling ray theory is usually applied to anisotropic common reference rays, but it is more accurate if it is applied to reference rays which are closer to the actual wave paths. If we know that a medium is close to uniaxial (transversely isotropic), it may be advantageous to trace reference rays which resemble the SH–wave and SV–wave rays. This paper is devoted to defining and tracing these SH and SV reference rays of elastic S waves in a heterogeneous generally anisotropic medium which is approximately uniaxial (approximately transversely isotropic), and to the corresponding equations of geodesic deviation (dynamic ray tracing). All presented equations are simultaneously applicable to ordinary and extraordinary reference rays of electromagnetic waves in a generally bianisotropic medium which is approximately uniaxially anisotropic. The improvement of the coupling–ray–theory seismograms calculated along the proposed SH and SV reference rays, compared to the coupling–ray–theory seismograms calculated along the anisotropic common reference rays, has already been numerically demonstrated by the authors in four approximately uniaxial velocity models. 相似文献
10.
Petr Bulant 《Pure and Applied Geophysics》1996,148(3-4):421-447
Algorithm for determination of all two-point rays of a given elementary wave by means of the shooting method is presented. The algorithm is designed for general 3-D models composed of inhomogeneous geological blocks separated by curved interfaces. It is independent of the initial conditions for rays and of the initial-value ray tracer. The algorithm described has been coded in Fortran 77, using subroutine packages MODEL and CRT for model specification and for initial-value ray tracing. 相似文献
11.
Diffraction and anelasticity problems involving decaying, “evanescent” or “inhomogeneous” waves can be studied and modelled using the notion of “complex rays”. The wavefront or “eikonal” equation for such waves is in general complex and leads to rays in complex position-slowness space. Initial conditions must be specified in that domain: for example, even for a wave originating in a perfectly elastic region, the ray to a real receiver in a neighbouring anelastic region generally departs from a complex point on the initial-values surface. Complex ray theory is the formal extension of the usual Hamilton equations to complex domains. Liouville's phase-space-incompressibility theorem and Fermat's stationary-time principle are formally unchanged. However, an infinity of paths exists between two fixed points in complex space all of which give the same final slowness, travel time, amplitude, etc. This does not contradict the fact that for a given receiver position there is a unique point on the initial-values surface from which this infinite complex ray family emanates.In perfectly elastic media complex rays are associated with, for example, evanescent waves in the shadow of a caustic. More generally, caustics in anelastic media may lie just outside the real coordinate subspace and one must trace complex rays around the complex caustic in order to obtain accurate waveforms nearby or the turning waves at greater distances into the lit region. The complex extension of the Maslov method for computing such waveforms is described. It uses the complex extension of the Legendre transformation and the extra freedom of complex rays makes pseudocaustics avoidable. There is no need to introduce a Maslov/KMAH index to account for caustics in the geometrical ray approximation, the complex amplitude being generally continuous. Other singular ray problems, such as the strong coupling around acoustic axes in anisotropic media, may also be addressed using complex rays.Complex rays are insightful and practical for simple models (e.g. homogeneous layers). For more complicated numerical work, though, it would be desirable to confine attention to real position coordinates. Furthermore, anelasticity implies dispersion so that complex rays are generally frequency dependent. The concept of group velocity as the velocity of a spatial or temporal maximum of a narrow-band wave packet does lead to real ray/Hamilton equations. However, envelope-maximum tracking does not itself yield enough information to compute synthetic seismogramsFor anelasticity which is weak in certain precise senses, one can set up a theory of real, dispersive wave-packet tracking suitable for synthetic seismogram calculations in linearly visco-elastic media. The seismologically-accepiable constant-Q rheology of Liu et al. (1976), for example, satisfies the requirements of this wave-packet theory, which is adapted from electromagnetics and presented as a reasonable physical and mathematical basis for ray modelling in inhomogeneous, anisotropic, anelastic media. Dispersion means that one may need to do more work than for elastic media. However, one can envisage perturbation analyses based on the ray theory presented here, as well as extensions like Maslov's which are based on the Hamiltonian properties. 相似文献
12.
当前的三维VSP地震数据偏移成像实现都是在共炮点道集或共检波点道集中逐个道集循环进行的,计算效率相对较低.根据三维VSP观测系统中炮点和检波点布置的特殊性和地震波场满足线性叠加的特性,本文提出了一种三维VSP数据的高效偏移成像方法,即首先通过对三维VSP共接收点道集进行地震数据的广义合成得到一种超道集,然后在共接收点道集的波场深度外推过程中逐步应用多震源波场对超道集进行偏移成像,即利用一次波场深度外推循环完成对所有共检波点道集数据的偏移成像.通过三维VSP模型数据与实际地震数据的偏移成像试验验证了这种高效的超道集偏移成像方法可取得与常规共检波点道集相当的偏移成像效果,还具有极高的计算效率,其计算量与单个共检波点道集的偏移成像计算量相当. 相似文献
13.
B. ARNTSEN 《Geophysical Prospecting》1988,36(5):478-503
Ray theories are a class of methods often chosen to compute synthetic seismograms due to their efficiency and ability to deal with complex, three-dimensional inhomogeneous media. To deal with the large number of rays needed to compute synthetic seismograms, a ray generation algorithm is given which is capable of generating a numerical code describing each ray. The code describes a subset of all possible rays by considering only pre-critical reflections. In a horizontally plane-layered medium the generation of rays and computation of amplitudes and traveltimes can be efficiently accomplished by grouping the rays into reflection order and dynamic analogue groups. Expressions summing all unconverted rays and rays with a single mode conversion are given for source and receiver located at arbitrary positions within the medium. Examples of zero-offset synthetic VSPs obtained by this method are given. 相似文献
14.
Emil Blias 《Geophysical Prospecting》2013,61(3):574-581
I introduce a new explicit form of vertical seismic profile (VSP) traveltime approximation for a 2D model with non‐horizontal boundaries and anisotropic layers. The goal of the new approximation is to dramatically decrease the cost of time calculations by reducing the number of calculated rays in a complex multi‐layered anisotropic model for VSP walkaway data with many sources. This traveltime approximation extends the generalized moveout approximation proposed by Fomel and Stovas. The new equation is designed for borehole seismic geometry where the receivers are placed in a well while the sources are on the surface. For this, the time‐offset function is presented as a sum of odd and even functions. Coefficients in this approximation are determined by calculating the traveltime and its first‐ and second‐order derivatives at five specific rays. Once these coefficients are determined, the traveltimes at other rays are calculated by this approximation. Testing this new approximation on a 2D anisotropic model with dipping boundaries shows its very high accuracy for offsets three times the reflector depths. The new approximation can be used for 2D anisotropic models with tilted symmetry axes for practical VSP geometry calculations. The new explicit approximation eliminates the need of massive ray tracing in a complicated velocity model for multi‐source VSP surveys. This method is designed not for NMO correction but for replacing conventional ray tracing for time calculations. 相似文献
15.
Seismic ray path variations in a 3D global velocity model 总被引：2，自引：0，他引：2
A three-dimensional (3D) ray tracing technique is used to investigate ray path variations of P, PcP, pP and PP phases in a global tomographic model with P wave velocity changing in three dimensions and with lateral depth variations of the Moho, 410 and 660 km discontinuities. The results show that ray paths in the 3D velocity model deviate considerably from those in the average 1D model. For a PcP wave in Western Pacific to East Asia where the high-velocity (1-2%) Pacific slab is subducting beneath the Eurasian continent, the ray path change amounts to 27 km. For a PcP ray in South Pacific where very slow (−2%) velocity anomalies (the Pacific superplume) exist in the whole mantle, the maximum ray path deviation amounts to 77 km. Ray paths of other phases (P, pP, PP) are also displaced by tens of kilometers. Changes in travel time are as large as 3.9 s. These results suggest that although the maximal velocity anomalies of the global tomographic model are only 1-2%, rays passing through regions with strong lateral heterogeneity (in velocity and/or discontinuity topography) can have significant deviations from those in a 1D model because rays have very long trajectories in the global case. If the blocks or grid nodes adopted for inversion are relatively large (3-5°) and only a low-resolution 3D model is estimated, 1D ray tracing may be feasible. But if fine blocks or grid nodes are used to determine a high-resolution model, 3D ray tracing becomes necessary and important for the global tomography. 相似文献
16.
To calculate the hydrodynamic interaction forces of the reservoir directly in the time-domain, the dynamic stiffness of each mode of the semi-infinite uniform fluid channel is either represented by a lumped-parameter model with frequency-independent real coefficients of the springs, dashpots and masses and with only a few additional internal degrees of freedom, or the interaction forces are calculated recursively. For each mode characterized by its eigenvalue, the coefficients of the lumped-parameter model and the recursive coefficients are specified, which can be used directly in a practical application. The procedures exhibit many advantages: the only approximation (replacing the rigorous dynamic stiffness by a ratio of two polynomials) can be evaluated visibly. No unfamiliar discrete-time manipulations such as the z-transformation are used. The stiffness, damping and mass matrices corresponding to the lumped-parameter model are automatically symmetrical. Stability of the procedures is also guaranteed. Combining the lumped-parameter model of the semi-infinite uniform channel with the finite-element discretization of the irregular fluid region or calculating the interaction forces recursively allows a reservoir of arbitrary shape to be analysed directly in the time domain. Non-linearities in the dam can, thus, be taken into consideration in a seismic analysis. 相似文献
17.
Offset continuation (OCO) is a seismic configuration transform designed to simulate a seismic section as if obtained with
a certain source-receiver offset using the data measured with another offset. Since OCO is dependent on the velocity model
used in the process, comparison of the simulated section to an acquired section allows for the extraction of velocity information.
An algorithm for such a horizon-oriented velocity analysis is based on so-called OCO rays. These OCO rays describe the output
point of an OCO as a function of the Root Mean Square (RMS) velocity. The intersection point of an OCO ray with the picked
traveltime curve in the acquired data corresponding to the output half-offset defines the RMS velocity at that position. We
theoretically relate the OCO rays to the kinematic properties of OCO image waves that describe the continuous transformation
of the common-offset reflection event from one offset to another. By applying the method of characteristics to the OCO image-wave
equation, we obtain a raytracing-like procedure that allows to construct OCO trajectories describing the position of the OCO
output point under varying offset. The endpoints of these OCO trajectories for a single input point and different values of
the RMS velocity form then the OCO rays. A numerical example demonstrates that the developed ray-tracing procedure leads to
reliable OCO rays, which in turn provide high-quality RMS velocities. The proposed procedure can be carried out fully automatically,
while conventional velocity analysis needs human intervention. Moreover, since velocities are extracted using offset sections,
more redundancy is available or, alternatively, OCO velocities can be studied as a function of offset. 相似文献
18.
Roberto de Franco 《Geophysical Prospecting》2001,49(4):395-404
A method to estimate interval velocities and thickness in a horizontal isotropic layered medium from wide-angle reflection traveltime curves is presented. The method is based on a relationship between the squared reflection traveltime differences and the squared offset differences relative to two adjacent reflectors. The envelope of the squared-time versus offset-difference curves, for rays with the same ray parameter, is a straight line, whose slope is the inverse of the square of the interval velocity and whose intercept is the square of the interval time. The method yields velocity and thickness estimates without any knowledge of the overlying stratification. It can be applied to wide-angle reflection data when either information on the upper crust and/or refraction control on the velocity is not available. Application to synthetic and real data shows that the method, used together with other methods, allows us to define a reliable 1D starting model for estimating a depth profile using either ray tracing or another technique. 相似文献
19.
M. ALI AK 《Geophysical Prospecting》1990,38(8):971-982
Seismic refraction surveying is still an important tool for determining the geometries and elastic wave propagation velocities of near-surface layers. Many analytical and graphical methods have been developed over the years for refraction interpretation, and these can be classified into two basic groups. The first group visualizes critically refracted rays converging on a common surface position, while the second group, which includes the wavefront methods, makes use of the critical rays emerging from a common point on the refractor. The method described in this paper is an analytical approach to the wavefront methods. The reverse refracted ray received by a geophone is intersected by the forward refracted rays received by subsequent geophones and a common critical refraction point on the refractor is estimated after a series of comparisons. This process is repeated for each geophone to yield the geometry and the velocity of the refractor. Several interpolations are performed to achieve a better accuracy. Palmer's models are used to test the efficiency of the algorithm. The results are presented together with those of other methods applied to the same models. 相似文献
20.
KLAUS HELBIG 《Geophysical Prospecting》1990,38(2):189-220
Wavefront charts in anisotropic gradient media are a useful tool in ray geometric constructions, particular in shear-wave exploration. They can be constructed by: (i) a family of wavefronts that contains a vertical plane as member - it is convenient to choose constant time increments; (ii) tracing one ray that makes everywhere the angle with the normal to the wavefront that is required by the anisotropy of the medium; (iii) scaling this ray to obtain a set of rays with different ray parameters; (iv) shifting these rays (with wavefront elements attached) so that they pass through a common source point; (v) interpolating the wavefronts between the elements. The construction is particularly simple in linear-gradient media, since here all members of the family of wavefronts are planes. Since the ray makes everywhere the angle prescribed by the anisotropy with the normal of the (plane) wavefronts, the ray has the shape of the slowness curve rotated by ?π/2. For isotropic media the slowness curve is a circle, and thus rays are circular arcs. The circles themselves intersect in the source point and in a second point above the surface of the earth. This provides a simple proof that wavefronts emanating from a point source in an isotropic linear-gradient medium are spheres: inversion of the set of circular rays with the source as centre maps the pencil of circular rays into a pencil of straight lines passing through a point. A pencil of concentric spheres around this point is perpendicular to the pencil of straight lines. On inverting back the pencil of spheres is mapped into another pencil of spheres that is perpendicular to the circular rays. 相似文献