高 穩(wěn)，喬 蘭，隋智力，李慶文

（1.北京科技大學(xué)土木與環(huán)境工程學(xué)院土木工程系，北京100083；2.北京科技大學(xué)金屬礦山高效開(kāi)采與安全教育部重點(diǎn)實(shí)驗(yàn)室，北京100083；3.北京城市學(xué)院，北京100083）

巖質(zhì)邊坡的穩(wěn)定性是保證大型露天鐵礦采掘作業(yè)安全的關(guān)鍵，采場(chǎng)邊坡內(nèi)部水力特性和滲流場(chǎng)的分布，對(duì)邊坡穩(wěn)定性有著重要的影響。坡體內(nèi)部的滲流特性，受巖體結(jié)構(gòu)面發(fā)育程度、巖石類(lèi)別、水頭高度等多方因素影響。土體中水力傳遞主要依靠的是土顆粒中的達(dá)西滲流作用，而巖體中的水力傳遞則由巖體節(jié)理發(fā)育程度決定。巖體內(nèi)存在的大量節(jié)理裂隙網(wǎng)絡(luò)，成為巖體中水力傳遞的主要通道［1］。在數(shù)值模擬分析中，常采用多孔連續(xù)介質(zhì)模型，將裂隙中的水等效平均到整個(gè)巖體中，再將其視為具有對(duì)稱(chēng)滲透張量的各向異性連續(xù)介質(zhì)體［2］。

水廠鐵礦是一座大型變質(zhì)巖型磁鐵礦床，礦區(qū)位于河北省遷安市境內(nèi)，西至北京200km，西南至唐山80km，東南至遷安市20km。礦區(qū)西幫與灤河河道緊鄰，灤河水系與此段邊坡內(nèi)地下水間的水力聯(lián)系密切。此外，礦區(qū)地下水還能接受大氣降水補(bǔ)給，西幫局部地段人工堆積廢石，透水性好，大氣降水就地入滲補(bǔ)給地下水，故對(duì)礦區(qū)西幫邊坡進(jìn)行滲流穩(wěn)定性分析是十分必要的。

1 水文地質(zhì)條件

北采場(chǎng)西幫有六組最優(yōu)勢(shì)節(jié)理裂隙，傾向分別為 NE60°～70°、SE120°～130°、SE160°～170°、SW200°～210°、NW340°～350°、NW270°～280°。其主要結(jié)構(gòu)面的粗糙度以平直光滑為主，結(jié)構(gòu)面的間距較大，其充填物主要是綠泥石化方解石，偶見(jiàn)輕微蝕變、局部氧化膜，巖體結(jié)構(gòu)性較好。且結(jié)構(gòu)面發(fā)育程度各個(gè)方向較為均勻，未體現(xiàn)出很強(qiáng)的方向性。由此可知，這樣產(chǎn)生的滲流場(chǎng)不會(huì)有很大的方向上的差異，對(duì)邊坡的力學(xué)特性的影響也比較均勻。根據(jù)Hoek－Brown提出的巖石力學(xué)理論，這種節(jié)理組合不會(huì)讓坡體呈現(xiàn)很明顯的各向異性，故其水力分析可以按照均質(zhì)體考慮。

根據(jù)地質(zhì)勘查報(bào)告可知，水廠鐵礦含水層特征如下所述：

1）第四系洪坡積弱～中等富水性孔隙潛水層：含水層不甚發(fā)育，厚度一般約10m左右，地下水埋深0.21～3.80m。

2）侏羅系中統(tǒng)后城組微弱富水性裂隙潛水含水巖層：巖石比較致密，節(jié)理不發(fā)育，含水微弱。

3）長(zhǎng)城系常州溝組弱～中等富水性裂隙潛水含水層：巖石堅(jiān)硬致密，形成陡峻的山脊。

4）太古界三屯營(yíng)組二段Ⅳ及Ⅵ層弱～中等富水性裂隙潛水承壓含水巖系：Ⅵ層含水層厚度19.67～79.00m，Ⅳ層含水層厚度20～62m。

5）太古界三屯營(yíng)組二段Ⅰ、Ⅱ、Ⅲ和Ⅴ層微弱～弱富水性裂隙潛水含水巖系：含水層厚度一般約80m左右。

2 采場(chǎng)邊坡三維滲流場(chǎng)分析

2.1 滲流數(shù)學(xué)模型

三維多孔連續(xù)介質(zhì)模型的數(shù)學(xué)表達(dá)式為［3］：

式中：Ss為巖體的單位儲(chǔ)存量，單位m3；kij為滲透系數(shù)，其雙腳標(biāo)的物理意義為兩次投影量。

對(duì)于三維穩(wěn)定滲流問(wèn)題，式（1）可以改寫(xiě)為［4］：

且在域Ω中：式（3）在S1上。

式（4）在S1上。

式（5）在S3和S4上。

式中：Ki（i＝x，y，z）為滲透系數(shù)；H為滲透水頭；Ω為滲流域；Z（X，Y）為浸潤(rùn)曲面或者溢出段曲面上相應(yīng)點(diǎn)（X，Y）上的位置水頭；q為邊界上單位面積流的流量空間滲流區(qū)域；S1為已知水頭值的邊界曲面；S2為給定流量曲面；S3為浸潤(rùn)面；S4為溢出段。

2.2 滲透系數(shù)的確定

2.2.1 通過(guò)單孔壓水試驗(yàn)確定滲透系數(shù)

壓水試驗(yàn)是以單位吸水量ω來(lái)表示的。在實(shí)驗(yàn)規(guī)程中，建議等效滲透系數(shù)按下式計(jì)算［5］：

式中：l為壓水試驗(yàn)段長(zhǎng)，一般為5m；r為壓水鉆孔半徑，單位m；ω為單位吸水量，單位為L(zhǎng)／（min·m2）。

2.2.2 通過(guò)注水試驗(yàn)確定滲透系數(shù)

當(dāng)鉆孔中地下水埋藏很深或者試驗(yàn)層為透水不含水層時(shí)，采用注水試驗(yàn)代替抽水試驗(yàn)，近似地測(cè)定巖層的滲透系數(shù)［6］，公式見(jiàn)式（7）。

式中：l為注水試驗(yàn)段長(zhǎng)，單位為m；r為注水鉆孔半徑，單位m；Q為注水量，單位為m3；S為水位變化量，單位m。

2.2.3 計(jì)算選用的滲透系數(shù)

現(xiàn)場(chǎng)部分壓水試驗(yàn)成果及注水試驗(yàn)成果曲線，分別如表1和圖1所示［7］。

三維滲流數(shù)值計(jì)算中所涉及到的滲透系數(shù)，如表2所示［8］。

表1 KE鉆孔壓水試驗(yàn)成果表

圖1 KE鉆孔159.60m深度注水試驗(yàn)綜合成果圖

表2 數(shù)值計(jì)算中涉及巖石的滲透系數(shù)表

2.3 數(shù)值模型的建立

根據(jù)工程巖組特征、巖體結(jié)構(gòu)特征、巖體不連續(xù)面特征、采礦設(shè)計(jì)及邊坡方位特征，水廠鐵礦北采場(chǎng)研究區(qū)共被劃分為Ⅰ～Ⅴ五個(gè)工程地質(zhì)區(qū)。本文對(duì)采場(chǎng)西幫邊坡的局部坡體進(jìn)行了三維滲流場(chǎng)計(jì)算和分析，采用FLAC3D有限差分軟件進(jìn)行建模計(jì)算。選取1剖面到7剖面之間的坡體建立計(jì)算模型，模型長(zhǎng)為2000m，寬和高均為900m，模型的開(kāi)挖底板標(biāo)高與開(kāi)挖設(shè)計(jì)深度一致，分別為－440m、－410m、－335m和－180m，模型及其在礦區(qū)的位置詳見(jiàn)圖2。

圖2 滲流計(jì)算模型及其在礦區(qū)中的位置

現(xiàn)場(chǎng)勘查資料表明，礦區(qū)巖體在滲透性上未顯示出各向異性的特點(diǎn)，故此次滲流計(jì)算采用各向同性的滲流模型進(jìn)行計(jì)算。計(jì)算中，取灤河水位按一般水位取值，為67m。

2.4 滲流模擬結(jié)果及分析

圖3給出了模型達(dá)到計(jì)算平衡后，西幫邊坡整體的三維滲流情況分布。分析圖3可知，邊坡出水點(diǎn)的平均標(biāo)高為－120m，出水點(diǎn)標(biāo)高隨坡高、坡角的減小而增高；坡體滲流場(chǎng)在遠(yuǎn)離坡面的坡體內(nèi)部受坡面狀況影響較小，基本穩(wěn)定，但在近坡面的坡體內(nèi)部影響大，變化劇烈。

分析圖4可知，剖面1的出水點(diǎn)標(biāo)高在－150m左右，水位線陡峭，與坡面呈小角度相交。

分析圖5可知，剖面3的出水點(diǎn)標(biāo)高在－110左右，水位線相對(duì)于剖面1來(lái)說(shuō)平緩了一些，與坡面呈小角度相交。

分析圖6可知，剖面5的出水點(diǎn)標(biāo)高在－100m左右，水位線相對(duì)于剖面1和剖面3來(lái)說(shuō)平緩了許多，已經(jīng)快接近水平。

圖3 西幫邊坡三維滲流分布圖

圖4 剖面1三維滲流分布圖

圖5 剖面3三維滲流分布圖

圖6 剖面5三維滲流分布圖

分析圖7可知，剖面7的出水點(diǎn)標(biāo)高在－80m左右，水位線幾乎接近水平。另外，還對(duì)剖面7豐水位時(shí)的滲流場(chǎng)進(jìn)行了計(jì)算，結(jié)果顯示此時(shí)的出水點(diǎn)標(biāo)高在－50m左右。

圖7 剖面7三維滲流分布圖

3 結(jié)論

本文在對(duì)水廠鐵礦西幫邊坡三維滲流場(chǎng)分析的基礎(chǔ)上，得出如下結(jié)論：

1）邊坡滲流場(chǎng)受邊坡高度影響較大，坡度越高，水位線越陡峭，出水點(diǎn)也越低。

2）此次計(jì)算并未考慮地表蒸發(fā)量變化對(duì)水位產(chǎn)生的影響，故整體計(jì)算結(jié)果相對(duì)保守，且未考慮斷層裂隙水的影響，不能代表局部裂隙滲水的情況。

3）剖面7相對(duì)灤河距離較近，而且豐水期出水點(diǎn)高程較高，雖然整體不會(huì)出現(xiàn)動(dòng)水破壞的情況，但是高水位對(duì)邊坡穩(wěn)定性有一定影響。建議在豐水期加強(qiáng)現(xiàn)場(chǎng)觀測(cè)，必要時(shí)可采取降排水措施。

4）開(kāi)挖對(duì)邊坡坡面附近一定范圍內(nèi)的滲流場(chǎng)影響較大，在評(píng)估邊坡穩(wěn)定時(shí)不能忽視其影響。

［1］王鵬，喬蘭，李長(zhǎng)洪，等.巖質(zhì)邊坡滲流場(chǎng)中等效連續(xù)介質(zhì)模型的應(yīng)用［J］.北京科技大學(xué)學(xué)報(bào)，2003，25（2）：99－130.

［2］張有天，周維垣.巖石高邊坡的變形與穩(wěn)定［M］.北京：中國(guó)水利水電出版社，1999.

［3］周維垣.高等巖石力學(xué)［M］.北京：水力電力出版社，1990.

［4］孫訥正.地下水流的數(shù)學(xué)模型和數(shù)值方法［M］.北京：地質(zhì)出版社，1981.

［5］東北勘測(cè)設(shè)計(jì)院修訂.SL25－92水利水電工程鉆孔壓水試驗(yàn)規(guī)程［S］.北京：中華人民共和國(guó)水利部、能源部，1992.

［6］蔡美峰，何滿潮，劉東燕.巖石力學(xué)與工程［M］.北京：科學(xué)出版社，2002.

［7］北京科技大學(xué)，首鋼礦業(yè)公司.露天礦深部安全強(qiáng)化開(kāi)采綜合技術(shù)研究專(zhuān)題研究報(bào)告［R］.2004.

［8］首鋼地質(zhì)勘查院地質(zhì)研究所.水廠鐵礦修改設(shè)計(jì)邊坡穩(wěn)定性研究補(bǔ)充地質(zhì)勘查鉆孔壓水試驗(yàn)整理成果表圖冊(cè)［R］.2010.