Numerical simulation of a stratigraphic model

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.     

求解变分不等式的Newton迭代的半局部收敛性分析

分析了求解变分不等式Newton方法的半局部收敛性，建立了类似于Kantorovich定理的收敛性结果。该结果不仅为判断Newton方法的收敛性提供了可计算的充分条件，也给出了Newton方法的收敛域以及问题解的存在区域。同时，文章也得到了Newton方法的若干收敛性质，包含收敛阶以及可计算的误差估计式等。     

GPS精密星历的轨道内插方法比较

在高精度的GPS测量数据处理过程中，获取高精度卫星轨道是重要的环节，内插是获取任意历元的精密轨道信息的重要手段，对拉格朗日插值法、牛顿插值法、Neville逐次线性内插法和切比雪夫插值法的效果做了细致的比较，并在此基础上进一步探讨了短时间外推的效果。     

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

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

不等式约束Partial EIV模型的WHP拟牛顿修正解法及其精度评定的SUT法

提出一种确定不等式约束Partial EIV模型解及精度评定的新方法，在总体最小二乘准则下，将附有不等式约束的Partial EIV模型转换为标准最优化问题。采取WHP拟牛顿修正的SQP方法求解，并利用SUT法对参数估值进行精度评定，可以减小迭代次数、提高收敛速度，且精度评定方法简单有效。     

一种快速有效的大数据区域网平差方法

针对摄影测量影像来源多样化、复杂化、大数据化等趋势,传统区域网平差算法在应对当前复杂多变的数据来源,矩阵排列毫无规律的法方程结构以及大数据量带来的高内存需求和低计算效率等问题上,遇到了前所未有的挑战,为了解决上述难题,本文引入了预条件共轭梯度法以及不精确牛顿解法求解区域网平差过程中的法方程,同时使用一种块状法方程系数矩阵压缩存储格式,构建了全新的区域网平差技术流程。本文方法避免了直接对法方程系数矩阵的求逆,压缩了法方程系数矩阵所需的内存空间,使得本文算法比传统算法所需计算机内存空间大幅减少,平差计算速度明显提升,同时保证了计算精度与传统方法相当。初步试验证明,本文方法对4500张影像、近900万像点数据的平差计算在普通电脑上仅需要约15min,且计算精度达到子像素级。     

二维视电阻率断面的快速最小二乘反演

快速最小二乘反演是以平滑限定的最小二乘方法为基础,是对二维视电阻率断面进行反演的一种方法。反演过程不需要提供初始模型,在首次迭代时使用一均匀介质地下模型作为初始模型,该模型的视电阻率偏导数值可以用解析法得到。在后面的迭代中,使用了拟牛顿法去修改每一次迭代的偏导数矩阵,避免了偏导数矩阵的直接计算,从而减少了计算时间和存储空间。同时运用牛顿矩阵校正技术解最小二乘方程组也减少了大量的计算时间。总之,该方     

南极洲纳尔逊冰帽的某些动力学特征*

纳尔逊冰帽是南设得兰群岛的一个小冰帽,受海洋性气候影响冰帽上降水丰富,冰温较高。冰帽表面运动速度完全由冰帽表面形态和冰帽底床形态所控制,其中E剖面更为复杂。冰帽驱动应力基本小于100kPa.由V/Z～Za曲线得到在V/7为0.6～3×10^-9范围内的流动参数:n约为1;B≈4.3×10¹⁰dynescm^-2sec.这表明在低应力区,冰体流动更接近牛顿流体。     

基于截断牛顿法的VTI介质声波多参数全波形反演

不同类别参数间的相互耦合使多参数地震全波形反演的非线性程度显著增加,地震波速度与各向异性参数取值数量级的巨大差异也会使反演问题的性态变差.合理使用Hessian逆算子可以减弱这两类问题对反演的影响,提高多参数反演的精度,而截断牛顿法是一种可以比较准确地估计Hessian逆算子的优化方法.本文采用截断牛顿法在时间域进行了VTI介质的声波双参数同时反演的研究.不同模型的反演试验表明,在VTI介质声波双参数同时反演中,截断牛顿法比有限内存BFGS(Limited-memory Broyden-Fletcher-Goldfarb-Shanno,L-BFGS)法能更准确地估计Hessian逆算子,进而较好地平衡两类不同参数的同时更新,得到了比较精确的反演结果.     

航空重力测量中牛顿中心差分器对确定载体加速度的影响

讨论不同点数的牛顿中心差分器对航空重力测量中确定载体加速度的影响。理论分析表明,当航空重力测量的有效频段位于牛顿中心差分器的通带时,不同点数的牛顿中心差分器所得结果经低通滤波后几乎没有差异。实验表明,当采样间隔为1 s、载体平均速度为155 m/s、所需重力数据的分辨率为10 km时,低通滤波后不同点数的牛顿中心差分器输出结果之间的较差优于2.1 mGal。