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

    考慮分時電價下需求響應(yīng)的用戶側(cè)儲能優(yōu)化配置

    2024-01-01 00:00:00陳景單茜王金鋒孫曉晨姜炎君
    陜西科技大學(xué)學(xué)報 2024年4期

    摘 要:針對用戶側(cè)儲能在投資與運行方面經(jīng)濟性不足的問題,提出一種考慮分時電價下需求響應(yīng)的用戶側(cè)儲能規(guī)劃與運行兩階段優(yōu)化策略.第一階段通過建立分時電價優(yōu)化模型,采用粒子群(PSO)和序列二次規(guī)劃(SQP)混合算法求解得出分時電價策略,第二階段基于第一階段得出的分時電價策略建立考慮兩個時間尺度的儲能雙層優(yōu)化模型,上層模型以最大化全壽命周期凈收益為目標(biāo)優(yōu)化長時間尺度的儲能配置,下層模型以最大化日運行收益為目標(biāo)優(yōu)化短時間尺度的儲能運行,結(jié)合混合整數(shù)線性規(guī)劃方法和KKT條件,調(diào)用CPLEX求解器進行求解,通過兩階段的聯(lián)合優(yōu)化,實現(xiàn)用戶側(cè)儲能配置和運行收益綜合最優(yōu).最后,結(jié)合實例進行仿真,結(jié)果表明所提方法可以有效提高用戶側(cè)儲能系統(tǒng)的經(jīng)濟性.

    關(guān)鍵詞:用戶側(cè)儲能; 分時電價; 需求響應(yīng); 兩階段優(yōu)化

    中圖分類號:TM721

    文獻標(biāo)志碼: A

    Optimal configuration of user-side energy storage considering demand response under time of use price

    CHEN Jing-wen1*, SHAN Xi1, WANG Jin-feng2, SUN Xiao-chen2, JIANG Yan-jun2

    (1.School of Electrical and Control Engineering, Shaanxi University of Science amp; Technology, Xi′an 710021, China;" 2.Economic Research Institute, State Grid Shaanxi Electric Power Company Limited, Xi′an 710065, China)

    Abstract:In view of the lack of economic benefits in investment and operation of user-side energy storage,a two-stage optimal strategy for user-side energy storage planning and operation considering demand response under time of use (TOU) price is proposed.In the first stage,the optimization model of TOU price is established,and the optimization of TOU price strategy is obtained by combining particle swarm optimization (PSO) and sequential quadratic programming (SQP) to guide users to use electricity and optimize the load curve.In the second stage,based on the TOU price strategy obtained in the first stage,a two-level optimal model of energy storage is established,which considers two time scales,the upper model optimizes the long-time storage allocation with the goal of maximizing the life-cycle net benefit,while the lower model optimizes the short-time storage operation with the goal of maximizing the daily operating benefit,combining mixed integer linear programming (MILP) method and KKT condition,and calling CPLEX solver to solve,the optimal configuration and operation benefit of the user-side energy storage can be achieved by the joint optimization of two-stage model.Ultimately,the simulation results demonstrate the efficacy of the proposed approach in enhancing the economic performance of the user-side energy storage system.

    Key words:user-side energy storage; time-of-use price; demand response; two-stage optimization

    0 引言

    在我國積極建設(shè)新型電力系統(tǒng)的進程中,各地區(qū)不斷完善電價政策,為用戶側(cè)儲能提供了廣闊前景[1,2].用戶側(cè)儲能主要依賴峰谷價差進行套利,多地政府已經(jīng)發(fā)布文件支持通過拉大峰谷電價差提升儲能應(yīng)用的比例.用戶側(cè)儲能可以節(jié)約用電成本,提升供電可靠性和用電質(zhì)量,還可以參與電網(wǎng)削峰填谷,平抑新能源并網(wǎng)波動性和需求響應(yīng)等服務(wù).然而,高昂的儲能電池成本很大程度上限制了用戶側(cè)儲能的應(yīng)用和發(fā)展[3].因此,在新型電力系統(tǒng)中,合理規(guī)劃用戶側(cè)儲能與需求響應(yīng)負荷等分布式資源至關(guān)重要.

    在電力市場競爭中,電價信號和激勵機制可以有效激勵用戶參與需求響應(yīng)[4,5],分時電價(Time of Use,TOU)是一種價格型需求響應(yīng)[6],通過設(shè)定不同時間段的電價,用戶能夠靈活調(diào)整用電行為.分時電價不僅能提高電力市場的運行效率和供需平衡,還能鼓勵用戶節(jié)約能源和優(yōu)化負荷曲線.因此,分時電價作為電力市場中的重要舉措,具有巨大的潛力和價值.文獻[7]基于支持向量機的模型,分析了分時電價與用戶響應(yīng)之間的相互作用.文獻[8]通過優(yōu)化分時電價對電網(wǎng)側(cè)變電站容量進行規(guī)劃,但在成本建模中未充分考慮供電方的利益損失.文獻[9]通過綜合考慮供電方和用電方的成本效益問題,開發(fā)了一種基于主從博弈模型的分時電價優(yōu)化策略,最大限度地實現(xiàn)供需平衡,降低電力系統(tǒng)的運營成本和用戶的用電成本.文獻[10]建立了一種基于改進帝國主義競爭算法的分時電價優(yōu)化模型,通過優(yōu)化分時電價,提高配電網(wǎng)光伏接納能力.上述文獻中,分時電價下需求響應(yīng)均沒有考慮用戶側(cè)儲能的作用.

    儲能優(yōu)化主要集中在電網(wǎng)側(cè)和用戶側(cè),利用儲能“低儲高發(fā)”的特性獲取電費利潤,改善新能源并網(wǎng)帶來的隨機性.文獻[11,12]從電網(wǎng)側(cè)出發(fā),通過考慮新能源接入電網(wǎng)的不確定性和儲能參與市場輔助服務(wù),對儲能的容量配置和充放電策略進行優(yōu)化.對于用戶側(cè)儲能,文獻[13]對用戶側(cè)儲能建立了全壽命周期模型,對其進行了容量配置,但沒有涉及儲能運行策略.文獻[14]基于層次分析法建立了儲能配置評估模型,針對適合配置儲能系統(tǒng)的用戶進行儲能運行優(yōu)化,未提出儲能配置的方案.文獻[15,16]同時考慮用戶側(cè)儲能的配置和運行問題,通過考慮兩種不同時間尺度下的儲能收益,建立了用戶側(cè)儲能規(guī)劃和調(diào)度綜合優(yōu)化模型.文獻[17]將需求響應(yīng)與儲能優(yōu)化同時考慮,以負荷聚合商為研究對象分析其運營模式,通過配置儲能系統(tǒng)提高需求響應(yīng)的可靠性,但未充分考慮儲能系統(tǒng)的經(jīng)濟收益問題.文獻[18]建立了考慮需求響應(yīng)的電熱氣云儲能配置和運行優(yōu)化模型,但未充分研究分時電價下需求響應(yīng)機制.上述研究對用戶側(cè)儲能優(yōu)化進行了研究,但未與需求響應(yīng)策略較好地相結(jié)合.

    綜上所述,分時電價策略和儲能配置的單獨研究較為成熟,但目前結(jié)合需求響應(yīng)以提升用戶側(cè)儲能經(jīng)濟性研究仍不充分.本文從該角度出發(fā)提出了一種考慮分時電價下需求響應(yīng)的用戶側(cè)儲能規(guī)劃和運行兩階段優(yōu)化策略.在第一階段,以兼顧負荷特性和用戶滿意度為目標(biāo)建立分時電價優(yōu)化模型,結(jié)合粒子群算法(Particle Swarm Optimization,PSO)和序列二次規(guī)劃算法(Sequential Quadratic Programming,SQP)進行尋優(yōu)求解.在此基礎(chǔ)上;在第二階段建立考慮不同時間尺度的配置運行聯(lián)合優(yōu)化雙層模型,基于KKT條件將雙層模型轉(zhuǎn)換為單層模型,結(jié)合混合整數(shù)線性規(guī)劃方法并調(diào)用CPLEX求解器進行求解.最后,通過算例仿真分析,驗證所提模型和策略的有效性.

    1 分時電價優(yōu)化

    分時電價優(yōu)化包括需求響應(yīng)建模、時段劃分和電價制定三個方面.

    1.1 分時電價下的需求響應(yīng)

    電力負荷彈性系數(shù)是指負荷對電價變化的敏感程度,即當(dāng)電價變化時用電量的相對變化率.用戶在i時段的負荷和j時段的價格之間的彈性系數(shù)定義為:

    e(i,j)=ΔL(i)Δρ(j)ρ0(j)L0(i)

    (1)

    式(1)中:e為電力負荷彈性系數(shù);ΔL為用電量的變化量;Δρ為電價的變化量;ρ0和L0為實施分時電價前的電價和用電量;當(dāng)i=j時,彈性系數(shù)為自彈性系數(shù),表示一個時段負荷變化相對于該時段電價變化的敏感程度;當(dāng)i≠j時,彈性系數(shù)為互彈性系數(shù),表示一個時段負荷變化相對于另一個時段電價變化的敏感程度.

    分時電價將一天的負荷分為谷、平、峰三種時段,各時段的電價可表示為:

    ρv=ρ0+δv,i∈Tv

    ρf=ρ0+δf,i∈Tf

    ρp=ρ0+δp,i∈Tp

    (2)

    式(2)中:ρv、ρf、ρp分別為谷、平、峰時段電價;δv、δf、δp分別為谷、平、峰時段電價波動量;Tv、Tf、Tp分別為谷、平、峰時段.

    各時段電價波動比可表示為:

    kv=δv/ρ0,i∈Tv

    kf=δf/ρ0,i∈Tf

    kp=δp/ρ0,i∈Tp

    (3)

    根據(jù)電力負荷彈性系數(shù)和分時電價環(huán)境下用戶需求響應(yīng)矩陣,可以得出分時電價下的負荷響應(yīng)行為:

    L=L0+λ·L0

    (4)

    其中,L0=Lv

    Lf

    Lp,λ=λvv-λvf-λvp

    λfvλffλfp

    λpvλpfλpp,L=L0v

    L0f

    L0p

    式(4)中:向量L表示實施分時電價后用戶負荷,λ為用戶需求響應(yīng)矩陣[19],其元素的值與負荷的彈性系數(shù)、分時電價的價格波動程度有關(guān),表達式為:

    λmm=kmE(i,i)""" i∈Tm

    λmn=(km-kn)E(i,j) i∈Tm,j∈Tn

    (5)

    式(5)中:m,n為不同時段,m,n∈(v,f,p)且m≠n,E(i,j)為負荷彈性系數(shù)矩陣.本文選用24數(shù)據(jù)點的日負荷曲線,則λ矩陣為24×24矩陣.

    1.2 分時電價策略優(yōu)化

    1.2.1 改進移動邊界方法的時段劃分

    移動邊界方法[20]用負荷均方差評估時段劃分結(jié)果,無法體現(xiàn)各時段的類間與類內(nèi)距離[21],本文提出添加輪廓系數(shù)這一指標(biāo),與均方差共同評估時段劃分效果.均方差計算表達式為:

    c(Bfv,Bpf)=124∑i∈[1,24](Li-L-m)2

    L-m=1Nm∑i∈mLi

    (6)

    式(6)中:Bfv、Bpf分別為平谷、峰平的邊界變量;Li是第i小時負荷;L-m為m時段負荷的平均值;Nm為m時段的小時數(shù).

    負荷輪廓系數(shù)的計算表達式為:

    s(Bfv,Bpf)=13∑sm,m∈(v,f,p)

    sm=1Nm∑i∈mb(i)-a(i)maxa(i),b(i)

    (7)

    式(7)中:sm為m時段的負荷輪廓系數(shù);a(i),b(i)分別表示類內(nèi)緊密度和類間分離度,表達式為:

    a(i)=1Nm-1∑i≠jLi-Lj,i,j∈m

    b(i)=minLi-Lk,i,km

    (8)

    綜上,時段劃分的目標(biāo)函數(shù)為:

    min[Fpp=c·(1-s)]

    (9)

    改進移動邊界方法的時段劃分步驟如下,求解流程如圖1所示.

    步驟1 輸入日負荷序列L=L1,L2,...,L24,對其升序排序,得到新的序列L′=L′1,L′2,...,L′24;

    步驟2 初始化邊界變量.令i=1,j=2,Bfv=L′1,Bpf=L′2,計算目標(biāo)函數(shù)的初始值;

    步驟3 更新邊界變量.如果jlt;24,則通過j=j+1更新Bpf,否則通過i=i+1,j=i+1更新邊界變量Bfv,Bpf,然后計算目標(biāo)函數(shù)值Fpp;

    步驟4 迭代停止.如果j=24,則停止迭代;否則,返回步驟3;

    步驟5 輸出結(jié)果.輸出目標(biāo)函數(shù)值Fpp及其對應(yīng)的邊界變量Bfv,Bpf.

    1.2.2 分時電價制定

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

    實施分時電價目的包括削峰填谷和提高用戶滿意度兩個方面[22].降低負荷峰谷差一般選取最小化峰谷差和最小化峰值負荷為目標(biāo)函數(shù),如式(10)所示.

    minf1=Lmax-LminL0max-L0min

    minf2=LmaxL0max

    (10)

    式(10)中:L0min,L0max分別為實施分時電價前日負荷曲線最大,最小值;Lmin,Lmax分別為實施分時電價后日負荷曲線的最大,最小值.

    用戶滿意度是從用電成本和用電習(xí)慣兩個方面綜合評估的結(jié)果[23],會對用戶參與需求響應(yīng)的積極性產(chǎn)生直接影響.用電成本改變引起的用戶滿意度變化可以表示為:

    f3=1-E-E0E0

    E=∑24i=1(Li·ρi)

    E0=∑24i=1(L0i·ρ0i)

    (11)

    式(11)中:E0,E分別為用戶實施分時電價前后的用電成本;i為小時數(shù),ρ0i,ρi分別為實施分時電價前后第i小時的電價;L0i,Li分別為實施分時電價前后第i小時的負荷.f3表示用電成本滿意度,取值等于1表示用戶自由用電,取值大于1表示用戶滿意度有所提高,小于1則相反.

    由于用電習(xí)慣改變引起的用戶滿意度變化可以用轉(zhuǎn)移負荷與初始負荷的比值表示:

    f4=1-ΔLL0

    ΔL=∑24i=1Li-L0i

    L0=∑24i=1L0i

    (12)

    式(12)中:ΔL,L0分別為轉(zhuǎn)移負荷與初始負荷.f4表示用電習(xí)慣滿意度,取值范圍為[0,1],1表示滿意度最高,0表示滿意度最低.

    最大化用戶滿意度的模型可以表示為:

    max(θ1f3+θ2f4)

    (13)

    式(13)中:θ1,θ2為權(quán)重系數(shù),滿足θ1+θ2=1,θ1≥0,θ2≥0,對于不同用戶權(quán)重系數(shù)會有所不同,對于電費支出較高的用戶(如工業(yè)、居民用戶),可以設(shè)置較大的θ1;而對于對用電時間要求較嚴格的用戶(如商業(yè)用戶),可以設(shè)置較大的θ2.具體取值參照表1.本文算例選用工業(yè)用戶,應(yīng)該滿足θ1≥θ2.

    由于削峰填谷和提高用戶滿意度的函數(shù)求解方向相互沖突,無法找到同時滿足它們的最優(yōu)解,故采用設(shè)定權(quán)重法構(gòu)建兼顧兩者,因此,分時電價優(yōu)化模型的目標(biāo)函數(shù)為:

    minFpo=μ1(f1+f2)-μ2(θ1f3+θ2f4)

    (14)

    式(14)中:μ1,μ2為權(quán)重系數(shù),要求與θ1,θ2相同,取值依賴于電網(wǎng)公司對于兩個方面的重視程度.

    (2)約束條件

    不同時期電價約束:實施分時電價后,峰期電價大于平期電價,平期電價大于谷期電價.

    ρp-ρfgt;0

    ρf-ρvgt;0

    (15)

    防止峰谷負荷反轉(zhuǎn)約束:實施分時電價后,避免出現(xiàn)負荷峰谷反轉(zhuǎn)的情況.

    Lp,min-Lv,maxgt;0

    (16)

    式(16)中:Lp,min,Lv,max分別為實施分時電價后峰時負荷最小值和谷時負荷最大值.

    用戶獲益約束:實施峰谷電價前后,保證用戶的利益不受損失,即實施分時電價后用戶單位用電成本不高于分時電價前的電價.

    EL≤E0L0,L=∑24i=1Li

    (17)

    式(17)中:L為實施分時電價后的負荷.

    2 用戶側(cè)儲能配置運行協(xié)同優(yōu)化

    2.1 用戶側(cè)儲能經(jīng)濟模型及運行約束

    2.1.1 成本模型

    用戶側(cè)儲能系統(tǒng)的全壽命周期投資成本主要由固定投資成本和運行維護成本構(gòu)成.

    固定投資成本Csys包括容量成本和功率成本兩部分,分別由額定容量和額定充/放電功率決定,可表示為:

    Csys=ceEmax+cpPmax

    (18)

    式(18)中:Emax為儲能系統(tǒng)的額定容量,Pmax為儲能系統(tǒng)的額定功率,ce為單位容量成本,cp為單位充放電功率成本.

    全壽命周期的運行維護成本Cope與額定功率有關(guān),可表示為:

    Cope=∑Tt=1(copPmax)1+ir1+drt

    (19)

    式(19)中:T為電池壽命,ir為通貨膨脹率,dr為貼現(xiàn)率,cop為單位年運行維護成本.

    2.1.2 收益模型

    用戶側(cè)儲能系統(tǒng)的全壽命周期收益包括削峰填谷收益、需量削減收益、降低用戶側(cè)變壓器成本和儲能電池的回收收益.用電成本包括需量電費和電量電費,需量電費由每月最大負荷值決定,電量電費由用電量決定,用戶側(cè)儲能對兩者都有削減效果,進而產(chǎn)生削峰填谷收益Bele和需量削減收益Bbas.

    一是儲能利用峰谷電價差,在峰時段放電,在谷時段充電,進而獲得削峰填谷收益:

    Bele=∑Tt=1DB1(1+ir1+dr)t

    B1=∑nDi=1ρ(i)[Pdis(i)-Pch(i)]Δt

    (20)

    式(20)中:B1為儲能一天的削峰填谷收益,D為年運行天數(shù),i為時段,Δt為時段的持續(xù)時間,ρ(i)為第i個時段的分時電價,nD為日運行周期內(nèi)總時段數(shù),Pdis(i)為儲能在第i個時段的放電功率,Pch(i)為儲能在第i個時段的充電功率.

    二是儲能通過削減每月最大負荷需量,產(chǎn)生需量削減收益:

    Bbas=∑Tt=1MB21+ir1+drt

    B2=a(Ppeak,m-P′peak,m)

    (21)

    式(21)中:B2為儲能一月的需量削減收益,M為年運行月數(shù),a為基本電價,Ppeak,m為原始月負荷峰值,P′peak,m為加裝儲能后的月負荷峰值.

    用戶加裝儲能后年負荷峰值降低可以減少變壓器成本[24],變壓器的壽命周期一般比儲能系統(tǒng)長,考慮資金的時間性,儲能全壽命周期內(nèi)降低變壓器成本Btran的收益可以表示為:

    Btran=∑Tt=1ftran∑Ttrant=1(1+ir1+dr)t(1+ir1+dr)t

    ftran=(1+k)Cdδy·Ppeak,yωcosφ

    (22)

    式(22)中:ftran為年節(jié)約的變壓器成本,ω為負載率,cosφ為功率因數(shù),Ttran為變壓器的壽命周期,k為變壓器安裝成本占其造價的比值,Cd變壓器的單位成本,δy為年負荷削峰率,可表示為:

    δy=Ppeak,y-P′peak,yPpeak×100%

    (23)

    式(23)中:Ppeak,y為用戶原始年負荷峰值,P′peak,y為加裝儲能后的年負荷峰值.

    儲能電池的回收收益Brec與固定投資成本Cinv有關(guān)[13],可表示為:

    Brec=γCinv

    (24)

    式(24)中:γ為儲能電池的回收系數(shù).

    2.1.3 儲能運行約束

    (1)功率平衡約束

    Pgrid=Pload(i)-Pdis(i)+Pch(i)

    (25)

    式(25)中:Pgrid為第i個時段與電網(wǎng)交換功率,Pload為第i個時段用戶負荷功率.

    (2)充放電功率約束

    0≤Pdis(i)≤PmaxBdis(i)

    0≤Pch(i)≤PmaxBch(i)

    (26)

    式(26)中:Bdis(i),Bch(i)為0-1變量表示儲能充放電狀態(tài),且滿足任意時刻Bdis(i)+Bcha(i)≤1.

    (3)SOC及其連續(xù)性約束

    儲能容量應(yīng)保持在一定范圍內(nèi),SOCmax和SOCmin分別是儲能荷電狀態(tài)(State of Charge,SOC)的上下限,分別取0.9和0.1.為了確保儲能系統(tǒng)連續(xù)運行的周期性,每個調(diào)度周期的初始荷電狀態(tài)要與最終時刻相等,初始荷電狀態(tài)SOC(0)為0.5,約束條件如下:

    SOCmin≤SOC(i)≤SOCmax

    SOC(i)=SOC(i-1)+[ηchPch(i)-Pdis(i)/ηdis]ΔtEmax

    SOC(0)=SOC(nD)

    (27)

    式(27)中:SOC(i)為儲能系統(tǒng)在第i個時段的荷電狀態(tài),ηch,ηdis分別為儲能系統(tǒng)的充、放電效率.

    (4)日吞吐量約束

    儲能電池的充放電狀態(tài)的轉(zhuǎn)換次數(shù)會影響電池壽命和運行成本,而約束儲能日吞吐量可以改善電池狀態(tài)的轉(zhuǎn)換情況.日吞吐量是指日運行周期內(nèi)電池充電量和放電量之和,約束條件如下:

    Qd≤2NEmax

    (28)

    式(28)中:Qd為第d天的日吞吐量,N為最大日循環(huán)次數(shù),本文取N=2,即每日儲能可以充放電兩次.

    (5)削峰填谷約束

    為了需量削減收益有效,用戶母線側(cè)合成負荷不得超過月負荷峰值,同時,月負荷峰值不得超過年負荷峰值,約束條件如下:

    Pload(i)-Pdis(i)+Pch(i)≤P′peak,m

    P′peak,m≤P′peak,y

    (29)

    2.2 用戶側(cè)儲能配置運行協(xié)同優(yōu)化雙層模型

    儲能的規(guī)劃和運行問題是相互聯(lián)系的,本文考慮多時間尺度下的用戶側(cè)儲能經(jīng)濟性,建立了儲能配置運行協(xié)同優(yōu)化雙層模型,框架如圖2所示.上層模型考慮長時間尺度的全壽命周期凈收益對儲能配置問題進行決策,決策變量包括額定功率、額定容量和年削峰率;下層模型考慮短時間尺度的日運行收益對儲能運行策略進行優(yōu)化,決策變量為日運行周期內(nèi)的充放電功率.上層模型將決策變量傳遞給下層模型問題,下層模型在基礎(chǔ)上優(yōu)化找到最優(yōu)解,將結(jié)果返回上層,不斷迭代找到最優(yōu)值[25].

    2.2.1 上層模型

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

    上層模型的優(yōu)化目標(biāo)是儲能系統(tǒng)全壽命周期內(nèi)的凈收益最大化:

    maxF1=Bele+Bbas+Btran+Brec-Cinv-Cope

    (30)

    式(30)中:F1為儲能全壽命周期凈收益.

    (2)約束條件

    上層模型的約束條件包括兩方面,分別是功率平衡約束和投資成本約束,功率平衡約束如式(25)所示,投資成本約束表示為:

    Csys≤Cmaxsys

    (31)

    式(31)中:Cmaxsys為投資成本最大限值.

    2.2.2 下層模型

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

    下層模型的優(yōu)化目標(biāo)是儲能系統(tǒng)日運行周期內(nèi)收益最大化,須將月需量電費平均到每日,最大化儲能系統(tǒng)日運行周期內(nèi)收益可表示為:

    maxF2=B1+B230

    (32)

    式(32)中:F2為儲能系統(tǒng)日運行收益.

    (2)約束條件

    下層模型的約束考慮儲能運行約束,主要包括充放電功率約束、荷電狀態(tài)(SOC)及其連續(xù)性約束、日吞吐量約束和削峰填谷約束,如式(26)~(29)所示.

    2.3 評估指標(biāo)

    本文選用四個經(jīng)濟指標(biāo)對儲能配置優(yōu)化結(jié)果進行評估,分別為固定投資成本Csys、全壽命周期凈收益F1、投資回報率R和投資回收期Ny.Csys和F1上文已定義,R和Ny的表達式如式(33)、(34)所示.

    R=F1(1+ir1+dr)tCsys+Cope×100%

    (33)

    Ny=Csys+CopeCinc+Ctran+Crec×T

    (34)

    3 模型求解

    兩階段聯(lián)合優(yōu)化模型的求解需要通過分階段進行決策[26].在第一階段,采用PSO-SQP混合算法求解分時電價優(yōu)化模型,生成最佳分時電價制定策略,將參數(shù)和優(yōu)化后的負荷曲線傳遞給第二階段,然后基于KKT條件將雙層模型轉(zhuǎn)換為單層模型求解,獲得最優(yōu)儲能配置方案和運行策略,求解過程如圖3所示.

    在第一階段,將PSO算法和SQP算法結(jié)合,PSO算法具有結(jié)構(gòu)簡單,計算效率高和全局搜索能力強等特點,但其局部搜索能力弱,容易為陷入局部最優(yōu);而SQP算法具有很強的局部搜索能力,計算精度高及穩(wěn)定性好,但其性能優(yōu)劣取決于初始值的選擇[27].因此本文結(jié)合兩種算法的優(yōu)點,采用PSO-SQP混合算法求解分時電價優(yōu)化問題,計算流程如圖4所示,主要步驟包括:

    步驟1 輸入負荷數(shù)據(jù),初始化系統(tǒng)參數(shù);

    步驟2 利用PSO算法對目標(biāo)函數(shù)進行優(yōu)化,若兩次迭代最優(yōu)解之差小于某個特定值ε,則認為達到局部最優(yōu),根據(jù)式(14)計算得到最優(yōu)目標(biāo)函數(shù)F(X1)和最優(yōu)變量X1=δv,δf,δp;

    步驟3 將步驟2得到的最優(yōu)變量X1作為SQP算法的初始值,計算得出目標(biāo)函數(shù)最優(yōu)值f(X2)和優(yōu)化變量X2=δ′v,δ′f,δ′p;

    步驟4 將步驟3得出的優(yōu)化變量X2返回至PSO算法中,計算得出目標(biāo)函數(shù)值F(X2);

    步驟5 比較F(X1)和F(X2)函數(shù)值,選擇其中較小的作為最終最優(yōu)解;

    步驟6 判斷是否滿足迭代停止條件,滿足則輸出最優(yōu)解結(jié)束算法,否則,更新粒子速度和位置,返回步驟2.

    在第二階段,所構(gòu)建的儲能配置運行協(xié)同優(yōu)化雙層模型中,上下層之間存在耦合關(guān)系,模型難以直接求解,為此,采用KKT條件通過構(gòu)造拉格朗日函數(shù),如式(35)所示.把下層模型轉(zhuǎn)換成上層模型的約束,進而將雙層模型簡化為等價的單層模型,然后在Matlab2018a環(huán)境中,利用工具箱YALMIP調(diào)用求解器CPLEX對單層模型進行求解,從而得到最優(yōu)解[28].

    L=-∑nDi=1ρ(i)[Pdis(i)-Pch(i)]Δt-a(Ppeak,m-P′peak,m)/30+

    ∑nDi=1λ1(i)SOC(i)-SOC(i-1)-[ηchPch(i)-Pdis(i)ηdis]Δt/Emax+

    λ2SOC(0)-SOC(nD)+∑nDi=1μ11(i)-Pdis(i)+

    ∑nDi=1μ12(i)Pdis(i)-Bdis(i)×Mdis+∑nDi=1μ13(i)Pdis(i)-Pmax+

    ∑nDi=1μ21(i)-Pch(i)+∑nDi=1μ22(i)Pch(i)-Bch(i)×Mch+

    ∑nDi=1μ23(i)Pch(i)-Bch(i)×Mch+∑nDi=1μ3(i)Bdis(i)+Bch(i)-1+

    ∑nDi=1μ41(i)SOCmin-SOC(i)+∑nDi=1μ42(i)SOC(i)-SOCmax+

    ∑nDi=1μ5(i)Pload(i)-Pdis(i)+Pch(i)-P′peak,m+μ6(P′peak,m-P′peak,y)+

    μ7∑nDi=1Pdis(i)+Pch(i)Δt-2NEmax

    (35)

    式(35)中:λ是等式約束對應(yīng)的拉格朗日乘子,μ是不等式約束對應(yīng)的拉格朗日乘子,下標(biāo)編號無實義,Mdis,Mch是很大的常數(shù).

    4 算例分析

    4.1 算例參數(shù)

    為驗證所提模型和策略的有效性,本文以某工業(yè)用戶典型日負荷曲線為算例,工業(yè)用戶滿意度模型的權(quán)重系數(shù)設(shè)置為θ1=0.9,θ2=0.1,初始電價為0.65元/(kW·h),通貨膨脹率為2%,貼現(xiàn)率為8%,需量電費按用戶側(cè)變壓器容量計費,需量電價為34元/(kW·月),變壓器的功率因數(shù)為0.85,負載率為0.75,壽命為20年.儲能電池選用適用性好的鋰離子電池,其回收價值[17]取回收價8 000元/t,參數(shù)如表2所示.

    4.2 分時電價策略優(yōu)化結(jié)果分析

    4.2.1 峰谷時段劃分結(jié)果對比分析

    電力負荷彈性系數(shù)如表3所示.兩種時段劃分方法的結(jié)果對比如表4和表5所示,把輪廓系數(shù)引入評估不同時段間的類內(nèi)和類間距離,輪廓系數(shù)在[-1,1]之間取值,若其值越靠近1,則聚類效果越好,與僅用均方差評估,劃分精度更高,對獲得合理的時段劃分結(jié)果有重要意義.

    4.2.2 分時電價優(yōu)化結(jié)果對比分析

    考慮到電網(wǎng)公司對分時電價優(yōu)化目標(biāo)有不同的側(cè)重點,本文設(shè)置三組權(quán)重系數(shù)賦值方案,方案一:μ1=0.9,μ2=0.1;方案二:μ1=0.5,μ2=0.5;方案三:μ1=0.8,μ2=0.2.采用PSO-SQP算法對模型進行求解,分時電價優(yōu)化結(jié)果如圖5所示.

    由圖5可知,采用分時電價可以有效減少負荷峰值以及峰谷差,降低負荷波動,實現(xiàn)負荷的削峰填谷,而用戶滿意度在不影響正常生產(chǎn)生活的前提下會稍微降低.因此,目標(biāo)函數(shù)權(quán)重系數(shù)的選擇關(guān)鍵在于電價優(yōu)化的目的側(cè)重于削峰填谷還是提高用戶滿意度,對于該工業(yè)用戶以方案一為例進行后續(xù)分析說明.

    分時電價優(yōu)化結(jié)果如表6所示.結(jié)果顯示,實施分時電價后,峰時段的電價高于初始電價,而谷時段的電價則低于初始電價,在電價彈性系數(shù)矩陣的作用下,峰谷電價差可以實現(xiàn)負荷從峰時段到谷時段的轉(zhuǎn)移,進而達到削峰填谷的目的,同時,平時段的電價略低于初始電價,旨在吸引少量負荷從高峰時段轉(zhuǎn)移到平時段.在兩種分時電價下,用戶分別參與需求響應(yīng)的負荷曲線如圖6所示,實施分時電價可以優(yōu)化用戶負荷曲線,為后續(xù)儲能配置優(yōu)化奠定了基礎(chǔ).

    4.3 儲能優(yōu)化結(jié)果分析

    4.3.1 經(jīng)濟性對比分析

    基于上述雙層優(yōu)化模型和求解算法,設(shè)置以下三種情景進行對比.情景一:不考慮分時電價下需求響應(yīng);情景二:考慮分時電價下需求響應(yīng),采用分時電價一策略;情景三:考慮分時電價下需求響應(yīng),采用分時電價二策略.

    儲能優(yōu)化配置結(jié)果如表7所示,經(jīng)濟評估指標(biāo)結(jié)果如表8所示.對比情景一與情景二、三的結(jié)果可知,情景二較情景一的固定投資成本增加了5.20%,但全生命周期凈收益增加了46.89%,投資回報率和投資回收期分別提升了22.81%和12.79%,情景三較情景一的固定投資成本減少了6.40%,全生命周期凈收益增加了36.19%,投資回報率和投資回收期分別提升了25.51%和14.20%,由此可知,考慮分時電價下的需求響應(yīng)對儲能配置的經(jīng)濟性有顯著的提升作用.

    再對比情景二、三的結(jié)果分析可知,情景三較情景二的固定投資成本減少了11.02%,全生命周期凈收益減少了7.28%,而投資回報率和投資回收期分別提升了2.7%和1.6%,情景三下凈收益減少是由于儲能額定功率和額定容量都低于情景二,但收益的減少量要比投資成本節(jié)約量稍微高一些,并且投資回報率和投資回收期都優(yōu)于情景二.因此,在分時電價二策略下配置儲能較分時電價一有更優(yōu)的經(jīng)濟性,綜上所述,驗證了本文所提兩階段優(yōu)化策略的合理性和有效性.

    4.3.2 配置儲能后負荷曲線對比分析

    在用戶側(cè)配置儲能系統(tǒng)前后日負荷曲線如圖7所示.由圖7可以看出,儲能系統(tǒng)可以有效削減峰時段負荷,減少基本電費,同時充放電策略與負荷曲線和分時電價曲線具有強跟隨性,儲能系統(tǒng)在峰時段放電,在谷時段充電,可以同時實現(xiàn)需量削減和峰谷套利兩種收益方式,為用戶節(jié)約用電成本與能源管理成本.

    對比情景二和情景三合成負荷曲線,分析可知分時電價一下儲能的負荷削峰效果比分時電價二下削峰效果更優(yōu),這是由于情景二配置的儲能系統(tǒng)容量大于情景三,情景二的合成負荷曲線在16:00~17:00時段的負荷值幾乎達到了全天的峰值負荷,這是由于情景二的時段劃分結(jié)果導(dǎo)致的,配置儲能后改變了原負荷曲線的雙峰值特性,并在18:00到峰時段,負荷曲線有大幅度削減,可以看出合成負荷曲線平滑度比情景三差,因此,情景三的模型使得儲能系統(tǒng)運行成本更低,也改善負荷曲線波動程度.

    4.3.3 日吞吐量約束驗證

    儲能運行模型中的約束加入了日吞吐量的約束,為驗證該約束條件的合理性,利用典型日負荷曲線數(shù)據(jù),對有無日吞吐量約束條件下的運行優(yōu)化分別仿真,結(jié)果如圖8所示.由圖8可以看出,加入日吞吐量約束可以減少儲能系統(tǒng)狀態(tài)(充電、放電和待機)轉(zhuǎn)換的次數(shù),能有效降低儲能系統(tǒng)的損耗,延長電池壽命,提升系統(tǒng)經(jīng)濟性.

    5 結(jié)論

    本文提出了一種考慮分時電價下需求響應(yīng)的用戶側(cè)儲能兩階段優(yōu)化策略,結(jié)合算例分析,得出以下結(jié)論:

    (1)分時電價策略引導(dǎo)用電,使用戶在峰時減少用電,在谷時增加用電,有效平抑了負荷曲線的波動性,提高了用戶用電的可靠性和經(jīng)濟性.

    (2)所提出的改進移動邊界方法對電價時段(谷、平、峰)進行劃分,引入輪廓系數(shù)作為目標(biāo)函數(shù)的一個因子,充分考慮時段劃分的類內(nèi)緊密度和類間分離度,提高了時段劃分精度,更準確地反映負荷波動性和用電需求變化,進而更好地指導(dǎo)了儲能配置和運行優(yōu)化.

    (3)所建立的儲能配置運行協(xié)同優(yōu)化雙層模型綜合考慮不同時間尺度下的經(jīng)濟性,最大化短期收益的同時考慮到長期效益和能源利用率,實現(xiàn)儲能配置和運行綜合最優(yōu),同時,儲能運行優(yōu)化約束中引入日吞吐量約束,有利于延長儲能壽命,實現(xiàn)更高的收益和可持續(xù)發(fā)展.

    參考文獻

    [1] 張智剛,康重慶.碳中和目標(biāo)下構(gòu)建新型電力系統(tǒng)的挑戰(zhàn)與展望[J].中國電機工程學(xué)報,2022,42(8):2 806-2 819.

    [2] Nadeem F,Hussain S M S,Tiwarip K,et al.Comparative review of energy storage systems, their roles, and impacts on future power systems[J].IEEE Access,2018,7:4 555-4 585.

    [3] 史林軍,楊 帆,劉 英,等.計及社會發(fā)展的多場景用戶側(cè)儲能容量優(yōu)化配置[J].電力系統(tǒng)保護與控制,2021,49(22):59-66.

    [4] Wang L,Hou C,Ye B,et al.Optimal operation analysis of integrated community energy system considering the uncertainty of demand response[J].IEEE Transactions on Power Systems,2021,36(4):3 681-3 692.

    [5] 范 帥,危怡涵,何光宇,等.面向新型電力系統(tǒng)的需求響應(yīng)機制探討[J].電力系統(tǒng)自動化,2022,46(7):1-12.

    [6] 王守相,張善濤,王 凱,等.計及分時電價下用戶需求響應(yīng)的分布式儲能多目標(biāo)優(yōu)化運行[J].電力自動化設(shè)備,2020,40(1):125-132.

    [7] 劉繼東,韓學(xué)山,韓偉吉,等.分時電價下用戶響應(yīng)行為的模型與算法[J].電網(wǎng)技術(shù),2013,37(10):2 973-2 978.

    [8] 姚建剛,付 強,葉 倫,等.考慮峰谷分時電價影響的變電站容量規(guī)劃方法[J].電力系統(tǒng)自動化,2017,41(13):53-61.

    [9] 胡 鵬,艾 欣,張 朔,等.基于需求響應(yīng)的分時電價主從博弈建模與仿真研究[J].電網(wǎng)技術(shù),2020,44(2):585-592.

    [10] 王守相,王瀚樟,趙倩宇,等.面向配電網(wǎng)光伏接納能力提升的分時電價優(yōu)化方法[J].電力系統(tǒng)自動化,2023,47(10):38-46.

    [11] 黃際元,李欣然,常 敏,等.考慮儲能電池參與一次調(diào)頻技術(shù)經(jīng)濟模型的容量配置方法[J].電工技術(shù)學(xué)報,2017,32(21):112-121.

    [12] 田 蓓,王朝暉,張 爽,等.面向風(fēng)光綜合消納的電力系統(tǒng)廣域儲能容量優(yōu)化配置研究[J].智慧電力,2020,48(6):67-72.

    [13] 薛金花,葉季蕾,陶 瓊,等.采用全壽命周期成本模型的用戶側(cè)電池儲能經(jīng)濟可行性研究[J].電網(wǎng)技術(shù),2016,40(8):2 471-2 476.

    [14] 郭久億,劉 洋,郭焱林,等.不同典型用戶側(cè)儲能配置評估與運行優(yōu)化模型[J].電網(wǎng)技術(shù),2020,44(11):4 245-4 254.

    [15] 趙乙潼,王慧芳,何奔騰,等.面向用戶側(cè)的電池儲能配置與運行優(yōu)化策略[J].電力系統(tǒng)自動化,2020,44(6):121-128.

    [16] 丁逸行,徐青山,呂亞娟,等.考慮需量管理的用戶側(cè)儲能優(yōu)化配置[J].電網(wǎng)技術(shù),2019,43(4):1 179-1 186.

    [17] 周欣怡,齊先軍,吳紅斌,等.提高需求響應(yīng)可靠性的儲能優(yōu)化配置研究[J].電力系統(tǒng)保護與控制,2021,49(2):107-114.

    [18] 丁 曦,姜 威,郭創(chuàng)新,等.考慮需求響應(yīng)的電/熱/氣云儲能優(yōu)化配置策略[J].電力建設(shè),2022,43(3):83-99.

    [19] 孔祥玉,楊 群,穆云飛,等.分時電價環(huán)境下用戶負荷需求響應(yīng)分析方法[J].電力系統(tǒng)及其自動化學(xué)報,2015,27(10):75-80.

    [20] Yang H,Wang L,Ma Y.Optimal time of use electricity pricing model and its application to electrical distribution system[J].IEEE Access,2019,99(3):1 991-2 002.

    [21] 高 原,楊賀鈞,郭凱軍,等.考慮負荷季節(jié)特性的電價型需求響應(yīng)最優(yōu)定價策略[J].電力建設(shè),2023,44(1):55-63.

    [22] 黃劍平,陳皓勇,林鎮(zhèn)佳,等.需求側(cè)響應(yīng)背景下分時電價研究與實踐綜述[J].電力系統(tǒng)保護與控制,2021,49(9):178-187.

    [23] 李國榮.基于用戶需求響應(yīng)的分時電價優(yōu)化研究[D].南京:南京郵電大學(xué),2023.

    [24] 麻秀范,陳 靜,余思雨,等.計及容量市場的用戶側(cè)儲能優(yōu)化配置研究[J].電工技術(shù)學(xué)報,2020,35(19):4 028-4 037.

    [25] 吳盛軍,李 群,劉建坤,等.基于儲能電站服務(wù)的冷熱電多微網(wǎng)系統(tǒng)雙層優(yōu)化配置[J].電網(wǎng)技術(shù),2021,45(10):3 822-3 832.

    [26] 王澤爽,陳嘉俊,朱建全,等.計及循環(huán)壽命的儲能優(yōu)化配置與運營策略[J].電力自動化設(shè)備,2021,41(10):75-81.

    [27] Yang Q,Huang Y,Zhang Q,et al.A bi-level optimization and scheduling strategy for charging stations considering battery degradation[J].Energies,2023,16(13):1-15.

    [28] Ding Y,Xu Q,Huang Y .Optimal sizing of user-side energy storage considering demand management and scheduling cycle[J].Electric Power Systems Research,2020,184:106 284-106 294.

    【責(zé)任編輯:陳 佳】

    午夜精品一区二区三区免费看| 欧美bdsm另类| 日韩 亚洲 欧美在线| 又粗又硬又长又爽又黄的视频| 寂寞人妻少妇视频99o| 99久久精品一区二区三区| 亚洲一区高清亚洲精品| 欧美日韩视频高清一区二区三区二| 亚洲欧美日韩无卡精品| 不卡视频在线观看欧美| 免费人成在线观看视频色| 亚洲一区高清亚洲精品| 亚洲综合精品二区| 成人漫画全彩无遮挡| 亚洲国产欧美在线一区| 日本与韩国留学比较| 三级国产精品欧美在线观看| 天天一区二区日本电影三级| a级一级毛片免费在线观看| 丝瓜视频免费看黄片| 狂野欧美白嫩少妇大欣赏| 丰满乱子伦码专区| 国产爱豆传媒在线观看| 看免费成人av毛片| 午夜福利成人在线免费观看| a级毛色黄片| 国产精品美女特级片免费视频播放器| av国产久精品久网站免费入址| 搡老妇女老女人老熟妇| 3wmmmm亚洲av在线观看| 日产精品乱码卡一卡2卡三| 亚洲经典国产精华液单| 欧美激情国产日韩精品一区| 中文在线观看免费www的网站| 秋霞在线观看毛片| 国产精品久久久久久av不卡| 一级a做视频免费观看| 汤姆久久久久久久影院中文字幕 | 久久久久性生活片| 成人午夜高清在线视频| 日本猛色少妇xxxxx猛交久久| 欧美日韩在线观看h| 国产精品一区二区性色av| 尾随美女入室| 真实男女啪啪啪动态图| 干丝袜人妻中文字幕| 99久久精品国产国产毛片| 亚洲精品久久午夜乱码| 日韩精品有码人妻一区| 午夜福利高清视频| 日韩欧美国产在线观看| 在线免费观看不下载黄p国产| 天天躁日日操中文字幕| 亚洲内射少妇av| 自拍偷自拍亚洲精品老妇| 你懂的网址亚洲精品在线观看| 搞女人的毛片| 麻豆成人av视频| 观看免费一级毛片| 最近2019中文字幕mv第一页| 日日摸夜夜添夜夜爱| 国产视频首页在线观看| 日日撸夜夜添| 日本午夜av视频| 亚洲国产高清在线一区二区三| 日韩国内少妇激情av| 久久精品夜色国产| 欧美日韩视频高清一区二区三区二| 亚洲国产精品专区欧美| 国产成人91sexporn| 亚洲成人中文字幕在线播放| 国产极品天堂在线| 美女脱内裤让男人舔精品视频| 免费黄频网站在线观看国产| 精品亚洲乱码少妇综合久久| 欧美成人精品欧美一级黄| 午夜福利在线观看吧| www.色视频.com| 日韩不卡一区二区三区视频在线| 精品人妻视频免费看| 波野结衣二区三区在线| 国产精品1区2区在线观看.| 免费看不卡的av| 欧美+日韩+精品| 秋霞伦理黄片| 秋霞伦理黄片| 女人久久www免费人成看片| 天堂网av新在线| 韩国高清视频一区二区三区| 成年人午夜在线观看视频 | 我的女老师完整版在线观看| 日韩欧美 国产精品| 网址你懂的国产日韩在线| 久久鲁丝午夜福利片| www.av在线官网国产| 特级一级黄色大片| 日韩欧美国产在线观看| 干丝袜人妻中文字幕| 中文字幕制服av| 在现免费观看毛片| 两个人视频免费观看高清| 日日摸夜夜添夜夜爱| 国产精品三级大全| 99久久中文字幕三级久久日本| 精品国产三级普通话版| 国产高清不卡午夜福利| 欧美丝袜亚洲另类| 久久99热这里只频精品6学生| 亚洲国产欧美在线一区| 亚洲成人精品中文字幕电影| 久久精品人妻少妇| 日韩欧美 国产精品| 久久久欧美国产精品| 欧美变态另类bdsm刘玥| 免费观看精品视频网站| 精品99又大又爽又粗少妇毛片| av国产久精品久网站免费入址| 一个人免费在线观看电影| 日日干狠狠操夜夜爽| 69av精品久久久久久| 欧美xxxx黑人xx丫x性爽| 久久久精品欧美日韩精品| 热99在线观看视频| 日韩精品有码人妻一区| 夜夜看夜夜爽夜夜摸| 国产日韩欧美在线精品| 国产成人a区在线观看| 伊人久久国产一区二区| 亚洲av电影不卡..在线观看| 亚洲va在线va天堂va国产| 亚洲精品国产成人久久av| 欧美性猛交╳xxx乱大交人| 精品一区在线观看国产| 听说在线观看完整版免费高清| 日韩成人av中文字幕在线观看| 精华霜和精华液先用哪个| 免费高清在线观看视频在线观看| 国产精品一区二区性色av| 成人亚洲欧美一区二区av| av在线亚洲专区| 日韩在线高清观看一区二区三区| 欧美潮喷喷水| 秋霞伦理黄片| av专区在线播放| 久久久精品免费免费高清| 国产淫语在线视频| 午夜精品一区二区三区免费看| 成人综合一区亚洲| 日韩视频在线欧美| 精品久久久久久久久久久久久| 亚洲国产精品国产精品| 又爽又黄a免费视频| 最近最新中文字幕大全电影3| 乱码一卡2卡4卡精品| 亚洲不卡免费看| 极品少妇高潮喷水抽搐| 日日啪夜夜爽| 中文在线观看免费www的网站| 青春草视频在线免费观看| 国产日韩欧美在线精品| 亚洲人成网站高清观看| 午夜精品一区二区三区免费看| 国产淫语在线视频| 搞女人的毛片| 欧美激情国产日韩精品一区| 精品酒店卫生间| 肉色欧美久久久久久久蜜桃 | 久久精品国产亚洲av天美| 成年av动漫网址| 丝瓜视频免费看黄片| 久久这里有精品视频免费| 久久亚洲国产成人精品v| 亚洲,欧美,日韩| 搞女人的毛片| 欧美不卡视频在线免费观看| 欧美xxⅹ黑人| 亚洲在线自拍视频| 欧美3d第一页| 成人鲁丝片一二三区免费| 欧美成人午夜免费资源| 欧美高清性xxxxhd video| videos熟女内射| 日韩视频在线欧美| 国产免费又黄又爽又色| 精品久久久久久久人妻蜜臀av| 午夜福利成人在线免费观看| 色视频www国产| 国产久久久一区二区三区| 直男gayav资源| 波多野结衣巨乳人妻| 女人久久www免费人成看片| 天堂中文最新版在线下载 | 精品人妻熟女av久视频| 在线免费观看不下载黄p国产| 九九爱精品视频在线观看| 国产一区有黄有色的免费视频 | 免费无遮挡裸体视频| 日本-黄色视频高清免费观看| 亚洲精品亚洲一区二区| 亚洲欧洲国产日韩| 亚洲欧洲国产日韩| 日本免费在线观看一区| h日本视频在线播放| 国产毛片a区久久久久| 国产精品综合久久久久久久免费| av国产免费在线观看| 亚洲精品色激情综合| av在线播放精品| av卡一久久| 午夜福利视频1000在线观看| 狂野欧美白嫩少妇大欣赏| 国产精品伦人一区二区| 久久精品国产鲁丝片午夜精品| 国产精品麻豆人妻色哟哟久久 | 国产真实伦视频高清在线观看| 亚洲怡红院男人天堂| 亚洲精品色激情综合| 亚洲18禁久久av| 夫妻午夜视频| av在线播放精品| 大片免费播放器 马上看| 男人爽女人下面视频在线观看| 免费在线观看成人毛片| 国产高清三级在线| 国产又色又爽无遮挡免| a级一级毛片免费在线观看| 亚洲人成网站在线观看播放| 可以在线观看毛片的网站| 亚洲人成网站在线观看播放| 久久久久免费精品人妻一区二区| 老司机影院毛片| 中文欧美无线码| 爱豆传媒免费全集在线观看| 亚洲av一区综合| 亚洲人与动物交配视频| 国产成人aa在线观看| 国产69精品久久久久777片| 国产成人a区在线观看| 日日摸夜夜添夜夜添av毛片| 国产精品久久久久久久电影| 精华霜和精华液先用哪个| 午夜福利成人在线免费观看| 如何舔出高潮| 国产成人精品一,二区| 欧美+日韩+精品| 中文天堂在线官网| 久久这里只有精品中国| 欧美一区二区亚洲| 久久久久精品性色| 极品教师在线视频| 美女被艹到高潮喷水动态| 国产男女超爽视频在线观看| 日韩精品青青久久久久久| 久久99热这里只有精品18| 亚洲欧洲日产国产| 午夜福利在线观看吧| 三级国产精品片| 乱码一卡2卡4卡精品| 国产伦一二天堂av在线观看| 插逼视频在线观看| 午夜日本视频在线| 色吧在线观看| 国产高清国产精品国产三级 | 99久久人妻综合| 欧美变态另类bdsm刘玥| 亚洲成色77777| 国产黄频视频在线观看| 亚洲自偷自拍三级| 亚洲欧美日韩卡通动漫| 简卡轻食公司| 精品久久久久久久久久久久久| 国产精品不卡视频一区二区| 在线a可以看的网站| 日产精品乱码卡一卡2卡三| 免费电影在线观看免费观看| 91在线精品国自产拍蜜月| 在线免费十八禁| 高清在线视频一区二区三区| 91精品一卡2卡3卡4卡| 国产午夜精品论理片| 伊人久久精品亚洲午夜| 纵有疾风起免费观看全集完整版 | 国产综合精华液| 99久久精品国产国产毛片| 99久国产av精品| 午夜免费观看性视频| 亚洲第一区二区三区不卡| av国产免费在线观看| 少妇的逼好多水| 在现免费观看毛片| 中文字幕免费在线视频6| 日韩欧美一区视频在线观看 | 网址你懂的国产日韩在线| 中文字幕av成人在线电影| 国产一区二区三区av在线| 97在线视频观看| 一级毛片我不卡| videos熟女内射| 成人亚洲欧美一区二区av| 亚洲精品久久久久久婷婷小说| 我的女老师完整版在线观看| 久久久久网色| 久久久久久久亚洲中文字幕| 人人妻人人澡人人爽人人夜夜 | 亚洲熟妇中文字幕五十中出| 91在线精品国自产拍蜜月| 美女高潮的动态| 80岁老熟妇乱子伦牲交| av在线观看视频网站免费| 毛片一级片免费看久久久久| 女人被狂操c到高潮| 亚洲丝袜综合中文字幕| 边亲边吃奶的免费视频| 777米奇影视久久| 成人av在线播放网站| 欧美一级a爱片免费观看看| 国产精品久久久久久久电影| 精品国产三级普通话版| 最近中文字幕高清免费大全6| 五月天丁香电影| 在线观看人妻少妇| 国产一区二区三区综合在线观看 | av福利片在线观看| 97人妻精品一区二区三区麻豆| 亚洲精品国产成人久久av| 国产视频首页在线观看| 国产亚洲精品久久久com| 国产精品久久久久久久久免| 日韩成人伦理影院| 91午夜精品亚洲一区二区三区| 日本猛色少妇xxxxx猛交久久| 男插女下体视频免费在线播放| 乱码一卡2卡4卡精品| 亚洲成人av在线免费| 亚洲av男天堂| 免费观看无遮挡的男女| 99热全是精品| 中文字幕av在线有码专区| 亚洲av日韩在线播放| 免费不卡的大黄色大毛片视频在线观看 | 一边亲一边摸免费视频| 亚洲欧美清纯卡通| 国产v大片淫在线免费观看| 在线播放无遮挡| av网站免费在线观看视频 | 亚洲av一区综合| 韩国av在线不卡| 久久亚洲国产成人精品v| 两个人的视频大全免费| 亚洲国产精品专区欧美| 国产av国产精品国产| 久久久久久久大尺度免费视频| 亚洲欧美成人精品一区二区| freevideosex欧美| 2018国产大陆天天弄谢| 亚洲欧美一区二区三区国产| 亚洲av成人精品一区久久| 高清视频免费观看一区二区 | 色哟哟·www| 欧美成人一区二区免费高清观看| 男女那种视频在线观看| 在线天堂最新版资源| 国产男人的电影天堂91| 国产精品美女特级片免费视频播放器| 日日摸夜夜添夜夜添av毛片| av网站免费在线观看视频 | 国产精品女同一区二区软件| 国产69精品久久久久777片| 伦精品一区二区三区| 久久国产乱子免费精品| 深爱激情五月婷婷| 精品久久久久久成人av| av福利片在线观看| 99久久中文字幕三级久久日本| 亚洲精品日本国产第一区| 看黄色毛片网站| 日本色播在线视频| 女人十人毛片免费观看3o分钟| 成人鲁丝片一二三区免费| 国产色爽女视频免费观看| 成人亚洲精品一区在线观看 | 精品午夜福利在线看| 国产老妇伦熟女老妇高清| 亚洲人成网站在线观看播放| 五月伊人婷婷丁香| 国产高清有码在线观看视频| 亚洲人与动物交配视频| 亚洲性久久影院| 91av网一区二区| 午夜福利视频精品| 亚洲天堂国产精品一区在线| 激情 狠狠 欧美| 能在线免费看毛片的网站| 午夜免费观看性视频| 亚洲av中文av极速乱| 国产精品女同一区二区软件| 婷婷色av中文字幕| 少妇熟女aⅴ在线视频| 久久这里只有精品中国| 性插视频无遮挡在线免费观看| 中文在线观看免费www的网站| 精品国产一区二区三区久久久樱花 | 日韩视频在线欧美| 亚洲精品aⅴ在线观看| 国产综合精华液| 亚洲美女搞黄在线观看| 国产午夜福利久久久久久| 成人特级av手机在线观看| 如何舔出高潮| 99re6热这里在线精品视频| 日本与韩国留学比较| 国产综合精华液| 亚洲天堂国产精品一区在线| 只有这里有精品99| 久久综合国产亚洲精品| 亚洲av成人精品一区久久| 22中文网久久字幕| 男女国产视频网站| 啦啦啦啦在线视频资源| 久久久久久久久久久免费av| 日韩三级伦理在线观看| 少妇人妻精品综合一区二区| 亚洲在久久综合| 男女边吃奶边做爰视频| 青春草国产在线视频| 黄色一级大片看看| 亚洲精品乱久久久久久| 午夜激情福利司机影院| 成人高潮视频无遮挡免费网站| 国产精品人妻久久久影院| av女优亚洲男人天堂| 麻豆成人av视频| 成人国产麻豆网| 成年人午夜在线观看视频 | 在线a可以看的网站| 少妇的逼好多水| 免费看av在线观看网站| 老司机影院毛片| 久久久久久国产a免费观看| 欧美高清性xxxxhd video| 亚洲第一区二区三区不卡| 久久久精品欧美日韩精品| 99热这里只有精品一区| 日本黄大片高清| 日本欧美国产在线视频| 精品亚洲乱码少妇综合久久| 神马国产精品三级电影在线观看| 69av精品久久久久久| 看黄色毛片网站| 国产亚洲91精品色在线| 国产单亲对白刺激| 成人毛片60女人毛片免费| 视频中文字幕在线观看| 成人综合一区亚洲| 国产不卡一卡二| 狂野欧美白嫩少妇大欣赏| 国产成人精品福利久久| 精品久久久久久久久av| 99热这里只有精品一区| ponron亚洲| 丰满少妇做爰视频| 久久久久久久久中文| 久久99热这里只有精品18| 国产美女午夜福利| 亚洲电影在线观看av| 少妇裸体淫交视频免费看高清| 国产高清三级在线| 97人妻精品一区二区三区麻豆| 亚洲国产色片| 老司机影院毛片| 国产精品女同一区二区软件| 日韩欧美国产在线观看| 亚洲欧美精品专区久久| 国产精品嫩草影院av在线观看| 日本熟妇午夜| 天堂√8在线中文| 久久久久久久国产电影| 久久久亚洲精品成人影院| 国产精品久久久久久av不卡| 人人妻人人看人人澡| 97超视频在线观看视频| 麻豆av噜噜一区二区三区| 天堂中文最新版在线下载 | 中文欧美无线码| 精品一区二区三区视频在线| 国产精品一区二区在线观看99 | 成人二区视频| 亚洲欧美精品自产自拍| 一二三四中文在线观看免费高清| av在线亚洲专区| 国产成人精品婷婷| 中文字幕亚洲精品专区| 国产乱人视频| 国产亚洲av嫩草精品影院| 黄色配什么色好看| 欧美一区二区亚洲| 成人一区二区视频在线观看| 亚洲欧洲国产日韩| 麻豆乱淫一区二区| 永久免费av网站大全| 国产精品久久久久久av不卡| 国产午夜精品论理片| 一边亲一边摸免费视频| 高清在线视频一区二区三区| 热99在线观看视频| 麻豆av噜噜一区二区三区| 一区二区三区高清视频在线| 成人无遮挡网站| 亚洲在久久综合| 中文精品一卡2卡3卡4更新| 国产在线一区二区三区精| 精品久久久久久久久亚洲| 非洲黑人性xxxx精品又粗又长| 亚洲成人av在线免费| 女人十人毛片免费观看3o分钟| 欧美变态另类bdsm刘玥| 亚洲精品色激情综合| 成年av动漫网址| 国产午夜精品久久久久久一区二区三区| 久久精品夜色国产| 欧美精品一区二区大全| 一本一本综合久久| 免费大片黄手机在线观看| 欧美性感艳星| 天堂av国产一区二区熟女人妻| 一级毛片 在线播放| 欧美 日韩 精品 国产| 午夜福利在线在线| 午夜精品在线福利| 老司机影院毛片| 少妇熟女欧美另类| 国产精品一区二区性色av| 亚洲婷婷狠狠爱综合网| 搞女人的毛片| 综合色av麻豆| 啦啦啦韩国在线观看视频| 九九久久精品国产亚洲av麻豆| 男人狂女人下面高潮的视频| 午夜日本视频在线| 国产淫片久久久久久久久| av福利片在线观看| 免费观看无遮挡的男女| 久久久久久久久久成人| 久久久久久久午夜电影| 日韩视频在线欧美| av免费在线看不卡| 亚洲美女搞黄在线观看| 日本色播在线视频| 99热这里只有是精品在线观看| av网站免费在线观看视频 | 亚洲va在线va天堂va国产| 国产老妇女一区| 女人十人毛片免费观看3o分钟| 高清视频免费观看一区二区 | 国产精品人妻久久久久久| 欧美另类一区| a级一级毛片免费在线观看| 一边亲一边摸免费视频| 午夜精品一区二区三区免费看| 一级a做视频免费观看| 久久精品熟女亚洲av麻豆精品 | 毛片一级片免费看久久久久| 国产成人一区二区在线| 中国国产av一级| 亚洲av中文字字幕乱码综合| 午夜激情久久久久久久| 中文在线观看免费www的网站| 精品一区二区免费观看| 国产精品麻豆人妻色哟哟久久 | 最近最新中文字幕大全电影3| 久久这里只有精品中国| 亚洲不卡免费看| 婷婷六月久久综合丁香| 日韩欧美国产在线观看| 欧美bdsm另类| 狠狠精品人妻久久久久久综合| 自拍偷自拍亚洲精品老妇| 国产精品一二三区在线看| av在线观看视频网站免费| 国产中年淑女户外野战色| 国产永久视频网站| 精品人妻熟女av久视频| 两个人的视频大全免费| 久久久久精品性色| 22中文网久久字幕| 99热全是精品| 日本与韩国留学比较| 人妻系列 视频| 国产av在哪里看| 亚洲成人一二三区av| 91av网一区二区| 极品教师在线视频| 国内揄拍国产精品人妻在线| 国产大屁股一区二区在线视频| 日韩制服骚丝袜av| 直男gayav资源| 日韩大片免费观看网站| 国产伦精品一区二区三区四那| 国产国拍精品亚洲av在线观看| 国产成人a区在线观看| 2018国产大陆天天弄谢| 亚洲欧美成人精品一区二区| 激情五月婷婷亚洲| av在线亚洲专区| 嫩草影院新地址| 看十八女毛片水多多多| 丰满少妇做爰视频| 国产av码专区亚洲av| 久久国内精品自在自线图片| 日本免费在线观看一区| 久久久久网色| 午夜日本视频在线| 久久99热6这里只有精品| 日本午夜av视频| 日本三级黄在线观看|