2008年10月西藏当雄MW6.3地震震后形变提取与余滑反演

2008年10月6日西藏当雄发生M_W6.3地震.本文利用震后2008年10月26日至2010年8月22日的16期ENVISAT ASAR数据,通过小基线集干涉测量、误差校正与MInTS(Multiscale InSAR Time Series)技术提取高精度的震后形变场,利用SDM(Steepest Descent Method)方法反演断层震后余滑演化过程,并分析震后余滑与同震滑动的关系.结果表明：当雄M_W6.3地震的近场震后形变场主要位于断层西侧,在时间演化上具有明显的对数函数衰减规律;震后余滑主要集中于断层中南段深0~15 km区间,最大的余滑量约0.07 m,位于断层深约9.28 km处,滑动角约-103°;震后余滑引起的地震矩能量M₀与矩震级M_W在时间演化上具有指数函数递增规律;当剪切模量μ=32 GPa,震后665天余滑释放的地震矩能量约为1.92×10¹⁷N·m,约占同震滑动释放地震矩的4.8%,相当于矩震级M_W5.46;虽然震后余滑已经延伸到断层浅部0~5 km区间,但由于余滑量相对较小,没有改变同震滑动在断层浅部区域的滑动亏损现象,这可能是2010年11月30日该区域又发生M_W5.3级余震的主要原因之一.

     

基于DFT模型的大场景InSAR图像配准

图像配准是实现干涉合成孔径雷达(InSAR)高精度相位提取及地形高程反演的关键,大场景图像的高效高精度配准成为近年高分宽幅InSAR成像应用研究的难点问题之一。由于大场景图像中不同区域偏移量及变化规律差异较大,传统最大相干系数配准方法需多分块及插值处理,面临计算量大且配准精度低等问题。针对此问题,本文提出一种基于DFT模型的大场景InSAR高效高精度图像配准算法。该方法利用最小均方差准则构建InSAR复图像配准的DFT模型,采用四叉树自适应分块及矩阵相乘DFT快速重采样配准方法,实现大场景InSAR图像各子块区域的高效高精度亚像素配准。仿真和实测数据验证本文算法的有效性,结果表明该算法不仅可实现大场景InSAR复图像亚像素级配准,还具有较高的运算效率,其运算效率相对于传统FFT配准方法通常可提升3倍以上。     

2017年四川九寨沟Ms7.0地震InSAR同震形变场及发震构造探讨

2017年8月8日四川省九寨沟县发生M_s7.0级地震,构造部位处于青藏高原东缘的巴颜喀拉地块东北角,震中位置是岷江断裂、塔藏断裂、虎牙断裂和雪山梁子断裂围闭的空震区。哪条断裂发震,如何界定其与周边活动断裂的关系,与青藏高原东缘近年来发生的大地震是否有成因联系等问题对于理解该区域现今构造活动模式、预判地震发展趋势和部署地震地质灾害防控等工作具有重要意义。利用地震前后两期Sentinel-1合成孔径雷达数据对地表同震形变场进行了InSAR测量,获取了极震区约2000 km²范围内的雷达视线向变形（-13~28 cm）和运动方向,呈现为主动盘单侧走滑兼逆冲的变形模式,结合震源机制、断裂展布、构造背景和近年地震迁移的分析,揭示了控震构造是巴颜喀拉地块北缘边界断裂弧形旋转体系的尾端构造,发震断层是该断裂系中塔藏断裂的南段,并有与虎牙断裂贯通的趋势,因此,应重视本次地震与虎牙断裂之间的空震区未来的强震危险性问题;从区域上看,此次九寨沟地震可能与汶川地震具有一定的时空成因联系,因在巴颜喀拉地块南北边界断裂破裂基本贯通的条件下,2008年汶川地震诱发的东缘中部锁固破裂导致块体加速向东挤出,2013年鲁甸地震又释放了东缘南段挤压构造应力,从而进一步加剧了东北角的应力集中,促使九寨沟地震的发生。     

基于SAR影像的喜马拉雅山珠穆朗玛峰地区冰川运动速度特征及其影响因素分析

合成孔径雷达(synthetic aperture radar, SAR)具有其全天候、全天时、穿云透雾的工作能力, 广泛应用于山地冰川动态监测中. 利用2006年6-9月三期ALOS/PALSAR雷达影像, 采用偏移量跟踪技术, 提取了喜马拉雅山珠穆朗玛峰(珠峰)区域的冰川运动速度, 分析了区域内冰川运动速度空间差异及其影响因素. 结果表明: 研究区31条山谷冰川平均运动速度为9.3 cm·d^-1, 总体上以珠峰-洛子峰南北向山脊线为界限, 东侧和东南侧冰川日均运动速度(11.1 cm·d^-1)普遍高于北部和西北部冰川日均运动速度(5.4 cm·d^-1). 冰川消融区非表碛区冰川平均运动速度为表碛覆盖区平均运动速度的2.2倍, 冰面湖的发育在一定程度上加剧冰川运动速度波动. 在气候与非气候因子共同作用和相互间的此消彼长中, 研究区65%的冰川的运动速度自中值高度往下显著减小, 16%的冰川自中值高度往下呈显著增大趋势, 19%冰川消融区运动速度无显著变化趋势.     

九寨沟M_s7.0地震的InSAR观测及发震构造分析

2017年8月8日四川省九寨沟县发生M_s7.0地震.本文基于Sentinel-1 SAR影像,利用InSAR技术获取了此次地震的同震形变场,反演获得同震滑动分布,计算了同震位错对余震分布和周边断层的静态库仑应力变化,并对发震构造进行了分析讨论.结果表明:①InSAR同震形变场显示,九寨沟地震造成地表形变最大量级约为20 cm(雷达视线方向),同震形变存在非对称性分布特征.②同震位错以左旋走滑为主,主要发生在4~16 km深度,最大滑动量约为77 cm,位于9 km深处.反演得到的矩震级为Mw6.46.同震错动未破裂到地表.③大部分余震发生在库仑应力增加区.此次地震增加了震中周边地区一些断裂的库仑应力,如东昆仑断裂带东段、龙日坝断裂、虎牙断裂等.④东昆仑断裂东段的未来地震危险性值得关注.⑤九寨沟地震的发震断层为树正断裂,可能是虎牙断裂的北西延伸隐伏部分,此次地震是巴颜喀拉块体南东向运动受到华南块体的强烈阻挡过程中发生的一次典型构造事件.     

现代地壳运动与地震监测预报研究的现状和发展趋势

综述了近年来美国、日本的GPS连续观测系统和INSAR空间测地技术在现代地壳运动和地震监测预报研究中的应用现状和发展趋势。     

北京某山坡湿陷性土层地基处理可行性分析与建议

北京市密云区某拟建别墅工程,属于山体边坡地基基础工程,任务要求对湿陷性地基土层进行湿陷性消除,并将其地基土承载力特征值提高到fak≥130 kPa。此外,还需要考虑边坡稳定性及其永久性加固问题。对该项目进行了深入分析、研究,并采用概念设计方法,对其进行设计预演,发现该项目存在工序复杂、不易实施、在永久性使用方面存在一定隐患、并具有责任不当等问题。对该地基基础工程进行了重新分析,建议取消地基处理措施,改用嵌岩桩作为拟建物的基础。     

轨道误差对InSAR数据处理的影响

轨道参数是InSAR数据处理中一个重要的参数,对从最初的图像配准到最后的高程值或形变值图像生成都有着重要的影响.含有误差的轨道参数造成基线误差以残差条纹的形式存在于干涉图中.完全去除残差条纹要求轨道绝对精度低于1 mm,目前的定轨精度远不能满足要求.这里推导了轨道误差和残差干涉条纹的关系,分析了轨道误差对高程值和形变值影响的特点,提出采用二次多项式拟合的方法去除残差干涉条纹,并以巴姆地区的ENVISAT数据证明了提出方法的有效性.     

Heterogeneous surface displacement pattern at the Hatchobaru geothermal field inferred from SAR interferometry time-series

We estimated surface displacements using persistent scatterer SAR interferometry (PS-InSAR) around the Hatchobaru geothermal field, Japan, from 18 ALOS/PALSAR images acquired from July 2007 to December 2010. Generally, geothermal fields, covered with natural targets such as rocky terrain and vegetation, have been one of the difficult targets for PS-InSAR analysis. However, we applied space adaptive filtering to increase the number of pixels for measuring surface displacement. The results of our analysis demonstrate ground subsidence with decaying velocity over the observation period around the geothermal field. The spatial pattern of ground subsidence includes sharp boundaries of subsidence that can be interpreted as fault traces. We demonstrated the usefulness of PS-InSAR analysis with the space adaptive filtering to estimate surface displacements with high spatial resolution and high spatial density around a geothermal field.     

利用DS-InSAR技术监测黄河三角洲地表形变

针对黄河三角洲地区湿地及农田多、范围大,导致PS-InSAR技术难以获取高密度地表形变信息的问题,提出一种基于分布式目标InSAR（DS-InSAR）的黄河三角洲地表形变监测方法。该方法通过置信区间估计选取同质像元点,利用特征值分解方法计算主散射体对应相位值以达到相位优化的目的,再根据时空相干性确定分布式目标,最后建模解算时序地表形变信息。以26景Sentinel-1A影像为数据源,提取2019-12～2020-12期间黄河三角洲地区的地表沉降信息,与PS-InSAR方法结果相比,点位密度提高5.56倍;两种方法获取的同名点对形变速率的相关系数为0.727,说明两者具有很好的一致性。实验结果表明,研究区内存在4处明显沉降区域,最大沉降速率达-238 mm/a,经分析及实地调查验证,其主要影响因素为地下卤水及油气开采。