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

    突發(fā)事件對我國鐵路客運量的干預特征研究

    2022-02-18 14:03:24汪志紅
    鐵道運輸與經濟 2022年1期
    關鍵詞:客運量離群傳遞函數(shù)

    汪志紅

    (廣東金融學院 金融數(shù)據(jù)與統(tǒng)計學院,廣東 廣州 510521)

    0 引言

    近年來,SARS、新冠肺炎疫情和各種自然災害等突發(fā)事件對人們生活和社會經濟造成不同程度的沖擊,階段性地影響著人們日常生活和社會發(fā)展。鐵路作為我國重要交通工具,其客運量變化能夠較為直觀地體現(xiàn)某些突發(fā)事件對人民生活的影響,研究分析其沖擊程度和發(fā)展趨勢有利于更好地預防和應對將來可能引起的影響,這也是各級應急管理機構重點關注的問題之一。

    突發(fā)事件對鐵路客運量的影響實際上是一個突發(fā)事件的干預問題,這種干預可以運用干預模型來進行定量分析。學者們通常從2 個方面基于時間序列構建突發(fā)事件干預模型并進行干預分析:一類是將某些影響因素和政策措施作為時間序列的干預因素[1-2],分析其干預常規(guī)變量對社會經濟影響程度的大小,對其干預趨勢進行研究。該類方法將鐵路客運量分解為ARIMA 模型擬合值與突發(fā)事件干預因素,一定程度上忽視了干預因素與ARIMA 模型內部可能存在的關系。另一類是用脈沖函數(shù)或階梯函數(shù)將某一突發(fā)事件作為干預因素,研究突發(fā)事件或政策干預的程度大小,同時研究其干預趨勢,其研究過程為在識別時間序列異常值和最優(yōu)ARIMA 模型基礎上,對異常值進行基于脈沖函數(shù)或階梯函數(shù)的虛擬變量針對性設計,確定干預過程的傳遞函數(shù),與最優(yōu)ARIMA 模型結合,構建時間序列的干預模型,這是當前國外時間序列干預研究的主流方法[3-5],如Habib[5]等以現(xiàn)代避孕措施的三階段作為階段干預因素,對加納和坦桑尼亞的公共衛(wèi)生事件進行試驗研究,得出帶控制組的時間序列干預模型能夠較為全面地分析公共衛(wèi)生事件干預的變化過程。在該類方法的實際應用中,學者們更多地將干預因素作為時間序列干預模型整體的一部分進行研究,相對第一類方法更能體現(xiàn)模型設計的系統(tǒng)性,但對干預因素類型選擇缺乏理論依據(jù),同時缺乏對趨勢干預過程的比較分析。時間序列模型在鐵路客運量實際應用研究方面,X-13ARIMASEATS 或ARIMA 等時間序列模型更多地應用于鐵路客運量的預測過程[6-9],如錢名軍等[6]運用SARIMA-GARCH 融合模型、SARIMA、ARIMA和NAR 動態(tài)神經網絡模型對鐵路月度客運量進行預測比較研究,認為SARIMA-GARCH 融合模型進行鐵路客運量預測精度更高,但基于時間序列干預模型對鐵路客運量的突發(fā)事件干預分析略顯不足,可以在干預分析的基礎上進一步細化鐵路客運量預測過程,提高預測精度。

    在此,針對突發(fā)事件對鐵路客運量干預分析過程的不完善問題,應用X-13ARIMA-SEATS 季節(jié)調整方法中的自動檢測突發(fā)事件離群值(異常值)功能,將離群值與脈沖函數(shù)或階梯函數(shù)進行特征比較分析,科學確定突發(fā)事件虛擬變量類型,運用設計的虛擬變量在確定突發(fā)事件影響程度大小的同時,基于ARMA 干預傳遞函數(shù)分析突發(fā)事件對鐵路客運量的沖擊趨勢,在總體上提升時間序列干預分析的科學性和系統(tǒng)性的同時,拓展時間序列干預分析在工程上的應用。

    1 基于離群值沖擊的干預傳遞函數(shù)設計

    在X-13ARIMA-SEATS 季節(jié)調整中,通常將離群值分為4 種類型,第1 類為突發(fā)可加(Additive Outlier,AO)離群值點,是在某一時刻T受到異常突發(fā)事件沖擊;第2 類為水平漂移(Level Shift,LS)離群值點,是從某時刻T瞬間變化到一個新的水平,并一直保持這一水平;第3 類為暫時變化(Temporary Change,TC)離群值點,為干預經過指數(shù)衰減回到原來水平;第4 類為斜線上升(Ramp Effect,RP)離群值點,為從T1時刻開始以線性速度逐漸變化到新的水平T2[10]。各離群值點表示形式如下。

    式中:AOt為t時刻的突發(fā)可加離群值;t為突發(fā)事件干預發(fā)生的時間變量;T為突發(fā)事件干預時間點;LSt為t時刻的水平漂移離群值;TCt為t時刻的暫時變化離群值;α為經指數(shù)衰減回到原始水平的速率;RPt為t時刻的斜線上升離群值;T1為線性速率逐漸變化起點;T2為線性速率逐漸變化終點。

    在時間序列干預模型中,將干預函數(shù)分為2 種類型,一種是在T時刻短期影響的暫時的脈沖函數(shù),一種是從T時刻長期影響的持續(xù)的階梯函數(shù),干預函數(shù)各類型公式如下。

    根據(jù)離群值和干預函數(shù)表示形式,并結合兩者表現(xiàn)形式特點,使用干預函數(shù)構建離群值干預變量,此時,Pt(T)和St(T)被稱為虛擬變量。根據(jù)干預過程影響時間的持續(xù)性特征,使用ARMA 模型表示干預傳遞函數(shù)。

    假設干預為AO 離群值類型的暫時性干預,即只在T時刻(或滯后q階)影響時間序列均值函數(shù),則干預大小表示為

    式中:Xt為干預大?。粀為干預函數(shù)系數(shù);B為滯后算子;q為算子滯后階數(shù)。

    以公式 ⑼ 為干預一般表示形式,結合離群值特點,4類離群值干預變化過程的傳遞函數(shù)表示如下[1-2]。

    (1)AO 離群值導致干預發(fā)生,干預傳遞函數(shù)顯示在T時刻(或干預滯后q階)產生w的干預沖擊。AO 離群值干預傳遞函數(shù)XtAO表示為

    (2)LS 離群值導致的干預突然發(fā)生,使時間序列均值發(fā)生永久性變化,其干預傳遞函數(shù)表示在T時刻(或干預滯后q階)處突然產生w的永久性干預。LS 離群值干預傳遞函數(shù)XtLS表示為

    (3)TC 離群值導致干預開始,使時間序列均值產生短暫影響,干預傳遞函數(shù)顯示在T時刻(或干預滯后q階)處從干預w以指數(shù)級δ p逐漸回歸到原來水平。TC 離群值干預傳遞函數(shù)XtTC表示為

    式中:δ為自回歸多項式系數(shù)。

    (4)RP 離群值導致干預發(fā)生,并產生永久性新水平變化,其干預傳遞函數(shù)顯示在沖擊的滯后q階處從w1+w2以指數(shù)級δ p逐漸變化到新水平w2水平。RP 離群值干預傳遞函數(shù)XtRP表示為式中:w1與w2分別為滑動平均多項式系數(shù)。

    將以上各類干預傳遞函數(shù)作為季節(jié)調整識別的最優(yōu)ARIMA 模型的外生傳遞函數(shù)變量,即可得到如下時間序列干預模型。

    式中:?d為d階非季節(jié)差分;?SD為季節(jié)周期為S(此處周期為12)的D階差分;Yt為鐵路客運量時間序列;θ1,θ2,…,θq為殘差非季節(jié)移動平均系數(shù);?1,?2,…,?Q為殘差季節(jié)移動平均系數(shù);QS為季節(jié)移動平均滯后階數(shù);φ1,φ2,…,φp為殘差非季節(jié)自回歸系數(shù);φ1,φ2,…,φP為殘差季節(jié)自回歸系數(shù);PS為季節(jié)自回歸滯后階數(shù);εt為零均值白噪音序列。

    公式 ⒁ 中,前一部分為干預傳遞函數(shù);后一部分為ARMA 模型部分,分子為q階移動平均(MA)和QS階移動平均季節(jié)多項式乘積,分母為p階自回歸(AR)和PS階自回歸季節(jié)多項式乘積。

    2 鐵路客運量干預模型構建

    2.1 數(shù)據(jù)來源與特征分析

    研究分析2002 年1 月—2020 年8 月鐵路客運量數(shù)據(jù),鐵路客運量時間序列分布曲線如圖1 所示。根據(jù)圖1,鐵路客運量表現(xiàn)出較強的季節(jié)性特征,且在2002 年與2020 年出現(xiàn)急劇下滑現(xiàn)象,說明有突發(fā)事件對鐵路客運量造成重大激冷沖擊。

    圖1 鐵路客運量時間序列分布曲線Fig.1 Time-series distribution curve of railway passenger volume

    根據(jù)數(shù)據(jù)季節(jié)性特征顯著、異常值凸顯和月度數(shù)據(jù)特點,在考慮季節(jié)性、周期性、發(fā)展趨勢和異常值影響等條件下,采用主要分析月度尺度和季度尺度時間序列的X-13ARIMA-SEATS 模型對其進行季節(jié)調整、最優(yōu)ARIMA 模型識別和離群值檢測。X-13ARIMA-SEATS 模型是X-12ARIMA 季節(jié)調整模型和TRAMO-SEATS 季節(jié)調整模型的結合,其優(yōu)點一是增加X-12ARIMA 季節(jié)調整中Reg-ARIMA 模塊的延展時間序列、自動識別ARIMA模型和檢測離群值能力,增加移動假日(如中國春節(jié)等)和交易日效應及預測能力;二是使用SEATS程序和X-11 非參數(shù)調整方法,產生基于季節(jié)因素建模的ARIMA 預調整和預測方法[11-12]。

    研究主要應用X-13ARIMA-SEATS 模型的最優(yōu)ARIMA 模型自動識別、離群值檢測和預測功能進行干預模型的設計與分析,得到季節(jié)調整結果。鐵路客運量的X-13ARIMA-SEATS 季節(jié)調整結果如表1 所示。

    表1 鐵路客運量的X-13ARIMA-SEATS 季節(jié)調整結果Tab.1 X-13ARIMA-SEATS seasonal adjustment results of railway passenger volume

    根據(jù)表1 季節(jié)調整結果顯示,去除季節(jié)效應和離群值的最優(yōu)ARIMA 模型為ARIMA (0,1,1)×(0,1,1)12,殘差季節(jié)性(QS= 1.021)、相關性(Box-Ljung= 35.00)和正態(tài)性檢驗(Shapiro= 0.975)通過無季節(jié)性、無自相關性和正態(tài)性檢驗,模型各個統(tǒng)計量的綜合Q統(tǒng)計量為0.740,值小于1,模型有效,該ARIMA 模型可用。模型結果顯示,非季節(jié)與季節(jié)移動平均系數(shù)顯示強顯著性,體現(xiàn)外部環(huán)境對鐵路客運量影響顯著,季節(jié)效應顯示鐵路客運量具有顯著的閏年正效應和工作日負效應,檢測到影響顯著的5 個離群值。根據(jù)歷年影響鐵路客運量的重大突發(fā)事件信息,2003 年5 月離群值(AO2003.May)為SARS 疫情對鐵路客運量的顯著激冷沖擊,2014 年12 月離群值(AO2014.Dec)為客票實名制全面實施的政策性顯著正向沖擊,2020 年2 月離群值(AO2020.Feb,LS2020.Feb)為新冠肺炎疫情初期引起的顯著負向沖擊,2020 年5 月離群值(LS2020.May)為新冠肺炎疫情解禁令對鐵路客運量引起的顯著正向沖擊。造成沖擊的離群值主要分為2 類,一類為突然發(fā)生且瞬間變化到某一水平的可加AO離群值,另一類為從某時刻瞬間變化到一個新水平的LS 離群值。在檢測到各異常值顯著性沖擊特征存在的條件下,需要對離群值干預趨勢作進一步研究,以分析其沖擊變化特點。

    2.2 鐵路客運量干預模型設計

    根據(jù)表1 中檢測的離群值特點,AO2003.May離群值為2003 年5 月SARS 疫情的暫時沖擊,其干預傳遞函數(shù)可以使用公式 ⑽ 表示,考慮其干預減弱趨勢,采用具有變化趨勢的公式 ⑿ 體現(xiàn)其干預和變化趨勢;同理,AO2014.Dec 離群值為2015 年1 月實施的鐵路客票實名制政策的前置短暫干預影響,其干預傳遞函數(shù)可以用公式 ⑿ 表示;AO2020.Feb 與LS2020.Feb 離群值為2020 年2 月新冠肺炎疫情對鐵路客運量的干預,該干預存在短暫強沖擊和后期持續(xù)弱沖擊,可以采用公式 ⑾ 和公式 ⑿ 的綜合沖擊函數(shù)表示;LS2020.May 離群值干預為2020 年4 月新冠肺炎疫情解禁令實施的長期影響,可以使用公式 ⑾ 表示。將離群值的干預傳遞函數(shù)替換表1 季節(jié)調整模型中檢測出的離群值變量,作為干預模型中突發(fā)事件干預部分。

    對于閏年效應和工作日效應,將其作為不需傳遞變換變量放入模型,其閏年和工作日效應變化借鑒X-13ARIMA-SEATS 閏年效應和工作日效應虛擬變量的構造方式[13]。

    相對于表1 中X-13ARIMA-SEATS 季節(jié)調整基礎模型,該模型不僅能體現(xiàn)離群值代表的突發(fā)事件的影響程度大小,更能通過自回歸模式顯示其突發(fā)事件的影響變化過程[14-15]。干預模型算法使用R語言實現(xiàn)[16-17]。

    3 鐵路客運量干預特征分析

    3.1 鐵路客運量干預模型分析

    根據(jù)構建的鐵路客運量干預模型(公式 ⒄),對鐵路客運量進行基于表1 最優(yōu)季節(jié)調整模型(ARIMA (0,1,1)×(0,1,1)12) 的干預分析,得到季節(jié)調整的干預模型分析結果。鐵路客運量干預模型分析結果如表2 所示。根據(jù)表2,在模型ARIMA (0,1,1)×(0,1,1)12基礎上,添加的各干預變量都表現(xiàn)出強顯著性影響,相對模型未添加干預因素的(0,1,1)×(0,1,1)12,該模型AIC值和極大似然估計值分別為-387.93 和205.97,顯示出更優(yōu)的模型擬合效果。

    表2 鐵路客運量干預模型分析結果Tab.2 Analysis results of intervention model for railway passenger volume

    根據(jù)表2,閏年效應(LeapYear)和工作日效應(Workday)顯著性影響程度變化不大,與季節(jié)調整模型結果基本一致;季節(jié)(sma1)和非季節(jié)(ma1)滑動平均部分系數(shù)由表1 中的正顯著影響變?yōu)樨擄@著影響,體現(xiàn)出客運量本身的外在隨機影響因素影響變小;離群值作為干預變量后,以激冷沖擊為主的離群值對模型隨機部分的負向影響加大,模型體現(xiàn)離群值表征的突發(fā)事件對鐵路客運量的影響更顯著。

    3.2 鐵路客運量離群值干預影響分析

    根據(jù)表2,SARS 疫情對鐵路客運量表現(xiàn)出顯著性負影響,并以0.437p的指數(shù)形式逐漸減弱,SARS 疫情在滯后的第2 期顯示出較強的恢復趨勢,季節(jié)與非季節(jié)一階差分鐵路客運量離群值干預與恢復趨勢比較如表3 所示,在沖擊滯后的第5 ~ 6 期鐵路客運量基本得到恢復,8 ~ 9期后基本完全恢復,這與智研咨詢調查結果基本一致[18]。從2002 年12月SARS 疫情出現(xiàn)初期到2003 年4 月疫情基本控制,鐵路客運量未受到很大的沖擊,直到2003 年4 月,一些省市(如廣東省)取消一系列團體旅行,同時世界衛(wèi)生組織發(fā)出SARS 疫情的全球警告,我國鐵路客運量在2003 年5 月(突發(fā)期)開始驟然下滑,5 月同比下降61.18%,7 月開始回穩(wěn),8 月增速回正,2004 年1 月(滯后8 期)完全恢復。

    表3 季節(jié)與非季節(jié)一階差分鐵路客運量離群值干預與恢復趨勢比較Tab.3 Comparison of intervention and recovery trend of outliers in seasonal and non-seasonal first-order differential railway passenger volumes

    新冠肺炎疫情初期首先對鐵路客運量以w3=-1.784 的強永久性顯著負影響,同時疊加以w4=-0.996 的顯著性負影響,之后以0.550p的指數(shù)形式逐漸減弱。根據(jù)其恢復趨勢,在2020 年12 月(滯后10 期),疫情初期帶來的巨大激冷沖擊只恢復到35.73%,疫情解禁令對鐵路客運量帶來的w5= 0.289的正向永久的顯著影響持續(xù)存在,與疫情初期影響合并,在2020 年12 月(滯后10 期),疫情初期帶來的巨大激冷沖擊恢復到46.12%。

    與疫情帶來的負向顯著影響不同,鐵路客票實名制政策對鐵路客運量顯著性影響是正負波動性逐漸減弱過程,體現(xiàn)出政策具有較強的社會性,在實施過程中還會出現(xiàn)一定的不穩(wěn)定性,且該政策對客運量影響時間較短。從2015 年1 月1 日全國全面實施鐵路客票實名制政策,受政策利空影響,2014 年12 月提前出現(xiàn)鐵路客運量大幅增加的激熱效應,具體實施實名制后,客運量出現(xiàn)小幅波動,在2015 年3 月(滯后3 期)基本完全恢復,體現(xiàn)出鐵路部門前期宣傳、試點工作很深入,廣大旅客能夠很快適應政策的變化。

    根據(jù)以上分析結果,得到的各離群值干預影響的趨勢變化,各種離群值對鐵路客運量干預比較圖如圖2 所示。圖2 直觀顯示出SARS 疫情與鐵路客票實名制政策對鐵路客運量影響時間比較短,分別在6 個月和4 個月后波動基本趨向為0;新冠肺炎疫情初期影響波動較大,2020 年2 月為沖擊最大點,后期開始緩慢恢復,在2020 年4 月解禁令和復工復產等利好政策的推動下,新冠肺炎疫情初期突發(fā)性帶來的激冷沖擊基本完全釋放,當前,處于新冠肺炎疫情防控措施和解禁令帶來的永久效應影響,這種效應會在較長時間內持續(xù)存在。

    圖2 各種離群值對鐵路客運量干預比較圖Fig.2 Comparison chart of various outliers intervening railway passenger volume

    3.3 鐵路客運量趨勢預測分析

    對季節(jié)調整模型與基于離群值特征的干預模型擬合值進行相對誤差比較分析,季節(jié)調整模型與干預模型絕對相對誤差比較如表4所示。分析結果顯示,干預模型絕對相對誤差相對季節(jié)調整模型低10%,顯示干預模型擬合結果與實際值更接近,對數(shù)據(jù)的擬合更有效。因此,將離群值作為干預變量進行分析能夠提高鐵路客運量的擬合精度。季節(jié)調整曲線、干預模型擬合曲線與原始曲線比較如圖3 所示。

    根據(jù)圖3,季節(jié)調整曲線消除各類季節(jié)效應的影響,調整曲線相對原始曲線光滑很多,可以較為客觀地體現(xiàn)鐵路客運量的整體發(fā)展趨勢;干預模型擬合曲線與原始曲線基本重合,可以更好地體現(xiàn)干預模型對鐵路客運量的信息解釋程度。在當前新冠肺炎疫情管控條件下,疫情對鐵路客運量沖擊在逐漸減小,鐵路客運量逐漸回暖,總體處于上升趨勢。

    4 結束語

    研究基于SARS、新冠肺炎疫情和鐵路客票實名制政策等主要干預事件構建我國鐵路客運量干預模型,從干預理論和可視化過程2 個方面分析突發(fā)事件的沖擊過程,在一定程度上能夠厘清突發(fā)事件對我國鐵路客運量干預的發(fā)展趨勢,該方法可以拓展到我國交通系統(tǒng)突發(fā)事件干預的多個方面。在實際研究過程中,由于數(shù)據(jù)的不可獲得性,沒有考慮鐵路線路增加和區(qū)域特點對鐵路客運量的影響,在不影響整體發(fā)展趨勢的條件下一定程度上會弱化或強化突發(fā)事件干預因素的干預強度,這也是當前不能深入進行這方面干預分析的難題所在。

    猜你喜歡
    客運量離群傳遞函數(shù)
    交通運輸部:3 月城市軌道交通客運量環(huán)比增長16.6%
    2018年北京市城市公共交通運行特征分析
    交通工程(2020年5期)2020-10-21 08:45:44
    2018年北京市軌道交通運行特征分析
    交通工程(2020年2期)2020-06-03 01:10:58
    基于LabVIEW的火焰?zhèn)鬟f函數(shù)測量系統(tǒng)
    測控技術(2018年7期)2018-12-09 08:58:46
    紅外成像系統(tǒng)的調制傳遞函數(shù)測試
    離群數(shù)據(jù)挖掘在發(fā)現(xiàn)房產銷售潛在客戶中的應用
    離群的小雞
    應用相似度測量的圖離群點檢測方法
    基于傳遞函數(shù)自我優(yōu)化的BP網絡算法改進
    一種基于核空間局部離群因子的離群點挖掘方法
    久久久久久国产a免费观看| 99久国产av精品| 99在线视频只有这里精品首页| 国产免费一级a男人的天堂| 在现免费观看毛片| 日本 av在线| 欧美日本视频| 麻豆成人午夜福利视频| 亚洲自偷自拍三级| 免费看光身美女| 久久精品国产亚洲网站| 亚洲自偷自拍三级| 婷婷精品国产亚洲av| 美女高潮的动态| 亚洲av一区综合| 身体一侧抽搐| 高清在线国产一区| 毛片一级片免费看久久久久 | 日韩欧美国产一区二区入口| 国产色爽女视频免费观看| 久久午夜福利片| 一进一出好大好爽视频| 亚洲电影在线观看av| 九九爱精品视频在线观看| 偷拍熟女少妇极品色| 人人妻人人澡欧美一区二区| 99久久成人亚洲精品观看| 波多野结衣高清作品| 久久欧美精品欧美久久欧美| 中文字幕高清在线视频| 亚洲一区高清亚洲精品| 成人无遮挡网站| 免费不卡的大黄色大毛片视频在线观看 | 亚洲国产精品成人综合色| avwww免费| 日本成人三级电影网站| 真人做人爱边吃奶动态| xxxwww97欧美| 一个人观看的视频www高清免费观看| 男人舔女人下体高潮全视频| 在线a可以看的网站| 日本一本二区三区精品| 久99久视频精品免费| 欧美性猛交╳xxx乱大交人| 亚洲精品一区av在线观看| 亚洲色图av天堂| 一进一出抽搐动态| 成年女人永久免费观看视频| 熟女电影av网| 精品人妻一区二区三区麻豆 | 日韩一区二区视频免费看| 国内毛片毛片毛片毛片毛片| 亚洲专区国产一区二区| 国产探花在线观看一区二区| 亚洲欧美日韩卡通动漫| 午夜福利高清视频| 亚洲成av人片在线播放无| 欧美色视频一区免费| 床上黄色一级片| 亚洲自拍偷在线| 国产精品乱码一区二三区的特点| 成人性生交大片免费视频hd| 国内毛片毛片毛片毛片毛片| 国产蜜桃级精品一区二区三区| 在线观看免费视频日本深夜| 91狼人影院| 观看免费一级毛片| 无人区码免费观看不卡| 69av精品久久久久久| 国产精品国产三级国产av玫瑰| 2021天堂中文幕一二区在线观| 中文亚洲av片在线观看爽| 成人高潮视频无遮挡免费网站| 亚洲人成网站高清观看| 狂野欧美白嫩少妇大欣赏| 成人午夜高清在线视频| 18禁黄网站禁片午夜丰满| 亚洲黑人精品在线| 午夜精品久久久久久毛片777| 很黄的视频免费| 美女cb高潮喷水在线观看| 久99久视频精品免费| 国产精品久久久久久久久免| 国产综合懂色| 看片在线看免费视频| 国产伦一二天堂av在线观看| av天堂中文字幕网| 成年女人看的毛片在线观看| 日日啪夜夜撸| 久久精品久久久久久噜噜老黄 | 男人的好看免费观看在线视频| 久久久久性生活片| 精品人妻一区二区三区麻豆 | 人妻制服诱惑在线中文字幕| 久久这里只有精品中国| 毛片女人毛片| 国产在线男女| 无遮挡黄片免费观看| 亚洲最大成人av| 夜夜夜夜夜久久久久| 亚洲国产日韩欧美精品在线观看| 欧美日韩精品成人综合77777| 日本a在线网址| 97超级碰碰碰精品色视频在线观看| 午夜免费男女啪啪视频观看 | 亚洲最大成人av| 99热这里只有是精品在线观看| av在线亚洲专区| 亚洲美女黄片视频| 啪啪无遮挡十八禁网站| 十八禁国产超污无遮挡网站| 男人舔奶头视频| 亚洲最大成人av| 麻豆av噜噜一区二区三区| 国产伦在线观看视频一区| 欧美日本亚洲视频在线播放| 丝袜美腿在线中文| 国产又黄又爽又无遮挡在线| 全区人妻精品视频| 丰满乱子伦码专区| 亚洲av不卡在线观看| 看免费成人av毛片| 亚洲一级一片aⅴ在线观看| 男插女下体视频免费在线播放| 精品一区二区三区av网在线观看| 免费观看在线日韩| 乱人视频在线观看| 精品久久久久久久久久久久久| 亚洲第一区二区三区不卡| 欧美黑人欧美精品刺激| 少妇人妻精品综合一区二区 | 国产淫片久久久久久久久| 啦啦啦韩国在线观看视频| 国产精品综合久久久久久久免费| 国产精品精品国产色婷婷| 国产亚洲欧美98| 国产免费男女视频| 色噜噜av男人的天堂激情| 夜夜爽天天搞| 国产色婷婷99| 午夜福利18| а√天堂www在线а√下载| 精品日产1卡2卡| 网址你懂的国产日韩在线| 久久久久久久久大av| 真人一进一出gif抽搐免费| 国产在线男女| 国产高清激情床上av| 日韩av在线大香蕉| 免费一级毛片在线播放高清视频| 国产三级在线视频| 久久天躁狠狠躁夜夜2o2o| 国产又黄又爽又无遮挡在线| 99国产极品粉嫩在线观看| 两个人视频免费观看高清| 别揉我奶头 嗯啊视频| 国产熟女欧美一区二区| 国产精品精品国产色婷婷| 极品教师在线视频| 日本-黄色视频高清免费观看| 国产av不卡久久| 97超级碰碰碰精品色视频在线观看| 一个人看视频在线观看www免费| 国产精品av视频在线免费观看| 嫁个100分男人电影在线观看| 可以在线观看毛片的网站| 最好的美女福利视频网| 男插女下体视频免费在线播放| 国产精品野战在线观看| 亚洲美女视频黄频| 日韩 亚洲 欧美在线| 韩国av一区二区三区四区| 国产真实伦视频高清在线观看 | 两个人的视频大全免费| 久久久久久大精品| 精品99又大又爽又粗少妇毛片 | 51国产日韩欧美| 在线国产一区二区在线| 国产精品一及| 亚洲成人久久爱视频| 免费一级毛片在线播放高清视频| 少妇人妻一区二区三区视频| 亚洲四区av| 欧美高清成人免费视频www| www日本黄色视频网| 一本久久中文字幕| 又黄又爽又免费观看的视频| 久久精品影院6| 成人二区视频| 欧美区成人在线视频| 99热网站在线观看| 欧美性猛交黑人性爽| 91午夜精品亚洲一区二区三区 | 国产视频内射| 日日摸夜夜添夜夜添小说| 免费看光身美女| 亚洲自拍偷在线| 嫩草影院精品99| 1024手机看黄色片| 乱码一卡2卡4卡精品| 国产精品1区2区在线观看.| 一级a爱片免费观看的视频| 国产精品一区二区免费欧美| 成人高潮视频无遮挡免费网站| 一a级毛片在线观看| 女人十人毛片免费观看3o分钟| 久久久久性生活片| 国产伦精品一区二区三区四那| 12—13女人毛片做爰片一| 精品乱码久久久久久99久播| videossex国产| 日本三级黄在线观看| 熟妇人妻久久中文字幕3abv| 国产男人的电影天堂91| 久久精品国产亚洲av天美| 97碰自拍视频| 麻豆国产av国片精品| 内射极品少妇av片p| 91麻豆精品激情在线观看国产| 97超级碰碰碰精品色视频在线观看| 欧美区成人在线视频| 极品教师在线免费播放| 99久久精品一区二区三区| 国产一区二区在线观看日韩| 中文字幕精品亚洲无线码一区| 在线观看66精品国产| 丰满乱子伦码专区| 午夜精品一区二区三区免费看| 少妇人妻一区二区三区视频| 一级黄片播放器| 麻豆成人av在线观看| 极品教师在线免费播放| 国产伦精品一区二区三区四那| 久久这里只有精品中国| 国产精品久久久久久久久免| 此物有八面人人有两片| 久久午夜亚洲精品久久| 国产精品亚洲一级av第二区| 欧美性猛交╳xxx乱大交人| 又黄又爽又刺激的免费视频.| 亚洲欧美日韩无卡精品| 美女 人体艺术 gogo| av天堂中文字幕网| 久久6这里有精品| 窝窝影院91人妻| 黄色女人牲交| 亚洲人成网站在线播| 免费在线观看影片大全网站| 成人无遮挡网站| 中文资源天堂在线| 国产欧美日韩精品亚洲av| 精品久久久久久成人av| 日日啪夜夜撸| 99久久无色码亚洲精品果冻| 国产主播在线观看一区二区| 亚洲经典国产精华液单| 国内揄拍国产精品人妻在线| 欧美3d第一页| 亚洲真实伦在线观看| 午夜免费激情av| 熟妇人妻久久中文字幕3abv| 一级黄色大片毛片| 两个人视频免费观看高清| 中文字幕av在线有码专区| 他把我摸到了高潮在线观看| 51国产日韩欧美| 日本黄色视频三级网站网址| a级毛片免费高清观看在线播放| 高清在线国产一区| 色噜噜av男人的天堂激情| 一边摸一边抽搐一进一小说| 亚洲最大成人中文| 亚洲一级一片aⅴ在线观看| 人妻久久中文字幕网| 国语自产精品视频在线第100页| 九九热线精品视视频播放| 可以在线观看的亚洲视频| 亚洲成人免费电影在线观看| 日日干狠狠操夜夜爽| 免费黄网站久久成人精品| 人人妻,人人澡人人爽秒播| 在线观看美女被高潮喷水网站| 黄色欧美视频在线观看| 日日摸夜夜添夜夜添av毛片 | 乱系列少妇在线播放| 亚洲性久久影院| АⅤ资源中文在线天堂| 中文在线观看免费www的网站| 我的老师免费观看完整版| 中国美女看黄片| 99热这里只有是精品50| 91麻豆精品激情在线观看国产| 麻豆成人午夜福利视频| 亚洲最大成人中文| 97人妻精品一区二区三区麻豆| 亚洲精品影视一区二区三区av| 成人永久免费在线观看视频| 久久精品久久久久久噜噜老黄 | 日本黄色视频三级网站网址| 免费电影在线观看免费观看| 老司机福利观看| 国产av在哪里看| av视频在线观看入口| 看黄色毛片网站| 日韩中字成人| 久久草成人影院| 免费无遮挡裸体视频| 美女xxoo啪啪120秒动态图| 两人在一起打扑克的视频| 又紧又爽又黄一区二区| 亚洲精品国产成人久久av| 色av中文字幕| 亚洲美女黄片视频| 可以在线观看的亚洲视频| 国内精品一区二区在线观看| 我要搜黄色片| 久久久久精品国产欧美久久久| 久久国产乱子免费精品| 性色avwww在线观看| 看免费成人av毛片| 精品久久国产蜜桃| a在线观看视频网站| 可以在线观看毛片的网站| 天堂影院成人在线观看| 一卡2卡三卡四卡精品乱码亚洲| 人人妻人人看人人澡| 99久久精品热视频| 国产精品久久视频播放| 九九爱精品视频在线观看| 国语自产精品视频在线第100页| 久久中文看片网| 亚洲国产精品成人综合色| 免费在线观看日本一区| 亚洲人成网站在线播放欧美日韩| 最新中文字幕久久久久| 精品人妻熟女av久视频| 波多野结衣巨乳人妻| 成人毛片a级毛片在线播放| 国内揄拍国产精品人妻在线| bbb黄色大片| 91狼人影院| 天堂动漫精品| 久久人妻av系列| 欧美三级亚洲精品| 精品一区二区免费观看| 国产美女午夜福利| 精品久久久久久久人妻蜜臀av| 国产av在哪里看| 乱人视频在线观看| 校园春色视频在线观看| 麻豆一二三区av精品| 91狼人影院| 亚洲av免费在线观看| 久久国内精品自在自线图片| 日韩一区二区视频免费看| 精品不卡国产一区二区三区| 两个人的视频大全免费| 成人特级av手机在线观看| 99在线视频只有这里精品首页| 在线播放无遮挡| 国产三级在线视频| 亚洲国产精品成人综合色| 中出人妻视频一区二区| 麻豆一二三区av精品| eeuss影院久久| 制服丝袜大香蕉在线| 内地一区二区视频在线| 99久久九九国产精品国产免费| 久久久久久伊人网av| 午夜亚洲福利在线播放| 一本精品99久久精品77| av福利片在线观看| 久久久精品欧美日韩精品| 啦啦啦啦在线视频资源| 日本五十路高清| 非洲黑人性xxxx精品又粗又长| 国产精品久久久久久亚洲av鲁大| 国产精品一区www在线观看 | 亚洲精品粉嫩美女一区| 天堂av国产一区二区熟女人妻| 国产精品日韩av在线免费观看| 国产人妻一区二区三区在| 国内少妇人妻偷人精品xxx网站| 日韩 亚洲 欧美在线| 午夜福利高清视频| 国产亚洲av嫩草精品影院| 999久久久精品免费观看国产| АⅤ资源中文在线天堂| 国产成人a区在线观看| 88av欧美| 亚洲avbb在线观看| 久久精品夜夜夜夜夜久久蜜豆| 99精品久久久久人妻精品| 日韩人妻高清精品专区| x7x7x7水蜜桃| 亚洲成人久久性| 国产麻豆成人av免费视频| 亚洲国产精品合色在线| 日本 欧美在线| 国产精品亚洲一级av第二区| 欧美zozozo另类| 婷婷亚洲欧美| 国产一区二区亚洲精品在线观看| 小蜜桃在线观看免费完整版高清| 亚洲av美国av| 男女那种视频在线观看| 变态另类丝袜制服| 成年女人毛片免费观看观看9| 97人妻精品一区二区三区麻豆| 成人亚洲精品av一区二区| 91久久精品电影网| 免费在线观看成人毛片| 欧美区成人在线视频| 亚洲四区av| 午夜福利高清视频| 毛片女人毛片| 日韩欧美国产一区二区入口| 国内精品久久久久久久电影| 国产乱人视频| 亚洲精品456在线播放app | 熟女人妻精品中文字幕| av天堂在线播放| 男插女下体视频免费在线播放| 精品一区二区三区人妻视频| 亚洲国产精品成人综合色| 小说图片视频综合网站| 给我免费播放毛片高清在线观看| 伦理电影大哥的女人| 国产成人影院久久av| 国产亚洲精品av在线| 国产男人的电影天堂91| 91在线观看av| 少妇的逼好多水| 国产在线男女| 偷拍熟女少妇极品色| 成人精品一区二区免费| 一个人免费在线观看电影| 人妻丰满熟妇av一区二区三区| 亚洲av日韩精品久久久久久密| 国产亚洲精品av在线| 看免费成人av毛片| 有码 亚洲区| 久久国产精品人妻蜜桃| 色综合站精品国产| 免费看光身美女| 偷拍熟女少妇极品色| 少妇高潮的动态图| 日韩欧美精品免费久久| 欧美不卡视频在线免费观看| 国内精品宾馆在线| 不卡视频在线观看欧美| 一本久久中文字幕| 日韩欧美国产一区二区入口| 成年人黄色毛片网站| 日韩精品青青久久久久久| 美女高潮的动态| 此物有八面人人有两片| 亚洲人成网站在线播| 免费不卡的大黄色大毛片视频在线观看 | 少妇猛男粗大的猛烈进出视频 | 性色avwww在线观看| 亚洲中文字幕一区二区三区有码在线看| 国产精品一区www在线观看 | 韩国av在线不卡| 精品久久久久久久久亚洲 | av国产免费在线观看| 国产在视频线在精品| 亚洲成人精品中文字幕电影| 午夜久久久久精精品| 亚洲国产欧美人成| 成人精品一区二区免费| 欧美在线一区亚洲| 国产一区二区三区视频了| 精品久久久久久久末码| 人妻夜夜爽99麻豆av| 免费在线观看成人毛片| 亚洲精品一区av在线观看| 亚洲国产精品久久男人天堂| 在线观看av片永久免费下载| 午夜福利成人在线免费观看| 久久久久久久久久黄片| 久久国产精品人妻蜜桃| 听说在线观看完整版免费高清| 亚洲国产精品合色在线| 成人午夜高清在线视频| 三级国产精品欧美在线观看| x7x7x7水蜜桃| 久久久午夜欧美精品| 亚洲人成网站在线播| 国产一级毛片七仙女欲春2| 午夜激情福利司机影院| 亚洲美女视频黄频| 精华霜和精华液先用哪个| 成熟少妇高潮喷水视频| 88av欧美| 国产精品99久久久久久久久| 噜噜噜噜噜久久久久久91| 国语自产精品视频在线第100页| 99在线视频只有这里精品首页| 亚洲第一区二区三区不卡| 97超级碰碰碰精品色视频在线观看| 91午夜精品亚洲一区二区三区 | 无人区码免费观看不卡| 国产精品亚洲美女久久久| 久久婷婷人人爽人人干人人爱| 又爽又黄无遮挡网站| 亚洲人成网站在线播放欧美日韩| 日韩av在线大香蕉| 观看免费一级毛片| 我的女老师完整版在线观看| 国产精品伦人一区二区| 亚洲精品在线观看二区| 夜夜爽天天搞| 国产白丝娇喘喷水9色精品| 久久久久久久午夜电影| 中文字幕久久专区| 亚洲精品成人久久久久久| 欧美丝袜亚洲另类 | 国内毛片毛片毛片毛片毛片| 国产在线精品亚洲第一网站| 久久久久久久久久黄片| 久久草成人影院| 婷婷色综合大香蕉| 美女cb高潮喷水在线观看| 免费观看在线日韩| 99九九线精品视频在线观看视频| 99热这里只有是精品50| 亚洲av免费在线观看| 99热这里只有是精品50| 国产av不卡久久| 天堂av国产一区二区熟女人妻| 女生性感内裤真人,穿戴方法视频| 亚洲精品456在线播放app | 亚洲美女搞黄在线观看 | 色哟哟·www| 天堂网av新在线| 色哟哟·www| 真人做人爱边吃奶动态| 91久久精品国产一区二区三区| 一级黄片播放器| .国产精品久久| 男人和女人高潮做爰伦理| 久久国产乱子免费精品| 日韩大尺度精品在线看网址| 99热只有精品国产| 日本一本二区三区精品| а√天堂www在线а√下载| 国产一区二区三区视频了| 亚洲成人久久爱视频| 久久精品国产亚洲av香蕉五月| 免费一级毛片在线播放高清视频| 最近视频中文字幕2019在线8| 熟女电影av网| 国产亚洲欧美98| 少妇被粗大猛烈的视频| 麻豆一二三区av精品| 午夜福利在线观看吧| 日本成人三级电影网站| 欧美丝袜亚洲另类 | 成人午夜高清在线视频| 国产伦在线观看视频一区| 一夜夜www| 亚洲va日本ⅴa欧美va伊人久久| 麻豆成人午夜福利视频| 久久这里只有精品中国| 亚洲天堂国产精品一区在线| 国产精品爽爽va在线观看网站| 九九在线视频观看精品| 一个人观看的视频www高清免费观看| .国产精品久久| 夜夜夜夜夜久久久久| 国产高清视频在线观看网站| 亚洲avbb在线观看| 国产精品一及| 亚洲三级黄色毛片| 色吧在线观看| 国产在线精品亚洲第一网站| 少妇丰满av| 一卡2卡三卡四卡精品乱码亚洲| 国产老妇女一区| 乱系列少妇在线播放| 国产三级在线视频| 能在线免费观看的黄片| 久久久久久九九精品二区国产| 男女视频在线观看网站免费| 国产精品1区2区在线观看.| 欧美性猛交╳xxx乱大交人| 熟妇人妻久久中文字幕3abv| 有码 亚洲区| 国产午夜精品论理片| 97超视频在线观看视频| 国内久久婷婷六月综合欲色啪| 久久精品国产清高在天天线| 极品教师在线视频| 热99在线观看视频| 国产精品永久免费网站| 美女大奶头视频| 欧美极品一区二区三区四区| 91午夜精品亚洲一区二区三区 | 国产激情偷乱视频一区二区| 超碰av人人做人人爽久久| 免费无遮挡裸体视频| 嫩草影院新地址| 在线观看66精品国产| 亚洲av第一区精品v没综合| a级一级毛片免费在线观看| 国产一区二区三区在线臀色熟女| 欧美一级a爱片免费观看看| 久久久久久大精品| 午夜精品一区二区三区免费看| 俄罗斯特黄特色一大片| 国产乱人伦免费视频| 欧美激情在线99| 久久香蕉精品热| 啪啪无遮挡十八禁网站|