半月刊

ISSN 1000-1026

CN 32-1180/TP

+高级检索 English
网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

基于藤Copula理论的海上风电建模及电力市场运行分析

  • 张帅龙 1,2
  • 郑可迪 1,2
  • 刘学 1,2
  • 程兰芬 3
  • 唐翀 3
  • 陈启鑫 1,2
1. 清华大学电机工程与应用电子技术系,北京市 100084; 2. 新型电力系统运行与控制全国重点实验室(清华大学),北京市100084; 3. 南方电网科学研究院有限责任公司,广东省广州市 510663

最近更新:2024-05-30

DOI:10.7500/AEPS20230804004

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

海上风电因具有更高的发电利用小时数,将成为中国实现“碳中和·碳达峰”的重要技术选项。随着海上风电装机容量的不断提升,其出力的间歇性、波动性将对电力系统的安全稳定和电力市场的经济运行造成显著影响。为刻画不同站点间风速的相关性,基于藤Copula理论对海上风电的出力进行建模,并利用C藤Copula函数构建多风电场出力的联合分布模型。在此基础上,提出基于共识聚类的电力市场典型阻塞场景构建方法,并采用简化的中国广东省电网架构数据开展算例仿真计算,对大规模海上风电并网前后电力系统的发电结构、碳排放量、市场电价等关键指标的变化进行实证分析。

0 引言

随着全球气候变化的加剧,能源电力部门朝着绿色低碳方向转型成为各个国家可持续发展的必经之路。截至2023年第一季度末,中国可再生能源装机容量已经达到了1 258 GW

1。其中,海上风电作为一种重要的可再生能源,大规模的并网和消纳将会加剧电网局部潮流的多样性和复杂性,直接影响电力市场的供需平衡。中国自2017年9月启动了以广东、蒙西、浙江等地区为代表的8个现货市场建设试点以来,已经起到了良好的保供稳价作用,但随着市场建设不断发展,将面临新型参与主体纳入、清洁能源消纳机制设计等问2。因此,有必要分析海上风电接入对电力市场经济性的影响。

目前,已有围绕海上风电对系统规划和电力市场影响的分析研究。文献[

3]建立了风光互补出力的Copula函数模型,并提出风光互补电网的扩展规划模型。文献[4]使用Copula函数和模糊C均值聚类法构建了2个风电场的风电功率概率模型,讨论了风电相关性对源网规划的影响。文献[5]提出将海上风电视为负的负荷,利用Nataf变换确定风电场间相关性,并分析并网后市场电价变化情况。文献[6]利用剩余需求曲线对风电参与爱尔兰电力市场造成的电价波动展开分析;文献[7]利用3种时间序列模型对比分析海上风电与陆上风电对批发电价的影响差别;文献[8]利用泛欧洲电力市场模型分析海上风电对区域电价、碳排放的影响。上述研究仅针对海上风电接入前后电价的变化进行分析,并且使用了较为简单的市场模型。

为了分析海上风电对市场的影响,首先,需要构建准确的海上风电模型。对于接入同一区域的不同风电场,其风速受到风活动的影响会呈现显著的随机相关性,加剧海上风电出力的波动性,直接影响系统的运行状态。目前,对风电场的相关性主要采用矩阵变换法和Copula函数进行研究。文献[

9-10]通过Chloesky分解对多维变量进行线性变换,应用矩阵变换法建立风速相关性模型;文献[11-12]将相关性数据转换为独立的正态分布变量,再通过逆Nataf变换建立风速相关性模型,解决了矩阵变换法仅适用于正态分布的问题,但仅计及变量间的线性相关关系。

为描述变量间的非线性相关性,文献[

13]使用Copula函数对4个风电场的风速数据进行联合建模,针对具有非对称上厚尾特性的2个风电场,通过混合Copula函数构建了其出力的相依结构。Copula函数可以很好地描述二维变量间的相关性,但随着变量维数增加,会面临维数灾问题。为克服这一问题,文献[14]基于Pair-Copula函数对多元联合概率密度函数进行分解,建立了多元风速相依模型;文献[15]利用C藤(canonical vine)结构Copula函数对多风电场功率进行建模;文献[16]利用D藤(drawable vine)结构和混合Copula函数对高维数据进行相关性建模,可以更好地拟合出数据的隐藏依赖结构。然而,上述研究只采用单一藤结构将高维变量转化为多个二维变量的Copula函数进行拟合,在实际应用中,因变量间的相关性复杂多变,准确度不高。

此外,由于节点负荷需求与海上风电出力波动性相互叠加,系统在不同时段下的阻塞情况不尽相同,这将对电力市场的经济运行产生直接影响。为了刻画市场长期运行中的典型阻塞分区情况,需要考虑线路阻塞状态实现电网分区划分。然而,目前相关研究主要面向电网功率平衡和电压控制进行分区划分。例如,文献[

17]以电气距离为权重进行无功分区的划分;文献[18]定义电耦合强度进行无功分区的划分;文献[19]利用5种电气距离指标进行分区划分。

针对上述研究现状,本文提出了一种基于藤结构Copula函数的海上风电建模及电力市场经济性分析方法。首先,针对多个海上风电场风速相关性的问题,通过藤结构Copula函数构建海上风电联合分布模型,利用拟合优度指标确定最佳分布模型,获取考虑风速相关性的联合出力数据;然后,提出了基于共识聚类的市场典型阻塞场景构建方法;最后,构建中国广东省电网拓扑,利用真实测算数据对海上风电并网消纳问题进行模拟仿真,对比分析海上风电接入前后系统发电结构、阻塞状态及市场经济性的变化。

1 基于藤结构Copula函数的联合分布建模

Copula函数可用于刻画变量间的非线性关系,本文借助藤结构Copula函数将高维变量转化为二维变量的相关性问题进行求解,并构建多个海上风电场风速的联合分布模型。

1.1 Copula函数

1959年,Sklar提出可以用一个Copula函数将n个变量各自的边缘分布联结起来得到一个n维的联合概率分布函

20,如式(1)所示。

F(x1,x2,,xn)=C(F1(x1),F2(x2),,Fn(xn)) (1)

式中:F(x1,x2,,xn)为各变量的联合分布函数;C(F1(x1),F2(x2),,Fn(xn))为Copula函数;F1(x1),F2(x2),,Fn(xn)分别为变量x1,x2,,xn对应的边缘分布函数。

式(1)等号两侧求导得到联合概率密度函数:

f(x1,x2,,xn)=c(F1(x1),F2(x2),,Fn(xn))o=1nfo(xo) (2)

式中:fx1,x2,,xn为联合概率密度函数;c(F1(x1),F2(x2),,Fn(xn))为Copula函数的概率密度函数;fo(xo)为各变量边缘分布的概率密度函数。

若各个随机变量的边缘分布均为连续函数,则Copula函数具有唯一性。因此,求解多维随机变量的联合概率分布函数只需求取各变量的边缘分布函数和Copula函数。

1.2 藤结构Copula函数

藤结构Copula函数是一种基于分位数的高维概率分布建模方

21,按照一定的树结构将高维随机变量之间的相关性转化为两两变量之间的相关性。以6维随机变量为例,不同藤结构下的连接示意图如附录A图A1所示。C藤结构适用于多个变量对某一随机变量的相关性均显著强于其余变量的情况;而各个随机变量之间的相关性较为平均时,D藤结构呈现出更好的拟合准确度。但对于高维变量而言,其变量间的相关性更为复杂,R藤结构往往具有更好的普适性。

对于d维随机变量x1,x2,,xd,具有F1(x1),F2(x2),,Fn(xn)的边缘分布函数,其R藤结构包括d-1层树,记作To,o=1,2,,d-1,每层树为一个包含No个节点、Eo条边连接的无向图,To=No,Eo满足:

No=d-o+1Eo=d-o (3)

式中:o=1,2,,d-1

此外,当树To+1中2个节点由同一边相连时,须保证这2个顶点在树To中对应2条边之间有一个公共节点。确定R藤的树结构后,其联合概率密度函数为:

fR(x1,x2,,xd)=q=1dfq(xq)r=1d-1eEocϕ(e),φ(e)ψe(Fϕ(e)ψe(xϕ(e)|xψe),Fφ(e)ψe(xφ(e)|xψe)) (4)

式中:r为树Tr的层数序号;ψe为边e条件变量编号的集合;ϕ(e)φ(e)为边e关系变量编号的集合;Fϕ(e)ψe(xϕ(e)|xψe)Fφ(e)ψe(xφ(e)|xψe)分别为变量xϕ(e)xφ(e)的条件分布函数;cϕ(e),φ(e)ψe(Fϕ(e)ψe(xϕ(e)|xψe), Fφ(e)ψe(xφ(e)|xψe))为对应变量条件分布的Copula函数的概率密度函

22

对于C藤结构下的每一层树To,均包含一个根节点Noroot,对第1层树的第1个根节点使用二元Copula函数对每对节点进行建模。对于树To,o2,以No-1root变量为条件变量,对其余变量进行相关性建模,C藤结构具有星形结

23-24。确定C藤的树结构后,其联合概率密度函数为:

fC(x1,x2,,xd)=q=1dfq(xq)r=1d-1h=1d-rcr,r+h1,2,,(r-1)(F(xh|x1,x2,,xh-1),F(xh+r|x1,x2,,xh-1)) (5)

对于D藤结构下的每一层树To,通过变量的固定顺序进行构建,若变量顺序为1,2,,d,则对第1层树中节点对(1,2),(2,3),,(d-1,d)进行二元Copula函数建模。对于树To,o2,采用同样顺序进行相关性建模,D藤结构具有路径结

24。确定D藤的树结构后,其联合概率密度函数为:

fD(x1,x2,,xd)=q=1dfq(xq)r=1d-1h=1d-rch,h+r(h+1),(h+2),,(h+r-1)(F(xh|xh+1,xh+2,,xh+r-1),F(xh+r|xh+1,xh+2,,xh+r-1)) (6)

1.3 模型构建与数据采样

本文根据序贯估计

22,通过Kendall秩相关系数γ确定各藤结构下各级树结构,计算各边Copula函数,最终确定最优的藤Copula函数模型,采样得到多维变量考虑相关性的模拟数据。具体求解过程如下:

步骤1:通过式(7)利用威布尔分布拟合各个风电场风速vw数据的概率密度函数fW,w(vw),计算得到其边缘分布函数FW,w(vw)

fW,w(vw)=abvwba-1e-vwba (7)

式中:a为形状参数;b为尺度参数。

步骤2:对第1层树所有节点对通过式(8)进行Kendall秩相关系数γ的计算,结合藤结构类型确定根节点及树结构。

γ=Nc-Nd12n(n-1) (8)

式中:Nc为两随机变量间的同序对数;Nd为两随机变量间的异序对数。|γ|的值越接近于1,表示随机变量之间相关性越强;|γ|的值越接近于0,表示随机变量之间相关性越弱。

步骤3:对第1层树的每条边估算二元Copula函数的参数值,利用赤池信息准则(Akaike information criterion,AIC)辅助确定最优二元Copula函数。

步骤4:利用步骤3所得Copula函数采样得到第2层树对应节点值,重复步骤2、3确定第2层树结构及Copula函数。

步骤5:重复上述步骤直至求得藤结构下所有树结构及Copula函数,得到藤Copula函数模型。

步骤6:基于蒙特卡洛法生成取值范围在[0,1]间独立均匀分布的随机数u1,u2,,ud,由随机数组成的向量u=[u1,u2,,ud]

步骤7:根据式(9)式(10)迭代求解得到采样数

25F1(v1),F2(v2),,Fd(vd)

F(vw|v)=Cvw,vkv-k(F(vw|v-k),F(vk|v-k))F(vk|v-k) (9)
u1=v1u2=F2|1(v2|v1)u3=F3|1,2(v3|v1,v2)ud=Fd|1,2,,d-1(vd|v1,v2,,vd-1) (10)

式中:v为包含vk的模拟点向量;v-k为向量v中去除vk的剩余向量;F(vw|v)F(vw|v-k)F(vk|v-k)为相应变量条件分布的联合分布函数;Cvw,vkv-k(F(vwv-k),F(vkv-k))为相应变量条件分布的Copula函数。

步骤8:基于vk=Fk-1(Fk(vk)),可得考虑相关性的风速数据,重复K次步骤6、7,可得K组风速数据v1,y,v2,y,,vd,y,y=1,2,,K

步骤9:确定最优藤结构Copula模型。根据式(11)计算各模型的Cramér-Von Mises(CvM)距

26以评估各模型拟合效果优劣。

DnCvM=n[0,1]d(Cn(u)-Ce(u))2dCe(u) (11)

式中:Cn(u)为构建的Copula函数;Ce(u)为经验Copula函数。CvM距离为检验的目标Copula函数与经验Copula函数之间的距离,该距离越小,模型的拟合准确度越高。

2 基于共识聚类的典型阻塞场景构建

考虑到海上风电的间歇性与波动性,在并网消纳的同时会对系统的潮流分布产生影响,造成不同时段输电阻塞线路的变化,进而影响系统各个节点的电价分布。不同时段的阻塞情况本质为节点的划分结果,共识聚类则可以将多个聚类划分结果进行再聚

27,以得到系统典型阻塞状态。

2.1 电网阻塞场景构建

电网由节点与线路构成,电网的分区划分本质为各节点的分区划分,即一个包含所有节点的集合U划分为I个包含各个节点的子集U1,U2,,UI,满足:

U1U2UI=UUi1Ui2=        i1=1,2,,I;i2=1,2,,I;                                   i1i2 (12)

本文利用节点边际阻塞成本(marginal congestion cost,MCC)作为输入特征,构建M个阻塞场景集合V1,V2,,VM。MCC为线路出现阻塞时,由线路的阻塞价格按照功率转移分布因子的比例所形成的价格信号:

λj,tMCC=σlGl,j (13)

式中:λj,tMCC为线路l阻塞引起节点j的MCC,λj,tMCC>0时,线路l的阻塞使得节点j电价升高,电能“稀缺”,λj,tMCC<0时,线路l的阻塞使得节点j电价降低,电能“充裕”;σl为线路l的阻塞价格;Gl,j为节点j对线路l的输出功率转移分布因子,是功率转移分布因子矩阵的第l行、第j列的元素。

通过阻塞分量划分得到系统各节点在各时段下的阻塞场景集合V1,V2,,VM,基于共识聚类找到最优的子集Vopt,使得其与其他划分集合间综合相似度最大,即:

Vopt=argmax s(Vg,Vm) (14)

式中:s(Vg,Vν)为划分子集间的相似度;VgVm分别为第gm个划分子集。

2.2 共识聚类求解方法

基于共识聚类对式(14)的划分求解包含以下步骤。

步骤1:获取节点相似度矩

28。对于目标区域的M个划分子集V1,V2,,VM,定义对应的节点隶属矩阵分别为H1,H2,,HM,第m个矩阵的元素Hμ,ν,m用以描述每种划分子集中节点μ是否隶属于分区Uνν=1,2,,I):

Hμ,ν,m=1μUν0μUν (15)

为获得最大的子集划分相似度,须选择2个节点划分至同一子集的方法。定义在所有的子集划分中,2个节点被划分到同一子集的比例为该节点对的节点相似度矩阵Snode,可通过式(16)求解。

Snode=1MHHT (16)

式中:H=[H1,H2,,HM],为目标区域所有节点隶属矩阵的级联矩阵。

步骤2:构建拉普拉斯矩阵。定义度矩阵SD为节点相似度矩阵的对角矩阵,则拉普拉斯矩阵可通过式(17)求解。

L=SD-Snode (17)

步骤3:特征值分解。拉普拉斯矩阵L为实对称矩阵,可以进行对角化和特征值分解(谱分解),并基于特征向量构建特征矩阵Λ(谱特征)。

步骤4:聚类。使用k-均值聚类算法将特征矩阵Λ进行聚类,将每个样本划分到其所属的分区中,从而得到最终的分区划分结构。

3 现货市场出清模型

为实现市场运行的准确模拟,本文采用日前安全约束经济调度(security constrained economic dispatch,SCED)模型模拟市场出清过程,并合理设置不同类型主体的报价行为。

3.1 日前SCED模型

SCED模型以最小化发电成本为目标,确定满足负荷需求的机组实时出力计划,其目标函数为:

miniΩGtOi,t(Pi,tG) (18)

式中:ΩG为机组编号集合;Pi,tG为机组i在时段t的有功功率;Oi,t(Pi,tG)为机组i在时段t的运行成本。

约束条件包括机组出力的上下限约束、系统负荷平衡约束及线路潮流约束,如式(19)所示。

Pi,tGminPi,tGPi,tGmaxiΩGPi,tG=PtD-PlmaxiΩGGl,iPi,tG-jGl,jPj,tDPlmax (19)

式中:Pi,tGminPi,tGmax分别为机组i在时段t允许的最小有功功率和最大有功功率;PtD为时段t的系统负荷;Plmax为线路l的潮流传输极限;Gl,i为机组i所在节点对线路l的输出功率转移分布因子;Pj,tD为节点j在时段t的母线负荷值。

3.2 市场主体报价行为设置

市场模型采用发电侧的单边竞价模型,用户侧只需考虑其用电计划,无须申报电价信息,并采用刚性约束满足所有负荷需求。

发电侧主体包含火电、核电、水电(含抽水蓄能)、高压直流输入功率和海上风电,其中,火电、核电、水电采用块报价的报价策略,机组i的块报价需要提交B组电量-电价对:(Qi,1,ρi,1), (Qi,2,ρi,2),,(Qi,b,ρi,b),,(Qi,B,ρi,B)。其中,1bBQi,b0ρi,b+1ρi,b,即机组ib块电量Qi,b的出清价格不低于ρi,b

高压直流输入功率作为省间电量交易在现货市场中固定为低报价ρHVDC,以满足基础电量需求。

考虑海上风电的发电边际成本几乎为0,在现货市场中,海上风电应尽可能申报低价以避免弃风,获得最大利润。因此,本文将海上风电视为负需求的负荷参与现货市

5

4 算例分析

4.1 算例概述

本文基于真实数据,结合文献[

29]构建了简化的中国广东电网拓扑,如附录A图A2所示。

具体而言,本文所构建的广东电网拓扑共包括113个节点,其中,机组节点60个、负荷节点53个,以及输电线路132条。除12台火电机组、9个直流逆变站、7台海上风电机组、4台核电机组、3台抽水蓄能机组外,还包含25台本地机组,对应编号82,83,…,106的节点,均为220 kV网络下的小型火电机组。各类型机组装机容量见附录A表A1。其中,海上风电根据《广东省海上风电发展规划(2017—2030年)(编修)

30中粤东、粤西、珠三角近海浅水区海上风电场规划,设置湛江市、阳江市、珠海市、惠州市、汕尾市、揭阳市、汕头市7个城市的海上风电接入站点。各个站点的装机容量、地理位置、电力接入站点信息见附录A表A2。为简化起见,将海上风电场总出力等效为多个海上风电机组出力的代数和。各站点风速数据取自哥白尼气候变化服务商(Copernicus Climate Change Service,C3S)提供的ERA5数31,选取2000至2022年共计201 624 h的风速数据;各海上风电场的风速基本特性见附录A表A3。

4.2 多风电场藤Copula模型构建

对7处海上风电场的风速分别构建R藤、C藤、D藤Copula函数联合分布模型。为验证各模型拟合效果优劣,利用式(11)计算各模型的CvM距离,结果如表1所示。

表1  藤Copula函数模型的CvM距离
Table 1  CvM distance of vine copula function models
模型类型CvM距离
R藤 Copula函数 0.094 18
C藤 Copula函数 0.068 01
D藤 Copula函数 0.125 42

表1可见,C藤结构Copula函数模型的CvM距离最小,D藤结构Copula函数模型的CvM距离最大,R藤结构Copula函数模型的CvM距离稍小。因此,C藤结构Copula函数模型的拟合相关性最佳,选用该结构进行海上风电场站点风速的联合建模与采样,所得模型各层树结构及采样数据和原始数据的分布图如附录A图A3所示,每层树内各变量间的二元Copula函数类型及相关性结果见附录A图A4。

4.3 结果分析

为分析海上风电并网消纳对广东电力系统及市场的影响,以明阳智能的MySE11型海上风电机组的参数为各海上风电场的典型机组,其切入、切出风速分别为3、25 m/s,风机轮毂高度为140 m,额定运行功率为11 MW,各站点尾流效应系数取0.9,风电机组可利用率取0.95。

算例模拟广东电网在海上风电场并网前后全年8 760 h的运行情况,市场出清间隔设为1 h,对比分析海上风电场接入前后系统发电结构、阻塞状态及市场经济性的变化。

4.3.1 系统发电结构分析

风电场实际发电利用小时数是指场内风电机组在额定功率下的运行小时数,用于表征风电场内发电设备利用率,可通过式(20)计算。

τW=Qτ,WPNW (20)

式中:τW为风电场的发电利用小时数;Qτ,W为风电场在时间τ内的实际发电量;PNW为风电场额定输出功率。

中国海上风电年平均发电利用小时数约2 500 h,陆上风电年平均发电利用小时数约2 000 h

32。计算广东省海上风电场模拟出力,结果如表2所示。各站点的全年实际发电利用小时数在2 130~3 736 h之间。

表2  广东省海上风电场模拟出力
Table 2  Simulated power output of offshore wind farms in Guangdong province, China
接入站点装机容量/MW有效出力/MW发电利用小时数/h
湛江市 1 200 340.6 2 486
阳江市 3 000 1 037.7 3 030
珠海市 500 121.6 2 130
惠州市 1 000 329.8 2 889
汕尾市 1 400 597.1 3 736
揭阳市 900 358.6 3 490
汕头市 1 850 773.9 3 664
合计 9 850 3 559.2 3 165

全省海上风电全年平均有效出力为3 559.2 MW,实际发电利用小时数为3 165 h,超出中国海上风电年平均发电利用小时数665 h,具有较高的发电设备利用率。图1为海上风电并网消纳前后各类型机组的实际出力变化量,其中,海上风电的并网消纳降低了接入站点当地火电机组的出力水平,如揭阳市、阳江市。得益于风能较低的边际成本,接入站点周边区域,例如,深圳市、潮州市,其火力机组出力水平也出现明显降低。部分接入站点城市火电机组出力水平下降幅度小于非接入站点,例如,湛江市、惠州市,这是由于城市原有的火电装机容量基数较小。

图1  海上风电场并网前后机组出力变化

Fig.1  Power output changes of generation units before and after grid-integration of offshore wind farms

海上风电场并网前后,系统火电的有效出力占比降低约3.58%,共计约28.9 TW∙h。设定水电、核电、风电机组的碳排放系数为0 g/(kW∙h),火电机组的平均碳排放系数为873 g/(kW∙h),系统的总碳排放量降低约2.5×107 t。根据《广东省海上风电发展规划(2017—2030年)(编修)

30,广东省到2030年底预计建成投产海上风电装机容量约30.00 GW,在现有算例基础上提高海上风电装机容量,计算各装机容量下系统发电结构的变化,所得结果如表3所示。

表3  海上风电装机容量对系统发电结构的影响
Table 3  Impact of offshore wind power installed capacity on power system generation structure
装机容量/GW火电发电量/(TW∙h)核电发电量/(TW∙h)水电发电量/(TW∙h)海上风电发电量/(TW∙h)
0 4 33.49 1 68.54 66.576 0
9.85 4 04.57 1 66.39 66.472 31.179
15.00 3 89.76 1 65.14 66.232 47.479
20.00 3 75.56 1 64.00 65.836 63.208
25.00 3 61.52 1 62.96 65.361 78.765
30.00 3 47.95 1 61.85 64.990 93.812

海上风电场的并网有效降低了火电的发电结构占比,有助于促进电力产业的绿色转型。但是,由于海上风电出力的随机性,提高其装机容量,将在负荷低谷时直接影响其他电源的出力,增加系统对灵活性资源的需求。

4.3.2 系统阻塞状态分析

计算海上风电场并网前后系统全年发生阻塞的总时长,并网前系统共发生阻塞5 397 h,并网后发生阻塞5 362 h,这是由于海上风电的并网消纳缓解了部分单线路阻塞情况,减少了系统阻塞时长。然而,在部分多线路阻塞情况下,额外的线路阻塞尽管不增加系统阻塞时长,但会加剧系统的阻塞程度。计算具体阻塞线路在海上风电场并网前后的阻塞情况,所得结果如表4所示。

表4  海上风电场并网前后系统阻塞情况对比
Table 4  Comparison of power system congestion before and after grid-integration of offshore wind farms
海上风电场并网前海上风电场并网后
输电线路注入节点流出节点阻塞时长/h输电线路注入节点流出节点阻塞时长/h
7 4 7 70 7 4 7 12
15 10 64 4 309 15 10 64 4 325
29 20 51 33 29 20 51 90
33 25 26 696 33 25 26 1 763
73 54 70 4 312 73 54 70 3 924
79 58 60 3 512 79 58 60 3 651
86 65 66 161 86 65 66 804
93 70 72 33 4 4 5 32

考虑海上风电的并网消纳对附录A图A2中广东省西翼、粤西南送出、东西电力交换、东翼南部沿海送出4个关键断

29造成的影响,其中,西翼送出断面中,线路7阻塞时长减少36 h;粤西南送出断面中,线路15、29、33阻塞时长共增加1 140 h;东西送出断面中,线路29、79阻塞时长共增加196 h;东翼南部沿海送出断面中,线路73阻塞时长减少388 h,线路86阻塞时长增加643 h。

利用基于共识聚类构建的广东电网典型阻塞分区,计算不同分区的节点边际电价,由高至低依次记为分区1至分区6,不同的阻塞分区反映了网络出现阻塞时各节点的电能稀缺度,所得结果见附录A图A5。

结合表4数据,在海上风电场并网后,虽然西翼送出关键断面的用电传输压力得以缓解,但整体上海上风电场并网会进一步增加电网输电断面的传输压力。面对不断创高的电力需求,需进一步完善网络拓扑以服务大规模清洁能源的发展与省外电力的输入。

4.3.3 系统经济性分析

设置发电侧主体的报价如附录A表A4所示,模拟测算海上风电场并网后的节点边际电价,所得广东省各市全年8 760 h的节点边际电价分布如图2所示。图中:a、b、c、d、e、f、g、h、i、j、k、l、m、n、o、p、q、r、s、t、u分别代表湛江市、茂名市、阳江市、云浮市、江门市、肇庆市、佛山市、中山市、珠海市、清远市、广州市、东莞市、深圳市、韶关市、惠州市、河源市、汕尾市、梅州市、揭阳市、汕头市、潮州市。

图2  海上风电场并网后的节点边际电价分布

Fig.2  Distribution of locational marginal prices after grid-integration of offshore wind farms

广东省各市的节点边际电价介于286.00~624.39元/(MW∙h),其中,最低电价与最高电价分别对应第773、4 983 h,系统平均节点边际电价为510.16元/(MW∙h),对比海上风电场未并网时515.04元/(MW∙h)下降了4.88元/(MW∙h)。

考虑到目前广东省海上风电装机容量占比较小,其全年平均有效出力占比约为3.8%,对市场电价影响较小,提高海上风电装机容量至30.00 GW,系统节点边际电价变化如表5所示。

表5  海上风电装机容量对节点边际电价的影响
Table 5  Impact of offshore wind power installed capacity on locational marginal prices
装机容量/GW装机容量占比/%平均节点边际电价/(元∙(MW∙h-1)峰值节点边际电价/(元∙(MW∙h-1)
0 0 515.04 629.79
9.85 4.90 510.10 653.00
15.00 7.47 507.47 657.39
20.00 9.96 503.41 658.70
25.00 12.45 498.38 665.55
30.00 14.94 491.26 661.97

由于现货市场的边际出清原则,海上风电的低边际成本使得火电等高成本机组出清受限,系统平均节点边际电价随装机容量的增加而降低。节点边际电价由电能量价格、阻塞价格、网损价格组成,受到海上风电接入站点均沿海分布的影响,大规模的风电并网消纳会加剧系统供需在空间分布的不平衡,导致部分节点阻塞价格涨幅高于电能量价格的下降幅度,表现出峰值电价的升高结果。

5 结语

在广东省面临大规模海上风电并网和消纳的工程背景下,本文提出了一种仿真分析方法,针对多个海上风电场彼此间存在出力相关性的问题,基于广东省海域风力资源的实际数据,利用藤结构Copula函数构建联合分布模型,其中,C藤结构对于广东省海上风电场站点的风速拟合具有更好的拟合效果。本文构建了220、500 kV电压等级下的广东电网拓扑结构,模拟仿真海上风电场并网前后全年8 760 h出清结果,利用共识聚类构建系统典型阻塞分区。算例结果表明,海上风电的并网消纳有利于降低广东省各地的节点边际电价,但机组出力的间歇性与波动性也将进一步造成电价的波动。此外,海上风电并网会缓解西翼送出断面的输电压力,但粤西南、东翼南部沿海送出断面的输电压力会显著增加,需进一步改进、完善网络拓扑以适应大规模清洁能源的发展。

由于目前海上风电的装机规模较小,本文设置其参与现货市场方式为完全消纳,在未来市场占比提高后可逐步放开报价,直至完全参与市场竞争,研究其在不同报价策略下对现货市场运行的影响。

附录

附录A

图A1  随机变量藤结构示意图

Fig.A1  Schematic diagrams of vine trees of random variables

图A2  简化的广东电网拓扑结构

Fig.A2  Simplified topology of Guangdong power grid

表A1  广东电网各类型机组装机容量
Table A1  Installed capacity of different generation units in Guangdong power grid
机组类型装机容量/MW
火电 1 234 25
核电 1 977 2
水电/抽水蓄能 7 600
特高压、超高压直流接入 4 020 0
海上风电 9 850
总容量 2 008 47
表A2  广东省海上风电场信息
Table A2  Informations of offshore wind farms in Gunagdong province
接入站点装机容量/MW北纬位置东经位置
湛江市 1 200 20°40′ 110°40′
阳江市 3 000 21°25′ 112°25′
珠海市 500 22°00′ 113°30′
惠州市 1 000 22°25′ 115°00′
汕尾市 1 400 22°39′ 116°17′
揭阳市 900 22°44′ 116°29′
汕头市 1 850 23°00′ 117°00′
表A3  海上风电场风速特性
Table A3  Wind speed characteristics of offshore wind farms
接入站点最大风速(m/s)平均风速(m/s)方差(m/s)偏度(m/s)峰度(m/s)
湛江市 33.73 6.75 8.00 0.56 4.98
阳江市 36.44 7.45 10.09 0.53 4.59
珠海市 33.94 6.27 7.68 0.73 5.01
惠州市 40.08 7.25 9.97 0.35 3.88
汕尾市 39.98 8.53 15.27 0.22 2.70
揭阳市 36.71 8.11 13.26 0.22 2.88
汕头市 35.62 8.42 15.11 0.20 2.62

图A3  C藤Copula模型拟合结果

Fig.A3  C-vine copula model fitting results

图A3中,上三角部分为采样数据(灰)与原始数据(黑)的散点分布图,采样数据较好的保有原始数据间的相关性。下三角部分为拟合的C藤Copula模型树结构,从下至上分别为第1至6层树结构,其中揭阳市对应站点与剩余六站点的风速相关性较高,作为第1层树的根节点进行二元Copula函数拟合,拟合结果在第2层树中选取珠海市对应变量作为根节点进行Copula拟合,第3至5层依次选取汕尾市、阳江市、惠州市作为根节点进行拟合,第6层树对剩余两变量进行分布拟合。

图A4  C藤Copula各层树结构及相关系数矩阵

Fig.A4  Correlation coefficients matrix and tree structure of C-vine copula model

图A5  海上风电并网后的广东电网典型阻塞分区

Fig.A5  Typical congestion partition of Guangdong power grid after grid-integration of offshore wind farms

表A4  发电侧主体报价信息
Table A4  Bidding information of entities on generating side
机组类型报价方式机组报价范围(元/(MW∙h))
500kV火电机组 块报价 453.75~703.13
220kV火电机组 块报价 507.50~601.25
核电 块报价 305.90~483.26
水电/抽水蓄能 块报价 254.00~317.50
特高压输入 固定报价 0
海上风电 固定报价 0

参 考 文 献

1

张智刚康重庆.碳中和目标下构建新型电力系统的挑战与展望[J].中国电机工程学报2022428):2806-2819. [百度学术] 

ZHANG ZhigangKANG Chongqing. Challenges and prospects for constructing the new-type power system towards a carbon neutrality future[J]. Proceedings of the CSEE2022428): 2806-2819. [百度学术] 

2

陈启鑫房曦晨郭鸿业.电力现货市场建设进展与关键问题[J].电力系统自动化2021456):3-15. [百度学术] 

CHEN QixinFANG XichenGUO Hongyeet al. Progress and key issues for construction of electricity spot market[J]. Automation of Electric Power Systems2021456): 3-15. [百度学术] 

3

李晖高涵宇张艳.考虑相关性的大规模风光互补电网扩展规划[J].电网技术2018427):2120-2127. [百度学术] 

LI HuiGAO HanyuZHANG Yanet al. Expansion planning of large scale hybrid wind-photovoltaic transmission network considering correlation[J]. Power System Technology2018427): 2120-2127. [百度学术] 

4

张玥王秀丽曾平良.基于Copula理论考虑风电相关性的源网协调规划[J].电力系统自动化2017419):102-108. [百度学术] 

ZHANG YueWANG XiuliZENG Pinglianget al. Copula theory based generation and transmission co-planning considering wind power correlation[J]. Automation of Electric Power Systems2017419): 102-108. [百度学术] 

5

MORALES J MCONEJO A JPÉREZ-RUIZ J. Simulating the impact of wind production on locational marginal prices[J]. IEEE Transactions on Power Systems2010262): 820-828. [百度学术] 

6

DE LA NIETA A A SCONTRERAS JMUÑOZ J Iet al. Modeling the impact of a wind power producer as a price-maker[J]. IEEE Transactions on Power Systems2014296): 2723-2732. [百度学术] 

7

HOSIUS ESEEBAß J VWACKER Bet al. The impact of offshore wind energy on Northern European wholesale electricity prices[J]. Applied Energy2023341120910. [百度学术] 

8

JANSEN MDUFFY CGREEN T Cet al. Island in the sea: the prospects and impacts of an offshore wind power hub in the North Sea[J]. Advances in Applied Energy20226100090. [百度学术] 

9

蔡德福陈金富石东源.风速相关性对配电网运行特性的影响[J].电网技术2013371):150-155. [百度学术] 

CAI DefuCHEN JinfuSHI Dongyuanet al. Impact of wind speed correlation on operation characteristics of distribution network[J]. Power System Technology2013371): 150-155. [百度学术] 

10

邓威李欣然徐振华.考虑风速相关性的概率潮流计算及影响分析[J].电网技术2012364):45-50. [百度学术] 

DENG WeiLI XinranXU Zhenhuaet al. Calculation of probabilistic load flow considering wind speed correlation and analysis on influence of wind speed correlation[J]. Power System Technology2012364): 45-50. [百度学术] 

11

QIN Z LLI W YXIONG X F. Generation system reliability evaluation incorporating correlations of wind speeds with different distributions[J]. IEEE Transactions on Power Systems2013281): 551-558. [百度学术] 

12

潘雄周明孔晓民.风速相关性对最优潮流的影响[J].电力系统自动化2013376):37-41. [百度学术] 

PAN XiongZHOU MingKONG Xiaominet al. Impact of wind speed correlation on optimal power flow[J]. Automation of Electric Power Systems2013376): 37-41. [百度学术] 

13

蔡菲严正赵静波.基于Copula理论的风电场间风速及输出功率相依结构建模[J].电力系统自动化20133717):9-16. [百度学术] 

CAI FeiYAN ZhengZHAO Jingboet al. Dependence structure models for wind speed and wind power among different wind farms based on copula theory[J]. Automation of Electric Power Systems20133717): 9-16. [百度学术] 

14

WANG ZWANG W SLIU Cet al. Probabilistic forecast for multiple wind farms based on regular vine copulas[J]. IEEE Transactions on Power Systems2018331): 578-589. [百度学术] 

15

KIM DKIM J MLIAO S Met al. Mixture of D-vine copulas for modeling dependence[J]. Computational Statistics & Data Analysis2013641-19. [百度学术] 

16

吴巍汪可友李国杰.基于Pair Copula的多维风电功率相关性分析及建模[J].电力系统自动化20153916):37-42. [百度学术] 

WU WeiWANG KeyouLI Guojieet al. Correlation analysis and modeling of multiple wind power based on pair copula[J]. Automation of Electric Power Systems20153916): 37-42. [百度学术] 

17

于琳孙莹徐然.改进粒子群优化算法及其在电网无功分区中的应用[J].电力系统自动化2017413):89-95. [百度学术] 

YU LinSUN YingXU Ranet al. Improved particle swarm optimization algorithm and its application in reactive power partitioning of power grid[J]. Automation of Electric Power Systems2017413): 89-95. [百度学术] 

18

赵晶晶贾然陈凌汉.基于深度学习和改进K-means聚类算法的电网无功电压快速分区研究[J].电力系统保护与控制20214914):89-95. [百度学术] 

ZHAO JingjingJIA RanCHEN Linghanet al. Research on fast partition of reactive power and voltage based on deep learning and an improved K-means clustering algorithm[J]. Power System Protection and Control20214914): 89-95. [百度学术] 

19

CUFFE PKEANE A. Visualizing the electrical structure of power systems[J]. IEEE Systems Journal2017113): 1810-1821. [百度学术] 

20

NELSEN R B. An introduction to copulas[M]. Berlin,GermanySpringer Science & Business Media2007. [百度学术] 

21

AAS KCZADO CFRIGESSI Aet al. Pair-copula constructions of multiple dependence[J]. Insurance: Mathematics and Economics2009442): 182-198. [百度学术] 

22

DIßMANN JBRECHMANN E CCZADO Cet al. Selecting and estimating regular vine copulae and application to financial returns[J]. Computational Statistics & Data Analysis20135952-69. [百度学术] 

23

FANG CXU Y LLI Y L. Optimized C-vine copula and environmental contour of joint wind-wave environment for sea-crossing bridges[J]. Journal of Wind Engineering and Industrial Aerodynamics2022225104989. [百度学术] 

24

QIU Y BLI QPAN Y Ret al. A scenario generation method based on the mixture vine copula and its application in the power system with wind/hydrogen production[J]. International Journal of Hydrogen Energy20194411): 5162-5170. [百度学术] 

25

CZADO CNAGLER T. Vine copula based modeling[J]. Annual Review of Statistics and Its Application20229453-477. [百度学术] 

26

SCHEPSMEIER U. Efficient information based goodness-of-fit tests for vine copula models with fixed margins: a comprehensive review[J]. Journal of Multivariate Analysis201513834-52. [百度学术] 

27

NGUYEN NCARUANA R. Consensus clusterings[C]//Seventh IEEE International Conference on Data Mining (ICDM 2007)October 28-312007Omaha, USA607-612. [百度学术] 

28

WANG YLEBOVITZ LZHENG K Det al. Consensus clustering for bi-objective power network partition[J]. CSEE Journal of Power and Energy Systems202284): 973-982. [百度学术] 

29

林勇陈允鹏王志勇.广东电网目标网架方案论证与建议[J].南方电网技术2020143):42-48. [百度学术] 

LIN YongCHEN YunpengWANG Zhiyonget al. Demonstration and suggestion on network scheme of Guangdong power system[J]. Southern Power System Technology2020143): 42-48. [百度学术] 

30

广东省发展改革委关于印发广东省海上风电发展规划(20172030年)(修编)的通知[EB/OL].(2018-04-11)[2023-06-21].http://drc.gd.gov.cn/gkmlpt/content/1/1060/mmpost_1060661.html#876. [百度学术] 

Issuing of Guangdong Development and Reform Commission on offshore wind power development plan of Guangdong province (from 2017 to 2030) (revised)[EB/OL]. (2018-04-11)[2023-06-21]. http://drc.gd.gov.cn/gkmlpt/content/1/1060/mmpost_1060661.html#876. [百度学术] 

31

HERSBACH HBELL BBERRISFORD Pet al. ERA5 hourly data on single levels from 1940 to present: copernicus climate change service (C3S) climate data store (CDS)[EB/OL]. [2023-07-07]. https://doi.org/10.24381/cds.adbb2d47. [百度学术] 

32

国家能源局.我国海上风电并网容量突破千万千瓦[EB/OL].(2021-05-26)[2023-06-21].http://www.nea.gov.cn/202105/26/c_139970677.Htm. [百度学术] 

National Energy Administration. China’s offshore wind power grid-connected capacity exceeds 10 million kilowatts[EB/OL]. (2021-05-26)[2023-06-21]. http://www.nea.gov.cn/202105/26/c_139970677.Htm. [百度学术]