A method for investigating regional photovoltaic development potential based on drones

Abstract

The present invention relates to the field of photovoltaic power generation technology, and specifically to a method for investigating regional photovoltaic development potential based on drones, comprising the following steps: conducting drone aerial surveys, planning routes and deploying image control points; establishing a digital model of regional terrain; analyzing and modeling regional lighting conditions; determining a suitable scope for the construction of a photovoltaic power station in the region; and determining the potential photovoltaic resource development capacity of the region. The present invention solves the problems of single means, long cycle, and inaccuracy in investigating the potential photovoltaic development capacity in a large region, especially in a county, provincial or larger region, and achieves the goal of improving the quality and efficiency of the regional photovoltaic development potential investigation work. The drone technology is used to obtain real data of natural environment mapping with outstanding efficiency and advantages such as mobility, speed, digitization and intelligence, and improves the efficiency of the regional photovoltaic development potential investigation.

Description

Technical Field

[0001] The invention relates to the technical field of photovoltaic power generation, and in particular to a method for investigating regional photovoltaic development potential based on unmanned aerial vehicles. Background Art

[0002] Since the country announced the "3060 target", a new energy development boom has been set off, and mastering high-quality resources is one of the hot topics that people are most concerned about. For photovoltaic resources, on the one hand, everyone is concerned about solar energy resources, and on the other hand, everyone is concerned about how to minimize its development costs. The main factors affecting the technical cost are terrain conditions, access conditions, etc. Due to the advantages of high mobility, speed, digitization and intelligence of drone technology, it has outstanding efficiency and advantages in quickly obtaining real data for natural environment mapping. Therefore, the introduction of drone technology in the process of photovoltaic development potential investigation will greatly improve the efficiency and reliability of the investigation, and bring convenience to the exploration of high-quality resources. Summary of the invention

[0003] In order to achieve the goal of improving the quality and efficiency of regional photovoltaic development potential investigation work, and at the same time to solve the problems of single means, long cycle and inaccuracy in the investigation of potential photovoltaic development capacity in a vast area, especially in county, provincial or larger areas, the present invention provides a regional photovoltaic development potential investigation method based on drone.

[0004] To achieve the above object, the present invention adopts the following technical solutions:

[0005] A method for investigating regional photovoltaic development potential based on drones, comprising the following steps:

[0006] S1: UAV aerial survey, route planning and image control point layout;

[0007] S2: Establish a digital model of regional terrain;

[0008] S3: Regional lighting condition analysis and modeling;

[0009] S4: Determine the area suitable for the construction of photovoltaic power stations;

[0010] S5: Determine the potential photovoltaic resource development capacity in the region.

[0011] Furthermore, the step S2 comprises the following steps:

[0012] S201: Create a new project in the 3D real scene modeling software, add images, and submit the aerial triangulation calculation. Modify the relevant parameters according to the actual situation on site, and add image control points after the aerial triangulation calculation is completed;

[0013] S202: According to the image control points added in step S1, define their correct coordinate system and attributes and complete matching with the image;

[0014] S203: re-perform aerial triangulation based on steps S201 and S202, check the quality report and image control point accuracy, filter out image control points with large errors and delete them;

[0015] S204: Generate a digital three-dimensional model for the new reconstruction project;

[0016] S205: Extracting elevation information from the three-dimensional model to generate a digital elevation model.

[0017] Furthermore, step S3 includes the following steps:

[0018] Step S3 is calculated based on the horizontal coordinate system. According to the requirements of the "Design Specifications for Photovoltaic Power Stations" (GB50797-2012), the calculation time point is considered to be the winter solstice of the year. Therefore, the shadow of the object is the longest on this day, and the declination angle δ in the northern hemisphere on this day is -23.45. The specific calculation is as follows:

[0019] S301: Discretize the regional digital elevation model into N continuous surfaces (the degree of discretization can be adjusted according to the design depth), each surface center point corresponds to an elevation Hi (0≤i≤N), each surface center point corresponds to a latitude Each face corresponds to an azimuth νi (0≤i≤N), each face corresponds to a slope θi (0≤i≤N), and each face corresponds to an area Si (0≤i≤N);

[0020] S302: Set the area shadow non-blocking time in the morning to T am , the unblocked time in the afternoon is T pm , the solar hour angles corresponding to the unobstructed time in the morning and the unobstructed time in the afternoon are ω am =15°×(T am -12),ω pm =15°×(T pm -12);

[0021] S303: Calculate the solar altitude angle ɑ during the morning period am And the solar altitude angle ɑ in the afternoon pm , the expressions are:

[0022]

[0023]

[0024] Calculate the solar azimuth angle β during the morning period am and the solar azimuth angle β in the afternoon pm , the expressions are:

[0025]

[0026]

[0027] S304: Obtain the shadow length ζ of each surface vertically facing north corresponding to the digital elevation model i , the digital elevation model corresponds to the shadow length of each face in the morning to the west ζ W和 The digital elevation model corresponds to the shadow length of each surface in the afternoon to the east ζ E , the expressions are:

[0028] ζi=Hi×cos(ɑ am -νi)(0≤i≤N);

[0029] ζ W =Hi×sin(ɑ am -νi) / tan(ɑ am )(0≤i≤N);

[0030] ζ E =Hi×sin(ɑ pm -νi) / tan(ɑ pm )(0≤i≤N);

[0031] S305: By ζi, ζ W and E The composed surface is the shadow surface Ssdi (0≤i≤N) at any point of the elevation model on the winter solstice. Then the corresponding shadow surface of the elevation model on the winter solstice is Ssd=∑Ssdi (0≤i≤N).

[0032] Furthermore, the step S4 comprises the following steps:

[0033] S401: Azimuth ν of any surface of the elevation model i (0≤i≤N) and slope θ i (0≤i≤N) can calculate the slope θ in the north-south direction at any point in the elevation model SNi =tan -1 (tanθ i ÷cos∣ν i ∣)(0≤i≤N) and the slope θ in the east-west direction at any point in the elevation model EWi =tan -1 (tanθ i ÷sin∣ν i ∣)(0≤i≤N);

[0034] S402: Considering the construction cost and power generation efficiency of the photovoltaic power station, the area with too large azimuth angle or large slope is identified as an unsuitable development area. The azimuth angle and slope here can be formulated in combination with the actual situation of the project. Generally speaking, the azimuth angle exceeds ±8° and the slope exceeds 35°. The area is identified as an unsuitable development area Suni (0≤i≤N) in the present invention. The area of ​​the unsuitable development area is calculated as Sun=∑Suni(0≤i≤N);

[0035] S403: The elevation model area of ​​the entire planning area is S = ∑Si (00≤i≤N). Combined with steps S3, S401 and S402, the shadow surface Ssd corresponding to the elevation model on the winter solstice and the area Sun that is not suitable for development are calculated, and the available photovoltaic power station area in the area is Sav = S-Ssd-Sun, which determines the area suitable for photovoltaic power station construction in the region.

[0036] Furthermore, the step S5 comprises the following steps:

[0037] S501: Based on step S3, the elevation Hi (0≤i≤N) of any center point of a surface in the region is calculated. Each center point of a surface corresponds to a latitude. Each face corresponds to a north-south slope θ SNi =tan -1 (tanθ i ÷cos∣ν i ∣) (0≤i≤N) calculation results, divide the available area of ​​the planning area into Class I terrain area, Class II terrain area and Class III terrain area according to elevation and slope;

[0038] S502: Assume that there are J usable areas in the terrain area of ​​type I, and the area of ​​the terrain area of ​​type I is S Ⅰj (0≤j≤J); there are M usable areas in the Class II terrain area, and the area of ​​the Class II terrain area is S Ⅱm (0≤m≤M); There are N usable areas in the Class III terrain area, and the area of ​​the Class III terrain area is S Ⅲn (0≤n≤N);

[0039] S503: Calculate the representative latitude value of the planning area Given the adopted component efficiency η, the planned voltage level of the photovoltaic power station, and the planned support form of the photovoltaic power station, the overall land use index S of each 10MW photovoltaic power station under three types of terrain conditions and various types of support installation types is calculated. 1ac , S 2ac , S 3ac , the present invention 1ac , S 2ac , S 3acThe calculation adopts the software "Design and Calculation Software for Reasonable Floor Space of Photovoltaic Power Station" (Software Registration No.: 2019SR0058318) independently developed by our institute;

[0040] S504: Based on the available area S of three types of terrain areas Ⅰj , S Ⅱm , S Ⅲn Calculate the photovoltaic power station capacity R in the usable area of ​​Class I terrain area Ⅰj , the photovoltaic power station capacity R in the area of ​​Class II terrain area can be used Ⅱm and the photovoltaic power station capacity R in the Class III terrain area Ⅲn , the calculation formulas are:

[0041]

[0042]

[0043]

[0044] S505: Calculate the regional potential photovoltaic resource development capacity R, the calculation formula is:

[0045]

[0046] Furthermore, the basis for the division of the planning area in step S501 is as follows: Class I terrain area refers to a plain area with no obvious undulations and a natural ground slope of less than or equal to 3°; Class II terrain area refers to a micro-hill area with little terrain undulations, a natural ground slope greater than 3° but less than or equal to 20°, and a relative height difference of less than 200m; Class III terrain area refers to a heavy hill or mountain area with large terrain undulations, a natural ground slope greater than 20°, and a relative height difference of more than 200m.

[0047] Compared with the prior art, the survey method of the present invention has the following advantages:

[0048] The present invention solves the problems of single means, long cycle and inaccuracy in investigating potential photovoltaic development capacity in a vast area, especially in a county, provincial or larger area.

[0049] The present invention achieves the goal of improving the quality and efficiency of regional photovoltaic development potential survey work, and uses drone technology to obtain natural environment mapping real data with outstanding efficiency and advantages such as mobility, speed, digitization and intelligence, which improves the efficiency of regional photovoltaic development potential survey. In addition, through a reasonable model, the regional boundary conditions (geographic location, terrain, layout plan, etc.) are objectively described, which improves the reliability of the regional photovoltaic development potential survey results, and brings convenience to the exploration of high-quality resources on the basis of making solar renewable energy planning meet objective and reasonable requirements.

[0050] The method of the invention is novel and highly operational. In view of the differences in environmental conditions in different regions, combined with UAV technology and a reasonable mathematical calculation model, a technical route for regional photovoltaic development potential investigation is proposed, aiming to achieve the goal of improving the quality and efficiency of regional photovoltaic development potential investigation work, which has important guiding significance for the development of solar energy projects in various regions. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 It is a flow chart of the regional photovoltaic development potential investigation method based on drone of the present invention;

[0052] Figure 2 This is a schematic diagram of the planning scope of a photovoltaic power station in Embodiment 2 of the present invention;

[0053] Figure 3 This is a schematic diagram of the available area range in Example 2 of the present invention;

[0054] Figure 4 It is a schematic diagram of the usable area range under different terrain conditions in Example 2 of the present invention. DETAILED DESCRIPTION

[0055] In order to enable those skilled in the art to better understand the technical solution of the present invention, the present invention is further described in detail below in conjunction with the accompanying drawings and specific embodiments. The following embodiments are only used to illustrate the present invention but are not intended to limit the scope of the present invention.

[0056] Example 1

[0057] This embodiment provides a method for investigating regional photovoltaic development potential based on a drone, comprising the following steps:

[0058] S1: UAV aerial survey, route planning and image control point layout;

[0059] S2: Establish a digital model of regional terrain. The specific steps are as follows:

[0060] S201: Create a new project in the 3D real scene modeling software, add images, and submit the aerial triangulation calculation. Modify the relevant parameters according to the actual situation on site, and add image control points after the aerial triangulation calculation is completed;

[0061] S202: According to the image control points added in step S1, define their correct coordinate system and attributes and complete matching with the image;

[0062] S203: re-perform aerial triangulation based on steps S201 and S202, check the quality report and image control point accuracy, filter out image control points with large errors and delete them;

[0063] S204: Generate a digital three-dimensional model for the new reconstruction project;

[0064] S205: Extracting elevation information from the three-dimensional model to generate a digital elevation model.

[0065] S3: Regional lighting condition analysis and modeling. The specific steps are as follows:

[0066] Step S3 is calculated based on the horizontal coordinate system. According to the requirements of the "Design Specifications for Photovoltaic Power Stations" (GB50797-2012), the calculation time point is considered to be the winter solstice of the year. Therefore, the shadow of the object is the longest on this day, and the declination angle δ in the northern hemisphere on this day is -23.45. The specific calculation is as follows:

[0067] S301: Discretize the regional digital elevation model into N continuous surfaces (the degree of discretization can be adjusted according to the design depth), each surface center point corresponds to an elevation Hi (0≤i≤N), each surface center point corresponds to a latitude Each face corresponds to an azimuth νi (0≤i≤N), each face corresponds to a slope θi (0≤i≤N), and each face corresponds to an area Si (0≤i≤N);

[0068] S302: Set the area shadow non-blocking time in the morning to T am , the unblocked time in the afternoon is T pm , the solar hour angles corresponding to the unobstructed time in the morning and the unobstructed time in the afternoon are ω am =15°×(T am -12),ω pm =15°×(T pm -12);

[0069] S303: Calculate the solar altitude angle ɑ during the morning period am And the solar altitude angle ɑ in the afternoon pm , the expressions are:

[0070]

[0071]

[0072] Calculate the solar azimuth angle β during the morning period am and the solar azimuth angle β in the afternoon pm , the expressions are:

[0073]

[0074]

[0075] S304: Obtain the shadow length ζ of each surface vertically facing north corresponding to the digital elevation model i , the digital elevation model corresponds to the shadow length of each face in the morning to the west ζ W和The digital elevation model corresponds to the shadow length of each surface in the afternoon to the east ζ E , the expressions are:

[0076] ζi=Hi×cos(ɑ am -νi)(0≤i≤N);

[0077] ζ W =Hi×sin(ɑ am -νi) / tan(ɑ am )(0≤i≤N);

[0078] ζ E =Hi×sin(ɑ pm -νi) / tan(ɑ pm )(0≤i≤N);

[0079] S305: By ζi, ζ W and E The composed surface is the shadow surface Ssdi (0≤i≤N) at any point of the elevation model on the winter solstice. Then the corresponding shadow surface of the elevation model on the winter solstice is Ssd=∑Ssdi (0≤i≤N).

[0080] S4: Determine the area suitable for the construction of photovoltaic power stations. The specific steps are as follows:

[0081] S401: Azimuth ν of any surface of the elevation model i (0≤i≤N) and slope θ i (0≤i≤N) can calculate the slope θ in the north-south direction at any point in the elevation model SNi =tan -1 (tanθ i ÷cos∣ν i ∣)(0≤i≤N) and the slope θ in the east-west direction at any point in the elevation model EWi =tan -1 (tanθ i ÷sin∣ν i ∣)(0≤i≤N);

[0082] S402: Considering the construction cost and power generation efficiency of the photovoltaic power station, the area with too large azimuth angle or large slope is identified as an unsuitable development area. The azimuth angle and slope here can be formulated in combination with the actual situation of the project. Generally speaking, the azimuth angle exceeds ±8° and the slope exceeds 35°. The area is identified as an unsuitable development area Suni (0≤i≤N) in the present invention. The area of ​​the unsuitable development area is calculated as Sun=∑Suni(0≤i≤N);

[0083] S403: The elevation model area of ​​the entire planning area is S = ∑Si (00≤i≤N). Combined with steps S3, S401 and S402, the shadow surface Ssd corresponding to the elevation model on the winter solstice and the area Sun that is not suitable for development are calculated, and the available photovoltaic power station area in the area is Sav = S-Ssd-Sun, which determines the area suitable for photovoltaic power station construction in the region.

[0084] S5: Determine the potential photovoltaic resource development capacity of the region. The specific steps are as follows:

[0085] S501: Based on step S3, the elevation Hi (0≤i≤N) of any center point of a surface in the region is calculated. Each center point of a surface corresponds to a latitude. Each face corresponds to a north-south slope θ SNi =tan -1 (tanθ i ÷cos∣ν i ∣) (0≤i≤N) calculation results, according to the elevation and slope, the available area of ​​the planning area is divided into Class I terrain area, Class II terrain area and Class III terrain area. Class I terrain area refers to the plain area with no obvious terrain undulation and the natural ground slope is less than or equal to 3°; Class II terrain area refers to the micro-hill area with small terrain undulation, the natural ground slope is greater than 3° but less than or equal to 20°, and the relative height difference is within 200m; Class III terrain area refers to the heavy hill or mountain area with large terrain undulation, the natural ground slope is greater than 20°, and the relative height difference is more than 200m;

[0086] S502: Assume that there are J usable areas in the terrain area of ​​type I, and the area of ​​the terrain area of ​​type I is S Ⅰj (0≤j≤J); there are M usable areas in the Class II terrain area, and the area of ​​the Class II terrain area is S Ⅱm (0≤m≤M); There are N usable areas in the Class III terrain area, and the area of ​​the Class III terrain area is S Ⅲn (0≤n≤N);

[0087] S503: Calculate the representative latitude value of the planning area Given the adopted component efficiency η, the planned voltage level of the photovoltaic power station, and the planned support form of the photovoltaic power station, the overall land use index S of each 10MW photovoltaic power station under three types of terrain conditions and various types of support installation types is calculated. 1ac , S 2ac , S 3ac , the present invention 1ac , S 2ac , S 3ac The calculation adopts the software "Design and Calculation Software for Reasonable Floor Space of Photovoltaic Power Station" (Software Registration No.: 2019SR0058318) independently developed by our institute;

[0088] S504: Based on the available area S of three types of terrain areas Ⅰj , S Ⅱm , S Ⅲn Calculate the photovoltaic power station capacity R in the usable area of ​​Class I terrain area Ⅰj , the photovoltaic power station capacity R in the area of ​​Class II terrain area can be used Ⅱm and the photovoltaic power station capacity R in the Class III terrain area Ⅲn , the calculation formulas are:

[0089]

[0090]

[0091]

[0092] S505: Calculate the regional potential photovoltaic resource development capacity R, the calculation formula is:

[0093]

[0094] Example 2

[0095] Figure 2 This is a schematic diagram of the planning range of a photovoltaic power station. Taking this as an example, the method for investigating regional photovoltaic development potential based on a drone provided by the present invention includes:

[0096] (1) UAV aerial survey, route planning and image control point layout, determine the planning scope of the photovoltaic power station, such as Figure 2 As shown. Use drone technology to form a digital elevation model within the planning area;

[0097] (2) According to the specific methods of step S3 and step S4 of the present invention, the available area within the planning range is analyzed to determine the area suitable for the construction of a photovoltaic power station, such as Figure 3 ;

[0098] (3) According to step S5 of the present invention, the land classification of the available area within the planning range is analyzed and obtained, such as Figure 4 ;

[0099] (4) Assume that the representative latitude of the planning scope in this planning scheme is 35.8°, the photovoltaic bracket installation method is fixed, the photovoltaic module efficiency is 22%, and the 330kV collection station scheme is adopted.

[0100] Based on the overall indicators of land use for fixed 10MW photovoltaic power stations under different terrain conditions in Table 1-3, the developable photovoltaic capacity within the scope of this planning is obtained, as shown in Table 4.

[0101] Table 1-General indicators of land use for fixed 10MW photovoltaic power stations under Class I terrain conditions

[0102] Unit: 10MW / hectare

[0103]

[0104]

[0105]

[0106] Table 2-General indicators of land use for fixed 10MW photovoltaic power stations under Class II terrain conditions

[0107] Unit: 10MW / hectare

[0108]

[0109]

[0110]

[0111] Table 3- Overall indicators of land use for fixed 10MW photovoltaic power stations under Class III terrain conditions

[0112] Unit: 10MW / hectare

[0113]

[0114]

[0115]

[0116] Table 4 Planned capacity of photovoltaic power station in Example 2

[0117]

[0118] The photovoltaic bracket installation method of the present invention can adopt other methods. If the photovoltaic bracket installation method adopts the flat single-axis tracking method, the photovoltaic developable capacity under the corresponding planning scheme within the planning scope can be calculated according to Table 5-7.

[0119] Table 5-1 General indicators of land use for flat single-axis tracking 10MW photovoltaic power station under Class I terrain conditions

[0120] Unit: 10MW / hectare

[0121]

[0122]

[0123]

[0124] Table 6-General indicators of land use for flat single-axis tracking 10MW photovoltaic power station under Class II terrain conditions

[0125] Unit: 10MW / hectare

[0126]

[0127]

[0128]

[0129] Table 7- Overall indicators of land use for flat single-axis tracking 10MW photovoltaic power station under Class III terrain conditions

[0130] Unit: 10MW / hectare

[0131]

[0132]

[0133]

[0134] If the photovoltaic bracket is installed in an inclined single-axis tracking manner, the photovoltaic developable capacity under the corresponding planning scheme within the planning scope can be calculated according to Table 8-10, as follows:

[0135] Table 8- Overall indicators of land use for inclined single-axis tracking 10MW photovoltaic power station under Class I terrain conditions

[0136] Unit: 10MW / hectare

[0137]

[0138]

[0139]

[0140] Table 9-Overall indicators of land use for inclined single-axis tracking 10MW photovoltaic power station under Class II terrain conditions

[0141] Unit: 10MW / hectare

[0142]

[0143]

[0144]

[0145] Table 10- Overall indicators of land use for inclined single-axis tracking 10MW photovoltaic power station under Class III terrain conditions

[0146] Unit: 10MW / hectare

[0147]

[0148]

[0149]

[0150] The preferred embodiments of the present invention are described in detail above. However, the present invention is not limited to the specific details in the above embodiments. Within the technical concept of the present invention, various changes can be made to the technical scheme of the present invention. These simple changes all fall within the protection scope of the present invention.

Claims

1. A method for investigating regional photovoltaic development potential based on drones, characterized in that: The following steps are involved: S1: UAV aerial survey, route planning and image control point layout; S2: Establish a digital model of regional terrain: S201: Create a new project in the 3D real scene modeling software, add images, and submit the aerial triangulation calculation. Modify the relevant parameters according to the actual situation on site. After the aerial triangulation calculation is completed, add image control points. S202: Based on the image control points added in step S1, define their correct coordinate system and attributes and match them with the image. S203: Re-perform aerial triangulation based on steps S201 and S202, check the quality report and image control point accuracy, filter out image control points with large errors and delete them. S204: Generate a digital 3D model for the new reconstruction project. S205: extracting elevation information from the three-dimensional model and generating a digital elevation model; S3: Regional lighting condition analysis modeling: S301: Discretize the regional digital elevation model into N continuous surfaces, each surface center point corresponds to an elevation Hi (0≤i≤N), each surface center point corresponds to a latitude φi (0≤i≤N), each surface corresponds to an azimuth νi (0≤i≤N), each surface corresponds to a slope θi (0≤i≤N), each surface corresponds to an area Si (0≤i≤N), S302: Set the area shadow non-blocking time in the morning to T am , the unblocked time in the afternoon is T pm , the solar hour angles corresponding to the unobstructed time in the morning and the unobstructed time in the afternoon are ω am =15°×(T am -12),ω pm =15°×(T pm -12), S303: Calculate the solar altitude angle ɑ during the morning period am and the solar altitude angle ɑ in the afternoon pm， The expressions are: ɑ am =arcsin(sinφ·sinδ+cosφi·cosδ·cosω am )(0≤i≤N), ɑ pm =arcsin(sinφ·sinδ+cosφi·cosδ·cosω pm )(0≤i≤N), Calculate the solar azimuth angle β during the morning period am and the solar azimuth angle β in the afternoon pm , the expressions are: b am =arcos[(sinɑ am ·sinφi-sinδ) / cosɑ am ·cosφi](0≤i≤N), βpm=arcos[(sinɑ pm ·sinφi-sinδ) / cosɑ pm ·cosφi](0≤i≤N), S304: Obtain the shadow length ζ of each surface vertically facing north corresponding to the digital elevation model i , the digital elevation model corresponds to the shadow length of each face in the morning to the west ζ W和 The digital elevation model corresponds to the shadow length of each surface in the afternoon to the east ζ E， The expressions are: ζi=Hi×cos(ɑ am -νi)(0≤i≤N), g W =Hi×sin(ɑ am -νi) / tan(ɑ am )(0≤i≤N), g E =Hi×sin(ɑ pm -νi) / tan(ɑ pm )(0≤i≤N), S305: By ζi, ζ W and E The composed surface is the shadow surface Ssdi (0≤i≤N) at any point of the elevation model on the winter solstice. Then the shadow surface corresponding to the elevation model on the winter solstice is Ssd = ∑Ssdi (0≤i≤N). S4: Determine the scope of suitable photovoltaic power station construction in the region: S401: Azimuth ν of any surface of the elevation model i (0≤i≤N) and slope θ i (0≤i≤N) can calculate the slope θ in the north-south direction at any point in the elevation model SNi =tan -1 (tanθ i ÷cos∣ν i ∣)(0≤i≤N) and the slope θ in the east-west direction at any point in the elevation model EWi =tan -1 (tanθ i ÷sin∣ν i ∣)(0≤i≤N), S402: Considering the construction cost and power generation efficiency of photovoltaic power stations, the area with an azimuth angle exceeding ±8° and a slope exceeding 35° is identified as an unsuitable development area Suni (0≤i≤N), and the area of ​​the unsuitable development area is calculated as Sun=∑Suni(0≤i≤N). S403: The elevation model area of ​​the entire planning area is S=∑Si(00≤i≤N). Combined with steps S3, S401 and S402, the shadow surface Ssd corresponding to the elevation model on the winter solstice and the unsuitable development area Sun are calculated, and the available photovoltaic power station area in the area is Sav=S-Ssd-Sun, and the suitable photovoltaic power station construction range in the area is determined; S5: Determine the potential photovoltaic resource development capacity in the region: S501: Based on step S3, the elevation Hi (0≤i≤N) of the center point of any surface in the region is calculated. Each center point of the surface corresponds to a latitude φi (0≤i≤N), and each surface corresponds to a slope θ in the north-south direction. SNi =tan -1 (tanθ i ÷cos∣ν i ∣)(0≤i≤N) calculation results, according to the elevation and slope, the available area of ​​the planning area is divided into Class I terrain area, Class II terrain area and Class III terrain area. Class I terrain area refers to the plain area with no obvious terrain undulation and the natural ground slope is less than or equal to 3°. Class II terrain area refers to the micro-hill area with small terrain undulation, the natural ground slope is greater than 3° but less than or equal to 20°, and the relative height difference is within 200m. Class III terrain area refers to the heavy hill or mountain area with large terrain undulation, the natural ground slope is greater than 20°, and the relative height difference is more than 200m. S502: Assume that there are J usable areas in the terrain area of ​​type I, and the area of ​​the terrain area of ​​type I is S Ⅰj (0≤j≤J); there are M usable areas in the Class II terrain area, and the area of ​​the Class II terrain area is S Ⅱm (0≤m≤M); there are N usable areas in the Class III terrain area, and the area of ​​the Class III terrain area is S Ⅲn (0≤n≤N), S503: Calculate the representative latitude value φ=max(φi)(0≤i≤N) of the planning area. Given the adopted component efficiency η, the voltage level of the planned photovoltaic power station, and the support form of the planned photovoltaic power station, calculate the overall land use index S of each 10MW photovoltaic power station under the three types of terrain conditions and various types of support installation types. 1ac , S 2ac , S 3ac , S504: Based on the available area S of three types of terrain areas Ⅰj , S Ⅱm、 S Ⅲn， Calculate the photovoltaic power station capacity R in the usable area of ​​Class I terrain area Ⅰj , the photovoltaic power station capacity R in the area of ​​Class II terrain area can be used Ⅱm and the photovoltaic power station capacity R in the Class III terrain area Ⅲn， The calculation formulas are: S505: Calculate the regional potential photovoltaic resource development capacity R, the calculation formula is:

Images