周鵬 王霜 王愷

摘 要：利用有限元分析軟件模擬雙層堤基、三層堤基以及多層堤基內(nèi)部的水力梯度分布，將堤基內(nèi)部各區(qū)域的水力梯度與相應(yīng)的允許水力梯度進(jìn)行比較，最終確定堤基管涌發(fā)生的臨界上游水頭。三層堤基和多層堤基中，堤基內(nèi)部存在細(xì)砂夾層，隨著細(xì)砂夾層埋深的增加，管涌發(fā)生的臨界上游水頭逐漸增大，通道深度不斷增加，管涌破壞發(fā)生的標(biāo)志是細(xì)砂夾層被沖破，發(fā)生豎向流土破壞;當(dāng)細(xì)砂夾層埋深達(dá)到一定深度時(shí)，細(xì)砂上部砂礫石層的滲水量增加，增加到足以將砂礫石表面的顆粒帶出，管涌破壞過(guò)程與雙層堤基類(lèi)似，通道主要形成在黏土下部的砂礫石表面，臨界水頭有所減小。將該數(shù)值模擬結(jié)果與前人管涌試驗(yàn)結(jié)果進(jìn)行比較，兩者吻合度較好。

關(guān)鍵詞：堤基結(jié)構(gòu);管涌;數(shù)值模擬;臨界水頭

Abstract：A finite element software was used to simulate the distribution of hydraulic gradient in two-stratum， three-stratum and multi-stratum foundations. The hydraulic gradient in different areas was compared to the critical hydraulic gradient of soil erosion in the corresponding area and then the critical hydraulic head of piping in the dike foundation was determined. In three-stratum and multi-stratum foundations， the critical hydraulic head gradually increased with the burial depth of the fine sand layer and the depth of the piping passage was gradually increased. Piping happened when the sand layer broke through with the vertical soil flow. When the burial depth of the fine sand layer came to a certain value， the piping passage formed on the surface of the sandy gravel layer below the clay layer， which was similar with that in two-stratum foundation. The critical hydraulic head decreased. Comparisons between the simulation results and the experimental value were conducted and the results showed that the goodness of fit was high.

Key words： structure of dike foundation; piping; numerical modeling; critical hydraulic head

歷史上的河流沖刷、淤積和改道，以及沉積條件的千變?nèi)f化，導(dǎo)致不同時(shí)期、不同河段所形成的堤基結(jié)構(gòu)有很大差異。堤基大多由黏土、粉質(zhì)黏土、壤土、淤泥質(zhì)土、砂壤土、粉細(xì)砂和砂礫石等組成，有多種組合類(lèi)型，大體上可以概括為單層堤基、雙層堤基、三層堤基和多層堤基4種類(lèi)型[1-3]。對(duì)于不同結(jié)構(gòu)的堤基，滲流場(chǎng)也是不同的，應(yīng)該對(duì)具體條件下的滲流場(chǎng)進(jìn)行分析。隨著計(jì)算機(jī)技術(shù)的迅速發(fā)展，數(shù)值模擬已經(jīng)漸漸成為管涌研究的一個(gè)重要手段和方法。殷建華[4]在國(guó)內(nèi)首次采用有限元方法來(lái)研究管涌問(wèn)題，計(jì)算了管涌發(fā)生區(qū)域的長(zhǎng)度和滲透系數(shù)對(duì)堤身滲流的影響。張家發(fā)等[5]針對(duì)長(zhǎng)江中下游堤防工程現(xiàn)狀概化出典型條件，采用有限元數(shù)值模擬方法對(duì)堤身堤基滲流場(chǎng)進(jìn)行了模擬計(jì)算，分析了作用水頭、堤身滲透性、弱透水覆蓋層厚度、強(qiáng)透水層厚度、外灘寬度和河泓切割程度等因素對(duì)滲流規(guī)律的影響。周曉杰等[6-7]將計(jì)算區(qū)域分為管流域與滲流域兩部分，對(duì)兩域分別采用管流理論與達(dá)西定律求解，通過(guò)對(duì)兩域共同邊界進(jìn)行耦合計(jì)算，從而獲得整個(gè)流場(chǎng)的信息。

本文針對(duì)上述不同結(jié)構(gòu)堤基，采用有限元數(shù)值模擬方法對(duì)不同作用水頭下堤基滲流場(chǎng)進(jìn)行分析，通過(guò)比較堤基內(nèi)部各區(qū)域的水力梯度與對(duì)應(yīng)的允許水力梯度的大小，判斷堤基滲透穩(wěn)定狀態(tài)，最終確定管涌發(fā)生時(shí)的臨界上游水頭。

1 堤基內(nèi)各區(qū)域臨界水力梯度計(jì)算

通過(guò)有限元數(shù)值模擬得到不同上游水頭下堤基內(nèi)部的水力梯度分布情況，將堤基內(nèi)部各區(qū)域的水力梯度與允許水力梯度進(jìn)行比較，從而確定在該上游水頭下堤基內(nèi)部是否有顆粒侵蝕。若無(wú)侵蝕，繼續(xù)逐級(jí)增大上游水頭，最終得到堤基內(nèi)部土體發(fā)生侵蝕時(shí)的上游水頭，也就是管涌發(fā)生的臨界上游水頭，同時(shí)可以得到堤基內(nèi)部顆粒侵蝕的位置與區(qū)域。砂礫石層以及細(xì)砂夾層的顆粒級(jí)配曲線見(jiàn)圖1、物理參數(shù)見(jiàn)表1，上覆黏土的滲透系數(shù)為5×10-6 m/s。

多層堤基中，存在滲透系數(shù)較小的細(xì)砂夾層，將堤基的砂礫石層分成夾層上部與夾層下部?jī)刹糠?。堤基發(fā)生滲透破壞存在兩種形式：夾層上部，在黏土層與砂礫石層界面發(fā)生侵蝕破壞;夾層下部，細(xì)砂夾層被沖破，發(fā)生豎向流土破壞，然后侵蝕逐漸向上游發(fā)展。不同堤基結(jié)構(gòu)發(fā)生滲透破壞的形式不同，主要取決于哪種滲透破壞形式的臨界水力條件先得到滿足，因此首先需要對(duì)這兩種滲透破壞的臨界條件進(jìn)行確定。

對(duì)于黏土層與砂礫石層界面上發(fā)生的水平方向侵蝕破壞，管涌口處的顆粒需首先滿足豎向運(yùn)動(dòng)的條件，顆粒沿豎向發(fā)生侵蝕，在管涌口處形成空洞，才能為水平方向的侵蝕提供足夠的空間，侵蝕才能沿水平方向逐漸向上游發(fā)展。根據(jù)康德拉且夫抗?jié)B比降公式[8]計(jì)算允許水力梯度Jcr：

對(duì)于細(xì)砂夾層被沖破而發(fā)生的豎向流土破壞，以管涌口下部單位底面積的夾砂層為研究對(duì)象，所受力主要有自身的浮重G′=γ′s2h2（1-n2），上部砂礫石的壓力F=γ′s1h1（1-n1），滲透力f=γwJh2，根據(jù)靜力平衡條件可得：

將堤基各層土體的物理參數(shù)代入式（5）可得，在不同的細(xì)砂層埋深下，細(xì)砂層及上部土體被沖破時(shí)細(xì)砂層底部所需的豎向臨界水力梯度（見(jiàn)表2）。

本文采用逐級(jí)提升水頭的方式模擬管涌集中滲漏通道的發(fā)生、發(fā)展。當(dāng)?shù)袒苡靠谏暗[石內(nèi)的水力梯度小于0.36或者細(xì)砂夾層底部水力梯度小于表2中所對(duì)應(yīng)的值時(shí)，兩種滲透破壞均不發(fā)生，顆粒不發(fā)生流失，限于文章篇幅，未在文中分析。

2 不同堤基結(jié)構(gòu)管涌發(fā)生的臨界上游水頭

2.1 雙層堤基管涌發(fā)生的臨界上游水頭

首先對(duì)典型的雙層堤基上覆黏土層出現(xiàn)管涌口后堤基內(nèi)部滲流場(chǎng)進(jìn)行模擬。模型以向堤外方向?yàn)閄軸正方向、向上為Y軸正方向建立坐標(biāo)系，堤基向兩邊無(wú)限延伸。模擬的是半無(wú)限邊界，上游側(cè)水頭的邊界條件作用于模型左側(cè)，下游側(cè)水頭的邊界條件作用于模型右側(cè)的黏土層表面。不考慮上部堤防滲流對(duì)堤基內(nèi)部滲流場(chǎng)產(chǎn)生的影響，假定堤防為不透水體，同時(shí)假定距堤腳26 m處為黏土層薄弱處，薄弱區(qū)域直徑為4 m，即在模型中將X軸坐標(biāo)為84～88 m的黏土層定為管涌口，滲透系數(shù)設(shè)定為0.5 m/s。

圖2為管涌口形成后、上游水頭H=10 m時(shí)堤基內(nèi)部的水頭等值線圖，可以看出，管涌口形成后，壓力水頭由砂礫石層承擔(dān)，在X軸向上，隨著X值的增大，壓力水頭呈逐漸下降趨勢(shì)，在管涌口附近壓力水頭發(fā)生突變，即管涌口附近的水力梯度最大。

雙層堤基的滲透破壞主要發(fā)生在黏土層與砂礫石層界面上，故土體的允許水力梯度是由式（1）所確定的0.36。將模擬的堤基內(nèi)部各區(qū)域水力梯度分布情況與該值進(jìn)行比較，可知當(dāng)H<19 m時(shí)不會(huì)出現(xiàn)顆粒流失，所以初始水頭定為19 m。圖3為該水頭下管涌口處Y方向水力梯度分布，管涌侵蝕的區(qū)域根據(jù)水力梯度分布圖和允許水力條件共同確定，將水力梯度分布圖上水力梯度超過(guò)顆粒允許水力梯度的區(qū)域，認(rèn)定為顆粒發(fā)生流失的區(qū)域。由圖3可以看出，顆粒發(fā)生流失的區(qū)域?yàn)閄=84.0～84.5 m。因此，在雙層堤基中，發(fā)生管涌的臨界上游水頭為19 m。

2.2 三層堤基管涌發(fā)生的臨界上游水頭

三層堤基中，細(xì)砂層在砂礫石層表面，也就是細(xì)砂層在砂礫石層中的埋深為0 m，當(dāng)H<4 m時(shí)不會(huì)出現(xiàn)顆粒流失，初始水頭定為4 m。該種堤基結(jié)構(gòu)發(fā)生的滲透破壞形式主要是細(xì)砂層被沖破后，侵蝕逐漸向上游發(fā)展。圖4為上游水頭4 m時(shí)Y方向水力梯度等值線圖，圖中堤基的各層土體從上到下依次為黏土層、細(xì)砂層、砂礫石層，可以看出Y方向水力梯度最大值發(fā)生在管涌口下端的細(xì)砂夾層底面。同時(shí)，由圖5可知管涌口下端細(xì)砂層底面（即A—A截面，見(jiàn)圖4）的最大水力梯度達(dá)到0.93，超過(guò)了表2中計(jì)算所得的細(xì)砂層埋深為0 m時(shí)對(duì)應(yīng)的允許水力梯度0.87。因此，在三層堤基結(jié)構(gòu)中，發(fā)生滲透破壞的臨界上游水頭為4 m。

圖6為細(xì)砂層底部（A—A截面）壓力水頭分布情況，細(xì)砂層的滲透系數(shù)較小，在細(xì)砂層的阻擋下，細(xì)砂層底部承擔(dān)了一定的壓力水頭，在管涌口處細(xì)砂層承擔(dān)的壓力水頭占上游水頭的44%左右，同時(shí)因管涌口處細(xì)砂層上部沒(méi)有泥沙顆粒的阻擋，故管涌口處的細(xì)砂層很容易被沖破，發(fā)生破壞的臨界水頭較小。

2.3 多層堤基管涌發(fā)生的臨界上游水頭

2.3.1 細(xì)砂夾層埋深1 m

多層堤基結(jié)構(gòu)中，當(dāng)細(xì)砂層在砂礫石層中的埋深為1 m、上游水頭為14 m時(shí)Y方向水力梯度的等值線見(jiàn)圖7，圖中堤基從上到下依次為黏土層、砂礫石層、細(xì)砂層、砂礫石層，水力梯度分布情況與圖4類(lèi)似，Y方向水力梯度最大值發(fā)生在管涌口下端的細(xì)砂夾層底面（A—A截面），最大水力梯度達(dá)到2.14（見(jiàn)圖8），超過(guò)了表2中計(jì)算所得的細(xì)砂層埋深為1 m時(shí)對(duì)應(yīng)的允許臨界水力梯度2.00，因此在上游水頭為14 m時(shí)，細(xì)砂層被沖破，發(fā)生管涌破壞。

多層堤基結(jié)構(gòu)中，因細(xì)砂層的分隔，還存在另一種形式的滲透破壞：在細(xì)砂層上部，黏土層與砂礫石層界面上可能發(fā)生沿著接觸面逐漸向上游發(fā)展的侵蝕。這種形式破壞發(fā)生的標(biāo)志是管涌口處細(xì)顆粒發(fā)生豎向運(yùn)動(dòng)并被帶出管涌口，是否發(fā)生主要取決于滲透破壞的臨界水力條件是否得到滿足。圖9為管涌口處（C—C截面，見(jiàn)圖7）Y方向水力梯度分布，由圖9可知，管涌口處的水力梯度小于顆粒允許水力梯度，從而認(rèn)定該種破壞并未發(fā)生。因此，對(duì)于夾砂層埋深1 m的多層堤基結(jié)構(gòu)，滲透破壞發(fā)生的標(biāo)志為細(xì)砂夾層被沖破，滲透破壞發(fā)生的臨界上游水頭為14 m。

2.3.2 細(xì)砂夾層埋深2 m

多層堤基結(jié)構(gòu)中，當(dāng)細(xì)砂層在砂礫石層中的埋深為2 m時(shí)，模擬結(jié)果與細(xì)砂層埋深為1 m時(shí)類(lèi)似。上游水頭為25 m時(shí)，細(xì)砂層底部Y方向的水力梯度分布情況見(jiàn)圖10，管涌口處Y方向水力梯度分布情況見(jiàn)圖11。

由圖10、圖11可知，Y方向水力梯度最大值發(fā)生在管涌口下端的細(xì)砂夾層底面，最大水力梯度達(dá)到3.22，超過(guò)了表2中計(jì)算所得的細(xì)砂層沖破的允許水力梯度3.12;管涌口處的水力梯度最大值發(fā)生在管涌口邊界處，最大為0.366，超過(guò)了顆粒發(fā)生豎向運(yùn)動(dòng)的允許水力梯度，然而管涌口邊界處超過(guò)允許水力梯度的區(qū)域很小，少量顆粒被帶出管涌口。對(duì)于多層堤基，細(xì)砂層埋深為2 m時(shí)，上游水頭升高到25 m，細(xì)砂夾層被沖破，同時(shí)管涌口邊界少量顆粒被帶出管涌口，兩種滲透破壞的臨界條件同時(shí)得到滿足。因此，上游水頭25 m為夾砂層埋深2 m的多層堤基結(jié)構(gòu)發(fā)生滲透破壞的臨界上游水頭。

2.3.3 細(xì)砂夾層埋深5 m

當(dāng)多層堤基中細(xì)砂層埋深增大至5 m、上游水頭為21 m時(shí)，由細(xì)砂層底部的水力梯度分布（見(jiàn)圖12）可知，管涌口下端夾砂層底面的水力梯度最大，達(dá)到1.87，小于表2中計(jì)算所得的細(xì)砂層埋深為5 m時(shí)對(duì)應(yīng)的允許水力梯度6.48，因此細(xì)砂層未被沖破。圖13為管涌口處的水力梯度分布情況，管涌口處的水力梯度最大值為0.395，超過(guò)了顆粒發(fā)生豎向運(yùn)動(dòng)的允許水力梯度，顆粒流失的區(qū)域?yàn)閄=84.0～84.5 m。因此，在細(xì)砂層埋深為5 m的多層堤基中，發(fā)生管涌的臨界上游水頭為21 m，且破壞主要沿著堤基黏土層與砂礫石層界面，細(xì)砂夾層未被沖破。

由上述模擬分析可知，在多層堤基中夾砂層埋深存在一個(gè)臨界埋深值，當(dāng)夾砂層埋深小于臨界埋深時(shí)，堤基主要發(fā)生豎向流土破壞;而夾砂層埋深超過(guò)臨界埋深時(shí)，堤基主要發(fā)生沿黏土層與砂礫石層界面的顆粒流失。模擬結(jié)果表明，夾砂層埋深2 m可以看作多層堤基中細(xì)砂夾層的臨界埋深。

此外，對(duì)于不同的堤基結(jié)構(gòu)，發(fā)生滲透破壞的形式不同，且發(fā)生滲透破壞的臨界上游水頭有很大差異。表3總結(jié)了不同堤基結(jié)構(gòu)發(fā)生管涌的臨界上游水頭，可知三層堤基結(jié)構(gòu)（細(xì)砂層埋深為0 m）發(fā)生管涌的臨界上游水頭最小，為4 m;對(duì)于多層堤基結(jié)構(gòu)，隨著細(xì)砂夾層埋深的增大，發(fā)生管涌的臨界上游水頭逐漸增大，直到細(xì)砂層埋深為5 m時(shí)，臨界上游水頭有所減小。對(duì)于含有細(xì)砂夾層，且?jiàn)A砂層埋深不超過(guò)2 m的堤基，發(fā)生滲透破壞的標(biāo)志是細(xì)砂夾層被沖破，發(fā)生豎向流土破壞;當(dāng)夾砂層埋深超過(guò)2 m時(shí)，發(fā)生滲透破壞的標(biāo)志是管涌口處夾砂層上部的細(xì)顆粒發(fā)生豎向運(yùn)動(dòng)并被帶出，破壞沿著黏土層與砂礫石層界面發(fā)展，破壞與雙層堤基類(lèi)似。

隨著細(xì)砂層埋深的增大，管涌破壞發(fā)生的臨界水頭會(huì)不斷增大，形成的管涌通道深度也不斷增加。當(dāng)細(xì)砂埋深達(dá)到一定深度時(shí)，細(xì)砂上部砂礫石層的滲水量增加，增加到足以將砂礫石表面的顆粒帶出，此時(shí)管涌破壞的過(guò)程與雙層堤基類(lèi)似，通道主要形成在黏土下部的砂礫石表面，臨界水頭將會(huì)有所減小，然而由于細(xì)砂夾層的存在，上部砂礫石層中的滲水量仍小于雙層堤基中砂礫石層中的滲水量，因此臨界水頭仍然大于雙層堤基中發(fā)生管涌的臨界上游水頭。

為了進(jìn)一步研究細(xì)砂夾層對(duì)上游水頭的阻隔作用，對(duì)相同水位、不同細(xì)砂層埋深下堤基內(nèi)部的壓力水頭以及管涌口處的流量進(jìn)行模擬分析。上游水頭為10 m時(shí)，雙層堤基管涌口處的流量為2.2×10-3 m3/s;當(dāng)細(xì)砂夾層埋深分別為0、1、2、5 m時(shí)，管涌口處的流量分別為9.4×10-4、1.5×10-3、1.7×10-3、1.9×10-3 m3/s，且細(xì)砂夾層底部壓力水頭分別占上游水頭的43.6%、30.4%、27.8%、26.2%。隨著細(xì)砂層埋深的增大，管涌口處的流量逐漸增大，細(xì)砂層底部壓力水頭占上游水頭的百分比逐漸減小，也就是說(shuō)隨著細(xì)砂層埋深的增大，細(xì)砂層上部砂礫石層中的流量逐漸增大，但仍小于雙層堤基中砂礫石層的流量。

3 算例分析

王霜等[12]通過(guò)自行設(shè)計(jì)的室內(nèi)試驗(yàn)裝置對(duì)雙層、三層和多層結(jié)構(gòu)堤基進(jìn)行了管涌發(fā)展的砂槽模型試驗(yàn)，試驗(yàn)裝置見(jiàn)圖14。模型槽長(zhǎng)100 cm，寬、高均為30 cm，模型槽左側(cè)為一個(gè)長(zhǎng)20 cm的進(jìn)水室，進(jìn)水室與砂槽之間由透水板隔開(kāi)，砂槽上游側(cè)設(shè)置7 cm長(zhǎng)的礫石層，使水流均勻地流入試樣內(nèi)。模型槽頂面為剛性玻璃蓋板，蓋板下面設(shè)置一定厚度的黏土層，模擬弱透水的上覆層。在玻璃蓋板上，距進(jìn)水室60 cm處預(yù)設(shè)直徑為5 cm的出水口，模擬上覆層被沖破后所形成的管涌口。模型槽內(nèi)部布置多根測(cè)壓管，測(cè)量試驗(yàn)過(guò)程中試樣內(nèi)部孔隙水壓力的變化情況。3組試驗(yàn)涉及3種不同的土層結(jié)構(gòu)，各試驗(yàn)方案的條件和試驗(yàn)結(jié)果見(jiàn)表4。管涌破壞的臨界水力梯度為試樣開(kāi)始涌砂時(shí)的水平表觀水力梯度i，i=H/L ，其中：H為上游水頭，L為管涌口到進(jìn)水室的水平距離。

試驗(yàn)結(jié)果表明，不同結(jié)構(gòu)的堤基滲透變形的臨界水力梯度有很大差異。雙層堤基結(jié)構(gòu)，管涌破壞主要沿著黏土層下部砂礫層表面發(fā)展，管涌發(fā)生的標(biāo)志為管涌口下部的砂礫石沿豎向被帶出管涌口，管涌發(fā)生的臨界水力梯度為0.57;細(xì)砂夾層在砂礫石表面的三層堤基，水壓力主要集中在細(xì)砂夾層底部，管涌發(fā)生的1.22標(biāo)志為細(xì)砂夾層被沖破，且管涌口處的細(xì)砂夾層上部沒(méi)有覆蓋，故管涌發(fā)生的臨界水力梯度較小，為0.39;細(xì)砂夾層埋深為1 cm的多層堤基，管涌發(fā)生的標(biāo)志同樣為細(xì)砂夾層被沖破，管涌發(fā)生的臨界水力梯度為1.22。該試驗(yàn)結(jié)果與本研究數(shù)值模擬結(jié)果完全吻合，當(dāng)?shù)袒写嬖诩?xì)砂夾層，且細(xì)砂夾層在砂礫石內(nèi)有一定埋深時(shí)，隨著細(xì)砂夾層埋深的增大，管涌發(fā)生的臨界水力梯度逐漸增大。文獻(xiàn)[12]推斷，當(dāng)細(xì)砂夾層埋深達(dá)到一定深度，即達(dá)到“臨界埋深”時(shí)，細(xì)砂夾層上部砂礫石層的水量將會(huì)很大，大到足以將砂礫石層表面的顆粒帶出，此時(shí)管涌破壞的過(guò)程與雙層堤基類(lèi)似，臨界水力梯度將會(huì)有所減小，通道主要形成在黏土下部的砂礫石層表面。本研究數(shù)值模擬結(jié)果驗(yàn)證了該推斷的正確性與合理性。

4 結(jié) 論

不同結(jié)構(gòu)的堤基，管涌發(fā)生破壞的機(jī)理不同。雙層堤基管涌破壞發(fā)生的標(biāo)志是黏土層與砂礫石界面土體發(fā)生侵蝕。三層堤基管涌破壞發(fā)生的標(biāo)志是細(xì)砂夾層被沖破。多層堤基中，細(xì)砂夾層的存在將堤基的砂礫石層分成了夾層上部與夾層下部?jī)刹糠?，堤基發(fā)生管涌時(shí)存在兩種形式的破壞：夾層上部，在黏土層與砂礫石層界面上發(fā)生侵蝕破壞，破壞以管涌口下部顆粒發(fā)生豎向運(yùn)動(dòng)為標(biāo)志;夾層下部，細(xì)砂夾層被沖破，發(fā)生豎向流土破壞。

當(dāng)?shù)袒鶅?nèi)部存在細(xì)砂夾層時(shí)，隨著細(xì)砂夾層埋深的增大，發(fā)生管涌的臨界上游水頭逐漸增大，管涌通道深度會(huì)不斷增加，發(fā)生管涌破壞的標(biāo)志是細(xì)砂夾層被沖破，發(fā)生豎向流土破壞。當(dāng)細(xì)砂夾層埋深達(dá)到一定深度時(shí)，細(xì)砂上部砂礫石層的水量增加，增加到足以將砂礫石層表面的顆粒帶出，此時(shí)管涌破壞過(guò)程與雙層堤基類(lèi)似，通道主要形成在黏土下部的砂礫石層表面，臨界水頭將會(huì)有所減小，然而因細(xì)砂夾層的存在，上部砂礫石層中的水量仍小于雙層堤基中砂礫石層中的水量，臨界水頭仍然大于雙層堤基中發(fā)生管涌的臨界上游水頭。

參考文獻(xiàn)：

[1] 劉川順，彭幼平，李莫.長(zhǎng)江堤防地基類(lèi)型與防滲方案[J].水文地質(zhì)工程地質(zhì)，2001，28（6）：20-22.

[2] 馬貴生，王造根.長(zhǎng)江重要堤防工程地質(zhì)勘察[J].人民長(zhǎng)江，2002，33（8）：60-62.

[3] 李廣誠(chéng)，司富安，白曉民.中國(guó)堤防工程地質(zhì)[M].北京：中國(guó)水利水電出版社，2005：9-12.

[4] 殷建華.土堤管涌區(qū)滲流的有限元模擬[J].巖石力學(xué)與工程學(xué)報(bào)，1998，17（6）：679-686.

[5] 張家發(fā)，王滿星，丁金華.典型條件下堤身堤基滲流規(guī)律分析[J].長(zhǎng)江科學(xué)院院報(bào)，2000，17（5）：23-27.

[6] 周曉杰，介玉新，李廣信.基于滲流和管流耦合的管涌數(shù)值模擬[J].巖土力學(xué)，2009，30（10）：3154-3158.

[7] ZHOU X J， JIE Y X， LI G X. Numerical Simulation of the Developing Course of Piping[J].Computers and Geotechnics， 2012，44：104-108.

[8] 劉杰.土的滲透穩(wěn)定與滲流控制[M].北京：水利電力出版社，1992：32-33.

[9] 劉忠玉，樂(lè)金朝，苗天德.無(wú)黏性土中管涌的毛管模型及其應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào)，2004，23（22）：3871-3876.

[10] SKEMPTON A W， BROGAN J M. Experiments on Piping in Sandy Gravels[J]. Geotechnique， 1994，44（3）：449-460.

[11] KOVACS G. Seepage Hydraulics[M]. New York： Elsevier Scientific Publishing Company， 1981：53-55.

[12] 王霜，陳建生，黃德文，等.土層結(jié)構(gòu)對(duì)管涌發(fā)展影響的試驗(yàn)研究[J].巖土工程學(xué)報(bào)，2013，35（12）：2334-2341.

【責(zé)任編輯 許立新】