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

    阿勒泰地區(qū)降水時(shí)間序列的多重分形特征分析

    2021-04-29 02:44:14周天祥
    關(guān)鍵詞:布爾津福海阿勒泰

    周天祥 陳 晶

    (1.紹興文理學(xué)院 機(jī)械與電氣工程學(xué)院,浙江 紹興 312000;2.紹興文理學(xué)院 土木工程學(xué)院,浙江 紹興 312000)

    自然界存在著許多復(fù)雜的動(dòng)力學(xué)過(guò)程,這些復(fù)雜的相互作用通過(guò)一系列看似嘈雜且不穩(wěn)定的信號(hào)表現(xiàn)出來(lái),通??梢允褂梅蔷€性微分方程來(lái)描述這些過(guò)程.一些非線性系統(tǒng)的時(shí)間序列信號(hào)存在著顯著的分形特征[1],即局部時(shí)間序列與整個(gè)時(shí)間序列完全相同或近似.經(jīng)過(guò)許多不同領(lǐng)域?qū)W者的努力,這種現(xiàn)象目前在經(jīng)濟(jì)學(xué)、生理學(xué)、地質(zhì)學(xué)以及氣象學(xué)等學(xué)科的復(fù)雜系統(tǒng)中均有發(fā)現(xiàn)[2-5].大氣中存在很多不規(guī)則的隨機(jī)運(yùn)動(dòng)(如大氣湍流),可將其視為在各種尺度的天氣系統(tǒng)相互作用下的復(fù)雜系統(tǒng),這些相互作用將導(dǎo)致在許多時(shí)空尺度上氣象參數(shù)以近似隨機(jī)的方式波動(dòng).Pielke等人[6]使用非線性微分方程對(duì)氣象過(guò)程進(jìn)行建模,但是數(shù)值方法與現(xiàn)實(shí)情況仍存在差距.

    降水的時(shí)空變化對(duì)農(nóng)業(yè)生產(chǎn)、水資源管理以及生態(tài)系統(tǒng)有著深刻的影響,同時(shí)也是氣候變化的重要參考指標(biāo).降雨量依賴(lài)于大氣過(guò)程的多樣性變化,并由于自然過(guò)程的復(fù)雜性在不同的時(shí)空尺度上呈現(xiàn)出隨機(jī)波動(dòng).一些學(xué)者基于統(tǒng)計(jì)學(xué)如Mann-Kendall秩次相關(guān)檢驗(yàn)、Pettitt檢驗(yàn)、Buishand檢驗(yàn)等方法對(duì)降水時(shí)間序列的變化趨勢(shì)進(jìn)行評(píng)價(jià)[7-9],然而這些方法無(wú)法捕捉降水時(shí)間序列復(fù)雜非線性行為中的潛在特征.隨著非線性動(dòng)力學(xué)的發(fā)展,多重分形級(jí)聯(lián)模型逐漸被應(yīng)用于研究天氣與氣候變化.Lovejoy和Mandel-brot[10]發(fā)現(xiàn)降雨區(qū)域的云層邊界投影到地面的形狀為分形圖形,并且降水序列的時(shí)空結(jié)構(gòu)服從雙曲線分布,從而確立了分形理論在氣象學(xué)的適用性.Telesca等人[11]應(yīng)用去趨勢(shì)波動(dòng)分析方法(DFA)對(duì)阿根廷中部地區(qū)近60年的降水時(shí)間序列進(jìn)行研究,發(fā)現(xiàn)了降水時(shí)間序列中的標(biāo)度不變性以及長(zhǎng)程相關(guān)性.然而DFA方法適用于描述單個(gè)標(biāo)度指數(shù),許多復(fù)雜系統(tǒng)中的信號(hào)并非表現(xiàn)出簡(jiǎn)單的單分形特征,而是在不同的時(shí)間尺度上具有不同的分形結(jié)構(gòu),即多重分形行為.多重分形去趨勢(shì)波動(dòng)分析方法(MF-DFA)能更加精確地描述復(fù)雜系統(tǒng)可能存在的多個(gè)分形特征[12].Kalamaras等人[13]應(yīng)用MF-DFA方法對(duì)希臘克里特島沿海氣象站的氣溫觀測(cè)數(shù)據(jù)進(jìn)行了研究,發(fā)現(xiàn)該地區(qū)的氣溫時(shí)間序列表現(xiàn)出多重分形特征,并且氣溫時(shí)間序列為長(zhǎng)期正相關(guān).Olsson[14]運(yùn)用概率密度函數(shù)以及功率譜分析了降水時(shí)間序列所服從的冪律分布,選取的不同統(tǒng)計(jì)矩階數(shù)隨尺度而變化證明該地區(qū)的降水時(shí)間序列具有多重分形特征.新疆北部地區(qū)水資源稀缺,降雨量的變化對(duì)該地區(qū)農(nóng)牧業(yè)的生產(chǎn)至關(guān)重要.本文基于新疆北部阿勒泰地區(qū)1961至2015年的逐月降雨量數(shù)據(jù),對(duì)該降水時(shí)間序列的分形特征和變化趨勢(shì)運(yùn)用MF-DFA方法進(jìn)行了詳細(xì)分析.

    1 研究區(qū)域概況及數(shù)據(jù)源

    本文研究資料來(lái)源于額爾齊斯河流域標(biāo)準(zhǔn)氣象站觀測(cè)數(shù)據(jù)[15],研究的三個(gè)地點(diǎn)分別位于新疆北部阿爾泰地區(qū)的阿勒泰氣象站(47.73 °N,88.08 °E)、福海氣象站(47.12 °N,87.47 °E)、布爾津氣象站(47.70 °N,86.87 °E),所在位置已在圖1中標(biāo)明.該區(qū)域深居歐亞大陸腹地,屬于溫帶大陸性氣候,冬季在強(qiáng)大的蒙古高壓控制下氣候寒冷干燥,年平均降水量在100 mm左右,降水多集中在夏季,冬季降雨量較春秋季節(jié)更多.三個(gè)觀測(cè)站的降水資料數(shù)據(jù)的時(shí)間范圍為1961年1月至2015年12月,觀測(cè)數(shù)據(jù)的時(shí)間尺度為月.

    2 MF-DFA原理

    分形是指局部與整體之間具有某種方式相似的形狀,即局部片段與整體之間具有自相似性,目前分形理論已從狹義的幾何對(duì)象擴(kuò)展至廣義分形.Mandelbrot[16]在研究美國(guó)棉花價(jià)格的變化時(shí)發(fā)現(xiàn)價(jià)格時(shí)間序列的短期波動(dòng)與長(zhǎng)期的變化趨勢(shì)存在標(biāo)度不變性,因而發(fā)現(xiàn)了時(shí)間序列中的自相似現(xiàn)象.若分形維隨時(shí)間尺度的不同而發(fā)生變化,則該時(shí)間序列具備多重分形特性.多重分形去趨勢(shì)波動(dòng)分析方法(MF-DFA)由Kandelhardt等人[12]在消除趨勢(shì)波動(dòng)分析(DFA)方法的基礎(chǔ)上拓展得到,能更加精確地描述復(fù)雜系統(tǒng)的動(dòng)力學(xué)行為,具體計(jì)算方法如下:

    第一步,首先對(duì)降水時(shí)間序列xi(i=1,2,…,n)求取累積離差序列:

    (1)

    第二步,將Yi使用長(zhǎng)度為s的Ns=int(N/s)個(gè)互不重疊的區(qū)間,對(duì)于每個(gè)長(zhǎng)度為s的區(qū)間進(jìn)行最小二乘擬合,得到該區(qū)間的局部趨勢(shì)多項(xiàng)式y(tǒng)v(s).對(duì)Yi上的每個(gè)子區(qū)間v計(jì)算局部均方誤差F2(v,s):

    圖1 研究區(qū)域所在位置

    (2)

    第三步,第qth階波動(dòng)函數(shù)可表示為:

    (3)

    第四步,對(duì)指定的q階值可使用下式確定波動(dòng)函數(shù)的標(biāo)度系數(shù):

    Fq(s)∝sH(q)

    (4)

    式中H(q)為廣義Hurst指數(shù),其取值范圍為[0,1],H(q)越接近1,表明時(shí)間序列存在的噪聲越少,趨勢(shì)越明顯.當(dāng)H=0.5時(shí),序列變化為標(biāo)準(zhǔn)布朗運(yùn)動(dòng),序列之間不相關(guān),表明當(dāng)前降水時(shí)間序列的趨勢(shì)不會(huì)影響未來(lái)的情況;當(dāng)H(q)>0.5時(shí),表明當(dāng)前降水時(shí)間序列與未來(lái)增量之間為正相關(guān);當(dāng)H(q)<0.5時(shí),表明當(dāng)前的降水時(shí)間序列與未來(lái)增量之間為負(fù)相關(guān).若H(q)指為常數(shù)時(shí)表明時(shí)間序列為單分形的情況,當(dāng)H(q)指數(shù)隨q值的增大而減小時(shí)表明時(shí)間序列屬于多重分形情況[17].

    第五步,計(jì)算多重分形標(biāo)度系數(shù)τ(q)滿(mǎn)足下列關(guān)系式:

    τ(q)=qH(q)-1

    (5)

    第六步,根據(jù)Legendre變換,多重分形奇異譜函數(shù)f(α)滿(mǎn)足下式:

    f(α)=qα-τ(q)

    (6)

    其中α=?τ(q)/?(q)為奇異性指數(shù),奇異譜的譜寬(Δα=αmax-αmin)描述了多重分形的程度,譜寬越大表明多重分形程度越高,降水時(shí)間序列波動(dòng)越劇烈.若時(shí)間序列具有多重分形屬性,其奇異譜將呈現(xiàn)為倒拋物線形[18].Δf=f(αmin)-f(αmax)表示最大概率子集與最小概率子集的分形差異,Δf可以反映時(shí)間序列的變化趨勢(shì),Δf<0表明未來(lái)時(shí)間序列數(shù)據(jù)的向低數(shù)值方向發(fā)展的概率較大[19].

    3 實(shí)驗(yàn)結(jié)果分析

    本文選取新疆阿勒泰地區(qū)的阿勒泰、福海、布爾津三個(gè)氣象觀測(cè)站的逐月觀測(cè)數(shù)據(jù),每個(gè)降水觀測(cè)數(shù)據(jù)的時(shí)間序列長(zhǎng)度都是660.由圖2可知,三個(gè)地區(qū)的歷年逐月降雨量基本都不足100 mm,降水時(shí)間序列全都表現(xiàn)出非線性、非平穩(wěn)性的特點(diǎn).福海地區(qū)降雨量最少,且在1984年出現(xiàn)一次大降水天氣,降水量超過(guò)了100 mm.阿勒泰地區(qū)降雨量趨勢(shì)受季節(jié)變化影響明顯,夏季降雨量要明顯高于其他季節(jié),因此應(yīng)先將原始的降水時(shí)間序列轉(zhuǎn)化為隨機(jī)游走結(jié)構(gòu)[20]以剔除季節(jié)更替帶來(lái)的影響.通過(guò)對(duì)三個(gè)氣象站的降水時(shí)間序列xi計(jì)算累積離差(式2)可得到隨機(jī)游走的降水時(shí)間序列Yi,如圖3所示,總體上降水累積離差序列在20世紀(jì)60年代至90年代呈下降趨勢(shì),而在60年代末和80年代出現(xiàn)了兩個(gè)小幅震蕩的時(shí)期,90年代后又呈增加趨勢(shì).

    (a)阿勒泰 (b)福海 (c)布爾津

    (a)阿勒泰 (b)福海 (c)布爾津

    對(duì)三個(gè)地區(qū)降水隨機(jī)游走時(shí)間序列Yi按照時(shí)間尺度s進(jìn)行分割,為防止多項(xiàng)式過(guò)擬合降水時(shí)間序列從而導(dǎo)致剩余波動(dòng)接近于0,時(shí)間尺度s應(yīng)滿(mǎn)足不小于多項(xiàng)式擬合階數(shù),同時(shí)最大時(shí)間尺度應(yīng)是時(shí)間序列劃分區(qū)間的個(gè)數(shù)大于10個(gè)[21].統(tǒng)計(jì)矩階數(shù)q的取值范圍對(duì)多重分形譜的計(jì)算十分關(guān)鍵,q值過(guò)大會(huì)導(dǎo)致計(jì)算時(shí)間呈指數(shù)增大,q值過(guò)小時(shí)多重分形譜的波動(dòng)會(huì)很大[22].由此本文選取Yi的子區(qū)間長(zhǎng)度s=(6∶66),按照經(jīng)驗(yàn)取q=(-10∶10),利用式(3)對(duì)不同取值的統(tǒng)計(jì)矩階數(shù)q計(jì)算波動(dòng)函數(shù)Fq(s).波動(dòng)函數(shù)Fq(s)關(guān)于時(shí)間尺度s在特定統(tǒng)計(jì)矩階數(shù)q=(-10,0,10)下的雙對(duì)數(shù)圖像如圖4所示,所有q值情況下的波動(dòng)函數(shù)Fq(s)關(guān)于時(shí)間尺度s都近似為線性關(guān)系,這表明三個(gè)地區(qū)的降水時(shí)間序列都符合冪律分布,同時(shí)也證明了降水時(shí)間序列的內(nèi)部波動(dòng)演化并非隨機(jī)行為,而是由遠(yuǎn)距離的某些特征有序引發(fā).另外,從圖4可以明顯看出,回歸線的斜率明顯隨q值的增大而減小,由式(4)可知該斜率即為廣義Hurst指數(shù),圖4中的擬合方程以及可決系數(shù)見(jiàn)表1.q值描述了局部波動(dòng)均方根占整體的比重[21],當(dāng)q值為正時(shí),廣義Hurst指數(shù)表示大波動(dòng)的標(biāo)度行為,q值為負(fù)時(shí),廣義Hurst指數(shù)表示小波動(dòng)的標(biāo)度行為,對(duì)于較小的時(shí)間尺度s受降水時(shí)間序列波動(dòng)幅度影響較大,而在較大的時(shí)間尺度上降水時(shí)間序列的波動(dòng)幅度將被削弱[22].由圖2(b)可以看出,福海地區(qū)的降水時(shí)間序列出現(xiàn)過(guò)一次大降水天氣情況,同時(shí)其累積離差分布也存在多個(gè)波動(dòng)趨勢(shì),降水時(shí)間序列受大的波動(dòng)控制,因此其廣義Hurst指數(shù)變化范圍最大(ΔH(q)=0.57).

    (a)阿勒泰 (b)福海 (c)布爾津

    廣義Hurst指數(shù)可用于判斷降水時(shí)間序列的分形特征以及長(zhǎng)程相關(guān)性.如圖5所示,三個(gè)地區(qū)的Hurst指數(shù)關(guān)于q為非線性關(guān)系,說(shuō)明三個(gè)地區(qū)的降水時(shí)間序列都存在多重分形特征.當(dāng)q值為2時(shí),廣義Hurst指數(shù)與經(jīng)典的Hurst指數(shù)相同,此時(shí)阿勒泰、福海、布爾津地區(qū)降水時(shí)間序列的廣義Hurst指數(shù)分別為0.56、0.50、0.52,因此當(dāng)q∈[-10,2]時(shí),H(q)>0.5,三個(gè)地區(qū)的降水時(shí)間序列都具有長(zhǎng)期正相關(guān)性.由圖5可知,阿勒泰與布爾津的降水時(shí)間序列的長(zhǎng)程相關(guān)性明顯強(qiáng)于其他福海地區(qū),表明未來(lái)的該地區(qū)降水趨勢(shì)將與當(dāng)前保持一致.

    表1 降水時(shí)間序列廣義Hurst指數(shù)的計(jì)算

    圖5 降水時(shí)間序列的廣義Hurst指數(shù)

    根據(jù)(式11)計(jì)算三個(gè)地區(qū)降水時(shí)間序列的多重分形譜如圖6所示,其多重分形譜均呈倒拋物線形,分形譜的最大值對(duì)應(yīng)為q=0時(shí).由表2可知,福海地區(qū)的譜寬最大,Δα=0.96,因此該地區(qū)降水時(shí)間序列的多重分形程度要強(qiáng)于阿勒泰以及布爾津地區(qū),并且其降水時(shí)間序列波動(dòng)幅度最劇烈.阿勒泰氣象站和布爾津氣象站的譜寬比較接近,因此其多重分形程度也相似.阿勒泰和福海的降水時(shí)間序列的譜高度Δf<0,表明這兩個(gè)地區(qū)未來(lái)降雨量可能有降低的趨勢(shì).布爾津地區(qū)降水時(shí)間序列的譜高度Δf>0,因此未來(lái)該地區(qū)降雨量增加的概率較大.

    表2 降水時(shí)間序列的多重分形參數(shù)

    (a)阿勒泰 (b)福海 (c)布爾津

    4 結(jié)論

    本文通過(guò)對(duì)新疆北部阿勒泰地區(qū)的阿勒泰、福海、布爾津三個(gè)氣象站的1961年至2015年的逐月降雨量時(shí)間序列的變化進(jìn)行多重分形檢驗(yàn),結(jié)論如下:

    (1)根據(jù)MF-DFA方法的檢驗(yàn)結(jié)果,三個(gè)地區(qū)的降水資料時(shí)間序列都表現(xiàn)出顯著的多重分形特征,符合冪律相關(guān)特性;

    (2)計(jì)算廣義Hurst指數(shù)發(fā)現(xiàn)降水時(shí)間序列在q<2時(shí)存在明顯的長(zhǎng)期正相關(guān),表明未來(lái)的降水趨勢(shì)具備可預(yù)測(cè)性;

    (3)福海地區(qū)的降水時(shí)間序列相比另外兩個(gè)地區(qū)的多重分形譜寬更大、波動(dòng)更強(qiáng)烈,表現(xiàn)出更強(qiáng)的多重分形程度;

    (4)對(duì)于長(zhǎng)期的降水時(shí)間序列使用MF-DFA方法,可以獲取其統(tǒng)計(jì)學(xué)的表征以便更好地理解其動(dòng)力學(xué)機(jī)制,為未來(lái)當(dāng)?shù)亟邓淖兓厔?shì)提供參考.

    猜你喜歡
    布爾津福海阿勒泰
    童話城堡新來(lái)的朋友們
    延河(2022年6期)2022-07-24 21:03:50
    倉(cāng)央嘉措
    孫福海書(shū)法作品欣賞
    印象·阿勒泰
    寶藏(2019年9期)2019-01-15 02:01:16
    我與福海有個(gè)約會(huì)
    收藏界(2018年5期)2018-10-08 09:10:32
    論阿勒泰烏梁海部弓箭的崇拜
    任福海:烏鴉反哺 心懷感恩
    公民與法治(2016年9期)2016-05-17 04:12:35
    阿勒泰水電站“以閥代井”研究
    劉志:三點(diǎn)導(dǎo)讀
    阿勒泰考古之行
    大眾考古(2015年7期)2015-06-26 08:41:00
    嫩草影院入口| 免费看日本二区| 最近在线观看免费完整版| 一区二区三区免费毛片| 男人和女人高潮做爰伦理| 在线观看日韩欧美| 可以在线观看的亚洲视频| 久久久久性生活片| 久久草成人影院| 婷婷丁香在线五月| 狂野欧美白嫩少妇大欣赏| 国产97色在线日韩免费| 桃色一区二区三区在线观看| 亚洲精品456在线播放app | 日本黄大片高清| 小说图片视频综合网站| 九九在线视频观看精品| 亚洲精品在线美女| 91字幕亚洲| 好男人电影高清在线观看| 午夜免费成人在线视频| 亚洲美女黄片视频| 搡女人真爽免费视频火全软件 | 国内精品美女久久久久久| 中文亚洲av片在线观看爽| 熟女电影av网| 免费无遮挡裸体视频| 久久久久国产精品人妻aⅴ院| 国产亚洲精品av在线| 欧美最黄视频在线播放免费| 日本一本二区三区精品| 亚洲人成网站高清观看| 亚洲不卡免费看| 国产成人欧美在线观看| 一级黄片播放器| 性色avwww在线观看| 国产精品 国内视频| 性色av乱码一区二区三区2| 国产精品爽爽va在线观看网站| www.色视频.com| 人人妻,人人澡人人爽秒播| 欧美日韩综合久久久久久 | 18禁黄网站禁片免费观看直播| 亚洲一区二区三区色噜噜| 欧美+日韩+精品| 搡老熟女国产l中国老女人| 香蕉久久夜色| 中文字幕久久专区| 黄色丝袜av网址大全| 午夜福利成人在线免费观看| 国产色婷婷99| 狂野欧美激情性xxxx| 国产日本99.免费观看| 禁无遮挡网站| 性色avwww在线观看| 男女视频在线观看网站免费| 国产精品精品国产色婷婷| www.999成人在线观看| 欧美黄色片欧美黄色片| 久久久久久久久大av| 悠悠久久av| 精品国产亚洲在线| 亚洲av成人不卡在线观看播放网| 亚洲美女黄片视频| 久久精品国产亚洲av香蕉五月| 国产欧美日韩一区二区三| 欧美性感艳星| 亚洲精品久久国产高清桃花| 美女 人体艺术 gogo| 亚洲av日韩精品久久久久久密| 黄片小视频在线播放| 热99在线观看视频| 色综合站精品国产| 国产高清视频在线观看网站| 国产精品一及| 成人三级黄色视频| 国产高清视频在线播放一区| 女同久久另类99精品国产91| 亚洲无线观看免费| 日日夜夜操网爽| 国产aⅴ精品一区二区三区波| 老司机福利观看| 国产麻豆成人av免费视频| 99久久无色码亚洲精品果冻| 一边摸一边抽搐一进一小说| 国产精品 国内视频| 亚洲精华国产精华精| 嫩草影院精品99| 久久久久久久午夜电影| 久久精品国产清高在天天线| 在线视频色国产色| 色噜噜av男人的天堂激情| 亚洲精品日韩av片在线观看 | eeuss影院久久| 男女下面进入的视频免费午夜| 亚洲精品在线观看二区| 少妇人妻一区二区三区视频| 日韩成人在线观看一区二区三区| 久久久久久九九精品二区国产| 身体一侧抽搐| 欧美日韩综合久久久久久 | 波多野结衣高清无吗| 在线观看免费视频日本深夜| 夜夜看夜夜爽夜夜摸| 99久久九九国产精品国产免费| av专区在线播放| 最近最新中文字幕大全电影3| 精品久久久久久久久久久久久| 观看免费一级毛片| 天堂网av新在线| 两个人看的免费小视频| 久久人妻av系列| 精品久久久久久成人av| 熟女少妇亚洲综合色aaa.| 人妻久久中文字幕网| 国产成人欧美在线观看| 亚洲男人的天堂狠狠| 久久草成人影院| 亚洲久久久久久中文字幕| 国产又黄又爽又无遮挡在线| 最近最新中文字幕大全电影3| 亚洲欧美日韩高清在线视频| 国产欧美日韩精品一区二区| 在线观看舔阴道视频| 色老头精品视频在线观看| 国产老妇女一区| 日本免费一区二区三区高清不卡| 神马国产精品三级电影在线观看| 国产野战对白在线观看| 亚洲国产精品999在线| 欧美bdsm另类| 久久久久久久午夜电影| 国产一区二区激情短视频| 欧美日韩中文字幕国产精品一区二区三区| 噜噜噜噜噜久久久久久91| 九九热线精品视视频播放| 搡老妇女老女人老熟妇| 女人十人毛片免费观看3o分钟| 免费观看的影片在线观看| 国产精品久久久久久久久免 | 亚洲精品粉嫩美女一区| 淫秽高清视频在线观看| 香蕉丝袜av| 嫩草影院入口| 黄片大片在线免费观看| 精品欧美国产一区二区三| 十八禁人妻一区二区| 久久久久国内视频| 欧美乱妇无乱码| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美激情综合另类| 欧洲精品卡2卡3卡4卡5卡区| 久久久精品大字幕| 99久久99久久久精品蜜桃| 亚洲国产精品合色在线| 黄色日韩在线| 老汉色∧v一级毛片| 嫩草影视91久久| 欧美bdsm另类| 高清日韩中文字幕在线| 久久久久久久久大av| 国产精品女同一区二区软件 | 少妇人妻一区二区三区视频| 午夜影院日韩av| 97超视频在线观看视频| 日韩av在线大香蕉| 少妇的丰满在线观看| 日本在线视频免费播放| 国产精品野战在线观看| 欧美不卡视频在线免费观看| 动漫黄色视频在线观看| 在线观看av片永久免费下载| 村上凉子中文字幕在线| 色吧在线观看| 亚洲性夜色夜夜综合| 69av精品久久久久久| 日韩欧美精品免费久久 | 18禁裸乳无遮挡免费网站照片| 精品欧美国产一区二区三| 欧美色欧美亚洲另类二区| 亚洲国产中文字幕在线视频| 亚洲av成人精品一区久久| 欧美高清成人免费视频www| 亚洲无线观看免费| 在线观看免费视频日本深夜| 女人高潮潮喷娇喘18禁视频| 天堂√8在线中文| 在线观看日韩欧美| 国产av麻豆久久久久久久| 日本黄大片高清| 长腿黑丝高跟| 可以在线观看的亚洲视频| 日日摸夜夜添夜夜添小说| 男人舔奶头视频| 五月玫瑰六月丁香| 午夜福利视频1000在线观看| 欧美一级毛片孕妇| 在线播放国产精品三级| 欧洲精品卡2卡3卡4卡5卡区| 99热精品在线国产| 日韩国内少妇激情av| 国内精品美女久久久久久| 18禁黄网站禁片免费观看直播| 欧美黄色片欧美黄色片| 99国产综合亚洲精品| 十八禁人妻一区二区| 亚洲人成伊人成综合网2020| 高清在线国产一区| 欧美高清成人免费视频www| 免费av不卡在线播放| 午夜福利高清视频| 美女黄网站色视频| 国产精品av视频在线免费观看| 国产高清videossex| 小说图片视频综合网站| 丰满乱子伦码专区| 午夜福利成人在线免费观看| 日本a在线网址| 国产69精品久久久久777片| 日本在线视频免费播放| 精品国内亚洲2022精品成人| 亚洲最大成人中文| 国产极品精品免费视频能看的| 色在线成人网| 黄色女人牲交| 久久久久久大精品| 国产精品 国内视频| 国产成人aa在线观看| 久久精品国产99精品国产亚洲性色| 法律面前人人平等表现在哪些方面| 免费观看精品视频网站| 欧美乱色亚洲激情| 天堂网av新在线| 看黄色毛片网站| 国产精品一区二区三区四区免费观看 | 国产色爽女视频免费观看| 精品一区二区三区av网在线观看| 搡老熟女国产l中国老女人| 国产亚洲精品久久久久久毛片| 精品国产三级普通话版| 99久久无色码亚洲精品果冻| x7x7x7水蜜桃| 在线播放国产精品三级| 看片在线看免费视频| 99久久九九国产精品国产免费| 中文字幕av成人在线电影| 3wmmmm亚洲av在线观看| 免费高清视频大片| 亚洲av熟女| 国产成人aa在线观看| 亚洲男人的天堂狠狠| 欧美最黄视频在线播放免费| 日韩国内少妇激情av| 国产极品精品免费视频能看的| 日韩成人在线观看一区二区三区| 美女高潮喷水抽搐中文字幕| 天天添夜夜摸| 久久久色成人| 少妇的丰满在线观看| 高潮久久久久久久久久久不卡| 国产亚洲精品一区二区www| 亚洲精品一卡2卡三卡4卡5卡| 免费在线观看影片大全网站| 一区福利在线观看| 国产一区二区亚洲精品在线观看| 最近最新免费中文字幕在线| 三级毛片av免费| 十八禁人妻一区二区| 国产综合懂色| 成人永久免费在线观看视频| h日本视频在线播放| 日韩欧美 国产精品| 亚洲av成人不卡在线观看播放网| 国产精品av视频在线免费观看| 久久久久久大精品| 国产精品 国内视频| 在线免费观看不下载黄p国产 | 亚洲精品一区av在线观看| 国产蜜桃级精品一区二区三区| 久久久久久久午夜电影| 亚洲五月婷婷丁香| 国产精品久久电影中文字幕| 亚洲成av人片在线播放无| 女同久久另类99精品国产91| 国产真实伦视频高清在线观看 | 欧美中文日本在线观看视频| 最近最新中文字幕大全电影3| 又黄又爽又免费观看的视频| 精品人妻1区二区| 亚洲国产欧洲综合997久久,| 亚洲欧美一区二区三区黑人| 操出白浆在线播放| 日本三级黄在线观看| 天堂av国产一区二区熟女人妻| 在线a可以看的网站| 在线视频色国产色| 欧美日韩中文字幕国产精品一区二区三区| 国产成人影院久久av| 成人亚洲精品av一区二区| 欧美乱码精品一区二区三区| 亚洲精品在线美女| 精品国产三级普通话版| 超碰av人人做人人爽久久 | 高清在线国产一区| 成人永久免费在线观看视频| 三级国产精品欧美在线观看| 国产精品一及| 国产精品久久久久久人妻精品电影| 欧美日韩福利视频一区二区| 色综合欧美亚洲国产小说| 国产精品电影一区二区三区| 日本黄色视频三级网站网址| 国产主播在线观看一区二区| 久9热在线精品视频| 免费人成在线观看视频色| 女人高潮潮喷娇喘18禁视频| 亚洲黑人精品在线| 国产综合懂色| 精品福利观看| 精华霜和精华液先用哪个| 亚洲av一区综合| 99热6这里只有精品| 真人做人爱边吃奶动态| 两个人的视频大全免费| 香蕉av资源在线| 69av精品久久久久久| 欧美黄色淫秽网站| 日日夜夜操网爽| av福利片在线观看| 美女免费视频网站| 乱人视频在线观看| 麻豆成人午夜福利视频| 日本三级黄在线观看| 日本熟妇午夜| 午夜福利在线在线| 欧美黄色片欧美黄色片| 久久天躁狠狠躁夜夜2o2o| 国产精品久久久久久久久免 | 成人午夜高清在线视频| 少妇的丰满在线观看| 无限看片的www在线观看| 国产三级在线视频| 嫩草影院精品99| 在线免费观看不下载黄p国产 | 免费电影在线观看免费观看| 日本免费a在线| 黄色片一级片一级黄色片| 国产精品一区二区免费欧美| av福利片在线观看| 精品国产三级普通话版| 亚洲av熟女| 亚洲一区高清亚洲精品| 亚洲色图av天堂| 国产精品,欧美在线| 午夜免费激情av| 欧美另类亚洲清纯唯美| 99国产精品一区二区三区| 国产av在哪里看| 听说在线观看完整版免费高清| 两人在一起打扑克的视频| 亚洲午夜理论影院| 女人高潮潮喷娇喘18禁视频| 欧美日韩中文字幕国产精品一区二区三区| 在线观看免费午夜福利视频| 亚洲在线自拍视频| 国产午夜福利久久久久久| 中文字幕人妻丝袜一区二区| 男人和女人高潮做爰伦理| 怎么达到女性高潮| 每晚都被弄得嗷嗷叫到高潮| 国产黄a三级三级三级人| 欧美日韩福利视频一区二区| 免费人成在线观看视频色| 日韩免费av在线播放| 久久香蕉国产精品| 在线观看日韩欧美| 国产色爽女视频免费观看| 一区二区三区激情视频| 久久久久久久久大av| 十八禁人妻一区二区| 亚洲在线观看片| 中文字幕av成人在线电影| 舔av片在线| 老司机午夜十八禁免费视频| 欧美一区二区精品小视频在线| 国产私拍福利视频在线观看| 日本一二三区视频观看| 久久久久精品国产欧美久久久| 全区人妻精品视频| 国产av不卡久久| 国产成人aa在线观看| 久久精品亚洲精品国产色婷小说| 欧美av亚洲av综合av国产av| 最近最新中文字幕大全免费视频| 成人亚洲精品av一区二区| 国产黄片美女视频| 久久精品国产清高在天天线| 国产亚洲av嫩草精品影院| 中文字幕久久专区| 国产三级中文精品| 99国产综合亚洲精品| 少妇的丰满在线观看| 亚洲国产精品合色在线| 美女免费视频网站| 亚洲国产精品成人综合色| 法律面前人人平等表现在哪些方面| 免费无遮挡裸体视频| 国产精品,欧美在线| 好男人在线观看高清免费视频| 美女被艹到高潮喷水动态| 无遮挡黄片免费观看| 亚洲国产精品成人综合色| 午夜久久久久精精品| 精品一区二区三区视频在线观看免费| 天堂动漫精品| 午夜激情福利司机影院| 搡老熟女国产l中国老女人| 亚洲精品一区av在线观看| 午夜亚洲福利在线播放| 婷婷丁香在线五月| 在线观看免费视频日本深夜| 亚洲精品国产精品久久久不卡| 操出白浆在线播放| 一本综合久久免费| 一a级毛片在线观看| 欧美成人一区二区免费高清观看| 九色国产91popny在线| 精品人妻偷拍中文字幕| 亚洲性夜色夜夜综合| 午夜福利成人在线免费观看| 亚洲av二区三区四区| 亚洲av免费在线观看| 熟女人妻精品中文字幕| 搡老熟女国产l中国老女人| 国产真人三级小视频在线观看| 日韩av在线大香蕉| 夜夜躁狠狠躁天天躁| 男人舔奶头视频| 日本免费a在线| 国产在视频线在精品| 老熟妇仑乱视频hdxx| 69人妻影院| 久久伊人香网站| 99久久久亚洲精品蜜臀av| 黄色视频,在线免费观看| 内地一区二区视频在线| 欧美日韩黄片免| 淫秽高清视频在线观看| 亚洲欧美精品综合久久99| 日本成人三级电影网站| 九色成人免费人妻av| 亚洲精品456在线播放app | 亚洲av成人av| 国产视频内射| 欧美区成人在线视频| 亚洲av一区综合| 三级毛片av免费| 国产精品1区2区在线观看.| 哪里可以看免费的av片| 国产69精品久久久久777片| av天堂中文字幕网| 国产精品爽爽va在线观看网站| 日本三级黄在线观看| 悠悠久久av| 88av欧美| 成年人黄色毛片网站| 国产高清三级在线| 亚洲av二区三区四区| 久久精品国产自在天天线| 俄罗斯特黄特色一大片| 国产日本99.免费观看| a级一级毛片免费在线观看| 1000部很黄的大片| 十八禁网站免费在线| 国产精品 国内视频| 99热这里只有精品一区| 黄色日韩在线| 亚洲av成人av| 一级黄片播放器| 激情在线观看视频在线高清| 日韩 欧美 亚洲 中文字幕| 69人妻影院| 内地一区二区视频在线| 网址你懂的国产日韩在线| 欧美在线黄色| 欧美性感艳星| 一级黄片播放器| 人人妻人人看人人澡| 久久久久亚洲av毛片大全| 亚洲无线在线观看| 操出白浆在线播放| 看黄色毛片网站| 波多野结衣高清无吗| 久久久久久大精品| 免费观看的影片在线观看| 国产av一区在线观看免费| 免费av不卡在线播放| 欧美三级亚洲精品| 精华霜和精华液先用哪个| 日韩欧美三级三区| 久久久久久大精品| 青草久久国产| 三级毛片av免费| 夜夜看夜夜爽夜夜摸| 男女之事视频高清在线观看| 亚洲在线自拍视频| 国产精品av视频在线免费观看| 一级黄片播放器| 亚洲一区高清亚洲精品| 三级男女做爰猛烈吃奶摸视频| 乱人视频在线观看| 精品久久久久久久人妻蜜臀av| 久久久久久国产a免费观看| 天堂网av新在线| 中文字幕人妻丝袜一区二区| 国产高清videossex| 国产免费av片在线观看野外av| 少妇人妻精品综合一区二区 | 美女免费视频网站| 欧美乱妇无乱码| 老熟妇乱子伦视频在线观看| 久久久久久大精品| 男女那种视频在线观看| 亚洲精品日韩av片在线观看 | 日韩欧美在线二视频| 男女午夜视频在线观看| 1024手机看黄色片| 网址你懂的国产日韩在线| 99国产综合亚洲精品| 国产欧美日韩精品一区二区| 丰满人妻一区二区三区视频av | 国内精品美女久久久久久| 午夜视频国产福利| 1000部很黄的大片| 国产国拍精品亚洲av在线观看 | 国产精品一区二区免费欧美| 亚洲电影在线观看av| or卡值多少钱| 成人一区二区视频在线观看| 他把我摸到了高潮在线观看| e午夜精品久久久久久久| 日韩欧美 国产精品| 少妇裸体淫交视频免费看高清| 夜夜躁狠狠躁天天躁| 啦啦啦免费观看视频1| 国产中年淑女户外野战色| 欧美日韩瑟瑟在线播放| 国产精品精品国产色婷婷| 国产不卡一卡二| 日本一二三区视频观看| 欧美激情久久久久久爽电影| 久久久久久大精品| 夜夜躁狠狠躁天天躁| 久久中文看片网| 校园春色视频在线观看| 午夜免费成人在线视频| 熟女人妻精品中文字幕| 在线a可以看的网站| а√天堂www在线а√下载| 国产亚洲精品综合一区在线观看| 日本精品一区二区三区蜜桃| 丰满人妻熟妇乱又伦精品不卡| 亚洲成人中文字幕在线播放| 美女 人体艺术 gogo| 悠悠久久av| 国产精品嫩草影院av在线观看 | 欧美在线黄色| 国产97色在线日韩免费| 国产精品影院久久| av女优亚洲男人天堂| 成人特级黄色片久久久久久久| 亚洲欧美一区二区三区黑人| 免费av观看视频| 欧美性猛交黑人性爽| 啦啦啦韩国在线观看视频| 亚洲国产欧洲综合997久久,| 亚洲国产精品久久男人天堂| 91av网一区二区| av片东京热男人的天堂| 国产乱人视频| 国产精品美女特级片免费视频播放器| 欧美三级亚洲精品| 欧美xxxx黑人xx丫x性爽| 婷婷精品国产亚洲av| 欧美激情在线99| 日本熟妇午夜| 欧美日韩乱码在线| 国产精品乱码一区二三区的特点| 免费看美女性在线毛片视频| 9191精品国产免费久久| 天天添夜夜摸| netflix在线观看网站| 国产精品爽爽va在线观看网站| 久久久久久久久中文| 欧美性猛交黑人性爽| 亚洲成av人片在线播放无| 中文字幕熟女人妻在线| 真人做人爱边吃奶动态| 麻豆国产97在线/欧美| 一级毛片女人18水好多| 欧美黑人欧美精品刺激| 日本在线视频免费播放| 亚洲美女视频黄频| 欧洲精品卡2卡3卡4卡5卡区| 久久精品国产自在天天线| 99久久综合精品五月天人人| 午夜福利18| 久久精品国产99精品国产亚洲性色| 成人三级黄色视频| www国产在线视频色| 亚洲国产精品sss在线观看| 特大巨黑吊av在线直播| 国产野战对白在线观看| 午夜激情欧美在线| 国产野战对白在线观看|