留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

利用TPCA分析北京平原区地面沉降的时空演化特征

赵亚丽 王彦兵 王新雨 田秀秀 李小娟 余洁

赵亚丽, 王彦兵, 王新雨, 田秀秀, 李小娟, 余洁. 利用TPCA分析北京平原区地面沉降的时空演化特征[J]. 武汉大学学报 ● 信息科学版, 2022, 47(9): 1498-1506. doi: 10.13203/j.whugis20200721
引用本文: 赵亚丽, 王彦兵, 王新雨, 田秀秀, 李小娟, 余洁. 利用TPCA分析北京平原区地面沉降的时空演化特征[J]. 武汉大学学报 ● 信息科学版, 2022, 47(9): 1498-1506. doi: 10.13203/j.whugis20200721
ZHAO Yali, WANG Yanbing, WANG Xinyu, TIAN Xiuxiu, LI Xiaojuan, YU Jie. Temporal and Spatial Analysis of Land Subsidence in Beijing Plain Based on TPCA[J]. Geomatics and Information Science of Wuhan University, 2022, 47(9): 1498-1506. doi: 10.13203/j.whugis20200721
Citation: ZHAO Yali, WANG Yanbing, WANG Xinyu, TIAN Xiuxiu, LI Xiaojuan, YU Jie. Temporal and Spatial Analysis of Land Subsidence in Beijing Plain Based on TPCA[J]. Geomatics and Information Science of Wuhan University, 2022, 47(9): 1498-1506. doi: 10.13203/j.whugis20200721

利用TPCA分析北京平原区地面沉降的时空演化特征

doi: 10.13203/j.whugis20200721
基金项目: 

北京市自然科学基金 8202009

虚拟现实技术与系统国家重点实验室开放基金 01117220010020

详细信息
    作者简介:

    赵亚丽,硕士,主要从事地理信息科学与InSAR地面沉降监测及分析研究。zhaoyali@cnu.edu.cn

    通讯作者: 王彦兵,博士,副教授。wyb@cnu.edu.cn
  • 中图分类号: P208;P642

Temporal and Spatial Analysis of Land Subsidence in Beijing Plain Based on TPCA

Funds: 

The Beijing Natural Science Foundation 8202009

the Open Project Program of the State Key Laboratory of Virtual Reality Technology and Systems, Beihang University, China 01117220010020

More Information
    Author Bio:

    ZHAO Yali, master, specializes in geographic information science and land subsidence monitoring and analysis research based on InSAR. E-mail: zhaoyali@cnu.edu.cn

    Corresponding author: WANG Yanbing, PhD, associate professor. E-mail: wyb@cnu.edu.cn
  • 摘要: 时间主成分分析(temporal principal component analysis,TPCA)可用于地学领域中提取时空数据的时序特征和空间分布特征,北京平原区的地面沉降具有典型的时序和空间特征。在利用永久散射体干涉测量技术获取的北京平原区2003—2010年地面沉降数据的基础上,采用TPCA方法分析了北京平原区地面沉降时空演化特征。经分析发现:(1)TPCA分析得到的第一主成分反映了地面沉降在该长时序阶段的空间分布特征。(2)第二主成分得分为正的空间点与可压缩层厚度在130 m以上的区域在空间分布上有一致性和相关性。(3)在空间上,第一主成分为负值与第二主成分为正值的永久散射体点分布在年均沉降速率30 mm/a以上的严重沉降区域。严重沉降区具有明显的南北沉降分类现象和季节性差异,具体表现为:北部沉降区在春夏季节的沉降量大于秋冬季节;南部沉降区则与之相反。总之,基于时间主成分分析方法可分析得到研究区的地面沉降时空演化规律,为城市安全监测提供数据支撑。
  • 图  1  TPCA数据结构

    Figure  1.  Data Structure of TPCA

    图  2  TPCA分析过程中的矩阵

    Figure  2.  Matrix from TPCA

    图  3  TPCA分析地面沉降时空特征技术路线

    Figure  3.  Technical Route to Analyze Land Subsidence Characteristics Using TPCA

    图  4  地面沉降与水准数据的精度验证

    Figure  4.  Accuracy Verification of Land Subsidence and Lev‍eling Data

    图  5  PC1~PC3时间序列载荷图

    Figure  5.  Time Series Loading Diagram of PC1-‍PC3

    图  6  PC1得分与年均沉降速率关系

    Figure  6.  Relationship Between PC1 Score and Average Annual Subsidence Rate

    图  7  PC2与可压缩厚度的关系(大于130 m)

    Figure  7.  Relationship Between PC2 Component and Compressible Thickness(> 130 m)

    图  8  两个主成分之间的关系以及沉降区域分类图

    Figure  8.  Relationship Between Two Principal Components

    图  9  2004年和2009年的差分沉降量图

    Figure  9.  Seasonal Variation of Cumulative Land Subsidence in 2004 and 2009

    图  10  季节性形变

    Figure  10.  Seasonal Land Subsidence

    图  11  历年的沉降速率剖面线

    Figure  11.  Subsidence Rate Profile Over the Years

    表  1  地面沉降的主成分特征值

    Table  1.   Principal Component Eigenvalues of Land Subsid‍ence

    成分 特征值 解释占比/% 累计占比/%
    PC1 43.95 86.18 86.18
    PC2 4.42 8.66 94.84
    PC3 1.21 2.37 97.21
    下载: 导出CSV
  • [1] 杨艳, 贾三满, 王海刚. 北京平原区地面沉降现状及发展趋势分析[J]. 上海地质, 2010, 31(4): 23-28 doi:  10.3969/j.issn.2095-1329.2010.04.005

    Yang Yan, Jia Sanman, Wang Haigang. The Status and Development of Land Subsidence in Beijing Plain[J]. Shanghai Geology, 2010, 31(4): 23-28 doi:  10.3969/j.issn.2095-1329.2010.04.005
    [2] 刘凯斯. 北京地铁M1/M6沿线区地面沉降演化特征及风险评价[D]. 北京: 首都师范大学, 2018

    Liu Kaisi. Evolution Characteristics and Risk Assessment of Land Subsidence in the Area along Beijing Subway M1/M6[D]. Beijing: Capital Normal University, 2018
    [3] 段光耀, 刘欢欢, 宫辉力, 等. 京津城际铁路沿线不均匀地面沉降演化特征[J]. 武汉大学学报·信息科学版, 2017, 42(12): 1847-1853 doi:  10.13203/j.whugis20150537

    Duan Guangyao, Liu Huanhuan, Gong Huili, et al. Evolution Characteristics of Uneven Land Subsid‍ence Along Beijing-Tianjin Inter-City Railway[J]. Geomatics and Information Science of Wuhan University, 2017, 42(12): 1847-1853 doi:  10.13203/j.whugis20150537
    [4] 罗三明, 杜凯夫, 万文妮, 等. 利用PSInSAR方法反演大时空尺度地表沉降速率[J]. 武汉大学学报·信息科学版, 2014, 39(9): 1128-1134 doi:  10.13203/j.whugis20130670

    Luo Sanming, Du Kaifu, Wan Wenni, et al. Ground Subsidence Rate Inversion of Large Temporal and Spatial Scales Based on Extended PSInSAR Method[J]. Geomatics and Information Science of Wuhan University, 2014, 39(9): 1128-1134 doi:  10.13203/j.whugis20130670
    [5] 朱邦彦, 姚冯宇, 孙静雯, 等. 利用InSAR与地质数据综合分析南京河西地面沉降的演化特征和成因[J]. 武汉大学学报·信息科学版, 2020, 45(3): 442-450 doi:  10.13203/j.whugis20190081

    Zhu Bangyan, Yao Fengyu, Sun Jingwen, et al. Attribution Analysis on Land Subsidence Feature in Hexi Area of Nanjing by InSAR and Geological Data[J]. Geomatics and Information Science of Wuhan University, 2020, 45(3): 442-450 doi:  10.13203/j.whugis20190081
    [6] Guo L, Gong H L, Zhu F, et al. Analysis of the Spatiotemporal Variation in Land Subsidence on the Beijing Plain, China[J]. Remote Sensing, 2019, 11‍(10): 1170-1189 doi:  10.3390/rs11101170
    [7] Zhou C D, Lan H X, Gong H L, et al. Reduced Rate of Land Subsidence Since 2016 in Beijing, China: Evidence from Tomo-PSInSAR Using RadarSAT‍-‍2 and Sentinel‍-‍1 Datasets[J]. International Journal of Remote Sensing, 2020, 41(4): 1259-1285 doi:  10.1080/01431161.2019.1662967
    [8] Zuo J J, Gong H L, Chen B B, et al. Time-Series Evolution Patterns of Land Subsidence in the East‍ern Beijing Plain, China[J]. Remote Sensing, 2019, 11(5): 539 doi:  10.3390/rs11050539
    [9] Richman M B. Rotation of Principal Components[J]. Journal of Climatology, 1986, 6(3): 293-335 doi:  10.1002/joc.3370060305
    [10] Lin Y N N, Kositsky A P, Avouac J P. PCAIM Joint Inversion of InSAR and Ground‍-‍Based Geodet‍ic Time Series: Application to Monitoring Magmatic Inflation Beneath the Long Valley Caldera[J]. Geophysical Research Letters, 2010, 37(23): 23301-23305
    [11] Ji K H, Herring T A. Transient Signal Detection Using GPS Measurements: Transient Inflation at Akutan Volcano, Alaska, During Early 2008[J]. Geophysical Research Letters, 2011, 38(6): 6307-6312
    [12] Zhang J P, Zhu T, Zhang Q H, et al. The Impact of Circulation Patterns on Regional Transport Pathways and Air Quality over Beijing and Its Surround‍ings[J]. Atmospheric Chemistry and Physics, 2012, 12(11): 5031-5053
    [13] 朱飙, 王振会, 李春华, 等. 江苏雷暴时空变化的气候特征分析[J]. 气象科学, 2009, 29(6): 849-852 doi:  10.3969/j.issn.1009-0827.2009.06.023

    Zhu Biao, Wang Zhenhui, Li Chunhua, et al. Anal‍ysis of Climate Spatial‍-‍Temporal Character of Thunderstorm over Jiangsu Province[J]. Scientia Meteorologica Sinica, 2009, 29(6): 849-852 doi:  10.3969/j.issn.1009-0827.2009.06.023
    [14] Neeti N, Eastman J R. Novel Approaches in Extended Principal Component Analysis to Compare Spatio‍-‍Temporal Patterns Among Multiple Image Time Series[J]. Remote Sensing of Environment, 2014, 148: 84-96 doi:  10.1016/j.rse.2014.03.015
    [15] Rudolph M L, Shirzaei M, Manga M, et al. Evolution and Future of the Lusi Mud Eruption Inferred from Ground Deformation[J]. Geophysical Research Letters, 2013, 40(6): 1089-1092 doi:  10.1002/grl.50189
    [16] Lipovsky B. Physical and Statistical Models in Deformation Geodesy [D]. Riverside, USA: University of California, Riverside, 2011
    [17] Chaussard E, Bürgmann R, Shirzaei M, et al. Predictability of Hydraulic Head Changes and Character‍ization of Aquifer‍-‍System and Fault Properties from InSAR-Derived Ground Deformation[J]. Journal of Geophysical Research: Solid Earth, 2014, 119(8): 6572-6590 doi:  10.1002/2014JB011266
    [18] 吴玉苗. 基于EOF与神经网络的隧道变形监测方法研究[D]. 成都: 西南交通大学, 2014

    Wu Yumiao. Investigation on Tunnel Deformation Monitoring Methods Based on the EOF and Neural Network[D]. Chengdu: Southwest Jiaotong University, 2014
    [19] 邹正波, 李辉, 吴云龙, 等. 日本Mw 9.0地震震区及其周缘2002-2015年卫星重力变化时空特征[J]. 地震学报, 2016, 38(3): 417-428 https://www.cnki.com.cn/Article/CJFDTOTAL-DZXB201603009.htm

    Zou Zhengbo, Li Hui, Wu Yunlong, et al. Spatial and Temporal Characteristics of Long‍-‍Term Satellite Gravity Change in the Epicenter of Mw 9.0 Japan Earthquake and Its Surrounding Regions[J]. Acta Seismologica Sinica, 2016, 38(3): 417-428 https://www.cnki.com.cn/Article/CJFDTOTAL-DZXB201603009.htm
    [20] Jiang L, Bai L, Zhao Y, et al. Combining InSAR and Hydraulic Head Measurements to Estimate Aquifer Parameters and Storage Variations of Confined Aquifer System in Cangzhou, North China Plain[J]. Water Resources Research, 2018, 54(10): 8234-8252
  • [1] 杨宏山, 聂建亮.  陕西地区垂直形变分析 . 武汉大学学报 ● 信息科学版, 2022, 47(11): 1854-1859. doi: 10.13203/j.whugis20200149
    [2] 张扬, 刘艳芳, 刘莹, 刘耀林, 陈雨露, 王征禹.  武汉市地面沉降时空分异特征及地理探测机制 . 武汉大学学报 ● 信息科学版, 2022, 47(9): 1486-1497. doi: 10.13203/j.whugis20210143
    [3] 谢萍, 张双喜, 周吕, 李庆隆, 肖家豪, 蔡剑锋.  武汉市中心城区地表形变与洪涝灾害防治新策略 . 武汉大学学报 ● 信息科学版, 2021, 46(7): 1015-1024. doi: 10.13203/j.whugis20190439
    [4] 许强, 蒲川豪, 赵宽耀, 何攀, 张含悦, 刘佳良.  延安新区地面沉降时空演化特征时序InSAR监测与分析 . 武汉大学学报 ● 信息科学版, 2021, 46(7): 957-969. doi: 10.13203/j.whugis20200146
    [5] 朱邦彦, 姚冯宇, 孙静雯, 王晓.  利用InSAR与地质数据综合分析南京河西地面沉降的演化特征和成因 . 武汉大学学报 ● 信息科学版, 2020, 45(3): 442-450. doi: 10.13203/j.whugis20190081
    [6] 蒲川豪, 许强, 蒋亚楠, 赵宽耀, 何攀, 张含悦, 李骅锦.  延安新区地面沉降分布及影响因素的时序InSAR监测分析 . 武汉大学学报 ● 信息科学版, 2020, 45(11): 1728-1738. doi: 10.13203/j.whugis20190372
    [7] 杨梦诗, 廖明生, 秦晓琼, 史绪国.  C和L波段SAR数据在填海新区的应用及特性分析 . 武汉大学学报 ● 信息科学版, 2017, 42(9): 1300-1305. doi: 10.13203/j.whugis20150356
    [8] 袁鹏, 孙宏飞, 秦昌威, 张理想.  安徽CORS参考站三维速度场分析 . 武汉大学学报 ● 信息科学版, 2016, 41(4): 535-540. doi: 10.13203/j.whugis20140552
    [9] 张金芝, 黄海军, 毕海波, 王权.  SBAS时序分析技术监测现代黄河三角洲地面沉降 . 武汉大学学报 ● 信息科学版, 2016, 41(2): 242-248. doi: 10.13203/j.whugis20140067
    [10] 钱 闯, 刘 晖, 丁志刚, 郑洪艳.  顾及非构造形变的参考站长期稳定性分析 . 武汉大学学报 ● 信息科学版, 2015, 40(9): 1259-1265. doi: 10.13203/j .whu g is20130715
    [11] 杨成生, 张勤, 赵超英, 季灵运.  短基线集InSAR技术用于大同盆地地面沉降、地裂缝及断裂活动监测 . 武汉大学学报 ● 信息科学版, 2014, 39(8): 945-950. doi: 10.13203/j.whugis20130656
    [12] 李永生, 张景发, 李振洪, 罗毅.  利用短基线集干涉测量时序分析方法监测北京市地面沉降 . 武汉大学学报 ● 信息科学版, 2013, 38(11): 1374-1377.
    [13] 李永生, 张景发, 罗毅, 姜文亮.  利用高分辨率聚束模式TerraSAR-X影像的PSInSAR监测地表变形 . 武汉大学学报 ● 信息科学版, 2012, 37(12): 1452-1455.
    [14] 何敏, 何秀凤.  利用时间序列干涉图叠加法监测江苏盐城地区地面沉降 . 武汉大学学报 ● 信息科学版, 2011, 36(12): 1461-1465.
    [15] 张诗玉, 李陶, 夏耶.  基于InSAR技术的城市地面沉降灾害监测研究 . 武汉大学学报 ● 信息科学版, 2008, 33(8): 850-853.
    [16] 董玉森, GeLinlin, ChangHsingchun, 张志.  基于差分雷达干涉测量的矿区地面沉降监测研究 . 武汉大学学报 ● 信息科学版, 2007, 32(10): 888-891.
    [17] 董国凤, 张蕾, 赵全, 赵新华.  天津市塘沽区地面沉降监测网络优化 . 武汉大学学报 ● 信息科学版, 2006, 31(11): 1015-1017.
    [18] 周东旭, 冯义楷, 张化疑, 付延光, 唐秋华.  联合卫星测高和GNSS观测的天津沿海近25年相对海平面变化分析 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20210532
    [19] 龚翔, 许强, 蒲川豪, 陈婉琳, 修德皓, 纪续.  兰州新区平山造城工后地表形变InSAR时序监测与分析 . 武汉大学学报 ● 信息科学版, 0, 0(0): -. doi: 10.13203/j.whugis20210553
    [20] 张双成, 张雅斐, 司锦钊, 罗勇, 余静, 雷坤超, 许强.  南水进京后升降轨InSAR解译北京地面沉降发展态势 . 武汉大学学报 ● 信息科学版, 0, 0(0): -. doi: 10.13203/j.whugis20210554
  • 加载中
图(11) / 表(1)
计量
  • 文章访问数:  499
  • HTML全文浏览量:  132
  • PDF下载量:  79
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-03-04
  • 刊出日期:  2022-09-05

利用TPCA分析北京平原区地面沉降的时空演化特征

doi: 10.13203/j.whugis20200721
    基金项目:

    北京市自然科学基金 8202009

    虚拟现实技术与系统国家重点实验室开放基金 01117220010020

    作者简介:

    赵亚丽,硕士,主要从事地理信息科学与InSAR地面沉降监测及分析研究。zhaoyali@cnu.edu.cn

    通讯作者: 王彦兵,博士,副教授。wyb@cnu.edu.cn
  • 中图分类号: P208;P642

摘要: 时间主成分分析(temporal principal component analysis,TPCA)可用于地学领域中提取时空数据的时序特征和空间分布特征,北京平原区的地面沉降具有典型的时序和空间特征。在利用永久散射体干涉测量技术获取的北京平原区2003—2010年地面沉降数据的基础上,采用TPCA方法分析了北京平原区地面沉降时空演化特征。经分析发现:(1)TPCA分析得到的第一主成分反映了地面沉降在该长时序阶段的空间分布特征。(2)第二主成分得分为正的空间点与可压缩层厚度在130 m以上的区域在空间分布上有一致性和相关性。(3)在空间上,第一主成分为负值与第二主成分为正值的永久散射体点分布在年均沉降速率30 mm/a以上的严重沉降区域。严重沉降区具有明显的南北沉降分类现象和季节性差异,具体表现为:北部沉降区在春夏季节的沉降量大于秋冬季节;南部沉降区则与之相反。总之,基于时间主成分分析方法可分析得到研究区的地面沉降时空演化规律,为城市安全监测提供数据支撑。

English Abstract

赵亚丽, 王彦兵, 王新雨, 田秀秀, 李小娟, 余洁. 利用TPCA分析北京平原区地面沉降的时空演化特征[J]. 武汉大学学报 ● 信息科学版, 2022, 47(9): 1498-1506. doi: 10.13203/j.whugis20200721
引用本文: 赵亚丽, 王彦兵, 王新雨, 田秀秀, 李小娟, 余洁. 利用TPCA分析北京平原区地面沉降的时空演化特征[J]. 武汉大学学报 ● 信息科学版, 2022, 47(9): 1498-1506. doi: 10.13203/j.whugis20200721
ZHAO Yali, WANG Yanbing, WANG Xinyu, TIAN Xiuxiu, LI Xiaojuan, YU Jie. Temporal and Spatial Analysis of Land Subsidence in Beijing Plain Based on TPCA[J]. Geomatics and Information Science of Wuhan University, 2022, 47(9): 1498-1506. doi: 10.13203/j.whugis20200721
Citation: ZHAO Yali, WANG Yanbing, WANG Xinyu, TIAN Xiuxiu, LI Xiaojuan, YU Jie. Temporal and Spatial Analysis of Land Subsidence in Beijing Plain Based on TPCA[J]. Geomatics and Information Science of Wuhan University, 2022, 47(9): 1498-1506. doi: 10.13203/j.whugis20200721
  • 北京地区最早发生的地面沉降出现在20世纪30年代,位于西单-东单区域。近年来,北京的快速发展需求使得地下水长期超采,导致地面沉降的幅度和范围逐年扩大。2003—2010年的最大年沉降速率达到110 mm/a,最大累计沉降量达到了723 mm,年均沉降速率达到30 mm/a的覆盖区域面积为480 000 000 m2。北京平原区形成海淀苏家坨、昌平沙河-八仙庄、顺义、朝阳来广营、东郊八里庄、大兴榆垡6大沉降区[1]。由于地面沉降严重威胁城市安全,因此需要分析该地区的地面沉降时空演化特征,预测演化趋势。

    有关地面沉降时空演化特征方面的研究方法分为时序分析与空间分析两类。时序演化特征分析通常采用典型点时序图法,从原始数据中观察地面沉降在时间序列上的变化特征或者通过统计某区域的年度沉降量来分析时序特征。刘凯斯[2]采用时序排列熵法分析北京地铁1号和6号线的地面沉降时间序列演化特征。段光耀等[3]用Mann-Kendall检验对北京平原区时空变化进行分析,研究历年发生突变现象的机理。空间演化特征分析一般采用剖面分析、梯度分析来分析空间上的不均匀地面沉降特征[4-6]。Zhou等[7]用空间分析方法等扇分析探究2012—2018年北京平原区地面沉降的扩张趋势,发现地面沉降由向东扩张变为向东和向北扩张。Zuo等[8]运用标准差椭圆方法发现北京市地面沉降漏斗的移动,揭示不均匀地面沉降。上述方法分别从时间或空间的角度研究沉降演化特征,在时间和空间上相分离,不能从时空角度发现数据中隐藏的信息和可能存在的规律。本文采用高维数据分析中的主成分分析(principal component analysis,PCA)方法研究地面沉降的时空演化特征,充分利用合成孔径雷达干涉测量(interferometric synthetic aperture radar,InSAR)所得的沉降信息具有长时序、覆盖范围广的优势。

    PCA常用于数据降维,本文的应用目的是通过降维挖掘地面沉降的主要时空特征。PCA应用于地学领域中,能够有效地从时空数据中提取出某信号的时间序列与空间分布,主成分分析模式分为6种,其中T模式(temporal mode)[9]已应用在GPS站点数据、电磁测距和潮汐计数据上,用来分离出瞬态形变事件[10-11]。气象领域中,可对得出的多种环流模式进行解释[12],分析气象雷暴日的规律[13]

    T模式时间序列分析能识别多个时间序列之间的相似空间模式[14]。Rudolph等[15]通过时间主成分分析(temporal principal component analysis,TPCA)从InSAR时序数据中提取主要的时间行为模式;Lipovsky[16]运用TPCA提取长时序形变的季节信号;Chaussard等[17]对小范围且量级小的InSAR监测结果进行TPCA分析,得到第一成分为沉降主趋势,第二主成分表征为季节性形变,与承压水空间覆盖范围较一致,第三主成分表现的空间特征与断裂带位置相关;吴玉苗[18]用类似于PCA的经验正交函数得到隧道两个方向的变形时空特征;邹正波等[19]基于重力场数据识别日本地震,并研究2002—2015年的重力场时空变化特征;Jiang等[20]在对沧州中部承压含水层系统的含水层参数和地下水储量变化定量研究过程中,通过多通道奇异谱分析对地表形变和地下水数据中的季节信号进行分离,推算出弹性骨架存储率,此外,分别重构出总地下水储量、可恢复地下水储量和不可逆地下水储量。

    综上所述,TPCA可应用于地学领域中,在未知先验知识条件下,提取时空数据中的时间序列和空间分布特征。本文使用TPCA方法来分析2003—2010年的北京平原区地面沉降,定量提取时空特征并进行合理解释。

    • TPCA是一种处理多因素复杂问题的多元统计方法,其主要用途就是进行数据降维和数据解释。作为统计分析“由表及里”的数学手段,常用于指标评价研究中,近些年也应用于地理数据的模式识别与数据挖掘中。

      利用合成孔径雷达(synthetic aperture radar,SAR)数据获取的地面沉降数据具有大量的空间点位、长时间的监测序列等特点。本文将TPCA方法应用到高维度的地面沉降数据中,可理解为对众多观测点的多维时间数据集以时间为变量进行PCA分析,得到数据本身所蕴含的时空特征和规律。

    • 通过TPCA构建M×N的矩阵X,如图 1(a)‍所示,M指观测点个数,代表永久散射体干涉测量技术(permanent scatterer interferometric synthetic aperture radar,PS-InSAR)识别的永久散射体(permanent scatterer,PS)点;N指日期,代表SAR影像的获取时间。本文的数据即图 1(b)中每一个地理位置的PS点,共近10万个点;每个PS点有51个监测数据。

      图  1  TPCA数据结构

      Figure 1.  Data Structure of TPCA

      XM个样本、N个变量构成的原始矩阵,其表达式为:

      X=X1,X2XN=x11x12x1Nx21x22x2NxM1xM2xMN

      i个变量可表示为Xi=(x1ix2ixMiT,通过TPCA对X线性变换后,得到新变量矩阵Z

      Z1=u11X1+u21X2++uN1XNZ2=u12X1+u22X2++uN2XN ZN=u1NX1+u2NX2++uNNXN

      其中,

      U=U1,U2UN=u11u12u1Nu21u22u2NuN1uN2uNN

      为了用Z1代替原来的N个变量X1X2XN,要求Z1尽可能多地反映原来N个变量的信息,信息量大小用方差来描述。方差Var(Z1)越大,表示Z1包含的信息越多。

      当变换矩阵满足UU'=I,存在满足条件的U1使得Var(Z1)达到最大,Z1就称为第一主成分(或者主分量);使得方差Var(Z2)达到第二最大,则Z2成为第二主成分,以此类推求得剩余成分。具体过程如下:(1)对矩阵X均值化:求矩阵X的均值E‍(X)‍=μX减去每一时间的平均值;(2)对均值化的矩阵进行标准化,得到标准化矩阵S;(3)基于标准化矩阵求得相关系数矩阵R为:

      R=S'S/(M-1)

      求式(2)中Z1=XU1,也就是求特征向量U1‍=(u11u21uN1',使U1U1=I时,Var(Z1)达到最大;等价于求相关矩阵R的特征值和特征向量。设λ=λ1R的最大特征值,对应的特征向量即为所求U1X的第i主成分得分Zi可通过R的第i特征向量Ui得到,即:

      Zi=XUi(i=1, 2N)

      式中,XM×N原始矩阵;UN×N特征向量矩阵,也称为载荷矩阵;Zi是新变量得分矩阵,称为主成分得分或者分量,代表原始数据X的线性组合,图 2为TPCA分析过程中的矩阵。

      图  2  TPCA分析过程中的矩阵

      Figure 2.  Matrix from TPCA

    • TPCA解算可以得到特征值、特征向量、主成分得分3个变量,其意义分别为:(1)特征值是矩阵的一种数学特性,用于决定提取分量的个数,测量主成分分量所解释的方差大小。特征值越大,主成分分量的方差越大,分量解释的信息占比(即方差贡献率=主成分特征值/特征值总和)越多。(2)特征向量是特征值对应的向量。时间主成分分析的特征向量可以表征主成分的时间特征,从而揭示可能存在的形变规律。(3)依据方差大小选择特征向量即主分量方向,将原始数据投影到主分量方向上,得到主成分得分。表征每个PS点的主要信息,这里主要是指PS点的空间分布特征。经TPCA分解出的特征向量和主成分得分能够反映数据变化的时间空间结构,按照特征值大小排列,量化可解释原信息的大小。从定义上来看,PCA目的是线性组合找到新变量代替原变量即时间变量,因而TPCA的新变量侧重描述空间特征。

    • 本文以北京平原区为研究对象,技术路线如下:(1)借助Sarproz平台,对2003—2010年的Envisat ASAR数据及辅助数据数字高程模型(digital elevation model,DEM)进行PS-InSAR处理,经过预处理、生成干涉图、解算形变量,获取北京平原区地面沉降信息,得到约10万个PS点及其时序沉降结果。(2)将沉降结果与水准测量结果进行对比验证。(3)精度验证合格后,利用TPCA方法分析并提取北京平原区地面沉降的时空演化特征,如图3所示。

      图  3  TPCA分析地面沉降时空特征技术路线

      Figure 3.  Technical Route to Analyze Land Subsidence Characteristics Using TPCA

    • 从Envisat ASAR数据集的地面沉降监测结果中提取水准点附近的监测点,将水准点附近的InSAR地面沉降数据与2005—2010年的水准点监测数据进行对比验证(见图 4),发现两者年均沉降速率的相关系数达到0.97。以水准监测值为基准,PS-InSAR监测值相对于水准监测值的平均绝对值误差为3.5 mm/a,在轻微沉降区(年均沉降量小于10 mm/a),二者相差在2.5 ‍mm/‍a以内;在一般沉降区和严重沉降区(年均沉降量大于10 ‍mm/a),二者相差10 mm/a以内。由此可知,2003—2010年的沉降结果精度相对可靠。

      图  4  地面沉降与水准数据的精度验证

      Figure 4.  Accuracy Verification of Land Subsidence and Lev‍eling Data

    • 对10万个永久散射体的时序累计沉降数据进行时间主成分分析,得到2003—2010年北京平原区地面沉降的主成分特征值(见表 1)和主成分特征向量(见图 5)。除了第1、第2、第3成分,其他各成分的方差贡献率不到3%,因此本文取前3个成分作为主成分进行研究。

      表 1  地面沉降的主成分特征值

      Table 1.  Principal Component Eigenvalues of Land Subsid‍ence

      成分 特征值 解释占比/% 累计占比/%
      PC1 43.95 86.18 86.18
      PC2 4.42 8.66 94.84
      PC3 1.21 2.37 97.21

      图  5  PC1~PC3时间序列载荷图

      Figure 5.  Time Series Loading Diagram of PC1-‍PC3

      特征值越大,对应主成分分量解释的信息占比越多。由表 1可知,PC1特征值的方差贡献率达到86.18%,提取了原数据的绝大部分信息和变化特征,说明PC1代表北京平原区的主要地面沉降变化特征。

      由于主成分Zj=u1jX1+u2jX2+…+uNjXN,cov‍(XiZj)=cov(ui1Z1+ui2Z2+…+uiNZNZj)=uijλj,由此可得Xi(所有原变量)与Zj(某一主成分得分)的相关系数为:

      ρXi,Zj=covXi,Zjσiλj=uijλjσi

      式中,σiXi方差;λjZj方差。

      可见,第j个主成分得分Zj与变量Xi的相关程度与对应主成分线性组合系数(特征向量)的大小有关,特征向量代表了主成分得分与原始变量的相关程度。

      特征向量表示线性组合后的主成分得分与原始变量的相关程度,也表示主成分特征的时间变化趋势。图 5中PC1的特征向量在0.15附近保持稳定,说明地面沉降的发展趋势在该段时间内保持一致。PC2和PC3的特征向量整体变化幅度相对PC1要大一些,PC2和PC3主成分通过计算可以揭示地面沉降的部分季节性变化特征。

      2005年之前,前3个特征向量的系数有明显的变化幅度,这是因为2003—2005年的采集时间稀疏,处理数据是累计沉降数据,前段时间的差异变化相对大,后期细微变化表现得不明显,导致前段时间的系数变化较大。

    • 特征向量代表最大的变化方向。将标准化数据投影到特征向量上,得到的投影值是新变量,即主成分得分。不同主成分得分代表不同的空间分布特征。PC1得分揭示地面沉降的空间变化主趋势。由图 6所知,PC1得分(图 6(b))与年平均沉降速率(图 6(a))空间分布一致,均体现出北京平原区的主要沉降漏斗和不均匀沉降,且相关系数达0.98(图 6(c)),PC1方向上,地面沉降速率分区明显(图 6(d));PC1特征向量的系数差异不大,变化趋势稳定;PC1的方差贡献率最大,代表最大的变化特征。PC1表示地面沉降空间特征在2003‍—2010年一直发展,不均匀沉降显著。

      图  6  PC1得分与年均沉降速率关系

      Figure 6.  Relationship Between PC1 Score and Average Annual Subsidence Rate

    • PC2与可压缩层厚度相关性较强。观察PC2得分的空间分布(图 7(b)),黑色虚线范围为-‍30 ‍mm/a沉降漏斗区(下文同),发现平原区北部沉降区和东部漏斗沉降区存在不同的空间特征。PC2得分的正值分布范围与130 m以深的可压缩层空间分布相似(图 7(a));且两者相关系数达0.61,相关性较大(图 7(c));PC1和PC2得分的散点图中,PC2方向上,较厚(> 130 m)与较薄(60~130 m)的PS点有较为明显的分区现象(图 7(d)),说明可压缩层厚度在PC2方向上差异明显。因此PC2与可压缩层厚度相关性较强。

      图  7  PC2与可压缩厚度的关系(大于130 m)

      Figure 7.  Relationship Between PC2 Component and Compressible Thickness(> 130 m)

    • PC1得分绝对值远大于PC2,将两者叠加会使PC2的空间特征被PC1的空间主趋势掩盖。通过比较两分量的相位关系,研究两者的变化关系,发现PC1得分为负、且PC2为正的PS点覆盖了-30 mm/a沉降漏斗范围的严重沉降区域,如图 8(a)的Ⅱ象限。其余两分量的关系分别为:图 8(a)的I象限代表PC1、PC2均为正值,覆盖了非沉降区域;III象限代表PC1、PC2均为负值,覆盖一部分沉降区域和非沉降区域;IV象限代表PC1为正值,PC2为负值,覆盖一部分非沉降区域。

      图  8  两个主成分之间的关系以及沉降区域分类图

      Figure 8.  Relationship Between Two Principal Components

      由于严重沉降区域对经济社会生活有重大影响,因此着重分析严重沉降区域PS点的规律。通过空间自相关分析发现严重沉降区域的分类情况。对图 8(a)II象限的监测点的PC1得分进行空间自相关分析,得到该PS点的聚集分类图,如图‍图 8(b)所示。研究区的北部沉降区域(海淀苏家坨、昌平沙河-八仙庄、顺义)呈现高高集聚(见图 8(b)橙色部分),南部沉降区域(朝阳来广营、金盏-东坝、黑庄户)呈现低低集聚(见图 8(b)绿色部分)。

    • 研究区北部和南部沉降区域体现出不同的季节性沉降。通过统计一年内前后半年的累计沉降量来分析季节性沉降,将作差之后的沉降量(即春夏季节沉降量-秋冬沉降量)称为半年差分沉降量。半年差分沉降量为负值,则春夏季节沉降量多于秋冬季节沉降量;为正值,则秋冬季节沉降量多于春夏季节沉降量。从半年差分沉降量(见图 9)的统计中发现,北部沉降区域的春夏季节沉降量较多,南部沉降区域的秋冬季节沉降量较多。严重沉降区域的高高集聚、低低集聚分类图(见图 8(b))与严重沉降区域北部、南部季节性沉降的空间分布相似(见图 9)。

      图  9  2004年和2009年的差分沉降量图

      Figure 9.  Seasonal Variation of Cumulative Land Subsidence in 2004 and 2009

      进一步分析研究区北部沉降区和南部沉降区的差异。沿北部沉降区域的红虚线AA'、南部沉降区域的黑实线BB'(见图 9(a))做剖面分析,得到沿线2004—2009年的半年地面沉降量剖面图(见图10)。由图 9图 10发现:(1)北部沉降区域整体上在春夏季节的沉降更多。2005年数据缺失严重,未体现春夏季节的沉降多;2006年、2008年、2009年均有部分季节性形变不明显的地方,其沉降速率相较前一年变化大(图 11(a)),导致季节性不明显,沉降速率变化大可能与该位置地下水的水位异常变化有关系。(2)南部沉降区域整体上秋冬季节沉降更多。2005年同样未体现秋冬季节的沉降多;2004年、2006年、2008年、2009年,范围约100~8 ‍000‍ m部分对应来广营沉降漏斗,该沉降漏斗区别于其他两个沉降漏斗,可能与其地下水开采有较大差异有关;2007年相对前一年的沉降速率变化最大(见图 11(b)),其季节性沉降不明显。

      图  10  季节性形变

      Figure 10.  Seasonal Land Subsidence

      图  11  历年的沉降速率剖面线

      Figure 11.  Subsidence Rate Profile Over the Years

    • 本文利用TPCA对地面沉降数据进行分解,PS点作为观测点,时间作为变量。TPCA实质上是从时间维度上进行降维,将具有相同空间变化特征的观测点分离,线性组合的新变量代替原始的时间变量,得到空间特征。得到以下结论:(1)‍TPCA应用于具有空间位置和时序属性的数据集,可以识别出该空间下的主要特征。特征值决定了各分量解释的信息量大小。本文主成分特征向量为一个时间序列,代表PC空间模式与沉降之间的相关关系。(2)主成分得分是TPCA分解得到的空间模式,表示的是不同的空间特征;本文的PC1得分与年平均沉降速率空间分布一致,揭示地面沉降的空间变化主趋势;PC2得分与可压缩层厚度相关性较强。(3)从主成分得分关系中得出严重沉降区域的分类情况,主要体现南北沉降区域的不同聚类现象,该现象与季节性沉降相关。同时研究区的北部、南部沉降区域存在较明显的季节性差异,与前面所描述的集聚分类空间分布相似。

      TPCA属于数据驱动方法,不需要先验模型做出假设,主要利用多元统计原理,仅根据数据内在结构提取具有统计意义的区域信息。所以相对假设驱动方法具有更好的自适应性和信号发掘能力,具有比较强的客观性,客观揭示数据内在信息。缺陷在于提取的分量不一定与现实吻合,缺乏一定的现实意义。由于TPCA是线性组合,通过找到方差最大的方向进行投影,得到的变量仅仅是不相关,并不相互独立。PCA只用到了原始数据的二阶统计信息,而忽略了其高阶统计信息。因此,需要借助旋转主成分进行优化,找到更多主成分蕴含的物理意义。

参考文献 (20)

目录

    /

    返回文章
    返回