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

    改進POT模型及其在邊坡安全監(jiān)測預警中的應用

    2016-07-09 14:37:07周志杰沈振中徐力群任杰
    南水北調(diào)與水利科技 2016年4期
    關(guān)鍵詞:閾值邊坡

    周志杰 沈振中 徐力群 任杰

    摘要:基于極值理論的BMM(Block Maximum Method)和POT模型是近來分析邊坡安全監(jiān)測資料、評估邊坡安全狀況的新興方法之一。相對簡便的BMM模型在數(shù)據(jù)取樣時往往忽略區(qū)間次極大值,在資料年限較短時樣本容量偏小,可能導致所得結(jié)果誤差較大。本文利用改進的Hill估計方法得到閾值,通過極大似然估計確定廣義帕累托分布參數(shù),從而利用超限數(shù)據(jù)序列來確定測值序列的整體分布,提出了改進POT(Modified Peaks over Threshold)模型,并應用于某邊坡工程的安全監(jiān)測預警指標分析。結(jié)果表明,在同一置信水平下利用超限值應用廣義帕累托分布擬合得到的預警指標小于利用塊極大值應用正態(tài)分布得到的預警指標,表明基于超限數(shù)據(jù)的改進POT模型得到的預警指標更能有效規(guī)避極端情況發(fā)生的風險,更有利于邊坡安全監(jiān)測和預警。

    關(guān)鍵詞:邊坡;極值理論;BMM模型;改進POT模型;預警指標;閾值;極大似然估計

    中圖分類號:U416 文獻標志碼:A 文章編號:1672-1683(2016)04-0192-06

    Abstract:BMM model and POT model which based on extreme value theory are one of the newly-developing methods to assess the safety conditions of slope in recent years.The ignorant of secondary maximum points and the relatively small size of sample in the case of short-time data may lead to a big error in BMM model.A modified POT model was proposed in the paper to get the threshold with improved Hill estimator method and obtain the Pareto parameters with the maximum likelihood estimation so that the whole distribution could be determined.With the analysis of a slope engineering,the early-warning index which used the generalized Pareto distribution with the data exceeding the threshold was less than that using the normal distribution with the block maximum values in the case of same confidence level.It could be drawn from the results that the early-warning index based on modified POT model was more effective in slope safety monitoring,which could be made use of reducing the risk when extreme conditions happened.

    Key words: slope;extreme value theory;BMM model;modified POT model;early-warning index;threshold;maximum likelihood estimation

    依據(jù)邊坡工程施工、運行的原型監(jiān)測資料,應用統(tǒng)計數(shù)學、力學等方法建立監(jiān)控模型,擬定不同置信水平下的預警指標,展開對邊坡的實時監(jiān)控和預警,是應對復雜工作條件下預防邊坡發(fā)生極端事件的有效措施[1-3]。而極值理論不僅提供了建立模型描述極端事件的理論基礎,且具有超越樣本數(shù)據(jù)的能力,有效地對隨機序列的最大(小)值的概率分布和數(shù)據(jù)序列的邊際概率分布尾部進行建模,進而建立有效的風險防范和預警系統(tǒng)[4],在實際工程監(jiān)控領域得到了廣泛的應用。

    傳統(tǒng)的分塊樣本極大值方法(BMM模型)將一獨立隨機觀測序列依據(jù)時間或者其它的標準分隔為幾個非重疊區(qū)域,在每個區(qū)域中選極大值,以這些極大值構(gòu)成的極值樣本數(shù)據(jù)序列進行分布擬合從而估計預警指標[5],主要用于處理具有明顯季節(jié)性數(shù)據(jù)的極值問題。但受限于取樣數(shù)據(jù)方法,忽略了區(qū)域次大值,往往造成區(qū)間內(nèi)除極值外大量數(shù)據(jù)的浪費,且在原型觀測資料年限較短的情況下,使無法滿足參數(shù)估計的樣本量,增加估計的誤差,顯然依據(jù)這些樣本而擬定的預警指標存在一定的偏差。并且鑒于實測數(shù)據(jù)分布可能存在的厚尾性,即極端事件引發(fā)的極端風險的真實值要比正態(tài)分布的大且更頻繁[6],所以BMM模型中基于正態(tài)分布的假設也同樣低估了尾部極端風險[7]。與之相對,POT模型是一種基于超限值的建模方法,即基于廣義帕累托分布(GPD)對超過某一閾值的所有觀測數(shù)據(jù)進行建模,漸進地刻畫分布的尾部特征。該模型更能充分利用樣本數(shù)據(jù)的極值序列包含的信息,理論上更加符合實際情況,適用于工程實際。因此,本文基于POT模型,以超過某一閾值的數(shù)據(jù)序列作為研究對象,利用廣義帕累托分布擬合超限分布,提出了改進的閾值估計方法,并利用極大似然估計確定分布參數(shù),建立了改進POT模型,并應用于邊坡變形監(jiān)測資料分析及預警,結(jié)合邊坡工程的實例對兩種模型所確定的預警指標進行比較,驗證了方法的合理性。

    1 BMM模型

    傳統(tǒng)的分塊樣本極大值方法(BMM)是基于廣義極值分布的參數(shù)模型,它是對某一同分布樣本分塊后的極大值進行建模。Fisher和Tippett[8]證明了如果獨立同分布序列的標準極大值收斂于某分布,那么其極限分布就是廣義極值分布(GEV)。

    以上基于超限值數(shù)據(jù)序列的改進POT模型旨在針對傳統(tǒng)BMM模型在取樣數(shù)據(jù)方法上僅選取極值的缺陷,將超過某一閾值的較大測值均納入考慮,可使所擬定的預警指標更貼近實際,有效規(guī)避風險。下面以某邊坡工程為例,利用其變形監(jiān)測資料建立以上兩種模型,并擬定預警指標進行分析比較,驗證改進POT模型的合理性。

    3 實例分析

    某水電站位于我國西南部,工程以發(fā)電為主,樞紐主要由混凝土面板堆石壩、右岸開敞式溢洪道、左岸引水發(fā)電系統(tǒng)、地面廠房等組成。壩頂高程643.0 m,最大壩高約135 m,壩頂長度315 m,水庫正常蓄水位639.0 m,可調(diào)節(jié)庫容3.341億m3。樞紐工程位于高山峽谷、陡峻岸坡地段,工程開挖形成了多處高開挖邊坡,主要有:壩址左岸邊坡、左岸下游邊坡、右岸高程643.0 m以上邊坡(主要屬于溢洪道引渠段邊坡)和溢洪道邊坡。

    高開挖邊坡主要集中在溢洪道所在的右岸邊坡。該范圍邊坡呈弧形,坡高一般60~80 m,末端高達114 m,主要為土質(zhì)邊坡和部分巖質(zhì)邊坡。巖質(zhì)邊坡存在不利結(jié)構(gòu)面組合構(gòu)成的不穩(wěn)定楔形塊體,施工過程中產(chǎn)生了兩次塌滑,穩(wěn)定性較差。土質(zhì)邊坡采用混凝土框架梁、節(jié)點加錨索固坡形式,由于支護滯后等原因,沿覆蓋層與基巖接觸面產(chǎn)生了滑塌。之后設計方案調(diào)整為:放緩開挖邊坡為1∶1.2~1∶1.0;上部采用1∶1.4,并用M10貼坡漿砌石護坡,厚度40 cm;下部邊坡及變坡部位采用C20貼坡鋼筋混凝土加錨索固坡;在坡面上均按梅花型布設排水孔。為了確保右岸邊坡的穩(wěn)定性,在右岸高程643.0 m以上邊坡布置了錨索測力計、表面變形監(jiān)測點、多點位移計及測斜孔等邊坡監(jiān)測設施。

    鑒于右岸邊坡工程的重要性,基于極值理論分析邊坡安全監(jiān)測資料,擬定某一置信水平下的預警指標,以便展開對右岸邊坡的實時監(jiān)控和預警。取右岸邊坡有效表面變形監(jiān)測點BPB-TP-07,各測點位置見圖1。自2005年1月16日至2014年7月17日順河向表面位移的測值序列,樣本總數(shù)為n=127,分別利用BMM模型和改進POT模型擬定右岸邊坡在95%置信水平下的水平位移預警值。

    3.1 利用BMM模型擬定預警值

    將所取測值序列以年為區(qū)間進行分塊處理,得到m=10的年極大值樣本序列{M1,M2,…,Mm},利用式(3)對分布參數(shù)ξ,μ,σ通過matlab編程得到極大似然估計值(結(jié)果見表1),代入GEV分布的概率密度函數(shù),易得在α=95%的置信水平下邊坡水平位移的預警值為201.247 mm。

    3.2 利用改進POT模型擬定預警值

    基于改進的Hill估計方法確定閾值,然后對分布參數(shù)采用極大似然估計,從而得到超限值的條件分布函數(shù)來表示總體分布函數(shù)。

    首先,將測值序列{x1,x2,…,xn}按升序排列得到次序統(tǒng)計量x(i),由式(8)通過matlab編程獲得k及對應的估計量γ(k),繪制點(k,γ(k)-1)集合的Hill圖見圖2。

    3.3 預警指標的比較分析

    通過算例發(fā)現(xiàn),在相同置信水平下利用超限值應用廣義帕累托分布擬合得到的預警指標小于利用塊極大值應用正態(tài)分布得到的預警指標。為排除單個測點測值序列所具有的偶然性,用同樣方法隨機取右岸邊坡若干有效測點進行計算同樣得到類似結(jié)論。由此可知,在實際監(jiān)測邊坡安全狀況的過程中,基于傳統(tǒng)BMM模型所得的預警指標往往可能偏于開放,尤其是在處理數(shù)據(jù)序列年限較短、有效區(qū)間次極大值被忽略的情況下,所得的預警指標存在較大偏差。而利用改進POT模型,可以充分考慮測值序列中的較大測值,所得的預警指標也更具參考性,使有效規(guī)避極端情況發(fā)生的風險。

    4 結(jié)語

    (1)借助極值理論和BMM模型提出了改進POT模型,研究了邊坡預警指標的擬定原理和算法。并結(jié)合具體邊坡工程的實例,分別利用區(qū)間極大值和超限值對邊坡測值序列進行整體分布擬合,從而比較同一置信水平下的預警指標,得到改進POT模型相較BMM模型所得預警指標更偏于保守,對于邊坡工程的風險防范和極端情況預警更具指導性。

    (2)改進POT模型通過引入斜率變點來界定Hill圖中序列的平穩(wěn)性,有效避免了主觀判定閾值產(chǎn)生的誤差,并通過編程使模型適用于大樣本數(shù)據(jù)的自動處理,為其他類似工程預警系統(tǒng)的建立和預警指標的擬定提供了參考。

    (3)基于一維情況對測值序列的極值進行建模處理來擬定預警指標,僅考慮了超越閾值的大小,沒有考慮隨機變量之間的相關(guān)性,因此可考慮用多維極值理論來研究其可行性。

    參考文獻(References):

    [1] 吳中如.水工建筑物安全監(jiān)控理論及其應用[M].北京:高等教育出版社,2003.(WU Zhong-ru.The safety of hydraulic structure monitoring theory and its application[M].Beijing:Higher Education Press,2003.(in Chinese))

    [2] Su H,Wen Z,Wu Z.Study on an intelligent inference engine in early-warning system of dam health[J].Water Resources Management,2011,25(6):1545-1563.

    [3] 蘇懷智,王鋒,劉紅萍.基于POT模型建立大壩服役性態(tài)預警指標[J].水利學報,2012(08):974-978,986.(SU Huai-zhi,WANG Feng,LIU Hong-ping.Early-warning index for dam service behavior based on POT model[J].Journal of Hydraulic Engineering,2012(8):974-978,986.(in Chinese))

    [4] 高娟.極值分布參數(shù)估計方法的研究[D].北京:華北電力大學,2009.(GAO Juan.The method study of extreme value distribution parameter estimation[D].Beijing:North China Electric Power University,2009.(in Chinese))

    [5] 劉莎莎.POT模型閾值的選取及應用[D].長春:吉林大學,2014.(LIU Sha-sha.The selection and application of the threshold of POT model[D].Changchun: Jilin University,2014.(in Chinese))

    [6] 李賀,葉中行.極值理論與VaR計算[J].寧夏大學學報:自然科學版,2007(02):124-127.(LI He,YE Zhong-xing.Extreme value theory and VaR computation[J].Journal of Ningxia University:Natural Science Edition,2007(02):124-127.(in Chinese))

    [7] 周敏娟,蘇鑫.極值理論POT模型與閾值選取研究[J].中國證券期貨,2011(6):199-200.(ZHOU Min-juan,SU Xin.Extreme value theory POT model and study on threshold[J].Securities and Futures of China,2011(6):199-200.(in Chinese))

    [8] Fisher R A,Tippett L H C.Limiting forms of the frequency distribution of the largest or smallest member of a sample[C].//Mathematical Proceedings of the Cambridge Philosophical Society.Cambridge University Press,1928,24(02):180-190.

    [9] Jenkinson A F.The frequency distribution of the annual maximum(or minimum)values of meteorological elements[J].Quarterly Journal of the Royal Meteorological Society,1955,81(348):158-171.

    [10] Leadbetter M R,Lindgren G,Rootzén H.Extremes and related properties of random sequences and processes[M].New York:Springer-Verlag,1983.

    [11] Pickands III J.Statistical inference using extreme order statistics[J].The Annals of Statistics,1975:119-131.

    [12] McNeil A J,F(xiàn)rey R.Estimation of tail-related risk measures for heteroscedastic financial time series:an extreme value approach[J].Journal of empirical finance,2000,7(3):271-300.

    [13] 桂文林.基于極端值理論(EVT)的金融風險度量[D].廣州:暨南大學,2005.(GUI Wen-lin.Calculate the financial risk according to Extreme Value Theory[D].Guangzhou:Jinan University,2005.(in Chinese))

    [14] Qi Y E W Y M B.Using the Improved Hill Estimator Model to Evaluate VaR[J].Journal of the Graduate School of the Chinese Academy of Science,2004,3:002.

    [15] Dacorogna M M,Müller U A,Pictet O V,et al.The distribution of extremal foreign exchange rate returns in extremely large data sets[M].Netherlands:Tinbergen Institute,1995.

    [16] 宋加山,李勇,彭誠,等.極值理論中閾值選取的Hill估計方法改進[J].中國科學技術(shù)大學學報,2008(9):1104-1108.(SONG Jia-shan,LI Yong,PENG Cheng,et al.Improving the traditional Hill estimator to select the threshold in extreme value theory[J].Journal of University of Science and Technology of China,2008(9):1104-1108.(in Chinese))

    [17] 葉五一,繆柏其,譚常春.基于分位點回歸模型變點檢測的金融傳染分析[J].數(shù)量經(jīng)濟技術(shù)經(jīng)濟研究,2007(10):151-161.(YE Wu-yi,LIAO Bo-qi,TAN Chang-chun.Analysis of financial contagion based on change point testing of quartile regression model[J].The Journal of Quantitative and Technical Economics,2007(10):151-161.(in Chinese))

    [18] Smith R L.Estimating tails of probability distributions[J].The annals of Statistics,1987:1174-1207.

    猜你喜歡
    閾值邊坡
    這邊山 那邊坡
    小波閾值去噪在深小孔鉆削聲發(fā)射信號處理中的應用
    基于自適應閾值和連通域的隧道裂縫提取
    比值遙感蝕變信息提取及閾值確定(插圖)
    河北遙感(2017年2期)2017-08-07 14:49:00
    水利水電工程高邊坡的治理與加固探討
    室內(nèi)表面平均氡析出率閾值探討
    SNS主動柔性防護網(wǎng)在高邊坡防護中的應用
    巖土高邊坡的設計方案
    基于SLOPE/W的邊坡穩(wěn)定分析
    基于不同軟件對高邊坡穩(wěn)定計算
    日韩欧美 国产精品| 十分钟在线观看高清视频www| 成人永久免费在线观看视频| 日韩欧美国产在线观看| 婷婷丁香在线五月| 国产真实乱freesex| 18禁裸乳无遮挡免费网站照片 | 亚洲欧美一区二区三区黑人| 国产成人精品久久二区二区91| 中国美女看黄片| 最近在线观看免费完整版| 操出白浆在线播放| 老司机午夜福利在线观看视频| 亚洲精品国产一区二区精华液| 欧美大码av| 亚洲午夜理论影院| 99久久99久久久精品蜜桃| 日本在线视频免费播放| 满18在线观看网站| 国产精品日韩av在线免费观看| 久久精品影院6| 97人妻精品一区二区三区麻豆 | 一区福利在线观看| 久久精品夜夜夜夜夜久久蜜豆 | 精品人妻1区二区| 男女下面进入的视频免费午夜 | 国产色视频综合| 极品教师在线免费播放| 亚洲九九香蕉| 国产不卡一卡二| 少妇裸体淫交视频免费看高清 | 校园春色视频在线观看| 老司机在亚洲福利影院| 午夜福利一区二区在线看| 亚洲七黄色美女视频| www国产在线视频色| 在线观看舔阴道视频| 欧美又色又爽又黄视频| 国产精品av久久久久免费| 久久这里只有精品19| 中文字幕久久专区| 亚洲精品在线美女| 国产精品日韩av在线免费观看| 欧美日韩福利视频一区二区| 午夜日韩欧美国产| 好男人电影高清在线观看| 哪里可以看免费的av片| 国产精品美女特级片免费视频播放器 | 变态另类丝袜制服| 亚洲一卡2卡3卡4卡5卡精品中文| 中文字幕高清在线视频| 亚洲欧美日韩无卡精品| 18禁美女被吸乳视频| 首页视频小说图片口味搜索| 757午夜福利合集在线观看| 亚洲自拍偷在线| 国产精品av久久久久免费| 久久久久久亚洲精品国产蜜桃av| 欧美性长视频在线观看| 不卡av一区二区三区| 精品国产国语对白av| 91成年电影在线观看| 少妇熟女aⅴ在线视频| 亚洲第一电影网av| 男男h啪啪无遮挡| 亚洲熟女毛片儿| 嫩草影视91久久| 久久狼人影院| 亚洲片人在线观看| 国内毛片毛片毛片毛片毛片| 真人一进一出gif抽搐免费| 国产亚洲欧美精品永久| 国产熟女xx| 欧美日韩精品网址| 少妇 在线观看| 色综合婷婷激情| 国产精品乱码一区二三区的特点| 欧美日本亚洲视频在线播放| 在线天堂中文资源库| 精品久久久久久久毛片微露脸| 亚洲欧美一区二区三区黑人| 18美女黄网站色大片免费观看| 国产麻豆成人av免费视频| 欧美国产日韩亚洲一区| 午夜日韩欧美国产| 亚洲人成网站在线播放欧美日韩| 18美女黄网站色大片免费观看| 国产精品久久视频播放| 中文字幕精品亚洲无线码一区 | 亚洲欧美精品综合久久99| 脱女人内裤的视频| 一二三四在线观看免费中文在| 久久精品国产清高在天天线| 久久九九热精品免费| 中文字幕久久专区| 日韩成人在线观看一区二区三区| 午夜福利在线观看吧| 级片在线观看| 91字幕亚洲| 国产精品爽爽va在线观看网站 | 麻豆久久精品国产亚洲av| 久久精品影院6| 无人区码免费观看不卡| 中文资源天堂在线| 国产精品久久久人人做人人爽| 亚洲一码二码三码区别大吗| 国产单亲对白刺激| 香蕉久久夜色| 久久人人精品亚洲av| 91av网站免费观看| 最近在线观看免费完整版| 黄片播放在线免费| 欧美精品啪啪一区二区三区| 国产av在哪里看| 啦啦啦 在线观看视频| 成人手机av| 99热这里只有精品一区 | 国产免费av片在线观看野外av| 香蕉丝袜av| x7x7x7水蜜桃| 亚洲av片天天在线观看| 大型黄色视频在线免费观看| 一进一出好大好爽视频| 成人午夜高清在线视频 | 一区福利在线观看| 亚洲一码二码三码区别大吗| 欧美一级毛片孕妇| 好看av亚洲va欧美ⅴa在| 悠悠久久av| 国产欧美日韩一区二区三| 91av网站免费观看| 亚洲成人久久爱视频| 国产av一区在线观看免费| 一本大道久久a久久精品| 久久久久久人人人人人| 香蕉av资源在线| 欧美成人性av电影在线观看| 日本一区二区免费在线视频| 观看免费一级毛片| 亚洲国产精品合色在线| 97碰自拍视频| 亚洲熟妇熟女久久| 久久性视频一级片| 国产久久久一区二区三区| 日本黄色视频三级网站网址| 亚洲性夜色夜夜综合| 色av中文字幕| 亚洲精品在线美女| 国产黄a三级三级三级人| 十分钟在线观看高清视频www| 欧美日韩精品网址| 精品久久久久久久毛片微露脸| 99久久综合精品五月天人人| 国产成人欧美| 久久精品国产亚洲av高清一级| www.熟女人妻精品国产| 波多野结衣高清无吗| 国产高清有码在线观看视频 | 久久久久久久久中文| 18禁黄网站禁片免费观看直播| 午夜福利在线观看吧| 99精品在免费线老司机午夜| 亚洲片人在线观看| 亚洲人成伊人成综合网2020| 黑人巨大精品欧美一区二区mp4| 老汉色av国产亚洲站长工具| 午夜久久久在线观看| 精品一区二区三区视频在线观看免费| 搡老妇女老女人老熟妇| 啦啦啦韩国在线观看视频| av中文乱码字幕在线| av视频在线观看入口| 91大片在线观看| 久久久久亚洲av毛片大全| av视频在线观看入口| 精品一区二区三区av网在线观看| 日韩成人在线观看一区二区三区| 久久久国产成人免费| 一本大道久久a久久精品| 亚洲激情在线av| 在线观看免费午夜福利视频| 亚洲色图av天堂| 亚洲av成人av| www.精华液| 免费在线观看影片大全网站| 免费人成视频x8x8入口观看| 久久久久久久久中文| 午夜a级毛片| 香蕉国产在线看| 91在线观看av| 十八禁网站免费在线| 首页视频小说图片口味搜索| 日本免费一区二区三区高清不卡| 一边摸一边抽搐一进一小说| 欧美绝顶高潮抽搐喷水| 亚洲av美国av| 校园春色视频在线观看| 国产男靠女视频免费网站| 亚洲国产欧美一区二区综合| 香蕉久久夜色| 国产精品日韩av在线免费观看| 国产亚洲欧美精品永久| 人人妻,人人澡人人爽秒播| 1024视频免费在线观看| 人人妻人人看人人澡| 国产精品香港三级国产av潘金莲| 在线看三级毛片| 亚洲五月色婷婷综合| or卡值多少钱| √禁漫天堂资源中文www| 日韩免费av在线播放| 亚洲一码二码三码区别大吗| 日日摸夜夜添夜夜添小说| 啦啦啦免费观看视频1| 精品欧美国产一区二区三| 超碰成人久久| 午夜福利免费观看在线| 久久伊人香网站| 一级片免费观看大全| 久久久国产精品麻豆| 亚洲激情在线av| 香蕉国产在线看| 国产熟女xx| 男女下面进入的视频免费午夜 | 日本一区二区免费在线视频| 国产在线精品亚洲第一网站| 动漫黄色视频在线观看| 99热只有精品国产| 后天国语完整版免费观看| 成人一区二区视频在线观看| 老司机靠b影院| 亚洲精品美女久久av网站| 欧美日韩一级在线毛片| 可以在线观看毛片的网站| 999久久久国产精品视频| 一个人免费在线观看的高清视频| 99热6这里只有精品| 欧美日本视频| 亚洲欧美日韩无卡精品| 最近最新中文字幕大全免费视频| 成人亚洲精品av一区二区| 熟女电影av网| 制服丝袜大香蕉在线| 黄色a级毛片大全视频| 亚洲国产精品999在线| 久9热在线精品视频| 日韩精品中文字幕看吧| 欧美性猛交╳xxx乱大交人| 九色国产91popny在线| 亚洲成人精品中文字幕电影| 热99re8久久精品国产| 久久国产精品影院| 两个人视频免费观看高清| 看片在线看免费视频| 在线永久观看黄色视频| 十分钟在线观看高清视频www| 不卡一级毛片| 久久精品成人免费网站| 国产黄片美女视频| 亚洲一码二码三码区别大吗| 老司机深夜福利视频在线观看| 在线天堂中文资源库| 国内久久婷婷六月综合欲色啪| 亚洲美女黄片视频| aaaaa片日本免费| 久久久久久免费高清国产稀缺| 欧美成人性av电影在线观看| 欧美乱色亚洲激情| 俄罗斯特黄特色一大片| 看片在线看免费视频| 国产精品免费一区二区三区在线| 精品欧美一区二区三区在线| 好男人电影高清在线观看| 高清毛片免费观看视频网站| 精品福利观看| 亚洲精品中文字幕一二三四区| 一级毛片精品| 波多野结衣高清作品| 亚洲男人的天堂狠狠| 欧美又色又爽又黄视频| tocl精华| 夜夜看夜夜爽夜夜摸| 啪啪无遮挡十八禁网站| 搞女人的毛片| 日本免费a在线| 中文字幕av电影在线播放| 色综合亚洲欧美另类图片| 久久久久国产一级毛片高清牌| 日本免费一区二区三区高清不卡| 国产色视频综合| 亚洲色图av天堂| 日本免费a在线| 看黄色毛片网站| 免费高清视频大片| 91老司机精品| 国产私拍福利视频在线观看| 伦理电影免费视频| 国产精品久久久久久精品电影 | 88av欧美| 久久伊人香网站| 亚洲第一av免费看| 十分钟在线观看高清视频www| 精品电影一区二区在线| 久99久视频精品免费| 国产亚洲精品久久久久久毛片| cao死你这个sao货| 亚洲最大成人中文| 最近最新中文字幕大全免费视频| 搡老熟女国产l中国老女人| 久久久水蜜桃国产精品网| 观看免费一级毛片| 免费在线观看完整版高清| 色尼玛亚洲综合影院| 一级a爱片免费观看的视频| 99久久精品国产亚洲精品| 悠悠久久av| 国产精华一区二区三区| 久久草成人影院| 国产伦一二天堂av在线观看| 久久狼人影院| 久久精品成人免费网站| 成人一区二区视频在线观看| 午夜福利18| 欧美成人一区二区免费高清观看 | 久久久久久久精品吃奶| 淫妇啪啪啪对白视频| 久久人人精品亚洲av| 亚洲精品久久国产高清桃花| 日日摸夜夜添夜夜添小说| 一级a爱片免费观看的视频| 黄片大片在线免费观看| 国产成人精品久久二区二区免费| 亚洲第一av免费看| 青草久久国产| 又黄又粗又硬又大视频| 亚洲熟妇熟女久久| 日日摸夜夜添夜夜添小说| 天堂√8在线中文| 精品一区二区三区av网在线观看| 亚洲熟妇中文字幕五十中出| 欧美最黄视频在线播放免费| 哪里可以看免费的av片| 黄色毛片三级朝国网站| av在线天堂中文字幕| 一区二区三区精品91| 怎么达到女性高潮| 亚洲第一av免费看| 色精品久久人妻99蜜桃| 亚洲天堂国产精品一区在线| 国产99白浆流出| 免费看美女性在线毛片视频| 日本三级黄在线观看| 老熟妇仑乱视频hdxx| 欧美zozozo另类| 中文字幕av电影在线播放| 国产高清激情床上av| 可以在线观看的亚洲视频| 成人国产综合亚洲| 国产一区二区激情短视频| 成熟少妇高潮喷水视频| 999久久久国产精品视频| 久久亚洲真实| 黄色毛片三级朝国网站| 国产免费av片在线观看野外av| 夜夜夜夜夜久久久久| 男女午夜视频在线观看| 久久久国产成人精品二区| 日韩欧美国产一区二区入口| 欧美丝袜亚洲另类 | 精品一区二区三区视频在线观看免费| 1024手机看黄色片| 18禁观看日本| 国产日本99.免费观看| 婷婷精品国产亚洲av在线| 无人区码免费观看不卡| 波多野结衣高清无吗| 中文字幕精品亚洲无线码一区 | 午夜福利在线观看吧| 欧美国产精品va在线观看不卡| 国产伦人伦偷精品视频| 婷婷亚洲欧美| 欧美成人一区二区免费高清观看 | www日本在线高清视频| 亚洲一码二码三码区别大吗| 侵犯人妻中文字幕一二三四区| 久久精品成人免费网站| 一边摸一边做爽爽视频免费| 久久午夜亚洲精品久久| 亚洲 国产 在线| 国产成人av教育| 韩国精品一区二区三区| 国产亚洲欧美98| 91成人精品电影| www国产在线视频色| 精华霜和精华液先用哪个| 国产v大片淫在线免费观看| 亚洲第一av免费看| 国产av一区在线观看免费| 宅男免费午夜| 可以在线观看毛片的网站| 久久久久久久午夜电影| 国产熟女xx| 亚洲欧美一区二区三区黑人| 亚洲精品一卡2卡三卡4卡5卡| 又大又爽又粗| 999精品在线视频| 在线观看免费午夜福利视频| 国产亚洲精品第一综合不卡| 校园春色视频在线观看| 免费电影在线观看免费观看| 国产激情欧美一区二区| 国内毛片毛片毛片毛片毛片| 老司机在亚洲福利影院| 给我免费播放毛片高清在线观看| 久久久久久九九精品二区国产 | 男男h啪啪无遮挡| 欧美性猛交╳xxx乱大交人| 亚洲自偷自拍图片 自拍| 亚洲aⅴ乱码一区二区在线播放 | 国产精品二区激情视频| 一区二区三区国产精品乱码| 精品免费久久久久久久清纯| 国产精品国产高清国产av| 一级片免费观看大全| 国产亚洲精品久久久久久毛片| 一本一本综合久久| av天堂在线播放| 性欧美人与动物交配| 国产精品香港三级国产av潘金莲| 成人18禁高潮啪啪吃奶动态图| 中文字幕人妻丝袜一区二区| 一区二区三区国产精品乱码| 日本免费a在线| 成人国产综合亚洲| 大香蕉久久成人网| 麻豆久久精品国产亚洲av| 国产真实乱freesex| 精品人妻1区二区| 成年免费大片在线观看| 男女下面进入的视频免费午夜 | 一级a爱片免费观看的视频| 国产精品久久电影中文字幕| 亚洲久久久国产精品| 中文字幕久久专区| 欧美日韩中文字幕国产精品一区二区三区| 美国免费a级毛片| 色播亚洲综合网| 麻豆成人av在线观看| 久久久久国内视频| 日韩欧美 国产精品| 桃色一区二区三区在线观看| 久久精品国产亚洲av香蕉五月| 最近最新中文字幕大全电影3 | 成年人黄色毛片网站| 精品一区二区三区视频在线观看免费| 美女 人体艺术 gogo| 久久国产乱子伦精品免费另类| 成人亚洲精品av一区二区| 国产91精品成人一区二区三区| 窝窝影院91人妻| 精品国产乱码久久久久久男人| 日韩精品中文字幕看吧| 香蕉国产在线看| 亚洲精品在线观看二区| 精品日产1卡2卡| 黄片小视频在线播放| 国产精品野战在线观看| 满18在线观看网站| 国产99白浆流出| 免费在线观看成人毛片| 日韩欧美一区视频在线观看| 我的亚洲天堂| 天天添夜夜摸| 一级a爱视频在线免费观看| 两性夫妻黄色片| 婷婷精品国产亚洲av在线| 99热这里只有精品一区 | 国产精品99久久99久久久不卡| 午夜精品在线福利| 色在线成人网| 最近最新中文字幕大全电影3 | 欧美一级a爱片免费观看看 | 中出人妻视频一区二区| 亚洲av美国av| 久久 成人 亚洲| 美女高潮到喷水免费观看| 亚洲精品av麻豆狂野| АⅤ资源中文在线天堂| 精品久久久久久久人妻蜜臀av| 亚洲无线在线观看| 午夜日韩欧美国产| 露出奶头的视频| 在线观看舔阴道视频| 免费看美女性在线毛片视频| 色婷婷久久久亚洲欧美| 99国产精品一区二区蜜桃av| 欧美日韩精品网址| 久久精品成人免费网站| 亚洲 欧美 日韩 在线 免费| 午夜精品在线福利| 欧美一区二区精品小视频在线| 午夜激情福利司机影院| 久久久久亚洲av毛片大全| 99热6这里只有精品| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av电影不卡..在线观看| 两个人免费观看高清视频| 亚洲精品国产一区二区精华液| 久久狼人影院| 国产高清videossex| 亚洲精品久久成人aⅴ小说| 亚洲av电影不卡..在线观看| 一a级毛片在线观看| 在线播放国产精品三级| 88av欧美| 久久久久久久久中文| 久久人妻福利社区极品人妻图片| 国产伦一二天堂av在线观看| 国产伦人伦偷精品视频| 久久久久国产精品人妻aⅴ院| 国产成人精品无人区| 亚洲专区国产一区二区| 久久久久九九精品影院| 不卡一级毛片| 日本一区二区免费在线视频| 97碰自拍视频| 免费在线观看成人毛片| 一进一出好大好爽视频| 亚洲国产日韩欧美精品在线观看 | 日韩欧美 国产精品| 特大巨黑吊av在线直播 | 日韩国内少妇激情av| 亚洲成av人片免费观看| 精品久久久久久,| 亚洲人成网站高清观看| 精品国产乱码久久久久久男人| 亚洲国产高清在线一区二区三 | 99久久久亚洲精品蜜臀av| 老司机福利观看| 黄色视频不卡| 欧美日韩精品网址| 国产成人av教育| 国产精品久久久人人做人人爽| 精品高清国产在线一区| 视频在线观看一区二区三区| 午夜福利在线在线| 国产精品永久免费网站| 欧美日韩中文字幕国产精品一区二区三区| 国产黄片美女视频| 午夜激情福利司机影院| 久久天堂一区二区三区四区| 韩国av一区二区三区四区| 十八禁人妻一区二区| 村上凉子中文字幕在线| 91老司机精品| 亚洲一码二码三码区别大吗| 午夜福利高清视频| 色尼玛亚洲综合影院| 免费高清在线观看日韩| 亚洲成人免费电影在线观看| 女人被狂操c到高潮| 亚洲精品色激情综合| 色婷婷久久久亚洲欧美| 日韩有码中文字幕| 久久99热这里只有精品18| 久久久久久久久中文| a在线观看视频网站| 在线视频色国产色| 脱女人内裤的视频| 在线视频色国产色| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品中文字幕一二三四区| 日韩欧美 国产精品| 亚洲国产精品久久男人天堂| 成在线人永久免费视频| 国产午夜福利久久久久久| 亚洲欧美精品综合久久99| 国产精品日韩av在线免费观看| 久久国产精品影院| 女人被狂操c到高潮| 久久久水蜜桃国产精品网| 欧美精品亚洲一区二区| 很黄的视频免费| 88av欧美| 老熟妇乱子伦视频在线观看| 免费在线观看成人毛片| 日本精品一区二区三区蜜桃| 亚洲午夜精品一区,二区,三区| 非洲黑人性xxxx精品又粗又长| 国产精品免费一区二区三区在线| 日本熟妇午夜| 国产成人av教育| 99热只有精品国产| 色老头精品视频在线观看| 精品无人区乱码1区二区| 免费电影在线观看免费观看| 欧美精品啪啪一区二区三区| 热99re8久久精品国产| 日韩欧美三级三区| 美女国产高潮福利片在线看| 精品一区二区三区四区五区乱码| 亚洲精品色激情综合| 白带黄色成豆腐渣| xxxwww97欧美| 哪里可以看免费的av片| 麻豆成人av在线观看| 精品久久蜜臀av无| 久久久国产精品麻豆| 妹子高潮喷水视频| 亚洲男人的天堂狠狠| 亚洲九九香蕉| 国产爱豆传媒在线观看 | 欧美亚洲日本最大视频资源| 久久伊人香网站| 制服丝袜大香蕉在线|