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

    基于Copula熵-隨機(jī)森林的中長期徑流預(yù)報研究

    2021-12-09 23:12:01黃朝君賈建偉秦赫王棟
    人民長江 2021年11期

    黃朝君 賈建偉 秦赫 王棟

    摘要:預(yù)測因子作為中長期預(yù)報模型的輸入項(xiàng),是影響預(yù)報結(jié)果精度的關(guān)鍵要素。為進(jìn)一步提高預(yù)報精度,提出了一種Copula熵與隨機(jī)森林模型相結(jié)合的中長期徑流預(yù)報方法。該方法首先采用Copula熵指標(biāo)對預(yù)測因子進(jìn)行篩選,然后將選取的預(yù)測因子作為輸入項(xiàng),導(dǎo)入隨機(jī)森林模型中對月徑流進(jìn)行相應(yīng)預(yù)測。將該方法應(yīng)用于漢江流域丹江口水庫的逐月入庫徑流預(yù)報中,并與相關(guān)系數(shù)篩選法進(jìn)行對比。結(jié)果表明:基于Copula熵指標(biāo)篩選出的預(yù)測因子對應(yīng)的模擬結(jié)果具有更高的精度,尤其對于汛期而言,其模擬值與實(shí)測值的擬合優(yōu)度顯著優(yōu)于比選方法,說明其篩選出的預(yù)測因子具有更好的合理性。

    關(guān)鍵詞:中長期徑流預(yù)報; 預(yù)測因子; 大氣環(huán)流因子; Copula熵; 隨機(jī)森林模型; 丹江口水庫

    中圖法分類號: TV211.1+2

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

    DOI:10.16232/j.cnki.1001-4179.2021.11.013

    0引 言

    中長期徑流預(yù)報作為水文預(yù)報的重要組成部分,是指依據(jù)已知信息對預(yù)見期為月、季、年尺度的徑流要素進(jìn)行定性或定量的預(yù)測[1-5]。與短期徑流預(yù)報相比,其預(yù)報結(jié)果可為水利管理部門提供更為充足的決策時間,在水庫調(diào)度、防洪減災(zāi)及水資源綜合管理等工作中發(fā)揮重要作用[6-8]。

    目前,針對中長期徑流預(yù)報已有大量研究,尤其隨著人工智能與數(shù)據(jù)挖掘技術(shù)的快速發(fā)展,大量機(jī)器學(xué)習(xí)算法被應(yīng)用于該領(lǐng)域。如酈于杰等[9]采用支持向量機(jī)對漢江流域皇莊站的長期徑流過程進(jìn)行了預(yù)報,并對預(yù)報結(jié)果的不確定性進(jìn)行了相應(yīng)分析;許斌等[10]以丹江口水庫為例,比較了隨機(jī)森林與梯度提升樹兩種機(jī)器學(xué)習(xí)模型的預(yù)報精度;謝帥等[11]將LASSO回歸與支持向量機(jī)相耦合,并應(yīng)用于龍羊峽水庫入庫徑流預(yù)報研究中;Huang等[12]將多種機(jī)器學(xué)習(xí)算法與BMA方法相結(jié)合,開展了基于多模型耦合的漢江流域中長期徑流預(yù)報研究;鄭炎輝等[13]采用隨機(jī)森林與RBF人工神經(jīng)網(wǎng)絡(luò)模型開展了新豐江水庫枯季入庫徑流中長期預(yù)報研究;王棟等[14]將經(jīng)驗(yàn)?zāi)B(tài)分解方法與BP神經(jīng)網(wǎng)絡(luò)相結(jié)合,對青海省湟水河流域月徑流過程進(jìn)行了預(yù)測;仕玉治等[15]將相關(guān)向量機(jī)、支持向量機(jī)及自動回歸滑動平均模型應(yīng)用于南方兩水庫入庫月徑流中長期預(yù)報研究中,并比較了三者的精度差異。然而,上述研究大都側(cè)重于預(yù)報模型的選取,對預(yù)報因子篩選普遍采用相關(guān)系數(shù)法,只考慮了預(yù)測因子與預(yù)測變量間的線性關(guān)系,缺乏對兩者間非線性相關(guān)性的探索。

    本文提出了一種基于Copula熵與隨機(jī)森林模型相結(jié)合的中長期徑流預(yù)報方法,并將其應(yīng)用于漢江流域丹江口水庫的入庫徑流預(yù)報研究中。該方法首先采用Copula熵指標(biāo)對預(yù)報因子進(jìn)行篩選,在此基礎(chǔ)上采用隨機(jī)森林模型對長預(yù)見期月徑流過程進(jìn)行相應(yīng)預(yù)報。該預(yù)報方法結(jié)果精度更高,可為水利管理部門提供更為可靠的決策依據(jù)。

    1計(jì)算方法

    1.1Copula熵理論

    1.1.1互信息(MI)

    作為描述不同變量間相互聯(lián)系的互信息指標(biāo)源于信息理論學(xué),其反映了不同變量間共有信息的大小,即互信息越大,兩者相關(guān)性越強(qiáng),反之則越弱[16]。假定隨機(jī)變量X與Y間存在一定聯(lián)系,則兩者間互信息可通過公式(1)計(jì)算。

    式中:xi及yi分別為隨機(jī)變量X與Y的樣本,N為樣本容量。

    由式(1)不難發(fā)現(xiàn),當(dāng)X與Y互不相關(guān)時,MI值接近于0,而當(dāng)X與Y呈現(xiàn)函數(shù)關(guān)系時,MI值將接近于正無窮。與其他相似性指標(biāo)相比,互信息指標(biāo)不僅可以反映變量間的非線性相關(guān)性,且隨機(jī)變量的任意可逆變換不改變互信息的大小。

    1.1.2Copula熵

    互信息指標(biāo)具有較強(qiáng)的信息挖掘能力,但實(shí)際研究中不同隨機(jī)變量間的聯(lián)合分布往往是偏態(tài)且非均質(zhì)性的,因此難以找到合適的分布類型對其進(jìn)行擬合。為解決該問題,MA等[17]引入Copula理論,結(jié)合互信息原理提出了Copula熵指標(biāo)。記u=F(x),v=F(y)分別為隨機(jī)變量X與Y的邊緣分布條件,則相應(yīng)Copula熵值可通過公式(2)計(jì)算。

    1.1.3預(yù)報因子篩選

    基于Copula熵理論,在相關(guān)系數(shù)顯著性檢驗(yàn)的基礎(chǔ)上,充分考慮不同變量與預(yù)測值間的非線性關(guān)系,來篩選合理的預(yù)測因子,具體步驟如下:

    (1) 通過理論分析,搜集會對預(yù)測變量產(chǎn)生影響的所有可能變量。

    (2) 計(jì)算所有初選因子與預(yù)測變量間的相關(guān)系數(shù),并對兩者的相關(guān)性進(jìn)行顯著性檢驗(yàn),挑選出顯著性變量作為備選因子。

    (3) 利用公式(2)逐一計(jì)算備選因子與月徑流間的Copula熵值,并計(jì)算各因子對應(yīng)的Hampel值。

    (4) 依據(jù)Hampel準(zhǔn)則,篩選合格的預(yù)測因子;若合格因子超過10個,則僅選取熵值排序前10的預(yù)測因子;若不足10個,則所有變量即為最終預(yù)測因子。

    1.2隨機(jī)森林模型

    隨機(jī)森林是Breiman于2001年提出的一種袋裝法與分類回歸樹(CART)相結(jié)合的并行增強(qiáng)機(jī)器學(xué)習(xí)算法[20]。作為一種白箱模型,隨機(jī)森林通過采用自舉法(Bootstrapping)對原始數(shù)據(jù)進(jìn)行大量采樣,并對各生成樣本分別建立決策樹,通過改變變量取值來判斷決策樹的變化,在此基礎(chǔ)上綜合考慮整個決策森林的變化結(jié)果。最終,該算法通過對所有決策樹的預(yù)測進(jìn)行組合來獲取模型模擬結(jié)果。具體步驟如下:

    (1) 采用自助采樣法對含m個樣本的數(shù)據(jù)集進(jìn)行T輪采樣。

    (2) 對每一輪采樣結(jié)果進(jìn)行訓(xùn)練時引入隨機(jī)屬性選擇,即先從回歸樹的每個節(jié)點(diǎn)的屬性集合(假設(shè)共有d個屬性)中隨機(jī)選擇一個包含k個屬性的子集,再從該子集中根據(jù)最優(yōu)屬性劃分,從而擬合各個分類回歸樹。

    2實(shí)例驗(yàn)證

    2.1研究區(qū)域概況

    漢江是長江中游最大的支流,干流流經(jīng)陜西、湖北兩省,于武漢市注入長江,干流全長1 577 km,流域面積約15.9萬km2。流域地勢西高東低,由西部的中低山區(qū)向東逐漸降至丘陵平原區(qū)。

    丹江口水庫位于漢江干流上游,是南水北調(diào)中線工程的水源地,具有防洪、供水、發(fā)電、灌溉、航運(yùn)、養(yǎng)殖等綜合功能。水庫以上流域面積約9.52萬km2,占漢江流域的60%,多為高山丘陵區(qū),屬北亞熱帶季風(fēng)氣候區(qū),年降水量700~1 800 mm,年徑流深300~900 mm。降水是徑流的主要來源,兩者年內(nèi)分配不均勻、年際變化較大。

    研究選取丹江口水庫1956~2016年天然徑流全年12個月數(shù)據(jù),對各月流量過程進(jìn)行模擬。丹江口水庫天然入庫徑流由長江水利委員會水文局提供:丹江口建庫前1956~1966年采用黃家港實(shí)測資料;建庫后的1967~2016年,基于黃家港實(shí)測資料,考慮上游水庫調(diào)蓄、引水量、庫面蒸發(fā)、上游工農(nóng)業(yè)及生活耗水等因素,根據(jù)水量平衡原理進(jìn)行還原。

    2.2預(yù)測因子篩選

    考慮到各月流量過程影響成因存在一定差異,研究對丹江口水庫入庫天然徑流全年12個月分別篩選預(yù)測因子,以選取的預(yù)測因子作為模型輸入項(xiàng),對各月流量過程進(jìn)行模擬。

    已有研究表明,大氣環(huán)流與長期水文過程間存在密切關(guān)系,其異?,F(xiàn)象將導(dǎo)致水文過程發(fā)生不同程度的變化。因此,本次研究將國家氣象局氣候中心提供的130項(xiàng)環(huán)流指數(shù)納入預(yù)測因子的初選范疇??紤]大氣環(huán)流因子與徑流過程間的遙相關(guān)性,以徑流過程發(fā)生前一年各月變量作為初選因子。其中,為驗(yàn)證基于Copula熵理論的預(yù)測因子篩選方法,研究設(shè)置了相應(yīng)比選方案,該方案采用相關(guān)系數(shù)指標(biāo)進(jìn)行篩選,挑選出相關(guān)系數(shù)最大的10個變量作為預(yù)報因子。

    以汛期7月平均流量為例,兩種方案選定的預(yù)報因子如表1所列,其中方案1為比選方案,方案2為基于Copula熵理論的預(yù)報因子篩選方案。由表1可以看出,對于7月而言,2種方案所選因子主要差異在于第9項(xiàng)與第10項(xiàng)。基于Copula熵理論的方案2挑選了南海副高脊線位置指數(shù)和印度洋偶極子指數(shù)。這兩項(xiàng)變量雖然與7月徑流過程的線性相關(guān)系數(shù)并不突出,但大量學(xué)者從物理成因角度對兩者與長江流域夏季洪澇災(zāi)害間的關(guān)系進(jìn)行了深入剖析。如印度洋偶極子,已有文獻(xiàn)認(rèn)為印度洋偶極子事件與弱厄爾尼諾現(xiàn)象共同促使了海洋羅斯貝波與溫躍層的下沉現(xiàn)象,引發(fā)了西北太平洋地區(qū)反氣旋異常,并加劇了東亞上空的西風(fēng)急流,從而導(dǎo)致長江中下游夏季暴雨頻發(fā)現(xiàn)象的產(chǎn)生[21];而南海副高活動則是通過增大西太平洋副熱帶高壓“北抬西升”效應(yīng),在冷暖空氣相互作用下,導(dǎo)致漢江上游易產(chǎn)生持續(xù)的強(qiáng)降雨,從而形成大洪水[22]。考慮文章篇幅問題,其他月份對應(yīng)預(yù)報因子篩選結(jié)果及對比并未展示在本文中。

    2.3模型結(jié)果分析

    本次研究將1956~2006年作為模型率定期,2007~2016年作為驗(yàn)證期,將2.2節(jié)中不同方案對應(yīng)預(yù)測因子導(dǎo)入隨機(jī)森林模型中,對丹江口水庫逐月徑流過程進(jìn)行模擬預(yù)測。為充分利用有限的資料并防止出現(xiàn)參數(shù)過擬合現(xiàn)象,研究采用留一交叉驗(yàn)證方法,將逐次模擬值與實(shí)測值間MAPE值最小作為優(yōu)化目標(biāo),確定了最優(yōu)參數(shù),最終不同方案對應(yīng)模的擬結(jié)果及相應(yīng)精度指標(biāo)如圖1~2所示及表2所列。

    由表2結(jié)果可知,基于兩種方案篩選出的預(yù)測因子,采用隨機(jī)森林模型對丹江口水庫月徑流過程進(jìn)行模擬,均取得較好的模擬效果。具體來看,率定期內(nèi)兩種方案Symbol值均大于80%,MAPE值小于30%,DC值大于0.8,而驗(yàn)證期的結(jié)果精度則略差于率定期,其中Symbol值較之率定期下降10%~14%,MAPE值變幅最大,整體提升20%~25%,DC值較之率定期下降0.20~0.21。對比2種方案,3種指標(biāo)均表明方案2對應(yīng)模擬精度顯著優(yōu)于方案1,這一現(xiàn)象說明方案2篩選的預(yù)報因子組合更加合理。

    圖1給出了丹江口水庫實(shí)測與模擬月流量過程對比圖,其結(jié)果表明在率定期及驗(yàn)證期內(nèi),豐水年份的模擬值略微偏小,而平水年及枯水年的模擬值略微偏大,對比圖1(a)與圖1(b),不難看出方案2對應(yīng)的模擬結(jié)果對汛期徑流的模擬性能要顯著優(yōu)于方案1。此外,通過對比率定期與驗(yàn)證期實(shí)測值與模擬值的散點(diǎn)圖(見圖2)可以看出,率定期及驗(yàn)證期的下半部分散點(diǎn)會略高于45°線,上半部分布散點(diǎn)則會略低于45°線,但所有散點(diǎn)都較為均勻分布在45°線附近,說明對應(yīng)精度尚可。對比2種方案模擬結(jié)果可知,方案2對應(yīng)結(jié)果整體更優(yōu),尤其對于汛期而言,其模擬結(jié)果與實(shí)測值更為接近。

    為進(jìn)一步比較不同方案模擬精度在年內(nèi)的變化特征,本次研究還計(jì)算了驗(yàn)證期內(nèi)各月份對應(yīng)模擬值與實(shí)測值間的相對誤差,如圖3所示。結(jié)果表明:除少數(shù)極端值外,方案1對應(yīng)模擬結(jié)果的相對誤差值絕大多數(shù)處于-65%~226%,而方案2對應(yīng)模擬結(jié)果的相對誤差則整體處于-57%~161%之間。同時,兩種方案對應(yīng)模擬結(jié)果精度的差異主要處于汛期內(nèi)(6~10月),這與上述結(jié)果得到的結(jié)論相互印證,說明方案2篩選出的預(yù)測因子組合對于預(yù)測來水量較大的汛期月份精度更高。

    3結(jié) 論

    本文提出了一種Copula熵與隨機(jī)森林模型相結(jié)合的中長期徑流預(yù)報方法,為了評估該方法的預(yù)測精度,研究將其應(yīng)用于漢江流域丹江口水庫逐月徑流預(yù)報中,并將相關(guān)系數(shù)篩選法作為比選方案進(jìn)行對比,結(jié)果表明:

    (1) 以不同方案篩選的預(yù)測因子作為輸入項(xiàng),隨機(jī)森林模型對于丹江口水庫逐月徑流的預(yù)測結(jié)果均取得了較高的精度,說明該模型能較好地應(yīng)用于丹江口水庫的中長期預(yù)報研究中。

    (2) 對比2種方案的模擬結(jié)果,不難看出基于Copula熵篩選的預(yù)報因子對應(yīng)預(yù)測精度顯著高于比選方案,尤其對于汛期而言,前者的預(yù)測值與實(shí)測值擬合優(yōu)度更好,說明采用Copula熵指標(biāo)篩選的預(yù)測因子組合具有更好的合理性。

    參考文獻(xiàn):

    [1]林劍藝,程春田.支持向量機(jī)在中長期徑流預(yù)報中的應(yīng)用[J].水利學(xué)報,2006,37(6):681-686.

    [2]楊旭,欒繼虹,馮國章.中長期水文預(yù)報研究評述與展望[J].西北農(nóng)業(yè)大學(xué)學(xué)報,2000,(6):203-207.

    [3]王文,馬駿.若干水文預(yù)報方法綜述[J].水利水電科技進(jìn)展,2005,25(1):56-60.

    [4]王富強(qiáng),霍風(fēng)霖.中長期水文預(yù)報方法研究綜述[J].人民黃河,2010,32(3):25-28.

    [5]張利平,王德智,夏軍,等.基于氣象因子的中長期水文預(yù)報方法研究[J].水電能源科學(xué),2003,21(3):4-6.

    [6]冉篤奎,李敏,武晟,等.丹江口水庫中長期徑流量的多模型預(yù)報結(jié)果分析及綜合研究[J].水利學(xué)報,2010,41(9):1069-1073.

    [7]李紅波,夏潮軍,王淑英.中長期徑流預(yù)報研究進(jìn)展及發(fā)展趨勢[J].人民黃河,2012,34(8):36-38,40.

    [8]陳柯兵,郭生練,王俊,等.長江上游ECMWF降水和徑流預(yù)報產(chǎn)品評估[J].人民長江,2020,51(3):73-80.

    [9]酈于杰,梁忠民,唐甜甜.基于支持向量回歸機(jī)的長期徑流預(yù)報及不確定性分析[J].南水北調(diào)與水利科技,2018,16(3):45-50.

    [10]許斌,楊鳳根,酈于杰.兩類集成學(xué)習(xí)算法在中長期徑流預(yù)報中的應(yīng)用[J].水力發(fā)電,2020,46(4):21-24,34.

    [11]謝帥,黃躍飛,李鐵鍵,等.LASSO回歸和支持向量回歸耦合的中長期徑流預(yù)報[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報,2018,26(4):709-722.

    [12]HUANG H,LIANG Z,LI B,et al.Combination of Multiple Data-Driven Models for long-term monthly runoff predictions based on Bayesian Model Averaging[J].Water Resources Management,2019,33(9):3321-3338.

    [13]鄭炎輝,張力瀾,田兆偉,等.基于隨機(jī)森林和RBF人工神經(jīng)網(wǎng)絡(luò)模型的新豐江水庫枯季入庫徑流中長期預(yù)報[J].人民珠江,2020,41(5):33-37.

    [14]王棟,魏加華,章四龍,等.基于CEEMD-BP模型的水文時間序列月徑流預(yù)測[J].北京師范大學(xué)學(xué)報(自然科學(xué)版),2020,56(3):376-386.

    [15]仕玉治,彭勇,周惠成.基于相關(guān)向量機(jī)的中長期徑流預(yù)報模型研究[J].大連理工大學(xué)學(xué)報,2012,52(1):79-84.

    [16]趙銅鐵鋼,楊大文.神經(jīng)網(wǎng)絡(luò)徑流預(yù)報模型中基于互信息的預(yù)報因子選擇方法[J].水力發(fā)電學(xué)報,2011,30(1):24-30.

    [17]MA J,SUN Z Q.Mutual information is Copula Entropy[J].Tsinghua Science and Technology,2011,16(1):51-54.

    [18]溫云亮,李艷玲,黃春艷,等.基于Copula熵理論的干旱驅(qū)動因子選擇[J].華北水利水電大學(xué)學(xué)報(自然科學(xué)版),2019,40(4):51-56.

    [19]陳璐,葉磊,盧韋偉,等.基于Copula熵的神經(jīng)網(wǎng)絡(luò)徑流預(yù)報模型預(yù)報因子選擇[J].水力發(fā)電學(xué)報,2014,33(6):25-29,90.

    [20]BREIMAN L.Random forests[J].Machine Learning,2001,45(1):5-32.

    [21]ZHOU Z Q,XIE S P,ZHANG R.Historic Yangtze flooding of 2020 tied to extreme Indian Ocean conditions[J].Proceedings of the National Academy of Sciences,2021,118(12):e2022255118.

    [22]陶玫,蔣薇,項(xiàng)瑛,等.1998和2010年長江流域汛期洪澇成因?qū)Ρ确治鯷J].氣象科學(xué),2012,32(3):282-287.

    (編輯:謝玲嫻)

    Abstract:As the key input of hydrological model for medium and long-term runoff forecast,forecast factors play an important role in improving the forecast accuracy.In order to further improve the accuracy of forecast results,we proposed a medium and long-term runoff forecast method combined with the Copula entropy and random forest model.For this method,the forecast factors were first selected based on the Copula entropy index,and the selected factors were used as input items and imported into the random forest model to simulate and forecast monthly runoff series.Finally,this method was applied to predict the monthly runoff series of the Danjiangkou Reservoir in Hanjiang River Basin,and compared with the correlation coefficient screening selection method.The results showed that the forecast results corresponding to the Copula entropy theory had a higher accuracy in forecasting monthly runoff series.Especially for flood season,the fitting effect of the simulated value and the measured value of this method was significantly better than that of the correlation coefficient method,indicating that the screened forecast factors are more reasonable.

    Key words:medium and long-term runoff forecast;forecast factors;global circulation factor;Copula entropy;random forest model;Danjiangkou Reservoir

    最近最新中文字幕大全免费视频| 亚洲国产欧洲综合997久久,| 欧美av亚洲av综合av国产av| 精品电影一区二区在线| 999精品在线视频| 麻豆成人av在线观看| 白带黄色成豆腐渣| 欧美日本视频| 岛国在线观看网站| 一二三四在线观看免费中文在| 中文亚洲av片在线观看爽| 嫩草影视91久久| 午夜a级毛片| 久久久久久国产a免费观看| 极品教师在线免费播放| cao死你这个sao货| 久久精品91无色码中文字幕| 亚洲成人久久爱视频| www.精华液| 丁香六月欧美| 亚洲午夜精品一区,二区,三区| 夜夜看夜夜爽夜夜摸| www.自偷自拍.com| 国产不卡一卡二| 亚洲av免费在线观看| 亚洲精品色激情综合| svipshipincom国产片| 窝窝影院91人妻| 在线十欧美十亚洲十日本专区| av在线蜜桃| 亚洲片人在线观看| 97超视频在线观看视频| 国产一区二区三区在线臀色熟女| 性色av乱码一区二区三区2| 99国产精品一区二区蜜桃av| 亚洲成av人片免费观看| av女优亚洲男人天堂 | 国产1区2区3区精品| 欧美乱码精品一区二区三区| 久久国产精品影院| 午夜福利高清视频| 国产不卡一卡二| 国产三级黄色录像| 欧美激情在线99| 国产亚洲av嫩草精品影院| av片东京热男人的天堂| 日本在线视频免费播放| 国产激情偷乱视频一区二区| 成人国产综合亚洲| 人妻夜夜爽99麻豆av| 一a级毛片在线观看| 女生性感内裤真人,穿戴方法视频| 亚洲aⅴ乱码一区二区在线播放| 禁无遮挡网站| 欧美黄色淫秽网站| 天堂网av新在线| 亚洲午夜理论影院| 三级男女做爰猛烈吃奶摸视频| 国产乱人视频| 久久婷婷人人爽人人干人人爱| 欧美一区二区国产精品久久精品| 国产黄色小视频在线观看| 别揉我奶头~嗯~啊~动态视频| 老司机福利观看| 国产1区2区3区精品| 精品不卡国产一区二区三区| 国产精品久久久久久亚洲av鲁大| 久久精品国产亚洲av香蕉五月| 波多野结衣高清作品| АⅤ资源中文在线天堂| 国产成人精品无人区| 99在线视频只有这里精品首页| 最近在线观看免费完整版| 国内毛片毛片毛片毛片毛片| 琪琪午夜伦伦电影理论片6080| 精华霜和精华液先用哪个| 成年女人看的毛片在线观看| 激情在线观看视频在线高清| 老司机在亚洲福利影院| 三级男女做爰猛烈吃奶摸视频| 免费无遮挡裸体视频| 成年女人看的毛片在线观看| 在线视频色国产色| 国产高清视频在线观看网站| 欧美在线一区亚洲| 女人被狂操c到高潮| 亚洲成av人片在线播放无| 12—13女人毛片做爰片一| 88av欧美| 美女cb高潮喷水在线观看 | 波多野结衣高清无吗| 亚洲国产欧美人成| 亚洲欧美日韩高清专用| 欧美大码av| 老司机在亚洲福利影院| www日本黄色视频网| 一个人免费在线观看的高清视频| 1024手机看黄色片| 国产又黄又爽又无遮挡在线| 成人高潮视频无遮挡免费网站| 亚洲专区中文字幕在线| 亚洲av免费在线观看| 亚洲av熟女| bbb黄色大片| 国产 一区 欧美 日韩| 婷婷亚洲欧美| 午夜福利18| 成人特级av手机在线观看| 嫁个100分男人电影在线观看| 成人三级做爰电影| 久久国产精品人妻蜜桃| 国产亚洲精品一区二区www| 国产亚洲av高清不卡| 2021天堂中文幕一二区在线观| 在线视频色国产色| 国产精品久久久人人做人人爽| 国产精品av视频在线免费观看| 亚洲激情在线av| 国产欧美日韩精品亚洲av| 国产精品av久久久久免费| 91久久精品国产一区二区成人 | 亚洲av日韩精品久久久久久密| 国产精品九九99| 18禁美女被吸乳视频| 国产精品久久视频播放| 无人区码免费观看不卡| 欧美另类亚洲清纯唯美| 中文字幕av在线有码专区| 亚洲电影在线观看av| 日本 av在线| 国产亚洲欧美在线一区二区| 国模一区二区三区四区视频 | 国产精品亚洲一级av第二区| 男人舔女人的私密视频| 久久香蕉国产精品| 观看美女的网站| 亚洲七黄色美女视频| 色播亚洲综合网| 男女视频在线观看网站免费| 搡老妇女老女人老熟妇| 亚洲欧美日韩高清专用| 日韩欧美国产在线观看| 搞女人的毛片| av中文乱码字幕在线| 18禁国产床啪视频网站| 观看美女的网站| 久久久久性生活片| 九色国产91popny在线| 日日夜夜操网爽| 免费看光身美女| 久久香蕉精品热| 久久国产乱子伦精品免费另类| 美女大奶头视频| 香蕉av资源在线| 亚洲无线在线观看| 色精品久久人妻99蜜桃| 国产成人av教育| 亚洲成av人片免费观看| 99久久综合精品五月天人人| 怎么达到女性高潮| 身体一侧抽搐| 欧美黑人欧美精品刺激| 性欧美人与动物交配| 在线观看66精品国产| 十八禁网站免费在线| 两人在一起打扑克的视频| 国语自产精品视频在线第100页| 欧美黑人欧美精品刺激| 成人亚洲精品av一区二区| 成年女人看的毛片在线观看| 亚洲欧美日韩卡通动漫| 亚洲av中文字字幕乱码综合| 他把我摸到了高潮在线观看| 久久精品国产99精品国产亚洲性色| 一二三四在线观看免费中文在| 久久久水蜜桃国产精品网| 国产精品亚洲av一区麻豆| 国产亚洲av高清不卡| 九九热线精品视视频播放| 国产成人啪精品午夜网站| 一进一出抽搐gif免费好疼| 一个人免费在线观看电影 | 成人三级黄色视频| 久9热在线精品视频| 好看av亚洲va欧美ⅴa在| 草草在线视频免费看| 欧美日韩一级在线毛片| 亚洲精品一卡2卡三卡4卡5卡| 宅男免费午夜| www国产在线视频色| 成年女人看的毛片在线观看| 亚洲国产精品成人综合色| 久久久精品大字幕| 99re在线观看精品视频| aaaaa片日本免费| 亚洲专区国产一区二区| 一夜夜www| 久久久精品大字幕| 黄色女人牲交| 国产日本99.免费观看| 曰老女人黄片| 免费大片18禁| 国产激情偷乱视频一区二区| 99热只有精品国产| 国产精品久久电影中文字幕| 日本五十路高清| 最近最新免费中文字幕在线| 一本综合久久免费| 日韩三级视频一区二区三区| 法律面前人人平等表现在哪些方面| 日韩欧美 国产精品| 亚洲av熟女| 欧美最黄视频在线播放免费| 丰满的人妻完整版| 中出人妻视频一区二区| 伊人久久大香线蕉亚洲五| 麻豆av在线久日| 91在线观看av| 亚洲国产欧美人成| 亚洲人成网站高清观看| 日韩欧美一区二区三区在线观看| 精品久久蜜臀av无| 日韩精品青青久久久久久| 久久天躁狠狠躁夜夜2o2o| 看黄色毛片网站| 久久久成人免费电影| 免费在线观看成人毛片| 老熟妇仑乱视频hdxx| 久久亚洲真实| 一本综合久久免费| 九九热线精品视视频播放| 美女扒开内裤让男人捅视频| 少妇人妻一区二区三区视频| 两性夫妻黄色片| 少妇熟女aⅴ在线视频| www.999成人在线观看| 久久九九热精品免费| 国产在线精品亚洲第一网站| 亚洲美女黄片视频| 一个人免费在线观看电影 | 久久久精品欧美日韩精品| 久久午夜综合久久蜜桃| 最新中文字幕久久久久 | 亚洲电影在线观看av| 狠狠狠狠99中文字幕| 可以在线观看毛片的网站| 久久精品夜夜夜夜夜久久蜜豆| 亚洲成人久久爱视频| 小说图片视频综合网站| 男女视频在线观看网站免费| 天天一区二区日本电影三级| 香蕉丝袜av| 午夜精品一区二区三区免费看| 男女下面进入的视频免费午夜| 国产成人aa在线观看| av天堂中文字幕网| 成人午夜高清在线视频| 亚洲自偷自拍图片 自拍| 亚洲精品粉嫩美女一区| 亚洲人成网站高清观看| av欧美777| x7x7x7水蜜桃| 亚洲专区中文字幕在线| av女优亚洲男人天堂 | 成人av在线播放网站| 久久国产精品影院| 国产免费男女视频| 色吧在线观看| 欧美极品一区二区三区四区| 精华霜和精华液先用哪个| bbb黄色大片| 久久精品91蜜桃| 日日干狠狠操夜夜爽| 国产1区2区3区精品| 久久中文看片网| 校园春色视频在线观看| 丁香欧美五月| 精品一区二区三区四区五区乱码| 日韩免费av在线播放| 日韩高清综合在线| 三级男女做爰猛烈吃奶摸视频| 亚洲av第一区精品v没综合| 一级毛片高清免费大全| 床上黄色一级片| 国产 一区 欧美 日韩| 最新中文字幕久久久久 | 又粗又爽又猛毛片免费看| a在线观看视频网站| 国产精品乱码一区二三区的特点| 一个人免费在线观看电影 | 夜夜爽天天搞| 日韩精品中文字幕看吧| 亚洲欧美日韩卡通动漫| 欧美日韩瑟瑟在线播放| 精品国产乱码久久久久久男人| 啦啦啦韩国在线观看视频| 香蕉av资源在线| 亚洲avbb在线观看| 午夜福利欧美成人| 成人三级做爰电影| 美女午夜性视频免费| 亚洲成a人片在线一区二区| 精品电影一区二区在线| 可以在线观看的亚洲视频| 国产极品精品免费视频能看的| 亚洲五月天丁香| 欧美绝顶高潮抽搐喷水| 一夜夜www| 国产熟女xx| 久久精品综合一区二区三区| 波多野结衣巨乳人妻| 99热只有精品国产| 91老司机精品| 亚洲,欧美精品.| 亚洲第一电影网av| 国产精品久久久久久精品电影| 欧美日韩亚洲国产一区二区在线观看| 啦啦啦免费观看视频1| 97超级碰碰碰精品色视频在线观看| 悠悠久久av| 两性夫妻黄色片| 九九在线视频观看精品| 丝袜人妻中文字幕| 熟女电影av网| 青草久久国产| 国产一级毛片七仙女欲春2| 一边摸一边抽搐一进一小说| 国产视频一区二区在线看| 日本 欧美在线| 白带黄色成豆腐渣| 午夜福利在线观看吧| 亚洲色图av天堂| 一夜夜www| 欧美乱色亚洲激情| 日韩欧美免费精品| 国产精品99久久久久久久久| 在线观看美女被高潮喷水网站 | 成年免费大片在线观看| 最好的美女福利视频网| e午夜精品久久久久久久| 又粗又爽又猛毛片免费看| 成人性生交大片免费视频hd| ponron亚洲| 好男人在线观看高清免费视频| 成年女人看的毛片在线观看| 久久久国产成人免费| 成年女人看的毛片在线观看| 久久久国产成人免费| 不卡一级毛片| 在线观看日韩欧美| 亚洲av成人不卡在线观看播放网| 亚洲成a人片在线一区二区| 亚洲精品在线观看二区| 国产成人一区二区三区免费视频网站| 国产av在哪里看| 99久久国产精品久久久| 精品人妻1区二区| 国产成人av教育| 又粗又爽又猛毛片免费看| 女警被强在线播放| 亚洲av成人不卡在线观看播放网| 两人在一起打扑克的视频| 女人高潮潮喷娇喘18禁视频| www国产在线视频色| 在线视频色国产色| 精品久久久久久成人av| 久久精品国产清高在天天线| 亚洲第一电影网av| 88av欧美| 国产1区2区3区精品| 琪琪午夜伦伦电影理论片6080| 亚洲av中文字字幕乱码综合| 欧美性猛交╳xxx乱大交人| 亚洲成av人片在线播放无| 最近视频中文字幕2019在线8| 人人妻人人看人人澡| 久久久精品欧美日韩精品| 国产亚洲精品av在线| 舔av片在线| 久久久久免费精品人妻一区二区| а√天堂www在线а√下载| 日本熟妇午夜| 亚洲国产欧美一区二区综合| 国产精品久久久久久久电影 | 久久久久久久久免费视频了| 亚洲精品粉嫩美女一区| 国产亚洲精品av在线| 免费大片18禁| 免费无遮挡裸体视频| 精品久久久久久久末码| 全区人妻精品视频| 最近在线观看免费完整版| 老熟妇仑乱视频hdxx| 热99re8久久精品国产| 此物有八面人人有两片| 免费看日本二区| 国产精品香港三级国产av潘金莲| 亚洲第一电影网av| 一二三四社区在线视频社区8| 成人一区二区视频在线观看| 99国产综合亚洲精品| 久久精品亚洲精品国产色婷小说| 欧美在线一区亚洲| 午夜两性在线视频| 在线观看日韩欧美| 久久人人精品亚洲av| 亚洲欧美精品综合久久99| 国产一区二区在线观看日韩 | 99久久精品一区二区三区| 琪琪午夜伦伦电影理论片6080| 9191精品国产免费久久| 搡老妇女老女人老熟妇| 在线观看66精品国产| 国产精华一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 久久这里只有精品19| 亚洲美女黄片视频| 国产久久久一区二区三区| 18禁国产床啪视频网站| 国产精品野战在线观看| 极品教师在线免费播放| 手机成人av网站| 麻豆av在线久日| 国产午夜福利久久久久久| 欧美日韩一级在线毛片| 精品久久久久久,| 亚洲av电影不卡..在线观看| 黄片小视频在线播放| 精品一区二区三区av网在线观看| 国产精品日韩av在线免费观看| 国产亚洲精品综合一区在线观看| 亚洲欧美日韩卡通动漫| 亚洲国产高清在线一区二区三| 亚洲色图 男人天堂 中文字幕| 成人国产一区最新在线观看| 99久久综合精品五月天人人| 香蕉av资源在线| 免费在线观看日本一区| 国产精品精品国产色婷婷| www日本在线高清视频| 久久天躁狠狠躁夜夜2o2o| 成年女人看的毛片在线观看| 精品无人区乱码1区二区| 国产精品久久久av美女十八| 一本精品99久久精品77| 日韩三级视频一区二区三区| 欧美乱色亚洲激情| 这个男人来自地球电影免费观看| 悠悠久久av| 夜夜夜夜夜久久久久| av欧美777| 男女那种视频在线观看| 国产一区二区三区在线臀色熟女| 手机成人av网站| 青草久久国产| 99国产极品粉嫩在线观看| 欧美一级毛片孕妇| 亚洲激情在线av| 国产欧美日韩精品一区二区| 一夜夜www| 黑人巨大精品欧美一区二区mp4| 99在线视频只有这里精品首页| av中文乱码字幕在线| 欧美另类亚洲清纯唯美| 一级毛片女人18水好多| or卡值多少钱| 婷婷亚洲欧美| 久9热在线精品视频| 亚洲真实伦在线观看| 亚洲欧美激情综合另类| 草草在线视频免费看| 怎么达到女性高潮| 91av网站免费观看| 在线观看66精品国产| 亚洲国产色片| 18禁黄网站禁片免费观看直播| 亚洲欧洲精品一区二区精品久久久| 亚洲精品美女久久久久99蜜臀| 亚洲av成人一区二区三| 亚洲精华国产精华精| 不卡一级毛片| 女人被狂操c到高潮| 色综合亚洲欧美另类图片| 成人欧美大片| 亚洲av熟女| 国产一级毛片七仙女欲春2| 日本撒尿小便嘘嘘汇集6| 亚洲片人在线观看| 成人欧美大片| 亚洲aⅴ乱码一区二区在线播放| 精品一区二区三区视频在线观看免费| 亚洲成a人片在线一区二区| 婷婷六月久久综合丁香| 丁香六月欧美| 色播亚洲综合网| 人妻夜夜爽99麻豆av| 成人鲁丝片一二三区免费| 欧美午夜高清在线| 波多野结衣巨乳人妻| 国产成人精品无人区| 身体一侧抽搐| 亚洲第一欧美日韩一区二区三区| 丰满人妻一区二区三区视频av | 国产精品一区二区精品视频观看| 久久久国产成人免费| 91av网一区二区| 国产 一区 欧美 日韩| 最近视频中文字幕2019在线8| 国产精品一区二区精品视频观看| 可以在线观看的亚洲视频| 在线观看一区二区三区| 99久久精品一区二区三区| 欧美zozozo另类| 欧美日韩黄片免| 又爽又黄无遮挡网站| 99热只有精品国产| 亚洲专区字幕在线| 桃红色精品国产亚洲av| av中文乱码字幕在线| 天堂影院成人在线观看| a级毛片a级免费在线| 亚洲色图av天堂| 可以在线观看的亚洲视频| av欧美777| 99久久精品一区二区三区| 91久久精品国产一区二区成人 | 日韩 欧美 亚洲 中文字幕| 人妻丰满熟妇av一区二区三区| 精品国产乱码久久久久久男人| 夜夜看夜夜爽夜夜摸| 日本黄色片子视频| 亚洲,欧美精品.| 国产伦一二天堂av在线观看| 久久九九热精品免费| 亚洲国产精品sss在线观看| 国产高清三级在线| 制服丝袜大香蕉在线| 夜夜看夜夜爽夜夜摸| 欧美黄色淫秽网站| 国产高清视频在线观看网站| 亚洲欧美日韩东京热| 真人一进一出gif抽搐免费| 精华霜和精华液先用哪个| 免费在线观看日本一区| 免费人成视频x8x8入口观看| 999久久久国产精品视频| www.精华液| 综合色av麻豆| 中文亚洲av片在线观看爽| 国产69精品久久久久777片 | 亚洲专区字幕在线| 精品欧美国产一区二区三| 99re在线观看精品视频| 国产1区2区3区精品| 色精品久久人妻99蜜桃| 人人妻人人澡欧美一区二区| 淫妇啪啪啪对白视频| 两个人视频免费观看高清| 国产精品综合久久久久久久免费| 亚洲男人的天堂狠狠| 亚洲中文av在线| 麻豆成人av在线观看| 美女高潮的动态| www.999成人在线观看| 久久久久国产精品人妻aⅴ院| 操出白浆在线播放| 丰满的人妻完整版| 国产精品99久久久久久久久| 麻豆成人午夜福利视频| 热99在线观看视频| 韩国av一区二区三区四区| 俺也久久电影网| 国产精品久久久av美女十八| 丁香六月欧美| 日韩高清综合在线| 97人妻精品一区二区三区麻豆| 国产激情欧美一区二区| 欧美成人一区二区免费高清观看 | 午夜久久久久精精品| 制服丝袜大香蕉在线| 日韩欧美 国产精品| 久久伊人香网站| 亚洲九九香蕉| 一区二区三区高清视频在线| 国内久久婷婷六月综合欲色啪| 视频区欧美日本亚洲| 全区人妻精品视频| 老司机福利观看| 久久久久久大精品| 亚洲人成网站在线播放欧美日韩| 亚洲第一欧美日韩一区二区三区| 亚洲中文av在线| 全区人妻精品视频| 老司机福利观看| 亚洲av电影在线进入| 国产精品久久久久久久电影 | h日本视频在线播放| e午夜精品久久久久久久| 国产欧美日韩精品亚洲av| 国产精品1区2区在线观看.| 午夜激情福利司机影院| 日本精品一区二区三区蜜桃| 精品福利观看| 国产97色在线日韩免费| 哪里可以看免费的av片| 在线观看66精品国产| 亚洲精品456在线播放app | 老熟妇乱子伦视频在线观看| 99久久精品国产亚洲精品| 国产成年人精品一区二区| 网址你懂的国产日韩在线| 色视频www国产| 国产一区二区在线av高清观看| 国产亚洲精品久久久久久毛片| 欧美一级毛片孕妇|