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

    基于動(dòng)態(tài)碳排放因子的電力低碳優(yōu)化調(diào)度策略

    2024-01-02 00:00:00劉洋廣
    電器與能效管理技術(shù) 2024年12期
    關(guān)鍵詞:儲(chǔ)能

    摘 要: 基于利用動(dòng)態(tài)碳排放因子變化規(guī)律引導(dǎo)用戶改變用能曲線,實(shí)現(xiàn)低碳減排,提出一種基于動(dòng)態(tài)碳排放因子的電力低碳優(yōu)化調(diào)度策略。以動(dòng)態(tài)碳排放因子為引導(dǎo)信號,通過儲(chǔ)能系統(tǒng)和可調(diào)度負(fù)荷的電力低碳特性引導(dǎo)用戶側(cè)進(jìn)行低碳激勵(lì)型響應(yīng),改變儲(chǔ)能系統(tǒng)的充放電和可調(diào)度負(fù)荷的運(yùn)行狀態(tài);有效利用動(dòng)態(tài)碳排放因子變化規(guī)律,將實(shí)時(shí)的碳排放量守恒原則擴(kuò)展到調(diào)度周期內(nèi),有序引導(dǎo)用戶側(cè)改變用能曲線。最后,通過算例分析對比驗(yàn)證所提策略的有效性和先進(jìn)性。

    關(guān)鍵詞: 碳排放因子; 電力低碳; 儲(chǔ)能; 可調(diào)度負(fù)荷; 碳勢

    中圖分類號: TM73

    文獻(xiàn)標(biāo)志碼: A

    文章編號: 2095-8188(2024)12-0026-07

    DOI: 10.16628/j.cnki.2095-8188.2024.12.005

    Low Carbon Optimal Dispatch Strategy for Electricity Based on Dynamic Carbon Emission Factors

    LIU Yangguang

    (Dating Renewable Energy Test and Research Institute Co.,Ltd., Beijing 100052, China)

    Abstract: A low carbon optimisation scheduling strategy based on the dynamic carbon emission factor is proposed to guide the user to change the energy consumption curves to achieve low carbon emission reduction by utilizing the changing laws of dynamic carbonemission factors.Taking the dynamic carbon emission factor as the guiding signal,the low carbon characteristics of the energy storage system and dispatchable load are used to guide the user side to make low carbon incentive responses and change the charging and discharging of the energy storage system and the operation state of the dispatchable load.The dynamic carbon emission factor changing rule is effectively utilized to extend the real-time carbon emission conservation principle to the scheduling cycle,so as to guide the user side to change the energy consumption curve in an orderly manner.Finally,the effectiveness and advancement of the proposed strategy are verified by the comparison with example analysis.

    Key words: carbon emission factor; electricity low carbon; energy storage; dispatchable load; carbon potential

    0 引 言

    當(dāng)前,隨著我國“碳達(dá)峰、碳中和”目標(biāo)[1的提出和各項(xiàng)相關(guān)工作的快速推進(jìn),國內(nèi)開始鼓勵(lì)大型企業(yè)積極發(fā)揮示范帶動(dòng)作用,逐步提升綠色能源的消費(fèi)比例[2。各地通過加強(qiáng)耗電量較高的企業(yè)使用綠色能源的硬性約束條件,推動(dòng)電力能源公司定期整理、發(fā)布本地綠色能源時(shí)段分布狀況,從而有序引導(dǎo)用電用戶積極消費(fèi)綠色能源[3。同時(shí),在電力供應(yīng)能力允許的時(shí)間段,在進(jìn)行需求側(cè)管理時(shí),優(yōu)先保障綠色能源消費(fèi)比例較高的用戶。

    將非零碳勢發(fā)電機(jī)組(物理碳排放不為0)產(chǎn)生的二氧化碳(CO2)按照電力碳排放分配規(guī)則分?jǐn)偨o用戶,是明確用戶碳排放量的重要手段[4-6。綠色電力是由零碳勢發(fā)電機(jī)組(物理碳排放為0)產(chǎn)生的,節(jié)點(diǎn)碳勢隨系統(tǒng)中各發(fā)電機(jī)組流入該節(jié)點(diǎn)的功率不同而變化,節(jié)點(diǎn)碳勢越低,表明該節(jié)點(diǎn)流入的低碳勢發(fā)電機(jī)組功率越多。有效利用動(dòng)態(tài)碳排放因子[7-8變化規(guī)律,有序引導(dǎo)用戶改變用能曲線,不但能有效降低用戶碳排放量,還能提高綠色電力消費(fèi)占比,進(jìn)而從整體上提高電力系統(tǒng)新能源消納率。

    文獻(xiàn)[9]通過引入碳排放流(CEF)理論,建立以最低碳排放分?jǐn)偭繛槟繕?biāo)函數(shù)的電力系統(tǒng)低碳優(yōu)化運(yùn)行數(shù)學(xué)模型,但僅考慮可轉(zhuǎn)移負(fù)荷在電力系統(tǒng)低碳優(yōu)化中的作用。文獻(xiàn)[10]進(jìn)一步給出電力系統(tǒng)計(jì)及碳排放流和用戶側(cè)需求響應(yīng)的雙層規(guī)劃模型,結(jié)果表明通過用戶側(cè)需求響應(yīng)對負(fù)荷曲線進(jìn)行調(diào)節(jié),可以提高綠色能源的利用率,但未探究不同類型的用戶側(cè)需求響應(yīng)負(fù)荷對低碳優(yōu)化的能力。文獻(xiàn)[11]提出一種含電轉(zhuǎn)氣、燃?xì)廨啓C(jī)、碳捕集電廠、風(fēng)電場聯(lián)合運(yùn)行的低碳經(jīng)濟(jì)調(diào)度模型,并驗(yàn)證所提模型能夠有效提高風(fēng)電消納能力。文獻(xiàn)[12]提出一種考慮電轉(zhuǎn)氣熱回收與蓄熱罐協(xié)調(diào)運(yùn)行的風(fēng)電消納低碳經(jīng)濟(jì)調(diào)度模型,有效提升風(fēng)電的消納能力,降低系統(tǒng)的總成本和碳排放量。文獻(xiàn)[13]將碳排放流理論和夏普利值法引入碳交易模型中,通過價(jià)格型的需求響應(yīng)對用戶進(jìn)行碳排放控制,從而實(shí)現(xiàn)用戶對高碳排放強(qiáng)度能源的消耗。文獻(xiàn)[14]強(qiáng)調(diào)儲(chǔ)能系統(tǒng)通過優(yōu)化調(diào)節(jié)充放電特性,能夠在一定程度上影響所在區(qū)域的對外碳強(qiáng)度特性,從而有效改善新能源和負(fù)荷的功率不匹配情況。上述研究表明,從“碳視角”出發(fā)的電力系統(tǒng)能量調(diào)度為實(shí)現(xiàn)電力低碳提供了行之有效的路徑。但是上述研究均未考慮到電力系統(tǒng)中可能存在的多種類型可調(diào)度負(fù)荷在電力低碳優(yōu)化中的作用。研究可調(diào)度負(fù)荷的電力低碳特性,對促進(jìn)電力系統(tǒng)低碳優(yōu)化具有重要的意義。

    以動(dòng)態(tài)碳排放因子為引導(dǎo)信號,優(yōu)化調(diào)度儲(chǔ)能系統(tǒng)和可調(diào)度負(fù)荷,對電力低碳發(fā)展具有重要的現(xiàn)實(shí)意義,本文提出一種基于動(dòng)態(tài)碳排放因子的電力低碳優(yōu)化調(diào)度策略。

    1 儲(chǔ)能碳排放特性

    碳排放流是表征系統(tǒng)潮流所攜帶的碳排放量的一種虛擬網(wǎng)絡(luò)流,相當(dāng)于給流過系統(tǒng)每個(gè)節(jié)點(diǎn)的潮流添加一個(gè)碳排放量的標(biāo)簽[15-16。碳排放從物理碳排放值不為0的發(fā)電機(jī)組出發(fā),跟隨潮流在網(wǎng)絡(luò)中流動(dòng),最終由用戶消費(fèi)。

    為了更清晰地描述電力系統(tǒng)網(wǎng)絡(luò)節(jié)點(diǎn)的碳排放特性,碳排放流理論中提出采用節(jié)點(diǎn)碳勢[17的概念進(jìn)行碳計(jì)量,其物理意義表述為在該節(jié)點(diǎn)消費(fèi)單位電量所造成的等效于發(fā)電側(cè)的碳排放值。當(dāng)考慮網(wǎng)損時(shí),電力網(wǎng)絡(luò)中各節(jié)點(diǎn)碳勢表示為

    EnN=∑Nn=1(Pi→nN+λPi→n,lossN)Ei→nN+∑Kk=1Pn,kGEn,kG∑Ni=1Pi→nN+∑Kk=1Pn,kG(1)

    式中: EnN——節(jié)點(diǎn)n的節(jié)點(diǎn)碳勢;

    N——電力網(wǎng)絡(luò)節(jié)點(diǎn)個(gè)數(shù);

    K——電力網(wǎng)絡(luò)中存在機(jī)組注入的節(jié)點(diǎn)個(gè)數(shù);

    Pi→nN——支路i→n流過的正向有功功率;

    Pi→n,lossN——支路i→n的網(wǎng)損;

    Ei→nN——支路i→n的碳流密度;

    Pn,kG——發(fā)電機(jī)組k在節(jié)點(diǎn)n注入的有功功率;

    En,kG——發(fā)電機(jī)組k的碳流密度;

    λ——網(wǎng)損分?jǐn)傁禂?shù),λ∈[0,1]。

    根據(jù)碳排放流特性,任意支路碳流密度都等于該支路始端節(jié)點(diǎn)碳勢,即Ei→nN=EiN。

    2 可調(diào)度負(fù)荷的電力低碳特性研究

    2.1 儲(chǔ)能低碳特性

    儲(chǔ)能系統(tǒng)在節(jié)點(diǎn)碳勢的變化中扮演著重要的角色。利用動(dòng)態(tài)碳排放因子變化規(guī)律,優(yōu)化調(diào)節(jié)儲(chǔ)能系統(tǒng)的充放電特性,將實(shí)時(shí)的碳排放量守恒原則擴(kuò)展到調(diào)度周期內(nèi),本質(zhì)上是將儲(chǔ)能系統(tǒng)在能量層面的平移作用映射到碳排放層面的平移作用。

    儲(chǔ)能系統(tǒng)充電時(shí),是將分?jǐn)偨o用戶的CO2進(jìn)行存儲(chǔ)的過程;儲(chǔ)能系統(tǒng)放電時(shí),是將存儲(chǔ)的CO2重新分?jǐn)偨o用戶的過程。儲(chǔ)能系統(tǒng)的碳排放特性示意如圖1所示。儲(chǔ)能系統(tǒng)充放電過程滿足實(shí)時(shí)的碳排放量守恒原則。

    假設(shè)儲(chǔ)能系統(tǒng)在t1~t2時(shí)間段內(nèi)充電,在t2時(shí)刻開始充電。由于儲(chǔ)能系統(tǒng)充放電存在損耗,因此可將儲(chǔ)能系統(tǒng)充放電損耗帶來的碳排放反映到儲(chǔ)能系統(tǒng)內(nèi)外碳勢的差異上,儲(chǔ)能系統(tǒng)的內(nèi)外碳勢可表示為

    Eint,t2ES="""""""""nbsp;"" Pint,t1ESEint,t1ES+∑Ni=1Pi→ESNEi→ESNΔt+∑Kk=1PES,kGEES,kGΔtPint,t1ES+∑Ni=1Pi→ESNΔt+∑Kk=1PES,kGΔt·ηCES(2)

    Eext,t2ES=Eint,t2ES/ηDES(3)

    式中: Eint,t2ES——t2時(shí)刻儲(chǔ)能系統(tǒng)的內(nèi)碳勢;

    Eext,t2ES——t2時(shí)刻儲(chǔ)能系統(tǒng)的外碳勢;

    Pint,t1ES——t1時(shí)刻前儲(chǔ)能系統(tǒng)存儲(chǔ)的電量;

    Eint,t1ES——t1時(shí)刻儲(chǔ)能系統(tǒng)的內(nèi)碳勢;

    Pi→ESN——支路i→ES流過的正向有功功率;

    Ei→ESN——支路i→n的碳流密度;

    PES,kG——發(fā)電機(jī)組k在節(jié)點(diǎn)ES注入的有功功率;

    EES,kG——發(fā)電機(jī)組k的碳流密度;

    Δt——時(shí)段長度;

    ηCES——儲(chǔ)能系統(tǒng)充電效率;

    ηDES——儲(chǔ)能系統(tǒng)放電效率。

    2.2 負(fù)荷低碳特性

    根據(jù)動(dòng)態(tài)碳排放因子變化對負(fù)荷側(cè)可調(diào)度資源進(jìn)行激勵(lì)型響應(yīng)調(diào)度,引導(dǎo)用戶改變用能曲線,本質(zhì)上也是將可調(diào)度負(fù)荷在能量層面的平移作用映射到碳排放層面的平移作用。

    可調(diào)度負(fù)荷根據(jù)負(fù)荷響應(yīng)特性一般可分為可平移負(fù)荷、可轉(zhuǎn)移負(fù)荷和可削減負(fù)荷。其中,可平移負(fù)荷和可轉(zhuǎn)移負(fù)荷本質(zhì)上均是錯(cuò)峰調(diào)節(jié)特性,即通過錯(cuò)開電網(wǎng)負(fù)荷高峰期和用電高峰期實(shí)現(xiàn)電力調(diào)節(jié);可削減負(fù)荷本質(zhì)上是避峰調(diào)節(jié)特性,即在高峰時(shí)段降低總用電量。

    2.2.1 可平移負(fù)荷

    可平移負(fù)荷的具體表現(xiàn)是調(diào)整用戶用電的時(shí)段,但可能會(huì)受生產(chǎn)工序等固定時(shí)段順序的限制,調(diào)整時(shí)只能將用電曲線進(jìn)行整體性的平移,可平移負(fù)荷的數(shù)學(xué)模型可表示為

    [tTrans,begin,tTrans,end](4)

    PtTrans=βtTransPtTrans(5)

    ∑tTrans,endt=tTrans,beginPtTrans=∑Tt=1PTrans(6)

    ∑t1+tTrans,min-1t=t1βtTrans≥tTrans,min(βt1Trans-βt1-1Trans)(7)

    CTrans=RTrans∑Tt=1βtTrans(PTrans-PtTrans)(8)

    式中:tTrans,begin——可平移負(fù)荷的可接受平移開始時(shí)間;

    tTrans,end——可平移負(fù)荷的可接受平移結(jié)束時(shí)間;

    βtTrans——可平移負(fù)荷的平移狀態(tài)系數(shù),為0~1的變量,βtTrans=1表示t時(shí)段可平移負(fù)荷發(fā)生負(fù)荷平移;

    PtTrans——可平移負(fù)荷參與調(diào)節(jié)后的負(fù)荷功率;

    PtTrans,max——可平移負(fù)荷參與平移調(diào)節(jié)的最大功率;

    PtTrans,min——可平移負(fù)荷參與平移調(diào)節(jié)的最小功率;

    PTrans——t時(shí)段可平移負(fù)荷參與調(diào)節(jié)前的負(fù)荷功率;

    tTrans,min——可平移負(fù)荷參與平移調(diào)節(jié)需滿足的最小平移時(shí)間;

    CTrans——平移補(bǔ)償費(fèi)用;

    RTrans——可平移負(fù)荷的單位功率平移補(bǔ)償成本。

    2.2.2 可轉(zhuǎn)移負(fù)荷

    可轉(zhuǎn)移負(fù)荷的具體表現(xiàn)是用戶用電總量保持不變,而用電功率可在一定時(shí)間范圍內(nèi)靈活調(diào)節(jié),一般采用多種調(diào)整激勵(lì)措施,激勵(lì)用戶自發(fā)調(diào)整用電時(shí)段和用電量??赊D(zhuǎn)移負(fù)荷的數(shù)學(xué)模型可表示為

    [tShift,begin,tShift,end](9)

    βtShiftPtShift,max≤PtShift≤βtShiftPtShift,min(10)

    ∑tShift,endt=tShift,beginPtShift=∑Tt=1PShift(11)

    ∑t1+tShift,min-1t=t1βtShift≥tShift,min(βt1Shift-βt1-1Shift)(12)

    CShift=RShift∑Tt=1βtShift(PShift-PtShift)(13)

    式中:tShift,begin——可轉(zhuǎn)移負(fù)荷的可接受轉(zhuǎn)移開始時(shí)間;

    tShift,end——可轉(zhuǎn)移負(fù)荷的可接受轉(zhuǎn)移結(jié)束時(shí)間;

    βtShift——可轉(zhuǎn)移負(fù)荷的轉(zhuǎn)移狀態(tài)系數(shù),為0~1的變量,βtShift=1表示t時(shí)段可轉(zhuǎn)移負(fù)荷發(fā)生負(fù)荷轉(zhuǎn)移;

    PtShift——可轉(zhuǎn)移負(fù)荷參與調(diào)節(jié)后的負(fù)荷功率;

    PtShift,max——可轉(zhuǎn)移負(fù)荷參與轉(zhuǎn)移調(diào)節(jié)的最大功率;

    PtShift,min——可轉(zhuǎn)移負(fù)荷參與轉(zhuǎn)移調(diào)節(jié)的最小功率;

    PShift——t時(shí)段可轉(zhuǎn)移負(fù)荷參與調(diào)節(jié)前的負(fù)荷功率;

    tShift,min——熱可轉(zhuǎn)移負(fù)荷參與轉(zhuǎn)移調(diào)節(jié)需滿足的最小轉(zhuǎn)移時(shí)間;

    CShift——轉(zhuǎn)移補(bǔ)償費(fèi)用;

    RShift——可轉(zhuǎn)移負(fù)荷的單位功率轉(zhuǎn)移補(bǔ)償成本。

    2.2.3 可削減負(fù)荷

    可削減負(fù)荷的特點(diǎn)是用電功率可部分或全部削減,允許一定程度的負(fù)荷中斷、功率下降或減少時(shí)間運(yùn)行??上鳒p負(fù)荷的數(shù)學(xué)模型可表示為

    PtCut=(1-βtCutθtCut)PCut(14)

    CCut=RCut∑Tt=1βtCut(PCut-PtCut)(15)

    式中: PtCut——t時(shí)段電可削減負(fù)荷參與削減調(diào)節(jié)后的負(fù)荷功率;

    βtCut——可削減負(fù)荷的削減狀態(tài)系數(shù),為0~1的變量,βtCut=1表示t時(shí)段可削減負(fù)荷發(fā)生削減調(diào)節(jié);

    θtCut——t時(shí)段可削減負(fù)荷的削減比例,0≤θtCut≤1,θtCut=1表示可削減負(fù)荷被全部削減;

    PCut——可削減負(fù)荷參與削減調(diào)節(jié)前的功率;

    CCut——削減補(bǔ)償費(fèi)用;

    RCut——可削減負(fù)荷的單位功率削減補(bǔ)償成本。

    3 基于動(dòng)態(tài)碳排放因子的電力低碳優(yōu)化

    根據(jù)碳排放理論可知,節(jié)點(diǎn)碳勢隨系統(tǒng)中各發(fā)電機(jī)組流入該節(jié)點(diǎn)的功率不同而變化。有效利用動(dòng)態(tài)碳排放因子變化規(guī)律,優(yōu)化調(diào)度儲(chǔ)能系統(tǒng)和可調(diào)度負(fù)荷,從而有序引導(dǎo)用戶改變用能曲線,可以有效降低用戶碳排放。電力低碳優(yōu)化調(diào)度特性曲線如圖2所示。其實(shí)質(zhì)是將高碳排放因子時(shí)段的電力能源通過儲(chǔ)能系統(tǒng)和可調(diào)度負(fù)荷轉(zhuǎn)移至低碳排放因子時(shí)段。

    3.1 目標(biāo)函數(shù)

    本文的電力低碳優(yōu)化調(diào)度策略以碳排放量最少和系統(tǒng)運(yùn)行成本最低為多目標(biāo)優(yōu)化模型,系統(tǒng)運(yùn)行成本考慮外購能源成本、設(shè)備啟停成本、系統(tǒng)運(yùn)維成本、可調(diào)度負(fù)荷功率轉(zhuǎn)移補(bǔ)償成本和減碳補(bǔ)償費(fèi)用。目標(biāo)函數(shù)表示為

    Obj(F1)=∑Tt=1EnN(t)[PnN(t)-PextES(t)]+EextES(t)PextES(t)(16)

    Obj(F2)=F2,1+F2,2+F2,3+F2,4+F2,5(17)

    式中: F1——碳排放量最低目標(biāo)函數(shù);

    F2——系統(tǒng)運(yùn)行成本最低目標(biāo)函數(shù);

    F2,1——外購能源成本;

    F2,2——設(shè)備啟停成本;

    F2,3——系統(tǒng)運(yùn)維成本;

    F2,4——可調(diào)度負(fù)荷功率轉(zhuǎn)移補(bǔ)償成本;

    F2,5——減碳補(bǔ)償費(fèi)用。

    F2,1=∑Tt=1Cbuy(t)Pbuy(t)Δt(18)

    F2,2=∑Tt=1[δ(t)-δ(t-1)]Cst(t)Δt(19)

    F2,3=∑Mm=1∑Tt=1Cm(t)Pm(t)Δt(20)

    F2,4=∑Tt=1[CTrans(t)+CShift(t)+CCut(t)]Δt(21)

    F2,5=∑Tt=1CCarbon(t)RCarbon(t)(22)

    式中: Cbuy(t)——t時(shí)段單位購能成本;

    Pbuy(t)——t時(shí)段購能量;

    Cst(t)——t時(shí)段單位啟停成本;

    δ(t)——t時(shí)段設(shè)備啟停狀態(tài);

    δ(t-1)——t-1時(shí)段設(shè)備啟停狀態(tài);

    Cm(t)——t時(shí)段m設(shè)備的單位運(yùn)維成本;

    Pm(t)——t時(shí)段m設(shè)備的運(yùn)行功率;

    CCarbon(t)——t時(shí)段減碳補(bǔ)償費(fèi)用;

    RCarbon(t)——t時(shí)段減碳量。

    3.2 約束條件

    3.2.1 功率平衡約束

    系統(tǒng)滿足實(shí)時(shí)的功率平衡約束,表達(dá)式為

    ∑Nn=1Pi→nN(t)+∑Kk=1Pn,kG(t)=PES(t)+PLoad(t)+βTrans(t)PTrans(t)+βShift(t)PShift(t)+βCut(t)PCut(t)(23)

    式中: Pi→nN(t)——t時(shí)段節(jié)點(diǎn)i流向節(jié)點(diǎn)n的功率;

    Pn,kG(t)——t時(shí)段發(fā)電機(jī)組k在節(jié)點(diǎn)n注入的功率;

    PES(t)——t時(shí)段儲(chǔ)能系統(tǒng)吸收或釋放的功率;

    PLoad(t)——t時(shí)段常規(guī)負(fù)荷功率;

    PTrans(t)——t時(shí)段平移的功率;

    PShift(t)——t時(shí)段轉(zhuǎn)移的功率;

    PCut(t)——t時(shí)段削減的功率。

    3.2.2 儲(chǔ)能充放電功率約束

    儲(chǔ)能充放電過程中滿足充放電功率PES約束,即

    PminES≤PES≤PmaxES(24)

    式中: PmaxES——儲(chǔ)能系統(tǒng)的充放電功率上限;

    PminES——儲(chǔ)能系統(tǒng)的充放電功率下限。

    3.2.3 儲(chǔ)能容量約束

    儲(chǔ)能容量SES不能超過儲(chǔ)能可利用的最大容量,即

    SminES≤SES≤SmaxES(25)

    式中: SmaxES——儲(chǔ)能設(shè)備的容量上限;

    SminES——儲(chǔ)能設(shè)備的容量下限。

    3.2.4 儲(chǔ)能荷電狀態(tài)約束

    假設(shè)充電起始荷電狀態(tài)為10%,充電結(jié)束荷電狀態(tài)為90%,放電結(jié)束時(shí)荷電狀態(tài)和充電起始荷電狀態(tài)相等。

    SOC(t)=SOC(t-1)+ηCES·H(t)·PES(t)SES(t)+

    [H(t)-1]·PES(t)/ηDESSES(t)(26)

    式中: SOC(t)——t時(shí)刻儲(chǔ)能系統(tǒng)的荷電狀態(tài);

    SOC(t-1)——t-1時(shí)刻儲(chǔ)能系統(tǒng)的荷電狀態(tài);

    H(t)——t時(shí)刻儲(chǔ)能充/放電的0~1函數(shù),充電時(shí)為1,放電時(shí)為0。

    3.3 求解方法

    將多目標(biāo)優(yōu)化問題用數(shù)學(xué)模型描述為

    s.t.Dd(t)=0,d=1,2,…,nBminb(t)≤Bb(t)≤Bmaxb(t),b=1,2,…,m(27)

    minObj(F1),Obj(F2)(28)

    式中:""" Dd(t)——等式約束;

    Bb(t)——不等式約束;

    Bminb(t)、Bmaxb(t)——不等式約束上、下限。

    相較于單目標(biāo)優(yōu)化求解,多目標(biāo)優(yōu)化求解無法找到唯一的最優(yōu)解,而是得到一組Pareto前沿解。傳統(tǒng)的多目標(biāo)優(yōu)化算法是將多目標(biāo)按照一定的權(quán)重系數(shù)重新將多目標(biāo)組合成單目標(biāo)函數(shù),但是目標(biāo)的權(quán)重系數(shù)往往是通過經(jīng)驗(yàn)值確定的,其運(yùn)算結(jié)果并不理想。

    因此,本文采用魯棒性好、搜索能力強(qiáng)的帶精英策略的非支配遺傳算法(NSGA-II)對電力低碳多目標(biāo)優(yōu)化問題進(jìn)行求解。

    得到Pareto最優(yōu)解集后需要在Pareto解集中選擇合適的解。常規(guī)的方法是利用趨勢改變平緩的區(qū)間作為加工建議的參考,這種簡單的方式較為粗糙,且無法解釋如何獲得最優(yōu)解。最優(yōu)解的選擇可以根據(jù)決策者的需求或偏好進(jìn)行選擇,也可以通過綜合評價(jià)來進(jìn)一步對最優(yōu)解進(jìn)行篩選。

    本文利用模糊隸屬度函數(shù)評價(jià)決策者對每個(gè)非劣解的綜合滿意度,并選擇綜合隸屬度作為最優(yōu)折衷解。將模糊隸屬函數(shù)定義為

    Φ=∑Ii=1φi(29)

    φi=1,mi≤mi,minmi,max-mimi,max-mi,min,mi,min≤mi≤mi,max0,mi≥mi,max(30)

    式中: Φ——目標(biāo)函數(shù)的綜合隸屬度;

    φi——第i個(gè)目標(biāo)的隸屬度;

    mi——第i個(gè)目標(biāo)的函數(shù)值;

    mi,min——第i個(gè)目標(biāo)函數(shù)的最小值;

    mj,max——第i個(gè)目標(biāo)函數(shù)的最大值。

    4 算例分析

    4.1 算例系統(tǒng)架構(gòu)及參數(shù)

    本文采用改進(jìn)后的IEEE 14節(jié)點(diǎn)系統(tǒng)進(jìn)行算例分析。IEEE 14節(jié)點(diǎn)系統(tǒng)如圖3所示。本算例中設(shè)定λ=1,即系統(tǒng)網(wǎng)絡(luò)損耗全部分?jǐn)偨o用戶,并假設(shè)G1發(fā)電機(jī)組為系統(tǒng)功率平衡節(jié)點(diǎn)。IEEE 14節(jié)點(diǎn)系統(tǒng)的發(fā)電機(jī)組碳勢如表1所示。

    設(shè)定節(jié)點(diǎn)13儲(chǔ)能系統(tǒng)ES的最大充放電功率均為400 kW,最大可存儲(chǔ)電量為2 MWh,充放電轉(zhuǎn)換效率均為90%。初始狀態(tài)下,儲(chǔ)能系統(tǒng)SOC為10%,外碳勢為0.5 kgCO2/kWh。

    系統(tǒng)中各節(jié)點(diǎn)負(fù)荷和發(fā)電機(jī)組出力在每個(gè)時(shí)段的取值已給定。IEEE 14節(jié)點(diǎn)系統(tǒng)各節(jié)點(diǎn)的碳勢變化曲線如圖4所示。

    4.2 優(yōu)化調(diào)度結(jié)果分析

    本算例的優(yōu)化運(yùn)行目標(biāo)為全系統(tǒng)全天的總碳排放量最少和總運(yùn)行成本最低;優(yōu)化決策變量是1天24 h內(nèi)的儲(chǔ)能系統(tǒng)和可調(diào)度負(fù)荷的運(yùn)行狀態(tài),共24個(gè)時(shí)段。

    本算例使用MATLAB R2020a,通過帶精英策略的非支配遺傳算法對本文所述的電力低碳多目標(biāo)優(yōu)化問題進(jìn)行求解。本文算例的多目標(biāo)優(yōu)化Pareto前沿解如圖5所示;儲(chǔ)能系統(tǒng)和可調(diào)度負(fù)荷的優(yōu)化運(yùn)行結(jié)果如圖6所示。

    由圖6可見,在節(jié)點(diǎn)3低碳勢的11:00~13:00時(shí)段,儲(chǔ)能系統(tǒng)為充電狀態(tài),存儲(chǔ)電力網(wǎng)絡(luò)中的低碳勢能量;在節(jié)點(diǎn)3高碳勢的18:00~21:00時(shí)段,儲(chǔ)能系統(tǒng)為放電狀態(tài),釋放儲(chǔ)能系統(tǒng)已存儲(chǔ)的低碳勢能量??赊D(zhuǎn)移負(fù)荷由10節(jié)點(diǎn)高碳勢的05:00~06:00時(shí)段轉(zhuǎn)移至低碳勢的10:00和14:00時(shí)段??上鳒p負(fù)荷在節(jié)點(diǎn)14高碳勢的01:00~08:00時(shí)段和16:00~24:00時(shí)段進(jìn)行削減??善揭曝?fù)荷由節(jié)點(diǎn)9的節(jié)點(diǎn)碳勢變化趨勢將負(fù)荷整體向右平移了8 h。通過上述算例的電力低碳優(yōu)化結(jié)果,降低了系統(tǒng)高碳勢能量的消耗,更多地利用了系統(tǒng)中的低碳勢能量,實(shí)現(xiàn)了低碳減排。

    5 結(jié) 語

    有效利用動(dòng)態(tài)碳排放因子變化規(guī)律,能夠有序引導(dǎo)用戶改變用能曲線,對電力低碳發(fā)展具有重要的現(xiàn)實(shí)意義。本文提出一種基于動(dòng)態(tài)碳排放因子的電力低碳優(yōu)化調(diào)度策略,該策略以動(dòng)態(tài)碳排放因子為引導(dǎo)信號,通過儲(chǔ)能系統(tǒng)和可調(diào)度負(fù)荷的電力低碳特性引導(dǎo)用戶側(cè)進(jìn)行低碳激勵(lì)型響應(yīng),改變儲(chǔ)能系統(tǒng)的充放電和可調(diào)度負(fù)荷的運(yùn)行狀態(tài)。

    通過算例優(yōu)化,實(shí)現(xiàn)了對儲(chǔ)能系統(tǒng)和可調(diào)度負(fù)荷的低碳優(yōu)化調(diào)度,有序引導(dǎo)用戶用能曲線根據(jù)動(dòng)態(tài)碳排放因子進(jìn)行改變。期望本文研究內(nèi)容能夠?yàn)殡娏Φ吞純?yōu)化領(lǐng)域的研究者提供一些新的研究思路,推動(dòng)電力系統(tǒng)低碳優(yōu)化調(diào)度相關(guān)研究的進(jìn)步。

    【參 考 文 獻(xiàn)】

    [1] 辛保安.為實(shí)現(xiàn)“碳達(dá)峰、碳中和”目標(biāo)貢獻(xiàn)智慧和力量[J].國家電網(wǎng),2021(3):2-3.

    [2] 佚名.七部門印發(fā)《促進(jìn)綠色消費(fèi)實(shí)施方案》[J].再生資源與循環(huán)經(jīng)濟(jì),2022,15(2):6.

    [3] 張智剛,康重慶.碳中和目標(biāo)下構(gòu)建新型電力系統(tǒng)的挑戰(zhàn)與展望[J].中國電機(jī)工程學(xué)報(bào),2022,42(8):2806-2819.

    [4] CHENG Y, ZHANG N, WANG Y, et al. Modeling carbon emission flow in multiple energy systems[J].IEEE Transactions on Smart Grid,2019,10(4):3562-3574.

    [5] CHENG Y, ZHANG N, ZHANG B, et al. Low-carbon operation of multiple energy systems based on energy-carbon integrated prices[J].IEEE Transactions on Smart Grid,2019,11(2):1307-1318.

    [6] 康重慶,杜爾順,李姚旺,等.新型電力系統(tǒng)的“碳視角”:科學(xué)問題與研究框架[J].電網(wǎng)技術(shù),2022,46(3):821-833.

    [7] 張寧,李姚旺,黃俊輝,等.電力系統(tǒng)全環(huán)節(jié)碳計(jì)量方法與碳表系統(tǒng)[J].電力系統(tǒng)自動(dòng)化,2023,47(9):2-12.

    [8] 翁格平,任嬌蓉,姚艷,等.考慮時(shí)變電碳因子的園區(qū)綜合能源系統(tǒng)低碳經(jīng)濟(jì)調(diào)度[J].浙江電力,2022,41(10):106-114.

    [9] 周天睿,康重慶.基于碳排放流的配電系統(tǒng)低碳優(yōu)化運(yùn)行方法研究[J].全球能源互聯(lián)網(wǎng),2019,2(3):241-247.

    [10] 趙偉,熊正勇,潘艷,等.計(jì)及碳排放流的電力系統(tǒng)低碳規(guī)劃[J].電力系統(tǒng)自動(dòng)化,2023,47(9):23-33.

    [11] 周洋,王進(jìn),任勝男,等.計(jì)及碳捕集設(shè)備的電-氣互聯(lián)綜合能源系統(tǒng)低碳經(jīng)濟(jì)調(diào)度[J].電器與能效管理技術(shù),2020(4):106-112.

    [12] 任勝男,王進(jìn),龔羅文,等.考慮電轉(zhuǎn)氣熱回收與蓄熱罐的風(fēng)電消納低碳經(jīng)濟(jì)調(diào)度[J].電器與能效管理技術(shù),2020(11):103-110.

    [13] 劉哲遠(yuǎn),邢海軍,程浩忠,等.考慮碳排放流及需求響應(yīng)的綜合能源系統(tǒng)雙層優(yōu)化調(diào)度[J].高電壓技術(shù),2023,49(1):169-178.

    [14] 李鵬,王瑞,冀浩然,等.低碳化智能配電網(wǎng)規(guī)劃研究與展望[J].電力系統(tǒng)自動(dòng)化,2021,45(24):10-21.

    [15] 汪超群,陳懿,文福拴,等.電力系統(tǒng)碳排放流理論改進(jìn)與完善[J].電網(wǎng)技術(shù),2022,46(5):1683-1693.

    [16] 周天睿,康重慶,徐乾耀,等.碳排放流在電力網(wǎng)絡(luò)中分布的特性與機(jī)理分析[J].電力系統(tǒng)自動(dòng)化,2012,36(15):39-44.

    [17] 周天睿,康重慶,徐乾耀,等.電力系統(tǒng)碳排放流分析理論初探[J].電力系統(tǒng)自動(dòng)化,2012,36(7):38-43,85.

    收稿日期: 20240710

    猜你喜歡
    儲(chǔ)能
    基于DCS平臺(tái)的儲(chǔ)能調(diào)頻主要控制系統(tǒng)的開發(fā)與應(yīng)用
    吉林電力(2021年4期)2021-11-20 06:23:36
    相變儲(chǔ)能材料的應(yīng)用
    煤氣與熱力(2021年6期)2021-07-28 07:21:24
    相變儲(chǔ)能材料研究進(jìn)展
    煤氣與熱力(2021年4期)2021-06-09 06:16:56
    Al3+摻雜MnO2空心納米球的制備及其在儲(chǔ)能上的應(yīng)用
    儲(chǔ)能技術(shù)在電力系統(tǒng)中的應(yīng)用
    儲(chǔ)能真要起飛了?
    能源(2017年12期)2018-01-31 01:42:59
    中壓一體化儲(chǔ)能變流器的設(shè)計(jì)
    電測與儀表(2016年8期)2016-04-15 00:30:16
    各類重要儲(chǔ)能系統(tǒng)綜述
    直流儲(chǔ)能型準(zhǔn)Z源光伏并網(wǎng)逆變器
    直流系統(tǒng)混合儲(chǔ)能與蓄電池單獨(dú)儲(chǔ)能的對比
    身体一侧抽搐| 丰满的人妻完整版| 18禁在线播放成人免费| 最近2019中文字幕mv第一页| 免费一级毛片在线播放高清视频| 亚洲一区二区三区色噜噜| avwww免费| 观看美女的网站| 国产美女午夜福利| 久久久欧美国产精品| 精品一区二区三区视频在线观看免费| 国产精品免费一区二区三区在线| 老司机影院成人| 国产综合懂色| 精品一区二区三区人妻视频| 免费高清视频大片| 日本一二三区视频观看| 成人欧美大片| 人妻夜夜爽99麻豆av| 久久99热这里只有精品18| 欧美性猛交╳xxx乱大交人| 国产人妻一区二区三区在| 久久人妻av系列| 亚洲成人av在线免费| 少妇丰满av| 日韩欧美免费精品| 久久草成人影院| 久久久久免费精品人妻一区二区| 亚洲精品日韩av片在线观看| 少妇猛男粗大的猛烈进出视频 | 在线国产一区二区在线| 午夜福利高清视频| 乱人视频在线观看| 欧美成人精品欧美一级黄| 91狼人影院| 此物有八面人人有两片| 中国国产av一级| 中文字幕精品亚洲无线码一区| 一级毛片aaaaaa免费看小| 国产午夜福利久久久久久| 国产精品,欧美在线| 美女内射精品一级片tv| 99热这里只有是精品50| 12—13女人毛片做爰片一| 国产欧美日韩精品亚洲av| 久久久色成人| 午夜免费男女啪啪视频观看 | 欧美成人精品欧美一级黄| 不卡视频在线观看欧美| 国产伦精品一区二区三区视频9| 日韩av在线大香蕉| 校园春色视频在线观看| 三级男女做爰猛烈吃奶摸视频| 国产高清激情床上av| 亚洲欧美日韩卡通动漫| 久久久a久久爽久久v久久| 国产老妇女一区| 国产精品一及| 久久久久九九精品影院| 可以在线观看毛片的网站| 免费在线观看影片大全网站| 露出奶头的视频| 午夜爱爱视频在线播放| 国产精品亚洲美女久久久| 99国产极品粉嫩在线观看| 色综合站精品国产| 中出人妻视频一区二区| 黄色一级大片看看| 色吧在线观看| 国产精品美女特级片免费视频播放器| 成人一区二区视频在线观看| 国产av在哪里看| 亚洲无线在线观看| 少妇熟女欧美另类| 亚洲人成网站在线观看播放| 深夜a级毛片| 狠狠狠狠99中文字幕| 女的被弄到高潮叫床怎么办| 老熟妇乱子伦视频在线观看| 国产精品国产三级国产av玫瑰| 成人无遮挡网站| 级片在线观看| 麻豆久久精品国产亚洲av| 天堂动漫精品| aaaaa片日本免费| 亚洲性久久影院| 亚洲美女黄片视频| 午夜福利在线观看吧| 久久久久精品国产欧美久久久| 国语自产精品视频在线第100页| 亚洲成人中文字幕在线播放| 亚洲五月天丁香| 国产aⅴ精品一区二区三区波| 不卡一级毛片| 久久久久久大精品| 亚洲精品456在线播放app| 少妇熟女aⅴ在线视频| 亚洲精品成人久久久久久| 观看美女的网站| 性欧美人与动物交配| 成人高潮视频无遮挡免费网站| 成人三级黄色视频| 午夜老司机福利剧场| 真实男女啪啪啪动态图| 日日啪夜夜撸| 少妇猛男粗大的猛烈进出视频 | 婷婷精品国产亚洲av| 国产精品免费一区二区三区在线| 美女 人体艺术 gogo| 中国美白少妇内射xxxbb| 日本色播在线视频| 国产精华一区二区三区| 在线播放无遮挡| 亚洲自拍偷在线| 日韩av在线大香蕉| 免费一级毛片在线播放高清视频| 国产精品福利在线免费观看| 免费搜索国产男女视频| 啦啦啦韩国在线观看视频| 成人av一区二区三区在线看| 亚洲无线在线观看| 日本三级黄在线观看| 精品不卡国产一区二区三区| 国产欧美日韩精品亚洲av| 午夜激情福利司机影院| 一个人看的www免费观看视频| 色吧在线观看| 97超级碰碰碰精品色视频在线观看| 日韩国内少妇激情av| 久久婷婷人人爽人人干人人爱| 97热精品久久久久久| 又爽又黄无遮挡网站| 亚洲精品乱码久久久v下载方式| 中文字幕免费在线视频6| 精品久久久久久久久久免费视频| 亚洲欧美成人精品一区二区| 在线a可以看的网站| 蜜臀久久99精品久久宅男| 免费电影在线观看免费观看| 天美传媒精品一区二区| 国产麻豆成人av免费视频| 日本在线视频免费播放| 日韩欧美在线乱码| 精品熟女少妇av免费看| 午夜a级毛片| 精品免费久久久久久久清纯| 欧美xxxx性猛交bbbb| 中文字幕av在线有码专区| 国产乱人视频| 少妇猛男粗大的猛烈进出视频 | 日本-黄色视频高清免费观看| 亚洲国产欧美人成| 亚洲图色成人| 大型黄色视频在线免费观看| 亚洲五月天丁香| 成人综合一区亚洲| 国产 一区精品| 欧美成人一区二区免费高清观看| 午夜免费男女啪啪视频观看 | 午夜精品一区二区三区免费看| 亚洲精品日韩在线中文字幕 | 亚洲人成网站在线播| av.在线天堂| 国产一区二区亚洲精品在线观看| 精品99又大又爽又粗少妇毛片| 夜夜看夜夜爽夜夜摸| 久久久久久久午夜电影| 晚上一个人看的免费电影| 欧美日韩一区二区视频在线观看视频在线 | 欧美又色又爽又黄视频| 性色avwww在线观看| 麻豆一二三区av精品| 看非洲黑人一级黄片| 久久国内精品自在自线图片| 欧美日韩国产亚洲二区| 深夜精品福利| 中文字幕免费在线视频6| 身体一侧抽搐| 免费搜索国产男女视频| 免费av观看视频| 少妇丰满av| 日韩亚洲欧美综合| 国产精品乱码一区二三区的特点| 国产成人福利小说| 精品久久久久久久久久久久久| 国产成人精品久久久久久| 成人漫画全彩无遮挡| 亚洲精品乱码久久久v下载方式| 悠悠久久av| 免费观看人在逋| 国产在视频线在精品| 久久热精品热| 亚洲va在线va天堂va国产| 亚洲成av人片在线播放无| 成人国产麻豆网| 99在线视频只有这里精品首页| 亚洲性久久影院| 91麻豆精品激情在线观看国产| 亚洲人与动物交配视频| 在线观看免费视频日本深夜| 男女之事视频高清在线观看| 中文字幕av在线有码专区| 婷婷六月久久综合丁香| 久久精品国产99精品国产亚洲性色| 伦理电影大哥的女人| 99在线人妻在线中文字幕| 国产精品伦人一区二区| 亚洲精品亚洲一区二区| 欧美xxxx黑人xx丫x性爽| 久久久欧美国产精品| 免费av毛片视频| 久久精品国产亚洲网站| 国内精品美女久久久久久| 一夜夜www| 人妻丰满熟妇av一区二区三区| 搡老岳熟女国产| 99久国产av精品国产电影| 国产精品99久久久久久久久| 长腿黑丝高跟| 高清毛片免费看| 亚洲av中文av极速乱| 夜夜爽天天搞| 91在线观看av| 国产精品三级大全| 校园春色视频在线观看| 国产午夜精品久久久久久一区二区三区 | 日本三级黄在线观看| 日韩av不卡免费在线播放| 嫩草影院入口| 亚洲欧美日韩东京热| АⅤ资源中文在线天堂| 亚洲天堂国产精品一区在线| 精品久久久久久久久亚洲| 成人av在线播放网站| 欧美不卡视频在线免费观看| 俄罗斯特黄特色一大片| 免费一级毛片在线播放高清视频| 九九久久精品国产亚洲av麻豆| 免费在线观看成人毛片| 91久久精品电影网| 美女cb高潮喷水在线观看| 成人永久免费在线观看视频| 女生性感内裤真人,穿戴方法视频| 精品99又大又爽又粗少妇毛片| 中文字幕av成人在线电影| 97热精品久久久久久| 国产午夜精品论理片| 亚洲一级一片aⅴ在线观看| 成年女人毛片免费观看观看9| 身体一侧抽搐| 免费观看的影片在线观看| 国产探花极品一区二区| 久久久精品欧美日韩精品| 久久人妻av系列| 18禁黄网站禁片免费观看直播| 亚洲国产精品成人久久小说 | 色吧在线观看| 欧美成人免费av一区二区三区| 久久热精品热| 国产淫片久久久久久久久| 不卡视频在线观看欧美| 久久综合国产亚洲精品| 免费人成在线观看视频色| 国内揄拍国产精品人妻在线| 晚上一个人看的免费电影| 嫩草影视91久久| 亚洲精品久久国产高清桃花| 精品福利观看| 成人午夜高清在线视频| 国产午夜精品论理片| 床上黄色一级片| 又爽又黄a免费视频| 搡老岳熟女国产| 狂野欧美白嫩少妇大欣赏| 成人亚洲欧美一区二区av| 日本爱情动作片www.在线观看 | av在线天堂中文字幕| 精品久久久噜噜| 婷婷六月久久综合丁香| 久久草成人影院| 欧美日韩一区二区视频在线观看视频在线 | 亚洲成人久久性| 一区福利在线观看| 偷拍熟女少妇极品色| 秋霞在线观看毛片| 国产大屁股一区二区在线视频| av天堂中文字幕网| 亚洲真实伦在线观看| 亚洲最大成人手机在线| 亚洲人成网站在线播放欧美日韩| 一个人观看的视频www高清免费观看| 人妻夜夜爽99麻豆av| 亚洲18禁久久av| 欧美高清性xxxxhd video| 99国产极品粉嫩在线观看| 欧美色视频一区免费| 欧美日韩综合久久久久久| 一夜夜www| 亚洲欧美成人精品一区二区| 日韩强制内射视频| 日韩欧美在线乱码| 床上黄色一级片| av在线老鸭窝| 搡女人真爽免费视频火全软件 | 国产视频一区二区在线看| av在线播放精品| 国产精品人妻久久久久久| 日本精品一区二区三区蜜桃| 亚洲人成网站高清观看| 深夜精品福利| 精品久久久噜噜| 91麻豆精品激情在线观看国产| 亚洲人成网站在线播| 亚洲欧美日韩无卡精品| 欧美中文日本在线观看视频| 村上凉子中文字幕在线| 美女高潮的动态| 日韩精品青青久久久久久| av女优亚洲男人天堂| 免费观看在线日韩| 日本一二三区视频观看| 亚洲第一区二区三区不卡| 亚洲精品影视一区二区三区av| 少妇熟女aⅴ在线视频| 深爱激情五月婷婷| 亚洲性久久影院| 噜噜噜噜噜久久久久久91| 99久久精品一区二区三区| 日韩三级伦理在线观看| 成人美女网站在线观看视频| 亚洲精品色激情综合| 国产乱人偷精品视频| 久久精品夜夜夜夜夜久久蜜豆| 99久国产av精品| 村上凉子中文字幕在线| 蜜桃亚洲精品一区二区三区| 日本色播在线视频| 午夜福利成人在线免费观看| 在线观看一区二区三区| 日本三级黄在线观看| 欧美bdsm另类| 欧美另类亚洲清纯唯美| 久久精品国产亚洲av香蕉五月| 午夜精品国产一区二区电影 | 亚洲熟妇中文字幕五十中出| 国内精品久久久久精免费| 国产一区亚洲一区在线观看| 精品久久久久久久久久免费视频| 亚洲av成人av| 99精品在免费线老司机午夜| 麻豆国产97在线/欧美| 国产色婷婷99| 少妇的逼好多水| 人妻丰满熟妇av一区二区三区| 黑人高潮一二区| 日韩欧美国产在线观看| 99久久中文字幕三级久久日本| av黄色大香蕉| 国产成年人精品一区二区| 国产精品一区二区免费欧美| 亚洲熟妇中文字幕五十中出| av专区在线播放| 哪里可以看免费的av片| 麻豆一二三区av精品| 久久久久久久午夜电影| 搡老熟女国产l中国老女人| 俄罗斯特黄特色一大片| 亚洲欧美日韩东京热| 岛国在线免费视频观看| 国产精品1区2区在线观看.| 六月丁香七月| 免费av不卡在线播放| 久久韩国三级中文字幕| 国产老妇女一区| 男人的好看免费观看在线视频| 在线观看免费视频日本深夜| 在现免费观看毛片| 日韩av不卡免费在线播放| 久久精品夜色国产| 国产午夜精品久久久久久一区二区三区 | 看黄色毛片网站| 一级毛片电影观看 | 在线看三级毛片| av在线观看视频网站免费| 麻豆av噜噜一区二区三区| 亚洲自拍偷在线| 久久草成人影院| 少妇丰满av| 日本三级黄在线观看| 丰满的人妻完整版| 国内精品久久久久精免费| 亚洲熟妇熟女久久| 在线播放无遮挡| 久久综合国产亚洲精品| 又黄又爽又免费观看的视频| 18禁在线无遮挡免费观看视频 | 非洲黑人性xxxx精品又粗又长| 波多野结衣高清无吗| 国产视频内射| 欧美bdsm另类| 国产综合懂色| 亚洲精品粉嫩美女一区| 亚洲国产日韩欧美精品在线观看| 日本爱情动作片www.在线观看 | 在线免费十八禁| 一进一出抽搐gif免费好疼| 可以在线观看的亚洲视频| 日韩一本色道免费dvd| 国产日本99.免费观看| 久久人人精品亚洲av| 特大巨黑吊av在线直播| www.色视频.com| 国产又黄又爽又无遮挡在线| 国产精品三级大全| 国产伦在线观看视频一区| 两个人视频免费观看高清| 亚洲av二区三区四区| 欧美成人一区二区免费高清观看| 色哟哟哟哟哟哟| 最近2019中文字幕mv第一页| 国产 一区精品| 蜜桃久久精品国产亚洲av| 国产精品一及| a级一级毛片免费在线观看| 中文字幕熟女人妻在线| 国产精品国产高清国产av| 欧美成人一区二区免费高清观看| 麻豆久久精品国产亚洲av| 欧美xxxx性猛交bbbb| 午夜影院日韩av| 欧美成人精品欧美一级黄| 久久精品国产亚洲av香蕉五月| 国产又黄又爽又无遮挡在线| 久久精品综合一区二区三区| 一个人看视频在线观看www免费| av视频在线观看入口| 三级男女做爰猛烈吃奶摸视频| 国产亚洲精品久久久com| ponron亚洲| 最新在线观看一区二区三区| 男女做爰动态图高潮gif福利片| 国产综合懂色| 丰满乱子伦码专区| 欧美激情久久久久久爽电影| 午夜福利视频1000在线观看| 午夜福利高清视频| 老司机午夜福利在线观看视频| 男女做爰动态图高潮gif福利片| 亚洲中文字幕日韩| 成人性生交大片免费视频hd| 此物有八面人人有两片| 成熟少妇高潮喷水视频| 日韩三级伦理在线观看| 日韩欧美在线乱码| 免费搜索国产男女视频| 国内精品美女久久久久久| 非洲黑人性xxxx精品又粗又长| 久久国产乱子免费精品| 中文字幕人妻熟人妻熟丝袜美| 久久精品影院6| 亚洲国产日韩欧美精品在线观看| 毛片一级片免费看久久久久| av在线蜜桃| 国产视频内射| 久久99热这里只有精品18| 中文字幕久久专区| 偷拍熟女少妇极品色| 国产黄片美女视频| 99在线视频只有这里精品首页| 亚洲真实伦在线观看| 国产精品久久电影中文字幕| 久久精品影院6| 搡老熟女国产l中国老女人| 久久韩国三级中文字幕| 1000部很黄的大片| 男人狂女人下面高潮的视频| 变态另类丝袜制服| 不卡视频在线观看欧美| 男人的好看免费观看在线视频| 大型黄色视频在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 国产精品免费一区二区三区在线| 国产亚洲精品综合一区在线观看| 91久久精品国产一区二区成人| av.在线天堂| 99热全是精品| 亚洲美女视频黄频| 联通29元200g的流量卡| 一个人免费在线观看电影| 99视频精品全部免费 在线| 日韩 亚洲 欧美在线| 亚洲精品日韩av片在线观看| 国产在线精品亚洲第一网站| 深爱激情五月婷婷| 中文亚洲av片在线观看爽| 狂野欧美白嫩少妇大欣赏| 少妇被粗大猛烈的视频| 中文字幕久久专区| 狠狠狠狠99中文字幕| 精品熟女少妇av免费看| 欧美丝袜亚洲另类| www.色视频.com| a级一级毛片免费在线观看| 亚洲av熟女| 国产精品一区二区性色av| 性欧美人与动物交配| 国产精品,欧美在线| 亚洲欧美成人精品一区二区| 一区二区三区免费毛片| 神马国产精品三级电影在线观看| 亚洲精品乱码久久久v下载方式| 变态另类丝袜制服| 国产毛片a区久久久久| 99久久精品国产国产毛片| 国产三级在线视频| 久久久a久久爽久久v久久| av在线天堂中文字幕| 女同久久另类99精品国产91| 在线看三级毛片| 国产精品久久久久久精品电影| 国产伦精品一区二区三区四那| 69人妻影院| 不卡一级毛片| 精品午夜福利视频在线观看一区| 亚洲av成人av| 黄色一级大片看看| 久99久视频精品免费| 少妇被粗大猛烈的视频| 一级毛片aaaaaa免费看小| av视频在线观看入口| 免费在线观看影片大全网站| 欧美极品一区二区三区四区| 国产一区二区激情短视频| 伦理电影大哥的女人| 人妻少妇偷人精品九色| 搡女人真爽免费视频火全软件 | 国产精品一区二区三区四区久久| 两性午夜刺激爽爽歪歪视频在线观看| 毛片女人毛片| 美女免费视频网站| 日韩大尺度精品在线看网址| 欧美zozozo另类| 国内精品美女久久久久久| 极品教师在线视频| 亚洲精品粉嫩美女一区| 一个人看的www免费观看视频| 欧美xxxx性猛交bbbb| 久久99热这里只有精品18| 18+在线观看网站| 97热精品久久久久久| 亚洲精品一区av在线观看| 中文字幕人妻熟人妻熟丝袜美| 99九九线精品视频在线观看视频| 亚洲电影在线观看av| 国产69精品久久久久777片| 亚洲精品色激情综合| 久久精品国产亚洲av涩爱 | 日日摸夜夜添夜夜添小说| 99国产精品一区二区蜜桃av| 91午夜精品亚洲一区二区三区| 内地一区二区视频在线| 欧美一区二区国产精品久久精品| 国产精品1区2区在线观看.| 久久久久性生活片| 国产精品国产三级国产av玫瑰| 欧美激情国产日韩精品一区| 麻豆国产av国片精品| 成人国产麻豆网| 国产精品人妻久久久久久| 高清毛片免费观看视频网站| 久久国产乱子免费精品| 国产欧美日韩精品一区二区| 99久久无色码亚洲精品果冻| 此物有八面人人有两片| 久久久a久久爽久久v久久| 嫩草影院新地址| 国内少妇人妻偷人精品xxx网站| 日本在线视频免费播放| 欧美一区二区亚洲| av在线蜜桃| 亚洲精品乱码久久久v下载方式| 舔av片在线| 欧美日韩一区二区视频在线观看视频在线 | 亚洲成人久久爱视频| 亚洲综合色惰| 伦精品一区二区三区| 国产精品嫩草影院av在线观看| 男女下面进入的视频免费午夜| 日韩欧美在线乱码| 久久九九热精品免费| 日韩精品有码人妻一区| 亚洲国产日韩欧美精品在线观看| 久久午夜亚洲精品久久| 亚洲欧美成人综合另类久久久 | 国产伦精品一区二区三区四那| 97热精品久久久久久| 国产 一区精品| 天堂av国产一区二区熟女人妻| 成人无遮挡网站| 搡老岳熟女国产| 亚洲中文日韩欧美视频| 国产亚洲av嫩草精品影院| 在线观看一区二区三区| 中文资源天堂在线| 好男人在线观看高清免费视频| 久久亚洲精品不卡| 中文资源天堂在线| 亚洲国产色片| 日韩制服骚丝袜av| 伦精品一区二区三区| 国产成人aa在线观看| 免费不卡的大黄色大毛片视频在线观看 | 亚洲综合色惰| 国产精品野战在线观看| 日韩中字成人|