CN109409420A - A kind of photovoltaic group string method for diagnosing faults under non-uniform irradiation degree - Google Patents
A kind of photovoltaic group string method for diagnosing faults under non-uniform irradiation degree Download PDFInfo
- Publication number
- CN109409420A CN109409420A CN201811168547.7A CN201811168547A CN109409420A CN 109409420 A CN109409420 A CN 109409420A CN 201811168547 A CN201811168547 A CN 201811168547A CN 109409420 A CN109409420 A CN 109409420A
- Authority
- CN
- China
- Prior art keywords
- photovoltaic
- spark
- group string
- fireworks
- component
- 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
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/2433—Single-class perspective, e.g. one-against-all classification; Novelty detection; Outlier detection
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Photovoltaic Devices (AREA)
- Testing Of Individual Semiconductor Devices (AREA)
Abstract
A kind of photovoltaic group string method for diagnosing faults under non-uniform irradiation degree of the present invention, in the case where unsettled causes that irradiation level is inconsistent in photovoltaic group string, the photovoltaic module environmental parameter of the irradiance sensor and temperature sensor measurement different location laid first by multiple spot carries out theoretical calculation to photovoltaic electric current, voltage and output power using photovoltaic cell equivalent-circuit model.Meanwhile at the header box of photovoltaic group string, actual photovoltaic electric current, voltage and output power are measured.Then open-circuit failure and short trouble carry out malfunction coding, obtain the theoretical output power function about malfunction coding.Finally using the absolute value of the difference of photovoltaic group string theory output power and real output as optimization aim, objective function is optimized using the reinforcing fireworks algorithm with Fast Algorithm for Mining and prospecting performance, realizes the fault diagnosis and location of photovoltaic group string.The present invention enables to the staff of photovoltaic plant operation and maintenance to overhaul in time to the photovoltaic module arrangement broken down.
Description
Technical field
The present invention relates to a kind of photovoltaic group string method for diagnosing faults, and in particular to a kind of photovoltaic under non-uniform irradiation degree
Group string method for diagnosing faults.
Background technique
As the increasingly depleted and problem of environmental pollution of the traditional energies such as fossil fuel is more and more prominent, the mankind are to new
The demand of the energy is more next urgent.Solar energy has many advantages, such as clean and environmental protection, safe and reliable, sustainable development, is more satisfactory
Renewable new energy.Currently, photovoltaic power generation is a kind of technology effectively converted solar energy into electrical energy, can effectively alleviate
Energy crisis and improvement problem of environmental pollution.Important component of the photovoltaic array as photovoltaic generating system, theoretically it is flat
The equal service life, photovoltaic module can encounter such as packaging technology, use environment evil up to 20~30 years, but in actual engineer application
The problems such as bad, hot spot phenomenon, mechanical damage, fault (open circuit and short circuit), so that photovoltaic module normal service life is significantly
Shorten.Wherein, short circuit (SC) failure of photovoltaic group string and open circuit (OC) failure influence the normal operation of photovoltaic power generation serious, very
To the safety for threatening entire photovoltaic plant.Therefore, it monitors the operating status of photovoltaic array in real time and finds failure, positioning in time
Failure be very it is necessary to.
In traditional photovoltaic array fault detection approach, there are offline inspection and on-line checking.Offline inspection is due to needing
It could implement after photovoltaic generating system is integrally stopped transport, but photovoltaic plant needs persistently to put into operation to ensure the stability of power supply
And continuity.Therefore, offline fault detection seldom large-scale use in practice.On-line checking generally has based on mathematical model
Photovoltaic array fault detection method and diagnosing failure of photovoltaic array system based on intelligence, not only mathematical model is multiple for these methods
It is miscellaneous, computationally intensive, real-time performance is poor, and under changeable such environmental effects, there are fault diagnosis accuracy rate is low and failure
The high disadvantage of rate of false alarm to aggravate the working strength of live operation maintenance personnel, or even occurs serious failure and fails to report situation, directly
Connect the benefit and security and stability for affecting photovoltaic generating system.
How engineering in practice, is efficiently in time analyzed the working condition of photovoltaic group string, quickly determines photovoltaic event
Hinder type and position is a urgent problem to be solved, and the important content of photovoltaic plant operation and maintenance from now on.
Summary of the invention
In order to overcome the shortcomings of above-mentioned existing fault diagnosis technology, the purpose of the present invention is to provide one kind based on reinforcing cigarette
The photovoltaic group string method for diagnosing faults under non-uniform irradiation degree of flower algorithm.
The present invention adopts the following technical scheme that realize:
A kind of photovoltaic group string method for diagnosing faults under non-uniform irradiation degree, comprising the following steps:
Step 1: for the photovoltaic group string of shared N block assembly, acquiring the output power P of photovoltaic group stringMeasured, voltage
VMeasuredWith electric current IMeasuredAnd in photovoltaic group string component i irradiation level GiWith component temperature Ti;
Step 2: according to photovoltaic equivalent-circuit model and the irradiation level G of acquisitioniWith component temperature TiPhotovoltaic group is calculated
The theoretical maximum output power of component i in string;
Step 3: calculating photovoltaic group string theory output power PArray;I.e. photovoltaic group string output power is equal to light in photovoltaic group string
Lie prostrate the summation of component output power;
Step 4: determining fault threshold ε, and judge whether photovoltaic group string has an exception;If PMeasured+ ε < PArray, then
Illustrate that a group string is operating abnormally, carries out step 5;Otherwise illustrate that group string is without exception;
Step 5: establishing fault equivalence model;For the short trouble and open-circuit fault in photovoltaic generating system, both
It is approximately a kind of switch models, short trouble is equivalent to photovoltaic panel parallel diode when occurring and is short-circuited, that is, is equivalent to photovoltaic panel
The switch K1 of a closure in parallel;It is then equivalent to the switch K2 of photovoltaic output one disconnection of series connection when open circuit, has cut off photovoltaic panel
External output channel;And normal operation is then equivalent to switch K1 disconnection, switch K2 closure;
Step 6: building fitness function, for characterizing the fitness size of each fireworks, the value of the fitness is determined
The spark quantity and burst radius that fireworks give out;
Fitness=f (NSC,NOC, SC, OC)=| PMeasured-PArray|
It is using the target for strengthening fireworks algorithm optimization then:
minf(NSC,NOC, SC, OC)=| PMeasured-PArray|
Wherein NSCFor short trouble number, NOCFor open-circuit fault number;
Step 7: initialization fireworks population randomly chooses a fireworks of M ' as preliminary iteration fireworks kind in feasible domain space
Group, and initialize maximum number of iterations iter;
Step 8: calculating the fitness size Fitness of fireworks i in fireworks populationi;Then it is normalized, obtainsIt will retain in populationM group carry out it is next
Step operation;
Step 9: design explosion operator, is (N according to each positionSC,NOC,SC,OC)iThe fitness of fireworks individual is big
It is small, the explosion initial radium R=ρ Fitness of the fireworks is calculated, wherein ρ is burst radius coefficient of dilatation, and the coefficient is by feasible
The specific size range in domain determines, formula show fitness it is smaller be more preferably fireworks, then its explosive range will reduce, progress
Search refinement;And fitness it is big fireworks burst radius it is big, carry out range exploitation;Similarly according to fireworks individual adaptation degree size,
Calculate explosive spark number Wherein θ is explosive spark number coefficient of dilatation, and the coefficient is by fireworks kind
The size of group determines, while adding the normal random number of 0 mean value again;To which the radius of each spark is Ri=R+N (0,
σR), increasing random number is the diversity in order to improve spark, while generating the deflection θ for being uniformly distributed generationi=U (0,
2);This is arrived, the position of spark i is obtained are as follows:
(NSC,NOC,SC,OC)′i=(NSC,NOC,SC,OC)i
+ (NSC,NOC,(Ricosθi)bin,(Risinθi)bin)i
Wherein, (Ricosθi)binIt indicates real number RicosθiSwitch to N binary forms, and is extracted as the position group of " 1 "
At set SC;Similarly (Risinθi)binBy real number RisinθiSwitch to N binary forms, and is extracted as the position number composition collection of " 1 "
Close OC;
Step 10: design Gaussian mutation operator carries out the spark that fireworks are broken into high to increase spark longitudinal direction diversity
This mutation operation selects n ' < n spark to carry out Gaussian mutation as variation population, that is,
Ri=N (R, σR);Wherein here i-th spark variation population in;To which the spark position of variation is
(NSC,NOC,SC,OC)′i=N (R, σR)×(NSC,NOC,SC,OC)i;
Step 11: design spark Collision operator strengthens the diversity of spark in fireworks onset process, i.e. fireworks exploded
Cheng Zhong, different spark collisions will generate new spark;That is former generation's spark (NSC,NOC,SC,OC)1(NSC,NOC,SC,OC)2It touches
It hits and generates offspring's spark (NSC,NOC, SC, OC) '1(NSC,NOC,SC,OC)′2;Be randomly generated first between one 0 to 1 with
Machine number n recycles multinomial probability distribution to generate two random number γ1And γ2,
Wherein
ηcIt is a nonnegative number;
To obtain offspring's spark are as follows:
(NSC,NOC,SC,OC)p1
=0.5 [((NSC,NOC,SC,OC)2+(NSC,NOC,SC,OC)1)
-γ1(|(NSC,NOC,SC,OC)2-(NSC,NOC,SC,OC)1|)]
(NSC,NOC,SC,OC)p2
=0.5 [((NSC,NOC,SC,OC)2+(NSC,NOC,SC,OC)1)
-γ2(|(NSC,NOC,SC,OC)2-(NSC,NOC,SC,OC)1|)]
Step 12: selecting fireworks of the next-generation spark as next step iteration;
Step 13: judging termination condition;If f ((NSC,NOC,SC,OC)k)≤ε obtains VArrayAnd IArrayIf
VMeasure=VArrayAnd IMeasure=IArray, then terminate, provide malfunction coding and decode, i.e., generate phase under non-uniform irradiation degree
In the case where same power, electric current, voltage, illustrates that the calculated fault type of model is consistent with actual conditions, that is, diagnose fault
For open-circuit fault or short trouble;Otherwise, next step iteration is carried out, goes to step 8;If iteration exceeds maximum number of iterations limit herein
Iter processed, then terminate, and illustrates that fault type is non-open-circuit fault or short trouble, as high loss is abnormal.
A further improvement of the present invention lies in that the photovoltaic group string fault method using Si-01TC-T irradiance sensor and
Photovoltaic monitoring system collects related data, and open-circuit failure and short trouble carry out binary equivalent coding, utilizes enhancing cigarette
Flower algorithm optimizes.
A further improvement of the present invention lies in that calculation formula is as follows in step 2:
Wherein,It is the output electric current of component i,It is the output power of component i,For in component i
Number of modules in parallel,For concatenated number of modules in component i,For module photogenerated current in component i,For in component i
Wired in parallel diode reverse saturation current,For the equivalent series resistance of component i,For the equivalent parallel electricity of component i
Resistance;For environment temperature at component i, NOCT is battery nominal operating temperature.
A further improvement of the present invention lies in that calculation formula is as follows in step 3:
Wherein, Pi(Ti,Gi) indicating i-th of component theory Maximum Power Output, it is the letter about irradiation level and temperature
Number;The set of short trouble in SC expression group string, OC indicate the set of open-circuit fault;For open circuit component i's and union II pole
The loss that pipe generates, To organize the branch current gone here and there where component i,For the parallel diode conducting voltage of component i.
The present invention has following beneficial technical effect:
The present invention by measure and analysis photovoltaic array actual current, voltage, power, photovoltaic plate temperature and irradiation level this
A little electric parameters and environmental parameter are established fault equivalence model, and are carried out using based on the optimal model for strengthening fireworks algorithm
Photovoltaic group string fault diagnosis.This method can not only carry out short trouble and the open-circuit fault diagnosis of photovoltaic group string, and can be with
Short trouble and open-circuit fault are positioned and distinguished under non-uniform irradiation degree, with model is simple, robustness is good, economical and practical
The good advantage with real-time performance;So that photovoltaic generating system ensures optimal operational condition under the premise of safety and stability, and it is
Operation maintenance personnel provides reliable fault early warning, and situation is failed to report and reported by mistake in reduction.
In conclusion the present invention enables to the staff of photovoltaic plant operation and maintenance in time to the photovoltaic to break down
Component arrangement maintenance ensure that the economic benefit of photovoltaic power plant so that realizing photovoltaic module continues optimum state operation.
Specific embodiment
Below with reference to photovoltaic group string, fault diagnosis is that the present invention will be described in more detail for example under non-uniform irradiation degree.
A kind of photovoltaic group string method for diagnosing faults under non-uniform irradiation degree provided by the invention, comprising the following steps:
Step 1: for the photovoltaic group string for sharing 10 block assemblies, acquiring photovoltaic group string output power PMeasured=
1070.50W, voltage VMeasured=155.20V and electric current IMeasuredThe irradiation level G of component i in=6.90A and photovoltaic group stringi
With component temperature Ti.Such as table 1.
1 environmental parameter of table
Step 2: according to i-th piece of photovoltaic module theory of irradiation level and temperature computation of photovoltaic equivalent-circuit model and acquisition
Maximum Power Output.
Wherein,It is the output electric current of component i,It is the output power of component i,For in component i
Number of modules in parallel,For concatenated number of modules in component i,For module photogenerated current in component i,For in component i
Wired in parallel diode reverse saturation current,For the equivalent series resistance of component i,For the equivalent parallel electricity of component i
Resistance.For environment temperature at component i, NOCT is battery nominal operating temperature.
Here, NOCT=48 DEG C,
Step 3: calculating photovoltaic group string theory output power PArray.I.e. photovoltaic group string output power is equal to light in photovoltaic group string
Lie prostrate the summation of component output power.
Wherein, Pi(Ti,Gi) indicating i-th of component theory Maximum Power Output, it is the letter about irradiation level and temperature
Number.The set of short trouble in SC expression group string, OC indicate the set of open-circuit fault.For open circuit component i's and union II pole
The loss that pipe generates, To organize the branch current gone here and there where component i,For the parallel diode conducting voltage of component i.
Step 4: determining fault threshold ε=0.1%PMeasured, and judge whether photovoltaic group string has an exception.If
PMeasured+ ε < PArray, then illustrate that a group string is operating abnormally, carry out step 5;Otherwise illustrate that group string is without exception.
Step 5: establishing fault equivalence model.For the short trouble and open-circuit fault in photovoltaic generating system, both
It can be approximated to be a kind of switch models, short trouble is equivalent to photovoltaic panel parallel diode when occurring and is short-circuited, that is, is equivalent to light
Lie prostrate the switch K1 of a plate closure in parallel;It is then equivalent to the switch K2 of photovoltaic output one disconnection of series connection when open circuit, has cut off light
Lie prostrate the external output channel of plate.And normal operation is then equivalent to switch K1 disconnection, switch K2 closure.
Step 6: building fitness function, for characterizing the fitness size of each fireworks, the value of the fitness is determined
The spark quantity and burst radius that fireworks give out.
Fitness=f (NSC,NOC, SC, OC)=| PMeasured-PArray|
It is based on the optimization object function for strengthening fireworks algorithm then:
minf(NSC,NOC, SC, OC)=| PMeasured-PArray|
Wherein NSCFor short trouble number, NOCFor open-circuit fault number.
Step 7: initialization fireworks population randomly chooses M '=20 fireworks as preliminary iteration cigarette in feasible domain space
Seeds of flowering plants group, and initialize maximum number of iterations iter=50.
Step 8: calculating the fitness size Fitness of fireworks i in fireworks populationi.Then it is normalized, obtainsIt will retain in populationM of (taking λ=0.6)
Group carries out next step operation.
Step 9: design explosion operator, is (N according to each positionSC,NOC,SC,OC)iThe fitness of fireworks individual is big
It is small, calculate the explosion initial radium R=ρ Fitness of the fireworks, wherein ρ=5 be burst radius coefficient of dilatation, the coefficient by
The specific size range of feasible zone determines, formula show fitness it is smaller be more preferably fireworks, then its explosive range will reduction,
Carry out search refinement;And fitness it is big fireworks burst radius it is big, carry out range exploitation.It is similarly big according to fireworks individual adaptation degree
It is small, calculate explosive spark numberTake σn=30, wherein θ=8 are the flexible system of explosive spark number
Number, which is determined by the size of fireworks population, while adding the normal random number of 0 mean value again.To each spark
Radius be Ri=R+N (0, σR), take σR=12 increase random numbers are the diversity in order to improve spark, while generating one
The deflection θ that even distribution generatesi=U (0,2 π).This is arrived, the position of spark i can be obtained are as follows:
(NSC,NOC,SC,OC)ii=(NSC,NOC,SC,OC)i
+(NSC,NOC,(Ricosθi)bin,(Risinθi)bin)i
Wherein, (Ricosθi)binIt indicates real number RicosθiSwitch to N binary forms, and is extracted as the position group of " 1 "
At set SC;Similarly (Risinθi)binBy real number RisinθiSwitch to N binary forms, and is extracted as the position number composition collection of " 1 "
Close OC.
Step 10: design Gaussian mutation operator carries out the spark that fireworks are broken into high to increase spark longitudinal direction diversity
This mutation operation selects n ' < n spark to carry out Gaussian mutation, that is, R as variation populationi=N (R, σR).Wherein this
In i-th spark variation population in.To the spark position of variation are as follows:
(NSC,NOC,SC,OC)′i=N (R, σR)×(NSC,NOC,SC,OC)i。
Step 11: design spark Collision operator strengthens the diversity of spark in fireworks onset process, i.e. fireworks exploded
Cheng Zhong, different spark collisions will generate new spark.That is former generation's spark (NSC,NOC,SC,OC)1(NSC,NOC,SC,OC)2It touches
It hits and generates offspring's spark (NSC,NOC,SC,OC)′1(NSC,NOC,SC,OC)′2。
The random number n between one 0 to 1 is randomly generated first, multinomial probability distribution is recycled to generate two random numbers
γ1And γ2,
Wherein
ηcIt is a nonnegative number, takes 0.8.
To obtain offspring's spark are as follows:
(NSC,NOC,SC,OC)p1
=0.5 [((NSC,NOC,SC,OC)2+(NSC,NOC,SC,OC)1)
-γ1(|(NSC,NOC,SC,OC)2-(NSC,NOC,SC,OC)1|)]
(NSC,NOC, SC, OC)p2
=0.5 [((NSC,NOC,SC,OC)2+(NSC,NOC,SC,OC)1)
-γ2(|(NSC,NOC,SC,OC)2-(NSC,NOC,SC,OC)1|)]
Step 12: selecting fireworks of the next-generation spark as next step iteration.
Step 13: judging termination condition.If f ((NSC,NOC,SC,OC)k)≤ε obtains VArrayAnd IArrayIf
VMeasure=VArrayAnd IMeasure=IArray, then terminate, provide malfunction coding and decode, i.e., generate phase under non-uniform irradiation degree
In the case where same power, electric current, voltage, illustrates that the calculated fault type of model is consistent with actual conditions, that is, diagnose fault
For open-circuit fault or short trouble.Otherwise, next step iteration is carried out, goes to step 8.If iteration exceeds maximum number of iterations limit herein
Iter processed, then terminate, and illustrates that fault type is non-open-circuit fault or short trouble, as high loss is abnormal.
Photovoltaic group string fault diagnosis is carried out through the above process, by 31 iteration, it can be deduced that fault message are as follows: the 1st, the
Open-circuit fault occurs for 8 and No. 9 components;Short trouble occurs for the 2nd and No. 5 component.The diagnostic result can be supplied to photovoltaic
The staff of power station operation and maintenance in time overhauls the photovoltaic module arrangement broken down, holds to realize photovoltaic module
Continuous optimum state operation, ensure that the economic benefit of photovoltaic power plant.
Claims (4)
1. a kind of photovoltaic group string method for diagnosing faults under non-uniform irradiation degree, which comprises the following steps:
Step 1: for the photovoltaic group string of shared N block assembly, acquiring the output power P of photovoltaic group stringMeasured, voltage VMeasuredWith
Electric current IMeasuredAnd in photovoltaic group string component i irradiation level GiWith component temperature Ti;
Step 2: according to photovoltaic equivalent-circuit model and the irradiation level G of acquisitioniWith component temperature TiIt is calculated in photovoltaic group string
The theoretical maximum output power of component i;
Step 3: calculating photovoltaic group string theory output power PArray;I.e. photovoltaic group string output power is equal to photovoltaic group in photovoltaic group string
The summation of part output power;
Step 4: determining fault threshold ε, and judge whether photovoltaic group string has an exception;If PMeasured+ ε < PArray, then illustrate
Group string is operating abnormally, and carries out step 5;Otherwise illustrate that group string is without exception;
Step 5: establishing fault equivalence model;For the short trouble and open-circuit fault in photovoltaic generating system, both are approximate
For a kind of switch models, short trouble is equivalent to photovoltaic panel parallel diode when occurring and is short-circuited, that is, is equivalent to photovoltaic panel parallel connection
The switch K1 of one closure;It is then equivalent to the switch K2 of photovoltaic output one disconnection of series connection when open circuit, it is external to have cut off photovoltaic panel
Output channel;And normal operation is then equivalent to switch K1 disconnection, switch K2 closure;
Step 6: building fitness function, for characterizing the fitness size of each fireworks, the value of the fitness determines cigarette
The spark quantity and burst radius that flower gives out;
Fitness=f (NSC, NOC, SC, OC)=| PMeasured-PArray|
It is using the target for strengthening fireworks algorithm optimization then:
minf(NSC, NOC, SC, OC)=| PMeasured-PArray|
Wherein NSCFor short trouble number, NOCFor open-circuit fault number;
Step 7: initialization fireworks population randomly chooses a fireworks of M ' as preliminary iteration fireworks population in feasible domain space,
And initialize maximum number of iterations iter;
Step 8: calculating the fitness size Fitness of fireworks i in fireworks populationi;Then it is normalized, obtainsIt will retain in populationM group carry out down
Single stepping;
Step 9: design explosion operator, is (N according to each positionSC, NOC, SC, OC)iThe fitness size of fireworks individual, meter
The explosion initial radium R=ρ Fitness of the fireworks is calculated, wherein ρ is burst radius coefficient of dilatation, and the coefficient is by feasible zone
Specific size range determines, formula show fitness it is smaller be more preferably fireworks, then its explosive range will reduction, refined
Search;And fitness it is big fireworks burst radius it is big, carry out range exploitation;Similarly according to fireworks individual adaptation degree size, calculate
Explosive spark number Wherein θ is explosive spark number coefficient of dilatation, and the coefficient is by fireworks population
Size determine, while adding the normal random number of 0 mean value again;To which the radius of each spark is Ri=R+N (0,
σR), increasing random number is the diversity in order to improve spark, while generating the deflection θ for being uniformly distributed generationi=U (0,2
π);This is arrived, the position of spark i is obtained are as follows:
(NSC, NOC, SC, OC) 'i=(NSC, NOC, SC, OC)i
+(NSC, NOC, (Ricosθi)bin, (Risinθi)bin)i
Wherein, (Ricosθi)binIt indicates real number RicosθiSwitch to N binary forms, and is extracted as the position number composition collection of " 1 "
Close SC;Similarly (Risinθi)binBy real number RisinθiSwitch to N binary forms, and is extracted as the position number composition set of " 1 "
OC;
Step 10: design Gaussian mutation operator carries out Gauss change to the spark that fireworks are broken into increase spark longitudinal direction diversity
ETTHER-OR operation selects n ' < n spark to carry out Gaussian mutation, that is, R as variation populationi=N (R, σR);Wherein here i-th
In spark variation population;To which the spark position of variation is (NSC, NOC, SC, OC) 'i=N (R, σR)×(NSC, NOC, SC, OC)i;
Step 11: design spark Collision operator strengthens the diversity of spark in fireworks onset process, i.e., fireworks are in blast process
In, different spark collisions will generate new spark;That is former generation's spark (NSC, NOC, SC, OC)1(NSC, NOC, SC, OC)2Collision
Generate offspring's spark (NSC, NOC, SC, OC) '1(NSC, NOC, SC, OC) '2;It is randomly generated first random between one 0 to 1
Number n, recycles multinomial probability distribution to generate two random number γ1And γ2,
Wherein
ηcIt is a nonnegative number;
To obtain offspring's spark are as follows:
(NSC, NOC, SC, OC)p1
=0.5 [((NSC, NOC, SC, OC)2+(NSC, NOC, SC, OC)1)
-γ1(|(NSC, NOC, SC, OC)2-(NSC, NOC, SC, OC)1|)]
(NSC, NOC, SC, OC)p2
=0.5 [((NSC, NOC, SC, OC)2+(NSC, NOC, SC, OC)1)
-γ2(|(NSC, NOC, SC, OC)2-(NSC, NOC, SC, OC)1|)]
Step 12: selecting fireworks of the next-generation spark as next step iteration;
Step 13: judging termination condition;If f ((NSC, NOC, SC, OC)k)≤ε obtains VArrayAnd IArrayIf VMeasure=
VArrayAnd IMeasure=IArray, then terminate, provide malfunction coding and decode, i.e., under non-uniform irradiation degree generate equal-wattage,
In the case where electric current, voltage, illustrate that the calculated fault type of model is consistent with actual conditions, i.e. troubleshooting is open circuit event
Barrier or short trouble;Otherwise, next step iteration is carried out, goes to step 8;If iteration exceeds maximum number of iterations limitation iter herein,
It then terminates, illustrates that fault type is non-open-circuit fault or short trouble, as high loss is abnormal.
2. a kind of photovoltaic group string method for diagnosing faults under non-uniform irradiation degree according to claim 1, feature exist
In the photovoltaic group string fault method collects related data using Si-01TC-T irradiance sensor and photovoltaic monitoring system, and right
Open-circuit fault and short trouble carry out binary equivalent coding, are optimized using enhancing fireworks algorithm.
3. a kind of photovoltaic group string method for diagnosing faults under non-uniform irradiation degree according to claim 1 or 2, feature
It is, in step 2, calculation formula is as follows:
Wherein,It is the output electric current of component i,It is the output power of component i,It is in parallel in component i
Number of modules,For concatenated number of modules in component i,For module photogenerated current in component i,For module in component i
Parallel diode reverse saturation current,For the equivalent series resistance of component i,For the equivalent parallel resistance of component i;For environment temperature at component i, NOCT is battery nominal operating temperature.
4. a kind of photovoltaic group string method for diagnosing faults under non-uniform irradiation degree according to claim 3, feature exist
In in step 3, calculation formula is as follows:
Wherein, Pi(Ti, Gi) indicating i-th of component theory Maximum Power Output, it is the function about irradiation level and temperature;SC
The set of short trouble in expression group string, OC indicate the set of open-circuit fault;Parallel diode for the component i that opens a way produces
Raw loss, To organize the branch current gone here and there where component i,For the parallel diode conducting voltage of component i.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811168547.7A CN109409420B (en) | 2018-10-08 | 2018-10-08 | Photovoltaic string fault diagnosis method under non-uniform irradiance |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811168547.7A CN109409420B (en) | 2018-10-08 | 2018-10-08 | Photovoltaic string fault diagnosis method under non-uniform irradiance |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109409420A true CN109409420A (en) | 2019-03-01 |
CN109409420B CN109409420B (en) | 2022-05-03 |
Family
ID=65466089
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811168547.7A Active CN109409420B (en) | 2018-10-08 | 2018-10-08 | Photovoltaic string fault diagnosis method under non-uniform irradiance |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109409420B (en) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150094967A1 (en) * | 2013-09-27 | 2015-04-02 | Hitachi, Ltd. | Failure diagnosis method and failure diagnosis system for photovoltaic system |
CN104767482A (en) * | 2014-01-02 | 2015-07-08 | 上海岩芯电子科技有限公司 | Method for diagnosing aging and short circuit of photovoltaic module online |
CN105652952A (en) * | 2016-04-18 | 2016-06-08 | 中国矿业大学 | Maximum power point tracking method for photovoltaic power generation system based on fireworks algorithm |
CN106776469A (en) * | 2017-03-17 | 2017-05-31 | 沈阳航空航天大学 | A kind of fireworks algorithm based on simulated annealing Yu Gauss disturbance |
CN107102223A (en) * | 2017-03-29 | 2017-08-29 | 江苏大学 | NPC photovoltaic DC-to-AC converter method for diagnosing faults based on improved hidden Markov model GHMM |
CN107579707A (en) * | 2017-10-13 | 2018-01-12 | 江苏大学 | A kind of diagnosing failure of photovoltaic array method based on parameter identification |
CN107658840A (en) * | 2017-09-30 | 2018-02-02 | 山东科技大学 | Distribution network failure restoration methods based on A* algorithms Yu fireworks algorithm |
CN108062571A (en) * | 2017-12-27 | 2018-05-22 | 福州大学 | Diagnosing failure of photovoltaic array method based on differential evolution random forest grader |
US20180238951A1 (en) * | 2016-09-07 | 2018-08-23 | Jiangnan University | Decision Tree SVM Fault Diagnosis Method of Photovoltaic Diode-Clamped Three-Level Inverter |
-
2018
- 2018-10-08 CN CN201811168547.7A patent/CN109409420B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150094967A1 (en) * | 2013-09-27 | 2015-04-02 | Hitachi, Ltd. | Failure diagnosis method and failure diagnosis system for photovoltaic system |
CN104767482A (en) * | 2014-01-02 | 2015-07-08 | 上海岩芯电子科技有限公司 | Method for diagnosing aging and short circuit of photovoltaic module online |
CN105652952A (en) * | 2016-04-18 | 2016-06-08 | 中国矿业大学 | Maximum power point tracking method for photovoltaic power generation system based on fireworks algorithm |
US20180238951A1 (en) * | 2016-09-07 | 2018-08-23 | Jiangnan University | Decision Tree SVM Fault Diagnosis Method of Photovoltaic Diode-Clamped Three-Level Inverter |
CN106776469A (en) * | 2017-03-17 | 2017-05-31 | 沈阳航空航天大学 | A kind of fireworks algorithm based on simulated annealing Yu Gauss disturbance |
CN107102223A (en) * | 2017-03-29 | 2017-08-29 | 江苏大学 | NPC photovoltaic DC-to-AC converter method for diagnosing faults based on improved hidden Markov model GHMM |
CN107658840A (en) * | 2017-09-30 | 2018-02-02 | 山东科技大学 | Distribution network failure restoration methods based on A* algorithms Yu fireworks algorithm |
CN107579707A (en) * | 2017-10-13 | 2018-01-12 | 江苏大学 | A kind of diagnosing failure of photovoltaic array method based on parameter identification |
CN108062571A (en) * | 2017-12-27 | 2018-05-22 | 福州大学 | Diagnosing failure of photovoltaic array method based on differential evolution random forest grader |
Non-Patent Citations (5)
Title |
---|
GUOJIANG XIONG等: "Parameter extraction of solar photovoltaic models using an improved whale optimization algorithm", 《ENERGY CONVERSION AND MANAGEMENT》 * |
MOHAMED HASSAN ALI等: "Real Time Fault Detection in Photovoltaic Systems", 《ENERGY PROCEDIA》 * |
叶雯雯等: "自适应爆炸半径和差分变异的增强烟花算法", 《应用数学进展》 * |
周仕豪: "考虑降低暂态电压失稳风险的动态无功优化配置方法", 《电力系统保护与控制》 * |
谭营等: "烟花算法研究进展", 《智能系统学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109409420B (en) | 2022-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102288856B (en) | Photovoltaic polar plate fault damage detection device based on wireless communication and method of same | |
CN104504607A (en) | Method for diagnosing photovoltaic power station faults on the basis of fuzzy clustering algorithm | |
CN104767481B (en) | Method and system for monitoring working state of solar photovoltaic power station | |
CN110096764B (en) | Method for identifying and optimizing fragile line of electric-gas coupling system | |
CN107978807A (en) | A kind of battery detecting and maintaining method and system | |
CN102520369A (en) | Online monitoring system for storage battery of fan variable pitch propeller | |
CN112083270B (en) | Wind power plant current collection line single-phase earth fault line selection method based on correlation coefficient | |
CN107766596B (en) | Low-voltage ride-through capability simulation evaluation method based on typical fault working condition set | |
CN106788243B (en) | The method for realizing the string state recognition of photovoltaic plant group based on monitor supervision platform | |
Chen et al. | A novel fault diagnosis method of PV based-on power loss and IV characteristics | |
CN105117591A (en) | Power device running status analysis method | |
Djalab et al. | Robust method for diagnosis and detection of faults in photovoltaic systems using artificial neural networks | |
Hare et al. | A review of faults and fault diagnosis in micro-grids electrical energy infrastructure | |
Yun et al. | Research on fault diagnosis of photovoltaic array based on random forest algorithm | |
Zaki et al. | Detection and localization the open and short circuit faults in PV system: A MILP approach | |
CN104779909B (en) | Method and system for monitoring working state of solar photovoltaic power station | |
CN109298228A (en) | A kind of Intelligence Diagnosis method and system based on photovoltaic group string current anomaly | |
Alajmi et al. | Fault detection and localization in solar photovoltaic arrays framework: hybrid methods of data-analysis and a network of voltage-current sensors | |
CN102608492A (en) | System and method for detecting fuse state of combiner box | |
CN109409420A (en) | A kind of photovoltaic group string method for diagnosing faults under non-uniform irradiation degree | |
CN110635767A (en) | Photovoltaic system and method for inhibiting potential attenuation and monitoring insulation | |
CN211350740U (en) | Cylindrical battery PACK box series conductor with alarm function | |
CN105158697B (en) | The method and device of energy-storage units status assessment | |
CN107271916A (en) | A kind of cell plate group string method for detecting health status | |
Wang et al. | Adaptability analysis for fault component protection under large scale wind power penetration |
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 |