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

    智慧流域水文預(yù)報(bào)技術(shù)研究進(jìn)展與開發(fā)前景

    2023-08-26 19:30:40王俊程海云郭生練張俊崔震
    人民長(zhǎng)江 2023年8期
    關(guān)鍵詞:水文洪水流域

    王俊 程海云 郭生練 張俊 崔震

    摘要:

    當(dāng)前處于傳統(tǒng)水利快速邁向智慧水利的高質(zhì)量發(fā)展新階段,開展智慧流域水文模擬和洪水預(yù)報(bào)系統(tǒng)研發(fā),是建設(shè)數(shù)字孿生流域的重要組成部分,是強(qiáng)化“四預(yù)”措施和推進(jìn)智慧水利建設(shè)的必然要求。綜述了流域水文模擬預(yù)報(bào)研究進(jìn)展,剖析中國(guó)洪水預(yù)報(bào)業(yè)務(wù)現(xiàn)狀和不足,從數(shù)據(jù)底板、規(guī)律機(jī)理、預(yù)報(bào)模型、方法技術(shù)和業(yè)務(wù)平臺(tái)等方面對(duì)水文預(yù)報(bào)技術(shù)進(jìn)行了探討展望。研究結(jié)果表明,智慧流域水文預(yù)報(bào)的關(guān)鍵技術(shù)與發(fā)展前景為:① 多源異構(gòu)數(shù)據(jù)匯集,夯實(shí)流域數(shù)據(jù)底板;② 加強(qiáng)河庫(kù)系統(tǒng)水汽循環(huán)研究,揭示成因機(jī)理;③ 聚焦強(qiáng)人類活動(dòng)影響特征,升級(jí)專業(yè)模型;④ 深化預(yù)報(bào)新技術(shù)研究,探索基于智能預(yù)報(bào)調(diào)度一體化業(yè)務(wù)體系;⑤ 打造通用化業(yè)務(wù)化平臺(tái),支撐多要素全周期“四預(yù)”功能。研究成果可為提高流域洪水預(yù)報(bào)水平,強(qiáng)化“四預(yù)”措施,建成數(shù)字孿生流域提供理論基礎(chǔ)和技術(shù)支撐。

    關(guān) 鍵 詞:

    智慧流域; 數(shù)字孿生; 水文預(yù)報(bào); 智能模擬; 概率預(yù)報(bào)

    中圖法分類號(hào): TV124

    文獻(xiàn)標(biāo)志碼: A

    DOI:10.16232/j.cnki.1001-4179.2023.08.001

    0 引 言

    國(guó)家“十四五”規(guī)劃綱要明確提出“構(gòu)建智慧水利體系,以流域?yàn)閱卧嵘闇y(cè)報(bào)和智能調(diào)度能力”。水利部高度重視智慧水利建設(shè),將其作為推動(dòng)新階段水利高質(zhì)量發(fā)展的六條實(shí)施路徑之一[1]。2021年12月,水利部召開推進(jìn)數(shù)字孿生流域建設(shè)工作會(huì)議,要求大力推進(jìn)數(shù)字孿生流域建設(shè)。2022年3月,水利部印發(fā)《數(shù)字孿生流域建設(shè)技術(shù)大綱(試行)》[2],指出數(shù)字孿生流域是智慧水利的核心與關(guān)鍵,其建設(shè)總體目標(biāo)之一是建成大江大河大湖及主要支流、重點(diǎn)流域和重點(diǎn)區(qū)域的數(shù)字孿生流域,實(shí)現(xiàn)與物理流域同步仿真運(yùn)行、虛實(shí)交互、迭代優(yōu)化,支撐“四預(yù)”(預(yù)報(bào)、預(yù)警、預(yù)演、預(yù)案)功能實(shí)現(xiàn)和“2+N”(流域防洪、水資源管理與調(diào)配+N項(xiàng)業(yè)務(wù)應(yīng)用)智能應(yīng)用運(yùn)行,加快構(gòu)建智慧水利體系,提升水利決策與管理的科學(xué)化、精準(zhǔn)化、高效化能力和水平,為新階段水利高質(zhì)量發(fā)展提供有力支撐和強(qiáng)力驅(qū)動(dòng)??梢?,在新時(shí)期新形勢(shì)下,開展智慧流域水文模擬和洪水預(yù)報(bào)系統(tǒng)研發(fā),是建設(shè)數(shù)字孿生流域的重要組成部分,是強(qiáng)化“四預(yù)”措施和推進(jìn)智慧水利建設(shè)的必然要求。

    中國(guó)幅員遼闊,氣候復(fù)雜,自然地理多樣,水旱災(zāi)害頻繁,嚴(yán)重影響著流域經(jīng)濟(jì)社會(huì)發(fā)展[3]。為了防治水旱災(zāi)害,各流域內(nèi)建設(shè)了大規(guī)模的水利工程群。這些水工程群的聯(lián)合調(diào)度運(yùn)用事關(guān)流域防洪、能源、航運(yùn)、供水、生態(tài)安全,對(duì)流域綜合管理和經(jīng)濟(jì)社會(huì)發(fā)展意義重大。然而,水工程群的建成使得流域內(nèi)水文循環(huán)、暴雨洪水產(chǎn)匯流特性受影響的程度有了新的變化,增加了洪水預(yù)報(bào)的難度和不確定性。本文在綜述流域水文模擬和洪水預(yù)報(bào)業(yè)務(wù)的研究進(jìn)展基礎(chǔ)上,從數(shù)據(jù)底板、規(guī)律機(jī)理、預(yù)報(bào)模型、方法技術(shù)和業(yè)務(wù)平臺(tái)等方面,探討流域水文預(yù)報(bào)技術(shù)研發(fā)的前景,為提高流域洪水預(yù)報(bào)精度,強(qiáng)化“四預(yù)”措施,建成數(shù)字孿生流域提供理論基礎(chǔ)和技術(shù)支撐[4]。

    1 流域水文模擬預(yù)報(bào)研究進(jìn)展

    1.1 流域水文模型

    縱觀國(guó)內(nèi)外洪水預(yù)報(bào)模型和技術(shù)的發(fā)展歷程,大體可分為產(chǎn)匯流機(jī)理研究、概念性水文模型、系統(tǒng)(黑箱)模型、分布式水文模型等發(fā)展階段。20世紀(jì)30~50年代,洪水預(yù)報(bào)的發(fā)展處于萌芽和起步階段,國(guó)內(nèi)外學(xué)者的研究多聚焦在水文循環(huán)各環(huán)節(jié)的演化機(jī)理方面,代表性成果有謝爾曼經(jīng)驗(yàn)單位線、霍頓下滲曲線、馬斯京根法、彭曼蒸發(fā)公式、Nash瞬時(shí)單位線等;20世紀(jì)60~80年代,開展了大量的野外山坡水文機(jī)理實(shí)驗(yàn),揭示了非飽和側(cè)向流、壤中流和飽和地表徑流多種產(chǎn)流機(jī)制的存在[5]。這些研究結(jié)果促進(jìn)了水文模型的發(fā)展。20世紀(jì)60年代末,世界上第一個(gè)水文模型(Stanford-IV)模型在美國(guó)誕生,這標(biāo)志著水文學(xué)家開始結(jié)合系統(tǒng)理論的思想,將流域水文循環(huán)各要素視為一個(gè)有機(jī)的整體進(jìn)行研究,具有里程碑意義;20世紀(jì)60~90年代,由于系統(tǒng)理論應(yīng)用的逐步深入以及計(jì)算機(jī)技術(shù)的快速發(fā)展,流域水文模型的研發(fā)和應(yīng)用進(jìn)入蓬勃發(fā)展期,代表性模型有美國(guó)Sacramento和API模型、中國(guó)新安江模型[6]、日本TANK模型、瑞典HBV模型、丹麥NAM模型、英國(guó)TOPMODEL模型、法國(guó)GR4J模型、愛爾蘭LPM模型、意大利CLS模型等。

    概念性模型和系統(tǒng)(黑箱)模型一般都是集總式的,為了表征流域在輸入、下墊面條件及模型參數(shù)的時(shí)空差異性,英國(guó)水文學(xué)家Freeze和Harlan[7]提出了分布式水文模型的構(gòu)造框架。20世紀(jì)90年代以后,由于地面站點(diǎn)建設(shè)不斷推進(jìn),尤其是3S技術(shù)的不斷革新與推廣和數(shù)值高程模型的建立,分布式水文模型快速發(fā)展[8],代表性成果包括歐洲SHE模型、意大利TOPKAPI模型、美國(guó)SWAT和VIC模型、中國(guó)GBHM模型[9]等。分布式水文模型耦合了流域物理特征、邊界條件以及水文過(guò)程的時(shí)空變化,可幫助人們深入認(rèn)識(shí)流域水文過(guò)程機(jī)理[10],在無(wú)資料地區(qū)水文預(yù)報(bào)和水資源量評(píng)估等方面已開展了廣泛應(yīng)用,但因?yàn)橄聣|面信息支撐不足等原因,其在洪水預(yù)報(bào)領(lǐng)域還處于試驗(yàn)應(yīng)用階段[11]。

    經(jīng)過(guò)90余年的發(fā)展,流域水文模型已臻成熟,近20 a來(lái),水文模擬技術(shù)的發(fā)展聚焦在與其他學(xué)科或技術(shù)的交叉融合上,如深化定量降水預(yù)測(cè)預(yù)報(bào)技術(shù)、利用遙感技術(shù)提取下墊面信息、改進(jìn)搜索技術(shù)率定模型參數(shù)等。實(shí)況輸入、降雨預(yù)報(bào)、模型、參數(shù)是制約洪水預(yù)報(bào)精度和預(yù)見期的關(guān)鍵所在。

    1.2 人工智能預(yù)報(bào)

    近年來(lái),人工智能和數(shù)據(jù)挖掘技術(shù)被引入至水文預(yù)報(bào)領(lǐng)域,機(jī)器學(xué)習(xí)是實(shí)現(xiàn)人工智能的基本途徑[12],因其無(wú)需考慮水量平衡原理等物理機(jī)制、僅從數(shù)據(jù)中尋求輸入輸出間的關(guān)系,因此可以更好地?cái)M合非線性的徑流序列。其中人工神經(jīng)網(wǎng)絡(luò)(ANN)是最具代表性的機(jī)器學(xué)習(xí)算法之一,根據(jù)網(wǎng)絡(luò)架構(gòu)是否有回饋?lái)?xiàng)可劃分為靜態(tài)神經(jīng)網(wǎng)絡(luò)和動(dòng)態(tài)神經(jīng)網(wǎng)絡(luò),前者的網(wǎng)絡(luò)架構(gòu)沒有回饋?lái)?xiàng),對(duì)資料結(jié)構(gòu)的長(zhǎng)期映射關(guān)系具有較好的解析效果,而后者具有回饋?lái)?xiàng),對(duì)資料結(jié)構(gòu)的短期映射關(guān)系具有較好的解析效果。在水文預(yù)報(bào)領(lǐng)域應(yīng)用較廣的靜態(tài)類神經(jīng)網(wǎng)絡(luò)包括反向傳播神經(jīng)網(wǎng)絡(luò)、徑向基神經(jīng)網(wǎng)絡(luò)、自組織映射網(wǎng)絡(luò)、時(shí)滯神經(jīng)網(wǎng)絡(luò)、適應(yīng)性網(wǎng)絡(luò)模糊推論系統(tǒng)等;應(yīng)用較廣的動(dòng)態(tài)類神經(jīng)網(wǎng)絡(luò)包括Elman回饋式神經(jīng)網(wǎng)絡(luò)、即時(shí)回饋式神經(jīng)網(wǎng)絡(luò)、非線性自回歸外因輸入模式等[13]。除ANN外,支持向量機(jī)(SVM)、決策樹(DT)、多層感知機(jī)(MLP)和隨機(jī)森林(RF)等許多機(jī)器學(xué)習(xí)算法同樣具有處理非線性和不確定性的強(qiáng)大能力,在水文模擬預(yù)報(bào)和水庫(kù)調(diào)度領(lǐng)域得到廣泛的應(yīng)用。張珂等[13]對(duì)比分析了DT、MLP、RF和SVM四種模型在陜西省3個(gè)半干旱半濕潤(rùn)區(qū)典型流域洪水預(yù)報(bào)的適用性,發(fā)現(xiàn)4種模型在半濕潤(rùn)區(qū)典型流域預(yù)報(bào)精度較高,適用性較優(yōu),在半干旱流域精度偏低。同時(shí),隨著預(yù)見期增加,SVM整體表現(xiàn)穩(wěn)定,RF和DT預(yù)報(bào)精度下降緩慢,MLP穩(wěn)定性較差。胡義明等[14]首先采用置換準(zhǔn)確度重要性度量法篩選關(guān)鍵預(yù)報(bào)因子,其次采用隨機(jī)搜索技術(shù)結(jié)合交叉驗(yàn)證方式確定模型參數(shù)和網(wǎng)絡(luò)結(jié)構(gòu),最后分析了AdaBoost模型(AdB)、RF和SVM的中長(zhǎng)期徑流預(yù)報(bào)精度,結(jié)果發(fā)現(xiàn)AdB模型適用于淮河流域的中長(zhǎng)期徑流預(yù)報(bào),預(yù)報(bào)性能優(yōu)于RF和SVM模型。

    盡管眾多機(jī)器學(xué)習(xí)模型算法廣泛應(yīng)用于時(shí)間序列預(yù)測(cè)問題,但仍存在泛化性能弱、預(yù)見期較短等問題,為解決該難題,出現(xiàn)了較多基于ANN的深度學(xué)習(xí)方法,如循環(huán)神經(jīng)網(wǎng)絡(luò)、卷積神經(jīng)網(wǎng)絡(luò)等。長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)(LSTM)則通過(guò)在神經(jīng)網(wǎng)絡(luò)的隱藏層中引入存儲(chǔ)單元,來(lái)選擇記憶當(dāng)前信息或遺忘過(guò)去記憶信息(如降雨-徑流映射關(guān)系),增強(qiáng)了神經(jīng)網(wǎng)絡(luò)的長(zhǎng)期記憶能力[15-18]。Hu等[19]將LSTM模型應(yīng)用于汾河流域洪水事件模擬研究中,發(fā)現(xiàn)LSTM模型的模擬性能優(yōu)于ANN模型,為洪水預(yù)報(bào)提出了新的數(shù)據(jù)驅(qū)動(dòng)方法。徐源浩等[20]探討了LSTM模型在汾河流域洪水預(yù)報(bào)中的應(yīng)用,并研究了超參數(shù)對(duì)預(yù)報(bào)精度的影響。結(jié)果表明,LSTM模型在6 h以上預(yù)見期預(yù)報(bào)效果相對(duì)較差,預(yù)報(bào)精度隨神經(jīng)元數(shù)量和訓(xùn)練次數(shù)的增加呈上升趨勢(shì),當(dāng)達(dá)到一定值時(shí),預(yù)報(bào)精度趨于穩(wěn)定。Xu等[21]采用粒子群優(yōu)化算法優(yōu)選LSTM模型的超參數(shù),建立PSO-LSTM深度學(xué)習(xí)模型,并在汾河和漯河流域進(jìn)行實(shí)驗(yàn),結(jié)果表明,PSO-LSTM模型學(xué)習(xí)數(shù)據(jù)特征的能力優(yōu)于ANN和LSTM等模型,提高了短期洪水預(yù)報(bào)的準(zhǔn)確性。李步等[22]采用主成分分析法提取氣象要素的空間特征,并與LSTM神經(jīng)網(wǎng)絡(luò)相結(jié)合,構(gòu)建了PCA-LSTM深度學(xué)習(xí)模型,以研究黃河源區(qū)氣象要素空間特征對(duì)預(yù)報(bào)精度的影響,結(jié)果表明所提模型具有較高的預(yù)報(bào)精度。張海榮等[23]采用偏互信息法篩選關(guān)鍵遙相關(guān)因子,將其作為L(zhǎng)STM神經(jīng)網(wǎng)絡(luò)輸入,在宜昌站進(jìn)行月徑流預(yù)報(bào),結(jié)果表明篩選因子使得LSTM神經(jīng)網(wǎng)絡(luò)更具有物理意義,可有效提高預(yù)報(bào)精度。

    隨著人工智能技術(shù)快速發(fā)展,采用編碼-解碼(ED)結(jié)構(gòu)的LSTM深度學(xué)習(xí)模型,得到了多時(shí)段和高精度的洪水過(guò)程預(yù)報(bào)[24-29]。Kao等[24]、Han等[25]、林康聆等[27]在石門水庫(kù)流域、美國(guó)Russian河流域、建溪流域分別構(gòu)建了LSTM-ED模型,結(jié)果表明LSTM-ED模型預(yù)報(bào)性能優(yōu)于LSTM模型,但當(dāng)預(yù)見期大于流域最大匯流時(shí)間時(shí),模型預(yù)報(bào)性能逐漸變差。Cui等[26]提出一種耦合XAJ模型預(yù)報(bào)流量的外源輸入編碼-解碼結(jié)構(gòu),并耦合至LSTM神經(jīng)網(wǎng)絡(luò)(LSTM-EDE),在陸水和建溪流域進(jìn)行3~12 h預(yù)見期的洪水預(yù)報(bào)。結(jié)果表明,所提模型可以克服遞歸編碼-解碼結(jié)構(gòu)的訓(xùn)練過(guò)程與驗(yàn)證過(guò)程不一致的問題,提高了洪水過(guò)程預(yù)報(bào)精度。Girihagama等[28]在加拿大10個(gè)流域構(gòu)建了基于注意力機(jī)制和編碼-解碼結(jié)構(gòu)的LSTM模型用于1~5 d預(yù)見期的徑流預(yù)報(bào),并與標(biāo)準(zhǔn)的LSTM-ED模型進(jìn)行對(duì)比,結(jié)果表明所提模型具有更優(yōu)的預(yù)報(bào)性能。崔震等[29]將ED結(jié)構(gòu)耦合至LSTM神經(jīng)網(wǎng)絡(luò),構(gòu)建了LSTM-ED深度學(xué)習(xí)模型,并采用貝葉斯預(yù)報(bào)處理器量化預(yù)報(bào)不確定性,用于三峽水庫(kù)1~7 d預(yù)見期的入庫(kù)洪水預(yù)報(bào),結(jié)果表明LSTM-ED模型的納什效率系數(shù)在0.92以上,同時(shí)考慮預(yù)報(bào)降雨信息可提高概率預(yù)報(bào)性能。

    綜上可知,深度學(xué)習(xí)算法的持續(xù)發(fā)展促使了其在預(yù)報(bào)業(yè)務(wù)中的需求逐漸增長(zhǎng),特征可視化或反映降雨徑流響應(yīng)規(guī)律的可解釋性深度學(xué)習(xí)、能夠量化預(yù)報(bào)不確定性的概率深度學(xué)習(xí)、用于缺資料流域的區(qū)域深度學(xué)習(xí)和多任務(wù)集成深度學(xué)習(xí)等研究均得到水文學(xué)者的廣泛關(guān)注,將是人工智能預(yù)報(bào)極具潛力的方向。然而,智能洪水預(yù)報(bào)仍處于研究探索階段,在國(guó)內(nèi)尚未投入生產(chǎn)業(yè)務(wù),制約了洪水預(yù)報(bào)客觀化的程度,亟待開展智能洪水預(yù)報(bào)的業(yè)務(wù)示范應(yīng)用。

    1.3 洪水概率預(yù)報(bào)

    傳統(tǒng)的水文預(yù)報(bào)是依據(jù)實(shí)測(cè)雨量進(jìn)行預(yù)報(bào),預(yù)見期較短,并僅給用戶提供一個(gè)確定的預(yù)報(bào)值,難以滿足實(shí)際工作的需要。隨著數(shù)值天氣預(yù)報(bào)技術(shù)的發(fā)展,數(shù)值降水預(yù)報(bào)為延長(zhǎng)預(yù)見期、提高洪水預(yù)報(bào)精度提供了可能。然而,數(shù)值降水預(yù)報(bào)存在誤差,降水作為水文模型最重要的輸入,它直接影響預(yù)報(bào)精度。此外,水文預(yù)報(bào)的不確定性還主要源于模型結(jié)構(gòu)、模型參數(shù)的誤差。隨著對(duì)于預(yù)報(bào)理論認(rèn)識(shí)的逐漸深入,人們發(fā)現(xiàn)基于特定水文模型和參數(shù)獲得的水文預(yù)報(bào)結(jié)果不可避免存在著預(yù)報(bào)誤差,可通過(guò)概率預(yù)報(bào)方法,即使用概率分布函數(shù)定量描述預(yù)報(bào)結(jié)果的不確定性,從而讓決策者更好地評(píng)估風(fēng)險(xiǎn),因此概率預(yù)報(bào)已經(jīng)成為當(dāng)前水文預(yù)報(bào)領(lǐng)域的研究熱點(diǎn)之一[30]。

    縱觀國(guó)內(nèi)外學(xué)者相關(guān)研究,目前洪水概率預(yù)報(bào)主要分為兩類途徑[31]:一類是集合概率分析途徑,另一類是總誤差分析途徑。集合概率預(yù)報(bào)是當(dāng)前刻畫降水預(yù)報(bào)和模型結(jié)構(gòu)不確定性的重要手段,然而集合預(yù)報(bào)結(jié)果對(duì)于不確定性的描述常常是不準(zhǔn)確的,必須通過(guò)統(tǒng)計(jì)后處理方法予以校正[32]。目前廣泛使用的統(tǒng)計(jì)后處理方法有兩大類,一類是集合模型輸出統(tǒng)計(jì)法(Ensemble Model Output Statistics,EMOS)[33],另一類是貝葉斯模型平均法(Bayesian Model Averaging,BMA) [34]。集合概率預(yù)報(bào)能甄別預(yù)報(bào)的不確定性來(lái)源,但計(jì)算量較大,且對(duì)輸入要求較高,在滿足實(shí)時(shí)作業(yè)預(yù)報(bào)的要求方面較困難。

    總誤差分析途徑則是基于確定性預(yù)報(bào)結(jié)果/預(yù)報(bào)誤差與預(yù)報(bào)量的關(guān)系進(jìn)行不確定性分析的概率預(yù)報(bào)方法,該法的輸入是確定性預(yù)報(bào)過(guò)程,重點(diǎn)是預(yù)報(bào)量的后驗(yàn)密度函數(shù)推求,其缺點(diǎn)是對(duì)歷史樣本依賴性強(qiáng),但計(jì)算快捷,使用方便。Krzysztofowicz[35]提出的貝葉斯概率預(yù)報(bào)系統(tǒng)(BFS)是該途徑最具代表性的成果,是公認(rèn)的通過(guò)確定性水文模型進(jìn)行概率預(yù)報(bào)解決水文預(yù)報(bào)不確定性的理論框架。針對(duì)BFS需進(jìn)行線性-正態(tài)假設(shè)的不足,劉章君等[36]利用Copula函數(shù)推導(dǎo)了流量先驗(yàn)分布及似然函數(shù)的解析表達(dá)式,構(gòu)建了Copula-BFS模型,不僅簡(jiǎn)化了計(jì)算難度,還可提高概率預(yù)報(bào)精度。近年來(lái)一些學(xué)者聚焦通過(guò)分析不同預(yù)見期、流量級(jí)下預(yù)報(bào)誤差的分布規(guī)律來(lái)量化預(yù)報(bào)的不確定性[37]。

    總體而言,概率預(yù)報(bào)計(jì)算的模型和方法已較為成熟,但在與調(diào)度風(fēng)險(xiǎn)、調(diào)度決策耦合方面的研究和應(yīng)用還比較缺乏,在當(dāng)前國(guó)內(nèi)以基于確定性理念的預(yù)報(bào)調(diào)度業(yè)務(wù)為主的背景下,概率預(yù)報(bào)仍未被生產(chǎn)決策部門接納,未真正實(shí)現(xiàn)落地推廣。

    2 中國(guó)洪水預(yù)報(bào)業(yè)務(wù)進(jìn)展

    2.1 洪水預(yù)報(bào)系統(tǒng)

    洪水預(yù)報(bào)系統(tǒng)是水文預(yù)報(bào)領(lǐng)域具有代表性的先進(jìn)技術(shù),作為實(shí)時(shí)作業(yè)預(yù)報(bào)的業(yè)務(wù)工具,洪水預(yù)報(bào)系統(tǒng)通過(guò)水雨情監(jiān)視預(yù)警、防洪形勢(shì)分析以及預(yù)報(bào)調(diào)度計(jì)算等,可為預(yù)報(bào)調(diào)度提供重要技術(shù)支撐,并極大地提高預(yù)報(bào)員的工作效率[38-39]。洪水預(yù)報(bào)系統(tǒng)在國(guó)內(nèi)外研究開發(fā)已有30~40 a的歷史,主要經(jīng)歷了4個(gè)發(fā)展階段:第一階段是聯(lián)機(jī)預(yù)報(bào)作業(yè)階段,其主要特點(diǎn)是集水情信息采集、傳輸、處理和洪水預(yù)報(bào)計(jì)算為一體,以便快速完成洪水預(yù)報(bào)作業(yè),如英國(guó)河川徑流預(yù)報(bào)系統(tǒng)、美國(guó)NOAA洪水預(yù)報(bào)系統(tǒng)(NWSRFS);第二階段是實(shí)時(shí)預(yù)報(bào)校正階段,在洪水預(yù)報(bào)系統(tǒng)中引入現(xiàn)代控制理論,實(shí)現(xiàn)實(shí)時(shí)信息和預(yù)報(bào)結(jié)果的實(shí)時(shí)校正,如長(zhǎng)江中下游洪水預(yù)報(bào)調(diào)度系統(tǒng)[40];第三階段是交互式洪水預(yù)報(bào)階段,利用圖形交互處理技術(shù)對(duì)洪水預(yù)報(bào)中間環(huán)節(jié)進(jìn)行人工干預(yù),充分利用專家、預(yù)報(bào)員的知識(shí)和經(jīng)驗(yàn),以有效地提高洪水預(yù)報(bào)水平,如國(guó)家防汛抗旱指揮系統(tǒng)一期工程——中國(guó)洪水預(yù)報(bào)系統(tǒng)[41]、長(zhǎng)江水利委員會(huì)推出的通用型水文預(yù)報(bào)平臺(tái)[42];目前,隨著計(jì)算機(jī)技術(shù)和網(wǎng)絡(luò)技術(shù)的進(jìn)步,Web服務(wù)模式下洪水預(yù)報(bào)系統(tǒng)的開發(fā)成為可能,洪水預(yù)報(bào)系統(tǒng)真正從模型開發(fā)走向應(yīng)用服務(wù),在業(yè)務(wù)功能處理上逐步實(shí)現(xiàn)了預(yù)報(bào)調(diào)度一體化、會(huì)商決策信息化和科學(xué)化支撐等功能,基于Web服務(wù)模式下構(gòu)建預(yù)報(bào)調(diào)度一體化系統(tǒng),正向洪水預(yù)報(bào)系統(tǒng)發(fā)展的第四階段邁進(jìn),如國(guó)家防汛抗旱指揮系統(tǒng)二期工程——長(zhǎng)江防洪預(yù)報(bào)調(diào)度系統(tǒng)[43]。

    當(dāng)前,從事洪水預(yù)報(bào)業(yè)務(wù)的各級(jí)水情部門和發(fā)電企業(yè)均開發(fā)了各自的洪水預(yù)報(bào)系統(tǒng),水平總體處在系統(tǒng)發(fā)展的第三或者第四階段,但大多數(shù)系統(tǒng)的定制化程度高,且自動(dòng)化、國(guó)產(chǎn)化、標(biāo)準(zhǔn)化程度不高,導(dǎo)致移植性較差,業(yè)務(wù)推廣不足。

    2.2 洪水預(yù)報(bào)業(yè)務(wù)應(yīng)用現(xiàn)狀及不足

    新中國(guó)成立以來(lái),中國(guó)水文情報(bào)預(yù)報(bào)工作從無(wú)到有、由點(diǎn)到面迅速發(fā)展。傳統(tǒng)的洪水預(yù)報(bào)方法是從工程水文學(xué)計(jì)算方法直接移植的,相關(guān)圖、謝爾曼單位線、馬斯京根法等是先期預(yù)報(bào)方法的主流,隨著計(jì)算機(jī)技術(shù)、信息處理技術(shù)、作業(yè)預(yù)報(bào)技術(shù)手段不斷創(chuàng)新發(fā)展,水情自動(dòng)測(cè)報(bào)系統(tǒng)、水文模擬技術(shù)、洪水預(yù)報(bào)模型和方法、洪水預(yù)報(bào)系統(tǒng)建設(shè)取得了長(zhǎng)足進(jìn)步[44]。但隨著經(jīng)濟(jì)社會(huì)的快速發(fā)展,一方面,人類活動(dòng)深刻影響著流域水文循環(huán),預(yù)報(bào)的不確定性影響因素更加復(fù)雜,應(yīng)推進(jìn)量化預(yù)報(bào)不確定性和調(diào)度風(fēng)險(xiǎn)的業(yè)務(wù)實(shí)踐,以及時(shí)做出相應(yīng)的優(yōu)化決策和防范措施。另一方面,跨區(qū)域、大規(guī)模水利工程群的多目標(biāo)聯(lián)合調(diào)度運(yùn)用是流域綜合管理的重要手段[45],相應(yīng)地,高精度的中、長(zhǎng)預(yù)見期徑流預(yù)報(bào)成為指導(dǎo)水資源管理的重要參考。同時(shí),傳統(tǒng)水利正快速邁向智慧水利階段,數(shù)字孿生流域建設(shè)亟需提升支撐能力、創(chuàng)新服務(wù)手段。這些都需要以水文預(yù)測(cè)預(yù)報(bào)為前提保障,因而對(duì)水文預(yù)報(bào)的發(fā)展提出了新的需求:① 提高預(yù)報(bào)精度、延長(zhǎng)預(yù)見期;② 客觀量化預(yù)報(bào)不確定性,更好支撐風(fēng)險(xiǎn)調(diào)度決策;③ 提高智慧化水平,解放生產(chǎn)力。面對(duì)新時(shí)期新形勢(shì)下的新需求,就流域洪水預(yù)報(bào)現(xiàn)狀而言,在算據(jù)獲取、洪水預(yù)報(bào)模型適應(yīng)性、新技術(shù)落地推廣及業(yè)務(wù)系統(tǒng)通用化等方面仍面臨較大挑戰(zhàn),是亟待解決的短板[46]。

    2.2.1 算據(jù)獲取能力不足

    豐富的數(shù)據(jù)底板是準(zhǔn)確、客觀描述流域水文循環(huán)的基礎(chǔ)。當(dāng)前,流域站網(wǎng)密度較優(yōu)的區(qū)域多集中在發(fā)達(dá)地區(qū)或平原區(qū),西部偏遠(yuǎn)山區(qū)、河源、中小河流等區(qū)域仍顯稀少,且站點(diǎn)維護(hù)不足,而高時(shí)空分辨率的降雨分布、下墊面信息、河道地形、大范圍土壤墑情、高精度傾斜攝影數(shù)據(jù)短缺,以及空天地一體化監(jiān)測(cè)的衛(wèi)星遙感、雷達(dá)監(jiān)測(cè)反演數(shù)據(jù)利用度不高,一定程度上制約了多時(shí)空尺度的水文模擬水平提升。上述均是流域數(shù)據(jù)底板構(gòu)建的短板,使得流域“算據(jù)”獲取能力不足。

    2.2.2 預(yù)報(bào)模型的適應(yīng)性面臨挑戰(zhàn)

    氣候變化導(dǎo)致陸氣互饋影響更復(fù)雜,局地暴雨特性改變,氣象水文預(yù)報(bào)耦合不緊密,降水預(yù)報(bào)難度加大,難以滿足洪水預(yù)報(bào)對(duì)定量化、定點(diǎn)化、定時(shí)化的精準(zhǔn)要求;掣肘于站網(wǎng)密度不足和預(yù)報(bào)模型未深入應(yīng)用等原因,未控區(qū)間產(chǎn)匯流計(jì)算并不精準(zhǔn);現(xiàn)有預(yù)報(bào)方案中產(chǎn)匯流模型的率定檢驗(yàn)多滿足于平均狀況下的目標(biāo)函數(shù)最優(yōu),對(duì)參數(shù)的異參同效等不夠重視;面臨全生命周期的水文預(yù)報(bào),即由主要為汛期洪水預(yù)報(bào)向全年的水文水資源預(yù)報(bào)轉(zhuǎn)變,需要考慮模型結(jié)構(gòu)和參數(shù)對(duì)季節(jié)性變化規(guī)律的響應(yīng)調(diào)整;人類活動(dòng)影響下的下墊面變化導(dǎo)致流域產(chǎn)匯流特性和機(jī)理改變,原有洪水預(yù)報(bào)模型的適應(yīng)性面臨挑戰(zhàn),既要考慮水庫(kù)調(diào)度對(duì)洪水預(yù)報(bào)互饋影響下的預(yù)報(bào)調(diào)度一體化,又要建立水文水力學(xué)相結(jié)合的河道(庫(kù)區(qū))匯流方法,以滿足變化條件下大流域、長(zhǎng)河系、多阻斷的流域水工程聯(lián)合調(diào)度和流域綜合管理的現(xiàn)實(shí)需要。

    2.2.3 新技術(shù)有待業(yè)務(wù)化落地和推廣

    面臨精準(zhǔn)預(yù)報(bào)、智能調(diào)度、風(fēng)險(xiǎn)決策等新需求,分布式模擬、智能模擬、集合預(yù)報(bào)、概率預(yù)報(bào)等新技術(shù)已成為行業(yè)發(fā)展的前沿方向,相關(guān)研究業(yè)已如火如荼地開展,但目前還多處于研究探索階段,尚未構(gòu)建成熟的技術(shù)方法體系,科研與應(yīng)用存在脫節(jié),與生產(chǎn)部門的實(shí)際需求和基礎(chǔ)結(jié)合不夠,未針對(duì)性開展業(yè)務(wù)化設(shè)計(jì)和典型示范,導(dǎo)致仍未實(shí)現(xiàn)落地推廣。例如,多模型業(yè)務(wù)預(yù)報(bào)仍然是點(diǎn)綴,成功應(yīng)用于業(yè)務(wù)預(yù)報(bào)的模型單一,開展集合預(yù)報(bào)困難;分布式模型因未能體現(xiàn)出相比集總式模型的模擬預(yù)報(bào)精度優(yōu)勢(shì),尚未在業(yè)務(wù)應(yīng)用展開;機(jī)器學(xué)習(xí)還是以神經(jīng)網(wǎng)絡(luò)等為初級(jí)表現(xiàn);特別要指出的是,引入風(fēng)險(xiǎn)概念的洪水概率預(yù)報(bào),尚未與調(diào)度決策耦合,仍未被生產(chǎn)決策部門接納應(yīng)用。

    2.2.4 業(yè)務(wù)系統(tǒng)的通用化和產(chǎn)品化程度有待提升

    縱觀科研和生產(chǎn)實(shí)踐現(xiàn)狀,關(guān)于水文建模和實(shí)時(shí)作業(yè)預(yù)報(bào),國(guó)內(nèi)已研發(fā)了種類齊全、數(shù)目繁多的業(yè)務(wù)系統(tǒng),但這些系統(tǒng)多基于特定模型或小流域構(gòu)建,定制化程度高,研發(fā)與生產(chǎn)應(yīng)用脫節(jié),普遍存在“四預(yù)”功能覆蓋不全,業(yè)務(wù)流程設(shè)計(jì)不完善,國(guó)產(chǎn)化、客觀化、自動(dòng)化、智能化程度不高,標(biāo)準(zhǔn)化、通用化、產(chǎn)品化不夠等問題,導(dǎo)致共享性、移植性差,推廣程度低、應(yīng)用場(chǎng)景有限[47]。

    3 流域水文預(yù)報(bào)關(guān)鍵技術(shù)與研發(fā)前景

    3.1 多源異構(gòu)數(shù)據(jù)匯集,夯實(shí)流域數(shù)據(jù)底板

    3.1.1 構(gòu)建基于空天地一體化的流域立體監(jiān)測(cè)體系

    數(shù)據(jù)源的多寡、質(zhì)量高低是制約洪水預(yù)報(bào)精度的主要因素之一。為夯實(shí)流域數(shù)據(jù)底板,需納入多維度信息監(jiān)測(cè)裝備,構(gòu)建基于空天地一體化的流域立體監(jiān)測(cè)體系,在無(wú)資料和少資料地區(qū)加密建設(shè)站點(diǎn),并提升維護(hù)能力,利用衛(wèi)星、多普勒雷達(dá)等遙感手段反演面雨量,采用高分衛(wèi)星遙感影像解譯、激光雷達(dá)、無(wú)人機(jī)航拍及高精度傾斜攝影獲取河道地形與大斷面、地形地貌和土地利用等資料。

    3.1.2 加強(qiáng)多主體信息匯集及共享平臺(tái)建設(shè)

    不同主體的監(jiān)測(cè)數(shù)據(jù)各有側(cè)重,匯集更全面、更多維的數(shù)據(jù)信息,有利于深入挖掘數(shù)據(jù)隱含的知識(shí)規(guī)則,提高數(shù)據(jù)使用效力。因此,流域水文預(yù)報(bào)需要加強(qiáng)與氣象、國(guó)土、環(huán)保、發(fā)電企業(yè)等不同行業(yè)、部門的信息共享,通過(guò)構(gòu)建流域?qū)用娑嘀黧w聯(lián)合的信息共享平臺(tái)和數(shù)據(jù)存儲(chǔ)中心,實(shí)現(xiàn)多源異構(gòu)數(shù)據(jù)的共享和存儲(chǔ)。

    3.1.3 開展多源數(shù)據(jù)融合和同化技術(shù)研究

    數(shù)據(jù)融合和同化是一個(gè)多源信息集成的過(guò)程。通過(guò)數(shù)據(jù)同化算法不斷融合時(shí)空上離散分布的不同來(lái)源和不同分辨率的直接或間接監(jiān)測(cè)信息,以動(dòng)態(tài)改善模型擬合軌跡、提高計(jì)算精度,是多源數(shù)據(jù)信息對(duì)水文預(yù)報(bào)模型預(yù)測(cè)精度的重要貢獻(xiàn)。因此,需持續(xù)開展多源數(shù)據(jù)融合和同化技術(shù)研究,使流域數(shù)據(jù)底板初始場(chǎng)數(shù)據(jù)更有效、更智慧。

    3.2 加強(qiáng)河庫(kù)系統(tǒng)水汽循環(huán)研究,揭示成因機(jī)理

    3.2.1 變化環(huán)境下流域暴雨特性演變成因分析

    大規(guī)模水庫(kù)群的建設(shè)運(yùn)行,導(dǎo)致流域陸面植被等下墊面條件發(fā)生顯著變化,引起蒸散發(fā)、土壤含水量等特性指標(biāo)發(fā)生變異,而水庫(kù)群超大庫(kù)面的形成和年內(nèi)、年際變化,又進(jìn)一步影響了陸氣能量交換和局地水汽循環(huán)規(guī)律,導(dǎo)致流域內(nèi)局部區(qū)域的暴雨過(guò)程、落區(qū)和強(qiáng)度發(fā)生改變,亟待研究水庫(kù)庫(kù)區(qū)暴雨特性變化及陸氣互饋影響機(jī)理。

    3.2.2 河庫(kù)系統(tǒng)影響下洪水傳播機(jī)理研究

    天然河道在不同類型水工程的作用下,被迫演變成長(zhǎng)河系、多阻斷格局,連續(xù)的洪水過(guò)程因多個(gè)串并聯(lián)相接的河庫(kù)系統(tǒng)單元呈現(xiàn)破碎化特征,而水庫(kù)蓄泄變化、蓄滯洪區(qū)分/滯洪過(guò)程的雙向互饋?zhàn)饔?,?dǎo)致洪水波動(dòng)特性改變,演進(jìn)規(guī)律呈現(xiàn)傳播時(shí)間縮短、坦化特性改變等新特征,迫切需要開展水工程阻斷斷面上、下游洪水波形態(tài)及關(guān)鍵區(qū)域分流變化特征和演變機(jī)理研究,以重構(gòu)河庫(kù)系統(tǒng)影響下洪水傳播特性變異的新認(rèn)識(shí),探明多阻斷條件下河道洪水傳播機(jī)理。

    3.3 聚焦強(qiáng)人類活動(dòng)影響特征,升級(jí)專業(yè)模型

    3.3.1 未控區(qū)間流域預(yù)報(bào)模型

    水工程建成以后,水工程之間、水工程與防洪保護(hù)點(diǎn)之間,形成了數(shù)目可觀的未控區(qū)間流域,上游區(qū)間的預(yù)報(bào)關(guān)乎水工程入流,進(jìn)而影響其調(diào)度運(yùn)用,而下游區(qū)間的預(yù)報(bào)則直接影響對(duì)下游防洪點(diǎn)的補(bǔ)償調(diào)度是否有效,因此,未控區(qū)間流域的預(yù)報(bào)一直是“河庫(kù)系統(tǒng)”洪水預(yù)報(bào)的難點(diǎn)痛點(diǎn)之一。一方面,未控區(qū)間的河流多屬中小河流,實(shí)時(shí)監(jiān)測(cè)站網(wǎng)和降雨預(yù)報(bào)分區(qū)的細(xì)化程度不足,區(qū)間來(lái)水的真值難以估計(jì),導(dǎo)致預(yù)報(bào)模型一般采用參數(shù)臨近移植法構(gòu)建,亟待深入開展率定檢驗(yàn)和參數(shù)修正研究;另一方面,在部分未控區(qū)間流域,尤其像長(zhǎng)江中下游干流區(qū)間覆蓋著數(shù)目眾多的排澇泵站、涵閘、洲灘民垸和中小水庫(kù),其在防洪緊張期的集中啟用,對(duì)預(yù)報(bào)的累積影響甚大,如2016年、2020年,長(zhǎng)江中游連續(xù)出現(xiàn)“中洪高水位”現(xiàn)象。加強(qiáng)實(shí)時(shí)信息獲取力度、動(dòng)態(tài)檢驗(yàn)參數(shù)、細(xì)化降雨預(yù)報(bào)分區(qū),是提高未控區(qū)間流域預(yù)報(bào)精度的主要途徑。

    3.3.2 多阻斷條件下的河道洪水演進(jìn)模型

    水工程的建成運(yùn)用改變了壩上、壩下洪水波形態(tài),原有河道匯流模型的適用性面臨挑戰(zhàn),需針對(duì)性開展庫(kù)區(qū)及壩下河道洪水演進(jìn)模型研究。對(duì)于河道型水庫(kù)庫(kù)區(qū)而言,入庫(kù)點(diǎn)因庫(kù)水位和入庫(kù)來(lái)水不同而不同,平水區(qū)、過(guò)渡區(qū)、河道區(qū)的范圍動(dòng)態(tài)變化,入庫(kù)流量的計(jì)算需分類適配,而由于動(dòng)庫(kù)容影響的存在,調(diào)洪計(jì)算宜采用水力學(xué)模型為主、靜庫(kù)容調(diào)洪模型為輔的多模型并行求解方法;對(duì)于壩下臨近河段,水庫(kù)出庫(kù)的激增或急減,使得河道洪水波呈現(xiàn)斷波特征,基于新樣本空間重構(gòu)動(dòng)態(tài)相關(guān)圖、分類馬斯京根等水文學(xué)模型,并探索基于黎曼間斷解的水力學(xué)求解方法,是提升下游河道洪水演進(jìn)預(yù)報(bào)精度的努力方向,對(duì)解決下游防洪保護(hù)點(diǎn)(如長(zhǎng)江中游沙市、城陵磯)的精準(zhǔn)補(bǔ)償調(diào)度具有重要價(jià)值;對(duì)于平原河網(wǎng)區(qū),研究構(gòu)建集水庫(kù)、蓄滯洪區(qū)、洲灘民垸、涵閘泵站、水文站、河湖于一體的一、二維水動(dòng)力學(xué)聯(lián)解模型,量化各類水工程對(duì)洪水的調(diào)蓄和滯排作用,可進(jìn)一步提升洪水精準(zhǔn)模擬預(yù)報(bào)能力。

    3.3.3 多元驅(qū)動(dòng)的中長(zhǎng)期徑流預(yù)測(cè)模型

    隨著流域內(nèi)水工程群的建成運(yùn)用,對(duì)中長(zhǎng)期預(yù)測(cè)的需求由定性轉(zhuǎn)向定量,而中長(zhǎng)期尺度的天氣系統(tǒng)變化機(jī)理尚不明晰,且徑流分布受水工程累積的影響越來(lái)越顯著,因此中長(zhǎng)期來(lái)水量預(yù)測(cè)成為當(dāng)前水文預(yù)報(bào)領(lǐng)域的又一重點(diǎn)、難點(diǎn)。建立穩(wěn)定可靠的長(zhǎng)系列資料庫(kù),開展中長(zhǎng)期徑流預(yù)測(cè)關(guān)鍵影響因子和水庫(kù)蓄泄因子分析,利用人工智能、機(jī)器學(xué)習(xí)、集合預(yù)報(bào)等新方法,耦合機(jī)理模型、數(shù)據(jù)挖掘、知識(shí)規(guī)則等多種模型實(shí)現(xiàn)中長(zhǎng)期趨勢(shì)分析和徑流預(yù)測(cè),進(jìn)而形成一套多元驅(qū)動(dòng)的流域中長(zhǎng)期徑流預(yù)測(cè)技術(shù),可望進(jìn)一步提高中長(zhǎng)期徑流預(yù)報(bào)精度。

    3.4 深化預(yù)報(bào)新技術(shù)研究,探索基于智能預(yù)報(bào)調(diào)度的一體化業(yè)務(wù)體系

    3.4.1 基于相似雨洪和參數(shù)適配的智能預(yù)報(bào)體系

    在實(shí)時(shí)作業(yè)預(yù)報(bào)中,雨洪規(guī)律和參數(shù)適應(yīng)性的認(rèn)知對(duì)預(yù)報(bào)結(jié)果產(chǎn)生重要影響。為了提高預(yù)報(bào)的客觀化水平,實(shí)現(xiàn)智能化模擬預(yù)報(bào),開展基于相似雨洪和參數(shù)智適配預(yù)報(bào)是兩條有效的解決途徑。一方面,提出暴雨洪水特征表征的指標(biāo)體系和規(guī)范化計(jì)量方法,準(zhǔn)確刻畫暴雨洪水特征指標(biāo),引入機(jī)器學(xué)習(xí)和圖像識(shí)別等算法,實(shí)現(xiàn)暴雨洪水關(guān)系挖掘綜合分析,即可辨識(shí)相似雨洪,由預(yù)見期降雨智能預(yù)報(bào)出相似洪水過(guò)程;另一方面,考慮不同特點(diǎn)的水文過(guò)程應(yīng)有不同適應(yīng)參數(shù),構(gòu)建動(dòng)態(tài)成長(zhǎng)的分類模型參數(shù)庫(kù),利用人工智能新技術(shù)創(chuàng)新流域水文模型參數(shù)的智適配方法,根據(jù)降雨預(yù)報(bào)智能匹配水文模型參數(shù),計(jì)算出精準(zhǔn)的洪水過(guò)程。

    3.4.2 基于調(diào)度規(guī)則智能模擬的精準(zhǔn)化決策體系

    水工程調(diào)度當(dāng)前已成為制約預(yù)報(bào)精度和預(yù)見期的主要影響因素之一,亟須研究調(diào)度規(guī)則的智能化模擬技術(shù),實(shí)現(xiàn)水工程調(diào)度的客觀化、智能化,以有效延長(zhǎng)預(yù)見期?;趯?shí)時(shí)調(diào)度規(guī)則庫(kù)、歷史洪水調(diào)度案例庫(kù),采用知識(shí)圖譜技術(shù),將預(yù)報(bào)調(diào)度相關(guān)的原理規(guī)律、規(guī)則、經(jīng)驗(yàn)、策略、方法等信息結(jié)構(gòu)化,并結(jié)合神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)、深度學(xué)習(xí)等智能技術(shù),構(gòu)建形成可累積式發(fā)展的實(shí)時(shí)洪水調(diào)度知識(shí)庫(kù)。根據(jù)水文氣象預(yù)報(bào)成果和流域防洪形勢(shì)分析,基于知識(shí)圖譜推理方法,智能推薦調(diào)度方案集,并進(jìn)行專家交互,形式可支撐調(diào)度的精準(zhǔn)化決策。

    3.4.3 基于概率的實(shí)時(shí)預(yù)報(bào)調(diào)度風(fēng)險(xiǎn)決策體系

    初始狀態(tài)誤差、輸入誤差、模型結(jié)構(gòu)和參數(shù)誤差使得洪水預(yù)報(bào)難免存在不確定性,確定性預(yù)報(bào)能提供的風(fēng)險(xiǎn)信息有限,為客觀描述預(yù)報(bào)的不確定性,從而為優(yōu)化決策、精細(xì)調(diào)度提供更好支撐,需從以下幾個(gè)方面開展預(yù)報(bào)不確定性量化及風(fēng)險(xiǎn)決策的關(guān)鍵技術(shù)攻關(guān):從總誤差途徑和集合途徑,研究可靠的概率預(yù)報(bào)模型;探索耦合概率預(yù)報(bào)的調(diào)度風(fēng)險(xiǎn)量化方法及風(fēng)險(xiǎn)決策技術(shù);完成概率預(yù)報(bào)和風(fēng)險(xiǎn)決策的業(yè)務(wù)化設(shè)計(jì),構(gòu)建基于風(fēng)險(xiǎn)的實(shí)時(shí)預(yù)報(bào)調(diào)度決策業(yè)務(wù)體系。

    3.5 打造通用化業(yè)務(wù)化平臺(tái),支撐多要素全周期“四預(yù)”功能

    為克服業(yè)界水文模型結(jié)構(gòu)不統(tǒng)一,建模流程不規(guī)范、客觀化程度不高等問題,亟待制定標(biāo)準(zhǔn)化的模型封裝、集成技術(shù),形成模型服務(wù)開發(fā)標(biāo)準(zhǔn)體系,構(gòu)建通用、共享的洪水預(yù)報(bào)模型庫(kù);融合流域自動(dòng)圈化、資料校核處理、傳播時(shí)間分析、參數(shù)智能率定、方案評(píng)價(jià)、預(yù)報(bào)體系敏捷搭建等功能,研發(fā)通用化、向?qū)降乃念A(yù)報(bào)建模平臺(tái),為“四預(yù)”提供良好的基礎(chǔ)支撐。

    順應(yīng)數(shù)字孿生流域需求,基于短中長(zhǎng)期相結(jié)合、水文氣象耦合、預(yù)報(bào)調(diào)度互饋的技術(shù)思路,構(gòu)建數(shù)值模式、水文模型、水力學(xué)模型、數(shù)理統(tǒng)計(jì)模型、智能模型、概率預(yù)報(bào)模型等多元模型驅(qū)動(dòng)和“分布式模擬-水工程調(diào)度-河道演進(jìn)”多層嵌套的流域水文模擬預(yù)報(bào)引擎,可實(shí)現(xiàn)大江大河洪水預(yù)報(bào)、山洪和中小河流預(yù)報(bào)、旱情預(yù)測(cè)、水資源預(yù)測(cè)、突發(fā)水事件水文應(yīng)急預(yù)報(bào)等多要素預(yù)報(bào)預(yù)警。研發(fā)集全周期預(yù)測(cè)預(yù)報(bào)、多要素監(jiān)視預(yù)警、水工程綜合調(diào)度預(yù)演、全鏈條預(yù)案生成與管理為一體的“四預(yù)”數(shù)字孿生平臺(tái)。智慧流域水文預(yù)報(bào)關(guān)鍵技術(shù)見圖1。

    4 結(jié) 語(yǔ)

    中國(guó)已進(jìn)入智慧水利、數(shù)字孿生流域建設(shè)的水利高質(zhì)量發(fā)展新階段,流域面臨人類活動(dòng)和氣候變化雙重影響的變化環(huán)境。本文通過(guò)綜述國(guó)內(nèi)外洪水預(yù)報(bào)研究和應(yīng)用的進(jìn)展,分析了目前流域水文預(yù)報(bào)業(yè)務(wù)發(fā)展的需求和面臨的挑戰(zhàn),并對(duì)智慧流域水文預(yù)報(bào)技術(shù)研究進(jìn)行了探討和展望,得到主要認(rèn)識(shí)和建議如下:

    (1) 制約流域產(chǎn)匯流及洪水預(yù)報(bào)模型研究和預(yù)報(bào)精度的關(guān)鍵是現(xiàn)有可用的資料十分有限,雨量和流量觀測(cè)站點(diǎn)稀少,僅依靠流域出口控制斷面流量系列來(lái)率定和檢驗(yàn)?zāi)P?,預(yù)報(bào)精度有待提高。特別是中西部欠發(fā)達(dá)地區(qū)和未控區(qū)間流域的產(chǎn)匯流計(jì)算,由于站網(wǎng)密度不足或無(wú)資料,直接導(dǎo)致水文模擬預(yù)報(bào)精度不高。因此,智慧流域建設(shè)首先應(yīng)加強(qiáng)空天地一體化監(jiān)測(cè)工作,收集整編流域氣象、水文、地理地貌、土地利用、植被覆蓋、水利交通等基礎(chǔ)數(shù)據(jù),并建立標(biāo)準(zhǔn)統(tǒng)一的國(guó)家基礎(chǔ)數(shù)據(jù)庫(kù)。

    (2) 新時(shí)期新形勢(shì)對(duì)流域水文預(yù)報(bào)的發(fā)展提出了新的需求:① 提高預(yù)報(bào)精度延長(zhǎng)預(yù)見期;② 客觀量化預(yù)報(bào)不確定性,更好支撐風(fēng)險(xiǎn)調(diào)度決策;③ 提高智慧化水平,解放生產(chǎn)力。面對(duì)新需求,流域洪水預(yù)報(bào)在算據(jù)獲取、預(yù)報(bào)模型適應(yīng)性、新技術(shù)落地推廣和業(yè)務(wù)系統(tǒng)通用化等方面仍面臨較大挑戰(zhàn)。

    (3) 開展智慧流域產(chǎn)匯流及洪水預(yù)報(bào)模型研發(fā),首先要聚焦數(shù)字孿生流域建設(shè)的“算據(jù)”支撐需求,實(shí)現(xiàn)多源異構(gòu)數(shù)據(jù)匯集,夯實(shí)流域數(shù)據(jù)底板;針對(duì)氣候變化和人類活動(dòng)深刻影響流域的背景,深入研究流域水汽循環(huán)和產(chǎn)匯流基礎(chǔ)規(guī)律,提升對(duì)流域暴雨洪水特性的認(rèn)知水平。

    (4) 抓住水工程群建成后的“大流域、長(zhǎng)河系、多阻斷”河庫(kù)系統(tǒng)新格局,開展水工程互饋影響下的水文預(yù)報(bào)模型及方法研究,升級(jí)專業(yè)模型平臺(tái),為數(shù)字孿生流域提供“算法”支撐;圍繞數(shù)字孿生流域建設(shè)總體目標(biāo)的智慧化模擬、精準(zhǔn)化決策路線,深化預(yù)報(bào)新技術(shù)研究,探索基于智能、風(fēng)險(xiǎn)的預(yù)報(bào)調(diào)度業(yè)務(wù)體系;以標(biāo)準(zhǔn)化、國(guó)產(chǎn)化、客觀化、智能化為目標(biāo),構(gòu)建通用型業(yè)務(wù)平臺(tái),提供多要素全周期水文預(yù)測(cè)預(yù)報(bào)預(yù)警服務(wù),有力支撐“四預(yù)”孿生功能。

    參考文獻(xiàn):

    [1] 李國(guó)英.建設(shè)數(shù)字孿生流域 推動(dòng)新階段水利高質(zhì)量發(fā)展[N].學(xué)習(xí)時(shí)報(bào),2022-06-29(1).

    [2] 水利部.數(shù)字孿生流域建設(shè)技術(shù)大綱(試行)[Z].2022.

    [3] 夏軍,石衛(wèi).變化環(huán)境下中國(guó)水安全問題研究與展望[J].水利學(xué)報(bào),2016,47(3):292-301.

    [4] 王俊.面向水庫(kù)群調(diào)度的水文數(shù)值模擬與預(yù)測(cè)技術(shù)[M].北京:中國(guó)水利水電出版社,2021.

    [5] 吳雷,許有朋,王躍峰,等.水文實(shí)驗(yàn)研究進(jìn)展[J].水科學(xué)進(jìn)展,2017,28(4):632-640.

    [6] 趙人俊.流域水文模擬:新安江模型與陜北模型[M].北京:水利電力出版社,1984.

    [7] FREEZE R A,HARLAN R.Blueprint for a physically based,digitally simulated hydrologic response model [J].Journal of Hydrology,1969,9(3):237-258.

    [8] 任立良,劉新仁.數(shù)字高程模型信息提取與數(shù)字水文模型研究進(jìn)展[J].水科學(xué)進(jìn)展,2000,11(4):463-469.

    [9] 賈仰文,王浩,倪廣恒,等.分布式流域水文模型原理與實(shí)踐[M].北京:中國(guó)水利水電出版社,2005.

    [10] 熊立華,郭生練.分布式流域水文模型[M].北京:中國(guó)水利水電出版社,2004.

    [11] 徐宗學(xué),程磊.分布式水文模型研究與應(yīng)用進(jìn)展[J].水利學(xué)報(bào),2010,41(9):1009-1017.

    [12] 周研來(lái),郭生練,張斐章,等.人工智能在水文預(yù)報(bào)中的應(yīng)用研究[J].水資源研究,2019,8(1):1-12.

    [13] 張珂,牛杰帆,李曦,等.洪水預(yù)報(bào)智能模型在中國(guó)半干旱半濕潤(rùn)區(qū)的應(yīng)用對(duì)比[J].水資源保護(hù),2021,37(1):28-35,60.

    [14] 胡義明,陳騰,羅序義,等.基于機(jī)器學(xué)習(xí)模型的淮河流域中長(zhǎng)期徑流預(yù)報(bào)研究[J].地學(xué)前緣,2022,29(3):284-291.

    [15] RUMELHART D E,HINTON G E,WILLIAMS R J.Learning representations by back-propagating errors[J].Nature,1986,323(6088):533-536.

    [16] HOCHREITER S,SCHMIDHUBER J.Long short-term memory[J].Neural Computation,1997,9(8):1735-1780.

    [17] KRATZERT F,KLOTZ D,BRENNER C,et al.Rainfall-runoff modelling using long short-term memory(LSTM) networks[J].Hydrology and Earth System Sciences,2018,22(11):6005-6022.

    [18] NEARING G S,KRATZERT F,SAMPSON A K,et al.What role does hydrological science play in the age of machine learning [J].Water Resources Research,2021,57(3):e2020WR028091.

    [19] HU C,WU Q,LI H,et al.Deep learning with a long short-term memory networks approach for rainfall-runoff simulation[J].Water,2018,10(11):1543.

    [20] 徐源浩,鄔強(qiáng),李常青,等.基于長(zhǎng)短時(shí)記憶(LSTM)神經(jīng)網(wǎng)絡(luò)的黃河中游洪水過(guò)程模擬及預(yù)報(bào)[J].北京師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2020,56(3):387-393.

    [21] XU Y,HU C,WU Q,et al.Research on particle swarm optimization in LSTM neural networks for rainfall-runoff simulation[J].Journal of Hydrology,2022,608:127553.

    [22] 李步,田富強(qiáng),李鈺坤,等.融合氣象要素時(shí)空特征的深度學(xué)習(xí)水文模型[J].水科學(xué)進(jìn)展,2022,33(6):904-913.

    [23] 張海榮,湯正陽(yáng),曹輝.基于偏互信息篩選的中長(zhǎng)期徑流預(yù)報(bào)方案比選[J].人民長(zhǎng)江,2022,53(10):71-75.

    [24] KAO I,ZHOU Y,CHANG L,et al.Exploring along short-term memory based encoder-decoder framework for multi-step-ahead flood forecasting[J].Journal of Hydrology,2020,583.

    [25] HAN H,CHOI C,JUNG J,et al.Deeplearning with long short term memory based sequence-to-sequence model for rainfall-runoff simulation[J].Water,2021,13(4).

    [26] CUI Z,ZHOU Y L,GUO S L,et al.Effective improvement of multi-step-ahead flood forecasting accuracy through encoder-decoder with an exogenous input structure[J].Journal of Hydrology,2022,609(4),127764.

    [27] 林康聆,陳華,陳清勇,等.耦合Encoder-Decoder的LSTM徑流預(yù)報(bào)模型研究[J].武漢大學(xué)學(xué)報(bào)(工學(xué)版),2022,55(8):755-761.

    [28] GIRIHAGAMA L,NAVEED KHALIQ M,et al.Streamflow modelling and forecasting for Canadian watersheds using LSTM networks with attention mechanism[J].Neural Computing and Applications,2022,34(22):19995-20015.

    [29] 崔震,郭生練,汪蕓,等.三峽入庫(kù)洪水概率預(yù)報(bào)的深度學(xué)習(xí)模型[J].水力發(fā)電學(xué)報(bào),2023,42(4):1-10.

    [30] 陸桂華,吳娟,吳志勇.水文集合預(yù)報(bào)試驗(yàn)及其研究進(jìn)展[J].水科學(xué)進(jìn)展,2012,23(5):728-734.

    [31] 張俊,楊文發(fā),牛文靜,等.長(zhǎng)江流域洪水概率預(yù)報(bào)業(yè)務(wù)化方法探討[J].人民長(zhǎng)江,2021,52(9):11-15.

    [32] 巴歡歡,郭生練,鐘逸軒,等.考慮降水預(yù)報(bào)的三峽入庫(kù)洪水集合概率預(yù)報(bào)方法比較[J].水科學(xué)進(jìn)展,2019,30(2):187-197.

    [33] RAFTERY A E,GNEITING T,BALABDAOUI F,et al.Using Bayesian model averaging to calibrate forecast ensembles[J].Monthly Weather Review,2005,133(5):1155-1174.

    [34] GNEITING T,RAFTERY A E,WESTVELD A H,et al.Calibrated probabilistic forecasting using ensemble model output statistics and minimum CRPS estimation [J].Monthly Weather Review,2005,133(5):1098-1119.

    [35] KRZYSZTOFOWICZ R.Bayesian theory of probabilistic forecasting via deterministic hydrologic model [J].Water Resources Research,1999,35(9):2739-2750.

    [36] 劉章君,郭生練,李天元,等.貝葉斯概率洪水預(yù)報(bào)模型及其比較應(yīng)用研究[J].水利學(xué)報(bào),2014,45(9):1019-1028.

    [37] 梁忠民,蔣曉蕾,錢名開,等.考慮誤差異分布的洪水概率預(yù)報(bào)方法研究[J].水力發(fā)電學(xué)報(bào),2017,36(4):18-25.

    [38] 張建云.中國(guó)水文預(yù)報(bào)技術(shù)發(fā)展的回顧與思考[J].水科學(xué)進(jìn)展,2010,21(4):435-443.

    [39] 葛守西.現(xiàn)代洪水預(yù)報(bào)技術(shù)[M].北京:中國(guó)水利水電出版社,1999.

    [40] 王威,施勇,陳敏.淺談長(zhǎng)江中下游實(shí)時(shí)洪水調(diào)度系統(tǒng)開發(fā)[J].人民長(zhǎng)江,2006,37(9):58-62.

    [41] 章四龍.中國(guó)洪水預(yù)報(bào)系統(tǒng)設(shè)計(jì)建設(shè)研究[J].水文,2002(1):32-34,16.

    [42] 鄒冰玉,高珺,李玉榮.通用型水文預(yù)報(bào)平臺(tái)開發(fā)與應(yīng)用[J].人民長(zhǎng)江,2011,42(6):101-105.

    [43] 秦昊,陳瑜彬.長(zhǎng)江洪水預(yù)報(bào)調(diào)度系統(tǒng)建設(shè)及應(yīng)用[J].人民長(zhǎng)江,2017,48(4):16-21.

    [44] 程海云.洪水預(yù)報(bào)系統(tǒng)及其在工程調(diào)度中的作用[J].中國(guó)水利,2020(17):11-13.

    [45] 芮孝芳.對(duì)流域水文模型的再認(rèn)識(shí)[J].水利水電科技進(jìn)展,2018,38(2):1-7.

    [46] 李致家.現(xiàn)代水文模擬與預(yù)報(bào)技術(shù)[M].南京:河海大學(xué)出版社,2010.

    [47] 黃艷,喻杉,羅斌,等.面向流域水工程防災(zāi)聯(lián)合智能調(diào)度的數(shù)字孿生長(zhǎng)江探索[J].水利學(xué)報(bào),2022,53(3):253-269.

    (編輯:江 文)

    Abstract:

    At present,traditional water conservancy is rapidly entering a new stage of high-quality development for intelligent water conservancy.It is an important component of building a digital twin watershed to develop an intelligent watershed hydrological simulation and flood forecasting system.It is also an inevitable requirement for strengthening the "four precautions"(forecasting,alarming,preview and plan) and promoting the construction of intelligent water conservancy.The research progress of hydrological simulation forecasting is reviewed,and the current situation and shortcomings of the application of flood forecasting systems in China are analyzed.Meanwhile,we discuss the key technologies and development prospect of intelligent basin hydrological simulation and flood forecasting system from the aspects of database,regular mechanism,forecast model,methodology and technology,and application platform,including: pooling multi-source heterogeneous data and consolidating basin database;strengthening the water vapor cycle study of river reservoir system and revealing the cause mechanism;focusing on characteristics of strong human activities and upgrading professional models;deepening the research of new forecasting technology and exploring the integrated business system based on intelligent forecasting and dispatching;building a general business platform to support the multi-factor full cycle "four precautions" function.The research results can provide a theoretical basis and technical support for improving the flood forecasting level in river basins,strengthening the "four precautions",and constructing digital dual river basin.

    Key words:

    intelligent watershed;digital twin;hydrological forecasting;intelligent simulation;probabilistic forecasting

    猜你喜歡
    水文洪水流域
    2022年《中國(guó)水文年報(bào)》發(fā)布
    壓油溝小流域
    堡子溝流域綜合治理
    羅堰小流域
    水文
    水文水資源管理
    洪水時(shí)遇到電線低垂或折斷該怎么辦
    水利規(guī)劃與設(shè)計(jì)(2018年1期)2018-01-31 01:53:37
    又見洪水(外二首)
    水文
    亚洲精品影视一区二区三区av| 国内久久婷婷六月综合欲色啪| 一卡2卡三卡四卡精品乱码亚洲| 日本在线视频免费播放| 制服人妻中文乱码| 国产在线精品亚洲第一网站| 国产亚洲精品综合一区在线观看| 国产精品99久久久久久久久| 亚洲国产精品合色在线| 色哟哟哟哟哟哟| 在线观看av片永久免费下载| 人人妻人人澡欧美一区二区| 18禁黄网站禁片免费观看直播| 日韩 欧美 亚洲 中文字幕| 3wmmmm亚洲av在线观看| 久久99热这里只有精品18| www.色视频.com| 十八禁网站免费在线| 成年女人永久免费观看视频| 我要搜黄色片| 午夜精品在线福利| 国产69精品久久久久777片| 少妇的丰满在线观看| 成人18禁在线播放| 国产高清三级在线| 国产伦一二天堂av在线观看| 少妇高潮的动态图| av在线蜜桃| 人人妻人人澡欧美一区二区| 欧美在线黄色| 18禁美女被吸乳视频| 中文字幕精品亚洲无线码一区| 国产激情欧美一区二区| 国内少妇人妻偷人精品xxx网站| 国产成年人精品一区二区| 成人无遮挡网站| 蜜桃久久精品国产亚洲av| 757午夜福利合集在线观看| 国产午夜精品久久久久久一区二区三区 | 香蕉av资源在线| 男女视频在线观看网站免费| 制服人妻中文乱码| 亚洲成人久久性| 亚洲精品色激情综合| 在线观看午夜福利视频| 人妻丰满熟妇av一区二区三区| 成人特级黄色片久久久久久久| 欧美午夜高清在线| 精品一区二区三区视频在线 | or卡值多少钱| 亚洲欧美激情综合另类| 中文资源天堂在线| 中文字幕人成人乱码亚洲影| 老鸭窝网址在线观看| 熟女少妇亚洲综合色aaa.| 国产毛片a区久久久久| 亚洲成人中文字幕在线播放| 亚洲精品久久国产高清桃花| 99热精品在线国产| 人人妻,人人澡人人爽秒播| 亚洲无线在线观看| 久久草成人影院| 又黄又粗又硬又大视频| 国产高清激情床上av| 欧美在线一区亚洲| 久久伊人香网站| 婷婷精品国产亚洲av| 久久精品夜夜夜夜夜久久蜜豆| 欧美在线一区亚洲| 夜夜躁狠狠躁天天躁| 国内精品美女久久久久久| 日本精品一区二区三区蜜桃| 在线观看免费午夜福利视频| 中文字幕av在线有码专区| 男人舔奶头视频| 九九热线精品视视频播放| 亚洲人与动物交配视频| 一个人看视频在线观看www免费 | 可以在线观看毛片的网站| 日韩人妻高清精品专区| 深爱激情五月婷婷| 色综合欧美亚洲国产小说| 日韩精品中文字幕看吧| 色综合亚洲欧美另类图片| 高清毛片免费观看视频网站| 国产欧美日韩精品一区二区| 免费看a级黄色片| 欧美一区二区亚洲| 人妻久久中文字幕网| 久久精品国产亚洲av涩爱 | 免费搜索国产男女视频| 日韩欧美在线二视频| 三级男女做爰猛烈吃奶摸视频| 国产欧美日韩一区二区三| 久久久久九九精品影院| 啦啦啦观看免费观看视频高清| 人妻夜夜爽99麻豆av| 久久久久久国产a免费观看| 热99re8久久精品国产| 亚洲午夜理论影院| 国产主播在线观看一区二区| www.www免费av| 两个人看的免费小视频| 国产三级中文精品| 女人十人毛片免费观看3o分钟| 久久伊人香网站| 亚洲成人免费电影在线观看| 综合色av麻豆| 午夜福利成人在线免费观看| 女同久久另类99精品国产91| 19禁男女啪啪无遮挡网站| 中文字幕精品亚洲无线码一区| 色吧在线观看| 69人妻影院| 人妻久久中文字幕网| 三级男女做爰猛烈吃奶摸视频| 在线免费观看的www视频| 嫩草影院入口| 最新美女视频免费是黄的| 精品国产超薄肉色丝袜足j| 久久香蕉国产精品| 亚洲欧美精品综合久久99| h日本视频在线播放| 免费电影在线观看免费观看| 欧美日韩福利视频一区二区| 亚洲最大成人中文| 欧美黑人欧美精品刺激| av片东京热男人的天堂| 在线a可以看的网站| 岛国视频午夜一区免费看| 操出白浆在线播放| 少妇高潮的动态图| 老汉色av国产亚洲站长工具| 性欧美人与动物交配| 国内精品一区二区在线观看| 国产欧美日韩一区二区三| 真人一进一出gif抽搐免费| 中文字幕高清在线视频| 草草在线视频免费看| 午夜福利视频1000在线观看| av专区在线播放| 人人妻人人澡欧美一区二区| 十八禁人妻一区二区| 午夜精品一区二区三区免费看| 国产爱豆传媒在线观看| 叶爱在线成人免费视频播放| 色尼玛亚洲综合影院| 在线天堂最新版资源| 久久香蕉国产精品| 精品午夜福利视频在线观看一区| 久久香蕉精品热| 床上黄色一级片| 欧美中文日本在线观看视频| 精品国产超薄肉色丝袜足j| 亚洲av日韩精品久久久久久密| 久久精品国产99精品国产亚洲性色| 亚洲专区国产一区二区| 久久精品综合一区二区三区| 欧美极品一区二区三区四区| 国产高潮美女av| 欧美日韩乱码在线| 人妻丰满熟妇av一区二区三区| 女人高潮潮喷娇喘18禁视频| 国产成人影院久久av| 欧美最黄视频在线播放免费| eeuss影院久久| 嫁个100分男人电影在线观看| 丁香六月欧美| 国产免费av片在线观看野外av| 精品一区二区三区视频在线 | 午夜老司机福利剧场| av专区在线播放| 91av网一区二区| 俺也久久电影网| 免费一级毛片在线播放高清视频| www.色视频.com| bbb黄色大片| 亚洲最大成人中文| 欧美一区二区国产精品久久精品| 精品久久久久久成人av| 欧美又色又爽又黄视频| 午夜福利视频1000在线观看| 欧美日韩瑟瑟在线播放| 哪里可以看免费的av片| svipshipincom国产片| 动漫黄色视频在线观看| 欧美成人性av电影在线观看| 国内精品久久久久精免费| 无遮挡黄片免费观看| 18美女黄网站色大片免费观看| 免费人成在线观看视频色| 国产午夜精品久久久久久一区二区三区 | 99热这里只有精品一区| 不卡一级毛片| 免费人成视频x8x8入口观看| 一级毛片高清免费大全| 精品国产美女av久久久久小说| 国产精品美女特级片免费视频播放器| 黄色日韩在线| av在线天堂中文字幕| 丰满人妻熟妇乱又伦精品不卡| 我的老师免费观看完整版| 在线观看美女被高潮喷水网站 | 在线国产一区二区在线| 美女高潮喷水抽搐中文字幕| 国产单亲对白刺激| 国产爱豆传媒在线观看| 91久久精品国产一区二区成人 | 国产欧美日韩一区二区三| 午夜激情福利司机影院| 亚洲电影在线观看av| 叶爱在线成人免费视频播放| 成人一区二区视频在线观看| 黄色女人牲交| 成人鲁丝片一二三区免费| 欧美三级亚洲精品| 日日干狠狠操夜夜爽| 欧美最黄视频在线播放免费| 国产精品一及| 夜夜夜夜夜久久久久| 女同久久另类99精品国产91| 听说在线观看完整版免费高清| 国产高清videossex| 日韩欧美国产一区二区入口| 日韩高清综合在线| 国产亚洲精品综合一区在线观看| 长腿黑丝高跟| 国产久久久一区二区三区| 国产乱人视频| 日韩亚洲欧美综合| 日本成人三级电影网站| 国产高清有码在线观看视频| 免费一级毛片在线播放高清视频| 久久99热这里只有精品18| 免费在线观看亚洲国产| 成人av一区二区三区在线看| 午夜老司机福利剧场| 成人永久免费在线观看视频| 久久久国产精品麻豆| svipshipincom国产片| 精品人妻1区二区| 男女视频在线观看网站免费| 亚洲国产欧美网| 国产私拍福利视频在线观看| 色av中文字幕| 亚洲欧美日韩无卡精品| x7x7x7水蜜桃| 日韩免费av在线播放| 国产精品久久视频播放| 色老头精品视频在线观看| 久久精品综合一区二区三区| 国产aⅴ精品一区二区三区波| 午夜福利18| 午夜精品在线福利| av福利片在线观看| 日日干狠狠操夜夜爽| 一级毛片高清免费大全| 变态另类丝袜制服| 久久久精品大字幕| 黄片小视频在线播放| 国产精品久久久久久久久免 | 精华霜和精华液先用哪个| 男插女下体视频免费在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 一个人免费在线观看的高清视频| 看免费av毛片| 国产伦精品一区二区三区视频9 | 人人妻,人人澡人人爽秒播| 中亚洲国语对白在线视频| 国产亚洲精品一区二区www| 99精品欧美一区二区三区四区| 婷婷丁香在线五月| 久久国产乱子伦精品免费另类| 一级作爱视频免费观看| 身体一侧抽搐| 国产伦人伦偷精品视频| 丰满人妻熟妇乱又伦精品不卡| 变态另类成人亚洲欧美熟女| 黄片小视频在线播放| 国产日本99.免费观看| www日本在线高清视频| 欧美午夜高清在线| 欧美乱妇无乱码| 日韩欧美在线乱码| 国产亚洲欧美在线一区二区| 午夜精品一区二区三区免费看| 日本 av在线| 性欧美人与动物交配| 中文字幕av成人在线电影| 俺也久久电影网| 国产精品国产高清国产av| 精品久久久久久久毛片微露脸| 亚洲在线自拍视频| 国产伦精品一区二区三区视频9 | 久久精品人妻少妇| 日本一二三区视频观看| 亚洲成人免费电影在线观看| 99riav亚洲国产免费| 手机成人av网站| 日韩欧美免费精品| 国产成+人综合+亚洲专区| 有码 亚洲区| 黑人欧美特级aaaaaa片| 久久中文看片网| 人妻久久中文字幕网| 亚洲美女黄片视频| 亚洲一区高清亚洲精品| 午夜免费成人在线视频| 女警被强在线播放| 国产高清视频在线播放一区| 午夜精品久久久久久毛片777| 日韩中文字幕欧美一区二区| 99riav亚洲国产免费| 成年版毛片免费区| 一个人观看的视频www高清免费观看| 搡老熟女国产l中国老女人| 琪琪午夜伦伦电影理论片6080| 亚洲在线观看片| 国产私拍福利视频在线观看| 最近最新中文字幕大全电影3| 久久亚洲精品不卡| 内地一区二区视频在线| 丁香欧美五月| 成年女人毛片免费观看观看9| 国产男靠女视频免费网站| 成人av在线播放网站| 日本黄色视频三级网站网址| 久久性视频一级片| 亚洲成人免费电影在线观看| 免费观看的影片在线观看| h日本视频在线播放| 91在线精品国自产拍蜜月 | 一区二区三区免费毛片| 在线观看免费视频日本深夜| www.www免费av| 波野结衣二区三区在线 | 好看av亚洲va欧美ⅴa在| 久久久久久大精品| 中文字幕高清在线视频| 最近视频中文字幕2019在线8| 国产精品野战在线观看| 九九热线精品视视频播放| 丁香六月欧美| 淫秽高清视频在线观看| 可以在线观看毛片的网站| www.999成人在线观看| 久久国产精品人妻蜜桃| 亚洲国产欧洲综合997久久,| 国产探花在线观看一区二区| 国产一区二区激情短视频| 一级毛片女人18水好多| 麻豆国产97在线/欧美| 国产色婷婷99| 国产黄a三级三级三级人| 免费观看精品视频网站| 少妇的逼好多水| 高潮久久久久久久久久久不卡| 88av欧美| 日韩大尺度精品在线看网址| 99精品欧美一区二区三区四区| 真实男女啪啪啪动态图| 99精品欧美一区二区三区四区| 男女那种视频在线观看| 国产国拍精品亚洲av在线观看 | 他把我摸到了高潮在线观看| 欧美大码av| 丰满的人妻完整版| 国产精品一区二区免费欧美| 哪里可以看免费的av片| 小蜜桃在线观看免费完整版高清| 亚洲性夜色夜夜综合| 在线播放无遮挡| 听说在线观看完整版免费高清| 久久精品国产综合久久久| or卡值多少钱| 欧美绝顶高潮抽搐喷水| 久久久久久久精品吃奶| 一级黄色大片毛片| 青草久久国产| 99久久久亚洲精品蜜臀av| tocl精华| 男人舔奶头视频| 很黄的视频免费| 人人妻,人人澡人人爽秒播| xxxwww97欧美| 97人妻精品一区二区三区麻豆| 国产精品久久久久久久久免 | 国产高清视频在线观看网站| 脱女人内裤的视频| 久久精品综合一区二区三区| 亚洲一区二区三区不卡视频| 久久6这里有精品| 白带黄色成豆腐渣| 嫁个100分男人电影在线观看| 国产99白浆流出| 最近在线观看免费完整版| 欧美一区二区精品小视频在线| 国产精品美女特级片免费视频播放器| 人人妻人人澡欧美一区二区| 亚洲av不卡在线观看| 色噜噜av男人的天堂激情| 亚洲精品乱码久久久v下载方式 | 搡老熟女国产l中国老女人| 超碰av人人做人人爽久久 | or卡值多少钱| 国产精品乱码一区二三区的特点| 中亚洲国语对白在线视频| 亚洲人成网站在线播放欧美日韩| 久久久国产成人免费| 大型黄色视频在线免费观看| 日本黄色片子视频| 精品国产三级普通话版| 两个人看的免费小视频| 床上黄色一级片| 色综合婷婷激情| 香蕉丝袜av| 色老头精品视频在线观看| 一区二区三区免费毛片| av在线蜜桃| 午夜激情福利司机影院| av国产免费在线观看| 亚洲精品乱码久久久v下载方式 | 亚洲成av人片在线播放无| 亚洲av一区综合| 欧美不卡视频在线免费观看| 欧美成人一区二区免费高清观看| 999久久久精品免费观看国产| 一级作爱视频免费观看| 国产野战对白在线观看| 欧洲精品卡2卡3卡4卡5卡区| 免费电影在线观看免费观看| 免费一级毛片在线播放高清视频| 少妇裸体淫交视频免费看高清| 99久久九九国产精品国产免费| 日韩高清综合在线| 精品人妻偷拍中文字幕| 少妇的丰满在线观看| 日本免费一区二区三区高清不卡| 可以在线观看的亚洲视频| 免费看a级黄色片| 人妻久久中文字幕网| 国产精品久久久久久精品电影| 黑人欧美特级aaaaaa片| 成年女人永久免费观看视频| 久久精品国产亚洲av涩爱 | 一边摸一边抽搐一进一小说| 12—13女人毛片做爰片一| 三级毛片av免费| 亚洲自拍偷在线| 亚洲av电影不卡..在线观看| 亚洲av美国av| 亚洲午夜理论影院| 欧美成人性av电影在线观看| 日韩中文字幕欧美一区二区| 在线免费观看不下载黄p国产 | 国产精品影院久久| 国产成人欧美在线观看| 久久婷婷人人爽人人干人人爱| 久久人妻av系列| 国产精品日韩av在线免费观看| 亚洲一区二区三区不卡视频| 中亚洲国语对白在线视频| 亚洲一区二区三区色噜噜| 国产中年淑女户外野战色| av国产免费在线观看| 中文字幕人妻熟人妻熟丝袜美 | 亚洲最大成人中文| 日韩高清综合在线| 精品人妻偷拍中文字幕| 少妇的丰满在线观看| 精品国产美女av久久久久小说| 99热精品在线国产| 免费看美女性在线毛片视频| 一级黄色大片毛片| 欧美日韩福利视频一区二区| 18禁黄网站禁片免费观看直播| 天堂影院成人在线观看| 成年女人永久免费观看视频| 亚洲国产精品合色在线| 午夜福利视频1000在线观看| 午夜免费成人在线视频| 亚洲国产中文字幕在线视频| 一二三四社区在线视频社区8| 亚洲午夜理论影院| 男女视频在线观看网站免费| 91久久精品电影网| 国产精品亚洲av一区麻豆| 国产精品久久久久久精品电影| 可以在线观看毛片的网站| 性色av乱码一区二区三区2| 在线观看免费视频日本深夜| 午夜激情福利司机影院| 最新中文字幕久久久久| 啦啦啦观看免费观看视频高清| 国产精品女同一区二区软件 | 欧美精品啪啪一区二区三区| 欧美乱妇无乱码| 中文字幕av成人在线电影| 女警被强在线播放| 久久精品人妻少妇| 一进一出抽搐动态| 亚洲中文字幕一区二区三区有码在线看| 中亚洲国语对白在线视频| 欧美日韩精品网址| 色播亚洲综合网| 网址你懂的国产日韩在线| 丝袜美腿在线中文| 少妇高潮的动态图| 欧美zozozo另类| 少妇熟女aⅴ在线视频| 一夜夜www| www日本在线高清视频| 久久6这里有精品| 老司机在亚洲福利影院| av视频在线观看入口| 婷婷亚洲欧美| 免费人成视频x8x8入口观看| 老司机在亚洲福利影院| 久久天躁狠狠躁夜夜2o2o| 亚洲精品一卡2卡三卡4卡5卡| 欧美成人a在线观看| 在线观看66精品国产| 国模一区二区三区四区视频| 国产视频一区二区在线看| 桃红色精品国产亚洲av| 久久精品91蜜桃| 午夜精品在线福利| 91久久精品电影网| 非洲黑人性xxxx精品又粗又长| 午夜福利在线在线| 久久久国产成人精品二区| 国产欧美日韩一区二区三| 9191精品国产免费久久| www.www免费av| 国产男靠女视频免费网站| 制服人妻中文乱码| 99国产综合亚洲精品| 亚洲成人免费电影在线观看| 首页视频小说图片口味搜索| 国产高清视频在线观看网站| 在线十欧美十亚洲十日本专区| 国产精品嫩草影院av在线观看 | www日本在线高清视频| 麻豆成人午夜福利视频| 一a级毛片在线观看| 免费观看精品视频网站| 久久久久久久精品吃奶| 在线观看av片永久免费下载| 日本在线视频免费播放| 国产一区二区三区在线臀色熟女| 在线观看免费视频日本深夜| 欧美黄色淫秽网站| 亚洲熟妇中文字幕五十中出| 国产精品乱码一区二三区的特点| 一级作爱视频免费观看| 国产高清激情床上av| 在线国产一区二区在线| www.www免费av| 亚洲成人精品中文字幕电影| 啦啦啦观看免费观看视频高清| 亚洲欧美日韩卡通动漫| 大型黄色视频在线免费观看| 级片在线观看| 午夜亚洲福利在线播放| 国产高潮美女av| 亚洲不卡免费看| 综合色av麻豆| 亚洲国产高清在线一区二区三| 少妇的丰满在线观看| 久久精品影院6| 精品福利观看| 在线免费观看不下载黄p国产 | 亚洲片人在线观看| 欧美高清成人免费视频www| 日韩高清综合在线| 少妇人妻一区二区三区视频| 国产精品电影一区二区三区| 色精品久久人妻99蜜桃| 亚洲乱码一区二区免费版| 夜夜躁狠狠躁天天躁| 日韩av在线大香蕉| 久久欧美精品欧美久久欧美| 两个人看的免费小视频| av视频在线观看入口| 欧美+亚洲+日韩+国产| 欧美中文综合在线视频| 乱人视频在线观看| 麻豆久久精品国产亚洲av| 亚洲国产精品久久男人天堂| 欧美性感艳星| 99久久久亚洲精品蜜臀av| 丰满人妻熟妇乱又伦精品不卡| 国内毛片毛片毛片毛片毛片| 成人av一区二区三区在线看| 国产精品永久免费网站| 俄罗斯特黄特色一大片| 99久久久亚洲精品蜜臀av| 高清在线国产一区| 国产精品女同一区二区软件 | 亚洲 欧美 日韩 在线 免费| 老司机午夜福利在线观看视频| 国产精品日韩av在线免费观看| 午夜两性在线视频| 欧美日韩亚洲国产一区二区在线观看| 亚洲av电影在线进入| 国产一区二区三区视频了| 日韩欧美在线二视频| 我要搜黄色片| 亚洲专区中文字幕在线| 观看免费一级毛片| 日韩av在线大香蕉| 日本与韩国留学比较| 啦啦啦观看免费观看视频高清| 久久精品91无色码中文字幕| 欧美av亚洲av综合av国产av|