留言板

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

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

利用锁眼卫星影像提取历史地震同震位移——以1976年土耳其Chaldiran地震为例

卢乐浚 周宇

卢乐浚, 周宇. 利用锁眼卫星影像提取历史地震同震位移——以1976年土耳其Chaldiran地震为例[J]. 武汉大学学报 ● 信息科学版, 2021, 46(2): 289-295. doi: 10.13203/j.whugis20190071
引用本文: 卢乐浚, 周宇. 利用锁眼卫星影像提取历史地震同震位移——以1976年土耳其Chaldiran地震为例[J]. 武汉大学学报 ● 信息科学版, 2021, 46(2): 289-295. doi: 10.13203/j.whugis20190071
LU Lejun, ZHOU Yu. Extracting Surface Displacements of Historical Earthquakes Using KH-9 Satellite Images: A Case Example of 1976 Chaldiran Earthquake, Turkey[J]. Geomatics and Information Science of Wuhan University, 2021, 46(2): 289-295. doi: 10.13203/j.whugis20190071
Citation: LU Lejun, ZHOU Yu. Extracting Surface Displacements of Historical Earthquakes Using KH-9 Satellite Images: A Case Example of 1976 Chaldiran Earthquake, Turkey[J]. Geomatics and Information Science of Wuhan University, 2021, 46(2): 289-295. doi: 10.13203/j.whugis20190071

利用锁眼卫星影像提取历史地震同震位移——以1976年土耳其Chaldiran地震为例

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

第二次青藏高原综合科学考察研究项目 2019QZKK0901

国家自然科学基金 41874020

国家重点研发计划 2017YFC1500101

广东省引进创新科研团队——环南海地质过程与灾害创新团队 2016ZT06N331

详细信息

Extracting Surface Displacements of Historical Earthquakes Using KH-9 Satellite Images: A Case Example of 1976 Chaldiran Earthquake, Turkey

Funds: 

The Second Tibetan Scientific Expedition and Research Program (STEP) 2019QZKK0901

the National Natural Science Foundation of China 41874020

the National Key Research and Development Program of China 2017YFC1500101

Guangdong Province Introduced Innovative Research and Development Team of Geological Processes and Natural Disasters Around the South China Sea 2016ZT06N331

More Information
  • 摘要: 大地震的地表破裂及形变特征对研究地震成因机制、断层运动及大陆变形具有重要意义。随着卫星技术的日益成熟,自1992年美国加州Landers地震以来,光学和雷达影像被广泛应用于同震破裂、震后形变等定量研究。然而,由于缺乏更早之前的影像资料,90年代以前的历史地震研究无法深入。美国锁眼卫星KH-9影像的公开解决了震前卫星影像匮乏的问题,为研究70~90年代历史地震提供了可行的条件。首先借助伊朗Tabas-e-Golshan地震和Khuli-Boniabad地震的研究实例,概述了KH-9卫星影像测量历史同震位移的方法与进展;然后利用1976年(震前)和1980年(震后)的KH-9影像对1976年土耳其Chaldiran地震的同震位移进行了计算,测得该地震东西向同震位移量为(3.1±0.7)m,与实地测量的地表位移相符,表明该走滑地震没有明显的浅部滑动亏损现象;最后对KH-9卫星影像定量研究历史地震的未来应用和限制进行了讨论与总结。
  • 图  1  利用光学影像获取同震位移场流程图

    Figure  1.  Flowchart of Estimating Earthquake Displacements Using Opotical Imagery

    图  2  基于KH-9影像研究全球历史地震的分布图

    Figure  2.  Distribution Map of Global Historical Earthquakes that Have Been Studied Using KH-9 Imagery

    图  3  1978年伊朗Tabas-e-Golshan地震东西向位移场图[3]

    Figure  3.  E-W Displacement of the 1978 Tabas-e-Golshan Earthquake, Iran[3]

    图  4  1979年伊朗Khuli-Boniabad地震东西向位移场图[4]

    Figure  4.  E-W Displacement of the 1979 Khuli-Boniabad Earthquake, Iran[4]

    图  5  中东地区活动构造图[23]

    Figure  5.  Active Tectonic Map of the Middle East[23]

    图  6  1976年土耳其Chaldiran地震东西向位移图

    Figure  6.  E-W Displacement of the 1976 Chaldiran Earthquake, Turkey

    图  7  1978年伊朗Tabas-e-Golshan地震的南北向位移图[3]

    Figure  7.  S-N Displacement of the 1978 Tabas-e-Golshan Earthquake, Iran[3]

    表  1  Chaldiran地震研究所使用的影像数据

    Table  1.   Images Data Used for the 1976 Chaldiran Earthquake

    统计项 震前影像KH-9 震后影像
    KH-9 Landsat 8
    时间 1976-07-27 1980-07-27 2014-08-13
    相机类型 框幅式 框幅式 推扫式
    分辨率/m 6~9 6~9 15
    焦距/cm 30.5 30.5 243.8
    下载: 导出CSV
  • [1] Zhang P Z. Beware of Slowly Slipping Faults[J]. Nature Geoscience, 2013, 6(5): 323-324 doi:  10.1038/ngeo1811
    [2] Massonnet D, Rossi M, Carmona C, et al. The Displacement Field of the Landers Earthquake Mapped by Radar Interferometry[J]. Nature, 1993, 364(6 433): 138-142
    [3] Zhou Y, Walker R T, Hollingsworth J, et al. Coseismic and Postseismic Displacements from the 1978 Mw 7.3 Tabas-e-Golshan Earthquake in Eastern Iran[J]. Earth and Planetary Science Letters, 2016, 452: 185-196 doi:  10.1016/j.epsl.2016.07.038
    [4] Marchandon M, Vergnolle M, Cavalié O, et al. Earthquake Sequence in the NE Lut, Iran: Observations from Multiple Space Geodetic Techniques[J]. Geophysical Journal International, 2018, 215(3): 1 604-1 621 doi:  10.1093/gji/ggy364
    [5] 赵秋艳.美国成像侦查卫星的发展[J].光机电信息, 2001(10): 15-23 doi:  10.3969/j.issn.1007-1180.2001.10.003

    Zhao Qiuyan. Development of American Imaging Reconnaissance Satellites[J]. OME Information, 2001(10): 15-23 doi:  10.3969/j.issn.1007-1180.2001.10.003
    [6] 王鑫, 张景发, 姜文亮, 等.美国锁眼侦查卫星遥感数据在活动断层研究中的应用——以郯庐断裂带江苏段为例[J].遥感学报, 2018, 22(S1): 233-246

    Wang Xin, Zhang Jingfa, Jiang Wenliang, et al. Application of Keyhole Satellite Data in Active Fault Study: A Case Example of Jiangsu Segment of Tan-Lu Fault Zone[J]. Journal of Remote Sensing, 2018, 22(S1): 233-246
    [7] Pressel P. Spy in the Sky: The KH-9 Hexagon[J]. Optics and Photonics News, 2013, 24(10): 28-35 doi:  10.1364/OPN.24.10.000028
    [8] Surazakov A, Aizen V. Positional Accuracy Evaluation of Declassified Hexagon KH-9 Mapping Camera Imagery[J]. Photogrammetric Engineering & Remote Sensing, 2010, 76(5): 603-608 doi:  10.14358/PERS.76.5.603
    [9] Burnett M G. Hexagon (KH-9)–Mapping Camera Program and Evolution [M]. Chantilly: Center for the Study of National Reconnaissance(CSNR), 2012
    [10] Guillemette R. Declassified US Spy Satellites Reveal Rare Look at Secret Cold War Space Program[EB/OL]. http://www.space.com/12996-secret-spy-satellites-declassified-nro.html, 2011
    [11] Hanies G K. Critical to US Security: The Development of the GAMBIT and HEXAGON Satellite Reconnaissance System[M]. Chantilly: National Reconnaissance Office, 2012
    [12] Fowler M J. The Archaeological Potential of Declassified Hexagon KH-9 Panoramic Camera Satellite Photographs[J]. AARGnews, 2016, 53: 30-36 https://www.researchgate.net/publication/308684304_The_archaeological_potential_of_declassified_HEXAGON_KH-9_panoramic_camera_satellite_photographs
    [13] 杨士超.子像素偏移追踪算法及其在实际地震同震形变场监测中的应用[D].合肥: 中国科学技术大学, 2015

    Yang Shichao. Sub-pixel Offset Tracking Algorithm and Its Application to Detection the Co-seismic Deformation of Real Earthquakes[D]. Hefei: University of Science and Technology of China, 2015
    [14] 孙家柄.遥感原理与应用[M].武汉:武汉大学出版社, 2009

    Sun Jiabing. Principles and Applications of Remote Sensing[M]. Wuhan:Wuhan University Press, 2009
    [15] Leprince S, Barbot S, Ayoub F, et al. Automatic and Precise Orthorectification, Coregistration, and Subpixel Correlation of Satellite Images, Application to Ground Deformation Measurements[J]. IEEE Transactions on Geoscience and Remote Sensing, 2007, 45(6): 1 529-1 558 doi:  10.1109/TGRS.2006.888937
    [16] Copley A. Postseismic Afterslip 30 Years After the 1978 Tabas-e-Golshan (Iran) Earthquake: Observations and Implications for the Geological Evolution of Thrust Belts[J]. Geophysical Journal International, 2014, 197(2): 665-679 doi:  10.1093/gji/ggu023
    [17] 李怀英.外电外刊对1978年9月16日伊朗塔巴斯地震的有关报道及评论(综合汇编)[J].国际地震动态, 1979, 5: 4-6 http://www.cnki.com.cn/Article/CJFDTotal-GJZT197905001.htm

    Li Huaiying. The Reports and Comments About the September 16th, 1978 Tabas, Iran Earthquake from Foreign News[J]. Recent Developments in World Seismology, 1979, 5: 4-6 http://www.cnki.com.cn/Article/CJFDTotal-GJZT197905001.htm
    [18] Walker R, Jackson J, Baker C. Surface Expression of Thrust Faulting in Eastern Iran: Source Parameters and Surface Deformation of the 1978 Tabas and 1968 Ferdows Earthquake Sequences[J]. Geophysical Journal International, 2003, 152(3): 749-765 doi:  10.1046/j.1365-246X.2003.01886.x
    [19] Ibrion M, Mokhtari M, Nadim F. Earthquake Disaster Risk Reduction in Iran: Lessons and "Lessons Learned" from Three Large Earthquake Disasters—Tabas 1978, Rudbar 1990, and Bam 2003[J]. International Journal of Disaster Risk Science, 2015, 6(4): 415-427 doi:  10.1007/s13753-015-0074-1
    [20] Berberian M. Earthquake Faulting and Bedding Thrust Associated with the Tabas-e-Golshan (Iran) Earthquake of September 16, 1978[J]. Bulletin of the Seismological Society of America, 1979, 69(6): 1 861-1 887
    [21] Nowroozi A A, Mohajer‐Ashjai A. Fault Movements and Tectonics of Eastern Iran: Boundaries of the Lut Plate[J]. Geophysical Journal of the Royal Astronomical Society, 1985, 83(1): 215-237 doi:  10.1111/j.1365-246X.1985.tb05164.x
    [22] Gulkan P, Gurpinar A, Celebi M, et al. Engineering Report on the Muradiye-Çaldiran, Turkey, Earthquake of 24 November 1976[M]. Washington D C: National Academy of Sciences, 1978
    [23] Selçuk A S, Erturaç M K, Nomade S. Geology of the Çaldıran Fault, Eastern Turkey: Age, Slip Rate and Implications on the Characteristic Slip Behaviour[J]. Tectonophysics, 2016, 680: 155-173 doi:  10.1016/j.tecto.2016.05.019
    [24] Reilinger R, McClusky S, Vernant P, et al. GPS Constraints on Continental Deformation in the Africa‐Arabia‐Eurasia Continental Collision Zone and Implications for the Dynamics of Plate Interactions[J]. Journal of Geophysical Research: Solid Earth, 2006, 111(B5): B05411
    [25] Toksöz M N, Arpat E, Şaroglu F. East Anatolian Earthquake of 24 November 1976[J]. Nature, 1977, 270(5 636): 423-425 https://www.nature.com/articles/270423b0
    [26] Toksöz M N, Nábělek J, Arpat E. Source Properties of the 1976 Earthquake in East Turkey: A Comparison of Field Data and Teleseismic Results[J]. Tectonophysics, 1978, 49(3-4): 199-205 doi:  10.1016/0040-1951(78)90178-6
    [27] Wells D L, Coppersmith K J. New Empirical Relationships Among Magnitude, Rupture Length, Rupture Width, Rupture Area, and Surface Displacement[J]. Bulletin of the Seismological Society of America, 1994, 84(4): 974-1 002
  • [1] 李欣, 杨宇辉, 杨博, 尹峰.  利用方向相位特征进行多源遥感影像匹配 . 武汉大学学报 ● 信息科学版, 2020, 45(4): 488-494. doi: 10.13203/j.whugis20180445
    [2] 顾国华, 王武星.  GPS测得的2018年夏威夷6.9级地震与火山喷发地壳运动 . 武汉大学学报 ● 信息科学版, 2019, 44(8): 1191-1197, 1204. doi: 10.13203/j.whugis20180463
    [3] 袁修孝, 袁巍, 陈时雨.  基于图论的遥感影像误匹配点自动探测方法 . 武汉大学学报 ● 信息科学版, 2018, 43(12): 1854-1860. doi: 10.13203/j.whugis20180154
    [4] 詹总谦, 李一挥, 王陈东, 郑莉.  顾及局部相对几何变形改正的影像匹配和空三逐步精化方法 . 武汉大学学报 ● 信息科学版, 2018, 43(11): 1620-1627. doi: 10.13203/j.whugis20170342
    [5] 叶沅鑫, 慎利, 陈敏, 王继成.  局部相位特征描述的多源遥感影像自动匹配 . 武汉大学学报 ● 信息科学版, 2017, 42(9): 1278-1284. doi: 10.13203/j.whugis20150576
    [6] 王奋飞, 任鑫, 刘建军, 施俊, 李春来.  嫦娥二号影像密集匹配及DEM制作方法 . 武汉大学学报 ● 信息科学版, 2017, 42(7): 918-923. doi: 10.13203/j.whugis20150147
    [7] 顾国华, 王武星.  2016年新西兰7.8级大地震GPS观测结果与弹性回跳模型 . 武汉大学学报 ● 信息科学版, 2017, 42(11): 1673-1680. doi: 10.13203/j.whugis20170286
    [8] 徐亚明, 石娟, 安动动, 马旭东.  利用影像分割和匹配特征进行无人机影像变化检测 . 武汉大学学报 ● 信息科学版, 2016, 41(10): 1286-1291. doi: 10.13203/j.whugis20140873
    [9] 徐秋辉, 佘江峰, 宋晓群, 肖鹏峰.  利用改进的DCCD和SIFT描述符的影像匹配方法 . 武汉大学学报 ● 信息科学版, 2015, 40(12): 1613-1617,1645. doi: 10.13203/j.whugis20130753
    [10] 张春森, 吕佩育, 郭丙轩.  基于控制点约束影像的密集匹配及其在考古发掘中的应用 . 武汉大学学报 ● 信息科学版, 2015, 40(12): 1575-1581. doi: 10.13203/j.whugis20130510
    [11] 张永军, 熊金鑫, 余磊, 凌霄.  严密定位模型辅助的国产卫星影像匹配 . 武汉大学学报 ● 信息科学版, 2014, 39(8): 897-900. doi: 10.13203/j.whugis201200178
    [12] 闫利, 陈林.  一种改进的SURF及其在遥感影像匹配中的应用 . 武汉大学学报 ● 信息科学版, 2013, 38(7): 770-773.
    [13] 张谦, 贾永红, 胡忠文.  多源遥感影像配准中的SIFT特征匹配改进 . 武汉大学学报 ● 信息科学版, 2013, 38(4): 455-459.
    [14] 叶沅鑫, 单杰, 熊金鑫, 董来根.  一种结合SIFT和边缘信息的多源遥感影像匹配方法 . 武汉大学学报 ● 信息科学版, 2013, 38(10): 1148-1151.
    [15] 蔡华, 孙汉荣, 赵齐乐, 李鹏.  GPS测定的2011年日本9.0级地震的中国大陆地区同震位移场 . 武汉大学学报 ● 信息科学版, 2012, 37(8): 953-955.
    [16] 宋妍, 田玉刚, 贾小霞.  基于Zernike矩的多时相遥感影像匹配方法研究 . 武汉大学学报 ● 信息科学版, 2009, 34(10): 1226-1230.
    [17] 江万寿, 郑顺义, 张祖勋, 张剑清.  航空影像特征匹配研究 . 武汉大学学报 ● 信息科学版, 2003, 28(5): 510-513.
    [18] 张祖勋, 张剑清, 张力.  一种利用影像匹配的城区三维景观图生成途径 . 武汉大学学报 ● 信息科学版, 1999, 24(4): 283-289.
    [19] 曾湧, 廖明生, 张剑清, 沈未名.  保持匹配一致性的遥感影像压缩编码 . 武汉大学学报 ● 信息科学版, 1997, 22(3): 245-247,251.
    [20] 李德仁, 王树根.  数字影像匹配质量的一种自动诊断方法 . 武汉大学学报 ● 信息科学版, 1995, 20(1): 1-6.
  • 加载中
图(7) / 表(1)
计量
  • 文章访问数:  29
  • HTML全文浏览量:  3
  • PDF下载量:  14
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-09-04
  • 刊出日期:  2021-02-05

利用锁眼卫星影像提取历史地震同震位移——以1976年土耳其Chaldiran地震为例

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

    第二次青藏高原综合科学考察研究项目 2019QZKK0901

    国家自然科学基金 41874020

    国家重点研发计划 2017YFC1500101

    广东省引进创新科研团队——环南海地质过程与灾害创新团队 2016ZT06N331

    作者简介:

    卢乐浚,博士生,研究方向为构造大地测量。lulj8@mail2.sysu.edu.cn

    通讯作者: 周宇,博士,教授。zhouyu36@mail.sysu.edu.cn
  • 中图分类号: P237;P315

摘要: 大地震的地表破裂及形变特征对研究地震成因机制、断层运动及大陆变形具有重要意义。随着卫星技术的日益成熟,自1992年美国加州Landers地震以来,光学和雷达影像被广泛应用于同震破裂、震后形变等定量研究。然而,由于缺乏更早之前的影像资料,90年代以前的历史地震研究无法深入。美国锁眼卫星KH-9影像的公开解决了震前卫星影像匮乏的问题,为研究70~90年代历史地震提供了可行的条件。首先借助伊朗Tabas-e-Golshan地震和Khuli-Boniabad地震的研究实例,概述了KH-9卫星影像测量历史同震位移的方法与进展;然后利用1976年(震前)和1980年(震后)的KH-9影像对1976年土耳其Chaldiran地震的同震位移进行了计算,测得该地震东西向同震位移量为(3.1±0.7)m,与实地测量的地表位移相符,表明该走滑地震没有明显的浅部滑动亏损现象;最后对KH-9卫星影像定量研究历史地震的未来应用和限制进行了讨论与总结。

English Abstract

卢乐浚, 周宇. 利用锁眼卫星影像提取历史地震同震位移——以1976年土耳其Chaldiran地震为例[J]. 武汉大学学报 ● 信息科学版, 2021, 46(2): 289-295. doi: 10.13203/j.whugis20190071
引用本文: 卢乐浚, 周宇. 利用锁眼卫星影像提取历史地震同震位移——以1976年土耳其Chaldiran地震为例[J]. 武汉大学学报 ● 信息科学版, 2021, 46(2): 289-295. doi: 10.13203/j.whugis20190071
LU Lejun, ZHOU Yu. Extracting Surface Displacements of Historical Earthquakes Using KH-9 Satellite Images: A Case Example of 1976 Chaldiran Earthquake, Turkey[J]. Geomatics and Information Science of Wuhan University, 2021, 46(2): 289-295. doi: 10.13203/j.whugis20190071
Citation: LU Lejun, ZHOU Yu. Extracting Surface Displacements of Historical Earthquakes Using KH-9 Satellite Images: A Case Example of 1976 Chaldiran Earthquake, Turkey[J]. Geomatics and Information Science of Wuhan University, 2021, 46(2): 289-295. doi: 10.13203/j.whugis20190071
  • 地震是一种对人类安全威胁很大的自然灾害,但人们对大地震(Mw≥7)的地震周期研究进展缓慢,其原因之一是大地震是一种小概率事件,可供研究的大地震震例很少。尤其是在大陆内部,断层应变积累缓慢,大地震周期往往很漫长,难以对其进行深入的研究[1]。目前,利用大地测量与遥感技术对已发生的大地震的定量研究仅能追溯到1992年美国加州Landers地震[2],亟需一种技术能够将大地测量与遥感技术用于对1992年以前的地震进行研究,但由于缺乏震前的影像数据,历史地震的研究也受到限制。美国锁眼(KeyHole,KH)卫星影像的公开解决了此前震前卫星影像匮乏的难题,为研究70~90年代的历史地震提供了可能。近年来,一些学者基于KH-9影像成功提取了全球部分70~90年代历史地震的同震位移,打开了定量研究历史地震的新局面[3-4]。针对历史地震的研究,本文首先介绍了KH-9影像和数据处理方法,对目前已取得的研究成果进行了概述;然后利用KH-9和Landsat 8影像提取了1976年土耳其Mw 7.3级地震的同震位移;最后讨论了该方法的局限性以及未来的可能应用。

    • KH系列卫星最开始发射于20世纪60年代,总发射次数达270余次。美国将此前的“日冕(CORONA)”“萨莫斯(SAMOS)”“发现者(DISCOVER)”等早期军事侦查卫星计划相合并,发展建立起KH卫星计划。KH系列卫星型号从KH-1发展到KH-12,根据装载的仪器类型可分为3代侦察卫星。第1代照相侦察卫星KH-1~6,类型为胶片返回式照像侦察卫星,装载有全景相机,工作寿命短,分辨率低。第2代照相侦察卫星KH-7~10,在第1代的基础上额外装载了红外扫描仪和多光谱扫描仪,工作寿命及分辨率均有提高。第3代照相侦察卫星KH-11和KH-12,类型为数字成像无线电传输卫星,具备数据实时传输功能,工作寿命及分辨率进一步提升[5-6]

      本文使用的KH-9卫星是由美国中央情报局于1965年计划研发的,用于替代将要退役的KH-4卫星。KH-9卫星于1971—1986年间共发射了20次,轨道高度为171 km,其中最后一次发射失败。该卫星是一颗可回收胶卷舱的卫星,卫星将相机拍摄到的胶卷影像存储到回收舱中,最终由回收舱将卫星影像运送回地球[7]

      在1973—1980年服役期间,KH-9卫星共完成12次制图任务(任务编号1 205~1 216)[8]。KH-9卫星使用框幅式相机进行拍摄,总共收集到29 000幅地面分辨率为6~9 m的黑白图像用于制图,每一幅图像覆盖的地面面积为125 km×250 km,拍摄区域基本覆盖了整个地球[9]。2011-09-17,美国国家侦察局将受严密保护的KH-9卫星影像资料公开,并将KH-9卫星在美国航空航天博物馆停车场向公众展示一天,以此纪念美国国家侦察局建局50周年[10]。在KH-9卫星的公开资料中,用于制图的相机参数并没有详细给出。经美国地质调查局(United States Geological Survey,USGS)扫描影像后对比实验发现,该摄影机与美国航空航天局的大像幅相机相似,使用的胶卷尺寸为23 cm×46 cm,镜头焦距为30.5 cm[11]

    • 基于KH-9影像估算地震同震位移的步骤包括影像获取、影像配准、正射校正、影像匹配,其流程如图 1所示。

      图  1  利用光学影像获取同震位移场流程图

      Figure 1.  Flowchart of Estimating Earthquake Displacements Using Opotical Imagery

      1)影像获取。震前的KH-9影像可以通过USGS EarthExplorer网站(https://earthexplorer.usgs.gov/)以30美元每幅的价格购买[12]。震后影像可根据需求使用免费的Landsat 8影像或购买SPOT、WorldView等高分辨率影像。

      2)影像配准。由于不同影像的成像条件不同(如成像时间、成像视角、传感器等),影像间的坐标系统存在差异,尤其是KH-9影像未经过几何校正和地理编码,因此在进行影像匹配前,需要利用控制点对影像进行匹配。由于靠近震中位置地表形变大,控制点选取应远离震中。通过在震前和震后影像中选取足够的同名点,计算影像间的映射函数[13],进而对KH-9影像进行校正,完成对震前KH-9影像与震后影像的配准。

      3)正射校正。卫星影像在成像过程中由于地面高程、大气折射以及地球曲率等因素影响,会造成影像发生不同程度的几何畸变失真。正射校正的作用是消除成像过程中由于地形起伏、高程变化引起的畸变,使得输出的正射校正影像为平面真实影像[14]

      4)影像匹配。由于大部分地震造成的形变量小于KH-9影像的分辨率,利用KH-9影像测量同震位移时需要达到亚像素级。为了保证精度,影像匹配通常在频率域进行[15]

      首先提取同一区域震前i1及震后i2两幅正射校正后的卫星影像。以卫星的质心为原点(O)建立空间坐标系(XYZ),Y轴为卫星运动方向,Z轴由地球质心指向卫星方向,X轴为YOZ平面的垂线方向,计算两幅影像的空间位移量:

      $$ {i_2}\left( {x, y} \right) = {i_1}\left( {x - {\varDelta _x}, y - {\varDelta _y}} \right) $$ (1)

      式中,$\left( {x, y} \right)$表示影像像素坐标;${{\varDelta _x}}$表示X轴上的位移量;${{\varDelta _y}}$表示Y轴上的位移量。代入以下等式关系:

      $$ {I_2}\left( {{\omega _x}, {\omega _y}} \right) = {I_1}\left( {{\omega _x}, {\omega _y}} \right){{\rm{e}}^{ - {\rm{j}}\left( {{\omega _x}{\varDelta _x} + {\omega _y}{\varDelta _y}} \right)}} $$ (2)

      式中,${I_1}$和${I_2}$分别表示震前和震后影像的傅里叶变换函数;${{\omega _x}}$和${{\omega _y}}$分别表示列与行的频率变量。将正射校正后的震前和震后影像进行傅里叶变换,从空间域转换到频率域,并代入下式计算:

      $$ {C_{{i_1}{i_2}}}\left( {{\omega _x}, {\omega _y}} \right) = \frac{{{I_1}\left( {{\omega _x}, {\omega _y}} \right)I_2^*\left( {{\omega _x}, {\omega _y}} \right)}}{{\left| {{I_1}\left( {{\omega _x}, {\omega _y}} \right)I_2^*\left( {{\omega _x}, {\omega _y}} \right)} \right|}} = {{\rm{e}}^{{\rm{j}}\left( {{\omega _x}{\varDelta _x} + {\omega _y}{\varDelta _y}} \right)}}{\rm{}} $$ (3)

      式中,${C_{{i_1}{i_2}}}\left( {{\omega _x}, {\omega _y}} \right)$表示交叉谱;*表示复共轭运算。通过计算交叉谱的相位得到像素的相对位移,运用傅里叶反函数${F^{ - 1}}$得到相关匹配方程:

      $$ \delta \left( {x + {\varDelta _x}, y + {\varDelta _y}} \right) = {F^{ - 1}}\left\{ {{{\rm{e}}^{{\rm{j}}\left( {{\omega _x}{\varDelta _x} + {\omega _y}{\varDelta _y}} \right)}}} \right\} $$ (4)

      式中,$\delta $表示待匹配窗口与匹配窗口的相关系数,用于衡量影像匹配的可靠性。若得到互相关匹配系数$\delta $大于设定的阈值,相应的${{\varDelta _x}}$和${{\varDelta _y}}$即为估算的地表形变。

    • 目前,有学者利用KH-9影像对全球两例70年代地震进行了定量研究,即1978年伊朗Tabas-e-Golshan地震[3]和1979年伊朗Khuli-Boniabad地震[4],如图 2所示。

      图  2  基于KH-9影像研究全球历史地震的分布图

      Figure 2.  Distribution Map of Global Historical Earthquakes that Have Been Studied Using KH-9 Imagery

    • 1978-09-16,伊朗东部Tabas地区发生Mw 7.3级地震,震中位置为(54.7°E,32.4°N),是由盲逆断层导致的[16-19]。震后野外调查发现85 km长的破裂带,测量到的单一破裂带上的最大同震滑移量仅为35 cm[20]。此次Mw 7.3级Tabas-e-Golshan地震的同震滑移量到底有多大,是一个谜题。

      为了估算此次地震的实际滑移量,Zhou等[3]获取了覆盖Tabas区域的1974年KH-9影像(震前)和2013年SPOT-6影像(震后)。基于§2描述的方法,Zhou等[3]获取了地震的形变场(图 3(a)),计算得到1974—2013年间Tabas逆断层东西向缩短约3.4 m(图 3(b)),反演结果表明,这一缩短量反映了断层面上约7 m的滑动。Tabas地区同时发育了两级河流阶地,通过河流阶地的高程剖面反演,Zhou等[3]发现Tabas断层显示出特征地震的行为,每次地震滑移量均为7 m。考虑到Tabas断层的地质滑动速率为1 mm/a,该地区地震复发周期约为3 500 a。

      图  3  1978年伊朗Tabas-e-Golshan地震东西向位移场图[3]

      Figure 3.  E-W Displacement of the 1978 Tabas-e-Golshan Earthquake, Iran[3]

    • 1979-11-27,伊朗东北部南霍拉桑省Koli村庄发生Mw 7.1级地震[4]。Khuli-Boniabad地震形成了震区中65 km长的东西向地表破裂,据野外调查,观测到震区内最大水平位移量为2.55 m,最大垂直位移量为3.80 m[21]

      Marchandon等[4]利用1976年KH-9影像(震前)和2017年Sentinel-2影像(震后)的匹配,获取了Khuli-Boniabad地震的形变场。结果显示,该地震在地表处形成了一条东西向的破裂痕迹,图 4A-B段、B-C段、C-D段的断层走向(断层迹线与北方向的夹角)分别为65°N、85°N、81°N,D-F段破裂不明显,而F处西端弯曲至160°N。P1东西向位移剖面显示,断层位移量为(4.92±0.73)m。

      图  4  1979年伊朗Khuli-Boniabad地震东西向位移场图[4]

      Figure 4.  E-W Displacement of the 1979 Khuli-Boniabad Earthquake, Iran[4]

    • 1976-11-24 14:00,土耳其东部城市范市(Van)发生Mw 7.3级地震,其中损害最严重的地区为Chaldiran城镇。据USGS的震源数据分析显示,该次地震震中位置为(44.02°E,39.10°N)。美国国家科学院的地震报告分析指出,该地震是由NW-SE右旋走滑断层运动所引起的,而该断层在地震前并未被发现[22]

    • 土耳其东部与伊朗西北部相接壤,位于东安纳托利亚高原,处于欧亚板块和阿拉伯板块持续碰撞的复杂交接地带,地质构造复杂。该交接地带分布着多个子块体,其中伊朗板块、高加索板块夹持于两大板块之间,且与相接板块交互作用,构成Zagros褶皱俯冲带。该地带的形成主要与中新世中期新特地斯洋的向北俯冲及消亡有关[23]

      Chaldiran走滑断层位于高加索板块的南缘,与伊朗板块的交接位置也处于Zagros褶皱俯冲带上,在地震发生前未被发现。基于GPS数据和弹性板块模型,Reilinger等[24]计算得到了伊朗板块、高加索板块、阿拉伯板块分别以相对于欧亚板块约16 mm/a、13 mm/a、15 mm/a的速度向北运动(图 5),速度差异导致两个板块交接处发生NW-SE向伸展,形成右旋走滑断层,从而引发本次地震。

      图  5  中东地区活动构造图[23]

      Figure 5.  Active Tectonic Map of the Middle East[23]

      据统计,Chaldiran地震共造成约3 800人死亡,9 000余间建筑严重坍塌,10 000余间建筑不同程度的破坏,破坏性极强[22]。该走滑断层的破裂在地表延伸长为55 km,顺破裂轨迹中西侧走向为N70°W,在Chaldiran东侧破裂有一定弯曲,走向变为N45°W。破裂轨迹西端延伸至西侧山脉处突然消失,东端地表位移变化逐渐减小,至伊朗边界处消失,可能在底下有继续延伸趋势。通过测量道路、河床、车辆轨迹等得到横向位移约200~250 cm,部分地区横向位移达370 cm[22, 25-26]

    • 为获取Chaldiran地震的同震位移,从USGS购买了1976-07-27(震前)和1980-07-27(震后)的KH-9影像,同时下载了1幅2014年的Landsat 8影像用于影像匹配(表 1)。首先,分别在两幅KH-9和Landsat 8影像上选取15个同名点,利用同名点在Landsat 8影像上已知的空间坐标计算两幅KH-9影像的外方位元素;然后利用30 m分辨率的ALOS World 3-D DEM对两幅KH-9影像进行正射校正,恢复影像真实的地理坐标;最后使用COSI-Corr软件对两幅正射校正后的KH-9影像进行子像素匹配,匹配窗口定为32×32像素,窗口移动步长定为8×8像素,得到位移场,如图 6所示。

      表 1  Chaldiran地震研究所使用的影像数据

      Table 1.  Images Data Used for the 1976 Chaldiran Earthquake

      统计项 震前影像KH-9 震后影像
      KH-9 Landsat 8
      时间 1976-07-27 1980-07-27 2014-08-13
      相机类型 框幅式 框幅式 推扫式
      分辨率/m 6~9 6~9 15
      焦距/cm 30.5 30.5 243.8

      图  6  1976年土耳其Chaldiran地震东西向位移图

      Figure 6.  E-W Displacement of the 1976 Chaldiran Earthquake, Turkey

      在靠近震中位置取剖面A-A',跨断层处位移剖面出现了明显的错动;用直线拟合断层两侧的位移,得到断层上东西向的位错约为(3.1±0.7)m。断层走向近东西,因此走滑量基本上表现为东西向位移。利用KH-9影像测得的断层位错(3.1±0.7)m反映的是地表以下的断层滑动,与Toksöz等[25]野外测量的地表走滑位移量(最大370 cm)相符,说明该地震不存在同震浅层滑动亏损的现象。

    • 利用KH-9影像为研究70~90年代历史地震带来了新的机遇。本文对KH-9影像处理方法和目前已取得的研究成果(包括1978年Tabas-e-Golshan地震和1979年Khuli-Boniabad地震)进行了梳理,同时将KH-9影像应用到1976年土耳其Chaldiran地震的研究中,获取了同震位移(约为(3.1±0.7)m),与野外调查的地表位错吻合,验证了该方法的可靠性。

      然而,KH-9影像用于历史地震的研究也具有一定的局限性,例如:①部分KH-9影像在扫描过程中产生了非常严重的扫描畸变,这一畸变主要在南北向,表现为规则的跳变[3],如图 7所示。因此目前对历史地震的研究都以东西方向为主,南北方向的形变分量基本可忽略。②KH-9影像分辨率不高,约为6~9 m(取决于数字化过程中的扫描分辨率),现有影像匹配算法的精度约为1/10个像素,所以基于KH-9影像的匹配精度约为0.9 m,KH-9影像仅能用于同震位移超过此量级的地震。根据Wells等[27]的地震定标律:

      图  7  1978年伊朗Tabas-e-Golshan地震的南北向位移图[3]

      Figure 7.  S-N Displacement of the 1978 Tabas-e-Golshan Earthquake, Iran[3]

      $$ M = 6.93 + 0.82\lg D $$ (5)

      可以计算出KH-9影像适用的最小地震震级为Mw 6.9。式(5)中,M表示地震矩震级;D表示平均地表位移,单位为m。考虑到扫描畸变的影响,KH-9影像可能仅适用于Mw 7级以上的大地震。

      对于全球70~90年代Mw 7级以上的大地震(例如1976年唐山地震),可以利用KH-9影像定量估算同震位移,同时结合干涉合成孔径雷达数据观测可能的震后效应或震间应变积累过程,从而完善对地震周期形变的认识和理解。

参考文献 (27)

目录

    /

    返回文章
    返回