999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

考慮電動汽車充電站選址定容的配電網規劃

2013-07-05 15:14:51陳婷衛志農吳霜孫國強韋延方許曉慧
電力系統及其自動化學報 2013年3期
關鍵詞:配電網規劃優化

陳婷,衛志農,吳霜,孫國強,韋延方,許曉慧

(1.河海大學能源與電氣學院,南京 210098;2.中國電力科學研究院,南京 210003)

考慮電動汽車充電站選址定容的配電網規劃

陳婷1,衛志農1,吳霜1,孫國強1,韋延方1,許曉慧2

(1.河海大學能源與電氣學院,南京 210098;2.中國電力科學研究院,南京 210003)

針對電動汽車充電站選址定容的配電網規劃問題,建立了以配電網投資及運行成本、電動汽車充電站投資及運行成本最小為優化目標的模型。采用帶精英策略的改進樹形結構編碼單親遺傳算法ITSE-PGA(improved tree structure encoding partheno-genetic algorithm),在配電網規劃的同時進行了電動汽車充電站位置和容量的優化。不同于用于傳統的配電網規劃的樹形結構編碼單親遺傳算法,該方法在含電動汽車充電站的優化過程中,配電網絡不一定連通,所以要對不可行解做特殊處理。通過算例仿真,驗證了該方法的可行性和有效性。

電動汽車充電站;配電網規劃;單親遺傳算法;樹形結構編碼

近年來,隨著全球環境的日益惡化以及石油資源的日趨緊張,電動汽車作為新型能源交通工具越來越受到關注,其發展前景廣闊[1]。充電站、充電樁及電池等充電設施,是發展電動汽車的重要配套基礎設施[2,3]。電動汽車充電站的接入將會引起配電網負荷的增長,以及配電網投資與運行成本的提高。為以較低成本滿足電動汽車的充電需求,亟需研究考慮充電站位置和容量優化的配電網規劃問題[4]。全面地規劃優化配電網架結構,能夠有效地縮小容量要求、降低網絡損耗、減少施工投入,提高電力公司人力物力資源的利用率,有效降低建設投資和維護費用,為國家和電力公司帶來可觀的經濟效益。

目前已有文獻對電動汽車充電站的規劃問題進行了研究[5~11]。文獻[5]提出用2步篩選法來確定充電站的候選站址,構造了電動汽車充電站最優規劃的數學模型,并采用改進的原對偶內點法來求解。文獻[6]分析了影響電動汽車充電站規劃的諸多因素,并就充電站的布局規劃問題,提出了應滿足充電半徑要求,與電動汽車交通密度、充電需求分布、城市總體規劃相配合等的原則性建議。文獻[7]應用區域交通流量守恒定理,將小區內的電動汽車總量視為充電站的負荷點,在此基礎上做出電動汽車的總量和分布預測,提出了充電站的年最小費用模型,并運用遺傳算法對實例進行求解。文獻[8]運用動態交通網絡思想建立了基于硬時間窗約束下的充電站布局及最佳規模確定的多目標優化模型,該模型將充電者充電成本和充電站投資者成本最小化作為優化目標。綜上所述,對于電動汽車充電站的選址定容研究,主要考慮了影響充電站本身經濟運行的諸多因素,而充電站與配電網結合進行規劃的文獻目前較少。文獻[4]考慮了電動汽車用戶的方便性和充電站接入電網的成本,文中提出的輻射狀約束可以保證規劃方案的輻射狀結構。文獻[9]針對集中型充電站接入的電網規劃問題,綜合考慮電力網絡和交通網絡因素,建立了集中型充電站的定址分容模型。電動汽車充電站作為電網的一部分,它的接入必然對電網的規劃造成影響,比如網絡投資和網損的增加,線路和變壓器過載,帶來新的負荷高峰,惡化電能質量等。如何使充電站成本和電網建設及運行成本最小是本文研究的問題。

在配電網規劃的求解算法中,簡單遺傳算法、禁忌搜索、蟻群算法、粒子群算法等隨機優化算法在生成個體時會產生大量的不可行解,導致計算速度慢,因此需要構造合理的編碼及進化策略,以提高算法效率[12~18]。樹形結構編碼單親遺傳算法應用于一般的配電網規劃中產生的解均是可行解[17],本文采用改進的樹形結構編碼單親遺傳算法ITSE-PGA(improved tree structure encoding partheno-genetic algorithm),在配電網規劃的同時進行了電動汽車充電站位置和容量的優化,在優化過程中產生的解不全是可行解,但只需在解碼過程中做出簡單的判斷,計算速度快,具有良好的尋優效果。

1 數學模型

1.1 問題描述

本文的研究對象為一個包含供電首端、充電站備選站址節點和傳統負荷節點的配電網絡。另外為了規劃充電站備選站址的容量,本文引入了單位容量節點,使其與所有的充電站備選站址虛擬相連。

1.2 目標函數

以配電網中線路的投資費用f1、運行費用(主要指網損)f2以及電動汽車充電站的投資及運行費用f3總和最小為目標建立模型,其模型為

式中:f1為線路投資費用,萬元;f2為規劃年限內的網絡損耗(主要指線損)費用,萬元;f3為在規劃年限內充電站付出的總成本,萬元;D為新建線路集合;li為第i條線路的長度,km;αi為第i條線路單位長度的建設費用,萬元/km;T為規劃年限,a;τ為年最大負荷利用小時數,h;γ為單位電量的費用系數,元/(kW·h);Ii為流過第i條線路的電流,A;ri為第i條線路的電阻,Ω;n為充電站備選站址的個數;Ki,f為充電站備選站址i新建的固定成本,萬元;xi為充電站i是否建設的二元決策變量,取值為1表示該待選站址將建設充電站,取值為0表示不在該待選站址建設充電站;ki,v為折算到一年,單位容量對應的平均可變投資,萬元/MW;ki,o為折算到一年,單位容量對應的平均運行成本,萬元/ MW;Qi為充電站備選站址i的容量,MW;Pi為頂點i的年最大負荷平均值,kW;UN為配電網的額定電壓,kV;cos φ為功率因數;ρi為電阻率,Ω·mm2/ km;αi為截面積,mm2。

1.3 約束條件

1)輻射狀結構和連通性約束

網架結構規劃首先要解決的就是網絡的連通性和輻射性約束問題,即避免通常所說的“孤島”與“環島”情形。

2)功率平衡約束

本文采用直流潮流作為網絡功率平衡的約束條件[9],負荷節點有功功率平衡約束為

式中:Pij為線路ij的有功傳輸功率;Pj,load為節點j處除充電站外的負荷;Pj,station為充電站的負荷;xj,station為節點j處是否建設充電站的二元決策變量;A為負荷節點集合;N為網絡總結點數。

3)不等式約束

1.4 約束條件的處理

本文輻射性約束自動滿足。由于充電站單位容量節點的加入,在隨機優化過程中,可能會將充電站單位容量節點當作交叉節點進行規劃,從而破壞網絡的連通性。本文在編寫程序時將此種情況予以排除,保證了網絡的連通性。線路功率約束在選線型時加以約束,線路中流過的電流不超過線路能承受的最大電流。

2 樹形結構編碼單親遺傳算法及其改進

2.1 染色體表達方式

本文中單親遺傳算法采用樹形結構編碼,將配電網中各負荷點及電源點之間的連接關系和充電站單位容量節點與充電站備選站址節點之間的連接關系用多個樹形結構(森林)表示。用一個節點類的對象數組Node代表染色體,采用節點類對象的屬性變量保存各節點及其父節點、子節點的編號,據此保存圖的樹形結構信息。并通過改變各節點的父節點及對應子節點,來達到改變圖的拓撲結構的目的。線路的修建與否、充電站單位容量節點與充電站備選站址節點之間連接與否由各節點間是否存在父子關系來決定。若兩節點間存在父子關系,則修建對應的線路或連接對應的充電站單位容量節點與充電站節點,否則不修建或不連接。

2.2 單親遺傳操作及其優點

單親遺傳算法PGA(partheno genetic algorithm)與傳統遺傳算法TGA(traditional genetic algorithm)類似,也是一種采用隨機搜索方式的種群算法。其基本原理[17]是:通過遺傳算子作用于父代種群,從而產生適應性更強的后代種群;通過反復上述過程,進而達到種群進化的目的。由于本文采用樹形結構編碼的單親遺傳算法,考慮到改變樹的結構所產生的效果主要是改變了某個或某些節點的父節點,故本文的單親遺傳算子采用移位算子和重分配算子[17]。這兩種算子的優點在于它們在對個體進行隨機優化時不改變網絡的輻射性和連通性結構。

2.3 適應度計算

適應度函數由數據模型中的目標函數所決定的,數據優化模型中的目標函數是以經濟效益為參考的。在種群產生之后,通過譯碼操作,將種群中每個個體的染色體轉化為可行解對應的數據變量,然后代入求解模型的目標函數中,獲得每組解所對應的經濟效益。經濟效益即為每個個體對應的適應度,經濟效益差的個體被淘汰,經濟效益好的個體優先生存,通過交叉、變異,最終產生最優解。

算法搜索到的最優解要使目標函數值F盡量小,而適應度作為評價優劣的標準總是尋求最大解,且適應度大于零。因此本文取適應度

式中,M為足夠大的常數,以保證適應度恒為正值。

2.4 精英保留策略

精英保留策略是為了解決由于隨機因素導致優化過程中優秀個體丟失的問題。算法開始時隨機產生一個初始群體A0,在此基礎上采用移位、重分配等遺傳操作產生一個新群體B0,A0和B0的群體規模均為N。將At和Bt并入到Ct中(初始時t= 0),針對Ct計算所有個體的適應度并從大到小排序,從中選取前N個個體進入At+1,直到t達到最大迭代次數。

2.5 電動汽車充電站位置和容量的優化

文獻[17]將樹形結構編碼單親遺傳算法應用于單一的配電網規劃問題中獲得了很好的效果。而本文在配電網規劃的同時進行電動汽車充電站位置和容量的優化,以期達到兩者經濟性最佳。因此本文引入充電站單位容量節點(1 MW)虛擬地與各個充電站備選站址節點相連,將充電站單位容量節點當作負荷節點,一起參與網絡拓撲的規劃。在個體解碼時,首先檢查是否將充電站單位容量節點當作中間交叉節點規劃到網絡拓撲中。當單位充電容量節點5同時與充電站備選站址節點3和節點4相連(如圖1所示),此為不可行解。從整體看,網絡呈輻射連通狀,但去除單位充電容量節點,網絡不連通,因此排除這一個體。由充電站備選站址節點連接的單位容量節點數,可以知道充電站備選站址節點容量。如圖2所示,充電站備選站址3連接了一個單位容量節點5,故充電站備選站址節點容量為1 MW,充電站備選站址4連接了兩個單位容量節點6和節點7,故充電站備選站址節點容量為2 MW。然后去除虛擬連接的充電站單位容量節點,計算目標函數以及適應度。將容量規劃的連續問題離散化后,采用改進的樹形結構編碼單親遺傳算法,快速解決了電動汽車充電站選址定容和配電網規劃聯合優化問題。

圖1 不可行解示例Fig.1One example of infeasible solutions

圖2 可行解確定容量Fig.2Determining capacity with a feasible solution

2.6 初始方案的產生

本文結合kruskal最小生成樹算法產生初始群體,以獲得可行(輻射及連通)且相比隨機產生的網絡較優的初始方案,步驟如下。

步驟1將充電站單位容量節點與充電站備選站址節點之間虛擬支路長度設為一個相同的大值,對所有支路(包括充電站單位容量節點與充電站節點之間虛擬支路)按長度從小到大進行排序,存入一個二維數組Q[n][2]中,Q[i][0]存支路首節點編號,Q[i][1]存支路末節點編號。

步驟2標記各節點標識signal為各節點編號i。

步驟3從Q中取出一條支路加入T[m][2]中,看這條支路的首節點與末節點標識是否相同,如果相同則刪除該支路;如果不同,則該邊加入T,且遍歷已形成的樹,將所有與末節點標識相同的節點標識修改為首節點標識。

步驟4檢查是否所有的節點都已連接,如果是則遍歷整棵樹,在Node[]數組中相應對象的屬性變量中保存此二節點的父子關系,結束該初始方案,否則轉步驟3。

3 算法流程

樹形結構編碼單親遺傳算法的具體步驟如下:

步驟1產生初始化種群A0,迭代次數t=0;

步驟2解碼,適應度值計算,進行個體評價;

步驟3判斷t是否大于最大迭代次數,如果是,輸出最優解;否則,進入下一步;

步驟4移位、重分配后解碼看是否為不可行解,如果是,則去除;否則由充電站備選站址節點連接了幾個單位容量節點,則可以知道充電站備選站址節點容量。并計算目標函數以及適應度值。如此循環往復直到得到相應個數的新種群Bt;

步驟5Ct=At∪Bt,針對Ct中的每個個體按適應度值從大到小排序;

步驟6采用精英保留策略,選取前N個個體進入At+1;

步驟7t=t+1,轉步驟3。

4 算例分析

算例如圖3所示[19],區域總充電需求為6 MW,數字為各節點編號,充電站備選站址為節點4~6,節點0~9間的虛線為可以架線的走廊,節點4~6與節點10~15間的虛線為虛擬連接的支路,以便確定充電站備選站址節點的容量。表1為網絡節點負荷情況,可以架設線路的長度如表2所示,充電站備選站址參數如表3所示,模型求解中采用的單親遺傳操作參數如表4所示。補充如下數據:取單位電量費用系數為0.5元/(kW·h),規劃周期為10 a,年最大負荷利用小時數為4 000 h,配電網額定電壓為10 kV,各個負荷點的功率因數cos φ為0.85。導線提供LGJ-50,LGJ-70,LGJ-95,LGJ-120,LGJ-150,LGJ-180和LGJ-240等7種型號,且只能單回路架線。

圖3 待規劃的配電網絡Fig.3Initial distribution network to be planned

表1 網絡節點負荷Tab.1Loads of network nodes

表2 線路長度Tab.2Length of the line

表3 充電站備選站址參數Tab.3Parameters of alternative sites for charging station

表4 單親遺傳算法的操作參數Tab.4Operation parameters of PGA

分2種方案進行優化規劃計算。方案1采用表3的參數計算,方案2將充電站備選站址節點4~6的容量設為固定值,分別為2、1、3MW。

在Visual C++環境下,在CPU主頻2.00 GHz、操作系統為Windows XP的計算機上上述算例的計算時間均不到1s。

圖4為方案1(即本文算法)的規劃結果,其中表5為方案1網架規劃結果,表6為方案1的充電站定址分容結果。利用改進的樹形結構編碼單親遺傳算法得出最佳個體生成代數為16,最小費用為1 338.617萬元。

圖4 方案1規劃結果Fig.4Planning results of scheme 1

圖5 方案2規劃結果Fig.5Planning results of scheme 2

表5 方案1網架規劃結果Tab.5 Results of distributionnet work planning of scheme1

表6 方案1充電站定址分容結果Tab.6 Results of site and capacity of stations of scheme 1

表7 方案2網架規劃結果Tab.7 Results of distribution net work planning of scheme2

表8 方案2充電站定址分容結果Tab.8 Results of site and capacity of stations of scheme 2

表92 種方案的成本比較Tab.9Cost comparison between two schemes

圖5為方案2的規劃結果,其中表7為方案2網架規劃結果,表8為方案2的充電站定址分容結果。最小費用為1 482.481萬元。

方案2指定各充電站備選站址的容量,而方案1利用改進的樹形結構編碼單親遺傳算法將電動汽車充電站的選址和定容考慮在配電網規劃模型中。由表9可知,方案1的規劃結果比方案2的規劃方案節省了大約10%的總成本。由規劃結果對比可知,將電動汽車充電站的選址和定容考慮到配電網優化模型中可有效減少包含充電站在內的配電網的總投資與運行成本。

本文提出的改進樹形結構編碼單親遺傳算法應用于配電網規劃和電動汽車充電站定址分容的綜合,不僅有效地解決了問題,而且能快速找到最優解,計算速度快。

5 結語

本文將改進的樹形結構編碼單親遺傳算法應用于配電網規劃和電動汽車充電站定址分容,有效地解決了電動汽車充電站接入電網的規劃問題。該算法在優化過程中配電網絡始終自然呈輻射狀,無需輻射性檢驗,提高了求解問題的效率。然而,由于樹形結構編碼的特殊性,電動汽車充電站容量優化這一連續性問題被離散化,因此,若要得到更精確的結果,所取單位容量需更為縮小。本文中,單位容量取1 MW,所得解為相對最優解。若將單位容量取為0.1 MW甚至更小,這樣從理論上來說可提高求解的精度,然而因為大大增大了解空間,對于隨機優化的智能算法而言,要想得到最優解將需要比較多的迭代次數,也可能陷入局部最優,所以選取合適的單位容量是比較重要的。

[1]Parks K,Denholm P,Markel T.Costs and emissions associated with plug-in hybrid electric vehicle charging in the Xcel energy Colorado service territory[EB/OL].http:// www.nrel.gov/vehiclesandfuels/pdfs/41410.pdf,2007.

[2]吳春陽,黎燦兵,杜力,等(Wu Chunyang,Li Canbing,Du Li,et al).電動汽車充電設施規劃方法(A method for electric vehicle charging infrastructure planning)[J].電力系統自動化(Automation of Electric Power Systems),2010,34(24):36-39,45.

[3]陳加盛,張建華,林建業,等(Chen Jiasheng,Zhang Jianhua,Lin Jianye,et al).以降低電網損耗為目標的電動汽車充電策略(Strategies for electric vehicle charging with aiming at reducing network losses)[J].電力系統及其自動化學報(Proceedings of the CSU-EPSA),2012,24(3):139-144.

[4]杜愛虎,胡澤春,宋永華,等(Du Aihui,Hu Zechun,Song Yonghua,et al).考慮電動汽車充電站布局優化的配電網規劃(Distribution network planning considering layout optimization of electric vehicle charging stations)[J].電網技術(Power System Technology),2011,35(11):35-42.

[5]劉志鵬,文福栓,薛禹勝,等(Liu Zhipeng,Wen Fushuan,Xue Yusheng,et al).電動汽車充電站的最優選址和定容(Optimal siting and sizing of electric vehicle charging stations)[J].電力系統自動化(Automation of Electric Power Systems),2012,36(3):54-59.

[6]徐凡,俞國勤,顧臨峰,等(Xu Fan,Yu Guoqin,Gu Lin-feng,et al).電動汽車充電站布局規劃淺析(Tentative analysis of layout of electrical vehicle charging stations)[J].華東電力(East China Electric Power),2009,37(10):1678-1682.

[7]李菱,李燕青,姚玉海,等(Li Ling,Li Yanqing,Yao Yuhai,et al).基于遺傳算法的電動汽車充電站的布局規劃(Layout planning of electric vehicle charging stations based on genetic algorithm)[J].華東電力(East China Electric Power),2011,39(6):1004-1006.

[8]任玉瓏,史樂峰,張謙,等(Ren Yulong,Shi Lefeng,Zhang Qian,et al).電動汽車充電站最優分布和規模研究(Optimal distribution and scale of charging stations for electric vehicles)[J].電力系統自動化(Automation of Electric Power Systems),2011,35(14):53-57.

[9]高賜威,張亮,薛飛,等(Gao Ciwei,Zhang Liang,Xue Fei,et al).考慮集中型充電站定址分容的電網規劃研究(Grid planning considering capacity and site of largescale centralized charging stations)[J].中國電機工程學報(Proceedings of the CSEE),2012,32(7):40-46.

[10]盧艷霞,張秀敏,蒲孝文(Lu Yanxia,Zhang Xiumin,Pu Xiaowen).電動汽車充電站諧波分析(Harmonic study of electric vehicle chargers)[J].電力系統及其自動化學報(Proceedings of the CSU-EPSA),2006,18(3):51-54.

[11]Wang Hengsong,Huang Qi,Zhang Changhua,et al.A novel approach for the layout of electric vehicle charging station[C]//International Conference on Apperceiving Computing and Intelligence Analysis.Chengdu,China:2010.

[12]孔濤,程浩忠,李鋼,等(Kong Tao,Cheng Haozhong,Li Gang,et al).配電網規劃研究綜述(Review of power distribution network planning)[J].電網技術(Power System Technology),2009,33(19):92-99.

[13]王雷,顧潔(Wang Lei,Gu Jie).中壓配電網優化規劃的改進單親遺傳算法(Improved partheno-genetic algorithm for medium voltage distribution network optimal planning)[J].電力系統及其自動化學報(Proceedings of the CSU-EPSA),2006,18(3):72-76.

[14]陳根軍,李繼洸,王磊,等(Chen Genjun,Li Jiguang,Wang Lei,et al).基于Tabu搜索的配電網絡規劃(Distribution system planning by Tabu search approach)[J].電力系統自動化(Automation of Electric Power Systems),2001,25(7):40-44.

[15]王志剛,楊麗徙,陳根永(Wang Zhigang,Yang Lixi,Chen Genyong).基于蟻群算法的配電網網架優化規劃方法(Ant colony algorithm for distribution network planning)[J].電力系統及其自動化學報(Proceedings of the CSUEPSA),2002,14(6):73-76.

[16]張憲,宛津莎,楊薛明,等(Zhang Xian,Yuan Jinsha,Yang Xueming,et al).基于粒子群算法的配電網網架規劃(Distribution network planning based on particle swarm optimization)[J].華北電力大學學報(Journal of North China Electric Power University),2006,33(3):14-17.

[17]章文俊,程浩忠,王一,等(Zhang Wenjun,Cheng Haozhong,Wang Yi,et al).基于樹形結構編碼單親遺傳算法的配電網優化規劃(Distribution network optimal planning based on tree structure encoding partheno-genetic algorithm)[J].電工技術學報(Transactions of China Electrotechnical Society),2009,24(5):154-160.

[18]Jonnavithula S,Billinton R.Minimum cost analysis of feeder routing in distribution system planning[J].IEEE Trans on Power Delivery,1996,11(4):1935-1940.

[19]戈東方,鐘大文.電力工程電氣設計手冊[M].北京:中國電力出版社,1989.

Distribution Network Planning by Considering Siting and Sizing of Electric Vehicle Charging Stations

CHEN Ting1,WEI Zhi-nong1,WU Shuang1,SUN Guo-qiang1,WEI Yan-fang1,XU Xiao-hui2
(1.College of Energy and Electrical Engineering,Hohai University,Nanjing 210098,China;
2.China Electric Power Research Institute(Nanjing),Nanjing 210003,China)

The problem of distribution network planning by considering siting and sizing of electric vehicle charging stations has recently

the increasing research attention.In this paper,an optimal model is proposed,which aims to minimize the investment and operation costs of the distribution network and electric vehicle charging stations. By using the improved tree structure encoding partheno-genetic algorithm(ITSE-PGA)with elitist strategy,the problems of distribution network planning and optimization of the location and capacity of electric vehicle charging stations can be tackled simultaneously.Different from the ITSE-PGA for traditional distribution network planning,the considered distribution network here is not necessarily connected in the optimization process when considering electric vehicle charging stations,so the infeasible solutions should be specially handled.Finally,simulation results are provided to show the feasibility and effectiveness of the proposed method.

electric vehicle charging station;distribution network planning;partheno-genetic algorithm;tree structure encoding

TM71

A

1003-8930(2013)03-0001-07

陳婷(1989—),女,碩士研究生,研究方向為配電網規劃、電動汽車。Email:ctjr@163.com

2012-11-29;

2013-01-14

國家自然科學基金項目(51277052,51107032,61104045);國網公司總部科技項目:智能配用電的技術體系及仿真基礎性問題研究

衛志農(1962—),男,教授,博士生導師,研究方向為電力系統運行分析與控制、輸配電系統自動化等。Email:wzn_nj@263.net

吳霜(1989—),女,碩士研究生,研究方向為配電網規劃與自動化。Email:wushuanghehai@sina.com

猜你喜歡
配電網規劃優化
超限高層建筑結構設計與優化思考
房地產導刊(2022年5期)2022-06-01 06:20:14
民用建筑防煙排煙設計優化探討
關于優化消防安全告知承諾的一些思考
一道優化題的幾何解法
配電網自動化的應用與發展趨勢
規劃引領把握未來
快遞業十三五規劃發布
商周刊(2017年5期)2017-08-22 03:35:26
多管齊下落實規劃
中國衛生(2016年2期)2016-11-12 13:22:16
基于IEC61850的配電網數據傳輸保護機制
電測與儀表(2016年5期)2016-04-22 01:14:14
迎接“十三五”規劃
主站蜘蛛池模板: 日韩免费无码人妻系列| 扒开粉嫩的小缝隙喷白浆视频| 波多野结衣在线se| 粗大猛烈进出高潮视频无码| 扒开粉嫩的小缝隙喷白浆视频| 成人在线天堂| 97国产精品视频人人做人人爱| 性喷潮久久久久久久久| 国产成人综合网| 毛片在线看网站| 亚洲av成人无码网站在线观看| 国产天天色| 免费日韩在线视频| 免费看的一级毛片| 欧美综合区自拍亚洲综合天堂| 日韩毛片免费| 91久久国产成人免费观看| 19国产精品麻豆免费观看| 91网红精品在线观看| 欧美一区二区啪啪| 亚洲人成日本在线观看| 青青久在线视频免费观看| 任我操在线视频| 欧美色99| 午夜啪啪网| 久久人体视频| 福利视频99| 亚洲男人的天堂久久香蕉| 亚洲视频四区| 亚洲精品777| 精品人妻一区二区三区蜜桃AⅤ| 久久99国产乱子伦精品免| 精品伊人久久大香线蕉网站| 国产高潮流白浆视频| 亚洲综合香蕉| 一级毛片在线播放免费| 国产超碰一区二区三区| 2020精品极品国产色在线观看 | 美女无遮挡免费网站| 97久久超碰极品视觉盛宴| 午夜a视频| 日本亚洲国产一区二区三区| 国产黑丝视频在线观看| 久久无码免费束人妻| 99无码熟妇丰满人妻啪啪| 久久精品国产在热久久2019 | 白丝美女办公室高潮喷水视频| 亚洲国产成人久久77| 99精品热视频这里只有精品7 | 97超碰精品成人国产| 青青热久免费精品视频6| 日日摸夜夜爽无码| 91九色国产在线| 精品视频91| 伊人久综合| 国产99热| 女人毛片a级大学毛片免费| 国产精品黑色丝袜的老师| 鲁鲁鲁爽爽爽在线视频观看| 日韩不卡高清视频| 久久永久精品免费视频| 婷婷色狠狠干| 国产白浆在线观看| 亚洲看片网| a级毛片免费看| 国产麻豆精品在线观看| 中文字幕 91| 国产人人射| 在线网站18禁| 久久这里只有精品8| 国产成人福利在线视老湿机| 在线免费看片a| 亚洲一级毛片免费观看| 国产精品嫩草影院av| 亚洲日本中文字幕天堂网| 国产在线观看成人91| 成人韩免费网站| 在线观看91精品国产剧情免费| 中文字幕在线欧美| 四虎精品免费久久| 99热这里只有免费国产精品| 国产电话自拍伊人|