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

    基于EMD與BP神經(jīng)網(wǎng)絡(luò)的北京市采暖季天然氣消耗量預(yù)測

    2021-11-15 09:01:42閔晶晶
    干旱氣象 2021年5期
    關(guān)鍵詞:消耗量分量氣象

    閔晶晶,王 華,董 顏

    (北京市氣象服務(wù)中心,北京 100097)

    引 言

    近年來,我國經(jīng)濟的迅速發(fā)展導(dǎo)致能源消費量快速增加,能源安全與環(huán)境問題日益顯著。天然氣作為優(yōu)質(zhì)高效的清潔能源極大改善了空氣環(huán)境質(zhì)量。伴隨北京煤改氣的全面推進,天然氣消費量呈現(xiàn)爆發(fā)式增長[1-4]。合理預(yù)測天然氣消耗量,不僅有利于城市管網(wǎng)的優(yōu)化調(diào)度,保證管網(wǎng)供氣量,而且在實現(xiàn)節(jié)能減排、國家政策分析與制定以及國民經(jīng)濟核算等方面具有深遠意義。據(jù)統(tǒng)計,天然氣作為北京市主要的供暖燃料,采暖季的消耗量占全年消耗總量的80%,具有明顯的季節(jié)性特征。如2009—2010、2012—2013年冬季北京遭遇了罕見的大雪、低溫天氣,天然氣出現(xiàn)供需不平衡情況,直接威脅到城市居民的冬季供暖。因此,災(zāi)害性天氣影響下冬季采暖季天然氣預(yù)測顯得尤為重要。

    天然氣消耗量常用的預(yù)測方法主要分為統(tǒng)計方法、經(jīng)濟數(shù)學模型和智能算法三類[5]。統(tǒng)計學方法以時間序列分析為主,廣泛應(yīng)用于能源預(yù)測[6-8],但統(tǒng)計模型需要長期數(shù)據(jù)積累,且模型中影響因子的貢獻率較為固定,無法反映隨時間推移某些特征因子貢獻率的衰減。經(jīng)濟數(shù)學模型如灰色系統(tǒng)理論,可量化不確定因素和小樣本信息,常用于構(gòu)建短期天然氣消耗量預(yù)測模型[9-11]。近年來,為了解決非線性決策問題,陸續(xù)采用神經(jīng)網(wǎng)絡(luò)、遺傳算法等智能算法構(gòu)建天然氣消耗量預(yù)測模型[12-15]。經(jīng)驗?zāi)B(tài)分解是處理非平穩(wěn)信號的一種時頻分析方法,可依據(jù)輸入信號的自身特點,將原始信號分解成若干個本征模態(tài)函數(shù)與殘差,該方法是對傅里葉分析和小波變換的重大突破[16]。BP神經(jīng)網(wǎng)絡(luò)方法在處理非線性問題上具有極大的靈活性和自適應(yīng)性,可通過訓(xùn)練數(shù)據(jù)集合不斷刺激系統(tǒng)來調(diào)整網(wǎng)絡(luò)權(quán)值和閾值,最終達到網(wǎng)絡(luò)誤差平方和最小的標準[17]。

    從天然氣消耗量的影響因素來看,一方面受國民經(jīng)濟、城市人口、能源政策等社會經(jīng)濟發(fā)展要素的影響,燃氣負荷表現(xiàn)出長期變化趨勢[18];另一方面受天氣影響,燃氣負荷呈現(xiàn)季節(jié)變化和波動特征[19-21]。然而,以往研究天然氣消耗量的影響因素大多局限在溫度、濕度、天氣狀況、降雨量、風等基本要素,且多以溫度作為唯一氣象參考因子[22-23],沒有過多考慮持續(xù)低溫或強降雪等災(zāi)害性天氣過程造成天然氣消耗量大幅增加的情況,同時也缺乏社會經(jīng)濟發(fā)展要素影響的考量。為此,本文以北京市為例,采用經(jīng)驗?zāi)B(tài)分解、BP神經(jīng)網(wǎng)絡(luò)和相關(guān)分析等方法,根據(jù)天然氣消耗量的非線性和非平穩(wěn)的特點,綜合分析經(jīng)濟、人口、供熱面積和氣象等多種影響要素,構(gòu)建北京市采暖季天然氣消耗量預(yù)測模型,并進行評估檢驗,以期為北京能源規(guī)劃和調(diào)控提供一定參考。

    1 資料和方法

    1.1 資 料

    天然氣資料為北京市2002—2018年采暖季(11月至次年3月)天然氣月消耗量,累加后得到采暖季天然氣總消耗量,其中2002—2015年作為預(yù)測模型構(gòu)建樣本,2016—2018年作為模型評估樣本。

    北京市有20個國家級氣象觀測站,其中佛爺頂、湯河口、上甸子、齋堂、霞云嶺5站海拔較高且地處遠郊,人類活動較少,不具代表性,因此選用采暖季其他15站最高、最低氣溫及平均氣溫和降水量、風速等氣象要素逐日觀測資料,并用上述15站平均值代表北京市氣候特征,同時還引入了表征冬季冷暖程度的負積溫和低溫日數(shù)2個參數(shù)[24]。其中,負積溫是指日平均氣溫低于0 ℃的累積值;低溫日數(shù)是指日最低氣溫小于某個閾值的累積天數(shù),本文將北京供暖部門關(guān)注的-5 ℃和-8 ℃作為閾值。

    社會因素資料選用了2002—2018年GDP、集中供熱面積、常住人口。

    1.2 經(jīng)驗?zāi)B(tài)分解方法

    經(jīng)驗?zāi)B(tài)分解(empirical mode decomposition, EMD)作為數(shù)據(jù)驅(qū)動的自適應(yīng)非線性時變信號分解方法,是以傅里葉分析和小波分析為基礎(chǔ),適用于對非線性、非平穩(wěn)信號逐級進行平穩(wěn)化處理。EMD算法的篩選過程是將復(fù)雜時序數(shù)據(jù)分解為有限個本征模函數(shù)(intrinsic mode function, IMF),分解得到的IMF分量包含了原數(shù)據(jù)在不同時間尺度上的波動信息[25]。具體處理過程如下:

    (1)對任一處理信號x(t),確定其局部最大值、最小值,將信號數(shù)據(jù)x(t)與上、下包絡(luò)線均值m1(t)的差值記為h1(t),其表達式如下:

    h1(t)=x(t)-m1(t)

    (1)

    (2)重復(fù)以上過程,從原始信號中篩選出的第一階段IMF分量h1(t)通常包含信號的最高頻分量。將h1(t)從x(t)中分離出來,得到去除高頻分量的差值信號r1(t),把r1(t)作為新信號重復(fù)以上篩分步驟,直到第n階段的殘余信號為單調(diào)函數(shù)不再篩分出IMF分量的rn(t),其表達式如下:

    (2)

    根據(jù)分解算法,x(t)可表示為n個IMF分量和1個殘差的和,表達式如下:

    (3)

    式中:rn(t)為殘差,表示信號中的平均趨勢;hj(t)為第j個IMF分量,j=1,2,…,n,分別表示信號從高到低不同頻率的成分。

    篩分過程的結(jié)束主要依據(jù)類似Cauchy收斂準則[26],標準偏差SD通常設(shè)置為0.2~0.3,表達式如下:

    (4)

    1.3 BP神經(jīng)網(wǎng)絡(luò)方法

    采用BP(back propagation)神經(jīng)網(wǎng)絡(luò)方法構(gòu)建預(yù)測模型。BP神經(jīng)網(wǎng)絡(luò)適用于一些高度復(fù)雜和非線性問題的建模,主要包含輸入層、輸出層和隱含層。一般采用梯度下降法作為權(quán)值修正,常利用Sigmoid函數(shù)計算各節(jié)點神經(jīng)元的激勵函數(shù),并歸一化到-1~1之間。設(shè)W(k+1)和W(k)為BP網(wǎng)絡(luò)各層節(jié)點之間相鄰訓(xùn)練數(shù)的連接權(quán)值[27],△W(k)為權(quán)值的修正值,k為訓(xùn)練次數(shù),權(quán)值修正公式如下:

    W(k+1)=W(k)+△W(k)

    (5)

    設(shè)有N組訓(xùn)練樣本和M組檢測樣本,Ti和Yi(k)分別為第i組訓(xùn)練樣本的目標值和網(wǎng)絡(luò)實際輸出值,而Tj和Yj(k)分別為第j組檢驗樣本的目標值和網(wǎng)絡(luò)實際輸出值,則訓(xùn)練誤差函數(shù)Etrain(k)和檢驗誤差函數(shù)Etest(k)可表示為:

    (6)

    (7)

    2 天然氣消耗量特征及其影響要素

    2.1 采暖季天然氣消耗量年際變化特征

    2002—2018年,北京地區(qū)采暖季天然氣消耗量呈波動增加趨勢,2008年以前增加緩慢,之后迅速增加[圖1(a)];采暖季天然氣消耗量的IMF1高頻分量呈波動變化,無明顯變化趨勢,2008年以前變化較平穩(wěn),之后波動幅度明顯增大,最高峰值和次高峰值分別出現(xiàn)在2012年和2009年,谷值出現(xiàn)在2013年和2016年[圖1(b)]。與采暖季天然氣消耗量的變化趨勢相似,天然氣消耗量的IMF2趨勢項呈持續(xù)增加趨勢,2002—2012年消耗量平穩(wěn)增加,2014—2015年急劇增加[圖1(c)],增量超過20%。

    圖1 北京地區(qū)2002—2018年采暖季天然氣消耗量(a)及分解量IMF1(b)和IMF2(c)的年際變化Fig.1 The inter-annual changes of natural gas consumption (a) and its decomposed IMF1 (b) and IMF2 (c) components in Beijing in heating season from 2002 to 2018

    2.2 影響要素及年際變化特征

    天然氣消耗量受社會因素和氣象條件的影響,具有典型的非平穩(wěn)特征,其中氣象條件是造成天然氣負荷短期波動的主要原因。利用EMD方法分析北京市2002—2018年采暖季天然氣消耗量,得到不同頻率的本征模態(tài)分量IMF1和趨勢分量IMF2,并與各影響因子作相關(guān)分析(表1)。可以看出,IMF1與平均氣溫、最低氣溫、最高氣溫、負積溫均呈顯著負相關(guān),而與持續(xù)低溫日數(shù)和降水量呈顯著正相關(guān);IMF2與GDP、集中供熱面積、常住人口數(shù)量呈顯著正相關(guān),而與氣象因子相關(guān)性不大,表明IMF2分量綜合反映了GDP、人口、城市化、能源結(jié)構(gòu)等社會因素的影響。因此,設(shè)定高頻分量IMF1為氣象消耗量,趨勢分量IMF2為社會消耗量。

    表1 北京市2002—2018年采暖季天然氣消耗量的IMF1和IMF2分量與各因子的相關(guān)系數(shù)Tab.1 Correlation coefficients between IMF1, IMF2 components of natural gas consumption and each factor in Beijing in heating season from 2002 to 2018

    2.2.1 氣象消耗量年際變化特征

    近17 a北京地區(qū)采暖季平均氣溫大多在0 ℃以上,2009年和2012年采暖季因出現(xiàn)多次強降雪過程導(dǎo)致平均氣溫較常年同期(0.7 ℃)明顯偏低,分別為-1.3、-0.9 ℃[圖2(a)],而降水量較常年同期(27.7 mm)異常偏多,分別為74.6、106 mm[圖2(b)],致使低溫持續(xù)時間較長,低溫日數(shù)多于常年,負積溫低于-400 ℃(常年平均為-295.6 ℃);2013、2014、2016年采暖季平均氣溫分別為2.4、2.3、2.0 ℃,均明顯高于常年同期,而降水量相對較少,低于-8 ℃的低溫日數(shù)在20 d左右,較常年同期(42 d)明顯偏少,負積溫僅在-150 ℃左右。因此,當采暖季氣溫明顯偏低或出現(xiàn)較強降雪、持續(xù)低溫等天氣過程時,天然氣的氣象消耗量會大幅增加,而當采暖季氣溫明顯偏高、負積溫偏少時,天然氣的氣象消耗量明顯減少。

    圖2 北京地區(qū)2002—2018年采暖季平均氣溫(a)、降水量(b)年際變化Fig.2 The inter-annual changes of average temperature (a) and precipitation (b) in Beijing in heating season from 2002 to 2018

    2.2.2 社會消耗量年際變化特征

    近17 a北京市GDP、集中供熱面積持續(xù)平穩(wěn)增長,而常住人口數(shù)量增長相對緩慢,特別是2015年以后人口略有減少(圖3),其中2013—2015年GDP、集中供熱面積及人口數(shù)量均線性平穩(wěn)增長,未出現(xiàn)與天然氣社會消耗量IMF2相對應(yīng)的明顯增長現(xiàn)象。2013年9月北京市政府發(fā)布了《北京市2013—2017年清潔空氣行動計劃》,這可能是導(dǎo)致2013—2015年采暖季天然氣用量明顯增加的主要原因。

    圖3 北京地區(qū)2002—2018年采暖季GDP、集中供熱面積及常住人口數(shù)量的年際變化Fig.3 The inter-annual changes of GDP, intensive heating supply area and permanent population number in Beijing in heating season from 2002 to 2018

    3 天然氣消耗量預(yù)測模型與檢驗

    3.1 預(yù)測模型的建立

    基于EMD方法在處理非平穩(wěn)信號和BP神經(jīng)網(wǎng)絡(luò)方法在處理非線性問題的優(yōu)勢,提出了基于EMD_BP方法的采暖季天然氣消耗量預(yù)測模型(圖4)。選取平均氣溫Tave、最低氣溫Tmin、最高氣溫Tmax、風速V、降水量R、負積溫NAT、小于-8 ℃的低溫日數(shù)LTD-8等作為IMF1預(yù)測模型的輸入因子,GDP、集中供熱面積、常住人口數(shù)量、時間(年份)等作為IMF2預(yù)測模型的輸入因子,分別建立天然氣氣象消耗量IMF1和社會消耗量IMF2的預(yù)測模型,兩分量預(yù)測值求和得到天然氣總消耗量預(yù)測值。其中,2002—2015年采暖季天然氣消耗量數(shù)據(jù)作為訓(xùn)練樣本,2016—2018年數(shù)據(jù)作為檢驗樣本。

    3.2 預(yù)測結(jié)果檢驗

    利用采暖季天然氣消耗量EMD_BP預(yù)測模型,分別對2002—2015年建模樣本進行擬合、2016—2018年檢驗樣本進行預(yù)測??梢钥闯?,EMD_BP預(yù)測模型對天然氣消耗量IMF1、IMF2分量和總消耗量均具有很好的擬合與預(yù)測效果(圖5),尤其是氣象消耗量IMF1預(yù)測模型能夠準確預(yù)測出陡增點(2009、2012、2015、2017年)和陡減點(2010、2013、2016年)。

    圖5 北京地區(qū)2002—2018年采暖季天然氣消耗量IMF1(a)、IMF2(b)分量和總消耗量(c)的擬合值、預(yù)測值與真實值對比Fig.5 Comparison of fitted and forecasted IMF1 (a), IMF2 (b) components of natural gas consumption and total consumption (c) with actual value in Beijing in heating season from 2002 to 2018

    通過擬合、預(yù)測與真實值的結(jié)果進行檢驗,發(fā)現(xiàn)EMD_BP預(yù)測模型對2002—2015年數(shù)據(jù)樣本擬合的平均相對誤差為4.5%,對2016—2018年檢驗樣本預(yù)測的平均相對誤差為5.6%??傮w來看,采暖季天然氣消耗量EMD_BP預(yù)測模型具有較好的預(yù)測效果,能夠準確預(yù)測出天氣變化帶來的影響。

    4 結(jié)論和討論

    本文引入對非線性非平穩(wěn)信號進行平穩(wěn)化處理的EMD方法,對北京市采暖季天然氣消耗量的時間序列進行分解及特征分析,探討天然氣消耗量的社會和氣象影響因子,采用BP神經(jīng)網(wǎng)絡(luò)方法構(gòu)建預(yù)測模型并檢驗,主要結(jié)論如下:

    (1)EMD方法能夠較好地提取北京市采暖季天然氣消耗量線性增加的社會消耗量和波動變化的氣象消耗量。

    (2)北京市采暖季天然氣的社會消耗量與GDP、集中供熱面積、常住人口數(shù)量呈顯著正相關(guān),但受政策影響,2014—2015年天然氣的社會消耗量出現(xiàn)劇增,此后增量有所減緩。

    (3)北京市采暖季天然氣的氣象消耗量與氣溫、負積溫呈顯著負相關(guān),而與降水量、小于-8 ℃的低溫日數(shù)呈顯著正相關(guān)。當采暖季氣溫明顯偏低或出現(xiàn)較強降雪、持續(xù)低溫等天氣過程時,天然氣的氣象消耗量會大幅增加,而在氣溫明顯偏高、負積溫偏少時,氣象消耗量明顯減少。

    (4)采暖季天然氣消耗量EMD_BP預(yù)測模型具有較好的預(yù)測效果,2002—2015年數(shù)據(jù)樣本擬合的平均相對誤差為4.5%,2016—2018年檢驗樣本預(yù)測的平均相對誤差為5.6%,能夠準確預(yù)測因天氣變化而引起的天然氣氣象消耗量的峰谷變化。

    猜你喜歡
    消耗量分量氣象
    氣象
    氣象樹
    路基石方爆破降低炸藥消耗量研究
    《內(nèi)蒙古氣象》征稿簡則
    帽子的分量
    一物千斤
    智族GQ(2019年9期)2019-10-28 08:16:21
    論《哈姆雷特》中良心的分量
    大國氣象
    分量
    有機化學反應(yīng)中試劑最大消耗量問題探析
    老司机福利观看| 国产精品99久久久久久久久| 精品人妻偷拍中文字幕| 国产av不卡久久| videossex国产| 麻豆精品久久久久久蜜桃| 成人综合一区亚洲| 欧美日韩综合久久久久久 | 亚洲最大成人手机在线| 级片在线观看| 欧美日韩综合久久久久久 | 桃色一区二区三区在线观看| 亚洲精华国产精华液的使用体验 | 国产高清不卡午夜福利| 国产精品一及| 亚洲精华国产精华精| 欧美成人免费av一区二区三区| 中文资源天堂在线| 国产真实伦视频高清在线观看 | 少妇高潮的动态图| 国产高清激情床上av| 国产高清有码在线观看视频| 欧美丝袜亚洲另类 | 亚洲av熟女| 97超级碰碰碰精品色视频在线观看| 美女被艹到高潮喷水动态| 国产高清不卡午夜福利| 深夜精品福利| 99久久精品国产国产毛片| 国产在视频线在精品| 亚洲av二区三区四区| 一级a爱片免费观看的视频| 久久久精品欧美日韩精品| 午夜精品久久久久久毛片777| 国产久久久一区二区三区| 亚洲欧美日韩东京热| 国产毛片a区久久久久| 久久热精品热| 亚洲成人精品中文字幕电影| 国产一区二区三区在线臀色熟女| x7x7x7水蜜桃| 久久精品国产亚洲av涩爱 | 成人毛片a级毛片在线播放| 亚洲av美国av| 精品一区二区三区视频在线| 午夜视频国产福利| 亚洲自拍偷在线| 欧美性猛交╳xxx乱大交人| 国产成人av教育| 亚洲av中文字字幕乱码综合| 自拍偷自拍亚洲精品老妇| 久久久精品欧美日韩精品| 亚洲人成伊人成综合网2020| 伊人久久精品亚洲午夜| 午夜福利欧美成人| 欧美成人a在线观看| 国产久久久一区二区三区| 国产精品98久久久久久宅男小说| 精品乱码久久久久久99久播| 久久久久九九精品影院| 中文字幕久久专区| 国产精品久久久久久亚洲av鲁大| 欧美性猛交黑人性爽| 国产美女午夜福利| 国产欧美日韩精品一区二区| 国产高潮美女av| 长腿黑丝高跟| 99热这里只有是精品在线观看| 国产精品不卡视频一区二区| 色在线成人网| 一个人观看的视频www高清免费观看| 嫩草影院精品99| x7x7x7水蜜桃| 伦理电影大哥的女人| www日本黄色视频网| 一个人免费在线观看电影| 女人被狂操c到高潮| 亚洲最大成人中文| 久久精品国产清高在天天线| 在线免费十八禁| 国产探花在线观看一区二区| 亚洲精品一区av在线观看| 免费大片18禁| 黄色丝袜av网址大全| 国语自产精品视频在线第100页| 成人国产麻豆网| 夜夜爽天天搞| 免费观看精品视频网站| 久久人人精品亚洲av| 我要搜黄色片| 啦啦啦啦在线视频资源| 欧美绝顶高潮抽搐喷水| 国产精品久久久久久精品电影| 国模一区二区三区四区视频| a级一级毛片免费在线观看| 不卡一级毛片| 亚洲狠狠婷婷综合久久图片| 内地一区二区视频在线| 国内精品美女久久久久久| 久久人人爽人人爽人人片va| 美女高潮的动态| 一个人看视频在线观看www免费| 搡老岳熟女国产| 亚洲欧美日韩卡通动漫| 国产色婷婷99| 最好的美女福利视频网| 一区二区三区免费毛片| 赤兔流量卡办理| 2021天堂中文幕一二区在线观| 欧美日韩乱码在线| 日韩欧美国产在线观看| 欧美一区二区精品小视频在线| 欧美性猛交黑人性爽| 国产亚洲欧美98| 特级一级黄色大片| 国产精品国产三级国产av玫瑰| 欧美色视频一区免费| 少妇的逼水好多| av在线观看视频网站免费| 搡老妇女老女人老熟妇| 国产色爽女视频免费观看| 悠悠久久av| 亚洲无线在线观看| 蜜桃久久精品国产亚洲av| 在线a可以看的网站| 久久久久久久久久黄片| 性欧美人与动物交配| 精品一区二区三区人妻视频| 精品久久国产蜜桃| 人妻久久中文字幕网| av天堂中文字幕网| 国产成人福利小说| 中出人妻视频一区二区| 色哟哟哟哟哟哟| 欧美日韩中文字幕国产精品一区二区三区| 成人二区视频| 免费观看在线日韩| 久久精品综合一区二区三区| 精品久久久久久成人av| 国产爱豆传媒在线观看| 夜夜爽天天搞| 中文字幕久久专区| 在线观看一区二区三区| 久9热在线精品视频| 在线国产一区二区在线| 热99在线观看视频| 热99re8久久精品国产| 老司机深夜福利视频在线观看| 久久精品久久久久久噜噜老黄 | 国产乱人伦免费视频| 一进一出好大好爽视频| 久久久久性生活片| 国产精品亚洲一级av第二区| 亚洲成人中文字幕在线播放| 亚洲精品一区av在线观看| av天堂在线播放| 此物有八面人人有两片| 99久久精品热视频| 久久6这里有精品| 国产成人aa在线观看| 欧美成人免费av一区二区三区| 日本三级黄在线观看| 成人精品一区二区免费| 网址你懂的国产日韩在线| 国产免费一级a男人的天堂| 国产真实伦视频高清在线观看 | 午夜福利成人在线免费观看| 亚洲在线观看片| h日本视频在线播放| 亚洲欧美日韩卡通动漫| 99riav亚洲国产免费| x7x7x7水蜜桃| 别揉我奶头 嗯啊视频| 日韩一区二区视频免费看| 国内精品一区二区在线观看| 国产极品精品免费视频能看的| 亚洲乱码一区二区免费版| 国产成年人精品一区二区| 国产精品久久久久久av不卡| bbb黄色大片| 精品久久久噜噜| 日韩国内少妇激情av| 免费在线观看影片大全网站| 精品久久久久久久久亚洲 | av在线老鸭窝| 一个人看视频在线观看www免费| 999久久久精品免费观看国产| 深夜a级毛片| 国产男人的电影天堂91| а√天堂www在线а√下载| 狂野欧美白嫩少妇大欣赏| 国内精品久久久久精免费| 亚洲成人免费电影在线观看| 日韩中字成人| 亚洲性夜色夜夜综合| 精品无人区乱码1区二区| 日韩欧美精品免费久久| 国产真实乱freesex| 国产在线精品亚洲第一网站| 在线a可以看的网站| 免费av毛片视频| 色精品久久人妻99蜜桃| 亚洲av电影不卡..在线观看| 黄色配什么色好看| 97人妻精品一区二区三区麻豆| 天美传媒精品一区二区| 久久午夜亚洲精品久久| 99久久成人亚洲精品观看| 亚洲va日本ⅴa欧美va伊人久久| 我的女老师完整版在线观看| 搡老岳熟女国产| 精品欧美国产一区二区三| 国产伦一二天堂av在线观看| 久久精品影院6| 色av中文字幕| 日本与韩国留学比较| 高清毛片免费观看视频网站| 日日摸夜夜添夜夜添av毛片 | 可以在线观看的亚洲视频| 高清在线国产一区| 国产亚洲欧美98| 老师上课跳d突然被开到最大视频| 成人特级黄色片久久久久久久| 色吧在线观看| 乱人视频在线观看| 亚洲性夜色夜夜综合| 亚洲欧美激情综合另类| 亚洲av五月六月丁香网| 国产91精品成人一区二区三区| 麻豆国产av国片精品| 在线国产一区二区在线| 22中文网久久字幕| 国产私拍福利视频在线观看| 亚洲精品影视一区二区三区av| 亚洲图色成人| 国产高清不卡午夜福利| 国产精品久久久久久精品电影| 最后的刺客免费高清国语| 中文字幕人妻熟人妻熟丝袜美| 女同久久另类99精品国产91| 亚洲av美国av| 长腿黑丝高跟| 日韩高清综合在线| 亚洲人成网站高清观看| 国产女主播在线喷水免费视频网站 | 久久精品国产99精品国产亚洲性色| 国产精品乱码一区二三区的特点| 日韩精品有码人妻一区| 日本一二三区视频观看| 亚洲中文日韩欧美视频| 人妻夜夜爽99麻豆av| 变态另类丝袜制服| xxxwww97欧美| 欧美绝顶高潮抽搐喷水| av国产免费在线观看| 少妇猛男粗大的猛烈进出视频 | 国产三级在线视频| 最近中文字幕高清免费大全6 | 国语自产精品视频在线第100页| 他把我摸到了高潮在线观看| 国产亚洲精品久久久com| 天堂av国产一区二区熟女人妻| 国产精品久久久久久亚洲av鲁大| 国产黄色小视频在线观看| 亚洲国产精品久久男人天堂| 亚洲专区国产一区二区| 免费高清视频大片| 一进一出抽搐动态| 欧美日韩国产亚洲二区| 亚洲自偷自拍三级| 国内久久婷婷六月综合欲色啪| 久久欧美精品欧美久久欧美| 成人av在线播放网站| 国产精品亚洲一级av第二区| 男女视频在线观看网站免费| 国产探花极品一区二区| 最近视频中文字幕2019在线8| 精品久久久噜噜| 久久这里只有精品中国| 国产精品久久久久久久久免| 国产乱人视频| 乱系列少妇在线播放| 国产亚洲精品久久久久久毛片| 国产黄色小视频在线观看| 看片在线看免费视频| 18禁黄网站禁片免费观看直播| 国产精品美女特级片免费视频播放器| 国内精品美女久久久久久| 久久天躁狠狠躁夜夜2o2o| 国产精品人妻久久久久久| 亚洲自偷自拍三级| 久久久久久久亚洲中文字幕| 色综合亚洲欧美另类图片| 久久精品国产99精品国产亚洲性色| 成人毛片a级毛片在线播放| 亚洲va在线va天堂va国产| 国产精品乱码一区二三区的特点| 欧美性感艳星| 色5月婷婷丁香| 国产精品一区二区三区四区久久| 国产av一区在线观看免费| 精品久久久久久久久亚洲 | 色噜噜av男人的天堂激情| 亚洲成a人片在线一区二区| 1000部很黄的大片| 国产淫片久久久久久久久| 变态另类成人亚洲欧美熟女| 一区二区三区免费毛片| 又紧又爽又黄一区二区| 天天躁日日操中文字幕| 国产毛片a区久久久久| 超碰av人人做人人爽久久| 麻豆一二三区av精品| 精品午夜福利在线看| 九九爱精品视频在线观看| 精品久久久久久,| 午夜激情福利司机影院| 久久久午夜欧美精品| 久久久色成人| 国内精品美女久久久久久| 俺也久久电影网| 99久久九九国产精品国产免费| 午夜亚洲福利在线播放| 在线播放无遮挡| 国国产精品蜜臀av免费| 中文字幕久久专区| 国产成人福利小说| 啦啦啦韩国在线观看视频| 我的女老师完整版在线观看| 蜜桃亚洲精品一区二区三区| 无遮挡黄片免费观看| 国产视频内射| 久久人人精品亚洲av| 国产精品乱码一区二三区的特点| 99热只有精品国产| 精品一区二区三区av网在线观看| 十八禁网站免费在线| 亚洲天堂国产精品一区在线| 久久精品夜夜夜夜夜久久蜜豆| 啦啦啦啦在线视频资源| 国产 一区精品| 夜夜夜夜夜久久久久| 亚洲 国产 在线| 成人性生交大片免费视频hd| 精品人妻偷拍中文字幕| bbb黄色大片| 国内精品美女久久久久久| 久久久久国内视频| h日本视频在线播放| 一区二区三区四区激情视频 | 啦啦啦韩国在线观看视频| 天堂动漫精品| 国内精品久久久久精免费| 88av欧美| 十八禁国产超污无遮挡网站| 国产成年人精品一区二区| 午夜久久久久精精品| 婷婷六月久久综合丁香| 老女人水多毛片| 色综合亚洲欧美另类图片| 亚洲av不卡在线观看| 3wmmmm亚洲av在线观看| 99热这里只有是精品在线观看| 熟女电影av网| 18禁在线播放成人免费| 成熟少妇高潮喷水视频| 国产亚洲精品av在线| 久久中文看片网| 精品无人区乱码1区二区| 看十八女毛片水多多多| 深夜精品福利| 蜜桃亚洲精品一区二区三区| 免费观看在线日韩| 校园人妻丝袜中文字幕| 成人综合一区亚洲| 午夜a级毛片| 精品无人区乱码1区二区| 一级黄色大片毛片| 久久精品综合一区二区三区| 直男gayav资源| 国产男靠女视频免费网站| 午夜福利视频1000在线观看| 一个人看视频在线观看www免费| 亚洲内射少妇av| 欧美高清性xxxxhd video| 日韩一本色道免费dvd| 中亚洲国语对白在线视频| 国产伦一二天堂av在线观看| 亚洲国产欧洲综合997久久,| 人妻久久中文字幕网| 18禁黄网站禁片免费观看直播| 女的被弄到高潮叫床怎么办 | 直男gayav资源| 久久精品国产清高在天天线| 日韩大尺度精品在线看网址| 国产一区二区亚洲精品在线观看| 中国美白少妇内射xxxbb| 一级毛片久久久久久久久女| 黄色一级大片看看| 18禁裸乳无遮挡免费网站照片| 欧美三级亚洲精品| 国产成人aa在线观看| 人人妻,人人澡人人爽秒播| 日本撒尿小便嘘嘘汇集6| 国产精品福利在线免费观看| 老司机福利观看| 久久亚洲精品不卡| 日本一本二区三区精品| 国产一区二区在线av高清观看| 久久久久久久久中文| 99九九线精品视频在线观看视频| 欧美精品啪啪一区二区三区| 久久久久国内视频| 亚洲成人久久性| 18禁黄网站禁片午夜丰满| 亚洲精品在线观看二区| 午夜福利成人在线免费观看| 免费电影在线观看免费观看| 国产伦精品一区二区三区四那| a级毛片a级免费在线| 韩国av一区二区三区四区| 国内毛片毛片毛片毛片毛片| 精品福利观看| 桃红色精品国产亚洲av| 亚洲专区中文字幕在线| 性欧美人与动物交配| 久久精品久久久久久噜噜老黄 | 久久久精品欧美日韩精品| 特级一级黄色大片| videossex国产| 永久网站在线| 51国产日韩欧美| 欧美日本亚洲视频在线播放| 人妻久久中文字幕网| 久久99热这里只有精品18| 精品一区二区三区av网在线观看| 日韩强制内射视频| 国产毛片a区久久久久| 在线观看午夜福利视频| 我的老师免费观看完整版| 精品无人区乱码1区二区| 久久国产乱子免费精品| 日韩欧美国产在线观看| 国产一区二区亚洲精品在线观看| 内射极品少妇av片p| 一级a爱片免费观看的视频| 成年女人毛片免费观看观看9| 久久亚洲精品不卡| 国内精品久久久久精免费| 成人特级av手机在线观看| 欧美+日韩+精品| 一个人观看的视频www高清免费观看| 九九在线视频观看精品| 国产美女午夜福利| 精品久久久久久久久亚洲 | 亚洲最大成人av| 久久精品国产亚洲av涩爱 | 制服丝袜大香蕉在线| 赤兔流量卡办理| 国产国拍精品亚洲av在线观看| 亚洲中文日韩欧美视频| 嫩草影院新地址| 尾随美女入室| 乱人视频在线观看| 欧美不卡视频在线免费观看| 日本一二三区视频观看| 婷婷精品国产亚洲av在线| 欧美3d第一页| 热99在线观看视频| 成人一区二区视频在线观看| 国内毛片毛片毛片毛片毛片| 日本成人三级电影网站| 国产精品一区www在线观看 | 国产 一区 欧美 日韩| 成人国产麻豆网| 精品99又大又爽又粗少妇毛片 | 久久久久性生活片| 亚洲电影在线观看av| 日韩av在线大香蕉| 成人永久免费在线观看视频| 亚洲最大成人av| 午夜日韩欧美国产| 欧美一区二区国产精品久久精品| 国内精品久久久久精免费| 99热只有精品国产| 国产 一区 欧美 日韩| 久久久久久久久大av| 亚洲美女搞黄在线观看 | 亚洲最大成人中文| 18禁在线播放成人免费| 91在线观看av| 丰满乱子伦码专区| 欧美最黄视频在线播放免费| 亚州av有码| 国产色婷婷99| 搡老妇女老女人老熟妇| 免费av不卡在线播放| 久久午夜亚洲精品久久| 免费看a级黄色片| 久久久久久久久久久丰满 | 国内久久婷婷六月综合欲色啪| 亚洲国产欧洲综合997久久,| 亚洲,欧美,日韩| videossex国产| 国产高清视频在线播放一区| 亚洲精品456在线播放app | 麻豆一二三区av精品| 日本爱情动作片www.在线观看 | АⅤ资源中文在线天堂| 国产精品乱码一区二三区的特点| 男人和女人高潮做爰伦理| 一个人看的www免费观看视频| 一本一本综合久久| 亚洲精品日韩av片在线观看| 日本一本二区三区精品| 少妇裸体淫交视频免费看高清| 变态另类成人亚洲欧美熟女| 非洲黑人性xxxx精品又粗又长| 搡女人真爽免费视频火全软件 | 又黄又爽又免费观看的视频| 亚洲最大成人手机在线| 国产 一区精品| 亚洲在线自拍视频| 日本爱情动作片www.在线观看 | 国产黄a三级三级三级人| 不卡视频在线观看欧美| 久久午夜亚洲精品久久| 搡老岳熟女国产| 亚洲欧美日韩高清在线视频| 日韩欧美免费精品| 亚洲电影在线观看av| 亚洲一区高清亚洲精品| 国产精品国产三级国产av玫瑰| 日日摸夜夜添夜夜添av毛片 | 午夜福利在线在线| 天堂影院成人在线观看| 18禁在线播放成人免费| 免费观看在线日韩| 国产单亲对白刺激| 免费高清视频大片| 国产毛片a区久久久久| 别揉我奶头~嗯~啊~动态视频| eeuss影院久久| 黄色欧美视频在线观看| 亚洲av熟女| 午夜福利高清视频| 亚洲男人的天堂狠狠| 欧美一区二区国产精品久久精品| 波野结衣二区三区在线| 日本撒尿小便嘘嘘汇集6| 免费电影在线观看免费观看| 国产免费男女视频| 久9热在线精品视频| 美女免费视频网站| 亚洲三级黄色毛片| 久久亚洲真实| 欧美日韩黄片免| 韩国av在线不卡| 99久久中文字幕三级久久日本| 国产精品亚洲美女久久久| 日本欧美国产在线视频| 精品久久久久久,| 老师上课跳d突然被开到最大视频| 深爱激情五月婷婷| 日本-黄色视频高清免费观看| 亚洲精品456在线播放app | 内射极品少妇av片p| 欧美色视频一区免费| 午夜a级毛片| 国产精品人妻久久久久久| 高清毛片免费观看视频网站| 成人三级黄色视频| 亚洲三级黄色毛片| 国产精品久久久久久av不卡| 99国产极品粉嫩在线观看| 真实男女啪啪啪动态图| 欧美成人性av电影在线观看| 中文字幕av成人在线电影| 国产久久久一区二区三区| 舔av片在线| 不卡一级毛片| 国产精品久久电影中文字幕| 国产一区二区三区在线臀色熟女| 免费观看在线日韩| 午夜福利成人在线免费观看| 综合色av麻豆| 国产熟女欧美一区二区| 极品教师在线视频| 91久久精品国产一区二区三区| 在线a可以看的网站| 欧美3d第一页| 国产黄a三级三级三级人| 在线观看66精品国产| 欧美3d第一页| 麻豆国产97在线/欧美| 国产精品98久久久久久宅男小说| 亚洲成a人片在线一区二区| 亚洲av成人av| 国产黄片美女视频| 日本黄色片子视频| 精品久久久久久久久av| 女同久久另类99精品国产91| 成年女人毛片免费观看观看9| 国产一区二区亚洲精品在线观看| 国产精品亚洲一级av第二区| 91麻豆av在线| 亚洲精品国产成人久久av| 午夜福利在线观看免费完整高清在 | 偷拍熟女少妇极品色| 国产一级毛片七仙女欲春2| 九九爱精品视频在线观看| 18禁在线播放成人免费| 一区二区三区高清视频在线| 精品人妻一区二区三区麻豆 | 国内精品宾馆在线|