CN116388190B - Distributed photovoltaic harmonic quantitative evaluation method, system, terminal equipment and medium - Google Patents
Distributed photovoltaic harmonic quantitative evaluation method, system, terminal equipment and medium Download PDFInfo
- Publication number
- CN116388190B CN116388190B CN202310658384.5A CN202310658384A CN116388190B CN 116388190 B CN116388190 B CN 116388190B CN 202310658384 A CN202310658384 A CN 202310658384A CN 116388190 B CN116388190 B CN 116388190B
- Authority
- CN
- China
- Prior art keywords
- feeder line
- harmonic
- feeder
- current
- quantitative evaluation
- 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
- 238000011158 quantitative evaluation Methods 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000011159 matrix material Substances 0.000 claims abstract description 70
- 230000008878 coupling Effects 0.000 claims abstract description 22
- 238000010168 coupling process Methods 0.000 claims abstract description 22
- 238000005859 coupling reaction Methods 0.000 claims abstract description 22
- 238000000926 separation method Methods 0.000 claims abstract description 22
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000011002 quantification Methods 0.000 claims description 7
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 claims 1
- 238000009826 distribution Methods 0.000 description 7
- 230000002087 whitening effect Effects 0.000 description 3
- 238000010248 power generation Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012880 independent component analysis Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/01—Arrangements for reducing harmonics or ripples
-
- 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
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/22—The renewable source being solar energy
- H02J2300/24—The renewable source being solar energy of photovoltaic origin
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Power Engineering (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computing Systems (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Photovoltaic Devices (AREA)
Abstract
The application relates to the technical field of distributed photovoltaic systems, and particularly discloses a distributed photovoltaic harmonic quantitative evaluation method, a system, terminal equipment and a medium; the method comprises the following steps: acquiring voltage at a public coupling point of a concerned time point and current data of each feeder line to form a feeder line current matrix; calculating the mutual information value of each feeder line current; selecting a feeder line with a mutual information value larger than a set value, adding current data of the feeder line to be used as a new column of data in a feeder line current matrix, and deleting original current data of the feeder line in the feeder line current matrix to form a new feeder line current matrix; calculating a new feeder current matrix by using an independent vector analysis algorithm to obtain a decoupled separation matrix; based on the voltage at the point of common coupling and the separation matrix, a result is obtainedhHarmonic quantitative evaluation results of subharmonics. The application can accurately evaluate the harmonic responsibility of a plurality of feeder lines with higher harmonic current correlation under the distributed photovoltaic access.
Description
Technical Field
The application relates to the technical field of distributed photovoltaic systems, in particular to a distributed photovoltaic harmonic quantitative evaluation method, a system, terminal equipment and a medium.
Background
The statements in this section merely provide background information related to the present disclosure and may not necessarily constitute prior art.
The photovoltaic power generation has intermittence, randomness and volatility, so that the state of the power distribution network is complex and changeable. Meanwhile, a large number of power electronic devices are connected into the power distribution network through photovoltaic power generation, a large number of low-frequency harmonics and switching order harmonics are introduced, so that the power quality problem of the power distribution network is more serious, and the method has important significance in harmonic quantitative evaluation research on the premise of improving the power quality.
The centralized harmonic quantitative evaluation is an algorithm for dividing the contribution of harmonic source loads connected by a plurality of feeder lines to harmonic voltages at a public coupling point, and because the power distribution network mostly presents a distributed structure, the centralized harmonic quantitative evaluation method is more common, however, in actual situations, current between the feeder lines has a coupling phenomenon, and the harmonic quantitative evaluation by adopting current data measured on the feeder lines is not necessarily in direct proportion to an actual value.
In order to obtain the theoretical current of each feeder line, a blind source separation algorithm is mainly adopted to separate the independent quantity at present, so that the theoretical current of each feeder line is obtained. However, the blind source separation algorithm needs to ensure the independence of each data, and under the scene of large-scale photovoltaic access, the current of each feeder line with photovoltaic access generally has higher correlation, and at the moment, the current cannot be directly calculated by using the blind source separation algorithm; in addition, in the centralized harmonic quantitative evaluation process, independence is also required between data, and the prior art needs to perform complex processing on the data to obtain a final harmonic quantitative evaluation result.
Disclosure of Invention
In order to solve the problems, the application provides a distributed photovoltaic harmonic quantitative evaluation method, a system, terminal equipment and a medium, independent component analysis is carried out through mutual information value verification to obtain theoretical currents of all feeder lines, and then harmonic quantitative evaluation results of all feeder lines are obtained through linear regression, weighted correlation coefficients and current fluctuation coefficients, so that the harmonic quantitative evaluation results are closer to the true level.
In some embodiments, the following technical scheme is adopted:
a distributed photovoltaic harmonic quantitative evaluation method, comprising:
acquiring voltage at a public coupling point of a concerned time point and current data of each feeder line to form a feeder line current matrix;
calculating mutual information values of currents of all feeder lines based on current data of all feeder lines;
selecting a feeder line with a mutual information value larger than a set value, adding current data of the feeder line to be used as a new column of data in a feeder line current matrix, and deleting original current data of the feeder line in the feeder line current matrix to form a new feeder line current matrix;
calculating a new feeder current matrix by using an independent vector analysis algorithm to obtain a decoupled separation matrix;
based on the voltage at the point of common coupling and the separation matrix, a result is obtainedhHarmonic quantitative evaluation results of subharmonics.
Wherein the saidhIn the harmonic quantitative evaluation result of subharmonic, the firstkColumn harmonic quantitative evaluation resultsThe sum of harmonic quantitative evaluation results of feeder lines with mutual information values larger than a set value is obtained;
obtaininghAfter harmonic quantitative evaluation results of subharmonic, supposing that the sum of feeder currents with mutual information values larger than a set value is decoupled, the obtained ith current data is thatIn the feeder line with mutual information value greater than the set value, calculating the firstmOf feed linesFeeder current and->Is a correlation coefficient of (2);
and carrying out harmonic quantitative evaluation on the feeder line with the mutual information value larger than a set value based on the correlation coefficient and the current fluctuation coefficient.
In other embodiments, the following technical solutions are adopted:
a distributed photovoltaic harmonic quantitative evaluation system comprising:
the data acquisition module is used for acquiring the voltage at the point of public coupling at the point of attention time and the current data of each feeder line to form a feeder line current matrix;
the mutual information value calculation module is used for calculating the mutual information value of the current of each feeder line based on the current data of each feeder line;
the theoretical current calculation module is used for selecting a feeder line with a mutual information value larger than a set value, adding current data of the feeder line to be used as a new column of data in a feeder line current matrix, and deleting original current data of the feeder line in the feeder line current matrix to form a new feeder line current matrix; calculating a new feeder current matrix by using an independent vector analysis algorithm to obtain a decoupled separation matrix;
a harmonic quantitative evaluation module for obtaining based on the voltage at the point of common coupling and the separation matrixhHarmonic quantitative evaluation results of subharmonics.
In other embodiments, the following technical solutions are adopted:
a terminal device comprising a processor and a memory, the processor for implementing instructions; the memory is used to store a plurality of instructions adapted to be loaded by the processor and to perform the distributed photovoltaic harmonic quantification assessment method described above.
In other embodiments, the following technical solutions are adopted:
a computer readable storage medium having stored therein a plurality of instructions adapted to be loaded by a processor of a terminal device and to perform the above-described distributed photovoltaic harmonic quantification assessment method.
Compared with the prior art, the application has the beneficial effects that:
(1) On the basis of traditional multi-feeder harmonic quantitative evaluation, the application considers the situation that the relevance of partial feeder harmonic sources is stronger under the condition of photovoltaic access and the traditional blind source separation algorithm is not applicable, restores the real current of each feeder line through the mutual information value and the independent vector analysis algorithm, and simultaneously considers the fluctuation factor of the harmonic current so that the harmonic quantitative evaluation result of each feeder line is more approximate to the real value; the power distribution system is particularly suitable for a power distribution system under high-proportion distributed photovoltaic access.
(2) Firstly, calculating the sum of harmonic quantitative evaluation results of feeder lines with lower correlation and harmonic quantitative evaluation results of several feeder lines with higher correlation; and then respectively calculating harmonic quantitative evaluation results of the feeder lines with higher correlation. The harmonic responsibility of a plurality of feeder lines with higher harmonic current correlation under the distributed photovoltaic access can be accurately evaluated.
Additional features and advantages of the application will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the application.
Drawings
FIG. 1 is a schematic diagram of a large-scale photovoltaic parallel-grid multi-feeder structure in an embodiment of the present application;
fig. 2 is a flowchart of a method for quantitatively evaluating distributed photovoltaic harmonics in an embodiment of the present application.
Detailed Description
It should be noted that the following detailed description is illustrative and is intended to provide further explanation of the application. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the present application. As used herein, the singular is also intended to include the plural unless the context clearly indicates otherwise, and furthermore, it is to be understood that the terms "comprises" and/or "comprising" when used in this specification are taken to specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof.
Example 1
In one or more embodiments, a distributed photovoltaic harmonic quantitative evaluation method is disclosed, and in combination with fig. 2, the method specifically includes the following steps:
s101: acquiring voltage at the point of common coupling at the point of interest and current data of each feeder line to form a feeder line current matrix;
As a specific example, fig. 1 shows a large-scale photovoltaic co-grid multi-feeder structure, with reference to a point of common coupling (Point of common coupling, PCC), dividing the system side and the user side; wherein, byhSubharmonic is exemplified by the system sidehSub-equivalent harmonic impedanceAndhsubequivalent harmonic current source->The parallel connection is formed, and the user side is formed by connecting a plurality of feeder lines in parallel and comprises a user feeder line 1, a photovoltaic feeder line 2, … and light Fu Kuixian n; wherein the photovoltaic feed lines comprise a relatively high proportion of the feed lines, each feed line comprisinghSub-equivalent harmonic impedance->Andhsubequivalent harmonic current source->Is formed by connecting the two parts in parallel,irepresenting the first of the usersiAnd (5) a feeder line.
Feeder current matrix of this embodimentThe method comprises the following steps:
wherein ,represent the firstnStrip feederhSubharmonic at time pointt n Current value at the time.
S102: calculating mutual information values of currents of all feeder lines based on current data of all feeder lines; the specific calculation formula is as follows:
wherein ,representing feederiAnd feeder linejA kind of electronic devicehMutual information value of the secondary current.
Representing feederiA kind of electronic devicehThe entropy of the secondary current data is calculated by the following formula:
representing feederiA kind of electronic devicehProbability value of occurrence of a certain number of secondary current data.
Representing feederiAnd feeder linejA kind of electronic devicehThe calculation formula of the joint entropy value of the secondary current data is as follows:
representing feederiAnd feeder linejA kind of electronic devicehJoint probability of secondary current data.
S103: selecting a feeder line with the mutual information value larger than a set value, adding current data of the feeder line to be used as a new column of data in a feeder line current matrix, and deleting original current data of the feeder line in the feeder line current matrix to form a new feeder line current matrix.
In this embodiment, for feeders having mutual information values greater than a set value, the correlation degree is considered to be relatively high; the set value can be selected according to actual needs; such as: in this embodiment, the set value is selected to be 0.3.
The present embodiment forms a new feeder current matrixThe method comprises the following steps:
wherein ,/>Is thatt 1 The time mutual information value is larger than the sum of feeder currents of the set values,lrepresenting a commonlAnd the mutual information value is larger than the feeder line of the set value.
S104: and calculating the new feeder current matrix by using an independent vector analysis algorithm to obtain a decoupled separation matrix.
The independent vector analysis separates signals by utilizing the independence in the data, and on the premise that the signals have different spatial distributions, the mixed signals are separated by performing matrix transformation on the data on the assumption that the signals are irrelevant in time and frequency.
Independent vector analysis algorithm model:
in the formula ,for observing signals, +.>For the source signal, a is the coefficient matrix.
The observation signals represent a new feeder current matrix formed after mutual information value calculation; the source signal represents the harmonic current emitted when each harmonic source acts alone.
And (3) carrying out averaging treatment on the observation signals to obtain:
in the formula ,is the mean signal of the observed signals.
Whitening the data, removing correlation between the metrology signals:
in the formula ,for whitening matrix +.>,/>Is->Eigenvalues of matrix +.>Is thatEigenvalue matrix of matrix,/>Matrix after data whitening, ++>And->Is a matrix of the same dimension.
Finally, obtaining a separation matrix (i.e. a theoretical current matrix) after iterative solution by a Niuzhan method:
in the formula ,Wfor an iteratively solved parameter matrix, whereinIs the first of (2)kThe sum of feeder currents with mutual information values larger than a set value is listed as the current after decoupling.
S105: based on the voltage at the point of common coupling and the separation matrix, a result is obtainedhHarmonic quantitative evaluation results of subharmonics.
In the embodiment, based on voltage data at the public coupling point, regression coefficients of all feeder lines are obtained by using a regression method;
wherein ,is n 1hA subpublic coupling point harmonic voltage matrix, E is an identity matrix of n 1,Z n is a feeder linenIs used to determine the regression coefficients of (a),U s is a harmonic voltage on the system side.
Then, based on the separation matrix, obtainhSub-resonanceThe harmonic quantitative evaluation result of the wave is:
in the formula ,is->Is the first of (2)lColumn (1)kColumn harmonic quantitative evaluation results->The sum of harmonic quantitative evaluation results of feeder lines with higher mutual information values.
S106: obtaininghAnd after the harmonic quantitative evaluation result of the subharmonic, carrying out harmonic quantitative evaluation on the feeder line with the mutual information value larger than the set value.
Specifically, after decoupling the sum of feeder currents with mutual information values greater than a set value, the obtained ith current data is assumed to beIn the feeder line with mutual information value greater than the set value, the firstmFeeder current and +.>Correlation coefficient of->The method specifically comprises the following steps:
wherein ,/>In the feeder line with mutual information value larger than the set value, the firstmFirst of the feeder linesi-individual current data; />In the feeder line with mutual information value larger than the set value, the firstmThe average value of the currents of the feeder lines; />And the average value of all current data is obtained after decoupling the sum of feeder currents with mutual information values larger than a set value.
Based on the correlation coefficient and the current fluctuation coefficient, carrying out harmonic quantitative evaluation on the feeder line with the mutual information value larger than a set value, wherein the method specifically comprises the following steps:
wherein ,is->Is a cumulative value of (a); />Is the current fluctuation coefficient; />In the feeder line with mutual information value larger than the set value, the firstmThe average value of the currents of the feeder lines; />In the feeder line with mutual information value larger than the set value, the firstmThe maximum current value of the strip feed line. />In the feeder line representing that the mutual information value is larger than the set value, the firstmOf feed lineshAnd quantitatively evaluating the result of subharmonic.
The process finally realizes the quantitative evaluation of the harmonic wave of each feeder line.
Example two
In one or more embodiments, a distributed photovoltaic harmonic quantitative evaluation system is disclosed, comprising in particular:
the data acquisition module is used for acquiring the voltage at the point of public coupling at the point of attention time and the current data of each feeder line to form a feeder line current matrix;
the mutual information value calculation module is used for calculating the mutual information value of the current of each feeder line based on the current data of each feeder line;
the theoretical current calculation module is used for selecting a feeder line with a mutual information value larger than a set value, adding current data of the feeder line to be used as a new column of data in a feeder line current matrix, and deleting original current data of the feeder line in the feeder line current matrix to form a new feeder line current matrix; calculating a new feeder current matrix by using an independent vector analysis algorithm to obtain a decoupled separation matrix;
a harmonic quantitative evaluation module for obtaining based on the voltage at the point of common coupling and the separation matrixhHarmonic quantitative evaluation results of subharmonics.
The specific implementation of each module has been described in detail in the first embodiment, and will not be described in detail here.
Example III
In one or more embodiments, a terminal device is disclosed that includes a server including a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing the distributed photovoltaic harmonic quantitative evaluation method of embodiment one when executing the program. For brevity, the description is omitted here.
It should be understood that in this embodiment, the processor may be a central processing unit CPU, and the processor may also be other general purpose processors, digital signal processors DSP, application specific integrated circuits ASIC, off-the-shelf programmable gate array FPGA or other programmable logic device, discrete gate or transistor logic devices, discrete hardware components, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory may include read only memory and random access memory and provide instructions and data to the processor, and a portion of the memory may also include non-volatile random access memory. For example, the memory may also store information of the device type.
In implementation, the steps of the above method may be performed by integrated logic circuits of hardware in a processor or by instructions in the form of software.
Example IV
In one or more embodiments, a computer-readable storage medium is disclosed, in which a plurality of instructions are stored, the instructions being adapted to be loaded by a processor of a terminal device and to perform the distributed photovoltaic harmonic quantification assessment method described in embodiment one.
While the foregoing description of the embodiments of the present application has been presented in conjunction with the drawings, it should be understood that it is not intended to limit the scope of the application, but rather, it is intended to cover all modifications or variations within the scope of the application as defined by the claims of the present application.
Claims (9)
1. A distributed photovoltaic harmonic quantitative evaluation method, comprising:
acquiring voltage at a public coupling point of a concerned time point and current data of each feeder line to form a feeder line current matrix;
calculating mutual information values of currents of all feeder lines based on current data of all feeder lines;
selecting a feeder line with a mutual information value larger than a set value, adding current data of the feeder line to be used as a new column of data in a feeder line current matrix, and deleting original current data of the feeder line in the feeder line current matrix to form a new feeder line current matrix;
calculating a new feeder current matrix by using an independent vector analysis algorithm to obtain a decoupled separation matrix;
based on the voltage at the point of common coupling and the separation matrix, a result is obtainedhHarmonic quantitative evaluation results of subharmonics;
the saidhHarmonic quantitative evaluation result of subharmonicIn (1)kColumn harmonic quantitative evaluation resultsThe sum of harmonic quantitative evaluation results of feeder lines with mutual information values larger than a set value is obtained;
obtaininghAfter harmonic quantitative evaluation results of subharmonic, supposing that the sum of feeder currents with mutual information values larger than a set value is decoupled, the obtained ith current data is thatIn the feeder line with mutual information value greater than the set value, calculating the firstmFeeder current and +.>Is a correlation coefficient of (2);
and carrying out harmonic quantitative evaluation on the feeder line with the mutual information value larger than a set value based on the correlation coefficient and the current fluctuation coefficient.
2. A distributed photovoltaic harmonic quantification assessment method according to claim 1, wherein the first calculation is performedmFeeder current of strip feederCorrelation coefficient of->The method specifically comprises the following steps:
wherein ,in the feeder line with mutual information value larger than the set value, the firstmFirst of the feeder linesi-individual current data; />In the feeder line with mutual information value larger than the set value, the firstmThe average value of the currents of the feeder lines; />And the average value of all current data is obtained after decoupling the sum of feeder currents with mutual information values larger than a set value.
3. The distributed photovoltaic harmonic quantitative evaluation method according to claim 2, wherein harmonic quantitative evaluation is performed on a feeder line with a mutual information value larger than a set value, specifically:
wherein ,is->Is a cumulative value of (a); />Is the current fluctuation coefficient; />In the feeder line with mutual information value larger than the set value, the firstmThe average value of the currents of the feeder lines; />In the feeder line with mutual information value larger than the set value, the firstmA current maximum value of the strip feeder; />In the feeder line representing that the mutual information value is larger than the set value, the firstmOf feed lineshAnd quantitatively evaluating the result of subharmonic.
4. The distributed photovoltaic harmonic quantitative evaluation method according to claim 1, wherein the calculating of the mutual information value of each feeder line current based on the current data of each feeder line is specifically as follows:
wherein ,representing feederiA kind of electronic devicehEntropy of subharmonic current; />Representing feederjA kind of electronic devicehEntropy of subharmonic current>Representing feederiAnd feeder linejA kind of electronic devicehA joint entropy value of the subharmonic current; />Representing feederiAnd feeder linejA kind of electronic devicehMutual information value of subharmonic current.
5. A distributed photovoltaic harmonic quantification evaluation method according to claim 1, wherein the voltage at the point of common coupling and the separation matrix are based onhThe harmonic quantitative evaluation result of subharmonic is specifically:
wherein ,is the firstlHarmonic quantitative evaluation results of the strip feeder lines; />Is a feeder linelRegression coefficients of (a); />Is thatIs the first of (2)lColumn (S)/(S)>Is a separation matrix; />Represents n 1hA secondary common coupling point harmonic voltage matrix.
6. A distributed photovoltaic harmonic quantitative evaluation method as claimed in claim 5 wherein the regression coefficients for each feed line are solved using:
wherein ,Z 1 、Z 2 、…、Z n the feed lines 1,2 …,ne is an identity matrix of n 1,U s is a harmonic voltage on the system side.
7. A distributed photovoltaic harmonic quantitative evaluation system, comprising:
the data acquisition module is used for acquiring the voltage at the point of public coupling at the point of attention time and the current data of each feeder line to form a feeder line current matrix;
the mutual information value calculation module is used for calculating the mutual information value of the current of each feeder line based on the current data of each feeder line;
the theoretical current calculation module is used for selecting a feeder line with a mutual information value larger than a set value, adding current data of the feeder line to be used as a new column of data in a feeder line current matrix, and deleting original current data of the feeder line in the feeder line current matrix to form a new feeder line current matrix; calculating a new feeder current matrix by using an independent vector analysis algorithm to obtain a decoupled separation matrix;
a harmonic quantitative evaluation module for obtaining based on the voltage at the point of common coupling and the separation matrixhHarmonic quantitative evaluation results of subharmonics;
the saidhIn the harmonic quantitative evaluation result of subharmonic, the firstkColumn harmonic quantitative evaluation resultsThe sum of harmonic quantitative evaluation results of feeder lines with mutual information values larger than a set value is obtained;
obtaininghAfter harmonic quantitative evaluation results of subharmonic, supposing that the sum of feeder currents with mutual information values larger than a set value is decoupled, the obtained ith current data is thatIn the feeder line with mutual information value greater than the set value, calculating the firstmFeeder current and +.>Is a correlation coefficient of (2);
and carrying out harmonic quantitative evaluation on the feeder line with the mutual information value larger than a set value based on the correlation coefficient and the current fluctuation coefficient.
8. A terminal device comprising a processor and a memory, the processor for implementing instructions; a memory for storing a plurality of instructions, wherein the instructions are adapted to be loaded by a processor and to perform the distributed photovoltaic harmonic quantification method of any of claims 1-6.
9. A computer readable storage medium having stored therein a plurality of instructions adapted to be loaded by a processor of a terminal device and to perform the distributed photovoltaic harmonic quantification assessment method of any of claims 1-6.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310658384.5A CN116388190B (en) | 2023-06-06 | 2023-06-06 | Distributed photovoltaic harmonic quantitative evaluation method, system, terminal equipment and medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310658384.5A CN116388190B (en) | 2023-06-06 | 2023-06-06 | Distributed photovoltaic harmonic quantitative evaluation method, system, terminal equipment and medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116388190A CN116388190A (en) | 2023-07-04 |
CN116388190B true CN116388190B (en) | 2023-08-18 |
Family
ID=86977241
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310658384.5A Active CN116388190B (en) | 2023-06-06 | 2023-06-06 | Distributed photovoltaic harmonic quantitative evaluation method, system, terminal equipment and medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116388190B (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412242A (en) * | 2013-08-01 | 2013-11-27 | 西南交通大学 | Method for locating harmonic source based on independent rapid component analysis and mutual information |
CN103424620A (en) * | 2013-08-08 | 2013-12-04 | 华南理工大学 | 10kV distribution network harmonic source identification method based on time sequence relevancy |
WO2014075581A1 (en) * | 2012-11-16 | 2014-05-22 | 天津大学 | Current-mode matching method for non-invasive power load monitoring and decomposition |
WO2014206007A1 (en) * | 2013-06-26 | 2014-12-31 | 国家电网公司 | Wind power plant harmonic suppression method, wind power generator sets, and wind power station |
CN104502704A (en) * | 2015-01-20 | 2015-04-08 | 西南交通大学 | Harmonic contribution calculating method applicable to background harmonic voltage change |
CN106208041A (en) * | 2016-07-18 | 2016-12-07 | 西南交通大学 | A kind of Distribution Network Harmonics current estimation method based on piecewise constant independent component analysis |
CN108089059A (en) * | 2017-12-19 | 2018-05-29 | 西南交通大学 | A kind of harmonic impedance evaluation method based on sparse independent component analysis |
CN114660362A (en) * | 2022-05-19 | 2022-06-24 | 四川大学 | System side harmonic impedance estimation method and device based on mutual information data optimization |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10110976B2 (en) * | 2016-10-10 | 2018-10-23 | Utilidata, Inc. | Systems and methods for system measurements integrity determination |
CN111239466B (en) * | 2020-01-20 | 2021-02-02 | 同济大学 | Synchronous motor harmonic current detection method and system |
-
2023
- 2023-06-06 CN CN202310658384.5A patent/CN116388190B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014075581A1 (en) * | 2012-11-16 | 2014-05-22 | 天津大学 | Current-mode matching method for non-invasive power load monitoring and decomposition |
WO2014206007A1 (en) * | 2013-06-26 | 2014-12-31 | 国家电网公司 | Wind power plant harmonic suppression method, wind power generator sets, and wind power station |
CN103412242A (en) * | 2013-08-01 | 2013-11-27 | 西南交通大学 | Method for locating harmonic source based on independent rapid component analysis and mutual information |
CN103424620A (en) * | 2013-08-08 | 2013-12-04 | 华南理工大学 | 10kV distribution network harmonic source identification method based on time sequence relevancy |
CN104502704A (en) * | 2015-01-20 | 2015-04-08 | 西南交通大学 | Harmonic contribution calculating method applicable to background harmonic voltage change |
CN106208041A (en) * | 2016-07-18 | 2016-12-07 | 西南交通大学 | A kind of Distribution Network Harmonics current estimation method based on piecewise constant independent component analysis |
CN108089059A (en) * | 2017-12-19 | 2018-05-29 | 西南交通大学 | A kind of harmonic impedance evaluation method based on sparse independent component analysis |
CN114660362A (en) * | 2022-05-19 | 2022-06-24 | 四川大学 | System side harmonic impedance estimation method and device based on mutual information data optimization |
Non-Patent Citations (1)
Title |
---|
基于互信息数据优选的***侧谐波阻抗估计;徐方维等;《电力工程技术》;第67-74页 * |
Also Published As
Publication number | Publication date |
---|---|
CN116388190A (en) | 2023-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sun et al. | Probabilistic peak load estimation in smart cities using smart meter data | |
Anderson et al. | Nonparametric estimation and inference about the overlap of two distributions | |
Bücher et al. | A note on bootstrap approximations for the empirical copula process | |
Zang et al. | Adaptive method for harmonic contribution assessment based on hierarchical K‐means clustering and Bayesian partial least squares regression | |
Yu et al. | Fractional-order modeling of lithium-ion batteries using additive noise assisted modeling and correlative information criterion | |
Ma et al. | Topology identification of distribution networks using a split-EM based data-driven approach | |
Li et al. | Classification of power quality disturbance based on S-transform and convolution neural network | |
Morzfeld et al. | Gaussian approximations in filters and smoothers for data assimilation | |
CN111537884A (en) | Method and device for acquiring service life data of power battery, computer equipment and medium | |
Warin | Variations on branching methods for non linear PDEs | |
CN115436812B (en) | Electrochemical model parameter identification method and system based on splicing factors | |
CN116388190B (en) | Distributed photovoltaic harmonic quantitative evaluation method, system, terminal equipment and medium | |
CN115688017A (en) | FRCMDE-based transformer core fault voiceprint diagnosis method and device | |
CN116505653A (en) | Transmission line monitoring system, method, device, computer equipment and storage medium | |
Yusuff et al. | Statistical method for identification of weak nodes in power system based on voltage magnitude deviation | |
Parpaei et al. | Rational random walk‐based optimal placement of phasor measurement units to enhance the initializing and guiding the optimization processes | |
He et al. | SOH estimation for lithium-ion batteries: an improved GPR optimization method based on the developed feature extraction | |
Tang et al. | Novelty method for the utility harmonic impedance estimation based on Gaussian mixed model | |
McMichael et al. | Optbayesexpt: sequential Bayesian experiment design for adaptive measurements | |
Li et al. | Effect of data enhancement on state-of-charge estimation of lithium-ion battery based on deep learning methods | |
CN106897328A (en) | A kind of image search method and device | |
CN110703038B (en) | Harmonic impedance estimation method suitable for fan access power distribution network | |
JP2008152619A (en) | Data processor and data processing program | |
CN116799879A (en) | AC/DC hybrid power grid state estimation method, device and computer equipment | |
CN108198173B (en) | A kind of online test method, device and the terminal device in distress in concrete region |
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 |