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

    確定直購電交易分?jǐn)傠娏W(wǎng)絡(luò)固定成本的新方法

    2015-07-18 11:06:35周盈文福拴朱炳銓項中明徐立中趙良浙江大學(xué)電氣工程學(xué)院杭州0027文萊科技大學(xué)電機(jī)與電子工程系斯里巴加灣BE40文萊國網(wǎng)浙江省電力公司杭州0027
    關(guān)鍵詞:分?jǐn)?/a>約束權(quán)重

    周盈,文福拴,2,朱炳銓,項中明,徐立中,趙良(.浙江大學(xué)電氣工程學(xué)院,杭州0027;2.文萊科技大學(xué)電機(jī)與電子工程系,斯里巴加灣BE40,文萊;.國網(wǎng)浙江省電力公司,杭州0027)

    確定直購電交易分?jǐn)傠娏W(wǎng)絡(luò)固定成本的新方法

    周盈1,文福拴1,2,朱炳銓3,項中明3,徐立中3,趙良3
    (1.浙江大學(xué)電氣工程學(xué)院,杭州310027;2.文萊科技大學(xué)電機(jī)與電子工程系,斯里巴加灣BE1410,文萊;3.國網(wǎng)浙江省電力公司,杭州310027)

    與大用戶直購電相關(guān)的如何確定其應(yīng)該承擔(dān)的電力網(wǎng)絡(luò)固定成本這一個重要問題,仍未得到很好解決。在此背景下,對大用戶直購電交易所需分?jǐn)偟碾娏W(wǎng)絡(luò)固定成本進(jìn)行了系統(tǒng)分析,并考慮了交易的多方面屬性。首先從用戶側(cè)、發(fā)電側(cè)和電網(wǎng)環(huán)節(jié)3個方面構(gòu)建了與直購電交易相關(guān)的電力網(wǎng)絡(luò)固定成本分?jǐn)偟闹笜?biāo)體系;之后,在數(shù)據(jù)包絡(luò)分析框架下建立了聯(lián)盟博弈模型,并利用區(qū)間層次分析法在該模型中引入了反映決策者偏好的權(quán)重約束,進(jìn)而獲得了一種能夠計及指標(biāo)權(quán)重約束的電力網(wǎng)絡(luò)固定成本分?jǐn)偰P?;最后,以IEEE14節(jié)點系統(tǒng)為例說明了所提方法的基本特征。算例結(jié)果表明,這種方法能夠在滿足決策者偏好的前提下保證分?jǐn)偨Y(jié)果的公平合理。

    直購電;固定成本分?jǐn)?;指?biāo)體系;數(shù)據(jù)包絡(luò)分析;聯(lián)盟博弈;權(quán)重約束

    大用戶直接交易是電力市場改革的重要組成部分。到目前為止,已有吉林、廣東、福建、安徽和遼寧等省份開展了大用戶直接交易試點工作,內(nèi)蒙古電力多邊交易市場中也包含大用戶直接交易[1]。然而,與大用戶直購電相關(guān)的輸配電價機(jī)制尚未形成,成為相關(guān)各方目前爭論的焦點之一,這已經(jīng)在相當(dāng)程度上影響了這項工作的推廣。在電力網(wǎng)絡(luò)成本中,固定成本占相當(dāng)大的比重,研究如何將電力網(wǎng)絡(luò)固定成本在各交易之間進(jìn)行分?jǐn)偩哂泻軓?qiáng)的應(yīng)用背景。

    在電力網(wǎng)絡(luò)固定成本分?jǐn)偡矫?,國?nèi)外已有相當(dāng)多的研究報道,所提出的方法從總體上可以分為綜合成本法、邊際成本法和基于博弈論的分?jǐn)偡椒ǖ取>C合成本法是目前使用最為廣泛的分?jǐn)偡椒?,其中郵票法、合同路徑法[2]等分?jǐn)偡椒ㄔ肀容^簡單,易于實現(xiàn),但其沒有考慮實際潮流狀況和各交易對電力網(wǎng)絡(luò)的使用情況和對系統(tǒng)產(chǎn)生的影響,因此分?jǐn)偨Y(jié)果難以公平合理地體現(xiàn)基于使用的分?jǐn)傇瓌t。潮流追蹤法、發(fā)電轉(zhuǎn)移分布因子法、廣義發(fā)電分布因子法和廣義負(fù)荷分布因子法等綜合成本分?jǐn)偡ㄊ峭ㄟ^將線路上的總潮流分解為各個“源”或“流”所引起的潮流分量之和,從而評估用戶(或交易)對電力網(wǎng)絡(luò)的使用程度,此類方法又稱為源流分析法[3-4]。實際上很難確定每個負(fù)荷所用電力來自哪個發(fā)電機(jī),所以這類方法均基于某種主觀假設(shè)。例如,潮流追蹤法假設(shè)電力在任一節(jié)點處的分配滿足比例共享原則。邊際成本法[5]分為短期邊際成本法和長期邊際成本法,短期邊際成本法不能有效回收固定成本,而長期邊際成本法需要對長期成本和負(fù)荷進(jìn)行預(yù)測,準(zhǔn)確性難以保證,計算也比較復(fù)雜。將博弈理論應(yīng)用于電力網(wǎng)絡(luò)成本分?jǐn)傋钤缬扇毡緦W(xué)者Tsukamoto等[6]提出。這類方法考慮了不同交易之間客觀存在的相互博弈關(guān)系,最終達(dá)成一個各方都愿意接受的分?jǐn)偡桨?。到目前為止,這類方法主要基于合作博弈,將參與電力交易的各方或電力市場成員視為局中人,將潮流節(jié)省量折算到線路固定成本或?qū)⒊杀竟?jié)省量視為特征函數(shù)[7]。

    上述固定成本分?jǐn)偡椒ㄋ紤]的因素一般側(cè)重于某個或某些方面,如依據(jù)潮流情況或用電量分?jǐn)偂kS著電力市場的深化和發(fā)展,需要考慮的因素會更多,例如交易發(fā)生的時間、交易電能的質(zhì)量等均會影響系統(tǒng)成本,根據(jù)“誰使用,誰承擔(dān)”的原則這些因素也需要適當(dāng)考慮。一種好的固定成本分?jǐn)偡椒☉?yīng)該能夠綜合考慮交易各方面的特性,并力求分?jǐn)偨Y(jié)果公平。

    數(shù)據(jù)包絡(luò)分析DEA(dataenvelopmentanalysis)是一種評價多屬性決策單元之間相對效率的方法。文獻(xiàn)[8]首次嘗試把DEA應(yīng)用于成本分?jǐn)倖栴},計及了分?jǐn)偳昂笏袇⑴c的成員效率不變性和帕雷托最小性。之后,許多學(xué)者開始將DEA與博弈理論方法相結(jié)合來解決公共平臺系統(tǒng)建設(shè)費用的分?jǐn)倖栴}[9]。近年來,有些學(xué)者提出了在數(shù)據(jù)包絡(luò)分析中體現(xiàn)決策者指標(biāo)或決策單元偏好的思路與方法。文獻(xiàn)[10]將反映決策者意愿的“偏好錐”引入DEA理論,建立了與帶產(chǎn)出錐和投入錐相應(yīng)的兩階段生產(chǎn)可能集;文獻(xiàn)[11]分別考慮了對各權(quán)重之比進(jìn)行約束和對虛擬輸出變量(輸出變量與相應(yīng)權(quán)重的乘積)之比進(jìn)行約束這兩種形式;文獻(xiàn)[12]將輸入輸出視為模糊數(shù),通過引入Campos平均指數(shù),將模糊集合之間的關(guān)系轉(zhuǎn)化為一般線性關(guān)系,建立了基于Campos指數(shù)的模糊DEA;文獻(xiàn)[13]采用層次分析法AHP(analytic hierarchy process)構(gòu)建權(quán)重約束,并在此基礎(chǔ)上運用DEA對電網(wǎng)規(guī)劃方案進(jìn)行了評價。上述方法有的過于復(fù)雜,限制了其應(yīng)用范圍(如文獻(xiàn)[10-12]);有的則過于依賴決策者的主觀判斷(如文獻(xiàn)[13])。區(qū)間層次分析法IAHP(interval analytic hierarchy process)[14]由AHP發(fā)展而來,是一種實用的多維決策方法,其將AHP與區(qū)間數(shù)相結(jié)合,更加適合解決主觀判斷和決策屬性的不確定性問題。

    本文將基于IAHP-DEA的聯(lián)盟博弈方法應(yīng)用于解決電力網(wǎng)絡(luò)固定成本分?jǐn)倖栴}。首先建立了較為全面的電力網(wǎng)絡(luò)固定成本分?jǐn)傊笜?biāo)體系,其中不僅考慮了電壓等級、電能質(zhì)量、交易電量,還考慮了發(fā)電側(cè)供電煤耗和排污水平、交易發(fā)生時間、地點等因素對固定成本分?jǐn)偟挠绊憽V?,利用基于DEA的聯(lián)盟博弈模型得到各交易的分?jǐn)偨Y(jié)果,借助了IAHP將決策者偏好信息轉(zhuǎn)換成約束條件加入到傳統(tǒng)DEA模型中,即對固定成本分?jǐn)傊笜?biāo)體系中各指標(biāo)權(quán)重進(jìn)行約束,并對無考慮權(quán)重約束時的分?jǐn)偨Y(jié)果進(jìn)行了比較。

    1 電力網(wǎng)絡(luò)固定成本分?jǐn)偟脑u估指標(biāo)體系

    建立科學(xué)而系統(tǒng)的指標(biāo)體系對于與大用戶直接交易相關(guān)的電力網(wǎng)絡(luò)固定成本分?jǐn)偩哂兄匾饬x。到目前為止,專門研究與電力網(wǎng)絡(luò)固定成本分?jǐn)傁嚓P(guān)的指標(biāo)體系的文獻(xiàn)很少。隨著電力交易的種類逐步增多,文獻(xiàn)[15]采用分時電價和電能質(zhì)量指標(biāo)表征交易負(fù)荷的峰谷特性和電能質(zhì)量要求對電力網(wǎng)絡(luò)成本的影響,在分?jǐn)傔^程中綜合考慮了分時電價、電能質(zhì)量和用電量3個指標(biāo)。在上述背景下,本文考慮了參與直購電的各利益相關(guān)方面,從用電側(cè)、發(fā)電側(cè)和對電力網(wǎng)絡(luò)影響3個方面出發(fā),建立了與大用戶直購電相關(guān)的電力網(wǎng)絡(luò)固定成本分?jǐn)傊笜?biāo)體系,如圖1所示。

    在圖1中,采用了交易電量、電壓等級和電能質(zhì)量這3個指標(biāo)描述電能交易在用戶側(cè)的屬性;引入供電煤耗和污染物排放水平指標(biāo)來描述電能交易在發(fā)電側(cè)的節(jié)能環(huán)保特性,以期在公平分?jǐn)偣潭ǔ杀镜耐瑫r合理引導(dǎo)大用戶的購電與用電行為;引入峰谷電價和線路負(fù)荷率不均衡度變化這2個指標(biāo)來表征交易發(fā)生時間和交易雙方所處位置對電力網(wǎng)絡(luò)成本的影響。

    圖1 電力網(wǎng)絡(luò)固定成本分?jǐn)傊笜?biāo)體系Tab.1 Index systeMofallocating power network fixed cost

    2 各指標(biāo)的量化與處理

    2.1 用電側(cè)指標(biāo)

    1)交易電量指標(biāo)Y1

    大用戶直購電的年交易電量越大,該筆交易應(yīng)分?jǐn)偟墓潭ǔ杀揪驮酱?。顯然,若只考慮交易電量這一指標(biāo),就蛻化為郵票法了。Y1的單位為萬kW·h。

    2)電壓等級指標(biāo)Y2

    用戶所在電壓等級不同,電力網(wǎng)絡(luò)送電所經(jīng)過的設(shè)備數(shù)量和距離可能不同,導(dǎo)致的網(wǎng)絡(luò)損耗一般也不同,則相關(guān)的電力網(wǎng)絡(luò)供電成本也不同。我國電力用戶的電壓等級一般分為不滿1 kV、1~10 kV、35~110 kV、110 kV和220 kV及以上5類,把Y2的量化值依次記為1、2、3、4和5。

    3)電能質(zhì)量指標(biāo)Y3

    在電力市場發(fā)展比較成熟的國家和地區(qū),電網(wǎng)企業(yè)可以根據(jù)用戶需要提供不同水平的電能質(zhì)量服務(wù),電能質(zhì)量要求高的用戶需要承擔(dān)較高或額外的輸電費用,這體現(xiàn)了公平負(fù)擔(dān)原則。利用電能質(zhì)量指標(biāo)Y3對電能質(zhì)量國家標(biāo)準(zhǔn)中所涉及的4個指標(biāo),即供電電壓允許偏差、公用電網(wǎng)諧波、電壓允許波動、三相電壓允許不平衡度,進(jìn)行綜合評估,具體計算公式為式中:SVD為電壓偏差指標(biāo);SHD為諧波畸變率指標(biāo);SVF為電壓波動指標(biāo);SVU為三相不平衡度指標(biāo)。各指標(biāo)計算方法參照文獻(xiàn)[15]。Y3的值越小,表示用戶對電能質(zhì)量的要求越高,需要分?jǐn)偟木W(wǎng)絡(luò)固定成本也就越大。

    2.2 發(fā)電側(cè)指標(biāo)

    考慮到不同發(fā)電廠生產(chǎn)單位電能所產(chǎn)生的污染排放量不同,大用戶選擇不同的發(fā)電廠進(jìn)行電力交易時產(chǎn)生的環(huán)境效益也不同。政府相關(guān)部門作為分?jǐn)倶?biāo)準(zhǔn)的制定者,為促進(jìn)節(jié)能減排,有必要引入相關(guān)指標(biāo)。

    發(fā)電側(cè)指標(biāo)主要考察發(fā)電廠的節(jié)能減排情況。從社會責(zé)任的角度考慮,某筆交易的發(fā)電側(cè)的供電煤耗及污染物排放水平越低,該筆交易分?jǐn)偟妮旊姽潭ǔ杀究梢栽降?,即通過價格引導(dǎo)大用戶選擇節(jié)能減排情況較好的機(jī)組,從而激勵電廠降低能耗及排放水平。

    1)供電煤耗指標(biāo)Y4

    參照文獻(xiàn)[16]中的機(jī)組性能考核系數(shù),將實際供電煤耗與機(jī)組優(yōu)化基準(zhǔn)供電煤耗之差占機(jī)組優(yōu)化基準(zhǔn)供電煤耗的比重作為供電煤耗指標(biāo),以避免機(jī)組基準(zhǔn)供電煤耗的差異對指標(biāo)的影響。具體計算模型為

    式中:Bh為機(jī)組運行實際供電煤耗,g/(kW·h),即總發(fā)電量除去廠用電量后向系統(tǒng)供電的實際煤耗;Bj-h為機(jī)組優(yōu)化基準(zhǔn)供電煤耗,g/(kW·h),Bj-h是指機(jī)組實際運行于某一工況時所能達(dá)到的最低煤耗,可借助機(jī)組性能在線監(jiān)測,通過改變機(jī)組的運行參數(shù),進(jìn)行大量優(yōu)化實驗及耗差分析確定。

    2)污染物排放指標(biāo)Y5

    污染物排放指標(biāo)包括3個子指標(biāo):煙塵排放指標(biāo)Y51、二氧化硫排放指標(biāo)Y52和氮氧化合物排放指標(biāo)Y53。則污染物排放指標(biāo)可表示為

    這3個子指標(biāo)的計算可參見文獻(xiàn)[17]。基本思路是按污染物排放限值將污染物的排放量分不同的等級,并且不同等級對應(yīng)不同的分值,分值都取0~1之間的數(shù);分值越高說明對應(yīng)污染物的排放量越高。

    2.3 電網(wǎng)影響指標(biāo)

    1)峰谷電價Y6

    電力負(fù)荷一般具有明顯的峰谷特性,同樣一筆交易發(fā)生在負(fù)荷高峰時期和低谷時期對電網(wǎng)運行的影響是不同的。為了應(yīng)對負(fù)荷高峰,電網(wǎng)公司不得不投入資金來增加輸電容量以保證電力供應(yīng)和系統(tǒng)運行的安全性。在一些國家或地區(qū),對用戶側(cè)的銷售電價實行峰谷電價也是出于這個考慮。因此,本文將峰谷電價指標(biāo)作為衡量交易發(fā)生時間對電力網(wǎng)絡(luò)影響的指標(biāo),其具體計算公式[15]為

    式中:n為參與直購電的大用戶數(shù)目;Ti為大用戶i在典型日的用電小時數(shù),h;α∶1∶β為典型日的峰、平、谷時段的電價比;Tif、Tip和Tig分別為大用戶i在典型日的峰、平、谷時段用電小時數(shù),h。顯然,谷時段利用小時數(shù)占全天開工小時數(shù)的比例越高,峰時段占的比例越小,該大用戶分?jǐn)偟降墓潭ǔ杀疽簿驮缴?,這有利于引導(dǎo)大用戶避開用電高峰期,以緩解電網(wǎng)在高峰時段的輸電壓力。

    2)線路負(fù)載率不均衡度變化指標(biāo)Y7

    線路負(fù)載率是指線路中實際流過的有功潮流與線路可輸送的極限有功潮流之比,計算公式為

    式中:fi為第i條線路的負(fù)載率;Pi為第i條線路上實際輸送的有功功率;Pimax為第i條線路可輸送的有功功率最大值。輸電線路的負(fù)載率越低,則相應(yīng)的線路容量裕度越大,系統(tǒng)調(diào)度的靈活性越大。

    對于實際電力系統(tǒng)而言,某一條或幾條輸電線路的負(fù)載率不能全面反映整個電力網(wǎng)絡(luò)的整體運行狀態(tài)。平均負(fù)載率表示系統(tǒng)中所有線路負(fù)載率的平均值,雖然能從整體上反映出系統(tǒng)的負(fù)荷水平,但不能反映各條輸電線路負(fù)載率和平均負(fù)載率之間的差,有可能出現(xiàn)某些遠(yuǎn)遠(yuǎn)大于或小于平均負(fù)載率的線路。本文采用文獻(xiàn)[18]提出的基于標(biāo)準(zhǔn)差的負(fù)載率不均衡度來衡量電力網(wǎng)絡(luò)的運行狀態(tài),即

    式中:E為線路負(fù)載率不均衡度;b為線路數(shù);fˉ為線路的加權(quán)平均負(fù)載率。實際上,E就是線路負(fù)載率的標(biāo)準(zhǔn)差。當(dāng)所有線路的負(fù)載率集中在同一水平線上下時,負(fù)載率不均衡度就會較小。當(dāng)有的線路負(fù)載率很高,而有的很低時,負(fù)載率不均衡度就會比較大。

    這里將交易j實施前后,線路不均衡度的變化作為固定成本分?jǐn)偟闹笜?biāo)之一,即

    式中:Eall和E-j分別為包含所有交易時和不包含第j筆交易時的電力網(wǎng)絡(luò)線路負(fù)載率不均衡度;t為交易數(shù)目??紤]第j筆交易后,線路負(fù)載率不均衡度增加的越多,該交易應(yīng)該分?jǐn)偟墓潭ǔ杀揪驮蕉啵蝗艨紤]該筆交易后線路負(fù)載率不均衡度下降,則該交易對系統(tǒng)的安全和經(jīng)濟(jì)運行是有利的,此時Y7取0。Y7的值與交易雙方在電力網(wǎng)絡(luò)中所處位置有關(guān),用Y7表征交易發(fā)生地點對成本分?jǐn)偟挠绊憽?/p>

    2.4 指標(biāo)類型的一致化處理

    在上述指標(biāo)中,有些指標(biāo)值與分?jǐn)偟墓潭ǔ杀臼钦嚓P(guān)的,即指標(biāo)值越小,應(yīng)分?jǐn)偟墓潭ǔ杀疽苍缴伲环粗畡t稱為負(fù)相關(guān)指標(biāo)。顯然,除了電壓等級和電能質(zhì)量指標(biāo)之外,其余均屬于正相關(guān)指標(biāo)。在進(jìn)行綜合評價之前,需對評價指標(biāo)類型進(jìn)行一致化處理。對于負(fù)相關(guān)指標(biāo),令

    式中:Y為負(fù)相關(guān)指標(biāo)的指標(biāo)值;Y*為經(jīng)一致化處理后的指標(biāo)值。

    2.5 指標(biāo)的無量綱化

    為盡可能地反映實際情況,避免由于各項指標(biāo)的量綱不同以及其數(shù)量級間的懸殊差別所帶來的負(fù)面影響,需要對指標(biāo)進(jìn)行無量綱化處理。這里采用歸一化處理方法,即

    式中,yrj和yrj′分別為交易j的第r個指標(biāo)歸一化處理前與處理后的值。

    3 基于IAHP-DEA的聯(lián)盟博弈方法

    3.1 滿足DEA有效的分?jǐn)偡桨讣?/p>

    將DEA引入到固定成本分?jǐn)倖栴}時,每筆交易就相對于一個決策單元DMU(decision making unit),而前文提出的指標(biāo)則相當(dāng)于系統(tǒng)的輸出指標(biāo),每筆交易分?jǐn)偟降墓潭ǔ杀颈闶瞧湮ㄒ坏妮斎胫笜?biāo)。這樣,第k個DMU的DEA效率可由DEA中的Charnes-Cooper-Rhodes(CCR)優(yōu)化模型[19]求取。具體模型可描述為

    式中:hk和yrk分別為第k個DMU的DEA效率和其第r個輸出指標(biāo)值;ur為第r個輸出指標(biāo)的權(quán)重;Rk和Rj分別為第k個和第j個DMU所分?jǐn)偟墓潭ǔ杀颈壤?;s為輸出指標(biāo)個數(shù)。式(12)~式(15)為約束條件。式(12)表明DEA效率并非絕對效率,而是衡量該DMU的生產(chǎn)能力與其對應(yīng)的前沿生產(chǎn)能力之間差距的相對效率,故效率值必然在[0,1]之間;式(13)表明各DMU承擔(dān)的固定成本之和應(yīng)等于待分?jǐn)偟目偣潭ǔ杀?;式?4)和式(15)均為模型中變量的非負(fù)約束。

    同理,DMU整體的DEA效率優(yōu)化模型可描述為

    式中:hall為DMU整體的DEA效率;hj為第j個DMU的DEA效率;約束條件同式(12)~式(15)。

    可以證明[19],所有DMU和整體均為DEA有效,即DEA效率均能達(dá)到1。綜合以上2個模型,同時滿足DMU個體和整體為DEA有效的分?jǐn)偡桨讣峡梢员硎緸?/p>

    3.2 DEA聯(lián)盟博弈模型

    假設(shè)S為決策集合N={1,2,…,t}的任一子集構(gòu)成的聯(lián)盟,定義S的各指標(biāo)值為

    式(20)表示在使得自身聯(lián)盟的分?jǐn)傋钚』耐瑫r,還要保證分?jǐn)偨Y(jié)果滿足DMU個體和整體為DEA有效。易得,V(φ)=0,V(N)=1,且具次加性,即V(S∪W)≤V(S)+V(W),因此可將V(S)作為聯(lián)盟S的特征函數(shù)。

    3.3 聯(lián)盟博弈的解

    多人合作博弈有多種不同的解的表達(dá)形式,常用的包括核心(core)、核仁(nucleolus)和Shapley值。其中只有核仁解必定存在且唯一[20]。為此,這里采用核仁作為聯(lián)盟博弈的解。

    假設(shè)R=(R1,R2,…,Rt)是t筆大用戶直接交易合作博弈的一個分?jǐn)偡桨福敲幢厝挥?/p>

    式中,|S|為S中包含的交易個數(shù)。該式說明了對于任何S,其成員單獨行動所付出的成本之和必然不小于形成聯(lián)盟后的總成本,否則就沒有成員愿意加入該聯(lián)盟。根據(jù)式(21)定義S對分?jǐn)偡桨窻的不滿意度為

    核仁算法的基本原理[20]是:優(yōu)先考慮最不滿意的聯(lián)盟,選擇分配使這種聯(lián)盟的不滿意度達(dá)到最??;再考慮次不滿意的聯(lián)盟,并選擇分配使其不滿意程度盡可能小。如此循環(huán),直到所有聯(lián)盟的不滿意度都盡可能小為止。因此,基于核仁解的固定成本分?jǐn)偰P蜑?/p>

    在最后的分?jǐn)偡桨钢懈鱾€DMU使用相同的權(quán)重向量u,表示了各指標(biāo)對分?jǐn)偨Y(jié)果的影響,其能夠給市場參與者提供一定的經(jīng)濟(jì)信號。

    3.4 權(quán)重約束的確定

    3.4.1 引入權(quán)重約束的必要性

    在各項指標(biāo)值不變的情況下,不同的權(quán)重值會導(dǎo)致不同的分?jǐn)偳闆r。各個DMU在博弈中勢必選擇對自身最有利的權(quán)重方案。因此,從某種意義上講,權(quán)重的靈活性為各成員提供了相互博弈的空間。但過于靈活的權(quán)重也是此方法的軟肋所在。從式(19)和式(20)可以看出,如果沒有權(quán)重約束,則最優(yōu)解為最有利于該DMU的指標(biāo)的權(quán)重為1,其余為0。然而,這種情況是不現(xiàn)實的,因為每個DMU在選擇各自的最優(yōu)權(quán)重時,不僅要使自己的分?jǐn)傤~盡可能小,還要考慮以下因素:

    (1)上層決策者對于各個指標(biāo)的偏好;

    (2)指標(biāo)體系的有效性。例如,各指標(biāo)的權(quán)重不能太高;所有指標(biāo)的權(quán)重不能為0。這是為了防止各DMU過度強(qiáng)調(diào)自己的強(qiáng)項,從而造成權(quán)重分布過于懸殊。

    鑒于此,本文采用上下限約束的形式對權(quán)重進(jìn)行了約束,即為每個指標(biāo)設(shè)置不同的權(quán)重變化范圍。

    3.4.2 用區(qū)間層次分析法求解權(quán)重約束

    考慮到人們對一些事物的比較判斷有時難以用準(zhǔn)確的數(shù)值來表達(dá),采用區(qū)間表示更為合適。因此,這里采用IAHP求解指標(biāo)的權(quán)重約束。IAHP用區(qū)間數(shù)代替AHP所用的點值來構(gòu)造判斷矩陣,并通過區(qū)間運算得到權(quán)重向量,原始數(shù)據(jù)和計算結(jié)果都以區(qū)間數(shù)形式表達(dá),從而可得到更具“柔性”的結(jié)果。這樣既能保證指標(biāo)體系的合理性,又能體現(xiàn)指標(biāo)之間的相對重要性。

    1)構(gòu)造判斷矩陣

    與AHP一樣,IAHP也要首先建立各指標(biāo)兩兩比較的判斷矩陣。但對于IAHP,其判斷矩陣中的元素不再是確定數(shù),而是區(qū)間數(shù),表示決策者對于某指標(biāo)與另一指標(biāo)之間相對重要性的估計范圍。可用如表1所示的1~9標(biāo)度方法將決策判斷定量化。如果認(rèn)為因素i與因素j相比介于“明顯”、“強(qiáng)烈”到“極端”之間,可以將結(jié)果表示為區(qū)間數(shù)[6,8]。根據(jù)互反性,因素j與因素i的比較結(jié)果為[1/8,1/6]。對角線元素aii表示元素i與自身比較的結(jié)果,顯然為“同等”,即[1,1]。

    表1 1~9標(biāo)度系統(tǒng)Tab.1 1~9 scalesystem

    2)一致性校驗

    設(shè)A=(aij)m×m為區(qū)間判斷矩陣,其中的元素為aij=[,。記=(a)m×m=(和分別為矩陣A的上邊界和下邊界點值矩陣。參照AHP的矩陣一致性檢驗方法,定義區(qū)間矩陣一致性的指標(biāo)為

    式中:CR為一致性比率,一般認(rèn)為當(dāng)CR≤0.1時,A的一致性程度滿足要求[21];CI為一致性指標(biāo),m為判斷矩陣的階數(shù),λ-和λ+分別為A-和A+的最大特征值;CS為隨機(jī)一致性指標(biāo),取值與m有關(guān)。

    3)求解判斷矩陣權(quán)重

    由于判斷矩陣元素為區(qū)間數(shù),就必須采用區(qū)間數(shù)判斷矩陣的權(quán)重求解算法。目前常用的算法包括區(qū)間特征根法IEM(intervaleigenvaluemethod)、隨機(jī)模擬法、迭代法等[21-22]。其中,隨機(jī)模擬法計算量大,迭代法誤差較大,而IEM在精度和計算量方面的綜合性能較好,本文采用IEM求解權(quán)重。

    假設(shè)x-、x+分別為A-、A+對應(yīng)于最大特征根λ-、λ+的具有正分量的歸一化特征向量。根據(jù)區(qū)間數(shù)及特征根運算基本定理[21],區(qū)間數(shù)權(quán)重向量u為

    式中,u-和u+為權(quán)重的下限約束和上限約束。在原始模型中加入權(quán)重約束u-≤u≤u+,即得到基于IAHP-DEA的聯(lián)盟博弈模型。

    4 算例

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

    以圖2所示IEEE 14節(jié)點系統(tǒng)[23]為例來說明所提出的方法,各支路參數(shù)參見文獻(xiàn)[23]。

    圖2 14節(jié)點系統(tǒng)拓?fù)浣Y(jié)構(gòu)Fig.2 14-bussysteMtopology structure

    假設(shè)有3個獨立發(fā)電機(jī)組G1、G2和G3和4個大用戶參與直接交易。具體交易信息見表2。

    表2 系統(tǒng)中大用戶直接交易數(shù)據(jù)Tab.2 Direct trade data in system

    表3和表4分別給出了參與直購電的大用戶及發(fā)電機(jī)組的具體信息。為簡化起見,這里直接給出了大用戶電能質(zhì)量的指標(biāo)值和污染物排放子指標(biāo)的分值。

    表3 系統(tǒng)中大用戶的基本情況Tab.3 Basic information of large consumers in the system

    給定劃分的峰谷電價時段:①峰時:10∶00—12∶00,18∶00—22∶00;②平時:8∶00—10∶00,12∶00—18∶00,22∶00—24∶00;③谷時:0∶00—8∶00。峰、平、谷電價之比為:2∶1∶0.5。由此可以計算出4個大用戶的分時電價指標(biāo)Y6=[1.083 3,0.500 0,1.454 5,1.454 5]。

    表4 系統(tǒng)中發(fā)電機(jī)組的基本情況Tab.4 Basic information ofgenerators in the system

    4.2 算法流程

    步驟1根據(jù)第2.3節(jié)的計算方法以及各支路的電抗值,計算出線路負(fù)載率不均衡度變化指標(biāo):Y7=[0.003 6,0.010 4,0.001 7,0.012 5,0.053 5]。然后對各交易的指標(biāo)數(shù)據(jù)作一致化及歸一化處理,可得到表5。

    表5 各交易的歸一化指標(biāo)Tab.5 Normalized indexmatrix

    步驟2采用IAHP方法計算各指標(biāo)的權(quán)重約束。首先形成如表6所示的區(qū)間判斷矩陣;然后,求取其一致性指標(biāo)CR=0.052 0,其小于0.1,因此滿足一致性要求;最后,根據(jù)IEM得到權(quán)重的上下限約束分別為:u+=[0.296 1;0.152 7;0.095 7;0.059 7;0.048 8;0.234 5;0.234 5],u-=[0.247 1; 0.100 8;0.062 7;0.049 6;0.038 2;0.189 6;0.189 6]。

    步驟3利用式(19)和式(20)分別計算未考慮權(quán)重約束和考慮權(quán)重約束時,各聯(lián)盟的特征函數(shù)及權(quán)重分配。

    表6 電力網(wǎng)絡(luò)固定成本分?jǐn)傊笜?biāo)區(qū)間判斷矩陣Tab.6 Interval judgmentmatrix for power network fixed costallocation

    步驟4根據(jù)式(23)~式(25)分別求取無權(quán)重約束和有權(quán)重約束條件下聯(lián)盟博弈的核仁解,即固定成本分?jǐn)偡桨浮?/p>

    4.3 分?jǐn)偨Y(jié)果及分析

    表7和表8分別顯示了無權(quán)重約束和有權(quán)重約束時固定成本分?jǐn)偡桨讣跋鄳?yīng)權(quán)重。圖3顯示了郵票法與基于DEA的聯(lián)盟博弈方法的對比。

    表7 聯(lián)盟博弈的核仁解Tab.7 Nucleolus solution of the cooperative game

    圖3 3種方法得到的固定成本分?jǐn)偙壤龑Ρ菷ig.3 Comparisonsofallocation proportionsof fixed costw ith threemethods

    由圖3可見,基于DEA的聯(lián)盟博弈算法由于綜合考慮了電壓、電能質(zhì)量、交易電量、線路負(fù)荷率不均衡度以及機(jī)組的節(jié)能減排等因素,得到的分?jǐn)偨Y(jié)果比較合理。對比交易1和交易3,在郵票法下由于交易3的交易電量遠(yuǎn)大于交易1,因此分?jǐn)偟墓潭ǔ杀疽脖冉灰?大。實際上,雖然交易1的交易電量小且電能質(zhì)量要求低,但其發(fā)電側(cè)節(jié)能減排效果不如交易3,且其峰谷電價指標(biāo)和線路負(fù)載率不均衡度變化指標(biāo)都比交易3大,這說明交易1對電力網(wǎng)絡(luò)的影響相對較大。采用基于DEA的聯(lián)盟博弈方法時交易1和交易3分?jǐn)偟降墓潭ǔ杀鞠嘟_@體現(xiàn)了相互博弈形成的結(jié)果比較公平合理的特征。

    表8 核仁解中各指標(biāo)的公共權(quán)重Tab.8 Common weightsofnucleolussolution

    從表8可知,在不考慮權(quán)重約束時,各交易間博弈的結(jié)果是:線路負(fù)載率不均衡度變化指標(biāo)的權(quán)重最大,沒有考慮電能質(zhì)量指標(biāo)、機(jī)組供電煤耗指標(biāo)及峰谷電價指標(biāo)。這樣的權(quán)重分配不夠合理,也不能給市場參與者提供正確的經(jīng)濟(jì)信號。引入權(quán)重約束一方面可保證所有指標(biāo)的權(quán)重差異不會過大,另一方面可使得各指標(biāo)的重要性符合決策者的主觀偏好。

    5 結(jié)語

    電力網(wǎng)絡(luò)固定成本分?jǐn)倢ν苿哟笥脩糁苯咏灰字陵P(guān)重要。本文采用基于IAHP-DEA的聯(lián)盟博弈模型來求解這一問題,并通過求取核仁解獲得分?jǐn)偙壤?。由算例計算結(jié)果可得到如下結(jié)論。

    (1)采用所提方法得到的分?jǐn)偨Y(jié)果不僅考慮了交易電量,還考慮了用戶所在電壓等級、要求的電能質(zhì)量,發(fā)電側(cè)機(jī)組的節(jié)能減排情況以及不同交易的發(fā)生時段、發(fā)生位置對固定成本分?jǐn)偟挠绊?,考慮的因素比較全面,得到的分?jǐn)偨Y(jié)果更為合理。

    (2)通過DEA與聯(lián)盟博弈相結(jié)合,實現(xiàn)了在滿足整體及各筆交易效率最大的前提下,各方對分?jǐn)偨Y(jié)果的不滿意度最小,這兼顧了效率與公平。

    (3)通過采用IAHP,引入了指標(biāo)權(quán)重的上下限約束,這不僅可以在相當(dāng)程度上避免過分依賴部分有利于自身的指標(biāo)的缺陷,而且能夠使得權(quán)重選取反映上層決策者的意愿,這體現(xiàn)了主客觀相結(jié)合的思路。

    (4)采用核仁解能夠得出一個確定的指標(biāo)權(quán)重,其表征了各指標(biāo)對于分?jǐn)偨Y(jié)果的重要性。這樣可使固定成本的分?jǐn)偨Y(jié)果透明化,并能夠給市場主體提供一定的經(jīng)濟(jì)信號。

    本文只是提出了電力網(wǎng)絡(luò)固定成本分?jǐn)偟囊环N新的思路,如要應(yīng)用于實際電力交易則還有一些具體問題需要研究。例如,如何協(xié)調(diào)跨省區(qū)大用戶直接交易、如何處理相關(guān)的輔助服務(wù)成本等問題。

    [1]夏清,白楊,鐘海旺,等(Xia Qing,BaiYang,Zhong Haiwang,etal).中國推廣大用戶直購電交易的制度設(shè)計與建議(Institutionaldesign and suggestions for promotion ofdirectelectricity purchase by large consumers in China)[J].電力系統(tǒng)自動化(Automation of Electric Power Systems),2013,37(20):1-7.

    [2]Kristiansen T.Comparison oftransmission pricingmodels[J]. International Journal of Electrical Power&Energy Systems,2011,33(4):947-953.

    [3]李傳棟,王建明,楊桂鐘(LiChuandong,Wang Jianming,Yang Guizhong).潮流追蹤解析算法(Analytical algorithMfor tracing power flow)[J].電力系統(tǒng)及其自動化學(xué)報(Proceedings of the CSU-EPSA),2013,25(3):119-123,146.

    [4]荊朝霞,廖名洋,聶更生,等(Jing Zhaoxia,Liao Mingyang,Nie Gengsheng,etal).輸電固定成本分?jǐn)傊凶畲蠡C合效益的源流分析方法(A source-sink analysis method formaximizing socialsurplus in transmission fixed costallocation)[J].電力系統(tǒng)自動化(Automation of Electric Power Systems),2013,37(23):71-77.

    [5]傅書逷,張銳,王海寧(Fu Shuti,Zhang Rui,Wang Haining).輸電費的計算與分?jǐn)偅–omputation and allocation of transmission charges)[J].電力系統(tǒng)自動化(Automation of Electric Power Systems),2004,28(3):5-9.

    [6]Tsukamoto Y,Iyoda I.Allocation of fixed transmission cost to wheeling transactions by cooperative game theory[J]. IEEETranson Power Systems,1996,11(2):620-629.

    [7]龔演平(Gong Yanping).電力市場中輸電固定成本分配方法對比研究(Studyonmethodology of fixed transmission costallocation in powermarket)[J].電力技術(shù)經(jīng)濟(jì)(Electric Power TechnologicEconomics),2009,21(5):10-17,31.

    [8]CookW D,KressM.Characterizing an equitable allocation of shared costs:a DEA approach[J].European Journalof OperationalResearch,1999,119(3):652-661.

    [9]孫玉華,曾慶鐸(Sun Yuhua,Zeng Qingduo).基于DEA的二階段網(wǎng)絡(luò)系統(tǒng)的固定成本分?jǐn)偡椒ǎ‵ixed costallocation of two-stage network systems based on DEA)[J].經(jīng)濟(jì)數(shù)學(xué)(Journal of Quantitative Economics),2013,30(1):33-36.

    [10]魏權(quán)齡,許余潔(WeiQuanling,Xu Yujie).具有偏好錐的鏈?zhǔn)骄W(wǎng)絡(luò)DEA模型(Chain network DEAmodelwith preference cone)[J].數(shù)學(xué)的實踐與認(rèn)識(Mathematics in Practiceand Theory),2010,40(24):228-240.

    [11]Estellita Lins MP,Moreira da Silva A C,Lovell C A K. Avoiding infeasibility in DEA modelswith weight restrictions[J].European JournalofOperationalResearch,2007,181(2):956-966.

    [12]張茂勤,李光金,尚文娟(Zhang Maoqin,LiGuangjin,Shang Wenjuan).基于Campos指數(shù)的模糊DEA(A fuzzy DEAmodelbased on Campos index)[J].系統(tǒng)工程理論與實踐(Systems Engineering—Theory&Practice),2004,24(4):41-48.

    [13]陳奎,韋曉廣,??∑迹龋–hen Kui,WeiXiaoguang,Niu Junping,etal).基于AHP-DEA模型的電網(wǎng)規(guī)劃方案綜合評判決策(Comprehensive judgment for power systeMplanning alternatives based on AHP-DEA)[J].電力系統(tǒng)保護(hù)與控制(PowerSysteMProtection and Control),2014,42(21):40-46.

    [14]阿遼沙·葉,祝恩國,成倩,等(Aliaosha·Ye,Zhu Enguo,Cheng Qian,etal).用電設(shè)備安全評估的改進(jìn)區(qū)間層次分析法(Improved intervalanalytic hierarchy processmethod for electrical equipment safety assessment)[J].電力系統(tǒng)及其自動化學(xué)報(Proceedingsof the CSU-EPSA),2015,27(1):32-36.

    [15]周明,鄭雅楠,李庚銀,等(Zhou Ming,Zheng Yanan,Li Gengyin,etal).考慮分時電價和電能質(zhì)量的大用戶轉(zhuǎn)運費用固定成本的綜合分?jǐn)偡椒ǎ⊿ynthetical method forallocatingwheeling fixed cost for large consumers considering time-of-use pricingand powerquality)[J].中國電機(jī)工程學(xué)報(Proceedingsof the CSEE),2008,28(19):125-130.

    [16]黃新,王健,陳海峰(Huang Xin,Wang Jian,Chen Haifeng).火電機(jī)組能耗及排放綜合考核管理方法研究(Comprehensive assessmentmanagementmethod for energy consumption and emissions of thermal power units)[J].華東電力(East China Electric Power),2011,39(10):1585-1587.

    [17]王曉璐(Wang Xiaolu).火電機(jī)組能效評價體系探究(Preliminary Research on Energy EfficientEvaluation systeMfor SupercriticalUnits)[D].北京:華北電力大學(xué)能源動力與機(jī)械工程學(xué)院(Beijing:Schoolof Energy,Power and Mechanical Engineering,North China Electric Power University),2012.

    [18]于海波(Yu Haibo).基于負(fù)載率均衡度的電力系統(tǒng)調(diào)度策略與風(fēng)電規(guī)劃研究(Research of Power SysteMDispatching Strategy and Wind Power Planning Based on Load Rate Balance Degree)[D].哈爾濱∶哈爾濱工業(yè)大學(xué)電氣工程及自動化學(xué)院(Harbin∶Schoolof Electrical Engineering and Automation,Harbin Institute of Technology),2013.

    [19]李勇軍(LiYongjun).基于DEA理論的固定成本分?jǐn)偡椒ㄑ芯浚≧esearch on the Fixed Cost Allocation based on DEA)[D].合肥:中國科學(xué)技術(shù)大學(xué)管理學(xué)院(Hefei:School of Management,University of Science and Technology ofChina),2008.

    [20]高效,彭建春,羅安(Gao Xiao,Peng Jianchun,Luo An).多種交易模式下核仁解分?jǐn)傒旊娋W(wǎng)固定成本(A nucleolus solution for fixed transmission cost allocation under hybrid transactionsmodes)[J].中國電機(jī)工程學(xué)報(Proceedingsof the CSEE),2007,27(10):120-124.

    [21]劉自發(fā),龐鋮鋮,魏建煒,等(Liu Zifa,Pang Chengcheng,Wei Jianwei,etal).基于IAHP和TOPSIS方法的負(fù)荷密度指標(biāo)計算(Index calculation of load density based on IAHPand TOPSISmethods)[J].電力系統(tǒng)自動化(Automation of Electric Power Systems),2012,36(13):56-60.

    [22]肖峻,羅鳳章,王成山,等(Xiao Jun,Luo Fengzhang,Wang Chengshan,etal).區(qū)間層次分析法的權(quán)重求解方法比較研究(Comparison amongmethods ofweights calculation in the interval-based AHP)[J].電力系統(tǒng)及其自動化學(xué)報(Proceedings of the CSU-EPSA),2004,16(3):12-16.

    [23]張伯明,陳壽孫.高等電力網(wǎng)絡(luò)分析[M].北京:清華大學(xué)出版社,1994.

    New Method for Allocating Power Network Fixed Costamong Direct Trades

    ZHOUYing1,WEN Fushuan1,2,ZHUBingquan3,XIANGZhongming3,XU Lizhong3,ZHAO Liang3
    (1.SchoolofElectricalEngineering,Zhejiang University,Hangzhou 310027,China;2.Department of Electrical&Electronic Engineering,Institut TeknologiBrunei,Bandar SeriBegawan BE1410,Brunei;3.State Grid Zhejiang Electric PowerCorporation,Hangzhou 310027,China)

    How to allocate powernetwork fixed cost isan importantbutnotyetwell-solved issue related to direct trades between large consumers and generation companies.Given this background,this issue is systematically analyzed considering variousattributes of direct trades.Specifically,an evaluation systeMfor allocating power network fixed cost is structured first froMthree aspects:the user side,generator side and the impactson the power system.Then,a cooperative gamemodel in the framework of data envelopment analysis is presented to address this issue.In thismodel,weightconstraintsare introduced by employing the intervalanalytic hierarchy process(IAHP)to reflect the preference of decision-makers,and a power network fixed costallocationmodel considering the index weight constraints is next attained.Finally,the IEEE 14-bus systeMisemployed to demonstrate the basic characteristicsof the developedmodel. Simulation results showed that the proposed approach can lead to fair and reasonable allocationswith decisionmakers' preference considered.

    direct trades;fixed costallocation;index system;data envelopmentanalysis;cooperative game;weight constraint

    TM73

    A

    1003-8930(2015)09-0012-09

    10.3969/j.issn.1003-8930.2015.09.003

    周盈(1991—),女,碩士研究生,主要從事電價機(jī)制、大用戶直購電和市場策略方面的研究。Email:zhouying7988@ gmail.com

    2015-03-03;

    2015-04-10

    國家重點基礎(chǔ)研究發(fā)展規(guī)劃(973計劃)資助項目(2013CB228202);國網(wǎng)浙江省電力公司科研項目(5211ZD13000R);國家電網(wǎng)公司總部科技項目(52110115009Q)

    文福拴(1965—),男,博士,教授,博士生導(dǎo)師,主要從事電力系統(tǒng)故障診斷與系統(tǒng)恢復(fù)、電力經(jīng)濟(jì)與電力市場、智能電

    網(wǎng)與電動汽車等方面的研究。Email:fushuan.wen@gmail.com朱炳銓(1967—),男,本科,高級工程師,主要從事電網(wǎng)調(diào)度運行管理工作。Email:hzzbq@sina.com

    猜你喜歡
    分?jǐn)?/a>約束權(quán)重
    系統(tǒng)工程學(xué)報(2021年4期)2021-12-21 06:21:16
    “碳中和”約束下的路徑選擇
    權(quán)重常思“浮名輕”
    約束離散KP方程族的完全Virasoro對稱
    為黨督政勤履職 代民行權(quán)重?fù)?dān)當(dāng)
    基于公約式權(quán)重的截短線性分組碼盲識別方法
    適當(dāng)放手能讓孩子更好地自我約束
    人生十六七(2015年6期)2015-02-28 13:08:38
    無形資產(chǎn)轉(zhuǎn)移定價視角下的獨立交易法與公式分?jǐn)偡ㄖ疇?發(fā)展中國家何去何從?
    電測與儀表(2014年3期)2014-04-04 09:08:14
    層次分析法權(quán)重的計算:基于Lingo的數(shù)學(xué)模型
    河南科技(2014年15期)2014-02-27 14:12:51
    99国产精品免费福利视频| 熟女电影av网| 久久99一区二区三区| 99热这里只有是精品在线观看| 成人国语在线视频| 五月开心婷婷网| 国产午夜精品一二区理论片| 少妇熟女欧美另类| 五月玫瑰六月丁香| 最黄视频免费看| 日本-黄色视频高清免费观看| 男女边摸边吃奶| 青春草视频在线免费观看| 欧美国产精品一级二级三级| 久久热精品热| 久久这里有精品视频免费| 久久精品国产鲁丝片午夜精品| 国产 精品1| 欧美另类一区| 女人精品久久久久毛片| 人体艺术视频欧美日本| 国产日韩欧美亚洲二区| 中文字幕精品免费在线观看视频 | 人妻 亚洲 视频| 精品久久久久久久久av| 国产极品粉嫩免费观看在线 | 国产成人a∨麻豆精品| 18禁裸乳无遮挡动漫免费视频| 寂寞人妻少妇视频99o| 国产精品久久久久久久久免| 国产69精品久久久久777片| 插阴视频在线观看视频| 成人国语在线视频| 少妇丰满av| 免费看不卡的av| 亚洲成人一二三区av| 国产亚洲最大av| 久久久午夜欧美精品| 成人免费观看视频高清| 日韩一区二区视频免费看| 九九爱精品视频在线观看| 精品少妇久久久久久888优播| 久久久久视频综合| 久久久久久久久久久免费av| 九草在线视频观看| 亚洲三级黄色毛片| 成年av动漫网址| 在线观看免费高清a一片| 亚洲精品色激情综合| 成年美女黄网站色视频大全免费 | 亚洲色图综合在线观看| 国产免费又黄又爽又色| 久久ye,这里只有精品| 婷婷色麻豆天堂久久| 成年人午夜在线观看视频| 精品午夜福利在线看| 国产永久视频网站| 国产成人精品一,二区| videossex国产| 国产极品天堂在线| 亚洲第一区二区三区不卡| 国产毛片在线视频| 国产亚洲av片在线观看秒播厂| 亚洲国产精品成人久久小说| 亚洲国产精品一区三区| 国产日韩一区二区三区精品不卡 | 少妇被粗大的猛进出69影院 | 在线精品无人区一区二区三| 亚洲国产色片| 在线观看国产h片| 欧美日韩视频精品一区| 成人手机av| 成人毛片60女人毛片免费| 欧美老熟妇乱子伦牲交| 18+在线观看网站| 日本欧美国产在线视频| 99热网站在线观看| 男的添女的下面高潮视频| 久久鲁丝午夜福利片| 性高湖久久久久久久久免费观看| 中文字幕人妻丝袜制服| 超色免费av| 亚洲精品亚洲一区二区| 欧美精品人与动牲交sv欧美| 亚洲精品,欧美精品| 成人免费观看视频高清| 在线观看一区二区三区激情| 啦啦啦啦在线视频资源| 在线 av 中文字幕| 国产亚洲av片在线观看秒播厂| 精品亚洲成国产av| tube8黄色片| 国产av一区二区精品久久| av在线老鸭窝| 亚洲av成人精品一区久久| 看免费成人av毛片| av在线观看视频网站免费| 亚洲国产av影院在线观看| a级毛片免费高清观看在线播放| 99九九在线精品视频| 欧美精品亚洲一区二区| 国产精品麻豆人妻色哟哟久久| 色5月婷婷丁香| 日韩熟女老妇一区二区性免费视频| 欧美精品人与动牲交sv欧美| 大码成人一级视频| 热re99久久国产66热| 热re99久久精品国产66热6| 成年女人在线观看亚洲视频| 亚洲精品亚洲一区二区| 啦啦啦中文免费视频观看日本| 亚洲精品国产av蜜桃| 麻豆成人av视频| 五月伊人婷婷丁香| 国模一区二区三区四区视频| 国产色婷婷99| 中文欧美无线码| 久久99热这里只频精品6学生| 9色porny在线观看| 中文字幕精品免费在线观看视频 | 最黄视频免费看| 在线观看一区二区三区激情| 亚洲第一区二区三区不卡| av免费观看日本| 国产国语露脸激情在线看| 少妇的逼好多水| 春色校园在线视频观看| 色婷婷久久久亚洲欧美| 国产成人精品一,二区| 免费大片黄手机在线观看| 99国产综合亚洲精品| av播播在线观看一区| 男女国产视频网站| 99热全是精品| 人人妻人人爽人人添夜夜欢视频| 蜜臀久久99精品久久宅男| 日韩大片免费观看网站| 日韩免费高清中文字幕av| 黄色视频在线播放观看不卡| 另类精品久久| 久久久a久久爽久久v久久| 久久久a久久爽久久v久久| 国产在线视频一区二区| 国产一区二区三区av在线| 国产欧美另类精品又又久久亚洲欧美| 久久人人爽人人爽人人片va| 春色校园在线视频观看| 黑人高潮一二区| 热re99久久国产66热| 久久精品久久久久久久性| 免费黄频网站在线观看国产| 久久亚洲国产成人精品v| 交换朋友夫妻互换小说| 欧美精品一区二区免费开放| 国产极品粉嫩免费观看在线 | 在线观看免费视频网站a站| www.av在线官网国产| 91久久精品电影网| 我的女老师完整版在线观看| 校园人妻丝袜中文字幕| 国产精品一区www在线观看| 午夜激情av网站| 超碰97精品在线观看| av专区在线播放| www.av在线官网国产| 一区二区三区精品91| 亚洲综合色网址| 成人无遮挡网站| 久热久热在线精品观看| 精品少妇黑人巨大在线播放| 日本午夜av视频| 国产精品欧美亚洲77777| 日韩免费高清中文字幕av| av不卡在线播放| 亚洲人成网站在线播| 国产精品99久久久久久久久| 一区二区三区精品91| 最近最新中文字幕免费大全7| 在线观看美女被高潮喷水网站| 亚洲精品成人av观看孕妇| 男女无遮挡免费网站观看| 色94色欧美一区二区| 国产精品成人在线| 亚洲精品成人av观看孕妇| 天天影视国产精品| 欧美日韩一区二区视频在线观看视频在线| 色婷婷av一区二区三区视频| 国产成人免费无遮挡视频| 毛片一级片免费看久久久久| 在线亚洲精品国产二区图片欧美 | 永久免费av网站大全| 亚洲国产av影院在线观看| 丰满迷人的少妇在线观看| 久久久久网色| 男女无遮挡免费网站观看| 你懂的网址亚洲精品在线观看| 菩萨蛮人人尽说江南好唐韦庄| 免费大片18禁| 国产精品99久久99久久久不卡 | 精品人妻熟女毛片av久久网站| 午夜av观看不卡| 啦啦啦中文免费视频观看日本| 日本91视频免费播放| 日韩av不卡免费在线播放| 亚洲精品av麻豆狂野| 韩国av在线不卡| 赤兔流量卡办理| 街头女战士在线观看网站| 亚洲一级一片aⅴ在线观看| 精品卡一卡二卡四卡免费| 中文字幕人妻熟人妻熟丝袜美| 欧美最新免费一区二区三区| 成人国产麻豆网| 亚洲国产最新在线播放| 日本欧美视频一区| 边亲边吃奶的免费视频| 午夜激情久久久久久久| 亚洲欧美日韩卡通动漫| a 毛片基地| 久久久久久久国产电影| av视频免费观看在线观看| 插逼视频在线观看| 国产午夜精品久久久久久一区二区三区| 亚洲综合色网址| 毛片一级片免费看久久久久| 热re99久久精品国产66热6| 欧美变态另类bdsm刘玥| 插逼视频在线观看| 日韩不卡一区二区三区视频在线| 国产高清国产精品国产三级| 下体分泌物呈黄色| 国产亚洲精品第一综合不卡 | 制服诱惑二区| 简卡轻食公司| 人体艺术视频欧美日本| 一二三四中文在线观看免费高清| 免费看不卡的av| 99热网站在线观看| 中文字幕制服av| 天堂中文最新版在线下载| 亚洲精品久久午夜乱码| 久久精品熟女亚洲av麻豆精品| 26uuu在线亚洲综合色| 精品一区二区三区视频在线| 熟女电影av网| 精品少妇内射三级| 日日摸夜夜添夜夜添av毛片| 自拍欧美九色日韩亚洲蝌蚪91| 飞空精品影院首页| 寂寞人妻少妇视频99o| 在线观看三级黄色| 如何舔出高潮| 男女边摸边吃奶| 黄片播放在线免费| 久久国产精品男人的天堂亚洲 | 一级片'在线观看视频| 欧美 日韩 精品 国产| 我的老师免费观看完整版| 91精品三级在线观看| 91午夜精品亚洲一区二区三区| 青春草亚洲视频在线观看| av在线老鸭窝| 国产黄频视频在线观看| 3wmmmm亚洲av在线观看| 久久国内精品自在自线图片| 欧美老熟妇乱子伦牲交| 亚洲av成人精品一二三区| 亚洲综合精品二区| 蜜桃在线观看..| 亚洲色图 男人天堂 中文字幕 | 国产日韩一区二区三区精品不卡 | 亚洲五月色婷婷综合| 99国产精品免费福利视频| 精品一区二区免费观看| 熟女人妻精品中文字幕| 日本-黄色视频高清免费观看| av在线播放精品| 精品久久久噜噜| 婷婷色综合www| 中文字幕亚洲精品专区| 91国产中文字幕| 国内精品宾馆在线| 欧美日韩综合久久久久久| 丰满迷人的少妇在线观看| 精品一区二区三卡| 爱豆传媒免费全集在线观看| 丁香六月天网| 草草在线视频免费看| 啦啦啦啦在线视频资源| 一区二区三区精品91| a级毛片黄视频| 亚洲av二区三区四区| 青春草亚洲视频在线观看| 青春草视频在线免费观看| 一区在线观看完整版| 狂野欧美白嫩少妇大欣赏| 亚洲第一区二区三区不卡| 嫩草影院入口| 人人妻人人爽人人添夜夜欢视频| 中国国产av一级| 男女无遮挡免费网站观看| 如日韩欧美国产精品一区二区三区 | 中文字幕av电影在线播放| 赤兔流量卡办理| 自拍欧美九色日韩亚洲蝌蚪91| 国产深夜福利视频在线观看| 亚洲欧美成人综合另类久久久| 久久久精品区二区三区| 亚洲精品日韩在线中文字幕| 97精品久久久久久久久久精品| 国产精品久久久久久久电影| 精品国产国语对白av| 十八禁网站网址无遮挡| 国产精品女同一区二区软件| 国产毛片在线视频| 免费大片18禁| 人妻一区二区av| 男女边吃奶边做爰视频| 久久狼人影院| 中文精品一卡2卡3卡4更新| 国产亚洲精品第一综合不卡 | 亚洲av电影在线观看一区二区三区| 日韩精品免费视频一区二区三区 | av天堂久久9| 日韩一区二区视频免费看| 各种免费的搞黄视频| 内地一区二区视频在线| 99热国产这里只有精品6| 韩国高清视频一区二区三区| 国产精品国产三级专区第一集| 五月天丁香电影| 99国产综合亚洲精品| 99久久综合免费| 亚洲精品美女久久av网站| 亚洲国产精品一区三区| 国产成人一区二区在线| 久久影院123| 一区二区日韩欧美中文字幕 | 亚洲精品国产色婷婷电影| 成年女人在线观看亚洲视频| 男女边吃奶边做爰视频| 另类亚洲欧美激情| 多毛熟女@视频| 欧美一级a爱片免费观看看| 久久久久久伊人网av| 老熟女久久久| 国产成人91sexporn| 99热这里只有是精品在线观看| 视频在线观看一区二区三区| 少妇人妻精品综合一区二区| 视频在线观看一区二区三区| 国产视频内射| 亚洲熟女精品中文字幕| 亚洲精品456在线播放app| 国产精品无大码| 欧美丝袜亚洲另类| 69精品国产乱码久久久| 成年av动漫网址| 欧美日韩视频精品一区| 纯流量卡能插随身wifi吗| 如日韩欧美国产精品一区二区三区 | 久久毛片免费看一区二区三区| 亚洲欧洲日产国产| 人妻 亚洲 视频| 国产亚洲午夜精品一区二区久久| 热99久久久久精品小说推荐| 18禁在线播放成人免费| 国产成人freesex在线| 黄色怎么调成土黄色| 亚洲久久久国产精品| 国产成人av激情在线播放 | 男女边吃奶边做爰视频| 18+在线观看网站| 久久ye,这里只有精品| 午夜免费男女啪啪视频观看| 美女xxoo啪啪120秒动态图| 国产不卡av网站在线观看| 狂野欧美激情性xxxx在线观看| 欧美97在线视频| 国产亚洲欧美精品永久| 国产av一区二区精品久久| a级毛片免费高清观看在线播放| 成人午夜精彩视频在线观看| 国产片内射在线| 蜜桃在线观看..| 亚洲av电影在线观看一区二区三区| 熟女av电影| 久久综合国产亚洲精品| 十八禁网站网址无遮挡| 毛片一级片免费看久久久久| 在线天堂最新版资源| 国产亚洲精品久久久com| 永久网站在线| 狂野欧美激情性bbbbbb| 18在线观看网站| 亚洲欧美日韩另类电影网站| 如日韩欧美国产精品一区二区三区 | 国产精品不卡视频一区二区| 黑人欧美特级aaaaaa片| 五月开心婷婷网| 美女内射精品一级片tv| 中文字幕人妻丝袜制服| 黄色视频在线播放观看不卡| 蜜臀久久99精品久久宅男| 久久久精品区二区三区| 26uuu在线亚洲综合色| av线在线观看网站| 中国三级夫妇交换| 精品一区二区三卡| 搡老乐熟女国产| a级片在线免费高清观看视频| 制服人妻中文乱码| 十八禁网站网址无遮挡| 久久午夜福利片| 精品国产一区二区三区久久久樱花| 春色校园在线视频观看| av网站免费在线观看视频| 水蜜桃什么品种好| 只有这里有精品99| 欧美成人精品欧美一级黄| 国产永久视频网站| 性色av一级| 日本av免费视频播放| 欧美亚洲日本最大视频资源| 亚洲国产av影院在线观看| 午夜福利视频在线观看免费| 色哟哟·www| 在线播放无遮挡| 精品少妇久久久久久888优播| 少妇高潮的动态图| 精品久久久精品久久久| 亚洲av日韩在线播放| 午夜激情av网站| 超色免费av| 特大巨黑吊av在线直播| 桃花免费在线播放| 国产淫语在线视频| 久久人人爽av亚洲精品天堂| 岛国毛片在线播放| 夜夜骑夜夜射夜夜干| 亚洲国产精品一区二区三区在线| a级毛片黄视频| freevideosex欧美| 亚洲精品国产av蜜桃| 亚洲av.av天堂| 亚洲av成人精品一区久久| 久久久国产一区二区| 两个人免费观看高清视频| 搡老乐熟女国产| 人妻制服诱惑在线中文字幕| 午夜91福利影院| 我的老师免费观看完整版| 亚洲精品第二区| 欧美精品一区二区免费开放| 免费黄网站久久成人精品| 亚洲av男天堂| 日韩电影二区| 国产成人av激情在线播放 | 在线观看美女被高潮喷水网站| 性高湖久久久久久久久免费观看| 亚洲av免费高清在线观看| av网站免费在线观看视频| 国产69精品久久久久777片| 亚洲欧洲国产日韩| 欧美精品高潮呻吟av久久| a级毛片免费高清观看在线播放| 99国产精品免费福利视频| 又大又黄又爽视频免费| 999精品在线视频| 人人妻人人澡人人看| 国产日韩欧美亚洲二区| av一本久久久久| 看免费成人av毛片| 国产亚洲精品久久久com| 亚洲精品乱码久久久久久按摩| 国产69精品久久久久777片| 草草在线视频免费看| 日韩视频在线欧美| 天堂8中文在线网| 黑丝袜美女国产一区| 丝袜脚勾引网站| 欧美日韩在线观看h| 亚洲精品中文字幕在线视频| 精品久久久久久久久亚洲| 亚洲少妇的诱惑av| 男女边吃奶边做爰视频| 免费大片黄手机在线观看| 我要看黄色一级片免费的| av网站免费在线观看视频| 成人二区视频| 国产精品国产三级国产av玫瑰| 亚洲精品久久成人aⅴ小说 | 999精品在线视频| 国产午夜精品久久久久久一区二区三区| 人妻系列 视频| 日韩一本色道免费dvd| 夜夜骑夜夜射夜夜干| 五月开心婷婷网| 黄色怎么调成土黄色| 伦精品一区二区三区| 中文字幕最新亚洲高清| 久久久精品免费免费高清| 国产亚洲欧美精品永久| 国产成人精品一,二区| 我的女老师完整版在线观看| 亚洲精品乱码久久久v下载方式| 欧美人与性动交α欧美精品济南到 | 国产色婷婷99| 一本色道久久久久久精品综合| 新久久久久国产一级毛片| 亚洲伊人久久精品综合| 亚洲人成网站在线播| 精品久久久久久久久亚洲| 欧美成人精品欧美一级黄| 成人国产麻豆网| 久久精品熟女亚洲av麻豆精品| 日本av手机在线免费观看| 中文字幕人妻丝袜制服| 热re99久久精品国产66热6| 七月丁香在线播放| a级片在线免费高清观看视频| 精品卡一卡二卡四卡免费| 天天躁夜夜躁狠狠久久av| 国产欧美日韩一区二区三区在线 | 少妇的逼水好多| 亚洲欧洲日产国产| 狂野欧美白嫩少妇大欣赏| 国产一级毛片在线| 亚洲精品乱码久久久久久按摩| 99久久精品国产国产毛片| 欧美三级亚洲精品| 国语对白做爰xxxⅹ性视频网站| 99国产精品免费福利视频| 亚洲欧美日韩卡通动漫| 国产成人freesex在线| 久久99热这里只频精品6学生| 中国美白少妇内射xxxbb| av一本久久久久| 国产不卡av网站在线观看| 中国国产av一级| 亚洲,欧美,日韩| 国产成人av激情在线播放 | 精品人妻偷拍中文字幕| 99热6这里只有精品| 国产精品国产三级专区第一集| 天堂俺去俺来也www色官网| 街头女战士在线观看网站| 亚洲成人一二三区av| 国产一区二区三区av在线| 国产欧美日韩一区二区三区在线 | av线在线观看网站| 女人精品久久久久毛片| 国产高清三级在线| 午夜激情久久久久久久| 另类亚洲欧美激情| 一本久久精品| 国产黄频视频在线观看| 亚洲av成人精品一区久久| 久久精品久久精品一区二区三区| 中文天堂在线官网| av线在线观看网站| 国产日韩一区二区三区精品不卡 | 男男h啪啪无遮挡| 22中文网久久字幕| 高清在线视频一区二区三区| 18禁在线播放成人免费| 国产免费福利视频在线观看| 亚洲国产欧美日韩在线播放| 久久久国产一区二区| 国产精品.久久久| 一级黄片播放器| 美女国产视频在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | av福利片在线| 赤兔流量卡办理| 一本久久精品| 观看av在线不卡| 两个人免费观看高清视频| 国模一区二区三区四区视频| 欧美精品一区二区免费开放| 男女无遮挡免费网站观看| 国产成人a∨麻豆精品| 汤姆久久久久久久影院中文字幕| 国产日韩欧美亚洲二区| 人人妻人人添人人爽欧美一区卜| 日本黄大片高清| 熟妇人妻不卡中文字幕| 青春草亚洲视频在线观看| 国产免费福利视频在线观看| 免费黄网站久久成人精品| 在线观看国产h片| 大话2 男鬼变身卡| 欧美成人午夜免费资源| 美女xxoo啪啪120秒动态图| 久久99一区二区三区| 久久精品国产亚洲av天美| 国产精品一区二区在线不卡| 欧美xxⅹ黑人| 九九爱精品视频在线观看| 丰满迷人的少妇在线观看| 成人毛片60女人毛片免费| 亚洲国产精品国产精品| 成人漫画全彩无遮挡| 午夜福利在线观看免费完整高清在| 欧美日韩亚洲高清精品| 欧美国产精品一级二级三级| 欧美日韩一区二区视频在线观看视频在线| 欧美日韩亚洲高清精品| 狂野欧美激情性xxxx在线观看| 寂寞人妻少妇视频99o| 老熟女久久久| 日本av手机在线免费观看| 亚洲精品亚洲一区二区| 一级爰片在线观看| 久久 成人 亚洲|