Distributed photovoltaic energy storage point selection layout optimization method

Abstract

The invention discloses a distributed photovoltaic energy storage point selection layout optimization method, relates to the technical field of new energy power systems, and solves the technical problem that point selection layout is not optimized enough due to the fact that factors considered during point selection are not comprehensive enough and systematic analysis is lacked. The limitation that a traditional method only pays attention to single or few factors is overcome, through comprehensive consideration, the area suitable for construction of the distributed photovoltaic energy storage system can be more accurately selected, the adaptability of the system to the surrounding environment and the power grid is improved, and the construction efficiency of the distributed photovoltaic energy storage system is improved. A GIS technology is fully utilized to carry out systematic analysis from macroscopically dividing a to-be-analyzed area to microcosmically dividing a standard layout area on a research area, the key factor of substation power grid node voltage deviation is considered in the point selection process, and optimization analysis is carried out on a point selection area according to the parameters to generate optimization information. The method is helpful for selecting an area with small influence on the stability of the power grid, and guarantees the safe and stable operation of the power grid.

Description

Technical Field

[0001] The present invention relates to the technical field of new energy power systems, and in particular to a method for optimizing the layout of distributed photovoltaic energy storage site selection. Background Art

[0002] As the demand for clean energy continues to increase, distributed photovoltaic energy storage systems have received widespread attention as a sustainable energy solution. In the construction of distributed photovoltaic energy storage systems, site selection and layout is a key issue, which directly affects the system's performance, cost, and impact on the surrounding environment and power grid.

[0003] According to the publication number CN113904374A, a distributed photovoltaic energy storage site selection and layout optimization method is disclosed, which includes: S1, considering the number of layout points and the economic efficiency of energy storage configuration, establishing an objective function, and the objective function includes a node voltage fluctuation index and a comprehensive benefit index; S2, considering the capacity of each point, constructing constraints that at least include energy storage layout point constraints and energy storage power and capacity constraints; S3, initializing the system model, calculating the power flow distribution of the system model, and substituting the power flow calculation results into the objective function and the constraints to perform optimization and obtain the optimal solution.

[0004] Traditional site selection methods may have problems such as incomplete considerations and lack of systematic analysis. They often only focus on individual factors such as solar energy resources or grid access convenience, while ignoring the comprehensive consideration of multiple factors such as topography, residential density, load demand matching, and the impact on grid stability.

[0005] At the same time, in the process of analyzing geographic data and regional characteristics, there is a lack of effective quantitative calculation methods and comprehensive evaluation means, resulting in less than optimized site selection and layout. Summary of the invention

[0006] In view of the shortcomings of the prior art, the present invention provides a method for optimizing the site selection and layout of distributed photovoltaic energy storage, which solves the problem that the factors considered when selecting sites are not comprehensive enough and there is a lack of systematic analysis, resulting in a less optimized site selection layout.

[0007] To achieve the above objectives, the present invention is implemented through the following technical solutions: a distributed photovoltaic energy storage site selection layout optimization method, the method specifically includes the following steps:

[0008] Step S001, obtaining geographic data of the study area based on GIS technology, and dividing the study area into multiple groups of areas to be analyzed;

[0009] Step S002, comprehensively analyzing the geographic data of the area to be analyzed, and selecting a standard layout area by calculating the comprehensive area value of the area to be analyzed;

[0010] Step S003, segmenting the obtained standard layout area to obtain multiple groups of segmented areas, and obtaining area information corresponding to the segmented areas, analyzing the area information and calculating the selected point values ​​of the segmented areas;

[0011] Step S004, screening the segmented areas according to the calculated selected point values ​​of the segmented areas to obtain selected point areas, and performing optimization analysis based on node voltage deviations of the substations in the selected point areas to generate optimization information;

[0012] Step S005: display the generated optimization information to the corresponding operator.

[0013] As a further solution of the present invention, the method of selecting the standard layout area in step S002 is:

[0014] All the areas to be analyzed are obtained and labeled as i, and i=1, 2, ..., j, where j represents the number of the areas to be analyzed. Then, the terrain value Di, solar energy resource value Yi and site selection distance value Li corresponding to the area to be analyzed i are calculated respectively, and the obtained terrain value Di, solar energy resource value Yi and site selection distance value Li are substituted into the formula Qi=(Di+Yi+Li)×θ to calculate the comprehensive area value Qi of the area to be analyzed i, where θ is the preset proportional coefficient, and the specific value is set by the operator. Then, the area to be analyzed corresponding to the maximum comprehensive area value Qi is selected as the standard layout area.

[0015] As a further solution of the present invention, the terrain value in step S002 is calculated as follows:

[0016] Based on GIS technology, the topography of the area i to be analyzed is obtained. At the same time, the proportion of the high mountain area in the area i to be analyzed is obtained according to the topography information, recorded as Ri, and the slope value of the area i to be analyzed is calculated and recorded as Pi. Then, the terrain undulation of the area i to be analyzed is obtained and recorded as Si. At the same time, the obtained parameters are substituted into the calculation formula The terrain value D i of the area to be analyzed is calculated, where a1, a2 and a3 are weight coefficients, and a1+a2+a3=1.

[0017] As a further solution of the present invention, the calculation method of the solar energy resource value Yi in step S002 is:

[0018] Get the average sunshine duration per unit time of the area to be analyzed, record it as Ti, and get the average sunshine weather Gi of the area to be analyzed based on historical data, then get the annual average solar radiation intensity corresponding to the area to be analyzed, record it as F i, and substitute the obtained parameters into the formula The solar energy resource value Yi corresponding to the area to be analyzed is calculated, where b1, b2 and b3 are weight coefficients, and b1+b2+b3=1.

[0019] As a further solution of the present invention, the calculation method of the site selection distance value in step S002 is:

[0020] Obtain all substations in the study area, and calculate the distance between the substation and the area to be analyzed i respectively, and calculate the mean of all distances as the site selection distance value Li of the area to be analyzed.

[0021] As a further solution of the present invention, the specific method of calculating the selected point value of the segmentation area in step S003 is:

[0022] Obtain a standard layout area, and divide the standard layout area into multiple groups of divided areas and record them as o, where o=1, 2, ..., p, where p represents the number of divided areas, then obtain the residential density of the divided area o and record it as Jo, and obtain the contact line fluctuation rate Xo of the divided area o;

[0023] Calculate the load demand matching value Mo of the segmented area o, obtain the load data of the segmented area o, and establish the corresponding load curve, and extract the peak load P from the load curve peak , valley load P valley , average load P avg and load factor LR, where the load factor is calculated according to the formula Calculate and substitute the obtained parameters into the formula The load demand matching degree Mo is calculated, where α is the preset proportional coefficient;

[0024] Substitute the obtained parameters into the formula Calculate the selected point value ω corresponding to the segmented area o o , where μ is the preset proportional coefficient.

[0025] As a further solution of the present invention, the specific method of obtaining the selected point area in step S004 is:

[0026] Sort the selected point values ​​of the segmentation areas from large to small, and select the segmentation area corresponding to the largest selected point value as the selected point area.

[0027] As a further solution of the present invention, the specific method of generating the optimization information in step S004 is:

[0028] Get the substation in the selected area, then get all the grid nodes of the substation, denoted as r, and r = 1, 2, ..., h, where h represents the number of grid nodes, and calculate the active network loss corresponding to the substation. The specific calculation formula is: The active network loss P of the power grid node is calculated according to the formula loss , where P lossrepresents the active network loss of the system, Vr represents the voltage amplitude of node r, G rh , B rh , δ rh Respectively represent the conductance, susceptance and phase angle difference between nodes r and h, h is the total number of nodes; ∑ represents the set of nodes connected to node r;

[0029] Then substitute the obtained parameters into the formula Calculate the node voltage deviation Ur of the substation grid node r grid , where η is the preset proportional coefficient, the specific value is set by the operator, Ur is the voltage amplitude of the r node; U N is the standard voltage amplitude, and the node voltage deviation Ur is calculated in order from small to large. grid Sorting generates optimization information.

[0030] The present invention provides a method for optimizing the layout of distributed photovoltaic energy storage sites. Compared with the prior art, it has the following beneficial effects:

[0031] The present invention overcomes the limitation of traditional methods that only focus on single or a few factors by considering multiple factors such as topography, solar energy resources, site selection distance, residential density, tie line fluctuation rate, load demand matching degree, etc. Through comprehensive consideration, it can more accurately select areas suitable for the construction of distributed photovoltaic energy storage systems, improve the adaptability of the system to the surrounding environment and the power grid, make full use of GIS technology to conduct a systematic analysis of the research area from macro-division of the area to be analyzed to micro-division of the standard layout area, consider the key factor of the voltage deviation of the substation power grid node in the site selection process, calculate the active network loss and node voltage deviation, and optimize the site selection area according to these parameters to generate optimization information. This helps to select areas with less impact on the stability of the power grid and ensure the safe and stable operation of the power grid. BRIEF DESCRIPTION OF THE DRAWINGS

[0032] Figure 1 It is a step method diagram of the present invention;

[0033] Figure 2 It is a flow chart of the present invention. DETAILED DESCRIPTION

[0034] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0035] See also Figure 1 and Figure 2The present application provides a method for optimizing the layout of distributed photovoltaic energy storage points, which specifically includes the following steps:

[0036] Step S001, based on GIS technology, the geographical data of the study area is obtained, and the study area is divided into multiple groups of areas to be analyzed. The study area here is the site selection point of photovoltaic energy storage, and the geographical data includes topography, land use type, solar energy resource distribution and other information.

[0037] Step S002, comprehensively analyzing the geographic data of the area to be analyzed, and selecting a standard layout area by calculating the comprehensive area value of the area to be analyzed.

[0038] Get all the areas to be analyzed and label them as i, where i=1, 2, ..., j, where j represents the number of areas to be analyzed, and then calculate the terrain value D i, solar energy resource value Y i and site selection distance value L i corresponding to the area to be analyzed i respectively;

[0039] The specific method for calculating the terrain value Di of the area to be analyzed is: based on GIS technology, the topography of the area to be analyzed is obtained, and the proportion of the high mountain area of ​​the area to be analyzed is obtained according to the topography information, which is recorded as Ri. The proportion of the high mountain area refers to the ratio of the area of ​​the high mountain area to the area to be analyzed. The slope value of the area to be analyzed is calculated and recorded as Pi. The specific calculation process is: using trigonometric functions in and are the rates of change of the digital elevation model in the x and y directions respectively. Then, the terrain relief of the area i to be analyzed is obtained and recorded as S i. The specific calculation process is: determine a 3×3 analysis window. If the resolution of the digital elevation model is 30 meters, then a cell represents a ground range of 30 meters × 30 meters. In each selected window, find the maximum value Zmax and the minimum value Zmin in the digital elevation model data, and calculate the terrain relief corresponding to a single window according to the formula RD=Zmax-Zmin. Similarly, calculate the terrain relief corresponding to all windows, and calculate the average of the terrain relief of all windows as the terrain relief of the area i to be analyzed. At the same time, substitute the obtained parameters into the calculation formula Calculate the terrain value Di of the area to be analyzed i, where a1, a2 and a3 are weight coefficients, and a1+a2+a3=1;

[0040] The specific method for calculating the solar energy resource value Yi of the area to be analyzed is as follows: obtain the average sunshine duration per unit time of the area to be analyzed, record it as Ti, and obtain the average sunshine weather Gi of the area to be analyzed based on historical data. The sunshine weather here means that the weather conditions are good. The specific good weather conditions are obtained through the corresponding weather software. At the same time, the average sunshine weather is the average value within a year. Then obtain the annual average solar radiation intensity corresponding to the area to be analyzed, record it as Fi, and substitute the obtained parameters into the formula Calculate the solar energy resource value Yi corresponding to the area i to be analyzed, where b1, b2 and b3 are weight coefficients, and b1+b2+b3=1;

[0041] The specific method of calculating the site selection distance value Li of the area to be analyzed is: obtain all substations in the study area, and calculate the distance between the substation and the area to be analyzed i respectively, and calculate the average of all distances as the site selection distance value Li of the area to be analyzed;

[0042] In a specific example, for example, assuming that there are 5 substations in the study area, their coordinates are A (10, 20), B (30, 40), C (50, 60), D (70, 80), E (90, 100), and the coordinates of the area to be analyzed is (45, 55). The distance between the analyzed area i and the substation A is calculated. Similarly, the above calculation process is used to calculate the distances of B, C, D and E in turn, and finally the mean distance is calculated.

[0043] Substitute the acquired terrain value Di, solar energy resource value Yi and site selection distance value Li into the formula Qi=(Di+Yi+Li)×θ to calculate the comprehensive regional value Qi of the area to be analyzed i, where θ is the preset proportional coefficient, and the specific value is set by the operator. Then select the area to be analyzed corresponding to the maximum comprehensive regional value Qi as the standard layout area.

[0044] Step S003, segmenting the obtained standard layout area to obtain a plurality of groups of segmented areas, and acquiring area information corresponding to the segmented areas, analyzing the area information and calculating the selected point values ​​of the segmented areas.

[0045] The standard layout area is obtained, and the standard layout area is divided into multiple groups of divided areas and recorded as o, where o=1, 2, ..., p, where p represents the number of divided areas. Then the residential density of the divided area o is obtained and recorded as Jo, and the residential density is expressed as the ratio of the number of people in the divided area to the area of ​​the divided area. For example, it is known through statistics that the number of people in the divided area o is 5000, and the area of ​​the area is measured to be 2 square kilometers, then the residential density Jo=5000÷2=2500 people / square kilometer. At the same time, the contact line fluctuation rate Xo of the divided area o is obtained. The specific contact line fluctuation rate is reflected by the contact line power change rate and the contact line power fluctuation rate, where the contact line power change rate is defined as: Where D grid is the power change rate of the tie line, P grid max and P grid min are the maximum and minimum values ​​of the tie line power in the Δt time interval, and the tie line power standard deviation is used to represent the fluctuation rate of the tie line power. The specific expression is: Where X grid is the power change rate of the tie line, P grid av is the average power of the tie line, P grid c is the power of the tie line at the cth time point, and the load demand matching value Mo of the segmented area o is calculated as follows: obtain the load data of the segmented area o, including the power load demand at different time intervals (such as every 15 minutes, every hour, etc.), and establish the corresponding load curve, and extract the peak load P from the load curve peak , valley load P valley , average load P avg and load factor LR, where the load factor is calculated according to the formula Calculate and substitute the obtained parameters into the formula The load demand matching degree Mo is calculated, where α is the preset proportional coefficient;

[0046] Substitute the obtained parameters into the formula Calculate the selected point value ω corresponding to the segmented area o o , where μ is the preset proportional coefficient.

[0047] Step S004, screening the segmented areas according to the calculated selected point values ​​of the segmented areas to obtain selected point areas, and performing optimization analysis based on node voltage deviations of the substations in the selected point areas to generate optimization information.

[0048] Sort the selected point values ​​of the segmented areas from large to small, and select the segmented area corresponding to the largest selected point value as the selected point area, and obtain the substation in the selected point area. Then obtain all the grid nodes of the substation and record them as r, and r = 1, 2, ..., h, where h represents the number of grid nodes. At the same time, calculate the active network loss corresponding to the substation. The specific calculation formula is: According to the above formula, the active network loss P of the power grid node is calculated loss , where P loss represents the active network loss of the system, Vr represents the voltage amplitude of node r, G rh , B rh , δ rh Respectively represent the conductance, susceptance and phase angle difference between nodes r and h, h is the total number of nodes; ∑ represents the set of nodes connected to node r;

[0049] Then substitute the obtained parameters into the formula Calculate the node voltage deviation Ur of the substation grid node r grid , where η is the preset proportional coefficient, the specific value is set by the operator, Ur is the voltage amplitude of the r node; U N is the standard voltage amplitude, and the node voltage deviation Ur is calculated in order from small to large. grid Sorting generates optimization information.

[0050] Step S005: display the generated optimization information to the corresponding operator.

[0051] Some of the data in the above formulas are numerically calculated by removing their dimensions. Meanwhile, the contents not described in detail in this specification belong to the prior art known to those skilled in the art.

[0052] The above embodiments are only used to illustrate the technical method of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical method of the present invention may be modified or replaced by equivalents without departing from the spirit and scope of the technical method of the present invention.

Claims

1. A distributed photovoltaic energy storage site selection layout optimization method, characterized in that: The method specifically comprises the following steps: Step S001, obtaining geographic data of the study area based on GIS technology, and dividing the study area into multiple groups of areas to be analyzed; Step S002, comprehensively analyzing the geographic data of the area to be analyzed, and selecting a standard layout area by calculating the comprehensive area value of the area to be analyzed; Step S003, segmenting the obtained standard layout area to obtain multiple groups of segmented areas, and obtaining area information corresponding to the segmented areas, analyzing the area information and calculating the selected point values ​​of the segmented areas; Step S004, screening the segmented areas according to the calculated selected point values ​​of the segmented areas to obtain selected point areas, and performing optimization analysis based on node voltage deviations of the substations in the selected point areas to generate optimization information; Step S005: display the generated optimization information to the corresponding operator.

2. A distributed photovoltaic energy storage site selection and layout optimization method according to claim 1, characterized in that: The method of selecting the standard layout area in step S002 is: All the areas to be analyzed are obtained and labeled as i, and i=1, 2, ..., j, where j represents the number of the areas to be analyzed. Then, the terrain value Di, solar energy resource value Yi and site selection distance value Li corresponding to the area to be analyzed i are calculated respectively, and the obtained terrain value Di, solar energy resource value Yi and site selection distance value Li are substituted into the formula Qi=(Di+Yi+Li)×θ to calculate the comprehensive area value Qi of the area to be analyzed i, where θ is the preset proportional coefficient, and the specific value is set by the operator. Then, the area to be analyzed corresponding to the maximum comprehensive area value Qi is selected as the standard layout area.

3. A distributed photovoltaic energy storage site selection and layout optimization method according to claim 2, characterized in that: The terrain value in step S002 is calculated as follows: Based on GIS technology, the topography of the area i to be analyzed is obtained. At the same time, the proportion of the high mountain area in the area i to be analyzed is obtained according to the topography information, recorded as Ri, and the slope value of the area i to be analyzed is calculated and recorded as Pi. Then, the terrain undulation of the area i to be analyzed is obtained and recorded as Si. At the same time, the obtained parameters are substituted into the calculation formula The terrain value Di of the area to be analyzed i is calculated, where a1, a2 and a3 are weight coefficients, and a1+a2+a3=1.

4. A distributed photovoltaic energy storage site selection and layout optimization method according to claim 2, characterized in that: The calculation method of the solar energy resource value Yi in step S002 is: Get the average sunshine duration per unit time of the area to be analyzed, record it as Ti, and get the average sunshine weather Gi of the area to be analyzed based on historical data, then get the annual average solar radiation intensity corresponding to the area to be analyzed, record it as Fi, and substitute the obtained parameters into the formula The solar energy resource value Yi corresponding to the area to be analyzed is calculated, where b1, b2 and b3 are weight coefficients, and b1+b2+b3=1.

5. A distributed photovoltaic energy storage site selection and layout optimization method according to claim 2, characterized in that: The calculation method of the site selection distance value in step S002 is: Obtain all substations in the study area, and calculate the distance between the substation and the area to be analyzed i respectively, and calculate the mean of all distances as the site selection distance value Li of the area to be analyzed.

6. A distributed photovoltaic energy storage site selection and layout optimization method according to claim 1, characterized in that: The specific method of calculating the selected point value of the segmentation area in step S003 is: Obtain a standard layout area, and divide the standard layout area into multiple groups of divided areas and record them as o, where o=1, 2, ..., p, where p represents the number of divided areas, then obtain the residential density of the divided area o and record it as Jo, and obtain the contact line fluctuation rate Xo of the divided area o; Calculate the load demand matching value Mo of the segmented area o, obtain the load data of the segmented area o, and establish the corresponding load curve, and extract the peak load P from the load curve peak , valley load P valley , average load P avg and load factor LR, where the load factor is calculated according to the formula Calculate and substitute the obtained parameters into the formula The load demand matching degree Mo is calculated, where α is the preset proportional coefficient; Substitute the obtained parameters into the formula Calculate the selected point value ω corresponding to the segmented area o o , where μ is the preset proportional coefficient.

7. A distributed photovoltaic energy storage site selection and layout optimization method according to claim 1, characterized in that: The specific method of obtaining the selected point area in step S004 is: Sort the point selection values ​​of the segmentation areas from large to small, and select the segmentation area corresponding to the largest point selection value as the point selection area.

8. A distributed photovoltaic energy storage site selection and layout optimization method according to claim 1, characterized in that: The specific method of generating the optimization information in step S004 is: Get the substation in the selected area, then get all the grid nodes of the substation, denoted as r, and r = 1, 2, ..., h, where h represents the number of grid nodes, and calculate the active network loss corresponding to the substation. The specific calculation formula is: The active network loss P of the power grid node is calculated according to the formula loss , where P loss represents the active network loss of the system, Vr represents the voltage amplitude of node r, G rh , B rh , δ rh Respectively represent the conductance, susceptance and phase angle difference between nodes r and h, h is the total number of nodes; ∑ represents the set of nodes connected to node r; Then substitute the obtained parameters into the formula Calculate the node voltage deviation Ur of the substation grid node r grid , where η is the preset proportional coefficient, the specific value is set by the operator, Ur is the voltage amplitude of the r node; U N is the standard voltage amplitude, and the node voltage deviation Ur is calculated in order from small to large. grid Sorting generates optimization information.

Images