基于VB的测量坐标正算与反算研究

坐标正反算在测量实际工作中经常用到,本文介绍了测量坐标正算与反算的数学模型,采用VB开发语言,设计了测量坐标正算与反算程序,为测绘工作中坐标转换提供了思路。     

Representations of inverse covariances by differential operators

In the cost function of three- or four-dimensional variational data assimilation, each term is weighted by the inverse of its associated error covariance matrix and the background error covariance matrix is usually much larger than the other covariance matrices. Although the background error covariances are traditionally normalized and parameterized by simple smooth homogeneous correlation functions, the covariance matrices constructed from these correlation functions are often too large to be inverted or even manipulated. It is thus desirable to find direct representations of the inverses of background error correlations. This problem is studied in this paper. In particular, it is shown that the background term can be written into f dx|Dv(x)|^2, that is, a squared L2 norm of a vector differential operator D, called the D-operator, applied to the field of analysis increment v(x). For autoregressive correlation functions, the D-operators are of finite orders. For Gaussian correlation functions, the D-operators are of infinite order. For practical applications, the Gaussian D-operators must be truncated to finite orders. The truncation errors are found to be small even when the Gaussian D-operators are truncated to low orders. With a truncated D-operator, the background term can be easily constructed with neither inversion nor direct calculation of the covariance matrix. D-operators are also derived for non-Gaussian correlations and transformed into non-isotropic forms.     

FURTHER STUDY ON THE PROPERTIES OF OPERATORS OF ATMOSPHERIC EQUATIONS AND THE EXISTENCE OF ATTRACTOR

The equivalent operator equation is derived from the full primitive nonlinear equations of theatmospheric motion and the properties and physical senses of the operators are studied.In theinfinite dimensional Hilbert space,the global asymptotic behavior of the atmosphere system withthe non-stationary external forcing is studied under the assumption of the bounded externalforcing.The existence theorems of the global absorbing set and the global attractor are obtained.Thus,the conclusions deduced from the large-scale atmosphere(Li and Chou 1996 a;1996 b)areextended to the general atmosphere.     

一种基于角点特征的遥感影像自动配准方法

本文通过对Harris算子进行改进,自适应地提取角点特征,然后基于角点特征进行由粗到精的匹配,完成影像的配准;并利用MODIS影像、TM影像和HJ-1卫星影像3种不同分辨率的遥感影像进行配准实验,结果表明该自动配准方法能够达到亚像素级的配准精度,是一种高精度的影像配准方法。     

基于深部位移的蠕滑型滑坡预警指标及预警预报研究

基于边坡监测数据建立数学模型,是边坡变形和稳定性分析的重要方法。但是单一预测模型的形式和应用范围具有一定的确定性,不同模型对数据的利用程度也有所差别,往往不能充分运用已知信息,导致模型精度不高,适用性不强。针对单一预测模型存在的问题,提出一种基于熵权法的PSO-SVR-NGM优化组合模型。该模型结合高精度变权缓冲NGM(1,1,k,c)模型和PSO-SVR模型,能够减小单一预测模型的误差,大幅度提高预测精度。首先通过引入变权缓冲算子λ和背景值权重系数η、κ改进无偏NGM(1,1,k,c)模型,构建新的3参数变权缓冲NGM(1,1,k,c)模型。结合最大灰色关联度和最小平均相对拟合误差重新构造粒子群算法的适应度函数,利用改进的粒子群算法对提出的变权缓冲模型进行搜索寻优,确定最佳的参数组合。然后通过熵权法对改进的变权缓冲NGM(1,1,k,c)模型和PSO-SVR模型进行赋权建立优化组合模型。最后,将该组合模型应用于3个不同变形特征的边坡工程中,并与其他单一模型进行对比分析。结果表明,相对于单一模型,本文所提出的组合模型的拟合和预测误差较小,与原始位移数据的相关性较好,能够更真实地反映边坡变形规律,具有较强的工程适应性。同时组合模型的提出与发展也促进了单一模型的优化改进,为解决实际工程问题提供了良好的思路。     

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

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

基于Shearlet变换的非局部均值地震噪声压制

在地震勘探中,由于野外地震数据采集环境及仪器性能本身的限制,采集到地震信号中不可避免地会混入较强的噪声,极大影响后续处理、解释工作。而近几年,多尺度几何分析方法以其独特优势成为压制噪声的研究热点,本文提出在Shearlet域中引入非局部均值算法对地震噪声进行压制,该算法首先对地震信号进行非下采样Shearlet变换,然后采用非局部均值法对分解后系数子集进一步处理,并采用8个Sobel算子近似表示全方向结构,对权重函数进行改进,最后对系数进行Shearlet反变换,得到去噪后的地震信号。实验结果表明相比于传统非局部均值法,该联合算法能有效地压制随机噪声,同时对弱同相轴具有更好的保护作用,在地震资料处理中具有良好的实用性。     

基于逆算子估计的AVO反演方法研究

传统反演算法以优化算法为主,而基于逆算子估计的AVO反演算法则利用了直接求逆的思路.算法的关键在于寻找存在逆函数的子域,进而可以在子域内直接求逆,这种解决反问题的思路不同于一般的优化类算法所采用的直接搜索解的方式,具有更高的效率.AVO反演利用了振幅随着偏移距的变化特征,反演的精度受到地震资料质量的影响,通过加入L1范数约束以及合理的初始模型有助于提高反演的稳定性以及准确度.模型测算和实际应用表明,基于逆算子估计的AVO反演方法具有较高的精确程度和可靠性.     

基于偏微分方程的模型光滑处理在偏移成像中的应用

本文针对射线类偏移成像当中的速度模型光滑处理问题,借鉴数字图像处理当中的偏微分方程法,基于能量泛函,应用变分方法导出基于速度模型的偏微分方程实现射线类偏移成像当中的速度模型的光滑处理.由于偏微分方程法具有线性叠加特性、模型解的唯一性和局部特征保持性,因此,应用该算法可以实现基于原始速度模型空间结构的模型光滑处理.通过在原始速度模型以及光滑处理后的速度模型上计算速度的空间分布以及地震波走时、射线路径可以得出,偏微分方程法对速度模型的光滑处理能够很好地保持原始模型的空间结构,偏移成像结果也证明了该方法的实用性.

     

基于非稳态匹配的角度域叠前道集去调谐方法

地层调谐效应改变了不同偏移距处反射波的干涉模式,使得远偏移距数据振幅和频率信息发生畸变,将会降低AVO/AVA分析及反演的置信度和分辨率.文中首先给出了调谐作用下反射波合成记录,通过正演模拟说明了调谐对叠前道集的影响.考虑到不同地层厚度下调谐效应的差异性,借助局域Lamoureux窗实现地震数据的自适应分解,利用角度数据之间的差异,构建了角度域地震数据的非平稳匹配目标函数,形成了非平稳匹配去调谐方法,实现了叠前道集的振幅和波形拉伸校正.本文方法在实际应用中取得了较好的效果,能够有效地改善大角度地震数据品质,为储层预测与流体识别奠定了数据基础.