共查询到20条相似文献，搜索用时 125 毫秒
1.
2.
Seismic diffracted waves carry valuable information for identifying geological discontinuities. Unfortunately, the diffraction energy is generally too weak, and standard seismic processing is biased to imaging reflection. In this paper, we present a dynamic diffraction imaging method with the aim of enhancing diffraction and increasing the signal‐to‐noise ratio. The correlation between diffraction amplitudes and their traveltimes generally exists in two forms, with one form based on the Kirchhoff integral formulation, and the other on the uniform asymptotic theory. However, the former will encounter singularities at geometrical shadow boundaries, and the latter requires the computation of a Fresnel integral. Therefore, neither of these methods is appropriate for practical applications. Noting the special form of the Fresnel integral, we propose a least‐squares fitting method based on double exponential functions to study the amplitude function of diffracted waves. The simple form of the fitting function has no singularities and can accelerate the calculation of diffraction amplitude weakening coefficients. By considering both the fitting weakening function and the polarity reversal property of the diffracted waves, we modify the conventional Kirchhoff imaging conditions and formulate a diffraction imaging formula. The mechanism of the proposed diffraction imaging procedure is based on the edge diffractor, instead of the idealized point diffractor. The polarity reversal property can eliminate the background of strong reflection and enhance the diffraction by same‐phase summation. Moreover,the fitting weakening function of diffraction amplitudes behaves like an inherent window to optimize the diffraction imaging aperture by its decaying trend. Synthetic and field data examples reveal that the proposed diffraction imaging method can meet the requirement of high‐resolution imaging, with the edge diffraction fully reinforced and the strong reflection mostly eliminated. 相似文献
3.
G.D. HUTTON 《Geophysical Prospecting》1987,35(6):681-699
The numerical modelling of seismic diffraction, e.g., at faults and other discontinuities, generally requires the use of fast approximate methods. The geophysicist responsible for the development of such numerical methods has a real need of exact solutions to certain ideal geometries to check the accuracy of his calculations. One such exact solution, which is available, is the acoustic wave solution to the perfectly reflecting wedge. The solution is three-dimensional and the source is an explosive point source. This model is ideal for seismic diffraction; the solution has the advantage of being exact, truly three-dimensional and of being in the convenient form of the temporal and spatial impulse response. More complicated sources which are extended in either space or time can, therefore, be modelled exactly by numerical integration. This paper presents some examples of the use of the perfectly reflecting wedge as a control model for an asymptotic high frequency diffraction modelling method. This control model has revealed that certain survey and wedge configurations can yield significant disagreement with, e.g., the Kirchhoff approximation. Such configurations could occur during VSP modelling when the survey lies in the near field or in the shadow zone of a high contrast fault. This control model has also been instructive in demonstrating why the high frequency, asymptotic, approximation is generally very good and has indicated a possible improvement to the Kirchhoff approximation for edge diffraction. 相似文献
4.
Numerical modeling of zero-offset laboratory data in a strong topographic environment: results for a spectral-element method and a discretized Kirchhoff integral method 
下载免费PDF全文
下载免费PDF全文 Nathalie Favretto-Cristini Anastasiya Tantsereva Paul Cristini Bjørn Ursin Dimitri Komatitsch Arkady M. Aizenberg 《地震科学（英文版）》2014,27(4):391-399
Accurate simulation of seismic wave propagation in complex geological structures is of particular interest nowadays. However conventional methods may fail to simulate realistic wavefields in environments with great and rapid structural changes, due for instance to the presence of shadow zones, diffractions and/or edge effects. Different methods, developed to improve seismic modeling, are typically tested on synthetic configurations against analytical solutions for simple canonical problems or reference methods, or via direct comparison with real data acquired in situ. Such approaches have limitations, especially if the propagation occurs in a complex environment with strong-contrast reflectors and surface irregularities, as it can be difficult to determine the method which gives the best approximation of the “real” solution, or to interpret the results obtained without an a priori knowledge of the geologic environment. An alternative approach for seismics consists in comparing the synthetic data with high-quality data collected in laboratory experiments under controlled conditions for a known configuration. In contrast with numerical experiments, laboratory data possess many of the characteristics of field data, as real waves propagate through models with no numerical approximations. We thus present a comparison of laboratory-scaled measurements of 3D zero-offset wave reflection of broadband pulses from a strong topographic environment immersed in a water tank with numerical data simulated by means of a spectral-element method and a discretized Kirchhoff integral method. The results indicate a good quantitative fit in terms of time arrivals and acceptable fit in amplitudes for all datasets. 相似文献
5.
Camilo Valencia Juan Gomez Juan Jaramillo Mario Saenz Juan Vergara 《Studia Geophysica et Geodaetica》2017,61(1):93-114
The problem of diffraction of cylindrical and plane horizontally polarized shear waves (SH waves) by a finite crack embedded in a plane bidimensional elastic full-space is revisited. Particularly, we construct an approximate solution by the addition of independent diffracted terms. In our method the derivation of the fundamental case of a semi-infinite crack obtained as a degenerate case of a generalized wedge is first considered. This result is then used as a building block to compute the diffraction of the main incident waves. The interaction between the opposite edges of the crack is later considered in terms of a series, one term at a time until a desired tolerance is reached. Moreover, we propose a procedure to determine the number of required interactions as a function of frequency. The solution derived with the superposition technique is shown to be effective at low and high frequencies and as shown by comparisons with a direct boundary element method software, highly accurate solutions are obtained after retaining just a few terms of the infinite series. 相似文献
6.
《Advances in water resources》1996,19(4):207-213
It is shown that from any solution of the linear diffusion equation, we may construct a solution of a realistic form of the Richards equation for unsaturated flow. Compared to the usual direct linearization method, our inverse approach involves a quite different sequence of transformations. This opens the possibility of exact solutions with a wider variety of continuously varying flux boundary conditions. Closed-form solutions are presented for two examples. In these, the varying water flux boundary conditions resemble (i) the passage of a peaking storm and (ii) the continuous opening of a valve preceding a steady water supply. Unlike earlier more systematic approaches to this problem, our method does not require the numerical solution of an integral equation. 相似文献
7.
本文从衍射波的物理定义出发,简化了Trorey提出的断层上Kirchhoff-Helmholtz衍射方程的解法;计算了不同深度、不同测线方向下的断层衍射波理论地震图;计算了衍射波振幅谱与相位谱,给出了利用衍射波求断层位置的公式;得出了一些新的结论,通过人工地震测深与地震勘探的实例,对断层衍射波的特性有了更明确的认识。 本文的结果表明:断层衍射波发生在地球介质剧烈变化处;衍射点两侧的衍射波走时曲线呈双曲线状;波初动清晰且半周期小;其优势频率振幅谱与反射波的相同,在衍射波与反射波走时曲线相切处附近,记录图中出现衍射波最大振幅,且波反相;视断点与真实断点一般不重合。上述特点可能为判定衍射波并确定断层位置提供判据。 相似文献
8.
A new method is presented for solving the 2D problem of diffraction of a plane wave by a wedge of arbitrary angle in a purely acoustic, constant-density medium with different constant compressional wave speeds inside and outside the wedge. The diffraction problem is formulated as integral equations, and a wavenumber–frequency representation of the scattered field is obtained. With the aid of the Cagniard–de Hoop method, exact analytical expressions in the space–time domain are obtained for the different wave constituents, i.e. geometric optical scattered waves and edge diffracted waves including head waves. These expressions can be computed to any degree of accuracy within reasonable computation times on a computer, and the semi-analytical method of solution presented thus constitutes a means of constructing reference solutions for wedge configurations. Such highly accurate reference solutions are of importance for verification of results that include diffraction phenomena modelled by general numerical approximate methods, e.g. finite differences, finite elements and spectral methods. Examples of such applications of the method of solution are given. 相似文献
9.
S.C. Mohapatra 《地球物理与天体物理流体动力学》2014,108(6):615-638
Expansion formulae associated with the interaction of oblique surface gravity waves with a floating flexible plate in the presence of a submerged horizontal flexible structure are derived using Green’s integral theorem in water of finite and infinite water depths. The associated Green’s functions are derived using the fundamental solution associated with the reduced wave equation. The integral forms of the Green’s functions and the velocity potentials are advantageous over the eigenfunction expansion method in situation when the roots of the dispersion relation coalesce. As an application of the expansion formulae, diffraction of oblique waves by a finite floating elastic plate in the presence of a submerged horizontal flexible membrane is investigated in water of finite depth. The accuracy of the numerical computation is demonstrated by analysing the convergence of the complex amplitude of the reflected waves and the energy relation. Effect of the submerged membrane on the diffraction of surface waves is studied by analysing the reflection and transmission coefficients for various parametric values. Further, the derivation of long wave equation under shallow water approximation is derived in a direct manner in the appendix. The concept and methodology can be easily extended to deal with acoustic wave interaction with flexible structures and related problems of mathematical physics and engineering. 相似文献
10.
11.
平面SV波在层状半空间中沉积谷地周围的散射 总被引：2，自引：0，他引：2
采用间接边界元方法求解了入射平面SV波在层状半空间中沉积谷地周围的散射问题.问题的解答包含自由场和散射场两部分.自由场可由直接刚度法求得,散射场由层状半空间中斜线荷载动力格林函数来模拟.文中以入射平面SV波在基岩上单一土层中沉积谷地周围散射为例研究了土层和沉积谷地周围的影响.结果表明,由于考虑了土层的动力特性,平面SV... 相似文献
12.
Acoustic plane wave scattering at a vertical fault structure represents the simplest two-dimensional model of geophysical exploration that can be investigated by analytical techniques. The exact and complete solution, in the time domain, for the scattering of the pressure field of an acoustic plane wave normally incident on a vertical fault structure is determined adapting previous results given for the frequency domain. The wave form of the pressure field of the incident plane wave is expressed by a causal time function that decays exponentially with time at every point above the fault (z<0). The zero-order term of the scattered pressure field has been computed above the fault. This zero-order term consists of an inverse Fourier transform which reduces to a closed expression forx=0, and contains an integral of a Hankel function forx#0. The high frequency part of the inverse Fourier transform forx#0 is computed employing asymptotic expressions for the Hankel function. The integral of the asymptotic expression of the Hankel function reduces to: (i) a Fresnel integral which contains a plane wave term for |x||z|; and (ii) a stationary point plane wave term plus an upper limit term for |x|=O(|z|). For the latter case the plane wave term cancels, leaving a cylindrical wave emanated from the edge of the fault. The wave front is well defined in shape, in phase and in amplitude. The amplitude of the scattered field is discontinuous atx=0, presents a jump and is well defined for |x| small and is rather smooth for |x| large. 相似文献
13.
B. I. Birger 《Izvestiya Physics of the Solid Earth》2011,47(6):541-554
Stationary solutions including wave solutions with constant amplitudes are found for nonlinear equations of thermal convection
in a layer with nonlinear rheology. The solution is based on the Fourier expansion of unknown velocities and temperatures
with only the first and first two terms retained in the velocity and temperature series, respectively. This method, which
can be regarded as the Lorenz method, yields the Lorenz equations that fairly well describe the thermal convection in a layer
with Newtonian rheology if the Rayleigh number is not very large. The obtained generalization of the Lorenz equations to the
case of an integral (having a memory) nonlinear rheology implies that only the first term is retained in the Fourier series
for the stress components, i.e., the nonlinear rheological equation is harmonically linearized. However, in the Fourier series
of temperature, it is essential to keep the second term: this term, which is independent of the horizontal coordinate, models
the thermal boundary layer that characterizes the developed convection. We constructed the bifurcation curves that describe
the stationary convection in the nonlinear integral medium simulating the rheology of the mantle, and analyzed the stability
of stationary convective flows. The Lorenz method is applied to study small-scale thermal convection in the lithosphere of
the Earth. 相似文献
14.
A. TARANTOLA 《Geophysical Prospecting》1984,32(6):998-1015
This is the first of a series of papers giving the solution of the inverse problem in seismic exploration. The acoustic approximation is used together with the assumption that the velocity field has the form . The forward problem is then linearized (thus neglecting multiple reflected waves) and the inverse problem of estimating δ is set up. Its rigorous solution can be obtained using an iterative algorithm, each step consisting of a classical Kirchhoff migration (hyperbola summation) plus a classical forward modeling step (circle summation). 相似文献
15.
Jianguo Sun 《Geophysical Prospecting》1996,44(3):351-374
Conventionally, the Fresnel zone and the geometrical spreading factor are investigated separately, because they belong to different theories of wave propagation. However, if the paraxial ray method is used for establishing the Fresnel–Kirchhoff diffraction formula for a laterally inhomogeneous multilayered medium, it can be shown that the normalized geometrical spreading factor is inversely proportional to the area of the first Fresnel zone associated with the reflection point. Therefore, if no diffracting edge cuts the first Fresnel zone, the geometrical optics approximation represents the principal part of the wavefield obtained by Fresnel–Kirchhoff diffraction theory. Otherwise, the geometrical optics approximation has to be corrected by adding edge diffractions. It is also shown that Kirchhoff-type migration and geometrical spreading factor correction both reduce the first Fresnel zone to a zone with unit area. 相似文献
16.
单程波波场延拓算子在地震偏移成像中有重要应用.单程波波场延拓算子按其实现方式可分为Kirchhoff积分、空间隐式有限差分和Fourier变换方法,他们代表了算子的不同表示方法,当截断使用这些方法时会得到不同的精度.象征表示对这些方法的导出和精度分析有重要作用.算子作用于正弦波函数所得函数称为算子的象征.算子的象征是褶积算子Fourier变换的推广.Fourier变换方法则直接用象征函数的可分表示求出.空间隐式有限差分则可以用象征函数的Padè近似或部分分式导出.单程波算子在深度域的积分称为单程波算子积分解.本文推导了单程波算子积分解的象征表达式,给出了算子象征的代数运算的头几阶表达式,这些表达式还未在前人文献中发现.Kirchhoff积分所需格林函数可以通过象征函数和鞍点法导出.基于积分解的象征表达式给出了非对称走时公式,对改善Kirchhoff积分的聚焦性能有重要意义. 相似文献
17.
18.
Morteza Eskandari-Ghadi Azizollah Ardeshir-Behrestaghi 《Soil Dynamics and Earthquake Engineering》2010
This paper is concerned with the investigation of the vertical vibration of a rigid circular disc buried at an arbitrary depth in a transversely isotropic half space in such a way the axis of material symmetry of the half space is normal to the surface of it and parallel to the vibration direction. By using the Hankel integral transforms, the mixed boundary-value problem is transformed to a pair of integral equations called dual integral equations, which generally can be reduced to a Fredholm integral equation of the second kind. With the aid of complex variable or contour integration, the governing integral equation is numerically solved in the general dynamic case. Two degenerated cases (i) the disc is buried in a transversely isotropic full space, and (ii) rigid circular disc is attached on the surface of the half space are discussed. The reduced static case of the dual integral equations is solved analytically and the vertical displacement, the contact pressure and the static impedance/compliance function are explicitly found. It is shown that the vertical pressure and the compliance function reduced for isotropic half space are identical to the previous solutions reported in the literature. The dynamic contact pressure under the disc and the impedance function are numerically evaluated in general dynamic case and graphically shown that the singularity exists in the contact pressure at the edge of the disc is the same as the static case. In addition, the impedance functions evaluated here for the isotropic domain are collapsed on the solution given by Luco and Mita. To show the effect of different material anisotropy, the numerical evaluations are given for some different transversely isotropic materials and compared. 相似文献
19.
Vladimir Priklonsky Anvar Shukurov Dmitry Sokoloff Andrew Soward 《地球物理与天体物理流体动力学》2013,107(1-2):97-114
Abstract We reconsider thin-disc global asymptotics for kinematic, axisymmetric mean-field dynamos with vacuum boundary conditions. Non-local terms arising from a small but finite radial field component at the disc surface are consistently taken into account for quadrupole modes. As in earlier approaches, the solution splits into a local part describing the field distribution along the vertical direction and a radial part describing the radial (global) variation of the eigenfunction. However, the radial part of the eigenfunction is now governed by an integro-differential equation whose kernel has a weak (logarithmic) singularity. The integral term arises from non-local interactions of magnetic fields at different radii through vacuum outside the disc. The non-local interaction can have a stronger effect on the solution than the local radial diffusion in a thin disc, however the effect of the integral term is still qualitatively similar to magnetic diffusion. 相似文献
20.
ABSTRACT Computation of the wavefield due to reflection from an irregular surface is carried out for subsurfaces with large radii of curvature. The Kirchhoff approximation is proved to be sufficiently accurate provided that the acoustic wavelength is sufficiently small with respect to the asperities of the rough surface. For cases where the irregular surface does not fulfil this condition, a series solution is proposed. The first term of this series appears to be the result obtained by conventional Kirchhoff approximation. The series, initially developed in the space–wavenumber domain by Meecham, is transformed into the space–time domain, and the general expression for the series is obtained by calculation of the normal derivative of the field function. The series solution, restricted to the first two terms, is illustrated by application to three synthetic examples. Applications show that the series approximation obtained by the Kirchhoff method contributes significantly to the modelling of narrow, steep and deep structures and consequently it appears that the second term in the series cannot be ignored in the computation of the wavefields arising from a rough surface. 相似文献