CN114692064B - Photovoltaic module service life calculation method based on environmental parameters - Google Patents
Photovoltaic module service life calculation method based on environmental parameters Download PDFInfo
- Publication number
- CN114692064B CN114692064B CN202210363024.8A CN202210363024A CN114692064B CN 114692064 B CN114692064 B CN 114692064B CN 202210363024 A CN202210363024 A CN 202210363024A CN 114692064 B CN114692064 B CN 114692064B
- Authority
- CN
- China
- Prior art keywords
- photovoltaic module
- photovoltaic
- output power
- annual
- value
- 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.)
- Active
Links
- 230000007613 environmental effect Effects 0.000 title claims abstract description 56
- 238000004364 calculation method Methods 0.000 title claims abstract description 34
- 230000002159 abnormal effect Effects 0.000 claims abstract description 46
- 230000005855 radiation Effects 0.000 claims abstract description 22
- 238000005070 sampling Methods 0.000 claims description 38
- 230000008030 elimination Effects 0.000 claims description 29
- 238000003379 elimination reaction Methods 0.000 claims description 29
- 238000003491 array Methods 0.000 claims description 4
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 230000001351 cycling effect Effects 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 19
- 230000032683 aging Effects 0.000 abstract description 9
- 238000010606 normalization Methods 0.000 abstract 2
- 238000012360 testing method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 2
- 239000000306 component Substances 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000013082 photovoltaic technology Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02S—GENERATION OF ELECTRIC POWER BY CONVERSION OF INFRARED RADIATION, VISIBLE LIGHT OR ULTRAVIOLET LIGHT, e.g. USING PHOTOVOLTAIC [PV] MODULES
- H02S40/00—Components or accessories in combination with PV modules, not provided for in groups H02S10/00 - H02S30/00
- H02S40/30—Electrical components
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02S—GENERATION OF ELECTRIC POWER BY CONVERSION OF INFRARED RADIATION, VISIBLE LIGHT OR ULTRAVIOLET LIGHT, e.g. USING PHOTOVOLTAIC [PV] MODULES
- H02S50/00—Monitoring or testing of PV systems, e.g. load balancing or fault identification
- H02S50/10—Testing of PV devices, e.g. of PV modules or single PV cells
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Photovoltaic Devices (AREA)
Abstract
The invention discloses a photovoltaic module service life calculation method based on environmental parameters, and belongs to the technical field of photovoltaic modules. Firstly, collecting the environmental temperature, the environmental humidity, the instantaneous irradiation and the ultraviolet radiation integral value of a photovoltaic module during operation within one year, and removing abnormal values; calculating the voltage and the current of the photovoltaic module according to the string current of the photovoltaic module and the voltage of the direct current bus, and calculating the output power normalization value of the photovoltaic module at different time by using environmental parameters; calculating the annual attenuation rate of the actual output power of the photovoltaic module; and then calculating parameters by using the output power normalization values at different times, and finally obtaining the service life of the photovoltaic module according to the calculated parameters. The method can accurately calculate the operation life of the photovoltaic module, has the advantages of less data volume requirement, economy, rapidness, suitability for various operation environments and the like, and the calculation result can be used for evaluating the aging degree of the photovoltaic module, so that the economy and the reliability of a photovoltaic system are improved.
Description
Technical Field
The invention relates to a photovoltaic module service life calculation method, in particular to a photovoltaic module service life calculation method based on environmental parameters, and belongs to the technical field of photovoltaic modules.
Background
The photovoltaic module generally operates outdoors, and a series of attenuation modes are caused by environmental factors such as temperature, humidity, ultraviolet radiation and the like, so that the output power of the photovoltaic module is attenuated. The photovoltaic modules are core components of a photovoltaic system, and the reliability and the service life of the photovoltaic modules are directly related to the power generation efficiency of a photovoltaic power station, the investment and recovery of the power station, the power generation cost and the like. It is therefore desirable to evaluate the lifetime of photovoltaic modules operating under different environmental conditions.
The document "Ismail Kaaya, sascha Lindig, karl-Anders Weiss, et al. Photoviral life time for echo model based on deletion mapping". Progress in Photovo1 tics: research applications, 2020, 28 (10): 978-992 ("photovoltaic life prediction model based on degradation mode," photovoltaic technology development: research and application ", 2020, volume 28, page 10-page 992) provides a method for calculating the life of a photovoltaic module based on different degradation modes, but the method requires actual measurement of power points of the photovoltaic module under different environmental conditions for a long time, and since the operating life of the photovoltaic module is generally dozens of years, actual measurement power points of the photovoltaic module having the whole life cycle in China are very few, and the life calculation is difficult to be performed only through the actual measurement power points of the power of the photovoltaic module.
The literature "Sun-Ho Park, jae-Hoon Kim. Application of gamma process models to estimate the lifetime of photovoltaic modules" solar energy ", 2017, 147:390-398. ("estimating the lifetime of a photovoltaic module using a gamma process model," Solar Energy, "2017, volume 147, pages 390-398) shows a method for calculating the lifetime of a photovoltaic module based on an accelerated aging test under a damp-heat condition, but the method requires an accelerated aging test for several months or more, and the environmental conditions under the accelerated aging test do not strictly correspond to the process of power decay under actual operating conditions, so the method predicts the lifetime inaccurately and is time-consuming and labor-consuming.
In summary, in the prior art of calculating the service life of the photovoltaic module, the following problems still exist:
1. the operation time limit of the photovoltaic module is generally dozens of years, and the current performance data set of the whole life cycle is short, so that the life can not be calculated only through the performance data;
2. if the service life of the photovoltaic module is calculated by combining the accelerated aging test, the accelerated aging test time is required to be carried out for several months, and the process of the output power attenuation of the photovoltaic module in the accelerated environment is not completely corresponding to the process of the output power attenuation in the actual operation process, so that the requirements of the calculation time and the calculation precision cannot be met.
Disclosure of Invention
The invention aims to solve the problems of long time consumption, inaccurate calculation result and the like of the service life calculation of a photovoltaic module in the prior art, and particularly, a determined relational expression can be established between the output power of the photovoltaic module and environmental factors only by the performance parameters and the environmental parameters of the photovoltaic module within one year, the service life of the photovoltaic module can be calculated through the relational expression, the service life of the photovoltaic module can be quickly and accurately calculated, the aging degree of the photovoltaic module can be evaluated, and the economy and the reliability of a photovoltaic system can be improved.
The purpose of the invention is realized by that the invention provides a method for calculating the service life of a photovoltaic component based on environmental parameters, wherein the photovoltaic component is one of photovoltaic arrays, and the photovoltaic array is composed of N 1 Each photovoltaic group is connected in series-parallel and connected into a combiner box, the combiner box is connected with a direct current bus, and each photovoltaic group string is composed of N 2 The photovoltaic modules with the same structure are connected in series, wherein N is 1 ，N 2 Are all positive integers; the environmental parameters comprise environmental temperature, environmental humidity, instantaneous irradiation and ultraviolet radiation integral value;
the photovoltaic module service life calculation method based on the environmental parameters comprises the following steps:
Step 1.1, setting the current year as the (Y + 1) th year of operation of the photovoltaic module, and acquiring environmental parameters of the photovoltaic module in the whole year of the Y th year of operation through an environmental monitor of a photovoltaic power station to obtain four groups of sampling data: n is 1 Ambient temperature sampling data, denoted as T i1 ，i1＝1，2，...，n 1 ；n 2 The sampling data of the environmental humidity is recorded as RH i2 ，i2＝1，2，...，n 2 ；n 3 Instantaneous irradiation sampling data, denoted as S i3 ，i3＝1，2，...，n 3 ；n 4 Sampling data of integral value of ultraviolet radiation, recorded as UV i4 ，i4＝1，2，...，n 4 ；
Calculating n 1 Ambient temperature sampling data T i1 Average value of (2)
And standard deviation σ 1 ，n 2 Ambient humidity sampling data RH i2 Average value of (2)
And standard deviation σ 2 ，n 3 Instantaneous irradiation sampling data S i3 Average value of (2)
And standard deviation σ 3 ，n 4 Sampling data UV of ultraviolet radiation integral value i4 Average value of (2)
And standard deviation σ 4 The calculation formulas are respectively as follows:
step 1.2, four groups of samples obtained in step 1.1 are sampledAnd eliminating abnormal values exceeding the upper and lower triple standard deviation ranges of the respective average values in the data to obtain four groups of environment parameter sampling data after eliminating the abnormal values: m is 1 Ambient temperature T 'after elimination of abnormal value' j1 ，j1＝1，2，...，m 1 ；m 2 Environmental humidity RH 'after elimination of abnormal values' j2 ，j2＝1，2，...，m 2 ；m 3 Instantaneous irradiation S 'after elimination of abnormal value' j3 ，j3＝1，2，...，m 3 ；m 4 Ultraviolet radiation integral value UV 'after eliminating abnormal value' j4 ，j4＝1，2，...，m 4 (ii) a Wherein, the four groups of data without abnormal values all comprise 1 month and 1 day sampling data and 12 months and 31 days sampling data;
from m 1 Ambient temperature T 'after elimination of abnormal value' j1 The highest ambient temperature of each day is extracted to form a group of data and is recorded as the highest ambient temperature T aj5 J5=1, 2., 365; from m 1 Ambient temperature T 'after eliminating abnormal value' j1 The lowest ambient temperature of each day is extracted to form a group of data and is recorded as the lowest ambient temperature T αj6 J6=1,2,. 365; from m 2 Environmental humidity RH 'after elimination of abnormal values' j2 Extracting the highest ambient humidity of each day to form a group of data, and recording as the highest ambient humidity RH βj7 J7=1,2,. 365; from m 2 Environmental humidity RH 'after elimination of abnormal values' j2 The lowest ambient humidity of each day is extracted to form a group of data and is recorded as the lowest ambient humidity RH βj8 ，j8＝1，2，...，365；
Calculating the annual average value T of the ambient temperature P Annual average value of ambient maximum temperature T U Annual average value of lowest ambient temperature T L Annual average ambient humidity RH P Annual average RH of maximum ambient humidity U And the annual average RH of the minimum humidity of the environment L Integral annual average value UV of ultraviolet radiation P The calculation formula is as follows:
step 1.3, at m 1 Ambient temperature T 'after eliminating abnormal value' j1 And m 3 Instantaneous irradiation S 'after elimination of abnormal value' j3 Take 1 month and 1 day t a The time data are recorded as the first environmental temperature T 01 And a first instant irradiation S 01 ，t a At any time within 1 month and 1 day, but the first environment temperature T at the time is ensured 01 And a first instant irradiation S 01 All exist at the same time; at m 1 Ambient temperature T 'after elimination of abnormal value' j1 And m 3 Instantaneous irradiation S 'after elimination of abnormal value' j3 Taking 12 months and 31 days t b The time data are recorded as the second environmental temperature T 02 And a second instantaneous irradiation S 02 ，t b At any time within 31 days of 12 months, but the second ambient temperature T is ensured at the time 02 And a second instantaneous irradiation S 02 Are all at the same time(ii) present;
step 1.4, acquiring t through a combiner box a Total current I of photovoltaic string at any moment 1 And DC bus voltage U 1 Obtaining t through the combiner box b Total current I of photovoltaic string at any moment 2 And DC bus voltage U 2 ；
Note t a At the moment, the current of the photovoltaic module is I' 1 ，
Note t a The voltage of the photovoltaic module is U 'at any moment' 1 ，
Note t b The current of the photovoltaic module is I at the moment' 2 ，
Note t b The voltage of the photovoltaic module is U' 2 ，
Obtaining an open-circuit voltage temperature coefficient omega, a short-circuit current temperature coefficient xi and a maximum power point power P of the photovoltaic module by inquiring the product specification of the photovoltaic module m ；
Calculating the 1 st, 1 st of the Yth year under the standard condition state of the photovoltaic module a Maximum power point voltage U at time m1 Maximum power point current I m1 And the normalized value P of the output power 1 (ii) a Calculating the 12 th month, 31 th day and t day of the Y year under the standard condition state of the photovoltaic module b Maximum power point voltage U at time m2 Maximum power point current I m2 And the normalized value P of the output power 2 The calculation formula is respectively:
step 3, calculating the service life of the photovoltaic module
Step 3.1, note k 1 Annual decay Rate, k, of output power decay of a photovoltaic Module caused by Damp and Heat modes 2 Annual decay Rate, k, for the output power decay of a photovoltaic Module caused by the ultraviolet radiation mode 3 Annual decay rate, k, for the temperature cycling mode causing the photovoltaic module to experience an output power decay 4 For the annual attenuation rate of the output power attenuation of the photovoltaic module caused by the humidity circulation mode, the calculation formula is respectively as follows:
annual decay rate k for causing output power decay of photovoltaic module caused by damp and hot mode 1 Annual attenuation rate k of photovoltaic module output power attenuation caused by ultraviolet radiation mode 2 Annual attenuation rate k of photovoltaic module with output power attenuation caused by temperature circulation mode 3 Annual attenuation rate k of photovoltaic module output power attenuation caused by humidity circulation mode 4 Obtaining the annual attenuation rate k of the actual output power of the photovoltaic module after coupling 0 The calculation formula is as follows:
k 0 ＝(1+k 1 ) 0.3 (1+k 2 ) 0.49 (1+k 3 ) 0.42 (1+k 4 ) 0.28 -1
wherein k is 1 ，k 2 ，k 3 ，k 4 ，k 0 The units of (A) are all%/year;
step 3.2, giving a normalized value P corresponding to the output power 1 And the normalized value P of the output power 2 The power attenuation formula of (1):
in the formula, B 0 Is the output power sensitivity coefficient mu of the photovoltaic module 0 The shape parameter of the output power attenuation curve of the photovoltaic module is obtained;
calculating the operating life t of the photovoltaic module 0 The calculation formula is as follows:
compared with the prior art, the invention has the beneficial effects that:
1. the operation life of the photovoltaic module can be accurately calculated, and a long-term accelerated aging test is not needed;
2. the method can be applied to photovoltaic power stations with short operation time, and operation data of the whole life cycle are not needed;
3. the service life of the photovoltaic module under any environmental condition can be calculated, and the method is economical and rapid;
4. the calculation result can be used for evaluating the aging degree of the photovoltaic assembly, and the economy and the reliability of the photovoltaic system are improved.
Drawings
Fig. 1 is a schematic connection diagram of a photovoltaic module according to the present invention.
Fig. 2 is a block diagram of a method for calculating the service life of a photovoltaic module according to the present invention.
Fig. 3 is a general flowchart of the photovoltaic module life calculating method according to the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
Fig. 1 is a schematic connection diagram of a photovoltaic module according to the present invention. As can be seen from FIG. 1, the photovoltaic module is one of the photovoltaic arrays consisting of N 1 Each photovoltaic group is connected in series-parallel and connected into a combiner box, the combiner box is connected with a direct current bus, and each photovoltaic group string is composed of N photovoltaic groups 2 The photovoltaic modules with the same structure are connected in series, wherein N is 1 ，N 2 Are all positive integers.
In this embodiment, N 1 ＝16，N 2 ＝22。
Fig. 2 is a block diagram of a method for calculating the lifetime of a photovoltaic module according to the present invention, and fig. 3 is a general flowchart of the method for calculating the lifetime of a photovoltaic module according to the present invention. As can be seen from fig. 2 and 3, the environmental parameters include an ambient temperature, an ambient humidity, an instantaneous irradiance and an ultraviolet radiation integral value. The photovoltaic module service life calculation method based on the environmental parameters comprises the following steps:
Step 1.1, setting the current year as the (Y + 1) th year of the operation of the photovoltaic module, passing through a ring of a photovoltaic power stationThe environment monitor acquires environmental parameters of a whole year in the Y-th year of operation of the photovoltaic module, and four groups of sampling data are obtained: n is 1 Sampling data of the ambient temperature, denoted as T i1 ，i1＝1，2，...，n 1 ；n 2 The sampling data of the environmental humidity is recorded as RH i2 ，i2＝1，2，...，n 2 ；n 3 Instantaneous irradiation sampling data, denoted as S i3 ，i3＝1，2，...，n 3 ；n 4 Sampling data of integral value of ultraviolet radiation, recorded as UV i4 ，i4＝1，2，...，n 4 。
Calculating n 1 Ambient temperature sampling data T i1 Average value of (2)
And standard deviation σ 1 ，n 2 Ambient humidity sampling data RH i2 Average value of (2)
Sum standard deviation σ 2 ，n 3 Instantaneous irradiation sampling data S i3 Average value of (2)
Sum standard deviation σ 3 ，n 4 Ultraviolet radiation integral value sampling data UV i4 Average value of (2)
And standard deviation σ 4 The calculation formulas are respectively as follows:
step 1.2, eliminating abnormal values in the standard deviation ranges exceeding the upper and lower triples of the respective average values in the four groups of sampling data obtained in step 1.1, and obtaining four groups of environmental parameter sampling data after eliminating the abnormal values: m is a unit of 1 Ambient temperature T 'after elimination of abnormal value' j1 ，j1＝1，2，...，m 1 ；m 2 Environmental humidity RH 'after elimination of abnormal values' j2 ，j2＝1，2，...，m 2 ；m 3 Instantaneous irradiation S 'after elimination of abnormal value' j3 ，j3＝1，2，...，m 3 ；m 4 Ultraviolet radiation integral value UV 'after eliminating abnormal value' j4 ，j4＝1，2，...，m 4 (ii) a The four groups of data without abnormal values all comprise 1 month and 1 day sampling data and 12 months and 31 days sampling data.
Specifically, in this embodiment, the environmental temperature sampling data of one year is obtained by the environmental detector, and the annual average value of the environmental temperatures is calculated to be 18.36 ℃, and the standard deviation is 8.35 ℃, so that the environmental temperatures outside the range from-6.69 ℃ to 43.41 ℃ in the sampled environmental temperatures are excluded.
From m 1 Ambient temperature T 'after elimination of abnormal value' j1 The highest ambient temperature of each day is extracted to form a group of data and is recorded as the highest ambient temperature T αj5 J5=1, 2., 365; from m 1 Ambient temperature T 'after elimination of abnormal value' j1 The lowest ambient temperature of each day is extracted to form a group of data and is recorded as the lowest ambient temperature T αj6 J6=1,2,. 365; from m 2 Environmental humidity RH 'after elimination of abnormal values' j2 The highest ambient humidity of each day is extracted to form a group of data and is recorded as the highest ambient humidity RH βj7 J7=1, 2., 365; from m 2 The environmental humidity RH' j after eliminating abnormal values 2 Extracting the lowest ambient humidity of each day to form a group of data, and recording the data as the lowest ambient humidity RH βj8 ，j8＝1，2，...，365。
Calculating the annual average value T of the ambient temperature P Annual average value of ambient maximum temperature T U Annual average value of lowest ambient temperature T L And annual average value of ambient humidity RH P Annual average RH of maximum ambient humidity U And the annual average value RH of the minimum ambient humidity L Integral annual average value UV of ultraviolet radiation P The calculation formula is as follows:
specifically, in this embodiment, T P ＝289.32Kelvin，T U ＝297.15Kelvin，T L ＝282.09Kelvin，RH P ＝70.11％，RH U ＝92.30％，RH L ＝36.30％，UV P ＝4.2Kwh/m 2 。
Step 1.3, at m 1 Ambient temperature T 'after elimination of abnormal value' j1 And m 3 Instantaneous irradiation S 'after elimination of abnormal value' j3 Take 1 month and 1 day t a The time data are recorded as the first environmental temperature T 01 And a first instant irradiation S 01 ，t a At any time within 1 month and 1 day, but the first environment temperature T at that time should be ensured 01 And a first instant irradiation S 01 All exist at the same time; at m 1 Ambient temperature T 'after elimination of abnormal value' j1 And m 3 Instantaneous irradiation S 'after elimination of abnormal value' j3 Taking 12 months and 31 days t b The data of the time are recorded as the second ambient temperature T 02 And a second instantaneous irradiation S 02 ，t b At any time within 31 days of 12 months, but the second ambient temperature T is ensured at the time 02 And a second instantaneous irradiation S 02 All exist simultaneously.
Specifically, in this embodiment, t a 1/month 1, 11 hours, 0 minutes, T in 2021 01 ＝17.5℃，S 01 ＝620W/m 2 ；t b 0 min at 12, 31 and 15 days in 2021, T 02 ＝12.0℃，S 01 ＝320W/m 2 。
Step 1.4, obtaining t through a combiner box a Total current I of photovoltaic string at any moment 1 And DC bus voltage U 1 Obtaining t through the combiner box b Total current I of photovoltaic string at any moment 2 And DC bus voltage U 2 。
Note t a At the moment, the current of the photovoltaic module is I' 1 ，
Note t a The voltage of the photovoltaic module is U 'at any moment' 1 ，
Note t b At the moment, the current of the photovoltaic module is I' 2 ，
Note t b The voltage of the photovoltaic module is U 'at any moment' 2 ，
In this example, I 1 ＝60.784A，U 1 ＝617.782V，I 2 ＝24.16A，U 2 ＝585.288V，I′ 1 ＝3.799A，U′ 1 ＝28.081V，I′ 2 ＝1.510A，U′ 2 ＝26.604V。
Obtaining an open-circuit voltage temperature coefficient omega, a short-circuit current temperature coefficient xi and a maximum power point power P of the photovoltaic module by inquiring the product specification of the photovoltaic module m ；
Calculating the 1 st, 1 st of the Yth year under the standard condition state of the photovoltaic module a Maximum power point voltage U at time m1 Maximum power point current I m1 And the normalized value P of the output power 1 (ii) a Calculating the 12 th month, 31 th day and t day of the Y year under the standard condition state of the photovoltaic module b Maximum power point voltage U at time m2 Maximum power point current I m2 And the normalized value P of the output power 2 The calculation formula is respectively:
in this embodiment, ω =0.330%/° c, ξ =0.058%/° c, P m =260W, the standard condition state of the photovoltaic module is 1000W/m 2 At 25 ℃. I is m1 ＝7.860A，U m1 ＝29.543V，P1＝0.8931，I m2 ＝7.792A，U m2 ＝29.444V，P 2 ＝0.8824。
Step 3, calculating the service life of the photovoltaic module
Step 3.1, note k 1 Annual decay Rate, k, for Damp and Hot modes causing output Power decay of photovoltaic modules 2 Annual decay rate, k, for the output power decay of a photovoltaic module caused by the ultraviolet radiation pattern 3 Annual decay rate, k, for the temperature cycling mode causing the photovoltaic module to experience an output power decay 4 For the annual attenuation rate of the output power attenuation of the photovoltaic module caused by the humidity circulation mode, the calculation formula is respectively as follows:
annual decay rate k for causing output power decay of photovoltaic module by damp and hot mode 1 Annual attenuation rate k of photovoltaic module output power attenuation caused by ultraviolet radiation mode 2 Annual attenuation rate k of photovoltaic module output power attenuation caused by temperature circulation mode 3 Annual attenuation rate k of photovoltaic module output power attenuation caused by humidity circulation mode 4 Obtaining the annual attenuation rate k of the actual output power of the photovoltaic module after coupling 0 The calculation formula is as follows:
k 0 ＝(1+k 1 ) 0.3 (1+k 2 ) 0.49 (1+k 3 ) 0.42 (1+k 4 ) 0.28 -1
wherein k is 1 ，k 2 ，k 3 ，k 4 And k 0 The units of (A) are all%/year.
In this embodiment, k 1 K of 1.6497%/year 2 K of 1.4795%/year 3 K of 1.1145%/year 4 K of 1.4060%/year 0 = 2.6606%/year.
Step 3.2, giving a normalized value P corresponding to the output power 1 And the normalized value P of the output power 2 The power attenuation formula of (1):
in the formula, B 0 Is the output power sensitivity coefficient mu of the photovoltaic module 0 Outputting shape parameters of a power attenuation curve for the photovoltaic module;
calculating the operating life t of the photovoltaic module 0 The calculation formula is as follows:
in this embodiment, B 0 ＝480，μ 0 ＝0.218，t 0 =20.33 years.
Claims (1)
1. A photovoltaic module service life calculation method based on environmental parameters is characterized in that the photovoltaic module is one of photovoltaic arrays, and the number of the photovoltaic arrays is N 1 Each photovoltaic group is connected in series-parallel and connected into a combiner box, the combiner box is connected with a direct current bus, and each photovoltaic group string is composed of N photovoltaic groups 2 The photovoltaic modules with the same structure are connected in series, wherein N 1 ，N 2 Are all positive integers; the environmental parameters comprise an environmental temperature, an environmental humidity, an instantaneous irradiation and an ultraviolet radiation integral value;
the photovoltaic module service life calculation method based on the environmental parameters is characterized by comprising the following steps of:
step 1, preprocessing environmental parameters and photovoltaic module performance parameters
Step 1.1, setting the current year as the (Y + 1) th year of the operation of the photovoltaic module, and acquiring environmental parameters of the photovoltaic module in the whole year in the Y-th year of the operation through an environmental monitor of a photovoltaic power station to obtain four groups of sampling data: n is a radical of an alkyl radical 1 Ambient temperature sampling data, denoted as T i1 ，i1＝1，2，...，n 1 ；n 2 The sampling data of the environmental humidity is recorded as RH i2 ，i2＝1，2，...，n 2 ；n 3 Instantaneous irradiation sampling data, denoted as S i3 ，i3＝1，2，...，n 3 ；n 4 Integral value of ultraviolet radiationSample data, denoted UV i4 ，i4＝1，2，...，n 4 ；
Calculating n 1 Ambient temperature sampling data T i1 Average value of (2)
Sum standard deviation σ 1 ，n 2 Ambient humidity sampling data RH i2 Average value of (2)
And standard deviation σ 2 ，n 3 Instantaneous irradiation sampling data S i3 Average value of (2)
And standard deviation σ 3 ，n 4 Ultraviolet radiation integral value sampling data UV i4 Average value of (2)
Sum standard deviation σ 4 The calculation formulas are respectively as follows:
step 1.2, eliminating abnormal values exceeding the upper and lower triple standard deviation ranges of the respective average values in the four groups of sampling data obtained in step 1.1 to obtain four groups of environmental parameter sampling data after eliminating the abnormal values: m is a unit of 1 Ambient temperature T 'after elimination of abnormal value' j1 ，j1＝1，2，...，m 1 ；m 2 Ambient humidity RH 'after elimination of abnormal values' j2 ，j2＝1，2，...，m 2 ；m 3 Instantaneous irradiation S 'after elimination of abnormal value' j3 ，j3＝1，2，...，m 3 ；m 4 Ultraviolet radiation integral value UV 'after eliminating abnormal value' j4 ，j4＝1，2，...，m 4 (ii) a Wherein, the four groups of data without abnormal values all comprise 1 month and 1 day sampling data and 12 months and 31 days sampling data;
from m 1 Ambient temperature T 'after elimination of abnormal value' j1 The highest ambient temperature of each day is extracted to form a group of data and is recorded as the highest ambient temperature T αj5 J5=1,2,. 365; from m 1 Ambient temperature T 'after elimination of abnormal value' j1 The lowest ambient temperature of each day is extracted to form a group of data and is recorded as the lowest ambient temperature T αj6 J6=1, 2.., 365; from m 2 Environmental humidity RH 'after elimination of abnormal values' j2 The highest ambient humidity of each day is extracted to form a group of data and is recorded as the highest environmentHumidity RH βj7 J7=1, 2., 365; from m 2 Environmental humidity RH 'after elimination of abnormal values' j2 Extracting the lowest ambient humidity of each day to form a group of data, and recording the data as the lowest ambient humidity RH βj8 ，j8＝1，2，...，365；
Calculating the annual average value T of the ambient temperature P Annual average value of ambient maximum temperature T U Annual average value of ambient minimum temperature T L And annual average value of ambient humidity RH P Annual average RH of maximum ambient humidity U And the annual average RH of the minimum humidity of the environment L Integral annual average value UV of ultraviolet radiation P The calculation formula is as follows:
step 1.3, at m 1 Ambient temperature T 'after eliminating abnormal value' j1 And m 3 Instantaneous irradiation S 'after elimination of abnormal value' j3 Take 1 month and 1 day t a The time data are recorded as the first environmental temperature T 01 And a first instant irradiation S 01 ，t a At any time within 1 month and 1 day, but the first environment temperature T at that time should be ensured 01 And a first instant irradiation S 01 All exist at the same time; at m 1 Ambient temperature T 'after elimination of abnormal value' j1 And m 3 Instantaneous irradiation S 'after elimination of abnormal value' j3 Taking 12 months and 31 days t b The data of the time are recorded as the second ambient temperature T 02 And a second instantaneous irradiation S 02 ，t b At any time within 31 days of 12 months, but the second ambient temperature T is ensured at the time 02 And a second instantaneous irradiation S 02 All exist at the same time;
step 1.4, obtaining t through a combiner box a Total current I of photovoltaic string at any moment 1 And DC bus voltage U 1 Obtaining t through the combiner box b Total current I of photovoltaic string at any moment 2 And DC bus voltage U 2 ；
Note t a The current of the photovoltaic module is I at the moment' 1 ，
Note t a The voltage of the photovoltaic module is U 'at any moment' 1 ，
Note t b At the moment, the current of the photovoltaic module is I' 2 ，
Note t b The voltage of the photovoltaic module is U' 2 ，
Step 2, normalizing the output power of the photovoltaic module
Obtaining an open-circuit voltage temperature coefficient omega, a short-circuit current temperature coefficient xi and a maximum power point power P of the photovoltaic module by inquiring the product specification of the photovoltaic module m ；
Calculating the 1 st, the Yth year and the 1 st day of the photovoltaic module under the standard condition state a Maximum power point voltage U at time m1 Maximum power point current I m1 And the normalized value P of the output power 1 (ii) a Calculating the 12 th month, 31 th day and t day of the Y year under the standard condition state of the photovoltaic module b Maximum power point voltage U at time m2 Maximum power point current I m2 And the normalized value P of the output power 2 The calculation formula is respectively:
step 3, calculating the service life of the photovoltaic module
Step 3.1, note k 1 Annual decay Rate, k, for Damp and Hot modes causing output Power decay of photovoltaic modules 2 Annual decay rate, k, for the output power decay of a photovoltaic module caused by the ultraviolet radiation pattern 3 Annual decay rate, k, for the temperature cycling mode causing the photovoltaic module to experience an output power decay 4 For the annual attenuation rate of the output power attenuation of the photovoltaic module caused by the humidity circulation mode, the calculation formula is respectively as follows:
annual decay rate k for causing output power decay of photovoltaic module caused by damp and hot mode 1 Annual attenuation rate k of photovoltaic module output power attenuation caused by ultraviolet radiation mode 2 Annual attenuation rate k of photovoltaic module output power attenuation caused by temperature circulation mode 3 Annual attenuation rate k of photovoltaic module output power attenuation caused by humidity circulation mode 4 Obtaining the annual attenuation rate k of the actual output power of the photovoltaic module after coupling 0 The calculation formula is as follows:
k 0 ＝(1+k 1 ) 0.3 (1+k 2 ) 0.49 (1+k 3 ) 0.42 (1+k 4 ) 0.28 -1
wherein k is 1 ，k 2 ，k 3 ，k 4 ，k 0 Are all in% >, andyear;
step 3.2, giving a normalized value P corresponding to the output power 1 And the normalized value P of the output power 2 The power attenuation formula (2):
in the formula, B 0 Is the output power sensitivity coefficient mu of the photovoltaic module 0 The shape parameter of the output power attenuation curve of the photovoltaic module is obtained;
calculating the operating life t of the photovoltaic module 0 The calculation formula is as follows:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210363024.8A CN114692064B (en) | 2022-04-07 | 2022-04-07 | Photovoltaic module service life calculation method based on environmental parameters |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210363024.8A CN114692064B (en) | 2022-04-07 | 2022-04-07 | Photovoltaic module service life calculation method based on environmental parameters |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114692064A CN114692064A (en) | 2022-07-01 |
| CN114692064B true CN114692064B (en) | 2022-11-29 |
Family
ID=82142460
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210363024.8A Active CN114692064B (en) | 2022-04-07 | 2022-04-07 | Photovoltaic module service life calculation method based on environmental parameters |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114692064B (en) |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160162616A1 (en) * | 2014-03-27 | 2016-06-09 | King Fahd University Of Petroleum And Minerals | Performance and life prediction model for photovoltaic module: effect of encapsulant constitutive behavior |
| CN105610401B (en) * | 2016-03-18 | 2018-09-11 | 天合光能股份有限公司 | Simulated photoinduced attenuation device and method for solar crystalline silicon photovoltaic cell |
| CN106204332B (en) * | 2016-07-19 | 2019-11-08 | 河海大学常州校区 | A kind of photovoltaic plant efficiency decaying prediction technique |
| CN110414026B (en) * | 2018-04-28 | 2023-08-22 | 东君新能源有限公司 | Method for evaluating service life of photovoltaic module |
| CN108663486A (en) * | 2018-06-11 | 2018-10-16 | 河海大学常州校区 | A kind of determination of the weather-proof life factor in photovoltaic module different regions and verification method |
| EP3624332B1 (en) * | 2018-09-17 | 2021-05-12 | IMEC vzw | Operational lifetime prediction of photovoltaic devices |
| CN109885856B (en) * | 2018-12-20 | 2022-05-31 | 合肥工业大学 | Photovoltaic module hot spot power loss calculation method based on I-V curve |
| CN112019158A (en) * | 2019-05-28 | 2020-12-01 | 汉能移动能源控股集团有限公司 | Outdoor life evaluation method and device for photovoltaic module |
-
2022
- 2022-04-07 CN CN202210363024.8A patent/CN114692064B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN114692064A (en) | 2022-07-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109873610B (en) | Photovoltaic array fault diagnosis method based on IV characteristic and depth residual error network | |
| CN107229824B (en) | Photovoltaic power station power generation unit power curve modeling method and device | |
| Whitaker et al. | Application and validation of a new PV performance characterization method | |
| CN107358335B (en) | Distributed photovoltaic efficiency evaluation method and system based on Internet | |
| CN116054740A (en) | Fault detection and fault detection model processing method and device for photovoltaic power generation equipment | |
| CN114139614B (en) | Fisher photovoltaic module hot spot diagnosis method and system based on typical correlation analysis feature extraction | |
| CN106597151B (en) | Photovoltaic power station switched inverter tracking efficiency testing method and system | |
| CN114692064B (en) | Photovoltaic module service life calculation method based on environmental parameters | |
| CN110336322B (en) | Photovoltaic power generation access capacity determination method based on daily minimum load confidence interval | |
| CN109885856B (en) | Photovoltaic module hot spot power loss calculation method based on I-V curve | |
| CN111799843A (en) | Multi-dimensional all-state data modeling method for probability fusion type photovoltaic power station | |
| Walters et al. | PV module loss analysis using system in-situ monitoring data | |
| CN116879672A (en) | Intelligent identification method for total fault of area based on time convolution network | |
| CN110796292A (en) | Photovoltaic power short-term prediction method considering haze influence | |
| CN114491389A (en) | Method for extracting and estimating equivalent circuit parameters of solar cell module | |
| CN116015209A (en) | Photovoltaic power station component performance management method, system, device and storage medium | |
| Almomani et al. | Application of artificial intelligence techniques for modeling and simulation of photovoltaic arrays | |
| Rodrigues et al. | Theoretical–experimental evaluation of partially shaded solar photovoltaic arrays through Methodological Framework: a Case Study Involving two identical 1.5 kWp photovoltaic generators | |
| Lodhi et al. | Modelling and experimental characteristics of photovoltaic modules in typical days at an actual photovoltaic power station | |
| JP6393732B2 (en) | Solar power generation monitoring device and solar power generation monitoring method | |
| CN113435038A (en) | Photovoltaic power generation system loss online analysis system and method | |
| Aboelmagd et al. | Failure Analysis in Photovoltaic Power Systems Using an Artificial Neural Network | |
| CN111027816A (en) | Photovoltaic power generation efficiency calculation method based on data envelope analysis | |
| CN111881565B (en) | Photovoltaic module process data reconstruction method | |
| CN117396733B (en) | System, computing device, and method for estimating solar irradiance falling on a photovoltaic string |
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 |