测距定位方程非线性平差的封闭牛顿迭代公式

测距方程;非线性强度;最小二乘;高斯-牛顿法;牛顿法;病态方程     

GPS测码伪距绝对定位的几种算法

GPS定位方程是非线性的,一般处理方法是按泰勒级数展开取至一次项进行线性化,再利用最小二乘原理求解,如果所取观测站坐标的初始值具有较大的偏差,略去二次微小量的模型误差,对解算结果将产生不能忽略的影响。本文研究分析了GPS测码伪距绝对定位的传统算法,并提出了一种通过求差法将GPS绝对定位的非线性观测方程转化成线性方程直接求解测站坐标的新算法,通过实例计算表明该方法计算简单,不需要测站的初始坐标信息,不需要求导计算和迭代计算,对于提高GPS测码伪距单点定位的解算速度和精度具有重要的意义。     

频率域海洋可控源电磁垂直各向异性三维反演

地层宏观电性各向异性会对可控源电磁响应产生重要影响.由于海底地层电性结构常表现为电导率各向异性,若仅对海洋可控源电磁(MCSEM)数据进行常规各向同性反演,有可能无法获得准确的反演解释结果,从而削弱MCSEM技术的可靠性.本文实现了电导率垂直各向异性(VTI)条件下频率域海洋可控源电磁数据三维反演算法.其中,三维正演采用基于二次场控制方程的交错网格有限体积法,并利用直接矩阵分解技术来求解离散所得的大型线性方程组,有利于快速计算多场源的响应.反演采用具有近似二次收敛性的高斯牛顿算法对目标函数进行最优化.最后,对具有VTI电性各向异性特征的盐丘构造模型的MCSEM合成数据分别进行了电导率各向同性和垂直各向异性三维反演,结果表明:各向同性三维反演算法无法对受VTI介质影响的MCSEM数据进行正确的反演解释,而垂直各向异性三维反演能够获得更为可靠的地下电阻率结构和异常体分布,展现出对海底电性各向异性结构更为优良的反演解释能力.     

混凝土结构缺陷检测的探地雷达资料波场反演方法

本文针对工程建筑中的混凝土无损探测问题,提出了一种新的反演成像方法.它避免了传统方法观察偏移探测数据的情况.这种方法以探地雷达为探测手段,以电磁场原理为理论依据,其结果更具定量化和可视性.具体算法是将阻尼高斯牛顿法和GCV方法(Generalized Cross Validation)结合起来,其中阻尼牛顿法可以避免传统牛顿法计算Hessian矩阵的困难,GCV方法可以自适应地选择正则化参数,既减小了计算量又克服了噪声及不适定性的影响.数值算例验证了算法的有效性.     

基于Born敏感核函数的速度、密度双参数全波形反演

速度、密度之间的相互耦合使得密度在多参数全波形反演中较难获得.本文将截断高斯-牛顿法用于声介质速度、密度双参数全波形反演,通过考虑近似Hessian矩阵中反映速度、密度相互作用的非主对角块元素,有效解决了多参数全波形反演中速度、密度之间的耦合问题,在不采用反演策略的情况下,仍能够获得精度较高的速度、密度反演结果.常规的截断牛顿类全波形反演通常利用一阶伴随状态法求取目标函数对模型参数的梯度,利用二阶伴随状态法或有限差分法求解Hessian-向量乘,在每一步内循环迭代过程中需要额外求解两次正演问题,计算量较大.本文基于Born近似,将梯度计算中的核函数-向量乘表示为具有明确物理意义的向量-标量乘的累加运算,同时将Hessian-向量乘转化为两次核函数-向量乘,无需额外求解正演问题,有效降低了计算量.数值实验证明了本文提出的方法的有效性.     

基于改进高斯-牛顿法的位场向下延拓

本文针对位场向下延拓的不适定问题,在分析最优化算法中高斯-牛顿法基本原理及滤波函数滤波特性的基础之上,提出基于正则参数指数递增计算方法和残差最小步长准则的改进高斯-牛顿法。基于理论重力模型和航磁实测数据的对比实验表明,改进后得到的自适应迭代法具有相对较高的位场向下延拓精度和很好的收敛性。     

水平层状介质中双侧向测井资料的迭代Tikhonov正则化反演

根据非线性反演理论与Morozov偏差原理研究建立从双侧向测井(DLL)资料中同时重构地层原状电阻率、侵入带电阻率、侵入半径、层界面位置以及井眼泥浆电阻率的迭代正则化算法.首先利用Tikhonov正则化反演理论将双侧向测井资料的反演问题转化为含有稳定泛函的非线性目标函数的极小化问题,并利用Gauss-Newton算法确定极小化解.为得到稳定的反演结果并有效实现测井资料的最佳拟合,在迭代过程中将Morozov偏差原理和Cholesky分解技术相结合,建立了一套后验选择正则化因子的方法.最后通过理论模型和大庆油田实际测井资料的处理结果,验证了该算法能够取得更为满意的反演效果.     

基于高斯牛顿法的频率域可控源电磁三维反演研究

三维反演解释是电磁法勘探发展的重要趋势,而如何提高三维反演的可靠性、稳定性和计算效率是算法开发者们目前的研究重点.本文实现了一种频率域可控源电磁(CSEM)三维反演算法.其中正演基于拟态有限体积法离散化,利用直接矩阵分解技术来求解大型线性系统方程,不仅准确、稳定,而且特别有利于含有大量发射场源位置的CSEM勘探情况;对目标函数的最优化采用高斯牛顿法(GN),具有近似二次的收敛性;使用预条件共轭梯度法(PCG)求解每次GN迭代所得到的法方程,避免了显式求解和存储灵敏度矩阵,减小了计算量.以上这些方法的结合应用,使得本文的三维反演算法准确、稳定且高效.通过陆地和海洋CSEM勘探场景中的典型理论模型的反演测试,验证了本文算法的有效性.     

A spectral approach to estimation and smoothing of continuous spatial processes

This paper is concerned with developing computational methods and approximations for maximum likelihood estimation and minimum mean square error smoothing of irregularly observed two-dimensional stationary spatial processes. The approximations are based on various Fourier expansions of the covariance function of the spatial process, expressed in terms of the inverse discrete Fourier transform of the spectral density function of the underlying spatial process. We assume that the underlying spatial process is governed by elliptic stochastic partial differential equations (SPDE's) driven by a Gaussian white noise process. SPDE's have often been used to model the underlying physical phenomenon and the elliptic SPDE's are generally associated with steady-state problems.A central problem in estimation of underlying model parameters is to identify the covariance function of the process. The cumbersome exact analytical calculation of the covariance function by inverting the spectral density function of the process, has commonly been used in the literature. The present work develops various Fourier approximations for the covariance function of the underlying process which are in easily computable form and allow easy application of Newton-type algorithms for maximum likelihood estimation of the model parameters. This work also develops an iterative search algorithm which combines the Gauss-Newton algorithm and a type of generalized expectation-maximization (EM) algorithm, namely expectation-conditional maximization (ECM) algorithm, for maximum likelihood estimation of the parameters.We analyze the accuracy of the covariance function approximations for the spatial autoregressive-moving average (ARMA) models analyzed in Vecchia (1988) and illustrate the performance of our iterative search algorithm in obtaining the maximum likelihood estimation of the model parameters on simulated and actual data.     

A GCV based method for nonlinear ill-posed problems

This paper discusses the inversion of nonlinear ill-posed problems. Such problems are solved through regularization and iteration and a major computational problem arises because the regularization parameter is not known a priori. In this paper we show that the regularization should be made up of two parts. A global regularization parameter is required to deal with the measurement noise, and a local regularization is needed to deal with the nonlinearity. We suggest the generalized cross validation (GCV) as a method to estimate the global regularization parameter and the damped Gauss-Newton to impose local regularization. Our algorithm is tested on the magnetotelluric problem. In the second part of this paper we develop a methodology to implement our algorithm on large-scale problems. We show that hybrid regularization methods can successfully estimate the global regularization parameter. Our algorithm is tested on a large gravimetric problem. This revised version was published online in July 2006 with corrections to the Cover Date.