共查询到20条相似文献，搜索用时 161 毫秒
1.
The dispersion relation for Love waves in a layer on a half-space is modified by introducing the wave number and its square instead of the phase velocity. The implicit function theorem is then used to derive the analytical formulae for the group velocity and for the phase- and group-velocity partial derivatives with respect to the parameters of the medium. The formulae are compared with those obtained by Novotný (1971) where the traditional formulation of the dispersion relation was used. 相似文献
2.
The dispersion relation for Rayleigh waves in a layer on a half-space is modified by introducing the quadratic wave number instead of the phase velocity. The implicit function theorem is then used to derive analytical formulae for the group velocity and for the phase- and group velocity partial derivatives with respect to the parameters of the medium. As two examples, the method is applied to the interpretation of dispersion curves for short-period Rayleigh waves observed in a sedimentary basin of the West Carpathians, and in the Moravo-Silesian region, Czech Republic. 相似文献
3.
Tariq Alkhalifah 《Geophysical Prospecting》2011,59(4):627-634
The wavefield dependence on a virtual shift in the source location can provide information helpful in velocity estimation and interpolation. However, the second‐order partial differential equation (PDE) that relates changes in the wavefield form (or shape) to lateral perturbations in the source location depends explicitly on lateral derivatives of the velocity field. For velocity models that include lateral velocity discontinuities this is problematic as such derivatives in their classical definition do not exist. As a result, I derive perturbation partial differential wave equations that are independent of direct velocity derivatives and thus, provide possibilities for wavefield shape extrapolation in complex media. These PDEs have the same structure as the wave equation with a source function that depends on the background (original source) wavefield. The solutions of the perturbation equations provide the coefficients of a Taylor's series type expansion for the wavefield. The new formulas introduce changes to the background wavefield only in the presence of lateral velocity variation or in general terms velocity variations in the perturbation direction. The accuracy of the representation, as demonstrated on the Marmousi model, is generally good. 相似文献
4.
Ｉｎｔｒｏｄｕｃｔｉｏｎ３－Ｄｓｅｉｓｍｉｃｔｏｍｏｇｒａｐｈｙｈａｓｂｅｅｎａｐｐｌｉｅｄｔｏｖａｒｉｏｕｓｇｅｏｐｈｙｓｉｃａｌｐｒｏｂｌｅｍｓ．ＡｋｉａｎｄＬｅｅ（１９７６）ａｎｄＨａｗｌｅｙｅｔａｌ．（１９８１）ｉｎｖｅｒｔｅｄ３－Ｄｍｏｄｅ... 相似文献
5.
Pantelis Soupios Irfan Akca Petros Mpogiatzis Ahmet T. Basokur Constantinos Papazachos 《Journal of Applied Geophysics》2011,75(3):479-489
Almost all earth sciences inverse problems are nonlinear and involve a large number of unknown parameters, making the application of analytical inversion methods quite restrictive. In practice, most analytical methods are local in nature and rely on a linearized form of the problem equations, adopting an iterative procedure which typically employs partial derivatives in order to optimize the starting (initial) model by minimizing a misfit (penalty) function. Unfortunately, especially for highly non-linear cases, the final model strongly depends on the initial model, hence it is prone to solution-entrapment in local minima of the misfit function, while the derivative calculation is often computationally inefficient and creates instabilities when numerical approximations are used. An alternative is to employ global techniques which do not rely on partial derivatives, are independent of the misfit form and are computationally robust. Such methods employ pseudo-randomly generated models (sampling an appropriately selected section of the model space) which are assessed in terms of their data-fit. A typical example is the class of methods known as genetic algorithms (GA), which achieves the aforementioned approximation through model representation and manipulations, and has attracted the attention of the earth sciences community during the last decade, with several applications already presented for several geophysical problems.In this paper, we examine the efficiency of the combination of the typical regularized least-squares and genetic methods for a typical seismic tomography problem. The proposed approach combines a local (LOM) and a global (GOM) optimization method, in an attempt to overcome the limitations of each individual approach, such as local minima and slow convergence, respectively. The potential of both optimization methods is tested and compared, both independently and jointly, using the several test models and synthetic refraction travel-time date sets that employ the same experimental geometry, wavelength and geometrical characteristics of the model anomalies. Moreover, real data from a crosswell tomographic project for the subsurface mapping of an ancient wall foundation are used for testing the efficiency of the proposed algorithm. The results show that the combined use of both methods can exploit the benefits of each approach, leading to improved final models and producing realistic velocity models, without significantly increasing the required computation time. 相似文献
6.
Summary The partial derivatives of the dispersion curves have been calculated by means of a method founded on the theorem of implicit functions, which was described in[1]. The partial derivatives of the dispersion curves of Love waves are studied in detail for a model of a single-layered medium. Analytical formulae for the partial derivatives are used for deriving the limits of these partial derivatives when the period approaches zero or the critical period. 相似文献
7.
From the partial differential equations of hydrodynamics governing the movements in the Earth's mantle of a Newtonian fluid with a pressure- and temperature-dependent viscosity, considering the bilateral symmetry of velocity and temperature distributions at the mid-plane of the plume, an analytical solution of the governing equations near the mid-plane of the plume was found by the method of asymptotic analysis. The vertical distribution of the upward velocity, viscosity and temperature at the mid-plane, and the temperature excess at the centre of the plume above the ambient mantle temperature were then calculated for two sets of Newtonian rheological parameters. The results obtained show that the temperature at the mid-plane and the temperature excess are nearly independent of the rheological parameters. The upward velocity at the mid-plane, however, is strongly dependent on the rheological parameters. 相似文献
8.
9.
J. GAZDAG 《Geophysical Prospecting》1980,28(1):60-70
A stacked seismic section represents a wave-field recorded at regularly spaced points on the surface. The seismic migration process transforms this recorded data into a reflectivity display. In recent years, Jon F. Claerbout and his co-workers developed migration techniques based on the numerical approximation of the wave equation by finite difference methods. This paper describes an alternative method, termed ASD (for Accurate Space Derivative), and its application to the wave equation migration problem. In this approach to the numerical solution of partial differential equations, partial derivatives are computed by finite Fourier transform methods. This migration method can accommodate media with vertical as well as horizontal velocity variations. 相似文献
10.
为了对多条二维剖面资料进行统一分析和解释,本文考虑了各剖面交点处的约束条件,通过建立统一的偏导矩阵,将各剖面资料联立求解.利用该方法对中国地震局地球物理勘探中心在首都圈内6条人工地震宽角反射/折射剖面资料进行了重新处理和解释,得到该地区的地壳三维速度结构和构造特征.在此基础上,结合其他地球物理研究成果,确定了壳内存在的4处低速层（体）的空间范围及断裂的走向. 相似文献
11.
本文以基于改进BISQ模型的二维双相各向同性介质一阶速度-应力方程为基础,推导出了曲线坐标系下对应的方程,然后采用低频散、低耗散的同位网格MacCormack有限差分法来离散方程,并采用紧致的单边MacCormack差分格式结合牵引力镜像法来施加自由地表边界条件,实现了地震波场数值模拟.曲线网格有限差分法采用贴体网格来描述自由表面,地表的网格线紧贴地形,避免了台阶近似造成的数值散射.数值模拟结果表明,在双相介质起伏自由地表和分界面处,各类波型复杂的反射透射规律可以清晰展现,曲线网格有限差分法可以精确地解决地震波在含起伏地表的双相各向同性介质中的传播问题. 相似文献
12.
In this paper, we deduced the corresponding first-order velocity–stress equation for curvilinear coordinates from the first-order velocity–stress equation based on the modified Biot/squirt model for a two-dimensional two-phase medium. The equations are then numerically solved by an optimized high-order non-staggered finite difference scheme, that is, the dispersion relation preserving/optimization MacCormack scheme. To implement undulating free-surface topography, we derive an analytical relationship between the derivatives of the particle velocity components and use the compact finite-difference scheme plus a traction-image method. In the undulating free surface and the undulating subsurface interface of two-phase medium, the complex reflected wave and transmitted wave can be clearly recognized in the numerical simulation results. The simulation results show that the curvilinear-grid finite-difference method, which uses a body-conforming grid to describe the undulating surface, can accurately reduce the numerical scattering effect of seismic wave propagation caused by the use of ladder-shaped grid to fit the surfaces when undulating topography is present in a two-phase isotropic medium. 相似文献
13.
Summary This paper is a continuation of[1]. It is mainly devoted to problems connected with the application of the method of determination of geometrical spreading in laterally inhomogeneous media with curved interfaces based on the solution of eight (in a three-dimensional medium) or two (in a two-dimensional medium) linear ordinary differential equations of the first order. The method of determination of the partial derivatives of velocity with respect to the special coordinates, connected with the ray under investigation, and the methods of determination of the initial values for the system of differential equations at the source and at the interfaces are proposed. 相似文献
14.
通过对球坐标梯度速度模型的角距离和走时公式的推导,获取了射线参数驱动的可解析计算的全球多震相时距及其对于模型,震中和地表速度的导数表达式.与TAUP软件进行对比发现除了散射震相外,所有震相的时距曲线误差均可忽略,该方法考虑了高程,在性能和应用方面也可以完全取代TAUP相应的模块.在此基础上升级了耦合一维模型的定位系统,并应用到龙门山地区的近震重新定位,通过对远震的局部化处理,增加了近震层析的射线覆盖范围和密度,使下地壳速度模型在中上地壳模型的约束下有所呈现.通过近震层析处理发现松潘地区浅部存在大量的高低速相间分布,四川盆地深部为扬子板块高速物质控制,浅部则以低速为主.龙门山断裂带浅部呈现高速为主的特征,而中地壳广泛存在低速物质的分布,龙门山断裂控制的彭灌杂岩体存在根部,收束在20 km深度中地壳滑脱带上. 相似文献
15.
中国大陆中西部乃至全球造山带普遍具有复杂地壳结构.随着矿产资源勘探和深部探测研究的深入,探测造山带及盆山耦合区下方地壳精细结构正逐渐成为当前面临的巨大挑战.人工源深地震测深方法正越来越清晰地揭示出不同构造域地壳速度结构的基本特征,然而传统的层状结构模型参数化方法难以准确描述复杂地质模型,通常情况下多忽略速度结构的精细间断面且采用层边界平滑处理,难以满足地壳精细结构成像的发展要求.针对上述困难,本文采用最近发展的块状结构建模方案构建三维复杂地壳模型,基于逐段迭代射线追踪正演走时计算方法,推导了走时对三角形界面深度以及网格速度的偏导数,开展了非线性共轭梯度走时反演方法研究.发展了利用直达波和反射波等多震相走时数据对界面深度和网格速度的多参数联合反演方法,并引人不同种类震相数据的权系数和不同类型参数偏导数归一化的方法.数值算例表明,基于块状结构的非线性共轭梯度走时反演方法适用于复杂地壳结构模型,在利用人工源走时数据反演复杂地壳精细结构领域具有良好的应用前景. 相似文献
16.
Emil Blias 《Geophysical Prospecting》2009,57(3):323-341
Lateral velocity changes (velocity anomalies) in the overburden may cause significant oscillations in normal moveout velocities. Explicit analytical moveout formulas are presented and provide a direct explanation of these lateral fluctuations and other phenomena for a subsurface with gentle deep structures and shallow overburden anomalies. The analytical conditions for this have been derived for a depth-velocity model with gentle structures with dips not exceeding 12°. The influence of lateral interval velocity changes and curvilinear overburden velocity boundaries can be estimated and analysed using these formulas. An analytical approach to normal moveout velocity analysis in a laterally inhomogeneous medium provides an understanding of the connection between lateral interval velocity changes and normal moveout velocities. In the presence of uncorrected shallow velocity anomalies, the difference between root-mean-square and stacking velocity can be arbitrarily large to the extent of reversing the normal moveout function around normal incidence traveltimes. The main reason for anomalous stacking velocity behaviour is non-linear lateral variations in the shallow overburden interval velocities or the velocity boundaries.
A special technique has been developed to determine and remove shallow velocity anomaly effects. This technique includes automatic continuous velocity picking, an inversion method for the determination of shallow velocity anomalies, improving the depth-velocity model by an optimization approach to traveltime inversion (layered reflection tomography) and shallow velocity anomaly replacement. Model and field data examples are used to illustrate this technique. 相似文献
A special technique has been developed to determine and remove shallow velocity anomaly effects. This technique includes automatic continuous velocity picking, an inversion method for the determination of shallow velocity anomalies, improving the depth-velocity model by an optimization approach to traveltime inversion (layered reflection tomography) and shallow velocity anomaly replacement. Model and field data examples are used to illustrate this technique. 相似文献
17.
A simple and exact numerical scheme for long-term simulations of 3D potential fully nonlinear periodic gravity waves is suggested. The scheme is based on the surface-following nonorthogonal curvilinear coordinate system. Velocity potential is represented as a sum of analytical and nonlinear components. The Poisson equation for the nonlinear component of velocity potential is solved iteratively. Fourier transform method, the second-order accuracy approximation of vertical derivatives on a stretched vertical grid and the fourth-order Runge–Kutta time stepping are used. The scheme is validated by simulation of steep Stokes waves. A one-processor version of the model for PC allows us to simulate evolution of a wave field with thousands degrees of freedom for hundreds of wave periods. The scheme is designed for investigation of nonlinear 2D surface waves, generation of extreme waves, and direct calculations of nonlinear interactions. 相似文献
18.
Analytic expressions for the velocity sensitivity to the elastic moduli for the most general anisotropic media 总被引：1，自引：0，他引：1
For non‐linear kinematic inversion of elastic anisotropy parameters and related investigations of the sensitivity of seismic data, the derivatives of the wavespeed (phase velocity and group velocity) with respect to the individual elastic moduli are required. This paper presents two analytic methods, called the eigenvalue and eigenvector methods, to compute the derivatives of the wavespeeds for wave propagation in a general anisotropic medium, which may be defined by up to 21 density‐normalized elastic moduli. The first method employs a simple and compact form of the eigenvalue (phase velocity) and a general form of the group velocity, and directly yields general expressions of the derivatives for the three wave modes (qP, qS1, qS2). The second method applies simple eigenvector solutions of the three wave modes and leads to other general forms of the derivatives. These analytic formulae show that the derivatives are, in general, functions of the 21 elastic moduli as well as the wave propagation direction, and they reflect the sensitivity of the wavespeeds to the individual elastic moduli. Meanwhile, we give results of numerical investigations with some examples for particular simplified forms of anisotropy. They show that the eigenvalue method is suitable for the qP‐, qS1‐ and qS2‐wave computations and mitigates the singularity problem for the two quasi‐shear waves. The eigenvector method is preferable to the eigenvalue method for the group velocity and the derivative of the phase velocity because it involves simpler expressions and independent computations, but for the derivative of the group velocity the derivative of the eigenvector is required. Both methods tackle the singularity problem and are applicable to any degree of seismic anisotropy for all three wave modes. 相似文献
19.
本文介绍了多震相的层析成像的思路和算法，通过穿透和反射走时可以同时作出２维和３维慢度（速度的倒数）重建。我们分析了在穿透和反射数据中确定速度和深度的不确定性，并认识到深度扰动对反射走时异常比慢度扰动更敏感。由不同波类型所提供的对速度和深度的约束，这个算法实际上减少了在一般反射层析成像在速度和深度之间的不确定性，并且也避免了在穿透层析成像中的不确定问题。线性化反演是通过从反射界面深度由分离速度参数迭代进行的。使用一个快速的２－Ｄ和３－Ｄ射线跟踪算法来计算穿透和反射走时和对幔度及反射界面深度的偏导数。深度和速度都用立方Ｂ样条函数来进行参数化。合成例子表明，当同时考虑穿透和反射时间，层析成像的结果得到改进。这个方法也应用到英国煤炭测量局（ＢｒｉｔｉｓｈＣｏａｌＭｅａｓｕｒｅｓ）沿跨线排列所记录的逆ＶＳＰ数据组。通过使用波形配合技术，用同时确定时间延迟和叠加权，可以自动拾取旅行时间。所观察到的逆ＶＳＰ层析成像可比周围介质具有较低速度的两个断层区域成像。断层的位置由附近的反射测线所确定。本文还讨论了在复杂２－Ｄ和３－Ｄ非均匀各向同性介质中地震射线跟踪方法。界面的几何形状和水平速度场都通过使用非均匀步长立方Ｂ－样条节点 相似文献
20.
本文从能量泛函的角度着手,基于速度梯度,采用最速下降法推导了基于偏微分方程的速度模型平滑公式,用于解决射线类偏移成像过程当中速度模型的平滑处理问题.同时针对偏微分方程速度模型平滑系数中阈值k的选取对速度模型空间结构的影响,在简单的高速体速度模型上分析了不同k值的选取对原始速度模型空间结构的改变,并通过射线路径和时间场的对比分析证明偏微分方程速度模型平滑处理相对于卷积算子平滑的优越性.最后通过在Marmousi、Sigsbee2A原始速度模型以及平滑处理后的速度模型上的成像结果分析验证偏微分方程模型平滑的有效性. 相似文献