當(dāng)前位置:湖南陽(yáng)光電子技術(shù)學(xué)校文章資訊技術(shù)園地電工焊工技術(shù)
招生辦公室電話(huà):13308461099 13807313137 全國(guó)免費(fèi)電話(huà):13807313137

電力系統(tǒng)短期負(fù)荷預(yù)測(cè)的模糊神經(jīng)網(wǎng)絡(luò)方法

減小字體 增大字體 作者:佚名  來(lái)源:本站整理  發(fā)布時(shí)間:2010-07-17 22:31:55
 

電力系統(tǒng)短期負(fù)荷預(yù)測(cè)的模糊神經(jīng)網(wǎng)絡(luò)方法

姜 勇 盧 毅  東南大學(xué)(南京210096)

    摘 要 針對(duì)電力系統(tǒng)短期負(fù)荷預(yù)測(cè)問(wèn)題,考慮到氣象因素對(duì)負(fù)荷的影響,提出了一種模糊神經(jīng)網(wǎng)絡(luò)的短期負(fù)荷預(yù)測(cè)方法,首先根據(jù)評(píng)價(jià)函數(shù)選取相似日學(xué)習(xí)樣本,然后利用隸屬函數(shù)對(duì)影響負(fù)荷的特征因素向量的分量進(jìn)行模糊處理,采用反向傳播算法,對(duì)24點(diǎn)每點(diǎn)建立一個(gè)預(yù)測(cè)模型,提高了學(xué)習(xí)效能。本方法適合在短期負(fù)荷預(yù)測(cè)中使用,具有較好的預(yù)測(cè)精度。
  關(guān)鍵詞 短期負(fù)荷預(yù)測(cè) 隸屬函數(shù) 人工神經(jīng)網(wǎng)絡(luò)


0 引 言
 
  短期負(fù)荷預(yù)報(bào)是電力系統(tǒng)管理現(xiàn)代化的重要內(nèi)容之一,是對(duì)發(fā)電、輸電和電能分配等合理介紹的必要前提,對(duì)提高電力系統(tǒng)的經(jīng)濟(jì)效益和社會(huì)效益,保障電力系統(tǒng)的安全經(jīng)濟(jì)運(yùn)行與國(guó)民經(jīng)濟(jì)的發(fā)展具有非常重要的影響。因此,尋求有效的負(fù)荷預(yù)報(bào)方法,提高預(yù)測(cè)結(jié)果的準(zhǔn)確度具有重要意義。
  傳統(tǒng)的預(yù)測(cè)模型是用顯式的數(shù)學(xué)表達(dá)式加以描述,這就決定了它的局限性。事實(shí)上,電力負(fù)荷變化受天氣情況和人們的社會(huì)活動(dòng)等因素的強(qiáng)烈影響,存在大量非線(xiàn)性關(guān)系,其發(fā)展規(guī)律很難用一個(gè)顯式的數(shù)學(xué)公式予以表示。因此,將具備模糊數(shù)據(jù)處理能力的模糊理論與擅長(zhǎng)擬合非線(xiàn)性映射的神經(jīng)網(wǎng)絡(luò)方法結(jié)合起來(lái),是一種比較有效的預(yù)測(cè)技術(shù)。

1 模糊理論與人工神經(jīng)網(wǎng)絡(luò)方法簡(jiǎn)介〔1,2〕 
 
  模糊神經(jīng)網(wǎng)絡(luò)是通過(guò)將模糊理論和人工神經(jīng)網(wǎng)絡(luò)理論有機(jī)結(jié)合來(lái)實(shí)現(xiàn)的,常見(jiàn)的結(jié)合方式有:a.仍采用普通神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu),但將普通非線(xiàn)性神經(jīng)元用模糊運(yùn)算神經(jīng)元代替;b.采用普通神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)和神經(jīng)元作為信息處理工具,而網(wǎng)絡(luò)的輸入量、輸出量等則采用輸入、輸出信息的模糊隸屬度。本文采用的模糊神經(jīng)網(wǎng)絡(luò)為后一種類(lèi)型,即將輸入量經(jīng)過(guò)隸屬度函數(shù)轉(zhuǎn)化為模糊量后,再交給神經(jīng)網(wǎng)絡(luò)進(jìn)行處理,以提高負(fù)荷預(yù)報(bào)的精度。
1.1 模糊集合論的概念
  客觀(guān)事物的差異在中介過(guò)渡時(shí)所呈現(xiàn)的亦此亦彼的現(xiàn)象稱(chēng)為模糊性,它體現(xiàn)了事物變化的連續(xù)過(guò)程。模糊集合論使用隸屬度來(lái)描述中介過(guò)渡,是以精確的數(shù)學(xué)語(yǔ)言對(duì)模糊性的一種表述。
  設(shè)論域u={x},u到閉區(qū)間〔0,1〕的任一映射uA(x)∶u→〔0,1〕,x→uA(x)確定了u的一個(gè)模糊子集,簡(jiǎn)稱(chēng)模糊集,記作A,該映射稱(chēng)為A的隸屬函數(shù)。uA(x)的大小反映了x對(duì)模糊集A的隸屬程度,簡(jiǎn)稱(chēng)為隸屬度。實(shí)數(shù)集合上常用的隸屬函數(shù)為F分布,主要有矩形分布,梯形分布,拋物形分布等,在實(shí)際應(yīng)用中可根據(jù)對(duì)象特點(diǎn)加以選擇。
1.2 人工神經(jīng)網(wǎng)絡(luò)和反向傳播算法的原理
  神經(jīng)網(wǎng)絡(luò)是由處理單元組成的一種并行、分布式信息處理結(jié)構(gòu),處理單元之間由單向信道相互連接。人工神經(jīng)元是神經(jīng)網(wǎng)絡(luò)的基本計(jì)算單元,模擬了人腦中神經(jīng)元的基本特征,一般是多輸入/單輸出的非線(xiàn)性單元,可以有一定的內(nèi)部狀態(tài)和閾值。
  反向傳播(Error Back Propagation-BP)算法是多層感知器的一種有效學(xué)習(xí)算法,它的模型為前向多層網(wǎng)絡(luò),如圖1所示。
  網(wǎng)絡(luò)不僅有輸入層節(jié)點(diǎn)、輸出層節(jié)點(diǎn),而且有隱含層節(jié)點(diǎn),經(jīng)過(guò)作用函數(shù)后,再把隱節(jié)點(diǎn)的輸出信息傳播到輸出節(jié)點(diǎn), 后給出結(jié)果。節(jié)點(diǎn)的作用函數(shù)通常選取s型函數(shù),如http://www.hnygpx.net/article/UploadPic/2009-5/200952114183721.jpg" width=103 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>這個(gè)算法的學(xué)習(xí)過(guò)程由正向傳播和反向傳播組成。在正向傳播過(guò)程中,輸入信號(hào)從輸入層經(jīng)隱含層逐層處理,并傳向輸出層,每一層神經(jīng)元的狀態(tài)只影響下一層神經(jīng)元的狀態(tài)。如果在輸出層不能得到期望的輸出,則轉(zhuǎn)向反向傳播,將誤差信號(hào)沿原來(lái)的連接通道返回,通過(guò)修改各層神經(jīng)元的權(quán)值,使得誤差信號(hào) 小。
    BP算法可描述為:
  a.權(quán)值和閾值初始化:隨機(jī)地給全部權(quán)值和神經(jīng)元的閾值賦以較小的初始值;
    b.給定輸入xk和目標(biāo)輸出yk;

http://www.hnygpx.net/article/UploadPic/2009-5/200952114183144.jpg" width=361 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>

http://www.hnygpx.net/article/UploadPic/2009-5/200952114183902.jpg" width=285 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>

其中http://www.hnygpx.net/article/UploadPic/2009-5/200952114183694.jpg" width=29 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>為第l層第j個(gè)神經(jīng)元到第l+1層的第i個(gè)神經(jīng)元的權(quán)系數(shù),η為增益項(xiàng),δ(l)ik為第l層i節(jié)點(diǎn)的k模式的誤差項(xiàng),且有
    若i為輸出節(jié)點(diǎn),則

http://www.hnygpx.net/article/UploadPic/2009-5/200952114183415.jpg" width=354 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>

  e.達(dá)到誤差精度要求則輸出結(jié)果,否則回到3。

2 預(yù)測(cè)方法〔3〕

2.1 BP網(wǎng)絡(luò)學(xué)習(xí)樣本的選擇
  根據(jù)差異評(píng)價(jià)函數(shù)選擇一批比較相似的負(fù)荷日,選擇方法如下:
  a.首先將影響因素向量的各分量數(shù)值化,這些因素括日類(lèi)型,工作日取1,雙休日取2;光照,晴天取1,少云取2,多云取3,陰天取4,下雨取5;雨量,無(wú)雨取0,小雨取1,中雨取2,大雨取3,暴雨取4; 高溫度、 低溫度、平均溫度等可取實(shí)際值。
    b.建立評(píng)價(jià)函數(shù)

http://www.hnygpx.net/article/UploadPic/2009-5/200952114183592.jpg" width=320 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>

  其中是α-β的范數(shù),α為已知日的影響因素向量,β為預(yù)測(cè)日的影響因素向量(預(yù)測(cè)值),δ為設(shè)定的偏離值,要根據(jù)不同地區(qū)的實(shí)際情況由經(jīng)驗(yàn)決定。凡是滿(mǎn)足該評(píng)價(jià)函數(shù)的已知日均可加入學(xué)習(xí)樣本。
2.2 預(yù)測(cè)方法
  為了提高BP網(wǎng)絡(luò)的學(xué)習(xí)速度和預(yù)測(cè)精度,對(duì)網(wǎng)絡(luò)輸出的實(shí)際負(fù)荷值進(jìn)行規(guī)一化處理,同時(shí)使用隸屬度函數(shù)對(duì)輸入網(wǎng)絡(luò)的影響因素向量進(jìn)行模糊化處理,使得輸入向量的不同分量如工作日、光照、雨量、溫度等都轉(zhuǎn)換為模糊量。
  日類(lèi)型的隸屬度函數(shù)采用半矩形分布,形式如下。
    對(duì)工作日的隸屬度函數(shù)為:

http://www.hnygpx.net/article/UploadPic/2009-5/200952114183888.jpg" width=358 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>
    溫度的隸屬度函數(shù)分別如下。
    對(duì)低溫的隸屬度函數(shù)采用偏小型梯形分布:

http://www.hnygpx.net/article/UploadPic/2009-5/200952114183955.jpg" width=350 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>
    對(duì)中溫的隸屬度函數(shù)采用中間型梯形分布:

http://www.hnygpx.net/article/UploadPic/2009-5/200952114184841.jpg" width=342 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>
  對(duì)高溫的隸屬度函數(shù)采用偏大型梯形分布:

http://www.hnygpx.net/article/UploadPic/2009-5/200952114184891.jpg" width=342 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>
  將 高溫度th代入以上3個(gè)公式,可分別求出對(duì)于低溫、中溫、高溫3個(gè)狀態(tài)的隸屬度。對(duì)于 低溫度t1、平均溫度tα,可采用同樣方法求出它們的3個(gè)狀態(tài)隸屬度。
  影響因素向量的其它分量如光照、雨量等的隸屬度函數(shù)與溫度相似,根據(jù)當(dāng)?shù)貙?shí)際情況選擇建立相應(yīng)的分布函數(shù),從而求出光照的5個(gè)狀態(tài)隸屬度以及雨量的5個(gè)狀態(tài)隸屬度。
  另外,為了降低求解規(guī)模,對(duì)1d的24點(diǎn)負(fù)荷每點(diǎn)建立1個(gè)預(yù)測(cè)模型,第i點(diǎn)的輸入輸出函數(shù)為:

http://www.hnygpx.net/article/UploadPic/2009-5/200952114184838.jpg" width=356 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>

  其中Lji是指第j個(gè)學(xué)習(xí)樣本第i點(diǎn)負(fù)荷的實(shí)際值,Ljimin是學(xué)習(xí)樣本中第i點(diǎn)負(fù)荷的 小值,Ljimax是學(xué)習(xí)樣本中第i點(diǎn)負(fù)荷的 大值,經(jīng)過(guò)規(guī)一化處理后,http://www.hnygpx.net/article/UploadPic/2009-5/200952114184773.jpg" width=87 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>
  fi是指第i點(diǎn)的輸入輸出函數(shù),α′j1,…α′jl是第j個(gè)樣本影響因素的隸屬度向量,括2個(gè)代表日類(lèi)型的隸屬度,3個(gè)代表 高溫度的隸屬度,3個(gè)代表 低溫度的隸屬度,3個(gè)代表平均溫度的隸屬度,5個(gè)代表光照的隸屬度,5個(gè)代表雨量的隸屬度。當(dāng)學(xué)習(xí)完成后,將預(yù)測(cè)日當(dāng)天影響因素的隸屬度向量http://www.hnygpx.net/article/UploadPic/2009-5/200952114184512.jpg" width=52 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>代入(10)式,則可得第i點(diǎn)的負(fù)荷預(yù)測(cè)值。
2.3 一些注意事項(xiàng)
    a.偽數(shù)據(jù)的處理
  因?yàn)樯窠?jīng)網(wǎng)絡(luò)所用的負(fù)荷數(shù)據(jù)來(lái)自電力部門(mén)的SCADA系統(tǒng),由于各種原因會(huì)造成一定數(shù)量的異常數(shù)據(jù)。考慮到負(fù)荷前后小時(shí)的自然變化,如果出現(xiàn)超常規(guī)值,必須將其剔除,代之以正常比例范圍內(nèi)的估計(jì)值。
    b.待選的相似日范圍
  因?yàn)殡S著時(shí)間的推移,系統(tǒng)負(fù)荷結(jié)構(gòu)會(huì)發(fā)生緩慢的變化,當(dāng)已知日和預(yù)測(cè)日相隔較遠(yuǎn)時(shí),即使它們的天氣情況等因素很相似,預(yù)測(cè)精度也不會(huì)高,因而取前3個(gè)星期的已知日作為待選范圍,同時(shí)還可以縮短程序選取樣本所花費(fèi)的時(shí)間。

3 計(jì)算實(shí)例分析與結(jié)論

3.1 實(shí)例計(jì)算
  對(duì)南京市某日24點(diǎn)的負(fù)荷進(jìn)行預(yù)測(cè),所得結(jié)果如表1。

http://www.hnygpx.net/article/UploadPic/2009-5/200952114184787.jpg" width=345 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>
http://www.hnygpx.net/article/UploadPic/2009-5/200952114184441.jpg" width=345 border=0 onload="return imgzoom(this,600);" onclick="javascript:window.open(this.src);" style="cursor:pointer;"/>

  預(yù)測(cè)的平均絕對(duì)百分誤差為2.56%, 大預(yù)測(cè)誤差為5.1%, 小誤差為0.1%,誤差超過(guò)5%的預(yù)測(cè)點(diǎn)有1個(gè),小于3%的點(diǎn)有15個(gè),預(yù)測(cè)效果良好。
3.2 結(jié) 論
  準(zhǔn)確進(jìn)行短期負(fù)荷預(yù)測(cè)是電力行業(yè)所企盼的,本文提出一種模糊神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)方法,利用模糊技術(shù)和神經(jīng)網(wǎng)絡(luò)各自的特長(zhǎng),充分發(fā)揮了ANN處理非線(xiàn)性問(wèn)題的能力,具有訓(xùn)練速度快,學(xué)習(xí)精度高,數(shù)值穩(wěn)定等優(yōu)點(diǎn),算例也表明這是一種行之有效的短期日負(fù)荷預(yù)測(cè)方法。

參考文獻(xiàn)
1 楊綸標(biāo),高英儀.模糊數(shù)學(xué)原理及應(yīng)用.華南理工大學(xué)出版社.2001,3
2 胡守仁,余少波,戴葵.神經(jīng)網(wǎng)絡(luò)導(dǎo)論.長(zhǎng)沙.國(guó)防科技大學(xué)出版社.1993,10
3 牛東曉,曹樹(shù)華,趙磊,張文文.電力負(fù)荷預(yù)測(cè)技術(shù)及其應(yīng)用.北京.中國(guó)電力出版社.1998,10

職稱(chēng)論文代寫(xiě)代發(fā)表(國(guó)家級(jí)發(fā)表3個(gè)月,省級(jí)期刊2個(gè)月)...

【javascript:history.back();">返回上頁(yè)】【javascript:window.print();">打印此頁(yè)】【查看參與評(píng)論】

Tags:

作者:佚名
[]
電工培訓(xùn)學(xué)校 電動(dòng)車(chē)維修學(xué)校 摩托車(chē)維修學(xué)校 摩托車(chē)維修培訓(xùn) 手機(jī)維修培訓(xùn) 家電維修培訓(xùn) 電腦維修培訓(xùn) 電動(dòng)工具維修培訓(xùn) 液晶電視維修培訓(xùn) 安防監(jiān)控培訓(xùn) 空調(diào)維修培訓(xùn) 網(wǎng)絡(luò)營(yíng)銷(xiāo)培訓(xùn) 網(wǎng)站設(shè)計(jì)培訓(xùn) 淘寶網(wǎng)店培訓(xùn) 電器維修培訓(xùn) 家電維修學(xué)校 電工培訓(xùn) 焊工培訓(xùn) 電工學(xué)校 電工培訓(xùn)學(xué)校 電動(dòng)車(chē)維修學(xué)校 摩托車(chē)維修學(xué)校 摩托車(chē)維修培訓(xùn) 手機(jī)維修培訓(xùn) 家電維修培訓(xùn) 電腦維修培訓(xùn) 電動(dòng)工具維修培訓(xùn) 液晶電視維修培訓(xùn) 安防監(jiān)控培訓(xùn) 空調(diào)維修培訓(xùn) 網(wǎng)絡(luò)營(yíng)銷(xiāo)培訓(xùn) 網(wǎng)站設(shè)計(jì)培訓(xùn) 淘寶網(wǎng)店培訓(xùn) 電器維修培訓(xùn) 家電維修學(xué)校 電工培訓(xùn) 焊工培訓(xùn) 電工學(xué)校 電工培訓(xùn)學(xué)校 電動(dòng)車(chē)維修學(xué)校 摩托車(chē)維修學(xué)校 摩托車(chē)維修培訓(xùn) 手機(jī)維修培訓(xùn) 家電維修培訓(xùn) 電腦維修培訓(xùn) 電動(dòng)工具維修培訓(xùn) 液晶電視維修培訓(xùn) 安防監(jiān)控培訓(xùn) 空調(diào)維修培訓(xùn) 網(wǎng)絡(luò)營(yíng)銷(xiāo)培訓(xùn) 網(wǎng)站設(shè)計(jì)培訓(xùn) 淘寶網(wǎng)店培訓(xùn) 電器維修培訓(xùn) 家電維修學(xué)校 電工培訓(xùn) 焊工培訓(xùn) 電工學(xué)校 電工培訓(xùn)學(xué)校 電動(dòng)車(chē)維修學(xué)校 摩托車(chē)維修學(xué)校 摩托車(chē)維修培訓(xùn) 手機(jī)維修培訓(xùn) 家電維修培訓(xùn) 電腦維修培訓(xùn) 電動(dòng)工具維修培訓(xùn) 液晶電視維修培訓(xùn) 安防監(jiān)控培訓(xùn) 空調(diào)維修培訓(xùn) 網(wǎng)絡(luò)營(yíng)銷(xiāo)培訓(xùn) 網(wǎng)站設(shè)計(jì)培訓(xùn) 淘寶網(wǎng)店培訓(xùn) 電器維修培訓(xùn) 家電維修學(xué)校 電工培訓(xùn) 焊工培訓(xùn) 電工學(xué)校
中山市,固原市,銀川市,玉樹(shù),海東,隴南市,酒泉市,張掖市,天水市,金昌市,蘭州市,榆林市,延安市,渭南市,銅川市,阿里,山南,拉薩市,怒江,文山州,楚雄州,普洱市,昭通市,玉溪市,昆明市,畢節(jié),銅仁,遵義市,貴陽(yáng)市,甘孜州,資陽(yáng)市,達(dá)州市,宜賓市,南充市,遂寧市,綿陽(yáng)市,瀘州市,自貢市,三亞市,崇左市,河池市,玉林市,欽州市,梧州市,柳州市,梅州市,肇慶市,湛江市,佛山市,珠海市,韶關(guān)市,湘西州,懷化市,郴州市,張家界市,邵陽(yáng)市,株洲市,仙桃市,隨州市,荊州市,荊門(mén)市,襄樊市,黃石市,駐馬店市,信陽(yáng)市,南陽(yáng)市,漯河市,中衛(wèi)市,石嘴山市,海西,海南藏州,黃南州,海北,甘南,慶陽(yáng)市,平?jīng)鍪?武威市,白銀市,嘉峪關(guān)市,安康市,漢中市,咸陽(yáng)市,寶雞市,林芝,日喀則,昌都,迪慶,德宏,大理,西雙版納,紅河州,臨滄市,麗江市,保山市,曲靖市,黔東州,黔西州,安順市,六盤(pán)水市,涼山州,阿壩州,雅安市,廣安市,眉山市,內(nèi)江市,廣元市,德陽(yáng)市,攀枝花市,成都市,海口市,來(lái)賓市,百色市,貴港市,北海市,桂林市,南寧市,云浮市,揭陽(yáng)市,潮州市,清遠(yuǎn)市,陽(yáng)江市,汕尾市,惠州市,茂名市,江門(mén)市,汕頭市,深圳市,廣州市,婁底市,永州市,益陽(yáng)市,岳陽(yáng)市,湘潭市,長(zhǎng)沙市,恩施州,黃岡市,孝感市,鄂州市,十堰市,武漢市,周口市,商丘市,三門(mén)峽市,許昌市,焦作市,安陽(yáng)市,鶴壁市,平頂山市,開(kāi)封市,鄭州市,聊城市,濱州市,德州市,萊蕪市,日照市,泰安市,煙臺(tái)市,濰坊市,東營(yíng)市,淄博市,上饒市,濟(jì)南市,撫州市,宜春市,贛州市,新余市,九江市,景德鎮(zhèn)市,寧德市,南平市,泉州市,莆田市,廈門(mén)市,宣城市,亳州市,六安市,宿州市,黃山市,滁州市,安慶市,淮北市,馬鞍山市,蚌埠市,蕪湖市,合肥市,麗水市,舟山市,衢州市,金華市,湖州市,嘉興市,寧波市,宿遷市,鎮(zhèn)江市,鹽城市,連云港市,蘇州市,徐州市,南京市,綏化市,牡丹江市,佳木斯市,大慶市,鶴崗市,哈爾濱市,白城市,白山市,遼源市,吉林市,葫蘆島市,鐵嶺市,盤(pán)錦市,阜新市,錦州市,本溪市,鞍山市,沈陽(yáng)市,錫林郭勒盟,通遼市,烏海市,呂梁市,忻州市,晉中市,晉城市,陽(yáng)泉市,太原市,廊坊市,承德市,保定市,邯鄲市,唐山市,寧夏,甘肅省,西藏,貴州省,重慶市,廣西,湖南省,河南省,江西省,安徽省,江蘇省,黑龍江省,遼寧省,山西省,天津市,四平市,內(nèi)蒙古,吳忠市,果洛,西寧市,定西市,商洛市,西安市,那曲,黔南州,巴中市,樂(lè)山市,賀州市,防城港市,東莞市,河源市,常德市,衡陽(yáng)市,咸寧市,宜昌市,濮陽(yáng)市,新鄉(xiāng)市,洛陽(yáng)市,菏澤市,臨沂市,威海市,濟(jì)寧市,棗莊市,青島市,吉安市,鷹潭市,萍鄉(xiāng)市,南昌市,龍巖市,漳州市,三明市,福州市,池州市,巢湖市,阜陽(yáng)市,銅陵市,淮南市,臺(tái)州市,紹興市,溫州市,杭州市,泰州市,揚(yáng)州市,淮安市,南通市,常州市,無(wú)錫市,大興安嶺,黑河市,七臺(tái)河市,伊春市,雙鴨山市,雞西市,齊齊哈爾市,延邊,松原市,通化市,長(zhǎng)春市,朝陽(yáng)市,遼陽(yáng)市,營(yíng)口市,丹東市,撫順市,大連市,阿拉善盟,興安盟,烏蘭察布市,巴彥淖爾市,呼倫貝爾市,鄂爾多斯市,赤峰市,頭市,呼和浩特市,臨汾市,運(yùn)城市,朔州市,長(zhǎng)治市,大同市,衡水市,滄州市,張家口市,邢臺(tái)市,秦皇島市,石家莊市,青海省,陜西省,云南省,四川省,海南省,廣東省,湖北省,山東省,福建省,浙江省,上海市,吉林省,河北省,北京市 主站蜘蛛池模板: 微拍秒拍99福利精品小视频 | 在线观看日本污污ww网站 | 欧美国产成人在线 | 美女扒开屁股让男人桶 | 中文字幕久久熟女人妻AV免费 | 久久久久琪琪精品色 | 欧美精品高潮呻吟AV久久 | 国语自产一区视频 | 大乳牛奶女在线观看 | 国产免费啪嗒啪嗒视频看看 | 3d在线看小舞被躁视频 | 18美女腿打开无遮软件 | 免费国产成人高清在线观看视频 | 99午夜视频 | 国产精品爽爽久久久久久竹菊 | 久久99re66热这里只有精品 | 一一本之道高清视频在线观看中文字幕 | 天天看片视频免费观看 | 中文无码第3页不卡av | 中文字幕爆乳JULIA女教师 | 无码任你躁久久久久久老妇双奶 | 好大的太粗好深BL | 五花大绑esebdsm国产 | 在线免费观看日本 | 18未满不能进的福利社 | 国产69精品久久久久乱码免费 | 西施打开双腿下面好紧 | 迅雷哥在线观看高清 | 揉抓捏打抽插射免费视频 | 纯肉高H种马艳遇风流多 | 一本大道无码AV天堂欧美 | 久久综合九色综合国产 | SM脚奴调教丨踩踏贱奴 | 国产综合欧美区在线 | 久久这里只有精品国产99 | 男人天堂黄色 | 精品久久久久中文字幕 | 99久久无码一区人妻A片竹菊 | 无码日韩人妻精品久久蜜桃免费 | 一道精品视频一区二区 | 美女被免费喷白浆视频 |