Patent Application
20130226537
The present application claims priority under 35 U.S.C 119(a) to Korean Application No. 10-2011-0086810, filed on Aug. 30, 2011, in the Korean Intellectual Property Office, which is incorporated herein by reference in its entirety set forth in full.
Exemplary embodiments relate to a method and an apparatus for providing a solar radiation model, and more particularly, a method and an apparatus for providing a one layer solar radiation model for calculation of insolation.
Solar radiation (insolation) is the primary motive power for moving the atmosphere of the earth and a necessary energy source for all the living things on the earth. The solar radiation energy and a change in the energy are used in various fields, such as not only atmospheric science, but agricultural, energy, biological, medical, and architectural researches and industrial activities.
Equipment for measuring the solar radiation energy is called a pyrheliometer, which measures insolation. The pyrheliometer, however, has a problem in that maintenance and calibration are very difficult and the measurable area is limited in comparison to other atmospheric equipment. Therefore, a solar radiation model is used to check distribution of solar radiation over a wide area.
However, a solar radiation model that performs calculation in consideration of several layers of the atmosphere and the input data of the solar radiation model is performed at limited observatories, such that the availability is very low.
An aspect of the present invention is to calculate the amount of solar energy reaching the surface of the earth with reference to the information on absorbed gases, aerosol, and cloud in the atmosphere by providing a one layer solar radiation model.
Another aspect of the present invention is to quickly and accurately calculate solar energy of a selected spot or a wide area using the data about the ground or the cloud measured by a satellite in consideration of limited features of the input data and based on the features of input data.
In accordance with an exemplary embodiment of the present invention, there is provided a method of providing a one layer solar radiation model for calculation of insolation including: calculating direct solar radiation at a selected position, based on extraterrestrial radiation, eccentricity of the earth from the sun, and transmittance of gases in the atmosphere; calculating final diffuse solar radiation at the selected position, by using diffuse solar radiation due to air molecules, diffuse solar radiation due to aerosol, and multi-diffuse solar radiation with the atmosphere and the surface of the earth; accumulating the calculated direct solar radiation and final diffuse solar radiation for a predetermined wavelength region; calculating the direct solar radiation reaching the surface of the earth and the final diffuse solar radiation reaching the surface of the earth by applying the amount of cloud to the direct solar radiation and the final diffuse solar radiation; and calculating global solar radiation at the selected positions by using the direct solar radiation reaching the surface of the earth and the final diffuse solar radiation reaching the surface of the earth.
In accordance with another exemplary embodiment of the present invention, there is provided an apparatus of providing a one layer solar radiation model for calculation of insolation including: a direct solar radiation unit configured to calculate direct solar radiation at a selected position, based on extraterrestrial radiation, eccentricity of the earth from the sun, and transmittance of gases in the atmosphere; a diffuse solar radiation calculating unit configured to calculate final diffuse solar radiation at the selected position, by using diffuse solar radiation due to air molecules, diffuse solar radiation due to aerosol, and multi-diffuse solar radiation with the atmosphere and the surface of the earth; an accumulating solar radiation unit configured to accumulate the calculated direct solar radiation and final diffuse solar radiation for a predetermined wavelength region; a cloud amount applying unit configured to calculate the direct solar radiation reaching the surface of the earth and the final diffuse solar radiation reaching the surface of the earth by applying the amount of cloud to the direct solar radiation and the final diffuse solar radiation; and a global solar radiation calculating unit configured to calculate global solar radiation at the selected positions by using the direct solar radiation reaching the surface of the earth and the final diffuse solar radiation reaching the surface of the earth.
The above and other aspects, features and other advantages will be more clearly understood from the following detailed description taken in conjunction with the accompanying drawings, in which:
Hereinafter, exemplary embodiments of the present invention will be described in detail with reference to the accompanying drawings. However, the present invention is not limited to the exemplary embodiment. Like reference numerals given in the drawings indicate like components.
As illustrated in
The eccentricity E0 is a change in distance between the earth and the sun generated due to the elliptical revolution of the earth around the sun, which can be calculated depending on the following Equation 1.
E
0=(r0/r)2=1.000110+0.034221 cos Γ+0.001280 sin Γ+0.000719 cos 2Γ+0.000077 sin 2Γ
Γ=2π(dn−1)/365 [Equation 1]
In the above Equation 1, r0 represents an average distance (1 AU=1.496×108 km) between the earth and the sun, r represents an actual distance, and dn represents a day number for one year. The eccentricity is the smallest (about 0.9666) at the summer solstice and the largest (about 1.046) at the winter solstice.
The solar zenith angle is a factor that has the most influence on the intensity of the sunlight on the surface of the earth, under the clear atmosphere without a cloud, which can be calculated depending on the following Equation 2.
cos θ=sin δ sin φ+cos δ cos φ cos α [Equation 2]
In the above Equation 2, δ represents the declination of the sun, φ represents the latitude, and α represents the hour angle. That is, on a sunny day without the cloud, the solar zenith angle or the altitude of the sun is the most important factor of a change of the sunlight that reaches the surface of the earth, which is determined based on the declination and the latitude of the sun and the local apparent time.
The declination δ means the angle between an equatorial plane of the earth which is generated due to the revolution of the earth at an angle of 23.5° and the sun and can be expressed by the following Equation 3.
δ=0.006918−0.399912 cos Γ+0.070257 sin Γ−0.006768 cos 2Γ+0.000907 sin 2Γ−0.002697 cos 3Γ+0.001480 sin 3Γ [Equation 3]
The local apparent time (LAT) means the actual solar time at the selected positions and is expressed by the following Equation 4.
LAT=LST+4(Ls−Le)+Et
E
t=(0.000075+0.001868 cos Γ−0.032077 sin Γ−0.014615 cos 2Γ−0.040849 sin 2Γ)*229.18 [Equation 4]
In the above Equation 4, Ls represents a standard longitude (135° for Korea), Le represents a longitude, and Et represents a time equation. Further, in accordance with an exemplary embodiment of the present invention, the LAT may be converted into universal time coordinated (UTC) and used.
Thereafter, the direct solar radiation is calculated for each wavelength by using the input data information (S120) and the direct solar radiation can be calculated by the following Equation 5.
Idλ=IodE0 cos θtλ [Equation 5]
In the above Equation 5, Idλ represents the direct solar radiation for each wavelength, Ioλ represents extraterrestrial radiation, E0 represents eccentricity showing a distance change between the earth and the sun, θ represents the zenith angle, and tλ represents the transmittance of the absorbed gases.
The transmittance of the absorbed gases can be calculated based on the optical thickness Tλ of each wavelength of the atmospheric components of tλ=exp(−τλ) and can be expressed by multiplying of the transmittance of each component in the atmosphere, based on the following Equation.
tλ=trλtaλtoλtwaλtgλ [Equation 6]
In the above Equation 6, trλ, taλ, toλ, twaλ, and tgλ are the transmittance of air molecules, aerosol, ozone, water vapor, and gas mixture, respectively. The transmittance can be calculated based on the following Equation 7.
t
oλ=exp(−koλlmo)
t
waλ=exp[−0.2385 kwaλwmr/(1+20.07 kwaλwmr)0.45]
t
rλ=exp(−0.008735×ma×λ−4.08)
t
aλ=exp(−β×(λ/0.55)−αmr)
t
gλ=exp[(−1.41kgλma)/(1+118.93kgλma)0.45] [Equation 7]
In the above Equation 7, koλ represents an absorption coefficient for each wavelength for ozone, kgλ represents the absorption coefficient for each wavelength for the gas mixture, kwaλ represents the absorption coefficient for each wavelength for the water vapor, a that is a calibrating constant of the aerosol represents is 1.027 when the wavelength is smaller than 0.5 μm and 1.206 for other wavelengths, l represents the vertical amount of ozone, and ω is the perceptible water.
Further, mr is relative optical mass of the air molecules, ma represents a value with the atmospheric pressure calibrated, and mo represents relative optical mass for the ozone, which are expressed by the following Equation 8.
m
r=[cos θ+0.15(93.885−θ)−1.253]−1
m
a
=m
r(P/1013.25)
m
o=(1+h0/6370)/(cos2θ+2h0/6370)0.5 [Equation 8]
In the above Equation 8, P represents the local atmospheric pressure and h0 represents the altitude for the largest ozone density, which is 22 km herein.
Thereafter, the final diffuse solar radiationdiffuse solar radiation at selected positions is calculated for each wavelength by using diffuse solar radiationdiffuse solar radiation due to the air molecules, diffuse solar radiationdiffuse solar radiation due to the aerosol, and multi-diffuse solar radiationdiffuse solar radiation with the surface of the earth (S130).
In detail, the diffusion insolation Isλ that reaches the surface of the earth represents the sum of the diffusion Irλ due to the air molecules, the diffusion Iaλ due to the aerosol, and the multi-diffusion Igλ with the atmosphere and the surface of the earth, and is expressed by the following Equation 9.
I
sλ
=I
rλ
+I
aλ
+I
gλ [Equation 9]
In the above Equation 9, Isλ represents the final diffusion insolation reaching the surface of the earth, Irλ represents the diffusion insolation due to the air molecules, Iaλ represents the diffusion insolation due to the aerosol, and Igλ represents the multi-diffuse solar radiationdiffuse solar radiation with the surface of the earth.
Further, the diffusion insolation Irλ due to the air molecules, the diffusion insolation Iaλ due to the aerosol, and the multi-diffuse solar radiationdiffuse solar radiation Igλ with the surface of the earth are calculated depending on the following Equation 10.
I
rλ
=I
0λ
E
0 cos θtoλtgλtwaλ×[0.5(1−toλ)taλ]
I
aλ
=I
0λ
E
0 cos θtoλtgλtwaλ×[Fcω0(1−taλ)trλ]
I
gλ
=Q
λ(ρgλρaλ/(1−ρgλρ′aλ))
ρ′aλ=toλtgλtwaλ[0.5(1−trλ)taλ+(1−Fc)ω0(1−taλ)trλ] [Equation 10]
In the above Equation 10, Fc represents a front diffusion ratio in the diffusion direction of the entire energy, 1−Fc represents a rear diffusion ratio, ω0 represents single diffusion albedo, which changes to about 0.6 in a city and about 0.9 in a country, Qλ represents direct and primary diffuse solar radiationdiffuse solar radiation reaching the surface of the earth, ρgλ represents surface albedo, and ρ′gλ represents atmospheric albedo.
Thereafter, the direct solar radiation or the final diffuse solar radiationdiffuse solar radiation, which is calculated for each wavelength, is subjected to accumulated calculation for a predetermined wavelength region (S140).
Thereafter, the calculation of applying the amount of cloud to the direct solar radiation or the final diffuse solar radiation which are subjected to the accumulated calculation is performed (S150). In detail, the insolation on a cloudy day by multiplying the insolation by the amount of cloud may be calculated, and the calculation of applying the amount of cloud to the insolation can be performed depending on the following Equation 11.
Fcloud=(1−cloud amount)×Fclear [Equation 11]
In the above Equation 11, Fcloud represents the insolation in the cloudy state, cloud amount represents the amount of cloud, and Fclear represents the direct solar radiation or the final diffuse solar radiation in the sunny state without cloud.
Thereafter, the global solar radiation is calculated from the direct solar radiation and the final diffuse solar radiation, which are calculated as described above (S160).
The global solar radiation is calculated depending on the following Equation 12.
I=I
dλcos θ+Isλ [Equation 12]
In the above Equation 12, λ represents the wavelength, I represents the global solar radiation, Idλ represents the direct solar radiation for each wavelength, Isλ represents the final diffuse solar radiation, and θ represents the solar zenith angle.
The insolation calculated as described above may be output, if necessary, as the insolation for each wavelength, the direct solar radiation, the final diffuse solar radiation, or the global solar radiation, in accordance with the purpose of use (S170).
As illustrated in
The direct solar radiation calculating unit 210 calculates the direct solar radiation at a selected position, using the extraterrestrial radiation, the eccentricity of the earth from the sun, and the transmittance of the gases in the atmosphere.
The direct solar radiation calculating unit 210 may calculate the direct solar radiation at a plurality of selected positions or may calculate the direct solar radiation at the selected positions for a predetermined period of time.
In this case, the transmittance of the gases in the atmosphere is calculated based on the transmittance of the air molecules, the aerosol, the ozone, the water vapor, and the gas mixture. Further, the extraterrestrial radiation is calculated based on the eccentricity of the earth from the sun, the declination that is an angle between the equatorial plane of the earth and the sun, the local apparent time that is the actual solar time at the selected position, and the solar zenith angle that is an angle between the rotation axis of the earth and the sun.
The diffuse solar radiation calculating unit 220 calculates the final diffuse solar radiation at the selected position, using the diffuse solar radiation due to the air molecules, the diffuse solar radiation due to the aerosol, and the multi-diffuse solar radiation with the atmosphere and the surface of the earth.
The diffuse solar radiation calculating unit 220 may calculate the final diffuse solar radiation at a plurality of selected positions or may calculate the final diffuse solar radiation at the selected positions for a predetermined period of time.
The accumulating solar radiation unit 230 accumulates the calculated direct solar radiation and the final diffuse solar radiation for the predetermined wavelength region.
The cloud amount applying unit 240 calculates the direct solar radiation reaching the surface of the earth and the final diffuse solar radiation reaching the surface of the earth, by applying the amount of cloud to the direct solar radiation and the final diffuse solar radiation. In this case, the cloud amount applying unit 240 may calculate the direct solar radiation and the final diffuse solar radiation that reach the surface of the earth depending on the above Equation 11.
The global solar radiation calculating unit 250 calculates the global solar radiation at the selected position, using the direct solar radiation reaching the surface of the earth and the final diffuse solar radiation reaching the surface of the earth.
In more detail, the global solar radiation calculating unit 250 can calculate the global solar radiation by applying the solar zenith angle, which is the angle between the rotation axis of the earth and the sun, to the direct solar radiation, and then by summing the final diffuse solar radiation for each wavelength.
Meanwhile, a variety of input data such as the ozone, the aerosol, the perceptible water, the surface albedo, the atmospheric pressure, and the atmospheric temperature is required to drive the solar radiation model in accordance with an exemplary embodiment of the present invention. In the data, the data of the atmospheric pressure and the perceptible water is obtained from a regional data assimilation and prediction system (RDAPS) having a resolution of 10 km, the data of the ozone is obtained from the data of an ozone monitoring instrument (OMI) sensor having a resolution of 1°×1°, and the data of the aerosol was obtained from the data generated by a moderate resolution imaging spectroradiometer (MODIS) sensor. The data of the albedo was obtained from the data of 0.05°×0.05° generated by the MODIS sensor and the data of the altitude of the calculated positions was obtained from the 3 seconds (90 m×90 m) data of the shuttle radar topographic mission (SRTM). Further, the insolation with a resolution of 4 km was calculated, when the sky of the Korean peninsula was sunny and cloudy, by using a GWNU model and data produced on the basis of the solar radiation theory described above.
In accordance with the exemplary embodiments of the present invention, it is possible to calculate the amount of solar energy reaching the surface of the earth with reference to the information on absorbed gases, aerosol, and cloud in the atmosphere by providing the one layer solar radiation model.
Further, in accordance with the exemplary embodiments of the present invention, it is possible to quickly and accurately calculate the solar energy of a selected spot or a wide area using the data about the ground or the cloud measured by a satellite in consideration of limited features of the input data and based on the features of input data and use the calculated solar energy for the research of insolation and remote sensing.
In addition, in accordance with the exemplary embodiments of the present invention, it is possible to use the one layer solar radiation model for the basis research result of the studies, such as researching into energy development, radiation measurement, climate data processing, and the like.
Although the present invention is described by specific matters such as concrete components, and the like, embodiments, and drawings, they are provided only for assisting in the entire understanding of the present invention. Therefore, the present invention is not limited to the embodiments. Various modifications and changes may be made by those skilled in the art to which the present invention pertains from this description. Therefore, the spirit of the present invention should not be limited to the above-described embodiments, and the following claims as well as all modified equally or equivalently to the claims are intended to fall within the scope and spirit of the invention.
| Number | Date | Country | Kind |
|---|---|---|---|
| 10-2011-0086810 | Aug 2011 | KR | national |
