CN110378028A - A kind of random solar irradiation intensity distribution calculation method of horizontal plane - Google Patents
A kind of random solar irradiation intensity distribution calculation method of horizontal plane Download PDFInfo
- Publication number
- CN110378028A CN110378028A CN201910660559.XA CN201910660559A CN110378028A CN 110378028 A CN110378028 A CN 110378028A CN 201910660559 A CN201910660559 A CN 201910660559A CN 110378028 A CN110378028 A CN 110378028A
- Authority
- CN
- China
- Prior art keywords
- irradiation
- horizontal plane
- day
- random
- average
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Abstract
The invention discloses a kind of random solar irradiation intensity distribution calculation methods of horizontal plane, according to horizontal plane day sun total amount of irradiation, establish horizontal plane monthly average scattering irradiation and the per day sun total amount of irradiation linear relationship of each weather condition lower horizontal plane, it is final to generate minute rank monthly mean level face always irradiation, random scatter irradiation and random direct projection irradiation profile curve at random in conjunction with sunny condition lower horizontal plane solar irradiation distribution calculation method and probability distribution function.The present invention accurate can calculate the random process that solar irradiation reaches ground under the influence of multiple factors, realize reasonable effective assessment to photovoltaic plant year irradiation.
Description
Technical field
The present invention relates to a kind of random solar irradiation intensity distribution calculation methods of horizontal plane, belong to solar energy photovoltaic system hair
Electricity assessment technology field.
Background technique
During planning and designing photovoltaic plant, photovoltaic plant annual electricity generating capacity is very important an index, comments
The quality of valence depends on the year irradiation for being irradiated to photovoltaic module surface.Therefore, it is necessary to one kind rationally can effectively assess irradiation
To photovoltaic module surface year irradiation method.
Currently, with good grounds typical meteorological year method generates typical day irradiation profile, or straight in the method for assessment year irradiation
It connects and irradiation profile is obtained according to year horizontal irradiation and transformation ratio product, these methods have ignored solar irradiation distribution by day
The multiple factors such as text, physics, meteorology influence, and solar irradiation is caused to eventually arrive at a case where random process is presented in ground.And
Common Metonorm meteorology software can only also provide main year, month, day and hour irradiation data, still not for minute rank
Data can be provided.
Summary of the invention
The technical problem to be solved by the present invention is to overcome the deficiencies of the prior art and provide a kind of random sun spoke of horizontal plane
According to the calculation method of intensity distribution, it is based on horizontal plane day sun total amount of irradiation, establishes horizontal plane monthly average scattering irradiation and each
The per day sun total amount of irradiation linear relationship of weather condition lower horizontal plane, and sunny condition lower horizontal plane solar irradiation is combined to be distributed
Calculation method and probability distribution function, it is final to generate minute rank monthly mean level face always irradiation, random scatter irradiation at random
And random direct projection irradiation profile curve.
In order to solve the above technical problems, the present invention provides a kind of random solar irradiation intensity distribution calculation method of horizontal plane,
The following steps are included:
(1) it collects to calculate and place horizontal plane day adds up total amount of irradiation historical data, calculate the horizontal plane of place in the neighbourhood
Day accumulative total amount of irradiation historical data and horizontal plane day accumulative scattering irradiation historical data;
(2) weather typing that total amount of irradiation carries out fine day, broken sky and rainy days is added up according to horizontal plane day;
(3) place accumulative average total irradiation of each weather typing condition lower horizontal plane day per month in the neighbourhood is calculated
Amount, and horizontal plane day accumulative average scattering irradiation monthly;Calculating place, each weather typing condition is lauched per month
Add up average total amount of irradiation plane day;
(4) it establishes about horizontal plane day accumulative average total amount of irradiation and horizontal plane day adds up the linear of average scattering irradiation
Regression equation, and solve to calculate and add up average scattering irradiation the horizontal plane day of place monthly;
(5) the ideal irradiation profile of typical day is calculated;
(6) handled using gaussian random and sliding average, obtain minute rank monthly average day horizontal plane at random always irradiation with
Random scatter irradiation profile curve.
The time span of the historical data of collection above-mentioned was at 1 year or more.
The acquisition interval of historical data above-mentioned keeps the frequency of minute grade.
Weather typing above-mentioned includes:
It extracts and monthly horizontal plane day adds up total amount of irradiation maximum value Gm,maxWith minimum value Gm,min;
Third by stages [Gm,min,Gm,max];
Horizontal plane day add up total amount of irradiation and is located at section [Gm,min,Gm,min+(Gm,max-Gm,min)/3] in be rainy days;
Horizontal plane day add up total amount of irradiation and is located at section [Gm,min+(Gm,max-Gm,min)/3,Gm,min+(Gm,max-Gm,min)2/
It 3] is broken sky in;
Horizontal plane day add up total amount of irradiation and is located at section [Gm,min+(Gm,max-Gm,min)2/3,Gm,max] in be fine day.
If historical data time span above-mentioned is greater than 1 year, selected monthly horizontal plane day adds up total amount of irradiation
Maximum value and minimum value are maximum value and minimum value in all same Januaries.
Equation of linear regression above-mentioned are as follows:
Ax+c=B2
Wherein, A is parameter matrix, and x is metewand phasor, and c is constant term, B2For parameter phasor;
Place will be calculated, and accumulative average total amount of irradiation of each weather typing condition lower horizontal plane day substitutes into parameter square in the neighbourhood
Battle array A will be calculated and be added up average scattering irradiation substitution parameter phasor B the horizontal plane day of place in the neighbourhood monthly2;
It is solved by Moore-Penrose inverse matrix, obtains the least square solution of equation,
Obtain metewand phasor x:
X=A+(B2-c);
Place each weather typing condition lower horizontal plane day accumulative average total amount of irradiation will be calculated and substitute into parameter matrix A, in conjunction with
The metewand phasor x of solution is solved and the monthly horizontal plane day for obtaining calculating place is added up average scattering irradiation.
The ideal irradiation profile above-mentioned for calculating typical day, comprising:
Monthly typical day as the monthly average day proposed using Klein, according to Haurwitz empirical model as ideal weather
Situation lower horizontal plane always irradiates computation model, calculation formula after amendment are as follows:
Wherein, R is correction factor, θZFor solar zenith angle, IGIt is always irradiated by force for the ideal weather condition lower horizontal plane sun
Degree;
Correction factor is solved by following formula integral calculation according to the ideal total irradiation intensity of the weather condition lower horizontal plane sun:
Typical day maximum, average and minimum irradiation item are arrived by continually changing solar zenith angle based on formula (1)
Ideal irradiation profile under part.
Monthly average day horizontal plane above-mentioned always irradiates randomization, comprising:
It is proposed three assumed conditions: (a1) is inscribed in t0, irradiation value Normal Distribution;(a2) normal distribution mean value takes
The value inscribed for C2 in t0;(a3) 3 times of standard deviation ranges of normal distribution mean value ± are taken, standard deviation is taken as inscribing (C1- when t0
C2)/3 with (C2-C1)/3 in absolute value smaller value;
Utilize taken standard deviation range, mean value and standard deviation, using the method for gaussian random obtain corresponding Gauss with
Machine number ultimately generates the irradiation profile curve generated after the total irradiation profile curve in typical day Random Level face and sliding average;
Wherein, C1 is the ideal irradiation profile curve under the conditions of typical day maximum permissible exposure, and C2 is typical day monthly average spoke
According to the ideal irradiation profile curve under the conditions of amount, C3 is the ideal irradiation profile curve under the conditions of typical day minimum irradiation.
Horizontal area scattering irradiation randomization above-mentioned, comprising:
The horizontal plane day of the calculating place calculated according to step (4) monthly adds up average scattering irradiation, in conjunction with
Haurwitz empirical model obtains revised typical day horizontal area scattering irradiation profile curve D1, and assumes: (b1) is in t0When
It inscribes, irradiation value Normal Distribution;(b2) normal distribution mean value is taken as in t0When inscribe the irradiation value of D1;(b3) standard deviation takes
For t0When inscribe absolute value smaller value in (C1-D1)/3 and (D1-C3)/3;
According to above-mentioned it is assumed that the typical day Random Level area scattering irradiation profile curve and cunning that output gaussian random generates
The irradiation profile curve generated after moving averagely.
Sliding average above-mentioned, comprising: use different data volumes as exposure mask size respectively, carry out sliding average and obtain
Random irradiation profile curve.
The beneficial effects obtained by the present invention are as follows are as follows:
The present invention accurate can calculate the random process that solar irradiation reaches ground under the influence of multiple factors, real
Now to reasonable effective assessment of photovoltaic plant year irradiation.
Detailed description of the invention
Fig. 1 is flow chart of the method for the present invention.
Fig. 2 is that certain in the embodiment of the present invention calculates the typical total irradiation profile curve in day Random Level face of place in August, 2018
Figure;
Fig. 3 is the total irradiation profile curve graph in typical day Random Level face in the embodiment of the present invention after sliding average;
Fig. 4 is that certain in the embodiment of the present invention calculates the typical day Random Level area scattering irradiation profile song of place in August, 2018
Line.
Specific embodiment
The invention will be further described below in conjunction with the accompanying drawings.Following instance is only used for clearly illustrating of the invention
Technical solution, herein, these embodiments of the invention individually or can be indicated generally with term " invention ", this
Just for the sake of convenient, and not intended to limit the protection scope of the present invention.
In the present invention, it is assumed that the solar irradiation for reaching level ground only scatters solar irradiation and direct projection solar irradiation, no
The incident irradiance of other forms is considered, so monthly mean level total amount of irradiation scatters irradiation and monthly average water by monthly mean level
It is straight to penetrate irradiation composition, and either fine day, have skies gas or a rainy weather, scattering irradiation all exists with certain proportion.
The present invention provides a kind of random solar irradiation intensity distribution calculation method of horizontal plane, referring to Fig. 1, including following step
It is rapid:
(1) the horizontal plane day accumulative total amount of irradiation historical data for calculating place is collected from weather site or experiment porch, with
And it calculates and adds up total amount of irradiation historical data and horizontal plane day accumulative scattering irradiation history the horizontal plane day of place in the neighbourhood
Data;
In the embodiment of the present invention, the time span Ying Yi or more of the historical data of collection, day acquisition time Ying Yougu
Determine range, and the interval of collection point should keep the frequency of minute grade, guarantee the integrality of data as far as possible.This method will utilize number
Calculating analysis is carried out according to more complete number of days.
(2) interval division is carried out to horizontal plane day accumulative total amount of irradiation, weather is divided into fine day, broken sky and rainy days.
Interval division is carried out to total amount of irradiation is added up horizontal plane day, comprising:
It extracts and middle of each month horizontal plane day adds up total amount of irradiation maximum value Gm,maxWith minimum value Gm,min;
By [Gm,min,Gm,max] section progress trisection, i.e., all horizontal plane days monthly add up total amount of irradiation and is all being spaced
Within the scope of;
It will add up total amount of irradiation horizontal plane day and be located at section [Gm,min,Gm,min+(Gm,max-Gm,min)/3] in be divided into yin
Rainy day;
It will add up total amount of irradiation horizontal plane day and be located at section [Gm,min+(Gm,max-Gm,min)/3,Gm,min+(Gm,max-Gm,min)
2/3] broken sky is divided into;
It will add up total amount of irradiation horizontal plane day and be located at section [Gm,min+(Gm,max-Gm,min)2/3,Gm,max] in be divided into it is fine
It.
If historical data time span is greater than 1 year, the maximum value and minimum value in all January are chosen at this time.
(3) place accumulative average total irradiation of each weather typing condition lower horizontal plane day per month in the neighbourhood is calculated
Amount, and horizontal plane day accumulative average scattering irradiation monthly;Calculating place, each weather typing condition is lauched per month
Add up average total amount of irradiation plane day.
Horizontal plane day accumulative average total amount of irradiation under the conditions of fine day, broken sky and rainy days weather typing is denoted as respectively
a1j、a2jWith a3j;Wherein, a1jIndicate the horizontal plane day accumulative average total amount of irradiation under j-th month sunny weather, a2jIndicate jth
Horizontal plane day accumulative average total amount of irradiation under a month broken sky weather, a3jIndicate the horizontal plane under j-th month rainy days of weather
Add up average total amount of irradiation day.If historical data time span is greater than 1 year, j-th month day, accumulative average total amount of irradiation was taken
All j-th month statistical averages.
It calculates the horizontal plane day monthly of place in the neighbourhood accumulative average scattering irradiation and is denoted as b2j, indicate j-th month
Horizontal plane day adds up average scattering irradiation.
(4) it establishes horizontal plane day accumulative average total amount of irradiation and adds up the linear regression of average scattering irradiation horizontal plane day
Equation solves and calculates the horizontal plane day of place monthly accumulative average scattering irradiation.
Linear regression equation is as follows:
Ax+c=B2
Wherein, A is parameter matrix, and x is metewand phasor, and c is constant term, B2For parameter phasor.
Add up average total amount of irradiation each weather typing condition lower horizontal plane day per month in the neighbourhood according to place is calculated, with
And horizontal plane day monthly adds up average scattering irradiation, solves metewand phasor, it may be assumed that
The horizontal plane day of place in the neighbourhood accumulative average total amount of irradiation a will be calculated1j、a2jWith a3jParameter matrix A is substituted into,
Parameter matrix A is made of 3 column, 12 row, and 3 column are respectively a1, a2, a3, the data of 12 behaviors 12 months.Place will be calculated nearby
Add up average scattering irradiation the horizontal plane day of area monthly and substitutes into parameter phasor B2;
It is solved by Moore-Penrose inverse matrix, obtains the least square solution of equation to get metewand phasor x is arrived:
X=A+(B2-c)。
By the horizontal plane day for calculating place accumulative average total amount of irradiation a1j、a2jWith a3jBring equation Ax+c=B into2Parameter
In matrix A, in conjunction with the metewand phasor x of solution, solves and add up average scattering irradiation the horizontal plane day moon for obtaining calculating place
Amount, is denoted as B1。
(5) the ideal irradiation profile of typical day is calculated;
Monthly typical day as the monthly average day proposed using Klein, typical day are in every month closest to average declination angle
One day, and computation model is always irradiated as ideal weather condition lower horizontal plane according to Haurwitz empirical model, after amendment
Calculation formula are as follows:
Wherein, R is correction factor, θZFor solar zenith angle, IGIt is always irradiated by force for the ideal weather condition lower horizontal plane sun
Degree.
Correction factor is solved by following formula integral calculation according to the ideal total irradiation intensity of the weather condition lower horizontal plane sun:
Based on formula (1) by continually changing solar zenith angle, typical day maximum, average and minimum irradiation can be obtained
Under the conditions of ideal irradiation profile.Enabling the ideal irradiation profile curve under the conditions of typical day maximum permissible exposure is C1, monthly average spoke
It is C3 under the conditions of minimum irradiation according to being C2 under the conditions of amount.
(6) handled using gaussian random and sliding average, obtain minute rank monthly average day horizontal plane at random always irradiation with
Random scatter irradiation profile curve.
It is near random when reaching ground because solar irradiation is influenced by various factors under physical condition, the present invention
To the randomization that horizontal plane always irradiates, three assumed conditions are proposed: (1) in t0When inscribe, irradiating the value that is likely to occur and obeying just
State distribution;(2) normal distribution mean value is taken as C2 in t0When the value inscribed, make to irradiate overall trend and average lamp phase respectively at random
Together;(3) according to 3 times of standard deviation ranges of normal distribution mean value ±, for probability close to 100%, standard deviation is taken as t0When inscribe
(C1-C2)/3 with the absolute value smaller value in (C2-C1)/3, a possibility that overflowing maximum Distribution value and minimum value is reduced.
According to above 3 points it is assumed that utilizing taken data area (± 3 times of standard deviation ranges), mean value, standard deviation, use
The method of gaussian random obtains corresponding Gauss number, ultimately generates the total irradiation profile curve in typical day Random Level face, bent
Line has fluctuating characteristic and curvilinear trend is consistent with Average profiles.
Because the curve data fluctuation that the method for gaussian random generates is larger, the present invention continues the method pair using sliding average
Random irradiation profile curve smoothing processing is obtained, i.e., uses different data volumes as exposure mask size respectively, carries out sliding average
Random irradiation profile curve is obtained, the curve generated is gentler while guaranteeing certain fluctuating characteristic.
The present invention irradiates horizontal area scattering be randomized simultaneously, according to the moon horizontal plane in the calculating place that step (4) calculate
Accumulative average scattering irradiation obtained the horizontal area scattering irradiation of revised typical day and divided in conjunction with Haurwitz empirical model day
Cloth curve D1.And its basic assumption is revised are as follows: (1) in t0When inscribe, irradiate the value Normal Distribution being likely to occur;(2)
Normal distribution mean value is taken as in t0When inscribe the irradiation value of D1;(3) standard deviation is taken as t0When inscribe (C1-D1)/3 and (D1-C3)/3
In absolute value smaller value.
According to above 3 points it is assumed that output gaussian random generate typical day Random Level area scattering irradiation profile curve with
And to the irradiation profile curve generated after its sliding average.
Embodiment
The present invention is set to 32.49 degree of north latitude by contraposition, and 119.9 degree of east longitude, the time is from September, 2017 in August, 2018
7 points of daily morning at 16 points in afternoon, and with the frequency of a collection point per minute, the total irradiation intensity of horizontal plane, scattering are irradiated
Intensity, direction direct projection irradiation intensity are monitored.
The present invention, which extracts to calculate, adds up total amount of irradiation historical data place horizontal plane day, calculates the level of place in the neighbourhood
The accumulative total amount of irradiation historical data in face day and horizontal plane day accumulative scattering irradiation historical data, horizontal plane day monthly is added up
Total amount of irradiation maximum value is denoted as Gm,max, horizontal plane day monthly add up total amount of irradiation minimum value and be denoted as Gm,min, third by stages
[Gm,min,Gm,max], wherein [Gm,min,Gm,min+(Gm,max-Gm,min)/3] it is divided into rainy days;[Gm,min+(Gm,max-Gm,min)/
3,Gm,min+(Gm,max-Gm,min) 2/3] it is divided into broken sky;[Gm,min+(Gm,max-Gm,min)2/3,Gm,max] it is divided into fine day.
Day under the conditions of fine day, broken sky, rainy days weather typing monthly is added up total amount of irradiation average value to be denoted as respectively
a1j、a2jWith a3j;It calculates the horizontal plane day of place in the neighbourhood accumulative scattering irradiation monthly average value and is denoted as b1j.Building is shaped like ' Ax
The matrix equation of+c=b1 ', and contain undetermined constant item c and metewand phasor x in equation.Present invention combination Moore-
Penrose inverse matrix method for solving obtains the mould least square solution of equation, as ' x=A+(b1- c) ', solution obtains a1~a3
Metewand phasor x be -0.0659,0.3573,0.1591, constant term coefficient be 0.002912, horizontal plane day is established with this
Accumulative scattering irradiation monthly average value and each weather condition lower horizontal plane day add up the linear relationship of total amount of irradiation average value.It will meter
Calculate a in place1j、a2jWith a3jIt substitutes into parameter matrix A, in conjunction with metewand phasor x and constant term coefficient c that solution obtains, obtains
To the horizontal plane day accumulative scattering irradiation monthly average value for calculating place.
The present invention is using calculating simple Haurwitz empirical model, for calculating the ideal spoke under the conditions of given irradiation
According to distribution, correction factor R, final formula (1) amendment are obtained by formula (1), (2) integral calculation are as follows:
The monthly typical day that the present invention is proposed according to Klein in 1977 et al. as monthly average day, is in August, 2018
Example, typical day are August 16th, and day issue is 228, test in the total irradiance data of horizontal plane of acquisition missing 7,8,10,20 and
Data on the 27th, the data of remaining 26 days of monitoring, the average value of day horizontal total amount of irradiation is 4.63, maximum value 6.80, minimum
Value is 0.65.By formula (1) and (2), calculate the irradiation profile of typical one day day, can obtain typical day as shown in Figure 2 it is maximum,
Ideal irradiation profile curve under the conditions of average and minimum irradiation, the ideal irradiation in figure under the conditions of typical day maximum permissible exposure
Distribution curve is C1, is C2 under the conditions of monthly average irradiation, is C3 under the conditions of minimum irradiation.
Solar irradiation is influenced by various factors under physical condition, near random when reaching ground, the present invention in order to
This random process is generated, for horizontal total irradiation randomization, three assumed conditions is proposed: (1) being inscribed in t0, irradiation may
The value Normal Distribution of appearance;(2) normal distribution mean value is taken as C2 in t0When the value inscribed, make to irradiate at random respectively it is whole become
Gesture is identical as average lamp;(3) according to 3 times of standard deviation ranges of normal distribution mean value ±, probability is close to 100%, standard deviation
It is taken as t0When inscribe absolute value smaller value in (C1-C2)/3 and (C2-C1)/3, reduce and overflow maximum Distribution value and minimum value
Possibility.
According to above-mentioned it is assumed that the available typical total irradiation profile in day Random Level face of in Augusts, 2018 as shown in Figure 2
Figure, the distribution curve trend is consistent with Average profiles, but its data fluctuations is larger, therefore this method passes through sliding average
Curve is carried out smooth.It is illustrated in figure 3 and 7,9,11 data volumes is used to obtain as exposure mask size progress sliding average respectively
Random irradiation profile curve, three curves are more gentle compared with the curve that gaussian random directly generates, while ensure that certain
Fluctuating characteristic.
The present invention uses method same as described above, the average scattering irradiation of typical day to the random irradiation profile of scattering is generated
For the curve D1 corrected according to the linear relationship of above-mentioned foundation and Haurwitz empirical model, basic assumption revision are as follows:
(1) in t0When inscribe, irradiate the value Normal Distribution being likely to occur;(2) normal distribution mean value is taken as in t0When inscribe the spoke of D1
According to value;(3) standard deviation is taken as t0When inscribe absolute value smaller value in (C1-D1)/3 and (D1-C3)/3.
According to above-mentioned it is assumed that the Random Level that the available typical day gaussian random of in Augusts, 2018 as shown in Figure 4 generates
Area scattering irradiation profile and to the irradiation profile generated after its sliding average.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, without departing from the technical principles of the invention, several improvement and deformations can also be made, these improvement and deformations
Also it should be regarded as protection scope of the present invention.
Claims (10)
1. a kind of random solar irradiation intensity distribution calculation method of horizontal plane, which comprises the following steps:
(1) it collects to calculate and place horizontal plane day adds up total amount of irradiation historical data, it is tired to calculate the horizontal plane day of place in the neighbourhood
Count total amount of irradiation historical data and horizontal plane day accumulative scattering irradiation historical data;
(2) weather typing that total amount of irradiation carries out fine day, broken sky and rainy days is added up according to horizontal plane day;
(3) it calculates place and adds up average total amount of irradiation each weather typing condition lower horizontal plane day per month in the neighbourhood, with
And horizontal plane day monthly adds up average scattering irradiation;Calculate place middle of each month each weather typing condition lower horizontal plane day
Accumulative average total amount of irradiation;
(4) establish about horizontal plane day accumulative average total amount of irradiation and add up horizontal plane day the linear regression of average scattering irradiation
Equation, and solve to calculate and add up average scattering irradiation the horizontal plane day of place monthly;
(5) the ideal irradiation profile of typical day is calculated;
(6) handled using gaussian random and sliding average, obtain minute rank monthly average day horizontal plane at random always irradiation with it is random
Scatter irradiation profile curve.
2. the random solar irradiation intensity distribution calculation method of a kind of horizontal plane according to claim 1, which is characterized in that institute
The time span of the historical data of collection is stated at 1 year or more.
3. the random solar irradiation intensity distribution calculation method of a kind of horizontal plane according to claim 2, which is characterized in that institute
The acquisition interval for stating historical data keeps the frequency of minute grade.
4. the random solar irradiation intensity distribution calculation method of a kind of horizontal plane according to claim 1, which is characterized in that institute
Stating weather typing includes:
It extracts and monthly horizontal plane day adds up total amount of irradiation maximum value Gm,maxWith minimum value Gm,min;
Third by stages [Gm,min,Gm,max];
Horizontal plane day add up total amount of irradiation and is located at section [Gm,min,Gm,min+(Gm,max-Gm,min)/3] in be rainy days;
Horizontal plane day add up total amount of irradiation and is located at section [Gm,min+(Gm,max-Gm,min)/3,Gm,min+(Gm,max-Gm,min) 2/3] in
It is broken sky;
Horizontal plane day add up total amount of irradiation and is located at section [Gm,min+(Gm,max-Gm,min)2/3,Gm,max] in be fine day.
5. the random solar irradiation intensity distribution calculation method of a kind of horizontal plane according to claim 4, which is characterized in that such as
Fruit historical data time span is greater than 1 year, then selected monthly horizontal plane day adds up total amount of irradiation maximum value and minimum value is
Maximum value and minimum value in all same Januaries.
6. the random solar irradiation intensity distribution calculation method of a kind of horizontal plane according to claim 1, which is characterized in that institute
State equation of linear regression are as follows:
Ax+c=B2
Wherein, A is parameter matrix, and x is metewand phasor, and c is constant term, B2For parameter phasor;
Place will be calculated, and accumulative average total amount of irradiation of each weather typing condition lower horizontal plane day substitutes into parameter matrix A in the neighbourhood,
It will calculate and add up average scattering irradiation substitution parameter phasor B the horizontal plane day of place in the neighbourhood monthly2;
It is solved by Moore-Penrose inverse matrix, obtains the least square solution of equation,
Obtain metewand phasor x:
X=A+(B2-c);
Place each weather typing condition lower horizontal plane day accumulative average total amount of irradiation will be calculated and substitute into parameter matrix A, in conjunction with solution
Metewand phasor x, solve obtain calculate place monthly horizontal plane day add up average scattering irradiation.
7. the random solar irradiation intensity distribution calculation method of a kind of horizontal plane according to claim 1, which is characterized in that institute
State the ideal irradiation profile for calculating typical day, comprising:
Monthly typical day as the monthly average day proposed using Klein, according to Haurwitz empirical model as ideal weather condition
Lower horizontal plane always irradiates computation model, calculation formula after amendment are as follows:
Wherein, R is correction factor, θZFor solar zenith angle, IGFor the ideal total irradiation intensity of the weather condition lower horizontal plane sun;
Correction factor is solved by following formula integral calculation according to the ideal total irradiation intensity of the weather condition lower horizontal plane sun:
Based on formula (1) by continually changing solar zenith angle to get to typical day maximum, averagely and under the conditions of minimum irradiation
Ideal irradiation profile.
8. the random solar irradiation intensity distribution calculation method of a kind of horizontal plane according to claim 1, which is characterized in that the moon
Average day horizontal plane always irradiates randomization, comprising:
It is proposed three assumed conditions: (a1) is inscribed in t0, irradiation value Normal Distribution;(a2) normal distribution mean value is taken as C2
The value inscribed in t0;(a3) 3 times of standard deviation ranges of normal distribution mean value ± are taken, standard deviation is taken as inscribing (C1-C2)/3 when t0
With the absolute value smaller value in (C2-C1)/3;
Taken standard deviation range, mean value and standard deviation are utilized, corresponding Gauss number is obtained using the method for gaussian random,
Ultimately generate the irradiation profile curve generated after the total irradiation profile curve in typical day Random Level face and sliding average;
Wherein, C1 is the ideal irradiation profile curve under the conditions of typical day maximum permissible exposure, and C2 is typical life average lamp amount
Under the conditions of ideal irradiation profile curve, C3 is the ideal irradiation profile curve under the conditions of typical day minimum irradiation.
9. the random solar irradiation intensity distribution calculation method of a kind of horizontal plane according to claim 1, which is characterized in that water
In-plane scatter irradiation randomization, comprising:
The horizontal plane day of the calculating place calculated according to step (4) monthly adds up average scattering irradiation, passes through in conjunction with Haurwitz
Model is tested, obtain revised typical day horizontal area scattering irradiation profile curve D1, and assume: (b1) is in t0When inscribe, irradiate
It is worth Normal Distribution;(b2) normal distribution mean value is taken as in t0When inscribe the irradiation value of D1;(b3) standard deviation is taken as t0When inscribe
(C1-D1)/3 with (D1-C3)/3 in absolute value smaller value;
According to above-mentioned it is assumed that typical day Random Level area scattering irradiation profile curve and sliding that output gaussian random generates are flat
The irradiation profile curve generated afterwards.
10. a kind of random solar irradiation intensity distribution calculation method of horizontal plane according to claim 8 or claim 9, feature exist
In the sliding average, comprising: use different data volumes as exposure mask size respectively, carry out sliding average and obtain random spoke
According to distribution curve.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910660559.XA CN110378028B (en) | 2019-07-22 | 2019-07-22 | Horizontal plane random solar irradiation intensity distribution calculation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910660559.XA CN110378028B (en) | 2019-07-22 | 2019-07-22 | Horizontal plane random solar irradiation intensity distribution calculation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110378028A true CN110378028A (en) | 2019-10-25 |
CN110378028B CN110378028B (en) | 2022-08-16 |
Family
ID=68254675
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910660559.XA Active CN110378028B (en) | 2019-07-22 | 2019-07-22 | Horizontal plane random solar irradiation intensity distribution calculation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110378028B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113916369A (en) * | 2021-10-09 | 2022-01-11 | 中国科学院西北生态环境资源研究院 | Snow pile solar radiation acquisition method and system |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763649A (en) * | 2018-04-27 | 2018-11-06 | 河海大学常州校区 | A kind of method of the received irradiation of optimized evaluation photovoltaic module cell piece |
CN109344491A (en) * | 2018-09-27 | 2019-02-15 | 河北工业大学 | A kind of solar irradiance modeling method considering weather state change and cloud cover |
-
2019
- 2019-07-22 CN CN201910660559.XA patent/CN110378028B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763649A (en) * | 2018-04-27 | 2018-11-06 | 河海大学常州校区 | A kind of method of the received irradiation of optimized evaluation photovoltaic module cell piece |
CN109344491A (en) * | 2018-09-27 | 2019-02-15 | 河北工业大学 | A kind of solar irradiance modeling method considering weather state change and cloud cover |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113916369A (en) * | 2021-10-09 | 2022-01-11 | 中国科学院西北生态环境资源研究院 | Snow pile solar radiation acquisition method and system |
Also Published As
Publication number | Publication date |
---|---|
CN110378028B (en) | 2022-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103955749B (en) | Direct solar radiation value Forecasting Methodology and system | |
CN113496311A (en) | Photovoltaic power station generated power prediction method and system | |
CN106021934B (en) | A kind of region can utilize solar energy resources appraisal procedure | |
CN106600060A (en) | Method for predicting amount of solar radiation based on similar day sunny coefficient correction | |
CN108763649B (en) | Method for optimizing and evaluating irradiation quantity received by photovoltaic module cell | |
Safaraliev et al. | Energy potential estimation of the region’s solar radiation using a solar tracker | |
Picotti et al. | Modelling the soiling of heliostats: Assessment of the optical efficiency and impact of cleaning operations | |
Myers et al. | Proposed reference spectral irradiance standards to improve concentrating photovoltaic system design and performance evaluation | |
CN110378028A (en) | A kind of random solar irradiation intensity distribution calculation method of horizontal plane | |
CN107688006B (en) | Quantification method for spectral distribution difference of different climate areas | |
CN104794545A (en) | Optimal planning method of photovoltaic grid-connecting scheme power distribution network | |
Abidin et al. | Optimal efficient energy production by PV module tilt-orientation prediction without compromising crop-light demands in Agrivoltaic systems | |
Chicco et al. | Characterization of solar irradiance profiles for photovoltaic system studies through data rescaling in time and amplitude | |
Xu et al. | Measurement of surface ozone and its precursors in urban and rural sites in Beijing | |
Belkilani et al. | Estimation and experimental evaluation of the shortfall of photovoltaic plants in Tunisia: case study of the use of titled surfaces | |
Castellano et al. | Photosynthetic photon flux density distribution inside photovoltaic greenhouses, numerical simulation, and experimental results | |
Wang et al. | Method for short-term photovoltaic generation power prediction base on weather patterns | |
Xiao et al. | A statistical photovoltaic power forecast model (spf) based on historical power and weather data | |
Meshyk et al. | Climate Resource Potential to Develop Solar Power in Belarus | |
CN106156453B (en) | Solar energy resource assessment method based on numerical weather forecast data | |
CN110390481A (en) | A kind of horizontal plane sun scattering irradiation appraisal procedure and device | |
Maru et al. | Model based optimization of tilt angle for solar PV panels in Jodhpur | |
Debrin et al. | Determining the factors affecting the characteristics of photovoltaic modules in the conditions of Krasnoyarsk Krai | |
CN105956685A (en) | Photovoltaic power factor table prediction method | |
Medykovskyy et al. | Processing of data on the intensity of solar radiation for solar power plant management systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |