• <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
    亚洲精品视频女| 女人精品久久久久毛片| 国产国拍精品亚洲av在线观看| 九色成人免费人妻av| 久久ye,这里只有精品| 9色porny在线观看| 日本欧美国产在线视频| 国产精品一二三区在线看| 熟妇人妻不卡中文字幕| 成人午夜精彩视频在线观看| 岛国毛片在线播放| 日韩人妻高清精品专区| 日日撸夜夜添| 久久精品久久精品一区二区三区| 国产国拍精品亚洲av在线观看| 中国国产av一级| av国产精品久久久久影院| 久久 成人 亚洲| 99久久综合免费| 高清欧美精品videossex| 男人爽女人下面视频在线观看| 国产亚洲一区二区精品| 十八禁高潮呻吟视频| 久久婷婷青草| 夫妻性生交免费视频一级片| 超色免费av| 精品人妻偷拍中文字幕| 久久人妻熟女aⅴ| 性高湖久久久久久久久免费观看| 国产免费又黄又爽又色| 国产精品99久久99久久久不卡 | 中文精品一卡2卡3卡4更新| 有码 亚洲区| 男女国产视频网站| 久久99热这里只频精品6学生| 人妻人人澡人人爽人人| 王馨瑶露胸无遮挡在线观看| 黄色一级大片看看| 韩国av在线不卡| 麻豆精品久久久久久蜜桃| 丰满少妇做爰视频| 国产成人精品婷婷| 少妇人妻 视频| 精品久久久久久电影网| 久久精品国产自在天天线| 欧美精品一区二区大全| 久久久精品免费免费高清| 国产精品人妻久久久久久| 亚洲精品中文字幕在线视频| 一级毛片我不卡| 纵有疾风起免费观看全集完整版| 岛国毛片在线播放| 男女无遮挡免费网站观看| 九草在线视频观看| 免费高清在线观看视频在线观看| 亚洲av福利一区| 国产高清国产精品国产三级| 91精品国产国语对白视频| 日韩亚洲欧美综合| 久久精品国产a三级三级三级| 久久久精品免费免费高清| 99精国产麻豆久久婷婷| 亚洲美女黄色视频免费看| 婷婷色麻豆天堂久久| 亚洲国产av影院在线观看| 黄片无遮挡物在线观看| 亚洲第一区二区三区不卡| 亚洲国产精品一区二区三区在线| 爱豆传媒免费全集在线观看| 好男人视频免费观看在线| 久久久久人妻精品一区果冻| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产精品999| 十八禁网站网址无遮挡| 男的添女的下面高潮视频| 少妇丰满av| 国产国拍精品亚洲av在线观看| 91久久精品国产一区二区成人| 最近的中文字幕免费完整| 亚洲不卡免费看| 免费看不卡的av| 美女脱内裤让男人舔精品视频| 如日韩欧美国产精品一区二区三区 | 国产探花极品一区二区| 日韩av不卡免费在线播放| 国产国语露脸激情在线看| 欧美另类一区| 精品少妇内射三级| 自线自在国产av| 免费看不卡的av| 亚洲欧美精品自产自拍| 少妇人妻久久综合中文| av国产精品久久久久影院| 亚洲精品日本国产第一区| 午夜激情福利司机影院| 久久久a久久爽久久v久久| 另类亚洲欧美激情| 亚洲欧美精品自产自拍| 国产一区亚洲一区在线观看| 一区二区三区免费毛片| 免费黄网站久久成人精品| 中国美白少妇内射xxxbb| 免费大片黄手机在线观看| 91久久精品电影网| 免费不卡的大黄色大毛片视频在线观看| 人人妻人人添人人爽欧美一区卜| 欧美精品亚洲一区二区| 极品少妇高潮喷水抽搐| h视频一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 成人国产麻豆网| 精品一区在线观看国产| 在现免费观看毛片| 午夜影院在线不卡| 欧美日韩综合久久久久久| 免费人妻精品一区二区三区视频| 高清av免费在线| 欧美少妇被猛烈插入视频| 精品卡一卡二卡四卡免费| 亚洲五月色婷婷综合| av不卡在线播放| 欧美亚洲日本最大视频资源| 精品人妻在线不人妻| 在线天堂最新版资源| 日韩伦理黄色片| 亚洲精品国产色婷婷电影| 夜夜骑夜夜射夜夜干| 丝袜脚勾引网站| 久久久精品区二区三区| 国产精品免费大片| 久久女婷五月综合色啪小说| 欧美人与善性xxx| 国产探花极品一区二区| 国产白丝娇喘喷水9色精品| 亚洲av二区三区四区| 熟女电影av网| 大话2 男鬼变身卡| 纵有疾风起免费观看全集完整版| 十八禁网站网址无遮挡| 大香蕉久久网| 色94色欧美一区二区| 成人毛片a级毛片在线播放| 亚洲国产精品专区欧美| 午夜91福利影院| 婷婷成人精品国产| 免费看不卡的av| 伊人亚洲综合成人网| 亚洲怡红院男人天堂| 日本av手机在线免费观看| 欧美xxxx性猛交bbbb| 久久久a久久爽久久v久久| www.av在线官网国产| 亚洲性久久影院| 永久网站在线| 亚洲av男天堂| 久久99热这里只频精品6学生| 久久精品夜色国产| 亚洲第一av免费看| 国产精品成人在线| 综合色丁香网| 成人二区视频| 久久毛片免费看一区二区三区| 免费观看在线日韩| 亚洲欧洲国产日韩| 久久99热这里只频精品6学生| 亚洲精品成人av观看孕妇| 成人综合一区亚洲| 国产av精品麻豆| 桃花免费在线播放| 一级毛片我不卡| 亚洲图色成人| 久久国产精品男人的天堂亚洲 | 成年人免费黄色播放视频| 午夜日本视频在线| 青春草亚洲视频在线观看| 久久久久国产精品人妻一区二区| 青春草视频在线免费观看| 成年人免费黄色播放视频| 寂寞人妻少妇视频99o| 久久久久久伊人网av| 王馨瑶露胸无遮挡在线观看| 99热这里只有精品一区| 国产成人精品在线电影| 桃花免费在线播放| 简卡轻食公司| 日韩人妻高清精品专区| 国产精品秋霞免费鲁丝片| 夜夜骑夜夜射夜夜干| 涩涩av久久男人的天堂| 黄色一级大片看看| 亚洲av福利一区| 午夜91福利影院| 久久这里有精品视频免费| 天天操日日干夜夜撸| 亚洲美女黄色视频免费看| 国产一区二区在线观看日韩| 亚洲内射少妇av| 欧美日韩视频高清一区二区三区二| 最后的刺客免费高清国语| 97在线人人人人妻| 国产女主播在线喷水免费视频网站| av又黄又爽大尺度在线免费看| 久久国产亚洲av麻豆专区| 妹子高潮喷水视频| 欧美激情国产日韩精品一区| 久久99精品国语久久久| 中文字幕av电影在线播放| 如何舔出高潮| 性高湖久久久久久久久免费观看| 国产精品一区二区在线不卡| 午夜福利视频在线观看免费| 日本av免费视频播放| 99九九在线精品视频| 亚洲丝袜综合中文字幕| 视频区图区小说| 色婷婷久久久亚洲欧美| 久久久国产精品麻豆| 一级片'在线观看视频| 国产精品一二三区在线看| 色网站视频免费| 黑人猛操日本美女一级片| 午夜精品国产一区二区电影| 日本免费在线观看一区| 自拍欧美九色日韩亚洲蝌蚪91| 五月天丁香电影| 精品视频人人做人人爽| 国产黄色免费在线视频| 亚洲怡红院男人天堂| 国产成人aa在线观看| 日韩在线高清观看一区二区三区| 亚洲美女黄色视频免费看| 高清不卡的av网站| videos熟女内射| 久久国产精品大桥未久av| 3wmmmm亚洲av在线观看| 亚洲精品av麻豆狂野| 狂野欧美激情性bbbbbb| 两个人免费观看高清视频| 亚洲精品久久久久久婷婷小说| 狂野欧美激情性xxxx在线观看| 国产精品一二三区在线看| 精品少妇黑人巨大在线播放| 国产精品秋霞免费鲁丝片| 亚洲综合精品二区| 一级毛片aaaaaa免费看小| 三级国产精品欧美在线观看| av一本久久久久| 亚洲欧美一区二区三区国产| 亚洲在久久综合| 交换朋友夫妻互换小说| 黑人猛操日本美女一级片| 99热这里只有是精品在线观看| 国产精品人妻久久久久久| 日本黄大片高清| 国产一区二区三区av在线| 亚洲色图 男人天堂 中文字幕 | 在线 av 中文字幕| 久久亚洲国产成人精品v| 国产黄频视频在线观看| 91精品伊人久久大香线蕉| 飞空精品影院首页| 国产精品一区二区在线观看99| 女人精品久久久久毛片| 久久青草综合色| 晚上一个人看的免费电影| 精品酒店卫生间| 搡女人真爽免费视频火全软件| 亚洲欧美清纯卡通| 夜夜爽夜夜爽视频| 伊人亚洲综合成人网| 99九九线精品视频在线观看视频| 如何舔出高潮| 欧美日韩综合久久久久久| 五月玫瑰六月丁香| 国产 精品1| 尾随美女入室| 少妇人妻久久综合中文| 久久久久国产精品人妻一区二区| 欧美激情国产日韩精品一区| 女人精品久久久久毛片| 国精品久久久久久国模美| 亚洲精品第二区| 老司机影院成人| 国产一区有黄有色的免费视频| 久久人妻熟女aⅴ| 又粗又硬又长又爽又黄的视频| 日本黄色片子视频| 亚洲国产日韩一区二区| 亚洲精品乱久久久久久| 简卡轻食公司| 日韩成人伦理影院| 亚洲综合色惰| 国产精品一区www在线观看| 人妻一区二区av| 亚洲怡红院男人天堂| 久久国内精品自在自线图片| 亚洲丝袜综合中文字幕| 制服诱惑二区| 欧美 亚洲 国产 日韩一| 男女无遮挡免费网站观看| 亚洲av国产av综合av卡| 亚洲欧美日韩卡通动漫| 午夜影院在线不卡| 中国三级夫妇交换| 韩国av在线不卡| 成人亚洲精品一区在线观看| 久久人人爽人人爽人人片va| 建设人人有责人人尽责人人享有的| 免费大片黄手机在线观看| 国产欧美另类精品又又久久亚洲欧美| 99热全是精品| 少妇的逼水好多| 国产极品天堂在线| 91精品伊人久久大香线蕉| 欧美一级a爱片免费观看看| 大香蕉久久网| 男女国产视频网站| 黑人巨大精品欧美一区二区蜜桃 | 美女中出高潮动态图| 夫妻午夜视频| 下体分泌物呈黄色| 国产成人freesex在线| 国产精品国产三级国产av玫瑰| 欧美精品亚洲一区二区| .国产精品久久| 精品久久蜜臀av无| 天堂俺去俺来也www色官网| 人妻制服诱惑在线中文字幕| 亚州av有码| 我要看黄色一级片免费的| 国产高清三级在线| 亚洲国产精品一区二区三区在线| 校园人妻丝袜中文字幕| 又黄又爽又刺激的免费视频.| 精品酒店卫生间| 亚洲欧洲日产国产| 久热久热在线精品观看| 大话2 男鬼变身卡| 啦啦啦视频在线资源免费观看| 秋霞在线观看毛片| 亚洲国产成人一精品久久久| 久久久精品免费免费高清| 国产欧美日韩综合在线一区二区| 制服人妻中文乱码| 97在线人人人人妻| 18在线观看网站| 亚洲美女视频黄频| 国精品久久久久久国模美| 精品亚洲成国产av| 亚洲av欧美aⅴ国产| 国产成人精品婷婷| 99久久精品国产国产毛片| 尾随美女入室| av又黄又爽大尺度在线免费看| 久久精品夜色国产| 水蜜桃什么品种好| 大又大粗又爽又黄少妇毛片口| 看十八女毛片水多多多| 免费高清在线观看视频在线观看| 久久人人爽av亚洲精品天堂| 99热这里只有精品一区| 观看av在线不卡| 午夜激情av网站| 精品一品国产午夜福利视频| 26uuu在线亚洲综合色| 男人爽女人下面视频在线观看| 国产色爽女视频免费观看| 国产黄频视频在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲人成网站在线播| 精品亚洲乱码少妇综合久久| 免费少妇av软件| 婷婷色麻豆天堂久久| av卡一久久| 久久这里有精品视频免费| 免费观看av网站的网址| 人人妻人人澡人人看| 一本久久精品| 久久久久久人妻| 18禁动态无遮挡网站| 久久女婷五月综合色啪小说| 汤姆久久久久久久影院中文字幕| 日韩精品有码人妻一区| 国产亚洲欧美精品永久| av国产久精品久网站免费入址| 97精品久久久久久久久久精品| 涩涩av久久男人的天堂| 日本黄色日本黄色录像| 成人手机av| √禁漫天堂资源中文www| 五月天丁香电影| 日韩精品免费视频一区二区三区 | 日韩中字成人| 免费观看性生交大片5| 99热全是精品| 午夜视频国产福利| 亚洲美女视频黄频| 蜜桃久久精品国产亚洲av| 亚洲,欧美,日韩| 3wmmmm亚洲av在线观看| 国产一区亚洲一区在线观看| 99久国产av精品国产电影| 国产精品麻豆人妻色哟哟久久| 国产一区二区三区av在线| 人妻人人澡人人爽人人| 色5月婷婷丁香| 国产精品秋霞免费鲁丝片| 国产一区二区在线观看日韩| 熟女人妻精品中文字幕| 欧美+日韩+精品| 日韩大片免费观看网站| 欧美国产精品一级二级三级| 青春草亚洲视频在线观看| 乱码一卡2卡4卡精品| 97超碰精品成人国产| 久久久久精品久久久久真实原创| 狂野欧美白嫩少妇大欣赏| 日韩av在线免费看完整版不卡| 亚洲综合精品二区| 午夜福利视频精品| 极品少妇高潮喷水抽搐| 亚州av有码| 久久久久久人妻| 日韩制服骚丝袜av| 999精品在线视频| 成人手机av| a 毛片基地| 人人妻人人澡人人爽人人夜夜| 啦啦啦视频在线资源免费观看| 国产不卡av网站在线观看| 国产男人的电影天堂91| 亚洲丝袜综合中文字幕| 777米奇影视久久| 婷婷色av中文字幕| 3wmmmm亚洲av在线观看| 超碰97精品在线观看| 欧美xxⅹ黑人| 简卡轻食公司| 国产免费福利视频在线观看| av卡一久久| 最近中文字幕2019免费版| 日本爱情动作片www.在线观看| 日韩成人伦理影院| 人人澡人人妻人| 国产精品.久久久| 黄片播放在线免费| 两个人免费观看高清视频| 伦理电影免费视频| 免费人成在线观看视频色| 3wmmmm亚洲av在线观看| 亚洲不卡免费看| 日本av免费视频播放| 在线 av 中文字幕| 久久久亚洲精品成人影院| 亚洲综合精品二区| 伊人久久精品亚洲午夜| 国产极品粉嫩免费观看在线 | 日韩中文字幕视频在线看片| 婷婷色av中文字幕| 日韩精品有码人妻一区| 99久久精品国产国产毛片| 国产成人精品在线电影| 一区二区三区四区激情视频| 国产av精品麻豆| 91久久精品国产一区二区成人| 午夜视频国产福利| 精品久久国产蜜桃| 蜜桃久久精品国产亚洲av| 桃花免费在线播放| 国产在线一区二区三区精| av又黄又爽大尺度在线免费看| av在线老鸭窝| 国产视频首页在线观看| 人人妻人人爽人人添夜夜欢视频| 自拍欧美九色日韩亚洲蝌蚪91| h视频一区二区三区| 精品久久久久久久久av| 99热这里只有精品一区| 日韩精品有码人妻一区| 97在线人人人人妻| 又大又黄又爽视频免费| 国产在线一区二区三区精| 亚洲美女黄色视频免费看| 亚洲精品456在线播放app| av网站免费在线观看视频| 国产精品久久久久久精品电影小说| 性高湖久久久久久久久免费观看| 久久午夜福利片| 国产成人freesex在线| 精品少妇内射三级| 老司机影院毛片| 亚洲成人手机| 亚洲第一区二区三区不卡| av国产久精品久网站免费入址| 夫妻性生交免费视频一级片| av在线播放精品| 人人妻人人澡人人爽人人夜夜| 黄色怎么调成土黄色| 亚洲中文av在线| 国产精品欧美亚洲77777| 国产精品人妻久久久久久| 大片电影免费在线观看免费| 考比视频在线观看| 高清av免费在线| 午夜免费观看性视频| 久久久久久久久大av| av在线老鸭窝| 日本vs欧美在线观看视频| 国产高清三级在线| 亚洲精品日韩av片在线观看| 久久精品国产鲁丝片午夜精品| av网站免费在线观看视频| 丰满饥渴人妻一区二区三| 哪个播放器可以免费观看大片| 日本黄色片子视频| 久久99一区二区三区| 麻豆乱淫一区二区| 中文字幕人妻丝袜制服| 热re99久久国产66热| 亚洲国产精品999| 国产精品人妻久久久影院| 美女视频免费永久观看网站| 精品一区二区免费观看| 久久久久久久国产电影| 搡女人真爽免费视频火全软件| 亚洲av综合色区一区| 国产男女超爽视频在线观看| 欧美日韩视频高清一区二区三区二| 亚洲美女搞黄在线观看| 全区人妻精品视频| www.色视频.com| 日本猛色少妇xxxxx猛交久久| 麻豆乱淫一区二区| 国产不卡av网站在线观看| 婷婷色综合大香蕉| 亚洲情色 制服丝袜| 在线亚洲精品国产二区图片欧美 | 亚洲久久久国产精品| 只有这里有精品99| 色视频在线一区二区三区| 在线观看www视频免费| 国产午夜精品一二区理论片| 国产永久视频网站| 成人国产av品久久久| 欧美丝袜亚洲另类| 欧美激情国产日韩精品一区| 少妇丰满av| 欧美性感艳星| 亚洲久久久国产精品| 国产亚洲最大av| 黄色一级大片看看| 国产黄片视频在线免费观看| 欧美成人精品欧美一级黄| 国产成人免费观看mmmm| 狠狠精品人妻久久久久久综合| 天堂中文最新版在线下载| 国产视频内射| 这个男人来自地球电影免费观看 | 少妇的逼好多水| 久久精品久久久久久久性| 少妇人妻 视频| 国产亚洲一区二区精品| 亚洲国产av影院在线观看| 国产精品蜜桃在线观看| 爱豆传媒免费全集在线观看| 日本91视频免费播放| 国产一区二区在线观看av| 亚洲丝袜综合中文字幕| 九色成人免费人妻av| 欧美+日韩+精品| 黑人欧美特级aaaaaa片| 大陆偷拍与自拍| 毛片一级片免费看久久久久| 人妻夜夜爽99麻豆av| 亚洲精品日本国产第一区| 亚洲天堂av无毛| 久久午夜福利片| 插逼视频在线观看| 一边摸一边做爽爽视频免费| 在线天堂最新版资源| 久久精品夜色国产| 成人国产av品久久久| 校园人妻丝袜中文字幕| 日日摸夜夜添夜夜添av毛片| 啦啦啦视频在线资源免费观看| 日韩av不卡免费在线播放| 涩涩av久久男人的天堂| xxxhd国产人妻xxx| 两个人的视频大全免费| 乱人伦中国视频| 成人漫画全彩无遮挡| 91精品三级在线观看| 国产精品一区www在线观看| 精品人妻偷拍中文字幕| 成年av动漫网址| 亚洲精品国产色婷婷电影| 亚洲国产精品999| 成年人午夜在线观看视频| 亚洲av日韩在线播放| 欧美 日韩 精品 国产| 26uuu在线亚洲综合色| kizo精华| 一区二区三区四区激情视频| 22中文网久久字幕| 欧美精品亚洲一区二区| 人妻 亚洲 视频| 18禁在线播放成人免费| 久久久欧美国产精品| 国产伦精品一区二区三区视频9| 午夜福利视频在线观看免费| 美女主播在线视频| 久久久精品免费免费高清| 99热网站在线观看| 国产午夜精品一二区理论片| 国内精品宾馆在线|