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

    基于模糊C-均值聚類(lèi)算法的臺(tái)區(qū)電壓與用戶(hù)關(guān)系辨識(shí)

    2021-07-21 14:07:02曾順奇吳杰康蔡志宏
    四川電力技術(shù) 2021年3期
    關(guān)鍵詞:皮爾遜臺(tái)區(qū)特征值

    曾順奇,吳杰康,李 欣,蔡志宏

    (1.廣東電網(wǎng)有限責(zé)任公司廣州供電局,廣東 廣州 510000;2.廣東工業(yè)大學(xué)自動(dòng)化學(xué)院,廣東 廣州 510006)

    0 引 言

    近年來(lái),隨著中國(guó)城鄉(xiāng)經(jīng)濟(jì)建設(shè)迅猛發(fā)展,用電負(fù)荷激增,有些小用戶(hù)發(fā)展成為大功率用戶(hù),超過(guò)了線(xiàn)路本身預(yù)留的空間,導(dǎo)致臺(tái)區(qū)電壓在用電高峰出現(xiàn)欠壓現(xiàn)象[1]。配電網(wǎng)網(wǎng)架結(jié)構(gòu)的相對(duì)滯后和迅速增長(zhǎng)的用電需求之間的矛盾越來(lái)越突出[2],導(dǎo)致臺(tái)區(qū)電壓波動(dòng)較大,嚴(yán)重影響用戶(hù)用電體驗(yàn)并降低配電網(wǎng)供電可靠性,給供電公司帶來(lái)了嚴(yán)峻的考驗(yàn)。針對(duì)以上問(wèn)題,業(yè)界學(xué)者開(kāi)展了針對(duì)性研究,例如:從原理上對(duì)低壓?jiǎn)栴}產(chǎn)生的原因深入分析,得到電壓偏低的原理[3];構(gòu)建多指標(biāo)的臺(tái)區(qū)低電壓成因識(shí)別指標(biāo)體系,建立基于支持向量機(jī)多分類(lèi)器的臺(tái)區(qū)低壓成因識(shí)別模型[2];通過(guò)人工排查、登記、數(shù)據(jù)預(yù)測(cè)和定期分析低電壓,動(dòng)態(tài)識(shí)別低電壓[4];采用改進(jìn)粒子群優(yōu)化(particle swarm optimization, PSO)算法優(yōu)化電容配置,改善負(fù)荷接入配電網(wǎng)對(duì)電網(wǎng)造成的沖擊[5]。

    以上文獻(xiàn)從低電壓成因以及治理方案進(jìn)行相關(guān)研究,但是充分利用大數(shù)據(jù)技術(shù)識(shí)別臺(tái)區(qū)電壓與用戶(hù)用電行為之間關(guān)系相關(guān)研究較少。因此,提出一種新的方法辨識(shí)臺(tái)區(qū)電壓與用戶(hù)的關(guān)系:通過(guò)智能電能表提供的海量用電數(shù)據(jù),首先,利用插值法對(duì)數(shù)據(jù)采樣過(guò)程中產(chǎn)生的缺值進(jìn)行補(bǔ)全;然后,對(duì)其進(jìn)行數(shù)據(jù)特征提取,通過(guò)主成分分析法(principal component analysis, PCA)進(jìn)行數(shù)據(jù)降維,大大減少計(jì)算量,提高辨識(shí)的時(shí)效性和聚類(lèi)的收斂性;最后,采用模糊C均值(fuzzy C-means, FCM)算法對(duì)用戶(hù)進(jìn)行分類(lèi),按照周期內(nèi)用戶(hù)用電有功值分為大、中、小3個(gè)等級(jí)用戶(hù),并結(jié)合皮爾遜相關(guān)系數(shù)分析各個(gè)等級(jí)用戶(hù)在峰-谷-平區(qū)間對(duì)臺(tái)區(qū)電壓的影響水平,進(jìn)而確定不同負(fù)荷等級(jí)用戶(hù)與配電臺(tái)區(qū)電壓之間的關(guān)系。

    1 臺(tái)戶(hù)關(guān)系

    1.1 臺(tái)戶(hù)拓?fù)浣Y(jié)構(gòu)

    配電網(wǎng)終端由變壓器、配電箱和配有智能電能表的用戶(hù)端組成。一般每一個(gè)智能電能表箱都有三相進(jìn)線(xiàn),由智能電能表箱分出A、B、C三相連接各個(gè)用戶(hù),一個(gè)智能電能表箱下可能存在單相用戶(hù),也有三相用戶(hù),這里只討論單相下的所有用戶(hù)。相與相之間電壓互不影響,每一相只受該相所接用戶(hù)的用電情況影響,具體的臺(tái)戶(hù)拓?fù)浣Y(jié)構(gòu)如圖1所示。

    圖2為臺(tái)戶(hù)等值電路,圖中Μ3點(diǎn)為變壓器A相低壓側(cè)電壓。通常低壓側(cè)的電壓與用戶(hù)用電功率、高壓側(cè)擋位和負(fù)荷功率有關(guān)。高壓側(cè)的擋位分別為±5%和0。由潮流計(jì)算可得Μ3點(diǎn)的電壓:

    圖2 臺(tái)戶(hù)等值電路

    VM3=VM1-ΔV

    (1)

    (2)

    下面主要研究在臺(tái)戶(hù)關(guān)系中,不同用戶(hù)的用電行為對(duì)應(yīng)的峰-谷-平區(qū)間與導(dǎo)致變壓器低壓側(cè)電壓發(fā)生偏移二者之間的關(guān)系。

    1.2 臺(tái)戶(hù)關(guān)系

    一個(gè)臺(tái)區(qū)下某一臺(tái)變壓器中一條單相線(xiàn)路連接了不同用電等級(jí)的用戶(hù),其中某一些用戶(hù)的用電行為在峰-谷-平期影響變壓器低壓側(cè)電壓,使其產(chǎn)生波動(dòng)。研究用戶(hù)用電與臺(tái)區(qū)變壓器低壓側(cè)電壓波動(dòng)之間的關(guān)系是基于臺(tái)區(qū)和用戶(hù)的歷史數(shù)據(jù)。用戶(hù)有功功率數(shù)據(jù)由智能電能表獲取,采集終端的時(shí)間分辨率為15 min,采集時(shí)長(zhǎng)為7 d,共有672個(gè)用戶(hù)側(cè)有功功率數(shù)據(jù)點(diǎn)。從用電采集系統(tǒng)歷史數(shù)據(jù)中獲取并構(gòu)建臺(tái)區(qū)用戶(hù)有功功率的數(shù)據(jù)矩陣X∈RN×D:

    (3)

    式中:X為N個(gè)用戶(hù)在采集區(qū)間內(nèi)D個(gè)采樣點(diǎn)所組成的矩陣;xi,tj為臺(tái)區(qū)用戶(hù)i在tj時(shí)刻的有功功率測(cè)量值;N為臺(tái)區(qū)下的用戶(hù)總數(shù);D為采樣時(shí)段內(nèi)臺(tái)區(qū)用戶(hù)有功功率的采樣點(diǎn)數(shù),即臺(tái)區(qū)用戶(hù)原始有功功率數(shù)據(jù)集的特征維度。X的列向量Xtj為臺(tái)區(qū)用戶(hù)在tj時(shí)刻的有功測(cè)量向量;X的行向量Xi為單個(gè)臺(tái)區(qū)用戶(hù)i在采集時(shí)段內(nèi)的有功測(cè)量向量,其單位為kW。

    假設(shè)以上N個(gè)用戶(hù)均屬于同一個(gè)臺(tái)區(qū)下某個(gè)變壓器的A相,按照D=96×7獲取A相電壓數(shù)據(jù),構(gòu)成變壓器低壓側(cè)電壓時(shí)間序列矩陣:

    V=[vt1,vt2,…,vtD]

    (4)

    式中:V為變壓器低壓側(cè)A相在采集周期內(nèi)的電壓向量;vtD為該相D時(shí)刻的測(cè)量電壓值。

    2 臺(tái)戶(hù)數(shù)據(jù)處理與修補(bǔ)

    2.1 功率缺值修補(bǔ)的插值法

    智能電能表采集數(shù)據(jù)時(shí)間跨度大,用戶(hù)的原件安裝、調(diào)試可能存在差異,加之停電或者采集失敗等因素,往往容易導(dǎo)致采集到的電力數(shù)據(jù)出現(xiàn)缺失現(xiàn)象。由于同一用戶(hù)有功數(shù)據(jù)與時(shí)間呈線(xiàn)性關(guān)系,可采用相鄰點(diǎn)(tt-1,xi,t-1)與(tt+1,xi,t+1)估計(jì)中間點(diǎn)的值。假設(shè)xi,t=f(t)呈線(xiàn)性關(guān)系,通過(guò)已知的兩點(diǎn)代入線(xiàn)性方程估算缺值xi,t。

    2.2 臺(tái)戶(hù)數(shù)據(jù)的皮爾遜相關(guān)系數(shù)

    應(yīng)用皮爾遜相關(guān)系數(shù)分析負(fù)載與臺(tái)區(qū)間電壓的關(guān)系?,F(xiàn)有臺(tái)區(qū)D個(gè)時(shí)間點(diǎn)的電壓數(shù)據(jù):V=[v1,v2,…,vD]。由RN×D通過(guò)皮爾遜相關(guān)系數(shù)計(jì)算公式計(jì)算各個(gè)用戶(hù)與臺(tái)區(qū)電壓之間的關(guān)聯(lián)程度。

    皮爾遜相關(guān)系數(shù)計(jì)算式為

    (5)

    式中:Ε(·)為該向量的期望值;cov(·)為兩個(gè)向量的協(xié)方差;ρ(V,Xi)取值范圍為[-1,1],小于0時(shí)為負(fù)相關(guān),大于0時(shí)為正相關(guān),當(dāng)且僅當(dāng)V與Xi有嚴(yán)格線(xiàn)性關(guān)系時(shí)取±1。

    表1為基于皮爾遜相關(guān)系數(shù)的相關(guān)強(qiáng)度估計(jì),通過(guò)表1相關(guān)系數(shù)范圍判斷變量之間的相關(guān)強(qiáng)度。

    表1 基于皮爾遜相關(guān)系數(shù)的相關(guān)強(qiáng)度估計(jì)

    2.3 數(shù)據(jù)降維處理的主成分分析方法

    在采集設(shè)備覆蓋率和采集成功率較高的臺(tái)區(qū),采集的有功數(shù)據(jù)往往分布較為集中,但是采集數(shù)據(jù)大多精度不高而且存在噪聲或者冗余。因此,不能簡(jiǎn)單使用加權(quán)平均法或者求和的方法進(jìn)行用戶(hù)等級(jí)區(qū)分。需要對(duì)采集到的數(shù)據(jù)特征進(jìn)行進(jìn)一步挖掘??梢圆捎肞CA方法進(jìn)行數(shù)據(jù)的特征提取和降維處理。經(jīng)過(guò)降維處理后的數(shù)據(jù),仍保留原始數(shù)據(jù)的分布特性。高維數(shù)據(jù)在數(shù)據(jù)分析中計(jì)算量大,時(shí)效性能不好,利用模糊聚類(lèi)算法對(duì)用戶(hù)等級(jí)分類(lèi)時(shí)可能會(huì)出現(xiàn)維度災(zāi)難、計(jì)算量大、計(jì)算時(shí)間長(zhǎng)等問(wèn)題。通過(guò)PCA法緩解這些問(wèn)題,盡可能解釋變量具有相關(guān)性的高維數(shù)據(jù)集。通過(guò)計(jì)算數(shù)據(jù)的協(xié)方差矩陣cov(X),尋找一個(gè)能反映出原有數(shù)據(jù)特征95%的特征矩陣,將原始高維空間過(guò)渡到低維空間,實(shí)際將數(shù)據(jù)映射到一個(gè)低維子空間,降維成一個(gè)線(xiàn)性無(wú)關(guān)的低維數(shù)據(jù)集。PCA具體分析過(guò)程可參考文獻(xiàn)[10]。

    PCA數(shù)據(jù)降維的具體步驟如下:

    1)對(duì)原始數(shù)據(jù)矩陣X進(jìn)行轉(zhuǎn)置處理得到XT,得到的矩陣行表示每一個(gè)用戶(hù)tj時(shí)刻的特征,共有D維特征。對(duì)XT零均值處理得到矩陣X′為

    (6)

    2)使用線(xiàn)性變換得到協(xié)方差矩陣為

    (7)

    3)通過(guò)求解|λI-R|=0得到特征值λ。最后計(jì)算特征值的累計(jì)貢獻(xiàn)率[10]:

    (8)

    式中,βi為累計(jì)大于95%特征值組成成分(λ1,λ2,…,λk)。對(duì)數(shù)據(jù)進(jìn)行分析后得到k值。同時(shí)得到特征值對(duì)應(yīng)的特征矩陣UD×D=[u1,u2,…,uD 1]。提取前k個(gè)特征值對(duì)應(yīng)的特征向量組成降維矩陣UD×k。降維矩陣Zn×k可表示為

    Zn×k=Xn×DUD×k=

    (9)

    式中,Zn×k為用k維的數(shù)據(jù)表示原始D維數(shù)據(jù)。

    Zn×k不僅降低了數(shù)據(jù)的冗余度,對(duì)提高用戶(hù)分類(lèi)速度有很大的作用。

    3 用戶(hù)分類(lèi)

    在臺(tái)區(qū)電壓與用戶(hù)用電關(guān)系辨識(shí)中,不同用電等級(jí)的用戶(hù)與臺(tái)區(qū)電壓具有不同的關(guān)聯(lián)程度。假設(shè)有數(shù)據(jù)集X,把數(shù)據(jù)分為c類(lèi),對(duì)應(yīng)就有c個(gè)類(lèi)中心C1,C2,…,Cc。把每個(gè)樣本j按照預(yù)設(shè)的類(lèi)簇歸類(lèi),得到用電等級(jí)分別為大、中、小3個(gè)等級(jí)的用戶(hù)。根據(jù)皮爾遜相關(guān)系數(shù)的物理含義,對(duì)不同等級(jí)的用戶(hù)分析其峰-谷-平期間對(duì)臺(tái)區(qū)電壓波動(dòng)的影響,得出哪些用戶(hù)在特定的用電期間與臺(tái)區(qū)電壓波動(dòng)有較強(qiáng)的關(guān)聯(lián)性。

    3.1 模糊聚類(lèi)分析

    模糊聚類(lèi)是在傳統(tǒng)k均值硬聚類(lèi)的基礎(chǔ)上加上隸屬度模糊概念的一種聚類(lèi)分析方法。假設(shè)每個(gè)樣本j屬于某一類(lèi)i存在隸屬度uij的關(guān)系,由此構(gòu)建FCM目標(biāo)函數(shù)及其約束條件:

    (10)

    (11)

    式中:m為一個(gè)隸屬度的因子;n為所有樣本個(gè)數(shù)。

    式(11)為約束條件,說(shuō)明一個(gè)樣本屬于所有類(lèi)的隸屬度之和恒為1。

    3.2 模糊聚類(lèi)中心

    采用模糊聚類(lèi)分析的方法,選取聚類(lèi)中心個(gè)數(shù)c,把所有用戶(hù)分為3類(lèi):大用戶(hù)、中等用戶(hù)及小用戶(hù)。構(gòu)建數(shù)據(jù)集的聚類(lèi)中心:Ci={ci,1,ci,2,…,ci,k}(i=1,2,3),其中ci,k為第i個(gè)聚類(lèi)中心第k維的聚類(lèi)中心特征值。

    3.3 最優(yōu)模糊聚類(lèi)中心矩陣迭代計(jì)算

    采用模糊聚類(lèi)分析方法,通過(guò)迭代計(jì)算,確定模糊C-均值聚類(lèi)矩陣的最優(yōu)聚類(lèi)中心,確定用戶(hù)所屬類(lèi)簇。具體迭代步驟如下:

    1)初始化:確定類(lèi)別數(shù)C、隸屬度因子m、迭代停止誤差ε以及最大迭代次數(shù)(LOOP)。

    2)初始化聚類(lèi)中心P。

    3)計(jì)算初始距離矩陣D。

    4)按式(12)更新用戶(hù)與聚類(lèi)中心的隸屬度。

    (12)

    式中,d(·)為樣本點(diǎn)到樣本中心的距離函數(shù)。若樣本距離中心距離為0,隸屬度設(shè)為1,即完全屬于該中心,否則按照推導(dǎo)式子確定隸屬度。

    5)更新聚類(lèi)中心。

    (13)

    6)重新計(jì)算距離矩陣,并計(jì)算目標(biāo)函數(shù)值。

    (14)

    7)若達(dá)到最大迭代次數(shù)或者前后兩次J的絕對(duì)差小于預(yù)設(shè)誤差ε則結(jié)束,否則轉(zhuǎn)步驟4。具體的模糊聚類(lèi)迭代流程如圖3所示。

    圖3 模糊聚類(lèi)流程

    4 實(shí)例仿真與分析

    以廣州某小區(qū)為實(shí)例,利用其配電系統(tǒng)用戶(hù)用電歷史數(shù)據(jù),選取其所屬臺(tái)區(qū)某個(gè)變壓器低壓側(cè)A相下91戶(hù)在6月21日至6月27日間的歷史用電有功數(shù)據(jù),分辨率為15 min,共672個(gè)有功功率初值,相應(yīng)地選取該變壓器低壓側(cè)A相672個(gè)電壓值。由于用電數(shù)據(jù)在采集過(guò)程中有缺值,利用插值法對(duì)空值進(jìn)行估算補(bǔ)全。

    4.1 數(shù)據(jù)特征分析

    采用PCA對(duì)數(shù)據(jù)進(jìn)行降維處理。經(jīng)過(guò)多次數(shù)據(jù)仿真估算,驗(yàn)證了當(dāng)采用63維PCA特征數(shù)據(jù)可以表達(dá)原始數(shù)據(jù)95%的特征貢獻(xiàn)度,而且可以將聚類(lèi)分析聚類(lèi)結(jié)果偏差控制在1%以?xún)?nèi)。因此,選擇把原始數(shù)據(jù)維度降為63維度,這樣既可保持?jǐn)?shù)據(jù)原始特征,又可以減少聚類(lèi)整體計(jì)算時(shí)間。

    維度選取為

    (15)

    式中:Sk,fron為前k個(gè)特征值組成的面積;Sall為所有特征值組成的總面積;ηk,ctb為前k個(gè)特征值特征貢獻(xiàn)度。特征貢獻(xiàn)度如圖4所示。

    圖4 PCA的降維貢獻(xiàn)度

    4.2 用戶(hù)的分類(lèi)

    采用插值法補(bǔ)全原始數(shù)據(jù),共91戶(hù),672維特征值。采用模糊C-均值聚類(lèi)結(jié)果如圖5所示。由于用戶(hù)之間存在用電行為相似度較高,其數(shù)據(jù)特征不明顯,同一個(gè)用戶(hù)在3個(gè)聚類(lèi)中心隸屬度比較時(shí)容易陷入某兩個(gè)中心最優(yōu)解,如圖6所示。聚類(lèi)迭代進(jìn)行到第8代至第10代,目標(biāo)函數(shù)在第8代陷入早熟、收斂,如圖7所示。

    圖5 基于插值法的模糊聚類(lèi)結(jié)果

    圖6 基于插值法的模糊聚類(lèi)隸屬度

    圖7 基于插值法的模糊聚類(lèi)迭代次數(shù)

    采用原始數(shù)據(jù)進(jìn)行聚類(lèi)的方法只能得到兩類(lèi)用戶(hù),幾乎不能分辨出第3類(lèi)用戶(hù)。只有在用戶(hù)數(shù)據(jù)量足夠多才能分辨出第3類(lèi)用戶(hù),但此時(shí)存在數(shù)據(jù)量冗余度高、計(jì)算量大、計(jì)算時(shí)間長(zhǎng)等問(wèn)題。

    采用皮爾遜相關(guān)系數(shù)聚類(lèi)方法,對(duì)用戶(hù)用電與臺(tái)區(qū)電壓的關(guān)聯(lián)程度進(jìn)行聚類(lèi),能夠區(qū)分出用戶(hù)對(duì)于臺(tái)區(qū)電壓影響程度。

    但是,皮爾遜相關(guān)系數(shù)僅僅反映的是某一戶(hù)與臺(tái)區(qū)的關(guān)聯(lián)程度。由于數(shù)據(jù)具有高維特性,可能是某一些維度對(duì)此造成了較大的影響,并不能區(qū)分該用戶(hù)等級(jí),不能體現(xiàn)整體性,只能作為部分參考對(duì)象。皮爾遜相關(guān)系數(shù)仿真結(jié)果如圖8至圖9所示。

    圖8 基于皮爾遜法的模糊聚類(lèi)結(jié)果

    圖9 基于皮爾遜法的模糊聚類(lèi)迭代次數(shù)

    考慮到數(shù)據(jù)維數(shù)較高、數(shù)據(jù)特征不突出容易陷入局部最優(yōu)、聚類(lèi)計(jì)算量大等原因,運(yùn)用PCA法可以解決上述問(wèn)題。采用PCA法可以提取原數(shù)據(jù)95%特征貢獻(xiàn)度,用63維數(shù)據(jù)表示原始數(shù)據(jù)672維特征。

    仿真結(jié)果說(shuō)明,運(yùn)用PCA法,迭代次數(shù)在35次左右,收斂時(shí)間一般可以控制在0.568 0~0.589 1 s。對(duì)比加權(quán)后從小到大排序,采用PCA聚類(lèi)得到的用戶(hù)等級(jí),聚類(lèi)偏差在4%左右,對(duì)應(yīng)上面按照95%特征貢獻(xiàn)度,負(fù)荷誤差可接受范圍整體誤差在1%左右。PCA法仿真結(jié)果如圖10至圖12所示。

    圖10 基于PCA法的模糊聚類(lèi)結(jié)果

    圖11 基于PCA法的模糊聚類(lèi)隸屬度

    圖12 基于PCA法的模糊聚類(lèi)迭代次數(shù)

    采用PCA法,確定大、中、小等級(jí)用戶(hù)數(shù)目分別是28、33、30;相對(duì)誤差率約為1%。而采用皮爾遜相關(guān)系數(shù)法,確定大、中、小等級(jí)用戶(hù)的數(shù)目分別是33、15、43。

    相對(duì)誤差率計(jì)算式為

    (16)

    式中:ηre為聚類(lèi)用戶(hù)等級(jí)相對(duì)誤差率;ncls為聚類(lèi)類(lèi)簇包含用戶(hù)數(shù);nset為實(shí)際應(yīng)分類(lèi)簇用戶(hù)數(shù);nz為總的用戶(hù)數(shù);ηpca為PCA特征貢獻(xiàn)率。

    表2為不同聚類(lèi)方法的比較結(jié)果。從表2可見(jiàn):1)采用基于插值法的模糊聚類(lèi),因用戶(hù)之間數(shù)值相似度較高,發(fā)現(xiàn)某些粒子只能陷于某一個(gè)或者兩個(gè)聚類(lèi)中心之間,容易陷入早熟,聚類(lèi)效果不好,小數(shù)據(jù)容易陷入早熟,大量數(shù)據(jù)聚類(lèi)消耗時(shí)間長(zhǎng);2)采用基于皮爾遜相關(guān)系數(shù)的模糊聚類(lèi),聚類(lèi)花費(fèi)時(shí)間短,分類(lèi)效果好,但是只能反映相關(guān)程度,體現(xiàn)某一個(gè)粒子與聚類(lèi)中心的整體關(guān)聯(lián)程度,不能體現(xiàn)各個(gè)時(shí)間點(diǎn)聯(lián)系緊密程度,且對(duì)于用戶(hù)等級(jí)的分類(lèi)誤差率較高,皮爾遜特征值聚類(lèi)不適合用于用戶(hù)等級(jí)分類(lèi),只適用于用戶(hù)對(duì)于變壓器低壓側(cè)電壓影響度辨識(shí);3)采用基于PCA法的模糊聚類(lèi),提取特征值聚類(lèi)收斂時(shí)間短,分類(lèi)效果好,相對(duì)誤差率達(dá)1%,適合用于用戶(hù)等級(jí)分類(lèi)。

    表2 不同聚類(lèi)方法的比較

    綜上,在臺(tái)區(qū)電壓與用戶(hù)用電關(guān)系辨識(shí)中,結(jié)合PCA特征值聚類(lèi)法與皮爾遜相關(guān)系數(shù),能找出既具有等級(jí)分化且有較高影響因子的用戶(hù),此類(lèi)用戶(hù)即是供電管理部門(mén)需要重點(diǎn)關(guān)注的對(duì)象。

    4.3 用戶(hù)用電行為分析

    4.3.1 峰-谷-平荷期大用戶(hù)用電行為分析

    由用戶(hù)聚類(lèi)結(jié)果與皮爾遜相關(guān)系數(shù),從大用戶(hù)類(lèi)簇中選取第85號(hào)、第86號(hào)用戶(hù),影響因子分別為-0.215 3和0.107 6。該用戶(hù)的峰-谷-平負(fù)荷曲線(xiàn)如圖13所示,第85號(hào)用戶(hù)15:00 — 16:00處于用電低谷區(qū)間,此時(shí)變壓器低壓端電壓基本保持在238 V以上;16:15以后用電功率逐漸到達(dá)峰值,變壓器電壓逐漸下降到235 V,該用戶(hù)高峰區(qū)間一直持續(xù)到20:15,變壓器低壓側(cè)電壓保持235 V,電壓波動(dòng)范圍3~4 V;22:15 — 00:00用電功率處于平荷區(qū),此時(shí)變壓器電壓穩(wěn)定在238 V。綜合皮爾遜相關(guān)系數(shù)表與功率-電壓仿真圖,該用戶(hù)的用電行為與變壓器電壓波動(dòng)具有較強(qiáng)關(guān)系。大用戶(hù)用電負(fù)荷曲線(xiàn)如圖13所示。

    圖13 大用戶(hù)在峰-谷-平期用電負(fù)荷特性

    4.3.2 峰-谷-平荷期小用戶(hù)用電行為分析

    對(duì)于小用戶(hù),選取第26戶(hù)為分析對(duì)象,皮爾遜相關(guān)系數(shù)為0.14且為小用戶(hù)中相關(guān)系數(shù)最大的用戶(hù),對(duì)其進(jìn)行關(guān)聯(lián)性分析:結(jié)果顯示該用戶(hù)在17:15 —18:45處于用電低谷區(qū),對(duì)比同時(shí)間段的第85號(hào)大用戶(hù)用電功率正處于高峰,電壓處于較低235 V;該用戶(hù)在22:00 — 22:30處于用電高峰,但此區(qū)間變壓器電壓呈現(xiàn)上升趨勢(shì)。綜合皮爾遜系數(shù)及臺(tái)戶(hù)功率-電壓仿真圖,第26號(hào)小用戶(hù)并不如第85號(hào)大用戶(hù)與變壓器電壓波動(dòng)具有較強(qiáng)的相關(guān)性,小用戶(hù)用電負(fù)荷曲線(xiàn)如圖14所示。

    圖14 小用戶(hù)在峰-谷-平期用電負(fù)荷特性

    在小用戶(hù)類(lèi)簇當(dāng)中,依據(jù)皮爾遜相關(guān)系數(shù)選取用電行為相近的戶(hù)號(hào):5、7、18、26、27、37,這些都屬于該臺(tái)區(qū)下A相小用戶(hù)。其皮爾遜相關(guān)系數(shù)保持在0.1左右,重新計(jì)算以上用戶(hù)連成一片用電區(qū)域與變壓器電壓之間的皮爾遜系數(shù)為-0.203 2。由負(fù)荷曲線(xiàn)與電壓曲線(xiàn)圖15可知,在這一片用戶(hù)中,15:00—16:00保持在用電低谷區(qū),此時(shí)變壓器電壓保持240 V;18:00—18:45逐漸達(dá)到用電峰值,變壓器電壓同步下降到235 V;22:45—23:15為該片區(qū)域用電平荷期,變壓器電壓同時(shí)保持238 V。結(jié)合皮爾遜相關(guān)系數(shù)與該片區(qū)域峰-谷-平期有功功率、變壓器電壓的變化趨勢(shì),該區(qū)域用戶(hù)用電行為與變壓器電壓波動(dòng)具有較強(qiáng)的相關(guān)性。多用戶(hù)用電負(fù)荷曲線(xiàn)如圖15所示。

    圖15 多用戶(hù)在峰-谷-平期用電負(fù)荷特性

    5 結(jié) 語(yǔ)

    提出了一種基于PCA改進(jìn)模糊聚類(lèi)的臺(tái)區(qū)電壓與用戶(hù)用電關(guān)系辨識(shí)方法,并分析了用戶(hù)在用電峰-谷-平荷用電區(qū)間用電行為與臺(tái)區(qū)電壓之間的關(guān)系。通過(guò)對(duì)廣州某個(gè)小區(qū)連續(xù)7天的數(shù)據(jù)進(jìn)行仿真,結(jié)果表明:

    1)皮爾遜相關(guān)系數(shù)達(dá)到-0.2以上且聚類(lèi)為大用戶(hù),其用電高峰區(qū)間導(dǎo)致電壓下降1.282%,在該相下對(duì)于變壓器電壓波動(dòng)具有較強(qiáng)的相關(guān)性。

    2)對(duì)于單個(gè)小用戶(hù),其用電特征為:?jiǎn)螀^(qū)段用電功率小,用電間歇時(shí)間長(zhǎng),相關(guān)系數(shù)較低。但是,某一相下存在一片用電行為高度相似的小用戶(hù),該片區(qū)域在同一時(shí)間段相關(guān)系數(shù)達(dá)到-0.2以上,其用電規(guī)律同樣對(duì)變壓器電壓波動(dòng)造成一定的影響,需要引起重視。

    3)單個(gè)小用戶(hù)的用電特性:區(qū)間功率小,用電間歇長(zhǎng),在臺(tái)區(qū)電壓與用戶(hù)用電之間影響因子較弱。若非存在多個(gè)用電行為相似的皮爾遜系數(shù)用戶(hù),其用電行為與變壓器電壓波動(dòng)可忽略。

    猜你喜歡
    皮爾遜臺(tái)區(qū)特征值
    一類(lèi)帶強(qiáng)制位勢(shì)的p-Laplace特征值問(wèn)題
    單圈圖關(guān)聯(lián)矩陣的特征值
    現(xiàn)代統(tǒng)計(jì)學(xué)之父:卡爾·皮爾遜
    現(xiàn)代統(tǒng)計(jì)學(xué)之父:卡爾·皮爾遜
    Excel在水文學(xué)教學(xué)中的應(yīng)用
    卡方分布的探源
    降低臺(tái)區(qū)實(shí)時(shí)線(xiàn)損整治工作方法
    電子制作(2017年2期)2017-05-17 03:55:22
    基于商奇異值分解的一類(lèi)二次特征值反問(wèn)題
    三合一集中器作為臺(tái)區(qū)線(xiàn)損考核表計(jì)的探討
    關(guān)于兩個(gè)M-矩陣Hadamard積的特征值的新估計(jì)
    久久影院123| 午夜福利视频在线观看免费| 国产三级黄色录像| 欧美日韩视频精品一区| 中文字幕最新亚洲高清| 欧美激情久久久久久爽电影 | 国产1区2区3区精品| 亚洲欧美激情在线| 国产不卡av网站在线观看| 久久99一区二区三区| 免费在线观看亚洲国产| netflix在线观看网站| 一进一出抽搐gif免费好疼 | 视频区图区小说| 9色porny在线观看| 亚洲欧美日韩另类电影网站| 久久精品aⅴ一区二区三区四区| 亚洲精品中文字幕一二三四区| 免费在线观看黄色视频的| 亚洲情色 制服丝袜| 久热爱精品视频在线9| 国产xxxxx性猛交| 日韩视频一区二区在线观看| 国产成人欧美在线观看 | 夜夜爽天天搞| 男女高潮啪啪啪动态图| 亚洲专区国产一区二区| avwww免费| 免费观看人在逋| 亚洲人成伊人成综合网2020| 国产1区2区3区精品| 久久久久久人人人人人| 久久精品人人爽人人爽视色| 99热网站在线观看| 美国免费a级毛片| 亚洲av成人不卡在线观看播放网| 国产av精品麻豆| 三级毛片av免费| 丁香六月欧美| 村上凉子中文字幕在线| 久久久久久久精品吃奶| 高清在线国产一区| 欧美亚洲日本最大视频资源| 免费在线观看完整版高清| 中国美女看黄片| 亚洲国产中文字幕在线视频| 国产视频一区二区在线看| 精品人妻在线不人妻| cao死你这个sao货| 亚洲欧洲精品一区二区精品久久久| 欧美日韩精品网址| 国产91精品成人一区二区三区| 久久久久国内视频| 久久人妻av系列| 九色亚洲精品在线播放| 国产蜜桃级精品一区二区三区 | 成人特级黄色片久久久久久久| 热99久久久久精品小说推荐| 国产精华一区二区三区| 人妻丰满熟妇av一区二区三区 | 国产伦人伦偷精品视频| 亚洲成a人片在线一区二区| 美女福利国产在线| 日韩欧美国产一区二区入口| 免费女性裸体啪啪无遮挡网站| 国产区一区二久久| 亚洲一区中文字幕在线| 国产在线观看jvid| 最近最新免费中文字幕在线| 精品久久久久久久久久免费视频 | 欧美av亚洲av综合av国产av| 一级片免费观看大全| 美女扒开内裤让男人捅视频| 一级黄色大片毛片| 欧美黄色片欧美黄色片| 欧美乱妇无乱码| 黑人欧美特级aaaaaa片| 自线自在国产av| 女人久久www免费人成看片| av福利片在线| 日韩免费高清中文字幕av| 久久国产乱子伦精品免费另类| 丝袜美腿诱惑在线| 国产精品欧美亚洲77777| 女人被躁到高潮嗷嗷叫费观| 91麻豆精品激情在线观看国产 | 中文字幕最新亚洲高清| 国产精品久久久久久精品古装| 亚洲成av片中文字幕在线观看| 麻豆国产av国片精品| 精品国内亚洲2022精品成人 | 亚洲人成电影观看| 亚洲精品久久午夜乱码| 免费看a级黄色片| 欧美日韩亚洲综合一区二区三区_| 亚洲免费av在线视频| 亚洲欧美一区二区三区久久| 日韩 欧美 亚洲 中文字幕| 亚洲av熟女| 亚洲国产中文字幕在线视频| 国产亚洲av高清不卡| 精品乱码久久久久久99久播| 热99re8久久精品国产| 国产成人精品久久二区二区免费| 国产成人系列免费观看| 亚洲精品自拍成人| 久久 成人 亚洲| 亚洲精品在线观看二区| 最近最新免费中文字幕在线| 免费观看精品视频网站| 亚洲美女黄片视频| 免费看十八禁软件| 黄色毛片三级朝国网站| 女人高潮潮喷娇喘18禁视频| 男女之事视频高清在线观看| 国产精品欧美亚洲77777| 黄网站色视频无遮挡免费观看| 美女高潮喷水抽搐中文字幕| 性少妇av在线| 麻豆成人av在线观看| 黄片小视频在线播放| 麻豆国产av国片精品| 国产成人精品久久二区二区免费| 午夜福利乱码中文字幕| 麻豆国产av国片精品| 99久久99久久久精品蜜桃| 天天影视国产精品| 亚洲五月婷婷丁香| netflix在线观看网站| 日本欧美视频一区| 成年人黄色毛片网站| 亚洲精品av麻豆狂野| 99久久人妻综合| а√天堂www在线а√下载 | 国产有黄有色有爽视频| 亚洲欧洲精品一区二区精品久久久| 国产一区二区三区视频了| 亚洲精品国产精品久久久不卡| 国产麻豆69| 美国免费a级毛片| av在线播放免费不卡| 亚洲精品自拍成人| 欧美乱码精品一区二区三区| 他把我摸到了高潮在线观看| 色婷婷av一区二区三区视频| 亚洲精品成人av观看孕妇| 人人妻人人爽人人添夜夜欢视频| 久久久久久亚洲精品国产蜜桃av| 亚洲成人免费电影在线观看| 99精国产麻豆久久婷婷| 两个人看的免费小视频| 新久久久久国产一级毛片| 满18在线观看网站| 看免费av毛片| 亚洲欧美色中文字幕在线| 99国产精品一区二区蜜桃av | 久久久久久久午夜电影 | 久久国产亚洲av麻豆专区| 国产精品影院久久| 99国产精品一区二区蜜桃av | 叶爱在线成人免费视频播放| 亚洲国产中文字幕在线视频| 日韩欧美国产一区二区入口| 十八禁人妻一区二区| 中文字幕av电影在线播放| 女人高潮潮喷娇喘18禁视频| 成人手机av| 色播在线永久视频| 美女高潮喷水抽搐中文字幕| 人妻久久中文字幕网| 久久国产精品影院| 国产精品自产拍在线观看55亚洲 | 久久久久精品人妻al黑| 自拍欧美九色日韩亚洲蝌蚪91| 侵犯人妻中文字幕一二三四区| 天堂√8在线中文| 日本撒尿小便嘘嘘汇集6| 亚洲欧美精品综合一区二区三区| 十八禁网站免费在线| 国产亚洲欧美98| 午夜激情av网站| 99re6热这里在线精品视频| aaaaa片日本免费| 色老头精品视频在线观看| 中文字幕人妻丝袜一区二区| 成熟少妇高潮喷水视频| 成人国产一区最新在线观看| 在线观看午夜福利视频| 久久国产亚洲av麻豆专区| av在线播放免费不卡| 一区在线观看完整版| 美国免费a级毛片| 国产一卡二卡三卡精品| 国产区一区二久久| 欧美日韩亚洲高清精品| a级毛片在线看网站| 窝窝影院91人妻| 国产一区有黄有色的免费视频| 超色免费av| 99久久99久久久精品蜜桃| 在线永久观看黄色视频| 欧美黄色淫秽网站| 黄色 视频免费看| а√天堂www在线а√下载 | 亚洲欧美一区二区三区黑人| 欧美精品av麻豆av| 亚洲国产毛片av蜜桃av| 在线观看免费午夜福利视频| 亚洲五月色婷婷综合| av有码第一页| 国产精品免费大片| 国产在线一区二区三区精| 国产成人精品久久二区二区91| 久久青草综合色| 精品国产一区二区三区四区第35| 午夜精品久久久久久毛片777| 两性午夜刺激爽爽歪歪视频在线观看 | 久久国产精品男人的天堂亚洲| 十八禁人妻一区二区| 亚洲午夜理论影院| 亚洲成人国产一区在线观看| 欧美另类亚洲清纯唯美| 午夜视频精品福利| 男人的好看免费观看在线视频 | 精品久久久久久,| 一边摸一边做爽爽视频免费| 丝袜美腿诱惑在线| 一级毛片高清免费大全| 男女床上黄色一级片免费看| 黄色丝袜av网址大全| av线在线观看网站| 欧美精品一区二区免费开放| 又紧又爽又黄一区二区| 亚洲精品中文字幕一二三四区| 国产99久久九九免费精品| 久久精品熟女亚洲av麻豆精品| 亚洲五月天丁香| 91精品国产国语对白视频| 国产免费av片在线观看野外av| 日本撒尿小便嘘嘘汇集6| 一区二区三区激情视频| 免费人成视频x8x8入口观看| 美国免费a级毛片| 亚洲成av片中文字幕在线观看| 精品午夜福利视频在线观看一区| 久久婷婷成人综合色麻豆| 久久人妻福利社区极品人妻图片| 黄色毛片三级朝国网站| 人妻一区二区av| 欧美丝袜亚洲另类 | 国产亚洲精品一区二区www | 美女国产高潮福利片在线看| 黄色视频,在线免费观看| 最新美女视频免费是黄的| 亚洲美女黄片视频| 日韩一卡2卡3卡4卡2021年| 亚洲精品在线美女| 最新的欧美精品一区二区| 女人爽到高潮嗷嗷叫在线视频| 午夜福利,免费看| 亚洲综合色网址| 自线自在国产av| 嫁个100分男人电影在线观看| 中文字幕精品免费在线观看视频| 国产精品自产拍在线观看55亚洲 | 国产亚洲精品久久久久5区| 欧美最黄视频在线播放免费 | 亚洲精品粉嫩美女一区| 一级片'在线观看视频| 成人av一区二区三区在线看| 母亲3免费完整高清在线观看| 一区二区三区国产精品乱码| 一本一本久久a久久精品综合妖精| 国产亚洲欧美精品永久| 99久久人妻综合| 午夜精品在线福利| 国产精品一区二区精品视频观看| 亚洲国产精品一区二区三区在线| 天天躁日日躁夜夜躁夜夜| 老汉色av国产亚洲站长工具| 成人18禁高潮啪啪吃奶动态图| 国产真人三级小视频在线观看| 国产欧美日韩综合在线一区二区| 黄网站色视频无遮挡免费观看| 黄频高清免费视频| 国产成人精品久久二区二区免费| 亚洲中文av在线| 国产人伦9x9x在线观看| 国产精品久久久人人做人人爽| 在线av久久热| 十八禁高潮呻吟视频| 国产av又大| av免费在线观看网站| 99在线人妻在线中文字幕 | 99热只有精品国产| 正在播放国产对白刺激| 国产蜜桃级精品一区二区三区 | 天堂动漫精品| www.熟女人妻精品国产| 村上凉子中文字幕在线| 一进一出抽搐动态| 天天躁夜夜躁狠狠躁躁| 国产国语露脸激情在线看| 动漫黄色视频在线观看| 最新的欧美精品一区二区| 国产精品国产高清国产av | 一级黄色大片毛片| 午夜福利在线观看吧| 两性午夜刺激爽爽歪歪视频在线观看 | 51午夜福利影视在线观看| 国产成人免费无遮挡视频| 中文字幕人妻丝袜一区二区| 午夜免费鲁丝| 捣出白浆h1v1| 黄色怎么调成土黄色| 国产精品国产av在线观看| 色综合欧美亚洲国产小说| 乱人伦中国视频| 免费观看a级毛片全部| 女同久久另类99精品国产91| 色精品久久人妻99蜜桃| a级毛片黄视频| 欧美日韩亚洲综合一区二区三区_| 成人精品一区二区免费| 90打野战视频偷拍视频| netflix在线观看网站| 美女高潮到喷水免费观看| 国产在线精品亚洲第一网站| 国产免费av片在线观看野外av| 成人三级做爰电影| 啦啦啦免费观看视频1| 国产精品亚洲av一区麻豆| 久久人人97超碰香蕉20202| 久久午夜亚洲精品久久| 久久婷婷成人综合色麻豆| 91国产中文字幕| 黄色成人免费大全| 中国美女看黄片| 777久久人妻少妇嫩草av网站| 99国产精品一区二区蜜桃av | 91国产中文字幕| 一边摸一边做爽爽视频免费| 91国产中文字幕| 亚洲三区欧美一区| 亚洲av欧美aⅴ国产| 免费在线观看视频国产中文字幕亚洲| а√天堂www在线а√下载 | 久久精品亚洲av国产电影网| 黄网站色视频无遮挡免费观看| 一边摸一边抽搐一进一小说 | 麻豆国产av国片精品| 热re99久久精品国产66热6| 一级片免费观看大全| 757午夜福利合集在线观看| 久久中文字幕一级| 超碰成人久久| 俄罗斯特黄特色一大片| 免费在线观看视频国产中文字幕亚洲| 叶爱在线成人免费视频播放| 欧美性长视频在线观看| 国产精品久久久久久人妻精品电影| 久久精品亚洲精品国产色婷小说| 国产亚洲精品久久久久5区| 日韩免费高清中文字幕av| 国内毛片毛片毛片毛片毛片| 一本大道久久a久久精品| 欧美日韩乱码在线| 热99国产精品久久久久久7| 久久人妻福利社区极品人妻图片| 国产精品乱码一区二三区的特点 | 国产欧美日韩精品亚洲av| 国产高清videossex| av视频免费观看在线观看| 成人18禁在线播放| 久久中文看片网| 成人影院久久| 亚洲视频免费观看视频| 成人影院久久| 一二三四社区在线视频社区8| 免费观看精品视频网站| 久久久国产精品麻豆| 国产成人精品久久二区二区免费| 久久人妻福利社区极品人妻图片| 免费在线观看视频国产中文字幕亚洲| 最近最新中文字幕大全免费视频| 下体分泌物呈黄色| 国产精品久久久久久人妻精品电影| 欧美日韩精品网址| 欧美日韩成人在线一区二区| 成人av一区二区三区在线看| x7x7x7水蜜桃| 91成人精品电影| 久久国产亚洲av麻豆专区| 女性生殖器流出的白浆| 精品久久久精品久久久| 久久99一区二区三区| 老熟妇乱子伦视频在线观看| 1024视频免费在线观看| 黄网站色视频无遮挡免费观看| 午夜福利在线观看吧| 一区二区三区国产精品乱码| 精品久久久久久久毛片微露脸| 欧美乱码精品一区二区三区| 满18在线观看网站| 亚洲精品成人av观看孕妇| 在线观看免费日韩欧美大片| 黄色 视频免费看| 日本一区二区免费在线视频| 十八禁高潮呻吟视频| 女人被狂操c到高潮| 久久国产乱子伦精品免费另类| 777米奇影视久久| 亚洲片人在线观看| 搡老乐熟女国产| 他把我摸到了高潮在线观看| 怎么达到女性高潮| videosex国产| 久久精品亚洲熟妇少妇任你| 夜夜爽天天搞| 国产99白浆流出| 99热国产这里只有精品6| 天天躁夜夜躁狠狠躁躁| 欧美黄色淫秽网站| 国产主播在线观看一区二区| 淫妇啪啪啪对白视频| 老汉色av国产亚洲站长工具| 欧美日韩国产mv在线观看视频| 性色av乱码一区二区三区2| 丰满人妻熟妇乱又伦精品不卡| 在线永久观看黄色视频| 一级毛片女人18水好多| 99精国产麻豆久久婷婷| www日本在线高清视频| 亚洲aⅴ乱码一区二区在线播放 | 欧美黑人欧美精品刺激| 午夜日韩欧美国产| 免费女性裸体啪啪无遮挡网站| 亚洲午夜精品一区,二区,三区| 亚洲一区二区三区欧美精品| 国产高清视频在线播放一区| 欧美黄色淫秽网站| 午夜福利一区二区在线看| 国产aⅴ精品一区二区三区波| 老汉色av国产亚洲站长工具| 精品久久久精品久久久| av超薄肉色丝袜交足视频| 一区二区三区激情视频| 日本一区二区免费在线视频| 国产在线观看jvid| av天堂久久9| 亚洲熟妇中文字幕五十中出 | 精品亚洲成国产av| 国产精品久久电影中文字幕 | 国产精品 欧美亚洲| 国产不卡av网站在线观看| 淫妇啪啪啪对白视频| 国产成人精品在线电影| 亚洲午夜精品一区,二区,三区| 欧美亚洲 丝袜 人妻 在线| 亚洲精品中文字幕在线视频| x7x7x7水蜜桃| 国产xxxxx性猛交| 免费av中文字幕在线| 欧美黄色片欧美黄色片| 亚洲av成人一区二区三| 国产成人精品久久二区二区91| 中文亚洲av片在线观看爽 | www.自偷自拍.com| 狂野欧美激情性xxxx| 一区二区日韩欧美中文字幕| 国产99久久九九免费精品| 亚洲九九香蕉| 国产精品98久久久久久宅男小说| 大香蕉久久成人网| 亚洲精品在线观看二区| 午夜成年电影在线免费观看| 亚洲成av片中文字幕在线观看| 精品高清国产在线一区| 欧美在线一区亚洲| 精品国产亚洲在线| 久久这里只有精品19| 亚洲国产看品久久| 怎么达到女性高潮| 久久中文字幕人妻熟女| 精品少妇一区二区三区视频日本电影| bbb黄色大片| 国产野战对白在线观看| 国产麻豆69| 日韩欧美国产一区二区入口| 国产精品美女特级片免费视频播放器 | 一级黄色大片毛片| 男女免费视频国产| 99国产精品一区二区三区| 国产av精品麻豆| 亚洲av成人一区二区三| 女性生殖器流出的白浆| 国产高清国产精品国产三级| 久久精品国产a三级三级三级| 母亲3免费完整高清在线观看| 看片在线看免费视频| www.999成人在线观看| 亚洲avbb在线观看| 国产深夜福利视频在线观看| 一级作爱视频免费观看| 亚洲成人免费电影在线观看| 中文字幕人妻丝袜制服| 王馨瑶露胸无遮挡在线观看| 一二三四社区在线视频社区8| xxx96com| 久久天堂一区二区三区四区| 精品欧美一区二区三区在线| 国产一区有黄有色的免费视频| 精品无人区乱码1区二区| 国产视频一区二区在线看| 18禁美女被吸乳视频| 亚洲自偷自拍图片 自拍| 中文字幕人妻丝袜制服| 好看av亚洲va欧美ⅴa在| 免费高清在线观看日韩| 亚洲午夜理论影院| 午夜福利在线免费观看网站| 热re99久久国产66热| 色94色欧美一区二区| 久久国产亚洲av麻豆专区| 捣出白浆h1v1| 两个人看的免费小视频| 两个人免费观看高清视频| 在线av久久热| 久久久久久久国产电影| 久久久久国内视频| 在线观看66精品国产| 在线看a的网站| 一本综合久久免费| 天天添夜夜摸| 啦啦啦在线免费观看视频4| 少妇被粗大的猛进出69影院| 可以免费在线观看a视频的电影网站| 亚洲伊人色综图| 国产欧美日韩综合在线一区二区| 身体一侧抽搐| 黄片大片在线免费观看| 国产99白浆流出| tocl精华| 身体一侧抽搐| 一进一出抽搐动态| 日本黄色视频三级网站网址 | netflix在线观看网站| 久久久国产精品麻豆| 丝袜美足系列| 亚洲第一青青草原| 国产男女内射视频| 一a级毛片在线观看| 久久精品国产99精品国产亚洲性色 | 亚洲成人免费电影在线观看| 亚洲国产精品一区二区三区在线| 亚洲av片天天在线观看| 国产一区二区激情短视频| 久久青草综合色| 国产免费现黄频在线看| 亚洲在线自拍视频| 人妻一区二区av| 国产亚洲欧美精品永久| av欧美777| 欧美中文综合在线视频| 午夜福利一区二区在线看| 午夜免费鲁丝| 在线观看一区二区三区激情| 精品午夜福利视频在线观看一区| 日韩成人在线观看一区二区三区| 日韩欧美在线二视频 | e午夜精品久久久久久久| 日韩精品免费视频一区二区三区| 国产亚洲欧美精品永久| 国产高清激情床上av| 美女扒开内裤让男人捅视频| 久久精品国产亚洲av香蕉五月 | 91麻豆精品激情在线观看国产 | 亚洲色图综合在线观看| 成在线人永久免费视频| 午夜福利影视在线免费观看| 夜夜爽天天搞| 亚洲精品在线美女| 欧美色视频一区免费| 伊人久久大香线蕉亚洲五| 日本wwww免费看| 日韩 欧美 亚洲 中文字幕| 丰满人妻熟妇乱又伦精品不卡| 亚洲av美国av| 真人做人爱边吃奶动态| 国内毛片毛片毛片毛片毛片| 女性被躁到高潮视频| 国产av精品麻豆| 成人18禁高潮啪啪吃奶动态图| 91麻豆精品激情在线观看国产 | 亚洲精品国产色婷婷电影| 99精品欧美一区二区三区四区| 看免费av毛片| 久久久国产精品麻豆| 日韩免费高清中文字幕av| 欧美在线一区亚洲| 麻豆乱淫一区二区| 亚洲五月婷婷丁香| 国产精品综合久久久久久久免费 | 王馨瑶露胸无遮挡在线观看| 国产精品一区二区在线观看99| 韩国av一区二区三区四区| 满18在线观看网站| 下体分泌物呈黄色| 男人操女人黄网站| 欧美+亚洲+日韩+国产| 精品欧美一区二区三区在线| 午夜免费观看网址| 美女福利国产在线| 欧美日韩亚洲综合一区二区三区_| 成年版毛片免费区|