導語：裂隙巖體滲流管理論文一文來源于網友上傳，不代表本站觀點，若需要原創文章可咨詢客服老師，歡迎參考。

1、引言

隨著人類工程和開采活動的日益增加，其影響越來越多地涉及到地下裂隙巖體，并且人類工程和開采活動與裂隙巖體之間的相互作用有著日益加劇的趨勢，這使得人們對裂隙巖體的工程及開采特性備加關注。隨著工程實踐和科學研究的深入，人們已逐漸認識到裂隙巖體所賦存地質環境的復雜性及其所誘發地質災害的多變性，于是劉繼山、仵彥卿、柴軍瑞、黃濤等國內一些專家學者在不同程度、不同角度對裂隙巖體賦存地質環境中各個因素之間相互影響作用的課題進行了有意義的探索和研究。裂隙巖體中存有地下水，地下水在巖體中會產生滲流，在復雜的地質環境中，不同的因素之間隨著時間、空間而發生著復雜的動態變化，其中裂隙巖體中滲流的熱學效應（滲流對溫度的影響機理）是一個比較重要的影響方面。

2、裂隙巖體地下水滲流的基本理論

對于裂隙巖體滲流的熱學效應的分析，首先應先了解巖體滲流的基本理論。一般情況下，巖體的各個裂隙中均含有水，而地下水在巖體中會產生滲流，所謂滲流是指含空隙（孔隙、裂隙等）介質中流體（液體、氣體）通過空隙的流動[1]。地下水的流動是最典型的滲流現象，裂隙巖體中地下水的滲流特性體現了滲流理論的主體特征。

2.1達西（Darcy）定律

達西定律的最初表達式[2]（Darcy，1855年提出）：

式中，為滲流速度，K為滲透系數，J為水力坡降,為滲透水頭H沿滲流方向的梯度。上式表明，滲流速度與水力坡降成線性關系，因此達西定律也被稱為現行滲透規律。

假設巖體中的滲流以裂隙滲流為主，忽略其間巖塊的滲透作用，那么裂隙的分布就相對比較密集，表征單元體比較小，視裂隙巖體為等效連續介質，用連續介質的方法描述巖體的滲流問題，并考慮巖體的非均質性和各向異性，以滲透系數張量來描述巖體裂隙的滲透性能。

本人根據矢量及張量的定義和連續性方程將達西定律進行了推廣，可得下式：

式中，為滲流速度矢量，為介質的滲流系數張量，為水力坡降矢量。上式的展開式為：

==－

2.2費克（FICK）定律

對溶質或壓縮性較大的氣體在介質空隙中的運動，可采用張蔚榛(1996年)關于費克定律[3]公式：

=

式中，為擴散通量或氣體流量,為擴散系數,為濃度(或密度)沿擴散方向的梯度。

2.3立方定律

單直平滑裂隙中的滲流滿足立方定律[4]，為了更好的描述裂隙介質的滲流規律，本文引用了田開銘、萬力（1989年）著作中的立方定律的公式：

==

式中，為通過單直平滑裂隙斷面的單寬流量，γ和分別為流體的容量和動力粘滯系數，

為裂隙的隙寬，為沿裂隙方向的水力坡降，=為裂隙的等效滲透系數，為滲流水頭沿裂隙方向的梯度，=－。

對于實際的粗糙不規則的裂隙及充填裂隙，須對立方定律進行修正[5]。

3、巖體溫度場分布

通常把人類工程活動所涉及的巖體內的溫度分布和熱狀態稱為巖體溫度場，巖體溫度場的分布一方面受制于地殼淺層溫度分布和熱狀態的控制，；另一方面受巖體工程和外界條件(包括地下水滲流)的影響。理論和實驗研究表明，熱力的傳遞（即傳熱）一般有三種方式：傳導、對流和輻射。工程巖體中溫度場的分布是以傳導和對流為主實現的，而對流則是以地下水的滲流運動為基礎進行的。

為了使巖體溫度場分布的描述更加具體化、形象化，本人根據熱量平衡原理，把傳導型傳熱方程、因地下水參與引起的對流項及通過地下熱源區因輻射增加的內熱源項進行了重新組合，引用了黃濤等（1999年）關于巖體溫度場的分布模型[6]，如下式：

Ca+cvi=+c(－)q(1)

式中，為巖體的熱動力彌散系數張量；xi、xj為笛卡爾坐標，,=1，2，3并遵循求和約定；為時間；為巖石的熱容量；c為地下水的熱容量；vi為滲透速度；q為單位體積含水巖體中源匯流量；為地下水溫;為源匯項溫度。

4、裂隙巖體中滲流的熱學效應分析

裂隙巖體滲流的熱學效應即巖體中滲流對溫度的影響機理。當裂隙巖體中有滲流發生時，一方面，地下水的滲流運動促成了巖體與地下水體之間發生熱傳導的熱量傳遞與交換。另一方面，地下水作為裂隙巖體中熱量交換的載體，通過地下水本身的滲流運動而產生熱對流的熱能轉移過程。以上兩個方面的熱能轉移，為達到地下水體與裂隙巖體二者溫度分布的均衡，最終導致裂隙巖體溫度場分布的總體衰減，即巖體溫度值的總體降低。

在一向導熱的情況下，當裂隙巖體中存在滲流時，熱量包括兩部分：一部分是由于巖體本身的熱傳導作用，等于－；另一部分是由滲流夾帶的熱量，等于，因此在總結朱伯芳等（1976年）關于熱傳導原理的基礎上柴軍瑞（2001年）在其著作中總結了熱流量的公式為[7][8]：

－(2)

式中，為沿一維坐標軸x方向的熱流量;為水的比熱；為水的密度；為巖體介質的導熱系數；為滲流速度；為溫度。

本人對該公式進行了分析，認為吸熱的熱源可作為負的熱源，將公式兩邊取負號，并兩邊同時作微分處理，得到了在單位時間內流入單位體積的凈熱量為：

－(3)

根據熱量平衡原理，這個熱量必須等于單位時間內裂隙巖體介質溫度升高所吸收的熱量，而單位時間內裂隙巖體介質溫度升高所吸收的熱量，故，即：

(4)

式中，為巖體介質的比熱，為巖體介質的密度。

為了更全面的對滲流影響下的溫度場進行分析，本人對式（4）進行了推廣：

如果溫度只隨平面內的兩個位置坐標而變，即溫度沿z方向是常數，，則溫度場是兩向的，推得滲流影響下的二維導熱方程：

(5)

同理，三向導熱的情況下，在不考慮源（匯）項的情況下，可推得考慮滲流影響的三維導熱方程：

(6)

由（6）式可以看出，溫度場的分布與滲流速度場的分布有密切的關系；滲流速度越大，對溫度場的影響也就越大。而滲流速度場的分布又由滲流場水頭的分布決定，即。所以式（6）定量地反映了滲流對溫度的影響機理。

5、結束語

通過對地下裂隙巖體滲流的熱學效應的分析可以看出，地下水在巖體中會產生滲流，地下水的滲流運動直接影響了巖體與地下水體之間發生熱傳導的熱量傳遞與交換，同時地下水本身的滲流運動通過熱對流的方式進行了熱能轉移，這些因素影響了裂隙巖體的溫度場分布，而且地下水在裂隙巖體中的滲流運動大大降低了巖體中的溫度分布值，地下水在裂隙巖體中的滲流速度的大小直接控制了巖體溫度的變化幅度，因此，滲流對溫度的影響是非常明顯的，對裂隙巖體滲流的熱學效應的分析研究具有很實際的意義，這還有利于加深對滲流影響溫度分布的因素的認識。

參考文獻

[1]朱詩鰲，壩工技術史，北京：水利電力出版社，1995

[2]祁慶和，水工建筑物，北京：水利水電出版社，1986

[3]張蔚榛，地下水與土壤水動力學，北京：中國水利電力出版社，1996

[4]田開銘、萬力，各向異性裂隙介質滲流性的研究與評價，北京：學苑出版社，1989

[5]忤彥卿、李俊亭、張倬元，地下水動態觀測網的優化設計，成都：成都科技大學出版社，1993

[6]黃濤、楊立中，隧道裂隙巖體溫度－滲流耦合數學模型研究，巖土工程學報，1999，21（5）：554～558

[7]朱伯芳、王同生、丁寶瑛、郭之章，水工混凝土的溫度應力與溫度控制，北京：水利電力出版社，1976

[8]柴軍瑞，大壩工程滲流力學，西藏：西藏人民出版社，2001