共查询到20条相似文献，搜索用时 140 毫秒
1.
Kinematical characteristics of reflected waves in anisotropic elastic media play an important role in the seismic imaging workflow. Considering compressional and converted waves, we derive new, azimuthally dependent, slowness-domain approximations for the kinematical characteristics of reflected waves (radial and transverse offsets, intercept time and traveltime) for layered orthorhombic media with varying azimuth of the vertical symmetry planes. The proposed method can be considered an extension of the well-known ‘generalized moveout approximation’ in the slowness domain, from azimuthally isotropic to azimuthally anisotropic models. For each slowness azimuth, the approximations hold for a wide angle range, combining power series coefficients in the vicinity of both the normal-incidence ray and an additional wide-angle ray. We consider two cases for the wide-angle ray: a ‘critical slowness match’ and a ‘pre-critical slowness match’ studied in Parts I and II of this work, respectively. For the critical slowness match, the approximations are valid within the entire slowness range, up to the critical slowness. For the ‘pre-critical slowness match’, the approximations are valid only within the bounded slowness range; however, the accuracy within the defined range is higher. The critical slowness match is particularly effective when the subsurface model includes a dominant high-velocity layer where, for nearly critical slowness values, the propagation in this layer is almost horizontal. Comparing the approximated kinematical characteristics with those computed by numerical ray tracing, we demonstrate high accuracy. 相似文献
2.
3.
The well‐known asymptotic fractional four‐parameter traveltime approximation and the five‐parameter generalised traveltime approximation in stratified multi‐layer transversely isotropic elastic media with a vertical axis of symmetry have been widely used for pure‐mode and converted waves. The first three parameters of these traveltime expansions are zero‐offset traveltime, normal moveout velocity, and quartic coefficient, ensuring high accuracy of traveltimes at short offsets. The additional parameter within the four‐parameter approximation is an effective horizontal velocity accounting for large offsets, which is important to avoid traveltime divergence at large offsets. The two additional parameters in the above‐mentioned five‐parameter approximation ensure higher accuracy up to a given large finite offset with an exact match at this offset. In this paper, we propose two alternative five‐parameter traveltime approximations, which can be considered extensions of the four‐parameter approximation and an alternative to the five‐parameter approximation previously mentioned. The first three short‐offset parameters are the same as before, but the two additional long‐offset parameters are different and have specific physical meaning. One of them describes the propagation in the high‐velocity layer of the overburden (nearly horizontal propagation in the case of very large offsets), and the other characterises the intercept time corresponding to the critical slowness that includes contributions of the lower velocity layers only. Unlike the above‐mentioned approximations, both of the proposed traveltime approximations converge to the theoretical (asymptotic) linear traveltime at the limit case of very large (“infinite”) offsets. Their accuracy for moderate to very large offsets, for quasi‐compressional waves, converted waves, and shear waves polarised in the horizontal plane, is extremely high in cases where the overburden model contains at least one layer with a dominant higher velocity compared with the other layers. We consider the implementation of the proposed traveltime approximations in all classes of problems in which the above‐mentioned approximations are used, such as reflection and diffraction analysis and imaging. 相似文献
4.
Slowness domain offset and traveltime approximations in layered vertical transversely isotropic media 
下载免费PDF全文
下载免费PDF全文 Considering horizontally layered transversely isotropic media with vertical symmetry axis and all types of pure‐mode and converted waves we present a new wide‐angle series approximation for the kinematical characteristics of reflected waves: horizontal offset, intercept time, and total reflection traveltime as functions of horizontal slowness. The method is based on combining (gluing) both zero‐offset and (large) finite‐offset series coefficients. The horizontal slowness is bounded by the critical value, characterised by nearly horizontal propagation within the layer with the highest horizontal velocity. The suggested approximation uses five parameters to approximate the offset, six parameters to approximate the intercept time or the traveltime, and seven parameters to approximate any two or all three kinematical characteristics. Overall, the method is very accurate for pure‐mode compressional waves and shear waves polarised in the horizontal plane and for converted waves. The application of the method to pure‐mode shear waves polarised in the vertical plane is limited due to cusps and triplications. To demonstrate the high accuracy of the method, we consider a synthetic, multi‐layer model, and we plot the normalised errors with respect to numerical ray tracing. 相似文献
5.
6.
Large-offset approximation to seismic reflection traveltimes 总被引：4，自引：0，他引：4
Conventional approximations of reflection traveltimes assume a small offset-to-depth ratio, and their accuracy decreases with increasing offset-to-depth ratio. Hence, they are not suitable for velocity analysis and stacking of long-offset reflection seismic data. Assuming that the offset is large, rather than small, we present a new traveltime approximation which is exact at infinite offset and has a decreasing accuracy with decreasing offset-to-depth ratio. This approximation has the form of a series containing powers of the offset from 1 to −∞. It is particularly accurate in the presence of a thin high-velocity layer above the reflector, i.e. in a situation where the accuracy of the Taner and Koehler series is poor. This new series can be used to gain insight into the velocity information contained in reflection traveltimes at large offsets, and possibly to improve velocity analysis and stacking of long-offset reflection seismic data. 相似文献
7.
Anisotropy in subsurface geological models is primarily caused by two factors: sedimentation in shale/sand layers and fractures. The sedimentation factor is mainly modelled by vertical transverse isotropy (VTI), whereas the fractures are modelled by a horizontal transversely isotropic medium (HTI). In this paper we study hyperbolic and non‐hyperbolic normal reflection moveout for a package of HTI/VTI layers, considering arbitrary azimuthal orientation of the symmetry axis at each HTI layer. We consider a local 1D medium, whose properties change vertically, with flat interfaces between the layers. In this case, the horizontal slowness is preserved; thus, the azimuth of the phase velocity is the same for all layers of the package. In general, however, the azimuth of the ray velocity differs from the azimuth of the phase velocity. The ray azimuth depends on the layer properties and may be different for each layer. In this case, the use of the Dix equation requires projection of the moveout velocity of each layer on the phase plane. We derive an accurate equation for hyperbolic and high‐order terms of the normal moveout, relating the traveltime to the surface offset, or alternatively, to the subsurface reflection angle. We relate the azimuth of the surface offset to its magnitude (or to the reflection angle), considering short and long offsets. We compare the derived approximations with analytical ray tracing. 相似文献
8.
A fast and robust method for two-point ray tracing in one-dimensional layered media is presented. This method is applicable to layered models with constant or linearly varying isotropic layer velocity. For given model properties and source and receiver positions, a ray path can be uniquely determined once its ray parameter (i.e. horizontal slowness) is known. The ray parameter can be obtained by numerically solving the nonlinear offset (i.e. source–receiver horizontal distance) equation using Newton's method, which generally works well at near and mid offsets. However, Newton's method becomes hard to converge at large offsets due to the oversensitivity of offset to ray parameter. Based on the analysis of the characteristic of the offset equation, a modified ray parameter is proposed and used to replace the generic ray parameter in numerical calculation. Numerical experiments show that the iteration process becomes stable and converges rapidly with the modified ray parameter. Moreover, a rational function that asymptotically approximates the shape of the offset equation is introduced for obtaining good initial estimates of the modified ray parameter. Numerical tests show that this method is robust in any situation, and an accurate ray parameter can be obtained within two or three iterations for a wide range of model velocity structure and source–receiver distance. Furthermore, the proposed two-point ray tracing method is easy to implement. 相似文献
9.
Seismic amplitude variations with offset contain information about the elastic parameters. Prestack amplitude analysis seeks to extract this information by using the variations of the reflection coefficients as functions of angle of incidence. Normally, an approximate formula is used for the reflection coefficients, and variations with offset of the geometrical spreading and the anelastic attenuation are often ignored. Using angle of incidence as the dependent variable is also computationally inefficient since the data are recorded as a function of offset. Improved approximations have been derived for the elastic reflection and transmission coefficients, the geometrical spreading and the complex travel-time (including anelastic attenuation). For a 1 D medium, these approximations are combined to produce seismic reflection amplitudes (P-wave, S-wave or converted wave) as a Taylor series in the offset coordinate. The coefficients of the Taylor series are computed directly from the parameters of the medium, without using the ray parameter. For primary reflected P-waves, dynamic ray tracing has been used to compute the offset variations of the transmission coefficients, the reflection coefficient, the geometrical spreading and the anelastic attenuation. The offset variation of the transmission factor is small, while the variations in the geometrical spreading, absorption and reflection coefficient are all significant. The new approximations have been used for seismic modelling without ray tracing. The amplitude was approximated by a fourth-order polynomial in offset, the traveltime by the normal square-root approximation and the absorption factor by a similar expression. This approximate modelling was compared to dynamic ray tracing, and the results are the same for zero offset and very close for offsets less than the reflector depth. 相似文献
10.
Using higher-order ray theory, we derived exact elastodynamic Green functions for three simple types of homogeneous anisotropy. The first type displays an orthorhombic symmetry, the other two types display transverse isotropy. In all cases, the slowness surfaces of waves are either ellipsoids, spheroids or spheres. All three Green functions are expressed by a ray series with a finite number of terms. The Green functions can be written in explicit and elementary form similar to the Stokes solution for isotropy. In two Green functions, the higher-order ray approximations form a near-singularity term, which is significant near a kiss singularity. In the third Green function, the higher-order ray approximations also form a near-field term, which is significant near the point source. No effect connected with the line singularity was observed. 相似文献
11.
OBC地震物理模型实验中的广角反射现象(英文) 总被引：1，自引：1，他引：0
广角地震采集作为获取野外高质量地震数据的一种重要方法,愈来愈受到人们的关注。因此,对广角反射的深入研究具有重要实际意义。本文通过层状物理模型海底电缆(0BC)地震勘探实验研究,比较详细地讨论了广角反射现象。一些实验结果并不支持理论数值模拟结论。主要实验结论有:1)广角反射波的振幅能量较非广角反射波强(约为1倍多),但和理论计算的15倍关系差异甚大,过临界角时反射能量缓慢增大,而不是理论计算的那样急剧增大;2)大炮检距时,反射同相轴仍保持近似双曲线趋势;3)广角反射波的主频较非广角反射波低(降低20-30%),且偏移距增大有所降低,而非广角反射在小偏移距段主频变化不明显;4)在临界角前后反射波波形没有发生突变,极性无变化;5)反射波组特征在临界角前后也未发生变化。6)水中直达波、多次波和水底折射波对海底层广角反射有一定的影响。 相似文献
12.
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. 相似文献
13.
Einar Iversen 《Pure and Applied Geophysics》1996,148(1-2):287-317
The motivation for this paper is to provide expressions for first-order partial derivatives of reflection point coordinates, taken with respect to model parameters. Such derivatives are expected to be useful for processes dealing with the problem of estimating velocities for depth migration of seismic data.The subject of the paper is a particular aspect of ray perturbation theory, where observed parameters—two-way reflection time and horizontal components of slowness, are constraining the ray path when parameters of the reference velocity model are perturbed. The methodology described here is applicable to general rays in a 3D isotropic, heterogeneous medium. Each ray is divided into a shot ray and a receiver ray, i.e., the ray portions between the shot/receiver and the reflection point, respectively. Furthermore, by freezing the initial horizontal slowness of these subrays as the model is perturbed,elementary perturbation quantities may be obtained, comprising derivatives of ray hit positions within theisochrone tangent plane, as well as corresponding time derivatives. The elementary quantities may be estimated numerically, by use of ray perturbation theory, or in some cases, analytically. In particular, when the layer above the reflection point is homogeneous, explicit formulas can be derived. When the elementary quantities are known,reflection point derivatives can be obtained efficiently from a set of linear expressions.The method is applicable for a common shot, receiver or offset data sorting. For these gather types, reflection point perturbationlaterally with respect to the isochrone is essentially different. However, in theperpendicular direction, a first-order perturbation is shown to beindependent of gather type.To evaluate the theory, reflection point derivatives were estimated analytically and numerically. I also compared first-order approximations to true reflection point curves, obtained by retracing rays for a number of model perturbations. The results are promising, especially with respect to applications in sensitivity analysis for prestack depth migration and in velocity model updating. 相似文献
14.
Refracted arrivals are analysed to estimate the near‐surface anisotropy of marine sediments using a vertical‐cable (VC) configuration. In the presence of dip, the horizontal and vertical ray‐slownesses are obtained from the observed apparent slownesses in the up‐ and downdip directions using a sum or difference at each azimuth. The multiple azimuths generated by a VC geometry permit the ray‐slowness distribution of the marine sediments to be determined. An inversion procedure is developed to provide dip and anisotropy parameters for refractive layers from the measured refraction traveltimes in multilayered azimuthally isotropic and anisotropic media. Two sets of transversely isotropic models are used to analyse the azimuthal variations of apparent and ray slownesses. In the first set, we fix the anisotropic parameters of the models but vary the dip (0°, 5° and 10°) to test the effects of the presence of dip. In the second set, we vary the P‐wave anisotropy strength (5.2%, 10.3%, 15.8% and 22.0%) to examine the sensitivity and accuracy of ray‐slowness approximations which are independent of dip. We test this inversion procedure on synthetic P‐wave VC data calculated for six different models by a finite‐difference method. The results of applications to real VC data acquired from the North Sea are also presented. 相似文献
15.
We study the azimuthally dependent hyperbolic moveout approximation for small angles (or offsets) for quasi‐compressional, quasi‐shear, and converted waves in one‐dimensional multi‐layer orthorhombic media. The vertical orthorhombic axis is the same for all layers, but the azimuthal orientation of the horizontal orthorhombic axes at each layer may be different. By starting with the known equation for normal moveout velocity with respect to the surface‐offset azimuth and applying our derived relationship between the surface‐offset azimuth and phase‐velocity azimuth, we obtain the normal moveout velocity versus the phase‐velocity azimuth. As the surface offset/azimuth moveout dependence is required for analysing azimuthally dependent moveout parameters directly from time‐domain rich azimuth gathers, our phase angle/azimuth formulas are required for analysing azimuthally dependent residual moveout along the migrated local‐angle‐domain common image gathers. The angle and azimuth parameters of the local‐angle‐domain gathers represent the opening angle between the incidence and reflection slowness vectors and the azimuth of the phase velocity ψphs at the image points in the specular direction. Our derivation of the effective velocity parameters for a multi‐layer structure is based on the fact that, for a one‐dimensional model assumption, the horizontal slowness and the azimuth of the phase velocity ψphs remain constant along the entire ray (wave) path. We introduce a special set of auxiliary parameters that allow us to establish equivalent effective model parameters in a simple summation manner. We then transform this set of parameters into three widely used effective parameters: fast and slow normal moveout velocities and azimuth of the slow one. For completeness, we show that these three effective normal moveout velocity parameters can be equivalently obtained in both surface‐offset azimuth and phase‐velocity azimuth domains. 相似文献
16.
17.
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.
The conventional acoustic logging interpretation method, which is based on vertical wells that penetrate isotropic formations, is not suitable for horizontal and deviated wells penetrating anisotropic formations. This unsuitability is because during horizontal and deviated well drilling, cuttings will splash on the well wall or fall into the borehole bottom and form a thin bed of cuttings. In addition, the high velocity layers at different depths and intrinsic anisotropy may affect acoustic logging measurements. In this study, we examine how these factors affect the acoustic wave slowness measured in horizontal and deviated wells that are surrounded by an anisotropic medium using numerical simulation. We use the staggered-grid finite difference method in time domain (FDTD) combined with hybrid-PML. First, we acquire the acoustic slowness using a simulated array logging system, and then, we analyze how various factors affect acoustic slowness measurements and the differences between the effects of these factors. The factors considered are high-velocity layers, thin beds of cuttings, dipping angle, formation thickness, and anisotropy. The simulation results show that these factors affect acoustic wave slowness measurements differently. We observe that when the wavelength is much smaller than the distance between the borehole wall and high velocity layer, the true slowness of the formation could be acquired. When the wavelengths are of the same order (i.e., in the near-field scenarios), the geometrical acoustics theory is no longer applicable. Furthermore, when a thin bed of cuttings exists at the bottom of the borehole, Fermat's principle is still applicable, and true slowness can be acquired. In anisotropic formations, the measured slowness changes with increments in the dipping angle. Finally, for a measurement system with specific spacing, the slowness of a thin target layer can be acquired when the distance covered by the logging tool is sufficiently long. Based on systematical simulations with different dipping angles and anisotropy in homogenous TI media, slowness estimation charts are established to quantitatively determine the slowness at any dipping angle and for any value of the anisotropic ratio. Synthetic examples with different acoustic logging tools and different elastic parameters demonstrate that the acoustic slowness estimation method can be conveniently applied to horizontal and deviated wells in TI formations with high accuracy. 相似文献
20.
The existence of‘*-waves’has, in recent years, prompted a renewed interest in these non-geometrical arrivals, which are generated by point sources located adajcent to plane interfaces. It has led to the re-evaluation of seismic data aquisition techniques and to the question of how to use this real phenomena in enhancing existing seismic interpretation methods. This paper considers a non-geometrical SH-arrival which is generated by a point torque source unrealistically buried within a half-space. The method of solution is essentially the same as presented in an earlier paper, with the modification that the limitation placed on the distance of the source from the interface has been removed in the saddle point method used to obtain a high-frequency approximate solution. In the earlier paper, a preliminary assumption forced the saddle point, which corresponded to the *-wave arrival, to be real when it is generally complex. However, for offsets removed from the distinct ray, the imaginary part of this complex quantity is negligible. A problem which arose when comparing exact synthetic traces with those obtained using zero-order saddle point methods, was the inability to match either the amplitude or phase of the geometrical arrival in the range of offsets when the *-wave and this corresponding geometrical ray were well separated. For this range of offsets the geometrical arrival was approaching grazing incidence and another term in the saddle point expansion of the integral was necessary to rectify this error. This method is also being used to validate the results for higher order terms obtained using asymptotic ray theory. Analytical formulae are given for both the *-wave and the higher order expansion of the geometrical event, together with a comparison of synthetic seismograms using the method developed here and a numerical integration algorithm. 相似文献