精密重力测量中相对重力仪格值系数的贝叶斯估计方法

相对重力仪的格值系数随时间会发生微小的变化,是影响精密重力测量精度的重要因素。通常需定期对相对重力仪进行专门的基线标定来评估仪器格值系数的变化。本文提出了一种利用重力观测数据进行格值系数评估的新方法,原理是利用测网中已知的多个绝对重力基准点作为先验约束,同时考虑仪器的非线性漂移变化,将格值系数作为超参数,基于贝叶斯原理和赤池贝叶斯信息准则(ABIC)估计最优值。通过对模拟数据的测试,该方法在高斯噪声和仪器非线性漂移等不确定性存在的情况下,可以获得格值系数的准确估计结果。对实测重力数据的测试表明:估计的格值系数与测量前在基线场标定的格值系数差值在5×10^-5以内;而且相较于采用标定不准确的格值系数,该方法可以获得与绝对重力测量结果差异更小的平差重力值。本文研究结果为有效提高精密重力测量的效率和精度提供了方法保障。     

重庆地震应急资料产出系统研发与应用

地震应急资料具有较强的时效性、专业性和多学科融合的特点。中强地震发生后启动应急响应,工作人员需在短时间内产出多学科、跨领域、多手段的专业地震会商资料。以地震会商系统为平台,利用其灵活开发、自动计算和智能分析的优势,建立地震应急资料产出系统（包括会商资料处理模块、应急图件产出模块和震情报告生成模块）,通过地震触发,根据预设流程和报告模板,自动产出专业报告,为震情形势研判和震后应急救援提供专业可靠的应急资料,最大限度地发挥应急救援和灾损评估的效能。     

满都拉—卡楚加的长剖面地壳介质密度分布及深部结构特征探榷

通过对南起中国内蒙古满都拉、向北经蒙古国、直抵北端俄罗斯贝加尔湖北卡楚加全长1320km长重力剖面的布格重力异常数据进行处理和解析, 依其构建了沿剖面的二维地壳密度结构模型, 并详细分析了沿剖面的地壳底界面(Moho界面)展布的深部构造特征, 及剖面辖域内的5个次级局部构造单元的重力异常场、地壳密度结构、界面起伏及断裂构造分布的特征.期望此研究结果能为该剖面跨越的中亚造山带东部地区的地壳结构、各次级构造单元的界域与关联方面的进一步研究提供相关的重力场依据.

     

塔北隆起X型走滑断裂成因机制的新解释

塔里木盆地塔北隆起发育两组呈小角度相交(40°)的透入性X型走滑断裂,分别沿着NNE走向和NNW延伸。在对塔北哈拉哈塘地区三维地震资料解释的基础上,对走滑断层的几何展布特征以及断层的剖面变形特征进行研究;同时重点解析了RP6断裂和HA13断裂,分析比较NNW向与NNE向断层的变形及发育特征差异;结合盆地重磁资料以及周缘造山带的活动特征,对塔北隆起小角度的X型走滑断层的发育机制以及演化进行了分析。研究表明,塔北隆起走滑断层在垂向上具有明显的分层变形特征,分为三个构造层:震旦系-中寒武统下构造层(TH₃界面以下)、上寒武统-中奥陶统中构造层(TH₃-TO₃t界面)和上奥陶统-石炭系上构造层(TO₃t-TP界面)。断层在下构造层和中构造层中整体处于压扭环境,多发育正花状构造;上构造层中断层主要发育负花状构造或正断层,整体处于张扭环境。两组断裂比较,NNW向断裂活动性强,在各构造层中均有显著的断裂特征发育,垂向连通性强,发育先存基底断裂,而NE向断层主要发育在中构造层,在下构造层和上构造层中断层发育不明显。活动性分析表明,断层的形成与演化具有多期性,走滑断层的形成经历了三期主要活动:中寒武世末、中晚奥陶世和志留纪-石炭纪。塔北隆起X型走滑断裂的形成受到了NNW向基底断裂和薄弱带的控制,NNW向先存基底断裂带或薄弱带优先发育走滑断裂,基底断裂与主挤压应力方向的夹角小于45°-Φ/2,NNE断层的发育受NNW向先存断裂限制,最终形成小角度相交的X型断裂。     

GRACE卫星观测到的与汶川Ms8.0地震有关的重力变化

利用GRACE卫星重力资料,计算了中国大陆及周边的卫星重力时变场和地表密度变化分布,获取了具有代表性的点位区域的每月重力变化时间序列.同时获得了WUSH、LHAS、KUNM、LUZH站相对于区域参考框架的GPS位移时间序列.卫星重力观测结果显示喜马拉雅弧形带的重力在2004年苏门答腊M_w9.3地震后快速下降, 2006~2008年尤为明显,西域地块西北边界带上震后重力下降也较为显著;而沿青藏高原北至东边界2007年出现明显的重力上升沿构造边界的弧形分布,且2008年南北地震带中南段重力上升变化显著.这些苏门答腊地震后的重力变化趋势到汶川地震发生后才开始改变.GPS位移结果显示四个台站均记录到苏门答腊大地震的同震信号,震后WUSH、LHAS、KUNM站水平位移向量出现明显的运动趋势改变,且一直持续到2008年汶川M_s8.0地震的发生.GRACE卫星揭示的青藏高原及周边地表质量的变化为解释汶川地震的动力机制提供了新的观测途径和资料.本文结合区域构造运动的特点和GPS位移,对GRACE观测的时变重力场特征及汶川地震的动力机制进行了初步解释和讨论.     

InSAR数据约束的2020年西藏定日Mw5.7地震源模型及构造意义

2020年3月20日,处在北喜马拉雅断裂与申扎一定结断裂交会区的西藏定日发生Mw5.7地震,此次地震发震断层源模型的研究对于认识该地区的复杂地质构造具有重要意义.本文利用升、降轨Sentinel-1A星载SAR数据,通过合成孔径雷达干涉测量(InSAR)处理获得了地震引起的地表位移场,其中,雷达视线方向(LOS)最大地...     

海南东寨港海草资源分布特征及影响因素

东寨港海草种类主要有2科2属3种,分别是单脉二药草Halodule uninervis、卵叶喜盐草Halophila ovalis以及贝克喜盐草Halophila beccarii,卵叶喜盐草与贝克喜盐草为东寨港优势种。东寨港海草资源面积约264.00 hm^2,平均盖度在52.13%,平均密度9901.00株/m^2,平均生物量37.46 g/m^2。东寨港内海草资源处于退化状态,主要受底质类型、水环境、潮汐汊道、渔业活动、陆源污染、非法渔业以及自然因素等影响。通过对东寨港海草分布特征及变化分析,提高公众对潟湖海草保护意识,为海南岛潟湖海草床的保护与修复提供科学依据。     

时变重力场球面模型反演算法和模拟实验

时变重力场是研究地球内部介质物性变化的重要手段。本文提出了一种适用于地面流动重力测量获得的时变重力信号的场源反演方法,该方法采用球坐标系下的六面体单元来模拟场源介质,适合大尺度地震流动重力测量数据的等效源模型构建。通过引入重力时变信号的一阶光滑先验条件,压制了时变重力信号中的短周期高频分量,可用于提取与地震孕育相关的长周期信号。通过理论和模型实验证明了本文算法的可靠性和稳定性,并使用南北地震带南段2014—2017年的流动重力实测数据进行了反演解释,获得了地壳内部等效场源的视密度时变信号,变化量级在正常地壳密度的±0.7‰之间,其空间形态受川滇菱形块体边界控制。研究成果可用于时变重力场模型解释和深部场源特征提取,可为地震重力前兆信号分析和相关研究提供完备的方法保障。     

基于等效源的重力异常空间导数计算方法及其在川滇地区的应用

在综合分析前人研究成果的基础上,采用基于贝叶斯准则的等效源模型滤波增强方法来识别地质体边界。通过理论模型实验讨论了在相同的噪声水平情况下,不同算法对于计算结果的影响,模型实验结果表明,与直接计算和传统高斯滤波增强技术相比,新算法可以有效消除传统算法中普遍存在的对于噪声的过度敏感问题,所得到的地质体边界形态失真更小。随后,利用此方法对川滇地区的重力数据进行了处理,并对该地区的重力水平总梯度(THD)、重力垂直梯度模(AVG)、重力解析信号模量(AAS)、重力水平梯度垂向导数(VGHD)应用效果进行了对比分析。研究结果提供了川滇地区的构造分区信息,对该地区构造演化、边界定位及动力学过程分析具有重要的指导意义。     

应用Slepian局部谱方法解算中国大陆重力场球谐模型

根据经典的球谐函数方法,为满足正交化要求,观测数据需要覆盖整个球面,而对于地表局部测量数据,则无法应用球谐方法解算重力场模型。针对此问题,采用Slepian局部谱分析方法解算中国大陆范围内的实测重力场变化数据,并以GOCE卫星球谐函数解作为已知模型,评估由于实际陆地重力测点的非均匀分布对球谐函数解的误差影响。通过计算多个阶次中国大陆局部范围的Slepian基函数分布;采用GOCE卫星获得重力场模型的前72阶球谐系数作为已知结果,评价实际测点非均匀分布的解算有效性,并针对中国大陆地区采用Slepian基函数进行解算,通过模型对比选择最优截段项数;针对2005—2008年中国大陆地区流动重力测量获得的重力场变化信号进行解算,获得了72阶重力场变化模型。