第一講 地下水?dāng)?shù)值模擬基礎(chǔ)知識(shí).ppt

1、第一講 地下水?dāng)?shù)值模擬基礎(chǔ)知識(shí),一、水循環(huán)與地下水,地下水賦存于地下介質(zhì),地下水是水循環(huán)的一個(gè)階段,土壤持水帶、包氣帶、飽和帶,The common assumption that the rate of groundwater withdrawal is “safe” or “sustainable” if it does not exceed the natural rate of recharge is not correct, because it ignores these changes in discharge from and recharge to the groundwat

2、er system,地下水的賦存,(1)巖石中的空隙與水 (2)包氣帶與飽水帶 (3)含水層、隔水層、弱透水層,包氣帶,特點(diǎn)： （1）巖石空隙未被水充滿； （2）固、液、氣三相介質(zhì)并存； （3）水的存在形式多樣：結(jié)合水、毛細(xì)水、重力水、氣態(tài)水。 包氣帶水的垂直分帶： （1）土壤水帶 （2）中間帶（過渡帶） （3）毛細(xì)水帶（支持毛細(xì)水帶、飽和毛細(xì)水帶）,飽水帶,巖石空隙被水完全充滿 是二相介質(zhì)（固相+液相水） 空隙中水的存在形式： 重力水 重力水：連續(xù)分布（孔隙是連邊）傳遞壓力在水頭差 作用下，地下水（空隙中的水）可以連續(xù)運(yùn)動(dòng)。 地下開挖，坑道，巷道，基坑，打井在此帶均有重力水涌出來！ 結(jié)合水,

3、含水層,隔水層與弱透水層：基本概念,飽水巖層中，根據(jù)巖層給水與透水能力而進(jìn)行的劃分： （1）含水層（Aquifer）：是能夠透過并給出相當(dāng)數(shù)量水的巖層，例如，各類砂土，砂巖等； （2）隔水層（Aquifuge)：不能透過與給出水或透過與給出的水量微不足道的巖層，例如，裂隙不發(fā)育的基巖、頁巖、板巖、粘土（致密）； （3）弱透水層（Aquitard）：滲透性很差，給出的水量微不足道，但在較大水力梯度作用下，具有一定的透水能力的巖層，例如，各種粘土，泥質(zhì)粉砂巖。,定義中的“相當(dāng)水量，微不足道，較大水力梯度”是模糊的；含水層與隔水層的劃分是相對(duì)的。從實(shí)際應(yīng)用來看，區(qū)分含水層與隔水層應(yīng)考慮巖層給出的水量

4、是否具有實(shí)際意義。從理論意義來看，巖層是否透水還取決于時(shí)間尺度。,含水層,隔水層與弱透水層：概念相對(duì)性,如華北平原早期地下水開采，深層水與淺層水之間水位差別不大，深層水與淺層水之間的粘土可作為隔水層；隨著深層水的不斷開采，水位大幅度降低，淺層水向深層“越流”，粘土層成為“透水層”。,含水層,隔水層與弱透水層：時(shí)間相對(duì)性,隔水層與弱透水層（諾曼與威瑟斯龐）,自然界中不存在絕對(duì)的隔水層，巖層是否透水取決于時(shí)間尺度。,含水層,隔水層與弱透水層: 應(yīng)用的相對(duì)性,巖性相同的地層根據(jù)不同研究目的可劃分為含水層或隔水層。 修水庫時(shí)，要考慮建庫后水庫是否滲漏？ 供水時(shí)，考慮水量是否足夠，是否為含水層？ 某組地

5、層是含水層還是隔水層？其界定要靈活運(yùn)用！ 現(xiàn)代水文地質(zhì)模擬計(jì)算，不再簡(jiǎn)單地劃分為含水層、隔水層，而是把不同巖層附于不同滲透參數(shù)。,地下水分類,廣義地下水（subsurface water）：地表以下巖石空隙中的水(包氣帶、飽水帶中的水) 狹義地下水（groundwater）：地表以下飽水帶巖石空隙中的水（重力水）,地下水分類：依據(jù)含水介質(zhì)類型、埋藏條件劃分含水介質(zhì)三類，埋藏條件分三類，組合共分為9類。,潛水:概念,潛水：地表以下第一個(gè)具有自由表面的穩(wěn)定含水層中的水 自由表面：即設(shè)有隔水層限制，與大氣直接相通，除大氣壓強(qiáng)外不受其它力 穩(wěn)定：具有一定的空間連續(xù)性（范圍）以示區(qū)分上層滯水 潛水含水層

6、：賦存潛水的巖層。 屋建筑時(shí)的基坑排水，大堤堤角處的散浸滲漏（潛水）,潛水:基本要素,潛水面 潛水位 潛水含水層 含水層厚度 潛水埋深,1 潛水含水層 2 隔水層 3、4 潛水面 潛水位 M 含水層厚度 D 潛水埋深 5 大氣降水入滲 6 蒸發(fā) 7 流向 8 泉,2,3,4,5,7,8,D,1,M,6,4,潛水:主要特征,補(bǔ)給：降水入滲，河湖入滲 排泄：泉，（河）泄流，蒸發(fā)補(bǔ)給或排泄通過含水層厚度變 化而儲(chǔ)水與釋水！ 動(dòng)態(tài)：受氣象，水文因素影響明顯， 變化快（水量、水位季節(jié)性變化） 受人為因素影響也顯著，易污染 水循環(huán)交替迅速：水循環(huán)周期短，更新恢復(fù)快,潛水等水位線圖,？一張潛水等水位線圖可以

7、提供哪些信息,承壓水:概念,1）定義： 充滿于2個(gè)隔水層（弱透水層）之間的含水層中的水。 2）基本要素：,承壓含水層 隔水頂板 隔水底板 承壓含水層厚 測(cè)壓水位線 承壓高度-H 補(bǔ)給區(qū) 承壓區(qū) 排泄區(qū) 自溢區(qū),基巖自流盆地中的承壓水,承壓水:主要特征,補(bǔ)給與排泄：有限區(qū)域與外界聯(lián)系，水循環(huán)交替慢，平均滯留時(shí)間長(zhǎng)（年齡老）可恢復(fù)性差。 水化學(xué)：變化較大，可以是淡水，也可能是鹵水。不易污染，一旦污染，很難凈化！ 動(dòng)態(tài)：較穩(wěn)定，如果分布面積大，厚度穩(wěn)定，則調(diào)節(jié)能力很強(qiáng)。,承壓水:含水層的儲(chǔ)水與釋水,彈性給水度Ss：承壓含水層中當(dāng)測(cè)壓水位下降1個(gè)單位，單位水平面積含水層柱體所釋放的水量。測(cè)區(qū)水位降低導(dǎo)

8、致： （1）含水層孔隙中水的壓力降低水體積膨脹釋水 （2）孔隙水壓力降低，含水層顆粒間有效應(yīng)力增加骨架被壓縮（顆粒不變骨架壓縮=空隙體積減小)發(fā)生釋水。 這兩部水是很有限的，所以e很?。ㄝ^d小10-210-3）,二、為什么需要開展地下水模擬工作,地下水相關(guān)問題,過度抽水、采礦疏干、地面沉降 地下水位上升、鹽漬化、工程問題 地下水污染 海水入侵與地下水排放,為了有效管理地下水，需要知道以下內(nèi)容,各含水層每年應(yīng)該開采的量 開采井位置、人工注射補(bǔ)給井位置，開采量 查明地下水質(zhì)量影響因素 地下水污染方面（1）工業(yè)廢棄物（2）垃圾填埋場(chǎng)（3）農(nóng)業(yè)活動(dòng)中化肥、殺蟲劑,地下水管理需要知道地下水系統(tǒng)對(duì)管理行為

9、的反應(yīng),管理就意味著做出的決定考慮一些特殊的制約因素 良好的管理需要知道政策執(zhí)行后管理對(duì)象的反應(yīng) 通過比較系統(tǒng)反應(yīng)，政策制定者就可以在各種策略之間比較 一旦知道含水層某處發(fā)現(xiàn)污染，地下水管理人員就要預(yù)測(cè)其路徑和歸宿，提出治理和修復(fù)計(jì)劃 地下水觀測(cè)網(wǎng)的布置，必須是在對(duì)地下水系統(tǒng)全面了解的基礎(chǔ)之上。,需要一個(gè)模型來提供地下水管理的信息,地下水運(yùn)動(dòng)的數(shù)學(xué)模型 化學(xué)組分遷移的數(shù)學(xué)模型 因?yàn)榉匠讨荒茉诤芤?guī)則、很簡(jiǎn)單的條件下得到解析解，所以用數(shù)值模擬的方法解方程，就是地下水?dāng)?shù)值模擬,數(shù)值模擬的類型,預(yù)測(cè)型模擬：預(yù)測(cè)人類活動(dòng)對(duì)地下水的影響，例如水位、水質(zhì)變化 研究型模擬：通過模擬認(rèn)識(shí)一些現(xiàn)象，例如通過考慮不

10、同化學(xué)作用條件下數(shù)值模擬，研究各種作用的影響大小 設(shè)計(jì)型模擬：設(shè)計(jì)監(jiān)測(cè)方案等 除了數(shù)值模擬，還有物理模擬模型，功能一樣,二、數(shù)學(xué)模型,描述地下水和化學(xué)組分的數(shù)學(xué)方程，加上初始條件、邊界條件，就是數(shù)學(xué)模型,32,微分方程,定解條件,邊界條件,初始條件,已知t=0時(shí)的因變量， H(x,y,z,0)=H0(x,y,z),已知水頭邊界(I類邊界) H(x,y,z,t)=f(x,y,z,t) (x,y,z)B1 特例：定水頭邊界 H(x,y,z,t)=C 已知流量邊界 特例：隔水邊界,數(shù)學(xué)模型,地下水運(yùn)動(dòng)的數(shù)學(xué)模型結(jié)構(gòu),1.地下水運(yùn)動(dòng)的基本方程,（1）連續(xù)性假設(shè) （2）達(dá)西定律 （3）水均衡原理,連續(xù)性

11、假設(shè),連續(xù)介質(zhì)力學(xué)的基本前提，假定整個(gè)物體的梯級(jí)都被組成物體的介質(zhì)充填，不留任何孔隙。這樣，物體內(nèi)的一些物理量，才可能是連續(xù)的，因而才可以用坐標(biāo)的連續(xù)函數(shù)表示他們的規(guī)律。 代表性單元體積：只要研究的對(duì)象體積足夠大，相鄰微粒之間的距離，都比物體尺寸小得多，連續(xù)性假設(shè)就不會(huì)產(chǎn)生顯著誤差,達(dá)西試驗(yàn),1856年，法國(guó)水力學(xué)家達(dá)西 H.Darcy 通過大量的室內(nèi)實(shí)驗(yàn)得 到線性滲流定律. 實(shí)驗(yàn)條件： 1）等徑圓筒裝入均勻砂樣，斷面為 ； 2）上（下各）置一個(gè)穩(wěn)定的溢水裝置穩(wěn)定水流 ； 3）實(shí)驗(yàn)時(shí)上端進(jìn)水，下端出水測(cè)出水量Q示意流線 ； 4）砂筒中安裝了2個(gè)測(cè)壓管。,試驗(yàn)裝置圖,達(dá)西定律,通過變水頭，多次實(shí)

12、驗(yàn)得出：出水端的流量Q與砂柱面積、測(cè)壓管水頭差 h 之間的關(guān)系為：,Q 滲透流量； 砂柱斷面面積； h 水頭損失；L 滲流途徑； K與試樣有關(guān)的比例常數(shù)。,水頭表達(dá)式,總水頭 測(cè)壓水頭,速度項(xiàng)很小,達(dá)西定律:滲透流速v,由于流速與流量的關(guān)系為：Q = v 與（2）式比較得: v = Kh/L=KI - （*） 式中 v稱為滲透流速。(*)式也稱為單位面積上的流量，即比流量。 達(dá)西定律的三個(gè)表達(dá)式中看出：滲透流速與水力梯度是一次方正比關(guān)系 ，故達(dá)西定律又稱為線性滲透定律。,達(dá)西定律,參數(shù)討論,1）滲透流速（v）和過水?dāng)嗝妫ǎ?Q = KI = V 過水?dāng)嗝媾c水流過斷面是否一致？,過水?dāng)嗝?是假想

13、的斷面 實(shí)際孔隙斷面 應(yīng)為n（孔隙度） 實(shí)際水流斷面 為ne (有效孔隙度),否！,達(dá)西定律,參數(shù)討論,滲透流速v=Q/ ; 實(shí)際流速u=Q/ 地下水滲透流速與平均實(shí)際流速的關(guān)系： v u ne 滲透流速v：是假設(shè)水流通過整個(gè)巖層斷面（骨架+空隙）時(shí)所具有的虛擬的平均流速。 意義：研究水量時(shí)，只考慮水流通過的總量與平均流速，而不去追蹤實(shí)際水質(zhì)點(diǎn)的運(yùn)移軌跡簡(jiǎn)化的研究,達(dá)西定律,參數(shù)討論:滲透系數(shù)K,有些教科書中也稱為水力傳導(dǎo)率（hydraulic conductivity） 定義：水力梯度為I1 時(shí)的滲透流速 （VKI） 滲透系數(shù)具有速度量綱L2T-1 由公式 VKI 分析： 當(dāng) I 一定時(shí)，巖

14、層的 K 愈大，則 V 也愈大， Q 也大。 因此，滲透系數(shù) K 是表征巖石透水性的定量指標(biāo)。,達(dá)西定律,達(dá)西定律,地下水水均衡原理 地下水的水量和鹽分在收支方面的數(shù)量關(guān)系，稱為地下水均衡。其中水量均衡為水均衡，鹽分均衡的話稱為鹽均衡。 在均衡期中，均衡區(qū)的補(bǔ)給量大于排泄量-正均衡。 在均衡期中，均衡區(qū)的補(bǔ)給量小于排泄量-負(fù)均衡。 地下水的均衡狀況是通過建立地下水均衡方程實(shí)現(xiàn)的。其原理就是水量平衡原理，一般: W=X + W1 + Z1 + Y1 + Z2 + W2 + Y2 其中：W均衡期內(nèi)地下水量的變化量，X：大氣降水的入滲補(bǔ)給量，W1：地下水流入量，Z1：凝結(jié)水補(bǔ)給量 ，Y1：地表水入滲

15、補(bǔ)給量，W2：地下水流出量，Z2：地下水蒸發(fā)量， Y2：地下水補(bǔ)給地表水量。,水均衡原理,44,取右圖所示得微小六面體。設(shè)與x, y, z,方向?qū)?yīng)得主滲透系數(shù)分別為Kx, Ky,Kz；建立均衡期t時(shí)段內(nèi)，微小均衡六面體的水量守恒方程。,45,同理，y、z-方向流入流出分別為:,x方向流入流出分別為:,t時(shí)段內(nèi)，六面體水量變化量為：,46,六面體內(nèi)地下水儲(chǔ)存量的變化為,由水均衡原理，得,方程兩端除以t，并取,和，則,47,一般密度的空間變化率很小，故,于是有,由達(dá)西定律,有,（2）水流連續(xù)性方程左端項(xiàng),都很小，可以忽略。,48,上式為非均質(zhì)各向異性承壓含水層的偏微分方程。,均質(zhì)各向異性非穩(wěn)定流

16、,均質(zhì)各向異性穩(wěn)定流,得到地下水三維流動(dòng)微分方程,1/L,53,溶質(zhì)運(yùn)移數(shù)學(xué)模型：緒論,隨著經(jīng)濟(jì)的快速發(fā)展，地下水被污染的程度日益嚴(yán)重，并引起了人們的廣泛關(guān)注，目前仍然存在很多問題題，迫切需要解決：,海（咸）水入侵：地下淡水的過量開采導(dǎo)致沿海地區(qū)的 地表污（廢）水排放和農(nóng)耕污染造成的硝酸鹽污染 石油和石油化工產(chǎn)品的污染 垃圾填埋場(chǎng)滲漏污染,54,早在1805年，F(xiàn)ick就提出了分子擴(kuò)散定律。1905年，Slichter報(bào)道了土壤中溶質(zhì)并不是以相同的速度運(yùn)移的現(xiàn)象。此后，人們逐漸提出并逐步形成了溶質(zhì)運(yùn)移的基本理論水動(dòng)力彌散理論。,溶質(zhì)運(yùn)移方程=Fick定律+質(zhì)量守恒原理,2、溶質(zhì)運(yùn)移數(shù)學(xué)模型,5

17、5,水動(dòng)力彌散理論,可混溶流體 兩種或兩種以上的流體在同一儲(chǔ)集空間中不存在明顯的突變界面，見下圖。如濱海含水層中海水入侵地下淡水。（示蹤劑） 不可混溶流體 多種（兩種或兩種以上）的流體在同一儲(chǔ)集空間中存在著明顯的突變界面，見下圖。如油、氣、水或其它有機(jī)物流體。（多相流體）,溶質(zhì)運(yùn)移數(shù)學(xué)模型,56,水動(dòng)力彌散理論（續(xù)1）,可混溶流體,石油污染物在含水層中運(yùn)移,不可混溶流體,不同性質(zhì)溶體之間無明顯的突變界,不同性質(zhì)溶體之間有明顯的突變界,油,水,污 染 物,水,溶質(zhì)運(yùn)移數(shù)學(xué)模型,57,水動(dòng)力彌散,分子擴(kuò)散,機(jī)械彌散,由濃度高的方向向濃度底 的方向運(yùn)動(dòng)，趨于均一,由于微觀多孔介質(zhì)中流速分布的不均一而

18、引起的示蹤劑（水質(zhì)點(diǎn)）濃度在地下水含水層中不均勻分布的現(xiàn)象。,1、水動(dòng)力彌散理論,兩部分,二、溶質(zhì)運(yùn)移數(shù)學(xué)模型,58,水動(dòng)力彌散理論：機(jī)械彌散原因,同一空隙中不同部位的流速分布不均勻 不同空隙的流速大小不同 固體骨架導(dǎo)致流速分布的不均勻,（1）,（2）,（3）,地下水質(zhì)點(diǎn)運(yùn)動(dòng)速度的差異是產(chǎn)生水動(dòng)力彌散的根本原因,二、溶質(zhì)運(yùn)移數(shù)學(xué)模型,59,平行于平均流速方向上的彌散,垂直于平均流速方向上的彌散,縱向彌散,橫向彌散,Fick定律,水動(dòng)力彌散理論：機(jī)械彌散原因,二、溶質(zhì)運(yùn)移數(shù)學(xué)模型,描述溶質(zhì)彌散規(guī)律,Fick定律,60,控制方程：質(zhì)量守恒定理,在多組分組成的流體體系中任取一點(diǎn)P (x, y, z)

19、,以P為中心取一微小的質(zhì)量平衡體，其側(cè)面分別平行與3個(gè)坐標(biāo)面，邊長(zhǎng)分別為x、y、 z。,質(zhì)量守恒原理:在時(shí)間t內(nèi)，組分在這個(gè)單元體中的凈流出（或流出）量（暫不考慮起內(nèi)部有質(zhì)量產(chǎn)生和消失），應(yīng)等于這個(gè)單元中組分的質(zhì)量變化,二、溶質(zhì)運(yùn)移數(shù)學(xué)模型,化學(xué)組分?jǐn)?shù)學(xué)模型,任意時(shí)間段內(nèi)化學(xué)組分儲(chǔ)存量的變化（包括溶解相和吸附相）等于 彌散、對(duì)流、源匯和化學(xué)反應(yīng)造成變化的總和,三、數(shù)值模擬方法,數(shù)值方法很多，但是最簡(jiǎn)單實(shí)用的是有限差分法：,有限差分法 有限單元法 積分有限差分法 半解析半數(shù)值法 邊界元法 有限體積法,Numerical Methods,All numerical methods involve

20、representing the flow domain by a limited number of discrete points called nodes. A set of equations are then derived to relate the nodal values of the dependent variable such that they satisfy the governing PDE, either approximately or exactly. Numerical Solutions Discrete solution of head at selec

21、ted nodal points. Involves numerical solution of a set of algebraic equations.,有限差分法的基本原理,有兩種方法建立差分方程 方法一以地下水流基本微分方程及其定解條件為基礎(chǔ)， 在滲流區(qū)剖分基礎(chǔ)上，用差商代替微商，將地下水流微分方程的求解轉(zhuǎn)化為差分方程（代數(shù)方程）求解。適用于二維矩形網(wǎng)格剖分、三維長(zhǎng)方體網(wǎng)格剖分。 方法二在滲流區(qū)剖分的基礎(chǔ)上，直接由達(dá)西定律和水均衡原理，建立各個(gè)均衡區(qū)的水均衡方程，從而得到差分方程。適用于矩形網(wǎng)格、三角形網(wǎng)格。,矩形網(wǎng)格,多邊形網(wǎng)格,1、網(wǎng)格劃分的基本類型,（1）先劃格線，格點(diǎn)位于網(wǎng)格中

22、心,均衡網(wǎng)格,節(jié)點(diǎn)網(wǎng)格,（2）先規(guī)定格點(diǎn)位置，再垂直平分兩相鄰結(jié)點(diǎn)的連線作格線，形成的網(wǎng)格即為水均衡區(qū),方法一：差商代替微商,MODFLOW網(wǎng)格系統(tǒng),方法一：差商代替微商,四、地下水?dāng)?shù)值模型建模步驟,模擬步驟,建立概念模型 建立數(shù)學(xué)模型 數(shù)值方法及軟件（編程） 參數(shù)準(zhǔn)備以及計(jì)算 模型校正與檢驗(yàn) 參數(shù)敏感性分析 預(yù)測(cè)模擬,軟件,一、概念模型（模型概化）,根據(jù)詳細(xì)的地形地貌、地質(zhì)、水文地質(zhì)、構(gòu)造地質(zhì)、水文地球化學(xué)、巖石礦物、水文、氣象、工農(nóng)業(yè)利用情況等 模擬的區(qū)域: 含水層類型: 潛水（無壓）、承壓、混合、多層 維數(shù):一維、二維、三維 水流狀態(tài)：穩(wěn)定流/非穩(wěn)定流、飽和流/非飽和流 介質(zhì)狀況: 均質(zhì)

23、和非均質(zhì)/各向同性和各向異性 孔隙/裂隙/雙重介質(zhì) 流體的密度差 邊界條件和初始條件 必要時(shí)需進(jìn)行一系列的室內(nèi)試驗(yàn)與野外試驗(yàn), 以獲取有關(guān)參數(shù), 如滲透系數(shù)、彌散系數(shù)、分配系數(shù)、反應(yīng)速率常數(shù)等。,（二）數(shù)學(xué)模型,三維地下水流動(dòng)問題控制方程,第二類邊界條件,第一類邊界條件,初始條件,絕大部分?jǐn)?shù)學(xué)模型是無法用解析法求解的，數(shù)值化就是將數(shù)學(xué)模型轉(zhuǎn)化為可解的數(shù)值模型。,三、數(shù)值方法及軟件（或編程）,有限差分法 有限單元法 積分有限差分法 半解析半數(shù)值法 邊界元法,有限差分法: MOFLOW系列 GMS 中MODFLOW Visual MODFLOW Processing MODFLOW 有限單元法:

24、FEFLOW 積分有限差分法: TOUGH2，TOUGH REACT,三、數(shù)值方法及軟件（或編程）,四、模型參數(shù)準(zhǔn)備以及計(jì)算,含水層參數(shù)：滲透系數(shù)，彈性釋水系數(shù)，孔隙度等 源匯項(xiàng)： 大氣降水入滲系數(shù)（分區(qū)、數(shù)值） 蒸發(fā)排泄系數(shù) 地表水體水位、底面高程、底面巖性特征 渠系灌溉入滲系數(shù) 人工開采（點(diǎn)狀、面狀） 邊界條件 初始條件,計(jì)算,資料,結(jié)果,輸入模塊,運(yùn)行模塊,輸出模塊,Visual MODFLOW,參數(shù)不確定性,鉆孔太少，地層資料少，鉆孔多，含水層結(jié)構(gòu)會(huì)發(fā)生變化,四、模型校正與檢驗(yàn),（1）鉆孔資料（2）抽水試驗(yàn),參數(shù)不確定性,泰斯井流抽水試驗(yàn)假設(shè)： （1） （2） （3） （4） （5）

25、（6）,（四）模型校正與檢驗(yàn),（1）鉆孔資料（2）抽水試驗(yàn),模擬：模型校正（參數(shù)識(shí)別）,將模擬結(jié)果與實(shí)測(cè)結(jié)果比較,進(jìn)行參數(shù)調(diào)整, 使模擬結(jié)果在給定的誤差范圍內(nèi)與實(shí)測(cè)結(jié)果吻合。 調(diào)參過程是一個(gè)復(fù)雜而辛苦的工作, 所調(diào)整的參數(shù)必須符合模擬區(qū)的具體情況。 人機(jī)交互與自動(dòng)調(diào)參相結(jié)合。盡管自動(dòng)調(diào)參程序(如PEST ) , 也不能代替人的工作。,模擬：模型檢驗(yàn),模型驗(yàn)證是在模型校正的基礎(chǔ)上, 進(jìn)一步調(diào)整參數(shù), 使模擬結(jié)果與第二次實(shí)測(cè)結(jié)果吻合, 以進(jìn)一步提高模型的置信度。,靈敏度分析,校正后的模型受參數(shù)值的時(shí)空分布、邊界條件、水流狀態(tài)等不確定度的影響。靈敏度分析就是為了確定不確定度對(duì)校正模型的影響程度。,預(yù)測(cè),用校正的參數(shù)值進(jìn)行預(yù)測(cè), 預(yù)測(cè)時(shí)需估算未來的水流狀態(tài)。,后續(xù)檢查與模型的重建（完善）,后續(xù)檢查在模擬研究結(jié)束數(shù)年后進(jìn)行。收集新的野外數(shù)據(jù)以確定預(yù)測(cè)結(jié)果是否