CN115411739A - Self-adaptive voltage control method of flexible power distribution system driven by data-physical fusion - Google Patents

Self-adaptive voltage control method of flexible power distribution system driven by data-physical fusion Download PDF

Info

Publication number
CN115411739A
CN115411739A CN202211155987.5A CN202211155987A CN115411739A CN 115411739 A CN115411739 A CN 115411739A CN 202211155987 A CN202211155987 A CN 202211155987A CN 115411739 A CN115411739 A CN 115411739A
Authority
CN
China
Prior art keywords
boundary
measurement
region
representing
distribution system
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.)
Pending
Application number
CN202211155987.5A
Other languages
Chinese (zh)
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.)
Guangdong Power Grid Co Ltd
Original Assignee
Guangdong Power Grid 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 Guangdong Power Grid Co Ltd filed Critical Guangdong Power Grid Co Ltd
Priority to CN202211155987.5A priority Critical patent/CN115411739A/en
Priority to PCT/CN2022/126385 priority patent/WO2024060344A1/en
Publication of CN115411739A publication Critical patent/CN115411739A/en
Pending legal-status Critical Current

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
    • H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • 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/004—Generation forecast, e.g. methods or systems for forecasting future energy generation
    • 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/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • H02J3/16—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by adjustment of reactive power
    • 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)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

A data-physical fusion driven flexible power distribution system adaptive voltage control method comprises the following steps: dividing the selected flexible power distribution system into a complete measurement region and an incomplete measurement region according to the spatial distribution condition of the measurement device, and inputting basic parameter information; initializing a physical guide coefficient, boundary information of a complete measurement region, boundary active power and reactive power requirements, and boundary information of an incomplete measurement region s, and boundary active power and reactive power requirements; calculating an initial value of a pseudo Jacobian matrix of the measured complete region r and an initial value of a distributed power supply operation strategy; establishing a voltage control model of an incomplete measurement area driven by a physical model according to basic parameter information; solving the model; establishing a data-driven self-adaptive voltage control model of a measurement complete region; and solving the self-adaptive voltage control model of the measurement complete region driven by the data. The invention effectively improves the voltage out-of-limit condition of the active power distribution system and completes the complementary support of the inter-area regulation and control resources.

Description

Self-adaptive voltage control method of flexible power distribution system driven by data-physical fusion
Technical Field
The invention relates to a self-adaptive voltage control method for a flexible power distribution system. In particular to a self-adaptive voltage control method of a flexible power distribution system driven by data-physics fusion.
Background
Diversified power supply equipment such as large-scale distributed energy sources are connected into the flexible power distribution system, so that the structural form of the flexible power distribution system evolves, the operation state changes frequently, and higher requirements are provided for the operation control method of the flexible power distribution system. With the improvement of informatization and digitization levels of power distribution systems, the flexible power distribution system accumulates massive multi-source heterogeneous operation data and historical information. Data-driven methods including statistical analysis methods, artificial intelligence methods, iterative learning methods, and the like have been widely used in power distribution system operation optimization problems. The data driving method establishes a data model based on the operation data, establishes an incidence relation between characteristics and research problems, does not depend on accurate physical model parameters, and has the advantages of simple form, high solving speed in a complex scene, strong adaptability and the like.
However, the data-driven optimal control method faces two types of problems: 1) The data model generated by the data driving method has the problems of unclear physical significance and poor interpretability, and can affect the reasonability of results and impact the flexible power distribution system under the condition that the initial value of a pseudo Jacobian matrix or a controlled strategy is unreasonable. 2) The situation of uneven measurement and resource distribution exists in an actual flexible power distribution system, so that a data driving model is difficult to establish due to lack of measurement data in a part of regions, or the regulation and control capability is insufficient due to limited regulation and control resources in a measurement complete region. These two types of problems make practical application of the data-driven control method challenging.
The flexible power distribution system optimization control algorithm based on the physical model establishes a model based on prior information, does not depend on system operation data, and has the advantages of globality, interpretability and the like. When a physical model or an approximate model of a part of flexible power distribution system can be obtained, a data-physical fusion driving mode can be considered to be adopted to construct a fusion driving model so as to overcome various problems of the data driving model and obtain a better application effect. Aiming at the self-adaptive control problem of the flexible power distribution system, the advantages of data and physical models can be combined, and the voltage control model of the flexible power distribution system driven by data-physical fusion is adopted to meet the operation optimization requirements of the flexible power distribution system in different operation scenes.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a flexible power distribution system self-adaptive voltage control method driven by data-physical fusion, which can obtain a distributed power supply operation control strategy under the conditions of lack of accurate parameters and uneven space distribution of a measuring device.
The technical scheme adopted by the invention is as follows: a data-physical fusion driven flexible power distribution system adaptive voltage control method comprises the following steps:
1) According to the space distribution condition of the measuring device, dividing the selected flexible power distribution system into a complete measuring region r and an incomplete measuring region s, and inputting basic parameter information, wherein the basic parameter information comprises the following steps: line parameter information, load, distributed power supply position and parameter, node voltage reference value, and setting start time to be t = t 0 + delta T, optimizing the total control time length to be T, and controlling the step length delta T;
2) Initializing a physical boot coefficient σ t 0 ]Measuring boundary information omega of the complete region r r [t 0 ]And boundary active and reactive power requirements
Figure BDA0003858634510000021
And
Figure BDA0003858634510000022
boundary information omega of incomplete measurement area s s [t 0 ]And boundary active and reactive power requirements
Figure BDA0003858634510000023
And
Figure BDA0003858634510000024
calculating the initial value of the pseudo Jacobian matrix of the measured complete region r based on the physical model according to the basic parameter information input in the step 1)
Figure BDA0003858634510000025
Initial value x of distributed power supply operation strategy p [t 0 ]；
3) The incomplete region s at the time of measuring t-delta t receives the boundary information omega of the adjacent complete region r r [t-Δt]And boundary active and reactive power requirements
Figure BDA0003858634510000026
And
Figure BDA0003858634510000027
according to the basic parameter information, establishing a physical model-driven measuring incomplete area s voltage control model, which comprises the following steps: the method comprises the following steps of taking the minimum deviation of the node voltage of the region, the minimum deviation of the boundary information of adjacent regions and the minimum deviation of the boundary power requirement as objective functions, and considering the operation constraint of a flexible power distribution system and the capacity constraint of a distributed power supply;
4) Solving the voltage control model of the incomplete measurement area s driven by the physical model to obtain a distributed power supply operation control strategy in the incomplete measurement area s, issuing and executing the strategy, and sending the boundary information omega s [t]And boundary active and reactive power requirements
Figure BDA0003858634510000028
Figure BDA0003858634510000029
Transmitting to the measurement completion region r;
5) The measurement complete region r obtains voltage measurement information of each measurement node, and receives boundary information omega of adjacent measurement incomplete regions s s [t-Δt]And boundary active and reactive power requirements
Figure BDA00038586345100000210
And
Figure BDA00038586345100000211
updating the physical guiding coefficient σ t]Establishing a data-driven self-adaptive voltage control model of the measurement completion region r, comprising the following steps: taking the minimum voltage deviation of the local area measurement node, the minimum boundary information deviation of adjacent areas and the minimum boundary power demand deviation as objective functions, and considering the capacity constraint of the distributed power supply;
6) Solving the data-driven measurement complete region r self-adaptive voltage control model in the step 5), obtaining a distributed power supply operation control strategy in the measurement complete region r, issuing and executing the strategy, and sending boundary information omega r [t]And boundary active and reactive power requirements
Figure BDA00038586345100000212
Transmitting to the incomplete measurement area s;
7) And updating the control time T = T + Δ T, judging whether T is greater than the optimization duration T, if not, turning to the step 3), and if so, ending.
The self-adaptive voltage control method of the data-physical fusion driven flexible power distribution system comprehensively considers the unknown line parameters of the flexible power distribution system, the uncertainty of the output condition of the distributed power supply and the uneven spatial distribution condition of the measuring device, and realizes the self-adaptive voltage control of the data-physical fusion driven flexible power distribution system by establishing a self-adaptive predictive voltage control strategy of the data-driven distributed energy storage system. A new idea is provided for the voltage optimization problem of the power distribution network, and the safety of the power distribution side and the user experience are improved.
The invention considers the adoption of a data-physical fusion driving mode to construct a fusion driving model, effectively improves the rationality and the interpretability of a data driving control effect, relieves the impact on a power distribution system caused by overlarge change of a control strategy, realizes the mutual assistance of inter-region target functions, fully calls the voltage regulation potential of a distributed power converter, effectively improves the voltage out-of-limit condition of an active power distribution system, completes the complementary support of inter-region regulation and control resources and improves the flexible and efficient operation level of the active power distribution system.
Drawings
FIG. 1 is a flow chart of the adaptive voltage control method of the flexible power distribution system driven by data-physical fusion according to the present invention;
FIG. 2 is a flexible power distribution system topology as employed in the present invention;
FIG. 3 is a predicted curve of distributed power output and load change;
FIG. 4 is a comparison of the voltages at node 18;
FIG. 5 is a distributed power reactive power take off strategy at node 18;
FIG. 6 is a scene one and a scene two 24 hour voltage control effect;
fig. 7 is a scene two and a scene three 24 hours voltage control effect.
Detailed Description
The data-physical fusion driven flexible power distribution system adaptive voltage control method of the invention is described in detail below with reference to the embodiments and the accompanying drawings.
As shown in fig. 1, the method for controlling adaptive voltage of a flexible power distribution system driven by data-physical fusion of the present invention comprises the following steps:
1) According to the space distribution condition of the measuring device, dividing the selected flexible power distribution system into a complete measuring region r and an incomplete measuring region s, and inputting basic parameter information, wherein the basic parameter information comprises the following steps: line parameter information, load, distributed power supply position and parameter, node voltage reference value, and setting start time to be t = t 0 + delta T, optimizing the total control time length to be T, and controlling the step length delta T;
2) Initializing the physical guiding coefficient σ t 0 ]Measuring boundary information omega of the complete region r r [t 0 ]And boundary active and reactive power requirements
Figure BDA0003858634510000031
And
Figure BDA0003858634510000032
measuring boundary information omega of incomplete area s s [t 0 ]And boundary active and reactive power requirements
Figure BDA0003858634510000033
And
Figure BDA0003858634510000034
calculating an initial value of a pseudo Jacobian matrix of the measured complete region r based on a physical model according to the basic parameter information input in the step 1)
Figure BDA0003858634510000035
Initial value x of distributed power supply operation strategy p [t 0 ](ii) a Wherein the content of the first and second substances,
the initial value of the pseudo Jacobian matrix of the measured complete region r is calculated based on the physical model
Figure BDA0003858634510000036
The following were used:
Figure BDA0003858634510000037
in the formula, x p Representing the controlled equipment strategy of the flexible power distribution system calculated by a physical model;
Figure BDA0003858634510000038
representing typical scenes, representing a typical scene set of the flexible power distribution system, N Λ Representing a typical number of scenes generated from historical data;
said baseCalculating and measuring initial value x of distributed power supply operation strategy in complete region r in physical model p [t 0 ]The following are:
x p [t 0 ]＝argmin(f) (2)
in the formula, x p [t 0 ]Expressing a strategy initial value of controlled equipment of the flexible power distribution system, considering the operation constraint of the flexible power distribution system, and solving through a flexible power distribution system operation optimization control model based on a physical model:
g＝min(x p ，y)
Figure BDA0003858634510000039
in the formula, x p Represents a control variable; y represents a flow variable matrix; g represents an objective function; m (x) p Y) represents inequality constraints in the flexible power distribution system operating constraints, including system safety constraints, distributed power supply operating constraints; n (x) p Y) represents an equality constraint in the flexible power distribution system operating constraints, including a power flow constraint; in the formula x max And x min For controlling the upper and lower limits of the variables, y max And y min The upper and lower limits of the power flow variable are shown.
3) The incomplete region s at the time of measuring t-delta t receives the boundary information omega of the adjacent complete region r r [t-Δt]And boundary active and reactive power requirements
Figure BDA00038586345100000310
And
Figure BDA00038586345100000311
according to the basic parameter information, establishing a physical model-driven measuring incomplete area s voltage control model, which comprises the following steps: the method comprises the following steps of taking the minimum deviation of the node voltage of the region, the minimum deviation of the boundary information of adjacent regions and the minimum deviation of the boundary power requirement as objective functions, and considering the operation constraint of a flexible power distribution system and the capacity constraint of a distributed power supply; wherein the content of the first and second substances,
the boundary information omega r [t-Δt]The method comprises the following steps of boundary link transmission power, boundary link current information and boundary voltage information of an adjacent measurement complete region r, wherein the boundary link transmission power, the boundary link current information and the boundary voltage information are expressed as follows:
Figure BDA0003858634510000041
in the formula (I), the compound is shown in the specification,
Figure BDA0003858634510000042
and
Figure BDA0003858634510000043
respectively representing the measurement values of active power and reactive power on a boundary transmission tie line l of the measurement completion region r at the time of t-delta t;
Figure BDA0003858634510000044
the voltage measurement value of the boundary node n of the complete measurement region r at the time of t-delta t is represented;
Figure BDA0003858634510000045
the current measurement value on the boundary link l of the measured complete region r at the time point t- Δ t is shown.
Considering the operation constraint of the flexible power distribution system, the s voltage control model of the incomplete measurement area driven by the physical model is expressed as follows:
Figure BDA0003858634510000046
Figure BDA0003858634510000047
wherein f represents an objective function; f. of 1 Representing the voltage deviation of a node in an s area of an incomplete measurement area; f. of 2 An iterative error penalty item representing the boundary value of the incomplete measurement region s and the adjacent region; f. of 3 Representing the deviation of the boundary power of the measured incomplete area s and the boundary requirement of the adjacent area; x is the number of s Indicates that the measurement is incompleteA controlled equipment operation strategy of the standby area s; x' s Representing the operation strategy of the controlled object of the measurement incomplete area s;
Figure BDA0003858634510000048
and
Figure BDA0003858634510000049
respectively representing boundary active power and reactive power requirements sent to adjacent regions by the complete measurement region r at the time t as boundary information; n is a radical of s Measuring the number of nodes in the incomplete area s; v. of s，i [t]Representing the voltage value of the s node i of the incomplete measurement area at the time t;
Figure BDA00038586345100000410
a node voltage reference value representing a measurement incomplete area s; y represents a flow variable matrix; omega s [t]Representing boundary information obtained by calculating through an s physical model of a measured incomplete area at the time t; omega s Representing an internal node set of the incomplete measurement area s;
Figure BDA00038586345100000411
and
Figure BDA00038586345100000412
respectively representing auxiliary variables of the incomplete measurement area s at the time t and the time t-delta t; omega s [t-Δt]Represents boundary state information m (x ') obtained by calculating a measured incomplete area s at the time of t-delta t' s Y) represents inequality constraints in the flexible power distribution system operation constraints, including system safety constraints and distributed power supply operation constraints; n (x' s And y) represents an equality constraint in the flexible power distribution system operating constraints, including power flow constraints. (given an explanation of all letters)
4) Solving the voltage control model of the incomplete measurement area s driven by the physical model to obtain a distributed power supply operation control strategy in the incomplete measurement area s, issuing and executing the strategy, and sending the boundary information omega s [t]And boundary active and reactive power requirements
Figure BDA0003858634510000051
Figure BDA0003858634510000052
Transmitting to the measurement completion region r;
the boundary information omega s [t]The method comprises the following steps of calculating the transmission power of the boundary connecting line, the current information of the boundary connecting line and the boundary voltage information of the incomplete measurement area s, and expressing the transmission power, the current information and the boundary voltage information of the boundary connecting line as follows:
Figure BDA0003858634510000053
in the formula (I), the compound is shown in the specification,
Figure BDA0003858634510000054
and
Figure BDA0003858634510000055
respectively representing boundary active power and reactive power requirements on a boundary connecting line l of an incomplete measurement region s at the moment t;
Figure BDA0003858634510000056
boundary voltage information of an s node n of an incomplete measurement area at the moment t is represented;
Figure BDA0003858634510000057
the current measurement value on the boundary link l of the complete region s is measured at time t.
5) The measurement complete region r obtains voltage measurement information of each measurement node, and receives boundary information omega of adjacent measurement incomplete regions s s [t-Δt]And boundary active and reactive power requirements
Figure BDA0003858634510000058
And
Figure BDA0003858634510000059
updating the physical guiding coefficient σ t]The method for establishing the data-driven measurement complete region r self-adaptive voltage control model comprises the following steps: to be provided withThe minimum voltage deviation of the regional measurement node, the minimum boundary information deviation of adjacent regions and the minimum boundary power demand deviation are objective functions, and distributed power supply capacity constraint is considered; wherein, the first and the second end of the pipe are connected with each other,
the updated physical guiding coefficient sigma [ t ] is as follows:
Figure BDA00038586345100000510
in the formula, σ [ t []And σ [ t- Δ t]Representing the physical guiding coefficient at time t and at-time t, Y r [t]The measured value of the r node voltage of the complete measurement region at the time t is shown, and gamma represents the voltage guiding range.
The self-adaptive voltage control model of the data-driven measurement complete region r is expressed as follows:
Figure BDA00038586345100000511
wherein J represents an objective function, J 1 Indicating the sum of the node voltage deviation and the variation amplitude of the controlled variable in the r region of the measurement completion region; j. the design is a square 2 An iterative error penalty term representing the boundary value between the measured complete region r and the adjacent region; j. the design is a square 3 Representing the deviation between the boundary power of the measured complete region r and the boundary power requirement of the adjacent region;
Figure BDA00038586345100000512
indicating the reference value of r node voltage in the measurement completion region; x r [t]Measuring a controlled equipment operation strategy of the complete region r for the time t; y is r [t]Representing a voltage measurement value of an r node of a complete measurement region at the time t; x' r [t]Representing the operation strategy of the controlled object of the complete measurement region r at the time t; delta X' r [t]＝X′ r [t]-X′ r [t-Δt]Representing the operation strategy variation of the controlled object of the complete measurement region r at the time t; delta X' r [t-Δt]＝X′ r [t-Δt]-X′ r [t-2Δt]Representing the operation strategy variation of the controlled object in the complete measurement region r at the time of t-delta t;
Figure BDA0003858634510000061
and two
Figure BDA0003858634510000062
Respectively representing the active power and reactive power requirements sent to adjacent regions by the complete measurement region r at the time t as boundary information;
Figure BDA0003858634510000063
and
Figure BDA0003858634510000064
auxiliary variables respectively representing the measurement completion region r at the time t and the time t-delta t;
Figure BDA0003858634510000065
the estimation value of the pseudo Jacobian matrix of the measured complete region r at the time t is represented by the following method:
Figure BDA0003858634510000066
Figure BDA0003858634510000067
in the formula, delta X' r [t-Δt]＝X′ r [t-Δt]-X′ r [t-2Δt]Representing the operation strategy variation of the controlled equipment at the time t-delta t; e represents a parameter; η and μ represent weight coefficients;
Figure BDA0003858634510000068
representing an initial value of a pseudo Jacobian matrix; (explanation of all letters is given)
The consideration of the distributed power capacity constraint is expressed as follows:
Figure BDA0003858634510000069
in formula (II), X' r [t]Measuring a controlled equipment operation strategy of the complete region r for the time t; pr [ t ]]Representing the active power output, C, of the distributed power supply in the complete measurement region r at the moment t r [t]And the capacity of the distributed power supply in the measurement completion region r at the time t is shown.
6) Solving the data-driven measurement complete region r self-adaptive voltage control model in the step 5), obtaining a distributed power supply operation control strategy in the measurement complete region r, issuing and executing the strategy, and sending boundary information omega r [t]And boundary active and reactive power requirements
Figure BDA00038586345100000610
Transmitting to the incomplete measurement area s;
7) And updating the control time T = T + delta T, judging whether T is greater than the optimization time length T, if not, turning to the step 3), and if yes, ending.
Specific examples are given below:
for the embodiment of the invention, the power distribution network comprises 33 nodes, the topological connection situation is as shown in fig. 2, a fan system with the capacity of 500kVA is accessed at the system nodes 13, 15, 16, 29 and 30, and a photovoltaic system with the capacity of 100kWp is accessed at the nodes 11, 12, 17, 18, 20, 21, 22, 23, 24, 25, 32 and 33. The distributed power output fluctuation curves are shown in fig. 3-4. The system voltage is 12.66kV, the reference power is 1MVA, and the active load and the reactive load of the system are 3715kW and 2300kvar respectively. Controlling the step length delta T =0.5min, and optimizing the time length T =24h; the voltage reference value of the distribution network was set to 1.0p.u. The values of the weight coefficients λ, ρ, η, μ, and δ are all 1.0, and the voltage guiding range γ = [0.97p.u.,1.03p.u ]). And optimizing by adopting a data-physical fusion driven voltage control method of the flexible power distribution system, and obtaining a distributed power output strategy through the steps. In order to verify the effectiveness of the method, the following three control scenes are adopted for the flexible power distribution system for comparison:
scene one: the distributed energy storage system is not controlled, and the initial running state of the flexible power distribution system is obtained;
scene two: a voltage control method of a flexible power distribution system driven by data-physical fusion is adopted.
Scene three: a data-driven voltage control method for a flexible power distribution system is adopted.
The computer hardware environment for executing the optimized calculation is Intel (R) Core (TM) CPU i5-10210U, the main frequency is 1.6GHz, and the memory is 16GB; the software environment is a Windows10 operating system.
An example topology used by an embodiment of the present invention is shown in FIG. 2. The effect pairs of the all-day optimization control under the three scenes are shown in table 1. The variation of the prediction curve of the distributed power output and load information is shown in fig. 3. The voltage at node 18 is plotted as shown in fig. 4. The distributed power reactive power out strategy at node 18 is shown in fig. 5. The effect of scene one and scene two 24 hour voltage control is shown in figure 6. The effects of scene two and scene three 24 hour voltage control are shown in figure 7.
As can be seen from table 1 and fig. 4 to 7, the second scenario can effectively adjust the voltage level of the flexible power distribution system in this embodiment, and effectively reduce the global voltage deviation, and the voltage deviation index is reduced by 64.13% compared with the first scenario. Compared with a data driving method in a third scene, the data-physical fusion driving method in the second scene can improve the global optimization effect, and the voltage deviation index of the data-physical fusion driving method is reduced by 43.10% compared with that of the third scene.
TABLE 1 comparison of the effects of the optimization control throughout the day
Figure BDA0003858634510000071

Claims (8)

1. A self-adaptive voltage control method of a data-physical fusion driven flexible power distribution system is characterized by comprising the following steps:
1) According to the space distribution condition of the measuring device, dividing the selected flexible power distribution system into a complete measuring region r and an incomplete measuring region s, and inputting basic parameter information, wherein the basic parameter information comprises the following steps: line parameter information, load, distributed power supply position and parameter, node voltage reference value, and setting start time to be t = t 0 + delta T, optimizing the total control time length to be T, and controlling the step length delta T;
2) Initializing a physical boot coefficient σ t 0 ]Measuring boundary information omega of the complete region r r [t 0 ]And boundary active and reactive power requirements
Figure FDA0003858634500000011
And
Figure FDA0003858634500000012
boundary information omega of incomplete measurement area s s [t 0 ]And boundary active and reactive power requirements
Figure FDA0003858634500000013
And
Figure FDA0003858634500000014
calculating the initial value of the pseudo Jacobian matrix of the measured complete region r based on the physical model according to the basic parameter information input in the step 1)
Figure FDA0003858634500000015
Initial value x of distributed power supply operation strategy p [t 0 ]；
3) The incomplete region s at the time of measuring t-delta t receives the boundary information omega of the adjacent complete region r r [t-Δt]And boundary active and reactive power requirements
Figure FDA0003858634500000016
And
Figure FDA0003858634500000017
according to the basic parameter information, establishing a physical model-driven measuring incomplete area s voltage control model, which comprises the following steps: the method comprises the following steps of taking the minimum deviation of the node voltage of the region, the minimum deviation of the boundary information of adjacent regions and the minimum deviation of the boundary power requirement as objective functions, and considering the operation constraint of a flexible power distribution system and the capacity constraint of a distributed power supply;
4) Driven by solving said physical modelMeasuring the voltage control model of the incomplete area s, obtaining the distributed power supply operation control strategy in the incomplete area s, issuing and executing the strategy, and sending the boundary information omega s [t]And boundary active and reactive power requirements
Figure FDA0003858634500000018
Figure FDA0003858634500000019
Transmitting to the measurement completion region r;
5) The measurement complete region r obtains voltage measurement information of each measurement node, and receives boundary information omega of adjacent measurement incomplete regions s s [t-Δt]And boundary active and reactive power requirements
Figure FDA00038586345000000110
And
Figure FDA00038586345000000111
updating the physical guiding coefficient σ t]Establishing a data-driven self-adaptive voltage control model of the measurement completion region r, comprising the following steps: taking the minimum voltage deviation of the measuring node in the region, the minimum boundary information deviation of adjacent regions and the minimum boundary power demand deviation as objective functions, and considering the capacity constraint of the distributed power supply;
6) Solving the data-driven measurement complete region r self-adaptive voltage control model in the step 5), obtaining a distributed power supply operation control strategy in the measurement complete region r, issuing and executing the strategy, and sending boundary information omega r [t]And boundary active and reactive power requirements
Figure FDA00038586345000000112
Figure FDA00038586345000000113
Transmitting to the incomplete measurement area s;
7) And updating the control time T = T + delta T, judging whether T is greater than the optimization time length T, if not, turning to the step 3), and if yes, ending.
2. The adaptive voltage control method of the data-physical fusion driven flexible power distribution system according to claim 1, wherein in step 2),
the initial value of the pseudo Jacobian matrix of the measured complete region r is calculated based on the physical model
Figure FDA00038586345000000114
The following were used:
Figure FDA00038586345000000115
in the formula, x p Representing the controlled equipment strategy of the flexible power distribution system calculated by a physical model; zeta represents a typical scenario, lambda represents a typical scenario set of the flexible power distribution system, N Λ Representing a typical number of scenes generated from historical data;
the initial value x of the distributed power supply operation strategy of the measurement complete region r is calculated based on the physical model p [t 0 ]The following are:
x p [t 0 ]＝argmin(f) (2)
in the formula, x p [t 0 ]Expressing a strategy initial value of controlled equipment of the flexible power distribution system, considering the operation constraint of the flexible power distribution system, and solving through a flexible power distribution system operation optimization control model based on a physical model:
g＝min(x p ,y)
Figure FDA0003858634500000021
in the formula, x p Represents a control variable; y represents a flow variable matrix; g represents an objective function; m (x) p Y) represents inequality constraints in the flexible power distribution system operation constraints, including system safety constraints and distributed power supply operation constraints; n (x) p And y) represents a flexible power distribution systemAn equality constraint in the operational constraint, including a power flow constraint; in the formula x max And x min For upper and lower limits of the control variable, y max And y min The upper and lower limits of the flow variable are shown.
3. The adaptive voltage control method for the flexible power distribution system driven by data-physics fusion as claimed in claim 1, wherein the boundary information ω in step 3) is r [t-Δt]The method comprises the following steps of boundary link transmission power, boundary link current information and boundary voltage information of an adjacent measurement complete region r, wherein the boundary link transmission power, the boundary link current information and the boundary voltage information are expressed as follows:
Figure FDA0003858634500000022
in the formula (I), the compound is shown in the specification,
Figure FDA0003858634500000023
and
Figure FDA0003858634500000024
respectively representing the measurement values of active power and reactive power on a boundary transmission tie line l of the measurement completion region r at the time of t-delta t;
Figure FDA0003858634500000025
the voltage measurement value of the boundary node n of the complete measurement region r at the time of t-delta t is represented;
Figure FDA0003858634500000026
the current measurement value on the boundary link l of the measured completion region r at the time point t- Δ t is shown.
4. The adaptive voltage control method of the data-physical fusion driven flexible power distribution system according to claim 1, wherein the flexible power distribution system operation constraint is considered, and the physical model driven measurement incomplete region s voltage control model in step 3) is represented as:
Figure FDA0003858634500000027
Figure FDA0003858634500000031
wherein f represents an objective function; f. of 1 Representing the voltage deviation of a node in an s area of an incomplete measurement area; f. of 2 An iterative error penalty item representing the boundary value of the incomplete measurement region s and the adjacent region; f. of 3 Representing the deviation of the boundary power of the measured incomplete area s and the boundary requirement of the adjacent area; x is the number of s Representing a controlled equipment operation strategy of the measurement incomplete area s; x' s Representing the operation strategy of the controlled object of the measurement incomplete area s;
Figure FDA0003858634500000032
and
Figure FDA0003858634500000033
respectively representing boundary active power and reactive power requirements sent to adjacent regions by the complete measurement region r at the time t as boundary information; n is a radical of s Measuring the number of nodes in the incomplete area s; v. of s,i [t]Representing the voltage value of the s node i of the incomplete measurement area at the time t;
Figure FDA0003858634500000034
a node voltage reference value representing a measurement incomplete area s; y represents a power flow variable matrix; omega s [t]Representing boundary information obtained by calculating through an s physical model of a measured incomplete area at the time t; omega s Representing an internal node set of the incomplete measurement area s;
Figure FDA0003858634500000035
and
Figure FDA0003858634500000036
respectively representing auxiliary variables of the incomplete measurement area s at the time t and the time t-delta t; omega s [t-Δt]And the boundary state information obtained by calculating the incomplete area s measured at the time of t-delta t is represented. m (x' s Y) represents inequality constraints in the flexible power distribution system operating constraints, including system safety constraints, distributed power supply operating constraints; n (x' s And y) represents an equality constraint in the flexible power distribution system operating constraints, including power flow constraints.
5. The adaptive voltage control method of the data-physical fusion driven flexible power distribution system according to claim 1, wherein the boundary information ω in step 4) is s [t]The method comprises the following steps of calculating the transmission power of the boundary connecting wire, the current information of the boundary connecting wire and the boundary voltage information of the boundary connecting wire in an incomplete measurement area s, and expressing the transmission power, the current information of the boundary connecting wire and the boundary voltage information as follows:
Figure FDA0003858634500000037
in the formula (I), the compound is shown in the specification,
Figure FDA0003858634500000038
and
Figure FDA0003858634500000039
respectively representing the boundary active power and reactive power requirements on the boundary connecting line l of the incomplete area s measured at the time t;
Figure FDA00038586345000000310
representing boundary voltage information of an s node n of an incomplete measurement area at the time t;
Figure FDA00038586345000000311
the current measurement value on the boundary link l of the complete region s is measured at time t.
6. The adaptive voltage control method of the data-physical fusion driven flexible power distribution system according to claim 1, wherein the updated physical guidance coefficient σ [ t ] in step 5) is:
Figure FDA00038586345000000312
in the formula, σ [ t []And σ [ t- Δ t]Representing the physical guiding coefficient, Y, at times t and t- Δ t, respectively r [t]The voltage measurement value of the r node of the complete measurement region at the time t is shown, and gamma represents the voltage guiding range.
7. The adaptive voltage control method of the data-physical fusion driven flexible power distribution system according to claim 1, wherein the adaptive voltage control model of the data-driven measurement completion region r in step 5) is represented as:
Figure FDA00038586345000000313
Figure FDA0003858634500000041
wherein J represents an objective function, J 1 Indicating the sum of the node voltage deviation and the variation amplitude of the controlled variable in the r region of the measurement completion region; j. the design is a square 2 An iterative error penalty term representing the boundary value between the measured complete region r and the adjacent region; j. the design is a square 3 Representing the deviation between the boundary power of the measured complete region r and the boundary power requirement of the adjacent region;
Figure FDA0003858634500000042
indicating a voltage reference value of an r node of a measurement completion region; x r [t]Measuring the operation strategy of the controlled equipment in the complete region r for the time t; y is r [t]Representing a voltage measurement value of an r node of a complete measurement region at the time t; x' r [t]Representing the operation strategy of the controlled object of the complete measurement region r at the time t; delta X' r [t]＝X′ r [t]-X′ r [t-Δt]Representing the operation strategy variation of the controlled object in the complete measurement region r at the time t; Δ X' r [t-Δt]＝X′ r [t-Δt]-X′ r [t-2Δt]Representing the operation strategy variation of the controlled object in the complete measurement region r at the time of t-delta t;
Figure FDA0003858634500000043
and
Figure FDA0003858634500000044
respectively representing the active power and reactive power requirements sent to adjacent regions by the complete measurement region r at the time t as boundary information;
Figure FDA0003858634500000045
and
Figure FDA0003858634500000046
auxiliary variables respectively representing the measurement completion region r at the time t and the time t-delta t;
Figure FDA0003858634500000047
the estimation value of the pseudo Jacobian matrix of the measured complete region r at the time t is represented by the following method:
Figure FDA0003858634500000048
Figure FDA0003858634500000049
in the formula, delta X' r [t-Δt]＝X′ r [t-Δt]-X′ r [t-2Δt]Representing the operation strategy variation of the controlled object at the time t-delta t; e represents a parameter; η and μ represent weight coefficients;
Figure FDA00038586345000000410
representing a pseudo JacobianAnd (5) initial values of the matrix.
8. The adaptive voltage control method of the data-physical fusion driven flexible power distribution system according to claim 1, wherein the consideration of the distributed power capacity constraint in step 5) is represented as follows:
Figure FDA00038586345000000411
in the formula, X r [t]Measuring a controlled equipment operation strategy of the complete region r for the time t; p is r [t]Representing the active power output, C, of the distributed power supply in the complete measurement region r at the moment t r [t]And the capacity of the distributed power supply in the measurement completion region r at the time t is shown.
CN202211155987.5A 2022-09-22 2022-09-22 Self-adaptive voltage control method of flexible power distribution system driven by data-physical fusion Pending CN115411739A (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202211155987.5A CN115411739A (en) 2022-09-22 2022-09-22 Self-adaptive voltage control method of flexible power distribution system driven by data-physical fusion
PCT/CN2022/126385 WO2024060344A1 (en) 2022-09-22 2022-10-20 Data-physics fusion-driven adaptive voltage control system for flexible power distribution system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211155987.5A CN115411739A (en) 2022-09-22 2022-09-22 Self-adaptive voltage control method of flexible power distribution system driven by data-physical fusion

Publications (1)

Publication Number Publication Date
CN115411739A true CN115411739A (en) 2022-11-29

Family

ID=84164933

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211155987.5A Pending CN115411739A (en) 2022-09-22 2022-09-22 Self-adaptive voltage control method of flexible power distribution system driven by data-physical fusion

Country Status (2)

Country Link
CN (1) CN115411739A (en)
WO (1) WO2024060344A1 (en)

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105226653B (en) * 2015-11-06 2017-12-26 贵州电网公司信息通信分公司 The exact linearization method of transformer model in a kind of active distribution network optimal load flow
CN106972539B (en) * 2017-05-13 2019-05-14 天津大学 A kind of distributed generation resource voltage control strategy setting method on the spot based on cone planning
CN111682594B (en) * 2020-06-15 2023-02-21 天津大学 Data-driven model-free adaptive voltage control method for flexible substation of power distribution network
CN113890016B (en) * 2021-09-26 2022-09-23 天津大学 Data-driven multi-time scale voltage coordination control method for power distribution network
CN114884063A (en) * 2022-04-11 2022-08-09 国网天津市电力公司 Distributed power supply and energy storage voltage control method and device considering interval coordination

Also Published As

Publication number Publication date
WO2024060344A1 (en) 2024-03-28

Similar Documents

Publication Publication Date Title
CN114142498B (en) Data-driven distributed energy storage self-adaptive prediction control voltage regulation method
CN109412202B (en) Distributed power supply coordination control method and system
CN113900380B (en) Robust output formation tracking control method and system for heterogeneous cluster system
CN111682594B (en) Data-driven model-free adaptive voltage control method for flexible substation of power distribution network
US10804702B2 (en) Self-organizing demand-response system
CN110858718B (en) Alternating current micro-grid distributed event-driven frequency control method considering economy
CN113315164B (en) Reactive voltage control method and device, medium and computing device
CN115345380A (en) New energy consumption electric power scheduling method based on artificial intelligence
CN114884063A (en) Distributed power supply and energy storage voltage control method and device considering interval coordination
CN113890017B (en) Power distribution network voltage self-adaptive control method based on key measurement
CN111969662B (en) Data-driven multi-intelligent soft switch partition cooperative adaptive voltage control method
CN114722693A (en) Optimization method of two-type fuzzy control parameter of water turbine regulating system
CN112952831B (en) Daily optimization operation strategy for providing stacking service by load side energy storage
CN114566971A (en) Real-time optimal power flow calculation method based on near-end strategy optimization algorithm
Mu et al. Adaptive composite frequency control of power systems using reinforcement learning
CN114444256A (en) Virtual power plant load prediction method and tracking control method based on big data
CN115411739A (en) Self-adaptive voltage control method of flexible power distribution system driven by data-physical fusion
CN115001002B (en) Optimal scheduling method and system for solving problem of energy storage participation peak clipping and valley filling
CN116976523A (en) Distributed economic dispatching method based on partially observable reinforcement learning
CN113269420B (en) Distributed event-driven power economy scheduling method based on communication noise
CN115360768A (en) Power scheduling method and device based on muzero and deep reinforcement learning and storage medium
CN115659792A (en) Multi-period multi-scene SCUC decoupling method, system, equipment and storage medium
CN115204056A (en) Power grid prospective optimization scheduling method and system
US20220373989A1 (en) System for configuring demand response for energy grid assets
CN111525547A (en) Low-voltage intelligent management method based on optimal reactive compensation

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