摘要：由于三峽庫區(qū)地勢高陡、植被茂密，傳統(tǒng)人工調(diào)查難以識別危巖崩塌地質(zhì)災(zāi)害隱患。以2023年10月7日龔家坊危巖崩塌為例，現(xiàn)場采集崩塌前后機載激光雷達數(shù)據(jù)，過濾植被生成地面點云數(shù)據(jù)，建立了點云對比的高程差分模型，從而識別出危巖崩塌空間位置及變形量。結(jié)果表明：該高程差分模型可以清晰識別危巖崩塌空間位置及變形量，通過與光學(xué)三維影像、崩塌前后數(shù)字高程模型（DEM）對比，驗證了高程差分模型的可行性。研究成果可為三峽庫區(qū)地質(zhì)災(zāi)害隱患識別、地災(zāi)監(jiān)測提供重要依據(jù)。

關(guān) 鍵 詞：危巖變形識別；機載激光雷達；點云對比；高程差分模型；龔家坊危巖；三峽庫區(qū)

中圖法分類號：P642.2 文獻標(biāo)志碼：ADOI：10.16232/j.cnki.1001-4179.2024.S2.021

0 引言

三峽水庫是中國重要水利樞紐工程，庫區(qū)兩岸山高坡陡、褶皺斷層發(fā)育，在歷史上形成了許多崩塌地質(zhì)災(zāi)害，其中巫山縣巫峽段就有龔家坊危巖^［1-2］、望霞危巖^［3-4］、箭穿洞危巖^［5-6］等地質(zhì)災(zāi)害，威脅航道運行及周邊群眾生命財產(chǎn)安全。由于區(qū)域植被發(fā)育、危巖高陡，傳統(tǒng)人工地面調(diào)查很難識別危巖地質(zhì)災(zāi)害隱患。

目前，研究學(xué)者通過兩期或者多期點云數(shù)據(jù)對比捕捉地質(zhì)災(zāi)害變形信息，可以判斷地質(zhì)災(zāi)害空間位置及變形趨勢，是地質(zhì)災(zāi)害隱患識別的有效途徑。點云主要獲取方式包括無人機貼近攝影測量（Nap-of-the-object photogrammetry）、地面三維激光掃描技術(shù)（Terrestrial Laser Scanning，TLS）及無人機激光雷達技術(shù)（Light Detection and Ranging，LiDAR），3種技術(shù)手段原理不同，在地質(zhì)災(zāi)害識別應(yīng)用中有各自優(yōu)勢和局限。貼近攝影測量是通過貼近被測物體表面而不斷調(diào)整拍攝角度的攝影測量技術(shù)，可以獲取毫米級超高分辨率影像^［7］。一方面，高分辨率影像可以用于危巖結(jié)構(gòu)調(diào)查^［8］；另一方面，在貼近攝影測量三維模型的基礎(chǔ)上采樣獲取密集點云數(shù)據(jù)，通過兩期數(shù)據(jù)對比獲取毫米級變形，可以定量分析庫區(qū)消落帶岸坡變化情況^［9］。但貼近攝影測量無法穿透植被，一般適用于庫區(qū)裸露危巖及庫岸消落帶變形分析。地面三維激光掃描是通過發(fā)射激光獲取高精度和高密度的三維點云，常被應(yīng)用于邊坡及危巖調(diào)查和地質(zhì)災(zāi)害監(jiān)測中。董秀軍等^［10］通過擬合三維點云并提取結(jié)構(gòu)產(chǎn)狀信息，實現(xiàn)了高陡邊坡及危巖體的結(jié)構(gòu)面調(diào)查。薛強等^［11］利用多期三維激光掃描數(shù)據(jù)，生成多期數(shù)字高程模型（DEM）及點云剖面，可以分析滑坡局部變形特征。高梓成^［12］針對箭穿洞危巖體開展多期三維激光掃描工作，通過多期危巖體數(shù)據(jù)切片比較和整體比較得到了形變量，驗證了利用三維激光掃描技術(shù)進行實際崩滑監(jiān)測與調(diào)查方法的可行性。但實踐表明：三維激光掃描架設(shè)受觀測角度、觀測距離等因素的制約，在庫區(qū)高陡且狹小空間內(nèi)使用存在一定的局限性。機載LiDAR是一種主動式遙感測量手段，通過多次回波技術(shù)可以獲得高密集的點云（origin point），經(jīng)過濾波、分類可以過濾植被獲取真實地面點云（ground point），彌補了貼近攝影測量無法穿透植被和地貌三維激光掃描架設(shè)限制的不足。研究學(xué)者一般將地面點云生成高精度數(shù)值高程模型（DEM），利用紋理和陰影等差異可以人工目視解譯出地質(zhì)災(zāi)害要素，從而進行地質(zhì)災(zāi)害調(diào)查^［13-15］、隱患識別^［16-19］和風(fēng)險評價^［20-21］。在兩期或多期機載LiDAR點云對比方面，姚頑強等^［22］采集多時相機載LiDAR點云數(shù)據(jù)，提出了點云直接比較的算法（cloud to cloud，C2C），驗證了其在礦山開采沉陷監(jiān)測中的應(yīng)用。

鑒于機載雷達可以過濾植被獲取地面點云及點云對比的能力，本文以2023年10月7日三峽庫區(qū)龔家坊岸坡危巖崩塌為研究對象，利用崩塌前后機載LiDAR點云數(shù)據(jù)，建立了一種通過兩期點云數(shù)據(jù)識別危巖變形空間位置及變形量的技術(shù)方法，以期為三峽庫區(qū)地質(zhì)災(zāi)害隱患識別、地災(zāi)監(jiān)測提供重要依據(jù)。

1 研究區(qū)概況

三峽庫區(qū)龔家坊岸坡危巖崩塌位于望霞背斜北西翼、巫山向斜南東翼之間（圖1～2）。望霞背斜形態(tài)近似箱狀，軸部平緩開闊，往兩翼變陡；巫山向斜南西段核部緩，翼部陡，呈屜形，往北東漸次變狹窄，軸面傾向北，樞紐起伏，兩翼基本對稱，次級褶皺發(fā)育。由此，區(qū)內(nèi)地質(zhì)構(gòu)造條件復(fù)雜，為崩塌地質(zhì)孕育提供了有利條件。

2 崩塌結(jié)構(gòu)特征

根據(jù)現(xiàn)場地質(zhì)調(diào)查，危巖體地層巖性主要為三疊系下統(tǒng)嘉陵江組四段（T₁j⁴）巖溶角礫巖及三疊系下統(tǒng)嘉陵江組三段（T₁j³）中—厚層狀白云質(zhì)灰?guī)r，巖體溶蝕明顯、結(jié)構(gòu)裂隙、坡體地表裂縫發(fā)育（圖3）。采用赤平投影分析可知崩塌塊體由坡面、層面、多個裂隙面相互切割形成，結(jié)構(gòu)破碎。其中，巖層產(chǎn)狀342°∠44°，坡向180°，結(jié)構(gòu)類型為反向坡。結(jié)構(gòu)裂隙主要發(fā)育有3組，危巖體的主控裂隙（LX1）優(yōu)勢產(chǎn)狀175°∠45°，張開度2～5 mm，含泥質(zhì)充填物，裂隙面平直，延伸長3～10 m，為危巖體的主控裂隙；第2組裂隙（LX2）優(yōu)勢產(chǎn)狀80°∠65°，裂隙張開度1～15 mm，含泥質(zhì)充填物，裂隙面鋸齒狀，延伸長0.5～3.0 m；第3組裂隙（LX3）優(yōu)勢產(chǎn)狀195°∠59°，裂隙張開度2～13 mm，局部少量泥質(zhì)充填。

3 高程差分模型構(gòu)建方法

采用華測大黃蜂BB4無人機搭載AS1300HL激光雷達系統(tǒng)采集激光雷達數(shù)據(jù)，采集時間分別是2022年5月10日和2024年3月30日，平均地面點云密度達到29.8點/m²。高程差分模型構(gòu)建方法如下（圖4）：首先，采用CoPre軟件將采集數(shù)據(jù)、華測云基站進行點云預(yù)處理，得到原始點云數(shù)據(jù)；隨后，采用CoProcess 2.6進行濾波、分類，過濾掉建筑物點云和植被點云，得到地面點云數(shù)據(jù)；最后，將兩期點云數(shù)據(jù)導(dǎo)入Polyworks，將第1期（崩塌前）地面點云數(shù)據(jù)生成三角化網(wǎng)格模型，并作為參考對象，將第2期（崩塌后）地面點云數(shù)據(jù)作為測量對象，按照高程方向進行差分，進而獲得高程差分模型。

4 高程差分模型分析與驗證

通過崩塌前后機載LiDAR點云數(shù)據(jù)對比，建立了危巖的高程差分模型。模型清晰顯示危巖崩塌空間位置及崩塌后坡體削方治理范圍（圖5（a））。其中，藍色區(qū)域代表減少，紅色區(qū)域代表增加。藍色最深區(qū)域表示此次位移崩塌位置，崩塌前后高程變形量達到13.3 m；藍色相對較淺的區(qū)域表示崩塌后坡體削方范圍，削方深度0.9～1.5 m。將高程差分模型與光學(xué)影像三維模型（圖5（b））、崩塌前后機載雷達DEM（圖5（c）、（d））進行對比驗證，結(jié)果表明：危巖崩塌范圍及變形量基本一致，驗證了本文高程差分模型構(gòu)建方法在庫區(qū)危巖崩塌變形識別中的可行性。

5 結(jié)論

（1）本文高程差分模型識別了危巖崩塌變形空間位置及變形量，顯示出“10·7”龔家坊崩塌變形量達到13.3 m，與現(xiàn)場結(jié)果基本一致，驗證了本文高程差分模型構(gòu)建方法的可行性。

（2）建議定期開展庫區(qū)大面域機載雷達點云數(shù)據(jù)變形對比，提高變形識別精度，這是庫區(qū)地質(zhì)災(zāi)害隱患識別的重要研究方向。

參考文獻：

［1］ GU D，HUANG D.A complex rock topple-rock slide failure of an anaclinal rock slope in the Wu Gorge，Yangtze River，China［J］.Engineering Geology，2016，208：165-180.

［2］ HUANG B，YIN Y，LIU G，et al.Analysis of waves generated by Gongjiafang landslide in Wu Gorge，Three Gorges Reservoir，on November 23，2008［J］.Landslides，2012，9（3）：395-405.

［3］ 高買燕，唐紅梅，曾云松，等.崩塌災(zāi)害危險性評價方法及應(yīng)用［J］.重慶交通大學(xué)學(xué)報（自然科學(xué)版），2013，32（3）：446-450，537.

［4］ 黃達，楊偉東，陳智強.考慮軟弱基座風(fēng)化效應(yīng)的望霞危巖崩塌機制分析［J］.人民長江，2018，49（6）：64-69.

［5］ 蔣文明，王魯琦，趙鵬，等.三峽庫區(qū)箭穿洞危巖體變形破壞模式與防治效果分析［J］.中國地質(zhì)災(zāi)害與防治學(xué)報，2021，32（5）：105-112.

［6］ 趙志男，李濱，賀凱，等.西南山區(qū)柱狀危巖體基座壓裂潰屈斷裂力學(xué)解析［J］.巖石力學(xué)與工程學(xué)報，2022，41（增2）：3132-3140.

［7］ 言司.獨辟蹊徑，不斷創(chuàng)新貼近攝影測量：第三種攝影測量方式的誕生：專訪武漢大學(xué)遙感信息工程學(xué)院張祖勛院士［J］.中國測繪，2019（10）：31-37.

［8］ 姚富潭，吳明堂，董秀軍，等.基于貼近攝影測量技術(shù)的高陡危巖體結(jié)構(gòu)面調(diào)查方法［J］.成都理工大學(xué)學(xué)報（自然科學(xué)版），2023，50（2）：218-228.

［9］ 王勛，蔣欽偉，姜正容，等.一種三峽庫區(qū)消落帶岸坡變化檢測方法［J］.遙感信息，2022，37（4）：21-28.

［10］董秀軍，黃潤秋.三維激光掃描技術(shù)在高陡邊坡地質(zhì)調(diào)查中的應(yīng)用［J］.巖石力學(xué)與工程學(xué)報，2006，25（增2）：3629-3635.

［11］薛強，畢俊擘，李政國，等.三維激光掃描技術(shù)在閻家溝滑坡變形監(jiān)測中的應(yīng)用［J］.中國地質(zhì)災(zāi)害與防治學(xué)報，2018，29（3）：108-112.

［12］高梓成.基于TLS的三峽庫區(qū)典型試點崩滑災(zāi)害監(jiān)測調(diào)查研究［D］.北京：北京建筑大學(xué)，2020.

［13］SHE J，MABI A，LIU Z，et al.Analysis using high-precision airborne LiDAR data to survey potential collapse geological hazards［J］.Advances in Civil Engineering，2021（24）：6475942.

［14］杜磊，陳潔，李敏敏，等.機載激光雷達技術(shù)在滑坡調(diào)查中的應(yīng)用：以三峽庫區(qū)張家灣滑坡為例［J］.國土資源遙感，2019，31（1）：180-186.

［15］劉圣偉，郭大海，陳偉濤，等.機載激光雷達技術(shù)在長江三峽工程庫區(qū)滑坡災(zāi)害調(diào)查和監(jiān)測中的應(yīng)用研究［J］.中國地質(zhì)，2012，39（2）：507-517.

［16］董遠峰，李偉，鄧玖林，等.松坎河流域地質(zhì)災(zāi)害識別與發(fā)育分布規(guī)律研究［J］.人民長江，2023，54（增1）：54-57.

［17］彭藝偉，董琦，田沖，等.基于機載激光雷達的地質(zhì)災(zāi)害識別關(guān)鍵技術(shù)及應(yīng)用研究［J］.安全與環(huán)境工程，2021，28（6）：100-108.

［18］郭晨，許強，董秀軍，等.基于SVF地形可視化方法的滑坡識別：以四川省丹巴縣城典型滑坡為例［J］.成都理工大學(xué)學(xué)報（自然科學(xué)版），2021，48（6）：705-713.

［19］GORUM T.Landslide recognition and mapping in amixed forest environment from airborne LiDAR data［J］.Engineering Geology，2019，258（14）：105155.

［20］甄俊偉，黃智煒，章桂芳，等.基于LiDAR數(shù)據(jù)的汕尾火山嶂地質(zhì)災(zāi)害風(fēng)險評價［J］.中山大學(xué)學(xué)報（自然科學(xué)版）（中英文），2024，63（1）：10-23.

［21］ABDULWAHID WM，PRADHAN B.Landslide vulnerability and risk assessment for multi-hazard scenarios using airborne laser scanning data（LiDAR）［J］.Landslides，2016（14）：1057-1076.

［22］姚頑強，高康洲，鄭俊良，等.基于機載LiDAR點云C2C算法的礦山沉陷監(jiān)測研究［J］.煤炭工程，2023，55（4）：162-167.

（編輯：劉媛）