， ，

(1．長江科學(xué)院 a．工程安全與災(zāi)害防治研究所；b．水利部水工程安全與病害防治工程技術(shù)研究中心；c．國家大壩安全工程技術(shù)研究中心，武漢 430010； 2.中南電力設(shè)計(jì)院 發(fā)電公司，武漢 430071)

1 研究背景

土石壩滲漏探測與滲流監(jiān)測研究是水利工程界長期以來的重要研究課題。滲流熱監(jiān)測技術(shù)是一種通過監(jiān)測大壩溫度分布來間接查找滲漏點(diǎn)、獲得滲流場分布、反演滲流場參數(shù)的方法，被稱為溫度示蹤滲流監(jiān)測技術(shù)。利用測壓管、滲壓計(jì)等點(diǎn)式儀器的傳統(tǒng)滲流監(jiān)測方法主要存在以下缺陷：①精度不高，分辨率較低；②點(diǎn)式測量，漏檢率高；③儀器的長期穩(wěn)定性差；④大量埋設(shè)儀器，成本較高。滲流熱監(jiān)測技術(shù)漏檢率較低，能較全面反映工程滲流活動(dòng)的狀態(tài)與防滲效果。較其他的新型滲流監(jiān)測技術(shù)，如同位素示蹤法、高密度電法、地質(zhì)雷達(dá)法、水質(zhì)分析法等，滲流熱監(jiān)測技術(shù)具有對環(huán)境無污染、對集中滲漏敏感度較高、對經(jīng)費(fèi)與人力資源要求較低的優(yōu)勢[1]。隨著分布式光纖測溫技術(shù)在土石壩溫度監(jiān)測中的應(yīng)用，人們可以獲得大壩溫度場時(shí)空連續(xù)分布，更準(zhǔn)確、有效地了解壩體內(nèi)部的滲流狀態(tài)，有助于加深對滲流狀態(tài)的發(fā)展變化過程的認(rèn)識(shí)。

2 滲流熱監(jiān)測技術(shù)理論發(fā)展過程

滲流熱監(jiān)測技術(shù)的理論基礎(chǔ)是土體中滲流場與溫度場的互相影響。對滲流熱監(jiān)測技術(shù)的理論研究經(jīng)歷了定性分析至定量分析的過程。

2.1 定性分析階段

滲流熱監(jiān)測技術(shù)起步于國外，從地下水地溫研究獲得啟發(fā)。1965年，美國加州Occidental大學(xué)地質(zhì)系的Joesph H.Birman將這一技術(shù)應(yīng)用于大壩的漏水檢查中，結(jié)果發(fā)現(xiàn)壩體低溫區(qū)正是大壩下游觀察到的滲漏區(qū)域，并由抽水實(shí)驗(yàn)得到了滲透系數(shù)和溫度測值的負(fù)相關(guān)關(guān)系[2]。隨后美國墾務(wù)局將此技術(shù)應(yīng)用于一些病險(xiǎn)土石壩的處理。蘇聯(lián)等國也開始采用溫度場探測大壩滲漏，通過溫度異常發(fā)現(xiàn)壩體、壩肩、壩址的集中滲漏通道[3]。定性分析階段的研究[4-5]基于觀測到大壩壩體溫度異常特征，通過比較大壩溫度曲線與氣溫曲線的差異，來判斷可能存在的滲漏位置，或是比較滲漏水流速變化、水頭變化與溫度測值曲線，建立滲漏水流速、水頭變化與溫度變化的定性關(guān)系。其目的是判斷滲漏點(diǎn)位置或?qū)α鲌鰠?shù)作定性分析，常作為其他探漏方法的輔助手段，或與其他物探結(jié)果作對比驗(yàn)證。定性分析階段國內(nèi)的相關(guān)研究較少，文獻(xiàn)[6]利用丹江口水利樞紐的溫度觀測資料，將壩基溫度場等溫線與等勢線進(jìn)行對比，以說明滲流場和溫度場具有變位一致性，具有較強(qiáng)的代表性。

2.2 定量分析階段

定量分析分為正演分析與反演分析。正演分析根據(jù)邊界條件與滲流類型不同，可分為集中滲漏研究和均勻滲漏研究2類。

集中滲漏問題研究多采用解析方法分析，這是由于集中滲漏研究模型多基于理想狀態(tài)下建立，可以求得精確的解析解[7]，而均勻滲流多采用數(shù)值模擬方法計(jì)算。集中滲漏問題研究的模型基礎(chǔ)是“虛擬熱源法”[8]，該方法將集中滲漏通道看作一個(gè)虛擬的熱源，滲漏通道外溫度分布僅與滲流帶最外層有關(guān)，滲漏通道內(nèi)部溫度由滲流水控制，外部溫度則由熱傳導(dǎo)支配，通過熱傳導(dǎo)方程推導(dǎo)溫度在滲漏通道外的分布，最后由溫度分布確定滲漏范圍。文獻(xiàn)[9-11]分別在線狀、面狀和圓柱狀熱源形狀假設(shè)條件下，或是根據(jù)疊加原理[12]，利用能量守恒定律和熱傳理論建立了管涌滲漏通道的持續(xù)線熱源模型，計(jì)算了集中滲流條件下的溫度場分布。各種假設(shè)的邊界形狀與實(shí)際滲流存在一定的差別，對復(fù)雜邊界條件的問題難以求出解析解。

滲流熱監(jiān)測技術(shù)需要將土石壩溫度場與滲流場計(jì)算結(jié)合起來。滲流熱監(jiān)測定量分析要真實(shí)模擬滲流場和溫度場互相影響規(guī)律，需建立包含影響項(xiàng)的溫度場、滲流場計(jì)算模型。朱伯芳[13]以考慮滲流的一維導(dǎo)熱方程的解析解為基礎(chǔ)，分析了壩體滲流隨壩高、水頭與滲透系數(shù)的變化對壩體穩(wěn)定溫度場的影響。由于可以求出解析解的均為邊界條件、定解條件比較簡單的情況，基于有限單元法、有限差分法等數(shù)值方法建立的模型更具有適用性。賴遠(yuǎn)明等[14]針對寒區(qū)大壩凍融問題，給出了帶相變溫度場和滲流場耦合模型，模型中的耦合參數(shù)是巖土體的比熱容與導(dǎo)熱系數(shù)，滲透系數(shù)也與含水率(未凍水率)有關(guān)，并推導(dǎo)了有限元計(jì)算弱形式。柴軍瑞等[15-16]考慮水的運(yùn)動(dòng)黏度的溫變效應(yīng)與溫度梯度勢流動(dòng)以及水的對流換熱效應(yīng)，建立并計(jì)算了混凝土壩的溫度場與滲流場耦合分析模型，此模型同樣適用于土石壩計(jì)算。陳建余等[17-19]考慮滲流場的熱效應(yīng)，對混凝土壩或碾壓混凝土壩的穩(wěn)定或非穩(wěn)定滲流問題建立了二維有限元耦合模型，定量分析了滲透場對溫度場的影響程度。崔皓東和朱岳明[20]對大壩蓄水初期的非穩(wěn)定滲流場建立了三維滲流場與溫度場耦合模型并進(jìn)行計(jì)算，結(jié)果表明混凝土壩滲流場溫度效應(yīng)明顯，滲透系數(shù)受溫度場直接影響。李端有等[21]對土石壩熱力學(xué)機(jī)理做了比較詳盡的分析，提出利用二維溫度場與滲流場耦合方程進(jìn)行雙場迭代求解的方法。孫曼[22]利用滲流場與溫度場在微分方程上的相似性，利用ANSYS軟件的溫度模塊來模擬滲流場，先計(jì)算單物理場，再利用中間傳遞變量進(jìn)行間接耦合，并與甘油試驗(yàn)成果進(jìn)行了比較，驗(yàn)證了計(jì)算的正確性。熊健[23]分析了水的物理化學(xué)參數(shù)的溫變效應(yīng)，分析了均質(zhì)土滲流場與溫度場互相影響程度，結(jié)果表明，隨著滲透系數(shù)的增大，滲流場對溫度場的影響增大，而溫度場對滲流場的影響減弱，當(dāng)滲透系數(shù)小于10-9m/s時(shí)，滲流對于溫度的影響非常小，可以忽略，溫度是由熱傳導(dǎo)控制；當(dāng)滲透系數(shù)大于10-6m/s時(shí)，溫度基本上由滲流水控制，對流換熱遠(yuǎn)遠(yuǎn)超過了熱傳導(dǎo)，當(dāng)滲透系數(shù)大于10-5m/s時(shí)，溫度完全由滲流水控制，壩體的溫度等于水體的溫度。吳志偉和宋漢周[24]考慮常態(tài)水的黏度與密度的溫變效應(yīng)，建立了溫度場和變物性滲流場全耦合模型，在Femab軟件中進(jìn)行了求解，結(jié)果表明考慮水的溫變效應(yīng)后，溫度場普遍增大，滲流流速減小，且滲流場變動(dòng)幅度大于溫度場變動(dòng)幅度，體現(xiàn)了滲流場對溫度場的影響程度較反過來要大。

利用實(shí)測溫度資料、已知滲漏點(diǎn)位置等信息的土石壩滲流場反演分析研究價(jià)值較重大，目前國內(nèi)利用溫度場的滲流反演分析還比較少。針對各種類型的集中滲漏問題，董漢洲等[25]通過集中滲漏傳熱模型進(jìn)行反演分析，取得了不少研究成果；熊健[23]在大壩內(nèi)選取一些計(jì)算點(diǎn)，將反演轉(zhuǎn)化為最優(yōu)化求解問題，采用網(wǎng)絡(luò)搜索法對滲透系數(shù)進(jìn)行了反演，采用后驗(yàn)差法對反分析結(jié)果進(jìn)行統(tǒng)計(jì)檢驗(yàn)，得到了有效的反分析結(jié)果，證明其研究路線是可行的。

可以看出，溫度-滲流雙場定量分析模型經(jīng)歷了由簡單到完備的過程。雙場耦合模型考慮的耦合參數(shù)除了最初的對流傳熱項(xiàng)和溫度梯度流，還考慮了水的物理化學(xué)特性的變溫特性，以及滲透系數(shù)的溫變效應(yīng)，說明耦合模型更加符合天然滲流場與溫度場的耦合過程。同時(shí)，從數(shù)量上來看：①計(jì)算滲流場對溫度場的影響研究較多，由溫度場計(jì)算滲流場針對土石壩的研究較少，這是由于雙場互相影響程度不同，由溫度場分析滲流場屬于強(qiáng)非線性問題，其結(jié)果不唯一；② 針對土石壩的定量分析研究以“熱源法”為主，數(shù)值模擬較少，雖然對混凝土壩雙場耦合研究在基本原理上與土石壩大致相同，但土石壩的三相性特點(diǎn)與混凝土壩有較大區(qū)別，對滲流場溫度場均有較大影響，對土石壩的耦合數(shù)值模擬研究存在欠缺。

從滲流熱監(jiān)測技術(shù)的發(fā)展歷程可以看出，研究初期的重點(diǎn)是弄清雙場間作用機(jī)理，建立符合實(shí)際情況的耦合參數(shù)表達(dá)式，隨著對土石壩滲流場與溫度場互相影響規(guī)律的研究深入，研究重點(diǎn)轉(zhuǎn)移到對土石壩滲流場和溫度場基本理論、特有性質(zhì)的發(fā)掘，使模型趨于完善，能更準(zhǔn)確地描述土石壩滲流場與溫度場的復(fù)雜動(dòng)態(tài)平衡關(guān)系。

3 土石壩滲流計(jì)算與溫度計(jì)算研究現(xiàn)狀

3.1 土石壩滲流計(jì)算研究現(xiàn)狀

人們對土石壩滲流場計(jì)算方法的研究由來已久。1856年Darcy提出線性滲流率定理，奠定了滲流計(jì)算的基礎(chǔ)，1889年H.E茹可夫斯基首先推導(dǎo)了滲流場微分方程，許多研究者基于此模型建立了各種解析方法，但僅適用于少數(shù)的簡單邊界條件與均質(zhì)土壩滲流問題。隨著電子計(jì)算機(jī)技術(shù)飛速發(fā)展，各種數(shù)值模擬方法得到重視，1965年O.C.Zienkiewiz最早將有限單元法應(yīng)用到流體力學(xué)中，從此該方法在滲流計(jì)算中廣泛使用。1963年Jaswon首先提出了邊界元法的概念，并在隨后得到完善與發(fā)展。目前在滲流分析中使用的計(jì)算程序和各種計(jì)算方法，仍是由上述方法發(fā)展、完善而來的。我國對土石壩滲流數(shù)值模擬研究起步于20世紀(jì)70年代，取得了許多研究成果，并編制了一些比較成熟的土壩滲流與岸坡穩(wěn)定分析程序等[26]。可見，對滲流計(jì)算基本方法的研究已經(jīng)較為成熟，目前研究重點(diǎn)集中在一些滲流模擬的難點(diǎn)問題及其程序化上。文獻(xiàn)[27]對當(dāng)前巖土、水工專業(yè)滲流計(jì)算研究重點(diǎn)進(jìn)行了總結(jié)，其中包括“對復(fù)雜條件下滲流計(jì)算的研究，特別是非飽和滲流場計(jì)算”，和滲流熱監(jiān)測技術(shù)關(guān)系較為密切。

以往的土壩滲流計(jì)算較多是在飽和滲流假設(shè)下進(jìn)行的，即計(jì)算滲流場不考慮自由水面以上的土體，僅計(jì)算自由面以下飽和區(qū)的滲流，而實(shí)際上自由面以上的近乎飽和的毛管水帶以及更上的非飽和區(qū)，或稱為包氣帶中也存在滲流，由文獻(xiàn)[28]可知非飽和滲流的存在一定程度上加大了滲流量，在降雨、庫水位升降等易引發(fā)土石壩垮塌條件下，應(yīng)把飽和土與非飽和土考慮成一個(gè)整體來分析。建立土石壩飽和-非飽和模型進(jìn)行滲流場單場分析的研究較多，但基于飽和-非飽和滲流模型進(jìn)行雙場耦合數(shù)值模擬的研究較少，韋立德和楊春和[29]基于飽和-非飽和滲流模型，考慮了滲流場、溫度場對應(yīng)力場的單向影響，編制了三維有限元計(jì)算程序，其成果偏重應(yīng)力場研究，沒有考慮滲流場與溫度場間的耦合作用。吳志偉和宋漢周[30]對典型心墻土石壩建立了滲流場與溫度場耦合數(shù)學(xué)模型，考慮了黏度的變溫效應(yīng)，采用VG模型說明滲透系數(shù)與含水率的關(guān)系，能模擬復(fù)雜的非穩(wěn)態(tài)滲流問題及受滲流影響的熱流問題，可作為后續(xù)研究的參考，由于模型未考慮密度的溫變效應(yīng)與溫度勢流，其關(guān)于溫度場對滲流場的計(jì)算與真實(shí)情況有所差別。

3.2 土石壩溫度計(jì)算研究現(xiàn)狀

滲流區(qū)土石壩內(nèi)的溫度計(jì)算屬于含水多孔介質(zhì)傳熱問題，其中包含輻射、導(dǎo)熱、對流、擴(kuò)散、相變等多種熱學(xué)過程，實(shí)際土石壩中各點(diǎn)的熱流運(yùn)動(dòng)為上述各過程中的一種或多種耦合。多孔介質(zhì)傳熱問題是一個(gè)涉及領(lǐng)域非常多，研究歷史悠久的課題，文獻(xiàn)[31]認(rèn)為在探求遷移機(jī)理方面，該問題涉及的理論包括能量理論、液體擴(kuò)散理論、毛細(xì)流動(dòng)理論和蒸發(fā)凝結(jié)理論。由于大壩空間尺度較大，在進(jìn)行計(jì)算時(shí)須對這些理論進(jìn)行簡化。目前國內(nèi)的溫度場計(jì)算模型多由能量平衡方程建立，從土石壩中選取一個(gè)單元體，考慮上述熱學(xué)過程單位時(shí)間內(nèi)流入的熱量。根據(jù)能量守恒定律，流入熱量應(yīng)等于單元體單位時(shí)間的溫度變化，由此建立溫度場計(jì)算微分方程，求解溫度場分布?？梢姕囟葓瞿Ｐ托问饺Q于滲流場形式?，F(xiàn)有研究中絕大多數(shù)滲流場、溫度場耦合數(shù)學(xué)模型是在飽和滲流假定下建立的，因此在溫度場計(jì)算時(shí)默認(rèn)單元體的熱學(xué)性質(zhì)是常數(shù)，當(dāng)考慮土壤非飽和滲流問題時(shí)，其導(dǎo)熱系數(shù)受到含水率的影響。文獻(xiàn)[30]采用S.Lu等[32]的半理論模型對其進(jìn)行了描述。

目前溫度場也主要采用基于變分原理的有限元法計(jì)算，通用的大型有限元軟件如ANSYS，ABAQUS中均有專門的溫度模塊，在給定的邊界條件與定解條件下能方便地進(jìn)行計(jì)算。土石壩具有比較復(fù)雜的溫度邊界，通常難以獲得全面的實(shí)測資料，文獻(xiàn)[21]對各溫度邊界的計(jì)算方法有較詳細(xì)的說明。

4 研究展望

滲流熱監(jiān)測理論研究已經(jīng)取得不少成果，但為了提高滲流熱監(jiān)測技術(shù)理論的實(shí)用性與精度，有必要對滲流場與溫度場的耦合規(guī)律作進(jìn)一步探討。

(1) 鑒于多數(shù)耦合模型耦合參數(shù)不完備，考慮開展參數(shù)敏感度研究，分析模型對參數(shù)的影響程度，確定建立耦合參數(shù)完備的雙場全耦合模型的必要性，若敏感度較高，則應(yīng)在模型中考慮水的溫度特性、滲透系數(shù)的溫度變化等。

(2) 基于飽和-非飽和滲流建立的數(shù)值模型尚存在較多的問題，對降雨、庫水位驟變等需考慮非飽和滲流的情況缺乏考慮，應(yīng)借鑒非飽和土的研究成果，開展飽和-非飽和雙場耦合模型研究。

(3) 目前的雙場耦合模擬對某些問題的滯后效應(yīng)缺乏考慮，如滲流場與水位變化的滯后性、非飽和滲流中毛細(xì)水的滯后效應(yīng)等，這些問題需要進(jìn)一步研究。

(4) 由溫度場求解滲流場的研究比較欠缺，一方面應(yīng)借鑒其他領(lǐng)域中對非線性問題的解決方法，加強(qiáng)反演分析研究；另一方面，為驗(yàn)證數(shù)值模擬與真實(shí)值的差距，應(yīng)加強(qiáng)溫度觀測新技術(shù)(如分布式光纖)的應(yīng)用研究，以獲取更加實(shí)用的溫度時(shí)空分布資料。

5 結(jié) 語

土石壩的溫度場分布是某種形式熱源的作用結(jié)果，其變化是連續(xù)的，不同類型熱源的分布和作用方式在壩體中將引起不同的溫度場分布，因此，通過對土石壩溫度的監(jiān)測，開展集中滲漏或壩體均勻滲流而引起的溫度場溫度異常的研究分析，能夠建立溫度異常與土石壩滲透特性的對應(yīng)聯(lián)系?；谙嚓P(guān)研究，得到以下2條結(jié)論：

(1) 土石壩滲流場與溫度場的互相影響規(guī)律是滲流熱監(jiān)測技術(shù)的理論基礎(chǔ)，目前的雙場耦合模型滲流場對滲流對流換熱、熱傳效應(yīng)已有較充分的考慮，可以對穩(wěn)定/非穩(wěn)定滲流場對溫度場的影響進(jìn)行數(shù)值模擬。而溫度場對滲流場的影響如水的物化性質(zhì)的溫變性、滲透參數(shù)的溫變特性考慮不足，這會(huì)增加由溫度場求解滲流場的非線性強(qiáng)度，不利于反演收斂，降低反演精度。

(2) 在未來的研究中進(jìn)一步加強(qiáng)模型對邊界條件、介質(zhì)特性的仿真度，真實(shí)描述土石壩滲流場、溫度場互相影響的動(dòng)態(tài)平衡過程，須從土石壩多孔介質(zhì)流的機(jī)理出發(fā)，引入土壤中水分運(yùn)動(dòng)、熱量輸移的一些研究成果，提高模型的計(jì)算精度。

參考文獻(xiàn)：

[1] 陳 江,張少杰,閔少鑫,等. 熱脈沖法監(jiān)測集中滲漏的模型試驗(yàn)[J].巖土力學(xué)，2010,31(9):2729-2735. (CHEN Jiang，ZHANG Shao-jie，MIN Shao-xin,etal. Model Experiment for Concentrated Leakage Monitoring by Heat Pulse Method[J]. Rock and Soil Mechanics，2010,31(9):2729-2735. (in Chinese))

[2] 王新建,許寶田,陳建生. 溫度場探測堤壩集中滲漏綜述[J].水利水電科學(xué)進(jìn)展，2008,28(3):89-94. (WANG Xin-jian, XU Bao-tian，CHEN Jian-sheng. Status and Prospects of Research on Temperature Field Tracer Method Used in Leakage Detection[J]. Advances in Science and Technology of Water Resources, 2008,28(3):89-94. (in Chinese))

[3] MYXER H A,ДИНОВ. 將溫度作為土石壩性態(tài)的指示因子[J].蘇伯林，譯. 大壩觀測技術(shù),1994, (2):22-24. (MYXER H A, ДИНОВ. Temperature as Indicators of Earth Rock Dam Behavior[J]. Translated by SU Bo-lin. Dam Observation Technology, 1994, (2):22-24. (in Chinese))

[4] ARMBRUSTER H, TROGER .溫度量測法探測土中的滲流[J]. 林叔忠，譯. 廣東水利水電,1994, (3):15-17. (ARMBRUSTER H, TROGER I. Temperature Measurement in Detecting Seepage in the Soil[J]. Translated by LIN Shu-zhong. Guangdong Hydropower, 1994, (3):15-17. (in Chinese))

[5] 王志遠(yuǎn),王占銳. 一項(xiàng)滲流監(jiān)測新技術(shù)——排水孔測溫法[J].大壩觀測與土工測試,1997,21(4):5-7. (WANG Zhi-yuan，WANG Zhan-rui. A New Technology of Leakage Monitoring: Measuring Temperature by Drain Hole[J]. Dam Observation and Geotechnical Tests, 1997,21(4):5-7. (in Chinese))

[6] 肖才忠,潘文昌. 由溫度場研究壩基滲流初探[J].人民長江,1999,30(5):21-23. (XIAO Cai-zhong, PAN Wen-chang. Preliminary Analysis on Dam Foundation Seepage through Temperature Field[J]. Yangtze River, 1999,30(5):21-23. (in Chinese))

[7] 吳 勇. 利用溫度場探測土石壩滲漏問題的正反演研究——以南方某水庫為例[D].南京：河海大學(xué),2008. (WU Yong. Forward Analysis and Inverse Analysis on the Seepage Problem of Earth-rock Dam in Consideration of Temperature Field: Case Study on a Reservoir in South China[D]. Nanjing: Hohai University, 2008. (in Chinese))

[8] 陳建生,董海洲,吳慶林,等. 虛擬熱源法研究壩基裂隙巖體滲漏通道[J].巖石力學(xué)與工程學(xué)報(bào),2005,24(22):4019-4024. (CHEN Jian-sheng, DONG Hai-zhou, WU Qing-lin，etal. Detection of Leakage Passage in Fissure Rock with Assumptive Heat Source Method[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(22): 4019-4024.(in Chinese))

[9] 陳建生,董海洲,余 波,等. 利用線熱源法研究堤防集中滲漏通道[J].地球物理學(xué)進(jìn)展，2003,18(3):401-403. (CHEN Jian-sheng, DONG Hai-zhou, YU Bo,etal. Using Line-heat Resource Method to Study Concentration Leakage Passage of Dike[J]. Progress in Geophysics, 2003,18(3):401-403. (in Chinese))

[10] 王新建,曾長女,陳建生. 利用穩(wěn)定溫度場探測堤壩面狀集中滲漏[J].人民黃河,2008,30(9):1800-1805. (WANG Xin-jian, ZENG Zhang-nv, CHEN Jian-sheng. Detecting Concentrated Seepage of Dam Surface by Using Stable Temperature Field[J]. Yellow River, 2008, 30(9): 1800 -1805. (in Chinese))

[11] 王新建,陳建生. 溫度探測土壩圓柱狀集中滲漏模型研究[J].水文地質(zhì)與工程地質(zhì)，2006, (4):31-36. (WANG Xin-jian, CHEN Jian-sheng. A Study of Temperature Field Tracer Method Used to Detect Tubular Leakage in Dam[J]. Hydrogeology & Engineering Geology, 2006, (4):31-36. (in Chinese))

[12] 王新建,潘紀(jì)順. 堤壩多集中滲漏通道位置溫度探測研究[J].巖土工程學(xué)報(bào),2010,32(11):1800-1805. (WANG Xin-jian，PAN Ji-shun. Location Detection of Concentrated-leakage Passages in Dam by Groundwater Temperature[J]. Chinese Journal of Geotechnical Engineering, 2010,32(11):1800-1805. (in Chinese))

[13] 朱伯芳. 大體積混凝土溫度應(yīng)力與溫度控制[M].北京:中國電力出版社,1999. (ZHU Bo-fang. Temperature Stress and Temperature Control of Mass Concrete [M]. Beijing: China Electric Power Press, 1999. (in Chinese))

[14] 賴遠(yuǎn)明,吳紫江,朱元林,等. 寒區(qū)隧洞溫度場、滲流場和應(yīng)力場耦合問題的非線性分析[J].巖土工程學(xué)報(bào),1999,21(5): 529-533. (LAI Yuan-ming, WU Zi-jiang, ZHU Yuan-lin,etal. Nonlinear Analyses for the Coupling Problem of Temperature, Seepage and Stress Fields in Cold Region Tunnels[J]. Chinese Journal of Geotechnical Engineering, 1999,21(5): 529-533. (in Chinese))

[15] 柴軍瑞. 混凝土壩滲流場與穩(wěn)定溫度場耦合分析的數(shù)學(xué)模型[J].水力發(fā)電學(xué)報(bào),2000, (1):27-35. (CHAI Jun-rui. On Mathematical Model for Coupled Seepage and Temperature Field in Concrete Dam[J]. Journal of hydroelectric engineering, 2000, (1):27-35. (in Chinese))

[16] 柴軍瑞,韓 群,仵彥卿. 巖體一維滲流場與溫度場耦合模型的解析演算[J].地下水,1999,21(4):180-181. (CHAI Jun-rui, HAN Qun, WU Yan-qin. Analytical Calculation of One-dimensional Coupled Model of Seepage Field and Temperature Field in Rock Mass[J]. Underground Water, 1999,21(4):180-181. (in Chinese))

[17] 陳建余,朱岳明,張建斌. 考慮滲流場影響的混凝土壩溫度場分析[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2003,31(2):119-123. (CHEN Jian-yu, ZHU Yue-ming, ZHANG Jian-bin. Temperature Fields in Concrete Dams in Consideration of Seepage Fields[J]. Journal of Hohai University (Natural Sciences), 2003,31(2):119-123. (in Chinese))

[18] 許增光,李康宏,柴軍瑞,等. 考慮滲流熱學(xué)效應(yīng)的大壩穩(wěn)定溫度場有限元分析[J].紅水河,2006,25(2):112-115. (XU Zeng-guang, LI Kang-hong, CHAI Jun-rui,etal. Finite Element Numerical Analysis of Steady Temperature Field of Dam Considering Thermal Effect of Seepage[J]. Hongshui River,2006,25(2): 112-115. (in Chinese))

[19] 陳向軍,李志敏. 考慮滲流場影響的大壩非穩(wěn)定溫度場[J].長沙理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2008,5(1):44-47. (CHEN Xiang-jun, LI Zhi-min. Unsteady Temperature Field in Dam in Consideration of Seepage Field[J]. Journal of Changsha University of Science and Technology(Natural Science), 2008,5(1):44-47. (in Chinese))

[20] 崔浩東,朱岳明. 蓄水初期的壩體非穩(wěn)定滲流場與溫度場耦合的理論模型及數(shù)值模擬[J].水利學(xué)報(bào),2009, 40(2):238- 242. (CUI Hao-dong, ZHU Yue-ming. Theoretical Model and Numerical Simulation by Coupling Unsteady Seepage Field with Temperature Field During Initial Impounding Period of Concrete Dams[J]. Journal of Hydraulic Engineering, 2009, 40(2):238-242. (in Chinese))

[21] 李端有,熊 健,於三大,等. 土石壩滲流熱監(jiān)測技術(shù)研究[J].長江科學(xué)院院報(bào),2005,22(6):29-33. (LI Duan-you, XIONG Jian, YU San-da,etal. Study on Seepage Monitoring Technology with Temperature in Embankment Dam[J]. Journal of Yangtze River Scientific Research Institute, 2005, 22(6): 29-33. (in Chinese))

[22] 孫 曼. 光纖Bragg光柵傳感技術(shù)用于工程結(jié)構(gòu)安全監(jiān)測的研究[D].成都：四川大學(xué),2005. (SUN Man. Safety Monitoring of Engineering Structures Using Fiber Bragg Grating Sensing Technology[D]. Chengdu： Sichuan University,2005. (in Chinese))

[23] 熊 健. 土石壩滲流熱監(jiān)測理論研究[D].武漢：長江科學(xué)院,2006. (XIONG Jian. Study on the Seepage Monitoring Theory with Temperature in Embankment Dam[D]. Wuhan: Yangtze River Scientific Research Institute,2006. (in Chinese))

[24] 吳志偉,宋漢周. 壩址溫度場與變物性滲流場全耦合分析[J].水利學(xué)報(bào),2010,41(6):703-710. (WU Zhi-wei, SONG Han-zhou. Fully-coupled Analysis of Temperature Field and Variable Properties Seepage Field around Dam-site[J]. Journal of Hydraulic Engineering, 2010,41(6):703-710. (in Chinese))

[25] 董漢洲,羅日洪,張 令,等. 壩基雙集中滲流通道傳熱模型及流速反演研究[J].四川大學(xué)學(xué)報(bào)(工程科學(xué)版),2012,44(3): 36-41. (DONG Han-zhou, LUO Ri-hong, ZHANG Ling,etal. Study on Double-concentrated-leakage Passage Heat Conduction Model and Retrieval of Velocity in Dam Basement[J]. Journal of Sichuan University (Engineering Science), 2012,44(3): 36-41. (in Chinese))

[26] 姜 帆,苾永寧,張 茹. 土石壩滲流研究發(fā)展綜述[J].水利與建筑工程學(xué)報(bào),2006,4(4):94-97. (JIANG Fan, MI Yong-ning, ZHANG Ru. Development of Research on Seepage of Embankment Dam[J]. Journal of Water Resources and Architectural Engineering, 2006, 4(4):94-97. (in Chinese))

[27] 黃站峰. 堤壩飽和-非飽和滲流研究及其程序化[D].武漢：武漢大學(xué),2005. (HUANG Zhan-feng. Study and Programming of Saturate & Saturate-Unsaturated Seepage of Embankment[D]. Wuhan: Wuhan University,2005. (in Chinese))

[28] 王增欣. 土石壩飽和-非飽和滲流的有限元分析及應(yīng)用[D].大連：大連理工大學(xué),2010. (WANG Zeng-xin. Finite Element Analysis of Saturated-Unsaturated Seepage for Embankment Dams and Application[D]. Dalian: Dalian University of Technology, 2010. (in Chinese))

[29] 韋立德,楊春和. 考慮飽和-非飽和滲流、溫度和應(yīng)力耦合的三維有限元程序研制[J].巖土力學(xué),2005,26(6):1000-1004. (WEI Li-de, YANG Chun-he. Program Development for 3D FEM of Coupling Saturated-Unsaturated Seepage, Temperature and Stress[J]. Rock and Soil Mechanics, 2005,26(6):1000-1004. (in Chinese))

[30] 吳志偉,宋漢周. 堤壩飽和-非飽和滲流場與溫度場耦合數(shù)值模擬[C]∥2013年中國水力發(fā)電工程學(xué)會(huì)大壩安全監(jiān)測專委會(huì)年會(huì)暨大壩安全評價(jià)技術(shù)交流會(huì)論文集.丹東,2013：273-279. (WU Zhi-wei, SONG Han-zhou. Numerical Simulation of Coupling Saturated-Unsaturated Seepage Field and Temperature Field[C]∥Proceedings of the 2013 Annual Meeting of China Hydropower Engineering Society. Dandong, 2013: 273-279. (in Chinese))

[31] 施明恒,虞維平,王補(bǔ)宣. 多孔介質(zhì)傳熱傳質(zhì)研究的現(xiàn)狀與展望[J].東南大學(xué)學(xué)報(bào),1994,24(增刊):1-5. (SHI Ming-heng, YU Wei-ping, WANG Bu-xuan. Research Status and Prospect of Heat and Mass Transfer in Porous Media[J]. Journal of Southeast University，1994,24(Sup.):1-5. (in Chinese))

[32] LU S, REN T, GONG Y S. An Improved Model for Predicting the Hydraulic Conductivity of Unsaturated Soils[J].Soil Science Society of America Journal, 1980, (44):892-898.