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

    考慮新能源爬坡的風光火耦合系統(tǒng)源荷匹配性分析及容量優(yōu)化配置

    2024-01-15 11:58:40夏芹芹羅永捷王榮茂羅桓桓李金燦周念成王強鋼
    上海交通大學學報 2024年1期
    關鍵詞:爬坡火電風光

    夏芹芹, 羅永捷, 王榮茂, 鄒 堯, 羅桓桓 李金燦, 周念成, 王強鋼

    (1.重慶大學 輸配電裝備及系統(tǒng)安全與新技術國家重點實驗室,重慶 400044;2.國網(wǎng)遼寧省電力有限公司,沈陽 110002; 3.廣西電網(wǎng)有限責任公司,南寧 530012)

    當前,國務院出臺關于碳達峰、碳中和的重大戰(zhàn)略決策[1-2],全面推進風力發(fā)電(風電)、光伏發(fā)電(光伏)大規(guī)模開發(fā)和高質(zhì)量發(fā)展,推動建立風、光互補調(diào)節(jié)的新能源發(fā)電基地.國家發(fā)展改革委員會和國家能源局指導相關企事業(yè)單位積極探索“風光水火儲一體化”實施,以提升各類能源互補協(xié)調(diào)能力[3].因此,研究風電、光伏的高效利用對推動新能源高質(zhì)量發(fā)展、能源互補協(xié)調(diào)、降低火電機組碳排放具有現(xiàn)實意義.

    然而,風光出力波動性和不確定性給電力系統(tǒng)帶來挑戰(zhàn).一方面,風光出力波動性表現(xiàn)在預測與實際出力的誤差,這些誤差需要電力系統(tǒng)的可控靈活資源進行調(diào)控,如果調(diào)控能力不足可能引發(fā)棄風棄光或供電不足的問題[4];另一方面,風光出力波動性也表現(xiàn)在短時間內(nèi)輸出功率發(fā)生大幅度單方向變化,即爬坡事件[5].風光爬坡事件主要與區(qū)域氣象資源有關[6-7].爬坡事件可能導致電力系統(tǒng)難以消納新能源,并可能威脅系統(tǒng)安全運行.因此,研究風電、光伏爬坡事件,減少其影響是新能源發(fā)電研究的重要課題之一.

    風光不確定性首先可通過同地區(qū)的風光互補和源荷匹配解決.文獻[8]中提出描述風光出力與負荷貼近程度的負荷匹配和互補指標,對新能源互補進行評價;文獻[9]中考慮風光2種新能源之間的相關性,提出大規(guī)模風光互補的多目標電網(wǎng)擴展規(guī)劃方法;文獻[10]中針對風光互補發(fā)電系統(tǒng),提出一種改進的優(yōu)化規(guī)模確定方法,以實現(xiàn)注入電網(wǎng)的功率波動更小的目標.

    除風光互補和源荷匹配外,風電、光伏不確定性主要由電力系統(tǒng)的可控靈活性發(fā)電資源進行消納.在我國北方地區(qū),火力發(fā)電(火電)是發(fā)電側主要靈活性資源[11].當前,通過機組改造和政策支持,火電已具備更強的調(diào)控能力以更好地消納新能源[12].與僅具備常規(guī)調(diào)峰(RPR)能力的機組相比,具備深度調(diào)峰(DPR)能力的火電機組能夠低于正常最小出力運行[13],以提供更多調(diào)峰能力.此外,我國東北地區(qū)已有調(diào)峰輔助服務政策,使火電可通過該服務獲得額外利益[14].

    火電與新能源的配合可采用耦合系統(tǒng)的形式.耦合系統(tǒng)是一種將物理距離較近的新能源和火電機組整合組成的集成系統(tǒng).因其通過同一并網(wǎng)點連接到主電網(wǎng),對外可視作單個電源由調(diào)度部門進行調(diào)控,系統(tǒng)運行者根據(jù)總出力對其內(nèi)部電源出力進行分配,以提高各電源的協(xié)同性,并通過輔助服務獲得更多經(jīng)濟效益.文獻[15]中考慮新能源與火電耦合系統(tǒng)的耦合機制,建立多維協(xié)同性能評價體系,并提出耦合潛力辨識方法,驗證耦合系統(tǒng)在電源協(xié)同互補方面的優(yōu)勢.文獻[16]中提出一種通過階梯式爬坡率進行表征的火電機組深度調(diào)峰運行模型,應用于含風電和火電耦合系統(tǒng)的優(yōu)化調(diào)度模型中,并通過仿真驗證含風電和火電的耦合系統(tǒng)中,通過充分挖掘火電的深度調(diào)峰特性可提高可再生能源的利用率.

    在新能源與常規(guī)電源耦合的規(guī)劃方面,相關研究主要考慮綜合能源系統(tǒng)或電源打捆的形式.文獻[17]中針對基于冷熱電聯(lián)供系統(tǒng)耦合的電力、天然氣區(qū)域綜合能源系統(tǒng),采用多階段規(guī)劃和多場景分析,獲得線路、管網(wǎng)規(guī)劃及冷熱電聯(lián)供系統(tǒng)的定容選址方案.文獻[18]中構建了一個考慮可靠性與成本的多目標協(xié)調(diào)電-氣耦合系統(tǒng)容量規(guī)劃模型,以實現(xiàn)發(fā)電機組與儲能裝置的選址定容.文獻[19]中提出一種基于時序運行模擬的風火打捆容量最優(yōu)配比整定方法,對比大范圍區(qū)域中不同風火配比下的風電消納能力及運行經(jīng)濟性數(shù)據(jù),給出兼顧環(huán)保性、安全性與經(jīng)濟性的風火配比.文獻[20]中以保證輸電工程的安全運行為目標,構建了一個考慮大區(qū)域新能源穿透功率極限、新能源出力反調(diào)峰及大波動風險的區(qū)域風光火打捆多目標優(yōu)化配置模型.

    總結上述文獻可知,目前已有的相關研究主要集中在多能源耦合系統(tǒng)容量規(guī)劃或發(fā)電側大范圍新能源與傳統(tǒng)電源打捆外送的宏觀層面容量規(guī)劃,缺少針對發(fā)電側含風、光、火且距離較近的多電源耦合系統(tǒng)的精細化源荷協(xié)同配合分析,缺少考慮風電、光伏爬坡特性和源荷匹配特性與耦合系統(tǒng)的相互影響,無法精確反映風光火資源在耦合系統(tǒng)內(nèi)部電源間、耦合系統(tǒng)和外部負荷需求間的協(xié)同作用效果.此外,未有研究針對給定火電機組如何充分挖掘火電調(diào)峰能力及利用源荷特性,合理配置風電、光伏容量,組成高效協(xié)同風光火耦合系統(tǒng)作定量分析.

    針對上述不足,以含風電、光伏、火電機組的耦合系統(tǒng)為基礎,考慮新能源爬坡特性和負荷關鍵指標,結合相關文件中控制新增火電、加快現(xiàn)役機組靈活性改造、推動火電向基礎保障性和調(diào)節(jié)性電源轉型的政策導向[1],研究在已有火電機組的情況下,風電、光伏的容量優(yōu)化配置方法.首先,分析風光火耦合系統(tǒng)的運行模式,對相關不確定性變量進行表征;其次,提取風電、光伏及其耦合后的年爬坡事件總發(fā)生時間,結合關鍵負荷特性指標進行源荷匹配性評價;再次,以綜合收益最優(yōu)為目標,建立考慮源荷特性、源荷匹配和多種成本、收益的不確定性風光火耦合系統(tǒng)容量優(yōu)化配置模型;最后,通過我國遼寧地區(qū)的實際算例,驗證所提容量優(yōu)化配置模型的有效性,并分析容量優(yōu)化配置結果、爬坡事件及源荷匹配性的影響因素,為耦合系統(tǒng)規(guī)劃人員提供不同情形中風、光裝機容量配置的參考.本文的創(chuàng)新性主要在于:針對發(fā)電側風光火耦合系統(tǒng),充分考慮源荷關鍵特性,確定已有火電機組如何配置風電、光伏的最優(yōu)裝機容量,解決發(fā)電基地工程實際中電源協(xié)調(diào)配置問題;針對區(qū)域特點,從源、荷兩個方面挖掘和利用風光火耦合系統(tǒng)的規(guī)律,為新能源最優(yōu)容量配置提供參考與提升建議.

    1 風光火耦合系統(tǒng)的不確定性建模及其關鍵特征分析

    1.1 風光火耦合系統(tǒng)的運行模式

    風光火耦合系統(tǒng)由物理距離相近的火電、風電和光伏電源組成,并通過同一點接入主電網(wǎng),其主要結構如圖1所示.風光火耦合系統(tǒng)運行模式為:①電力系統(tǒng)調(diào)度部門綜合考慮風、光、負荷預測數(shù)據(jù)和火電特性,制訂負荷需求計劃曲線;②耦合系統(tǒng)以自身經(jīng)濟性最優(yōu)原則,對外根據(jù)負荷需求計劃曲線發(fā)電,對內(nèi)協(xié)調(diào)各電源出力;③耦合系統(tǒng)可從外部電力系統(tǒng)購買一定電量,也可進行一定的棄風、棄光;④耦合系統(tǒng)可參與調(diào)峰輔助服務獲得收益.在容量配置階段,調(diào)度部門可與耦合系統(tǒng)進行配合,針對性調(diào)整負荷需求曲線以實現(xiàn)源側和負荷側的匹配及雙向高效利用.

    圖1 風光火耦合系統(tǒng)示意圖

    1.2 不確定性的表征

    本文考慮的不確定性變量包括風電、光伏發(fā)電功率和負荷需求,采用場景法表征不確定性.場景法的主要思想如下:①對每種類型的相互獨立不確定性變量或具有相關性的不確定性變量,根據(jù)其相應概率密度函數(shù)或聯(lián)合概率密度函數(shù),獲得代表性場景及場景的相應概率[21-22];②根據(jù)不同類型的相互獨立不確定性變量概率密度函數(shù)及具有相關性的不確定性變量聯(lián)合概率密度函數(shù),通過場景生成、聚合和削減,得到多個包含所有類型不確定性變量的場景及其相應概率,以此將不確定性問題轉換為多個確定性問題;由于得到的不確定性變量場景概率之和為1,最終可通過期望值體現(xiàn)在待求問題中[21].為表征風、光出力和負荷的不確定性,采用場景法的步驟如下:通過風、光發(fā)電功率和負荷1年12個月的月度數(shù)據(jù),采用削減方法,每月生成一條具有代表性的典型日每小時期望風、光出力和負荷需求曲線,如圖2所示.考慮風、光出力、負荷需求之間相互獨立,采用正態(tài)分布表征每小時實際風、光出力、負荷需求與期望值的誤差,其中風、光出力、負荷需求對應正態(tài)分布的均值為其每小時期望值,即圖2中的值,方差與期望值和最大出力有關[16].使用上文所述場景法獲得同時考慮風、光出力和負荷整體曲線的代表性場景及其相應概率,用于耦合系統(tǒng)容量優(yōu)化配置仿真,每小時最終場景數(shù)為5個,如圖3所示.使用MATLAB 進行場景生成和聚合,GAMS軟件提供的Scenred[23]模塊用于場景削減.結合12個典型日場景的風、光出力、負荷曲線及每個典型日中每個時刻的5個場景計算,每種不確定性變量的總場景數(shù)量為1 440個,根據(jù)文獻[21]中的研究,場景削減后的精度已經(jīng)過驗證,故綜合考慮計算量與容量優(yōu)化配置的目的,經(jīng)削減后的場景數(shù)量能夠滿足精度需求.場景法的更多詳情可參考文獻[21-23].

    圖2 每月典型日的期望風電、光伏出力和負荷需求曲線

    圖3 場景法得到的每月典型日的風電、光伏出力和負荷需求具體場景

    1.3 風電/光伏的爬坡事件分析

    爬坡事件可表征風光出力的大尺度波動,是分析耦合系統(tǒng)所在地區(qū)風光資源,研究其出力互補性的重要手段.根據(jù)最優(yōu)爬坡時段,同時考慮不確定性,對每個具體場景的風電、光伏及其耦合后的爬坡事件進行辨識分析.

    目前對于爬坡事件主要有4種定義形式[6],其核心內(nèi)容均為表征給定觀測時間段內(nèi)功率變化是否超過閾值.在給定火電裝機容量配置風電、光伏裝機容量的耦合系統(tǒng)中,因靈活性調(diào)節(jié)資源固定,所關心的爬坡事件為觀測時間內(nèi)新能源發(fā)電功率在單位時間內(nèi)是否超過給定的實際數(shù)值.因此,風電、光伏爬坡事件的表征如下:給定時間和功率的序列{(t1,P1),(t2,P2),…,(tn,Pn)},如果在時間段[ts,te]內(nèi),存在

    (Pe-Ps)/(te-ts)>|Ra|

    (1)

    則認為在起始時刻ts至結束時刻te時間段內(nèi)發(fā)生爬坡事件.式中:Pe、Ps分別為te、ts時刻的功率;|Ra|為爬坡率閾值,當(Pe-Ps)/(te-ts)>Ra時,發(fā)生上爬坡;當(Pe-Ps)/(te-ts)<-Ra時,發(fā)生下爬坡.

    (2)

    式中:J(ti,tj)是在時間間隔[ti,tj]內(nèi)獲得的最大分數(shù),其值可通過tj-ti個子問題的最大值計算;W(ti,tk)是時間間隔[ti,tk]的權重,其值與發(fā)生爬坡事件的時間段長度有關,時間段長度越長,其值越大.更多關于最優(yōu)爬坡時段的計算可參考文獻[24].

    1.4 負荷的關鍵特征選取

    風光火耦合系統(tǒng)的供電能力需要與耦合點外部的負荷需求相匹配才能夠更好地發(fā)揮耦合系統(tǒng)中發(fā)電資源的作用.因耦合系統(tǒng)同時為本地負荷和上級電力系統(tǒng)供電,針對其負荷需求曲線的制訂雖具有一定靈活調(diào)整空間,但負荷需求曲線的部分特性仍主要與耦合系統(tǒng)所在地區(qū)有關,難以改變.針對上述問題,從我國常用負荷特性指標[25]中,選擇年最大負荷、年峰谷差率及其對風光火耦合系統(tǒng)容量優(yōu)化配置的結果影響進行研究.這些指標可隨負荷需求曲線的制訂直接改變,且在耦合系統(tǒng)容量優(yōu)化配置研究中具有代表性.其中,年最大負荷涉及耦合系統(tǒng)最大供電能力,年峰谷差率與耦合系統(tǒng)調(diào)節(jié)能力相關.

    1.5 源荷匹配性評價

    為分析耦合系統(tǒng)所在區(qū)域的源荷耦合程度及影響因素,選取已有研究中的2個典型源荷評價指標,再針對其局限性提出額外2個源荷評價指標,共同對風光火耦合系統(tǒng)源荷匹配性進行評價.其中,選取的指標及其計算方法如下:

    (1) 源荷追蹤系數(shù)[26].表征風、光出力曲線與負荷需求曲線的貼合程度,其值越大表示貼合程度越高.其公式表示為

    (3)

    (2) 源荷變化系數(shù)[26].表征風、光出力曲線與負荷需求曲線差值相對于其平均值的波動水平,其值越大表示曲線整體較平均值的波動水平越大.其公式表示為

    (4)

    上述指標主要表征源荷匹配后曲線的整體趨勢,無法反映源荷匹配后曲線在相鄰時刻的波動水平以分析火電機組爬坡壓力和調(diào)節(jié)能力.因此,提出源荷平均絕對值變化率和源荷動態(tài)時間彎曲距離指標,作為上述指標的補充,具體如下:

    (1) 源荷平均絕對值變化率.表征風、光出力曲線與負荷需求曲線差值在相鄰時刻的波動率,其值越大表示曲線相鄰時刻波動越大.其公式為

    (5)

    (2) 源荷動態(tài)時間彎曲距離.表征風、光出力曲線與負荷需求曲線形態(tài)的相似程度,其值越小表示兩條曲線越相似[27],即相鄰時刻曲線差值變化越小.該指標通過迭代計算為

    (6)

    式中:Li,j為i時刻和j時刻的動態(tài)時間彎曲距離,L0,0=0,L0,j=Li,0=+∞;di,j為i時刻和j時刻風光出力與負荷需求的歐氏距離;I4為源荷動態(tài)時間彎曲距離.

    所選取的源荷匹配評價指標具有顯著針對性,從源荷曲線數(shù)值差異、整體波動、相鄰時間波動和曲線時序相似度4個具體方面直觀反映源荷匹配中的不同關鍵特征.這些評價指標均可獨立使用,不依賴其指標之間的協(xié)助判斷,不受指標之間相關性的影響.

    2 風光火耦合發(fā)電系統(tǒng)容量優(yōu)化配置模型

    從耦合系統(tǒng)運營者角度,風光火耦合系統(tǒng)容量規(guī)劃目標是最大化經(jīng)濟效益,結合上文所述提倡減少新增火電機組的原則,風光火耦合系統(tǒng)容量配置主要考慮我國東北缺少水電地區(qū)中已投產(chǎn)的火電機組配置風光機組容量的場景,并以最大化經(jīng)濟效益為目標.

    2.1 目標函數(shù)

    考慮風光出力和負荷需求的不確定性,建立一個兩階段不確定性風光火耦合系統(tǒng)容量規(guī)劃模型,在考慮耦合系統(tǒng)發(fā)電收益、調(diào)峰輔助服務收益、投資成本、運行維護成本和污染排放成本的條件下,實現(xiàn)年化總收入(FA)最大,目標函數(shù)如下:

    maxFA=FGI-(FOP,d+FOP,u)-(FTE,d+FTE,u)+

    FWTI,u+FPVI,u-FWTC,u-FPVC,u-Fgs,u

    (7)

    目標函數(shù)含有兩個決策階段,第1階段使用下標d表示,第2階段使用下標u表示.其中第1階段是確定性的,該階段的決策變量在不確定性被考慮前就已被確定;第2階段是不確定性的,該階段的決策變量是實時可變動資源對不確定因素在第1階段中產(chǎn)生影響后進行的修正.使用不同場景法表征模型中的不確定性因素,并在最終通過計算期望值將第2階段中的場景由不確定性轉化為確定性.式(7)中各部分計算的均為年化收益期望值.

    式(7)中FGI為火電機組售電收益,與耦合系統(tǒng)參與調(diào)峰輔助服務的情況有關.根據(jù)東北電力市場輔助服務相關規(guī)定[14],當運行的火電機組發(fā)電功率小于補償基準時,可獲調(diào)峰輔助服務補償.風光火耦合系統(tǒng)可充分利用風光資源獲得調(diào)峰輔助服務補償,取得競爭優(yōu)勢.火電機組的發(fā)電收益由下式計算:

    (8)

    (9)

    pTH_EXP(g,tM)=PTH(g,tM)+

    FGI(g,tM)=

    (10)

    式(7)中FOP,d+FOP,u為火電機組運行成本,考慮了火電機組常規(guī)調(diào)峰和深度調(diào)峰兩種狀態(tài).常規(guī)調(diào)峰狀態(tài)運行費用考慮了機組單位發(fā)電量的煤耗和維護費用;深度調(diào)峰狀態(tài)由于火電機組工況較差,運行費用還需考慮此狀態(tài)的額外損耗費用,運行成本通過下式計算:

    (11)

    (12)

    式中:aOP、bOP、cOP、aDPR、bDPR為擬合函數(shù)中的計算系數(shù).

    式(7)中FTE,d+FTE,u為火電機組污染排放處理費用:

    (13)

    式中:CTE(n)、ρTE(n)、τTE(n)分別為火電機組第n種污染物的單位當量處理費用、單位發(fā)電量排放系數(shù)和污染當量值;N為污染物總數(shù),考慮的污染物為煙塵、SO2和NOx.

    式(7)中FWTI,u、FPVI,u分別為風電、光伏的售電收益:

    (14)

    (15)

    式中:PWT(tM,s)、PPV(tM,s)分別為考慮不同場景不確定性波動的風、光發(fā)電功率;CWT、CPV分別為風、光單位發(fā)電量售電收益.

    式(7)中FWTC,u、FPVC,u分別為風、光的投資和維護成本.風電、光伏在壽命周期內(nèi)維護成本相對固定且費用較低,同投資成本類似,可折算至單位裝機容量進行表征[29]:

    (16)

    (17)

    式中:SWT、SPV分別為風電、光伏的規(guī)劃裝機容量;IWT、IPV分別為風電、光伏單位裝機容量的投資和維護成本;yWT、yPV分別為風電、光伏的使用年限;r為年利率.

    式(7)中Fgs,u為耦合系統(tǒng)無法滿足負荷需求時從外部電力系統(tǒng)的購電費用:

    Fgs,u=

    (18)

    式中:Cgs為從外部電力系統(tǒng)購買單位電量的費用;Pgs(tM,s)為不同場景下耦合系統(tǒng)從外部電力系統(tǒng)購買的發(fā)電功率,即耦合系統(tǒng)發(fā)電功率缺額.

    2.2 約束條件

    風光火耦合系統(tǒng)容量優(yōu)化配置的約束條件包括功率平衡約束,電源相關發(fā)電約束,以及源荷匹配限制爬坡事件約束.約束條件針對每個典型日獨立列寫,不同月份典型日之間的約束無關聯(lián).

    2.2.1功率平衡約束 功率平衡約束反映耦合系統(tǒng)內(nèi)部的發(fā)電功率應滿足耦合點的負荷需求:

    PWT(t,s)=PL(t,s)-Pgs(t,s)

    (19)

    Pgs(t,s)應小于該時刻負荷需求預測值PL_pre(t)的一定比例,即Pgs(t,s)≤φLPL_pre(t),其中φL為比例系數(shù).

    2.2.2火電機組發(fā)電約束 火電機組發(fā)電約束主要包括運行狀態(tài)約束和爬坡約束:

    (20)

    火電機組的實時靈活性調(diào)節(jié)資源PTHR(g,t,s)會受到火電機組的爬坡率,當前出力狀態(tài)、出力范圍限制:

    (21)

    2.2.3風電/光伏發(fā)電約束 耦合系統(tǒng)中的風電和光伏在發(fā)電時可通過降低利用率以滿足功率平衡:

    (22)

    2.2.4風光爬坡事件限制約束 為限制耦合系統(tǒng)中的風電和光伏安裝容量過大引起爬坡事件增加,利用源荷耦合對相鄰時刻的凈功率差進行限制:

    |(PWT(t″,s)+PPV(t″,s)-PL(t″,s))-

    (PWT(t″,s)+PPV(t′,s)-PL(t′,s))|/

    (t″-t′)≥σR

    (23)

    式中:1≤t′

    2.3 決策變量

    建立的風光火耦合系統(tǒng)容量優(yōu)化配置模型中包括第1階段和第2階段決策變量.第1階段決策變量包括火電確定性發(fā)電功率PTH(g,t),風電、光伏的裝機容量SPV、SWT.第2階段決策變量包括火電靈活性發(fā)電功率PTHR(g,t,s),風電、光伏的實時出力PWT(t,s)、PPV(t,s),耦合系統(tǒng)供電缺額Pgs(t,s).

    3 算例分析

    3.1 算例基礎數(shù)據(jù)

    以我國遼寧某區(qū)域具備形成圖1所示耦合系統(tǒng)條件的2臺600 MW火電機組為算例進行風光火耦合系統(tǒng)容量優(yōu)化配置仿真.風、光出力及負荷需求曲線使用圖3所示的12個典型日場景及其相應概率,作為容量優(yōu)化配置模型的輸入之一.本算例的調(diào)峰輔助服務相關政策參考文獻[14],火電機組的污染排放相關參數(shù)參考文獻[16],火電機組常規(guī)運行狀態(tài)和深度調(diào)峰狀態(tài)的最小出力為機組容量的50%和30%,爬坡率分別為9 MW/min和6 MW/min;源荷耦合后相鄰時刻的凈功率差限值σR為 100 MW;風電、光伏的裝機與運維綜合成本分別為 1.615 8×107元/MW和4.001×106元/MW,售電價格分別為850元/(MW·h)和740元/(MW·h),壽命分別為30 a和25 a;年利率為5%.為避免耦合系統(tǒng)中電源因各類情況無法供電可能導致的較大負荷缺額,最大負荷值僅考慮略大于或等于火電機組總裝機容量的情形.

    建立的風光火耦合系統(tǒng)容量優(yōu)化配置模型為混合整數(shù)非線性規(guī)劃模型.為方便求解,通過分段線性化和大M法將優(yōu)化模型中整數(shù)變量松弛后[30]使用商業(yè)求解器Gurobi[31]進行求解,展示的結果均為期望值.

    3.2 容量優(yōu)化配置結果的影響因素分析

    為研究負荷需求曲線最大值對容量優(yōu)化配置結果的影響,基于區(qū)域原始負荷需求曲線,結合1.4節(jié)中選取的兩個負荷特性指標,考慮該區(qū)域的源荷特性,調(diào)整負荷需求曲線,研究上述因素對耦合系統(tǒng)容量優(yōu)化配置結果的影響,結果如表1所示.需要注意的是,負荷曲線基于源荷特性僅作略微調(diào)整,保證容量優(yōu)化配置有解且符合實際情況.

    表1 不同情形風光火耦合系統(tǒng)容量優(yōu)化配置結果

    在峰谷差率和新能源利用率相等時,最大負荷減小會導致風光裝機容量和總綜合收益減小,且每減小100 MW最大負荷可導致新能源滲透率減少4%左右,如情形1、7、11和5、9,這是因為負荷整體減小后,更少的風光裝機容量即可滿足負荷需求.

    在最大負荷和新能源利用率相等時,峰谷差率減小可導致風光裝機容量增大,總綜合收益增大,如情形1、5、7、9和3、6、8、10.這是由于負荷需求曲線的峰谷差率減小,使火電可投入更多的調(diào)節(jié)資源給風光,而風光在獲得更多調(diào)節(jié)資源后,能增加裝機容量使總收益變大.雖然負荷峰谷差率受區(qū)域負荷特性影響無法大幅調(diào)整,但可略微減小以提升耦合系統(tǒng)綜合收益.

    在最大負荷和峰谷差率相等時,新能源利用率減小可導致風光總裝機容量和總收益變大,如情形1~4和情形7~8.這是由于新能源最小利用率的降低表示耦合系統(tǒng)可棄更多風光發(fā)電負荷.由此,在一些風光出力較極端且耦合系統(tǒng)的新能源最小利用率較低的場景下,通過棄風棄光即可滿足功率平衡需求;此外,更多的棄風棄光意味著火電投入的調(diào)節(jié)資源更少.上述因素綜合作用導致風光裝機容量及綜合收益增大,仍需指出,新能源最小利用率減小雖利于提升綜合收益,但棄風棄光過多與保證新能源最小利用率的政策不符.因此耦合系統(tǒng)容量優(yōu)化配置中,應考慮新能源最小利用率達一定標準.

    計算不同情形下風電、光伏的裝機容量比可知,風電裝機容量大于光伏裝機容量,新能源最小利用率對風光裝機容量比影響較小,如情形1~4,這是由于新能源中風電在全時間段中具有更好的發(fā)電能力,且爬坡事件發(fā)生較少,新能源最小利用率主要與運行情況有關;最大負荷相同時,負荷峰谷差率減小可導致風光裝機容量比增大,如情形1和5,這是由于負荷峰谷差率減小可使火電具有更多的靈活性資源平抑新能源波動,而新能源中風電在全時間段中具有更好的發(fā)電能力;負荷峰谷差率相同時,最大負荷減小可導致風光裝機容量比減小,如情形1、7、11,這是由于最大負荷減小后,新能源裝機容量需求減小,而風電裝機成本較光伏高.

    綜合上述情形,本文研究區(qū)域的耦合系統(tǒng)中,考慮最大負荷略大于火電機組總裝機容量的情況,新能源的最佳滲透率主要在30%~40%左右,其中風電裝機容量大于光伏裝機容量,且隨著新能源滲透率增大,風光裝機容量比增加;當最大負荷較大、峰谷差率和新能源最小利用率均較小時,耦合系統(tǒng)的新能源滲透率可達到40%以上,如情形6所示.

    3.3 爬坡事件的影響因素分析

    采用1.3節(jié)中最優(yōu)爬坡區(qū)間計算方法,分析耦合系統(tǒng)中風、光年爬坡事件,風光耦合后爬坡事件的影響因素.

    首先分析不同爬坡閾值對上述爬坡事件發(fā)生時間的影響,設置最大負荷為 1 400 MW,峰谷差率為0.55,新能源最小利用率為0.9,爬坡閾值設置為30~120 MW,間隔1 MW,分別計算風、光和風光耦合后年爬坡事件總發(fā)生時間,結果如圖4所示.由圖可知,在相同爬坡閾值下,光伏比風電的年爬坡事件總發(fā)生時間高出很多,結合圖2中風光期望出力曲線可知,該地區(qū)風電整體波動較小,而光伏在每天日出和日落時功率變化很大,容易發(fā)生爬坡事件.爬坡閾值較小和較大,即30~50 MW和75~120 MW時,風光耦合后的年爬坡事件總發(fā)生時間小于光伏的年爬坡事件總發(fā)生時間;爬坡閾值處于中等值,即50~75 MW時,風光耦合后的年爬坡事件總發(fā)生時間大于光伏的年爬坡事件總發(fā)生時間;在本節(jié)研究的爬坡閾值下,風光耦合后的年爬坡事件總發(fā)生時間幾乎均小于風電、光伏的年爬坡事件總發(fā)生時間直接相加.上述現(xiàn)象說明該地區(qū)的風光互補性較好,風電可以減少光伏爬坡事件,尤其在抑制較大或者較小的爬坡事件時效果較好.

    圖4 風電、光伏及其耦合后年爬坡事件總發(fā)生時間與爬坡閾值的關系

    其次分析給定爬坡閾值下,不同容量優(yōu)化配置結果對風、光年爬坡事件總發(fā)生時間和風光耦合后年爬坡事件總發(fā)生時間的影響,爬坡閾值設為80 MW,結果如圖5所示,圖中10組數(shù)據(jù)從左至右分別對應情形1~10.分析可得,雖然新能源裝機容量大會引起其波動變大,更易達到給定爬坡閾值,但年爬坡事件總發(fā)生時間與風光裝機容量并非完全正相關,如情形4和5.這是因為最大負荷、峰谷差率較大,新能源利用率較小時,風電、光伏會為了追求更多收益進行更多棄風、棄光,從而引起更大波動,導致爬坡事件增加.綜上,耦合系統(tǒng)的新能源最小利用率不應過小,否則會導致出現(xiàn)更多爬坡事件,增加火電機組的調(diào)峰壓力,并可能進一步引起其他問題.

    圖5 80 MW爬坡閾值下不同情形風電、光伏及其耦合后爬坡事件影響分析

    綜上所述,該地區(qū)風電是較為優(yōu)質(zhì)的新能源,其爬坡事件發(fā)生概率較小,同時與光伏相比全天均具備發(fā)電能力,因此風電可較光伏配置更多裝機容量.

    3.4 源荷匹配性的影響因素分析

    根據(jù)容量優(yōu)化配置結果計算源荷匹配評價指標并分析其影響因素,結果如圖6所示.由圖可見,最大負荷減小,源荷追蹤系數(shù)增大,而源荷變化系數(shù)、源荷平均絕對值變化和源荷動態(tài)時間彎曲距離均減小,如情形1、7、11,這是因為最大負荷減小導致負荷需求曲線整體波動和變化變小,同時風光出力可更好地追蹤負荷需求,并更好地與負荷需求曲線保持一致;峰谷差率減小可導致源荷追蹤系數(shù)、源荷平均絕對值變化和源荷動態(tài)時間彎曲距離減小,源荷變化系數(shù)增大,如情形3、6,這是由于峰谷差率減小時的負荷需求曲線更平緩,但負荷需求曲線的平緩會導致新能源出力曲線與其更難貼合和近似;新能源最小利用率減小會導致4個源荷匹配指標增大,如情形1~4,這是由于新能源最小利用率減小后,風光可靈活棄風、棄光以更好地貼近負荷曲線并保持形狀相似.這有利于源荷匹配且可減小其匹配后的波動,但因新能源最小利用率減小導致風光裝機容量增大,在裝機容量增大的同時為保證耦合系統(tǒng)利益最大化,會在風光充足時優(yōu)先出力,上述原因共同作用使源荷變化系數(shù)、源荷平均絕對值變化和源荷動態(tài)時間彎曲距離略微提升.

    圖6 源荷匹配評價指標分析

    從源荷匹配評價可以看出,在源荷匹配較好時,適當降低可再生能源最小利用率可進一步增加風、光裝機容量,即盡可能合適的負荷需求曲線制定對耦合系統(tǒng)的收益及新能源的充分利用具有較大作用.

    4 結語

    研究風光火耦合系統(tǒng)的容量優(yōu)化配置,針對已有火電機組,考慮多種源荷特性、成本和收益,配置風電、光伏的裝機容量,并分析多種源荷因素與容量優(yōu)化配置結果之間的相互影響.以遼寧地區(qū)含兩臺600 MW 火電機組的風光火耦合系統(tǒng)容量優(yōu)化配置為算例進行仿真,結果表明:減小新能源最小利用率、峰谷差率和發(fā)電缺額率可增大裝機容量;該地區(qū)風電發(fā)生爬坡事件的概率較光伏小,風、光互補后,可較好地平抑較大和較小的爬坡事件;新能源最小利用率小會導致更多爬坡事件發(fā)生,增加火電的調(diào)峰壓力;最大負荷和峰谷差率對源荷匹配性有相反的影響,減小新能源最小利用率有利于更好的源荷匹配.綜上所述,本文研究區(qū)域的耦合系統(tǒng)中,考慮最大負荷略大于火電機組總裝機容量的情況,新能源的最佳滲透率在30%~40%左右;其中風電可較光伏配置更多裝機容量,以實現(xiàn)更好的爬坡事件平抑及源荷匹配;若源荷匹配較好,適當減小新能源最小利用率可進一步增加新能源的裝機容量.

    猜你喜歡
    爬坡火電風光
    陜西農(nóng)民收入:一路爬坡過坎
    當代陜西(2021年13期)2021-08-06 09:24:32
    風光新580
    汽車觀察(2021年11期)2021-04-24 20:47:38
    爬坡過坎
    當代工人(2020年1期)2020-05-11 11:47:32
    風光如畫
    海峽姐妹(2019年12期)2020-01-14 03:25:02
    風光ix5:當轎跑邂逅SUV
    汽車觀察(2018年12期)2018-12-26 01:05:36
    火電施工EPC項目管理探討
    我國防腐木進口又現(xiàn)平穩(wěn)“爬坡”
    各國首都風光
    向下的火電
    能源(2015年8期)2015-05-26 09:15:36
    火電脫硝“大限”將至
    自動化博覽(2014年4期)2014-02-28 22:31:18
    成人特级av手机在线观看| 精品久久久久久久久av| av卡一久久| 最近的中文字幕免费完整| 尤物成人国产欧美一区二区三区| 日韩大尺度精品在线看网址| 欧美不卡视频在线免费观看| 日本黄大片高清| avwww免费| 99九九线精品视频在线观看视频| 能在线免费观看的黄片| 精品一区二区三区av网在线观看| 黄色配什么色好看| 欧美日韩在线观看h| 亚洲精品粉嫩美女一区| 禁无遮挡网站| 老师上课跳d突然被开到最大视频| 不卡视频在线观看欧美| 国产三级中文精品| 成年女人永久免费观看视频| 在现免费观看毛片| 中文在线观看免费www的网站| 看免费成人av毛片| 午夜福利18| 亚洲av.av天堂| 日本爱情动作片www.在线观看 | 久久国产乱子免费精品| 97在线视频观看| 亚洲最大成人手机在线| 国产色爽女视频免费观看| 看十八女毛片水多多多| 少妇人妻一区二区三区视频| 欧美性感艳星| 97超级碰碰碰精品色视频在线观看| 99精品在免费线老司机午夜| 国产女主播在线喷水免费视频网站 | 亚洲精品色激情综合| 69av精品久久久久久| 色视频www国产| 精品久久久久久久人妻蜜臀av| 不卡一级毛片| 免费高清视频大片| 麻豆国产av国片精品| 国产 一区 欧美 日韩| 少妇熟女aⅴ在线视频| 插逼视频在线观看| 婷婷六月久久综合丁香| 性插视频无遮挡在线免费观看| 极品教师在线视频| 国产午夜精品论理片| 欧美+亚洲+日韩+国产| 亚洲熟妇中文字幕五十中出| 18禁裸乳无遮挡免费网站照片| 午夜福利成人在线免费观看| 少妇熟女aⅴ在线视频| 国产男靠女视频免费网站| 国产私拍福利视频在线观看| 中国美女看黄片| 如何舔出高潮| 看免费成人av毛片| 久久婷婷人人爽人人干人人爱| av在线播放精品| 99热这里只有是精品在线观看| 国产伦精品一区二区三区四那| 身体一侧抽搐| 国产成年人精品一区二区| 久久午夜福利片| 国产久久久一区二区三区| 91狼人影院| 日日啪夜夜撸| 久久久久久大精品| 少妇熟女欧美另类| 99精品在免费线老司机午夜| 国产三级中文精品| 在线观看美女被高潮喷水网站| 99热全是精品| 麻豆一二三区av精品| 久久精品人妻少妇| av在线天堂中文字幕| 最近的中文字幕免费完整| 麻豆一二三区av精品| 嫩草影视91久久| 一级毛片电影观看 | 美女黄网站色视频| 亚洲四区av| 亚洲av成人精品一区久久| 国产高清不卡午夜福利| 亚州av有码| 免费人成视频x8x8入口观看| 在线免费十八禁| 精品日产1卡2卡| 日韩大尺度精品在线看网址| 色视频www国产| 性欧美人与动物交配| 欧洲精品卡2卡3卡4卡5卡区| 午夜日韩欧美国产| 你懂的网址亚洲精品在线观看 | 午夜激情福利司机影院| 一个人观看的视频www高清免费观看| 一级毛片久久久久久久久女| 日本a在线网址| 少妇猛男粗大的猛烈进出视频 | 可以在线观看的亚洲视频| 亚洲精品久久国产高清桃花| 好男人在线观看高清免费视频| 久久精品国产99精品国产亚洲性色| 晚上一个人看的免费电影| 午夜精品国产一区二区电影 | 老熟妇仑乱视频hdxx| 少妇高潮的动态图| 精品熟女少妇av免费看| 亚洲真实伦在线观看| 欧美bdsm另类| 国产一区亚洲一区在线观看| 97在线视频观看| 欧美激情在线99| 亚洲欧美日韩高清专用| 少妇熟女欧美另类| 国产乱人偷精品视频| 18+在线观看网站| 又黄又爽又免费观看的视频| 国产亚洲精品久久久久久毛片| 99久久精品热视频| 毛片女人毛片| 久久亚洲国产成人精品v| 亚洲精品日韩av片在线观看| videossex国产| 又黄又爽又刺激的免费视频.| 全区人妻精品视频| 免费高清视频大片| 干丝袜人妻中文字幕| 日日摸夜夜添夜夜添av毛片| 久久精品国产自在天天线| 女人十人毛片免费观看3o分钟| 一卡2卡三卡四卡精品乱码亚洲| 白带黄色成豆腐渣| 久久中文看片网| 毛片女人毛片| 国产亚洲精品综合一区在线观看| 亚洲国产精品sss在线观看| 日日干狠狠操夜夜爽| 欧美高清性xxxxhd video| 黄色一级大片看看| 深夜精品福利| 中文字幕精品亚洲无线码一区| 欧美人与善性xxx| 搡老岳熟女国产| 搡老妇女老女人老熟妇| 日韩欧美国产在线观看| 久久热精品热| 男女之事视频高清在线观看| 精品人妻一区二区三区麻豆 | 国产一区二区三区在线臀色熟女| 美女高潮的动态| 看免费成人av毛片| 国产一区亚洲一区在线观看| 桃色一区二区三区在线观看| 大香蕉久久网| 日本与韩国留学比较| 免费人成视频x8x8入口观看| 丝袜喷水一区| 一区二区三区高清视频在线| 丰满乱子伦码专区| 亚洲av不卡在线观看| 国产精品亚洲一级av第二区| 亚洲18禁久久av| 可以在线观看的亚洲视频| 伦精品一区二区三区| 日韩制服骚丝袜av| 最近在线观看免费完整版| 少妇的逼水好多| 国产午夜精品久久久久久一区二区三区 | 乱系列少妇在线播放| 中文资源天堂在线| 草草在线视频免费看| 国内精品宾馆在线| 少妇的逼水好多| 精品久久久久久久久亚洲| 欧美成人一区二区免费高清观看| 亚洲人成网站在线观看播放| 成人三级黄色视频| 免费观看的影片在线观看| 天堂av国产一区二区熟女人妻| 国产三级在线视频| 最近手机中文字幕大全| 亚洲性久久影院| 日本 av在线| 国产亚洲精品久久久com| 波多野结衣巨乳人妻| 精品久久久久久久久亚洲| 亚洲四区av| 欧美日韩国产亚洲二区| 国产私拍福利视频在线观看| 成年免费大片在线观看| 欧美日韩综合久久久久久| 亚洲精品粉嫩美女一区| 亚洲美女搞黄在线观看 | 免费看av在线观看网站| 日本黄色视频三级网站网址| 欧美性猛交╳xxx乱大交人| 99热这里只有精品一区| 国产免费一级a男人的天堂| 舔av片在线| 久久久久精品国产欧美久久久| 91久久精品电影网| 女人十人毛片免费观看3o分钟| 精品一区二区三区av网在线观看| 97人妻精品一区二区三区麻豆| 亚洲婷婷狠狠爱综合网| 午夜精品一区二区三区免费看| 亚洲中文字幕一区二区三区有码在线看| 听说在线观看完整版免费高清| 亚洲久久久久久中文字幕| 12—13女人毛片做爰片一| 美女xxoo啪啪120秒动态图| 久久精品国产亚洲av天美| 日本黄大片高清| 少妇裸体淫交视频免费看高清| 亚洲人成网站在线播| 人妻夜夜爽99麻豆av| 美女xxoo啪啪120秒动态图| 尤物成人国产欧美一区二区三区| 国产三级中文精品| 午夜免费激情av| 国产精品亚洲美女久久久| 99久久中文字幕三级久久日本| 少妇人妻精品综合一区二区 | 99久久九九国产精品国产免费| АⅤ资源中文在线天堂| 校园春色视频在线观看| 又粗又爽又猛毛片免费看| 久久天躁狠狠躁夜夜2o2o| 波多野结衣巨乳人妻| 久久精品国产亚洲av香蕉五月| 99久国产av精品国产电影| 日日摸夜夜添夜夜添av毛片| 亚洲熟妇中文字幕五十中出| 国产精品永久免费网站| 91午夜精品亚洲一区二区三区| 毛片一级片免费看久久久久| 尾随美女入室| 亚洲国产欧洲综合997久久,| 高清毛片免费看| 国产爱豆传媒在线观看| 天天躁日日操中文字幕| 欧美日本视频| 国产男人的电影天堂91| 亚洲精品乱码久久久v下载方式| 亚洲激情五月婷婷啪啪| 男人的好看免费观看在线视频| 97超碰精品成人国产| 久久久久久久久大av| 色视频www国产| 人妻夜夜爽99麻豆av| 女人十人毛片免费观看3o分钟| 欧美潮喷喷水| 国产真实伦视频高清在线观看| 久久久午夜欧美精品| 一级毛片电影观看 | 日本一本二区三区精品| 变态另类成人亚洲欧美熟女| 精品人妻视频免费看| 亚洲熟妇熟女久久| 别揉我奶头~嗯~啊~动态视频| 欧洲精品卡2卡3卡4卡5卡区| 少妇熟女aⅴ在线视频| 欧美日韩精品成人综合77777| а√天堂www在线а√下载| 欧美成人免费av一区二区三区| 久久久久国产网址| 无遮挡黄片免费观看| 亚洲七黄色美女视频| 欧美成人免费av一区二区三区| 伦理电影大哥的女人| 亚洲美女视频黄频| 亚洲欧美清纯卡通| 人妻丰满熟妇av一区二区三区| 日韩人妻高清精品专区| 一级毛片久久久久久久久女| 成年女人毛片免费观看观看9| 内射极品少妇av片p| 欧美日韩在线观看h| 中国美女看黄片| 插阴视频在线观看视频| 久久精品夜色国产| 美女被艹到高潮喷水动态| 香蕉av资源在线| 亚洲国产精品成人综合色| 晚上一个人看的免费电影| 亚洲av中文字字幕乱码综合| 精品熟女少妇av免费看| 免费高清视频大片| 最近2019中文字幕mv第一页| av在线亚洲专区| 日日撸夜夜添| 日韩高清综合在线| 啦啦啦观看免费观看视频高清| 国内少妇人妻偷人精品xxx网站| 听说在线观看完整版免费高清| 免费看日本二区| 老熟妇乱子伦视频在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 国产亚洲精品av在线| 91av网一区二区| 午夜福利在线观看吧| 小蜜桃在线观看免费完整版高清| 欧美bdsm另类| 国产精品久久视频播放| 99热精品在线国产| 人妻夜夜爽99麻豆av| 波多野结衣高清作品| 精品人妻熟女av久视频| 少妇的逼水好多| 在线观看av片永久免费下载| 国产精品人妻久久久影院| 久99久视频精品免费| 在线国产一区二区在线| 村上凉子中文字幕在线| 久久久久九九精品影院| 久久久精品大字幕| 两个人的视频大全免费| 久久久久久大精品| 亚洲欧美日韩高清专用| 国产精品野战在线观看| 日日啪夜夜撸| 深爱激情五月婷婷| 亚洲美女搞黄在线观看 | 日韩欧美免费精品| 在线播放国产精品三级| 夜夜看夜夜爽夜夜摸| АⅤ资源中文在线天堂| 麻豆av噜噜一区二区三区| 婷婷亚洲欧美| 亚洲18禁久久av| 麻豆成人午夜福利视频| 国产亚洲精品综合一区在线观看| 男女下面进入的视频免费午夜| 亚洲av一区综合| 亚洲天堂国产精品一区在线| 男女啪啪激烈高潮av片| 少妇裸体淫交视频免费看高清| 欧美成人免费av一区二区三区| 国产成人精品久久久久久| 男女视频在线观看网站免费| 亚洲精品一卡2卡三卡4卡5卡| ponron亚洲| 国产精品久久视频播放| 一级黄色大片毛片| 六月丁香七月| 国产视频一区二区在线看| 男女视频在线观看网站免费| 亚洲经典国产精华液单| 免费搜索国产男女视频| 婷婷精品国产亚洲av| 日本爱情动作片www.在线观看 | av视频在线观看入口| 久久久久久久久久成人| 深爱激情五月婷婷| 熟妇人妻久久中文字幕3abv| 男女之事视频高清在线观看| 久久久久国产网址| 在线观看美女被高潮喷水网站| 91在线观看av| 一个人免费在线观看电影| 色5月婷婷丁香| 久久久久久国产a免费观看| 联通29元200g的流量卡| 国产一区二区亚洲精品在线观看| 淫妇啪啪啪对白视频| 亚洲精品成人久久久久久| 国产亚洲欧美98| 日韩精品青青久久久久久| 国产黄色视频一区二区在线观看 | 久久久久国产精品人妻aⅴ院| 日韩成人av中文字幕在线观看 | 久久韩国三级中文字幕| 亚洲av五月六月丁香网| 欧美xxxx黑人xx丫x性爽| 色综合色国产| 久久久久国内视频| 精品熟女少妇av免费看| 国产亚洲精品久久久com| 免费观看在线日韩| 精品一区二区三区视频在线| 熟女电影av网| 99久国产av精品国产电影| 精品一区二区三区视频在线观看免费| 日本欧美国产在线视频| 国产高清视频在线播放一区| 最近的中文字幕免费完整| 亚洲av中文av极速乱| 非洲黑人性xxxx精品又粗又长| 22中文网久久字幕| 一进一出抽搐动态| av在线观看视频网站免费| 三级毛片av免费| av天堂在线播放| 秋霞在线观看毛片| 亚洲真实伦在线观看| 桃色一区二区三区在线观看| 村上凉子中文字幕在线| 亚洲av不卡在线观看| 久久精品久久久久久噜噜老黄 | 久99久视频精品免费| 午夜福利高清视频| 麻豆精品久久久久久蜜桃| 午夜影院日韩av| 色视频www国产| 国产真实伦视频高清在线观看| 老司机午夜福利在线观看视频| 一级a爱片免费观看的视频| 老司机福利观看| 国产一区二区在线观看日韩| 欧美国产日韩亚洲一区| 国产精品一及| 精品一区二区三区av网在线观看| 色哟哟·www| 亚洲欧美成人精品一区二区| 午夜亚洲福利在线播放| 国语自产精品视频在线第100页| 丝袜美腿在线中文| 国产精品久久久久久久电影| 亚洲国产色片| 蜜臀久久99精品久久宅男| 桃色一区二区三区在线观看| 少妇被粗大猛烈的视频| 老女人水多毛片| 蜜桃亚洲精品一区二区三区| 美女大奶头视频| 国产91av在线免费观看| 免费av毛片视频| av黄色大香蕉| 亚洲内射少妇av| 国产精品不卡视频一区二区| 日本免费a在线| 日本熟妇午夜| 在线观看午夜福利视频| 麻豆久久精品国产亚洲av| 国产精品亚洲美女久久久| 精品乱码久久久久久99久播| 国产色婷婷99| 哪里可以看免费的av片| 久久久久免费精品人妻一区二区| 日韩av在线大香蕉| 国产老妇女一区| 麻豆国产av国片精品| 日韩精品青青久久久久久| 狂野欧美白嫩少妇大欣赏| 日本五十路高清| 亚洲精品456在线播放app| av国产免费在线观看| 成年av动漫网址| 久久精品国产清高在天天线| 国产探花在线观看一区二区| 我要看日韩黄色一级片| 欧美激情久久久久久爽电影| 成人无遮挡网站| 三级国产精品欧美在线观看| 麻豆乱淫一区二区| 国产一区二区在线观看日韩| 国产色爽女视频免费观看| 国产黄a三级三级三级人| 国产精品不卡视频一区二区| 俄罗斯特黄特色一大片| 少妇的逼水好多| 床上黄色一级片| 亚洲欧美日韩高清在线视频| avwww免费| 99热全是精品| 成人特级黄色片久久久久久久| 日日干狠狠操夜夜爽| 国产精品99久久久久久久久| 床上黄色一级片| 国产aⅴ精品一区二区三区波| 少妇被粗大猛烈的视频| 免费看av在线观看网站| 国产精品一区二区性色av| 国产精品久久视频播放| 大香蕉久久网| 日本免费a在线| 国产探花在线观看一区二区| 国产精品一及| 久久精品国产鲁丝片午夜精品| 97人妻精品一区二区三区麻豆| 麻豆精品久久久久久蜜桃| 日韩一本色道免费dvd| 在现免费观看毛片| 日本-黄色视频高清免费观看| 一级毛片我不卡| 亚洲成人av在线免费| 18禁在线无遮挡免费观看视频 | 69人妻影院| 免费看光身美女| 一个人免费在线观看电影| 欧美人与善性xxx| 欧美一区二区亚洲| 亚洲国产日韩欧美精品在线观看| 91久久精品国产一区二区三区| 青春草视频在线免费观看| 亚洲国产欧洲综合997久久,| 最近最新中文字幕大全电影3| 久久精品国产清高在天天线| 国内久久婷婷六月综合欲色啪| 婷婷色综合大香蕉| 国产精品伦人一区二区| 久久99热这里只有精品18| 日日摸夜夜添夜夜添小说| АⅤ资源中文在线天堂| 99精品在免费线老司机午夜| a级毛片a级免费在线| 成年版毛片免费区| 国模一区二区三区四区视频| 熟女电影av网| 一本一本综合久久| 久久国产乱子免费精品| 精品一区二区三区av网在线观看| 观看美女的网站| 在线a可以看的网站| av女优亚洲男人天堂| 日本色播在线视频| 天堂网av新在线| 观看美女的网站| 免费观看在线日韩| 欧美丝袜亚洲另类| 国内揄拍国产精品人妻在线| 老女人水多毛片| 国产色婷婷99| 校园人妻丝袜中文字幕| 亚洲中文字幕日韩| 啦啦啦啦在线视频资源| 97超视频在线观看视频| 亚洲欧美成人精品一区二区| 一边摸一边抽搐一进一小说| 日日啪夜夜撸| 午夜久久久久精精品| 亚洲熟妇熟女久久| 色综合色国产| 一进一出抽搐动态| 99久久精品热视频| 看黄色毛片网站| 卡戴珊不雅视频在线播放| 久久精品91蜜桃| 免费在线观看成人毛片| 国产成人91sexporn| 看片在线看免费视频| 最近中文字幕高清免费大全6| 不卡一级毛片| 成熟少妇高潮喷水视频| 精品久久国产蜜桃| 婷婷精品国产亚洲av| 久久精品夜色国产| 一进一出好大好爽视频| 亚洲无线观看免费| 中国美女看黄片| 99国产极品粉嫩在线观看| 国产在视频线在精品| 成人特级av手机在线观看| .国产精品久久| 国产成人aa在线观看| 少妇丰满av| 看片在线看免费视频| 美女内射精品一级片tv| 亚洲激情五月婷婷啪啪| 尤物成人国产欧美一区二区三区| 俄罗斯特黄特色一大片| 永久网站在线| 午夜福利成人在线免费观看| 国产精品日韩av在线免费观看| 黄色日韩在线| 国产日本99.免费观看| 亚洲综合色惰| 国产精品亚洲美女久久久| 久久久久九九精品影院| 亚洲精品日韩av片在线观看| 亚洲三级黄色毛片| 夜夜爽天天搞| 成人毛片a级毛片在线播放| 精品久久久久久成人av| 国产高清视频在线观看网站| 在现免费观看毛片| 成人永久免费在线观看视频| 俄罗斯特黄特色一大片| 国产亚洲精品久久久久久毛片| 性欧美人与动物交配| 长腿黑丝高跟| 亚洲av五月六月丁香网| 免费在线观看成人毛片| 深夜精品福利| 亚洲国产精品合色在线| 99久国产av精品国产电影| 看非洲黑人一级黄片| 看黄色毛片网站| 亚洲欧美清纯卡通| 99视频精品全部免费 在线| 国产日本99.免费观看| 日韩三级伦理在线观看| 伊人久久精品亚洲午夜| 成人无遮挡网站| 少妇人妻一区二区三区视频| 可以在线观看毛片的网站| videossex国产| 精品熟女少妇av免费看| 精品99又大又爽又粗少妇毛片| 国国产精品蜜臀av免费| 看片在线看免费视频| 欧美成人a在线观看| 午夜影院日韩av| 国产精品电影一区二区三区| 乱人视频在线观看| 人妻夜夜爽99麻豆av| 色尼玛亚洲综合影院| 黄色视频,在线免费观看| 国产白丝娇喘喷水9色精品| 成年免费大片在线观看| 亚洲熟妇熟女久久|