2021年青海玛多7 4级地震的同震变形分析

2021年5月22日,青海省玛多县发生了 7.4级地震,该地震发生在巴颜喀拉地块北部边界东昆仑断裂带以南约70 km,属于块体内部断裂带地震.根据中国大陆构造环境监测网络提供资料,距离震源30多公里的玛多台站记录到东西向永久位移约25 cm.同时,InSAR也观测到明显的形变场,升轨和降轨的最大相对形变量分别约1.87 m和2.32 m.为了解释这些大地测量观测数据,本文利用该地震的三个断层滑动模型,基于不同地球模型的地震位错理论,计算同震变形场,并分别与GNSS观测数据和InSAR视线向形变量对比分析,结果显示基于InSAR数据反演的断层滑动模型产生的位移场与球形地球模型的理论计算结果最为吻合.进一步,利用较优断层模型计算2021年青海玛多7.4级地震的理论同震位移、大地水准面、重力和应变等变化,该结果为玛多地震的GNSS和重力观测的解释提供理论参考依据.     

三轴磁通门传感器水平和定向对地磁日变化观测数据的影响

针对多套三轴磁通门传感器在同一观测环境下,日变化观测曲线不是完全重合,存在观测数据的不一致性问题,通过理论计算和地磁台站测试实验,分析仪器的定向、底座水平与地磁各分量观测数据之间的变化规律并给出其定量关系.实验结果表明:若存在定向误差角α,则表现为D分量的观测曲线中包含了一定比例的H分量;在水平面内,若存在磁传感器以H分量为轴在ZOD平面内的旋转夹角φ,则表现为D分量的观测曲线中包含了一定比例的Z分量;若存在磁传感器以D分量为轴在ZOH平面内的旋转夹角θ,则表现为H和Z两分量的记录数据中分别包含了一定比例的对方分量成分.      

蓄水池塘放水对昆明地震台形变观测影响

监控昆明地震台观测山洞附近蓄水池塘的放水过程,掌握放水对水管倾斜、洞体应变观测造成的干扰形态,计算观测数据的卸载响应比,根据干扰类型,选用合理方法处理数据,更好地提升水管倾斜、洞体应变潮汐观测质量;对于未来地震异常落实,可用于类似干扰因素的排除及干扰量的计算。     

地球南北半球的非对称性

依据新的计算分析和空间观测数据,进一步论述了地球南北半球的非对称性. 全球热散失量的计算得出,南半球高出北半球33髎;南半球地幔热散失量是北半球的2倍. 比较南北半球S波速度分布,得出南半球的上地幔为低速、高温,北半球的上地幔为高速、低温. 计算地幔各层的质心位置发现,地球的质心偏于北半球. 计算地球经、纬圈长度的年变化率表明,南半球在扩张,北半球在收缩. 用空间大地测量数据的检测结果证实,南半球处于扩张状态,北半球处于压缩状态. 对地球的非对称性作了初步的动力学解释.     

基于重力卫星数据监测地表质量变化的三维点质量模型法

利用卫星重力测量手段监测全球质量变化取得了巨大成功,本文基于牛顿万有引力定律在三维空间直角坐标系中导出利用重力卫星观测数据监测全球质量变化的三维点质量模型法,该方法可直接利用重力卫星的轨道和星间观测数据或时变重力场模型计算全球质量变化,由于利用卫星观测数据计算地表质量变化的向下延拓过程以及观测数据噪声的影响,需要采用合适的空间约束方程或正则化技术对解算结果进行约束或平滑处理.利用合成全球质量变化模型模拟一个月的GRACE双星轨道和星间距离变率数据计算全球质量变化,对三维点质量模型法进行分析验证,采用零阶Tikhonov正则化技术处理病态问题.结果表明,三维点质量模型法可有效用于重力卫星观测数据监测全球质量变化,为利用重力卫星观测数据监测全球质量变化提供一种可选的途径.     

利用 Yabuki＆Matsu'ura 反演方法计算2011年日本东北地区太平洋海域Mw9.0级地震同震滑动分布

Yabuki & Matsu'ura反演方法是利用ABIC最佳模型参数选取方法和平滑的滑动分布作为约束条件,由形变观测数据计算发震断层滑动分布.本文基于日本列岛同震GPS观测数据和发震断层曲面构造模型,利用Yabuki&Matsu'ura反演方法计算2011年日本东北地区太平洋海域Mw9.0级地震的发震断层同震滑动分布.反演结果表明,断层面上的最大滑动量为35 m,较大滑动分布在浅于30 km的震源中心上部,最大破裂集中在20 km深度的地方,其地震矩约为3.63×1022N·m,对应的矩震级为Mw9.0.模拟结果显示Yabuki&Matsu'ura反演方法更适用于倾角低于40°的断层模型反演.最后,本文基于上述方法获得的发震断层滑动模型,利用地球体位错理论正演计算该地震在中国及其邻区产生的远场形变,正演计算结果基本可以解释由中国GPS陆态网络观测到的同震形变.     

基于Monte Carlo方法的由GPS观测计算地应变率的误差分析

在利用GPS观测计算地应变率中，应变率的误差分布在地球动力学研究中非常重要. 为解决这一问题，文中提出基于Monte Carlo技术，利用GPS观测的速度计算地应变率进行误差分析的方法；以青藏高原地区为例，利用带有误差的GPS观测数据进行了大量、独立的应变率重复计算，然后运用统计方法给出了应变率分量及其误差的分布. 结果显示青藏高原的喜马拉雅及中东部地区地应变率的东西向、南北向及剪应变率分量的误差相对较大，但整个青藏高原及附近地区应变率结果稳定可靠；研究结果进一步表明，由GPS观测计算应变率的误差主要决定于GPS观测速度的误差， 在变形小的区域（如中国东部地区），目前GPS的观测精度还不能完全满足地应变率计算的要求.     

台站地磁观测数据精度提高方法

在当前多种地磁日常数据处理方法基础上,探讨提高地磁观测数据精度的问题.通过使用一种改进的地磁观测基线值计算方法,能有效减小外磁场变化对观测结果的影响,考虑到基线值日变化等因素,使用高阶多项式拟合基线值日变化曲线.利用该方法处理的地磁数据与实测数据进行结果比对后,发现此方法能有效提高地磁观测数据的精度和质量.     

网式大地电磁阻抗张量及其影响因素分析

网式大地电磁（Network-MT,N-MT）法采用长数公里至数十公里的电话线为电极线测量电场,很难形成两条笔直且相互垂直的电极线,因此阻抗张量的计算不如大地电磁法中直接.本文依据阻抗张量的旋转规则提出了一种计算N-MT阻抗张量的简便算法.依据该算法计算了中国东北地区5个N-MT测站的阻抗张量,获得了基于阻抗张量的视电阻率、相位曲线和最佳主轴方位角分布图像,为我国东北N-MT资料的进一步处理和解释提供了基础性数据.此外,本文对比分析了朝阳测站中6条N-MT观测电极线上阻抗张量旋转值与观测值之间的差异,重点讨论了产生这种系统性偏差的各种因素,提出其主要因素可能来自“电场等效各向异性”效应,即测站附近的地壳内部存在与观测电极线尺度相比拟的横向非均匀构造,而测站各电极线沿不同方向跨越不均匀构造,此时各电极线上的电场分量不遵循同一电场矢量的分解准则,导致地表观测三角形内阻抗分量不满足统一的阻抗张量旋转规则.     

2000年日本鸟取县7.3级地震断层滑动分布的研究

利用Hori(2001)的反演方法,由同震GPS观测数据和水准观测数据,计算了2000年日本鸟取县西部地区7.3级地震断层滑动分布,并与余震分布进行了比较。结果显示,其最大滑动量达到5m,断层滑动分布主要集中在浅部,与余震分布并不一致。     

利用北斗GEO卫星观测数据分析中低纬地区小尺度VTEC时变率

利用北斗GEO卫星的静地特性,研究了中低纬地区小尺度垂直总电子含量(VTEC)时变率.给出了利用GEO卫星双频载波相位观测数据计算VTEC时变率(RVTEC)的公式,引入VTEC时变率指数(RVTECI)作为反映电离层活动程度的指标.通过对不同卫星观测结果的比较,确定30°为合适的RVTEC观测截止仰角.对处于南北纬30°之间5个测站共12天的观测数据进行了计算,结果表明:中纬地区VTEC变化最为平缓,RVTECI一般不超过0.01 TECu/s;赤道附近地区某些夜间时段电离层变化较为剧烈,RVTECI可达0.5 TECu/s,1 min内VTEC的变化量可达8 TECu.     

航线上重力基准构建的单点延拓方法

本文利用联合低阶重力位模型和地形数据的单点延拓方法将地面重力观测数据延拓到空中,构建了航线上重力基准,实现了两套航空重力观测数据的质量评价,实验结果表明:(1)单点延拓方法完全可以应用于航空重力观测数据的外部精度检核,效果良好,性价比高;(2)相比仅采用重力位模型或地形数据进行单点延拓计算,联合重力位模型和地形数据的延拓计算对航空重力观测数据的评估结果更加可靠;(3)联合低阶重力位模型和地形数据的单点延拓方法计算速度快、效率高.本文研究方法,也可适用于新研航空重力仪的标度因数、零偏等参数的在线标定.     

利用F10.7和Mg II构建太阳极紫外辐射长时间序列

利用SOHO/SEM在1996-2008年的太阳EUV观测数据,比较和评估F10.7和Mg II作为EUV代理参数的代表性,不能支持Viereck等关于Mg II是比F10.7更好的代理参数的结论.通过比较对两种参数的多种回归计算结果,确立双因子极大似然估计方法构建EUV计算模式,通过模型计算结果与SEM观测数据比较,表明该模型能够很好地重建EUV数据系列.利用该模式,构建了1978年11月以来的太阳极紫外辐射数据序列.     

利用Yabuki & Matsu'ura反演方法计算 2011年日本东北地区太平洋海域Mw9.0级地震同震滑动分布

Yabuki & Matsu'ura反演方法是利用ABIC最佳模型参数选取方法和平滑的滑动分布作为约束条件,由形变观测数据计算发震断层滑动分布.本文基于日本列岛同震GPS观测数据和发震断层曲面构造模型,利用Yabuki & Matsu'ura反演方法计算2011年日本东北地区太平洋海域M_w9.0级地震的发震断层同震滑动分布.反演结果表明,断层面上的最大滑动量为35 m,较大滑动分布在浅于30 km的震源中心上部,最大破裂集中在20 km深度的地方.其地震矩约为3.63×10²²N·m,对应的矩震级为M_w9.0.模拟结果显示Yabuki & Matsu'ura反演方法更适用于倾角低于40°的断层模型反演.最后,本文基于上述方法获得的发震断层滑动模型,利用地球体位错理论正演计算该地震在中国及其邻区产生的远场形变,正演计算结果基本可以解释由中国GPS陆态网络观测到的同震形变.     

水质量分布变化对重力观测的影响

利用全球7000多个陆地气象台站观测资料、TOPEX/Poseidon卫星测高资料、多层海水温盐度数据和负荷格林函数,采用数值积分法,计算了陆地水储量及海水质量分布变化引起我国20个测站的重力变化.结果表明,我国一些测站观测的重力受水质量迁移影响超过3μGal,可见在应用重力仪观测资料研究地球动力学时,有必要作水负荷改正.另外,文中还对武汉测站的超导重力仪观测的残差和水负荷的计算结果作了比较,并讨论了不同的地球模型及气象数据对计算测站重力的影响.     

冕宁地震台深井地电阻率观测系统简述及典型干扰数据特征分析

通过对冕宁台新旧观测场地电阻率数据曲线形态、完整率、连续率、月精度等指标进行对比分析,认为新场地三个测向的深井地电阻率观测精度、信噪比和相对变化量都优于旧场地浅地表观测数据质量.对新场地深井地电观测系统试运行期间的典型干扰数据进行分析,发现该台地电阻率的干扰源主要有自然环境、场地环境和观测系统故障.     

基于非全张量卫星重力梯度数据的张量不变量法

在非全张量卫星重力梯度观测数据的处理过程中,由于卫星姿态角误差、梯度观测数据误差和非全张量观测等原因,重力梯度值从卫星重力梯度仪系转换到地固系后,精度损失严重.本文研究了张量不变量法以解决上述问题.首先在重力梯度张量不变量线性化的基础上,建立了基于卫星轨道面的不变量观测模型,完整地推导了两类重力梯度张量不变量的球近似和...     

"多维"多尺度重磁位场数据融合方法

随着多种地球物理观测手段的应用和发展,不同观测手段采集的位场数据融合方法的研究越来越重要.根据观测维度、观测比例尺的不同,本文将重、磁力异常数据融合分维度(单维和多维)和尺度(单尺度和多尺度)进行讨论分析.针对多维多尺度重、磁位场数据观测位置、观测精度、数据基准的差异,采用空间域迭代法延拓、加权平均、回归分析等方法将不同观测手段获得的位场数据在同一标准下归算至同一平面或曲面,基于此提出了一套适用于多维多尺度重、磁位场数据融合方法;模型测试和实际资料处理结果显示,本文提出的融合方法在多维多尺度重、磁数据的融合中效果良好,计算误差小且符合位场数据的特点.本文所提出的多维多尺度重、磁位场数据融合方法适用于航空、地面和海洋观测的重、磁位场数据融合,具有良好的实用意义和推广价值.     

基于Oracle数据库的电磁扰动数据转换

为提供丰富的观测数据产品,编制基于Oracle数据库的电磁扰动产品数据计算软件.本软件能够将秒采样观测数据转换为分钟值、小时值、日均值数据,为应用人员提供丰富的研究数据.     

基于G856F的地磁野外测量处理程序

介绍了基于G856F的地磁野外测量处理程序的功能、使用方法及在山西地磁野外观测中的应用.认为该程序可在野外观测结束后立即计算观测结果,如果数据出现异常或误差超限,可现场查找原因或进行重测,从而提高地磁野外观测的效率和质量.