亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        梯級(jí)泵站輸水系統(tǒng)旬優(yōu)化調(diào)度及經(jīng)濟(jì)運(yùn)行研究

        2018-09-10 01:19:12吳怡李智王京晶曹悅雷曉輝
        人民黃河 2018年5期
        關(guān)鍵詞:優(yōu)化調(diào)度經(jīng)濟(jì)運(yùn)行

        吳怡 李智 王京晶 曹悅 雷曉輝

        摘要:南水北調(diào)來水調(diào)入密云水庫調(diào)蓄工程是北京南水北調(diào)配套工程的重要組成部分,該工程后半段從懷柔水庫通過3級(jí)泵站提升輸水至密云水庫。為實(shí)現(xiàn)該工程高效經(jīng)濟(jì)運(yùn)行,基于北京市分時(shí)電價(jià),以懷柔水庫至密云水庫梯級(jí)泵站輸水系統(tǒng)運(yùn)行單位輸水成本最小為目標(biāo),考慮懷柔水庫蓄水量、輸水線路水頭損失、流量平衡和總揚(yáng)程等約束條件,采用大系統(tǒng)分解協(xié)調(diào)模型求解該工程旬內(nèi)優(yōu)化調(diào)度及經(jīng)濟(jì)運(yùn)行問題。將該大系統(tǒng)分解協(xié)調(diào)模型分解為梯級(jí)泵站揚(yáng)程優(yōu)化分配模型和旬優(yōu)化調(diào)度模型2個(gè)子系統(tǒng),采用動(dòng)態(tài)規(guī)劃法實(shí)現(xiàn)各子系統(tǒng)的優(yōu)化,從而實(shí)現(xiàn)大系統(tǒng)的全局優(yōu)化。應(yīng)用情況表明:模型具有較好的適用性,可為工程優(yōu)化調(diào)度提供依據(jù),指導(dǎo)工程高效經(jīng)濟(jì)運(yùn)行。

        關(guān)鍵詞:優(yōu)化調(diào)度;經(jīng)濟(jì)運(yùn)行;梯級(jí)泵站;大系統(tǒng)協(xié)調(diào)分解模型;動(dòng)態(tài)規(guī)劃法

        中圖分類號(hào):TV122+.5 文獻(xiàn)標(biāo)志碼:A doi:10.3969/j.issn.1000-1379.2018.05.031

        南水北調(diào)來水調(diào)入密云水庫調(diào)蓄工程是北京南水北調(diào)配套工程的重要組成部分,對消納南水北調(diào)來水,彌補(bǔ)水資源虧空,實(shí)現(xiàn)北京水資源優(yōu)化配置具有重要作用。該工程分為團(tuán)城湖至懷柔水庫段(前半段)和懷柔水庫至密云水庫段(后半段)兩段,分別通過6級(jí)和3級(jí)泵站串聯(lián)輸水。后半段梯級(jí)泵站輸水系統(tǒng)中間沒有調(diào)蓄工程,總揚(yáng)程大,是由泵站、輸水線路、閘門等構(gòu)成的復(fù)雜輸水系統(tǒng),其優(yōu)化調(diào)度與經(jīng)濟(jì)運(yùn)行面臨著懷柔水庫(水源)蓄量、渠道水頭損失、流量平衡、總揚(yáng)程等約束條件。

        對于梯級(jí)泵站輸水系統(tǒng)優(yōu)化調(diào)度與經(jīng)濟(jì)運(yùn)行問題,國內(nèi)外學(xué)者通常將泵站和輸水管渠作為整體研究對象,根據(jù)系統(tǒng)分析理論,將泵站實(shí)際運(yùn)行參數(shù)(或泵站機(jī)組性能測算結(jié)果)與渠道(或管道)運(yùn)行水力學(xué)特性參數(shù)結(jié)合起來[1-5],運(yùn)用遺傳算法[6]、動(dòng)態(tài)規(guī)劃法[7]、粒子群優(yōu)化算法[8-10]、大系統(tǒng)分解協(xié)調(diào)法[11]等優(yōu)化算法進(jìn)行求解。桑國慶[12]以梯級(jí)泵站輸水系統(tǒng)為研究對象、中長期運(yùn)行費(fèi)用最小為目標(biāo)函數(shù),綜合考慮泵站站前站后水位、輸水總量、分時(shí)電價(jià)等因素,建立梯級(jí)泵站系統(tǒng)中長期運(yùn)行費(fèi)用優(yōu)化模型,確定梯級(jí)泵站系統(tǒng)中長期運(yùn)行優(yōu)化方案。這些研究沒有全面考慮分時(shí)電價(jià)、水源蓄量、輸水線路水頭損失、流量平衡、總揚(yáng)程等約束條件,無法直接應(yīng)用于后半段梯級(jí)泵站輸水系統(tǒng)優(yōu)化調(diào)配過程。筆者通過建立旬優(yōu)化調(diào)度與經(jīng)濟(jì)運(yùn)行模型,將上述約束條件考慮在內(nèi),以后半段梯級(jí)泵站輸水系統(tǒng)單位輸水成本最小為目標(biāo),實(shí)現(xiàn)該工程的經(jīng)濟(jì)運(yùn)行,最大程度發(fā)揮工程效益。

        1 研究區(qū)域

        南水北調(diào)來水調(diào)人密云水庫調(diào)蓄工程后半段為本次研究區(qū)域。懷柔水庫在南水北調(diào)來水調(diào)人密云水庫調(diào)蓄工程中發(fā)揮旬調(diào)節(jié)功能,最小蓄水量為850萬m3,最大蓄水量為1.44億m3,懷柔水庫在調(diào)度過程中不允許產(chǎn)生棄水。懷柔水庫至密云水庫輸水線路總長約30km,總揚(yáng)程114.54m,泵站和閘門等建筑物布置見圖1。經(jīng)懷柔水庫調(diào)節(jié)部分水回補(bǔ)密懷順?biāo)吹兀渲幸徊糠謱Τ卑缀舆M(jìn)行生態(tài)補(bǔ)水及向水廠供水;另一部分通過懷柔水庫進(jìn)水閘旁新建的郭家塢泵站提升,經(jīng)京密引水渠反向輸水至北臺(tái)上倒虹吸,再經(jīng)新建的雁棲泵站加壓,由京密引水渠側(cè)新建DN2600PCCP輸水管道人白河電站下游調(diào)節(jié)池,最后由新建溪翁莊泵站加壓后通過白河發(fā)電洞將水加壓送人密云水庫,設(shè)計(jì)輸水流量10m3/s,在實(shí)際運(yùn)行中,輸水流量只能為10m3/s或0。

        2 數(shù)學(xué)模型

        該工程旬優(yōu)化調(diào)度與經(jīng)濟(jì)運(yùn)行問題屬于大系統(tǒng)優(yōu)化問題,采用大系統(tǒng)分解協(xié)調(diào)模型求解。把該大系統(tǒng)分解成2個(gè)子系統(tǒng),采用動(dòng)態(tài)規(guī)劃法實(shí)現(xiàn)各子系統(tǒng)的優(yōu)化。根據(jù)整個(gè)大系統(tǒng)的總目標(biāo),考慮各子系統(tǒng)之間的關(guān)聯(lián),實(shí)現(xiàn)大系統(tǒng)的全局優(yōu)化。大系統(tǒng)分解協(xié)調(diào)模型見圖2,包括2層:第一層為梯級(jí)泵站揚(yáng)程優(yōu)化分配模型,將總揚(yáng)程在梯級(jí)泵站之間進(jìn)行分配,使得梯級(jí)泵站運(yùn)行費(fèi)用最小,其中水頭損失通過建立一維水力學(xué)模型求得;第二層為旬優(yōu)化調(diào)度模型,優(yōu)化懷柔水庫調(diào)度過程,使得梯級(jí)泵站輸水系統(tǒng)單位輸水成本最小。

        2.1 一維水力學(xué)仿真模型

        針對梯級(jí)泵站輸水系統(tǒng)的特點(diǎn),構(gòu)建一維非恒定流水力學(xué)仿真模型,基于仿真模型,對泵站、倒虹吸、漸變段等復(fù)雜的內(nèi)部構(gòu)筑物進(jìn)行概化處理,并將概化的內(nèi)部構(gòu)筑物與圣維南方程組進(jìn)行耦合,同時(shí)采用穩(wěn)定性好、計(jì)算精度高的四點(diǎn)時(shí)空偏心Preissmann格式對方程組進(jìn)行離散,用高效率的雙掃描法求解。在此基咄上采用恒定非均勻流計(jì)算模塊對研究區(qū)恒定流和非恒定流進(jìn)行模擬,從而進(jìn)行渠道水力學(xué)特性分析和梯級(jí)泵站優(yōu)化[13]。輸水線路非恒定流計(jì)算采用一維圣維南方程組,該方程組由連續(xù)性方程和動(dòng)量方程組成:式中;A、Q為斷面面積和流量;x、t為空間和時(shí)間坐標(biāo);q為單位長度渠道上的側(cè)向入流流量;α為動(dòng)量修正系數(shù);Z為水位;g為重力加速度;Sf為水力坡度;S0為渠道底坡;K為流量模數(shù)。

        2.2 梯級(jí)泵站揚(yáng)程優(yōu)化分配模型

        最后一級(jí)泵站站后水位減去第一級(jí)泵站站前水位可得到梯級(jí)泵站系統(tǒng)凈揚(yáng)程,加上由水力學(xué)模型計(jì)算得到的各渠段水頭損失即可得到總揚(yáng)程[14]。在滿足各渠段水力聯(lián)系、各級(jí)泵站站前站后水位約束等條件下,進(jìn)行梯級(jí)泵站揚(yáng)程分配,使梯級(jí)泵站耗能最小。這是一個(gè)總揚(yáng)程在梯級(jí)泵站之間分配的空間優(yōu)化問題。將所研究的問題賦予時(shí)間特性,即將泵站按順序編號(hào),把每級(jí)泵站作為一個(gè)階段,該優(yōu)化問題就變成了一個(gè)多階段決策過程的優(yōu)化問題,可以采用動(dòng)態(tài)規(guī)劃技術(shù)來求解。

        (1)階段變量。將泵站編號(hào)作為階段變量j(j=1,2,…,m),其中m為投入運(yùn)行的泵站數(shù)(本研究中m取3)。

        (2)狀態(tài)變量。將第j階段至最末階段m的累計(jì)揚(yáng)程作為狀態(tài)變量。式中:sj為第j階段至最末階段m的累計(jì)揚(yáng)程(狀態(tài)變量);Hi為第i階段的揚(yáng)程;i為編號(hào)。

        (3)決策變量。將每級(jí)泵站的揚(yáng)程Hj作為決策變量。對于決策變量的離散化,離散的步長越小,計(jì)算精度越高,但計(jì)算量顯著增加。

        (4)狀態(tài)轉(zhuǎn)移方程。狀態(tài)轉(zhuǎn)移方程表示梯級(jí)泵站中第j+1階段(即第j+1個(gè)泵站)狀態(tài)變量Sj+1和第j階段(即第j個(gè)泵站)狀態(tài)變量Sj與決策變量Hj之間的關(guān)系。

        Sj+1=Sj-Hj(5)式中:S1=H*(H*為梯級(jí)泵站總揚(yáng)程);Sm+1=0。

        (5)目標(biāo)函數(shù)。對于梯級(jí)泵站,考慮各泵站的電費(fèi),其目標(biāo)函數(shù)為式中:Qj為第j個(gè)泵站的出水流量;ηj為第j個(gè)泵站在Qj、Hj組合工況下的效率;△t為梯級(jí)泵站運(yùn)行時(shí)間,取24h;c為日平均電價(jià);ρ為水的密度。

        (6)約束條件。總揚(yáng)程約束:

        單級(jí)泵站揚(yáng)程約束:

        Hjmin≤Hj≤Hjmax(8)式中;Zm為最后一級(jí)泵站站后水位;Z0為第一級(jí)泵站站前水位;hj,j+1為相鄰兩級(jí)泵站間輸水線路的水頭損失,由水力學(xué)模型計(jì)算得出;Hjmin、Hjmax為第j級(jí)泵站的最小、最大揚(yáng)程。

        站前站后水位約束:每一級(jí)泵站的站前站后水位都要在控制范圍內(nèi)。

        (7)求解方法。采用逆向遞推、正向決策法計(jì)算,遞推方程為

        Fm+1(Sm+1)=0(9)

        Fj*(Si)=min[Lj(Si,Hj)+Fj+1*(Sj+1)](10)式中:Fm+1(Sm+1)為最末階段結(jié)束時(shí)的電費(fèi);Lj(Si,Hj)為第j階段的費(fèi)用函數(shù);Fj*(sj)為第j階段的最小費(fèi)用函數(shù)。

        2.3 旬優(yōu)化調(diào)度

        根據(jù)懷柔水庫來水和回補(bǔ)地下水、供水、生態(tài)補(bǔ)水需求,在工程各約束條件下,優(yōu)化旬調(diào)度過程,使得梯級(jí)泵站輸水系統(tǒng)單位輸水成本最小。采用動(dòng)態(tài)規(guī)劃法描述此模型。

        (1)階段變量。階段變量采用k(k=1,2,…,N)表示,將1旬分為N個(gè)時(shí)段,即N天。

        (2)狀態(tài)變量。將懷柔水庫每時(shí)段初的蓄水量Vk作為狀態(tài)變量。

        (3)決策變量。將每時(shí)段懷柔水庫的調(diào)水流量Qk作為決策變量。

        (4)狀態(tài)轉(zhuǎn)移方程。由水量平衡可得狀態(tài)轉(zhuǎn)移方程:

        Vk+1=Vk+(Q1In+Q2In-Q1Out-Q2Out-Q3Out)×86400-Qk×Δtk(11)式中:Q1In為天然來水流量;Q2In為團(tuán)城湖來水流量;Q1Out為懷柔水庫向水廠供水流量;Q2Out為回補(bǔ)地下水流量(密懷順地下水源地,設(shè)計(jì)要求為10m3/s);Q3Out為生態(tài)補(bǔ)水流量; Qk為日平均調(diào)水流量Δtk為第k天的調(diào)水時(shí)間。

        此外,可通過Vk查水位庫容曲線得到第k天的懷柔水庫水位。

        (5)目標(biāo)函數(shù)。目標(biāo)函數(shù)為式中:Qk,j,p為第k天第p時(shí)段第j座泵站的流量,實(shí)際運(yùn)行中為10m3/s或0;Hk,j為第k天第j座泵站的揚(yáng)程;η(k,j)為第k天第j座泵站的效率;△tp為某天第p

        (6)約束條件。懷柔水庫蓄水量約束:Vmin≤Vk≤Vmax,Vmin和Vmax分別為懷柔水庫的最小和最大蓄水量;流量平衡約束:Qk,j,p=Qk;效率約束:η(k,j)≠0;揚(yáng)程約束和水位約束同前所述。

        (7)求解方法。采用逆向遞推、正向決策法計(jì)算,遞推方程為

        FN+1(VN+1)=0(13)

        Fk*(Vk)=min[Lk(Vk,Qk)+Fk+1*(Vk+1)](14)式中:FN+1(VN+1)為最末階段結(jié)束時(shí)的電費(fèi);Lk(Vk,Qk)為k階段的費(fèi)用函數(shù);Fk*(Vk)為k階段的最小費(fèi)用函數(shù)。

        2.4 優(yōu)化結(jié)果

        第一級(jí)泵站(郭家塢)站前水位為懷柔水庫水位,由水量平衡求出旬調(diào)度過程逐日蓄水量,反查水位庫容曲線得到懷柔水庫逐日水位且認(rèn)為一天內(nèi)不變;最后一級(jí)泵站(溪翁莊)站后水位為密云水庫水位,認(rèn)為一旬內(nèi)不變。根據(jù)工程實(shí)際運(yùn)行需求,基于北京市現(xiàn)行分時(shí)電價(jià),將1d概化為3個(gè)調(diào)度時(shí)段,即7:00至15:00、15:00至23;00和23;00至次日7;00,電價(jià)分別為1.1412、1.0205、0.3818元/(kW·h)。

        假設(shè)懷柔水庫在某一旬(10d)的調(diào)度過程中,團(tuán)城湖向懷柔水庫輸水流量為20m3/s,懷柔水庫天然來水流量為20m3/s,回補(bǔ)密懷順?biāo)吹亓髁繛?0m3/s,供水流量為18m3/s,生態(tài)補(bǔ)水流量為5m3/s,懷柔水庫初、末水位均為58m,密云水庫水位為132m。為方便計(jì)算,以及使計(jì)算結(jié)果便于運(yùn)行管理,泵站揚(yáng)程離散步長為0.1m,假設(shè)每天調(diào)水小時(shí)數(shù)為整數(shù),一旬中每天調(diào)水時(shí)間前9d為16h,第10d為24h,前9d的調(diào)水時(shí)間為15:00至次日7:00、第10d全天輸水。

        懷柔水庫水位見圖3,可知旬內(nèi)懷柔水庫水位變化很小,最大不足4cm。由水力學(xué)仿真模型計(jì)算出輸水線路水頭損失為0.35m,可認(rèn)為旬內(nèi)總揚(yáng)程為定值,從而計(jì)算出郭家塢、雁棲和溪翁莊泵站旬內(nèi)的揚(yáng)程分別為2.3、33.8、48.4m,3個(gè)泵站運(yùn)行效率分別為0.698、0.900、0.900。由此得出旬調(diào)水總量為604.8萬m3,總電費(fèi)為942605元,單位輸水成本為0.1559元/m3。

        3 結(jié)論

        針對南水北調(diào)來水調(diào)入密云水庫調(diào)蓄工程懷柔水庫至密云水庫段工程特點(diǎn),考慮懷柔水庫蓄水量約束、輸水損失、流量平衡、梯級(jí)泵站總揚(yáng)程等約束條件,以梯級(jí)泵站輸水系統(tǒng)單位輸水成本最小為優(yōu)化目標(biāo),基于北京市分時(shí)電價(jià),采用大系統(tǒng)分解協(xié)調(diào)模型分層建立旬優(yōu)化調(diào)度與經(jīng)濟(jì)運(yùn)行模型,包括梯級(jí)泵站揚(yáng)程優(yōu)化分配模型和旬優(yōu)化調(diào)度模型,并運(yùn)用動(dòng)態(tài)規(guī)劃法求解,其中輸水線路水頭損失由一維水力學(xué)模型計(jì)算。分層結(jié)構(gòu)模型具有較大的靈活性和適應(yīng)性,各子模型單獨(dú)優(yōu)化計(jì)算,層層遞進(jìn),可以保證達(dá)到全局最優(yōu)。從優(yōu)化結(jié)果可以看出,一旬內(nèi)懷柔水庫水位變化較小,各泵站每天運(yùn)行時(shí)間內(nèi)揚(yáng)程和效率均不變,且梯級(jí)泵站優(yōu)先在電價(jià)低的時(shí)候運(yùn)行,保證單位輸水成本最小。

        參考文獻(xiàn):

        [1]嚴(yán)登豐.泵與泵裝置特性預(yù)測[J].排灌機(jī)械工程學(xué)報(bào),2012,30(3):315-323.

        [2]賈仁甫,王紅,金明宇,等.調(diào)水工程中梯級(jí)泵站的優(yōu)化調(diào)度研究[J].揚(yáng)州大學(xué)學(xué)報(bào)(自然科學(xué)版),2006,9(2):69-73.

        [3]MORADI-JALAL M,MARINO M A,AFSHAR A.OptimalDesign and Operation of Irrigation Pumping Stations[J].Journal of Irrigation and Drainage Engineering,2003,129(3):149-154.

        [4]RODIN S I.Use of Genetic Algorithms for Optimal Control ofBulk Water Supply[J].Journal of Irrigation and DrainageEngineering,2004,130(5):357-365.

        [5]PULIDO-CALVO I,ROLDANJ,L6PEZ-LUQUE R,et al.Demand Forecasting for Irrigation Water Distribution Systems[J].Journal of Irrigation and Drainage Engineering,2003,129(6):422-431.

        [6]SAVIC D A,WALTERS G A,SMITH M R,et al.Cost Sav-ings on Large Water Distribution Systems:Design ThroughGenetic Algorithm Optimization[C]//Proceedings of JointConference on Water Resource Engineering and Water Re-sources Planning and Management.Minneapolis,Minnesota:ASCE,2000:1-10.

        [7]COSTA B E,ALMEIDA R A,VIANA A N C.Optimizationof Parallel Variable-Speed-Driven Centrifugal Pumps Operation[J].Energy Efficiency,2008(3):167-173.

        [8]OSTFELD A,TUBALTZEV A.Ant Colony Optimization forLeast-Cost Design and Operation of Pumping Water Distribu-tion Systems[1].Journal of Water Resources Planning andManagement,2008,134(2):107-118.

        [9]吳月秋,紀(jì)昌明,王麗萍,等.基于混沌粒子群算法的水電站水庫優(yōu)化調(diào)度[J].人民黃河,2008,30(11):96-112.

        [10]趙曉軍,田富強(qiáng),胡和平.粒子群優(yōu)化算法在水量調(diào)度方案優(yōu)化中的應(yīng)用[J].人民黃河,2005,27(11)26-27.

        [11]朱勁木,龍新平,劉梅清,等.東深供水工程梯級(jí)泵站的優(yōu)化調(diào)度[J].水力發(fā)電學(xué)報(bào),2005,24(3):123-127.

        [12]桑國慶.基于動(dòng)態(tài)平衡的梯級(jí)泵站輸水系統(tǒng)優(yōu)化運(yùn)行及控制研究[D].濟(jì)南:山東大學(xué),2012:142-149.

        [13]張大偉.南水北調(diào)中線干線水質(zhì)水量聯(lián)合調(diào)控關(guān)鍵技術(shù)研究[D].上海:東華大學(xué),2014:11-32.

        [14]劉波波.梯級(jí)泵站調(diào)水工程優(yōu)化運(yùn)行模擬研究[D].北京:北京工業(yè)大學(xué),2015:59-61.

        猜你喜歡
        優(yōu)化調(diào)度經(jīng)濟(jì)運(yùn)行
        一季度農(nóng)業(yè)農(nóng)村經(jīng)濟(jì)運(yùn)行良好
        2019全球經(jīng)濟(jì)運(yùn)行的三個(gè)特點(diǎn)
        中國外匯(2019年6期)2019-07-13 05:44:02
        2018年1—12月成都市經(jīng)濟(jì)運(yùn)行情況
        先鋒(2019年2期)2019-03-27 09:31:22
        對新醫(yī)改下的醫(yī)院經(jīng)濟(jì)運(yùn)行管理探討
        電力系統(tǒng)優(yōu)化調(diào)度及決策方法研究
        基于多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度的微電網(wǎng)模型及方法研究
        淺談含風(fēng)電裝機(jī)電力系統(tǒng)的發(fā)電調(diào)度
        智慧燃?xì)饩C合管理平臺(tái)的探討
        2016年1—11月全市經(jīng)濟(jì)運(yùn)行簡況
        投資北京(2017年1期)2017-02-13 20:44:17
        供水工作節(jié)能降耗思路探析
        北極光(2016年4期)2016-06-06 14:55:24
        北条麻妃在线视频观看| 风韵人妻丰满熟妇老熟| 亚洲精品久久国产精品| 国产女人高潮叫床免费视频| 亚洲欧美日韩国产综合一区二区| 亚洲高清有码在线观看| 97中文字幕一区二区| 国产三级视频不卡在线观看| 未满十八勿入av网免费| 又黄又爽的成人免费视频 | 国产成人精品久久一区二区三区 | 巨爆乳中文字幕爆乳区| 国产一区二区精品久久呦| 日韩精品乱码中文字幕| 狠狠噜狠狠狠狠丁香五月| 精品国产18久久久久久| 人妖精品视频在线观看| 青青草视频免费在线播放| 狂猛欧美激情性xxxx大豆行情| 尤物网址在线观看| 亚洲无码一二专区| 国产av乳头久久一区| 国产亚洲精品在线视频| 东京道一本热中文字幕| 免费啪啪视频一区| 精品国产自拍在线视频| 日本美女在线一区二区| 亚洲成a∨人片在线观看不卡| 中文字幕高清在线一区二区三区| 国产V亚洲V天堂A无码| 精品国产色哟av一区二区三区| 亚无码乱人伦一区二区| 精品麻豆国产色欲色欲色欲www| 国产av天堂成人网| 午夜精品一区二区三区视频免费看| 亚洲高清在线天堂精品| 国产精品_国产精品_k频道w| 无码日韩人妻AV一区免费| 日本人妻三级在线观看| 成人欧美一区二区三区黑人| 亚洲精华国产精华液的福利|