殷 翔 劉 煒 姚生海 蓋海龍

（青海省地震局, 西寧 810001）

引言

活動(dòng)構(gòu)造的定量研究需獲取一系列表征活動(dòng)構(gòu)造特征的幾何參數(shù)，如斷裂長(zhǎng)度、同震位移、滑動(dòng)速率等（畢海蕓等，2017）。突發(fā)的斷裂運(yùn)動(dòng)會(huì)在近斷層處保留豐富的斷錯(cuò)地貌，獲取地貌參數(shù)是研究斷裂性質(zhì)、位錯(cuò)量、活動(dòng)速率、發(fā)震機(jī)理的重要基礎(chǔ)（潘家偉等，2007）。

為更精準(zhǔn)地獲取斷錯(cuò)地貌變形數(shù)據(jù)，開(kāi)展了大量活動(dòng)構(gòu)造研究，研究技術(shù)從早期的航衛(wèi)片解譯、皮尺測(cè)量、全站儀測(cè)量等，發(fā)展到目前廣泛使用的差分GPS、LiDAR 及無(wú)人機(jī)低空攝影等，為地貌參數(shù)獲取提供了更準(zhǔn)確可靠的活動(dòng)構(gòu)造定量參數(shù)。相比于傳統(tǒng)測(cè)量方法，近年來(lái)興起的無(wú)人機(jī)攝影測(cè)量技術(shù)具有高清晰度、大比例尺、小面積、高現(xiàn)勢(shì)性的優(yōu)點(diǎn)，特別適合獲取條狀區(qū)域航拍影像。同時(shí)，該技術(shù)克服了傳統(tǒng)測(cè)量方法工作量大、效率低、受自然條件限制等缺點(diǎn)，采用超低空飛行進(jìn)行多視角拍攝，提供厘米級(jí)定位數(shù)據(jù)，從而顯著提升圖像元數(shù)據(jù)的絕對(duì)精度（羅國(guó)文等，2012），后期通過(guò)影像生成高精度的數(shù)字高程模型（DEM）和正射影像圖（DOM），并提取復(fù)雜地貌面上的構(gòu)造信息，為斷裂研究提供定量化、精準(zhǔn)化的基礎(chǔ)數(shù)據(jù)。

昆中斷裂位于東昆侖隆升區(qū)內(nèi)，研究程度較低，筆者發(fā)現(xiàn)昆中斷裂在東段（巴隆-溫泉段）存在長(zhǎng)約110 km、線性特征明顯的斷層跡線，而在我國(guó)現(xiàn)行的第五代區(qū)劃圖中，將昆中斷裂劃分為7.0 級(jí)淺源，但缺乏基礎(chǔ)資料作為潛源劃分依據(jù)。因此，選擇昆中斷裂巴隆?溫泉段作為研究區(qū)域，通過(guò)無(wú)人機(jī)航空攝影測(cè)量技術(shù)生成的高精度DEM 與DOM 數(shù)據(jù)，提取水平及垂直位錯(cuò)量等數(shù)據(jù)，為斷裂定量研究提供數(shù)據(jù)基礎(chǔ)。

1 數(shù)據(jù)采集原理及流程

無(wú)人機(jī)攝影測(cè)量即通過(guò)無(wú)人機(jī)低空攝影獲取高清晰影像數(shù)據(jù)，在具有多角度重疊的影像中，選擇1 張影像作為主糾正影像，而其他影像作為從屬影像，用于補(bǔ)償主糾正影像上被遮擋部分的信息（馬建等，2019），即在數(shù)字微分糾正過(guò)程中，以數(shù)字表面模型（DSM）為基礎(chǔ)，進(jìn)行數(shù)字微分糾正，進(jìn)而生成三維點(diǎn)云與數(shù)字高程模型，實(shí)現(xiàn)地理信息的快速獲取，并在計(jì)算機(jī)上使用相應(yīng)軟件可采集模型內(nèi)任意需求的坐標(biāo)數(shù)據(jù)（魏占玉等，2015）。

1.1 數(shù)據(jù)采集原理

選取大疆精靈4 RTK 小型多旋翼高精度航測(cè)無(wú)人機(jī)，采用后差分定位技術(shù)（PPK）對(duì)山前發(fā)育的斷層陡坎及微地貌構(gòu)造進(jìn)行圖像采集。采用PPK 技術(shù)時(shí)，基于1 臺(tái)同步觀測(cè)的基準(zhǔn)站和至少1 臺(tái)流動(dòng)觀測(cè)站，對(duì)GNSS 衛(wèi)星進(jìn)行同步觀測(cè)，即基準(zhǔn)站保持連續(xù)觀測(cè)，流動(dòng)站在測(cè)量過(guò)程中需保持對(duì)衛(wèi)星的連續(xù)跟蹤?；鶞?zhǔn)站和流動(dòng)站同步接收的數(shù)據(jù)在計(jì)算機(jī)中進(jìn)行線性組合，形成虛擬的載波相位觀測(cè)量，確定接收機(jī)之間的相對(duì)位置，最后引入基準(zhǔn)站的已知坐標(biāo)，從而獲得流動(dòng)站的三維坐標(biāo)（劉衛(wèi)剛等，2015；朱逍賢，2018）。

1.2 數(shù)據(jù)處理技術(shù)

在無(wú)人機(jī)低空攝影獲得特定地點(diǎn)照片后，采用Photoscan 軟件生成研究點(diǎn)DOM 和DEM 數(shù)據(jù)。處理過(guò)程主要分為以下步驟：

（1）剔除原數(shù)據(jù)中的異常照片，避免生成的正射影像出現(xiàn)色彩異常；

（2）導(dǎo)入照片后，Photoscan 軟件將自動(dòng)對(duì)齊照片并優(yōu)化，在優(yōu)化過(guò)程中自動(dòng)估算校正參數(shù)；

（3）獲取相鄰照片之間的重復(fù)點(diǎn)，并構(gòu)建稀疏的點(diǎn)云數(shù)據(jù)；

（4）創(chuàng)建網(wǎng)格，并重構(gòu)幾何圖形；

（5）根據(jù)預(yù)估的相機(jī)位置，分析其深度信息，并生成密集點(diǎn)云數(shù)據(jù)；

（6）生成網(wǎng)格及紋理，并生成DEM 數(shù)據(jù)。

2 研究區(qū)地震構(gòu)造綜述

研究區(qū)位于青藏高原中北部，區(qū)內(nèi)活動(dòng)斷裂眾多，周邊發(fā)育多條區(qū)域性斷裂，其中，以北部的昆中斷裂帶、東昆侖斷裂帶，中部的西藏大溝-昌馬河斷裂帶、瑪多-甘德斷裂帶、江錯(cuò)斷裂，南部的巴顏喀拉主峰斷裂帶等。對(duì)于這些斷裂，研究程度不一（姚生海等，2021）。該地區(qū)也是地震頻發(fā)區(qū)，歷史上曾發(fā)生過(guò)1902 年秀溝6.9 級(jí)地震、1963 年阿蘭克湖7.0 級(jí)地震、1937 年花石峽7?級(jí)地震、1947 年達(dá)日7?級(jí)地震、2021 年5 月22 日瑪多7.4 級(jí)地震，表明巴顏喀拉塊體仍為我國(guó)大陸強(qiáng)震的主體活動(dòng)地區(qū)。

昆中斷裂西起昆侖山昆侖河附近，東延經(jīng)大干溝、清水泉、青根河至鄂拉山溫泉附近，被NNW 向鄂拉山斷裂切割后，呈隱伏狀態(tài)繼續(xù)東延，與東昆侖斷裂帶近乎平行。鄂拉山以西，斷裂構(gòu)成東昆侖北坡斷隆與柴達(dá)木南緣臺(tái)緣褶帶的分界。同時(shí)，昆中斷裂是分割青藏高原甚至中國(guó)南北陸殼不同基底的分界線。斷裂沿區(qū)內(nèi)主要重力梯級(jí)帶和磁力梯級(jí)帶展布，其東段磁異常呈串珠狀分布，與北秦嶺南緣斷裂相連，莫霍面沿?cái)嗔殉霈F(xiàn)高達(dá)7～10 km 的陡坎和臺(tái)階，這是其他斷裂沒(méi)有的（曾秋生，1999）。

研究區(qū)為昆中斷裂東段，地理位置位于青海省瑪多縣苦海以北（圖1）。斷裂在山前形成了明顯的斷層陡坎，并發(fā)現(xiàn)了大量水系左旋，斷層埡口、擠壓隆起等構(gòu)造地貌特征，由于該地區(qū)位于三江源保護(hù)區(qū)內(nèi)，受人為因素的影響小，這些地貌特征保存完整，因此采用無(wú)人機(jī)攝影測(cè)量技術(shù)，有針對(duì)性地進(jìn)行采集，快速獲取高精度的地形地貌數(shù)據(jù)。

圖1 區(qū)域構(gòu)造Fig. 1 Regional structure map

3 研究點(diǎn)數(shù)據(jù)提取

在前期無(wú)人機(jī)低空攝影和后期數(shù)據(jù)處理的基礎(chǔ)上，獲取了研究區(qū)域內(nèi)的正射影像。同時(shí)，基于高密度點(diǎn)云生成的DEM 數(shù)據(jù)分辨率較高，與斷層相關(guān)的微地貌特征較明顯。為清晰展示斷裂的活動(dòng)跡象，本文選取2 個(gè)典型地貌點(diǎn)進(jìn)行詳細(xì)解譯，并通過(guò)實(shí)地野外調(diào)查驗(yàn)證此次無(wú)人機(jī)數(shù)據(jù)采集結(jié)果。

本次飛行作業(yè)分別選取了龍通村北、乎勒塘村北作為本次工作的2 個(gè)試驗(yàn)點(diǎn)（以下簡(jiǎn)稱(chēng)為點(diǎn)1 和點(diǎn)2），飛行高度設(shè)定為80 m，分別獲取419，489 張照片，照片重疊率達(dá)70%以上。其中，點(diǎn)1 處共有227 262 個(gè)連接點(diǎn)，生成密集點(diǎn)云11 784 833 個(gè)點(diǎn)，生成的DEM 數(shù)據(jù)分辨率為6 985 pix×4 552 pix，覆蓋面積為0.549 km2；點(diǎn)2 處共有253 024 個(gè)連接點(diǎn)，生成密集點(diǎn)云12 423 887 個(gè)點(diǎn)，生成的DEM 數(shù)據(jù)分辨率為7 004 pix×4 233 pix，覆蓋面積為0.325 km2。

航高為航攝時(shí)飛機(jī)的飛行高度，根據(jù)起算基準(zhǔn)的不同可分為絕對(duì)航高與相對(duì)航高。相對(duì)航高為無(wú)人機(jī)在飛行時(shí)相機(jī)相對(duì)于某一基準(zhǔn)面的高度，是相對(duì)于作業(yè)區(qū)域內(nèi)地面平均高程基準(zhǔn)面的設(shè)計(jì)航高。一般來(lái)說(shuō)，在其他因素不變的情況下，航高越低（飛機(jī)飛行越低），地面分辨率越高（朱海斌等，2018）。

根據(jù)高建彩（2022）的規(guī)定，相對(duì)航高計(jì)算公式如下：

式中，H為相對(duì)航高；f為攝影鏡頭的焦距；G為影像的地面分辨率；a為像元尺寸。

通過(guò)換算可得到點(diǎn)1 的數(shù)據(jù)精度為2.07 cm/pix，點(diǎn)2 的數(shù)據(jù)精度為2.15 cm/pix，滿足斷裂位錯(cuò)量的計(jì)算要求。

傳統(tǒng)測(cè)量斷裂的垂直位移方法主要為差分GPS，該方法獲得的數(shù)據(jù)與DEM 數(shù)據(jù)相差較小，可滿足工作需要，而測(cè)量斷裂水平位移時(shí)主要使用皮尺等工具，受人為因素影響較大，野外地形相對(duì)復(fù)雜時(shí)測(cè)量數(shù)據(jù)無(wú)法真實(shí)反映斷層位錯(cuò)量。為排除個(gè)人主觀因素的影響，本文從高精度DEM 數(shù)據(jù)上直接提取陡坎高程數(shù)據(jù)，同時(shí)，剖面線的選取能夠反映斷層陡坎兩側(cè)的地形變化。

3.1 龍通村北試驗(yàn)點(diǎn)（點(diǎn)1）

點(diǎn)1 位于興?？h溫泉鄉(xiāng)西北約28 km 處，從遙感影像中可看出斷裂在此處的線性特征明顯（圖2），呈NWW 向展布，長(zhǎng)約700 m。斷層陡坎主要發(fā)育在山前沖洪積扇體后緣（圖3），水系有明顯的左旋位錯(cuò)量（劉超等，2021）。

圖2 點(diǎn)1 遙感影像Fig. 2 Remote sensing image of point one

圖3 點(diǎn)1 斷層陡坎Fig. 3 Fault scarp of point one

通過(guò)野外地質(zhì)調(diào)查發(fā)現(xiàn)，該處山前形成多期洪積扇，斷裂在洪積扇后緣形成明顯的斷層陡坎、水系扭錯(cuò)及地震形成的跌水。斷層陡坎連續(xù)性較好，多形成長(zhǎng)約50 m、高0.5～1.1 m 的斷層陡坎。在一處沖溝內(nèi)，形成長(zhǎng)約30 m、高2.5 m 的斷層跌水。

由于該地區(qū)海拔達(dá)4 500 m，人工測(cè)量微地貌難度較大，為更好地獲取該試驗(yàn)點(diǎn)微地貌特征，利用無(wú)人機(jī)航拍及數(shù)字化處理技術(shù)，對(duì)該試驗(yàn)點(diǎn)進(jìn)行高精度低空攝影，獲取大量無(wú)人機(jī)數(shù)據(jù)，通過(guò)后期處理，生成DEM/DOM 矢量圖（圖4（a））。分析DEM/DOM 矢量圖可知該處發(fā)育3 期洪積扇（圖4（b）），斷層橫切Fan2 和Fan3 洪積扇。同時(shí)，調(diào)查發(fā)現(xiàn)，F(xiàn)an2 洪積扇上的沖溝被左旋水平斷錯(cuò)，形成了斷頭溝。

基于生成的高分辨率DEM 數(shù)據(jù)，對(duì)3 條沖溝進(jìn)行了位移重塑，其中，東側(cè)斷層水平位移僅2.1 m（圖4（c）），并在沖溝后緣形成了高達(dá)2.5 m 的跌水；西側(cè)斷層水平位移可劃分為2 期，最大左旋位移為15.4 m，最小位移為4.85 m（圖4（d））。依據(jù)水平位移推測(cè)，認(rèn)為西側(cè)的沖溝也應(yīng)存在1 期2 m 左右的水平位移，該位錯(cuò)量可能由最新一次斷裂活動(dòng)所致。

圖4 點(diǎn)1 基于DEM 數(shù)據(jù)的水系位錯(cuò)恢復(fù)Fig. 4 Water system dislocation recovery based on DEM data of point one

昆中斷裂在研究段不僅存在水平位錯(cuò)量，同時(shí)還具備一定的逆沖分量。根據(jù)此次獲取的DEM 數(shù)據(jù)可知，在山前的沖洪積扇上，有1 條明顯的延伸線（圖5），其空間展布未受地形影響，通過(guò)實(shí)地調(diào)查，確認(rèn)該延伸線為斷裂在山前形成的斷層陡坎。

圖5 點(diǎn)1 DEM 提取剖面位置示意Fig. 5 Schematic diagram of DEM extraction section position of point one

利用本次生成的DEM 數(shù)據(jù)共獲取4 條陡坎高度（圖6～圖9），陡坎高度基本為0.6～0.9 m，且4 條剖面線在斷層陡坎兩側(cè)形態(tài)并無(wú)較大差異，原始地貌保存較好，筆者認(rèn)為這4 條陡坎高度基本一致，應(yīng)是同次斷層活動(dòng)形成，且陡坎高度可代表此次斷裂活動(dòng)的實(shí)際垂直位錯(cuò)量，而陡坎高度差異可能是斷層陡坎形成后受水系侵蝕作用導(dǎo)致的。

圖6 點(diǎn)1 AA′地形剖面Fig. 6 AA' topographic profile of point one

圖7 點(diǎn)1 BB′地形剖面Fig. 7 BB' topographic profile of point one

圖8 點(diǎn)1 CC′地形剖面Fig. 8 CC' topographic profile of point one

圖9 點(diǎn)1 DD′地形剖面Fig. 9 DD' topographic profile of point one

3.2 乎勒塘村北試驗(yàn)點(diǎn)（點(diǎn)2）

點(diǎn)2 位于溫泉鄉(xiāng)西北約33 km 處，從遙感影像中可看出斷裂在此處的線性特征較明顯（圖10），斷裂呈NWW 向展布，陡坎長(zhǎng)度約610 m，斷層陡坎主要發(fā)育在山前沖洪積扇后緣，水系有明顯的左旋位錯(cuò)量（圖11（a））。點(diǎn)2 處的斷層陡坎地貌特征明顯（圖11（b）），由圖12（a）和圖12（b）可知，山前至少發(fā)育了3 期洪積扇，斷層橫切Fan2 和Fan3 洪積扇。

圖10 點(diǎn)2 遙感影像Fig. 10 Remote sensing image of point two

圖11 點(diǎn)2 斷層地貌特征Fig. 11 Geomorphological characteristics of fault in point two

圖12 點(diǎn)2 基于DEM 數(shù)據(jù)的水系位錯(cuò)恢復(fù)Fig. 12 Water system dislocation recovery based on DEM data of point two

通過(guò)解譯無(wú)人機(jī)數(shù)據(jù)發(fā)現(xiàn)，F(xiàn)an2 洪積扇上的沖溝被左旋水平斷錯(cuò)，并形成了2 處斷頭溝。通過(guò)對(duì)沖溝位錯(cuò)的恢復(fù)，圖12（c）中的沖溝位錯(cuò)共發(fā)育2 期，2 期的左旋位錯(cuò)量分別為4.8 m 和12.1 m；圖12（d）中的沖溝左旋位錯(cuò)量為5.2 m；與點(diǎn)1 處的水系左旋位錯(cuò)量基本一致。

利用本次生成的DEM 數(shù)據(jù)共獲取4 條陡坎高度（圖13），陡坎高度基本為0.45～0.75 m（圖14～圖17），且3 條剖面線在斷層陡坎兩側(cè)形態(tài)并無(wú)較大差異，原始地貌保存較好，筆者認(rèn)為這3 條陡坎高度基本一致，應(yīng)為同次斷層活動(dòng)形成，但與點(diǎn)1 相比，該處的垂直位錯(cuò)量明顯較小，應(yīng)是由于此處水流順著陡坎方向匯集，使陡坎上側(cè)的地貌侵蝕更劇烈，因此陡坎高度較小，因此，該處陡坎高度不能代表此次斷裂活動(dòng)的實(shí)際垂直位錯(cuò)量，點(diǎn)1 點(diǎn)處的陡坎高度更能代表斷裂的實(shí)際位錯(cuò)量。

圖13 點(diǎn)2 DEM 提取剖面位置示意Fig. 13 Schematic diagram of DEM extraction section position

圖14 點(diǎn)2 AA′地形剖面Fig. 14 AA' topographic profile

圖15 點(diǎn)2 BB′地形剖面Fig. 15 BB' topographic profile

圖16 點(diǎn)2 CC′地形剖面Fig. 16 CC' topographic profile

圖17 點(diǎn)2 DD′地形剖面Fig. 17 DD' topographic profile

4 結(jié)論

本文采用大疆精靈4 RTK 無(wú)人機(jī)測(cè)繪影像進(jìn)行3 D 模型構(gòu)建，能夠較準(zhǔn)確地還原實(shí)際地貌場(chǎng)景，獲得較高分辨率的DEM 數(shù)據(jù)。通過(guò)對(duì)采集圖像進(jìn)行處理，可進(jìn)行斷層地貌變形的精確測(cè)量，為研究微地貌特征、斷裂幾何展布等提供方便有效的方法。

本文利用無(wú)人機(jī)攝影測(cè)量技術(shù)獲取了斷裂局部地區(qū)的微地貌數(shù)據(jù)?斷層陡坎及水系左旋的高精度DEM，對(duì)研究區(qū)域內(nèi)斷裂水平、垂直位錯(cuò)量進(jìn)行精細(xì)研究，確定了昆中斷裂（巴隆-溫泉段）活動(dòng)斷層陡坎位置及斷裂最新活動(dòng)量，得出以下結(jié)論：

（1）無(wú)人機(jī)航拍技術(shù)相較于傳統(tǒng)測(cè)量方式，具有低成本、易操作、獲取數(shù)據(jù)方便的特點(diǎn)。同時(shí)，利用PhotoScan 進(jìn)行數(shù)據(jù)處理時(shí)，數(shù)據(jù)處理速度快、誤差小，能夠滿足斷層陡坎精細(xì)化測(cè)量的數(shù)據(jù)要求，凸顯無(wú)人機(jī)航拍技術(shù)的應(yīng)用潛力。但該技術(shù)也受天氣、飛行高度及續(xù)航等因素限制，飛行數(shù)據(jù)易產(chǎn)生異常，不利于后續(xù)工作的開(kāi)展。

（2）利用無(wú)人機(jī)獲取的高精度DEM 數(shù)據(jù)，分別得到昆中斷裂（巴隆-溫泉段）水平位錯(cuò)量為2.1、4.85、15.4 m，位錯(cuò)量之間近似存在著倍數(shù)關(guān)系，而15.4 m 的位錯(cuò)量可能是由4 次斷層活動(dòng)形成的累計(jì)位錯(cuò)量。因此，在下一步工作中，還應(yīng)探明是否存在1 期8 m 左右的斷層累計(jì)活動(dòng)位錯(cuò)量。

（3）通過(guò)對(duì)比2 處試驗(yàn)點(diǎn)DEM 數(shù)據(jù)提取的8 條陡坎數(shù)據(jù)，陡坎兩側(cè)的地貌形態(tài)基本相似，無(wú)明顯差異，乎勒塘村北試驗(yàn)點(diǎn)處的陡坎高度受水流侵蝕較嚴(yán)重，而龍通村北試驗(yàn)點(diǎn)處的陡坎高度均為0.6～0.9 m，可初步認(rèn)為這是最新一次斷層運(yùn)動(dòng)形成的陡坎，可代表斷裂實(shí)際位錯(cuò)量。