SEEPV3.0

或就火豫淤嘲水痘簪微

二00四年一月

1滲流基本理論

1.1水工滲流的危害及滲流分析計算的任務

流體在多孔介質中的運動稱為滲流。水是最為常見的流體，水利水電工程中由于廣泛建

造堤、壩、圍堰、水閘等擋水建筑物形成了水頭差，這些建筑物或其地基一般是透水口勺多孔

介質，因此水工滲流現象十分普遍。

水工滲流導致多方面的危害。滲流導致水庫、渠道水量損失；滲流使堤壩、圍堰土體飽

和，減少壩體的有效容重和抗剪強度，也許導致壩坡失穩(wěn)：建筑物地基滲流對建筑物底部產

生揚壓力，也不利于建筑物的穩(wěn)定：滲流流速過大時，還也許導致壩體或建筑物地基日勺土體

顆粒流失，發(fā)生滲透變形，從而使堤壩倒塌或建筑物滑移、傾覆：水庫滲流還也許引起卜游

地下水位升高，導致農田冷浸漬害、鹽堿化，使作物減產；攔污壩滲流導致地下水環(huán)境污染。

水工滲流分析計算的任務就是要研究水在滲流區(qū)域H勺滲流流速、流星、水頭分布及浸

潤線等，從而為采用合理的滲流控制措施提供根據，以防止或減緩滲流危害。

1?2達西定律

19世紀50年代，法國工程師亨利?達西（H.Darcy）通過對裝在圓筒中的均質砂土進

行滲透試驗發(fā)現，通過兩個滲流斷面間日勺平均滲流流速，正比于兩斷面間口勺水頭差Ah,反

比于滲徑長度L,且與土粒構造及流體性質有關。這就是著名的達西定律，可用公式體現為:

(1.2.1)

式中h-測壓管水頭,z為位置高度，p為壓強，丫為水冰、J容垂。由于

丫2g

滲流H勺流速?般很小,流速水頭?？珊鲆?，故h=z+20

2gY

k一反應土粒構造及流體性質的系數，即滲透系數，對于某一詳細的流體(例如水)而

言，k值僅與土粒構造有關。

J一滲透坡降，J=—.

ds

式中的負號“-”表達水總是流向水頭減小的方向。

應當注意，達西定律中口勺流速是全斷面上的平均流速V,而不是上體孔隙中H勺流速V,

這兩種流速存在如下關系：

v=nv'(1.2.2)

式中n為體積孔隙率，可見達西流速不大于上體孔隙中的流速.

還應注意，達西定律只能合用于層流狀態(tài)的滲流運動。在水利工程中，除了堆石壩、堆

石排水體等大孔隙介質中的滲流為流之外，絕大多數滲流都屬于層流，達西定律都可合用。

對于非層流滲流，其流動規(guī)律可用如下公式形式體現：

1

v=k,Jm()

上式中當m=l時，為層流滲流；當m=2時，為完全紊流滲流：當l〈m<2時，為層流到紊流

的過渡區(qū)。

將式()等號兩邊向X、？、z軸投影，便得到空間直角坐標系中的達西公式：

滲流運動持續(xù)性方程

y

圖1.2.1微元體進出流量示意

從滲流場中取某一微分單元體(如圖)，其體積V=4rdydz,流入左側面的水體質量的

速率為流出右側面的水體質量H勺速率為(pv.+Ji，Mz)c￡r小，，則左右面進出流

*oz"

量之差為

^-pvxdxdydz

ox

同樣，對于前背面和上卜面作流進流出的流入量計算，最終累加各凈有流入量，得到微

元體內總得流入量為

da

--pv,+—

辦OZ

將上式展開為

+生+生

V-

/dz糕+襟+吟上

式中后一項與前一項相比，小得可以忽視，故改寫成

o

上式即為水體質量在微元體內積累的速率，根據質量守恒原理，它應等于微元體內水體

質量V隨時間口勺變化速率:

cM_d(npV)=3絲+即包+亞農

dtdtdtdtdt

式中，n為土體的孔隙率，P為水的密度，V為微元體的體積。式（）右邊三項分別代表土

體骨架、孔隙體積及流體密度H勺變化速率。

引入彈性壓縮理論，可導出（）式中得等價體現為:

6M,/,,\1Zdh

—=P-^+,；/y）v-

式（）中a——土顆粒骨架的壓縮性（即壓縮模量）

B一—水的壓縮性

h——滲流測壓管水頭（m）,h=z+E

根據質量守恒原理，式（）與式（）相等，則得到

(d\\8八dv.

--+--+―-=pg(a+邛虺

(及dy3z)

假定水和土不可壓縮時，上式變?yōu)?

dxdydz

式（）為不可壓縮流體在剛體介質中流動的持續(xù)性方程。

滲流微分方程

將達西公式代入滲流持續(xù)性方程（）可得：

=汕+〃畔=渭

dx)dz\'dzJ

式（）為非穩(wěn)定滲流微分方程的?般形式，既適合于承壓含水層，也適合于無壓滲流。式中

S,=必（。+〃夕）稱為單位儲存量（尺度為1/L）,其含義是：單位體積得飽和土體在水頭

下降1m時，由于土體壓縮（pga）和水體膨脹（pgnp）所釋放出來得儲存水量。

當介質為均質各向同性（即幻=勺.=￡=攵）時，式（）變?yōu)?/p>

222

a/?dhdhsvdh

dx2*dy2+dz2kdt

當假定水和土為不可壓縮時，式（）、（）分別變成

式（）、（）分別為各向異性、各向同性時的三維穩(wěn)定滲流微分方程。

對于堤壩滲流及建筑物地基滲流，在許多狀況下可簡化為垂直剖面上的二維滲流問

題，則式（）、（）可變成

a（dh\d（,dh\八

—kt<—+—k.—=0（）

dx\Jdz\'dz）

d2hd2h.、，

——+———0（1.2.15）

dx~dz~

有自由面變動的滲流微分方程

對于土壩（或堤防）在水庫（或河道）水位下降時H勺非穩(wěn)定滲流狀況，自由面下降

引起的土體壓縮或彈性釋放水量與自由面下降時所排出的水量相比很小，故可令式（）中的

&=0,因此，有自由面日勺三維非穩(wěn)定滲流微分方程為

a（.%d（劭）ardh\,、

—k——+—kt—+—k.——=0n（）

dx\xdx）dy\vdy）V'dz）

其二維非穩(wěn)定滲流微分方程為

上面二式雖然在形式上與穩(wěn)定滲流方程式（）、（）完全相似，但結合自由面變動的邊界

條件所得到H勺水頭分布是空間座標與時間H勺函數，而不象穩(wěn)定滲流方程式I燈解答，只是空間

座標的函數。

水工滲流問題的定解條件

發(fā)生在有限空間流場內H勺滲流運動，不僅受滲流微分方程支配，也受流場邊界條件和初

始滲流狀態(tài)條件支配。下面以土壩非穩(wěn)定滲流狀態(tài)簡介水工滲流定解條件的幾種類型。

圖土壩非穩(wěn)定滲流定解條件

初始條件：〃|go=%（x,z"）（僅對非穩(wěn)定滲流）（）

邊界條件：

⑴水頭邊界胃口=/（兀2,力（）

在圖1.2.2中，水頭邊界包括：

上游已知水頭邊界：hABC=H1（t）

下游已知水頭邊界：hKFC=H2

下游坡滲出段已知水頭邊界：hDE=2

⑵流量邊界：夕二%“一|「2=/2（X，Z/）（1.2.20）

on“

在圖中，自由而AD下降時，由自由而流進的單寬流量為：

ah

q-Li—cos。（A）

dt

式中。一自由面的法線與鉛直線的夾角

口一給水度，即單位體積土體在飽和含水狀況下，當自由面下降1m后排出I向空隙水

量，乂稱土體的排水空隙率或有效空隙率。

又因

q=-k?—=-k?—cos。（B）

"dn“dz

由式（A）、（B）,自由面流量邊界條件改寫成：

4ohdh

-----=----（C）

kndtdz

當壩體土各向異性時，幻取心值。

因滲流自由面與大氣相通，壓力水頭為零，故在自由面上還應滿足：

h=z"0

2SEEPV3.0的功能和特色

（1）合用于多種工況（如水位固定或上游水位快降）、復雜邊界（包括多種防滲、排水設施）

和復雜土層條件下的：堤防、土石壩滲流計算；各類不透水材料壩（如混凝土壩和漿砌石壩）、

水閘、泵房口勺地基滲流計算。

（2）采用有限元法計算，可得到滲流量、浸潤線、等水頭線等計算成果。

（3）采用可視化界面建立原始數據文獻,也可在記事本中建立原始數據文獻。

（4）計算成果既可打印成表格，也可白動繪制出浸潤線、等水頭線圖。

3SEEPV3.0解題要點

3.1坐標系的選用及計算范圍確實定

如圖3.1.1所示，水平坐標x取向右為正，豎直坐標與當地高程體系一致。

對于上游入滲段長度LI及下游滲出段長度L2,一?般取0.，~1倍的水頭（或取實測值）。

圖3.1.1

3.2單元剖分

a）單元自動剖分的分組、分塊與結點編號

根據滲流場的幾何形狀、土層分布等，將整個滲流場劃分為某些四邊形的子域，如

圖321所示。這些子域稱為塊，為使結點優(yōu)化并滿足自由面調整迭代H勺需要，將所劃分的

塊又提成若干組，每組內由若干塊構成。詳細規(guī)定；每組內各坎首尾相接，上下貫穿整個滲

流場：每組內的塊由上而下逐塊編號，且從左到右一組一組地進行，但塊的結點標號是任意

的（提議按照從左到右，從上到下H勺次序）。圖3.2.1中由結點3,7,11,12,15連成的折

線為事先假定的浸潤線。

VH130⑵訂

6)⑶/⑶

⑶1418回12(5)15

———.[3]

⑶調期因小8（心乙、二也

1w<4)9:⑶13⑤16(4)21

③（*國盤》⑶[Q⑶

2⑷610(2)14<5)17(4)20W22

圖3.2.1

b）四條邊剖分點布置

本軟件進行滲流有限元計完的基本微小單元為三角形單元,為此還需將上述四邊形塊單

元剖分為三角形單元.如圖3.2.2所示，從滲流場中取出任意一塊，圖中nl,n2,n3,n4

分別表達對應邊上剖分點數（剖分點數決定三角形單元打勺數目或大?。?，且按順時針方向編

號，nl與n3、n2與n4分別為兩組對邊，兩對邊上的剖分點數可相等，也可不相等，但為

防止出錯，最佳使各對邊H勺點數相等或相差小某些。剖分點的疏密應按水力條件來確定，在

下列部位的劑分單元應加密：水力坡降大的部位、壩軸線附近、截水槽及上游鋪蓋進水處，

下游出口及排水溝附近。此外，剖分點數的確定還必須滿足如下三點：

①三角形單元總數＜500：

②一〃JV〃2，%|〃4一〃21y3：

③浸潤線變動范圍內々=〃3

O剖分次序

剖分次序依次為：劃分塊單元、分組并編組號、編塊號并查對塊單元的結點數、編結點

號、在各塊邊上布置剖分點數。

同一塊內的滲誘系數必先相似，伸同一滲誘系數土層可劃分為若干塊。

每個塊都必須是四邊形塊。在滲流區(qū)域邊界上有時不可防止地出現三角形塊，此時在該

塊的邊界邊上任取一點作為結點，這樣?來，三角形塊可以當;乍是蛻化的四邊形塊。

4SEEPV3.0的操作環(huán)節(jié)與文獻內容簡介

4.1DQB文獻的產生

工程有關參數:

(1)填入“分塊組數”和“土層種類數”；

(2)等勢線間距有5%和10%兩個數可供選擇，分別表達在滲流場分布圖中相鄰兩條等

勢線之間的水頭勢能相差5%和10%；

(3)滲流狀況：假如計算非穩(wěn)定滲流，需要填寫水位降落經歷口勺時間、初始水位和降后

水位：假如計和穩(wěn)定滲流，需要選擇與否計冕滲流量。

視圖有關參數：

(1)“最大坐標”和“最小坐標”：是您在繪圖過程中所能輸入的最大最小坐標值，為了

以便繪圖，提議將斷面的最大X、Y坐標值分別加上30,將斷面的最小X、Y坐標

值分別減去30后輸入：

(2)“l(fā)m=()P”表達繪圖時一米長度所顯示日勺象素個數，假如圖形中有比較短的線條，

可填入較大值，例如10：假如圖形中線條都較長，提議采用默認值4；

(3)”點的選用半徑”默認為3個像素，設置大些可以更好地捕捉點，但過大會影響

線條的選用，提議采用默認值；“字體大小”一般默認為10。

第二步：單擊“工程及視圖參數界面”中確煙按鈕，進入圖形繪制界面(如圖4.1.2)：

1攵141CI攵1必〃八―15而

默認狀態(tài)a下進行圖形繪制：

(1)圖形繪制：在右下角“LineFrom”的背面直接輸入節(jié)點坐標并按回車鍵:

(2)假如節(jié)點坐標有錯誤，可以點擊工進行修改；

(3)點擊聞，用左鍵選中線條，再單擊右鍵可以刪除線條

切換到狀態(tài)里進行單元剖分；

(1)節(jié)點編號：按編號次序直接點擊該節(jié)點，假如不能選口該節(jié)點，可以點擊_L加大

“點口勺選用半徑”：

(2)剖分控制：先選中塊的四條邊(假如不能選中該邊，可以點擊減小”點的選用

半徑”，或者將圖形放大)，然后點擊■進行剖分填寫，最終按塊的編號次序進行

其他塊的剖分。

數據填寫和文獻保留:

(1)填寫土層滲透系數：點擊■輸入每層土的水平和垂直滲透系數標、Kv(cm/s),

注意：每輸入一層土的滲透系數，都要點擊邈按鈕進行數據保留：

(2)剖分塊信息修改：假如剖分信息有錯誤，點擊旦進行修改(先選中Reg下某一塊

日勺編號，再單擊酗按鈕)：

(3)設置有關參數：點擊上打開有關參數設置對話框，檢查“分組數”、“等勢線百分

數間距”和“滲流狀況”與否對的，注意：當計算工況為穩(wěn)定滲流且計算滲流量的

時，需填寫“滲流斷面穿過的塊號”;

(4)保留dqb.dat文獻:點擊正彈出dqb文獻保留對話框,選擇保留途徑后點擊保留⑤)

按鈕，注意：計算時文獻名必須為dqb.dat：

(5)生成斷面輪廓圖：假如需要斷面輪廓圖.有接在圖中點擊鼠標片彼選中所有要保留

gjju

的輪廓線后，點擊一彈出保留對話框，填入文獻名并保留。

4.2運行計算

點擊合一按鈕，彈出計和操作界面(如圖4.2)：

函4個：.|?8隰/任田而

(1)點擊|查看DQB］按鈕,檢查產生的DQB文獻與否對H勺(可以參照闡明書中4.3DQB

.DAT文獻的內容和格式)；

(2)點擊［剖分單元|按鈕,生成滲流計算需要的原始數據文獻G-DAT和剖分圖NETO.BLN：

(3)點擊|顯示剖分圖|按鈕,進行單元剖分同步生成網格圖NET.BLN；

(4)點擊|編輯數據|按鈕,打開計算原始數據文獻G-DAT,注意：在該文獻末尾另起一行

輸入0并保留；

(5)點擊旺回按鈕，生成成果文獻SEEP?OCT和滲流場圖形文獻GS?BLN；

(6)點擊|等勢線圖|啟動SIRFER,繪制滲流場分布圖：

a、單擊MAP菜單中晌LOADBASEMAP出現IMPORTFILE視窗，雙擊BOUND?BLN.

單擊OK,隨即出現土壩橫斷面圖；

b、單擊MAP菜單中的LOADBASEMAP出現IMPORTFILE視窗，雙擊GS?BLN,單

擊OK,隨即出現滲流場圖形；

c、將土壩橫斷面圖與滲流場圖形進行登加：用鼠標將滲流場圖形拖至堤防橫斷面

圖，用鼠標同步框選堤防橫斷面圖和單個滲流場圖形，單擊MAP出現MAP菜單，

單擊MAP菜單中的OVERLAYMAPS即完畢圖形疊加；

d、進行文字編輯、坐標軸設定等，最終得到完美口勺滲流場分布圖。

（7）若想查看計算成果SEEP?OUT則點擊計算匯報|按鈕。

4.3DQB.DAT文獻的內容和格式

產生的數據文獻DQB.DAT內容格式如下：

第一行：SEEPLSEEP2,SEEP3

第二行：ZU,KUALDIAN,BIAN,TUCENG,SLDM,SLK,SD,BFS

第三行：KXX（1：KUAI,1：9）

第四行：ZB（1：DIAN,I：2）

第五行：SLCS（1：TUCENG,1：3）

第六行：BJXX（1：BIAN,I：3）

第七行：SLKH（1：SLK）

第八行：SLDMXX（1：SLDM,1：3）

第九行:SJSW（1：2,I：2）

數聽闡明：

（1）SEEPI—鑒別與否保留滲流繪圖數據，'Y'為肯定，'N'為否認。

SEEP2—鑒別與否自動剖分數據，'Y'為肯定，'N'為否認。

SEEP3一設定為，N'。

（2）ZU一分塊組數（計算區(qū)域內剖分的塊依次從左到右分為若干組，同一組間應為上下關

系，即第I塊口勺下邊（由）與第1+1塊的上邊（m）應為同一邊，每組內由上而下次序編號）;

KUAI一四邊形塊H勺數目：

DIAN一塊角點數：

BIAN一含各類已知邊界的邊數;

TUCENG一不一樣滲透系數的土層種類數：

SLDM一需要計算滲流量的滲流斷面?zhèn)€數，假如是零表達不計算滲流量；

SLK一需要計算滲流量的斷面所穿過的J總計塊數，假如是零表達不計算滲流量；

SD—計算非穩(wěn)定滲流時水庫水位降落H勺時段數，假如是零表達計算穩(wěn)定滲流；

BFS一計算等勢線的百分數詞距值，默認為I0o

(3)KXX(1:ME,1:9)塊信息數組

每塊9個分量，按塊號從小到大次序填寫。第1—4分量表達塊從左上角起按順時針方

向的四個角點點號；第5—8分量依次表達塊從上邊起按順時針方向4條邊口勺剖分點數：第

9分量為塊的滲透系數區(qū)號。

(4)ZB(I：DIAN,1：2)塊口勺結點坐標數組

按塊結點序號次序表達每一結點破JX,Y坐標。

(5)SLCS(I：TUCENG,I：3)

依次為每層土的水平滲透系數Kx(單位：cm/s),垂直滲透系數Ky(單位：cm/s),給

水度二°117四(該式中滲透系數k的單位為：m/day).

(6)BJXX(1:BN,1:3)已知之界信息數組

每條已知邊界有三個分量,第一分量為邊界邊所在塊的塊號：第二分量為該邊界邊在塊

中的邊序號；第三分量為邊的邊界類型號，規(guī)定為：

1一一滲流自由面邊界(斜、心墻壩僅為上部自由而)

3一—半有東流邊界

5——斜、心墻壩下游壩殼自由面邊界

6一—均質上壩下游壩坡或排水體滲出段邊界

7一一上游已知水頭邊界

8一—下游已知水頭邊界

9——斜、心墻土壩斜墻心墻后坡滲出段邊界

（7）SLKH（1：SLK）

需計算滲流量的過流斷面上的塊號。

（8）SLDMXX（1：SLDM,1：3）滲流斷面歐J塊單元信息

第?分量為某過流斷面所通過的塊數；第二分量為該過流斷面與塊邊界相交的信息，

水工滲流計算一般為1：第三分量為該過流斷面在塊中位置的節(jié)點鏈序號，一般為1或2。

（9）SJSW（1：2,1：2）起降時間（天）、水位；降后時間（天）、水位。

4.4G.DAT數據內容

本軟件可將四邊形塊單元自動剖分為三角形單元，司步由DQB.DAT自動產生

G.DA1;G.DAT是本軟件進行滲流計算的原始數據文獻，其數據內容格式如卜.：

第一行：SEEPI,SEEP2,SEEP3

第二行：NE,NN,NNU,NND,NNM1,NNM2,NNF,TICENG,AA,AA1,AA2,

BB1,BB,BB2,NNP,NQ,SD,NW,NG,NNM3,LU,LD,MA

第三行；Hmax，HI,H2,H3,EPS,EPSI,DT,YAA，XAA

aooooo

第十四行：SW（1：SD）

第十五行：SJJG（I：SD）

第十六行：SLCS（1：TUCENG,1：3）

第十七行：QTE（I：NQ）

第十八行：QEN（1：QE）

解釋：

（1）SEEP1,SEEP2,SEEP3-RDQB.DAT

（2）NE一三角形單元總數

aooooo

<3）He一上游最高水位

Hl—上游水位H2一下游水位

H3一井水位

EPS一自由面迭代容許誤差（0.1~0.001）

EPS1一井附近迭代容許誤差（0.5~0.005）

DT—水位變化總歷時（單位：天）

。0。。。。。

（14）SW（1：SD）

各時段末上游水位。

（15）SJJG（1：SD）

各時段時間間隔天數

（16）SLCS（I：TUCENG,I：3）

依次填寫每層土的水平滲透系數K,（單位：m/day）：垂直滲透系數Kv（單位：m/day）；

給水度〃，〃=0/17探（式中k為滲透系數,單位：m/day）.

（17）QTE（I：NQ）

計算滲流量斷面日勺單元總數

(18)QEN(1：QE)

計算滲流量斷面的所有單元號。

5算例

例一均質土壩穩(wěn)定滲流計算

某均質土壩橫斷面尺寸、上下游水位及滲透系數如圖5.1所示，規(guī)定計算穩(wěn)定滲流場

解題環(huán)節(jié)如下:

(1)運行CCED.EXE文獻，彈出工程及視圖有關參數界面，在“分塊組數”中填入5,在

“土層種類數”中填入2,選擇“穩(wěn)定滲流”，并選擇“計算滲流量”：

(2)單擊|視圖有關參?！白畲笞鴺恕敝蠿填寫170,Y填寫100,“最小坐標”中X填

寫-30,Y填寫10,其他的按照默認值:

(3)單擊雁|按鈕,進入圖形繪制界面，參照“4.1DQB文獻的產生”中詳細的操作環(huán)節(jié)

建立均質壩橫斷面剖分圖(如圖5.2)：

㈤匚Of居癡出曲百加八㈤

(4)保留dqb.dal文獻，內容如下：

，、.，，丫，，“

58167212010

145232321

367423442

478542421

6910744432

71011842421

912131084842

1013141182821

1315161442421

0.00045.000

0.00040.000

40.00055.000

10.00045.000

10.00040.000

42.00055.000

42.00045.000

42.00()40.00()

67.00065.000

67.00045.000

67.00040.000

121.00048.000

118.00045.000

118.00040.000

140.00045.000

140.00040.000

0.0.0.

0.0.0.

117

247

217

417

61I

628

818

67

212

（5）點擊圖形繪制界面中"勺R按鈕，彈出計算操作界面：

a、點擊|杳看DQB|按鈕,檢杳產生的DQB文獻與否對的（可以參照闡明竹中4.3DQB.DAT文

獻的內容和格式）：

b、點擊|剖分單元|按鈕,生成滲流計算需要H勺原始數據文獻G-DAT和剖分圖NETO.BLN；

c、點擊|顯示剖分圖|按鈕,進行剖分同步生成網格圖NET.BLN；

（1、點擊解輯數據|按鈕,打開L算原始數據文獻G-DAT,注意：在該文獻末尾另起一行輸入

0并保留；

e、點擊肝回按鈕,生成成果文獻SEEP-OUT和滲流場圖形文獻GS-BLN：

八點擊|等勢線圖|啟動SURFER,察看滲流場分布圖：

選擇MAP菜單中的LOADBASEMAP出現IMPORTFILE視窗,雙擊BOUND-BLN,單擊OK,

隨即出現上壩橫斷面圖：

選擇MAP菜單中的LOADBASEMAP出現IMPORTFILE視窗，雙擊GS?BLN,單擊OK,

隨即出現滲流場圖形：

將土壩橫斷面圖與滲流場圖形進行疊加：用鼠標將滲流場圖形拖至堤防橫斷面圖，同步

框選堤防橫斷面圖和單個滲流場圖形，選擇MAP菜單中的OVERLAYMAPS即完畢圖形疊加；

進行文字編輯、坐標軸設定等，最終得到完美的滲流場分布圖（如圖5.3）。

圖5.3均質壩滲流場分布圖

打開計算匯報（SEEP-OUT）,得到單位壩長口勺日滲漏流量1=0.0313（m3/day/m）＜,

例二堤防非穩(wěn)定滲流計算

堤防、堤基橫斷面尺寸及滲透系數如圖5,4所示，下游水位為31.0m,河水位由38.0m經

4.5天降至36.0m,乂經3.5天降至34.0m。規(guī)定計算河水位降落過程中H勺堤防滲流場。

圖5.4堤防橫斷面

解題環(huán)節(jié)可以參照例一，建立堤防橫斷面剖分圖（如圖5.5）:

圖5.5堤防橫斷面剖分圖

產生本問題日勺DQB-DAT文獻內容如下:

Y''Y''N'

35127200210

156245452

37841331331

48951331531

591061551552

911121045452

0.00030.000

0.00020.000

30.00038.000

20.()0034.000

10.00030.000

10.00020.000

53.00033.500

58.00031.000

60.00030.000

60.00020.000

70.00030.000

70.00020.()00

0.0.