基于多尺度迭代求解三维频率-空间域声波方程的全波形速度反演研究-理论模型

使用广义最小残量方法迭代求解三维频率-空间域声波方程,反演时使用多尺度、多重网格的策略,探讨了如何快速实现高分辨率的三维频率-空间域迭代法声波全波形速度反演.通过对理论模型进行三维频率-空间域迭代法声波全波形反演数值试验,证实该方法的计算速度快、计算效率高,反演所得速度的分辨率高.从而为基于多尺度迭代求解三维频率-空间域声波方程的全波形速度反演成像打下方法基础.     

一维声波方程联合反演的分步正则化法

本文用一个纵波信息，对一维声波方程的速度和源函数进行联合反演，并考虑到声波方程的反问题是一个不适应问题，对源函数和波速分别和正则化法分步迭代求解，减少反问题的计算工作量，改善该问题的计算稳定性，为计算实际工程和岩性学问题供了一种方法。文中给出只用一个反问题补充条件同时进行多参数反演的公式，并对相应的数值算例进行分析和比较。     

小波变换多尺度地震波形反演

完全非线性地震波形反演问题是石油地球物理勘探领域中一个非常重要而又难度很大的问题。本文提出了多尺度地震波形反演的小波变换方法,对于一维非线性地震波形反演问题,将该方法和已有的简单迭代法及多重网格法作了比较,数值实验结果表明,本方法效果较好。     

二维大地电磁数据的整体反演

从麦克斯韦第一、二方程出发，经过变量替换，将二维电磁问题线性化，得到不同测点、不同浓度下关于导电率微小变化的线性方程组，用Gauss-Seidel迭代法求解。实现了整体对二维大地电磁数据的反演，由于对各个测点下的数据同时反演，克服了以往用一维方法来反演二维问题时固有的缺点，讨论了不同的初始模型对反演结果的影响，提出了一种简便实用的初始值的给定方法，数值模拟结果表明该方法有一定的抗噪能力。     

一维波动方程波阻抗反演的同伦方法

文中从地震勘探一维波动方程反问题出发,研究了一种反演地层参数的同伦方法,该方法把非线性方程组的求解转化成常微分方程初值问题的数值求解,从而给出一种稳定的计算速度快、抗噪能力强的全局收敛的反演方法．理论模型和实例试算的结果表明了同伦方法是一种有效的反演算法,特别适用于非线性、多极值的地球物理反演问题,在地球物理非线性反演中具有广泛的应用前景．     

介质密度反演偏导矩阵的精确计算方法

实现反演偏导矩阵的计算是基于导数最优化反演方法的关键,然而目前的地震反演几乎都是基于Zoeppritz方程近似实现的,使计算精度和适应范围受到限制.本文利用Zoeppritz方程建立了反射系数对地层介质密度比偏导方程,导出了Zoeppritz方程矩阵元对介质密度比的导数.通过求解偏导方程获得了反射系数对介质密度比偏导数的精确计算(考虑了速度中含介质密度的问题).利用数值算例分析了反射系数对介质密度比偏导数的变化特点.本文采用直接解法求解偏导矩阵方程组,获得了快的计算速度和高的计算精度,为实现地层介质密度反演(包括大角度反演)提供了偏导矩阵的计算方法.     

散射波积分方程的Adomian分解解法

在背景模型基础上,求解模型扰动后的地震波散射场,这是目前地震反演中的一个关键步骤.本文将计算数学中求解非线性积分方程的Adomian分解方法,应用到求解标量波散射场的Lippmann-Schwinger积分方程和Ricatti积分方程中,分别得到了散射场的Born序列解和Rytov序列解.通过一维和二维数值算例说明:在满足一定的条件下,散射场的这两种序列解稳定收敛,与传统的Born和Rytov近似解相比,引入散射序列中的高阶项可以更精确地描述地震波散射场.     

基于Zoeppritz偏导方程精确解的地层密度多角度反演

基于Zoeppritz方程对介质密度偏导数所建立的偏导方程的精确解,构造了多角度反演地层介质密度的反演方程,在偏导数求解过程中考虑了介质密度对波速度的影响因素,并由此实现了利用反射系数梯度精确解计算地层密度的多角度联合反演.通过数值算例考察了计算方法,结果显示:反演方法对层状地层模型不论反射波是否存在相干现象均获得了较好的反演结果,反演迭代10次后计算结果的最大相对误差能够收敛到1%之内;随着反演角度的增加地层介质密度反演的精度逐步提高,反演具有自动校正能力,有快的计算速度.本方法克服了传统AVO(Amplitude Versus Offset)基于Zoeppritz方程近似所遇到的困难,不受反演角度大小及反射界面对波反射强弱的限制,为地层介质密度的多角度包括大角度反演提供了一种新的快速有效的计算方法.     

Multi-angle inversion of formation densities based on the accurate solutions of Zoeppritz's partial derivative equations

基于Zoeppritz方程对介质密度偏导数所建立的偏导方程的精确解,构造了多角度反演地层介质密度的反演方程,在偏导数求解过程中考虑了介质密度对波速度的影响因素,并由此实现了利用反射系数梯度精确解计算地层密度的多角度联合反演.通过数值算例考察了计算方法,结果显示:反演方法对层状地层模型不论反射波是否存在相干现象均获得了较好的反演结果,反演迭代10次后计算结果的最大相对误差能够收敛到1％之内;随着反演角度的增加地层介质密度反演的精度逐步提高,反演具有自动校正能力,有快的计算速度.本方法克服了传统AVO(Amplitude Versus Offset)基于Zoeppritz方程近似所遇到的困难,不受反演角度大小及反射界面对波反射强弱的限制,为地层介质密度的多角度包括大角度反演提供了一种新的快速有效的计算方法.     

频率域波动方程的双参数识别方法

利用反射波地震记录恢复波动方程的两个系数以研究地球内部构造，是当前应用地球物理研究的热点之一，本文在波动方程反演方法的理论基础上，导出了频率域识别介质的密度波速两个重要参数的新算法，从而把当前的波动方程单参数反演双参数反演方法，通过理论模型及实际测井资料的反演试算对比，均藜得良好效果，在算法中还省去了Ｇｒｅｅｎ函数的繁杂计算，节省了计算量，增加实用性。     

Ground penetrating radar inversion in 1-D: an approach for the estimation of electrical conductivity, dielectric permittivity and magnetic permeability

This paper presents a method for inverting ground penetrating radargrams in terms of one-dimensional profiles. We resort to a special type of linearization of the damped E-field wave equation to solve the inverse problem. The numerical algorithm for the inversion is iterative and requires the solution of several forward problems, which we evaluate using the matrix propagation approach. Analytical expressions for the derivatives with respect to physical properties are obtained using the self-adjoint Green's function method. We consider three physical properties of materials; namely dielectrical permittivity, magnetic permeability and electrical conductivity. The inverse problem is solved minimizing the quadratic norm of the residuals using quadratic programming optimization. In the iterative process to speed up convergence we use the Levenberg–Mardquardt method. The special type of linearization is based on an integral equation that involves derivatives of the electric field with respect to magnetic permeability, electrical conductivity and dielectric permittivity; this equation is the result of analyzing the implication of the scaling properties of the electromagnetic field. The ground is modeled using thin horizontal layers to approximate general variations of the physical properties. We show that standard synthetic radargrams due to dielectric permittivity contrasts can be matched using electrical conductivity or magnetic permeability variations. The results indicate that it is impossible to differentiate one property from the other using GPR data.     

Iterative convergence of boundary-volume integral equation method

The boundary-volume integral equation numerical technique can be a powerful tool for piecewise heterogeneous media, but it is limited to small problems or low frequencies because of great computational cost. Therefore, a restarted GMRES method is applied to solve large-scale boundary-volume scattering problems in this paper to overcome the computational barrier. The iterative method is firstly applied to responses of dimensionless frequency to a semicircular alluvial valley filled with sediments, compared with the standard Gaussian elimination method. Then the method is tested by a heterogeneous multilayered model to show its applicability. Numerical experiments indicate that the preconditioned GMRES method can significantly improve computational efficiency especially for large Earth models and high frequencies, but with a faster convergence for the left diagonal preconditioning.     

A central difference method with low numerical dispersion for solving the scalar wave equation

In this paper, we propose a nearly‐analytic central difference method, which is an improved version of the central difference method. The new method is fourth‐order accurate with respect to both space and time but uses only three grid points in spatial directions. The stability criteria and numerical dispersion for the new scheme are analysed in detail. We also apply the nearly‐analytic central difference method to 1D and 2D cases to compute synthetic seismograms. For comparison, the fourth‐order Lax‐Wendroff correction scheme and the fourth‐order staggered‐grid finite‐difference method are used to model acoustic wavefields. Numerical results indicate that the nearly‐analytic central difference method can be used to solve large‐scale problems because it effectively suppresses numerical dispersion caused by discretizing the scalar wave equation when too coarse grids are used. Meanwhile, numerical results show that the minimum sampling rate of the nearly‐analytic central difference method is about 2.5 points per minimal wavelength for eliminating numerical dispersion, resulting that the nearly‐analytic central difference method can save greatly both computational costs and storage space as contrasted to other high‐order finite‐difference methods such as the fourth‐order Lax‐Wendroff correction scheme and the fourth‐order staggered‐grid finite‐difference method.     

利用积分方程法计算双侧向测井仪的井眼校正曲线

对电测井积分方程中的各积分项进行了理论推导，利用任意闭合曲面对曲面外一点所张立体角为零的特点将电势积分项由对柱体侧面的积分转化为对柱体底面的积分，从而使在纵向采用线性插值条件下给出了电势积分项的解析解．对电流积分项进行了化简，使之由曲面积分化为一维数值积分．利用所得结果计算了双侧向测井仪的井眼校正曲线．     

Lippmann-Schwinger积分方程广义超松弛迭代解法及其收敛特性

为克服Born级数只适用于弱散射和短程传播的弱点,从光学散射理论中引入了针对强光学散射问题的Lippmann-Schwinger方程广义超松弛迭代解法,并对这种方法在反射地震条件下的算法表现、收敛特点、计算效率以及计算精度等问题进行了分析和讨论.理论分析和数值计算结果均表明:(1)与散射级数重整化相比,超松弛法具有更坚实的数学基础;(2)超松弛法可以有效地克服Born级数的弱点,得到与有限差分相当的数值模拟结果;(3)在反射地震条件下LS方程超松弛法表现优异,完全可以解决反射地震数值模拟中的强散射问题.     

Fast and accurate three‐dimensional controlled source electromagnetic modelling†

We present a fast approximate method for three‐dimensional low frequency controlled source electro‐magnetic modeling. We apply the method to a synthetic model in a typical marine controlled source electromagnetic scenario, where conductivity and permittivity are different from the known background medium. For 3D configurations, fast computational methods are relevant for both forward and inverse modelling studies. Since this problem involves a large number of unknowns, it has to be solved efficiently to obtain results in a timely manner, without compromising accuracy. For this reason, the Born approximation, extended Born approximation and iterative extended Born approximation are implemented and compared with the full solution of the conjugate gradient fast Fourier transformation method. These methods are based on an electric field domain integral equation formulation. It is shown here how well the iterative extended Born approximation method performs in terms of both accuracy and speed with different configurations and different source positions. The improved accuracy comes at virtually no additional computational cost. With the help of this method, it is now possible to perform sensitivity analysis using 3D modelling in a timely manner, which is vital for controlled source electromagnetic applications. For forward modeling the solution at the sea‐bottom is of interest, because that is where the receivers are usually located. For inverse modeling, the accuracy of the solution in the target zone is important to obtain reasonably accurate conductivity values from the inversion using this approximate solution method. Our modelling studies show that the iterative extended Born approximation method is fast and accurate for both forward and inverse modelling. Sensitivity analysis as a function of the source position and different reservoir sizes validate the accuracy of the iterative extended Born approximation.     

具有滞迟耗能装置的高层建筑结构的随机响应分析

本文提出了一个近似分析具有滞迟耗能装置的高层建筑的方法。结果被模拟为一个有滞迟非线性支座的悬壁梁。其控制方程由Ｒｉｔｚ法导出，用等效线性化方法分析。     

广义WKB法在波场模拟中的应用研究

针对地震勘探中的波场模拟和反演困难等问题,本文从二维波动方程出发,采用Fourier变换得到了赫姆霍兹方程,并创造性地提出了一种全新的近似解法－广义WKB法;该法可以更加近似赫姆霍兹方程,且它的解中包含WKB的近似值.数值模拟显示了WKB法与广义WKB求解的异同.结果表明：广义WKB法是WKB法的一种推广,在求解非线性波数问题时具有更好的优势.     

用于弹性波方程数值模拟的有限差分系数确定方法

有限差分方法是波场数值模拟的一个重要方法,交错网格差分格式比规则网格差分格式稳定性更好,但方法本身都存在因网格化而形成的数值频散效应,这会降低波场模拟的精度与分辨率.为了缓解有限差分算子的数值频散效应,精确求解空间偏导数,本文把求解波动方程的线性化方法推广到用于求解弹性波方程交错网格有限差分系数;同时应用最大最小准则作为模拟退火(SA)优化算法求解差分系数的数值频散误差判定标准来求解有限差分系数.通过上述两种方法,分别利用均匀各向同性介质和复杂构造模型进行了数值正演模拟和数值频散分析,并与传统泰勒展开算法、最小二乘算法进行比较,验证了线性化方法和模拟退火方法都能有效压制数值频散,并比较了各个算法的特点.