几种不同矩阵算法的遥感图像几何精纠正效果比较

分别采用Gauss、Householder及Givens三种变换形式进行遥感图像几何精纠正试验，并将试验结果与商业遥感软件PCI的纠正结果进行比较。结果表明，Givens算法在一至五次多项式纠正中，其总体精度优于其它算法。     

加权总体最小二乘在重心化Bursa模型中的应用

研究了基于加权总体最小二乘的重心化布尔沙模型的坐标转换算法,针对坐标转换中原始坐标和目标坐标均存在误差的问题,根据误差的影响程度不同而给予不同的权值,利用加权总体最小二乘方法求解转换参数。坐标重心化的方法可以解决布尔沙模型在局部地区容易出现病态的问题,将两种方法结合可以很好地提高坐标转换的精度,通过实例验证了该方法的优越性。     

中国地震台网监测能力评估和台站检测能力评分(2008—2015年)

针对中国地震台网"十五"项目建成后的地震监测能力科学评估的需求,为进一步优化台网布局、提升边疆海域等重点地区监测能力,本研究利用"基于概率的完整性震级"（PMC）方法,对中国地震台网1001个台站以及2008-10-01-2015-09-17期间实际产出的地震观测资料进行了研究,分析了指定震级档下的检测概率P_E和最小完整性震级M_P的分布.除台网整体监测能力分布外,为直观地用单分值表述逐个台站的地震检测能力,本文发展了基于等振幅曲线的"地震检测能力评分表",给出了国家台和区域台每个台站的地震检测能力评分统计特征和空间分布特征.此外,研究中还采用设定"最佳"地震监测能力目标函数的方式,模拟了通过改进观测条件可获得的地震台网监测能力提升的理论结果.研究结果表明,我国华北和东南沿海等东部地区地震监测能力较高,西部尤其是青藏高原南部地区M_p仅约为4.5,近海海域M_p仅约为3.5;从单个地震台站的运行效益角度,台网运行水平和地震观测资料的分析程度对台站的实际的地震检测能力影响显著,新疆等部分台站稀疏地区地震检测能力较高,而中等台站密度的贵州等部分区域相对较低;国家台的地震检测能力评分D_score系统优于区域台,新疆等西部边疆地区,以及福建等东南沿海地区的D_score明显高于台站密集的东部地区;模拟结果显示,在现有台站布局条件下,通过台站优化改造和提升运维管理水平,可显著提升对内蒙古西部、四川西部、甘肃-青海的北部交界地区、鄂尔多斯地块内部、贵州大部分地区,以及我国近海海域、朝鲜半岛北部和中南半岛北部地区的地震监测能力.     

地壳形变分析的抗差最小二乘配置迭代解法

利用最小二乘配置进行地壳形变分析,其结果的合理性关键在于经验协方差函数的拟合.考虑到观测数据存在粗差的情况,提出基于观测值中位数初值的抗差最小二乘配置方法和基于中位参数法的抗差最小二乘配置方法.两种方法首先分别利用观测值中位数给出观测值初始权阵以及利用中位参数法给出最小二乘配置初始解,然后均在给定协方差函数参数初始值的情况下,应用合适的等价权进行抗差估计并通过迭代计算,最终获得稳健的协方差函数参数估值及最小二乘配置解.利用本文提出的两种方法以及传统方法分别对庐山地震的GPS垂直位移数据和意大利L'Aquila地震的InSAR同震位移数据进行处理分析.结果表明:相对传统方法,基于观测值中位数初值的抗差最小二乘配置方法效果更好,更具稳健性.     

基于热点参考系的板块绝对运动模型

本文在最新一代板块相对运动模型(MORVEL模型)的基础上,通过最小二乘反演观测热点方向数据建立板块绝对运动模型.经过对前人用于约束板块绝对运动模型的所有热点数据与MORVEL模型的一致性作系统的统计学检验,发现存在离群数据.相应地,提出两种系统剔除离群数据的新方法:逐一剔除法和全局搜索法.结合对总残差的卡方检验和对残差频率分布的正态检验,逐步筛选最优模型,最终得到一个基于热点参考系的板块绝对运动新模型(T87模型).该模型能够合理拟合全球分布的87个热点方向数据,而模型预测的板块绝对运动速率比观测热点火山的迁移速率系统偏小,最大偏差达到4 cm·a~(-1).这样的偏差可由观测热点速率的系统误差或地幔返回流造成的热点间运动解释.不论是哪种解释成立,本文结果表明热点方向数据能够独立、有效地定义全球热点参考系.这样定义的热点参考系可方便地应用于板块绝对运动、地幔对流及真极移的研究.     

基于非结构化网格的最小二乘逆时偏移

相对于传统的逆时偏移,最小二乘逆时偏移具有更高的成像质量,这种改善是通过迭代反演来获得的,另外其精度与效率依赖于求解波动方程算法的精度与效率.本文给出了基于非结构化网格的最小二乘逆时偏移的方法,该方法能够充分结合最小二乘算法与非结构化网格精细刻画地下界面以及随速度自适应剖分的优点;并采用带补偿的拉普拉斯滤波算法,来消除梯度计算中的低频噪声,从而加速目标函数的收敛速度.通过简单倾斜模型以及复杂Marmousi模型测试,显示了该方法的有效性和潜力.     

基于平面波加速的VTI介质最小二乘逆时偏移

地震各向异性集中表现为速度各向异性,势必影响地震波运动学特征.传统声波逆时偏移(RTM)和最小二乘逆时偏移(LSRTM)没有考虑介质各向异性特征,导致反射波不能正确归位、同相轴出现扭曲及寻优速度慢或不收敛等,VTI介质逆时偏移(VTI-RTM)矫正了声波成像的不足,但仍存在低频干扰严重、中深部成像不佳、振幅保持差等缺陷.为此,本文首先实现了VTI介质最小二乘逆时偏移(VTI-LSRTM)方法,为了节省I/O及内存需求并提高效率,进一步引入平面波编码技术,提出了一种基于平面波加速的VTI介质最小二乘逆时偏移(VTI-PLSRTM)策略.在此基础上开展了简单模型及复杂Marmousi模型成像试验,并与标准逆时偏移剖面对比表明:本方法能够校正各向异性造成的相位畸变,且在迭代中自动压制串扰及低频噪声、补偿中深部能量,是一种兼具质量与效率的保幅成像策略;对速度误差的敏感性测试说明该方法需要相对正确的偏移速度及Thomsen参数模型.     

基于不等式约束的三维电阻率探测混合反演方法

三维电阻率探测的线性反演和非线性反演中均存在着多解性的固有难题.电阻率线性反演方法的效率较高,但反演结果对初始模型的依赖性较强,易陷入局部极小;而非线性反演方法不依赖初始模型,但搜索效率极低,尚未见到关于三维电阻率非线性反演的文献.针对上述问题,融合线性与非线性反演方法的互补优势,提出了最小二乘法(线性方法)与改进遗传算法(非线性方法)相结合的混合反演方法的概念和思想.首先,提出了将介质电阻率变化范围作为不等式约束引入反演方程的思路,以实现压制多解性、提高可靠性的目标.提出了宽松不等式约束和基于钻孔推断的局部严格不等式约束的获取及定义方法.在此基础上,分别提出了基于不等式约束的最小二乘线性反演方法和遗传算法非线性反演方法.其次,对于遗传算法在变异搜索方向控制、初始群体产生等方面进行了改进,优化了其搜索方向和初始群体多样性.然后,提出了混合反演方法及其实现方案,利用改进遗传算法进行第一阶段反演,发挥其对初始模型的依赖程度低的优势,搜索到最优解附近的空间,输出当前最优个体;利用最小二乘法进行第二阶段反演,将遗传算法得到的当前最优个体作为初始模型,在最优解附近空间执行高效率的局部线性搜索,最终实现地电结构的三维成像.最后,开展了合成数据与实际工程算例验证,与传统最小二乘方法进行了对比,发现混合反演方法在压制多解性、摆脱初始模型依赖和提高反演效果方面有较好效果.     

GA辅助NLS的GNSS-IR土壤湿度反演方法北大核心CSCD

采用具有阻尼因子的函数模型,使用遗传算法(genetic algorithm,GA)辅助非线性最小二乘(nonlinear least squares,NLS)方法对相位参数进行求解。结果表明:1)相较于标准余弦函数模型,该方法的反演相位与土壤湿度的相关系数有较为明显的提升,反演结果也更加稳定,在5°~15°、5°~20°、5°~25°三个高度角范围内的相关系数均大于0.68,不同高度角之间的相关系数差值小于0.07;2)反演精度有不同程度提高,R 2提高5.72%~76.06%,RMSE减小6.12%~24.24%,MAE减小2.7%~28.3%,将该方案所求相位用于多星线性回归模型后平均RMSE减小10%。     

基于三维层状介质模型地震层析成像正演网格研究

地震层析成像技术是研究地层速度结构的一个重要手段,该技术是利用震源到接收点的地震波走时或波形的观测数据,结合建立的数学物理模型来重建地层速度结构。而层析成像正演模拟的精度将直接影响反演出的速度与地层真实速度的拟合度,因此找出一种精度较高的层析成像正演模拟算法是非常重要的。本文针对三维层状介质,通过网格化模型,采用最小走时射线追踪方法展开层析成像正演研究,分别采用三角形网格、矩形网格和六边形网格进行模型参数化,并对比分析各自正演结果的精度,总结出针对三维层状介质模型的最佳网格划分方式。