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

基于无人机测量的微地貌提取−以昆中断裂(巴隆-温泉段)为例

殷翔 刘炜 姚生海 盖海龙

马兴全, 王志铄, 张扬. 基于OLI和DEM数据初步研究盘谷寺断裂新构造运动特征[J]. 震灾防御技术, 2020, 15(2): 359-368. doi: 10.11899/zzfy20200213
引用本文: 殷翔,刘炜,姚生海,盖海龙,2022. 基于无人机测量的微地貌提取−以昆中断裂(巴隆-温泉段)为例. 震灾防御技术,17(4):701−709. doi:10.11899/zzfy20220410. doi: 10.11899/zzfy20220410
Ma Xingquan, Wang Zhishuo, Zhang Yang. Preliminary Study on Neotectonic Movement Characteristics of Pangusi Fault Based on OLI and DEM Data[J]. Technology for Earthquake Disaster Prevention, 2020, 15(2): 359-368. doi: 10.11899/zzfy20200213
Citation: Yin Xiang, Liu Wei, Yao Shenghai, Gai Hailong. Micro Landform Extraction Based on UAV Photography Technology−Taking Kunzhong Fault (Balong Wenquan Section) as An Example[J]. Technology for Earthquake Disaster Prevention, 2022, 17(4): 701-709. doi: 10.11899/zzfy20220410

基于无人机测量的微地貌提取−以昆中断裂(巴隆-温泉段)为例

doi: 10.11899/zzfy20220410
基金项目: 中国地震局星火计划项目(XH21038Y);青海省地震科学基金(2019A07);青海省地震灾害风险普查项目
详细信息
    作者简介:

    殷翔,男,生于1989年。工程师。主要从事地震构造等研究。E-mail:544261740@qq.com

Micro Landform Extraction Based on UAV Photography Technology−Taking Kunzhong Fault (Balong Wenquan Section) as An Example

  • 摘要: 无人机测量具有高清晰度、大比例尺、小面积、高现势性的优点,为地貌参数获取提供更准确可靠的活动构造定量参数,克服传统测量方法工作量大、效率低、受自然条件限制等缺点,可提供厘米级定位数据,从而显著提升图像元数据的绝对精度。利用大疆精灵4 RTK小型多旋翼高精度航测无人机,获取昆中断裂(巴隆-温泉段)在龙通村北的高精度DEM地貌数据,通过对微地貌的提取,初步确定断裂在该处的水平位错量为2.1~15.4 m。分析获取的8条陡坎剖面,认为其中5条陡坎形成后受到水流侵蚀作用较小,陡坎高度基本相似,断裂实际垂直位错量为0.6~0.9 m。研究结果表明,无人机航测技术是识别复杂地貌构造信息并提取相关活动构造参数的有效手段,可为断裂的定量研究提供可靠的数据基础。
  • 近年来,由于遥感技术的不断进步,遥感地质学在矿产普查、区域地质和水文地质等方面开始发挥日益重要的作用(赵鸿燕等,2010齐信等,2012魏永明等,2015段瑞琪等,2017)。通过遥感影像识别断裂各类构造地貌信息,进而讨论断裂新构造运动特征的研究方法已颇为成熟(贾营营等,2010郑颖平等,2017)。目前,随着对构造地貌理解的不断深入,众多学者已通过不同构造地貌类型的组合和发育现状反演构造活动历史,开展更全面的运动学研究(刘华国等,2011文力等,2018)。

    断裂构造地貌是由断裂活动在地表塑造的地貌,由于断裂活动方式多样,在地表塑造的地貌类型也不尽相同(杨景春等,2011)。在各类地貌中,线性特征是断裂最明显的构造标志,也是断裂遥感解译中最重要的地物之一。通常情况下,走滑断裂的线性特征较倾滑断裂更加典型。此外,断裂垂向错动还可形成断层陡坎(中屠炳明等,1991)、台地(师亚芹等,2009)及河流裂点(毕丽思等,2011)等地貌;断裂水平扭动错断山脊线和水系,可形成眉脊和断尾沟等地貌。这些地貌的形成与发育和断裂活动方式、强度、频度、岩性特征、地形、气候等因素密切相关。通过构造地貌的识别与分析,可开展断裂空间几何学、运动学等方面的研究。随着高分辨率遥感影像的逐步应用和GIS技术的不断提高,结合断裂构造地貌特征与现有地质资料,有助于提高断裂解译的精度和深度。

    盘谷寺断裂为太行山南缘盘谷寺-新乡断裂的西段,是我国东部1条重要的岩石圈构造带(邓起东等,1980马杏垣,1989张培震,1999)、地壳厚度陡变带(江娃利等,1984徐杰等,2000)和地震活动带(李清武,2015),同时也是太行山断块隆起和南华北坳陷的分界断裂。盘谷寺断裂活动始于燕山运动,于燕山运动晚期—喜马拉雅运动中期强烈活动,最终止于中更新世(中国地震局地球物理勘探中心,2015)。目前针对盘谷寺断裂新构造运动特征的深入分析仍较少,分析该断裂新构造运动特征对研究太行山东南缘构造地貌演化,乃至我国东部新构造活动格局具有重要意义。本文综合利用多波段与全色波段遥感数据,在ArcGIS平台上,采用陆地成像仪(Operational Land Imager,OLI)数据及先进星载热发射和反射辐射仪全球数字高程模型(Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model,ASTER GDEM)数据,通过地物判读、提取坡向坡度及构建三维地貌模型等手段,对盘谷寺断裂进行遥感解译,确定该断裂空间展布形态,并利用地壳掀斜规律和剖面分析,初步总结断裂活动特征,取得良好效果。

    研究区地处河南省西北部,与山西省相邻,构造上位于太行山断块隆起与南华北坳陷的交界部位(图 1)。南太行基底由太古界登封群和与之不整合接触的下元古界嵩山群组成,表现为紧闭线型复式褶皱;盖层为中元古界汝阳群和寒武系-中奥陶统,多呈单斜构造和开阔褶皱(河南省地质矿产局,1989)。

    图 1  研究区地震构造(红框为研究区)
    Figure 1.  Seismotectonic map of the study area(the red box denothing the study area)
    F1盘谷寺断裂;F2武陟断裂;F3封门口断裂;F4郑州—开封断裂

    根据已有成果,在太行山南缘发育的1组近东西走向、陡倾的阶梯状正断层控制了本地区构造地貌的形成和发育(马寅生等,2007)。自新生代以来,本地区共经历3次快速的构造隆升运动,在间歇期于海拔1800m、1100—1300m和350—500m处分别形成了北台期、太行期和唐县期三级夷平面(程绍平等,1989吴忱等,1999)。进入第四纪后,断裂活动性减弱,仅山前断裂——盘谷寺断裂错断了第四纪沉积物(马寅生等,2007)。

    盘谷寺-新乡断裂为太行山隆起与济源-开封凹陷的分界断裂,该断裂大致以柏山、大高村为界,由西向东划分为3个次级段落。西段位于柏山以西,为盘谷寺断裂,中段和东段均为隐伏断裂。盘谷寺断裂南部济源凹陷的西界和东界分别为封门口断裂和武陟断裂(图 1),与凹陷内部发育的多条近东西向次级正断层共同控制了凹陷整体的形成和发展过程。该凹陷在古近纪时期强烈断陷下沉,南北向伸展率可达35%—60%(潘澄雨,2013),内部堆积了巨厚的古近系。进入新近纪和第四纪后,凹陷伸展断陷作用显著减弱(马寅生等,2007)。

    盘谷寺-新乡断裂西段为中更新世活动段,中段和东段为晚更新世活动段(中国地震局地球物理勘探中心,2015)。断裂附近曾于1587年发生2次破坏性地震,震级分别为5.5级和6.0级。现今小震活动主要集中于南部的济源-开封凹陷内。

    遥感影像数据选用Landsat 8 OLI陆地成像仪和ASTER GDEM数据。Landsat 8 OLI陆地成像仪发射于2013年,共9个波段,包括1个分辨率15m的全色波段,空间分辨率为30m,成像宽幅185km×185km,选取的影像云量均小于1%。

    高程数据选用ASTER GDEM数据,该数据发布于2011年,其空间分辨率为30m,垂直精度约20m,每幅影像覆盖面积约60km×60km(ASTER GDEM Validation Team,2009)。

    本文所用遥感影像数据和高程数据均来源于地理空间数据云,根据该网站提供的产品信息,Landsat 8 OLI影像已经系统辐射校正和几何校正,故本次数据处理主要包括影像镶嵌、融合与裁剪等(图 2)。

    图 2  数据处理流程图
    Figure 2.  Flow chart of data processing

    遥感影像的处理在ENVI环境下进行,首先将选取的2幅遥感影像各波段数据进行拼接,然后将OLI数据多光谱波段和全色波段进行融合处理,通过重采样生成分辨率为15m的多光谱数据,采用6、5、4波段进行假彩色合成,通过线性拉伸、边缘增强等处理突出纹理特征,地质信息丰富,视觉效果得到明显提升。ASTER GDEM数据的镶嵌和裁剪处理均在ArcGIS下进行。最终将OLI遥感影像与DEM数据叠加,构建地区三维地貌模型。

    流水型构造地貌对构造活动的响应十分灵敏(史兴民等,2003),构造运动通过改变地表地形高差、地壳掀斜特征等影响冲沟和河谷的发育及形态。在气候、岩性相近的条件下,构造运动是影响冲沟形态最重要的因素。正断裂发生差异升降运动后,上盘相对下降,局部侵蚀基准面随之降低,由下盘而来的水流将在沟谷中产生强烈的溯源侵蚀作用,最终在下盘发育一系列与断裂近垂直的冲沟。断裂作为溯源侵蚀作用的起点,其差异升降运动的强弱往往决定着溯源侵蚀作用的强弱。反之,一系列近平行的冲沟可作为断裂解译的标志。

    通过判读盘谷寺断裂的遥感影像,发现沿盘谷寺断裂走向,在断裂下盘发育的冲沟形态差别十分明显,大致以济源市河口村为界,在河口村以西,山前发育一系列密集、近平行、长约数百米的深切冲沟,冲沟沟尾线性相连,如图 3(a)所示,是判断断裂位置的良好标志;在河口村以东,断裂线性特征较清晰,影像中显示清晰的断层三角面及大型山间沟谷,如图 3(b)中红框所示;沿走向继续向东,断层三角面形态逐渐消失,断裂下盘表现为经大面积强烈侵蚀作用后的低矮残丘,如图 3(c)所示。残丘与断裂上盘仍存在一定高差,以致断裂在影像上虽依稀表现出线性特征,但较断裂西段已明显弱化。

    图 3  盘谷寺断裂下盘地貌表现
    Figure 3.  Landscape map on the lower walls of Pangusi fault

    基于以上解译认识,基本可确定断裂的空间展布方式。由西向东,在断裂下盘依次发育深切冲沟→断层三角面→侵蚀残丘,客观反映了溯源侵蚀作用由西向东逐渐减弱的规律,据此初步判断断裂差异升降运动由西向东逐渐减弱。

    坡状地貌也是构造地貌学研究的重要内容之一,坡状地貌的形成演化是构造运动和侵蚀作用长期共同作用的结果(Davis,1899杨景春等,2011)。构造运动是坡状地貌的诱发因素,具有脉冲式特点,而侵蚀作用则不断将坡状地貌夷平,具有长期性。通过评价坡状地貌形态可反演地区历经的构造运动强度和规模。坡向和坡度是定量描述坡状地貌的2个重要参数。在ArcGIS平台上,采用ASTER GDEM数据提取了盘谷寺断裂沿线坡向、坡度图,如图 4所示。由图 4可知,由于构造运动引起的地壳掀斜和后期侵蚀堆积作用,盘谷寺断裂上盘地表普遍向南倾斜,倾斜区域呈楔形分布,宽度由西向东逐渐变小,至博爱县北附近,地表倾斜现象已不明显,如图 4(a)所示。由坡度图可知,由于断裂的存在,太行山脉和济源凹陷间存在1条地形坡度分界线,该界线以北为太行山脉,地形坡度大,起伏多变;以南为济源凹陷,地形平缓。地形坡度的差异由断裂活动造成,该地形坡度分界线即为盘谷寺断裂的构造位置,断裂活动造成两侧地形坡度产生差异。

    图 4  盘谷寺断裂两侧坡向、坡度图
    Figure 4.  Aspect and slope on both sides of Pagusi fault

    通过对比坡向图和坡度图可知,断裂西段地壳掀斜范围、地形倾角变化率均宽于和大于断裂东段,盘谷寺断裂西段对地表地形的影响强于东段,这是由于断裂西段的活动性较强。

    分析地形剖面可获得断裂垂向断错特征,地形剖面形态主要与构造活动、侵蚀及堆积过程有关,通过识别地形剖面中的夷平面、断面及洪积扇堆积方式等要素,可研究断裂的活动期次与活动强度等特征(杨源源等,2012)。

    笔者提取了4条近垂直穿过盘谷寺断裂的地形剖面(图 5),这4条剖面由西向东依次穿过了该断裂的西段和东段,充分反映了断裂整体的垂向断错特征。

    图 5  盘谷寺断裂地形剖面
    Figure 5.  Topographic section across Pangusi fault

    A—A'剖面位于济源市以西,剖面显示,断裂下盘的太行山脉海拔普遍为700—1000m,沟壑纵横,地形起伏大;断裂活动在地表形成了数百米高且陡直的断层面;在断裂上盘,由于靠近断裂的位置存在洪积扇体堆积,造成局部地形坡度较大,向南地形平缓,海拔较低,普遍不超过200m。

    B—B'剖面位于济源市以东,与A—A'剖面类似,断裂两侧地形迥异,下盘海拔普遍为700—1000m,上盘地形平坦且海拔不超过200m,断裂活动在地表形成了数百米高且陡直的断层面。

    C—C'剖面位于沁阳市西侧,跨越盘谷寺断裂东段,剖面北部为太行山脉,海拔由900m向南逐渐降至300m,断裂活动未在地表形成高耸的断层面,而是表现为斜坡状地形,斜坡顶部与底部高差仅数十米,该剖面反映断裂在该处的活动性变弱,大面积侵蚀和堆积作用开始在塑造地形的过程中占据主导地位。

    D—D'剖面位于盘谷寺断裂东端,剖面整体形态与C—C'剖面类似,由于断裂活动性弱,海拔高度由北部太行山向济源凹陷连续过渡,地形剖面几乎已无法表现出断面形态。

    结合地形剖面与前文所述冲沟发育、坡向坡度特征,认为盘谷寺断裂差异升降运动由西向东逐渐减弱,因此济源凹陷西部的沉降幅度较东部大,西部新近系底界埋深也大于东部,如图 1所示(河南省地质矿产局,1989王志铄,2017)。

    将OLI数据和高程数据叠加,可获得盘谷寺断裂三维地貌模型。叠加后的三维地貌模型较原遥感影像更立体。通过增加阴影、扩大Z值因子等手段,对盘谷寺断裂地貌细节进行突出显示(图 6)。

    图 6  盘谷寺断裂三维地貌模型
    Figure 6.  Three-dimensional geomorphic model of Pangusi fault

    图 6可知,经过三维地貌模型突出显示后,断裂两盘地形地貌反差明显,通过突出细节,使地表迹线更清晰,便于研究断裂的空间展布方式和表现形式。该模型同时直观地展示出盘谷寺断裂下盘由西向东依次发育了不同的构造地貌:断裂西段下盘表现为深切侵蚀沟,断裂由沟尾附近穿过;向东发展,断裂下盘表现为较完整的断层三角面;至断裂东段,由于断裂活动性减弱,以大面积侵蚀作用为主,断裂下盘主要表现为低矮的侵蚀残丘。

    本文采用Landsat 8 OLI数据和ASTER GDEM数据,以ArcGIS为平台,对太行山南缘盘谷寺断裂开展了构造地貌提取和分析,通过分析断裂下盘冲沟发育特征、断裂两盘坡向和坡度变化规律及跨断层地形剖面等手段,判定断裂空间展布形态,得到以下结论:

    (1) 盘谷寺断裂下盘由西向东依次发育了深切侵蚀(冲)沟、断层三角面和低矮的侵蚀残丘3类构造地貌,这是伴随构造活动强度的减弱,侵蚀方式由强烈的溯源侵蚀逐渐向大面积侵蚀作用转变的结果。

    (2) 不同类型的构造地貌表现与地形剖面图共同说明盘谷寺断裂为1条倾向南的正断层,断裂差异升降运动由西向东逐渐减弱,导致新近系底界在断裂西部的埋深大于东部,至中更新世时期断裂停止活动。

    (3) 断裂西部的沉降幅度大于东部,说明该断裂新构造运动以来的差异升降运动不均匀,西部差异升降运动强于东部,因此推断太行山在隆升剥露的过程中并非同步抬升,值得进一步深入研究。

  • 图  1  区域构造

    Figure  1.  Regional structure map

    图  2  点1遥感影像

    Figure  2.  Remote sensing image of point one

    图  3  点1断层陡坎

    Figure  3.  Fault scarp of point one

    图  4  点1基于DEM数据的水系位错恢复

    Figure  4.  Water system dislocation recovery based on DEM data of point one

    图  5  点1 DEM提取剖面位置示意

    Figure  5.  Schematic diagram of DEM extraction section position of point one

    图  6  点1 AA′地形剖面

    Figure  6.  AA' topographic profile of point one

    图  7  点1 BB′地形剖面

    Figure  7.  BB' topographic profile of point one

    图  8  点1 CC′地形剖面

    Figure  8.  CC' topographic profile of point one

    图  9  点1 DD′地形剖面

    Figure  9.  DD' topographic profile of point one

    图  10  点2遥感影像

    Figure  10.  Remote sensing image of point two

    图  11  点2断层地貌特征

    Figure  11.  Geomorphological characteristics of fault in point two

    12  点2基于DEM数据的水系位错恢复

    12.  Water system dislocation recovery based on DEM data of point two

    图  13  点2 DEM提取剖面位置示意

    Figure  13.  Schematic diagram of DEM extraction section position

    图  14  点2 AA′地形剖面

    Figure  14.  AA' topographic profile

    图  15  点2 BB′地形剖面

    Figure  15.  BB' topographic profile

    图  16  点2 CC′地形剖面

    Figure  16.  CC' topographic profile

    图  17  点2 DD′地形剖面

    Figure  17.  DD' topographic profile

  • 毕海芸, 郑文俊, 曾江源等, 2017. SfM摄影测量方法在活动构造定量研究中的应用. 地震地质, 39(4): 656—674 doi: 10.3969/j.issn.0253-4967.2017.04.003

    Bi H. Y. , Zheng W. J. , Zeng J. Y. , et al. , 2017. Application of SfM photogrammetry method to the quantitative study of active tectonics. Seismology and Geology, 39(4): 656—674. (in Chinese) doi: 10.3969/j.issn.0253-4967.2017.04.003
    高建彩, 2022. 倾斜摄影测量在智慧城市建设中的应用. 黑龙江科学, 13(6): 91—93 doi: 10.3969/j.issn.1674-8646.2022.06.027

    Gao J. C. , 2022. Application of tilt photogrammetry in smart city construction. Heilongjiang Science, 13(6): 91—93. (in Chinese) doi: 10.3969/j.issn.1674-8646.2022.06.027
    刘超, 雷启云, 余思汗等, 2021. 基于无人机摄影测量技术的地震地表破裂带定量参数提取——以1709年中卫南M7½地震为例. 地震学报, 43(1): 113—123 doi: 10.11939/jass.20200039

    Liu C. , Lei Q. Y. , Yu S. H. , et al. , 2021. Using UAV photogrammetry technology to extract the Quanti-Tative parameters of earthquake surface rupture zone: a case study of the southern Zhongwei M7½ earthquake in 1709. Acta Seismologica Sinica, 43(1): 113—123. (in Chinese) doi: 10.11939/jass.20200039
    刘卫刚, 马静, 2015. 浅谈GPS-PPK技术在1: 1万基础测绘中的应用. 测绘技术装备, 17(2): 80—81. doi: 10.3969/j.issn.1674-4950.2015.02.027
    罗国文, 阴志宏, 杨树文, 2012. 断裂构造遥感识别和提取方法的现状与展望. 山东国土资源, 28(2): 29—33 doi: 10.3969/j.issn.1672-6979.2012.02.010

    Luo G. W. , Yin Z. H. , Yang S. W. , 2012. Present condition and prospect of ertracton and identification methods by using remote sensing technology in faults. Shandong Land and Resources, 28(2): 29—33. (in Chinese) doi: 10.3969/j.issn.1672-6979.2012.02.010
    马建, 黄帅堂, 吴国栋, 2019. 基于微型无人机摄影测量技术的微构造信息提取研究——以博-阿断裂乌苏通沟东岸为例. 中国地震, 35(3): 550—557 doi: 10.3969/j.issn.1001-4683.2019.03.013

    Ma J. , Huang S. T. , Wu G. D. , 2019. Microstructure information extraction based on the UAV photogrammetry technology-a case study of the bolokenu-aqikekuduke fault on the east bank of Wusutonggou. Earthquake Research in China, 35(3): 550—557. (in Chinese) doi: 10.3969/j.issn.1001-4683.2019.03.013
    潘家伟, 李海兵, Van Der Woerd J. 等, 2007. 西昆仑山前冲断带晚新生代构造地貌特征. 地质通报, 26(10): 1368—1379 doi: 10.3969/j.issn.1671-2552.2007.10.014

    Pan J. W. , Li H. B. , Van Der Woerd J. , et al. , 2007. Late cenozoic morphotectonic features of the thrust belt in the front of the west Kunlun mountains. Geological Bulletin of China, 26(10): 1368—1379. (in Chinese) doi: 10.3969/j.issn.1671-2552.2007.10.014
    魏占玉, Ramon A. , 何宏林等, 2015. 基于SfM方法的高密度点云数据生成及精度分析. 地震地质, 37(2): 636—648 doi: 10.3969/j.issn.0253-4967.2015.02.024

    Wei Z. Y. , Ramon A. , He H. L. , et al. , 2015. Accuracy analysis of terrain point cloud acquired by "Structure from Motion" using aerial photos. Seismology and Geology, 37(2): 636—648. (in Chinese) doi: 10.3969/j.issn.0253-4967.2015.02.024
    姚生海, 盖海龙, 殷翔等, 2021. 青海玛多MS7.4地震地表破裂带的基本特征和典型现象. 地震地质, 43(5): 1060—1072

    Yao S. H. , Gai H. L. , Yin X. , et al. , 2021. The basic characteristics and typical phenomena of the surface rupture zone of the Maduo MS7.4 earthquake in Qinghai. Seismology and Geology, 43(5): 1060—1072. (in Chinese)
    曾秋生, 1999. 青海地震综合研究. 北京: 地震出版社.
    朱海斌, 王妍, 李亚梅, 2018. 基于无人机的露天矿区测绘研究. 煤炭工程, 50(10): 162—166

    Zhu H. B. , Wang Y. , Li Y. M. , 2018. Surveying and mapping of open-pit mining area using unmanned aerial vehicle. Coal Engineering, 50(10): 162—166. (in Chinese)
    朱逍贤, 2018. 基于无人机倾斜摄影制作真正射影像研究. 成都: 成都理工大学.

    Zhu X. X., 2018. The study of the real projective image based on TDOM. Chengdu: Chengdu University of Technology. (in Chinese)
  • 加载中
图(18)
计量
  • 文章访问数:  162
  • HTML全文浏览量:  41
  • PDF下载量:  32
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-03-28
  • 刊出日期:  2022-12-31

目录

/

返回文章
返回