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

    計(jì)及風(fēng)電-負(fù)荷不確定性的風(fēng)-火-核-碳捕集多源協(xié)調(diào)優(yōu)化調(diào)度

    2022-02-21 07:49:06趙冬梅王浩翔
    電工技術(shù)學(xué)報(bào) 2022年3期
    關(guān)鍵詞:存儲(chǔ)單元核電機(jī)組調(diào)峰

    趙冬梅 王浩翔 陶 然

    計(jì)及風(fēng)電-負(fù)荷不確定性的風(fēng)-火-核-碳捕集多源協(xié)調(diào)優(yōu)化調(diào)度

    趙冬梅 王浩翔 陶 然

    (華北電力大學(xué)電氣與電子工程學(xué)院 北京 102206)

    隨著核電技術(shù)的飛速發(fā)展,大容量核電機(jī)組接入電網(wǎng)是實(shí)現(xiàn)電力低碳環(huán)保的必然趨勢(shì),而負(fù)荷峰谷差的增大和間歇式電源滲透率的提高,迫切需要核電機(jī)組以靈活運(yùn)行方式分擔(dān)電網(wǎng)調(diào)峰壓力。首先,基于核電機(jī)組實(shí)際運(yùn)行特點(diǎn),線性化核電調(diào)峰深度;然后,分析綜合靈活運(yùn)行方式碳捕集電廠的“削峰填谷”特性,實(shí)現(xiàn)核電-碳捕集機(jī)組捆綁調(diào)峰;接著,引入模糊參數(shù)表征風(fēng)電和負(fù)荷的不確定性,在計(jì)及核電調(diào)峰安全性和經(jīng)濟(jì)性的前提下,以各類機(jī)組聯(lián)合總運(yùn)行成本最低為目標(biāo)函數(shù),建立考慮風(fēng)電消納效益的風(fēng)、火、核、碳捕集多源協(xié)調(diào)模糊優(yōu)化調(diào)度模型;最后,通過仿真算例驗(yàn)證所提模型和方法的有效性。結(jié)果表明,該優(yōu)化調(diào)度模型在保證核電安全運(yùn)行的基礎(chǔ)上,提高了電網(wǎng)調(diào)度柔性,實(shí)現(xiàn)了運(yùn)行的經(jīng)濟(jì)性和低碳性。

    核電機(jī)組 碳捕集機(jī)組 調(diào)峰 削峰填谷 優(yōu)化調(diào)度

    0 引言

    隨著經(jīng)濟(jì)社會(huì)的不斷進(jìn)步,溫室效應(yīng)與能源問題日益凸顯,大力發(fā)展核電、風(fēng)電、光伏等清潔能源成為解決環(huán)境問題和能源危機(jī)的重要方式[1-2]。但風(fēng)電、光伏具有間歇性和強(qiáng)波動(dòng)性,目前仍處在小規(guī)模的應(yīng)用階段,而核電因其低碳性與環(huán)保性在電力結(jié)構(gòu)中的比重逐年增大[3]。截至2019年底,全球共有443臺(tái)在運(yùn)核電機(jī)組,總凈裝機(jī)容量為392GWe,發(fā)電量占全球發(fā)電量約10.4%[4]。在核電占比較高地區(qū),核電站因保障其安全性要求以基本負(fù)荷運(yùn)行為主,年利用小時(shí)長(zhǎng),運(yùn)行負(fù)荷率高,基本不參與電網(wǎng)調(diào)峰,對(duì)電網(wǎng)安全運(yùn)行影響的負(fù)面性隨著核電裝機(jī)容量增大而上升,這無疑加劇了電網(wǎng)的調(diào)峰壓力,加重了燃煤機(jī)組深度調(diào)峰任務(wù),對(duì)電網(wǎng)的安全性和經(jīng)濟(jì)性均造成重大影響。因此,考慮到電網(wǎng)調(diào)峰形勢(shì)越加嚴(yán)峻[5-6],使核電積極參與電網(wǎng)日調(diào)峰的需求變得越發(fā)強(qiáng)烈。

    目前有關(guān)核電調(diào)峰的研究成果較少,且多以多電源聯(lián)合運(yùn)行緩解電網(wǎng)的調(diào)峰需求為主[7-9]。文獻(xiàn)[10]介紹了核電機(jī)組參與電網(wǎng)調(diào)峰的運(yùn)行方式。文獻(xiàn)[11]在考慮核電機(jī)組出力特性的前提下,將核電設(shè)置為三擋固定調(diào)峰深度,制定了核電和抽水蓄能聯(lián)合調(diào)峰調(diào)度策略。文獻(xiàn)[12]針對(duì)高比例風(fēng)電的反調(diào)峰特性,構(gòu)建了風(fēng)-火-核多源協(xié)調(diào)優(yōu)化調(diào)度模型以緩解風(fēng)電消納困難問題。文獻(xiàn)[13]提出計(jì)及核電安全可調(diào)節(jié)域和調(diào)峰安全約束的兩階段調(diào)度決策方法,對(duì)比分析了兩種典型場(chǎng)景下各機(jī)組的開停機(jī)及出力情況。上述文獻(xiàn)所建立的核電調(diào)峰數(shù)學(xué)模型相對(duì)簡(jiǎn)單,核電機(jī)組調(diào)峰深度固定,運(yùn)行方式不夠靈活。同時(shí),文獻(xiàn)[11]中,由于抽水蓄能電站受到地理位置的制約,其與位于負(fù)荷中心的核電站捆綁調(diào)峰存在一定的限制。所以,探索一種易于和核電進(jìn)行調(diào)峰配合且工況轉(zhuǎn)換靈活,相較其他電源在負(fù)荷跟蹤能力上具有明顯優(yōu)勢(shì)的調(diào)峰電源受到了專家學(xué)者的廣泛關(guān)注。

    近年來,我國(guó)正在積極調(diào)整電力結(jié)構(gòu),但火電作為最大的碳排放源,仍在我國(guó)電力結(jié)構(gòu)中占據(jù)主導(dǎo)地位。為實(shí)現(xiàn)電力低碳化,將常規(guī)燃煤電廠改造成碳捕集電廠可以有效降低機(jī)組的碳排放量[14],避免負(fù)荷峰谷差過大。此外,碳捕集電廠與傳統(tǒng)燃煤電廠相比,具有良好的功率調(diào)整特性,可以與不同類型電源協(xié)調(diào)運(yùn)行[15]。目前有關(guān)碳捕集電廠配合其他電源參與電網(wǎng)調(diào)峰的研究較少。文獻(xiàn)[16]分析了碳捕集電廠靈活運(yùn)行的基本原理,指出碳捕集電廠具有較好的負(fù)荷跟蹤能力和更深的調(diào)峰深度,但只定性闡述了其靈活運(yùn)行方式的實(shí)現(xiàn)方法和技術(shù)特性。文獻(xiàn)[17]定義碳捕集電廠的電碳特性,建立了靈活運(yùn)行方式碳捕集電廠的數(shù)學(xué)模型,但靈活運(yùn)行方式?jīng)]有涉及電力系統(tǒng)層面。文獻(xiàn)[18]在改進(jìn)傳統(tǒng)單碳量模型的基礎(chǔ)上引入雙碳量模型,建立包含虛擬電廠的經(jīng)濟(jì)調(diào)度模型,實(shí)現(xiàn)了碳捕集系統(tǒng)吸收環(huán)節(jié)與解析環(huán)節(jié)的解耦。文獻(xiàn)[19]提出了風(fēng)-光-水-碳捕集多區(qū)域虛擬電廠協(xié)調(diào)優(yōu)化調(diào)度模型,通過優(yōu)化不同類型機(jī)組出力時(shí)序,有效降低了電網(wǎng)的碳排放強(qiáng)度。上述文獻(xiàn)對(duì)碳捕集電廠的運(yùn)行原理、經(jīng)濟(jì)調(diào)度進(jìn)行了深入分析,但對(duì)于碳捕集電廠的調(diào)峰能力及“削峰填谷”特性缺少理論研究。

    由此可見,以核電機(jī)組和碳捕集機(jī)組等低碳和環(huán)境友好型電源代替碳排放強(qiáng)度較大的燃煤機(jī)組作為可控發(fā)電單元是未來的重點(diǎn)研究方向。同時(shí),由于風(fēng)電、光伏等清潔能源在電網(wǎng)中的滲透率不斷提高,其預(yù)測(cè)出力的不確定性使模型計(jì)算結(jié)果與實(shí)際值之間存在偏差[20-21],分析間歇性電源和負(fù)荷預(yù)測(cè)的不確定性顯得尤為重要。

    綜上所述,本文在現(xiàn)有研究成果的基礎(chǔ)上,為進(jìn)一步降低電網(wǎng)的運(yùn)行成本與碳排放量,充分挖掘核電機(jī)組、碳捕集機(jī)組與可再生能源聯(lián)合運(yùn)行的潛力,構(gòu)建了考慮風(fēng)電消納效益的風(fēng)、火、核、碳捕集多源協(xié)調(diào)的模糊優(yōu)化調(diào)度模型。首先,分析核電機(jī)組的可調(diào)度性,線性化核電調(diào)峰深度,使核電調(diào)峰模型更加精確化;其次,建立綜合靈活運(yùn)行方式碳捕集電廠的數(shù)學(xué)模型,對(duì)碳捕集電廠的“削峰填谷”特性進(jìn)行理論分析;接著引入模糊參數(shù)表征風(fēng)電和負(fù)荷的不確定性,對(duì)模糊期望約束和模糊機(jī)會(huì)約束進(jìn)行等效處理;最后,通過算例分析風(fēng)、核、碳捕集協(xié)調(diào)調(diào)度的經(jīng)濟(jì)性和低碳性,驗(yàn)證了本文模型的適用性和有效性。

    1 核電機(jī)組調(diào)峰能力分析及其數(shù)學(xué)模型

    1.1 核電機(jī)組調(diào)峰能力分析

    核電在滿足自身安全、經(jīng)濟(jì)的前提下參與電網(wǎng)調(diào)峰一般分為日負(fù)荷跟蹤運(yùn)行、長(zhǎng)期降功率運(yùn)行和季節(jié)性調(diào)峰三種運(yùn)行方式,但受電價(jià)機(jī)制影響,主要采用日負(fù)荷跟蹤運(yùn)行方式。在設(shè)計(jì)上,我國(guó)目前已投運(yùn)核電機(jī)組均具備日負(fù)荷跟蹤能力。以CNP650為代表的二代核電機(jī)組可以在30%~100%N功率調(diào)整范圍內(nèi)以每分鐘1%~3%N的速度實(shí)現(xiàn)功率調(diào)節(jié)。對(duì)于大力發(fā)展的三代核電機(jī)組,其調(diào)節(jié)范圍更寬,最低下調(diào)峰深度可達(dá)額定功率的85%,具有每分鐘5%N最大功率調(diào)節(jié)速度和10%階躍功率變化幅度,日負(fù)荷跟蹤能力更加優(yōu)越,可以滿足電網(wǎng)的調(diào)峰需求。

    在實(shí)際運(yùn)行方面,圖1給出了法國(guó)Golfech核電廠1號(hào)機(jī)組參與調(diào)峰運(yùn)行數(shù)據(jù)[22]。在一個(gè)換料周期內(nèi),核電機(jī)組進(jìn)行了155次負(fù)荷跟蹤操作,最低功率水平從燃料循環(huán)的65%之后隨著燃耗的加深逐漸提高至85%N,在燃料循環(huán)的70%之后,反應(yīng)堆的靈活運(yùn)行能力有所下降,基本不參與調(diào)峰。

    另外,核電機(jī)組參與調(diào)峰,存在導(dǎo)致機(jī)組運(yùn)行可靠性降低的風(fēng)險(xiǎn)(如功率調(diào)節(jié)棒動(dòng)作頻繁、燃料包殼破損概率增大等問題),對(duì)運(yùn)行人員的操作能力和實(shí)踐經(jīng)驗(yàn)有更高的要求。因此,需要從滿足安全性的角度建立核電日調(diào)峰數(shù)學(xué)模型。

    1.2 核電機(jī)組調(diào)峰數(shù)學(xué)模型

    圖2 核電機(jī)組調(diào)峰運(yùn)行原理圖

    核電機(jī)組的運(yùn)行約束為

    1)最大調(diào)峰深度約束

    2)狀態(tài)標(biāo)志耦合約束

    當(dāng)功率過渡時(shí)間為3h時(shí),有

    當(dāng)功率過渡時(shí)間為2h時(shí),有

    功率線性上升和線性下降時(shí)間相等約束為

    3)額定運(yùn)行狀態(tài)和低功率運(yùn)行狀態(tài)時(shí)間約束

    需要注意的是,國(guó)內(nèi)大部分核電機(jī)組的技術(shù)源自法國(guó)核電機(jī)組采用的技術(shù)。根據(jù)法國(guó)核電豐富的調(diào)峰經(jīng)驗(yàn)及目前我國(guó)核電機(jī)組的調(diào)峰現(xiàn)狀,尚沒有研究和報(bào)告表明核電機(jī)組頻繁參與調(diào)峰會(huì)導(dǎo)致放射性排放顯著增加,即核電調(diào)峰可以滿足放射性要求。

    2 碳捕集電廠數(shù)學(xué)模型及其“削峰填谷”特性分析

    碳捕集電廠的靈活運(yùn)行方式包括煙氣分流、溶液存儲(chǔ)及綜合靈活運(yùn)行方式。以煙氣分流方式為基礎(chǔ),在吸收塔和解析塔之間配置儲(chǔ)液?jiǎn)卧礊榫C合靈活運(yùn)行方式。

    碳捕集電廠綜合靈活運(yùn)行方式運(yùn)行原理如圖3a所示[18]。該流程實(shí)現(xiàn)了發(fā)電循環(huán)和吸收環(huán)節(jié)及CO2吸收與解析的解耦。

    圖3 碳捕集電廠運(yùn)行原理及其功率流

    2.1 碳捕集電廠數(shù)學(xué)模型

    由圖3b可知,電廠對(duì)外表現(xiàn)出的輸出功率(凈輸出功率)為總輸出功率扣除碳捕集系統(tǒng)的能量損耗,即

    根據(jù)圖3b,碳捕集系統(tǒng)CO2捕集量為吸收塔和存儲(chǔ)單元共同提供的CO2質(zhì)量,即

    由于吸收環(huán)節(jié)的能量損耗只占碳捕集損耗的2%~10%,本文忽略其影響,則碳捕集系統(tǒng)的運(yùn)行損耗可以表示為

    由此,得到電廠的凈輸出功率為

    由式(15)可知:加入存儲(chǔ)單元后,吸收與解析環(huán)節(jié)相互獨(dú)立,碳捕集電廠可以通過調(diào)整對(duì)碳捕集系統(tǒng)的能量供應(yīng)來迅速改變凈輸出功率。

    碳捕集電廠的運(yùn)行約束為

    此外,碳捕集電廠出力上下限約束、最小起停時(shí)間約束和爬坡約束等與燃煤機(jī)組類似,本文不再贅述。

    2.2 碳捕集電廠“削峰填谷”特性分析

    碳捕集電廠功率運(yùn)行區(qū)間如圖4所示。傳統(tǒng)燃煤電廠由于沒有加裝碳捕集系統(tǒng),運(yùn)行方式不夠靈活,凈輸出功率運(yùn)行區(qū)間最小。采用綜合靈活運(yùn)行方式后,碳捕集電廠的凈輸出功率運(yùn)行區(qū)間最大,當(dāng)富液存儲(chǔ)單元儲(chǔ)存的CO2溶液接近飽和時(shí),解析塔工作在最大運(yùn)行狀態(tài),碳捕集系統(tǒng)運(yùn)行損耗最大,凈輸出功率最小;當(dāng)富液存儲(chǔ)單元儲(chǔ)存的CO2溶液逐漸減少為零時(shí),碳捕集系統(tǒng)運(yùn)行損耗逐漸降低,最小凈輸出功率隨之提高,最終與煙氣分流式碳捕集電廠的最小凈輸出功率相等[24]。

    圖4 碳捕集電廠功率運(yùn)行區(qū)間

    與傳統(tǒng)燃煤電廠相比,碳捕集電廠同時(shí)引入煙氣旁路系統(tǒng)和溶液存儲(chǔ)單元,擴(kuò)大了電廠的凈輸出功率運(yùn)行范圍,從調(diào)度側(cè)看,等效于降低了電廠的出力下限,賦予其更大的調(diào)峰深度。

    另外,在負(fù)荷高峰時(shí)段,碳捕集機(jī)組跟蹤負(fù)荷需提高凈輸出功率,CO2排放量較大,此時(shí)若利用碳捕集系統(tǒng)捕集CO2,會(huì)降低機(jī)組凈輸出功率,使上調(diào)峰能力不足。在引入溶液存儲(chǔ)單元后,富液存儲(chǔ)單元儲(chǔ)碳但不釋放,減少碳排放量的同時(shí)增加了凈輸出功率(存儲(chǔ)單元釋放能量)。在負(fù)荷低谷時(shí)段,碳捕集機(jī)組跟蹤負(fù)荷需降低凈輸出功率,此時(shí)富液存儲(chǔ)單元釋放CO2增大碳捕集系統(tǒng)損耗,以此來減小凈輸出功率(存儲(chǔ)單元儲(chǔ)存能量),實(shí)現(xiàn)碳捕集電廠“削峰填谷”功能。

    3 考慮風(fēng)電消納效益的風(fēng)、火、核、碳捕集多源協(xié)調(diào)的模糊優(yōu)化調(diào)度模型

    對(duì)于風(fēng)、火、核、碳捕集多電源共同參與的電力系統(tǒng)優(yōu)化調(diào)度問題,應(yīng)綜合考慮各電源的運(yùn)行約束和系統(tǒng)的安全約束。假設(shè)存在一個(gè)中央集權(quán)的管理者,能夠以整體社會(huì)效益最大化為目標(biāo)建立調(diào)度決策方案,統(tǒng)一管轄各類發(fā)電廠,本文在經(jīng)濟(jì)調(diào)度的基礎(chǔ)上引入碳交易機(jī)制,同時(shí)考慮源荷雙側(cè)的不確定性,以多電源運(yùn)行成本最小為目標(biāo),通過優(yōu)化方法確定各電源的出力模式,實(shí)現(xiàn)電網(wǎng)運(yùn)行的低碳性和經(jīng)濟(jì)性。

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

    考慮到實(shí)際系統(tǒng)中源-荷雙側(cè)數(shù)據(jù)信息獲取不充分、收集困難的情況,本文采用模糊數(shù)來表征風(fēng)電出力和負(fù)荷的不確定性,并用三角隸屬度函數(shù)描述其不確定性程度。受限于篇幅,對(duì)應(yīng)的風(fēng)電和負(fù)荷模糊數(shù)表征方法參見文獻(xiàn)[25]。

    在日調(diào)峰優(yōu)化模型中,當(dāng)電網(wǎng)調(diào)峰能力不足時(shí),考慮將核電機(jī)組滿功率運(yùn)行方式轉(zhuǎn)換為低功率運(yùn)行方式參與電網(wǎng)調(diào)峰。以風(fēng)電機(jī)組、燃煤機(jī)組、核電機(jī)組及碳捕集機(jī)組多源系統(tǒng)總運(yùn)行成本費(fèi)用最小為目標(biāo)函數(shù),即

    1)燃煤機(jī)組運(yùn)行成本

    燃煤發(fā)電成本為

    2)碳捕集機(jī)組運(yùn)行成本

    碳捕集機(jī)組發(fā)電運(yùn)行成本為

    設(shè)備折舊成本為

    3)核電機(jī)組運(yùn)行成本

    核電發(fā)電成本為

    核電調(diào)峰成本為

    4)風(fēng)電棄風(fēng)成本

    5)碳交易成本

    核電和風(fēng)電均為清潔能源,系統(tǒng)的碳排放主要來源于燃煤機(jī)組和碳捕集機(jī)組,則全天總碳排放量為

    碳交易成本為

    3.2 約束條件

    1)功率平衡方程

    2)正負(fù)旋轉(zhuǎn)備用約束

    正負(fù)旋轉(zhuǎn)備用約束參考文獻(xiàn)[25],具體不再贅述。

    3)燃煤機(jī)組約束

    燃煤機(jī)組約束包括常規(guī)機(jī)組出力上下限、爬坡和最小開停機(jī)時(shí)間約束,本文不再給出。

    4)核電機(jī)組約束

    核電機(jī)組約束由式(1)~式(9)表示。

    5)碳捕集機(jī)組約束

    碳捕集組約束由式(10)~式(17)表示。

    6)網(wǎng)絡(luò)潮流約束

    3.3 模型求解

    本文將考慮風(fēng)電消納效益的風(fēng)、火、核、碳捕集多源協(xié)調(diào)的模糊優(yōu)化調(diào)度模型中功率平衡方程和碳交易成本采用模糊期望約束表示,將正負(fù)旋轉(zhuǎn)備用和棄風(fēng)悲觀值采用模糊機(jī)會(huì)約束表示,則功率平衡方程通過模糊期望計(jì)算公式可以轉(zhuǎn)換為確定性等價(jià)模型,即

    碳交易成本經(jīng)轉(zhuǎn)換后的確定性等價(jià)模型為

    棄風(fēng)悲觀值、正負(fù)旋轉(zhuǎn)備用經(jīng)轉(zhuǎn)換后的確定性等價(jià)模型類似于文獻(xiàn)[25, 29]。轉(zhuǎn)換后的確定性等價(jià)模型可以在YALMIP環(huán)境下借助商用優(yōu)化軟件Gurobi高效求解。模糊優(yōu)化調(diào)度模型的求解流程如圖5所示。

    4 算例分析

    為驗(yàn)證本文所提風(fēng)、火、核、碳捕集多源協(xié)調(diào)模糊優(yōu)化調(diào)度模型的有效性,選取改進(jìn)的IEEE-39節(jié)點(diǎn)測(cè)試系統(tǒng)進(jìn)行算例分析。

    4.1 算例基礎(chǔ)數(shù)據(jù)

    根據(jù)是否引入碳捕集機(jī)組及是否考慮風(fēng)電和負(fù)荷的不確定性設(shè)置以下四個(gè)場(chǎng)景進(jìn)行對(duì)比:

    場(chǎng)景1:不考慮碳捕集電廠,核電機(jī)組參與調(diào)峰,不計(jì)及風(fēng)電出力和負(fù)荷的不確定性。

    場(chǎng)景2:考慮分流式碳捕集電廠,核電機(jī)組參與調(diào)峰,不計(jì)及風(fēng)電出力和負(fù)荷的不確定性。

    場(chǎng)景3:考慮綜合靈活運(yùn)行方式碳捕集電廠(含溶液存儲(chǔ)單元),核電機(jī)組參與調(diào)峰,不計(jì)及風(fēng)電出力和負(fù)荷的不確定性。

    場(chǎng)景4:考慮綜合靈活運(yùn)行方式碳捕集電廠,核電機(jī)組參與調(diào)峰,計(jì)及風(fēng)電出力和負(fù)荷的不確定性。

    4.2 算例結(jié)果分析

    圖6、圖7給出了核電采用日負(fù)荷跟蹤運(yùn)行方式,不同場(chǎng)景下所有機(jī)組和核電機(jī)組出力情況。圖6中四個(gè)場(chǎng)景的不同機(jī)組出力從下至上依次為燃煤機(jī)組4、燃煤機(jī)組5、燃煤機(jī)組6、燃煤機(jī)組1、燃煤機(jī)組2、燃煤機(jī)組3、碳捕集電廠1、碳捕集電廠2、核電機(jī)組1、核電機(jī)組2、風(fēng)電。

    圖6 不同場(chǎng)景所有機(jī)組出力

    圖7 不同場(chǎng)景核電機(jī)組出力

    由圖6a、圖7a可知:在負(fù)荷低谷時(shí)期,核電機(jī)組1、2分別在2~8h與1~7h參與電網(wǎng)調(diào)峰,調(diào)峰深度為15.6%和21.6%,共計(jì)調(diào)峰量為2 451.9MW?h。與調(diào)峰深度采用固定檔位運(yùn)行方式相比,線性化調(diào)峰深度將核電低功率狀態(tài)刻畫得更加準(zhǔn)確,避免了調(diào)峰深度不精確導(dǎo)致核電調(diào)峰成本的增加,但部分燃煤機(jī)組以起停方式參與電網(wǎng)調(diào)峰。為進(jìn)一步提高運(yùn)行經(jīng)濟(jì)性,降低電網(wǎng)碳排放水平,場(chǎng)景2中將部分燃煤電廠改造為分流式碳捕集電廠,并考慮碳捕集電廠的運(yùn)行機(jī)制,圖6b、圖7b為引入分流式碳捕集電廠后所有機(jī)組和核電機(jī)組出力情況。

    分析圖6b、圖7b可知:相比于場(chǎng)景1,引入分流式碳捕集電廠后,核電機(jī)組1、2在2~8h和1~7h參與調(diào)峰的調(diào)峰深度為10.0%和17.8%,調(diào)峰深度減小,原因在于負(fù)荷低谷時(shí)段,分流式碳捕集電廠煙氣旁路系統(tǒng)通過控制流入吸收塔的CO2量,調(diào)節(jié)碳捕集能耗,降低機(jī)組出力,使碳捕集電廠具有較大的調(diào)峰深度,從而減小核電機(jī)組的調(diào)峰壓力,提高核電機(jī)組調(diào)峰的安全性。同時(shí),電網(wǎng)中起停調(diào)峰燃煤機(jī)組數(shù)量為零,棄風(fēng)電量降低。分流式碳捕集電廠的引入,提高了電網(wǎng)運(yùn)行的經(jīng)濟(jì)性和調(diào)度柔性,降低了電網(wǎng)碳排放量(后文介紹)。但分流式碳捕集電廠的調(diào)峰能力有限,因此場(chǎng)景3中在分流式碳捕集電廠基礎(chǔ)上增加溶液存儲(chǔ)單元,并考慮其綜合靈活運(yùn)行方式,圖6c、圖7c為引入綜合靈活運(yùn)行方式碳捕集電廠后所有機(jī)組和核電機(jī)組出力。

    由圖6c、圖7c可知:引入綜合靈活運(yùn)行方式碳捕集電廠后,僅核電機(jī)組1參與電網(wǎng)調(diào)峰,核電機(jī)組2滿出力運(yùn)行,減少的調(diào)峰電量由碳捕集電廠承擔(dān)。與場(chǎng)景2相比,場(chǎng)景3核電調(diào)峰電量更少,為1 260MW?h,這是因?yàn)閳?chǎng)景2碳捕集損耗在負(fù)荷低谷時(shí)段較小,在負(fù)荷高峰時(shí)段較大,而場(chǎng)景3碳捕集機(jī)組運(yùn)行更加靈活,在負(fù)荷低谷時(shí)段煙氣旁路系統(tǒng)提供的CO2與富液存儲(chǔ)單元增大提供的CO2相結(jié)合,使碳捕集設(shè)備損耗增大,凈輸出功率更小(相當(dāng)于存儲(chǔ)單元儲(chǔ)存能量即“填谷”),從而實(shí)現(xiàn)輔助核電機(jī)組調(diào)峰的任務(wù)。在負(fù)荷高峰時(shí)段,富液存儲(chǔ)單元減小所提供的CO2量,碳捕集設(shè)備損耗降低,凈輸出功率提高(相當(dāng)于存儲(chǔ)單元釋放能量即“削峰”),減輕了其他機(jī)組的上調(diào)峰壓力。同時(shí),場(chǎng)景3運(yùn)行成本相比場(chǎng)景2下降1.5%。碳捕集機(jī)組和燃煤機(jī)組總出力如圖8所示。分析可知:綜合靈活運(yùn)行方式碳捕集電廠可以降低燃煤機(jī)組的出力峰谷差,大小為418MW,其“削峰填谷”機(jī)制的有效性得到驗(yàn)證。

    圖8 兩場(chǎng)景碳捕集機(jī)組和燃煤機(jī)組總輸出功率

    圖9為場(chǎng)景3碳捕集電廠中富液存儲(chǔ)單元與貧液存儲(chǔ)單元總體積變化量。由圖9可知:富液存儲(chǔ)單元在負(fù)荷低谷時(shí)段釋放CO2來增大碳捕集設(shè)備的損耗,即富液存儲(chǔ)單元液體體積下降、貧液存儲(chǔ)單元液體體積上升;在負(fù)荷高峰時(shí)段儲(chǔ)存CO2來減小碳捕集設(shè)備損耗,即富液存儲(chǔ)單元液體體積上升、貧液存儲(chǔ)單元液體體積下降,印證了碳捕集電廠的“削峰填谷”特性。

    圖9 溶液存儲(chǔ)單元總?cè)芤后w積

    在場(chǎng)景4中,所有機(jī)組出力如圖6d所示,核電出力情況如圖7d所示。由圖6d、圖7d可知:通過采用模糊期望約束和模糊機(jī)會(huì)約束來考慮風(fēng)電預(yù)測(cè)出力和負(fù)荷預(yù)測(cè)值的不確定性后,電網(wǎng)總運(yùn)行成本降低1.1%,且棄風(fēng)電量為零。不同場(chǎng)景下電網(wǎng)風(fēng)電調(diào)度情況如圖10所示。圖中,場(chǎng)景4風(fēng)電調(diào)度出力最大,而場(chǎng)景1風(fēng)電調(diào)度出力最小,證明本文所提模型可以實(shí)現(xiàn)風(fēng)電消納的最大化。

    圖10 不同場(chǎng)景風(fēng)電調(diào)度出力

    四種場(chǎng)景成本比較見表1。由表1可知:場(chǎng)景2相比于場(chǎng)景1,核電機(jī)組日負(fù)荷跟蹤運(yùn)行時(shí)調(diào)峰成本下降,雖然出現(xiàn)折舊成本,但由于其棄風(fēng)成本小、燃煤機(jī)組起停調(diào)峰數(shù)量少,使得總運(yùn)行成本小于場(chǎng)景1;場(chǎng)景3相比于場(chǎng)景2,碳排放量降低14.7%,雖然總折舊成本增大,但存儲(chǔ)單元的折舊成本較低,其帶來的效益可以抵消增大的折舊成本,使得總運(yùn)行成本下降,體現(xiàn)了綜合靈活運(yùn)行方式碳捕集機(jī)組的經(jīng)濟(jì)性和低碳性。

    表1 各場(chǎng)景成本比較

    Tab.1 Cost comparison of each scenario

    5 結(jié)論

    本文在核電參與日調(diào)峰的基礎(chǔ)上,將常規(guī)燃煤電廠改造成為碳捕集電廠,并計(jì)及風(fēng)電預(yù)測(cè)出力和負(fù)荷預(yù)測(cè)值的不確定性,構(gòu)建了考慮風(fēng)電消納效益的風(fēng)、火、核、碳捕集多源協(xié)調(diào)的模糊優(yōu)化調(diào)度模型,通過測(cè)試算例驗(yàn)證了所提模型和方法的有效性和適用性。得到如下結(jié)論:

    1)線性化調(diào)峰深度后,核電機(jī)組調(diào)峰電量更加精確。同時(shí),核電以日負(fù)荷跟蹤運(yùn)行方式參與電網(wǎng)調(diào)峰,有效緩解了燃煤機(jī)組調(diào)峰壓力,提高了電網(wǎng)運(yùn)行的經(jīng)濟(jì)性。

    2)含有存儲(chǔ)單元的分流式碳捕集機(jī)組以綜合靈活運(yùn)行方式參與電網(wǎng)運(yùn)行,可以代替抽水蓄能電站實(shí)現(xiàn)“削峰填谷”,并輔助核電機(jī)組完成調(diào)峰任務(wù)。相比于不含溶液存儲(chǔ)單元情況,電網(wǎng)總運(yùn)行成本下降1.5%,碳排放量下降14.7%,驗(yàn)證了綜合靈活運(yùn)行碳捕集機(jī)組的經(jīng)濟(jì)性和低碳性。

    3)考慮風(fēng)電和負(fù)荷的不確定性后,可以實(shí)現(xiàn)風(fēng)電的完全消納,且總運(yùn)行成本相比確定性降低1.1%。

    本文所建立的核電機(jī)組連續(xù)調(diào)峰模型并未深入研究放射性污染物的增加量。雖然核電調(diào)峰可以滿足放射性要求,不會(huì)導(dǎo)致放射性排放顯著增加,但對(duì)調(diào)峰運(yùn)行所產(chǎn)生的放射性廢物進(jìn)行量化仍是至關(guān)重要的。下一步研究將討論核電機(jī)組連續(xù)調(diào)峰運(yùn)行對(duì)放射性三廢、維護(hù)維修成本等因素的影響。

    附錄 仿真系統(tǒng)參數(shù)

    附圖1 改進(jìn)的IEEE 39節(jié)點(diǎn)拓?fù)鋱D

    App.Fig.1 Improved IEEE 39 node topology

    附圖2 負(fù)荷功率和風(fēng)電功率預(yù)測(cè)曲線

    App.Fig.2 Load and wind forecast profile

    附表1 風(fēng)電和負(fù)荷三角模糊參數(shù)

    App.Tab.1 Triangular fuzzy parameters

    模糊參數(shù)風(fēng)電負(fù)荷 數(shù)值0.51.50.91.1

    [1] 李爭(zhēng), 張蕊, 孫鶴旭, 等. 可再生能源多能互補(bǔ)制-儲(chǔ)-運(yùn)氫關(guān)鍵技術(shù)綜述[J]. 電工技術(shù)學(xué)報(bào), 2021, 36(3): 446-462.

    Li Zheng, Zhang Rui, Sun Hexu, et al. Review on key technologies of hydrogen generation, storage and transportation based on multi-energy complementary renewable energy[J]. Transactions of China Electrotechnical Society, 2021, 36(3): 446-462.

    [2] 孫鶴旭, 李爭(zhēng), 陳愛兵, 等. 風(fēng)電制氫技術(shù)現(xiàn)狀及發(fā)展趨勢(shì)[J]. 電工技術(shù)學(xué)報(bào), 2019, 34(19): 4071-4083.

    Sun Hexu, Li Zheng, Chen Aibing, et al. Current status and development trend of hydrogen production technology by wind power[J]. Transactions of China Electrotechnical Society, 2019, 34(19): 4071-4083.

    [3] 宋翔宇. 世界核電發(fā)展現(xiàn)狀[J]. 中國(guó)核電, 2017, 10(3): 439-443.

    Song Xiangyu. World nuclear power status[J]. China Nuclear Power, 2017, 10(3): 439-443.

    [4] 王政, 伍浩松. IAEA公布2019年全球核電數(shù)據(jù)[J]. 國(guó)外核新聞, 2020(8): 22-24.

    Wang Zheng, Wu Haosong. IAEA releases global nuclear power data for 2019[J]. Foreign Nuclear News, 2020(8): 22-24.

    [5] 羅桓桓, 王昊, 葛維春, 等. 考慮報(bào)價(jià)監(jiān)管的動(dòng)態(tài)調(diào)峰輔助服務(wù)市場(chǎng)競(jìng)價(jià)機(jī)制設(shè)計(jì)[J]. 電工技術(shù)學(xué)報(bào), 2021, 36(9): 1935-1947.

    Luo Huanhuan, Wang Hao, Ge Weichun, et al. Design of dynamic peak regulation ancillary service market bidding mechanism considering quotation supervision[J]. Transactions of China Electrotechnical Society, 2021, 36(9): 1935-1947.

    [6] 麻秀范, 王戈, 朱思嘉, 等. 計(jì)及風(fēng)電消納與發(fā)電集團(tuán)利益的日前協(xié)調(diào)優(yōu)化調(diào)度[J]. 電工技術(shù)學(xué)報(bào), 2021, 36(3): 579-587.

    Ma Xiufan, Wang Ge, Zhu Sijia, et al. Coordinated day-ahead optimal dispatch considering wind power consumption and the benefits of power generation group[J]. Transactions of China Electrotechnical Society, 2021, 36(3): 579-587.

    [7] 王駿, 趙潔, 劉滌塵, 等. 考慮核電參與的調(diào)峰優(yōu)化運(yùn)行模型[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2018, 38(6): 1665-1674, 1903.

    Wang Jun, Zhao Jie, Liu Dichen, et al. Optimal scheduling model of peak load regulation considering participation of nuclear power plant[J]. Proceedings of the CSEE, 2018, 38(6): 1665-1674, 1903.

    [8] 李旭東, 艾欣, 胡俊杰, 等. 計(jì)及碳交易機(jī)制的核-火-虛擬電廠三階段聯(lián)合調(diào)峰策略研究[J]. 電網(wǎng)技術(shù), 2019, 43(7): 2460-2470.

    Li Xudong, Ai Xin, Hu Junjie, et al. Three-stage combined peak regulation strategy for nuclear-thermal-virtual power plant considering carbon trading mechanism[J]. Power System Technology, 2019, 43(7): 2460-2470.

    [9] 林毅, 潘璽安, 林章歲, 等. 含高比例核電的電力系統(tǒng)調(diào)峰技術(shù)綜述[J]. 現(xiàn)代電力, 2020, 37(1): 51-58.

    Lin Yi, Pan Xian, Lin Zhangsui, et al. Review of peak regulation for power system with a high penetration of nuclear power[J]. Modern Electric Power, 2020, 37(1): 51-58.

    [10] 周剛. 考慮核電參與的多類型電源聯(lián)合調(diào)峰優(yōu)化方法[J]. 電網(wǎng)技術(shù), 2019, 43(3): 928-935.

    Zhou Gang. An approach on coordinative peaking optimization for power system with nuclear power plant participation[J]. Power System Technology, 2019, 43(3): 928-935.

    [11] Wang Jun,Zhao Jie, Ye Xiaoli, et al.Safety constraints and optimal operation of large-scale nuclear power plant participating in peak load regulation of power system[J].IET Generation,Transmission & Distribution,2017,11(13): 3332-3340.

    [12] 胡弘, 韋化, 李昭昱. 風(fēng)電接入下核電參與電力系統(tǒng)調(diào)峰的協(xié)調(diào)優(yōu)化模型[J]. 電力自動(dòng)化設(shè)備, 2020, 40(5): 31-39.

    Hu Hong, Wei Hua, Li Zhaoyu. Coordinated optimization model considering nuclear power participating in peak load regulation of power system with wind power[J]. Electric Power Automation Equipment, 2020, 40(5): 31-39.

    [13] 侯文庭, 韋化. 考慮核電可調(diào)度性的風(fēng)-光-核-水-火多源協(xié)調(diào)短期優(yōu)化調(diào)度[J]. 電工技術(shù)學(xué)報(bào), 2018, 33(12): 2873-2882.

    Hou Wenting, Wei Hua. A multi-source coordinated short-term dispatch model considering the dispatchability of nuclear power plants[J]. Transactions of China Electrotechnical Society, 2018, 33(12): 2873-2882.

    [14] Li Xue, Zhang Rufeng, Bai Linquan, et al. Stochastic low-carbon scheduling withcarbon capture power plants and coupon-based demand response[J]. Applied Energy, 2018(210):1219-1228.

    [15] 盧志剛, 隋玉珊, 馮濤, 等. 考慮儲(chǔ)熱裝置與碳捕集設(shè)備的風(fēng)電消納低碳經(jīng)濟(jì)調(diào)度[J]. 電工技術(shù)學(xué)報(bào), 2016, 31(17): 41-51.

    Lu Zhigang, Sui Yushan, Feng Tao, et al. Wind power accommodation low-carbon economic dispatch considering heat accumulator and carbon capture devices[J]. Transactions of China Electrotechnical Society, 2016, 31(17): 41-51.

    [16] 康重慶, 陳啟鑫, 夏清. 應(yīng)用于電力系統(tǒng)的碳捕集技術(shù)及其帶來的變革[J]. 電力系統(tǒng)自動(dòng)化, 2010, 34(1): 1-7.

    Kang Chongqing, Chen Qixin, Xia Qing. Innovation incurred by carbon capture technologies utilized in power systems[J]. Automation of Electric Power Systems, 2010, 34(1): 1-7.

    [17] Chen Qixin, Kang Chongqing, Xia Qing. Modeling flexible operation mechanism of CO2capture power plant and its effects on power-system operation[J]. IEEE Transactions on Energy Conversion, 2010, 25(3): 853-861.

    [18] 周任軍, 孫洪, 唐夏菲, 等. 雙碳量約束下風(fēng)電-碳捕集虛擬電廠低碳經(jīng)濟(jì)調(diào)度[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2018, 38(6): 1675-1683, 1904.

    Zhou Renjun, Sun Hong, Tang Xiafei, et al. Low-Carbon economic dispatch based on virtual power plant made up of carbon capture unit and wind power under double carbon constraint[J]. Proceedings of the CSEE, 2018, 38(6): 1675-1683, 1904.

    [19] 孫惠娟, 蒙錦輝, 彭春華. 風(fēng)-光-水-碳捕集多區(qū)域虛擬電廠協(xié)調(diào)優(yōu)化調(diào)度[J]. 電網(wǎng)技術(shù), 2019, 43(11): 4040-4051.

    Sun Huijuan, Meng Jinhui, Peng Chunhua. Coordinated optimization scheduling of multi-region virtual power plant with wind-power /photovoltaic /hydropower /carbon-capture units[J]. Power System Technology, 2019, 43(11): 4040-4051.

    [20] Fang Xin, Hu Qinran, Li Fangxing, et al. Coupon-based demand response considering wind power uncertainty: a strategic bidding model for load serving entities[J]. IEEE Transactions on Power System, 2016, 31(2): 1025-1037.

    [21] Yin Jiafu, Zhao Dongmei. Economic dispatch coordinated with information granule chance constraint goal programming under the manifold uncertainties[J]. IET Renewable Power Generation, 2019, 13(8): 1329-1337.

    [22] Cany C, Mansilla C, Mathonnière G, et al. Nuclear power supply: going against the misconceptions. Evidence of nuclear flexibility from the French experience[J]. Energy, 2018(151): 289-296.

    [23] 邵尤國(guó), 趙潔, 劉滌塵, 等. 考慮核電調(diào)峰的風(fēng)電-核電協(xié)調(diào)優(yōu)化調(diào)度[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2019, 39(4): 1018-1029.

    Shao Youguo, Zhao Jie, Liu Dichen, et al. Coordinated optimal dispatch of wind-nuclear power considering peak load shaving of nuclear power plant[J]. Proceedings of the CSEE, 2019, 39(4): 1018-1029.

    [24] 崔楊, 曾鵬, 惠鑫欣, 等. 考慮碳捕集電廠綜合靈活運(yùn)行方式的低碳經(jīng)濟(jì)調(diào)度[J]. 電網(wǎng)技術(shù), 2021, 45(5): 1877-1886.

    Cui Yang, Zeng Peng, Hui Xinxin, et al. Low-carbon economic dispatch considering the integrated flexible operation mode of carbon capture power plant[J]. Power System Technology, 2021, 45(5): 1877-1886.

    [25] 熊虎, 向鐵元, 陳紅坤, 等. 含大規(guī)模間歇式電源的模糊機(jī)會(huì)約束機(jī)組組合研究[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2013, 33(13): 36-44.

    Xiong Hu, Xiang Tieyuan, Chen Hongkun, et al. Research of fuzzy chance constrained unit commitmentcontaining large-scale intermittent power[J]. Proceedings of the CSEE, 2013, 33(13): 36-44.

    [26] IRENA. Renewable power generation costs in 2017[R]. Abu Dhabi: International Renewable Energy Agency, 2018.

    [27] 中華人民共和國(guó)國(guó)家發(fā)展和改革委員會(huì). 國(guó)家發(fā)展改革委關(guān)于調(diào)整光伏發(fā)電陸上風(fēng)電標(biāo)桿上網(wǎng)電價(jià)的通知[EB/OL]. 2016. http://www.ndrc.gov.cn/fzgggz/jggl/ zcfg/201612/t20161228_833049.html.

    [28] 崔楊, 曾鵬, 王錚, 等. 考慮碳捕集電廠能量轉(zhuǎn)移特性的棄風(fēng)消納多時(shí)間尺度調(diào)度策略[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2021, 41(3): 946-961.

    Cui Yang, Zeng Peng, Wang Zheng, et al. Multiple time scales scheduling strategy of wind power dissipation considering energy transfer characteristics of carbon capture power plant[J]. Proceedings of the CSEE, 2021, 41(3): 946-961.

    [29] 孫宇軍, 王巖, 王蓓蓓, 等. 考慮需求響應(yīng)不確定性的多時(shí)間尺度源荷互動(dòng)決策方法[J]. 電力系統(tǒng)自動(dòng)化, 2018, 42(2): 106-113, 159.

    Sun Yujun, Wang Yan, Wang Beibei, et al. Multi-time scale decision method for source-load interaction considering demand response uncertainty[J]. Automation of Electric Power Systems, 2018, 42(2): 106-113, 159.

    A Multi-Source Coordinated Optimal Scheduling Model Considering Wind-Load Uncertainty

    Zhao Dongmei Wang Haoxiang Tao Ran

    (School of Electrical and Electronic Engineering North China Electric Power University Beijing 102206 China)

    With the rapid development of nuclear power technology, it is an inevitable trend for large capacity nuclear power units to be connected to the power grid to achieve low carbon and environmental protection. However, with the increase of load peak-valley difference and the increase of intermittent power supply permeability, it is urgent for nuclear power units to share the peak pressure of the power grid in a flexible operation mode. Based on the actual operation characteristics of nuclear power units, the peak regulation depth of nuclear power units was linearized. The characteristics of peak shaving and valley filling of carbon capture power plants with integrated flexible operation mode were analyzed to realize binding peak regulation of nuclear-carbon capture units. Then, fuzzy parameters were introduced to characterize the uncertainty of wind power and load. On the premise of considering the safety and economy of nuclear power peak regulation, a multi-source coordinated fuzzy optimal scheduling model considering the benefit of wind power consumption was established by minimizing the total cost of joint operation as the objective function. The simulation example verifies the validity of the proposed model and method. The results show that the optimal dispatching model can improve the flexibility of power network dispatching and realize the economic and low-carbon operation on the basis of ensuring the safe operation of nuclear power.

    Nuclear power unit, carbon capture unit, peak load regulation, peak shaving and valley filling, optimal scheduling

    10.19595/j.cnki.1000-6753.tces.201654

    TM73

    國(guó)家重點(diǎn)研發(fā)計(jì)劃(2017YFB0902600)和國(guó)家電網(wǎng)公司科技項(xiàng)目(SGJS0000DKJS1700840)資助。

    2020-12-16

    2021-03-25

    趙冬梅 女,1965年生,教授,博士生導(dǎo)師,研究方向?yàn)殡娏ο到y(tǒng)分析與控制、新能源發(fā)電與智能電網(wǎng)。E-mail:zhao-dm@ncepu.edu.cn

    王浩翔 男,1994年生,博士研究生,研究方向?yàn)殡娏ο到y(tǒng)分析與控制。E-mail:18645750858@163.com(通信作者)

    (編輯 赫蕾)

    猜你喜歡
    存儲(chǔ)單元核電機(jī)組調(diào)峰
    新常態(tài)下電站鍋爐深度調(diào)峰改造與調(diào)試實(shí)踐
    一種28 nm工藝下抗單粒子翻轉(zhuǎn)SRAM的12T存儲(chǔ)單元設(shè)計(jì)
    我國(guó)第二臺(tái)“華龍一號(hào)”核電機(jī)組并網(wǎng)發(fā)電
    調(diào)峰保供型和普通型LNG接收站罐容計(jì)算
    煤氣與熱力(2021年5期)2021-07-22 09:02:14
    重慶市天然氣調(diào)峰儲(chǔ)氣建設(shè)的分析
    煤氣與熱力(2021年2期)2021-03-19 08:56:04
    數(shù)據(jù)在計(jì)算機(jī)內(nèi)存中的存儲(chǔ)形式及實(shí)驗(yàn)驗(yàn)證
    “華龍一號(hào)”核電機(jī)組穹頂?shù)跹b成功
    一種成本更低的全新靜態(tài)DRAM存儲(chǔ)單元
    MiR-125a-5p is Upregulated in Plasma of Residents from An Electronic Waste Recycling Site
    中國(guó)力爭(zhēng)2030年前在“一帶一路”沿線建30臺(tái)核電機(jī)組
    中亞信息(2016年3期)2016-12-01 06:08:25
    欧美 亚洲 国产 日韩一| 人妻人人澡人人爽人人| 91久久精品国产一区二区三区| 久久久久视频综合| 性少妇av在线| 亚洲精品日韩在线中文字幕| 亚洲经典国产精华液单| 精品少妇黑人巨大在线播放| 精品人妻熟女毛片av久久网站| 国产毛片在线视频| 免费看不卡的av| 男女午夜视频在线观看| 久久精品国产亚洲av天美| 亚洲精品成人av观看孕妇| av电影中文网址| 91精品伊人久久大香线蕉| 亚洲国产精品一区三区| 国产片内射在线| 精品一区在线观看国产| 可以免费在线观看a视频的电影网站 | 激情五月婷婷亚洲| 爱豆传媒免费全集在线观看| 日韩大片免费观看网站| 免费观看a级毛片全部| 交换朋友夫妻互换小说| 在线天堂最新版资源| 香蕉丝袜av| 亚洲中文av在线| 1024香蕉在线观看| 美女午夜性视频免费| 熟妇人妻不卡中文字幕| 超碰成人久久| 9色porny在线观看| √禁漫天堂资源中文www| 亚洲成国产人片在线观看| 街头女战士在线观看网站| 国产免费福利视频在线观看| 十八禁网站网址无遮挡| 久久亚洲国产成人精品v| 男女午夜视频在线观看| 久久久国产欧美日韩av| 中文字幕制服av| 韩国av在线不卡| 侵犯人妻中文字幕一二三四区| 国产成人aa在线观看| 亚洲国产色片| 亚洲精品日本国产第一区| 咕卡用的链子| 国产1区2区3区精品| 嫩草影院入口| 毛片一级片免费看久久久久| 毛片一级片免费看久久久久| 久久韩国三级中文字幕| 9色porny在线观看| 精品第一国产精品| 国产激情久久老熟女| xxx大片免费视频| 亚洲国产欧美在线一区| 黄片无遮挡物在线观看| 男女边摸边吃奶| 天堂中文最新版在线下载| 性高湖久久久久久久久免费观看| 亚洲视频免费观看视频| 丝袜喷水一区| 亚洲国产成人一精品久久久| 午夜精品国产一区二区电影| 曰老女人黄片| 一级黄片播放器| 在线观看www视频免费| 九九爱精品视频在线观看| 极品少妇高潮喷水抽搐| 国产精品香港三级国产av潘金莲 | 亚洲精品自拍成人| 国产精品国产三级专区第一集| 国产精品成人在线| 少妇被粗大猛烈的视频| 亚洲综合色惰| 人人妻人人澡人人看| 精品久久蜜臀av无| 一级毛片电影观看| 天堂8中文在线网| 久久国产亚洲av麻豆专区| 哪个播放器可以免费观看大片| 国产精品成人在线| 国产精品亚洲av一区麻豆 | 欧美日韩一级在线毛片| 久久精品久久久久久噜噜老黄| 亚洲精品国产一区二区精华液| 欧美中文综合在线视频| 亚洲欧洲精品一区二区精品久久久 | 自线自在国产av| 免费黄频网站在线观看国产| 成年女人在线观看亚洲视频| 精品人妻在线不人妻| 国产又爽黄色视频| 国产色婷婷99| 亚洲人成网站在线观看播放| 另类亚洲欧美激情| 亚洲视频免费观看视频| 深夜精品福利| 欧美最新免费一区二区三区| 在线天堂最新版资源| 久久久亚洲精品成人影院| 欧美亚洲 丝袜 人妻 在线| 青草久久国产| 日本午夜av视频| 久久久精品国产亚洲av高清涩受| 中文字幕色久视频| 成年女人毛片免费观看观看9 | 街头女战士在线观看网站| 欧美人与性动交α欧美精品济南到 | 国产精品偷伦视频观看了| 美女主播在线视频| 又大又黄又爽视频免费| 久久久久视频综合| 熟女少妇亚洲综合色aaa.| www.熟女人妻精品国产| 国产激情久久老熟女| 日日爽夜夜爽网站| 午夜福利视频精品| 热re99久久国产66热| 日日摸夜夜添夜夜爱| 国产免费视频播放在线视频| 久久 成人 亚洲| 黄色毛片三级朝国网站| 亚洲国产av新网站| 久久久久精品久久久久真实原创| 欧美人与性动交α欧美精品济南到 | 人成视频在线观看免费观看| 黄频高清免费视频| 老鸭窝网址在线观看| 成人毛片60女人毛片免费| 男的添女的下面高潮视频| 只有这里有精品99| 69精品国产乱码久久久| 国产 一区精品| 欧美日韩精品成人综合77777| 亚洲伊人久久精品综合| 少妇人妻 视频| 最黄视频免费看| 久久精品久久久久久噜噜老黄| 中文字幕人妻熟女乱码| 精品国产乱码久久久久久小说| av福利片在线| 亚洲色图 男人天堂 中文字幕| 日产精品乱码卡一卡2卡三| 国产黄色免费在线视频| 91午夜精品亚洲一区二区三区| 老鸭窝网址在线观看| 日韩中字成人| 女人精品久久久久毛片| 亚洲内射少妇av| 制服人妻中文乱码| 久久精品久久精品一区二区三区| 伦理电影免费视频| 中国国产av一级| 亚洲国产精品一区三区| 人人澡人人妻人| 校园人妻丝袜中文字幕| av卡一久久| 亚洲综合精品二区| 欧美日韩国产mv在线观看视频| 亚洲国产色片| 美女国产高潮福利片在线看| 欧美精品亚洲一区二区| 欧美亚洲日本最大视频资源| 秋霞伦理黄片| 人妻 亚洲 视频| 午夜免费观看性视频| 亚洲精华国产精华液的使用体验| 国产精品二区激情视频| 免费观看在线日韩| 午夜福利网站1000一区二区三区| 日韩人妻精品一区2区三区| 少妇被粗大猛烈的视频| 亚洲国产欧美在线一区| 国产精品欧美亚洲77777| av国产久精品久网站免费入址| 国产不卡av网站在线观看| 波野结衣二区三区在线| 日本-黄色视频高清免费观看| 日韩人妻精品一区2区三区| 看非洲黑人一级黄片| av免费在线看不卡| 国产精品女同一区二区软件| 久久国产亚洲av麻豆专区| 欧美人与性动交α欧美精品济南到 | 秋霞在线观看毛片| 亚洲精品久久午夜乱码| 人人澡人人妻人| 99久久综合免费| 久久久久国产精品人妻一区二区| 成年av动漫网址| 亚洲欧美成人精品一区二区| 99国产综合亚洲精品| 国产 一区精品| 在线观看三级黄色| xxxhd国产人妻xxx| 亚洲国产色片| 黑人巨大精品欧美一区二区蜜桃| 少妇的逼水好多| 制服丝袜香蕉在线| 街头女战士在线观看网站| 欧美精品国产亚洲| av一本久久久久| 国产精品一区二区在线不卡| a 毛片基地| 国产色婷婷99| 亚洲av日韩在线播放| 亚洲色图 男人天堂 中文字幕| 国产亚洲午夜精品一区二区久久| 亚洲国产最新在线播放| 妹子高潮喷水视频| 午夜福利乱码中文字幕| 国产免费福利视频在线观看| 国产97色在线日韩免费| 咕卡用的链子| 熟妇人妻不卡中文字幕| 国产毛片在线视频| 免费观看av网站的网址| 免费播放大片免费观看视频在线观看| 国产人伦9x9x在线观看 | 捣出白浆h1v1| 欧美av亚洲av综合av国产av | 考比视频在线观看| 免费观看在线日韩| 免费不卡的大黄色大毛片视频在线观看| 色视频在线一区二区三区| 亚洲精品中文字幕在线视频| 日韩精品免费视频一区二区三区| 成人免费观看视频高清| 精品福利永久在线观看| 国产精品三级大全| 又黄又粗又硬又大视频| 丰满饥渴人妻一区二区三| 亚洲男人天堂网一区| 免费观看av网站的网址| 99久久中文字幕三级久久日本| 国产av国产精品国产| 一级黄片播放器| 亚洲欧洲精品一区二区精品久久久 | 国产白丝娇喘喷水9色精品| 只有这里有精品99| 国产精品免费大片| 国产精品久久久久久精品电影小说| 18禁国产床啪视频网站| 欧美精品人与动牲交sv欧美| 黄片播放在线免费| 纯流量卡能插随身wifi吗| 日韩伦理黄色片| 水蜜桃什么品种好| 纯流量卡能插随身wifi吗| 国产在线免费精品| 日韩一卡2卡3卡4卡2021年| 久久午夜福利片| 纯流量卡能插随身wifi吗| 国产免费现黄频在线看| 免费人妻精品一区二区三区视频| 久久青草综合色| av电影中文网址| 亚洲精品成人av观看孕妇| 成年人免费黄色播放视频| 日韩制服骚丝袜av| 9191精品国产免费久久| h视频一区二区三区| 最新的欧美精品一区二区| 久久99热这里只频精品6学生| 男人舔女人的私密视频| 大话2 男鬼变身卡| 激情五月婷婷亚洲| 午夜激情av网站| 熟女av电影| 亚洲国产av新网站| 色婷婷久久久亚洲欧美| 国产一区有黄有色的免费视频| 免费高清在线观看日韩| 国产亚洲午夜精品一区二区久久| 少妇熟女欧美另类| 国产av一区二区精品久久| 91午夜精品亚洲一区二区三区| 97在线视频观看| 国产精品熟女久久久久浪| av在线app专区| 成人亚洲精品一区在线观看| 精品一区二区免费观看| 性色avwww在线观看| 乱人伦中国视频| av在线观看视频网站免费| 在线 av 中文字幕| 国产免费福利视频在线观看| 最近手机中文字幕大全| 久久久久国产网址| 人人澡人人妻人| 精品人妻在线不人妻| 国产精品秋霞免费鲁丝片| 999精品在线视频| 母亲3免费完整高清在线观看 | 久久精品国产亚洲av天美| 只有这里有精品99| 国产日韩欧美视频二区| 视频在线观看一区二区三区| 九色亚洲精品在线播放| 电影成人av| 高清欧美精品videossex| 国产白丝娇喘喷水9色精品| 欧美日韩视频高清一区二区三区二| 日韩精品免费视频一区二区三区| 久久久国产一区二区| 十八禁高潮呻吟视频| 看非洲黑人一级黄片| 久久久精品区二区三区| 91国产中文字幕| 777久久人妻少妇嫩草av网站| 人体艺术视频欧美日本| 性色avwww在线观看| videosex国产| 天堂8中文在线网| 热99久久久久精品小说推荐| 多毛熟女@视频| 欧美97在线视频| 国产精品秋霞免费鲁丝片| 黄色 视频免费看| 免费高清在线观看视频在线观看| 99热全是精品| 久久99精品国语久久久| 性色av一级| 亚洲av国产av综合av卡| 欧美成人午夜精品| 免费高清在线观看视频在线观看| 国产精品国产三级国产专区5o| 电影成人av| 国产片特级美女逼逼视频| 99精国产麻豆久久婷婷| 一本久久精品| 天堂中文最新版在线下载| 午夜老司机福利剧场| 亚洲经典国产精华液单| 9191精品国产免费久久| 久久精品人人爽人人爽视色| 一二三四中文在线观看免费高清| 国产精品.久久久| 亚洲色图综合在线观看| 捣出白浆h1v1| 免费观看无遮挡的男女| 性高湖久久久久久久久免费观看| av卡一久久| 国产欧美亚洲国产| 久久人人爽av亚洲精品天堂| 欧美激情 高清一区二区三区| 激情五月婷婷亚洲| 亚洲精品一二三| 男女高潮啪啪啪动态图| 18禁国产床啪视频网站| 人人澡人人妻人| 国产伦理片在线播放av一区| 国产一区二区 视频在线| 国产又色又爽无遮挡免| 国产一区二区 视频在线| 国产又色又爽无遮挡免| 激情五月婷婷亚洲| 黄片小视频在线播放| 成人18禁高潮啪啪吃奶动态图| av视频免费观看在线观看| 亚洲欧洲日产国产| 成人午夜精彩视频在线观看| 国产亚洲最大av| 午夜免费观看性视频| 日韩三级伦理在线观看| 超碰97精品在线观看| 国产有黄有色有爽视频| 精品久久久久久电影网| 亚洲欧洲精品一区二区精品久久久 | 成人黄色视频免费在线看| 高清不卡的av网站| 亚洲图色成人| 亚洲av欧美aⅴ国产| 最近中文字幕高清免费大全6| 亚洲av在线观看美女高潮| 中文字幕人妻丝袜制服| 人人澡人人妻人| 秋霞伦理黄片| 91久久精品国产一区二区三区| www日本在线高清视频| 久久久国产精品麻豆| 亚洲成人手机| 肉色欧美久久久久久久蜜桃| 久久久久精品性色| 1024视频免费在线观看| 日韩中文字幕欧美一区二区 | 亚洲第一区二区三区不卡| 99久久精品国产国产毛片| 久久人妻熟女aⅴ| 热99国产精品久久久久久7| 制服丝袜香蕉在线| 国产精品99久久99久久久不卡 | 亚洲精品国产色婷婷电影| 80岁老熟妇乱子伦牲交| 黑人欧美特级aaaaaa片| 成年人免费黄色播放视频| 啦啦啦视频在线资源免费观看| 日本黄色日本黄色录像| 亚洲国产欧美网| 国产男女超爽视频在线观看| 日本欧美国产在线视频| 91午夜精品亚洲一区二区三区| 久久人人爽人人片av| av在线观看视频网站免费| 另类精品久久| 亚洲欧美成人综合另类久久久| 丝瓜视频免费看黄片| 日日撸夜夜添| 青春草视频在线免费观看| 男女边吃奶边做爰视频| 熟女av电影| 午夜免费鲁丝| 亚洲成色77777| 国语对白做爰xxxⅹ性视频网站| 免费观看在线日韩| 日韩 亚洲 欧美在线| 亚洲人成电影观看| 久久久久久人人人人人| 99久久中文字幕三级久久日本| 一级爰片在线观看| 亚洲欧洲日产国产| kizo精华| 啦啦啦在线免费观看视频4| 亚洲成色77777| 午夜福利影视在线免费观看| av网站免费在线观看视频| 侵犯人妻中文字幕一二三四区| 免费人妻精品一区二区三区视频| 三上悠亚av全集在线观看| 国产麻豆69| 久久久国产欧美日韩av| 夫妻午夜视频| 十八禁网站网址无遮挡| 精品国产国语对白av| 久久人妻熟女aⅴ| 下体分泌物呈黄色| 青草久久国产| 韩国精品一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 午夜免费观看性视频| 国产一区二区 视频在线| 黄色配什么色好看| 精品久久久精品久久久| 一级片免费观看大全| 日韩中文字幕欧美一区二区 | 亚洲av中文av极速乱| 久久久精品区二区三区| 免费在线观看黄色视频的| 涩涩av久久男人的天堂| 男女边吃奶边做爰视频| 捣出白浆h1v1| 日韩免费高清中文字幕av| 中文天堂在线官网| 亚洲精品在线美女| 少妇人妻 视频| tube8黄色片| 国产精品三级大全| 日本黄色日本黄色录像| 亚洲欧美清纯卡通| 一级,二级,三级黄色视频| 国产精品久久久久久精品古装| 国产国语露脸激情在线看| 精品人妻一区二区三区麻豆| 夫妻午夜视频| 不卡视频在线观看欧美| 精品久久久精品久久久| 1024视频免费在线观看| 久久久久久久国产电影| 中文精品一卡2卡3卡4更新| 日韩在线高清观看一区二区三区| 亚洲国产av影院在线观看| 久久人人爽人人片av| 欧美日韩精品成人综合77777| videos熟女内射| 欧美中文综合在线视频| 日本黄色日本黄色录像| 一本大道久久a久久精品| 一级片免费观看大全| 婷婷成人精品国产| av免费在线看不卡| 久久毛片免费看一区二区三区| 精品少妇久久久久久888优播| 黑丝袜美女国产一区| a 毛片基地| 波多野结衣一区麻豆| 涩涩av久久男人的天堂| 熟女少妇亚洲综合色aaa.| 国产精品一二三区在线看| 免费人妻精品一区二区三区视频| 激情视频va一区二区三区| 哪个播放器可以免费观看大片| 久热久热在线精品观看| 亚洲熟女精品中文字幕| 国产精品久久久久久精品古装| 亚洲熟女精品中文字幕| 午夜精品国产一区二区电影| 免费在线观看黄色视频的| 亚洲精品成人av观看孕妇| 美女大奶头黄色视频| 在现免费观看毛片| av.在线天堂| 一本大道久久a久久精品| 国产有黄有色有爽视频| 人妻人人澡人人爽人人| 波野结衣二区三区在线| 日日摸夜夜添夜夜爱| 肉色欧美久久久久久久蜜桃| 麻豆精品久久久久久蜜桃| 又大又黄又爽视频免费| 精品一品国产午夜福利视频| 啦啦啦中文免费视频观看日本| 欧美日韩一区二区视频在线观看视频在线| 97精品久久久久久久久久精品| 最新中文字幕久久久久| h视频一区二区三区| 亚洲色图综合在线观看| 最近最新中文字幕大全免费视频 | 久久久久网色| 亚洲美女黄色视频免费看| 91成人精品电影| 91久久精品国产一区二区三区| 久久久久久久久久久免费av| 国产精品久久久久久精品电影小说| 免费观看a级毛片全部| 久久久国产欧美日韩av| 久久综合国产亚洲精品| 中文字幕人妻丝袜制服| 亚洲成人一二三区av| 日日撸夜夜添| 久久av网站| 丰满少妇做爰视频| 亚洲精品日韩在线中文字幕| 国产在线一区二区三区精| 免费久久久久久久精品成人欧美视频| 九草在线视频观看| 在线 av 中文字幕| 欧美日韩一级在线毛片| 日韩视频在线欧美| 日本av手机在线免费观看| 狠狠婷婷综合久久久久久88av| 韩国高清视频一区二区三区| 亚洲精品乱久久久久久| 久久精品国产亚洲av涩爱| 免费日韩欧美在线观看| 男女啪啪激烈高潮av片| 水蜜桃什么品种好| 亚洲精品日韩在线中文字幕| 国产成人aa在线观看| 永久网站在线| 99久久精品国产国产毛片| 黄频高清免费视频| 热99国产精品久久久久久7| 欧美精品一区二区大全| 狂野欧美激情性bbbbbb| 欧美亚洲 丝袜 人妻 在线| 婷婷色av中文字幕| 日本色播在线视频| 亚洲国产成人一精品久久久| 亚洲国产av新网站| 亚洲精品久久成人aⅴ小说| 在线观看免费高清a一片| 亚洲激情五月婷婷啪啪| 极品少妇高潮喷水抽搐| 国产精品女同一区二区软件| 日本-黄色视频高清免费观看| 国产成人欧美| 熟女少妇亚洲综合色aaa.| 女人久久www免费人成看片| 亚洲欧美精品综合一区二区三区 | 99热全是精品| 女人精品久久久久毛片| 亚洲成av片中文字幕在线观看 | 夫妻性生交免费视频一级片| 人妻系列 视频| 国产精品亚洲av一区麻豆 | 最黄视频免费看| 国产精品香港三级国产av潘金莲 | 在线观看美女被高潮喷水网站| 中文字幕人妻丝袜制服| 女人被躁到高潮嗷嗷叫费观| 一级毛片电影观看| 国产熟女欧美一区二区| 国产成人欧美| 国产精品 欧美亚洲| 99国产综合亚洲精品| 久久久久久免费高清国产稀缺| 国产精品人妻久久久影院| 99国产精品免费福利视频| 国产片特级美女逼逼视频| 久久久精品国产亚洲av高清涩受| 欧美激情 高清一区二区三区| 99久国产av精品国产电影| 最黄视频免费看| 香蕉国产在线看| 一个人免费看片子| 亚洲色图综合在线观看| 日韩一本色道免费dvd| 亚洲伊人色综图| 我的亚洲天堂| 亚洲精品aⅴ在线观看| 丝瓜视频免费看黄片| 久久热在线av| 天天操日日干夜夜撸| 亚洲图色成人| 观看美女的网站| 亚洲精品aⅴ在线观看| 亚洲国产精品国产精品| 9色porny在线观看| 国产福利在线免费观看视频| 国产精品一二三区在线看| 欧美日韩亚洲国产一区二区在线观看 | 毛片一级片免费看久久久久|