CN111952968A - 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
CN111952968A
CN111952968A CN202010818367.XA CN202010818367A CN111952968A CN 111952968 A CN111952968 A CN 111952968A CN 202010818367 A CN202010818367 A CN 202010818367A CN 111952968 A CN111952968 A CN 111952968A
Authority
CN
China
Prior art keywords
node
function
constraint
power supply
representing
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.)
Granted
Application number
CN202010818367.XA
Other languages
Chinese (zh)
Other versions
CN111952968B (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)
  • Computer Hardware Design (AREA)
  • Power Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computer Networks & Wireless Communication (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 in China and the development of the economic society, an incremental power distribution network is indispensable, the planning and operation mode of the traditional power distribution network needs to be changed urgently, and the incremental power distribution network is combined with the latest power distribution network operation technology, so that the goals of reducing cost, improving quality and efficiency, saving energy and reducing emission are favorably realized. 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 for 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 for 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 for constraining the third target function according to the third target data, wherein the first function constraint comprises constraint for limiting the number of distributed power supply equipment accessed by preset nodes in a distribution network node 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 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 capacity of the capacitor.
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)2From the auxiliary variables, the following target equation is derived:
Figure BDA0002633565000000021
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 BDA0002633565000000031
reconstructing the relaxation inequality through a second order cone to obtain the following target inequality:
Figure BDA0002633565000000032
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 BDA0002633565000000033
Figure BDA0002633565000000034
Figure BDA0002633565000000035
Figure BDA0002633565000000036
(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 BDA0002633565000000037
represents node i anda voltage phase angle difference between nodes 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.
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 BDA0002633565000000041
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.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.
Preferably, in the step S200,
the expression of each line tide constraint is as follows:
Figure BDA0002633565000000042
Figure BDA0002633565000000043
Figure BDA0002633565000000044
Figure BDA0002633565000000045
Figure BDA0002633565000000046
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 BDA0002633565000000051
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 power between node j and node l at time tElectric quantity, 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.
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 BDA0002633565000000052
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.
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 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.
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 sub-distributed power supply nodes and a first function constraint constraining the first objective function according to the first objective data, build a second objective function corresponding to the incremental power distribution network nodes and a second function constraint constraining 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 constraining the third objective function according to the third objective 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 capacitance 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.
[ PROBLEMS ] the present invention
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 for 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 for 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 for constraining the third target function according to the third target data, wherein the first function constraint comprises constraint for limiting the number of distributed power supply equipment accessed by preset nodes in a distribution network node 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.
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 BDA0002633565000000091
w in formula (1)I.DGIs formula (3):
Figure BDA0002633565000000092
w in formula (1)OM.DGIs formula (4):
Figure BDA0002633565000000101
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.ssgRepresenting a distributionConstruction costs of the power supply unit, in particular 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 nodes, and the early-stage construction investment cost and the operation and maintenance expenditure cost of the distributed power equipment are reduced by the first objective function, so as to improve the electricity selling income of the distributed power nodes, and maximize the net income of the distributed power nodes.
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 BDA0002633565000000102
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, takes a valueIs 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 historical data of the distributed power nodes and is stored in a preset database in advance, and when needed, the first objective data can be directly obtained 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 BDA0002633565000000111
w in formula (8)L.DNIs expressed as formula (10):
Figure BDA0002633565000000112
w in formula (8)E.DNIs expressed by formula (11):
formula (11)
Figure BDA0002633565000000113
W in formula (8)B1.DNIs formula (12):
formula (12)
Figure BDA0002633565000000114
W in formula (8)B2.DNIs formula (13):
Figure BDA0002633565000000115
w in formula (8)B3.DGIs expressed as formula (14):
Figure BDA0002633565000000116
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; 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; and the uncertainty of the distributed power supply equipment access node influences the operation safety of the incremental power distribution network system more, so that the economy of the incremental power distribution network node optimization planning scheme is reduced, and the loss caused by the output fluctuation of the distributed power supply equipment is expected to be reduced. 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 BDA0002633565000000121
Formula (16)
Figure BDA0002633565000000122
Formula (17)
Figure BDA0002633565000000123
Formula (18)
Figure BDA0002633565000000124
Formula (19)
Figure BDA0002633565000000125
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 BDA0002633565000000131
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.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.
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 BDA0002633565000000132
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, it is desirable 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 shift out loads during the peak electricity price period according to the time-of-use electricity price informationThe load is transferred during the low valley period of the electricity price.
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 BDA0002633565000000141
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. 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 one "max" in each of the first, second and third objective functions, so in the specific calculation in equation (24), W isDG、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 supply 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 capacity of the capacitor.
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 a node 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 formula (25):
equation (25)
Figure BDA0002633565000000161
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.tIndicating that branch ij is transmitting at time tActive electric quantity of, Qij.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 BDA0002633565000000162
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 BDA0002633565000000163
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 BDA0002633565000000164
Formula (29)
Figure BDA0002633565000000165
Formula (30)
Figure BDA0002633565000000166
Formula (31)
Figure BDA0002633565000000167
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 BDA0002633565000000171
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 BDA0002633565000000172
it should be noted that after determining the power capacity of each node of the distributed power supply devices to be accessed, the power 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 capacity of the node of the distributed power supply devices to be accessed and the power 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 capacity of the capacitor. 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 BDA0002633565000000181
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 BDA0002633565000000182
Figure BDA0002633565000000191
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 6:00-8:00&11:00-18:00 period of each 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 sub-distributed power supply nodes 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 nodes 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 on the number of distributed power supply devices accessed to 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 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 capacitance 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 BDA0002633565000000211
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 BDA0002633565000000212
reconstructing the relaxation inequality through a second order cone to obtain the following target inequality:
Figure BDA0002633565000000213
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 BDA0002633565000000214
Figure BDA0002633565000000215
Figure BDA0002633565000000216
Figure BDA0002633565000000221
(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 BDA0002633565000000222
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 representing number of devices accessing distributed power supply at node i，Ni.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 BDA0002633565000000223
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.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 BDA0002633565000000224
Figure BDA0002633565000000231
Figure BDA0002633565000000232
Figure BDA0002633565000000233
Figure BDA0002633565000000234
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 BDA0002633565000000235
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 BDA0002633565000000236
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;
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 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 capacity of the capacitor.
2. The method of configuring incremental distributed power supply equipment according to claim 1, wherein in the 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 FDA0002633564990000021
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 FDA0002633564990000022
reconstructing the relaxation inequality through a second order cone to obtain the following target inequality:
Figure FDA0002633564990000023
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 FDA0002633564990000024
Figure FDA0002633564990000025
Figure FDA0002633564990000026
Figure FDA0002633564990000027
(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 FDA0002633564990000031
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 FDA0002633564990000032
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 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 FDA0002633564990000033
Figure FDA0002633564990000034
Figure FDA0002633564990000041
Figure FDA0002633564990000042
Figure FDA0002633564990000043
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 FDA0002633564990000044
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 FDA0002633564990000045
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 any one of claims 1 to 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 sub-distributed power supply nodes and a first function constraint constraining the first objective function according to the first objective data, build a second objective function corresponding to the incremental power distribution network nodes and a second function constraint constraining 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 constraining the third objective function according to the third objective 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 capacitance 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 FDA0002633564990000061
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 FDA0002633564990000062
reconstructing the relaxation inequality through a second order cone to obtain the following target inequality:
Figure FDA0002633564990000063
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 FDA0002633564990000071
Figure FDA0002633564990000072
Figure FDA0002633564990000073
Figure FDA0002633564990000074
(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 FDA0002633564990000076
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.
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 FDA0002633564990000075
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 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.
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 true CN111952968A (en) 2020-11-17
CN111952968B 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)

Citations (2)

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

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2497177A2 (en) * 2009-11-03 2012-09-12 Spirae, Inc. Dynamic distributed power grid control system
CN103094900A (en) * 2012-12-20 2013-05-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
CN111952968B (en) 2022-05-03

Similar Documents

Publication Publication Date Title
Zhong et al. Cooperative P2P energy trading in active distribution networks: An MILP-based Nash bargaining solution
CN109034563B (en) Multi-subject game incremental power distribution network source-network-load collaborative planning method
Chen et al. Research on day-ahead transactions between multi-microgrid based on cooperative game model
Yang et al. A multi-objective stochastic optimization model for electricity retailers with energy storage system considering uncertainty and demand response
CN111082451B (en) Incremental distribution network multi-objective optimization scheduling model based on scene method
Yuan et al. Real-time pricing for smart grid with multi-energy microgrids and uncertain loads: a bilevel programming method
CN109919450B (en) Game optimization method for solving scheduling of comprehensive intelligent energy system
CN112101607A (en) Active power distribution network rolling optimization scheduling method considering demand response time effect
CN110429653A (en) Consider energy storage and the rural power grids distributed photovoltaic consumption method and terminal device of DR
CN110378729A (en) A kind of integration requirement response method based on dynamic energy price strategy
CN109272353B (en) System dynamic probability energy flow analysis method considering comprehensive demand response uncertainty
Yin et al. Equilibrium stability of asymmetric evolutionary games of multi-agent systems with multiple groups in open electricity market
Sun et al. A coalitional game theoretic energy transaction algorithm for networked microgrids
Zhao et al. A secure intra-regional-inter-regional peer-to-peer electricity trading system for electric vehicles
Zhou et al. A fully-decentralized transactive energy management under distribution network constraints via peer-to-peer energy trading
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
CN115940299B (en) Active and reactive collaborative optimization method for power distribution system source network load storage
CN111952968B (en) Configuration method and device of incremental distributed power supply equipment and readable storage medium
Fan et al. Optimization of power supply capacity of distribution network considering the participation of power sales companies in spot power trading
Zhang et al. Expanding the role of DSOs in establishing flexible local markets: A brokerage perspective
Yan et al. Consensus alternating direction multiplier method based fully distributed peer-to-peer energy transactions considering the network transmission distance
Zhang et al. Coordinated Optimal Scheduling of WPHTNS Power System Based on Adaptive Improved Genetic Algorithm (July 2023)
Chen et al. Stochastic economic dispatch based optimal market clearing strategy considering flexible ramping products under wind power uncertainties
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