Calculation of the Inverse of the Covariance

In reservoir characterization, the covariance is often used to describe the spatial correlation and variation in rock properties or the uncertainty in rock properties. The inverse of the covariance, on the other hand, is seldom discussed in geostatistics. In this paper, I show that the inverse is required for simulation and estimation of Gaussian random fields, and that it can be identified with the differential operator in regularized inverse theory. Unfortunately, because the covariance matrix for parameters in reservoir models can be extremely large, calculation of the inverse can be a problem. In this paper, I discuss four methods of calculating the inverse of the covariance, two of which are analytical, and two of which are purely numerical. By taking advantage of the assumed stationarity of the covariance, none of the methods require inversion of the full covariance matrix.     

利用GCV优化的正则化方法对硬X射线望远镜图象重建研究

图象重建是一个求逆过程具有不适定性。本文介绍的是一种从少数投影数据进行图象重建的线性代数法－Ｐｈｉｌｌｉｐｓ－Ｔｉｋｈｏｎｏｖ正则化方法,同时选择优化工具ＧＣＶ确定最优的正则参数。为验证该方法,我们对日本Ｙｏｈｋｏｈ卫星上塔载的硬Ｘ射线望远镜（ＨＸＴ）拍摄的太阳图象进行重建。     

一种高频雷达回波谱反演海浪的新算法

在高频地波雷达海浪谱反演问题中,广泛采用的Barrick后向散射公式属于第一类非线性Fredholm积分方程。此类积分方程的解在本质上是不适定的,加之高频雷达二阶回波信号信噪比较低,使得反演海浪谱存在解不稳定的问题。本文提出一种稳定且低复杂度的反演算法,此算法首先根据高频雷达一阶回波谱测量海浪方向,并将其引入积分方程求解过程,减少求解变量的个数,降低反演算法的复杂度。为解决反演结果不稳定的问题,使用Tikhonov正则化方法并利用广义交叉验证法(GCV)确定其正则化系数。通过在不同测试条件下对反演算法的仿真测试分析,表明此方法具有运算量小、稳定性好的特点。     

基于最小支持的2.5维直流电阻率正则化反演研究

地球物理反问题存在多解性和不稳定性,正则化是得到稳定解的重要手段之一,正则化方法中稳定因子的设计与正则化因子的选择是两项重要研究内容。稳定因子的主要功能是对模型的解空间进行限制,以减少多解性,求得稳定解;为更好的定位地下目标体,研究工作采用了最小支持稳定因子,并给出了统一表达形式。正则化因子在模型误差与数据误差之间起到平衡作用,研究工作在＂L曲线法＂自动选取正则化因子的基础上,提出一种修正的计算正则化因子算法。以2.5维直流电阻率为例,运用非线性共轭梯度法,进行了试算研究,结果表明了修正的＂L曲线法＂算法快速、稳定,最小支持稳定因子具有较好的聚焦特性,有利于实现陡变边界反演。     

机载激光雷达数据的建筑物三维模型重建

机载激光雷达技术可以快速获取建筑物高精度高密度的三维空间信息.文章提出一种建筑物三维模型重建方法,该方法利用栅格数据使用RANSAC算法提取建筑边界上的长线段,采用规则化策略获取建筑轮廓;在模型生成阶段,基于原始点云和屋顶拓扑关系约束,获取屋顶间的交线和交点,并结合已得到的建筑边界、屋顶间交线和交点获取三维建筑模型.实验结果表明,该方法思路清晰、科学可行且适用性强,提取的建筑边界和模型垂直精度高.     

The Kepler problem: A unifying view

This paper is an attempt to bring unity in the study of the classical Kepler problem by combining, through simple vectorial and quaternionic techniques, its two peculiar aspects: the determination of the constants of the motion and the regularization at the origin.Research supported by the Consiglio Nazionale delle Ricerche of Italy (C.N.R.-G.N.F.M.).     

一种基于流形正则化的半监督指纹定位算法

针对传统指纹定位算法采集带标签训练数据成本高的问题,本文提出了一种基于流形正则化的半监督指纹定位算法。首先以流形假设为依据,利用批量输入的带标签数据与无标签数据之间的相似度构建图拉普拉斯算子;然后与极限学习机算法相结合,通过随机特征映射建立隐含层;最后在流形正则化框架下,求解隐含层和输出层之间的权值矩阵,从而建立位置估计模型。仿真结果表明,与INN、SVR、ELM 3种算法相比,该算法的训练和测试时间相对较短,且在带标签训练数据稀疏的前提下仍能保持较高的准确率与稳定性。     

基于机载LiDAR数据的建筑物轮廓规则化方法

在点云数据分类和建筑物轮廓提取的基础上,提出一种基于最小外包矩形的建筑物轮廓规则化方法。以任意两相邻轮廓点构建初始外包矩形,以迭代的方式将内点数最多者作为最终外包矩形,继而以面积比率、距离和均方根作为指标对规则化结果给出评判,最后以实际数据进行实验。实验结果表明,文中提出的方法能够有效对分类后的建筑物轮廓点进行规则化。     

LSQR法在位场反演中的分析与评价

LSQR法具有计算效率高、对计算机内存要求低的优点,适合于大规模问题的求解。为探讨其应用于位场反演的稳定性和可靠性,笔者以加入不同噪声的两个合成模型数据为实验对象,比较分析了Tikhonov正则化与LSQR法求解结果,显示直接利用LSQR法求解位场反问题能够得到满意的正则化解,其解模型相对Tikhonov正则化,最大相对误差仅为0.36%,说明直接利用LSQR法求解位场反问题是可行的。将其应用于四川盆地雅安地区重力三维反演,极大地降低计算成本,获取了区内沉积盆及主要断裂分布情况,为页岩气靶区优选提供了有力支撑。     

扰动重力梯度的球冠谐分析建模

从球冠谐理论出发,详细推导了球冠坐标系下扰动重力梯度的无奇异性计算公式。基于Tikhonov正则化方法,利用GOCE卫星实际观测数据解算局部重力场球冠谐模型。数值计算表明,基于扰动重力梯度的球冠谐分析建模方法能够有效地恢复局部重力场中的短波信号,与GO_CONS_GCF_2_DIR_R5模型的差异在±0.3×10~(-5) m/s~2水平。