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

    考慮負(fù)荷裕度的區(qū)域綜合能源系統(tǒng)儲(chǔ)能雙層優(yōu)化配置

    2022-07-20 01:44:26劉書琪賴柏希金之儉
    電力自動(dòng)化設(shè)備 2022年7期
    關(guān)鍵詞:裕度熱網(wǎng)時(shí)段

    劉書琪,顧 潔,賴柏希,金之儉

    (上海交通大學(xué)電子信息與電氣工程學(xué)院,上海 200240)

    0 引言

    綜合能源系統(tǒng)IES(Integrated Energy System)作為能源互聯(lián)網(wǎng)的重要實(shí)現(xiàn)形式,將是未來(lái)電力系統(tǒng)發(fā)展的重要方向之一。在運(yùn)行時(shí),IES需要預(yù)留一定的負(fù)荷裕度以應(yīng)對(duì)系統(tǒng)供需的快速變化,但由于耦合了多種能源,IES在供能側(cè)與用能側(cè)都存在更多的不確定性;同時(shí),負(fù)荷和分布式電源DG(Distributed Generation)出力的波動(dòng)性導(dǎo)致系統(tǒng)在不同時(shí)刻所擁有的負(fù)荷裕度差異巨大,運(yùn)行過(guò)程中面臨的不確定性因素也更多,這為IES 的運(yùn)行帶來(lái)了較大的風(fēng)險(xiǎn)。作為IES 融合的紐帶,儲(chǔ)能及其優(yōu)化配置和運(yùn)行為IES抵御風(fēng)險(xiǎn)、實(shí)現(xiàn)安全經(jīng)濟(jì)運(yùn)行提供了重要保障。

    目前,國(guó)內(nèi)外對(duì)儲(chǔ)能優(yōu)化配置的研究主要存在以下問(wèn)題。①以投資主體的單一利益為目標(biāo)進(jìn)行規(guī)劃,可能忽略了儲(chǔ)能的其他隱性價(jià)值。一方面,部分文獻(xiàn)以某一類型儲(chǔ)能帶來(lái)的社會(huì)效益最佳為目標(biāo)進(jìn)行規(guī)劃:在電力系統(tǒng)范疇,??紤]儲(chǔ)能投資和運(yùn)行成本最低[1-2]、年綜合失負(fù)荷成本最?。?],或儲(chǔ)能對(duì)配電網(wǎng)網(wǎng)損和節(jié)點(diǎn)電壓波動(dòng)的改善程度[4];在IES 范疇,儲(chǔ)能配置拓展為電/熱/冷/氣等多能存儲(chǔ)的規(guī)劃研究,如文獻(xiàn)[5-7]均從實(shí)現(xiàn)系統(tǒng)經(jīng)濟(jì)性最優(yōu)的角度進(jìn)行了IES 儲(chǔ)電/儲(chǔ)熱設(shè)備的容量配置。但上述文獻(xiàn)僅考慮了配置儲(chǔ)能后傳統(tǒng)的系統(tǒng)經(jīng)濟(jì)性,忽略了儲(chǔ)能在平抑負(fù)荷和DG 出力波動(dòng)、均衡系統(tǒng)負(fù)荷裕度方面所帶來(lái)的社會(huì)效益。另一方面,以提升負(fù)荷裕度為研究重點(diǎn)的文獻(xiàn)多采用切負(fù)荷、進(jìn)行無(wú)功補(bǔ)償和優(yōu)化機(jī)組出力等策略[8-10],忽略了儲(chǔ)能在提高IES 負(fù)荷裕度方面的作用。文獻(xiàn)[11]通過(guò)控制儲(chǔ)能進(jìn)行有序充放電來(lái)緩解DG 出力與負(fù)荷需求的錯(cuò)峰問(wèn)題,實(shí)現(xiàn)了配電網(wǎng)負(fù)荷裕度的均衡分布;文獻(xiàn)[12]研究了電動(dòng)汽車充電、儲(chǔ)能和可中斷負(fù)荷調(diào)度對(duì)配電網(wǎng)負(fù)荷峰值和谷值裕度的提升能力,但并沒(méi)有考慮系統(tǒng)全時(shí)段負(fù)荷裕度的均衡性,上述研究主要針對(duì)電力系統(tǒng),不能完全適用于含多能耦合的IES。②忽略了IES 中儲(chǔ)能的規(guī)劃和運(yùn)行問(wèn)題之間的相互影響。文獻(xiàn)[13]將儲(chǔ)能的規(guī)劃與運(yùn)行進(jìn)行協(xié)同考慮,但研究對(duì)象為單一供能網(wǎng)絡(luò),并未涉及IES。文獻(xiàn)[14-15]均考慮了IES 中多元儲(chǔ)能的互補(bǔ)協(xié)調(diào)特性,分別從儲(chǔ)能配置和調(diào)度兩方面進(jìn)行了研究,但在儲(chǔ)能規(guī)劃的過(guò)程中也應(yīng)考慮儲(chǔ)能建成之后的運(yùn)行問(wèn)題,用運(yùn)行來(lái)指導(dǎo)規(guī)劃,這樣雖然增加了優(yōu)化的難度,但能使儲(chǔ)能的規(guī)劃更加合理,同時(shí)使IES 的運(yùn)行更具經(jīng)濟(jì)性。

    綜上,本文將儲(chǔ)能規(guī)劃和運(yùn)行問(wèn)題相結(jié)合,并將其作為改善IES 負(fù)荷裕度的手段,提出一種考慮負(fù)荷裕度的IES 儲(chǔ)能優(yōu)化配置方法。首先,明確包含電/熱/氣3種能源形式的區(qū)域綜合能源系統(tǒng)RIES(Regional Integrated Energy System)結(jié)構(gòu),并建立含儲(chǔ)能的RIES 模型。綜合考慮RIES 能量耦合、儲(chǔ)能運(yùn)行特性、負(fù)荷裕度的時(shí)空均衡分布、DG 出力和負(fù)荷需求的波動(dòng)性,建立了一種考慮負(fù)荷裕度的RIES儲(chǔ)能配置雙層優(yōu)化模型,通過(guò)上/下層模型之間的相互作用與反饋,在最小化社會(huì)綜合成本的同時(shí),也實(shí)現(xiàn)了儲(chǔ)能容量配置及其充、放能計(jì)劃的最優(yōu)化,提高了RIES負(fù)荷裕度均衡度。

    1 含儲(chǔ)能的RIES模型

    1.1 RIES結(jié)構(gòu)

    本文所研究的對(duì)象為RIES,在RIES 中,主要存在以下3 個(gè)利益主體:能量樞紐運(yùn)營(yíng)商EHO(Energy Hub Operator)、終端用戶(部分含光伏)和儲(chǔ)能運(yùn)營(yíng)商ESO(Energy Storage Operator)。

    1)EHO。在RIES的市場(chǎng)機(jī)制下,EHO 作為能源網(wǎng)絡(luò)與用戶之間的媒介,融合了電、熱、氣為用戶提供能源服務(wù),同時(shí)投資建設(shè)熱電聯(lián)供CHP(Combined Heat and Power)系統(tǒng)以豐富能量來(lái)源并降低購(gòu)能成本,以此獲得更高的利潤(rùn)[16]。

    2)終端用戶。在RIES 中,除了傳統(tǒng)用戶,部分用戶擁有光伏發(fā)電的控制權(quán),其可利用光伏發(fā)電以及自主用能管理降低用能成本。

    3)ESO。ESO 同時(shí)經(jīng)營(yíng)著儲(chǔ)電設(shè)備和儲(chǔ)熱設(shè)備,兩者與其他設(shè)備(如CHP 和DG 等)相配合,可以提高RIES 對(duì)新能源的消納能力,發(fā)揮減少能源浪費(fèi)并降低用能成本的作用。

    1.2 RIES建模

    1.2.1 供能網(wǎng)絡(luò)模型

    1)電網(wǎng)模型。

    輻射形配電網(wǎng)的潮流模型可用如下線性化后的DistFlow模型來(lái)表述:

    式(1)、(2)分別為配電網(wǎng)的有功、無(wú)功平衡模型;式(3)為線路(i,j)的壓降模型;式(4)為功率、電壓和電流之間的關(guān)系式;式(5)—(7)分別為電流、電壓和發(fā)電機(jī)出力的上下限約束。

    2)熱網(wǎng)模型。

    本文忽略二次熱網(wǎng)熱能傳輸損耗及延時(shí),將換熱站與二次熱網(wǎng)等效為熱網(wǎng)負(fù)荷節(jié)點(diǎn),根據(jù)質(zhì)調(diào)節(jié)方式,建立計(jì)及供/回水管道的熱網(wǎng)熱力模型如下:

    3)天然氣網(wǎng)絡(luò)模型。

    天然氣網(wǎng)絡(luò)包括氣源中心、天然氣管道、壓縮機(jī)和天然氣負(fù)荷,對(duì)天然氣管道的建模如下:

    式中:fij為管道(i,j)流過(guò)的天然氣流量;ε為管道傳輸效率;sij=sign(pi-pj)為管道(i,j)潮流方向,即天然氣從高氣壓側(cè)流向低氣壓側(cè),其中sign(·)為符號(hào)函數(shù);Mk為連接節(jié)點(diǎn)i和j的管道k的計(jì)算參數(shù);pi和pj分別為節(jié)點(diǎn)i和節(jié)點(diǎn)j的氣壓;T0和p0分別為標(biāo)準(zhǔn)溫度(15 ℃)和標(biāo)準(zhǔn)大氣壓(101.325 kPa);Dk為管道k的內(nèi)直徑;G為天然氣的相對(duì)密度,取0.6;Lk為管道k的長(zhǎng)度;Tka為流過(guò)管道k的天然氣的平均溫度,取環(huán)境溫度25 ℃;Za為天然氣的平均壓縮系數(shù),取1。

    天然氣在傳輸過(guò)程中會(huì)與管道發(fā)生摩擦并出現(xiàn)氣壓降低的現(xiàn)象,為保證配氣網(wǎng)絡(luò)的正常運(yùn)行,在部分管道接入天然氣壓縮機(jī),對(duì)其的建模如下:

    式中:fim為壓縮機(jī)出口節(jié)點(diǎn)m流入網(wǎng)絡(luò)節(jié)點(diǎn)i的流量;fc為壓縮機(jī)消耗的天然氣流量;fnj為網(wǎng)絡(luò)節(jié)點(diǎn)j流入壓縮機(jī)入口節(jié)點(diǎn)n的天然氣流量;Bk為管道k壓縮機(jī)參數(shù);fk為管道k流過(guò)的天然氣流量;α為天然氣的比熱容,取1.3 J/(kg·K);ηk為管道k壓縮效率;pn/pm為壓縮比。

    1.2.2 電/熱耦合環(huán)節(jié)模型

    在本文所研究的RIES 中,電/熱耦合環(huán)節(jié)是通過(guò)CHP 實(shí)現(xiàn)的,可用能量樞紐來(lái)描述其中的能源耦合關(guān)系,其基本結(jié)構(gòu)如圖1所示[17]。

    圖1 能量樞紐結(jié)構(gòu)Fig.1 Structure of energy hub

    本文所采用的CHP 由電力變壓器、微型燃?xì)廨啓C(jī)MT(MicroTurbine)和燃?xì)忮仩tGB(Gas Boiler)共同構(gòu)成。輸入環(huán)節(jié)包括電能和天然氣,前者直接輸入變壓器,而后者同時(shí)輸入到MT和GB;輸出環(huán)節(jié)包含了電能和冷熱能兩部分,前者由變壓器和MT 供給,而后者則由GB 和MT 共同產(chǎn)生。由此可以得到輸入/輸出能源耦合關(guān)系式如下:

    1.2.3 儲(chǔ)能能量狀態(tài)模型

    儲(chǔ)電設(shè)備和儲(chǔ)熱設(shè)備的運(yùn)行特性基本相同,為更好地體現(xiàn)兩者間的協(xié)調(diào)關(guān)系,可建立如下統(tǒng)一的電/熱儲(chǔ)能能量狀態(tài)SOE(State Of Energy)模型,均以電能及電功率來(lái)進(jìn)行表述:

    1.2.4 RIES負(fù)荷裕度模型

    在電力系統(tǒng)中,最大供電能力LSC(Load Supply Capacity)與供電裕度LSM(Load Supply Margin)等指標(biāo)被廣泛應(yīng)用于負(fù)荷裕度水平的評(píng)價(jià)。然而,這2 個(gè)指標(biāo)不能直接適用于含多能耦合的RIES,因?yàn)橹苯訉?duì)RIES 各子網(wǎng)的能源進(jìn)行疊加是不合理的。因此,本文一方面通過(guò)上文建立的RIES 的耦合環(huán)節(jié)模型實(shí)現(xiàn)子網(wǎng)間能源的耦合,另一方面將在下文考慮儲(chǔ)電和儲(chǔ)熱設(shè)備的能量耦合關(guān)系,并將LSC 和LSM 拓展為最大供能能力ESC(Energy Supply Capacity)和供能裕度ESM(Energy Supply Margin),用以評(píng)價(jià)RIES 各子網(wǎng)的負(fù)荷裕度水平。若各子網(wǎng)在研究時(shí)段內(nèi)均能實(shí)現(xiàn)負(fù)荷裕度均衡度最優(yōu),則表示RIES整體運(yùn)行風(fēng)險(xiǎn)得到了降低。

    1)ESC。

    RIES 中,各子網(wǎng)的ESC 是指在滿足供能、用能平衡約束,且所有設(shè)備均不發(fā)生過(guò)負(fù)荷現(xiàn)象時(shí)網(wǎng)絡(luò)所能供應(yīng)負(fù)荷的最大值[18-20]。在滿足文獻(xiàn)[21]所提約束條件下,其計(jì)算方法如下:

    式中:N為所研究的子網(wǎng)負(fù)荷節(jié)點(diǎn)數(shù);fESCi為所研究的子網(wǎng)的負(fù)荷節(jié)點(diǎn)i在最大運(yùn)行方式下的供能負(fù)荷大小,其既可以是電負(fù)荷,也可以是熱負(fù)荷;Li為負(fù)荷節(jié)點(diǎn)i消耗的能量;ηi為負(fù)荷節(jié)點(diǎn)i供能設(shè)備的供能能力系數(shù),即該節(jié)點(diǎn)額定負(fù)荷Li,n與實(shí)際消耗能量Li的比值。

    2)ESM。

    ESM 是指在子網(wǎng)網(wǎng)絡(luò)結(jié)構(gòu)不變、負(fù)荷不變且滿足安全運(yùn)行約束的情況下,網(wǎng)絡(luò)可以額外接入的負(fù)荷量,其計(jì)算方法為:

    式中:fESMi為所研究的子網(wǎng)的負(fù)荷節(jié)點(diǎn)i在最大運(yùn)行方式下的負(fù)荷裕度大小。

    式(22)表明,Li的大小將直接影響RIES 負(fù)荷裕度的大小和分布。負(fù)荷裕度的分布不均可能導(dǎo)致RIES 總體負(fù)荷裕度充足,但部分時(shí)段部分節(jié)點(diǎn)負(fù)荷裕度不足的問(wèn)題。因此,對(duì)負(fù)荷裕度的優(yōu)化目標(biāo)應(yīng)是對(duì)其分布進(jìn)行均衡化處理[11],即通過(guò)控制儲(chǔ)能設(shè)備進(jìn)行有序充、放能,以應(yīng)對(duì)負(fù)荷需求的驟變,同時(shí)調(diào)和DG出力同負(fù)荷需求的錯(cuò)峰問(wèn)題。

    2 考慮負(fù)荷裕度的RIES 儲(chǔ)能雙層優(yōu)化配置模型

    優(yōu)化配置與優(yōu)化運(yùn)行兩者聯(lián)系緊密:一方面,前者是后者的基礎(chǔ);另一方面,后者給前者方案的形成提供了數(shù)據(jù)參考,這符合雙層規(guī)劃的特征。在上層模型中,對(duì)于儲(chǔ)能優(yōu)化配置問(wèn)題,考慮規(guī)劃層面配置方案的社會(huì)綜合成本盡可能最優(yōu);在下層模型中,考慮運(yùn)行層面的負(fù)荷裕度均衡度最優(yōu),并用于指導(dǎo)儲(chǔ)能的充、放能計(jì)劃。為了充分發(fā)揮RIES 不同能源轉(zhuǎn)換、互補(bǔ)等耦合特性的優(yōu)勢(shì),本文建立了考慮負(fù)荷裕度的RIES儲(chǔ)能雙層優(yōu)化配置模型。

    2.1 上層規(guī)劃的數(shù)學(xué)模型

    上層規(guī)劃的數(shù)學(xué)模型以儲(chǔ)能接入容量為決策變量,以RIES 社會(huì)年綜合成本C最小為目標(biāo),綜合考慮RIES 電/氣能耗成本CEG、儲(chǔ)能的投資成本CINV、設(shè)備運(yùn)維成本COM以及碳交易成本CCAB,建立上層規(guī)劃的目標(biāo)函數(shù)如下:

    式中:r為貼現(xiàn)率,取6%;Ness為儲(chǔ)能安裝數(shù)量;cessi為節(jié)點(diǎn)i處儲(chǔ)能的安裝費(fèi)用,與其類型有關(guān);n為儲(chǔ)能使用年限;nessi為節(jié)點(diǎn)i處安裝的儲(chǔ)能數(shù)量。

    3)設(shè)備運(yùn)維成本COM。

    式中:gcoal為配電網(wǎng)供電耗煤量;ucoal為排放因子;ωCAB為碳排放稅。

    上層規(guī)劃模型的約束條件如下。

    1)儲(chǔ)能接入容量與功率約束。

    2)儲(chǔ)能安裝位置數(shù)量約束。

    2.2 下層規(guī)劃的數(shù)學(xué)模型

    在下層規(guī)劃的數(shù)學(xué)模型中,本文以各時(shí)段儲(chǔ)能充、放能功率為決策變量,以RIES 負(fù)荷裕度最大為目標(biāo),同時(shí)考慮其均衡分布,建立如下的負(fù)荷裕度均衡度優(yōu)化模型:

    由于受到潮流和DG 出力波動(dòng)的影響,僅以負(fù)荷裕度最大為目標(biāo)時(shí)會(huì)導(dǎo)致ESM 在時(shí)間和空間上的分布不均衡,為此引入均衡項(xiàng)λt Abal,通過(guò)調(diào)整該項(xiàng)在式(30)中所占的比例,可達(dá)到調(diào)控ESM 在時(shí)間尺度上的均衡性的效果。

    理論上,若fESCi/Li在某時(shí)段偏小,則表明此時(shí)系統(tǒng)運(yùn)行越限將面臨較大的切負(fù)荷風(fēng)險(xiǎn),因此賦給該時(shí)段一個(gè)較大的λt值,增大均衡項(xiàng)的作用效果,并鼓勵(lì)儲(chǔ)能在此時(shí)段釋放能量;相反,若處于谷時(shí)段時(shí)系統(tǒng)越限的概率低,則賦給該時(shí)段一個(gè)較小的λt,降低均衡項(xiàng)的作用效果,并鼓勵(lì)儲(chǔ)能在此時(shí)段儲(chǔ)存能量。綜上所述,λt的取值方法如下:

    下層規(guī)劃模型的約束條件如下。

    1)拓?fù)浼s束。

    電力、天然氣和熱力系統(tǒng)之間的生產(chǎn)、轉(zhuǎn)換及需求環(huán)節(jié)存在的耦合關(guān)系,以及DG 和儲(chǔ)能設(shè)備的接入,使得RIES 的拓?fù)浣Y(jié)構(gòu)和特性與傳統(tǒng)電網(wǎng)有所不同。因此,參考文獻(xiàn)[11],以生成樹原理為基礎(chǔ)形成RIES拓?fù)浼s束如下:

    式中:βl,ij為表示線路l潮流方向的0-1 整數(shù)變量,當(dāng)βl,ij=1 時(shí)線路l的功率從節(jié)點(diǎn)i流向節(jié)點(diǎn)j,反之則從節(jié)點(diǎn)j流向節(jié)點(diǎn)i;A(i,l)為網(wǎng)絡(luò)的無(wú)向關(guān)聯(lián)矩陣的元素,若節(jié)點(diǎn)i與支路l相關(guān)聯(lián),則A(i,l)=1,否則A(i,l)=0;B為RIES的線路總數(shù)。

    2)供能網(wǎng)絡(luò)安全約束。

    同時(shí),還需要滿足式(5)—(7)分別所示的電流、電壓和發(fā)電機(jī)出力約束,以及式(12)所示的節(jié)點(diǎn)溫度約束。

    3)聯(lián)絡(luò)線功率約束。

    4)節(jié)點(diǎn)功率平衡約束。

    RIES 的節(jié)點(diǎn)功率平衡約束包括配電網(wǎng)功率平衡約束式(1)和式(2)、配熱網(wǎng)功率平衡約束式(8)和配氣網(wǎng)功率平衡約束式(15)。

    3 模型簡(jiǎn)化與求解

    3.1 非線性項(xiàng)處理

    由于供能網(wǎng)絡(luò)的非線性和儲(chǔ)能充、放能狀態(tài)的0-1約束,本文所建立的模型是混合整數(shù)非線性規(guī)劃問(wèn)題。對(duì)此,本文參考文獻(xiàn)[15]所提方法對(duì)模型中的非線性項(xiàng)進(jìn)行處理,以降低求解難度。

    3.2 模型求解思路

    經(jīng)過(guò)3.1節(jié)的處理后,上層規(guī)劃模型為一個(gè)僅含線性不等式約束的混合整數(shù)非線性規(guī)劃問(wèn)題,下層屬于混合整數(shù)二階錐規(guī)劃問(wèn)題。求解時(shí)通常采用上下層交替優(yōu)化,上層模型考慮下層模型返回的各時(shí)段RIES 負(fù)荷裕度優(yōu)化結(jié)果,以解決多場(chǎng)景下儲(chǔ)能容量?jī)?yōu)化配置問(wèn)題,下層依據(jù)上層給出的儲(chǔ)能配置容量,優(yōu)化儲(chǔ)能的充、放能計(jì)劃,最大化RIES 負(fù)荷裕度均衡度,以解決每個(gè)場(chǎng)景下RIES 的負(fù)荷裕度優(yōu)化問(wèn)題。本文采用改進(jìn)遺傳算法與Gurobi求解器的混合策略求解上述雙層優(yōu)化模型,具體求解思路如附錄A圖A1所示。

    4 算例分析

    為了驗(yàn)證本文所提出的考慮儲(chǔ)能配置的RIES負(fù)荷裕度雙層優(yōu)化模型的有效性,本節(jié)將基于由IEEE 33節(jié)點(diǎn)系統(tǒng)、比利時(shí)14節(jié)點(diǎn)配氣系統(tǒng)、12節(jié)點(diǎn)配熱系統(tǒng)和4 個(gè)能量樞紐(CHP1—CHP4)所組成的RIES進(jìn)行仿真,系統(tǒng)結(jié)構(gòu)見(jiàn)附錄A圖A2[22]。為探索ESO 儲(chǔ)能設(shè)備容量規(guī)劃的最優(yōu)值,以CHP4的接入點(diǎn)為例進(jìn)行負(fù)荷裕度優(yōu)化。

    4.1 算例參數(shù)

    如1.1 節(jié)所述,RIES 的電/氣/熱耦合環(huán)節(jié)是通過(guò)CHP 實(shí)現(xiàn)的,故為了盈利,EHO 必須充分利用CHP,將低價(jià)購(gòu)入的天然氣生產(chǎn)為電能與熱能高價(jià)售出。本文采用的CHP 的各個(gè)設(shè)備參數(shù)見(jiàn)附錄B表B1。

    配電網(wǎng)、配熱網(wǎng)與配氣網(wǎng)的聯(lián)絡(luò)線或聯(lián)絡(luò)管道的容量決定了EHO 可以向各個(gè)網(wǎng)絡(luò)購(gòu)能或售能的最大功率,其參數(shù)見(jiàn)附錄B表B2。

    通過(guò)對(duì)某地區(qū)負(fù)荷/光伏年時(shí)序數(shù)據(jù)進(jìn)行場(chǎng)景削減,得到春夏秋冬中4 個(gè)典型日用于表征規(guī)劃年。4 個(gè)規(guī)劃典型日的電、熱負(fù)荷以及光伏出力變化見(jiàn)附錄C 圖C1—C3,該系統(tǒng)總電、熱負(fù)荷值分別為5.5 MW+j3.4 Mvar和4 MW。

    在本文所配置的儲(chǔ)能設(shè)備中,儲(chǔ)電設(shè)備接入DG所在電節(jié)點(diǎn),作為DG 的輔助設(shè)備;儲(chǔ)熱設(shè)備接入CHP 設(shè)備所在熱節(jié)點(diǎn),作為熱源的輔助設(shè)備。待配置儲(chǔ)能的技術(shù)和經(jīng)濟(jì)參數(shù)分別見(jiàn)附錄C 表C1 和表C2。

    ESO 向EHO 的購(gòu)能或售能價(jià)格見(jiàn)附錄C 圖C4[22]。上述因素(儲(chǔ)能設(shè)備的技術(shù)/經(jīng)濟(jì)參數(shù)和購(gòu)、售能價(jià)格等)將為本文后續(xù)優(yōu)化方案的選擇提供依據(jù)。

    4.2 場(chǎng)景設(shè)置

    為驗(yàn)證本文所提模型的有效性,設(shè)定以下不同的規(guī)劃場(chǎng)景進(jìn)行分析對(duì)比。

    1)場(chǎng)景1:不考慮儲(chǔ)能配置,電網(wǎng)負(fù)荷由大電網(wǎng)、DG和CHP設(shè)備承擔(dān),熱網(wǎng)負(fù)荷由CHP設(shè)備承擔(dān)。

    2)場(chǎng)景2:在場(chǎng)景1 的基礎(chǔ)上考慮電儲(chǔ)能的配置,不考慮熱儲(chǔ)能,將儲(chǔ)電設(shè)備作為輔助設(shè)備參與負(fù)荷裕度優(yōu)化。

    3)場(chǎng)景3:在場(chǎng)景1 的基礎(chǔ)上考慮熱儲(chǔ)能的配置,不考慮電儲(chǔ)能,將儲(chǔ)熱設(shè)備作為輔助設(shè)備參與負(fù)荷裕度優(yōu)化。

    4)場(chǎng)景4:在前3種場(chǎng)景的基礎(chǔ)上加入電/熱儲(chǔ)能,考慮儲(chǔ)電設(shè)備與儲(chǔ)熱設(shè)備的協(xié)調(diào)運(yùn)行,并采用本文所提模型進(jìn)行負(fù)荷裕度優(yōu)化。

    4.3 仿真結(jié)果及分析

    4.3.1 規(guī)劃方案經(jīng)濟(jì)性分析

    表1 為4 種規(guī)劃場(chǎng)景下所得到的儲(chǔ)能配置及各項(xiàng)成本的比較結(jié)果。在表中所示的儲(chǔ)能配置結(jié)果中,括號(hào)內(nèi)的數(shù)字表示儲(chǔ)電設(shè)備(儲(chǔ)熱設(shè)備)在配電網(wǎng)(配熱網(wǎng))中的安裝位置,括號(hào)外的數(shù)字表示安裝容量或功率,如(e21)6.71表示在配電網(wǎng)節(jié)點(diǎn)21處配置6.71 MW·h 的儲(chǔ)電設(shè)備;(h5)8.54 表示在配熱網(wǎng)節(jié)點(diǎn)5處配置8.54 MW·h的儲(chǔ)熱設(shè)備。

    表1 規(guī)劃方案結(jié)果Table 1 Results of planning schemes

    對(duì)表1 所得的規(guī)劃結(jié)果進(jìn)行分析對(duì)比可得以下結(jié)論。

    1)與場(chǎng)景1 相比,后3 種場(chǎng)景所需運(yùn)維成本更高,原因在于這3 種場(chǎng)景下的儲(chǔ)能裝置有較高的運(yùn)行維護(hù)費(fèi)用。

    2)與場(chǎng)景1 相比,后3 種場(chǎng)景所需能耗費(fèi)用更低,原因在于后3 種場(chǎng)景下通過(guò)更合理的儲(chǔ)能配置及其相應(yīng)的充、放能計(jì)劃,緩解了ESC與負(fù)荷需求的錯(cuò)峰問(wèn)題,減小了RIES 對(duì)上級(jí)電網(wǎng)的依賴程度。另外,場(chǎng)景4 由于考慮了RIES 負(fù)荷裕度均衡分布和儲(chǔ)能的協(xié)同運(yùn)行,突破了單一儲(chǔ)能的運(yùn)行限制,提高了RIES 對(duì)DG 的利用效率和靈活性,進(jìn)一步降低了系統(tǒng)的能耗成本。

    3)場(chǎng)景4的碳交易成本僅為14.8萬(wàn)元,相比前3種場(chǎng)景分別降低了11.4%、3.9%和5.1%,原因在于在儲(chǔ)電設(shè)備和儲(chǔ)熱設(shè)備的協(xié)同作用下,場(chǎng)景4 中DG利用率得到了提高,且向上級(jí)電網(wǎng)的購(gòu)電量也有所降低。

    4)場(chǎng)景4 下的社會(huì)年綜合成本比場(chǎng)景1—3 分別降低了104.3、10.2、116.1萬(wàn)元。由此可見(jiàn),雖然場(chǎng)景4 下的儲(chǔ)能投資成本較高,但在綜合考慮多項(xiàng)成本、效益指標(biāo)后,該方案實(shí)現(xiàn)了社會(huì)綜合成本的最小化。

    在經(jīng)濟(jì)性分析的基礎(chǔ)上,下文將首先分析場(chǎng)景1 下(即未進(jìn)行負(fù)荷裕度優(yōu)化)的RIES 供能能力,然后對(duì)優(yōu)化前后的RIES負(fù)荷裕度進(jìn)行對(duì)比分析。

    4.3.2 RIES供能能力分析

    由于RIES 的負(fù)荷需求隨時(shí)間變化而變化,且各子網(wǎng)變化趨勢(shì)有差異,故本文選取某典型日分別對(duì)場(chǎng)景1 下RIES 的配電網(wǎng)和配熱網(wǎng)進(jìn)行供能能力分析,結(jié)果如圖2所示。

    圖2 配電網(wǎng)、配熱網(wǎng)一天內(nèi)的ESC變化曲線Fig.2 Curve of ESC variation in whole day for distributed power network and distributed heat network

    由圖2 可以看出,無(wú)論是配電網(wǎng)還是配熱網(wǎng),其ESC 的變化都與一天內(nèi)的用戶側(cè)負(fù)荷變化相反:在峰時(shí)段,允許額外接入的負(fù)荷較小,子網(wǎng)的負(fù)荷裕度較??;在平時(shí)段或谷時(shí)段,子網(wǎng)的負(fù)荷裕度較大。

    以上結(jié)論說(shuō)明當(dāng)不考慮負(fù)荷裕度均衡分布時(shí),RIES 的ESC 變化與負(fù)荷需求存在錯(cuò)峰問(wèn)題,在負(fù)荷需求上升時(shí),其ESC 出現(xiàn)下降,而在負(fù)荷需求下降時(shí),其ESC反而上升。因此,有必要通過(guò)合理的儲(chǔ)能配置提高負(fù)荷裕度均衡度。

    4.3.3 RIES負(fù)荷裕度優(yōu)化結(jié)果分析

    本節(jié)選取場(chǎng)景1、2 和場(chǎng)景4 下某典型日對(duì)配電網(wǎng)負(fù)荷裕度優(yōu)化結(jié)果進(jìn)行對(duì)比分析,選取場(chǎng)景1、3和場(chǎng)景4 下某典型日對(duì)配熱網(wǎng)負(fù)荷裕度優(yōu)化結(jié)果進(jìn)行對(duì)比分析。

    1)配電網(wǎng)負(fù)荷裕度優(yōu)化。

    圖3 展示了優(yōu)化前后配電網(wǎng)負(fù)荷裕度的對(duì)比,可以得到以下結(jié)論。

    圖3 配電網(wǎng)負(fù)荷裕度優(yōu)化前后對(duì)比Fig.3 Comparison of load margin for distributed power network between before and after optimization

    1)在優(yōu)化前(即場(chǎng)景1),10:00—18:00時(shí)段系統(tǒng)的負(fù)荷裕度處于較低水平,且16:00—18:00 時(shí)段負(fù)荷裕度小于0,表明RIES已經(jīng)運(yùn)行在不安全狀態(tài)。

    2)根據(jù)場(chǎng)景2、4 的負(fù)荷裕度優(yōu)化結(jié)果可得,在儲(chǔ)電設(shè)備的參與下,21:00—24:00、01:00—08:00時(shí)段負(fù)荷裕度有所削減(儲(chǔ)電設(shè)備充能);09:00—18:00時(shí)段負(fù)荷裕度有所提高(儲(chǔ)電設(shè)備放能),且此時(shí)段恰好是負(fù)荷高峰期,失負(fù)荷風(fēng)險(xiǎn)較大,證明了儲(chǔ)能設(shè)備的參與對(duì)調(diào)整負(fù)荷裕度的有效性。

    3)由于在下層模型的目標(biāo)函數(shù)中考慮了負(fù)荷裕度均衡度,場(chǎng)景4 下RIES 的負(fù)荷裕度變化相較于其他場(chǎng)景更加平緩,其方差由場(chǎng)景1 下的0.034 縮小為0.011,負(fù)荷裕度與負(fù)荷需求錯(cuò)峰問(wèn)題得到緩解。

    4)場(chǎng)景4對(duì)負(fù)荷裕度的均衡度優(yōu)于場(chǎng)景2,這是由于場(chǎng)景2下儲(chǔ)電設(shè)備僅用于緩解DG 出力波動(dòng),故即使在部分時(shí)段(如03:00—06:00時(shí)段)其負(fù)荷裕度值大于場(chǎng)景4,但其總體波動(dòng)性也更大。

    2)配熱網(wǎng)負(fù)荷裕度優(yōu)化。

    圖4 展示了優(yōu)化前后配熱網(wǎng)負(fù)荷裕度的對(duì)比,可以得到以下結(jié)論。

    圖4 配熱網(wǎng)負(fù)荷裕度優(yōu)化前后對(duì)比圖Fig.4 Comparison of load margin for distributed heat network between before and after optimization

    1)對(duì)于場(chǎng)景1,由于沒(méi)有儲(chǔ)能設(shè)備的參與,多個(gè)時(shí)段均出現(xiàn)了負(fù)荷裕度小于0,即配熱網(wǎng)運(yùn)行在不安全狀態(tài)的情況。

    2)對(duì)于場(chǎng)景3,由于01:00—05:00 時(shí)段EHO的售熱價(jià)格非常低,僅為0.3 元/(kW·h)左右,而06:00—12:00 時(shí)段則為0.75 元/(kW·h)左右。因此,場(chǎng)景3 下ESO 的儲(chǔ)熱設(shè)備會(huì)選擇在01:00—05:00 時(shí)段以較大功率進(jìn)行儲(chǔ)熱,并在06:00 時(shí)開始放熱。由于RIES 中CHP 系統(tǒng)容量較小,EHO 不得不從配熱網(wǎng)直接購(gòu)買大量熱能來(lái)滿足ESO的儲(chǔ)熱需要,導(dǎo)致配熱網(wǎng)處于不安全運(yùn)行狀態(tài)。

    3)對(duì)于場(chǎng)景4,隨著負(fù)荷水平的降低,雖然其在部分時(shí)段的負(fù)荷裕度相較于場(chǎng)景1、3 有所下降,但縱觀一天24 h,場(chǎng)景4下均不存在負(fù)荷裕度小于0的情況,即使是在負(fù)荷高峰期,場(chǎng)景4 下RIES 的配熱網(wǎng)也能保有400 kW左右的裕度。

    4.3.4 儲(chǔ)能設(shè)備充、放能分析

    為了分析考慮RIES 負(fù)荷裕度對(duì)儲(chǔ)能充、放能計(jì)劃的調(diào)節(jié)作用,本節(jié)選取場(chǎng)景2和場(chǎng)景4下的ESO儲(chǔ)電設(shè)備優(yōu)化充、放電計(jì)劃,以及場(chǎng)景3 和場(chǎng)景4 下的ESO儲(chǔ)熱設(shè)備優(yōu)化充、放熱計(jì)劃進(jìn)行對(duì)比分析。

    圖5 展示了儲(chǔ)電和儲(chǔ)熱設(shè)備充、放能計(jì)劃優(yōu)化結(jié)果,圖中功率為正表示放能,功率為負(fù)表示儲(chǔ)能,由此得出以下結(jié)論。

    圖5 儲(chǔ)電和儲(chǔ)熱設(shè)備充、放能計(jì)劃優(yōu)化結(jié)果Fig.5 Optimal results of charging and discharging energy plan of electric storage and heat storage equipment

    1)在考慮負(fù)荷裕度均衡后,儲(chǔ)電設(shè)備在負(fù)荷需求較低且負(fù)荷裕度較充足的23:00—24:00、01:00—05:00 時(shí)段進(jìn)行充電,在負(fù)荷需求攀升但DG 出力不足的07:00—10:00、17:00—23:00 時(shí)段進(jìn)行放電,在DG 出力充足的10:00—15:00 時(shí)段進(jìn)行充電??梢?jiàn)充放電控制合理,較好地均衡了負(fù)荷裕度的時(shí)間分布。

    2)在場(chǎng)景2 下,由于ESO 會(huì)以最大化經(jīng)濟(jì)利益為目標(biāo)制定充、放電計(jì)劃,即使配電網(wǎng)的負(fù)荷裕度在17:00—18:00時(shí)段已小于0,即系統(tǒng)運(yùn)行在不安全的狀態(tài)下,ESO 依然會(huì)避免在此時(shí)大幅放電;相比之下,場(chǎng)景4 下的儲(chǔ)電設(shè)備在17:00—18:00 時(shí)段的放電功率有所提高,雖然此時(shí)放電的收益更低,但卻保證了配電網(wǎng)運(yùn)行在安全范圍內(nèi),且獲得了1 000 kW左右的負(fù)荷裕度。

    3)在電、熱負(fù)荷晚高峰時(shí)段(19:00—23:00 時(shí)段),場(chǎng)景4 下的CHP4將增加供熱功率和發(fā)電功率,此時(shí)儲(chǔ)熱設(shè)備蓄熱,儲(chǔ)電設(shè)備放電,兩者協(xié)調(diào)運(yùn)行對(duì)RIES起到了削峰填谷的作用。

    4)在場(chǎng)景4 下,由于考慮了負(fù)荷裕度均衡,ESO放棄了在熱價(jià)最低時(shí)段(01:00—05:00時(shí)段)以最大功率儲(chǔ)熱,而是選擇在時(shí)間更長(zhǎng)時(shí)段以較平緩的功率蓄熱,使得場(chǎng)景3 下凌晨的負(fù)荷裕度驟降的現(xiàn)象消失。

    以上結(jié)果證明由于考慮了負(fù)荷裕度均衡度,本文所建立的雙層優(yōu)化模型對(duì)合理調(diào)整ESO儲(chǔ)電以及儲(chǔ)熱設(shè)備充、放能計(jì)劃具有指導(dǎo)意義。

    5 結(jié)論

    本文首先建立了RIES 模型及其耦合環(huán)節(jié)模型,定義了ESC 與ESM 這2 個(gè)指標(biāo)以刻畫RIES 的負(fù)荷裕度,最后建立了考慮負(fù)荷裕度的RIES 儲(chǔ)能雙層優(yōu)化配置模型。為了驗(yàn)證所提模型的有效性,基于典型RIES進(jìn)行了仿真,并得出以下結(jié)論。

    1)本文的雙層優(yōu)化模型將規(guī)劃與運(yùn)行問(wèn)題相結(jié)合,在優(yōu)化儲(chǔ)能配置的同時(shí),也優(yōu)化了其充、放能計(jì)劃,并實(shí)現(xiàn)了RIES 負(fù)荷裕度均衡化。與其他規(guī)劃方案相比,在本文配置方案下儲(chǔ)能協(xié)同配置減少了單類型儲(chǔ)能的配置需求,有效降低了社會(huì)綜合成本,使RIES的運(yùn)行具有更高的經(jīng)濟(jì)效益和環(huán)境效益。

    2)電儲(chǔ)能和熱儲(chǔ)能協(xié)調(diào)運(yùn)行能有效削峰填谷和緩解DG出力波動(dòng)性,提高了DG消納能力和RIES負(fù)荷裕度均衡度,并實(shí)現(xiàn)了多能存儲(chǔ)方面的優(yōu)勢(shì)互補(bǔ)。

    3)本文在優(yōu)化配置儲(chǔ)能容量的基礎(chǔ)上,通過(guò)考慮負(fù)荷裕度均衡性,進(jìn)一步優(yōu)化了儲(chǔ)能設(shè)備的充、放能計(jì)劃,由于各類儲(chǔ)能充、放能時(shí)段不同,使RIES 在各時(shí)段均具有一定能量?jī)?chǔ)備,提高了系統(tǒng)可靠性。

    附錄見(jiàn)本刊網(wǎng)絡(luò)版(http://www.epae.cn)。

    猜你喜歡
    裕度熱網(wǎng)時(shí)段
    熱網(wǎng)異常工況的辨識(shí)
    煤氣與熱力(2022年4期)2022-05-23 12:44:44
    新建熱網(wǎng)與現(xiàn)狀熱網(wǎng)并網(wǎng)升溫方案
    基于動(dòng)態(tài)三維交互的二級(jí)熱網(wǎng)仿真系統(tǒng)
    煤氣與熱力(2021年9期)2021-11-06 05:22:46
    關(guān)于熱網(wǎng)換熱站節(jié)能監(jiān)管與應(yīng)用探究
    四個(gè)養(yǎng)生黃金時(shí)段,你抓住了嗎
    基于DFIG可用無(wú)功裕度的風(fēng)電場(chǎng)無(wú)功電壓控制方法
    三環(huán)路核電廠的抗震裕度評(píng)價(jià)
    基于ANN模型的在線電壓穩(wěn)定裕度評(píng)估
    傍晚是交通事故高發(fā)時(shí)段
    分時(shí)段預(yù)約在PICC門診維護(hù)中的應(yīng)用與探討
    国产麻豆69| 国产成人一区二区在线| 免费看不卡的av| 伦理电影大哥的女人| 纵有疾风起免费观看全集完整版| 国产麻豆69| www.熟女人妻精品国产 | 免费少妇av软件| 高清视频免费观看一区二区| 边亲边吃奶的免费视频| 18禁裸乳无遮挡动漫免费视频| 人人妻人人添人人爽欧美一区卜| 在线免费观看不下载黄p国产| 女人被躁到高潮嗷嗷叫费观| 在线精品无人区一区二区三| 亚洲精品久久久久久婷婷小说| 一级毛片 在线播放| 黄色怎么调成土黄色| 热99国产精品久久久久久7| 久久久精品区二区三区| 少妇的逼水好多| 久久久国产精品麻豆| 街头女战士在线观看网站| 久久人人爽人人爽人人片va| 校园人妻丝袜中文字幕| 激情视频va一区二区三区| 午夜福利视频在线观看免费| 人体艺术视频欧美日本| 亚洲美女视频黄频| 黄色毛片三级朝国网站| 国产爽快片一区二区三区| 国产精品人妻久久久久久| av.在线天堂| 少妇被粗大猛烈的视频| 国产精品女同一区二区软件| 老司机亚洲免费影院| 成年av动漫网址| 亚洲,欧美精品.| 午夜福利在线观看免费完整高清在| 观看av在线不卡| 男女高潮啪啪啪动态图| 免费黄色在线免费观看| 久热这里只有精品99| 在线天堂最新版资源| 最黄视频免费看| 国产视频首页在线观看| 久久精品夜色国产| 国产淫语在线视频| 亚洲精品aⅴ在线观看| 国产一区二区在线观看av| 日本av手机在线免费观看| 九色亚洲精品在线播放| 久久这里有精品视频免费| 插逼视频在线观看| 国产欧美另类精品又又久久亚洲欧美| av视频免费观看在线观看| 久久国产精品大桥未久av| 亚洲色图 男人天堂 中文字幕 | 亚洲av在线观看美女高潮| 国产在线视频一区二区| 草草在线视频免费看| 久久青草综合色| 性色av一级| 亚洲欧美成人精品一区二区| 看非洲黑人一级黄片| 成人影院久久| 中文字幕精品免费在线观看视频 | 多毛熟女@视频| 亚洲精品日本国产第一区| 亚洲欧美色中文字幕在线| 中国美白少妇内射xxxbb| 成年人免费黄色播放视频| 人人妻人人爽人人添夜夜欢视频| 亚洲第一av免费看| 欧美97在线视频| 伦精品一区二区三区| 成人亚洲欧美一区二区av| 亚洲av中文av极速乱| 日本wwww免费看| 夫妻午夜视频| 亚洲精品av麻豆狂野| 欧美日韩国产mv在线观看视频| 在线免费观看不下载黄p国产| 高清在线视频一区二区三区| 97超碰精品成人国产| 国产在线免费精品| 国产免费一级a男人的天堂| 大片电影免费在线观看免费| 日韩 亚洲 欧美在线| 亚洲三级黄色毛片| 亚洲精品,欧美精品| 伦理电影免费视频| 久久精品久久久久久久性| 97超碰精品成人国产| 母亲3免费完整高清在线观看 | 国产精品国产三级国产av玫瑰| 精品少妇内射三级| 精品视频人人做人人爽| 只有这里有精品99| 久久人人爽av亚洲精品天堂| 性色av一级| 捣出白浆h1v1| 一区二区三区精品91| 观看美女的网站| 国产成人欧美| 捣出白浆h1v1| 日韩精品有码人妻一区| 高清欧美精品videossex| 美女中出高潮动态图| av网站免费在线观看视频| 日本猛色少妇xxxxx猛交久久| 下体分泌物呈黄色| 亚洲久久久国产精品| 日本色播在线视频| 欧美日韩国产mv在线观看视频| 亚洲人与动物交配视频| 免费看不卡的av| 丝袜人妻中文字幕| 久久国产亚洲av麻豆专区| 日韩制服丝袜自拍偷拍| 男的添女的下面高潮视频| 午夜久久久在线观看| 最近的中文字幕免费完整| 热99久久久久精品小说推荐| 亚洲av成人精品一二三区| 亚洲婷婷狠狠爱综合网| 午夜福利乱码中文字幕| 午夜av观看不卡| 国产精品久久久久久av不卡| 一级爰片在线观看| 丝袜人妻中文字幕| 亚洲婷婷狠狠爱综合网| 另类亚洲欧美激情| 免费播放大片免费观看视频在线观看| 亚洲av免费高清在线观看| 亚洲精品久久成人aⅴ小说| 亚洲精品自拍成人| 大片免费播放器 马上看| 国产乱人偷精品视频| 日本免费在线观看一区| 国产高清国产精品国产三级| 熟妇人妻不卡中文字幕| 交换朋友夫妻互换小说| 啦啦啦视频在线资源免费观看| 久久精品国产a三级三级三级| 中文字幕另类日韩欧美亚洲嫩草| 伦理电影免费视频| 99久久人妻综合| 日本免费在线观看一区| 久久综合国产亚洲精品| 久久久国产欧美日韩av| 国产黄频视频在线观看| 丰满少妇做爰视频| 亚洲欧美成人精品一区二区| 精品久久久久久电影网| 国产 精品1| 日韩成人伦理影院| 一级毛片 在线播放| 黄色视频在线播放观看不卡| 欧美激情极品国产一区二区三区 | 亚洲精品久久成人aⅴ小说| 婷婷色麻豆天堂久久| 91在线精品国自产拍蜜月| 国产在线一区二区三区精| 国产成人精品久久久久久| 搡老乐熟女国产| 国产精品久久久久久久久免| 日本免费在线观看一区| 校园人妻丝袜中文字幕| 日本91视频免费播放| a级毛片在线看网站| av国产精品久久久久影院| 十八禁高潮呻吟视频| 大话2 男鬼变身卡| 校园人妻丝袜中文字幕| 国产爽快片一区二区三区| 另类亚洲欧美激情| 久久久久久久精品精品| 亚洲欧洲日产国产| 欧美少妇被猛烈插入视频| av天堂久久9| 欧美日韩国产mv在线观看视频| 欧美另类一区| 色婷婷av一区二区三区视频| 99热网站在线观看| 午夜日本视频在线| 日韩,欧美,国产一区二区三区| 国产成人午夜福利电影在线观看| 黄色配什么色好看| 亚洲第一区二区三区不卡| 日本色播在线视频| 免费大片18禁| 99热这里只有是精品在线观看| 80岁老熟妇乱子伦牲交| 亚洲欧美精品自产自拍| 菩萨蛮人人尽说江南好唐韦庄| av在线播放精品| 免费看av在线观看网站| a级片在线免费高清观看视频| 成人毛片a级毛片在线播放| 五月玫瑰六月丁香| 欧美日韩国产mv在线观看视频| 亚洲伊人色综图| 女人久久www免费人成看片| 久久这里只有精品19| 成年人午夜在线观看视频| 午夜福利,免费看| 国产麻豆69| 亚洲人与动物交配视频| 肉色欧美久久久久久久蜜桃| 久久影院123| 国产亚洲精品第一综合不卡 | 日本欧美国产在线视频| videosex国产| 欧美少妇被猛烈插入视频| 成年av动漫网址| 一级片免费观看大全| 亚洲精品av麻豆狂野| 成年动漫av网址| 婷婷色综合大香蕉| 人妻系列 视频| 国产精品久久久久久久电影| 国产日韩一区二区三区精品不卡| 精品亚洲成国产av| 1024视频免费在线观看| 国产精品熟女久久久久浪| 熟女av电影| 一本色道久久久久久精品综合| 热re99久久精品国产66热6| 亚洲第一区二区三区不卡| 亚洲精品久久午夜乱码| 青春草视频在线免费观看| 欧美人与性动交α欧美精品济南到 | 又大又黄又爽视频免费| 久久人人爽人人片av| 国产av精品麻豆| 亚洲av电影在线观看一区二区三区| 亚洲少妇的诱惑av| 波多野结衣一区麻豆| 2021少妇久久久久久久久久久| 欧美激情国产日韩精品一区| 在线观看三级黄色| av线在线观看网站| 婷婷色综合大香蕉| 在线观看www视频免费| 国产精品久久久久久久电影| 亚洲高清免费不卡视频| 不卡视频在线观看欧美| 国产视频首页在线观看| 亚洲精品久久成人aⅴ小说| 亚洲国产av新网站| 午夜精品国产一区二区电影| 欧美精品一区二区免费开放| 久久人人爽人人片av| 全区人妻精品视频| 丰满迷人的少妇在线观看| 免费人妻精品一区二区三区视频| 国产精品久久久久久av不卡| 欧美+日韩+精品| 又大又黄又爽视频免费| 久久久久国产精品人妻一区二区| 中文字幕另类日韩欧美亚洲嫩草| 久久午夜福利片| 亚洲欧美成人精品一区二区| 97在线视频观看| 十八禁高潮呻吟视频| 国产精品一国产av| 97精品久久久久久久久久精品| 全区人妻精品视频| 国产成人a∨麻豆精品| 亚洲精华国产精华液的使用体验| 少妇人妻 视频| 熟女人妻精品中文字幕| 国产成人精品久久久久久| 午夜福利视频精品| 不卡视频在线观看欧美| 国产成人av激情在线播放| 亚洲国产看品久久| 日韩大片免费观看网站| 亚洲精品国产av蜜桃| 制服丝袜香蕉在线| 高清欧美精品videossex| 免费看光身美女| 免费久久久久久久精品成人欧美视频 | 热99久久久久精品小说推荐| 男女边吃奶边做爰视频| 免费黄频网站在线观看国产| 亚洲国产精品专区欧美| 乱人伦中国视频| 一区在线观看完整版| 国产成人aa在线观看| 美女脱内裤让男人舔精品视频| 日本爱情动作片www.在线观看| 80岁老熟妇乱子伦牲交| 色哟哟·www| 男女午夜视频在线观看 | tube8黄色片| 五月伊人婷婷丁香| 热re99久久精品国产66热6| 天天躁夜夜躁狠狠久久av| 咕卡用的链子| 国产成人精品福利久久| 啦啦啦视频在线资源免费观看| 亚洲国产精品999| 亚洲欧洲日产国产| 秋霞伦理黄片| 爱豆传媒免费全集在线观看| 日日啪夜夜爽| 中文字幕人妻熟女乱码| 国产精品久久久久久久久免| 成人黄色视频免费在线看| 最近2019中文字幕mv第一页| 国产免费一级a男人的天堂| 国语对白做爰xxxⅹ性视频网站| 日本黄色日本黄色录像| 午夜免费鲁丝| 国产精品国产av在线观看| 少妇的丰满在线观看| 亚洲美女搞黄在线观看| 国产麻豆69| 国产欧美日韩一区二区三区在线| 婷婷色综合大香蕉| 久久久久久人人人人人| av福利片在线| 亚洲成人一二三区av| 国产精品国产三级专区第一集| 亚洲欧美清纯卡通| 亚洲人成77777在线视频| 久久热在线av| 午夜免费男女啪啪视频观看| 日本猛色少妇xxxxx猛交久久| 日本欧美视频一区| 午夜激情久久久久久久| 日韩电影二区| 欧美精品国产亚洲| 久久精品熟女亚洲av麻豆精品| 日韩中文字幕视频在线看片| 老女人水多毛片| 日本wwww免费看| 天天躁夜夜躁狠狠躁躁| 日本色播在线视频| 国产欧美日韩一区二区三区在线| 免费播放大片免费观看视频在线观看| 丝瓜视频免费看黄片| 日本与韩国留学比较| 一本久久精品| 亚洲丝袜综合中文字幕| 国产欧美另类精品又又久久亚洲欧美| 国产乱人偷精品视频| 欧美日韩精品成人综合77777| 建设人人有责人人尽责人人享有的| 女人久久www免费人成看片| 亚洲国产欧美在线一区| 一个人免费看片子| 亚洲国产精品专区欧美| 国产日韩欧美视频二区| 寂寞人妻少妇视频99o| 人妻人人澡人人爽人人| 中文字幕另类日韩欧美亚洲嫩草| 日本欧美国产在线视频| 交换朋友夫妻互换小说| 久久精品夜色国产| 超色免费av| 久久久久精品性色| 亚洲精品视频女| 97精品久久久久久久久久精品| 精品人妻偷拍中文字幕| 99热全是精品| 久久99精品国语久久久| 久久99热6这里只有精品| 欧美日韩成人在线一区二区| 亚洲精品中文字幕在线视频| av免费观看日本| 在线观看www视频免费| 大香蕉久久成人网| 另类精品久久| 一区在线观看完整版| 黄色 视频免费看| 中文欧美无线码| 激情视频va一区二区三区| 国产黄频视频在线观看| 欧美成人精品欧美一级黄| 观看av在线不卡| 在线亚洲精品国产二区图片欧美| 一区二区av电影网| 97在线人人人人妻| 精品午夜福利在线看| 日本欧美视频一区| 97在线视频观看| 热99久久久久精品小说推荐| 国产成人av激情在线播放| 亚洲在久久综合| 视频中文字幕在线观看| 女性生殖器流出的白浆| 国语对白做爰xxxⅹ性视频网站| 国产有黄有色有爽视频| 欧美国产精品va在线观看不卡| 老司机影院毛片| 亚洲精品国产色婷婷电影| 中文字幕免费在线视频6| 免费av不卡在线播放| av免费在线看不卡| 久久青草综合色| 少妇 在线观看| 久久精品国产鲁丝片午夜精品| 久久av网站| 亚洲一码二码三码区别大吗| 亚洲熟女精品中文字幕| 91久久精品国产一区二区三区| 国产1区2区3区精品| 精品一区二区三区视频在线| 国产精品久久久av美女十八| 大话2 男鬼变身卡| 久久人人爽av亚洲精品天堂| 精品视频人人做人人爽| 精品一区二区三区四区五区乱码 | 一区二区三区乱码不卡18| 久久久久久久久久成人| 九色成人免费人妻av| 国产深夜福利视频在线观看| 国产黄频视频在线观看| 精品一区二区三区视频在线| 国产高清不卡午夜福利| 亚洲情色 制服丝袜| 成人综合一区亚洲| 精品一区二区三区四区五区乱码 | 熟女电影av网| 一二三四中文在线观看免费高清| 人妻 亚洲 视频| 久久久久久伊人网av| 三上悠亚av全集在线观看| 精品福利永久在线观看| 在线观看免费高清a一片| 久久久精品免费免费高清| 一本—道久久a久久精品蜜桃钙片| 亚洲国产成人一精品久久久| 91午夜精品亚洲一区二区三区| 久久国产亚洲av麻豆专区| 丝袜喷水一区| 成人黄色视频免费在线看| 欧美日韩av久久| 精品熟女少妇av免费看| 一二三四中文在线观看免费高清| 精品一品国产午夜福利视频| 中文字幕人妻熟女乱码| 少妇人妻精品综合一区二区| 五月开心婷婷网| 两个人看的免费小视频| 久久久欧美国产精品| 大码成人一级视频| 麻豆精品久久久久久蜜桃| 亚洲av日韩在线播放| 国语对白做爰xxxⅹ性视频网站| 亚洲国产日韩一区二区| 好男人视频免费观看在线| 免费观看av网站的网址| 丝袜美足系列| 黄色一级大片看看| 青青草视频在线视频观看| 国产日韩一区二区三区精品不卡| av线在线观看网站| 91aial.com中文字幕在线观看| 久久久久网色| 美女脱内裤让男人舔精品视频| 久久青草综合色| 蜜桃在线观看..| √禁漫天堂资源中文www| 日本91视频免费播放| 搡女人真爽免费视频火全软件| 国产无遮挡羞羞视频在线观看| 国产永久视频网站| 亚洲国产最新在线播放| 人妻一区二区av| 国产精品免费大片| 国产高清三级在线| 草草在线视频免费看| 18禁动态无遮挡网站| 大香蕉久久成人网| 秋霞在线观看毛片| 中文字幕另类日韩欧美亚洲嫩草| 一区二区三区四区激情视频| 国产毛片在线视频| 国产亚洲午夜精品一区二区久久| 国产精品麻豆人妻色哟哟久久| 国产麻豆69| 久久久久久久国产电影| 哪个播放器可以免费观看大片| 熟女人妻精品中文字幕| 成人毛片60女人毛片免费| 美女内射精品一级片tv| 国产在线免费精品| 国产av国产精品国产| 亚洲国产精品999| 欧美精品av麻豆av| 成人午夜精彩视频在线观看| 中文精品一卡2卡3卡4更新| 国产精品久久久av美女十八| 少妇人妻 视频| 久久久久久久大尺度免费视频| 日产精品乱码卡一卡2卡三| 亚洲欧美一区二区三区国产| 欧美精品人与动牲交sv欧美| 久久国内精品自在自线图片| 中文字幕另类日韩欧美亚洲嫩草| 精品卡一卡二卡四卡免费| 黄色一级大片看看| xxxhd国产人妻xxx| 久久久久久久久久成人| 成人18禁高潮啪啪吃奶动态图| 国产成人午夜福利电影在线观看| 国产片内射在线| 国产又色又爽无遮挡免| 国产片内射在线| 另类亚洲欧美激情| 激情五月婷婷亚洲| 成年动漫av网址| 又粗又硬又长又爽又黄的视频| 午夜影院在线不卡| 啦啦啦在线观看免费高清www| 国产精品麻豆人妻色哟哟久久| 97在线视频观看| 伦理电影大哥的女人| 老司机影院成人| 成人黄色视频免费在线看| 91成人精品电影| 狂野欧美激情性xxxx在线观看| 日本黄色日本黄色录像| 亚洲精品日本国产第一区| 亚洲天堂av无毛| 欧美精品高潮呻吟av久久| 欧美性感艳星| 最近中文字幕2019免费版| 久久久久国产精品人妻一区二区| 久久这里只有精品19| 日韩欧美一区视频在线观看| 久久国产亚洲av麻豆专区| 边亲边吃奶的免费视频| 狠狠精品人妻久久久久久综合| 免费久久久久久久精品成人欧美视频 | 永久网站在线| 欧美亚洲日本最大视频资源| 日韩中文字幕视频在线看片| 不卡视频在线观看欧美| 色婷婷久久久亚洲欧美| 欧美日韩精品成人综合77777| 国产又爽黄色视频| 国产av国产精品国产| 国产福利在线免费观看视频| 久久久久久人人人人人| 国产亚洲午夜精品一区二区久久| 国产精品.久久久| 婷婷色综合www| 天天操日日干夜夜撸| 国产精品久久久久久av不卡| 三级国产精品片| 亚洲av电影在线进入| 夜夜骑夜夜射夜夜干| 91午夜精品亚洲一区二区三区| 啦啦啦啦在线视频资源| 国产国语露脸激情在线看| videos熟女内射| 国产男女内射视频| 国产精品人妻久久久久久| 九色成人免费人妻av| 日韩一区二区视频免费看| 搡老乐熟女国产| 免费大片黄手机在线观看| 最近中文字幕高清免费大全6| videosex国产| 亚洲国产最新在线播放| 亚洲精品aⅴ在线观看| 久久久久国产精品人妻一区二区| 亚洲精品日韩在线中文字幕| 欧美精品亚洲一区二区| 国产福利在线免费观看视频| 亚洲人成网站在线观看播放| 亚洲国产精品一区三区| 免费看光身美女| 18禁在线无遮挡免费观看视频| 最近中文字幕2019免费版| 国产精品熟女久久久久浪| 男女边摸边吃奶| 丰满饥渴人妻一区二区三| 成人黄色视频免费在线看| 亚洲美女搞黄在线观看| 黑人高潮一二区| 日韩欧美一区视频在线观看| 视频中文字幕在线观看| 国产白丝娇喘喷水9色精品| 亚洲国产最新在线播放| 国产免费一级a男人的天堂| 一级黄片播放器| 亚洲天堂av无毛| av在线app专区| 久久久国产一区二区| 尾随美女入室| 美女福利国产在线| 永久网站在线| 97精品久久久久久久久久精品| 成人毛片60女人毛片免费| 免费播放大片免费观看视频在线观看| 韩国av在线不卡| 69精品国产乱码久久久| 国产日韩欧美亚洲二区| 啦啦啦在线观看免费高清www| 国产成人一区二区在线| 最近的中文字幕免费完整| 久久99蜜桃精品久久| 亚洲精品成人av观看孕妇| 国产女主播在线喷水免费视频网站| 丰满饥渴人妻一区二区三| 大香蕉久久成人网| 丝瓜视频免费看黄片| 一本色道久久久久久精品综合| 成年女人在线观看亚洲视频|