基于D-InSAR的樱岛火山喷发形变场解算与反演分析

采用D-InSAR技术对覆盖2020年日本鹿儿岛市樱岛火山喷发前后时段的2景降轨Sentinel-1A数据进行处理,获取日本鹿儿岛市樱岛火山南岳火山口喷发事件引起的地表形变场,并在此基础上结合点源Mogi模型对樱岛火山喷发时的岩浆源进行反演分析。结果表明,2020-07-28~08-09期间沉降主要集中在火山中心地区,最大沉降区域为火山中心地带,最大沉降量与平均沉降量分别为5.5 cm与2.85 cm;抬升主要发生在火山边缘区域,最大抬升量与平均抬升量分别为5 cm与2.24 cm,抬升可能与艾拉火山口下方的岩浆活动有关。根据点源Mogi模型反演得到岩浆源深度和体积变化量分别为1.016 km和-0.139×10⁶ m³,岩浆源位于南岳下方,与南岳火山口的喷发活动有关。     

函数模型与随机模型双约束的GPS-InSAR数据融合方法建立三维形变场

针对InSAR技术研究地表三维形变时监测信息不足的问题，以GPS监测信息为先验信息，建立附有随机模型约束的地表三维形变模型。考虑到SAR卫星极轨方式运行导致LOS向观测量对南北向形变不敏感的问题，以GPS南北向形变观测值作为强约束，构建三维形变解算的函数约束条件。模拟数据与西安地区实测数据的计算结果表明，基于随机模型与函数模型共同约束的地表三维形变参数最小二乘解的精度优于仅有函数模型约束或仅有随机模型约束及无任何约束的参数解精度。     

利用SBAS-DInSAR技术提取腾冲火山区形变时间序列

基于6景JERS-1 L波段SAR影像,利用小基线集-合成孔径雷达干涉测量技术,通过线性形变相位、非线性形变相位、大气延迟相位以及地形残差相位的分离,提取了腾冲火山地区1995—1997年间地表形变时间序列（雷达视线向）,与2003—2004年的GPS观测结果对比表明,SBAS-DInSAR技术提取地壳形变的精度可达亚厘米级。时间序列形变显示胆札-高田断裂两侧形变差异性显著,可能与其下方存在的地壳岩浆囊的活动有关。打鹰山地区地表形变揭示其下方可能存在隐伏断裂。     

基于GNSS观测研究2015年尼泊尔MW7.8地震震后地壳形变特征及其机制

收集及处理尼泊尔境内的GPS连续观测站和中国藏南地区的GPS基准站数据,获得2015年尼泊尔MW7.8地震震后3 a的GPS水平形变场。结果显示,尼泊尔地震的震后形变主要分布于尼泊尔北部及中尼边境区域,且东西方向形变较小,南北方向形变较大,整体继续向南运动,最大震后位移约为10.93 cm。采用孔隙弹性回弹模型计算的理论地表位移远小于GPS观测值,无法解释GPS观测到的震后形变。采用震后余滑模型反演的结果表明,震后余滑主要集中在断层的下倾延伸部分,且空间分布较广,余滑释放的地震矩为1.09×1020 Nm。采用PSGRN/PSCMP程序计算粘弹性引起的理论地表形变结果显示,粘弹性松弛模型不能解释近场GPS观测值,但在远场区域的运动方向与GPS观测值一致。采用粘弹性松弛和震后余滑组合机制模型进行反演,余滑释放的地震矩降为1.08×1020 Nm,且空间分布更加集中。研究结果表明,组合机制模型在保证了模型拟合精度的基础上,反演结果与应力驱动模型反演结果更接近。     

采用GPS观测值进行三维形变分析的构想

在已知地表度量张量的基础上,给出了根据测站的GPS复测值计算地表形变张量、主方向、线应变、面应变、高斯曲率与平均曲率的有关公式,并着重对地表的内蕴与外蕴形变概念作了阐述,最后指出利用地表内蕴与外蕴形变量能较好地描述地表的三维形变。     

2021年云南漾濞地震InSAR监测及断层滑动分布反演

采用双轨D-InSAR技术获取2021-05-21云南漾濞6.1级地震的升降轨同震形变场,并以此为约束,基于Okada均匀弹性半空间位错模型反演发震断层的几何参数与同震滑动分布,最后通过正演模型预测地震的同震三维形变场,分析地震的形变特征及发震断层的几何形态和构造运动。结果显示,在InSAR形变场中,升轨LOS向最大形变量约10 cm,降轨LOS向最大形变量约13 cm;利用Okada模型反演得出,发震断层长约12.53 km,宽约6.44 km,位于地下3～9 km处,断层走向约136.34°,倾角约83.79°,滑动角约-175.46°;断层的同震滑动主要分布在沿走向9～21 km与沿倾向向下3～9 km范围内,最大滑动位移出现在地下4 km深度处,滑移量达到0.64 m;通过正演模型的地表三维形变预测发现,垂直方向最大抬升约3.22 cm,最大沉降约4.04 cm。总体而言,此次地震的发震断层类型为右旋走滑,且带有正断分量,初步判断地震发生在一条未知的NW-SE向断裂带上,该断裂可能属于维西-乔后-巍山断裂和红河断裂北段的次级或分支断裂。     

基于GPS观测值研究日本M_W9.0地震震后形变机制

基于GPS观测的震后水平位移对2011年日本MW9.0地震的震后形变特征进行研究。震后近5a(截至2015-12),震后水平位移累积达到东向60~165cm,南向20~65cm,距离震中较远处的G104、G105及J192站点观测到的震后位移累积变化已超过同震位移,且震后形变还在持续。联合震后余滑和粘弹性位错理论模拟震后形变,利用这2种作用机制对震后GPS水平形变进行解释。研究表明,震后余滑在震后形变最初阶段起主要作用,但随着时间增长而逐渐衰减,粘滞性松弛作用的贡献随着时间增长而变大,GPS观测到的震后地表形变可由这2种机制结合得到较合理的解释。利用震后GPS形变模拟估算地震区域的地幔粘滞系数在1.5×1019 Pa·s量级。     

新疆呼图壁地下储气库的InSAR形变监测与模拟

利用TerraSAR-X卫星2013-08~2014-08的17景雷达影像,采用小基线集（small baseline subset,SBAS）InSAR技术获取呼图壁(HTB)地下储气库（underground gas storage,UGS）运行期间的地表形变序列,并结合UGS注（采）气井口的压力数据,采用多点源Mogi模型,对HTB UGS的形变场进行模拟。结果表明,整个UGS区域的形变特征为非连续分布,形变与注（采）气压力变化具有较好的相关性;注（采）气期间沿卫星视线向(LOS)的形变峰值分别为10 mm和-8 mm;采用自适应前向搜索法,基于多点源Mogi模型初步模拟注（采）气期间的形变过程,当UGS的注（采）气平均气压为18 MPa和15 MPa时,LOS的形变可达7 mm和-4 mm,地表形变的大小与注（采）气井口密度有关;UGS的储气分布呈非均匀状态,即地下气库结构复杂多变。     

合成孔径雷达信号子带干涉技术及其应用

分别选用两景Radarsat-2和TerraSAR-X数据,应用子带干涉技术获取夏威夷火山国家公园及其附近区域的沉降量,并将其与偏移跟踪技术结果进行对比,证明了子带干涉技术的地表形变监测能力,且分辨率越高、带宽越宽的TerraSAR-X数据得到的监测精度越高。用两景PALSAR数据,应用子带干涉技术获得唐洞煤矿地表方位向位移,与同期GPS监测数据进行对比,进一步验证了子带干涉技术的形变监测能力。     

日本Mw9.0强震的高频GPS形变波特征初探

利用高频GPS数据长基线双差精密定位技术,获取2011年3月日本Mw9.0强震近场和远场GPS连续参考站的瞬时动态形变。分析表明：距震源较近的IGS站（USUD）的同震动态形变水平方向幅度超过60 cm,震后水平形变达到30 cm;位于远场的GPS连续参考站（JLYJ）的动态形变,清晰地给出了地震面波引起的地表形变波形。对JLYJ站的谱分析表明,地震面波引起的动态形变周期为10～17 s,并出现了两次峰值,东西向和垂向上的动态位移量远大于北南向,说明JLYJ站的形变主要由Rayleigh波及其谐波所主导,Rayleigh波的传播速度为3.5 km/s,波长约50 km。     

����GPS��InSAR����ͬ��ر���άλ�Ƴ�

???GPS???λ????InSAR??LOS?α???????????GPS??λ?????LOS????????GPS??????LOS??λ????????Delaunay???????????????????????壬??InSAR??LOS?α??????о?????????????????????λ??????????GPS???????????????????????????????λ?????????Χ???????????????????????????????λ???????????????????????InSAR??LOSλ?????????????????????????????????????????????α????????????????????????á?     

����GPS��InSAR�۲��������봨������ά�ر��α�

?????????????????????InSAR???????????????????α??????????????????λ???????λ??????????????????????????????GPS??InSAR??????????????λ????????????????????????м????????????????????????????????????????????     

BDS/GPS单历元阻尼LAMBDA算法及其在边坡变形监测中的应用效果分析

论述了BDS/GPS双系统组合下的单频单历元阻尼LAMBDA算法原理,并基于该算法实现一种实时远程变形监测系统。以理县某边坡变形监测应用为例,比较了GPS、BDS及其组合系统下的观测条件和阻尼LAMBDA算法定位效果：BDS的观测卫星数多于GPS,PDOP值相较于GPS也更加稳定;BDS单历元解标准差在N、E、U方向上分别为0.40 cm、0.31 cm、1.00 cm,优于GPS的0.61 cm、0.40 cm、1.74 cm。长期实验结果表明,在我国南部边坡监测中,BDS相比GPS具有一定优势。理县边坡23个月内位移量显著,3个监测站平均位移在N、E、U方向上分别达到8.70 cm、43.63 cm、18.03 cm。强降雨是引起土质山体滑坡的主要因素,使用累积降雨量和累积位移量建立线型回归模型,线性相关系数在0.98以上。在边坡监测中,可将实时的位移数据和降雨数据作为滑坡预警的重要依据。     

Absolute calibration of HY-2A and Jason-2 altimeters for sea surface height using GPS buoy in Qinglan,China

GPS buoy methodology is one of the main calibration methodologies for altimeter sea surface height calibration.This study introduces the results of the Qinglan calibration campaign for the HY-2 A and Jason-2 altimeters.It took place in two time slices;one was from August to September 2014,and the other was in July 2015.One GPS buoy and two GPS reference stations were used in this campaign.The GPS data were processed using the real-time kinematic(RTK)technique.The final error budget estimate when measuring the sea surface height(SSH)with a GPS buoy was better than 3.5 cm.Using the GPS buoy,the altimeter bias estimate was about-2.3 cm for the Jason-2 Geophysical Data Record(GDR)Version 'D' and from-53.5 cm to-75.6 cm for the HY-2 A Interim Geophysical Data Record(IGDR).The bias estimates for Jason-2 GDR-D are similar to the estimates from dedicated calibration sites such as the Harvest Platform,the Crete Site and the Bass Strait site.The bias estimates for HY-2 A IGDR agree well with the results from the Crete calibration site.The results for the HY-2 A altimeter bias estimated by the GPS buoy were verified by cross-calibration,and they agreed well with the results from the global analysis method.     

应用GPS监测城市地表形变的初步分析

根据实例,总体介绍了用ＧＰＳ监测城市地表形变的观测及数据处理方法。首先根据ＧＰＳ获得的基线结果进行３维自由网平差,由定期观测获得了各个监测点的水平位移和垂直位移。然后针对城市监测范围较小,在将其按照矩形格网划分后,给出了两种拟合方法：线性内插法和多面函数拟合法。据此,计算获得了城市地表的水平位移场和垂直位移（沉降）面,并对两种方法获得的结果作了比较。根据计算结果绘制了水平位移场、垂直沉降面及其等值线图,可以直观地了解地表形变信息。同时对于用ＧＰＳ监测城市形变的观测和数据处理中需要注意的问题作了简要的分析。