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

    新型聯(lián)合模型在多站點(diǎn)同步降水預(yù)報(bào)中的應(yīng)用

    2018-11-28 09:40:24吳金冉李維德孔德萌
    節(jié)水灌溉 2018年11期
    關(guān)鍵詞:果洛月度協(xié)整

    吳金冉,李維德,孔德萌

    (蘭州大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,蘭州 730000)

    淡水資源危機(jī)現(xiàn)已成為嚴(yán)重制約西北半旱區(qū)可持續(xù)發(fā)展的首要因素。而其中,半旱區(qū)降水作為該地區(qū)主要淡水來源,其的預(yù)測研究分析對半旱區(qū)的居民農(nóng)牧活動及野生物種生存等各方面具有重要指導(dǎo)意義。而本文以生態(tài)穩(wěn)定性較為脆弱的青海省海東地區(qū)和果洛州的月度降水為研究對象,通過構(gòu)造具有高精度的降水預(yù)報(bào)模型以對該地區(qū)的經(jīng)濟(jì)農(nóng)業(yè)生活提供有效的指引。

    一般而言,依據(jù)不同的模型構(gòu)造原理,降水預(yù)報(bào)模型通??煞譃閮纱箢愋停簲?shù)據(jù)驅(qū)動模型(Data Driven Models)和物理過程模型(Physical Process Models)[1,2]。第一大類方法主要依據(jù)所獲得的各種氣象數(shù)據(jù)集,通過“黑箱”建立輸入和輸處變量間的數(shù)量關(guān)系,這里需要指明,該方法僅僅依賴一條序列信息也可以獲得較為突出預(yù)報(bào)效果;而第二大類方法,通常要求較高的空氣動力學(xué)系統(tǒng)基礎(chǔ),并依賴于大量的變量信息(濕度,日照,溫度等),通過極為復(fù)雜數(shù)學(xué)公式來有效的表示降水的物理過程,從而獲得較為準(zhǔn)確的降水預(yù)報(bào)值[3,4]。而近些年來,以機(jī)器學(xué)習(xí)為代表的數(shù)據(jù)驅(qū)動模型,憑借其出色的獲取歷史經(jīng)驗(yàn)數(shù)據(jù)集特征的能力,被廣泛運(yùn)用于水文數(shù)據(jù)預(yù)報(bào)研究中,其中,又以支持向量機(jī)模型(Support Vector Machine, SVM)和人工神經(jīng)網(wǎng)絡(luò)模型(artificial neural network, ANN)最為常見[5-7]。雖然該類建模方法可以通過觀測點(diǎn)數(shù)據(jù)集來建立非線性系統(tǒng),但大量該方法研究僅僅局限于單點(diǎn)(觀測點(diǎn))降水預(yù)報(bào)上,而并未將空間相關(guān)站(多觀測)點(diǎn)信息運(yùn)用于模型構(gòu)造中。這里需要提及,機(jī)器學(xué)習(xí)學(xué)家Vapnik于1995年創(chuàng)立的支持向量機(jī)模型以其出色的非線性擬合能力和魯棒性,該模型現(xiàn)已廣泛的運(yùn)用于河流徑流量、臺風(fēng)降水、月度降水和地下水井深等水文變量預(yù)報(bào)中來[8-11];統(tǒng)計(jì)學(xué)家Engle和Granger為了解決多條時間序列的非平穩(wěn)性所導(dǎo)致的“偽回歸”現(xiàn)象的發(fā)生,提出了協(xié)整分析(Co-integration Analysis)[12];而當(dāng)前,該方法已初步被應(yīng)用于降水預(yù)報(bào)和空氣顆粒預(yù)報(bào)等環(huán)境工程研究中[13,14]。而本文以西北半旱區(qū)海東地區(qū)和果洛州月度降水為研究對象,將上述兩種方法進(jìn)行結(jié)構(gòu)以構(gòu)建聯(lián)合降水預(yù)報(bào)模型。具體的來看,將協(xié)整分析運(yùn)用于降水預(yù)報(bào)研究中,打破了傳統(tǒng)單一觀測點(diǎn)的降水預(yù)報(bào)建模方法,且以經(jīng)典的支持向量機(jī)模型作為主算法來獲得輸入變量(input variables)和目標(biāo)(targets)之間的“黑箱”(非線性系統(tǒng))來對多觀測點(diǎn)降水進(jìn)行同步降水預(yù)報(bào)。

    本文主要介紹降水預(yù)報(bào)的意義、常用的建模技術(shù)以及文章所提出的方法技術(shù)的發(fā)展;本文采用的預(yù)報(bào)技術(shù)詳情,主要包括聯(lián)合預(yù)報(bào)模型所運(yùn)用到的協(xié)整理論、支持向量機(jī)原理和聯(lián)合模型構(gòu)造的過程;主要包含驗(yàn)證模型所運(yùn)用的研究對象簡介、數(shù)據(jù)集、模型仿真結(jié)果以及相關(guān)對比模型的比較分析結(jié)果;總結(jié)所構(gòu)建的基于協(xié)整關(guān)系—支持向量機(jī)模型(CI-SVM)在半干旱區(qū)降水預(yù)報(bào)研究中的突出的效果,以及該預(yù)報(bào)模型的結(jié)果在海東地區(qū)和果洛州的降水資源管理規(guī)劃中的實(shí)際意義。

    1 預(yù)報(bào)技術(shù)

    1.1 支持向量機(jī)

    為了獲得高效的機(jī)器學(xué)習(xí)模型來解決分類和回歸問題,1997年Vapnik教授改進(jìn)了支持向量機(jī)模型[14]。該方法主要是尋求模型的經(jīng)驗(yàn)誤差與模型復(fù)雜度之間的最佳權(quán)衡關(guān)系,而此關(guān)系是將支持向量機(jī)回歸(SVR)的表達(dá)函數(shù)f(·)約束至分類的超平面所得到的。這里可以認(rèn)為,f(·)只決定于減少的訓(xùn)練集的支持向量(support vectors),但支持向量對于優(yōu)化問題的約束空間。具體的來說,存在m個樣本(sample)數(shù)據(jù)集(xi,yi)∈RN×R,其中N為輸入變量的維度,可以得到,支持向量回歸優(yōu)化問題的數(shù)學(xué)表達(dá)式為:

    maxW(α,α*)=

    (1)

    (2)

    表達(dá)式(2)中,C為模型的負(fù)責(zé)度懲罰參數(shù),α和α*是約束的二變量。需要說明的是,支持向量機(jī)回歸引入了核函數(shù)(kernel function)k(xi,x)來將非線性回歸問題轉(zhuǎn)化為線性回歸問題來解決,進(jìn)一步,結(jié)合公式(1)可得到模型的最終模型為:

    (3)

    1.2 協(xié)整理論

    由于時間序列變量間的回歸容易導(dǎo)致“偽回歸”的發(fā)生;因此,為了描述非平穩(wěn)時間序列間的長期均衡關(guān)系的而提出了協(xié)整的概念。若存在某非平穩(wěn)時間序列經(jīng)過d次差分后為平穩(wěn)序列,則稱該序列為d階單整,也記為I(d)[12]。若存在Xt和Yt為同階單整序列,且其的線性組合ut=Yt-βXt為平穩(wěn)序列,即uy~I(xiàn)(d),則稱Xt和Yt為協(xié)整關(guān)系。這里需要重點(diǎn)指出,為了判斷兩條序列間存在顯著的因果關(guān)系,1987年Engle和Granger提出了E-G兩步檢驗(yàn)法。該主要想法是,第一步,運(yùn)用最小二乘法進(jìn)行用Xt對Yt做回歸模型,得到殘差εt;第二步,檢驗(yàn)殘差εt是否為平穩(wěn)序列,若是則認(rèn)為二變量間為協(xié)整關(guān)系,若否則認(rèn)為不存在協(xié)整關(guān)系。

    1.3 聯(lián)合模型的構(gòu)建

    本文所建立的聯(lián)合同步降水預(yù)報(bào)模型----協(xié)整關(guān)系的支持向量機(jī)模型(CI-SVM)的主要想法是,利用協(xié)整理論檢驗(yàn)多個站點(diǎn)間的空間相關(guān)性,并以此來增加在構(gòu)建降水預(yù)報(bào)模型中可運(yùn)用到的相關(guān)信息,然后,采用支持向量機(jī)模型來建立輸入和目標(biāo)變量間的非線性系統(tǒng),進(jìn)而得到多空間同步聯(lián)合降水預(yù)測結(jié)果。具體的模型流程如下:

    Step 1:輸入兩地區(qū)降水時間序列信息,記為{Xt}和{Yt}。

    Step 2: 對{Xt}和{Yt}的是否存在偽回歸進(jìn)行協(xié)整檢驗(yàn)。

    Step 3: 用協(xié)整檢驗(yàn)結(jié)果確定兩條序列間的輸入變量和目標(biāo)之間的映射關(guān)系。

    (4)

    式中:p和q分別為{Xt}和{Yt}的滯后期數(shù);f為支持向量機(jī)所獲得的“黑箱”。

    Step 4: 通過所建立的非線性映射系統(tǒng)進(jìn)行多空間聯(lián)合降水預(yù)報(bào)。

    2 仿真及其結(jié)果

    2.1 研究對象

    為了驗(yàn)證所提出的聯(lián)合降水預(yù)報(bào)模型的效果,本文選取了位于西北半干旱區(qū)的海東地區(qū)和果洛州的月度降水量為預(yù)報(bào)對象。海東地區(qū)是青海省重要的農(nóng)牧業(yè)生產(chǎn)基地和鄉(xiāng)鎮(zhèn)企業(yè)發(fā)達(dá)區(qū)域之一,而果洛州位于青藏高原的腹地,且二地區(qū)皆具有高寒、缺氧、低溫、干旱和日照較長等特征,二者的年平均降水量取值分別為319~531和400~760 mm。

    2.2 數(shù)據(jù)集的選擇

    本研究對象為青海省海東地區(qū)和果洛州的月度降水h,而本次試驗(yàn)數(shù)據(jù)集主要源于青海省統(tǒng)計(jì)年鑒(http:∥www.qhtjj.gov.cn/tjData/qhtjnj/)。其具體月度降水量趨勢如圖1所示。

    海東地區(qū)和果洛州的月度降水量如圖1 所示,可以得出,二地區(qū)的月度降水量存在明顯的周期性變化規(guī)律;具體的看,夏季(6-8月)的降水量遠(yuǎn)遠(yuǎn)高于其他三季,而冬季的降水量最低;并且,對比兩折線趨勢可發(fā)現(xiàn),海東地區(qū)和果洛州具有較高的同步性,而該特征將被運(yùn)用于協(xié)整關(guān)系的建模分析中。通過SPSS軟件對兩地降水量進(jìn)行描述統(tǒng)計(jì)分析得到表1。

    表1 海東和果洛月度降水量描述分析 mm

    圖1 海東和果洛月度降水量

    這里需要說明,為了驗(yàn)證所提出聯(lián)合預(yù)報(bào)模型的高效,本文所提到的模型所采用的訓(xùn)練數(shù)據(jù)集和測試數(shù)據(jù)集分別為前96個月降水序列數(shù)據(jù)和后12個月度降水序列數(shù)據(jù)。

    2.3 協(xié)整檢驗(yàn)及分析

    海東和果洛的月度降水量平穩(wěn)性和協(xié)整檢驗(yàn)結(jié)果見表2。

    表2 海東和果洛的協(xié)整檢驗(yàn)

    通過表2可得,序列平穩(wěn)性ADF單位根檢驗(yàn)的原假設(shè)為存在單位根,通過計(jì)算得到海東和果洛的月度降水量序列的p值分別為0和0.083,在置信水平10%下,拒絕原假設(shè),則認(rèn)為該兩條降水序列皆為0階單整的。然后,通過建立回歸模型得出對應(yīng)的殘差序列,在置信水平1%下,該序列的平穩(wěn)性檢驗(yàn)的p值為0,拒絕原假設(shè),則可認(rèn)為其是平穩(wěn)的,進(jìn)一步得到海東地區(qū)和果洛州的月降水量間存在顯著的協(xié)整關(guān)系。經(jīng)過多次試驗(yàn)嘗試,本次仿真以海東地區(qū)和果洛州前36月降水量為預(yù)報(bào)模型的輸入變量,兩地下一月的降水量為預(yù)測模型的目標(biāo)變量進(jìn)行滾動預(yù)測。

    2.4 預(yù)報(bào)結(jié)果

    在Matlab2018a版本下建立支持向量機(jī)回歸模型,結(jié)果見表3。

    表3 CI-SVM 模型對海東和果洛月度降水測試指標(biāo)

    表3為提出的CI-SVM模型分別對海東地區(qū)和果洛州月度降水量的預(yù)報(bào)結(jié)果。① 從絕對水平上,二地區(qū)的所有平均絕對誤差MAE皆小于12.5;② 從相對水平上,標(biāo)準(zhǔn)化平均方差誤NMSE全都小于0.13;③從相關(guān)水平上,納什系數(shù)NS都大于0.86;可認(rèn)為所提出的聯(lián)合模型CI-SVM對海東地區(qū)及果洛州的月降水量具有較好的預(yù)報(bào)效果。具體來看,海東地區(qū)和果洛州的月降水預(yù)報(bào)的三指標(biāo)分別為7.97、0.1、0.89、12.37、0.12、0.87。其對海東地區(qū)和果洛州的月度降水預(yù)報(bào)的結(jié)果具體見圖2。

    圖2 CI-SVM模型的海東和果洛聯(lián)合預(yù)報(bào)結(jié)果

    2.5 模型比較

    在此構(gòu)建兩種模型(即SVM和CI-NN)作為參照模型,和CI-SVM模型就海東地區(qū)和果洛州月降水量的預(yù)測精確度進(jìn)行對比分析。SVM即持向量回歸模型,不考慮序列相關(guān)性,除去CI-SVM的協(xié)整關(guān)系;CI-NN即BP神經(jīng)網(wǎng)絡(luò)模型,建立在協(xié)整關(guān)系基礎(chǔ)上,用ANN替換CI-SVM模型中的SVM。具體對比分析見表4。

    參照表4,對比CI-SVM、CI-NN、SVM 3種模型預(yù)測結(jié)果可得,CI-SVM模型的預(yù)測精確度最好,該模型中平均MAE10.17、平均NMSE0.11、平均NS系數(shù)0.88,此三項(xiàng)指標(biāo)優(yōu)于CI-NN模型的相應(yīng)指標(biāo)數(shù)據(jù)10.23、0.15、0.85和SVM模型的相應(yīng)指標(biāo)13.51、0.17、0.83。

    表4 模型指標(biāo)對比分析

    由表4得出,CI-SVM模型優(yōu)于CI-NN模型,由此可得,在聯(lián)合模型中SVM貢獻(xiàn)值優(yōu)于ANN,即SVM作為主體算法所得的效果更好。從數(shù)據(jù)分析看來,CI-SVM指標(biāo)優(yōu)于SVM,由此得出,多觀測點(diǎn)同步預(yù)報(bào)的精度優(yōu)于單觀測點(diǎn)預(yù)報(bào)的精度,即由協(xié)整理論引入的多觀測點(diǎn)信息提高了模型的精準(zhǔn)性。綜上所述,相較于其他模型,經(jīng)由CI-SVM模型所得出的預(yù)測結(jié)果是非常精確、科學(xué)且合理的。

    3 結(jié) 語

    本文所構(gòu)建的二空間聯(lián)合同步降水預(yù)報(bào)模型CI-SVM在海東地區(qū)和果洛州月度降水?dāng)?shù)據(jù)集的驗(yàn)證中獲得較為突出的效果,其的平均MAE、平均NMSE和平均NS系數(shù)分別為10.17、0.11和0.88。進(jìn)一步可以得到,所提出的CI-SVM模型利用協(xié)整關(guān)系將傳統(tǒng)的單觀測點(diǎn)降水預(yù)報(bào)拓展到了多空間同步聯(lián)合預(yù)測,通過提高建模中所使用的相關(guān)空間信息,進(jìn)而取得了更優(yōu)的預(yù)報(bào)精度。并且,通過其他對比模型,該聯(lián)合模型驗(yàn)證了機(jī)器學(xué)習(xí)中的支持向量回歸模型可以提高月度降水預(yù)報(bào)的效果。最終,可以將該聯(lián)合模型的預(yù)報(bào)結(jié)果運(yùn)用于半旱區(qū)的生態(tài)監(jiān)管和規(guī)劃中,且運(yùn)用計(jì)算出來的結(jié)果給農(nóng)牧民的產(chǎn)品種植以及家畜養(yǎng)殖提供些指導(dǎo)意見,進(jìn)而,促進(jìn)海東地區(qū)和果洛州的可持續(xù)發(fā)展。

    猜你喜歡
    果洛月度協(xié)整
    故鄉(xiāng)果洛(外一首)
    果洛州草原建設(shè)綜合措施初探
    外商直接投資對我國進(jìn)出口貿(mào)易影響的協(xié)整分析
    智富時代(2019年2期)2019-04-18 07:44:42
    河南金融發(fā)展和城鄉(xiāng)居民收入差距的協(xié)整分析
    智富時代(2018年3期)2018-06-11 16:10:44
    果洛瑪沁的九天九夜——青海果洛瑪沁“8 ·30”山體滑坡?lián)岆U(xiǎn)紀(jì)實(shí)
    中國公路(2017年19期)2018-01-23 03:06:32
    果洛州瑪沁縣大武鎮(zhèn)商品豬旋毛蟲感染情況調(diào)查報(bào)告
    月度聚焦
    中國記者(2015年8期)2015-05-09 08:30:35
    月度聚焦
    中國記者(2014年4期)2014-05-14 06:04:39
    月度聚焦
    中國記者(2014年9期)2014-03-01 01:44:22
    月度聚焦
    中國記者(2014年6期)2014-03-01 01:39:52
    如日韩欧美国产精品一区二区三区| 国产av精品麻豆| 国产91精品成人一区二区三区 | 久久久久久久久免费视频了| 大码成人一级视频| 精品久久久久久电影网| 亚洲人成电影免费在线| a级毛片在线看网站| 免费女性裸体啪啪无遮挡网站| 欧美日韩亚洲综合一区二区三区_| 美女高潮到喷水免费观看| 亚洲第一青青草原| 99国产极品粉嫩在线观看| 精品福利观看| 欧美日韩精品网址| 午夜福利影视在线免费观看| 欧美另类一区| 老司机福利观看| 亚洲综合色网址| 黄色片一级片一级黄色片| 免费看十八禁软件| 又大又爽又粗| 亚洲av国产av综合av卡| 99国产精品一区二区蜜桃av | 黄色a级毛片大全视频| www.熟女人妻精品国产| 男女国产视频网站| 9色porny在线观看| 精品国内亚洲2022精品成人 | 成人国产一区最新在线观看| 人人妻,人人澡人人爽秒播| 自线自在国产av| 欧美国产精品一级二级三级| 一进一出抽搐动态| 亚洲人成电影免费在线| 一进一出抽搐动态| 日本一区二区免费在线视频| 亚洲av男天堂| 婷婷丁香在线五月| 亚洲欧洲日产国产| 久久久久国产一级毛片高清牌| 亚洲人成电影免费在线| 国产伦人伦偷精品视频| 国产伦理片在线播放av一区| 久久久国产欧美日韩av| 国产亚洲av片在线观看秒播厂| 精品亚洲成国产av| 在线观看人妻少妇| 国产亚洲精品第一综合不卡| 国产日韩欧美在线精品| 久久精品aⅴ一区二区三区四区| 久热这里只有精品99| 在线永久观看黄色视频| 一区福利在线观看| 亚洲视频免费观看视频| 精品一品国产午夜福利视频| 国产精品亚洲av一区麻豆| 成人av一区二区三区在线看 | 中文欧美无线码| 一区二区三区激情视频| 国产日韩欧美视频二区| 男女免费视频国产| 午夜福利一区二区在线看| 丁香六月天网| 国产精品久久久av美女十八| 午夜福利乱码中文字幕| 久久精品熟女亚洲av麻豆精品| 巨乳人妻的诱惑在线观看| 午夜影院在线不卡| 成年动漫av网址| 一区二区日韩欧美中文字幕| www.999成人在线观看| 成人免费观看视频高清| 亚洲一码二码三码区别大吗| 狠狠狠狠99中文字幕| 永久免费av网站大全| 免费看十八禁软件| 亚洲一卡2卡3卡4卡5卡精品中文| 国产欧美日韩综合在线一区二区| 99热网站在线观看| 2018国产大陆天天弄谢| 青春草亚洲视频在线观看| 一级毛片女人18水好多| 91精品国产国语对白视频| 欧美日韩国产mv在线观看视频| 又黄又粗又硬又大视频| 精品福利观看| 丰满饥渴人妻一区二区三| 欧美日韩黄片免| 日韩欧美国产一区二区入口| 亚洲第一av免费看| 久久狼人影院| 日本精品一区二区三区蜜桃| 亚洲视频免费观看视频| 最近中文字幕2019免费版| 久久99一区二区三区| 丝袜人妻中文字幕| 午夜精品久久久久久毛片777| 少妇人妻久久综合中文| 久久久久国产一级毛片高清牌| 黑人猛操日本美女一级片| 丁香六月欧美| 亚洲中文日韩欧美视频| 亚洲中文日韩欧美视频| 亚洲午夜精品一区,二区,三区| 亚洲熟女毛片儿| 亚洲国产欧美日韩在线播放| 免费黄频网站在线观看国产| 人妻 亚洲 视频| 日韩三级视频一区二区三区| 90打野战视频偷拍视频| 久久中文字幕一级| 国产成人精品无人区| 欧美精品一区二区免费开放| 国产成人av教育| 欧美另类亚洲清纯唯美| 人人妻人人添人人爽欧美一区卜| 色视频在线一区二区三区| 真人做人爱边吃奶动态| 一边摸一边抽搐一进一出视频| 国产亚洲av高清不卡| 在线观看www视频免费| 12—13女人毛片做爰片一| 嫩草影视91久久| 日韩精品免费视频一区二区三区| 99re6热这里在线精品视频| 国内毛片毛片毛片毛片毛片| 亚洲一码二码三码区别大吗| 色精品久久人妻99蜜桃| 国产熟女午夜一区二区三区| e午夜精品久久久久久久| 国产男女内射视频| 欧美精品av麻豆av| 亚洲精品av麻豆狂野| 99久久综合免费| 黄片大片在线免费观看| 三级毛片av免费| 各种免费的搞黄视频| 高清av免费在线| 国产99久久九九免费精品| 91精品国产国语对白视频| netflix在线观看网站| 色播在线永久视频| 岛国在线观看网站| 欧美日韩黄片免| 99久久人妻综合| 国产精品久久久人人做人人爽| www.999成人在线观看| 免费观看人在逋| 欧美在线一区亚洲| 国产熟女午夜一区二区三区| 久久久久久亚洲精品国产蜜桃av| 欧美大码av| 精品久久久精品久久久| 久久久精品区二区三区| 嫩草影视91久久| 久久久精品免费免费高清| 国产男女超爽视频在线观看| 视频在线观看一区二区三区| 18禁裸乳无遮挡动漫免费视频| 天堂8中文在线网| 中文字幕精品免费在线观看视频| 亚洲欧美日韩另类电影网站| 99久久人妻综合| 久久久久精品人妻al黑| 成年美女黄网站色视频大全免费| 国产色视频综合| 国产亚洲欧美精品永久| 99精品欧美一区二区三区四区| 亚洲精品美女久久久久99蜜臀| 母亲3免费完整高清在线观看| 亚洲欧美成人综合另类久久久| 夫妻午夜视频| 午夜免费鲁丝| 老熟妇乱子伦视频在线观看 | 欧美+亚洲+日韩+国产| h视频一区二区三区| 五月开心婷婷网| 日韩人妻精品一区2区三区| 国产成人一区二区三区免费视频网站| 国产日韩一区二区三区精品不卡| 99久久综合免费| 麻豆乱淫一区二区| 国产又爽黄色视频| 久久久国产精品麻豆| 在线观看免费日韩欧美大片| 91麻豆av在线| 王馨瑶露胸无遮挡在线观看| 动漫黄色视频在线观看| 日韩视频在线欧美| 日本黄色日本黄色录像| 国产福利在线免费观看视频| 国产一区二区在线观看av| 久久久久久久精品精品| 亚洲午夜精品一区,二区,三区| 亚洲av成人不卡在线观看播放网 | 日本猛色少妇xxxxx猛交久久| 亚洲av成人一区二区三| 少妇粗大呻吟视频| 他把我摸到了高潮在线观看 | 国产有黄有色有爽视频| 别揉我奶头~嗯~啊~动态视频 | 精品视频人人做人人爽| 国产有黄有色有爽视频| 热99re8久久精品国产| 久久人人爽av亚洲精品天堂| 曰老女人黄片| 丁香六月天网| 人人妻人人澡人人爽人人夜夜| 免费观看a级毛片全部| 中国国产av一级| 人妻人人澡人人爽人人| 美女视频免费永久观看网站| 国产日韩欧美在线精品| 成人免费观看视频高清| 亚洲一区二区三区欧美精品| 国产麻豆69| 久久热在线av| 9热在线视频观看99| 99国产极品粉嫩在线观看| 久久亚洲国产成人精品v| 叶爱在线成人免费视频播放| 欧美日韩视频精品一区| 国产成人啪精品午夜网站| 无限看片的www在线观看| 久久久久久人人人人人| 两个人免费观看高清视频| 亚洲av电影在线进入| 久久精品亚洲av国产电影网| 亚洲精品国产一区二区精华液| 亚洲国产欧美网| 国产无遮挡羞羞视频在线观看| 人妻 亚洲 视频| 最近最新中文字幕大全免费视频| 日本一区二区免费在线视频| 亚洲 欧美一区二区三区| 国产精品99久久99久久久不卡| 国产又爽黄色视频| 国产有黄有色有爽视频| 国产精品秋霞免费鲁丝片| 啦啦啦中文免费视频观看日本| 亚洲三区欧美一区| 国产又爽黄色视频| 欧美激情高清一区二区三区| 国产极品粉嫩免费观看在线| 一区二区三区四区激情视频| 国产精品99久久99久久久不卡| 国产野战对白在线观看| 大香蕉久久网| a在线观看视频网站| 欧美少妇被猛烈插入视频| 男人爽女人下面视频在线观看| 亚洲少妇的诱惑av| 日韩熟女老妇一区二区性免费视频| 如日韩欧美国产精品一区二区三区| 国产高清videossex| 欧美日韩精品网址| 欧美在线黄色| 51午夜福利影视在线观看| 久久久久国产一级毛片高清牌| 人妻人人澡人人爽人人| 亚洲少妇的诱惑av| 性高湖久久久久久久久免费观看| 日本av免费视频播放| 精品少妇久久久久久888优播| 亚洲精品国产区一区二| 国产伦理片在线播放av一区| 久久精品亚洲av国产电影网| 午夜久久久在线观看| 欧美国产精品va在线观看不卡| 国产在线观看jvid| 亚洲人成77777在线视频| 女警被强在线播放| 日日爽夜夜爽网站| 18禁观看日本| 国产欧美日韩精品亚洲av| 欧美亚洲 丝袜 人妻 在线| 黑人欧美特级aaaaaa片| 亚洲中文字幕日韩| 深夜精品福利| 如日韩欧美国产精品一区二区三区| 黄片播放在线免费| 久热爱精品视频在线9| 午夜福利乱码中文字幕| 美女大奶头黄色视频| 国产精品欧美亚洲77777| av免费在线观看网站| 考比视频在线观看| 国产精品国产三级国产专区5o| 国产又色又爽无遮挡免| 中亚洲国语对白在线视频| 精品一区在线观看国产| 精品少妇一区二区三区视频日本电影| 欧美黑人精品巨大| 高清黄色对白视频在线免费看| 午夜成年电影在线免费观看| 水蜜桃什么品种好| 亚洲精品一二三| 国产一区二区三区av在线| xxxhd国产人妻xxx| 午夜精品国产一区二区电影| 国产精品一区二区在线观看99| 国产成人一区二区三区免费视频网站| 中文字幕高清在线视频| 新久久久久国产一级毛片| 中文字幕另类日韩欧美亚洲嫩草| 18禁裸乳无遮挡动漫免费视频| 欧美激情久久久久久爽电影 | 夫妻午夜视频| 亚洲一区中文字幕在线| h视频一区二区三区| 亚洲国产欧美一区二区综合| 久久国产精品人妻蜜桃| 精品国产超薄肉色丝袜足j| 欧美午夜高清在线| 欧美成狂野欧美在线观看| 我的亚洲天堂| 成人国产一区最新在线观看| 久久国产精品大桥未久av| 中文字幕色久视频| 久久精品亚洲熟妇少妇任你| 一区福利在线观看| 十八禁人妻一区二区| 看免费av毛片| 91字幕亚洲| 91av网站免费观看| 午夜久久久在线观看| 一级黄色大片毛片| 一本一本久久a久久精品综合妖精| 别揉我奶头~嗯~啊~动态视频 | 国产在线视频一区二区| 69精品国产乱码久久久| 天天躁日日躁夜夜躁夜夜| 男女无遮挡免费网站观看| 日韩人妻精品一区2区三区| 中国国产av一级| 日韩欧美国产一区二区入口| 在线看a的网站| 精品一区在线观看国产| 激情视频va一区二区三区| 久久久国产精品麻豆| 亚洲免费av在线视频| 欧美精品亚洲一区二区| 国产欧美日韩精品亚洲av| 国产成人av教育| 亚洲专区字幕在线| 国产成人系列免费观看| 欧美精品人与动牲交sv欧美| 男女高潮啪啪啪动态图| 超碰97精品在线观看| 国产免费现黄频在线看| e午夜精品久久久久久久| a 毛片基地| 黄片小视频在线播放| 99久久国产精品久久久| 免费日韩欧美在线观看| 久久天堂一区二区三区四区| 如日韩欧美国产精品一区二区三区| 大片免费播放器 马上看| 欧美97在线视频| 制服诱惑二区| 久久精品国产综合久久久| 日韩电影二区| 嫁个100分男人电影在线观看| 丁香六月欧美| 国产精品久久久av美女十八| 汤姆久久久久久久影院中文字幕| 亚洲欧美日韩高清在线视频 | 午夜精品国产一区二区电影| 国产伦人伦偷精品视频| 国产高清国产精品国产三级| 成人国产一区最新在线观看| 久久精品国产亚洲av高清一级| 国产高清videossex| 91字幕亚洲| 亚洲一码二码三码区别大吗| 黄色a级毛片大全视频| 免费黄频网站在线观看国产| 免费看十八禁软件| 久久国产亚洲av麻豆专区| 一个人免费看片子| 国产免费福利视频在线观看| 色视频在线一区二区三区| 精品人妻一区二区三区麻豆| 伊人亚洲综合成人网| 窝窝影院91人妻| 欧美在线黄色| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品国产av成人精品| 另类亚洲欧美激情| 视频区图区小说| 午夜福利影视在线免费观看| 国产不卡av网站在线观看| 久久久精品国产亚洲av高清涩受| 欧美激情久久久久久爽电影 | 黄色片一级片一级黄色片| 久久午夜综合久久蜜桃| www日本在线高清视频| 午夜福利,免费看| 少妇的丰满在线观看| 黄网站色视频无遮挡免费观看| 国产精品av久久久久免费| 亚洲精品国产色婷婷电影| 亚洲 国产 在线| 满18在线观看网站| 高潮久久久久久久久久久不卡| 丁香六月天网| 亚洲五月婷婷丁香| 一本综合久久免费| 两个人免费观看高清视频| 成年人黄色毛片网站| 美国免费a级毛片| 精品乱码久久久久久99久播| 99久久国产精品久久久| 水蜜桃什么品种好| 国产精品一区二区在线观看99| 欧美另类一区| www.自偷自拍.com| 男女边摸边吃奶| 日韩免费高清中文字幕av| 日本精品一区二区三区蜜桃| 首页视频小说图片口味搜索| 亚洲精品国产区一区二| 亚洲欧美精品自产自拍| 啦啦啦中文免费视频观看日本| 青春草视频在线免费观看| 97精品久久久久久久久久精品| 他把我摸到了高潮在线观看 | 亚洲熟女精品中文字幕| 女性被躁到高潮视频| 日韩,欧美,国产一区二区三区| www.精华液| 18禁观看日本| 男女边摸边吃奶| 精品久久久精品久久久| 另类亚洲欧美激情| 精品国产乱码久久久久久小说| 国产精品免费视频内射| 国产精品一区二区在线观看99| 亚洲精品美女久久久久99蜜臀| 电影成人av| av在线播放精品| 久久毛片免费看一区二区三区| 免费不卡黄色视频| 啦啦啦啦在线视频资源| 亚洲国产精品一区二区三区在线| 亚洲 欧美一区二区三区| 欧美黑人精品巨大| 最黄视频免费看| 亚洲熟女毛片儿| 亚洲av日韩在线播放| 免费一级毛片在线播放高清视频 | 91精品国产国语对白视频| 一区二区日韩欧美中文字幕| 操出白浆在线播放| 日韩 亚洲 欧美在线| 亚洲av电影在线进入| 久久综合国产亚洲精品| 日韩 欧美 亚洲 中文字幕| 精品国产一区二区久久| 亚洲欧美日韩另类电影网站| 性色av乱码一区二区三区2| 啦啦啦中文免费视频观看日本| 精品人妻在线不人妻| 日本vs欧美在线观看视频| 久久人人爽人人片av| 99国产综合亚洲精品| av在线播放精品| 午夜福利乱码中文字幕| 国产成人免费观看mmmm| 国产无遮挡羞羞视频在线观看| 免费女性裸体啪啪无遮挡网站| 少妇猛男粗大的猛烈进出视频| 考比视频在线观看| 成人国语在线视频| 国产精品久久久久久精品古装| 国产精品久久久久久人妻精品电影 | 久久av网站| 大片电影免费在线观看免费| 欧美 日韩 精品 国产| 最近最新中文字幕大全免费视频| 日本av手机在线免费观看| 波多野结衣av一区二区av| 美女主播在线视频| 色婷婷久久久亚洲欧美| kizo精华| 免费观看a级毛片全部| 欧美精品av麻豆av| 美女福利国产在线| 中文字幕精品免费在线观看视频| 777久久人妻少妇嫩草av网站| av不卡在线播放| 国产亚洲av高清不卡| 老汉色av国产亚洲站长工具| 亚洲国产欧美在线一区| 国产深夜福利视频在线观看| 国产精品偷伦视频观看了| 日韩 亚洲 欧美在线| 满18在线观看网站| 天天躁日日躁夜夜躁夜夜| 亚洲中文av在线| 9色porny在线观看| 日韩视频在线欧美| 十分钟在线观看高清视频www| 亚洲精品在线美女| a级毛片在线看网站| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品 欧美亚洲| 黑人巨大精品欧美一区二区蜜桃| 欧美成狂野欧美在线观看| 一级a爱视频在线免费观看| 99国产精品一区二区蜜桃av | 老鸭窝网址在线观看| 99国产精品一区二区三区| 黑丝袜美女国产一区| 窝窝影院91人妻| 午夜福利,免费看| 国产精品免费大片| 在线亚洲精品国产二区图片欧美| 最近最新中文字幕大全免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 黑人欧美特级aaaaaa片| h视频一区二区三区| 91麻豆av在线| 亚洲欧美色中文字幕在线| 少妇的丰满在线观看| 国产在线观看jvid| 色婷婷av一区二区三区视频| 极品少妇高潮喷水抽搐| 女人精品久久久久毛片| 老熟女久久久| 一本—道久久a久久精品蜜桃钙片| 极品少妇高潮喷水抽搐| 日本撒尿小便嘘嘘汇集6| 日韩欧美国产一区二区入口| 黄色片一级片一级黄色片| 极品少妇高潮喷水抽搐| 亚洲国产av影院在线观看| 欧美成人午夜精品| 一区二区三区精品91| 亚洲成人免费电影在线观看| 久久九九热精品免费| 性色av一级| 亚洲av美国av| 欧美另类亚洲清纯唯美| 波多野结衣av一区二区av| 精品久久蜜臀av无| 久久精品国产a三级三级三级| av超薄肉色丝袜交足视频| 成年人午夜在线观看视频| 天天操日日干夜夜撸| 九色亚洲精品在线播放| 在线看a的网站| av不卡在线播放| 伊人亚洲综合成人网| 亚洲精品一卡2卡三卡4卡5卡 | 国产精品一区二区在线不卡| www.av在线官网国产| 色综合欧美亚洲国产小说| 丝瓜视频免费看黄片| 热99re8久久精品国产| 成年av动漫网址| 亚洲av男天堂| 国产又爽黄色视频| 可以免费在线观看a视频的电影网站| 巨乳人妻的诱惑在线观看| 少妇的丰满在线观看| av在线老鸭窝| 亚洲第一欧美日韩一区二区三区 | 日韩 亚洲 欧美在线| 亚洲精品久久久久久婷婷小说| 亚洲精品一卡2卡三卡4卡5卡 | av一本久久久久| 亚洲,欧美精品.| 日日摸夜夜添夜夜添小说| 国产免费现黄频在线看| 黑人操中国人逼视频| 成年美女黄网站色视频大全免费| 色婷婷久久久亚洲欧美| 中亚洲国语对白在线视频| 亚洲熟女毛片儿| 国产欧美日韩一区二区精品| 精品人妻1区二区| 精品欧美一区二区三区在线| 蜜桃在线观看..| 久久青草综合色| 久久久欧美国产精品| 亚洲情色 制服丝袜| 成人免费观看视频高清| 国产免费一区二区三区四区乱码| 国产又爽黄色视频| 亚洲精品国产一区二区精华液| 久久久欧美国产精品| 午夜免费鲁丝| 亚洲人成77777在线视频| 免费在线观看视频国产中文字幕亚洲 | 后天国语完整版免费观看| 久久久久久久国产电影| 性少妇av在线| 波多野结衣av一区二区av| www.精华液| 国产伦理片在线播放av一区| 少妇猛男粗大的猛烈进出视频| 国产亚洲精品久久久久5区| 国产精品熟女久久久久浪| 亚洲一区二区三区欧美精品| 国产一级毛片在线| 99久久精品国产亚洲精品| 久久狼人影院| 久久午夜综合久久蜜桃| 性高湖久久久久久久久免费观看| 国产成+人综合+亚洲专区| 69av精品久久久久久 | 动漫黄色视频在线观看| 日本一区二区免费在线视频|