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

    考慮辨識(shí)結(jié)果連續(xù)性的Type-Ⅲ型工商業(yè)負(fù)荷辨識(shí)方法

    2021-12-29 07:27:00張振宇李巍巍
    電力系統(tǒng)自動(dòng)化 2021年24期
    關(guān)鍵詞:互信息特征選擇準(zhǔn)確率

    段 晶,李 勇,張振宇,李巍巍,蔣 林,李 磊

    (1. 湖南大學(xué)電氣與信息工程學(xué)院,湖南省長(zhǎng)沙市 410082;2. 長(zhǎng)沙新奧長(zhǎng)燃能源發(fā)展有限公司,湖南省長(zhǎng)沙市 410005;3. 利物浦大學(xué)電氣工程與電子系,利物浦L69 3BX,英國(guó);4. 長(zhǎng)沙新奧湘江新能源發(fā)展有限公司,湖南省長(zhǎng)沙市 410005)

    0 引言

    為提高電力用戶終端負(fù)荷的能源利用效率、減少電量消耗、促進(jìn)用戶參與需求側(cè)管理或?qū)崿F(xiàn)遠(yuǎn)程優(yōu)化控制,電力負(fù)荷的在線監(jiān)測(cè)必不可少[1-2]。電力負(fù)荷的在線監(jiān)測(cè)方法主要有侵入式負(fù)荷監(jiān)測(cè)(ILM)和非侵入式負(fù)荷監(jiān)測(cè)(NILM)兩種。ILM 需要為每個(gè)監(jiān)測(cè)設(shè)備配置一個(gè)或多個(gè)傳感器,雖然監(jiān)測(cè)結(jié)果更為準(zhǔn)確,但成本高、配置復(fù)雜、維護(hù)困難。NILM的概念最早由Hart 于1992 年提出[3],它只需在被監(jiān)測(cè)系統(tǒng)的電力入口處安裝一個(gè)電力信息采集裝置,再通過信號(hào)分析和處理實(shí)現(xiàn)負(fù)荷的分解與辨識(shí),雖然精度稍低,但是硬件結(jié)構(gòu)簡(jiǎn)單,便于安裝和維護(hù),更切合智能電網(wǎng)建設(shè)的實(shí)際需求[4-5]。因此,近年來NILM 方法的研究備受國(guó)內(nèi)外學(xué)者的關(guān)注。

    負(fù)荷辨識(shí)是負(fù)荷監(jiān)測(cè)的核心內(nèi)容,負(fù)荷辨識(shí)方法可分為三大類:優(yōu)化方法、監(jiān)督學(xué)習(xí)方法和無(wú)監(jiān)督學(xué)習(xí)方法。優(yōu)化方法是負(fù)荷辨識(shí)最直接的方法,將提取的負(fù)荷特征數(shù)據(jù)與模板庫(kù)中的特征數(shù)據(jù)進(jìn)行對(duì)比,采用優(yōu)化算法最小化誤差,找到匹配度最高的電氣設(shè)備狀態(tài)組合,包括動(dòng)態(tài)時(shí)間彎曲(DTW)算法[6]、遺傳算法[7]、粒子群優(yōu)化算法[8]等,但優(yōu)化方法進(jìn)行負(fù)荷辨識(shí)普遍需要根據(jù)歷史信息提前建立負(fù)荷特征模板庫(kù)。監(jiān)督學(xué)習(xí)方法需要已知電氣設(shè)備的特定運(yùn)行信息進(jìn)行模型的初始訓(xùn)練,且算法的性能很大程度取決于模型學(xué)習(xí)的數(shù)據(jù)標(biāo)簽,基本上是將負(fù)荷辨識(shí)問題轉(zhuǎn)換為分類學(xué)習(xí)問題進(jìn)行解決,主要方法包括貝葉斯分類器[9]、支持向量機(jī)(SVM)[10]和人工神經(jīng)網(wǎng)絡(luò)(ANN)[11-13]。無(wú)監(jiān)督學(xué)習(xí)方法可以在沒有任何先驗(yàn)信息的情況下運(yùn)行,相較于監(jiān)督學(xué)習(xí)方法減少了訓(xùn)練步驟的干擾性,主要算法包括隱馬爾可夫模型(HMM)[14]及其擴(kuò)展,如文獻(xiàn)[15]中的因子隱馬爾可夫模型(FHMM)和文獻(xiàn)[16]中的隱式半馬爾可夫模型(HSMM)均被應(yīng)用到負(fù)荷辨識(shí)中。

    不同類型的負(fù)荷具有不同的負(fù)荷特征,按照其運(yùn)行狀態(tài)可大致分為3 類[17]。對(duì)于連續(xù)變化(Type-Ⅲ)型的負(fù)荷,由于其工作狀態(tài)的功率值保持連續(xù)變化的特性,導(dǎo)致其狀態(tài)空間極其巨大。文獻(xiàn)[18]指出,對(duì)于常用的基于HMM 的負(fù)荷辨識(shí)方法,僅對(duì)Type-Ⅰ型的開關(guān)負(fù)荷表現(xiàn)較好,在面對(duì)狀態(tài)更為復(fù)雜的其他類型負(fù)荷時(shí)效果將被極大的削弱。因此Type-Ⅲ型負(fù)荷的辨識(shí)一直是NILM 中研究的重點(diǎn)和難點(diǎn)之一。

    為有效辨識(shí)Type-Ⅲ型負(fù)荷,滿足負(fù)荷辨識(shí)算法對(duì)辨識(shí)準(zhǔn)確率和算法魯棒性的要求,本文提出了一種基于深度卷積神經(jīng)網(wǎng)絡(luò)(CNN)和HMM 的非侵入式負(fù)荷辨識(shí)算法。通過深度CNN 提取負(fù)荷的多維特征,進(jìn)行負(fù)荷的初次辨識(shí),將連續(xù)變化的負(fù)荷狀態(tài)辨識(shí)問題重構(gòu)為二分類狀態(tài)辨識(shí);根據(jù)負(fù)荷狀態(tài)的時(shí)間依賴性,利用HMM 解決負(fù)荷的二分類狀態(tài)再辨識(shí),減少CNN 辨識(shí)結(jié)果中出現(xiàn)的狀態(tài)斷點(diǎn)現(xiàn)象,提高辨識(shí)準(zhǔn)確率。本文以工商業(yè)電力負(fù)荷為例,在具有代表性的Type-Ⅲ型負(fù)荷數(shù)據(jù)上進(jìn)行訓(xùn)練和測(cè)試,證明了本文算法的有效性,并且可達(dá)到97.9%的辨識(shí)準(zhǔn)確率。

    1 負(fù)荷數(shù)據(jù)處理與特征選擇

    附錄A 圖A1 展示了3 種典型的負(fù)荷曲線,其中Type-Ⅲ型負(fù)荷的狀態(tài)連續(xù)可變性特點(diǎn)非常明顯。工商業(yè)用電在總負(fù)荷中占有很大的比重,研究工商業(yè)用電特性,對(duì)電力系統(tǒng)安全經(jīng)濟(jì)運(yùn)行和電力市場(chǎng)的發(fā)展有著重要的作用。對(duì)工商業(yè)電力負(fù)荷進(jìn)行監(jiān)測(cè)和辨識(shí),可以幫助電力部門了解客戶的用電行為,進(jìn)而為不同的客戶群制訂相應(yīng)的市場(chǎng)策略,使雙方都受益。因此,本文使用Type-Ⅲ型工商業(yè)負(fù)荷來驗(yàn)證非侵入式負(fù)荷辨識(shí)算法的性能。

    1.1 數(shù)據(jù)處理

    本文所用的數(shù)據(jù)集為實(shí)測(cè)電力數(shù)據(jù)。測(cè)量?jī)x器為日置電能質(zhì)量分析儀PQ3198,可完成無(wú)間隙的連續(xù)運(yùn)算、事件檢測(cè)、GPS 時(shí)間同步等標(biāo)準(zhǔn)要求的測(cè)量,測(cè)量項(xiàng)目包括有功功率、無(wú)功功率、頻率、電壓、電流及其有效值、最值和高次諧波成分等多維負(fù)荷特征數(shù)據(jù),并具備波形記錄和存儲(chǔ)功能。在電力入口處安裝數(shù)據(jù)采集裝備,采樣頻率為1 Hz,采樣時(shí)長(zhǎng)為7 d。為了驗(yàn)證辨識(shí)效果,在待辨識(shí)負(fù)荷處也安裝采集設(shè)備以準(zhǔn)確獲取其運(yùn)行狀態(tài),從而為辨識(shí)模型提供參考數(shù)據(jù)。

    對(duì)于電力入口處采集的任一負(fù)荷特征數(shù)據(jù),可用時(shí)間序列R=[R0,R1,…,RT]表示,其中T為負(fù)荷數(shù)據(jù)的時(shí)間長(zhǎng)度。數(shù)據(jù)處理中可通過差分的方式來增強(qiáng)特征的表達(dá)能力,因此可以得到負(fù)荷的差值特征序列x=[x1,x2,…,xT],其中任一元素xt皆滿足:

    給定時(shí)間窗口長(zhǎng)度w,并通過特征選擇確定n種負(fù)荷特征,由此可生成CNN 模型能學(xué)習(xí)的輸入負(fù)荷數(shù)據(jù){H1,H2,…,HL},其中L為輸入數(shù)據(jù)的總長(zhǎng)度。

    式中:Xi,j為第j種負(fù)荷特征的第i個(gè)窗口的差值時(shí)間序列,如式(3)所示。

    式中:i=1,2,…,T-w+1;j=1,2,…,n。

    對(duì)于采集的待辨識(shí)負(fù)荷數(shù)據(jù),以其電流有效值Ii+w/2的大小為標(biāo)準(zhǔn),生成設(shè)備的標(biāo)簽狀態(tài)集{y1,y2,…,yL},其中yi為第i個(gè)時(shí)間窗口內(nèi)的i+w/2 時(shí)刻點(diǎn)的負(fù)荷的運(yùn)行狀態(tài),i=1,2,…,L。

    其中,將1 A 作為狀態(tài)標(biāo)記的基準(zhǔn)是根據(jù)待辨識(shí)負(fù)荷運(yùn)行時(shí)的最小電流值所選擇的。

    1.2 特征選擇

    由電能質(zhì)量分析儀采集的電力數(shù)據(jù)共包含3 537 種負(fù)荷特征,其中包括能反映待辨識(shí)負(fù)荷運(yùn)行狀態(tài)變化的相關(guān)特征和與待辨識(shí)負(fù)荷狀態(tài)變化不相關(guān)的無(wú)關(guān)特征。除此之外,負(fù)荷特征中還包括一類冗余特征,即這類特征中所包含的信息可從其他特征中提?。?9]。冗余特征在負(fù)荷辨識(shí)模型學(xué)習(xí)中幾乎不起作用,而無(wú)關(guān)特征反而會(huì)對(duì)模型學(xué)習(xí)造成干擾。因此,在進(jìn)行模型訓(xùn)練之前要進(jìn)行特征選擇,去除無(wú)關(guān)特征和冗余特征,從而提高負(fù)荷辨識(shí)的準(zhǔn)確率,并且特征選擇可以在很大程度上降低維數(shù),減少了模型的運(yùn)算量。

    時(shí)間窗口長(zhǎng)度、特征維度與CNN 的乘-加操作(multiply-accumulate operations,MAC)次數(shù)之間的關(guān)系如附錄A 圖A2 所示。隨著特征數(shù)目的增加,神經(jīng)網(wǎng)絡(luò)的操作次數(shù)逐漸變得不可接受,另外,由于特征之間的組合狀態(tài)爆炸增長(zhǎng),導(dǎo)致需要大量的數(shù)據(jù)才能對(duì)網(wǎng)絡(luò)進(jìn)行合適的訓(xùn)練。因此,為了降低網(wǎng)絡(luò)運(yùn)算復(fù)雜度、縮短模型訓(xùn)練時(shí)間,本文采用互信息算法進(jìn)行特征選擇。

    在概率論和信息論中,2 個(gè)隨機(jī)變量的互信息是變量間相互依賴性的量度,互信息熵越大表明2 個(gè)隨機(jī)變量之間的相關(guān)性越強(qiáng),反之互信息熵越小,2 個(gè)隨機(jī)變量之間的相關(guān)性就越弱[20]?;バ畔⒅礗(x,y)的計(jì)算如下:

    式中:p(x)和p(y)分別為隨機(jī)變量x和y的邊緣概率密度函數(shù);p(x,y)為x和y的聯(lián)合概率密度函數(shù)。

    本文采用互信息進(jìn)行特征選擇的具體步驟為:根據(jù)式(5)計(jì)算每個(gè)負(fù)荷特征與待辨識(shí)負(fù)荷運(yùn)行狀態(tài)之間的互信息值;按照互信息的大小對(duì)所有特征按從大到小進(jìn)行排序;然后選擇前n個(gè)互信息最大的特征,構(gòu)成負(fù)荷辨識(shí)的初始特征集合{f1,f2,…,fn};計(jì)算fi與fj(i<j)之間的互信息大小,若fi與某個(gè)特征fj的互信息值大于閾值,則將fj移除特征子集,并按照順序添加新的特征;重新計(jì)算新的特征子集中特征之間的互信息值,重復(fù)篩選直到兩兩特征之間的互信息值皆小于閾值,并構(gòu)建新的特征集合。

    附錄A 圖A3 以熱圖的形式展示了負(fù)荷特征之間的相關(guān)性。圖A3(a)中互信息值很小,不同特征之間的相關(guān)性很弱,而圖A3(b)則相反,說明不同負(fù)荷特征之間存在冗余性,特征提取時(shí)可從不同的特征中提取重疊的信息,以致于增加了許多不必要的運(yùn)算。因此,去除冗余特征可使得模型在相同的運(yùn)算維度內(nèi),提取到更多的隱含特征。

    2 深度CNN-HMM 負(fù)荷辨識(shí)算法原理

    對(duì)于Type-Ⅲ型負(fù)荷,本文引入了具備強(qiáng)特征提取能力的CNN 模型,用于提取連續(xù)的負(fù)荷功率變化中隱含的特征。但是,純CNN 類的方法具有另一種弊端,在本文中將其表述為辨識(shí)結(jié)果的非連續(xù)性問題,如附錄A 圖A4 所示。由于CNN 本身不具備時(shí)序建模能力,所以其對(duì)時(shí)序數(shù)據(jù)的感知范圍僅能停留在所選取的時(shí)間窗長(zhǎng)度內(nèi),這導(dǎo)致CNN 的辨識(shí)結(jié)果中經(jīng)常出現(xiàn)局部性的辨識(shí)錯(cuò)誤。值得注意的是,使用CNN 重新提取負(fù)荷數(shù)據(jù)中的特征并將其重構(gòu)為運(yùn)行或非運(yùn)行的二分類狀態(tài),可以解決HMM面臨的狀態(tài)空間過大的問題,使其重新適用于HMM 對(duì)負(fù)荷狀態(tài)的表達(dá)。因此,本文最終確立了使用CNN-HMM 架構(gòu)的Type-Ⅲ型負(fù)荷運(yùn)行狀態(tài)辨識(shí)的解決方案。

    2.1 算法流程

    基于深度CNN-HMM 的非侵入式負(fù)荷辨識(shí)算法流程如圖1 所示。

    圖1 基于深度CNN-HMM 的負(fù)荷辨識(shí)算法流程圖Fig.1 Flow chart of load identification algorithm based on deep CNN-HMM

    第1 階段為數(shù)據(jù)處理:采用互信息算法對(duì)電力入口處采集的高維負(fù)荷數(shù)據(jù)進(jìn)行特征選擇,并采用差分處理的方式強(qiáng)化特征表達(dá)能力;對(duì)于待辨識(shí)負(fù)荷數(shù)據(jù),根據(jù)電流有效值大小確定其運(yùn)行狀態(tài),生成設(shè)備的標(biāo)簽數(shù)據(jù)集。

    第2 階段為模型訓(xùn)練:將負(fù)荷特征數(shù)據(jù)和標(biāo)簽數(shù)據(jù)合成模型的訓(xùn)練數(shù)據(jù)集,訓(xùn)練深度CNN,使其學(xué)習(xí)擬合特征數(shù)據(jù)與標(biāo)簽之間的非線性映射關(guān)系;通過深度CNN 得到負(fù)荷的初始辨識(shí)結(jié)果,并將其作為HMM 的已知觀測(cè)序列,訓(xùn)練HMM 的參數(shù)。

    第3 階段為模型評(píng)估:將負(fù)荷特征數(shù)據(jù)和標(biāo)簽數(shù)據(jù)合成模型的測(cè)試數(shù)據(jù)集,采用深度CNN 模型提取負(fù)荷的隱含特征并進(jìn)行初次辨識(shí),然后采用HMM 優(yōu)化CNN 辨識(shí)結(jié)果中的狀態(tài)斷點(diǎn)問題,由此得到負(fù)荷的最終狀態(tài)序列,并根據(jù)負(fù)荷標(biāo)簽數(shù)據(jù)對(duì)CNN-HMM 辨識(shí)結(jié)果的準(zhǔn)確率進(jìn)行評(píng)估。

    2.2 深度CNN 結(jié)構(gòu)

    殘差神經(jīng)網(wǎng)絡(luò)(ResNet)因其特殊的捷徑連接結(jié)構(gòu),不僅能通過增加網(wǎng)絡(luò)的深度來提高準(zhǔn)確率,而且可緩解由深度增加所帶來的梯度消失問題[21-22]。因此,本文采用ResNet-50 網(wǎng)絡(luò)作為深度CNN 的基本架構(gòu),通過ResNet-50 的學(xué)習(xí),逐漸增加輸入特征向量的深度,減小其數(shù)據(jù)維度,提取數(shù)據(jù)的深層隱含特征,由此優(yōu)化神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)性能。附錄A 表A1為ResNet-50 網(wǎng)絡(luò)結(jié)構(gòu)的超參數(shù)。

    2.3 HMM 參數(shù)計(jì)算

    HMM 是關(guān)于時(shí)間序列的概率模型,可用一個(gè)五元組參數(shù)λ={U,V,A,B,π}來表示,其參數(shù)定義如下。

    1)U={u1,u2,…,uN}是一組隱含狀態(tài)的集合,其中N為隱含狀態(tài)數(shù)目,系統(tǒng)在t時(shí)刻的隱含狀態(tài)Qt∈U。本文采用CNN 進(jìn)行負(fù)荷特征提取,并將其重構(gòu)為二分類狀態(tài)負(fù)荷辨識(shí)問題,所以本文實(shí)驗(yàn)所選擇的待辨識(shí)負(fù)荷的隱含狀態(tài)數(shù)目N=2,只有工作狀態(tài)和非工作狀態(tài),即U={0,1}。

    2)V={v1,v2,…,vM}是一組觀測(cè)狀態(tài)的集合,其中M為每一個(gè)隱含狀態(tài)可能輸出的不同觀測(cè)值的數(shù)目,系統(tǒng)在t時(shí)刻的觀測(cè)值St∈V。對(duì)于待辨識(shí)負(fù)荷,觀測(cè)狀態(tài)數(shù)目M=2,即V={0,1}。

    3)狀態(tài)轉(zhuǎn)移概率矩陣A=(aij),其中aij為在t時(shí)刻狀態(tài)為ui的條件下,在t+1 時(shí)刻狀態(tài)為uj的概率,計(jì)算公式如下。

    4)觀測(cè)狀態(tài)轉(zhuǎn)移概率矩陣B=(bi(k)),其中bi(k)為在t時(shí)刻隱含狀態(tài)為ui的條件下,觀測(cè)狀態(tài)為vk的概率。

    5)初始狀態(tài)概率矩陣π=(πi)為t=1 時(shí)刻初始狀態(tài)的概率分布。

    本文提出的非侵入式負(fù)荷辨識(shí)算法采用ResNet-50 模型對(duì)負(fù)荷數(shù)據(jù)進(jìn)行初辨識(shí),并將其重構(gòu)為運(yùn)行或非運(yùn)行的二分類問題,再采用HMM 對(duì)初始辨識(shí)結(jié)果進(jìn)行優(yōu)化。HMM 算法的具體步驟為:①初始化HMM 可更新參數(shù)λ′0={A,B,π};②根據(jù)深度CNN 的初辨識(shí)結(jié)果生成觀測(cè)序列S=[S1,S2,…,SL];③迭代更新的參數(shù)λ′,直到觀測(cè)序列出現(xiàn)的概率P(S|λ′)達(dá)最大值;④通過最優(yōu)參數(shù)λ′和觀測(cè)序列S,采用維特比算法求使得可能性P(Q|S,λ′)最大的狀態(tài)序列Q=[Q1,Q2,…,QL],即為待辨識(shí)負(fù)荷的最終辨識(shí)結(jié)果。

    3 算例分析

    3.1 評(píng)價(jià)指標(biāo)

    負(fù)荷辨識(shí)算法的性能可用準(zhǔn)確率Acc、精確度P′、召回率R以及調(diào)和平均評(píng)價(jià)指標(biāo)Fscore進(jìn)行評(píng)估。準(zhǔn)確率Acc的計(jì)算如下:

    式中:m為經(jīng)過算法辨識(shí)預(yù)測(cè)正確的數(shù)量。

    精確度P′用來表示算法抗誤檢的能力,召回率R用來表示算法抗漏檢的能力。

    式中:TP為被辨識(shí)為正樣本、實(shí)際為正樣本的數(shù)量;FP為被辨識(shí)為正樣本、實(shí)際為負(fù)樣本的數(shù)量;FN為被辨識(shí)為負(fù)樣本、實(shí)際為正樣本的數(shù)量。

    Fscore是用來平衡精確度和召回率的評(píng)價(jià)指標(biāo)。

    此外,本文采用Wasserstein 距離衡量負(fù)荷辨識(shí)結(jié)果的運(yùn)行時(shí)間概率分布與實(shí)際負(fù)荷運(yùn)行時(shí)間概率分布之間的差異[23],距離可表征算法解決識(shí)別錯(cuò)誤斷點(diǎn)的能力,該值越小表明2 個(gè)概率分布之間的擬合程度越好。Wasserstein 距離的計(jì)算公式如下。

    式中:pr和pg為需要計(jì)算二者之間距離的2 個(gè)概率分布;inf(·)表示取最小值;Π(pr,pg)為以pr和pg為邊緣分布的所有可能的聯(lián)合概率分布γ的集合。從γ中采樣(x,y)~γ,計(jì)算其范數(shù)距離||x-y||,由此計(jì)算該聯(lián)合分布下的距離期望值E(x,y)~γ(||x-y||),所有可能的聯(lián)合分布中取該期望值的下確界為Wasserstein 距離。

    3.2 算例數(shù)據(jù)集介紹

    本文實(shí)驗(yàn)所用數(shù)據(jù)集為實(shí)測(cè)電力數(shù)據(jù),數(shù)據(jù)采集地點(diǎn)分別為湖南大學(xué)研究生教學(xué)樓和湘江歡樂城能源站。其中,教學(xué)樓背景負(fù)荷較為簡(jiǎn)單,主要用能設(shè)備包括電梯與中央空調(diào)。而能源站用能結(jié)構(gòu)較為復(fù)雜,它是一個(gè)冷熱電綜合能源站,包括溴化鋰機(jī)組、蒸汽鍋爐、燃?xì)鉄崴仩t、熱水鍋爐循環(huán)泵等多種用能設(shè)備以及內(nèi)燃發(fā)電機(jī)等發(fā)電設(shè)備。除去大型設(shè)備外,也存在數(shù)量眾多的中小型負(fù)荷,由于供能覆蓋范圍較大,難以窮舉。為滿足2 所酒店與冰雪世界的用能需求,設(shè)備機(jī)組經(jīng)常處于復(fù)雜的動(dòng)態(tài)調(diào)控過程中,用電情況多變、辨識(shí)難度較高,是一個(gè)典型的工商業(yè)復(fù)雜負(fù)荷。

    附錄A 圖A5 為實(shí)驗(yàn)數(shù)據(jù)采集地點(diǎn)的電網(wǎng)拓?fù)浣Y(jié)構(gòu)圖,數(shù)據(jù)采集裝置的安裝位置已在圖中標(biāo)記。圖A5(a)為教學(xué)樓低壓側(cè)的拓?fù)?,包括相互?dú)立的1 號(hào)與2 號(hào)電梯、中央空調(diào)、照明設(shè)備等,本文選擇1 號(hào)電梯負(fù)荷為待辨識(shí)負(fù)荷,采集的部分電梯負(fù)荷數(shù)據(jù)見附錄B。圖A5(b)為能源站低壓側(cè)的拓?fù)渑c其部分負(fù)荷,為保障能源站供電的安全可靠性,其中包含兩支獨(dú)立的進(jìn)線,電源進(jìn)線為消防負(fù)荷供電,非消防負(fù)荷進(jìn)線為其他電力負(fù)荷供電,本文選擇鍋爐負(fù)荷作為待辨識(shí)負(fù)荷。

    圖2 為待辨識(shí)負(fù)荷的有功功率和無(wú)功功率的聯(lián)合分布圖??梢钥闯?,電梯負(fù)荷和鍋爐負(fù)荷的有功和無(wú)功功率變化范圍都非常廣泛,并且分布較隨機(jī),不能簡(jiǎn)單地將其聚類處理劃分為有限個(gè)狀態(tài),且其有功功率并非維持在一個(gè)或多個(gè)穩(wěn)定值附近。電梯負(fù)荷和鍋爐負(fù)荷的狀態(tài)多變復(fù)雜,是典型的Type-Ⅲ型負(fù)荷,因此可以作為驗(yàn)證本文算法有效性的待辨識(shí)負(fù)荷。

    圖2 待辨識(shí)負(fù)荷的有功、無(wú)功功率聯(lián)合分布圖Fig.2 Diagrams of active and reactive power joint distribution for load to be identified

    3.3 時(shí)間窗口長(zhǎng)度和特征數(shù)目對(duì)負(fù)荷辨識(shí)的影響

    對(duì)采集的數(shù)據(jù)集進(jìn)行預(yù)處理和特征選擇時(shí),選擇不同的時(shí)間窗口長(zhǎng)度w和特征數(shù)量n會(huì)產(chǎn)生不同的辨識(shí)結(jié)果。表1 為固定w=50 s 時(shí),不同的特征數(shù)目對(duì)ResNet-50 網(wǎng)絡(luò)每輪訓(xùn)練時(shí)間的影響。由表1 可知,隨著特征數(shù)目的增加,輸入數(shù)據(jù)維度變大,網(wǎng)絡(luò)的訓(xùn)練時(shí)間也會(huì)隨之變長(zhǎng),因此采用互信息進(jìn)行特征選擇可以縮短模型的訓(xùn)練時(shí)間,從而提高模型運(yùn)算效率。

    表1 ResNet-50 網(wǎng)絡(luò)每輪的訓(xùn)練時(shí)間Table 1 Training time of each epoch in ResNet-50 network

    圖3 展示了算法辨識(shí)準(zhǔn)確率與時(shí)間長(zhǎng)度和特征數(shù)目的關(guān)系。由圖可見,隨著n的增加,負(fù)荷辨識(shí)模型可從輸入數(shù)據(jù)中提取更多的隱含特征,但當(dāng)n達(dá)到某一數(shù)值后,選擇的負(fù)荷特征中開始出現(xiàn)重疊信息或是無(wú)關(guān)信息,反而會(huì)對(duì)模型辨識(shí)產(chǎn)生干擾,降低辨識(shí)的準(zhǔn)確率;隨著w的增加,模型可從一個(gè)時(shí)間窗口內(nèi)學(xué)習(xí)待辨識(shí)負(fù)荷狀態(tài)之間的時(shí)間連續(xù)性,從而減少辨識(shí)結(jié)果非連續(xù)性的現(xiàn)象,但當(dāng)w覆蓋負(fù)荷的最大運(yùn)行周期值時(shí),模型逐漸達(dá)到飽和,負(fù)荷辨識(shí)準(zhǔn)確率便不會(huì)隨著w的增大而增大。綜合考慮算法的運(yùn)算量和準(zhǔn)確率,本文選擇時(shí)間窗口長(zhǎng)度w=60 s、特征數(shù)量n=40 進(jìn)行特征選擇和模型訓(xùn)練。

    圖3 時(shí)間長(zhǎng)度和特征數(shù)量與辨識(shí)準(zhǔn)確率的關(guān)系Fig.3 Relationship of time length, number of characteristics and identification accuracy

    3.4 ResNet-50 網(wǎng)絡(luò)訓(xùn)練過程

    從實(shí)驗(yàn)數(shù)據(jù)集中選擇周六的數(shù)據(jù)作為網(wǎng)絡(luò)的訓(xùn)練數(shù)據(jù),選擇周一的數(shù)據(jù)作為測(cè)試數(shù)據(jù),并且在訓(xùn)練數(shù)據(jù)中將80%用于訓(xùn)練,其余20%用于模型的評(píng)估,從而選擇出效果最佳的模型。在訓(xùn)練過程中,采用Adam 優(yōu)化器,設(shè)置初始學(xué)習(xí)率為0.001,并使用線性整流函數(shù)(ReLu)以加快模型的收斂速度,用交叉熵函數(shù)計(jì)算分類誤差。附錄A 圖A6 為損失函數(shù)值隨訓(xùn)練過程迭代次數(shù)的變化曲線圖。可以看出,在迭代次數(shù)達(dá)到30 左右時(shí),無(wú)論是訓(xùn)練誤差還是評(píng)估誤差,模型的損失值都趨于平緩,達(dá)到收斂狀態(tài)。

    3.5 與其他辨識(shí)算法的對(duì)比

    本文將幾種常用的非侵入式負(fù)荷辨識(shí)算法作為對(duì)比算法進(jìn)行測(cè)試,包括k最近鄰(KNN)[24]、隨機(jī)森林(RF)[25]、HMM[14]、長(zhǎng)短期記憶(LSTM)網(wǎng)絡(luò)[12]、CNN[11]以及CNN-LSTM 模型。其中,選擇HMM和CNN 算法是為了對(duì)比本文算法相較于單獨(dú)模型的優(yōu)化效果。KNN 和RF 算法是經(jīng)典的分類識(shí)別算法,而LSTM 網(wǎng)絡(luò)模型是處理時(shí)間序列的神經(jīng)網(wǎng)絡(luò),這3 種算法均被普遍應(yīng)用于非侵入式負(fù)荷辨識(shí)領(lǐng)域中。

    待辨識(shí)負(fù)荷的辨識(shí)結(jié)果如表2 所示,圖4 展示了幾種負(fù)荷辨識(shí)算法在相同時(shí)間內(nèi)的辨識(shí)結(jié)果。由圖4 可見,LSTM 和CNN 算法的辨識(shí)結(jié)果中都存在一定的狀態(tài)斷點(diǎn)問題,導(dǎo)致負(fù)荷辨識(shí)的準(zhǔn)確率降低;而2 層模型CNN-LSTM 和CNN-HMM 相較于單獨(dú)模型在一定程度上優(yōu)化了辨識(shí)結(jié)果的非連續(xù)性,提高了負(fù)荷辨識(shí)的準(zhǔn)確率,且HMM 比LSTM 的優(yōu)化效果更為顯著,本文提出的CNN-HMM 算法具有更高的精度。結(jié)合表2 中的評(píng)價(jià)指標(biāo)對(duì)比,可在數(shù)據(jù)上證實(shí)這一觀點(diǎn)。雖然在精確度P′和召回率R上,本文算法并非是最優(yōu)的,即本文算法在所有對(duì)比算法中并不具有最強(qiáng)的抗誤檢和抗漏檢的能力,但是對(duì)于能綜合反映辨識(shí)結(jié)果的可靠程度指標(biāo)Fscore,提出的算法取得了0.987 9 的最高值;對(duì)于反映負(fù)荷辨識(shí)結(jié)果與實(shí)際運(yùn)行狀態(tài)之間差異的Wasserstein 距離,提出的算法取得了最小值,而d值越小表明算法解決識(shí)別錯(cuò)誤斷點(diǎn)的能力越強(qiáng)。綜上所述,本文算法對(duì)于Type-Ⅲ型工商業(yè)負(fù)荷的辨識(shí)效果優(yōu)于其他對(duì)比算法,并且具有最高的辨識(shí)準(zhǔn)確率,能達(dá)到97.9%。

    表2 待辨識(shí)負(fù)荷的辨識(shí)結(jié)果Table 2 Identification results of load to be identified

    圖4 辨識(shí)結(jié)果對(duì)比Fig.4 Comparison of identification results

    3.6 用戶用電行為分析

    用戶的用電行為是指在什么時(shí)間用了什么電器產(chǎn)生了多大的負(fù)荷,通過NILM 技術(shù)可以方便地獲得各電器的使用情況,進(jìn)而可以對(duì)各電器的用電信息如啟停時(shí)間、持續(xù)時(shí)間及用電量等進(jìn)行統(tǒng)計(jì),以得到用戶準(zhǔn)確的用電行為,從而幫助電網(wǎng)公司制定供需互動(dòng)策略。

    附錄A 圖A7 為根據(jù)負(fù)荷辨識(shí)結(jié)果得到的電梯負(fù)荷和鍋爐負(fù)荷一天的用能統(tǒng)計(jì)。由圖A7(a)可知,電梯負(fù)荷的運(yùn)行具有明顯的時(shí)間特性,主要用電行為集中在09:00—23:00,而在00:00—06:00 時(shí)段內(nèi)幾乎沒有運(yùn)作,且工作日的運(yùn)行時(shí)長(zhǎng)較周末的略長(zhǎng),這與電梯負(fù)荷所處的學(xué)校環(huán)境有關(guān),電梯的運(yùn)行時(shí)段與學(xué)生上下課時(shí)段重合。由圖A7(b)可知,鍋爐負(fù)荷的運(yùn)行幾乎與時(shí)間無(wú)關(guān),因其處于供能站內(nèi),要滿足站內(nèi)所有的供熱需求,所以鍋爐負(fù)荷每天的工作時(shí)間都很長(zhǎng),用能量很大,且工作日與周末的用能并沒有很大區(qū)別。

    4 結(jié)語(yǔ)

    為有效辨識(shí)Type-Ⅲ型負(fù)荷,本文提出了基于深度CNN-HMM 的非侵入式負(fù)荷辨識(shí)算法。首先,根據(jù)互信息原理進(jìn)行特征選擇并去除冗余特征,采用參數(shù)搜索確定最優(yōu)時(shí)間窗口長(zhǎng)度和特征數(shù)目;由于ResNet-50 網(wǎng)絡(luò)可有效解決深度網(wǎng)絡(luò)訓(xùn)練過程中梯度消失的問題,將其作為深度CNN 模型的架構(gòu),提取負(fù)荷的多維特征,并完成初次辨識(shí);根據(jù)初次辨識(shí)結(jié)果,計(jì)算HMM 的參數(shù),并完成負(fù)荷的二次辨識(shí),從而優(yōu)化CNN 模型辨識(shí)中存在的辨識(shí)結(jié)果非連續(xù)性問題;最后,與其他負(fù)荷辨識(shí)算法的對(duì)比測(cè)試結(jié)果表明,本文算法能有效辨識(shí)Type-Ⅲ型工商業(yè)電力負(fù)荷,具有更好的準(zhǔn)確率和辨識(shí)精度。

    下一步的工作重點(diǎn)是研究Type-Ⅲ型負(fù)荷運(yùn)行時(shí)的典型狀態(tài)辨識(shí)問題。

    附錄見本刊網(wǎng)絡(luò)版(http://www.aeps-info.com/aeps/ch/index.aspx),掃英文摘要后二維碼可以閱讀網(wǎng)絡(luò)全文。

    猜你喜歡
    互信息特征選擇準(zhǔn)確率
    乳腺超聲檢查診斷乳腺腫瘤的特異度及準(zhǔn)確率分析
    健康之家(2021年19期)2021-05-23 11:17:39
    不同序列磁共振成像診斷脊柱損傷的臨床準(zhǔn)確率比較探討
    2015—2017 年寧夏各天氣預(yù)報(bào)參考產(chǎn)品質(zhì)量檢驗(yàn)分析
    高速公路車牌識(shí)別標(biāo)識(shí)站準(zhǔn)確率驗(yàn)證法
    Kmeans 應(yīng)用與特征選擇
    電子制作(2017年23期)2017-02-02 07:17:06
    基于互信息的貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)學(xué)習(xí)
    聯(lián)合互信息水下目標(biāo)特征選擇算法
    改進(jìn)的互信息最小化非線性盲源分離算法
    基于增量式互信息的圖像快速匹配方法
    基于特征選擇和RRVPMCD的滾動(dòng)軸承故障診斷方法
    久久久久精品性色| 超碰av人人做人人爽久久| 国产高清有码在线观看视频| 国产淫片久久久久久久久| 日本黄色片子视频| a级一级毛片免费在线观看| 蜜桃在线观看..| 欧美区成人在线视频| 男人舔奶头视频| 99久久精品热视频| 亚洲国产精品一区三区| 丝袜脚勾引网站| 蜜桃亚洲精品一区二区三区| 大又大粗又爽又黄少妇毛片口| 国产在线视频一区二区| 亚洲av国产av综合av卡| 五月玫瑰六月丁香| 国模一区二区三区四区视频| 特大巨黑吊av在线直播| 国产在线免费精品| 成人二区视频| 日本vs欧美在线观看视频 | 国产精品无大码| 日日摸夜夜添夜夜添av毛片| 大香蕉97超碰在线| 亚洲欧美日韩卡通动漫| 亚洲精品久久久久久婷婷小说| 精品一品国产午夜福利视频| 性色avwww在线观看| 亚洲av成人精品一二三区| 国产在线视频一区二区| 最近中文字幕高清免费大全6| 亚洲精品一二三| 亚洲国产av新网站| 欧美 日韩 精品 国产| 亚洲真实伦在线观看| 少妇 在线观看| 日韩av在线免费看完整版不卡| 欧美丝袜亚洲另类| 永久免费av网站大全| 日韩欧美精品免费久久| 久久国产精品大桥未久av | 日韩大片免费观看网站| 美女中出高潮动态图| 噜噜噜噜噜久久久久久91| 亚洲国产精品999| 国产精品一区二区在线不卡| 欧美变态另类bdsm刘玥| 日韩不卡一区二区三区视频在线| 欧美3d第一页| 亚洲精品,欧美精品| 在线免费观看不下载黄p国产| 亚洲国产高清在线一区二区三| 欧美xxxx黑人xx丫x性爽| 性高湖久久久久久久久免费观看| 我要看日韩黄色一级片| 国产av一区二区精品久久 | 国产成人freesex在线| 18禁在线播放成人免费| 国产av一区二区精品久久 | 身体一侧抽搐| 亚洲av免费高清在线观看| 丝袜脚勾引网站| 我的老师免费观看完整版| 女性被躁到高潮视频| 高清黄色对白视频在线免费看 | 国内精品宾馆在线| 夜夜看夜夜爽夜夜摸| 日本一二三区视频观看| 欧美日韩一区二区视频在线观看视频在线| av卡一久久| 国产白丝娇喘喷水9色精品| 日本爱情动作片www.在线观看| 亚洲在久久综合| 青春草亚洲视频在线观看| 91精品一卡2卡3卡4卡| 激情 狠狠 欧美| 蜜桃在线观看..| 综合色丁香网| 一级爰片在线观看| 久久精品久久久久久噜噜老黄| 亚洲不卡免费看| 国产精品av视频在线免费观看| 国产精品久久久久久精品古装| 亚洲熟女精品中文字幕| 欧美极品一区二区三区四区| 国产精品av视频在线免费观看| 欧美成人午夜免费资源| 2018国产大陆天天弄谢| 亚洲国产日韩一区二区| 久久久久网色| 精品国产乱码久久久久久小说| 国产淫语在线视频| 久久 成人 亚洲| 国产乱人视频| 97在线人人人人妻| 天堂8中文在线网| 国产 一区精品| 精品国产乱码久久久久久小说| 亚洲精品自拍成人| 老司机影院毛片| 亚洲在久久综合| 久久这里有精品视频免费| 国产精品嫩草影院av在线观看| 99视频精品全部免费 在线| 久久97久久精品| 亚洲精品久久久久久婷婷小说| 日韩视频在线欧美| 亚洲高清免费不卡视频| 欧美精品一区二区免费开放| av专区在线播放| 性色avwww在线观看| 3wmmmm亚洲av在线观看| 免费观看av网站的网址| 女性生殖器流出的白浆| 人人妻人人爽人人添夜夜欢视频 | 18禁在线无遮挡免费观看视频| 久久人人爽av亚洲精品天堂 | 卡戴珊不雅视频在线播放| 日本与韩国留学比较| 蜜桃在线观看..| 国产黄色视频一区二区在线观看| 最近最新中文字幕免费大全7| 人妻制服诱惑在线中文字幕| 麻豆精品久久久久久蜜桃| 2018国产大陆天天弄谢| av.在线天堂| 国产深夜福利视频在线观看| 青青草视频在线视频观看| 91午夜精品亚洲一区二区三区| 国产成人精品一,二区| 欧美成人一区二区免费高清观看| 午夜福利在线观看免费完整高清在| 亚洲不卡免费看| 国产高清有码在线观看视频| 亚洲精品,欧美精品| 少妇的逼水好多| 中文字幕精品免费在线观看视频 | 久久亚洲国产成人精品v| 国产美女午夜福利| 国产欧美另类精品又又久久亚洲欧美| 国产高清不卡午夜福利| 国产爽快片一区二区三区| 啦啦啦中文免费视频观看日本| av播播在线观看一区| 日韩成人伦理影院| 青春草国产在线视频| 国产精品一区二区在线不卡| 国产黄色免费在线视频| av视频免费观看在线观看| 亚洲精品日本国产第一区| 久久精品国产亚洲av涩爱| 国产亚洲精品久久久com| 色吧在线观看| 成人亚洲欧美一区二区av| 成年av动漫网址| 亚洲,一卡二卡三卡| 亚洲熟女精品中文字幕| 国产探花极品一区二区| 国产视频内射| 亚洲精品国产色婷婷电影| 国产黄色免费在线视频| 尾随美女入室| 国产精品国产三级国产专区5o| 成人影院久久| 汤姆久久久久久久影院中文字幕| 少妇的逼水好多| 建设人人有责人人尽责人人享有的 | 日韩一区二区三区影片| 久久久精品94久久精品| 午夜日本视频在线| 国产精品国产av在线观看| 少妇熟女欧美另类| 免费大片18禁| 最黄视频免费看| 三级经典国产精品| 新久久久久国产一级毛片| 日韩欧美精品免费久久| 久久久久视频综合| 人妻夜夜爽99麻豆av| 亚洲丝袜综合中文字幕| 黄色视频在线播放观看不卡| 国产亚洲最大av| av播播在线观看一区| 日韩,欧美,国产一区二区三区| 成年av动漫网址| 亚洲精品国产av成人精品| 777米奇影视久久| 男人舔奶头视频| 国产探花极品一区二区| av线在线观看网站| 热99国产精品久久久久久7| 九色成人免费人妻av| 国产极品天堂在线| 一个人看视频在线观看www免费| 国产爽快片一区二区三区| 婷婷色综合www| 国产男女内射视频| 自拍偷自拍亚洲精品老妇| 国产精品一及| 免费久久久久久久精品成人欧美视频 | 国产69精品久久久久777片| 91精品伊人久久大香线蕉| 亚洲av免费高清在线观看| 亚洲不卡免费看| 色5月婷婷丁香| 99热这里只有精品一区| 国产精品嫩草影院av在线观看| 在线观看国产h片| 男女下面进入的视频免费午夜| 亚洲电影在线观看av| 小蜜桃在线观看免费完整版高清| 国产无遮挡羞羞视频在线观看| 大片免费播放器 马上看| 一边亲一边摸免费视频| 久久99蜜桃精品久久| av免费在线看不卡| 欧美精品国产亚洲| 久久久久人妻精品一区果冻| 免费观看无遮挡的男女| 亚洲婷婷狠狠爱综合网| 亚洲欧洲日产国产| 国产精品秋霞免费鲁丝片| 纵有疾风起免费观看全集完整版| 超碰97精品在线观看| 国产成人a∨麻豆精品| 九色成人免费人妻av| 最黄视频免费看| av网站免费在线观看视频| 日本黄色日本黄色录像| 国产日韩欧美在线精品| 免费看av在线观看网站| 久久午夜福利片| 成人国产av品久久久| 人妻 亚洲 视频| 国产一级毛片在线| 日本欧美国产在线视频| 观看美女的网站| 亚洲自偷自拍三级| 免费观看a级毛片全部| 亚洲第一区二区三区不卡| 尾随美女入室| 又黄又爽又刺激的免费视频.| 在线免费观看不下载黄p国产| 欧美成人一区二区免费高清观看| 大片电影免费在线观看免费| 久久国产乱子免费精品| 制服丝袜香蕉在线| 国产成人精品福利久久| 人妻 亚洲 视频| 国产深夜福利视频在线观看| 亚洲欧美一区二区三区国产| 六月丁香七月| 成人漫画全彩无遮挡| 看免费成人av毛片| 夜夜爽夜夜爽视频| 国产精品蜜桃在线观看| 日韩 亚洲 欧美在线| 丰满少妇做爰视频| 综合色丁香网| 日韩av免费高清视频| 日韩伦理黄色片| 亚洲精品aⅴ在线观看| 亚洲成色77777| 伦精品一区二区三区| 丝袜喷水一区| 亚洲精品456在线播放app| 男人爽女人下面视频在线观看| 日日摸夜夜添夜夜爱| 久久久精品94久久精品| 日本午夜av视频| 午夜福利网站1000一区二区三区| 日韩,欧美,国产一区二区三区| tube8黄色片| 啦啦啦啦在线视频资源| 一本一本综合久久| 午夜福利视频精品| 综合色丁香网| 欧美另类一区| 又粗又硬又长又爽又黄的视频| 日韩不卡一区二区三区视频在线| 2022亚洲国产成人精品| 亚洲欧美一区二区三区黑人 | 色综合色国产| 久久亚洲国产成人精品v| 一级毛片aaaaaa免费看小| 国产 一区 欧美 日韩| 18禁动态无遮挡网站| 水蜜桃什么品种好| 午夜福利影视在线免费观看| 啦啦啦中文免费视频观看日本| 免费av中文字幕在线| 中文精品一卡2卡3卡4更新| 婷婷色综合大香蕉| 久久热精品热| 国产精品一区二区性色av| 熟女电影av网| 成人免费观看视频高清| 国产成人freesex在线| 人妻一区二区av| 免费看不卡的av| 日韩一区二区视频免费看| 欧美97在线视频| 最近最新中文字幕大全电影3| 精品99又大又爽又粗少妇毛片| 免费看日本二区| 超碰97精品在线观看| 建设人人有责人人尽责人人享有的 | 国产黄色免费在线视频| 国产精品国产三级国产av玫瑰| 亚洲,一卡二卡三卡| 亚洲aⅴ乱码一区二区在线播放| 免费大片18禁| 少妇裸体淫交视频免费看高清| 一级二级三级毛片免费看| 美女xxoo啪啪120秒动态图| 少妇的逼好多水| 久久久国产一区二区| 欧美+日韩+精品| 日韩成人伦理影院| 亚洲精华国产精华液的使用体验| 久久久精品免费免费高清| 最近中文字幕2019免费版| 亚洲真实伦在线观看| 高清日韩中文字幕在线| 91久久精品国产一区二区成人| 免费在线观看成人毛片| 久久精品熟女亚洲av麻豆精品| 黑丝袜美女国产一区| 免费观看a级毛片全部| 欧美日韩精品成人综合77777| 搡女人真爽免费视频火全软件| 色婷婷av一区二区三区视频| 亚洲自偷自拍三级| 一级a做视频免费观看| 国产视频首页在线观看| 天堂8中文在线网| 日日摸夜夜添夜夜爱| 成人无遮挡网站| 美女cb高潮喷水在线观看| 国产老妇伦熟女老妇高清| 色婷婷久久久亚洲欧美| 夜夜骑夜夜射夜夜干| 国产精品欧美亚洲77777| 97超视频在线观看视频| 欧美日韩在线观看h| 最新中文字幕久久久久| 精品久久久久久电影网| 97超视频在线观看视频| 亚洲内射少妇av| 99九九线精品视频在线观看视频| 日本av手机在线免费观看| 美女福利国产在线 | 久久精品夜色国产| 亚洲精品中文字幕在线视频 | 丰满人妻一区二区三区视频av| 亚洲自偷自拍三级| 精品一区在线观看国产| 天美传媒精品一区二区| 国产成人a∨麻豆精品| 中文字幕精品免费在线观看视频 | 中文精品一卡2卡3卡4更新| av一本久久久久| 日本av手机在线免费观看| 一区二区三区免费毛片| 精品一区二区三区视频在线| 麻豆精品久久久久久蜜桃| 国产欧美日韩一区二区三区在线 | 天堂中文最新版在线下载| 91午夜精品亚洲一区二区三区| 伊人久久精品亚洲午夜| 日韩,欧美,国产一区二区三区| 黑丝袜美女国产一区| 22中文网久久字幕| 亚洲熟女精品中文字幕| 成人国产av品久久久| 夜夜看夜夜爽夜夜摸| 亚洲精品亚洲一区二区| 少妇丰满av| 少妇的逼好多水| 国产黄片视频在线免费观看| 久久 成人 亚洲| 国产精品嫩草影院av在线观看| a 毛片基地| 国产视频内射| 尤物成人国产欧美一区二区三区| 国产亚洲av片在线观看秒播厂| 亚洲图色成人| 精品久久久噜噜| 高清视频免费观看一区二区| 精品国产露脸久久av麻豆| 国产精品国产三级专区第一集| 国产精品av视频在线免费观看| 九色成人免费人妻av| 91在线精品国自产拍蜜月| 欧美日韩视频精品一区| 99久久中文字幕三级久久日本| 噜噜噜噜噜久久久久久91| 亚洲精品自拍成人| 日本vs欧美在线观看视频 | 大陆偷拍与自拍| 久久精品人妻少妇| 久久99热这里只有精品18| 亚洲精品国产色婷婷电影| 王馨瑶露胸无遮挡在线观看| 51国产日韩欧美| 99久久人妻综合| 国产精品国产av在线观看| 欧美日韩综合久久久久久| 国产乱人视频| 久久99热这里只有精品18| 啦啦啦在线观看免费高清www| 欧美性感艳星| 日本一二三区视频观看| 久久精品国产亚洲网站| 日日啪夜夜爽| 国产精品99久久久久久久久| 精品国产一区二区三区久久久樱花 | 亚洲美女视频黄频| 国产精品精品国产色婷婷| 99久国产av精品国产电影| 欧美成人精品欧美一级黄| 在线观看一区二区三区激情| 免费看光身美女| 成人毛片60女人毛片免费| 视频区图区小说| 一本一本综合久久| 91狼人影院| 3wmmmm亚洲av在线观看| 国模一区二区三区四区视频| 欧美97在线视频| 免费人妻精品一区二区三区视频| 色视频在线一区二区三区| 久久精品国产亚洲av天美| 欧美老熟妇乱子伦牲交| 在线看a的网站| 精品一区二区免费观看| 亚洲欧美日韩卡通动漫| 国产人妻一区二区三区在| av线在线观看网站| 日日摸夜夜添夜夜添av毛片| av网站免费在线观看视频| 久久久久久伊人网av| 青春草亚洲视频在线观看| 男人爽女人下面视频在线观看| 国产在线一区二区三区精| 久久精品国产亚洲av涩爱| 亚洲精品久久久久久婷婷小说| 国产精品人妻久久久影院| 久久韩国三级中文字幕| 晚上一个人看的免费电影| 精品国产一区二区三区久久久樱花 | 亚洲精品久久久久久婷婷小说| 国产精品女同一区二区软件| 一区二区三区乱码不卡18| 中文字幕制服av| 日韩国内少妇激情av| 国产午夜精品久久久久久一区二区三区| 免费久久久久久久精品成人欧美视频 | 亚洲欧美一区二区三区国产| 国产成人精品一,二区| 一级毛片黄色毛片免费观看视频| 九九久久精品国产亚洲av麻豆| 久久97久久精品| 国产v大片淫在线免费观看| 亚洲精品一二三| 午夜福利在线在线| 噜噜噜噜噜久久久久久91| 简卡轻食公司| 国产欧美日韩精品一区二区| 1000部很黄的大片| 啦啦啦在线观看免费高清www| 国产精品久久久久成人av| 又黄又爽又刺激的免费视频.| 日本免费在线观看一区| 亚洲一级一片aⅴ在线观看| 久久午夜福利片| videos熟女内射| 插逼视频在线观看| 国产成人91sexporn| 日韩欧美精品免费久久| videos熟女内射| 免费观看av网站的网址| 亚洲自偷自拍三级| 男人爽女人下面视频在线观看| 国产淫语在线视频| 午夜福利高清视频| 91久久精品国产一区二区三区| 欧美xxxx性猛交bbbb| 日韩一区二区视频免费看| 久久久久久九九精品二区国产| 久久久国产一区二区| 亚洲人与动物交配视频| 熟女av电影| 一区二区三区精品91| 日韩成人av中文字幕在线观看| 色视频www国产| 国产精品99久久久久久久久| 我要看黄色一级片免费的| 亚洲三级黄色毛片| 你懂的网址亚洲精品在线观看| 爱豆传媒免费全集在线观看| 好男人视频免费观看在线| 国产一级毛片在线| 十八禁网站网址无遮挡 | 黑丝袜美女国产一区| 国产一区亚洲一区在线观看| 日本欧美国产在线视频| 成人一区二区视频在线观看| 国产黄色免费在线视频| 综合色丁香网| 99久久人妻综合| 国产成人精品久久久久久| 黑人猛操日本美女一级片| 超碰97精品在线观看| 亚洲,欧美,日韩| 日韩欧美精品免费久久| 久久国产精品男人的天堂亚洲 | 国产成人freesex在线| 亚洲怡红院男人天堂| 我要看黄色一级片免费的| 久久这里有精品视频免费| 国产精品成人在线| 天天躁日日操中文字幕| 又黄又爽又刺激的免费视频.| 男女边摸边吃奶| 黄色怎么调成土黄色| 丝瓜视频免费看黄片| 亚洲av日韩在线播放| 国产视频首页在线观看| 纯流量卡能插随身wifi吗| 欧美三级亚洲精品| 一级毛片电影观看| 亚洲av不卡在线观看| 成年美女黄网站色视频大全免费 | 最近2019中文字幕mv第一页| 不卡视频在线观看欧美| 我的女老师完整版在线观看| 最近的中文字幕免费完整| 国产成人精品福利久久| 亚洲精品国产av成人精品| 久久久欧美国产精品| 亚洲欧美成人综合另类久久久| 中文乱码字字幕精品一区二区三区| 少妇的逼水好多| 少妇人妻 视频| 日韩 亚洲 欧美在线| 久久久久久久久久久免费av| 少妇的逼好多水| 国产欧美日韩一区二区三区在线 | 黑人高潮一二区| 青春草国产在线视频| 黄色日韩在线| 国产精品爽爽va在线观看网站| 丝瓜视频免费看黄片| 22中文网久久字幕| av卡一久久| 人人妻人人看人人澡| 久久99热这里只有精品18| 少妇人妻精品综合一区二区| 国产精品一区二区在线观看99| 黑人猛操日本美女一级片| 国产久久久一区二区三区| 亚洲精品日韩在线中文字幕| 成人无遮挡网站| 亚洲久久久国产精品| 日本av手机在线免费观看| 99国产精品免费福利视频| 国产精品精品国产色婷婷| 国产又色又爽无遮挡免| 九色成人免费人妻av| 欧美一级a爱片免费观看看| 日韩大片免费观看网站| av在线app专区| 久久久久久伊人网av| 中文在线观看免费www的网站| 3wmmmm亚洲av在线观看| 午夜激情福利司机影院| 欧美三级亚洲精品| 久久久久性生活片| 永久网站在线| 亚洲aⅴ乱码一区二区在线播放| 国产午夜精品一二区理论片| 丰满人妻一区二区三区视频av| 蜜桃亚洲精品一区二区三区| 在线看a的网站| 国模一区二区三区四区视频| 亚洲在久久综合| 91久久精品国产一区二区三区| 又黄又爽又刺激的免费视频.| 看非洲黑人一级黄片| 最近手机中文字幕大全| 久久久久久久亚洲中文字幕| 精品一区二区免费观看| 久久久久久人妻| 亚洲精品视频女| 中文乱码字字幕精品一区二区三区| 天堂中文最新版在线下载| 国产精品久久久久成人av| 寂寞人妻少妇视频99o| 亚洲精品日本国产第一区| 永久免费av网站大全| 国产 精品1| 亚洲精品视频女| 国产视频内射| 哪个播放器可以免费观看大片| 高清毛片免费看| 亚洲欧美成人综合另类久久久| 免费看日本二区| 街头女战士在线观看网站| 一级黄片播放器| av网站免费在线观看视频| 在线 av 中文字幕| 肉色欧美久久久久久久蜜桃| 欧美日韩国产mv在线观看视频 |