SWAT培訓(xùn)文檔-數(shù)據(jù)庫建立等

第一屆BNU-SWAT&SWAT-CUP培訓(xùn)2011年10月14-17日第一部分sWAT模型的建模及應(yīng)用框架下午1600-18:00第一講培訓(xùn)框架介紹Karim教授/魚京善教授上午8:30-11:30魚京善教授第三講DEM的利用魚京善教授魚京善教授魚京善教授上機(jī)操作魚京善教授孫文超博士下午13:30-17:00孫文超博士孫文超博士孫文超博士孫文超博士魚京善教授上機(jī)操作孫文超博士李占杰老師第三部分BNU-SWAT的應(yīng)用上午8:30-11:30魚京善教授第十二講BNU-SWAT模型率定魚京善教授魚京善教授第十四講降水及氣溫變化響應(yīng)模擬魚京善教授第十五講土地利用及作物變化模擬魚京善教授第十六講BNU-SWAT地下水埋深模擬魚京善教授第十七講WEC&WQS模型介紹魚京善教授上機(jī)操作孫文超博士李占杰老師下午13:30-17:00晚上19:00-21:30第五部分SWAT-CUP的應(yīng)用第一部分SWAT模型的建模及應(yīng)用框架BNU-SWAT使用手冊(cè)BNU-SWAT&SWAT-CUP培訓(xùn)教材1.新建BNU-SWAT工程使用BNU-SWATV1.0進(jìn)行模型移植時(shí)，必須確認(rèn)ArcSWAT工程是使用2.1(Beta)版2.TxtInOut文件的復(fù)制首先，刪除BNU-SWAT新工程的【Bin/Scenarios/Default/TxtInOut】文件夾內(nèi)的所有3.數(shù)據(jù)庫文件的復(fù)制制到BNU-SWAT新工程【Bin/Scenarios/Default】文件夾中，改名為input.mdb,并覆蓋原ArcSWAT工程用的SWAT2005.mdb文件復(fù)制到BNU-SWAT新工程的【Bin】文件夾4.確認(rèn)工程路徑雙擊新文件夾下的SWAT_CN.EXE文件啟動(dòng)BNU-SWAT工程，單擊按鈕進(jìn)入路徑設(shè)置菜單，確認(rèn)工程路徑是否指向新工程文件夾。確認(rèn)之后，關(guān)閉BNU-SWAT工程。1BNU-SWAT&SWAT-CUP培訓(xùn)教材用戶及畫面設(shè)置用戶及畫面設(shè)置使用人&預(yù)五 北京師范大學(xué)水科學(xué)研究院數(shù)字流域?qū)嶒?yàn)室2BNU-SWAT&SWAT-CUP培訓(xùn)教材Format:FLogtoFiF:u?ers\ingshanlappdatal\lecal\ta(注意：格式轉(zhuǎn)換時(shí)，文件夾名稱中不能帶有中文)(3)將調(diào)整好的TAB格式文件復(fù)制到BNU-SWAT新工程的TAB文件夾內(nèi)，覆蓋原有文件；啟動(dòng)BNU-SWAT工程，【流域分區(qū)】模塊中，打開利用平移工具調(diào)整GIS畫面北京師范大學(xué)水科學(xué)研究院數(shù)字流域?qū)嶒?yàn)室3 圖6BNU-SWAT中畫面位置的調(diào)整和保存(4)BNU-SWATGIS畫面的右上角單擊鼠標(biāo)右鍵，打開GIS屬性設(shè)置窗口，【CurrentTool】依次選為ZoomIn,ZoomOut及Pan等，單擊【確定按鈕】,移動(dòng)畫面到適當(dāng)位置后，再次單擊【保存】工具，關(guān)閉BNU-SWAT工程。顏色應(yīng)用A) 圖7BNU-SWAT中右上角GIS畫面位置的調(diào)整4 BNU-SWAT&SWAT-CUP培訓(xùn)教材流域分區(qū)數(shù)據(jù)輸入運(yùn)行模型結(jié)果分析刻數(shù)率定數(shù)據(jù)庫路徑D:ABNUSWATTest貴陽SW目標(biāo)數(shù)據(jù)庫土境屬性化肥使用條蟲劑使用城鎮(zhèn)用地屬性擬系統(tǒng)設(shè)置系統(tǒng)數(shù)程庫土值水設(shè)置溫懂結(jié)構(gòu)槽被名稱農(nóng)業(yè)用地代碼AGRL植被分類暖季一年生消號(hào)1DLAl0eaunit/heaurit] 槽被名稱農(nóng)業(yè)用地代碼AGRL植被分類暖季一年生消號(hào)1DLAl0eaunit/heaurit] FgPkgbiomsVPDfRBMJBEESlioethal05題附脫表5CS徑流曲故費(fèi)寧系數(shù)率)4BD23455565789756666604高日蒙6s5田差5高35得麥5高稻2g型運(yùn)行正常。北京師范大學(xué)水科學(xué)研究院數(shù)字流域?qū)嶒?yàn)室51系統(tǒng)啟動(dòng) 12用戶及畫面設(shè)置 13流域分區(qū) 2 54.1系統(tǒng)數(shù)據(jù)庫設(shè)置功能 5 84.3地下水設(shè)置功能 95數(shù)據(jù)輸入模塊 5.3排污口 5.4用水?dāng)?shù)據(jù) 5.5水庫數(shù)據(jù) 6.1運(yùn)行模型功能 6.2情景設(shè)置功能 7結(jié)果分析 7.1子流域按月 27.3河段按月 8參數(shù)率定模塊 8.1待率定參數(shù) 8.3子流域、土地利用、土壤類型 268.4待率定參數(shù)的意義及取值范圍 9.1降水氣溫模塊 9.3農(nóng)業(yè)管理模塊 附錄A不同計(jì)算結(jié)果的對(duì)比方法 附錄B雨量站降水?dāng)?shù)據(jù)和統(tǒng)計(jì)年鑒降水?dāng)?shù)據(jù)的統(tǒng)一方法 40附錄C氣象數(shù)據(jù)的編輯和導(dǎo)入方法 SWAT-CN啟動(dòng)界面如圖1-1所示，主要功能依次為流域分區(qū)、數(shù)據(jù)輸入、北京市縣級(jí)SWAT模型應(yīng)用系統(tǒng)流域分區(qū)數(shù)據(jù)輸入云行模結(jié)果分析系統(tǒng)設(shè)置模型啟動(dòng)后，點(diǎn)擊主菜單上的笠按鈕，可進(jìn)行SWAT-CN模型的用戶及界圖2-1所示。用戶及畫面設(shè)置用戶及畫面設(shè)置透明效果不透明圖2-1用戶及畫面透明度設(shè)置選中?透明效果2控制顯示選項(xiàng)控制顯示選項(xiàng)3打開降水(氣溫)數(shù)據(jù)輸出圖片(B)通州區(qū)大興區(qū)平谷區(qū)房山區(qū)密云縣專題圖(I)工具欄包含多個(gè)選項(xiàng)可以對(duì)地圖進(jìn)行自由的放大，縮小，移動(dòng)，測(cè)量及保存。具有Mapinfo的圖層顯示功能，單擊可以對(duì)圖層進(jìn)行添加，刪除，顯示，可選，可編輯等多種操作。區(qū)縣名稱區(qū)縣界線區(qū)縣界面r窗體界面右側(cè)選項(xiàng)用于顯示SWAT圖層數(shù)據(jù)。包括3大選項(xiàng)：基本及專題圖層選項(xiàng)，數(shù)據(jù)選擇及時(shí)間選擇。單擊基本按鈕，用戶可以從下拉框中選擇基本圖層在窗體中進(jìn)行顯示?；緢D層如圖3-4。圖3-5從下往上依次為土壤類型圖，土地利用圖，高程圖和不顯示圖層。4單擊專題選項(xiàng)，用戶可以使用數(shù)據(jù)選擇及時(shí)間選擇選項(xiàng)對(duì)SWAT模擬結(jié)果一專題圖的形式在窗體中顯示出來。數(shù)據(jù)選項(xiàng)包括融雪，ET,潛在ET等8項(xiàng)。時(shí)間選項(xiàng)包括年選項(xiàng)和月選項(xiàng)兩項(xiàng)。融雪土壤含水根系層滲漏坡面流地下徑流徑流量年2001年2002年2004年2005年2月3月4月5月6月7月時(shí)間選擇=5圖3-7專題圖層樣例4系統(tǒng)設(shè)置模塊應(yīng)北京市縣應(yīng)系統(tǒng)設(shè)置流域分區(qū)系統(tǒng)設(shè)置圖4-1系統(tǒng)設(shè)置菜單界面6編編輯STAI數(shù)據(jù)庫文件確定圖4-2系統(tǒng)數(shù)據(jù)庫設(shè)置界面本模塊允許用戶新建，保存，導(dǎo)入和導(dǎo)出氣象站，并在導(dǎo)入的情況下對(duì)里面的具7數(shù)據(jù)庫寫入按鈕數(shù)據(jù)庫寫入按鈕植被類型數(shù)目：95植被類型數(shù)目：95IC四2345B967889696596994999200053949399B44亮箱5梯牧草陸尾草)693植被頑蓋植被類型參數(shù)植被名稱農(nóng)業(yè)用地33.5455071212unts46是寧系數(shù)斷率)ABD圖4-3土地利用數(shù)據(jù)庫設(shè)置界面完成最終的更新之后，退出本界面再重新加載，之后就可以看到新的氣象站點(diǎn)已SVAT數(shù)據(jù)庫的編輯-氣象站點(diǎn)參鼓 經(jīng)度度)緯度度)高程(m)數(shù)據(jù)年數(shù)月4月5月8月9月10月12月日均降水降水標(biāo)準(zhǔn)差降水協(xié)方差降水天數(shù)日均太陽輻射圖4-4氣象站點(diǎn)新建界面82)氣象站的數(shù)據(jù)修改氣象站TongZ導(dǎo)入導(dǎo)出經(jīng)度度度)高程(而11月12月日均降水降水標(biāo)準(zhǔn)差降水協(xié)方差降水天數(shù)日均太陽輻射3)子流域關(guān)聯(lián)功能9系統(tǒng)參數(shù)設(shè)置系統(tǒng)參數(shù)設(shè)置確定取消圖4-6土壤參數(shù)設(shè)置4.3地下水設(shè)置功能這時(shí)可以看到在左邊的文本框中出現(xiàn)了兩個(gè)數(shù)字，均為1000這表示當(dāng)前的參數(shù)取值為1000,單位為毫米(mm),修改“取值”文本框中的值，改為1001,點(diǎn)擊賦值按鈕，退出然后再打開就會(huì)發(fā)現(xiàn)，值已經(jīng)被修改，若選擇列表框中的“全部”選項(xiàng)，則表示對(duì)相關(guān)數(shù)據(jù)類型的全部執(zhí)行操作，建議在模擬前就設(shè)定好該功能，盡量采用精細(xì)設(shè)置的方法逐一確定，其他參數(shù)的操作與之類似。取值范圍：1005000mm子流域3參數(shù)修改區(qū)域上壤類型賦值退出圖4-7地下水設(shè)置模塊5數(shù)據(jù)輸入模塊數(shù)據(jù)輸入模塊是SWAT_CN為了方便用戶進(jìn)行操作而設(shè)計(jì)的數(shù)據(jù)輸入接口，據(jù)、氣溫?cái)?shù)據(jù)及用水資料。有關(guān)本章使用的數(shù)據(jù)結(jié)構(gòu)CSV文件結(jié)構(gòu)，請(qǐng)參照附錄A中相關(guān)說明。水圖5-1系統(tǒng)界面圖數(shù)據(jù)。點(diǎn)擊“查找”按鈕，系統(tǒng)會(huì)彈出對(duì)話框，如圖5-2所示，用戶可以選擇降溫?cái)?shù)據(jù)后面的“查找”按鈕來導(dǎo)入氣溫?cái)?shù)據(jù)，最終結(jié)果如圖5-4所示。目目設(shè)置氣象數(shù)據(jù)降水?dāng)?shù)據(jù)查找圖5-2氣象數(shù)據(jù)設(shè)置界面一當(dāng)用戶最終完成所有設(shè)置之后，點(diǎn)擊“導(dǎo)入”按鈕即可完成最終設(shè)置結(jié)果的開始1996-01-01結(jié)束：2005-12-31站點(diǎn)數(shù)：1小工具pOp>Cvstmp>Cvs按統(tǒng)計(jì)年鑒修改降水查找5.2入境水量5.3排污口5.4用水?dāng)?shù)據(jù)導(dǎo)出的數(shù)據(jù)格式與氣象數(shù)據(jù)相同，均為CSV格式，相關(guān)說明參見附錄A。萬m3子流域2子流域3子流域4子流域5子流域7子流域8導(dǎo)入導(dǎo)出賦值退出10月11月12月圖5-5用水?dāng)?shù)據(jù)模塊數(shù)據(jù)的修改工作非常便捷，用戶只需要單擊要修改的數(shù)值，利用鍵盤上的件名就可以完成數(shù)據(jù)的導(dǎo)入導(dǎo)出工作。(注：所有設(shè)置完成之后需要重新運(yùn)行模導(dǎo)入用水?dāng)?shù)據(jù)導(dǎo)入用水?dāng)?shù)據(jù)2004年用水?dāng)?shù)據(jù).esvpep.csv我的電腦文件名Q廠以只讀方式打開B)圖5-6用水?dāng)?shù)據(jù)導(dǎo)入模塊6運(yùn)行模型模塊輸入之后，只有最終完成模型的運(yùn)行才可以生成相應(yīng)的結(jié)果，該模塊包含兩個(gè)功莫型應(yīng)思莫型應(yīng)思退系統(tǒng)設(shè)置6.1運(yùn)行模型功能運(yùn)行SWAT運(yùn)行時(shí)間運(yùn)行運(yùn)行進(jìn)度?輸出HRU運(yùn)行進(jìn)度圖6-2運(yùn)行模型設(shè)置界面據(jù)，調(diào)整模擬的時(shí)間(注意必須在降水范圍內(nèi)調(diào)整)當(dāng)運(yùn)行百分率從0%達(dá)到100%之后，點(diǎn)擊生成數(shù)據(jù)庫，將相關(guān)數(shù)據(jù)輸出到保存情景功能：是指當(dāng)模型運(yùn)行完1次之后，點(diǎn)擊按鈕，系統(tǒng)將記憶相關(guān)參6.2情景設(shè)置功能圖6-3情景功能設(shè)置界面情景設(shè)置功能可以將整個(gè)運(yùn)行文件夾設(shè)置為當(dāng)前文件夾，默認(rèn)是Default文7結(jié)果分析京市京市模型應(yīng)用系子流域按年圖7-1結(jié)果分析界面圖系統(tǒng)設(shè)置11月11月二子謙域ET按月計(jì)算結(jié)果開始2000年1月輸出數(shù)據(jù)保存圖片圖7-2子流域按月計(jì)算結(jié)果段，點(diǎn)擊“更新圖表”,圖表將更新，如圖7-3所示為2002年8月至2005年5單位(mm單位(mm開始2002年8月結(jié)束2005年12月原新自0零四需器需mur寫四需輸出數(shù)據(jù)當(dāng)前子流城?顯示降水添加數(shù)據(jù)保存圖片圍圖7-4當(dāng)前子流域下拉列表框01年開始2000年1月結(jié)束[2005年12月更新圖表01年開始2000年1月結(jié)束[2005年12月更新圖表0(讀取數(shù)據(jù)?顯示降水=保存圖片圖表右上角下拉列表框?yàn)橹饕M水文要素(圖7-6),可通過“添加數(shù)據(jù)”按鈕添加到圖表中，如圖7-7所示。降融降融水滲滲地下水回歸圖7-6主要模擬水文要素圖7-7添加土壤含水量和地表徑流計(jì)算結(jié)果為便于比較，單擊計(jì)算結(jié)果柱狀圖，圖表以折線圖顯示，如圖7-8所示。雙05年謂05年謂0器器需器器器輸出數(shù)據(jù)01年身11月降水(mm)01年身11月降水(mm)11月05年明更新圖表0四票?顯示降0圖7-8子流域計(jì)算結(jié)果折線圖測(cè)數(shù)據(jù),雙擊后返回計(jì)算結(jié)果表，如圖7-9所示，以便觀測(cè)和模擬數(shù)據(jù)進(jìn)行對(duì)比。雙擊柱狀圖變?yōu)檎劬€圖，如圖7-10所示。單位(mm)10月12月單位(mm)10月12月0四0圖7-9模擬ET與遙感ET對(duì)比單位(mm)單位(mm)00r圖7-10模擬ET與遙感ET對(duì)比(折線圖)8保存類型T)ET-土壤水-地表徑流.esvET-土壤水-地表徑流取消保存(四)圖7-11輸出結(jié)果文本格式和圖片格式保存文件)編輯@)格式@)查看V)幫助0022000年2月，33.268,73.25,4.2212000年4月,12.9,64.655,.0092000年6月,12.882,58.164,.2832000年8月,20.026,59.84,2.3652000年9月，21.366,50.267,1.1462000年10月，37.898,32.761,1.912文件便)編輯@)格式(Q)查看V幫助QH)32000年1謂，62.54,62.898,19.5562000年2月,33.268,73.25,4.22000年3月,29.952,77.001,4.12000年4月,12.9,64.655,.02000年5月,5.669,59.352000年6月,12.882,58.164,.2832000年7月,6.816,51.911,02000年8月,20.026,59.84,2.3652000年9月，21.366,50.267,1.1462000年10月，37.898,32.761,1.912圖7-12讀取、輸出數(shù)據(jù)文本格式00年1日2000年1月開始]2000年1月結(jié)束2005年12月a譬年當(dāng)前子流域?顯示降水要3輸出數(shù)據(jù)營圖7-14不變的縱坐標(biāo)范圍輸出數(shù)據(jù)保存圖片結(jié)束2005年F顯示降水001年11月11月03年01年11月11月03年各河段流量計(jì)算結(jié)果添加數(shù)據(jù)讀取數(shù)據(jù)F圖7-16河段按月計(jì)算結(jié)果自各河段按年計(jì)算結(jié)果添加數(shù)據(jù)流出7.5HRU按月單擊“HRU按月”按鈕，彈出“水文響應(yīng)單元按月計(jì)算結(jié)果”對(duì)話框，如圖5-18所示，此對(duì)話框操作基本同7-1子流域按月。,才會(huì)輸出所有需要注意的是，在運(yùn)行SWAT過程中，選中,才會(huì)輸出所有HRU相關(guān)的信息，默認(rèn)情況下只輸出第一個(gè)HRU的相關(guān)信息，因一個(gè)HRU一個(gè)HRU的相關(guān)模擬數(shù)據(jù)，選其它HRU,彈出_對(duì)話框。HSTAT_CMxHRU輸出已被限制。定隆水(mm單位(mm)11月隆水(mm單位(mm)11月11月水文響應(yīng)單元按月計(jì)算結(jié)果開始2000年1月P顯示降水更新圖表0輸出數(shù)據(jù)添加數(shù)據(jù)保存圖片生地利用子流域主00圖7-18水文響應(yīng)單元按月計(jì)算結(jié)果地下水補(bǔ)給地下水補(bǔ)給承壓水儲(chǔ)量地表徑流圖7-19水文響應(yīng)單元模擬水文要素進(jìn)行率定。點(diǎn)擊“參數(shù)率定”,彈出參數(shù)率定對(duì)話框，如圖8-1所示。參數(shù)率定參數(shù)率定圖8-1參數(shù)率定8.1待率定參數(shù)根據(jù)經(jīng)驗(yàn)，結(jié)合模擬區(qū)域的實(shí)際情況，對(duì)SWAT模型模率定參數(shù)選擇率定參數(shù)選擇基流回歸常王壤水植被可和用量王壤飽和水力傳導(dǎo)度量大葉冠指數(shù)圖8-2率定參數(shù)選擇當(dāng)前值如圖8-4所示，通過交叉選取得到不同HRU,在不同級(jí)(流域級(jí)、子流域參數(shù)率定參數(shù)率定控制地表徑流取值范圍：0-100當(dāng)前值取值8.4.1徑流曲線數(shù)(CN2)硬質(zhì)裸地CN值較大，在80左右，松軟耕地和草地CN值較低，在30-60圖8-5CN值與降雨-徑流關(guān)系將CN值調(diào)為原理的1.2倍，模擬后結(jié)果見圖8-6,紅色保存圖片圖8-6相同降水不同CN值徑流值越小，就能從下層土壤補(bǔ)充更多的水分用于蒸發(fā)，蒸發(fā)會(huì)有所增加，如圖8-7所示，綠色部分“第一次模擬ET”采用的ESCO為0.1。8.4.4植物攝取補(bǔ)償系數(shù)(EPCO)植物攝取補(bǔ)償系數(shù)(EPCO),控制植物的水分?jǐn)z取，取值8.4.5土壤水植被可利用量(SOL_AWC)土壤水植被可利用量(SOL_AWC),是土壤水含量的指標(biāo)性參數(shù)，取值范土壤各層的植被可利用量不同，不同HRU,不同土層采用不同的值，賦值全部參數(shù)率定控制土壤水含量取值范圍：0-當(dāng)前值p15,0.14,0.08,0圖8-8土壤水植被可利用量賦值9情景分析模塊型應(yīng)用型應(yīng)用降水氣溫圖9-1系統(tǒng)菜單界面農(nóng)業(yè)管理北京流域分區(qū)數(shù)據(jù)輸入模擬系統(tǒng)設(shè)置降水?dāng)?shù)據(jù)[%百分率120氣溫?cái)?shù)據(jù)度)年份變化值確定取消降水?dāng)?shù)據(jù)(2年份變化值全部取消圖9-2降水氣溫模擬設(shè)置界面則表示當(dāng)年降水變?yōu)樵瓉淼?20%,若修改氣溫“變化值”文本框?yàn)?.1則表示支流域2支流域3支流域4支流域5支流域6支流域7支流域8支流域9支流域10支流域11n33363若在實(shí)際應(yīng)用中需要對(duì)土地利用類型進(jìn)行調(diào)整，如對(duì)子流域1進(jìn)行調(diào)整，則 相關(guān)操作如圖所示。(注：需要重新運(yùn)行模型才可得到相關(guān)結(jié)果)文件名面積(Km2)高程[m]雨量站支流域163346支流域11n土地利用1調(diào)整類型類型面積比例調(diào)整類整比例調(diào)整比例圖9-4地塊調(diào)整示例圖文件名面積(Km2)高程(m)支流域1支流域2支流域3支流域563支流域7支流域8支流域9支流域11去000080000.gub1.511739.921832結(jié)果調(diào)整調(diào)整后土地利用類型土地利用1土地利用2圖9-5地塊調(diào)整結(jié)果圖農(nóng)管理設(shè)置農(nóng)管理設(shè)置循環(huán)周期C單擊顯示類型05/09玉米8單擊顯示類型是單擊顯示類型月日壤類型O無灌溉C河道水C土地利用全部單擊顯示類型是否殘留3465圖9-6農(nóng)業(yè)調(diào)整示例圖一數(shù)據(jù)列，然后手動(dòng)敲入作物的名稱，但這樣做的后果往往是數(shù)據(jù)的不正確，造成無法與后臺(tái)數(shù)據(jù)庫匹配，因此，建議用戶單擊“播種”這兩個(gè)字，然后就會(huì)彈出如上圖所示的下拉列表框，在其中選擇作物類型即可，系統(tǒng)會(huì)在數(shù)據(jù)表中自耕作類型參數(shù)圖9-7耕作類型設(shè)置圖肥料名稱氤肥23456氨氮比例8細(xì)菌分配系數(shù)7)殺蟲劑模塊：本模塊功能類似于上述的施肥模塊兒，但無自動(dòng)殺蟲的功肥料類型數(shù)目：54有機(jī)氮比例00000寫入數(shù)據(jù)庫09附錄A不同計(jì)算結(jié)果的對(duì)比方法0開頭是不留空，以行為單位。2一行數(shù)據(jù)不垮行，無空行。4列內(nèi)容如存在半角逗號(hào)(即，)則用半角引號(hào)(即"")將該字段值包5列內(nèi)容如存在半角引號(hào)(即")則應(yīng)替換成半角雙引號(hào)("")轉(zhuǎn)義。6文件讀寫時(shí)引號(hào)，逗號(hào)操作規(guī)則互逆。文件P)編輯Q)格式Q)查看V幫助Q0高程，26.82,21.94,20.11,20.42,25.60,23.16,17.67,17.06,19.50,24.99,18.592000-01-01,000.0,000.0,000.0,000.0,000.0,000.0,000.0.000.0,000.0,2000-01-02,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-01-03,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,2000-01-04,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-01-05,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,2000-01-06,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-01-07,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-01-08,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,2000-01-09,000.0,000.0,000.0,000.8,000.0,000.0,000.0,000.0,000.02000-01-10,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,2000-01-11,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-01-12,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,2000-01-13,000.0,009.0,000.0,000.8,000.0,000.0,000.0,000.0,090.02000-01-14,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,2000-01-15,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-01-16,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-01-17,000.0,000.0,000.0,000.8,090.0,000.0,000.0,000.0,000.0,2000-01-18,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-01-19,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,2000-01-20,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-01-21,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,2000-01-22,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-01-23,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,2000-01-24,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-01-25,000.0.000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0.2000-01-26,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-01-27,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-01-28,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,2000-01-29,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-01-30,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,2000-01-31,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0.000.02000-02-01,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,2000-82-02,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-02-03,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,2000-02-04,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.02000-02-05,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.01H2000-2可以在Excel中方便的對(duì)其進(jìn)行編輯和存儲(chǔ)，而不需要關(guān)心SWAT原來的數(shù)據(jù)結(jié)構(gòu)。只需要遵守逗號(hào)隔開這一原則即可。2)關(guān)于情景對(duì)比的數(shù)據(jù)輸出和使用由于系統(tǒng)沒有提供在同一運(yùn)行窗口中比較兩種不同情景的功能，因此需要用戶首先利用結(jié)果分析模塊中的輸出數(shù)據(jù)功能將數(shù)據(jù)輸出，具體操作如下圖所示。在本次示例中共顯示了三組數(shù)據(jù)即：ET、潛在ET、和三組數(shù)據(jù)。點(diǎn)擊“保存數(shù)據(jù)”輸出之后，生成的CSV文件中依次包含了這當(dāng)重新模擬情景之后需要將其中一個(gè)數(shù)據(jù)加入，但是表中有3列數(shù)據(jù)，系統(tǒng)依據(jù)CSV文件第一列的數(shù)字進(jìn)行選擇，如圖4中所示，當(dāng)前值為“2”也就是說要對(duì)第二列數(shù)據(jù)進(jìn)行顯示即“ET”數(shù)據(jù)將被加入對(duì)比結(jié)果當(dāng)中。如圖5所示，在新的運(yùn)行情景中，點(diǎn)擊“讀取數(shù)據(jù)”按鈕，將已輸出的文件加入則系統(tǒng)自動(dòng)默認(rèn)加入“ET”數(shù)據(jù)進(jìn)行對(duì)比。用戶可以在Excel中手動(dòng)修改相關(guān)數(shù)值，以選擇對(duì)比的數(shù)據(jù)列。單位(mm)降水(mm)11月結(jié)束2005年12月單位(mm)降水(mm)11月結(jié)束2005年12月0輸出數(shù)據(jù)讀取數(shù)據(jù)?顯示降水0皿四器u器器圖3數(shù)據(jù)對(duì)比顯示窗口圖4導(dǎo)出數(shù)據(jù)對(duì)話窗口回2000年2000年2000年6潛潛在土51.4944.337圖5數(shù)據(jù)導(dǎo)出后顯示方式圖圖6數(shù)據(jù)導(dǎo)入后結(jié)果圖附錄B雨量站降水?dāng)?shù)據(jù)和統(tǒng)計(jì)年鑒降水?dāng)?shù)據(jù)的統(tǒng)一方法年\月，1,2,3,4,5,6,7,8,9,12000,9.4,0,9.9,13.2,37.2,12.5,55.8,96.6,23.9,28.2,8.2001,11.4,1.5,,9.6,5.1,100.5,133.4,65.9,8.9,50.4,20.2,2002,,0,5.4,32.5,29.7,92.4,54.4,120.6,54.7,18.2,0,5.2003,8.7,0.6,31.9,33.9,20.7,43,41.4,64.3,92.7,71.6,35.3,2004,0,7.7,0.1,23.5,54.9,98.1,140.2,59.4,79.8,11.7,8.8,1.2005,0.9,9.4,0.6,18.8,56.7,144.4,58.1,184.7,18.9,0.6,0,新的日降水?dāng)?shù)據(jù)文件原日降水?dāng)?shù)據(jù)csv文件原月降水?dāng)?shù)據(jù)csv文件降水?dāng)?shù)據(jù)校正系數(shù)csv文件本系統(tǒng)中采用的輸入數(shù)據(jù)主要是氣溫?cái)?shù)據(jù)和降水?dāng)?shù)據(jù)，在SWAT的運(yùn)行過程中，需要輸入日單位的降水和氣溫?cái)?shù)據(jù)這兩個(gè)數(shù)據(jù)的正確輸入與否，直接影響到模擬的結(jié)果。由于原有的SWAT模型中要求采用dbf格式輸入相關(guān)數(shù)據(jù)，其格式難度大，操作復(fù)雜，易出錯(cuò)且不易發(fā)現(xiàn)，因此根據(jù)SWAT實(shí)際應(yīng)用情況SWAT_CN采用新的數(shù)據(jù)讀入方式，以方便用戶快速準(zhǔn)確地輸入降水及氣溫相關(guān)數(shù)據(jù)。下面對(duì)兩種數(shù)據(jù)的輸入文件結(jié)構(gòu)進(jìn)行說明，二者結(jié)構(gòu)類似，以降水?dāng)?shù)據(jù)為例。降水?dāng)?shù)據(jù)文件共包含三個(gè)部分：分別為文件信息部分，測(cè)站信息部分和數(shù)據(jù)信息部分。圖1表結(jié)構(gòu)概覽(1)文件信息：文件信息是指該輸入文件的基本信息包括：數(shù)據(jù)的起始時(shí)間和終止時(shí)間，測(cè)站的個(gè)數(shù)，信息來源及說明。數(shù)據(jù)的起止時(shí)間代表下面第三部分輸入的時(shí)間序列長度，測(cè)站的個(gè)數(shù)則表示數(shù)據(jù)輸入的列數(shù)，文件信息是指數(shù)據(jù)的提供單位，而說明是指對(duì)文件的必要說明。ABCDEF34經(jīng)度20.42入一個(gè)數(shù)據(jù)即日降水量數(shù)據(jù)，但是對(duì)于一個(gè)氣溫00000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000第三列輸入測(cè)站個(gè)數(shù)，以阿拉伯字母輸入，第ABCDEF使用說明第1行：第一列輸入“經(jīng)度”二字第2行：第一列輸入“緯度”二字第3行：第一列輸入“高程”二字。相關(guān)測(cè)站的數(shù)據(jù)是以北京當(dāng)?shù)刈鴺?biāo)為單位進(jìn)行換算的，其值以米為單位 (m)。依次輸入完所有測(cè)站信息即可(第2個(gè)測(cè)站的信息在第3列，第3個(gè)測(cè)站的信息在第4列，以此類推),測(cè)站輸入的數(shù)目要與第一部分文件234始就要輸入實(shí)測(cè)的數(shù)據(jù)，從第5行開始每1行的第1列都要輸入日期，日列輸入第1個(gè)測(cè)站當(dāng)日的數(shù)據(jù)，第3列輸入第2個(gè)測(cè)站當(dāng)日的數(shù)據(jù)，以此類推，從中我們可以看出每一行表示當(dāng)日18個(gè)000000000000000第二個(gè)測(cè)站當(dāng)日的數(shù)據(jù)00000000000000000000000000000000000000000000000000000000000000000000000圖7第三部分輸入說明2列開始，每隔1列為一個(gè)測(cè)站的數(shù)據(jù)，即一個(gè)測(cè)站包括兩列數(shù)據(jù)，若有3個(gè)測(cè)站，則應(yīng)該包括6列數(shù)據(jù)，每個(gè)測(cè)站的這兩列數(shù)據(jù)中第1列是最高氣的順序順次輸入。A2緯度B3.6 2.7-4.42圖8氣象數(shù)據(jù)輸入說明示意圖在Excel中編輯好的數(shù)據(jù)，點(diǎn)擊“另存為”選擇存為“CSV(逗號(hào)分隔)”格式即可完成存儲(chǔ)即可。4)SWAT_CN系統(tǒng)中氣溫?cái)?shù)據(jù)及降水?dāng)?shù)據(jù)的編輯(記事本)在一般情況下使用Excel進(jìn)行編輯即方便又可滿足應(yīng)用需要，但是CSV文件的另一個(gè)特點(diǎn)是應(yīng)用的廣泛性，因?yàn)樗梢员挥浭卤敬蜷_，所以即使計(jì)算機(jī)上沒有安裝Excel也可進(jìn)行相關(guān)的編輯。高程，26.82,21.94,20.11,20.42,25.60,23.16,17.67,17.06,19.52000-01-03,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,02000-01-05,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,02000-01-08,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,02000-01-11,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,02000-01-13,000.0,000.0,000.0,000.0,000.0,000.0,000.0,000.0,0圖9在記事本中打開降水?dāng)?shù)據(jù)的CSV文件從上圖中可以看到其排列組織方式與Excel中打開的類似，沒有了×系統(tǒng)設(shè)置氣溫?cái)?shù)據(jù)。點(diǎn)擊“查找”按鈕，系統(tǒng)會(huì)彈出對(duì)話框，如圖2所示，用戶可13所示。查找選擇降水?dāng)?shù)據(jù)文件選擇降水?dāng)?shù)據(jù)文件查找范圍(Importpep1.csv桌面網(wǎng)上部居打開@)為方便用戶導(dǎo)入數(shù)據(jù)，系統(tǒng)提供了幾個(gè)方便的小工具，包括PCP->CSV,tmp->CSV這兩個(gè)按鈕的功能是將SWAT運(yùn)行時(shí)的降水?dāng)?shù)降水?dāng)?shù)據(jù)(參看附錄B)。查找退出導(dǎo)入6)SWAT_CN系統(tǒng)中氣溫?cái)?shù)據(jù)及降水?dāng)?shù)據(jù)空間位置顯示為了方便用戶校驗(yàn)雨量站和氣溫站的位置是否與實(shí)際位置相匹配，在打開系統(tǒng)的GIS模塊，點(diǎn)擊菜單中的“文件”選項(xiàng)，選擇“打開降水(氣溫)數(shù)據(jù)”,彈出如下對(duì)話框。選擇降水(氣溫)數(shù)據(jù)文件選擇降水(氣溫)數(shù)據(jù)文件查找范圍)打開@)如果用戶先前編輯的數(shù)據(jù)有效，系統(tǒng)會(huì)自動(dòng)根據(jù)CSV文件中第2、3、4行的數(shù)據(jù)顯示相關(guān)測(cè)站的空間位置。圖16顯示結(jié)果附錄D投影坐標(biāo)的轉(zhuǎn)換方法軸為極軸，所有通過地球南北極的平面，均稱例如：北京位于北緯40度和東經(jīng)116度的交點(diǎn)附近，昆明近似橢球體投影得到的，它對(duì)應(yīng)于某個(gè)地理坐標(biāo)系。投影坐標(biāo)系由以下參數(shù)確定：投影方法(比如高斯一克呂格、Lambert投影、Mercator投影)要對(duì)其設(shè)置統(tǒng)一的坐標(biāo)系統(tǒng)，對(duì)于雨量站、參數(shù)數(shù)值Standard_Parallel_1Standard_Parallel_20表2地理坐標(biāo)系參數(shù)參數(shù)值0D_Krasovsky_1940Krasovsky_19403)地理坐標(biāo)向投影坐標(biāo)轉(zhuǎn)換的方法(以Arcgis9.2為例)先設(shè)置圖層組的投影信息，后將該圖層組內(nèi)的各圖層輸出為具有正確信息的數(shù)據(jù)。行政界面如果在ArcGIS9.2中載入一個(gè)沒有投影坐標(biāo)，只有北京當(dāng)?shù)氐乩碜鴺?biāo)的數(shù)據(jù)，則ArcGIS9.2會(huì)提示，這個(gè)圖像沒有相關(guān)的投影信息，是否要繼圖層組的屬性信息b▼(注：不是圖層)2口X圖2選擇圖層組屬性“Properties”選項(xiàng)可彈出下圖中的對(duì)話框：AnmotationGroups|ExtGeneral|DataFrame[CoordinaCurentcoordnatesystem:UnknownClearJransformations,..Nodiy,中Layers沒有設(shè)置相圖3投影信息對(duì)話框在這個(gè)對(duì)話框中，顯示了該圖層組的屬性，選中標(biāo)簽“CoodinateSystem”可以看到該圖層組的投影信息為“unknown”,需要設(shè)置相關(guān)信AmnotationGroups|ExtentRectaGeneral|DataPraneCoordinateUnknownJransformatiors..Nodfy.選擇投影信息圖4創(chuàng)建新的投影坐標(biāo)點(diǎn)擊“new”按鈕，選擇ProjectedCoordinateSystem會(huì)彈出如下對(duì)話框，可以進(jìn)行具體的參數(shù)設(shè)置。IeyIeyProjectedCoordin輸入名稱選擇投影Name!0.000000000000000000選擇度量單位(默認(rèn)即可)perunit:選擇地理坐標(biāo)Select..New.輸入表2中的數(shù)據(jù)Unit-圖5投影參數(shù)的設(shè)置在上面的對(duì)話框中，首先輸入自己定義的投影坐標(biāo)的名稱，如，“swat_cn”,設(shè)置內(nèi)容包括三個(gè)部分，分別為：Projection,LinearUnit,二個(gè)選項(xiàng)),然后在下面的參數(shù)輸入部分中將表1的參數(shù)填入即可。LinearUnit:在本次應(yīng)用中，采用系統(tǒng)的默認(rèn)值即可。NorthAmerieaIndonesianNational.prttModfy...確定圖8完成設(shè)定此時(shí)地圖的顯示形式已經(jīng)跟隨投影發(fā)生變化，需要將設(shè)置好的數(shù)據(jù)進(jìn)行輸出，產(chǎn)生另一份信息完全一樣，但投影方式不同的數(shù)據(jù)。shapefile:F:5WAT_BJ|SWATBeij確定取消應(yīng)用(6此處為圖層的信息XBenovekAnalystToolsConvertSymlolozytoBap渡圖10數(shù)據(jù)輸出菜單選項(xiàng)氧0sforcaiEspotAlesussloriopploaitgougapothoafpohurodalasotmagoodyta2圖11輸出對(duì)話框圖12新的圖層投影信息NH3(x10)北京師范大學(xué)水科學(xué)研究院2009年7月1水量方面(河道自產(chǎn)流量、河段間取水量分配、沿河點(diǎn)源排污口廢水、污水量)和水質(zhì)方面(沿河點(diǎn)源與非點(diǎn)源污染負(fù)荷),空間上沿著河流縱向方向上的水量圖1系統(tǒng)啟動(dòng)界面2河流定義模塊圖2河流設(shè)置模塊該模塊可以設(shè)置河流編號(hào)(注：編號(hào)作為計(jì)算標(biāo)識(shí)符，應(yīng)具有唯一性),河23系統(tǒng)設(shè)置模塊設(shè)置河段編號(hào)及匯流關(guān)系設(shè)置河段編號(hào)及匯流關(guān)系444計(jì)算單元計(jì)算單元長度米圖3系統(tǒng)設(shè)置模塊3入口入口42679出口處(即河源處)向干流的出口處排序。遇到的支流依次尋找支流的河源，從河流4由于1號(hào)河(通惠河下段)及2號(hào)河(小中河)匯入3號(hào)河(北運(yùn)河北段)。若河段B為0則表示僅有A匯入，如第7行表示只有6號(hào)河(涼水河西段)匯入7號(hào)河(涼水河中段);若A和B相等則表示，無河流匯入，僅為自身河段進(jìn)行計(jì)4月7月3500調(diào)水水源名稱隔榆河河流號(hào)名稱北運(yùn)河中段涼水河西段河段號(hào)圖5入流設(shè)置模塊5單擊輸入相關(guān)值即可，若設(shè)置年均值則雙擊年均值文本框系統(tǒng)即可自動(dòng)將值平均至各月，若有詳細(xì)資料也可對(duì)月均值進(jìn)行單獨(dú)設(shè)置。單擊“保存”按鈕即可完成數(shù)據(jù)設(shè)置。取水設(shè)置取水設(shè)置1月2月12月保存D回流號(hào)名稱圖6取水設(shè)置模塊66點(diǎn)源設(shè)置河遠(yuǎn)號(hào)名稱河段名稱周3點(diǎn)源名稱點(diǎn)進(jìn)112月5圖7點(diǎn)源設(shè)置07月9月非點(diǎn)源名稱果圖保存河段名稱河段號(hào)河段長度回流號(hào)名稱選擇河段000000078結(jié)果分析模塊0圖9結(jié)果分析模塊1)一級(jí)衰變動(dòng)力學(xué)方程(i)流向下一個(gè)(i+1)計(jì)算單元前的濃度，即第i個(gè)計(jì)算單元的斷尾濃度，第i+1個(gè)計(jì)算單元的段首濃度(第i+1個(gè)計(jì)算單元沒有排污口的情況下)按照下面式(3-2)中，k——污染物綜合自凈系數(shù)(1/d);u——平均流速(km/d);82)混合方程Q.——累計(jì)流量(m3/s);P??——排污口污染物負(fù)荷(t/a);3)混合方程J;=LQL+——污染物濃度(mg/L);Q?——累計(jì)流量(m3/s);JA——支流污染物負(fù)荷(g/s)。n表示匯入該干流的支流終點(diǎn)計(jì)算單元位置，對(duì)于不同長度和單元長度的河9操作說明單一河段查看：若需要查看單個(gè)河段值的變化，在添加河段按鈕的右側(cè)的下拉框中選擇相應(yīng)河段即可，所選河段的流量、COD、及氨氮的變化情況即會(huì)顯示在圖中。(注：為使坐標(biāo)軸顯示范圍盡量一致，氨氮的值在顯示過程中是乘以10得到的，實(shí)際應(yīng)用中在背后數(shù)據(jù)庫中所存儲(chǔ)值為真實(shí)值)用戶可以通過“數(shù)據(jù)輸出”按鈕將其輸出為可供Excel和記事本打開的CSV文件格式，同時(shí)也可以將當(dāng)前圖片進(jìn)行輸出以便直觀顯示結(jié)果。匯流關(guān)系查看：在實(shí)際模擬過程中，需要查看具有復(fù)雜匯流關(guān)系的全流域的模擬結(jié)果，在本模塊中，若按照匯流河段的順序添加相關(guān)河段，即可得到相關(guān)河段的水質(zhì)變化圖。如圖9所示即可小中河進(jìn)入北運(yùn)河北段后其水質(zhì)的變化。門門施加河北運(yùn)河北段j50圖10具有匯流關(guān)系的多河段結(jié)果顯示監(jiān)測(cè)數(shù)據(jù)的添加；若用戶有河段某一具體位置的監(jiān)測(cè)數(shù)據(jù)，則可以通過“讀取數(shù)據(jù)”按鈕進(jìn)行添加，點(diǎn)擊讀取數(shù)據(jù)按鈕導(dǎo)入相關(guān)的CSV文件，就可以在相應(yīng)位置上添加河段監(jiān)測(cè)點(diǎn)的信息，如下圖所示：②圖11數(shù)據(jù)讀入窗口北運(yùn)河北段北運(yùn)河中段北運(yùn)河下段流量CODNH3(x10)0距離(km)【北運(yùn)河北段~北運(yùn)河中段>北運(yùn)河下段】北運(yùn)河下段全年讀取數(shù)據(jù)輸出數(shù)據(jù)保存圖片圖12監(jiān)測(cè)數(shù)據(jù)顯示所采用的CSV格式，在示例數(shù)據(jù)中，有顯示，第一個(gè)字段為監(jiān)測(cè)點(diǎn)距離河段上游的距離，接下來所填寫數(shù)據(jù)依次為COD氨氮的實(shí)際監(jiān)測(cè)值，以CSV格式進(jìn)行保存即可，數(shù)據(jù)輸入應(yīng)從第二行開始。北運(yùn)河2004年水質(zhì)數(shù)據(jù)C 1 1 1 2 42.4容量計(jì)算功能 52.5河段水質(zhì)預(yù)測(cè)功能 7 9WaterSoftWECToolVers1WECToolVer1.0是為河段水環(huán)境容量(水體納污能力)計(jì)算而設(shè)計(jì)的，包括河段數(shù)據(jù)瀏覽、排污口數(shù)據(jù)瀏覽、河段水環(huán)境容量計(jì)算、河段水質(zhì)預(yù)測(cè)、水環(huán)境容量標(biāo)準(zhǔn)查詢修改等功能，可滿足河段污染物排放削減目標(biāo)計(jì)算的需要。系統(tǒng)界面如圖1所示。河段屬性-梨園鎮(zhèn)東、張家灣鎮(zhèn)東、河段添加河段田除數(shù)據(jù)導(dǎo)出水質(zhì)分析圣廟排口東方排口上店污水口氨氮降解系數(shù)起始斷面終止斷面99比運(yùn)河地北運(yùn)河比325mg/L氮顏現(xiàn)狀濃度081/日COD容量氨氮容量立方米散氮氮總負(fù)荷578.19噸容量計(jì)水環(huán)境容量計(jì)算模塊提供河段及其相關(guān)排污口的數(shù)據(jù)瀏覽的功能，使用者可以通過系統(tǒng)上方標(biāo)有“前”“后”的圖標(biāo)，對(duì)記錄進(jìn)行移動(dòng)，依次瀏覽各河段信息及每個(gè)河段所包含的排污口的基本信息，由于河段信息和排污口信息具有對(duì)應(yīng)關(guān)系，因此二者的顯示方式采用聯(lián)動(dòng)顯示的方式，即每顯示一個(gè)不同的河段，就顯示所在河段排污口的信息。WaterSoftWECToolVers2排污口排污口數(shù)日CDD總負(fù)荷7406.34噸氨氮總負(fù)荷578.19噸添排污口起始斷面終止斷面圖2河段基本信息瀏覽圖3排污口數(shù)據(jù)瀏覽2.2數(shù)據(jù)編輯功能算過程中需要對(duì)已有信息進(jìn)行修改，添加或刪除，模塊可以有效的幫助使用者完成上述圖4數(shù)據(jù)操作相關(guān)按鈕圖河段長度設(shè)計(jì)流速COD降解系數(shù)氨氮降解系數(shù)3立方米渺米終止斷面所屬行政區(qū)潮白河圖5數(shù)據(jù)填寫方式WaterSoftWECToolVers32)應(yīng)用數(shù)據(jù)添加3)取消數(shù)據(jù)添加4)河段數(shù)據(jù)刪除5)河段數(shù)據(jù)修改2.2.2排污口數(shù)據(jù)操作添加排污口添加排污口1)數(shù)據(jù)添加2)刪除WaterSoftWECToolVersion43)修改需要對(duì)排污口數(shù)據(jù)進(jìn)行修改時(shí)，點(diǎn)擊窗體下方“編輯數(shù)據(jù)”的按鈕，用鼠標(biāo)左鍵點(diǎn)擊排污口記錄所在行，然后點(diǎn)輸入相關(guān)數(shù)據(jù)即可完成操作。4)數(shù)據(jù)導(dǎo)出功能利用該功能將計(jì)算結(jié)果導(dǎo)出為可供Excel(csv格式)或文本文件瀏覽程序打開的文件，只需要指定文件即可。導(dǎo)出之后的數(shù)據(jù)格式如圖7所示。ABCBEFGMJK2 3C0201北運(yùn)河北關(guān)間4北運(yùn)河下端C0202北運(yùn)河12VMd205C0300運(yùn)湖減河北關(guān)分洪閘C0400港溝河許各莊6708通惠河溫榆河新河閘9000801原水河馬橋村西張家清閘小中河下端C0902故道間北關(guān)分洪閘VV 0.737G5中中壩河中壩河安里玉帶河西門篩子莊北關(guān)閘于府間軍屯間小屯 VVVVV0.3140.4060.036U00861MN0P0RSUY置0299.389 29.838梨園鎮(zhèn)東、張7406.3451.984部縣鎮(zhèn)北西47.60.916潮縣鎮(zhèn)、永樂486.9202 025.543臺(tái)湖鎮(zhèn)南、馬 200 221.294馬駒橋鎮(zhèn)南2.016潮縣鎮(zhèn)東部000.28永順鎮(zhèn)、梨園022圖7導(dǎo)出數(shù)據(jù)的格式5WaterSoftWECToolVers水質(zhì)標(biāo)準(zhǔn)水質(zhì)標(biāo)準(zhǔn)1類IV類區(qū)氨氨圖9環(huán)境容量計(jì)算標(biāo)準(zhǔn)設(shè)置模塊環(huán)境容量圖10規(guī)劃水質(zhì)目標(biāo)的設(shè)定系統(tǒng)將自動(dòng)調(diào)用相應(yīng)的水質(zhì)標(biāo)準(zhǔn)進(jìn)行水環(huán)境容量計(jì)算。62.5河段水質(zhì)預(yù)測(cè)功能圖11河段水質(zhì)預(yù)測(cè)功能污口污染排放統(tǒng)計(jì)和降解系數(shù)選擇的合理性。WaterSoftWECToo73.附錄計(jì)算模型及其參數(shù)的選擇3.1水環(huán)境容量計(jì)算模型水環(huán)境容量(河段納污能力)計(jì)算采用基于一維對(duì)流推移方程的段尾控制模型，假設(shè)集中點(diǎn)源的自凈長M——納污能力(kg/d);C?——下段水質(zhì)標(biāo)準(zhǔn)(mg/L);Q,——設(shè)計(jì)流量(m3/s);k——污染物綜合自凈系數(shù)(1/d);L——河流間距(km);u——平均流速(km/d)。3.2模型參數(shù)的設(shè)定根據(jù)所選定的模型，對(duì)兩種重要水環(huán)境指標(biāo)COD和NH?-N做納污能力計(jì)算，計(jì)算河流的納污能力需要知(1)水質(zhì)標(biāo)準(zhǔn)C?、河流間距L的確定根據(jù)我國《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)基本項(xiàng)目標(biāo)準(zhǔn)限值》(GB3838-2002),本計(jì)算中涉及的水質(zhì)標(biāo)準(zhǔn)C數(shù)值如下表7-3所示。河流間距L由劃分河段時(shí)根據(jù)統(tǒng)計(jì)資料獲得。表1水質(zhì)標(biāo)準(zhǔn)Cs、Co數(shù)值IV類水質(zhì)標(biāo)準(zhǔn)Cs、Co(2)水質(zhì)降解系數(shù)k確定污能力核定中推薦的降解系數(shù)值(如下表2)。WaterSoftWECToolVersi8表2一般河道水質(zhì)降解系數(shù)參考值表水質(zhì)降解系數(shù)參考值(1/d)優(yōu)(相應(yīng)水質(zhì)為Ⅱ~Ⅲ類)中(相應(yīng)水質(zhì)為Ⅲ~IV類)劣(相應(yīng)水質(zhì)為V類或劣V類)(3)設(shè)計(jì)流量Q的確定(4)設(shè)計(jì)流速u確定WaterSoftWECToolVersi91)水質(zhì)降解系數(shù)k確定2)設(shè)計(jì)流量Q的確定分段起止點(diǎn)自產(chǎn)流量北運(yùn)河(2段)北關(guān)閘一榆林閘通惠河(1段)馬橋村西一新河閘涼水河(3段)新河閘一張家灣閘張家灣閘一榆林莊房辛店一于府閘鳳港減河(3段)于府閘一軍屯閘蕭太后河(1段)口子村一張家灣涼水河 3)設(shè)計(jì)流速u的確定分段起止點(diǎn)設(shè)計(jì)流量(m3/s)設(shè)計(jì)流速(m/s)北運(yùn)河(2段)北關(guān)閘一榆林閘通惠河(1段)馬橋村西一新河閘涼水河(3段)新河閘一張家灣閘張家灣閘一榆林莊，入北運(yùn)河房辛店一于府閘鳳港減河(3段)于府閘一軍屯閘蕭太后河(1段)口子村一張家灣，入涼水河 步驟使用數(shù)據(jù)使用軟件子流域數(shù)據(jù)準(zhǔn)備DEM、河流GIS圖層ArcGIS/GoogleEarth土壤、土地?cái)?shù)據(jù)準(zhǔn)備土壤、土地利用屬性數(shù)據(jù)Excel計(jì)算表/1stOpt/SPAW氣象數(shù)據(jù)的準(zhǔn)備氣象數(shù)據(jù)SwatWeather水文響應(yīng)單元的劃分土壤、土地利用GIS圖層ArcSWATArcSWAT初步運(yùn)行ArcSWATArcSWAT到BNU-SWAT的移植MapInfo步驟使用數(shù)據(jù)使用軟件BNU-SWAT的率定徑流監(jiān)測(cè)數(shù)據(jù)BNU-SWAT/SWAT-CUPBNU-SWAT結(jié)果分析BNU-SWATBNU-SWAT氣候模擬BNU-SWATBNU-SWAT土地利用及作物變化模擬BNU-SWATBNU-SWAT地下水模擬BNU-SWAT事COD、氨氮濃度模擬WQS、WEC1BNU-SWAT&SWAT-CUP培訓(xùn)教材第二講培訓(xùn)相關(guān)軟件1.培訓(xùn)相關(guān)軟件1)培訓(xùn)用計(jì)算機(jī)中的軟件軟件名稱用途河流數(shù)字化土壤粒徑轉(zhuǎn)換計(jì)算Excel計(jì)算表生壤參數(shù)計(jì)算工具氣象數(shù)據(jù)提取及氣象發(fā)生器參數(shù)生成工具GIS數(shù)據(jù)shp到tab格式的轉(zhuǎn)換水環(huán)境容量計(jì)算工具wQsV1.0SWAT模型的參數(shù)率定工具2)培訓(xùn)光盤中的軟件及數(shù)據(jù)(1)GIS軟件1BNU-SWAT&SWAT-CUP培訓(xùn)教材方正綜藝簡體BNU-SGDEM數(shù)據(jù)土地利用分布圖水系圖北京市通州區(qū)BNU-SWAT樣例2.ArcSWATV2.1的安裝進(jìn)入ArcSWAT文件夾執(zhí)行安裝文件SWATInstall.msi,按提示完成安裝。安裝之前，系統(tǒng)必須裝有ArcGIS9.2(SP4)。進(jìn)入GoogleEarth文件夾執(zhí)行安裝文件，按提示完成安裝。4.BNU-SWATV1.0的安裝進(jìn)入培訓(xùn)光盤的MapXTrial文件夾，雙擊文件SETUPEXE,單擊【Next】按紐進(jìn)入2BNU-SWAT&SWAT-CUP培訓(xùn)教材這一步將顯示安裝可選項(xiàng)，此處必須將所有選項(xiàng)全部安裝，如下圖所示，單擊【Next】按鈕繼續(xù)；Selectthecomponentyouwanttoinstal,clearthe按默認(rèn)選項(xiàng)，完成其余安裝步驟。(2)字體的安裝進(jìn)入【控制面板】-【字體】菜單，安裝培訓(xùn)光盤Font文件夾中的字體。(3)BNU-SWAT的安裝進(jìn)入培訓(xùn)光盤的BNU-SWAT文件夾，雙擊SETUPEXE文件，按照提示安裝軟件。安裝BNU-SWAT時(shí)自動(dòng)安裝。安裝BNU-SWAT時(shí)自動(dòng)安裝。5.SWAT-CUPV2.15(Beta)的安裝進(jìn)入培訓(xùn)光盤的SWAT-CUP文件夾，雙擊SwatCupSetup2.1.5.exe文件，按照提示安裝北京師范大學(xué)水科學(xué)研究院數(shù)字流域?qū)嶒?yàn)室3BNU-SWAT&SWAT-CUP培訓(xùn)教材第三講DEM數(shù)據(jù)的利用國際科學(xué)數(shù)據(jù)服務(wù)平臺(tái)InternationalScient軟軟件/文檔用戶空間關(guān)于我們圖1中科院國際數(shù)據(jù)服務(wù)平臺(tái)1)SRTM數(shù)據(jù)的獲取SRTM的數(shù)據(jù)組織方式為：每5度經(jīng)緯度方格劃分一個(gè)文件，共分為24行(-60至60度)和72列(-180至180度)。文件命名規(guī)則為srtm_XX_YY.zip,XX表示列數(shù)(01-72),圓BNU-SWAT&SWAT-CUP培訓(xùn)教材數(shù)據(jù)的空間范圍。國際科學(xué)數(shù)據(jù)服務(wù)平臺(tái)InternationalScientincDataServ網(wǎng)站首頁數(shù)據(jù)目錄數(shù)據(jù)檢索數(shù)據(jù)服務(wù)共身成軟件/文檔用戶空間關(guān)于我們提文后調(diào)圖2SRTM數(shù)據(jù)查詢界面在輸入研究區(qū)域范圍后點(diǎn)擊右下角的提交按鈕后即進(jìn)入數(shù)據(jù)下載界面，如圖3所示。選擇屏幕右側(cè)SRTMV4.1(tiff)文件進(jìn)行下載。將下載的文件解壓縮后即為ArcGIS9.2可打開的DEM數(shù)據(jù)。國際科學(xué)數(shù)據(jù)服務(wù)平臺(tái)enno16m注料字運(yùn)距InternationalScientincDataServicePl網(wǎng)站首頁數(shù)需目錄數(shù)需檢索數(shù)據(jù)服務(wù)共享成果軟件/文檔用戶空間關(guān)于我們神1◆刺I15口數(shù)字高程(0mg)坡度(mg)境向(mg)坡位0mg)圖3SRTM數(shù)據(jù)下載界面北京師范大學(xué)水科學(xué)研究院數(shù)字流域?qū)嶒?yàn)室2BNU-SWAT&SWAT-CUP培訓(xùn)教材該數(shù)據(jù)集利用ASTERGDEM第一版本(V1)的數(shù)據(jù)進(jìn)行加工得來，是全球空間分辨率為30米的數(shù)字高程數(shù)據(jù)產(chǎn)品。數(shù)據(jù)以1度經(jīng)緯度網(wǎng)格為一個(gè)數(shù)據(jù)文件，采用UTM平面坐標(biāo)系統(tǒng)進(jìn)行存儲(chǔ)。數(shù)據(jù)查詢界面如圖4所示，需要輸入研究區(qū)域經(jīng)緯度確定數(shù)據(jù)下載范圍。點(diǎn)擊篩選按鈕后即獲得數(shù)據(jù)文件。點(diǎn)擊右側(cè)操作按鈕，選擇數(shù)字高程數(shù)據(jù)下當(dāng)研究區(qū)域的DEM數(shù)據(jù)存在于多個(gè)數(shù)據(jù)文件時(shí)(如圖5),需要對(duì)多景DEM數(shù)據(jù)進(jìn)圖5拼接前的兩景數(shù)據(jù)以下部分以拼接left_piece和right_piece兩景DEM數(shù)據(jù)為例。在ArcMap中將兩個(gè)文件打開。使用ArcToolbox中DataManagementTools中Raster工具的Mos1)InputRater:選擇要拼接的數(shù)據(jù)，即left_piece和right_piece兩個(gè)文件3)Raterdatasetnamewithextension:確定生成文件的名稱4)Pixel_