水文模型薩克拉門(mén)托模型

薩克拉門(mén)托流域水文模型現(xiàn)在是1頁(yè)\一共有34頁(yè)\編輯于星期三薩克拉門(mén)托流域水文模型薩克拉門(mén)托流域水文模型是美國(guó)天氣局水文辦公室薩克拉門(mén)托預(yù)報(bào)中心、和于70年代初期在第IV號(hào)斯坦福模型基礎(chǔ)上改進(jìn)和發(fā)展的。1973年研制成功了日流量模擬程序，1975年又進(jìn)一步提出了6h時(shí)段模擬程序，因它始用于薩克拉門(mén)托河而得名該模型功能較完善，能應(yīng)用于大、中流域，又能適應(yīng)濕潤(rùn)地區(qū)和干旱地區(qū)，它在美國(guó)的水文預(yù)報(bào)中得到了較廣泛的應(yīng)用，也是我國(guó)引進(jìn)的模型中人們較為熟悉的模型之一?，F(xiàn)在是2頁(yè)\一共有34頁(yè)\編輯于星期三薩克拉門(mén)托河薩克拉門(mén)托河位于美國(guó)加利福尼亞州中部，是該州最長(zhǎng)的河流，它發(fā)源于卡斯克德山脈的沙斯特山附近，最后注入太平洋舊金山灣北支舒森灣。河流全長(zhǎng)615km，總流域面積約為7.02萬(wàn)平方公里，多年平均流量856立方米/秒，多年平均徑流量270億立方米。薩克拉門(mén)托河流域所在地區(qū)為薩克拉門(mén)托河谷與圣華金河谷，均是中央河谷的一部分，其北部為喀斯喀特山脈，東部為內(nèi)華達(dá)山脈，西部為沿海山嶺。中央河谷南北長(zhǎng)約700km，東西寬約90km，面積約4.5萬(wàn)平方公里，是美國(guó)著名的農(nóng)業(yè)地帶。河谷北部年降水量約為1000mm，南部?jī)H150mm，流域內(nèi)年徑流量3/4集中在北部。現(xiàn)在是3頁(yè)\一共有34頁(yè)\編輯于星期三薩克拉門(mén)托河現(xiàn)在是4頁(yè)\一共有34頁(yè)\編輯于星期三基本概念薩克拉門(mén)托模型是集總式連續(xù)運(yùn)算的確定性流域水文模型，是以土壤含水量貯存、滲透、排水和蒸散發(fā)特性為基礎(chǔ)來(lái)模擬水溫循環(huán)的綜合的河川徑流流域模型。薩克拉門(mén)托模型主要由降雨、融雪、蒸散發(fā)、入滲、河網(wǎng)總?cè)肓骷捌浜泳W(wǎng)匯流六個(gè)部分組成，其中，土壤含水量計(jì)算模型是該系統(tǒng)的中心組件。現(xiàn)在是5頁(yè)\一共有34頁(yè)\編輯于星期三通用水文模型現(xiàn)在是6頁(yè)\一共有34頁(yè)\編輯于星期三模型結(jié)構(gòu)圖現(xiàn)在是7頁(yè)\一共有34頁(yè)\編輯于星期三土壤含水量模型降雨不透水面透水面上層下層滲透方程上層土壤極其貯水量下層土壤極其貯水量充滿上層張力水，產(chǎn)生上層自由水滲透到下層，形成下層張力水、淺層以及深層自由水現(xiàn)在是8頁(yè)\一共有34頁(yè)\編輯于星期三模型參數(shù)薩克拉門(mén)托模型的參數(shù)，除蒸發(fā)和匯流部分之外，土壤含水量模型中還有17個(gè)參數(shù)?，F(xiàn)在是9頁(yè)\一共有34頁(yè)\編輯于星期三模型參數(shù)現(xiàn)在是10頁(yè)\一共有34頁(yè)\編輯于星期三水源劃分及徑流計(jì)算壤中流淺層基流深層基流直接徑流地面徑流現(xiàn)在是11頁(yè)\一共有34頁(yè)\編輯于星期三水源劃分及徑流計(jì)算1、直接徑流（DRO）又分永久性不透水面積上產(chǎn)生的直接徑流（ROIMP）和暫時(shí)性不透水面積上產(chǎn)生的直接徑流(SDRO）。1）永久不透水面積上形成的直接徑流模型設(shè)想雨水降落在湖面、沼澤或河網(wǎng)等以及其毗連的永久性不透水面積上時(shí)，都產(chǎn)生直接徑流，所以ROIMP=PX.PCTIMPX——降水量PCTIM——永久性不透水面積比重2）可變不透水面積上形成得直接徑流因?yàn)樵诳勺儾煌杆娣e上也分為上下兩個(gè)土層，各層的張力水容量與透水面積上的一樣，但不設(shè)自由水。總的張力水蓄量由上下兩層的張力水蓄量組成，上層張力水蓄量等于透水面積上的上土層張力水蓄量,下土層張力水蓄量為ADIMC-UZTWC?，F(xiàn)在是12頁(yè)\一共有34頁(yè)\編輯于星期三水源劃分及徑流計(jì)算現(xiàn)在是13頁(yè)\一共有34頁(yè)\編輯于星期三水源劃分及徑流計(jì)算現(xiàn)在是14頁(yè)\一共有34頁(yè)\編輯于星期三滲透計(jì)算下層水分飽和時(shí)水分從上層滲透到下層的滲透率為飽和滲透率（PBASE）它是兩種下層自由水貯存量與其相應(yīng)的排水率乘積之和：現(xiàn)在是15頁(yè)\一共有34頁(yè)\編輯于星期三滲透計(jì)算現(xiàn)在是16頁(yè)\一共有34頁(yè)\編輯于星期三滲透計(jì)算現(xiàn)在是17頁(yè)\一共有34頁(yè)\編輯于星期三滲透計(jì)算滲透公式中，滲透量是按一定的比例（PFREE）分配給下層張力水和下層自由水的。而分配給下層自由水的滲透量，又按淺層自由水和深層自由水的含水比（RATLS,RATLP）分配給淺層自由水和深層自由水。現(xiàn)在是18頁(yè)\一共有34頁(yè)\編輯于星期三土壤水的水平交換現(xiàn)在是19頁(yè)\一共有34頁(yè)\編輯于星期三蒸散發(fā)計(jì)算(1)上土層張力水蒸散發(fā)量降落在透水面積上的時(shí)段雨量,首先補(bǔ)充上土層張力水蓄量，當(dāng)滿足上土層張力水的缺水量后其余的雨量成為有效降雨。UZTWM～為上土層張力水容量（mm）；UZTWC～為上土層張力水蓄量（mm）。流域的蒸散發(fā)能力EP。(2)上土層自由水蒸散發(fā)量UZFWC～為上層自由水蓄量(mm)現(xiàn)在是20頁(yè)\一共有34頁(yè)\編輯于星期三蒸散發(fā)計(jì)算(3)下土層張力水蒸散發(fā)量

UZTWM～為上土層張力水容量(mm)；LZTWM～為下土層張力水容量(mm)。LZTWC~為下層張力水蓄量。（4）水面蒸發(fā)量

SARVA～為河網(wǎng)、湖泊及水生植物的面積占全流域面積的百分?jǐn)?shù)；PCTIM～為永久不透水面積占全流域面積的百分?jǐn)?shù)。現(xiàn)在是21頁(yè)\一共有34頁(yè)\編輯于星期三蒸散發(fā)計(jì)算

(5)可變不透水面積上的蒸散發(fā)量UZTWM～為上土層張力水容量(mm)；LZTWM～為下土層張力水容量(mm)?？偟膹埩λ盍緼DIMC由上下兩層的張力水蓄量組成，流域的蒸散發(fā)能力EP?，F(xiàn)在是22頁(yè)\一共有34頁(yè)\編輯于星期三模型計(jì)算框圖開(kāi)始輸入?yún)?shù)、各種初值、降水量、蒸發(fā)能力各時(shí)段降水、蒸發(fā)取值計(jì)算ROIMP，各種蒸發(fā)量、調(diào)整各層土壤含水量劃分步長(zhǎng)，得6h時(shí)段內(nèi)步長(zhǎng)NINC,步長(zhǎng)內(nèi)有效降雨量PINC將日出流系數(shù)換為步長(zhǎng)出流系數(shù)按每步長(zhǎng)計(jì)算滲透量及產(chǎn)流量計(jì)算1.SDRO2.地下徑流3.滲透量4.壤中流5.地面徑流累加得徑流時(shí)段量輸出結(jié)果現(xiàn)在是23頁(yè)\一共有34頁(yè)\編輯于星期三模型優(yōu)缺點(diǎn)及應(yīng)用情況優(yōu)點(diǎn)：模型的參數(shù)、所用的變量和模擬的過(guò)程具有一定物理意義，易于理解和便于根據(jù)實(shí)測(cè)的降水、流量及流域特征資料估算其初值。缺點(diǎn)：參數(shù)眾多，調(diào)試?yán)щy該模型在美國(guó)應(yīng)用較多，在我國(guó)應(yīng)用不多。該模型經(jīng)過(guò)改進(jìn)可用于寒冷濕潤(rùn)地區(qū)的徑流模擬現(xiàn)在是24頁(yè)\一共有34頁(yè)\編輯于星期三現(xiàn)在是25頁(yè)\一共有34頁(yè)\編輯于星期三模型結(jié)構(gòu)圖現(xiàn)在是26頁(yè)\一共有34頁(yè)\編輯于星期三土壤含水量模型降雨不透水面透水面上層下層滲透方程上層土壤極其貯水量下層土壤極其貯水量充滿上層張力水，產(chǎn)生上層自由水滲透到下層，形成下層張力水、淺層以及深層自由水現(xiàn)在是27頁(yè)\一共有34頁(yè)\編輯于星期三徑流直接徑流壤中流地面徑流基流永久不透水面可變不透水面1.不透水面積與下層張力水相對(duì)蓄水量平方成正比2.產(chǎn)生的直接徑流透水面積暫時(shí)不透水面透水部分現(xiàn)在是28頁(yè)\一共有34頁(yè)\編輯于星期三下層水分飽和時(shí)水分從上層滲透到下層的滲透率為飽和滲透率（PBASE）它是兩種下層自由水貯存量與其相應(yīng)的排水率乘積之和：滲透計(jì)算2.考慮下層土壤含水量缺水率(DEFR)影響3.考慮有效供水量影響:水量不充足時(shí)為上層自由水含水量現(xiàn)在是29頁(yè)\一共有34頁(yè)\編輯于星期三滲透計(jì)算現(xiàn)在是30頁(yè)\一共有34頁(yè)\編輯于星期三滲透計(jì)算現(xiàn)在是31頁(yè)\一共有34頁(yè)\編輯于星期三蒸散發(fā)計(jì)算土壤各部分蒸散發(fā)是隨蒸散發(fā)能力及張力水貯存的分布而不同，上層蒸發(fā)的水分，其速率依張力水貯存量而定，下層速率取決于剩余蒸散發(fā)能力乘以下層張力水含水量與張力水總?cè)萘恐?。若含水量不足則全部蒸發(fā)。蒸散發(fā)上層張力水蒸散發(fā)量上層自由水蒸散發(fā)量下層張力水蒸散發(fā)量河湖面積上蒸發(fā)量和水生植物面積蒸散發(fā)量可變不透水面積上的蒸散發(fā)量現(xiàn)在是32頁(yè)\一共有34頁(yè)\編輯于星期三土壤水平交換蒸發(fā)后，如果自由水含量高于張力水，則由自由水補(bǔ)充張力水現(xiàn)在是33頁(yè)\一共有34頁(yè)\編輯于星期三模型計(jì)算框圖開(kāi)始輸入?yún)?shù)、各種初值、降水量、蒸發(fā)能力