Distributed power supply locating and sizing method applied to expressway energy system

Abstract

The invention discloses a distributed power supply locating and sizing method applied to an expressway energy system, and belongs to the field of distributed power supply locating and sizing. The method comprises the following steps: establishing output constraint of a distributed power supply, and processing historical data of the distributed power supply according to the output constraint of the distributed power supply to obtain a plurality of output scenes; clustering each output scene by using a clustering algorithm to obtain a typical output scene; according to the typical output scene, the construction position and capacity of the distributed power supply are obtained by taking the minimum cost of the highway energy system under the normal condition and the minimum load loss of the highway energy system under the fault condition as targets, and locating and sizing of the distributed power supply are completed. The method solves the problem that the prior art cannot give consideration to normal and fault operation states of an expressway energy system to perform wind and light power station locating and sizing.

Description

Technical Field

[0001] The present invention belongs to the field of distributed power source site selection and capacity determination, and in particular relates to a distributed power source site selection and capacity determination method applied to a highway energy system. Background Art

[0002] The deployment of wind and solar power stations before disasters can effectively improve the resilience of highway energy systems to cope with extreme disasters. However, few existing technologies consider both normal and faulty operating states of highway energy systems when selecting wind and solar power station sites and determining their capacity. Summary of the invention

[0003] In view of the above-mentioned deficiencies in the prior art, the present invention provides a distributed power source site selection and sizing method applied to a highway energy system, which solves the problem that the prior art cannot take into account both normal and fault operating states of the highway energy system when selecting the site and sizing the wind and solar power stations.

[0004] In order to achieve the above-mentioned invention object, the technical solution adopted by the present invention is: a distributed power source site selection and capacity determination method applied to a highway energy system, comprising the following steps:

[0005] Establish output constraints of distributed power sources, and process historical data of distributed power sources according to the output constraints of distributed power sources to obtain several output scenarios;

[0006] Clustering algorithms are used to cluster various output scenarios to obtain typical output scenarios;

[0007] According to typical output scenarios, with the goal of minimizing the cost of the highway energy system under normal circumstances and minimizing the load loss of the highway energy system under fault conditions, the construction location and capacity of the distributed power source are obtained, and the site selection and sizing of the distributed power source are completed.

[0008] The beneficial effects of the present invention are: by rationally planning the capacity and location of distributed power sources, the economy of the highway energy system is improved; by configuring wind and solar power stations before a disaster, the resilience of the highway energy system is effectively improved.

[0009] Furthermore, the expression of the output constraint of the distributed power source is:

[0010]

[0011]

[0012]

[0013]

[0014] T C (t) = Ttemp (t)+(N oct -20)I(t) / I stc

[0015]

[0016]

[0017] in, is the maximum output power of wind farm i1 at time t; v(t) is the wind speed at the hub of the wind turbine generator at time t; v in is the cut-in wind speed of the wind turbine; v out is the cut-out wind speed of the wind turbine; is the rated power of wind farm i1; v r is the rated wind speed of the wind turbine; is the actual output active power of wind farm i1 at time t; is the minimum output reactive power of wind farm i1; is the actual output reactive power of wind farm i1 at time t; is the maximum output reactive power of wind farm i1; is the maximum output power of photovoltaic power station i2 at time t; I(t) is the light intensity at time t; φ is the derating coefficient of photovoltaic cell; is the rated power of the photovoltaic power station i2; β is the power temperature coefficient; T C (t) is the actual operating temperature of the photovoltaic cell at time t; T stc (t) is the standard operating temperature of the photovoltaic cell at time t; stc is the light intensity under standard test conditions; T temp (t) is the ambient temperature at time t; N oct is the nominal temperature of the photovoltaic cell; is the actual output active power of photovoltaic power station i2 at time t; is the minimum output reactive power of photovoltaic power station i2; is the actual output reactive power of photovoltaic power station i2 at time t; is the maximum output reactive power of photovoltaic power station i2.

[0018] The beneficial effect of the above further scheme is that expressing the wind turbine output constraints and photovoltaic output constraints in the form of a mathematical model can significantly improve the efficiency of calculating wind and solar output scenarios.

[0019] Furthermore, the typical output scenario includes a typical output scenario under normal conditions and a typical output scenario under fault conditions.

[0020] The beneficial effect of the above further scheme is that it targets the output scenarios of the highway energy system under different conditions and is more comprehensive in analysis.

[0021] Furthermore, according to the typical output scenario, the location and capacity determination strategy of distributed power sources is obtained with the goal of minimizing the cost of the highway energy system under normal conditions and minimizing the load loss of the highway energy system under fault conditions, which is specifically:

[0022] A multi-objective optimization model is constructed with the goal of minimizing the cost of the highway energy system under normal circumstances and minimizing the load loss of the highway energy system under fault conditions; the feasible construction locations of distributed power sources are obtained, and the distribution network nodes are obtained; the multi-objective optimization model is used to determine the specific construction locations of distributed power sources in the distribution network nodes, and the distributed power source capacity is determined to complete the site selection and capacity determination of distributed power sources.

[0023] The beneficial effect of the above further scheme is that it can obtain the site selection and capacity determination results of wind and solar power stations that maximize the benefits and resilience improvements.

[0024] Furthermore, the expression of the multi-objective optimization model is:

[0025] C=C1+C2+C3

[0026]

[0027]

[0028]

[0029] Among them, C is a multi-objective optimization model; C1 is the construction cost of the wind and solar power station; C2 is the economic benefit of the operation stage under normal conditions; C3 is the annual load loss cost under fault conditions; f is the depreciation coefficient; e is the cost of the wind and solar power station construction ... w is the unit installed capacity construction cost of wind power; u w N is the annual operation and maintenance cost per unit installed capacity of wind power; W Assemble for wind farms; is the rated power of wind farm i1; e s is the construction cost of photovoltaic unit installed capacity; u s N is the annual operation and maintenance cost per unit installed capacity of photovoltaic power generation; S Assemble for photovoltaic power station; is the rated power of the photovoltaic power station i2; W is the set of typical output scenarios under normal conditions; Γ w is the probability of the typical output scenario w under normal conditions; C w,t is the distribution network electricity price at time t under the typical output scenario w under normal conditions; T is the number of divided time stages; N is the set of distribution network nodes; is the active power output of the upper distribution network node i at time t under the typical output scenario w under normal conditions; ij is the distribution network line from node i to node j; l is the set of distribution network lines; is the square of the current amplitude on the distribution network line ij at time t under the typical output scenario w under normal conditions; is the resistance of the distribution network line ij at time t in the typical output scenario w under normal conditions; N d is the number of extreme disasters that occur in a year; S is the set of typical output scenarios under fault conditions; π s is the probability of occurrence of typical output scenario s under fault conditions; T z is the duration of the fault condition; c s is the load loss cost coefficient of ordinary load; w i is the importance weight of the critical load at the distribution network node i; is the electricity price at time t1 under the typical output scenario s under fault conditions; is the load reduction percentage at the distribution network node i at time t1 in the typical output scenario s under fault conditions; is the active load at the node i of the distribution network at time t1 in the typical output scenario s under fault conditions; Δt1 is the time period length.

[0030] The beneficial effect of the above further scheme is that by determining the objective function of the site selection and sizing scheme, the goals of minimizing economic costs and improving resilience are guaranteed.

[0031] Furthermore, the rated power of the photovoltaic power station is and wind farm rated power The constraints are:

[0032]

[0033]

[0034]

[0035]

[0036] Among them, m k is a non-negative integer variable; P PV is the rated power of the unit photovoltaic power station; P PV,max The maximum rated power installed for the photovoltaic power station; P WT is the rated power of the unit wind farm; P WT,max Install the maximum power rating for the wind farm.

[0037] The beneficial effect of the above further scheme is that by setting the unit wind and solar power station rated power, the rationality of the wind and solar power station rated power configuration is guaranteed.

[0038] Furthermore, the normal operating constraints of the highway energy system are:

[0039]

[0040]

[0041]

[0042]

[0043]

[0044]

[0045]

[0046]

[0047] in, is the PV active output of the upper distribution network node i at time t under the typical output scenario w under normal conditions; is the wind farm active output of the upper distribution network node i at time t under the typical output scenario w under normal conditions; is the active load of the upper distribution network node i at time t under the typical output scenario w under normal conditions; is the line transmission active power of the distribution network line ij at time t under the typical output scenario w under normal conditions; ki is the distribution network line from node k to node i; is the line transmission active power of distribution network line ki at time t under the typical output scenario w under normal conditions; is the resistance of the distribution network line ki at time t in the typical output scenario w under normal conditions; is the square of the current amplitude of the distribution network line ki at time t under the typical output scenario w under normal conditions; is the transmission reactive power of node i in the upper distribution network at time t under the typical output scenario w under normal conditions; is the photovoltaic reactive power output of distribution network node i at time t under the typical output scenario w under normal conditions; is the reactive power output of the wind farm at the distribution network node i at time t under the typical output scenario w under normal conditions; is the reactive load of the distribution network node i at time t under the typical output scenario w under normal conditions; is the transmission reactive power of distribution network line ij at time t under the typical output scenario w under normal conditions; is the transmission reactive power of distribution network line ki at time t under the typical output scenario w under normal conditions; is the reactance of distribution network line ki at time t under the typical output scenario w under normal conditions; is the square of the voltage amplitude of the distribution network node j at time t in the typical output scenario w under normal conditions; is the square of the voltage amplitude of the distribution network node i at time t in the typical output scenario w under normal conditions; is the reactance of the distribution network line ij at time t under the typical output scenario w under normal conditions; is the minimum value of the square of the node voltage amplitude under the typical output scenario w under normal conditions; is the maximum value of the square of the node voltage amplitude under the typical output scenario w under normal conditions; is the maximum square value of the line current under the typical output scenario w under normal conditions; P min,w is the minimum value of the active power transmitted by the line under the typical output scenario w under normal conditions; P max,w is the maximum value of the active power transmitted by the line under the typical output scenario w under normal conditions; Q min,w is the minimum value of the line transmission reactive power under the typical output scenario w under normal conditions; Q max,w It is the maximum value of the reactive power transmitted by the line under the typical output scenario w under normal conditions.

[0048] The beneficial effect of the above further scheme is: by setting the operation constraints of the highway energy system under normal circumstances, the stable operation of the system under normal circumstances is guaranteed.

[0049] Furthermore, the operation constraints under the condition of highway energy system failure are:

[0050]

[0051]

[0052]

[0053]

[0054]

[0055]

[0056]

[0057]

[0058]

[0059]

[0060]

[0061]

[0062]

[0063] in, is the transmitted active power of node i in the upper distribution network at time t1 under the typical output scenario s under fault conditions; is the PV active output of distribution network node i at time t1 under the typical output scenario s under fault conditions; is the active output of the wind farm at the distribution network node i at time t1 under the typical output scenario s under fault conditions; is the line transmission active power of distribution network line ij at time t1 under the typical output scenario s under fault conditions; ki is the distribution network line from node k to node i; is the line transmission active power of distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the resistance of the distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the square of the current amplitude of the distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the transmitted reactive power of distribution network node i at time t1 under the typical output scenario s under fault conditions; is the PV reactive output of distribution network node i at time t1 under the typical output scenario s under fault conditions; is the reactive power output of the wind farm at the node i of the distribution network at time t1 under the typical output scenario s under fault conditions; is the reactive load of the distribution network node i at time t1 under the typical output scenario s under fault conditions; is the transmission reactive power of distribution network line ij at time t1 under the typical output scenario s under fault conditions; is the transmission reactive power of distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the reactance of distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the square of the voltage amplitude of the distribution network node j at time t1 in the typical output scenario s under fault conditions; M is a constant; is a 0-1 variable of the distribution network line ij at time t1 under the typical output scenario s under fault conditions, indicating the on-off status of the line; is the square of the voltage amplitude of the distribution network node i at time t1 in the typical output scenario s under fault conditions; is the resistance of the distribution network line ij at time t1 under the typical output scenario s under fault conditions; is the reactance of distribution network line ij at time t1 under the typical output scenario s under fault conditions; is the square of the current amplitude of the distribution network line ij at time t1 under the typical output scenario s under fault conditions; is the minimum value of the square of the node voltage amplitude under the typical output scenario s under fault conditions; is the maximum value of the square of the node voltage amplitude under the typical output scenario s under fault conditions; is the maximum square value of the line current under the typical output scenario s under fault conditions; P min,s is the minimum value of the active power transmitted by the line under the typical output scenario s under fault conditions; P max,s is the maximum value of the active power transmitted by the line under the typical output scenario s under fault conditions; Q min,s is the minimum value of the reactive power transmitted by the line under the typical output scenario s under fault conditions; Q max,s is the maximum value of the reactive power transmitted by the line under the typical output scenario s under fault conditions; is a radial constraint auxiliary 0-1 variable. When the distribution network line ij is connected and j is the parent node of i, otherwise, is a radial constraint auxiliary 0-1 variable. When the distribution network line ji is connected and i is the parent node of j, otherwise, N i is the set of nodes directly connected to i; Auxiliary 0-1 variable for the radial constraint of the distribution network line from distribution network node 1 to distribution network node j.

[0064] The beneficial effect of the above further scheme is: by setting the operation constraints in the event of a highway energy system failure, the stable operation of the system in the event of a failure is ensured. BRIEF DESCRIPTION OF THE DRAWINGS

[0065] Figure 1 The figure is a flow chart of the method of the present invention.

[0066] Figure 2 This is a topological diagram of an improved IEEE 14-node system in an embodiment of the present invention.

[0067] Figure 3 1 is a topological diagram of a highway network in an embodiment of the present invention.

[0068] Figure 4 It is a schematic diagram of the planning results of the wind and solar power station in an embodiment of the present invention. DETAILED DESCRIPTION

[0069] The specific implementation modes of the present invention are described below so that those skilled in the art can understand the present invention. However, it should be clear that the present invention is not limited to the scope of the specific implementation modes. For those of ordinary skill in the art, as long as various changes are within the spirit and scope of the present invention as defined and determined by the attached claims, these changes are obvious, and all inventions and creations utilizing the concept of the present invention are protected.

[0070] like Figure 1 As shown, in one embodiment of the present invention, a distributed power source site selection and capacity determination method applied to a highway energy system includes the following steps:

[0071] Establish output constraints of distributed power sources, and process historical data of distributed power sources according to the output constraints of distributed power sources to obtain several output scenarios;

[0072] Clustering algorithms are used to cluster various output scenarios to obtain typical output scenarios;

[0073] According to typical output scenarios, with the goal of minimizing the cost of the highway energy system under normal circumstances and minimizing the load loss of the highway energy system under fault conditions, the construction location and capacity of the distributed power source are obtained, and the site selection and sizing of the distributed power source are completed.

[0074] In this embodiment, two improved IEEE14 node systems and a highway network in a certain place in the southwest are taken as examples. The topological structure diagram of the two improved IEEE14 node systems is as follows: Figure 2 As shown in the figure, the topology of the original system is represented by solid lines, and the dotted lines represent the unclosed tie switches. The reference voltage of the system is 35kV, the reference power is 100MW, and the voltage range is 0.85~1.05pu. The topology of the highway network in a certain place in the southwest is as follows Figure 3 As shown in Figure 1, the highway network is divided into different energy-consuming facilities, the energy characteristics and needs of different energy-consuming facilities are considered, the main energy consumption scenarios of highways are constructed, and the load parameters under different energy consumption scenarios are calculated. The correspondence between highway network nodes and distribution network nodes is shown in Table 1. Table 1 is a schematic table of the correspondence between highway network nodes and distribution network nodes.

[0075] Table 1

[0076] Highway network node Corresponding distribution network node 1 2 2 1 3 11 4 14 5 13 6 4 7 12 8 16 9 19 10 21 11 7 12 5 13 3 14 8 15 9 16 26 17 15 18 17 19 22 20 23 21 18 22 27 23 28 24 24 25 20

[0077] The expression of the output constraint of the distributed power source is:

[0078]

[0079]

[0080]

[0081]

[0082] T C (t) = T temp (t)+(N oct -20)I(t) / I stc

[0083]

[0084]

[0085] in, is the maximum output power of wind farm i1 at time t; v(t) is the wind speed at the hub of the wind turbine generator at time t; v in is the cut-in wind speed of the wind turbine; v out is the cut-out wind speed of the wind turbine; is the rated power of wind farm i1; v r is the rated wind speed of the wind turbine; is the actual output active power of wind farm i1 at time t; is the minimum output reactive power of wind farm i1; is the actual output reactive power of wind farm i1 at time t; is the maximum output reactive power of wind farm i1; is the maximum output power of photovoltaic power station i2 at time t; I(t) is the light intensity at time t; φ is the derating coefficient of photovoltaic cell; is the rated power of the photovoltaic power station i2; β is the power temperature coefficient; T C (t) is the actual operating temperature of the photovoltaic cell at time t; T stc (t) is the standard operating temperature of the photovoltaic cell at time t; stc is the light intensity under standard test conditions; T temp (t) is the ambient temperature at time t; N oct is the nominal temperature of the photovoltaic cell; is the actual output active power of photovoltaic power station i2 at time t; is the minimum output reactive power of photovoltaic power station i2; is the actual output reactive power of photovoltaic power station i2 at time t; is the maximum output reactive power of photovoltaic power station i2.

[0086] The typical output scenarios include typical output scenarios under normal conditions and typical output scenarios under fault conditions.

[0087] In this embodiment, the output constraints of wind and solar power are taken into consideration, and different output scenarios are generated based on historical data of wind and solar power. The specific steps are as follows:

[0088] ① Collect historical data of wind and solar power, including wind speed, temperature and other data;

[0089] ② Based on the output constraints of wind and solar power, the historical data of wind and solar power, wind speed, temperature, etc. are substituted into the output constraint model, and the wind and solar power output data under different scenarios are generated through calculation.

[0090] In this embodiment, the improved k-means method is used to cluster and reduce a large number of generated scenarios to obtain typical output scenarios. The specific steps are as follows:

[0091] ① Consider the natural distribution of a large amount of generated wind and solar power output scene data and determine the initial cluster center;

[0092] ② Using the k-means method, a large amount of wind and solar scene data with initial cluster centers are clustered and reduced to construct typical wind and solar output scenarios under normal and fault conditions.

[0093] According to the typical output scenario, the location and capacity determination strategy of distributed power sources is obtained with the goal of minimizing the cost of the highway energy system under normal conditions and minimizing the load loss of the highway energy system under fault conditions, which is specifically:

[0094] A multi-objective optimization model is constructed with the goal of minimizing the cost of the highway energy system under normal circumstances and minimizing the load loss of the highway energy system under fault conditions; the feasible construction locations of distributed power sources are obtained, and the distribution network nodes are obtained; the multi-objective optimization model is used to determine the specific construction locations of distributed power sources in the distribution network nodes, and the distributed power source capacity is determined to complete the site selection and capacity determination of distributed power sources.

[0095] The expression of the multi-objective optimization model is:

[0096] C=C1+C2+C3

[0097]

[0098]

[0099]

[0100] Among them, C is a multi-objective optimization model; C1 is the construction cost of the wind and solar power station; C2 is the economic benefit of the operation stage under normal conditions; C3 is the annual load loss cost under fault conditions; f is the depreciation coefficient; e is the cost of the wind and solar power station construction ... w is the unit installed capacity construction cost of wind power; u w N is the annual operation and maintenance cost per unit installed capacity of wind power;W Assemble for wind farms; is the rated power of wind farm i1; e s is the construction cost of photovoltaic unit installed capacity; u s N is the annual operation and maintenance cost per unit installed capacity of photovoltaic power generation; S Assemble for photovoltaic power station; is the rated power of the photovoltaic power station i2; W is the set of typical output scenarios under normal conditions; Γ w is the probability of the typical output scenario w under normal conditions; C w,t is the distribution network electricity price at time t under the typical output scenario w under normal conditions; T is the number of divided time stages; N is the set of distribution network nodes; is the active power output of the upper distribution network node i at time t under the typical output scenario w under normal conditions; ij is the distribution network line from node i to node j; l is the set of distribution network lines; is the square of the current amplitude on the distribution network line ij at time t under the typical output scenario w under normal conditions; is the resistance of the distribution network line ij at time t in the typical output scenario w under normal conditions; N d is the number of extreme disasters that occur in a year; S is the set of typical output scenarios under fault conditions; π s is the probability of occurrence of typical output scenario s under fault conditions; T z is the duration of the fault condition; c s is the load loss cost coefficient of ordinary load; w i is the importance weight of the critical load at the distribution network node i; is the electricity price at time t1 under the typical output scenario s under fault conditions; is the load reduction percentage at the distribution network node i at time t1 in the typical output scenario s under fault conditions; is the active load at the node i of the distribution network at time t1 in the typical output scenario s under fault conditions; Δt1 is the time period length.

[0101] The rated power of the photovoltaic power station and wind farm rated power The constraints are:

[0102]

[0103]

[0104]

[0105]

[0106] Among them, mk is a non-negative integer variable; P PV is the rated power of the unit photovoltaic power station; P PV,max The maximum rated power installed for the photovoltaic power station; P WT is the rated power of the unit wind farm; P WT,max Install the maximum power rating for the wind farm.

[0107] The operating constraints of the highway energy system under normal circumstances are:

[0108]

[0109]

[0110]

[0111]

[0112]

[0113]

[0114]

[0115]

[0116] in, is the PV active output of the upper distribution network node i at time t under the typical output scenario w under normal conditions; is the wind farm active output of the upper distribution network node i at time t under the typical output scenario w under normal conditions; is the active load of the upper distribution network node i at time t under the typical output scenario w under normal conditions; is the line transmission active power of the distribution network line ij at time t under the typical output scenario w under normal conditions; ki is the distribution network line from node k to node i; is the line transmission active power of distribution network line ki at time t under the typical output scenario w under normal conditions; is the resistance of the distribution network line ki at time t in the typical output scenario w under normal conditions; is the square of the current amplitude of the distribution network line ki at time t under the typical output scenario w under normal conditions; is the transmission reactive power of node i in the upper distribution network at time t under the typical output scenario w under normal conditions; is the photovoltaic reactive power output of distribution network node i at time t under the typical output scenario w under normal conditions; is the reactive power output of the wind farm at the distribution network node i at time t under the typical output scenario w under normal conditions; is the reactive load of the distribution network node i at time t under the typical output scenario w under normal conditions; is the transmission reactive power of distribution network line ij at time t under the typical output scenario w under normal conditions; is the transmission reactive power of distribution network line ki at time t under the typical output scenario w under normal conditions; is the reactance of distribution network line ki at time t under the typical output scenario w under normal conditions; is the square of the voltage amplitude of the distribution network node j at time t in the typical output scenario w under normal conditions; is the square of the voltage amplitude of the distribution network node i at time t in the typical output scenario w under normal conditions; is the reactance of the distribution network line ij at time t under the typical output scenario w under normal conditions; is the minimum value of the square of the node voltage amplitude under the typical output scenario w under normal conditions; is the maximum value of the square of the node voltage amplitude under the typical output scenario w under normal conditions; is the maximum square value of the line current under the typical output scenario w under normal conditions; P min,w is the minimum value of the active power transmitted by the line under the typical output scenario w under normal conditions; P max,w is the maximum value of the active power transmitted by the line under the typical output scenario w under normal conditions; Q min,w is the minimum value of the line transmission reactive power under the typical output scenario w under normal conditions; Q max,w It is the maximum value of the reactive power transmitted by the line under the typical output scenario w under normal conditions.

[0117] The operating constraints under highway energy system failure conditions are:

[0118]

[0119]

[0120]

[0121]

[0122]

[0123]

[0124]

[0125]

[0126]

[0127]

[0128]

[0129]

[0130]

[0131] in, is the transmitted active power of node i in the upper distribution network at time t1 under the typical output scenario s under fault conditions; is the PV active output of distribution network node i at time t1 under the typical output scenario s under fault conditions; is the active output of the wind farm at the distribution network node i at time t1 under the typical output scenario s under fault conditions; is the line transmission active power of distribution network line ij at time t1 under the typical output scenario s under fault conditions; ki is the distribution network line from node k to node i; is the line transmission active power of distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the resistance of the distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the square of the current amplitude of the distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the transmitted reactive power of distribution network node i at time t1 under the typical output scenario s under fault conditions; is the PV reactive output of distribution network node i at time t1 under the typical output scenario s under fault conditions; is the reactive power output of the wind farm at the node i of the distribution network at time t1 under the typical output scenario s under fault conditions; is the reactive load of the distribution network node i at time t1 under the typical output scenario s under fault conditions; is the transmission reactive power of distribution network line ij at time t1 under the typical output scenario s under fault conditions; is the transmission reactive power of distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the reactance of distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the square of the voltage amplitude of the distribution network node j at time t1 in the typical output scenario s under fault conditions; M is a constant; is a 0-1 variable of the distribution network line ij at time t1 under the typical output scenario s under fault conditions, indicating the on-off status of the line; is the square of the voltage amplitude of the distribution network node i at time t1 in the typical output scenario s under fault conditions; is the resistance of the distribution network line ij at time t1 under the typical output scenario s under fault conditions; is the reactance of distribution network line ij at time t1 under the typical output scenario s under fault conditions; is the square of the current amplitude of the distribution network line ij at time t1 under the typical output scenario s under fault conditions; is the minimum value of the square of the node voltage amplitude under the typical output scenario s under fault conditions; is the maximum value of the square of the node voltage amplitude under the typical output scenario s under fault conditions; is the maximum square value of the line current under the typical output scenario s under fault conditions; P min,s is the minimum value of the active power transmitted by the line under the typical output scenario s under fault conditions; P max,s is the maximum value of the active power transmitted by the line under the typical output scenario s under fault conditions; Q min,s is the minimum value of the reactive power transmitted by the line under the typical output scenario s under fault conditions; Q max,s is the maximum value of the reactive power transmitted by the line under the typical output scenario s under fault conditions; is a radial constraint auxiliary 0-1 variable. When the distribution network line ij is connected and j is the parent node of i, otherwise, is a radial constraint auxiliary 0-1 variable. When the distribution network line ji is connected and i is the parent node of j, otherwise, N i is the set of nodes directly connected to i; Auxiliary 0-1 variable for the radial constraint of the distribution network line from distribution network node 1 to distribution network node j.

[0132] In this embodiment, it is assumed that in a normal scenario, one typical day each in spring, summer, autumn and winter is taken as an example, and the switching of the contact switch is not considered; it is assumed that in a fault scenario, Line 1 and Line 17 are broken, the fault occurs on a typical day in spring, the fault period is 10:00-13:00, and lasts for three hours.

[0133] Table 2 is a schematic diagram of relevant parameters for wind and solar power station planning.

[0134] Table 2

[0135] Object parameter Numeric Construction cost (yuan / kW) 7780 Wind Farm Annual operation and maintenance cost (yuan / kW) 78 Depreciation factor 0.01 Construction cost (yuan / kW) 6900 Photovoltaic power station Annual operation and maintenance cost (yuan / kW) 69 Depreciation factor 0.01 other Planning period 1

[0136] The schematic diagram of wind and solar power station planning results is as follows Figure 4As shown, in order to verify that the present invention can effectively improve the economy and resilience of the highway energy system, the economic benefits before and after planning are compared. Table 3 is a comparison table of the results before and after planning.

[0137] Table 3

[0138]

[0139]

[0140] It can be seen from Table 3 that by planning the wind and solar power stations, the construction cost and economic cost of the normal operation stage of the wind and solar power stations after planning were reduced by 88.80028 million yuan, the annual load loss cost was reduced by 11,100 yuan, and the total cost was reduced by 88.80139 million yuan. This proves that the present invention can effectively improve the economy and resilience of the highway energy system.

Claims

1. A distributed power source site selection and capacity determination method applied to a highway energy system, characterized in that: The following steps are involved: Establish output constraints of distributed power sources, and process historical data of distributed power sources according to the output constraints of distributed power sources to obtain several output scenarios; Clustering algorithms are used to cluster various output scenarios to obtain typical output scenarios; According to typical output scenarios, with the goal of minimizing the cost of the highway energy system under normal circumstances and minimizing the load loss of the highway energy system under fault conditions, the construction location and capacity of the distributed power source are obtained, and the site selection and sizing of the distributed power source are completed.

2. The distributed power source site selection and capacity determination method applied to highway energy system according to claim 1 is characterized in that: The expression of the output constraint of the distributed power source is: T C (t)=T temp (t)+(N oct -20)I(t) / I stc in, is the maximum output power of wind farm i1 at time t; v(t) is the wind speed at the hub of the wind turbine generator at time t; v in is the cut-in wind speed of the wind turbine; v out is the cut-out wind speed of the wind turbine; is the rated power of wind farm i1; v r is the rated wind speed of the wind turbine; is the actual output active power of wind farm i1 at time t; is the minimum output reactive power of wind farm i1; is the actual output reactive power of wind farm i1 at time t; is the maximum output reactive power of wind farm i1; is the maximum output power of photovoltaic power station i2 at time t; I(t) is the light intensity at time t; φ is the derating coefficient of photovoltaic cell; is the rated power of the photovoltaic power station i2; β is the power temperature coefficient; T C (t) is the actual operating temperature of the photovoltaic cell at time t; T stc (t) is the standard operating temperature of the photovoltaic cell at time t; stc is the light intensity under standard test conditions; T temp (t) is the ambient temperature at time t; N oct is the nominal temperature of the photovoltaic cell; is the actual output active power of photovoltaic power station i2 at time t; is the minimum output reactive power of photovoltaic power station i2; is the actual output reactive power of photovoltaic power station i2 at time t; is the maximum output reactive power of photovoltaic power station i2.

3. The distributed power source site selection and capacity determination method applied to highway energy system according to claim 1 is characterized in that: The typical output scenarios include typical output scenarios under normal conditions and typical output scenarios under fault conditions.

4. The method for site selection and capacity determination of distributed power sources applied to highway energy systems according to claim 1 is characterized in that: According to the typical output scenario, the location and capacity determination strategy of distributed power sources is obtained with the goal of minimizing the cost of the highway energy system under normal conditions and minimizing the load loss of the highway energy system under fault conditions, which is specifically: A multi-objective optimization model is constructed with the goal of minimizing the cost of the highway energy system under normal circumstances and minimizing the load loss of the highway energy system under fault conditions; the feasible construction locations of distributed power sources are obtained, and the distribution network nodes are obtained; the multi-objective optimization model is used to determine the specific construction locations of distributed power sources in the distribution network nodes, and the distributed power source capacity is determined to complete the site selection and capacity determination of distributed power sources.

5. The method for selecting and sizing the distributed power source for a highway energy system according to claim 4 is characterized in that: The expression of the multi-objective optimization model is: C=C1+C2+C3 Among them, C is a multi-objective optimization model; C1 is the construction cost of the wind and solar power station; C2 is the economic benefit of the operation stage under normal conditions; C3 is the annual load loss cost under fault conditions; f is the depreciation coefficient; e is the cost of the wind and solar power station construction ... w is the unit installed capacity construction cost of wind power; u w N is the annual operation and maintenance cost per unit installed capacity of wind power; W Assemble for wind farms; is the rated power of wind farm i1; e s is the construction cost of photovoltaic unit installed capacity; u s N is the annual operation and maintenance cost per unit installed capacity of photovoltaic power generation; S Assemble for photovoltaic power station; is the rated power of the photovoltaic power station i2; W is the set of typical output scenarios under normal conditions; Γ w is the probability of the typical output scenario w under normal conditions; C w,t is the distribution network electricity price at time t under the typical output scenario w under normal conditions; T is the number of divided time stages; N is the set of distribution network nodes; is the active power output of the upper distribution network node i at time t under the typical output scenario w under normal conditions; ij is the distribution network line from node i to node j; l is the set of distribution network lines; is the square of the current amplitude on the distribution network line ij at time t under the typical output scenario w under normal conditions; is the resistance of the distribution network line ij at time t in the typical output scenario w under normal conditions; N d is the number of extreme disasters that occur in a year; S is the set of typical output scenarios under fault conditions; π s is the probability of occurrence of typical output scenario s under fault conditions; T z is the duration of the fault condition; c s is the load loss cost coefficient of ordinary load; w i is the importance weight of the critical load at the distribution network node i; is the electricity price at time t1 under the typical output scenario s under fault conditions; is the load reduction percentage at the distribution network node i at time t1 in the typical output scenario s under fault conditions; is the active load at the node i of the distribution network at time t1 in the typical output scenario s under fault conditions; Δt1 is the time period length.

6. The method for selecting and sizing the distributed power source for use in a highway energy system according to claim 5 is characterized in that: The rated power of the photovoltaic power station and wind farm rated power The constraints are: Among them, m k is a non-negative integer variable; P PV is the rated power of the unit photovoltaic power station; P PV,max The maximum rated power installed for the photovoltaic power station; P WT is the rated power of the unit wind farm; P WT,max Install the maximum power rating for the wind farm.

7. The method for site selection and capacity determination of distributed power sources applied to highway energy systems according to claim 5 is characterized in that: The operating constraints of the highway energy system under normal circumstances are: in, is the PV active output of the upper distribution network node i at time t under the typical output scenario w under normal conditions; is the wind farm active output of the upper distribution network node i at time t under the typical output scenario w under normal conditions; is the active load of the upper distribution network node i at time t under the typical output scenario w under normal conditions; is the line transmission active power of the distribution network line ij at time t under the typical output scenario w under normal conditions; ki is the distribution network line from node k to node i; is the line transmission active power of distribution network line ki at time t under the typical output scenario w under normal conditions; is the resistance of the distribution network line ki at time t in the typical output scenario w under normal conditions; is the square of the current amplitude of the distribution network line ki at time t under the typical output scenario w under normal conditions; is the transmission reactive power of node i in the upper distribution network at time t under the typical output scenario w under normal conditions; is the photovoltaic reactive power output of distribution network node i at time t under the typical output scenario w under normal conditions; is the reactive power output of the wind farm at the distribution network node i at time t under the typical output scenario w under normal conditions; is the reactive load of the distribution network node i at time t under the typical output scenario w under normal conditions; is the transmission reactive power of distribution network line ij at time t under the typical output scenario w under normal conditions; is the transmission reactive power of distribution network line ki at time t under the typical output scenario w under normal conditions; is the reactance of distribution network line ki at time t under the typical output scenario w under normal conditions; is the square of the voltage amplitude of the distribution network node j at time t in the typical output scenario w under normal conditions; is the square of the voltage amplitude of the distribution network node i at time t in the typical output scenario w under normal conditions; is the reactance of the distribution network line ij at time t under the typical output scenario w under normal conditions; is the minimum value of the square of the node voltage amplitude under the typical output scenario w under normal conditions; is the maximum value of the square of the node voltage amplitude under the typical output scenario w under normal conditions; is the maximum square value of the line current under the typical output scenario w under normal conditions; P min,w is the minimum value of the active power transmitted by the line under the typical output scenario w under normal conditions; P max,w is the maximum value of the active power transmitted by the line under the typical output scenario w under normal conditions; Q min,w is the minimum value of the line transmission reactive power under the typical output scenario w under normal conditions; Q max,w It is the maximum value of the reactive power transmitted by the line under the typical output scenario w under normal conditions.

8. The method for selecting and sizing the distributed power source for use in a highway energy system according to claim 5 is characterized in that: The operating constraints under highway energy system failure conditions are: in, is the transmitted active power of node i in the upper distribution network at time t1 under the typical output scenario s under fault conditions; is the PV active output of distribution network node i at time t1 under the typical output scenario s under fault conditions; is the active output of the wind farm at the distribution network node i at time t1 under the typical output scenario s under fault conditions; is the line transmission active power of distribution network line ij at time t1 under the typical output scenario s under fault conditions; ki is the distribution network line from node k to node i; is the line transmission active power of distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the resistance of the distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the square of the current amplitude of the distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the transmitted reactive power of distribution network node i at time t1 under the typical output scenario s under fault conditions; is the PV reactive output of distribution network node i at time t1 under the typical output scenario s under fault conditions; is the reactive power output of the wind farm at the node i of the distribution network at time t1 under the typical output scenario s under fault conditions; is the reactive load of the distribution network node i at time t1 under the typical output scenario s under fault conditions; is the transmission reactive power of distribution network line ij at time t1 under the typical output scenario s under fault conditions; is the transmission reactive power of distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the reactance of distribution network line ki at time t1 under the typical output scenario s under fault conditions; is the square of the voltage amplitude of the distribution network node j at time t1 in the typical output scenario s under fault conditions; M is a constant; is a 0-1 variable of the distribution network line ij at time t1 under the typical output scenario s under fault conditions, indicating the on-off status of the line; is the square of the voltage amplitude of the distribution network node i at time t1 in the typical output scenario s under fault conditions; is the resistance of the distribution network line ij at time t1 under the typical output scenario s under fault conditions; is the reactance of distribution network line ij at time t1 under the typical output scenario s under fault conditions; is the square of the current amplitude of the distribution network line ij at time t1 under the typical output scenario s under fault conditions; is the minimum value of the square of the node voltage amplitude under the typical output scenario s under fault conditions; is the maximum value of the square of the node voltage amplitude under the typical output scenario s under fault conditions; is the maximum square value of the line current under the typical output scenario s under fault conditions; P min,s is the minimum value of the active power transmitted by the line under the typical output scenario s under fault conditions; P max,s is the maximum value of the active power transmitted by the line under the typical output scenario s under fault conditions; Q min,s is the minimum value of the reactive power transmitted by the line under the typical output scenario s under fault conditions; Q max,s is the maximum value of the reactive power transmitted by the line under the typical output scenario s under fault conditions; is a radial constraint auxiliary 0-1 variable. When the distribution network line ij is connected and j is the parent node of i, otherwise, is a radial constraint auxiliary 0-1 variable. When the distribution network line ji is connected and i is the parent node of j, otherwise, N i is the set of nodes directly connected to i; Auxiliary 0-1 variable for the radial constraint of the distribution network line from distribution network node 1 to distribution network node j.

Images