• <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à)變換和高維不變子空間
    欧美日本视频| 免费无遮挡裸体视频| 日韩欧美在线二视频| 久久香蕉激情| 最近最新免费中文字幕在线| 久久久久久久久久黄片| 悠悠久久av| 中文字幕av在线有码专区| 久久久国产成人精品二区| 精品欧美国产一区二区三| 免费观看人在逋| 久久天堂一区二区三区四区| 99久久久亚洲精品蜜臀av| 国产亚洲精品一区二区www| 精品久久久久久久人妻蜜臀av| 国产区一区二久久| 欧美黑人巨大hd| 国产亚洲欧美在线一区二区| 午夜福利欧美成人| 999精品在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美午夜高清在线| 18禁观看日本| 777久久人妻少妇嫩草av网站| 高潮久久久久久久久久久不卡| 亚洲成a人片在线一区二区| 1024视频免费在线观看| 精品久久久久久久毛片微露脸| 国产成人系列免费观看| 国产精品美女特级片免费视频播放器 | 国产av一区二区精品久久| 19禁男女啪啪无遮挡网站| 90打野战视频偷拍视频| 岛国在线免费视频观看| 桃红色精品国产亚洲av| 老司机午夜十八禁免费视频| 1024手机看黄色片| 最近最新免费中文字幕在线| 久久精品国产亚洲av高清一级| 亚洲精品中文字幕一二三四区| 久久中文字幕一级| 国产午夜精品久久久久久| 久久亚洲精品不卡| 村上凉子中文字幕在线| 日本成人三级电影网站| 日韩精品中文字幕看吧| 91九色精品人成在线观看| 久久久精品大字幕| 男女那种视频在线观看| 国产精品永久免费网站| 女人被狂操c到高潮| 国产麻豆成人av免费视频| 毛片女人毛片| 桃红色精品国产亚洲av| 亚洲免费av在线视频| 久久精品91无色码中文字幕| 国内毛片毛片毛片毛片毛片| 在线观看美女被高潮喷水网站 | 国产精品日韩av在线免费观看| 极品教师在线免费播放| av在线天堂中文字幕| 香蕉丝袜av| 叶爱在线成人免费视频播放| 亚洲av第一区精品v没综合| 欧美成人性av电影在线观看| 男女之事视频高清在线观看| 午夜亚洲福利在线播放| а√天堂www在线а√下载| 亚洲人成77777在线视频| 91成年电影在线观看| 床上黄色一级片| 国产精品,欧美在线| 婷婷六月久久综合丁香| 美女大奶头视频| 午夜福利在线在线| 老司机靠b影院| 人妻夜夜爽99麻豆av| 久热爱精品视频在线9| 国内精品一区二区在线观看| 久久国产精品人妻蜜桃| 欧美丝袜亚洲另类 | 真人做人爱边吃奶动态| 午夜福利成人在线免费观看| 国产麻豆成人av免费视频| 国产熟女午夜一区二区三区| 一个人免费在线观看电影 | 午夜免费观看网址| 国产在线观看jvid| 午夜福利高清视频| 日韩中文字幕欧美一区二区| 成人高潮视频无遮挡免费网站| 白带黄色成豆腐渣| 亚洲自偷自拍图片 自拍| 人妻丰满熟妇av一区二区三区| 午夜久久久久精精品| 亚洲欧美精品综合久久99| 亚洲欧美激情综合另类| 欧美色欧美亚洲另类二区| 别揉我奶头~嗯~啊~动态视频| 亚洲国产欧美一区二区综合| 狠狠狠狠99中文字幕| 一个人免费在线观看电影 | 美女免费视频网站| 亚洲一区二区三区色噜噜| 黑人巨大精品欧美一区二区mp4| 黄色成人免费大全| 老司机在亚洲福利影院| 国产精品1区2区在线观看.| 老熟妇乱子伦视频在线观看| 在线国产一区二区在线| 黑人欧美特级aaaaaa片| 国产人伦9x9x在线观看| 狂野欧美激情性xxxx| 欧美一级毛片孕妇| 久久精品国产99精品国产亚洲性色| av片东京热男人的天堂| 亚洲人与动物交配视频| 女警被强在线播放| av在线天堂中文字幕| 大型黄色视频在线免费观看| 亚洲国产高清在线一区二区三| 黄色成人免费大全| 三级国产精品欧美在线观看 | 国产精品亚洲美女久久久| 国产精品野战在线观看| 欧美大码av| 级片在线观看| 日本在线视频免费播放| 精品久久久久久成人av| 精品久久久久久成人av| 国产精品永久免费网站| 丰满人妻熟妇乱又伦精品不卡| 久久午夜综合久久蜜桃| 亚洲欧洲精品一区二区精品久久久| 亚洲在线自拍视频| 国产激情欧美一区二区| 中文字幕高清在线视频| 久久精品亚洲精品国产色婷小说| 亚洲一区中文字幕在线| 婷婷精品国产亚洲av| 女同久久另类99精品国产91| 三级国产精品欧美在线观看 | 在线观看66精品国产| 身体一侧抽搐| 亚洲欧美日韩无卡精品| 婷婷六月久久综合丁香| 黄片小视频在线播放| 国产黄a三级三级三级人| 色综合站精品国产| 黄频高清免费视频| 国产精品久久久久久人妻精品电影| 亚洲精华国产精华精| 中国美女看黄片| 又大又爽又粗| 老司机午夜十八禁免费视频| 美女扒开内裤让男人捅视频| 观看免费一级毛片| 看片在线看免费视频| 好看av亚洲va欧美ⅴa在| 日韩精品中文字幕看吧| 亚洲 欧美 日韩 在线 免费| 18禁观看日本| 婷婷精品国产亚洲av在线| 中文字幕人成人乱码亚洲影| 欧美黑人精品巨大| 五月玫瑰六月丁香| 亚洲熟妇熟女久久| 俄罗斯特黄特色一大片| 免费人成视频x8x8入口观看| 亚洲国产欧美一区二区综合| 黄色片一级片一级黄色片| www国产在线视频色| 国产高清videossex| 人妻丰满熟妇av一区二区三区| 看黄色毛片网站| 国产精品一区二区三区四区免费观看 | 18禁黄网站禁片午夜丰满| 亚洲人成网站高清观看| 色综合婷婷激情| 岛国在线观看网站| 啦啦啦韩国在线观看视频| 国产成+人综合+亚洲专区| 久久精品亚洲精品国产色婷小说| 又黄又爽又免费观看的视频| 亚洲精品在线美女| 国产真人三级小视频在线观看| 日韩欧美三级三区| 亚洲国产欧美人成| 一进一出好大好爽视频| 久久这里只有精品中国| 国产高清视频在线播放一区| 欧美+亚洲+日韩+国产| 国产男靠女视频免费网站| 欧美黑人精品巨大| 中文字幕人妻丝袜一区二区| 午夜福利免费观看在线| 国产成人精品无人区| 一级作爱视频免费观看| 国产视频内射| 午夜久久久久精精品| 日韩精品中文字幕看吧| 午夜视频精品福利| 五月伊人婷婷丁香| 五月玫瑰六月丁香| 视频区欧美日本亚洲| 99国产精品一区二区蜜桃av| 在线观看免费日韩欧美大片| 老司机靠b影院| 婷婷精品国产亚洲av在线| 久久国产乱子伦精品免费另类| 18禁国产床啪视频网站| 久久中文字幕人妻熟女| 99在线视频只有这里精品首页| 国产精品 欧美亚洲| 欧美激情久久久久久爽电影| 人人妻,人人澡人人爽秒播| 亚洲片人在线观看| 亚洲av电影在线进入| 中国美女看黄片| 熟女少妇亚洲综合色aaa.| 欧美高清成人免费视频www| 熟女电影av网| 在线观看一区二区三区| 国产精品日韩av在线免费观看| 免费搜索国产男女视频| 亚洲成人久久爱视频| 成人18禁高潮啪啪吃奶动态图| 91在线观看av| 在线观看www视频免费| 国产真人三级小视频在线观看| 国产欧美日韩精品亚洲av| 岛国在线免费视频观看| 欧美黄色淫秽网站| 亚洲人成网站高清观看| 不卡av一区二区三区| 久久婷婷成人综合色麻豆| 国产一区二区在线观看日韩 | 99精品欧美一区二区三区四区| 给我免费播放毛片高清在线观看| 美女 人体艺术 gogo| 丰满人妻熟妇乱又伦精品不卡| 欧美精品啪啪一区二区三区| 国产男靠女视频免费网站| 国产高清视频在线观看网站| 最近最新免费中文字幕在线| 国产精品永久免费网站| 免费在线观看视频国产中文字幕亚洲| 亚洲成人中文字幕在线播放| 国产亚洲精品久久久久5区| 给我免费播放毛片高清在线观看| 亚洲av成人精品一区久久| 亚洲 国产 在线| 国产精品av久久久久免费| 色av中文字幕| 在线a可以看的网站| 99精品欧美一区二区三区四区| 可以在线观看的亚洲视频| 夜夜躁狠狠躁天天躁| 麻豆成人av在线观看| 色噜噜av男人的天堂激情| 亚洲av电影在线进入| 午夜激情福利司机影院| 精品久久久久久久人妻蜜臀av| 国产av在哪里看| 露出奶头的视频| 亚洲成av人片免费观看| 校园春色视频在线观看| 久久精品国产综合久久久| 妹子高潮喷水视频| 99re在线观看精品视频| 亚洲精品国产精品久久久不卡| 麻豆av在线久日| 99热这里只有精品一区 | 国产成人欧美在线观看| 国产精品乱码一区二三区的特点| av在线播放免费不卡| 亚洲成av人片免费观看| 国产精品永久免费网站| 日韩国内少妇激情av| 成人国产一区最新在线观看| 国产高清视频在线观看网站| 十八禁人妻一区二区| 可以在线观看毛片的网站| 亚洲av成人不卡在线观看播放网| 免费观看人在逋| 91大片在线观看| 少妇粗大呻吟视频| 麻豆国产97在线/欧美 | 久久伊人香网站| 精品福利观看| 日韩 欧美 亚洲 中文字幕| 听说在线观看完整版免费高清| 国模一区二区三区四区视频 | 女人被狂操c到高潮| 亚洲午夜理论影院| 黄色 视频免费看| 深夜精品福利| 成人手机av| 欧美大码av| 欧美激情久久久久久爽电影| 在线视频色国产色| 视频区欧美日本亚洲| 日韩有码中文字幕| 日韩大尺度精品在线看网址| 欧美日韩精品网址| 亚洲人成电影免费在线| 久久久久九九精品影院| 色综合欧美亚洲国产小说| 好看av亚洲va欧美ⅴa在| 神马国产精品三级电影在线观看 | 成年免费大片在线观看| netflix在线观看网站| 精品久久久久久,| 亚洲黑人精品在线| 成人三级做爰电影| 在线观看舔阴道视频| 亚洲精品一区av在线观看| 国产欧美日韩精品亚洲av| 毛片女人毛片| svipshipincom国产片| 亚洲av五月六月丁香网| 日韩精品免费视频一区二区三区| 午夜亚洲福利在线播放| 中出人妻视频一区二区| 成人国产综合亚洲| 国产又色又爽无遮挡免费看| 两个人视频免费观看高清| 欧美中文日本在线观看视频| 一本大道久久a久久精品| 亚洲国产精品sss在线观看| 精品国产乱码久久久久久男人| 夜夜夜夜夜久久久久| 欧美乱色亚洲激情| 2021天堂中文幕一二区在线观| 一个人免费在线观看的高清视频| 久久精品国产99精品国产亚洲性色| 精品福利观看| 国产精品日韩av在线免费观看| 国产av又大| 淫妇啪啪啪对白视频| 国产精品综合久久久久久久免费| 欧美性猛交黑人性爽| 黄色片一级片一级黄色片| 久久精品亚洲精品国产色婷小说| 97人妻精品一区二区三区麻豆| 亚洲性夜色夜夜综合| 久久久精品国产亚洲av高清涩受| 中文字幕高清在线视频| 亚洲av成人不卡在线观看播放网| videosex国产| 国产69精品久久久久777片 | 免费在线观看日本一区| 12—13女人毛片做爰片一| 国产黄a三级三级三级人| bbb黄色大片| 在线观看舔阴道视频| 亚洲成人中文字幕在线播放| 欧美在线一区亚洲| 黑人欧美特级aaaaaa片| 日韩欧美精品v在线| 岛国在线免费视频观看| 成人亚洲精品av一区二区| 亚洲自偷自拍图片 自拍| 香蕉久久夜色| 999久久久国产精品视频| av视频在线观看入口| 久久伊人香网站| 三级男女做爰猛烈吃奶摸视频| 精品国产美女av久久久久小说| 国产精品久久久久久人妻精品电影| 亚洲欧洲精品一区二区精品久久久| 久久久久久免费高清国产稀缺| 丰满的人妻完整版| 成人手机av| 一区二区三区国产精品乱码| 日韩免费av在线播放| 精品国产乱码久久久久久男人| 99久久综合精品五月天人人| 天天一区二区日本电影三级| 国产一区二区在线观看日韩 | 久久精品国产99精品国产亚洲性色| 一本大道久久a久久精品| 国产蜜桃级精品一区二区三区| 18禁裸乳无遮挡免费网站照片| 亚洲 欧美一区二区三区| av国产免费在线观看| 国产精华一区二区三区| 可以在线观看毛片的网站| 男人舔奶头视频| 两个人视频免费观看高清| 亚洲欧洲精品一区二区精品久久久| 悠悠久久av| 久久久久久大精品| 亚洲天堂国产精品一区在线| 欧美一区二区精品小视频在线| 人成视频在线观看免费观看| 人妻久久中文字幕网| 妹子高潮喷水视频| АⅤ资源中文在线天堂| 可以免费在线观看a视频的电影网站| 久久国产精品影院| 男女床上黄色一级片免费看| 在线观看免费午夜福利视频| 国产精品乱码一区二三区的特点| 1024香蕉在线观看| 最近在线观看免费完整版| 黑人操中国人逼视频| 亚洲美女视频黄频| 麻豆成人午夜福利视频| 熟妇人妻久久中文字幕3abv| 男女做爰动态图高潮gif福利片| a在线观看视频网站| 日韩国内少妇激情av| 香蕉av资源在线| 国产麻豆成人av免费视频| 亚洲成av人片在线播放无| 美女午夜性视频免费| 久久性视频一级片| 亚洲国产欧美网| 欧美黑人精品巨大| 欧美人与性动交α欧美精品济南到| 亚洲国产日韩欧美精品在线观看 | 99久久精品国产亚洲精品| 婷婷亚洲欧美| 男女床上黄色一级片免费看| 人妻久久中文字幕网| 国产精品影院久久| 夜夜爽天天搞| 欧美一级a爱片免费观看看 | 精品熟女少妇八av免费久了| 人人妻人人看人人澡| 亚洲国产精品sss在线观看| a在线观看视频网站| 午夜日韩欧美国产| 国产欧美日韩一区二区三| 日本一区二区免费在线视频| 成人国产综合亚洲| 久久久久久久久久黄片| 精品无人区乱码1区二区| 在线观看66精品国产| 久久香蕉国产精品| 亚洲男人的天堂狠狠| 97碰自拍视频| 一边摸一边抽搐一进一小说| 在线观看舔阴道视频| 国产精品久久久人人做人人爽| bbb黄色大片| 成人国产综合亚洲| 日韩三级视频一区二区三区| 免费看日本二区| 久久亚洲精品不卡| av国产免费在线观看| 一个人观看的视频www高清免费观看 | 日韩欧美 国产精品| 免费观看人在逋| 国产成人av激情在线播放| 亚洲全国av大片| 又大又爽又粗| 欧美日本视频| 国产精品,欧美在线| 成人国语在线视频| 999精品在线视频| 亚洲欧美一区二区三区黑人| 怎么达到女性高潮| 一级a爱片免费观看的视频| 制服诱惑二区| 亚洲最大成人中文| 国产欧美日韩精品亚洲av| 校园春色视频在线观看| 伊人久久大香线蕉亚洲五| 在线观看免费日韩欧美大片| 欧美激情久久久久久爽电影| 亚洲精品美女久久av网站| 亚洲国产中文字幕在线视频| 三级国产精品欧美在线观看 | 国产精品九九99| 欧洲精品卡2卡3卡4卡5卡区| 国产三级中文精品| 日韩大码丰满熟妇| 无限看片的www在线观看| 亚洲国产欧洲综合997久久,| 成人国语在线视频| 亚洲精品中文字幕在线视频| 亚洲av熟女| 麻豆一二三区av精品| 久久精品国产清高在天天线| 免费在线观看亚洲国产| 欧美成人午夜精品| 欧美性猛交╳xxx乱大交人| 久久人妻av系列| 久久中文看片网| 精品久久久久久久久久久久久| 日韩欧美精品v在线| 18禁黄网站禁片午夜丰满| 一二三四在线观看免费中文在| 国产精品av久久久久免费| 精品久久久久久久久久免费视频| 巨乳人妻的诱惑在线观看| 在线观看日韩欧美| 亚洲专区国产一区二区| 国产精品自产拍在线观看55亚洲| 欧美久久黑人一区二区| 中文字幕人成人乱码亚洲影| 久久这里只有精品19| 动漫黄色视频在线观看| 欧美在线一区亚洲| 在线观看舔阴道视频| 久久香蕉国产精品| 日韩欧美国产在线观看| 超碰成人久久| 90打野战视频偷拍视频| 中文字幕最新亚洲高清| 国产精品一及| 黄色丝袜av网址大全| av国产免费在线观看| 别揉我奶头~嗯~啊~动态视频| 一边摸一边做爽爽视频免费| 日韩欧美免费精品| 成人18禁在线播放| 亚洲专区中文字幕在线| 国产欧美日韩一区二区精品| 国产探花在线观看一区二区| 国产精品影院久久| 久9热在线精品视频| 欧美三级亚洲精品| av超薄肉色丝袜交足视频| 亚洲最大成人中文| 我要搜黄色片| 一级毛片高清免费大全| 亚洲精品美女久久av网站| 中文字幕久久专区| 午夜免费观看网址| 色尼玛亚洲综合影院| av欧美777| 亚洲欧美精品综合久久99| 精品人妻1区二区| 无人区码免费观看不卡| 久久久久久九九精品二区国产 | 成人高潮视频无遮挡免费网站| 国产精品爽爽va在线观看网站| 看黄色毛片网站| 最近视频中文字幕2019在线8| 亚洲国产欧美网| 90打野战视频偷拍视频| 真人做人爱边吃奶动态| 国产午夜精品久久久久久| 黑人欧美特级aaaaaa片| 国产精品免费一区二区三区在线| 色尼玛亚洲综合影院| 久久精品夜夜夜夜夜久久蜜豆 | 特级一级黄色大片| 中文字幕熟女人妻在线| 一级毛片高清免费大全| 日韩三级视频一区二区三区| 黄色a级毛片大全视频| 18禁国产床啪视频网站| 正在播放国产对白刺激| 国产av一区二区精品久久| 国产精品久久电影中文字幕| 丰满的人妻完整版| 国产成人啪精品午夜网站| 中文字幕人成人乱码亚洲影| 亚洲第一欧美日韩一区二区三区| 日韩成人在线观看一区二区三区| 成年免费大片在线观看| 午夜免费成人在线视频| 精品一区二区三区四区五区乱码| av天堂在线播放| 女人被狂操c到高潮| 在线视频色国产色| 欧美成人免费av一区二区三区| av在线播放免费不卡| 国内久久婷婷六月综合欲色啪| 国产不卡一卡二| svipshipincom国产片| 制服人妻中文乱码| 久久99热这里只有精品18| 精品久久久久久久久久久久久| 久久精品国产亚洲av香蕉五月| √禁漫天堂资源中文www| 一区福利在线观看| 午夜成年电影在线免费观看| 国内揄拍国产精品人妻在线| 级片在线观看| 免费在线观看完整版高清| 久久人妻av系列| 美女黄网站色视频| 淫妇啪啪啪对白视频| 哪里可以看免费的av片| 精品久久久久久久人妻蜜臀av| 制服丝袜大香蕉在线| 变态另类成人亚洲欧美熟女| 国产欧美日韩一区二区精品| 精品免费久久久久久久清纯| 无限看片的www在线观看| 一夜夜www| 亚洲av五月六月丁香网| 精品久久久久久久末码| 在线观看免费视频日本深夜| 国产精华一区二区三区| 999精品在线视频| 无人区码免费观看不卡| 亚洲真实伦在线观看| 日本撒尿小便嘘嘘汇集6| 久久国产乱子伦精品免费另类| 午夜福利成人在线免费观看| 99久久精品热视频| 制服人妻中文乱码| 国产欧美日韩精品亚洲av| 一区二区三区国产精品乱码| 久久久精品欧美日韩精品| 欧美乱色亚洲激情| 狠狠狠狠99中文字幕| 久久国产精品影院| 久久午夜亚洲精品久久|