SWAT水文模型之歐陽科創(chuàng)編

歐陽科創(chuàng)編 2021.02.05SWAT水文模型介紹時間：2021.02.05創(chuàng)作：歐陽科1概述SWAT（SoilandWaterAssessment丁001）模型是美國農(nóng)業(yè)部（USDA）農(nóng)業(yè)研究局（ARS）開發(fā)的基于流域尺度的一個長時段的分布式流域水文模型。它主要基于SWRRB模型，并吸取了CREAMS、GLEAMS、EPIC和ROTO的主要特征。SWAT具有很強的物理基礎，能夠利用GIS和RS提供的空間數(shù)據(jù)信息模擬地表水和地下水的水量和水質(zhì)，用來協(xié)助水資源管理，即預測和評估流域內(nèi)水、泥沙和農(nóng)業(yè)化學品管理所產(chǎn)生的影響。該模型主要用于長期預測，對單一洪水事件的演算能力不強，模型主要由 8個部分組成：水文、氣象、泥沙、土壤溫度、作物生長、營養(yǎng)物、農(nóng)業(yè)管理和殺蟲劑。SWAT模型擁有參數(shù)自動率定模塊，其采用的是Q.Y.Duan等在1992年提出的SCE-UA算法。模型采用模塊化編程，由各水文計算模塊實現(xiàn)各水文過程模擬功能，其源代碼公開，方便用戶對模型的改進和維護。2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.052模型原理SWAT模型在進行模擬時，首先根據(jù)DEM把流域劃分為一定數(shù)目的子流域，子流域劃分的大小可以根據(jù)定義形成河流所需要的最小集水區(qū)面積來調(diào)整，還可以通過增減子流域出口數(shù)量進行進一步調(diào)整。然后在每一個子流域內(nèi)再劃分為水文響應單元HRU。HRU是同一個子流域內(nèi)有著相同土地利用類型和土壤類型的區(qū)域。每一個水文響應單元內(nèi)的水平衡是基于降水、地表徑流、蒸散發(fā)、壤中流、滲透、地下水回流和河道運移損失來計算的。地表徑流估算一般采用SCS徑流曲線法。滲透模塊采用存儲演算方法，并結(jié)合裂隙流模型來預測通過每一個土壤層的流量，一旦水滲透到根區(qū)底層以下則成為地下水或產(chǎn)生回流。在土壤剖面中壤中流的計算與滲透同時進行。每一層土壤中的壤中流采用動力蓄水水庫來模擬。河道中流量演算采用變動存儲系數(shù)法或馬斯金根演算法。模型中提供了三種估算潛在蒸散發(fā)量的計算方法一Hargreaves、PriestleyTaylor和PenmanMonteith。每一個子流域內(nèi)侵蝕和泥沙量的估算采用改進的USLE方程，河道內(nèi)泥沙演算采用改進的Bagnold泥沙運移方程。植物吸收的氮采用供需方法計算，植物的氮日需求量是植物與生物量中氮濃度的函數(shù)。土壤中向植物供給氮，當需求超過供給時，出現(xiàn)營養(yǎng)物壓力。2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05地表徑流、壤中流和滲透過程運移的硝態(tài)氮量由水量和土壤層中的平均硝態(tài)氮濃度來估計。泥沙中運移的有機氮采用McElroyetal.開發(fā)的負荷方程，后經(jīng)進一步改進。該負荷方程基于土壤表層的有機氮濃度、泥沙量和富集率來估計徑流中的有機氮損失。植物吸引I的磷采用與氮相似的供需方法。徑流中帶走的可溶解磷采用土壤表層中的不穩(wěn)定磷、徑流量和磷土分離系數(shù)來計算。泥沙運移的磷采用與有機氮運移相同的方程。河道中營養(yǎng)物的動態(tài)模擬采用QUAL2E模型。模型中采用的水量平衡表達式為：式中：“”為土壤最終含水量，mm；"'為土壤前期含水量，mm；t為時間步長，d；，”.為第i天降水量，mm；為第i天的地表徑流，mm；e為第i天的蒸發(fā)量，mm；h「為第i天存在于土壤剖面底層的滲透量和側(cè)流量，mm；Q『為第i天地下水含量，mm。SWAT模型水文循環(huán)陸地階段主要有水文、天氣、沉積、土壤溫度、作物產(chǎn)量、營養(yǎng)物質(zhì)和農(nóng)業(yè)管理等部分組成。模型徑流量產(chǎn)生計算流程圖如下圖2-1所示。圖2-1SWAT模型產(chǎn)流計算流程圖2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.052.1.1地表徑流當落到地表的降水量多余入滲量時產(chǎn)生地表徑流。SWAT模型采用SCS徑流曲線法計算。SCS曲線方程自上世紀50年代逐漸得到廣泛使用，屬于經(jīng)驗模型，是對全美小流域降水與徑流關系20多年的研究成果。模型能反應不同土壤類型和土地利用方式及前期土壤含水量對降雨徑流的影響，它是基于流域的實際入滲量（F）與實際徑流量（Q）之比等于流域該場降雨前的最大可能入滲量（S）與最大可能徑流量（:［；）之比的假定基礎上建立的5?5模型的降雨徑流基本關系表達式如下：式中：假定潛在徑流量（：1）為降水量（P）與由徑流產(chǎn)生前植物截留、初滲和填洼蓄水構(gòu)成的流域初損（）的差值。由此推到上式有：初損受土地利用、耕作方式、灌溉條件、冠層截留、下滲、填洼等因素的影響，它與土壤最大可能入滲量S呈一定的正比關系，美國農(nóng)業(yè)部土壤保持局在分析了大量長期的實驗結(jié)果基礎上，提出了二者最合適的比例系數(shù)為0.2，即：由此可得SCS方程為：流域當時最大可能滯留量::在空間上與土地利用方式、土壤類型和坡度等下墊面因素密切相關，模型弓1入的；n值可歐陽科創(chuàng)編 2021.02.05歐陽科創(chuàng)編 2021.02.05較好地確定:，公式如下：CN是一個無量綱參數(shù)，CN值是反映降雨前期流域特征的一個綜合參數(shù)，它是前期土壤濕度、坡度、土地利用方式和土壤類型狀況等因素的綜合。2.1.2蒸散發(fā)模型考慮的蒸散發(fā)是指所有地表水轉(zhuǎn)化為水蒸氣的過程，包括樹冠截留的水分蒸發(fā)、蒸騰和升華及土壤水的蒸發(fā)。蒸散發(fā)是水分轉(zhuǎn)移出流域的主要途徑，在許多江河流域，蒸發(fā)量都大于徑流量。準確地評價蒸散發(fā)量是估算水資源量的關鍵，也是研究氣候和土地覆蓋變化對河川徑流影響的關鍵。（1）潛在蒸散發(fā)模型提供了PenmanMonteith、PriestleyTaylor和Hargreaves三種計算潛在蒸散發(fā)的方法，另外還可以使用實測資料或已經(jīng)計算好的逐日潛在蒸散發(fā)資料。一般采用PenmanMonteith方法來計算流域的潛在蒸散發(fā)。（2）實際蒸散發(fā)實際蒸散發(fā)以潛在蒸散發(fā)為計算基礎。在計算流域?qū)嶋H蒸散發(fā)量的時候，模型首先計算植物冠層截留水分的蒸發(fā)，然后計算最大蒸騰量、最大升華量和最大土壤蒸發(fā)量，最后計算實際的升華量和土壤水分蒸發(fā)量。2021.02.05歐陽科創(chuàng)編2021.02.05歐陽科創(chuàng)編 2021.02.05(3)冠層截留蒸發(fā)量模型在計算實際蒸發(fā)時假定盡可能蒸發(fā)冠層截留的水分，如果潛在蒸發(fā)量二小于冠層截留的自由水量1.，貝卜式中：三為某日流域的實際蒸發(fā)量，mm；'為某日冠層自由水蒸發(fā)量，mm；巳為某日的潛在蒸發(fā)量，mm;為某日植被冠層自由水初始含量，mm；,：為某日植被冠層自由水終止含量，mm。如果潛在蒸發(fā)量二大于冠層截留的自由水含量1.貝：當植被冠層截留的自由水被全部蒸發(fā)掉，繼續(xù)蒸發(fā)所需的水分就會從植被和土壤中得到。(4)植物蒸騰假設植物生長在一個理想的條件下，植物蒸騰可用以下表達式計算：當0這LAI這3.0時當LAI>3,0時式中：巳為某日最大蒸騰量，mm/?e：為植被冠層自由水蒸發(fā)調(diào)整后的潛在蒸發(fā)匚一二6％mm；LAI為葉面積指數(shù)。因為沒有考慮到植物下面圖層的含水量問題，由此公式計算處的蒸騰量可能比實際蒸騰量要大一些。(5)土壤水分蒸發(fā)在計算土壤水分蒸發(fā)時，首先區(qū)分出不同深度土壤層所需要的蒸發(fā)量，土壤深度層次的劃分決定土壤允許的最2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05大蒸發(fā)量，可由下式計算：式中：1…為z深度處蒸發(fā)需要的水量，mm；z為地表以下土壤的深度，mm。表達式中的系數(shù)是為了滿足50%的蒸發(fā)所需水分來自土壤表層10mm,以及95%的蒸發(fā)所需水分來自0?100mm土壤深度范圍內(nèi)。土壤水分蒸發(fā)所需要的水量是有土壤上層蒸發(fā)需水量與土壤下層蒸發(fā)需水量決定的：式中：卜為ly層的蒸發(fā)需水量，mm;Z"二:為土壤下層的蒸發(fā)需水量，mm；「*二一為土壤上層的蒸發(fā)需水量，mm。土壤深度的劃分假設50%的蒸發(fā)需水量由0?10mm內(nèi)土壤上層的含水量提供，因此100mm的蒸發(fā)需水量中50mm都要由10mm的上層土壤提供，顯然上層無法滿足需要，這就需要建立一個系數(shù)來調(diào)整土壤層深度的劃分，以滿足蒸發(fā)需水量，調(diào)整后的公式可以表示為：式中：esco為土壤蒸發(fā)調(diào)節(jié)系數(shù)，該系數(shù)是SWAT為調(diào)整土壤因毛細作用和土壤裂隙等因素對不同土層蒸發(fā)量二提出的，對于不同的esco值對應著相應的土壤層劃分深度。2.1.3土壤水滲入到土壤中的水有多種不同運動方式。土壤水可以被植物吸收或蒸騰而損耗，可以滲透到土壤底層最終補給2021.02.05歐陽科創(chuàng)編

2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05地下水，也可以在地表形成徑流，即壤中流。由于主要考慮徑流量的多少，因此對壤中流的計算簡要概括。模型采用動力儲水方法計算壤中流。相對飽和區(qū)厚度：計算公式為：式中：為土壤飽和區(qū)內(nèi)可流出的水量，mm；；血為山坡坡長，m；中,為土壤可出流的孔隙率；］表示土壤層總孔隙度，即山.與土壤層水分含量達到田間持水量的孔隙度中:.?之差。山坡出口斷面的凈水量為：式中：「為出口斷面處的流速，mmh。其表達式為：式中：k為土壤飽和導水率，mmh；山為坡度。總結(jié)上面表達式，模型中壤中流最終計算公式為：??地下水模型采用以下表達式來計算流域地下水:式中：：為第i天進入河道的地下水補給量， mm；?！笧榈冢╥）天進入河道的地下水補給量，mm；”為時間步長，d；時間步長，d；土壤侵蝕模型2021.02.05歐陽科創(chuàng)編2021.02.05歐陽科創(chuàng)編 2021.02.05根據(jù)天氣條件，泥沙輸移量可以高于或者低于輸移能力，導致沉積過量的泥沙通過渠道侵蝕再懸浮輸移泥沙。流速方程為：式中：?為流量，m3/s；為渠道寬度，m；.」為徑流深，m。對于低于齊岸深度的徑流，徑流深使用Manning方程來計算，假定渠道寬度遠大于深度：式中：”為渠道曼寧系數(shù)；「為渠道坡度，m/m。由于降水和徑流產(chǎn)生的土壤侵蝕是用MUSLE方程來計算的，MUSLE是修正的通用土壤流失方程（USLE）。USLE方程是通過降水動能函數(shù)預測年均侵蝕量，而在 MUSLE中，用徑流因子代替降水動能，改善了泥沙產(chǎn)量的預測，這樣就不需要泥沙輸移系數(shù)，并且可以將方程用于單次暴雨事件，因為徑流因子是先行濕度和降水動能的函數(shù)。USLE中需要輸移系數(shù)是因為降水動能因子表示的能量只在作用流域內(nèi)起作用。修正的通用土壤流失方程為；式中："「為土壤侵蝕量，t；Q為地表徑流，mm/h；」」為洪峰流量，m3/s；7丁為水文響應單元的面積，hm2；1,ME為土壤侵蝕因子；「le為植被覆蓋和管理因子；「UE為水土保持措施因子；Ffle為地形因子；口內(nèi)。為粗碎屑因子。2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05土壤侵蝕當其他影響侵蝕的因子不變時，K因子反映不同類型土壤抵抗侵蝕力的高低。它與土壤物理性質(zhì)的影響，如機械組成、有機質(zhì)含量、土壤結(jié)構(gòu)、土壤滲透性等有關。當土壤顆粒粗、滲透性大時，K值就低，反之則高；一般情況下k值得變幅在0.02?0.75之間。K值得直接測定方法是：在標準小區(qū)（坡長為22.1m，寬為1.83m，坡度為9%）上沒有任何植被，完全休閑，無水土保持措施，降水后收集由于坡面徑流而沖蝕到集流槽內(nèi)的土壤，烘干、稱重，由公式計算得到卜：值。試驗測算KK式中：記為顆粒尺度參數(shù)；。式中：&償為土壤有機質(zhì)因子；L....式中：⑺為粒徑在0.05?2.00mm沙粒的百分含量；…:為粒徑在0.002?0.05mm的淤泥、細砂百分含量；,為粒徑小于0.002mm的粘土百分含量；；.為各土壤層中有機碳含量，％。植被覆蓋因子（7』植被覆蓋和管理因子c：、：.表示植物覆蓋和作物栽培措施對防止土壤侵蝕的綜合效益，其含義是在地形、土壤、降水條件相同的情況下，種植作物或林草地的土地與連續(xù)休2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05養(yǎng)地土壤流失量的比值，最大取值為1.0.由于植被覆蓋受植物生長期的影響，SWAT模型通過下面的方程調(diào)整植被覆蓋和管理因子式中：，「常是最小植被覆蓋和管理因子值；是地表植物殘留量，kg/hm2。最小C因子可以由已知年平均C值，通過以下方程計算。式中：』皿…表示不同植被覆蓋的年均C值。水土保持措施因子 』水土保持因子LtLE是指有保持措施的地表土壤流失與不采取任何措施的地表土壤流失的比值，這里的水土保持措施包括等高耕作、帶狀種植和梯田。等高耕作對于中低強度的降水侵蝕具有保護水土流失的作用，但對于高強度的降水其保護作用則很小，等高耕作對坡度為3%?8%之間的土地非常有效。地形因子地形因子匚g』：的計算公式如下：式中："為坡長；由為坡長指數(shù)；胴為坡度。坡長指數(shù)「.的計算公式如下：式中：口為水文相應單元的坡度，R-二⑴」2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05225CFR(■因子□R。因子是通過下面公式計算的：式中：為第一層土壤中礫石的百分比，％。污染負荷模型SWAT模型可以模擬不同形態(tài)氮的遷移轉(zhuǎn)化過程，地表徑流流失、入滲淋失、化肥輸入等物理過程，有機氮礦化、反硝化等化學過程以及作物吸收等生物過程，氮可以分為有機氮、作物氮和硝酸鹽氮三種化學狀態(tài)，氮的生物固定、有機氮向無機氮的轉(zhuǎn)化以及溶解性氮隨側(cè)向壤中流的遷移等過程，有機氮又被劃分為活潑有機氮和惰性有機氮兩種狀態(tài)，以及銨態(tài)氮揮發(fā)過程的模擬。硝態(tài)氮污染負荷模型硝態(tài)氮主要隨地表徑流、側(cè)向流或滲流在水體中遷移，要計算隨水體遷移的硝態(tài)氮量必須先計算自由水中硝態(tài)氮濃度，用這個濃度乘以各個水路流動水的總量，即可得到從土壤中流失的硝態(tài)氮總量。自由水部分的硝態(tài)氮濃度可用下面公式計算：式中：…㈠■為自由水中硝態(tài)氮濃度，kg/mm:,..為土壤中硝態(tài)氮的量，kg/hm2；省血於.為土壤中自由水的量，mm「I為孔隙度；"「為土壤飽和含水量。2021.02.05歐陽科創(chuàng)編2021.02.05歐陽科創(chuàng)編 2021.02.05（1）通過地表徑流流失的溶解態(tài)氮計算公式：式中：：'.為通過地表徑流流失的硝態(tài)氮，kg/hm2；3.為硝態(tài)氮滲流系數(shù)；Q…為地表徑流，mm。（2）通過側(cè)向流流失的溶解態(tài)氮的量計算公式：對于地表10mm圖層：對于10mm以下的土層：式中："二”也為通過側(cè)向流流失的硝態(tài)氮，kg/hm2”,為硝態(tài)氮滲流系數(shù)；為自由水的硝態(tài)氮濃度，kg/mm；1"：i為側(cè)向流，mm。（3）通過滲流流失的溶解態(tài)氮量計算公式：式中：入二」…為通過滲流流失的硝態(tài)氮，kg/hm2；為自由水的硝態(tài)氮濃度，kg/mm；'［為滲流,mm。有機氮污染負荷模型有機氮通常是吸附在土壤顆粒上隨徑流遷移的，這種形式的氮負荷與土壤流失量密切相關，土壤流失量直接反映了有機氮負荷。有機氮隨土壤流失的輸移量計算公式為：式中：p…為有機氮流失量，kg/hm2；、「為有機氮在表層土壤中的濃度，kg/t；⑴為土壤流失量，t；。為水文相應單元的面積，hm2；為氮富集系數(shù)，氮富集系數(shù)是隨土壤流失的有機氮濃度和土壤表層有機氮濃度的比值。2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05計算富集系數(shù)的公式如下：式中：：:“為地表徑流中泥沙含量。：「的計算公式如下：式中：田為土壤流失量，t；，'「為水文相應單元面積，hm2；Q”為地表徑流,mm。溶解態(tài)磷污染負荷模型溶解態(tài)磷在土壤中的遷移主要是通過擴散作用實現(xiàn)的，擴散是指離子在微小尺度下由于濃度梯度而弓I起的溶質(zhì)遷移，由于溶解態(tài)磷不很活躍，所以由地表徑流以溶解態(tài)形式帶走的土壤表層的磷很少，地表徑流輸移的溶解態(tài)磷可由下面公式計算：式中： 為通過地表徑流流失的溶解態(tài)磷， kg/hm2；p：士-.…為土壤中溶解態(tài)磷，kg/hm2；為土壤溶質(zhì)密度，mg/m3；h…為表層土壤深度，mm/?上；二為土壤磷分配系數(shù)，表層土壤中溶解態(tài)磷的濃度和地表徑流中溶解態(tài)磷濃度的比值。有機磷和礦物質(zhì)磷污染負荷模型有機磷和礦物質(zhì)磷通常是吸附在土壤顆粒上通過徑流遷移的，這種形式的磷負荷與土壤流失量密切相關，土壤流失量直接反映了有機磷和礦物質(zhì)磷負荷，有機磷和礦物質(zhì)磷隨土壤流失輸移量計算公式為：2021.02.05歐陽科創(chuàng)編2021.02.05歐陽科創(chuàng)編 2021.02.05式中：m.為有機磷流失量，kg/hm2；,為有機磷在表層土壤中的濃度，kg/t；田為土壤流失量，t；'為水文相應單元的面積，hm2；；二為磷富集系數(shù)。河道演算模型模型中定義的河道均是明渠流，SWAT用曼寧公式來定義河道糙率和水流流速。水流在河道中演進過程使用變儲量演算法或馬斯金根法，兩種方法都是動力波方程。在模擬中一般使用馬斯金根法來計算。（1）馬斯金根法馬斯金根法假設河道內(nèi)水體形狀是由一個楔形蓄水體和一個棱形蓄水體組成，如下圖 所示。圖 ？?河道水體示意圖當洪水波行進到某個河段槽，入流量大于出流量便形成了楔形蓄水體。當洪水波退去，在河段槽便出現(xiàn)了出流量大于入流量的負楔蓄體。另外對于楔蓄水體，河段槽內(nèi)始終包含一個體積為流域長度上橫截面不變的棱柱狀水體?？偟男钏萘繛椋菏街校盒樾钏萘?，m??；，爭為入流量，m??s；：「川為出流量，m??s;、為穩(wěn)定流情況下的河段傳播時間；X為流量比重因素。2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05該公式可以重新整理為如下形式：流量比重因素的下限為0.0,上限為0.5。這個因子是楔蓄量的函數(shù)。對于水庫式蓄水，沒有楔蓄，'.：」「；而對于一個完整的楔蓄，.、心/?對于河流，X落在0.0和0.3之間，其平均值接近0.2。對于蓄水容量的定義可以加入連續(xù)公式并簡化為：式中：」為該時間段開始時的入流量，m3/s/」為該時間段結(jié)束時的入流量，m3/s；」 為該時間段結(jié)束時的出流量，m3/s；「為該時間段結(jié)束時的出流量,m3/s。其中，「"小」1。為了用體積單位表示所有值，在蓄水容量公式兩端乘以該時間段得到：為了保持數(shù)值穩(wěn)定和避免出現(xiàn)負出流量的計算，必須滿足以下條件：流量比重因素X的值由使用者輸入，蓄水時間常數(shù)的值估算如下：式中：k為穩(wěn)定流情況下的河段傳播時間，s；z和川為權(quán)重系數(shù)，由使用者輸入；廠心訕為穩(wěn)定劉情況下渠道蓄滿水的河段傳播時間，s；，；面山為渠道蓄滿1/10水量時河段傳播時間，s。要計算匕心和一；面山，Cunge于1969年提出一個公式：式中：廠為渠道長度，km1為指定深度處的波速，m/s。2021.02.05歐陽科創(chuàng)編

2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05波速的計算公式為：其中流速”.由曼寧公式求解，將曼寧公式代入上式得：式中：R二為濕周，m「ip二為河段坡度，%;為曼寧系數(shù)；、為流速，m/s。(2)變儲量演算法對于一個給定的河段，儲量演算基于連續(xù)方程' ，' ，,：、可寫為：式中：%為穩(wěn)定流情況下的河段傳播時間，s；.為該時間段開始時的入流量，m3/s；l為該時間段結(jié)束時的入流量，m3/s；「為該時間段開始時的出流量，m3/s；「.為該時間段結(jié)束時的出流量，m3/s；」.為該時間段開始時的蓄水容量，m3；'「:：為該時間段結(jié)束時的蓄水容量，m3。演進時間是由渠道中的水容量除以水流流量：式中：、一為蓄水容量，m3；“為出流量，m3/s。聯(lián)合以上兩式，簡化為：2A其中，y為蓄水系數(shù)，乩= ;為入流平均流時,1+勺運2里,“?二，m3/s。整理上式得：為用體積單位來表達所有變量，公式兩邊都乘以時間段，則有：2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.053^件操作輸入數(shù)據(jù)要創(chuàng)建SWAT數(shù)據(jù)集，ArcSWAT需要訪問提供流域信息的準確類型的ArcGIS的柵格、矢量和數(shù)據(jù)文件。這些必需的地圖主題和數(shù)據(jù)文件必須在運行SWAT之前準備好。空間數(shù)據(jù)空間數(shù)據(jù)包括必需的ARCSWAT空間數(shù)據(jù)和可選的ARCSWAT空間數(shù)據(jù)。SWAT空間數(shù)據(jù)集可以以任何投影類型來創(chuàng)建，但所有的地圖必須用同一投影。（1）必需的ARCSWAT空間數(shù)據(jù)DEM，ESRIGRIDFormat。DEM的高程值可以用整型或者實數(shù)型。決定地圖分辨率的單位不必與高程的單位保持一致。例如地圖分辨率可以是米（m），而高程可以是英尺（feet）。地圖分辨率單位必須定義為以下幾種：米（m），公里（km），英尺（feet），碼（yards），英里（miles）。高程單位必須定義為以下幾種：米（m），厘米（cm），碼（yards），英尺（feet），英寸（inches）。LandCover/LandUse，ESRIGRID，Shapefile，orFeatureClassFormat。土地利用圖中的類別需要重新2021.02.05歐陽科創(chuàng)編2021.02.05歐陽科創(chuàng)編 2021.02.05分類為SWAT需要的土地利用類型。用戶可以用三種方法重新分類土地利用，第一種是創(chuàng)建地圖時用美國地質(zhì)勘探局的分類代碼；第二種是加載土地利用圖時為每一分類選定SWAT土地利用類型；第三種是為土地利用圖的不同分類創(chuàng)建一個能夠識別4位SWAT代碼的lookup表。Soil，ESRIGRID，Shapefile，orFeatureClassFormat。用戶需要輸入SWAT土壤文件，在創(chuàng)建工程之前把圖中每一類的土壤數(shù)據(jù)輸入UserSoildatabase中。要對土壤圖的土壤類型進行重分類，信息必須在界面中手動輸入，列出土壤信息100kup表（2）可選的ARCSWAT空間數(shù)據(jù)集DEMMask，ESRIGRID，ShapefileorFeatureClassFormat。模型可以加載一個mask在DEM之上。STREAMS，ShapefileorFeatureClassFormat。在那些地勢很低的地方，DEM地圖網(wǎng)格不能精確地判斷河流的位置，就需要有河網(wǎng)描繪的 shapefile線性文件。User-DefinedWatersheds，ShapefileorFeatureClassFormat。用戶可以自定義子流域劃分加載到模型中去，如果選擇自定義子流域也必須自定義河網(wǎng)。User-DefinedStreams，ShapefileorFeatureClass2021.02.05歐陽科創(chuàng)編2021.02.05歐陽科創(chuàng)編 2021.02.05Format。用戶自定義子流域后，模型將不會進行河網(wǎng)提取，需要用戶自行定義。屬性數(shù)據(jù)屬性數(shù)據(jù)包括：土壤屬性數(shù)據(jù)、氣象資料數(shù)據(jù)和水文數(shù)據(jù)。SWAT模型通過三個個數(shù)據(jù)庫來存儲屬性數(shù)據(jù)，其中氣象資料數(shù)據(jù)通過模型自帶的數(shù)據(jù)庫進行存儲，另外自定義的土壤屬性數(shù)據(jù)和氣象站參數(shù)數(shù)據(jù)通過兩個附加的數(shù)據(jù)庫進行存儲，水文數(shù)據(jù)用于模擬結(jié)果率定。（1）土壤屬性數(shù)據(jù)模型需要的土壤屬性數(shù)據(jù)包括兩大類：土壤物理屬性和土壤化學屬性。土壤的物理屬性對土壤剖面中水和氣的運動狀況起決定作用，并顯著影響水文響應單元（HRU）中水文過程的模擬。土壤物理屬性參數(shù)主要包括土層厚度、土壤層數(shù)、密度、土壤顆粒大小分布和土壤飽和水力傳導率等。由于獲得的土壤屬性數(shù)據(jù)難以完全滿足構(gòu)建SWAT模型土壤庫的要求，因此，部分無法直接獲取的參數(shù)如土壤可利用有效水、飽和水力傳到率等采用美國農(nóng)業(yè)部開發(fā)的土壤水特性程序SPAW進行估算；部分參數(shù)通過對前人的研究成果和經(jīng)驗公式獲得。估算結(jié)果僅簡略反映土壤參數(shù)特征，在模型參數(shù)校準的過程中還需做進一步調(diào)整。2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05(2)氣象、水文資料數(shù)據(jù)SWAT模型所需要的氣象數(shù)據(jù)包括氣象站點的位置以及各氣象站點的實測數(shù)據(jù)。氣象站點包括雨量站、濕度站、溫度站、風速站、太陽輻射站，其實測數(shù)據(jù)分別為日降水量、相對濕度、日最高/最低氣溫、風速和太陽輻射。本文采用北京站的氣象資料來構(gòu)建模型“天氣發(fā)生器”所需的數(shù)據(jù)。SWAT模型自帶的天氣發(fā)生器WXGEN，其功能有二，一是用于生成氣候數(shù)據(jù)，通過大量基礎氣象數(shù)據(jù)構(gòu)建完成之后，可以用來生成任意年份的氣象數(shù)據(jù)，二是對缺失的數(shù)據(jù)進行填補。輸出數(shù)據(jù)SWAT模型每次模擬都會生成五個輸出文件：輸入?yún)R總文件(input.std)、輸出匯總文件(output.std)、HRU輸出文件(output.hru)、子流域輸出文件(output.sub)、河道輸出文件(output.rch)。3.2.1output.hruHRU輸出文件包含了流域中每一個水文響應單元的信息，文件中輸出數(shù)據(jù)介紹如下：MON:時間步長，根據(jù)用戶選擇可以為日、月、年；AREA：HRU的面積，km2；2021.02.05歐陽科創(chuàng)編2021.02.05歐陽科創(chuàng)編 2021.02.05PRECIP:時間步長內(nèi)HRU的降雨量，mm；SNOFALL：時間步長內(nèi)HRU的降雪量，mm；SNOMELT：時間步長內(nèi)HRU的冰雪融化量，mm；IRR:灌溉用水量，mm；PET：潛在蒸散發(fā)量，mm；ET：實際蒸散發(fā)量，mm；SWINIT:土壤初始含水量，mm；SWEND:時間步長后土壤含水量，mm；PERC：植物蒸騰水量，mm；GWRCHG：入滲量，mm；DARCHG:深層地下水補給量，mm；REVAP:淺層地下水向表層土壤回流量，mm；SAIRR:淺層地下水灌溉量，mm；DAIRR:深層地下水灌溉量，mm；SAST:淺層地下水儲水量，mm；DAST:深層地下水儲水量，mm；SURQGEN:地表徑流量，mm；SURQCNT:地表徑流進入河道水量，mm；TLOSS：輸移損失水量，mm；LATQ:側(cè)向流水量，mm；GWQ:地下水補充河道水量，mm；WYLD:總產(chǎn)流量，包括地表徑流、地下水、側(cè)向流,2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05扣除損失水量，mm；DAILYCN:日徑流曲線值；TMP-AV：日平均溫度，攝氏度；TMP-MX：日最高溫度，攝氏度；TMP-MN：日最低溫度，攝氏度；SOL-TMP：土壤溫度，攝氏度；SOLAR：日輻射量，MJ/m2；SYLD:產(chǎn)沙量，t/ha；USLE:土壤流失量，t/ha；N-APP：氮肥施用量，kg/ha；P-APP:磷肥施用量，kg/ha；NRAIN:氮素隨降雨進入土壤量，kg/ha；NFIX:植物固氮量，kg/ha；ORGN:有機氮產(chǎn)生量，kg/ha；ORGP:有機磷產(chǎn)生量，kg/ha；SEDP:隨泥沙磷流失量，kg/ha；NSURQ:地表徑流中硝態(tài)氮量，kg/ha；NLATQ:側(cè)向流中硝態(tài)氮含量，kg/ha；NO3L：土壤中被沖刷出的硝態(tài)氮，kg/ha；NO3GW:地下水流失硝態(tài)氮量，kg/ha；SOLP:溶解態(tài)磷產(chǎn)生量，kg/ha；P-GW:地下水溶解態(tài)磷流失量，kg/ha。2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.053.2.2output.subMON:時間步長，據(jù)用戶選擇可以為日、月、年;AREA:子流域面積，km2；PRECIP:降水量，mm；SNOMELT：冰雪融化量，mm；PET：潛在蒸發(fā)量，mm；ET：實際蒸發(fā)量，mm；SW:土壤含水量，mm；PERC：入滲量，mm；SURQ:地表徑流進入河道水量，mm；GWQ:地下水補給河道水量，mm；WYLD:流域產(chǎn)水量，mm；SYLD:流域產(chǎn)沙量，t/ha；ORGN:流域有機氮產(chǎn)生量，kg/ha；ORGP:流域有機磷產(chǎn)生量，kg/ha；NSURQ:地表徑流中硝態(tài)氮含量，kg/ha；SOLP:流域溶解態(tài)磷產(chǎn)生量，kg/ha；SEDP:流域礦化磷產(chǎn)生量，kg/ha；3.2.3output.rchMON:時間步長，據(jù)用戶選擇可以為日、月、年;AREA:河道集水面積，km2；2021.02.05歐陽科創(chuàng)編2021.02.05歐陽科創(chuàng)編 2021.02.05FLOWIN:日平均入流量，m3/s；FLOWOUT：日平均出流量，m3/s；EVAP:河道蒸發(fā)量，m3/s；TLOSS:河道輸移損失水量，m3/s；SEDIN:泥沙入流量，t；SEDOUT:泥沙出流量，tSEDCONC：泥沙淤積量，mg/L；ORGNIN:有機氮輸入量，kg；ORGNOUT:有機氮輸出量，kg；ORGPIN:有機磷輸入量，kg；ORGPOUT:有機磷輸出量，kg；NO3IN:硝態(tài)氮輸入量，kg；NO3OUT:硝態(tài)氮輸出量，kg；NH4IN:氨氮輸入量，kg；NH4OUT:氨氮輸出量，kg；NO2IN:亞硝態(tài)氮輸入量，kg；NO2OUT:亞硝態(tài)氮輸出量，kg；MINPIN:礦化磷輸入量，kg；MINPOUT:礦化磷輸出量，kg；CBODIN:生物需氧量輸入量，kg；CBODOUT:生化需氧量輸出量，kg；DISOXIN:溶解氧輸入量，kg2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05DISOXOUT:溶解氧輸出量，kg；流域模擬基于DEM的水文參數(shù)計算DEM包括地理平面坐標和高程數(shù)據(jù)兩類信息，由DEM自動提取水系或子流域會遇到兩個問題，即凹陷與平坦處水流方向的確定和水道起始點位置的確定。水道起始點位置處理采用最小水道集水面積閾值的方法形成河流水系的基本水文過程。運用GIS軟件，對流域DEM進行分析和處理，可得到柵格格式的水流流向，流域分水線，自動生成河網(wǎng)及其子流域，河道和子流域編碼、面積，河網(wǎng)結(jié)構(gòu)的拓撲關系等。（1）DEM填洼采用墊高添平方法對流域DEM模型中存在的洼地進行凹陷處理，將洼地劃分為阻擋型和凹陷型分別處理，并運用高程增量疊加算法來設定平坦部分的流向，由此得到一個與原數(shù)字高程模型相對應的無凹陷的口￡片模型。（2）水流流向確定計算DEM模型每個網(wǎng)格單元與其相嶺的八個網(wǎng)格單元之間的坡度，然后按最陡坡度原則確定單元水流流向，最終形成包括所有網(wǎng)格單元在內(nèi)的水流方向新矩陣。2021.02.05歐陽科創(chuàng)編2021.02.05歐陽科創(chuàng)編 2021.02.05（3）集水面積計算沿坡度最陡原則確定的水流路徑可以計算任何網(wǎng)格單元上的坡面集水面積，其集水面積的量值以網(wǎng)格數(shù)目的多少表示，從而產(chǎn)生包含每個網(wǎng)格單元上坡面匯水面積的新矩陣。（4）流域分水線劃分首先給出閉合流域出口斷面的準確地理位置，即DEM模型網(wǎng)格單元上的行列坐標值；流域出口點位置確定后，可以按照確定網(wǎng)格流向的方法來勾繪流域邊界。（5）河網(wǎng)生成與編碼定義最小河道集水面積閾值，當集水面積超過此閾值時，則這些網(wǎng)格點就確定為河道。進行河網(wǎng)編碼時，首先從流域出口斷面向上游靠左搜索，對每個節(jié)點依次遞增編碼，直至河道起始點，再從河道起點反方向搜尋，直到找到新的節(jié)點或河道起點位置，同時按遞增的順序規(guī)則編碼，依次類推，最終搜尋至流域出口斷面。（6）特征參數(shù)提取生成完整的河網(wǎng)后，可以確定每一河段左右岸集水面積、河道上游末端節(jié)點及相應的子流域分水線，從而建立河網(wǎng)節(jié)點、河段和子流域之間的關系，包括河段坡度、高程、上游集水面積及其他信息。流域勾繪過程中，依據(jù)DEM計算流域坡長、河流流向和流域水系。2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05子流域劃分建立流域模擬的第一步是將流域劃分為子流域?？臻g數(shù)據(jù)的輸入時基于網(wǎng)格，網(wǎng)格劃分是子流域勾繪的基礎。SWAT模型基于DEM數(shù)據(jù)，在每個子流域上應用概念性集總模型進行計算，再進行匯流演算，最后求得出口斷面流量和污染物質(zhì)量。每個子流域包含無限制的HRU（至少一個）、河段（一個）、主渠道網(wǎng)絡中的池塘（可選）和點源排放（可選）等。子流域保持流域的地理位置并同其他的子流域具有空間聯(lián)系。子流域勾繪可以從地表地形定義的子流域邊界中得到，從而子流域的整個面積都流向子流域出口，或者子流域勾繪可以從網(wǎng)格邊界中獲得，河段或者主河道同流域中的每個子流域想聯(lián)系。應用GIS可以將流域劃分為多個水文學上相聯(lián)系的子流域以便在SWAT流域模擬中使用。SWAT流域工具應用并擴展了ArcView的SpatialAnalystExtension的功能，可以在DEM數(shù)據(jù)基礎上進行子流域勾繪。定義的參數(shù)范圍會影響所生成的子流域大小和數(shù)目，勾繪過程要求 grid格式的DEM數(shù)據(jù)。應用基于柵格DEM進行流域勾繪的八流點法則，產(chǎn)生水流矢量柵格來填充凹陷點；通過計算柵格中到達任一單元的水流貢獻單元數(shù)產(chǎn)生匯流柵格。那些潛在的座位河網(wǎng)部分的單元會有較大的匯流值，而那些接近流域邊界和坡面漫流控制的單元則會具有較低的匯流值。2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05選擇用來模擬子流域的數(shù)目取決于整個流域的大小，輸入數(shù)據(jù)的空間詳細程度以及達到研究目標所要求的詳細度。在劃分子流域時，要注意流域的地貌特征（坡度、坡長、渠道長度和渠道寬度等）在子流域水平上計算和匯總。河流的支流數(shù)目受到形成支流的匯流網(wǎng)格單元數(shù)目閾值的限制，通過對河流支流下游邊界的定位可以對子流域出口點進行定義。SWAT提供兩種附加的設置工具：DEM屬性和閾值面積。這兩種設置工具是相互聯(lián)系的，第一種設置DEM的水平和垂直單位，這對計算地貌系數(shù)以及從以公頃為單位的數(shù)值中計算匯流單元數(shù)目閾值有重要作用。模型中流域勾繪的其他部分是一套獨立的工具按鈕，可以單擊流域視圖上相應的位置或者輸入點位置來增加出口或者入口點，也可以去除選中的出口點。當出水口確定后，定義主流域出水口。出水口的類型有7種，分別為：L連接河流增加的出水口、O表格增加的出水口、T手動增加的出水口、P手動增加的點源排放口、D表格增加的點源排放口、I表格增加的排水流域入口、W手動增加的排水流域入口。然后跟蹤每個柵格單元的流向直到遇到流域出口單元或者到達DEM網(wǎng)格的邊界，完成子流域勾繪?；诹饔虻牡匦翁卣鲗⒘饔騽澐譃樽恿饔?，保留了河水、泥沙和化學物質(zhì)的實際演算所要求的自然河流2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05的路徑、邊界和渠道。勾繪完成后，會在當前的項目中加入詳細的地形報告，隨后結(jié)果專題圖也會加入到流域視圖中，包括：子流域、河流、出水口和水庫（可選）。地形報告描述了流域內(nèi)合子流域單元的高程分布；結(jié)果專題圖含有流域的特征參數(shù)。水文響應單元當子流域勾繪完成后，可以選擇模擬每個子流域內(nèi)單個的土壤類型、土地利用將子流域劃分為多個水文響應單元（HRUs）。水文響應單元是子流域的一部分，含有唯一的土地利用和土壤屬性，被假定為在子流域中有統(tǒng)一的水文行為。HRU并不同于田間小區(qū)，而是子流域內(nèi)特定的土地利用和土壤類型的總面積，而單個田間小區(qū)具有特殊的土地利用和土壤類型，在整個子流域內(nèi)可能是離散的，這些面積聚集在一起形成HRU。HRU應用在模型中，通過聚集所有相似的土壤類型和土地利用面積構(gòu)成耽擱的響應單元，從而簡化了模型運行。HRU的好處是增加了對子流域負荷預測的準確性。不同種類的植被生長和發(fā)育有很大的不同，當考慮了子流域內(nèi)植被覆蓋的多樣性時，從子流域中進入主渠道的徑流量會預測的更準確，對給定的子流域應該有1?10個HRU。當2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05對子流域內(nèi)多個HRU進行模擬時，子流域內(nèi)水文循環(huán)的陸地階段也進行了模擬，并對子流域內(nèi)所有的HRU的負荷進行了總匯，子流域的凈負荷通過流域河網(wǎng)進行演算。HRU的概念表明，假定在子流域的HRU之間不存在相互影響，每個HRU的負荷分別進行計算，然后匯總在一起確定子流域的總負荷。可以對每個子流域分配單個HRU或者多個HRU，如果選擇一個子流域單個HRU，則HRU有子流域內(nèi)的主要土地利用分類和土壤類型來決定；如果選擇多個HRU，可以設定土地利用和土壤數(shù)據(jù)的敏感性，以用來確定每個子流域內(nèi)HRU的數(shù)目和種類，對多個HRU，要定義土地利用和土壤類型面積閾值。如果一種土地利用類型同另一種土地利用的相互作用是重要的，可以將這些土地利用面積定義為子流域而不是HRU，只有在子流域水平上，空間關系才能被詳細描述。參數(shù)敏感度分析在分布式水文模型的模擬計算中，對現(xiàn)有的分布式水文模型與環(huán)境模擬模型的精度要求也越來越高。模型發(fā)展趨勢表現(xiàn)為兩方面：一是模型本身的適用范圍逐漸擴展，空間上擴展到流域尺度，時間上精確計算到日；二是模型模擬范圍逐漸拓展，在陸地循環(huán)和河道循環(huán)中計算水文要素和水質(zhì)數(shù)據(jù)。但這需要模型校正為前提，一般的模型校2021.02.05歐陽科創(chuàng)編2021.02.052021.02.05歐陽科創(chuàng)編2021.02.05正需要大量的數(shù)據(jù)，如水文站、監(jiān)測站的實測數(shù)據(jù)或者是實驗室試驗數(shù)據(jù)等。事實上，模型的準確性和可靠性是有限的。模擬精確度較高的是模擬小區(qū)域不透水的集水面上徑流的水文模型，誤差在百分之幾左右；可靠性最差的是模擬大流域的水質(zhì)模型，誤差可能達到一個數(shù)據(jù)量級甚至更大。根據(jù)非點源污染產(chǎn)出的特點，首先率定的是水文過程，然后是泥沙和污染物的遷移。在水文或者侵蝕過程出現(xiàn)的誤差，會轉(zhuǎn)移并且擴大到其他過程，致使非點源污染負荷模擬誤差增大。對分布式模型的大量參數(shù)進行分析。一是運用摩爾斯分類篩選法，即“微擾動”分析方法，考察大量的輸入?yún)?shù)與輸出結(jié)果之間的相關性，對模型大量的參數(shù)進行定性分級篩選，繼而計算經(jīng)驗敏感指標值；二是在篩選分類的基礎上運用付利葉敏感型檢驗法計算相關性較好的一些參數(shù)計算敏感度指標值，分析其對輸出結(jié)果的影響程度，即模型輸出結(jié)果對參數(shù)的敏感性。SWAT模型采用后者進行參數(shù)敏感性分析計算。模型校準與驗證當模型的結(jié)構(gòu)和輸入?yún)?shù)初步確定后，需要對模型進行參數(shù)校準和驗證。通常將所使用的資料系列分為兩部2021.02.05歐陽科創(chuàng)編2021.02.05歐陽科創(chuàng)編 2021.02.05分：其中一部分用于模型參數(shù)校準，而另一部分則用于模型驗證。參數(shù)校準是模型驗證的重要一步，他能夠揭示模型在設計和執(zhí)行過程中的缺陷，在不能或者難以獲得必要的參數(shù)值時，參數(shù)校準是相當有用的。標準的參數(shù)校準過程是尋找實測和預測狀態(tài)的細微差別，并通過統(tǒng)計的泥河度來衡量模型的適用性。當模型參數(shù)校準完成后，應用參數(shù)校準數(shù)據(jù)集意外的實驗數(shù)據(jù)或者現(xiàn)場觀測數(shù)據(jù)對模型模擬值進行對比分析與驗證，以評價模型的適用性。為證明一個經(jīng)過校準了得模型是否可以使用，可以用不同于校準的數(shù)據(jù)來進行驗證。在SWAT模擬中，可以選擇相對誤差Re、相關系數(shù)R2和NashSuttcliffe系數(shù)來評價模型的適用性。相對誤差計算公式為：式中：3為模型模擬相對誤差；2為模擬值；。為實測值。若八為正值，說明模型預測或模擬值偏大；若吐為負值，模型預測或模擬值偏??；若此、則說明模型模擬結(jié)果與實測值正好吻合。相關系數(shù)計算公式為：式中：Q為實測值；。為模擬值；。，為平均值；?為實測數(shù)據(jù)個數(shù)。決定系數(shù)是實測值與模擬值做線性回歸得到，表示兩者之間的線性關系，用于評價兩者之間的吻合程度。如果=1，表示觀測值與模擬值呈完全線性關系，完全吻合；如果0<<1，表示模擬值與實測值有一定的線性相關2021.02.05歐陽科創(chuàng)編2021.02.05歐陽科創(chuàng)編 2021.02.05性，其值越大線性相關性越好；如果=0,說明模擬值與實測值無任何線性關系。NashSuttcliffe系數(shù)1的計算公式為：式中：q為實測值；。為模擬值；。，為平均值；n為實測數(shù)據(jù)個數(shù)。當時，J=i；若E為負值，說明模型模擬平均值比直接使用實測平均值的可信度更低。SWAT模型的校準一般分為水平衡和河流流量、泥沙、營養(yǎng)物質(zhì)校準等。在SWAT模型中，產(chǎn)水量是指離開子流域，并在子流域出口匯入河流的水量，是地表徑流、根區(qū)側(cè)向流和地下水徑流減去子流域渠道輸移損失和水池分離后的水量。泥沙產(chǎn)量是指從子流域到達出口的泥沙量，是子流域出口處測量到的泥沙量。（1）徑流驗證徑流校準采用數(shù)字濾波技術(shù)將實測徑流分割為直接徑流和基流兩部分，然后分別對直接徑流和基流進行校準。先對基流進行校準，對于基流校準調(diào)整的參數(shù)包括：地下水再蒸發(fā)系數(shù)R卜一」、土壤蒸發(fā)補償系數(shù)卜”。、土壤可利用水量"Lkc和基流."血系數(shù)。在基流校準達到標準時，對直接徑流進行校準，調(diào)整的參數(shù)為徑流曲線數(shù)八-，其調(diào)整范圍為初始值的?、，根據(jù)不同HRU的八值，計算得其面積平均初始值，最終使總徑流量得到較好的擬合。（2）泥沙負荷驗證2021.02.05歐陽科創(chuàng)編2021