• <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)中試劑最大消耗量問題探析
    黄色视频,在线免费观看| 在线天堂最新版资源| 欧美人与善性xxx| 内射极品少妇av片p| 日本成人三级电影网站| 2021天堂中文幕一二区在线观| 婷婷色综合大香蕉| av又黄又爽大尺度在线免费看 | 日韩欧美一区二区三区在线观看| 亚洲精品国产成人久久av| 少妇的逼好多水| 国产一区二区三区av在线 | 青春草亚洲视频在线观看| 亚洲精品456在线播放app| 美女黄网站色视频| 精品少妇黑人巨大在线播放 | 国产乱人偷精品视频| 99热只有精品国产| 亚洲美女视频黄频| 九九久久精品国产亚洲av麻豆| 久久精品91蜜桃| 国产一级毛片七仙女欲春2| 日韩成人伦理影院| 久久亚洲国产成人精品v| 又爽又黄无遮挡网站| 我的女老师完整版在线观看| 99久久久亚洲精品蜜臀av| 最近最新中文字幕大全电影3| 一进一出抽搐动态| 亚洲五月天丁香| 成年女人永久免费观看视频| 国内揄拍国产精品人妻在线| 你懂的网址亚洲精品在线观看 | 国产精品99久久久久久久久| 99视频精品全部免费 在线| 国产精品综合久久久久久久免费| 国产亚洲5aaaaa淫片| 欧美人与善性xxx| 99久久久亚洲精品蜜臀av| 精品少妇黑人巨大在线播放 | 一边亲一边摸免费视频| 欧美+日韩+精品| 成人鲁丝片一二三区免费| 一本一本综合久久| 国产精品久久视频播放| 日韩精品青青久久久久久| 久久久欧美国产精品| 婷婷色av中文字幕| 在线观看午夜福利视频| 97在线视频观看| 美女黄网站色视频| av又黄又爽大尺度在线免费看 | 综合色av麻豆| 麻豆国产av国片精品| 老司机福利观看| 国产精品电影一区二区三区| 亚洲成人av在线免费| 亚洲乱码一区二区免费版| 最新中文字幕久久久久| 一级毛片久久久久久久久女| 六月丁香七月| 国产v大片淫在线免费观看| 联通29元200g的流量卡| 欧美三级亚洲精品| 日日撸夜夜添| 欧美成人a在线观看| 午夜激情福利司机影院| 日韩一区二区视频免费看| 99视频精品全部免费 在线| 国产精品无大码| 床上黄色一级片| 日日干狠狠操夜夜爽| 欧美xxxx性猛交bbbb| 国产久久久一区二区三区| 成人综合一区亚洲| 波野结衣二区三区在线| 蜜臀久久99精品久久宅男| 久久婷婷人人爽人人干人人爱| 一边亲一边摸免费视频| 一本一本综合久久| 国产日韩欧美在线精品| 人体艺术视频欧美日本| 亚洲av.av天堂| 综合色丁香网| 欧美丝袜亚洲另类| 国产伦在线观看视频一区| 中文亚洲av片在线观看爽| 自拍偷自拍亚洲精品老妇| 国产成人精品婷婷| 日韩欧美精品免费久久| 欧美3d第一页| a级一级毛片免费在线观看| 国产高清不卡午夜福利| 国产一区二区三区av在线 | 人妻制服诱惑在线中文字幕| 只有这里有精品99| 一个人免费在线观看电影| 精品无人区乱码1区二区| 国产伦精品一区二区三区四那| 欧美三级亚洲精品| 亚洲自偷自拍三级| 亚洲国产色片| 国产成人精品一,二区 | 亚洲一区高清亚洲精品| 赤兔流量卡办理| 91久久精品电影网| 成年女人永久免费观看视频| 亚洲欧美中文字幕日韩二区| 青春草视频在线免费观看| 插逼视频在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲国产精品国产精品| 一个人看视频在线观看www免费| 自拍偷自拍亚洲精品老妇| 听说在线观看完整版免费高清| 国产美女午夜福利| 国产伦理片在线播放av一区 | 午夜精品国产一区二区电影 | 久久这里只有精品中国| 一个人观看的视频www高清免费观看| 欧美3d第一页| 免费看光身美女| 国产精品久久久久久亚洲av鲁大| 国产成年人精品一区二区| 欧美xxxx黑人xx丫x性爽| 在线播放国产精品三级| 成人亚洲欧美一区二区av| 国产一区二区在线观看日韩| 超碰av人人做人人爽久久| 简卡轻食公司| 久久久欧美国产精品| 国产亚洲5aaaaa淫片| 舔av片在线| 国产视频内射| 国产精品.久久久| 国产精品av视频在线免费观看| av在线播放精品| 晚上一个人看的免费电影| 日本熟妇午夜| 丰满人妻一区二区三区视频av| 久久久久九九精品影院| 久久久久久久久中文| 久久久午夜欧美精品| 狂野欧美激情性xxxx在线观看| 中文字幕精品亚洲无线码一区| 久久久欧美国产精品| 国产三级在线视频| 91久久精品国产一区二区成人| 成人亚洲精品av一区二区| 免费无遮挡裸体视频| 九九久久精品国产亚洲av麻豆| 在线观看午夜福利视频| 波多野结衣巨乳人妻| 免费看a级黄色片| 我要看日韩黄色一级片| 亚洲人成网站在线播| 伦理电影大哥的女人| 婷婷六月久久综合丁香| 亚洲av免费高清在线观看| 免费一级毛片在线播放高清视频| 日韩欧美精品v在线| 美女xxoo啪啪120秒动态图| 超碰av人人做人人爽久久| 黄色视频,在线免费观看| 97超视频在线观看视频| 精品一区二区三区人妻视频| 少妇熟女aⅴ在线视频| 免费一级毛片在线播放高清视频| 欧美+日韩+精品| 久久这里有精品视频免费| 久久亚洲国产成人精品v| 国产日韩欧美在线精品| 亚洲欧美精品综合久久99| 嘟嘟电影网在线观看| 99riav亚洲国产免费| 亚洲成人久久爱视频| 99热网站在线观看| av在线天堂中文字幕| 亚洲成人久久爱视频| 九九爱精品视频在线观看| 成人午夜精彩视频在线观看| 欧美日韩精品成人综合77777| 一级毛片电影观看 | 久久久久久国产a免费观看| 色播亚洲综合网| 国产精品一二三区在线看| 真实男女啪啪啪动态图| 国产黄片美女视频| 99久久久亚洲精品蜜臀av| 精华霜和精华液先用哪个| 久久久久久久久大av| 国产一区二区三区在线臀色熟女| 亚洲精品国产av成人精品| 一进一出抽搐动态| 国产 一区 欧美 日韩| 黄色配什么色好看| 免费不卡的大黄色大毛片视频在线观看 | 激情 狠狠 欧美| 国产麻豆成人av免费视频| 啦啦啦韩国在线观看视频| 国产熟女欧美一区二区| 国产一区二区激情短视频| 亚洲最大成人手机在线| 国产 一区 欧美 日韩| 天天躁夜夜躁狠狠久久av| 99热只有精品国产| 国产男人的电影天堂91| 性欧美人与动物交配| 天堂av国产一区二区熟女人妻| 网址你懂的国产日韩在线| 丰满人妻一区二区三区视频av| 国产成人aa在线观看| 亚洲性久久影院| 日韩欧美 国产精品| 国产精品久久久久久精品电影| 九九热线精品视视频播放| 夜夜夜夜夜久久久久| 久久人妻av系列| 亚洲最大成人中文| 波多野结衣高清无吗| 日韩强制内射视频| 热99re8久久精品国产| 又黄又爽又刺激的免费视频.| 日韩在线高清观看一区二区三区| av免费观看日本| 日韩强制内射视频| 熟女人妻精品中文字幕| 国产高清不卡午夜福利| 亚洲av免费在线观看| 麻豆国产av国片精品| 国产高清三级在线| 亚洲第一区二区三区不卡| 国产av不卡久久| 男插女下体视频免费在线播放| 久久久久久久久大av| 乱系列少妇在线播放| 精品国内亚洲2022精品成人| 中文字幕av成人在线电影| 久久久久久久久大av| 能在线免费看毛片的网站| 51国产日韩欧美| 99久久精品一区二区三区| 欧美区成人在线视频| 变态另类丝袜制服| 久久精品夜色国产| 中国美白少妇内射xxxbb| 亚洲va在线va天堂va国产| 亚洲精品456在线播放app| 国产免费男女视频| 草草在线视频免费看| 男人和女人高潮做爰伦理| 老师上课跳d突然被开到最大视频| 在线观看66精品国产| 热99在线观看视频| 日韩欧美精品免费久久| 国内精品久久久久精免费| 亚洲精品乱码久久久久久按摩| 亚洲国产欧洲综合997久久,| 51国产日韩欧美| 国内精品久久久久精免费| 午夜福利在线在线| 天美传媒精品一区二区| 极品教师在线视频| 人妻夜夜爽99麻豆av| 国产高清视频在线观看网站| 午夜爱爱视频在线播放| 久久精品国产亚洲av涩爱 | 久久久久网色| 亚洲在线观看片| 在线a可以看的网站| 亚洲国产欧美在线一区| 国模一区二区三区四区视频| 老女人水多毛片| 欧美激情国产日韩精品一区| 国产一区二区在线观看日韩| 久久6这里有精品| 中文在线观看免费www的网站| 在线天堂最新版资源| 久久亚洲精品不卡| 床上黄色一级片| 亚洲人成网站在线播| 亚洲经典国产精华液单| 欧美日韩在线观看h| 一级毛片久久久久久久久女| 亚洲久久久久久中文字幕| 欧美精品一区二区大全| 九色成人免费人妻av| 国产精品爽爽va在线观看网站| 美女cb高潮喷水在线观看| 六月丁香七月| 欧美激情在线99| 一本久久中文字幕| 女人被狂操c到高潮| 亚洲不卡免费看| 国产在线男女| 日韩成人伦理影院| 国产伦理片在线播放av一区 | 男人舔奶头视频| 国产真实伦视频高清在线观看| 99久久精品热视频| АⅤ资源中文在线天堂| 国产一区二区三区av在线 | 日日摸夜夜添夜夜爱| 久久久精品欧美日韩精品| 久久6这里有精品| 久久久久免费精品人妻一区二区| 国产精品一区二区性色av| av又黄又爽大尺度在线免费看 | 亚洲真实伦在线观看| 色播亚洲综合网| 久久草成人影院| 久久久久久久久久久免费av| 婷婷六月久久综合丁香| 国产成人福利小说| 一本精品99久久精品77| 国产成人freesex在线| 精品人妻视频免费看| 深爱激情五月婷婷| 日韩视频在线欧美| 只有这里有精品99| 麻豆成人午夜福利视频| av在线亚洲专区| 久久人人爽人人爽人人片va| 日韩一区二区三区影片| 伊人久久精品亚洲午夜| 国产精品三级大全| 成人二区视频| 国产精品久久久久久亚洲av鲁大| 卡戴珊不雅视频在线播放| 国产伦精品一区二区三区四那| 国产一区二区三区av在线 | 搞女人的毛片| 婷婷六月久久综合丁香| 又爽又黄无遮挡网站| 久久鲁丝午夜福利片| 久久久久久国产a免费观看| 美女 人体艺术 gogo| .国产精品久久| 日韩欧美三级三区| 性欧美人与动物交配| 欧美激情久久久久久爽电影| 只有这里有精品99| av视频在线观看入口| 少妇人妻精品综合一区二区 | 日韩精品有码人妻一区| 国产精品一二三区在线看| 久久久欧美国产精品| 国产高清视频在线观看网站| 51国产日韩欧美| 亚洲第一区二区三区不卡| 丰满乱子伦码专区| 一本精品99久久精品77| 亚洲精品久久久久久婷婷小说 | 最近中文字幕高清免费大全6| 亚洲美女搞黄在线观看| 少妇高潮的动态图| 婷婷色综合大香蕉| 亚洲av.av天堂| 97超视频在线观看视频| 91久久精品国产一区二区成人| 哪里可以看免费的av片| 一级毛片我不卡| 精品人妻一区二区三区麻豆| 成人高潮视频无遮挡免费网站| 日韩人妻高清精品专区| 免费av毛片视频| 69av精品久久久久久| 亚洲精品国产av成人精品| 在线播放国产精品三级| 国产精品.久久久| 九九久久精品国产亚洲av麻豆| 偷拍熟女少妇极品色| 一级av片app| 亚洲在线观看片| 免费看光身美女| 一级黄色大片毛片| 22中文网久久字幕| 深夜a级毛片| 九色成人免费人妻av| 午夜福利视频1000在线观看| 晚上一个人看的免费电影| 中文字幕av在线有码专区| 亚洲色图av天堂| av在线观看视频网站免费| 国产高清视频在线观看网站| 中文字幕熟女人妻在线| 最近手机中文字幕大全| 欧美极品一区二区三区四区| 久久亚洲国产成人精品v| 日日啪夜夜撸| 好男人视频免费观看在线| 看片在线看免费视频| 深夜a级毛片| 日本一二三区视频观看| 看片在线看免费视频| 亚洲国产欧美在线一区| 少妇熟女欧美另类| h日本视频在线播放| 有码 亚洲区| 国产精品久久久久久精品电影小说 | 青春草视频在线免费观看| 黄片无遮挡物在线观看| 国产真实乱freesex| 在线播放无遮挡| av国产免费在线观看| 亚洲国产精品成人综合色| 中文字幕制服av| 人体艺术视频欧美日本| 国产一区二区激情短视频| 可以在线观看的亚洲视频| 两性午夜刺激爽爽歪歪视频在线观看| 免费av毛片视频| 久久久久九九精品影院| 97超碰精品成人国产| 此物有八面人人有两片| 两个人视频免费观看高清| 久久人妻av系列| 变态另类成人亚洲欧美熟女| 在线观看午夜福利视频| 日韩欧美在线乱码| 草草在线视频免费看| 一级黄片播放器| 国产精品久久久久久亚洲av鲁大| 色播亚洲综合网| 国产在视频线在精品| 国产乱人视频| 久久韩国三级中文字幕| 亚洲丝袜综合中文字幕| 九九久久精品国产亚洲av麻豆| 久久久国产成人免费| 欧美+日韩+精品| 精品国产三级普通话版| 国产久久久一区二区三区| 日韩欧美 国产精品| av视频在线观看入口| 亚洲三级黄色毛片| 婷婷色综合大香蕉| 欧美激情在线99| 免费一级毛片在线播放高清视频| 久久久久久伊人网av| 少妇猛男粗大的猛烈进出视频 | 精品熟女少妇av免费看| 在线观看美女被高潮喷水网站| 热99re8久久精品国产| 亚洲av中文av极速乱| 日本爱情动作片www.在线观看| 国产亚洲av片在线观看秒播厂 | 99久久精品国产国产毛片| 亚洲人成网站在线观看播放| 久久精品久久久久久噜噜老黄 | 国产亚洲精品久久久com| 亚洲欧美成人精品一区二区| 中国国产av一级| 亚洲美女搞黄在线观看| 国产三级中文精品| 最近手机中文字幕大全| 国产精品电影一区二区三区| 久久精品综合一区二区三区| 99久久人妻综合| 亚洲欧美精品自产自拍| 国产黄片视频在线免费观看| 国产精品一区二区三区四区免费观看| 天美传媒精品一区二区| 高清毛片免费看| 国产三级在线视频| 免费观看的影片在线观看| 日本-黄色视频高清免费观看| 久久综合国产亚洲精品| 免费人成在线观看视频色| 搡女人真爽免费视频火全软件| 国产v大片淫在线免费观看| 日韩人妻高清精品专区| 欧美成人一区二区免费高清观看| 国产成人午夜福利电影在线观看| 99久久精品热视频| 99久久人妻综合| 最好的美女福利视频网| 久久韩国三级中文字幕| 联通29元200g的流量卡| 男女下面进入的视频免费午夜| 国产精品一区二区三区四区免费观看| 嫩草影院新地址| 青春草视频在线免费观看| 亚洲欧洲国产日韩| 99在线视频只有这里精品首页| 永久网站在线| 狠狠狠狠99中文字幕| 啦啦啦韩国在线观看视频| 欧美日韩在线观看h| 国产黄片视频在线免费观看| 久久精品国产99精品国产亚洲性色| 日韩国内少妇激情av| 插逼视频在线观看| 最新中文字幕久久久久| 嫩草影院入口| 亚洲精品乱码久久久v下载方式| 人妻制服诱惑在线中文字幕| 精品国产三级普通话版| 日韩av在线大香蕉| 亚洲成a人片在线一区二区| 国产欧美日韩精品一区二区| 日本黄色视频三级网站网址| 18禁黄网站禁片免费观看直播| 亚洲七黄色美女视频| 蜜桃久久精品国产亚洲av| 精品少妇黑人巨大在线播放 | 亚洲三级黄色毛片| ponron亚洲| 一个人免费在线观看电影| 精品久久久久久成人av| 少妇熟女aⅴ在线视频| 99热网站在线观看| 午夜福利成人在线免费观看| 免费看a级黄色片| 99在线视频只有这里精品首页| 网址你懂的国产日韩在线| 欧美性感艳星| 免费人成在线观看视频色| 91麻豆精品激情在线观看国产| 午夜福利视频1000在线观看| 亚洲激情五月婷婷啪啪| 午夜精品在线福利| 在线a可以看的网站| 精品免费久久久久久久清纯| 久久久精品大字幕| 亚洲av男天堂| 亚洲图色成人| 自拍偷自拍亚洲精品老妇| 国产麻豆成人av免费视频| 久久国产乱子免费精品| 最近的中文字幕免费完整| av.在线天堂| 亚洲五月天丁香| 国产伦一二天堂av在线观看| 毛片一级片免费看久久久久| 国产毛片a区久久久久| 日韩成人av中文字幕在线观看| 日本一本二区三区精品| 精品久久国产蜜桃| av在线老鸭窝| 丰满的人妻完整版| 日本黄色片子视频| 女同久久另类99精品国产91| 亚洲av熟女| 国产探花在线观看一区二区| 丰满的人妻完整版| 亚洲av不卡在线观看| 亚洲一区二区三区色噜噜| 中文字幕熟女人妻在线| 国产精品一区二区三区四区免费观看| 嫩草影院入口| 一个人看的www免费观看视频| 亚洲欧美日韩高清专用| 99视频精品全部免费 在线| 1024手机看黄色片| 欧美日本亚洲视频在线播放| 欧美色欧美亚洲另类二区| 最新中文字幕久久久久| 99久国产av精品| 国产蜜桃级精品一区二区三区| 精品日产1卡2卡| 久久精品国产清高在天天线| 国产综合懂色| 亚洲一级一片aⅴ在线观看| 精品久久久久久久久久久久久| 亚洲成人中文字幕在线播放| 国产成人a区在线观看| 久久婷婷人人爽人人干人人爱| 老师上课跳d突然被开到最大视频| 少妇被粗大猛烈的视频| 直男gayav资源| av免费观看日本| 成人鲁丝片一二三区免费| 亚洲乱码一区二区免费版| 一区二区三区四区激情视频 | 国产一区二区亚洲精品在线观看| 久久久久久久久久久丰满| 久久久久久久久大av| 天堂影院成人在线观看| 久久久精品欧美日韩精品| 可以在线观看毛片的网站| 搡女人真爽免费视频火全软件| 日韩一区二区视频免费看| 国产一区二区亚洲精品在线观看| 日本熟妇午夜| 亚洲av熟女| 中文字幕av在线有码专区| 亚洲国产精品成人综合色| 免费无遮挡裸体视频| 熟女人妻精品中文字幕| 麻豆久久精品国产亚洲av| 亚洲va在线va天堂va国产| 日产精品乱码卡一卡2卡三| 夜夜夜夜夜久久久久| www.色视频.com| 国产综合懂色| 亚洲美女视频黄频| 十八禁国产超污无遮挡网站| 国产精品久久久久久精品电影| 精品国内亚洲2022精品成人| 国产私拍福利视频在线观看| 成人二区视频| 丰满的人妻完整版| 久久草成人影院| 长腿黑丝高跟| 亚洲国产精品国产精品| 国产精品爽爽va在线观看网站| 麻豆一二三区av精品| 欧美性猛交黑人性爽| 亚洲在久久综合| 99久久无色码亚洲精品果冻| 国产精品综合久久久久久久免费| 国产黄色视频一区二区在线观看 | kizo精华| 六月丁香七月|