共查询到20条相似文献，搜索用时 406 毫秒
1.
单程波近似实际上是一种多次前向散射和单次后向散射近似.利用单程波近似来描述波传播可以极大地节省地震数值模拟的计算时间和内存,实现地震波长距离传播模拟和三维地震模拟快速计算.本文基于单程波近似和波动积分方程的分离变量逼近,从广义Lippmann-Schwinger波动积分方程推导出耦合反射/透射系数的单程波传播算子.该算子由两部分构成：分离变量Fourier单程波传播算子和薄板间的反射/透射系数表达.前者将常规的Fourier分裂步单程波传播算子(SSF)推广适应横向强速度变化介质和大角度传播波场.后者是利用垂直波数来表示反射/透射系数,自然耦合到波场传播的计算过程中,其为地质界面倾角的隐式表达,精确描述振幅随入射角的变化,能适应任意复杂的模型.通过两个数值算例和一个实际地质模型的计算,本文将该方法和边界元法进行了比较,结果表明：在算例给出的介质横向速度变化情况下,本文提出的方法在相位和振幅方面与全波数值方法基本吻合. 相似文献
2.
Wave Propagation, Scattering and Imaging Using Dual-domain One-way and One-return Propagators 总被引：4，自引：0，他引：4
R.-S. Wu 《Pure and Applied Geophysics》2003,160(3-4):509-539
— Dual-domain one-way propagators implement wave propagation in heterogeneous media in mixed domains (space-wavenumber domains). One-way propagators neglect wave reverberations between heterogeneities but correctly handle the forward multiple-scattering including focusing/defocusing, diffraction, refraction and interference of waves. The algorithm shuttles between space-domain and wavenumber-domain using FFT, and the operations in the two domains are self-adaptive to the complexity of the media. The method makes the best use of the operations in each domain, resulting in efficient and accurate propagators. Due to recent progress, new versions of dual-domain methods overcame some limitations of the classical dual-domain methods (phase-screen or split-step Fourier methods) and can propagate large-angle waves quite accurately in media with strong velocity contrasts. These methods can deliver superior image quality (high resolution/high fidelity) for complex subsurface structures. One-way and one-return (De Wolf approximation) propagators can be also applied to wave-field modeling and simulations for some geophysical problems. In the article, a historical review and theoretical analysis of the Born, Rytov, and De Wolf approximations are given. A review on classical phase-screen or split-step Fourier methods is also given, followed by a summary and analysis of the new dual-domain propagators. The applications of the new propagators to seismic imaging and modeling are reviewed with several examples. For seismic imaging, the advantages and limitations of the traditional Kirchhoff migration and time-space domain finite-difference migration, when applied to 3-D complicated structures, are first analyzed. Then the special features, and applications of the new dual-domain methods are presented. Three versions of GSP (generalized screen propagators), the hybrid pseudo-screen, the wide-angle Padé-screen, and the higher-order generalized screen propagators are discussed. Recent progress also makes it possible to use the dual-domain propagators for modeling elastic reflections for complex structures and long-range propagations of crustal guided waves. Examples of 2-D and 3-D imaging and modeling using GSP methods are given. 相似文献
3.
Sensitivity of seismic waves to structure 总被引：2，自引：0，他引：2
Luděk Klime? 《Studia Geophysica et Geodaetica》2012,56(2):483-520
We study how the perturbations of a generally heterogeneous isotropic or anisotropic structure manifest themselves in the wavefield, and which perturbations can be detected within a limited aperture and a limited frequency band. A short-duration broad-band incident wavefield with a smooth frequency spectrum is considered. In-finitesimally small perturbations of elastic moduli and density are decomposed into Gabor functions. The wavefield scattered by the perturbations is then composed of waves scattered by the individual Gabor functions. The scattered waves are estimated using the first-order Born approximation with the paraxial ray approximation. For each incident wave, each Gabor function generates at most 5 scattered waves, propagating in specific directions and having specific polarisations. A Gabor function corresponding to a low wavenumber may generate a single broad-band unconverted wave scattered in forward or narrow-angle directions. A Gabor function corresponding to a high wavenumber usually generates 0 to 5 narrow-band Gaussian packets scattered in wide angles, but may also occasionally generate a narrow-band P to S or S to P converted Gaussian packet scattered in a forward direction, or a broad-band S to P (and even S to S in a strongly anisotropic background) converted wave scattered in wide angles. In this paper, we concentrate on the Gaussian packets caused by narrow-band scattering. For a particular source, each Gaussian packet scattered by a Gabor function at a given spatial location is sensitive to just a single linear combination of 22 values of the elastic moduli and density corresponding to the Gabor function. This information about the Gabor function is lost if the scattered wave does not fall into the aperture covered by the receivers and into the legible frequency band. 相似文献
4.
利用地震波正向传播方程对属于波形线性反演问题近似求解方法的地震数据偏移成像进行重新推导,得到了适合散射地震数据的散射偏移成像方法和适合反射地震数据的反射偏移成像方法.以地震波传播的散射理论为出发点,首先根据描述一次散射波正向传播的线性方程研究建立散射地震数据的偏移成像方法理论;利用高频近似对产生散射波场的地下速度扰动函数的空间变化进行近似,推导出地下反射率函数,再由散射波传播方程推导出基于反射率函数的反射波传播方程,然后根据描述一次反射波正向传播的线性方程研究建立反射地震数据的偏移成像方法理论.本文指出和修正了Claerbout偏移成像方法中的不足,提出的地震数据偏移成像方法是对当前偏移成像方法理论的完善,使反射地震数据偏移成像具有了更坚实的数学物理理论基础,得到的偏移成像结果相位正确、位置准确、分辨率提高. 相似文献
5.
6.
傅里叶有限差分法三维波动方程正演模拟 总被引：4，自引：6，他引：4
傅里叶有限差分(FFD)法兼有相位屏法和隐式有限差分法二者的优势，能够处理复杂地质构造中的波传播问题，但在三维情形下，算子的双向分裂会引起明显的方位各向异性误差.本文用Fourier变换计算双向分裂过程中的高阶交叉项，消除了方位各向异性误差.该方法充分利用了FFD法在双域实现的算法结构，明显减少了由于引入误差校正所带来的计算量.将该方法应用于修改后的三维French模型的地震正演问题，并将得到的叠后记录、单炮记录同全波有限差分法的模拟结果进行对比，结果证实了该方法对一次反射波具有较高的模拟精度，在内存需求和计算效率方面则具有更大的优势. 相似文献
7.
For wave propagation simulation in piecewise heterogeneous media, Gaussian-elimination-based full-waveform solutions to the generalized Lippmann–Schwinger integral equation (GLSIE) are highly accurate, but involved with extremely time-consuming computations because of the very large size of the resulting boundary–volume integral equation matrix to be inverted. Several flexible approximations to the GLSIE are scaled in an iterative way to adapt numerical solutions to the smoothness of heterogeneous media in terms of incident wavelengths, with a great saving of computing time and memory. Among various typical iterative schemes to the GLSIE matrix, the generalized minimal residual method (GMRES) is an efficient approach to reduce the computational intensity to some degree. The most efficient approximation can be obtained using a Born series, as an alternative iterative solution, to both the boundary-scattering and volume-scattering waves, leading to the Born-series approximation (BSA) scheme and the improved Born-series approximation (IBSA) scheme. These iteration schemes are validated by dimensionless frequency responses to a heterogeneous semicircular alluvial valley, and then applied to a heterogeneous multilayered model by calculating synthetic seismograms to evaluate approximation accuracies. Numerical experiments, compared with the full-waveform numerical solution, indicate that the convergence rates of these methods decrease gradually with increasing velocity perturbations. The comparison also shows that the BSA scheme has a faster convergence than the GMRES method for velocity perturbations less than 10 percent, but converges slowly and even hardly achieves convergence for velocity perturbations greater than 15 percent. The IBSA scheme gives a superior performance over the other methods, with the least iterations to achieve the necessary convergence. 相似文献
8.
9.
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. 相似文献
10.
Morten Jakobsen 《Studia Geophysica et Geodaetica》2012,56(1):1-20
Forward seismic modelling in the acoustic approximation, for variable velocity but constant density, is dealt with. The wave
equation and the boundary conditions are represented by a volume integral equation of the Lippmann-Schwinger (LS) or Fredholm
type. A T-matrix (or transition operator) approach from quantum mechanical potential scattering theory is used to derive a
family of linear and nonlinear approximations (cluster expansions), as well as an exact numerical solution of the LS equation.
For models of 4D anomalies involving small or moderate contrasts, the Born approximation gives identical numerical results
as the first-order t-matrix approximation, but the predictions of an exact T-matrix solution can be quite different (depending
on spatial extention of the perturbations). For models of fluid-saturated cavities involving large or huge contrasts, the
first-order t-matrix approximation is much more accurate than the Born approximation, although it does not lead to significantly
more time-consuming computations. If the spatial extention of the perturbations is not too large, it is practical to use the
exact T-matrix solution which allows for arbitrary contrasts and includes all the effects of multiple scattering. 相似文献
11.
Broadband constant-coefficient propagators 总被引：4，自引：1，他引：4
Li-Yun Fu 《Geophysical Prospecting》2005,53(3):299-310
The phase error between the real phase shift and the Gazdag background phase shift, due to lateral velocity variations about a reference velocity, can be decomposed into axial and paraxial phase errors. The axial phase error depends only on velocity perturbations and hence can be completely removed by the split‐step Fourier method. The paraxial phase error is a cross function of velocity perturbations and propagation angles. The cross function can be approximated with various differential operators by allowing the coefficients to vary with velocity perturbations and propagation angles. These variable‐coefficient operators require finite‐difference numerical implementation. Broadband constant‐coefficient operators may provide an efficient alternative that approximates the cross function within the split‐step framework and allows implementation using Fourier transforms alone. The resulting migration accuracy depends on the localization of the constant‐coefficient operators. A simple broadband constant‐coefficient operator has been designed and is tested with the SEG/EAEG salt model. Compared with the split‐step Fourier method that applies to either weak‐contrast media or at small propagation angles, this operator improves wavefield extrapolation for large to strong lateral heterogeneities, except within the weak‐contrast region. Incorporating the split‐step Fourier operator into a hybrid implementation can eliminate the poor performance of the broadband constant‐coefficient operator in the weak‐contrast region. This study may indicate a direction of improving the split‐step Fourier method, with little loss of efficiency, while allowing it to remain faster than more precise methods such as the Fourier finite‐difference method. 相似文献
12.
13.
研究了三维弱各向异性近似下,利用伪P波(伪纵波)模拟弹性波场P分量在倾斜对称轴的横向各向同性(TTI)介质中的传播过程,并对比了分别基于弹性Hooke定律、弹性波投影和运动学色散方程所建立的三种二阶差分伪P波方程的正演特点.目前这些伪P波方程数值计算主要采用规则网格差分,但是规则网格在TTI模拟中有低效率、低精度以及不稳定的缺点.为了提高计算的精度,本文构建出相应方程的交错网格有限差分格式.通过对比伪P波方程在三维TTI介质中不同的数值模拟的表达形式,本文认为基于色散方程所建立的伪P波方程在模拟弹性波中P波传播的过程中具有最小的噪声.本文分析不同的各向同性对称轴空间角度的频散特征,并引入适当的横波速度维持计算的稳定.二维模型算例表明,本文提出的交错网格正演算法可以得到稳定光滑的伪P波正演波场.使用本文交错网格算法对二维BP TTI模型的逆时偏移也具有较稳定的偏移结果. 相似文献
14.
—We investigate a new nonlinear inversion method for low frequencies to determine the bulk and shear modulus as well as the material density and the location of subsurface inhomogeneities. The solution is a direct exact nonlinear inversion of single scattered waves containing near- and far-field terms for incident P and scattered P and S waves, allowing for inversion of parameters in the vicinity and at distance from the sources and receivers. Because the approach is based on single scattering theory, the range of application includes single strong scattering anomalies of various sizes like magma chambers, gas- or fluid-filled cavities, or buried near-surface obstacles. The replacement of the material properties by a new set of parameters, referred to as scattering factors, allows the inverse problem to be solved analytically. The nonlinear nature of the scattering problem is investigated and implications for the inversion process are discussed. The deviations in the elastic parameters as a function of the scattering factors show a strong asymmetry about zero, and therefore linearized approximations will perform differently, depending on the sign of the perturbation. Based on the low frequency (Rayleigh) approximation, we introduce and evaluate a pair of approximations (Mie) derived by numerical and analytical integration of the Rayleigh approximation. Both approximations are based on the underlying principle of subdividing the inhomogeneities into a number of small noninteracting parts and subsequent integration over the total volume, thus increasing the Rayleigh limit and producing better resolution of the parameter estimates during the inversion. The two Mie approximations, when evaluated as a function of scattering angle and distance, produce similar results in the mid- and far-field of the inhomogeneity and reveal better resolution than the Rayleigh approximation. For three anomalies of ± 50% in bulk modulus, shear modulus, and density, the relative error between the exact solution and the two Mie approximations remains below 10%, 20%, and 30%, respectively, for values of k p R < 3.0, where R is the radius of the heterogeneity. However, smaller errors for individual cases are found for values up to k p R≈ 4.5. The performance of the inversion based on the analytically and the numerically integrated Mie approximation is tested for single parameter perturbations, revealing reliable and stable inversion results for the bulk and the shear modulus, reasonable results for the density, and crosstalk between the shear modulus and the density. The results show well-defined locations of the anomalies and slight deviations in the estimates of their magnitudes, which can be explained by amplitude and phase deviations between the analytical solution used for forward modeling and the approximations used for the inversion. The analytical Mie approximation provides a fast means to estimate elastic parameters compared to the more time consuming numerically integrated approximation, while the latter can be applied to more arbitrarily shaped inhomogeneities. 相似文献
15.
16.
本文以基于改进BISQ模型的二维双相各向同性介质一阶速度-应力方程为基础,推导出了曲线坐标系下对应的方程,然后采用低频散、低耗散的同位网格MacCormack有限差分法来离散方程,并采用紧致的单边MacCormack差分格式结合牵引力镜像法来施加自由地表边界条件,实现了地震波场数值模拟.曲线网格有限差分法采用贴体网格来描述自由表面,地表的网格线紧贴地形,避免了台阶近似造成的数值散射.数值模拟结果表明,在双相介质起伏自由地表和分界面处,各类波型复杂的反射透射规律可以清晰展现,曲线网格有限差分法可以精确地解决地震波在含起伏地表的双相各向同性介质中的传播问题. 相似文献
17.
An efficient approximate scheme is presented for wave-propagation simulation in piecewise heterogeneous media by applying the Born-series approximation to volume-scattering waves. The numerical scheme is tested for dimensionless frequency responses to a heterogeneous alluvial valley where the velocity is perturbed randomly in the range of 5 %–25 %, compared with the full-waveform numerical solution. Then, the scheme is extended to a heterogeneous multilayered model by calculating synthetic seismograms to evaluate approximation accuracies. Numerical experiments indicate that the convergence rate of this method decreases gradually with increasing velocity perturbations. The method has a fast convergence for velocity perturbations less than 15 %. However, the convergence becomes slow drastically when the velocity perturbation increases to 20 %. The method can hardly converge for the velocity perturbation up to 25 %. 相似文献
18.
The perturbation method in elastic wave scattering 总被引：3，自引：0，他引：3
Ru-Shan Wu 《Pure and Applied Geophysics》1989,131(4):605-637
Methods of theoretical study in seismic wave scattering are reviewed with the emphasis on the perturbation method. Detailed analysis for weak scattering using Born approximation is given. For elastic random media, the mean square amplitudes of scattered waves are derived using a new approach by working directly in the spectrum domain. The conditions for the scalar wave approximation are obtained. The problem of sensitivity of fore- and backscattering to heterogeneities with different scales and properties (velocity or impedance) is discussed. 相似文献
19.
模拟弹性波在孔隙介质中传播，需要稳定有效的吸收边界来消除或尽可能的减小由人工边界引起的虚假反射. 本文在前人工作基础上，首次建立了弹性孔隙介质情况下完全匹配层吸收边界的高阶速度-应力交错网格有限差分算法，并详细讨论了完全匹配层的构建及其有限差分算法实现. 首先，本文通过均匀孔隙模型的数值解与解析解的对比，验证所提出的数值方法的正确性；然后，本文考察了完全匹配层对不同入射角度入射波和自由表面上的瑞利波的吸收性能，将完全匹配层与廖氏和阻尼吸收边界进行了对比，研究了这三种吸收边界在不同吸收厚度情况下对弹性波吸收能力. 数值结果表明，在孔隙介质中，完全匹配层作为吸收边界能十分有效地吸收衰减外行波，无论对体波还是面波，是一种高效边界吸收算法. 相似文献
20.
This paper is concerned with numerical tests of several rock physical relationships. The focus is on effective velocities and scattering attenuation in 3D fractured media. We apply the so‐called rotated staggered finite‐difference grid (RSG) technique for numerical experiments. Using this modified grid, it is possible to simulate the propagation of elastic waves in a 3D medium containing cracks, pores or free surfaces without applying explicit boundary conditions and without averaging the elastic moduli. We simulate the propagation of plane waves through a set of randomly cracked 3D media. In these numerical experiments we vary the number and the distribution of cracks. The synthetic results are compared with several (most popular) theories predicting the effective elastic properties of fractured materials. We find that, for randomly distributed and randomly orientated non‐intersecting thin penny‐shaped dry cracks, the numerical simulations of P‐ and S‐wave velocities are in good agreement with the predictions of the self‐consistent approximation. We observe similar results for fluid‐filled cracks. The standard Gassmann equation cannot be applied to our 3D fractured media, although we have very low porosity in our models. This is explained by the absence of a connected porosity. There is only a slight difference in effective velocities between the cases of intersecting and non‐intersecting cracks. This can be clearly demonstrated up to a crack density that is close to the connectivity percolation threshold. For crack densities beyond this threshold, we observe that the differential effective‐medium (DEM) theory gives the best fit with numerical results for intersecting cracks. Additionally, it is shown that the scattering attenuation coefficient (of the mean field) predicted by the classical Hudson approach is in excellent agreement with our numerical results. 相似文献