A method and device for quantifying the effect of mountain terrain on solar radiation shielding

Abstract

The present invention discloses a method and device for quantifying the shielding effect of mountain terrain on solar radiation, and belongs to the field of terrain shielding data processing technology. The method comprises: obtaining the skyline and visible field elevation angle of the target research point under the influence of terrain shielding, calculating the solar altitude angle and azimuth; calculating the first direct radiation intensity and the first scattered radiation intensity of any surface per minute when the shielding effect is not considered; calculating the second direct radiation intensity and the second scattered radiation intensity of any surface per minute under the influence of terrain shielding; obtaining the solar direct radiation and scattered radiation of any surface under the influence of mountain terrain according to the second direct radiation intensity and the second scattered radiation intensity, and then combining the direct radiation and scattered radiation of any surface when the shielding effect is not considered, calculating the direct radiation shielding rate and scattered radiation shielding rate of any surface under the influence of mountain terrain. The present invention can quickly predict the shielding effect of mountain terrain on solar radiation, and improve the accuracy of shielding data prediction.

Description

Technical Field

[0001] The present invention relates to the technical field of terrain shielding data processing, and in particular to a method and device for quantifying the shielding effect of mountain terrain on solar radiation. Background Art

[0002] Mountainous areas are rich in solar energy. Fully developing and utilizing mountain solar energy can not only solve the winter heating problem in mountainous areas and improve the living environment, but also reduce carbon emissions and energy waste. Maximizing the utilization of solar energy in mountainous areas and accurately predicting the amount of solar radiation received by each surface are the basic prerequisites. However, the terrain in mountainous areas is complex and changeable, and the solar radiation received by different planes is affected by different degrees of terrain shielding, and the shielding effect is very significant. Therefore, in order to accurately predict the amount of solar radiation in mountainous areas, it is first necessary to quantify the degree of shielding effect of mountainous terrain on the solar radiation received by each plane.

[0003] In mountainous areas with complex and changeable terrain, there are few meteorological stations that can monitor solar radiation data, resulting in a lack of solar radiation data in most mountainous areas; most existing solar radiation models are based on unobstructed plain areas, without considering the shielding effect of mountainous terrain on solar radiation. It is difficult to accurately quantify and analyze the shielding effect of mountainous terrain on solar radiation, which restricts the effective utilization of solar energy resources in mountainous areas. In addition, in mountainous urban environments, it is also very important to quantify the impact of terrain on solar radiation shielding on any surface for the optimization of solar equipment layout and installation in mountainous towns or the solar energy application design of building facades in various directions. Although GIS (Geographic Information System)-based research based on DEM (Digital Elevation Model) can be applied to the study of solar radiation under mountainous terrain shielding, these existing studies only consider the shielding effect of mountainous terrain on horizontal solar radiation, and cannot quantify the degree of shielding effect of mountainous terrain on solar radiation on any surface. Summary of the invention

[0004] In order to overcome the technical problems existing in the related art at least to a certain extent, the present invention provides a method and device for quantifying the shielding effect of mountainous terrain on solar radiation.

[0005] In a first aspect, an embodiment of the present invention provides a method for quantifying the effect of mountain terrain on solar radiation shielding, comprising the following steps:

[0006] Based on the digital elevation model data, the skyline and visible field elevation angle of the target research point under the influence of terrain occlusion are obtained;

[0007] Analyze the sun's motion trajectory according to the latitude, longitude and altitude of the target research point, and calculate the sun's altitude angle and azimuth angle based on the full-year time dimension with a calculation step of minutes;

[0008] By using the solar altitude angle and azimuth as well as the inclination and orientation of the target research surface, the first direct radiation intensity and the first scattered radiation intensity of any surface are calculated every minute without considering the influence of shielding;

[0009] Calculate the second direct radiation intensity and the second scattered radiation intensity of any surface per minute under the influence of terrain shielding according to the visible field elevation angle, the solar altitude angle and the azimuth angle, and the inclination and orientation of the target research surface;

[0010] By integrating the second direct radiation intensity and the second scattered radiation intensity over time, the direct solar radiation and scattered radiation of any surface under the influence of the mountainous terrain within any time interval are obtained. Then, combined with the direct radiation and scattered radiation of any surface without considering the shielding effect, the direct radiation shielding rate and scattered radiation shielding rate of any surface under the influence of the mountainous terrain within any time interval are calculated.

[0011] Preferably, the step of obtaining the skyline and visible field elevation angle of the target research point under the influence of terrain occlusion based on the digital elevation model data specifically includes:

[0012] Converting the digital elevation model data into an irregular triangular grid, and determining the specific location of the target research point in the irregular triangular grid according to the longitude and latitude of the target research point;

[0013] Obtaining the skyline of the target research point under the influence of terrain occlusion;

[0014] Determine the visible field of view of the horizontal plane research point under the influence of terrain occlusion and obtain the visible field elevation angle on the visible field boundary corresponding to each degree of azimuth.

[0015] Preferably, the step of calculating the first direct radiation intensity and the first scattered radiation intensity of any surface minute by minute without considering the influence of shielding specifically includes:

[0016] The ASHRAE and Perez model algorithms are used to calculate the first direct radiation intensity and the first scattered radiation intensity of any surface every minute without shielding.

[0017] Preferably, the step of calculating the second direct radiation intensity and the second scattered radiation intensity of any surface minute by minute under the influence of terrain occlusion specifically includes:

[0018] Based on the sun's motion trajectory and ray tracing method, the visible field is used to calculate the second direct radiation intensity and the second scattered radiation intensity of any surface minute by minute under the influence of terrain occlusion.

[0019] In a second aspect, an embodiment of the present invention provides a device for quantifying the effect of mountain terrain on solar radiation shielding, including:

[0020] A data acquisition module is used to obtain the skyline and visible field elevation angle of the target research point under the influence of terrain occlusion based on the digital elevation model data;

[0021] An angle calculation module, used to analyze the sun's motion trajectory according to the latitude and longitude and altitude of the target research point, and calculate the sun's altitude angle and azimuth angle based on the full-year time dimension and in minutes as the calculation step;

[0022] An unobstructed radiation intensity calculation module is used to calculate the first direct radiation intensity and the first scattered radiation intensity of any surface per minute without considering the influence of obstruction by using the solar altitude angle and azimuth as well as the inclination and orientation of the target research surface;

[0023] There is an obstruction radiation intensity calculation module, which is used to calculate the second direct radiation intensity and the second scattered radiation intensity of any surface minute by minute under the influence of terrain obstruction according to the visible field elevation angle, solar altitude angle and azimuth, and the inclination and orientation of the target research surface;

[0024] The radiation shielding rate calculation module is used to obtain the solar direct radiation and scattered radiation of any surface under the influence of mountainous terrain within any time interval by integrating the second direct radiation intensity and the second scattered radiation intensity over time, and then calculate the direct radiation shielding rate and scattered radiation shielding rate of any surface under the influence of mountainous terrain within any time interval by combining the direct radiation and scattered radiation of any surface without considering the shielding effect.

[0025] Preferably, the data acquisition module is specifically used for:

[0026] Converting the digital elevation model data into an irregular triangular grid, and determining the specific location of the target research point in the irregular triangular grid according to the longitude and latitude of the target research point;

[0027] Obtaining the skyline of the target research point under the influence of terrain occlusion;

[0028] Determine the visible field of view of the horizontal plane research point under the influence of terrain occlusion and obtain the visible field elevation angle on the visible field boundary corresponding to each degree of azimuth.

[0029] Preferably, the unobstructed radiation intensity calculation module is specifically used to:

[0030] The ASHRAE and Perez model algorithms are used to calculate the first direct radiation intensity and the first scattered radiation intensity of any surface every minute without shielding.

[0031] Preferably, the shielded radiation intensity calculation module is specifically used for:

[0032] Based on the sun's motion trajectory and ray tracing method, the visible field is used to calculate the second direct radiation intensity and the second scattered radiation intensity of any surface minute by minute under the influence of terrain occlusion.

[0033] The technical solution provided by the embodiments of the present invention has the following beneficial effects:

[0034] The embodiment of the present invention first obtains the skyline and visible field elevation of the target research point under the influence of terrain shielding, analyzes the movement trajectory of the sun position according to the longitude and latitude and the altitude, and calculates the solar altitude angle and azimuth; then calculates the first direct radiation intensity and the first scattered radiation intensity of any surface every minute when the shielding effect is not considered, and calculates the second direct radiation intensity and the second scattered radiation intensity of any surface every minute when the terrain shielding effect is considered; then performs an integral operation to obtain the solar direct radiation and scattered radiation of any surface under the influence of mountainous terrain in any time interval; finally, the direct radiation and scattered radiation of any surface when the shielding effect is not considered are combined to calculate the direct radiation shielding rate and scattered radiation shielding rate of any surface under the influence of mountainous terrain in any time interval, thereby quantifying the solar radiation data under the shielding of complex terrain in mountainous areas. The embodiment of the present invention can be used for the correction calculation of solar radiation data of any surface under the influence of mountainous terrain, and effectively improves the accuracy of the prediction of solar radiation shielding data of any receiving surface by mountainous terrain.

[0035] It should be understood that the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present application and, together with the description, serve to explain the principles of the present application.

[0037] Figure 1 It is a flow chart of a method for quantifying the effect of mountain terrain on solar radiation shielding provided by an embodiment of the present invention.

[0038] Figure 2 It is a block diagram of a device for quantifying the shielding effect of mountainous terrain on solar radiation provided by an embodiment of the present invention.

[0039] The accompanying drawings in the figure are marked as follows:

[0040] Data acquisition module 100; Angle calculation module 200; Unobstructed radiation intensity calculation module 300;

[0041] There is a shielding radiation intensity calculation module 400 and a radiation shielding rate calculation module 500. DETAILED DESCRIPTION

[0042] Exemplary embodiments will be described in detail herein, examples of which are shown in the accompanying drawings. When the following description refers to the drawings, unless otherwise indicated, the same numbers in different drawings represent the same or similar elements. The implementations described in the following exemplary embodiments do not represent all implementations consistent with the present application. Instead, they are merely examples of methods and devices consistent with some aspects of the present application as detailed in the appended claims.

[0043] In order to make the objectives, technical solutions and advantages of the present invention more clear, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings.

[0044] Example 1

[0045] like Figure 1 FIG. 1 is a flow chart of a method for quantifying the effect of mountain terrain on solar radiation shielding provided by an embodiment of the present invention. The method may include the following steps:

[0046] Step S100: Based on the digital elevation model data, the skyline and visible field elevation angle of the target research point under the influence of terrain occlusion are obtained.

[0047] Among them, the Digital Elevation Model (DEM) is a physical ground model that represents the ground elevation in the form of a set of ordered numerical arrays. Preferably, the DEM data used comes from the ASTER GDEM V2 version data in the geospatial data cloud platform of the Computer Network Information Center of the Chinese Academy of Sciences. It is currently the only high-resolution elevation image data covering the global land surface and has high accuracy.

[0048] In this embodiment, the digital elevation model data can be processed in the ESRI ArcGIS platform.

[0049] During implementation, the digital elevation model data is first converted into an irregular triangular grid to accurately express the mutation phenomenon of complex terrain in mountainous areas and improve the calculation accuracy. Then, the specific location of the target research point is determined in the irregular triangular grid according to the longitude and latitude of the target research point. Next, the skyline tool of the spatial analyst module in ESRI ArcGIS is used to obtain the skyline of the target research point under the influence of terrain occlusion. Finally, the viewshed analysis tool is used to perform a detailed analysis of the skyline to determine the visible viewshed of the horizontal plane research point under terrain occlusion, and to obtain the elevation data on the visible viewshed boundary corresponding to each degree of azimuth. The above data are easier to obtain and the database is more complete.

[0050] In another embodiment, fisheye camera shooting and image processing technology may also be used to obtain skyline and elevation angle data.

[0051] Step S200: Analyze the sun's motion trajectory according to the latitude, longitude and altitude of the target research point, and calculate the sun's altitude angle and azimuth angle based on the annual time dimension with minutes as the calculation step.

[0052] Specifically, the position of the sun in the sky can usually be expressed as the solar altitude angle ( β sun ) and azimuth ( ). The solar altitude angle is defined as the angle between the direct sun beam and the horizontal plane. The solar altitude angle is 0° at sunrise and sunset and reaches its maximum at noon. The solar azimuth angle is the angle between the projection of the direct sun beam on the horizontal plane and due south, which is negative in the morning, positive in the afternoon, and 0° at noon.

[0053] The solar azimuth and altitude are affected by the local latitude ( ), solar declination angle (δ) and solar hour angle ( ). The calculation formula of the solar altitude angle (βsun) is as follows:

[0054] Solar azimuth ( ) can be expressed by the following formula: Step S300: Calculate the first direct radiation intensity and the first scattered radiation intensity of any surface per minute without considering the effect of shielding by using the solar altitude angle and azimuth angle as well as the inclination and orientation of the target research surface.

[0055] In this embodiment, the ASHRAE and Perez model algorithms can be used in the MATLAB environment to calculate the first direct radiation intensity and the first scattered radiation intensity of any surface minute by minute in the absence of shielding.

[0056] Specifically, the normal direct radiation intensity (W / m 2 ) and horizontal diffuse radiation intensity (W / m 2 ), the calculation formula is as follows: In the formula, I dn is the unobstructed normal direct radiation intensity, in w / m 2 ; I hd is the scattered radiation intensity of the unobstructed horizontal surface, in w / m 2 , , and is a constant on the 21st day of each month, obtained by fitting the measured data. When there is no terrain shielding, the direct radiation intensity of any slope ( I direct-0 ) can be used to calculate the normal direct radiation intensity ( I dn ) and the incident angle of direct solar radiation on the inclined surface ( α ), the specific formula is as follows: In the formula, The calculation formula is as follows: in, θ A is the azimuth of the study plane; θ S is the inclination of the study plane.

[0057] For the scattered radiation intensity of any surface without obstruction, since the scattered radiation is anisotropic and the propagation process is extremely complicated, in order to estimate the scattered radiation data of any surface, it is necessary to use a conversion model to decompose it into several isotropic parts and then perform conversion calculations.

[0058] Among them, the Perez model decomposes the solar diffuse radiation in space into three parts, namely, the sky dome diffuse radiation, the disc diffuse radiation and the horizontal diffuse radiation. The dome diffuse radiation is defined as the isotropically distributed diffuse radiation covering the entire dome. The disc diffuse radiation is considered to be the diffuse radiation around the sun with the sun as the center, and its intensity decreases with the distance from the center of the circle. The horizontal diffuse radiation is regarded as a diffuse radiation around the horizon, and its intensity is independent of the orientation.

[0059] According to the Perez model, the scattered radiation intensity (W / m 2 ) is calculated as follows: In the formula, The dome diffuse radiation intensity when there is no shielding, in W / m 2 ; The scattered radiation intensity of the optical disc when there is no shielding, in W / m 2 ; The horizontal scattered radiation intensity when there is no shielding, in W / m 2 ; The calculation formula is as follows: In the formula, To study the plane tilt angle; , It is a parameter describing the position relationship between the sun and the inclined surface; is the solar disk emissivity; is the horizontal radiation coefficient; is the horizontal scattered radiation, in W / m 2.

[0060] Since scattered radiation is anisotropic and unevenly distributed in actual space, it is difficult to quantify the shielding effect of terrain on scattered radiation. In this embodiment, based on the Perez model, the scattered radiation is decomposed into three isotropically distributed radiations, namely, solar disc scattered radiation, dome scattered radiation, and horizontal scattered radiation. Then, based on the principle of ray tracing, the degree of shielding of the three types of scattered radiation by terrain is quantified respectively, and finally the shielding effect of complex mountain terrain on scattered radiation of any surface is calculated.

[0061] Step S400: Calculate the second direct radiation intensity and the second scattered radiation intensity of any surface per minute under the influence of terrain shielding according to the visible field elevation angle, solar altitude angle and azimuth, and the inclination and orientation of the target research surface.

[0062] In this embodiment, based on the sun's motion trajectory and visible field, a ray tracing method is used to calculate the second direct radiation intensity and the second scattered radiation intensity of any surface minute by minute under the influence of terrain occlusion.

[0063] In order to quantify the effect of terrain on direct solar radiation and determine whether direct radiation will be blocked by the mountains around the study point, it is necessary to combine the trajectory of the sun's movement and compare the sun's altitude angle with the altitude angle of any visible field under terrain blocking.

[0064] If the solar altitude angle

[0065] Less than or equal to the altitude angle of any visible field of view under terrain occlusion

[0066] , then the direct sunlight will be blocked. For this, the direct radiation reception factor is introduced

[0067] To describe whether direct radiation is blocked.

[0068] It is defined as the ratio of the direct radiation intensity of the inclined surface under terrain shielding to the direct radiation intensity of the inclined surface without shielding.

[0069] Less than or equal to

[0070] When there is no direct radiation at the research point,

[0071] is 0, otherwise it is 1. The specific calculation formula is as follows: According to the classic model, the direct radiation intensity of the inclined surface without terrain shielding can be calculated ( ) is calculated, which is the normal direct radiation intensity when there is no shielding Cosine of the angle of incidence of direct solar radiation on the inclined surface Then, the direct radiation intensity of the inclined surface under the shielding of mountainous terrain (W / m2 ) can be calculated using the direct radiation intensity of the inclined surface without terrain shielding and the direct radiation receiving factor. The specific formula is as follows: For the scattered radiation intensity of any surface under terrain occlusion, in order to study the occlusion effect of terrain on scattered radiation, the Perez model is first used to convert the scattered radiation intensity of the inclined surface when it is not occluded into Decomposed into three parts, namely, the unobstructed isotropic dome diffuse radiation intensity , the scattered radiation intensity of the unobstructed solar disk , and the unobstructed horizontal diffuse radiation intensity , and then analyze the occlusion effect of the terrain on it.

[0072] The intensity of scattered radiation from the solar disk The occlusion judgment is the same as direct radiation, and its calculation is shown as follows: When calculating the diffuse radiation of the sky dome, the usual method is to divide the sky dome into many small blocks, calculate the radiation intensity of each small block separately, and then accumulate them. If the sky dome is divided more, a more accurate dome radiation intensity can be obtained. Therefore, in order to improve the simulation accuracy of the model, it is assumed that the number of small blocks is infinite, and each small block can be regarded as a point light source. Then the double integral method is used to calculate the diffuse radiation intensity of the sky dome under the occlusion of mountainous terrain. Dome diffuse reception factor It is defined as the ratio of the dome diffuse radiation intensity of any surface under terrain occlusion to the dome diffuse radiation intensity of any surface without terrain occlusion. After dividing the visible field into 360 sectors at intervals of 1° in azimuth, double integration is performed to obtain The calculation formula is as follows: in is the visible field elevation angle of the nth sector area; The elevation angle of the intersection line between the study plane and the sky hemisphere for the nth sector area.

[0073] In summary, the scattered radiation intensity of the inclined dome under terrain shading (W / m 2 ) is calculated as: Horizontal diffuse radiation reception factor ( ) is defined as the horizontal diffuse radiation intensity received by the inclined surface under terrain shading ( ) and the horizontal scattered radiation intensity of the inclined surface when there is no shielding ( ) ratio. Analysis of visible field It can be seen that the mountainous research site is surrounded by mountains, and the horizontal scattered radiation from around the horizon is usually completely blocked. Therefore, for mountainous areas with undulating terrain, this paper will It is set to a constant of 0, which means that the horizontal diffuse radiation received by the research point is completely blocked by the terrain.

[0074] Horizontal scattered radiation intensity of the inclined surface under the shielding of mountainous terrain (W / m 2 ) is calculated as: Since the position of the sun is constantly changing, the radiation received by the receiving surface is different at different times. The simulation model in this embodiment establishes a dynamic mathematical model for the position of the sun, the position of terrain obstructions, the position of the receiving surface, the orientation and inclination angle of the receiving surface, and can predict the exact position of the sun in the sky and the dynamic motion trajectory of the sun every minute, so as to analyze the radiation intensity and incident angle of the direct and scattered beams of the sun entering the receiving surface every minute, and can determine whether there are any obstructions on the incident path based on the principle of ray tracing, so as to calculate the solar radiation data under real obstruction conditions.

[0075] Step S500: By integrating the second direct radiation intensity and the second scattered radiation intensity over time, the direct solar radiation and scattered radiation of any surface under the influence of the mountainous terrain within any time interval are obtained, and then combined with the direct radiation and scattered radiation of any surface without considering the shielding effect, the direct radiation shielding rate and scattered radiation shielding rate of any surface under the influence of the mountainous terrain within any time interval are calculated.

[0076] Taking into account both the amount of calculation and the calculation accuracy, this embodiment sets 1 minute as the calculation step when calculating the solar radiation, calculates the solar radiation once every 1 minute, and assumes that the solar radiation at any time within the minute is a constant.

[0077] The definition of direct radiation intensity of any surface under mountain shading is the cumulative intensity of direct radiation intensity of any surface under mountain terrain shading within a fixed time period, in units of kW·h / m 2 The specific formula is: In the formula, For the i The direct solar radiation intensity under the mountain shade at the moment of one minute, in W / m 2 .

[0078] Similarly, the amount of solar scattered radiation on any surface under mountain shading (kW·h / m 2 ) is calculated as follows: In the formula, For the i The intensity of scattered solar radiation under mountain cover at the moment of one minute, in W / m 2 For the calculation of solar radiation shielding rate, in order to obtain the degree of shielding of horizontal solar radiation by mountain terrain in a specific time period, this embodiment introduces the solar radiation shielding rate in a specific time period. In a certain time period, the direct radiation shielding rate is defined as the amount of direct radiation when there is no terrain shielding. Direct radiation with terrain occlusion The ratio of the difference between the value of and the direct radiation when there is no terrain shielding. The calculation formula is as follows: Similarly, the calculation formula for the scattered radiation shielding rate in a certain period of time is as follows: in, is the amount of scattered radiation without shading, It is the amount of scattered radiation affected by terrain occlusion.

[0079] Through the shielding rate calculation formula, the direct radiation shielding rate and scattered radiation shielding rate can be calculated in any time period such as hourly, daily, monthly, seasonal or throughout the year.

[0080] Example 2

[0081] The embodiment of the present invention also provides a device for quantifying the effect of mountain terrain on solar radiation shielding, such as Figure 2 As shown, the device may include a data acquisition module 100 , an angle calculation module 200 , an unobstructed radiation intensity calculation module 300 , an obstructed radiation intensity calculation module 400 and a radiation obstruction rate calculation module 500 .

[0082] Specifically, the data acquisition module 100 is used to acquire the skyline and visible field elevation angle of the target research point under the influence of terrain occlusion based on the digital elevation model data.

[0083] The angle calculation module 200 is used to analyze the sun's motion trajectory according to the latitude, longitude and altitude of the target research point, and calculate the sun's altitude angle and azimuth angle based on the full-year time dimension with minutes as the calculation step.

[0084] The unobstructed radiation intensity calculation module 300 is used to calculate the first direct radiation intensity and the first scattered radiation intensity of any surface every minute without considering the influence of obstruction through the solar altitude angle and azimuth angle as well as the inclination and orientation of the target research surface.

[0085] The obstructed radiation intensity calculation module 400 is used to calculate the second direct radiation intensity and the second scattered radiation intensity of any surface minute by minute under the influence of terrain obstruction according to the visible field elevation angle, solar altitude angle and azimuth, and the inclination and orientation of the target research surface.

[0086] The radiation shielding rate calculation module 500 is used to obtain the solar direct radiation and scattered radiation of any surface under the influence of mountainous terrain within any time interval by integrating the second direct radiation intensity and the second scattered radiation intensity over time, and then calculate the direct radiation shielding rate and scattered radiation shielding rate of any surface under the influence of mountainous terrain within any time interval by combining the direct radiation and scattered radiation of any surface without considering the shielding effect.

[0087] In one implementation, the data acquisition module 100 is specifically used to:

[0088] Converting the digital elevation model data into an irregular triangular grid, and determining the specific location of the target research point in the irregular triangular grid according to the longitude and latitude of the target research point;

[0089] Obtaining the skyline of the target research point under the influence of terrain occlusion;

[0090] Determine the visible field of view of the horizontal plane research point under the influence of terrain occlusion and obtain the visible field elevation angle on the visible field boundary corresponding to each degree of azimuth.

[0091] In one implementation, the unobstructed radiation intensity calculation module 300 is specifically used to:

[0092] The ASHRAE and Perez model algorithms are used to calculate the first direct radiation intensity and the first scattered radiation intensity of any surface every minute without shielding.

[0093] In one implementation, the shielded radiation intensity calculation module 400 is specifically used to:

[0094] Based on the sun's motion trajectory and ray tracing method, the visible field is used to calculate the second direct radiation intensity and the second scattered radiation intensity of any surface minute by minute under the influence of terrain occlusion.

[0095] In summary, the embodiment of the present invention first obtains the skyline and visible field elevation angle of the target research point under the influence of terrain occlusion, analyzes the movement trajectory of the sun position according to the longitude and latitude and the altitude, and calculates the solar altitude angle and azimuth; then calculates the first direct radiation intensity and the first scattered radiation intensity of any surface every minute when the occlusion effect is not considered, and calculates the second direct radiation intensity and the second scattered radiation intensity of any surface every minute when the terrain occlusion effect is considered; then performs an integral operation to obtain the solar direct radiation and scattered radiation of any surface under the influence of mountainous terrain in any time interval; finally, the direct radiation and scattered radiation of any surface when the occlusion effect is not considered are combined to calculate the direct radiation occlusion rate and scattered radiation occlusion rate of any surface under the influence of mountainous terrain in any time interval, thereby quantifying the solar radiation data under the occlusion of complex terrain in mountainous areas. The embodiment of the present invention can be used for the correction calculation of solar radiation data of any surface under the influence of mountainous terrain, and effectively improves the accuracy of the prediction of solar radiation occlusion data of any receiving surface by mountainous terrain.

[0096] It can be understood that the same or similar parts of the above embodiments can be referenced to each other, and the contents not described in detail in some embodiments can refer to the same or similar contents in other embodiments.

[0097] It should be noted that, in the description of this application, the terms "first", "second", etc. are only used for descriptive purposes and cannot be understood as indicating or implying relative importance. In addition, in the description of this application, unless otherwise specified, the meaning of "plurality" refers to at least two.

[0098] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, segment or portion of code that includes one or more executable instructions for implementing the steps of a specific logical function or process, and the scope of the preferred embodiments of the present application includes alternative implementations in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in the reverse order depending on the functions involved, which should be understood by technicians in the technical field to which the embodiments of the present application belong.

[0099] It should be understood that the various parts of the present application can be implemented with hardware, software, firmware or a combination thereof. In the above-mentioned embodiments, multiple steps or methods can be implemented with software or firmware stored in a memory and executed by a suitable instruction execution system. A person of ordinary skill in the art can understand that all or part of the steps carried by the above-mentioned embodiment method can be completed by instructing the relevant hardware through a program, and the program can be stored in a computer-readable storage medium, which, when executed, includes one of the steps of the method embodiment or a combination thereof.

[0100] In addition, each functional unit in each embodiment of the present application may be integrated into a processing module, or each unit may exist physically separately, or two or more units may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.

[0101] The storage medium mentioned above can be a read-only memory, a magnetic disk or an optical disk, etc.

[0102] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representation of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner.

[0103] Although the embodiments of the present application have been shown and described above, it can be understood that the above embodiments are exemplary and cannot be understood as limitations on the present application. Ordinary technicians in this field can change, modify, replace and modify the above embodiments within the scope of the present application.

Claims

1. A method for quantifying the effect of mountain terrain on solar radiation shielding, characterized in that: The following steps are involved: Based on the digital elevation model data, the skyline and visible field elevation angle of the target research point under the influence of terrain occlusion are obtained; Analyze the sun's motion trajectory according to the latitude, longitude and altitude of the target research point, and calculate the sun's altitude angle and azimuth angle based on the full-year time dimension with a calculation step of minutes; By using the solar altitude angle and azimuth as well as the inclination and orientation of the target research surface, the first direct radiation intensity and the first scattered radiation intensity of any surface are calculated every minute without considering the influence of shielding; Calculate the second direct radiation intensity and the second scattered radiation intensity of any surface per minute under the influence of terrain shielding according to the visible field elevation angle, the solar altitude angle and the azimuth angle, and the inclination and orientation of the target research surface; By integrating the second direct radiation intensity and the second scattered radiation intensity over time, the direct solar radiation and scattered radiation of any surface under the influence of the mountainous terrain in any time interval are obtained, and then the direct radiation and scattered radiation of any surface without considering the shielding effect are combined to calculate the direct radiation shielding rate and scattered radiation shielding rate of any surface under the influence of the mountainous terrain in any time interval; wherein, The direct radiation shielding rate of any surface under the influence of mountain terrain in any time interval = (direct radiation without terrain shielding - direct radiation with terrain shielding) / direct radiation without terrain shielding; The scattered radiation shielding rate of any surface under the influence of mountain terrain in any time interval = (scattered radiation without shielding - scattered radiation under the influence of terrain shielding) / scattered radiation without shielding; When there is terrain occlusion The unit is W / m 2 , I direct-i is the direct solar radiation intensity under the mountain shade at the i-th minute; Under the influence of terrain occlusion The unit is W / m 2 , E diffuse-i is the intensity of solar scattered radiation under mountain shade at i minute; The step of obtaining the skyline and visible field elevation angle of the target research point under the influence of terrain occlusion based on the digital elevation model data specifically includes: Converting the digital elevation model data into an irregular triangular grid, and determining the specific location of the target research point in the irregular triangular grid according to the longitude and latitude of the target research point; Obtaining the skyline of the target research point under the influence of terrain occlusion; Determine the visible field of view of the horizontal plane research point under the influence of terrain occlusion and obtain the visible field elevation angle at the visible field boundary corresponding to each degree of azimuth; The step of calculating the first direct radiation intensity and the first scattered radiation intensity of any surface every minute without considering the influence of occlusion specifically includes: The ASHRAE and Perez model algorithms are used to calculate the first direct radiation intensity and the first scattered radiation intensity of any surface every minute without shielding; The step of calculating the second direct radiation intensity and the second scattered radiation intensity of any surface minute by minute under the influence of terrain occlusion specifically includes: Based on the sun's motion trajectory and ray tracing method, the visible field is used to calculate the second direct radiation intensity and the second scattered radiation intensity of any surface minute by minute under the influence of terrain occlusion.

2. A device for quantifying the effect of mountain terrain on solar radiation shielding, characterized in that: The device comprises: A data acquisition module is used to obtain the skyline and visible field elevation angle of the target research point under the influence of terrain occlusion based on digital elevation model data; An angle calculation module, used to analyze the sun's motion trajectory according to the latitude and longitude and altitude of the target research point, and calculate the sun's altitude angle and azimuth angle based on the full-year time dimension and in minutes as the calculation step; An unobstructed radiation intensity calculation module is used to calculate the first direct radiation intensity and the first scattered radiation intensity of any surface per minute without considering the influence of obstruction by using the solar altitude angle and azimuth as well as the inclination and orientation of the target research surface; There is an obstruction radiation intensity calculation module, which is used to calculate the second direct radiation intensity and the second scattered radiation intensity of any surface minute by minute under the influence of terrain obstruction according to the visible field elevation angle, solar altitude angle and azimuth, and the inclination and orientation of the target research surface; The radiation shielding rate calculation module is used to obtain the solar direct radiation and scattered radiation of any surface under the influence of mountainous terrain in any time interval by integrating the second direct radiation intensity and the second scattered radiation intensity over time, and then calculate the direct radiation shielding rate and scattered radiation shielding rate of any surface under the influence of mountainous terrain in any time interval by combining the direct radiation and scattered radiation of any surface without considering the shielding effect; wherein, The direct radiation shielding rate of any surface under the influence of mountain terrain in any time interval = (direct radiation without terrain shielding - direct radiation with terrain shielding) / direct radiation without terrain shielding; The scattered radiation shielding rate of any surface under the influence of mountain terrain in any time interval = (scattered radiation without shielding - scattered radiation under the influence of terrain shielding) / scattered radiation without shielding; Straightness when there is terrain occlusion The unit is W / m 2 , I direct-i is the direct solar radiation intensity under the mountain shade at the i-th minute; Under the influence of terrain occlusion The unit is W / m 2 , E diffuse-i is the intensity of solar scattered radiation under mountain shade at i minute; The data acquisition module is specifically used for: Converting the digital elevation model data into an irregular triangular grid, and determining the specific location of the target research point in the irregular triangular grid according to the longitude and latitude of the target research point; Obtaining the skyline of the target research point under the influence of terrain occlusion; Determine the visible field of view of the horizontal plane research point under the influence of terrain occlusion and obtain the visible field elevation angle at the visible field boundary corresponding to each degree of azimuth; The unobstructed radiation intensity calculation module is specifically used for: The ASHRAE and Perez model algorithms are used to calculate the first direct radiation intensity and the first scattered radiation intensity of any surface every minute without shielding; The shielded radiation intensity calculation module is specifically used for: Based on the sun's motion trajectory and ray tracing method, the visible field is used to calculate the second direct radiation intensity and the second scattered radiation intensity of any surface minute by minute under the influence of terrain occlusion.

Images