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

    天津地區(qū)用電量與氣溫的非線性響應(yīng)關(guān)系研究

    2021-09-17 08:34:40蘭輝于佳卉曹經(jīng)福劉玉坤孫玫玲黃鶴
    氣象與環(huán)境學報 2021年4期
    關(guān)鍵詞:用電量天津市氣溫

    蘭輝 于佳卉 曹經(jīng)福 劉玉坤 孫玫玲 黃鶴

    (1.天津市氣象服務(wù)中心,天津 300074; 2.天津市氣象科學研究所,天津 300074;3.國網(wǎng)天津市電力公司,天津 300010; 4.天津市氣候中心,天津 300074)

    引言

    準確預(yù)測用電量需求對電力能源供應(yīng)規(guī)劃和經(jīng)濟評估工作具有重要意義,而氣溫作為直接影響夏季降溫和冬季取暖兩大用電需求的氣象要素,是預(yù)測模型中必須考慮的因素[1-3]。在以往的研究中,不同學者利用多種方法分析了氣溫和用電量之間的非線性特征[4-6],確定了用電量與氣溫之間的響應(yīng)關(guān)系,為建立合適的預(yù)測模型奠定基礎(chǔ)。早期研究對于用電需求與氣溫之間的響應(yīng)關(guān)系一般設(shè)定為“V”型[7],即假設(shè)存在一個閾值氣溫,當氣溫超過該值,電能主要用于冷卻降溫,用電量隨氣溫升高而升高(冷卻區(qū)),而氣溫低于該值,電能主要用于加熱取暖,用電量隨氣溫降低而升高(加熱區(qū))。有學者研究發(fā)現(xiàn)分段后中間的閾值氣溫實際上是一個溫度區(qū)間[4,8],在這個溫度區(qū)間內(nèi),電力需求隨氣溫的改變沒有明顯變化,并定義該溫度區(qū)間為舒適區(qū),從而形成了氣溫和用電量之間的“U”型響應(yīng)關(guān)系模型。其中,用電需求與氣溫不同分段區(qū)域的關(guān)系曲線最多的為線性不對稱模型和非線性模型[9]。

    根據(jù)用電量和氣溫的不同響應(yīng)關(guān)系,在用電量預(yù)測時廣泛采用分段預(yù)測的思路,如常用的度日值預(yù)測方法[10-11]。由于分段閾值氣溫的選取會直接影響用電量需求模型的建立,因此,很多學者開展了專項研究。在“V”型響應(yīng)關(guān)系中,通常將閾值氣溫選為18.3 ℃,Bessec和Fouquau[5]研究發(fā)現(xiàn),歐洲15個國家的閾值氣溫存在明顯差異;而在“U”型響應(yīng)關(guān)系中,受社會環(huán)境、經(jīng)濟結(jié)構(gòu)、建筑節(jié)能等因素的可能影響,不同地區(qū)的閾值氣溫也存在差異[10-13]。

    雖然以往研究取得了大量成果,但大部分關(guān)于用電量與氣象條件的非線性關(guān)系研究僅根據(jù)用電量預(yù)測的需要,對單一模型進行分析。國內(nèi)則一般直接引入已有的閾值氣溫研究結(jié)果,而缺乏對本地適用性的分析。在閾值氣溫的計算中,多僅根據(jù)氣溫與用電量的響應(yīng)關(guān)系確定,而未考慮其他氣象要素的影響,也未考慮同一地區(qū)不同時段閾值氣溫的差異,從而導(dǎo)致閾值氣溫的計算結(jié)果產(chǎn)生偏差,進而影響用電量需求預(yù)測的準確率。針對以上問題,本文利用計量經(jīng)濟學上普遍應(yīng)用的非線性滑動回歸(STR)模型和門限回歸(TR)模型,分析了“V”型非線性、“U”型線性和“U”型非線性3種模型下氣溫與用電量的響應(yīng)關(guān)系,計算不同模型對應(yīng)的閾值氣溫,并在此基礎(chǔ)上,分析相對濕度對閾值氣溫的影響及同一區(qū)域閾值氣溫隨時間的變化,為進一步開展氣溫累積效應(yīng)和用電量需求分析等提供參考。

    1 資料與方法

    1.1 資料來源及處理

    所用數(shù)據(jù)包括電力數(shù)據(jù)和氣象數(shù)據(jù)。電力數(shù)據(jù)來自國網(wǎng)天津市電力公司,為天津市2002—2005年和2014—2019年逐日用電量數(shù)據(jù);氣象數(shù)據(jù)來自天津市氣象局,為2002—2005年和2014—2019年天津市城市氣候監(jiān)測站(站號54517)的逐日平均氣溫、相對濕度及1960—2019年年平均氣溫數(shù)據(jù),其中2014—2019年的電力和氣象數(shù)據(jù)用于開展用電量與氣象要素分析,2002—2005年電力和氣象數(shù)據(jù)用于評估閾值氣溫隨時間變化,1960—2019年年平均氣溫主要用于分析閾值氣溫改變的原因。由于國外多個研究[4-5,8]已論證非線性滑動回歸(STR)模型和門限回歸(TR)模型在確定閾值氣溫方面的適用性,因此,本文未選取檢驗樣本開展評估工作。

    為評估用電量與氣象要素的非線性關(guān)系,需要剔除非氣象因素對用電量產(chǎn)生的影響。根據(jù)Scapin等[11]的模型,用電量除受氣象條件影響外,還表現(xiàn)出明顯的周周期特點和長期趨勢特征,模型具體公式為:

    (1)

    (2)

    式(1)中,D為日用電量;t為日數(shù);I為周周期的虛擬變量;V為氣象因子;α、β和γ分別為系數(shù);ε為殘差。式(1)右側(cè)第1項為用電量時間序列的趨勢項,主要反映經(jīng)濟形勢和消費習慣的長期變化對用電量的影響;右側(cè)第2項代表用電量的周期性變化。有研究表明[4],周二至周四用電量比較接近,周五略偏少,周一次之,周末最少,將I1在周一時設(shè)置為1,其他時間設(shè)置為0,以此類推,I2、I3、I4、I5分別在周二至周四、周五、周六、周日時設(shè)置為1,其他時間設(shè)置為0;右側(cè)第3項代表氣象條件影響下的用電量變化,主要受氣溫條件影響,由于最高、最低氣溫更適合預(yù)測日用電量的峰值和低谷[14],因此本文采用平均氣溫作為預(yù)測因子用于預(yù)測日總用電量;右側(cè)第4項為殘差項。

    采用二次曲線擬合用電量隨時間變化趨勢(右側(cè)第1項),并僅選取周二至周四的數(shù)據(jù),以消除周周期變化(右側(cè)第2項)對分析結(jié)果的影響。

    1.2 用電量與氣溫關(guān)系的模型

    1.2.1 TR模型

    Hansen在門限回歸模型(Threshold Regression Model,TR)上進行了大量研究[15-16]。先后提出時間序列門限自回歸模型的估計和檢驗方法、具有個體效應(yīng)的面板門限模型的計量分析方法、多元門限回歸模型的思想等一系列創(chuàng)新性研究,并發(fā)展出一套完備的門限回歸模型,其中Hansen的雙閾值單元線性回歸模型的計算公式為:

    (3)

    式(3)中,DF為氣象用電量;MT為平均氣溫;μ為常數(shù);ω為系數(shù);ε為殘差;t為時間;Th1和Th2分別為加熱區(qū)與舒適區(qū)、舒適區(qū)與冷卻區(qū)的閾值氣溫。

    Moral-Carcedo等[4]利用式(3)計算得到用電量隨氣溫變化的雙閾值氣溫,通過建立3個線性模型,計算不同模型之間的交點(Th1和Th2),即上下限的閾值氣溫。式(3)中假設(shè)加熱區(qū)和冷卻區(qū)的用電量與氣溫之間呈線性變化,對其進行調(diào)整,可以得到用電量和氣溫非線性響應(yīng)關(guān)系下的雙閾值氣溫,公式如下:

    (4)

    式(4)中,k為系數(shù);其他各變量與式(3)一致。

    1.2.2 STR模型

    平滑轉(zhuǎn)換自回歸模型(Smooth Transition Auto-regression Model,STAR)是由GRANGER等提出的,目前已經(jīng)形成了一套較為成熟的建模程序,包括模型的設(shè)定、估計和診斷測試,STAR模型的提出為非線性滑動回歸模型(Smooth Transition Regression Model,STR)奠定了理論基礎(chǔ)[17]。STR模型是對STAR模型的改進和發(fā)展,為模擬電力需求對溫度的響應(yīng)而提出的STR模型[4]可以表示為:

    DFt=(α1+β1MTt)[1-G(MTt;γ,c)]+(α2+β2MTt)G(MTt;γ,c)+εt

    (5)

    式(5)中,εt為獨立分布的誤差,假設(shè)其服從正態(tài)分布;轉(zhuǎn)換函數(shù)G(MTt;γ,c) 是一個取值范圍為[0,1]的連續(xù)函數(shù),常見的形式是logistic函數(shù)形式和指數(shù)函數(shù)形式,其中 logistic 函數(shù)形式的表達式為:

    G(MTt;γ,c)=[1+exp{-γ(MTt-c)}]-1

    (6)

    式(6)中,MT表示平均氣溫;斜率γ反映了模型描述的狀態(tài)“0”過渡到狀態(tài)“1”的速度大小及轉(zhuǎn)換的平滑性;c表示常數(shù);t表示日數(shù)。式(5)和式(6)共同構(gòu)成了logistic平滑回歸模型,即LSTR模型。具體計算過程可參考昌春艷等[17]的研究成果。

    1.2.3 模型評價方法

    Hansen[15]以殘差平方和最小化為條件確定門限值,本文采用類似的均方根誤差確定最優(yōu)閾值,均方根誤差被廣泛應(yīng)用于模型檢驗[18-19],公式如下:

    (7)

    分析結(jié)果時,出現(xiàn)兩個臨近閾值均方差相同的情況,本文采用決定系數(shù)輔助選取最優(yōu)閾值,同時決定系數(shù)可用于對比不同時間段均方差較大的模型效果評估,公式如下[20-21]:

    (8)

    2 結(jié)果分析

    2.1 用電量與氣溫關(guān)系模型

    利用STR模型計算的2014—2019年天津市逐日用電量與氣溫“V”型非線性擬合曲線(圖1),利用TR模型計算的“U”型線性擬合和非線性擬合曲線(圖2),由于數(shù)據(jù)特征明顯不符合“V”型線性擬合曲線,因此未計算該種情況下的閾值氣溫。分析可以發(fā)現(xiàn),“V”型非線性響應(yīng)關(guān)系下的閾值氣溫為18.8 ℃;“U”型響應(yīng)關(guān)系下,線性和非線性模型的舒適區(qū)范圍存在明顯差異,其中線性曲線模型下的舒適區(qū)范圍為12.3—23.4 ℃,非線性模型下的舒適區(qū)范圍則明顯變小,為13.7—21.7 ℃,這主要是由于非線性曲線可以捕捉用電量隨氣溫增長幅度較小時的變化區(qū)間,而采用線性曲線,該區(qū)間被全部并入舒適區(qū)。

    圖1 2014—2019年天津市逐日用電量與氣溫“V”型非線性擬合曲線以及閾值氣溫Fig.1 Daily electricity consumption and temperature "V"-shaped nonlinear fitting curve as well as threshold temperature in Tianjin area from 2014 to 2019

    圖2 2014—2019年天津市逐日用電量與氣溫“U”型擬合曲線以及閾值氣溫Fig.2 Daily electricity consumption and temperature "U"-shaped fitting curve as well as threshold temperature in Tianjin area from 2014 to 2019

    由圖2可知,冷卻區(qū)用電量隨氣溫升高的速度明顯大于加熱區(qū),在“U”型線性曲線中,加熱區(qū)用電量隨氣溫的增長率為-2.6 GWh/℃,冷卻區(qū)用電量隨氣溫的增長率為8.0 GWh/℃,同時冷卻區(qū)的用電量峰值也明顯偏大,這主要是由于夏季制冷主要由空調(diào)用電提供,而冬季加熱則是由電力和天然氣組合提供。除此之外,可以看出相較于舒適區(qū)與冷卻區(qū)之間比較連貫的過渡,擬合曲線在“U”型線性關(guān)系模型下,加熱區(qū)與舒適區(qū)之間的過渡并不連貫,這是由于天津市采用集體供暖的方式,在供暖季啟停日期的前后,室外溫度的變化不大,但用電量會出現(xiàn)明顯變化,而該段時期的平均氣溫一般位于12.3 ℃左右,從而造成擬合曲線在舒適區(qū)和加熱區(qū)之間的不連貫,而舒適區(qū)和冷卻區(qū)之間的用電量主要受工業(yè)、商業(yè)和居民自發(fā)的用電行為影響,總體變化比較規(guī)律,因此擬合曲線在閾值氣溫處的用電量基本一致。

    表1為各模型計算的2014—2019年天津市逐日用電量閾值氣溫及對應(yīng)的預(yù)測效果,對比3種模型的預(yù)測效果發(fā)現(xiàn),單閾值的非線性模型優(yōu)于雙閾值的線性模型,而雙閾值的非線性模型優(yōu)于單閾值的非線性模型,表明雙閾值優(yōu)于單閾值,非線性優(yōu)于線性。

    表1 2014—2019年天津市逐日用電量與氣溫在不同響應(yīng)關(guān)系下的閾值氣溫及預(yù)測效果Table 1 Threshold temperature and prediction effect under different response relationships between daily electricity consumption and temperature in Tianjin area from 2014 to 2019

    2.2 閾值氣溫的影響要素

    2.2.1 相對濕度對閾值氣溫的影響

    有研究表明,除了氣溫外,相對濕度對用電量的影響也十分顯著[22-23],有研究從體感溫度或相對濕度角度分析閾值氣溫變化[7-8]。本文將所有數(shù)據(jù)按相對濕度范圍分為3組,分別為30%—50%、50%—70%和其他,計算在“U”非線型模型(利用TR模型)中前兩組的閾值氣溫,以評估相對濕度對閾值氣溫的可能影響。

    2014—2019年天津市逐日用電量與氣溫擬合曲線在不同相對濕度范圍內(nèi)對應(yīng)的閾值氣溫見圖3,從圖3可以看出,不同相對濕度下的舒適區(qū)存在明顯差異,相對濕度在30%—50%時的舒適區(qū)范圍明顯大于相對濕度在50%—70%時的舒適區(qū)范圍,特別是舒適區(qū)與冷卻區(qū)的閾值氣溫,相差2.2 ℃,表明相對濕度對舒適區(qū)和冷卻區(qū)的閾值氣溫影響較大,對舒適區(qū)和加熱區(qū)閾值氣溫的影響則并不明顯。

    圖3 2014—2019年天津市逐日用電量與氣溫擬合曲線在不同相對濕度范圍內(nèi)對應(yīng)的閾值氣溫Fig.3 The threshold temperature corresponding to the daily electricity consumption and temperature fitting curve in different relative humidity ranges in Tianjin area from 2014 to 2019

    值得注意的是,在冷卻區(qū)中,當氣溫相同時,相對濕度為50%—70%時的用電量基本大于相對濕度為30%—50%時的用電量,而這一特點在舒適區(qū)和加熱區(qū)不存在,用電量最大的日期也基本出現(xiàn)在相對濕度為50%—70%時,兩種相對濕度范圍內(nèi),冷卻區(qū)內(nèi)的用電量隨氣溫變化曲線基本一致,這主要是二者在舒適區(qū)與冷卻區(qū)的閾值氣溫差異較大,從而造成相同氣溫時,相對濕度為50%—70%時的用電量均高于相對濕度為30%—50%時。

    由于日平均相對濕度低于30%及大于70%的樣本數(shù)量偏少,本文未做詳細分析。以往研究表明[22-24],隨著相對濕度增加,冷卻用電量(或負荷)將明顯上升,結(jié)合本文的研究成果,可以推測在相對濕度超過70%情況下,冷卻區(qū)與舒適區(qū)的閾值氣溫將低于19.4 ℃,反之在相對濕度低于30%情況下,冷卻區(qū)與舒適區(qū)的閾值氣溫將高于21.6 ℃。

    2.2.2 閾值氣溫隨時間變化特征

    有研究表明,氣候變化會對建筑采暖降溫負荷產(chǎn)生影響[25],為了探究氣候變化可能對閾值氣溫的影響程度,選取2002—2005年的天津市逐日數(shù)據(jù)做相同的統(tǒng)計分析,以評估閾值氣溫隨時間的變化,如表2所示。

    表2 2002—2005年和2014—2019年天津市逐日用電量與氣溫在不同響應(yīng)關(guān)系下閾值氣溫和預(yù)測效果對比Table 2 Comparison of threshold air temperature and prediction effect under different response relationships between daily electricity consumption and air temperature in Tianjin area between 2002-2005 and 2014-2019

    由表2可知,在不同響應(yīng)關(guān)系中,閾值氣溫均隨時間發(fā)生了變化,利用STR模型計算“V”型非線性模型下的閾值氣溫從18.2 ℃增長到18.8 ℃,表明居民采取加熱和冷卻措施的氣溫向高溫方向偏移。由圖4可知,氣溫整體的變暖趨勢非常明顯,決定系數(shù)為0.825,每年增溫幅度高達0.0689 ℃,按此增幅計算2002—2005年和2014—2019年的增幅約為0.9 ℃,與閾值氣溫變化幅度比較接近,這可能是由于人體適應(yīng)環(huán)境的能力,造成舒適中值隨全球變暖有所提升。

    圖4 1960—2019年天津市年平均氣溫變化曲線及趨勢Fig.4 Annual average temperature change curve and trend in Tianjin area from 1960 to 2019

    “U”型非線性模型在不同時間下的計算結(jié)果表明,2014—2019年的舒適區(qū)范圍(13.7—21.7 ℃)較2002—2005年的范圍(13.6—23.0 ℃)明顯縮小,這主要是因為天津經(jīng)濟高速發(fā)展過程中,居民和商業(yè)調(diào)溫設(shè)備大量增加,且人們對環(huán)境舒適的要求更高,而基本電價未發(fā)生明顯改變,人們普遍會過早的采取加熱或者降溫措施,特別是降溫措施(降低了1.3 ℃)。

    同時,天津經(jīng)濟的快速發(fā)展、氣候的持續(xù)變暖和居民用電習慣的改變還造成每年用電量的明顯提升(2019年用電量為2002年用電量的3倍多),而節(jié)能政策的推出與成果應(yīng)用及相關(guān)產(chǎn)業(yè)的規(guī)?;l(fā)展,使得GDP的增速(2019年GDP為2002年的6倍多)明顯大于用電量的增長,單位GDP需要的能耗明顯減少。

    3 結(jié)論與討論

    (1)本文運用計量經(jīng)濟學的兩種方法計算了2014—2019年天津市用電量和氣溫在不同響應(yīng)關(guān)系下的閾值氣溫,發(fā)現(xiàn)閾值氣溫存在明顯差異,就日用電量預(yù)測效果而言,“U”型優(yōu)于“V”型模型,非線性模型優(yōu)于線性模型。

    (2)相對濕度可對閾值氣溫產(chǎn)生較大影響,舒適區(qū)范圍會隨相對濕度增大明顯減小,特別是舒適區(qū)和冷卻區(qū)交界處的閾值氣溫,會隨相對濕度增大向低溫區(qū)偏移。

    (3)對比2002—2005年和2014—2019年兩段時間天津市用電量與氣溫不同響應(yīng)關(guān)系下的閾值氣溫,發(fā)現(xiàn)閾值氣溫存在明顯差異,這可能是由于經(jīng)濟增長和氣候變化影響了居民用電習慣,從而引起閾值氣溫發(fā)生改變。

    (4)在實際開展用電量和氣象條件相關(guān)分析和預(yù)測時,應(yīng)根據(jù)相對濕度范圍分別計算閾值氣溫;當分析或者預(yù)測的數(shù)據(jù)時段較長,應(yīng)分段選取不同的閾值氣溫。本文僅以天津為例,確定了用電量和氣溫在不同模型下的閾值氣溫,并分析了相對濕度和時間變化對閾值氣溫產(chǎn)生的影響,并未開展其他地區(qū)的對比研究。因此,需更多地區(qū)共同開展類似研究,從而發(fā)現(xiàn)相對濕度和時間變化對閾值氣溫影響更為普適性的規(guī)律。

    猜你喜歡
    用電量天津市氣溫
    02 國家能源局:1~7月全社會用電量同比增長3.4%
    01 國家能源局:3月份全社會用電量同比增長3.5%
    基于FY-3D和FY-4A的氣溫時空融合
    有趣的云
    深冬氣溫多變 蔬菜管理要隨機應(yīng)變
    如果畫筆會說話?
    天津市第三屆“未來之星”初中數(shù)學邀請賽
    1~10月全社會用電量累計56552億千瓦時同比增長8.7%
    與氣溫成反比的東西
    2014年全社會用電量
    露出奶头的视频| 在线观看午夜福利视频| 一级作爱视频免费观看| 18禁裸乳无遮挡动漫免费视频| 熟女少妇亚洲综合色aaa.| 91精品国产国语对白视频| 变态另类成人亚洲欧美熟女 | 高潮久久久久久久久久久不卡| 高清欧美精品videossex| 一进一出抽搐gif免费好疼 | 十八禁人妻一区二区| 18禁黄网站禁片午夜丰满| 久久影院123| 高清欧美精品videossex| 亚洲av成人不卡在线观看播放网| 欧美黑人欧美精品刺激| 亚洲国产欧美一区二区综合| 99久久综合精品五月天人人| 国产亚洲精品一区二区www | 久久久国产一区二区| 黄色a级毛片大全视频| 美国免费a级毛片| 国产欧美日韩一区二区三区在线| 高清视频免费观看一区二区| 麻豆av在线久日| 成年女人毛片免费观看观看9 | 日本精品一区二区三区蜜桃| 又大又爽又粗| 人人澡人人妻人| 欧美激情久久久久久爽电影 | 一二三四在线观看免费中文在| 欧美av亚洲av综合av国产av| 亚洲av电影在线进入| 91国产中文字幕| 丰满人妻熟妇乱又伦精品不卡| 久久精品亚洲精品国产色婷小说| 99国产极品粉嫩在线观看| 亚洲久久久国产精品| 久久精品成人免费网站| 伦理电影免费视频| 国产精品99久久99久久久不卡| 欧美av亚洲av综合av国产av| 精品国产乱码久久久久久男人| 99久久国产精品久久久| 久久青草综合色| 久热爱精品视频在线9| 欧美在线黄色| 久久精品国产亚洲av高清一级| 一进一出抽搐动态| 精品国产乱子伦一区二区三区| 亚洲国产精品合色在线| 国产一区在线观看成人免费| 午夜亚洲福利在线播放| 很黄的视频免费| 日日爽夜夜爽网站| 午夜精品国产一区二区电影| 日本精品一区二区三区蜜桃| 欧美成狂野欧美在线观看| 久久午夜综合久久蜜桃| 精品少妇一区二区三区视频日本电影| 久久精品国产清高在天天线| 欧美在线一区亚洲| 正在播放国产对白刺激| 老汉色∧v一级毛片| 91麻豆精品激情在线观看国产 | 午夜成年电影在线免费观看| 色婷婷久久久亚洲欧美| 村上凉子中文字幕在线| 精品国产一区二区三区久久久樱花| 中文字幕人妻熟女乱码| 色综合婷婷激情| 午夜精品国产一区二区电影| 亚洲av成人av| 国产色视频综合| 欧美久久黑人一区二区| 成人av一区二区三区在线看| 99国产精品免费福利视频| 国产色视频综合| 欧美大码av| 侵犯人妻中文字幕一二三四区| 精品人妻在线不人妻| 最近最新中文字幕大全免费视频| 亚洲综合色网址| 中文字幕精品免费在线观看视频| 巨乳人妻的诱惑在线观看| 99精国产麻豆久久婷婷| 777米奇影视久久| 国产精品久久久久成人av| 国产日韩欧美亚洲二区| 亚洲av熟女| 女人被狂操c到高潮| 丝袜美足系列| 757午夜福利合集在线观看| 丝瓜视频免费看黄片| 国产欧美亚洲国产| 99久久精品国产亚洲精品| 午夜91福利影院| 亚洲伊人色综图| 国产av精品麻豆| 女人精品久久久久毛片| 一边摸一边抽搐一进一出视频| 在线观看一区二区三区激情| 性色av乱码一区二区三区2| 99久久精品国产亚洲精品| 日韩免费高清中文字幕av| av天堂在线播放| 国产精品自产拍在线观看55亚洲 | 少妇被粗大的猛进出69影院| 日韩大码丰满熟妇| 丁香六月欧美| 无遮挡黄片免费观看| 叶爱在线成人免费视频播放| 亚洲第一欧美日韩一区二区三区| 一边摸一边抽搐一进一出视频| 亚洲精品久久成人aⅴ小说| 一级毛片女人18水好多| 老熟女久久久| 中文字幕另类日韩欧美亚洲嫩草| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品国产高清国产av | 黄色丝袜av网址大全| 国产精品一区二区精品视频观看| 黄片小视频在线播放| 麻豆国产av国片精品| 久久青草综合色| 大型黄色视频在线免费观看| 黑丝袜美女国产一区| 十分钟在线观看高清视频www| 午夜福利影视在线免费观看| 久久狼人影院| 别揉我奶头~嗯~啊~动态视频| 窝窝影院91人妻| 纯流量卡能插随身wifi吗| 成年人免费黄色播放视频| 操出白浆在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 下体分泌物呈黄色| 99国产精品免费福利视频| bbb黄色大片| 午夜福利影视在线免费观看| 少妇被粗大的猛进出69影院| 欧美日韩成人在线一区二区| 久久国产亚洲av麻豆专区| 乱人伦中国视频| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲欧美精品永久| 不卡一级毛片| 99国产综合亚洲精品| 久久精品国产99精品国产亚洲性色 | 男女免费视频国产| 国内久久婷婷六月综合欲色啪| 中文亚洲av片在线观看爽 | 一本综合久久免费| 国产一区在线观看成人免费| 久久亚洲精品不卡| 狠狠狠狠99中文字幕| 人妻久久中文字幕网| 亚洲色图av天堂| 欧美黄色片欧美黄色片| 日本vs欧美在线观看视频| 欧美国产精品va在线观看不卡| 80岁老熟妇乱子伦牲交| 美女 人体艺术 gogo| 国产精品久久久久成人av| 热99久久久久精品小说推荐| bbb黄色大片| 高清欧美精品videossex| 国产精品国产av在线观看| 久久九九热精品免费| 亚洲午夜理论影院| 另类亚洲欧美激情| 久久久精品国产亚洲av高清涩受| 老熟妇乱子伦视频在线观看| 天堂动漫精品| 自拍欧美九色日韩亚洲蝌蚪91| 精品少妇一区二区三区视频日本电影| 午夜福利免费观看在线| 中出人妻视频一区二区| 999久久久精品免费观看国产| 国产精品九九99| 亚洲成a人片在线一区二区| 无限看片的www在线观看| 啦啦啦免费观看视频1| av超薄肉色丝袜交足视频| 9色porny在线观看| 欧美精品高潮呻吟av久久| 精品福利永久在线观看| 91精品国产国语对白视频| videos熟女内射| 久久亚洲真实| 91精品国产国语对白视频| 亚洲欧美色中文字幕在线| 亚洲精品自拍成人| 精品卡一卡二卡四卡免费| 巨乳人妻的诱惑在线观看| 欧美亚洲 丝袜 人妻 在线| 国精品久久久久久国模美| 国产日韩欧美亚洲二区| 男人操女人黄网站| 亚洲色图av天堂| 中文字幕人妻丝袜一区二区| 成人国语在线视频| 悠悠久久av| 中文字幕人妻丝袜制服| 久久人人爽av亚洲精品天堂| 嫩草影视91久久| 叶爱在线成人免费视频播放| av一本久久久久| 国产精品久久久久久精品古装| 俄罗斯特黄特色一大片| 久久久久久久久免费视频了| 午夜影院日韩av| 国产成人影院久久av| 在线国产一区二区在线| 18禁国产床啪视频网站| 成人国语在线视频| 国产高清视频在线播放一区| а√天堂www在线а√下载 | 性色av乱码一区二区三区2| 久久午夜综合久久蜜桃| 欧美在线黄色| 国产精品秋霞免费鲁丝片| 国产精品久久久久成人av| 亚洲精品一二三| 免费观看a级毛片全部| 色播在线永久视频| 免费在线观看完整版高清| 免费在线观看黄色视频的| 国产成+人综合+亚洲专区| 国产精品免费一区二区三区在线 | 悠悠久久av| 大型av网站在线播放| 国产精品一区二区免费欧美| 国产精品久久久av美女十八| 欧美乱妇无乱码| 午夜福利影视在线免费观看| www.999成人在线观看| 18禁裸乳无遮挡动漫免费视频| 黄网站色视频无遮挡免费观看| 国产一区二区三区在线臀色熟女 | 亚洲欧美日韩另类电影网站| 免费高清在线观看日韩| 免费观看精品视频网站| 丁香欧美五月| 成人黄色视频免费在线看| 少妇的丰满在线观看| 国产精品一区二区精品视频观看| 亚洲性夜色夜夜综合| 丰满人妻熟妇乱又伦精品不卡| av中文乱码字幕在线| 这个男人来自地球电影免费观看| 精品国产国语对白av| 久久中文看片网| 成人国语在线视频| 国产高清视频在线播放一区| 中亚洲国语对白在线视频| 亚洲中文日韩欧美视频| 一本一本久久a久久精品综合妖精| 精品无人区乱码1区二区| 变态另类成人亚洲欧美熟女 | 美女高潮到喷水免费观看| 777米奇影视久久| 国产精品综合久久久久久久免费 | 亚洲欧美色中文字幕在线| 18在线观看网站| 国产在视频线精品| 国产精品九九99| 男女高潮啪啪啪动态图| 亚洲欧美色中文字幕在线| 久久香蕉激情| 曰老女人黄片| 欧美激情高清一区二区三区| 中文字幕av电影在线播放| 亚洲片人在线观看| 国产国语露脸激情在线看| 国产精品免费视频内射| 国产亚洲精品久久久久5区| 黄片播放在线免费| 一进一出好大好爽视频| 国产亚洲欧美精品永久| 好男人电影高清在线观看| 亚洲精华国产精华精| 99久久99久久久精品蜜桃| 两个人免费观看高清视频| 99热网站在线观看| 伊人久久大香线蕉亚洲五| 在线观看午夜福利视频| 黄色成人免费大全| 亚洲精品中文字幕一二三四区| 看片在线看免费视频| 他把我摸到了高潮在线观看| 亚洲综合色网址| 99久久人妻综合| 波多野结衣一区麻豆| 欧美精品亚洲一区二区| 18禁美女被吸乳视频| 亚洲成人手机| 亚洲精品国产精品久久久不卡| av一本久久久久| 天天影视国产精品| 亚洲免费av在线视频| 老司机靠b影院| av网站免费在线观看视频| 91在线观看av| av有码第一页| 精品少妇久久久久久888优播| 国产蜜桃级精品一区二区三区 | 18禁国产床啪视频网站| 男女高潮啪啪啪动态图| 午夜福利乱码中文字幕| 久久久久视频综合| 99久久人妻综合| 精品熟女少妇八av免费久了| 欧美 日韩 精品 国产| 又大又爽又粗| 中文亚洲av片在线观看爽 | 老司机靠b影院| 大香蕉久久成人网| 国产精品一区二区精品视频观看| 亚洲av欧美aⅴ国产| 成人永久免费在线观看视频| 成人亚洲精品一区在线观看| 亚洲精品av麻豆狂野| 中文字幕人妻丝袜一区二区| 夜夜躁狠狠躁天天躁| 国产高清激情床上av| 欧美丝袜亚洲另类 | 丝袜在线中文字幕| 如日韩欧美国产精品一区二区三区| 欧美精品亚洲一区二区| 久久久久久久久久久久大奶| www日本在线高清视频| 日韩视频一区二区在线观看| 男男h啪啪无遮挡| 99在线人妻在线中文字幕 | 国产男女超爽视频在线观看| xxxhd国产人妻xxx| 亚洲精品粉嫩美女一区| 免费一级毛片在线播放高清视频 | 18禁国产床啪视频网站| 两个人看的免费小视频| 99re6热这里在线精品视频| 久久久久精品国产欧美久久久| 99re在线观看精品视频| 亚洲国产中文字幕在线视频| 国产xxxxx性猛交| 少妇的丰满在线观看| 午夜福利视频在线观看免费| 女人被狂操c到高潮| 新久久久久国产一级毛片| 99久久国产精品久久久| 婷婷精品国产亚洲av在线 | 久热这里只有精品99| 777米奇影视久久| 久久精品国产a三级三级三级| 黄色片一级片一级黄色片| 国产欧美日韩一区二区三区在线| 欧美日韩一级在线毛片| 自线自在国产av| 欧美日韩一级在线毛片| www.自偷自拍.com| 免费在线观看亚洲国产| 波多野结衣av一区二区av| av片东京热男人的天堂| 国产亚洲精品第一综合不卡| 日本wwww免费看| 国产精品影院久久| 国产成人免费观看mmmm| 天天影视国产精品| 国产成人啪精品午夜网站| 国产淫语在线视频| 多毛熟女@视频| av欧美777| 成熟少妇高潮喷水视频| 精品人妻在线不人妻| 男女高潮啪啪啪动态图| 亚洲av美国av| 亚洲在线自拍视频| 国产一区二区三区综合在线观看| 日日摸夜夜添夜夜添小说| 99re6热这里在线精品视频| 高潮久久久久久久久久久不卡| 亚洲久久久国产精品| 国产在线精品亚洲第一网站| 亚洲熟妇中文字幕五十中出 | 精品国产亚洲在线| 国产色视频综合| 激情视频va一区二区三区| 久久久久国内视频| 热re99久久精品国产66热6| 国产av又大| 免费在线观看影片大全网站| 精品国产一区二区久久| 免费黄频网站在线观看国产| 亚洲精品国产一区二区精华液| 啦啦啦免费观看视频1| 欧美一级毛片孕妇| 国产淫语在线视频| 久久影院123| 男女床上黄色一级片免费看| 美女福利国产在线| 叶爱在线成人免费视频播放| a级毛片黄视频| 精品福利永久在线观看| 亚洲三区欧美一区| 别揉我奶头~嗯~啊~动态视频| 午夜精品国产一区二区电影| 女人精品久久久久毛片| avwww免费| 欧美国产精品一级二级三级| 99热只有精品国产| 99精品在免费线老司机午夜| 大陆偷拍与自拍| 老司机深夜福利视频在线观看| 亚洲av日韩在线播放| 国产精品一区二区精品视频观看| 欧美乱妇无乱码| 国产精品综合久久久久久久免费 | 成人影院久久| 免费在线观看完整版高清| xxxhd国产人妻xxx| 天天添夜夜摸| 亚洲午夜精品一区,二区,三区| 麻豆av在线久日| 动漫黄色视频在线观看| 无人区码免费观看不卡| 一级a爱片免费观看的视频| www.精华液| 久久中文看片网| 国产99久久九九免费精品| 国产三级黄色录像| 国产精品久久久久成人av| 啪啪无遮挡十八禁网站| 久久久精品区二区三区| 一区在线观看完整版| 老熟妇仑乱视频hdxx| 99热网站在线观看| 99re在线观看精品视频| 国产日韩一区二区三区精品不卡| 国产精品.久久久| 女人爽到高潮嗷嗷叫在线视频| 满18在线观看网站| 精品少妇一区二区三区视频日本电影| 亚洲一卡2卡3卡4卡5卡精品中文| 精品亚洲成国产av| 成人18禁在线播放| 国产精品一区二区免费欧美| 亚洲久久久国产精品| 在线观看免费高清a一片| 大陆偷拍与自拍| 国产精品国产高清国产av | av天堂在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 丰满的人妻完整版| 国产激情欧美一区二区| 黄色视频,在线免费观看| 欧美激情久久久久久爽电影 | 夜夜爽天天搞| 美女视频免费永久观看网站| 国产97色在线日韩免费| 两个人免费观看高清视频| 国产成人av激情在线播放| 一个人免费在线观看的高清视频| 老司机亚洲免费影院| 在线观看午夜福利视频| 亚洲 国产 在线| 亚洲第一欧美日韩一区二区三区| 久久精品国产清高在天天线| 成人精品一区二区免费| 天天操日日干夜夜撸| 动漫黄色视频在线观看| 亚洲成a人片在线一区二区| 精品一品国产午夜福利视频| 免费观看精品视频网站| 狠狠婷婷综合久久久久久88av| 下体分泌物呈黄色| 久久香蕉激情| 国产在线观看jvid| 久久国产精品人妻蜜桃| 午夜福利在线免费观看网站| 激情在线观看视频在线高清 | av天堂在线播放| 免费久久久久久久精品成人欧美视频| 亚洲精品美女久久久久99蜜臀| 精品国产国语对白av| 国产精品久久久久久人妻精品电影| 丰满饥渴人妻一区二区三| 亚洲人成伊人成综合网2020| 欧美黑人精品巨大| 亚洲精品av麻豆狂野| 男女之事视频高清在线观看| 国产精品1区2区在线观看. | 国产99白浆流出| 亚洲人成伊人成综合网2020| 国产av一区二区精品久久| 亚洲精品久久午夜乱码| 午夜亚洲福利在线播放| 久久久精品区二区三区| 最新在线观看一区二区三区| 他把我摸到了高潮在线观看| 亚洲中文日韩欧美视频| a级毛片在线看网站| 中文字幕制服av| 精品午夜福利视频在线观看一区| 亚洲成av片中文字幕在线观看| 99久久99久久久精品蜜桃| 日韩人妻精品一区2区三区| 少妇粗大呻吟视频| 色综合欧美亚洲国产小说| 日本精品一区二区三区蜜桃| 久久精品亚洲精品国产色婷小说| 黄片大片在线免费观看| 国产精品影院久久| 亚洲精品国产色婷婷电影| 久久久国产精品麻豆| 久久久久国产精品人妻aⅴ院 | 无遮挡黄片免费观看| 制服人妻中文乱码| 俄罗斯特黄特色一大片| 建设人人有责人人尽责人人享有的| 老司机午夜十八禁免费视频| 另类亚洲欧美激情| 亚洲国产欧美网| 亚洲av欧美aⅴ国产| 中文字幕制服av| 高清av免费在线| 欧美黑人欧美精品刺激| 不卡av一区二区三区| 制服诱惑二区| 搡老岳熟女国产| 亚洲欧美日韩高清在线视频| 久久精品国产亚洲av高清一级| 欧美日韩黄片免| 日本vs欧美在线观看视频| 午夜免费观看网址| 女人爽到高潮嗷嗷叫在线视频| 十八禁人妻一区二区| 男女午夜视频在线观看| 久久人人爽av亚洲精品天堂| 天堂中文最新版在线下载| 999久久久精品免费观看国产| 女人精品久久久久毛片| 午夜福利,免费看| 免费一级毛片在线播放高清视频 | 免费女性裸体啪啪无遮挡网站| 精品国产超薄肉色丝袜足j| 国产成人欧美在线观看 | 午夜91福利影院| 十八禁人妻一区二区| 91成人精品电影| 黑人欧美特级aaaaaa片| 久久人妻熟女aⅴ| 国产成人欧美| 国产激情久久老熟女| 亚洲中文av在线| 亚洲专区国产一区二区| 很黄的视频免费| 亚洲成人国产一区在线观看| 又黄又粗又硬又大视频| 69精品国产乱码久久久| av在线播放免费不卡| 在线观看免费高清a一片| 午夜福利在线观看吧| 老熟女久久久| 久久人妻av系列| 一夜夜www| 午夜视频精品福利| 免费av中文字幕在线| 免费高清在线观看日韩| 久久久精品区二区三区| 亚洲综合色网址| 国产精品.久久久| 一级黄色大片毛片| 久久久久久亚洲精品国产蜜桃av| 天堂动漫精品| 亚洲色图av天堂| 日韩一卡2卡3卡4卡2021年| 亚洲成人国产一区在线观看| 99国产精品一区二区三区| 免费观看a级毛片全部| 久久亚洲精品不卡| 午夜91福利影院| 欧美乱妇无乱码| 精品视频人人做人人爽| 亚洲国产欧美网| 中文欧美无线码| 巨乳人妻的诱惑在线观看| 国产1区2区3区精品| 99riav亚洲国产免费| 午夜福利视频在线观看免费| av国产精品久久久久影院| 国产精品久久视频播放| 免费黄频网站在线观看国产| 国产99久久九九免费精品| 久久天躁狠狠躁夜夜2o2o| 日本五十路高清| 人妻丰满熟妇av一区二区三区 | 99国产精品一区二区三区| 操出白浆在线播放| 国产不卡av网站在线观看| 久久精品亚洲精品国产色婷小说| 别揉我奶头~嗯~啊~动态视频| 少妇被粗大的猛进出69影院| 高清欧美精品videossex| 欧美日韩福利视频一区二区| 久热这里只有精品99| 国产成+人综合+亚洲专区| 亚洲人成电影免费在线| 国产免费av片在线观看野外av| 欧美乱色亚洲激情| 久99久视频精品免费| 日韩欧美一区二区三区在线观看 | 久久久国产精品麻豆| 亚洲精品在线美女| 午夜成年电影在线免费观看|