陳志杰,馮 曦,李傲松

(四川大唐國(guó)際甘孜水電開(kāi)發(fā)有限公司,四川康定 626001)

三維巖體裂隙網(wǎng)絡(luò)模擬研究及應(yīng)用

陳志杰,馮 曦,李傲松

(四川大唐國(guó)際甘孜水電開(kāi)發(fā)有限公司,四川康定 626001)

為合理評(píng)價(jià)某水庫(kù)庫(kù)岸巖體結(jié)構(gòu)特征及工程整治提供依據(jù)，對(duì)區(qū)內(nèi)巖體進(jìn)行了裂隙調(diào)查及三維網(wǎng)絡(luò)模擬研究。研究基于Monte-Carlo模擬方法，由走向玫瑰圖、傾向和傾角分布的直方圖對(duì)裂隙的產(chǎn)狀分布特征分別進(jìn)行了分析，結(jié)合裂隙等密度圖分析得出了結(jié)構(gòu)面優(yōu)勢(shì)產(chǎn)狀、用直接法產(chǎn)生隨機(jī)數(shù)并利用計(jì)算機(jī)模擬得出符合實(shí)際的三維裂隙網(wǎng)絡(luò)圖。引入廣義計(jì)算，討論了值的變化可能影響庫(kù)岸巖體在不同方向的穩(wěn)定性。

巖體裂隙；三維網(wǎng)絡(luò)；模擬研究

1 基本原理

從20世紀(jì)80年代初開(kāi)始研究裂隙網(wǎng)絡(luò)模擬技術(shù)[1]以來(lái)，裂隙網(wǎng)絡(luò)模擬技術(shù)已逐漸成熟，其成果已在工程上得到了廣泛的應(yīng)用[2，5]。裂隙網(wǎng)絡(luò)模擬研究過(guò)程一般包括：1）在野外采樣的基礎(chǔ)上對(duì)裂隙樣本進(jìn)行統(tǒng)計(jì)分析，包括對(duì)樣本進(jìn)行分組和各組樣本隨機(jī)變量（如走向、傾向、傾角、間距、跡長(zhǎng)等）的統(tǒng)計(jì)；2）對(duì)樣本分布形式進(jìn)行擬合優(yōu)度檢驗(yàn)，判斷各隨機(jī)變量的統(tǒng)計(jì)分布形式和分布參數(shù)；3）根據(jù)裂隙各隨機(jī)變量的統(tǒng)計(jì)分布形式，生成符合裂隙分布規(guī)律的隨機(jī)數(shù)，并以此生成裂隙網(wǎng)絡(luò)圖。

2 裂隙網(wǎng)絡(luò)統(tǒng)計(jì)與模擬

2.1 工程背景

某庫(kù)段位于興山縣中部偏南的香溪河谷，為興山縣原縣城所在地。區(qū)域上位于大巴山與秭歸盆地過(guò)渡地帶，屬于構(gòu)造剝蝕中低山區(qū)，山高坡陡，相對(duì)高差800～1 000 m。研究區(qū)內(nèi)地層出露齊全，從古元界至新生界均有出露，該地區(qū)分布地層有三疊系，侏羅系及第四系。區(qū)域地質(zhì)構(gòu)造上高陽(yáng)鎮(zhèn)處于3個(gè)構(gòu)造單元的交接部位，區(qū)內(nèi)地質(zhì)構(gòu)造復(fù)雜，次級(jí)褶皺斷層發(fā)育。

按該水庫(kù)設(shè)計(jì)水位，該地區(qū)將大部分被淹沒(méi)，將極大地改變?cè)摰貐^(qū)的工程地質(zhì)條件，進(jìn)一步惡化其環(huán)境地質(zhì)條件。為研究該庫(kù)段巖體地質(zhì)結(jié)構(gòu)特征，進(jìn)一步為地質(zhì)災(zāi)害防治工程設(shè)計(jì)施工提供依據(jù)，進(jìn)行此庫(kù)段巖體裂隙網(wǎng)絡(luò)模擬研究。

2.2 裂隙野外采樣

要進(jìn)行巖體裂隙網(wǎng)絡(luò)的計(jì)算機(jī)模擬，首先必須通過(guò)采樣由統(tǒng)計(jì)分析確定巖體中裂隙空間幾何參數(shù)的分布形式，需要的參數(shù)有：裂隙的空間產(chǎn)狀（走向、傾向和傾角）、裂隙的間距、裂隙的跡長(zhǎng)等；然后根據(jù)幾何參數(shù)的分布形式進(jìn)行裂隙網(wǎng)絡(luò)模擬。對(duì)裂隙的幾何參數(shù)進(jìn)行采樣，常見(jiàn)的測(cè)量方法有統(tǒng)計(jì)窗法和測(cè)線測(cè)量法，本次研究采用測(cè)線測(cè)量法。實(shí)際測(cè)量是選擇在研究區(qū)塌岸帶內(nèi)出露最為廣泛的香溪組、中—下統(tǒng)聶家山組地層內(nèi)布置了12條測(cè)線，共1 235條節(jié)理。

2.3 裂隙幾何參數(shù)分析統(tǒng)計(jì)

裂隙統(tǒng)計(jì)分析的方法很多，有玫瑰花圖法、直方圖法、極點(diǎn)圖法和等密度圖法等。本次研究中，首先由走向玫瑰圖、傾向和傾角分布的直方圖對(duì)裂隙的產(chǎn)狀分布特征分別進(jìn)行了分析，然后結(jié)合裂隙等密度圖分析得出了結(jié)構(gòu)面優(yōu)勢(shì)產(chǎn)狀。

1）裂隙走向分布特征。由圖1可以看出：NW組較NE組發(fā)育，且裂隙在NW象限內(nèi)的分布是比較均勻的；而對(duì)于NE組則普遍出現(xiàn)裂隙發(fā)育較分散的現(xiàn)象，NE組內(nèi)的裂隙主要集中在0°～40°之間。

2）傾向和傾角。由傾向和傾角分布直方圖可以看出：傾向大致呈均勻分布，但多在110°～130°之間和190°～240°之間；傾角大多為陡傾，一般在45°～80°之間。

3）裂隙產(chǎn)狀分布特征。根據(jù)結(jié)構(gòu)面的統(tǒng)計(jì)數(shù)據(jù)，可以作出結(jié)構(gòu)面的等密度圖。研究區(qū)內(nèi)普遍發(fā)育的結(jié)構(gòu)面有以下 5 組：15°∠80°，120°∠45°，120°∠70°，195°∠80°和 220°∠65°。

4）間距和跡長(zhǎng)。根據(jù)前人的工作[1]，裂隙的間距一般服從對(duì)數(shù)正態(tài)分布或負(fù)指數(shù)分布。在本研究區(qū)域內(nèi)，裂隙的間距服從對(duì)數(shù)正態(tài)分布，裂隙的平均間距較小，一般為0.25～0.30 m。裂隙的跡長(zhǎng)分布規(guī)律均不與四種常見(jiàn)分布（均勻分布、負(fù)指數(shù)分布、正態(tài)分布和對(duì)數(shù)正態(tài)分布）相吻合，且裂隙的平均半跡長(zhǎng)在1.3～1.6 m之間，規(guī)模相對(duì)較小，如表1所示。

表1 裂隙分布規(guī)律及模擬結(jié)果統(tǒng)計(jì)表

2.3 直接法及裂隙網(wǎng)絡(luò)模擬過(guò)程

本次研究的裂隙有相當(dāng)一部分隨機(jī)變量并不服從常見(jiàn)的分布形式（見(jiàn)表1），此時(shí)如仍采用常規(guī)Monte-Carlo法，則無(wú)法產(chǎn)生滿足精度要求的隨機(jī)數(shù)或無(wú)法繼續(xù)進(jìn)行裂隙模擬，而直接法則較好地解決了這個(gè)問(wèn)題。

2.3.1 隨機(jī)數(shù)的產(chǎn)生

從數(shù)學(xué)上來(lái)說(shuō)如只進(jìn)行裂隙網(wǎng)絡(luò)模擬直接法可以取代常規(guī)法，其省去了對(duì)隨機(jī)變量進(jìn)行是否服從某些分布形式的檢驗(yàn)和生成相應(yīng)隨機(jī)數(shù)的數(shù)學(xué)變換過(guò)程，方法更明了，且精度更高，其均值和方差更接近樣本，x2檢驗(yàn)值更小。下面重點(diǎn)介紹生成裂隙隨機(jī)數(shù)的直接法，具體過(guò)程為：

將樣本分布的區(qū)間 [a，b]劃分為n個(gè)小區(qū)間[a，a1]，[a1，a2],…，[an-1，b]，記樣本落入第 j個(gè)小區(qū)間的概率為pj（1≤j≤N）。

在樣本分布區(qū)間[a，b]內(nèi)用蒙特—卡洛法抽取出足夠多的均勻分布的隨機(jī)數(shù)，設(shè)有M個(gè)，其中M為足夠大的整數(shù)。這樣就能產(chǎn)生由M個(gè)數(shù)構(gòu)成的隨機(jī)數(shù)庫(kù)。在此基礎(chǔ)上，可以從產(chǎn)生的隨機(jī)數(shù)庫(kù)中抽取所需要的隨機(jī)數(shù)，使得抽取的隨機(jī)數(shù)在[a，b]范圍內(nèi)的任一小區(qū)間內(nèi)出現(xiàn)的概率與樣本出現(xiàn)的概率充分接近，具體方法為：

其中rj為在[0，1]服從均勻分布的隨機(jī)數(shù)。

對(duì)任一隨機(jī)數(shù)tj（a≤tj≤b），在規(guī)定a0=a，an=b時(shí)，若：

成立，則將落該小區(qū)間[aj-1，aj]內(nèi)的數(shù)目加1，即：

其中，k為抽取的隨機(jī)數(shù)tj落入該區(qū)間的次數(shù)，nj為第j個(gè)小區(qū)間內(nèi)的隨機(jī)數(shù)的個(gè)數(shù)。這樣，在總共所要抽取的 N隨機(jī)數(shù)中，出現(xiàn)在該區(qū)間的隨機(jī)數(shù)的概率為：

對(duì)給定的充分接近于1.0的正數(shù)ε，若：

成立，表明落入該隨機(jī)區(qū)間的隨機(jī)數(shù)已滿足要求，不再需要抽取該區(qū)間的隨機(jī)數(shù)，若上式不成立，則重復(fù)上述過(guò)程，直到每一個(gè)小區(qū)間內(nèi)的隨機(jī)數(shù)都滿足要求為止。可見(jiàn)，無(wú)論樣本服從的隨機(jī)分布多么復(fù)雜，都可以完成抽取滿足要求隨機(jī)數(shù)的過(guò)程。

2.3.2 裂隙的三維網(wǎng)絡(luò)模擬

在三維網(wǎng)絡(luò)的模型中，不連續(xù)面被模擬成圓盤(pán)面。這是因?yàn)椋鶕?jù)研究表明，不連續(xù)面在其走向方向的長(zhǎng)度和其傾向方向的長(zhǎng)度大致相等，出于數(shù)學(xué)方面的原因可簡(jiǎn)化成為圓盤(pán)[2]。圓盤(pán)由6個(gè)幾何參數(shù)來(lái)表達(dá)——坐標(biāo)系為右手坐標(biāo)系，x為正北方向，y為正東方向，圓盤(pán)的圓心坐標(biāo)，圓盤(pán)的直徑d,傾向（x軸正向順時(shí)針到圓盤(pán)面指向上的法向量在水平面的投影向量的夾角），傾角（節(jié)理圓盤(pán)的法向與z軸的夾角）。由這 6個(gè)參數(shù)可以得到下面的兩個(gè)方程：

式（5）為圓盤(pán)所在平面的空間解析幾何方程；式（6）為以點(diǎn)（x0，y0，z0）為圓心，以 d 為直徑的球體的球面方程，二者相交即為圓盤(pán)的圓周。

三維裂隙網(wǎng)絡(luò)模擬的具體過(guò)程如下：

1）確定模擬范圍及范圍內(nèi)不連續(xù)面的個(gè)數(shù)。

2）確定每條不連續(xù)面的空間位置：假定節(jié)理圓盤(pán)面服從均勻分布規(guī)則，通過(guò)蒙特卡羅模擬方法，隨機(jī)產(chǎn)生不連續(xù)面的中心點(diǎn)坐標(biāo)。

3）用蒙特卡羅法模擬生成不連續(xù)面直徑并使其服從已知的最佳概率分布。

4）用蒙特卡羅法模擬生成各不連續(xù)面產(chǎn)狀。

5）把上述各步驟模擬的結(jié)果進(jìn)行組合，從而形成一個(gè)完整的模型。

這樣，只要通過(guò)在生成域內(nèi)設(shè)置檢驗(yàn)域，每生成一條裂隙即對(duì)檢驗(yàn)域進(jìn)行一次檢驗(yàn)，直到得出的檢驗(yàn)域內(nèi)的裂隙平均面密度與實(shí)測(cè)結(jié)果充分接近來(lái)保證生成裂隙數(shù)目合理，裂隙間距的分布形式和分布參數(shù)就可滿足要求[5]。

3 有關(guān)問(wèn)題的探討

3.1 巖體各向異性廣義RQDt值計(jì)算

RQD指標(biāo)被應(yīng)用于各種巖體工程，在當(dāng)今最流行的RMR巖體工程分類系統(tǒng)和Q分類系統(tǒng)中，RQD都作為最重要的評(píng)定參數(shù)。根據(jù)Deere的定義，RQD為長(zhǎng)度大于0.1 m整段巖芯的累計(jì)長(zhǎng)度與巖芯總長(zhǎng)度之比的百分?jǐn)?shù)。人們發(fā)現(xiàn)，在鉆探中對(duì)于所有巖芯長(zhǎng)度都小于0.1 m的巖體，該方法不適合評(píng)價(jià)巖體的性質(zhì)，同樣的，如果所取得的巖芯長(zhǎng)度都大于0.1 m，RQD將為100%，這也不符合實(shí)際；通常鉆孔深度有限，巖體也是各向異性的，通過(guò)鉆孔法統(tǒng)計(jì)的RQD值只能反映局部的巖體情況，有很大的局限性。

要合理全面地利用RQD 評(píng)價(jià)巖體的質(zhì)量應(yīng)考慮巖體中不同方向及在不同界限長(zhǎng)度時(shí)RQD的變化規(guī)律。在裂隙網(wǎng)絡(luò)模擬中可以通過(guò)裂隙的間距計(jì)算出RQD值。

對(duì)于任一間距閾值，沿某一測(cè)線方向大于間距的巖體與總測(cè)線長(zhǎng)度之比的百分?jǐn)?shù)定義為廣義RQD[3]，用 RQDt表示（在本論文里，t分別取 0.1，0.2，0.3，0.4，0.5，0.6，0.7，0.8 m）。即：

式中：n代表長(zhǎng)度大于t的間距數(shù)目，為大于t的間距長(zhǎng)度；L為測(cè)線的總長(zhǎng)度。

因結(jié)構(gòu)面間距具有隨機(jī)分布特征，所以RQDt也是隨機(jī)分布的。RQDt的隨機(jī)分布規(guī)律可用間距的概率函數(shù)來(lái)表達(dá)，其數(shù)學(xué)期望E（RQDt）描述了巖石質(zhì)量指標(biāo)。根據(jù)Sen[3]建立了不同間距分布的E（RQDt）模型（結(jié)構(gòu)面的間距一般服從負(fù)指數(shù)分布或?qū)?shù)正態(tài)分布）。當(dāng)間距呈均勻分布時(shí)：

當(dāng)間距呈負(fù)指數(shù)分布時(shí)：

當(dāng)間距呈對(duì)數(shù)正態(tài)分布時(shí)：

式中：F[x]表示從－∞到x這一區(qū)間內(nèi)標(biāo)準(zhǔn)正態(tài)分布曲線下的面積；qnx為間距方差；λ代表結(jié)構(gòu)面的線密度，λ=1/（為結(jié)構(gòu)面平均間距）。

在三維裂隙網(wǎng)絡(luò)模擬圖中，在臨空面處進(jìn)行切割，獲得了臨空面的跡線模擬圖。然后再根據(jù)間距分布形式，進(jìn)行隨機(jī)抽樣產(chǎn)生各種概率分布的隨機(jī)變量。利用（8）、（9）或（10）式可以較方便地計(jì)算出不同方向和界限值所對(duì)應(yīng)的廣義RQDt值。

由圖2可以看出：在各個(gè)閥值所對(duì)應(yīng)的值中,最大的方向是NE24°，最小的方向SE98°，如表2所示。在 NE24°，SE98°方向的RQDt與t所對(duì)應(yīng)的曲線如圖3所示。

表2 NE24°方向的RQDt的均值

各塊徑的RQD在不同方位值變化很大，說(shuō)明該區(qū)結(jié)構(gòu)面的發(fā)育情況及巖體的完整性都呈現(xiàn)明顯的各向異性，這將必然導(dǎo)致巖體在不同方向上的物理、力學(xué)性質(zhì)差別較大。沿某一方向的RQD值越小，則巖體越破碎；反之，則越完整。在本研究區(qū)，巖體為塊狀結(jié)構(gòu)和層間結(jié)合較好的中厚層或厚層狀結(jié)構(gòu)，根據(jù)GB1350287-99《水利水電工程地質(zhì)勘查規(guī)范》，巖體的RQD大致范圍在0.5～0.75，在勘察設(shè)計(jì)時(shí)，NE24°方向閥值 t=0.4 m，SE98°方向閥值t=0.2 m比較合理。在RQD值最小的方向與河谷走向的交角不大，從而使得該區(qū)滑坡現(xiàn)象嚴(yán)重[2]。

4 結(jié)語(yǔ)

根據(jù)以上研究結(jié)果及研究區(qū)地質(zhì)調(diào)查情況，以下基本結(jié)論可供工程應(yīng)用時(shí)參考：

1）研究區(qū)內(nèi)普遍發(fā)育的結(jié)構(gòu)面有3～5組，其傾角較陡，一般介于40°～80°之間，傾向以110°～135°，190°～235°，300°～340°方向發(fā)育的最為集中。

2）結(jié)構(gòu)面半跡長(zhǎng)多在1.0～1.5 m之間，規(guī)模相對(duì)較小。結(jié)構(gòu)面的平均間距較小，一般為20～49 cm，總體上具有各向異性。

3）該區(qū)0.1 m塊徑時(shí)的RQD值很大，一般都達(dá)90%左右，而0.3 m塊徑時(shí)為40%～75%，0.5 m塊徑時(shí)為15%～50%左右，0.7 m塊徑時(shí)為6.5%～35%。

4）在通過(guò)對(duì)巖體出露平面的統(tǒng)計(jì)，在各個(gè)閥值所對(duì)應(yīng)的RQDt值中,最大的方向是NE24°，最小的方向SE98°。根據(jù)規(guī)范所提供的巖體質(zhì)量指標(biāo)，NE24°方向閥值 t=0.4 m，SE98°方向閥值 t=0.2 m測(cè)得的RQDt比較合理。

5）本次研究中采用了直接法和廣義RQD值等方法，取得了一定成果，可為評(píng)價(jià)研究區(qū)巖體結(jié)構(gòu)特征和工程整治提供依據(jù)，但其研究仍然是初步的，有待進(jìn)一步深入。

［1］徐鐘濟(jì).蒙特卡羅方法［M］.上海：上海科學(xué)技術(shù)出版社，1985：1-144.

［2］孫廣忠.巖體結(jié)構(gòu)力學(xué)［M］.北京：科學(xué)出版社，1988.

［3］杜時(shí)貴，何芳象，王思敬.RQD研究的幾個(gè)理論問(wèn)題［J］.現(xiàn)代地質(zhì)，1998（1）.

［4］段蔚平.邊坡巖體結(jié)構(gòu)面的模擬［J］.金屬礦山，1994（6）.

［5］魏安.巖體裂隙網(wǎng)絡(luò)的計(jì)算機(jī)模擬及其應(yīng)用［J］.西南交通大學(xué)學(xué)報(bào)，1995（4）.

［6］楊天鴻，尚富行.露天礦邊坡巖體結(jié)構(gòu)面調(diào)查及滲透特性分析［J］.勘察科學(xué)技術(shù)，1998（3）.

Simulation and application of 3D rock crack network

CHEN Zhi-jie,FENG Xi,LI Ao-song

In order to evaluate the structural characteristics of reservoir bank rock reasonably and supply the basis for project regulation,the rock crack has been investigated in the area and a 3D network simulation is studied.The study analyzes the distribution characteristics of crack occurrence by trend rose and distribution histograms of dip and dip angle based on the Monte-Carlo simulation method and obtained the dominant texture of occurrence by the crack density map,the random number by direct method and the 3D crack network chart by computer simulation conforming to reality.The study introduces the calculation of RQD2 and discusses the possible influence of its value change on reservoir bank rock stability in different directions.

rock crack;3D network;simulation study

X171.3

B

1002－0624（2011）08－0049－04

2011-05-05