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

    基于熱需求響應(yīng)的綜合能源系統(tǒng)可靠供能能力評(píng)估

    2022-05-06 12:00:08張治魏振華鄧子琦
    科學(xué)技術(shù)與工程 2022年11期
    關(guān)鍵詞:區(qū)間峰值可靠性

    張治, 魏振華*, 鄧子琦

    (1.國(guó)核電力規(guī)劃設(shè)計(jì)研究院有限公司, 北京 100095; 2.華北電力大學(xué)控制與計(jì)算機(jī)工程學(xué)院, 北京 102206)

    隨著社會(huì)的持續(xù)發(fā)展,能源結(jié)構(gòu)正在逐步向清潔發(fā)電轉(zhuǎn)變,能源生產(chǎn)和消費(fèi)模式正在發(fā)生重大轉(zhuǎn)變。綜合能源系統(tǒng)(integrated energy system, IES)是一種新型的混合動(dòng)力系統(tǒng)[1],通過先進(jìn)的能源轉(zhuǎn)換和能源傳輸儲(chǔ)存技術(shù)整合多種不同能源形式(如電能、天然氣、熱能等),實(shí)現(xiàn)多個(gè)能源子系統(tǒng)的協(xié)調(diào)運(yùn)行,梯級(jí)利用能源,從而提升能源效率[2]。在新一輪工業(yè)革命的條件下,綜合能源系統(tǒng)是滿足多元供應(yīng)體系的具體實(shí)現(xiàn)方式,已經(jīng)成為世界各國(guó)能源轉(zhuǎn)型爭(zhēng)相發(fā)展的重點(diǎn)[3]。

    目前,已有許多學(xué)者在綜合能源系統(tǒng)的規(guī)劃分析、優(yōu)化調(diào)度上方面進(jìn)行了相關(guān)研究。文獻(xiàn)[4]從站-網(wǎng)整體規(guī)劃的角度出發(fā),對(duì)能源站數(shù)量、位置和設(shè)備容量配置,供能網(wǎng)絡(luò)布局進(jìn)行統(tǒng)一規(guī)劃研究。文獻(xiàn)[5]在綜合能源配置優(yōu)化問題中,引入可靠性約束,形成兼顧經(jīng)濟(jì)性與可靠性的綜合能源系統(tǒng)優(yōu)化配置框架。文獻(xiàn)[6-7]考慮了需求響應(yīng)對(duì)綜合能源系統(tǒng)優(yōu)化調(diào)度的影響。文獻(xiàn)[8]綜合考慮了能效、環(huán)境因素,提出了基于能源價(jià)格響應(yīng)的綜合能源系統(tǒng)優(yōu)化調(diào)度模型。文獻(xiàn)[9]考慮了源荷不確定因素的日前經(jīng)濟(jì)調(diào)度問題,構(gòu)建了基于可信性理論的模糊機(jī)會(huì)約束規(guī)劃優(yōu)化調(diào)度模型。

    現(xiàn)有文獻(xiàn)已經(jīng)對(duì)綜合能源系統(tǒng)的理論建模、規(guī)劃運(yùn)行方面等方面進(jìn)行了研究,但綜合能源系統(tǒng)的系統(tǒng)評(píng)估這一方面的研究仍處于起步階段。文獻(xiàn)[10-11]結(jié)合馬爾可夫和蒙特卡洛對(duì)設(shè)備故障情況建模,分析了設(shè)備故障對(duì)系統(tǒng)可靠性的影響,沿用電力系統(tǒng)的可靠性評(píng)估方法對(duì)系統(tǒng)的可靠性進(jìn)行了評(píng)估。文獻(xiàn)[12]考慮不同子系統(tǒng)的能量損失等級(jí)差異,提出了基于能量質(zhì)量特征的綜合能源系統(tǒng)評(píng)價(jià)指標(biāo)。同樣考慮了能源品質(zhì)差異,文獻(xiàn)[13]將不同品質(zhì)的能源統(tǒng)一折算,構(gòu)建了多能源系統(tǒng)綜合能效的評(píng)估模型。文獻(xiàn)[14]首次將配電網(wǎng)的最大供電能力評(píng)估指標(biāo)引入綜合能源系統(tǒng),為綜合能源系統(tǒng)的供能能力評(píng)估提供了思路。文獻(xiàn)[15]設(shè)計(jì)了綜合能源系統(tǒng)的效益評(píng)價(jià)指標(biāo)體系,根據(jù)層次分析法建立了綜合效益評(píng)價(jià)模型,并基于模糊理論對(duì)系統(tǒng)分級(jí)。

    綜上,現(xiàn)有文獻(xiàn)對(duì)綜合能源系統(tǒng)的研究主要集中在優(yōu)化調(diào)度、綜合評(píng)估方面,但對(duì)綜合能源系統(tǒng)的多種能源綜合可靠供能能力評(píng)估的研究不夠深入;需求響應(yīng)作為調(diào)節(jié)負(fù)荷的重要手段,鮮有文獻(xiàn)對(duì)需求響應(yīng)在高負(fù)荷運(yùn)行下對(duì)綜合能源系統(tǒng)的多能供能能力的提高和作為后備資源保障系統(tǒng)可靠供能的能力進(jìn)行研究。

    鑒于此,現(xiàn)提出一種基于多目標(biāo)差分算法的協(xié)同進(jìn)化算法(co-evolution algorithm based on multi-objective differential evolution algorithm, Coevolution-NSDE)來評(píng)估綜合能源系統(tǒng)的供能能力。利用多目標(biāo)優(yōu)化算法自動(dòng)尋優(yōu)的優(yōu)點(diǎn)解決不同子系統(tǒng)間的耦合問題,在系統(tǒng)設(shè)備約束條件內(nèi)對(duì)可供能峰值進(jìn)行尋優(yōu),來計(jì)算系統(tǒng)的供能區(qū)間。其次,將設(shè)備故障情況納入考慮,在考慮系統(tǒng)運(yùn)行約束的基礎(chǔ)上,對(duì)系統(tǒng)的可靠性需求與系統(tǒng)的供能能力的關(guān)系進(jìn)行研究;并通過對(duì)供暖型熱負(fù)荷參與需求響應(yīng)進(jìn)行分析來研究熱需求響應(yīng)的實(shí)施對(duì)系統(tǒng)可靠供能區(qū)間的影響。最后,通過仿真算例證明,本文所提算法對(duì)于綜合能源系統(tǒng)供能區(qū)間的求解是有效的,且熱需求響應(yīng)作為用戶側(cè)參與優(yōu)化負(fù)荷、合理調(diào)配資源的手段,能夠進(jìn)一步擴(kuò)大系統(tǒng)的可靠供能能力。

    1 綜合能源系統(tǒng)建模

    1.1 綜合能源系統(tǒng)結(jié)構(gòu)組成

    結(jié)合能源集線器思想,本文中綜合能源系統(tǒng)結(jié)構(gòu)如圖1所示,包括供能側(cè),能源轉(zhuǎn)換部分、能源存儲(chǔ)部分以及負(fù)荷側(cè)。

    圖1 綜合能源系統(tǒng)結(jié)構(gòu)圖Fig.1 Structure diagram of the integrated energy system

    本文中所提綜合能源系統(tǒng)模型構(gòu)造如下:①供能側(cè):包括外部電網(wǎng)提供的電能、風(fēng)電、光伏以及外部氣網(wǎng)提供的天然氣;②能量存儲(chǔ)部分:包括儲(chǔ)氫設(shè)備和蓄熱設(shè)備;③能量轉(zhuǎn)換部分:包括燃?xì)怆姛崧?lián)產(chǎn)機(jī)組(cogeneration, combined heat and power,CHP)、燃?xì)鉄岜?、燃料電池、電制氫設(shè)備、氫氣甲烷化設(shè)備。負(fù)荷側(cè)包括熱負(fù)荷、電負(fù)荷、氣負(fù)荷。各設(shè)備具體建模公式及運(yùn)行約束公式見文獻(xiàn)[16-18]。

    1.2 運(yùn)行策略

    綜合能源系統(tǒng)的運(yùn)行需要滿足功率平衡約束,即

    (1)

    式(1)中:LE,buy表示從電網(wǎng)購電;LE,CHP、LE,fc、LE,WT、LE,PV表示CHP、燃料電池、風(fēng)電、光伏產(chǎn)生的電量;LQ,buy表示從電網(wǎng)購氣;LQ,CH4表示甲烷化設(shè)備產(chǎn)生的氣體;LH,CHP、LH,eb、LH,HP、LH,hs表示熱電聯(lián)產(chǎn)、電鍋爐、儲(chǔ)氫罐產(chǎn)生的熱量。

    對(duì)上述綜合能源系統(tǒng),其運(yùn)行過程中,各能源優(yōu)先滿足同類負(fù)荷需求,若不足或有余再通過能源轉(zhuǎn)換設(shè)備轉(zhuǎn)換成其他種類能源滿足負(fù)荷需求或進(jìn)行儲(chǔ)存。

    即電負(fù)荷優(yōu)先由光伏、風(fēng)電滿足,不夠部分由熱電聯(lián)產(chǎn)機(jī)組、燃料電池等滿足,若仍不滿足則從外部電網(wǎng)買電補(bǔ)足;熱負(fù)荷由燃?xì)鉄岜霉?yīng)、不足部分由熱電聯(lián)產(chǎn)機(jī)組、儲(chǔ)熱、電鍋爐供應(yīng);氣負(fù)荷優(yōu)先由外部氣網(wǎng)供給,不足部分由電轉(zhuǎn)氣裝置供給。由于相比于電轉(zhuǎn)氫過程,電轉(zhuǎn)氣過程損耗較小,轉(zhuǎn)化效率較高,故儲(chǔ)能過程中的策略具體為:當(dāng)可再生能源及熱電聯(lián)產(chǎn)機(jī)組所發(fā)電量大于當(dāng)前時(shí)刻電負(fù)荷,多余電量通過電解槽轉(zhuǎn)化為氫氣。當(dāng)儲(chǔ)氫罐儲(chǔ)滿不能繼續(xù)工作或電解槽以最大功率運(yùn)行但系統(tǒng)仍有余電時(shí),通過儲(chǔ)熱消納多余電能。只有當(dāng)供氣不足或儲(chǔ)氫罐及蓄熱裝置不能繼續(xù)工作時(shí),轉(zhuǎn)化的氫氣才通過甲烷化設(shè)備轉(zhuǎn)化為甲烷以滿足氣負(fù)荷;大部分時(shí)間內(nèi),系統(tǒng)將氫氣儲(chǔ)存在儲(chǔ)氫罐中。

    2 基于Coevolution-NSDE的綜合能源可靠供能能力評(píng)估方法

    2.1 系統(tǒng)供能區(qū)間

    綜合能源系統(tǒng)的供能區(qū)間是由多維數(shù)組構(gòu)成的邊界與坐標(biāo)軸圍成的區(qū)間,本文中通過求取系統(tǒng)的多維供能區(qū)間來評(píng)估系統(tǒng)的供能能力。在本文中,這個(gè)多維數(shù)組為三維數(shù)組x=[x1,x2,x3],分別代表電、熱、氣三種供能需求。即供電區(qū)間為[0,x1];供熱區(qū)間為[0,x2];供氣區(qū)間為[0,x3]。

    綜合能源系統(tǒng)負(fù)荷種類多樣,系統(tǒng)電、熱、氣供能情況不僅具有耦合特性,能夠互濟(jì)互聯(lián),而且各能源間的轉(zhuǎn)換關(guān)系并非是簡(jiǎn)單線性關(guān)系。從內(nèi)部設(shè)備層面來看,設(shè)備自身具有復(fù)雜的約束條件,并且儲(chǔ)氫、儲(chǔ)熱設(shè)備具有時(shí)序特性,難以通過數(shù)學(xué)解析的方法得到最大供能區(qū)間邊界;從系統(tǒng)層面來看,某種能源供應(yīng)能力大小受其余種類能源供給能力的影響,系統(tǒng)供能能力不存在唯一最優(yōu)解。

    為了解決上述問題,在滿足系統(tǒng)配置的設(shè)備的約束條件的基礎(chǔ)上,將系統(tǒng)能夠供應(yīng)的最大負(fù)荷能力作為系統(tǒng)的供能能力進(jìn)行系統(tǒng)的最大供能區(qū)間求解。以綜合能源系統(tǒng)的供電、熱、氣能力作為優(yōu)化變量,采用多目標(biāo)優(yōu)化算法對(duì)進(jìn)行求解。

    將優(yōu)化變量xi乘以典型日的時(shí)序系數(shù)得到全年的負(fù)荷曲線,進(jìn)行系統(tǒng)的運(yùn)行仿真。即優(yōu)化變量xi代表系統(tǒng)可供應(yīng)范圍內(nèi)的負(fù)荷峰值個(gè)體,通過不斷的優(yōu)化迭代,得到系統(tǒng)可以滿足的最大負(fù)荷峰值即所求的目標(biāo),系統(tǒng)供能能力。具體年負(fù)荷曲線生成公式為

    Li(t)=xiLi,i=1, 2, 3,…

    (2)

    式(2)中:Li(t)表示全年各時(shí)刻負(fù)荷,xi表示優(yōu)化變量;Li表示典型日各時(shí)刻負(fù)荷與負(fù)荷峰值比值;i分別表示三種負(fù)荷。由于熱負(fù)荷具有季節(jié)特性,為了貼近實(shí)際,不同季節(jié)應(yīng)使用不同負(fù)荷比值曲線表示。

    優(yōu)化目標(biāo)是使系統(tǒng)可供應(yīng)電、熱、氣負(fù)荷能力最大。參變量包括迭代次數(shù)、個(gè)體數(shù)目、綜合能源系統(tǒng)配置等。優(yōu)化結(jié)果所得個(gè)體均滿足系統(tǒng)設(shè)備約束以及可靠性約束,其中每個(gè)個(gè)體都代表了系統(tǒng)一種最大可供應(yīng)負(fù)荷能力方案。但由于各能源彼此耦合相互約束關(guān)系,最終得到的pareto最優(yōu)解集為一個(gè)多維數(shù)組的集合,代表了系統(tǒng)能夠供應(yīng)的最大負(fù)荷邊界。

    2.2 求解算法Coevolution-NSDE

    為了求得系統(tǒng)在功率平衡約束及設(shè)備運(yùn)行約束條件下的最大供能區(qū)間,將系統(tǒng)能夠供應(yīng)的最大電、熱、氣負(fù)荷作為尋優(yōu)目標(biāo),這些尋優(yōu)目標(biāo)是互相影響、互相競(jìng)爭(zhēng)、耦合在一起的,它們的競(jìng)爭(zhēng)性和復(fù)雜性使得其優(yōu)化求解變得困難。采用多目標(biāo)優(yōu)化算法,利用群智能算法的自動(dòng)尋優(yōu)特性解決各個(gè)能源子系統(tǒng)的耦合問題進(jìn)行求解,流程如圖2所示。

    圖2 Coevolution-NSDE算法求解流程圖Fig.2 Coevolution-NSDE algorithm solving flow chart

    為了保持種群的多樣性,得到更為完整的pareto前沿,本文中采用基于多目標(biāo)差分進(jìn)化算法的協(xié)同進(jìn)化算法進(jìn)行求解。

    非支配排序差分進(jìn)化算法(non-dominated sorting differential evolution algorithm, NSDE)將差分進(jìn)化算法(differential evolution algorithm, DE)與快速非支配排序操作相融合,具有收斂速度快、pareto最優(yōu)解分布特性好等特點(diǎn),在電力系統(tǒng)相關(guān)研究領(lǐng)域得到了大量應(yīng)用。然而,NSDE 算法仍存在易于早熟、全局收斂能力差等缺點(diǎn)。故本文中采用協(xié)同進(jìn)化提高其全局搜索能力以得到更全面的pareto解集。

    選取非支配排序遺傳算法(non-dominated sorting genetic algorithm, NSGA-II)生成輔助種群進(jìn)行協(xié)同進(jìn)化,為NSDE的子種群增加多樣性,同時(shí)也提供了NSDE子種群跳出局部最優(yōu)的可能性,避免過早收斂。

    具體優(yōu)化求解過程如下:①NSDE與NSGA-Ⅱ算法設(shè)置相同的種群數(shù)量N、變量范圍等初始化參數(shù)。各自生成初始化種群;②二者各自進(jìn)化生成子代種群;③將二者子代種群混合產(chǎn)生2N個(gè)子代,根據(jù)非支配排序?qū)旌献哟N群進(jìn)行排序;④選取最優(yōu)N個(gè)子代作為NSDE的子代,NSGA-Ⅱ的種群子代仍為NSGA-Ⅱ的原始子代;⑤重復(fù)步驟②~步驟④直至到達(dá)迭代終止條件。

    2.3 系統(tǒng)可靠供能區(qū)間

    實(shí)際情況中,評(píng)估綜合能源系統(tǒng)供能能力必須要考慮系統(tǒng)設(shè)備的狀態(tài)。當(dāng)系統(tǒng)設(shè)備出現(xiàn)故障時(shí),可能使系統(tǒng)無法如理想狀態(tài)供能,導(dǎo)致供能能力達(dá)不到預(yù)期。但負(fù)荷峰值及系統(tǒng)設(shè)備故障時(shí)間在整個(gè)仿真時(shí)間中只占很小的一部分,若要求系統(tǒng)時(shí)時(shí)滿足負(fù)荷需求要求過于苛刻。故在一定的可靠指標(biāo)約束下求取的系統(tǒng)供能區(qū)間更能挖掘系統(tǒng)的供能能力,也更具有實(shí)際意義。

    計(jì)算系統(tǒng)可靠性一般用解析法和模擬法。由于本文所討論的綜合能源系統(tǒng)包含儲(chǔ)氫、蓄熱設(shè)備,綜合考慮系統(tǒng)的時(shí)序特性和計(jì)算的準(zhǔn)確性,采用序貫蒙特卡洛方法進(jìn)行可靠性指標(biāo)計(jì)算。負(fù)荷峰值越大,當(dāng)系統(tǒng)設(shè)備出現(xiàn)故障時(shí),其造成的系統(tǒng)期望缺供量也會(huì)隨之增大。故以系統(tǒng)期望缺供量為約束條件會(huì)在一定程度上對(duì)系統(tǒng)最大供應(yīng)負(fù)荷能力起到限制作用。以馬爾可夫兩狀態(tài)模型描述設(shè)備狀態(tài)的不確定性,將可靠性指標(biāo)作為約束條件之一進(jìn)行供能區(qū)間求解。

    假設(shè)設(shè)備正常工作時(shí)間與故障持續(xù)時(shí)間均服從指數(shù)分布,其蒙特卡洛抽樣過程公式為

    (3)

    (4)

    式中:下標(biāo)x代表設(shè)備類型;Tg,x為設(shè)備正常持續(xù)時(shí)間;Tr,x為設(shè)備修復(fù)時(shí)間;λx為設(shè)備故障率;δx為設(shè)備修復(fù)率;σx為0~1的隨機(jī)數(shù)。通過式(3)、式(4)可以得出仿真時(shí)間內(nèi)系統(tǒng)各設(shè)備的時(shí)序狀態(tài)。

    為求解系統(tǒng)可靠供能區(qū)間,以系統(tǒng)期望缺供量(expected demand not supplied,EENS)這一可靠性指標(biāo)作為約束條件,系統(tǒng)期望缺供量是指統(tǒng)計(jì)時(shí)間內(nèi)系統(tǒng)由于能量供應(yīng)不足而削減的部分負(fù)荷量。在本文系統(tǒng)中具體為系統(tǒng)期望缺供電量、系統(tǒng)缺供熱量和系統(tǒng)缺供氣量。計(jì)算公式為

    (5)

    式(5)中:T為仿真時(shí)長(zhǎng);E為系統(tǒng)的缺負(fù)荷量;i為能源種類。

    2.4 熱需求響應(yīng)參與下的系統(tǒng)可靠供能區(qū)間

    在綜合能源系統(tǒng)中,負(fù)荷具有多種形式。由于熱傳導(dǎo)的延遲特性,部分供暖類型熱負(fù)荷也具有可調(diào)控潛力。與電負(fù)荷不同,當(dāng)系統(tǒng)供熱設(shè)備發(fā)生故障,或供熱能力不足時(shí),溫度作為熱負(fù)荷的呈現(xiàn)形式并不會(huì)突變,而是在一定時(shí)間內(nèi)呈非線性下降趨勢(shì)。對(duì)負(fù)荷側(cè)而言,一定范圍內(nèi)的溫度波動(dòng)雖然對(duì)用戶體驗(yàn)有所影響,但是可接受的,即用戶側(cè)供能故障發(fā)生時(shí)間滯后于設(shè)備故障發(fā)生時(shí)間。

    如圖3所示,若供熱設(shè)備在t1時(shí)刻發(fā)生故障,對(duì)于用戶而言,在t2時(shí)刻介質(zhì)溫度才降低到允許范圍之下,超出允許范圍。若不考慮介質(zhì)熱傳導(dǎo)延遲特性,實(shí)時(shí)調(diào)整供熱量以維持介質(zhì)溫度處于某一恒定值,則會(huì)造成能源的浪費(fèi),降低系統(tǒng)的整體供能能力。故在評(píng)估系統(tǒng)可靠性時(shí),將這種允許范圍內(nèi)的溫度波動(dòng)列入考慮是更符合實(shí)際情況的。對(duì)于綜合能源系統(tǒng),由于各能源之間相互耦合,考慮熱負(fù)荷響應(yīng)不僅影響系統(tǒng)供熱能力,當(dāng)系統(tǒng)設(shè)備有故障發(fā)生時(shí),系統(tǒng)產(chǎn)生其他類型負(fù)荷缺供情況,允許系統(tǒng)供熱量在一定范圍內(nèi)偏離熱負(fù)荷,優(yōu)先滿足其他負(fù)荷的供給,即作為柔性熱負(fù)荷使用。

    圖3 熱傳導(dǎo)延遲特性Fig.3 Thermal conduction delay characteristics

    以熱水系統(tǒng)作為典型的熱負(fù)荷模型,即

    (6)

    (7)

    式中:T(t-1)、T(t)分別為上一時(shí)刻、這一時(shí)刻的熱水溫度;Qmax,H為熱容器額定熱容量;R為熱容器熱阻;Tout(t)、Tin為環(huán)境溫度和進(jìn)水溫度;Te為額定溫度;Ph為輸入熱功率。

    在考慮熱負(fù)荷響應(yīng)時(shí),系統(tǒng)的運(yùn)行策略具體為:當(dāng)系統(tǒng)產(chǎn)生缺供時(shí),判斷負(fù)荷側(cè)供暖型熱負(fù)荷供應(yīng)溫度是否在正常范圍內(nèi),若處于可接受范圍,計(jì)算若在可接受范圍內(nèi)降低熱負(fù)荷供應(yīng),能否滿足缺供量。即以溫度作為評(píng)判缺供與否的標(biāo)志。

    考慮熱需求響應(yīng)參與的綜合能源系統(tǒng)可靠供能區(qū)間求取流程如圖4所示。

    圖4 考慮熱需求響應(yīng)的系統(tǒng)可靠供能區(qū)間求解流程圖Fig.4 Flow chart for solving system’s reliable energy supply interval considering heat demand response

    如圖4所示,求取可靠供能區(qū)間步驟如下。

    (1)更新三維個(gè)體,分別代表系統(tǒng)可供應(yīng)的電、熱、氣負(fù)荷年峰值。

    (2)通過負(fù)荷比例曲線將個(gè)體代表的可供負(fù)荷峰值轉(zhuǎn)化為年負(fù)荷曲線。

    (3)對(duì)系統(tǒng)設(shè)備進(jìn)行序貫蒙特卡洛模擬抽樣,得到設(shè)備的時(shí)序狀態(tài)。根據(jù)事先設(shè)定好的系統(tǒng)配置、運(yùn)行策略,在仿真年限內(nèi)進(jìn)行運(yùn)行仿真。

    (4)運(yùn)行仿真結(jié)束,計(jì)算可靠性指標(biāo)。

    (5)判斷生成的可靠性指標(biāo)是否滿足設(shè)定的約束條件。

    (6)若不滿足,重新生成個(gè)體,返回步驟(1)。

    (7)比較進(jìn)化結(jié)果,擇優(yōu)選取。

    3 算例分析

    3.1 算法性能對(duì)比

    為測(cè)試所提Coevolution-NSDE算法性能,以不考慮需求響應(yīng)和可靠性約束的基礎(chǔ)最大供能區(qū)間進(jìn)行仿真對(duì)比,采用解集覆蓋指標(biāo)(C-metric)和間距評(píng)估指標(biāo)(Sa)[19]對(duì)基礎(chǔ)NSDE和Coevolution-NSDE二者按圖4所述流程進(jìn)行仿真。其中,設(shè)備具體參數(shù)如表1所示。

    表1 綜合能源系統(tǒng)設(shè)備參數(shù)

    典型日負(fù)荷比例曲線如圖5所示,全年風(fēng)光出力如圖6所示。

    圖5 典型日負(fù)荷比例曲線Fig.5 Typical daily load ratio curve

    圖6 全年風(fēng)光出力圖Fig.6 Annual wind power and photovoltaic output

    (1)解集覆蓋指標(biāo)(C-metric)。用來比較兩個(gè)解集哪個(gè)更接近真實(shí)pareto解集。解集覆蓋指標(biāo)C(A,B)代表B中解被A中解支配的個(gè)數(shù)占B非劣解集中的比例。若C(A,B)>C(B,A)則說明解集A表現(xiàn)優(yōu)于解集B,反之亦然。計(jì)算公式為

    (8)

    式(8)中:A、B表示兩個(gè)非劣解集;分子表示B中解被A中解支配的個(gè)數(shù);分母表示B中解的個(gè)數(shù)。

    對(duì)兩個(gè)算法采取同樣的種群數(shù)量200和迭代次數(shù)50次,得到結(jié)果如表2所示??梢姳疚姆椒ㄋ媒饧瘍?yōu)于NSDE解集。

    表2 算法的解集覆蓋指標(biāo)對(duì)比

    (2)間距評(píng)估指標(biāo)(Sa)。間距評(píng)估指標(biāo)用來評(píng)估指標(biāo)中的解的均勻性,指標(biāo)越小代表解集均勻性越好。計(jì)算公式為

    (9)

    式(9)中:I為非劣解集中解的個(gè)數(shù);d為解集的相鄰歐氏距離;dave為相鄰歐氏距離的平均值。

    對(duì)兩個(gè)算法采取同樣的種群數(shù)量200和迭代次數(shù)50次,得到結(jié)果如表3所示。

    表3 算法的間距評(píng)估指標(biāo)指標(biāo)對(duì)比

    本文所得的Sa指標(biāo)小于NSDE所得指標(biāo),可見本文方法所得解集有著更好的均勻性。

    3.2 系統(tǒng)可靠供能區(qū)間

    以不同的系統(tǒng)年期望缺供量約束條件代表系統(tǒng)的不同可靠性需求程度。系統(tǒng)期望缺供越小,代表系統(tǒng)對(duì)供能可靠性的需求越高。以年期望缺供電、熱、氣量分別為1 500 kW、1 500 kW和1 500 m3;年期望缺供電、熱、氣量為1 000 kW、1 000 kW和1 000 m3和年期望缺供電、熱、氣量為500 kW、500 kW和500 m3為約束條件對(duì)系統(tǒng)進(jìn)行求解。每種約束條件下各100個(gè)個(gè)體。各設(shè)備故障率及修復(fù)率如表4所示,仿真時(shí)間為200年。

    表4 綜合能源系統(tǒng)設(shè)備可靠性參數(shù)

    采用第2節(jié)所述的Coevolution-NSDE算法對(duì)約束條件的最大供能區(qū)間進(jìn)行求解。得到結(jié)果如圖7所示。

    圖7 不同可靠性需求約束的系統(tǒng)可靠供能區(qū)間Fig.7 Reliable energy supply interval of the system with different reliability demand constraints

    其中,在年期望缺供電量為1 500 kW,年期望缺供熱量1 500 kW,年期望缺供氣量1 500 m3以內(nèi)的約束條件下,得到的最大可供電能峰值為3 010.5 kW,最大可供熱能峰值為2 101.6 kW,最大可供氣能負(fù)荷為1 100.7 m3;在年期望缺供電量1 000 kW,年期望缺供熱量1 000 kW,年期望缺供1 000 m3以內(nèi)的約束條件下,得到的最大可供電能峰值為2 992.6 kW,最大可供熱能峰值為2 079.1 kW,最大可供氣能負(fù)荷為1 006.5 m3;在年期望缺供電量500 kW,年期望缺供熱量500 kW,年期望缺供500 m3以內(nèi)的約束條件下,得到的最大可供電能峰值為2 980.5 kW,最大可供熱能峰值為1 993.9 kW,最大可供氣能負(fù)荷為981.0 m3??梢姡S著系統(tǒng)的可靠性需求提高,系統(tǒng)的最大供能能力有所下降。在不改變系統(tǒng)配置的條件下,系統(tǒng)的可靠性與供能能力呈負(fù)相關(guān)關(guān)系。

    為了更清楚地得到各種類能源最大供能能力的關(guān)系及對(duì)比不同的可靠性需求約束對(duì)系統(tǒng)供能能力的影響,固定年供電能峰值1 500 kW,得到熱-氣供能能力在不同可靠性需求約束下的差異,如圖8(a)所示。同樣地,固定年供熱能峰值為1 500 kW,得到電-氣供能能力在不同可靠性需求約束下的差異,如圖8(b)所示;固定年供氣能峰值為700 m3時(shí),得到電-熱供能能力在不同可靠性需求約束下的差異,如圖8(c)所示。用第2節(jié)的多目標(biāo)求解方法進(jìn)行仿真求解,結(jié)果如圖8所示。

    圖8 不同可靠性需求約束下的系統(tǒng)最大供能區(qū)間邊界Fig.8 The boundary of the maximum energy supply interval of the system under the constraints of different reliability requirements

    由圖8可以看出,隨著可靠性需求的提高,系統(tǒng)的供電、熱、氣能力都有所下降,且從圖8(a)中可以看出系統(tǒng)的最大供熱、氣關(guān)系呈負(fù)相關(guān)關(guān)系。這是由于系統(tǒng)供熱沒有直接熱源,熱能的供給需要通過各能源轉(zhuǎn)化設(shè)備將電、氣能轉(zhuǎn)化為熱能。熱能需求量越高,則需要越多的能源轉(zhuǎn)化設(shè)備運(yùn)作,其對(duì)系統(tǒng)能源轉(zhuǎn)換設(shè)備的可靠性要求越高。圖8(b)所示的系統(tǒng)最大供電、氣能力其相關(guān)性則不明顯,僅在最大供電能力小于500 kW時(shí),存在負(fù)相關(guān)特性,這是由于本文系統(tǒng)電氣能源轉(zhuǎn)換僅配備了熱電聯(lián)產(chǎn)和電轉(zhuǎn)氣設(shè)備,而熱電聯(lián)產(chǎn)的運(yùn)行策略為以熱定電,故在系統(tǒng)最大供熱能力固定的情況下,電氣能源轉(zhuǎn)換主要為電轉(zhuǎn)氣設(shè)備,僅在系統(tǒng)供電能源較小時(shí)電轉(zhuǎn)氣設(shè)備才會(huì)作用。圖8(c)所示的系統(tǒng)最大供電、熱能力在電能小于1 500 kW、熱能大于1 450 kW 時(shí),呈現(xiàn)較強(qiáng)的負(fù)相關(guān)關(guān)系,這是由于在本文配置的系統(tǒng)中,熱電能源的轉(zhuǎn)換主要是依靠電鍋爐完成的,根據(jù)優(yōu)先供給本類型能源的原則,僅在熱能需求高且電能需求低時(shí)電鍋爐工作進(jìn)行能源轉(zhuǎn)換,符合實(shí)際情況。

    3.3 需求響應(yīng)參與下的可靠供能區(qū)間

    本文中假設(shè)供暖型熱負(fù)荷占總熱負(fù)荷比例取70%。期望溫度為75 ℃,允許溫度范圍55~75 ℃??紤]熱負(fù)荷響應(yīng)時(shí),以溫度為缺供標(biāo)志進(jìn)行判斷;未考慮熱負(fù)荷響應(yīng)時(shí),以熱負(fù)荷與系統(tǒng)實(shí)際供能差值為缺供標(biāo)志進(jìn)行判斷。當(dāng)設(shè)備故障、或可再生能源出力低時(shí),供暖部分熱負(fù)荷參與需求響應(yīng),通過在溫度變化允許范圍內(nèi)短時(shí)間的降低熱負(fù)荷供給,以溫度作為缺供與否的評(píng)判標(biāo)準(zhǔn)更能準(zhǔn)確反映熱負(fù)荷的用本質(zhì)。

    對(duì)不同供暖型熱負(fù)荷占比的系統(tǒng)進(jìn)行最大供能區(qū)間的求解,當(dāng)系統(tǒng)不能滿足供能需求時(shí),可通過降低供暖型負(fù)荷的供熱,以滿足其余負(fù)荷需求。在年期望缺供電量為1 000 kW,年期望缺供熱量1 000 kW,年期望缺供氣量1 000 m3以內(nèi)的約束條件下,供暖型負(fù)荷占比w分別為0.7、0.5、0.3,得到仿真結(jié)果如圖9所示。

    圖9 不同供暖型熱負(fù)荷占比的系統(tǒng)可靠供能區(qū)間Fig.9 Reliable energy supply interval of the system with different heating load proportions

    其中,在當(dāng)供暖型熱負(fù)荷占比0.7時(shí),最大可供電峰值為3 040.2 kW,最大可供熱能峰值為2 123.6 kW,最大可供氣能負(fù)荷為1 153.7 m3;在當(dāng)供暖型熱負(fù)荷占比0.5時(shí),得到的最大可供電能峰值為3 021.8 kW,最大可供熱能峰值為2 101.4 kW,最大可供氣能負(fù)荷為1 127.9 m3;在當(dāng)供暖型熱負(fù)荷占比0.3時(shí),得到的最大可供電能峰值為2 999.7 kW,最大可供熱能峰值為2 086.9 kW,最大可供氣能負(fù)荷為1 099.4 m3??梢?,隨著系統(tǒng)的供暖型熱負(fù)荷比例增大,系統(tǒng)最大供能能力也有所增加。通過合理調(diào)用需求側(cè)資源,利用需求側(cè)響應(yīng)的參與使得系統(tǒng)資源分配效率提高,使得系統(tǒng)的最大供能邊界增大。

    同樣地,為了更清楚地對(duì)比不同的供暖型負(fù)荷占比對(duì)系統(tǒng)供能能力的影響,固定年供電能峰值1 500 kW,得到熱-氣供能能力在不同供暖型負(fù)荷占比下的差異,如圖10(a)所示;固定年供熱能峰值為1 500 kW,得到電-氣供能能力在不同供暖型負(fù)荷占比下的差異,如圖10(b)所示;固定年供氣能峰值為700 m3時(shí),得到電-熱供能能力在不同供暖型負(fù)荷占比下的差異,如圖10(c)所示。用第2節(jié)的多目標(biāo)求解方法進(jìn)行仿真求解,結(jié)果如圖10所示。

    圖10 不同供暖型熱負(fù)荷占比系統(tǒng)最大可靠供能區(qū)間邊界Fig.10 The boundary of the maximum reliable energy supply interval of different heating systems with different heating load proportions

    由圖10可以看出,以溫度作為缺供判斷標(biāo)志,將供暖型熱負(fù)荷作為可替代備用能源,從需求側(cè)進(jìn)一步優(yōu)化系統(tǒng)的運(yùn)行能力,可以一定程度上提高系統(tǒng)的供熱能力。并且由于綜合能源系統(tǒng)各類型能源可相互轉(zhuǎn)換,系統(tǒng)供電、供氣能力也會(huì)隨之有所提高。

    由圖10可以看出,隨著供暖型負(fù)荷占比的增加,系統(tǒng)的最大供電、熱、氣能力均有所增加,且由圖10(a)、圖10(c)可以看出隨著系統(tǒng)供熱的增加,供暖型負(fù)荷占比對(duì)系統(tǒng)供電、氣能力的影響愈發(fā)明顯。這是由于系統(tǒng)供熱增加,供暖型負(fù)荷也隨之增加,即在系統(tǒng)缺供時(shí),系統(tǒng)將有更多的供暖能源在允許范圍內(nèi)支援其他類型缺供能源。

    4 結(jié)論

    本文提出了基于多目標(biāo)差分進(jìn)化算法的協(xié)同進(jìn)化算法對(duì)綜合能源系統(tǒng)的可靠供能區(qū)間進(jìn)行求解。分別對(duì)不同占比的供暖型熱負(fù)荷參與需求響應(yīng)的情況進(jìn)行仿真,分析熱需求響應(yīng)參與對(duì)系統(tǒng)的可靠供能能力的影響,可得下述結(jié)論。

    (1)采用基于非支配排序多目標(biāo)差分進(jìn)化算法的協(xié)同進(jìn)化算法對(duì)綜合能源系統(tǒng)的可靠供能區(qū)間進(jìn)行求解,通過多目標(biāo)優(yōu)化算法自動(dòng)尋優(yōu)的優(yōu)點(diǎn)解決了綜合能源系統(tǒng)復(fù)雜的耦合特性為求解帶來的困難。該算法在綜合能源系統(tǒng)的多能供能區(qū)間求解中是有效的。在綜合能源系統(tǒng)中,各類型能源的供能能力呈現(xiàn)負(fù)相關(guān)特性。

    (2)通過結(jié)合設(shè)備狀態(tài)的不確定性進(jìn)一步分析了綜合能源系統(tǒng)的供能能力,系統(tǒng)的供能能力隨著供能的可靠性需求的增高而下降。

    (3)供暖型熱負(fù)荷作為一種需求響應(yīng)的手段參與到能源的調(diào)控中,能夠進(jìn)一步擴(kuò)大系統(tǒng)的可靠供能區(qū)間,優(yōu)化綜合能源系統(tǒng)的運(yùn)行,提高系統(tǒng)合理配置資源的能力。

    在后續(xù)研究中,計(jì)劃進(jìn)一步考慮需求響應(yīng)的不確定性及風(fēng)光的不確定性對(duì)系統(tǒng)供能的影響。

    猜你喜歡
    區(qū)間峰值可靠性
    解兩類含參數(shù)的復(fù)合不等式有解與恒成立問題
    “四單”聯(lián)動(dòng)打造適齡兒童隊(duì)前教育峰值體驗(yàn)
    你學(xué)會(huì)“區(qū)間測(cè)速”了嗎
    可靠性管理體系創(chuàng)建與實(shí)踐
    電子制作(2017年2期)2017-05-17 03:55:06
    寬占空比峰值電流型準(zhǔn)PWM/PFM混合控制
    基于峰值反饋的電流型PFM控制方法
    區(qū)間對(duì)象族的可鎮(zhèn)定性分析
    基于可靠性跟蹤的薄弱環(huán)節(jié)辨識(shí)方法在省級(jí)電網(wǎng)可靠性改善中的應(yīng)用研究
    可靠性比一次采購成本更重要
    風(fēng)能(2015年9期)2015-02-27 10:15:24
    国产在线免费精品| 成年女人毛片免费观看观看9 | 我的亚洲天堂| 建设人人有责人人尽责人人享有的| 亚洲欧美色中文字幕在线| 欧美黄色片欧美黄色片| 国产成人精品久久久久久| 操出白浆在线播放| 啦啦啦中文免费视频观看日本| 国产精品99久久99久久久不卡 | 男女高潮啪啪啪动态图| 777久久人妻少妇嫩草av网站| 97在线人人人人妻| 老司机影院毛片| 国产一区二区三区av在线| 日韩中文字幕视频在线看片| 夜夜骑夜夜射夜夜干| 纯流量卡能插随身wifi吗| 亚洲欧美成人精品一区二区| 久久天躁狠狠躁夜夜2o2o | 色婷婷久久久亚洲欧美| tube8黄色片| 久久久久久久久久久免费av| 国产一区二区在线观看av| 色视频在线一区二区三区| 九色亚洲精品在线播放| 两个人免费观看高清视频| 免费观看人在逋| 又大又黄又爽视频免费| 精品少妇久久久久久888优播| 一二三四在线观看免费中文在| 波多野结衣av一区二区av| 久久久国产精品麻豆| 亚洲美女黄色视频免费看| 欧美另类一区| 一边摸一边做爽爽视频免费| 人妻一区二区av| 国产精品 国内视频| 精品午夜福利在线看| 一级毛片我不卡| 女人被躁到高潮嗷嗷叫费观| 国产探花极品一区二区| 免费在线观看黄色视频的| 人人妻人人添人人爽欧美一区卜| 美女视频免费永久观看网站| 亚洲国产日韩一区二区| 国产国语露脸激情在线看| av天堂久久9| 一级a爱视频在线免费观看| 国产一区有黄有色的免费视频| 久久久国产精品麻豆| av福利片在线| 热re99久久国产66热| 亚洲国产成人一精品久久久| 国产免费现黄频在线看| 成人手机av| 国产人伦9x9x在线观看| 青春草国产在线视频| 成人午夜精彩视频在线观看| 亚洲精品视频女| 久久久久精品国产欧美久久久 | 国产精品三级大全| 亚洲精品成人av观看孕妇| 精品一品国产午夜福利视频| 麻豆av在线久日| 国产精品免费大片| 99久国产av精品国产电影| 国产野战对白在线观看| 丝袜美腿诱惑在线| 中文精品一卡2卡3卡4更新| 国产在视频线精品| 精品人妻一区二区三区麻豆| 欧美激情高清一区二区三区 | 成年动漫av网址| 成年美女黄网站色视频大全免费| 欧美日韩亚洲国产一区二区在线观看 | 日韩不卡一区二区三区视频在线| 丝袜美足系列| 亚洲精品国产av蜜桃| 看免费成人av毛片| 国产精品秋霞免费鲁丝片| 91精品国产国语对白视频| 亚洲成人一二三区av| 不卡视频在线观看欧美| 一区二区av电影网| 欧美人与性动交α欧美软件| 日本欧美视频一区| 高清黄色对白视频在线免费看| 亚洲国产欧美在线一区| 少妇 在线观看| 少妇被粗大的猛进出69影院| 国产精品无大码| 久久综合国产亚洲精品| 建设人人有责人人尽责人人享有的| 亚洲精华国产精华液的使用体验| 国产熟女午夜一区二区三区| 欧美日韩国产mv在线观看视频| 乱人伦中国视频| 免费高清在线观看日韩| 午夜激情久久久久久久| 国产乱来视频区| a级片在线免费高清观看视频| 麻豆av在线久日| 久久ye,这里只有精品| 在线免费观看不下载黄p国产| 免费看不卡的av| 又大又爽又粗| 看免费av毛片| 极品少妇高潮喷水抽搐| 国产女主播在线喷水免费视频网站| 亚洲综合色网址| 下体分泌物呈黄色| 男女床上黄色一级片免费看| 精品国产一区二区久久| av电影中文网址| 搡老岳熟女国产| 国产成人精品无人区| 2018国产大陆天天弄谢| 热re99久久精品国产66热6| 不卡视频在线观看欧美| 日韩不卡一区二区三区视频在线| 老司机影院成人| 成人午夜精彩视频在线观看| 看非洲黑人一级黄片| 久久 成人 亚洲| 嫩草影视91久久| 人体艺术视频欧美日本| √禁漫天堂资源中文www| 一边摸一边做爽爽视频免费| 成人亚洲精品一区在线观看| 大片电影免费在线观看免费| 久久久久久免费高清国产稀缺| 深夜精品福利| 成人国产av品久久久| 国产精品.久久久| 成人国产麻豆网| 国产免费一区二区三区四区乱码| 午夜影院在线不卡| 51午夜福利影视在线观看| 精品午夜福利在线看| 曰老女人黄片| 久久天堂一区二区三区四区| 婷婷色麻豆天堂久久| 免费av中文字幕在线| 亚洲精品成人av观看孕妇| 亚洲,欧美精品.| 国产精品成人在线| av国产精品久久久久影院| 蜜桃国产av成人99| 在线天堂中文资源库| 少妇被粗大猛烈的视频| 卡戴珊不雅视频在线播放| 乱人伦中国视频| 香蕉丝袜av| 成年女人毛片免费观看观看9 | 黄色 视频免费看| 岛国毛片在线播放| 午夜激情久久久久久久| 国产成人欧美| 欧美日韩综合久久久久久| 久久久久国产一级毛片高清牌| 亚洲欧美一区二区三区国产| 精品酒店卫生间| 免费不卡黄色视频| 亚洲精品第二区| 亚洲精品美女久久av网站| 日韩一卡2卡3卡4卡2021年| 无限看片的www在线观看| 精品一品国产午夜福利视频| 欧美精品一区二区免费开放| 婷婷色综合www| 久久久精品免费免费高清| 一区二区三区精品91| 国产1区2区3区精品| 国产野战对白在线观看| 中文欧美无线码| 久久国产亚洲av麻豆专区| 午夜福利视频在线观看免费| 制服诱惑二区| 91精品伊人久久大香线蕉| av在线老鸭窝| 国产 一区精品| 好男人视频免费观看在线| 男人添女人高潮全过程视频| 操出白浆在线播放| 日韩 亚洲 欧美在线| 深夜精品福利| 又粗又硬又长又爽又黄的视频| 韩国高清视频一区二区三区| 你懂的网址亚洲精品在线观看| 国产男女内射视频| 久久久精品区二区三区| 免费观看人在逋| 看免费成人av毛片| 国产精品免费大片| 狠狠精品人妻久久久久久综合| 免费观看人在逋| 久久ye,这里只有精品| 午夜精品国产一区二区电影| 一级片免费观看大全| xxxhd国产人妻xxx| 日韩精品免费视频一区二区三区| 1024视频免费在线观看| 9191精品国产免费久久| 久久av网站| 亚洲欧美成人综合另类久久久| 老司机亚洲免费影院| 亚洲欧美成人精品一区二区| 可以免费在线观看a视频的电影网站 | 久久人妻熟女aⅴ| 久久久精品94久久精品| 十八禁高潮呻吟视频| 满18在线观看网站| 国产免费视频播放在线视频| 在线观看免费视频网站a站| 男人操女人黄网站| 下体分泌物呈黄色| 99国产综合亚洲精品| 中文精品一卡2卡3卡4更新| 国产成人欧美| 日韩大码丰满熟妇| 亚洲自偷自拍图片 自拍| 亚洲第一青青草原| 高清av免费在线| 女的被弄到高潮叫床怎么办| 日韩大码丰满熟妇| 五月天丁香电影| 黄色一级大片看看| 国产精品一区二区在线不卡| 免费少妇av软件| 别揉我奶头~嗯~啊~动态视频 | 亚洲婷婷狠狠爱综合网| 十八禁网站网址无遮挡| 欧美国产精品一级二级三级| 日本91视频免费播放| 91国产中文字幕| 老司机深夜福利视频在线观看 | 久热这里只有精品99| 精品国产一区二区久久| 亚洲av综合色区一区| 黄频高清免费视频| 国产成人a∨麻豆精品| 成人午夜精彩视频在线观看| 国产男人的电影天堂91| 国产亚洲精品第一综合不卡| 国产国语露脸激情在线看| 国产黄色视频一区二区在线观看| 欧美日韩一区二区视频在线观看视频在线| 日韩成人av中文字幕在线观看| 在线观看一区二区三区激情| 亚洲av成人精品一二三区| 精品少妇一区二区三区视频日本电影 | 精品国产一区二区久久| 亚洲熟女毛片儿| 在线天堂最新版资源| 黑人欧美特级aaaaaa片| 考比视频在线观看| 又粗又硬又长又爽又黄的视频| 狠狠精品人妻久久久久久综合| 亚洲激情五月婷婷啪啪| 超碰成人久久| 国产熟女欧美一区二区| 久久97久久精品| 另类亚洲欧美激情| xxx大片免费视频| 久久久久久久国产电影| 成人黄色视频免费在线看| 一边摸一边抽搐一进一出视频| 国产在线免费精品| 久久久久久久精品精品| 赤兔流量卡办理| 国产成人欧美在线观看 | 亚洲专区中文字幕在线 | 美女视频免费永久观看网站| 国产熟女欧美一区二区| 欧美另类一区| 国产亚洲av高清不卡| 人体艺术视频欧美日本| 观看美女的网站| 亚洲美女搞黄在线观看| 中文字幕人妻熟女乱码| av电影中文网址| avwww免费| 亚洲av中文av极速乱| 亚洲国产欧美网| 欧美精品高潮呻吟av久久| 日韩不卡一区二区三区视频在线| 在线 av 中文字幕| 色网站视频免费| 欧美日韩亚洲国产一区二区在线观看 | 成人亚洲欧美一区二区av| 9色porny在线观看| 日本91视频免费播放| 99热国产这里只有精品6| 日韩免费高清中文字幕av| 午夜福利乱码中文字幕| 色综合欧美亚洲国产小说| 日韩一区二区三区影片| 永久免费av网站大全| 制服人妻中文乱码| 精品亚洲乱码少妇综合久久| 免费高清在线观看视频在线观看| 日韩成人av中文字幕在线观看| 赤兔流量卡办理| 两性夫妻黄色片| 国产黄色视频一区二区在线观看| 日韩熟女老妇一区二区性免费视频| 久久影院123| 伊人久久国产一区二区| 男女下面插进去视频免费观看| 午夜福利视频精品| 国产精品麻豆人妻色哟哟久久| 人人妻人人澡人人看| 我要看黄色一级片免费的| 成人亚洲精品一区在线观看| 男女高潮啪啪啪动态图| 久久狼人影院| 又大又黄又爽视频免费| 欧美亚洲日本最大视频资源| 最近中文字幕2019免费版| 国产精品女同一区二区软件| 亚洲第一av免费看| 久久99精品国语久久久| 丰满饥渴人妻一区二区三| 丁香六月欧美| 女性生殖器流出的白浆| 热re99久久国产66热| 看免费成人av毛片| 我的亚洲天堂| 青春草视频在线免费观看| 免费不卡黄色视频| 日韩中文字幕欧美一区二区 | 男女边摸边吃奶| 久久青草综合色| 欧美老熟妇乱子伦牲交| 桃花免费在线播放| 9热在线视频观看99| 久久精品国产亚洲av高清一级| 国产片内射在线| 精品国产一区二区三区四区第35| 操出白浆在线播放| 午夜精品国产一区二区电影| a级毛片在线看网站| 免费日韩欧美在线观看| 19禁男女啪啪无遮挡网站| 国产精品女同一区二区软件| 一个人免费看片子| 这个男人来自地球电影免费观看 | 成人亚洲精品一区在线观看| 男女国产视频网站| 亚洲成人免费av在线播放| 少妇人妻久久综合中文| 在线观看www视频免费| 国产精品国产三级专区第一集| 大片免费播放器 马上看| 看免费成人av毛片| 国产亚洲欧美精品永久| 久久久久久人人人人人| 最近最新中文字幕大全免费视频 | 成年人免费黄色播放视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品麻豆人妻色哟哟久久| 亚洲av国产av综合av卡| 日韩 亚洲 欧美在线| 亚洲av日韩精品久久久久久密 | 又粗又硬又长又爽又黄的视频| 热re99久久精品国产66热6| 国产男女超爽视频在线观看| 色播在线永久视频| kizo精华| 亚洲国产精品一区三区| 日韩电影二区| 国产激情久久老熟女| 如何舔出高潮| av福利片在线| 秋霞在线观看毛片| 国产黄频视频在线观看| 中文乱码字字幕精品一区二区三区| 亚洲欧美精品自产自拍| 久久97久久精品| 美女高潮到喷水免费观看| 伦理电影免费视频| 宅男免费午夜| 国产一区亚洲一区在线观看| 久久国产精品男人的天堂亚洲| 超色免费av| 激情视频va一区二区三区| 国产av码专区亚洲av| 9色porny在线观看| 久久狼人影院| 极品人妻少妇av视频| 热re99久久精品国产66热6| 男女之事视频高清在线观看 | 欧美日韩亚洲国产一区二区在线观看 | 久久精品久久久久久噜噜老黄| 亚洲欧美清纯卡通| 男女免费视频国产| 亚洲国产欧美网| 丝袜美腿诱惑在线| 日日摸夜夜添夜夜爱| 女人久久www免费人成看片| www.熟女人妻精品国产| 久热这里只有精品99| 视频区图区小说| 日日啪夜夜爽| 最黄视频免费看| 丝瓜视频免费看黄片| 亚洲精品视频女| 丰满饥渴人妻一区二区三| 国产免费福利视频在线观看| 国产 精品1| 人人妻人人澡人人爽人人夜夜| 欧美乱码精品一区二区三区| 91精品三级在线观看| 日日爽夜夜爽网站| 亚洲国产精品999| xxxhd国产人妻xxx| 亚洲精华国产精华液的使用体验| 午夜免费鲁丝| 免费看不卡的av| 亚洲av电影在线观看一区二区三区| 国产一区有黄有色的免费视频| 午夜激情久久久久久久| 男女午夜视频在线观看| 国产男人的电影天堂91| 新久久久久国产一级毛片| 国产精品人妻久久久影院| 欧美日韩亚洲高清精品| 日日撸夜夜添| 久久久久久人妻| 色吧在线观看| 三上悠亚av全集在线观看| 菩萨蛮人人尽说江南好唐韦庄| 久久久精品区二区三区| 国产国语露脸激情在线看| 国产精品国产三级国产专区5o| 男人操女人黄网站| 亚洲欧美激情在线| 少妇被粗大猛烈的视频| 亚洲精品在线美女| 99热国产这里只有精品6| 一二三四在线观看免费中文在| 久久精品亚洲av国产电影网| av女优亚洲男人天堂| 亚洲人成电影观看| 女人爽到高潮嗷嗷叫在线视频| 久久99一区二区三区| 91精品国产国语对白视频| 叶爱在线成人免费视频播放| 久久久久久久精品精品| 美女脱内裤让男人舔精品视频| 久久综合国产亚洲精品| 蜜桃在线观看..| 大香蕉久久网| 成人三级做爰电影| 亚洲国产精品999| 日日摸夜夜添夜夜爱| 叶爱在线成人免费视频播放| 国产黄色视频一区二区在线观看| av免费观看日本| 国产乱人偷精品视频| 中文字幕色久视频| 少妇被粗大的猛进出69影院| 一区二区三区四区激情视频| 狠狠婷婷综合久久久久久88av| 9热在线视频观看99| 国产精品偷伦视频观看了| 18禁动态无遮挡网站| 女人精品久久久久毛片| 国产午夜精品一二区理论片| 日日摸夜夜添夜夜爱| 欧美日韩一区二区视频在线观看视频在线| 一级a爱视频在线免费观看| 国产 一区精品| 精品人妻一区二区三区麻豆| 欧美国产精品一级二级三级| 九草在线视频观看| 亚洲精品国产av成人精品| 天天操日日干夜夜撸| 性色av一级| 日韩精品有码人妻一区| 久久狼人影院| 久久99热这里只频精品6学生| 成人亚洲欧美一区二区av| 丝袜脚勾引网站| 黄片无遮挡物在线观看| 男女午夜视频在线观看| 9色porny在线观看| 午夜福利免费观看在线| 欧美日韩av久久| 日本色播在线视频| 精品人妻一区二区三区麻豆| 男女国产视频网站| 亚洲av日韩精品久久久久久密 | av女优亚洲男人天堂| 日本av免费视频播放| av在线app专区| 夫妻午夜视频| 女人久久www免费人成看片| 女性被躁到高潮视频| av在线观看视频网站免费| 午夜精品国产一区二区电影| 一级片'在线观看视频| 亚洲欧美成人精品一区二区| 国产免费一区二区三区四区乱码| 久久免费观看电影| 男人舔女人的私密视频| 亚洲精品自拍成人| 黑人巨大精品欧美一区二区蜜桃| 国产亚洲欧美精品永久| 亚洲精品美女久久久久99蜜臀 | 欧美日韩精品网址| 高清黄色对白视频在线免费看| 亚洲三区欧美一区| 久久久久精品久久久久真实原创| 狠狠精品人妻久久久久久综合| 啦啦啦在线观看免费高清www| 亚洲,欧美精品.| 黑人欧美特级aaaaaa片| 香蕉丝袜av| 欧美精品一区二区大全| 欧美国产精品va在线观看不卡| 一区二区三区精品91| 最近的中文字幕免费完整| 视频区图区小说| 777久久人妻少妇嫩草av网站| 交换朋友夫妻互换小说| 三上悠亚av全集在线观看| 一级片'在线观看视频| 国产一区二区三区综合在线观看| 老司机影院成人| 久久女婷五月综合色啪小说| 少妇人妻精品综合一区二区| 精品国产乱码久久久久久小说| 人人妻,人人澡人人爽秒播 | 亚洲国产av影院在线观看| 女的被弄到高潮叫床怎么办| 一边亲一边摸免费视频| 久久精品熟女亚洲av麻豆精品| 亚洲精品国产色婷婷电影| 亚洲人成77777在线视频| netflix在线观看网站| 亚洲精品日韩在线中文字幕| 九色亚洲精品在线播放| 国产一区二区在线观看av| 两个人免费观看高清视频| 久久久久人妻精品一区果冻| 精品国产乱码久久久久久男人| 9191精品国产免费久久| av在线观看视频网站免费| 女人高潮潮喷娇喘18禁视频| 少妇的丰满在线观看| 在线精品无人区一区二区三| 免费观看性生交大片5| 成人毛片60女人毛片免费| 我要看黄色一级片免费的| 在线看a的网站| 国产一级毛片在线| 操出白浆在线播放| 赤兔流量卡办理| 亚洲婷婷狠狠爱综合网| 90打野战视频偷拍视频| 久久久久国产精品人妻一区二区| 色精品久久人妻99蜜桃| 亚洲av成人不卡在线观看播放网 | 亚洲精品自拍成人| 久久久久国产精品人妻一区二区| 日韩一卡2卡3卡4卡2021年| 老熟女久久久| 丰满迷人的少妇在线观看| 人人妻人人澡人人爽人人夜夜| 日韩免费高清中文字幕av| 午夜av观看不卡| 国产精品久久久久久精品古装| 日本av免费视频播放| 亚洲专区中文字幕在线 | 亚洲精品成人av观看孕妇| 只有这里有精品99| 天天影视国产精品| 一级片'在线观看视频| 亚洲av日韩精品久久久久久密 | 欧美精品一区二区大全| 9热在线视频观看99| 日韩人妻精品一区2区三区| 纯流量卡能插随身wifi吗| 男女国产视频网站| 你懂的网址亚洲精品在线观看| 又黄又粗又硬又大视频| 成人黄色视频免费在线看| 精品一品国产午夜福利视频| 啦啦啦在线观看免费高清www| 国产av一区二区精品久久| 女人爽到高潮嗷嗷叫在线视频| 久久国产精品大桥未久av| 久久久久久免费高清国产稀缺| av网站在线播放免费| 亚洲国产看品久久| 亚洲国产欧美网| 中文欧美无线码| 老司机靠b影院| 国产精品欧美亚洲77777| 最近中文字幕2019免费版| www日本在线高清视频| 成人黄色视频免费在线看| 两个人免费观看高清视频| 久久ye,这里只有精品| 国产成人精品久久久久久| 自线自在国产av| 亚洲av男天堂| 亚洲国产日韩一区二区| 亚洲精品第二区| 99热网站在线观看| 日韩一区二区视频免费看| 女人爽到高潮嗷嗷叫在线视频| 久久99精品国语久久久|