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

    考慮階梯碳交易和最優(yōu)建設(shè)時(shí)序的園區(qū)綜合能源系統(tǒng)規(guī)劃

    2023-12-25 05:43:06江訓(xùn)譜呂施霖王健張瑩穎包哲靜于淼
    電測(cè)與儀表 2023年12期
    關(guān)鍵詞:碳達(dá)峰階梯排放量

    江訓(xùn)譜,呂施霖,王健,張瑩穎,包哲靜,于淼

    (1.浙江大學(xué)工程師學(xué)院,杭州 310015; 2.國(guó)網(wǎng)浙江省電力有限公司綜合服務(wù)分公司,杭州 310005; 3.浙江大學(xué)電氣工程學(xué)院,杭州 310027)

    0 引言

    近年來,隨著能源需求不斷地增長(zhǎng),如何實(shí)現(xiàn)能源的可持續(xù)發(fā)展得到了全社會(huì)的關(guān)注,我國(guó)也提出了2030 年實(shí)現(xiàn)碳達(dá)峰、2060 年實(shí)現(xiàn)碳中和的“雙碳”戰(zhàn)略目標(biāo)[1]。隨著能源需求的增長(zhǎng),熱電聯(lián)產(chǎn)( CHP) 等新的技術(shù)也在迅速地發(fā)展[2],電、熱、氣等多種能源得以相互耦合和協(xié)同,園區(qū)綜合能源系統(tǒng)( PIES) 可以實(shí)現(xiàn)多能源互補(bǔ),從而實(shí)現(xiàn)能源綜合利用,已經(jīng)成為了實(shí)現(xiàn)能源可持續(xù)發(fā)展的重要途徑之一。由于規(guī)劃的年限比較長(zhǎng),根據(jù)不同的能源需求和供給情況,合理的建設(shè)時(shí)序可以最大限度地提高能源利用效率,減少投資浪費(fèi)和降低環(huán)境污染。

    PIES 的規(guī)劃問題關(guān)系到綜合能源系統(tǒng)運(yùn)行的經(jīng)濟(jì)性、可靠性等多個(gè)方面,受到學(xué)者的廣泛關(guān)注。文獻(xiàn)[3]針對(duì)系統(tǒng)運(yùn)行的經(jīng)濟(jì)性建立了考慮激勵(lì)型需求響應(yīng)的多區(qū)域規(guī)劃模型;文獻(xiàn)[4]提出了通過多區(qū)域冷熱聯(lián)設(shè)備和儲(chǔ)能設(shè)備的協(xié)調(diào)運(yùn)行來改善系統(tǒng)運(yùn)行的經(jīng)濟(jì)性;文獻(xiàn)[5]提出了考慮經(jīng)濟(jì)性、可持續(xù)性和可靠性協(xié)同優(yōu)化的電氣系統(tǒng)聯(lián)合規(guī)劃模型;文獻(xiàn)[6-7]提出了在長(zhǎng)周期的背景下考慮經(jīng)濟(jì)性的多階段規(guī)劃模型; 文獻(xiàn)[8]提出了考慮源菏互動(dòng)和需求響應(yīng)的園區(qū)雙層規(guī)劃模型以提高系統(tǒng)運(yùn)行的自主性;文獻(xiàn)[9]提出了考慮風(fēng)電不確定性的混合魯棒區(qū)間優(yōu)化的模型以協(xié)調(diào)系統(tǒng)運(yùn)行的經(jīng)濟(jì)性與可靠性。但是上述研究鮮有考慮系統(tǒng)設(shè)備最優(yōu)的建設(shè)時(shí)序問題。

    PIES 集多種能源于一體,實(shí)現(xiàn)能源的多能互補(bǔ)、協(xié)同優(yōu)化,從而降低碳排放。在“雙碳”目標(biāo)的背景下,有必要將碳排放的影響納入園區(qū)規(guī)劃中進(jìn)行考慮。文獻(xiàn)[10]提出建立電-碳價(jià)格相關(guān)性和氫儲(chǔ)能的綜合能源系統(tǒng)以降低碳排放的模型;文獻(xiàn)[11-12]研究碳交易對(duì)電氣綜合系統(tǒng)運(yùn)行的影響; 文獻(xiàn)[13]研究碳交易在綜合能源系統(tǒng)中的應(yīng)用及未來的發(fā)展;文獻(xiàn)[14-15]建立了階梯碳交易模型,研究階梯碳交易機(jī)制對(duì)系統(tǒng)容量配置的影響;文獻(xiàn)[16]在中長(zhǎng)期的多階段規(guī)劃引入氫能和碳交易的減排手段來提升園區(qū)的經(jīng)濟(jì)性與減排效果;文獻(xiàn)[17-18]建立綠證與碳交易聯(lián)合交易的模型以促進(jìn)新能源的消納量和降低系統(tǒng)的碳排量。上述研究在園區(qū)中引入碳交易機(jī)制,有效地提高了系統(tǒng)的低碳性。

    文中在上述研究基礎(chǔ)上,綜合考慮階梯碳交易和最優(yōu)建設(shè)時(shí)序,建立了階梯碳交易成本模型,目標(biāo)函數(shù)取為全壽命周期總費(fèi)用折現(xiàn)到投資初年的現(xiàn)值,并結(jié)合碳達(dá)峰、碳中和的“雙碳”目標(biāo)設(shè)置碳排放約束條件,根據(jù)PIES 的實(shí)際發(fā)展情況去考慮負(fù)荷在不同階段的增長(zhǎng),進(jìn)一步考慮PIES 的建設(shè)時(shí)序,將投建年份設(shè)置為待尋優(yōu)量,求解得出最佳建設(shè)時(shí)序; 最后,討論了碳交易基準(zhǔn)價(jià)格和投建次數(shù)對(duì)園區(qū)規(guī)劃結(jié)果的影響。

    1 園區(qū)綜合能源系統(tǒng)建模

    1.1 PIES 的結(jié)構(gòu)示意圖

    PIES 以能源耦合、多種能源相互關(guān)聯(lián)和互補(bǔ)為特征,包含儲(chǔ)能設(shè)備、能源生產(chǎn)與轉(zhuǎn)換設(shè)備,滿足系統(tǒng)終端用戶的電/熱/氣的負(fù)荷要求。為表示所有設(shè)備的能量流動(dòng)關(guān)系,文中采用電、熱、氣統(tǒng)一母線[19]表示PIES的基本結(jié)構(gòu),如圖1 所示。通過該結(jié)構(gòu)示意圖,可以清晰直觀地看到各個(gè)設(shè)備的連接、耦合關(guān)系。

    圖1 園區(qū)綜合能源系統(tǒng)結(jié)構(gòu)示意圖Fig.1 Schematic diagram of PIES structure

    1.2 PIES 的設(shè)備建模

    文中考慮PIES 的候選規(guī)劃的能源生產(chǎn)設(shè)備有光伏( PV) ,能源轉(zhuǎn)換設(shè)備有熱泵( HP) 、熱電聯(lián)產(chǎn)、燃?xì)忮仩t( GB) ,儲(chǔ)能設(shè)備有儲(chǔ)電裝置( ES) 、蓄熱裝置( HS) ,各個(gè)候選設(shè)備的建模如下:

    1) 光伏。

    基于典型日的光照數(shù)據(jù),假設(shè)PV 以最大功率點(diǎn)追蹤[20]模式( MPPT) 工作,PV 在MPPT 的最大輸出功率曲線是PV 的預(yù)測(cè)輸出功率曲線,PV 在任意時(shí)刻的實(shí)際輸出功率均不能超過該最大輸出功率,如式( 1)所示:

    式中PPV(t) 表示t時(shí)段的PV 實(shí)際輸出功率;表示t時(shí)段PV 在MPPT 下的最大輸出功率。

    2) 能源轉(zhuǎn)換設(shè)備。

    如果將各種能源轉(zhuǎn)換設(shè)備的輸入功率、輸出功率、轉(zhuǎn)換效率分別表示為Ii(t) 、Oi(t) 和ηi,則能源轉(zhuǎn)換設(shè)備的出力特性及其出力約束如式(2) 和式(3) 所示:

    其中,式(2) 表示各能源轉(zhuǎn)換設(shè)備的輸入輸出功率之間的關(guān)系;式(3) 表示設(shè)備出力的上下限約束;指各設(shè)備的輸出功率的上限,由設(shè)備當(dāng)下的裝機(jī)容量決定。

    3) 儲(chǔ)能設(shè)備。

    儲(chǔ)電、蓄熱裝置的能量轉(zhuǎn)換關(guān)系相似,文中采用一個(gè)廣義儲(chǔ)能系統(tǒng)的通用模型[21]對(duì)二者的運(yùn)行特性進(jìn)行描述,儲(chǔ)能裝置的運(yùn)行特性模型如下:

    其中,Pcha(t) 和Pdis(t) 表示儲(chǔ)能設(shè)備的充、放能功率;和表示儲(chǔ)能設(shè)備的充、放能功率上限,由儲(chǔ)能設(shè)備當(dāng)下的裝機(jī)容量決定;SOC(t) 表示儲(chǔ)能設(shè)備在t時(shí)段的荷電狀態(tài)( State of Charge,SOC) ; ηcha和ηdis是儲(chǔ)能設(shè)備的充、放能效率;WESS表示儲(chǔ)能設(shè)備的容量;SOCmin和SOCmax表示儲(chǔ)能荷電狀態(tài)的上下限;Δt表示時(shí)間間隔,文中取Δt=1 h。式(4) 和式(5) 表示儲(chǔ)能設(shè)備的充、放能功率的上下限約束; 式( 6) 表示充、放能過程不能同時(shí)進(jìn)行; 式( 7) 表示儲(chǔ)能設(shè)備的SOC 隨充放電過程的變化;式(8) 表示儲(chǔ)能荷電狀態(tài)的上下限約束;式( 9) 表示儲(chǔ)能裝置的周期始、末荷能狀態(tài)相同,考慮儲(chǔ)能裝置以日為周期循環(huán),所以T0=24 。

    2 園區(qū)綜合能源系統(tǒng)的階梯碳交易成本計(jì)算模型

    2.1 PIES 的碳排放量和碳排放配額

    碳交易政策最初是聯(lián)合國(guó)為了應(yīng)對(duì)氣候問題而構(gòu)建的一種國(guó)際貿(mào)易體系,在我國(guó)處于逐漸發(fā)展的階段[22]。在碳交易市場(chǎng)中,碳市場(chǎng)將碳排放權(quán)分配給企業(yè)并允許其自由交易。在電力行業(yè),我國(guó)主要是選擇免費(fèi)分配的方法用以分配碳配額[23],文中選取基準(zhǔn)線法[15]以確定園區(qū)的無償碳配額,認(rèn)為PIES 的碳排放主要來自外購(gòu)電力、GB、CHP 這三個(gè)部分,并考慮外購(gòu)電力全部來自于火力發(fā)電。

    PIES 的實(shí)際碳排放量如式(10) ~式(13) 所示:

    其中,En、Egrid,n、EGB,n、ECHP,n分別表示第n年P(guān)IES、外購(gòu)電力、GB、CHP 的實(shí)際碳排放量; βe、βh分別表示單位電量、單位熱量的碳排放量;Pgrid,n(t) 表示第n年的第t時(shí)段PIES 向電網(wǎng)的購(gòu)電功率;HGB,n(t) 表示第n年的第t時(shí)段燃?xì)忮仩t的熱輸出功率;HCHP,n(t) 、PCHP,n(t) 表示第n年的第t時(shí)段CHP 機(jī)組的熱、電輸出功率;文中碳交易費(fèi)用以年為時(shí)間尺度進(jìn)行結(jié)算,所以T=8760; φeh是CHP 機(jī)組發(fā)電量折算到發(fā)熱量的折算系數(shù)[23]。

    PIES 的碳排放配額如式(14) 所示:

    PIES 第n年考慮碳配額后的交易碳排放量為Ef,n,表示為:

    2.2 階梯碳交易的計(jì)算模型

    文中建立了階梯碳交易的模型,規(guī)定了若干個(gè)交易碳排放區(qū)間,以及交易碳排放量Ef,n在各個(gè)交易碳排放區(qū)間上對(duì)應(yīng)的價(jià)格。當(dāng)Ef,n為負(fù)數(shù)時(shí),表示實(shí)際碳排放量小于碳配額,此時(shí)PIES 可以將多余的碳配額出售以獲取收益;當(dāng)Ef,n為正數(shù)時(shí),表示實(shí)際碳排放量大于碳配額,那么PIES 需要向碳交易市場(chǎng)購(gòu)買所需的碳配額,且所需購(gòu)買的碳配額越多,碳交易價(jià)格也會(huì)越高。

    階梯碳交易的交易區(qū)間及其個(gè)數(shù)通常由上級(jí)碳排放監(jiān)管部分來制定,園區(qū)是碳交易的執(zhí)行者。在制定過程中,通常會(huì)考慮氣候變化數(shù)據(jù)、溫室氣體排放量以及預(yù)測(cè)模型、國(guó)家或地區(qū)的碳減排目標(biāo),同時(shí)兼顧不同行業(yè)的差異和經(jīng)濟(jì)實(shí)際情況等因素。

    第n年的階梯碳交易費(fèi)用CCO2,n的計(jì)算如式(16)所示:

    式中c表示碳交易基準(zhǔn)價(jià)格;d表示碳排放區(qū)間長(zhǎng)度; α 表示碳交易價(jià)格的增長(zhǎng)幅度。碳交易價(jià)格和交易碳排放量Ef,n的關(guān)系由圖2 所示。

    圖2 碳交易價(jià)格與交易碳排放量的關(guān)系圖Fig.2 Relationship between carbon trading price and traded carbon emission

    3 考慮最優(yōu)建設(shè)時(shí)序的園區(qū)綜合能源系統(tǒng)的多階段規(guī)劃模型

    對(duì)于規(guī)劃周期年限比較長(zhǎng)的PIES,在制定規(guī)劃時(shí)常常采用分期進(jìn)行的多階段規(guī)劃方法。為規(guī)避設(shè)備出現(xiàn)壽命到期退役的情況發(fā)生,文中將園區(qū)的規(guī)劃周期設(shè)置為各類設(shè)備的壽命最小值,并將該規(guī)劃周期分成多個(gè)階段分期規(guī)劃,對(duì)規(guī)劃期內(nèi)負(fù)荷的增長(zhǎng)速度進(jìn)行預(yù)測(cè),將投建年份和設(shè)備的投資容量設(shè)置為待尋優(yōu)量,然后求解優(yōu)化模型得出優(yōu)化結(jié)果。

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

    文中的PIES 模型以計(jì)及全壽命周期內(nèi)的投資成本、運(yùn)行成本、維護(hù)成本和碳交易成本之和再減去規(guī)劃期末的殘值費(fèi)用,并折算到投資初年的現(xiàn)值總費(fèi)用最小作為目標(biāo)函數(shù)。投資成本在各規(guī)劃階段的首年年初結(jié)算,殘值費(fèi)用在規(guī)劃期末的末年年末結(jié)算,運(yùn)行成本、維護(hù)成本和碳交易成本在規(guī)劃期內(nèi)的每年年末結(jié)算一次。文中模型的目標(biāo)函數(shù)為:

    式中C表示園區(qū)的全生命周期總費(fèi)用現(xiàn)值;K是園區(qū)規(guī)劃階段的數(shù)量;N表示園區(qū)規(guī)劃的周期年數(shù);k表示處于第k個(gè)規(guī)劃階段;n表示處于規(guī)劃的第n年;nk表示第k個(gè)規(guī)劃階段處于園區(qū)規(guī)劃周期中的第nk年; γ 表示年折算率;Cinv,k表示在第k個(gè)規(guī)劃階段PIES 系統(tǒng)的設(shè)備投資費(fèi)用;Cope,n、Cmain,n、CCO2,n表示第n年P(guān)IES 系統(tǒng)的運(yùn)行、維護(hù)與碳交易費(fèi)用;FRV表示規(guī)劃期末設(shè)備的殘值費(fèi)用。碳交易費(fèi)用的計(jì)算如式(16) 所示,投資成本、運(yùn)行成本、維護(hù)成本與殘值費(fèi)用的計(jì)算如下:

    1) 投資費(fèi)用。

    PIES 在第k個(gè)階段的設(shè)備投資費(fèi)用的計(jì)算如式(18) 所示:

    2) 運(yùn)行費(fèi)用。

    PIES 在第n年的運(yùn)行費(fèi)用計(jì)算如式(19) 所示:

    式中ce(t) 表示t時(shí)段的購(gòu)電價(jià)格;cg表示單位折算后的購(gòu)氣價(jià)格;Ggrid,n(t) 表示在第n年的第t時(shí)段園區(qū)向氣網(wǎng)的購(gòu)氣功率。

    3) 維護(hù)費(fèi)用。

    PIES 在第n年的維護(hù)費(fèi)用計(jì)算如式(20) 所示:

    4) 殘值費(fèi)用。

    在規(guī)劃期末,無論設(shè)備是否壽命到期,都存在一定的殘值費(fèi)用,所以需要考慮它們折舊后的殘值費(fèi)用,PIES 在規(guī)劃末期的殘值費(fèi)用計(jì)算如式(21) 所示:

    式中Cinv,j表示第j類設(shè)備的投資費(fèi)用; δj是第j類設(shè)備的凈殘值率;Tj是第j類設(shè)備在規(guī)劃期末的已使用年數(shù);Nj是第j類設(shè)備的壽命年限。

    3.2 約束條件

    1) 功率平衡約束。

    PIES 需要滿足電、熱、氣的功率平衡約束,如式(22) ~式(24) 所示:

    其中,PL,n(t) 、HL,n(t) 、GL,n(t) 分別表示第n年的第t時(shí)段電、熱、氣負(fù)荷功率;PHP,n(t) 與HHP,n(t) 分別表示第n年的第t時(shí)段熱泵的電輸入功率與熱輸出功率;GCHP,n(t) 與GGB,n(t) 分別表示第n年的第t時(shí)段熱電聯(lián)產(chǎn)和燃?xì)忮仩t的氣輸入功率;PchaES,n(t) 與PdisES,n(t)表示第n年的第t時(shí)段儲(chǔ)電裝置的充電、放電功率;與表示第n年的第t時(shí)段蓄熱裝置的蓄熱、放熱功率。

    2) 設(shè)備運(yùn)行約束。

    PIES 的設(shè)備運(yùn)行約束由式(1) ~式(9) 所示。

    3) 上級(jí)網(wǎng)絡(luò)交互功率約束。

    PIES 與上級(jí)電網(wǎng)、氣網(wǎng)的交互功率需要被限制在一定范圍以內(nèi),如式(25) ~式(26) 所示:

    4) “雙碳”約束。

    碳達(dá)峰約束如式(27) 所示:

    式中x表示碳達(dá)峰在規(guī)劃的第x年實(shí)現(xiàn),碳達(dá)峰指的是PIES 考慮碳配額后的交易碳排量最多。

    根據(jù)“雙碳”目標(biāo),在碳排放達(dá)到峰值之后,碳排放逐步回落。文中設(shè)置在碳達(dá)峰后,交易碳排放量每年下降v%及以上,故碳中和約束如式(28) 所示:

    經(jīng)常性團(tuán)建聚餐、每周分享會(huì)、按員工意愿與特長(zhǎng)分配工作、老父親般慈祥地對(duì)待員工失職借口。我還親手在公司搭健身角,鼓勵(lì)大家多鍛煉身體,保持良好狀態(tài)。

    式中y表示達(dá)到碳達(dá)峰之后的年份。

    文中構(gòu)建的是混合整數(shù)線性規(guī)劃模型,通過Matlab R 2021a 中的Cplex 求解器進(jìn)行求解,得到園區(qū)多階段規(guī)劃的優(yōu)化結(jié)果。

    4 算例研究

    4.1 算例數(shù)據(jù)

    文中的算例來源于我國(guó)某地區(qū)的電/熱/氣耦合的PIES,為真實(shí)地反映PIES 的運(yùn)行情況,將一年分成三個(gè)典型日,分別為夏季典型日、冬季典型日、過渡季典型日,三個(gè)典型日的電負(fù)荷、熱負(fù)荷、氣負(fù)荷需求曲線和PV 的最大輸出功率曲線如圖3 所示。

    圖3 典型日下負(fù)荷、光伏出力預(yù)測(cè)曲線Fig.3 Predicted load and PV output curves of typical days

    PIES 的規(guī)劃周期N設(shè)置為15 年,折算率γ 取為8%,設(shè)備殘值率δj取為0.06,儲(chǔ)能裝置的周期始、末荷能狀態(tài)SOC 設(shè)置為0.3。PIES 采用分時(shí)電價(jià)機(jī)制,園區(qū)向上級(jí)電網(wǎng)購(gòu)電的日內(nèi)分時(shí)電價(jià)如表1 所示,購(gòu)氣價(jià)格為0.28 元/( kW·h) 。參考全國(guó)碳交易市場(chǎng)的最新成交數(shù)據(jù),碳交易基準(zhǔn)價(jià)格定為56 元/t,碳交易相關(guān)的其他參數(shù)參考文獻(xiàn)[14]。

    表1 電價(jià)Tab.1 Electricity price

    各類待規(guī)劃設(shè)備的投資、維護(hù)及壽命等參數(shù)參考了相關(guān)文獻(xiàn)[24-27],并根據(jù)市場(chǎng)價(jià)格波動(dòng)適當(dāng)?shù)剡M(jìn)行調(diào)整,相關(guān)參數(shù)參考文獻(xiàn)[14]和文獻(xiàn)[24]。根據(jù)PIES 的負(fù)荷增長(zhǎng)特點(diǎn),即: 規(guī)劃初期處于快速增長(zhǎng)期,負(fù)荷的增長(zhǎng)率較快;規(guī)劃中期處于平穩(wěn)增長(zhǎng)期,負(fù)荷的增長(zhǎng)率有一定降低;規(guī)劃后期達(dá)到成熟期,負(fù)荷增速進(jìn)一步降低;給出PIES 的負(fù)荷信息如表2 所示。運(yùn)行算例的計(jì)算機(jī)CPU 是Core i7-1195G7,內(nèi)存16G。

    表2 PIES 的負(fù)荷信息Tab.2 Load information of PIES

    4.2 各場(chǎng)景仿真結(jié)果分析

    對(duì)于文中建立的考慮最優(yōu)建設(shè)時(shí)序和階梯碳交易的規(guī)劃模型,設(shè)置了如下4 種場(chǎng)景進(jìn)行對(duì)比。其中,設(shè)置PIES 規(guī)劃的第8 年達(dá)峰; 在碳達(dá)峰后,交易碳排放量每年下降8%及以上。

    1) 場(chǎng)景1:單階段規(guī)劃,考慮階梯碳交易;

    3) 場(chǎng)景3:多階段規(guī)劃,規(guī)劃段數(shù)為3,在規(guī)劃期的第1、6、11 年投建設(shè)備,考慮階梯碳交易;

    4) 場(chǎng)景4:多階段規(guī)劃,規(guī)劃段數(shù)為3,投建年份待尋優(yōu),考慮階梯碳交易。

    各類場(chǎng)景下的投建年份及設(shè)備的配置容量見表3,各類場(chǎng)景折算到規(guī)劃初年的現(xiàn)值費(fèi)用及碳排放量見表4。

    表3 各類場(chǎng)景下的投建年份及設(shè)備的配置容量Tab.3 Construction year and equipment configuration capacity under different scenarios

    表4 各類場(chǎng)景折算到規(guī)劃初年的現(xiàn)值費(fèi)用及碳排放量Tab.4 Present value costs and carbon emissions of scenarios converted to the initial planning year

    對(duì)比場(chǎng)景1 和場(chǎng)景4 可以發(fā)現(xiàn),場(chǎng)景4 的PV、HP、CHP 這些設(shè)備的投建總?cè)萘勘葓?chǎng)景1 高,場(chǎng)景4 的GB設(shè)備的總?cè)萘勘葓?chǎng)景1 低,這是多階段規(guī)劃在規(guī)劃后期階段投入的設(shè)備的殘值較高的原因。場(chǎng)景1 與場(chǎng)景4 的投資費(fèi)用與殘值費(fèi)用之差比較接近,但是運(yùn)行費(fèi)用場(chǎng)景1 比場(chǎng)景4 高3.06%,這是因?yàn)樵谝?guī)劃前期階段,場(chǎng)景1 投建的設(shè)備冗余嚴(yán)重,雖然投資費(fèi)用低,但減去殘值費(fèi)用與場(chǎng)景4 接近;但在規(guī)劃后期階段,場(chǎng)景1 的設(shè)備投建容量相對(duì)場(chǎng)景4 稍低,更多地依賴于從外部購(gòu)氣、購(gòu)電,所以運(yùn)行費(fèi)用更高。場(chǎng)景1 的碳排放量比場(chǎng)景4 高27.11%,這是由于多階段規(guī)劃中設(shè)備規(guī)劃更加靈活,在規(guī)劃后期可以選擇投建更多的PV、CHP 這類碳排放較少的設(shè)備,所以多階段規(guī)劃中總碳排放量大大降低。這說明在進(jìn)行多階段規(guī)劃后,在滿足各類負(fù)荷需求的基礎(chǔ)上,可以減少單階段規(guī)劃因超前投建所導(dǎo)致的設(shè)備冗余和折舊損耗,設(shè)備投建更加靈活,也可以減少總碳排放,對(duì)規(guī)劃期內(nèi)的全壽命周期成本和低碳運(yùn)行都有較大的優(yōu)化效果。

    將場(chǎng)景2 和場(chǎng)景4 進(jìn)行對(duì)比可以得到,場(chǎng)景4 的PV、CHP 設(shè)備的投資總?cè)萘勘葓?chǎng)景2 稍高,場(chǎng)景4 的HP 設(shè)備的投資總?cè)萘勘葓?chǎng)景2 稍低,這是因?yàn)閳?chǎng)景2的碳交易以固定的基準(zhǔn)價(jià)格進(jìn)行,相當(dāng)于變相地降低了碳交易價(jià)格,所以場(chǎng)景2 的規(guī)劃傾向于選擇相對(duì)便宜但碳排放量高的HP 設(shè)備,場(chǎng)景2 的總碳排放量遠(yuǎn)大于場(chǎng)景4。然而,場(chǎng)景2 的全壽命周期費(fèi)用比場(chǎng)景4 略低,這是因?yàn)閳?chǎng)景2 的碳交易費(fèi)用的單價(jià)降低了。綜合來看,場(chǎng)景4 選擇了階梯碳交易雖然會(huì)導(dǎo)致總費(fèi)用稍微增加了0.03%,但是碳排放量相對(duì)場(chǎng)景2 降低了12.87%,這說明了階梯碳交易機(jī)制比以固定的價(jià)格進(jìn)行碳交易在實(shí)現(xiàn)低碳性的方面更有優(yōu)勢(shì)。

    對(duì)比場(chǎng)景3 和場(chǎng)景4 可以發(fā)現(xiàn),場(chǎng)景4 僅僅是將PIES 內(nèi)設(shè)備的投資年份設(shè)置為變量,通過求解得出設(shè)備的最佳投建年份,其他均與場(chǎng)景3 相同。場(chǎng)景4 的投建費(fèi)用與全壽命周期成本比場(chǎng)景3 略低,且總碳排放量比場(chǎng)景3 下降了1.02%。由此可知,將投建年份設(shè)為變量,各類規(guī)劃設(shè)備的配置會(huì)更加靈活,全壽命周期成本和低碳性可以進(jìn)一步得到優(yōu)化。

    場(chǎng)景1 - 場(chǎng)景4 的求解時(shí)間分別是10. 498 s、31.145 s、11.107 s、36.279 s,可見將規(guī)劃年份設(shè)置為待尋優(yōu)量后,求解時(shí)間明顯變長(zhǎng)。

    選取PIES 在場(chǎng)景4 下的規(guī)劃期最后一年的過渡季典型日的電、熱、氣功率平衡進(jìn)行分析,各設(shè)備出力及負(fù)荷結(jié)果如圖4 所示??梢园l(fā)現(xiàn),儲(chǔ)電裝置在電價(jià)低谷時(shí)段充電,在電價(jià)高峰時(shí)放電,符合“低儲(chǔ)高發(fā)”原則;由于電熱相互耦合,蓄熱裝置也是在電價(jià)低谷時(shí)段蓄熱,在電價(jià)高峰時(shí)放熱。PIES 中的各類能源設(shè)備協(xié)同調(diào)度,在滿足用戶各類負(fù)荷需求的同時(shí)盡量降低自身的運(yùn)行成本。

    圖4 規(guī)劃期末年過渡季典型日的功率平衡Fig.4 Supply and demand balance on typical day of transition season at the ending year of planning period

    4.3 碳交易基準(zhǔn)價(jià)格的影響分析

    基于場(chǎng)景4,探究分析碳交易基準(zhǔn)價(jià)格對(duì)園區(qū)規(guī)劃的影響,改變碳交易的基準(zhǔn)價(jià)格,得到對(duì)應(yīng)的總碳排放量和碳交易費(fèi)用現(xiàn)值,結(jié)果如圖5 所示。由圖5 可以得到,伴隨著碳交易基準(zhǔn)價(jià)格越來越高,總碳排放量在逐步地下降,碳交易費(fèi)用現(xiàn)值總體來說是先升高后降低的趨勢(shì)。當(dāng)碳交易費(fèi)用較低時(shí),園區(qū)出于考慮經(jīng)濟(jì)性的目標(biāo),更多地選擇外購(gòu)價(jià)格較低的電和投建價(jià)格相對(duì)較低但碳排放量較多的HP 設(shè)備,由于碳交易價(jià)格的升高,碳交易費(fèi)用對(duì)園區(qū)規(guī)劃的影響權(quán)重增加,PIES開始逐步地提高光伏、CHP 這類清潔設(shè)備的出力,改變園區(qū)的用能方式,從而降低系統(tǒng)的碳排放量。

    圖5 碳排放量、碳交易總成本、碳交易基準(zhǔn)價(jià)格的關(guān)系Fig.5 Relations between carbon emission,total cost of carbon trading and benchmark price of carbon trading

    4.4 PIES 規(guī)劃的投建次數(shù)的影響

    考慮階梯碳交易,將PIES 的投建年份設(shè)置為待尋優(yōu)量,同時(shí)將園區(qū)的投建次數(shù)K設(shè)置為變量。設(shè)置一個(gè)0-1 變量Xi用以表示PIES 在第i年是否投建了新設(shè)備,Xi為0 表示該年沒有投建任何設(shè)備,Xi為1 表示投建了新設(shè)備,即= K。

    研究投建次數(shù)K對(duì)全壽命周期費(fèi)用現(xiàn)值的影響,結(jié)果如表5 所示。根據(jù)表5 可知,隨著投建次數(shù)的增加,全壽命周期費(fèi)用現(xiàn)值有所下降。這是因?yàn)殡S著投建次數(shù)的增加,待規(guī)劃設(shè)備投建可以更加靈活,超前投建導(dǎo)致的冗余更少,設(shè)備的利用率更高,全壽命周期費(fèi)用現(xiàn)值越少。但是投建次數(shù)越多,投建所需的人工費(fèi)用越多,文中未將這部分費(fèi)用計(jì)入目標(biāo)函數(shù),所以在園區(qū)規(guī)劃時(shí)還應(yīng)根據(jù)實(shí)際情況選擇合適的投建次數(shù)。

    表5 投建次數(shù)與全壽命周期費(fèi)用現(xiàn)值的關(guān)系Tab.5 Relations between the number of construction phases and the present value of life cycle cost

    4.5 改變“雙碳”約束的影響分析

    在場(chǎng)景4 的基礎(chǔ)上設(shè)置在碳達(dá)峰后交易碳排放量每年下降8%及以上,同時(shí)改變園區(qū)碳達(dá)峰的年份,其他條件都保持不變。將碳達(dá)峰的年份分別設(shè)置為園區(qū)規(guī)劃的第5 年、第8 年、第11 年,不同的碳達(dá)峰年份的規(guī)劃結(jié)果的如表6 所示。根據(jù)表6 可知,隨著碳達(dá)峰年份的提前,園區(qū)的CHP 配置總?cè)萘繒?huì)增加,同時(shí)全壽命周期費(fèi)用增加,總碳排放量減少。這是因?yàn)樘歼_(dá)峰年份的提前,對(duì)碳排放的約束更高,總碳排放量會(huì)有所降低,園區(qū)傾向于選擇CHP 這類清潔設(shè)備,同時(shí)全壽命周期費(fèi)用也會(huì)有所升高。

    表6 不同的碳達(dá)峰年份的規(guī)劃結(jié)果Tab.6 Planning results for different carbon peak years

    在場(chǎng)景4 的基礎(chǔ)上設(shè)置PIES 規(guī)劃的第8 年達(dá)峰,同時(shí)改變園區(qū)在碳達(dá)峰后交易碳排放量每年下降的速度,其他條件都保持不變。在第8 年園區(qū)碳達(dá)峰后,交易碳排放量分別設(shè)置為每年至少下降4%、8%、12%,碳達(dá)峰后交易碳排放量不同下降速度所對(duì)應(yīng)的規(guī)劃結(jié)果如表7 所示。根據(jù)表7 可知,隨著碳達(dá)峰后交易碳排放量每年下降比例的增加,對(duì)碳排放的約束也會(huì)更嚴(yán)格,碳排放量會(huì)有所下降,園區(qū)也會(huì)傾向于選擇CHP這類清潔設(shè)備,但是會(huì)導(dǎo)致園區(qū)總費(fèi)用有一定的增加。

    表7 碳達(dá)峰后交易碳排放量不同下降速度的規(guī)劃結(jié)果Tab.7 Planning results of different decline rates of carbon emissions after reaching carbon peak

    5 結(jié)束語

    為了更好地實(shí)現(xiàn)國(guó)家“雙碳”戰(zhàn)略任務(wù),文中提出了考慮階梯碳交易和最優(yōu)建設(shè)時(shí)序的PIES 多階段規(guī)劃,以全壽命周期總費(fèi)用折算到投資初年的現(xiàn)值最小作為優(yōu)化目標(biāo),求解得到設(shè)備投建的最優(yōu)年份和最優(yōu)設(shè)備容量。通過仿真實(shí)驗(yàn),得到以下結(jié)論:

    1) 所提出的考慮最優(yōu)建設(shè)時(shí)序的PIES 多階段規(guī)劃方法可以使得設(shè)備投建更加靈活,設(shè)備冗余更少,對(duì)園區(qū)整體的經(jīng)濟(jì)性與低碳性有一定優(yōu)化效果;

    2) 基于階梯碳交易機(jī)制,同時(shí)與我國(guó)“雙碳”戰(zhàn)略目標(biāo)相結(jié)合,在模型中考慮碳達(dá)峰的年份以及碳達(dá)峰后碳排放量逐年下降的約束; 研究碳價(jià)對(duì)碳排放量和碳交易總成本的影響,從而分析碳交易價(jià)格對(duì)園區(qū)用能方式的影響;

    3) 考慮投建次數(shù)對(duì)園區(qū)規(guī)劃的影響,可以發(fā)現(xiàn)投建次數(shù)越多,全壽命周期成本越少,從而分析得到: 隨著投建次數(shù)增多,設(shè)備可以投建地更加靈活,設(shè)備的利用率更高;

    4) 研究“雙碳”約束對(duì)園區(qū)規(guī)劃的影響。隨著碳達(dá)峰年份的提前或碳達(dá)峰后交易碳排放量下降速度的增加,對(duì)園區(qū)碳約束更加嚴(yán)格,配置的CHP 清潔設(shè)備有所增加,總碳排放量下降,同時(shí)也會(huì)導(dǎo)致全壽命周期費(fèi)用有所增加。

    未來的研究可以深入探討光伏出力不確定性對(duì)規(guī)劃結(jié)果的影響,引入分布魯棒模型刻畫光伏出力不確定性,并構(gòu)建系統(tǒng)容量規(guī)劃的分布魯棒優(yōu)化模型,降低優(yōu)化結(jié)果的保守性。

    猜你喜歡
    碳達(dá)峰階梯排放量
    碳中和·碳達(dá)峰
    關(guān)于“碳達(dá)峰”和“碳中和”——2030碳達(dá)峰與2060碳中和
    新疆鋼鐵(2021年1期)2021-10-14 08:45:32
    實(shí)現(xiàn)“碳達(dá)峰碳中和”應(yīng)當(dāng)采用何種立法思路
    碳達(dá)峰碳中和要求下的設(shè)計(jì)思考
    天然氣輸配系統(tǒng)甲烷排放量化方法
    煤氣與熱力(2021年6期)2021-07-28 07:21:40
    黑龍江省碳排放量影響因素研究
    爬階梯
    時(shí)光階梯
    幸福(2016年9期)2016-12-01 03:08:50
    有趣的階梯
    全國(guó)機(jī)動(dòng)車污染物排放量
    ——《2013年中國(guó)機(jī)動(dòng)車污染防治年報(bào)》(第Ⅱ部分)
    成人av在线播放网站| 欧美成人性av电影在线观看| 99国产极品粉嫩在线观看| 色综合婷婷激情| 国产精品嫩草影院av在线观看 | av欧美777| 欧美日韩国产亚洲二区| 亚洲第一欧美日韩一区二区三区| а√天堂www在线а√下载| 一a级毛片在线观看| 亚洲专区国产一区二区| 欧美另类亚洲清纯唯美| 亚洲五月天丁香| 亚洲18禁久久av| 亚洲中文字幕一区二区三区有码在线看| 人人妻人人看人人澡| 亚洲一区高清亚洲精品| 嫩草影视91久久| 亚洲精品一区av在线观看| 在线天堂最新版资源| 1000部很黄的大片| 欧美黄色淫秽网站| 亚洲欧美激情综合另类| 女人被狂操c到高潮| 男插女下体视频免费在线播放| 噜噜噜噜噜久久久久久91| www.色视频.com| 久久欧美精品欧美久久欧美| 国产黄片美女视频| 国产激情偷乱视频一区二区| 美女高潮的动态| 亚洲五月天丁香| 我要搜黄色片| 国产国拍精品亚洲av在线观看 | 亚洲国产欧洲综合997久久,| 丁香欧美五月| 国产蜜桃级精品一区二区三区| 日韩欧美国产一区二区入口| 婷婷亚洲欧美| 午夜日韩欧美国产| 在线播放国产精品三级| 美女高潮的动态| 久99久视频精品免费| 日本免费a在线| 欧美性感艳星| 国产一区二区在线观看日韩 | 国产精品av视频在线免费观看| 久久久成人免费电影| 日韩欧美精品v在线| 日本黄大片高清| 老熟妇乱子伦视频在线观看| 又黄又粗又硬又大视频| 亚洲欧美日韩高清在线视频| 成熟少妇高潮喷水视频| 三级毛片av免费| 最近最新免费中文字幕在线| 九九热线精品视视频播放| 久久久成人免费电影| 人人妻人人澡欧美一区二区| 久久精品91蜜桃| 亚洲av五月六月丁香网| 欧美+日韩+精品| 精品国内亚洲2022精品成人| 亚洲国产精品久久男人天堂| 18禁黄网站禁片午夜丰满| 午夜免费激情av| 天天躁日日操中文字幕| 俄罗斯特黄特色一大片| 国产av不卡久久| 99久久精品热视频| 国产真实伦视频高清在线观看 | 丰满乱子伦码专区| 日本 av在线| 丰满的人妻完整版| 久久久久久久午夜电影| 99热只有精品国产| 深夜精品福利| 无遮挡黄片免费观看| 亚洲成人久久性| 国产黄片美女视频| 欧美日韩国产亚洲二区| 国产激情欧美一区二区| 天堂网av新在线| 九色国产91popny在线| 国产精品久久久久久久久免 | 欧美在线黄色| 久久久久亚洲av毛片大全| av天堂中文字幕网| 亚洲中文日韩欧美视频| 国产高清三级在线| 色尼玛亚洲综合影院| 午夜福利成人在线免费观看| 国产精品久久久久久久久免 | 女生性感内裤真人,穿戴方法视频| 精品乱码久久久久久99久播| 色综合亚洲欧美另类图片| 国产成人系列免费观看| 日本与韩国留学比较| 中文字幕av成人在线电影| 亚洲欧美日韩无卡精品| 老熟妇仑乱视频hdxx| 日本黄色视频三级网站网址| 老司机在亚洲福利影院| 日本熟妇午夜| 午夜日韩欧美国产| 精品久久久久久久人妻蜜臀av| 亚洲欧美日韩高清在线视频| 亚洲天堂国产精品一区在线| 欧美一级毛片孕妇| 精品无人区乱码1区二区| 99在线视频只有这里精品首页| 国产成人影院久久av| 观看美女的网站| 最近视频中文字幕2019在线8| 国产欧美日韩精品亚洲av| 国产91精品成人一区二区三区| 日韩 欧美 亚洲 中文字幕| 亚洲真实伦在线观看| 国产精品久久久久久亚洲av鲁大| 国产高潮美女av| 一进一出抽搐gif免费好疼| 日本黄大片高清| 三级国产精品欧美在线观看| 日本成人三级电影网站| 国产精品99久久久久久久久| 欧美高清成人免费视频www| 日韩欧美在线乱码| 中文字幕高清在线视频| av专区在线播放| av专区在线播放| 真人做人爱边吃奶动态| 窝窝影院91人妻| 国产在视频线在精品| www.熟女人妻精品国产| 国产高清videossex| 黄色日韩在线| 国产毛片a区久久久久| h日本视频在线播放| 免费av不卡在线播放| 久久人人精品亚洲av| 亚洲国产欧洲综合997久久,| 精品免费久久久久久久清纯| 中文字幕熟女人妻在线| 女生性感内裤真人,穿戴方法视频| 国产高清视频在线观看网站| 久久久久久久久久黄片| 精品国产美女av久久久久小说| 国内揄拍国产精品人妻在线| 精品国产三级普通话版| 国产91精品成人一区二区三区| 久久中文看片网| 亚洲av成人av| 国内精品一区二区在线观看| 国产精品99久久久久久久久| 在线视频色国产色| or卡值多少钱| 精品无人区乱码1区二区| 99热精品在线国产| 在线观看66精品国产| 欧美av亚洲av综合av国产av| 九九热线精品视视频播放| 叶爱在线成人免费视频播放| 超碰av人人做人人爽久久 | 久9热在线精品视频| 91av网一区二区| 国产一区在线观看成人免费| 国产真实伦视频高清在线观看 | 成人18禁在线播放| 99精品在免费线老司机午夜| 亚洲av二区三区四区| 国产精品1区2区在线观看.| 国产精品久久久久久久电影 | 久久亚洲真实| 禁无遮挡网站| x7x7x7水蜜桃| 桃色一区二区三区在线观看| 免费看日本二区| 亚洲国产中文字幕在线视频| 久久国产精品人妻蜜桃| 欧美激情久久久久久爽电影| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放| 欧美xxxx黑人xx丫x性爽| 黄片大片在线免费观看| 性色avwww在线观看| 久久久色成人| 熟女人妻精品中文字幕| aaaaa片日本免费| 国产成人啪精品午夜网站| 尤物成人国产欧美一区二区三区| 国产探花在线观看一区二区| 91九色精品人成在线观看| 亚洲国产精品合色在线| 精品国产亚洲在线| 神马国产精品三级电影在线观看| 亚洲成人精品中文字幕电影| 一进一出好大好爽视频| 国产av不卡久久| 日本a在线网址| 男插女下体视频免费在线播放| 操出白浆在线播放| 欧美丝袜亚洲另类 | 日韩成人在线观看一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 国产高清视频在线观看网站| 亚洲狠狠婷婷综合久久图片| 岛国在线免费视频观看| 久久久久九九精品影院| 18禁黄网站禁片午夜丰满| 成人亚洲精品av一区二区| 观看免费一级毛片| 久久久久久久久中文| 99久久九九国产精品国产免费| 亚洲人成网站在线播放欧美日韩| 每晚都被弄得嗷嗷叫到高潮| av天堂在线播放| 国产精品av视频在线免费观看| 长腿黑丝高跟| 一夜夜www| 精品欧美国产一区二区三| e午夜精品久久久久久久| 成人18禁在线播放| 免费无遮挡裸体视频| 精品国产超薄肉色丝袜足j| 久久久久久九九精品二区国产| 亚洲自拍偷在线| 宅男免费午夜| 久久精品影院6| 一夜夜www| 久久久国产成人免费| 最新在线观看一区二区三区| 亚洲久久久久久中文字幕| 国产在视频线在精品| 少妇的逼水好多| 一级黄色大片毛片| 久久欧美精品欧美久久欧美| 男人舔女人下体高潮全视频| 成年版毛片免费区| 午夜福利高清视频| 国产一区在线观看成人免费| 免费大片18禁| 久久精品91无色码中文字幕| 操出白浆在线播放| 国产成人啪精品午夜网站| 成人欧美大片| 精品不卡国产一区二区三区| 男人舔奶头视频| 久久亚洲精品不卡| 深夜精品福利| av视频在线观看入口| 麻豆国产av国片精品| 国产精品亚洲av一区麻豆| 亚洲天堂国产精品一区在线| 欧美日韩国产亚洲二区| 精品国产亚洲在线| 亚洲成人中文字幕在线播放| 国产精品av视频在线免费观看| 国产真实伦视频高清在线观看 | 国产午夜精品久久久久久一区二区三区 | 国产精品一区二区免费欧美| 免费一级毛片在线播放高清视频| 757午夜福利合集在线观看| 老熟妇仑乱视频hdxx| 美女黄网站色视频| 岛国在线免费视频观看| 美女 人体艺术 gogo| 亚洲精品乱码久久久v下载方式 | 国产一区二区在线av高清观看| 男女视频在线观看网站免费| 首页视频小说图片口味搜索| 欧美一级a爱片免费观看看| 色老头精品视频在线观看| 精品人妻偷拍中文字幕| 久久精品国产清高在天天线| 一个人免费在线观看电影| 国产真实乱freesex| 亚洲第一电影网av| 国产午夜精品论理片| 我的老师免费观看完整版| 国产主播在线观看一区二区| 欧美性猛交黑人性爽| 亚洲成av人片免费观看| 成人亚洲精品av一区二区| 中文字幕高清在线视频| 亚洲中文字幕一区二区三区有码在线看| 国内揄拍国产精品人妻在线| 偷拍熟女少妇极品色| 欧美日韩一级在线毛片| 中文字幕人妻熟人妻熟丝袜美 | 乱人视频在线观看| 精品熟女少妇八av免费久了| 国产 一区 欧美 日韩| 淫妇啪啪啪对白视频| 亚洲人成网站高清观看| 成人特级黄色片久久久久久久| 五月伊人婷婷丁香| 亚洲国产欧美人成| 久久性视频一级片| 亚洲av不卡在线观看| 免费av观看视频| 欧美日韩黄片免| 热99在线观看视频| 男女之事视频高清在线观看| 18禁黄网站禁片午夜丰满| 国产色爽女视频免费观看| 亚洲色图av天堂| 2021天堂中文幕一二区在线观| 国产一区二区在线av高清观看| 国产精华一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 搡老熟女国产l中国老女人| 在线观看av片永久免费下载| 99热只有精品国产| 中亚洲国语对白在线视频| 国产69精品久久久久777片| 麻豆一二三区av精品| 欧美一区二区精品小视频在线| 色吧在线观看| 亚洲精品久久国产高清桃花| 法律面前人人平等表现在哪些方面| 老司机深夜福利视频在线观看| 1000部很黄的大片| 亚洲精品色激情综合| 一个人观看的视频www高清免费观看| 亚洲精品456在线播放app | 变态另类丝袜制服| 青草久久国产| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 丝袜美腿在线中文| 啦啦啦观看免费观看视频高清| ponron亚洲| 国产国拍精品亚洲av在线观看 | 少妇的丰满在线观看| 国内精品久久久久精免费| 国产黄片美女视频| 精品国产三级普通话版| 国产精品 国内视频| 国产精品久久电影中文字幕| 成人国产综合亚洲| 最新中文字幕久久久久| 成人午夜高清在线视频| 久久久国产成人精品二区| 国产亚洲精品综合一区在线观看| 国产爱豆传媒在线观看| 欧美成狂野欧美在线观看| 亚洲国产日韩欧美精品在线观看 | 国产中年淑女户外野战色| 可以在线观看毛片的网站| 久久国产乱子伦精品免费另类| 亚洲av第一区精品v没综合| 精品99又大又爽又粗少妇毛片 | 精品久久久久久久末码| 成人亚洲精品av一区二区| 特级一级黄色大片| 人人妻人人澡欧美一区二区| 听说在线观看完整版免费高清| 亚洲欧美激情综合另类| 成人永久免费在线观看视频| 嫩草影视91久久| 欧美日韩中文字幕国产精品一区二区三区| 日韩亚洲欧美综合| 怎么达到女性高潮| 精品国产超薄肉色丝袜足j| 18禁国产床啪视频网站| 午夜福利视频1000在线观看| 男女那种视频在线观看| 大型黄色视频在线免费观看| 亚洲无线在线观看| 成年免费大片在线观看| 欧美+日韩+精品| 国产一区二区在线观看日韩 | 国产欧美日韩一区二区精品| 啦啦啦韩国在线观看视频| ponron亚洲| or卡值多少钱| 欧美黄色淫秽网站| 在线十欧美十亚洲十日本专区| 久久亚洲精品不卡| 午夜免费激情av| 看免费av毛片| 日韩欧美 国产精品| 国产三级黄色录像| 日本 欧美在线| 国产97色在线日韩免费| 久久久久久久久中文| 黄色日韩在线| 久久久成人免费电影| 午夜免费激情av| 国产伦精品一区二区三区四那| 国产精品影院久久| 成人特级av手机在线观看| 老司机午夜十八禁免费视频| 欧美日本视频| 美女大奶头视频| 免费搜索国产男女视频| 最近最新中文字幕大全电影3| 人妻夜夜爽99麻豆av| 色综合站精品国产| 久久久国产成人精品二区| 欧美日韩精品网址| 91麻豆精品激情在线观看国产| 欧美成狂野欧美在线观看| 免费av观看视频| 男女下面进入的视频免费午夜| 久久九九热精品免费| 99精品在免费线老司机午夜| 亚洲精品在线观看二区| 丝袜美腿在线中文| 两个人看的免费小视频| 午夜福利高清视频| 亚洲av第一区精品v没综合| 欧美国产日韩亚洲一区| 欧美日韩黄片免| 嫩草影院精品99| 色av中文字幕| 亚洲人成网站在线播放欧美日韩| 亚洲内射少妇av| 国产单亲对白刺激| 女同久久另类99精品国产91| 久久香蕉国产精品| 精品日产1卡2卡| 精品不卡国产一区二区三区| 国产精品三级大全| 亚洲av不卡在线观看| 精品人妻1区二区| 亚洲欧美日韩高清专用| 欧美+日韩+精品| 男女做爰动态图高潮gif福利片| 青草久久国产| 18+在线观看网站| 中文字幕人成人乱码亚洲影| 青草久久国产| 成人永久免费在线观看视频| 国产精品一及| 国产成人av教育| 网址你懂的国产日韩在线| 亚洲18禁久久av| 国产97色在线日韩免费| 麻豆成人av在线观看| 久久久国产精品麻豆| 国产精品久久电影中文字幕| 亚洲一区二区三区色噜噜| or卡值多少钱| 级片在线观看| 午夜福利在线在线| 99久久无色码亚洲精品果冻| 色综合欧美亚洲国产小说| 精品国产超薄肉色丝袜足j| 99国产综合亚洲精品| 丰满人妻熟妇乱又伦精品不卡| 精品电影一区二区在线| 欧美zozozo另类| 成人欧美大片| 一区福利在线观看| 久9热在线精品视频| 国产精品乱码一区二三区的特点| 欧美黄色片欧美黄色片| 美女高潮的动态| 动漫黄色视频在线观看| 久久久久久久午夜电影| 成人亚洲精品av一区二区| 免费在线观看影片大全网站| 日韩欧美免费精品| 国产又黄又爽又无遮挡在线| 午夜免费成人在线视频| 国模一区二区三区四区视频| 窝窝影院91人妻| 久久亚洲精品不卡| 亚洲成人久久爱视频| 激情在线观看视频在线高清| 亚洲国产欧美人成| 在线播放国产精品三级| 亚洲乱码一区二区免费版| 亚洲av成人精品一区久久| 91久久精品电影网| 日本免费a在线| 国产在线精品亚洲第一网站| 午夜福利18| 午夜福利在线观看吧| 99久久九九国产精品国产免费| 99在线人妻在线中文字幕| 熟女少妇亚洲综合色aaa.| 免费在线观看影片大全网站| 91久久精品国产一区二区成人 | 国产精品一区二区三区四区免费观看 | 亚洲无线观看免费| 精品乱码久久久久久99久播| 午夜福利成人在线免费观看| 亚洲av成人不卡在线观看播放网| 欧美成人一区二区免费高清观看| 好看av亚洲va欧美ⅴa在| 亚洲精品色激情综合| 欧美极品一区二区三区四区| 伊人久久大香线蕉亚洲五| www.www免费av| 成人高潮视频无遮挡免费网站| 久久久精品欧美日韩精品| 中文字幕久久专区| 黄色女人牲交| 国产97色在线日韩免费| 成人av一区二区三区在线看| 欧美成人a在线观看| 国产伦在线观看视频一区| 国产国拍精品亚洲av在线观看 | 日日摸夜夜添夜夜添小说| 热99在线观看视频| 国产97色在线日韩免费| 91久久精品电影网| 嫁个100分男人电影在线观看| 伊人久久精品亚洲午夜| 亚洲国产欧美人成| 国产aⅴ精品一区二区三区波| 午夜影院日韩av| 狂野欧美白嫩少妇大欣赏| 精品国产亚洲在线| 亚洲黑人精品在线| 伊人久久精品亚洲午夜| 综合色av麻豆| 99热精品在线国产| 免费av毛片视频| 亚洲欧美精品综合久久99| 又粗又爽又猛毛片免费看| 精品欧美国产一区二区三| 九九在线视频观看精品| 岛国在线免费视频观看| 搡老熟女国产l中国老女人| 国产国拍精品亚洲av在线观看 | 一本综合久久免费| 亚洲精品粉嫩美女一区| 国产高清三级在线| 亚洲精品成人久久久久久| av片东京热男人的天堂| 动漫黄色视频在线观看| 亚洲欧美一区二区三区黑人| 18+在线观看网站| 天堂动漫精品| 亚洲精品乱码久久久v下载方式 | 亚洲av免费高清在线观看| 激情在线观看视频在线高清| 九九热线精品视视频播放| 两个人看的免费小视频| www.色视频.com| 亚洲天堂国产精品一区在线| 高清日韩中文字幕在线| 亚洲国产中文字幕在线视频| 亚洲美女黄片视频| 51国产日韩欧美| 久久精品91蜜桃| 精品久久久久久久末码| 日韩国内少妇激情av| 99国产综合亚洲精品| 人人妻人人看人人澡| 亚洲精品国产精品久久久不卡| 欧美一级a爱片免费观看看| 欧美午夜高清在线| 欧洲精品卡2卡3卡4卡5卡区| 91av网一区二区| 女人被狂操c到高潮| 此物有八面人人有两片| 久9热在线精品视频| 精品国产超薄肉色丝袜足j| 在线观看午夜福利视频| 成年版毛片免费区| 色播亚洲综合网| 亚洲最大成人中文| а√天堂www在线а√下载| 国产熟女xx| 欧美性猛交黑人性爽| 级片在线观看| 中文字幕人成人乱码亚洲影| 国产精品美女特级片免费视频播放器| 成人高潮视频无遮挡免费网站| АⅤ资源中文在线天堂| 亚洲国产中文字幕在线视频| 最好的美女福利视频网| 国产精华一区二区三区| 国产精品99久久久久久久久| 国产高清视频在线观看网站| 亚洲国产精品久久男人天堂| 亚洲国产精品999在线| 亚洲成a人片在线一区二区| 久久久久亚洲av毛片大全| 亚洲人成伊人成综合网2020| 国产男靠女视频免费网站| 99久久精品国产亚洲精品| 午夜福利18| 精品日产1卡2卡| 男女做爰动态图高潮gif福利片| 99热这里只有精品一区| 久久婷婷人人爽人人干人人爱| 久久久久精品国产欧美久久久| 黄色女人牲交| 两个人看的免费小视频| 天堂动漫精品| 免费一级毛片在线播放高清视频| a级毛片a级免费在线| 欧美+日韩+精品| 欧美成人性av电影在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精华国产精华精| 亚洲人成网站在线播放欧美日韩| 久久性视频一级片| 成人av一区二区三区在线看| 国产成人av教育| 午夜福利18| 99热这里只有精品一区| 国产私拍福利视频在线观看| 国内精品久久久久精免费| 久久香蕉国产精品| 中文资源天堂在线| 国产黄色小视频在线观看| 手机成人av网站| 日本一二三区视频观看| 两个人的视频大全免费| 无人区码免费观看不卡|