共查询到19条相似文献,搜索用时 171 毫秒
1.
2.
《地球物理学进展》2017,(2)
传统的地震逆散射广义Radon变换(GRT)保幅反演方法是建立在散射场一阶Born近似(单散射)的基础上,仅仅适用于弱扰动介质模型.本文从散射场积分方程出发,通过研究二次散射的特征,讨论和验证了基于局部二阶Born近似的GRT非线性保幅反演方法,将传统GRT线性保幅反演算子的适用范围扩展至非均匀强扰动介质.数值测试结果表明:在散射场近似模拟方面,二阶Born近似比一阶Born近似更为准确,二次散射效应主要集中在主散射点周围的局部区域内,超过这一范围,二次散射强度趋于稳定;在保幅反演方面,本文基于局部二阶Born近似的GRT非线性反演算法,明显优于传统的GRT线性反演算法,可以准确重构强扰动介质模型,而计算效率与线性反演方法相当. 相似文献
3.
传统的波动方程走时核函数(或走时Fréchet导数)多基于互相关时差测量方式及地震波场的一阶Born近似导出,其成立条件非常苛刻.然而,地震波走时与大尺度的速度结构具有良好的线性关系,对于小角度的前向散射波场,Rytov近似优于Born近似.因此,本文基于Rytov近似和互相关时差测量方式,导出了基于Rytov近似的有限频走时敏感度核函数的两种等价形式:频率积分和时间积分表达式.在此基础之上,本文提出了一种隐式矩阵向量乘方法,可以直接计算Hessian矩阵或者核函数与向量的乘积,而无需显式计算和存储核函数及Hessian矩阵.基于隐式矩阵向量乘方法,本文利用共轭梯度法求解法方程实现了一种高效的Gauss-Newton反演算法求解走时层析反问题.与传统的敏感度核函数反演方法相比,本文方法在每次迭代过程中,无需显式计算和存储核函数,极大降低了存储需求.与基于Born近似的伴随状态方法走时层析相比,本文方法具有准二阶的收敛速度,且适用范围更广.数值试验证明了本文方法的有效性. 相似文献
4.
有限频率层析成像考虑了非均匀介质中波的散射、衍射、波前愈合等物理性质,使得其对速度异常体的分辨能力远大于射线层析成像.推导和计算有限频率敏感核是进行有限频率层析成像的关键,当前推导有限频率敏感核多借助一阶Born近似,但这只适用于弱散射介质的情况.本文基于二阶Born近似并利用傅里叶变换推导了三维均匀介质情况下有限频率敏感核的解析表达式,并将其推广到非均匀介质中得到了三维非均匀介质中有限频率敏感核.研究表明:当介质中速度扰动小于2%时,基于二阶Born近似的有限频率敏感核与基于一阶Born近似的有限频率敏感核差别很小,可近似认为相同;当介质中速度扰动大于5%时,基于二阶Born近似的有限频率敏感核与基于一阶Born近似的有限频率敏感核有较大不同,表明此时已不能忽略二次散射. 相似文献
5.
在背景模型基础上,求解模型扰动后的地震波散射场,这是目前地震反演中的一个关键步骤.本文将计算数学中求解非线性积分方程的Adomian分解方法,应用到求解标量波散射场的Lippmann-Schwinger积分方程和Ricatti积分方程中,分别得到了散射场的Born序列解和Rytov序列解.通过一维和二维数值算例说明:在满足一定的条件下,散射场的这两种序列解稳定收敛,与传统的Born和Rytov近似解相比,引入散射序列中的高阶项可以更精确地描述地震波散射场. 相似文献
6.
7.
本文将普遍声逆散射微扰论应用于弹性波层析成像问题,在Born变换下推出了以旋转角为补偿参数的各阶微扰重建公式,实现了对非均匀各向同性散射体内3个参数(质量密度ρ和两个Lamé系数λ,μ)的同时重建. 对于层析成像问题,在弹性波的传播过程中P波与SV波有耦合,但它们不会和SH波发生耦合,于是可以得到3个形式相对简单的标量方程. 在Born变换下,在散射波中引入微扰参数,将散射体的3个参数分别按该微扰参数展开,然后利用二维自由空间的Green函数分别得到散射的P波、SV波和SH波的积分表示. 最后,经一维傅氏变换后,得到Born变换下散射体3个参数的各阶微扰重建公式. 相似文献
8.
9.
10.
应用当前的数值方法求解动态海面与目标的复合电磁散射,由于海面的变化,在不同时刻需要对阻抗矩阵各元素及海面表面电流重新求解,因而要耗费大量内存和运算量.为了克服这一问题,本文应用物理光学(PO)近似求解了导体海面表面电流及导体平板的一阶散射场,应用基尔霍夫近似给出了海面的后向散射场,同时借助互易性定理降低了求解平板和海面之间二次耦合散射场的难度,讨论了平板尺寸、风速等对后向复合散射场的影响.另外,本文还推导出了耦合散射场多普勒谱频移的理论公式,详细分析了复合后向散射场的Doppler 频谱特性. 相似文献
11.
An analytical approximation for the calculation of the stationary reliability of linear dynamic systems with higher‐dimensional output under Gaussian excitation is presented. For systems with certain parameters theoretical and computational issues are discussed for two topics: (1) the correlation of failure events at different parts of the failure boundary and (2) the approximation of the conditional out‐crossing rate across the failure boundary by the unconditional one. The correlation in the first topic is approximated by a multivariate integral, which is evaluated numerically by an efficient algorithm. For the second topic some existing semi‐empirical approximations are discussed and a new one is introduced. The extension to systems with uncertain parameters requires the calculation of a multi‐dimensional reliability integral over the space of the uncertain parameters. An existing asymptotic approximation is used for this task and an efficient scheme for numerical calculation of the first‐ and second‐order derivatives of the integrand is presented. Stochastic simulation using an importance sampling approach is also considered as an alternative method, especially for cases where the dimension of the uncertain parameters is moderately large. Comparisons between the proposed approximations and Monte Carlo simulation for some examples related to earthquake excitation are made. It is suggested that the proposed analytical approximations are appropriate for problems that require a large number of consistent error estimates of the probability of failure, as occurs in reliability‐based design optimization. Numerical problems regarding computational efficiency may arise when the dimension of both the output and the uncertain parameters is large. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
12.
Theory of equivalent staggered‐grid schemes: application to rotated and standard grids in anisotropic media 下载免费PDF全文
The previous finite‐difference numerical schemes designed for direct application to second‐order elastic wave equations in terms of displacement components are strongly dependent on Poisson's ratio. This fact makes theses schemes useless for modelling in offshore regions or even in onshore regions where there is a high Poisson's ratio material. As is well known, the use of staggered‐grid formulations solves this drawback. The most common staggered‐grid algorithms apply central‐difference operators to the first‐order velocity–stress wave equations. They have been one of the most successfully applied numerical algorithms for seismic modelling, although these schemes require more computational memory than those mentioned based on second‐order wave equations. The goal of the present paper is to develop a general theory that enables one to formulate equivalent staggered‐grid schemes for direct application to hyperbolic second‐order wave equations. All the theory necessary to formulate these schemes is presented in detail, including issues regarding source application, providing a general method to construct staggered‐grid formulations to a wide range of cases. Afterwards, the equivalent staggered‐grid theory is applied to anisotropic elastic wave equations in terms of only velocity components (or similar displacements) for two important cases: general anisotropic media and vertical transverse isotropy media using, respectively, the rotated and the standard staggered‐grid configurations. For sake of simplicity, we present the schemes in terms of velocities in the second‐ and fourth‐order spatial approximations, with second‐order approximation in time for 2D media. However, the theory developed is general and can be applied to any set of second‐order equations (in terms of only displacement, velocity, or even stress components), using any staggered‐grid configuration with any spatial approximation order in 2D or 3D cases. Some of these equivalent staggered‐grid schemes require less computer memory than the corresponding standard staggered‐grid formulation, although the programming is more evolved. As will be shown in theory and practice, with numerical examples, the equivalent staggered‐grid schemes produce results equivalent to corresponding standard staggered‐grid schemes with computational advantages. Finally, it is important to emphasize that the equivalent staggered‐grid theory is general and can be applied to other modelling contexts, e.g., in electrodynamical and poroelastic wave propagation problems in a systematic and simple way. 相似文献
13.
本文首先解决了声波方程的非Born近似的正演计算问题,从而获得理论上不带近似的正演数据;然后,推导了井间(CBP)、垂直地震剖面(VSP)和地面反射(SRP)三种不同的数据采集方式下的衍射CT的重建公式;利用这些重建算法和正演数据,系统地研究了影响到地球物理CT成象质量的三种因素,即:(1)数据采集方式,(2)异常程度和(3)成象区域的尺寸,对重建图象的影响;并比较了衍射地震CT和射线地震CT的成象质量。 相似文献
14.
This paper discusses Born/Rytov approximation tomographic velocity inversion methods constrained by the Fresnel zone. Calculations of the sensitivity kernel function and traveltime residuals are critical in tomographic velocity inversion. Based on the Born/Rytov approximation of the frequency-domain wave equation, we derive the traveltime sensitivity kernels of the wave equation on the band-limited wave field and simultaneously obtain the traveltime residuals based on the Rytov approximation. In contrast to single-ray tomography, the modified velocity inversion method improves the inversion stability. Tests of the near-surface velocity model and field data prove that the proposed method has higher accuracy and Computational efficiency than ray theory tomography and full waveform inversion methods. 相似文献
15.
We present preconditioned non‐linear conjugate gradient algorithms as alternatives to the Gauss‐Newton method for frequency domain full‐waveform seismic inversion. We designed two preconditioning operators. For the first preconditioner, we introduce the inverse of an approximate sparse Hessian matrix. The approximate Hessian matrix, which is highly sparse, is constructed by judiciously truncating the Gauss‐Newton Hessian matrix based on examining the auto‐correlation and cross‐correlation of the Jacobian matrix. As the second preconditioner, we employ the approximation of the inverse of the Gauss‐Newton Hessian matrix. This preconditioner is constructed by terminating the iteration process of the conjugate gradient least‐squares method, which is used for inverting the Hessian matrix before it converges. In our preconditioned non‐linear conjugate gradient algorithms, the step‐length along the search direction, which is a crucial factor for the convergence, is carefully chosen to maximize the reduction of the cost function after each iteration. The numerical simulation results show that by including a very limited number of non‐zero elements in the approximate Hessian, the first preconditioned non‐linear conjugate gradient algorithm is able to yield comparable inversion results to the Gauss‐Newton method while maintaining the efficiency of the un‐preconditioned non‐linear conjugate gradient method. The only extra cost is the computation of the inverse of the approximate sparse Hessian matrix, which is less expensive than the computation of a forward simulation of one source at one frequency of operation. The second preconditioned non‐linear conjugate gradient algorithm also significantly saves the computational expense in comparison with the Gauss‐Newton method while maintaining the Gauss‐Newton reconstruction quality. However, this second preconditioned non‐linear conjugate gradient algorithm is more expensive than the first one. 相似文献
16.
17.
18.
针对传统射线层析存在的种种局限性,菲涅尔带走时层析成像摒弃了传统的数学射线,考虑到地震信号具有一定的频带宽度,中央射线附近的介质对地震波的传播产生不同程度的影响。本文提出了多频段组合菲涅尔带走时层析成像方法。该方法以频率域波动方程Born和Rytov近似为基础,推导出建立在带限地震波理论基础上的波动方程 Rytov 近似走时敏感核函数,实现第一菲涅尔带约束下的波动方程走时层析反演方法。同时由于多个频段的引入,充分利用低频段和高频段的特有优势,从而兼顾菲涅尔带层析的计算效率与分辨率。模型试算结果证明了本方法的有效性和稳定性。 相似文献
19.
M. Emin Candansayar 《Geophysical Prospecting》2008,56(1):141-157
I investigated the two‐dimensional magnetotelluric data inversion algorithms in studying two significant aspects within a linearized inversion approach. The first one is the method of minimization and second one is the type of stabilizing functional used in parametric functionals. The results of two well‐known inversion algorithms, namely conjugate gradient and the least‐squares solution with singular value decomposition, were compared in terms of accuracy and CPU time. In addition, magnetotelluric data inversion with various stabilizers, such as L2‐norm, smoothing, minimum support, minimum gradient support and first‐order minimum entropy, were examined. A new inversion algorithm named least‐squares solution with singular value decomposition and conjugate gradient is suggested in seeing the outcomes of the comparisons carried out on least‐squares solutions with singular value decomposition and conjugate gradient algorithms subject to a variety of stabilizers. Inversion results of synthetic data showed that the newly suggested algorithm yields better results than those of the individual implementations of conjugate gradient and least‐squares solution with singular value decomposition algorithms. The suggested algorithm and the above‐mentioned algorithms inversion results for the field data collected along a line crossing the North Anatolian Fault zone were also compared each other and results are discussed. 相似文献