Method and device for calculating total radiation quantity of horizontal plane and storage medium

Abstract

The invention provides a method and device for calculating the total radiation energy of a horizontal plane and a storage medium, and the method comprises the following steps: (1), obtaining a daily clearance index in a year through a random number algorithm according to an average clearness index in each month in the year, constructing a TAG model, and calculating a clearness index in each hour in the year; (2) calculating the extraterrestrial horizontal clear sky radiation quantity every day in one year according to the latitude position, and then calculating the clear sky radiation quantity of the extraterrestrial horizontal plane every hour in one year; and (3) multiplying the clearness index per hour in one year obtained in the step (1) by the extraterrestrial horizontal clearness radiation quantity per hour in one year obtained in the step (2) to obtain the total horizontal plane radiation quantity per hour in one year. According to the method, the irradiation data per day and per hour are obtained innovatively through the algorithm according to the monthly average irradiation data and the latitude position, and the problems that the prior art is only suitable for a small part of positions of each country or region, equipment is expensive, and time is delayed are solved.

Description

Technical Field

[0001] The invention discloses a method, a device and a storage medium for calculating the total radiation energy of a horizontal plane, and belongs to the field of illumination calculation. Background Art

[0002] Since many design codes and studies require hourly irradiance data, monthly average irradiance data are no longer sufficient for today's needs. However, since the measured radiation value series are only available for a small number of locations in each country or region, there are still problems in many cases, such as when irradiance data are required for inclined surfaces instead of horizontal surfaces and the recording time is not long enough or there are gaps in the recording process. Therefore, in order to generate hourly values ​​at any desired location, random models are used. The values ​​generated by random models have the same statistical properties as the measured data, namely the mean, variance and characteristic sequence (autocorrelation). And the generated data approximate the natural characteristics to a large extent. Recent studies have also shown that the data generated in this way can be a good substitute for long-term measured data. Generating synthetic radiation values ​​is often the only practical way to obtain radiation data on a daily or hourly time scale.

[0003] The radiation value series measured by traditional methods are only applicable to a small number of locations in each country or region, and there are problems in many cases, such as when the recording time is not long enough or there are intervals in the recording process. In addition, some high-precision measurement equipment may be relatively expensive and require regular maintenance and calibration to ensure the accuracy of the data. The terrain, obstructions and meteorological conditions at the same location may cause data deviations. There may be time delays in obtaining some data, which is not suitable for applications that require real-time data. All of these may lead to low accuracy of the results. Summary of the invention

[0004] Aiming at the fact that the radiation value sequence measured by traditional methods is only applicable to a small part of each country or region and the measurement value has low accuracy, the present invention proposes a random number algorithm and a TAG model, which can more conveniently calculate the daily clear sky index and the hourly clear sky index and radiation.

[0005] The technical solution of the present invention is as follows:

[0006] A method for calculating the total radiation energy of a horizontal surface comprises the following steps:

[0007] (1) The clear sky index for each day of the year is obtained by using a random number algorithm based on the average clear sky index for each month of the year. Then a TAG model is constructed to calculate the clear sky index for each hour of the year.

[0008] (2) Calculate the clear sky radiation outside the Earth's horizontal plane every day of the year from the latitude position, and then calculate the clear sky radiation outside the Earth's horizontal plane every hour of the year;

[0009] (3) Multiply the clear sky index for each hour of the year obtained in step (1) by the extraterrestrial horizontal clear sky radiation for each hour of the year obtained in step (2) to obtain the total horizontal radiation for each hour of the year.

[0010] Preferably, the random number algorithm of step (1) above adopts a Markov transfer matrix table, and is calculated using multiple sites, where the site areas cover the main climate zones. The result is a 10*10*10 matrix. The specific steps are as follows:

[0011] (2-1) Determine the matrix in the Markov transfer matrix table according to the interval of the average clear sky index of this month, as well as the maximum and minimum values ​​of the monthly average clear sky index; and divide the maximum and minimum values ​​into 10 intervals, the interval less than the minimum value is the first interval, and the interval greater than the maximum value is the 12th interval;

[0012] (2-2) The average clear sky index on the first day of the month is defined as the average clear sky index of the previous month;

[0013] (2-3) The i-th row of the determined matrix is ​​determined by the average clear sky index of the previous day in the i-th interval determined above, and if i>10, i=10;

[0014] (2-4) Generate a random number R between 0 and 1. The average clear sky index for the next day is the sum of the values ​​in the i-th row until it is greater than R.

[0015] (2-5) Repeat steps (2-3) and (2-4) to obtain the average clear sky index for each day of a month.

[0016] Preferably, the specific steps of calculating the clear sky index per hour in a year in the above step (1) are as follows:

[0017] (3-1) Calculate the sun angle and sunrise and sunset angles at the beginning, middle and end of each hour:

[0018]

[0019] ω=arccos[-tanδ*tan(lat)]

[0020] Where ω1 represents the sun angle at the beginning of each hour, ω2 represents the sun angle at the middle of each hour, ω3 represents the sun angle at the end of each hour, ω represents the sunrise and sunset angles, lat is the latitude in radians, and δ is the solar declination angle calculated as follows:

[0021]

[0022] Where n is the number of days in a year;

[0023] If ω1≤-ω or ω2≥ω, it is non-sunshine time, and the clear sky index for each hour is 0;

[0024] (3-2) If ω1>-ω and ω2<ω, the model consists of two parts: the first part calculates the daily average clear sky index; the second part simulates the intermittent hourly changes by superimposing a first-order autoregressive process;

[0025] k t =k tM +σy(h)

[0026] where k t is the hourly clear sky index, k tm is the average clear sky index; σ is the standard deviation of the clear sky index for each hour of the day corresponding to the clear sky index for each day.

[0027]

[0028] The algorithm constants λ, ε, k and the solar altitude angle h are s They are:

[0029]

[0030] ε=0.32-1.6(K t -0.5)

[0031]

[0032] h s =arcsin(cosω2*cosδ*coslat+sinδ*sinlat)

[0033] According to the study of the sequential characteristics of the y variable, its current hourly value is only significantly dependent on the value of the previous hour, so the autoregressive function y(h) is defined as:

[0034]

[0035] In the formula, is the autocorrelation coefficient: r is a randomly generated variable that conforms to the Gaussian distribution, and its expected value is zero:

[0036]

[0037] In the formula, z is a floating point number greater than or equal to 0 and less than 1, following a uniform distribution, that is, it randomly falls in the half-open interval [0.0, 1.0), σ′ where the standard deviation σ is: σWhere A and B are algorithm constants:

[0038] When calculating, if the hourly clear sky index k t <0 or k t >k cs , where k cs The clear sky index for clear skies is:

[0039]

[0040] Then regenerate a floating point number z, calculate r, y(h) and k t , and set the maximum number of iterations to 10. If the iterations exceed 10, if k t <0, then k t =0, if k t >k cs , then k t =k cs ;

[0041] In this way, a clear sky index for each hour of the year can be generated.

[0042] Preferably, the specific calculation steps of the above step (2) are as follows:

[0043] (4-1) Calculate the daily clear sky radiation B0d on the outer horizontal surface:

[0044]

[0045] Where ω is the sunrise and sunset angle, δ is the solar declination angle, lat is the latitude in radians, and ε is the correction factor for the earth's eccentricity:

[0046]

[0047] Where n is the day of the year;

[0048] (4-2) Calculate the hourly clear sky radiation G0c on the outer horizontal surface:

[0049] If ω1≤-ω or ω2≥ω, it is non-sunshine time, and its value is 0.

[0050] If ω1>-ω and ω2<ω,

[0051]

[0052] Where ω2 is the sun angle at the middle of each hour, and ω is the sunrise and sunset angle.

[0053] The present invention further provides a calculation device for calculating the total radiation energy of a horizontal surface, comprising a memory and a processor, wherein the memory is used to store a computer program, and the processor is used to read the computer program and execute the above-mentioned method for calculating the total radiation energy of a horizontal surface.

[0054] The present invention also provides a computer-readable storage medium, on which instructions are stored. When the instructions are run on a computer, the above-mentioned method for calculating the total radiation energy of a horizontal surface is executed.

[0055] Key Terms:

[0056] Clear sky index: refers to the ratio of the total solar radiation incident on the horizontal plane to the astronomical radiation, with a value between 0 and 1.

[0057] Radiation: It is the radiation energy per unit area per unit time on the solid earth surface after the direct, scattered and reflected effects of the atmosphere. Its unit is: Watt / square meter (W / ㎡).

[0058] Markov transfer matrix: It is a random process with the "Markov property", that is, given the current state, the probability distribution of the future state depends only on the current state and not on the past state.

[0059] TAG model: It is a method for time series modeling that combines the concepts of autoregressive model and Gaussian process. It is mainly used to process time series data, which contains time correlation and uncertainty modeling.

[0060] Clear sky radiation on the outer horizontal plane: refers to the radiation energy of solar radiation directly hitting the horizontal plane of the earth outside the earth's atmosphere (outside the earth). This value is usually used to describe the maximum intensity of solar radiation on the earth under ideal conditions, that is, without the influence of atmospheric absorption and scattering. Its unit is: Watt / square meter (W / ㎡).

[0061] The beneficial effects of the present invention are:

[0062] The present invention innovatively obtains daily and hourly irradiation data through an algorithm based on the monthly average irradiation data and the latitude position, overcoming the problems that existing technologies are only applicable to a small number of locations in each country or region, expensive equipment, and time delays. BRIEF DESCRIPTION OF THE DRAWINGS

[0063] Figure 1 It is a flow chart of the method of the present invention. DETAILED DESCRIPTION

[0064] The present invention will be further described below in conjunction with the accompanying drawings. The following embodiments are only used to more clearly illustrate the technical solution of the present invention, and cannot be used to limit the protection scope of the present invention.

[0065] like Figure 1 As shown, a method for calculating the total radiation energy of a horizontal surface comprises the following steps:

[0066] (1) The clear sky index for each day of the year is obtained by using a random number algorithm based on the average clear sky index for each month of the year. Then a TAG model is constructed to calculate the clear sky index for each hour of the year.

[0067] (2) Calculate the clear sky radiation outside the Earth's horizontal plane every day of the year from the latitude position, and then calculate the clear sky radiation outside the Earth's horizontal plane every hour of the year;

[0068] (3) Multiply the clear sky index for each hour of the year obtained in step (1) by the extraterrestrial horizontal clear sky radiation for each hour of the year obtained in step (2) to obtain the total horizontal radiation for each hour of the year.

[0069] Preferably, the random number algorithm of step (1) above adopts a Markov transfer matrix table and is calculated using multiple sites, where the site areas cover the main climate zones, resulting in a 10*10*10 matrix.

[0070] Firstly, a table of known Markov transfer matrices is given, which is a matrix of 10 rows and 10 columns, where Kt is the monthly average clear sky index; and the maximum and minimum values ​​of the clear sky index in the 10 matrices defined by the monthly average clear sky index.

[0071] Markov transition matrix table and maximum and minimum value table:

[0072] Matrix 1: Kt≤0.30

[0073] [[0.229,0.333,0.208,0.042,0.083,0.042,0.042,0.021,0.000,0.000],

[0074] [0.167,0.319,0.194,0.139,0.097,0.028,0.042,0.000,0.014,0.000],

[0075] [0.250,0.250,0.091,0.136,0.091,0.046,0.046,0.023,0.068,0.000],[0.158,0.237,0.158,0.263,0.026,0.053,0.079,0.026,0.000,0.000],[0.2 11,0.053,0.211,0.158,0.053,0.053,0.158,0.105,0.000,0.000],[0.125,0.125,0.250,0.188,0.063,0.125,0.000,0.125,0.000,0.000],[0.040,0 .240,0.080,0.120,0.080,0.080,0.120,0.120,0.080,0.040],[0.000,0.250,0.000,0.125,0.000,0.125,0.125,0.250,0.063,0.063],[0.000,0.250 ,0.000,0.125,0.250,0.000,0.250,0.000,0.000,0.125],[0.000,0.000,0.000,0.000,0.000,0.000,0.500,0.250,0.000,0.250]])Matrix 2: 0.30<Kt≤0.35

[0076] [[0.000,0.000,0.091,0.000,0.364,0.091,0.182,0.000,0.273,0.000],[0.118,0.118,0.176,0.118,0.059,0.118,0.176,0.059,0.059,0.000],[0.067,0.267,0.067,0.200,0.067,0.0 00,0.133,0.133,0.000,0.067],[0.118,0.235,0.000,0.235,0.059,0.176,0.118,0.000,0.059,0.000],[0.077,0.154,0.308,0.077,0.154,0.077,0.000,0.077,0.000],[0.083,0 .000,0.167,0.250,0.083,0.167,0.000,0.083,0.167,0.000],[0.222,0.222,0.000,0.111,0.111,0.000,0.111,0.222,0.000,0.000],[0.091,0.182,0.273,0.000,0.091,0.273,0.000, Matrix 3: 0.35<Kt≤0.40

[0077] [[0.206,0.088,0.176,0.176,0.088,0.029,0.176,0.029,0.029,0.000],[0.120,0.100,0.140,0.160,0.120,0.220,0.100,0.000,0.020,0.020],[0.077,0.123,0.185,0.123,0.077,0.1 39,0.092,0.123,0.061,0.000],[0.048,0.111,0.095,0.206,0.206,0.190,0.095,0.048,0.000,0.000],[0.059,0.137,0.118,0.137,0.098,0.118,0.118,0.157,0.059,0.000],[0.014,0 .097,0.139,0.153,0.125,0.139,0.208,0.056,0.042,0.028],[0.073,0.101,0.116,0.145,0.087,0.159,0.203,0.087,0.029,0.000],[0.019,0.037,0.111,0.056,0.074,0.111,0.185, Matrix 4: 0.40<Kt≤0.45

[0078] [[0.167,0.167,0.167,0.000,0.083,0.125,0.000,0.167,0.125,0.000],[0.117,0.117,0.150,0.117,0.083,0.117,0.200,0.067,0.017,0.017],[0.049,0.085,0.134,0.158,0.098,0.1 10,0.134,0.134,0.061,0.037],[0.039,0.090,0.141,0.141,0.167,0.141,0.090,0.141,0.039,0.013],[0.009,0.139,0.074,0.093,0.194,0.139,0.167,0.093,0.074,0.019],[0.036,0 .018,0.117,0.099,0.144,0.180,0.180,0.117,0.072,0.036],[0.000,0.046,0.061,0.061,0.136,0.159,0.273,0.167,0.098,0.000],[0.016,0.056,0.080,0.128,0.104,0.080,0.160, Matrix 5: 0.45<Kt≤0.50

[0079] [[0.120,0.200,0.160,0.120,0.120,0.120,0.080,0.000,0.040,0.040],[0.100,0.080,0.120,0.140,0.140,0.200,0.180,0.040,0.000,0.000],[0.046,0.114,0.068,0.171,0.125,0.1 71,0.080,0.159,0.057,0.011],[0.015,0.061,0.084,0.099,0.191,0.153,0.153,0.115,0.115,0.015],[0.024,0.030,0.098,0.098,0.165,0.195,0.195,0.140,0.043,0.012],[0.015,0 .026,0.062,0.124,0.144,0.170,0.170,0.222,0.062,0.005],[0.000,0.013,0.045,0.108,0.112,0.175,0.188,0.224,0.117,0.018],[0.008,0.023,0.054,0.066,0.093,0.125,0.191, Matrix 6: 0.50<Kt≤0.55

[0080] [[0.250,0.179,0.107,0.107,0.143,0.071,0.107,0.036,0.000,0.000],[0.133,0.022,0.089,0.111,0.156,0.178,0.111,0.133,0.067,0.000],[0.064,0.048,0.143,0.048,0.175,0.1 43,0.206,0.095,0.079,0.000],[0.000,0.022,0.078,0.111,0.156,0.156,0.244,0.167,0.044,0.022],[0.016,0.027,0.037,0.069,0.160,0.219,0.230,0.160,0.075,0.005],[0.013,0 .025,0.030,0.093,0.144,0.202,0.215,0.219,0.055,0.004],[0.006,0.041,0.035,0.064,0.090,0.180,0.337,0.192,0.049,0.006],[0.012,0.021,0.029,0.035,0.132,0.123,0.184, 0.371,0.082,0.012],[0.008,0.016,0.016,0.024,0.071,0.103,0.159,0.270,0.309,0.024],[0.000,0.000,0.000,0.000,0.059,0.000,0.059,0.294,0.412,0.176]])Matrix 7: 0.55<Kt≤0.60

[0081] [[0.217,0.087,0.000,0.174,0.130,0.087,0.087,0.130,0.087,0.000],[0.026,0.079,0.132,0.079,0.026,0.158,0.158,0.132,0.158,0.053],[0.020,0.020,0.020,0.040,0.160,0 .180,0.160,0.200,0.100,0.100],[0.025,0.013,0.038,0.076,0.076,0.139,0.139,0.266,0.215,0.013],[0.030,0.030,0.050,0.020,0.091,0.131,0.162,0.283,0.131,0.071],[0.0 06,0.006,0.013,0.057,0.057,0.121,0.204,0.287,0.185,0.064],[0.004,0.026,0.037,0.030,0.093,0.107,0.193,0.307,0.167,0.037],[0.011,0.009,0.014,0.042,0.041,0.071, Matrix 8: 0.60 <Kt≤0.65

[0082] [[0.067,0.133,0.133,0.067,0.067,0.200,0.133,0.133,0.067,0.000],[0.118,0.059,0.059,0.059,0.059,0.118,0.118,0.235,0.118,0.059],[0.000,0.024,0.024,0.049,0.146,0.0 73,0.195,0.244,0.195,0.049],[0.026,0.000,0.026,0.026,0.053,0.184,0.263,0.184,0.237,0.000],[0.014,0.000,0.042,0.056,0.069,0.097,0.139,0.306,0.278,0.000],[0.009,0 .009,0.052,0.069,0.052,0.112,0.215,0.285,0.138,0.060],[0.009,0.009,0.026,0.017,0.094,0.099,0.232,0.283,0.210,0.021],[0.010,0.014,0.016,0.019,0.027,0.062,0.163, 0.467,0.202,0.019],[0.004,0.007,0.031,0.017,0.033,0.050,0.086,0.252,0.469,0.050],[0.000,0.000,0.015,0.046,0.031,0.046,0.077,0.123,0.446,0.215]])Matrix 9: 0.65<Kt≤0.70

[0083] [[0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,1.000,0.000],[0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,1.000,0.000],[0.000,0.000,0.000,0.000,0.000,0.000,0. 000,0.250,0.250,0.500,0.000],[0.000,0.000,0.000,0.000,0.250,0.000,0.000,0.375,0.250,0.125],[0.000,0.000,0.000,0.083,0.000,0.167,0.167,0.250,0.333,0.000],[0.000 ,0.000,0.042,0.042,0.042,0.083,0.083,0.292,0.292,0.125],[0.000,0.000,0.032,0.000,0.000,0.032,0.129,0.387,0.355,0.065],[0.000,0.000,0.000,0.038,0.038,0.075,0.0 47,0.340,0.415,0.047],[0.004,0.004,0.007,0.007,0.011,0.030,0.052,0.141,0.654,0.089],[0.000,0.000,0.000,0.000,0.061,0.061,0.030,0.030,0.349,0.470]])Matrix 10: Kt>0.70

[0084] [[0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,1.000,0.000],[0.100,0.100,0.100,0.100,0.100,0.100,0.100,0.100,0.100,0.100,0.100],[0.000,0.000,0.000,0.250,0.000 ,0.000,0.000,0.500,0.250,0.000],[0.000,0.000,0.143,0.143,0.000,0.143,0.429,0.000,0.000],[0.000,0.000,0.000,0.200,0.000,0.000,0.200,0.400,0.200,0.000], [0.000,0.000,0.000,0.000,0.000,0.000,0.222,0.444,0.333,0.000],[0.000,0.000,0.000,0.000,0.080,0.080,0.480,0.240,0.040],[0.000,0.000,0.027,0.009,0.027, 0.018,0.135,0.523,0.252,0.009],[0.000,0.000,0.000,0.022,0.000,0.043,0.043,0.326,0.511,0.054],[0.000,0.000,0.000,0.143,0.000,0.000,0.000,0.143,0.714,0.000]])

[0085]

[0086] (2-1) Determine the matrix in the Markov transfer matrix table according to the interval of the average clear sky index of this month, as well as the maximum and minimum values ​​of the monthly average clear sky index; and divide the maximum and minimum values ​​into 10 intervals, the interval less than the minimum value is the first interval, and the interval greater than the maximum value is the 12th interval;

[0087] (2-2) The average clear sky index on the first day of the month is defined as the average clear sky index of the previous month;

[0088] (2-3) The i-th row of the determined matrix is ​​determined by the average clear sky index of the previous day in the i-th interval determined above, and if i>10, i=10;

[0089] (2-4) Generate a random number R between 0 and 1. The average clear sky index for the next day is the sum of the values ​​in the i-th row until it is greater than R.

[0090] (2-5) Repeat steps (2-3) and (2-4) to obtain the average clear sky index for each day of a month.

[0091] Preferably, the specific steps of calculating the clear sky index per hour in a year in the above step (1) are as follows:

[0092] (3-1) Calculate the sun angle and sunrise and sunset angles at the beginning, middle and end of each hour:

[0093]

[0094] ω=arccos[-tanδ*tan(lat)]

[0095] Where ω1 represents the sun angle at the beginning of each hour, ω2 represents the sun angle at the middle of each hour, ω3 represents the sun angle at the end of each hour, ω represents the sunrise and sunset angles, lat is the latitude in radians, and δ is the solar declination angle calculated as follows:

[0096]

[0097] Where n is the number of days in a year;

[0098] If ω1≤-ω or ω2≥ω, it is non-sunshine time, and the clear sky index for each hour is 0;

[0099] (3-2) If ω1>-ω and ω2<ω, the model consists of two parts: the first part calculates the daily average clear sky index; the second part simulates the intermittent hourly changes by superimposing a first-order autoregressive process;

[0100] k t =k tM +σy(h)

[0101] where k t is the hourly clear sky index, k tm is the average clear sky index; σ is the standard deviation of the clear sky index for each hour of the day corresponding to the clear sky index for each day.

[0102]

[0103] Wherein the algorithm constants λ, ∈, k and the solar altitude angle h s They are:

[0104]

[0105] ∈=0.32-1.6(K t -0.5)

[0106]

[0107] h s =arcsin(cosω2*cosδ*coslat+sinδ*sinlat)

[0108] According to the study of the sequential characteristics of the y variable, its current hourly value is only significantly dependent on the value of the previous hour, so the autoregressive function y(h) is defined as:

[0109]

[0110] In the formula, is the autocorrelation coefficient: r is a randomly generated variable that conforms to the Gaussian distribution, and its expected value is zero:

[0111]

[0112] In the formula, z is a floating point number greater than or equal to 0 and less than 1, following a uniform distribution, that is, it randomly falls in the half-open interval [0.0, 1.0), σ′ where the standard deviation σ is: σWhere A and B are algorithm constants:

[0113] When calculating, if the hourly clear sky index k t <0 or k t >k cs , where k cs The clear sky index for clear skies is:

[0114]

[0115] Then regenerate a floating point number z, calculate r, y(h) and k t , and set the maximum number of iterations to 10. If the iterations exceed 10, if k t <0, then k t =0, if k t >k cs , then k t =k cs ;

[0116] In this way, a clear sky index for each hour of the year can be generated.

[0117] Preferably, the specific calculation steps of the above step (2) are as follows:

[0118] (4-1) Calculate the daily clear sky radiation B0d on the outer horizontal surface:

[0119]

[0120] Where ω is the sunrise and sunset angle, δ is the solar declination angle, lat is the latitude in radians, and ε is the correction factor for the earth's eccentricity:

[0121]

[0122] Where n is the day of the year;

[0123] (4-2) Calculate the hourly clear sky radiation G0c on the outer horizontal surface:

[0124] If ω1≤-ω or ω2≥ω, it is non-sunshine time, and its value is 0.

[0125] If ω1>-ω and ω2<ω,

[0126]

[0127] Where ω2 is the sun angle at the middle of each hour, and ω is the sunrise and sunset angle.

[0128] Calculate the hourly radiation G0 on the horizontal surface during the year:

[0129] G0=G0c*k t

[0130] Where G0c is the hourly clear sky radiation on the outer horizontal surface, k t It is the clear sky index for each hour of the year.

[0131] The above is only a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the technical principles of the present invention. These improvements and modifications should also be regarded as the scope of protection of the present invention.

Claims

1. A method for calculating the total radiation energy of a horizontal surface, characterized in that: The steps include: (1) The clear sky index for each day of the year is obtained by using a random number algorithm based on the average clear sky index for each month of the year. Then a TAG model is constructed to calculate the clear sky index for each hour of the year. (2) Calculate the clear sky radiation outside the Earth's horizontal plane every day of the year from the latitude position, and then calculate the clear sky radiation outside the Earth's horizontal plane every hour of the year; (3) Multiply the clear sky index for each hour of the year obtained in step (1) by the extraterrestrial horizontal clear sky radiation for each hour of the year obtained in step (2) to obtain the total horizontal radiation for each hour of the year.

2. A method for calculating the total radiation energy of a horizontal surface according to claim 1, characterized in that: The random number algorithm of step (1) adopts a Markov transfer matrix table and is calculated using multiple sites, where the site areas cover the main climate zones. The result is a 10*10*10 matrix. The specific steps are as follows: (2-1) Determine the matrix in the Markov transfer matrix table according to the interval of the average clear sky index of this month, as well as the maximum and minimum values ​​of the monthly average clear sky index; and divide the maximum and minimum values ​​into 10 intervals, the interval less than the minimum value is the first interval, and the interval greater than the maximum value is the 12th interval; (2-2) The average clear sky index on the first day of the month is defined as the average clear sky index of the previous month; (2-3) The i-th row of the determined matrix is ​​determined by the average clear sky index of the previous day in the i-th interval determined above, and if i>10, i=10; (2-4) Generate a random number R between 0 and 1. The average clear sky index for the next day is the sum of the values ​​in the i-th row until it is greater than R. (2-5) Repeat steps (2-3) and (2-4) to obtain the average clear sky index for each day of a month.

3. The method for calculating the total radiation energy of a horizontal plane according to claim 1, characterized in that: The specific steps of calculating the clear sky index per hour in a year in step (1) are as follows: (3-1) Calculate the sun angle and sunrise and sunset angles at the beginning, middle and end of each hour: ω=arccos[-tanδ*tan(lat)] Where ω1 represents the sun angle at the beginning of each hour, ω2 represents the sun angle at the middle of each hour, ω3 represents the sun angle at the end of each hour, ω represents the sunrise and sunset angles, lat is the latitude in radians, and δ is the solar declination angle calculated as follows: Where n is the number of days in a year; If ω1≤-ω or ω2≥ω, it is non-sunshine time, and the clear sky index for each hour is 0; (3-2) If ω1>-ω and ω2<ω, the model consists of two parts: the first part calculates the daily average clear sky index; the second part simulates the intermittent hourly changes by superimposing a first-order autoregressive process; k t =k tM +σy(h) where k t is the hourly clear sky index, k tm is the average clear sky index; σ is the standard deviation of the clear sky index per hour on the day corresponding to the clear sky index per day; Wherein the algorithm constants λ, ∈, κ and the solar altitude angle h s They are: ∈＝0.32-1.6(K t -0.5)2 h s =arcsin(cosω2*cosδ*coslat+sinδ*sinlat) According to the study of the sequential characteristics of the y variable, its current hourly value is only significantly dependent on the value of the previous hour, so the autoregressive function y(h) is defined as: In the formula, is the autocorrelation coefficient: r is a randomly generated variable that conforms to the Gaussian distribution, and its expected value is zero: In the formula, z is a floating point number greater than or equal to 0 and less than 1, following a uniform distribution, that is, it randomly falls in the half-open interval [0.0, 1.0), σ′ where the standard deviation σ is: σWhere A and B are algorithm constants: When calculating, if the hourly clear sky index k t <0 or k t >k cs , where k cs The clear sky index for clear skies is: Then regenerate a floating point number z, calculate r, y(h) and k t , and set the maximum number of iterations to 10. If the iterations exceed 10, if k t <0, then k t =0, if k t >k cs , then k t =k cs ; In this way, a clear sky index for each hour of the year can be generated.

4. A method for calculating the total radiation energy of a horizontal plane according to claim 3, characterized in that: The specific calculation steps of step (2) are as follows: (4-1) Calculate the daily clear sky radiation B0d on the outer horizontal surface: Where ω is the sunrise and sunset angle, δ is the solar declination angle, lat is the latitude in radians, and ε is the correction factor for the earth's eccentricity: Where n is the day of the year; (4-2) Calculate the hourly clear sky radiation G0c on the outer horizontal surface: If ω1≤-ω or ω2≥ω, it is non-sunshine time, and its value is 0. If ω1>-ω and ω2<ω, Where ω2 is the sun angle at the middle of each hour, and ω is the sunrise and sunset angle.

5. A device for calculating the total radiation energy of a horizontal surface, characterized in that: The method comprises a memory and a processor, wherein the memory is used to store a computer program, and the processor is used to read the computer program and execute the method for calculating the total radiation energy of a horizontal surface as claimed in any one of claims 1 to 4.

6. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores instructions, and when the instructions are executed on a computer, the method for calculating the total radiation energy of a horizontal surface according to any one of claims 1 to 4 is executed.

Images