求解模型对欧拉矢量的影响

In the process of solving Euler vectors based on GNSS horizontal movement field, the number of estimated parameters can affect Euler vector results. This issue is analyzed through theoretical deduction and practical example in this paper. Firstly, the difference between the results of Euler vectors in different solving models is deduced. Meanwhile, based on GNSS horizontal movement field in the Chinese mainland from 2004 to 2007, two common models （ RRM and REHSM） are used to discuss the impact of solving models on Euler vectors and the follow-up study. The result shows that the maximum value of the difference in a block＇s entire rotation can reach 2. 6mm/a, and should not be ignored. Therefore, the results of horizontal movement are different using different kinematic block models, and this should be paid more attention in the analysis of crustal horizontal movement.     

不同数据用于欧拉方程的模型计算

在欧拉反褶积方法的研究中,过去常直接用原始位场异常值及其导数进行反演,由于数据的局限性,有时导致反演结果的发散.为改善反演结果的收敛性,本文中采用了垂直导数、水平导数及解析信号作为源数据进行欧拉反演.(1)推导了重、磁异常不同数据类型的欧拉方程;(2)找出了重、磁异常不同数据类型的欧拉方程之间的联系和差异;(3)通过模型计算比较,证明用导数或解析信号进行欧拉反演计算,可改善反演效果,提高反演结果的收敛性.     

基于ITRF2005的全球板块运动模型

利用最新的国际地球自转服务(IERS)发表的国际地球参考框架ITRF2005的速度场,建立了一个完全基于现代空间大地测量实测结果的现今全球板块运动模型ITRF2005VEL,并与地质模型NNR-NUVELlA及以往利用ITRF序列所建立的模型进行了比较.结果表明,从整体上看ITRF2005VEL与NNR-NUVELlA模型有较好的一致性,且较以往的ITRF序列的模型精度更高,但同时也存在着一定的差异性,这与观测台站、板块本身等相关因素都有关.本文还由ITRF2005VEL模型计算了构造板块的总角动量,结果不为零,说明ITRF2005框架有可能不满足无整体旋转的要求.     

基于欧拉矢量的中国大陆地壳水平运动自适应拟合推估模型

中国大陆地壳运动除具有整体趋势外,还具有十分复杂的局部变化特征,因此建立中国大陆地壳运动模型应当同时顾及板块的长期运动趋势和小范围的局部运动变化.本文将欧拉矢量模型作为地壳运动的趋势项模型,将局部变化信息作为信号项模型,采用趋势变化与随机变化相结合的拟合推估法进行解算.拟合推估模型的信号向量易受先验协方差函数的误差影响,本文采用方差分量估计构造自适应因子调整随机信号的先验协方差阵,从而自适应地调整信号向量与观测向量的先验权比,达到平衡随机观测误差和随机信号对形变参数估值影响的目的.计算结果表明,基于欧拉矢量的拟合推估模型更好地顾及了位移场的整体变化趋势和局部运动特征,提高了中国大陆地壳水平运动场的计算精度.     

复数矢量球函数的理论研究

给出了正交归一化的复数矢量球函数的定义和性质。复数矢量球面调和函数是球体位移的研究基础,采用复数矢量球面调和函数可把位移场的球形场和环形场一同表示出来,为各种动力学效应的统一研究提供了理论依据。同时还研究了在球坐标系下弹性就变张量与复数矢量球函数的一些关系,为地球弹性动力学方程的解算提供参考。     

古地磁欧拉极及对华北地块的应用

本文详细讨论了古地磁欧拉最佳小圆拟合法的基本原理和实际方法，并提出了提高拟合速度和优度的措施，同时利用此方法对前人取得的华北地块中石炭世至晚侏罗世这一时间的古地磁视极移曲线进行了最佳小圆拟合。     

球面几何构建欧拉极的方法

欧拉旋转可用来描绘刚性板块在球体上的运动,并用于古大陆的再造.一般构建欧拉极的方法是通过匹配等时线、海底磁异常、转换断层和古地磁数据等来确定的.但是在一些地质环境下,不能直接构建欧拉极.本文直接根据古板块位置图,从几何上构建欧拉极,从而为古大陆重建提供了一个新途径.本文主要介绍了在球面几何上构建欧拉极的方法原理,并编写了相应的程序功能.选取了印度大陆0 Ma和50 Ma的位置图作为实例计算了欧拉极.由本方法计算得到的欧拉极将板块旋转后的位置,与原始古板块位置图几乎完全重合.说明本文构建欧拉极的方法是正确而有效的.     

结合欧拉矢量的反演算法构建青藏高原东北缘地壳运动速度场模型

探讨地壳运动速度场模型的构建方法,提出结合欧拉矢量的维多样性动态权重粒子群算法构建地壳运动速度场模型。通过模拟算例验证该算法的稳定性和有效性,建立的速度场模型与线性权重粒子群算法和非线性权重粒子群算法的计算结果相比具有较高的精度,且收敛速度较快。利用青藏高原东北缘1999—2013年中国地壳运动观测网络观测到的GPS水平速率结果,在块体划分和模型辨识的基础上,建立青藏高原东北缘地壳运动速度场模型,并将其与最小二乘配置法的计算结果进行比较,结果表明改进的粒子群算法建立的地壳运动速度场模型具有较高的精度。     

华北地区近期地壳水平运动与应力应变场特征

利用华北GPS监测网 1 992年、1 995年、1 996年的观测资料 ,应用最小二乘配置给出了华北地区相对水平位移场、应变场的分布图像 .经初步研究表明 :华北地区 1 992—1 995年间的水平位移和应变场表现为整体性不均匀的压性运动 ,1 995— 1 996年测区东部仍以水平压性运动为主 ,但测区西部则主要表现为张性运动 .水平运动 （方向、大小 ）发生显著变化和应变高值区的地带主要位于块体边界带和主要断裂带附近 .燕山断块南边界的北东向断裂存在着较显著的左旋运动 .区内最大剪应变、面膨胀的高值区在天津、北京、唐山一带 .结合非连续变形数值分析方法 （DDA）初步分析认为 ,1 992— 1 995年GPS观测结果显示的华北地区存在东、西部构造应力作用的明显差别 ,华北东部以东西向压应力作用为主 ,而西部的南北向构造应力作用又明显大于东部 .     

华北地区近期地壳水平运动与应力应变场特征

利用华北GPS监测网 1 992年、1 995年、1 996年的观测资料 ,应用最小二乘配置给出了华北地区相对水平位移场、应变场的分布图像 .经初步研究表明 :华北地区 1 992-1 995年间的水平位移和应变场表现为整体性不均匀的压性运动 ,1 995- 1 996年测区东部仍以水平压性运动为主 ,但测区西部则主要表现为张性运动 .水平运动 （方向、大小 ）发生显著变化和应变高值区的地带主要位于块体边界带和主要断裂带附近 .燕山断块南边界的北东向断裂存在着较显著的左旋运动 .区内最大剪应变、面膨胀的高值区在天津、北京、唐山一带 .结合非连续变形数值分析方法 （DDA）初步分析认为 ,1 992- 1 995年GPS观测结果显示的华北地区存在东、西部构造应力作用的明显差别 ,华北东部以东西向压应力作用为主 ,而西部的南北向构造应力作用又明显大于东部 .     

顾及抗差性的GNSS多尺度地壳运动场构建方法

精准提取应变率场是利用GNSS(Global Navigation Satellite System)观测资料开展区域地震危险性分析工作的重要环节, 而其中的关键点是构建精确的地壳运动场模型.本文提出一种顾及抗差性的多尺度地壳运动场构建方法, 在多尺度球面小波方法的基础上, 考虑到测站观测数据质量不优或观测环境不好等因素产生的局部非构造变形对提取地壳形变场的影响.通过在实际观测资料中加入粗差的模拟分析验证本方法具有好的抗差性.通过对观测数据高精度统一解算, 获取川滇地区最新的震间GNSS速度场, 抗差迭代算法从中检测出20%的降权站点和2%的除权站点, 而这些测站的观测标墩均建设在土层上, 且数据解算的误差较大, 表明存在明显的局部非构造变形, 如果在解算应变率场时与其他测站作等权处理或基于解算的形式误差确定权重, 必然会在解算的应变率场产生显著粗差, 进而影响到地震趋势判定或危险性分析等重要工作.

     

An analysis on short—wave components of the global stress field

Introduction The intraplate stresses are affected by two types of forces: one is called the first-order stresses (Zoback, 1992) or the long-wave components of stresses, which are closely related to the driving mechanism of the plate motions such as ridge push, continental collision resistance, slab pull (suc-tion) and mantle drag force (Richardson, 1992; WEI, 2000). This type of stress field is character-ized by a uniform feature in very large horizontal scale with 50 times or more of the th…     

Construction of GPS MonitoringNetwork,Characteristics of RecentCrustal Movement,and Relationbetween Stress Field andEarthquake Activity in the NingxiaRegion

INT正ODU**!0厂Nln以la Is located In the nol’th segmentdthe North－southselsnuc Belt,Whlchls thejuncturedthe Qinghai－Xizang(Dbet)block,Alxa block and Odos block．In that region,the tectonic activity。svery strong;。y strong earthquakes occu。d In history,such as the 1739 ingluo M。8．0 earthqu拙eandl920 Halyuan M。8．5 eafthquake,etc．; Merately strongeaFthquakes are still veryactlve atpresent,onlyln the Wuzhong－un斟m area,there have been7 e血hquakes tvlth M。＞5．0 to occ…     

基于欧拉反褶积方法估算华北地区重力变化的等效源参数

利用华北地区2009—2014年绝对重力与相对重力多期重复观测资料, 得到不同时空尺度的华北区域重力场动态变化图像。 采用欧拉反褶积方法, 通过对理论模型试算, 获得最优反演参数, 对引起华北地区重力场变化的场源深度和空间分布规律进行了反演计算和分析。 结果表明: 当构造指数为1时, 适合对流动重力场变化数据进行反演。 实际资料的反演结果在2009—2014年间的场源位置集中于河套断裂带。 本文研究方法可以用于重力场反演和以场求源的定量研究。 本工作为流动重力变化信号的反演和解释提供一条新思路, 也为构建地震重力预报指标体系提供了定量依据。     

云南地区GNSS应变场时空演化与地震事件关联分析

采用GAMIT/GLOBK软件对云南境内及邻区近400个GNSS测点1999～2018年的观测数据进行解算,在各个测点时间序列和速度场的基础上,采用克里金插值方法分时段估计该区域在1999～2004年、2004～2007年、2009～2013年、2013～2015年、2015～2018年共计5个时间区域应变率场;根据区域地壳面应变率和最大剪应变率的空间变化以及相应时段之后3年内的MS≥5. 0地震事件分布特征,分析发现:绝大部分震例发生在面应变高梯度带的张压转换区和最大剪应变高值区,可见研究区各个观测时段GNSS应变率场对后期1～3年内的中强震发生区域有一定的指示意义;以2014年盈江6. 1级、鲁甸6. 5级和景谷6. 6级地震为样本,建立监视块体获取应变时间序列,分析发现:地震前三个月左右均出现震中附近短期应变趋势改变、快速增强、转折的现象,这些形变异常变化或许反映了发震区应力-应变积累在接近临界破裂状态时的非线性调整,为地震短临预测尤其是时间要素的判断提供参考。     

地壳受力的向量场及奇点的存在性

在地壳的每一点,由于受力定义了整个地壳上的向量场,然后利用同胚,将其转化为球面上的自映射,并根据球面的欧拉示性数等于2,说明该自映射至少有一个不动点,即向量场至少有一个奇点.这样,从理论上解释了地震、海啸、火山等自然现象发生的必然性.     

Estimating Field Source Parameters of Gravity Change in North China Based on the Euler Deconvolution Method

Based on the absolute and relative gravity observations in North China from 2009 to 2014, spatial dynamic variations of the regional gravity field are obtained. We employed the Euler deconvolution method and the theoretical model to get the best estimates of parameters. Gravity field change caused by the depth and distribution in North China is calculated by back analysis. The results show the structural index that equals 1 is suitable for inversion of the gravity variation data. The inversion results indicate that the depths of anomaly field sources are spread over the Hetao fault. The research method of this paper can be used in the quantitative study on the field source and may shed new light on the interpretations of gravity change, and also provide quantitative basis for earthquake prediction index criterions based on the gravity change.