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

    蒙特卡洛與馬爾可夫方法在降水預(yù)測中的應(yīng)用

    2020-10-20 06:37:11黃鶴張維江李娟
    人民黃河 2020年5期
    關(guān)鍵詞:原州區(qū)

    黃鶴 張維江 李娟

    摘 要:為了更深入地分析原州區(qū)的降水特征,為該區(qū)域水資源預(yù)測提供依據(jù),基于原州區(qū)1957—2016年60 a降水資料,采用蒙特卡洛方法推求降水分布,采用K-S檢驗(yàn)對模型進(jìn)行顯著性檢驗(yàn),采用基于歐氏距離的層次聚類方法進(jìn)行狀態(tài)劃分,確定了原州區(qū)的降水分布,建立了適用于原州區(qū)的滑動(dòng)平均加權(quán)馬爾可夫預(yù)測模型。根據(jù)已有數(shù)據(jù)驗(yàn)證了預(yù)測結(jié)果的有效性,再結(jié)合已確定的降水分布,通過K-S檢驗(yàn),檢驗(yàn)了未來5 a降水預(yù)測的準(zhǔn)確性。結(jié)果表明:原州區(qū)降水分布符合P-Ⅲ型分布;馬爾可夫模型適用于原州區(qū)降水預(yù)測,且未來5 a的降水預(yù)測結(jié)果是準(zhǔn)確的,分別為508.5、520.8、554.9、451.0、466.6 mm。

    關(guān)鍵詞:蒙特卡洛方法;馬爾可夫模型;隨機(jī)模擬;降水預(yù)測;原州區(qū);K-S檢驗(yàn)

    中圖分類號:TV11 文獻(xiàn)標(biāo)志碼:A

    doi:10.3969/j.issn.1000-1379.2020.05.004

    Abstract: In order to analyze the characteristics of precipitation in YuanzhouDistrict and provide a basis for water resources prediction in the region, based on the precipitation data from 1957 to 2016 of the district, the Monte Carlo method was used to estimate the precipitation distribution and the K-S test was used to model the significant test. It used the hierarchical clustering method based on Euclidean distance to divide the state, determined the precipitation distribution of the region and established a sliding average weighted Markov prediction model which was suitable for Yuanzhou District. Based on the existing data, the true validity of the prediction results was verified. Combined with the determined precipitation distribution, the accuracy of precipitation prediction for the next 5 years was tested by K-S test. Comparing with Pearson three-type distribution, the results show that the precipitation distribution in Yuanzhou District is more consistent with the log-normal distribution. The Markov model can be applied to the precipitation forecast in the region, and the prediction results of the sliding average precipitation in the next 5 years are real and effective, respectively 508.5, 520.8, 554.9, 451.0 and 466.6 mm.

    Key words: Monte Carlo method; Markov model; stochastic simulation; precipitation prediction; Yuanzhou District; K-S test

    1 引 言

    由于客觀世界中的一些現(xiàn)象可能與另一種現(xiàn)象存在著某種相似性,因此我們經(jīng)常從一種現(xiàn)象出發(fā)來研究另一種現(xiàn)象。當(dāng)某個(gè)概率模型可以描述隨機(jī)系統(tǒng)并可以基于此概率模型進(jìn)行實(shí)驗(yàn)時(shí),這種實(shí)驗(yàn)方法即為隨機(jī)模擬方法[1]。近年來,隨機(jī)模擬在水文系統(tǒng)預(yù)測中發(fā)展迅速,取得了顯著進(jìn)展[2]。劉新立[3]研究了隨機(jī)過程情況下隨機(jī)模擬和馬爾可夫鏈在水災(zāi)風(fēng)險(xiǎn)管理中所起的作用,研究表明兩者結(jié)合可以評估未來若干個(gè)時(shí)間單位水災(zāi)所造成的風(fēng)險(xiǎn);尹正杰等[4]提出了一個(gè)同時(shí)含有趨勢、季節(jié)、隨機(jī)3個(gè)分量的時(shí)間序列模型,并通過此時(shí)間序列模型對灌區(qū)灌溉需水量進(jìn)行隨機(jī)模擬;溫季等[5]探討了集中作物需水量的隨機(jī)模擬及預(yù)測模型,結(jié)果表明隨機(jī)模擬在作物灌溉管理中具有很強(qiáng)的應(yīng)用性及普適性。

    蒙特卡洛方法與馬爾可夫方法在水文系統(tǒng)方面的應(yīng)用發(fā)展十分迅速。ZHANG等[6]研究開發(fā)了基于馬爾可夫鏈-蒙特卡洛的多級因子分析方法,從而更好地對水文模型參數(shù)不確定性進(jìn)行評估;KNIGHTON等[7]通過蒙特卡洛方法將建立的概率分布應(yīng)用于已知模型,用來估計(jì)水文模型的不確定性;劉悅憶等[8]在淮河流域水動(dòng)力學(xué)-水質(zhì)模型的基礎(chǔ)上,使用蒙特卡洛方法隨機(jī)模擬了大量入流數(shù)據(jù)并將其用作模型的輸入條件進(jìn)行計(jì)算,建立了基于蒙特卡洛模擬的水質(zhì)概率預(yù)報(bào)模型;李娟等[9]應(yīng)用滑動(dòng)平均-馬爾可夫模型對固原市隆德縣的降水進(jìn)行預(yù)測,研究證明改進(jìn)后的馬爾可夫模型預(yù)測精度較高;王艷等[10]在傳統(tǒng)的馬爾可夫方法上用最優(yōu)分割法優(yōu)化了分級標(biāo)準(zhǔn);岳遙等[11]提出了一種基于投影距離的處理級別特征值的方法來代替?zhèn)鹘y(tǒng)的處理方法,并將其引入馬爾可夫模型,從而應(yīng)用于對水質(zhì)的定量預(yù)測;李亞斌等[12]用樣本均值-均方差方法對銅川地區(qū)降水量進(jìn)行分級,建立相應(yīng)的馬爾可夫模型進(jìn)行降水預(yù)測。

    降水受氣候和人類活動(dòng)兩者的共同影響,且由于氣候因素本身具有很強(qiáng)的變異性、復(fù)雜性以及多樣性,因此降水系統(tǒng)呈現(xiàn)出十分復(fù)雜的行為特征。過于復(fù)雜的系統(tǒng)難以建立精確的數(shù)學(xué)模型進(jìn)行準(zhǔn)確預(yù)測。筆者運(yùn)用隨機(jī)模擬原理,在確定地區(qū)降水的分布模型后,結(jié)合蒙特卡洛方法以及馬爾可夫方法對固原市原州區(qū)的降水進(jìn)行隨機(jī)模擬,以期預(yù)測未來原州區(qū)降水的發(fā)展變化趨勢,為水資源合理利用及調(diào)控提供依據(jù)。

    2 研究方法

    2.1 蒙特卡洛方法

    蒙特卡洛方法又稱為概率統(tǒng)計(jì)法,是一種基于概率論思想,對隨機(jī)變量進(jìn)行數(shù)理統(tǒng)計(jì)實(shí)驗(yàn)及分布概率模擬,從而近似求解得到預(yù)測值的方法[13]。只要能構(gòu)建出適當(dāng)?shù)哪P?,此方法都能夠進(jìn)行模擬應(yīng)用,其基本框架:①假設(shè)變量X服從某一概率分布;②用隨機(jī)抽樣的方法對概率分布進(jìn)行抽樣從而得到樣本值,一般選擇常用的均勻分布模型,產(chǎn)生(0,1)區(qū)間內(nèi)的隨機(jī)數(shù),再結(jié)合原始數(shù)據(jù)產(chǎn)生服從特定分布的隨機(jī)序列,即樣本值;③確定和選取統(tǒng)計(jì)值;④由統(tǒng)計(jì)量的算術(shù)平均值得到統(tǒng)計(jì)量的估計(jì)值,從而近似求解出預(yù)測值。

    2.2 馬爾可夫方法

    2.2.1 馬爾可夫基本原理

    馬爾可夫模型是基于馬爾可夫鏈建立起來的,馬爾可夫鏈?zhǔn)且环N特殊的隨機(jī)過程[14]。設(shè)有一隨機(jī)運(yùn)動(dòng)的系統(tǒng),它可能處的狀態(tài)記為E0、E1、…、En。這個(gè)系統(tǒng)只可能在時(shí)刻t(t=1,2,…,n)上改變它的狀態(tài)。隨著隨機(jī)運(yùn)動(dòng)的進(jìn)行,定義一列隨機(jī)變量Xn(n=0,1,2,…),其中Xn=k,表示在t=n時(shí),系統(tǒng)狀態(tài)為Ek。

    2.2.2 模糊集理論中的級別特征值

    使用傳統(tǒng)的馬爾可夫只能夠預(yù)測到某個(gè)區(qū)間,無法預(yù)測較為準(zhǔn)確的降水值,在實(shí)際應(yīng)用中作用有限,筆者選用模糊集理論中的級別特征值方法有效地解決了這個(gè)問題。

    假定最大概率的狀態(tài)為i,當(dāng)H>i時(shí),年降水量的預(yù)測值X預(yù)報(bào)=TiHi+0.5;當(dāng)H

    2.3 K-S分布檢驗(yàn)

    K-S檢驗(yàn)是用于檢測一組數(shù)據(jù)的分布與已知分布之間相似程度的一種檢驗(yàn)方法,通過顯著性水平來判別是否相似。

    設(shè)累積頻率曲線的理論分布形式為Fe(x),樣本的累積頻率F1(x)=k/n。其中:n為總觀測數(shù),k為不大于x的次數(shù)。

    3 原州區(qū)降水隨機(jī)模擬分析

    3.1 原州區(qū)降水量理論分布函數(shù)

    根據(jù)我國水文計(jì)算的相關(guān)規(guī)范,我國主要流域的降水量可假定服從P-Ⅲ型分布或?qū)?shù)正態(tài)分布[15]。通過計(jì)算機(jī)軟件用原始數(shù)據(jù)畫頻率曲線發(fā)現(xiàn),當(dāng)變差系數(shù)CV=0.26、偏態(tài)系數(shù)CS=0.52時(shí),曲線模擬效果很好。又由P-Ⅲ型分布特性可知,當(dāng)CS=2CV時(shí),參數(shù)a0=0,此時(shí)的P-Ⅲ型分布就是伽瑪分布,故筆者采用a0=0的P-Ⅲ型分布,即伽瑪分布進(jìn)行隨機(jī)模擬。因?yàn)楣P者主要研究降水的隨機(jī)模擬,需要建立最為適合的模型,所以對兩種分布進(jìn)行對比,選出最貼近原始數(shù)據(jù)的分布,從而建立最優(yōu)模型。

    3.1.1 伽瑪分布

    3.1.2 對數(shù)正態(tài)分布

    3.1.3 參數(shù)確定

    通過原始降水?dāng)?shù)據(jù)可求得式(7)和式(8)的未知參數(shù)。由原始降水?dāng)?shù)據(jù)可知,原州區(qū)近60 a降水量的期望值為442.0 mm,方差為11 659.2 mm,從而求得對數(shù)正態(tài)分布參數(shù)μ為6.064、σ為0.234,伽瑪分布中的形狀參數(shù)α為14.80、尺度參數(shù)β為0.03。

    3.1.4 對比分析

    將原州區(qū)近60 a降水量的經(jīng)驗(yàn)分布分別與對數(shù)正態(tài)函數(shù)理論分布及伽瑪函數(shù)理論分布作對比,如圖1所示??梢钥闯?,樣本數(shù)據(jù)與對數(shù)正態(tài)分布、伽瑪分布的擬合效果均較好,但無法進(jìn)一步判斷最適合的理論分布,故筆者應(yīng)用K-S檢驗(yàn)對兩者進(jìn)行概率分布判斷。

    為了使結(jié)果更加精確,取置信度α=0.01,n=60,其K-S臨界值Dα(n)為0.206 7。當(dāng)假設(shè)數(shù)據(jù)服從對數(shù)正態(tài)分布時(shí),最大偏差D(n)為0.094 1;假設(shè)數(shù)據(jù)服從伽瑪分布時(shí),最大偏差D(n)為0.107 6。通過對比分布的最大偏差與K-S檢驗(yàn)臨界值,發(fā)現(xiàn)兩者均不能拒絕原假設(shè),而且在99%的置信區(qū)間上符合對數(shù)正態(tài)分布及伽瑪分布。由此可知,原州區(qū)降水?dāng)?shù)據(jù)既服從對數(shù)正態(tài)分布,又服從P-Ⅲ型分布。此處選取一般水文規(guī)范所用到的P-Ⅲ型分布進(jìn)行隨機(jī)模擬及預(yù)測研究。

    3.2 降水量蒙特卡洛隨機(jī)模擬

    借助計(jì)算機(jī)編程,引用均勻分布模型,通過蒙特卡洛方法產(chǎn)生(0,1)區(qū)間的隨機(jī)數(shù),將產(chǎn)生的結(jié)果視為隨機(jī)變量概率。為了使模擬更具可靠性,使用伽瑪分布的反函數(shù),根據(jù)原始降水序列的均值和方差,產(chǎn)生1 000 000組隨機(jī)降水序列值。用平方誤差作為評判標(biāo)準(zhǔn),選取與原始降水序列數(shù)據(jù)平方誤差最小的一組作為隨機(jī)模擬降水量結(jié)果,見圖2。

    將原始降水序列依次與通過伽瑪分布函數(shù)隨機(jī)模擬出的降水序列進(jìn)行對比分析??梢钥闯?,在伽瑪分布隨機(jī)模擬出的新降水序列中,整體波動(dòng)較大,部分極值點(diǎn)在時(shí)間趨勢上保持一致。部分重要指標(biāo)的對比見表1,伽瑪分布所模擬出來的均值、標(biāo)準(zhǔn)差等指標(biāo)與原始數(shù)據(jù)相差不大,模擬效果良好,進(jìn)一步證明了原始數(shù)據(jù)服從P-Ⅲ型分布的結(jié)論。

    3.3 馬爾可夫降水預(yù)測驗(yàn)證

    當(dāng)前對馬爾可夫預(yù)測準(zhǔn)確度的評定大都是通過原始數(shù)據(jù)進(jìn)行對比驗(yàn)證,沒有考慮未來數(shù)據(jù)驗(yàn)證的準(zhǔn)確性。而從隨機(jī)模擬結(jié)果得到的P-Ⅲ型分布就可以很好地解決這個(gè)問題。筆者先通過對原始數(shù)據(jù)的預(yù)測及驗(yàn)證判別馬爾可夫方法的可行性,之后再進(jìn)行5 a預(yù)測,并進(jìn)行K-S檢驗(yàn),判別其是否依舊符合P-Ⅲ型分布,從而確定預(yù)測結(jié)果的準(zhǔn)確性。

    3.3.1 分級標(biāo)準(zhǔn)及狀態(tài)確定

    考慮到降水序列的連續(xù)性及降水特性,采取降水序列的3 a滑動(dòng)平均值進(jìn)行降水預(yù)測??紤]水文現(xiàn)象的本身特性及序列數(shù)據(jù)的結(jié)構(gòu)合理性,將降水量序列分為5類,即將降水量劃分為5個(gè)區(qū)間,分別為枯水年、偏枯年、平水年、偏豐年、豐水年,采用基于歐氏距離的層次聚類法進(jìn)行分級,見表2,其中x為年降水量。

    3.3.3 相關(guān)系數(shù)與權(quán)重的確定

    3.3.4 模型檢驗(yàn)及分析

    基于隨機(jī)模擬出的2007—2011年的年降水?dāng)?shù)據(jù),采用加權(quán)馬爾可夫模型對2012年降水量進(jìn)行預(yù)測,結(jié)果見表5,2012年的預(yù)測值為狀態(tài)3,根據(jù)模糊集理論,級別特征值為2.824,2012年的降水量預(yù)測值為500.8 mm,相對誤差為13.9%,在中長期水文預(yù)報(bào)允許誤差20%之內(nèi)。2013年、2014年預(yù)測值分別為521.4、499.5 mm,相對誤差分別為6.5%、5.8%,都在允許誤差以內(nèi)。

    通過以上對原始數(shù)據(jù)的驗(yàn)證,能在一定程度上證明加權(quán)馬爾可夫模型在降水預(yù)測應(yīng)用上的可行性。但是,其只能對未來降水進(jìn)行預(yù)測,并不能對未來降水預(yù)測的準(zhǔn)確性進(jìn)行檢驗(yàn),要確保對未來預(yù)測的準(zhǔn)確性,還需要進(jìn)一步的驗(yàn)證。由于前文已經(jīng)得到寧南地區(qū)降水符合P-Ⅲ型分布,因此筆者通過對寧南山區(qū)降水進(jìn)行接下來3 a滑動(dòng)平均的預(yù)測,再通過K-S檢驗(yàn),看是否滿足P-Ⅲ型分布,預(yù)測結(jié)果見表6。

    此時(shí),n=63,取置信度α=0.01,其K-S臨界值Dα(n)為0.201 8,假設(shè)數(shù)據(jù)服從P-Ⅲ型分布,最大偏差D(n)為0.080 6。通過對比分布的最大偏差與K-S檢驗(yàn)臨界值,發(fā)現(xiàn)不能拒絕原假設(shè),而且在99%的置信區(qū)間上符合P-Ⅲ型分布??梢姡瑢ξ磥? a的降水預(yù)測結(jié)果是可靠的。

    4 結(jié) 論

    (1)依據(jù)原州區(qū)1957—2016年實(shí)測降水資料,運(yùn)用蒙特卡洛方法進(jìn)行隨機(jī)模擬,通過K-S檢驗(yàn),得出原州區(qū)降水符合P-Ⅲ型分布。

    (2)通過對降水?dāng)?shù)據(jù)進(jìn)行3 a滑動(dòng)平均處理,采用基于歐氏距離的層次聚類法進(jìn)行狀態(tài)劃分,確定了枯水年、偏枯年、正常年、偏豐年和豐水年5個(gè)狀態(tài),建立了滑動(dòng)平均-加權(quán)馬爾可夫鏈預(yù)測模型。

    (3)通過對2012年、2013年、2014年滑動(dòng)平均降水量的預(yù)測,得出馬爾可夫降水預(yù)測模型可以應(yīng)用于原州區(qū)降水預(yù)測。

    (4)對未來5 a降水進(jìn)行預(yù)測,預(yù)測值分別為508.5、520.8,554.9、451.0、466.6 mm,通過K-S檢驗(yàn)進(jìn)行驗(yàn)證,證明了預(yù)測的準(zhǔn)確性。

    參考文獻(xiàn):

    [1] 魏艷華,王丙參.概率論與數(shù)理統(tǒng)計(jì)[M].成都:西南交通大學(xué)出版社,2013:264-269.

    [2] 丁晶,鄧育仁.隨機(jī)水文學(xué)[M].成都:成都科技大學(xué)出版社,1988:9-10.

    [3] 劉新立.隨機(jī)過程與隨機(jī)模擬在水災(zāi)風(fēng)險(xiǎn)管理中的應(yīng)用研究[J].經(jīng)濟(jì)科學(xué),2003(1):114-119.

    [4] 尹正杰,袁宏源,崔遠(yuǎn)來,等.灌區(qū)灌溉需水量的隨機(jī)模擬[J].中國農(nóng)村水利水電,2001(11):19-22.

    [5] 溫季,郭樹龍,盧聞航.作物灌溉隨機(jī)模擬技術(shù)研究[J].人民黃河,2004,26(5):39-41.

    [6] ZHANG J L, LI Y P, HUANG G H, et al. Assessment of Parameter Uncertainty in Hydrological Model Using a Markov-Chain-Monte-Carlo-Based Multilevel-Factorial-Analysis Method[J].Journal of Hydrology, 2016, 538(7):471-486.

    [7] KNIGHTON J, WHITE E, LENNON E, et al. Development of Probability Distributions for Urban Hydrologic Model Parameters and a Monte Carlo Analysis of Model Sensitivity[J].Hydrological Processes, 2015, 28(19):5131-5139.

    [8] 劉悅憶,趙建世,黃躍飛,等.基于蒙特卡洛模擬的水質(zhì)概率預(yù)報(bào)模型[J].水利學(xué)報(bào),2015,46(1):51-57.

    [9] 李娟,張維江,馬軼.滑動(dòng)平均-馬爾可夫模型在降水預(yù)測中的應(yīng)用[J].水土保持研究,2005,12(6):200-202,209.

    [10] 王艷,毛明志,范晶,等.最優(yōu)分割法確定的加權(quán)馬爾可夫鏈在降雨量預(yù)測中的應(yīng)用[J].統(tǒng)計(jì)與決策,2009,25(11):17-18.

    [11] 岳遙,李天宏.基于模糊集理論的馬爾可夫模型在水質(zhì)定量預(yù)測中的應(yīng)用[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào),2011,19(2):231-242.

    [12] 李亞斌,徐盼盼,錢會,等.加權(quán)馬爾可夫鏈在銅川地區(qū)降水量預(yù)測中的應(yīng)用[J].灌溉排水學(xué)報(bào),2017,36(5):96-102.

    [13] GARCA-ALONSO C R, ARENAS-ARROYO E, PREZ-ALCAL G M. A Macro-Economic Model to Forecast Remittances Based on Monte-Carlo Simulation and Artificial Intelligence[J].Expert Systems with Applications, 2012, 39(9):7929-7937.

    [14] 夏樂天.馬爾可夫鏈預(yù)測方法及其在水文序列中的應(yīng)用研究[D].南京:河海大學(xué),2005:11-12.

    [15] 夏樂天,朱元甡.馬爾可夫鏈預(yù)測方法的統(tǒng)計(jì)試驗(yàn)研究[J].水利學(xué)報(bào),2007,38(增刊1):372-378.

    【責(zé)任編輯 張 帥】

    猜你喜歡
    原州區(qū)
    淺談原州區(qū)草畜產(chǎn)業(yè)發(fā)展情況
    固原市原州區(qū)水土保持生態(tài)建設(shè)成效及經(jīng)驗(yàn)啟示
    固原市原州區(qū)河川鄉(xiāng)干旱半干旱區(qū)油用牡丹栽培技術(shù)
    固原市原州區(qū)“十三五”林業(yè)特色優(yōu)勢產(chǎn)業(yè)發(fā)展探討
    加快推進(jìn)美麗原州建設(shè)的主要措施
    固原市原州區(qū)馬鈴薯種薯產(chǎn)業(yè)現(xiàn)狀及發(fā)展對策
    固原市原州區(qū)林木種苗產(chǎn)業(yè)發(fā)展現(xiàn)狀及對策
    原州區(qū)“十二五”時(shí)期經(jīng)濟(jì)社會發(fā)展取得的成就
    固原市原州區(qū)“十三五”林業(yè)生態(tài)建設(shè)發(fā)展探討
    固原市原州區(qū)中河鄉(xiāng)“20160731”暴雨洪水分析
    国产91av在线免费观看| 男人和女人高潮做爰伦理| av不卡在线播放| 亚洲精品色激情综合| 校园人妻丝袜中文字幕| 丰满迷人的少妇在线观看| 一区二区av电影网| 国产乱人偷精品视频| 多毛熟女@视频| 欧美区成人在线视频| 日本欧美视频一区| 国内精品宾馆在线| 成人亚洲精品一区在线观看| 十八禁高潮呻吟视频 | 免费人妻精品一区二区三区视频| 久久人人爽人人片av| 成人18禁高潮啪啪吃奶动态图 | 国产亚洲午夜精品一区二区久久| 亚洲精品日韩在线中文字幕| 欧美最新免费一区二区三区| 最新的欧美精品一区二区| 亚洲精品亚洲一区二区| 十八禁网站网址无遮挡 | a级毛色黄片| 亚洲色图综合在线观看| 国产日韩欧美视频二区| 亚洲精品456在线播放app| 亚洲美女搞黄在线观看| 欧美精品高潮呻吟av久久| 色94色欧美一区二区| 水蜜桃什么品种好| 新久久久久国产一级毛片| 两个人免费观看高清视频 | 成人无遮挡网站| 欧美三级亚洲精品| 免费观看的影片在线观看| a级片在线免费高清观看视频| 三级经典国产精品| 如日韩欧美国产精品一区二区三区 | 热99国产精品久久久久久7| 久久久久久久国产电影| 欧美日本中文国产一区发布| h日本视频在线播放| 日本与韩国留学比较| 久久久国产一区二区| 18禁在线无遮挡免费观看视频| 99热全是精品| 六月丁香七月| 免费人成在线观看视频色| 婷婷色综合大香蕉| 国产高清有码在线观看视频| 两个人的视频大全免费| 国产一区二区在线观看av| 精品久久久久久久久av| 99九九在线精品视频 | 蜜臀久久99精品久久宅男| 亚洲精品一区蜜桃| 亚洲美女搞黄在线观看| 亚洲情色 制服丝袜| 国精品久久久久久国模美| 亚洲成人手机| 亚洲四区av| 麻豆成人午夜福利视频| 男女啪啪激烈高潮av片| 尾随美女入室| 久久这里有精品视频免费| 性高湖久久久久久久久免费观看| 欧美日韩在线观看h| 日本猛色少妇xxxxx猛交久久| 国模一区二区三区四区视频| 国产淫片久久久久久久久| 成人国产麻豆网| 夜夜看夜夜爽夜夜摸| 色视频www国产| 亚洲美女搞黄在线观看| 久久久久久久久久久免费av| 亚洲欧洲国产日韩| 日日摸夜夜添夜夜爱| 国产精品嫩草影院av在线观看| 午夜激情久久久久久久| 国产高清三级在线| 啦啦啦啦在线视频资源| 如日韩欧美国产精品一区二区三区 | 男人狂女人下面高潮的视频| 亚洲国产av新网站| 久久人人爽人人爽人人片va| 亚洲欧美清纯卡通| 又黄又爽又刺激的免费视频.| 黄色一级大片看看| 99热国产这里只有精品6| 国产午夜精品久久久久久一区二区三区| 日本91视频免费播放| 美女xxoo啪啪120秒动态图| 中文乱码字字幕精品一区二区三区| 亚洲国产精品国产精品| 黄色日韩在线| 自拍欧美九色日韩亚洲蝌蚪91 | 国产男女超爽视频在线观看| 在线看a的网站| 80岁老熟妇乱子伦牲交| 精品国产国语对白av| 偷拍熟女少妇极品色| 国产视频内射| 国产欧美亚洲国产| 精品一区二区免费观看| 高清欧美精品videossex| 极品少妇高潮喷水抽搐| 人人妻人人澡人人爽人人夜夜| 亚洲欧美一区二区三区黑人 | 欧美高清成人免费视频www| 日本91视频免费播放| 黑人巨大精品欧美一区二区蜜桃 | 中国国产av一级| 哪个播放器可以免费观看大片| 欧美精品一区二区大全| 精品久久久久久久久av| a 毛片基地| 最近中文字幕高清免费大全6| 免费看日本二区| 欧美激情极品国产一区二区三区 | 国产精品一区二区性色av| 如日韩欧美国产精品一区二区三区 | 又粗又硬又长又爽又黄的视频| 一区二区三区免费毛片| 久久国产精品男人的天堂亚洲 | 亚洲av二区三区四区| 日韩制服骚丝袜av| 色婷婷久久久亚洲欧美| 狂野欧美白嫩少妇大欣赏| 亚洲欧美日韩卡通动漫| 久久久久久久久大av| 一区二区三区四区激情视频| 99热这里只有是精品50| 男女边摸边吃奶| 国产精品一区二区在线不卡| 国产爽快片一区二区三区| 亚洲中文av在线| 日产精品乱码卡一卡2卡三| 日本色播在线视频| 少妇丰满av| 日韩av在线免费看完整版不卡| 黑丝袜美女国产一区| 日韩制服骚丝袜av| 国产黄色视频一区二区在线观看| 国产成人freesex在线| 日本色播在线视频| 久久国内精品自在自线图片| 秋霞伦理黄片| 男女边摸边吃奶| 亚洲天堂av无毛| 亚洲中文av在线| 久久久久国产网址| 国产极品粉嫩免费观看在线 | 丰满迷人的少妇在线观看| 少妇猛男粗大的猛烈进出视频| 观看免费一级毛片| 精品熟女少妇av免费看| av网站免费在线观看视频| 亚洲欧美日韩东京热| 久久精品国产亚洲av涩爱| 日韩在线高清观看一区二区三区| 日韩制服骚丝袜av| 亚洲精品第二区| 久久av网站| 日韩中字成人| 国产精品久久久久成人av| 午夜影院在线不卡| a级毛片在线看网站| 91久久精品电影网| 亚洲精品久久久久久婷婷小说| 只有这里有精品99| av在线播放精品| 乱人伦中国视频| 国产精品.久久久| 少妇人妻 视频| 99久久精品国产国产毛片| 一级爰片在线观看| 中国美白少妇内射xxxbb| 少妇被粗大的猛进出69影院 | 久久久久久伊人网av| 熟女人妻精品中文字幕| 老司机影院毛片| 少妇的逼好多水| 久久精品国产亚洲网站| 国产精品一区二区在线观看99| 69精品国产乱码久久久| h日本视频在线播放| 亚洲无线观看免费| 国产白丝娇喘喷水9色精品| 99热网站在线观看| 国产精品一二三区在线看| 久久女婷五月综合色啪小说| 久久人人爽av亚洲精品天堂| 大又大粗又爽又黄少妇毛片口| 国产黄频视频在线观看| 六月丁香七月| av在线app专区| 黄色欧美视频在线观看| 不卡视频在线观看欧美| h日本视频在线播放| 各种免费的搞黄视频| 97精品久久久久久久久久精品| 蜜桃在线观看..| 亚洲欧美成人综合另类久久久| 热99国产精品久久久久久7| 欧美+日韩+精品| 美女国产视频在线观看| 精品人妻熟女av久视频| 99久久精品一区二区三区| 中文字幕免费在线视频6| 内地一区二区视频在线| 欧美97在线视频| 免费人妻精品一区二区三区视频| 夜夜骑夜夜射夜夜干| 人人妻人人看人人澡| 精品一区二区免费观看| 嘟嘟电影网在线观看| 国产精品国产av在线观看| 伊人久久精品亚洲午夜| 国产熟女午夜一区二区三区 | 亚洲av成人精品一区久久| 免费高清在线观看视频在线观看| 热re99久久国产66热| 亚洲av成人精品一区久久| 久久国产乱子免费精品| 51国产日韩欧美| 久久毛片免费看一区二区三区| 日本黄色日本黄色录像| 欧美亚洲 丝袜 人妻 在线| 日韩大片免费观看网站| 97精品久久久久久久久久精品| 夫妻性生交免费视频一级片| 视频中文字幕在线观看| 在线观看免费视频网站a站| 成人国产麻豆网| 精品一区二区三区视频在线| 国产av一区二区精品久久| 秋霞伦理黄片| 免费看不卡的av| 中文欧美无线码| 国产视频首页在线观看| 中文字幕久久专区| 国产 一区精品| 国产精品熟女久久久久浪| 久久av网站| 在线观看三级黄色| 亚洲av在线观看美女高潮| 亚洲av欧美aⅴ国产| 国产一区二区在线观看日韩| 日本vs欧美在线观看视频 | 精品亚洲成国产av| 寂寞人妻少妇视频99o| 啦啦啦中文免费视频观看日本| 国产日韩欧美在线精品| 美女国产视频在线观看| 内地一区二区视频在线| 日本午夜av视频| 女性被躁到高潮视频| 久久久久国产网址| 视频区图区小说| 国语对白做爰xxxⅹ性视频网站| 男人爽女人下面视频在线观看| 18禁在线播放成人免费| 性色av一级| 色网站视频免费| 色哟哟·www| 伊人久久国产一区二区| 精华霜和精华液先用哪个| 欧美xxxx性猛交bbbb| 亚洲成人av在线免费| 少妇的逼好多水| 啦啦啦啦在线视频资源| 国产成人一区二区在线| 天堂中文最新版在线下载| av又黄又爽大尺度在线免费看| 亚洲欧美一区二区三区国产| 欧美日韩一区二区视频在线观看视频在线| 亚洲,一卡二卡三卡| 99久久精品一区二区三区| 我的女老师完整版在线观看| 久久久久久人妻| 伦精品一区二区三区| 一区在线观看完整版| 少妇猛男粗大的猛烈进出视频| 最近手机中文字幕大全| 久久久久久久久久久丰满| 国产精品人妻久久久久久| 观看av在线不卡| 丰满少妇做爰视频| 内射极品少妇av片p| 久久精品国产自在天天线| 三级国产精品片| 边亲边吃奶的免费视频| 18禁在线无遮挡免费观看视频| 一个人看视频在线观看www免费| 欧美精品国产亚洲| 日本黄色片子视频| 人妻一区二区av| 亚洲精品自拍成人| 成人漫画全彩无遮挡| 亚洲国产毛片av蜜桃av| 97超视频在线观看视频| 一个人免费看片子| 亚洲精品日韩av片在线观看| 亚洲精品456在线播放app| 婷婷色麻豆天堂久久| 热99国产精品久久久久久7| 国产高清不卡午夜福利| 毛片一级片免费看久久久久| 一级爰片在线观看| 91久久精品电影网| 国产成人免费观看mmmm| 能在线免费看毛片的网站| 精品亚洲乱码少妇综合久久| 精品国产一区二区久久| 成人二区视频| 哪个播放器可以免费观看大片| 国产日韩欧美视频二区| 免费人妻精品一区二区三区视频| 久久99热6这里只有精品| 日韩一本色道免费dvd| 免费大片18禁| 在线观看一区二区三区激情| 午夜福利网站1000一区二区三区| 久久久久视频综合| 日本猛色少妇xxxxx猛交久久| 精品人妻熟女av久视频| 久久久久久久久久久久大奶| 午夜91福利影院| 99久久人妻综合| 中文字幕人妻丝袜制服| 亚洲av免费高清在线观看| 久久鲁丝午夜福利片| 在线观看美女被高潮喷水网站| 在线播放无遮挡| 国产亚洲91精品色在线| av不卡在线播放| av在线app专区| 天天操日日干夜夜撸| 七月丁香在线播放| 丰满饥渴人妻一区二区三| av网站免费在线观看视频| 在线 av 中文字幕| 22中文网久久字幕| 精品国产国语对白av| 国产黄片美女视频| 午夜老司机福利剧场| 国产成人精品福利久久| 中文字幕av电影在线播放| 一个人免费看片子| 全区人妻精品视频| 亚洲自偷自拍三级| 人妻少妇偷人精品九色| 亚洲国产精品999| 精品一区在线观看国产| a级片在线免费高清观看视频| 国产视频首页在线观看| 午夜福利,免费看| 99九九线精品视频在线观看视频| 在线精品无人区一区二区三| 久久久久久久精品精品| 国产在线视频一区二区| 少妇猛男粗大的猛烈进出视频| a级片在线免费高清观看视频| 亚洲av男天堂| 久久久久国产精品人妻一区二区| 国产伦精品一区二区三区四那| 国产精品不卡视频一区二区| 一级,二级,三级黄色视频| 国产淫片久久久久久久久| 国产精品.久久久| 91精品国产国语对白视频| 国产精品一区二区三区四区免费观看| 精品少妇内射三级| 欧美日韩一区二区视频在线观看视频在线| 在现免费观看毛片| 边亲边吃奶的免费视频| 国产精品一区www在线观看| .国产精品久久| 国产精品女同一区二区软件| 久久精品国产a三级三级三级| 国产永久视频网站| 日韩成人av中文字幕在线观看| 亚洲不卡免费看| 男女边摸边吃奶| av女优亚洲男人天堂| 亚洲第一区二区三区不卡| 男人和女人高潮做爰伦理| 国产探花极品一区二区| av卡一久久| 久久久国产精品麻豆| 久久精品久久精品一区二区三区| 三级经典国产精品| 一本大道久久a久久精品| 国产av精品麻豆| 精品少妇久久久久久888优播| 亚洲人成网站在线观看播放| 街头女战士在线观看网站| 毛片一级片免费看久久久久| 美女中出高潮动态图| 国产精品秋霞免费鲁丝片| 国产淫片久久久久久久久| 久久久久久久久久久久大奶| 99热这里只有是精品在线观看| 欧美日韩在线观看h| 香蕉精品网在线| 美女主播在线视频| 日日啪夜夜爽| 日本av免费视频播放| 亚洲四区av| 一级毛片黄色毛片免费观看视频| 一级毛片久久久久久久久女| 精品久久久噜噜| 国产真实伦视频高清在线观看| 王馨瑶露胸无遮挡在线观看| a级片在线免费高清观看视频| 亚洲欧洲国产日韩| 欧美精品高潮呻吟av久久| 久久久国产一区二区| 少妇的逼水好多| 亚洲情色 制服丝袜| 九草在线视频观看| 午夜福利影视在线免费观看| 人妻夜夜爽99麻豆av| 51国产日韩欧美| 中文欧美无线码| 热re99久久精品国产66热6| 交换朋友夫妻互换小说| 啦啦啦视频在线资源免费观看| 91成人精品电影| 只有这里有精品99| 女性被躁到高潮视频| 一个人看视频在线观看www免费| 一二三四中文在线观看免费高清| 久久亚洲国产成人精品v| 久久99热这里只频精品6学生| 蜜桃在线观看..| 久久人人爽人人爽人人片va| 高清在线视频一区二区三区| 少妇人妻 视频| 美女视频免费永久观看网站| 高清黄色对白视频在线免费看 | 在线 av 中文字幕| 涩涩av久久男人的天堂| 看免费成人av毛片| 成人国产av品久久久| 99热这里只有是精品50| 国产午夜精品一二区理论片| 最新中文字幕久久久久| 中文在线观看免费www的网站| 久久人人爽人人爽人人片va| 国产成人精品一,二区| 中国三级夫妇交换| 夜夜爽夜夜爽视频| 亚洲欧洲精品一区二区精品久久久 | 菩萨蛮人人尽说江南好唐韦庄| 热re99久久国产66热| 在线观看免费视频网站a站| 国产男人的电影天堂91| 精品亚洲成国产av| 黑人猛操日本美女一级片| 内射极品少妇av片p| 好男人视频免费观看在线| 我的老师免费观看完整版| 久久99热6这里只有精品| 人人妻人人澡人人看| 中文在线观看免费www的网站| 午夜久久久在线观看| 久久av网站| 亚洲婷婷狠狠爱综合网| 成人免费观看视频高清| 日日爽夜夜爽网站| 人人妻人人爽人人添夜夜欢视频 | 日本欧美国产在线视频| 高清av免费在线| 欧美日韩亚洲高清精品| 日本色播在线视频| 丝袜喷水一区| 亚洲av电影在线观看一区二区三区| 国产在线免费精品| 在线 av 中文字幕| 国产真实伦视频高清在线观看| 18禁动态无遮挡网站| 一区在线观看完整版| 热re99久久精品国产66热6| 国产女主播在线喷水免费视频网站| 精品少妇内射三级| 国产精品一区二区三区四区免费观看| 最近中文字幕高清免费大全6| 亚洲成人av在线免费| 国产成人aa在线观看| 观看av在线不卡| 亚洲国产毛片av蜜桃av| 国产精品女同一区二区软件| 九九在线视频观看精品| 久久精品国产亚洲网站| 精品一区二区免费观看| 国产又色又爽无遮挡免| av卡一久久| av网站免费在线观看视频| 欧美精品亚洲一区二区| 男的添女的下面高潮视频| 性色avwww在线观看| 成人免费观看视频高清| 天堂8中文在线网| 午夜福利,免费看| av国产精品久久久久影院| 欧美三级亚洲精品| 免费看av在线观看网站| 丰满人妻一区二区三区视频av| av视频免费观看在线观看| 免费大片18禁| 欧美xxⅹ黑人| 久热久热在线精品观看| 日韩中文字幕视频在线看片| 99国产精品免费福利视频| 午夜老司机福利剧场| freevideosex欧美| 新久久久久国产一级毛片| 亚洲av男天堂| 免费黄频网站在线观看国产| 亚洲欧洲国产日韩| 亚洲精品456在线播放app| 国产黄频视频在线观看| 啦啦啦啦在线视频资源| 久久久久精品久久久久真实原创| 午夜免费观看性视频| 卡戴珊不雅视频在线播放| 三级国产精品欧美在线观看| 女的被弄到高潮叫床怎么办| 精品久久久久久久久av| 国产白丝娇喘喷水9色精品| 乱系列少妇在线播放| 免费播放大片免费观看视频在线观看| 午夜免费观看性视频| 久久韩国三级中文字幕| 永久免费av网站大全| 国产一级毛片在线| 日韩精品有码人妻一区| 亚洲精品一区蜜桃| 亚洲av.av天堂| 精品亚洲成a人片在线观看| 精品国产乱码久久久久久小说| 在线亚洲精品国产二区图片欧美 | 人人妻人人添人人爽欧美一区卜| 一级毛片黄色毛片免费观看视频| 妹子高潮喷水视频| 又粗又硬又长又爽又黄的视频| 人妻制服诱惑在线中文字幕| 麻豆成人午夜福利视频| 少妇人妻 视频| 18+在线观看网站| 插阴视频在线观看视频| 边亲边吃奶的免费视频| 黑人巨大精品欧美一区二区蜜桃 | 看十八女毛片水多多多| 91精品伊人久久大香线蕉| 免费看光身美女| 国产成人免费无遮挡视频| 精品人妻熟女毛片av久久网站| 国产老妇伦熟女老妇高清| 国产欧美日韩综合在线一区二区 | 精品午夜福利在线看| 久久久久国产网址| 欧美xxxx性猛交bbbb| 国产成人免费观看mmmm| 久久精品国产亚洲av涩爱| 亚洲国产色片| 亚洲伊人久久精品综合| 日韩成人伦理影院| 日韩亚洲欧美综合| 精品久久久久久久久亚洲| av在线app专区| 一级片'在线观看视频| av视频免费观看在线观看| 91久久精品国产一区二区三区| 日韩在线高清观看一区二区三区| 狂野欧美激情性xxxx在线观看| 在线精品无人区一区二区三| 国产男女内射视频| 熟女av电影| 一本大道久久a久久精品| 一区在线观看完整版| av黄色大香蕉| 精品久久国产蜜桃| 国产成人freesex在线| 中文字幕精品免费在线观看视频 | 亚洲av电影在线观看一区二区三区| 国产精品人妻久久久久久| 18禁在线播放成人免费| 久久久久久久久久久免费av| 国产日韩欧美视频二区| 亚洲国产欧美在线一区| 日韩熟女老妇一区二区性免费视频| 99热全是精品| 大码成人一级视频| 成人亚洲精品一区在线观看| 丰满饥渴人妻一区二区三| 2018国产大陆天天弄谢| 久久久久久久久久久免费av| 国产亚洲一区二区精品| 国产欧美另类精品又又久久亚洲欧美| 国产在线免费精品| 人妻系列 视频| 又粗又硬又长又爽又黄的视频| 亚洲久久久国产精品| 天美传媒精品一区二区| 国产高清不卡午夜福利| 一区二区三区免费毛片| 欧美xxⅹ黑人| 狂野欧美激情性bbbbbb| 黄片无遮挡物在线观看| 久久久久久人妻| 啦啦啦在线观看免费高清www| 亚洲性久久影院|