巖質(zhì)邊坡穩(wěn)定分析一站式解決方案

呂 超，孫旭良，廖慶龍，李 琦

（廣西桂禹工程咨詢有限公司，廣西 南寧 530023）

0 引言

巖體結(jié)構(gòu)面產(chǎn)狀的量測、優(yōu)勢結(jié)構(gòu)面統(tǒng)計、結(jié)構(gòu)面對邊坡穩(wěn)定性的影響分析是工程地質(zhì)勘察中重要的基礎工作。目前我國基礎建設迅猛發(fā)展，水電工程、鐵路工程、高速路工程等很多項目在勘察設計階段和施工階段均會涉及巖體工程，而巖質(zhì)邊坡穩(wěn)定性分析評價是其中的重要內(nèi)容。這部分工程關注巖體結(jié)構(gòu)面產(chǎn)狀的量取、結(jié)構(gòu)面產(chǎn)狀統(tǒng)計分析、結(jié)構(gòu)面穩(wěn)定性分析、地質(zhì)評價報告編制等內(nèi)容。

當前野外地質(zhì)勘察中，測量巖體結(jié)構(gòu)面信息的主要工具是地質(zhì)羅盤和皮尺。這種方法在所有勘察方法中最為普遍，但也存在一定的不足之處，工作量大、費時費力，調(diào)查人員無法接觸的區(qū)域?qū)y量造成阻礙，人員安全會受到威脅。另外一種接觸式測量方法是鉆探方法，這種方法成孔的質(zhì)量要求較高，范圍小，精度低。近年來興起的近景攝影測量法和三維激光掃描儀快速獲取點位坐標法的應用，一定程度上給地質(zhì)調(diào)查帶來了一種新的革命，但這兩種方法也存在弊病，操作復雜、需要其他輔助設備精確定位、搬運費時費力、設備昂貴等。

小型無人機的問世給地質(zhì)調(diào)查帶來了新希望，小型無人機自身帶有定位系統(tǒng)，價格低廉，操作簡單，數(shù)據(jù)精度可用，但數(shù)據(jù)為高程和坐標，需要通過空間幾何算法計算結(jié)構(gòu)面傾向和傾角。結(jié)構(gòu)面產(chǎn)狀統(tǒng)計分析包括產(chǎn)狀分組統(tǒng)計、各組百分比含量、節(jié)理玫瑰圖繪制、優(yōu)勢結(jié)構(gòu)面選取；結(jié)構(gòu)面穩(wěn)定性分析主要應用赤平投影方法對結(jié)構(gòu)面及其組合的交割線進行穩(wěn)定性分析評價。這兩種方法，目前理正巖土軟件、邁達斯GTS 等商用軟件都有獨立模塊進行統(tǒng)計與分析。

但目前軟件對三維實景模型結(jié)構(gòu)面產(chǎn)狀計算、結(jié)構(gòu)面產(chǎn)狀統(tǒng)計分析、結(jié)構(gòu)面穩(wěn)定性分析、地質(zhì)分析評價報告自動編制等工作之間的結(jié)合尚不緊密，工程建設隨著進度推進，結(jié)構(gòu)面數(shù)據(jù)的變化，相關統(tǒng)計分析工作都會變動，報告成果也會改變，使用不夠便捷。

因此開發(fā)一套針對巖質(zhì)邊坡穩(wěn)定性進行分析的軟件，以完成上述分析的全部內(nèi)容，優(yōu)化結(jié)構(gòu)面?zhèn)鹘y(tǒng)測量方法，使結(jié)構(gòu)面產(chǎn)狀計算、結(jié)構(gòu)面產(chǎn)狀統(tǒng)計、結(jié)構(gòu)面穩(wěn)定性分析、地質(zhì)評價報告自動編制等工作一站式完成，無論結(jié)構(gòu)面數(shù)據(jù)如何更新、分析邊坡部位如何變化，均可快速可靠的一鍵完成穩(wěn)定性分析，并生成地質(zhì)評價報告，使得各工序無縫銜接，促進成果質(zhì)量和工作效率的大幅提高，對地質(zhì)工作數(shù)字化標準化具有重要意義。

1 穩(wěn)定分析一站式解決方案

1.1 技術路線

穩(wěn)定分析工作需要完成巖層產(chǎn)狀量測、結(jié)構(gòu)面統(tǒng)計分析、邊坡穩(wěn)定性分析和生成報告幾部分工作。本次以上述分析工作為研究對象，以高效、高質(zhì)量完成設計分析工作為目標，以系統(tǒng)思維為統(tǒng)領，以標準化為準繩，提出標準化設計流程和成果范本，以計算機技術為手段，開發(fā)分析軟件，為供結(jié)構(gòu)面分析提供一站式解決方案。分析軟件按照其功能劃分為結(jié)構(gòu)面產(chǎn)狀采集模塊、產(chǎn)狀統(tǒng)計分析模塊、邊坡穩(wěn)定分析和報告生成4 個模塊。

該軟件可實現(xiàn)從三維實景模型提取巖體結(jié)構(gòu)面坐標，并完成結(jié)構(gòu)面產(chǎn)狀計算（傾向與傾角）；在此基礎上統(tǒng)計結(jié)構(gòu)面分布規(guī)律并分析優(yōu)勢結(jié)構(gòu)面，同時完成節(jié)理玫瑰花圖繪制；在上一步結(jié)構(gòu)面統(tǒng)計分析的基礎上，依據(jù)赤平投影方法分析結(jié)構(gòu)面及其交割線對邊坡穩(wěn)定性的影響，并自動繪制赤平投影圖；最后，整合上述3個過程，生成分析報告。

Excel是美國微軟公司開發(fā)的商業(yè)辦公軟件，具有強大的數(shù)據(jù)計算分析能力和豐富的圖表工具，同時提供了二次開發(fā)功能，可實現(xiàn)定制服務。Excel 的二次開發(fā)有多種方式，包括著名的VBA（VB for Application）、VSTO（Visual Studio Tools for Office）和DNA（DotNet for Application）等。本次開發(fā)采用DNA 技術，該技術是基于.Net 框架的開發(fā)方式，可以依托.Net 的優(yōu)勢，快速開發(fā)出功能強大的應用程序，與其他開發(fā)方式相比具有安全性好，不易破解，無需安裝注冊，程序的分發(fā)和部署非常便利等優(yōu)勢。

軟件以Excel 工作表為輸入輸出平臺，將工作表劃分為輸入?yún)^(qū)、輸出區(qū)和圖表區(qū)，開發(fā)相應的應用程序，完成上述分析計算、圖表繪制和報告生成工作。軟件開發(fā)技術路線如圖1所示。

圖1 軟件開發(fā)技術路線圖Fig.1 Technology roadmap of software development

1.2 結(jié)構(gòu)面產(chǎn)狀提取和計算

1.2.1 基本原理

在三維實景模型中提取巖體結(jié)構(gòu)面上不在同一直線的3個點的坐標P1（x1，y1，z1）、P2（x2，y2，z2）、P3（x3，y3，z3），可通過空間幾何計算出該巖體結(jié)構(gòu)面的產(chǎn)狀（傾向α和傾角A）。

巖體結(jié)構(gòu)面的傾向α需要確定該結(jié)構(gòu)面單位法向量n→在XY平面的投影n′所在的象限，具體計算如下：

（1）當nz≥0時：

巖體結(jié)構(gòu)面的傾角A為單位法向量n→（nx，ny，nz）與Z軸的夾角，具體計算如下：A=arccos［nz］。

其中，巖體結(jié)構(gòu)面的單位法向量n→（nx，ny，nz）計算如下：

1.2.2 坐標提取

（1）復制坐標。為從實景模型中提取坐標，本次編制了提取程序，具體步驟為：在實體模型中選擇巖體結(jié)構(gòu)面上不在同一直線上的3 個點，依次選擇點坐標，并將其復制到剪貼板，如圖2所示。

圖2 坐標點選擇及復制Fig.2 Selection and copy of coordinate points

（2）提取坐標。在選項卡中點擊“提取產(chǎn)狀”按鈕，可彈出“結(jié)構(gòu)面提取”對話框，完成上述復制坐標后，點擊對話框中的“提取點”按鈕，程序?qū)⒃擖c坐標填充到工作表中的指定位置，重復1）和點擊“提取點”按鈕可完成其他坐標點的提取。當完成結(jié)構(gòu)面提取后，點擊“計算”按鈕，程序自動完成傾向和傾角計算公式的填充，并完成傾角傾向計算，如圖3所示。

圖3 坐標提取Fig.3 Extraction of coordinates

1.2.3 產(chǎn)狀計算

結(jié)構(gòu)面計算需要輸入的條件是3 個不在一條直線上的3 個點的坐標，共9 個參數(shù)，計算結(jié)果為傾向和傾角。在Excel 工作表中編制標準表格，第A 列為序號，第B 列為結(jié)構(gòu)面名稱，第C～K 列為3個點的三維坐標，第L、M 列為傾向和傾角。A～K 列為輸入?yún)^(qū)，L、M為計算結(jié)果輸出區(qū)，表格形式見圖3。

為實現(xiàn)傾向和傾角的計算，按照上述計算原理，編制2個自定義函數(shù)，分別為傾向計算函數(shù)：QingXiang（C6：K6）和傾角計算函數(shù)QingJiao（C6：K6）。2個函數(shù)輸入?yún)?shù)為結(jié)構(gòu)面的3個點的9 個坐標參數(shù)。輸入?yún)?shù)后，程序自動完成結(jié)構(gòu)面傾向和傾角2個參數(shù)的計算。

Function QingXiang（P1ToP3 As Object（））As Single ′傾向計算

Dim xyz（9）As Single

For i=1 To 9

xyz（i）=P1ToP3（i）

Next i

Dim qj As Single

Dim qx As Single

CalQingJiao xyz，qj，qx

Return qx

End Function

Function QingJiao（P1ToP3 As Object（））As Single ′傾角計算

For i=1 To 9

xyz（i）=P1ToP3（i）

Next i

Dim qj As Single

Dim qx As Single

CalQingJiao xyz，qj，qx

Return qj

End Function

Private Function CalQingJiao（xy，QingJiao As Single，QingXiang As Single）

……

If nz ＞ 0 Then

QingJiao=Application.WorksheetFunction.Acos（nz）* 180 /pi

Else

QingJiao=180-Application.WorksheetFunction.Acos（nz）*180 / pi

End If

′傾向計算

ddd=（nx * nx+ny * ny）^ 0.5

If nz ＞=0 Then

If nx ＞ 0 Then

QingXiang=Application.WorksheetFunction.Acos（ny / ddd）* 180 / pi

Else

QingXiang=360-Application.WorksheetFunction.Acos（ny /ddd）* 180 / pi

End If

Else ′nz＜0

If nx ＜ 0 Then

QingXiang=Application.WorksheetFunction.Acos（-ny /ddd）* 180 / pi

Else

QingXiang=360-Application.WorksheetFunction.Acos（-ny/ ddd）* 180 / pi

End If

End If

End Function

1.3 結(jié)構(gòu)面產(chǎn)狀統(tǒng)計分析

當工程區(qū)巖層、節(jié)理、裂隙等結(jié)構(gòu)面發(fā)育，結(jié)構(gòu)面產(chǎn)狀數(shù)量較多時，必須統(tǒng)計結(jié)構(gòu)面分布規(guī)律，分析主要的優(yōu)勢結(jié)構(gòu)面，為巖質(zhì)邊坡穩(wěn)定性分析提供重要數(shù)據(jù)。程序?qū)?60°分成36組，每10°傾向為一組統(tǒng)計，分析統(tǒng)計每組內(nèi)的平均傾向、平均傾角、每組裂隙數(shù)量占總數(shù)量的百分比；同時為更加直觀反映結(jié)構(gòu)面分布規(guī)律，軟件可根據(jù)統(tǒng)計結(jié)果自動繪制傾向節(jié)理玫瑰花圖和傾角分布直方圖。

為實現(xiàn)上述功能，在Excel 工作表編制結(jié)構(gòu)面統(tǒng)計標準表格，第Q～X 列分別為編號、傾向分組、平均傾角、百分比、傾向、累計條數(shù)、改組條數(shù)等，表格形式如下圖。之后，在表格中輸入Excel 提供的內(nèi)部函數(shù)即可完成統(tǒng)計，其中累計條數(shù)列輸入“COUNTIF（$L$6：$L$500，V6）”函數(shù)，該函數(shù)第1 個參數(shù)為前面計算的傾向列，第2 個參數(shù)為判斷條件，位于V 列，其他項目計算公式略。標準表格及圖表如圖4所示。

圖4 結(jié)構(gòu)面統(tǒng)計成果圖表Fig.4 Results chart of statistics of discontinuities

1.4 巖質(zhì)邊坡穩(wěn)定性分析程序

1.4.1 原 理

根據(jù)優(yōu)勢結(jié)構(gòu)面分析結(jié)果，待評價邊坡的坡向和坡角，采用赤平投影方法，對巖質(zhì)邊坡的穩(wěn)定性進行初步分析。為與《建筑邊坡工程技術規(guī)范》（GB 50330-2013）和《滑坡防治工程勘查規(guī)范》（GB/T 32864-2016）的邊坡穩(wěn)定狀態(tài)劃分相適應，本程序?qū)⑦吰路€(wěn)定性劃分為四級：不穩(wěn)定、欠穩(wěn)定、基本穩(wěn)定和穩(wěn)定，穩(wěn)定性分級和評價指標見表1。

表1 邊坡穩(wěn)定性分級與評價指標Tab.1 Classification and evaluation index of slope stability

1.4.2 實現(xiàn)方法

赤平投影分析是在上步結(jié)構(gòu)面產(chǎn)狀統(tǒng)計分析的基礎上，針對選定的結(jié)構(gòu)面進行赤平投影圖繪制和分析，提出穩(wěn)定分析結(jié)論。其輸入條件坡面和巖層結(jié)構(gòu)面產(chǎn)狀。

（1）定標準表格。在Excel 工作表中編制標準表格，設置輸入?yún)^(qū)、圖表區(qū)和成果輸出區(qū)。本程序AL4：AP14區(qū)域，右側(cè)為圖表區(qū)，下部為成果輸出區(qū)。表格形式見圖6。

（2）繪制赤平投影圖。根據(jù)輸入的坡面和巖層產(chǎn)狀信息，按照赤平投影原理和計算公式，編制程序生成基圓、結(jié)構(gòu)面投影弧線的離散的數(shù)據(jù)點和交割線的坐標點，自動寫入工作表的指定區(qū)域并完成圖表繪制。

（3）優(yōu)勢節(jié)理面選擇。為銜接產(chǎn)狀統(tǒng)計和赤平投影分析，本次開發(fā)了赤平投影計算模塊，該模塊可以實現(xiàn)節(jié)理面組合快速選擇和參數(shù)輸入，之后完成赤平投影分析，界面如圖5所示。

圖5 赤平投影計算對話框Fig.5 Calculation dialog box of stereographic projection

（4）赤平投影主要成果。完成赤平投影分析，程序完成各分區(qū)結(jié)構(gòu)面（交割線）數(shù)量統(tǒng)計，形成分析結(jié)論。結(jié)構(gòu)面參數(shù)、分析成果及有關圖表如圖6所示。

圖7 程序加載界面Fig.7 Interface of program

1.5 報告生成

根據(jù)專業(yè)特點，提出了報告范本，開發(fā)相應程序，根據(jù)上述產(chǎn)狀提取、統(tǒng)計分析和赤平投影分析的參數(shù)、圖表等成果，自動生成專題報告。程序略。

1.6 Excel-DNA開發(fā)

1.6.1 開發(fā)環(huán)境

本次開發(fā)采用Excel-DNA 版本是ExcelDna-0.34.6，其開發(fā)包可從網(wǎng)上下載。將下載的ExcelDna-0.34.6.zip 壓縮包解壓到指定目錄下，本程序為“d：sxlExcelDna-0.34.6”目錄，開發(fā)包中包含ExcelDna.xll、ExcelDna64.xll、ExcelDnaPack.exe 和ExcelDna.Integration.dll這4個核心文件。各文件功能如表2所示。

表2 Excel-DNA開發(fā)包主要文件及其功能Tab.2 The Key Files and their Functions

1.6.2 加載項的生成

（1）文件準備。Excel-DNA 加載項的生成需要文件名相同，擴展名不同的3 個文件，本程序為：dizhiFenxi.dna、dizhi-Fenxi.dll、dizhiFenxi.xll。其中dizhiFenxi.dna 文件為文本文件，

dizhiFenxi.xll 是將開發(fā)包中的ExcelDna.xll 復制后修改文件名得到。（如果是64 位Excel，則將ExcelDna64.xll 復制后修改文件名得到。）

dizhiFenxi.dll動態(tài)鏈接庫，為程序的核心，由本次開發(fā)程序編譯后得到。

（2）dizhiFenxi.dll的生成。由于Visual Studio 開發(fā)平臺提供了豐富的功能，有利于開發(fā)和調(diào)試，本次選擇其進行dizhiFenxi.dll的開發(fā)和編譯。操作步驟簡述如下：①啟動Visual Studio，在“新建項目”對話框中，選擇Visual Basic 語言類型庫項目，并指定項目名稱；②在項目中添加ExcelDna.Integration 引用，位置為開發(fā)包目錄，本程序為“D：sxlExcelDna-0.34.6Distribution”。之后，打開Class1.vb 文件，在文件頂部導入ExcelDna.Integration指令。③完成上述操作后，可按照開發(fā)與常規(guī)VB.Net 相同，不在詳述。④完成程序編制后，選擇菜單“生成”→“生成解決方案”，Visual Studio會生成“dizhiFenxi.dll”文件。

（3）程序打包。準備好上述3 個文件后，即可進行打包，本次生成的最終加載項文件為“dizhiFenxi32.xll”。為方便操作，本次編制了打包批處理文件，內(nèi)容如下：

cd D：\sxl\dizhiFenxi\dizhiFenxi\bin\Debug\del dizhiFenxi.xll copy ExcelDna.xll dizhiFenxi.xll D：\sxl\ExcelDna-0.34.6\Distribution\ExcelDnaPack.exe dizhiFenxi.d na using bass addin dizhiFenxi.xll /y del dizhiFenxi32.xll copy dizhiFenxi-packed.xll dizhiFenxi32.xll pause

（4）程序加載。將加載項文件“dizhiFenxi32.xll”分發(fā)給用戶，啟動Excel 后，選擇“開發(fā)工具”選項卡中的“加載項”按鈕（如無“開發(fā)工具”選項卡，可右擊工具欄，然后選擇“自定義功能區(qū)”→勾選“開發(fā)工具”），可彈出如下加載宏對話框，瀏覽程序所在目錄，加載 “dizhiFenxi32.xll”。加載完成后，“功能區(qū)”中會自動生成本次添加的“地質(zhì)工具箱”選項卡，利用該選項卡中的按鈕即可進行分析計算。

2 應用與驗證

巖體結(jié)構(gòu)面產(chǎn)狀傳統(tǒng)量測方法以地質(zhì)羅盤為主，然而地質(zhì)羅盤量測也有很大的弊端，需要地質(zhì)調(diào)查人員現(xiàn)場對結(jié)構(gòu)面裂隙逐一量測，作業(yè)時間長、效率低下。而且只能獲得可以人為接觸的結(jié)構(gòu)面信息，從而導致結(jié)構(gòu)面量測范圍受地形限制，信息量有限、代表性不足，很大程度上影響巖體穩(wěn)定性分析評價結(jié)果的準確性。羅盤測量時指針擺動，讀數(shù)時系統(tǒng)誤差加上磁偏角誤差再加上識讀誤差一般會在5°～10°左右。地質(zhì)工作者多數(shù)情況下會碰到無法攀爬的高陡邊坡巖體，這時候往往需要有經(jīng)驗的地質(zhì)工作者手持羅盤比劃一個產(chǎn)狀信息，作為計算其穩(wěn)定性的基礎數(shù)據(jù)，這時的誤差在10°以上。而無人機三維實景模型上提取巖體結(jié)構(gòu)面數(shù)據(jù)進行產(chǎn)狀幾何計算不受高陡地形限制，量測范圍大、結(jié)構(gòu)面信息量足、代表性好。

本軟件以工程實例為背景，對結(jié)構(gòu)面產(chǎn)狀、結(jié)構(gòu)面統(tǒng)計及赤平投影分析進行了驗證，具體過程如下。

2.1 結(jié)構(gòu)面產(chǎn)狀計算驗證

采用精靈4 RTK多旋翼無人機，利用攝影測量3D規(guī)劃航線對工程區(qū)開展無人機航拍工作，獲得高清影像數(shù)據(jù)，并通過RTK 同步獲取野外布設的控制點坐標數(shù)據(jù)，利用控制點坐標數(shù)據(jù)進行空中三角測量之后，生成具有坐標和高程的三維實景模型。在三維實景模型中巖體結(jié)構(gòu)面上提取不在同一直線上的三點坐標，通過程序自動計算出該結(jié)構(gòu)面的傾向和傾角。通過本軟件自動計算和手算，結(jié)果驗證一致，說明軟件自動計算產(chǎn)狀結(jié)果精度良好，見表3。

表3 結(jié)構(gòu)面計算成果對比表Tab.3 Comparison table of discontinuities calculation result

三維實景模型中選取10組出露清晰的結(jié)構(gòu)面，計算其傾向和傾角，同羅盤測量獲取的傾向和傾角間的對比中誤差均在3°以內(nèi)，具有較好的精度，說明了無人機三維實景模型解譯產(chǎn)狀信息的可靠性，可以在三維實景模型中大量解譯巖體結(jié)構(gòu)面產(chǎn)狀信息應用于巖質(zhì)邊坡穩(wěn)定性分析中。

2.2 結(jié)構(gòu)面的統(tǒng)計

結(jié)構(gòu)面統(tǒng)計按每10°傾向為一組，共計36 組，自動統(tǒng)計每組內(nèi)的平均傾向、平均傾角、每組裂隙數(shù)量占總數(shù)量的百分比，通過本軟件計算和市場上某軟件計算，結(jié)果驗證一致，說明軟件計算成果精度好，程序可靠，對比表見表4。

表4 結(jié)構(gòu)面統(tǒng)計成果表Tab.4 The statistics result of discontinuities

2.3 巖質(zhì)邊坡穩(wěn)定性評價

根據(jù)表3 結(jié)構(gòu)面數(shù)據(jù)進行優(yōu)勢結(jié)構(gòu)面分析后，采用本軟件和市場上常用赤平投影分析軟件對擬開挖邊坡進行穩(wěn)定性分析，兩者穩(wěn)定性分析結(jié)論一致。說明本軟件具有較好的可靠性。

3 結(jié)論與建議

3.1 結(jié)論

（1）本軟件將巖層結(jié)構(gòu)面分析中的結(jié)構(gòu)面量測、結(jié)構(gòu)面統(tǒng)計、穩(wěn)定性分析及報告編制等工作進行了融合，提供了一站式解決方案，在系統(tǒng)化、標準化方面做了一些探索，促進了生產(chǎn)效率和成果質(zhì)量的大幅提升。

（2）巖層產(chǎn)狀采集由三維實景模型提取，改變了傳統(tǒng)工作方式，降低高陡邊坡產(chǎn)狀采集難度，有效避免了工作人員的危險性，同時工作效率也明顯提高。

（3）軟件采用Excel-DNA 技術開發(fā)，安全性較VBA 大幅提升，同時避免了VSTO 和VC 等開發(fā)技術對開發(fā)人員要求較高的問題，為非計算機專業(yè)人員提供了一個高效實用工具。

3.2 建 議

（1）巖質(zhì)邊坡結(jié)構(gòu)面數(shù)據(jù)的全面性和代表性直接影響穩(wěn)定性分析結(jié)論，如何進一步擴大地表以下巖體結(jié)構(gòu)面的采集范圍和數(shù)據(jù)量，是我們進一步研究的方向。

（2）在無人機航拍和三維實景模型生成過程中，如何有效剔除植被影響，高精度反映實際邊坡狀態(tài)有待進一步研究。