滲流計算操作指南

SEEPV3.0

或就火豫淤嘲水痘簪微

二00四年一月

1滲流基本理論

1.1水工滲流的危害及滲流分析計算的任務(wù)

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

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

介質(zhì)，因此水工滲流現(xiàn)象十分普遍。

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

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

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

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

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

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

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

1?2達(dá)西定律

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

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

比于滲徑長度L,且與土粒構(gòu)造及流體性質(zhì)有關(guān)。這就是著名的達(dá)西定律，可用公式體現(xiàn)為:

(1.2.1)

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

丫2g

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

2gY

k一反應(yīng)土粒構(gòu)造及流體性質(zhì)的系數(shù)，即滲透系數(shù)，對于某一詳細(xì)的流體(例如水)而

言，k值僅與土粒構(gòu)造有關(guān)。

J一滲透坡降，J=—.

ds

式中的負(fù)號“-”表達(dá)水總是流向水頭減小的方向。

應(yīng)當(dāng)注意，達(dá)西定律中口勺流速是全斷面上的平均流速V,而不是上體孔隙中H勺流速V,

這兩種流速存在如下關(guān)系：

v=nv'(1.2.2)

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

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

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

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

1

v=k,Jm()

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

的過渡區(qū)。

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

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

y

圖1.2.1微元體進(jìn)出流量示意

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

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

*oz"

量之差為

^-pvxdxdydz

ox

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

元體內(nèi)總得流入量為

da

--pv,+—

辦OZ

將上式展開為

+生+生

V-

/dz糕+襟+吟上

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

o

上式即為水體質(zhì)量在微元體內(nèi)積累的速率，根據(jù)質(zhì)量守恒原理，它應(yīng)等于微元體內(nèi)水體

質(zhì)量V隨時間口勺變化速率:

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

dtdtdtdtdt

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

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

引入彈性壓縮理論，可導(dǎo)出（）式中得等價體現(xiàn)為:

6M,/,,\1Zdh

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

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

B一—水的壓縮性

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

根據(jù)質(zhì)量守恒原理，式（）與式（）相等，則得到

(d\\8八dv.

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

(及dy3z)

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

dxdydz

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

滲流微分方程

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

=汕+〃畔=渭

dx)dz\'dzJ

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

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

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

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

222

a/?dhdhsvdh

dx2*dy2+dz2kdt

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

式（）、（）分別為各向異性、各向同性時的三維穩(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)定滲流方程式（）、（）完全相似，但結(jié)合自由面變動的邊界

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

座標(biāo)的函數(shù)。

水工滲流問題的定解條件

發(fā)生在有限空間流場內(nèi)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下降時，由自由而流進(jìn)的單寬流量為：

ah

q-Li—cos。（A）

dt

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

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

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

又因

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

"dn“dz

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

4ohdh

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

kndtdz

當(dāng)壩體土各向異性時，幻取心值。

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

h=z"0

2SEEPV3.0的功能和特色

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

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

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

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

（3）采用可視化界面建立原始數(shù)據(jù)文獻(xiàn),也可在記事本中建立原始數(shù)據(jù)文獻(xiàn)。

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

3SEEPV3.0解題要點

3.1坐標(biāo)系的選用及計算范圍確實定

如圖3.1.1所示，水平坐標(biāo)x取向右為正，豎直坐標(biāo)與當(dāng)?shù)馗叱腆w系一致。

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

圖3.1.1

3.2單元剖分

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

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

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

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

流場：每組內(nèi)的塊由上而下逐塊編號，且從左到右一組一組地進(jìn)行，但塊的結(jié)點標(biāo)號是任意

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

線為事先假定的浸潤線。

VH130⑵訂

6)⑶/⑶

⑶1418回12(5)15

———.[3]

⑶調(diào)期因小8（心乙、二也

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

③（*國盤》⑶[Q⑶

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

圖3.2.1

b）四條邊剖分點布置

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

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

分別表達(dá)對應(yīng)邊上剖分點數(shù)（剖分點數(shù)決定三角形單元打勺數(shù)目或大小），且按順時針方向編

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

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

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

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

①三角形單元總數(shù)＜500：

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

③浸潤線變動范圍內(nèi)々=〃3

O剖分次序

剖分次序依次為：劃分塊單元、分組并編組號、編塊號并查對塊單元的結(jié)點數(shù)、編結(jié)點

號、在各塊邊上布置剖分點數(shù)。

同一塊內(nèi)的滲誘系數(shù)必先相似，伸同一滲誘系數(shù)土層可劃分為若干塊。

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

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

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

4.1DQB文獻(xiàn)的產(chǎn)生

工程有關(guān)參數(shù):

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

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

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

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

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

視圖有關(guān)參數(shù)：

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

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

值分別減去30后輸入：

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

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

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

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

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

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

默認(rèn)狀態(tài)a下進(jìn)行圖形繪制：

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

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

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

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

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

“點口勺選用半徑”：

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

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

其他塊的剖分。

數(shù)據(jù)填寫和文獻(xiàn)保留:

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

注意：每輸入一層土的滲透系數(shù)，都要點擊邈按鈕進(jìn)行數(shù)據(jù)保留：

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

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

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

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

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

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

按鈕，注意：計算時文獻(xiàn)名必須為dqb.dat：

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

gjju

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

4.2運行計算

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

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

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

.DAT文獻(xiàn)的內(nèi)容和格式)；

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

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

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

輸入0并保留；

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

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

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

單擊OK,隨即出現(xiàn)土壩橫斷面圖；

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

擊OK,隨即出現(xiàn)滲流場圖形；

c、將土壩橫斷面圖與滲流場圖形進(jìn)行登加：用鼠標(biāo)將滲流場圖形拖至堤防橫斷面

圖，用鼠標(biāo)同步框選堤防橫斷面圖和單個滲流場圖形，單擊MAP出現(xiàn)MAP菜單，

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

d、進(jìn)行文字編輯、坐標(biāo)軸設(shè)定等，最終得到完美口勺滲流場分布圖。

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

4.3DQB.DAT文獻(xiàn)的內(nèi)容和格式

產(chǎn)生的數(shù)據(jù)文獻(xiàn)DQB.DAT內(nèi)容格式如下：

第一行：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）

數(shù)聽闡明：

（1）SEEPI—鑒別與否保留滲流繪圖數(shù)據(jù)，'Y'為肯定，'N'為否認(rèn)。

SEEP2—鑒別與否自動剖分?jǐn)?shù)據(jù)，'Y'為肯定，'N'為否認(rèn)。

SEEP3一設(shè)定為，N'。

（2）ZU一分塊組數(shù)（計算區(qū)域內(nèi)剖分的塊依次從左到右分為若干組，同一組間應(yīng)為上下關(guān)

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

KUAI一四邊形塊H勺數(shù)目：

DIAN一塊角點數(shù)：

BIAN一含各類已知邊界的邊數(shù);

TUCENG一不一樣滲透系數(shù)的土層種類數(shù)：

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

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

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

BFS一計算等勢線的百分?jǐn)?shù)詞距值，默認(rèn)為I0o

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

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

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

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

(4)ZB(I：DIAN,1：2)塊口勺結(jié)點坐標(biāo)數(shù)組

按塊結(jié)點序號次序表達(dá)每一結(jié)點破JX,Y坐標(biāo)。

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

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

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

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

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

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

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

3一—半有東流邊界

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

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

7一一上游已知水頭邊界

8一—下游已知水頭邊界

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

（7）SLKH（1：SLK）

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

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

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

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

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

4.4G.DAT數(shù)據(jù)內(nèi)容

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

G.DA1;G.DAT是本軟件進(jìn)行滲流計算的原始數(shù)據(jù)文獻(xiàn)，其數(shù)據(jù)內(nèi)容格式如卜.：

第一行：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一三角形單元總數(shù)

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）

各時段時間間隔天數(shù)

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

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

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

（17）QTE（I：NQ）

計算滲流量斷面日勺單元總數(shù)

(18)QEN(1：QE)

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

5算例

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

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

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

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

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

(2)單擊|視圖有關(guān)參矗“最大坐標(biāo)”中X填寫170,Y填寫100,“最小坐標(biāo)”中X填

寫-30,Y填寫10,其他的按照默認(rèn)值:

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

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

㈤匚Of居癡出曲百加八㈤

(4)保留dqb.dal文獻(xiàn)，內(nèi)容如下：

，、.，，丫，，“

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|按鈕,檢杳產(chǎn)生的DQB文獻(xiàn)與否對的（可以參照闡明竹中4.3DQB.DAT文

獻(xiàn)的內(nèi)容和格式）：

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

c、點擊|顯示剖分圖|按鈕,進(jìn)行剖分同步生成網(wǎng)格圖NET.BLN；

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

0并保留；

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

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

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

隨即出現(xiàn)上壩橫斷面圖：

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

隨即出現(xiàn)滲流場圖形：

將土壩橫斷面圖與滲流場圖形進(jìn)行疊加：用鼠標(biāo)將滲流場圖形拖至堤防橫斷面圖，同步

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

進(jìn)行文字編輯、坐標(biāo)軸設(shè)定等，最終得到完美的滲流場分布圖（如圖5.3）。

圖5.3均質(zhì)壩滲流場分布圖

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

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

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

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

圖5.4堤防橫斷面

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

圖5.5堤防橫斷面剖分圖

產(chǎn)生本問題日勺DQB-DAT文獻(xiàn)內(nèi)容如下:

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.