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

    基于聯(lián)邦學(xué)習(xí)的短期負(fù)荷預(yù)測模型協(xié)同訓(xùn)練方法

    2022-11-14 03:19:45車亮徐茂盛崔秋實(shí)
    關(guān)鍵詞:負(fù)荷預(yù)測機(jī)器學(xué)習(xí)

    車亮 徐茂盛 崔秋實(shí)

    摘要:針對機(jī)器學(xué)習(xí)方法在電力系統(tǒng)短期負(fù)荷預(yù)測領(lǐng)域的應(yīng)用過程中,存在數(shù)據(jù)樣本不足、模型泛化能力差以及數(shù)據(jù)隱私保護(hù)要求較高等問題,以氣象、日期以及歷史負(fù)荷數(shù)據(jù)為輸入特征,構(gòu)建基于長短期記憶(Long Short-Term Memory,LSTM)網(wǎng)絡(luò)的短期負(fù)荷預(yù)測模型,提出基于聯(lián)邦學(xué)習(xí)(Federated Learning,F(xiàn)L)的短期負(fù)荷預(yù)測模型協(xié)同訓(xùn)練方法.通過分散訓(xùn)練、中心聚合的方式對模型參數(shù)進(jìn)行迭代更新,實(shí)現(xiàn)各負(fù)荷運(yùn)營商在保證數(shù)據(jù)隱私的情況下協(xié)同構(gòu)建預(yù)測模型.在GEFCom2012比賽的多個(gè)地區(qū)負(fù)荷數(shù)據(jù)集上進(jìn)行仿真驗(yàn)證,結(jié)果表明,所提方法在保證各運(yùn)營商數(shù)據(jù)隱私的同時(shí),有效提升了短期負(fù)荷預(yù)測準(zhǔn)確率,所訓(xùn)練出的模型在多場景下具有優(yōu)秀的泛化能力.

    關(guān)鍵詞:數(shù)據(jù)隱私;機(jī)器學(xué)習(xí);負(fù)荷預(yù)測;聯(lián)邦學(xué)習(xí)

    中圖分類號:TM715文獻(xiàn)標(biāo)志碼:A

    Federated Learning Based Coordinated Training Method of a Short-term Load Forecasting Model

    CHE Liang1,XU Maosheng1,CUI Qiushi2

    (1. College of Electrical and Information Engineering,Hunan University,Changsha 410082,China;2. School of Electrical Engineering,Chongqing University,Chongqing 400044,China)

    Abstract:Machine learning methods have been widely used in the field of short-term load forecasting of power systems. However,it is difficult for load operators to obtain high-performance forecasting models due to insufficient data samples,poor model generalization ability,and high data privacy protection requirements in the application process. In this paper,meteorological,date,and historical load are used as input features to construct a short-term load forecasting model based on Long Short-Term Memory(LSTM). A federated learning(FL)based coordinated training method of a short-term load forecasting model is proposed. The proposed method mainly iteratively updates model parameters through decentralized training and aggregation of centers,so as to realize cooperative construction of the prediction model by all load operators under the condition of data privacy. The simulation results based on the GEFCom2012 dataset show that the proposed method not only ensures the data privacy of operators but also effectively improves the forecasting accuracy of the load forecasting model,and the trained model has satisfied generalization ability in multiple scenarios.

    Key words:data privacy;machine learning;load forecasting;federated learning

    電力系統(tǒng)短期負(fù)荷預(yù)測是指根據(jù)過去一段時(shí)間的負(fù)荷以及日期、氣象等數(shù)據(jù),對未來短時(shí)間的負(fù)荷變化進(jìn)行預(yù)測[1].精準(zhǔn)的負(fù)荷預(yù)測是電網(wǎng)實(shí)施精確調(diào)度、制訂發(fā)電計(jì)劃以及設(shè)置合理備用的基礎(chǔ),是保證電力系統(tǒng)可靠運(yùn)行與增加經(jīng)濟(jì)效益的必要條件[2].隨著電網(wǎng)數(shù)字化程度不斷加深,電力公司陸續(xù)安裝了多種智能儀表,獲得了更多精細(xì)的負(fù)荷與氣象等數(shù)據(jù),這使得采用人工智能方法進(jìn)行精準(zhǔn)負(fù)荷預(yù)測變得可行[3].

    近年來,機(jī)器學(xué)習(xí)算法被廣泛應(yīng)用在短期負(fù)荷預(yù)測領(lǐng)域,研究者們對各種算法進(jìn)行改進(jìn),提出了許多解決短期負(fù)荷預(yù)測領(lǐng)域相關(guān)問題的方法[4].研究者將支持向量機(jī)、BP(Back Propagation)神經(jīng)網(wǎng)絡(luò)(Neural Network,NN)、LSTM網(wǎng)絡(luò)等多種機(jī)器學(xué)習(xí)算法應(yīng)用于負(fù)荷預(yù)測,并對其進(jìn)行改造以提升預(yù)測精度.文獻(xiàn)[5]提出基于深度信念網(wǎng)絡(luò)的短期負(fù)荷預(yù)測方法,將無監(jiān)督與有監(jiān)督訓(xùn)練結(jié)合,進(jìn)而提升預(yù)測準(zhǔn)確度;文獻(xiàn)[6]將Attention機(jī)制加入LSTM網(wǎng)絡(luò)構(gòu)建預(yù)測模型,從輸入數(shù)據(jù)中提取對負(fù)荷預(yù)測起到關(guān)鍵作用的特征;文獻(xiàn)[7]基于Hadoop架構(gòu)中的MapReduce框架構(gòu)建BP神經(jīng)網(wǎng)絡(luò)預(yù)測模型;文獻(xiàn)[8]構(gòu)建基于雙向長短期網(wǎng)絡(luò)的預(yù)測模型,并融入特征選擇方法,選擇出最佳輸入特征.一般來說,單一類型的神經(jīng)網(wǎng)絡(luò)只擅長挖掘數(shù)據(jù)的某一類特征,而負(fù)荷預(yù)測數(shù)據(jù)包括負(fù)荷、氣象和日期等多類特征,利用不同的神經(jīng)網(wǎng)絡(luò)能更有效提取不同數(shù)據(jù)類型的特征[9];文獻(xiàn)[10]先運(yùn)用卷積神經(jīng)網(wǎng)絡(luò)(Convolutional Neural Network,CNN)提取特征與負(fù)荷的聯(lián)系,最后將結(jié)果輸入到GRU網(wǎng)絡(luò)中進(jìn)行負(fù)荷預(yù)測,提升預(yù)測準(zhǔn)確率;文獻(xiàn)[11]將卷積神經(jīng)網(wǎng)絡(luò)(CNN)和長短期記憶(LSTM)網(wǎng)絡(luò)相結(jié)合以構(gòu)建短期負(fù)荷預(yù)測模型;文獻(xiàn)[12]綜合集合經(jīng)驗(yàn)?zāi)B(tài)分解(Ensemble Empirical Mode Decomposition,EEMD)、門控循環(huán)單元神經(jīng)網(wǎng)絡(luò)(Gated Recurrent Unit,GRU)和多元線性回歸(Multiple Linear Regression,MLR)方法構(gòu)建預(yù)測模型;文獻(xiàn)[13]提出了一種基于GRU-NN模型的短期電力負(fù)荷預(yù)測方法;文獻(xiàn)[14]提出了基于LSTM和XGBoost的組合預(yù)測模型.

    由上可見,目前機(jī)器學(xué)習(xí)方法在短期負(fù)荷預(yù)測領(lǐng)域的應(yīng)用已取得了一定的成果,但在實(shí)際應(yīng)用中,現(xiàn)有方法還存在以下問題:

    1)未解決實(shí)際情況下負(fù)荷數(shù)據(jù)樣本數(shù)量與特征不足的問題.現(xiàn)有基于機(jī)器學(xué)習(xí)的負(fù)荷預(yù)測方法依賴于預(yù)測者擁有樣本數(shù)足夠多、特征種類足夠豐富的訓(xùn)練數(shù)據(jù)集,但在實(shí)際負(fù)荷預(yù)測工作中,部分預(yù)測者往往難以滿足該條件(例如小規(guī)模電力公司、新接入電網(wǎng)的用戶和負(fù)荷聚合商等通常難以建立滿足負(fù)荷預(yù)測需求的數(shù)據(jù)量和特征足夠的訓(xùn)練數(shù)據(jù)集[15-16]).

    2)預(yù)測模型泛化能力問題.僅使用單個(gè)地區(qū)數(shù)據(jù)訓(xùn)練出的模型,難以直接移植到新的地區(qū)進(jìn)行預(yù)測,預(yù)測模型的泛化能力亟待增強(qiáng).

    3)數(shù)據(jù)隱私問題.若能實(shí)現(xiàn)不同地區(qū)間共享數(shù)據(jù)訓(xùn)練模型,就可以有效提升模型預(yù)測能力[17].但目前各電力公司均具有較高的數(shù)據(jù)保密要求,歷史負(fù)荷等數(shù)據(jù)是嚴(yán)格保密的,這限制了數(shù)據(jù)共享訓(xùn)練,因此相關(guān)數(shù)據(jù)隱私保密問題亟待解決.

    因此,目前亟需一個(gè)在保護(hù)負(fù)荷等數(shù)據(jù)隱私的前提下實(shí)現(xiàn)協(xié)同構(gòu)建負(fù)荷預(yù)測模型的方法.針對這一需求,本文引入聯(lián)邦學(xué)習(xí)(Federated Learning,F(xiàn)L)方法.聯(lián)邦學(xué)習(xí)可以實(shí)現(xiàn)多個(gè)數(shù)據(jù)持有者在不共享其數(shù)據(jù)情況下協(xié)同訓(xùn)練機(jī)器學(xué)習(xí)模型,在保護(hù)數(shù)據(jù)隱私的同時(shí)降低數(shù)據(jù)通信量,提高通信效率[18].目前,聯(lián)邦學(xué)習(xí)已被研究者們應(yīng)用在一些需要考慮數(shù)據(jù)隱私及降低通信開銷的場景[19].文獻(xiàn)[20]將深度強(qiáng)化學(xué)習(xí)模型進(jìn)行聯(lián)邦,實(shí)現(xiàn)高準(zhǔn)確率、低時(shí)延的設(shè)備節(jié)點(diǎn)選擇;文獻(xiàn)[21]將聯(lián)邦學(xué)習(xí)應(yīng)用于物聯(lián)網(wǎng)框架下的邊緣計(jì)算過程;文獻(xiàn)[22]將聯(lián)邦學(xué)習(xí)應(yīng)用于霧計(jì)算過程;文獻(xiàn)[23]將聯(lián)邦學(xué)習(xí)應(yīng)用在電力計(jì)量系統(tǒng)中,分別應(yīng)用于計(jì)量裝置故障診斷、無線接入安全檢測以及電費(fèi)風(fēng)險(xiǎn)識別場景,均取得了較好的效果,證明了將聯(lián)邦學(xué)習(xí)方法應(yīng)用在電力系統(tǒng)領(lǐng)域具有實(shí)用性和可行性.

    針對上述負(fù)荷預(yù)測模型構(gòu)建所需數(shù)據(jù)樣本數(shù)量與特征不足、模型泛化能力差以及負(fù)荷相關(guān)數(shù)據(jù)高度保密無法共享訓(xùn)練的問題,本文引入聯(lián)邦學(xué)習(xí)方法,以LSTM網(wǎng)絡(luò)的短期負(fù)荷預(yù)測模型為基礎(chǔ),提出基于聯(lián)邦學(xué)習(xí)的短期負(fù)荷預(yù)測模型協(xié)同訓(xùn)練方法,結(jié)合短期負(fù)荷預(yù)測場景的特點(diǎn),對數(shù)據(jù)傳輸過程與模型聚合過程進(jìn)行優(yōu)化改進(jìn).該方法實(shí)現(xiàn)了各負(fù)荷運(yùn)營商在保護(hù)數(shù)據(jù)隱私的前提下對模型進(jìn)行協(xié)同訓(xùn)練,并且有效提升了短期負(fù)荷預(yù)測模型的預(yù)測準(zhǔn)確率和泛化能力.

    1基于LSTM的短期負(fù)荷預(yù)測模型

    1.1長短期記憶(LSTM)網(wǎng)絡(luò)

    LSTM網(wǎng)絡(luò)是由循環(huán)神經(jīng)網(wǎng)絡(luò)改進(jìn)而來的,解決了循環(huán)神經(jīng)網(wǎng)絡(luò)存在的長期依賴問題.相比常規(guī)循環(huán)神經(jīng)網(wǎng)絡(luò),LSTM添加了輸入門、遺忘門和輸出門[24].LSTM網(wǎng)絡(luò)基本單元如圖1所示.

    ft=σ(Wfht-1+Wfxt+bf)(1)

    式中:Wf、Wi、Wo分別為遺忘門、輸入門、輸出門的權(quán)重矩陣;bf、bi、bo分別為遺忘門、輸入門、輸出門的輸出偏置項(xiàng);*表示矩陣的Hadamard積.

    1.2基于LSTM網(wǎng)絡(luò)的短期負(fù)荷預(yù)測模型構(gòu)建

    由于負(fù)荷數(shù)據(jù)數(shù)量級較高,與溫度、日期等特征數(shù)量級差距較大,不利于模型訓(xùn)練.為了消除各類特征之間量綱不同的影響,使模型快速收斂,并且考慮到各負(fù)荷運(yùn)營商之間負(fù)荷數(shù)據(jù)數(shù)量級存在較大差異,故利用本文方法對輸入數(shù)據(jù)進(jìn)行歸一化處理,將所有輸入數(shù)據(jù)歸一化在[0,1]內(nèi),如下式所示:

    式中:x*t為歸一化后的數(shù)據(jù);xt為歸一化前的數(shù)據(jù);xmax和xmin分別表示該特征數(shù)據(jù)的最大值和最小值.LSTM網(wǎng)絡(luò)模型結(jié)構(gòu)如圖2所示.

    負(fù)荷預(yù)測效果評價(jià)指標(biāo)采用均方根誤差(Root Mean Square Error,RMSE)和平均絕對百分比誤差(Mean Absolute Percentage Error,MAPE),指標(biāo)值越小,表示預(yù)測值與真實(shí)值之間的誤差越小,模型的預(yù)測效果越好.指標(biāo)計(jì)算公式分別如式(8)和式(9)所示.

    2基于FL的短期負(fù)荷預(yù)測模型協(xié)同訓(xùn)練方法

    2.1方法概述

    基于聯(lián)邦學(xué)習(xí)的短期負(fù)荷預(yù)測模型協(xié)同訓(xùn)練方法能在保護(hù)各負(fù)荷運(yùn)營商數(shù)據(jù)隱私的同時(shí),利用各個(gè)運(yùn)營商所擁有的負(fù)荷數(shù)據(jù)資源協(xié)同訓(xùn)練短期負(fù)荷預(yù)測模型,有效提高模型的預(yù)測準(zhǔn)確率及多場景泛化能力.

    本文方法整體系統(tǒng)架構(gòu)如圖3所示.其中,本文第1節(jié)所構(gòu)建的負(fù)荷預(yù)測模型部署在服務(wù)器端與各個(gè)負(fù)荷運(yùn)營商端,負(fù)荷運(yùn)營商可以是某地區(qū)電力公司或某些電力用戶的集合.

    在結(jié)構(gòu)上,本文方法整體組成包括服務(wù)器端、通信網(wǎng)絡(luò)、本地負(fù)荷運(yùn)營商及其對應(yīng)的負(fù)荷、氣象和日期數(shù)據(jù).令L={l1,l2,…,ln}表示所有參與協(xié)同訓(xùn)練的負(fù)荷運(yùn)營商的集合,每一個(gè)lk(k=1,2,…,n)具有對應(yīng)的訓(xùn)練數(shù)據(jù)以及進(jìn)行本地模型訓(xùn)練、模型推理的能力.圖3中,令V表示訓(xùn)練數(shù)據(jù)(包括負(fù)荷數(shù)據(jù)、氣象數(shù)據(jù)、日期數(shù)據(jù))集合,則第k個(gè)負(fù)荷運(yùn)營商lk所擁有的訓(xùn)練數(shù)據(jù)表示為Vl,k;W表示服務(wù)器下發(fā)到本地的模型權(quán)重參數(shù);令Wk(k=1,2,…,n)表示第k個(gè)本地端lk從本地服務(wù)端上傳至服務(wù)器端的模型權(quán)重參數(shù).

    整體流程主要包括預(yù)測模型下發(fā)、負(fù)荷運(yùn)營商本地訓(xùn)練、各運(yùn)營商模型上傳以及模型聚合4個(gè)步驟.具體步驟如算法1所示.

    1)以LSTM網(wǎng)絡(luò)的短期負(fù)荷預(yù)測模型為基礎(chǔ),服務(wù)器端初始化各參數(shù),并下發(fā)模型到本地端.

    2)本地端接收下發(fā)的模型,用該模型參數(shù)更新自身模型參數(shù),并以自身數(shù)據(jù)訓(xùn)練模型.

    3)將已訓(xùn)練完成的模型上傳至服務(wù)器端,采用不傳輸梯度而是直接傳輸完整權(quán)重參數(shù)的方法降低數(shù)據(jù)泄露風(fēng)險(xiǎn),避免被針對模型參數(shù)進(jìn)行攻擊的算法竊取數(shù)據(jù).

    4)服務(wù)器端執(zhí)行模型聚合,考慮電力負(fù)荷預(yù)測場景的特點(diǎn),在模型聚合階段,根據(jù)各參與方數(shù)據(jù)集質(zhì)量形成適用于負(fù)荷預(yù)測場景的模型聚合策略.最終,形成完整的基于聯(lián)邦學(xué)習(xí)的短期負(fù)荷預(yù)測模型協(xié)同訓(xùn)練方法.

    選取算法1步驟中較為主要的10)~12)在2.2節(jié)詳細(xì)介紹;步驟13)和14)在2.3節(jié)詳細(xì)介紹;步驟16)和17)在2.4節(jié)詳細(xì)介紹.

    2.2模型下發(fā)與本地訓(xùn)練更新

    在服務(wù)器端進(jìn)行模型預(yù)訓(xùn)練,得到服務(wù)器端模型權(quán)重參數(shù)W.通過通信網(wǎng)絡(luò)將W下發(fā)至各個(gè)參與

    訓(xùn)練的聯(lián)邦參與本地端lk.各lk接收服務(wù)器下發(fā)的W,更新本地模型參數(shù)Wk.隨后以其所擁有的數(shù)據(jù)集Vl,k作為本地訓(xùn)練樣本集,進(jìn)行本地訓(xùn)練并更新Wk.

    模型訓(xùn)練過程采用Adam(Adaptive Moment Estimation)優(yōu)化算法.該算法在RMSprop算法的基礎(chǔ)上添加了動量項(xiàng),根據(jù)t時(shí)刻梯度gt的一階矩估計(jì)和二階矩估計(jì)來動態(tài)調(diào)整每個(gè)參數(shù)的學(xué)習(xí)率.該算法的

    主要優(yōu)點(diǎn)為:經(jīng)過偏置校正后,每一次迭代學(xué)習(xí)率都存在于一個(gè)確定的范圍,使得參數(shù)較為穩(wěn)定.

    基于Adam的模型訓(xùn)練過程描述如下:

    1)根據(jù)t時(shí)刻的梯度gt分別計(jì)算梯度和梯度平方的指數(shù)移動平均數(shù)mt和vt,β1與β2分別取0.9、0.999,m0與v0均設(shè)為0:

    4)模型權(quán)重參數(shù)更新如式(13)所示,式中Wk表示參數(shù)更新后的第k個(gè)本地模型權(quán)重,ΣΔWk,t表示第k個(gè)本地模型總的權(quán)重變化量:

    Wk=W+ΣΔWk,t(13)

    采用平均絕對誤差(Mean Absolute Error,MAE)作為模型的目標(biāo)損失函數(shù)yMAE.

    2.3模型上傳

    本地模型訓(xùn)練完成后,每個(gè)負(fù)荷運(yùn)營商lk經(jīng)過通信網(wǎng)絡(luò)上傳模型權(quán)重參數(shù)Wk至服務(wù)器端,但在參數(shù)傳輸過程中存在多種不同的傳輸方式.

    經(jīng)典的聯(lián)邦學(xué)習(xí)方法在訓(xùn)練迭代時(shí)傳輸梯度數(shù)據(jù)gt,每位本地參與者首先將當(dāng)前訓(xùn)練的模型梯度相關(guān)系數(shù)上傳至服務(wù)器,服務(wù)器通過求取每位參與者最新的梯度相關(guān)系數(shù)的平均值來更新全局模型;在下一次迭代中,參與者再下載最新的全局模型參數(shù)繼續(xù)訓(xùn)練-但是,此種方法已被指出存在隱私泄露風(fēng)險(xiǎn)-有學(xué)者發(fā)現(xiàn)在聯(lián)邦過程中,如果發(fā)生梯度參數(shù)泄露,則可用算法通過梯度反推出該次迭代對應(yīng)的訓(xùn)練數(shù)據(jù)-例如較為典型的Deep Leakage from Gradients(DLG)算法[25]:首先假設(shè)虛擬的輸入和對應(yīng)的標(biāo)簽;隨后,計(jì)算虛擬的梯度,在此基礎(chǔ)上優(yōu)化虛擬梯度與真實(shí)梯度之間的距離,通過梯度的匹配使得假設(shè)的虛擬數(shù)據(jù)不斷靠近原始數(shù)據(jù);最終,獲取參與方的真實(shí)訓(xùn)練數(shù)據(jù).

    為了避免由于梯度參數(shù)泄露導(dǎo)致反推出對應(yīng)數(shù)據(jù)樣本的風(fēng)險(xiǎn),本文采用不傳輸梯度而是直接傳輸完整權(quán)重參數(shù)的方法,在參與方本地模型訓(xùn)練完成后,上傳整體模型權(quán)重參數(shù)Wk,而非單個(gè)數(shù)據(jù)樣本訓(xùn)練時(shí)產(chǎn)生的梯度-權(quán)重參數(shù)經(jīng)過該參與方使用自身全部數(shù)據(jù)通過Adam算法進(jìn)行多次迭代訓(xùn)練,根據(jù)上述Adam訓(xùn)練過程及DLG算法原理,此時(shí)已無法推導(dǎo)出每一條數(shù)據(jù)樣本訓(xùn)練時(shí)所對應(yīng)的梯度計(jì)算值,從而無法推導(dǎo)出對應(yīng)的原始數(shù)據(jù)樣本-因此,在權(quán)重參數(shù)上傳過程中,即使遇到模型參數(shù)竊取行為,竊取者也無法推導(dǎo)出該參與方對應(yīng)的所有訓(xùn)練數(shù)據(jù)集樣本.

    同時(shí),為了服務(wù)于模型聚合階段的權(quán)值計(jì)算過程,在此階段還需上傳各本地負(fù)荷運(yùn)營商lk所擁有數(shù)據(jù)Vl,k的時(shí)間跨度信息Uk.

    2.4模型聚合

    當(dāng)所有本地負(fù)荷運(yùn)營商lk完成模型上傳流程后,服務(wù)器端對接收到的全部模型執(zhí)行聚合.以下介紹本文執(zhí)行模型聚合步驟時(shí)權(quán)值wB,k的計(jì)算方法.

    現(xiàn)有的聯(lián)邦學(xué)習(xí)模型聚合方法并不適應(yīng)負(fù)荷預(yù)測場景.例如,在經(jīng)典聯(lián)邦學(xué)習(xí)方法中,某聯(lián)邦參與方lk模型聚合的權(quán)值wB,k僅與其參與訓(xùn)練的樣本數(shù)量大小有關(guān).這種方法沒有結(jié)合實(shí)際的應(yīng)用場景,未充分考慮各個(gè)聯(lián)邦參與方數(shù)據(jù)樣本的質(zhì)量.這會導(dǎo)致在聚合時(shí)使用了不合理的聚合權(quán)重參數(shù),影響最終模型的預(yù)測準(zhǔn)確率和泛化能力.而在本文方法所針對的短期負(fù)荷預(yù)測場景中,擁有越多樣化、越全面的數(shù)據(jù)的負(fù)荷運(yùn)營商所訓(xùn)練出的本地負(fù)荷預(yù)測模型往往有著更優(yōu)秀的預(yù)測準(zhǔn)確率與泛化能力,在模型聚合時(shí)應(yīng)該給予更高的聚合權(quán)重.

    本文結(jié)合負(fù)荷預(yù)測場景的特點(diǎn),定義各個(gè)負(fù)荷運(yùn)營商的貢獻(xiàn)度G,進(jìn)而更新聚合權(quán)值wB,k.貢獻(xiàn)度Gk與運(yùn)營商lk所擁有的數(shù)據(jù)樣本數(shù)量大小以及對應(yīng)的數(shù)據(jù)日期跨度Uk均有關(guān);當(dāng)一個(gè)負(fù)荷運(yùn)營商自身數(shù)據(jù)樣本數(shù)量越大或其所擁有的相對其他運(yùn)營商數(shù)據(jù)的新時(shí)刻數(shù)據(jù)越多,其貢獻(xiàn)度G應(yīng)該越大;該模型對服務(wù)器模型的聚合過程也應(yīng)該有更大的影響,參數(shù)聚合時(shí)擁有較高的權(quán)重.

    對每一個(gè)存在數(shù)據(jù)的時(shí)刻賦予總量為1的貢獻(xiàn)量,對于i時(shí)刻,若在該時(shí)刻共有Qi個(gè)負(fù)荷運(yùn)營商存在數(shù)據(jù)樣本,則共同平分該時(shí)刻的貢獻(xiàn)量,每個(gè)運(yùn)營商得到的貢獻(xiàn)分量為Bk,i;若負(fù)荷運(yùn)營商lk在該時(shí)刻i沒有數(shù)據(jù)樣本,則其貢獻(xiàn)分量Bk,i=0.

    貢獻(xiàn)度確定后,按照各參與方對聯(lián)邦模型的貢獻(xiàn)度Gk更新wB,k,執(zhí)行模型聚合,更新模型參數(shù).更新方法如式(17)(18)所示.

    式中:W為已更新后模型權(quán)重參數(shù);Wk為參與方lk在經(jīng)過本地訓(xùn)練后所上傳的模型權(quán)重參數(shù).

    3仿真驗(yàn)證

    3.1試驗(yàn)環(huán)境

    本文采用多個(gè)相互獨(dú)立的計(jì)算機(jī)模擬參與聯(lián)邦的各個(gè)負(fù)荷運(yùn)營商,進(jìn)而驗(yàn)證本文所提出方法對預(yù)測準(zhǔn)確率與泛化能力的提升.采用5臺計(jì)算機(jī)模擬5個(gè)本地負(fù)荷運(yùn)營商(地區(qū)),各計(jì)算機(jī)配置均為英特爾至強(qiáng)6248R CPU,RTX 3080 GPU,512 GB內(nèi)存;同時(shí)采用一臺個(gè)人計(jì)算機(jī)作為服務(wù)器端,配置為Intel i7 CPU,16 GB 內(nèi)存,Intel UHD Graphics 630集成顯卡.上述設(shè)備之間可點(diǎn)對點(diǎn)通信.

    本文采用在負(fù)荷預(yù)測中具有代表性的GEF- Com2012比賽數(shù)據(jù)[26],作為試驗(yàn)樣本數(shù)據(jù)集.該數(shù)據(jù)集包含20個(gè)地區(qū)每1h記錄1次的電力負(fù)荷數(shù)據(jù)、對應(yīng)溫度數(shù)據(jù)以及日期信息,將其整理成可用于本文LSTM模型輸入的數(shù)據(jù)集.深度學(xué)習(xí)算法采用Python 3.6編程語言編程,負(fù)荷預(yù)測模型的深度學(xué)習(xí)框架采用TensorFlow2.

    3.2仿真結(jié)果

    1)為了驗(yàn)證本文方法對模型預(yù)測能力的提升,同時(shí)排除日期等因素對預(yù)測精度的影響,在3.2.1節(jié)中將參與協(xié)同訓(xùn)練的各地區(qū)所擁有數(shù)據(jù)的時(shí)間跨度均設(shè)置為同一時(shí)間段,樣本設(shè)為相同數(shù)量.

    2)為了驗(yàn)證本文方法對模型泛化能力的提升,在3.2.2節(jié)中將各參與方的數(shù)據(jù)時(shí)間跨度設(shè)置為不同時(shí)間段,樣本設(shè)為不同數(shù)量.設(shè)置3個(gè)算例驗(yàn)證其在3種場景下對模型預(yù)測精度和泛化能力的提升.

    3.2.1相同時(shí)間跨度與數(shù)據(jù)量的參與方進(jìn)行聯(lián)邦

    為了體現(xiàn)本文方法對模型預(yù)測準(zhǔn)確率的提升,選取5個(gè)地區(qū)作為相互獨(dú)立的負(fù)荷運(yùn)營商參與聯(lián)邦過程,分別記為l1、l2、l3、l4、l5.數(shù)據(jù)集選擇該5個(gè)地區(qū)2004年量測數(shù)據(jù)Vl,1~Vl,5,每1h有1次記錄,全年每個(gè)Vl,k有8 784個(gè)數(shù)據(jù)樣本.

    需要說明的是,本算例主要目的是驗(yàn)證本文方法對預(yù)測精度的提升,所以應(yīng)排除各地區(qū)數(shù)據(jù)量及選取日期差異的干擾.因此,本算例采用1月1日至7月27日的數(shù)據(jù)樣本作為模型訓(xùn)練樣本集,以7月28日的24個(gè)數(shù)據(jù)樣本作為測試集.

    依次將5個(gè)地區(qū)作為目標(biāo)對象,每一個(gè)地區(qū)均進(jìn)行以下工作:首先以地區(qū)lk(k=1、2、3、4、5)自身所擁有的數(shù)據(jù)集進(jìn)行模型訓(xùn)練,并將該訓(xùn)練好的模型作為服務(wù)器模型,記錄此時(shí)模型對該日負(fù)荷的預(yù)測結(jié)果;隨后,將該服務(wù)器模型下發(fā)至除lk外其余4個(gè)本地端,各lk根據(jù)自身所擁有的數(shù)據(jù)集進(jìn)行本地模型訓(xùn)練;依次進(jìn)行本地模型更新、本地模型上傳、服務(wù)器接收模型、服務(wù)器模型聚合更新流程;地區(qū)lk得到更新后的服務(wù)器模型,并在測試集進(jìn)行測試.最終,記錄5個(gè)地區(qū)單獨(dú)模型負(fù)荷預(yù)測及聯(lián)邦后預(yù)測的結(jié)果,并對比實(shí)際負(fù)荷、單一模型及本文方法的預(yù)測結(jié)果.各地區(qū)負(fù)荷預(yù)測結(jié)果對比分別如圖4~圖8所示,RMSE和MAPE結(jié)果對比分別如表1和表2所示.

    從圖4~圖8可以看出,使用單一負(fù)荷預(yù)測模型預(yù)測的負(fù)荷曲線與實(shí)際負(fù)荷曲線存在著明顯差距,但經(jīng)過本文方法訓(xùn)練的模型所預(yù)測的負(fù)荷曲線與真實(shí)曲線之間的差距明顯縮小.由表1和表2可知,在5次測試中,各地區(qū)單一模型的預(yù)測結(jié)果相對該地區(qū)實(shí)際負(fù)荷的誤差均較大,地區(qū)4的MAPE甚至超過了5%,這表明單獨(dú)使用單一負(fù)荷預(yù)測模型并不具備良好的誤差控制效果.使用本文所提方法,模型的誤差指標(biāo)RMSE與MAPE都有了明顯下降.地區(qū)1的誤差下降最為明顯,相比單一模型,RMSE與MAPE分別減少了32%和43%.由此可見,本文所提方法在保護(hù)數(shù)據(jù)隱私的情況下,豐富了參與到模型訓(xùn)練的數(shù)據(jù)特征,彌補(bǔ)了單一地區(qū)數(shù)據(jù)集特征不夠全面的缺陷,最終驗(yàn)證了經(jīng)過本文所提方法訓(xùn)練的負(fù)荷預(yù)測模型具有比單一數(shù)據(jù)集訓(xùn)練的模型更加優(yōu)秀的預(yù)測準(zhǔn)確度和更好的誤差控制能力.

    3.2.2多場景測試

    本文所提方法也適合更多復(fù)雜的場景.為了驗(yàn)證經(jīng)過本文方法訓(xùn)練后的模型在多場景下具有良好的預(yù)測準(zhǔn)確率和泛化能力,設(shè)置如下3個(gè)算例:算例1和算例2模擬利用本文方法實(shí)現(xiàn)擴(kuò)充數(shù)據(jù)集時(shí)間跨度的場景;算例3模擬將模型移植到新的地區(qū)進(jìn)行預(yù)測的場景.

    對仿真所用到的數(shù)據(jù)集進(jìn)行預(yù)處理,地區(qū)1~地區(qū)5的數(shù)據(jù)量及日期特征均不相同.其中地區(qū)4由于為新設(shè)立的區(qū)域,并無歷史數(shù)據(jù),但具有負(fù)荷、氣象等數(shù)據(jù)量測能力.詳細(xì)設(shè)置如表3所示.

    算例1以地區(qū)2的數(shù)據(jù)進(jìn)行訓(xùn)練,得到地區(qū)2 負(fù)荷預(yù)測模型,并用該模型預(yù)測地區(qū)2的負(fù)荷,記錄預(yù)測結(jié)果;再以地區(qū)1聯(lián)合地區(qū)2進(jìn)行協(xié)同訓(xùn)練,用訓(xùn)練出的模型預(yù)測地區(qū)2的負(fù)荷,記錄此時(shí)預(yù)測結(jié)果.算例1負(fù)荷預(yù)測結(jié)果對比如圖9所示,預(yù)測誤差對比如表4所示.

    由圖9和表4可知,采用本文所提方法后,在保證地區(qū)1數(shù)據(jù)隱私的同時(shí),訓(xùn)練出的模型降低了地區(qū)2因7月數(shù)據(jù)缺失造成7月份負(fù)荷預(yù)測準(zhǔn)確度下降的影響,預(yù)測出的負(fù)荷曲線更接近實(shí)際數(shù)據(jù),其誤差也大幅降低,RMSE與MAPE分別降低了21%和28%.這體現(xiàn)出本文方法有效地豐富了地區(qū)5模型訓(xùn)練可利用的樣本特征,驗(yàn)證了本文方法在數(shù)據(jù)量不足情況下對模型預(yù)測準(zhǔn)確率的提升效果.

    算例2以地區(qū)5的數(shù)據(jù)進(jìn)行訓(xùn)練,得到地區(qū)5負(fù)荷預(yù)測模型,并用該模型預(yù)測地區(qū)5的負(fù)荷,記錄預(yù)測結(jié)果;再以地區(qū)3聯(lián)合地區(qū)5進(jìn)行協(xié)同訓(xùn)練,用訓(xùn)練出的模型預(yù)測地區(qū)5的負(fù)荷,記錄此時(shí)預(yù)測結(jié)果.結(jié)果如表5所示,預(yù)測負(fù)荷曲線圖如圖10所示.

    從表5和圖10可以看出,相較地區(qū)5僅用自身數(shù)據(jù)訓(xùn)練的模型,采用本文方法訓(xùn)練后的模型預(yù)測曲線與真實(shí)曲線之間的差距大幅縮小,對應(yīng)的誤差RMSE和MAPE分別下降了74%和80%.在本算例中,地區(qū)5自身所擁有的數(shù)據(jù)極少,僅用該數(shù)據(jù)訓(xùn)練出的模型預(yù)測能力和泛化能力非常弱.由于本文方法對數(shù)據(jù)量進(jìn)行了大幅擴(kuò)充,且補(bǔ)充了更多不同日期的特征,對預(yù)測準(zhǔn)確率的提升十分顯著,使得誤差大幅降低.這也符合深度學(xué)習(xí)在樣本不夠全面的情況下,當(dāng)訓(xùn)練樣本增加時(shí)模型性能也提升的特點(diǎn).與此同時(shí),在協(xié)同訓(xùn)練過程中有效地保證了地區(qū)3、5 各自的數(shù)據(jù)隱私.綜上所述,本算例說明本文方法能有效提升模型的預(yù)測能力.

    算例3利用地區(qū)1、2、3、5的數(shù)據(jù)協(xié)同訓(xùn)練負(fù)荷預(yù)測模型,將該模型應(yīng)用于新加入聯(lián)邦的地區(qū)4的負(fù)荷預(yù)測;先量測地區(qū)4預(yù)測時(shí)刻前24 h時(shí)刻的歷史負(fù)荷數(shù)據(jù)以及預(yù)測時(shí)刻的氣象、日期數(shù)據(jù),作為聯(lián)合預(yù)測模型的輸入,進(jìn)而預(yù)測下一時(shí)刻的負(fù)荷變化,記錄下一日預(yù)測結(jié)果.預(yù)測負(fù)荷曲線圖如圖11所示,預(yù)測誤差對比如表6所示.從圖11和表6中可以看出,模型在新地區(qū)具有良好的預(yù)測準(zhǔn)確率,相比單一模型,本文方法預(yù)測曲線與真實(shí)負(fù)荷曲線差距較小,RMSE和MAPE分別降低了19%和23%.這是因?yàn)楸疚姆椒ㄓ?xùn)練出的短期負(fù)荷預(yù)測模型學(xué)習(xí)了多個(gè)地區(qū)數(shù)據(jù)的特征,適應(yīng)能力更強(qiáng),具有良好的泛化能力.

    4結(jié)論

    本文針對構(gòu)建負(fù)荷預(yù)測模型所需數(shù)據(jù)樣本數(shù)量與特征不足、模型泛化能力差以及負(fù)荷相關(guān)數(shù)據(jù)高度保密無法共享訓(xùn)練的問題,引入聯(lián)邦學(xué)習(xí)方法,提出基于聯(lián)邦學(xué)習(xí)的短期負(fù)荷預(yù)測模型協(xié)同訓(xùn)練方法,實(shí)現(xiàn)各負(fù)荷運(yùn)營商在滿足保證數(shù)據(jù)隱私的前提下協(xié)同訓(xùn)練短期負(fù)荷預(yù)測模型.在GEFCom2012比賽的多個(gè)地區(qū)負(fù)荷數(shù)據(jù)集上進(jìn)行了仿真測試,驗(yàn)證了所提方法對模型預(yù)測準(zhǔn)確率的提升,證明了通過本文方法訓(xùn)練出的短期負(fù)荷預(yù)測模型在多場景下具有優(yōu)秀的泛化能力.在本文基礎(chǔ)上,未來可進(jìn)一步研究配電網(wǎng)多主體協(xié)同的去中心化分布式訓(xùn)練機(jī)制.

    參考文獻(xiàn)

    [1]梁智,孫國強(qiáng),李虎成,等.基于VMD與PSO優(yōu)化深度信念網(wǎng)絡(luò)的短期負(fù)荷預(yù)測[J].電網(wǎng)技術(shù),2018,42(2):598-606.

    LIAN G Z,SUN G Q,LI H C,et al. Short-term load forecasting based on VMD and PSO optimized deep belief network [J]. Power System Technology,2018,42(2):598-606. (In Chinese)

    [2]陳宏義,李存斌,施立剛.基于聚類分析的短期負(fù)荷智能預(yù)測方法研究[J].湖南大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,41(5):94-98.

    CHEN H Y,LI C B,SHI L G. A new forecasting approach for short-term load intelligence based on cluster method[J]. Journal of Hunan University (Natural Sciences),2014,41(5):94-98. (In Chinese)

    [3]孔祥玉,李闖,鄭鋒,等.基于經(jīng)驗(yàn)?zāi)B(tài)分解與特征相關(guān)分析的短期負(fù)荷預(yù)測方法[J].電力系統(tǒng)自動化,2019,43(5):46-52.

    KONG X Y,LI C,ZHENG F,et al. Short-term load forecasting method based on empirical mode decomposition and feature correlation analysis[J]. Automation of Electric Power Systems,2019,43(5):46-52.(In Chinese)

    [4]彭顯剛,鄭偉欽,林利祥,等.考慮負(fù)荷自適應(yīng)檢測和修復(fù)的魯棒極限學(xué)習(xí)機(jī)短期負(fù)荷預(yù)測方法[J].中國電機(jī)工程學(xué)報(bào),2016,36(23):6409-6417.

    PENG X G,ZHENG W Q,LIN L X,et al. Short-term load forecasting method based on outlier robust extreme learning machine considering adaptive load detection and repair[J].Proceedings of the CSEE,2016,36(23):6409-6417.(In Chinese)

    [5]孔祥玉,鄭鋒,鄂志君,等.基于深度信念網(wǎng)絡(luò)的短期負(fù)荷預(yù)測方法[J].電力系統(tǒng)自動化,2018,42(5):133-139.

    KONG X Y,ZHENG F,E Z J,et al. Short-term load forecasting based on deep belief network[J]. Automation of Electric Power Systems,2018,42(5):133-139.(In Chinese)

    [6]彭文,王金睿,尹山青.電力市場中基于Attention-LSTM的短期負(fù)荷預(yù)測模型[J].電網(wǎng)技術(shù),2019,43(5):1745-1751.

    PENG W,WANG J R,YIN S Q. Short-term load forecasting model based on Attention-LSTM in electricity market[J]. Power System Technology,2019,43(5):1745-1751.(In Chinese)

    [7]蘇學(xué)能,劉天琪,曹鴻謙,等.基于Hadoop架構(gòu)的多重分布式BP神經(jīng)網(wǎng)絡(luò)的短期負(fù)荷預(yù)測方法[J].中國電機(jī)工程學(xué)報(bào),2017,37(17):4966-4973.

    SU X N,LIU T Q,CAO H Q,et al.A multiple distributed BP neural networks approach for short-term load forecasting based on Hadoop framework[J]. Proceedings of the CSEE,2017,37(17):4966-4973.(In Chinese)

    [8]楊龍,吳紅斌,丁明,等.新能源電網(wǎng)中考慮特征選擇的Bi-LSTM網(wǎng)絡(luò)短期負(fù)荷預(yù)測[J].電力系統(tǒng)自動化,2021,45 (3):166-173.

    YANG L,WU H B,DING M,et al. Short-term load forecasting in renewable energy grid based on Bi-directional long short-term memory network considering feature selection[J]. Automation of Electric Power Systems,2021,45(3):166-173.(In Chinese)

    [9]龐昊,高金峰,杜耀恒.基于多神經(jīng)網(wǎng)絡(luò)融合的短期負(fù)荷預(yù)測方法[J].電力自動化設(shè)備,2020,40力):37-43.

    PANG H,GAO J F,DU Y H.Short-term load forecasting method based on fusion of multiple neural networks[J]. Electric Power Automation Equipment,2020,40(6):37-43. (In Chinese)

    [10]姚程文,楊蘋,劉澤健.基于CNN-GRU混合神經(jīng)網(wǎng)絡(luò)的負(fù)荷預(yù)測方法[J].電網(wǎng)技術(shù),2020,44(9):3416-3424.

    YAO C W,YANG P,LIU Z J.Load forecasting method based on CNN-GRU hybrid neural network[J]. Power System Technology,2020,44(9):3416-3424.(In Chinese)

    [11]陸繼翔,張琪培,楊志宏,等.基于CNN-LSTM混合神經(jīng)網(wǎng)絡(luò)模型的短期負(fù)荷預(yù)測方法[J].電力系統(tǒng)自動化,2019,43(8):131-137.

    LU J X,ZHANG Q P,YANG Z H,et al.Short-term load forecasting method based on CNN-LSTM hybrid neural network model[J]. Automation of Electric Power Systems,2019,43(8):131-137.(In Chinese)

    [12]鄧帶雨,李堅(jiān),張真源,等.基于EEMD-GRU-MLR的短期電力負(fù)荷預(yù)測[J].電網(wǎng)技術(shù),2020,44(2):593-602.

    DENG D Y,LI J,ZHANG Z Y,et al. Short-term electric load forecasting based on EEMD-GRU-MLR[J].Power System Technology,2020,44(2):593-602.(In Chinese)

    [13]王增平,趙兵,紀(jì)維佳,等.基于GRU-NN模型的短期負(fù)荷預(yù)測方法[J].電力系統(tǒng)自動化,2019,43(5):53-58.

    WANG Z P,ZHAO B,JI W J,et al. Short-term load forecasting method based on GRU-NN model [J]. Automation of Electric Power Systems,2019,43(5):53-58. (In Chinese)

    [14]陳振宇,劉金波,李晨,等.基于LSTM與XGBoost組合模型的超短期電力負(fù)荷預(yù)測[J].電網(wǎng)技術(shù),2020,44(2):614-620.

    CHEN Z Y,LIU J B,LI C,et al. Ultra short-term power load forecasting based on combined LSTM-XGBoost model [J]. Power System Technology,2020,44(2):614-620. (In Chinese)

    [15]孫曉燕,李家釗,曾博,等.基于特征遷移學(xué)習(xí)的綜合能源系統(tǒng)小樣本日前電力負(fù)荷預(yù)測[J].控制理論與應(yīng)用,2021,38(1):63-72.

    SUN X Y,LI J Z,ZENG B,et al Small-sample day-ahead power load forecasting of integrated energy system based on feature transfer learning [J]. Control Theory & Applications,2021,38(1):63-72 (In Chinese)

    [16]張勇,陶一凡,鞏敦衛(wèi).遷移學(xué)習(xí)引導(dǎo)的變源域長短時(shí)記憶網(wǎng)絡(luò)建筑負(fù)荷預(yù)測[J].控制與決策,2021,36(10):2328-2338.

    ZHANG Y,TAO Y F,GONG D W. Load forecasting of buildings using LSTM based on transfer learning with variable source domain[J] Control and Decision,2021,36(10):2328-2338 (In Chinese)

    [17] GONQALVES C,BESSA R J,PINSON P. Privacy-preserving distributed learning for renewable energy forecasting[J]. IEEE Transactions on Sustainable Energy,2021,12(3):1777-1787

    [18] KONECNY J,MCMAHAN H B,YUF X,et al. Federated learning:strategies for improving communication efficiency [EB/OL]. https://arxiv.org/pdf/1610.05492.pdf.

    [19] YANG Q,LIU Y,CHEN T J,et al. Federated machine learning[J]. ACM Transactions on Intelligent Systems and Technology,2019,10(2):1-19

    [20]賀文晨,郭少勇,邱雪松,等.基于DRL的聯(lián)邦學(xué)習(xí)節(jié)點(diǎn)選擇方法[J].通信學(xué)報(bào),2021,42(6):62-71.

    HE W C,GUO S Y,QIU X S,et al. Node selection method in federated learning based on deep reinforcement learning[J]. Journal on Communications,2021,42(6):62-71 (In Chinese)

    [21] MILLS J,HU J,MIN G Y Communication-efficient federated learning for wireless edge intelligence in IoT[J]. IEEE Internet of Things Journal,2020,7(7):5986-5994

    [22] ZHOU C Y,F(xiàn)U A M,YU S,et al Privacy-preserving federated learning in fog computing[J] IEEE Internet of Things Journal,2020,7(11):10782-10793

    [23]鄭楷洪,肖勇,王鑫,等.一個(gè)面向電力計(jì)量系統(tǒng)的聯(lián)邦學(xué)習(xí)框架[J].中國電機(jī)工程學(xué)報(bào),2020,40(S1):122-133.

    ZHENG K H,XIAO Y,WANG X,et al A federated learning framework for power grid metering system [J]. Proceedings of the CSEE,2020,40(S1):122-133 (In Chinese)

    [24] HOCHREITER S,SCHMIDHUBER J. Long short-term memory[J]. Neural Computation,1997,9(8):1735-1780

    [25] ZHU L G,LIU Z J,HAN S. Deep leakage from gradients[EB/OL]. https://arxiv.org/pdf/1906.08935.pdf.

    [26] HONG T,PINSON P,F(xiàn)AN S Global energy forecasting competition 2012[J]. International Journal of Forecasting,2014,30(2):357-363

    猜你喜歡
    負(fù)荷預(yù)測機(jī)器學(xué)習(xí)
    基于高階BP神經(jīng)網(wǎng)絡(luò)的日最大負(fù)荷預(yù)測
    電力系統(tǒng)短期負(fù)荷研究綜述
    基于詞典與機(jī)器學(xué)習(xí)的中文微博情感分析
    電力系統(tǒng)負(fù)荷的一元線性模型預(yù)測方法
    面向云計(jì)算的智能電網(wǎng)負(fù)荷預(yù)測平臺模型研究
    基于電力系統(tǒng)負(fù)荷預(yù)測方法的配電網(wǎng)規(guī)劃中的實(shí)踐研究
    基于機(jī)器學(xué)習(xí)的圖像特征提取技術(shù)在圖像版權(quán)保護(hù)中的應(yīng)用
    基于網(wǎng)絡(luò)搜索數(shù)據(jù)的平遙旅游客流量預(yù)測分析
    前綴字母為特征在維吾爾語文本情感分類中的研究
    基于支持向量機(jī)的金融數(shù)據(jù)分析研究
    亚洲电影在线观看av| 亚洲欧美日韩高清在线视频| 好看av亚洲va欧美ⅴa在| 亚洲 欧美 日韩 在线 免费| 亚洲激情在线av| 最近最新中文字幕大全免费视频| 久久久久性生活片| 可以在线观看毛片的网站| 曰老女人黄片| 少妇丰满av| 黄色视频,在线免费观看| 久久性视频一级片| 亚洲第一电影网av| 在线观看66精品国产| 美女大奶头视频| 变态另类丝袜制服| 亚洲成人久久爱视频| 日本一本二区三区精品| 亚洲熟女毛片儿| 美女午夜性视频免费| 琪琪午夜伦伦电影理论片6080| 十八禁人妻一区二区| 国产精品久久电影中文字幕| 国产乱人伦免费视频| 精华霜和精华液先用哪个| 露出奶头的视频| 久久久国产成人免费| 久久热在线av| 国产精品1区2区在线观看.| 国产精品久久久av美女十八| 51午夜福利影视在线观看| 国产麻豆成人av免费视频| 国产在线精品亚洲第一网站| 身体一侧抽搐| 国产精品九九99| 毛片女人毛片| 人妻丰满熟妇av一区二区三区| 国产精品一区二区精品视频观看| 99久久精品热视频| 99riav亚洲国产免费| 亚洲国产精品久久男人天堂| 日本成人三级电影网站| 日韩欧美三级三区| 精品熟女少妇八av免费久了| 精品久久久久久久人妻蜜臀av| 999久久久精品免费观看国产| 久久久成人免费电影| av视频在线观看入口| 又爽又黄无遮挡网站| 午夜日韩欧美国产| 亚洲成人精品中文字幕电影| 51午夜福利影视在线观看| 亚洲精品456在线播放app | 国产主播在线观看一区二区| 嫩草影院入口| 久久午夜综合久久蜜桃| 亚洲中文字幕一区二区三区有码在线看 | 亚洲黑人精品在线| 特大巨黑吊av在线直播| 熟妇人妻久久中文字幕3abv| 国产野战对白在线观看| 日韩欧美免费精品| 老司机福利观看| www日本在线高清视频| 久久久国产欧美日韩av| 狠狠狠狠99中文字幕| 老司机午夜福利在线观看视频| 91九色精品人成在线观看| av天堂在线播放| 久久久国产精品麻豆| 国产麻豆成人av免费视频| 免费一级毛片在线播放高清视频| 亚洲真实伦在线观看| 亚洲无线在线观看| 国产97色在线日韩免费| h日本视频在线播放| 久久精品aⅴ一区二区三区四区| 两个人视频免费观看高清| 亚洲一区高清亚洲精品| 亚洲午夜精品一区,二区,三区| 亚洲成人中文字幕在线播放| 一边摸一边抽搐一进一小说| 天堂网av新在线| 人妻夜夜爽99麻豆av| 亚洲熟妇熟女久久| 99久久精品国产亚洲精品| 国产精品一区二区精品视频观看| 宅男免费午夜| 99精品久久久久人妻精品| 脱女人内裤的视频| 婷婷精品国产亚洲av在线| 亚洲av电影在线进入| 美女 人体艺术 gogo| 真人做人爱边吃奶动态| 色在线成人网| 欧美黑人巨大hd| 日韩欧美在线二视频| 国产黄片美女视频| 99热只有精品国产| 精品无人区乱码1区二区| 看黄色毛片网站| 好男人电影高清在线观看| 亚洲精品国产精品久久久不卡| av女优亚洲男人天堂 | 黄色视频,在线免费观看| 最新在线观看一区二区三区| 亚洲色图av天堂| 国产高清激情床上av| a级毛片在线看网站| 成人18禁在线播放| 麻豆久久精品国产亚洲av| 国产一区二区在线av高清观看| 亚洲人成电影免费在线| 看片在线看免费视频| 又粗又爽又猛毛片免费看| 久久久国产精品麻豆| 桃色一区二区三区在线观看| 欧美日韩精品网址| 深夜精品福利| 亚洲人成伊人成综合网2020| 欧美激情久久久久久爽电影| 精品乱码久久久久久99久播| 成年人黄色毛片网站| 人妻丰满熟妇av一区二区三区| 日韩免费av在线播放| 国产亚洲欧美在线一区二区| 久久草成人影院| 欧美不卡视频在线免费观看| 国产私拍福利视频在线观看| 欧美日韩一级在线毛片| 国产亚洲精品av在线| 村上凉子中文字幕在线| 久久久精品欧美日韩精品| 欧美色欧美亚洲另类二区| 久久精品国产综合久久久| 欧美+亚洲+日韩+国产| 久久精品综合一区二区三区| x7x7x7水蜜桃| 搡老熟女国产l中国老女人| 人人妻,人人澡人人爽秒播| 成在线人永久免费视频| 成年女人毛片免费观看观看9| 免费看光身美女| 国产真人三级小视频在线观看| 少妇裸体淫交视频免费看高清| 免费在线观看亚洲国产| 又黄又粗又硬又大视频| 少妇人妻一区二区三区视频| 啪啪无遮挡十八禁网站| 中文字幕最新亚洲高清| 九九在线视频观看精品| 国产精品一区二区三区四区久久| av女优亚洲男人天堂 | 欧美日韩中文字幕国产精品一区二区三区| 美女高潮的动态| 欧美日韩国产亚洲二区| 久久久成人免费电影| 亚洲精品久久国产高清桃花| 国内精品久久久久精免费| 日韩欧美 国产精品| 亚洲国产欧美人成| 99久国产av精品| 亚洲欧美日韩东京热| 亚洲成人中文字幕在线播放| 99riav亚洲国产免费| 精品久久久久久久久久免费视频| 老汉色∧v一级毛片| 欧美日韩中文字幕国产精品一区二区三区| 校园春色视频在线观看| 熟女电影av网| 国产午夜精品久久久久久| 精品一区二区三区四区五区乱码| 成人高潮视频无遮挡免费网站| 黄片小视频在线播放| 精品一区二区三区四区五区乱码| 网址你懂的国产日韩在线| 中文字幕人妻丝袜一区二区| 午夜福利成人在线免费观看| 一区二区三区高清视频在线| 两个人视频免费观看高清| 亚洲av成人精品一区久久| 午夜福利在线在线| 丝袜人妻中文字幕| 亚洲成av人片在线播放无| avwww免费| 美女cb高潮喷水在线观看 | 久久久久久久精品吃奶| 高潮久久久久久久久久久不卡| 午夜免费成人在线视频| 夜夜爽天天搞| 国产又黄又爽又无遮挡在线| 少妇人妻一区二区三区视频| 美女cb高潮喷水在线观看 | 日本 欧美在线| 啦啦啦观看免费观看视频高清| 欧美又色又爽又黄视频| 这个男人来自地球电影免费观看| 久久久久免费精品人妻一区二区| 99热只有精品国产| 国产亚洲欧美98| 一级黄色大片毛片| 最近视频中文字幕2019在线8| 欧美极品一区二区三区四区| 成人av一区二区三区在线看| 亚洲精华国产精华精| 美女被艹到高潮喷水动态| 午夜福利18| 欧美xxxx黑人xx丫x性爽| 国产av麻豆久久久久久久| 在线视频色国产色| 国产淫片久久久久久久久 | 女人被狂操c到高潮| 久久香蕉精品热| 亚洲狠狠婷婷综合久久图片| 色尼玛亚洲综合影院| 91字幕亚洲| 精品免费久久久久久久清纯| 精品久久久久久久人妻蜜臀av| 草草在线视频免费看| 久9热在线精品视频| 国产成+人综合+亚洲专区| 国产欧美日韩一区二区三| x7x7x7水蜜桃| 麻豆久久精品国产亚洲av| 一夜夜www| 九色国产91popny在线| 97碰自拍视频| 国产精品 国内视频| 波多野结衣巨乳人妻| 在线观看一区二区三区| 国产69精品久久久久777片 | 国产成人欧美在线观看| 国产精品久久久久久久电影 | 日本黄色片子视频| 三级国产精品欧美在线观看 | 久久精品91蜜桃| 91av网站免费观看| 日日夜夜操网爽| 97人妻精品一区二区三区麻豆| 一本久久中文字幕| 亚洲一区二区三区不卡视频| 久久久水蜜桃国产精品网| 成人亚洲精品av一区二区| 不卡一级毛片| 一个人免费在线观看电影 | 99国产综合亚洲精品| www.熟女人妻精品国产| 啦啦啦免费观看视频1| 青草久久国产| 中文字幕高清在线视频| 色综合亚洲欧美另类图片| 国产日本99.免费观看| 国产视频一区二区在线看| 国产一区二区三区在线臀色熟女| 天堂网av新在线| 日韩 欧美 亚洲 中文字幕| 久久久国产成人免费| 欧美不卡视频在线免费观看| 欧美性猛交╳xxx乱大交人| 亚洲人与动物交配视频| 欧美不卡视频在线免费观看| 噜噜噜噜噜久久久久久91| 国产美女午夜福利| 又黄又粗又硬又大视频| 国产精品爽爽va在线观看网站| 高清在线国产一区| 国产成人系列免费观看| 午夜激情欧美在线| 成人特级黄色片久久久久久久| 国产精品女同一区二区软件 | 噜噜噜噜噜久久久久久91| 午夜久久久久精精品| 少妇人妻一区二区三区视频| 国产成人aa在线观看| 黑人巨大精品欧美一区二区mp4| 香蕉丝袜av| 天堂网av新在线| avwww免费| 免费在线观看成人毛片| 日韩免费av在线播放| 久久久久国内视频| 久久久久国产精品人妻aⅴ院| 噜噜噜噜噜久久久久久91| 中文字幕久久专区| 波多野结衣高清作品| 国产av一区在线观看免费| 村上凉子中文字幕在线| 亚洲人成伊人成综合网2020| 午夜福利高清视频| 最近在线观看免费完整版| 国产精品一区二区精品视频观看| av视频在线观看入口| 人人妻人人看人人澡| 男女床上黄色一级片免费看| 特大巨黑吊av在线直播| 少妇人妻一区二区三区视频| 两人在一起打扑克的视频| 亚洲成人久久性| 18美女黄网站色大片免费观看| 熟妇人妻久久中文字幕3abv| av在线天堂中文字幕| 国产三级在线视频| 熟女人妻精品中文字幕| 亚洲成人久久性| 99riav亚洲国产免费| 岛国在线观看网站| 在线观看美女被高潮喷水网站 | 一本一本综合久久| 日本黄色视频三级网站网址| 亚洲五月天丁香| 91九色精品人成在线观看| 免费高清视频大片| 真实男女啪啪啪动态图| 一进一出好大好爽视频| 亚洲国产中文字幕在线视频| 国产不卡一卡二| 国产日本99.免费观看| 12—13女人毛片做爰片一| 精品久久久久久久末码| 国产精品国产高清国产av| 99国产综合亚洲精品| 国产精品野战在线观看| 日韩三级视频一区二区三区| 淫妇啪啪啪对白视频| 色综合亚洲欧美另类图片| 老司机福利观看| 午夜影院日韩av| 国产精品美女特级片免费视频播放器 | 久久精品综合一区二区三区| 国产精品电影一区二区三区| 哪里可以看免费的av片| 最近最新中文字幕大全电影3| av欧美777| 少妇的丰满在线观看| 国产精品影院久久| 老司机午夜十八禁免费视频| 亚洲 国产 在线| 欧美色视频一区免费| 熟女电影av网| 黄片大片在线免费观看| 欧美日韩综合久久久久久 | 91九色精品人成在线观看| 久久精品影院6| 天堂网av新在线| 亚洲av美国av| 久久午夜综合久久蜜桃| 美女黄网站色视频| 真实男女啪啪啪动态图| 亚洲精品色激情综合| 久久国产精品人妻蜜桃| 亚洲av美国av| 人妻丰满熟妇av一区二区三区| 深夜精品福利| 啦啦啦韩国在线观看视频| 亚洲av成人av| 国产伦精品一区二区三区四那| 波多野结衣高清作品| 亚洲在线自拍视频| 国产成人福利小说| 免费一级毛片在线播放高清视频| 一个人免费在线观看的高清视频| 国产亚洲精品一区二区www| 九九热线精品视视频播放| 亚洲中文日韩欧美视频| 亚洲九九香蕉| 又黄又粗又硬又大视频| 欧美激情在线99| 在线十欧美十亚洲十日本专区| 51午夜福利影视在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产精品女同一区二区软件 | a在线观看视频网站| 国产亚洲av高清不卡| 国产成人影院久久av| 亚洲成人免费电影在线观看| 男人和女人高潮做爰伦理| 亚洲 欧美 日韩 在线 免费| 好男人在线观看高清免费视频| 黄色 视频免费看| 精品欧美国产一区二区三| 我的老师免费观看完整版| 男人和女人高潮做爰伦理| 老汉色∧v一级毛片| 一级作爱视频免费观看| 波多野结衣高清无吗| 午夜激情欧美在线| 亚洲欧美日韩无卡精品| 美女扒开内裤让男人捅视频| 国产91精品成人一区二区三区| 久久国产精品人妻蜜桃| 亚洲精品粉嫩美女一区| 中文资源天堂在线| 中文字幕最新亚洲高清| 丰满人妻一区二区三区视频av | 9191精品国产免费久久| 久久中文字幕一级| 婷婷亚洲欧美| 天天添夜夜摸| 国产探花在线观看一区二区| 日韩欧美 国产精品| 久久久精品大字幕| 亚洲激情在线av| 国产久久久一区二区三区| 亚洲国产看品久久| 亚洲精品久久国产高清桃花| 欧美高清成人免费视频www| 国产精品日韩av在线免费观看| 两个人看的免费小视频| 久久这里只有精品中国| xxxwww97欧美| 国产乱人视频| 这个男人来自地球电影免费观看| 亚洲专区字幕在线| 国产精品九九99| 啦啦啦观看免费观看视频高清| 男人舔女人下体高潮全视频| 女人被狂操c到高潮| 日韩精品青青久久久久久| 久久久色成人| 757午夜福利合集在线观看| 一二三四社区在线视频社区8| 国产人伦9x9x在线观看| 美女大奶头视频| 亚洲成a人片在线一区二区| 国产亚洲av高清不卡| 一个人免费在线观看的高清视频| 成人精品一区二区免费| 日韩欧美免费精品| av福利片在线观看| 久久午夜亚洲精品久久| 身体一侧抽搐| 夜夜夜夜夜久久久久| 色噜噜av男人的天堂激情| 日韩欧美国产在线观看| 在线免费观看的www视频| 两个人看的免费小视频| 19禁男女啪啪无遮挡网站| 日韩大尺度精品在线看网址| 亚洲精品色激情综合| 成在线人永久免费视频| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品美女久久av网站| 最近最新中文字幕大全电影3| 变态另类丝袜制服| 无限看片的www在线观看| 久久久久久久精品吃奶| 一个人看视频在线观看www免费 | 国产私拍福利视频在线观看| 99热精品在线国产| 亚洲自拍偷在线| 午夜福利成人在线免费观看| 亚洲av美国av| 亚洲国产精品久久男人天堂| 亚洲,欧美精品.| 首页视频小说图片口味搜索| 人人妻,人人澡人人爽秒播| 欧美成人免费av一区二区三区| 国产欧美日韩一区二区三| 精品国内亚洲2022精品成人| 国产精品免费一区二区三区在线| 小蜜桃在线观看免费完整版高清| 亚洲av熟女| 成人特级黄色片久久久久久久| 床上黄色一级片| 亚洲中文日韩欧美视频| 欧美激情在线99| 久久精品91蜜桃| 免费看十八禁软件| 好男人电影高清在线观看| 中出人妻视频一区二区| 国产精品久久久人人做人人爽| 成人亚洲精品av一区二区| 日韩欧美三级三区| 韩国av一区二区三区四区| 免费搜索国产男女视频| 亚洲国产精品合色在线| 久久久色成人| 在线视频色国产色| 亚洲天堂国产精品一区在线| 在线永久观看黄色视频| 精品久久久久久久人妻蜜臀av| а√天堂www在线а√下载| 91九色精品人成在线观看| 国产精品 国内视频| 每晚都被弄得嗷嗷叫到高潮| 天天躁日日操中文字幕| 精品一区二区三区视频在线观看免费| 狂野欧美白嫩少妇大欣赏| 在线十欧美十亚洲十日本专区| 男人的好看免费观看在线视频| 嫩草影视91久久| 久久精品91无色码中文字幕| 亚洲专区中文字幕在线| 亚洲欧美精品综合久久99| 伊人久久大香线蕉亚洲五| 久久草成人影院| 在线免费观看的www视频| 久久伊人香网站| 国产精品九九99| 久久久久国产精品人妻aⅴ院| 偷拍熟女少妇极品色| 黄色日韩在线| 亚洲人成电影免费在线| 日本免费一区二区三区高清不卡| 欧美丝袜亚洲另类 | 亚洲成av人片在线播放无| 熟妇人妻久久中文字幕3abv| 色综合亚洲欧美另类图片| 一卡2卡三卡四卡精品乱码亚洲| 国内精品美女久久久久久| 久久天躁狠狠躁夜夜2o2o| 看黄色毛片网站| 久久香蕉国产精品| 精品久久久久久成人av| 久久人人精品亚洲av| 免费人成视频x8x8入口观看| 亚洲国产欧美网| 欧美又色又爽又黄视频| 精品熟女少妇八av免费久了| 婷婷精品国产亚洲av在线| 我要搜黄色片| 一个人看视频在线观看www免费 | 国产精品精品国产色婷婷| 国产精品一区二区精品视频观看| 亚洲成人久久爱视频| 18禁观看日本| 欧美另类亚洲清纯唯美| 国产一区二区激情短视频| 手机成人av网站| 老鸭窝网址在线观看| 又粗又爽又猛毛片免费看| 国产精品 欧美亚洲| 国产精品九九99| 色综合婷婷激情| 18美女黄网站色大片免费观看| 男人舔奶头视频| 女警被强在线播放| 少妇丰满av| 看片在线看免费视频| 老鸭窝网址在线观看| 久久香蕉精品热| 色视频www国产| 88av欧美| 亚洲在线自拍视频| 岛国在线免费视频观看| 久久久久久人人人人人| 一卡2卡三卡四卡精品乱码亚洲| 国产单亲对白刺激| 色噜噜av男人的天堂激情| 99精品在免费线老司机午夜| 国产精品日韩av在线免费观看| h日本视频在线播放| 久久这里只有精品19| 老熟妇仑乱视频hdxx| 日韩欧美在线二视频| 午夜福利免费观看在线| 丁香六月欧美| 国产精品 国内视频| 人人妻,人人澡人人爽秒播| 一夜夜www| 999精品在线视频| 精品99又大又爽又粗少妇毛片 | 村上凉子中文字幕在线| 精品久久蜜臀av无| 老汉色av国产亚洲站长工具| 757午夜福利合集在线观看| 熟女人妻精品中文字幕| 国产av在哪里看| 国产在线精品亚洲第一网站| 日日干狠狠操夜夜爽| 国产精品一区二区免费欧美| 最近最新中文字幕大全免费视频| 日本免费一区二区三区高清不卡| 成人亚洲精品av一区二区| 一a级毛片在线观看| 久久久久久久久中文| 一本久久中文字幕| 12—13女人毛片做爰片一| 午夜成年电影在线免费观看| 最近视频中文字幕2019在线8| 小蜜桃在线观看免费完整版高清| 丝袜人妻中文字幕| 91av网一区二区| 成人一区二区视频在线观看| 久久久久久大精品| 一区二区三区激情视频| 特大巨黑吊av在线直播| 在线永久观看黄色视频| 又黄又粗又硬又大视频| 精品熟女少妇八av免费久了| 麻豆久久精品国产亚洲av| 国产麻豆成人av免费视频| 美女黄网站色视频| 欧美日韩瑟瑟在线播放| 成人无遮挡网站| 久久久国产欧美日韩av| 可以在线观看毛片的网站| 在线观看午夜福利视频| 国内精品美女久久久久久| 亚洲精品乱码久久久v下载方式 | 最新中文字幕久久久久 | 草草在线视频免费看| 成人亚洲精品av一区二区| 国产一区二区三区在线臀色熟女| 国产99白浆流出| 午夜成年电影在线免费观看| 亚洲精品色激情综合| 老司机在亚洲福利影院| 成年免费大片在线观看| 国产黄a三级三级三级人| 熟女少妇亚洲综合色aaa.| 久久久国产欧美日韩av| 欧美另类亚洲清纯唯美| 国产91精品成人一区二区三区| 久久精品国产99精品国产亚洲性色|