网站首页 公文大全 个人文档 实用范文 讲话致辞 实用工具 心得体会 哲学范文 总结范文 范文大全 报告 合同 文书 信函 实用
  • 汇报体会
  • 节日庆典
  • 礼仪
  • 毕业论文
  • 评语寄语
  • 导游词
  • 口号大全
  • 其他范文
  • 百花范文网 > 实用范文 > 其他范文 > 基于无人机测量的微地貌提取−以昆中断裂(巴隆-温泉段)为例1

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

    时间:2023-02-24 11:16:16来源:百花范文网本文已影响

    殷 翔 刘 炜 姚生海 盖海龙

    (青海省地震局, 西宁 810001)

    活动构造的定量研究需获取一系列表征活动构造特征的几何参数,如断裂长度、同震位移、滑动速率等(毕海芸等,2017)。突发的断裂运动会在近断层处保留丰富的断错地貌,获取地貌参数是研究断裂性质、位错量、活动速率、发震机理的重要基础(潘家伟等,2007)。

    为更精准地获取断错地貌变形数据,开展了大量活动构造研究,研究技术从早期的航卫片解译、皮尺测量、全站仪测量等,发展到目前广泛使用的差分GPS、LiDAR 及无人机低空摄影等,为地貌参数获取提供了更准确可靠的活动构造定量参数。相比于传统测量方法,近年来兴起的无人机摄影测量技术具有高清晰度、大比例尺、小面积、高现势性的优点,特别适合获取条状区域航拍影像。同时,该技术克服了传统测量方法工作量大、效率低、受自然条件限制等缺点,采用超低空飞行进行多视角拍摄,提供厘米级定位数据,从而显著提升图像元数据的绝对精度(罗国文等,2012),后期通过影像生成高精度的数字高程模型(DEM)和正射影像图(DOM),并提取复杂地貌面上的构造信息,为断裂研究提供定量化、精准化的基础数据。

    昆中断裂位于东昆仑隆升区内,研究程度较低,笔者发现昆中断裂在东段(巴隆-温泉段)存在长约110 km、线性特征明显的断层迹线,而在我国现行的第五代区划图中,将昆中断裂划分为7.0 级浅源,但缺乏基础资料作为潜源划分依据。因此,选择昆中断裂巴隆−温泉段作为研究区域,通过无人机航空摄影测量技术生成的高精度DEM 与DOM 数据,提取水平及垂直位错量等数据,为断裂定量研究提供数据基础。

    无人机摄影测量即通过无人机低空摄影获取高清晰影像数据,在具有多角度重叠的影像中,选择1 张影像作为主纠正影像,而其他影像作为从属影像,用于补偿主纠正影像上被遮挡部分的信息(马建等,2019),即在数字微分纠正过程中,以数字表面模型(DSM)为基础,进行数字微分纠正,进而生成三维点云与数字高程模型,实现地理信息的快速获取,并在计算机上使用相应软件可采集模型内任意需求的坐标数据(魏占玉等,2015)。

    1.1 数据采集原理

    选取大疆精灵4 RTK 小型多旋翼高精度航测无人机,采用后差分定位技术(PPK)对山前发育的断层陡坎及微地貌构造进行图像采集。采用PPK 技术时,基于1 台同步观测的基准站和至少1 台流动观测站,对GNSS 卫星进行同步观测,即基准站保持连续观测,流动站在测量过程中需保持对卫星的连续跟踪。基准站和流动站同步接收的数据在计算机中进行线性组合,形成虚拟的载波相位观测量,确定接收机之间的相对位置,最后引入基准站的已知坐标,从而获得流动站的三维坐标(刘卫刚等,2015;
    朱逍贤,2018)。

    1.2 数据处理技术

    在无人机低空摄影获得特定地点照片后,采用Photoscan 软件生成研究点DOM 和DEM 数据。处理过程主要分为以下步骤:

    (1)剔除原数据中的异常照片,避免生成的正射影像出现色彩异常;

    (2)导入照片后,Photoscan 软件将自动对齐照片并优化,在优化过程中自动估算校正参数;

    (3)获取相邻照片之间的重复点,并构建稀疏的点云数据;

    (4)创建网格,并重构几何图形;

    (5)根据预估的相机位置,分析其深度信息,并生成密集点云数据;

    (6)生成网格及纹理,并生成DEM 数据。

    研究区位于青藏高原中北部,区内活动断裂众多,周边发育多条区域性断裂,其中,以北部的昆中断裂带、东昆仑断裂带,中部的西藏大沟-昌马河断裂带、玛多-甘德断裂带、江错断裂,南部的巴颜喀拉主峰断裂带等。对于这些断裂,研究程度不一(姚生海等,2021)。该地区也是地震频发区,历史上曾发生过1902 年秀沟6.9 级地震、1963 年阿兰克湖7.0 级地震、1937 年花石峡7½级地震、1947 年达日7¾级地震、2021 年5 月22 日玛多7.4 级地震,表明巴颜喀拉块体仍为我国大陆强震的主体活动地区。

    昆中断裂西起昆仑山昆仑河附近,东延经大干沟、清水泉、青根河至鄂拉山温泉附近,被NNW 向鄂拉山断裂切割后,呈隐伏状态继续东延,与东昆仑断裂带近乎平行。鄂拉山以西,断裂构成东昆仑北坡断隆与柴达木南缘台缘褶带的分界。同时,昆中断裂是分割青藏高原甚至中国南北陆壳不同基底的分界线。断裂沿区内主要重力梯级带和磁力梯级带展布,其东段磁异常呈串珠状分布,与北秦岭南缘断裂相连,莫霍面沿断裂出现高达7~10 km 的陡坎和台阶,这是其他断裂没有的(曾秋生,1999)。

    研究区为昆中断裂东段,地理位置位于青海省玛多县苦海以北(图1)。断裂在山前形成了明显的断层陡坎,并发现了大量水系左旋,断层垭口、挤压隆起等构造地貌特征,由于该地区位于三江源保护区内,受人为因素的影响小,这些地貌特征保存完整,因此采用无人机摄影测量技术,有针对性地进行采集,快速获取高精度的地形地貌数据。

    图1 区域构造Fig. 1 Regional structure map

    在前期无人机低空摄影和后期数据处理的基础上,获取了研究区域内的正射影像。同时,基于高密度点云生成的DEM 数据分辨率较高,与断层相关的微地貌特征较明显。为清晰展示断裂的活动迹象,本文选取2 个典型地貌点进行详细解译,并通过实地野外调查验证此次无人机数据采集结果。

    本次飞行作业分别选取了龙通村北、乎勒塘村北作为本次工作的2 个试验点(以下简称为点1 和点2),飞行高度设定为80 m,分别获取419,489 张照片,照片重叠率达70%以上。其中,点1 处共有227 262 个连接点,生成密集点云11 784 833 个点,生成的DEM 数据分辨率为6 985 pix×4 552 pix,覆盖面积为0.549 km2;
    点2 处共有253 024 个连接点,生成密集点云12 423 887 个点,生成的DEM 数据分辨率为7 004 pix×4 233 pix,覆盖面积为0.325 km2。

    航高为航摄时飞机的飞行高度,根据起算基准的不同可分为绝对航高与相对航高。相对航高为无人机在飞行时相机相对于某一基准面的高度,是相对于作业区域内地面平均高程基准面的设计航高。一般来说,在其他因素不变的情况下,航高越低(飞机飞行越低),地面分辨率越高(朱海斌等,2018)。

    根据高建彩(2022)的规定,相对航高计算公式如下:

    式中,H为相对航高;
    f为摄影镜头的焦距;
    G为影像的地面分辨率;
    a为像元尺寸。

    通过换算可得到点1 的数据精度为2.07 cm/pix,点2 的数据精度为2.15 cm/pix,满足断裂位错量的计算要求。

    传统测量断裂的垂直位移方法主要为差分GPS,该方法获得的数据与DEM 数据相差较小,可满足工作需要,而测量断裂水平位移时主要使用皮尺等工具,受人为因素影响较大,野外地形相对复杂时测量数据无法真实反映断层位错量。为排除个人主观因素的影响,本文从高精度DEM 数据上直接提取陡坎高程数据,同时,剖面线的选取能够反映断层陡坎两侧的地形变化。

    3.1 龙通村北试验点(点1)

    点1 位于兴海县温泉乡西北约28 km 处,从遥感影像中可看出断裂在此处的线性特征明显(图2),呈NWW 向展布,长约700 m。断层陡坎主要发育在山前冲洪积扇体后缘(图3),水系有明显的左旋位错量(刘超等,2021)。

    图2 点1 遥感影像Fig. 2 Remote sensing image of point one

    图3 点1 断层陡坎Fig. 3 Fault scarp of point one

    通过野外地质调查发现,该处山前形成多期洪积扇,断裂在洪积扇后缘形成明显的断层陡坎、水系扭错及地震形成的跌水。断层陡坎连续性较好,多形成长约50 m、高0.5~1.1 m 的断层陡坎。在一处冲沟内,形成长约30 m、高2.5 m 的断层跌水。

    由于该地区海拔达4 500 m,人工测量微地貌难度较大,为更好地获取该试验点微地貌特征,利用无人机航拍及数字化处理技术,对该试验点进行高精度低空摄影,获取大量无人机数据,通过后期处理,生成DEM/DOM 矢量图(图4(a))。分析DEM/DOM 矢量图可知该处发育3 期洪积扇(图4(b)),断层横切Fan2 和Fan3 洪积扇。同时,调查发现,Fan2 洪积扇上的冲沟被左旋水平断错,形成了断头沟。

    基于生成的高分辨率DEM 数据,对3 条冲沟进行了位移重塑,其中,东侧断层水平位移仅2.1 m(图4(c)),并在冲沟后缘形成了高达2.5 m 的跌水;
    西侧断层水平位移可划分为2 期,最大左旋位移为15.4 m,最小位移为4.85 m(图4(d))。依据水平位移推测,认为西侧的冲沟也应存在1 期2 m 左右的水平位移,该位错量可能由最新一次断裂活动所致。

    图4 点1 基于DEM 数据的水系位错恢复Fig. 4 Water system dislocation recovery based on DEM data of point one

    昆中断裂在研究段不仅存在水平位错量,同时还具备一定的逆冲分量。根据此次获取的DEM 数据可知,在山前的冲洪积扇上,有1 条明显的延伸线(图5),其空间展布未受地形影响,通过实地调查,确认该延伸线为断裂在山前形成的断层陡坎。

    图5 点1 DEM 提取剖面位置示意Fig. 5 Schematic diagram of DEM extraction section position of point one

    利用本次生成的DEM 数据共获取4 条陡坎高度(图6~图9),陡坎高度基本为0.6~0.9 m,且4 条剖面线在断层陡坎两侧形态并无较大差异,原始地貌保存较好,笔者认为这4 条陡坎高度基本一致,应是同次断层活动形成,且陡坎高度可代表此次断裂活动的实际垂直位错量,而陡坎高度差异可能是断层陡坎形成后受水系侵蚀作用导致的。

    图6 点1 AA′地形剖面Fig. 6 AA" topographic profile of point one

    图7 点1 BB′地形剖面Fig. 7 BB" topographic profile of point one

    图8 点1 CC′地形剖面Fig. 8 CC" topographic profile of point one

    图9 点1 DD′地形剖面Fig. 9 DD" topographic profile of point one

    3.2 乎勒塘村北试验点(点2)

    点2 位于温泉乡西北约33 km 处,从遥感影像中可看出断裂在此处的线性特征较明显(图10),断裂呈NWW 向展布,陡坎长度约610 m,断层陡坎主要发育在山前冲洪积扇后缘,水系有明显的左旋位错量(图11(a))。点2 处的断层陡坎地貌特征明显(图11(b)),由图12(a)和图12(b)可知,山前至少发育了3 期洪积扇,断层横切Fan2 和Fan3 洪积扇。

    图10 点2 遥感影像Fig. 10 Remote sensing image of point two

    图11 点2 断层地貌特征Fig. 11 Geomorphological characteristics of fault in point two

    图12 点2 基于DEM 数据的水系位错恢复Fig. 12 Water system dislocation recovery based on DEM data of point two

    通过解译无人机数据发现,Fan2 洪积扇上的冲沟被左旋水平断错,并形成了2 处断头沟。通过对冲沟位错的恢复,图12(c)中的冲沟位错共发育2 期,2 期的左旋位错量分别为4.8 m 和12.1 m;
    图12(d)中的冲沟左旋位错量为5.2 m;
    与点1 处的水系左旋位错量基本一致。

    利用本次生成的DEM 数据共获取4 条陡坎高度(图13),陡坎高度基本为0.45~0.75 m(图14~图17),且3 条剖面线在断层陡坎两侧形态并无较大差异,原始地貌保存较好,笔者认为这3 条陡坎高度基本一致,应为同次断层活动形成,但与点1 相比,该处的垂直位错量明显较小,应是由于此处水流顺着陡坎方向汇集,使陡坎上侧的地貌侵蚀更剧烈,因此陡坎高度较小,因此,该处陡坎高度不能代表此次断裂活动的实际垂直位错量,点1 点处的陡坎高度更能代表断裂的实际位错量。

    图13 点2 DEM 提取剖面位置示意Fig. 13 Schematic diagram of DEM extraction section position

    图14 点2 AA′地形剖面Fig. 14 AA" topographic profile

    图15 点2 BB′地形剖面Fig. 15 BB" topographic profile

    图16 点2 CC′地形剖面Fig. 16 CC" topographic profile

    图17 点2 DD′地形剖面Fig. 17 DD" topographic profile

    本文采用大疆精灵4 RTK 无人机测绘影像进行3 D 模型构建,能够较准确地还原实际地貌场景,获得较高分辨率的DEM 数据。通过对采集图像进行处理,可进行断层地貌变形的精确测量,为研究微地貌特征、断裂几何展布等提供方便有效的方法。

    本文利用无人机摄影测量技术获取了断裂局部地区的微地貌数据−断层陡坎及水系左旋的高精度DEM,对研究区域内断裂水平、垂直位错量进行精细研究,确定了昆中断裂(巴隆-温泉段)活动断层陡坎位置及断裂最新活动量,得出以下结论:

    (1)无人机航拍技术相较于传统测量方式,具有低成本、易操作、获取数据方便的特点。同时,利用PhotoScan 进行数据处理时,数据处理速度快、误差小,能够满足断层陡坎精细化测量的数据要求,凸显无人机航拍技术的应用潜力。但该技术也受天气、飞行高度及续航等因素限制,飞行数据易产生异常,不利于后续工作的开展。

    (2)利用无人机获取的高精度DEM 数据,分别得到昆中断裂(巴隆-温泉段)水平位错量为2.1、4.85、15.4 m,位错量之间近似存在着倍数关系,而15.4 m 的位错量可能是由4 次断层活动形成的累计位错量。因此,在下一步工作中,还应探明是否存在1 期8 m 左右的断层累计活动位错量。

    (3)通过对比2 处试验点DEM 数据提取的8 条陡坎数据,陡坎两侧的地貌形态基本相似,无明显差异,乎勒塘村北试验点处的陡坎高度受水流侵蚀较严重,而龙通村北试验点处的陡坎高度均为0.6~0.9 m,可初步认为这是最新一次断层运动形成的陡坎,可代表断裂实际位错量。

    猜你喜欢 冲沟左旋剖面 ATC系统处理FF-ICE四维剖面的分析中国科技纵横(2022年17期)2022-10-25基于历史卫星影像估算黄土丘陵区冲沟发育速率农业工程学报(2022年1期)2022-03-25左旋的柳布达拉(2019年3期)2019-06-11高原湿陷性黄土地区冲沟高填方路基施工工艺及质量控制价值工程(2017年35期)2018-02-08复杂多约束条件通航飞行垂直剖面规划方法北京航空航天大学学报(2016年6期)2016-11-16船体剖面剪流计算中闭室搜索算法船舶标准化工程师(2015年5期)2015-12-03左旋肉碱对肥胖患者血脂水平影响的meta分析郑州大学学报(医学版)(2015年2期)2015-02-27长输管道定向钻穿越高原地区冲沟、陡坡技术企业文化·中旬刊(2014年12期)2014-12-30罗哌卡因与左旋布比卡因复合舒芬太尼用于分娩镇痛效果比较中国药业(2014年17期)2014-05-26不同时间段服用左旋氨氯地平治疗老年非杓型高血压患者31例中国药业(2014年21期)2014-05-26

    相关热词搜索:无人机 地貌 为例

    • 范文大全
    • 说说大全
    • 学习资料
    • 语录
    • 生肖
    • 解梦
    • 十二星座