CN111952968B - Configuration method, device and readable storage medium for incremental distributed power supply equipment - Google Patents

Abstract

The present invention provides a configuration method, device and readable storage medium for incremental distributed power equipment, wherein the method includes the following steps: constructing a distributed power incremental distribution network model including at least a node distribution network structure; In the distribution network model, the corresponding objective functions and function constraints of distributed power nodes, incremental distribution network nodes and power nodes are multiplied by the corresponding weights to obtain the corresponding weight functions. The weight functions corresponding to the distribution network nodes and the power nodes are added to obtain the overall planning function; the function constraints of the incremental distribution network nodes are subjected to second-order cone relaxation processing, and based on the processed function constraints and the unprocessed function constraints , the overall planning function is solved by the solver, and the power capacity of the node to be connected to the distributed power equipment is obtained. The invention improves the rationality of the power supply capacity of each node in the incremental distribution network system.

Description

Configuration method, apparatus and readable storage medium of incremental distributed power supply equipment

technical field

The present invention relates to the technical field of power systems, and in particular, to a configuration method, apparatus and readable storage medium for incremental distributed power equipment.

Background technique

In principle, incremental distribution network refers to power grids with voltage levels of 110kV (kilovolt) and below and local power grids such as industrial parks (economic development zones) with voltage levels of 220 (330)kV and below, and does not involve the construction of transmission grids of 220kV and above. Incremental distribution networks include new incremental distribution networks, capacity expansion and expansion of distribution networks, and existing distribution networks outside the stock of power grid companies. With the continuous development of urbanization and economic and social development in China, incremental distribution network is indispensable. It is urgent to change the planning and operation mode of traditional distribution network, and to combine incremental distribution network with the latest distribution network operation technology. , which is conducive to achieving the goals of reducing costs, improving quality and efficiency, and saving energy and reducing emissions. In addition, the distributed generation (DG) equipment is integrated into the incremental distribution network system, and the direction of the power flow changes from one-way to two-way, which increases the uncertainty of operation. It also brings certain security and stability challenges. Therefore, under the premise of ensuring the safe and stable operation of the incremental distribution network system, it is urgent to provide an incremental distribution network optimization planning method to improve the rationality of the power supply capacity of each node in the incremental distribution network system.

SUMMARY OF THE INVENTION

Based on the above status quo, it is urgent to provide a method for improving the rationality of the power supply capacity of each node in the incremental distribution network system on the premise of ensuring the safe and stable operation of the incremental distribution network system. The main purpose of the present invention is to provide an incremental The configuration method, device and readable storage medium of distributed power supply equipment are provided to improve the rationality of the power supply capacity of each node in the incremental distribution network system.

For achieving the above object, the technical scheme adopted in the present invention is as follows:

A configuration method for incremental distributed power equipment, the configuration method for incremental distributed power equipment includes the following steps:

S100. Build an incremental distribution network model of distributed power sources in an incremental distribution network system, wherein the incremental distribution network model includes at least a node distribution network structure;

S200: Acquire, from a preset database of the incremental power distribution network system, first target data related to the distributed power supply node and second target data related to the incremental power distribution network node in the incremental power distribution network model data and third target data related to the power node, construct a first target function corresponding to the distributed power supply node and a first function constraint that constrains the first target function according to the first target data, according to the first target function The second objective function corresponding to the incremental distribution network node and the second function constraint constraining the second objective function are constructed with the two objective data, and the third objective corresponding to the power node is constructed according to the third objective data Functions and constraints on the third function constraint of the third objective function, wherein the first functional constraints include the constraints on the number of preset nodes connected to distributed power equipment in the node distribution network structure, and the upper limit constraints on the penetration rate of distributed power equipment and the total power supply capacity constraints of distributed power equipment, the second function constraints include line power flow constraints, voltage constraints and power constraints in the incremental distribution network model, and the third function constraints include the transfer out of the target time. The power load constraints and the power load constraints transferred out at the target time;

S300, multiply the first objective function by the first weight corresponding to the first objective function to obtain a first weight function, and multiply the second objective function by the second weight corresponding to the second objective function, Obtain the second weight function, multiply the third objective function by the third weight corresponding to the third objective function, and obtain the third weight function;

S400, adding the first weight function, the second weight function and the third weight function to obtain an overall planning function;

S500. Perform second-order cone relaxation processing on the second functional constraint to obtain a relaxed second functional constraint, and based on the first functional constraint, the third functional constraint, and the relaxed second functional constraint, pass The solver solves the overall planning function, obtains the power supply capacity of the node to be connected to the distributed power supply device in the node distribution network structure, and configures the distributed power supply device corresponding to the node according to the power supply capacity.

Preferably, in the step S500, during the second-order cone relaxation process, auxiliary variables U _it '=(U _it ) ² and I _ij.t '=(I _ij.t ) ² are defined, and according to the auxiliary variables The variables get the following objective equation:

I _ij.t represents the line current between node i and node j in the node distribution network structure, U _it represents the voltage amplitude of node i at time t, P _ij.t represents the active power transmitted by branch ij at time t, Q _ij.t represents the reactive power transmitted by branch ij at time t;

The objective equation is relaxed to obtain the following relaxed inequality:

Reconstructing the relaxation inequalities with second-order cones yields the following objective inequalities:

Perform second-order cone relaxation processing on each line power flow constraint in the second functional constraint according to the target inequality, to obtain the second functional constraint after relaxation processing;

The expression of each line power flow constraint in the second function constraint after relaxation is:

(U _jt ) ² =U _it ′-2(R _ij P _ij.t +X _ij Q _ij.t )+[(R _ij ) ² +(X _ij ) ² ]·I _ij.t ′;

表示节点i与节点j之间的电压相角差；u(j)表示与节点j相连的上游节点集合，d(j)表示与节点j相连的下游节点集合，P_jl.t表示节点j与节点l间在t时刻的有功电量，Q_jl.t表示节点j与节点l间在t时刻的无功电量。Among them, P _it represents the active energy of node i at time t, Q _it represents the reactive energy of node i at time t, P _jt represents the active energy of node j at time t, Q _jt represents the reactive energy of node j at time t Electricity, U _it represents the voltage amplitude of node i at time t, U _jt represents the voltage amplitude of node j at time t, R _ij represents the resistance of the line between node i and node j, X _ij represents node i and node j The reactance of the line between nodes, G _ij represents the conductance of the line between node i and node j, B _ij represents the susceptance of the line between node i and node j;

Represents the voltage phase angle difference between node i and node j; u(j) represents the upstream node set connected to node j, d(j) represents the downstream node set connected to node j, and P _jl.t represents node j and The active power between node l at time t, Q _jl.t represents the reactive power between node j and node l at time t.

Preferably, in the step S200,

The constraints on the number of nodes connected to the distributed power supply equipment are preset as: N _{i.min ≤N i ≤N i.max} ;

Among them, N _i.min represents the minimum value of the number of distributed power supply devices connected to node i, N _i.max represents the maximum number of distributed power supply devices connected to the preset node i, and N _i represents the ith node connected to the number of distributed power equipment;

The upper limit constraint on the penetration rate of distributed power equipment is expressed as:

Among them, β represents the penetration rate after the distributed power equipment is connected to the node, P _load represents the power load of the distributed power equipment at the preset node; n _i represents a variable, which takes a value of 0 or 1, and n _i =0 means In the node distribution network structure, the ith node is not connected to the distributed power equipment, and n _i =1 means that the ith node is connected to the distributed power equipment; P _sg.DG represents the rated active power of each distributed power equipment;

The total power capacity constraint of distributed power equipment is expressed as: P _{min.DG ≤P t.DG ≤P max.DG} ;

Among them, P _min.DG represents the lower limit of the total power supply capacity of the distributed power equipment, P _max.DG represents the upper limit of the total power supply capacity of the distributed power equipment; P _t.DG represents the power supply provided by the distributed power equipment at time t Active power.

Preferably, in the step S200,

The expression of each line power flow constraint is:

表示节点i与节点j之间的电压相角差；u(j)表示与节点j相连的上游节点集合，d(j)表示与节点j相连的下游节点集合，P_jl.t表示节点j与节点l间在t时刻的有功电量，Q_jl.t表示节点j与节点l间在t时刻的无功电量；Among them, P _it represents the active energy of node i at time t, Q _it represents the reactive energy of node i at time t, P _jt represents the active energy of node j at time t, Q _jt represents the active energy of node j at time t Reactive power, P _ij.t represents the active power transmitted by branch ij at time t, Q _ij.t represents the reactive power transmitted by branch ij at time t, U _it represents the voltage amplitude of node i at time t, U _jt represents the voltage amplitude of node j at time t, R _ij represents the resistance of the line between node i and node j, X _ij represents the reactance of the line between node i and node j, G _ij represents the difference between node i and node j The conductance of the line between nodes, B _ij represents the susceptance of the line between node i and node j;

Represents the voltage phase angle difference between node i and node j; u(j) represents the upstream node set connected to node j, d(j) represents the downstream node set connected to node j, and P _jl.t represents node j and Active power between node l at time t, Q _jl.t represents the reactive power between node j and node l at time t;

The voltage constraint is: U _i.min ≤U _it ≤U _i.max , wherein U _i.min represents the minimum value of the voltage amplitude of node i, and U _i.max represents the maximum value of the voltage amplitude of node i;

The power constraint is: P _{ij.t ≤P ij.max} , where P _ij.t represents the power flowing on the line between node i and node j at time t, and P _ij.max represents node i and node j The maximum value of power flowing on the inter-line.

Preferably, in the step S200,

The target time is time t, and the expressions of the power load constraint transferred from the target time and the power load constraint transferred from the target time are:

Among them, ρ _min P _t.load represents the minimum value of the load transfer out power coefficient at time t, ρ _max P _t.load represents the maximum value of the load transfer out power coefficient at time t; σ _min P _t.load represents the load transfer input at time t The minimum value of the power coefficient, σ _max P _t.load represents the maximum value of the power coefficient transferred in at time t; P _t.out represents the amount of power load transferred at time t; P _t.in represents the power load transferred in at time t quantity.

Preferably, in step S500, based on the first functional constraint, the third functional constraint and the relaxed second functional constraint, the overall planning function is solved by a solver to obtain the node distribution network structure in the The steps of the power supply capacity of the node to be connected to the distributed power supply equipment include:

Based on the first function constraint, the third function constraint, and the relaxed second function constraint, the overall planning function is solved by the Cplex solver, and the corresponding nodes of the distributed power equipment to be connected in the node distribution network structure are obtained. The number of devices of distributed power equipment;

Obtain the power supply capacity of a distributed power supply device, multiply the power supply capacity of the distributed power supply device by the number of devices, and obtain the power supply capacity of the node to be connected to the distributed power supply device, wherein the power supply capacity of each distributed power supply device is The power supply capacity is the same.

The present invention also provides a configuration device for incremental distributed power equipment, and the configuration device for incremental distributed power equipment includes:

a building module for building an incremental distribution network model of distributed power sources in an incremental distribution network system, wherein the incremental distribution network model includes at least a node distribution network structure;

The acquiring module is configured to acquire, from the preset database of the incremental distribution network system, the first target data related to the distributed power supply node in the incremental distribution network model, and the data related to the incremental distribution network node. second target data and third target data associated with the power node;

The construction module is further configured to construct a first objective function corresponding to the distributed power supply node and a first function constraint that constrains the first objective function according to the first target data, and construct the first objective function according to the second target data. The second objective function corresponding to the incremental distribution network node and the second function constraint that constrains the second objective function, and the third objective function corresponding to the power node is constructed according to the third objective data and the constraint is described. The third function constraint of the third objective function, wherein the first function constraint includes a preset number of nodes connected to the distributed power equipment in the node distribution network structure, the upper limit of the penetration rate of the distributed power equipment, and the distributed power equipment. Total power supply capacity constraints, the second function constraints include power flow constraints, voltage constraints and power constraints of each line in the incremental distribution network model, and the third function constraints include power load constraints transferred at the target time and Constraints on the amount of power load transferred out at the target time;

The calculation module is used for multiplying the first objective function by the first weight corresponding to the first objective function to obtain a first weight function, and multiplying the second objective function by the first weight corresponding to the second objective function. two weights to obtain a second weight function, multiply the third objective function by the third weight corresponding to the third objective function to obtain a third weight function; combine the first weight function, the second weight function and the third weight function The three weight functions are added to obtain the overall planning function;

a relaxation processing module, configured to perform second-order cone relaxation processing on the second functional constraint to obtain a relaxed second functional constraint;

A solving module is configured to solve the overall planning function through a solver based on the first function constraint, the third function constraint and the second function constraint after relaxation processing, and obtain the distributed distribution to be accessed in the node distribution network structure The power capacity of the node of the power supply device;

The configuration module is configured to configure the distributed power supply device corresponding to the node according to the power supply capacity.

The present invention also provides a computer-readable storage medium, where a detection program is stored on the computer-readable storage medium, and when the detection program is executed by a processor, the steps of the configuration method for an incremental distributed power supply device as described above are implemented. .

Through the technical scheme of the present invention, the mutual coordination of each subject is realized, the feasibility of each decision is improved, the value of the overall planning function is optimized, and each node of the distributed power equipment is connected to the incremental distribution network system. It improves the rationality of the power supply capacity of each node in the incremental distribution network system, that is, improves the rationality of the configuration of the distributed power equipment corresponding to the nodes in the node distribution network structure, and ensures the incremental power distribution. Safe and stable operation of the network system.

Other beneficial effects of the present invention will be illustrated in the specific embodiments through the introduction of specific technical features and technical solutions. Those skilled in the art should be able to understand the technical features and technical solutions through the introduction of these technical features and technical solutions. beneficial technical effects.

Description of drawings

Preferred embodiments according to the present invention will be described below with reference to the accompanying drawings. In the picture:

FIG. 1 is a block diagram of a configuration method of an incremental distributed power supply device according to an embodiment of the present invention.

2 is a schematic diagram of an IEEE33 node distribution network structure in an embodiment of the present invention;

FIG. 3 is a schematic structural diagram of a configuration device for incremental distributed power equipment according to the present invention.

Detailed ways

The preferred embodiments of the present invention will now be described in detail with reference to the accompanying drawings.

It should be noted that, in the present invention, step numbers (letters or numbers) are used to refer to some specific method steps, which are only for the purpose of convenience and brevity of description, and are not intended to limit these method steps with letters or numbers. Order. Those skilled in the art can understand that the sequence of related method steps should be determined by the technology itself, and should not be unduly limited due to the existence of step numbers.

As shown in FIG. 1 , it is a flowchart of a method for configuring an incremental distributed power supply device according to an embodiment of the present invention.

S100. Build an incremental distribution network model of distributed power sources in an incremental distribution network system, where the incremental distribution network model includes at least a node distribution network structure.

The incremental distribution network model of distributed power generation is constructed in the incremental distribution network system. Among them, the incremental distribution network system is a distributed system, and the incremental distribution network model includes at least a node distribution network structure, which can be IEEE (Institute of Electrical and Electronics Engineers, Institute of Electrical and Electronics Engineers) 33 node distribution network structure. network structure, IEEE14 node distribution network structure or IEEE30 node distribution network structure, etc. Specifically, referring to FIG. 2 , FIG. 2 is a schematic diagram of an IEEE33 node distribution network structure in an embodiment of the present invention. As can be seen from FIG. 2 , in the node distribution network structure, node 5 , node 8 , node 14 , node 16 , node 18 and node 31 are the access nodes of the distributed power equipment. The incremental distribution network model also includes branch data information and node data information. The node distribution network structure describes the connection relationship between the nodes in the incremental distribution network model; the branch data information specifically refers to the impedance of the network lines in the incremental distribution network model; the node data information includes the The load size of the node, the total power capacity of the distributed power equipment connected to each node, and the node where each distributed power equipment is located. It can be understood that the node where the distributed power equipment is located is the installation location of the distributed power equipment. The distributed power equipment is an electric power equipment, specifically, a photovoltaic power generation equipment. The incremental distribution network system includes three stakeholders, namely DG operation users, distribution network investment users and power users. These three stakeholders all exist in the incremental distribution network system in the form of nodes, that is, incremental In the distribution network system, there are distributed power nodes corresponding to DG operating users, incremental distribution network nodes corresponding to distribution network investment users, and power nodes corresponding to power users. For the convenience of description, distributed power nodes, incremental distribution network nodes and power nodes are used to describe various stakeholders in the following.

S200: Acquire, from a preset database of the incremental power distribution network system, first target data related to the distributed power supply node and second target data related to the incremental power distribution network node in the incremental power distribution network model data and third target data related to the power node, construct a first target function corresponding to the distributed power supply node and a first function constraint that constrains the first target function according to the first target data, according to the first target function The second objective function corresponding to the incremental distribution network node and the second function constraint constraining the second objective function are constructed with the two objective data, and the third objective corresponding to the power node is constructed according to the third objective data Functions and constraints on the third function constraint of the third objective function, wherein the first functional constraints include the constraints on the number of preset nodes connected to distributed power equipment in the node distribution network structure, and the upper limit constraints on the penetration rate of distributed power equipment and the total power supply capacity constraints of distributed power equipment, the second function constraints include line power flow constraints, voltage constraints and power constraints in the incremental distribution network model, and the third function constraints include the transfer out of the target time. The power load constraints and the power load constraints transferred out at the target time.

In the incremental distribution network system, a database is preset, and the database stores the target data corresponding to the objective functions and function constraints of the distributed power supply nodes, the incremental distribution network nodes, and the power nodes. In this embodiment, The target data related to the distributed power supply node is recorded as the first target data, the target data related to the incremental distribution network node is recorded as the second target data, and the target data related to the power node is recorded as the third target data . When the corresponding objective function and function constraints need to be constructed, the first target data, the second target data and the third target data are obtained from the database, and the first target data is used to construct the first target function corresponding to the distributed power node and the constraint first target data. A first function constraint of an objective function, construct a second objective function corresponding to the incremental distribution network node and a second function constraint constraining the second objective function according to the second objective data, construct a first function corresponding to the power node according to the third objective data Three objective functions and a third function constraint that constrains the third objective function. Wherein, the first function constraint includes the preset number of nodes connected to the distributed power equipment in the node distribution network structure, the upper limit of the penetration rate of the distributed power equipment, and the total power capacity of the distributed power equipment, and the second function constraint includes the increment The power flow constraints, voltage constraints and power constraints of each line in the distribution network model, and the third function constraints include the power load constraints transferred at the target time and the power load constraints transferred at the target time. It should be noted that, in the process of constructing the objective function and function constraints, from the perspective of individual rationality of distributed power nodes, incremental distribution network nodes and power nodes, due to distributed power nodes, incremental distribution network nodes and power nodes. Nodes participating in the incremental distribution network model correspond to different goals in the planning and decision-making of the incremental distribution network system. Therefore, it is necessary to establish distributed power nodes, incremental distribution network nodes, and The objective function corresponding to the power node.

Specifically, in the step S200, the first objective function expression is formula (1):

Formula (1) maxW _DG (n _i , N _i )=W _S.DG -W _I.DG -W _OM.DG ;

Among them, the expression of W _S.DG in formula (1) is formula (2):

The expression of W _I.DG in formula (1) is formula (3):

The expression of W _OM.DG in formula (1) is formula (4):

n _i represents a variable, taking the value of 0 or 1, n _i =0 indicates that the ith node in the node distribution network structure is not connected to the distributed power equipment, and n _i =1 indicates that the ith node is connected to the distributed power equipment ; _{Ni represents the number of distributed power equipment connected to the i-th node; W S.DG represents} the electricity sales revenue of all distributed power equipment in the node distribution network structure within the preset years, W _I.DG represents the node distribution network structure W _OM.DG represents the operation and maintenance cost of all distributed power equipment within the preset years, U _t represents the collection at all times in a day, which is 24 hours; q _es represents distributed power The unit selling price of the equipment, the unit selling price of each distributed power equipment is the same; P _t.DG represents the active power provided by the distributed power equipment at time t; q _sg represents the construction cost of the distributed power equipment, specifically, q _sg represents the construction investment of the distributed power equipment corresponding to the unit power supply capacity; U _i represents the set of nodes that are preset to be connected to the distributed power equipment in the node distribution network structure; P _sg.DG represents the rated value of each distributed power equipment Active power, the rated active power of each distributed power equipment is the same; Y represents the life cycle of each distributed power equipment, and the life cycle of each distributed power equipment is the same; a represents the discount rate; q _om represents the distributed power equipment The operation and maintenance cost of the unit power supply. The size of the preset years can be set according to specific needs, for example, it can be set as 5 years, 10 years, or 12 years. The preliminary construction costs before the use of distributed power equipment include but are not limited to design costs, purchase costs and installation costs of distributed power equipment. It should be noted that the first objective function is the optimization objective function of the distributed power node, and the first objective function is used to reduce the initial construction investment cost and operation and maintenance expenditure of the distributed power equipment, so as to increase the electricity sales revenue of the distributed power node. , in order to maximize the net income of distributed power nodes.

The constraint on the number of pre-set nodes accessing distributed power equipment is expressed as formula (5): N _{i.min ≤N i ≤N i.max} .

Among them, N _i.min represents the minimum value of the number of distributed power supply devices connected to node i, N _i.max represents the maximum number of distributed power supply devices connected to the preset IEEE33 node i, and N _i represents the i-th node. The number of connected distributed power equipment.

The upper limit constraint of the penetration rate of distributed power equipment is expressed as formula (6):

Among them, β represents the penetration rate after the distributed power equipment is connected to the node, and the penetration rate is the ratio of the distributed power capacity to the rated capacity of the distribution network. P _load represents the power load of the distributed power supply equipment at the preset node; n _i represents a variable, which takes the value of 0 or 1, and n _i =0 indicates that the ith node in the node distribution network structure is not connected to the distributed power supply equipment , n _i =1 indicates that the ith node is connected to the distributed power equipment; P _sg.DG indicates the rated active power of each distributed power equipment.

The total power capacity constraint of distributed power equipment is expressed as formula (7): P _{min.DG ≤P t.DG ≤P max.DG} .

Among them, P _min.DG represents the lower limit of the total power supply capacity of the distributed power equipment, P _max.DG represents the upper limit of the total power supply capacity of the distributed power equipment; P _t.DG represents the power supply provided by the distributed power equipment at time t Active power. Each data in the first objective function and the first function constraint expression is the first objective data, and the first objective data can be obtained by the incremental distribution network system according to the historical data of the distributed power node, and stored in advance in the preset data. In the database, when needed, it can be obtained directly from the local database.

Specifically, in the step S200, the expression of the second objective function is formula (8):

Formula (8) _{maxW DN} (pi )=W _S.DN -W _L.DN -W _E.DN -W _B1.DN -W _B2.DN -W _B3.DG ;

Among them, the expression of W _S.DG in formula (8) is formula (9):

Formula (9):

The expression of W _L.DN in formula (8) is formula (10):

The expression of W _E.DN in formula (8) is formula (11):

formula (11)

The expression of W _B1.DN in formula (8) is formula (12):

formula (12)

The expression of W _B2.DN in formula (8) is formula (13):

The expression of W _B3.DG in Equation (8) is Equation (14):

Among them, pi represents the objective function variable of the incremental distribution network node, W _S.DN represents the electricity sales revenue of the incremental distribution network node within the preset period, and W _L.DN represents the line loss cost _within the preset period , W _E.DN represents the fault repair cost within the preset year, W _B1.DN represents the electricity purchase cost from the superior within the preset year, W _B2.DN represents the electricity purchase cost from the distributed power node within the preset year, W _B3.DG represents the loss caused by the output fluctuation of distributed power equipment, U _t represents the set of all times in a day, f _es represents the unit electricity selling price to the power node; P _t.load represents the initial power load at time t; P _t.out represents the amount of power load transferred at time t; P _t.in represents the amount of power load transferred at time t; P _t.is represents the amount of power load that can be interrupted at time t; P _t.loss represents the amount of power load transferred at time t Active power loss; EENS _t represents the expected value of less energy supplied to the power node at time t; U _b represents the set of all lines in the power grid in the incremental distribution network model; λ _b represents the probability of failure of the bth line; U _n represents the set of nodes that are preset to be connected to distributed power equipment in the node distribution network structure; P _ntload represents the initial power load of node n at time t; f _eb1 represents the unit power purchase price from the upper-level power grid; f _eb2 represents from The unit power purchase price of the distributed power node; P _p represents the punitive fee to be paid due to random fluctuations within the preset period; ρ _s represents the probability of power loss due to the output fluctuation of the distributed power equipment within the preset period; ΔQ _s represents Random fluctuation amount when the output of distributed power equipment fluctuates. It should be noted that the preset years in the first objective function and the second objective function are the same. Establish an incremental distribution network node optimization objective function, that is, build a second objective function, expect to reduce transmission line losses, early construction, and power purchases from superiors, increase electricity sales revenue, and maximize the benefits of incremental distribution network nodes as much as possible. ; and the uncertainty of the access nodes of distributed power equipment will affect the safety of the incremental distribution network system operation, resulting in a reduction in the economy of the optimal planning scheme of incremental distribution network nodes, so it is also expected to reduce the distribution of distributed power. Losses caused by fluctuations in the output of power equipment. In this embodiment, the loss caused by the output fluctuation of the distributed power supply equipment can be represented by the cost corresponding to the loss. The loss caused by the output fluctuation of the distributed power equipment is the loss caused by the unstable power supply of the distributed power equipment.

The expressions of the power flow constraints of each line are formula (15)-formula (19):

formula (15)

formula (16)

formula (17)

formula (18)

formula (19)

表示节点i与节点j之间的电压相角差；u(j)表示与节点j相连的上游节点集合，d(j)表示与节点j相连的下游节点集合，P_jl.t表示节点j与节点l间在t时刻的有功电量，Q_jl.t表示节点j与节点l间在t时刻的无功电量。Among them, P _it represents the active energy of node i at time t, Q _it represents the reactive energy of node i at time t, P _jt represents the active energy of node j at time t, Q _jt represents the active energy of node j at time t Reactive power, P _ij.t represents the active power transmitted by branch ij at time t, Q _ij.t represents the reactive power transmitted by branch ij at time t, U _it represents the voltage amplitude of node i at time t, U _jt represents the voltage amplitude of node j at time t, R _ij represents the resistance of the line between node i and node j, X _ij represents the reactance of the line between node i and node j, G _ij represents the difference between node i and node j The conductance of the line between nodes, B _ij represents the susceptance of the line between node i and node j;

Represents the voltage phase angle difference between node i and node j; u(j) represents the upstream node set connected to node j, d(j) represents the downstream node set connected to node j, and P _jl.t represents node j and The active power between node l at time t, Q _jl.t represents the reactive power between node j and node l at time t.

The expression of the voltage constraint is formula (20): U _i.min ≤U _it ≤U _i.max .

Among them, U _i.min represents the minimum value of the voltage amplitude of node i, and U _i.max represents the maximum value of the voltage amplitude of node i.

The expression of the power constraint is formula (21): P _{ij.t ≤P ij.max} .

Among them, P _ij.t represents the power flowing on the line between node i and node j at time t, and P _ij.max represents the maximum value of power flowing on the line between node i and node j.

It should be noted that, in this embodiment, after the distributed power equipment is connected to the node, problems such as power flow direction change, active management mode change, and increase of uncertain factors in the incremental distribution network system are considered. The loss caused by output fluctuations is added to the second objective function of , which reduces the adverse effects of uncertain factors after the distributed power equipment is connected to the node, that is, after entering the incremental distribution network system, thereby reducing its impact on the overall planning function. interference, which improves the accuracy of the resulting overall planning model.

Each data in the second objective function and the second function constraint expression is the second objective data, and the second objective data can be obtained by the incremental distribution network system according to the historical data of the incremental distribution network nodes, and stored in advance in the preset In a good database, when needed, it can be obtained directly from the local database.

Specifically, in the step S200, the expression of the third objective function is formula (22):

formula (22)

Among them, γ _eb represents the unit price of electricity purchased by the power node from the incremental distribution network node, P _t.out represents the power load transferred at time t; P _t.in represents the power load transferred at time t; P _{t .is} represents the power load that can be interrupted at time t; P _t.loss represents the active power loss at time t, and U _t represents the set of all moments in a day. The power load transferred out is the power load that is used more by the power node relative to the set reference power load; the power load transferred in is the power node that is used less relative to the set reference power load. Load; the power load of the interruption is the power load of the power node that is not useful when the power is interrupted. It should be noted that the objective function of the power node is established, and it is expected to reduce the electricity bill by adjusting the electricity consumption plan. The considered demand side response method is the price-based demand side response (DSR) based on the time-of-use electricity price. Participate in DSR According to the time-of-use electricity price information, the power nodes of the node remove the load during the peak period of the electricity price, and transfer the load during the trough period of the electricity price.

The target time is time t, and the expression of the power load amount constraint transferred from the target time and the power load amount constraint transferred from the target time is formula (23):

formula (23)

Among them, ρ _min P _t.load represents the minimum value of the load transfer out power coefficient at time t, ρ _max P _t.load represents the maximum value of the load transfer out power coefficient at time t; σ _min P _t.load represents the load transfer input at time t The minimum value of the power coefficient, σ _max P _t.load represents the maximum value of the power coefficient transferred in at time t; P _t.out represents the amount of power load transferred at time t; P _t.in represents the power load transferred in at time t quantity. Each data in the third objective function and the constraint expression of the third function is the third objective data. The third objective data can be obtained by the incremental distribution network system according to the historical data of the power nodes, and stored in the preset database in advance. , when needed, it can be obtained directly from the local database.

In this embodiment, each of the first objective function, the second objective function and the third objective function has a "max", in order to ensure the final value obtained by the first objective function, the second objective function and the third objective function are the largest, even if the benefits of distributed power nodes, incremental distribution network nodes and power nodes are maximized.

S300, multiply the first objective function by the first weight corresponding to the first objective function to obtain a first weight function, and multiply the second objective function by the second weight corresponding to the second objective function, A second weight function is obtained, and the third objective function is multiplied by a third weight corresponding to the third objective function to obtain a third weight function.

S400. Add the first weight function, the second weight function and the third weight function to obtain an overall planning function.

It is understandable that the three market players, namely distributed power nodes, incremental distribution network nodes and power nodes, have different interest demands, and their goals in participating in planning are different, and they all hope to maximize their respective benefits. The maximization of the individual interests of each subject may lead to an unfavorable situation in which the overall interests are far from the optimum. Therefore, this embodiment sets a weight for each objective function. The weight proportion and importance of the power grid node and the power node in the incremental power distribution network system are determined, and this embodiment does not limit the size of each weight. For the convenience of description, the first weight is denoted as r ₁ , the second weight is denoted as r ₂ , and the third weight is denoted as r ₃ , where r ₁ +r ₂ +r ₃ =1, and the overall planning function is denoted as W, Then the overall planning function can be expressed as formula (24):

Formula (24) W=γ ₁ W _DG +γ ₂ W _DN +γ ₃ W _US .

It should be noted that there is a "max" in the first objective function, the second objective function and the third objective function, so in the specific calculation in formula (24), W _DG , W _DN and W _US are also The maximum value is only in the formula (24), for the convenience of description, the "max" in front of the first objective function, the second objective function and the third objective function is omitted. It can be known from formula (24) that r ₁ W _DG is the first weight function, r ₂ W _DN is the second weight function, and r ₃ W _US is the third weight function.

S500. Perform second-order cone relaxation processing on the second functional constraint to obtain a relaxed second functional constraint, and based on the first functional constraint, the third functional constraint, and the relaxed second functional constraint, pass The solver solves the overall planning function, obtains the power supply capacity of the node to be connected to the distributed power supply device in the node distribution network structure, and configures the distributed power supply device corresponding to the node according to the power supply capacity and capacity.

The second-order cone relaxation process is performed on the second function constraint to obtain the second function constraint after the relaxation process. Specifically, the second-order cone relaxation process is performed on each line flow constraint in the second function constraint to obtain the second function constraint after relaxation, and then the first function constraint, the third function constraint and the second relaxation process are Under the condition of function constraints, the overall planning function is solved by the solver, and the power supply capacity of the node to be connected to the distributed power equipment in the node distribution network structure is obtained. In each node of the node distribution network structure, the position of each node can be At least one distributed power supply device is connected, but not all nodes in the node distribution network structure must be connected to the distributed power supply device. Specifically, the solver is a Cplex solver. In this embodiment, in order to improve the solution speed of solving the overall planning function, not all nodes in the node distribution network structure are regarded as the nodes to be connected to the distributed power supply equipment, but only a few of them are determined as the nodes to be connected. When the power supply capacity of the node to be connected to the distributed power equipment is 0, it indicates that the node to be connected to the distributed power equipment does not need to be connected to the distributed power equipment. When the power supply capacity of the distributed power equipment node is not 0, it indicates that the node to be connected to the distributed power equipment needs to be connected to the distributed power equipment.

Specifically, in the step S500, during the second-order cone relaxation process, auxiliary variables U _it '=(U _it ) ² and I _ij.t '=(I _ij.t ) ² are defined, and according to the auxiliary variables The variable yields the following target equation, whose expression is Equation (25):

formula (25)

Among them, I _ij.t represents the line current between node i and node j in the node distribution network structure, U _it represents the voltage amplitude of node i at time t, and P _ij.t represents the active power transmitted by branch ij at time t Electricity, Q _ij.t represents the reactive energy transmitted by branch ij at time t;

The objective equation is relaxed, and the following relaxation inequality is obtained, and the expression of the relaxation inequality is formula (26):

formula (26)

The relaxation inequality is reconstructed by the second-order cone, and the following target inequality is obtained. The expression of the target inequality is formula (27):

formula (27)

Perform second-order cone relaxation processing on each line power flow constraint in the second functional constraint according to the target inequality, to obtain the second functional constraint after relaxation processing;

The expressions of each line power flow constraint in the second function constraint after relaxation processing are formula (28)-formula (32):

formula (28)

formula (29)

formula(30)

formula (31)

Formula (32)(U _jt ) ² =U _it ′-2(R _ij P _ij.t +X _ij Q _ij.t )+[(R _ij ) ² +(X _ij ) ² ]·I _ij.t ′ .

表示节点i与节点j之间的电压相角差；u(j)表示与节点j相连的上游节点集合，d(j)表示与节点j相连的下游节点集合，P_jl.t表示节点j与节点l间在t时刻的有功电量，Q_jl.t表示节点j与节点l间在t时刻的无功电量。Among them, P _it represents the active energy of node i at time t, Q _it represents the reactive energy of node i at time t, P _jt represents the active energy of node j at time t, Q _jt represents the reactive energy of node j at time t Electricity, U _it represents the voltage amplitude of node i at time t, U _jt represents the voltage amplitude of node j at time t, R _ij represents the resistance of the line between node i and node j, X _ij represents node i and node j The reactance of the line between nodes, G _ij represents the conductance of the line between node i and node j, B _ij represents the susceptance of the line between node i and node j;

Represents the voltage phase angle difference between node i and node j; u(j) represents the upstream node set connected to node j, d(j) represents the downstream node set connected to node j, and P _jl.t represents node j and The active power between node l at time t, Q _jl.t represents the reactive power between node j and node l at time t.

In the second function constraint after relaxation, formula (28) is the same as formula (15), and formula (29) is the same as formula (16), that is, during the relaxation process, formula (15) is not and Equation (16) for processing.

It should be noted that when solving the overall planning function through the solver, it is required to obtain the maximum value of the overall planning function, that is, maxW=r ₁ W _DG +r ₂ W _DN +r ₃ W _US . By solving the overall planning function with the solver, the constraints on the power supply capacity of the nodes to be connected to the distributed power equipment in the node distribution network structure can be expressed as:

It should be noted that, after determining the power supply capacity of each distributed power supply device node to be connected, the power supply capacity of each distributed power supply device is obtained, and the power supply capacity of each distributed power supply device node to be connected and each distributed power supply The power supply capacity of the device can determine the number of distributed power supply devices connected to the node, and then configure the corresponding distributed power supply device to the node, that is, configure the distributed power supply device corresponding to the node according to the power supply capacity. In this embodiment, the sum of the power supply capacity of the distributed power supply equipment connected to the node is equal to or greater than the power supply capacity of the node to be connected to the distributed power supply equipment. The power capacity of the distributed power equipment connected to the same node can be the same or different.

In order to facilitate understanding, for example, if the IEEE33 node distribution network example is used, the number of branches included in this incremental distribution network system is 33, and the simulation period is 10 years, that is, the preset period is 10 years. The nodes of the distributed power supply equipment are node 5 , node 8 , node 14 , node 16 , node 18 and node 31 . By solving the overall planning function, it can be determined that the distributed power equipment should be connected to the nodes among node 5, node 8, node 14, node 16, node 18 and node 31. , node 18 and node 31, not all nodes are necessarily connected to distributed power equipment, and the number of distributed power equipment connected to each node may or may not be the same.

Using the solver to calculate the maximum overall market revenue within the simulation period is 1,066,150 yuan, that is, the maximum value of the overall planning function is 1,066,150 yuan. When the overall planning function is equal to the maximum value, the node locations and power supply capacity of distributed power equipment are as follows The table shows:

The values of distributed power nodes, incremental distribution network nodes and power nodes (revenue of each subject) and the maximum value of the overall planning function (overall benefit) are shown in the following table:

Further, in order to improve the accuracy of the calculation of the maximum value of the overall planning function and the accuracy of the power node revenue, in the process of calculating the maximum value of the overall planning function, the unit electricity price in different time periods should be considered, such as during peak electricity consumption. During the period of 8:00-11:00&18:00-23:00 every day, the electricity purchase price is 800 (yuan/MW·h); in the low electricity consumption period, every day 6:00-8:00&11:00-18 During the period from 23:00 to 6:00, the electricity purchase price is 400 (yuan/MW·h); during the period of 23:00-6:00 every day, the electricity purchase price is 600 (yuan/MW·h).

This embodiment realizes the mutual coordination of each subject, improves the feasibility of their respective decision-making, optimizes the value of the overall planning function, and enables the power supply capacity configuration of each node connected to the distributed power supply equipment in the incremental distribution network system Optimal, which improves the rationality of the power supply capacity of each node in the incremental distribution network system, that is, improves the rationality of the configuration of the distributed power equipment corresponding to the nodes in the node distribution network structure, and ensures the safety of the incremental distribution network system. Stable operation.

Further, another embodiment of the configuration method of the incremental distributed power equipment of the present invention is proposed.

The difference between another embodiment of the method for configuring an incremental distributed power supply device and the above embodiment of the method for configuring an incremental distributed power supply device is that in step S500, the third The step of obtaining the power supply capacity of the node to be connected to the distributed power supply equipment in the node distribution network structure includes:

Step a, based on the first function constraint, the third function constraint and the relaxed second function constraint, solve the overall planning function through the Cplex solver, and obtain the distributed power supply to be connected in the node distribution network structure The device number of the device node corresponding to the distributed power supply device;

Step b: Obtain the power supply capacity of a distributed power supply device, multiply the power supply capacity of the distributed power supply device by the number of devices, and obtain the power supply capacity of the node to be connected to the distributed power supply device, wherein each distributed power supply device is The power supply capacity of the power supply unit is the same.

When the second function constraint after relaxation processing is obtained, based on the first function constraint, the third function constraint and the second function constraint after relaxation processing, the overall planning function is solved by the Cplex solver, and the node distribution network structure to be The number of devices connected to the distributed power supply equipment node corresponds to the number of distributed power supply equipment, and then obtains the power supply capacity of a distributed power supply device, and multiplies the power supply capacity of the distributed power supply equipment by the number of devices corresponding to each distributed power supply equipment node to be connected. , to obtain the power supply capacity of the node to be connected to the distributed power supply equipment, wherein the power supply capacity of each distributed power supply equipment is the same.

Further, if the power supply capacity of each distributed power supply device is different, the power supply capacity of each distributed power supply device can also be obtained to determine which distributed power supply device needs to be connected according to the power supply capacity of the node to be connected to the distributed power supply device. The power supply device is connected to the node to be connected to the distributed power supply device.

In this embodiment, the power supply capacity of each distributed power supply device is the same, and the speed of obtaining the power supply capacity of the node to be connected to the distributed power supply device in the node distribution network structure is improved.

The present invention also provides a configuration device for incremental distributed power equipment. Referring to FIG. 3 , the configuration device for incremental distributed power equipment includes:

A building module 10, configured to build an incremental distribution network model of distributed power sources in an incremental distribution network system, wherein the incremental distribution network model includes at least a node distribution network structure;

The obtaining module 20 is configured to obtain, from the preset database of the incremental power distribution network system, the first target data related to the distributed power supply nodes in the incremental power distribution network model, and the first target data related to the incremental power distribution network nodes The second target data and the third target data related to the power node;

The construction module 10 is further configured to construct, according to the first target data, a first objective function corresponding to the distributed power supply node and a first function constraint that constrains the first objective function, and construct a first objective function based on the second target data. The second objective function corresponding to the incremental distribution network node and the second function constraint that constrains the second objective function, and the third objective function corresponding to the power node and the constraint are constructed according to the third objective data. The third function constraints of the third objective function, wherein the first function constraints include the constraints on the number of preset nodes connected to distributed power equipment in the node distribution network structure, the upper limit constraint on the penetration rate of distributed power equipment, and the distributed power equipment. The total power supply capacity constraint of the equipment, the second function constraint includes each line power flow constraint, voltage constraint and power constraint in the incremental distribution network model, and the third function constraint includes the power load amount transferred out at the target time constraint. and the power load constraints transferred out at the target time;

The calculation module 30 is used to multiply the first objective function by the first weight corresponding to the first objective function to obtain a first weight function, and multiply the second objective function by the corresponding second objective function. The second weight is obtained to obtain a second weight function, and the third objective function is multiplied by the third weight corresponding to the third objective function to obtain a third weight function; the first weight function, the second weight function and the The third weight function is added to obtain the overall planning function;

A relaxation processing module 40, configured to perform second-order cone relaxation processing on the second functional constraint to obtain a relaxed second functional constraint;

A solving module 50 is configured to solve the overall planning function through a solver based on the first functional constraint, the third functional constraint and the relaxed second functional constraint, and obtain the distribution to be accessed in the node distribution network structure The power capacity of the node of the type power supply device;

The configuration module 60 is configured to configure the distributed power supply device corresponding to the node according to the power supply capacity.

Further, in the second-order cone relaxation process, auxiliary variables U _it ′=(U _it ) ² and I _ij.t ′=(I _ij.t ) ² are defined, and the following objective equation is obtained according to the auxiliary variables:

I _ij.t represents the line current between node i and node j in the node distribution network structure, U _it represents the voltage amplitude of node i at time t, P _ij.t represents the active power transmitted by branch ij at time t, Q _ij.t represents the reactive power transmitted by branch ij at time t;

The objective equation is relaxed to obtain the following relaxed inequality:

Reconstructing the relaxation inequalities with second-order cones yields the following objective inequalities:

Perform second-order cone relaxation processing on each line power flow constraint in the second functional constraint according to the target inequality, to obtain the second functional constraint after relaxation processing;

The expression of each line power flow constraint in the second function constraint after relaxation is:

表示节点i与节点j之间的电压相角差；u(j)表示与节点j相连的上游节点集合，d(j)表示与节点j相连的下游节点集合，P_jl.t表示节点j与节点l间在t时刻的有功电量，Q_jl.t表示节点j与节点l间在t时刻的无功电量。Among them, P _it represents the active energy of node i at time t, Q _it represents the reactive energy of node i at time t, P _jt represents the active energy of node j at time t, Q _jt represents the reactive energy of node j at time t Electricity, U _it represents the voltage amplitude of node i at time t, U _jt represents the voltage amplitude of node j at time t, R _ij represents the resistance of the line between node i and node j, X _ij represents node i and node j The reactance of the line between nodes, G _ij represents the conductance of the line between node i and node j, B _ij represents the susceptance of the line between node i and node j;

Represents the voltage phase angle difference between node i and node j; u(j) represents the upstream node set connected to node j, d(j) represents the downstream node set connected to node j, and P _jl.t represents node j and The active power between node l at time t, Q _jl.t represents the reactive power between node j and node l at time t.

Further, the preset constraints on the number of nodes accessing distributed power equipment are expressed as: N _{i.min ≤N i ≤N i.max} ;

Among them, N _i.min represents the minimum value of the number of distributed power supply devices connected to node i, N _i.max represents the maximum number of distributed power supply devices connected to the preset node i, and N _i represents the ith node connected to the number of distributed power equipment;

The upper limit constraint on the penetration rate of distributed power equipment is expressed as:

Among them, β represents the penetration rate after the distributed power equipment is connected to the node, P _load represents the power load of the distributed power equipment at the preset node; n _i represents a variable, which takes a value of 0 or 1, and n _i =0 means In the node distribution network structure, the ith node is not connected to the distributed power equipment, and n _i =1 means that the ith node is connected to the distributed power equipment; P _sg.DG represents the rated active power of each distributed power equipment;

The total power capacity constraint of distributed power equipment is expressed as: P _{min.DG ≤P t.DG ≤P max.DG} ;

Among them, P _min.DG represents the lower limit of the total power supply capacity of the distributed power equipment, P _max.DG represents the upper limit of the total power supply capacity of the distributed power equipment; P _t.DG represents the power supply provided by the distributed power equipment at time t Active power.

Further, the expression of each line power flow constraint is:

表示节点i与节点j之间的电压相角差；u(j)表示与节点j相连的上游节点集合，d(j)表示与节点j相连的下游节点集合，P_jl.t表示节点j与节点l间在t时刻的有功电量，Q_jl.t表示节点j与节点l间在t时刻的无功电量；Among them, P _it represents the active energy of node i at time t, Q _it represents the reactive energy of node i at time t, P _jt represents the active energy of node j at time t, Q _jt represents the active energy of node j at time t Reactive power, P _ij.t represents the active power transmitted by branch ij at time t, Q _ij.t represents the reactive power transmitted by branch ij at time t, U _it represents the voltage amplitude of node i at time t, U _jt represents the voltage amplitude of node j at time t, R _ij represents the resistance of the line between node i and node j, X _ij represents the reactance of the line between node i and node j, G _ij represents the difference between node i and node j The conductance of the line between nodes, B _ij represents the susceptance of the line between node i and node j;

Represents the voltage phase angle difference between node i and node j; u(j) represents the upstream node set connected to node j, d(j) represents the downstream node set connected to node j, and P _jl.t represents node j and Active power between node l at time t, Q _jl.t represents the reactive power between node j and node l at time t;

The voltage constraint is: U _i.min ≤U _it ≤U _i.max , wherein U _i.min represents the minimum value of the voltage amplitude of node i, and U _i.max represents the maximum value of the voltage amplitude of node i;

The power constraint is: P _{ij.t ≤P ij.max} , where P _ij.t represents the power flowing on the line between node i and node j at time t, and P _ij.max represents node i and node j The maximum value of power flowing on the inter-line.

Further, the target time is time t, and the expressions of the power load constraint transferred from the target time and the power load constraint transferred from the target time are:

Among them, ρ _min P _t.load represents the minimum value of the load transfer out power coefficient at time t, ρ _max P _t.load represents the maximum value of the load transfer out power coefficient at time t; σ _min P _t.load represents the load transfer input at time t The minimum value of the power coefficient, σ _max P _t.load represents the maximum value of the power coefficient transferred in at time t; P _t.out represents the amount of power load transferred at time t; P _t.in represents the power load transferred in at time t quantity.

Further, the solving module 50 includes:

A solving unit, configured to solve the overall planning function through the Cplex solver based on the first function constraint, the third function constraint and the second function constraint after relaxation processing, to obtain the distribution to be accessed in the node distribution network structure The number of distributed power equipment nodes corresponding to distributed power equipment;

The acquisition unit is used to acquire the power capacity of a distributed power supply device;

The computing unit is configured to multiply the power supply capacity of the distributed power supply device by the number of devices to obtain the power supply capacity of the node to be connected to the distributed power supply device, wherein the power supply capacity of each distributed power supply device is the same.

The specific implementations of the configuration apparatus for incremental distributed power equipment according to the present invention are basically the same as the above-mentioned embodiments of the configuration method for incremental distributed power equipment, and will not be repeated here.

The present invention also provides a computer-readable storage medium, where a detection program is stored on the computer-readable storage medium, and when the detection program is executed by a processor, the steps of the above-mentioned method for configuring an incremental distributed power supply device are implemented.

The specific implementation manner of the computer-readable storage medium of the present invention is basically the same as the above-mentioned embodiments of the configuration method of the incremental distributed power supply device, and will not be repeated here.

Those skilled in the art can understand that, under the premise of no conflict, the above preferred solutions can be freely combined and superimposed.

It should be understood that the above-mentioned embodiments are only exemplary rather than restrictive, and those skilled in the art can make various obvious or equivalent to the above-mentioned details without departing from the basic principles of the present invention. Modifications or substitutions will be included within the scope of the claims of the present invention.

Claims (10)

1. a configuration method of incremental distributed power supply equipment, is characterized in that, the configuration method of described incremental distributed power supply equipment comprises the following steps: S100. Build an incremental distribution network model of distributed power sources in an incremental distribution network system, wherein the incremental distribution network model includes at least a node distribution network structure; S200: Acquire, from a preset database of the incremental power distribution network system, first target data related to the distributed power supply node and second target data related to the incremental power distribution network node in the incremental power distribution network model data and third target data related to the power node, construct a first target function corresponding to the distributed power supply node and a first function constraint that constrains the first target function according to the first target data, according to the first target function The second objective function corresponding to the incremental distribution network node and the second function constraint constraining the second objective function are constructed with the two objective data, and the third objective corresponding to the power node is constructed according to the third objective data Functions and constraints on the third function constraint of the third objective function, wherein the first functional constraints include the constraints on the number of preset nodes connected to distributed power equipment in the node distribution network structure, and the upper limit constraints on the penetration rate of distributed power equipment and the total power supply capacity constraints of distributed power equipment, the second function constraints include line power flow constraints, voltage constraints and power constraints in the incremental distribution network model, and the third function constraints include the transfer out of the target time. The power load amount constraint and the power load amount constraint transferred at the target time; wherein, in the step S200, the expression of the first objective function is the formula: max W _DG (n _i , N _i )=W _S.DG -W _I.DG -W _OM.DG ; Among them, the expression of W _S.DG in the formula is the formula:

The expression of W _I.DG in the formula is the formula:

The expression of W _OM.DG in the formula is the formula:

n _i represents a variable, taking the value of 0 or 1, n _i =0 indicates that the ith node in the node distribution network structure is not connected to the distributed power equipment, and n _i =1 indicates that the ith node is connected to the distributed power equipment ; _{Ni represents the number of distributed power equipment connected to the i-th node; W S.DG represents} the electricity sales revenue of all distributed power equipment in the node distribution network structure within the preset years, W _I.DG represents the node distribution network structure W _OM.DG represents the operation and maintenance cost of all distributed power equipment within the preset years, U _t represents the collection at all times in a day, which is 24 hours; q _es represents distributed power The unit selling price of the equipment, the unit selling price of each distributed power equipment is the same; P _t.DG represents the active power provided by the distributed power equipment at time t; q _sg represents the construction cost of the distributed power equipment, specifically, q _sg represents the construction investment of the distributed power equipment corresponding to the unit power supply capacity; U _i represents the set of nodes that are preset to be connected to the distributed power equipment in the node distribution network structure; P _sg.DG represents the rated value of each distributed power equipment Active power, the rated active power of each distributed power equipment is the same; Y represents the life cycle of each distributed power equipment, and the life cycle of each distributed power equipment is the same; a represents the discount rate; q _om represents the distributed power equipment The operation and maintenance cost per unit of power supply; wherein, in the step S200, the expression of the second objective function is the formula: Formula max W _DN (pi )=W _S.DN -W _L.DN -W _E.DN -W _B1.DN -W _{B2.DN -W B3.DG ;} Among them, the expression of W _S.DG in the formula is the formula: formula:

The expression of W _L.DN in the formula is the formula:

Among them, pi represents the objective function variable of the incremental distribution network node, W _S.DN represents the electricity sales revenue of the incremental distribution network node within the preset period, and W _L.DN represents the line loss cost _within the preset period , W _E.DN represents the fault repair cost within the preset year, W _B1.DN represents the electricity purchase cost from the superior within the preset year, W _B2.DN represents the electricity purchase cost from the distributed power node within the preset year, W _B3.DG represents the loss caused by the output fluctuation of distributed power equipment, U _t represents the set of all times in a day, f _es represents the unit electricity selling price to the power node; P _t.load represents the initial power load at time t; P _t.out represents the amount of power load transferred at time t; P _t.in represents the amount of power load transferred at time t; P _t.is represents the amount of power load that can be interrupted at time t; P _t.loss represents the amount of power load transferred at time t Active power loss; EENS _t represents the expected value of less energy supplied to the power node at time t; U _b represents the set of all lines in the power grid in the incremental distribution network model; λ _b represents the probability of failure of the bth line; U _n represents the set of nodes that are preset to be connected to distributed power equipment in the node distribution network structure; P _ntload represents the initial power load of node n at time t; f _eb1 represents the unit power purchase price from the upper-level power grid; f _eb2 represents from The unit power purchase price of the distributed power node; P _p represents the penalty fee to be paid due to random fluctuations within the preset period; ρ _s represents the probability of power loss caused by the output fluctuation of the distributed power equipment within the preset period; ΔQ _s represents Random fluctuations when the output of distributed power equipment fluctuates; The expression of W _E.DN in the formula is the formula:

The expression of W _B1.DN in the formula is the formula:

The expression of W _B2.DN in the formula is the formula:

The expression of W _B3.DG in the formula is the formula:

Wherein, in the step S200, the expression of the third objective function is the formula:

Among them, γ _eb represents the unit price of electricity purchased by the power node from the incremental distribution network node, P _t.out represents the power load transferred at time t; P _t.in represents the power load transferred at time t; P _{t .is} represents the power load that can be interrupted at time t; P _t.loss represents the active power loss at time t, and U _t represents the set of all moments in a day; S300, multiply the first objective function by the first weight corresponding to the first objective function to obtain a first weight function, and multiply the second objective function by the second weight corresponding to the second objective function, Obtain the second weight function, multiply the third objective function by the third weight corresponding to the third objective function, and obtain the third weight function; S400, adding the first weight function, the second weight function and the third weight function to obtain an overall planning function; S500. Perform second-order cone relaxation processing on the second functional constraint to obtain a relaxed second functional constraint, and based on the first functional constraint, the third functional constraint, and the relaxed second functional constraint, pass The solver solves the overall planning function, obtains the power supply capacity of the node to be connected to the distributed power supply device in the node distribution network structure, and configures the distributed power supply device corresponding to the node according to the power supply capacity. 2 . The method for configuring an incremental distributed power supply device according to claim 1 , wherein, in the step S500 , in the second-order cone relaxation process, an auxiliary variable U _it ′=(U _it ) is defined. 3 . ² and I _ij.t ′=(I _ij.t ) ² , the following objective equation is obtained according to the auxiliary variable:

I _ij.t represents the line current between node i and node j in the node distribution network structure, U _it represents the voltage amplitude of node i at time t, P _ij.t represents the active power transmitted by branch ij at time t, Q _ij.t represents the reactive power transmitted by branch ij at time t; The objective equation is relaxed to obtain the following relaxed inequality:

Reconstructing the relaxation inequalities with second-order cones yields the following objective inequalities:

Perform second-order cone relaxation processing on each line power flow constraint in the second functional constraint according to the target inequality, to obtain the second functional constraint after relaxation processing; The expression of each line power flow constraint in the second function constraint after relaxation is:

(U _jt ) ² =U _it ′-2(R _ij P _ij.t +X _ij Q _ij.t )+[(R _ij ) ² +(X _ij ) ² ]·I _ij.t ′; Among them, P _it represents the active energy of node i at time t, Q _it represents the reactive energy of node i at time t, P _jt represents the active energy of node j at time t, Q _jt represents the reactive energy of node j at time t Electricity, U _it represents the voltage amplitude of node i at time t, U _jt represents the voltage amplitude of node j at time t, R _ij represents the resistance of the line between node i and node j, X _ij represents node i and node j The reactance of the line between nodes, G _ij represents the conductance of the line between node i and node j, B _ij represents the susceptance of the line between node i and node j;

Represents the voltage phase angle difference between node i and node j; u(j) represents the upstream node set connected to node j, d(j) represents the downstream node set connected to node j, and P _jl.t represents node j and The active power between node l at time t, Q _jl.t represents the reactive power between node j and node l at time t. 3. The method for configuring an incremental distributed power supply device according to claim 1, wherein in the step S200, The constraints on the number of nodes connected to the distributed power supply equipment are preset as: N _{i.min ≤N i ≤N i.max} ; Among them, N _i.min represents the minimum value of the number of distributed power supply devices connected to node i, N _i.max represents the maximum number of distributed power supply devices connected to the preset node i, and N _i represents the ith node connected to the number of distributed power equipment; The upper limit constraint on the penetration rate of distributed power equipment is expressed as:

Among them, β represents the penetration rate after the distributed power equipment is connected to the node, P _load represents the power load of the distributed power equipment at the preset node; n _i represents a variable, which takes a value of 0 or 1, and n _i =0 means In the node distribution network structure, the ith node is not connected to the distributed power equipment, and n _i = 1 means that the ith node is connected to the distributed power equipment; P _sg .DG represents the rated active power of each distributed power equipment; The total power capacity constraint of distributed power equipment is expressed as: P _{min.DG ≤P t.DG ≤P max.DG} ; Among them, P _min.DG represents the lower limit of the total power supply capacity of the distributed power equipment, P _max.DG represents the upper limit of the total power supply capacity of the distributed power equipment; P _t.DG represents the power supply provided by the distributed power equipment at time t Active power. 4. The method for configuring an incremental distributed power supply device according to claim 1, wherein in the step S200, The expression of each line power flow constraint is:

Among them, P _it represents the active energy of node i at time t, Q _it represents the reactive energy of node i at time t, P _jt represents the active energy of node j at time t, Q _jt represents the active energy of node j at time t Reactive power, P _ij.t represents the active power transmitted by branch ij at time t, Q _ij.t represents the reactive power transmitted by branch ij at time t, U _it represents the voltage amplitude of node i at time t, U _jt represents the voltage amplitude of node j at time t, R _ij represents the resistance of the line between node i and node j, X _ij represents the reactance of the line between node i and node j, G _ij represents the difference between node i and node j The conductance of the line between nodes, B _ij represents the susceptance of the line between node i and node j;

Represents the voltage phase angle difference between node i and node j; u(j) represents the upstream node set connected to node j, d(j) represents the downstream node set connected to node j, and P _jl.t represents node j and Active power between node l at time t, Q _jl.t represents the reactive power between node j and node l at time t; The voltage constraint is: U _i.min ≤U _it ≤U _i.max , wherein U _i.min represents the minimum value of the voltage amplitude of node i, and U _i.max represents the maximum value of the voltage amplitude of node i; The power constraint is: P _{ij.t ≤P ij.max} , where P _ij.t represents the power flowing on the line between node i and node j at time t, and P _ij.max represents node i and node j The maximum value of power flowing on the inter-line. 5. The method for configuring an incremental distributed power supply device according to any one of claims 1 to 4, wherein in the step S200, The target time is time t, and the expressions of the power load constraint transferred from the target time and the power load constraint transferred from the target time are:

Among them, ρ _min P _t.load represents the minimum value of the load transfer out power coefficient at time t, ρ _max P _t.load represents the maximum value of the load transfer out power coefficient at time t; σ _min P _t.load represents the load transfer input at time t The minimum value of the power coefficient, σ _max P _t.load represents the maximum value of the power coefficient transferred in at time t; P _t.out represents the amount of power load transferred at time t; P _t.in represents the power load transferred in at time t quantity. 6 . The method for configuring an incremental distributed power supply device according to claim 5 , wherein, in step S500 , the second function based on the first function constraint, the third function constraint and the relaxation processed Constraints, the steps of obtaining the power supply capacity of the nodes to be connected to the distributed power equipment in the node distribution network structure by solving the overall planning function by a solver include: Based on the first function constraint, the third function constraint, and the relaxed second function constraint, the overall planning function is solved by the Cplex solver, and the corresponding nodes of the distributed power equipment to be connected in the node distribution network structure are obtained. The number of devices of distributed power equipment; Obtain the power supply capacity of a distributed power supply device, multiply the power supply capacity of the distributed power supply device by the number of devices, and obtain the power supply capacity of the node to be connected to the distributed power supply device, wherein the power supply capacity of each distributed power supply device is The power supply capacity is the same. 7. A configuration device for incremental distributed power supply equipment, wherein the configuration device for incremental distributed power supply equipment comprises: a building module for building an incremental distribution network model of distributed power sources in an incremental distribution network system, wherein the incremental distribution network model includes at least a node distribution network structure; The acquiring module is configured to acquire, from the preset database of the incremental distribution network system, the first target data related to the distributed power supply node in the incremental distribution network model, and the data related to the incremental distribution network node. second target data and third target data associated with the power node; The construction module is further configured to construct a first objective function corresponding to the distributed power supply node and a first function constraint that constrains the first objective function according to the first target data, and construct the first objective function according to the second target data. The second objective function corresponding to the incremental distribution network node and the second function constraint that constrains the second objective function, and the third objective function corresponding to the power node is constructed according to the third objective data and the constraint is described. The third function constraint of the third objective function, wherein the first function constraint includes a preset number of nodes connected to the distributed power equipment in the node distribution network structure, the upper limit of the penetration rate of the distributed power equipment, and the distributed power equipment. Total power supply capacity constraints, the second function constraints include power flow constraints, voltage constraints and power constraints of each line in the incremental distribution network model, and the third function constraints include power load constraints transferred at the target time and The power load amount constraint transferred out at the target time; wherein, the expression of the first objective function is the formula: Formula max W _DG (n _i ,N _i )=W _S.DG -W _I.DG -W _OM.DG ; Among them, the expression of W _S.DG in the formula is the formula:

The expression of W _I.DG in the formula is the formula:

The expression of W _OM.DG in the formula is the formula:

n _i represents a variable, taking the value of 0 or 1, n _i =0 indicates that the ith node in the node distribution network structure is not connected to the distributed power equipment, and n _i =1 indicates that the ith node is connected to the distributed power equipment ; _{Ni represents the number of distributed power equipment connected to the i-th node; W S.DG represents} the electricity sales revenue of all distributed power equipment in the node distribution network structure within the preset years, W _I.DG represents the node distribution network structure W _OM.DG represents the operation and maintenance cost of all distributed power equipment within the preset years, U _t represents the collection at all times in a day, which is 24 hours; q _es represents distributed power The unit selling price of the equipment, the unit selling price of each distributed power equipment is the same; P _t.DG represents the active power provided by the distributed power equipment at time t; q _sg represents the construction cost of the distributed power equipment, specifically, q _sg represents the construction investment of the distributed power equipment corresponding to the unit power supply capacity; U _i represents the set of nodes that are preset to be connected to the distributed power equipment in the node distribution network structure; P _sg.DG represents the rated value of each distributed power equipment Active power, the rated active power of each distributed power equipment is the same; Y represents the life cycle of each distributed power equipment, and the life cycle of each distributed power equipment is the same; a represents the discount rate; q _om represents the distributed power equipment The operation and maintenance cost of the unit power supply; wherein, in the step S200, the expression of the second objective function is the formula: Formula max W _DN (pi )=W _S.DN -W _L.DN -W _E.DN -W _B1.DN -W _{B2.DN -W B3.DG ;} Among them, the expression of W _S.DG in the formula is the formula: formula:

The expression of W _L.DN in the formula is the formula:

Among them, pi represents the objective function variable of the incremental distribution network node, W _S.DN represents the electricity sales revenue of the incremental distribution network node within the preset period, and W _L.DN represents the line loss cost _within the preset period , W _E.DN represents the fault repair cost within the preset year, W _B1.DN represents the electricity purchase cost from the superior within the preset year, W _B2.DN represents the electricity purchase cost from the distributed power node within the preset year, W _B3.DG represents the loss caused by the output fluctuation of distributed power equipment, U _t represents the set of all times in a day, f _es represents the unit electricity selling price to the power node; P _t.load represents the initial power load at time t; P _t.out represents the amount of power load transferred at time t; P _t.in represents the amount of power load transferred at time t; P _t.is represents the amount of power load that can be interrupted at time t; P _t.loss represents the amount of power load transferred at time t Active power loss; EENS _t represents the expected value of less energy supplied to the power node at time t; U _b represents the set of all lines in the power grid in the incremental distribution network model; λ _b represents the probability of failure of the bth line; U _n represents the set of nodes that are preset to be connected to distributed power equipment in the node distribution network structure; P _ntload represents the initial power load of node n at time t; f _eb1 represents the unit power purchase price from the upper-level power grid; f _eb2 represents from The unit power purchase price of the distributed power node; P _p represents the penalty fee to be paid due to random fluctuations within the preset period; ρ _s represents the probability of power loss caused by the output fluctuation of the distributed power equipment within the preset period; ΔQ _s represents Random fluctuations when the output of distributed power equipment fluctuates; The expression of W _E.DN in the formula is the formula: formula

The expression of W _B1.DN in the formula is the formula: formula

The expression of W _B2.DN in the formula is the formula:

The expression of W _B3.DG in the formula is the formula:

Wherein, in the step S200, the expression of the third objective function is the formula: formula

Among them, γ _eb represents the unit price of electricity purchased by the power node from the incremental distribution network node, P _t.out represents the power load transferred at time t; P _t.in represents the power load transferred at time t; P _{t .is} represents the power load that can be interrupted at time t; P _t.loss represents the active power loss at time t, and U _t represents the set of all moments in a day; The calculation module is used for multiplying the first objective function by the first weight corresponding to the first objective function to obtain a first weight function, and multiplying the second objective function by the first weight corresponding to the second objective function. two weights to obtain a second weight function, multiply the third objective function by the third weight corresponding to the third objective function to obtain a third weight function; combine the first weight function, the second weight function and the third weight function The three weight functions are added to obtain the overall planning function; a relaxation processing module, configured to perform second-order cone relaxation processing on the second functional constraint to obtain a relaxed second functional constraint; A solving module is used to solve the overall planning function through a solver based on the first function constraint, the third function constraint and the second function constraint after relaxation processing, and obtain the distribution network structure to be accessed in the node distribution network structure. The power capacity of the node of the power supply device; The configuration module is configured to configure the distributed power supply device corresponding to the node according to the power supply capacity. 8 . The configuration device for incremental distributed power equipment according to claim 7 , wherein, in the second-order cone relaxation process, auxiliary variables U _it ′=(U _it ) ² and I _ij.t ′ are defined. 9 . =(I _ij.t ) ² , the following objective equation is obtained according to the auxiliary variable:

I _ij.t represents the line current between node i and node j in the node distribution network structure, U _it represents the voltage amplitude of node i at time t, P _ij.t represents the active power transmitted by branch ij at time t, Q _ij.t represents the reactive power transmitted by branch ij at time t; The objective equation is relaxed to obtain the following relaxed inequality:

Reconstructing the relaxation inequalities with second-order cones yields the following objective inequalities:

Perform second-order cone relaxation processing on each line power flow constraint in the second functional constraint according to the target inequality, to obtain the second functional constraint after relaxation processing; The expression of each line power flow constraint in the second function constraint after relaxation is:

(U _jt ) ² =U _it ′-2(R _ij P _ij.t +X _ij Q _ij.t )+[(R _ij ) ² +(X _ij ) ² ]·I _ij.t ′; Among them, P _it represents the active energy of node i at time t, Q _it represents the reactive energy of node i at time t, P _jt represents the active energy of node j at time t, Q _jt represents the reactive energy of node j at time t Electricity, U _it represents the voltage amplitude of node i at time t, U _jt represents the voltage amplitude of node j at time t, R _ij represents the resistance of the line between node i and node j, X _ij represents node i and node j The reactance of the line between nodes, G _ij represents the conductance of the line between node i and node j, B _ij represents the susceptance of the line between node i and node j;

Represents the voltage phase angle difference between node i and node j; u(j) represents the upstream node set connected to node j, d(j) represents the downstream node set connected to node j, and P _jl.t represents node j and The active power between node l at time t, Q _jl.t represents the reactive power between node j and node l at time t. 9 . The device for configuring incremental distributed power equipment according to claim 7 , wherein the constraint on the number of nodes connected to distributed power equipment is preset as: N _{i.min ≤N i ≤N i.max} ; Among them, N _i.min represents the minimum value of the number of distributed power supply devices connected to node i, N _i.max represents the maximum number of distributed power supply devices connected to the preset node i, and N _i represents the ith node connected to the number of distributed power equipment; The upper limit constraint on the penetration rate of distributed power equipment is expressed as:

Among them, β represents the penetration rate after the distributed power equipment is connected to the node, P _load represents the power load of the distributed power equipment at the preset node; n _i represents a variable, which takes a value of 0 or 1, and n _i =0 means In the node distribution network structure, the ith node is not connected to the distributed power equipment, and n _i = 1 means that the ith node is connected to the distributed power equipment; P _sg .DG represents the rated active power of each distributed power equipment; The total power capacity constraint of distributed power equipment is expressed as: P _{min.DG ≤P t.DG ≤P max.DG} ; Among them, P _min.DG represents the lower limit of the total power supply capacity of the distributed power equipment, P _max.DG represents the upper limit of the total power supply capacity of the distributed power equipment; P _t.DG represents the power supply provided by the distributed power equipment at time t Active power. 10. A computer-readable storage medium, characterized in that a detection program is stored on the computer-readable storage medium, and when the detection program is executed by a processor, the enhancement according to any one of claims 1-6 is realized. Steps of a configuration method of a distributed power supply device.

Images