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

Configuration method and device of incremental distributed power supply equipment and readable storage medium Download PDF

Info

Publication number
CN111952968B
CN111952968B CN202010818367.XA CN202010818367A CN111952968B CN 111952968 B CN111952968 B CN 111952968B CN 202010818367 A CN202010818367 A CN 202010818367A CN 111952968 B CN111952968 B CN 111952968B
Authority
CN
China
Prior art keywords
node
representing
power supply
distributed power
constraint
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010818367.XA
Other languages
Chinese (zh)
Other versions
CN111952968A (en
Inventor
张渊
李绮
严以臻
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Changzhou Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Changzhou Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Jiangsu Electric Power Co Ltd, Changzhou Power Supply Co of State Grid Jiangsu Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202010818367.XA priority Critical patent/CN111952968B/en
Publication of CN111952968A publication Critical patent/CN111952968A/en
Application granted granted Critical
Publication of CN111952968B publication Critical patent/CN111952968B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/04—Power grid distribution networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Power Engineering (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention provides a configuration method, a configuration device and a readable storage medium of incremental distributed power supply equipment, wherein the method comprises the following steps: constructing a distributed power incremental distribution network model at least comprising a node distribution network structure; constructing target functions and function constraints corresponding to distributed power supply nodes, incremental power distribution network nodes and power nodes in an incremental power distribution network model, multiplying the target functions by corresponding weights to obtain corresponding weight functions, and adding the weight functions corresponding to the distributed power supply nodes, the incremental power distribution network nodes and the power nodes to obtain an overall planning function; and performing second-order cone relaxation treatment on the function constraint of the incremental distribution network node, and solving an overall planning function through a solver based on the treated function constraint and the untreated function constraint to obtain the power supply capacity of the node to be accessed to the distributed power supply equipment. The invention improves the rationality of the power supply capacity of each node in the incremental power distribution network system.

Description

Configuration method and device of incremental distributed power supply equipment and readable storage medium
Technical Field
The invention relates to the technical field of power systems, in particular to a configuration method and device of incremental distributed power supply equipment and a readable storage medium.
Background
The incremental distribution network refers to a local area power grid of a voltage class of 110kV (kilovolt) and below and a voltage class of 220(330) kV and below, an industrial park (an economic development area) and the like in principle, and does not relate to the construction of a power transmission network of 220kV and above. The incremental power distribution network comprises a newly-built incremental power distribution network, a power distribution network capacity increase and expansion and a power distribution network stock outside the power distribution network enterprise stock. With the continuous development of the urbanization process and the development of the economic society in China, an incremental distribution network is indispensable, the planning and operation mode of the traditional distribution network needs to be changed urgently, and the incremental distribution network is combined with the latest distribution network operation technology, so that the aims of reducing cost, improving quality and efficiency, saving energy and reducing emission are favorably fulfilled. In addition, Distributed Generation (DG) devices are incorporated into an incremental power distribution network system, and the direction of power flow changes from one-way to two-way, so that the uncertainty factor of operation increases, and certain safety and stability challenges are brought to the power selling company while the cost is reduced and the benefit is increased. Therefore, it is urgently needed to provide an incremental distribution network optimization planning method on the premise of ensuring the safe and stable operation of an incremental distribution network system so as to improve the rationality of the power supply capacity of each node in the incremental distribution network system.
Disclosure of Invention
Based on the above current situation, it is urgently needed to provide a method for improving the rationality of the power supply capacity of each node in an incremental power distribution network system on the premise of ensuring the safe and stable operation of the incremental power distribution network system.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
a method of configuring an incremental distributed power supply device, the method comprising the steps of:
s100, constructing a distributed power incremental distribution network model in an incremental distribution network system, wherein the incremental distribution network model at least comprises a node distribution network structure;
s200, obtaining first target data related to distributed power supply nodes, second target data related to incremental power distribution network nodes and third target data related to power nodes in the incremental power distribution network model from a preset database of the incremental power distribution network system, constructing a first target function corresponding to the distributed power supply nodes and a first function constraint constraining the first target function according to the first target data, constructing a second target function corresponding to the incremental power distribution network nodes and a second function constraint constraining the second target function according to the second target data, and constructing a third target function corresponding to the power nodes and a third function constraint constraining the third target function according to the third target data, wherein the first function constraint comprises constraint, access number of distributed power supply equipment of preset nodes in a distribution network structure, The second function constraint comprises power flow constraint, voltage constraint and power constraint of each line in the incremental distribution network model, and the third function constraint comprises power load constraint transferred at a target moment and power load constraint transferred at a target moment;
s300, multiplying the first objective function by a first weight corresponding to the first objective function to obtain a first weight function, multiplying the second objective function by a second weight corresponding to the second objective function to obtain a second weight function, and multiplying the third objective function by a third weight corresponding to the third objective function to obtain a 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, performing second-order cone relaxation processing on the second function constraint to obtain a second function constraint after relaxation processing, solving 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 to obtain the power supply capacity of the node to be connected to the distributed power supply equipment in the node distribution network structure, and configuring the distributed power supply equipment corresponding to the node according to the power supply capacity.
Preferably, in the step S500, during the second-order cone relaxation process, the auxiliary variable U is definedi.t′＝(Ui.t)2And Iij.t′＝(Iij.t)2According to saidThe auxiliary variables yield the following target equation:
Figure GDA0003452711880000021
Iij.trepresenting the line current, U, between node i and node j in a node-distribution network architecturei.tRepresenting the voltage amplitude, P, of node i at time tij.tRepresenting the active electric quantity, Q, transmitted by branch ij at time tij.tThe reactive power transmitted by the branch ij at the moment t is represented;
and carrying out relaxation treatment on the target equation to obtain the following relaxation inequality:
Figure GDA0003452711880000031
reconstructing the relaxation inequality through a second order cone to obtain the following target inequality:
Figure GDA0003452711880000032
performing second-order cone relaxation processing on each line power flow constraint in the second function constraint according to the target inequality to obtain a second function constraint after relaxation processing;
the expression of each line power flow constraint in the relaxed second function constraint is as follows:
Figure GDA0003452711880000033
Figure GDA0003452711880000034
Figure GDA0003452711880000035
Figure GDA0003452711880000036
(Uj.t)2＝Ui.t′-2(RijPij.t+XijQij.t)+[(Rij)2+(Xij)2]·Iij.t′；
wherein, Pi.tRepresents the active electric quantity, Q, of the node i at the moment ti.tRepresenting the reactive power of the node i at time t, Pj.tRepresents the active electric quantity, Q, of the node j at the moment tj.tRepresenting the reactive electric quantity, U, of the node j at time ti.tRepresenting the voltage amplitude, U, of node i at time tj.tRepresenting the voltage magnitude, R, of node j at time tijRepresenting the resistance of the line between node i and node j, XijRepresenting the reactance of the line between node i and node j, GijRepresenting the conductance of the line between node i and node j, BijRepresenting susceptance of the line between node i and node j;
Figure GDA0003452711880000037
represents the voltage phase angle difference between node i and node j; u (j) represents an upstream set of nodes connected to node j, d (j) represents a downstream set of nodes connected to node j, Pjl.tRepresents the active electric quantity between the node j and the node l at the moment t, Qjl.tAnd the reactive power between the node j and the node l at the time t is represented.
Preferably, in the step S200,
the constraint of the number of the preset nodes accessing the distributed power supply devices is represented as follows: n is a radical ofi.min≤Ni≤Ni.max；
Wherein N isi.minMinimum value, N, representing the number of devices accessing the distributed power supply at node ii.maxMaximum value N representing the number of devices accessing the distributed power supply at preset node iiThe number of the distributed power supply devices accessed by the ith node is represented;
the upper limit constraint on the permeability of the distributed power equipment is expressed as:
Figure GDA0003452711880000041
wherein, beta represents the permeability of the distributed power supply equipment after being connected to the node, PloadRepresenting the electric load quantity of the distributed power supply equipment at a preset node; n isiRepresents a variable with a value of 0 or 1, niWhen the node is 0, the ith node in the node distribution network structure is not connected with the distributed power supply device, ni1 represents that the ith node is accessed into the distributed power supply equipment; psg.DGThe rated active power of each distributed power supply device is represented;
the total power capacity constraint for a distributed power plant is expressed as: pmin.DG≤Pt.DG≤Pmax.DG；
Wherein, Pmin.DGLower limit, P, representing total power supply capacity of distributed power supply equipmentmax.DGAn upper limit value representing a total power capacity of the distributed power equipment; pt.DGRepresenting the active power provided by the distributed power equipment at time t.
Preferably, in the step S200,
the expression of each line tide constraint is as follows:
Figure GDA0003452711880000042
Figure GDA0003452711880000043
Figure GDA0003452711880000044
Figure GDA0003452711880000045
Figure GDA0003452711880000046
wherein, Pi.tRepresents the active electric quantity, Q, of the node i at the moment ti.tRepresenting the reactive power of the node i at time t, Pj.tThe expression is the active electric quantity, Q, of the node j at the moment tj.tThe expression is the reactive power of the node j at the time t, Pij.tRepresenting the active electric quantity, Q, transmitted by branch ij at time tij.tRepresenting the reactive electric quantity, U, transmitted by branch ij at time ti.tRepresenting the voltage amplitude, U, of node i at time tj.tRepresenting the voltage magnitude, R, of node j at time tijRepresenting the resistance of the line between node i and node j, XijRepresenting the reactance of the line between node i and node j, GijRepresenting the conductance of the line between node i and node j, BijRepresenting susceptance of the line between node i and node j;
Figure GDA0003452711880000051
represents the voltage phase angle difference between node i and node j; u (j) represents an upstream set of nodes connected to node j, d (j) represents a downstream set of nodes connected to node j, Pjl.tRepresents the active electric quantity between the node j and the node l at the moment t, Qjl.tRepresenting the reactive power between the node j and the node l at the time t;
the voltage constraints are: u shapei.min≤Ui.t≤Ui.maxWherein, Ui.minRepresents the minimum value, U, of the voltage amplitude of node ii.maxRepresents the maximum value of the voltage amplitude of the node i;
the power constraint is: pij.t≤Pij.maxWherein P isij.tRepresenting the power flowing on the line between node i and node j at time t, Pij.maxWhich represents the maximum value of the power flowing on the line between node i and node j.
Preferably, in the step S200,
the target time is t time, and the expressions of the electric power load quantity constraint transferred at the target time and the electric power load quantity constraint transferred at the target time are as follows:
Figure GDA0003452711880000052
wherein ρminPt.loadRepresents the minimum value of the power coefficient of the load transfer at time t, rhomaxPt.loadThe maximum value of the power coefficient of the load transfer at the moment t is shown; sigmaminPt.loadRepresenting the minimum value of the power coefficient, σ, of the load transfer at time tmaxPt.loadRepresenting the maximum value of the power coefficient of load transfer at the moment t; pt.outRepresenting the electric load amount transferred out at the time t; pt.inIndicating the amount of electrical load transferred at time t.
Preferably, in step S500, the step of obtaining the power supply capacity of the node to be connected to the distributed power supply device in the node distribution network structure by solving the overall planning function through a solver based on the first function constraint, the third function constraint and the relaxed second function constraint includes:
based on the first function constraint, the third function constraint and the relaxed second function constraint, solving the overall planning function through a Cplex solver to obtain the equipment number of the distributed power equipment nodes to be accessed in the node distribution network structure, wherein the equipment number corresponds to the distributed power equipment;
the method comprises the steps of obtaining the power capacity of one piece of distributed power supply equipment, and multiplying the power capacity of the distributed power supply equipment by the number of the equipment to obtain the power capacity of nodes to be connected into the distributed power supply equipment, wherein the power capacity of each piece of distributed power supply equipment is the same.
The invention also provides a configuration device of the incremental distributed power supply equipment, which comprises the following steps:
the incremental distribution network model comprises a building module and a control module, wherein the building module is used for building a distributed power incremental distribution network model in an incremental distribution network system, and the incremental distribution network model at least comprises a node distribution network structure;
the acquisition module is used for acquiring first target data related to distributed power source nodes, second target data related to incremental power distribution network nodes and third target data related to power nodes in the incremental power distribution network model from a preset database of the incremental power distribution network system;
the building module is further configured to build a first objective function corresponding to the distributed power supply node and a first function constraint constraining the first objective function according to the first target data, build a second objective function corresponding to the incremental power distribution network node and a second function constraint constraining the second objective function according to the second target data, and build a third objective function corresponding to the power node and a third function constraint constraining the third objective function according to the third target data, where the first function constraint includes a constraint on the number of distributed power supply devices accessed by preset nodes in a node distribution network structure, a distributed power supply device permeability upper limit constraint and a distributed power supply device total power supply capacity constraint, and the second function constraint includes a power flow constraint, a voltage constraint and a power constraint of each line in the incremental power distribution network model, the third function constraint comprises a power load amount constraint transferred at a target moment and a power load amount constraint transferred at the target moment;
a calculating module, configured to multiply the first objective function by a first weight corresponding to the first objective function to obtain a first weight function, multiply the second objective function by a second weight corresponding to the second objective function to obtain a second weight function, and multiply the third objective function by a third weight corresponding to the third objective function to obtain a third weight function; adding the first weight function, the second weight function and the third weight function to obtain an overall planning function;
the relaxation processing module is used for performing second-order cone relaxation processing on the second function constraint to obtain a second function constraint after the relaxation processing;
the solving module is used for solving the overall planning function through a solver based on the first function constraint, the third function constraint and the relaxed second function constraint to obtain the power supply capacity of the node to be connected to the distributed power supply equipment in the node distribution network structure;
and the configuration module is used for configuring the distributed power supply equipment corresponding to the node according to the power supply capacity.
The invention also provides a computer-readable storage medium, on which a detection program is stored, which, when executed by a processor, implements the steps of the method for configuring an incremental distributed power supply device as described.
By the technical scheme, mutual coordination of the main bodies is realized, the feasibility of respective decision is improved, the value of the overall planning function is optimized, the power capacity configuration of each node connected with the distributed power supply equipment in the incremental power distribution network system is optimized, the rationality of the power capacity of each node in the incremental power distribution network system is improved, namely the rationality of the configuration of the distributed power supply equipment corresponding to the node in the node distribution network structure is improved, and the safe and stable operation of the incremental power distribution network system is ensured.
Other advantages of the present invention will be described in the detailed description, and those skilled in the art will understand the technical features and technical solutions presented in the description.
Drawings
Preferred embodiments according to the present invention will be described below with reference to the accompanying drawings. In the figure:
fig. 1 is a flow chart of a configuration method of an incremental distributed power supply device according to an embodiment of the present invention.
Fig. 2 is a schematic diagram of a distribution network structure of IEEE33 nodes according to an embodiment of the present invention;
fig. 3 is a schematic structural diagram of a configuration device of the incremental distributed power supply apparatus of the present invention.
Detailed Description
The preferred embodiments of the present invention will now be described in detail with reference to the accompanying drawings.
It should be noted that step numbers (letter or number numbers) are used to refer to some specific method steps in the present invention only for the purpose of convenience and brevity of description, and the order of the method steps is not limited by letters or numbers in any way. It will be clear to a person skilled in the art that the order of the steps of the method in question, as determined by the technology itself, should not be unduly limited by the presence of step numbers.
Fig. 1 is a flow chart of a configuration method of an incremental distributed power supply device according to an embodiment of the present invention.
S100, constructing a distributed power incremental distribution network model in an incremental distribution network system, wherein the incremental distribution network model at least comprises a node distribution network structure.
And constructing a distributed power supply incremental distribution network model in the incremental distribution network system, wherein the incremental distribution network model is a DG incremental distribution network model. The incremental power distribution network system is a distributed system, the incremental power distribution network model at least comprises a node distribution network structure, and the node distribution network structure can be an IEEE (Institute of Electrical and Electronics Engineers) 33 node distribution network structure, an IEEE14 node distribution network structure or an IEEE30 node distribution network structure and the like. Specifically, referring to fig. 2, fig. 2 is a schematic diagram of a distribution network structure of IEEE33 nodes in the embodiment of the present invention. As can be seen from fig. 2, in the node distribution network structure, the node 5, the node 8, the node 14, the node 16, the node 18, and the node 31 are access nodes of the distributed power supply device. The incremental power distribution network model also includes branch data information and node data information. The node distribution network structure describes the connection relation among nodes in the incremental distribution network model; the branch data information specifically refers to the impedance of a network line in the incremental power distribution network model; the node data information includes the load size of each node of the incremental distribution network, the total power capacity of the distributed power supply device accessed by each node, and the node where each distributed power supply device is located. The incremental power distribution network system comprises three benefit agents which are respectively DG operation users, power distribution network investment users and power users, and the three benefit agents are stored in the incremental power distribution network system in a node mode, namely distributed power source nodes corresponding to the DG operation users, incremental power distribution network nodes corresponding to the power distribution network investment users and power nodes corresponding to the power users are stored in the incremental power distribution network system. For convenience of description, distributed power nodes, incremental distribution network nodes and power nodes are used subsequently to describe various benefit subjects.
S200, obtaining first target data related to distributed power supply nodes, second target data related to incremental power distribution network nodes and third target data related to power nodes in the incremental power distribution network model from a preset database of the incremental power distribution network system, constructing a first target function corresponding to the distributed power supply nodes and a first function constraint constraining the first target function according to the first target data, constructing a second target function corresponding to the incremental power distribution network nodes and a second function constraint constraining the second target function according to the second target data, and constructing a third target function corresponding to the power nodes and a third function constraint constraining the third target function according to the third target data, wherein the first function constraint comprises constraint, access number of distributed power supply equipment of preset nodes in a distribution network structure, The second function constraint comprises power flow constraint, voltage constraint and power constraint of each line in the incremental distribution network model, and the third function constraint comprises power load constraint transferred at a target moment and power load constraint transferred at a target moment.
In the incremental power distribution network system, a database is preset and stores and constructs target data of corresponding target functions and function constraints of distributed power nodes, incremental power distribution network nodes and power nodes. When a corresponding target function and function constraint need to be constructed, acquiring first target data, second target data and third target data from a database, constructing a first target function corresponding to a distributed power supply node and a first function constraint constraining the first target function by the first target data, constructing a second target function corresponding to an incremental power distribution network node and a second function constraint constraining the second target function according to the second target data, and constructing a third target function corresponding to a power node and a third function constraint constraining the third target function according to the third target data. The first function constraint comprises the constraint of the number of preset nodes connected into the distributed power supply equipment in the node distribution network structure, the upper limit constraint of the permeability of the distributed power supply equipment and the total power supply capacity constraint of the distributed power supply equipment, the second function constraint comprises the power flow constraint, the voltage constraint and the power constraint of each line in the incremental distribution network model, and the third function constraint comprises the power load constraint transferred at the target moment and the power load constraint transferred at the target moment. It should be noted that in the process of constructing the objective function and the function constraint, the objective function corresponding to the distributed power source nodes, the incremental power distribution network nodes and the power nodes is set from the individual rationality of the distributed power source nodes, the incremental power distribution network nodes and the power nodes, and because the objective deviation is different when the distributed power source nodes, the incremental power distribution network nodes and the power nodes participate in the incremental power distribution network model and correspond to the incremental power distribution network system grid planning decision, the objective function corresponding to the distributed power source nodes, the incremental power distribution network nodes and the power nodes needs to be respectively established according to the investment and expected profit of each benefit principal.
Specifically, in the step S200, the first objective function expression is formula (1):
formula (1) maxWDG(ni,Ni)＝WS.DG-WI.DG-WOM.DG；
Wherein, W in the formula (1)S.DGIs expressed by formula (2):
Figure GDA0003452711880000091
w in formula (1)I.DGIs formula (3):
Figure GDA0003452711880000092
w in formula (1)OM.DGIs formula (4):
Figure GDA0003452711880000101
nirepresents a variable with a value of 0 or 1, niWhen the node is 0, the ith node in the node distribution network structure is not connected with the distributed power supply device, ni1 represents that the ith node is accessed into the distributed power supply equipment; n is a radical ofiThe number of the distributed power supply devices accessed by the ith node is represented; wS.DGRepresents the electricity sales income within the preset years of all distributed power supply equipment in the node distribution network structure, WI.DGRepresents the early construction cost, W, of all distributed power supply equipment in the node distribution network structure before useOM.DGRepresents the operation and maintenance cost, U, of all distributed power supply equipment within the preset agetRepresents the set of all the time in one day, which is 24 hours; q. q.sesThe unit electricity selling price of the distributed power supply equipment is represented, and the unit electricity selling price of each distributed power supply equipment is the same; pt.DGRepresenting the active power provided by the distributed power supply equipment at the moment t; q. q.ssgRepresents a distributed power supply equipment construction cost, specifically, qsgThe construction investment of the distributed power supply equipment corresponding to the unit power supply capacity is represented; u shapeiRepresenting that a node set to be accessed to distributed power equipment is preset in a node distribution network structure; psg.DGThe rated active power of each distributed power supply device is represented, and the rated active power of each distributed power supply device is the same; y represents the life cycle of each distributed power supply device, and the life cycle of each distributed power supply device is the same; a represents the discount rate; q. q.somAnd the operation and maintenance cost of unit power supply of the distributed power supply equipment is represented. The size of the preset year can be set according to specific needs, for example, the preset year can be set to 5 years, 10 years or 12 years. Upfront construction costs of the distributed power equipment prior to use include, but are not limited to, distributed power equipment design costs, purchase costs, and installation costs. It should be noted that the first objective function is an optimization objective function of the distributed power node, and the early-stage construction investment of the distributed power equipment is reduced by the first objective functionThe cost of the distributed power nodes is increased, and the cost of the operation and maintenance expenditure is increased, so that the net income of the distributed power nodes is maximized.
The constraint of the number of the preset nodes accessing the distributed power supply equipment is expressed as a formula (5): n is a radical ofi.min≤Ni≤Ni.max。
Wherein N isi.minMinimum value, N, representing the number of devices accessing the distributed power supply at node ii.maxMaximum value, N, representing the number of devices accessing the distributed power supply at preset IEEE33 node iiIndicating the number of distributed power supply devices accessed by the ith node.
The upper limit constraint on the permeability of the distributed power equipment is expressed by the formula (6):
Figure GDA0003452711880000102
wherein, beta represents the permeability of the distributed power supply equipment after being connected to the node, and the permeability is the proportion of the capacity of the distributed power supply to the rated capacity of the power distribution network. PloadRepresenting the electric load quantity of the distributed power supply equipment at a preset node; n isiRepresents a variable with a value of 0 or 1, niWhen the node is 0, the ith node in the node distribution network structure is not connected with the distributed power supply device, ni1 represents that the ith node is accessed into the distributed power supply equipment; psg.DGAnd the rated active power of each distributed power supply device is shown.
The total power capacity constraint of the distributed power equipment is expressed by the formula (7): pmin.DG≤Pt.DG≤Pmax.DG。
Wherein, Pmin.DGLower limit, P, representing total power supply capacity of distributed power supply equipmentmax.DGAn upper limit value representing a total power capacity of the distributed power equipment; pt.DGRepresenting the active power provided by the distributed power equipment at time t. Each data in the first objective function and the first function constraint expression is first objective data, the first objective data can be obtained by the incremental distribution network system according to the historical data of the distributed power nodes and is pre-stored in a pre-set database, and when needed,and directly obtaining the data from a local database.
Specifically, in the step S200, the expression of the second objective function is formula (8):
equation (8) maxWDN(pi)＝WS.DN-WL.DN-WE.DN-WB1.DN-WB2.DN-WB3.DG；
Wherein W in the formula (8)S.DGIs expressed by formula (9):
formula (9):
Figure GDA0003452711880000111
w in formula (8)L.DNIs expressed as formula (10):
Figure GDA0003452711880000112
w in formula (8)E.DNIs expressed by formula (11):
formula (11)
Figure GDA0003452711880000113
W in formula (8)B1.DNIs formula (12):
formula (12)
Figure GDA0003452711880000114
W in formula (8)B2.DNIs formula (13):
Figure GDA0003452711880000115
w in the formula (8)B3.DGIs expressed as formula (14):
Figure GDA0003452711880000116
wherein p isiVariable of an objective function, W, representing a node of an incremental distribution networkS.DNRepresenting the electricity selling and receiving of incremental distribution network nodes within a preset yearIn, WL.DNRepresenting the line loss cost, W, over a predetermined period of timeE.DNIndicating a fault repair cost, W, over a predetermined period of timeB1.DNRepresents the electricity purchase cost from the higher level within a preset age, WB2.DNRepresenting the purchase of electricity from distributed power nodes within a predetermined age, WB3.DGRepresents the loss, U, caused by the output fluctuation of the distributed power supply equipmenttRepresenting the set of all times of the day, fesRepresenting a unit price for electricity sold to the power node; pt.loadRepresenting the initial electric power load quantity at the time t; pt.outRepresenting the electric load amount transferred out at the time t; pt.inRepresenting the amount of electric power load transferred at time t; pt.isRepresenting the amount of electric load which can be interrupted at the moment t; pt.lossRepresenting the active loss at the moment t; EENStThe expected value of less power supply energy to the power node at the moment t is shown; u shapebRepresenting a set of all lines of a power grid in the incremental power distribution grid model; lambda [ alpha ]bA probability value representing the fault occurrence of the b-th line; u shapenRepresenting that a node set to be accessed to distributed power equipment is preset in a node distribution network structure; pn.t.loadRepresenting the initial electric load quantity of the node n at the time t; f. ofeb1Represents a unit purchase price from an upper-level power grid; f. ofeb2Represents a unit purchase price from the distributed power node; ppRepresenting punishment cost required to be paid due to random fluctuation within a preset age; rhosRepresenting the probability of power loss caused by the output fluctuation of the distributed power supply equipment within a preset period; delta QsRepresenting the random fluctuation amount when the distributed power supply equipment is fluctuated. It should be noted that the preset age of the first objective function and the second objective function is the same. Establishing an incremental distribution network node optimization objective function, namely establishing a second objective function, expecting to reduce the cost of transmission line loss, early-stage construction, upper-stage electricity purchasing and the like, improving the electricity selling income and maximizing the benefit of the incremental distribution network node as much as possible; furthermore, the uncertainty of the distributed power supply equipment access node more affects the operation safety of the incremental power distribution network system, so that the economy of the incremental power distribution network node optimization planning scheme is reduced, and the reduction of the distributed power supply equipment output fluctuation is expected to be reducedAnd (4) loss. In this embodiment, the loss due to the output fluctuation of the distributed power supply apparatus can be represented by the cost corresponding to the loss. The loss caused by the output fluctuation of the distributed power supply equipment is the loss caused by the unstable power supply of the distributed power supply equipment.
The expression of each line power flow constraint is formula (15) -formula (19):
formula (15)
Figure GDA0003452711880000121
Formula (16)
Figure GDA0003452711880000122
Formula (17)
Figure GDA0003452711880000123
Formula (18)
Figure GDA0003452711880000124
Formula (19)
Figure GDA0003452711880000125
Wherein, Pi.tRepresents the active electric quantity, Q, of the node i at the moment ti.tRepresenting the reactive power of the node i at time t, Pj.tThe expression is the active electric quantity, Q, of the node j at the moment tj.tThe expression is the reactive power of the node j at the time t, Pij.tRepresenting the active electric quantity, Q, transmitted by branch ij at time tij.tRepresenting the reactive electric quantity, U, transmitted by branch ij at time ti.tRepresenting the voltage amplitude, U, of node i at time tj.tRepresenting the voltage magnitude, R, of node j at time tijRepresenting the resistance of the line between node i and node j, XijRepresenting the reactance of the line between node i and node j, GijRepresenting the conductance of the line between node i and node j, BijRepresenting susceptance of the line between node i and node j;
Figure GDA0003452711880000131
represents the voltage phase angle difference between node i and node j; u (j) represents an upstream set of nodes connected to node j, d (j) represents a downstream set of nodes connected to node j, Pjl.tRepresents the active electric quantity between the node j and the node l at the moment t, Qjl.tAnd the reactive power between the node j and the node l at the moment t is represented.
The expression of the voltage constraint is formula (20): u shapei.min≤Ui.t≤Ui.max。
Wherein, Ui.minRepresenting the minimum value, U, of the voltage amplitude of node ii.maxRepresenting the maximum value of the voltage amplitude at node i.
The expression of the power constraint is formula (21): pij.t≤Pij.max。
Wherein, Pij.tIndicating the power flowing on the line between node i and node j at time t, Pij.maxWhich represents the maximum value of the power flowing on the line between node i and node j.
It should be noted that, in this embodiment, after the distributed power supply device is accessed to the node, problems of a trend change, an active management mode change, an increase in uncertain factors, and the like in the incremental power distribution network system are considered, and by adding a loss caused by an output fluctuation to the second objective function of the incremental power distribution network node, an adverse effect of the uncertain factors generated after the distributed power supply device is accessed to the node, that is, after the distributed power supply device enters the incremental power distribution network system, is reduced, so that interference of the distributed power supply device on the overall planning function is reduced, and the accuracy of the overall planning model obtained through the method is improved.
Each data in the second objective function and the second function constraint expression is second objective data, the second objective data can be obtained by the incremental distribution network system according to historical data of the incremental distribution network nodes, and is stored in a preset database in advance, and when needed, the second objective data can be directly obtained from the local database.
Specifically, in the step S200, the expression of the third objective function is formula (22):
formula (22)
Figure GDA0003452711880000132
Wherein, γebRepresenting the unit price, P, of electricity purchased by a power node to an incremental distribution network nodet.outRepresenting the electric load amount transferred out at the time t; pt.inRepresenting the amount of electric power load transferred at time t; pt.isRepresenting the amount of electric load which can be interrupted at the moment t; pt.lossRepresenting the active loss, U, at time ttRepresenting the set of all times of the day. The transferred electric power load amount is a multi-purpose electric power load amount of the electric power node relative to a set reference electric power load amount; the transferred power load amount is the power load amount which is less than the set reference power load amount of the power node; the amount of interrupted power load is the amount of power load that the power node is not using at the time of the power interruption. It is to be noted that, the power node objective function is established, and it is desired to reduce the electricity fee expenditure by adjusting the electricity usage plan, the demand-side response method to be considered is a price-based demand-side response (DSR) based on the time-of-use electricity price, and the power nodes participating in the DSR move out of the load during the peak period of the electricity price and move in the load during the valley period of the electricity price based on the time-of-use electricity price information.
The target time is time t, and the expressions of the electric power load amount constraint transferred at the target time and the electric power load amount constraint transferred at the target time are formula (23):
formula (23)
Figure GDA0003452711880000141
Where ρ isminPt.loadRepresents the minimum value of the power coefficient of the load transfer at time t, rhomaxPt.loadThe maximum value of the power coefficient of the load transfer at the moment t is shown; sigmaminPt.loadRepresenting the minimum value of the power coefficient, σ, of the load transfer at time tmaxPt.loadRepresenting the maximum value of the power coefficient of load transfer at the moment t; pt.outRepresenting the electric load amount transferred out at the time t; pt.inIndicating a shift at time tThe amount of electrical load. And each data in the third objective function and the third function constraint expression is 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 is stored in a preset database in advance, and when needed, the third objective data can be directly obtained from the local database.
In this embodiment, there is a "max" in the first, second, and third objective functions to ensure that the values obtained through the first, second, and third objective functions are all maximized, even though the benefits of the distributed power nodes, incremental distribution nodes, and power nodes are maximized.
And S300, multiplying the first objective function by a first weight corresponding to the first objective function to obtain a first weight function, multiplying the second objective function by a second weight corresponding to the second objective function to obtain a second weight function, and multiplying the third objective function by a third weight corresponding to the third objective function to obtain a third weight function.
And S400, adding the first weight function, the second weight function and the third weight function to obtain an overall planning function.
It can be understood that the benefit appeal of the three market subjects, namely the distributed power supply node, the incremental distribution network node and the power node, is different, the target bias is different when the distributed power supply node, the incremental distribution network node and the power node participate in planning, and the respective benefit is expected to be maximized. The maximization of the personal interests of each subject may cause a disadvantage that the overall interests are far from the optimal, therefore, in this embodiment, a weight is set for each objective function, the weight size corresponding to each objective function is determined according to the weight and importance of the distributed power nodes, the incremental distribution network nodes and the power nodes in the incremental distribution network system, and the size of each weight is not limited in this embodiment. For convenience of description, the first weight is denoted as r1And the second weight is denoted as r2And the third weight is denoted as r3Wherein r is1+r2+r3When the overall planning function is denoted as W, the overall planning function can be expressed as formula (24):
formula (24) W ═ γ1WDG+γ2WDN+γ3WUS。
It should be noted that there is a "max" in each of the first, second and third objective functions, and therefore in the specific calculation in equation (24), WDG、WDNAnd WUSAlso the maximum value, just in equation (24), "max" preceding the first, second and third objective functions is omitted for convenience of description. As can be seen from the formula (24), r1WDGIs a first weight function, r2WDNIs a second weight function, r3WUSIs a third weighting function.
S500, performing second-order cone relaxation processing on the second function constraint to obtain a second function constraint after relaxation processing, solving 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 to obtain the power capacity of a node to be connected to the distributed power supply equipment in the node distribution network structure, and configuring the distributed power supply equipment corresponding to the node according to the power capacity.
And performing second-order cone relaxation processing on the second function constraint to obtain a second function constraint after the relaxation processing. Specifically, second-order cone relaxation processing is performed on each line power flow constraint in the second function constraint to obtain a second function constraint after relaxation processing, then an overall planning function is solved through a solver under the conditions of the first function constraint, the third function constraint and the second function constraint after relaxation processing to obtain the power capacity of nodes to be connected with distributed power equipment in the node distribution network structure, in each node of the node distribution network structure, each node position can be connected with at least one distributed power equipment, but not all nodes in the node distribution network structure need to be connected with the distributed power equipment. Specifically, the solver is a Cplex solver. In this embodiment, in order to improve the solving speed of solving the overall planning function, all nodes in the node distribution network structure are not taken as nodes of the distributed power supply device to be accessed, but only a few of the nodes are determined as nodes of the distributed power supply device to be accessed, when the power capacity of the nodes of the distributed power supply device to be accessed is 0, it indicates that the nodes of the distributed power supply device to be accessed do not need to be accessed to the distributed power supply device, and only when the power capacity of the nodes of the distributed power supply device to be accessed is not 0, it indicates that the nodes of the distributed power supply device to be accessed need to be accessed to the distributed power supply device.
Specifically, in the step S500, during the second-order cone relaxation process, the auxiliary variable U is definedi.t′＝(Ui.t)2And Iij.t′＝(Iij.t)2And obtaining the following target equation according to the auxiliary variable, wherein the expression of the target equation is the formula (25):
equation (25)
Figure GDA0003452711880000161
Wherein, Iij.tRepresenting the line current, U, between node i and node j in a node-distribution network architecturei.tRepresenting the voltage amplitude, P, of node i at time tij.tRepresenting the active electric quantity, Q, transmitted by branch ij at time tij.tRepresenting the reactive power transmitted by the branch ij at the moment t;
and (3) carrying out relaxation treatment on the target equation to obtain the following relaxation inequality, wherein the expression of the relaxation inequality is a formula (26):
formula (26)
Figure GDA0003452711880000162
Reconstructing a relaxation inequality through a second-order cone to obtain the following target inequality, wherein the expression of the target inequality is formula (27):
formula (27)
Figure GDA0003452711880000163
Performing second-order cone relaxation processing on each line power flow constraint in the second function constraint according to the target inequality to obtain a second function constraint after relaxation processing;
the expression of each line power flow constraint in the second function constraint after the relaxation processing is formula (28) -formula (32):
formula (28)
Figure GDA0003452711880000164
Formula (29)
Figure GDA0003452711880000165
Formula (30)
Figure GDA0003452711880000166
Formula (31)
Figure GDA0003452711880000167
Equation (32) (U)j.t)2＝Ui.t′-2(RijPij.t+XijQij.t)+[(Rij)2+(Xij)2]·Iij.t′。
Wherein, Pi.tRepresents the active electric quantity, Q, of the node i at the moment ti.tRepresenting the reactive power of the node i at time t, Pj.tRepresents the active electric quantity, Q, of the node j at the moment tj.tRepresenting the reactive electric quantity, U, of the node j at time ti.tRepresenting the voltage amplitude, U, of node i at time tj.tRepresenting the voltage magnitude, R, of node j at time tijRepresenting the resistance of the line between node i and node j, XijRepresenting the reactance of the line between node i and node j, GijRepresenting the conductance of the line between node i and node j, BijRepresenting susceptance of the line between node i and node j;
Figure GDA0003452711880000171
represents the voltage phase angle difference between node i and node j; u (j) represents an upstream set of nodes connected to node j, d (j) represents a downstream set of nodes connected to node j, Pjl.tRepresents the active electric quantity between the node j and the node l at the moment t, Qjl.tAnd the reactive power between the node j and the node l at the moment t is represented.
In the second function constraint after the relaxation process, equation (28) is the same as equation (15), and equation (29) is the same as equation (16), i.e., equation (15) and equation (16) are not processed during the relaxation process.
When solving the overall planning function using the solver, it is required to obtain the maximum value of the overall planning function, that is, maxW ═ r1WDG+r2WDN+r3WUS. Solving the overall planning function through a solver to obtain constraint conditions of power supply capacity of the nodes of the distributed power supply equipment to be accessed in the node distribution network structure, wherein the constraint conditions can be expressed as follows:
Figure GDA0003452711880000172
it should be noted that, after determining the power supply capacity of each node of the distributed power supply devices to be accessed, the power supply capacity of each distributed power supply device is obtained, the number of the distributed power supply devices accessed by the node can be determined according to the power supply capacity of the node of the distributed power supply devices to be accessed and the power supply capacity of each distributed power supply device, and then the corresponding distributed power supply device is configured to the node, that is, the distributed power supply device corresponding to the node is configured according to the power supply capacity. In this embodiment, the sum of the power capacities of the distributed power supply devices of the access node is equal to or greater than the power capacity of the node to be accessed. The power capacities of the distributed power devices connected to the same node may be the same or different.
For convenience of understanding, for example, as an IEEE33 node power distribution network example is adopted, the incremental power distribution network system includes 33 branches, the simulation year is 10 years, that is, the preset year is 10 years, and the nodes to be connected to the distributed power supply device are node 5, node 8, node 14, node 16, node 18, and node 31, respectively. By solving the overall planning function, it can be determined which nodes of the nodes 5, 8, 14, 16, 18 and 31 the distributed power supply device is to be accessed to, in the nodes 5, 8, 14, 16, 18 and 31, not all the nodes are necessarily to be accessed to the distributed power supply device, and the number of the distributed power supply devices accessed to each node may be the same or different.
The highest market overall income within the simulation year obtained by calculation of a solver is 106.615 ten thousand yuan, namely the maximum value of the overall planning function is 106.615 ten thousand yuan, and the node positions and the power supply capacity of the distributed power supply equipment access when the overall planning function is equal to the maximum value are shown in the following table:
Figure GDA0003452711880000181
the values of the distributed power nodes, incremental distribution network nodes and power nodes (individual subject gains) and the maximum value of the overall planning function (overall benefit) are shown in the following table:
Figure GDA0003452711880000182
Figure GDA0003452711880000191
further, in order to improve the accuracy of the calculation of the maximum value of the overall planning function and the accuracy of the profit of the power node, the unit electricity price of different time periods is considered in the process of calculating the maximum value of the overall planning function, for example, the electricity purchase price is 800 (yuan/MW · h) during the peak period of electricity utilization and 8:00-11:00&18:00-23:00 each day; in the electricity consumption valley period, the electricity purchasing price is 400 (yuan/MW & h) in the period of 6:00-8:00&11:00-18:00 every day; in the flat period of electricity utilization, the electricity purchasing price is 600 (yuan/MW & h) in the 23:00-6:00 period of each day.
The embodiment realizes the mutual coordination of all the main bodies, improves the feasibility of respective decision, optimizes the value of the whole planning function, optimizes the power capacity configuration of all the nodes accessed with the distributed power supply equipment in the incremental power distribution network system, improves the rationality of the power capacity of all the nodes in the incremental power distribution network system, namely improves the rationality of the configuration of the distributed power supply equipment corresponding to the nodes in the node distribution network structure, and ensures the safe and stable operation of the incremental power distribution network system.
Further, another embodiment of the method for configuring an incremental distributed power supply device of the present invention is provided.
The difference between the another embodiment of the method for configuring incremental distributed power supply devices and the above embodiment of the method for configuring incremental distributed power supply devices is that, in step S500, the step of obtaining the power supply capacity of the node to be connected to the distributed power supply device in the node distribution network structure by solving the overall planning function through a solver based on the first function constraint, the third function constraint and the relaxed second function constraint includes:
step a, solving the overall planning function through a Cplex solver based on the first function constraint, the third function constraint and the relaxed second function constraint to obtain the number of equipment to be connected to the distributed power equipment nodes in the node distribution network structure, wherein the equipment corresponds to the distributed power equipment nodes;
and b, acquiring the power capacity of one distributed power supply device, and multiplying the power capacity of the distributed power supply device by the number of the devices to obtain the power capacity of the nodes to be accessed to the distributed power supply device, wherein the power capacity of each distributed power supply device is the same.
After the second function constraint after the relaxation processing is obtained, solving an overall planning function through a Cplex solver based on the first function constraint, the third function constraint and the second function constraint after the relaxation processing to obtain the number of the distributed power equipment corresponding to the distributed power equipment node to be accessed in the node distribution network structure, then obtaining the power capacity of one distributed power equipment, and multiplying the power capacity of the distributed power equipment by the number of the equipment corresponding to each distributed power equipment node to be accessed to obtain the power capacity of the distributed power equipment node to be accessed, wherein the power capacities of all the distributed power equipment are the same.
Further, if the power capacities of the distributed power devices are different, the power capacity of each distributed power device may also be obtained, so as to determine, according to the power capacity of the node to be accessed to the distributed power device, which distributed power device needs to be accessed to the node to be accessed to the distributed power device.
According to the embodiment, the power capacity of each distributed power supply device is the same, and the speed of obtaining the power capacity of the nodes of the distributed power supply devices to be connected in the node distribution network structure is improved.
The present invention also provides a configuration apparatus for an incremental distributed power supply device, and referring to fig. 3, the configuration apparatus for an incremental distributed power supply device includes:
the incremental distribution network model building method comprises a building module 10, a calculation module and a calculation module, wherein the building module is used for building a distributed power incremental distribution network model in an incremental distribution network system, and the incremental distribution network model at least comprises a node distribution network structure;
the obtaining module 20 is configured to obtain, from a preset database of the incremental power distribution network system, first target data related to a distributed power source node, second target data related to the incremental power distribution network node, and third target data related to a power node in the incremental power distribution network model;
the building module 10 is further configured to build a first objective function corresponding to the distributed power source node and a first function constraint that constrains the first objective function according to the first objective data, build a second objective function corresponding to the incremental power distribution network node and a second function constraint that constrains the second objective function according to the second objective data, and build a third objective function corresponding to the power node and a third function constraint that constrains the third objective function according to the third objective data, where the first function constraint includes a constraint that a preset node in a node distribution network structure accesses the number of distributed power source devices, a distributed power source device permeability upper limit constraint, and a distributed power source device total power source capacity constraint, and the second function constraint includes a power flow constraint, a voltage constraint, and a power constraint of each line in the incremental power distribution network model, the third function constraint comprises a power load amount constraint transferred at a target moment and a power load amount constraint transferred at the target moment;
a calculating module 30, configured to multiply the first objective function by a first weight corresponding to the first objective function to obtain a first weight function, multiply the second objective function by a second weight corresponding to the second objective function to obtain a second weight function, and multiply the third objective function by a third weight corresponding to the third objective function to obtain a third weight function; adding the first weight function, the second weight function and the third weight function to obtain an overall planning function;
the relaxation processing module 40 is configured to perform second-order cone relaxation processing on the second function constraint to obtain a second function constraint after the relaxation processing;
the solving module 50 is configured to solve the overall planning function through a solver based on the first function constraint, the third function constraint and the relaxed second function constraint to obtain a power supply capacity of a node to be connected to the distributed power supply device in the node distribution network structure;
and a configuration module 60, 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, an auxiliary variable U is definedi.t′＝(Ui.t)2And Iij.t′＝(Iij.t)2From the auxiliary variables, the following target equation is derived:
Figure GDA0003452711880000211
Iij.trepresenting the line current, U, between node i and node j in a node-distribution network architecturei.tRepresenting the voltage amplitude, P, of node i at time tij.tRepresenting the active electric quantity, Q, transmitted by branch ij at time tij.tRepresenting the reactive power transmitted by the branch ij at the moment t;
and carrying out relaxation treatment on the target equation to obtain the following relaxation inequality:
Figure GDA0003452711880000212
reconstructing the relaxation inequality through a second order cone to obtain the following target inequality:
Figure GDA0003452711880000213
performing second-order cone relaxation processing on each line power flow constraint in the second function constraint according to the target inequality to obtain a second function constraint after relaxation processing;
the expression of each line power flow constraint in the relaxed second function constraint is as follows:
Figure GDA0003452711880000214
Figure GDA0003452711880000215
Figure GDA0003452711880000216
Figure GDA0003452711880000221
Figure GDA0003452711880000222
wherein, Pi.tRepresents the active electric quantity, Q, of the node i at the moment ti.tRepresenting the reactive power of the node i at time t, Pj.tRepresents the active electric quantity, Q, of the node j at the moment tj.tRepresenting nodesj reactive power at time t, Ui.tRepresenting the voltage amplitude, U, of node i at time tj.tRepresenting the voltage magnitude, R, of node j at time tijRepresenting the resistance of the line between node i and node j, XijRepresenting the reactance of the line between node i and node j, GijRepresenting the conductance of the line between node i and node j, BijRepresenting susceptance of the line between node i and node j;
Figure GDA0003452711880000223
represents the voltage phase angle difference between node i and node j; u (j) represents an upstream set of nodes connected to node j, d (j) represents a downstream set of nodes connected to node j, Pjl.tRepresents the active electric quantity between the node j and the node l at the moment t, Qjl.tAnd the reactive power between the node j and the node l at the moment t is represented.
Further, the constraint of the number of the preset nodes accessing the distributed power supply devices is represented as: n is a radical ofi.min≤Ni≤Ni.max；
Wherein N isi.minMinimum value, N, representing the number of devices accessing the distributed power supply at node ii.maxMaximum value N representing the number of devices accessing the distributed power supply at preset node iiThe number of the distributed power supply devices accessed by the ith node is represented;
the upper limit constraint on the permeability of the distributed power equipment is expressed as:
Figure GDA0003452711880000224
wherein, beta represents the permeability of the distributed power supply equipment after being connected to the node, PloadRepresenting the electric load quantity of the distributed power supply equipment at a preset node; n isiRepresents a variable with a value of 0 or 1, ni0 represents that the ith node in the node distribution network structure is not connected with the distributed power supply equipment, and ni1 represents that the ith node is accessed into the distributed power supply equipment; psg.DGRepresenting the rated active power of each distributed power supply device;
the total power capacity constraint for a distributed power plant is expressed as: pmin.DG≤Pt.DG≤Pmax.DG；
Wherein, Pmin.DGLower limit, P, representing total power supply capacity of distributed power supply equipmentmax.DGAn upper limit value representing a total power capacity of the distributed power equipment; pt.DGRepresenting the active power provided by the distributed power equipment at time t.
Further, the expression of each line power flow constraint is as follows:
Figure GDA0003452711880000225
Figure GDA0003452711880000231
Figure GDA0003452711880000232
Figure GDA0003452711880000233
Figure GDA0003452711880000234
wherein, Pi.tRepresents the active electric quantity, Q, of the node i at the moment ti.tRepresenting the reactive power of the node i at time t, Pj.tThe expression is the active electric quantity, Q, of the node j at the moment tj.tThe expression is the reactive power of the node j at the time t, Pij.tRepresenting the active electric quantity, Q, transmitted by branch ij at time tij.tRepresenting the reactive electric quantity, U, transmitted by branch ij at time ti.tRepresenting the voltage amplitude, U, of node i at time tj.tRepresenting the voltage magnitude, R, of node j at time tijRepresenting the resistance of the line between node i and node j, XijRepresenting the reactance of the line between node i and node j, GijRepresenting node i and node jConductance of the line between, BijRepresenting susceptance of the line between node i and node j;
Figure GDA0003452711880000235
represents the voltage phase angle difference between node i and node j; u (j) represents an upstream set of nodes connected to node j, d (j) represents a downstream set of nodes connected to node j, Pjl.tRepresents the active electric quantity between the node j and the node l at the moment t, Qjl.tRepresenting the reactive power between the node j and the node l at the time t;
the voltage constraints are: u shapei.min≤Ui.t≤Ui.maxWherein, Ui.minRepresenting the minimum value, U, of the voltage amplitude of node ii.maxRepresents the maximum value of the voltage amplitude of the node i;
the power constraint is: pij.t≤Pij.maxWherein P isij.tRepresenting the power flowing on the line between node i and node j at time t, Pij.maxWhich represents the maximum value of the power flowing on the line between node i and node j.
Further, the target time is time t, and the expressions of the electric power load amount constraint shifted out at the target time and the electric power load amount constraint shifted out at the target time are as follows:
Figure GDA0003452711880000236
where ρ isminPt.loadRepresents the minimum value of the power coefficient of the load transfer at time t, rhomaxPt.loadThe maximum value of the power coefficient of the load transfer at the moment t is shown; sigmaminPt.loadRepresenting the minimum value of the power coefficient, σ, of the load transfer at time tmaxPt.loadRepresenting the maximum value of the power coefficient of load transfer at the moment t; pt.outRepresenting the electric load amount transferred out at the time t; pt.inIndicating the amount of electrical load transferred at time t.
Further, the solving module 50 includes:
the solving unit is used for solving the overall planning function through a Cplex solver based on the first function constraint, the third function constraint and the relaxed second function constraint to obtain the number of equipment to be connected to the distributed power equipment in the node distribution network structure, wherein the equipment corresponds to the distributed power equipment;
the device comprises an acquisition unit, a control unit and a control unit, wherein the acquisition unit is used for acquiring the power supply capacity of one distributed power supply device;
and the calculating unit is used for multiplying the power capacity of the distributed power supply equipment by the number of the equipment to obtain the power capacity of the nodes of the distributed power supply equipment to be accessed, wherein the power capacity of each distributed power supply equipment is the same.
The specific implementation of the configuration apparatus of incremental distributed power supply devices of the present invention is substantially the same as each embodiment of the configuration method of incremental distributed power supply devices, and will not be described again here.
The invention also proposes a computer-readable storage medium on which a detection program is stored, which, when executed by a processor, implements the steps of the method for configuring an incremental distributed power supply device as described above.
The specific implementation of the computer-readable storage medium of the present invention is substantially the same as the above-mentioned embodiments of the configuration method of the incremental distributed power supply device, and will not be described again here.
It will be appreciated by those skilled in the art that the above-described preferred embodiments may be freely combined, superimposed, without conflict.
It will be understood that the embodiments described above are illustrative only and not restrictive, and that various obvious and equivalent modifications and substitutions for details described herein may be made by those skilled in the art without departing from the basic principles of the invention.

Claims (10)

1. A configuration method of an incremental distributed power supply device is characterized by comprising the following steps:
s100, constructing a distributed power incremental distribution network model in an incremental distribution network system, wherein the incremental distribution network model at least comprises a node distribution network structure;
s200, obtaining first target data related to distributed power supply nodes, second target data related to incremental power distribution network nodes and third target data related to power nodes in the incremental power distribution network model from a preset database of the incremental power distribution network system, constructing a first target function corresponding to the distributed power supply nodes and a first function constraint constraining the first target function according to the first target data, constructing a second target function corresponding to the incremental power distribution network nodes and a second function constraint constraining the second target function according to the second target data, and constructing a third target function corresponding to the power nodes and a third function constraint constraining the third target function according to the third target data, wherein the first function constraint comprises constraint, access number of distributed power supply equipment of preset nodes in a distribution network structure, The second function constraint comprises power flow constraint, voltage constraint and power constraint of each line in the incremental distribution network model, and the third function constraint comprises power load constraint transferred at a target moment and power load constraint transferred at a target moment; in step S200, the first objective function expression is a formula:
max WDG(ni,Ni)＝WS.DG-WI.DG-WOM.DG；
wherein W in the formulaS.DGIs the formula:
Figure FDA0003452711870000011
w in the formulaI.DGIs the formula:
Figure FDA0003452711870000012
in the formulaW ofOM.DGIs the formula:
Figure FDA0003452711870000013
nirepresents a variable with a value of 0 or 1, niWhen the node is 0, the ith node in the node distribution network structure is not connected with the distributed power supply device, ni1 represents that the ith node is accessed into the distributed power supply equipment; n is a radical ofiThe number of the distributed power supply devices accessed by the ith node is represented; wS.DGRepresents the electricity sales income within the preset years of all distributed power supply equipment in the node distribution network structure, WI.DGRepresents the early construction cost, W, of all distributed power supply equipment in the node distribution network structure before useOM.DGRepresents the operation and maintenance cost, U, of all distributed power supply equipment within the preset agetRepresents the set of all the time in one day, which is 24 hours; q. q.sesThe unit electricity selling price of the distributed power supply equipment is represented, and the unit electricity selling price of each distributed power supply equipment is the same; pt.DGRepresenting the active power provided by the distributed power supply equipment at the moment t; q. q ofsgRepresents a distributed power supply equipment construction cost, specifically, qsgThe construction investment of the distributed power supply equipment corresponding to the unit power supply capacity is represented; u shapeiRepresenting that a node set to be accessed to distributed power equipment is preset in a node distribution network structure; psg.DGThe rated active power of each distributed power supply device is represented, and the rated active power of each distributed power supply device is the same; y represents the life cycle of each distributed power supply device, and the life cycle of each distributed power supply device is the same; a represents the discount rate; q. q.somThe operation and maintenance cost of unit power supply of the distributed power supply equipment is represented; in step S200, the expression of the second objective function is formula:
formula max WDN(pi)＝WS.DN-WL.DN-WE.DN-WB1.DN-WB2.DN-WB3.DG；
Wherein W in the formulaS.DGIs the formula:
the formula:
Figure FDA0003452711870000021
w in the formulaL.DNIs the formula:
Figure FDA0003452711870000022
wherein p isiVariable of an objective function, W, representing a node of an incremental distribution networkS.DNRepresents the electricity selling income of the incremental distribution network nodes within the preset year, WL.DNRepresenting the line loss cost, W, over a predetermined period of timeE.DNIndicating a fault repair cost, W, over a predetermined period of timeB1.DNRepresents the electricity purchase cost from the upper level within a preset age, WB2.DNRepresenting the purchase of electricity from distributed power nodes within a predetermined age, WB3.DGRepresents the loss, U, caused by the output fluctuation of the distributed power supply equipmenttRepresenting the set of all times of the day, fesRepresenting a unit price for electricity sold to the power node; pt.loadRepresenting the initial electric power load quantity at the time t; pt.outRepresenting the electric load amount transferred out at the time t; pt.inRepresenting the amount of electric power load transferred at time t; p ist.isRepresenting the amount of power load which can be interrupted at the time t; pt.lossRepresenting the active loss at the moment t; EENStThe expected value of less power supply energy to the power node at the moment t is shown; u shapebRepresenting a set of all lines of a power grid in the incremental power distribution grid model; lambda [ alpha ]bA probability value representing the fault occurrence of the b-th line; u shapenRepresenting that a node set to be accessed to distributed power equipment is preset in a node distribution network structure; pn.t.loadRepresenting the initial electric load quantity of the node n at the time t; f. ofeb1Represents a unit purchase price from an upper-level power grid; f. ofeb2Represents a unit purchase price from the distributed power node; ppRepresenting punishment cost required to be paid due to random fluctuation within a preset age limit; rhosRepresenting the probability of power loss caused by the output fluctuation of the distributed power supply equipment within a preset period; delta QsRepresenting the random fluctuation amount of the distributed power supply equipment during the output fluctuation;
w in the formulaE.DNIs the formula:
Figure FDA0003452711870000031
w in the formulaB1.DNIs the formula:
Figure FDA0003452711870000032
w in the formulaB2.DNIs the formula:
Figure FDA0003452711870000033
w in the formulaB3.DGIs the formula:
Figure FDA0003452711870000034
in step S200, the expression of the third objective function is formula:
Figure FDA0003452711870000035
wherein, γebRepresenting the unit price, P, of electricity purchased by a power node to an incremental distribution network nodet.outRepresenting the electric load amount transferred out at the time t; p ist.inRepresenting the amount of electric power load transferred at time t; pt.isRepresenting the amount of electric load which can be interrupted at the moment t; pt.lossRepresenting the active loss, U, at time ttRepresenting the set of all times of the day;
s300, multiplying the first objective function by a first weight corresponding to the first objective function to obtain a first weight function, multiplying the second objective function by a second weight corresponding to the second objective function to obtain a second weight function, and multiplying the third objective function by a third weight corresponding to the third objective function to obtain a 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, performing second-order cone relaxation processing on the second function constraint to obtain a second function constraint after relaxation processing, solving 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 to obtain the power supply capacity of the node to be connected to the distributed power supply equipment in the node distribution network structure, and configuring the distributed power supply equipment corresponding to the node according to the power supply capacity.
2. The method for configuring incremental distributed power equipment according to claim 1, wherein in step S500, an auxiliary variable U is defined during a second-order cone relaxation processi.t′＝(Ui.t)2And Iij.t′＝(Iij.t)2From the auxiliary variables, the following target equation is derived:
Figure FDA0003452711870000041
Iij.trepresenting the line current, U, between node i and node j in a node-distribution network architecturei.tRepresenting the voltage amplitude, P, of node i at time tij.tRepresenting the active electric quantity, Q, transmitted by branch ij at time tij.tRepresenting the reactive power transmitted by the branch ij at the moment t;
and carrying out relaxation treatment on the target equation to obtain the following relaxation inequality:
Figure FDA0003452711870000042
reconstructing the relaxation inequality through a second order cone to obtain the following target inequality:
Figure FDA0003452711870000043
performing second-order cone relaxation processing on each line power flow constraint in the second function constraint according to the target inequality to obtain a second function constraint after relaxation processing;
the expression of each line power flow constraint in the relaxed second function constraint is as follows:
Figure FDA0003452711870000044
Figure FDA0003452711870000045
Figure FDA0003452711870000051
Figure FDA0003452711870000052
(Uj.t)2＝Ui.t′-2(RijPij.t+XijQij.t)+[(Rij)2+(Xij)2]·Iij.t′；
wherein, Pi.tRepresents the active electric quantity, Q, of the node i at the moment ti.tRepresenting the reactive power of the node i at time t, Pj.tRepresents the active electric quantity, Q, of the node j at the moment tj.tRepresenting the reactive electric quantity, U, of the node j at time ti.tRepresenting the voltage amplitude, U, of node i at time tj.tRepresenting the voltage magnitude, R, of node j at time tijRepresenting the resistance of the line between node i and node j, XijRepresenting the reactance of the line between node i and node j, GijRepresenting lines between nodes i and jConductance, BijRepresenting susceptance of the line between node i and node j;
Figure FDA0003452711870000053
represents the voltage phase angle difference between node i and node j; u (j) represents an upstream set of nodes connected to node j, d (j) represents a downstream set of nodes connected to node j, Pjl.tRepresents the active electric quantity between the node j and the node l at the moment t, Qjl.tAnd the reactive power between the node j and the node l at the moment t is represented.
3. The method of configuring an incremental distributed power apparatus of claim 1, wherein in said step S200,
the constraint of the number of the preset nodes accessing the distributed power supply devices is represented as follows: n is a radical ofi.min≤Ni≤Ni.max；
Wherein N isi.minMinimum value, N, representing the number of devices accessing the distributed power supply at node ii.maxMaximum value N representing the number of devices accessing the distributed power supply at preset node iiThe number of the distributed power supply devices accessed by the ith node is represented;
the upper limit constraint on the permeability of the distributed power equipment is expressed as:
Figure FDA0003452711870000054
wherein, beta represents the permeability of the distributed power supply equipment after being connected to the node, PloadRepresenting the electric load quantity of the distributed power supply equipment at a preset node; n isiRepresents a variable with a value of 0 or 1, niWhen the node is 0, the ith node in the node distribution network structure is not connected with the distributed power supply device, ni1 represents that the ith node is accessed into the distributed power supply equipment; p issgDG represents the rated active power of each distributed power supply device;
the total power capacity constraint for a distributed power plant is expressed as: pmin.DG≤Pt.DG≤Pmax.DG；
Wherein, Pmin.DGLower limit, P, representing total power supply capacity of distributed power supply equipmentmax.DGAn upper limit value representing a total power capacity of the distributed power equipment; pt.DGRepresenting the active power provided by the distributed power equipment at time t.
4. The method of configuring an incremental distributed power apparatus of claim 1, wherein in said step S200,
the expression of each line tide constraint is as follows:
Figure FDA0003452711870000061
Figure FDA0003452711870000062
Figure FDA0003452711870000063
Figure FDA0003452711870000064
Figure FDA0003452711870000065
wherein, Pi.tRepresents the active electric quantity, Q, of the node i at the moment ti.tRepresenting the reactive power of the node i at time t, Pj.tThe expression is the active electric quantity, Q, of the node j at the moment tj.tThe expression is the reactive power of the node j at the time t, Pij.tRepresenting the active electric quantity, Q, transmitted by branch ij at time tij.tRepresenting the reactive electric quantity, U, transmitted by branch ij at time ti.tRepresenting the voltage magnitude, U, of node i at time tj.tRepresenting the voltage magnitude, R, of node j at time tijRepresenting the resistance of the line between node i and node j, XijRepresenting the reactance of the line between node i and node j, GijRepresenting the conductance of the line between node i and node j, BijRepresenting susceptance of the line between node i and node j;
Figure FDA0003452711870000066
represents the voltage phase angle difference between node i and node j; u (j) represents an upstream set of nodes connected to node j, d (j) represents a downstream set of nodes connected to node j, Pjl.tRepresents the active electric quantity between the node j and the node l at the moment t, Qjl.tRepresenting the reactive power between the node j and the node l at the time t;
the voltage constraints are: u shapei.min≤Ui.t≤Ui.maxWherein, Ui.minRepresenting the minimum value, U, of the voltage amplitude of node ii.maxRepresents the maximum value of the voltage amplitude of the node i;
the power constraint is: pij.t≤Pij.maxWherein P isij.tRepresenting the power flowing on the line between node i and node j at time t, Pij.maxWhich represents the maximum value of the power flowing on the line between node i and node j.
5. The method for configuring an incremental distributed power apparatus according to any one of claims 1 to 4, wherein in the step S200,
the target time is t time, and the expressions of the electric power load quantity constraint transferred at the target time and the electric power load quantity constraint transferred at the target time are as follows:
Figure FDA0003452711870000071
where ρ isminPt.loadRepresents the minimum value of the power coefficient of the load transfer at time t, rhomaxPt.loadThe maximum value of the power coefficient of the load transfer at the moment t is shown; sigmaminPt.loadRepresenting the minimum value of the power coefficient, σ, of the load transfer at time tmaxPt.loadRepresenting the maximum value of the power coefficient of load transfer at the moment t; pt.outRepresenting the electric load amount transferred out at the time t; pt.inIndicating the amount of electrical load transferred at time t.
6. The method according to claim 5, wherein in step S500, the step of obtaining the power capacity of the node to be connected to the distributed power supply device in the node distribution network structure by solving the overall planning function through a solver based on the first function constraint, the third function constraint, and the relaxed second function constraint includes:
based on the first function constraint, the third function constraint and the relaxed second function constraint, solving the overall planning function through a Cplex solver to obtain the number of equipment to be connected to the distributed power equipment in the node distribution network structure, wherein the equipment corresponds to the distributed power equipment;
the method comprises the steps of obtaining the power capacity of one distributed power supply device, and multiplying the power capacity of the distributed power supply device by the number of the devices to obtain the power capacity of nodes to be connected to the distributed power supply device, wherein the power capacity of each distributed power supply device is the same.
7. An incremental distributed power supply device configuration apparatus, characterized in that the incremental distributed power supply device configuration apparatus comprises:
the incremental distribution network model comprises a building module and a control module, wherein the building module is used for building a distributed power incremental distribution network model in an incremental distribution network system, and the incremental distribution network model at least comprises a node distribution network structure;
the acquisition module is used for acquiring first target data related to distributed power source nodes, second target data related to incremental power distribution network nodes and third target data related to power nodes in the incremental power distribution network model from a preset database of the incremental power distribution network system;
the building module is further configured to build a first objective function corresponding to the distributed power supply node and a first function constraint constraining the first objective function according to the first target data, build a second objective function corresponding to the incremental power distribution network node and a second function constraint constraining the second objective function according to the second target data, and build a third objective function corresponding to the power node and a third function constraint constraining the third objective function according to the third target data, where the first function constraint includes a constraint on the number of distributed power supply devices accessed by preset nodes in a node distribution network structure, a distributed power supply device permeability upper limit constraint and a distributed power supply device total power supply capacity constraint, and the second function constraint includes a power flow constraint, a voltage constraint and a power constraint of each line in the incremental power distribution network model, the third function constraint comprises a power load amount constraint transferred at a target moment and a power load amount constraint transferred at the target moment; wherein the first objective function expression is a formula:
formula max WDG(ni,Ni)＝WS.DG-WI.DG-WOM.DG；
Wherein W in the formulaS.DGIs the formula:
Figure FDA0003452711870000081
w in the formulaI.DGIs the formula:
Figure FDA0003452711870000082
w in the formulaOM.DGIs the formula:
Figure FDA0003452711870000083
nirepresents a variable with a value of 0 or 1, niWhen the node is 0, the ith node in the node distribution network structure is not connected with the distributed power supply device, ni1 represents that the ith node is accessed into the distributed power supply equipment; n is a radical ofiIs shown asThe number of the i nodes connected to the distributed power supply equipment; wS.DGRepresents the electricity sales income within the preset years of all distributed power supply equipment in the node distribution network structure, WI.DGRepresents the upfront construction cost, W, of all distributed power equipment in the node distribution network structure before useOM.DGRepresents the operation and maintenance cost, U, of all distributed power supply equipment within the preset agetRepresents the set of all the time in one day, which is 24 hours; q. q ofesThe unit electricity selling price of the distributed power supply equipment is represented, and the unit electricity selling price of each distributed power supply equipment is the same; pt.DGRepresenting the active power provided by the distributed power supply equipment at the moment t; q. q.ssgRepresents a distributed power supply equipment construction cost, specifically, qsgThe construction investment of the distributed power supply equipment corresponding to the unit power supply capacity is represented; u shapeiRepresenting that a node set to be accessed to distributed power equipment is preset in a node distribution network structure; psg.DGThe rated active power of each distributed power supply device is represented, and the rated active power of each distributed power supply device is the same; y represents the life cycle of each distributed power supply device, and the life cycle of each distributed power supply device is the same; a represents the discount rate; q. q.somThe operation and maintenance cost of unit power supply of the distributed power supply equipment is represented; in step S200, the expression of the second objective function is formula:
formula max WDN(pi)＝WS.DN-WL.DN-WE.DN-WB1.DN-WB2.DN-WB3.DG；
Wherein W in the formulaS.DGIs the formula:
the formula:
Figure FDA0003452711870000091
w in the formulaL.DNIs the formula:
Figure FDA0003452711870000092
wherein p isiObjective function variable representing incremental distribution network node，WS.DNRepresents the electricity selling income of the incremental distribution network nodes within the preset year, WL.DNRepresenting the line loss cost, W, over a predetermined period of timeE.DNIndicating a fault repair cost, W, over a predetermined period of timeB1.DNRepresents the electricity purchase cost from the upper level within a preset age, WB2.DNRepresenting the purchase of electricity from distributed power nodes within a predetermined age, WB3.DGRepresents the loss, U, caused by the output fluctuation of the distributed power supply equipmenttRepresenting the set of all times of the day, fesRepresenting a unit price for electricity sold to the power node; pt.loadRepresenting the initial electric power load quantity at the time t; pt.outRepresenting the electric load amount transferred out at the time t; pt.inRepresenting the amount of electric power load transferred at time t; pt.isRepresenting the amount of electric load which can be interrupted at the moment t; pt.lossRepresenting the active loss at the moment t; EENStThe expected value of less power supply energy to the power node at the moment t is shown; u shapebRepresenting a set of all lines of a power grid in the incremental power distribution grid model; lambda [ alpha ]bA probability value representing the fault occurrence of the b-th line; u shapenRepresenting that a node set to be accessed to distributed power equipment is preset in a node distribution network structure; pn.t.loadRepresenting the initial electric load quantity of the node n at the time t; f. ofeb1Represents a unit purchase price from an upper-level power grid; f. ofeb2Represents a unit purchase price from the distributed power node; ppRepresenting punishment cost required to be paid due to random fluctuation within a preset age; rhosRepresenting the probability of power loss caused by the output fluctuation of the distributed power supply equipment within a preset period; delta QsRepresenting the random fluctuation amount of the distributed power supply equipment during the output fluctuation;
w in the formulaE.DNIs the formula:
formula (la)
Figure FDA0003452711870000093
W in the formulaB1.DNThe expression of (b) is the formula:
formula (II)
Figure FDA0003452711870000101
W in the formulaB2.DNThe expression of (b) is the formula:
Figure FDA0003452711870000102
w in the formulaB3.DGIs the formula:
Figure FDA0003452711870000103
in step S200, the expression of the third objective function is formula:
formula (II)
Figure FDA0003452711870000104
Wherein, γebRepresenting the unit price, P, of electricity purchased by a power node to an incremental distribution network nodet.outRepresenting the electric load amount transferred out at the time t; pt.inRepresenting the amount of electric power load transferred at time t; p ist.isRepresenting the amount of electric load which can be interrupted at the moment t; pt.lossRepresenting the active loss, U, at time ttRepresenting the set of all times of the day;
a calculating module, configured to multiply the first objective function by a first weight corresponding to the first objective function to obtain a first weight function, multiply the second objective function by a second weight corresponding to the second objective function to obtain a second weight function, and multiply the third objective function by a third weight corresponding to the third objective function to obtain a third weight function; adding the first weight function, the second weight function and the third weight function to obtain an overall planning function;
the relaxation processing module is used for performing second-order cone relaxation processing on the second function constraint to obtain a second function constraint after the relaxation processing;
the solving module is used for solving the overall planning function through a solver based on the first function constraint, the third function constraint and the relaxed second function constraint to obtain the power supply capacity of the node to be connected to the distributed power supply equipment in the node distribution network structure;
and the configuration module is used for configuring the distributed power supply equipment corresponding to the node according to the power supply capacity.
8. The incremental distributed power apparatus configuration device according to claim 7, wherein an auxiliary variable U is defined during the second order cone relaxation processi.t′＝(Ui.t)2And Iij.t′＝(Iij.t)2From the auxiliary variables, the following target equation is derived:
Figure FDA0003452711870000111
Iij.trepresenting the line current, U, between node i and node j in a node-distribution network architecturei.tRepresenting the voltage amplitude, P, of node i at time tij.tRepresenting the active electric quantity, Q, transmitted by branch ij at time tij.tRepresenting the reactive power transmitted by the branch ij at the moment t;
and carrying out relaxation treatment on the target equation to obtain the following relaxation inequality:
Figure FDA0003452711870000112
reconstructing the relaxation inequality through a second order cone to obtain the following target inequality:
Figure FDA0003452711870000113
performing second-order cone relaxation processing on each line power flow constraint in the second function constraint according to the target inequality to obtain a second function constraint after relaxation processing;
the expression of each line power flow constraint in the relaxed second function constraint is as follows:
Figure FDA0003452711870000114
Figure FDA0003452711870000115
Figure FDA0003452711870000116
Figure FDA0003452711870000117
(Uj.t)2＝Ui.t′-2(RijPij.t+XijQij.t)+[(Rij)2+(Xij)2]·Iij.t′；
wherein, Pi.tRepresents the active electric quantity, Q, of the node i at the moment ti.tRepresenting the reactive power of the node i at time t, Pj.tRepresents the active electric quantity, Q, of the node j at the moment tj.tRepresenting the reactive electric quantity, U, of the node j at time ti.tRepresenting the voltage magnitude, U, of node i at time tj.tRepresenting the voltage magnitude, R, of node j at time tijRepresenting the resistance of the line between node i and node j, XijRepresenting the reactance of the line between node i and node j, GijRepresenting the conductance of the line between node i and node j, BijRepresenting susceptance of the line between node i and node j;
Figure FDA0003452711870000122
represents the voltage phase angle difference between node i and node j; u (j) represents an upstream set of nodes connected to node j, d (j) represents a downstream set of nodes connected to node j, Pjl.tRepresenting the active electric quantity between the node j and the node l at the moment t，Qjl.tAnd the reactive power between the node j and the node l at the moment t is represented.
9. The incremental distributed power supply device configuration apparatus according to claim 7, wherein the constraint on the number of preset nodes accessing the distributed power supply device is represented as: n is a radical ofi.min≤Ni≤Ni.max；
Wherein N isi.minMinimum value, N, representing the number of devices accessing the distributed power supply at node ii.maxMaximum value N representing the number of devices accessing the distributed power supply at preset node iiThe number of the distributed power supply devices accessed by the ith node is represented;
the upper limit constraint on the permeability of the distributed power equipment is expressed as:
Figure FDA0003452711870000121
wherein, beta represents the permeability of the distributed power supply equipment after being connected to the node, PloadRepresenting the electric load quantity of the distributed power supply equipment at a preset node; n isiRepresents a variable with a value of 0 or 1, niWhen the node is 0, the ith node in the node distribution network structure is not connected with the distributed power supply device, ni1 represents that the ith node is accessed into the distributed power supply equipment; psgDG represents the rated active power of each distributed power supply device;
the total power capacity constraint of a distributed power device is expressed as: pmin.DG≤Pt.DG≤Pmax.DG；
Wherein, Pmin.DGLower limit, P, representing total power supply capacity of distributed power supply equipmentmax.DGAn upper limit value representing a total power capacity of the distributed power equipment; pt.DGRepresenting the active power provided by the distributed power equipment at time t.
10. A computer-readable storage medium, characterized in that the computer-readable storage medium has stored thereon a detection program which, when executed by a processor, carries out the steps of the method of configuration of an incremental distributed power supply device according to any one of claims 1 to 6.
CN202010818367.XA 2020-08-14 2020-08-14 Configuration method and device of incremental distributed power supply equipment and readable storage medium Active CN111952968B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010818367.XA CN111952968B (en) 2020-08-14 2020-08-14 Configuration method and device of incremental distributed power supply equipment and readable storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010818367.XA CN111952968B (en) 2020-08-14 2020-08-14 Configuration method and device of incremental distributed power supply equipment and readable storage medium

Publications (2)

Publication Number Publication Date
CN111952968A CN111952968A (en) 2020-11-17
CN111952968B true CN111952968B (en) 2022-05-03

Family

ID=73343341

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010818367.XA Active CN111952968B (en) 2020-08-14 2020-08-14 Configuration method and device of incremental distributed power supply equipment and readable storage medium

Country Status (1)

Country Link
CN (1) CN111952968B (en)

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8401709B2 (en) * 2009-11-03 2013-03-19 Spirae, Inc. Dynamic distributed power grid control system
CN103094900B (en) * 2012-12-20 2014-10-08 河海大学 Distributed generation power distribution network three-phase load flow calculation method taking phase sequence mixing method into consideration

Also Published As

Publication number Publication date
CN111952968A (en) 2020-11-17

Similar Documents

Publication Publication Date Title
CN109034563B (en) Multi-subject game incremental power distribution network source-network-load collaborative planning method
CN110276698B (en) Distributed renewable energy transaction decision method based on multi-agent double-layer collaborative reinforcement learning
Yang et al. A multi-objective stochastic optimization model for electricity retailers with energy storage system considering uncertainty and demand response
Chen et al. Research on day-ahead transactions between multi-microgrid based on cooperative game model
CN111082451B (en) Incremental distribution network multi-objective optimization scheduling model based on scene method
CN110782281B (en) Day-ahead market clearing method for multi-owner cascade power station basin electric quantity transfer
Yuan et al. Real-time pricing for smart grid with multi-energy microgrids and uncertain loads: a bilevel programming method
CN110443398A (en) A kind of regional complex energy resource system optimizing operation method based on repeated game model
CN110378729A (en) A kind of integration requirement response method based on dynamic energy price strategy
CN110429653A (en) Consider energy storage and the rural power grids distributed photovoltaic consumption method and terminal device of DR
CN112101607A (en) Active power distribution network rolling optimization scheduling method considering demand response time effect
CN109272353B (en) System dynamic probability energy flow analysis method considering comprehensive demand response uncertainty
CN112290536A (en) Online scheduling method of electricity-heat comprehensive energy system based on near-end strategy optimization
CN117494993A (en) Ladder carbon transaction-based source network load distribution collaborative low-carbon scheduling method
Sun et al. A coalitional game theoretic energy transaction algorithm for networked microgrids
Tang et al. Unit maintenance strategy considering the uncertainty of energy intensive load and wind power under the carbon peak and carbon neutral target
Fan et al. Medium and long-term electricity trading considering renewable energy participation
CN110556821B (en) Multi-microgrid double-layer optimization scheduling method considering interactive power control and bilateral bidding transaction
CN111952968B (en) Configuration method and device of incremental distributed power supply equipment and readable storage medium
CN115940299B (en) Active and reactive collaborative optimization method for power distribution system source network load storage
Fan et al. Optimization of power supply capacity of distribution network considering the participation of power sales companies in spot power trading
CN114069613B (en) Method and system for regulating and controlling participation peak shaving of self-contained power plant based on enterprise energy consumption characteristics
Yan et al. Consensus alternating direction multiplier method based fully distributed peer-to-peer energy transactions considering the network transmission distance
CN110752599B (en) Distributed power supply grid-connected configuration method
CN112529610A (en) End-to-end electric energy trading market user decision method based on reinforcement learning

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant