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

電動汽車充電站最優(yōu)規(guī)模和布局

2014-02-14 02:21:22王澤黎
電力建設(shè) 2014年4期
關(guān)鍵詞:優(yōu)化

王澤黎,齊 燦

(華北電力大學(xué)電氣與電子工程學(xué)院,北京市102206)

0 引 言

為減少霧霾天氣、降低PM2.5 濃度,發(fā)展節(jié)能、環(huán)保的交通工具越來越成為各界關(guān)注的焦點(diǎn)[1]。在此背景下,電動汽車產(chǎn)業(yè)得到飛速發(fā)展,銷售市場已初具規(guī)模。但是,對于電動汽車采用充電還是換電模式補(bǔ)充能量,國內(nèi)外一直爭議不斷。國家電網(wǎng)公司和南方電網(wǎng)公司推薦換電池為主的電動汽車運(yùn)營方案,制定電動汽車電池充電站的最優(yōu)規(guī)劃迫在眉睫。

基于電池租賃的換電池模式有如下優(yōu)點(diǎn):(1)電池由電網(wǎng)公司或?qū)I(yè)租賃公司購買,消費(fèi)者只需要交納租賃價(jià)格,成本大降,激發(fā)用戶購車積極性;(2)換電池有利于統(tǒng)一規(guī)范電池標(biāo)準(zhǔn)以及實(shí)現(xiàn)電池生產(chǎn)、管理、運(yùn)輸、維護(hù)等過程的規(guī)模效益;(3)電池更換既能保證電能迅速補(bǔ)給又不會損害電池,有利于延長電池的使用壽命[2];(4)電池充電站在滿足用戶對充電電池需求的同時(shí),也可以協(xié)助電網(wǎng)實(shí)現(xiàn)“削峰填谷”,提高負(fù)荷利用率。

文獻(xiàn)[3]根據(jù)地理因素和充電站服務(wù)半徑確定候選站址,同時(shí)以充電站的投資、運(yùn)行成本和網(wǎng)損費(fèi)用之和最小為目標(biāo),構(gòu)造電動汽車充電站規(guī)劃的最優(yōu)數(shù)學(xué)模型,并采用改進(jìn)的原對偶內(nèi)點(diǎn)求解。文獻(xiàn)[4]將充電者的充電時(shí)間限制在1個(gè)硬時(shí)間窗內(nèi),以充電站投資者投入成本最小為目標(biāo),建立了計(jì)及充電路徑選擇、充電電價(jià)的充電站規(guī)模選擇和布局的最優(yōu)模型,以滿足不同地點(diǎn)電動汽車動態(tài)充電需求。文獻(xiàn)[5]引入了對電動汽車分層、分區(qū)調(diào)度的理念,通過優(yōu)化各電動汽車代理商在各時(shí)段對充電負(fù)荷和放電出力的調(diào)度,使系統(tǒng)在目標(biāo)時(shí)間區(qū)間內(nèi)總負(fù)荷水平的方差最小,建立了基于雙層優(yōu)化的可入網(wǎng)電動汽車充放電調(diào)度模型。文獻(xiàn)[6]對電動汽車充電的負(fù)荷特性進(jìn)行分析,在簡化鋰電池I-U 充電模型的基礎(chǔ)上,得到單個(gè)電動汽車的充電模型,提出2 階段泊松分布的電動汽車充電站集聚模型和基于充電站日充電負(fù)荷曲線的電動汽車充電站負(fù)荷集聚模型的建模方法。

本文將網(wǎng)絡(luò)損耗和電力線路建設(shè)費(fèi)用引入電池充電站投資模型,建立一種考慮電動汽車時(shí)間-空間分布、電能損耗和新建線路費(fèi)用的電動汽車電池充電站選址定容新模型。模型用高峰期的交通狀況模擬電動汽車,最大程度地滿足用戶需求。在線路傳輸功率、無功補(bǔ)償上下限、站內(nèi)電池?cái)?shù)量和服務(wù)半徑的約束下,由新建站初始投資、網(wǎng)損費(fèi)用、新建線路費(fèi)用和行駛耗電成本構(gòu)成目標(biāo)函數(shù)。針對建立的模型,采用混沌粒子群優(yōu)化算法優(yōu)化充電站布局。

1 電池充電站的經(jīng)濟(jì)性模型

電池充電站規(guī)劃取決于其投資周期費(fèi)用和用戶便利性。投資周期費(fèi)用是指從開始考慮建設(shè)電池充電站,到電池充電站投入正常使用整個(gè)過程的費(fèi)用;用戶便利性以用戶在電池充電站的等待時(shí)間和用戶到電池充電站的路程耗時(shí)最短為目標(biāo)。等待時(shí)間最短意味著電池充電站的數(shù)量、位置和規(guī)模都能滿足用戶更換電池最高峰時(shí)刻的需求,減少甚至消除用戶排隊(duì)概率;排除交通擁堵和道路維修等不確定因素的影響,行駛耗時(shí)與路程一一對應(yīng),行駛路程越長,路上耗時(shí)越長,反之亦然。因此,路上行駛距離折算到耗電成本即可以作為衡量用戶便利性的指標(biāo)。本研究以最小化電池充電站投資費(fèi)用及行駛耗電成本為目標(biāo),確定電池充電站的規(guī)模、位置、服務(wù)范圍及用戶換電耗時(shí),其數(shù)學(xué)模型為

式中:G 為電池充電站年費(fèi)用;Gbj為電池充電站j 的投資成本;Gcj為用戶到電池充電站j 的行駛耗電成本,j=1,2,3…,n 分別表示電池充電站的序號;Gbij為電池充電站j 的第i 類投資成本,i 的取值為1,2,3 時(shí)分別表示電池充電站初始投資成本、網(wǎng)損成本、新建線路成本。

1.1 投資成本

1.1.1 初始投資成本

新建站j 的初始投資成本Gb1j包括購買設(shè)備的費(fèi)用Ej和購地費(fèi)用Aj。Ej包括交流充電樁、分箱充電機(jī)、電池更換系統(tǒng)、配電監(jiān)控、充電監(jiān)控、計(jì)費(fèi)系統(tǒng)及應(yīng)急充電機(jī)等設(shè)備購買費(fèi)用,設(shè)備的個(gè)數(shù)和容量與該電池充電站的規(guī)模有關(guān)。Aj由電池充電站的位置和占地面積決定。Gb1j的計(jì)算公式為

式中:r0表示投資回收率;m 表示電池充電站的運(yùn)行年限。

1.1.2 網(wǎng)損費(fèi)用

網(wǎng)損費(fèi)用Gb2j是指新建電池充電站j 接入原有電網(wǎng)后,新增網(wǎng)絡(luò)損耗引起的費(fèi)用。該費(fèi)用與電池充電站接入電網(wǎng)的位置有關(guān),是衡量電池充電站選址對電網(wǎng)經(jīng)濟(jì)性影響的指標(biāo)之一。Gb2j的計(jì)算公式為

式中:rL表示網(wǎng)損費(fèi)用的年值折算系數(shù);ΔAj表示新建電池充電站j 接入原有電網(wǎng)后的網(wǎng)損增加量。

1.1.3 新建線路費(fèi)用

新建線路費(fèi)用Gb3j是指新建電池充電站j 到距離最近的變電站的線路投資,該費(fèi)用不僅與電池充電站位置有關(guān)而且與規(guī)劃區(qū)域電力負(fù)荷級別有關(guān)。Gb3j的計(jì)算公式為

式中:λ1表示單位長度雙回線路的投資費(fèi)用;lj表示電池充電站j 到距離最近的變電站的線路長度。

1.2 行駛耗電成本

行駛耗電成本Gcj的計(jì)算公式為

式中:α 表示道路曲折系數(shù);β 表示道路暢通系數(shù);k 表示平均每輛車的充電次數(shù);ω 表示距離對電價(jià)的折算系數(shù);Dij表示第i 輛電動汽車到電池充電站j 的直線距離;Ij表示到電池充電站j 換電池的電動汽車集合。

1.3 約束條件

為了避免影響原有電網(wǎng)的安全性、電能質(zhì)量和滿足用戶需求,電池充電站規(guī)劃模型還需加入以下約束條件:

(1)電池充電站線路傳輸功率約束。

式中:Pj(t)表示電池充電站j 和變電站之間的連接線路在t 時(shí)段的傳輸功率;Pjmax(t)表示電池充電站j和變電站之間的連接線路在t 時(shí)段的最大傳輸功率。

(2)無功補(bǔ)償上下限約束。

式中:Qj表示電池充電站j 的無功補(bǔ)償功率;Qjmax和Qjmin分別表示電池充電站j 的無功補(bǔ)償功率上、下限。

(3)電池?cái)?shù)量約束。

式中:nj表示電池充電站j 的電池更換系統(tǒng)單位小時(shí)服務(wù)電池?cái)?shù)量;Tj表示電池充電站j 的日平均工作時(shí)間;Nj表示電池充電站j 的充電電池?cái)?shù)量。

(4)服務(wù)半徑約束。

式中:rj表示第j個(gè)電池充電站的服務(wù)半徑;esoc表示電動汽車的平均荷電狀態(tài);LN表示電動汽車單次額定行駛里程。

2 混沌粒子群算法

粒子群優(yōu)化(particle swarm optimization,PSO)算法模擬鳥群覓食行為的一種通過迭代實(shí)現(xiàn)的矩陣隨機(jī)搜索算法[7]。粒子更新自己的速度和位置的公式為:

式中:c1、c2為學(xué)習(xí)因子;w 為慣性權(quán)重;為第k次循環(huán)中粒子的位置;為第k次循環(huán)中粒子的速度;為粒子的全局最優(yōu)坐標(biāo);xg,j為粒子的局部最優(yōu)坐標(biāo)。PSO 算法使用了速度-位置同時(shí)搜索的模型,全局尋優(yōu)速度和能力較強(qiáng);但基本PSO 算法存在搜索結(jié)果精度不高、容易陷入局部解的缺點(diǎn)。為克服這些不足,將混沌運(yùn)動引入粒子群優(yōu)化算法,利用混沌運(yùn)動的遍歷性、隨機(jī)性等特點(diǎn),當(dāng)粒子陷入早熟收斂時(shí),用混沌擾動來跳出局部最優(yōu),并快速搜尋到最優(yōu)解,提高解的精度和收斂速度。根據(jù)“多點(diǎn)隨機(jī)全局搜索策略+面向問題的局部尋優(yōu)算法=最有效的全局優(yōu)化算法”[8],運(yùn)用基于混沌序列的粒子群優(yōu)化(chaos particle swarm optimization,CPSO)算法,得到理想結(jié)果。

混沌運(yùn)動的數(shù)學(xué)過程可描述為:對于一個(gè)隨機(jī)初始向量的N 維向量,x0=[x0,1,x0,2,…,x0,N]T,x0,n∈[0,1],根據(jù)某種產(chǎn)生混沌序列的模型,開始迭代。根據(jù)參考文獻(xiàn)[9],本文采用邏輯自映射函數(shù),得到迭代序列xm,n。其數(shù)學(xué)表達(dá)式為

式中:i=0,1,2…,M;j=0,1,…,N。

混沌運(yùn)動的步驟為:

(1)根據(jù)式(13),在[-1,1]內(nèi)隨機(jī)初始化,生成1個(gè)N×2 的矩陣;

(2)根據(jù)下式將混沌運(yùn)動的變量映射成種群粒子;

混沌變量xi,j變換到以粒子當(dāng)前全局最優(yōu)位置xg,j為圓心,Ri,j為半徑的圓上。Ri,j為各鄰域內(nèi)混沌運(yùn)動的半徑,所以通過邏輯自映射函數(shù)優(yōu)化后的種群粒子yi,j的取值范圍如下:

(3)計(jì)算目標(biāo)函數(shù)的適應(yīng)值f(xm,n),并更新全局最優(yōu)適應(yīng)值和全局最優(yōu)位置,更新在混沌迭代過程中的歷史最優(yōu)適應(yīng)值f(gbest)和歷史最優(yōu)位置xibest。如果混沌迭代后的全局最優(yōu)適應(yīng)值優(yōu)于歷史全局最優(yōu)適應(yīng)Gbest,則用混沌迭代后的位置和速度替換原有的歷史最優(yōu)位置和速度。其中,速度計(jì)算公式為

綜上所述,基于CPSO 算法的電池充電站規(guī)劃流程圖如圖1 所示。

圖1 基于CPSO 算法的電池充電站規(guī)劃流程Fig.1 Planning process of charging station based on CPSO

需要說明的是:

(1)Ni表示滿足規(guī)劃區(qū)用戶需求時(shí),建設(shè)規(guī)模為i 的充電站數(shù)量,Nmax表示充電站最大需求數(shù)量;Ni的計(jì)算公式為

式中:Cmax表示1 天內(nèi)最大充電需求時(shí)刻電動汽車的數(shù)目;α1表示同時(shí)工作的充電站比例;α2表示充電機(jī)的同時(shí)率;β 表示充電機(jī)的充電效率;cosφ 表示功率因數(shù);根據(jù)《國家電網(wǎng)公司電動汽車充電設(shè)施建設(shè)指導(dǎo)意見》,將電動汽車充電站按服務(wù)能力分成不同等級,i=1、2、3;li表示等級為i 的充電站的服務(wù)能力。

(2)判定粒子運(yùn)動是否停滯的標(biāo)準(zhǔn)為

式中:δ,NC為根據(jù)實(shí)際情況設(shè)定的常數(shù);NΔG的初始值設(shè)為0。當(dāng)不滿足ΔGi<δ 時(shí),NΔG=NΔG+1,如果不滿足式(18),則判定粒子運(yùn)動停滯,重新進(jìn)入混沌運(yùn)動過程。

3 算例分析

某小區(qū)總面積為70.3 km2,常住人口10.7 萬人,東西跨度11.06 km,南北跨度9.12 km;該小區(qū)主干道2條,次干道4條,支路20條,被主干道劃分為4個(gè)分區(qū);日均車流總量1.9 萬輛。該小區(qū)的工業(yè)和旅游業(yè)發(fā)達(dá),電動公交車和出租車占電動汽車總量的比重較大,具體基礎(chǔ)數(shù)據(jù)如表1 所示,電動汽車的位置重心和數(shù)量見附表。該小區(qū)的居民住宅、旅游區(qū)、商業(yè)區(qū)集中,各用地類型的地價(jià)如表2,不同規(guī)模的充電站設(shè)備的詳細(xì)配置如表3。由于現(xiàn)實(shí)中電動汽車類型、服務(wù)耗時(shí)、道路曲折系數(shù)、道路暢通系數(shù)不唯一,為簡化問題,本算例取平均值進(jìn)行計(jì)算,其中各參數(shù)的取值如下:r0=0.1;rl=0.2;λ1=1.65萬元/km;α=1.05 ~1.45;β=1.0 ~1.3;m=20;k =242次/輛;ω=1.1元/km;T =10 h/a;esoc=50% ~75%;LN=150 ~200 km。

表1 基礎(chǔ)數(shù)據(jù)Tab.1 Basic data

表2 土地單價(jià)Tab.2 Land price 元/m2

表3 不同規(guī)模電池充電站的詳細(xì)配置Tab.3 Detailed configurations of charging station with different scales

根據(jù)表1 的換電電池組需求,待建的充電站必須保證每天可以服務(wù)5 400 組電池,結(jié)合表2 和公式(17),當(dāng)電池充電站按規(guī)模1 配置時(shí),最少需要6 座;當(dāng)電池充電站按規(guī)模3 配置時(shí),最多需要18 座電池充電站;以n(6 ~18)為循環(huán)變量,計(jì)算當(dāng)前數(shù)量的電池充電站年總費(fèi)用,取最小值對應(yīng)的結(jié)果作為規(guī)模選擇和布局的最優(yōu)規(guī)劃結(jié)果。

理論上,初始設(shè)備配置是造成不同規(guī)模電池充電站投資費(fèi)用差異的主要因素,隨著電池充電站的增多,其規(guī)模逐漸減小,單個(gè)電池更換系統(tǒng)占設(shè)備投資的比重升高。而電池更換系統(tǒng)對初始建設(shè)費(fèi)用的影響最大,因此,電池充電站越多,總投資越高。對算例獨(dú)立運(yùn)行100次,計(jì)算出建設(shè)6 ~18 座電池充電站的最小費(fèi)用分別為3 070.2,3 076.4,3 062.7,3 055.9,3 063.1,3 070.5,3 078.4,3 084.6,3 091.0,3 097.1,3 110.0,3 122.3,3 135.1萬元。可見,在給定的條件下,規(guī)劃區(qū)內(nèi)建設(shè)9 座電池充電站時(shí)結(jié)果最優(yōu),優(yōu)化后的各充電站位置和規(guī)模如表4,電池充電站的各項(xiàng)投資費(fèi)用、用戶路上耗電成本和時(shí)間如表5,其中,用戶路上耗電時(shí)間是在假設(shè)電動汽車的行駛速度保持在100 km/h的基礎(chǔ)上得出的。結(jié)果表明,該算法優(yōu)化得到的結(jié)果既保證投資費(fèi)用最少,又能夠滿足用戶的便利性。圖2 為采用CPSO 算法和PSO 算法的最優(yōu)規(guī)劃的動態(tài)演化過程,圖3 為采用CPSO算法和PSO 算法的相鄰兩點(diǎn)的距離變化軌跡。從圖2、3 可以看出:CPSO 算法在避免陷入早熟、收斂速度、最優(yōu)結(jié)果等三方面均優(yōu)于PSO 算法。結(jié)果表明,本文建立的電池充電站規(guī)模選擇和布局模型一方面降低了建設(shè)成本,縮小了投資;另一方面又能滿足不同交通密度的換電需求,符合預(yù)期設(shè)想。

表4 充電站位置的優(yōu)化結(jié)果Tab.4 Optimization results of charging station location

表5 電池充電站優(yōu)化的仿真結(jié)果Tab.5 Simulation results of charging station optimization

圖2 CPSO 和PSO 算法的最優(yōu)規(guī)劃的動態(tài)演化過程Fig.2 Optimal planning of dynamic evolution process of CPSO and PSO

圖3 CPSO 和PSO 算法的相鄰兩點(diǎn)的距離變化軌跡Fig.3 Change trajectory between two adjacent points in CPSO and PSO

4 結(jié) 論

(1)本文提出的數(shù)學(xué)模型中,既考慮了投資的費(fèi)用,又加入了用戶便利性,滿足了建設(shè)方和使用方的需求,從全面性和經(jīng)濟(jì)性上都反映了充電站規(guī)劃問題的實(shí)質(zhì)。

(2)CPSO 算法利用混沌對初值的敏感性和混沌遍歷性初始化種群;再通過邏輯自映射函數(shù),形成混沌序列和優(yōu)化變量取值范圍的對應(yīng)關(guān)系。通過混沌初始化可以使PSO 算法能從較好的初始值開始尋優(yōu),同時(shí),迭代更新過程的停滯標(biāo)準(zhǔn)使搜索更為精細(xì),克服了PSO 算法容易早熟、陷入局部極值的缺陷。

(3)對某小區(qū)電動汽車電池充電站進(jìn)行了規(guī)劃,證明了本文提出的優(yōu)化計(jì)算模型是科學(xué)、可行的。

[1]徐敬,丁國安,顏鵬,等. 北京地區(qū)PM_(2.5)的成分特征及來源分析[J].應(yīng)用氣象學(xué)報(bào),2007,18(5):645-654.

[2]高賜威,張亮,薛飛,等.考慮集中型充電站定址分容的電網(wǎng)規(guī)劃研究[J].中國電機(jī)工程學(xué)報(bào),2012,32(7):40-46.

[3]任玉瓏,史樂峰,張謙,等.電動汽車充電站最優(yōu)分布和規(guī)模研究[J].電力系統(tǒng)自動化,2011,35(14):53-57.

[4]何永秀,周波,熊威.北京市電動汽車用戶響應(yīng)度研究[J].電力建設(shè),2013,34(1):1-6.

[5]Liu Z P,Wen F S,Ledwich G. Optimal planning of electric-vehicle charging stations in distribution systems[J]. IEEE Transactions on Power Delivery,2013,28(1):102-110.

[6]胡澤春,宋永華,徐智威,等. 電動汽車接入電網(wǎng)的影響與利用[J].中國電機(jī)工程學(xué)報(bào),2012,32(4):1-10.

[7]姚偉鋒,趙俊華,文福拴,等. 基于雙層優(yōu)化的電動汽車充放電調(diào)度策略[J].電力系統(tǒng)自動化,2012,36(11):30-37.

[8]He J,Zhou B,F(xiàn)eng C,et al. Electric vehicle charging station planning based on multiple-population hybrid genetic algorithm[C]// Control Engineering and Communication Technology 2012 International Conference. Liaoning University,Shenyang:IEEE,2012:403-406.

[9]紀(jì)震,廖惠連,吳青華. 粒子群算法及應(yīng)用[M]. 北京:科學(xué)出版社,2009:28-45.

[10]劉長平,葉春明.基于邏輯自映射的變尺度混沌粒子群優(yōu)化算法[J].計(jì)算機(jī)應(yīng)用研究,2011,28(8):25-27.

[11]張樂平,許燕灝,胡紅,等. 電動汽車對電網(wǎng)影響的評價(jià)指標(biāo)與評估方法[J].電力建設(shè),2013,34(8):46-51.

[12]李菱,李燕青,姚玉海,等.基于遺傳算法的電動汽車充電站的布局規(guī)劃[J].華東電力,2011,39(6):1004-1006.

[13]劉自發(fā),張偉,王澤黎.基于量子粒子群優(yōu)化算法的城市電動汽車充電站優(yōu)化布局[J]. 中國電機(jī)工程學(xué)報(bào),2012,32(20):39-45.

猜你喜歡
優(yōu)化
超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
PEMFC流道的多目標(biāo)優(yōu)化
能源工程(2022年1期)2022-03-29 01:06:28
民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
關(guān)于優(yōu)化消防安全告知承諾的一些思考
一道優(yōu)化題的幾何解法
由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
圍繞“地、業(yè)、人”優(yōu)化產(chǎn)業(yè)扶貧
事業(yè)單位中固定資產(chǎn)會計(jì)處理的優(yōu)化
4K HDR性能大幅度優(yōu)化 JVC DLA-X8 18 BC
幾種常見的負(fù)載均衡算法的優(yōu)化
電子制作(2017年20期)2017-04-26 06:57:45
主站蜘蛛池模板: a网站在线观看| 亚洲欧美日韩中文字幕一区二区三区 | 国产精品免费久久久久影院无码| 久久96热在精品国产高清| 91丨九色丨首页在线播放 | 国产性生交xxxxx免费| 国内精品伊人久久久久7777人| 成人年鲁鲁在线观看视频| 精品一区国产精品| 国产传媒一区二区三区四区五区| 国产日本一线在线观看免费| 国产成人亚洲精品蜜芽影院| 欧美激情综合一区二区| 日韩成人午夜| 五月婷婷综合色| 国产18在线| 亚洲天堂高清| 97在线观看视频免费| 国内精品91| 玖玖精品视频在线观看| jizz在线观看| 日韩精品一区二区三区大桥未久 | 国产va视频| 91精品啪在线观看国产| 久久无码av三级| 欧美精品亚洲精品日韩专区va| 亚洲IV视频免费在线光看| 3D动漫精品啪啪一区二区下载| 亚洲男女在线| 制服丝袜一区| 欧美国产在线精品17p| 国产99热| 国产幂在线无码精品| 全裸无码专区| 久久免费视频播放| 先锋资源久久| 91一级片| 亚洲精品无码在线播放网站| 亚洲视频影院| 综合成人国产| 国产偷国产偷在线高清| 欧美精品一区在线看| 这里只有精品在线| 自拍偷拍一区| 超级碰免费视频91| 中文字幕亚洲另类天堂| 亚洲综合色吧| 美女免费黄网站| 国产自在线拍| 国产日韩精品欧美一区喷| 日本黄色不卡视频| 三上悠亚精品二区在线观看| 亚洲第一综合天堂另类专| 日韩精品一区二区三区大桥未久 | 日韩视频福利| 国产xx在线观看| 亚洲αv毛片| 久久亚洲黄色视频| 国产一区二区精品福利| 麻豆国产精品| 亚洲男女在线| 久久精品aⅴ无码中文字幕| 亚洲综合久久一本伊一区| 免费看的一级毛片| 中文字幕欧美日韩| 亚洲区一区| 9啪在线视频| 午夜精品久久久久久久99热下载| 99激情网| 91福利片| 超清无码熟妇人妻AV在线绿巨人| 好久久免费视频高清| 九九九国产| 18禁黄无遮挡网站| 久久精品国产在热久久2019| 欧美日韩动态图| 欧美成人区| 亚洲欧美成人在线视频| 亚洲一区国色天香| 就去吻亚洲精品国产欧美| 丁香婷婷综合激情| 日韩毛片免费观看|