程 琦

(塔城地區(qū)水利水電勘察設計院，新疆 塔城 834700)

0 引 言

排水孔是水利水電工程以及巖土工程中較為常見的排水降壓、提高工程穩(wěn)定性的控制措施之一[1]。在實際工程中，將一排或多排的排水孔組成排水孔幕，從而減小水荷載對工程的影響。那么如何合理布置排水孔則是一個關鍵的問題。因為排水孔對所在區(qū)域及附近滲流場水頭分布、水力坡降等影響巨大，其附近的水力梯度相對其它區(qū)域非常大，同時排水孔的孔徑尺寸較小且數(shù)量眾多，排列也較為密集，因此無論采用哪種方法對排水孔進行模擬計算都相對來說比較復雜。如果追求計算結(jié)果的精確性，則計算中網(wǎng)絡劃分工作量將非常大，計算過程將較為復雜；但如果粗略的模擬排水孔，那么計算結(jié)果精確度將有折扣[2]。因此尋找兩者之間的平衡點，不但能精確模擬排水孔或排水孔幕的排水效果，而且還要達到簡單易行的原則將是一個難題。同時能夠精確模擬出排水效果，又是合理設計排水孔間距、深度、布置等問題的基礎。因此，對排水孔的模擬計算，國內(nèi)外研究都是非常重視。目前，無論理論基礎還是數(shù)值模擬，排水孔模擬計算都取得了較大的成績[3-4]。

利用空氣單元法對實際工程方案進行模擬計算時，必須加密排水孔附近的網(wǎng)格,從而達到所要求的精確性[5]。但是模型區(qū)域比較大以及情況復雜等原因，加密網(wǎng)格后數(shù)目將非常巨大，這對模擬計算是很大的阻礙。而且實際工程中，排水孔常常以排水孔幕的形式布置，數(shù)量較多，加密后的網(wǎng)絡數(shù)量將更加巨大，這對此方法的實際應用會帶來一定的阻力[6-7]。利用空氣單元法計算時，對于包含有壓排水設施結(jié)構(gòu)的滲流模擬計算很有幫助，不但可以減輕數(shù)據(jù)前期的處理工作量，而且具用一定的精確度。

1 排水孔模擬的空氣單元法原理

空氣單元法本質(zhì)是將排水孔視為以空氣為滲流介質(zhì)的單元體。在模擬計算中，給定合適的滲透系數(shù)，與計算區(qū)域其它巖體、土體等單元體一樣直接引入到傳統(tǒng)有限元滲流分析計算方法中。那么如何正確給定滲透系數(shù)，則是此方法模擬排水孔排水效果精確性的前提。本文利用空氣單元法得到滿足計算誤差范圍的排水孔模擬計算滲透系數(shù)取值。

在實際工程中，排水孔按自身特點主要分為兩種形式。一種是孔口出流式：水從孔內(nèi)向上流動并從上部孔口流出，此時排水孔內(nèi)部處于有壓滲流狀態(tài)，孔口各個節(jié)點水位為相應的位置高程，如地下洞室的底板排水孔，見圖1(a)。另一種是孔底出流：這種排水孔水是從孔壁滴落并從孔底流出，此時孔內(nèi)壁即為地下水溢出面，孔壁上的各個節(jié)點水頭應與各自位置高程相等，如地下洞室頂拱上的排水孔，見圖1(b)。

圖1 排水孔示意圖

在實際工程模擬計算中，排水孔中的水頭是未知的，只有個別的特殊點水頭是可以事先得到，而特殊水頭部位以及數(shù)量與排水孔的形式有關?？卓诔隽魇降目卓谒^就是孔口的位置高程，孔口以下各點水頭應大于孔口水頭。而對于孔底排水式，其孔內(nèi)壁作為地下水的溢出點，水頭應等于位置高程。由滲流基本理論可得，當水流通過滲透性突變的界面時，水流會發(fā)生折射的現(xiàn)象，這種折射現(xiàn)象是由于界面上水流的連續(xù)性引起的。

大小裂隙相互交叉處，就會出現(xiàn)所謂的偏流現(xiàn)象。這種現(xiàn)象表明大裂隙(即滲透性大)起強導水的作用，細小裂隙的水流會被匯入到大裂隙中，并導向排泄點。因此，當巖體中出現(xiàn)強導水結(jié)構(gòu)，就會產(chǎn)生排水作用。而實際工程中的排水孔就是人為布置在巖體中有強導水特性的管狀結(jié)構(gòu)。巖體和土體滲流產(chǎn)生的原因就是由于存在裂隙和孔隙，它們?yōu)闈B流提供了貯存和流動的空間?？臻g的連通以及大小決定了介質(zhì)的滲透性能。為了便于分析，可以用滲透系數(shù)K來表征滲透性能。正常情況下，排水孔就是一個充滿空氣的中空結(jié)構(gòu)，所以可將它作為一種比一般滲流介質(zhì)的滲透性大很多的特殊介質(zhì)。因此，假若能夠用一個合適的滲透系數(shù)來表征排水孔的滲透性能，并和其他實體單元一樣進行模擬計算，這種方法是完全可行的。利用這種方法進行滲流計算時，數(shù)據(jù)的準備工作也較為簡單。這種方法模擬計算結(jié)果的正確與否依賴于滲透系數(shù)K的取值合理性。

2 排水孔滲透系數(shù)的確定

在實際工程中，因為周圍介質(zhì)是不同的，因此排水孔模擬時排水孔滲透系數(shù)取值也會發(fā)生相應的變化。因此假設D等于排水孔滲透系數(shù)K1與周圍介質(zhì)滲透系數(shù)K2比值，即D=K1/K2，定義D為相對滲透系數(shù)。對其進行敏感性分析，研究D的取值對模擬計算的影響。只需確定D取值范圍，就能夠得到排水孔滲透系數(shù)在不同的工程問題中正確的取值范圍，從而精確模擬排水孔排水效果。

利用兩種完全不同的方法進行排水孔模擬計算：一種是在計算中無需模擬排水孔單元，而是在布置排水孔的位置處給定排水孔各點的水頭值，這是一種較為精確的方法，稱之為棄單元法。另一種是利用本文提出的將排水孔模擬成實體單元法，確定巖體滲透系數(shù)后，通過改變D值，使以上兩種不同的計算方法之間誤差E滿足工程模擬允許誤差范圍，以此來確定排水孔滲透系數(shù)。

3 模型計算分析

3.1 孔底出流式

3.1.1 計算模型

圖2為一地下洞室，周圍為各向同性均質(zhì)巖體。在洞室的頂拱中心含一垂直排水孔。其中寬度T=2 m，即排水孔間距為4 m。模型上下游方向即與Z軸垂直，兩個面水頭分別取H1=500 m，H2=350 m。棄單元法中排水孔壁上各節(jié)點應預先給定水頭值，水頭值與各自的位置高程相等；而利用實體單元模擬時，排水孔無需給定各點的水頭值。其余各面為不透水邊界。排水孔半徑為0.2 m，由洞室頂板向上延至巖體上表面，排水孔深40 m。假設巖體的滲透系數(shù)為K1=1×10-7m/s，而排水孔滲透系數(shù)為K2。定義D為相對滲透系數(shù)，取D=K2/K1。模型剖分時，排水孔附近網(wǎng)格加密，遠離排水孔區(qū)域，為便于計算，網(wǎng)格相對疏松。將兩種模擬方法的計算結(jié)果進行對比，通過改變D值尋找使這兩種方法的計算誤差E滿足一定要求的D的解。

圖2 計算模型及邊界示意圖

3.1.2 計算結(jié)果

以通常采用的精度較高的棄單元法為基準，選取排水孔上距孔口10 m 處的一點計算進行兩種方法的對比，見圖3。

圖3 棄單元法和空氣單元法模擬排水孔的E-lgD關系圖

當2.12.7或lgD<2.1時，計算誤差E隨著D的變化而顯著變大。因此，用空氣單元法計算孔底溢流排水孔的問題，可取2.1

3.2 孔口出流式

3.2.1 計算模型

圖2所示的各項同性三維巖體計算模型，模型內(nèi)部含一排水孔。其中寬度T=2 m，即排水孔間距為4 m。模型上下游方向即與Z軸垂直，兩個面水頭分別取H1=500 m，H2=350 m。棄單元法中孔口各節(jié)點的水位即為相應的位置高程，排水孔上表面處為自由溢流邊界；而利用實體單元模擬時，排水孔無需給定各點的水頭值。其余各面為不透水邊界。排水孔半徑為0.2 m，由洞室頂板向上延至巖體上表面,排水孔深40 m。假設巖體的滲透系數(shù)為K1=1×10-7m/s，而排水孔滲透系數(shù)為K2。定義D為相對滲透系數(shù)，取D=K2/K1。模型剖分時，排水孔附近網(wǎng)格加密，遠離排水孔區(qū)域，為便于計算，網(wǎng)格相對疏松。將兩種模擬方法的計算結(jié)果進行對比，通過改變D值尋找使這兩種方法的計算誤差E滿足一定要求的D的解。

3.2.2 計算結(jié)果分析

棄單元法和空氣單元法兩者之間的計算誤差E，是隨排水孔滲透系數(shù)的增大而減小。當lgD<3時，E隨著D的增大有明顯的減小。當lgD>2.2時，誤差可以控制在5%。

但誤差依舊隨D的增大而減小，當lgD>3時，誤差減小速率明顯減弱。當lgD在3到5之間，隨著D的增大，計算誤差E在0.1%左右振蕩。而實際工程中，認為計算誤差在5%以內(nèi)是可行的。因此，利用空氣單元法計算含排水孔的滲流問題時，可取lgD=3即排水孔的滲透系數(shù)為周圍巖體或土體滲透性的1 000倍。

4 結(jié) 論

通過空氣單元法分析，可以得到以下結(jié)論：

因為排水孔分為兩種形式，一種是有壓滲流，另一種是無壓滲流。兩種排水形式的滲透系數(shù)K值相差較大,其原因是孔底出流式的排水孔處于非飽和流狀態(tài),且水流呈跌落式流動；而孔口排水式的排水孔處于飽和流狀態(tài)。從滲透介質(zhì)的條件看,兩者的飽和度不同,流態(tài)也不同,致使K值不同。利用空氣單元法進行模擬時，兩種排水孔滿足工程要求的滲透系數(shù)。

1) 孔底出流式排水孔，即無壓滲流情況時，排水孔的滲透系數(shù)為周圍巖體(土體)滲透系數(shù)的126～500倍來模擬排水孔的作用。當lgD=2.5即取D=300時，E最小。

2) 孔口出流式排水孔，即有壓滲流情況時，排水孔的滲透系數(shù)取大于周圍巖體或土體滲透性的1 000倍。