• <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網絡算法改進
    一種基于核空間局部離群因子的離群點挖掘方法
    你懂的网址亚洲精品在线观看 | 亚洲三级黄色毛片| 亚洲最大成人中文| av在线观看视频网站免费| 免费高清视频大片| 国模一区二区三区四区视频| 久99久视频精品免费| 国产在视频线在精品| 少妇高潮的动态图| 草草在线视频免费看| 69人妻影院| 麻豆精品久久久久久蜜桃| 免费一级毛片在线播放高清视频| 成人av一区二区三区在线看| 久久人人爽人人爽人人片va| 国产亚洲精品久久久com| 色尼玛亚洲综合影院| 又黄又爽又免费观看的视频| 男女啪啪激烈高潮av片| 欧美性猛交黑人性爽| 99在线视频只有这里精品首页| 老师上课跳d突然被开到最大视频| 99久久九九国产精品国产免费| 午夜福利18| 噜噜噜噜噜久久久久久91| 国产精品美女特级片免费视频播放器| 精品一区二区三区av网在线观看| 日韩欧美 国产精品| 国产成年人精品一区二区| 在线观看午夜福利视频| 级片在线观看| 亚洲国产高清在线一区二区三| 久久99热这里只有精品18| 99国产精品一区二区蜜桃av| 国产乱人偷精品视频| 精品少妇黑人巨大在线播放 | 欧美三级亚洲精品| 亚洲av成人av| 热99re8久久精品国产| 看片在线看免费视频| 免费在线观看影片大全网站| 成人特级av手机在线观看| 简卡轻食公司| 97超碰精品成人国产| 国产三级中文精品| 免费看美女性在线毛片视频| 蜜桃亚洲精品一区二区三区| 久久精品国产亚洲av涩爱 | 欧美一区二区亚洲| 精品日产1卡2卡| 免费看光身美女| 久久国产乱子免费精品| 六月丁香七月| 精品一区二区三区av网在线观看| 97碰自拍视频| 舔av片在线| 久久久久性生活片| 成人高潮视频无遮挡免费网站| 日韩欧美精品v在线| 国产亚洲91精品色在线| 51国产日韩欧美| 麻豆一二三区av精品| 精品一区二区三区视频在线观看免费| 亚洲精品成人久久久久久| 日韩制服骚丝袜av| 99久久精品热视频| 亚洲精品日韩av片在线观看| 草草在线视频免费看| 国产男人的电影天堂91| 亚洲精品国产成人久久av| 亚洲av熟女| 男人和女人高潮做爰伦理| a级毛色黄片| 午夜久久久久精精品| 99久久精品一区二区三区| 伦理电影大哥的女人| 亚洲中文字幕日韩| 欧美日本亚洲视频在线播放| 久久久成人免费电影| 女生性感内裤真人,穿戴方法视频| 免费观看人在逋| 久久久国产成人免费| 欧美潮喷喷水| 91av网一区二区| 亚洲av熟女| 日本三级黄在线观看| 亚洲经典国产精华液单| 有码 亚洲区| 午夜激情欧美在线| 亚洲精华国产精华液的使用体验 | 中国国产av一级| 最新在线观看一区二区三区| 熟妇人妻久久中文字幕3abv| 搡女人真爽免费视频火全软件 | 少妇熟女aⅴ在线视频| 日韩中字成人| 99riav亚洲国产免费| av在线老鸭窝| 22中文网久久字幕| av在线播放精品| 亚洲精品乱码久久久v下载方式| 国产亚洲精品久久久久久毛片| 丰满人妻一区二区三区视频av| 久久九九热精品免费| 日韩大尺度精品在线看网址| 国产av在哪里看| 日韩成人av中文字幕在线观看 | 十八禁网站免费在线| 麻豆成人午夜福利视频| 亚洲性久久影院| 黄片wwwwww| 一区二区三区免费毛片| 精品久久久久久久人妻蜜臀av| 国产精品久久久久久亚洲av鲁大| 欧美日韩乱码在线| av在线老鸭窝| 亚洲自拍偷在线| 日日摸夜夜添夜夜爱| 久久久久精品国产欧美久久久| 日本黄色视频三级网站网址| 国国产精品蜜臀av免费| 成熟少妇高潮喷水视频| 日本黄大片高清| 国产一区二区亚洲精品在线观看| 老司机福利观看| 国产精品亚洲一级av第二区| 综合色av麻豆| 国内久久婷婷六月综合欲色啪| 夜夜看夜夜爽夜夜摸| 国产黄a三级三级三级人| 麻豆乱淫一区二区| 在线a可以看的网站| 日本精品一区二区三区蜜桃| 成人二区视频| 麻豆久久精品国产亚洲av| 九九热线精品视视频播放| 啦啦啦啦在线视频资源| 国产成人91sexporn| 中国国产av一级| 国产91av在线免费观看| 你懂的网址亚洲精品在线观看 | 高清日韩中文字幕在线| 最近2019中文字幕mv第一页| 黄色日韩在线| 亚洲国产欧洲综合997久久,| 一级av片app| 三级男女做爰猛烈吃奶摸视频| 亚洲自偷自拍三级| 成熟少妇高潮喷水视频| 国产乱人视频| 午夜福利在线观看吧| 天堂网av新在线| 国产亚洲精品久久久久久毛片| 日韩欧美一区二区三区在线观看| 91午夜精品亚洲一区二区三区| 午夜影院日韩av| 一个人免费在线观看电影| 免费搜索国产男女视频| 亚洲av不卡在线观看| 亚洲欧美清纯卡通| 成人亚洲精品av一区二区| 久久久久免费精品人妻一区二区| 97人妻精品一区二区三区麻豆| 欧美一区二区国产精品久久精品| 欧美成人a在线观看| 能在线免费观看的黄片| 精品日产1卡2卡| av在线亚洲专区| 精品一区二区三区视频在线观看免费| 99热网站在线观看| 欧美激情在线99| 在线播放国产精品三级| 国产视频一区二区在线看| 婷婷六月久久综合丁香| 性插视频无遮挡在线免费观看| 午夜福利成人在线免费观看| 亚洲性久久影院| 亚洲国产精品sss在线观看| 黄色视频,在线免费观看| 国产精品一二三区在线看| 美女被艹到高潮喷水动态| 中文字幕av成人在线电影| 狠狠狠狠99中文字幕| 99九九线精品视频在线观看视频| 尤物成人国产欧美一区二区三区| 波多野结衣高清无吗| 99久久无色码亚洲精品果冻| 亚洲国产日韩欧美精品在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 99久久精品热视频| 欧美3d第一页| ponron亚洲| 女人十人毛片免费观看3o分钟| 国产精品一区二区免费欧美| 色噜噜av男人的天堂激情| 久久国内精品自在自线图片| 看十八女毛片水多多多| 我要看日韩黄色一级片| 免费观看精品视频网站| 亚洲成人精品中文字幕电影| 99久国产av精品国产电影| 日韩在线高清观看一区二区三区| 亚洲综合色惰| 亚洲精品456在线播放app| 亚洲第一区二区三区不卡| 男女边吃奶边做爰视频| 中文字幕精品亚洲无线码一区| av黄色大香蕉| 国产精华一区二区三区| 高清毛片免费看| 亚洲性久久影院| 日本五十路高清| 中文字幕熟女人妻在线| aaaaa片日本免费| 成人美女网站在线观看视频| 成年版毛片免费区| 日本与韩国留学比较| 亚洲无线在线观看| 国产亚洲精品久久久com| 又爽又黄a免费视频| 精品少妇黑人巨大在线播放 | 国产69精品久久久久777片| 亚洲欧美成人精品一区二区| 精品熟女少妇av免费看| 菩萨蛮人人尽说江南好唐韦庄 | 国产成人freesex在线 | 91精品国产九色| 插逼视频在线观看| 日日摸夜夜添夜夜爱| 亚洲av中文字字幕乱码综合| 99久久精品热视频| 人妻久久中文字幕网| 国产女主播在线喷水免费视频网站 | 国产单亲对白刺激| 97碰自拍视频| 亚洲国产色片| 日本-黄色视频高清免费观看| 一级毛片久久久久久久久女| 天堂av国产一区二区熟女人妻| 国产蜜桃级精品一区二区三区| 午夜激情福利司机影院| 两个人的视频大全免费| 色播亚洲综合网| 人人妻人人看人人澡| 一进一出好大好爽视频| 国产精品久久电影中文字幕| 美女免费视频网站| 夜夜夜夜夜久久久久| 三级国产精品欧美在线观看| 日韩欧美精品免费久久| 成年版毛片免费区| 69人妻影院| 天堂av国产一区二区熟女人妻| 菩萨蛮人人尽说江南好唐韦庄 | 日韩成人av中文字幕在线观看 | 岛国在线免费视频观看| av在线亚洲专区| 天堂√8在线中文| 蜜桃久久精品国产亚洲av| 久久欧美精品欧美久久欧美| 成人av在线播放网站| 免费av不卡在线播放| 美女xxoo啪啪120秒动态图| 国产在线精品亚洲第一网站| 久久热精品热| 一级毛片久久久久久久久女| 精品一区二区三区人妻视频| 美女黄网站色视频| 久久精品夜色国产| 国产精品1区2区在线观看.| 国产精品三级大全| 91精品国产九色| 色哟哟哟哟哟哟| 最近中文字幕高清免费大全6| 国产69精品久久久久777片| 变态另类成人亚洲欧美熟女| 日本 av在线| 给我免费播放毛片高清在线观看| 美女被艹到高潮喷水动态| 日本色播在线视频| 最后的刺客免费高清国语| 国内久久婷婷六月综合欲色啪| 亚洲成人av在线免费| 丰满人妻一区二区三区视频av| 国产中年淑女户外野战色| 精品少妇黑人巨大在线播放 | 级片在线观看| 日韩av不卡免费在线播放| 亚洲国产精品成人综合色| 老司机福利观看| 亚洲精品一卡2卡三卡4卡5卡| 日本精品一区二区三区蜜桃| 麻豆精品久久久久久蜜桃| 国语自产精品视频在线第100页| 天堂动漫精品| 欧美另类亚洲清纯唯美| 亚洲欧美精品自产自拍| 成人漫画全彩无遮挡| 国产 一区 欧美 日韩| 日韩精品青青久久久久久| 国产精品久久久久久av不卡| 中文字幕av在线有码专区| 欧美bdsm另类| 精品一区二区三区视频在线| 亚洲成av人片在线播放无| 欧美日韩精品成人综合77777| 丝袜美腿在线中文| 日日摸夜夜添夜夜添小说| 少妇丰满av| 男人和女人高潮做爰伦理| 欧美绝顶高潮抽搐喷水| 天堂√8在线中文| 成人亚洲精品av一区二区| 精品一区二区三区视频在线观看免费| 九九热线精品视视频播放| 秋霞在线观看毛片| 少妇熟女aⅴ在线视频| 一区福利在线观看| 久久国内精品自在自线图片| 欧美日韩精品成人综合77777| a级毛片免费高清观看在线播放| 女人十人毛片免费观看3o分钟| 一区二区三区四区激情视频 | 久久草成人影院| 欧美最新免费一区二区三区| 精品午夜福利在线看| 国产亚洲欧美98| 小说图片视频综合网站| 亚洲国产欧美人成| 色av中文字幕| 麻豆久久精品国产亚洲av| 欧美日韩乱码在线| 国产午夜精品论理片| 九九热线精品视视频播放| 可以在线观看毛片的网站| 亚州av有码| 精品人妻一区二区三区麻豆 | 国产黄色视频一区二区在线观看 | 亚洲成人av在线免费| 色播亚洲综合网| 人妻制服诱惑在线中文字幕| 最近中文字幕高清免费大全6| 亚洲性久久影院| 我要看日韩黄色一级片| 国产伦一二天堂av在线观看| 免费观看人在逋| 欧美日本视频| 欧美一区二区亚洲| 天天一区二区日本电影三级| 91在线精品国自产拍蜜月| 熟妇人妻久久中文字幕3abv| av在线播放精品| 国产成人freesex在线 | 欧美日韩一区二区视频在线观看视频在线 | 少妇的逼好多水| 久久久久久久亚洲中文字幕| 国产男人的电影天堂91| 成人一区二区视频在线观看| 色综合站精品国产| 国产亚洲精品av在线| 噜噜噜噜噜久久久久久91| 国产在线男女| 欧美性感艳星| 人妻制服诱惑在线中文字幕| 国产精品久久久久久av不卡| 国产精品亚洲一级av第二区| 此物有八面人人有两片| 99久国产av精品国产电影| 免费观看在线日韩| 亚洲国产精品成人久久小说 | 国产成年人精品一区二区| 久久久a久久爽久久v久久| 色吧在线观看| 国产成人a区在线观看| 国产精品人妻久久久影院| 国产成人a区在线观看| 亚洲18禁久久av| 日韩大尺度精品在线看网址| 国国产精品蜜臀av免费| 中文资源天堂在线| 一区二区三区四区激情视频 | 99久久精品热视频| 男女那种视频在线观看| 免费高清视频大片| 久久久欧美国产精品| 激情 狠狠 欧美| 91狼人影院| 在线观看美女被高潮喷水网站| 内射极品少妇av片p| 亚洲七黄色美女视频| 特级一级黄色大片| 69av精品久久久久久| 看十八女毛片水多多多| 91在线精品国自产拍蜜月| 看片在线看免费视频| 国产一区二区三区av在线 | 欧美日本亚洲视频在线播放| 精品久久国产蜜桃| 精品久久久久久成人av| av天堂在线播放| 男人的好看免费观看在线视频| 欧美激情国产日韩精品一区| 国产男人的电影天堂91| 亚洲不卡免费看| 国产一区二区在线观看日韩| 日本撒尿小便嘘嘘汇集6| 嫩草影视91久久| 1024手机看黄色片| 婷婷精品国产亚洲av| 亚洲国产欧洲综合997久久,| 成人亚洲精品av一区二区| 亚洲精品日韩在线中文字幕 | 国产精华一区二区三区| 国产精品99久久久久久久久| 天堂√8在线中文| 干丝袜人妻中文字幕| 日本黄色视频三级网站网址| 99久久无色码亚洲精品果冻| 国产伦在线观看视频一区| 国产亚洲91精品色在线| 亚洲美女搞黄在线观看 | 久久精品国产亚洲av天美| 97碰自拍视频| 国产亚洲精品久久久com| 免费无遮挡裸体视频| 嫩草影院精品99| 日韩欧美国产在线观看| 18禁黄网站禁片免费观看直播| 国产真实乱freesex| 波野结衣二区三区在线| 最近视频中文字幕2019在线8| av在线老鸭窝| 一a级毛片在线观看| 热99re8久久精品国产| 精品熟女少妇av免费看| 久久韩国三级中文字幕| 亚洲电影在线观看av| 亚洲国产色片| 欧美日韩在线观看h| 亚洲av不卡在线观看| 国产精品久久久久久精品电影| 观看免费一级毛片| 男女边吃奶边做爰视频| 欧美绝顶高潮抽搐喷水| 在线观看午夜福利视频| 22中文网久久字幕| 中出人妻视频一区二区| 色吧在线观看| 久久久久久久久中文| 最近中文字幕高清免费大全6| .国产精品久久| 日韩亚洲欧美综合| 亚洲美女黄片视频| 国产乱人偷精品视频| 亚洲av成人av| 亚洲,欧美,日韩| 国产成人a∨麻豆精品| 午夜福利在线观看免费完整高清在 | 99热6这里只有精品| 精品午夜福利视频在线观看一区| 在线观看一区二区三区| 99热这里只有是精品50| or卡值多少钱| 欧美不卡视频在线免费观看| 免费观看的影片在线观看| 在线观看av片永久免费下载| 91精品国产九色| 国产私拍福利视频在线观看| 亚洲人成网站在线观看播放| 亚洲aⅴ乱码一区二区在线播放| 亚洲中文字幕日韩| 99热全是精品| 男女做爰动态图高潮gif福利片| 十八禁国产超污无遮挡网站| 亚洲精品国产成人久久av| 久久久久国产网址| 国产亚洲精品久久久com| 神马国产精品三级电影在线观看| 中文字幕免费在线视频6| 黑人高潮一二区| 91在线观看av| 天天躁日日操中文字幕| 久久久久久九九精品二区国产| 免费看av在线观看网站| 国产成人91sexporn| 国产av麻豆久久久久久久| 一级毛片久久久久久久久女| 插逼视频在线观看| 免费在线观看成人毛片| 久久精品国产自在天天线| 国产 一区 欧美 日韩| 91久久精品国产一区二区三区| 欧美激情久久久久久爽电影| 日本欧美国产在线视频| 亚洲自拍偷在线| 中文字幕人妻熟人妻熟丝袜美| 色哟哟·www| 天天一区二区日本电影三级| 99久国产av精品国产电影| 免费看a级黄色片| 国产av一区在线观看免费| 国产一区二区激情短视频| 人妻制服诱惑在线中文字幕| 久久久色成人| 美女免费视频网站| 国产色婷婷99| 欧美日韩乱码在线| 日本-黄色视频高清免费观看| 在线看三级毛片| 亚洲美女视频黄频| 国产激情偷乱视频一区二区| 卡戴珊不雅视频在线播放| 亚洲一区二区三区色噜噜| 欧美丝袜亚洲另类| 久久久久久伊人网av| 国产在线男女| 亚洲五月天丁香| 日韩一本色道免费dvd| 亚洲av电影不卡..在线观看| 亚洲av成人av| 久久久久国产网址| 国产真实乱freesex| 国产一区二区在线观看日韩| 国产精品电影一区二区三区| 日日摸夜夜添夜夜爱| 亚洲成人久久性| 欧美成人精品欧美一级黄| 亚洲av成人av| 又黄又爽又免费观看的视频| 午夜爱爱视频在线播放| 午夜a级毛片| 日日摸夜夜添夜夜添小说| 亚洲国产精品成人综合色| 成人毛片a级毛片在线播放| 成人亚洲精品av一区二区| 99在线人妻在线中文字幕| 国产午夜精品论理片| 男人的好看免费观看在线视频| 亚洲乱码一区二区免费版| 真实男女啪啪啪动态图| 日本免费a在线| 在线天堂最新版资源| 国产精品一二三区在线看| 色吧在线观看| 欧美激情国产日韩精品一区| 日韩亚洲欧美综合| 欧美bdsm另类| 赤兔流量卡办理| videossex国产| 蜜臀久久99精品久久宅男| 国产大屁股一区二区在线视频| 欧美一区二区亚洲| 日本黄大片高清| 成熟少妇高潮喷水视频| 国国产精品蜜臀av免费| 在线观看66精品国产| 国产日本99.免费观看| 国产精品国产三级国产av玫瑰| 成人永久免费在线观看视频| 国内精品宾馆在线| 久久精品国产自在天天线| 国产激情偷乱视频一区二区| 欧美xxxx黑人xx丫x性爽| 欧美又色又爽又黄视频| 国产综合懂色| 一夜夜www| 12—13女人毛片做爰片一| 欧美高清性xxxxhd video| 久久久久久大精品| 九九在线视频观看精品| 免费看日本二区| 精品久久久久久久久av| 亚洲专区国产一区二区| a级毛色黄片| 天堂动漫精品| 国产白丝娇喘喷水9色精品| 亚洲人与动物交配视频| 香蕉av资源在线| 国产亚洲精品av在线| 国内精品宾馆在线| avwww免费| 久久久欧美国产精品| 免费观看的影片在线观看| 国产一区二区在线观看日韩| 中文字幕av在线有码专区| 最近手机中文字幕大全| 一级a爱片免费观看的视频| 成人一区二区视频在线观看| 久久午夜福利片| 亚洲自拍偷在线| 亚洲丝袜综合中文字幕| 欧美成人一区二区免费高清观看| 亚洲美女黄片视频| 国产色婷婷99| 亚洲高清免费不卡视频| 亚洲乱码一区二区免费版| 99热精品在线国产| 亚洲精品久久国产高清桃花| 国产真实伦视频高清在线观看| 美女被艹到高潮喷水动态| 一级毛片我不卡| 精品人妻视频免费看| 国产伦精品一区二区三区四那| 99riav亚洲国产免费| 亚洲精品成人久久久久久| 天堂av国产一区二区熟女人妻| 高清毛片免费观看视频网站| 18禁黄网站禁片免费观看直播| 少妇的逼水好多| 久久午夜亚洲精品久久| av国产免费在线观看| 午夜视频国产福利| 亚洲最大成人手机在线| 久久精品久久久久久噜噜老黄 | 精品日产1卡2卡|