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

    日徑流序列混沌識(shí)別中合理長度的確定

    2017-02-27 23:25:30周長讓陳元芳顧圣華黃琴
    南水北調(diào)與水利科技 2017年1期
    關(guān)鍵詞:北碚武隆延遲時(shí)間

    周長讓+陳元芳+顧圣華+黃琴

    摘要:混沌理論是進(jìn)行水文時(shí)間序列分析的重要手段。為保證分析結(jié)果的可靠,主張充分利用現(xiàn)有資料,但目前缺乏時(shí)間序列長度對(duì)混沌特性識(shí)別影響的研究。以長江上游武隆站和北碚站日徑流序列(1951年-2012年)為例,通過對(duì)二者進(jìn)行混沌分析,研究了最大Lyapunov指數(shù)對(duì)序列長度的響應(yīng)。結(jié)果表明,日徑流時(shí)間序列長度過小時(shí)會(huì)影響混沌識(shí)別結(jié)果,使結(jié)果缺乏可靠性;并不是樣本序列長度越長混沌識(shí)別結(jié)果越好;當(dāng)序列長度達(dá)到3 000左右時(shí),序列的混沌特性達(dá)到穩(wěn)定,結(jié)果可靠并縮短了計(jì)算時(shí)間。

    關(guān)鍵詞:水文時(shí)間序列分析;日徑流時(shí)間序列;時(shí)間序列長度;混沌識(shí)別;飽和關(guān)聯(lián)維數(shù);最大 Lyapunov指數(shù);混沌特性

    中圖分類號(hào):TV121文獻(xiàn)標(biāo)志碼:A文章編號(hào):1672-1683(2017)01-0055-05

    Abstract:Chaotic theory is an important means of hydrology time series analysis.In order to get reliable analysis results,it is recommended to make a full use of time series.But the research about how the length of time series affects the identification of chaotic characteristics is rare.In this paper,we carried out a study about the responding effect of the maximum Lyapunov exponent to the length of time series with the use of daily runoff time series of gauged stations named Wulong and Beibei in Yangtze River.The result suggested that short daily runoff time series would affect the result of chaotic identification and make the result unreliable;besides,when the length of daily runoff time series reached 3 000,the chaotic characteristics became stable and reliable,and it saved a lot of computing time at the same time.

    Key words:hydrology time series analysis;daily runoff time series;the length of time series;chaotic identification;saturated correlation dimensions;the maximal Lyapunov exponents;chaotic characteristics

    流域的徑流過程的研究方法和手段有很多:基于物理基礎(chǔ)的水文模型、隨機(jī)水文時(shí)間序列分析方法、人工神經(jīng)網(wǎng)絡(luò)、模糊集、模糊邏輯、模糊神經(jīng)、遺傳算法、以及混沌等[1]?;煦绗F(xiàn)象存在于諸多研究領(lǐng)域,包括天文學(xué)、生物學(xué)、化學(xué)、生態(tài)學(xué)、工程學(xué)以及物理學(xué)等[2]?;煦缈茖W(xué)是非線性科學(xué)的一個(gè)重要分支,同時(shí)是解決非線性問題的重要工具。水文現(xiàn)象是一種具有非線性規(guī)律的自然現(xiàn)象,以往的研究結(jié)果表明水文時(shí)間序列中存在混沌現(xiàn)象,因此探究水文序列的混沌特性對(duì)深入了解水文規(guī)律具有重要意義。關(guān)于利用混沌理論進(jìn)行水文時(shí)間序列分析,已經(jīng)有若干進(jìn)展。李國良[4]利用混沌理論,系統(tǒng)而全面地對(duì)時(shí)間序列的混沌分析方法進(jìn)行了介紹。樓玉[5]建立了騰龍橋站日徑流和蘭溪站月徑流的混沌徑流預(yù)報(bào)模型,模型預(yù)報(bào)結(jié)果良好。牟麗琴[6]等則對(duì)汀江流域的月降雨混沌序列進(jìn)行了混沌識(shí)別。袁鵬[7]等證明了四川省6個(gè)水文站的月降雨量時(shí)間序列可能存在混沌現(xiàn)象。周寅康[8]等對(duì)淮河流域的洪澇序列進(jìn)行了混沌識(shí)別,并確定了該時(shí)間序列的最大可預(yù)報(bào)時(shí)間長度。但是針對(duì)時(shí)間序列長度對(duì)序列混沌特性識(shí)別的影響,研究甚少。此外,在進(jìn)行相空間重構(gòu)時(shí),考慮到不同序列長度所需要的計(jì)算時(shí)長不同,序列的混沌特性隨著序列長度的變化的規(guī)律也尚不清晰。因此探究時(shí)間序列長度與最大Lyapunov指數(shù)之間的關(guān)系,對(duì)合理選擇序列長度、混沌識(shí)別可靠性及減少計(jì)算工作量具有重要意義。

    1 混沌識(shí)別與最大Lyapunov指數(shù)的計(jì)算

    1.1 混沌識(shí)別

    混沌是非線性動(dòng)力系統(tǒng)所特有的一種運(yùn)動(dòng)形式,它是既普遍存在又極具復(fù)雜性的現(xiàn)象[9]。混沌現(xiàn)象是貌似無序中的有序,是確定的映射中的類隨機(jī)性[10]。混沌理論在進(jìn)行時(shí)間序列分析時(shí)可分為混沌識(shí)別和混沌預(yù)測(cè)。混沌識(shí)別先對(duì)水文時(shí)間序列進(jìn)行相空間重構(gòu),繼而在高維空間中識(shí)別其混沌特性;混沌預(yù)測(cè)主要是結(jié)合混沌識(shí)別與水文預(yù)測(cè)模型,對(duì)水文序列未來可能的變化進(jìn)行預(yù)測(cè)。

    混沌識(shí)別的主要步驟為:(1)相空間重構(gòu)(自相關(guān)函數(shù)法確定延遲時(shí)間;飽和關(guān)聯(lián)維數(shù)法同時(shí)確定混沌序列的最小嵌入維數(shù)和飽和關(guān)聯(lián)維數(shù));(2)混沌特性的定量識(shí)別(Wolf方法確定序列的最大Lyapunov指數(shù))。

    1.2 方法介紹

    1.2.1 自相關(guān)函數(shù)法確定延遲時(shí)間

    自相關(guān)函數(shù)描述了序列在第i時(shí)刻和第i+τ時(shí)刻,運(yùn)動(dòng)的相似程度。一般而言,時(shí)間間隔τ越小,xi和xi+τ的相似程度越高,r(τ)也就越大。對(duì)于時(shí)間序列{x1,x2,…,xn},自相關(guān)函數(shù)可以表示為:

    實(shí)際工作中,延遲時(shí)間的選取不宜過大和過小[11],因此應(yīng)根據(jù)自相關(guān)函數(shù),選擇合適的延遲時(shí)間τ。一般選擇自相關(guān)函數(shù)首次經(jīng)過零點(diǎn)時(shí)對(duì)應(yīng)的時(shí)間為延遲時(shí)間;但當(dāng)自相關(guān)函數(shù)首次取零時(shí)的延遲時(shí)間過大,則可以減小相空間重構(gòu)的時(shí)間移動(dòng)值,選擇自相關(guān)函數(shù)下降趨于穩(wěn)定時(shí)對(duì)應(yīng)的時(shí)間為延遲時(shí)間。

    1.2.2 Wolf方法計(jì)算最大Lyapunov指數(shù)

    衡量動(dòng)力系統(tǒng)對(duì)初始條件敏感依賴性的一個(gè)重要指標(biāo)就是Lyapunov指數(shù)[12]。Lyapunov指數(shù)是混沌現(xiàn)象的定性和定量的特征參數(shù),它表現(xiàn)為相點(diǎn)在相空間中的相鄰軌道之間以指數(shù)形式快速分離或者收斂[13]。假設(shè)空間中有兩個(gè)相點(diǎn)xi和xj,兩點(diǎn)間的初始距離為d0=‖xi-xj‖,在系統(tǒng)演化τ時(shí)段后,兩點(diǎn)間的距離為dk=‖xi+τ-xj+τ‖,如果dk≈d0eλτ,即兩點(diǎn)間距離呈指數(shù)分離,其中λ就是最大Lyapunov指數(shù)[13]。對(duì)于最大Lyapunov指數(shù)的計(jì)算,常用Wolf方法估計(jì)。計(jì)算公式如下:

    式中:l為向量個(gè)數(shù);τ為延遲時(shí)間;Di為重構(gòu)的相空間里兩個(gè)相點(diǎn)之間的最小歐氏距離;λ為最大Lyapunov指數(shù),其他參數(shù)意義同前。當(dāng)最大Lyapunov指數(shù)大于0時(shí),可以判斷出序列可能為混沌序列,且最大Lyapunov指數(shù)的取值越大,序列的混沌特性越明顯。

    2 實(shí)例研究

    本文以長江流域的武隆站和北碚站歷年逐日徑流時(shí)間序列為例,對(duì)時(shí)間序列長度與水文序列的混沌特性進(jìn)行研究。為了進(jìn)行相似水文背景下的研究,選取兩個(gè)水文站1951年-2012年的日平均流量序列作為研究對(duì)象。

    2.1 自相關(guān)函數(shù)法確定延遲時(shí)間

    根據(jù)以上所介紹的自相關(guān)函數(shù)法,利用已有的日徑流資料繪制出給定時(shí)間移動(dòng)值τ下的自相關(guān)函數(shù)。在實(shí)際計(jì)算和繪圖過程中,為防止延遲時(shí)間選擇的不合理,分別采用τ=1 000與τ=100繪制了兩組自相關(guān)函數(shù)曲線,以便于對(duì)比分析。結(jié)果表明,兩個(gè)水文站在τ=1 000時(shí)自相關(guān)函數(shù)曲線首次經(jīng)過零點(diǎn)時(shí)所對(duì)應(yīng)的延遲時(shí)間均在80~100之間,明顯偏大,因此不予考慮。從而主要取τ=100時(shí)自相關(guān)曲線變化趨于平緩所對(duì)應(yīng)的時(shí)間,作為序列相空間重構(gòu)的延遲時(shí)間。兩個(gè)站日徑流時(shí)間序列的自相關(guān)曲線見圖1、圖2。計(jì)算結(jié)果表明,武隆站日徑流時(shí)間序列的延遲時(shí)間τ0為10,北碚站日徑流時(shí)間序列的延遲時(shí)間τ0為4。

    2.2 飽和關(guān)聯(lián)維數(shù)法確定日徑流時(shí)間序列的嵌入維數(shù)和飽和關(guān)聯(lián)維數(shù)

    實(shí)例研究中,根據(jù)以上所介紹的飽和關(guān)聯(lián)維數(shù)法,取不同的嵌入維數(shù)m對(duì)原序列進(jìn)行相空間重構(gòu),得出對(duì)應(yīng)維數(shù)下的飽和關(guān)聯(lián)維數(shù)D2(m),最后建立嵌入維數(shù)m與飽和關(guān)聯(lián)維數(shù)D2(m)之間的關(guān)系曲線,取曲線達(dá)到穩(wěn)定時(shí)對(duì)應(yīng)的嵌入維數(shù)和飽和關(guān)聯(lián)維數(shù)。圖3和圖4中的(a)圖,自上而下為m取2,3,…,12時(shí)對(duì)應(yīng)的飽和關(guān)聯(lián)維數(shù)曲線,(b)圖為飽和關(guān)聯(lián)維數(shù)隨嵌入維數(shù)變化曲線。

    圖3、圖4對(duì)應(yīng)的計(jì)算結(jié)果見表2。從表2可以得出,兩個(gè)水文站的關(guān)聯(lián)維數(shù)在有限的數(shù)據(jù)長度下,隨著嵌入維數(shù)的升高,均表現(xiàn)出了趨于飽和的現(xiàn)象,因此兩個(gè)水文站的水文時(shí)間序列可能為混沌序列。

    王文等[14]指出,在進(jìn)行混沌分析也不一定需要數(shù)千、數(shù)萬個(gè)點(diǎn)的時(shí)間序列,在維數(shù)不高的情況下,有500個(gè)點(diǎn)左右的長度基本滿足要求。由于研究對(duì)象為日徑流時(shí)間序列,且混沌識(shí)別時(shí)的嵌入維數(shù)均較高,因此序列長度(記為L)應(yīng)盡量長。以往研究已證明水文現(xiàn)象普遍具有混沌特性。但是當(dāng)L在100以內(nèi)時(shí),時(shí)間序列并不表現(xiàn)出混沌特性(如北碚站的最大Lyapunov指數(shù)為-123.8)。因此,使用較少的時(shí)間序列進(jìn)行混沌識(shí)別和預(yù)測(cè)缺乏準(zhǔn)確性、可靠性。

    在實(shí)際計(jì)算過程中,當(dāng)L從100增加到6 000的過程中,計(jì)算時(shí)間由不到1 s增加到20 min以上(1 200 s),計(jì)算時(shí)間隨著計(jì)算長度的增加呈近似指數(shù)的增長。而隨著L的增加,時(shí)間序列的混沌特性逐漸達(dá)到穩(wěn)定,從而選擇合適的L不僅可以保證結(jié)果的可靠性,而且可以大大縮短計(jì)算時(shí)間。

    根據(jù)不同時(shí)間步長所獲得的最大Lyapunov指數(shù)(簡稱Lyapunov),繪制Lyapunov與L相關(guān)曲線圖(見圖5)。其中,圖中的曲線為局部二階多項(xiàng)式回歸所對(duì)應(yīng)的趨勢(shì)線。從圖5中可以看出在L較小時(shí),樣本對(duì)總體的代表性較差,因此曲線波動(dòng)較大。之后,隨著L的增加,武隆站的Lyapunov先是表現(xiàn)出明顯的減小趨勢(shì),隨后趨于穩(wěn)定;北碚站的Lyapunov隨著L增加,曲線始終波動(dòng)較為平緩,最終達(dá)到穩(wěn)定,Lyapunov的取值也趨于穩(wěn)定。對(duì)Lyapunov進(jìn)行滑動(dòng)平均的均值,并繪制相關(guān)曲線,見圖6。

    由圖6(a)圖和(b)可知,當(dāng)L=3 000時(shí),時(shí)間序列的Lyapunov基本達(dá)到穩(wěn)定。

    為了比較時(shí)間序列混沌特性的強(qiáng)弱,對(duì)所選取的兩個(gè)水文站的Lyapunov序列穩(wěn)定前后的統(tǒng)計(jì)特性進(jìn)行了分析,結(jié)果見表2。從表2可以看出武隆站的Lyapunov大于北碚站的,換而言之,武隆站日徑流時(shí)間序列的混沌特性更強(qiáng)。此外,武隆站的Lyapunov序列隨時(shí)間序列長度的波動(dòng)相對(duì)較小,混沌特性變化較小,而北碚站的Lyapunov隨時(shí)間長度的波動(dòng)較大,混沌特性波動(dòng)較小。

    王亦農(nóng)[15]指出“水文水資源系統(tǒng)具有產(chǎn)生混沌的基本條件—對(duì)初始條件的敏感性和內(nèi)在的隨機(jī)性。對(duì)于同一水文系統(tǒng),不同的邊界條件會(huì)產(chǎn)生完全不同的水文過程?!?換而言之,對(duì)于同一場(chǎng)暴雨過程,由于兩個(gè)站點(diǎn)的位置不同,二者的水文過程也大相徑庭。借此,對(duì)以上結(jié)果(武隆站的Lyapunov指數(shù)大于北碚站)進(jìn)行分析可知,武隆站對(duì)初值的敏感性更強(qiáng),內(nèi)在的隨機(jī)性也更強(qiáng)。影響時(shí)間序列混沌特性的因素眾多,若從水文站點(diǎn)所在流域特征來說,嘉陵江流域的集水面積相對(duì)烏江更大,調(diào)蓄能力較強(qiáng),其對(duì)初始條件的敏感性相對(duì)于烏江更弱。從而武隆站日徑流時(shí)間序列的Lyapunov指數(shù)較大,因此混沌特性相對(duì)較強(qiáng),而北碚站則較弱。

    4 結(jié)論

    通過對(duì)實(shí)例中的武隆站和北碚站1951年-2012年日平均流量資料序列進(jìn)行混沌分析和相應(yīng)的計(jì)算結(jié)果可以得到如下結(jié)論。

    (1)從混沌特征角度。(a)飽和關(guān)聯(lián)維數(shù)隨著嵌入維數(shù)的升高呈現(xiàn)出趨于飽和的現(xiàn)象;(b)最大Lyapunov指數(shù)均大于0,體現(xiàn)了序列對(duì)初值具有很強(qiáng)的敏感性。因此這兩個(gè)序列具有混沌序列的基本特征,從而兩個(gè)時(shí)間序列可能為混沌序列。

    (2)從最大Lyapunov指數(shù)的統(tǒng)計(jì)特性角度。(a)武隆站日徑流時(shí)間序列的最大Lyapunov指數(shù)較大,因此混沌特性相對(duì)較強(qiáng),而北碚站則較弱;(b)武隆站日徑流時(shí)間序列的混沌特性隨時(shí)間序列長度的增加波動(dòng)較大,北碚站的變化較??;

    (3)從序列長度角度。日徑流時(shí)間序列長度過小時(shí)樣本的混沌分析結(jié)果的可靠性降低,從而在資料長度允許條件下,盡量使用較長的時(shí)間序列;對(duì)于長度較長的序列,進(jìn)行混沌分析時(shí)并非長度越長越好,資料長度過長往往計(jì)算需要較長的計(jì)算時(shí)間,甚至超出計(jì)算內(nèi)存。從實(shí)例研究的結(jié)果來看,當(dāng)序列長度達(dá)到3 000左右時(shí),序列的混沌特性達(dá)到穩(wěn)定,且結(jié)果可靠、計(jì)算時(shí)間相對(duì)較少。從而,選擇合適的序列長度可以在減少工作量的同時(shí),得到代表性較好而穩(wěn)定的結(jié)果。

    參考文獻(xiàn)(References):

    [1] Rahman Khatibi,Bellie Sivakumar,Mohammad Ali Ghorbani,et al.Investigating chaos in river stage and discharge time series[J].Journal of Hydrology,2012,414-415:108-117.

    [2] W W Ng,U S Panu,W C Lennox.Chaos based analytical techniques for daily extreme hydrological observations[J].Journal of Hydrology,2007,342:17-41.

    [3] 李國良.混沌理論及其在水文時(shí)間序列中的應(yīng)用研究[D].哈爾濱:東北農(nóng)業(yè)大學(xué),2007.12.(LI Guo-liang.Studies on chaos theory and its application in hydrology time series[D].Harbin :Northeast Agriculture University,2007.12.(in Chinese)

    [4] 樓玉.混沌時(shí)間序列方法在徑流預(yù)報(bào)中的應(yīng)用研究[D].杭州:浙江大學(xué),2005.8.18.(LOU Yu.Chaotic time series methods in the study of the application of the runoff forecast[D].Zhejiang Hangzhou Zhejiang University,2005.8.18.(in Chinese)

    [5] 牟麗琴,林煥新.汀江流域降水時(shí)間序列的混沌分析[J].人民珠江,2010(01):27-29,55.(MOU Li-qin,LI Huan-xin.Chaotic analysis of rainfall time series in Tingjiang catchment[J].Pearl River,2010(01):27-29,55.(in Chinese)

    [6] 袁鵬,李謂新,王文圣,等.月降雨量時(shí)間序列中的混沌現(xiàn)象[J].四川大學(xué)學(xué)報(bào):工程科學(xué)版,2002,34(1):16-19.(YUAN Peng,LI Wei-xin,WANG Wen-sheng,et al.Chaos analysis of monthly precipitation time series[J].Journal of Sichuan University:Engineering Science Edition,2002,34(1):16-19.(in Chinese)

    [7] 周寅康,王臘春,張捷.淮河流域洪澇變化的混沌特征[J].自然災(zāi)害學(xué)報(bào),1999,8(1):42-47.(ZHOU Yin-kang,WANG La-chun,ZHANG Jie.Some chaotic characteristics of the flood series in Huai River basin[J].Journal of Natural Disasters,1999,8(1),42-47.(in Chinese)

    [8] 李新杰.河川徑流時(shí)間序列的非線性特征識(shí)別與分析[D].武漢:武漢大學(xué),2013.3.(LI Xin-jie.Nonlinear characteristics identification and analysis for river runoff time series[D].Wuhan:Wuhan University ,2013.3 .(in Chinese)

    [9] 陳關(guān)榮,呂金虎.Lorenz系統(tǒng)族的動(dòng)力學(xué)分析、控制與同步[M].北京:科學(xué)出版社,2003.(CHEN Guan-rong,LYU Jin-hu.Dynamics analysis、control and synchronization of Lorenz system[M].Beijing:Science Press,2003.

    [10] 王文,許武成.混沌時(shí)間序列可預(yù)報(bào)時(shí)間長度分析[J].河海大學(xué)學(xué)報(bào):自然科學(xué)版,2004,32(4):367-371.(WANG Wen,XU Wu-cheng.Analysis of length of predictable time of chaotic time series[J].Journal of Hohai University :Natural Sciences,2004,32(4):367-371.(in Chinese)

    [11] 袁鵬,李謂新,王文圣,等.月降雨量時(shí)間序列中的混沌現(xiàn)象[J].四川大學(xué)學(xué)報(bào):工程科學(xué)版,2002,34(1):16-19.(YUAN Peng,LI Wei-xin,WANG Wen-sheng,et al.Chaos analysis of monthly precipitation time series[J].Journal of Sichuan University :Engineering Science Edition,2002,34(1),16-19.(in Chinese)

    [12] Alan Wolf,Jack B Swift,Harry L Swinney,et al.Determining Lyapunov exponents from a time series[J].Physica 16D,1985,285-317.

    [13] Lorenz,Edward N.Dimension of weather and climate attractors[J].Nature,1991,353(6341):241-244.

    [14] 王文,許武成.對(duì)水文時(shí)間序列混沌特征參數(shù)估計(jì)問題的討論[J].水科學(xué)進(jìn)展,2005,16(4):609-616.(WANG Wen,XU Wu-cheng.Some issues on the characteristics of chaotic properties of hydrologic time series[J].Advances in Water Science,2005,16(4):609-616.(in Chinese)

    [15] 王亦農(nóng).淺談混沌與水文[J].資源環(huán)境,2011,3(40):71-72.(WANG Yi-nong.A brief talk about chaotic and hydrology[J].Resources and Environment,2011,3(40):71-72.(in Chinese)

    猜你喜歡
    北碚武隆延遲時(shí)間
    “大美武隆”詩詞選登
    中華詩詞(2022年12期)2022-07-31 02:45:16
    北碚:讓綠色成為高質(zhì)量發(fā)展最亮底色
    抗美援朝運(yùn)動(dòng)中的武隆人
    紅巖春秋(2022年1期)2022-04-12 00:37:34
    二氧化碳對(duì)乙烷燃燒著火延遲時(shí)間的影響
    煤氣與熱力(2021年3期)2021-06-09 06:16:22
    中國新詩創(chuàng)研中心在重慶北碚揭牌
    鴨綠江(2021年35期)2021-04-19 12:24:26
    重慶武?。贺毨丝谝劳新糜螛I(yè)脫貧
    LTE 系統(tǒng)下行鏈路FDRX 節(jié)能機(jī)制研究
    基于分層COX模型的跟馳反應(yīng)延遲時(shí)間生存分析
    梁漱溟在重慶北碚的文化教育活動(dòng)
    延遲時(shí)間對(duì)氣輔注射成型氣體穿透行為影響的數(shù)值模擬和實(shí)驗(yàn)研究
    中國塑料(2016年8期)2016-06-27 06:35:02
    18禁美女被吸乳视频| 12—13女人毛片做爰片一| 久久久水蜜桃国产精品网| 少妇的丰满在线观看| 亚洲中文字幕一区二区三区有码在线看 | 麻豆久久精品国产亚洲av| 黑人巨大精品欧美一区二区蜜桃| 国产精品精品国产色婷婷| 亚洲人成电影免费在线| 法律面前人人平等表现在哪些方面| 757午夜福利合集在线观看| 国产麻豆成人av免费视频| 久久久久久人人人人人| 亚洲欧美激情综合另类| 一级a爱片免费观看的视频| av有码第一页| 在线观看免费午夜福利视频| av天堂久久9| 18禁美女被吸乳视频| 欧美绝顶高潮抽搐喷水| 19禁男女啪啪无遮挡网站| 变态另类成人亚洲欧美熟女 | 久久久久亚洲av毛片大全| 欧美在线一区亚洲| 天天一区二区日本电影三级 | 伊人久久大香线蕉亚洲五| 18禁黄网站禁片午夜丰满| 美女午夜性视频免费| 正在播放国产对白刺激| 久久久久久久午夜电影| 欧美性长视频在线观看| 国产亚洲av嫩草精品影院| 国产成人系列免费观看| 久久午夜亚洲精品久久| 999久久久国产精品视频| 久久久久精品国产欧美久久久| 久久香蕉国产精品| 免费无遮挡裸体视频| 国产av精品麻豆| 一级黄色大片毛片| 狂野欧美激情性xxxx| 黄色视频不卡| av免费在线观看网站| 国内毛片毛片毛片毛片毛片| 男男h啪啪无遮挡| 亚洲国产精品sss在线观看| 欧美激情久久久久久爽电影 | 无遮挡黄片免费观看| 国产极品粉嫩免费观看在线| 成人免费观看视频高清| 少妇被粗大的猛进出69影院| 成人永久免费在线观看视频| 夜夜看夜夜爽夜夜摸| 他把我摸到了高潮在线观看| 日日夜夜操网爽| 久9热在线精品视频| 国产精品久久久久久精品电影 | 国产日韩一区二区三区精品不卡| or卡值多少钱| 777久久人妻少妇嫩草av网站| e午夜精品久久久久久久| 岛国视频午夜一区免费看| 精品少妇一区二区三区视频日本电影| 性少妇av在线| 啪啪无遮挡十八禁网站| 男人舔女人的私密视频| 亚洲精品国产色婷婷电影| 亚洲精品美女久久av网站| 国产亚洲精品av在线| 嫩草影院精品99| 国产午夜福利久久久久久| 久久久久国产一级毛片高清牌| 国产成人精品久久二区二区91| 看片在线看免费视频| 国产亚洲精品久久久久5区| 亚洲欧美激情在线| 一区在线观看完整版| 亚洲一区二区三区色噜噜| 久久香蕉激情| 波多野结衣av一区二区av| 岛国在线观看网站| www.熟女人妻精品国产| 1024香蕉在线观看| 午夜成年电影在线免费观看| 欧美日韩乱码在线| 国产精品99久久99久久久不卡| 亚洲精品国产色婷婷电影| 国产熟女午夜一区二区三区| or卡值多少钱| 国产精品九九99| 一区二区日韩欧美中文字幕| 老汉色∧v一级毛片| 好男人电影高清在线观看| 午夜福利一区二区在线看| 国产成人精品在线电影| 夜夜爽天天搞| 三级毛片av免费| 我的亚洲天堂| 欧美绝顶高潮抽搐喷水| 精品福利观看| 亚洲专区国产一区二区| 一边摸一边做爽爽视频免费| 亚洲最大成人中文| 免费一级毛片在线播放高清视频 | 国产成人影院久久av| 黄色视频不卡| 亚洲精品国产区一区二| 女性被躁到高潮视频| 99精品在免费线老司机午夜| 国产激情欧美一区二区| 国产成人av激情在线播放| 中文字幕久久专区| 亚洲欧美日韩高清在线视频| 久久人人爽av亚洲精品天堂| av片东京热男人的天堂| 午夜老司机福利片| 亚洲五月天丁香| 亚洲av日韩精品久久久久久密| 久久久水蜜桃国产精品网| 午夜视频精品福利| 精品国产乱子伦一区二区三区| 国产区一区二久久| 成人手机av| 青草久久国产| av视频在线观看入口| 免费看美女性在线毛片视频| 天天躁夜夜躁狠狠躁躁| 亚洲av成人不卡在线观看播放网| 国产精品亚洲一级av第二区| 在线十欧美十亚洲十日本专区| 大陆偷拍与自拍| 国产一卡二卡三卡精品| 大型黄色视频在线免费观看| 乱人伦中国视频| 国产成人系列免费观看| 精品电影一区二区在线| 日本免费一区二区三区高清不卡 | 韩国av一区二区三区四区| 一级毛片女人18水好多| 最近最新中文字幕大全电影3 | 精品国产乱码久久久久久男人| 中文字幕人妻丝袜一区二区| 日韩精品青青久久久久久| 亚洲欧美日韩无卡精品| 午夜福利一区二区在线看| 亚洲精品在线观看二区| 国产一区二区三区视频了| 精品国产亚洲在线| 色老头精品视频在线观看| 亚洲一区中文字幕在线| 国产亚洲精品久久久久久毛片| 免费在线观看影片大全网站| 少妇粗大呻吟视频| 国产精品亚洲一级av第二区| 老熟妇仑乱视频hdxx| 亚洲一区高清亚洲精品| 亚洲免费av在线视频| 免费在线观看黄色视频的| 男人的好看免费观看在线视频 | 久久国产精品男人的天堂亚洲| 麻豆久久精品国产亚洲av| 国产精品综合久久久久久久免费 | 女警被强在线播放| 99久久综合精品五月天人人| 国产麻豆69| 99精品久久久久人妻精品| 国产精品一区二区免费欧美| 在线观看免费视频日本深夜| 亚洲精品久久成人aⅴ小说| 变态另类丝袜制服| 精品欧美国产一区二区三| 正在播放国产对白刺激| 国产精品一区二区精品视频观看| 亚洲少妇的诱惑av| 国产精品爽爽va在线观看网站 | 麻豆一二三区av精品| 欧美亚洲日本最大视频资源| 国产欧美日韩精品亚洲av| 嫁个100分男人电影在线观看| 一卡2卡三卡四卡精品乱码亚洲| www.熟女人妻精品国产| 精品不卡国产一区二区三区| 午夜激情av网站| 黄色 视频免费看| 久热这里只有精品99| 淫秽高清视频在线观看| 国产国语露脸激情在线看| 久久久久精品国产欧美久久久| 亚洲va日本ⅴa欧美va伊人久久| ponron亚洲| 男女床上黄色一级片免费看| 色老头精品视频在线观看| 久久人人精品亚洲av| 一级片免费观看大全| 欧美最黄视频在线播放免费| 亚洲自拍偷在线| 日韩大码丰满熟妇| 成人永久免费在线观看视频| 最好的美女福利视频网| 欧美老熟妇乱子伦牲交| 怎么达到女性高潮| 咕卡用的链子| 在线观看www视频免费| 免费av毛片视频| 欧美成人一区二区免费高清观看 | 多毛熟女@视频| 侵犯人妻中文字幕一二三四区| 在线观看一区二区三区| 久久精品国产亚洲av高清一级| 中文字幕色久视频| 亚洲自偷自拍图片 自拍| 欧美日本亚洲视频在线播放| 男女之事视频高清在线观看| 久久草成人影院| 激情视频va一区二区三区| 亚洲性夜色夜夜综合| 国产aⅴ精品一区二区三区波| 最好的美女福利视频网| 两性夫妻黄色片| 成人永久免费在线观看视频| 啦啦啦 在线观看视频| 99久久综合精品五月天人人| 亚洲伊人色综图| 1024视频免费在线观看| 757午夜福利合集在线观看| 啦啦啦 在线观看视频| 欧美另类亚洲清纯唯美| 精品久久久久久,| 高清黄色对白视频在线免费看| svipshipincom国产片| 18美女黄网站色大片免费观看| 国产精品秋霞免费鲁丝片| 可以在线观看毛片的网站| 12—13女人毛片做爰片一| 波多野结衣高清无吗| 一进一出好大好爽视频| 啪啪无遮挡十八禁网站| 精品久久蜜臀av无| 国产亚洲精品av在线| 波多野结衣av一区二区av| 日本a在线网址| 免费久久久久久久精品成人欧美视频| 国语自产精品视频在线第100页| 男女下面进入的视频免费午夜 | 亚洲精品久久成人aⅴ小说| 18美女黄网站色大片免费观看| 久久热在线av| 国产高清视频在线播放一区| 国产精品自产拍在线观看55亚洲| 99久久综合精品五月天人人| 女人高潮潮喷娇喘18禁视频| 国产精品美女特级片免费视频播放器 | 人人妻人人爽人人添夜夜欢视频| 高清在线国产一区| 国产一区在线观看成人免费| 国产欧美日韩综合在线一区二区| 色哟哟哟哟哟哟| 大型av网站在线播放| 波多野结衣巨乳人妻| 欧美午夜高清在线| 欧美日韩亚洲国产一区二区在线观看| 在线视频色国产色| 久99久视频精品免费| 亚洲色图av天堂| 一个人免费在线观看的高清视频| 免费在线观看日本一区| 国产av在哪里看| 男女下面进入的视频免费午夜 | 桃色一区二区三区在线观看| 日日摸夜夜添夜夜添小说| 乱人伦中国视频| 欧美激情高清一区二区三区| 大码成人一级视频| 国产精品一区二区精品视频观看| 亚洲av美国av| 在线观看免费午夜福利视频| 成人永久免费在线观看视频| avwww免费| 两性夫妻黄色片| 免费一级毛片在线播放高清视频 | www国产在线视频色| 97人妻天天添夜夜摸| 宅男免费午夜| 一进一出抽搐动态| 久久久久久人人人人人| 久久久久九九精品影院| 黑人巨大精品欧美一区二区mp4| 亚洲七黄色美女视频| 真人做人爱边吃奶动态| 亚洲中文av在线| 亚洲欧洲精品一区二区精品久久久| 亚洲自拍偷在线| 午夜老司机福利片| 韩国精品一区二区三区| 午夜老司机福利片| 日韩精品青青久久久久久| 好男人在线观看高清免费视频 | 国产免费av片在线观看野外av| 久久久国产成人精品二区| 国产成人欧美| 麻豆国产av国片精品| 9热在线视频观看99| av在线天堂中文字幕| 亚洲 欧美一区二区三区| 欧美午夜高清在线| 欧美国产日韩亚洲一区| 精品一品国产午夜福利视频| 免费久久久久久久精品成人欧美视频| 一a级毛片在线观看| 日本撒尿小便嘘嘘汇集6| 亚洲精品中文字幕在线视频| 亚洲人成电影观看| 夜夜看夜夜爽夜夜摸| 国产高清有码在线观看视频 | 欧美人与性动交α欧美精品济南到| 国产乱人伦免费视频| 搡老妇女老女人老熟妇| 伦理电影免费视频| 精品人妻1区二区| 可以在线观看的亚洲视频| 欧美不卡视频在线免费观看 | 日韩成人在线观看一区二区三区| 一区在线观看完整版| 精品国产超薄肉色丝袜足j| 又大又爽又粗| 亚洲欧美日韩另类电影网站| 亚洲人成电影观看| 日本黄色视频三级网站网址| 亚洲成av人片免费观看| 亚洲国产欧美日韩在线播放| 一区在线观看完整版| 精品国产亚洲在线| 男人操女人黄网站| 久久久久久人人人人人| 三级毛片av免费| 午夜福利高清视频| 国产极品粉嫩免费观看在线| 亚洲国产高清在线一区二区三 | 好男人在线观看高清免费视频 | 精品久久久久久成人av| 亚洲欧美激情在线| 99在线人妻在线中文字幕| 久久久国产成人精品二区| 亚洲 欧美 日韩 在线 免费| 女性生殖器流出的白浆| 国产成人啪精品午夜网站| 12—13女人毛片做爰片一| 日日干狠狠操夜夜爽| 久久久久久久久免费视频了| 国产一区在线观看成人免费| 欧美日韩福利视频一区二区| 老熟妇乱子伦视频在线观看| 国产精品日韩av在线免费观看 | 久热这里只有精品99| 免费在线观看亚洲国产| 国产视频一区二区在线看| 丝袜美腿诱惑在线| 最近最新中文字幕大全电影3 | 精品一区二区三区视频在线观看免费| 欧美成狂野欧美在线观看| 久久香蕉精品热| 高清毛片免费观看视频网站| 日本黄色视频三级网站网址| 国产熟女午夜一区二区三区| 好男人电影高清在线观看| 十八禁人妻一区二区| 一本综合久久免费| 国产欧美日韩精品亚洲av| 中文字幕久久专区| 久久精品人人爽人人爽视色| 亚洲av成人一区二区三| 亚洲av成人av| 久久 成人 亚洲| 久久久国产精品麻豆| 国产成人精品久久二区二区免费| 国内毛片毛片毛片毛片毛片| 久久久久国产精品人妻aⅴ院| 大陆偷拍与自拍| 亚洲一区二区三区不卡视频| av视频在线观看入口| 曰老女人黄片| 在线播放国产精品三级| 国产又爽黄色视频| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久久久久精品电影 | 777久久人妻少妇嫩草av网站| 老鸭窝网址在线观看| 国产亚洲精品一区二区www| 免费观看人在逋| 亚洲最大成人中文| 午夜福利在线观看吧| 久久香蕉精品热| 成年人黄色毛片网站| 久久婷婷成人综合色麻豆| 精品熟女少妇八av免费久了| 日韩欧美国产一区二区入口| aaaaa片日本免费| 村上凉子中文字幕在线| 精品午夜福利视频在线观看一区| 国产蜜桃级精品一区二区三区| 日韩欧美国产在线观看| 欧美日本中文国产一区发布| 欧美日本视频| 日本a在线网址| 一边摸一边抽搐一进一出视频| 亚洲五月天丁香| 中文亚洲av片在线观看爽| 又大又爽又粗| 丝袜在线中文字幕| 国产麻豆69| 村上凉子中文字幕在线| www.精华液| 宅男免费午夜| 国产激情欧美一区二区| 99国产综合亚洲精品| 又黄又粗又硬又大视频| 色精品久久人妻99蜜桃| 亚洲中文字幕一区二区三区有码在线看 | 国产伦一二天堂av在线观看| 人成视频在线观看免费观看| 亚洲欧美激情在线| 日本一区二区免费在线视频| 国产男靠女视频免费网站| 亚洲免费av在线视频| 久久国产精品男人的天堂亚洲| √禁漫天堂资源中文www| 亚洲第一青青草原| 91成年电影在线观看| 亚洲中文av在线| 天堂影院成人在线观看| av欧美777| 久久草成人影院| 人人妻人人澡欧美一区二区 | 黑人操中国人逼视频| 亚洲精品中文字幕一二三四区| 国产精品国产高清国产av| 国产精品久久视频播放| 丝袜在线中文字幕| 如日韩欧美国产精品一区二区三区| 纯流量卡能插随身wifi吗| 日韩一卡2卡3卡4卡2021年| 国产野战对白在线观看| √禁漫天堂资源中文www| 欧美老熟妇乱子伦牲交| 精品高清国产在线一区| 男女之事视频高清在线观看| 国产精品1区2区在线观看.| 在线观看免费日韩欧美大片| 波多野结衣高清无吗| 久久草成人影院| 色综合亚洲欧美另类图片| 怎么达到女性高潮| 日韩欧美免费精品| 黑人巨大精品欧美一区二区mp4| 麻豆国产av国片精品| 日本五十路高清| 美女扒开内裤让男人捅视频| 91精品三级在线观看| 9191精品国产免费久久| 国产av精品麻豆| 757午夜福利合集在线观看| 亚洲人成电影观看| 波多野结衣高清无吗| 久久精品国产99精品国产亚洲性色 | 精品人妻在线不人妻| 午夜久久久在线观看| 欧美激情 高清一区二区三区| 人人澡人人妻人| 亚洲人成77777在线视频| 亚洲va日本ⅴa欧美va伊人久久| 视频在线观看一区二区三区| 国产高清视频在线播放一区| 夜夜夜夜夜久久久久| 久久婷婷人人爽人人干人人爱 | 国产精品精品国产色婷婷| 亚洲va日本ⅴa欧美va伊人久久| 久久久久国内视频| 日韩免费av在线播放| 亚洲第一欧美日韩一区二区三区| 免费在线观看完整版高清| 色综合婷婷激情| 人成视频在线观看免费观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲一区高清亚洲精品| 18禁国产床啪视频网站| 91在线观看av| 又大又爽又粗| 婷婷精品国产亚洲av在线| 欧洲精品卡2卡3卡4卡5卡区| 一二三四社区在线视频社区8| 韩国av一区二区三区四区| av超薄肉色丝袜交足视频| 色哟哟哟哟哟哟| 国产精品亚洲av一区麻豆| 午夜免费观看网址| 亚洲午夜理论影院| 视频在线观看一区二区三区| 亚洲av成人不卡在线观看播放网| 日本三级黄在线观看| 欧美激情极品国产一区二区三区| 99久久99久久久精品蜜桃| 亚洲性夜色夜夜综合| 亚洲在线自拍视频| 免费在线观看视频国产中文字幕亚洲| 国产精品九九99| 啦啦啦韩国在线观看视频| 中文字幕人成人乱码亚洲影| 欧美日韩一级在线毛片| 成年版毛片免费区| 国产精品日韩av在线免费观看 | 伊人久久大香线蕉亚洲五| 操出白浆在线播放| 一区二区三区高清视频在线| a在线观看视频网站| 亚洲电影在线观看av| 久久久久九九精品影院| 一级a爱视频在线免费观看| 一级片免费观看大全| 黄色成人免费大全| 亚洲五月天丁香| 麻豆成人av在线观看| 午夜免费激情av| 亚洲中文字幕日韩| 欧美成狂野欧美在线观看| 国产乱人伦免费视频| 麻豆av在线久日| 夜夜看夜夜爽夜夜摸| 日韩有码中文字幕| 一区二区日韩欧美中文字幕| 热re99久久国产66热| 久久天躁狠狠躁夜夜2o2o| 18禁观看日本| 高清在线国产一区| 国产av精品麻豆| 亚洲欧美激情在线| 首页视频小说图片口味搜索| 在线免费观看的www视频| 一本综合久久免费| 大香蕉久久成人网| 村上凉子中文字幕在线| 欧美一级a爱片免费观看看 | 露出奶头的视频| 一级黄色大片毛片| 国产区一区二久久| 91精品三级在线观看| 欧美精品啪啪一区二区三区| 久久香蕉激情| 婷婷六月久久综合丁香| 一级毛片高清免费大全| 亚洲aⅴ乱码一区二区在线播放 | 女同久久另类99精品国产91| 国产成人av教育| 久热这里只有精品99| 亚洲国产欧美网| 十八禁人妻一区二区| 午夜影院日韩av| 午夜免费鲁丝| 久久久久久久精品吃奶| 伊人久久大香线蕉亚洲五| 国产精品永久免费网站| 免费高清在线观看日韩| 色综合站精品国产| 一边摸一边做爽爽视频免费| 十分钟在线观看高清视频www| 亚洲成av片中文字幕在线观看| 在线十欧美十亚洲十日本专区| 69精品国产乱码久久久| 一区二区日韩欧美中文字幕| 黑人操中国人逼视频| 成人手机av| 国产xxxxx性猛交| 久久人妻av系列| 国产精品av久久久久免费| 夜夜看夜夜爽夜夜摸| 91麻豆av在线| 色综合站精品国产| 久久香蕉精品热| 可以免费在线观看a视频的电影网站| 老汉色av国产亚洲站长工具| 欧美激情久久久久久爽电影 | 一区福利在线观看| 男男h啪啪无遮挡| 桃红色精品国产亚洲av| 国产真人三级小视频在线观看| 99精品在免费线老司机午夜| 色精品久久人妻99蜜桃| 亚洲成人精品中文字幕电影| 亚洲色图综合在线观看| 亚洲国产欧美一区二区综合| 国产亚洲精品第一综合不卡| 91老司机精品| 一级片免费观看大全| 午夜影院日韩av| 啦啦啦韩国在线观看视频| 1024视频免费在线观看| 国产成人系列免费观看| 久久 成人 亚洲| 老鸭窝网址在线观看| 女人高潮潮喷娇喘18禁视频| 人妻丰满熟妇av一区二区三区| 老鸭窝网址在线观看| 午夜影院日韩av| 亚洲精品av麻豆狂野| 18禁裸乳无遮挡免费网站照片 | 亚洲av美国av| 成年版毛片免费区| 最新在线观看一区二区三区| 黄色视频不卡| 中文亚洲av片在线观看爽| 成人亚洲精品av一区二区| 国产精品自产拍在线观看55亚洲| 久久婷婷成人综合色麻豆| 在线播放国产精品三级| 一夜夜www| 大型黄色视频在线免费观看|