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

    降雨空間分布的模糊熵聚類分析

    2016-08-17 01:19:25張繼國管耀宗朱永忠
    關(guān)鍵詞:雨量站河海大學(xué)降雨量

    張繼國,管耀宗,朱永忠

    (1. 河海大學(xué)水利信息統(tǒng)計與管理研究所,江蘇 常州 213022; 2. 河海大學(xué)理學(xué)院,江蘇 南京 210098)

    ?

    降雨空間分布的模糊熵聚類分析

    張繼國1,管耀宗2,朱永忠2

    (1. 河海大學(xué)水利信息統(tǒng)計與管理研究所,江蘇 常州213022; 2. 河海大學(xué)理學(xué)院,江蘇 南京210098)

    為了提高降雨量插值精度,充分挖掘降雨變量信息,利用模糊熵聚類分析算法,對流域內(nèi)雨量站進(jìn)行模糊熵聚類研究,通過基于可能性分布和距離判定的聚類有效性函數(shù)確定模糊熵系數(shù)和聚類數(shù),從而得到模糊聚類結(jié)果,改進(jìn)原有的插值方法。以淮河流域蚌埠站以上區(qū)域99個雨量站雨量數(shù)據(jù),分別在一般情況下和模糊熵聚類情況下做交叉驗證,結(jié)果顯示,模糊熵聚類分析在反距離平方插值法中對降雨精度有所提升。

    降雨空間分布;降雨數(shù)據(jù)精度;模糊熵聚類分析;聚類有效性分析;降雨量插值

    HohaiUniversity,Changzhou213022,China;2.CollegeofSciences,HohaiUniversity,Nanjing210098,China)

    降雨量是水文模型中徑流模擬最基本、最主要的一個輸入項,是研究其他水文問題的基礎(chǔ)。其空間分布特征是影響產(chǎn)匯流模擬及其他一系列水文問題的重要控制因素[1]。隨著研究的深入,水文模型對降雨數(shù)據(jù)精度和廣度的要求越來越高。理論上獲取高精度降雨數(shù)據(jù)的方法是建立高密度的雨量站網(wǎng),但是由于經(jīng)濟條件和技術(shù)手段的約束,大部分地區(qū)氣象觀測站點數(shù)量不足,分布密度有限。因此,利用現(xiàn)有氣象觀測站的數(shù)據(jù),通過空間插值對觀測數(shù)據(jù)進(jìn)行補充尤為重要,孔云峰等[2]通過多種插值方法探究了美國德州的空間雨量數(shù)據(jù)。然而,大尺度流域上的降雨空間具有很強的時空分布不均勻性和復(fù)雜性,對此,在區(qū)域內(nèi)對已有站點作聚類分區(qū)處理,即將復(fù)雜的降雨測量站點系統(tǒng)劃分成不同的子系統(tǒng),減少不確定性因素的影響,是一種切實有效的研究方法[3]。李生辰等[4]在2007年研究了青藏高原降雨分區(qū)問題。楊絢等[5]在2008年通過降雨變化特征對中國干旱地區(qū)進(jìn)行了聚類劃分;鄭永宏等[6]2012年研究了湖北省的降雨分區(qū)問題。

    聚類分區(qū)主要分硬聚類和模糊聚類。相對于硬聚類,模糊聚類方法能夠?qū)︻惻c類之間有交叉的數(shù)據(jù)樣本集進(jìn)行有效的聚類,所得的聚類結(jié)果明顯優(yōu)于硬聚類方法。由于模糊聚類建立了數(shù)據(jù)樣本對于類別的不確定性的描述,表達(dá)了樣本類屬的模糊性,因此能夠更客觀地反映實際情況[7],并被廣泛應(yīng)用于水文研究中[8-9]。本文根據(jù)模糊熵聚類算法,將淮河流域蚌埠站以上的99個雨量站進(jìn)行模糊劃分,并研究模糊聚類分析在降雨量插值精度中的應(yīng)用,為流域內(nèi)雨量建模分析、水文循環(huán)研究、災(zāi)害預(yù)報等提供理論依據(jù)。

    1 模糊熵聚類

    1.1模糊熵目標(biāo)函數(shù)

    Tran等提出的模糊熵聚類算法(fuzzy entropy clustering)[10]是在模糊C均值聚類算法(fuzzyC-means clustering)基礎(chǔ)上,引入熵的概念,對隸屬度值分布進(jìn)行了算法優(yōu)化。

    (1)

    (2)

    式中:C——聚類數(shù);T——樣本數(shù);uit——成員xt對聚類Ci的隸屬度;n——模糊熵系數(shù);θi——聚類中心;d(xt,θi)——樣本xt與θi的差異距離。

    利用拉格朗日算子對目標(biāo)函數(shù)求極值可以得到模糊熵聚類算法隸屬度矩陣和聚類中心的更新方程(推導(dǎo)過程見文獻(xiàn)[10]):

    uit={∑Cj=1ed2(xt,θi)ed2(xt,θj)[?1n}-1θi=∑Tt=1uitxt∑Tt=1uitì?í?????ù?úúúúú(3)(4)

    1.2聚類數(shù)的確定

    相應(yīng)地,每一個聚類中心對T個樣本構(gòu)成的可能性分布也有一個香農(nóng)信息熵,由此定義可能性劃分熵為

    (5)

    范九倫[12]根據(jù)Bezdek劃分熵和可能性劃分熵定義了基于可能性分布的聚類有效性函數(shù)(式(6)),并提出當(dāng)HP(U;C)取得最大值的時候,U為最佳聚類隸屬度矩陣、C為最佳聚類數(shù)。

    (6)

    1.3模糊熵系數(shù)的確定

    在已知數(shù)據(jù)樣本中選取Xi、Xj、Xk、Xl(i、j、k、l=1,2,…,T)4個樣本作為樣本空間Ω,其隸屬度向量為ui、uj、uk、ul,定義基于距離判定的模糊聚類有效函數(shù)為

    (7)

    式中:P{*}——事件*發(fā)生占總樣本空間Ω的比例。

    1.4研究步驟

    1.5聚類合理性驗證

    通過真實降雨量矩陣Z、一般情況下的插值雨量矩陣Z′和模糊熵聚類情況下的插值雨量矩陣Z*計算多個交叉統(tǒng)計量[13]。記Xo,Xe分別是已知降雨量數(shù)據(jù)和交叉檢驗的插值計算降雨量數(shù)據(jù),N為數(shù)據(jù)個數(shù),各統(tǒng)計量(相關(guān)系數(shù)R、平均相對誤差RMAE、均方根誤差RMSE、復(fù)合相對誤差CRE)計算公式如下:

    (8)

    相關(guān)系數(shù)除去了偏差和方差的影響,考慮了插值估計數(shù)據(jù)與實際數(shù)據(jù)變化的同步性,表示了插值估計序列替代實際觀測序列的潛在能力。平均相對誤差和均方根誤差反映了插值估計序列與實際序列比較得到的誤差平均情況。復(fù)合相對誤差是描述插值序列與實際序列的相似性指標(biāo),該統(tǒng)計量對大誤差數(shù)據(jù)十分敏感。

    2 研 究 實 例

    2.1數(shù)據(jù)來源及相似度計算

    研究數(shù)據(jù)來自淮河流域蚌埠站以上區(qū)域99個雨量站1953—2013年732個月的降雨量數(shù)據(jù),站點基本情況見文獻(xiàn)[14]。

    淮河流域位于東經(jīng)112°~118°、北緯31°~35°的區(qū)域內(nèi),介于長江和黃河兩大流域之間。在氣候上,它處于南北氣候過渡帶,降水時空分布嚴(yán)重不均,差異較大?;春佑质俏覈媳狈降囊粭l自然分界線。因此,研究淮河流域的降水時空不確定性具有較高的科學(xué)價值。

    (9)

    (10)

    2.2聚類結(jié)果

    研究發(fā)現(xiàn),當(dāng)聚類數(shù)過大時,出現(xiàn)聚類中心彼此靠近并有重合現(xiàn)象;模糊熵系數(shù)小于0.02時聚類劃分過于分明,近似硬聚類,當(dāng)系數(shù)大于0.03時聚類間過于模糊。因此,擬定聚類數(shù)C=2、3、4、5,模糊熵系數(shù)n=0.02、0.022、0.024、0.026、0.028、0.03和停止閾值ε=0.000 1。計算相應(yīng)的聚類有效性值HP,當(dāng)C=2、3、4、5時,HP=0.000 4、0.101 1、0.058 4、0.032 6。比較可得,當(dāng)C=3時,HP取得最大值,因此選取最優(yōu)聚類數(shù)C*=3。

    在C*=3的情況下分別計算各模糊熵系數(shù)對應(yīng)隸屬度矩陣的距離判定有效性值E,當(dāng)n=0.02、0.022、0.024、0.026、0.028、0.03時,E=0.761 8、0.763 2、0.764 0、0.764 2、0.762 5、0.764 1。比較可得,在C*=3的情況下,當(dāng)n*=0.026時E取得最大值。從而得到相應(yīng)的最優(yōu)隸屬度矩陣U*。聚類結(jié)果如圖1所示(圖中所示星點表示已知雨量站位置)。

    圖1 聚類隸屬度分布Fig. 1 Distribution of degree of membership in EFC

    圖1中3幅圖分別表示3個聚類的隸屬度分布,其中顏色越淺表示該區(qū)域?qū)Υ司垲惖碾`屬度越高,反之越低。中部雨量站集中在第一聚類中,南部雨量站集中在第二聚類中,北部雨量站集中在第三聚類中??梢钥闯?算法得到的聚類結(jié)果有明顯的地理位置聚攏性,十分符合實際情況,即地理位置較近的地方,降雨量差異性較小。此外,分區(qū)多以緯度劃分為主,說明降雨量在緯度變化中的差異變化大,符合我國淮河流域基本自然情況,即南北降雨量差異大、東西降雨量差異小。

    表1 交叉驗證統(tǒng)計量對比Table 1 Comparison of statistics in cross-validation

    2.3驗證結(jié)果

    模糊熵聚類情況下和一般情況下反距離平方加權(quán)插值法各交叉統(tǒng)計量如表1所示。

    可以看出,與一般情況相比,模糊熵聚類分析情況下,站點雨量插值與測量值相關(guān)系數(shù)有所提高,各類誤差都有所下降,表明模糊熵聚類方法具有一定的優(yōu)越性。

    3 結(jié)  語

    通過模糊熵聚類分析,可以深入挖掘降雨信息在流域內(nèi)的分布,有利于更加深入地探究降雨系統(tǒng)內(nèi)部的關(guān)系,以及多種不確定因素。本文對淮河流域蚌埠站以上區(qū)域降雨量數(shù)據(jù)進(jìn)行了模糊熵聚類分析,獲得模糊熵聚類結(jié)果。同時,通過交叉驗證法,說明了模糊熵聚類算法在反距離平方加權(quán)插值中的實用性。

    [1] 石朋,芮孝芳. 降雨空間插值方法的比較與改進(jìn)[J]. 河海大學(xué)學(xué)報(自然科學(xué)版),2005,33(4):361-365. (SHI Peng, RUI Xiaofang. Comparison and improvement of spatial rainfall interpolation methods. [J]. Journal of Hehai University (Natural Sciences), 2005,33(4):361-365. (in Chinese))

    [2] 孔云峰,仝文偉. 降雨量地面觀測數(shù)據(jù)空間探索與插值方法探討[J]. 地理研究,2008,27(5):1097-1108.(KONG Yunfeng, TONG Wenwei. Spatial exploration and interpolation of the surface precipitation data[J]. Geographical Research, 2008, 27(5):1097-1108. (in Chinese))

    [3] 張繼國, 謝平, 龔艷冰, 等.降雨信息空間插值研究評述與展望[J].水資源與水工程學(xué)報, 2012, 23(1):6-9. (ZHANG Jiguo, XIE Ping, GONG Yanbing, et al. Review and perspectives of the research on spatial interpolation of rainfall information [J]. Journal of Water Resources and Water Engineering, 2012,23(1):6-9.(in Chinese))

    [4] 李生辰, 徐亮, 郭英香, 等.近34a青藏高原年降水變化及其分區(qū)[J].中國沙漠, 2007, 27(2):307-314. (LI Shengchen, XU Liang, GUO Yingxiang, et al. Change of annual precipitation over QinghaiXizang Plateau and subregions in recent 34 years [J]. Journal of Desert Research, 2007, 27(2): 307-314. (in Chinese))

    [5] 楊絢, 李棟梁.中國干旱氣候分區(qū)及其降水量變化特征[J].干旱氣象, 2008, 26(2):17-24. (YANG Xuan, LI Dongliang. Precipitation variation characteristics and arid climate division in China [J]. Arid Meteorology, 2008, 26(2):17-24.(in Chinese))

    [6] 鄭永宏, 林愛文, 代偵勇.湖北省降水分區(qū)研究[J].長江流域資源與環(huán)境, 2012, 21(7):859-863. (ZHENG Yonghong, LIN Aiwen, DAI Zhenyong. Research on precipitation regionalization in Hubei Province [J]. Resources and Environment in the Yangtze Basin, 2012,21(7):859-863.(in Chinese))

    [7] 雷鳴. 模糊聚類新算法的研究[D].天津:天津大學(xué),2007.

    [8] 冀鴻蘭,卞雪軍,徐晶. 黃河內(nèi)蒙古段流凌預(yù)報可變模糊聚類循環(huán)迭代模型[J]. 水利水電科技進(jìn)展,2013,33(4):14-17. (JI Honglan, BIAN Xuejun, XU Jing. Variable fuzzy clustering loop iteration model for ice-run forecast in Inner Mongolia reach of Yellow River [J]. Advances in Science and Technology of Water Resources, 2013, 33(4):14-17.(in Chinese))

    [9] 樊哲超,陳建生,董海洲,等. 應(yīng)用環(huán)境同位素和模糊聚類方法研究堤防滲漏[J]. 水利水電科技進(jìn)展,2005,25(2):8-10,57. (FAN Zhechao, CHEN Jiansheng, DONG Haizhou, et al. Application of environmental isotope and fuzzy clustering method to study of seepage from dykes [J]. Advances In Science and Technology of Water Resources, 2005, 25(2):8-10,57.(in Chinese))

    [10] WU Xiaohong, ZHOU Jianjiang. Possibilistic fuzzy entropy clustering[J]. Journal of Computational Information Systems, 2007, 3(1):25-33.

    [11] BEZDEK J C. Pattem recognition with fuzzy objective function algorithms [M]. New York: Plenum, 1981.

    [12] 范九倫. 模糊聚類新算法與聚類有效性問題研究[D].西安:西安電子科技大學(xué),1998.

    [13] 熊秋芬,黃玫,熊敏詮,等. 基于國家氣象觀測站逐日降水格點數(shù)據(jù)的交叉檢驗誤差分析[J]. 高原氣象,2011,30(6):1615-1625.(XIONG Qiufen, HUANG Mei, XIONG Minquan, et al. Cross-validation error analysis of daily gridded precipitation based on China meteorological observation [J]. Plateau Meteorology,2011,30(6):1615-1625. (in Chinese))

    [14] 張繼國.降雨時空分布不均勻性信息熵研究[D].南京:河海大學(xué),2004.

    Fuzzy entropy clustering analysis of spatial distribution of precipitation

    ZHANG Jiguo1, GUAN Yaozong2, ZHU Yongzhong2

    (1.InstituteofInformationStatisticsandManagementofWaterResources,

    In order to improve the accuracy of precipitation interpolation and fully explore the information regarding precipitation variables, fuzzy entropy clustering (FEC) was carried out at rain gauge stations in a basin. A clustering validity function, based on possibility distribution and distance determination, was used to determine the fuzzy entropy coefficient and the number of clusters, so as to obtain the fuzzy clustering results and improve the original interpolation method. Based on data from 99 rain gauge stations located above the Bengbu Station in the Huaihe River Basin, cross validation was conducted under non-clustering and FEC conditions. The results demonstrate that FEC improves the precipitation accuracy in the inverse distance squared interpolation method.

    spatial distribution of precipitation; accuracy of precipitation data; fuzzy entropy clustering (FEC) analysis; clustering validity analysis; precipitation interpolation

    1000-1980(2016)04-0353-05

    10.3876/j.issn.1000-1980.2016.04.012

    2015-10-17

    江蘇省自然科學(xué)基金(BK20131135);江蘇省自然科學(xué)青年基金(BK20130242)

    張繼國(1956—),男,湖北漢川人,教授,主要從事水文不確定性分析、信息熵理論與方法研究。E-mail:zhangjg@hhuc.edu.cn

    P467

    A

    猜你喜歡
    雨量站河海大學(xué)降雨量
    遼西山洪災(zāi)害重點防治區(qū)遙測雨量站點優(yōu)化布設(shè)研究
    《河海大學(xué)學(xué)報(哲學(xué)社會科學(xué)版)》征稿簡則
    《河海大學(xué)學(xué)報(哲學(xué)社會科學(xué)版)》征稿簡則
    降雨量與面積的關(guān)系
    信息熵方法在遼寧省不同分區(qū)雨量站網(wǎng)布設(shè)的應(yīng)用研究
    河海大學(xué)學(xué)術(shù)期刊創(chuàng)辦百年
    中小流域雨量站網(wǎng)密度規(guī)劃與研究
    洞庭湖區(qū)降雨特性分析
    雨量站網(wǎng)測量精度的評估
    河海大學(xué)迎來建校百年華誕
    中國水利(2015年20期)2015-02-01 07:50:47
    一卡2卡三卡四卡精品乱码亚洲| 久久久久免费精品人妻一区二区| 51国产日韩欧美| 十八禁网站免费在线| 女同久久另类99精品国产91| 精品99又大又爽又粗少妇毛片 | 亚洲av第一区精品v没综合| 亚洲国产精品sss在线观看| 精品不卡国产一区二区三区| 欧美+日韩+精品| 啦啦啦韩国在线观看视频| 在线免费十八禁| 成年女人永久免费观看视频| 亚洲美女黄片视频| 91午夜精品亚洲一区二区三区 | 波多野结衣高清无吗| 日韩精品中文字幕看吧| 亚洲在线观看片| 精品午夜福利视频在线观看一区| 欧美色视频一区免费| 嫩草影院入口| 九九热线精品视视频播放| 国产精品亚洲一级av第二区| 韩国av一区二区三区四区| 国产精品野战在线观看| 精品久久久噜噜| 久久久久久久久久久丰满 | 国产高清三级在线| 午夜久久久久精精品| 亚洲第一电影网av| 欧美成人免费av一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 在线播放国产精品三级| 精品99又大又爽又粗少妇毛片 | 亚洲中文字幕日韩| 国产高清不卡午夜福利| 老司机深夜福利视频在线观看| 免费人成在线观看视频色| 天堂网av新在线| 国产综合懂色| 如何舔出高潮| 又爽又黄无遮挡网站| 干丝袜人妻中文字幕| 成年免费大片在线观看| 免费看日本二区| 99久国产av精品| 亚洲不卡免费看| 色尼玛亚洲综合影院| 在线播放国产精品三级| 中文字幕人妻熟人妻熟丝袜美| 女的被弄到高潮叫床怎么办 | 国产一区二区三区av在线 | 午夜精品久久久久久毛片777| 国产精品嫩草影院av在线观看 | 成熟少妇高潮喷水视频| 欧美成人免费av一区二区三区| 美女 人体艺术 gogo| 成人性生交大片免费视频hd| 久99久视频精品免费| 久久久色成人| 别揉我奶头~嗯~啊~动态视频| 亚洲最大成人手机在线| 亚洲欧美日韩高清专用| 日本一二三区视频观看| 禁无遮挡网站| av在线蜜桃| 在线观看66精品国产| 变态另类成人亚洲欧美熟女| 成人高潮视频无遮挡免费网站| 久久久久国产精品人妻aⅴ院| 久久久久久久久久成人| 色综合色国产| 国内精品一区二区在线观看| 日韩一区二区视频免费看| 亚洲三级黄色毛片| 久久久精品大字幕| av在线老鸭窝| 国产综合懂色| 99热这里只有是精品50| 久久99热6这里只有精品| 两个人视频免费观看高清| 久久精品国产亚洲av香蕉五月| 婷婷色综合大香蕉| 91av网一区二区| 国产午夜精品久久久久久一区二区三区 | 成人国产一区最新在线观看| 国产精品一区二区性色av| 中文亚洲av片在线观看爽| 欧美日韩黄片免| 不卡一级毛片| 国产精品国产高清国产av| 亚洲精品一卡2卡三卡4卡5卡| 中文字幕av在线有码专区| 老女人水多毛片| 精品久久久久久成人av| 精品人妻一区二区三区麻豆 | 噜噜噜噜噜久久久久久91| 免费看美女性在线毛片视频| 国产一区二区三区在线臀色熟女| 狂野欧美激情性xxxx在线观看| 听说在线观看完整版免费高清| 国产成人福利小说| 久久久久久大精品| 亚洲最大成人av| 久久久久久大精品| 国产精品自产拍在线观看55亚洲| 最近中文字幕高清免费大全6 | 在现免费观看毛片| 一级黄片播放器| 99久国产av精品| 精品久久久久久久久久久久久| 级片在线观看| 少妇的逼水好多| 18+在线观看网站| 亚洲午夜理论影院| 日韩精品中文字幕看吧| 在线a可以看的网站| 老熟妇仑乱视频hdxx| 欧美又色又爽又黄视频| eeuss影院久久| 丰满的人妻完整版| 综合色av麻豆| 99九九线精品视频在线观看视频| 性插视频无遮挡在线免费观看| 午夜a级毛片| 男女啪啪激烈高潮av片| 国产乱人视频| 日本与韩国留学比较| 国产午夜精品久久久久久一区二区三区 | 国产成人影院久久av| 在线观看午夜福利视频| 亚洲欧美日韩高清专用| 国产精品伦人一区二区| 亚洲第一电影网av| 国内精品久久久久精免费| 看片在线看免费视频| 日韩欧美在线乱码| 欧美日韩国产亚洲二区| av福利片在线观看| 精品一区二区三区视频在线观看免费| 一区二区三区高清视频在线| ponron亚洲| 欧美又色又爽又黄视频| 少妇被粗大猛烈的视频| 国产蜜桃级精品一区二区三区| 99在线人妻在线中文字幕| 69人妻影院| 久久久久九九精品影院| 最近在线观看免费完整版| 非洲黑人性xxxx精品又粗又长| 亚洲一区高清亚洲精品| 亚洲第一电影网av| 人人妻人人看人人澡| 两性午夜刺激爽爽歪歪视频在线观看| 欧美人与善性xxx| 在线a可以看的网站| 日韩精品青青久久久久久| 又粗又爽又猛毛片免费看| 一级a爱片免费观看的视频| 五月伊人婷婷丁香| 久久婷婷人人爽人人干人人爱| 亚洲熟妇熟女久久| 欧美日韩乱码在线| 1024手机看黄色片| 日本黄色片子视频| 99久久成人亚洲精品观看| 欧美高清成人免费视频www| 国产熟女欧美一区二区| 日韩欧美国产一区二区入口| 午夜a级毛片| 麻豆国产97在线/欧美| 在线观看美女被高潮喷水网站| 美女cb高潮喷水在线观看| 亚洲国产色片| 成人亚洲精品av一区二区| 天堂av国产一区二区熟女人妻| 中国美白少妇内射xxxbb| 全区人妻精品视频| 波多野结衣高清无吗| 色吧在线观看| 国产精品自产拍在线观看55亚洲| 欧美成人a在线观看| 欧美色欧美亚洲另类二区| 欧美日韩亚洲国产一区二区在线观看| 欧美高清性xxxxhd video| 国产v大片淫在线免费观看| 国产又黄又爽又无遮挡在线| 小蜜桃在线观看免费完整版高清| 亚洲精品一卡2卡三卡4卡5卡| 午夜免费成人在线视频| 我的女老师完整版在线观看| 露出奶头的视频| 岛国在线免费视频观看| 乱系列少妇在线播放| 99热这里只有精品一区| 人人妻人人看人人澡| 国产视频一区二区在线看| 人人妻,人人澡人人爽秒播| 91麻豆精品激情在线观看国产| 超碰av人人做人人爽久久| 99久久精品国产国产毛片| 少妇人妻精品综合一区二区 | 亚洲人成网站在线播放欧美日韩| 久久精品人妻少妇| 国产高清视频在线观看网站| 最近最新免费中文字幕在线| 亚洲欧美激情综合另类| 免费观看精品视频网站| av福利片在线观看| 午夜影院日韩av| 男女边吃奶边做爰视频| 一区二区三区免费毛片| 女人被狂操c到高潮| 精品国内亚洲2022精品成人| 白带黄色成豆腐渣| 大型黄色视频在线免费观看| 香蕉av资源在线| 色播亚洲综合网| 精品午夜福利视频在线观看一区| 亚洲精华国产精华液的使用体验 | 99久久无色码亚洲精品果冻| 狂野欧美白嫩少妇大欣赏| 老师上课跳d突然被开到最大视频| 伊人久久精品亚洲午夜| 亚洲中文字幕日韩| 精华霜和精华液先用哪个| 九九热线精品视视频播放| 久久国内精品自在自线图片| 精品日产1卡2卡| 欧美日本视频| 在线国产一区二区在线| 日韩大尺度精品在线看网址| 三级男女做爰猛烈吃奶摸视频| 自拍偷自拍亚洲精品老妇| 精品午夜福利在线看| 禁无遮挡网站| 色尼玛亚洲综合影院| 久久精品国产亚洲av涩爱 | 日韩高清综合在线| 亚洲图色成人| 白带黄色成豆腐渣| 国产三级在线视频| 色哟哟·www| 精品欧美国产一区二区三| 中文字幕人妻熟人妻熟丝袜美| 久久久久国产精品人妻aⅴ院| 精品人妻1区二区| 天堂√8在线中文| 毛片女人毛片| 高清在线国产一区| 又黄又爽又免费观看的视频| 亚洲自拍偷在线| 变态另类丝袜制服| 日韩亚洲欧美综合| 久久久久久久亚洲中文字幕| 全区人妻精品视频| 午夜福利欧美成人| 美女xxoo啪啪120秒动态图| 日本撒尿小便嘘嘘汇集6| 午夜日韩欧美国产| 日韩精品青青久久久久久| 国内毛片毛片毛片毛片毛片| 特大巨黑吊av在线直播| 熟妇人妻久久中文字幕3abv| 亚洲自拍偷在线| 国产黄色小视频在线观看| 麻豆成人午夜福利视频| 老师上课跳d突然被开到最大视频| 日韩,欧美,国产一区二区三区 | 国内少妇人妻偷人精品xxx网站| 十八禁网站免费在线| 免费观看的影片在线观看| 国产精品永久免费网站| 日韩,欧美,国产一区二区三区 | 亚洲 国产 在线| 国产91精品成人一区二区三区| videossex国产| 久久精品国产亚洲av天美| 白带黄色成豆腐渣| 亚洲精品成人久久久久久| 国产精品三级大全| 国产成人影院久久av| 又爽又黄a免费视频| 一级黄色大片毛片| 日本撒尿小便嘘嘘汇集6| 亚洲不卡免费看| 少妇高潮的动态图| 一区福利在线观看| 中亚洲国语对白在线视频| 国产av一区在线观看免费| 网址你懂的国产日韩在线| 久久精品久久久久久噜噜老黄 | av女优亚洲男人天堂| 成人欧美大片| 91av网一区二区| 久久久色成人| 日本在线视频免费播放| 免费无遮挡裸体视频| 色综合亚洲欧美另类图片| 天堂√8在线中文| 日韩欧美在线二视频| 老熟妇乱子伦视频在线观看| 婷婷精品国产亚洲av在线| 亚洲最大成人av| 偷拍熟女少妇极品色| 日韩av在线大香蕉| 极品教师在线视频| 国语自产精品视频在线第100页| 午夜a级毛片| 桃红色精品国产亚洲av| 久久久久久久亚洲中文字幕| 国产精品亚洲一级av第二区| 深夜精品福利| av女优亚洲男人天堂| 国产精品电影一区二区三区| 国产精品久久久久久久电影| www.色视频.com| 日韩欧美精品免费久久| 国产高清三级在线| 国产极品精品免费视频能看的| 精品午夜福利在线看| 狂野欧美白嫩少妇大欣赏| 日日摸夜夜添夜夜添av毛片 | 男插女下体视频免费在线播放| 日韩大尺度精品在线看网址| 不卡一级毛片| 久久久久久久午夜电影| 欧美色欧美亚洲另类二区| 亚洲国产精品合色在线| 久久久久久久久大av| 久久天躁狠狠躁夜夜2o2o| 欧美另类亚洲清纯唯美| 国语自产精品视频在线第100页| 国产91精品成人一区二区三区| 一本一本综合久久| 国产高清三级在线| 国产午夜福利久久久久久| 成年女人永久免费观看视频| 一级黄片播放器| 两性午夜刺激爽爽歪歪视频在线观看| 自拍偷自拍亚洲精品老妇| 国产男人的电影天堂91| 亚洲精品成人久久久久久| 丝袜美腿在线中文| 国产又黄又爽又无遮挡在线| 亚洲成人久久性| 久久久久久久午夜电影| 国产单亲对白刺激| 日韩精品中文字幕看吧| 九九久久精品国产亚洲av麻豆| 日韩高清综合在线| 成人美女网站在线观看视频| 亚洲精品一区av在线观看| 免费看美女性在线毛片视频| 亚洲久久久久久中文字幕| av福利片在线观看| 少妇的逼好多水| 美女高潮的动态| 一级黄色大片毛片| 美女高潮喷水抽搐中文字幕| 三级毛片av免费| 亚洲性久久影院| 久久久成人免费电影| 一边摸一边抽搐一进一小说| 人妻丰满熟妇av一区二区三区| 悠悠久久av| 乱人视频在线观看| 精品午夜福利视频在线观看一区| 黄色日韩在线| 97超视频在线观看视频| 99在线视频只有这里精品首页| 麻豆一二三区av精品| 亚洲美女黄片视频| 一个人观看的视频www高清免费观看| 中文字幕人妻熟人妻熟丝袜美| 国产爱豆传媒在线观看| 琪琪午夜伦伦电影理论片6080| 成人二区视频| 黄色日韩在线| 变态另类丝袜制服| 蜜桃亚洲精品一区二区三区| 麻豆成人av在线观看| 亚洲无线在线观看| 91麻豆av在线| 99视频精品全部免费 在线| 色吧在线观看| x7x7x7水蜜桃| 色5月婷婷丁香| 99久久中文字幕三级久久日本| 好男人在线观看高清免费视频| 波多野结衣巨乳人妻| 成人国产麻豆网| 国产人妻一区二区三区在| 成人高潮视频无遮挡免费网站| 日日撸夜夜添| 国内精品宾馆在线| 亚洲无线观看免费| 亚洲综合色惰| 很黄的视频免费| h日本视频在线播放| 男女边吃奶边做爰视频| 成人亚洲精品av一区二区| 欧美日韩黄片免| 免费看美女性在线毛片视频| 最近最新中文字幕大全电影3| 村上凉子中文字幕在线| 99riav亚洲国产免费| 在线国产一区二区在线| 99久久中文字幕三级久久日本| 欧美激情在线99| 又黄又爽又刺激的免费视频.| 大型黄色视频在线免费观看| 在线观看美女被高潮喷水网站| 高清在线国产一区| 亚洲精品亚洲一区二区| 久久天躁狠狠躁夜夜2o2o| 日韩一区二区视频免费看| 在线观看一区二区三区| 别揉我奶头~嗯~啊~动态视频| 亚洲va日本ⅴa欧美va伊人久久| 国产欧美日韩精品亚洲av| 一级av片app| 亚洲成人免费电影在线观看| 18禁在线播放成人免费| 国产大屁股一区二区在线视频| 男人舔奶头视频| 国模一区二区三区四区视频| 人人妻人人看人人澡| 国产91精品成人一区二区三区| 久久草成人影院| 国产女主播在线喷水免费视频网站 | 麻豆成人午夜福利视频| 91狼人影院| 久久久久久伊人网av| 精品无人区乱码1区二区| 欧美区成人在线视频| 免费大片18禁| 免费观看精品视频网站| 少妇被粗大猛烈的视频| 欧美xxxx性猛交bbbb| av福利片在线观看| 黄色视频,在线免费观看| 久久国产精品人妻蜜桃| 自拍偷自拍亚洲精品老妇| 亚洲av免费高清在线观看| 国产毛片a区久久久久| 精品久久国产蜜桃| 国产综合懂色| 亚洲成人精品中文字幕电影| 简卡轻食公司| 国产高清不卡午夜福利| 超碰av人人做人人爽久久| 国内精品久久久久精免费| 久久亚洲精品不卡| 成年免费大片在线观看| 成人特级黄色片久久久久久久| 男女之事视频高清在线观看| 99精品在免费线老司机午夜| www.色视频.com| 女人被狂操c到高潮| 国产私拍福利视频在线观看| 精品欧美国产一区二区三| 国产精品乱码一区二三区的特点| 久久久久久久久大av| 老司机午夜福利在线观看视频| 午夜福利高清视频| 18禁裸乳无遮挡免费网站照片| 少妇高潮的动态图| 色综合站精品国产| 啦啦啦观看免费观看视频高清| 亚洲av美国av| 别揉我奶头 嗯啊视频| 婷婷色综合大香蕉| 亚洲成a人片在线一区二区| 日韩一本色道免费dvd| 在线免费十八禁| 亚洲中文字幕日韩| 嫩草影院入口| 欧美不卡视频在线免费观看| 内射极品少妇av片p| 国产一区二区在线av高清观看| 久久久久精品国产欧美久久久| 男女啪啪激烈高潮av片| 久久亚洲真实| 亚洲最大成人中文| 国产精品爽爽va在线观看网站| 成人欧美大片| 国产男人的电影天堂91| 在线播放国产精品三级| 国内揄拍国产精品人妻在线| 成人鲁丝片一二三区免费| 亚洲av中文av极速乱 | 欧美区成人在线视频| 日本 欧美在线| 黄色一级大片看看| 舔av片在线| 精品国内亚洲2022精品成人| 亚洲av免费在线观看| 亚洲第一电影网av| 国产主播在线观看一区二区| 色av中文字幕| 美女大奶头视频| 最新在线观看一区二区三区| 亚洲av中文av极速乱 | 国产免费一级a男人的天堂| 成年人黄色毛片网站| 网址你懂的国产日韩在线| 成人午夜高清在线视频| 日韩国内少妇激情av| 成人特级av手机在线观看| 国产激情偷乱视频一区二区| 国产成人福利小说| 可以在线观看毛片的网站| 国产精品国产三级国产av玫瑰| 99热这里只有精品一区| 18禁黄网站禁片免费观看直播| 蜜桃亚洲精品一区二区三区| 色在线成人网| 国产日本99.免费观看| 精品人妻偷拍中文字幕| 一本一本综合久久| 亚洲一区二区三区色噜噜| 成年人黄色毛片网站| 日本-黄色视频高清免费观看| 国产精华一区二区三区| av.在线天堂| 毛片一级片免费看久久久久 | 久久精品国产亚洲av香蕉五月| 中出人妻视频一区二区| 美女大奶头视频| 精品日产1卡2卡| 一级a爱片免费观看的视频| 久久精品国产99精品国产亚洲性色| 国产三级在线视频| 狂野欧美白嫩少妇大欣赏| 97碰自拍视频| 97超级碰碰碰精品色视频在线观看| 一区二区三区高清视频在线| 久9热在线精品视频| 国产极品精品免费视频能看的| 一本久久中文字幕| 国产精品久久久久久久电影| 最后的刺客免费高清国语| 免费av不卡在线播放| 欧美xxxx黑人xx丫x性爽| 99视频精品全部免费 在线| 国产精品久久视频播放| 欧美最新免费一区二区三区| 一区二区三区免费毛片| 国语自产精品视频在线第100页| 久久亚洲真实| 在线a可以看的网站| 男插女下体视频免费在线播放| 亚洲国产高清在线一区二区三| 午夜福利在线观看吧| 久久香蕉精品热| 国产亚洲精品久久久com| 日本-黄色视频高清免费观看| 国产亚洲精品久久久com| 亚洲第一区二区三区不卡| 亚洲欧美日韩高清专用| 国产女主播在线喷水免费视频网站 | 中国美白少妇内射xxxbb| 深爱激情五月婷婷| 色5月婷婷丁香| 一个人免费在线观看电影| 婷婷亚洲欧美| 日韩高清综合在线| 身体一侧抽搐| 国产麻豆成人av免费视频| 欧美又色又爽又黄视频| 亚洲人成伊人成综合网2020| 此物有八面人人有两片| 两人在一起打扑克的视频| 搡老岳熟女国产| 亚洲国产精品久久男人天堂| 草草在线视频免费看| 一本精品99久久精品77| 波多野结衣高清作品| 免费在线观看日本一区| a级毛片a级免费在线| 亚洲无线观看免费| 91在线精品国自产拍蜜月| 又黄又爽又免费观看的视频| 99国产精品一区二区蜜桃av| 国语自产精品视频在线第100页| 亚洲精华国产精华精| 九九在线视频观看精品| 禁无遮挡网站| 国产黄片美女视频| 亚洲黑人精品在线| 三级男女做爰猛烈吃奶摸视频| 最后的刺客免费高清国语| 永久网站在线| 中文字幕精品亚洲无线码一区| 亚洲美女黄片视频| 成人亚洲精品av一区二区| 久久99热这里只有精品18| 国产在线精品亚洲第一网站| 亚洲av美国av| 久久久精品大字幕| 日韩强制内射视频| 日本熟妇午夜| 午夜老司机福利剧场| av在线蜜桃| aaaaa片日本免费| 成人三级黄色视频| 欧美高清性xxxxhd video| 国产精品国产三级国产av玫瑰| 联通29元200g的流量卡| 黄色视频,在线免费观看| 日韩国内少妇激情av| 欧美xxxx黑人xx丫x性爽| 永久网站在线| 亚洲欧美激情综合另类|