CN108039739B - A dynamic stochastic economic dispatch method for active distribution network - Google Patents

A dynamic stochastic economic dispatch method for active distribution network Download PDF

Info

Publication number
CN108039739B
CN108039739B CN201711205188.3A CN201711205188A CN108039739B CN 108039739 B CN108039739 B CN 108039739B CN 201711205188 A CN201711205188 A CN 201711205188A CN 108039739 B CN108039739 B CN 108039739B
Authority
CN
China
Prior art keywords
power
distribution network
period
constraints
node
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
Application number
CN201711205188.3A
Other languages
Chinese (zh)
Other versions
CN108039739A (en
Inventor
吴素农
吴文传
杨为群
张伯明
熊宁
于金镒
栗子豪
孙宏斌
李迎军
徐俊杰
李伟伟
翟亚军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Jiangxi Electric Power Co
Tsinghua University
State Grid Corp of China SGCC
Economic and Technological Research Institute of State Grid Jiangxi Electric Power Co Ltd
Original Assignee
State Grid Jiangxi Electric Power Co
Tsinghua University
State Grid Corp of China SGCC
Economic and Technological Research Institute of State Grid Jiangxi Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Jiangxi Electric Power Co, Tsinghua University, State Grid Corp of China SGCC, Economic and Technological Research Institute of State Grid Jiangxi Electric Power Co Ltd filed Critical State Grid Jiangxi Electric Power Co
Priority to CN201711205188.3A priority Critical patent/CN108039739B/en
Publication of CN108039739A publication Critical patent/CN108039739A/en
Application granted granted Critical
Publication of CN108039739B publication Critical patent/CN108039739B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for AC mains or AC distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48—Controlling the sharing of the in-phase component
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)

Abstract

本发明提出一种主动配电网动态随机经济调度方法，属于电力系统调度与控制领域。该方法首先建立由目标函数和约束条件构成的主动配电网动态随机经济调度模型；然后，对模型的约束条件进行转化，收集配电网中所有节点在每个时段的有功负荷的预测误差值集合和分布式电源功率预测误差值集合，并根据统计信息分别构建对应的不确定量的概率分布集合，构建包含旋转备用约束的机会约束，并利用凸松弛将其转化为确定性线性约束；最后，应用凸规划算法对模型求解，得到各发电机组最优经济调度功率计划。本发明考虑功率随机性，可以用于解决功率预测不确定性的主动配电网经济调度问题，调度策略具备鲁棒性和可靠性，有较高的应用价值。The invention proposes a dynamic random economic dispatch method for an active distribution network, which belongs to the field of power system dispatch and control. The method firstly establishes a dynamic stochastic economic dispatch model of active distribution network composed of objective function and constraints; then, the constraints of the model are transformed, and the prediction error value of active load of all nodes in the distribution network in each period is collected. Sets and sets of power prediction error values of distributed power sources, and respectively constructs the probability distribution sets of corresponding uncertain quantities according to statistical information, constructs chance constraints including spinning reserve constraints, and uses convex relaxation to convert them into deterministic linear constraints; finally , the convex programming algorithm is used to solve the model, and the optimal economic dispatch power plan for each generator set is obtained. The invention considers the randomness of power, and can be used to solve the economical dispatching problem of active distribution network with uncertainty of power prediction. The dispatching strategy has robustness and reliability, and has high application value.

Description

一种主动配电网动态随机经济调度方法A dynamic stochastic economic dispatch method for active distribution network

技术领域technical field

本发明属于电力系统调度与控制技术领域，特别涉及一种主动配电网动态随机经济调度方法。The invention belongs to the technical field of power system scheduling and control, in particular to a dynamic random economic scheduling method for an active distribution network.

背景技术Background technique

为了应对以光伏、风电为主的分布式电源在主动配电网中大规模接入所带来的技术问题，降低主动配电网整体全体运行发电成本，需要对主动配电网进行动态经济调度，并制定出最优的发电机组功率方案，以实现主动配电网中高效经济运行的目标。In order to cope with the technical problems brought about by the large-scale access of distributed power sources mainly photovoltaic and wind power in the active distribution network, and reduce the overall operating power generation cost of the active distribution network, it is necessary to carry out dynamic economic dispatching of the active distribution network. , and formulate the optimal generator set power scheme to achieve the goal of efficient and economical operation in the active distribution network.

由于分布式电源功率，包括全天分布式电源的有功、无功功率，受天气和环境因素的影响而具有显著的波动性和间歇性，现有的预测技术无法对分布式电源未来功率进行精准预测；同样的，现有预测技术也无法对配电网中的节点负荷进行准确预测。因此，分布式电源功率和负荷预测误差，为主动配电网中的动态经济调度问题引入了很强的不确定性。Because the power of distributed power, including the active and reactive power of distributed power throughout the day, is affected by weather and environmental factors and has significant fluctuations and intermittence, the existing prediction technology cannot accurately predict the future power of distributed power. Prediction; Similarly, the existing prediction technology cannot accurately predict the node load in the distribution network. Therefore, the power and load prediction errors of DGs introduce strong uncertainty to the dynamic economic dispatch problem in active distribution networks.

然而，现有的确定性配电网动态随机经济调度方法并未考虑上述不确定性的存在，在调度问题建模过程中仅采用分布式电源功率和负荷的预测值。另一方面，传统基于机会约束的配电网动态随机经济调度方法在实际应用中面临两大问题：(1)需要精确的随机变量概率密度函数，而该函数在现实中大多数难以获得；(2)该方法所建立的随机优化模型基本上基于抽样场景法，计算量过大。However, the existing deterministic distribution network dynamic stochastic economic dispatch methods do not consider the existence of the above-mentioned uncertainties, and only use the predicted values of distributed power and load in the process of dispatching problem modeling. On the other hand, the traditional dynamic stochastic economic dispatch methods for distribution networks based on chance constraints face two major problems in practical applications: (1) an accurate random variable probability density function is required, and most of this function is difficult to obtain in reality; ( 2) The stochastic optimization model established by this method is basically based on the sampling scene method, and the calculation amount is too large.

发明内容SUMMARY OF THE INVENTION

本发明的目的是为克服已有技术的不足之处，提出一种主动配电网动态随机经济调度方法。本发明在求解主动配电网动态随机经济调度问题时考虑功率随机性，可以用于解决功率预测不确定性的主动配电网经济调度问题，调度策略具备鲁棒性和可靠性，有较高的应用价值。The purpose of the present invention is to propose a dynamic random economic dispatch method for active distribution network in order to overcome the deficiencies of the prior art. The present invention considers the power randomness when solving the dynamic random economic dispatching problem of the active distribution network, and can be used to solve the economical dispatching problem of the active distribution network with the uncertainty of power prediction. The dispatching strategy has robustness and reliability, and has higher application value.

本发明提出一种主动配电网动态随机经济调度方法，其特征在于，该方法包括以下步骤：The present invention proposes a dynamic random economic dispatch method for an active distribution network, which is characterized in that the method comprises the following steps:

1)建立主动配电网动态随机经济调度模型，该模型由目标函数和约束条件构成；具体步骤如下：1) Establish a dynamic stochastic economic dispatch model of active distribution network, which consists of objective functions and constraints; the specific steps are as follows:

1-1)建立模型的目标函数；表达式如式(1)所示：1-1) Establish the objective function of the model; the expression is shown in formula (1):

Figure GDA0002562561940000021
Figure GDA0002562561940000021

其中，

Figure GDA0002562561940000022
为机组j在t时段的发电功率；ΨG为主动配电网中所有发电机组集合；Ai,2、Ai,1、Ai,0分别为发电机组j的耗量特性系数；Γ为调度周期总的时段数；in,
Figure GDA0002562561940000022
is the generating power of unit j at time t; Ψ G is the set of all generating units in the active distribution network; A i,2 , A i,1 , and A i,0 are the consumption characteristic coefficients of generating unit j respectively; Γ is The total number of periods in the scheduling period;

1-2)确定模型的约束条件；具体如下：1-2) Determine the constraints of the model; the details are as follows:

1-2-1)主动配电网的功率平衡约束，如式(2)所示：1-2-1) Power balance constraints of active distribution network, as shown in formula (2):

Figure GDA0002562561940000023
Figure GDA0002562561940000023

其中，

Figure GDA0002562561940000024
为节点i在t时段有功负荷实际功率，
Figure GDA0002562561940000025
为节点i在t时段分布式电源实际功率，PLoss,t为主动配电网在t时段的网损值，Ψn为主动配电网中所有节点的集合；in,
Figure GDA0002562561940000024
is the actual power of the active load of node i in period t,
Figure GDA0002562561940000025
is the actual power of the distributed power generation of node i in the period t, P Loss,t is the network loss value of the active distribution network in the period t, and Ψ n is the set of all nodes in the active distribution network;

1-2-2)主动配电网的网损约束，如式(3)所示：1-2-2) Network loss constraints of active distribution network, as shown in formula (3):

Figure GDA0002562561940000026
Figure GDA0002562561940000026

其中，wLoss为发电网损系数；Among them, w Loss is the power generation network loss coefficient;

1-2-3)发电机机组功率约束，如式(4)所示：1-2-3) The power constraint of the generator set, as shown in formula (4):

Figure GDA0002562561940000027
Figure GDA0002562561940000027

其中，

Figure GDA0002562561940000028
分别为发电机组j的功率下限和上限；in,
Figure GDA0002562561940000028
are the lower and upper power limits of generator set j, respectively;

1-2-4)发电机组的爬坡速率约束，如式(5)所示：1-2-4) The ramp rate constraint of the generator set, as shown in formula (5):

Figure GDA0002562561940000029
Figure GDA0002562561940000029

其中，

Figure GDA00025625619400000210
分别为发电机组j的下爬坡限值和上爬坡限值；in,
Figure GDA00025625619400000210
are the down-climbing limit and the up-climbing limit of generator set j, respectively;

1-2-5)主动配电网的上下旋转备用约束，分别如式(6)和式(7)所示：1-2-5) The upper and lower rotation reserve constraints of the active distribution network are shown in equations (6) and (7) respectively:

Figure GDA00025625619400000211
Figure GDA00025625619400000211

Figure GDA00025625619400000212
Figure GDA00025625619400000212

其中，

Figure GDA00025625619400000213
分别为主动配电网在t时段的上、下备用要求；in,
Figure GDA00025625619400000213
are the upper and lower standby requirements of the active distribution network in the period t, respectively;

2)对步骤1)的约束条件进行转化；具体步骤如下：2) Transform the constraints of step 1); the specific steps are as follows:

2-1)根据约束条件式(6)和(7)构建机会约束，分别如式(8)和(9)所示：2-1) Construct the chance constraint according to the constraint equations (6) and (7), as shown in equations (8) and (9) respectively:

Figure GDA00025625619400000214
Figure GDA00025625619400000214

Figure GDA0002562561940000031
Figure GDA0002562561940000031

其中，Pr()为事件发生的概率，ξ为给定的不等式约束式(8)和(9)被破坏的概率；Among them, Pr() is the probability of event occurrence, ξ is the probability that the given inequality constraints (8) and (9) are destroyed;

2-2)收集配电网中所有节点在每个时段的有功负荷的预测误差值集合记为

Figure GDA0002562561940000032
收集配电网中所有节点在每个时段的分布式电源功率预测误差值集合记为
Figure GDA0002562561940000033
其中
Figure GDA0002562561940000034
为节点i在t时段的有功负荷预测误差，
Figure GDA0002562561940000035
为节点i在t时段的分布式电源功率预测误差；2-2) Collect the prediction error value set of the active load of all nodes in the distribution network in each period, and denote it as
Figure GDA0002562561940000032
The set of distributed power prediction error values of all nodes in the distribution network at each time period is collected and denoted as
Figure GDA0002562561940000033
in
Figure GDA0002562561940000034
is the active load prediction error of node i at time t,
Figure GDA0002562561940000035
is the prediction error of distributed power supply power of node i in time period t;

分别对

Figure GDA0002562561940000036
求取对应的误差标幺参数，如式(10)和(11)所示：respectively
Figure GDA0002562561940000036
Find the corresponding error per-unit parameters, as shown in equations (10) and (11):

Figure GDA0002562561940000037
Figure GDA0002562561940000037

Figure GDA0002562561940000038
Figure GDA0002562561940000038

其中，max(||)为求集合中元素绝对值的最大值；

Figure GDA0002562561940000039
为节点i在t时段有功负荷误差标幺参数，
Figure GDA00025625619400000310
为节点i在t时段分布式电源功率误差标幺参数；Among them, max(||) is the maximum value of the absolute value of the elements in the set;
Figure GDA0002562561940000039
is the per-unit parameter of active load error of node i at time t,
Figure GDA00025625619400000310
is the per-unit parameter of the power error of the distributed power supply for node i in the t period;

设定

Figure GDA00025625619400000311
为节点i在t时段的有功负荷标幺化预测误差，
Figure GDA00025625619400000312
为节点i在t时段的分布式电源功率标幺化预测误差，
Figure GDA00025625619400000313
的概率分布集合分别为
Figure GDA00025625619400000314
为定义在[-1,1]上且均值为0的任意相互独立分布组成的集合；set up
Figure GDA00025625619400000311
is the per-unit prediction error of the active load of node i in the period t,
Figure GDA00025625619400000312
is the per-unit prediction error of the distributed power supply power of node i in time period t,
Figure GDA00025625619400000313
The set of probability distributions of
Figure GDA00025625619400000314
is a set of any mutually independent distributions defined on [-1,1] with mean 0;

2-3)将配电网中有功负荷实际功率与分布式电源实际功率分别表示为式(12)和(13)所示的形式：2-3) The actual power of the active load and the actual power of the distributed power supply in the distribution network are expressed as formulas (12) and (13), respectively:

Figure GDA00025625619400000315
Figure GDA00025625619400000315

Figure GDA00025625619400000316
Figure GDA00025625619400000316

Figure GDA00025625619400000317
Figure GDA00025625619400000317

其中，

Figure GDA00025625619400000318
为节点i在t时段有功负荷预测功率，
Figure GDA00025625619400000323
为节点i在t时段分布式电源预测功率；in,
Figure GDA00025625619400000318
is the predicted power for the active load of node i at time t,
Figure GDA00025625619400000323
Predict the power of distributed power generation for node i at time period t;

2-4)根据式(2)、(3)、(10)、(11)、(12)、(13)，将

Figure GDA00025625619400000319
表示为如式(14)所示：2-4) According to formula (2), (3), (10), (11), (12), (13), the
Figure GDA00025625619400000319
It is expressed as formula (14):

Figure GDA00025625619400000320
Figure GDA00025625619400000320

2-5)将式(14)代入约束条件式(8)和(9)中，根据机会约束凸松弛转化方法，则式(8)和(9)分别转化为式(15)和(16)所示：2-5) Substitute equation (14) into constraint equations (8) and (9), according to the opportunity constraint convex relaxation transformation method, equations (8) and (9) are transformed into equations (15) and (16) respectively shown:

Figure GDA00025625619400000321
Figure GDA00025625619400000321

Figure GDA00025625619400000322
Figure GDA00025625619400000322

其中，系数

Figure GDA0002562561940000041
定义如式(17)所示：Among them, the coefficient
Figure GDA0002562561940000041
The definition is shown in formula (17):

Figure GDA0002562561940000042
Figure GDA0002562561940000042

3)对模型求解；3) Solve the model;

根据目标函数式(1)，约束条件式(2)、(3)、(4)、(5)、(15)、(16)、(17)，应用凸规划算法对模型进行求解；最终得到

Figure GDA0002562561940000043
即是各发电机组t时段最优经济调度的发电功率，各时段的
Figure GDA0002562561940000044
汇集起来为各发电机组全天各时段最优经济调度的发电功率。According to the objective function formula (1), the constraint condition formulae (2), (3), (4), (5), (15), (16), (17), the convex programming algorithm is used to solve the model;
Figure GDA0002562561940000043
That is, the optimal economic dispatching power of each generating set in period t, and the power in each period is
Figure GDA0002562561940000044
Collected together to generate power for the optimal economic dispatch of each generator set throughout the day and time period.

本发明的特点及有益效果在于：The characteristics and beneficial effects of the present invention are:

本发明提出一种主动配电网动态随机经济调度算法，根据已知的统计信息构建出一个不确定量的概率分布集，构建包含旋转备用约束的机会约束，并利用凸松弛将其转化为确定性线性约束，使该主动配电网动态随机经济调度问题得到有效求解。本发明符合生产的实际情况，所得到的调度策略具备鲁棒性和可靠性，在考虑运行随机性的情况下计算效率高，有较高的应用价值。The invention proposes a dynamic stochastic economic dispatch algorithm for active distribution network, constructs a probability distribution set of uncertain quantities according to known statistical information, constructs chance constraints including rotating reserve constraints, and uses convex relaxation to convert them into deterministic Therefore, the dynamic stochastic economic dispatch problem of active distribution network can be effectively solved. The invention conforms to the actual situation of production, and the obtained scheduling strategy has robustness and reliability, and has high calculation efficiency and high application value when the randomness of operation is considered.

具体实施方式Detailed ways

本发明提出的一种主动配电网动态随机经济调度方法，下面结合具体实施例进一步详细说明如下。A dynamic random economic dispatch method for an active distribution network proposed by the present invention is further described in detail below with reference to specific embodiments.

本发明提出的一种主动配电网动态随机经济调度方法，包括以下步骤：A dynamic random economic dispatch method for an active distribution network proposed by the present invention includes the following steps:

1)建立主动配电网动态随机经济调度模型，该模型由目标函数和约束条件构成；具体步骤如下：1) Establish a dynamic stochastic economic dispatch model of active distribution network, which consists of objective functions and constraints; the specific steps are as follows:

1-1)建立模型的目标函数；表达式如式(1)所示：1-1) Establish the objective function of the model; the expression is shown in formula (1):

Figure GDA0002562561940000045
Figure GDA0002562561940000045

其中，

Figure GDA0002562561940000046
为机组j在t时段的发电功率；ΨG为主动配电网中所有发电机组集合；Ai,2、Ai,1、Ai,0分别为发电机组j(可为燃煤、燃气、水电等机组)的耗量特性系数(上述系数由配电网调度中心给出)；Γ为调度周期总的时段数，通常取为96；in,
Figure GDA0002562561940000046
is the generating power of unit j in the period t; Ψ G is the set of all generating units in the active distribution network; A i,2 , A i,1 , A i,0 are the generating units j (can be coal, gas, (the above coefficients are given by the distribution network dispatch center); Γ is the total number of time periods in the dispatch cycle, usually taken as 96;

1-2)确定模型的约束条件；具体如下：1-2) Determine the constraints of the model; the details are as follows:

1-2-1)主动配电网的功率平衡约束，如式(2)所示：1-2-1) Power balance constraints of active distribution network, as shown in formula (2):

Figure GDA0002562561940000051
Figure GDA0002562561940000051

其中，

Figure GDA0002562561940000052
为节点i在t时段有功负荷实际功率，
Figure GDA0002562561940000053
为节点i在t时段分布式电源实际功率，PLoss,t为主动配电网在t时段的网损值，Ψn为主动配电网中所有节点的集合；in,
Figure GDA0002562561940000052
is the actual power of the active load of node i in period t,
Figure GDA0002562561940000053
is the actual power of the distributed power generation of node i in the period t, P Loss,t is the network loss value of the active distribution network in the period t, and Ψ n is the set of all nodes in the active distribution network;

1-2-2)主动配电网的网损约束，如式(3)所示：1-2-2) Network loss constraints of active distribution network, as shown in formula (3):

Figure GDA0002562561940000054
Figure GDA0002562561940000054

其中，wLoss为发电网损系数(常规取值为0.01-0.05)；Among them, w Loss is the power generation network loss coefficient (the conventional value is 0.01-0.05);

1-2-3)发电机机组功率约束，如式(4)所示：1-2-3) The power constraint of the generator set, as shown in formula (4):

Figure GDA0002562561940000055
Figure GDA0002562561940000055

其中，

Figure GDA00025625619400000517
分别为发电机组j的功率下限和上限；in,
Figure GDA00025625619400000517
are the lower and upper power limits of generator set j, respectively;

1-2-4)发电机组的爬坡速率约束，如式(5)所示：1-2-4) The ramp rate constraint of the generator set, as shown in formula (5):

Figure GDA0002562561940000056
Figure GDA0002562561940000056

其中，

Figure GDA0002562561940000057
分别为发电机组j的下爬坡限值和上爬坡限值；in,
Figure GDA0002562561940000057
are the down-climbing limit and the up-climbing limit of generator set j, respectively;

1-2-5)主动配电网的上下旋转备用约束，分别如式(6)和式(7)所示：1-2-5) The upper and lower rotation reserve constraints of the active distribution network are shown in equations (6) and (7) respectively:

Figure GDA0002562561940000058
Figure GDA0002562561940000058

Figure GDA0002562561940000059
Figure GDA0002562561940000059

其中，

Figure GDA00025625619400000510
分别为主动配电网在t时段的上、下备用要求，一般取为系统总负荷的5％；in,
Figure GDA00025625619400000510
are respectively the upper and lower standby requirements of the active distribution network in the t period, generally taken as 5% of the total system load;

2)对步骤1)的约束条件进行转化；具体步骤如下：2) Transform the constraints of step 1); the specific steps are as follows:

2-1)根据约束条件式(6)和(7)构建机会约束，分别如式(8)和(9)所示：2-1) Construct the chance constraint according to the constraint equations (6) and (7), as shown in equations (8) and (9) respectively:

Figure GDA00025625619400000511
Figure GDA00025625619400000511

Figure GDA00025625619400000512
Figure GDA00025625619400000512

其中，Pr()为事件发生的概率，ξ为给定的不等式约束式(8)和(9)被破坏的概率，取值范围为[0,1]，本实例中取值为0.1；Among them, Pr() is the probability of event occurrence, ξ is the probability that the given inequality constraints (8) and (9) are destroyed, the value range is [0, 1], and in this example, the value is 0.1;

2-2)收集配电网中所有节点在每个时段的有功负荷的预测误差值集合记为

Figure GDA00025625619400000513
收集配电网中所有节点在每个时段的分布式电源功率预测误差值集合记为
Figure GDA00025625619400000514
其中
Figure GDA00025625619400000515
为节点i在t时段的有功负荷预测误差，
Figure GDA00025625619400000516
为节点i在t时段的分布式电源功率预测误差；收集全天各个时段(每隔15分钟一个时段，全天共96各时段)的误差数据，各时段数据量依据预测机构数据提供程度，越多越好。预测误差具体指为测量实际值与预测值之差(实际功率减去对应的预测功率)；2-2) Collect the prediction error value set of the active load of all nodes in the distribution network in each period, and denote it as
Figure GDA00025625619400000513
The set of distributed power prediction error values of all nodes in the distribution network at each time period is collected and denoted as
Figure GDA00025625619400000514
in
Figure GDA00025625619400000515
is the active load prediction error of node i at time t,
Figure GDA00025625619400000516
is the power prediction error of the distributed power supply of node i in period t; the error data of each period of the day (one period every 15 minutes, a total of 96 periods of the whole day) are collected. More is better. The prediction error specifically refers to the difference between the measured actual value and the predicted value (the actual power minus the corresponding predicted power);

分别对

Figure GDA0002562561940000061
求取对应的误差标幺参数，如式(10)和(11)所示：respectively
Figure GDA0002562561940000061
Find the corresponding error per-unit parameters, as shown in equations (10) and (11):

Figure GDA0002562561940000062
Figure GDA0002562561940000062

Figure GDA0002562561940000063
Figure GDA0002562561940000063

其中，max(||)为求集合中元素绝对值的最大值；

Figure GDA0002562561940000064
为节点i在t时段有功负荷误差标幺参数，
Figure GDA0002562561940000065
为节点i在t时段分布式电源功率误差标幺参数；Among them, max(||) is the maximum value of the absolute value of the elements in the set;
Figure GDA0002562561940000064
is the per-unit parameter of active load error of node i at time t,
Figure GDA0002562561940000065
is the per-unit parameter of the power error of the distributed power supply for node i in the t period;

设定

Figure GDA0002562561940000066
为节点i在t时段的有功负荷标幺化预测误差，
Figure GDA0002562561940000067
为节点i在t时段的分布式电源功率标幺化预测误差，
Figure GDA0002562561940000068
的概率分布集合分别为
Figure GDA0002562561940000069
为定义在[-1,1]上且均值为0的任意相互独立分布组成的集合；set up
Figure GDA0002562561940000066
is the per-unit prediction error of the active load of node i in the period t,
Figure GDA0002562561940000067
is the per-unit prediction error of the distributed power supply power of node i in time period t,
Figure GDA0002562561940000068
The set of probability distributions of
Figure GDA0002562561940000069
is a set of any mutually independent distributions defined on [-1,1] with mean 0;

2-3)将配电网中有功负荷实际功率与分布式电源实际功率分别表示为式(12)和(13)所示的形式：2-3) The actual power of the active load and the actual power of the distributed power supply in the distribution network are expressed as formulas (12) and (13), respectively:

Figure GDA00025625619400000610
Figure GDA00025625619400000610

Figure GDA00025625619400000611
Figure GDA00025625619400000611

Figure GDA00025625619400000612
Figure GDA00025625619400000612

其中，

Figure GDA00025625619400000613
为节点i在t时段有功负荷预测功率，
Figure GDA00025625619400000614
为节点i在t时段分布式电源预测功率，上述两种预测功率由专门的预测机构给出；
Figure GDA00025625619400000615
为节点i在t时段的有功负荷不确定性变量(标幺化预测误差)，
Figure GDA00025625619400000616
为节点i在t时段有功负荷误差标幺参数(大于等于0)，
Figure GDA00025625619400000617
为节点i在t时段的分布式电源功率不确定变量(标幺化预测误差)，
Figure GDA00025625619400000618
为节点i在t时段分布式电源功率误差标幺参数(大于等于0)。本发明内容中参数全部认为已经过标幺化，单位皆为1。in,
Figure GDA00025625619400000613
is the predicted power for the active load of node i at time t,
Figure GDA00025625619400000614
Predict the power of the distributed power supply for node i in the period t, the above two predicted powers are given by a special prediction agency;
Figure GDA00025625619400000615
is the active load uncertainty variable of node i at time t (per unit forecast error),
Figure GDA00025625619400000616
is the per-unit parameter (greater than or equal to 0) of the active load error of node i in the t period,
Figure GDA00025625619400000617
is the uncertain variable of distributed power supply power of node i in time period t (per unit forecast error),
Figure GDA00025625619400000618
is the per-unit parameter (greater than or equal to 0) of the power error of the distributed power supply for node i in the t period. All parameters in the content of the present invention are considered to have been per-unitized, and the unit is all 1.

2-4)根据式(2)、(3)、(10)、(11)、(12)、(13)，将

Figure GDA00025625619400000619
表示为如式(14)所示：2-4) According to formula (2), (3), (10), (11), (12), (13), the
Figure GDA00025625619400000619
It is expressed as formula (14):

Figure GDA00025625619400000620
Figure GDA00025625619400000620

2-5)将式(14)代入约束条件式(8)和(9)中，根据机会约束凸松弛转化方法，则式(8)和(9)分别转化为式(15)和(16)所示：2-5) Substitute equation (14) into constraint equations (8) and (9), according to the opportunity constraint convex relaxation transformation method, equations (8) and (9) are transformed into equations (15) and (16) respectively shown:

Figure GDA00025625619400000621
Figure GDA00025625619400000621

Figure GDA0002562561940000071
Figure GDA0002562561940000071

其中，系数

Figure GDA0002562561940000072
定义如式(17)所示：Among them, the coefficient
Figure GDA0002562561940000072
The definition is shown in formula (17):

Figure GDA0002562561940000073
Figure GDA0002562561940000073

3)对模型求解；3) Solve the model;

根据目标函数式(1)，约束条件式(2)、(3)、(4)、(5)、(15)、(16)、(17)，应用凸规划算法对模型进行求解；最终得到

Figure GDA0002562561940000074
即是各发电机组t时段最优经济调度的发电功率。各时段的
Figure GDA0002562561940000075
汇集起来为各发电机组全天各时段经济调度的发电功率。According to the objective function formula (1), the constraint condition formulae (2), (3), (4), (5), (15), (16), (17), the convex programming algorithm is used to solve the model;
Figure GDA0002562561940000074
That is, the optimal economic dispatching power of each generating set in period t. of each period
Figure GDA0002562561940000075
Collected together to generate power for economic dispatch of each generator set at all times of the day.

Claims (1)

1.一种主动配电网动态随机经济调度方法，其特征在于，该方法包括以下步骤：1. A dynamic random economic dispatch method for active distribution network, characterized in that, the method comprises the following steps: 1)建立主动配电网动态随机经济调度模型，该模型由目标函数和约束条件构成；具体步骤如下：1) Establish a dynamic stochastic economic dispatch model of active distribution network, which consists of objective functions and constraints; the specific steps are as follows: 1-1)建立模型的目标函数；表达式如式(1)所示：1-1) Establish the objective function of the model; the expression is shown in formula (1):
Figure FDA0002562561930000011
Figure FDA0002562561930000011
 其中，
Figure FDA0002562561930000012
为机组j在t时段的发电功率；ΨG为主动配电网中所有发电机组集合；Ai,2、Ai,1、Ai,0分别为发电机组j的耗量特性系数；Γ为调度周期总的时段数；in,
Figure FDA0002562561930000012
is the generating power of unit j at time t; Ψ G is the set of all generating units in the active distribution network; A i,2 , A i,1 , and A i,0 are the consumption characteristic coefficients of generating unit j respectively; Γ is The total number of periods in the scheduling period; 1-2)确定模型的约束条件；具体如下：1-2) Determine the constraints of the model; the details are as follows: 1-2-1)主动配电网的功率平衡约束，如式(2)所示：1-2-1) Power balance constraints of active distribution network, as shown in formula (2):
Figure FDA0002562561930000013
Figure FDA0002562561930000013
 其中，
Figure FDA0002562561930000014
为节点i在t时段有功负荷实际功率，
Figure FDA0002562561930000015
为节点i在t时段分布式电源实际功率，PLoss,t为主动配电网在t时段的网损值，Ψn为主动配电网中所有节点的集合；in,
Figure FDA0002562561930000014
is the actual power of the active load of node i in period t,
Figure FDA0002562561930000015
is the actual power of the distributed power generation of node i in the period t, P Loss,t is the network loss value of the active distribution network in the period t, and Ψ n is the set of all nodes in the active distribution network; 1-2-2)主动配电网的网损约束，如式(3)所示：1-2-2) Network loss constraints of active distribution network, as shown in formula (3):
Figure FDA0002562561930000016
Figure FDA0002562561930000016
 其中，wLoss为发电网损系数；Among them, w Loss is the power generation network loss coefficient; 1-2-3)发电机机组功率约束，如式(4)所示：1-2-3) The power constraint of the generator set, as shown in formula (4):
Figure FDA0002562561930000017
Figure FDA0002562561930000017
 其中，
Figure FDA0002562561930000018
分别为发电机组j的功率下限和上限；in,
Figure FDA0002562561930000018
are the lower and upper power limits of generator set j, respectively; 1-2-4)发电机组的爬坡速率约束，如式(5)所示：1-2-4) The ramp rate constraint of the generator set, as shown in formula (5):
Figure FDA0002562561930000019
Figure FDA0002562561930000019
 其中，
Figure FDA00025625619300000110
分别为发电机组j的下爬坡限值和上爬坡限值；in,
Figure FDA00025625619300000110
are the down-climbing limit and the up-climbing limit of generator set j, respectively; 1-2-5)主动配电网的上下旋转备用约束，分别如式(6)和式(7)所示：1-2-5) The upper and lower rotation reserve constraints of the active distribution network are shown in equations (6) and (7) respectively:
Figure FDA00025625619300000111
Figure FDA00025625619300000111
Figure FDA00025625619300000112
Figure FDA00025625619300000112
 其中，
Figure FDA00025625619300000113
分别为主动配电网在t时段的上、下备用要求；in,
Figure FDA00025625619300000113
are the upper and lower standby requirements of the active distribution network in the period t, respectively; 2)对步骤1)的约束条件进行转化；具体步骤如下：2) Transform the constraints of step 1); the specific steps are as follows: 2-1)根据约束条件式(6)和(7)构建机会约束，分别如式(8)和(9)所示：2-1) Construct the chance constraint according to the constraint equations (6) and (7), as shown in equations (8) and (9) respectively:
Figure FDA0002562561930000021
Figure FDA0002562561930000021
Figure FDA0002562561930000022
Figure FDA0002562561930000022
 其中，Pr()为事件发生的概率，ξ为不等式约束式(8)和(9)被破坏的概率；Among them, Pr() is the probability of event occurrence, ξ is the probability that inequality constraints (8) and (9) are broken; 2-2)收集配电网中所有节点在每个时段的有功负荷的预测误差值集合记为
Figure FDA0002562561930000023
收集配电网中所有节点在每个时段的分布式电源功率预测误差值集合记为
Figure FDA0002562561930000024
其中
Figure FDA0002562561930000025
为节点i在t时段的有功负荷预测误差，
Figure FDA0002562561930000026
为节点i在t时段的分布式电源功率预测误差；2-2) Collect the prediction error value set of the active load of all nodes in the distribution network in each period, and denote it as
Figure FDA0002562561930000023
The set of distributed power prediction error values of all nodes in the distribution network at each time period is collected and denoted as
Figure FDA0002562561930000024
in
Figure FDA0002562561930000025
is the active load prediction error of node i at time t,
Figure FDA0002562561930000026
is the prediction error of distributed power supply power of node i in time period t; 分别对
Figure FDA0002562561930000027
求取对应的误差标幺参数，如式(10)和(11)所示：respectively
Figure FDA0002562561930000027
Find the corresponding error per-unit parameters, as shown in equations (10) and (11):
Figure FDA0002562561930000028
Figure FDA0002562561930000028
Figure FDA0002562561930000029
Figure FDA0002562561930000029
 其中，max(||)为求集合中元素绝对值的最大值；
Figure FDA00025625619300000210
为节点i在t时段有功负荷误差标幺参数，
Figure FDA00025625619300000211
为节点i在t时段分布式电源功率误差标幺参数；Among them, max(||) is the maximum value of the absolute value of the elements in the set;
Figure FDA00025625619300000210
is the per-unit parameter of active load error of node i at time t,
Figure FDA00025625619300000211
is the per-unit parameter of the power error of the distributed power supply for node i in the t period; 设定
Figure FDA00025625619300000212
为节点i在t时段的有功负荷标幺化预测误差，
Figure FDA00025625619300000213
为节点i在t时段的分布式电源功率标幺化预测误差，
Figure FDA00025625619300000214
的概率分布集合分别为
Figure FDA00025625619300000215
为定义在[-1,1]上且均值为0的任意相互独立分布组成的集合；set up
Figure FDA00025625619300000212
is the per-unit prediction error of the active load of node i in the period t,
Figure FDA00025625619300000213
is the per-unit prediction error of the distributed power supply power of node i in time period t,
Figure FDA00025625619300000214
The set of probability distributions are
Figure FDA00025625619300000215
is a set of any mutually independent distributions defined on [-1,1] with mean 0; 2-3)将配电网中有功负荷实际功率与分布式电源实际功率分别表示为式(12)和(13)所示的形式：2-3) The actual power of the active load and the actual power of the distributed power supply in the distribution network are expressed as formulas (12) and (13), respectively:
Figure FDA00025625619300000216
Figure FDA00025625619300000216
Figure FDA00025625619300000217
Figure FDA00025625619300000217
Figure FDA00025625619300000218
Figure FDA00025625619300000218
 其中，
Figure FDA00025625619300000219
为节点i在t时段有功负荷预测功率，
Figure FDA00025625619300000220
为节点i在t时段分布式电源预测功率；in,
Figure FDA00025625619300000219
is the predicted power for the active load of node i at time t,
Figure FDA00025625619300000220
Predict the power of distributed power generation for node i at time period t; 2-4)根据式(2)、(3)、(10)、(11)、(12)、(13)，将
Figure FDA00025625619300000221
表示为如式(14)所示：2-4) According to formula (2), (3), (10), (11), (12), (13), the
Figure FDA00025625619300000221
It is expressed as formula (14):
Figure FDA00025625619300000222
Figure FDA00025625619300000222
 2-5)将式(14)代入约束条件式(8)和(9)中，根据机会约束凸松弛转化方法，则式(8)和(9)分别转化为式(15)和(16)所示：2-5) Substitute equation (14) into constraint equations (8) and (9), according to the opportunity constraint convex relaxation transformation method, equations (8) and (9) are transformed into equations (15) and (16) respectively shown:
Figure FDA0002562561930000031
Figure FDA0002562561930000031
Figure FDA0002562561930000032
Figure FDA0002562561930000032
 其中，系数
Figure FDA0002562561930000033
定义如式(17)所示：Among them, the coefficient
Figure FDA0002562561930000033
The definition is shown in formula (17):
Figure FDA0002562561930000034
Figure FDA0002562561930000034
 3)对模型求解；3) Solve the model; 根据目标函数式(1)，约束条件式(2)、(3)、(4)、(5)、(15)、(16)、(17)，应用凸规划算法对模型进行求解；最终得到
Figure FDA0002562561930000035
即是各发电机组t时段最优经济调度的发电功率，各时段的
Figure FDA0002562561930000036
汇集起来为各发电机组全天各时段最优经济调度的发电功率。According to the objective function formula (1), the constraint condition formulae (2), (3), (4), (5), (15), (16), (17), the convex programming algorithm is used to solve the model;
Figure FDA0002562561930000035
That is, the optimal economic dispatching power of each generating set in period t, and the power in each period is
Figure FDA0002562561930000036
Collected together to generate power for the optimal economic dispatch of each generator set throughout the day and time period.
CN201711205188.3A 2017-11-27 2017-11-27 A dynamic stochastic economic dispatch method for active distribution network Active CN108039739B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711205188.3A CN108039739B (en) 2017-11-27 2017-11-27 A dynamic stochastic economic dispatch method for active distribution network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711205188.3A CN108039739B (en) 2017-11-27 2017-11-27 A dynamic stochastic economic dispatch method for active distribution network

Publications (2)

Publication Number Publication Date
CN108039739A CN108039739A (en) 2018-05-15
CN108039739B true CN108039739B (en) 2020-09-25

Family

ID=62093453

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711205188.3A Active CN108039739B (en) 2017-11-27 2017-11-27 A dynamic stochastic economic dispatch method for active distribution network

Country Status (1)

Country Link
CN (1) CN108039739B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109713713B (en) * 2018-10-19 2021-03-02 云南电网有限责任公司 A stochastic optimization method for unit start-stop based on chance-constrained convex relaxation
CN109657898B (en) * 2018-10-19 2022-11-01 云南电网有限责任公司 Renewable energy random dynamic economic dispatching method based on convex relaxation

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104600747A (en) * 2015-01-21 2015-05-06 西安交通大学 Operation optimizing method capable of coordinating operation risk and wind energy consumption of power system
CN105207272A (en) * 2015-09-18 2015-12-30 武汉大学 Electric power system dynamic random economic dispatching method and device based on general distribution
CN105244869A (en) * 2015-10-13 2016-01-13 国网山东省电力公司电力科学研究院 Dynamic random scheduling control method for power distribution network containing micro-grid
CN105846425A (en) * 2016-04-08 2016-08-10 江苏省电力试验研究院有限公司 Economic dispatching method based on general wind power forecasting error model
CN106099984A (en) * 2016-07-29 2016-11-09 清华大学 A kind of active distribution network distributed power source heap(ed) capacity appraisal procedure of data-driven

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104600747A (en) * 2015-01-21 2015-05-06 西安交通大学 Operation optimizing method capable of coordinating operation risk and wind energy consumption of power system
CN105207272A (en) * 2015-09-18 2015-12-30 武汉大学 Electric power system dynamic random economic dispatching method and device based on general distribution
CN105244869A (en) * 2015-10-13 2016-01-13 国网山东省电力公司电力科学研究院 Dynamic random scheduling control method for power distribution network containing micro-grid
CN105846425A (en) * 2016-04-08 2016-08-10 江苏省电力试验研究院有限公司 Economic dispatching method based on general wind power forecasting error model
CN106099984A (en) * 2016-07-29 2016-11-09 清华大学 A kind of active distribution network distributed power source heap(ed) capacity appraisal procedure of data-driven

Also Published As

Publication number Publication date
CN108039739A (en) 2018-05-15

Similar Documents

Publication Publication Date Title
CN112186766B (en) Multistage decomposition method based on Benders decomposition
CN103296679B (en) The medium-term and long-term long-term wind power run that optimizes of power system is exerted oneself model modelling approach
CN113705962B (en) Virtual power plant day-ahead scheduling method based on distribution robust optimization
CN110400056B (en) Method and device for cascade hydropower day-ahead optimized scheduling based on wind-light space-time characteristics
CN103997039B (en) Method for predicting rotating standby interval with wind power acceptance considered based on probability interval prediction
CN107910863A (en) Consider the power distribution network dispatching method that photovoltaic is contributed with workload demand forecast interval
CN107292766A (en) Towards the power system peak regulation means economic evaluation method and system of wind electricity digestion
CN107977744A (en) A kind of electric system based on traditional Benders decomposition methods Robust Scheduling method a few days ago
CN107240933A (en) A kind of wind fire for considering wind power characteristic coordinates rolling scheduling method
CN104299173B (en) It is a kind of to optimize dispatching method a few days ago suitable for the robust that various energy resources are accessed
CN104682447A (en) Power system economic dispatching method containing multiple wind power plants
CN109299862B (en) A kind of convex loose appraisal procedure of wind-powered electricity generation maximum digestion capability
CN102904248B (en) Electric power system dispatching method based on wind electricity output uncertainty aggregation
CN108039739B (en) A dynamic stochastic economic dispatch method for active distribution network
CN112186734A (en) Medium-and-long-term operation simulation method for power system, storage medium and computing equipment
CN116111597A (en) A method, system, memory and device for constructing a power flow section scene set in a medium- and long-term dispatch plan
CN111327080A (en) Multi-category unit oriented medium-and-long-term electric quantity combined decomposition method for power system
CN115187018A (en) A double-layer optimal scheduling method and device for a comprehensive energy system in a park
CN106160009A (en) A kind of interconnected electric power system based on risk node electricity price dispersion coordinated dispatching method
CN107968397B (en) Distribution network transmission capacity calculation method considering operation randomness
CN109255102B (en) Photovoltaic power generation maximum grid-connected capacity evaluation method based on convex relaxation
CN104659818A (en) Optimal allocation method for positive and negative spinning reserve capacity in system comprising wind farm
CN104951650A (en) Method for evaluating outer power transmission trading capacity of power exchange point of large-scale wind power grid
Zhang et al. Two‐Stage Optimization Model of Centralized Energy Storage Participating in Peak Shaving with Maximum Reserve Capacity and Minimum Carbon Emission of the System
CN114692988B (en) CatBoost-based net load spare capacity demand prediction method, device and storage medium

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
GR01 Patent grant
GR01 Patent grant