2025-09-05 3524 0
高陡巖質(zhì)邊坡低傾角裂縫三維連通性計算方法
3D connectivity calculation method for low-dip fractures in high-steeprock slopes

摘要:
低傾角裂縫是控制高陡巖質(zhì)邊坡穩(wěn)定性的主要裂縫之一,直接影響拱壩壩肩的抗滑穩(wěn)定性。巖體中裂隙的空間復(fù)雜性和隨機性給三維連通性的精確測定帶來了挑戰(zhàn)。本文提出了一種計算三維連通性的新方法。通過修正裂縫與測量窗口夾角造成的采樣偏差,采用基于copula的二元分布模型對裂縫產(chǎn)狀進行建模,構(gòu)建了能準確反映裂縫空間分布的三維離散裂縫網(wǎng)絡(luò)(DFN)模型。采用最優(yōu)路徑搜索算法(A*算法)識別DFN模型中不同斷面的潛在破壞路徑,并采用離散插值分析預(yù)測不同高程下的潛在三維剪切面。然后使用布爾運算將涉及穿透破壞的裂縫投影到三維表面上,最終獲得反映工程巖體破壞模式的三維連通性。應(yīng)用該方法對西藏某水電工程壩肩巖體低傾裂縫連通性進行了評價。數(shù)值仿真結(jié)果驗證了該方法的準確性。此外,通過與傳統(tǒng)二維方法的對比分析,進一步驗證了所提出的三維連通率計算方法的可靠性。本文為水電工程壩肩巖體的抗滑穩(wěn)定性和地基面的選擇提供了參考,并對巖體裂隙的三維連通性特征有了新的認識。
引言:
裂縫,也稱為節(jié)理,是沿應(yīng)力導致巖石局部失去凝聚力的平面,通常兩側(cè)巖石都沒有明顯的位移。隨機分布的裂隙對控制巖體的基本力學行為和整體穩(wěn)定性起著至關(guān)重要的作用。它們不僅影響巖體的強度,而且是工程失穩(wěn)的潛在幾何和力學邊界。在實際工程中,巖體抗剪強度問題往往涉及節(jié)理和巖橋的聯(lián)合抗剪強度,這在很大程度上依賴于沿潛在破壞面斷裂連通性的準確計算。特別是在高陡巖質(zhì)邊坡的穩(wěn)定性分析中,低傾角裂縫的連通性特征對壩肩巖體的抗滑穩(wěn)定性評價起著關(guān)鍵作用
低傾角裂縫是構(gòu)成邊坡運動的主要構(gòu)造面之一,其連通性直接影響巖體的整體強度和穩(wěn)定性。然而,由于巖體中裂縫發(fā)育的空間復(fù)雜性和隨機性,以往的研究人員發(fā)現(xiàn),準確量化裂縫的連通性是一項挑戰(zhàn)。他們保守地估計了巖體的強度,假設(shè)裂縫已經(jīng)完全穿透。這種方法雖然簡單可行,但嚴重低估了巖體的強度。為了提高巖體的利用率,提出了二維裂隙連通性的概念。
研究人員最初將二維連接率定義為實測裂縫的平均軌跡長度等于貫通路徑長度之和(裂縫的平均軌跡長度加上巖橋的平均長度)。然而,由于現(xiàn)場測量窗口大小的限制,獲得裂縫的真實軌跡長度和巖橋之間的間距是極其困難的。窗口估計模型(H-H法),Wu等人根據(jù)測量窗口內(nèi)可見、部分可見和不可見裂縫痕跡的比例估計連通性,顯著解決了節(jié)理間距估計問題。此外,還提出了帶寬投影法計算連通性,獲得裂縫在剪切方向上的發(fā)育特征。帶寬投影法假定巖體通常不是沿單個裂縫破壞,而是由分布在特定帶寬內(nèi)的多個裂縫控制。因此,2D連通性速率可以定義為假設(shè)投影基線上指定帶寬內(nèi)所有裂縫投影長度之和與投影基線長度之比。在尚未被破壞的巖體中確定投影計算的帶寬仍然是一個挑戰(zhàn)。盡管H-H法和帶寬投影法都簡單而優(yōu)越,它們在三維分析計算中的適用性受到限制,因為它們沒有考慮裂縫的三維空間特征。
根據(jù)二維連通性的定義,三維連通性是指裂縫在參考平面上的投影面積與參考平面面積之比。然而,目前的技術(shù)仍然無法精確表征裂縫面復(fù)雜的3D形態(tài)。為了獲得裂縫的空間分布特征,利用物探、三維激光掃描、無人機攝影測量、圖像分析等多種野外調(diào)查技術(shù)獲取裂縫信息。然而,這些技術(shù)只捕獲了淺層表面上一些發(fā)育的大裂縫的形態(tài),并且所獲得的連通性率的準確性一直存在爭議。提出了FERM(巖體法醫(yī)開挖)方法,通過向巖體中注入可膨脹化學劑來重建裂縫的三維形態(tài),從而獲得三維連通性率。該方法可以獲得真實有效的裂縫連通性特征,但成本較高,不適用于大規(guī)模深部巖體裂縫。在過去的半個世紀里,基于統(tǒng)計的離散裂縫網(wǎng)絡(luò)(DFN)建模方法逐漸成為獲取裂縫三維連通性特征的主流方法。Priest等人首先提出使用計算機程序生成以數(shù)字或圖形表示的不連續(xù)裂縫,隨后被廣泛應(yīng)用于巖土工程。大量的實例研究表明,巖體中裂縫的幾何參數(shù)(如形狀、產(chǎn)狀、跡長和密度)往往遵循一定的概率分布特征。蒙特卡羅裂縫網(wǎng)模擬可以根據(jù)所測道參數(shù)的分布模式對不確定性裂縫進行預(yù)測和模擬,通過增加模擬次數(shù)來解決構(gòu)造真實裂縫網(wǎng)絡(luò)圖的問題。因此,確定裂縫的幾何分布特征對于建立三維DFN模型至關(guān)重要。許多學者研究了裂縫的幾何形狀,假設(shè)DFN模型中的裂縫幾何形狀可以表示為薄圓盤或固定長寬比的四邊形。在裂縫姿態(tài)表征方面,F(xiàn)isher引入了基于圍繞平均姿態(tài)旋轉(zhuǎn)對稱的二元Fisher分布。該分布只包含一個濃度參數(shù),并且是可積的,因此應(yīng)用廣泛。然而,它最適合具有明確單峰濃度的數(shù)據(jù)集。當濃度參數(shù)較大時,分布的尾部衰減迅速,限制了其表示離群態(tài)度的能力。相比之下,Bingham分布在表示具有多模態(tài)或橢圓特征的數(shù)據(jù)時提供了更大的靈活性,因為它由對稱矩陣控制。然而,它的應(yīng)用需要復(fù)雜的數(shù)據(jù)處理和龐大的數(shù)據(jù)集來準確估計參數(shù)。對于三維DFN中的跡線長度和密度估計,Priest等人提出了不同的計算方法來估計平均跡線長度、空間維度和裂縫密度。
巖體內(nèi)隨機離散分布的裂縫在三維空間中相交,形成高度復(fù)雜的網(wǎng)絡(luò)。盡管3D DFN建模技術(shù)已廣泛應(yīng)用于可視化裂縫的三維形態(tài),并提供了計算三維連通性的基礎(chǔ)模型,不同的研究人員提出了明顯不同的參數(shù)化方法。此外,許多研究忽略了抽樣偏差對模型精度的影響,使人們對所構(gòu)建的DFN模型的可靠性產(chǎn)生懷疑。此外,目前計算三維裂縫連通性的方法通常會將定義區(qū)域內(nèi)的所有裂縫不加選擇地投影到參考平面上。這種方法不可避免地導致計算出的連通性與現(xiàn)場觀察到的裂縫網(wǎng)絡(luò)的實際穿過特性之間存在差異。總的來說,對于表征深部巖體中裂縫的三維分布和連通性,目前還缺乏系統(tǒng)有效的解決方案。
針對這些問題,本研究提出了一種基于裂縫面與采樣面幾何相交的裂縫采樣偏差校正方法,從而構(gòu)建更能反映裂縫系統(tǒng)空間特征的三維DFN模型。在此改進的三維DFN模型的基礎(chǔ)上,首次應(yīng)用A*算法在三維空間中識別潛在的穿過剪切面。這種方法可以選擇不同海拔高度的關(guān)鍵裂縫,減少傳統(tǒng)的基于投影的方法的局限性。將離散點插值與布爾運算相結(jié)合,提出了一種計算三維裂縫的新方法連接建立。該方法在西藏某水電站壩肩巖體低傾裂縫評價中得到了應(yīng)用和驗證。該方法有效解決了傳統(tǒng)連通性計算的若干固有局限性,進一步拓展了三維DFN模型在剪切路徑識別和裂縫連通性定量評價方面的應(yīng)用潛力。為分析巖體裂隙的空間連通性特征提供了新的視角。
圖表:

圖1 二維連通性計算方法示意圖

圖2 三維裂縫連通性計算方法流程圖

圖3 裂縫與測量窗口的相交關(guān)系

圖4 圓盤的中心、半徑和軌跡長度之間的幾何關(guān)系示意圖

圖5 DFN模型生成過程:(a)蒙特卡羅法生成DFN模型;(b)確定性大型裂縫坐標計算

圖6 DFN模型的二維軌跡生成過程

圖7 A*算法搜索潛在故障路徑的實現(xiàn)過程

圖8 三維剪切面生成過程:(a)空間離散點插值分析;(b)位勢三維剪切面;(c)骨折投影面積;(d)剪切面剩余面積

圖9 研究領(lǐng)域的基本特征:(a)地理位置;(b)右岸壩肩現(xiàn)場照片

圖10 裂縫信息獲取流程:(a)裂縫面初始識別;(b)點云提取和平面擬合;(c)坡面三維激光掃描點云圖;(d)裂縫跡線圖
表1特殊情況下裂縫產(chǎn)狀計算方法


圖11 裂縫極點立體輪廓圖(Schmidt網(wǎng)):(a)現(xiàn)場測量結(jié)果,(b) K-means聚類結(jié)果,(c) copula模擬結(jié)果

圖12 裂縫采樣頻率校正結(jié)果:(a)校正頻率與傾角散點圖;(b)校正頻率和傾斜方向散點圖
表2 DFN模型低傾角裂縫基本參數(shù)


圖13 低角度裂縫三維連通性計算過程:(a)三維DFN模型及計算截面;(b)潛在剪切面交點裂縫選擇與投影計算
表3不同標高低角度三維裂縫三維連通性結(jié)果


圖14 野外露頭與DFN模型露頭的跡線和姿態(tài)比較

圖15 DFN模型跡線長度與現(xiàn)場跡線長度的比較:(a)集合1;(b)組2;(c)組3;(d)第4組
表4 DFN模型走線長度與現(xiàn)場走線長度的偏差


圖16 抗剪強度參數(shù)箱線圖:(a)斷裂面;(b)完整巖石
表5巖石橋梁與裂隙力學參數(shù)取值


圖17 數(shù)值模型及試驗加載方法:(a)數(shù)值模型;(b)裝載方式

圖18 不同圍壓和莫爾庫侖強度包絡(luò)線下的軸向應(yīng)力-應(yīng)變曲線:(a)不同圍壓下的軸向應(yīng)力-應(yīng)變曲線;(b)莫爾-庫侖強度包絡(luò)線
圖19 右岸壩肩坑道低傾裂縫示蹤圖(局部視圖)

圖20 坑道內(nèi)低傾裂縫二維連通性計算結(jié)果:(a)不同帶寬下二維連通性計算結(jié)果;(b)二維與三維連通性計算結(jié)果對比

圖21 裂縫連通性的方向敏感性分析:(a)切片方向示意圖。(b)不同剪切方向下連通性玫瑰圖
表6不同剪切方向連通性及穩(wěn)定性計算結(jié)果


圖22 連通性與穩(wěn)定性因子的關(guān)系
結(jié)論:
研究了巖體中裂縫的空間連通性特征,提出了一種計算三維裂縫連通性的新方法。本研究的主要貢獻總結(jié)如下:
(1)提出了基于裂縫與露頭面相交概率的裂縫采樣頻率修正模型,為后續(xù)三維連通性計算提供了可靠數(shù)據(jù)。
(2)引入聯(lián)結(jié)二元分布函數(shù)來表示裂縫走向與傾角之間的相關(guān)性,克服了傳統(tǒng)方法單獨處理這些參數(shù)的局限性。
(3)利用統(tǒng)計建模參數(shù)和修正后的裂縫頻率,建立了包含確定性大尺度裂縫和露頭的更精確的三維DFN模型。A*算法結(jié)合離散點插值分析預(yù)測潛在的三維剪切面,布爾運算進行三維投影計算裂縫涉及潛在的破壞路徑,從而產(chǎn)生更精確的三維連通性計算方法。
(4)將該方法與傳統(tǒng)的二維計算方法相結(jié)合,對西藏某水電站壩肩巖體低傾裂縫連通性進行了評價。連通性結(jié)果分別為37.3%和47.6%。與機械反分析得到的值(32.5%)相比,該方法產(chǎn)生的偏差僅為4.8%,顯示出更高的精度。
(5)三維連通性計算方法具有方向敏感性,連通性隨剪切方向與邊坡走向夾角的增大而增大。這種行為受到邊坡中主要裂縫產(chǎn)狀的影響,并突出了不僅根據(jù)幾何形狀,而且根據(jù)巖體的地應(yīng)力條件選擇剪切方向的重要性。
參考文獻
[1]Liang P, Wei Y, Guo H, et al. 3D connectivity calculation method for low-dip fractures in high-steep rock slopes[J]. Engineering Geology, 2025, 355: 108227.
Hot News
成功提示
錯誤提示
警告提示
評論 (0)