CN106712105B - 一种新能源接纳空间离散概率序列计算方法 - Google Patents
一种新能源接纳空间离散概率序列计算方法 Download PDFInfo
- Publication number
- CN106712105B CN106712105B CN201510447034.XA CN201510447034A CN106712105B CN 106712105 B CN106712105 B CN 106712105B CN 201510447034 A CN201510447034 A CN 201510447034A CN 106712105 B CN106712105 B CN 106712105B
- Authority
- CN
- China
- Prior art keywords
- load
- max
- thermal power
- day
- new energy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims description 12
- 238000000034 method Methods 0.000 claims abstract description 20
- 230000002354 daily effect Effects 0.000 claims description 92
- 230000003203 everyday effect Effects 0.000 claims description 14
- 238000010248 power generation Methods 0.000 claims description 12
- 238000009434 installation Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- ACWBQPMHZXGDFX-QFIPXVFZSA-N valsartan Chemical compound C1=CC(CN(C(=O)CCCC)[C@@H](C(C)C)C(O)=O)=CC=C1C1=CC=CC=C1C1=NN=NN1 ACWBQPMHZXGDFX-QFIPXVFZSA-N 0.000 claims description 3
- 230000007613 environmental effect Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种新能源接纳空间离散概率序列求解方法。首先制定火电机组开机方式。将新能源出力纳入火电机组开机电力平衡,得到每日的火电机组开机容量和最小技术出力数据。然后,考虑每日最大负荷与火电机组最小技术出力的相关性,求出两者的联合概率分布;同理求出每日最大负荷与负荷的联合概率分布。最后,将相同的日最大负荷离散值对应的负荷序列与火电最小出力序列卷差,求得这个日最大负荷离散值对应的接纳空间序列,依次求出每个日最大负荷离散值对应的接纳空间序列,整理后即得系统的接纳空间离散概率序列。
Description
技术领域
本发明涉及一种新能源并网仿真领域的离散概率序列计算方法,具体涉及一种新能源接纳空间离散概率序列计算方法。
背景技术
近年来,新能源在世界范围内迅猛发展,一些国家已将新能源的开发利用作为改善能源结构、推动环境保护、保持经济和社会可持续发展的重大举措,积极推动新能源的高效利用。目前我国的新能源利用主要是风电和光伏发电。
作为可再生能源中的重要一员,风能已成为发展最快、最具开发潜力的绿色能源。由于风电具有随机性、间歇性等特点,其大规模并网后对电网的稳定性、电能质量和调度机制等方面都带来巨大挑战,使得电网在接纳风电同时必须付出相应的代价。
光伏发电并网在一定程度上缓解了能源危机和环境压力,同时也给电力系统安全运行带来了新的挑战。与火电、水电等常规发电方式不同,由于太阳辐射的波动性和间歇性,光伏电站的出力具有随机性,这使得光伏电站难以调度和控制。除此之外,光伏发电具有天然的昼夜交替特性,夜间光伏电站出力为零,这与风电等其他可再生能源又具有明显的区别。
为保证能源的清洁高效利用,需要对新能源进行优先消纳。然而,新能源在多个时间尺度上具有与生俱来的波动性和随机性,电网消纳新能源的难度随之加大,部分地区已无法全额消纳新能源,必须研究新能源的接纳空间,以实现对新能源消纳电量和限电电量的定量计算。
新能源的接纳空间是系统最多可接纳的新能源电力,为保证尽可能多的接纳新能源,实现新能源的优先消纳,应在保证负荷需求和系统正备用的前提下尽可能少开火电机组,且尽量使火电机组以最小技术出力运行,则新能源的接纳空间就是负荷超出火电机组最小出力的部分。因为负荷和火电机组最小出力是随时间变化的,因此,新能源的接纳空间也是时变的。
发明内容
为解决上述现有技术中的不足,本发明的目的是提供一种新能源接纳空间离散概率序列计算方法。
本发明的目的是采用下述技术方案实现的:
本发明提供一种新能源接纳空间离散概率序列计算方法,其改进之处在于,所述方法基于联合概率分布,包括下述步骤:
步骤1:计算新能源出力参与火电开机的可替代容量;
步骤2:制定火电机组每日开机方式,确定每日火电最小技术出力Pmin;
步骤3:确定每日最大负荷loadmax与每日火电最小技术出力Pmin的联合概率分布;
步骤4:确定每日最大负荷loadmax与负荷load的联合概率分布;
步骤5:用卷差确定新能源接纳空间概率序列。
进一步地,所述步骤1包括下述步骤:
步骤1.1:获取并整理数据:获取一定时间周期的负荷、新能源(风电和光伏发电)出力和装机时序数据,数据总量都是T,研究周期的总天数为D,设负荷、新能源出力和装机每天的数据量为N,则T=N·D;若负荷等时序数据的采样间隔为1h,则N=24;若采样间隔为15min,则N=96;
步骤1.2:求解新能源参与火电开机的可替代容量:新能源的预测误差为20%,将新能源预测值减去其装机的20%参与火电开机;若新能源出力预测值大于其装机容量的20%,则超出部分纳入火电开机电力平衡,减少火电开机容量;否则,新能源不参与火电机组开机,则研究周期T内新能源参与火电开机的可替代容量;Crenewableenergy(t)满足下式:
Crenewableenergy(t)=max(Prenewableenergy(t)-Capacityrenewableenergy(t)×20%,0)
Prenewableenergy(t)=Pwind(t)+PPV(t)
Capacityrenewableenergy(t)=Capacitywind(t)+CapacityPV(t)
其中,t表示研究周期内的第t个时段,0≤t≤T,Prenewableenergy(t)是t时段的新能源出力预测值,Capacityrenewableenergy(t)是t时段的新能源装机容量,Pwind(t)是t时段的风力发电预测值,Capacitywind(t)是t时段的风电装机容量,PPV(t)是t时段的光伏发电预测值,CapacityPV(t)是t时段的光伏装机容量。
进一步地,所述步骤2中包括以下步骤,
步骤2.1:求等效负荷eqload和日最大等效负荷eqloadmax:Crenewableenergy(t)参与火电机组开机电力平衡后的等效负荷满足下式:
eqload(t)=load(t)-Crenewableenergy(t)
其中eqload(t)是t时段的等效负荷,load(t)是t时段的负荷,Crenewableenergy(t)是t时段的新能源可替代容量;
eqload的数据总量也为T,每天的数据量为N,求出eqload每天的最大值,即为日最大等效负荷eqloadmax,eqloadmax的数据总量是D;
步骤2.2:制定火电机组每日开机方式,求出每日火电最小技术出力:
根据每日最大等效负荷和系统正备用,结合系统中的火电机组参数,可依次求出每天的开机方式,包括:
将电力系统中的必开火电机组开机,供热期为保证供热,有必开机组,非供热期无必开机组;
若等效负荷的当日最大值加上电力系统正备用大于必开机组的最大出力,则增开其余机组;为使新能源接纳电量尽可能大,在最大出力相同的情况下,最小出力小的机组优先开机,直至电力系统中火电机组的最大出力满足日最大等效负荷和系统正备用,停止增开;
所述火电机组每日开机方式用混合整数规划模型表示的目标函数为:
minPmin(d)
其中,Pmin(d)表示火电机组第d天的最小出力功率;
火电机组开机方式混合整数规划模型的约束条件为:
(1)火电机组第d天的最小出力等于第d天开启的所有火电机组最小出力之和;
其中,Xj(d)示第j台机组第d天的运行状态,为二进制变量,0表示机组停机,1则表示机组正在运行;TPj,min(d)是第j台火电机组第d天的最小出力功率;
(2)火电机组每日开机的最大出力必须保证当日的负荷需求和系统备用,即:
其中,TPj,max(d)是第j台火电机组第d天的最大出力功率,eqloadmax(d)是第d天的日最大等效负荷,Pre是系统的正备用;
(3)每种火电机组每天的的启机台数在其最大台数和最小台数之间,即:
Sj,min(d)≤Sj(d)≤Sj,max(d)
其中,Sj(d)是第d天的第j台火电机组的启机台数,Sj,min(d)是第d天第j台火电机组的最小启机台数,Sj,max(d)是第d天第j台火电机组的最大启机台数;
根据所述模型依次求出火电机组每日开机方式,进而求出研究周期的火电机组日最小技术出力Pmin,Pmin的数据总量为D。
进一步地,所述步骤3包括下述步骤:
步骤3.1:将每日最大负荷loadmax与每日火电最小技术出力Pmin分别在其变化范围内分区:loadmax和Pmin的数据总量都为D;选取每日最大负荷loadmax的离散化区间长度为C1,每日火电最小技术出力Pmin的离散化区间长度为C2,loadmax和Pmin的区间个数分别为L1和L2;
步骤3.2:设置循环初值:d=1,nik=0,其中,d表示第d天,0≤d≤D,nik为loadmax(d)归入其第i个离散化区间,且Pmin(d)归入其第k个离散化区间的概数,i=1,2,...,L1,k=1,2,...,L2;
步骤3.3:读入第d天的每日最大负荷loadmax(d)和每日火电最小技术出力Pmin(d),并归入各自对应的离散化区间;
步骤3.4:d=d+1,nik=nik+1;
步骤3.5:若d≤D,则转至步骤3.3,否则,转至步骤3.6;
步骤3.6:确定每日最大负荷loadmax每个离散化区间内所有数据的平均值averageloadmax(i),i=1,2...L1;确定每日火电最小技术出力Pmin每个区间内所有数据的平均值averagePmin(k),k=1,2...L2;
步骤3.7:确定loadmax归入其第i个离散化区间,且Pmin归入其第k个离散化区间的概率Pik,loadmax和Pmin归入的离散化区间组合共有L1·L2个,依次求出每种离散化区间组合的概率,进而得到日最大负荷与日火电最小出力的联合概率分布。
进一步地,所述步骤4包括:
步骤4.1:将每日最大负荷loadmax与负荷load分别在其变化范围内分区:loadmax的数据总量为D,load的数据总量为T,选取loadmax的离散化区间长度为C1,选取负荷load的离散化区间长度为C3,loadmax和load区间个数分别为L1和L3;
步骤4.2:设置循环初值,d=1,nis=0,其中,d表示第d天,nis为loadmax(d)归入其第i个离散化区间,且t时段的负荷load(t)归入其第s个离散化区间的概数,i=1,2,...,L1,s=1,2,...,L3;
步骤4.3:读入第d天的每日最大负荷loadmax(d),归入其所对应的区间;
步骤4.4:设置初值,t=(d-1)·N+1,此时,t表示第d天的第1个时段。
步骤4.5:读入t时段的负荷load(t),并归入其对应的离散化区间;
步骤4.6:t=t+1,nis=nis+1;
步骤4.7:若t≤d·N,则转至步骤4.5,否则,转至步骤4.8;
步骤4.8:d=d+1;
步骤4.9:若d≤D,则转至步骤4.3,否则,转至步骤4.10;
步骤4.10:确定每日最大负荷loadmax每个离散化区间内所有数据的平均值averageloadmax(i),i=1,2...L1;确定负荷load每个区间内所有数据的平均值averageload(s),s=1,2...L3;
步骤4.11:确定每日最大负荷loadmax归入其第i个离散化区间,且每日火电最小技术出力Pmin归入其第s个离散化区间的概率Pis,loadmax和load归入的离散化区间组合共有L1·L3个,依次求出每种区间组合的概率,进而得到日最大负荷与负荷的联合概率分布。
进一步地,所述步骤5包括:日最大负荷loadmax每个离散化区间内所有数据的平均值averageloadmax(i)均对应每日火电最小技术出力的离散概率序列gPmini,gPmini的第一行是离散功率值,第二行是对应的概率,即:
gPmini(1,k)=averagePmin(k)
每日最大负荷loadmax每个离散化区间内所有数据的平均值averageloadmax(i)对应的负荷离散概率序列gloadi满足下式:
gloadi(1,s)=averageload(s)
对gPmini和gloadi,每日火电机组最小技术出力离散功率值与负荷离散功率值所有的组合数为L2·L3,若某一种组合中,负荷离散功率值是gloadi(1,s),每日火电机组最小技术出力离散功率值是gPmini(1,k),则此组合的新能源接纳空间功率gACCOMi(1,j)满足下式:
gACCOMi(1,j)=gloadi(1,s)-gPmini(1,k)
此组合的概率为:
gACCOMi(2,j)=gloadi(2,s)·gPmini(2,k)
确定所有组合的新能源接纳空间功率值及其概率,得到averageloadmax(i)对应的新能源接纳空间离散概率序列gACCOMi;
上述为将gloadi与gPmini做卷差运算,表示为,
由于gACCOMi是averageloadmax(i)对应的接纳空间概率序列,在整个概率空间中,当loadmax归入离散值averageloadmax(i)时的接纳空间离散概率序列gACCOM'i满足下式,
gACCOM'i(1,j)=gACCOMi(1,j)
本发明提供的技术方案具有的优异效果是:
(1)将新能源出力纳入火电机组开机电力平衡,可减少火电机组开机容量,能更好地消纳新能源。
(2)与采用负荷预测时序数据和对应的火电最小出力时序数据直接相减的方法求解接纳空间的方法不同,采用概率序列运算的方法求解接纳空间,减小了负荷和新能源时序预测数据的偶然性对接纳空间计算结果的影响,能更好地反映接纳空间的规律性。
(3)充分考虑日最大负荷与火电机组最小技术出力的相关性,日最大负荷与负荷的相关性,将相同日最大负荷离散值对应的负荷概率序列与火电最小出力概率序列卷差,求得日最大负荷离散值对应的接纳空间序列。避免了使用负荷和火电最小出力时序数据转换的离散概率序列直接卷差造成的误差。
附图说明
图1是本发明提供的新能源接纳空间离散概率序列计算方法的流程图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
以下描述和附图充分地示出本发明的具体实施方案,以使本领域的技术人员能够实践它们。其他实施方案可以包括结构的、逻辑的、电气的、过程的以及其他的改变。实施例仅代表可能的变化。除非明确要求,否则单独的组件和功能是可选的,并且操作的顺序可以变化。一些实施方案的部分和特征可以被包括在或替换其他实施方案的部分和特征。本发明的实施方案的范围包括权利要求书的整个范围,以及权利要求书的所有可获得的等同物。在本文中,本发明的这些实施方案可以被单独地或总地用术语“发明”来表示,这仅仅是为了方便,并且如果事实上公开了超过一个的发明,不是要自动地限制该应用的范围为任何单个发明或发明构思。
本发明提供了一种基于联合概率分布的新能源接纳空间离散概率序列求解方法,由5个步骤组成,其流程图如图1所示:
步骤1:求解新能源出力参与火电开机的可替代容量。
步骤1.1:整理并获取数据:获取一定时间周期(一般为一年或几个月)的负荷、新能源(风电和光伏发电)出力和装机时序数据,数据总量都是T,研究周期的总天数为D,设负荷、新能源出力和装机每天的数据量为N,T=N·D。若负荷等时序数据的采样间隔为1h,则N=24;若采样间隔为15min,则N=96。
步骤1.2:求解新能源参与火电开机的可替代容量:为了尽量多接纳新能源,将各时段新能源发电出力纳入火电机组开机容量计算范畴,可替代一部分火电,减小火电开机容量。目前,新能源的预测误差约为20%,因此,将新能源预测值减去其装机的20%参与火电开机。若新能源出力预测值大于其装机容量的20%,则超出部分纳入火电开机电力平衡,减少火电开机容量;否则,新能源不参与火电机组开机。则研究周期T内新能源参与火电开机的可替代容量Crenewableenergy(t)满足下式:
Crenewableenergy(t)=max(Prenewableenergy(t)-Capacityrenewableenergy(t)×20%,0)
Prenewableenergy(t)=Pwind(t)+PPV(t)
Capacityrenewableenergy(t)=Capacitywind(t)+CapacityPV(t)
其中,t表示研究周期内的第t个时段,0≤t≤T,Prenewableenergy(t)是t时段的新能源出力预测值,Capacityrenewableenergy(t)是t时段的新能源装机容量,Pwind(t)是t时段的风力发电预测值,Capacitywind(t)是t时段的风电装机容量,PPV(t)是t时段的光伏发电预测值,CapacityPV(t)是t时段的光伏装机容量。
步骤2:制定火电机组每日开机方式,确定每日火电最小技术出力Pmin。
步骤2.1:求等效负荷eqload和日最大等效负荷eqloadmax:Crenewableenergy(t)参与火电机组开机电力平衡后的等效负荷满足下式:
eqload(t)=load(t)-Crenewableenergy(t)
其中eqload(t)是t时段的等效负荷,load(t)是t时段的负荷,Crenewableenergy(t)是由步骤1求出的t时段的新能源可替代容量。
eqload的数据总量也是T,每天的数据量也是N,求出eqload每天的最大值,就是日最大等效负荷eqloadmax,eqloadmax的数据总量是D。
步骤2.2:制定火电机组每日开机方式,求出每日火电最小技术出力。
根据每日最大等效负荷和系统正备用,结合系统中的火电机组参数,可依次求出每天的开机方式。具体方法是:
将电力系统中的必开火电机组开机,供热期为保证供热,有必开机组,非供热期无必开机组;
若等效负荷的当日最大值加上电力系统正备用大于必开机组的最大出力,则增开其余机组;为使新能源接纳电量尽可能大,在最大出力相同的情况下,最小出力小的机组优先开机,直至电力系统中火电机组的最大出力满足日最大等效负荷和系统正备用,停止增开。
上述的日开机方式用混合整数规划模型表示的目标函数为:
minPmin(d)
其中,Pmin(d)表示火电机组第d天的最小出力功率。上式表示在制定每日的火电机组开机计划时,应使每日开启的火电机组最小出力尽可能小。
火电机组开机方式混合整数规划模型的约束条件为:
(1)火电机组第d天的最小出力等于第d天开启的所有火电机组最小出力之和。
其中,Xj(d)示第j台机组第d天的运行状态,为二进制变量,0表示机组停机,1则表示机组正在运行。TPj,min(d)是第j台火电机组第d天的最小出力功率。
(3)火电机组每日开机的最大出力必须保证当日的负荷需求和系统备用,
其中,TPj,max(d)是第j台火电机组第d天的最大出力功率,eqloadmax(d)是第d天的日最大等效负荷,Pre是系统的正备用。
(4)每种火电机组每天的的启机台数应在其最大台数和最小台数之间,即
Sj,min(d)≤Sj(d)≤Sj,max(d)
其中,Sj(d)是第d天的第j台火电机组的启机台数,Sj,min(d)是第d天第j台火电机组的最小启机台数,Sj,max(d)是第d天第j台火电机组的最大启机台数。
根据步骤2所述方法,依次求出每天的开机方式,就可求出研究周期的火电机组日最小技术出力Pmin,Pmin的数据总量为D。
步骤3:求每日最大负荷loadmax与Pmin的联合概率分布。
由于火电最小技术出力Pmin由每日最大等效负荷与系统正备用之和决定,而系统的正备用一般是定值,每日最大等效负荷又与每日最大负荷存在较强的相关性,因此,需通过联合概率分布的方法求取每日最大负荷loadmax与Pmin的联合概率分布。
步骤3.1:将每日最大负荷loadmax与每日火电最小技术出力Pmin分别在其变化范围内分区:loadmax和Pmin的数据总量都为D。选取每日最大负荷loadmax的离散化区间长度为C1,每日火电最小技术出力Pmin的离散化区间长度为C2,loadmax和Pmin的区间个数分别为L1、L2;
步骤3.2:设置循环初值:d=1,nik=0,其中,d表示第d天,0≤d≤D,nik为loadmax(d)归入其第i个离散化区间,且Pmin(d)归入其第k个离散化区间的概数,i=1,2,...,L1,k=1,2,...,L2;
步骤3.3:读入第d天的每日最大负荷loadmax(d)和每日火电最小技术出力Pmin(d),并归入各自对应的离散化区间;
步骤3.4:d=d+1,nik=nik+1;
步骤3.5:若d≤D,则转至步骤3.3,否则,转至步骤3.6;
步骤3.6:确定每日最大负荷loadmax每个离散化区间内所有数据的平均值averageloadmax(i),i=1,2...L1;确定每日火电最小技术出力Pmin每个区间内所有数据的平均值averagePmin(k),k=1,2...L2;
步骤3.7:确定loadmax归入其第i个离散化区间,且Pmin归入其第k个离散化区间的概率Pik,loadmax和Pmin归入的离散化区间组合共有L1·L2个,依次求出每种区间组合的概率,就得到了日最大负荷与日火电最小出力的联合概率分布。
求得的loadmax与Pmin的联合概率分布a2×1如表1所示。
表1 loadmax与Pmin的联合概率分布表
步骤4:求loadmax与负荷load的联合概率分布:
loadmax与load也存在较强相关性,也需通过联合概率分布的方法求取每日最大负荷loadmax与load的联合概率分布。
步骤4.1:将每日最大负荷loadmax与负荷load分别在其变化范围内分区:loadmax的数据总量为D,load的数据总量为T,每日最大负荷loadmax与步骤3.1中的分区方法相同,选取loadmax的离散化区间长度为C1,选取负荷load的离散化区间长度为C3,loadmax和load区间个数分别为L1和L3;
步骤4.2:设置循环初值,d=1,nis=0,其中,d表示第d天,nis为loadmax(d)归入其第i个离散化区间,且t时段的负荷load(t)归入其第s个离散化区间的概数,i=1,2,...,L1,s=1,2,...,L3;
步骤4.3:读入第d天的每日最大负荷loadmax(d),归入其所对应的区间;
步骤4.4:设置初值,t=(d-1)·N+1,此时,t表示第d天的第1个时段。
步骤4.5:读入t时段的负荷load(t),并归入其对应的离散化区间;
步骤4.6:t=t+1,nis=nis+1;
步骤4.7:若t≤d·N,则转至步骤4.5,否则,转至步骤4.8;
步骤4.8:d=d+1;
步骤4.9:若d≤D,则转至步骤4.3,否则,转至步骤4.10;
步骤4.10:确定每日最大负荷loadmax每个离散化区间内所有数据的平均值averageloadmax(i),i=1,2...L1;确定负荷load每个区间内所有数据的平均值averageload(s),s=1,2...L3;
步骤4.11:确定每日最大负荷loadmax归入其第i个离散化区间,且每日火电最小技术出力Pmin归入其第s个离散化区间的概率Pis,loadmax和load归入的离散化区间组合共有L1·L3个,依次求出每种区间组合的概率,进而得到日最大负荷与负荷的联合概率分布。求出loadmax与load的联合概率分布,如表2所示。
表2 loadmax与load的联合概率分布表
步骤5:卷差求接纳空间概率序列。
由表2可知,每个averageloadmax(i)都对应一个火电机组日最小出力的离散概率序列gPmini,gPmini的第一行是离散功率值,第二行是对应的概率,即,
gPmini(1,k)=averagePmin(k)
同理,每个averageloadmax(i)对应的负荷离散概率序列gloadi满足下式,
gloadi(1,s)=averageload(s)
对gPmini和gloadi,火电机组日最小技术出力离散功率值与负荷离散功率值所有可能的组合数为L2·L3,若某一种组合中,负荷离散功率值是gloadi(1,s),火电日最小出力离散功率值是gPmini(1,k),则此组合的新能源接纳空间功率gACCOMi(1,j)满足下式:
gACCOMi(1,j)=gloadi(1,s)-gPmini(1,k)
此组合的概率为:
gACCOMi(2,j)=gloadi(2,s)·gPmini(2,k)
求出所有组合的新能源接纳空间功率值及其概率,就可得到averageloadmax(i)对应的新能源接纳空间离散概率序列gACCOMi。
以上其实是将gloadi与gPmini做卷差运算,可表示为,
又由于gACCOMi是averageloadmax(i)对应的接纳空间概率序列,所以在整个概率空间中,当loadmax归入离散值averageloadmax(i)时的接纳空间离散概率序列gACCOM'i满足下式,
gACCOM'i(1,j)=gACCOMi(1,j)
以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员依然可以对本发明的具体实施方式进行修改或者等同替换,这些未脱离本发明精神和范围的任何修改或者等同替换,均在申请待批的本发明的权利要求保护范围之内。
Claims (6)
1.一种新能源接纳空间离散概率序列计算方法,其特征在于,所述方法基于联合概率分布,包括下述步骤:
步骤1:计算新能源出力参与火电开机的可替代容量;
步骤2:制定火电机组每日开机方式,确定每日火电最小技术出力Pmin;
步骤3:确定每日最大负荷loadmax与每日火电最小技术出力Pmin的联合概率分布;
步骤4:确定每日最大负荷loadmax与负荷load的联合概率分布;
步骤5:用卷差确定新能源接纳空间概率序列。
2.如权利要求1所述的新能源接纳空间离散概率序列计算方法,其特征在于,所述步骤1包括下述步骤:
步骤1.1:获取并整理数据:获取一定时间周期的负荷、新能源出力和装机时序数据,数据总量都是T,研究周期的总天数为D,设负荷、新能源出力和装机每天的数据量为N,则T=N·D;若负荷等时序数据的采样间隔为1h,则N=24;若采样间隔为15min,则N=96;
步骤1.2:求解新能源参与火电开机的可替代容量:新能源的预测误差为20%,将新能源预测值减去其装机的20%参与火电开机;若新能源出力预测值大于其装机容量的20%,则超出部分纳入火电开机电力平衡,减少火电开机容量;否则,新能源不参与火电机组开机,则研究周期T内新能源参与火电开机的可替代容量Crenewableenergy(t)满足下式:
Crenewableenergy(t)=max(Prenewableenergy(t)-Capacityrenewableenergy(t)×20%,0)
Prenewableenergy(t)=Pwind(t)+PPV(t)
Capacityrenewableenergy(t)=Capacitywind(t)+CapacityPV(t)
其中,t表示研究周期内的第t个时段,0≤t≤T,Prenewableenergy(t)是t时段的新能源出力预测值,Capacityrenewableenergy(t)是t时段的新能源装机容量,Pwind(t)是t时段的风力发电预测值,Capacitywind(t)是t时段的风电装机容量,PPV(t)是t时段的光伏发电预测值,CapacityPV(t)是t时段的光伏装机容量。
3.如权利要求1所述的新能源接纳空间离散概率序列计算方法,其特征在于,所述步骤2中包括以下步骤,
步骤2.1:求等效负荷eqload和日最大等效负荷eqloadmax:Crenewableenergy(t)参与火电机组开机电力平衡后的等效负荷满足下式:
eqload(t)=load(t)-Crenewableenergy(t)
其中eqload(t)是t时段的等效负荷,load(t)是t时段的负荷,Crenewableenergy(t)是t时段的新能源可替代容量;
eqload的数据总量也为T,每天的数据量为N,求出eqload每天的最大值,即为日最大等效负荷eqloadmax,eqloadmax的数据总量是D;
步骤2.2:制定火电机组每日开机方式,求出每日火电最小技术出力:
根据每日最大等效负荷和系统正备用,结合系统中的火电机组参数,可依次求出每天的开机方式,包括:
将电力系统中的必开火电机组开机,供热期为保证供热,有必开机组,非供热期无必开机组;
若等效负荷的当日最大值加上电力系统正备用大于必开机组的最大出力,则增开其余机组;为使新能源接纳电量尽可能大,在最大出力相同的情况下,最小出力小的机组优先开机,直至电力系统中火电机组的最大出力满足日最大等效负荷和系统正备用,停止增开;
所述火电机组每日开机方式用混合整数规划模型表示的目标函数为:
minPmin(d)
其中,Pmin(d)表示火电机组第d天的最小出力功率;
火电机组开机方式混合整数规划模型的约束条件为:
(1)火电机组第d天的最小出力等于第d天开启的所有火电机组最小出力之和;
其中,Xj(d)示第j台机组第d天的运行状态,为二进制变量,0表示机组停机,1则表示机组正在运行;TPj,min(d)是第j台火电机组第d天的最小出力功率;
(2)火电机组每日开机的最大出力必须保证当日的负荷需求和系统备用,即:
其中,TPj,max(d)是第j台火电机组第d天的最大出力功率,eqloadmax(d)是第d天的日最大等效负荷,Pre是系统的正备用;
(3)每种火电机组每天的的启机台数在其最大台数和最小台数之间,即:
Sj,min(d)≤Sj(d)≤Sj,max(d)
其中,Sj(d)是第d天的第j台火电机组的启机台数,Sj,min(d)是第d天第j台火电机组的最小启机台数,Sj,max(d)是第d天第j台火电机组的最大启机台数;
根据所述模型依次求出火电机组每日开机方式,进而求出研究周期的火电机组日最小技术出力Pmin,Pmin的数据总量为D。
4.如权利要求1所述的新能源接纳空间离散概率序列计算方法,其特征在于,所述步骤3包括下述步骤:
步骤3.1:将每日最大负荷loadmax与每日火电最小技术出力Pmin分别在其变化范围内分区:loadmax和Pmin的数据总量都为D;选取每日最大负荷loadmax的离散化区间长度为C1,每日火电最小技术出力Pmin的离散化区间长度为C2,loadmax和Pmin的区间个数分别为L1和L2;
步骤3.2:设置循环初值:d=1,nik=0,其中,d表示第d天,0≤d≤D,nik为loadmax(d)归入其第i个离散化区间,且Pmin(d)归入其第k个离散化区间的概数,i=1,2,...,L1,k=1,2,...,L2;
步骤3.3:读入第d天的每日最大负荷loadmax(d)和每日火电最小技术出力Pmin(d),并归入各自对应的离散化区间;
步骤3.4:d=d+1,nik=nik+1;
步骤3.5:若d≤D,则转至步骤3.3,否则,转至步骤3.6;
步骤3.6:确定每日最大负荷loadmax每个离散化区间内所有数据的平均值averageloadmax(i),i=1,2...L1;确定每日火电最小技术出力Pmin每个区间内所有数据的平均值averagePmin(k),k=1,2...L2;
5.如权利要求1所述的新能源接纳空间离散概率序列计算方法,其特征在于,所述步骤4包括:
步骤4.1:将每日最大负荷loadmax与负荷load分别在其变化范围内分区:loadmax的数据总量为D,load的数据总量为T,选取loadmax的离散化区间长度为C1,选取负荷load的离散化区间长度为C3,loadmax和load区间个数分别为L1和L3;
步骤4.2:设置循环初值,d=1,nis=0,其中,d表示第d天,nis为loadmax(d)归入其第i个离散化区间,且t时段的负荷load(t)归入其第s个离散化区间的概数,i=1,2,...,L1,s=1,2,...,L3;
步骤4.3:读入第d天的每日最大负荷loadmax(d),归入其所对应的区间;
步骤4.4:设置初值,t=(d-1)·N+1,此时,t表示第d天的第1个时段;
步骤4.5:读入t时段的负荷load(t),并归入其对应的离散化区间;
步骤4.6:t=t+1,nis=nis+1;
步骤4.7:若t≤d·N,则转至步骤4.5,否则,转至步骤4.8;
步骤4.8:d=d+1;
步骤4.9:若d≤D,则转至步骤4.3,否则,转至步骤4.10;
步骤4.10:确定每日最大负荷loadmax每个离散化区间内所有数据的平均值averageloadmax(i),i=1,2...L1;确定负荷load每个区间内所有数据的平均值averageload(s),s=1,2...L3;
6.如权利要求1所述的新能源接纳空间离散概率序列计算方法,其特征在于,所述步骤5包括:日最大负荷loadmax每个离散化区间内所有数据的平均值averageloadmax(i)均对应每日火电最小技术出力的离散概率序列gPmini,gPmini的第一行是离散功率值,第二行是对应的概率,即:
gPmini(1,k)=averagePmin(k)
每日最大负荷loadmax每个离散化区间内所有数据的平均值averageloadmax(i)对应的负荷离散概率序列gloadi满足下式:
gloadi(1,s)=averageload(s)
对gPmini和gloadi,每日火电机组最小技术出力离散功率值与负荷离散功率值所有的组合数为L2·L3,若某一种组合中,负荷离散功率值是gloadi(1,s),每日火电机组最小技术出力离散功率值是gPmini(1,k),则此组合的新能源接纳空间功率gACCOMi(1,j)满足下式:
gACCOMi(1,j)=gloadi(1,s)-gPmini(1,k)
此组合的概率为:
gACCOMi(2,j)=gloadi(2,s)·gPmini(2,k)
确定所有组合的新能源接纳空间功率值及其概率,得到averageloadmax(i)对应的新能源接纳空间离散概率序列gACCOMi;
上述为将gloadi与gPmini做卷差运算,表示为,
由于gACCOMi是averageloadmax(i)对应的接纳空间概率序列,在整个概率空间中,当loadmax归入离散值averageloadmax(i)时的接纳空间离散概率序列gACCOM'i满足下式,
gACCOM'i(1,j)=gACCOMi(1,j)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510447034.XA CN106712105B (zh) | 2015-07-27 | 2015-07-27 | 一种新能源接纳空间离散概率序列计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510447034.XA CN106712105B (zh) | 2015-07-27 | 2015-07-27 | 一种新能源接纳空间离散概率序列计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106712105A CN106712105A (zh) | 2017-05-24 |
CN106712105B true CN106712105B (zh) | 2021-01-15 |
Family
ID=58895350
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510447034.XA Active CN106712105B (zh) | 2015-07-27 | 2015-07-27 | 一种新能源接纳空间离散概率序列计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106712105B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106410780A (zh) * | 2015-07-28 | 2017-02-15 | 中国电力科学研究院 | 一种新能源接纳空间离散概率序列计算方法 |
CN107342592B (zh) * | 2017-08-31 | 2022-05-13 | 国电南瑞科技股份有限公司 | 一种基于故障触发的紧急协调电负荷平衡和快速切负荷措施的安全稳定控制方法 |
CN108805388A (zh) * | 2018-04-09 | 2018-11-13 | 中国电力科学研究院有限公司 | 一种未来年负荷时间序列场景的确定方法和装置 |
CN112134273A (zh) * | 2020-08-06 | 2020-12-25 | 国家电网公司西北分部 | 一种基于风险控制的新能源纳入电力系统备用方法及装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8285719B1 (en) * | 2008-08-08 | 2012-10-09 | The Research Foundation Of State University Of New York | System and method for probabilistic relational clustering |
CN104167765B (zh) * | 2014-07-11 | 2016-08-24 | 海南电网公司 | 一种基于接纳能力分布的最大风电装机容量计算方法 |
CN104463378B (zh) * | 2014-12-29 | 2017-08-29 | 国家电网公司 | 用于水电资源丰富地区省级电网风力发电消纳能力的算法 |
CN104617591B (zh) * | 2015-01-28 | 2017-01-25 | 西安交通大学 | 基于多场景新能源发电模拟的日运行方式安排和调峰方法 |
CN106099993B (zh) * | 2016-06-28 | 2018-06-08 | 中国电力工程顾问集团西南电力设计院有限公司 | 一种适应新能源大规模接入的电源规划方法 |
-
2015
- 2015-07-27 CN CN201510447034.XA patent/CN106712105B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106712105A (zh) | 2017-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109103926A (zh) | 基于多辐照特性年气象场景的光伏发电接纳能力计算方法 | |
CN106532764B (zh) | 一种就地消纳光伏发电的电动汽车充电负荷调控方法 | |
CN113095791B (zh) | 一种综合能源系统运行方法及系统 | |
Piacentino et al. | Sustainable and cost-efficient energy supply and utilisation through innovative concepts and technologies at regional, urban and single-user scales | |
Ceran | The concept of use of PV/WT/FC hybrid power generation system for smoothing the energy profile of the consumer | |
CN109858673A (zh) | 一种光伏发电系统功率预测方法 | |
Zhang et al. | Economical assessment of large-scale photovoltaic power development in China | |
Castañeda et al. | Sizing methods for stand-alone hybrid systems based on renewable energies and hydrogen | |
CN106712105B (zh) | 一种新能源接纳空间离散概率序列计算方法 | |
Xu et al. | The production efficiency of renewable energy generation and its influencing factors: Evidence from 20 countries | |
CN104268659A (zh) | 一种光伏电站发电功率超短期预测方法 | |
Tan et al. | A multi-objective optimization dispatching and adaptability analysis model for wind-PV-thermal-coordinated operations considering comprehensive forecasting error distribution | |
Tan et al. | Photovoltaic power generation in China: Development potential, benefits of energy conservation and emission reduction | |
CN105354637A (zh) | 一种使用弃风电加热的电采暖储热运行计划的制定方法 | |
CN107392791B (zh) | 多能互补系统的分布式光伏与气电混合容量规划方法和系统 | |
CN112134271A (zh) | 一种多能互补系统装机容量优化方法及系统 | |
CN107181272A (zh) | 一种利用储能系统提升调峰空间的风电消纳方法 | |
Ma et al. | Long-term coordination for hydro-thermal-wind-solar hybrid energy system of provincial power grid | |
CN107622331B (zh) | 一种发电机组与电力用户直接交易方式的优化方法和装置 | |
Zhao et al. | Optimal configuration of grid connected microgrid considering CCHP and analysis of energy saving and emission reduction | |
Mahoor et al. | Smart energy management for a micro-grid with consideration of demand response plans | |
Niaz et al. | Environment Friendly Hybrid Solar-Hydro Power Distribution Scheduling on Demand Side | |
CN105634005B (zh) | 一种用于评估光伏接纳能力的方法及系统 | |
CN114530882A (zh) | 一种多分布式能源的供用能协同自治方法及系统 | |
CN112580897A (zh) | 基于鹦鹉算法的多能源电力系统优化调度的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15 Applicant after: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd. Applicant after: STATE GRID CORPORATION OF CHINA Address before: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15 Applicant before: China Electric Power Research Institute Applicant before: State Grid Corporation of China |
|
GR01 | Patent grant | ||
GR01 | Patent grant |