• ISSN 1673-5722
  • CN 11-5429/P

DEM分辨率对提取东北地区坡度的影响研究

刘佳轩 高瑞 师黎静

文鑫涛,李华玥,段乙好,郑通彦,2021. 2020年中国大陆地震灾害损失述评. 震灾防御技术,16(4):651−656. doi:10.11899/zzfy20210406. doi: 10.11899/zzfy20210406
引用本文: 刘佳轩,高瑞,师黎静,2021. DEM分辨率对提取东北地区坡度的影响研究. 震灾防御技术,16(2):237−244. doi:10.11899/zzfy20210202. doi: 10.11899/zzfy20210202
Wen Xintao, Li Huayue, Duan Yihao, Zheng Tongyan. Earthquake Disasters Loss on Chinese Mainland in 2020[J]. Technology for Earthquake Disaster Prevention, 2021, 16(4): 651-656. doi: 10.11899/zzfy20210406
Citation: Liu Jiaxuan, Gao Rui, Shi Lijing. Research on the Influence of DEM Resolution on the Slope Extraction in Northeast, China[J]. Technology for Earthquake Disaster Prevention, 2021, 16(2): 237-244. doi: 10.11899/zzfy20210202

DEM分辨率对提取东北地区坡度的影响研究

doi: 10.11899/zzfy20210202
基金项目: 国家重点研发计划(2018YFC1504601);国家自然科学基金面上项目(51978635)
详细信息
    作者简介:

    刘佳轩,男,生于1995年。硕士研究生。主要从事工程抗震方面的研究。E-mail:liujiaxuanhit@163.com

    通讯作者:

    高瑞,男,生于1993年。硕士研究生。主要从事工程抗震方面的研究。E-mail:gaoruihy@163.com

Research on the Influence of DEM Resolution on the Slope Extraction in Northeast, China

  • 摘要: 目前国内外广泛采用的坡度-等效剪切波速关系模型中,坡度主要基于30″分辨率数字高程数据的提取。DEM分辨率对坡度的提取具有重要影响,且对于不同地区、地貌单元的影响不同。本文以我国东北地区1″、3″、30″分辨率DEM为基础,分别提取研究区域内公里网格点的百分比坡度,分析不同分辨率下提取坡度差异的原因。区分平原、丘陵和山地统计分析坡度随分辨率变化特征,通过比较多种拟合模型的拟合优度和坡度分级平均相对差值,给出不同分辨率下的最佳坡度转换关系。研究结果表明,坡度体现计算面积的平均特性,随着分辨率的增大,栅格计算面积减小,导致坡度有所差异,区域坡度变化越大差异越大;在我国东北地区,随着DEM分辨率的降低,坡度-频率分布曲线向低坡度段移动,且均值、方差及值域变小,不同地貌的坡度-频率分布曲线向低坡度方向移动,且峰值区域变窄,其中平原地貌变化趋势性显著,丘陵次之,山地最小;平原及山地地貌最佳坡度拟合为线性拟合,丘陵地貌下1″、3″与30″分辨率下最佳坡度拟合分别为线性拟合和多项式拟合。
  • 2020年,中国境内共发生5.0级及以上地震28次(中国大陆地区发生20次,台湾地区及海域地区发生8次),其中5.0~5.9级地震25次,6.0~6.9级地震3次,最大地震为2020年7月23日发生在西藏那曲市尼玛县的6.6级地震(表1图1)。2020年中国大陆发生的20次5.0级以上地震中有10次发生在新疆、6次发生在西藏,6.0级以上地震中有2次发生在新疆、1次发生在西藏。

    表 1  2020年中国5.0级及以上地震事件目录及成灾事件
    Table 1.  Catalogue of earthquakes (MS≥5.0) and disaster - causing events in China, 2020
    序号时间纬度/°经度/°震中位置震级/M成灾事件
    12020-01-16 16:32:3841.2183.6新疆阿克苏地区库车市5.6(1)
    22020-01-18 00:05:5039.8377.18新疆喀什地区伽师县5.4
    32020-01-19 21:27:5539.8377.21新疆喀什地区伽师县6.4(2)
    42020-01-19 22:23:0139.8977.46新疆克孜勒苏州阿图什市5.2
    52020-01-25 06:56:0531.9895.09西藏昌都市丁青县5.1
    62020-01-29 07:39:2927.16126.6东海海域5.3
    72020-02-03 00:05:4130.74104.46四川成都市青白江区5.1
    82020-02-15 19:00:0723.95121.49台湾花莲县5.4
    92020-02-21 02:01:4034.5685.68西藏阿里地区改则县5.0
    102020-02-21 23:39:1439.8777.47新疆喀什地区伽师县5.1
    112020-03-10 02:12:1132.8485.52西藏阿里地区改则县5.0
    122020-03-12 23:44:0332.8885.55西藏阿里地区改则县5.1
    132020-03-20 09:33:1528.6387.42西藏日喀则市定日县5.9
    142020-03-23 03:21:3941.7581.11新疆阿克苏地区拜城县5.0
    152020-04-01 20:23:2733.0498.92四川甘孜州石渠县5.6(3)
    162020-05-03 11:24:4023.29121.6台湾台东县海域5.4
    172020-05-06 18:51:0039.7174.1新疆克孜勒苏州乌恰县5.0
    182020-05-09 23:35:5940.7778.76新疆阿克苏地区柯坪县5.2
    192020-05-18 21:47:5927.18103.16云南昭通市巧家县5.0(4)
    202020-06-14 04:18:5924.29122.41台湾宜兰县海域5.5
    212020-06-26 05:05:2035.7382.33新疆和田地区于田县6.4(5)
    222020-07-12 06:38:2539.78118.44河北唐山市古冶区5.1
    232020-07-13 09:28:0244.4280.82新疆伊犁州霍城县5.0
    242020-07-23 04:07:2033.1986.81西藏那曲市尼玛县6.6
    252020-07-26 20:52:2724.27122.48台湾花莲县海域5.5
    262020-09-29 04:50:5322.29121.1台湾台东县海域5.0
    272020-09-30 12:37:1824.85122.14台湾宜兰县海域5.0
    282020-12-10 21:19:5824.74121.99台湾宜兰县海域5.8
    注:“()”中表示为地震灾害事件。
    下载: 导出CSV 
    | 显示表格
    图 1  2020年中国5.0级及以上地震事件及成灾事件空间分布图
    Figure 1.  Spatial distribution of earthquakes (MS≥5.0) in China, 2020

    2020年,中国大陆地区共发生地震灾害事件5次,造成5人死亡,30人受伤,直接经济损失约18.47亿元(表2)。其中,灾害损失最严重的地震为1月19日新疆伽师6.4级地震,造成1人死亡,2人受伤,直接经济损失15.26亿元。人员伤亡最严重的地震为5月18日云南巧家5.0级地震,造成4人死亡,28人受伤,直接经济损失约1.04亿元。

    表 2  2020年中国大陆地震灾害损失
    Table 2.  Losses caused by earthquake disasters in China mainland, 2020
    序号时间震中位置震级/M死亡人数/人受伤人数/人直接经济损失/万元
    12020-01-16 16:32新疆阿克苏地区库车市5.600712
    22020-01-19 21:27新疆喀什地区伽师县6.412152 642
    32020-04-01 20:23四川甘孜州石渠县5.60019 242.69
    42020-05-18 21:47云南昭通市巧家县5.042810 430
    52020-06-26 05:05新疆和田地区于田县6.4001 650
    下载: 导出CSV 
    | 显示表格

    (1)地震成灾事件总体偏少

    2020年共发生5次地震灾害事件,低于2000年以来平均水平。

    (2)地震灾害相对集中

    5次地震灾害事件中有3次发生在新疆维吾尔自治区,共造成1人死亡,2人受伤,直接经济损失约15.5亿元。地震造成的人员伤亡主要集中在云南省,其中云南省巧家县5.0级地震虽未造成巨大的直接经济损失,但造成4人死亡,28人受伤。

    发生在新疆维吾尔自治区的3次地震灾害事件分别为:①2020年1月16日16时32分阿克苏地区库车市发生的5.6级地震,震源深度16 km,未造成人员伤亡,直接经济损失712万元。震中位于库车市,极震区烈度为Ⅵ度,Ⅵ度区面积813.6 km2表3)。灾区场地位于塔里木盆地腹地,场地对地震动有显著的放大作用,易产生不均匀沉降,加重了建筑物震害。灾区农居房屋抗震能力普遍提高,震中“安居富民房”与农村安居房均未出现破坏,有效保护灾区群众生命财产安全,同时缓解了抗震救灾和转移安置压力,仅少数建造年代较早的抗震安居房出现破坏。②2020年1月19日21时27分喀什地区伽师县发生的6.4级地震,震源深度16 km,造成1人死亡,2人受伤,直接经济损失152 642万元。Ⅵ度及以上区总面积为7 599 km2,其中,Ⅵ度区面积4 945 km2;Ⅶ度区面积为2 397 km2;Ⅷ度区面积为257 km2,主要涉及伽师县西克尔库勒镇、古勒鲁克乡。极震区烈度为Ⅷ度,农村安居房设防烈度为8度,安居房主要承重构件未发生明显破坏,有效抵御了本次地震灾害,在保障人民群众生命财产安全及震后转移安置、余震防范中发挥了重要作用。本次地震为前震-主震-余震型,前震为5.4级,主震为6.4级,最大余震为5.2级,对震中区产生多次震害影响,自建砖木结构房屋墙体为黏土砖砌筑,震害叠加效应显著。当地群众防震减灾意识在历次地震灾害实践中不断得到加强,防震减灾意识、应急避险知识和自救互救能力等得到明显提升。③2020年6月26日5时5分,和田地区于田县发生的6.4级地震,震源深度10 km,未造成人员伤亡,直接经济损失1 650万元。于田县阿羌乡3.5 km道路出现塌方,1座中桥出现多处裂缝,1座自建木桥受损,1个涵洞严重受损,1座拱桥受损。本次地震有感范围较广,但破坏性较小。震区附近乡镇村庄居民普遍反映地震造成的晃动强度一般,大部分区域房屋抗震设防水平高,“富民安居工程”房屋和早期的抗震安居房屋基本完好,仅个别民居院内的自建砖木结构房屋出现轻微程度破坏。

    表 3  2020年中国大陆地震灾区范围统计
    Table 3.  Statistics of range of earthquake disaster area in China mainland, 2020
    序号时间震中位置震级/M极震区烈度/度震源深度/km乡镇
    /个
    人口
    /人
    各烈度区面积/km2
    12020-01-16 16:32新疆阿克苏地区库车市5.61612685813.600
    22020-01-19 21:27新疆喀什地区伽师县
    6.41612279 3884 9452 397257
    32020-04-01 20:23四川甘孜州石渠县5.61063603 1400
    42020-05-18 21:47云南昭通市巧家县5.08667 08333000
    52020-06-26 05:05新疆和田地区于田县6.410
    注:新疆伽师地震Ⅵ度区之外的部分地区也受到波及,个别老旧房屋出现破坏受损现象,为Ⅵ度异常点;四川甘孜州石渠县地震位于Ⅵ度区之外色达县泥朵镇和其他地区也受到波及,零星房屋有破坏现象。
    下载: 导出CSV 
    | 显示表格

    (3)地震次生灾害突出

    今年发生的地震灾害中,次生灾害多发,如新疆伽师地震造成震区1座水库出现险情,当地政府紧急疏散安置受影响群众,云南巧家地震次生地质灾害造成2人死亡,数人受伤,震区交通等基础设施受损,再次为各级政府敲响警钟,应对西部地区地震次生灾害引起足够重视,才能进一步减轻地震灾害。

    1991—2020年主要震害统计数据如表4所示,3个主要统计数据分布情况(许永江等,1999郑通彦等,20122015a2015b陈通等,2016文鑫涛等,2018林向洋等,20182020a2020b)如图2所示。由表4图2可知,2011—2020年共造成2 063.08亿元的经济损失,其中2013、2014年地震造成的灾害是2011年以来最严重的,这两年造成灾害最重的地震分别为2013年4月20日四川芦山7.0级地震和2014年8月3日云南鲁甸6.5级地震,这2年合计的经济损失、人员死亡(失踪)和人员受伤数量分别占2011—2020年总数的82.9%、81.8%和65.5%。

    表 4  1991—2020年主要震害统计数据
    Table 4.  The statistics of the annual damage caused by earthquakes from 1991 to 2020
    年份成灾次数/次死亡人数/人受伤人数/人直接经济损失/亿元
    2011153250660.11
    201211861 33182.88
    20131429415 671995.36
    2014107363 688355.64
    201512331 217180.00
    201616210366.80
    20171137638217.40
    20181108127.30
    2019131741159.12
    2020553018.47
    2011—2020118124223 6762 063.08
    2001—201010872372398 9178 984.90
    1991—200013056452 668120.18
    下载: 导出CSV 
    | 显示表格
    图 2  1991—2020年中国大陆成灾地震年度震害分布图
    Figure 2.  Distribution of annual earthquakes losses in China mainland from 1991 to 2020

    近10年中国大陆年均成灾地震频次较1991年以来的总体水平略低,近10年期间地震灾害造成的经济损失、人员死亡(失踪)和人员受伤数量分别占30年来的18.47%、1.67%和4.98%,可见地震灾害造成的经济损失、人员伤亡数量均低于1991年以来的平均水平。

    2020年我国大陆未发生特重大地震灾害事件,地震灾害损失总体偏轻,成灾地震次数较少,低于2000年以来的平均水平。地震灾害相对集中,5次成灾地震中3次发生在新疆维吾尔自治区,四川和云南各发生1次。地震次生灾害多发,小震致灾致亡现象突出,云南省巧家县5.0级地震造成4人死亡,28人受伤,其中次生地质灾害造成2人死亡。各级政府应对西部地区地震次生灾害引起足够重视,进而减轻人员伤亡。

  • 图  1  东北地区不同分辨率下高程数据模型图

    Figure  1.  Elevation maps of 1 arc second, 3 arc second, and 30 arc second resolution in the Northeast, China

    图  2  计算高程梯度变化率的像元窗口

    Figure  2.  3×3 unit schematic diagram of slope calculation point

    图  3  东北地区不同分辨率下坡度图

    Figure  3.  Slope maps of 1 arc second, 3 arc second, and 30 arc second resolution in the Northeast, China

    图  4  不同DEM分辨率下的坡度统计

    Figure  4.  Slope statistics of different DEM resolutions

    图  5  地貌分布图

    Figure  5.  Landform distribution map

    图  6  不同分辨率下不同地貌单元坡度-频率分布曲线

    Figure  6.  Slope-frequency distribution curves of different topography at different resolutions

    图  7  不同地貌单元不同分辨率下坡度转换结果

    Figure  7.  Conversion results of DEM slopes with different resolutions under different landforms

    表  1  不同地貌单元回归模型拟合参数及评价参数(1″与30″分辨率下)

    Table  1.   Fitting and evaluation parameters of regression models for different landforms (1 arc second and 30 arc second)

    地貌单元类型拟合模型拟合参数R2ΔmeanK
    abcd
    平原线性函数00.085 630.313 400.535 70.585 0
    多项式函数0.114 600−0.154 800.090 5000.322 200.763 90.421 7
    幂函数0.047 1600.658 5000.080 000.892 30.089 6
    丘陵线性函数00.132 570.552 500.320 81.722 2
    多项式函数0.259 270−0.736 800.711 4600.610 810.486 51.255 5
    幂函数0.043 2490.329 2300.077 80.842 30.092 3
    山地线性函数00.199 240.598 490.345 11.734 2
    多项式函数0.445 110−0.996 520.686 5900.687 300.714 80.961 5
    幂函数0.079 9900.287 7300.092 450.965 30.095 7
    下载: 导出CSV

    表  2  不同地貌单元回归模型拟合参数及评价参数(3″与30″分辨率下)

    Table  2.   Fitting and evaluation parameters of regression models for different landforms (3 arc second and 30 arc second)

    地貌拟合模型拟合参数R2ΔmeanK
    abcd
    平原线性函数00.180 00.313 400.106 72.937 20
    多项式函数0.485 0−1.476 01.64600.322 200.324 40.993 20
    幂函数0.160 00.713 000.080 000.525 90.152 10
    丘陵线性函数00.246 00.552 500.385 11.434 70
    多项式函数0.449 0−1.695 02.33700.610 810.410 11.489 40
    幂函数0.082 10.465 000.077 800.579 90.134 10
    山地线性函数00.279 80.598 490.354 11.690 10
    多项式函数0.635 0−2.104 02.36600.687 300.716 00.959 90
    幂函数0.103 60.360 200.092 450.782 70.118 11
    下载: 导出CSV

    表  3  不同地貌单元不同分辨率下坡度转换公式

    Table  3.   The slope conversion formula of different resolution DEM for different landforms

    地貌单元类型与30″的转换关系
    1″(i=1″)3″(i=3″)
    平原$ {S_{30}} = {\text{0}}{\text{.085 63}}{S_{{i}}} $$ {S_{30}} = {\text{0}}{\text{.18}}{S_{{i}}} $
    丘陵$ {S_{30}} = {\text{0}}{\text{.132 57}}{S_{{i}}} $$ {S_{30}} = {\text{0}}{\text{.449}}{S_{{i}}}{{ - 1}}{\text{.695}}{S_{{i}}}^{\text{2}}{\text{ + 2}}{\text{.337}}{S_{{i}}}^{\text{3}} $
    山地$ {S_{30}} = {\text{0}}{\text{.199 24}}{S_{{i}}} $$ {S_{30}} = {\text{0}}{\text{.279 8}}{S_{{i}}} $
    下载: 导出CSV
  • [1] 李昕蕾, 2020. 基于多因素的区域场地分类方法研究. 哈尔滨: 中国地震局工程力学研究所.

    Li X. L., 2020. Regional site classification method based on multiple factors. Harbin: Institute of Engineering Mechanics, China Earthquake Administration. (in Chinese)
    [2] 刘飞, 范建容, 崔兆岩等, 2019. 基于DEM分形特征的坡度尺度变换模型. 山地学报, 37(1): 129—136

    Liu F., Fan J. R., Cui Z. Y., et al., 2019. A model of re-scaling slope based on DEM fractal feature. Mountain Research, 37(1): 129—136. (in Chinese)
    [3] 刘利峰, 2018. DEM分辨率对白马小流域坡度与坡向提取的影响. 山西水土保持科技, (2): 10—12 doi: 10.3969/j.issn.1008-0120.2018.02.005

    Liu L. F., 2018. DEM resolving capability influence on extract slope and aspect in Baima minor watershed. Soil and Water Conservation Science and Technology in Shanxi, (2): 10—12. (in Chinese) doi: 10.3969/j.issn.1008-0120.2018.02.005
    [4] 刘晓, 赵荣, 梁勇等, 2017. 顾及地貌与DEM分辨率的坡度算法适应性研究. 测绘科学, 42(3): 29—34

    Liu X., Zhao R., Liang Y., et al., 2017. Research of algorithms for deriving slope based on different geomorphic types and DEM resolution. Science of Surveying and Mapping, 42(3): 29—34. (in Chinese)
    [5] 乔之轩, 2019. 基于DEM的多分辨率水文特征提取. 北京: 北京林业大学.

    Qiao Z. X., 2019. Multi-resolution hydrological feature extraction based on DEM. Beijing: Beijing Forestry University. (in Chinese)
    [6] 师动, 朱奇峰, 杨勤科等, 2020. DEM分辨率对坡度算法选择影响的分析. 山地学报, 38(6): 935—944

    Shi D., Zhu Q. F., Yang Q. K., et al., 2020. DEM resolution influence on slope algorithm selection. Mountain Research, 38(6): 935—944. (in Chinese)
    [7] 史大成, 2009. 基于GIS的场地分类新方法研究. 哈尔滨: 中国地震局工程力学研究所.

    Shi D. C., 2009. Study on new methods of site classification based on GIS. Harbin: Institute of Engineering Mechanics, China Earthquake Administration. (in Chinese)
    [8] 土祥, 王春梅, 庞国伟等, 2018. 黄土丘陵沟壑区坡度尺度效应空间分异分析. 山地学报, 36(6): 964—972

    Tu X., Wang C. M., Pang G. W., et al., 2018. On the spatial differentiation of scale effect of slopes in Loess Hill and Gully Region, China. Mountain Research, 36(6): 964—972. (in Chinese)
    [9] 王英, 龚家国, 贾仰文等, 2019. 基于不同分辨率DEM提取坡度值的转换关系研究. 水利水电技术, 50(8): 45—51

    Wang Y., Gong J. G., Jia Y. W., et al., 2019. Study on conversion relationship of slope information extracted from different resolution DEM. Water Resources and Hydropower Engineering, 50(8): 45—51. (in Chinese)
    [10] 吴效勇, 2019. 基于多因素的地震影响场研究. 北京: 中国地震局地震预测研究所.

    Wu X. Y., 2019. Study on seismic influence field based on multi-factors. Beijing: Institute of Earthquake Forecasting, China Earthquake Administration. (in Chinese)
    [11] Grohmann C. H., 2015. Effects of spatial resolution on slope and aspect derivation for regional-scale analysis. Computers & Geosciences, 77: 111—117.
    [12] Horn B. K. P., 1981. Hill shading and the reflectance map. Proceedings of the IEEE, 69(1): 14—47. doi: 10.1109/PROC.1981.11918
    [13] Mukherjee S., Mukherjee S., Garg R. D., et al., 2013. Evaluation of topographic index in relation to terrain roughness and DEM grid spacing. Journal of Earth System Science, 122(3): 869—886. doi: 10.1007/s12040-013-0292-0
    [14] Olaya V., 2009. Basic land-surface parameters. Developments in Soil Science, 33: 141—169.
    [15] Thompson E. M., Wald D. J., Worden C. B., 2014. A Vs30 map for California with geologic and topographic constraints. Bulletin of the Seismological Society of America, 104(5): 2313—2321. doi: 10.1785/0120130312
    [16] Thornley J. H. M., 1976. Mathematical models in plant physiology. London: Academic Press, 86—110.
    [17] Wald D. J., Earle P. E., Lin K. W., et al., 2006. Challenges in rapid ground motion estimation for the prompt assessment of global urban earthquakes. Bulletin of the Earthquake Research Institute, 81: 273—281.
    [18] Walker J. P., Willgoose G. R., 1999. On the effect of digital elevation model accuracy on hydrology and geomorphology. Water Resources Research, 35(7): 2259—2268. doi: 10.1029/1999WR900034
    [19] Wang C. M., Yang Q. K., Guo W. L., et al., 2012. Influence of resolution on slope in areas with different topographic characteristics. Computers & Geosciences, 41: 156—168.
  • 加载中
图(7) / 表(3)
计量
  • 文章访问数:  154
  • HTML全文浏览量:  107
  • PDF下载量:  12
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-04-26
  • 刊出日期:  2021-06-30

目录

/

返回文章
返回