摘要：

危巖墜落受多種因素影響，具有不確定性，研究危巖的潛在受災(zāi)區(qū)域并區(qū)分不同影響程度能夠為災(zāi)害防治提供數(shù)據(jù)支撐。以廣西南寧市內(nèi)莊屯危巖為例，基于傾斜與貼近攝影融合建模技術(shù)構(gòu)建三維模型，精準(zhǔn)獲取巖體結(jié)構(gòu)信息及危巖體特征信息；并通過三維數(shù)值模擬對危巖的運動路徑、運動過程、實時動能及彈跳高度進行研究，根據(jù)多次模擬結(jié)果提出了將潛在受災(zāi)區(qū)域進行影響區(qū)劃分的方法。研究表明：融合建模有效提高了陡立山體細(xì)節(jié)特征的可視性，更利于危巖運動模擬參數(shù)的精確提?。换诙啻文M結(jié)果分析，可以確定潛在受災(zāi)區(qū)域；影響區(qū)劃分方法能夠進一步將潛在受災(zāi)區(qū)進行簡易、快速的區(qū)域劃分。研究結(jié)果對高陡危巖體的調(diào)查具有參考意義。

關(guān)" 鍵" 詞：

高位危巖； 融合建模； 傾斜攝影； 貼近攝影； 數(shù)值模擬； 落石分區(qū)

中圖法分類號： P231;P642.2

文獻(xiàn)標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2024.07.018

收稿日期：

2023-12-29

；接受日期：

2024-04-19

基金項目：

國家自然科學(xué)基金項目（42072306）

作者簡介：

劉兵兵，男，碩士研究生，主要從事攝影測量與工程地質(zhì)方面的研究。E-mail：1749910575@qq.com

通信作者：

吳" 福，男，正高級工程師，碩士，主要從事地質(zhì)環(huán)境監(jiān)測、地質(zhì)災(zāi)害防治方面的研究。E-mail：1404630357@qq.com

Editorial Office of Yangtze River. This is an open access article under the CC BY-NC-ND 4.0 license.

文章編號：1001-4179（2024） 07-0138-07

引用本文：

劉兵兵，吳福，劉昶，等.

基于無人機融合建模的危巖運動模擬及影響區(qū)劃分

［J］.人民長江，2024，55（7）：138-144，159.

0" 引 言

斜坡上已有變形跡象、具有失穩(wěn)破壞征兆的巖塊體，稱為危巖［1］。近年來，中國發(fā)生了多起危巖失穩(wěn)導(dǎo)致人員傷亡、公路河道堵塞、堤壩損毀等事件，如2018年7月，湖南省張家界大峽谷景區(qū)發(fā)生巖體滾落，使4名游客受傷，2名游客死亡；2021年3月，四川省樂山市峨邊縣發(fā)生巨石滾落事故，導(dǎo)致村民3人受傷，1人被埋；2022年7月27日，三峽庫區(qū)奉節(jié)縣瞿塘峽吊嘴危巖發(fā)生崩塌墜入長江，濺起水花高約25 m，所幸無過往船只被擊中。高位危巖具有地面高差大、突發(fā)性強等特點［2］，通過人工難以接近進行精細(xì)化調(diào)查，如何查明危巖體特征信息并預(yù)測失穩(wěn)后的受災(zāi)區(qū)域一直是危巖調(diào)查研究中的重點。

隨著無人機技術(shù)的迅猛發(fā)展，國內(nèi)外采用無人機攝影測量進行三維模型構(gòu)建，基于模型對危巖開展調(diào)查的方式被眾多學(xué)者所采用［3］，如Yan等［4］利用無人機攝影測量開發(fā)了一種新的半自動危險巖塊的提取方法，能夠自動獲取識別塊體的位置、幾何參數(shù)和穩(wěn)定性結(jié)果；鄭紅霞等［5］基于貼近攝影測量技術(shù)提出了一種自動識別巖層剖面特征點進行巖層面重構(gòu)的方法；Vasuki等［6］利用無人機攝影測量技術(shù)采集地質(zhì)結(jié)構(gòu)體上的影像信息，通過采集的影像信息，基于圖像處理技術(shù)，半自動化提取出巖體表面結(jié)構(gòu)面的產(chǎn)狀。以往多數(shù)研究者采用傾斜攝影測量技術(shù)進行數(shù)據(jù)獲取，而貼近攝影測量技術(shù)提出之后，也有學(xué)者將貼近攝影測量技術(shù)應(yīng)用到巖體特征信息的調(diào)查中來，如孫琦等［7］基于貼近攝影技術(shù)研究了采樣尺寸、采樣間距對結(jié)構(gòu)面粗糙度評價結(jié)果的影響。經(jīng)實踐表明，無論是傾斜攝影或貼近攝影技術(shù)均存在一定的缺陷［8］，如傾斜攝影在陡立面山體容易出現(xiàn)拍攝盲區(qū)，導(dǎo)致模型出現(xiàn)空洞拉花等不良現(xiàn)象［9］；貼近攝影在植被發(fā)育區(qū)近距離飛行時易造成無人機墜機現(xiàn)象［10］。將兩種技術(shù)結(jié)合能夠彌補單一技術(shù)采集的不足，但利用傾斜與貼近技術(shù)進行融合建模的研究尚少。

基于三維模型的落石運動特征模擬結(jié)果更加準(zhǔn)確直觀。對于落石運動特征的研究，國內(nèi)學(xué)者常采用數(shù)值模擬的方法，其中二維模擬軟件有Rockfall、二維離散元UDEC等，三維軟件有Rockyfor3D、Unity3D等。王豪等［11］基于Unity3D平臺研發(fā)了可以輸出運動軌跡、速度、能量及彈跳高度的崩塌運動過程模擬軟件；呂權(quán)儒等［12］利用Rockfall軟件模擬崩落體失穩(wěn)后的滾落運動特征，并劃分出不同級別的崩塌影響區(qū)，為最終的防治提供依據(jù)。但以上多數(shù)研究止步于確定出危巖潛在受災(zāi)區(qū)域。為了能夠給災(zāi)害治理提供更加詳細(xì)的數(shù)據(jù)支持，對潛在受災(zāi)區(qū)域進行不同影響程度的區(qū)域劃分方法還需進一步研究。

本文以南寧市馬山縣內(nèi)莊屯房屋后陡立山體危巖為研究對象，通過傾斜與貼近攝影融合建模技術(shù)獲取了該地的精細(xì)化三維模型，基于三維模型對巖體結(jié)構(gòu)面特征及危巖體特征信息進行提取，利用崩塌模擬軟件MSRA-3D對其動力學(xué)進行分析，并根據(jù)多次模擬結(jié)果將危巖失穩(wěn)后可能造成的受災(zāi)范圍進行不同影響程度的劃分，以為后期危巖防治工程設(shè)計提供可靠的數(shù)據(jù)支撐。

1" 研究區(qū)概況

研究區(qū)位于南寧市馬山縣福蘭村內(nèi)莊屯房屋后山體（圖1），該地屬于峰叢洼地地貌區(qū)，峰叢基座相連，峰頂多呈錐狀，相對高差約300 m，山體中下部坡度為40°～50°，中上部為80°～90°，危巖多分布于山體陡立面處，調(diào)查人員難以靠近。出露地層巖性為泥盆系中統(tǒng)東崗嶺組（D2d）灰?guī)r和第四系溶余堆積（Qedl）。泥盆系中統(tǒng)東崗嶺組灰?guī)r呈灰色，中—厚層狀，雨量充沛，風(fēng)化溶蝕作用明顯，溶漕、溶隙、溶溝發(fā)育；第四系溶余堆積黏土分布于洼地及周邊，棕黃、灰褐色，厚度0～5 m，研究區(qū)南側(cè)200 m左右發(fā)育有一不明性質(zhì)斷層，長約25 km，走向北西。

2" 傾斜與貼近攝影技術(shù)融合建模

2.1" 原理概述

傾斜攝影測量技術(shù)是通過在同一攝影平臺上搭載多臺傳感器，從5個不同的角度進行數(shù)據(jù)采集，快速、高效地獲取豐富的數(shù)據(jù)信息［13-14］，但由于拍攝時飛行高度較高，在拍攝有巖體突出的陡崖數(shù)據(jù)時，在凹陷處容易出現(xiàn)拍攝盲區(qū)（圖2），從而在模型生產(chǎn)時會出現(xiàn)空洞、拉花、粘連等現(xiàn)象。而貼近攝影測量技術(shù)是正對物體表面進行拍攝［15］，能夠有效彌補傾斜攝影拍攝時的盲區(qū)空白，且其飛行時距離物體較近，能夠獲得更加清晰的影像數(shù)據(jù)。將貼近攝影對傾斜攝影進行拍攝補充，不僅在建立模型時能防止空洞等缺陷的出現(xiàn)，而且能提高重點區(qū)域的精細(xì)度（圖3），后續(xù)基于模型的信息提取結(jié)果也會更加完整、準(zhǔn)確。

2.2" 數(shù)據(jù)采集及模型建立

2.2.1" 影像數(shù)據(jù)采集

研究區(qū)數(shù)據(jù)采集選取大疆精靈4 RTK無人機，采集工作分為兩步：

（1） 對研究區(qū)整體范圍開展傾斜攝影，將傾斜攝影所獲得的高清影像照片利用大疆智圖軟件進行快速處理，生成粗模。此時的模型在陡崖處由于拍攝盲區(qū)可能會出現(xiàn)一些拉花、空洞等缺陷。該步驟的意義在于完成非重點區(qū)域的航測工作，找到模型缺陷位置，以及為貼近攝影提供基礎(chǔ)地形。

（2） 根據(jù)粗模確定基巖出露及有模型缺陷的重點區(qū)域，然后規(guī)劃貼近攝影飛行航線并開展正對坡面的攝影工作，貼近攝影的方式分為自動飛行和手動飛行，自動飛行是根據(jù)粗模及航線規(guī)劃軟件進行航線繪制，該種方法的優(yōu)點在于能夠保證無人機飛行期間的重疊度要求，避免出現(xiàn)漏拍現(xiàn)象，缺點是需要避免直立面雜生的植被影響，注意無人機的飛行安全。手動飛行則是飛手全程控制無人機飛行并定點拍攝，十分考驗飛手的專業(yè)能力，但能夠根據(jù)現(xiàn)場條件對飛行軌跡及時做出調(diào)整。本次貼近攝影由于重點區(qū)域數(shù)量多、面積小，采用手動飛行方式，兩次飛行共計獲取高清影像2 633張。

2.2.2" 三維模型構(gòu)建

導(dǎo)出無人機獲取的高清影像數(shù)據(jù)，通過影像數(shù)據(jù)內(nèi)置POS信息及高清影像信息，利用ContextCapture軟件進行空三解算、融合拼接等關(guān)鍵步驟，得到研究區(qū)的融合三維模型。圖4展示了該模型陡立面易出現(xiàn)拍攝盲區(qū)，未出現(xiàn)拉花、空洞等明顯模型缺陷，符合危巖信息提取的要求。

3" 巖體特征調(diào)查及穩(wěn)定性評價

3.1" 巖體結(jié)構(gòu)面特征識別

3.1.1" 結(jié)構(gòu)面產(chǎn)狀計算原理

結(jié)構(gòu)面特征包括結(jié)構(gòu)面的組數(shù)、產(chǎn)狀、跡長、間距及粗糙度等。結(jié)構(gòu)面的分布、產(chǎn)狀信息等參數(shù)能夠體現(xiàn)出巖體的裂隙特征，分析結(jié)構(gòu)面的分布及組合關(guān)系能夠了解巖體的斷裂模式、滑動模式等變形特征，從而評估巖體可能發(fā)生的破壞形式和規(guī)模。對于結(jié)構(gòu)面產(chǎn)狀信息提取，傳統(tǒng)方法是現(xiàn)場使用羅盤等工具進行直接測量，但對于高陡危巖體結(jié)構(gòu)面產(chǎn)狀信息由于地質(zhì)工作人員難以接近，無法精準(zhǔn)測量［16］，在實景三維模型上直接進行結(jié)構(gòu)面產(chǎn)狀信息的測量及統(tǒng)計，能夠有效解決人上不去的難題，為解決高陡地區(qū)對危巖穩(wěn)定性研究的阻礙提供新的思路和方法?；谌S模型的結(jié)構(gòu)面產(chǎn)狀信息量測，是根據(jù)三維模型上的每一點都具有坐標(biāo)信息，在平面上選取不在同一條直線上的3個特征點，來進行擬合平面，從而確定傾向傾角，簡稱為“三點法”，其原理如下：

首先設(shè)結(jié)構(gòu)面擬合平面方程式為

Z=AX+BY+C

（1）

式中：A、B、C為擬合平面方程參數(shù)，且A、B不同時為0，該平面法向量為n=［A，B，1］

然后采用最小二乘法解算得出A、B、C的值：

ABC=X1Y11X2Y21X3Y31-1Z1Z2Z3

（2）

根據(jù)結(jié)構(gòu)面產(chǎn)狀與平面方程參數(shù)的對應(yīng)關(guān)系可以求得傾角α和傾向β：

α=arctan（" A2+B2）

（3）

β0=arctanBA

（4）

當(dāng)A＜0時，

B≤0，β=β0

B＞0，β=β0+π

（5）

當(dāng)A＞0時，

β=β0+π

（6）

由此即可求出傾角α和傾向β［17］。

3.1.2" 巖體結(jié)構(gòu)面參數(shù)及精度驗證

研究區(qū)以灰?guī)r出露為主，溶蝕風(fēng)化強烈，巖層產(chǎn)狀為95°∠10°，主要發(fā)育兩組陡傾結(jié)構(gòu)面，第一組結(jié)構(gòu)面產(chǎn)狀為28°～45°或180°～210°∠79°～88°，該組結(jié)構(gòu)面走向與陡崖延伸方向近于平行，表面平直粗糙，貫通性好，發(fā)育較廣泛；第二組結(jié)構(gòu)面產(chǎn)狀為90°～110°或270°～300°∠80°～89°，該組結(jié)構(gòu)面表面平直光滑，貫通性較好。

為驗證模型量測產(chǎn)狀信息的準(zhǔn)確性，筆者在數(shù)據(jù)采集期間從研究區(qū)易于到達(dá)處選取4個點位分別按照傳統(tǒng)方法使用羅盤進行產(chǎn)狀信息量測，將羅盤量測結(jié)果與模型量測結(jié)果進行比較，模型測量結(jié)果與羅盤測量結(jié)果傾向傾角測量誤差均在3°以內(nèi)（表1），滿足工程精度要求。

3.2" 危巖識別及失穩(wěn)模式分析

基于三維模型，可以從各個角度直觀地對研究區(qū)域進行目視解譯，優(yōu)質(zhì)的三維模型，使得巖體的細(xì)節(jié)能夠更清晰地被觀察和發(fā)現(xiàn)，有助于發(fā)現(xiàn)更小的裂縫、節(jié)理等特征，通過分析巖層的形態(tài)、結(jié)構(gòu)、傾向、傾角等特征，并結(jié)合巖體臨空面、裂隙發(fā)育等條件能夠快速識別危巖，危巖編號及位置見圖5。

危巖體的失穩(wěn)模式即危巖的破壞模式，是其穩(wěn)定性研究及治理措施研究的基礎(chǔ)［18］。由于分類角度、調(diào)查區(qū)域不同，危巖體失穩(wěn)模式存在多種分類方案。陳洪凱等［1］根據(jù)危巖失穩(wěn)破壞的可能模式，將危巖分為墜落式、滑移式、傾倒式3種（圖6），其中墜落式危巖通常具有良好的臨空面，處于懸空或懸挑式的巖體由于發(fā)生拉斷、切斷而發(fā)生墜落；滑移式危巖通常在巖體與母巖之間發(fā)育有軟弱結(jié)構(gòu)面，危巖沿軟弱面滑移，于陡崖（坡）處塌落；傾倒式危巖多發(fā)育在黃土、直立或陡傾坡內(nèi)的巖層處，后緣發(fā)育裂隙逐漸擴張，危巖繞轉(zhuǎn)動點轉(zhuǎn)動發(fā)生傾倒塌落。

隨著三維模型在地質(zhì)災(zāi)害調(diào)查方面的應(yīng)用發(fā)展，利用EarthSurvey軟件，對危巖高度、空間尺寸等參數(shù)進行快速量測，依據(jù)危巖體裂隙及受力狀態(tài)，判斷危巖破壞特征及失穩(wěn)模式，提取參數(shù)及失穩(wěn)類型見表2。

3.3" 巖體穩(wěn)定性分析

3.3.1" 穩(wěn)定性分析方法

赤平投影能夠有效地表示空間結(jié)構(gòu)面之間的關(guān)系［19-20］，因此可以方便地用于邊坡結(jié)構(gòu)體的分析。大圓分析法是將各結(jié)構(gòu)面產(chǎn)狀以大圓表示的分析方法。該方法根據(jù)每組結(jié)構(gòu)面在赤平投影平面上的投影大圓或每兩組結(jié)構(gòu)面大圓組合形成的交線，與坡面在該投影平面上投影大圓和以及摩擦圓的相對位置來對其穩(wěn)定性做出判斷［21］。若代表結(jié)構(gòu)面的大圓或兩組結(jié)構(gòu)面大圓的交線小于邊坡面的傾角且都在摩擦圓中時，如圖7中結(jié)構(gòu)面1和結(jié)構(gòu)面2的交線，則認(rèn)為該結(jié)構(gòu)是潛在不穩(wěn)定的。

3.3.2" 危巖穩(wěn)定性分析

本文通過三維模型拾取坐標(biāo)，采用三點擬合平面法量測得WY1兩組主要控制節(jié)理產(chǎn)狀信息及坡面產(chǎn)狀信息，其中節(jié)理產(chǎn)狀為J1：189°∠81°，J2：283°∠83°，坡面產(chǎn)狀為192°∠83°。并根據(jù)南寧市相關(guān)地質(zhì)勘察資料，確定該處巖體內(nèi)摩擦角φ取值為30°，在赤平投影圖中進行結(jié)構(gòu)面組合分析，結(jié)果如圖8所示，J1、J2兩組節(jié)理的交點傾角小于坡面傾角且在摩擦圓中，易發(fā)生滑動，處于潛在不穩(wěn)定狀態(tài)。

4" 危巖運動特征模擬及風(fēng)險區(qū)域劃分

為研究危巖失穩(wěn)后的運動過程及影響，本文以WY1為例，采用崩塌運動全過程模擬軟件MSRA-3D對其動力學(xué)進行分析研究，該軟件能夠三維展示巖體失穩(wěn)后的運動路徑，輸出運動過程中巖體在不同位置上的速度、動能、彈跳高度等信息，能夠為災(zāi)害防治工程設(shè)計提供支持［12］。

4.1" 模型建立及參數(shù)設(shè)置

模擬所采用的坡體數(shù)據(jù)為前文通過無人機技術(shù)融合建模所得到的三維模型，通過數(shù)據(jù)轉(zhuǎn)換輸出FBX格式的模型，可直接導(dǎo)入使用。坡體表面特征參數(shù)參考廣西壯族自治區(qū)地方標(biāo)準(zhǔn)DB 45/T 1696-2018《危巖防治工程技術(shù)規(guī)范》，并通過已有落石情況進行反演測試，多次微調(diào)系數(shù)直至模擬結(jié)果與已有落石情況相符，有效保證后續(xù)模擬結(jié)果的準(zhǔn)確性。最終動摩擦系數(shù)、靜摩擦系數(shù)、恢復(fù)系數(shù)取值見表3。固定關(guān)節(jié)斷裂力閾值參考以往工程相同巖性實驗數(shù)據(jù)確定。危巖體數(shù)據(jù)則是根據(jù)前文所提取的巖體結(jié)構(gòu)面特征、危巖體特征等信息來確定危巖體模型的形狀、尺寸、空間位置、空間姿態(tài)。

4.2" 落石運動特征

在整個運動過程中，受結(jié)構(gòu)面切割的危巖脫離母巖后，首先在重力作用下發(fā)生墜落運動，再與坡面發(fā)生多次碰撞彈跳，最后進行滾落運動，運動軌跡及速度如圖9所示（此處僅展示W(wǎng)Y1的模擬結(jié)果）。輸出運動過程中的運動路徑、實時動能、彈跳高度等計算結(jié)果，如圖10所示（此處僅展示W(wǎng)Y1的計算結(jié)果）。根據(jù)模擬結(jié)果分析可知，該危巖在極端天氣下一旦發(fā)生墜落，若不加以防護，有擊中房屋的危險，對居民安全造成嚴(yán)重威脅。

4.3" 落石風(fēng)險區(qū)劃

通過多次模擬結(jié)果可以確定落石可能造成的災(zāi)害范圍，本文提出了一種將災(zāi)害范圍進行不同影響區(qū)域劃分的方法，分區(qū)結(jié)果在后期災(zāi)害防治等多個方面具有參考價值。

4.3.1" 影響區(qū)域劃分方法

不同影響區(qū)域劃分需要有明確的劃分規(guī)則，用來確定其比例或范圍，但如果根據(jù)某一坡形進行實驗分析，則缺乏了劃分規(guī)則的普適性。為此，本文提出用理

想化模型去研究落石?？课恢玫姆植家?guī)律，以求能從中獲取啟發(fā)。所謂理想化模型，即指巖體墜落區(qū)域是一片平地，無高大樹木或群集房屋等障礙物，不影響落石的運動軌跡，由于巖體的傾向、大小、高程等信息在多次實驗中不發(fā)生變化，且受上述條件控制，巖體墜落后一定有一個最佳理論位置，如圖11（a）中藍(lán)色的點。當(dāng)實驗的次數(shù)達(dá)到一定數(shù)量時，根據(jù)隨機試驗的原則，其分布則會呈現(xiàn)出距離理論位置越近越密集，越遠(yuǎn)越稀疏的規(guī)律。將該分布以垂直巖體主要滾落方向為橫坐標(biāo)，密集程度為縱坐標(biāo)進行統(tǒng)計，如圖11（b）所示。發(fā)現(xiàn)其分布近似于正態(tài)分布曲線，而正態(tài)分布可以根據(jù)標(biāo)準(zhǔn)差將其整體劃分為不同區(qū)域，標(biāo)準(zhǔn)差在正態(tài)分布中表征數(shù)據(jù)的離散程度，將此分類標(biāo)準(zhǔn)應(yīng)用于密集程度統(tǒng)計中，該點即表示為密集程度變化較大的分界位置，也正是分區(qū)較為符合的一個標(biāo)準(zhǔn)。因此將正態(tài)分布劃分比例應(yīng)用于樣本點劃分中，如圖11（c）所示，一個標(biāo)準(zhǔn)差之內(nèi)的區(qū)域占比約68％，兩個標(biāo)準(zhǔn)差之內(nèi)的區(qū)域占比約96％，兩個標(biāo)準(zhǔn)差之外的占比約4％，按照該區(qū)域比例與樣本點數(shù)目相結(jié)合，確定分區(qū)界限控制點，如圖11（a）所示，按此方法將影響區(qū)域劃分為重、中、輕度3個等級。

4.3.2" 影響區(qū)域劃分結(jié)果

危巖崩塌受多種因素影響，為了更加符合真實情況，MSRA-3D軟件在計算時引入了一定的隨機參數(shù)來減小確定參數(shù)數(shù)值模擬所造成的局限性，實驗過程中會隨機小幅度改變隨機參數(shù)，用來代替實際中遇到的降雨、地震等其他不確定因素影響。本文以WY1為例，通過多次模擬，得到100處危巖最終?？课恢脭?shù)據(jù)，選取處于最外圍4個位置作為可能造成的災(zāi)害范圍邊界，在垂直巖體主要滾落方向上，以最邊界的點為原點，距原點距離為左縱軸進行數(shù)據(jù)統(tǒng)計，且按此距離進行排序編號，并計算相鄰點位的差值為右縱軸。為方便制圖觀察，右縱軸設(shè)置向下為正方向，差值越小，則表明點位越密集，統(tǒng)計結(jié)果如圖12所示。相鄰點差值曲線呈現(xiàn)整體上先減小再增大的趨勢，說明點位分布特征為先稀疏再密集又稀疏，該特征證明采用正態(tài)分布規(guī)律劃分區(qū)域的方法是可行的。

將統(tǒng)計結(jié)果按照分布比例確定為不同影響區(qū)域，不同危巖體之間影響區(qū)的相互重疊區(qū)域采用風(fēng)險分級矩陣（圖13）確定影響等級；如圖14所示，研究區(qū)危巖落石共造成3處影響區(qū)域，其中WY1～4墜落位置形成區(qū)域一，WY5和WY6墜落位置分別形成區(qū)域二和區(qū)域三。本次分區(qū)結(jié)果能夠為危巖的治理與防范提供重要參考依據(jù)。

5" 結(jié) 論

（1）" 采用傾斜與貼近攝影技術(shù)融合構(gòu)建三維模型的方法，能夠有效避免傾斜攝影由拍攝盲區(qū)導(dǎo)致的模型空洞、拉花等現(xiàn)象，提高了模型的完整性，使信息提取結(jié)果更加準(zhǔn)確、真實。

（2）" 基于危巖運動模擬結(jié)果分析，房屋有被擊中的風(fēng)險，需要采取措施進行防護。

（3）" 首次提出采用正態(tài)分布規(guī)律將可能造成的災(zāi)害范圍進行不同影響程度的區(qū)域劃分方法，并在該示范區(qū)進行應(yīng)用，將落石范圍分為了重、中、輕3種不同程度的區(qū)域。

參考文獻(xiàn)：

［1］" 陳洪凱，唐紅梅，王蓉.三峽庫區(qū)危巖穩(wěn)定性計算方法及應(yīng)用［J］.巖石力學(xué)與工程學(xué)報，2004（4）：614-619.

［2］" 楊緒波，程強，袁進科，等.基于無人機傾斜攝影的危巖識別及數(shù)值模擬分析［J］.人民長江，2023，54（6）：112-119.

［3］" UYSAL M，TOPRAK A S，POLAT N.DEM generation with UAV photogrammetry and analysis in Sahitler hill［J］.Measurement，2015，73（1）：539-543.

［4］" YAN J H，CHEN J P，ZHANG Y S，et al.Semi-automatic extraction of dangerous rock blocks from jointed rock exposures based on a discontinuity trace map［J］.Computers and Geotechnics，2023，156：105265.

［5］" 鄭紅霞，李凌昊，盛輝，等.基于無人機貼近攝影測量的剖面數(shù)字露頭系統(tǒng)研發(fā)及巖層面重構(gòu)方法［J］.中國石油大學(xué)學(xué)報（自然科學(xué)版），2023，47（2）：53-63.

［6］" VASUKI Y，HOLDEN E J，KOVESI P，et al.Semi-automatic mapping of geological structures using UAV-based photogra mmetricdata：an image analysis approach［J］.Computers amp; Geosciences，2014，69：22-32.

［7］" 孫琦，張文，趙云鵬，等.高陡巖質(zhì)斜坡結(jié)構(gòu)面的無人機貼近采集與粗糙度分析［J］.工程地質(zhì)學(xué)報，2021，29（5）：1460-1468.

［8］" 張釗，吳鋒，尚海興，等.貼近攝影與傾斜攝影測量技術(shù)融合在水電站高壩精細(xì)建模中的應(yīng)用［J］.西北水電，2021（5）：47-50，55.

［9］" 張林杰，馮剛，龍超，等.無人機傾斜攝影與貼近攝影融合建模在高位危巖體調(diào)查中的應(yīng)用［J］.水電與抽水蓄能，2023，9（增1）：85-90.

［10］買小爭，馮學(xué)勝，顏振能，等.貼近攝影影像與傾斜影像融合精細(xì)化建模方法［J］.北京測繪，2022，36（9）：1188-1192.

［11］王豪，黃健，黃祥，等.一種利用Unity3D模擬崩塌三維運動全過程的方法［J］.武漢大學(xué)學(xué)報（信息科學(xué)版），2023，48（12）：1990-1998.

［12］呂權(quán)儒，曾斌，孟小軍，等.基于無人機傾斜攝影技術(shù)的崩塌隱患早期識別及影響區(qū)劃分方法［J］.地質(zhì)科技通報，2021，40（6）：313-325，334.

［13］鄭史芳，黎治坤.結(jié)合傾斜攝影技術(shù)的地質(zhì)災(zāi)害監(jiān)測［J］.測繪通報，2018（8）：88-92.

［14］吳明堂，姚富潭，楊建元，等.基于傾斜攝影測量與InSAR技術(shù)的庫區(qū)滑坡識別［J］.人民長江，2023，54（3）：130-137.

［15］何佳男.貼近攝影測量及其關(guān)鍵技術(shù)研究［D］.武漢：武漢大學(xué)，2022.

［16］葉震，許強，劉謙，等.無人機傾斜攝影測量在邊坡巖體結(jié)構(gòu)面調(diào)查中的應(yīng)用［J］.武漢大學(xué)學(xué)報（信息科學(xué)版），2020，45（11）：1739-1746.

［17］楊力龍.基于輕小型無人機的航空攝影測量技術(shù)在高陡邊坡幾何信息勘察中的應(yīng)用研究［D］.成都：西南交通大學(xué)，2017.

［18］徐剛，張志斌，范澤英.危巖失穩(wěn)模式分類研究［C］∥中國地質(zhì)學(xué)會，中國地質(zhì)學(xué)會工程地質(zhì)專業(yè)委員會.2017年全國工程地質(zhì)學(xué)術(shù)年會論文集.北京：科學(xué)出版社，2017.

［19］張倬元，王士天.工程地質(zhì)分析原理［M］.4版.北京：地質(zhì)出版社，2016.

［20］趙中強，王有林，王文革.金川水電站路基邊坡穩(wěn)定性分析與評價［J］.水利水電快報，2022，43（5）：58-62.

［21］陳祖煜.巖質(zhì)邊坡穩(wěn)定分析原理·方法·程序［M］.北京：中國水利水電出版社，2005.

（編輯：劉 媛）

Simulation of dangerous rock mass movement and delineation of impacted zone based on UAV fusion modeling

LIU Binging1，WU Fu2，LIU Chang2，HE Na2，HUANG Hai2

（1.State Key Laboratory of Geological Disaster Prevention and Environmental Protection，Chengdu University of Technology，Chengdu 610059，China；" 2.Geological Environment Monitoring Station of Guangxi Zhuang Autonomous Region，Nanning 530201，China）

Abstract：

Rock falls are influenced by a variety of factors and inherent uncertainties.Investigating potential disaster zones and distinguishing varying degrees of impact can provide critical data for disaster prevention and control.This study focused on the dangerous rock masses in Neizhuangtun，Nanning City，Guangxi.A 3-dimensional model was constructed using a fusion of oblique photography and nap-of-the-object photogrammetry modeling technology to accurately capture the structural and characteristic information of the rock mass.Through 3D numerical simulation，the movement path，process，real-time kinetic energy，and bouncing height of the dangerous rock masses were analyzed.Based on the simulation results，a method of dividing the potentially impacted areas was proposed.The findings indicated that the fusion modeling technique significantly enhanced the visibility of detailed features in steep mountainous regions，facilitating the accurate extraction of simulation parameters for rock movement.By analyzing multiple simulation results，the potential disaster risk zones were identified.The proposed method for dividing the affected areas enables a quick and easy division of potential disaster zones.The outcomes of this study provide insights for the investigation of high and steep dangerous rock masses.

Key words：

high dangerous rocks； fusion modeling； oblique photography； nap-of-the-object photogrammetry； numerical simulation； rockfall area zoning