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

    萬家寨水庫上游的冰情特征分析及預(yù)報(bào)

    2017-03-27 10:14:41熊運(yùn)阜
    水土保持通報(bào) 2017年1期
    關(guān)鍵詞:萬家寨開河冰情

    張 璐, 張 生, 李 超, 熊運(yùn)阜

    (1.內(nèi)蒙古農(nóng)業(yè)大學(xué) 水利與土木建筑工程學(xué)院, 內(nèi)蒙古 呼和浩特市 010018;2.黃河萬家寨水利樞紐有限公司, 山西 太原 036412)

    萬家寨水庫上游的冰情特征分析及預(yù)報(bào)

    張 璐1, 張 生1, 李 超1, 熊運(yùn)阜2

    (1.內(nèi)蒙古農(nóng)業(yè)大學(xué) 水利與土木建筑工程學(xué)院, 內(nèi)蒙古 呼和浩特市 010018;2.黃河萬家寨水利樞紐有限公司, 山西 太原 036412)

    [目的] 研究黃河萬家寨水庫上游至托克托縣喇嘛灣河段冰情特征及預(yù)報(bào)方法,為防凌實(shí)踐提供重要依據(jù)。[方法] 通過對1998—2015年度冰情分析,研究了該河段的封、開河過程及二者之間的聯(lián)系,以野外觀測為基礎(chǔ),分別對封、開河時間歷時與氣溫、流量、冰厚、封河天數(shù)進(jìn)行了相關(guān)性分析,選取了與封、開河時間相關(guān)性系數(shù)較高的影響因子,利用BP神經(jīng)網(wǎng)絡(luò)模型對萬家寨上游河段的封河歷時、開河歷時、開河日期進(jìn)行了預(yù)報(bào)。[結(jié)果] 通過與實(shí)際觀測日期對比,該預(yù)報(bào)方法合格率達(dá)到93%。[結(jié)論] 神經(jīng)網(wǎng)絡(luò)模型對于黃河萬家寨水庫上游段封河歷時、開河歷時、開河日期預(yù)報(bào)的適用性較強(qiáng),更適合解決受復(fù)雜因素影響的問題。

    冰凌; BP神經(jīng)網(wǎng)絡(luò); 封河; 開河; 萬家寨水庫

    文獻(xiàn)參數(shù): 張璐, 張生, 李超, 等.萬家寨水庫上游的冰情特征分析及預(yù)報(bào)[J].水土保持通報(bào),2017,37(1):196-200.DOI:10.13961/j.cnki.stbctb.2017.01.035; Zhang Lu, Zhang Sheng, Li Chao, et al. Analysis and prediction of ice characteristics in upstream of Wajiazhai dam on Yellow River[J]. Bulletin of Soil and Water Conservation, 2017,37(1):196-200.DOI:10.13961/j.cnki.stbctb.2017.01.035

    黃河萬家寨水利樞紐于1998年建成蓄水,在下閘蓄水以前,壩址上游喇嘛灣(距壩76.6 km)以下河段整個冬季常年都以淌凌為主,這主要是由于該段河槽比降大,上游卡冰形成冰蓋后,使得該河段上游來冰量減少,因而該河段由于水氣熱交換而形成的冰凌可以在此順利下泄。自1998年水利樞紐運(yùn)用后,抬升了上游河道水位,減緩了水流流速,使得冰凌在水庫回水末端冰蓋形成后,上游來冰在庫尾上游彎道處堆積,導(dǎo)致庫區(qū)由之前淌凌為主河段變?yōu)榉€(wěn)定封凍河段。其中,在距壩址58 km處的牛龍灣,冰凌容易形成堆積,導(dǎo)致上游63 km處的水泥廠水位壅高,因此,在開河期該處幾乎每年都會產(chǎn)生凌汛災(zāi)害,尤其2007—2008年度堆冰嚴(yán)重,最高水位達(dá)982.24 m,水面高出路面0.8~1.0 m,曹家灣耕地全部被淹沒,造成較大損失。因此,開展封開河冰情預(yù)報(bào)研究可對防凌實(shí)踐提供重要技術(shù)支持。

    河道冰情分析和研究多采用實(shí)地監(jiān)測和模型模擬的方法,在國外,利用熱交換原理和冰力學(xué)原理建立的模型已經(jīng)在冰情預(yù)報(bào)中被廣泛應(yīng)用[1-4],這一方法已用于加拿大圣勞倫斯河蒙特利爾河段,預(yù)測精度較高。在國內(nèi),黃河冰凌問題一直以來備受研究領(lǐng)域和管理部門的關(guān)注,黃河下游冰情預(yù)報(bào)模型自20世紀(jì)80年代末建立后[5]。模糊優(yōu)選神經(jīng)網(wǎng)絡(luò)方法成功地引入到黃河內(nèi)蒙段冰情預(yù)報(bào)[6-9]中,取得了較好的效果,許可誤差為6 d的情況下,預(yù)報(bào)精度達(dá)86.7%。而對于本文的研究區(qū)域的相關(guān)研究主要為冰塞、冰壩等冰情災(zāi)害的研究[10-11]。本文擬在此基礎(chǔ)上,選取1998—1999年至2014—2015年共17個冬季(當(dāng)年11月至次年3月)的實(shí)測數(shù)據(jù),分別對黃河萬家寨上游的溫度、流量、水位及河道地形變化進(jìn)行分析,以此來分析萬家寨上游冰情變化特點(diǎn),并選取關(guān)聯(lián)度大的預(yù)報(bào)因子,運(yùn)用BP神經(jīng)網(wǎng)絡(luò)模型,對萬家寨上游封河歷時、開河歷時、開河日期進(jìn)行預(yù)報(bào)分析,為防凌實(shí)踐提供重要依據(jù)。

    1 區(qū)域及研究方法

    1.1 研究區(qū)概況

    研究區(qū)域?yàn)閺娜f家寨水庫壩址至頭道拐河段,總長112.4 km。該河段屬于峽谷型河道,庫區(qū)河道呈“U”型,河寬約300~500 m,從頭道拐至拐上河段屬于天然河道,河面開闊且縱坡較緩,拐上是河道縱坡由緩變陡的轉(zhuǎn)折點(diǎn),拐上以下河道的比降達(dá)1.17‰。整個研究區(qū)域處于黃河內(nèi)蒙段下游,由于冬季受到蒙古冷高壓的控制,氣候寒冷且干燥,雨雪較少[12-14],冰期長4—5個月,凌汛期12月與1月氣溫最低,多年平均氣溫分別為-8,-10.6 ℃,最低氣溫分別為-24.8,-23.3 ℃。萬家寨庫區(qū)一般年在11月開始流凌,12月進(jìn)入穩(wěn)封期,3月解凍開河,多年穩(wěn)封期平均流量為387 m3/s,開河期平均流量為885 m3/s。距壩58 km處的牛龍灣為“S”型彎道,此處的河床比降與河面寬度變化較大,特殊的地形條件及橋墩、護(hù)堤等人工建筑物的阻擋,使得冰凌在此處不易下泄,極易卡冰結(jié)壩,造成嚴(yán)重的冰凌災(zāi)害。

    1.2 人工神經(jīng)網(wǎng)絡(luò)

    人工神經(jīng)網(wǎng)絡(luò)[15-16]是由大量神經(jīng)元組成的極其復(fù)雜的非線性網(wǎng)絡(luò)系統(tǒng)。其過程如下[17]:設(shè)輸入神經(jīng)元為h,隱含層神經(jīng)元為i,輸出神經(jīng)元為j,nh,ni,nj為3層神經(jīng)網(wǎng)絡(luò)的節(jié)點(diǎn)數(shù)目,θi,θj為隱含層節(jié)點(diǎn)i,輸出層節(jié)點(diǎn)j的閾值,whi,wij為輸入層節(jié)點(diǎn)h與隱含層節(jié)點(diǎn)i間、隱含層節(jié)點(diǎn)i和輸出層節(jié)點(diǎn)j間的連線的權(quán)值,各節(jié)點(diǎn)的輸入為x,輸出為y。

    (1) 初始化。設(shè)已歸一化的輸入、輸出樣本為

    {xk,h,dk,j|k=1,2,…,nk;h=1,2,…,nh;j=1,2,…,nj}

    nk為樣本容量。給各連接權(quán){whi},{wij} 和閾值{θi},{θj} 賦予(-0.1,0.1)區(qū)間上的隨機(jī)值。

    (2) 置k=1,把樣本對(xk,h,dk,j)提供給網(wǎng)絡(luò)(h=1,2,…,nh;j=1,2,…,nj)。

    (3) 計(jì)算隱含層各節(jié)點(diǎn)的輸入xi,輸出yi(i=1,2,…,ni)。

    yi=1/(1+e-xi)

    (4) 計(jì)算輸出層各節(jié)點(diǎn)的輸入xj,輸出yj(i=1,2,…,nj)

    (5) 計(jì)算輸出層各節(jié)點(diǎn)所收到的總輸入變化時單樣本點(diǎn)誤差Ek的變化率。

    (6) 計(jì)算隱含層各節(jié)點(diǎn)所收到的總輸入變化時單樣本點(diǎn)誤差的變化率。

    (7) 修正各鏈接的權(quán)值和閾值。

    式中:t——修正次數(shù),學(xué)習(xí)速率η∈(0,1),動量因子α∈(0,1)。η較大,則算法收斂快,但是不穩(wěn)定,可能出現(xiàn)震蕩,η較小則算法收斂緩慢;α的作用恰好與η相反。

    (8) 置k=k+1,取學(xué)習(xí)模式對(xk,h,dk,j)提供給網(wǎng)絡(luò),轉(zhuǎn)步驟3,直至全部nk個模式對訓(xùn)練完畢,轉(zhuǎn)步驟9。

    (9) 重復(fù)步驟2至步驟8,直到網(wǎng)絡(luò)全局誤差函數(shù)。

    小于預(yù)先設(shè)定的一個較小值或?qū)W習(xí)次數(shù)大于預(yù)先設(shè)定的值,結(jié)束學(xué)習(xí)。

    當(dāng)信號輸入到神經(jīng)網(wǎng)絡(luò)的時候,首先傳到隱含層節(jié)點(diǎn),經(jīng)過網(wǎng)絡(luò)作用函數(shù)轉(zhuǎn)化后,再把隱含層節(jié)點(diǎn)的輸出信號傳播到輸出層節(jié)點(diǎn),經(jīng)過處理后輸出結(jié)果。其中,節(jié)點(diǎn)的作用函數(shù)通常采用Sigmoid型函數(shù),其數(shù)學(xué)模型是:f(x)=1/(1+e-x)。

    1.3 模型參數(shù)設(shè)置

    在BP神經(jīng)網(wǎng)絡(luò)系統(tǒng)中,設(shè)定含有輸入層(接收外界輸入的信息,并傳遞給中間層各神經(jīng)元)、隱含層(處理內(nèi)部信息,負(fù)責(zé)信息變換)和輸出層(向外界進(jìn)行信息輸出)的3層網(wǎng)絡(luò)結(jié)構(gòu)。根據(jù)預(yù)報(bào)因子數(shù)不同輸入層的輸入節(jié)點(diǎn)數(shù)也不同,預(yù)報(bào)封河歷時選擇了3個影響因子,因此輸入節(jié)點(diǎn)數(shù)為3,預(yù)報(bào)開河歷選擇了4個影響因子,因此輸入節(jié)點(diǎn)數(shù)為4,預(yù)報(bào)開河日期選擇了6個影響因子,因此節(jié)點(diǎn)數(shù)為6;經(jīng)過反復(fù)數(shù)值訓(xùn)練表明,隱含層設(shè)5個隱節(jié)點(diǎn)比較合適;輸出層僅有一個單結(jié)點(diǎn)輸出即為預(yù)報(bào)對象。參考文獻(xiàn)[7]設(shè)置BP神經(jīng)網(wǎng)絡(luò)的運(yùn)行參數(shù)為:學(xué)習(xí)效率為0.9,動量系數(shù)為0.05,訓(xùn)練次數(shù)一般為2 000~100 000,經(jīng)過模型的反復(fù)訓(xùn)練,得出訓(xùn)練次數(shù)為10 000時預(yù)報(bào)精度達(dá)到最高。根據(jù)上述的模型算法,用1998—2010年資料對網(wǎng)絡(luò)進(jìn)行反復(fù)訓(xùn)練,用2010—2015年的數(shù)據(jù)進(jìn)行驗(yàn)證,直到滿足精度要求為止。

    1.4 預(yù)報(bào)因子的篩選

    影響凌汛的2大因素是自然因素和人為因素[18]。其中自然因素包括熱力因素、動力因素、河勢因素。熱力因素主要包括水溫、氣溫、降水等;動力因素主要有流量、流速、水位、風(fēng)速等;河勢因素主要包括河道比降、河面寬度、糙率、彎曲度等。人為因素主要包括人工建筑物和人類活動的影響。通過野外實(shí)地觀測,發(fā)現(xiàn)入冬時氣溫持續(xù)低于-5 ℃以下開始流凌,冰凌首先會在庫尾回水末端或彎道處卡冰堆積,溫度持續(xù)降低,封河距離增加,所以流凌封河時的主控因素為溫度,而日均氣溫對封河的影響不大,所以選取了11月至封河時的累積負(fù)氣溫,17a來累積負(fù)氣溫平均為-177.45 ℃,同時選取了負(fù)氣溫的持續(xù)天數(shù)及11月25至12月6日的流量,凌期流量平均在554.80m3/s左右。進(jìn)入開河期,隨著氣溫的回升,冰蓋變薄且脆性增大,冰下過流能力增強(qiáng),導(dǎo)致冰蓋破裂而開河,因此自然因素為開河的主要影響因素,本文選取了1月9至26日流量均值及最大冰厚作為影響因子,因?yàn)榉夂悠诹髁康拇笮『捅駥﹂_河有重要影響,近年來平均冰厚為0.98m,累積正氣溫與正氣溫累積天數(shù)也是重要影響因子。本文未考慮河勢因素,原因是研究時段內(nèi)河道每年的擺動比較小,地形變化不大,所以河勢因素相對變化比較小,故未予考慮。冰情數(shù)據(jù)源于萬家寨水庫上游岔河口冰情站沿線各測站的統(tǒng)計(jì)資料,基本能夠反映萬家寨水庫上游的冰情現(xiàn)象,且精度可靠。綜上所述,基于野外觀測與理論分析,選取了不同日期的溫度、流量、冰厚、封河天數(shù)及封開河時的流量與溫度,分別確定各影響因子與相應(yīng)預(yù)報(bào)對象的相關(guān)性系數(shù),如表1為各影響因子與各預(yù)報(bào)對象的相關(guān)性分析表。

    表1 各預(yù)報(bào)因子與封河歷時、開河歷時、開河日期的相關(guān)系數(shù)

    2 萬家寨水利樞紐冰情變化分析

    萬家寨水庫建立后,根據(jù)野外觀測發(fā)現(xiàn)該河段冰情主要變化特點(diǎn)為:喇嘛灣以下穩(wěn)定封凍河段增長;冰塞、冰壩生成幾率增大;水位擁高,造成冰凌災(zāi)害[19]。由黃河內(nèi)蒙古段下游典型斷面水位變化如圖1(所用的數(shù)據(jù)為2013—2014年度冰期實(shí)測數(shù)據(jù)),從圖1可以看出,水位從頭道拐、蒲灘拐、喇嘛灣、水泥廠、岔河口沿程降低,從頭道拐至岔河口河段水位降低15 m左右,幅度較大,由此也可以得出該河段比降較大??梢悦黠@看出封河時岔河口斷面首先封河,然后朔源而上各斷面依次封河,開河時則相反。其中,頭道拐斷面整個冰期水位變化不大,封開河時水位也沒有明顯漲落,說明該斷面地勢平緩,封開河較平穩(wěn)。蒲灘拐斷面水位在封開河時有較小的波動,但整體比較平穩(wěn),沒有太大的變化。喇嘛灣斷面水位在封開河時較蒲灘拐變化明顯,因?yàn)樵摂嗝嫣幱趶澋郎嫌危瑥澋揽ū?,極易壅高此處水位。水泥廠斷面的水位漲落最明顯,也是最易發(fā)生凌汛災(zāi)害的斷面,該斷面下游為“S”型彎道,11月23日時氣溫由正轉(zhuǎn)負(fù),到26日時降至11月最低氣溫-7.5 ℃,氣溫的突降導(dǎo)致下游卡冰堆積而壅高該斷面水位,之后進(jìn)入穩(wěn)封期,3月17日氣溫升高至7.9 ℃且持續(xù)上升時,水位第2次上漲,說明下游河段因?yàn)樯嫌魏佣谓鈨鲩_河的冰凌堆積而產(chǎn)生冰壩,壅高了上游水位。岔河口斷面水位漲落較明顯但是比較緩慢,該斷面地勢較低,流凌易在此堆積,但不會造成大的凌汛災(zāi)害。

    圖1 萬家寨水庫上游水位變化曲線

    繪制冬季冰情縱斷面圖如圖2所示。從圖2可以看出,整個河段明顯分為4部分。壩址-32 km河段因?yàn)槭芩畮熘苯佑绊?,水流流速非常小,雖然該段河床比降非常大,但是水面比降小,所以氣溫由正轉(zhuǎn)負(fù)后,該河段便開始結(jié)冰封河,屬于平封河段。距壩32~52.1 km河段主要受到溫度的影響,而導(dǎo)致封河形態(tài)的位置不同,11月氣溫偏低,立封河段會增加,氣溫偏高,平封河段會增加,屬于平立封交替河段。距壩52.1~67.6 km河段有明顯的高程落差,而且

    河床比降和水面比降都很大,水庫未建庫蓄水以前,該河段因?yàn)榈靥幱趰{谷型河段,河道比降、水流流速較大而從不封凍,建庫后,由于受到回水末端的影響,使得該河段流速減小,水位抬升,再加上大壩對流凌的攔截作用,導(dǎo)致冰凌極易在彎道處堆積,封河時此河段主要受動力因素影響而形成動力型冰蓋,即該河段為立封河段。距壩67.6 km以上河段,河槽比降小,封河時主要受熱力因素影響而形成熱力型冰蓋,即該河段為平封河段;整個河段封河時朔源而上,11月氣溫由正轉(zhuǎn)負(fù)后,由于壩前的流速幾乎為零,壩前至32 km河段開始結(jié)冰,且冰面平滑,氣溫持續(xù)降低,便開始流凌,冰凌在平封邊緣處開始堆積,受地形影響,到距壩67.6 km后,才由立封轉(zhuǎn)為平封河段。開河時上游首先開河,開河形成的流冰堆積到52.1~67.6 km河段內(nèi),極易生成冰壩,形成武開河局勢,擁高水位造成冰凌災(zāi)害。

    圖2 萬家寨水庫上游冬季冰情縱斷面

    3 模型的預(yù)測

    根據(jù)《水文情報(bào)預(yù)報(bào)規(guī)范》[20]中有關(guān)春訊和冰情預(yù)報(bào)的規(guī)定,并參照多年的防凌減災(zāi)實(shí)踐工作,預(yù)見期定為16 d,根據(jù)線性內(nèi)插法,得到許可誤差為6 d。根據(jù)表2預(yù)報(bào)檢驗(yàn)結(jié)果可知,封河歷時和開河歷時的預(yù)報(bào)結(jié)果均在許可誤差范圍內(nèi),且和實(shí)測值非常接近,而開河日期的預(yù)報(bào)除2011—2012年度的預(yù)報(bào)誤差峰值為9 d大于許可誤差外,其余的預(yù)報(bào)值都在許可誤差內(nèi)。

    表2 萬家寨上游封河歷時、開河歷時、開河日期預(yù)報(bào)模型結(jié)果

    分析其原因(如圖3所示),發(fā)現(xiàn)該年度雖然封河晚、冰層薄,但是封凍時間較長,該年度進(jìn)入3月以后氣溫偏低,氣溫遲遲達(dá)不到開河的條件,3月15日之前主要以融冰為主,并且河段內(nèi)發(fā)生許多新增的清溝與原有的清溝逐漸加寬加長且連通,因此導(dǎo)致開河過程減緩,較其它年份來說開河較晚,所以該年度的預(yù)測結(jié)果偏差較大。綜上所述,共15個預(yù)報(bào)結(jié)果中,預(yù)報(bào)合格率為93.3%,預(yù)報(bào)精度較高,滿足水文預(yù)報(bào)的精度要求。

    圖3 2010—2015年及3月份多年平均溫度對比

    4 結(jié) 論

    (1) 從所選取的預(yù)報(bào)因子來看,氣溫、封凍天數(shù)、冰厚和流量是影響冰情變化的主控因素。

    (2) 河段的封河形式主要取決于水面比降的大小。

    (3) 頭道拐以下庫尾河段總是先于頭道拐以上河段解凍開河,且形成庫尾交通水泥廠附近河段冰壩的冰源來自頭道拐以下,主要是距壩67~100 km河段的冰量。

    (4) 由于黃河封、開河受多種因素的綜合影響,且各影響因子雖單獨(dú)與預(yù)報(bào)因子呈線性關(guān)系,但綜合影響下關(guān)系較復(fù)雜,所以神經(jīng)網(wǎng)絡(luò)模型對于黃河萬家寨水庫上游段封河歷時、開河歷時、開河日期預(yù)報(bào)的適用性較強(qiáng),更適合解決上述受復(fù)雜因素影響的問題。經(jīng)過神經(jīng)網(wǎng)絡(luò)反復(fù)訓(xùn)練,最終得出預(yù)報(bào)合格率為93.3%,預(yù)報(bào)精度較高,滿足水文預(yù)報(bào)的精度要求。

    [1] Shen Hungtao, Yapa P D. A unified degree-day method for river ice cover thickness simulation[J]. Canadian Journal of Civil Engineering, 1985,12(1):54-62.

    [2] Foltyn E P, Shen H T. St. Lawrence River freeze-up forecast[J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 1986,112(4):467-481.

    [3] Shen Hungtao, Shen H, Tsai S M. Dynamic transport of river ice[J]. Journal of Hydraulic Research, 1990,28(6): 659-671.

    [4] Shen Hungtao. Under cover transport and accumulation of frazil granules[J]. Journal of Hydraulic Engineering,1995,121(2):184-195.

    [5] 冀鴻蘭,朝倫巴根,陳守煜.基于遺傳算法的模糊優(yōu)選BP網(wǎng)絡(luò)的冰凌預(yù)報(bào)模型[J].中國農(nóng)村水利水電,2009(1):5-7.

    [6] 陳守煜,冀鴻蘭.冰凌預(yù)報(bào)模糊優(yōu)選神經(jīng)網(wǎng)絡(luò)BP方法[J].水利學(xué)報(bào),2004(6):114-118.

    [7] 冀鴻蘭,張傲妲,高瑞忠,等.黃河內(nèi)蒙古段開河日期預(yù)報(bào)模型及應(yīng)用[J].水利水電科技進(jìn)展,2012,32(2):42-45.

    [8] 李鳳玲,高瑞忠,馮國華,等.兩種改進(jìn)的BP模型預(yù)報(bào)黃河內(nèi)蒙段冰情信息及效果分析[J].中國農(nóng)村水利水電,2011(7):69-72.

    [9] 冀鴻蘭,朝倫巴根,陳守煜.模糊優(yōu)選神經(jīng)網(wǎng)絡(luò)BP模型在黃河內(nèi)蒙段封開河預(yù)報(bào)中的應(yīng)用[J].水利水電科技進(jìn)展,2008,28(3):70-72.

    [10] 鄭春茂,裴群,閆新光,等.萬家寨水庫運(yùn)用對上下游河段防凌的影響[J].人民黃河,2002,24(3):34-35.

    [11] 可素娟,錢云平,郝守英,等.萬家寨水庫凌情特點(diǎn)及初始封河期調(diào)度方案研究[J].冰川凍土,2002,24(2):209-213.

    [12] 顧潤源,周偉燦,白美蘭,等.氣候變化對黃河內(nèi)蒙古段凌汛期的影響[J].中國沙漠,2012,32(6):1751-1756.

    [13] 丁永建,葉佰生,劉時銀.祁連山中部地區(qū)40 a來氣候變化及其對徑流的影響[J].冰川凍土,2002,22(3):193-199.

    [14] 王根緒,沈永平,程國棟.黃河河源區(qū)生態(tài)環(huán)境變化與成因分析[J].冰川凍土,2000,22(3):200-205.

    [15] 吳昌友.神經(jīng)網(wǎng)絡(luò)的研究及應(yīng)用[D].哈爾濱:東北林業(yè)大學(xué),2007:22-34.

    [16] 宋玉強(qiáng).人工神經(jīng)網(wǎng)絡(luò)在時間序列預(yù)測中的應(yīng)用研究[D].西安:西安建筑科技大學(xué),2005:14-27.

    [17] 金龍.神經(jīng)網(wǎng)絡(luò)氣象預(yù)報(bào)建模理論方法與應(yīng)用[M].北京:氣象出版社,2004:75-80.

    [18] 姚惠明,秦福興,沈國昌,等.黃河寧蒙河段凌情特性研究[J].水科學(xué)進(jìn)展,2007,18(6):893-899.

    [19] 馬喜祥,熊運(yùn)阜,徐偉.萬家寨水庫冰情淺議[J].泥沙研究,2003(4):66-72.

    [20] 水利部水利信息中心.水文情報(bào)預(yù)報(bào)規(guī)范(GB/T22482-2008)[S].北京:中國水利水電出版社,2008.

    Analysis and Prediction of Ice Characteristics in Upstream of Wajiazhai Dam on Yellow River

    ZHANG Lu1, ZHANG Sheng1, LI Chao1, XIONG Yunfu2

    (1.CollegeofWaterResourcesandCivilEngineering,InnerMongoliaAgriculturalUniversity,Hohhot,InnerMongolia010018,China; 2.YellowRiverWanjiazhaiWaterMulti-purposeDamProjectCo.,Ltd.,Taiyuan,Shanxi036412,China)

    [Objective] The objective of this paper is to explore the characteristics and forecasting methods of ice regime of upper reach of Wanjiazhai dam and Lamawan reach in Togtoh County after Wanjiazhai reservoir, and provide a basis for ice prevention and scheduling. [Methods] The processes of freezing-up and thawing of the river reach were analyzed based on the field observations from 1998 to 2015, including the dates and durations of both freezing-up and thawing and their relationships to the affecting factors like air temperature, flow rate, ice thickness, and number of days of frozen. The significant correlation coefficients were selected and BP neural network model was applied to predict the durations of the freezing-up, durations of thawing, and dates of thawing of the upstream reach respectively. [Results] The accuracy of the predicting results were up to 93% compared with the observed date. [Conclusion] Neural network model is suitable for forecasting freezing-up duration, breaking-up duration and breaking-date in upper reach of Wanjiazhai dam, and is more suitable for solving the problem of ice regime influenced by complex factors.

    ice; BP neural network approach; freeze-up; thawing; Wanjiazhai reservoir

    2016-03-15

    2016-06-10

    國家自然科學(xué)基金項(xiàng)目“凍融過程中湖泊污染物多介質(zhì)遷移轉(zhuǎn)化規(guī)律及機(jī)制研究”(51339002), “寒區(qū)湖泊冰封期營養(yǎng)物質(zhì)冰水多介質(zhì)環(huán)境過程及對富營養(yǎng)化影響效應(yīng)研究”(51369017)

    張璐(1990—),女(漢族),內(nèi)蒙古自治區(qū)巴彥淖爾市人,碩士研究生,研究方向?yàn)辄S河河冰研究。E-mail:1026093195@qq.com。

    張生(1960—),男(漢族),內(nèi)蒙古自治區(qū)包頭市,教授,主要從事水環(huán)境科學(xué)與工程研究。E-mail:shengzhang@imau.edu.cn。

    B

    1000-288X(2017)01-0196-05

    P338+.4

    猜你喜歡
    萬家寨開河冰情
    松花江干流哈爾濱江段封凍
    中國水利(2020年23期)2020-12-24 06:24:40
    成都院中標(biāo)智慧萬家寨建設(shè)規(guī)劃項(xiàng)目
    黃河萬家寨庫區(qū)水質(zhì)監(jiān)測評價(jià)及保護(hù)措施簡介
    詠菊致友人
    萬家寨引黃入晉工程清水及消防供水控制系統(tǒng)改造設(shè)計(jì)
    禪意
    黃河之聲(2016年24期)2016-04-22 02:39:44
    小開河流淌出一條美麗生態(tài)帶
    走向世界(2016年1期)2016-04-13 06:04:34
    破冰開河
    支點(diǎn)(2016年3期)2016-03-21 13:01:12
    海冰基礎(chǔ)知識及船舶冰區(qū)航行的注意事項(xiàng)
    黃河萬家寨庫區(qū) 水質(zhì)狀況分析
    国产日韩一区二区三区精品不卡| 两个人免费观看高清视频| 国产一卡二卡三卡精品| 亚洲av美国av| 大陆偷拍与自拍| 亚洲黑人精品在线| 视频区图区小说| 操美女的视频在线观看| 人妻久久中文字幕网| 国产人伦9x9x在线观看| 桃花免费在线播放| 久久影院123| 亚洲欧洲精品一区二区精品久久久| 国产成人a∨麻豆精品| 亚洲精品国产av成人精品| 亚洲va日本ⅴa欧美va伊人久久 | 国产精品一区二区精品视频观看| 飞空精品影院首页| 日本五十路高清| 国产不卡av网站在线观看| 亚洲第一av免费看| 国产精品久久久人人做人人爽| 纯流量卡能插随身wifi吗| 在线亚洲精品国产二区图片欧美| 12—13女人毛片做爰片一| 日本91视频免费播放| 考比视频在线观看| 在线天堂中文资源库| 亚洲性夜色夜夜综合| 18禁国产床啪视频网站| 亚洲av日韩精品久久久久久密| 亚洲久久久国产精品| 国产精品.久久久| 伊人亚洲综合成人网| 极品人妻少妇av视频| 正在播放国产对白刺激| 久久国产亚洲av麻豆专区| 在线 av 中文字幕| 久久久精品区二区三区| 亚洲少妇的诱惑av| 老汉色∧v一级毛片| 午夜福利影视在线免费观看| 高清黄色对白视频在线免费看| 久久久久国内视频| 麻豆国产av国片精品| 夜夜夜夜夜久久久久| 大片电影免费在线观看免费| 看免费av毛片| 免费日韩欧美在线观看| 国产精品免费视频内射| 韩国精品一区二区三区| 国产真人三级小视频在线观看| 三上悠亚av全集在线观看| 老汉色∧v一级毛片| 丁香六月天网| 国产欧美日韩一区二区三 | 国产精品九九99| 51午夜福利影视在线观看| 亚洲美女黄色视频免费看| 国产成人免费观看mmmm| 欧美中文综合在线视频| 国产99久久九九免费精品| 久久精品亚洲熟妇少妇任你| 另类精品久久| 午夜福利乱码中文字幕| videosex国产| 国产成人一区二区三区免费视频网站| 欧美 亚洲 国产 日韩一| 国产在线一区二区三区精| 搡老岳熟女国产| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲天堂av无毛| 最新在线观看一区二区三区| 亚洲成人手机| 久久人人爽人人片av| 久久国产精品影院| 欧美精品啪啪一区二区三区 | 丝袜喷水一区| 中文字幕人妻丝袜一区二区| 亚洲一码二码三码区别大吗| 久久久国产成人免费| 日本撒尿小便嘘嘘汇集6| svipshipincom国产片| av在线播放精品| 亚洲va日本ⅴa欧美va伊人久久 | 99re6热这里在线精品视频| 久久人妻熟女aⅴ| 美女主播在线视频| 高清av免费在线| 99久久综合免费| 亚洲精品第二区| 免费黄频网站在线观看国产| 岛国在线观看网站| 国产精品 欧美亚洲| 啦啦啦在线免费观看视频4| 在线观看舔阴道视频| 一级黄色大片毛片| 新久久久久国产一级毛片| 一个人免费看片子| 老汉色av国产亚洲站长工具| 最黄视频免费看| 啦啦啦在线免费观看视频4| 日韩欧美国产一区二区入口| 国产精品.久久久| 色婷婷av一区二区三区视频| 黄色视频在线播放观看不卡| 99国产精品免费福利视频| 亚洲一卡2卡3卡4卡5卡精品中文| 精品久久久精品久久久| 两个人免费观看高清视频| 在线观看免费日韩欧美大片| 在线观看舔阴道视频| 午夜影院在线不卡| 免费高清在线观看视频在线观看| 亚洲成人手机| 久热这里只有精品99| 最新的欧美精品一区二区| 黄色怎么调成土黄色| 久久人人97超碰香蕉20202| 国产精品一区二区在线观看99| 日韩大片免费观看网站| 国产精品1区2区在线观看. | 无限看片的www在线观看| 国产淫语在线视频| 国产色视频综合| 男男h啪啪无遮挡| 少妇精品久久久久久久| a在线观看视频网站| 高潮久久久久久久久久久不卡| 自线自在国产av| 久久人妻熟女aⅴ| 免费av中文字幕在线| 在线观看www视频免费| 日韩熟女老妇一区二区性免费视频| 日本wwww免费看| 亚洲自偷自拍图片 自拍| 久久精品aⅴ一区二区三区四区| 欧美日韩成人在线一区二区| 一本综合久久免费| 久久久欧美国产精品| 亚洲熟女精品中文字幕| 国产精品久久久久成人av| 久久精品aⅴ一区二区三区四区| 中文欧美无线码| 91精品国产国语对白视频| tube8黄色片| 母亲3免费完整高清在线观看| 色视频在线一区二区三区| 国产深夜福利视频在线观看| 国产福利在线免费观看视频| 一区二区三区四区激情视频| 成年动漫av网址| 成人手机av| 在线观看免费高清a一片| 免费黄频网站在线观看国产| 青草久久国产| 国产精品久久久久成人av| 亚洲视频免费观看视频| 国产国语露脸激情在线看| 老司机深夜福利视频在线观看 | 丁香六月欧美| 精品人妻在线不人妻| 香蕉国产在线看| 欧美日韩国产mv在线观看视频| 成年动漫av网址| 久9热在线精品视频| 午夜福利影视在线免费观看| 十八禁网站网址无遮挡| 人妻久久中文字幕网| 精品国产超薄肉色丝袜足j| 91成年电影在线观看| 亚洲综合色网址| 在线观看www视频免费| 久久国产精品大桥未久av| 国产精品免费大片| 午夜精品国产一区二区电影| 男女边摸边吃奶| 免费在线观看日本一区| 成年动漫av网址| 69精品国产乱码久久久| 欧美精品亚洲一区二区| 国产精品自产拍在线观看55亚洲 | 在线av久久热| 丝袜在线中文字幕| 国产成人精品无人区| 黄片播放在线免费| 日韩 欧美 亚洲 中文字幕| 久久九九热精品免费| netflix在线观看网站| 亚洲精品久久成人aⅴ小说| 国产av又大| 一本久久精品| 岛国毛片在线播放| 男女高潮啪啪啪动态图| 国产野战对白在线观看| 波多野结衣一区麻豆| 国产精品熟女久久久久浪| 亚洲成av片中文字幕在线观看| 久久久久久亚洲精品国产蜜桃av| 每晚都被弄得嗷嗷叫到高潮| 狠狠婷婷综合久久久久久88av| 正在播放国产对白刺激| 天天影视国产精品| 菩萨蛮人人尽说江南好唐韦庄| 老司机在亚洲福利影院| 成在线人永久免费视频| 亚洲 欧美一区二区三区| 国产xxxxx性猛交| 久9热在线精品视频| 老司机亚洲免费影院| 中文字幕人妻丝袜一区二区| 一级毛片电影观看| 日本91视频免费播放| 中文欧美无线码| 国产精品久久久久成人av| 亚洲久久久国产精品| 精品久久蜜臀av无| 国产在线一区二区三区精| 1024视频免费在线观看| 国产精品久久久久成人av| 91精品三级在线观看| 亚洲综合色网址| 午夜精品国产一区二区电影| 国产色视频综合| 欧美性长视频在线观看| 午夜福利在线免费观看网站| 久久国产精品人妻蜜桃| 精品国产一区二区三区四区第35| 日韩制服骚丝袜av| 久久久久久久大尺度免费视频| 满18在线观看网站| 色视频在线一区二区三区| 久久久国产欧美日韩av| 免费在线观看黄色视频的| 欧美日韩av久久| 一个人免费在线观看的高清视频 | 日韩视频在线欧美| 日韩制服骚丝袜av| 天堂俺去俺来也www色官网| 免费观看av网站的网址| 丝袜美足系列| 91精品三级在线观看| 美女扒开内裤让男人捅视频| 涩涩av久久男人的天堂| 国产成人av教育| 亚洲精品成人av观看孕妇| 国内毛片毛片毛片毛片毛片| 国产成人免费无遮挡视频| 日韩欧美一区视频在线观看| 亚洲第一欧美日韩一区二区三区 | 国产在线一区二区三区精| 黄片大片在线免费观看| 欧美精品亚洲一区二区| 满18在线观看网站| 欧美人与性动交α欧美软件| 老司机影院毛片| 久久久久久久久免费视频了| 欧美变态另类bdsm刘玥| 亚洲av男天堂| 国产精品影院久久| 视频区图区小说| 建设人人有责人人尽责人人享有的| 少妇裸体淫交视频免费看高清 | 欧美黑人欧美精品刺激| 一级片'在线观看视频| av线在线观看网站| 精品国产乱码久久久久久小说| 国产亚洲av高清不卡| 亚洲精品中文字幕一二三四区 | 亚洲av成人不卡在线观看播放网 | 精品亚洲乱码少妇综合久久| 美女主播在线视频| 欧美日韩黄片免| 最近最新免费中文字幕在线| 女性被躁到高潮视频| 成人免费观看视频高清| 淫妇啪啪啪对白视频 | 首页视频小说图片口味搜索| 十八禁高潮呻吟视频| 国产av又大| 国产精品免费视频内射| 母亲3免费完整高清在线观看| av电影中文网址| 亚洲精品粉嫩美女一区| 久久综合国产亚洲精品| 国产不卡av网站在线观看| 女人精品久久久久毛片| 亚洲精品第二区| www.av在线官网国产| 中文字幕人妻熟女乱码| 亚洲自偷自拍图片 自拍| 蜜桃在线观看..| 在线观看免费午夜福利视频| 动漫黄色视频在线观看| 91国产中文字幕| 水蜜桃什么品种好| 亚洲黑人精品在线| 女性被躁到高潮视频| 美女脱内裤让男人舔精品视频| 亚洲精品久久午夜乱码| 在线精品无人区一区二区三| 两个人看的免费小视频| 亚洲精品一二三| 欧美中文综合在线视频| 欧美黑人精品巨大| 免费久久久久久久精品成人欧美视频| 欧美在线黄色| 午夜激情久久久久久久| 国产日韩欧美在线精品| 国产成人欧美在线观看 | 亚洲精品在线美女| 一区二区三区精品91| 国产精品一区二区精品视频观看| 亚洲国产欧美一区二区综合| 手机成人av网站| 香蕉丝袜av| 91精品伊人久久大香线蕉| 欧美日本中文国产一区发布| 国产av国产精品国产| 午夜福利,免费看| 欧美少妇被猛烈插入视频| 亚洲欧美色中文字幕在线| 操美女的视频在线观看| 91成人精品电影| 国产成人一区二区三区免费视频网站| 丝瓜视频免费看黄片| 久久久国产成人免费| 亚洲精品在线美女| 亚洲天堂av无毛| 亚洲中文字幕日韩| 桃红色精品国产亚洲av| 丰满饥渴人妻一区二区三| 久久久久国内视频| 亚洲成人免费电影在线观看| 久久九九热精品免费| 亚洲熟女毛片儿| 亚洲av国产av综合av卡| 亚洲中文字幕日韩| 成年人黄色毛片网站| 天堂俺去俺来也www色官网| 又大又爽又粗| 建设人人有责人人尽责人人享有的| 1024视频免费在线观看| 亚洲成人免费电影在线观看| 国产精品成人在线| 亚洲国产欧美日韩在线播放| av福利片在线| 日韩制服丝袜自拍偷拍| 两人在一起打扑克的视频| 嫩草影视91久久| bbb黄色大片| 亚洲精品国产av蜜桃| 最黄视频免费看| 成人av一区二区三区在线看 | 免费黄频网站在线观看国产| 成年女人毛片免费观看观看9 | 久久精品亚洲熟妇少妇任你| 欧美在线一区亚洲| 少妇 在线观看| 亚洲精华国产精华精| av免费在线观看网站| 欧美精品一区二区大全| 日韩 亚洲 欧美在线| 日韩人妻精品一区2区三区| 九色亚洲精品在线播放| 这个男人来自地球电影免费观看| 精品国产超薄肉色丝袜足j| 国产精品成人在线| 中文字幕制服av| 国产欧美日韩综合在线一区二区| 亚洲精品一区蜜桃| 狂野欧美激情性xxxx| 国产精品麻豆人妻色哟哟久久| 成人国产av品久久久| 欧美老熟妇乱子伦牲交| 青青草视频在线视频观看| 宅男免费午夜| 中国国产av一级| 人人妻人人爽人人添夜夜欢视频| 亚洲精品av麻豆狂野| 亚洲精品国产精品久久久不卡| 丁香六月欧美| 国产av精品麻豆| 欧美乱码精品一区二区三区| 在线看a的网站| 手机成人av网站| 国产精品 国内视频| 侵犯人妻中文字幕一二三四区| 最近最新免费中文字幕在线| 久久久精品区二区三区| 欧美日韩福利视频一区二区| 国产深夜福利视频在线观看| 爱豆传媒免费全集在线观看| 成年动漫av网址| 男女国产视频网站| 十八禁网站免费在线| 精品少妇久久久久久888优播| 日韩,欧美,国产一区二区三区| 国产99久久九九免费精品| 日本黄色日本黄色录像| 午夜福利视频精品| 欧美激情高清一区二区三区| 精品福利观看| 桃红色精品国产亚洲av| 午夜精品国产一区二区电影| 久久ye,这里只有精品| 亚洲精品一二三| 精品一区二区三区av网在线观看 | 成年美女黄网站色视频大全免费| 下体分泌物呈黄色| 91国产中文字幕| 欧美日韩亚洲高清精品| 久久久精品国产亚洲av高清涩受| 国产精品久久久av美女十八| 一级片'在线观看视频| 久久精品国产亚洲av香蕉五月 | 国产亚洲欧美精品永久| 日韩欧美一区视频在线观看| 50天的宝宝边吃奶边哭怎么回事| 美女高潮到喷水免费观看| 99久久综合免费| 肉色欧美久久久久久久蜜桃| 欧美日韩av久久| 18在线观看网站| 丝袜在线中文字幕| 亚洲av男天堂| 日韩欧美一区视频在线观看| 午夜免费鲁丝| 国产一区二区三区在线臀色熟女 | 成年人黄色毛片网站| 纯流量卡能插随身wifi吗| 三级毛片av免费| 亚洲精品一区蜜桃| 亚洲第一欧美日韩一区二区三区 | 国产精品.久久久| 国产亚洲精品一区二区www | 高清黄色对白视频在线免费看| 亚洲精品一卡2卡三卡4卡5卡 | 97精品久久久久久久久久精品| 亚洲激情五月婷婷啪啪| 女人被躁到高潮嗷嗷叫费观| 首页视频小说图片口味搜索| 国产精品久久久人人做人人爽| 黄片播放在线免费| 老司机影院毛片| 成年人午夜在线观看视频| 国产成人影院久久av| 久久久久久人人人人人| 久久久久久亚洲精品国产蜜桃av| 考比视频在线观看| 国产成人系列免费观看| 国产视频一区二区在线看| 丝袜美腿诱惑在线| 一级毛片精品| 午夜福利视频在线观看免费| 汤姆久久久久久久影院中文字幕| 免费日韩欧美在线观看| 国产男女超爽视频在线观看| 美女福利国产在线| 亚洲精品美女久久av网站| 欧美大码av| 久久精品国产亚洲av香蕉五月 | 久久人人97超碰香蕉20202| 91国产中文字幕| 少妇 在线观看| 亚洲黑人精品在线| 亚洲午夜精品一区,二区,三区| 51午夜福利影视在线观看| 一本大道久久a久久精品| 亚洲人成77777在线视频| 啦啦啦在线免费观看视频4| 国产亚洲欧美在线一区二区| 王馨瑶露胸无遮挡在线观看| 国产三级黄色录像| 国产高清国产精品国产三级| 免费在线观看完整版高清| 法律面前人人平等表现在哪些方面 | 免费在线观看视频国产中文字幕亚洲 | 亚洲av欧美aⅴ国产| 国产成+人综合+亚洲专区| 99国产综合亚洲精品| 王馨瑶露胸无遮挡在线观看| 成人免费观看视频高清| 啦啦啦在线免费观看视频4| 老汉色av国产亚洲站长工具| 9色porny在线观看| 欧美黑人欧美精品刺激| 亚洲av电影在线观看一区二区三区| 国产主播在线观看一区二区| 美女扒开内裤让男人捅视频| 日韩电影二区| 婷婷丁香在线五月| 国产又爽黄色视频| 亚洲国产欧美在线一区| 搡老岳熟女国产| 国产男人的电影天堂91| tube8黄色片| 爱豆传媒免费全集在线观看| 日本91视频免费播放| 亚洲免费av在线视频| 免费黄频网站在线观看国产| 各种免费的搞黄视频| 天天影视国产精品| 老司机福利观看| h视频一区二区三区| svipshipincom国产片| 国产欧美亚洲国产| 成年人午夜在线观看视频| 亚洲色图 男人天堂 中文字幕| 日本五十路高清| 免费av中文字幕在线| 亚洲中文av在线| 欧美激情久久久久久爽电影 | 窝窝影院91人妻| 国产免费福利视频在线观看| 精品国产超薄肉色丝袜足j| 国产成人免费观看mmmm| a级片在线免费高清观看视频| 中文字幕高清在线视频| 国产精品久久久久久人妻精品电影 | 亚洲精品一区蜜桃| 伦理电影免费视频| 大陆偷拍与自拍| 久热爱精品视频在线9| bbb黄色大片| 久久精品aⅴ一区二区三区四区| 丁香六月欧美| 母亲3免费完整高清在线观看| 纯流量卡能插随身wifi吗| 2018国产大陆天天弄谢| 亚洲国产日韩一区二区| 女性被躁到高潮视频| 久久香蕉激情| 男人添女人高潮全过程视频| 看免费av毛片| 亚洲一码二码三码区别大吗| 国产欧美日韩一区二区三区在线| 成人三级做爰电影| 国产一区二区激情短视频 | 久久香蕉激情| avwww免费| 国产成人a∨麻豆精品| 亚洲av欧美aⅴ国产| 成年av动漫网址| 国产精品免费大片| 成人亚洲精品一区在线观看| 精品少妇一区二区三区视频日本电影| 国产成人一区二区三区免费视频网站| 日韩欧美一区二区三区在线观看 | 在线av久久热| 99re6热这里在线精品视频| 在线十欧美十亚洲十日本专区| 久9热在线精品视频| 国产淫语在线视频| 男男h啪啪无遮挡| 久久久久国内视频| 丝袜喷水一区| 精品国内亚洲2022精品成人 | 久久精品aⅴ一区二区三区四区| 真人做人爱边吃奶动态| 中国美女看黄片| 午夜免费成人在线视频| 久久影院123| 国内毛片毛片毛片毛片毛片| 国产一级毛片在线| 亚洲国产中文字幕在线视频| 女性生殖器流出的白浆| 一边摸一边做爽爽视频免费| 一区二区av电影网| 久久亚洲国产成人精品v| 老司机深夜福利视频在线观看 | 日本欧美视频一区| 啦啦啦中文免费视频观看日本| 欧美精品亚洲一区二区| 色婷婷av一区二区三区视频| 人妻一区二区av| 亚洲自偷自拍图片 自拍| 麻豆国产av国片精品| 欧美黑人精品巨大| 91九色精品人成在线观看| 国产一区二区三区av在线| 精品视频人人做人人爽| 亚洲国产欧美网| 免费高清在线观看视频在线观看| 在线观看一区二区三区激情| 啦啦啦中文免费视频观看日本| a级毛片在线看网站| 国产野战对白在线观看| 久久ye,这里只有精品| 狠狠婷婷综合久久久久久88av| 老熟妇仑乱视频hdxx| 国产精品欧美亚洲77777| 亚洲性夜色夜夜综合| 欧美亚洲日本最大视频资源| 国产精品免费大片| 国产av精品麻豆| 亚洲国产欧美日韩在线播放| 一级毛片电影观看| 国产亚洲精品久久久久5区| 在线亚洲精品国产二区图片欧美| 亚洲欧洲精品一区二区精品久久久| 中文字幕人妻丝袜制服| 亚洲精品国产av蜜桃| 中国国产av一级| 亚洲一码二码三码区别大吗| 中文字幕制服av| 欧美大码av| 99国产精品一区二区蜜桃av | 午夜精品国产一区二区电影| 久久ye,这里只有精品| 狠狠狠狠99中文字幕| 久久精品国产a三级三级三级| 天堂8中文在线网| 日韩人妻精品一区2区三区| 免费黄频网站在线观看国产|