福建省连续重力观测资料分析

借助Matlab平台,采用功率谱密度估计方法和傅里叶谱,分析2012年4月11日苏门答腊8.6级地震引起的地球自由振荡和2012年第5号台风"泰利"的卓越频率等特征。结果显示:厦门台和福州台能监测到地球自由振荡信号的基频振型;漳州台和福州台监测到的台风的卓越频率约为0.25 Hz。     

大连台连续重力观测质量分析

选取2015-07-01~2016-07-01大连台gPhone重力连续观测数据资料,进行数据降采样、预处理及调和分析,得到该台观测精度及噪声水平。结果表明,大连台连续重力数据在1 min内和1 h内2个周期的噪声水平分别为9.18×10－8 ms－2和0.28×10－8 ms－2,表明该台连续重力观测数据噪声水平较低;重力潮汐主波振幅观测频段低于重力台网中心观测基地超导重力的白噪声水平,说明台站外部观测环境较好。     

重力卫星观测俯冲带特大地震同震重力场变化

利用去相关的DDK3滤波器对得克萨斯空间中心（CSR）发布的GRACE RL06月重力场模型数据进行处理,采用多项式拟合时间序列分析方法成功提取了2004年苏门答腊M_W9.3、2010年智利M_W8.8和2011年日本东北M_W9.0三个地震事件的同震重力变化、大地水准面变化、垂线偏差变化及重力梯度变化。结果表明,3个震例的同震重力变化范围分别为-15.5~6.5 μGal、-9.1~2.1 μGal和-11.1~4.2 μGal;同震大地水准面变化范围为-5.9~0.8 mm、-3.0~0.8 mm 和 -3.2~0.5 mm;垂线偏差NS向变化范围为 -1.2~2.2 mas、-0.9~1.0 mas 和 -1.1~1.4 mas,EW向变化范围为 -1.8~1.0 mas、-0.8~0.8 mas和 -0.7~1.0 mas;重力梯度各分量数值中,rr分量同震变化最大,其次是rθ分量。3次地震的同震信号空间分布均表现为：同震大地水准面和重力变化信号呈非对称两极分布;垂线偏差呈负-正-负或正-负-正三极分布;重力梯度变化信号呈较复杂的多极分布。观测结果与现有的球体位错理论计算结果较一致,可作为重力卫星数据约束震源机制和地球粘度结构的可靠信号来源。     

济南地震台连续重力仪资料质量跟踪与映震能力研究

通过对济南地震台连续重力仪2016~2019年资料质量进行跟踪分析发现,连续重力仪观测精度较高,资料稳定性好,固体潮波形光滑、周期明显,记录到的地震波形态以阶跃和震荡型为主。分析研究其映震能力后得出,济南台连续重力仪对全球7级以上地震均有记录,记录到的最小地震为2.9级,最远震中距为19 114 km,有较强的映震能力。同时,推导出济南台连续重力仪映震能力公式为M=1.875lgD-1.543。     

基于绝对重力基准控制的流动重力观测资料动态平差方法研究

选取滇西重力测网1986～2009年流动重力观测资料及下关、丽江、昆明3个绝对重力基准点多年复测资料作基准控制,采用线性速率模型拟合观测过程中重力场的时变信息,建立绝对重力基准点约束下的动态平差模型,讨论段差观测值及绝对重力基准观测值协方差矩阵的构建、平差模型初值选取及迭代计算等细节问题。经分析,测区内测点的重力年变率范围为(-0.039 5~0.024 1)×10^-5ms^-2/a,精度范围为(0.001~0.029 8) ×10⁵ms^-2/a,与分期静态平差拟合结果的差异约为(0.003~0.020)×10^-5ms^-2/a。在处理全国网和较大面积的区域网多年重复观测的重力资料时,尤其是多网多期重力资料联合处理时,应采用动态平差模型来获得统一时间基准和空间基准的重力变化。     

2022年门源MS6.9地震同震应变观测

以2022年门源M_S6.9地震为例,将17个台站的四分量钻孔应变仪1 sps采样率的同震应变观测值与利用QSSP软件正演计算得到的同震应变模拟值进行拟合,分析二者的相关性。结果表明,17个台站NS向和EW向的P波波形拟合较好,相关系数平均值分别为0.835和0.842,达到高度相关,表明钻孔应变观测数据可较好地反映地震破裂信息,在理论上进一步验证钻孔应变观测数据可用于地震破裂过程的研究。     

福建区域重力台网同震响应特征分析

本文收集整理了福建区域重力台网2015年以来的观测资料,分析各台gPhone型重力仪的同震响应特征,对比并分析震级、震中距、震源深度对同震响应的影响。结果表明,同震响应的最大变形幅度与震级及震中距有较强的相关性,但没有严格的比例关系;面波延迟时间与震中距成正比;浅源地震在福建地区的危害性要大于深源地震;同震持续时间与仪器分辨率有关。     

陆态网络腾冲台连续重力观测质量分析

对陆态网络腾冲台2012-07-01~2013-06-30连续重力观测资料进行处理,获得该台站观测精度及噪声水平。结果表明,腾冲台1~100s噪声水平在9×10^-8ms^-2左右,100s~100min噪声水平在0.3×10^-8ms^-2左右,白噪声水平为0.509×10^-8ms^-2,对O₁波、M₂波等主要潮波振幅因子测量精度达到SG-053超导重力仪的观测量级,其中M₂波振幅因子测量精度最高（为0.000 18）。腾冲台可为川滇地区强震和腾冲火山活动监测、震前重力变化及相关研究提供高精度、可靠的连续重力观测资料。     

不同采样率下四分量钻孔应变观测同震变化特征

利用湟源台四分量钻孔应变观测的分钟值、1 sps、10 sps和100 sps四种不同采样率的观测数据,通过自检分析、同震应变阶分析和频谱分析等方法研究青海玛多M7.4地震同震变化特征。研究结果表明,分钟采样记录的地震波信息缺失严重,用分钟采样数据进行地震波初动、同震变化幅度等研究将会得到信度较低的结论;采样率越高,记录应变地震波信息的能力越强,但100 sps采样和10 sps采样结果相差不大,10 sps采样已能记录到比较全面的应变地震波信息;同震应变阶的变化性质和变化幅度与采样率无关;未来布设四分量钻孔应变仪时,建议应将采样率至少提高至1 sps。     

½̬�����ڳ�̨���������۲���������

????????????2012-07-01~2013-06-30?????????????????д????????????????????????????????????1~100??s????????9×10^-8ms^-2???????100??s~100??min????????0.3×10^-8ms^-2????????????????0.509×10^-8ms^-2??????O₁????M₂??????????????????????????SG-053???????????????????????M₂?????????????????????0.000 18????????????????????????????????????????仯??????о??????????????????????????????     

沈阳连续重力非潮汐变化异常分析

沈阳台定点重力非潮汐变化数据稳定,但自2017-06以来数据曲线连续5个月明显上升,且速率较大,幅度约为5μGal/月,打破了自2015年断电后重归的稳定状态。通过对仪器观测系统、区域地质特征、GPS资料及周边环境等分析发现,GPS垂直分量于同年5月开始下降,GPS垂向变化引起的重力量值-1.714μGal/mm远小于重力变化;抽水的干扰仅产生0.2181μGal的变化,远小于5μGal/月的量值变化。对GPS垂向引起的重力变化进行定量计算分析,并对水井抽水进行正演模型计算,排除抽水和垂向重力变化的影响,重力非潮汐变化量值约为14.7819μGal。研究发现,外界影响变化不足以引起重力的非潮汐异常。综合分析认为,该重力非潮汐变化异常为地震前兆异常,为后续地震预报提供了可靠依据。     

利用EOF分析方法研究康定MS6.3地震前后重力时空分布异常特征

对康定MS6.3地震前后的流动重力观测数据采用正交经验函数(EOF)分析方法进行处理,分析主要模式与资料时段的MS5.0以上地震震中位置的关系,探讨主要模式代表的物理含义。结果表明,川滇地区7期累积重力观测资料中大于200 μGal的数据是偶然误差的概率大于95%。第二模式规格向量和累积重力变化均具有相似的分区特征,2013~2014年6次MS6.0以上地震和第二模式规格向量正负转换带有超过83.3%的对应率,具备重力长期变化模式。第一模式地震前后规格值序列明显的同震信号、空间相似的4象限分布以及2014年4次MS6.0以上地震和第二模式规格向量正负转换带100%的对应率表明,该模式代表重力短期变化模式。康定MS6.3地震前第一模式包含4象限特点,震中位置规格值不断增大,且地震前后变化从负转为正。     

动态与静态平差方法在流动重力数据处理中的对比研究

为研究重力时变因素对流动重力网数据处理结果的影响,利用动态和静态2种平差方法对南北地震带南段流动重力观测数据进行处理,并对2种方法的重力变化结果进行对比分析。结果表明：1）重力时变因素会导致静态平差方法计算的重力变化结果存在误差,对0.5 a和1 a时间尺度下的重力变化具有较为明显的影响,但对2 a以上时间尺度下的重力变化影响较小,因此计算0.5 a和1 a时间尺度下的重力变化时宜采用动态平差方法;2）研究区内九寨沟7.0级、长宁6.0级和漾濞6.4级地震的发震地点与重力变化零值线具有较好的对应关系,重力变化图像可反映3次地震的发震背景。     

一种航空重力测线交叉点搜索新方法

针对目前航空重力测量数据处理中常规搜索方法存在搜索速度慢、正确率低、适用性差的不足,提出以主测线点为搜索中心,采用一定宽度的搜索框对主副测线点进行搜索,将搜索到的主副测线上相邻测点组成线段,运用行列式方法精确求出交叉点及不符值的方法,即滑动窗口搜索法。从搜索时间和正确率两方面对滑动窗口搜索法进行分析,结果表明本文所提方法简单有效,能明显提高测线交叉点搜索效率。     

2020年土耳其MW5.7地震InSAR同震形变特征与破裂滑动分布

为准确认识2020-06-14土耳其MW5.7地震的发震位置、构造特点以及地震危险性,利用D-InSAR技术对Sentinel-1A数据进行处理,基于GACOS进行大气校正,获得视线向（line of sight,LOS）同震形变场。降轨LOS向同震形变场显示,断层北侧抬升,最大抬升形变量约8.87 cm;南侧沉降,最大沉降量约-7.75 cm。以LOS向同震形变为约束,先采用贝叶斯自举优化法反演发震断层几何参数,然后使用有限断层方法反演地震破裂滑动分布。结果显示,断层走向约257.48°±0.65°、倾角约79.69°±0.98°、滑动角约154.2°±3.8°,震中位置为（40.754°E,39.389°N）,破裂区域长度约8 km、宽度约6 km,破裂的最浅埋深约0.8 km、最大埋深约8.9 km,最大滑动量约0.57 m,对应深度约4.278 km。地震释放地震矩约4.54×1017 Nm,对应矩震级MW5.7,与土耳其灾害和紧急情况管理局公布的结果一致。此次土耳其地震受近东西向的潜伏断层控制,以右旋走滑为主兼具少量的逆冲特性。地震造成部分地区的库伦破裂应力增量超过0.1 bar,这些区域未来的地震危险性值得关注。     

LCR-G型相对重力仪周期项参数的确定方法

基于国家重力基准网观测数据,采用正则化方法对LCR-G型相对重力仪周期项参数进行筛选,并分析平差解算的结果。结果表明,当惩罚因子λ取80时,正则化方法筛选相对重力仪周期项的效果良好,正则化参数α值的变化对周期项筛选结果影响甚微;周期项振幅量级最大超过10 μGal,且对重力基准网平差重力值结果的影响达到μGal量级;相对重力仪的周期会随着观测环境或者时间变化,导致解算结果不一致,这是由于周期函数模型与实际周期变化有差异,模型只能近似而不能真正模拟周期的变化。     

扩边尺度对重力异常分层分离处理的影响——以插值切割法为例

以插值切割法为例,对理想模型在地表的重力异常数据进行不同尺度的扩边,并对扩边后的数据进行不同深度层源的切割分离.基于分离后的结果,从不同扩边尺度、不同深度层源的重力异常差异、扩边尺度与重力场源分离误差的角度研究扩边尺度对重力异常分层分离处理的影响.结果表明:1)理想模型数据区域的边界效应与扩边尺度、切割分离深度、异常的...     

钻孔应变观测原位标定方法研究进展

钻孔应变观测的原位标定是制约钻孔应变定量观测理论与技术发展的基础性问题。为推进钻孔应变观测在地球动力学研究中的应用,从标定模型、参考应变信号等方面系统总结钻孔应变观测原位标定方法研究的最新进展,指出钻孔应变观测原位标定中存在的主要问题,并展望今后的发展方向。