附合重力网的严密平差

本文首先证明了传统的强制符合条件平差法已不适用于现代高精度附合重力网平差,而应根据相关分组平差原理将控制网起始数据作为具有先验精度的虚拟观测值与附合网观测值一起进行平差,其结果与对控制网和附合网实行不等精度一并平差等价.按长段差法组成相对观测方程能获得较高的重力值和仪器格值线性项求定精度.本文推导了有关数学模型,并通过算例对这种平差方法的结果和精度进行分析讨论.     

多期复测水准网和重力网的动态联合平差

在利用重力测量结果研究地壳垂直运动时，重力－高程因子的确定是一个重要研究的问题。本文提出利用多期复测水准网与重力网动态联合平差的方法，同时求解地壳垂直运动速率和重力－高程因子。利用此方法对北京地区１９８３－１９９１年的水准和重力复测结果进行了计算，结果表明：两类不同观测值的权比配置是恰当的，平差模型是正确的，重力－高程因子的分布符合一般规律。     

重复绝对重力测量观测的滇西地区和拉萨点的重力变化及其意义

根据近年来在滇西地区和西藏拉萨绝对重力点上的一些新的绝对重力重复观测结果 ,对1990年以来在这些点上绝对重力重复观测结果进行了分析. 滇西地区的5个绝对重力观测 点的近10年的重复观测结果表明,大部分观测时段没有出现与地球动力学事件有关的重力变 化,只有丽江和洱源2个绝对重力观测点的观测结果显示丽江地震前后有变化. 为了研究该 重力变化的原因,本文正演模拟了出现重力变化期间丽江地震同震位错引起的重力变化,模 拟结果与实际观测结果有较好的一致性. 西藏拉萨近10年以来的重复观测结果给出了拉萨点 的绝对重力值以-1.82±0.9×10-8m·s-2·a-1速率下降,这从重力学 的角度反映出青藏高原的隆升. 对拉萨点的重力变化机制进行了探讨,根据印度板块与欧亚 大陆俯冲模型计算的拉萨点重力变化速率与观测到的重力变化速率较一致,表明现今拉萨地 区重力变化是由于印度板块与欧亚大陆俯冲所引起. 根据印度板块向欧亚大陆俯冲的位错模 型计算的拉萨点的重力变化与地表垂直位移的关系,将重力变化转换为拉萨点的隆升速率为 8.7mm/a.      

磁法和重力测量基点网平差计算方法探讨

通过公式推导和举例说明,介绍一种高精度磁法测量和重力测量基点网平差的计算方法。此方法简便快捷,并且能更好地应用计算机进行平差计算。     

青藏高原GRACE卫星重力长期变化

本文采用最新的GRACE(Gravity Recovery and Climate Experiment)(RL05)数据,通过水文模型(Global Land Data Assimilation System, GLDAS与Climate Prediction Center, CPC)扣除土壤水及雪水的影响,利用Paulson提供的冰川模型结果扣除GIA(Glacial Isostatic Adjustment)的影响,采用尺度因子的方法减少数据处理过程中误差的影响,最终基于最小二乘计算方法得到2003—2013中国及周边地区长期性重力异常变化情况.结果发现青藏高原有较为明显的重力上升信号,我们认为该信号可能由印度板块俯冲欧亚板块导致青藏高原地壳增厚所引起.接着依据GPS观测结果和艾黎均衡假说构建了地壳形变模型并通过直立长方体模型予以正演模拟分析.以班公湖—怒江断裂带为界将青藏高原划分为南北两大区块,结果显示青藏高原重力异常大致以0.2 μGal·a-1的速率在递增,小于GRACE得到的0.3±0.08 μGal·a-1的增长速率(对应于地壳增厚速率约3 mm·a-1),剩余未解释部分可能与湖水、冰川因素、冻土因素等有关.该结果对于认识青藏高原隆升动力学有一定参考意义.     

基于欧拉反褶积方法计算辽宁地区重力变化场源特征

以2012—2018年辽宁地区的流动重力观测资料为基础,结合实际地下地质体构造特征构建理论模型,应用欧拉反褶积法对重力变化场源深度及空间分布特征进行反演和解释。结果表明:在辽宁西部辽蒙交界和海城—岫岩断裂带附近场源位置集中,地下介质活动增强;欧拉反褶积方法适用于对流动重力观测资料的重力场反演及以场至源的定量研究应用。     

利用小波分析重力的长期变化

运用小波滤波方法估算Chandler和周年项的潮汐因子.本文分析了四个台站(Brussels, Boulder, Membach以及Strasbourg)的观测记录,运用合成潮方法得到重力残差后,用Daubechies小波带通滤波器滤波残差,得到256～512 d时间尺度上的序列,根据标准差最小原则确定观测极潮周年和Chandler项的周期,然后利用最小二乘法估算它们的潮汐因子,同时给出未经模型改正的周年重力.由于高阶Daubechies小波构造的滤波器具有良好的频率响应,且能压制信号中的高阶异常成分,使滤波的信号更加光滑,因此计算结果具有更小的均方差,更加可靠.     

滇西动态重力网的分形特征及空间分辨力

介绍了用分形学描述地球物理观测网非均匀分布和联合应用分形特征与Ｓｈａｎｎｏｎ定理确定最佳网格距的方法。研究了滇西动态重力网的分形特征，它的分形维Ｄｆ＝１．２７，最佳网格距ｒ０＝２８ｋｍ。讨论了９１８６年至１９９４年间网区发生的８次Ｍｓ〉５．０级地震前重力场变化的可信性，指出这些变化的波长均大于１００ｋｍ，至少２倍于网 最小空间分辨力，因而是可信的。     

京津唐地区的重力变化

我们重力测量组在1981年7月用三台 LCR 重力仪建立了重力测网.测量精度优于100nm/s2.文章描述了京,津,唐地区重力随时间的变化.结果表明在天津地区增加了约1000nm/s2,天津是沉降中心,每年平均沉降约7——8cm.在其它台站,如堂二里,霸县等地重力值也逐年增加,但其幅度小于天津.重力变化的空间分布由东南向西北逐渐减小.     

改算重力段差的简便方法

重力仪的格值因温度变化或其它原因发生变化后,用原格值计算的重力段差值要进行改算或重新计算。本文推导了重力仪格值变动后,改算重力段差值的简单公式,方法比较简便,能保证所需精度。     

多网重力数据联合平差和系统误差改正

华北地区现行的地震重力监测网由于独立联测, 自成体系, 绝对控制点较少, 获得的重力场动态变化图像在绝对控制较弱的地区会产生畸变, 也不能有效的消除仪器标定系统所引起的测量误差。 本文以华北大网作为控制网, 结合绝对重力点进行联合平差计算, 使各省区地震重力监测网平差精度提高了(0~4)×10^-8 m/s², 平差起算基准相统一, 绝对控制较弱的地区畸变减小。 采用相关系数检验法, 对华北地区各省区地震重力监测网2010年二期重力平差计算结果进行检验, 部分重力点位系统误差高达60×10^-8 m/s², 经过系统误差改正后, 重力场动态变化趋于一致。 经过各项改正后, 华北地区2010年重力场动态差分图像对地震前兆具有更好的反映。     

基于GRACE RL05数据研究智利地震长期震后重力变化特征

利用GRACE卫星重力可对地震引起的大范围重力变化进行观测,并从重力数据中发现主要的变化特征.发生于2010年的M_W8.8智利地震震级较高,可观测到震中附近广泛的同震和震后长期重力变化.本文基于GRACE RL05 Level-2时变重力场数据,对2010年智利地震的同震和震后长期变化进行了计算.对同震变化的计算发现,智利地震引起的同震变化极值达-5 μGal,而本文为减小水文信号的干扰而采用的3年平均的方法可以获得良好的效果.在对震后重力变化的计算中发现,智利地震震后在2011—2016年间的重力变化存在先增大后逐渐衰减的过程.对震后变化的拟合表明,智利地震震中附近有约1μGal的震后重力变化,震后变化的特征时间约1.1年.同时,在智利地震中未出现较明显的两个震后变化阶段（短期、长期）.

     

滇西地区基于单体建筑的震后人员死亡评估方法研究

震后死亡人数是震后应急响应和救援工作的重要基础。为提高震后人员死亡评估精度,从云南县(区)城镇化率、建筑面积和建筑楼层数考虑,提出了一种基于单体建筑的人口分配方法;以建筑的致死性矩阵为基础,使用滇西地区建筑数据和历史地震死亡人数计算本地化修正系数,得到了该地区基于建筑类型的人员死亡评估方法;以2021年云南漾濞6.4级地震为例,对评估结果进行了检验和三维可视化。结果表明：(1)基于单体建筑的人口分布空间比传统的公里格网更加精细;(2)基于单体建筑评估的受灾人口和实际受灾人口之间的相对误差为3.7%,精度较好;(3)基于单体建筑致死性评估的漾濞地震人员死亡评估结果为0～3人,符合实际情况;(4)三维表达能够增加信息表达的直观性和丰富程度。     

辽南地区重力变化与地震活动的研究

本文根据近年来多期流动重力监测资料，对辽南地区重力场时空变化特征进行分析研究，认为该区重力场时空变化与地震活动存在着一定的关系。