2013年芦山Ms 7.0级地震断层参数模型反演

2013年4月20日四川芦山地区发生了Ms 7.0级地震。利用GPS三维同震形变数据获取地表形变场,基于位错模型反演芦山地震的断层几何参数及滑动分布。首先采用多峰值颗粒群算法（multiple peak particle swarm optimization,MPSO）得到断层几何参数,其中断层走向206.47°,倾角44.11°,长度21.94 km,离地表最浅处为7.66 km,最深处为17.84 km。为了反演断层面的精细滑动分布,分析地震所在的龙门山断裂西南段破裂面具有的铲状型特征,将芦山地震破裂面确立为铲状模型,即将断层的倾角预设为上陡下缓,倾角变化范围为21°~50°。结果显示,断层破裂面在不同深度区域出现了两个滑动峰值,其中最大滑动量为0.68 m,深度位于13 km。地震释放的能量为1.47×1019 N·m,对应的矩震级为Mw 6.74,与地震学的研究结果一致。     

利用光学影像相关获取2019年Ridgecrest地震序列同震形变

光学影像覆盖范围广,质量越来越高,可以应用于地表形变监测研究.文中以基于ENVI的COSI-Corr作为数据处理平台,通过光学影像获取地表形变.文中选取2019年加州Ridgecrest地震序列覆盖区域.首先,探究哨兵2号影像误差处理过程中横向条带处理方法,对其进行改进;然后,使用哨兵2号光学影像获取Ridgecres...     

系数矩阵误差对地壳应变参数反演的影响

针对地壳应变参数反演模型中系数矩阵含随机和非随机元素及观测数据存在相关性等情况,以部分变量误差（partial-errors-in-variables,PEIV）模型为基础,采用了地壳应变参数反演的加权总体最小二乘算法,该算法不受系数矩阵和权矩阵结构的限制,能够快速、有效解决系数矩阵含有随机误差的模型问题。结合推导得到的最小二乘改正项公式,对地壳反演模型中坐标点误差对反演参数求解的影响进行了分析。通过对模拟数据和川滇地区的实际数据进行处理,得出系数矩阵误差对地壳应变参数反演的影响主要受GPS站点坐标值量级以及应变参数量级的牵制。     

MATLAB在重力图制作中的应用

重力图是重力测量的主要成果之一,在地质、探矿等领域中有着广泛的用途。MATLAB是数据处理中应用广泛的优秀软件之一,具有编程简单、计算结果可视化等优点。本文在给出重力图相关知识的基础上,用MATLAB来实现重力图的制作,并通过实例显示了用MATLAB制作重力图的可行性和简单性。     

病态总体最小二乘问题的共轭梯度解法

针对测量数据处理中存在的病态问题,该文提出采用共轭梯度方法求解病态总体最小二乘问题。利用Tikhonov正则化的思想构造目标函数,将总体最小二乘问题转化为无约束最优化问题,然后利用共轭梯度法进行求解,避免了过程中的矩阵求逆运算。通过模拟数值算例表明了该方法在解决病态总体最小二乘问题中的有效性,并分析比较了该方法同正则化总体最小二乘方法之间的差异。     

利用稳健WTLS方法进行三维激光扫描标靶球定位

针对标靶球定位缺少有效的确定协因数阵方法的问题,顾及点到平面的距离反映了点与平面的相关性及入射角对点云数据点位精度的影响,将两者推广到三维激光扫描标靶球定位中。以距离、先验入射角确定各点协因数值,并给出观测向量协因数阵及系数矩阵协因数阵,利用稳健加权总体最小二乘方法进行标靶球定位。实例表明,以距离确定协因数阵的稳健加权总体最小二乘方法解算标靶球参数估计比其他方法更精确。     

变异函数模型参数的加权总体最小二乘回归法

顾及距离值的随机误差,提出用加权总体最小二乘回归法估计变异函数模型参数。通过协方差传播律发现,分组后的变异函数值和距离值是不等精度的。给出距离值的定权方法,结合熵权法和点对数法迭代解算模型参数。以幂函数模型为例,模拟数据和实测数据的结果表明,加权总体最小二乘回归法更加合理,参数的估计精度也更高。     

病态乘性误差模型的加权最小二乘正则化迭代解法及精度评定

针对乘性误差模型的病态问题,引入Tikhonov正则化方法,导出了病态乘性误差模型的加权最小二乘正则化解.顾及加权最小二乘正则化法在求解病态乘性误差模型时,参数估值与观测值之间存在复杂的非线性关系,本文利用一种无需求导、通过加权的方式便能够计算非线性函数的均值和均方误差阵的比例对称采样的无迹变换(scaled unscented transformation,SUT)法,对病态乘性误差模型进行精度评定.模拟算例和真实算例结果表明,本文提出的加权最小二乘正则化迭代解法可以有效减弱模型的病态性,基于SUT法的精度评定方法能够得到比已有方法更为合理的精度信息,具有较强的适用性.     

测角网坐标平差模型的Bootstrap参数估计方法

针对线性最小二乘法处理非线性模型产生模型误差的问题,文章将高斯牛顿迭代法引入测角网坐标平差模型中,给出测角网坐标平差模型的高斯牛顿迭代法计算过程.考虑到非线性平差模型的参数估计值是有偏估计,结合Bootstrap重采样方法对参数估值进行改善,提出测角网坐标平差模型的Bootstrap参数估计方法,并给出详细的迭代流程图.针对等精度与不等精度角度观测数据,设计两个测角网案例.实验结果表明,测角网坐标平差模型的高斯牛顿迭代解法能够削弱线性近似带来的模型误差影响,其参数估值优于传统的线性近似方法;而测角网坐标平差模型的Bootstrap参数估计方法比高斯牛顿迭代解法在提高测角网坐标平差参数估值质量方面更加有效.实验证明将高斯牛顿迭代解法应用于测角网坐标平差模型的必要性与实用性,也证明将Bootstrap重采样参数估计方法与高斯牛顿迭代解法结合并用于测角网坐标平差的可行性与有效性.     

火山Mogi模型形变反演的虚拟观测法及方差分量估计

针对Mogi模型垂直位移与水平位移联合反演中的病态问题,改进火山形变总体最小二乘(Total Least Squares,TLS)联合反演的虚拟观测法,并使用方差分量估计(Variance Components Estimation,VCE)方法确定病态问题的正则化参数.将附有先验信息的参数作为观测方程,与垂直位移和水平位移的观测方程联合解算,推导了三类观测方程联合反演的求解公式及基于总体最小二乘方差分量估计确定正则化参数的表达式,给出了算法的迭代流程.通过算例实验,研究了总体最小二乘联合反演的虚拟观测法在火山Mogi模型形变反演中的应用;算例结果表明,三类数据的联合平差及方差分量估计方法可以确定权比因子并得到修正后的压力源参数,具有一定的实际参考价值.