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

    情景分析在飽和負(fù)荷預(yù)測(cè)中的應(yīng)用

    2016-08-11 09:29:32陳屹東程浩忠黃錦華韓新陽(yáng)尤鐘曉

    陳屹東,程浩忠,黃錦華,韓新陽(yáng),尤鐘曉

    (1.上海交通大學(xué)電力傳輸與功率變換控制教育部重點(diǎn)實(shí)驗(yàn)室,上海 200240;2.浙江省電力公司經(jīng)濟(jì)技術(shù)研究院,杭州 310008;3.國(guó)網(wǎng)能源研究院,北京 100052)

    ?

    情景分析在飽和負(fù)荷預(yù)測(cè)中的應(yīng)用

    陳屹東1,程浩忠1,黃錦華2,韓新陽(yáng)3,尤鐘曉1

    (1.上海交通大學(xué)電力傳輸與功率變換控制教育部重點(diǎn)實(shí)驗(yàn)室,上海 200240;2.浙江省電力公司經(jīng)濟(jì)技術(shù)研究院,杭州 310008;3.國(guó)網(wǎng)能源研究院,北京 100052)

    摘 要:飽和負(fù)荷預(yù)測(cè)時(shí)間跨度長(zhǎng)、影響因素多,僅基于歷史數(shù)據(jù)進(jìn)行預(yù)測(cè)很可能造成較大誤差。因此,該文將情景分析思想引入飽和負(fù)荷預(yù)測(cè),情景分析法通過(guò)假設(shè)未來(lái)不同的發(fā)展情景,使預(yù)測(cè)結(jié)果能夠體現(xiàn)未來(lái)用電需求的可能出現(xiàn)的不同發(fā)展?fàn)顩r。針對(duì)負(fù)荷規(guī)模增長(zhǎng)呈“S”型的特點(diǎn),文中選取Logistic生長(zhǎng)曲線模型進(jìn)行預(yù)測(cè),分析了基于Logistic曲線的階段劃分方法,并且采用Levenberg-Marquardt算法確定模型的具體參數(shù)。通過(guò)引入情景分析的思想,也可以彌補(bǔ)Logistic曲線在基于歷史數(shù)據(jù)進(jìn)行擬合時(shí)只能得出一種預(yù)測(cè)方案的不足。實(shí)際算例驗(yàn)證了情景分析法在用于飽和負(fù)荷預(yù)測(cè)時(shí)的有效性。

    關(guān)鍵詞:飽和負(fù)荷;預(yù)測(cè);情景分析;負(fù)荷生長(zhǎng)曲線;發(fā)展階段劃分

    從發(fā)達(dá)國(guó)家的電力發(fā)展歷程來(lái)看,城市建設(shè)與發(fā)展初期,在經(jīng)濟(jì)高速發(fā)展與人口增長(zhǎng)的帶動(dòng)下,電力負(fù)荷呈快速增長(zhǎng)趨勢(shì)[1-2]。但當(dāng)城市發(fā)展到一定階段后,受到土地、交通、環(huán)境、自然資源及政策等因素的限制,電力負(fù)荷增長(zhǎng)速度也逐漸放緩,隨后負(fù)荷增長(zhǎng)率變得非常小甚至不再增長(zhǎng)或有所下降,最后達(dá)到飽和狀態(tài)。

    我國(guó)經(jīng)濟(jì)發(fā)展存在城鄉(xiāng)差異、地域差異,區(qū)域經(jīng)濟(jì)發(fā)展也存在很大的不平衡性。電力負(fù)荷與經(jīng)濟(jì)的發(fā)展有著密切的聯(lián)系,一些經(jīng)濟(jì)發(fā)展比較超前的地區(qū),其電力負(fù)荷增長(zhǎng)已開(kāi)始呈飽和趨勢(shì)。這些地區(qū)電網(wǎng)就需要綜合自身特點(diǎn)對(duì)其電力負(fù)荷的最終規(guī)模進(jìn)行預(yù)測(cè)來(lái)制定電網(wǎng)發(fā)展的遠(yuǎn)期規(guī)劃,因此電力飽和負(fù)荷的概念應(yīng)運(yùn)而生[3]。

    目前,對(duì)于飽和負(fù)荷的研究大多數(shù)是針對(duì)城市飽和負(fù)荷[4],對(duì)于區(qū)域級(jí)飽和負(fù)荷研究較少,而且對(duì)飽和負(fù)荷的概念也沒(méi)有明確、統(tǒng)一的定義。文獻(xiàn)[5]采用了一種改進(jìn)的Logistic參數(shù)估計(jì)方法,并驗(yàn)證了Logistic模型在用于飽和負(fù)荷預(yù)測(cè)時(shí)的良好預(yù)測(cè)效果。文獻(xiàn)[6]采用人均電量法預(yù)測(cè)了區(qū)域飽和負(fù)荷,該方法的優(yōu)點(diǎn)是模型簡(jiǎn)單,所需基礎(chǔ)數(shù)據(jù)少,但誤差也較大。文獻(xiàn)[7-9]根據(jù)對(duì)經(jīng)濟(jì)社會(huì)發(fā)展現(xiàn)狀的分析,采用負(fù)荷密度法對(duì)城市中心城區(qū)飽和負(fù)荷的總體水平及到達(dá)時(shí)間進(jìn)行預(yù)測(cè),這種方法能很好地反映出飽和時(shí)期城市各個(gè)區(qū)塊的負(fù)荷分布,但需要預(yù)測(cè)城市的具體發(fā)展定位及城區(qū)負(fù)荷密度數(shù)據(jù),不適用于以區(qū)域級(jí)負(fù)荷為對(duì)象的負(fù)荷預(yù)測(cè)。文獻(xiàn)[10]采用了系統(tǒng)動(dòng)力學(xué)結(jié)合計(jì)量經(jīng)濟(jì)學(xué)的方法來(lái)預(yù)測(cè)飽和負(fù)荷,但是這類方法由于預(yù)測(cè)所需基礎(chǔ)數(shù)據(jù)量大,在用于省級(jí)和區(qū)域級(jí)飽和負(fù)荷預(yù)測(cè)時(shí),很多基礎(chǔ)數(shù)據(jù)難以獲取。文獻(xiàn)[11]基于智能電網(wǎng)環(huán)境下的高級(jí)量測(cè)技術(shù),基于每小時(shí)的數(shù)據(jù)預(yù)測(cè)長(zhǎng)期負(fù)荷,能有效提高預(yù)測(cè)精度。文獻(xiàn)[12-13]將組合預(yù)測(cè)理論用于中長(zhǎng)期負(fù)荷預(yù)測(cè)之中,也為飽和負(fù)荷預(yù)測(cè)提供了新的思路。

    生長(zhǎng)曲線法預(yù)測(cè)所需原始數(shù)據(jù)少、計(jì)算過(guò)程簡(jiǎn)單,也能從整體上較好地把握電力負(fù)荷的發(fā)展趨勢(shì)。而負(fù)荷密度法或系統(tǒng)動(dòng)力學(xué)法需要收集人口、經(jīng)濟(jì)、氣候、產(chǎn)業(yè)結(jié)構(gòu)、負(fù)荷密度等大量基礎(chǔ)數(shù)據(jù),雖然預(yù)測(cè)結(jié)果能體現(xiàn)負(fù)荷的分布情況以及各隨機(jī)因素對(duì)負(fù)荷波動(dòng)的影響,但是預(yù)測(cè)所需的很多基礎(chǔ)數(shù)據(jù)難以準(zhǔn)確獲取,操作較為困難。對(duì)于區(qū)域級(jí)飽和負(fù)荷預(yù)測(cè)而言,首先也需要從整體上把握其發(fā)展趨勢(shì),因此生長(zhǎng)曲線法操作更為簡(jiǎn)便。本文分析了基于Logistic生長(zhǎng)曲線的發(fā)展階段劃分方法,并采用Levenberg-Marquardt算法求取曲線的待定參數(shù)和特征點(diǎn),該方法具有梯度下降法和高斯-牛頓法的優(yōu)點(diǎn),具有更好的收斂性,大大降低了陷入局部最優(yōu)的可能性。

    與傳統(tǒng)的針對(duì)具體年限的負(fù)荷預(yù)測(cè)不同,飽和負(fù)荷預(yù)測(cè)時(shí)間跨度大,涉及面廣,影響因素眾多。無(wú)論采用哪種預(yù)測(cè)方法,都只是根據(jù)歷史數(shù)據(jù)預(yù)測(cè)未來(lái)的最可能發(fā)生的某一種情況,如果未來(lái)外部條件如政策、資源等發(fā)生較大變化,很可能導(dǎo)致最終預(yù)測(cè)結(jié)果出現(xiàn)較大誤差。因此,本文引入了情景分析法,分析未來(lái)可能出現(xiàn)的不同狀況,并基于不同的未來(lái)情景分別進(jìn)行預(yù)測(cè),這樣便可彌補(bǔ)運(yùn)用生長(zhǎng)曲線法預(yù)測(cè)時(shí)只能得出一種預(yù)測(cè)方案的不足。而我國(guó)大部分區(qū)域電力都還處于高速發(fā)展期,未來(lái)發(fā)展不確定因素眾多,所以引入情景分析法對(duì)于中國(guó)未來(lái)飽和預(yù)測(cè)也具有一定的指導(dǎo)意義。

    1 負(fù)荷生長(zhǎng)曲線

    1.1飽和負(fù)荷基本概念

    一個(gè)地區(qū)的負(fù)荷增長(zhǎng)通常呈如下規(guī)律:在經(jīng)濟(jì)發(fā)展初期,用電需求增長(zhǎng)較慢;隨著經(jīng)濟(jì)的快速發(fā)展,用電需求出現(xiàn)快速增長(zhǎng)趨勢(shì);但是受土地面積、人口、資源等因素制約,用電需求不會(huì)無(wú)限制地增長(zhǎng),而是呈現(xiàn)出飽和增長(zhǎng)態(tài)勢(shì),即總體呈現(xiàn)“S”型。

    飽和負(fù)荷涵蓋了兩方面內(nèi)容,一是某地區(qū)全社會(huì)用電量達(dá)到飽和時(shí)的規(guī)模,二是某地區(qū)年最大負(fù)荷達(dá)到飽和時(shí)的規(guī)模。

    1.2基于Logistic曲線的發(fā)展階段劃分

    Logistic曲線是一種最常用的“S”型曲線模型。該模型首先源于對(duì)生物種群的研究,通過(guò)大量的科學(xué)觀察研究,發(fā)現(xiàn)經(jīng)濟(jì)、社會(huì)、人口以及科技領(lǐng)域中很多事物的成長(zhǎng)過(guò)程,以及事物成長(zhǎng)過(guò)程中的某個(gè)或某些定量化特性也符合生物成長(zhǎng)過(guò)程中的曲線規(guī)律。該曲線模型的基本方程為

    該曲線可以分成3部分,首先是比較緩慢的增長(zhǎng),而后快速的增長(zhǎng),最后增長(zhǎng)速度變得緩和而逐漸趨向于穩(wěn)定。曲線穩(wěn)定的極限水平漸近線是l

    ti→m∞y=c。

    對(duì)式(1)分別求取二階和三階導(dǎo)數(shù)得

    令式(2)為0可得

    令式(3)為0可得

    在T1、T2和T3點(diǎn)上,Logistic函數(shù)值分別為

    圖1 Logistic曲線階段劃分示意Fig.1 Stage division of Logistic curve

    由圖1可知,Logistic函數(shù)單調(diào)增長(zhǎng),其增長(zhǎng)速度的最大值對(duì)應(yīng)于特征點(diǎn)P2,),該點(diǎn)也是曲線的對(duì)稱中心點(diǎn)。特征點(diǎn)P1、P2、P3將Logistic曲線分為4段,可分別解釋為初始發(fā)展階段(0~T1)、快速發(fā)展階段(T1~T2)、后發(fā)展階段(T2~T3)、飽和發(fā)展階段(T3~+∞),且特征點(diǎn)P1、P2、P3對(duì)應(yīng)的函數(shù)值y1、y2、y3僅與曲線最終飽和值c有關(guān),與參數(shù)a、b均無(wú)關(guān)。

    如果預(yù)測(cè)模型可變參數(shù)太少,則模型自由度低,擬合效果可能很差;擴(kuò)展模型可變參數(shù)數(shù)目,可增加模型自由度,提高參數(shù)估計(jì)效果和歷史數(shù)據(jù)擬合精度[14]。擴(kuò)展后的Logistic函數(shù)方程為

    1.3飽和負(fù)荷判定指標(biāo)

    發(fā)達(dá)國(guó)家和地區(qū)的用電需求及對(duì)應(yīng)的社會(huì)經(jīng)濟(jì)背景可以對(duì)我國(guó)電力負(fù)荷飽和階段的界定提供重要參考和借鑒,表1列出了部分發(fā)達(dá)國(guó)家和地區(qū)用電需求飽和時(shí)的經(jīng)濟(jì)、電力相關(guān)數(shù)據(jù)。

    表1 部分發(fā)達(dá)國(guó)家和地區(qū)負(fù)荷飽和時(shí)的指標(biāo)Tab.1 Indicators on saturation stage in some developed countries and regions

    參考發(fā)達(dá)國(guó)家和地區(qū)的用電需求增長(zhǎng)率、人均用電量、人口增長(zhǎng)率、人均GDP、產(chǎn)業(yè)結(jié)構(gòu)等各項(xiàng)指標(biāo),得到負(fù)荷進(jìn)入飽和階段的參考判定指標(biāo),如表2所示。判定指標(biāo)可分為必要指標(biāo)和輔助指標(biāo)兩類,可根據(jù)這些指標(biāo)對(duì)飽和負(fù)荷預(yù)測(cè)結(jié)果進(jìn)行校驗(yàn),其中必要指標(biāo)為在校驗(yàn)時(shí)必須滿足的指標(biāo),也是各發(fā)達(dá)國(guó)家和地區(qū)在電力負(fù)荷達(dá)到飽和時(shí)共同滿足的指標(biāo)。而輔助指標(biāo)根據(jù)各個(gè)國(guó)家和地區(qū)的具體國(guó)情、資源稟賦、地理位置等因素不同而存在一定差異,在對(duì)預(yù)測(cè)結(jié)果進(jìn)行校驗(yàn)時(shí)某地區(qū)只需滿足其中2個(gè)指標(biāo)即可認(rèn)為電力負(fù)荷已經(jīng)達(dá)到飽和。

    表2 負(fù)荷進(jìn)入飽和階段的參考判定指標(biāo)Tab.2 Indicators on saturation stage

    2 參數(shù)估計(jì)的Levenberg-Marquardt算法

    非線性函數(shù)的一般表達(dá)形式為

    式中:f為非線性函數(shù);x1,x2,…,xm為m個(gè)自變量;a1,a2,…,ap為p個(gè)待估未知參數(shù);ε為隨機(jī)誤差項(xiàng)。對(duì)y和x1,x2,…,xm通過(guò)n次觀察得到n組觀察值:(xi1,xi2,…,xim,yi),i=1,2,…,n。

    將自變量的第i次觀察值代入函數(shù)f(x1,x2,…,xm;a1,a2,…,ap)=f(xi,a),給a賦初始值a(0)=(a,a,…,a),然后將 f(xi,a)在a(0)處按泰勒級(jí)數(shù)展開(kāi),略去二次及其以上項(xiàng),可得

    根據(jù)最小二乘法,令

    其中d≥0為阻尼因子,當(dāng)取d=0時(shí),就是高斯-牛頓法,高斯-牛頓法就是Levenberg-Marquadt算法的特殊形式,它對(duì)迭代初始值的選擇更為嚴(yán)格。

    令Q分別對(duì)a1,a2,…,ap的一階偏導(dǎo)數(shù)為0,即,可得

    式(11)可以轉(zhuǎn)化為

    其中:

    從而可得

    若解得aj與a)之差的絕對(duì)值很小,則可認(rèn)為估計(jì)成功。如果較大,則把上一步算得的aj作為新的a代入式(12),得到新的aj又作為a代入式(12),如此反復(fù)迭代直到aj與a之差可以忽略為止[15]。

    對(duì)非線性最小二乘問(wèn)題,Levenberg-Marquadt是一種經(jīng)典的優(yōu)化迭代求解算法,它通過(guò)在求解被估計(jì)的函數(shù)參數(shù)矩陣的迭代過(guò)程中增加阻尼項(xiàng),改善法方程的正則化條件,克服法方程系數(shù)矩陣出現(xiàn)病態(tài)或奇異的不利影響,使函數(shù)模型獲得穩(wěn)定解。

    3 情景分析法

    情景(scenario)一詞最早由英國(guó)的Kahn等[16]提出,他們認(rèn)為未來(lái)是多樣的,通向這種或那種未來(lái)結(jié)果的途徑的描述也不是唯一的,對(duì)可能出現(xiàn)的未來(lái)以及實(shí)現(xiàn)這種未來(lái)的途徑的假設(shè)就構(gòu)成了一個(gè)情景?;谇榫暗母拍睿榫胺治鼍褪欠治?、構(gòu)建未來(lái)情景的過(guò)程。情景分析法(scenario analysis method)是在推測(cè)的基礎(chǔ)之上,對(duì)未來(lái)可能的情景加以描述和預(yù)測(cè)的方法。情景分析法通過(guò)對(duì)環(huán)境的研究,識(shí)別影響研究主體發(fā)展的外部因素,對(duì)不同政策、規(guī)劃、技術(shù)等對(duì)未來(lái)發(fā)展趨勢(shì)產(chǎn)生的影響和效果進(jìn)行定性與定量相結(jié)合的系統(tǒng)分析,模擬外部因素可能發(fā)生的多種交叉情景和各種可能前景,并以此來(lái)分析與預(yù)測(cè)影響未來(lái)發(fā)展的各種不確定因素的相互作用[17]。

    情景分析法在進(jìn)行預(yù)測(cè)時(shí),不僅可根據(jù)預(yù)測(cè)對(duì)象的內(nèi)在產(chǎn)生機(jī)理從定量方法上進(jìn)行推理與歸納,還可對(duì)各不確定因素(自變量)的幾種典型的可能情況采取人為決策,從而可以更合理地模擬現(xiàn)實(shí)。

    此外,情景分析法與傳統(tǒng)預(yù)測(cè)方法相比還有一個(gè)顯著的不同點(diǎn)。傳統(tǒng)預(yù)測(cè)方法是試圖描繪被預(yù)測(cè)對(duì)象未來(lái)最可能發(fā)生的某一種狀況,及這種可能程度的大小。而情景分析法則采用的是一種多路徑式的預(yù)測(cè)方式,研究各種假設(shè)條件下的被預(yù)測(cè)對(duì)象未來(lái)可能出現(xiàn)何種情況。因此在某種意義上,傳統(tǒng)預(yù)測(cè)方法也可以看作是情景分析法的一種特例。在情景分析中,雖然各種假設(shè)條件不一定會(huì)出現(xiàn),但通過(guò)這樣的分析,可幫助人們了解如果要讓被研究對(duì)象出現(xiàn)某種結(jié)果需要采取的前提條件是什么。單一預(yù)測(cè)方法和情景分析法的區(qū)別如圖2所示[18]。

    圖2 單一預(yù)測(cè)和情景分析法的比較Fig.2 Comparison between single prediction and scenario analysis

    4 飽和負(fù)荷預(yù)測(cè)流程

    本文采用的基于情景分析法的飽和負(fù)荷預(yù)測(cè)基本步驟如下。

    步驟1基準(zhǔn)情景設(shè)定。

    基準(zhǔn)情景主要是基于歷史數(shù)據(jù)進(jìn)行外推,該情景作為與其他情景進(jìn)行對(duì)比的基礎(chǔ)。

    步驟2確定電力負(fù)荷發(fā)展階段劃分方案。

    采用參數(shù)擴(kuò)展后的Logistic模型,即式(7)對(duì)用電量及最高負(fù)荷序列進(jìn)行建模分析,用Levenberg-Marquardt算法對(duì)曲線待定參數(shù)進(jìn)行估計(jì)。根據(jù)得到的參數(shù)求取曲線的3個(gè)特征時(shí)間點(diǎn),并得出階段劃分方案。

    步驟3飽和負(fù)荷時(shí)間點(diǎn)和飽和規(guī)模預(yù)測(cè)。

    用Logistic曲線分別對(duì)用電量和最高負(fù)荷序列進(jìn)行分析預(yù)測(cè),首先取曲線極值c的95%對(duì)應(yīng)的值和年份作為飽和值,然后用判定指標(biāo)進(jìn)行校驗(yàn),若有指標(biāo)未達(dá)標(biāo),則將年份推后一年繼續(xù)計(jì)算,直到各項(xiàng)指標(biāo)都滿足要求為止。根據(jù)用電量和最高負(fù)荷的飽和年份得出達(dá)到飽和的時(shí)間范圍。

    步驟4其他情景設(shè)定。

    基于基準(zhǔn)情景下的預(yù)測(cè)結(jié)果,設(shè)定其他的發(fā)展情景,根據(jù)不同情景下的限制條件,調(diào)整模型參數(shù)值,再次進(jìn)行預(yù)測(cè)?;谇榫胺治龅娘柡拓?fù)荷分析思路和步驟如圖3所示。

    圖3 基于情景分析法的飽和負(fù)荷預(yù)測(cè)流程Fig.3 Flow chart of saturation load forecasting based on scenario analysis

    5 算例分析

    為了說(shuō)明基于情景分析的飽和預(yù)測(cè)思路和方法,以某省電網(wǎng)1990—2013年的實(shí)際全社會(huì)用電量、年最高負(fù)荷數(shù)據(jù)為例,通過(guò)設(shè)置不同發(fā)展情景,對(duì)該省的飽和負(fù)荷進(jìn)行預(yù)測(cè),并對(duì)預(yù)測(cè)結(jié)果進(jìn)行評(píng)價(jià)分析。

    設(shè)定3種不同的發(fā)展情景,即基準(zhǔn)情景、高速發(fā)展情景和低速發(fā)展情景?;鶞?zhǔn)情景主要是考慮過(guò)去的發(fā)展特點(diǎn),該情景作為與其他情景進(jìn)行對(duì)比的基礎(chǔ)。高速發(fā)展情景主要是考慮一些可以預(yù)見(jiàn)的因素,如人口、技術(shù)進(jìn)步、電氣化水平提高等變化可能導(dǎo)致未來(lái)電力負(fù)荷會(huì)保持更長(zhǎng)時(shí)間的快速增長(zhǎng),從而使得電力需求飽和規(guī)模也相應(yīng)增加,文獻(xiàn)[19]指出美國(guó)通過(guò)發(fā)展電動(dòng)汽車(chē),預(yù)計(jì)到2030年增加的電力需求比不考慮電動(dòng)汽車(chē)的情景要高5%~6%。低速發(fā)展情景則是考慮以后經(jīng)濟(jì)發(fā)展因受較多不利因素影響而面臨一定風(fēng)險(xiǎn),從而可能會(huì)導(dǎo)致電力發(fā)展速度變慢,電力需求飽和規(guī)模也會(huì)相應(yīng)降低。本文假設(shè)高速發(fā)展情景和低速發(fā)展情景分別是在基準(zhǔn)情景下飽和量的基礎(chǔ)上增加10%和降低10%,在不同發(fā)展情景下,分別得到全社會(huì)用電量和年最高負(fù)荷的擬合曲線如圖4和圖5所示。

    圖4 3種情景下全社會(huì)用電量擬合預(yù)測(cè)Fig.4 Electricity consumption fitting and forecasting under three different scenario

    圖5 3種情景下年最大負(fù)荷擬合預(yù)測(cè)Fig.5 Annual peak load fitting and forecasting under three different scenario

    在3種情景下分別得到的預(yù)測(cè)結(jié)果和發(fā)展階段劃分時(shí)間節(jié)點(diǎn)如表3和表4所示。其中情景1、情景2、情景3分別對(duì)應(yīng)基準(zhǔn)情景、高速發(fā)展情景和低速發(fā)展情景。

    從預(yù)測(cè)結(jié)果中可以看出,當(dāng)運(yùn)用Logistic模型進(jìn)行飽和負(fù)荷預(yù)測(cè)和階段劃分時(shí),雖然在不同情景下飽和規(guī)模不同,但飽和時(shí)間點(diǎn)和階段劃分時(shí)間節(jié)點(diǎn)都變化不大。說(shuō)明基于Logistic曲線進(jìn)行階段劃分時(shí),特征時(shí)間點(diǎn)對(duì)飽和量的變化敏感度低,即不同情景下雖然飽和規(guī)模發(fā)生變化,但是其對(duì)于合理推斷電力的發(fā)展階段影響不太大。

    表3 不同發(fā)展情景下的特征時(shí)間點(diǎn)和飽和規(guī)模Tab.3 Feature points and saturation scale under different scenarios

    不同情景下部分年份的預(yù)測(cè)值和相對(duì)誤差如表5所示。采用灰色殘差檢驗(yàn)法對(duì)建模精度進(jìn)行檢驗(yàn),設(shè)X為實(shí)際值,X?為預(yù)測(cè)值,首先計(jì)算相對(duì)誤差絕對(duì)值的平均值,即,則建模精度為 p=(1-εˉ)×100%。一般要求p>80%,最好是p>90%。對(duì)全社會(huì)用電量而言,3種情景下分別為1=0.059,ˉ2=0.058,ˉ3=0.058,其建模精度分別為p1=94.1%,p2=94.2%,p3=94.2%。對(duì)最高負(fù)荷而言,3種情景下相對(duì)誤差絕對(duì)值的平均值分別為1=0.061,ˉ2=0.063,ˉ3=0.059,建模精度分別為p1=93.9%,p2=93.7%,p3=94.1%。可見(jiàn)不同情景下模型對(duì)全社會(huì)用電量和最高負(fù)荷均具有較好的建模精度。

    表4 不同發(fā)展情景下的飽和負(fù)荷預(yù)測(cè)結(jié)果Tab.4 Forecasting results under different scenarios

    表5 不同情景下的預(yù)測(cè)誤差分析Tab.5 Error analysis under different scenarios

    6 結(jié)論

    (1)由于飽和負(fù)荷預(yù)測(cè)時(shí)間跨度大,影響因素多,單一的預(yù)測(cè)方法僅根據(jù)歷史數(shù)據(jù)進(jìn)行外推可能導(dǎo)致預(yù)測(cè)結(jié)果出現(xiàn)較大偏差。因此,本文將情景分析思想引入飽和負(fù)荷預(yù)測(cè)。情景分析法在進(jìn)行預(yù)測(cè)時(shí),不僅可根據(jù)預(yù)測(cè)對(duì)象的內(nèi)在產(chǎn)生機(jī)理從定量方法上進(jìn)行推理與歸納,還可對(duì)各不確定因素的幾種典型的可能情況采取人為決策,從而更合理地模擬現(xiàn)實(shí)。

    (2)雖然傳統(tǒng)的預(yù)測(cè)方法也可以通過(guò)采用不同預(yù)測(cè)模型得出高中低等不同方案,但是總體思路仍是以歷史數(shù)據(jù)為基礎(chǔ),沒(méi)有充分考慮未來(lái)可能發(fā)生的情況,而且采用的不同預(yù)測(cè)模型不一定都適合飽和負(fù)荷預(yù)測(cè),尤其是對(duì)于處在快速發(fā)展時(shí)期的地區(qū)來(lái)說(shuō),未來(lái)發(fā)展不確定因素眾多,因此情景分析法可以彌補(bǔ)傳統(tǒng)預(yù)測(cè)方法的不足。

    (3)從預(yù)測(cè)結(jié)果可以看出,當(dāng)用Logistic曲線模型進(jìn)行預(yù)測(cè)時(shí),在不同的發(fā)展情景下雖然飽和規(guī)模不同,但是達(dá)到飽和的年份以及劃分發(fā)展階段的時(shí)間節(jié)點(diǎn)都變化不大。

    (4)飽和負(fù)荷的預(yù)測(cè)是一項(xiàng)長(zhǎng)期工作,只有不斷更新用電需求與社會(huì)經(jīng)濟(jì)發(fā)展相關(guān)數(shù)據(jù),并根據(jù)最新數(shù)據(jù)對(duì)未來(lái)情景進(jìn)行滾動(dòng)修正,才能得到更加科學(xué)、合理的結(jié)論。

    參考文獻(xiàn):

    [1] Fumio A.History of power system development in Japan [C]//IEEE Conference on the History of Electric Power. Newark,USA,2007:1-9.

    [2]Elders I M,Ault G W,Galloway S.Identification of longterm scenarios of electricity network development[C]//In?ternational Conference on Future Power Systems.Amster?dam,Netherlands,2005.

    [3]崔凱,李敬如,劉海波,等(Cui Kai,Li Jingru,Liu Haibo,et al).城市負(fù)荷飽和階段電力規(guī)劃方法及其在濟(jì)南電網(wǎng)中的應(yīng)用(City′s power planning methods at the stage of load saturation and its application in Jinan power grid)[J].電網(wǎng)技術(shù)(Power System Technology),2007,31(S2):131-134.

    [4]肖欣,周渝慧,張寧,等(Xiao Xin,Zhou Yuhui,Zhang Ning,et al).城市電力飽和負(fù)荷分析技術(shù)及其應(yīng)用研究綜述(Survey of saturated load analysis technology for ur?ban power system and its application)[J].電力自動(dòng)化設(shè)備(Electric Power Automation Equipment),2014,34(6):146-152.

    [5]Jia Yudong,Li Shenghu,Tan Yun,et al.Improved para?metric estimation of logistic model for saturated load fore?cast[C]//Asia-Pacific Power and Energy Engineering Con?ference.Shanghai,China,2012.

    [6]王偉,房婷婷(Wang Wei,F(xiàn)ang Tingting).人均用電量法在區(qū)域飽和負(fù)荷預(yù)測(cè)中的應(yīng)用研究(The application of per-person electricity consumption method in saturation load forecasting)[J].電力需求側(cè)管理(Power Demand Side Management),2012,14(1):21-23.

    [7]崔凱,張麗娟,李敬如,等(Cui Kai,Zhang Lijuan,Li Jin?gru,et al).天津市中心城區(qū)飽和負(fù)荷分析與預(yù)測(cè)(Analysis and forecast of saturated load for the central city district of Tianjin)[J].電力技術(shù)經(jīng)濟(jì)(Electric Power Technologic Economics),2008,20(5):32-36.

    [8]王晶,馮顯時(shí),郭紅珍(Wang Jing,F(xiàn)eng Xianshi,Guo Hongzhen).基于蟻群元胞自動(dòng)機(jī)理論的城市飽和負(fù)荷預(yù)測(cè)(Urban load saturation forecast based on ant cellular automata theory)[J].中國(guó)電力(Electric Power),2011,44 (7):17-20.

    [9]Vasquez-Arnez R L,Jardini J A,Casolari R,et al.A meth?odology for electrical energy forecast and its spatial alloca?tion over developing boroughs[C]//IEEE/PES Transmis?sion and Distribution Conference&Exposition.Chicago,USA,2008:1346-1351.

    [10]何永秀,吳良器,戴愛(ài)英,等(He Yongxiu,Wu Liangqi,Dai Aiying,et al).基于系統(tǒng)動(dòng)力學(xué)與計(jì)量經(jīng)濟(jì)模型的城市飽和負(fù)荷綜合預(yù)測(cè)方法(Combined saturation load forecast model based on system dynamics and economet?rics)[J].電力需求側(cè)管理(Power Demand Side Manage? ment),2010,13(1):21-25.

    [11]Hong Tao,Wilson Jason,Xie Jingrui.Long term probabi?listic load forecasting and normalization with hourly infor?mation[J].IEEE Trans on Smart Grid,2014,5(1):456-462.

    [12]金鑫,羅滇生,孫廣強(qiáng),等(Jin Xin,Luo Diansheng,Sun Guangqiang,et al).中長(zhǎng)期電力負(fù)荷預(yù)測(cè)模型篩選與組合方法(Shifing and combination method of medium-andlong-term load forecasting model)[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào)(Proceedings of the CSU-EPSA),2012,24(4):150-156.

    [13]尚芳屹(Shang Fangyi).組合預(yù)測(cè)在區(qū)域級(jí)飽和負(fù)荷預(yù)測(cè)中的應(yīng)用(Application of a Combination Model in Re?gional Saturation Load Forecasting)[D].上海:上海交通大學(xué)電子信息與電氣工程學(xué)院(Shanghai:School of Electronic Information and Electrical Engineering,Shang?hai Jiao Tong University),2013.

    [14]寧波,康重慶,夏清(Ning Bo,Kang Chongqing,Xia Qing).中長(zhǎng)期負(fù)荷預(yù)測(cè)模型的擴(kuò)展策略(The expan?sion strategy for medium and long term load forecasting model)[J].中國(guó)電力(Electric Power),2000,33(10):36-38,108.

    [15]孫鵬,陳紹輝,張彩慶(Sun Peng,Chen Shaohui,Zhang Caiqing).基于Marquardt法參數(shù)估計(jì)的變電設(shè)備壽命周期故障率評(píng)估(Assessment of failure rate for substa?tion equipment life cycle based on Marquardt parameter estimation method)[J].電力系統(tǒng)保護(hù)與控制(Power Sys?tem Protection and Control),2012,40(1):85-90.

    [16]Kahn H,Wiener A J.The Year 2000:A Framework for Speculation on the Next Thirty-Three Years[M].New York:Macmillan,1967.

    [17]夏德建(Xia Dejian).基于情景分析的發(fā)電側(cè)碳排放生命周期計(jì)量研究(The Research on Life Cycle Carbon Emissions Measurement of Electric Power Generation Side Based on the Scenario Analysis Method)[D].重慶:重慶大學(xué)經(jīng)濟(jì)與工商管理學(xué)院(Chongqing:College of Economics and Business Administration,Chongqing Uni?versity),2010.

    [18]于紅霞(Yu Hongxia).情景分析在港口發(fā)展戰(zhàn)略中的應(yīng)用研究(Research on Scenario Analysis of Port Develop?ment Strategy)[D].天津:天津大學(xué)建筑工程學(xué)院(Tian?jin:School of Civil Engineering,Tianjin University),2004.

    [19]劉振亞.中國(guó)電力與能源[M].北京:中國(guó)電力出版社,2012.

    陳屹東(1990—),男,碩士研究生,研究方向?yàn)殡娋W(wǎng)規(guī)劃與優(yōu)化、電力系統(tǒng)負(fù)荷預(yù)測(cè)。Email:cydt2009@gmail.com

    程浩忠(1962—),男,博士,教授,研究方向?yàn)殡娋W(wǎng)規(guī)劃、電壓穩(wěn)定性、電能質(zhì)量。Email:hzcheng@sjtu.edu.cn

    黃錦華(1979—),男,碩士,高級(jí)工程師,研究方向?yàn)殡娏κ袌?chǎng)分析預(yù)測(cè)、電力系統(tǒng)規(guī)劃等。Email:huang_jinhua@sina. com

    中圖分類號(hào):TM715

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1003-8930(2016)06-0007-07

    DOI:10.3969/j.issn.1003-8930.2016.06.002

    作者簡(jiǎn)介:

    收稿日期:2014-09-22;修回日期:2015-11-25

    基金項(xiàng)目:國(guó)家自然科學(xué)基金重點(diǎn)資助項(xiàng)目(51337005)

    Application of Scenario Analysis in Saturation Load Forecasting

    CHEN Yidong1,CHENG Haozhong1,HUANG Jinhua2,HAN Xinyang3,YOU Zhongxiao1
    (1.Key Laboratory of Control of Power Transmission and Conversion,Ministry of Education,Shanghai Jiao Tong University,Shanghai 200240,China;2.Zhejiang Electric Power Company Economic Research Institute,Hangzhou 310008,China;3.State Grid Energy Research Institute,Beijing 100052,China)

    Abstract:Saturation load forecasting involves long span of time and many influencing factors,and predictions simply based on historical data may probably cause error.Thus scenario analysis is introduced into saturation load forecasting in this paper.The prediction results based on scenario analysis can reflect different possible future scenarios of the elec?tricity demand.As load growth presents‘S’type characteristic,Logistic growth curve model is selected for prediction and the phase division based on Logistic curve is analyzed in this paper.The specific parameters of the model can be de?termined by Levenberg-Marquardt algorithm.Through introducing scenario analysis,multiple prediction results can be acquired by Logistic model.Finally,the validity of the scenario analysis method is verified by a practical example.

    Key words:saturation load;forecasting;scenario analysis;load growth curve;development phase division

    久久久精品94久久精品| 中国国产av一级| 亚洲精品色激情综合| 久久久久久久久久黄片| 亚洲国产欧洲综合997久久,| 亚洲av电影在线观看一区二区三区 | 久久人妻av系列| 欧美日本视频| 少妇高潮的动态图| 日韩欧美精品免费久久| 高清日韩中文字幕在线| 久久99精品国语久久久| 国产成人a区在线观看| 亚洲综合精品二区| 国语对白做爰xxxⅹ性视频网站| 久久99热这里只频精品6学生 | 午夜福利视频1000在线观看| 亚洲精品国产av成人精品| 久久国内精品自在自线图片| 中文乱码字字幕精品一区二区三区 | 免费电影在线观看免费观看| 国产欧美另类精品又又久久亚洲欧美| 日韩成人av中文字幕在线观看| 热99re8久久精品国产| 国产色婷婷99| 国产极品精品免费视频能看的| 九九久久精品国产亚洲av麻豆| 久久精品久久精品一区二区三区| 亚洲av电影不卡..在线观看| 亚洲av二区三区四区| 最新中文字幕久久久久| 亚洲高清免费不卡视频| videos熟女内射| h日本视频在线播放| 国产免费一级a男人的天堂| 男女啪啪激烈高潮av片| 亚洲av成人av| 天美传媒精品一区二区| 亚洲久久久久久中文字幕| 中国美白少妇内射xxxbb| av在线观看视频网站免费| 国产精品爽爽va在线观看网站| 亚洲国产色片| 久久久欧美国产精品| 一区二区三区高清视频在线| 淫秽高清视频在线观看| 国产v大片淫在线免费观看| 国产伦理片在线播放av一区| 精华霜和精华液先用哪个| 我要搜黄色片| 亚洲精品乱码久久久v下载方式| 女人久久www免费人成看片 | 午夜精品一区二区三区免费看| 人体艺术视频欧美日本| 国产高清视频在线观看网站| 午夜精品在线福利| 成人午夜精彩视频在线观看| ponron亚洲| АⅤ资源中文在线天堂| 久99久视频精品免费| 色播亚洲综合网| 亚洲av成人av| 桃色一区二区三区在线观看| 日本熟妇午夜| 精品一区二区免费观看| av国产久精品久网站免费入址| 国产v大片淫在线免费观看| 人妻制服诱惑在线中文字幕| 色综合亚洲欧美另类图片| 六月丁香七月| 天堂网av新在线| 永久网站在线| 亚洲精品国产成人久久av| 最近2019中文字幕mv第一页| 国产午夜精品一二区理论片| 亚洲国产精品sss在线观看| 综合色丁香网| 欧美极品一区二区三区四区| 亚洲av二区三区四区| 亚洲av男天堂| 乱系列少妇在线播放| 久久这里有精品视频免费| 午夜福利高清视频| 国产视频首页在线观看| 女的被弄到高潮叫床怎么办| 九九在线视频观看精品| 亚洲av一区综合| 亚洲欧美中文字幕日韩二区| 久久久欧美国产精品| av在线观看视频网站免费| 少妇裸体淫交视频免费看高清| 久久99热6这里只有精品| 欧美日韩国产亚洲二区| 五月伊人婷婷丁香| 国模一区二区三区四区视频| 国产高清三级在线| 天天躁夜夜躁狠狠久久av| 欧美一区二区国产精品久久精品| 嘟嘟电影网在线观看| 内射极品少妇av片p| 国产三级在线视频| 久久精品久久精品一区二区三区| 亚洲av成人av| 美女cb高潮喷水在线观看| 成人毛片60女人毛片免费| 日韩亚洲欧美综合| 国产极品天堂在线| 2022亚洲国产成人精品| 久久久成人免费电影| 久久人人爽人人爽人人片va| 亚洲av成人av| 99九九线精品视频在线观看视频| 女人十人毛片免费观看3o分钟| 国产黄片视频在线免费观看| 精品久久久久久久末码| 秋霞在线观看毛片| 人妻系列 视频| 国产亚洲精品av在线| 国产 一区精品| 毛片一级片免费看久久久久| 久久久久久久亚洲中文字幕| 久久99精品国语久久久| 丰满少妇做爰视频| 午夜爱爱视频在线播放| 能在线免费观看的黄片| 精品久久久久久久人妻蜜臀av| 日日啪夜夜撸| 狠狠狠狠99中文字幕| 亚洲图色成人| 在线观看av片永久免费下载| av在线天堂中文字幕| 亚洲精品456在线播放app| 99久久九九国产精品国产免费| 国产伦精品一区二区三区视频9| 午夜福利在线观看吧| 天堂网av新在线| 嫩草影院入口| 久热久热在线精品观看| 日本-黄色视频高清免费观看| 亚洲av福利一区| 99久国产av精品国产电影| 亚洲丝袜综合中文字幕| www日本黄色视频网| 好男人视频免费观看在线| 日本一本二区三区精品| av在线观看视频网站免费| 久久精品久久精品一区二区三区| 日本一本二区三区精品| 久久99精品国语久久久| 69av精品久久久久久| 午夜福利在线观看吧| 少妇裸体淫交视频免费看高清| 在线免费观看的www视频| 免费观看a级毛片全部| 精品久久久久久久久亚洲| 成年女人永久免费观看视频| 性插视频无遮挡在线免费观看| 午夜爱爱视频在线播放| 免费大片18禁| 欧美日韩在线观看h| 久久99蜜桃精品久久| 国产男人的电影天堂91| 岛国在线免费视频观看| 国产亚洲精品av在线| 亚洲精品乱码久久久久久按摩| 国产精品精品国产色婷婷| 欧美成人精品欧美一级黄| av天堂中文字幕网| 丰满乱子伦码专区| 欧美成人午夜免费资源| 日韩精品青青久久久久久| 免费观看性生交大片5| 日韩欧美国产在线观看| 日本五十路高清| 欧美成人a在线观看| 国产成人a区在线观看| 久久99精品国语久久久| 久久国内精品自在自线图片| 国内精品一区二区在线观看| 国产成人福利小说| 超碰97精品在线观看| 国产精品一区二区在线观看99 | 色噜噜av男人的天堂激情| 久久午夜福利片| 99热精品在线国产| 久久精品久久久久久噜噜老黄 | 国产精品福利在线免费观看| 午夜福利视频1000在线观看| 欧美人与善性xxx| 久久婷婷人人爽人人干人人爱| 一本久久精品| 日韩国内少妇激情av| 国产精品一及| 国产精品伦人一区二区| 永久网站在线| 久久久精品94久久精品| 成年免费大片在线观看| 大话2 男鬼变身卡| 久久久精品94久久精品| 麻豆乱淫一区二区| 欧美三级亚洲精品| 国产成人a∨麻豆精品| 精品欧美国产一区二区三| 国产成人福利小说| 国产精品女同一区二区软件| 国国产精品蜜臀av免费| 少妇的逼好多水| 日韩强制内射视频| 欧美性感艳星| 最近的中文字幕免费完整| 可以在线观看毛片的网站| 黄色配什么色好看| 亚洲欧美精品专区久久| 久久久久久久午夜电影| 99久国产av精品国产电影| 内地一区二区视频在线| 少妇高潮的动态图| 成人三级黄色视频| 我要搜黄色片| 国产亚洲一区二区精品| 欧美高清成人免费视频www| 在线观看av片永久免费下载| 卡戴珊不雅视频在线播放| 99热精品在线国产| 我的女老师完整版在线观看| 晚上一个人看的免费电影| 国产免费一级a男人的天堂| 日本av手机在线免费观看| av线在线观看网站| 日本与韩国留学比较| 亚洲综合精品二区| 久久久久久久久久久免费av| 国产亚洲最大av| 直男gayav资源| 午夜日本视频在线| 免费观看的影片在线观看| 色5月婷婷丁香| 在线免费观看的www视频| 日本黄色片子视频| 精品人妻偷拍中文字幕| 99久久精品国产国产毛片| 亚洲欧美日韩东京热| 最近最新中文字幕免费大全7| 三级国产精品片| 国产熟女欧美一区二区| 亚洲av免费在线观看| 久久99热6这里只有精品| 国产免费福利视频在线观看| 最后的刺客免费高清国语| 欧美日本视频| 九色成人免费人妻av| 午夜精品国产一区二区电影 | 非洲黑人性xxxx精品又粗又长| 色尼玛亚洲综合影院| 中文字幕久久专区| 久久鲁丝午夜福利片| 国产成人91sexporn| 国产精品久久久久久久电影| 欧美bdsm另类| 国产精品爽爽va在线观看网站| 国产黄色视频一区二区在线观看 | 美女大奶头视频| 一边亲一边摸免费视频| 亚洲av免费高清在线观看| 国产精品麻豆人妻色哟哟久久 | 神马国产精品三级电影在线观看| 亚洲av男天堂| 最近2019中文字幕mv第一页| 国产伦精品一区二区三区四那| 欧美日韩在线观看h| 国产人妻一区二区三区在| 亚洲经典国产精华液单| 青青草视频在线视频观看| 午夜精品一区二区三区免费看| 最新中文字幕久久久久| 国产单亲对白刺激| 国产三级在线视频| 乱系列少妇在线播放| 在线观看66精品国产| 国产一区亚洲一区在线观看| 亚洲高清免费不卡视频| 成人毛片a级毛片在线播放| 大话2 男鬼变身卡| АⅤ资源中文在线天堂| 久久草成人影院| 99视频精品全部免费 在线| 国产国拍精品亚洲av在线观看| 一级av片app| 人体艺术视频欧美日本| 欧美性猛交╳xxx乱大交人| 99久久人妻综合| 麻豆久久精品国产亚洲av| 国产黄a三级三级三级人| 欧美成人精品欧美一级黄| 天堂中文最新版在线下载 | 蜜桃久久精品国产亚洲av| 国产精品一区二区性色av| 天美传媒精品一区二区| 久久久精品大字幕| 欧美日本亚洲视频在线播放| 国产在视频线精品| 欧美高清成人免费视频www| 99热网站在线观看| 亚洲va在线va天堂va国产| 久久这里有精品视频免费| 草草在线视频免费看| 中文乱码字字幕精品一区二区三区 | 国产精品伦人一区二区| 全区人妻精品视频| 国产精品精品国产色婷婷| 丝袜喷水一区| 亚洲精品乱码久久久久久按摩| 免费搜索国产男女视频| 欧美97在线视频| 欧美日韩精品成人综合77777| 成人漫画全彩无遮挡| 蜜桃久久精品国产亚洲av| 精品人妻偷拍中文字幕| 午夜亚洲福利在线播放| 日韩欧美国产在线观看| 欧美又色又爽又黄视频| 少妇猛男粗大的猛烈进出视频 | 青春草亚洲视频在线观看| 色尼玛亚洲综合影院| 亚洲精品亚洲一区二区| 日韩,欧美,国产一区二区三区 | 欧美极品一区二区三区四区| 国产高清不卡午夜福利| 成年女人永久免费观看视频| 亚洲美女搞黄在线观看| 亚洲av日韩在线播放| 亚洲av熟女| 久久精品国产99精品国产亚洲性色| 国产免费男女视频| www日本黄色视频网| 日本欧美国产在线视频| 欧美激情国产日韩精品一区| 国产亚洲91精品色在线| 少妇熟女欧美另类| 一个人看视频在线观看www免费| 美女黄网站色视频| 成年免费大片在线观看| 一区二区三区高清视频在线| 一边摸一边抽搐一进一小说| 国产午夜精品论理片| 热99在线观看视频| 特大巨黑吊av在线直播| 久久精品熟女亚洲av麻豆精品 | 亚洲成av人片在线播放无| 亚洲精品久久久久久婷婷小说 | 国产精品蜜桃在线观看| 欧美日本视频| 国内揄拍国产精品人妻在线| 精品久久久久久久末码| 美女被艹到高潮喷水动态| 亚洲精品日韩在线中文字幕| 岛国在线免费视频观看| 日韩欧美国产在线观看| 国产三级中文精品| 一区二区三区乱码不卡18| 亚洲国产日韩欧美精品在线观看| 深爱激情五月婷婷| 日本熟妇午夜| 性插视频无遮挡在线免费观看| 亚洲色图av天堂| av视频在线观看入口| 国内精品一区二区在线观看| 亚洲五月天丁香| 少妇高潮的动态图| 国模一区二区三区四区视频| 免费观看在线日韩| eeuss影院久久| 欧美成人精品欧美一级黄| 亚洲成人中文字幕在线播放| 亚洲av不卡在线观看| 久久国产乱子免费精品| 女的被弄到高潮叫床怎么办| 男女下面进入的视频免费午夜| 久久久久久久国产电影| 亚洲怡红院男人天堂| 国产av一区在线观看免费| 久久这里只有精品中国| 老师上课跳d突然被开到最大视频| 久久久成人免费电影| 午夜久久久久精精品| 日韩欧美三级三区| 欧美97在线视频| 精品久久久久久久久av| 天堂av国产一区二区熟女人妻| 99在线人妻在线中文字幕| 一级毛片aaaaaa免费看小| 成年女人永久免费观看视频| 日韩成人伦理影院| 成人午夜高清在线视频| 一级黄色大片毛片| 久久久色成人| 日本免费一区二区三区高清不卡| 非洲黑人性xxxx精品又粗又长| 日本五十路高清| 欧美3d第一页| 色噜噜av男人的天堂激情| 亚洲图色成人| 亚洲国产欧美在线一区| 免费一级毛片在线播放高清视频| 亚洲最大成人av| 波多野结衣巨乳人妻| a级毛片免费高清观看在线播放| av专区在线播放| 高清av免费在线| 亚洲三级黄色毛片| 亚洲内射少妇av| 99久久精品一区二区三区| 亚洲国产欧美人成| 中文资源天堂在线| 91在线精品国自产拍蜜月| 欧美激情久久久久久爽电影| 97热精品久久久久久| 亚洲电影在线观看av| 久久久国产成人免费| 日韩一本色道免费dvd| 免费黄色在线免费观看| 哪个播放器可以免费观看大片| 蜜桃久久精品国产亚洲av| 99热网站在线观看| 国产精品电影一区二区三区| 美女xxoo啪啪120秒动态图| 久久久久国产网址| 国模一区二区三区四区视频| 我的老师免费观看完整版| 精品人妻熟女av久视频| 2021天堂中文幕一二区在线观| 亚洲欧美精品综合久久99| 久久久久精品久久久久真实原创| 久久人妻av系列| 国产成人福利小说| 国产伦一二天堂av在线观看| 国产精品久久久久久av不卡| 精品久久久久久久久久久久久| 亚洲精品成人久久久久久| 欧美另类亚洲清纯唯美| 色综合亚洲欧美另类图片| 欧美一级a爱片免费观看看| 国产一区二区在线av高清观看| 国产成人午夜福利电影在线观看| 一区二区三区四区激情视频| 夜夜看夜夜爽夜夜摸| 能在线免费观看的黄片| 国产三级在线视频| 直男gayav资源| 亚洲18禁久久av| av专区在线播放| 99九九线精品视频在线观看视频| 白带黄色成豆腐渣| 国产精品嫩草影院av在线观看| 我的老师免费观看完整版| 韩国av在线不卡| 国产一区二区在线观看日韩| 国产一区有黄有色的免费视频 | 日本免费a在线| 女人久久www免费人成看片 | 国产真实伦视频高清在线观看| 成年av动漫网址| 七月丁香在线播放| 免费一级毛片在线播放高清视频| 免费看美女性在线毛片视频| h日本视频在线播放| 亚洲av免费高清在线观看| 久久久午夜欧美精品| 久热久热在线精品观看| 久久久久久久久久黄片| 国产亚洲精品av在线| 色播亚洲综合网| 国产一级毛片七仙女欲春2| 99热全是精品| 国产精品av视频在线免费观看| 久久人妻av系列| 色综合站精品国产| 99久久九九国产精品国产免费| 99久久精品国产国产毛片| 久久久午夜欧美精品| 精品99又大又爽又粗少妇毛片| 黑人高潮一二区| 免费av观看视频| 少妇猛男粗大的猛烈进出视频 | 色网站视频免费| 久99久视频精品免费| 91aial.com中文字幕在线观看| 久久久久久久久久久丰满| 国产男人的电影天堂91| 欧美+日韩+精品| 午夜老司机福利剧场| 91久久精品电影网| 女的被弄到高潮叫床怎么办| 亚洲国产欧洲综合997久久,| 欧美激情国产日韩精品一区| 男插女下体视频免费在线播放| 纵有疾风起免费观看全集完整版 | 观看免费一级毛片| 舔av片在线| 波野结衣二区三区在线| 国产麻豆成人av免费视频| 97人妻精品一区二区三区麻豆| 国产乱人偷精品视频| 亚洲内射少妇av| 两性午夜刺激爽爽歪歪视频在线观看| 最后的刺客免费高清国语| 亚洲精品亚洲一区二区| 久久精品国产亚洲av天美| 色哟哟·www| 伊人久久精品亚洲午夜| 国产精华一区二区三区| 亚洲自拍偷在线| 最近最新中文字幕大全电影3| 日韩强制内射视频| 老司机影院毛片| 人妻少妇偷人精品九色| www日本黄色视频网| 天堂网av新在线| 超碰av人人做人人爽久久| 三级男女做爰猛烈吃奶摸视频| 最近2019中文字幕mv第一页| 亚洲欧美成人综合另类久久久 | 一夜夜www| 国产午夜福利久久久久久| a级毛色黄片| 国产精品一及| 色5月婷婷丁香| 久久久久九九精品影院| av播播在线观看一区| 久久人妻av系列| 久久综合国产亚洲精品| 日本-黄色视频高清免费观看| 啦啦啦啦在线视频资源| 亚洲自偷自拍三级| 高清视频免费观看一区二区 | 国产成人精品一,二区| 观看美女的网站| 亚洲中文字幕一区二区三区有码在线看| a级毛色黄片| 日韩大片免费观看网站 | 亚洲av福利一区| 日本熟妇午夜| 成人av在线播放网站| www日本黄色视频网| 久久精品熟女亚洲av麻豆精品 | 99在线视频只有这里精品首页| 精华霜和精华液先用哪个| 青春草视频在线免费观看| 亚洲欧美日韩卡通动漫| 天天躁夜夜躁狠狠久久av| 高清毛片免费看| 91久久精品国产一区二区三区| 成人亚洲欧美一区二区av| 麻豆精品久久久久久蜜桃| 国产成人91sexporn| 日韩成人伦理影院| 国产亚洲精品久久久com| 美女大奶头视频| 国产淫语在线视频| 欧美zozozo另类| 男女啪啪激烈高潮av片| 国产在视频线在精品| 免费看美女性在线毛片视频| 亚洲欧美精品自产自拍| 成人亚洲精品av一区二区| 国产av不卡久久| 亚洲国产精品合色在线| 亚洲欧美中文字幕日韩二区| 国产精品av视频在线免费观看| 97超碰精品成人国产| 美女cb高潮喷水在线观看| 日本午夜av视频| 国产激情偷乱视频一区二区| 97热精品久久久久久| 中文欧美无线码| 亚洲人成网站高清观看| 欧美激情在线99| 亚洲av中文av极速乱| 久久人妻av系列| 免费大片18禁| 国产精品久久久久久久电影| 成人鲁丝片一二三区免费| 久久久久久久久久黄片| 内射极品少妇av片p| 国产人妻一区二区三区在| 久久久久久九九精品二区国产| 在线播放国产精品三级| 免费观看人在逋| 国产精品综合久久久久久久免费| 在线播放国产精品三级| 少妇裸体淫交视频免费看高清| 国产高潮美女av| 久久精品国产亚洲av天美| 中国国产av一级| 国产午夜福利久久久久久| 久久久久久久国产电影| 我要看日韩黄色一级片| 大话2 男鬼变身卡| 久久精品国产亚洲av天美| 青青草视频在线视频观看| 永久免费av网站大全| 亚洲国产精品久久男人天堂| 色综合站精品国产| 一级黄色大片毛片| 嘟嘟电影网在线观看| 欧美性猛交黑人性爽| 精品国产一区二区三区久久久樱花 | 亚洲人成网站在线观看播放| 春色校园在线视频观看| 哪个播放器可以免费观看大片| 麻豆一二三区av精品| 国产亚洲av嫩草精品影院| 成人国产麻豆网| 高清视频免费观看一区二区 | 久久久欧美国产精品|