利用VLBI观测量对月球天平动参数解算及着陆器定位和测速的改进

在嫦娥三号探月工程中,月球天平动参数是影响其着陆器定位精度的重要参数,研究月球天平动参数估计方法对于中国探月工程的发展具有重要意义。目前,月球激光测距（lunar laser ranging,LLR）是解算月球天平动参数的常用方法,该方法通过数值拟合将拟合系数放在历表文件中供用户读取。首先,基于差分甚长基线干涉测量技术（very long baseline interferometry, VLBI）解算月球物理天平动参数,其仿真结果表明,与DE421星历的插值结果对比,其欧拉角$\Omega$、$i$和$\mu$的改进值分别补偿至-0.692 4″、0.009 6″和-0.009 7″;然后,基于补偿后的月球物理天平动参数求解着陆器坐标和速度,结果显示,相比于补偿前, 着陆器在X、Y、Z方向的定位精度分别提高了24.204 m、0.405 m、1.996 m,速度误差的估计精度也分别提高了0.010 6 m/s、0.013 5 m/s、0.007 2 m/s。上述研究结果可为未来月球天平动参数解算的相关研究提供参考。     

全球降水观测计划多卫星联合反演降水产品的极端降水监测潜力研究

基于全国自动站与CMORPH融合降水数据集,综合评估全球降水计划（GPM）下的IMERG与GSMaP的4套纯卫星降水数据,从极端降水指标、极端降水事件探测能力、对于不同历时的极端降水事件精度评估3个方面研究纯卫星降水产品监测极端降水的能力。依据极端降水的阈值将研究区分为3个区域进行研究,结果表明：（1）在RX1指标中,IMERG、GSMaP数据均在受季风影响的复杂地形区域呈现出明显的高估状态,在其余地区存在不同程度的低估现象;在R95pTOT指标中,卫星数据均表现较好,与地面参考数据具有较高的相关性。（2）在极端降水事件探测方面,4套纯卫星产品在东北地区的表现均优于其他区域;GSMaP表现优于IMERG数据,具有较低的误报率,但对极端降水的反演精度较低。（3）在不同历时的极端降水事件精度评估中,IMERG、GSMaP卫星降水产品在历时较长的极端降水事件中的表现较好,具有更高的精度;对于极端日降水事件,卫星降水产品在高雨强下的误差显著性十分明显,远远高于复杂地形对卫星降水反演精度的影响,导致卫星降水在复杂地形区域（Ⅲ区）的表现优于其他区域。总体上,IMERG产品在研究区对极端降水的监测能力优于GSMaP产品,其中又以IMERG_Late表现最佳;4套卫星降水产品均能表现出研究区的极端降水区域特征,但在研究区大部分区域呈现出低估状态,卫星降水产品对于雨强的误差订正仍是未来极端降水反演的重点与难点之一。     

被动微波遥感雪深反演混合像元问题研究EI北大核心CSCD

积雪深度是反映地表积雪量变化的重要因子,是水文模型和气候模式中的重要参数之一。被动微波遥感以其穿云透雾、对雪层信息敏感等特点,被广泛应用于雪深的反演研究中。被动微波传感器的低空间分辨率(数千米至数十千米)及地表覆盖的复杂性,使星载被动微波遥感影像中的混合像元现象十分突出,严重制约了被动微波雪深遥感监测的应用和发展。目前,被动微波遥感雪深反演中的混合像元问题研究仍存在着诸多挑战:①理论支撑不足,对被动微波混合像元亮温响应特征及影响机制的研究尚不充分;②针对混合像元问题所发展的被动微波雪深反演算法,对地表异质性特征考虑不足。     

利用TSVD参数估值变化特性确定算法截断参数

TSVD是大地测量病态问题解算的常用有效方法。影响TSVD解算效果的关键因素是截断参数,现有截断参数确定方法可提供有效的截断参数,但仍难以给出最优截断参数。以均方误差最小为准则确定截断参数是一种理论依据较充分的截断参数确定方法,但均方误差计算所需的模型参数真值在实际应用中无法获得,导致该方法难以给出理论最优截断参数。鉴于此,本文研究了基于均方误差影响下(方差与偏差联合影响)参数估值变化特性的TSVD截断参数确定方法。通过TSVD依次截掉小奇异值,获得奇异值截掉前后的方差与参数估值变化,利用两者变化分析确定偏差影响,避免依赖参数真值计算偏差,从而确定出均方误差最小理论下的截断参数。数值与应用试验结果表明,本文方法确定的截断参数可有效改善TSVD解算效果,是一种行之有效的截断参数确定方法。     

加乘性混合误差模型精度评定的SUT法EI北大核心CSCD

已有加乘性混合误差模型参数估计方法能达到二阶精度,但精度评定方法只能达到一阶精度,若通过传统泰勒级数展开近似函数法来获取参数估值的二阶精度信息,由于加乘性混合误差模型中参数估值与观测值为一个复杂的非线性关系,必然需要复杂的求导运算。针对该问题,本文使用一种无须求导、无须了解非线性函数构成的比例无迹变换(scaled unscented transformation,SUT)法来计算参数估值的二阶精度信息。通过算例分析表明,利用SUT法求解加乘性混合误差模型能够有效避免复杂的求导运算,所求得的参数估值及其协方差阵均能达到二阶精度,从而验证了本文方法的可行性和优势。     

Re-Os同位素在沉积地层精确定年及古环境反演中的应用进展

Re-Os同位素亲有机质的性质,使得富有机质沉积岩在沉积的过程中能够吸附富集海水中的Re、Os,沉积岩的沉积压实过程也是其中Re-Os同位素体系封闭计时的过程,沉积岩Re-Os同位素等时线年龄代表地层沉积时代,Os同位素初始比值187Os/188Os反映沉积时海水的Os同位素比值。Re-Os同位素体系在富有机质沉积岩中的成功应用,能够直接确定地层沉积时代,从而对地层界线进行直接厘定,并且能够对一些沉积矿床形成时代、冰川事件发生时代进行厘定和限制。通过沉积岩Re-Os同位素特征,可以对古环境进行反演,有助于了解全球大气海洋的演化,气候的变化,对研究生物灭绝等重大地质事件发生的时限和机制以及金属矿床的成矿物质来源具有重要意义。文章阐述了Re-Os同位素在富有机质沉积岩中应用的原理,并列举了Re-Os同位素在沉积地层精确定年及古环境反演中应用实例,说明了富有机质沉积岩Re-Os同位素分析精度及其影响因素,以使Re-Os同位素体系作为一种强有力的工具,解决更多的沉积地层时代问题以及更好地对古环境进行反演。     

云南省江城地区盐盆地重力勘探与深部找钾启示

云南省江城勐野井钾盐矿是目前我国发现的唯一固体钾盐矿。该地区地质构造复杂,近几年通过野外地质、地球物理及钻井勘探 ,认为钾盐来自深部侏罗系,并沿断裂构造带挤压塑流到表层,可能在深部还存在大面积的“盐源”。为实现找矿突破,亟须了解控制盐矿形成的构造及围岩的三维结构特征。通过对江城地区高精度重力数据位场分离、边缘检测、3D物性反演的地球物理处理与解释,获得了勐野井地区江城湖盆的基底变化、断裂展布等地质信息,确定了研究区构造格架、地下密度体的三维分布与岩盐矿的关系,推断在勐野井矿区西北深部可能还有侏罗系盐岩甚至钾盐的存在,该地区将是未来找矿工作的重点。研究成果为江城地区下一步钾盐矿勘查提供了线索和依据。     

地形对不同反演的影响

[摘 要]根据观测数据和模型参数之间的关系,将求误差函数的最优问题的方法通常分为线性反演和非线性反演; 线性反演要求有相当准确的背景介质模型的条件,这和地球物理问题中的实际情况相矛盾。特殊情况下,线性反演甚至出现发散的情况。而非线性反演不用求目标函数的偏导数,全局随机搜索模型空间的最优解的办法,避免了问题的线性化,在山区复杂地形条件下,通过对比论证,非线性反演能够获取更合理的真实地电模型,从而说明由于非线性反演适应性较强、精度较高已被广泛的应用于地球物理的各个领域。     

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

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

Mohr-Coulomb准则下基于滑动面深度的边坡参数反分析方法

利用上限定理与强度折减法,提出了一种新的边坡反分析方法。该方法基于边坡的滑动面形状,不需要位移、应力等监测数据。分析边坡的临界滑动面位置与强度参数的关系,可以发现,当边坡的几何形状、土体重度以及孔隙水压力分布确定,则临界滑动面的形状只与土体黏聚力c、内摩擦角φ的关系式c/tanφ有关,并在此基础上引入边坡的无量纲系数c,λcφ。利用极限分析上限定理与强度折减法,得到了边坡临界滑动面方程,提出了其深度的计算公式,建立了边坡c,λcφ值与滑动面深度的关系。通过已有边坡的反分析实例,验证了方法的正确性和有效性。