• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    SWAT模型在我國流域水環(huán)境模擬應(yīng)用中的評估驗證過程評價

    2020-11-25 09:10:42秦成新杜鵬飛
    環(huán)境科學(xué)研究 2020年11期
    關(guān)鍵詞:庫中流域文獻(xiàn)

    榮 易, 秦成新, 孫 傅, 杜鵬飛

    清華大學(xué)環(huán)境學(xué)院, 北京 100084

    流域水環(huán)境模型是對水體及其匯水區(qū)受自然過程和人類活動影響呈現(xiàn)出的復(fù)雜響應(yīng)關(guān)系的抽象定量表達(dá),是評估氣候變化、土地管理、人類活動等對流域水環(huán)境影響的重要工具[1]. 近年來,流域水環(huán)境模型在我國得到廣泛應(yīng)用,在非點源污染負(fù)荷估算[2]、污染負(fù)荷源解析[3]、水環(huán)境容量核算[4]、流域規(guī)劃情景預(yù)測[5]等方面發(fā)揮了支撐作用. 由于真實環(huán)境系統(tǒng)的復(fù)雜性,流域水環(huán)境模型通常具有復(fù)雜的非線性結(jié)構(gòu),部分參數(shù)、輸入及邊界條件難以直接獲取,需要通過主觀估計確定,導(dǎo)致模型的模擬結(jié)果存在不確定性. 為了評估和管控模型不確定性對管理決策的影響,模型評估驗證是模型開發(fā)和應(yīng)用不可或缺的工作. 模型評估驗證是指為確定模型及其結(jié)果的質(zhì)量是否適宜支持管理決策而收集和分析信息的過程,具體包括模型遴選、收集輸入數(shù)據(jù)、確定邊界條件、建立模型并識別參數(shù)、模型結(jié)果不確定性分析等步驟.

    SWAT(soil and water assessment tool)模型是美國農(nóng)業(yè)部農(nóng)業(yè)研究所(USDA-ARS)開發(fā)設(shè)計、目前已得到廣泛應(yīng)用的成熟模型,它可在日尺度上連續(xù)模擬土壤蒸散發(fā)、地表徑流、地下水與地表交互及污染物遷移等水文水質(zhì)過程[6]. SWAT模型也是我國流域水環(huán)境模擬的常用模型之一,已有大量應(yīng)用案例[7]. 然而,由于該模型是基于美國相關(guān)研究開發(fā)的,也出現(xiàn)了不適用于我國流域的情況[8]. 因此,將SWAT模型應(yīng)用于我國流域水環(huán)境管理決策時,必須開展完整的模型評估驗證. 但是,利用SWAT模型開展應(yīng)用研究的現(xiàn)有文獻(xiàn)對模型評估驗證的重視程度參差不齊,并且通常重點關(guān)注模型參數(shù)識別、不確定性分析等步驟,對模型空間概化、輸入數(shù)據(jù)獲取等缺少評估[9-11],尚不足以為該模型的評估驗證提供規(guī)范性指導(dǎo).

    該研究針對流域水環(huán)境模型評估驗證過程的質(zhì)量問題,以SWAT模型為例,梳理總結(jié)模型評估驗證過程的評價方法,并利用該方法評價2015—2017年發(fā)表的SWAT模型在我國流域應(yīng)用研究中的模型評估驗證過程,分析模型評估驗證的現(xiàn)狀以及存在的問題和不足,以期為完善SWAT模型的評估驗證方法、促進(jìn)模型的標(biāo)準(zhǔn)化應(yīng)用提出建議.

    1 SWAT模型評估驗證過程的評價方法

    為了評價現(xiàn)有SWAT模型應(yīng)用研究中模型評估驗證過程的質(zhì)量,該研究結(jié)合SWAT模型的建模過程和模型評估驗證的基本步驟,梳理總結(jié)了SWAT模型評估驗證過程的評價方法,具體流程如圖1所示.

    圖1 SWAT模型評估驗證過程的評價方法

    1.1 確定研究或決策目標(biāo)階段

    本階段的評價重點是SWAT模型應(yīng)用中是否清晰識別研究區(qū)域、確認(rèn)研究目標(biāo)并分析模型應(yīng)用的期望產(chǎn)出. 識別研究區(qū)域需關(guān)注研究區(qū)域的地理位置、面積及流域劃分等. 根據(jù)模型應(yīng)用關(guān)注的主要問題,可將研究目標(biāo)分為解釋研究、預(yù)測研究和方法研究. 解釋研究側(cè)重利用模型解釋流域水文水質(zhì)過程的歷史和現(xiàn)狀,預(yù)測研究側(cè)重模擬輸入要素變化下流域的水文水質(zhì)響應(yīng),而方法研究側(cè)重提出模型構(gòu)建及評估驗證新方法. 模型的期望產(chǎn)出應(yīng)匹配研究目標(biāo),選定重點關(guān)注的水文水質(zhì)過程及變量,如河道徑流、泥沙及氮磷含量等.

    1.2 獲取輸入數(shù)據(jù)階段

    本階段的評價重點是模型輸入數(shù)據(jù)的可溯源性、準(zhǔn)確性、精密性和完備性,由此判斷其是否足以支持模型建立與分析. 可溯源性是指輸入數(shù)據(jù)應(yīng)從可靠來源獲取,并詳細(xì)說明來源;準(zhǔn)確性和精密性是指輸入數(shù)據(jù)的產(chǎn)生程序、質(zhì)量控制、誤差水平等應(yīng)符合相關(guān)標(biāo)準(zhǔn)和規(guī)范要求;完備性是指輸入數(shù)據(jù)的數(shù)量和代表性應(yīng)滿足既定研究目標(biāo)的需求.

    1.3 構(gòu)建模型及率定驗證階段

    本階段的評價內(nèi)容包括流域空間劃分、行為設(shè)定、參數(shù)選擇、參數(shù)識別方法、率定期和驗證期長度、模擬效果評價指標(biāo)等. 模型概化時,流域空間劃分方式和人類活動行為等設(shè)定應(yīng)符合流域?qū)嶋H并滿足研究需求. 模型參數(shù)識別時,應(yīng)重點關(guān)注對流域水文水質(zhì)響應(yīng)靈敏的參數(shù),并根據(jù)研究目標(biāo)選擇適宜的參數(shù)識別方法,如基于定向搜索的尋優(yōu)方法和基于采樣的貝葉斯方法. 模型參數(shù)識別和驗證時,還應(yīng)關(guān)注率定期和驗證期觀測數(shù)據(jù)的時間長度以及二者的相對長短,并選擇適當(dāng)?shù)脑u價指標(biāo)評判模型模擬效果.

    1.4 展示模型結(jié)果階段

    本階段的評價重點是模型模擬效果和模擬結(jié)果的不確定性,即評價模擬結(jié)果與觀測數(shù)據(jù)的吻合程度,并表征模擬結(jié)果的不確定性及其對決策的影響.

    2 我國SWAT模型評估驗證過程評價

    該研究評價的我國SWAT模型應(yīng)用研究文獻(xiàn)均來自中國知網(wǎng)(https:www.cnki.net). 數(shù)據(jù)檢索策略:主題為SWAT模型,并限制出版時間為2015—2017年,對文獻(xiàn)類型不做限制. 從檢索結(jié)果中排除綜述性文章、同一研究同時發(fā)表于期刊及作為學(xué)位論文主要內(nèi)容、研究區(qū)域位于我國境外等情況后,共保留428篇文獻(xiàn),建立我國SWAT模型評估驗證過程評價研究文獻(xiàn)庫(簡稱“研究庫”),其中2015年145篇、2016年158篇、2017年125篇.

    2.1 確定研究或決策目標(biāo)

    2.1.1識別研究區(qū)域

    圖2 研究流域的位置和面積分布

    將研究庫中文獻(xiàn)的研究區(qū)域按地理位置和流域面積分類,結(jié)果如圖2所示. 其中,地理位置按流域所在省份歸類至我國地理區(qū)劃,流域面積按全部研究流域面積的四分位數(shù)歸類,分為小型流域(0.62~1 911.5 km2)、中小型流域(1 911.5~5 394.5 km2)、大中型流域(5 394.5~19 032.75 km2)和大型流域(19 032.75~805 000 km2). 從地理位置看,SWAT模型在西北地區(qū)的應(yīng)用研究最多,占到總數(shù)的27.1%,其中,應(yīng)用于大中型及以上流域的文獻(xiàn)達(dá)到75篇,遠(yuǎn)超其他區(qū)域. 這與SWAT模型適合開展大尺度流域模擬,且干旱半干旱的西北地區(qū)更關(guān)注SWAT模型揭示的水資源變化趨勢有關(guān). 在水資源較為缺乏的華北和東北地區(qū),SWAT模型的應(yīng)用也較為廣泛,分別占到文獻(xiàn)總數(shù)的16.1%和15.7%. 從應(yīng)用的具體流域來看,海河流域[12]、渭河流域[13]、長三角水系[14-15]和三峽庫區(qū)[16-17]等是研究熱點.

    2.1.2確認(rèn)研究目標(biāo)

    將研究庫中文獻(xiàn)的研究目標(biāo)分類,結(jié)果如圖3所示. 由圖3可見:解釋研究的文獻(xiàn)有294篇,目標(biāo)包括模擬水文水質(zhì)過程(如蒸散發(fā)量[18]、融雪徑流[19]等)和輔助流域管理(如計算生態(tài)需水量[20]、水環(huán)境容量[21]等);預(yù)測研究的文獻(xiàn)有183篇,目標(biāo)主要是分析氣候、土地利用、管理措施等不同驅(qū)動條件下流域的水文水質(zhì)響應(yīng)[22-23];方法研究的文獻(xiàn)有61篇,目標(biāo)涉及模型改進(jìn)及與其他模型耦合[24]、模型輸入和結(jié)構(gòu)不確定性分析[25]、參數(shù)識別和靈敏度分析方法[26]等. 綜上,我國對SWAT模型的研究仍以應(yīng)用為主,對模型開發(fā)和應(yīng)用中涉及的方法學(xué)研究仍有探索空間.

    圖3 各類研究目標(biāo)對應(yīng)的文獻(xiàn)數(shù)量

    圖4 各類期望產(chǎn)出對應(yīng)的文獻(xiàn)數(shù)量

    2.1.3分析期望產(chǎn)出

    研究庫中文獻(xiàn)模擬或預(yù)測的流域水文水質(zhì)變量主要包括河道徑流量、其他水文過程量(如蒸發(fā)量等)、河道泥沙含量、河道污染物含量等,各類期望產(chǎn)出涉及的文獻(xiàn)數(shù)量如圖4所示. 應(yīng)用SWAT模型研究河道徑流量的文獻(xiàn)最多,占68.0%,其次是研究河道污染物含量,占25.9%. 研究河道泥沙含量和蒸發(fā)量等其他水文過程量的文獻(xiàn)相對較少,這可能與相應(yīng)的監(jiān)測數(shù)據(jù)獲取困難有關(guān). 隨著我國日益重視水生態(tài)保護(hù),SWAT模型也被用于模擬河道重金屬[27]、農(nóng)藥[28]等非常規(guī)污染物,但目前相關(guān)文獻(xiàn)數(shù)量極少,尚有較大發(fā)展空間.

    2.2 獲取輸入數(shù)據(jù)

    SWAT模型必需的輸入數(shù)據(jù)包括數(shù)字高程模型(DEM)、土地利用類型分布、土壤分布與屬性、氣象等數(shù)據(jù). 研究庫中文獻(xiàn)輸入數(shù)據(jù)的可溯源性總體較好,89.5%的文獻(xiàn)提供了所有數(shù)據(jù)的來源信息,主要來源如表1所示. 由表1可見,當(dāng)?shù)叵嚓P(guān)主管部門(如測繪局、氣象局、水文局等)是主要的數(shù)據(jù)來源. 當(dāng)研究區(qū)域較大時,研究者也會使用公開數(shù)據(jù)庫及開源數(shù)據(jù). 但是,不同文獻(xiàn)對開源數(shù)據(jù)的描述方式不一致,有的描述為數(shù)據(jù)名稱(如Landsat TM),有的描述為數(shù)據(jù)來源(如地理空間數(shù)據(jù)云等網(wǎng)絡(luò)數(shù)據(jù)庫),同一數(shù)據(jù)集在不同文獻(xiàn)中可能出現(xiàn)不同的描述. 為了準(zhǔn)確說明模型所用的數(shù)據(jù),應(yīng)同時報告數(shù)據(jù)來源、名稱、版本號等信息,如文獻(xiàn)[29].

    表1 文獻(xiàn)報告的輸入數(shù)據(jù)來源

    由于文獻(xiàn)使用的數(shù)據(jù)通常具有可靠來源,因此一般未報告數(shù)據(jù)的準(zhǔn)確性和精密性. 然而,研究者自行觀測的數(shù)據(jù)(如徑流量、土壤理化性質(zhì)等),則應(yīng)盡可能報告數(shù)據(jù)獲取的相關(guān)方法和依據(jù),以確保數(shù)據(jù)質(zhì)量. 研究庫中12篇文獻(xiàn)采用了自行觀測的數(shù)據(jù),其中6篇文獻(xiàn)報告了數(shù)據(jù)觀測相關(guān)內(nèi)容,其余則受篇幅所限未作報告.

    視研究目標(biāo)不同,SWAT模型可能還需要流域水文水質(zhì)、耕地管理措施等數(shù)據(jù)或信息,需要關(guān)注輸入數(shù)據(jù)的完備性. 例如,在耕地比重較大的流域,耕地管理措施(農(nóng)田耕種行為等)較為重要,研究庫中22.7%的文獻(xiàn)做了相關(guān)設(shè)定.

    2.3 構(gòu)建模型及率定驗證

    2.3.1構(gòu)建模型

    在流域空間劃分方面,研究庫中67.5%的文獻(xiàn)報告了子流域劃分結(jié)果,54.7%的文獻(xiàn)報告了水文響應(yīng)單元(hydrologic response unit, HRU)劃分結(jié)果. 在劃分HRU時,各土地利用類型和土壤類型的劃分閾值越大,計算速度越快,但同時會增加不確定性,因此需要確定合理的閾值. 研究庫中36.9%的文獻(xiàn)報告了劃分閾值,絕大部分不超過20%. 在研究土地利用變化的水文水質(zhì)響應(yīng)時,土地利用類型劃分閾值通常設(shè)為0%,以保留不同土地利用的空間分布信息.

    圖5 參與靈敏度分析以及被判定為靈敏參數(shù)的SWAT模型參數(shù)

    在行為設(shè)定方面,研究庫中文獻(xiàn)設(shè)定最多的行為是農(nóng)業(yè)耕種措施,包括種植、施肥、灌溉、收割等. 這些措施的實施時間通常采用已有文獻(xiàn)的調(diào)研結(jié)果或?qū)嵉卣{(diào)研獲取,在模擬中固定在每年的同一時間. 然而,由于施肥、灌溉等行為受降雨影響較大,固定耕種措施實施時間可能出現(xiàn)降雨時施肥等不合理行為,導(dǎo)致模型模擬與實際出現(xiàn)較大偏差[30]. 為降低這一問題引起的輸出不確定性,可對實施時間開展靈敏度分析,并結(jié)合降雨時間做出適當(dāng)調(diào)整.

    2.3.2模型率定驗證

    研究庫文獻(xiàn)中,SWAT模型涉及的95個參數(shù)參與了參數(shù)靈敏度分析和參數(shù)識別過程,其中水文、泥沙和水質(zhì)相關(guān)的參數(shù)分別為43、10和42個,如圖5所示. 在水文參數(shù)中,SCS徑流曲線數(shù)(CN2)、基流衰退系數(shù)(ALPHA_BF)、土壤蒸發(fā)補償系數(shù)(ESCO)、土壤持水度(SOL_AWC)、主河道水力系數(shù)(CH_K2)、曼寧系數(shù)(CH_N2)等在較多的應(yīng)用案例中表現(xiàn)為靈敏參數(shù),而影響氣象、降雪、匯流和側(cè)向流過程的參數(shù)通常為不靈敏參數(shù). 在泥沙參數(shù)中,影響產(chǎn)沙的通用土壤流失方程耕作管理因子(USLE_C)和水土保持因子(USLE_P),以及影響河道傳輸過程的河道泥沙運移方程的線性參數(shù)(SPCON)和指數(shù)參數(shù)(SPEXP)等,在較多應(yīng)用案例中被判定為靈敏參數(shù). 在水質(zhì)參數(shù)中,氮滲透系數(shù)(NPERCO)、磷滲透系數(shù)(PPERCO)、生物混合系數(shù)(BIOMIX)、磷分配系數(shù)(PHOSKD)等與物質(zhì)遷移相關(guān)的參數(shù)在較多應(yīng)用案例中是靈敏參數(shù),而與初始化、物質(zhì)形態(tài)轉(zhuǎn)化、河道反應(yīng)過程相關(guān)的參數(shù)通常不敏感. 總體而言,SWAT模型水文模擬的控制過程主要包括產(chǎn)流、蒸散發(fā)以及土壤、地下水和河道之間的交互,而水質(zhì)模擬的控制過程主要是污染物的遷移過程. 由于不同流域的水文條件不同,模型的主要控制過程也會出現(xiàn)變化,因此需要在確定參與靈敏度分析的模型參數(shù)前識別流域模擬的主要控制過程,并基于控制過程選擇參數(shù)[31-32].

    研究庫中75.2%的文獻(xiàn)提及了模型參數(shù)靈敏度和不確定性分析過程. 其中,49.3%的文獻(xiàn)使用了SWAT-CUP[33]軟件,9.6%使用了SWAT自帶的靈敏度分析及自動校準(zhǔn)模塊,6.1%采用手動方式調(diào)整參數(shù). SWAT-CUP是專為SWAT模型開發(fā)的,能夠自動完成SWAT模型參數(shù)的靈敏度分析、參數(shù)識別和不確定性分析. 在開展參數(shù)靈敏度分析時,使用SWAT-CUP的多元線性回歸法[34]的文獻(xiàn)最多,而LH-OAT法[35]、Sobol法[36]、Morris法[37]等方法也有文獻(xiàn)使用. 在參數(shù)識別和不確定性分析方法中,SWAT-CUP特有的序列不確定性擬合方法(SUFI-2)[38]使用最為廣泛,少量文獻(xiàn)使用了粒子群算法(PSO)[39]、復(fù)合型混合演化算法(SCE-UA)[40]、廣義似然不確定性估計方法(GLUE)[41]等. 由于SWAT-CUP軟件的廣泛應(yīng)用,我國SWAT模型應(yīng)用研究中采用的靈敏度和不確定性分析方法大部分為該軟件支持的方法.

    在率定和驗證期長度方面,研究庫中85.7%的文獻(xiàn)報告了率定期和驗證期觀測數(shù)據(jù)的時間長度. 現(xiàn)有文獻(xiàn)中率定期數(shù)據(jù)時間長度的中位數(shù)與驗證期接近,均為5年左右,但率定期數(shù)據(jù)時間長度總體上長于驗證期,50%的文獻(xiàn)中率定期數(shù)據(jù)時間長度為驗證期的1~2倍,25%的文獻(xiàn)中率定期數(shù)據(jù)時間長度為驗證期的2倍以上. 研究庫中也有17篇文獻(xiàn)在率定期和驗證期使用的數(shù)據(jù)時間長度均不足1年,這可能會制約模型在不同水文條件下的模擬能力.

    研究庫中文獻(xiàn)使用的模型模擬效果評價指標(biāo)主要為統(tǒng)計學(xué)指標(biāo),包括相關(guān)系數(shù)或決定系數(shù)(correlation coefficientcoefficient of determination),相對誤差或百分比偏差(relative errorpercent bias)以及納什效率系數(shù)(Nash-Sutcliffe efficiency coefficient)3種,它們在現(xiàn)有文獻(xiàn)中分別被使用了345、207和380次. 除此之外,也有少量文獻(xiàn)使用均方根誤差(root mean square error)、均方根誤差與觀測值標(biāo)準(zhǔn)偏差比(RMSE-observations standard deviation ratio)、克林-古普塔效率系數(shù)(Kling-Gupta efficiency)等評價指標(biāo). 國外SWAT模型研究中使用的流量歷時曲線坡度偏差、峰值時間及強度偏差等動力學(xué)評價指標(biāo),在研究庫文獻(xiàn)中尚未見采用. 已有研究表明,使用動力學(xué)評價指標(biāo)可更好地識別模型參數(shù)[42],與統(tǒng)計學(xué)指標(biāo)組合使用可以提高模型對極低流量的模擬精度[43],并且有助于發(fā)現(xiàn)模型參數(shù)或結(jié)構(gòu)問題[44].

    2.4 展示模型結(jié)果

    研究庫中文獻(xiàn)報告的決定系數(shù)和納什效率系數(shù)值較高,所有決定系數(shù)均大于0.6,絕大部分納什效率系數(shù)均大于0.5. 兩個評價指標(biāo)的中位數(shù)均超過0.8,高于擬合程度“相當(dāng)好”的標(biāo)準(zhǔn)0.75[45],說明大部分文獻(xiàn)中SWAT模型模擬結(jié)果與觀測數(shù)據(jù)的擬合效果較好. 但也有2篇文獻(xiàn)的納什效率系數(shù)小于0.5,表明這些文獻(xiàn)建立的模型對流域的表征能力弱,模擬結(jié)果不可靠.

    研究庫中文獻(xiàn)對模型模擬結(jié)果不確定性的討論相對較少,研究者通常只使用單一的模型“最優(yōu)”參數(shù)開展預(yù)測模擬,因而可能導(dǎo)致預(yù)測結(jié)果出現(xiàn)偏差. 部分文獻(xiàn)[46]通過統(tǒng)計所有接受的模擬結(jié)果得到模擬結(jié)果的95%置信區(qū)間,并分析了該置信區(qū)間涵蓋的監(jiān)測數(shù)據(jù)數(shù)量和模擬結(jié)果標(biāo)準(zhǔn)差,以評價模擬結(jié)果的覆蓋度和分散性. 國際上一些文獻(xiàn)將模型模擬結(jié)果的不確定性轉(zhuǎn)換為風(fēng)險損失費用,以更直觀地表現(xiàn)不同決策所需承受的風(fēng)險[47],但研究庫文獻(xiàn)中尚未見類似研究.

    3 SWAT模型評估驗證問題識別與建議

    3.1 模型模擬過程和參數(shù)識別的選擇

    SWAT模型的徑流模擬是泥沙模擬的基礎(chǔ),而水質(zhì)模擬又與泥沙模擬有關(guān). 因此,開展泥沙模擬應(yīng)首先完成徑流模擬,而開展水質(zhì)模擬則需同時完成徑流和泥沙模擬[33]. 然而,部分文獻(xiàn)在應(yīng)用SWAT模型時并未遵守該流程. 例如,在模擬研究區(qū)污染負(fù)荷時,未率定模型對流量和泥沙的模擬效果,而直接率定模型對總氮、總磷的模擬效果[48];或僅率定了流量和水質(zhì)的模擬效果,而未率定泥沙的模擬效果,也未將泥沙相關(guān)參數(shù)納入靈敏度分析和參數(shù)識別的范圍[49]. 相較于流量和水質(zhì)數(shù)據(jù),泥沙數(shù)據(jù)較難獲取,因此模型率定存在客觀限制. 但在識別水質(zhì)過程時,仍應(yīng)將泥沙相關(guān)參數(shù)納入?yún)?shù)識別范圍.

    從模型參數(shù)看,SWAT模型中使用的城市、植物數(shù)據(jù)庫以及徑流曲線數(shù)(CN2)、曼寧系數(shù)(OV_N)等均是基于美國環(huán)境特征設(shè)定的,直接移用到我國流域可能存在適用性的問題. 已有研究[50]表明,SWAT模型默認(rèn)的CN2對照表與我國實際情況存在較大偏差. 再例如,研究庫中文獻(xiàn)開展參數(shù)靈敏度分析涉及的95個參數(shù)中,僅包含植物數(shù)據(jù)庫中最大潛在葉面積指數(shù)(BLAI)這一個參數(shù),但是在其他被忽視的植物參數(shù)中,太陽輻射利用率(BIO_E)、作物最佳葉面積曲線第一參數(shù)點(FRGRW1)等參數(shù)可能對徑流模擬是靈敏的[51]. 因此,在我國開展SWAT模型評估驗證時,建議盡可能使用本地化數(shù)據(jù)庫或?qū)⑾嚓P(guān)模型參數(shù)納入靈敏度分析,以降低模型模擬的不確定性.

    3.2 參數(shù)取值超出合理范圍

    如前文所述,SUFI-2為SWAT模型最常用的參數(shù)識別算法,該算法識別得到的參數(shù)最優(yōu)值可能落在參數(shù)先驗分布范圍外[52]. 這一特點可糾正不合理的參數(shù)先驗分布范圍,但是當(dāng)參數(shù)的取值范圍具有明確物理意義時,參數(shù)取值落在設(shè)定的范圍外則出現(xiàn)錯誤. 常見的參數(shù)取值錯誤主要有3種形式:①參數(shù)取值過小,地下水相關(guān)參數(shù)(ALPHA_BF、GE_REVAP、GWQMN、RCHDG_DP)取值為負(fù)數(shù)的情況在文獻(xiàn)中頻繁出現(xiàn)[53-60],SCS徑流曲線數(shù)(CN2)[61-63]、冠層最大蓄水量(CANMX)[64]、土壤持水量(SOL_AWC)[53]、主河道曼寧系數(shù)(CH_N2)[65]等參數(shù)取值過小情況也時有發(fā)生. ②參數(shù)取值過大,主要體現(xiàn)在與河道(CH_N2)[66]、匯流(CH_N1)[66]、土壤(ESCO、SOL_AWC)[67-70]等相關(guān)的參數(shù)上. ③參數(shù)取值的相對大小不合理,例如一些文獻(xiàn)在識別降雪參數(shù)時得到了最小溶雪因子(SMFMN)取值大于最大溶雪因子(SMFMX)的結(jié)果[71-75]. 雖然SWAT模型會將這些不合理取值調(diào)整為合理數(shù)值,但這也將導(dǎo)致參數(shù)報告值與實際率定值不一致,容易引起研究者和模型用戶對參數(shù)取值的誤解. 因此,當(dāng)出現(xiàn)模型率定出現(xiàn)不合理參數(shù)取值時,研究者應(yīng)對這些參數(shù)取值進(jìn)行分析和調(diào)整.

    3.3 流域水文過程表征不合理

    研究庫中部分文獻(xiàn)應(yīng)用SWAT模型的模擬效果不佳,在模型率定期表現(xiàn)尚可但在驗證期表現(xiàn)較差[60]. 分析這些文獻(xiàn)的參數(shù)取值和時間序列模擬結(jié)果可以發(fā)現(xiàn),模擬效果較差不只是因為產(chǎn)流模式與模型預(yù)設(shè)不匹配、驗證期水文環(huán)境變化等[60,76],還與率定期模型水文過程表征不準(zhǔn)確有關(guān).

    模型對水文過程的表征與模型參數(shù)取值密切相關(guān). 例如,在北洛河流域的研究案例中,模型率定得到的參數(shù)ALPHA_BF取值為0.64,CN2為35,而驗證期徑流模擬的納什效率系數(shù)低至-32.125[76]. CN2值較低導(dǎo)致降雨幾乎不產(chǎn)生地表徑流而直接進(jìn)入土壤和地下水,而ALPHA_BF值過高又導(dǎo)致進(jìn)入淺層地下水的水量會立刻進(jìn)入河道,從而降低基流的流量. 因此,模型對北洛河水文過程不準(zhǔn)確的表征,造成了徑流量較低時模擬偏差較大的問題. 建議在模型評估驗證時,既注重提高模型率定期的模擬效果,也注重對水文過程的合理表征.

    4 結(jié)論與展望

    4.1 結(jié)論

    a) 我國現(xiàn)有的SWAT模型應(yīng)用研究總體比較重視模型評估驗證過程,部分文獻(xiàn)報告了輸入數(shù)據(jù)來源和質(zhì)量、模型概化、參數(shù)靈敏度分析、參數(shù)識別、不確定性分析等過程的技術(shù)細(xì)節(jié),模型模擬效果較好.

    b) SWAT模型評估驗證過程的主要薄弱環(huán)節(jié)是參數(shù)識別,具體表現(xiàn)為關(guān)鍵參數(shù)選取與研究區(qū)域?qū)嶋H特征存在偏差、參數(shù)識別結(jié)果不合理等,進(jìn)而影響模型模擬效果和水文水質(zhì)過程表征.

    c) 國際上SWAT模型評估驗證的一些新技術(shù),例如使用動力學(xué)指標(biāo)評價模擬效果、分時段分析模擬結(jié)果誤差、在情景模擬評價中討論模型輸出不確定性等,在我國的研究文獻(xiàn)中尚未見采用.

    4.2 展望

    a) 在SWAT模型應(yīng)用研究中應(yīng)重點關(guān)注模型參數(shù)靈敏度分析和參數(shù)識別,結(jié)合流域?qū)嶋H情況識別關(guān)鍵模擬過程、篩選潛在靈敏參數(shù)、選取合理參數(shù)取值,確保過程模擬結(jié)果與流域?qū)嶋H相符.

    b) 建立SWAT模型本地化應(yīng)用數(shù)據(jù)庫,更新模型涉及的非本地化設(shè)計的數(shù)據(jù),如植物數(shù)據(jù)庫、城市數(shù)據(jù)庫、CN值對照表等.

    c) SWAT模型應(yīng)用于流域水環(huán)境管理決策時,應(yīng)完整報告模型評估驗證過程相關(guān)信息,并開展不確定性分析,降低決策風(fēng)險.

    d) 加強SWAT模型評估驗證方法研究,引入和研發(fā)評估驗證新技術(shù),提高模擬結(jié)果的準(zhǔn)確度與可靠性,為模型應(yīng)用于流域水環(huán)境管理決策提供支撐.

    猜你喜歡
    庫中流域文獻(xiàn)
    壓油溝小流域
    動物城堡
    動物城堡
    Hostile takeovers in China and Japan
    速讀·下旬(2021年11期)2021-10-12 01:10:43
    堡子溝流域綜合治理
    羅堰小流域
    Cultural and Religious Context of the Two Ancient Egyptian Stelae An Opening Paragraph
    大東方(2019年12期)2019-10-20 13:12:49
    打造智慧流域的思路及構(gòu)想——以討賴河流域為例
    The Application of the Situational Teaching Method in English Classroom Teaching at Vocational Colleges
    智能盤庫在自動化立體庫中的探索和應(yīng)用
    国产真实乱freesex| 午夜激情福利司机影院| 国产一区二区激情短视频| 亚洲精品一卡2卡三卡4卡5卡| 国产午夜精品论理片| 亚洲av成人一区二区三| 免费一级毛片在线播放高清视频| 亚洲国产欧美网| 最近视频中文字幕2019在线8| 全区人妻精品视频| 岛国在线观看网站| 欧美激情久久久久久爽电影| 超碰成人久久| 欧美另类亚洲清纯唯美| 欧美日韩福利视频一区二区| 国产精品久久视频播放| 国产精品久久视频播放| 嫁个100分男人电影在线观看| 成人手机av| 一本一本综合久久| 国产真实乱freesex| 中文字幕高清在线视频| 色尼玛亚洲综合影院| 久久精品亚洲精品国产色婷小说| 日本五十路高清| 久久精品影院6| 一个人观看的视频www高清免费观看 | 久久精品夜夜夜夜夜久久蜜豆 | 黑人欧美特级aaaaaa片| 成人特级黄色片久久久久久久| 一a级毛片在线观看| 久久人妻福利社区极品人妻图片| 久久久国产成人精品二区| 精品国产乱子伦一区二区三区| 在线观看免费视频日本深夜| av欧美777| 男女之事视频高清在线观看| 精品欧美一区二区三区在线| 亚洲午夜精品一区,二区,三区| 在线视频色国产色| 成人午夜高清在线视频| 久9热在线精品视频| 天堂影院成人在线观看| 全区人妻精品视频| 哪里可以看免费的av片| 非洲黑人性xxxx精品又粗又长| 啦啦啦韩国在线观看视频| 国产99白浆流出| 丰满人妻一区二区三区视频av | 色噜噜av男人的天堂激情| 在线观看美女被高潮喷水网站 | 亚洲中文字幕一区二区三区有码在线看 | 国产精品,欧美在线| 亚洲五月天丁香| √禁漫天堂资源中文www| 国产精品一区二区三区四区免费观看 | 国产亚洲精品一区二区www| 伊人久久大香线蕉亚洲五| 变态另类丝袜制服| 免费在线观看亚洲国产| 一本综合久久免费| 欧美在线黄色| 欧美zozozo另类| 特大巨黑吊av在线直播| 变态另类成人亚洲欧美熟女| 男女床上黄色一级片免费看| 99精品久久久久人妻精品| 啦啦啦观看免费观看视频高清| 妹子高潮喷水视频| 国产三级黄色录像| 99久久精品国产亚洲精品| 黑人操中国人逼视频| 午夜视频精品福利| 成人av一区二区三区在线看| 欧美性猛交黑人性爽| 免费看十八禁软件| 91大片在线观看| 欧美三级亚洲精品| 久久久水蜜桃国产精品网| 别揉我奶头~嗯~啊~动态视频| 又爽又黄无遮挡网站| 久久草成人影院| 午夜福利18| 老司机靠b影院| 美女高潮喷水抽搐中文字幕| 一进一出抽搐动态| 一区二区三区国产精品乱码| 亚洲精品一区av在线观看| 俄罗斯特黄特色一大片| 成人国语在线视频| 久久国产精品影院| 亚洲精品粉嫩美女一区| 欧美3d第一页| 久久久久国产精品人妻aⅴ院| 久久亚洲精品不卡| 久久性视频一级片| 色综合婷婷激情| av中文乱码字幕在线| 听说在线观看完整版免费高清| 97超级碰碰碰精品色视频在线观看| 亚洲国产日韩欧美精品在线观看 | 欧美日本视频| 欧美在线黄色| 国产精品久久视频播放| 中文字幕久久专区| 伊人久久大香线蕉亚洲五| 最近视频中文字幕2019在线8| 午夜老司机福利片| 日韩欧美国产在线观看| 久久天躁狠狠躁夜夜2o2o| 国产精品久久久人人做人人爽| 中文字幕人成人乱码亚洲影| 成人18禁高潮啪啪吃奶动态图| 久久 成人 亚洲| 又紧又爽又黄一区二区| 好看av亚洲va欧美ⅴa在| 天天躁夜夜躁狠狠躁躁| 国产欧美日韩一区二区三| 亚洲精华国产精华精| 一进一出抽搐动态| 亚洲av片天天在线观看| 国产伦在线观看视频一区| 白带黄色成豆腐渣| 亚洲真实伦在线观看| 午夜免费成人在线视频| 国产蜜桃级精品一区二区三区| 亚洲av熟女| 1024视频免费在线观看| 国产区一区二久久| 午夜日韩欧美国产| 熟女电影av网| 日本 av在线| 久久久久久久精品吃奶| 又爽又黄无遮挡网站| 国产激情偷乱视频一区二区| 波多野结衣高清无吗| 嫩草影院精品99| e午夜精品久久久久久久| 日日摸夜夜添夜夜添小说| 大型黄色视频在线免费观看| 国产精品日韩av在线免费观看| 一级作爱视频免费观看| 欧美久久黑人一区二区| 成人特级黄色片久久久久久久| 亚洲,欧美精品.| 久久久久国内视频| 一本精品99久久精品77| 久久久久免费精品人妻一区二区| 好男人电影高清在线观看| 床上黄色一级片| av视频在线观看入口| 1024视频免费在线观看| 免费看日本二区| 波多野结衣高清无吗| 天天躁夜夜躁狠狠躁躁| 1024视频免费在线观看| 免费看日本二区| 看黄色毛片网站| 999久久久精品免费观看国产| 五月伊人婷婷丁香| 欧美日韩亚洲综合一区二区三区_| 一级作爱视频免费观看| 淫妇啪啪啪对白视频| 欧美成人午夜精品| 国产精品香港三级国产av潘金莲| 俄罗斯特黄特色一大片| 久久久精品国产亚洲av高清涩受| 日韩大尺度精品在线看网址| 国产精品久久久久久久电影 | 欧美日韩亚洲国产一区二区在线观看| 身体一侧抽搐| 日韩欧美免费精品| 91老司机精品| 午夜福利18| 亚洲精品久久国产高清桃花| 特级一级黄色大片| 丰满人妻一区二区三区视频av | 中文字幕熟女人妻在线| 校园春色视频在线观看| 亚洲av成人不卡在线观看播放网| 亚洲五月婷婷丁香| 亚洲成av人片在线播放无| 久久精品aⅴ一区二区三区四区| 1024视频免费在线观看| 特大巨黑吊av在线直播| 天堂影院成人在线观看| 在线观看舔阴道视频| 此物有八面人人有两片| 18禁美女被吸乳视频| www.精华液| 又黄又爽又免费观看的视频| 身体一侧抽搐| 丰满人妻一区二区三区视频av | 无限看片的www在线观看| 97碰自拍视频| a级毛片a级免费在线| 国产成人一区二区三区免费视频网站| 中文字幕久久专区| 亚洲一区二区三区色噜噜| 久久精品国产清高在天天线| 色播亚洲综合网| 亚洲avbb在线观看| 白带黄色成豆腐渣| 国产av一区二区精品久久| 一个人观看的视频www高清免费观看 | 狠狠狠狠99中文字幕| 在线永久观看黄色视频| 国内毛片毛片毛片毛片毛片| 美女 人体艺术 gogo| 怎么达到女性高潮| 亚洲国产高清在线一区二区三| 久久精品国产亚洲av香蕉五月| 欧美成人性av电影在线观看| 最新在线观看一区二区三区| 欧美日韩国产亚洲二区| 国产三级在线视频| 俄罗斯特黄特色一大片| 国产亚洲欧美98| 国产熟女xx| 午夜福利18| 日日干狠狠操夜夜爽| 精品欧美一区二区三区在线| 日韩精品青青久久久久久| 精品久久久久久久末码| 国产成+人综合+亚洲专区| 亚洲av美国av| 熟女少妇亚洲综合色aaa.| 大型黄色视频在线免费观看| 一边摸一边抽搐一进一小说| 欧美绝顶高潮抽搐喷水| 久久久久久久午夜电影| 男女那种视频在线观看| 久久这里只有精品中国| 中出人妻视频一区二区| 国产精品久久久久久亚洲av鲁大| tocl精华| 国产乱人伦免费视频| 亚洲中文av在线| 国产激情久久老熟女| 亚洲av成人不卡在线观看播放网| 亚洲成av人片免费观看| 曰老女人黄片| 99国产综合亚洲精品| 亚洲在线自拍视频| 久久亚洲精品不卡| 国产单亲对白刺激| 老熟妇仑乱视频hdxx| 亚洲国产精品成人综合色| 99精品在免费线老司机午夜| 91麻豆精品激情在线观看国产| 97超级碰碰碰精品色视频在线观看| 亚洲国产高清在线一区二区三| 男女之事视频高清在线观看| netflix在线观看网站| 欧美黄色片欧美黄色片| 99热这里只有精品一区 | 不卡一级毛片| 老鸭窝网址在线观看| 久久精品国产综合久久久| 亚洲成人国产一区在线观看| 一卡2卡三卡四卡精品乱码亚洲| 又粗又爽又猛毛片免费看| 十八禁人妻一区二区| 天天一区二区日本电影三级| 亚洲人成77777在线视频| 制服人妻中文乱码| 国产不卡一卡二| 亚洲成人久久爱视频| 欧美成人午夜精品| 国产精品 国内视频| 国产成+人综合+亚洲专区| 亚洲熟妇熟女久久| 国产成人av激情在线播放| 美女午夜性视频免费| 狂野欧美白嫩少妇大欣赏| 小说图片视频综合网站| 中亚洲国语对白在线视频| 国产伦在线观看视频一区| 国产黄a三级三级三级人| 日韩欧美国产在线观看| 久久婷婷人人爽人人干人人爱| 老司机午夜十八禁免费视频| www.自偷自拍.com| 18禁黄网站禁片午夜丰满| 日本黄色视频三级网站网址| 亚洲精华国产精华精| 欧美 亚洲 国产 日韩一| 久久久水蜜桃国产精品网| 伦理电影免费视频| 午夜a级毛片| 首页视频小说图片口味搜索| 亚洲乱码一区二区免费版| 免费在线观看黄色视频的| 精品少妇一区二区三区视频日本电影| 亚洲中文日韩欧美视频| 波多野结衣高清作品| 黄色丝袜av网址大全| www日本黄色视频网| 日本撒尿小便嘘嘘汇集6| 不卡一级毛片| 久久久久九九精品影院| 久久香蕉精品热| 亚洲aⅴ乱码一区二区在线播放 | 亚洲成av人片免费观看| 全区人妻精品视频| 黄片大片在线免费观看| 日本黄色视频三级网站网址| 男人的好看免费观看在线视频 | 毛片女人毛片| 中文亚洲av片在线观看爽| 69av精品久久久久久| 搡老熟女国产l中国老女人| 哪里可以看免费的av片| 一级毛片女人18水好多| 一进一出好大好爽视频| 亚洲欧洲精品一区二区精品久久久| 亚洲国产看品久久| 欧美成人免费av一区二区三区| 亚洲av成人不卡在线观看播放网| 狂野欧美白嫩少妇大欣赏| 黄色 视频免费看| 这个男人来自地球电影免费观看| 国产精品一区二区精品视频观看| 人成视频在线观看免费观看| 国产乱人伦免费视频| 成人午夜高清在线视频| 欧美不卡视频在线免费观看 | 亚洲人成网站在线播放欧美日韩| 日韩大码丰满熟妇| 精品少妇一区二区三区视频日本电影| 亚洲精品久久国产高清桃花| 午夜两性在线视频| 免费看十八禁软件| 夜夜看夜夜爽夜夜摸| 亚洲精华国产精华精| 曰老女人黄片| 免费人成视频x8x8入口观看| 国产精品国产高清国产av| 久久中文字幕人妻熟女| 国产亚洲精品av在线| 国产精品九九99| 99热这里只有是精品50| 长腿黑丝高跟| 高潮久久久久久久久久久不卡| 在线免费观看的www视频| 国内久久婷婷六月综合欲色啪| 国产一区在线观看成人免费| 成人亚洲精品av一区二区| 精华霜和精华液先用哪个| 成年免费大片在线观看| 人妻夜夜爽99麻豆av| 精品国产乱码久久久久久男人| 欧美日韩中文字幕国产精品一区二区三区| 久久久久久久精品吃奶| 亚洲精品av麻豆狂野| 伊人久久大香线蕉亚洲五| 在线观看美女被高潮喷水网站 | 精品久久久久久久末码| 久久精品国产综合久久久| 黄片小视频在线播放| 老司机午夜十八禁免费视频| 精品不卡国产一区二区三区| 黄色视频,在线免费观看| 久久久久精品国产欧美久久久| 免费在线观看影片大全网站| 舔av片在线| 久久午夜综合久久蜜桃| 伊人久久大香线蕉亚洲五| 母亲3免费完整高清在线观看| 免费一级毛片在线播放高清视频| av欧美777| 狂野欧美白嫩少妇大欣赏| av免费在线观看网站| 精品久久久久久久毛片微露脸| 亚洲欧洲精品一区二区精品久久久| 亚洲真实伦在线观看| 在线观看免费日韩欧美大片| 国产熟女xx| 深夜精品福利| 老司机福利观看| 嫩草影视91久久| 丝袜人妻中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 别揉我奶头~嗯~啊~动态视频| 两个人的视频大全免费| 午夜老司机福利片| 日韩成人在线观看一区二区三区| 欧美一级毛片孕妇| 亚洲片人在线观看| 美女扒开内裤让男人捅视频| 男女做爰动态图高潮gif福利片| 亚洲第一电影网av| 免费在线观看完整版高清| 淫妇啪啪啪对白视频| АⅤ资源中文在线天堂| 亚洲精品av麻豆狂野| 欧美性猛交黑人性爽| 波多野结衣巨乳人妻| 亚洲狠狠婷婷综合久久图片| 国产精品乱码一区二三区的特点| 久久精品国产亚洲av高清一级| 亚洲最大成人中文| 18美女黄网站色大片免费观看| 午夜两性在线视频| 99精品久久久久人妻精品| 久久精品人妻少妇| 桃色一区二区三区在线观看| av在线播放免费不卡| www.精华液| 欧美午夜高清在线| 国产精品电影一区二区三区| 一个人免费在线观看的高清视频| 日本熟妇午夜| 母亲3免费完整高清在线观看| 国产精品亚洲一级av第二区| 麻豆成人午夜福利视频| 人人妻人人澡欧美一区二区| 亚洲精品国产一区二区精华液| 国产精品av视频在线免费观看| 天堂影院成人在线观看| 变态另类丝袜制服| 一本综合久久免费| 亚洲五月天丁香| 动漫黄色视频在线观看| 无人区码免费观看不卡| 97人妻精品一区二区三区麻豆| 91麻豆av在线| 18禁国产床啪视频网站| 日韩大尺度精品在线看网址| 欧美性猛交╳xxx乱大交人| 久久精品亚洲精品国产色婷小说| 亚洲国产欧美一区二区综合| 亚洲一区二区三区色噜噜| 国产成人aa在线观看| 欧美日本亚洲视频在线播放| 亚洲av电影不卡..在线观看| 99久久国产精品久久久| 嫩草影视91久久| 久久久国产成人精品二区| 巨乳人妻的诱惑在线观看| 亚洲国产中文字幕在线视频| 国产精品,欧美在线| 99精品在免费线老司机午夜| 香蕉丝袜av| 亚洲av五月六月丁香网| 男人的好看免费观看在线视频 | 老司机午夜福利在线观看视频| 日本 av在线| 久久热在线av| 久久精品综合一区二区三区| 国产精品乱码一区二三区的特点| 俺也久久电影网| 免费在线观看影片大全网站| 久久久久久亚洲精品国产蜜桃av| 久久精品国产亚洲av香蕉五月| 亚洲在线自拍视频| 成在线人永久免费视频| 国产精品国产高清国产av| 俺也久久电影网| 欧美一区二区精品小视频在线| 成年免费大片在线观看| 国产蜜桃级精品一区二区三区| 怎么达到女性高潮| 免费看美女性在线毛片视频| 国产精品98久久久久久宅男小说| 毛片女人毛片| 国产三级在线视频| 国产片内射在线| 十八禁网站免费在线| 久久精品国产综合久久久| 免费人成视频x8x8入口观看| 三级国产精品欧美在线观看 | 亚洲精品美女久久久久99蜜臀| 国产激情偷乱视频一区二区| 亚洲七黄色美女视频| 久热爱精品视频在线9| 久久香蕉激情| 精华霜和精华液先用哪个| 国产aⅴ精品一区二区三区波| 国产一区二区激情短视频| 欧美绝顶高潮抽搐喷水| 久久久久久久精品吃奶| 三级国产精品欧美在线观看 | 一本大道久久a久久精品| a级毛片在线看网站| 不卡一级毛片| 欧美中文日本在线观看视频| 国产欧美日韩一区二区三| 黄色毛片三级朝国网站| 国产成人精品久久二区二区91| 欧美三级亚洲精品| 久久久久久亚洲精品国产蜜桃av| 91麻豆精品激情在线观看国产| 99国产精品99久久久久| а√天堂www在线а√下载| 全区人妻精品视频| 国产视频内射| 国产精品一及| 看黄色毛片网站| 亚洲aⅴ乱码一区二区在线播放 | 可以免费在线观看a视频的电影网站| 丰满的人妻完整版| 亚洲国产欧美网| 色精品久久人妻99蜜桃| 97碰自拍视频| 丰满人妻熟妇乱又伦精品不卡| 哪里可以看免费的av片| www日本在线高清视频| 国产精华一区二区三区| 女人高潮潮喷娇喘18禁视频| 天堂影院成人在线观看| 成在线人永久免费视频| 欧美成人午夜精品| 桃色一区二区三区在线观看| 宅男免费午夜| 一本综合久久免费| 久久性视频一级片| 日韩成人在线观看一区二区三区| 国产精品一区二区免费欧美| 制服丝袜大香蕉在线| 三级男女做爰猛烈吃奶摸视频| 99久久国产精品久久久| 少妇被粗大的猛进出69影院| 窝窝影院91人妻| 69av精品久久久久久| 日韩欧美 国产精品| 色av中文字幕| 女生性感内裤真人,穿戴方法视频| 超碰成人久久| 国产熟女午夜一区二区三区| 日韩成人在线观看一区二区三区| 久久久久久久久免费视频了| 国内毛片毛片毛片毛片毛片| 亚洲精品在线美女| 国产v大片淫在线免费观看| www.999成人在线观看| 亚洲va日本ⅴa欧美va伊人久久| 18禁裸乳无遮挡免费网站照片| www.自偷自拍.com| 亚洲欧洲精品一区二区精品久久久| 亚洲人成77777在线视频| 小说图片视频综合网站| 国产亚洲欧美98| 日韩欧美三级三区| 日韩国内少妇激情av| 桃红色精品国产亚洲av| 日韩成人在线观看一区二区三区| 淫妇啪啪啪对白视频| 亚洲七黄色美女视频| 一区二区三区高清视频在线| www.精华液| 淫秽高清视频在线观看| 国产精品98久久久久久宅男小说| 亚洲人成伊人成综合网2020| 亚洲成人国产一区在线观看| 又黄又粗又硬又大视频| 两性夫妻黄色片| 欧美在线黄色| 亚洲av中文字字幕乱码综合| 免费看十八禁软件| 久久人人精品亚洲av| 伦理电影免费视频| 精品国产乱码久久久久久男人| 国产一级毛片七仙女欲春2| 国产精品亚洲av一区麻豆| 88av欧美| 色老头精品视频在线观看| 女警被强在线播放| 香蕉av资源在线| 熟妇人妻久久中文字幕3abv| 精品欧美国产一区二区三| 男女下面进入的视频免费午夜| 成人精品一区二区免费| 日韩大码丰满熟妇| 丁香欧美五月| 国产1区2区3区精品| 在线观看午夜福利视频| 天天躁狠狠躁夜夜躁狠狠躁| 中文字幕久久专区| 日本 欧美在线| 欧美日本亚洲视频在线播放| 日韩大尺度精品在线看网址| 一个人观看的视频www高清免费观看 | 黄色成人免费大全| 久99久视频精品免费| 久久精品aⅴ一区二区三区四区| 亚洲国产欧美网| 青草久久国产| 国产精品香港三级国产av潘金莲| 日韩欧美在线乱码| 亚洲成av人片免费观看| 久久亚洲精品不卡| 人妻丰满熟妇av一区二区三区| 亚洲av日韩精品久久久久久密| 欧美乱码精品一区二区三区| 一夜夜www| 国产精品亚洲一级av第二区| 后天国语完整版免费观看| 亚洲精品在线观看二区| 亚洲av电影不卡..在线观看| 黄片小视频在线播放| 久久 成人 亚洲| 一夜夜www| 欧美成人午夜精品| 女警被强在线播放| 欧美在线黄色| 黑人欧美特级aaaaaa片| 国产一区二区三区在线臀色熟女| 亚洲欧美日韩东京热| 黄片大片在线免费观看| 国产三级中文精品| 成年女人毛片免费观看观看9| 午夜免费观看网址| 波多野结衣巨乳人妻| 大型av网站在线播放|