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

    省間電力現(xiàn)貨交易優(yōu)化設(shè)計(jì)與定價(jià)機(jī)制

    2024-03-06 08:21:12孫宏斌
    電力系統(tǒng)自動(dòng)化 2024年4期
    關(guān)鍵詞:區(qū)域

    劉 昊,郭 燁,孫宏斌,3,4

    (1.清華-伯克利深圳學(xué)院,清華大學(xué)深圳國(guó)際研究生院,廣東省深圳市 518000;2.電力規(guī)劃設(shè)計(jì)總院,北京市 100120;3.電力系統(tǒng)及大型發(fā)電設(shè)備安全控制和仿真國(guó)家重點(diǎn)實(shí)驗(yàn)室(清華大學(xué)),北京市 100084;4.清華大學(xué)電機(jī)工程與應(yīng)用電子技術(shù)系,北京市 100084)

    0 引言

    2021 年,中國(guó)提出“碳達(dá)峰·碳中和”戰(zhàn)略目標(biāo),不斷提高新能源發(fā)電占比,加快建設(shè)新型電力系統(tǒng)[1-3]??紤]新能源發(fā)電與負(fù)荷的逆向地理分布現(xiàn)狀,完善全國(guó)統(tǒng)一電力市場(chǎng)體系有利于電力資源在中國(guó)更大范圍內(nèi)實(shí)現(xiàn)優(yōu)化配置。而省間電力現(xiàn)貨市場(chǎng)作為全國(guó)統(tǒng)一電力市場(chǎng)體系的關(guān)鍵環(huán)節(jié),可以促進(jìn)中國(guó)西部地區(qū)新能源發(fā)電的消納,加強(qiáng)不同區(qū)域電網(wǎng)間的互聯(lián)互濟(jì),加快“碳達(dá)峰·碳中和”目標(biāo)的實(shí)現(xiàn)進(jìn)程[4-5]。

    對(duì)于跨區(qū)域電力現(xiàn)貨交易,國(guó)外已有諸多研究和應(yīng)用。在市場(chǎng)出清方面,美國(guó)中西部獨(dú)立系統(tǒng)運(yùn)營(yíng) 商(midcontinent independent system operator,MISO)和PJM 獨(dú)立系統(tǒng)運(yùn)營(yíng)商(Pennsylvania-New Jersey-Maryland independent system operator,PJMISO)、MISO 和 美 國(guó) 西 南 電 網(wǎng)(southwest power pool,SPP)、紐約獨(dú)立系統(tǒng)運(yùn)營(yíng)商(New York independent system operator,NYISO)和PJMISO 分別就區(qū)域市場(chǎng)間協(xié)調(diào)(market-to-market,M2M)達(dá)成了聯(lián)合運(yùn)營(yíng)協(xié)議,以管理實(shí)時(shí)運(yùn)營(yíng)中的輸電線路阻 塞[6]。 協(xié) 調(diào) 交 易 調(diào) 度(coordinated transaction scheduling,CTS)已經(jīng)在不同獨(dú)立系統(tǒng)運(yùn)營(yíng)商(independent system operator,ISO)間得到應(yīng)用,但其依賴于代理節(jié)點(diǎn)的設(shè)置,會(huì)產(chǎn)生模型誤差,造成環(huán)流效應(yīng)[7-9]。廣義協(xié)調(diào)交易調(diào)度(generalized CTS,GCTS)進(jìn)一步消除了對(duì)代理節(jié)點(diǎn)的依賴,在一定條件下可以與聯(lián)合經(jīng)濟(jì)調(diào)度(joint economic dispatch,JED)的 結(jié) 果 等 價(jià)[10]。在 歐 洲,可 用 傳 輸 容 量(available transfer capacity,ATC)模型僅考慮區(qū)域間交易量限制并分配聯(lián)絡(luò)線容量,出清結(jié)果較為保守[11]。基于流的模型采用功率轉(zhuǎn)移分布因子量化每個(gè)區(qū)域?qū)β?lián)絡(luò)線的影響,相較于ATC 模型可以提高區(qū)域間聯(lián)絡(luò)線的利用率,降低系統(tǒng)發(fā)電成本[12]。

    在市場(chǎng)結(jié)算方面,美國(guó)的CTS 機(jī)制中,跨區(qū)域投標(biāo)者根據(jù)ISO 預(yù)測(cè)的代理節(jié)點(diǎn)的節(jié)點(diǎn)邊際電價(jià)(locational marginal price,LMP)提交投標(biāo)向量。在實(shí)時(shí)市場(chǎng)中,ISO 通過(guò)增加交易頻率和降低交易費(fèi)用來(lái)激勵(lì)投標(biāo)人參與區(qū)域間電力交易,并采用代理節(jié)點(diǎn)的LMP 來(lái)進(jìn)行結(jié)算[7]。在歐洲,不同國(guó)家和地區(qū)根據(jù)線路擁堵情況靈活劃分報(bào)價(jià)區(qū),每個(gè)報(bào)價(jià)區(qū)采用統(tǒng)一價(jià)格[13-14]。歐洲各國(guó)聯(lián)合開(kāi)展了區(qū)域價(jià)格耦合(price coupling of regions,PCR)項(xiàng)目,開(kāi)發(fā)了單一價(jià)格耦合算法EUPHEMIA(pan-European hybrid electricity market integration algorithm),用于求解ATC 和Flow-based 模型并獲得隱式拍賣價(jià)格[15]。

    在中國(guó),北京電力交易中心和廣州電力交易中心已經(jīng)成立,用于組織中國(guó)跨區(qū)跨省電力交易的出清和結(jié)算[16-17]。目前,中國(guó)國(guó)內(nèi)的省間現(xiàn)貨市場(chǎng)已具有較大規(guī)模。以國(guó)家電網(wǎng)有限公司為例,2022年,超過(guò)6 000 家發(fā)電企業(yè)參與省間現(xiàn)貨交易。全年累計(jì)交易電量27.8 TW·h,最大單日交易電量超過(guò)19 GW·h,發(fā)電側(cè)平均交易電價(jià)0.87 元/(kW·h)[16,18]。在省間現(xiàn)貨交易的出清過(guò)程中,每個(gè)省份采用一個(gè)或多個(gè)代理節(jié)點(diǎn)供市場(chǎng)參與者競(jìng)價(jià)。每筆交易的買賣雙方都需要在出價(jià)時(shí)指定交易路徑,調(diào)度中心基于交易路徑來(lái)出清,根據(jù)報(bào)價(jià)高低依次扣除對(duì)應(yīng)路徑的相應(yīng)容量,直至節(jié)點(diǎn)間交易路徑無(wú)可用容量[19]。這種基于交易路徑的市場(chǎng)出清方法雖然易于理解,但會(huì)導(dǎo)致交易流和物理潮流之間存在差異,在線路容量不足時(shí),會(huì)造成實(shí)際物理潮流越限。

    本文針對(duì)省間電力現(xiàn)貨交易設(shè)計(jì)與定價(jià)機(jī)制開(kāi)展相關(guān)研究,給出了不依賴于跨區(qū)域投標(biāo)者的省間電力現(xiàn)貨交易流程和市場(chǎng)出清模型,采用分區(qū)邊際定價(jià)的方法為每個(gè)報(bào)價(jià)區(qū)計(jì)算出清價(jià)格,并計(jì)算區(qū)域間各交易所需承擔(dān)的輸電費(fèi)用。

    1 省間電力現(xiàn)貨交易流程設(shè)計(jì)

    本文基于分區(qū)邊際定價(jià)原則設(shè)計(jì)的不依賴于跨區(qū)域投標(biāo)者的省間電力現(xiàn)貨交易流程如下:

    1)按地理位置或行政因素將全局系統(tǒng)劃分為多個(gè)報(bào)價(jià)區(qū);

    2)基于多區(qū)域電力系統(tǒng)中的網(wǎng)絡(luò)等值方法,獲得用于省間現(xiàn)貨交易出清的簡(jiǎn)化公共網(wǎng)絡(luò);

    3)各區(qū)域運(yùn)營(yíng)商(如省級(jí)調(diào)度中心和電力交易中心)收集本區(qū)域投標(biāo),協(xié)同完成市場(chǎng)出清,計(jì)算得到各筆交易的出清數(shù)量;

    4)各區(qū)域運(yùn)營(yíng)商將該區(qū)域的邊際報(bào)價(jià)作為交易出清價(jià)格對(duì)本區(qū)域交易進(jìn)行結(jié)算,進(jìn)一步計(jì)算各筆交易的出清費(fèi)用和輸電費(fèi)用。

    本章將進(jìn)一步討論報(bào)價(jià)區(qū)劃分和得到簡(jiǎn)化公共網(wǎng)絡(luò)的方法,市場(chǎng)出清模型和結(jié)算方法將分別在第2 章和第3 章中進(jìn)行闡述。

    1.1 報(bào)價(jià)區(qū)劃分

    在歐洲,報(bào)價(jià)區(qū)通常按地理位置劃分,在開(kāi)展跨區(qū)域交易的過(guò)程中,各國(guó)家內(nèi)部的輸電線路一般視為容量充足。主要考慮聯(lián)絡(luò)線容量約束,這極大提高了出清效率,適用于線路容量較為充足的多區(qū)域電力系統(tǒng)[20]。類似地,可以將中國(guó)的每個(gè)省作為一個(gè)報(bào)價(jià)區(qū),在組織省間電力現(xiàn)貨交易的過(guò)程中,各報(bào)價(jià)區(qū)內(nèi)的傳輸線路一般視為非阻塞。值得注意的是,本文將區(qū)域內(nèi)易阻塞線路納入“關(guān)鍵設(shè)備”集合中進(jìn)行考慮,并將其容量約束納入市場(chǎng)出清模型。

    以圖1 所示的兩區(qū)域電力系統(tǒng)為例進(jìn)行說(shuō)明,假定區(qū)域i從區(qū)域j進(jìn)口電力。根據(jù)前文描述,將區(qū)域i和區(qū)域j均作為報(bào)價(jià)區(qū),每個(gè)報(bào)價(jià)區(qū)由省級(jí)調(diào)度中心和電力交易中心負(fù)責(zé)運(yùn)營(yíng)。在電力系統(tǒng)實(shí)時(shí)運(yùn)行時(shí),存在因各種不確定性因素導(dǎo)致省內(nèi)輸電線路的阻塞問(wèn)題。為解決此問(wèn)題,一方面,在省間現(xiàn)貨市場(chǎng)完成出清確定各省邊界狀態(tài)后,省內(nèi)市場(chǎng)出清過(guò)程中會(huì)計(jì)及省內(nèi)線路容量約束,進(jìn)一步降低線路潮流越限風(fēng)險(xiǎn)[18];另一方面,針對(duì)實(shí)時(shí)運(yùn)行中出現(xiàn)的阻塞或潮流越限問(wèn)題,采用機(jī)組再調(diào)度、改變線路拓?fù)?、調(diào)整靈活負(fù)荷水平、儲(chǔ)能出力等技術(shù)手段進(jìn)行處理,具體可參考文獻(xiàn)[6,21]。

    圖1 兩區(qū)域電力系統(tǒng)示例Fig.1 Example of two-region power system

    1.2 簡(jiǎn)化公共網(wǎng)絡(luò)形成

    在省間現(xiàn)貨交易過(guò)程中,考慮多區(qū)域電力系統(tǒng)網(wǎng)絡(luò)結(jié)構(gòu)一般較為復(fù)雜,如果對(duì)全系統(tǒng)進(jìn)行完整建模會(huì)降低計(jì)算效率,并增加通信負(fù)擔(dān)。本文采用網(wǎng)絡(luò)等值方法建立簡(jiǎn)化公共網(wǎng)絡(luò),在保證出清結(jié)果最優(yōu)的前提下提高計(jì)算效率。

    簡(jiǎn)化公共網(wǎng)絡(luò)由多區(qū)域電力系統(tǒng)中的“關(guān)鍵設(shè)備”組成,通常包含區(qū)域間的聯(lián)絡(luò)線和受區(qū)域間交易影響較大的每個(gè)區(qū)域的擁堵線路、變壓器和換流站等。值得注意的是,區(qū)域間電力交易對(duì)設(shè)備的影響通常由功率轉(zhuǎn)移分布因子進(jìn)行量度,其可以基于直流最優(yōu)潮流模型計(jì)算得到[22]。

    區(qū)域間電力交易與電力設(shè)備(如輸電線路等)的功率轉(zhuǎn)移分布因子高于區(qū)域運(yùn)營(yíng)商設(shè)定的閾值時(shí),則將相應(yīng)的設(shè)備納入“關(guān)鍵設(shè)備”集合,置入簡(jiǎn)化公共網(wǎng)絡(luò)。

    以圖1 所示的兩區(qū)域電力系統(tǒng)為例,區(qū)域i和區(qū)域j分別為省i和省j,此時(shí)“關(guān)鍵設(shè)備”包括兩條省間聯(lián)絡(luò)線1-3 和2-4。在將區(qū)域i和區(qū)域j分別作為報(bào)價(jià)區(qū)后,對(duì)每個(gè)省份網(wǎng)絡(luò)對(duì)外作網(wǎng)絡(luò)等值,得到兩區(qū)域電力系統(tǒng)的公共網(wǎng)絡(luò)如圖2 所示。圖中:虛線表示區(qū)域i和區(qū)域j的對(duì)外等價(jià)支路。此外,區(qū)域i的等效邊界功率注入如式(1)所示。

    圖2 兩區(qū)域電力系統(tǒng)簡(jiǎn)化公共網(wǎng)絡(luò)Fig.2 Simplified common network for two-region power system

    式中:qi為區(qū)域i的跨區(qū)域交易;ei為關(guān)聯(lián)矩陣,功率注入節(jié)點(diǎn)對(duì)應(yīng)元素為1,其他節(jié)點(diǎn)對(duì)應(yīng)元素為0;BBi-Ii為與區(qū)域i中的邊界節(jié)點(diǎn)和內(nèi)部節(jié)點(diǎn)相關(guān)的節(jié)點(diǎn)導(dǎo)納矩陣;BIi-Ii為區(qū)域i的自導(dǎo)納矩陣。BBi-Ii和B-1Ii-Ii包含的元素均為支路電抗的倒數(shù)。區(qū)域j的等效邊界功率注入也可以采用類似方法計(jì)算得到,表示為P?j。

    直流潮流模型在高壓輸電網(wǎng)中已有廣泛應(yīng)用,準(zhǔn)確度較高[23],本文基于直流潮流模型對(duì)潮流轉(zhuǎn)移因子進(jìn)行求解。此外,采用目前在多區(qū)域電力系統(tǒng)中已經(jīng)應(yīng)用的靜態(tài)網(wǎng)絡(luò)等值方法得到簡(jiǎn)化公共網(wǎng)絡(luò),該方法在穩(wěn)態(tài)情況下不引入額外的建模誤差,詳見(jiàn)文獻(xiàn)[23-24]。

    2 市場(chǎng)出清模型

    不依賴于跨區(qū)域投標(biāo)者直接參與的省間現(xiàn)貨機(jī)制的市場(chǎng)出清流程包括以下步驟:

    步驟1:與傳統(tǒng)發(fā)電負(fù)荷參與現(xiàn)貨市場(chǎng)的投標(biāo)方式相同,參與省間現(xiàn)貨交易的市場(chǎng)主體對(duì)未來(lái)某一時(shí)段進(jìn)行報(bào)量和報(bào)價(jià);

    步驟2:省級(jí)或區(qū)域電力交易中心收集投標(biāo)向量以形成該區(qū)域的供需曲線;

    步驟3:區(qū)域間協(xié)調(diào)者收集所有報(bào)價(jià)區(qū)的供需曲線并進(jìn)行聯(lián)合市場(chǎng)出清。

    下文將對(duì)以上3 個(gè)步驟分別進(jìn)行說(shuō)明。步驟1中,市場(chǎng)參與者k向省級(jí)電力交易中心提交分時(shí)電力-價(jià)格曲線。在投標(biāo)價(jià)格方面,賣方和買方在各交易出清時(shí)段申報(bào)多段投標(biāo)價(jià)格;在報(bào)量方面,要求各段報(bào)量出力區(qū)間首尾相接,各出力報(bào)價(jià)曲線隨出力增加應(yīng)滿足單調(diào)非遞減。在具體描述方式上,市場(chǎng)參與者k的電力-價(jià)格曲線的第s個(gè)區(qū)間可采用投標(biāo)向量Rk,s表示,該向量格式為:

    式中:tk、nk分別為投標(biāo)的目標(biāo)時(shí)間段、購(gòu)/售電的節(jié)點(diǎn)位置編號(hào);pk,s、-qk,s、qˉk,s分別為第s個(gè)區(qū)間的投標(biāo)價(jià)格、投標(biāo)數(shù)量下限、投標(biāo)數(shù)量上限。考慮信息披露的要求[19],在投標(biāo)環(huán)節(jié)各市場(chǎng)參與者無(wú)法掌握網(wǎng)架信息。市場(chǎng)出清過(guò)程中充分考慮目標(biāo)電網(wǎng)線路容量約束,各參與者投標(biāo)環(huán)節(jié)不考慮目標(biāo)電網(wǎng)架構(gòu),認(rèn)為其不影響市場(chǎng)活力。

    步驟2 中,省級(jí)或區(qū)域電力交易中心收集所有投標(biāo)向量,形成各個(gè)區(qū)域的供給或需求曲線,如圖3所示。此時(shí),省級(jí)電力交易中心代表其區(qū)域內(nèi)的市場(chǎng)參與者在省間現(xiàn)貨市場(chǎng)中進(jìn)行電力買賣。根據(jù)投標(biāo)位置,區(qū)域運(yùn)營(yíng)商可以通過(guò)式(1)計(jì)算區(qū)域i中所有投標(biāo)對(duì)邊界總線的等效功率注入。

    圖3 各區(qū)域供給和需求曲線Fig.3 Supply and demand curve of each region

    步驟3 中,考慮多區(qū)域電力市場(chǎng)中不同區(qū)域間的跨區(qū)域交易,在各區(qū)域匯總形成各自的電力-價(jià)格曲線后,建立多區(qū)域市場(chǎng)的省間現(xiàn)貨交易出清模型,如式(3)—式(9)所示。

    s.t.

    式中:R為購(gòu)電區(qū)域集合;C為售電區(qū)域集合;eR、eC分別為購(gòu)電區(qū)域、售電區(qū)域的關(guān)聯(lián)矩陣;決策變量qR、qC分別為購(gòu)電區(qū)域、售電區(qū)域的投標(biāo)出清數(shù)量向量;pR、pC分別為購(gòu)電區(qū)域、售電區(qū)域的報(bào)價(jià)向量;St,B為聯(lián)絡(luò)線和邊界節(jié)點(diǎn)之間的功率轉(zhuǎn)移分布因子矩陣;SR,B、SC,B分別為購(gòu)電區(qū)域、售電區(qū)域的“關(guān)鍵設(shè)備”和邊界節(jié)點(diǎn)之間的功率轉(zhuǎn)移分布因子矩陣;SR,R為購(gòu)電區(qū)域“關(guān)鍵設(shè)備”和購(gòu)電區(qū)域內(nèi)部節(jié)點(diǎn)之間的功率轉(zhuǎn)移分布因子矩陣;SC,C為售電區(qū)域“關(guān)鍵設(shè)備”和售電區(qū)域內(nèi)部節(jié)點(diǎn)之間的功率轉(zhuǎn)移分布因子矩陣;fˉt為聯(lián)絡(luò)線容量限制向量;fˉR、fˉC分別為購(gòu)電區(qū)域、售電區(qū)域“關(guān)鍵設(shè)備”容量限制向量;-qR、qˉR分別為購(gòu)電區(qū)域的投標(biāo)數(shù)量的下限、上限向量;-qC、qˉC分別為售電區(qū)域的投標(biāo)數(shù)量的下限、上限向量。

    在需求總量方面,跨區(qū)域電力現(xiàn)貨交易主要針對(duì)用電缺口或發(fā)電余量進(jìn)行優(yōu)化配置。在設(shè)定缺口或余量投標(biāo)上、下限值時(shí),如式(8)和式(9)所示,已經(jīng)考慮了有關(guān)電網(wǎng)的需求總量,不再需要對(duì)各區(qū)域電網(wǎng)的需求總量進(jìn)行設(shè)定。在需求彈性方面,各市場(chǎng)參與者提交的電力-價(jià)格曲線體現(xiàn)了供需各方出清量對(duì)價(jià)格的彈性。

    與引言中介紹的美國(guó)CTS 機(jī)制、GCTS 機(jī)制相比,本文設(shè)計(jì)的機(jī)制不依賴于跨區(qū)域投標(biāo)者,解決了因跨區(qū)域投標(biāo)者套利導(dǎo)致的跨區(qū)域交易不充分的問(wèn)題。與歐洲電力市場(chǎng)采用的ATC 機(jī)制相比,所提方法能夠充分利用輸電線路容量,提升經(jīng)濟(jì)性。與歐洲采用的Flow-based 機(jī)制相比,所提方法基于網(wǎng)絡(luò)等值方法后的簡(jiǎn)化公共網(wǎng)絡(luò)完成出清,簡(jiǎn)化了投標(biāo)流程,保護(hù)了各區(qū)域內(nèi)部網(wǎng)架信息隱私。

    3 定價(jià)機(jī)制

    本章對(duì)省間電力現(xiàn)貨交易的定價(jià)機(jī)制進(jìn)行說(shuō)明。對(duì)于省i中的區(qū)域間市場(chǎng)參與者k,它的支付金額Mi,k由下式計(jì)算得到:

    式中:Hi,k、Li,k分別為市場(chǎng)清算支付、所需承擔(dān)的輸電費(fèi)用。這兩部分費(fèi)用的計(jì)算分別在3.1 節(jié)和3.2節(jié)中進(jìn)行闡述。

    3.1 計(jì)算分區(qū)LMP

    在分區(qū)邊際定價(jià)的規(guī)則下,每個(gè)區(qū)域的市場(chǎng)出清價(jià)格等于該區(qū)域的邊際投標(biāo)價(jià)格[20]。以圖2 中的兩區(qū)域電力系統(tǒng)為例進(jìn)行說(shuō)明,其市場(chǎng)出清過(guò)程如圖4 所示。圖中:星號(hào)表示邊際投標(biāo)。

    圖4 區(qū)域i 和j 的市場(chǎng)出清價(jià)格Fig.4 Market clearing prices of region i and region j

    由出清模型式(3)—式(7)對(duì)應(yīng)的Karush-Kuhn-Tucker(KKT)條件,區(qū)域i和區(qū)域j的區(qū)域邊際價(jià)格用λi和λj表示,計(jì)算公式分別如式(11)和式(12)所示。

    式中:pi,m、pj,m分別為區(qū)域i、區(qū)域j的邊際投標(biāo)價(jià)格。

    區(qū)域i中的區(qū)域間市場(chǎng)參與者k所需支付的市場(chǎng)出清金額Hi,k可由式(13)計(jì)算得到。

    式中:qi,k為省i中的區(qū)域間市場(chǎng)參與者k的出清量。區(qū)域j的計(jì)算過(guò)程類似,此處省略。

    3.2 計(jì)算跨省跨區(qū)輸電價(jià)格

    本節(jié)計(jì)算各筆省間現(xiàn)貨交易所需支付的輸電費(fèi)用。

    首先,將簡(jiǎn)化公共網(wǎng)絡(luò)的傳輸價(jià)格向量lC表示如下:

    式中:lP、lR、lI分別為省、地區(qū)(如:華北電網(wǎng)、華東電網(wǎng)等)、區(qū)域間輸電價(jià)格。中國(guó)各省、地區(qū)和跨地區(qū)輸電價(jià)格詳見(jiàn)文獻(xiàn)[25-26],對(duì)每個(gè)省電網(wǎng)、區(qū)域電網(wǎng)、每條區(qū)域間聯(lián)絡(luò)線分別核定其輸電價(jià)格。

    其次,為了量化區(qū)域間交易對(duì)輸電線路的潮流貢獻(xiàn),需要明確每筆省間交易的源節(jié)點(diǎn)和匯節(jié)點(diǎn),進(jìn)而將出清的投標(biāo)數(shù)量乘以等效功率注入和線路之間的功率轉(zhuǎn)移分布因子以獲得潮流貢獻(xiàn)。在完成市場(chǎng)出清后,通過(guò)如圖5 所示的等效雙邊交易求解方法得到等效的雙邊交易,明確各省間現(xiàn)貨交易的源匯節(jié)點(diǎn)。圖中:矩陣S、矩陣D分別包含電力賣出區(qū)域和電力購(gòu)入?yún)^(qū)域的投標(biāo)信息,其中矩陣的1—4 列向量分別是價(jià)格、數(shù)量、節(jié)點(diǎn)編號(hào)、區(qū)域編號(hào);矩陣C包含等價(jià)雙邊交換信息,其中矩陣的1—6 列分別為數(shù)量、源匯節(jié)點(diǎn)價(jià)差、源節(jié)點(diǎn)數(shù)、匯節(jié)點(diǎn)數(shù)、電力購(gòu)入?yún)^(qū)域編號(hào)、電力賣出區(qū)域編號(hào)。

    圖5 等效雙邊交易求解方法Fig.5 Solving method of equivalent bilateral trading

    以圖1 所示的兩區(qū)域電力系統(tǒng)為例。根據(jù)圖5所示的算法流程,可以基于高低配對(duì)的原則形成5 筆 雙 邊 交 易,分 別 為qi,1-qj,1、qi,2-qj,1、qi,2-qj,2、qi,3-qj,2、qi,3-qj,3。假設(shè)區(qū)域i的市場(chǎng)參與者k(購(gòu)入電力)和區(qū)域j的市場(chǎng)參與者α(賣出電力)之間存在省間現(xiàn)貨交易,其數(shù)量用qk,α表示。值得注意的是,目前中國(guó)的跨區(qū)跨省輸電費(fèi)用均由電力購(gòu)買者一側(cè)承擔(dān)[19]。在此條件下,市場(chǎng)參與者k承擔(dān)的輸電費(fèi)用Li,k為:

    式中:li、lj分別為區(qū)域i、區(qū)域j的輸電價(jià)格;li,j為區(qū)域i與區(qū)域j間跨區(qū)域輸電線路的輸電價(jià)格向量;qk,α為區(qū)域i的市場(chǎng)參與者k(購(gòu)入電力)和區(qū)域j的市場(chǎng)參與者α(賣出電力)之間存在的省間現(xiàn)貨交易數(shù)量。式中向量B與e的含義與式(1)中相同,詳見(jiàn)1.2 節(jié)中對(duì)式(1)各變量的相關(guān)解釋說(shuō)明。

    市場(chǎng)參與者α承擔(dān)的輸電費(fèi)用Lj,α可由式(16)計(jì)算得到。

    在輸電費(fèi)用計(jì)算方面,相較于目前中國(guó)國(guó)內(nèi)省間現(xiàn)貨市場(chǎng)中基于合同路徑的輸電費(fèi)用計(jì)算方法[19],本文提出的方法充分考慮跨區(qū)域交易對(duì)區(qū)域間聯(lián)絡(luò)線和區(qū)域內(nèi)線路的物理潮流影響,采用潮流轉(zhuǎn)移分布因子對(duì)該影響進(jìn)行量化,從而按潮流貢獻(xiàn)合理收取輸電費(fèi)用。該方法解決了合同路徑法由于物理潮流對(duì)合同潮流的偏離無(wú)法實(shí)現(xiàn)“誰(shuí)受益,誰(shuí)承擔(dān)”的問(wèn)題[27]。

    4 算例驗(yàn)證

    采用兩省系統(tǒng)和三省系統(tǒng)對(duì)本文所提出的中國(guó)省間現(xiàn)貨交易優(yōu)化設(shè)計(jì)與定價(jià)機(jī)制進(jìn)行驗(yàn)證,相關(guān)算例在Intel Core i7 1.10 GHz 處理器和16 GB 內(nèi)存的PC 端完成。

    4.1 算例1:兩省(6 節(jié)點(diǎn))系統(tǒng)

    首先,采用兩省系統(tǒng)(6 節(jié)點(diǎn))來(lái)驗(yàn)證所提出的省間電力現(xiàn)貨市場(chǎng)出清與定價(jià)機(jī)制,如圖6 所示。選擇邊界上的5 號(hào)節(jié)點(diǎn)作為全局系統(tǒng)的相角參考節(jié)點(diǎn),所有內(nèi)部線路和聯(lián)絡(luò)線的電抗均設(shè)置為1.0 p.u.,假定所有輸電線路的容量為無(wú)窮大,省A、省B、聯(lián)絡(luò)線2-4、聯(lián)絡(luò)線3-5 的輸電價(jià)格分別為lA、lB、l2-4、l3-5。在此算例中,假設(shè)省B 從省A 進(jìn)口電力。省間電力現(xiàn)貨交易投標(biāo)情況如表1 所示。

    表1 算例1 中跨區(qū)域投標(biāo)參數(shù)Table 1 Cross-regional bidding parameters in case 1

    圖6 算例1 的兩省(6 節(jié)點(diǎn))系統(tǒng)Fig.6 Two-province (6-bus) power system in case 1

    其次,對(duì)本文提出的市場(chǎng)出清模型(式(3)—式(7))進(jìn)行求解,計(jì)算得到投標(biāo)1、投標(biāo)2、投標(biāo)3、投標(biāo)4 的出清量分別為10、20、20、10 MW,此時(shí),省A 與省B 之間的跨區(qū)域交易量為30 MW,兩省的邊際出清的交易分別為投標(biāo)2 和4。因此,根據(jù)式(11)—式(13),兩省跨區(qū)域交易的市場(chǎng)出清邊際λ和投標(biāo)1、投標(biāo)2、投標(biāo)3、投標(biāo)4 的支付費(fèi)用H在表2 中給出。

    表2 算例1 中市場(chǎng)出清價(jià)格和支付Table 2 Market clearing prices and payments in case 1

    然后,對(duì)投標(biāo)1~4 所需承擔(dān)的輸電費(fèi)用進(jìn)行計(jì)算。根據(jù)圖5 所示的等效雙邊交易求解方法,可以得到等效的雙邊交易分別為q1-q3、q2-q3、q2-q4,這三筆等效雙邊交易的出清量均為10 MW。表3 給出了等效雙邊交易對(duì)聯(lián)絡(luò)線的潮流貢獻(xiàn)δ,根據(jù)式(15)和式(16)可以計(jì)算得到區(qū)域間交易需承擔(dān)的輸電費(fèi)用(見(jiàn)表4)。

    表3 算例1 中等效雙邊交易對(duì)聯(lián)絡(luò)線的潮流貢獻(xiàn)Table 3 Contribution of equivalent bilateral transactions to power flow of tie lines in case 1

    表4 算例1 中跨區(qū)域輸電費(fèi)用Table 4 Cross-regional transmission tariffs in case 1

    根據(jù)3.2 節(jié),因?yàn)橥稑?biāo)1 和投標(biāo)2 是電力賣出方,這兩個(gè)投標(biāo)不承擔(dān)跨省跨區(qū)輸電費(fèi)用,而電力購(gòu)買方投標(biāo)3 和投標(biāo)4 支付每個(gè)等價(jià)雙邊交易中的輸電費(fèi)用。因此,投標(biāo)1、投標(biāo)2、投標(biāo)3、投標(biāo)4 承擔(dān)的輸電費(fèi)用分別為0、0、650、325 元/h。而省A運(yùn)營(yíng)商、省B 運(yùn)營(yíng)商和協(xié)調(diào)者的輸電收入分別為150、150、675 元/h。最后,投標(biāo)1、投標(biāo)2、投標(biāo)3、投標(biāo)4 的總支付金額分別為-2 000、-4 000、8 650、4 325 元/h,本算例的計(jì)算結(jié)果如圖7 所示。

    圖7 算例1 中市場(chǎng)出清支付和輸電費(fèi)用Fig.7 Market clearing payments and transmission tariffs in case 1

    此外,將本文所提的市場(chǎng)出清方法與JED、GCTS、ATC 機(jī)制在跨區(qū)域交易量方面進(jìn)行對(duì)比。采用本文所提的市場(chǎng)出清方法、JED、GCTS、ATC 機(jī)制得到的跨區(qū)域交易量分別為30、30、28、25 MW??梢钥闯?本文方法在簡(jiǎn)化出清流程、保護(hù)各區(qū)域內(nèi)部網(wǎng)架隱私的前提下,求解得到的跨區(qū)域交易量等于JED 求解的交易量,驗(yàn)證了第2 章的結(jié)論。

    將所提的輸電費(fèi)用計(jì)算方法所得結(jié)果與當(dāng)前基于合同路徑的計(jì)算方法所得結(jié)果進(jìn)行對(duì)比,結(jié)果如表5 所示。可以看出,本文所提輸電費(fèi)用計(jì)算方法準(zhǔn)確考慮了跨區(qū)域交易對(duì)省間聯(lián)絡(luò)線的物理潮流貢獻(xiàn),驗(yàn)證了第3 章的結(jié)論。

    表5 算例1 中不同輸電費(fèi)用計(jì)算方法對(duì)比Table 5 Comparison of transmission tariffs calculation methods in case 1

    最后,對(duì)比考慮完整網(wǎng)絡(luò)架構(gòu)和考慮網(wǎng)絡(luò)等值后的潮流計(jì)算結(jié)果(見(jiàn)表6)和交易出清結(jié)果(見(jiàn)表7)??梢钥闯?網(wǎng)絡(luò)等值后與基于完整網(wǎng)絡(luò)的計(jì)算結(jié)果相同,驗(yàn)證了所提方法的合理性。

    表6 算例1 中網(wǎng)絡(luò)等值前后輸電線路潮流對(duì)比Table 6 Comparison of power flow of transmission lines before and after network equivalence in case 1

    表7 算例1 中網(wǎng)絡(luò)等值前后交易量對(duì)比Table 7 Comparison of trading volumes before and after network equivalence in case 1

    4.2 算例2:三省系統(tǒng)(118 節(jié)點(diǎn))

    將本文所提方法擴(kuò)展到三省系統(tǒng)(118 節(jié)點(diǎn)),該系統(tǒng)如圖8(a)所示[24]。應(yīng)用網(wǎng)絡(luò)等值方法后得到的簡(jiǎn)化公共網(wǎng)絡(luò)如圖8(b)所示。圖中:紅色虛線表示各區(qū)域?qū)ν獾戎稻€路。

    圖8 三省系統(tǒng)(118 節(jié)點(diǎn))Fig.8 Three-province power system (118-bus)

    假設(shè)省1 賣出電力,省2 和3 有進(jìn)口電力需求,表8 給出了省間現(xiàn)貨投標(biāo)信息。聯(lián)絡(luò)線30-3 和聯(lián)絡(luò)線68-81 的容量分別設(shè)置為50 MW 和25 MW,其他輸電線路則設(shè)置為無(wú)窮大。省1、省2、省3 的輸電價(jià)格分別為200、300、400 元/(MW·h),表9 給出了算例2 中聯(lián)絡(luò)線的輸電價(jià)格。

    表8 算例2 中跨區(qū)域投標(biāo)參數(shù)Table 8 Cross-regional bidding parameters in case 2

    表9 算例2 中聯(lián)絡(luò)線輸電價(jià)格Table 9 Transmission prices of tie-lines in case 2

    首先,根據(jù)式(8)—式(12)計(jì)算得到市場(chǎng)出清結(jié)果如表10 所示?;诜謪^(qū)邊際定價(jià)的方法,通過(guò)列寫模型的KKT 條件計(jì)算每個(gè)省的邊際價(jià)格λ,并用其來(lái)結(jié)算該省的投標(biāo)。在此算例中,投標(biāo)1~10的邊際投標(biāo)價(jià)格為300 元/(MW·h),而投標(biāo)11~15的邊際投標(biāo)價(jià)格為650 元/(MW·h),省間現(xiàn)貨交易1~15 在市場(chǎng)出清過(guò)程中需要支付的費(fèi)用如表10 所示,用符號(hào)H表示。

    表10 算例2 中跨區(qū)域交易出清結(jié)果Table 10 Market clearing results of cross-regional trading in case 2

    其次,根據(jù)圖5 所示的算法流程,可以計(jì)算得到省間投標(biāo)的等效雙邊交易,結(jié)果如表11 所示。根據(jù)式(15)和式(16),計(jì)算得到等效雙邊交易需要承擔(dān)的輸電費(fèi)用。輸電費(fèi)用在中國(guó)主要由電力購(gòu)入端進(jìn)行分?jǐn)?投標(biāo)1~15 所需承擔(dān)的輸電費(fèi)用如表12 所示,圖9 給出了各筆省間現(xiàn)貨交易的總支付金額。綜上,本文所提出的市場(chǎng)清算和結(jié)算方式在此三省系統(tǒng)(118 節(jié)點(diǎn))中得到了驗(yàn)證。

    表11 算例2 中等效雙邊交易Table 11 Equivalent bilateral tradings in case 2

    表12 算例2 中跨區(qū)域投標(biāo)的輸電費(fèi)用Table 12 Transmission tariffs of cross-regional biddings in case 2

    圖9 算例2 中市場(chǎng)出清支付和輸電費(fèi)用Fig.9 Market clearing payments and transmission tariffs in case 2

    最后,對(duì)比算例2 中不同跨區(qū)域現(xiàn)貨交易機(jī)制出清量,本文所提方法跨區(qū)域交易量為92.1 MW,JED、GCTS、ATC 的 跨 區(qū) 域 交 易 量 分 別 為92.1、88.3、82.5 MW。輸電費(fèi)用計(jì)算方面的對(duì)比結(jié)果如表13 所示,分別驗(yàn)證了第2 章和第4 章分析結(jié)果的準(zhǔn)確性。

    表13 算例2 中不同輸電費(fèi)用計(jì)算方法對(duì)比Table 13 Comparison of different transmission tariff calculation methods in case 2

    5 結(jié)語(yǔ)

    在實(shí)現(xiàn)“碳達(dá)峰·碳中和”目標(biāo)的背景下,發(fā)展完善省間電力現(xiàn)貨交易有利于實(shí)現(xiàn)更大范圍內(nèi)的資源優(yōu)化配置,助推中國(guó)綠色低碳發(fā)展進(jìn)程。本文面向多區(qū)域互聯(lián)電力系統(tǒng),研究了適用于省間現(xiàn)貨交易的出清流程與定價(jià)機(jī)制。首先,根據(jù)地理或行政因素劃分報(bào)價(jià)區(qū),并基于網(wǎng)絡(luò)等值得到簡(jiǎn)化公共網(wǎng)絡(luò)。其次,設(shè)計(jì)省間現(xiàn)貨市場(chǎng)出清方法,通過(guò)求解出清模型得到各市場(chǎng)參與者投標(biāo)的出清數(shù)量。最后,采用區(qū)域邊際價(jià)格對(duì)每個(gè)區(qū)域進(jìn)行結(jié)算,并基于等效雙邊交易匹配原則計(jì)算省間現(xiàn)貨交易所應(yīng)承擔(dān)的輸電費(fèi)用。

    下一步,將在省間現(xiàn)貨交易中進(jìn)一步考慮不同形式輸電成本分?jǐn)倷C(jī)制(如容量費(fèi)用等)的應(yīng)用。

    猜你喜歡
    區(qū)域
    分割區(qū)域
    探尋區(qū)域創(chuàng)新的密碼
    科學(xué)(2020年5期)2020-11-26 08:19:22
    基于BM3D的復(fù)雜紋理區(qū)域圖像去噪
    軟件(2020年3期)2020-04-20 01:45:18
    小區(qū)域、大發(fā)展
    商周刊(2018年15期)2018-07-27 01:41:20
    論“戎”的活動(dòng)區(qū)域
    區(qū)域發(fā)展篇
    區(qū)域經(jīng)濟(jì)
    關(guān)于四色猜想
    分區(qū)域
    公司治理與技術(shù)創(chuàng)新:分區(qū)域比較
    汤姆久久久久久久影院中文字幕| 成人二区视频| av视频免费观看在线观看| 国产免费现黄频在线看| 精品国产一区二区久久| 久久99热这里只频精品6学生| 狂野欧美激情性bbbbbb| 精品一品国产午夜福利视频| 亚洲第一区二区三区不卡| 国产精品国产三级专区第一集| 日本vs欧美在线观看视频| 9色porny在线观看| 蜜桃国产av成人99| 欧美人与善性xxx| freevideosex欧美| 搡老乐熟女国产| 成年人午夜在线观看视频| 亚洲图色成人| 人体艺术视频欧美日本| 最近最新中文字幕免费大全7| 91精品国产国语对白视频| 成人二区视频| 看十八女毛片水多多多| 国产伦精品一区二区三区视频9| 日本av手机在线免费观看| 一边摸一边做爽爽视频免费| 久久av网站| 亚洲精品自拍成人| 一级二级三级毛片免费看| 国产黄色视频一区二区在线观看| 国产淫语在线视频| 午夜日本视频在线| 久久综合国产亚洲精品| 国产精品蜜桃在线观看| 亚洲人成77777在线视频| 精品国产乱码久久久久久小说| 青青草视频在线视频观看| 九色亚洲精品在线播放| 黄色欧美视频在线观看| 母亲3免费完整高清在线观看 | 久久久欧美国产精品| 边亲边吃奶的免费视频| 男女国产视频网站| 成人二区视频| 国产精品一区二区在线观看99| 精品久久久久久久久亚洲| 综合色丁香网| 午夜精品国产一区二区电影| 亚洲成人av在线免费| 日韩制服骚丝袜av| 我要看黄色一级片免费的| 精品人妻一区二区三区麻豆| 波野结衣二区三区在线| 爱豆传媒免费全集在线观看| 男人爽女人下面视频在线观看| tube8黄色片| 肉色欧美久久久久久久蜜桃| 九色成人免费人妻av| 女性被躁到高潮视频| 国产精品麻豆人妻色哟哟久久| 插逼视频在线观看| 亚州av有码| 高清欧美精品videossex| 丁香六月天网| 中文天堂在线官网| 久久久精品94久久精品| 免费黄色在线免费观看| 97超碰精品成人国产| 午夜日本视频在线| 色视频在线一区二区三区| 国产无遮挡羞羞视频在线观看| 亚洲性久久影院| av播播在线观看一区| 国产成人av激情在线播放 | 永久网站在线| 2018国产大陆天天弄谢| 99热全是精品| 看十八女毛片水多多多| 亚洲av福利一区| 国产老妇伦熟女老妇高清| 中文字幕久久专区| 日韩人妻高清精品专区| 观看美女的网站| 蜜桃国产av成人99| 免费av不卡在线播放| 精品久久久久久电影网| 能在线免费看毛片的网站| 高清在线视频一区二区三区| 日本免费在线观看一区| 男女高潮啪啪啪动态图| 韩国av在线不卡| 日韩成人伦理影院| 少妇人妻 视频| 精品卡一卡二卡四卡免费| 久久精品国产鲁丝片午夜精品| 黑丝袜美女国产一区| 中国三级夫妇交换| 国产成人aa在线观看| 国产午夜精品一二区理论片| 成人黄色视频免费在线看| av国产精品久久久久影院| 中文字幕制服av| 伦理电影大哥的女人| 草草在线视频免费看| 精品久久久精品久久久| 人妻 亚洲 视频| 免费高清在线观看视频在线观看| 精品卡一卡二卡四卡免费| 久久久a久久爽久久v久久| 精品熟女少妇av免费看| 国产精品一区二区在线不卡| 久热久热在线精品观看| 肉色欧美久久久久久久蜜桃| 欧美性感艳星| 久久久久久久久久久丰满| 国产精品一区二区在线观看99| 18禁裸乳无遮挡动漫免费视频| 国产精品嫩草影院av在线观看| 99久久精品一区二区三区| 婷婷色麻豆天堂久久| 国产综合精华液| 在线亚洲精品国产二区图片欧美 | 国产老妇伦熟女老妇高清| 午夜视频国产福利| 免费久久久久久久精品成人欧美视频 | 中国国产av一级| 妹子高潮喷水视频| 99视频精品全部免费 在线| 在线观看三级黄色| 欧美成人精品欧美一级黄| 亚洲成人av在线免费| 国产男人的电影天堂91| a级毛色黄片| 成年人免费黄色播放视频| 国产一区二区在线观看av| 男男h啪啪无遮挡| 国产亚洲精品第一综合不卡 | 免费av中文字幕在线| 免费看av在线观看网站| 最近中文字幕2019免费版| 国产一区亚洲一区在线观看| 人妻系列 视频| www.av在线官网国产| 人人妻人人澡人人看| 一本久久精品| 在线播放无遮挡| 日本色播在线视频| 中文乱码字字幕精品一区二区三区| 免费看不卡的av| av福利片在线| 亚洲欧美清纯卡通| 熟女av电影| 亚洲精品日韩av片在线观看| 精品国产国语对白av| 极品少妇高潮喷水抽搐| 亚洲av中文av极速乱| 久久久精品94久久精品| 国产精品久久久久久精品古装| 国产探花极品一区二区| 九九久久精品国产亚洲av麻豆| 91精品国产国语对白视频| 久久久a久久爽久久v久久| 伊人亚洲综合成人网| 大码成人一级视频| 狂野欧美白嫩少妇大欣赏| 国产成人一区二区在线| 美女大奶头黄色视频| 激情五月婷婷亚洲| 18禁观看日本| av网站免费在线观看视频| 国产亚洲午夜精品一区二区久久| 免费看av在线观看网站| 国产有黄有色有爽视频| 大片电影免费在线观看免费| 日韩 亚洲 欧美在线| 国产精品人妻久久久久久| 婷婷色麻豆天堂久久| 黄色怎么调成土黄色| 亚洲av不卡在线观看| 老司机影院毛片| 亚洲精品乱久久久久久| 免费黄色在线免费观看| 亚洲av二区三区四区| 亚洲精品乱久久久久久| 最近的中文字幕免费完整| 国产精品不卡视频一区二区| 国产伦精品一区二区三区视频9| 一区二区三区乱码不卡18| 大又大粗又爽又黄少妇毛片口| 婷婷色麻豆天堂久久| 中文字幕av电影在线播放| 日日啪夜夜爽| 国产精品国产三级专区第一集| 国产成人免费无遮挡视频| 乱码一卡2卡4卡精品| 国产免费福利视频在线观看| 欧美性感艳星| 欧美日韩亚洲高清精品| 亚洲欧美精品自产自拍| 狂野欧美激情性xxxx在线观看| 亚洲成人一二三区av| 成人免费观看视频高清| 亚洲精品色激情综合| 在线免费观看不下载黄p国产| 看免费成人av毛片| 中文精品一卡2卡3卡4更新| 精品人妻熟女av久视频| 欧美激情极品国产一区二区三区 | av一本久久久久| 亚洲精品中文字幕在线视频| 一级毛片黄色毛片免费观看视频| 国产片内射在线| 三上悠亚av全集在线观看| 欧美xxⅹ黑人| 观看美女的网站| 国产爽快片一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 天堂8中文在线网| 十八禁高潮呻吟视频| 肉色欧美久久久久久久蜜桃| 天美传媒精品一区二区| 国产高清三级在线| 哪个播放器可以免费观看大片| 好男人视频免费观看在线| 免费少妇av软件| 2022亚洲国产成人精品| 亚洲国产最新在线播放| 久久久精品免费免费高清| 国产精品成人在线| 日本欧美国产在线视频| 欧美+日韩+精品| 亚洲熟女精品中文字幕| 在现免费观看毛片| 欧美激情国产日韩精品一区| 大码成人一级视频| 亚洲人成网站在线播| 成人毛片a级毛片在线播放| 亚洲精品av麻豆狂野| 大又大粗又爽又黄少妇毛片口| 国产欧美日韩综合在线一区二区| 日韩一区二区三区影片| 纵有疾风起免费观看全集完整版| av电影中文网址| 不卡视频在线观看欧美| 国产国拍精品亚洲av在线观看| 中文字幕免费在线视频6| 日韩一本色道免费dvd| 成人无遮挡网站| 中国美白少妇内射xxxbb| 最近最新中文字幕免费大全7| 欧美亚洲日本最大视频资源| 在线看a的网站| 国产av国产精品国产| 在线观看免费视频网站a站| 亚洲欧美成人综合另类久久久| 最近的中文字幕免费完整| 国产有黄有色有爽视频| 一区二区三区免费毛片| 午夜免费男女啪啪视频观看| 精品午夜福利在线看| 亚洲av不卡在线观看| 亚洲激情五月婷婷啪啪| 国产老妇伦熟女老妇高清| 亚洲成色77777| 欧美日韩视频高清一区二区三区二| 欧美一级a爱片免费观看看| 久久精品国产自在天天线| 一级黄片播放器| 肉色欧美久久久久久久蜜桃| 丝袜美足系列| www.av在线官网国产| 日韩强制内射视频| 赤兔流量卡办理| 免费不卡的大黄色大毛片视频在线观看| 欧美亚洲日本最大视频资源| 国产精品人妻久久久影院| 蜜桃国产av成人99| 婷婷色麻豆天堂久久| 国产精品嫩草影院av在线观看| 色94色欧美一区二区| 99九九线精品视频在线观看视频| 午夜视频国产福利| 精品久久国产蜜桃| 国产色婷婷99| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av在线观看美女高潮| 国产精品成人在线| 韩国av在线不卡| 午夜福利网站1000一区二区三区| 夜夜骑夜夜射夜夜干| 黄色配什么色好看| 9色porny在线观看| 熟女av电影| 国产亚洲午夜精品一区二区久久| 免费人成在线观看视频色| 日本vs欧美在线观看视频| 欧美日韩亚洲高清精品| 国产 一区精品| 久久精品夜色国产| 亚洲一区二区三区欧美精品| 午夜激情av网站| 全区人妻精品视频| 日韩欧美精品免费久久| 国产成人一区二区在线| 9色porny在线观看| 国产成人一区二区在线| 国产精品女同一区二区软件| 高清在线视频一区二区三区| 男女啪啪激烈高潮av片| 成人国产麻豆网| 免费黄色在线免费观看| av国产久精品久网站免费入址| 18禁观看日本| 老司机影院毛片| 久久久久久久大尺度免费视频| 成人午夜精彩视频在线观看| 只有这里有精品99| a级毛片免费高清观看在线播放| 插阴视频在线观看视频| 大码成人一级视频| 大码成人一级视频| 大香蕉97超碰在线| 日本欧美视频一区| 91国产中文字幕| 午夜91福利影院| 91成人精品电影| 国产精品国产三级专区第一集| 国产成人免费无遮挡视频| 久久久a久久爽久久v久久| 91午夜精品亚洲一区二区三区| 中文字幕av电影在线播放| 日韩精品免费视频一区二区三区 | 亚洲av免费高清在线观看| 91aial.com中文字幕在线观看| 国语对白做爰xxxⅹ性视频网站| 久久狼人影院| 亚洲av日韩在线播放| 亚洲,一卡二卡三卡| 人妻 亚洲 视频| 成人国产av品久久久| 国模一区二区三区四区视频| 久久久欧美国产精品| 少妇人妻 视频| 韩国av在线不卡| 熟女电影av网| 在线观看免费高清a一片| 日韩欧美一区视频在线观看| 美女国产视频在线观看| 国产成人aa在线观看| 久久毛片免费看一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 特大巨黑吊av在线直播| 国产黄频视频在线观看| 成年人午夜在线观看视频| 亚洲成色77777| 人人澡人人妻人| 在线观看免费视频网站a站| 成人国产麻豆网| 久久99热6这里只有精品| 狠狠婷婷综合久久久久久88av| 亚洲欧洲精品一区二区精品久久久 | 国产精品一区二区在线不卡| 97在线人人人人妻| 国产精品偷伦视频观看了| 国产精品嫩草影院av在线观看| 成人免费观看视频高清| av天堂久久9| 国产精品麻豆人妻色哟哟久久| 久久久久精品久久久久真实原创| 韩国av在线不卡| 大话2 男鬼变身卡| av卡一久久| 亚洲无线观看免费| 各种免费的搞黄视频| 国产无遮挡羞羞视频在线观看| 日本-黄色视频高清免费观看| 啦啦啦中文免费视频观看日本| 一级毛片aaaaaa免费看小| 国产熟女午夜一区二区三区 | 99精国产麻豆久久婷婷| 大码成人一级视频| 99热全是精品| 国产男女超爽视频在线观看| 国产亚洲av片在线观看秒播厂| 国产亚洲一区二区精品| 国产av码专区亚洲av| 国产综合精华液| 欧美97在线视频| 最黄视频免费看| 久久精品夜色国产| 天天操日日干夜夜撸| 亚洲精品国产av成人精品| 看免费成人av毛片| 一级a做视频免费观看| 乱人伦中国视频| 日本黄大片高清| 国产一区亚洲一区在线观看| 99热全是精品| 男的添女的下面高潮视频| 中文字幕av电影在线播放| 免费人成在线观看视频色| 一区二区av电影网| 国产乱来视频区| 午夜福利在线观看免费完整高清在| 国产亚洲最大av| videos熟女内射| 五月伊人婷婷丁香| a 毛片基地| 精品国产一区二区三区久久久樱花| 成人影院久久| 男人添女人高潮全过程视频| 一区在线观看完整版| 亚洲欧美中文字幕日韩二区| 国产欧美亚洲国产| 男女高潮啪啪啪动态图| 中文字幕人妻熟人妻熟丝袜美| 免费观看a级毛片全部| 亚洲精品成人av观看孕妇| 夫妻性生交免费视频一级片| 日韩强制内射视频| 欧美精品高潮呻吟av久久| 精品久久久久久久久亚洲| 最近的中文字幕免费完整| 国产精品久久久久久久电影| 日本爱情动作片www.在线观看| 欧美精品高潮呻吟av久久| 观看av在线不卡| 中文字幕免费在线视频6| 国产成人91sexporn| 精品久久久噜噜| 午夜久久久在线观看| 国产精品不卡视频一区二区| 久久国产亚洲av麻豆专区| 精品国产一区二区三区久久久樱花| 日韩欧美精品免费久久| 国产成人精品福利久久| 久久久久视频综合| 国产精品国产三级国产av玫瑰| 99热6这里只有精品| 色94色欧美一区二区| 在线观看国产h片| 大又大粗又爽又黄少妇毛片口| 国产成人91sexporn| 嘟嘟电影网在线观看| 我的老师免费观看完整版| 美女内射精品一级片tv| 久久久久久久久久成人| 午夜久久久在线观看| 午夜福利视频在线观看免费| 国产 一区精品| 韩国av在线不卡| 交换朋友夫妻互换小说| 国产免费视频播放在线视频| 久久亚洲国产成人精品v| 男女边吃奶边做爰视频| videos熟女内射| 一区二区av电影网| 中文字幕亚洲精品专区| 91aial.com中文字幕在线观看| av在线播放精品| 看十八女毛片水多多多| 国产精品嫩草影院av在线观看| 婷婷色麻豆天堂久久| 18禁动态无遮挡网站| 色吧在线观看| 考比视频在线观看| 久久97久久精品| 欧美日韩一区二区视频在线观看视频在线| 久久精品夜色国产| 91国产中文字幕| 亚洲av在线观看美女高潮| 国产国拍精品亚洲av在线观看| 尾随美女入室| 能在线免费看毛片的网站| 欧美精品人与动牲交sv欧美| 亚洲精品中文字幕在线视频| 十分钟在线观看高清视频www| 一级毛片aaaaaa免费看小| 在线观看美女被高潮喷水网站| 波野结衣二区三区在线| 水蜜桃什么品种好| 最黄视频免费看| 欧美精品高潮呻吟av久久| 久久ye,这里只有精品| 日韩,欧美,国产一区二区三区| 国语对白做爰xxxⅹ性视频网站| 天天躁夜夜躁狠狠久久av| 如何舔出高潮| 各种免费的搞黄视频| 亚洲精华国产精华液的使用体验| 啦啦啦视频在线资源免费观看| 夜夜看夜夜爽夜夜摸| 少妇熟女欧美另类| 亚洲av成人精品一二三区| 久久国产精品男人的天堂亚洲 | 嫩草影院入口| 天天躁夜夜躁狠狠久久av| 亚洲国产欧美日韩在线播放| 国产精品人妻久久久影院| 欧美激情 高清一区二区三区| 美女内射精品一级片tv| 91午夜精品亚洲一区二区三区| 国产精品久久久久久久久免| 色吧在线观看| 人妻一区二区av| 国语对白做爰xxxⅹ性视频网站| 视频在线观看一区二区三区| 亚洲人成网站在线播| 国产国拍精品亚洲av在线观看| 在线精品无人区一区二区三| 一本久久精品| 国产精品成人在线| 少妇精品久久久久久久| 18禁在线播放成人免费| 国产精品久久久久久久电影| 久久久精品区二区三区| 亚洲av.av天堂| 赤兔流量卡办理| 精品少妇内射三级| 男女边摸边吃奶| 亚洲精品国产av蜜桃| 国产成人精品婷婷| 九草在线视频观看| a 毛片基地| 日韩精品有码人妻一区| 一区二区日韩欧美中文字幕 | 春色校园在线视频观看| 超色免费av| 国产亚洲av片在线观看秒播厂| 在线观看国产h片| 插逼视频在线观看| 国产69精品久久久久777片| 成人国语在线视频| 91精品三级在线观看| 国产精品熟女久久久久浪| 成人手机av| 少妇精品久久久久久久| 精品少妇久久久久久888优播| 寂寞人妻少妇视频99o| 妹子高潮喷水视频| 精品午夜福利在线看| 国产一区二区三区综合在线观看 | 久久精品熟女亚洲av麻豆精品| 日本爱情动作片www.在线观看| 三级国产精品欧美在线观看| 纯流量卡能插随身wifi吗| 日本猛色少妇xxxxx猛交久久| 免费大片黄手机在线观看| 国产精品嫩草影院av在线观看| 哪个播放器可以免费观看大片| 精品一区二区三卡| 少妇猛男粗大的猛烈进出视频| a级毛片在线看网站| 亚洲av在线观看美女高潮| 国产精品一区二区在线不卡| 你懂的网址亚洲精品在线观看| 热re99久久精品国产66热6| 国产精品人妻久久久影院| 国产在线视频一区二区| 丝袜脚勾引网站| 日韩精品有码人妻一区| 久久毛片免费看一区二区三区| 能在线免费看毛片的网站| 成人手机av| 男人添女人高潮全过程视频| 在线观看免费高清a一片| 视频在线观看一区二区三区| 精品亚洲成国产av| 精品少妇久久久久久888优播| 交换朋友夫妻互换小说| 乱人伦中国视频| 美女xxoo啪啪120秒动态图| 日本-黄色视频高清免费观看| 国产精品久久久久久久久免| 男女边吃奶边做爰视频| 最近中文字幕2019免费版| 免费看av在线观看网站| 九色成人免费人妻av| 久久女婷五月综合色啪小说| 天堂俺去俺来也www色官网| 97超视频在线观看视频| 黄色配什么色好看| 在线免费观看不下载黄p国产| 最近最新中文字幕免费大全7| 精品人妻熟女毛片av久久网站| 欧美日韩av久久| 亚洲av男天堂| 久久久久久久久久成人| 下体分泌物呈黄色| 亚洲国产av新网站| 国产精品99久久久久久久久| 久久久久精品性色| 亚洲综合色惰| 国产一区二区三区综合在线观看 | 日韩在线高清观看一区二区三区| 国产亚洲精品久久久com| 少妇丰满av| 欧美老熟妇乱子伦牲交| 五月天丁香电影| 免费播放大片免费观看视频在线观看| 少妇被粗大的猛进出69影院 | 国产69精品久久久久777片| 国产国拍精品亚洲av在线观看| 欧美最新免费一区二区三区| 婷婷色av中文字幕| 欧美少妇被猛烈插入视频| 欧美日韩一区二区视频在线观看视频在线| 如何舔出高潮| 中文字幕最新亚洲高清| 欧美三级亚洲精品| 高清毛片免费看| 人人妻人人添人人爽欧美一区卜| 国产女主播在线喷水免费视频网站| 高清在线视频一区二区三区| 国产黄片视频在线免费观看|