許孝臣，許小杰，占　川，李洪強(qiáng)，湯明禮

（1.浙江廣川工程咨詢有限公司，浙江　杭州　310020；2.浙江省水利河口研究院，浙江　杭州　310020）

?

土石壩滲流監(jiān)測資料分析方法研究

許孝臣1，2，許小杰1，2，占川1，2，李洪強(qiáng)1，2，湯明禮1，2

（1.浙江廣川工程咨詢有限公司，浙江杭州310020；2.浙江省水利河口研究院，浙江杭州310020）

為更好地發(fā)揮土石壩滲流監(jiān)測資料的作用，以浙江省現(xiàn)狀條件下的監(jiān)測狀況、已研究的滲流分析方法為基礎(chǔ)，結(jié)合浙江省水利工程標(biāo)準(zhǔn)化建設(shè)的要求，研究了滲流監(jiān)測資料的分析方法，通過分析滲流在時間、空間上的分布規(guī)律及其影響因素，對滲流穩(wěn)定和大壩安全做出判斷，有效發(fā)揮滲流監(jiān)測資料的作用。

滲流監(jiān)測；資料分析；土石壩；分析模型

1　問題的提出

土石壩是當(dāng)?shù)夭牧蠅?，有諸多優(yōu)點，如造價低、施工技術(shù)成熟、便于管理等，在國內(nèi)外的應(yīng)用較廣。據(jù)統(tǒng)計，土石壩在國外占所有大壩的比例超過62.0%［1］，國內(nèi)土石壩所占的比例更高，尤其是在小型水庫中，表現(xiàn)更為明顯，如浙江省擁有3 885座小型水庫，其中土石壩所占的比例高達(dá)90.5%［2］。

土石壩所面臨的問題主要有2個，即結(jié)構(gòu)問題和滲流問題。土顆粒在高水頭滲流力的作用下容易被帶走，進(jìn)而發(fā)展為滲透破壞，表現(xiàn)形式有管涌、流土等。據(jù)統(tǒng)計，在土石壩所有的安全問題中，由滲流引起的占37.1%［3］。因此，土石壩的滲流穩(wěn)定是影響大壩安全的重要因素，滲流分析是大壩安全分析的重要內(nèi)容。

目前大量的土石壩中埋設(shè)了滲壓計，可以實現(xiàn)對壩體內(nèi)滲流場的實時監(jiān)控，為分析和研究壩體內(nèi)滲流場的變化提供了寶貴資料，但目前對土石壩監(jiān)測資料分析的深度和廣度有待提高。

從2015年年底開始，為規(guī)范水利工程管理，浙江省迅速地推廣水利工程標(biāo)準(zhǔn)化建設(shè)，把大壩監(jiān)測工作作為水庫管理的重要內(nèi)容，對大壩監(jiān)測在數(shù)據(jù)采集、資料整編、數(shù)據(jù)分析方面都作出了明確規(guī)定。數(shù)據(jù)采集、資料整編可以通過自動化或軟件實現(xiàn)，但資料分析相對困難。因此，本文結(jié)合浙江省土石壩的滲流監(jiān)測現(xiàn)狀、以往滲流資料的分析方法研究，研究了基于土石壩滲流監(jiān)測資料的滲流分析方法，以期為浙江省水利工程標(biāo)準(zhǔn)化建設(shè)中的滲流分析提供參考。

2　滲流壓力分析

2.1過程線分析

過程線分析是規(guī)范中推薦的方法，以時間為橫軸，以各觀測項目的工程物理量為縱軸，繪制物理量隨時間變化的過程線。為分析各影響因子之間的關(guān)系，常把庫水位、降雨量、溫度等影響因子繪制在同一張圖上，以便分析不同物理量在時間和空間上的變化過程及變化趨勢。著重觀察不同變量之間的關(guān)系，分析測值和因變量之間的因果關(guān)系，查看有無持續(xù)增加或減小的趨勢，有無突變，結(jié)合土石壩滲流分布規(guī)律，判斷是否存在異常。

2.2浸潤線分析

在土石壩斷面圖上繪制壩體浸潤線，可直觀地看出斷面上滲流壓力的分布規(guī)律。將實測壩體浸潤線與設(shè)計壩體浸潤線進(jìn)行比較，當(dāng)實測壩體浸潤線高于設(shè)計壩體浸潤線，或出現(xiàn)其它異常時，應(yīng)進(jìn)行詳細(xì)分析。

2.3相關(guān)性分析

通過計算不同測點或不同時段內(nèi)滲流壓力與水位的相關(guān)系數(shù)，可分析判斷不同位置測點與庫水位的相關(guān)性大小，根據(jù)滲流壓力分布規(guī)律判斷其是否合理。計算不同時段的相關(guān)系數(shù)，可以分析判斷滲流發(fā)展趨勢，向惡化或者改善方向發(fā)展。

2.4滯后性分析

滲流壓力的滯后性有2種，一種是測量裝置的滯后性，一種是庫水位變化導(dǎo)致測點發(fā)生變化的滯后性，分析中應(yīng)加以區(qū)別。滲流測量裝置一般有測壓管或埋入式滲壓計，測壓管測量滲流壓力是通過水滲入和流出測壓管內(nèi)導(dǎo)致管內(nèi)水位變化，根據(jù)水位上升或下降，從而反映滲流壓力的變化。測壓管補(bǔ)充或排出水量需要的時間，導(dǎo)致測壓管出現(xiàn)滯后性，這種效應(yīng)一般在滲透性較小的土中表現(xiàn)得特別明顯，因此《土石壩安全監(jiān)測技術(shù)規(guī)范》規(guī)定，在土的滲透系數(shù)小于1×10-5cm/s的情況下，不適合埋設(shè)測壓管進(jìn)行觀測。埋入式滲壓計不存在測壓管的滯后性問題，在有壓飽和土體內(nèi)，滲流壓力的傳遞十分迅速，時間可以忽略不計，但土石壩為無壓滲流，常表現(xiàn)出一定的滯后性，分析中應(yīng)明確滯后性產(chǎn)生的原因。

采用單因子回歸模型可以計算滲流壓力與水位的滯后時間。通過分析滯后性的原因和時間長短，可以判定測點與庫水位的連通性。

2.5滲流分布等值線分析

通過繪制滲流分布等值線圖，可以觀察壩體不同部位的滲流分布情況，分析滲流分布規(guī)律，解釋滲流規(guī)律成因。繪制可以借助現(xiàn)有軟件，例如采用SUFER或其它程序進(jìn)行繪制。

2.6壩基抗?jié)B分析

計算各斷面相鄰測點的滲透坡降，與允許滲透坡降進(jìn)行比較，判斷滲透穩(wěn)定性。

2.7位勢分析

拉普拉斯方程可以用來描述穩(wěn)定的大壩滲流場。對于建在強(qiáng)透水地基上、采用鋪蓋進(jìn)行防滲的壩基滲流，可以簡化為有壓滲流。根據(jù)拉普拉斯方程，勢函數(shù)在滲透介質(zhì)固定的情況下，是空間位置的函數(shù)，不隨時間和庫水位的變化而變化。對于滲流測點，空間位置確定，作為無量綱量的勢函數(shù)也不隨時間變化，為常數(shù)，一旦發(fā)生變化，說明滲透介質(zhì)發(fā)生改變。不能簡化為有壓滲流的情況下，可以分析在同一庫水位情況下的位勢，也就是庫水位一定，滲透介質(zhì)一定，位勢僅是空間位置的函數(shù)，在一定庫水位情況下，位勢可視為常數(shù)，若發(fā)生變化，說明滲透介質(zhì)發(fā)生變化，可據(jù)此來分析滲流監(jiān)測資料。

勢函數(shù)可用位勢表示。位勢和測滲壓水位的關(guān)系見式（1）：

式中：φ為位勢，無量綱量；H為滲透壓力水位，m；H1為上游水位，m；H2為下游水位，m。

通過對過程線的觀測與分析，判斷滲流安全，對分析判斷準(zhǔn)則進(jìn)行歸納（見表1）［4］。

表1　位勢分析準(zhǔn)則表

2.8回歸模型分析

采用多元回歸方法建立滲流壓力統(tǒng)計模型［5］，統(tǒng)計模型的因子為上游水深、降水量、時間等，根據(jù)相關(guān)系數(shù)，就回歸效果給出定性結(jié)論。回歸模型見式（2）～（4）。

根據(jù)統(tǒng)計模型可以分析并預(yù)測滲流發(fā)展趨勢，同時對水位、雨量、時效線性項和時效對數(shù)項進(jìn)行分離變量分析。確定各個因素的影響程度，找出主要影響因素。當(dāng)主要影響因素為時效線性項或水位時，應(yīng)分析主要影響因素的逐年變化規(guī)律，確定大壩是否存在安全隱患。

3　滲流量分析

滲流量過程線分析一般以滲流量為縱坐標(biāo)，以時間為橫坐標(biāo)，作滲流量過程線圖，同時繪出上下游水位過程線和降雨量分布圖，這些是常用的資料整理方法。通過對比滲流量和庫水位、降雨量的變化過程和相關(guān)程度，判斷滲流量的影響因素及變化規(guī)律。對滲流量同樣可以進(jìn)行相關(guān)性分析和回歸模型分析，方法同滲流壓力分析。

3.1根據(jù)測壓管水位估算壩基滲流量

通常情況下，浙江省的土石壩下游無水，沒有條件安裝滲流量觀測設(shè)施時，無法直接觀測滲流量，可以將大壩壩基等效為有限均勻地基，通過壩下游的測壓管水位換算壩基滲流量。均勻有限地基上的無壓滲流量計算示意見圖1。

圖1　壩基滲流量計算示意圖

滲流速度可采用達(dá)西定律，方程如下：

單寬流量方程如下：

分離變量積分可以得到如下方程：

帶入邊界條件（x = 0時，y = h1；x = l時，y = h2），可以得到單寬流量如下計算方程：

3.2化引流量分析

化引流量可用如下方程進(jìn)行計算：

式中：Q為總滲流量，m3/s；B為過水寬度，m。

對具體工程，k和B是常數(shù)，因此，式（9）可變?yōu)椋?/p>

化引流量和位勢分別相當(dāng)于流函數(shù)和勢函數(shù)，根據(jù)式（10），化引流量可作為無量綱量，其分析方法與位勢基本相同。對于有壓滲透，化引流量不隨時間變化；對于無壓滲流，在同一庫水位下，化引流量不隨時間變化，若發(fā)生變化，則滲透介質(zhì)發(fā)生變化［6］。將化引流量與周圍滲流壓力結(jié)合進(jìn)行分析，具體分析準(zhǔn)則見表2［4］。

表2　化引流量分析準(zhǔn)則表

4　結(jié)　語

針對浙江省現(xiàn)狀條件下土石壩的監(jiān)測項目及監(jiān)測情況，在前人所研究的滲流分析方法的基礎(chǔ)上，研究了滲流監(jiān)測資料的分析方法，通過對滲流壓力和滲流量的過程線、分布圖、相關(guān)性、滯后性、位勢、化引流量及滲流模型進(jìn)行分析，結(jié)合滲流在時間、空間上的分布規(guī)律及影響因素的分析，可以清楚地了解滲流的分布規(guī)律、分布成因及發(fā)展趨勢，對滲流穩(wěn)定和大壩安全做出判斷，有效發(fā)揮滲流監(jiān)測資料的作用，為大壩安全運行提供技術(shù)支撐，以期為水利工程標(biāo)準(zhǔn)化建設(shè)的大壩監(jiān)測資料分析提供參考。

［1］孫少華，謝紅強(qiáng)，左林勇.基于有限元法的土石壩滲流特性及防滲效果研究［J］.礦業(yè)研究與開發(fā)，2009，29（3）：27 - 29.

［2］王良，崔英華.浙江省小型水庫土石壩病險分析及對策［J］.浙江水利科技，2009（1）：26 - 28.

［3］丁樹云，蔡正銀.土石壩滲流研究綜述［J］.人民長江，2008，39（2）：33 - 36.

［4］龐毅，王保澤，楊萬志，等.土石壩安全監(jiān)測分析評價技術(shù)與工程應(yīng)用［M］.北京：中國水利水電出版社，2011.

［5］吳中如.水工建筑物安全監(jiān)控理論及其應(yīng)用［M］.北京：高等教育出版社，2003.

［6］王韋.土壩滲流觀測資料的分析及應(yīng)用［J］.水利水運科學(xué)研究，1980（2）：10 - 19.

（責(zé)任編輯黃超）

Study on the Analysis Method of Embankment Dam Seepage Monitoring Data

XU Xiao - chen1，2，XU Xiao - jie1，2，ZHAN Chuan1，2，LI Hong - qiang1，2，TANG Ming - li1，2
（1.Zhejiang Guangchuan Engineering Consulting Co.，Ltd.，Hangzhou 310020，Zhejiang，China；2.Zhejiang Institute of Hydraulics & Estuary，Hangzhou 310020，Zhejiang，China）

To better play the role of embankment dam seepage monitoring data，on the basis of monitoring situation under present conditions in Zhejiang Province and previous seepage analysis methods，this paper studied on the analysis method of seepage monitoring data，combining with the requirements of water conservancy standardization construction in Zhejiang Province.By analyzing the distribution rule and its controlling factors of the seepage flow in time and space，it made judgments on seepage stability and dam safety，playing the role of seepage monitoring data effectively.

seepage monitoring；data analysis；embankment dam；analysis model

TV223.4

A

1008 - 701X（2016）05 - 0021 - 03

10.13641/j.cnki.33 - 1162/tv.2016.05.009

2016-07-22

浙江省科技計劃項目（2014F50016）。

許孝臣（1982 - ），男，工程師，碩士，主要從事大壩安全監(jiān)測、巖土滲流與控制工作。