• <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ì)狀況分析
    一级,二级,三级黄色视频| xxx大片免费视频| 午夜日本视频在线| 日韩人妻精品一区2区三区| 肉色欧美久久久久久久蜜桃| 精品国产露脸久久av麻豆| 色网站视频免费| 色网站视频免费| 老鸭窝网址在线观看| 久久精品亚洲av国产电影网| 亚洲欧美成人综合另类久久久| 国产一区亚洲一区在线观看| 国产人伦9x9x在线观看 | 99国产精品免费福利视频| 亚洲av日韩在线播放| 欧美日韩视频精品一区| 久久综合国产亚洲精品| 亚洲精品成人av观看孕妇| 天天躁狠狠躁夜夜躁狠狠躁| 叶爱在线成人免费视频播放| 日韩一卡2卡3卡4卡2021年| 黄色怎么调成土黄色| 久久人人爽人人片av| 国产精品香港三级国产av潘金莲 | 亚洲精品美女久久av网站| 视频在线观看一区二区三区| 免费av中文字幕在线| 国语对白做爰xxxⅹ性视频网站| 午夜精品国产一区二区电影| 99re6热这里在线精品视频| freevideosex欧美| 国产人伦9x9x在线观看 | 精品国产乱码久久久久久小说| 青春草亚洲视频在线观看| 国产成人aa在线观看| 美女午夜性视频免费| 午夜福利在线免费观看网站| 精品一区二区三区四区五区乱码 | 亚洲美女视频黄频| 男男h啪啪无遮挡| 18禁动态无遮挡网站| 亚洲精华国产精华液的使用体验| 婷婷色麻豆天堂久久| 亚洲欧洲日产国产| 最近的中文字幕免费完整| 一个人免费看片子| 在线精品无人区一区二区三| 欧美最新免费一区二区三区| 国产精品二区激情视频| 欧美日韩综合久久久久久| 国产女主播在线喷水免费视频网站| 中文欧美无线码| 亚洲一级一片aⅴ在线观看| 男女高潮啪啪啪动态图| 一二三四在线观看免费中文在| 男女无遮挡免费网站观看| 国产有黄有色有爽视频| 久久久久视频综合| av电影中文网址| 男女午夜视频在线观看| 欧美最新免费一区二区三区| 中国三级夫妇交换| 免费在线观看视频国产中文字幕亚洲 | 精品一区二区三卡| 男的添女的下面高潮视频| 午夜激情久久久久久久| av女优亚洲男人天堂| 免费看av在线观看网站| 亚洲国产精品一区三区| av片东京热男人的天堂| 成人黄色视频免费在线看| 日本vs欧美在线观看视频| 亚洲综合色网址| 亚洲欧洲日产国产| 亚洲色图综合在线观看| 黄色 视频免费看| 国产毛片在线视频| 69精品国产乱码久久久| 少妇猛男粗大的猛烈进出视频| 久久久久久久国产电影| 91精品三级在线观看| 亚洲成人av在线免费| 日韩一本色道免费dvd| av又黄又爽大尺度在线免费看| 国产成人免费无遮挡视频| 国产成人a∨麻豆精品| 有码 亚洲区| 色吧在线观看| 日本vs欧美在线观看视频| 一区在线观看完整版| 最近中文字幕高清免费大全6| 国产精品蜜桃在线观看| 亚洲人成77777在线视频| 女人精品久久久久毛片| 亚洲一区中文字幕在线| 久久久欧美国产精品| 男的添女的下面高潮视频| 91久久精品国产一区二区三区| 免费黄色在线免费观看| 欧美人与善性xxx| 波多野结衣一区麻豆| 啦啦啦中文免费视频观看日本| 中文字幕制服av| 看十八女毛片水多多多| 少妇的丰满在线观看| 国产精品久久久久久精品电影小说| 成年女人在线观看亚洲视频| 日产精品乱码卡一卡2卡三| 精品第一国产精品| 三上悠亚av全集在线观看| 日韩熟女老妇一区二区性免费视频| 亚洲av成人精品一二三区| 大片免费播放器 马上看| 女人高潮潮喷娇喘18禁视频| 亚洲一码二码三码区别大吗| 色哟哟·www| 精品久久久精品久久久| 欧美成人午夜精品| 国产无遮挡羞羞视频在线观看| 美女国产高潮福利片在线看| 久久人人爽av亚洲精品天堂| 久久青草综合色| a 毛片基地| 涩涩av久久男人的天堂| 精品人妻熟女毛片av久久网站| av网站免费在线观看视频| 精品一区在线观看国产| 一区二区三区激情视频| videosex国产| 国产激情久久老熟女| 亚洲欧美一区二区三区国产| 国产精品不卡视频一区二区| 精品一区二区三区四区五区乱码 | 日韩人妻精品一区2区三区| 久久热在线av| 久久精品国产亚洲av天美| 美女脱内裤让男人舔精品视频| 亚洲国产av影院在线观看| 春色校园在线视频观看| 欧美激情 高清一区二区三区| 亚洲熟女精品中文字幕| av卡一久久| 亚洲精品中文字幕在线视频| 国产麻豆69| 黄频高清免费视频| 国产成人精品久久久久久| 最新中文字幕久久久久| 午夜久久久在线观看| 国语对白做爰xxxⅹ性视频网站| av在线老鸭窝| 叶爱在线成人免费视频播放| 欧美精品一区二区免费开放| 国产亚洲欧美精品永久| 9色porny在线观看| 一级a爱视频在线免费观看| 亚洲欧美精品综合一区二区三区 | 在线观看人妻少妇| 亚洲精品美女久久久久99蜜臀 | 免费黄频网站在线观看国产| 成人免费观看视频高清| 亚洲色图 男人天堂 中文字幕| 一本—道久久a久久精品蜜桃钙片| 亚洲精品久久午夜乱码| 久久久久久久大尺度免费视频| 国产伦理片在线播放av一区| 电影成人av| 777米奇影视久久| 色视频在线一区二区三区| 国精品久久久久久国模美| 亚洲五月色婷婷综合| 久久狼人影院| av网站免费在线观看视频| 日韩一本色道免费dvd| √禁漫天堂资源中文www| 久久热在线av| 人妻系列 视频| 蜜桃在线观看..| 国产av精品麻豆| 交换朋友夫妻互换小说| 日韩中字成人| 大片电影免费在线观看免费| 亚洲国产欧美网| 亚洲成人一二三区av| 免费少妇av软件| 国产精品一区二区在线观看99| 免费高清在线观看视频在线观看| 午夜日韩欧美国产| 99久久精品国产国产毛片| 十八禁高潮呻吟视频| 亚洲人成77777在线视频| 欧美亚洲 丝袜 人妻 在线| 大片电影免费在线观看免费| 夜夜骑夜夜射夜夜干| 精品亚洲成国产av| 国产欧美亚洲国产| 亚洲激情五月婷婷啪啪| 桃花免费在线播放| 两个人看的免费小视频| 色婷婷久久久亚洲欧美| 一级爰片在线观看| 香蕉丝袜av| 观看av在线不卡| 成人国产av品久久久| 亚洲,欧美精品.| 永久网站在线| 可以免费在线观看a视频的电影网站 | 欧美最新免费一区二区三区| 国产精品偷伦视频观看了| 久久热在线av| 热re99久久精品国产66热6| 国产欧美亚洲国产| 国产成人91sexporn| 天美传媒精品一区二区| 亚洲欧美一区二区三区国产| 考比视频在线观看| 日韩一本色道免费dvd| 成人漫画全彩无遮挡| 国产精品99久久99久久久不卡 | 亚洲精品中文字幕在线视频| 久久av网站| 久久久亚洲精品成人影院| 99久久中文字幕三级久久日本| 一二三四中文在线观看免费高清| 亚洲一区中文字幕在线| 国产一级毛片在线| 国产成人一区二区在线| 国产男女超爽视频在线观看| 最黄视频免费看| 丰满迷人的少妇在线观看| 亚洲精品成人av观看孕妇| 免费黄频网站在线观看国产| 最近中文字幕2019免费版| 免费不卡的大黄色大毛片视频在线观看| 国产1区2区3区精品| 午夜精品国产一区二区电影| 女性生殖器流出的白浆| 又大又黄又爽视频免费| 精品少妇久久久久久888优播| 妹子高潮喷水视频| 国产xxxxx性猛交| 久久精品国产亚洲av高清一级| 成人国产av品久久久| 亚洲欧美精品自产自拍| 亚洲av男天堂| videossex国产| 另类精品久久| 人妻 亚洲 视频| 久久精品亚洲av国产电影网| 久久影院123| 亚洲av综合色区一区| 五月开心婷婷网| 少妇人妻 视频| 精品卡一卡二卡四卡免费| 国产精品久久久久久精品古装| 免费看不卡的av| 日本av免费视频播放| 久久久久久久久免费视频了| 久久人人爽人人片av| 人妻一区二区av| 中国三级夫妇交换| 三级国产精品片| 看十八女毛片水多多多| 免费在线观看黄色视频的| 日日爽夜夜爽网站| 女的被弄到高潮叫床怎么办| 十八禁网站网址无遮挡| 最新中文字幕久久久久| 2021少妇久久久久久久久久久| 欧美在线黄色| 天天躁夜夜躁狠狠久久av| 有码 亚洲区| 激情视频va一区二区三区| 赤兔流量卡办理| 一级片免费观看大全| www.精华液| 亚洲精品在线美女| 国产午夜精品一二区理论片| 久久99精品国语久久久| 精品少妇久久久久久888优播| 丁香六月天网| 日韩一本色道免费dvd| 亚洲精品国产色婷婷电影| 中文字幕制服av| 欧美人与善性xxx| 日本午夜av视频| 亚洲国产精品一区三区| 亚洲国产精品国产精品| 高清黄色对白视频在线免费看| 日韩大片免费观看网站| a 毛片基地| 精品人妻熟女毛片av久久网站| 视频区图区小说| 最近手机中文字幕大全| 精品视频人人做人人爽| 亚洲av中文av极速乱| 日韩成人av中文字幕在线观看| 亚洲av在线观看美女高潮| 熟女电影av网| 久久久久人妻精品一区果冻| 乱人伦中国视频| 欧美日韩综合久久久久久| av免费在线看不卡| 夜夜骑夜夜射夜夜干| 亚洲成av片中文字幕在线观看 | 精品少妇内射三级| 亚洲国产看品久久| 亚洲精品自拍成人| 免费观看无遮挡的男女| 18+在线观看网站| 一级片免费观看大全| 亚洲精品国产色婷婷电影| 999久久久国产精品视频| 狠狠婷婷综合久久久久久88av| 中文精品一卡2卡3卡4更新| 在线免费观看不下载黄p国产| 大香蕉久久网| 久久久久视频综合| 秋霞在线观看毛片| 国产在线视频一区二区| 99热国产这里只有精品6| 久久韩国三级中文字幕| av在线老鸭窝| 日韩av在线免费看完整版不卡| 国产精品一区二区在线不卡| 亚洲少妇的诱惑av| 亚洲一区二区三区欧美精品| 热re99久久国产66热| 一级爰片在线观看| 国产成人精品久久二区二区91 | 午夜av观看不卡| 一级毛片我不卡| 国产不卡av网站在线观看| 亚洲国产精品一区二区三区在线| 黑人欧美特级aaaaaa片| videos熟女内射| 人成视频在线观看免费观看| 美女xxoo啪啪120秒动态图| 纯流量卡能插随身wifi吗| 国产精品一区二区在线观看99| 极品人妻少妇av视频| 亚洲av免费高清在线观看| 午夜福利一区二区在线看| 美女主播在线视频| 超碰成人久久| 日韩在线高清观看一区二区三区| 观看av在线不卡| 不卡av一区二区三区| 亚洲,欧美,日韩| 十分钟在线观看高清视频www| av女优亚洲男人天堂| 亚洲,欧美,日韩| 亚洲成国产人片在线观看| 水蜜桃什么品种好| 夫妻性生交免费视频一级片| 十分钟在线观看高清视频www| 热99国产精品久久久久久7| 日本欧美视频一区| 欧美精品人与动牲交sv欧美| 在线观看美女被高潮喷水网站| 高清在线视频一区二区三区| 一级片免费观看大全| 国产1区2区3区精品| 久久国产精品男人的天堂亚洲| 国产精品亚洲av一区麻豆 | 熟女电影av网| 老司机亚洲免费影院| 亚洲欧美精品自产自拍| 亚洲av免费高清在线观看| 亚洲欧美精品自产自拍| 欧美日韩av久久| 电影成人av| 亚洲欧美一区二区三区国产| 久久99精品国语久久久| 国产精品久久久久久久久免| 中文字幕人妻熟女乱码| 国产爽快片一区二区三区| 日韩中字成人| 亚洲伊人色综图| 搡女人真爽免费视频火全软件| 18禁动态无遮挡网站| 少妇被粗大猛烈的视频| 精品99又大又爽又粗少妇毛片| 免费播放大片免费观看视频在线观看| 十八禁高潮呻吟视频| 国产成人91sexporn| 国产亚洲精品第一综合不卡| 欧美精品国产亚洲| 国产精品国产av在线观看| 久久午夜综合久久蜜桃| 捣出白浆h1v1| 久久韩国三级中文字幕| 熟女电影av网| 9热在线视频观看99| 日韩av在线免费看完整版不卡| 亚洲精品成人av观看孕妇| 欧美日韩成人在线一区二区| 欧美成人精品欧美一级黄| 99香蕉大伊视频| 中文字幕制服av| 少妇人妻 视频| 精品人妻一区二区三区麻豆| 中文精品一卡2卡3卡4更新| 国产精品香港三级国产av潘金莲 | 最近最新中文字幕免费大全7| 妹子高潮喷水视频| 伊人亚洲综合成人网| 夜夜骑夜夜射夜夜干| 成人免费观看视频高清| 欧美精品国产亚洲| 香蕉精品网在线| 在线看a的网站| 亚洲,一卡二卡三卡| 人妻 亚洲 视频| 99热网站在线观看| 婷婷色av中文字幕| 欧美日韩精品成人综合77777| 熟妇人妻不卡中文字幕| 亚洲av男天堂| 亚洲精品久久成人aⅴ小说| 亚洲欧洲精品一区二区精品久久久 | 老熟女久久久| 中文字幕色久视频| 精品人妻在线不人妻| 在线观看国产h片| 精品少妇黑人巨大在线播放| 久久午夜福利片| 欧美xxⅹ黑人| 麻豆乱淫一区二区| 国产极品粉嫩免费观看在线| 一个人免费看片子| av一本久久久久| 免费在线观看完整版高清| 亚洲色图综合在线观看| 香蕉丝袜av| 国产精品久久久久久精品电影小说| 亚洲av国产av综合av卡| 国产精品久久久久久久久免| av片东京热男人的天堂| 一区福利在线观看| 人人妻人人添人人爽欧美一区卜| 免费高清在线观看日韩| 老司机影院成人| 婷婷色麻豆天堂久久| 美女视频免费永久观看网站| 七月丁香在线播放| av视频免费观看在线观看| 国产高清国产精品国产三级| 女人久久www免费人成看片| 咕卡用的链子| 久久久亚洲精品成人影院| 婷婷色av中文字幕| 久久人人97超碰香蕉20202| 欧美日韩一级在线毛片| tube8黄色片| 中文字幕色久视频| 久久亚洲国产成人精品v| 高清不卡的av网站| 欧美激情高清一区二区三区 | 久久这里有精品视频免费| 伦理电影免费视频| 国产欧美亚洲国产| 不卡av一区二区三区| 2018国产大陆天天弄谢| 亚洲三级黄色毛片| 亚洲av福利一区| av在线老鸭窝| 制服人妻中文乱码| 成人免费观看视频高清| 午夜91福利影院| 男的添女的下面高潮视频| 热re99久久精品国产66热6| 亚洲美女视频黄频| 国产色婷婷99| 国产一区二区 视频在线| 看十八女毛片水多多多| 亚洲国产色片| 国产熟女午夜一区二区三区| 国产精品麻豆人妻色哟哟久久| 91久久精品国产一区二区三区| 尾随美女入室| 91精品伊人久久大香线蕉| 高清黄色对白视频在线免费看| 日韩中文字幕欧美一区二区 | 国产精品成人在线| 黑人巨大精品欧美一区二区蜜桃| 中文精品一卡2卡3卡4更新| 免费在线观看完整版高清| 嫩草影院入口| 久久免费观看电影| 欧美日韩亚洲高清精品| 自拍欧美九色日韩亚洲蝌蚪91| 免费在线观看完整版高清| 好男人视频免费观看在线| 久久99一区二区三区| 亚洲三区欧美一区| 一区二区av电影网| 人妻 亚洲 视频| 日日爽夜夜爽网站| 纵有疾风起免费观看全集完整版| 妹子高潮喷水视频| 18禁裸乳无遮挡动漫免费视频| 人妻 亚洲 视频| 日日摸夜夜添夜夜爱| 亚洲av日韩在线播放| 国产精品国产三级国产专区5o| 一区二区av电影网| 欧美变态另类bdsm刘玥| 亚洲成人手机| 电影成人av| 18禁观看日本| 熟女少妇亚洲综合色aaa.| 不卡视频在线观看欧美| 99久久精品国产国产毛片| 大香蕉久久成人网| 久久久a久久爽久久v久久| 在线观看国产h片| 国产男女超爽视频在线观看| 亚洲五月色婷婷综合| 久久这里只有精品19| 80岁老熟妇乱子伦牲交| 99久久综合免费| 久久免费观看电影| 久久精品国产亚洲av天美| 97在线视频观看| 久久久久久久久久久免费av| 曰老女人黄片| 最黄视频免费看| 久久狼人影院| 久久国产亚洲av麻豆专区| 伊人亚洲综合成人网| 国产爽快片一区二区三区| 日本猛色少妇xxxxx猛交久久| 九草在线视频观看| 亚洲一区中文字幕在线| 久久精品久久久久久噜噜老黄| 一区福利在线观看| 男男h啪啪无遮挡| 91午夜精品亚洲一区二区三区| 国产深夜福利视频在线观看| 精品一区在线观看国产| 色网站视频免费| 青春草国产在线视频| 成人国产av品久久久| 成人毛片60女人毛片免费| 18禁动态无遮挡网站| 亚洲精品国产一区二区精华液| 国产成人精品福利久久| 国产精品99久久99久久久不卡 | 精品国产国语对白av| 男女啪啪激烈高潮av片| 亚洲少妇的诱惑av| 两个人看的免费小视频| 制服丝袜香蕉在线| 久久久久久人人人人人| 欧美少妇被猛烈插入视频| 亚洲人成77777在线视频| 麻豆乱淫一区二区| 爱豆传媒免费全集在线观看| 国产一区二区三区综合在线观看| 五月天丁香电影| 观看av在线不卡| 午夜福利乱码中文字幕| 王馨瑶露胸无遮挡在线观看| 亚洲激情五月婷婷啪啪| 国产亚洲av片在线观看秒播厂| 国产成人av激情在线播放| 久久人妻熟女aⅴ| 久久久久久久大尺度免费视频| 欧美中文综合在线视频| 性高湖久久久久久久久免费观看| 国产黄色视频一区二区在线观看| 久久久精品国产亚洲av高清涩受| 18禁国产床啪视频网站| 国语对白做爰xxxⅹ性视频网站| 日本av免费视频播放| 男女无遮挡免费网站观看| 毛片一级片免费看久久久久| 亚洲人成电影观看| 久久久a久久爽久久v久久| 久久精品久久精品一区二区三区| 男女下面插进去视频免费观看| 中文字幕亚洲精品专区| 午夜久久久在线观看| 黑丝袜美女国产一区| 超碰97精品在线观看| 又粗又硬又长又爽又黄的视频| 看免费成人av毛片| 亚洲五月色婷婷综合| 可以免费在线观看a视频的电影网站 | 99精国产麻豆久久婷婷| 国产无遮挡羞羞视频在线观看| 2018国产大陆天天弄谢| 久久国产精品大桥未久av| 国产无遮挡羞羞视频在线观看| 亚洲内射少妇av| 色网站视频免费| 国产成人欧美| av一本久久久久| 亚洲第一区二区三区不卡| 欧美成人精品欧美一级黄| 激情五月婷婷亚洲| 成年美女黄网站色视频大全免费| 天天躁夜夜躁狠狠躁躁| 丰满迷人的少妇在线观看| 国产免费福利视频在线观看| 亚洲美女搞黄在线观看| 久久久久久久久免费视频了| 成人毛片a级毛片在线播放| av在线app专区| 我的亚洲天堂| 一级,二级,三级黄色视频| 欧美av亚洲av综合av国产av | 伊人久久国产一区二区| 国产精品女同一区二区软件| 成人毛片a级毛片在线播放|