黎照洪，胡 成，陳 剛，劉靜華

(中國(guó)地質(zhì)大學(xué)(武漢)環(huán)境學(xué)院，湖北武漢430074)

基于離散裂隙網(wǎng)絡(luò)的煙臺(tái)水封洞庫(kù)滲水點(diǎn)分析

黎照洪，胡 成，陳 剛，劉靜華

(中國(guó)地質(zhì)大學(xué)(武漢)環(huán)境學(xué)院，湖北武漢430074)

建立在致密結(jié)晶巖地區(qū)的地下水封洞庫(kù)開(kāi)挖過(guò)程中洞室及水幕巷道的滲水具有不均勻性，而現(xiàn)階段主要在開(kāi)挖后根據(jù)實(shí)際滲水量組織施工排水，若排水不及時(shí)容易造成洞室或巷道大面積浸水而形成安全風(fēng)險(xiǎn)。為在洞室開(kāi)挖前對(duì)滲水點(diǎn)分布及滲水量進(jìn)行一定預(yù)測(cè)，從而提前且合理安排排水工作，降低施工風(fēng)險(xiǎn)，以煙臺(tái)地下水封洞庫(kù)丁烷庫(kù)區(qū)為例，在對(duì)該庫(kù)區(qū)裂隙調(diào)查結(jié)果進(jìn)行一定的統(tǒng)計(jì)分區(qū)后，采用Monte-Carlo隨機(jī)模擬方法建立離散裂隙網(wǎng)絡(luò)模型，通過(guò)多次迭代確定裂隙模擬置信區(qū)間，并選取在置信區(qū)間范圍內(nèi)的裂隙模擬結(jié)果按照立方定律求取滲透系數(shù)，按區(qū)間進(jìn)行水封洞庫(kù)滲水點(diǎn)分布區(qū)劃分，同時(shí)與實(shí)測(cè)大流量滲水區(qū)域進(jìn)行比對(duì)，結(jié)果發(fā)現(xiàn)兩者具有較好的一致性，表明該方法能借助前期裂隙統(tǒng)計(jì)結(jié)果預(yù)測(cè)洞室開(kāi)挖后滲水情況，協(xié)助組織施工期洞室排水。

地下水封洞庫(kù);離散裂隙網(wǎng)絡(luò);隨機(jī)模擬;滲水預(yù)測(cè);煙臺(tái)市

地下工程施工過(guò)程中巷道滲水是影響工程進(jìn)度乃至引發(fā)事故的常見(jiàn)工程地質(zhì)問(wèn)題，而地下水封洞庫(kù)一般建于致密基巖地區(qū)地下水位以下，開(kāi)挖過(guò)程中洞室滲水問(wèn)題尤為普遍。由于裂隙介質(zhì)不均勻的導(dǎo)水性，滲水點(diǎn)的水量差異明顯，故采取一定手段合理預(yù)判地下洞庫(kù)涌水量，能幫助施工過(guò)程中合理安排排水工作或指導(dǎo)進(jìn)行其他工程措施。目前，預(yù)測(cè)地下洞室滲水量的方法主要包括水文地質(zhì)比擬法、經(jīng)驗(yàn)解析法或公式法和數(shù)值模擬法。水文地質(zhì)比擬法一般通過(guò)類比同類工程經(jīng)驗(yàn)，但目前我國(guó)已建成地下水封洞庫(kù)較少且地質(zhì)條件差異較大，適合作為類比的工程項(xiàng)目不多。經(jīng)驗(yàn)解析法以地下水動(dòng)力學(xué)為基礎(chǔ)，通過(guò)概化庫(kù)區(qū)的地下水滲流條件而對(duì)洞室中的涌水量進(jìn)行預(yù)測(cè)，如大島揚(yáng)志方法、佐藤邦明方法、我國(guó)鐵道工程關(guān)于洞室初期涌水量和穩(wěn)定涌水量的計(jì)算公式等。數(shù)值模擬法則根據(jù)概念模型對(duì)開(kāi)挖區(qū)地下水滲流場(chǎng)進(jìn)行計(jì)算機(jī)模擬，如王者超等［1］、何國(guó)富［2］綜合運(yùn)用多種方法分析了地下水封洞庫(kù)工程的涌水量;龐家偉等［3］運(yùn)用數(shù)值模擬手段分析煙臺(tái)地下水封洞庫(kù)在施工及運(yùn)營(yíng)期的涌水量;韓曼［4］分析了有無(wú)水幕條件下地下水封洞庫(kù)涌水量的差異。上述三種方法在粗略估算地下水滲水總量方面起到了良好效果，但對(duì)于巖體各向異性方面考慮不足，不能較好地預(yù)測(cè)滲水點(diǎn)具體分布區(qū)域。為了對(duì)滲水點(diǎn)具體區(qū)域進(jìn)行預(yù)測(cè)，需要對(duì)地下水封洞庫(kù)地下滲透系數(shù)張量進(jìn)行評(píng)估，如張龍等［5］引入試井理論，介紹了地下水封洞庫(kù)水平水幕鉆孔滲透張量的求解方法;季惠彬［6］則根據(jù)Hyung等［7］總結(jié)的地下水封洞庫(kù)巖體滲透系數(shù)張量的合理范圍，評(píng)價(jià)了黃島水封洞庫(kù)巖體滲透系數(shù)張量。此外，通過(guò)建立離散裂隙網(wǎng)絡(luò)的方法，可有效研究地下裂隙的不均勻分布，并借用立方定律，可求取巖體滲透系數(shù)［8-13］。但目前在地下水封洞庫(kù)領(lǐng)域應(yīng)用較少，僅有馬峰［14］根據(jù)黃島水封洞庫(kù)裂隙統(tǒng)計(jì)結(jié)果求解了裂隙巖體滲透系數(shù)張量，并分析了滲透系數(shù)張量的整體特點(diǎn)，尚缺乏根據(jù)滲透系數(shù)與實(shí)際滲水量進(jìn)行對(duì)比的研究。

本文以煙臺(tái)地下水封洞庫(kù)丁烷庫(kù)區(qū)為例，基于裂隙調(diào)查結(jié)果，按裂隙密度對(duì)研究區(qū)進(jìn)行劃分后，采用模擬軟件Fracman的Monte-Carlo隨機(jī)模擬方法生成離散裂隙網(wǎng)絡(luò)，通過(guò)多次迭代確定裂隙模擬置信區(qū)間，并選取在置信區(qū)間范圍內(nèi)的裂隙模擬結(jié)果按照立方定律求取裂隙巖體滲透系數(shù)，按區(qū)間進(jìn)行水封沿庫(kù)滲水點(diǎn)分布區(qū)統(tǒng)計(jì)劃分，同時(shí)與實(shí)測(cè)滲水點(diǎn)進(jìn)行對(duì)比，結(jié)果表明該方法能較好地預(yù)測(cè)洞庫(kù)滲水區(qū)域及滲水特點(diǎn)。

1 煙臺(tái)水封洞庫(kù)工程概況

工程項(xiàng)目為煙臺(tái)某公司商業(yè)儲(chǔ)備區(qū)，分為丁烷洞庫(kù)庫(kù)區(qū)、LPG洞庫(kù)庫(kù)區(qū)和丙烷洞庫(kù)庫(kù)區(qū)，其中丁烷洞庫(kù)庫(kù)區(qū)為本文的研究對(duì)象。該煙臺(tái)地下水封洞庫(kù)丁烷庫(kù)區(qū)截面呈馬蹄形，拱頂跨度約18 m，高度約26 m，包含3條水幕巷道、1個(gè)豎井和3個(gè)主洞室，分別編號(hào)為BC1、BC2、BC3，主洞室之間配有連接巷道，水幕巷道在洞室上方20 m，總長(zhǎng)度784 m。庫(kù)址區(qū)地形地貌上屬山間河谷沖洪積丘陵到平原地貌過(guò)渡帶，基巖為燕山早期的黑云母二長(zhǎng)花崗巖，主要巖脈為閃長(zhǎng)玢巖脈、煌斑巖脈、燕山期花崗斑巖脈等，庫(kù)址西側(cè)靠近丁烷庫(kù)區(qū)發(fā)育一局部斷裂段。

2 離散裂隙網(wǎng)絡(luò)的建立

2．1 裂隙測(cè)量統(tǒng)計(jì)

本次通過(guò)8組鉆孔電視成像、2組已開(kāi)挖洞室地質(zhì)素描，對(duì)研究區(qū)結(jié)構(gòu)面裂隙進(jìn)行了測(cè)量，得到水幕巷道裂隙314條、水幕巷道下方儲(chǔ)藏品主洞室裂隙 388條，裂隙優(yōu)勢(shì)產(chǎn)狀為:NE62°～74°∠79°～89°、SE111°～118°∠78°～84°、SW242°～253°∠83°～89°、NW289°～305°∠66°～79°。

2．2 模擬區(qū)劃分

將洞室施工區(qū)域適當(dāng)向外擴(kuò)展，選取長(zhǎng)為500 m、寬為250 m、高為200 m的空間區(qū)域作為模擬區(qū)，如圖1所示。

圖1 模擬區(qū)范圍平面圖Fig．1 Plan of simulation area

根據(jù)工程重要性在水平方向劃分為中部洞室主體區(qū)域以及左、右兩側(cè)外圍區(qū)域共3部分，并根據(jù)裂隙密度隨深度變化特征，將每個(gè)部分在垂直方向上各劃分為3個(gè)裂隙等密度區(qū)域，共計(jì)9個(gè)分區(qū)，具體分區(qū)位置見(jiàn)圖2，其中V區(qū)為洞庫(kù)工程的水幕巷道及油氣品洞室所在區(qū)域。各分區(qū)空間位置、裂隙密度及優(yōu)勢(shì)產(chǎn)狀統(tǒng)計(jì)見(jiàn)表1。

2．3 裂隙結(jié)構(gòu)面幾何要素統(tǒng)計(jì)及裂隙網(wǎng)絡(luò)生成

裂隙結(jié)構(gòu)面幾何要素統(tǒng)計(jì)包括產(chǎn)狀統(tǒng)計(jì)、裂隙跡長(zhǎng)統(tǒng)計(jì)、裂隙隙寬統(tǒng)計(jì)等，本次研究對(duì)洞室主體部分Ⅳ、Ⅴ、Ⅵ區(qū)裂隙結(jié)構(gòu)面幾何要素進(jìn)行了統(tǒng)計(jì)，結(jié)果顯示:裂隙跡長(zhǎng)服從對(duì)數(shù)正態(tài)分布;裂隙隙寬服從負(fù)指數(shù)分布;Ⅳ、Ⅴ區(qū)產(chǎn)狀分布類型為二元賓漢類型，Ⅵ區(qū)產(chǎn)狀分布類型為費(fèi)雪類型。

圖2 模擬區(qū)分區(qū)示意圖Fig．2 Zoning of the simulation area

表1 離散裂隙模擬分區(qū)情況統(tǒng)計(jì)表Table 1 Statistics of the simulation zone of the discrete fracture

在裂隙結(jié)構(gòu)面幾何要素統(tǒng)計(jì)結(jié)果基礎(chǔ)之上，采用Monte-Carlo法隨機(jī)生成各個(gè)分區(qū)內(nèi)的裂隙，研究區(qū)裂隙結(jié)構(gòu)面網(wǎng)絡(luò)模擬結(jié)果見(jiàn)圖3。

圖3 研究區(qū)離散裂隙網(wǎng)絡(luò)Fig．3 Discrete fracture network of the study area

圖3中每一種顏色的片狀六邊形分別表示子分區(qū)中某一優(yōu)勢(shì)產(chǎn)狀控制的裂隙集，六邊形的分布位置及密集程度表示裂隙集的空間發(fā)育密集位置及其密集程度。由圖3可見(jiàn)，垂向上(沿Z軸正方向)，自上而下裂隙發(fā)育程度逐步降低;水平方向上，自西向東(即X軸正方向)裂隙發(fā)育密集程度下降，裂隙發(fā)育特點(diǎn)受研究區(qū)西側(cè)的局部斷裂帶影響明顯。

3 滲透系數(shù)求取與滲水區(qū)域?qū)Ρ?/h2>

3．1 置信區(qū)間分析與滲透系數(shù)求取

為合理預(yù)測(cè)地下封洞庫(kù)滲水情況，需從模擬結(jié)果中選取與洞室開(kāi)挖空間相連通的裂隙。雖然Monte-Carlo方法［15］是一種統(tǒng)計(jì)意義上概率分布的隨機(jī)模擬方法，不同模擬結(jié)果間存在一定程度差異，但隨著模擬次數(shù)增加能表現(xiàn)其宏觀規(guī)律性，故利用Fracman軟件的宏命令功能產(chǎn)生100次連通性分析隨機(jī)模擬結(jié)果，并分別統(tǒng)計(jì)每次所得出的與洞室連通的裂隙數(shù)量，最終得到丁烷洞庫(kù)各個(gè)子洞室(BC1、BC2、BC3)的連通裂隙數(shù)量分布，詳見(jiàn)圖4至圖6。

圖4 子洞室BC1連通裂隙數(shù)量迭代統(tǒng)計(jì)結(jié)果Fig．4 Iterative statistical result of the numbers of fractures connecting to the sub-cavern BC1

圖5 子洞室BC2連通裂隙數(shù)量迭代統(tǒng)計(jì)結(jié)果Fig．5 Iterative statistical result of the numbers of fractures connecting to the sub-cavern BC2

本文利用SPSS軟件對(duì)裂隙數(shù)據(jù)進(jìn)行分析，采用95%置信率下的連通裂隙數(shù)數(shù)量范圍作為置信區(qū)間，用于滲透性分析。各子洞室對(duì)應(yīng)的連通裂隙數(shù)量置信區(qū)間見(jiàn)表2。

圖6 子洞室BC3連通裂隙數(shù)量迭代統(tǒng)計(jì)結(jié)果Fig．6 Iterative statistical result of the numbers of fractures connecting to the sub-cavern BC3

同樣，為保證裂隙模擬結(jié)果的可靠性，進(jìn)行了多次模擬，得到一個(gè)在置信區(qū)間范圍內(nèi)的連通裂隙數(shù)量模擬結(jié)果值，再根據(jù)立方定律求取對(duì)應(yīng)的滲透系數(shù)，并根據(jù)滲透系數(shù)數(shù)量級(jí)進(jìn)行分區(qū)間裂隙數(shù)量統(tǒng)計(jì)，其結(jié)果見(jiàn)表3。

表2 子洞室連通裂隙數(shù)量置信區(qū)間Table 2 Confidence interval of the numbers of fractures connecting to the sub-caverns

表3 子洞室連通裂隙數(shù)量模擬結(jié)果及對(duì)應(yīng)的滲透系數(shù)Table 3 Simulation result of the numbers of fractures connecting to the sub-carvens and the corresponding permeability coefficients

3．2 模擬結(jié)果與實(shí)測(cè)滲水區(qū)域?qū)Ρ?/p>

將連通裂隙數(shù)量多(取連通裂隙數(shù)量大于10條的部分)且滲透性相對(duì)較大的部分作為預(yù)測(cè)大流量滲水區(qū)，即BC1子洞室的K090～130 m及K0240～270 m區(qū)域、BC2子洞室的K0150～200 m區(qū)域以及BC3子洞室的K0150～230 m區(qū)域，并在洞室開(kāi)挖完成后對(duì)現(xiàn)場(chǎng)實(shí)測(cè)滲水點(diǎn)進(jìn)行統(tǒng)計(jì)，劃分實(shí)測(cè)大流量滲水區(qū)域，對(duì)3個(gè)子洞室共劃分出6個(gè)實(shí)測(cè)大流量滲水區(qū)域，詳見(jiàn)表4。

表4 丁烷洞庫(kù)各子洞室大流量滲水區(qū)域Table 4 High-flux seepage zones in each sub-carven of the butane reservoir

在丁烷洞庫(kù)平面圖上將實(shí)測(cè)和預(yù)測(cè)大流量滲水區(qū)域這兩部分區(qū)域進(jìn)行投影，可得到該丁烷洞庫(kù)實(shí)測(cè)與預(yù)測(cè)大流量滲水區(qū)域的對(duì)比圖，見(jiàn)圖7。從實(shí)測(cè)滲水區(qū)域與預(yù)測(cè)滲水區(qū)域在平面圖上的重疊情況來(lái)看，6個(gè)實(shí)測(cè)滲水區(qū)域中，有5個(gè)實(shí)現(xiàn)了較準(zhǔn)確的預(yù)測(cè)，其預(yù)測(cè)效果良好。

圖7 丁烷洞庫(kù)實(shí)測(cè)與預(yù)測(cè)大流量滲水區(qū)域?qū)Ρ葓DFig．7 Comparison of the measured and predicted high-flux seepage zone in butane reservoir

4 結(jié) 論

本文通過(guò)對(duì)煙臺(tái)地下水封洞庫(kù)丁烷庫(kù)區(qū)開(kāi)展離散裂隙網(wǎng)絡(luò)構(gòu)建工作與現(xiàn)場(chǎng)滲水統(tǒng)計(jì)工作，對(duì)比了通過(guò)數(shù)值模擬得到的滲透系數(shù)偏大區(qū)域與實(shí)測(cè)大流量滲水區(qū)域的分布差異，發(fā)現(xiàn)兩者具有較好的一致性。因此，在洞室未開(kāi)挖階段，通過(guò)收集整理前期的裂隙展布數(shù)據(jù)，建立離散裂隙網(wǎng)絡(luò)模型，在一定程度上可以實(shí)現(xiàn)對(duì)開(kāi)挖階段大流量滲水區(qū)域的預(yù)測(cè)，以便協(xié)助組織施工期洞室排水。

［1］王者超，李術(shù)才，梁建毅，等．地下水封石油洞庫(kù)滲水量預(yù)測(cè)與統(tǒng)計(jì)［J］．巖土工程學(xué)報(bào)，2014，36(8):1490-1497．

［2］何國(guó)富．湛江某地下水封洞庫(kù)涌水量估算與分析［J］．西部探礦工程，2011，23(11):19-24．

［3］龐家偉，丁國(guó)平，陳剛，等．地下水封洞庫(kù)水位降深與涌水量的數(shù)值分析［J］．油氣儲(chǔ)運(yùn)，2013，32(9):1007-1011．

［4］韓曼．地下水封石油洞庫(kù)滲流場(chǎng)及溶質(zhì)運(yùn)移模擬研究［D］．青島:中國(guó)海洋大學(xué)，2007．

［5］張龍，陳剛，胡成．試井理論在洞庫(kù)水幕效率試驗(yàn)中的應(yīng)用［J］．安全與環(huán)境工程，2015，22(1):7-10．

［6］季惠彬．水封地質(zhì)條件對(duì)地下水封洞庫(kù)選址的影響研究［J］．安全與環(huán)境工程，2013，20(2):138-141．

［7］Hyung S Y，Jae G K，Kyung S K，et al．Groundwater flow characterization in the vicinity of the underground caverns in fractured rockmasses by numerical modeling［J］．Geosciences Journal，2004，8(4): 401-413．

［8］賀少輝，廖國(guó)華．裂隙巖體裂隙網(wǎng)絡(luò)滲流模型研究［J］．礦冶工程，1997，17(2):11-14．

［9］王培濤，楊天鴻，于慶磊，等．基于離散裂隙網(wǎng)絡(luò)模型的節(jié)理巖體滲透張量及特性分析［J］．巖土力學(xué)，2013，34(2):448-455．

［10］Miller S M．A statistical method to evaluate homogeneity of structural populations［J］．Mathematical Geology，1983，15(2):317-328．

［11］Mahtab M A，Yegulalp T M．A similarity test for grouping Orientation data in rock mechanics［C］//Proceedings of the25th U．S．Symposium on Rock Mechanics．New York:American Institute of Mining，Metallurgical，and Petroleum Engineers，1984:495-502．

［12］孫蓉琳，梁杏，靳孟貴．裂隙巖體滲透系數(shù)確定方法綜述［J］．水文地質(zhì)與工程地質(zhì)，2006(6):120-123．

［13］任文明，梁久正．大型地下洞庫(kù)巖體力學(xué)參數(shù)取值及工程應(yīng)用［J］．石油工程建設(shè)，2012，38(1):14-15．

［14］馬峰．黃島地下水封洞庫(kù)裂隙巖體滲透性研究［D］．武漢:中國(guó)地質(zhì)大學(xué)武漢，2010．

［15］Oda M．Permeability tensor for discontinuous rock asses［J］．Geotechnique，1984，35:483-495．

Seepage Analysis of the Water-sealed Cavern in Yantai City Based on the Discrete Fracture Network Model

LI Zhaohong，HU Cheng，CHEN Gang，LIU Jinghua
(School of Environmental Studies，China University of Geosciences，Wuhan430074，China)

Seepage inhomogeneity is common characteristics in carvens and water curtain tunnels during excavation which are built in compact crystalline rock area．At presentthe site drainage is deployed after the excavation has begun according to the actual seepage quantity．Once the water is not discharged in time，the caves or the tunnels tend to immerse to a large extent with a hysteretic drainage，which may cause safety risk．To forecast the seepage point and its quantity before the caves are dug，so as to arrange reasonably the drainage in advance and lower the constructive risk，this paper calculates and classifies the fracture in the water-sealed cavern of Yantai city and builds the discrete fracture network model based on the stochastic method by Monte-Carlo．With multiple iterations，the paper confirms the confidence interval of the fracture simulation，then calculates the permeability coefficients using cubic law according to the simulated result of fractures which lies in confidence interval．Furthermore，the paper divides seepage points of water-sealed cavern into several parts by the confidence interval and compares with the virtual parts．The results exhibit fine consistency between the simulated and the measured parts，showing that the method can predict the seepage after caving，which helps in site drainage in construction period．

underground water-sealed cavern;discrete fracture network;stochastic simulation;seepage prediction; Yantai City

X948;TU926

ADOI:10．13578/j．cnki．issn．1671-1556．2016．05．029

1671-1556(2016)05-0170-04

胡 成(1976—)，男，博士，副教授，主要從事3S技術(shù)在水文地質(zhì)中的應(yīng)用以及工程水文地質(zhì)方面的研究。E-mail:hu-cheng @cug．edu．cn

2015-11-09

2016-08-24

黎照洪(1990—)，男，碩士研究生，主要研究方向?yàn)楣こ趟牡刭|(zhì)。E-mail:526152304@qq．com