摘要:利用2016年1月—2021年12月的89景哨兵一号雷达影像,采用相干点目标分析合成孔径雷达干涉测量技术(IPTA-InSAR)进行数据处理,获取了武汉地区的地面沉降信息,并联合GNSS基准站观测成果对InSAR监测的形变进行分析。结果表明2016—2021年期间,武汉市主城区形成了汉口、沙湖北和白沙洲3个较为明显的沉降中心,地面沉降呈连片化发展趋势。汉口地区的地面沉降速率最高,部分区域沉降超过10 mm/a,典型沉降区的地面沉降过程伴随有一定的波动特征,其中下沉趋势在不同年间会有区别,2016—2020年间特征点缓慢下沉,2018年开始加速下沉,至2020年下沉速度再次放缓。GNSS与InSAR特征点沉降分析表明2种技术的监测结果整体上具有很好的一致性。
地面沉降又称为地面下沉或地陷,是在人类工程经济活动影响下,由于地下松散地层固结压缩,导致地壳表面标高降低的一种局部的下降运动(或工程地质现象)。城市地面沉降是我国常见的城市地质灾害之一,它可能会导致建筑物、桥梁、道路等基础设施的破坏,危及人民生命财产安全[1]。常规地面沉降监测与分析方法(如水准测量、GNSS、加速度传感器等)受空间分辨率低、目标监测点布设困难、作业成本高等因素限制,难以实现大空间尺度与高空间分辨率的变形监测与分析[2,3,4]。合成孔径雷达干涉测量(interferometric synthetic aperture radar,In SAR)技术能够提取地表微小形变信息,且具有受天气影响小、覆盖范围大、成本低、测量精度高等特点,弥补了传统形变监测的不足,已经在众多城市地面沉降监测中进行了应用[5,6]。
本文选取覆盖武汉地区2016年1月—2021年12月的89景Sentinel-1A影像,采用相干点目标分析合成孔径雷达干涉测量技术(IPTA-In SAR)获取武汉地区的地面形变速率,对典型沉降区域进行时序分析,并结合GNSS观测成果对In SAR监测的精度进行了验证分析。
1、研究区概况与实验数据
1.1 研究区域
武汉市地处江汉平原东部、长江中游,是国际湿地城市,市内江河纵横、百湖密布,水域面积占全市总面积的1/4。武汉市地貌属鄂东南丘陵经汉江平原东缘向大别山南麓低山丘陵过渡地区,中间低平,南北丘陵、岗垄环抱,北部低山林立。
1.2 GNSS数据
本文同时收集了研究区域内的7个CORS基准站2016—2022年的观测数据,利用GAMIT10.71软件对观测期间的数据进行了基线解算,解算时将中国境内及周边的部分IGS基准站纳入组网解算,基线处理策略如表1所示。获取基线解文件后,利用GLOBK软件进行平差处理,进而获取各基准站统一参考框架下的坐标。
表1 基线处理策略
1.3 In SAR数据处理方法
本文收集覆盖武汉地区的2016年1月—2022年9月的101景哨兵一号SAR数据,成像模式为宽幅干涉模式(IW)的单视复数(SLC)影像,数据的基本信息如表2所示。本文采用相干点目标分析合成孔径雷达干涉测量技术(IPTA-In SAR)来进行数据处理。IPTA技术是GAMMA软件提出的干涉点目标分析技术,它将覆盖研究区的多景SAR影像进行统计回归分析得到相位和幅度信息,能够有效识别建筑,裸露的岩石、人工角反射器等大量的点目标[7,8]。本文经过公共主影像的优化选取、SAR影像配准、差分干涉处理、相干点目标提取、IPTA相位解缠、大气相位等残余相位分离、形变速度场解算等步骤,获取了基于InSAR监测的武汉市地表沉降信息。
表2 SAR数据基本信息
2、结果分析
2.1 研究区域总体变化特征
本文通过IPTA技术对哨兵一号数据进行处理,进而获取了影像覆盖地区的形变监测速度。整体而言,2016—2021年研究区域的地面沉降呈平稳态势,年均形变速率在-5~5 mm/a的地面形变点数量占总地面形变点数量的98.28%(负值代表地面沉降、正值代表地面抬升,下同),年均形变速率在-5 mm/a以下的地面形变点数量占总地面形变点数量的1.71%,表明监测区存在部分区域处于地面沉降状态。
进一步选择武汉主城区为研究对象,对城市地面沉降进行分析。武汉主城区分割后得到3 442 256个点位,其中形变速率在-10 mm/a到-5 mm/a的地面形变点数量占总地面形变点数量的1.4%;形变速率在-5 mm/a到0 mm/a的地面形变点数量占总地面形变点数量的33.80%。年均形变速率在0 mm/a以上的地面形变点数量占总地面形变点数量的64.39%;年均形变速率在+5 mm/a以上的地面形变点数量,还存在0.43%的点位形变速率在-10 mm/a以上。表明监测区内存在一定数量的地面沉降趋势,其中还有少部分区域沉降趋势极为明显。
在监测周期内,部分地区沉降量呈逐渐增加的趋势,沉降量较大的区域范围呈逐渐扩大并有连成一片的趋势。截至2021年12月,武汉市主城区形成3个较为明显的沉降中心,即汉口地区、沙湖北部地区和白沙洲地区,而地面抬升的区域主要分布在武汉市主城区的东南部和西南部。从沉降空间分布上看,地面沉降速率较高(<-5 mm/a)的区域主要是长江左岸的汉口地区和长江右岸的沙湖北部地区、白沙洲大道附近地区。其中,汉口地区的地面沉降速率最高,部分区域的地面沉降速率超过10 mm/a,沉降明显区域连片分布,沉降面积最大。其次是沙湖北地区,沉降明显区域由沙湖北部连绵至长江天兴洲段南岸;除部分沉降突出的区域外,白沙洲地区的地面沉降速率较汉口地区和沙湖北部地区偏低,其分布沉降明显区域集聚程度低于汉口地区和沙湖北部地区。
2.2 典型沉降区特征点分析
为了分析地表沉降随时间的演变,选取了分布在不同沉降区的沉降特征点对重点地区(后湖公园及周边地区、洪山区青菱河路白沙洲农副产品大市场、白沙州大道张家湾社区及周边地区、天河机场、汇丰企业天地、青山公园、北湖项目部、严东湖、玖通达物流园、中央水景公园)进行分析,具体位置如图1中红色圆点所示,各特征点沉降序列结果如表3所示。2016—2017年期间,汇丰企业、武商、后湖公园累计沉降量呈稳定波动;2017—2018年期间,除青山公园,汇丰企业园累计沉降量呈稳定波动外,其他特征点累计沉降量呈现波动下降趋势,其中后湖公园累计沉降量呈断崖式下降波动趋势;2018—2020年期间,所选特征点累计沉降量呈波动下降趋势,白沙洲市场特征点下降最为严重,在整个研究期间持续下沉;2020—2021年期间,除白沙洲市场特征点、张家湾特征点累计沉降量呈波动下降趋势外,其他特征点累计沉降量呈稳定波动。在2016—2021年期间,白沙市场特征点下降最为突出,累计沉降量达到-160 mm,与其同为洪山区的张家湾特征点累计沉降量也达到-80 mm,两地的沉降速率分别达-29.26 mm/a和-15.28 mm/a。总体而言,选择的沉降区特征点均具有明显的下降趋势,且下降过程伴随有一定的波动特征,其中下沉趋势在不同年间会有区别,2016—2020年间多数特征点下沉速度较慢,2018年前后开始加速下沉,至2020年下沉速度开始放缓。
表3 特征点位置及沉降信息
地质构造背景、地下水状况、地面荷载、地下空间开发、工程施工等因素均会对地面沉降产生影响。武汉地区存在碳酸盐岩条带,高危险区为第Ⅰ型地质结构分布区主要位于:白沙洲条带(L3)和汉南条带(L6)长江两岸全新统地层覆盖区,粉细砂直接位于碳酸盐岩地层之上。本区约占碳酸盐岩分布面积的3.6%。由于全新统粉细砂的凝聚力很小或等于零,很容易在重力和渗透力作用下进入下伏的岩溶通道中,从而产生岩溶地面塌陷。初步推测白沙洲大市场、张家湾社区所发生的地面沉降与白沙洲条带(L3)南岸地质构造碳酸盐岩条带有一定的相关性。
2.3 主城区地面沉降自相关与冷热点分析
利用Arc GIS空间分析工具,对主城区地面沉降空间集聚特征进行分析。武汉市主城区2016—2021年间地面沉降速率的莫兰指数值是0.763 3,通过了显著性检验(Z>1.96,P<0.05)。空间自相关分析结果表明2016—2021年武汉市主城区地面沉降存在空间正相关性,空间集聚分布明显。进一步对武汉主城区地面沉降进行冷热点分析,结果如图1所示。
图1 2016—2021年研究地区地面沉降热点分布图
武汉市主城区地面沉降热点区域主要是青山区、江岸区、江汉区、武昌区沙湖公园部分区域。2016—2021年武汉市主城区地面沉降热点区域面积较大,地面沉降呈连片化发展。汉口区域的地面沉降热点区域遍布汉口城区,沙湖北部的地面沉降区域也有连片趋势,青山区及江岸区三环线地区沉降明显增多。同时在洪山区白沙洲大道及汉阳四新地区沉降分布也具有连片化特征。
2.4 In SAR处理成果与GNSS站点结果的精度评价
为了对GNSS和In SAR两项监测技术进行相互验证,同时也为了更好地研究武汉地区的地表变形特征,本文进一步比较了GNSS监测点的垂直形变与对应In SAR IPT点结果的一致性,如表4所示。从表中可以看出,所有站点的GNSS与In SAR监测结果的互差均不大于5 mm/a,互差绝对值的平均值为2.35 mm/a,说明两者整体上具有很好的一致性。
通过功率谱分析可知,周年和半周年变化是GNSS基准站主要的周期特征,同时U方向的振幅要显著大于水平方向。垂向的周年振幅范围为1.63 (W005)至8.66 mm (W003),振幅平均值4.17 mm;半周年振幅相对较小,振幅范围为0.64 (W001)至1.72 mm(W005),平均振幅1.03 mm;相对于周年振幅而言,半周年振幅仅为周年振幅的25%。南北方向的周年振幅范围为0.40 (W006)至2.36 mm (W002),振幅平均值1.10 mm;半周年振幅相对较小,振幅范围为0.34(W001)至1.04 mm (W002),平均振幅0.56 mm;半周年振幅为周年振幅的51%,GNSS基准站所在地整体处于相对稳定的状态。
从GNSS基准站坐标时间序列与In SAR形变序列的对比可以看出(为了更好地比较2种监测技术得到形变结果的区别,2种序列进行了标准化处理),GNSS结果除了反映趋势项,还能反映出周期性的变化特征,In SAR相比较而言对于速度变化的结果可以较好地反映出来,但是周期特征不显著。
表4 GNSS和In SAR监测成果比较/(mm/a)
3、结语
本文利用IPTA技术对武汉地区2016—2021年间89景哨兵一号SAR数据进行了处理,获取了研究区域的地面沉降信息,并联合GNSS数据对结果进行了验证。结果表明从沉降空间分布上看,地面沉降速率较高(<-5 mm/a)的区域主要是长江左岸的汉口地区和长江右岸的沙湖北部地区、白沙洲大道附近地区。通过冷热点分析可知,武汉市主城区地面沉降热点区域主要是青山区、江岸区、江汉区、武昌区沙湖公园部分区域,2016—2021年武汉市主城区地面沉降热点区域面积较大,地面沉降呈连片化发展。通过沉降特征点时间序列变化分析表明,选择的沉降区特征点均具有明显的下降趋势,且下降过程伴随有一定的波动特征,其中下沉趋势在不同年间会有区别,2016—2020年间多数特征点下沉速度较慢,2018年前后开始加速下沉,至2020年下沉速度开始放缓。最后联合GNSS监测点的垂直形变对In SAR结果进行了验证,所有站点的GNSS与In SAR监测结果的互差均不大于5 mm/a,互差绝对值的平均值为2.35 mm/a,说明两者整体上具有很好的一致性。
参考文献:
[1]袁铭,白俊武,秦永宽.国内外地面沉降研究综述[J].苏州科技学院学报(自然科学版),2016,33(1):1-5
[2]刘国祥,陈强,罗小军.In SAR原理与应用[M].北京:科学出版社,2019
[5]潘超,江利明,孙奇石,等.基于Sentinel-1雷达影像的成都市地面沉降In SAR监测分析[J].大地测量与地球动力学,2020,40(2):198-203
[6]许才军,何平,温扬茂,等.In SAR技术及应用研究进展[J].测绘地理信息,2015,40(2):1-9
[7]张力彪,闫若鹏.基于IPTA-In SAR技术的太原市地面形变监测[J].山西建筑,2023,49(1):4-6
[8]王心雨.高分辨率数据用于西安市表面位移的In SAR监测研究[D].西安:长安大学,2015
文章来源:魏以宽,闵天,邹崇尧等.基于InSAR的武汉地区2016—2021年地面沉降监测[J].地理空间信息,2023,21(10):69-72.
分享:
随着我国交通基础设施的发展,我国道路建设的速度和规模早已居于世界第一[1]。然而,由于道路建设施工质量不佳、建设工期和线路选择不合理以及道路基础地质条件复杂等因素,造成近几年我国城市道路塌陷事故频发,严重威胁人民群众的生命安全和城市交通秩序[2]。
2024-01-05晚二叠世作为全球地质历史时期重要的成煤—聚煤期[1],在我国华南地区也形成了大量的可采煤层,尤以西南各省煤炭资源最为丰富,滇东、黔西、川南等含煤区已成为我国南方主要的煤炭生产基地[2]。滇东含煤区晚二叠世所形成的1套易受海进—海退变化影响的近海型含煤岩系,呈现出陆相→海陆交互(过渡)相→海相的沉积相变特征,为分析该时期沉积环境变迁和聚煤规律提供了较为理想的地层基础[3,4]。
2024-01-03中国高速铁路地震预警系统由高速铁路地震预警监测系统和车载地震紧急处置装置组成[1,2],是实现高速铁路沿线地震实时监测和地震台网信息接入、进行综合分析处理后生成及传输并发布地震警报信息和紧急处置信息、通过车地联动的方式对列车采取紧急处置措施的系统。高速铁路地震预警监测系统联调联试(简称地震联调联试)是通过对地震预警监控单元、铁路局集团公司中心系统及其与相关系统接口的功能以及整体性能的测试[3]。验证系统是否满足相关标准规定和设计要求,依据测试结果指导动态调试,为动态验收提供依据。
2023-12-29地质矿产是自然资源的重要组成部分,地质是指地球的成分构造和土壤的质地;矿产泛指一切埋藏在地下,可供人类利用的天然矿物或岩石资源[1]。《自然资源部信息化建设总体方案》提出:立足已有基础,统筹整合土地、地质、矿产、海洋、测绘地理信息的信息化资源,运用新一代信息技术,通过完善、优化和创新,建设自然资源“一张图”,但目前地质矿产信息化建设严重不足,大量地质矿产数据还采用纸质人工管理,即便搭建了系统也主要为了满足地质档案馆的数字化要求,未实现网络服务和自然资源管理融合[2]。
2023-12-22水害的防治一直是煤矿安全生产工作中的要点,需要矿企结合《煤矿防治水细则》中的要求,每三年进行矿井水文地质类型的修定,实现对水害事故的控制,降低经济、人员等方面的损失。在实际的水文地质分析中,要对矿井已有的地质及水文信息资料进行全面收集,有针对的设计地质勘测及调查工作,还要结合新的勘查数据进行整理分析,评价矿井当前的水文地质条件。
2023-12-06根据弹性回跳理论分析断层内部的应变积累、闭锁程度和滑动亏损分布,是研究地震的重要手段。根据现代大地测量资料反演断层深浅部的闭锁程度和滑动亏损分布特征是判断断裂带未来地震危险性的重要手段[1]。基于Okada模型[2]的Defnode负位错反演程序[3,4],可利用GPS、水准、地质数据等结果,采用网格搜索和模拟退火方法求出块体旋转运动的欧拉极、块体边界断层的闭锁程度和滑动亏损速率。
2023-11-30油田水是石油和天然气勘探开发过程中产生的最大废水,其来源较为复杂,包括天然存在的深层地下水以及在油气开采过程中用于储层处理、压力支撑和水力压裂的注入水[1]。随着油田开发的不断深入,油田水含量不断增加。近年来我国大部分油田已进入高产水期,一般综合含水率在80%以上,每年产生的油田水超过5×108t[2]。油田水组成十分复杂,包括溶解态以及颗粒态的有机及无机化合物、天然放射性物质、挥发性有机物、溶解性气体、微生物以及为了提高产油/气量而人为添加的表面活性剂、润滑剂、杀菌剂、除垢剂等[3,4]。
2023-11-24我国海陆过渡相页岩气资源丰富,约占页岩气总资源量的1/4(匡立春等, 2020)。二叠系海陆过渡相富有机质泥页岩层系分布广泛,主要有以山西组、太原组和龙潭组为代表的含煤建造型和以孤峰组、大隆组为代表的含硅质建造型(翟刚毅等, 2020; 曹磊等, 2020)。
2023-11-24随着勘探开发的深入,发现储层酸化效果极差,油气开采遇阻,究其原因为绿泥石遇酸会产生大量氢氧化铁,堵塞孔隙,严重破坏油层孔隙结构;不同含量及赋存状态的绿泥石造成相邻井产能差异较大,难以对区块产能进行有效性评价。由于前人针对研究区齐古组绿泥石缺乏系统性研究,面对现场问题难以进行有效地质指导。针对上述问题,本次研究利用野外剖面、岩心观察等宏观分析手段和薄片鉴定、扫描电镜等微观分析手段,详细分析了阜东斜坡区齐古组砂岩储层中自生绿泥石赋存状态、分布特征以及形成机制,重点加强了优质储层中绿泥石含量及其赋存状态的分析,
2023-11-24环境同位素是地下水流动系统中查找水资源循环的重要方法,在天然地表及地下水体的离子之间作用中,其物质组成和成分是记录水体的转换、循环、流动、富集的重要标志。特别是天然水体中的氢氧稳定同位素完整的体现在水循环演化过程中,是展现地表水和地下水及地下水各含水系统中水岩相互作用及循环的有效标志和示踪剂,国内各流域范围内各水体之间的补给、排泄关系是目前重要环境科学研究内容之一,受到各方广泛重视。
2023-11-08我要评论
期刊名称:地质学刊
期刊人气:1266
主管单位:江苏省自然资源厅
主办单位:江苏省地质调查研究院,江苏省地质学会,中国地质学会
出版地方:江苏
专业分类:科学
国际刊号:1674-3636
国内刊号:32-1796/P
创刊时间:1977年
发行周期:季刊
期刊开本:16开
见刊时间:4-6个月
影响因子:2.546
影响因子:0.659
影响因子:0.149
影响因子:1.209
影响因子:2.866
400-069-1609
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!