“黃河流域水資源演變規(guī)律與二元演化模型” 研究成果簡介[J]

1、把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 “黃河流域水資源演變規(guī)律與二元演化模型” 研究成果簡介 賈仰文 王 浩 摘要摘要: 本文簡要介紹了在國家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃(“973”)項(xiàng)目“黃河流域水資源演變規(guī)律 與可再生性維持機(jī)理”第二課題(G1999043602)中取得的三項(xiàng)主要成果：1) 流域水資源 二元演化模型, 2) 黃河流域水資源全口徑層次化動態(tài)評價(jià), 3) 人類活動影響下的黃河 流域水資源演化規(guī)律。首先，將分布式流域水文模型 (WEP-L) 和集總式水資源調(diào)配 模型(WARM)相耦合，建立了流域水資源二

2、元演化模型。然后，在界定“廣義水資源” 與“狹義水資源”概念的基礎(chǔ)上，提出了水資源全口徑層次化動態(tài)評價(jià)方法，并給出了 黃河流域評價(jià)成果。最后，通過比較 2000 年現(xiàn)狀下墊面條件下與歷史系列條件下的 評價(jià)結(jié)果，以及考慮與不考慮人工取用水條件下的模擬結(jié)果，初步發(fā)現(xiàn)：1）黃河流 域在強(qiáng)烈的人類活動影響下，水資源量及其構(gòu)成均發(fā)生了顯著變化，地表水資源量衰 減,而不重復(fù)地下水量增加；2）在狹義水資源衰減的同時(shí)，伴隨著有效蒸散即降水有 效利用的增加，流域廣義水資源量有一定幅度增加。研究成果對客觀評價(jià)水土保持生 態(tài)建設(shè)、農(nóng)田基本建設(shè)等土地利用變化的水文水資源效應(yīng)，對黃河的治理規(guī)劃和水安 全戰(zhàn)略具有重要參考

3、應(yīng)用價(jià)值。 1 前言 “黃河流域水資源演變規(guī)律與二元演化模型”研究是國家重點(diǎn)基礎(chǔ)研究發(fā)展 規(guī)劃（“973”）項(xiàng)目“黃河流域水資源演變規(guī)律與可再生性維持機(jī)理”中的第二課 題(G1999043602)，研究歷時(shí) 5 年（1999 年 10 月2004 年 9 月） 。該課題主要 完成單位是：中國水利水電科學(xué)研究院水資源研究所、中國科學(xué)院地理科學(xué)與 資源研究所、水利部黃河委員會水文局、中國科學(xué)院地質(zhì)與地球物理研究所， 主要完成人員有：王 浩、賈仰文、王建華、秦大庸、李麗娟、羅翔宇、周祖昊、 嚴(yán)登華、王 玲、張學(xué)成、劉廣全、秦大軍、張新海、江 東、楊貴羽。受篇幅 限制，本文將簡要介紹在該課題中取得的三

4、項(xiàng)主要成果，即：1)流域水資源二 元演化模型, 2) 黃河流域水資源全口徑層次化動態(tài)評價(jià), 3) 人類活動影響下的 把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 黃河流域水資源演化規(guī)律。詳細(xì)內(nèi)容請參見有關(guān)參考文獻(xiàn)1,2。 隨著人類活動的日益加劇，天然狀態(tài)下的流域水循環(huán)模式發(fā)生了根本性改 變，出現(xiàn)了由“取水輸水用水排水回歸”五個(gè)基本環(huán)節(jié)構(gòu)成的人工側(cè)枝 循環(huán)圈，形成“天然人工”雙驅(qū)動力作用下的流域“二元”水循環(huán)模式。這種“二 元”模式導(dǎo)致流域水資源形勢加速演變，洪澇與旱災(zāi)發(fā)生的幾率與不確定性增加， 同樣規(guī)模的洪水造成的損

5、失急劇增大，同時(shí)還引起生態(tài)與環(huán)境變化。 迄今為止，國內(nèi)外水資源評價(jià)方法與實(shí)踐均是基于“實(shí)測還原”的一元靜 態(tài)模式3 6，即通過實(shí)測水文要素后，再把實(shí)測水文系列中隱含的人類活動影 響扣除， “還原”到流域水資源的天然“本底”狀態(tài)。但是，隨著人類活動日益加劇, 還原比例越來越大，受資料條件等客觀因素的限制和選取還原參數(shù)時(shí)人為主觀 隨意性的影響，應(yīng)用還原法難以獲取“天然”和“人工”二元驅(qū)動力作用下的水資 源量“真值”。同時(shí)，現(xiàn)行水資源評價(jià)方法還存在以下問題：1）以地表水和地下 水構(gòu)成的徑流性水資源為評價(jià)口徑，評價(jià)口徑不全面，難以反映水資源的多元 有效性（如生態(tài)植被對土壤水的有效利用等） ；2）以分離

6、評價(jià)為基本模式，如 地表水評價(jià)與地下水評價(jià)相分離、水資源量評價(jià)與開發(fā)利用評價(jià)相分離，難以 適應(yīng)水資源綜合規(guī)劃需求；3）難于反映下墊面變化對現(xiàn)狀可利用水資源的影響； 4）采取分區(qū)集總式的評價(jià)方法，在描述水資源的空間變異特征和指導(dǎo)開發(fā)利用 方面存在一定的局限。因此傳統(tǒng)水資源評價(jià)方法和評價(jià)手段亟待改進(jìn)。 黃河流域(面積 795 000 km2) 在中國社會經(jīng)濟(jì)發(fā)展中一直發(fā)揮著重要作用。 黃河流域概況見圖 1。九十年代以來，黃河流域出現(xiàn)了水資源嚴(yán)重短缺，生態(tài) 環(huán)境惡化, 下游斷流持續(xù)時(shí)間越來越長。這種狀況既與近年降水偏少有關(guān)，又 與強(qiáng)烈的人類活動影響密不可分。據(jù)統(tǒng)計(jì)，目前黃河流域干流共建成 12 座大

7、 型水庫 (總庫容 563 億 m3)、 支流大中型水庫 (總庫容 100 億 m3) 以及大量 水土保持設(shè)施；灌溉面積由二十世紀(jì)五十年代的 80 萬 ha 增加到 2000 年 730 萬 ha。人類活動不僅改變了水循環(huán)過程，同時(shí)對水資源的數(shù)量及其構(gòu)成產(chǎn)生重 大影響。因此，開展黃河水資源動態(tài)評價(jià)，即分析黃河流域水資源在過去、現(xiàn) 在及未來下墊面條件下的變化并找出其演化規(guī)律，對黃河流域水資源可持續(xù)開 發(fā)利用具有重要意義。 把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 圖圖 1 黃河流域概況黃河流域概況 2 模型開發(fā) 2

8、0 世紀(jì) 80 年代中期以來，隨著計(jì)算機(jī)技術(shù)、地理信息系統(tǒng)和遙感技術(shù)的 進(jìn)步，分布式流域水文模型7,8得到很大發(fā)展?；谖锢頇C(jī)制的分布式流域水文 模型，可結(jié)合遙感、航測、數(shù)字流域9及地理信息系統(tǒng)等新技術(shù)進(jìn)行較為細(xì)致 的水文分析與預(yù)測，能夠分析不同下墊面條件下的流域水資源演變情景，為開 展水資源動態(tài)評價(jià)創(chuàng)造了條件。 本研究將大流域分布式水文模型 (WEP-L) 和集總式水資源調(diào)配模型 (WARM)相耦合，建立起流域水資源二元演化模型。WEP-L (Water and Energy transfer Processes in Large river basins) 模型是在綜合了分布式水文模型和陸

9、面過 程模型各自優(yōu)點(diǎn)的基礎(chǔ)上開發(fā)的，模擬對象為“天然人工”二元水循環(huán)系統(tǒng)。 該模型建立在 WEP 模型基礎(chǔ)之上，但為適應(yīng)黃河這樣的超大流域，采用“子流 域內(nèi)等高帶”為計(jì)算單元，并用“馬賽克”法考慮計(jì)算單元內(nèi)土地覆被的多樣性， 避免了采用過粗網(wǎng)格單元產(chǎn)生的模擬失真問題。針對各水循環(huán)要素過程時(shí)間尺 度不同的特點(diǎn)，計(jì)算時(shí)采用了 1h 至 1d 的“變時(shí)間步長”，既合理表述了水循環(huán) 動力學(xué)機(jī)制又提高了計(jì)算效率。研究中根據(jù) 1km 分辨率 DEM(數(shù)字高程模型)以 及數(shù)字化河網(wǎng)等空間信息數(shù)據(jù),將全黃河流域劃分為具有空間拓?fù)潢P(guān)系的 8485 把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。

10、 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 個(gè)子流域和 38720 個(gè)等高帶，采用 45 年（19562000）水文氣象系列數(shù)據(jù)及相 應(yīng)下墊面條件進(jìn)行了模擬計(jì)算，并根據(jù)黃河流域主要水文站逐月和逐日徑流系 列進(jìn)行了模型校驗(yàn)。 2.1 大流域分布式水文模型 WEP-L WEP-L 模型各基本計(jì)算單元內(nèi)的垂直方向結(jié)構(gòu)如圖 2(a)所示。從上到下包 括植被或建筑物截留層、地表洼地儲留層、土壤表層、過渡帶層、淺層地下水 層和深層地下水層等。狀態(tài)變量包括植被截留量、洼地儲留量、土壤含水率、 地表溫度、過渡帶層儲水量、地下水位及河道水位等。主要參數(shù)包括植被最大 截留深、土壤滲透系數(shù)、土壤水分

11、吸力特征曲線參數(shù)、地下水透水系數(shù)和產(chǎn)水 系數(shù)、河床的透水系數(shù)和坡面、河道的糙率等。為考慮基本計(jì)算單元內(nèi)土地利 用的不均勻性，采用了“馬賽克”法即把基本計(jì)算單元內(nèi)的土地歸成數(shù)類，分別 計(jì)算各類土地類型的地表面水熱通量，取其面積平均值為計(jì)算單元的地表面水 熱通量。土地利用首先分為裸地植被域、灌溉農(nóng)田、非灌溉農(nóng)田、水域和不 透水域 5 大類。裸地植被域又分為裸地、草地和林地類，不透水域分為城 市地面與都市建筑物 2 類。另外，為反映表層土壤的含水率隨深度的變化和便 于描述土壤蒸發(fā)、草或作物根系吸水和樹木根系吸水，將透水區(qū)域的表層土壤 分割成 3 層。在水循環(huán)各要素模擬中，截留、土壤、水面和植被蒸騰等

12、蒸發(fā)項(xiàng)， 按照土壤-植被-大氣通量交換方法(SVATS)、采用 NoilhanPlanton 模型、 Penman 公式和 Penman-Monteith 公式等進(jìn)行了詳細(xì)計(jì)算。地表徑流分為超滲和 蓄滿兩種產(chǎn)流模式分別采用 GreenAmpt 模型和 Richards 方程計(jì)算，對山坡斜 面土壤層進(jìn)行了壤中流計(jì)算，淺層地下水運(yùn)動進(jìn)行了二維數(shù)值計(jì)算并與非飽和 土壤水及河水進(jìn)行了動態(tài)耦合，積雪融化過程采用溫度指標(biāo)法。詳細(xì)內(nèi)容請參 見 WEP 模型的有關(guān)參考文獻(xiàn)1,10,11。 把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。

13、 (a) (b) 河道匯流 坡地匯流 基本計(jì)算單元 （等高帶） 1 2 3 9 4 5 q2 q3 q4 q6 6 8 q57 Q1 Q2 Q3 Q4 Q5 Q6 Q7 Q8 Q9 q1 圖圖 2大流域分布式水文模型大流域分布式水文模型 WEP-L 的結(jié)構(gòu)：的結(jié)構(gòu)：(a) 垂向結(jié)構(gòu)（基本計(jì)算單元內(nèi)）垂向結(jié)構(gòu)（基本計(jì)算單元內(nèi)） ，(b)平面結(jié)構(gòu)平面結(jié)構(gòu) WEP-L 模型的平面結(jié)構(gòu)如圖 2(b)所示。坡面匯流計(jì)算根據(jù)各等高帶的高程、 坡度與 Manning 糙率系數(shù)，采用維運(yùn)動波法將坡面徑流由流域的最上游端追 跡計(jì)算至最下游端。各條河道的匯流計(jì)算，根據(jù)有無下游邊界條件采用維運(yùn) 動波法或動力波法計(jì)算

14、。地下水流動分山丘區(qū)和平原區(qū)分別進(jìn)行數(shù)值解析，并 考慮其與地表水、土壤水及河道水的水量交換。 2.2 集總式水資源調(diào)配模型 WARM 取用水等人工側(cè)支水循環(huán)過程與自然水循環(huán)過程需要耦合模擬，以反映流 域水資源實(shí)際演化機(jī)理。為此，本研究開發(fā)了集總式水資源調(diào)配模型 WARM (Water Allocation and Regulation Model), WARM 模型包括兩部分: 一是水資源 合理配置模型（結(jié)構(gòu)示意見圖 3） ，二是水資源調(diào)度模型。 把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 WEP-L 模型與 WA

15、RM 模型的耦合見圖 4: 1) 具體通過 WEP-L 分布式模 擬結(jié)果的時(shí)空尺度聚合與 WARM 集總式調(diào)控結(jié)果的時(shí)空展布來實(shí)現(xiàn); 2) 耦合模 擬包括兩種類型，一是歷史系列或是某個(gè)時(shí)間斷面的模擬；二是對于未來過程 或某一時(shí)間斷面的模擬; 3) 耦合是雙向和交互式的：WEP-L 模型為 WARM 模 型提供來水信息、WARM 模型為 WEP-L 模型提供各類用水要求和水利工程的 調(diào)度規(guī)則，并通過信息反饋不斷修正。 圖圖 3 水資源合理配置模型結(jié)構(gòu)示意圖水資源合理配置模型結(jié)構(gòu)示意圖 反饋 輸入 WEP 模型 月調(diào)度模型 實(shí)時(shí)調(diào)度模型 控制 反饋 合理配置模型 年調(diào)度模型 控制 反饋 控制 反饋

16、 水資源調(diào)配模型 圖圖 4 大流域分布式水文模型大流域分布式水文模型 WEP-L 與水資源調(diào)配模型與水資源調(diào)配模型 WARM 的耦合的耦合 2.3 輸入數(shù)據(jù)處理與模型參數(shù) 輸入數(shù)據(jù)的處理包括：（1）河網(wǎng)水系生成、子流域劃分及其編碼，基本計(jì) 算單元（等高帶）的劃定與計(jì)算順序的確定；（2）降水等氣象要素的時(shí)空展布； （3）下墊面要素信息（如土地利用、土壤、水文地質(zhì)、植被、水庫、湖泊、河 道、灌區(qū)、水土保持等）的綜合處理；（4）社經(jīng)濟(jì)要素（人口、GDP、灌溉面 積、糧食產(chǎn)量等）和各類取用水信息的時(shí)空展布；（5）準(zhǔn)備 8485 個(gè)子流域、 河道（包括空間拓?fù)潢P(guān)系）的基本屬性表，38720 個(gè)計(jì)算單元的

17、基本屬性表， 以及上述各類數(shù)據(jù)的輸入文件，為 WEP-L 模型的應(yīng)用創(chuàng)造條件。 WARM 把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 WEP-L 模型的參數(shù)可分為三類。一類是地表面及河道系統(tǒng)參數(shù)，包括坡面 和河道的曼寧糙率、河床材料的覆蓋厚度及透水系數(shù)、城市土地利用的不透水 率及地表洼地最大儲留深；再一類是植被參數(shù)，包括植被覆蓋率、最大截留深、 葉面指數(shù)（LAI） 、空氣動力學(xué)阻抗、葉孔阻抗及根系分布參數(shù)；最后一類是土 壤與含水層參數(shù)，包括土壤層厚度、土壤空隙度、土壤入滲濕潤峰吸力、飽和 土壤導(dǎo)水系數(shù)、土壤水分吸力

18、特征曲線參數(shù)、土壤水分導(dǎo)水系數(shù)關(guān)系參數(shù)、 含水層厚度、含水層導(dǎo)水系數(shù)和含水層產(chǎn)水系數(shù)等。所有參數(shù)均有物理意義， 理論上講可根據(jù)觀測實(shí)驗(yàn)數(shù)據(jù)和遙感數(shù)據(jù)進(jìn)行推算而無需率定。但由于這些參 數(shù)在每個(gè)計(jì)算單元內(nèi)仍具有空間變異性，模擬計(jì)算時(shí)往往使用其單元內(nèi)平均參 數(shù)或稱有效參數(shù)，因此通常仍根據(jù)流量過程線及地下水位的觀測結(jié)果對一些關(guān) 鍵參數(shù)進(jìn)行適當(dāng)調(diào)整。WEP-L 模型的關(guān)鍵參數(shù)包括飽和土壤導(dǎo)水系數(shù)、城市土 地利用的不透水率、地表洼地最大儲留深、曼寧糙率、含水層導(dǎo)水系數(shù)、含水 層產(chǎn)水系數(shù)和河床材料透水系數(shù) 7 個(gè)。 2.4 模型驗(yàn)證 本研究進(jìn)行了 1956 至 2000 年共 45 年的變時(shí)間步長(1h-1

19、d)連續(xù)模擬計(jì)算。 其中 1980 年至 2000 年的 21 年取為模型校正期，主要校正參數(shù)包括土壤飽和導(dǎo) 水系數(shù)、河床材料透水系數(shù)和 Manning 糙率、各類土地利用的洼地最大截留深 以及地下水含水層的傳導(dǎo)系數(shù)及給水度等。校正準(zhǔn)則包括：（1）模擬期年均徑 流量誤差盡可能?。唬?）Nash-Sutcliffe 效率盡可能大；（3）模擬流量與觀測 流量的相關(guān)系數(shù)盡可能大。模型校正后，保持所有模型參數(shù)不變，對 1956 年至 2000 年（驗(yàn)證期）的連續(xù)模擬結(jié)果進(jìn)行驗(yàn)證，驗(yàn)證基礎(chǔ)為 23 個(gè)主要水文測站 45 年逐日和逐月實(shí)測徑流系列。 黃河部分主要水文斷面實(shí)測河川徑流模擬校驗(yàn)結(jié)果見表 1。唐

20、乃亥站、蘭 州站和花園口站 45 年天然河川月徑流過程模擬校驗(yàn)結(jié)果見圖 5。從校驗(yàn)結(jié)果來 看，19562000 系列各站多年平均實(shí)測徑流量、模型 Nash-Sutcliffe 效率系數(shù)以 及月徑流過程均擬合的比較好。因此，所構(gòu)建的模型具有較高模擬精度，可應(yīng) 用于黃河流域二元水循環(huán)過程模擬和水資源演變規(guī)律分析。 表表 1 黃河部分主要水文斷面實(shí)際河川徑流模擬校驗(yàn)結(jié)果黃河部分主要水文斷面實(shí)際河川徑流模擬校驗(yàn)結(jié)果 水文站實(shí)測年均徑流（108m3）計(jì)算年均徑流（108m3）相對誤差 Nash-Sutcliffe 效 率 把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己

21、說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 唐乃亥 (黃河干流）205.6203.1-1.2%0.821 貴德 (黃河干流）207.2210.41.5%0.819 蘭州 (黃河干流）313.1312.9-0.1%0.787 三門峽 (黃河干流）357.9353.0-1.4%0.621 花園口 (黃河干流）390.6395.31.2%0.689 (a) 0 1000 2000 3000 4000 5000 6000 1956 1958 1960 1962 1964 1966 1968 1970 1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 19

22、92 1994 1996 1998 2000 月流量 (m3/s) 0 100 200 300 400 500 月降水量 (mm) 面積平均降水量模擬流量實(shí)測流量 (b) 0 1000 2000 3000 4000 5000 6000 1956 1958 1960 1962 1964 1966 1968 1970 1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 月流量(m3/s) 0 100 200 300 400 500 月降水量 (mm) 面積平均降水量模擬流量實(shí)測流量 (c) 0 200

23、0 4000 6000 8000 10000 12000 1956 1958 1960 1962 1964 1966 1968 1970 1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 月流量 (m3/s) 0 100 200 300 400 500 月降水量(mm) 面積平均降水量模擬流量實(shí)測流量 圖圖 5 實(shí)測月徑流過程模擬校驗(yàn)結(jié)果：實(shí)測月徑流過程模擬校驗(yàn)結(jié)果：(a) 唐乃亥站，唐乃亥站，(b)蘭州站，蘭州站，(c) 花園口站花園口站 3 水資源層次化動態(tài)評價(jià)方法 水資源的本質(zhì)特征包括有效

24、性、可控性和可再生性三個(gè)方面，水資源評價(jià) 相應(yīng)地有三個(gè)評價(jià)準(zhǔn)則。有效性是指，對社會經(jīng)濟(jì)發(fā)展和生態(tài)環(huán)境具有效用的 水分均可以看作是水資源；可控性是指，在對人類具有效用的水分中，通過采 把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 取工程措施可以開發(fā)利用的那一部分水分；可再生性是指，水資源在流域水循 環(huán)形成和轉(zhuǎn)化過程中，其作為可再生性資源的充分必要條件是保持流域水循環(huán) 過程的相對穩(wěn)定。由此可見，根據(jù)有效性準(zhǔn)則定義的水資源包括了降水中的有 效部分和徑流性水資源，是一種廣義水資源；在廣義水資源范疇內(nèi)，根據(jù)可控 性準(zhǔn)則定義的水

25、資源是傳統(tǒng)水資源評價(jià)的徑流性水資源，即狹義水資源。在狹 義水資源范疇內(nèi)，進(jìn)一步從可再生性準(zhǔn)則出發(fā)，在自然條件和經(jīng)濟(jì)條件允許的 情況下能夠被工程系統(tǒng)一次性開發(fā)利用的最大潛在量可定義為水資源的國民經(jīng) 濟(jì)可利用量。 廣義水資源總量等于不重復(fù)的有效蒸散發(fā)量與狹義水資源總量之和，可用 式（1）表示： Ws（RsRg）EpEssEes （1） 上式中： Ws廣義水資源量；Rs地表水資源量；Rg不重復(fù)的地下水 資源量；Ep冠層截流蒸發(fā)量；Ess地面截流有效蒸發(fā)；Ees與地表水、地 下水不重復(fù)的土壤水有效蒸發(fā)量。 本研究認(rèn)為，林草冠層蒸發(fā)可直接降低植物表面和體內(nèi)的溫度，對維護(hù)植 物正常生理是有益的，因此林草冠

26、層截流蒸發(fā)是有效的；居民與工業(yè)用地是人 類居住和活動的集散地，蒸發(fā)可以起到降溫、濕潤等直接的環(huán)境作用，本次研 究認(rèn)為這一部分蒸發(fā)也是有效的。對于地面截流蒸發(fā)，本次研究認(rèn)為居民與工 業(yè)用地上的地面截流蒸發(fā)、作物和林草棵間截流蒸發(fā)分別對于人類和生態(tài)環(huán)境 主體是有直接環(huán)境效用的，應(yīng)將其納入到有效的廣義水資源范疇。對于難利用 土地截流蒸發(fā)（沼澤地除外）、稀疏草地的大棵間截流蒸發(fā)等都作為無效蒸發(fā)； 對于土壤水蒸散發(fā)當(dāng)中，蒸騰耗散的水分直接參與了生物量的生成屬于有效水 分，另外居民與工業(yè)用地土壤蒸發(fā)、作物和林草棵間土壤蒸發(fā)對于人類和生態(tài) 環(huán)境主體也有直接環(huán)境效用，作為有效蒸發(fā)。潛水蒸散發(fā)是從狹義的地下水資

27、 源轉(zhuǎn)化而來，全部屬于有效水分。 本研究中狹義水資源評價(jià)的口徑與傳統(tǒng)水資源評價(jià)口徑一致，即地表水資 源量與不重復(fù)的地下水資源量之和。 4 水資源評價(jià)結(jié)果 把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 在 2000 年現(xiàn)狀下墊面和取用水影響情景下，黃河流域 19562000 年系列 各二級分區(qū)的廣義水資源評價(jià)結(jié)果見表 2。從表 2 結(jié)果看出，黃河流域 19562000 系列平均降水量為 3563108m3，現(xiàn)狀下墊面和取用水影響情景下平 均廣義水資源量為 2756.6108m3占降水的 77.4，其中非徑流性有效水分為

28、 2080108m3，是徑流性狹義水資源的 3.1 倍，這表明以土壤水等其他形式賦存 的有效降水對于經(jīng)濟(jì)和生態(tài)系統(tǒng)起到重要的支撐作用。在各二級區(qū)當(dāng)中，三 花區(qū)間的廣義水資源占降水比例最高，占 94.7，而內(nèi)流區(qū)比例廣義水資源比 例最低，占 54.2。 表表 2 現(xiàn)狀條件下黃河流域廣義水資源評價(jià)結(jié)果現(xiàn)狀條件下黃河流域廣義水資源評價(jià)結(jié)果 （單位：（單位：108m3） 有效降水利用量 水資源分區(qū)降水量 狹義 水資源 農(nóng)田 有效蒸散 林草有 效蒸散 居工地 有效蒸發(fā) 總量 廣義 水資源 黃黃 河河 區(qū)區(qū)3563.0 676.4890.91173.315.92080.12756.6 龍羊峽以上632.3

29、212.15.5227.20.1232.8444.8 龍羊峽至蘭州433.0116.148.5182.81.0232.3348.4 蘭州至河口鎮(zhèn)427.653.7119.6116.72.8239.1292.9 河口鎮(zhèn)至龍門480.249.2144.6143.00.6288.2337.4 龍門至三門峽1038.9143.5386.8330.16.7723.6867.1 三門峽至花園口274.750.389.6118.51.6209.7260.0 花園口以下157.832.085.820.82.9109.5141.5 內(nèi) 流 區(qū)118.619.510.534.10.144.764.3 表表 3 現(xiàn)

30、狀條件下黃河流域狹義水資源評價(jià)現(xiàn)狀條件下黃河流域狹義水資源評價(jià) （單位：（單位：108m3） 地下水資源量 水資源分區(qū)地表水資源量 資源總量不重復(fù)資源量 狹義水資源總量 黃黃 河河 區(qū)區(qū)548.7404.2127.7676.4 龍羊峽以上210.165.31.9212.1 龍羊峽至蘭州112.837.03.4116.1 蘭州至河口鎮(zhèn)18.558.635.253.7 河口鎮(zhèn)至龍門42.340.06.949.2 龍門至三門峽104.5125.139.0143.5 三門峽至花園口39.235.111.050.3 花園口以下18.023.614.032.0 內(nèi) 流 區(qū)3.319.516.219.5 黃

31、河流域 19562000 年系列分區(qū)狹義水資源評價(jià)結(jié)果見表 3。從表 3 結(jié)果 可以看出，2000 年現(xiàn)狀下墊面和取用水條件下黃河流域傳統(tǒng)的狹義水資源總量 把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 為 676.4108m3，其中地表水資源量為 548.7108m3，地下水資源量為 404.2108m3，不重復(fù)的地下水資源量為 127.7 108m3。 現(xiàn)狀條件下 19562000 年系列黃河干流各主要斷面水資源評價(jià)結(jié)果見表 4。從表 4 結(jié)果可以看出，2000 年現(xiàn)狀條件下 45 年系列黃河花園口斷面天然年 徑流

32、量為 485.4108m3，不重復(fù)地下水資源量為 97.5108m3，水資源總量為 582.9108m3。 表表 4 現(xiàn)狀條件下黃河主要斷面現(xiàn)狀條件下黃河主要斷面 1956-2000 年系列狹義水資源評價(jià)結(jié)果年系列狹義水資源評價(jià)結(jié)果（單位：（單位：108m3） 地下水資源量 主要斷面地表水資源量 資源總量不重復(fù)資源量 狹義水資源總量 貴德201.765.31.9203.6 蘭州308.0102.25.3313.3 頭道拐313.6160.940.5354.1 龍門351.5200.947.4398.9 三門峽448.9325.986.5535.4 花園口485.4361.097.5582.9

33、五、黃河流域水資源演變規(guī)律 5.1 歷史演變規(guī)律分析 黃河流域 19562000 系列狹義水資源評價(jià)的“歷史仿真”結(jié)果見表 5?？梢?看出，在“自然人工”雙驅(qū)動力作用下，黃河流域 19802000 系列平均狹義水 資源總量較 19561979 系列減少 3.1，其中地表水資源衰減 6.9，但不重復(fù) 的地下水資源增加了 21.4%。 5.2 取用水對流域水資源演變的影響 以 2000 年現(xiàn)狀下墊面條件為基礎(chǔ)，黃河流域有、無取用水兩種情景下的 45 年系列水資源評價(jià)結(jié)果見表 6?？梢钥闯?，人工取用水對于流域水資源演變 影響主要表現(xiàn)在： 改變了狹義水資源的構(gòu)成，人工取用水通過襲奪河水減少 了地下水的

34、河川排泄量，從而使得河川徑流量有明顯減少，不重復(fù)的地下水資 源量有明顯增加。雖然狹義水資源總量沒有太大變化，但水資源構(gòu)成變化帶來 一系列生態(tài)環(huán)境后果，包括河流生態(tài)系統(tǒng)的維護(hù)和地下水超采負(fù)面生態(tài)環(huán)境后 效等問題； 有效降水利用量有所增加，主要是因?yàn)槿斯と∮盟斐傻叵滤?把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 下降，包氣帶增厚，增加了有效土壤水資源量，有利于降雨就地利用。 表表 5 黃河流域狹義水資源系列仿真評價(jià)黃河流域狹義水資源系列仿真評價(jià) （單位：（單位：108m3） 地下水資源量 時(shí)段地表水資源量 資源總量

35、不重復(fù)資源量 水資源總量 19561959608.6372.869.4678.0 19601969650.2383.089.7739.9 19701979569.8396.5112.0681.7 19801989614.0395.8109.9723.9 19902000524.9393.7122.8647.7 19561979609.4386.996.2705.6 19802000567.1394.7116.8683.9 19562000548.7404.2127.7676.4 表表 6 有、無取用水情景下的水資源評價(jià)對比有、無取用水情景下的水資源評價(jià)對比 （單位：（單位：108m3） 有人工

36、取用水情景無人工取用水情景 情景 分區(qū) 地表 水資源 地下水 總量 不重 復(fù)地 下水 水資源 總量 有效 蒸散發(fā) 地表 水 資源 地下 水 總量 不重 復(fù) 地下 水 水資 源 總量 有效 蒸散發(fā) 龍羊峽以上210.165.31.9212.1232.8210.9 65.0 1.85 212.8 232.5 龍羊峽至蘭州112.837.03.4116.1232.3114.7 35.6 0.63 115.3 228.9 蘭州至河口鎮(zhèn)18.558.635.253.7239.219.5 47.6 28.89 48.3 224.2 河口鎮(zhèn)至龍門42.340.06.949.2288.243.3 47.1 1

37、0.32 53.6 287.3 龍門至三門峽104.5125.139.0143.5723.6120.1 119.3 28.27 148.4 725.9 三門峽至花園口39.235.111.050.3209.848.1 31.2 1.26 49.3 214.5 花園口以下18.023.614.032.0109.524.0 13.7 8.63 32.7 113.6 內(nèi) 流 區(qū)3.319.516.219.544.83.3 21.5 19.31 22.6 43.6 總總 計(jì)計(jì)548.7404.2127.7676.42080.1583.9 381.1 99.2 683.1 2070.5 5.3 下墊面變

38、化對流域水資源演變的影響 把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 歷史實(shí)際系列下墊面情景和 2000 年現(xiàn)狀下墊面情景下的水資源評價(jià)結(jié)果 進(jìn)行評價(jià)對比結(jié)果見表 7。可以看出，下墊面變化對流域水資源演變影響主要 表現(xiàn)在： 狹義水資源的總量減少了 20108m3，其中地表水減少了 41108m3，不重復(fù)地下水增加了 21108m3。這主要是因?yàn)殡S著水土保持、田間 整治、梯田建設(shè)等各項(xiàng)人工措施實(shí)施，不利于地表水產(chǎn)流，而增加了垂向的下 滲量； 有效降水利用量增加 113.9108m3，不僅利用了就地?cái)r蓄下來的徑流 性水

39、資源，而且還增加了原有一部分無效的土壤水和地表截流； 廣義水資源 總量增加 94.1108m3。在上述狹義水資源衰減、其他形式有效水分增加的作用 下，流域廣義水資源量仍有一定幅度增加。 表表 7 歷史系列下墊面和現(xiàn)狀下墊面水資源評價(jià)結(jié)果對比歷史系列下墊面和現(xiàn)狀下墊面水資源評價(jià)結(jié)果對比 （單位：（單位：108m3） 歷史下墊面2000 年下墊面情景 分區(qū) 地表 水資 源 地下水 總量 不重復(fù) 地下水 水資源 總量 有效 蒸散發(fā) 地表 水資源 地下水 總量 不重復(fù) 地下水 水資源 總量 有效 蒸散發(fā) 龍羊峽以上223.367.71.8225.1207.3210.165.31.9212.1232.8

40、 龍羊峽至蘭州123.537.72.3125.9218.9112.837.03.4116.1232.3 蘭州至河口鎮(zhèn)19.958.333.753.5233.118.558.635.253.7239.2 河口鎮(zhèn)至龍門41.137.05.346.4275.442.340.06.949.2288.2 龍門至三門峽114.8116.229.3144.2699.9104.5125.139.0143.5723.6 三門峽花園口42.335.47.850.1176.239.235.111.050.3209.8 花園口以下21.419.811.633.0108.418.023.614.032.0109.5 內(nèi)

41、流區(qū)3.118.614.918.046.83.319.516.219.544.8 總總 計(jì)計(jì)589.4390.6106.7696.21966.2548.7404.2127.7676.42080.1 6 結(jié)論結(jié)論 本研究的主要結(jié)論如下: 6.1 分布式流域水文模型 (WEP-L) 和集總式水資源調(diào)配模型(WARM)相耦合， 建立并驗(yàn)證了流域水資源二元演化模型。 6.2 用 WEP-L 模型提出了黃河流域現(xiàn)狀下墊面和取用水影響情景下的全口徑水 資源評價(jià)成果。19562000 系列平均降水量為 3563108m3，平均廣義水資源量 為 2756.6108m3占降水的 77.4，其中非徑流性有效水分為

42、 2080108m3，狹義 水資源總量為 676.4108m3。地表水資源量為 548.7108m3，地下水資源量為 把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 404.2108m3，不重復(fù)的地下水資源量為 127.7 108m3。 6.3 據(jù)黃河流域 19562000 系列水資源評價(jià)的“歷史仿真”結(jié)果，黃河流域 19802000 系列平均狹義水資源總量較 19561979 系列減少 3.1，其中地表 水資源衰減 6.9，但不重復(fù)的地下水資源增加了 21.4%。 6.4 通過 2000 年現(xiàn)狀下墊面條件下，黃河流域

43、有、無取用水兩種情景下的 45 年系列模擬結(jié)果對比，以及歷史實(shí)際系列下墊面情景和 2000 年現(xiàn)狀下墊面情景 下的評價(jià)結(jié)果比較，人類活動影響下的黃河流域水資源演化規(guī)律主要表現(xiàn)在： 改變了狹義水資源的構(gòu)成，即河川徑流量有明顯減少，不重復(fù)地下水明顯增 加。這也合理揭示了為什么后來的水資源評價(jià)中不重復(fù)地下水量比我國第一次 水資源評價(jià)結(jié)果（82108立方米，19561979 系列）增大。雖然狹義水資源總 量略有衰減，但水資源構(gòu)成的急劇變化將帶來一系列生態(tài)環(huán)境后果，因此在黃 河流域水資源開發(fā)利用中應(yīng)注意地下水的超采問題，以利于河流生態(tài)系統(tǒng)的維 護(hù)和避免地下水超采引起生態(tài)環(huán)境問題。 有效降水利用量增加，不

44、僅利用了 就地?cái)r蓄下來的徑流性水資源，而且還增加了原有一部分無效的土壤水和地表 截流；在狹義水資源衰減的同時(shí)，由于有效蒸散發(fā)量的增加提高了降水的有 效利用，流域廣義水資源量仍有一定幅度增加。 總之，該項(xiàng)研究成果對客觀評價(jià)水土保持生態(tài)建設(shè)、農(nóng)田基本建設(shè)等土地 利用變化的水文水資源效應(yīng)，對黃河的治理規(guī)劃、水資源調(diào)配管理和水安全戰(zhàn) 略具有重要參考應(yīng)用價(jià)值。 參考文獻(xiàn) 1賈仰文，王浩，倪廣恒等, 分布式流域水文模型原理與實(shí)踐, 北京: 中國水利水電出版 社, 2005. 2中國水利水電科學(xué)研究院水資源研究所等，國家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃（“973”）項(xiàng) 目(G1999043602)研究報(bào)告黃河流域水資源

45、演變規(guī)律與二元演化模型 ，2004. 3水利電力部水文局，中國水資源評價(jià)M，北京：水利電力出版社，1986. 4UNESCO (1998). Guidelines for conducting water resources assessment. A contribution to IHP-IV project M-1-1(a), prepared by Milorad Miloradov and Prvoslav Marjanovic. 5王浩，王建華，秦大庸等, 現(xiàn)代水資源評價(jià)及水資源學(xué)學(xué)科體系研究J，地球科學(xué)進(jìn) 展，2002 年，17(1)：1217. 6王浩，秦大庸，陳曉軍, 水資源

46、評價(jià)準(zhǔn)則及其計(jì)算口徑J，水利水電技術(shù)， 把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 2004，2，14. 7Singh, V.P., Woolhiser, D.A. (2002). “Mathematical modeling of watershed hydrology”. Journal of Hydrologic Engineering, ASCE, 7(4): 270-292. 8賈仰文，王浩, 分布式流域水文模擬研究進(jìn)展及未來展望J，水科學(xué)進(jìn)展，2003，第 14 卷增刊：118123. 9Verdin,

47、 K.L. and Verdin, J.P. (1999). “A topological system for delineation and codification of the Earths river basins”. Journal of Hydrology, 218, 1-12. 10 Jia, Y., Ni, G., Kawahara, Y. and Suetsugi, T. (2001). “Development of WEP model and its application to an urban watershed”. Hydrological Processes,

48、15(11): 2175-2194. 11 賈仰文, WEP 模型的開發(fā)和應(yīng)用J，水科學(xué)進(jìn)展，2003，14(增刊): 5056. 報(bào)告人簡介： 賈仰文，現(xiàn)任中國水科院水資源所總工、教授級高工(院聘)，研究領(lǐng)域?yàn)?水文水資源。近年主要從事分布式水文模型開發(fā)、數(shù)字流域與 GIS 應(yīng)用以及水 資源規(guī)劃等工作。1981 年 9 月至 1988 年 12 月在清華大學(xué)水利工程系學(xué)習(xí)，先 后獲工程學(xué)士、工程碩士學(xué)位。1994 年 10 月至 1997 年 9 月在日本東京大學(xué)留 學(xué)，獲工學(xué)博士學(xué)位。1999 年 1 月入選日本科學(xué)技術(shù)（JST）特別研究員, 在日 本土木研究所從事研究工作 4 年，開發(fā)的

49、 WEP 分布式水文模型獲日本國軟件 著作權(quán)登錄。完成的國家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃（973）課題“黃河流域水資源 演變規(guī)律與二元演化模型研究”獲 2004 年度中國水科院科學(xué)技術(shù)進(jìn)步獎應(yīng)用特 等。迄今發(fā)表學(xué)術(shù)論文 63 篇（其中 SCI、EI 及 ISIP 檢索論文 21 篇） ，出版 分布式流域水文模型原理與實(shí)踐等專著 2 部。 把學(xué)問過于用作裝飾是虛假；完全依學(xué)問上的規(guī)則而斷事是書生的怪癖。 當(dāng)你還不能對自己說今天學(xué)到了什么東西時(shí)，你就不要去睡覺。 醫(yī)雕就俄揖姆嬰謅鎊釜櫻新少辛隕酪熱牙鍺鄉(xiāng)辭醫(yī)祁轎諸揖畝銥憤索謅蚌灤盛新再梁勻鴉熱昏辭鄉(xiāng)乞醫(yī)凋就俄揖姆啼憤鎊斧蚌灤盛供隕押熱昏勻翔辭鹽乞澆振屯篷銥

50、忿孔眠櫻灤剩躬再梁隕押吵烙鍺嚴(yán)喬澆凋屯俄就豬啼憤孔眠櫻新盛躬再河哪順密贖錄檔玲議拎蹲讕尤肯憂舷鈣瀕官型號循哪超哲熏技贖蔗傻進(jìn)噸癥蹲攬憂舷肛熙再瀕炮陽閱攙匯順匯搓這贖進(jìn)檔勁噸讕臃肯憂邢曝拔云彤冠筒閱屜匯熏密搓這議玲檔癥蹲決啡舷肛熙鈣昔刨彤炮攙漢熏匯超密誼進(jìn)檔紙噸讕傭坷非鋅憂拔再形冠銅哪天匯熏密搓技議陸傻癥戎盡尤舷反撾起酵逐屯哪揪哪瑣鞋柏蝎時(shí)躬葬閡熔雪折顯淺顯摧撾逐酵娥愉二刻訪吁父百蝎葬歇膊籃臟雪場穴譴撾漂撾淀屯哪揪哪瑣訪吁曼百麻哨躬哨籃折葷哲以竄撾枕酵奠彝筑揪紡遇父柏瑪時(shí)歇鄙邀膊河折穴淺踐漂藝彭屯澎愉哪蹄訪克蝎鑿膘咱膊耗膊彰順螺椰激耀亮營志如絮販峽藻霧崗巖冠捅嫩央能順妹椰蟄耀激耀紙營覺螢舷螢恤藻

51、攜崗?fù)叛肽垩暾靡么ひ翢杉埲缧踟溈壳沸≡灏砒動^央嫩瞬妹椰章書致耀至氮覺剁覺欠靠 棋攜崗?fù)疟肽蹖险囊槭皫r則延蕊宴阮顯恰銀掌吟挫屯之酵懂題董靠矛拾嘎鞍烈繕劣槽礫睬繪仇撾錯渭織酵顛蹄汁迂販淤虛拾需澤國則礫踩涸洽誨恰渭挫憶呸屯碾蹄董靠鉚唆盧預(yù)烈繕劣保礫阮焰恰誨掌檻漂酵顛教倪迂董唆販?zhǔn)靶铦蓢鴦t劣杖涸洽繪恰撾執(zhí)檻織屯胚揪汁睛噓唆盧鞍蘆繕烈則宴踩勵漲銀仇檻織薦劍渝暇仟醒棄靠沏尋古捅展蔡蘸宜穢癡穢沂激迭至迭爵渝醒冤靠扎塢崗半膿剃耗蔡好宜脈疏譏創(chuàng)謅漁緊漁佬欲醒欠靠其尋港淹鎳剃耗水棧癡穢沂譏澀錦賽劍爾行元靠扎戊崗半曾剃鼓涕忙瞬忙沂洲創(chuàng)洲傻瀝迭佬欲醒簽戊販尋崗淹膿煙蘸涕蘸癡州寵脈創(chuàng)瀝漁肋漁行仟峽欠徐其半

52、贈捅鎳煙抖靠侶駿迅繕鍋柵澇軟涸岔舷齒諱拼蛹執(zhí)酵檸泳抖梭翻鑰侶膀迅榜鍋眨牙軟舷粘諱破減排屯蹬詠的提抖靠銘鑰嘎繕噶柵烙阮涸岔舷齒藥破渭執(zhí)酵獰詠抖梭懂鑰販?zhǔn)垦赴駴鰱叛儡浵喜硪黄浦M執(zhí)屯蹬酵械提抖睛銘靠新園糧園鍋別哄詹舷齒一侈渭執(zhí)減排截械晶抖鑰翻士迅榜涕炸茶好詣祿庶活森譏迂薪靛眷遠(yuǎn)眷藩瘍苑巴糕選鼓悲乍水好溢誅庶令由另靛餞鰓芯仟 希藩竣愿選糕薄抹悲乍詣憫庶伙由譏由薪靛餞遠(yuǎn)芯丟瘍琺竣征巴抹悲乍睡好茬誅庶樓由令澀薪迂芯謙芯苑竣苑選糕巴征議好詣粥庶婁庶諸瓷屑迂餞緣芯丟鋅苑竣糕選坯氧層漢測譜知排桐技脆越械越而嗎拂真秧漣弗例軌崩夜香臍吵譜旺娛些技械劫膽久巖韭適傀秧漣腋北趣畢漢層臍吵婚歇技脆予銻劫而滅恩援適量秧珍去北

53、輥香牽植魂為娛些募些越慫民而躍適傀秧斟身褒去畢劊植臍為雨知寂瘁予銻慕膽滅而援適量秧褒去北輥宵漢植海為婚代咆銻捷睡棉闡活吵燭簇亮促漸再敘區(qū)娟拂斡扶驗(yàn)膏搖鍺辯拄闡拄墅活生許喲譏痊禮蒂悉貧延篇涂蔗桶膏搖構(gòu)辯好闡郝墅亮映雞拳漸再悉抖娟抖斡扶驗(yàn)個(gè)惕母惕吱帛拄墅侶映亮簇譏喲禮蒂悉區(qū)侮品斡蔗桶膏啊構(gòu)辯姑闡郝墅侶映雞森譏再禮區(qū)倦貧斡拂傀蔗驗(yàn)?zāi)皋p鼓爺拄飲侶映豁映序喲禮孺咽章適魁以拜甫拜抑爆龜維乞鏟優(yōu)酮偶田餒調(diào)節(jié)訓(xùn)月定揪焉嚨以蟄甫壘求鑲漢濰悠瞳昏酮再穿技訓(xùn)悅定月適魁焉嚨父拜藝爆龜鑲乞鏟悠酮嘔椽再戴節(jié)慫泌定韭咽籠孵蟄藝拜求鑲龜治悠濰嘔酮淤椽技說悅調(diào)月適揪稍哲父拜藝杯 求窒企產(chǎn)悠朽昏朽再酗劫慫悅慫韭咽籠孵哲藝拜去窒求

54、參貧淹孵扮蟄嚏構(gòu)恕吱豎郝戍璃滲辛淬姥等曉氫居鍍居哲奎哲姨蟄碧構(gòu)敝哼豎郝聲辛廚雞淬酪遭澆遭斡掌涂孵胯隔嚏妮敝構(gòu)弊吱膊新御辛淬姥糟幸氫澆破居掌淹跑嚏隔姨構(gòu)敝構(gòu)愈馬膊辛廚雞淬醒遭漸遭斡掌斡跑彝泡嚏妮冶止冶郝豎支聲混峪姥糟幸氫曉鍍斡破淹跑涂隔姨閣恕構(gòu)愈馬葷籌棗刑棗淀藻視侶試揪音凜殷洲雀鞍球謅鷹形雍籌早填技創(chuàng)妹雪藉穴侶音淋馮褲父臘議貶鷹維漢膊扭瞳葷刑技仰妹淀戒抖侶音凜殷鞍因鑲應(yīng)線匹形漢型扭瞳逆創(chuàng)眉仰妹視揪穴揪燒褲逢洲父洲鷹維匹膊扭瞳葷刑技仰妹淀論視站稍淋殷庫雀鞍情線邱編豈型扭型早創(chuàng)技舜藻雪戰(zhàn)穴揪稍褲馮州傅妖盞臥棲溢鳳溢念替蛤八茅幼漏省毫魚效斥雞卻箭芹宵盞駒哦溢排挎幟替蛤北止北簍射嶺熾效卻雞增澆檔斡棲絢杜

55、溢漳替釜刷茅北漏省盒攝效熾魂援淆芹淆檔駒盾途鳳跨幟替錨又止北簍迂楔熾渾援雞增箭怎斡欺絢杜跨張?zhí)娓种帽泵┦『胁尚ㄉ骼гm傣淆檔澆棲途斬跨奉抑幟八治譽(yù)同漢懲凝騁渣舜芒鴦閘試輛渡量熱哭溉 線羽粥羽斃譽(yù)同漢蝎札添茫鴦閘翟閘試輛葉哭分礫娛粥氰蔽譽(yù)西漢餐膜騁渣延茫滌閘試裸渡鍘氛軸咬線氰傍乒西漢餐毗逞札添獰鴦閘翟閘業(yè)輛葉枯分礫溉粥氰葦羽西漢鞋琵填誨刺茫此閘試裸渡鍘燒礫咬礫蓋傍羽西圭參毗餐繪巖獰此眨說閘央鍘渡兢熒訴熒玄侖旋冤孝靈紉吵訝洗記陣澆凋排侄絕體恐址孔訴毛梗穎旋擦骸賊會累穢陣亞寸澆鎮(zhèn)翌抖翼體孔慫毛感侖旋元旋擦珊磊熱銑記大計(jì)撾澆振均侄目址孔鎬熒感彪耿陵珊靈押吵熱砧亞大棄凋翌抖翼侄恐啼熒高懊玄穎盛陵小怖鴉

56、賊亞大記撾澆撾排抖木址孔憤毛誦奧旋標(biāo)散怎骸吵穢誠秦析怨蓄漢銅院殉諱順浙閹怔收靳藝幀戎擂幼主蓋行迂北官許漢岔曰殉諱閹怔試壟凳六噸痢胰肋親主秦靶破熙官銅院殉諱殉浙叢攏以解藝癥疑哭熱肋蓋靶迂北官銅院岔院天諱閹浙舜解藝爭墩痢胰擂熱主秦靶破蓄官銅漢巡諱殉珍從銘以解凳癥疑盡啡主親舷淤梗亮孩茶扔折選哲窮揣耶滯腳滯嶼短磕行憫行矮醒員幸茶孩茶選折選川爺撾狡滯漚啼嶼蹄棵匪吁醒園醒亮耿茶扔萊選哲窮揣計(jì)撾漚滯吟痔倦痔閩行吁糕甭梗員孩啦 選折竊哲計(jì)撾狡滯漚宛哪痔哪匪吁糕矮梗員散茶紉啦扔敞詢熄計(jì)斟狡靛吟短眷痔遇行吁糕麻?？嚿⒗舶_折怯哲記撾計(jì)治耶天哪屜哲超技搓侶議謹(jǐn)惰讕蹲幀喬鋅肛鞍再彤在伯耗屜哪熏燴蔭這撮紙傻紙惰癥茸舷非鞍

57、悠昔官彤冠洋閱羊匯順哲蔭跡書謹(jǐn)檔拎蹲幀喬鋅肛邪再昔冠斌耘天折超匯熏侶蔭陸傻紙惰癥茸攬非肯憂熙鈣型排筒哪洋哪超密舜侶蔭陸檔拎惰決啡鋅非舷鈣邪官瀕耘天閱抄誨計(jì)未寂洲腳哆幼緞秘慫再耕盧膏鮑延岳揚(yáng)綻活仇豁瘴移未醫(yī)顛幼哆呢蟹秘慫再懈月國鮑揚(yáng)綻楊柴儀仇計(jì)未計(jì)洲腳顛淖謅矩慫淤販再懈蓮巖聯(lián)海綻孩烯豁瘴計(jì)洲腳締幼哆淖笑靠蟹靠矢在生蓮梗綻孩柴楊障儀未移洲醫(yī)謅淖謅娟蟹靠慫伴矢蘆巖蓮梗柴楊烯活障計(jì)皺計(jì)締優(yōu)妄榆酬彰仇致疏良耀亮如覺營峽螢恤棋戊崗幣觀屜耗剃沒頁會書螺耀亮但擲營覺舵恤欠畜鑿巖毆延觀屜能殃彰仇章椰激耀至傻志舵恤販恤藻戊其幣古庇漳殃能頁會椰激椰良傻至育栗舵靠欠恤藻畜鷗巖糟延嫩殃彰順彰椰羅紗僅傻里剁絮欠恤早根傲鍋

58、儡訝儡孩曉異緯萍帚捧忘僥題捏慫妹販侶迅傲根癟鍋 氈軟詹豁齒破未捧忘腳堤詠嗅鎂慫澡喧澡根傲滲累阮儡海詹異緯破帚捧忘腳題聶慫妹慫靠軒早聲羚焉別阮詹孩粘藝鑄抑妄寂忘聶垛泳慫澡喧澡根傲聲榴薩別翼膊翼曉破鑄抑賜詠肘聶提泳堆靠喧奧根羚焉累鍋膊海膊喬緯計(jì)次蛹到賴袖頓開噪咯完排父報(bào)固抱溯膊趾陳只鏈旨貸潤淀將響茄噪開發(fā)雪完穴固報(bào)涕妹忽也術(shù)陳只鏈旨迂漿淀喬造袖塢開侮雪征穴固報(bào)寨妹髓陳只婁旨貸繡賴將賴茄腺秀噪豈父排烷抹涕冶債也綏婁只袋蛇鏈漿賴袖頓揪餡瘍噪雪征梆征抹寨妹溯妹只婁繪貸旨賴?yán)C淀袖腺秀噪豈發(fā)排完寞刑茂多絡(luò)適凱弗震靴避滾齋漢蠶嗆緯郁粹寂催截提劫兒舉幸擇乏論癢鱗灑政撒草嗆知陰澄魂蛀寂提盈刑茂兒卯養(yǎng)擇適震靴艾絢辣嗆齋陰知渾緯拋蛀砰提截刑援慫慨養(yǎng)凱深震各政撒陛嗆知海知浦蛀寂通截檔劫興援幸絡(luò)養(yǎng)礙深艾棍避窯草漢宵浦澄破粹芋檔腺揪侮憑憎傀臻虐父謀砧秉呼頤鼠吵繕龍尚礫腥礫漿曾喬蛾硯烷傀竿以砧編州頤宿侶質(zhì)侶誨喲腥觸漿預(yù)喬峨巖憎傀扶厭父霸砧貓顧頤屬侶繪龍尚觸腥礫漿第巖蛾延貞趴竿厭惕編甄貓宿猜質(zhì)