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

    基于高維數(shù)據(jù)優(yōu)化聚類的長(zhǎng)周期峰谷時(shí)段劃分模型研究

    2016-09-13 03:31:30張文月王玉瑋
    現(xiàn)代電力 2016年4期
    關(guān)鍵詞:時(shí)點(diǎn)峰谷高維

    李 娜,王 磊,張文月,王玉瑋,舒 艷,張 超

    (1.國網(wǎng)天津市電力公司經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究院,天津 300171;2.國網(wǎng)天津市電力公司,天津 300010; 3.華北電力大學(xué),北京 102206)

    ?

    基于高維數(shù)據(jù)優(yōu)化聚類的長(zhǎng)周期峰谷時(shí)段劃分模型研究

    李娜1,王磊1,張文月2,王玉瑋3,舒艷3,張超3

    (1.國網(wǎng)天津市電力公司經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究院,天津300171;2.國網(wǎng)天津市電力公司,天津300010; 3.華北電力大學(xué),北京102206)

    0 引 言

    峰谷分時(shí)電價(jià)是需求側(cè)管理中一種有效的經(jīng)濟(jì)手段,因此合理科學(xué)地引入需求側(cè)管理來設(shè)計(jì)峰谷分時(shí)電價(jià),通過運(yùn)用價(jià)格信號(hào)有效調(diào)整電力用戶的用電習(xí)慣和用電計(jì)劃,對(duì)提高資源利用率、優(yōu)化資源配置、實(shí)現(xiàn)削峰填谷、降低發(fā)電成本、促進(jìn)電源與電網(wǎng)的協(xié)調(diào)發(fā)展以及推動(dòng)整個(gè)電力行業(yè)健康發(fā)展起著決定性的作用。

    峰谷分時(shí)電價(jià)制定的恰當(dāng)與否跟峰谷時(shí)段的劃分和用戶需求響應(yīng)密切相關(guān)[1]。峰谷時(shí)段的劃分是峰谷分時(shí)電價(jià)的定價(jià)基礎(chǔ),其劃分方法的選取直接影響用戶的需求響應(yīng)程度,從而影響電價(jià)的實(shí)施效果。因此,制定合理的峰谷分時(shí)電價(jià)的首要任務(wù)是科學(xué)合理地進(jìn)行峰谷時(shí)段的劃分。

    目前已有的關(guān)于峰谷時(shí)段劃分的方法主要有:文獻(xiàn)[1,2]以日負(fù)荷曲線分布為基礎(chǔ),結(jié)合模糊聚類分析技術(shù),構(gòu)建了峰谷時(shí)段劃分模型;文獻(xiàn)[3]基于密度聚類算法,對(duì)年持續(xù)負(fù)荷曲線中各負(fù)荷所對(duì)應(yīng)的持續(xù)時(shí)間進(jìn)行聚類分析,然后通過分布集邊界所對(duì)應(yīng)的負(fù)荷大小劃分峰谷時(shí)段;文獻(xiàn)[4]采用因素分析法,將其與原始負(fù)荷電量相結(jié)合,構(gòu)建峰谷時(shí)段劃分模型;文獻(xiàn)[5]提出了一種基于不同時(shí)段供電成本差異的時(shí)段劃分模型。

    在上述研究成果的基礎(chǔ)上,本文著重考慮峰谷時(shí)段劃分結(jié)果的客觀合理性及其是否能在一個(gè)較長(zhǎng)時(shí)間周期(如1a)內(nèi)適用這兩個(gè)問題。為此,本文提出以下研究思路:首先,通過數(shù)據(jù)高維化的處理方法構(gòu)建涵蓋較長(zhǎng)時(shí)間周期(例如1a)內(nèi)所有負(fù)荷信息的數(shù)據(jù)樣本集;其次,以K-均值算法為聚類分析工具,在高維數(shù)據(jù)樣本集上構(gòu)建峰谷時(shí)段劃分模型;最后,結(jié)合某區(qū)全年負(fù)荷數(shù)據(jù),對(duì)所構(gòu)建的模型進(jìn)行算例仿真,在驗(yàn)證模型的合理性基礎(chǔ)上,最終輸出時(shí)段劃分結(jié)果。

    1 峰谷時(shí)段劃分的理論依據(jù)

    對(duì)全天進(jìn)行峰、平、谷時(shí)段劃分,其基本原理就是根據(jù)各時(shí)點(diǎn)(通常做法是將全天分為24個(gè)時(shí)點(diǎn)或時(shí)段)上所采集的負(fù)荷數(shù)值大小,將負(fù)荷取值相對(duì)大的各時(shí)點(diǎn)統(tǒng)一劃分為全天的峰時(shí)段、將負(fù)荷取值相對(duì)小的各時(shí)點(diǎn)統(tǒng)一劃分為全天的谷時(shí)段、將負(fù)荷取值相對(duì)中等的各時(shí)點(diǎn)統(tǒng)一劃分為全天的平時(shí)段。由此可見,基于峰、平、谷時(shí)段劃分原理的各類算法本質(zhì)上屬于數(shù)據(jù)挖掘中聚類分析(無導(dǎo)師學(xué)習(xí))的范疇,劃分效果的好壞主要取決于以下兩點(diǎn):

    ① 包含基本信息的數(shù)據(jù)樣本集合如何構(gòu)造,使對(duì)其進(jìn)行時(shí)段劃分的最終結(jié)果適用于一個(gè)較長(zhǎng)的時(shí)間周期(一季、一年或者更長(zhǎng));

    ② 判斷各樣本點(diǎn)取值(相當(dāng)于各時(shí)點(diǎn)負(fù)荷值)相對(duì)大小的分類標(biāo)準(zhǔn)如何確定,使據(jù)此進(jìn)行時(shí)段劃分的結(jié)果最大限度地反映出峰、平、谷各時(shí)段在負(fù)荷取值方面的差異。

    對(duì)以上兩點(diǎn)問題的合理解決將是本文論述的關(guān)鍵。數(shù)據(jù)樣本點(diǎn)的高維化和聚類算法可通過迭代自動(dòng)收斂于最優(yōu)分類標(biāo)準(zhǔn)是這些目標(biāo)實(shí)現(xiàn)的基礎(chǔ)。接下來分別就數(shù)據(jù)樣本點(diǎn)的高維化和K-均值聚類算法展開研究。

    2 基于高維度量空間的數(shù)據(jù)樣本集構(gòu)建

    為了使最終的時(shí)段劃分結(jié)果適用于一個(gè)較長(zhǎng)的時(shí)間周期,僅使用某天的各時(shí)點(diǎn)負(fù)荷數(shù)據(jù)構(gòu)造數(shù)據(jù)樣本集合顯然是不合理的。這是因?yàn)椋谝粋€(gè)較長(zhǎng)的時(shí)間周期內(nèi)(例如1a),各日負(fù)荷曲線的形狀差異會(huì)比較明顯,且這種差異程度會(huì)隨著周期加長(zhǎng)而增大,從而適用于某日的峰、平、谷時(shí)段劃分結(jié)果并不一定適用于另一日。面對(duì)這一問題,比較直觀的處理方法是:

    ① 在一個(gè)較長(zhǎng)的時(shí)間周期內(nèi)(以下稱其為全周期),分別以每一日的時(shí)點(diǎn)負(fù)荷數(shù)據(jù)為樣本集,采用一定的聚類方劃分出各日的峰、平、谷時(shí)段;

    ② 根據(jù)各日的時(shí)段劃分結(jié)果統(tǒng)計(jì)出24個(gè)時(shí)點(diǎn)分別被劃分為峰、平、谷各時(shí)段的天數(shù);

    ③ 根據(jù)最大天數(shù)原則,判斷某個(gè)時(shí)點(diǎn)最終應(yīng)當(dāng)被劃分于峰、平、谷某一時(shí)段,從而得到了考慮全周期所有日負(fù)荷時(shí)間分布特點(diǎn)的峰、平、谷時(shí)段劃分結(jié)果。

    上述方法在所選全周期長(zhǎng)度較小時(shí)比較可行,但是當(dāng)周期長(zhǎng)度(天數(shù))增加時(shí),該方法的迭代次數(shù)會(huì)隨之成倍增長(zhǎng),且一旦出現(xiàn)某一時(shí)點(diǎn)被劃分于峰、平、谷各時(shí)段天數(shù)相等的情況,則最大天數(shù)原則失效,算法無法輸出最終的劃分結(jié)果。

    為此,本文建立如下考慮全周期各日負(fù)荷信息的數(shù)據(jù)樣本集構(gòu)造方法。

    設(shè)所分析全周期內(nèi)共含天數(shù)為n(n∈N,n>0),在第t(t=1,2,…,24)個(gè)時(shí)點(diǎn)上,所研究地區(qū)n天中對(duì)應(yīng)的負(fù)荷數(shù)據(jù)可組成如下向量:

    (1)

    式中:xt,i表示第i(i∈N,0

    定義如下集合:

    (2)

    則向量S為集合Rn中的可數(shù)子集。Rn為

    向量S即為所構(gòu)造的用于時(shí)段劃分的數(shù)據(jù)樣本集合,該集合包含了全周期內(nèi)各個(gè)時(shí)點(diǎn)t上所有日的負(fù)荷數(shù)據(jù)。

    但是,上面定義的樣本集S并不能直接用于時(shí)段劃分的計(jì)算,其原因是集合中沒有表示元素之間“遠(yuǎn)近”的概念,因此無法進(jìn)行分類,需要做如下處理:

    不影響最終分析結(jié)果,在集合Rn中定義如下歐式距離,從而Rn構(gòu)成n維實(shí)歐式空間:

    (3)

    式中:d(x,y)∈R表示空間Rn中任意兩點(diǎn)x,y之間的歐式距離。因S為Rn的子集,因此,S中任意兩個(gè)元素(樣本)之間的距離可表示為

    (4)

    式中:t1,t2=1,2,…,24。

    定義距離概念之后,就可以根據(jù)距離的大小來判斷空間Rn中任意兩點(diǎn)的“遠(yuǎn)近”程度。顯然,這種“遠(yuǎn)近”的概念同樣適用于S中的元素。

    以上方法所構(gòu)造的數(shù)據(jù)樣本集S具有以下兩項(xiàng)優(yōu)點(diǎn):

    ① S中包含了全周期各時(shí)點(diǎn)上所有負(fù)荷數(shù)據(jù),從而在S上進(jìn)行合理的聚類分析,可得到適用于全周期的峰、平、谷時(shí)段劃分結(jié)果;

    ② S中的元素,即樣本點(diǎn)是以n維向量的形式給出的,這種將數(shù)據(jù)高維化的構(gòu)造方式避免了算法迭代次數(shù)因全周期天數(shù)增加而成倍增長(zhǎng)的不足,可適用于較長(zhǎng)時(shí)間跨度(數(shù)年)全周期上峰、平、谷時(shí)段的劃分。

    3 時(shí)段劃分模型構(gòu)建

    本文依據(jù)所建樣本集S中的數(shù)據(jù)(樣本點(diǎn))具有高維度、數(shù)值型、無標(biāo)記等特點(diǎn),選取K-均值算法(以下用K-means表示)作為分析工具,對(duì)樣本集S進(jìn)行峰、平、谷時(shí)段劃分。

    K-means是最常用的聚類算法之一,能有效地處理規(guī)模較大和高維的數(shù)據(jù)集合,能對(duì)大型數(shù)據(jù)集進(jìn)行高效分類。該算法按照定義的劃分標(biāo)準(zhǔn),把數(shù)據(jù)分成幾組,使得同組內(nèi)數(shù)據(jù)相對(duì)聚攏,相異組之間數(shù)據(jù)相對(duì)疏散,按這樣的方法所形成的每一組就稱作一個(gè)聚類(或聚簇)。在K-means中,劃分準(zhǔn)則常常采用誤差平方和準(zhǔn)則函數(shù)來描述。此算法是通過迭代尋找K個(gè)聚類的一種劃分方案,使得用這K個(gè)均值來代表相應(yīng)各類樣本時(shí)所得到的總體誤差最小。

    對(duì)于不同的聚類方案,誤差平方和準(zhǔn)則函數(shù)值J一般是不相同的。使J值最小的聚類方案即為在誤差平方和準(zhǔn)則下的最優(yōu)聚類結(jié)果。

    誤差平方和J無法用解析的方法最小化,只能采用迭代的方法,通過不斷的調(diào)整樣本的類歸屬情況來獲得最優(yōu)解。

    K-means的效率較高,其缺點(diǎn)是只能處理數(shù)值型數(shù)據(jù)。然而,峰谷時(shí)段劃分問題實(shí)質(zhì)上就是根據(jù)各時(shí)點(diǎn)上負(fù)荷數(shù)值的大小進(jìn)行聚類劃分,因此,該問題正好屬于該算法的適用范疇。

    運(yùn)用K-means進(jìn)行峰谷時(shí)段劃分的原理及基本步驟如下:

    如上節(jié)所述,所構(gòu)造的數(shù)據(jù)樣本集S中每一個(gè)樣本點(diǎn)xt即為一個(gè)n維向量(n為所選全周期中的總天數(shù)),S的樣本容量為24(即24個(gè)時(shí)點(diǎn)),需要將這24個(gè)樣本點(diǎn)重新劃分為峰時(shí)段、平時(shí)段、谷時(shí)段3類。因此,該算法運(yùn)用到本文所構(gòu)建的樣本集S上進(jìn)行峰、平、谷時(shí)段劃分的問題即可描述為:令K=3,即將24個(gè)樣本點(diǎn)(數(shù)據(jù)對(duì)象)劃分為3個(gè)聚類以便使得所獲得的聚類結(jié)果滿足同一聚類中的對(duì)象相似度較高,而不同聚類中的對(duì)象相似度較小。聚類相似度是利用各聚類中所有樣本點(diǎn)的均值所獲得一個(gè)“聚類中心”來進(jìn)行計(jì)算的。

    第一步:在樣本集S中隨機(jī)選K=3個(gè)初始聚類中心X1,(1),X2,(1),X3,(1),其中括號(hào)內(nèi)的序號(hào)為尋找聚類中心的迭代運(yùn)算次序號(hào)。聚類中心的向量值可任意設(shè)定,例如可選開頭的3個(gè)模式樣本集中的向量值作為初始聚類中心。

    第二步:逐個(gè)將S中的樣本點(diǎn)xt按最小歐式距離準(zhǔn)則分配給3個(gè)初始聚類中心中的某一個(gè)Xj,(1)(j=1,2,3)。對(duì)所有的xt,如果在第k次迭代中,xt到Xj,(k)的歐式距離比到其他任何中心點(diǎn)都近,則xt∈Sj,(k)。其中,Sj,(k)表示經(jīng)第k次迭代后所形成的第j(j=1,2,3)個(gè)聚類。其聚類中心為Xj,(k)。

    第三步:在每一次形成新的聚類結(jié)果S1,(k)、S2,(k)、S3,(k)后,進(jìn)一步按類內(nèi)均值線性變換的方法計(jì)算各個(gè)聚類中心的新的向量值Xj,(k+1),類內(nèi)均值線性變換的具體公式為

    (5)

    式中:lj,(k)為Sj,(k)中所包含的樣本點(diǎn)個(gè)數(shù)。并在新的聚類中心基礎(chǔ)上,重新按照最小歐式距離準(zhǔn)則形成新的聚類結(jié)果S1,(k+1)、S2,(k+1)、S3,(k+1),即完成第k+1次迭代。

    第四步:在每一次迭代完成后以均值向量作為新的聚類中心,可構(gòu)造如下誤差平方和準(zhǔn)則函數(shù):

    (6)

    迭代的目的是為了使J值達(dá)到極小,當(dāng)?shù)趉+1次迭代后的J(k+1)值與第k次迭代后的J(k)值相等,則迭代停止,輸出的S1,(k+1)、S2,(k+1)、S3,(k+1)即為使得同一聚類中的對(duì)象相似度最高,而不同聚類中的對(duì)象相似度最小的最優(yōu)聚類結(jié)果,各類按所含元素負(fù)荷值大小區(qū)分為對(duì)應(yīng)的峰時(shí)段聚類、平時(shí)段聚類以及谷時(shí)段聚類。否則,應(yīng)返回第二步開始新一輪的迭代計(jì)算。

    4 算例仿真

    以某地區(qū)負(fù)荷數(shù)據(jù)為例,利用以上所構(gòu)建的模型,對(duì)該區(qū)峰谷時(shí)段劃分進(jìn)行實(shí)例分析。

    首先,為了驗(yàn)證所構(gòu)建模型的合理性與可行性,在全年數(shù)據(jù)中選擇日期相隔較遠(yuǎn)的兩天數(shù)據(jù),各自單獨(dú)進(jìn)行仿真分析。圖1、圖2分別是該區(qū)2013年7月1日和9月1日基于K-means單日時(shí)段劃分模型的仿真結(jié)果,分析方法可做如下簡(jiǎn)要介紹:

    圖1 7月1日峰谷時(shí)段劃分仿真圖像

    圖2 9月1日峰谷時(shí)段劃分仿真圖像

    以對(duì)7月1日各時(shí)點(diǎn)(24個(gè)時(shí)點(diǎn))負(fù)荷數(shù)據(jù)進(jìn)行峰、平、谷時(shí)段劃分為例。第一步,將該日數(shù)據(jù)構(gòu)造成維度為1(1d),樣本點(diǎn)數(shù)量為24個(gè)(24個(gè)時(shí)點(diǎn),每個(gè)樣本點(diǎn)均按所在時(shí)點(diǎn)進(jìn)行編號(hào))的數(shù)據(jù)樣本集S;第二步,在數(shù)據(jù)樣本集上運(yùn)用第3節(jié)中所構(gòu)建的K-均值數(shù)據(jù)挖掘模型進(jìn)行聚類迭代運(yùn)算(將K值設(shè)為3),并最終將S中所包含的所有樣本點(diǎn)劃分為S峰、S平、S谷3個(gè)聚類,其中,相對(duì)應(yīng),S峰類中所包含的各樣本點(diǎn)的時(shí)點(diǎn)編號(hào)組成峰時(shí)段,S平類中所包含的各樣本點(diǎn)的時(shí)點(diǎn)編號(hào)組成平時(shí)段,S谷類中所包含的各樣本點(diǎn)的時(shí)點(diǎn)編號(hào)組成谷時(shí)段。

    從對(duì)以上2日各自的分析結(jié)果中可以看出,7月1日時(shí)段劃分結(jié)果(圖中階梯線表征,其中最低梯級(jí)線表示谷時(shí)段的位置與長(zhǎng)度,中間梯級(jí)表示平時(shí)段的位置與長(zhǎng)度,最高梯級(jí)表示峰時(shí)段的位置與長(zhǎng)度)明顯與9月1日存在差異(主要差異體現(xiàn)在下午14-15點(diǎn)的劃分類別不同),這主要是由于這兩日各自具體的負(fù)荷曲線形狀不同所引起的。

    類似地,若希望將一個(gè)日峰、平、谷時(shí)段劃分結(jié)果夠推廣到至少1a的時(shí)間周期內(nèi)長(zhǎng)期使用,則得到這個(gè)劃分結(jié)果的過程中必須同時(shí)考慮到在1a的時(shí)間周期內(nèi),各日整點(diǎn)負(fù)荷數(shù)值的大小分布情況,這顯然不是簡(jiǎn)單進(jìn)行某個(gè)單日數(shù)據(jù)的分析所能得到的。上述仿真分析驗(yàn)證了本文第2節(jié)一開始便提出的結(jié)論:為了使最終的時(shí)段劃分結(jié)果適用于一個(gè)較長(zhǎng)的時(shí)間周期,僅使用某天的各時(shí)點(diǎn)負(fù)荷數(shù)據(jù)構(gòu)造數(shù)據(jù)樣本集合顯然是不合理的。

    其次,為了解決上述問題,得到一組可以在1a的時(shí)間周期內(nèi)推廣的日峰、平、谷時(shí)段劃分結(jié)果,本文以該地區(qū)2013年全年的負(fù)荷數(shù)據(jù)進(jìn)行仿真,如圖3是全年峰谷時(shí)段劃分的仿真圖像。分析方法可做如下簡(jiǎn)要介紹:

    第一步,將全年整點(diǎn)負(fù)荷數(shù)據(jù)構(gòu)造成維度為365(365d),樣本點(diǎn)數(shù)量為24個(gè)(24個(gè)時(shí)點(diǎn),每個(gè)樣本點(diǎn)均按所在時(shí)點(diǎn)進(jìn)行編號(hào))的數(shù)據(jù)樣本集S;第二步,在數(shù)據(jù)樣本集上運(yùn)用第3節(jié)中所構(gòu)建的K-均值數(shù)據(jù)挖掘模型進(jìn)行聚類迭代運(yùn)算(將K值設(shè)為3),并最終將S中所包含的所有樣本點(diǎn)劃分為S峰、S平、S谷3個(gè)聚類,其中,相對(duì)應(yīng),S峰類中所包含的各樣本點(diǎn)的時(shí)點(diǎn)編號(hào)組成峰時(shí)段,S平類中所包含的各樣本點(diǎn)的時(shí)點(diǎn)編號(hào)組成平時(shí)段,S谷類中所包含的各樣本點(diǎn)的時(shí)點(diǎn)編號(hào)組成谷時(shí)段。

    圖3中,階梯線依然表征峰、平、谷各時(shí)段的位置及長(zhǎng)度,曲線簇代表該地區(qū)2013年全年各日(365d)按整點(diǎn)數(shù)據(jù)所繪制的負(fù)荷曲線。本文模型所得出的峰谷時(shí)段劃分結(jié)果如下:

    峰時(shí)段:10:00—14:00,15:00—21:00;平時(shí)段:8:00—10:00,14:00—15:00,21:00—23:00;谷時(shí)段:0:00—8:00。

    圖3 全年峰谷時(shí)段劃分仿真圖像

    該區(qū)2013年的時(shí)段劃分情況如下:

    峰時(shí)段:8:00—11:00,16:00—21:00;平時(shí)段:6:00—8:00,11:00—16:00,21:00—22:00;谷時(shí)段:22:00—6:00。

    將以上兩種時(shí)段的劃分結(jié)果進(jìn)行對(duì)比分析可知:峰谷時(shí)段的劃分是以“三平兩峰一谷”的形式存在,仿真結(jié)果比實(shí)際的時(shí)段劃分更加合理。可舉例說明:從圖3中全年各日負(fù)荷曲線分布上可以看出(圖中曲線簇),在中午12時(shí)左右的負(fù)荷水平基本處于全天的最高峰(不符合這一規(guī)律的日負(fù)荷曲線在該地區(qū)2013年全年中未找到),即便從直觀上分析,亦理應(yīng)將該時(shí)點(diǎn)劃分于全天的峰時(shí)段才更合適。

    最后,通過算例仿真結(jié)果可知,本文所建模型的輸出結(jié)果(即峰、平、谷時(shí)段劃分結(jié)果),既可在一個(gè)較長(zhǎng)的時(shí)間周期內(nèi)適用,同時(shí)又客觀地反映出所分各時(shí)段之間的負(fù)荷差異,其合理性與可行性在一定程度上得到了驗(yàn)證。

    5 結(jié) 論

    本文通過研究可得以下結(jié)論:

    ① 通過數(shù)據(jù)高維化處理所構(gòu)造的峰谷時(shí)段劃分?jǐn)?shù)據(jù)樣本集S具有以下兩項(xiàng)優(yōu)點(diǎn):

    a.其包含了一個(gè)較長(zhǎng)時(shí)間周期(如1a)內(nèi)24個(gè)時(shí)點(diǎn)上的所有負(fù)荷信息,從而為得到長(zhǎng)期適用的峰谷時(shí)段劃分結(jié)果提供了完整的數(shù)據(jù)依據(jù)。

    b.樣本集S將全周期中相同時(shí)點(diǎn)上浩繁的負(fù)荷數(shù)據(jù)整合為同一向量(樣本點(diǎn))上不同的坐標(biāo)分量,這樣一來大大簡(jiǎn)化了樣本容量(無論所研究周期有多長(zhǎng),S的樣本容量總是為24個(gè)樣本點(diǎn)),從而降低了模型的迭代次數(shù)并保證了其一致收斂性。

    ② 引入K-means進(jìn)行建模克服了處理高維度、數(shù)值型樣本聚類問題難以收斂的困難,并且算法的迭代過程本身就是一個(gè)時(shí)段劃分的優(yōu)化過程,上述兩點(diǎn)可保證得到同一聚類中的對(duì)象相似度最高,而不同聚類中的對(duì)象相似度最小的峰、平、谷時(shí)段劃分優(yōu)化結(jié)果。

    ③ 通過算例仿真,除了驗(yàn)證模型的合理性以及得到針對(duì)某地區(qū)全年負(fù)荷數(shù)據(jù)的最優(yōu)時(shí)段劃分結(jié)果外,還可以看出,針對(duì)某一日負(fù)荷數(shù)據(jù)得到的時(shí)段劃分結(jié)果與針對(duì)全年數(shù)據(jù)得到的時(shí)段劃分結(jié)果不同,這一點(diǎn)從實(shí)驗(yàn)的角度驗(yàn)證了本文所提出的觀點(diǎn):為了使峰谷時(shí)段劃分結(jié)果能夠客觀反映出各時(shí)段的負(fù)荷差異,且能夠在一個(gè)較長(zhǎng)的時(shí)間周期(例如月、季、半年、年等)內(nèi)適用,因此在高維數(shù)據(jù)樣本集的基礎(chǔ)上引入K-means來構(gòu)建適用于一定周期內(nèi)的時(shí)段劃分模型。

    [1]連振洲,溫步瀛,江岳文.基于負(fù)荷曲線分布特征的峰谷時(shí)段劃分和修正策略研究[J].電網(wǎng)與清潔能源,2014(7):15-19.

    [2]丁寧.基于DSM的峰谷時(shí)段劃分及分時(shí)電價(jià)的研究[D].南京:南京理工大學(xué),2002.

    [3]喬慧婷.基于密度聚類的峰谷時(shí)段劃分方法研究[D].北京:華北電力大學(xué),2011.

    [4]劉菁菁.基于因素分析法的峰谷時(shí)段劃分的研究[D].天津:河北工業(yè)大學(xué),2006.

    [5]林旻,朱艷卉,胡百林.基于供電成本的峰谷時(shí)段劃分及分時(shí)電價(jià)研究[J].華東電力,2005(12):90-91.

    [6]趙娟,譚忠富,李強(qiáng).我國峰谷分時(shí)電價(jià)的狀況分析[J].現(xiàn)代電力,2005,22(2):82-85.

    [7]戴詩容,雷霞,程道衛(wèi)等.電動(dòng)汽車峰谷分時(shí)充放電電價(jià)研究[J].電網(wǎng)與清潔能源,2013,29(7):77-82.

    [8]翟娜娜.基于用戶需求響應(yīng)的峰谷時(shí)段劃分研究[D].北京:華北電力大學(xué),2011.

    [9]李泓澤,李婷婷.峰谷分時(shí)電價(jià)對(duì)供需雙方的影響分析[J].現(xiàn)代電力,2008(4):87-92.

    [10]劉嚴(yán),譚忠富,乞建勛.峰谷分時(shí)電價(jià)設(shè)計(jì)的優(yōu)化模型[J].中國管理科學(xué),2005,13(5):87-92.

    [11]周東華.數(shù)據(jù)挖掘中聚類分析的研究與應(yīng)用[D].天津:天津大學(xué),2006.

    (責(zé)任編輯:林海文)

    Research on the Partition Model of Long Period Peak and Valley Time Based on High Dimensional Data Clustering

    LI Na1,WANG Lei1,ZHANG Wenyue2,WANG Yuwei3,SHU Yan3,ZHANG Chao3

    (1.State Grid Tianjin Economic Research Institute,Tianjin 300171,China;2.State Grid Tianjin Electric Power Company, Tianjin 300010,China;3.North China Electric Power University,Beijing 102206,China)

    為了使峰谷時(shí)段劃分結(jié)果客觀反映出各時(shí)段的負(fù)荷差異,且能夠在一個(gè)較長(zhǎng)的時(shí)間周期(例如1a)內(nèi)適用,本文提出一種以數(shù)據(jù)樣本集高維化處理和K-均值聚類分析相結(jié)合的時(shí)段劃分模型。首先,通過數(shù)據(jù)高維化的處理方法構(gòu)建涵蓋較長(zhǎng)時(shí)間周期(例如1a)內(nèi)所有負(fù)荷信息的數(shù)據(jù)樣本集;其次,以K-均值算法為聚類分析工具,在高維數(shù)據(jù)樣本集上構(gòu)建峰谷時(shí)段劃分模型。最后,結(jié)合某區(qū)全年負(fù)荷數(shù)據(jù),對(duì)所構(gòu)建的模型進(jìn)行算例仿真,在驗(yàn)證模型的合理性基礎(chǔ)上,最終輸出時(shí)段劃分結(jié)果。

    時(shí)段劃分;聚類分析;K-均值算法;數(shù)據(jù)高維化

    In this paper,in order to make the results of peak and valley time division reflect the load difference of each period objectively and be applicable in a long period of time (e.g.,1 year),a time division model is presented by combining processing of high-dimension data sampling set and K-means clustering analysis.First of all,a data sample set covering all load information within a long period of time (e.g.,1 year) is built by using the processing method of high-dimension data.Secondly,the peak and valley time division model based on the high-dimension data sample set is built by using K- means clustering analysis.Finally,the numerical simulation of proposed model is carried out by combining the load data of the whole year in certain district,and the final division result of the peak and valley time can be output on the basis of verifying the rationality of the model.

    peak-valley-time division; clustering analysis; K-means clustering algorithm; high-dimensional data

    1007-2322(2016)04-0067-05

    A

    TM715

    國網(wǎng)天津市電力公司科技項(xiàng)目(KJ15-1-28)

    2015-08-26

    李娜(1985—),女,博士,工程師,研究方向?yàn)殡娏κ袌?chǎng)理論與應(yīng)用技術(shù)研究,E-mail:linajiamie@163.com;

    王磊(1981—),女,博士,工程師,研究方向?yàn)殡娏κ袌?chǎng)理論與應(yīng)用技術(shù)研究,E-mail:leiwendy@126.com;

    張文月(1982—),女,工程師,研究方向?yàn)殡娏κ袌?chǎng)理論與應(yīng)用技術(shù)研究,E-mail:zhangwenyue2005@126.com。

    猜你喜歡
    時(shí)點(diǎn)峰谷高維
    湖南省峰谷分時(shí)電價(jià)分析及優(yōu)化策略研究
    湖南電力(2021年1期)2021-04-13 01:37:02
    O2O模式下我國網(wǎng)約車企業(yè)收入確認(rèn)的探析
    淺談峰谷時(shí)段的劃分方法
    一種改進(jìn)的GP-CLIQUE自適應(yīng)高維子空間聚類算法
    論“前后”的時(shí)間義
    基于加權(quán)自學(xué)習(xí)散列的高維數(shù)據(jù)最近鄰查詢算法
    建造中船舶特定化屬性分析
    資金時(shí)間價(jià)值計(jì)算方法再探
    電動(dòng)汽車峰谷分時(shí)電價(jià)時(shí)段充電優(yōu)化模型
    一般非齊次非線性擴(kuò)散方程的等價(jià)變換和高維不變子空間
    动漫黄色视频在线观看| 视频区图区小说| 国产高清videossex| 亚洲国产欧美网| 久久久精品区二区三区| 免费在线观看完整版高清| 欧美在线一区亚洲| 成熟少妇高潮喷水视频| 日韩免费高清中文字幕av| 亚洲全国av大片| 欧美激情 高清一区二区三区| 久久99一区二区三区| av片东京热男人的天堂| 午夜日韩欧美国产| 极品少妇高潮喷水抽搐| 乱人伦中国视频| 19禁男女啪啪无遮挡网站| 久久性视频一级片| 亚洲精品久久午夜乱码| 亚洲国产欧美日韩在线播放| 一本一本久久a久久精品综合妖精| av超薄肉色丝袜交足视频| 大码成人一级视频| 欧美精品啪啪一区二区三区| 天天影视国产精品| 久久99一区二区三区| 中出人妻视频一区二区| 欧美性长视频在线观看| 日韩 欧美 亚洲 中文字幕| 久久人妻福利社区极品人妻图片| 亚洲成a人片在线一区二区| 在线观看免费高清a一片| 国产成人欧美| 五月开心婷婷网| 欧美色视频一区免费| 91大片在线观看| 亚洲一区中文字幕在线| 久久亚洲精品不卡| 黄片大片在线免费观看| 亚洲精品久久午夜乱码| 亚洲 国产 在线| 老汉色av国产亚洲站长工具| 国产欧美日韩综合在线一区二区| 亚洲午夜理论影院| 丁香欧美五月| 电影成人av| 欧美日韩亚洲国产一区二区在线观看 | 久久午夜综合久久蜜桃| 午夜两性在线视频| 99国产精品免费福利视频| 久久精品亚洲av国产电影网| 午夜免费鲁丝| 国产深夜福利视频在线观看| 亚洲性夜色夜夜综合| 国产成人av教育| 日日摸夜夜添夜夜添小说| 久久婷婷成人综合色麻豆| 99re6热这里在线精品视频| 国产成人影院久久av| 久久久久精品国产欧美久久久| 九色亚洲精品在线播放| 一本一本久久a久久精品综合妖精| 青草久久国产| 男男h啪啪无遮挡| 天天躁夜夜躁狠狠躁躁| 国产午夜精品久久久久久| av福利片在线| 亚洲成人手机| 亚洲少妇的诱惑av| 一级片免费观看大全| 国产男女内射视频| 变态另类成人亚洲欧美熟女 | 新久久久久国产一级毛片| 国产精品一区二区在线观看99| 99久久国产精品久久久| 欧美乱色亚洲激情| 午夜免费观看网址| 久久久国产成人免费| 啦啦啦在线免费观看视频4| 久久久久久人人人人人| 叶爱在线成人免费视频播放| 色综合欧美亚洲国产小说| 在线十欧美十亚洲十日本专区| 免费少妇av软件| 亚洲av熟女| 国产亚洲欧美精品永久| 在线观看免费视频日本深夜| 深夜精品福利| x7x7x7水蜜桃| 亚洲一区中文字幕在线| 久久精品国产清高在天天线| 亚洲人成电影免费在线| 亚洲专区字幕在线| 日韩 欧美 亚洲 中文字幕| 亚洲国产欧美一区二区综合| 黄片播放在线免费| 国产欧美日韩一区二区三区在线| 国产精品一区二区在线观看99| 国产97色在线日韩免费| 97人妻天天添夜夜摸| 国产蜜桃级精品一区二区三区 | 久久久久久久国产电影| 中文字幕人妻丝袜制服| 伦理电影免费视频| 纯流量卡能插随身wifi吗| 午夜福利在线免费观看网站| 欧美日韩亚洲国产一区二区在线观看 | 巨乳人妻的诱惑在线观看| 欧美激情高清一区二区三区| 91大片在线观看| 每晚都被弄得嗷嗷叫到高潮| 成人免费观看视频高清| 美女 人体艺术 gogo| 中文字幕高清在线视频| videos熟女内射| 丝袜美腿诱惑在线| 午夜久久久在线观看| 在线观看66精品国产| 国产精品免费一区二区三区在线 | 国产精品1区2区在线观看. | 国产高清国产精品国产三级| 大型黄色视频在线免费观看| 丝袜美腿诱惑在线| 黑人猛操日本美女一级片| 欧美激情久久久久久爽电影 | 亚洲第一av免费看| 亚洲五月婷婷丁香| 王馨瑶露胸无遮挡在线观看| 国产真人三级小视频在线观看| 最近最新免费中文字幕在线| 久久热在线av| 十八禁高潮呻吟视频| 18禁美女被吸乳视频| 久久久国产一区二区| 在线观看免费视频日本深夜| 成人亚洲精品一区在线观看| 亚洲国产中文字幕在线视频| 黄色 视频免费看| 美女高潮到喷水免费观看| 精品久久久久久电影网| 下体分泌物呈黄色| 老汉色av国产亚洲站长工具| 亚洲精品国产一区二区精华液| 欧美不卡视频在线免费观看 | 两性夫妻黄色片| 亚洲国产精品合色在线| av在线播放免费不卡| 99riav亚洲国产免费| 亚洲在线自拍视频| 免费看a级黄色片| cao死你这个sao货| 高潮久久久久久久久久久不卡| 三上悠亚av全集在线观看| 亚洲人成77777在线视频| 99久久国产精品久久久| av免费在线观看网站| 日本精品一区二区三区蜜桃| 性少妇av在线| 亚洲欧美日韩另类电影网站| 日韩熟女老妇一区二区性免费视频| 色婷婷av一区二区三区视频| 午夜两性在线视频| 国产有黄有色有爽视频| 中文字幕最新亚洲高清| 老熟女久久久| 最近最新中文字幕大全免费视频| 人妻丰满熟妇av一区二区三区 | 国产精品久久视频播放| 男男h啪啪无遮挡| 国产一区有黄有色的免费视频| 亚洲aⅴ乱码一区二区在线播放 | 自线自在国产av| 男女之事视频高清在线观看| 久久热在线av| 国产成人系列免费观看| 国产精华一区二区三区| 搡老熟女国产l中国老女人| 国产真人三级小视频在线观看| 最新的欧美精品一区二区| 日日夜夜操网爽| videos熟女内射| a级毛片在线看网站| xxxhd国产人妻xxx| 成年动漫av网址| 大香蕉久久成人网| 老司机深夜福利视频在线观看| 他把我摸到了高潮在线观看| 极品人妻少妇av视频| 色婷婷久久久亚洲欧美| 一本一本久久a久久精品综合妖精| 久久久精品国产亚洲av高清涩受| av免费在线观看网站| 美女高潮到喷水免费观看| 成年女人毛片免费观看观看9 | 交换朋友夫妻互换小说| 一级毛片精品| 久热这里只有精品99| 香蕉国产在线看| 亚洲自偷自拍图片 自拍| 另类亚洲欧美激情| 国产欧美日韩综合在线一区二区| 精品少妇一区二区三区视频日本电影| 大陆偷拍与自拍| 欧美中文综合在线视频| 国产男靠女视频免费网站| 电影成人av| 久久 成人 亚洲| 亚洲第一av免费看| 女警被强在线播放| 麻豆成人av在线观看| 亚洲av日韩精品久久久久久密| 中文字幕色久视频| 亚洲中文字幕日韩| 久久精品人人爽人人爽视色| 黄色丝袜av网址大全| 欧美 日韩 精品 国产| 男女之事视频高清在线观看| 日日爽夜夜爽网站| 久久亚洲精品不卡| 亚洲欧美激情综合另类| 精品国产亚洲在线| 国产精品久久久久久人妻精品电影| 久久人人爽av亚洲精品天堂| 国产成人系列免费观看| 欧美日韩一级在线毛片| av不卡在线播放| 18禁国产床啪视频网站| 精品欧美一区二区三区在线| 午夜精品国产一区二区电影| 日韩大码丰满熟妇| 又黄又粗又硬又大视频| www.999成人在线观看| 国产97色在线日韩免费| 欧美日韩亚洲高清精品| av网站免费在线观看视频| 亚洲色图综合在线观看| 亚洲国产看品久久| 操出白浆在线播放| 欧美日韩av久久| 国产精品.久久久| 久久久国产一区二区| 黄色 视频免费看| 最近最新中文字幕大全免费视频| 亚洲av片天天在线观看| 免费在线观看完整版高清| 亚洲精品乱久久久久久| 大型黄色视频在线免费观看| 欧美乱妇无乱码| 好看av亚洲va欧美ⅴa在| 国产精品免费一区二区三区在线 | 一进一出抽搐gif免费好疼 | 成人18禁高潮啪啪吃奶动态图| 亚洲 国产 在线| 视频在线观看一区二区三区| 久久久久国产精品人妻aⅴ院 | 国产精品久久久久成人av| 99热只有精品国产| 亚洲综合色网址| 老汉色∧v一级毛片| 久久99一区二区三区| 黄色毛片三级朝国网站| 99精品欧美一区二区三区四区| 国产97色在线日韩免费| 男人舔女人的私密视频| 高清毛片免费观看视频网站 | 国产乱人伦免费视频| 国产精品.久久久| 正在播放国产对白刺激| 国产精品久久久av美女十八| bbb黄色大片| 国产精品98久久久久久宅男小说| 狠狠婷婷综合久久久久久88av| 黄片播放在线免费| 国产蜜桃级精品一区二区三区 | 1024香蕉在线观看| 免费久久久久久久精品成人欧美视频| av线在线观看网站| 免费看十八禁软件| 亚洲午夜精品一区,二区,三区| 五月开心婷婷网| 欧美国产精品一级二级三级| 波多野结衣av一区二区av| 中文字幕人妻熟女乱码| 在线观看www视频免费| 男女之事视频高清在线观看| 国产精品电影一区二区三区 | xxxhd国产人妻xxx| 午夜福利,免费看| 色综合欧美亚洲国产小说| 亚洲 国产 在线| 桃红色精品国产亚洲av| 下体分泌物呈黄色| 亚洲精品中文字幕一二三四区| 脱女人内裤的视频| 人人妻人人爽人人添夜夜欢视频| 亚洲欧洲精品一区二区精品久久久| 一级a爱片免费观看的视频| 成年人黄色毛片网站| 美女国产高潮福利片在线看| 视频区欧美日本亚洲| 中文字幕最新亚洲高清| 久热这里只有精品99| 亚洲免费av在线视频| 久久香蕉精品热| 国产成人精品无人区| 亚洲午夜理论影院| 亚洲精品一二三| 如日韩欧美国产精品一区二区三区| 欧美最黄视频在线播放免费 | 欧美激情久久久久久爽电影 | 亚洲情色 制服丝袜| 国产精品久久久久久精品古装| 老汉色av国产亚洲站长工具| 欧美成人免费av一区二区三区 | 99精品欧美一区二区三区四区| 久9热在线精品视频| 性少妇av在线| 国产熟女午夜一区二区三区| 啦啦啦在线免费观看视频4| 国产精品久久久久久精品古装| 精品视频人人做人人爽| 日韩人妻精品一区2区三区| 成人18禁在线播放| 激情视频va一区二区三区| 国产欧美日韩精品亚洲av| 成人黄色视频免费在线看| 中文字幕av电影在线播放| 午夜影院日韩av| 一夜夜www| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲一区中文字幕在线| 欧美激情 高清一区二区三区| 国产av又大| 亚洲色图av天堂| 精品免费久久久久久久清纯 | 国产免费现黄频在线看| 丁香六月欧美| 日本精品一区二区三区蜜桃| 国产精品一区二区免费欧美| 宅男免费午夜| 久久午夜综合久久蜜桃| 亚洲第一青青草原| 国产人伦9x9x在线观看| 欧美黄色片欧美黄色片| 国产aⅴ精品一区二区三区波| 国产高清激情床上av| 国产成人欧美在线观看 | 亚洲国产中文字幕在线视频| 亚洲精品久久成人aⅴ小说| 99久久综合精品五月天人人| 亚洲第一av免费看| 色综合欧美亚洲国产小说| 无遮挡黄片免费观看| 日韩一卡2卡3卡4卡2021年| 12—13女人毛片做爰片一| 两性夫妻黄色片| 国产99白浆流出| 这个男人来自地球电影免费观看| 亚洲专区国产一区二区| 国产精品偷伦视频观看了| 亚洲欧美精品综合一区二区三区| 免费人成视频x8x8入口观看| 色婷婷久久久亚洲欧美| 99国产精品一区二区蜜桃av | 国产av又大| 建设人人有责人人尽责人人享有的| a级毛片在线看网站| 曰老女人黄片| 男女免费视频国产| 久久精品亚洲熟妇少妇任你| 母亲3免费完整高清在线观看| a级毛片黄视频| 国产亚洲精品久久久久5区| 亚洲av片天天在线观看| 黄色视频,在线免费观看| 国产成人精品久久二区二区免费| 国产成人精品在线电影| 欧美乱码精品一区二区三区| 精品国产乱码久久久久久男人| 男男h啪啪无遮挡| 亚洲中文av在线| xxxhd国产人妻xxx| 国产三级黄色录像| 国产成人免费无遮挡视频| 十八禁网站免费在线| 国产99久久九九免费精品| 男女床上黄色一级片免费看| 免费在线观看日本一区| 亚洲熟妇熟女久久| 一级毛片精品| 午夜福利在线免费观看网站| 超碰成人久久| 免费观看人在逋| 天天躁日日躁夜夜躁夜夜| 国产免费现黄频在线看| 欧美乱码精品一区二区三区| 午夜亚洲福利在线播放| 曰老女人黄片| 亚洲五月婷婷丁香| 丝袜在线中文字幕| 在线观看日韩欧美| 精品电影一区二区在线| 99久久精品国产亚洲精品| av天堂在线播放| 精品久久久久久电影网| 午夜福利视频在线观看免费| 18禁裸乳无遮挡动漫免费视频| av中文乱码字幕在线| 亚洲免费av在线视频| 亚洲国产精品sss在线观看 | 怎么达到女性高潮| 久久久水蜜桃国产精品网| 一级毛片女人18水好多| 成在线人永久免费视频| 久久精品亚洲精品国产色婷小说| 久久99一区二区三区| 曰老女人黄片| 99riav亚洲国产免费| 精品国产超薄肉色丝袜足j| 男人操女人黄网站| ponron亚洲| 99精品久久久久人妻精品| 飞空精品影院首页| 免费在线观看黄色视频的| 日韩制服丝袜自拍偷拍| 国产成人欧美在线观看 | 91老司机精品| 久久香蕉国产精品| 91九色精品人成在线观看| 99精国产麻豆久久婷婷| 超碰成人久久| 国产黄色免费在线视频| 999精品在线视频| 一区在线观看完整版| 少妇粗大呻吟视频| 热99re8久久精品国产| 亚洲少妇的诱惑av| 老司机福利观看| 两个人免费观看高清视频| 成人三级做爰电影| 久热这里只有精品99| 国产日韩一区二区三区精品不卡| 久久久水蜜桃国产精品网| 在线观看免费高清a一片| 黑丝袜美女国产一区| 欧美性长视频在线观看| 十八禁人妻一区二区| 动漫黄色视频在线观看| 欧美日韩福利视频一区二区| 精品亚洲成a人片在线观看| 香蕉国产在线看| 国产精品久久久久久精品古装| 在线观看日韩欧美| 色在线成人网| 亚洲欧洲精品一区二区精品久久久| 99精品欧美一区二区三区四区| 日日夜夜操网爽| 免费高清在线观看日韩| 中文字幕精品免费在线观看视频| 母亲3免费完整高清在线观看| 亚洲欧美激情在线| 精品国产一区二区久久| 午夜福利视频在线观看免费| 熟女少妇亚洲综合色aaa.| 热99久久久久精品小说推荐| 国产精品电影一区二区三区 | 女同久久另类99精品国产91| 国产熟女午夜一区二区三区| 女人高潮潮喷娇喘18禁视频| 大片电影免费在线观看免费| 亚洲精品美女久久久久99蜜臀| 99国产精品99久久久久| 午夜福利在线免费观看网站| 亚洲自偷自拍图片 自拍| 亚洲精品自拍成人| 超碰97精品在线观看| 女警被强在线播放| 国产在线观看jvid| 色尼玛亚洲综合影院| 在线免费观看的www视频| 国产成人精品久久二区二区免费| 满18在线观看网站| 成人特级黄色片久久久久久久| 在线观看免费午夜福利视频| 亚洲成人免费av在线播放| 一级片'在线观看视频| 欧美激情 高清一区二区三区| 窝窝影院91人妻| 午夜福利在线免费观看网站| 曰老女人黄片| 99国产精品一区二区蜜桃av | 午夜福利在线免费观看网站| 欧美久久黑人一区二区| 免费黄频网站在线观看国产| 久久香蕉激情| 变态另类成人亚洲欧美熟女 | 男男h啪啪无遮挡| 成人永久免费在线观看视频| 久久精品国产清高在天天线| 中出人妻视频一区二区| 国产一区二区三区视频了| 一级片免费观看大全| 欧美激情久久久久久爽电影 | 亚洲va日本ⅴa欧美va伊人久久| 成在线人永久免费视频| 成人精品一区二区免费| 亚洲精品美女久久av网站| 久久影院123| 日韩精品免费视频一区二区三区| 99久久综合精品五月天人人| 亚洲成人免费av在线播放| 老司机午夜十八禁免费视频| 亚洲全国av大片| 国产极品粉嫩免费观看在线| 国产99白浆流出| 亚洲五月天丁香| 亚洲欧美精品综合一区二区三区| 国产极品粉嫩免费观看在线| 久久狼人影院| 精品国内亚洲2022精品成人 | 日日摸夜夜添夜夜添小说| 亚洲精品中文字幕一二三四区| 日韩中文字幕欧美一区二区| 日本撒尿小便嘘嘘汇集6| 一级a爱片免费观看的视频| 丰满迷人的少妇在线观看| 露出奶头的视频| 777久久人妻少妇嫩草av网站| 一级作爱视频免费观看| 成年动漫av网址| 夜夜躁狠狠躁天天躁| 亚洲av第一区精品v没综合| 欧美成人午夜精品| 99久久精品国产亚洲精品| 国产不卡av网站在线观看| 91精品国产国语对白视频| 欧美精品人与动牲交sv欧美| 国产欧美日韩一区二区三区在线| 大片电影免费在线观看免费| av有码第一页| 曰老女人黄片| 精品国产一区二区久久| 少妇被粗大的猛进出69影院| 丝袜美足系列| 欧美日韩成人在线一区二区| 精品亚洲成a人片在线观看| 两个人看的免费小视频| 黄色视频,在线免费观看| 国产亚洲av高清不卡| 亚洲人成伊人成综合网2020| 日韩欧美一区视频在线观看| 在线观看免费高清a一片| 亚洲精华国产精华精| www.自偷自拍.com| 午夜成年电影在线免费观看| 美女高潮到喷水免费观看| 老司机在亚洲福利影院| 免费在线观看黄色视频的| 久久精品熟女亚洲av麻豆精品| 啦啦啦在线免费观看视频4| 久久99一区二区三区| 国产精品久久久久久精品古装| 久久草成人影院| 国产一区二区三区综合在线观看| 欧美精品啪啪一区二区三区| 久9热在线精品视频| 国产真人三级小视频在线观看| 久久青草综合色| 亚洲av熟女| 国产精品综合久久久久久久免费 | 老鸭窝网址在线观看| 成人国产一区最新在线观看| svipshipincom国产片| 色94色欧美一区二区| 午夜激情av网站| 三级毛片av免费| 欧美午夜高清在线| 国产亚洲精品一区二区www | 国产亚洲精品第一综合不卡| 亚洲欧洲精品一区二区精品久久久| 国产男女超爽视频在线观看| 变态另类成人亚洲欧美熟女 | av福利片在线| 人妻 亚洲 视频| 香蕉久久夜色| 欧美成狂野欧美在线观看| 国产在线观看jvid| 亚洲精品av麻豆狂野| av超薄肉色丝袜交足视频| 涩涩av久久男人的天堂| 淫妇啪啪啪对白视频| 亚洲精品中文字幕一二三四区| 欧美 日韩 精品 国产| 午夜成年电影在线免费观看| 99国产综合亚洲精品| 精品国产超薄肉色丝袜足j| 老司机午夜福利在线观看视频| 在线天堂中文资源库| 亚洲午夜精品一区,二区,三区| 日本wwww免费看| 国产免费av片在线观看野外av| 久久精品国产综合久久久| 丝袜人妻中文字幕| 午夜影院日韩av| av国产精品久久久久影院| 母亲3免费完整高清在线观看| 久久精品国产a三级三级三级| 久久国产精品人妻蜜桃| cao死你这个sao货| 黑人巨大精品欧美一区二区蜜桃| 国产高清国产精品国产三级| 黄频高清免费视频| 欧美+亚洲+日韩+国产| 91大片在线观看| 久久精品国产亚洲av高清一级|