CN116388190A - 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 PDF

Info

Publication number
CN116388190A
CN116388190A CN202310658384.5A CN202310658384A CN116388190A CN 116388190 A CN116388190 A CN 116388190A CN 202310658384 A CN202310658384 A CN 202310658384A CN 116388190 A CN116388190 A CN 116388190A
Authority
CN
China
Prior art keywords
feeder line
feeder
harmonic
current
matrix
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
Application number
CN202310658384.5A
Other languages
Chinese (zh)
Other versions
CN116388190B (en
Inventor
孙媛媛
尹书林
许庆燊
丁磊
李亚辉
孙瑞泽
亓德民
单鹏博
李博文
孟广泽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University
Original Assignee
Shandong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shandong University filed Critical Shandong University
Priority to CN202310658384.5A priority Critical patent/CN116388190B/en
Publication of CN116388190A publication Critical patent/CN116388190A/en
Application granted granted Critical
Publication of CN116388190B publication Critical patent/CN116388190B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/01Arrangements for reducing harmonics or ripples
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381Dispersed generators
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10Power 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
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20The dispersed energy generation being of renewable origin
    • H02J2300/22The renewable source being solar energy
    • H02J2300/24The renewable source being solar energy of photovoltaic origin

Abstract

The invention 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 solutionA coupled 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 invention can accurately evaluate the harmonic responsibility of a plurality of feeder lines with higher harmonic current correlation under the distributed photovoltaic access.

Description

Distributed photovoltaic harmonic quantitative evaluation method, system, terminal equipment and medium
Technical Field
The invention 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 invention 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 results
Figure SMS_1
The 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 that
Figure SMS_2
In the feeder line with mutual information value greater than the set value, calculating the firstmFeeder current and +.>
Figure SMS_3
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 invention has the beneficial effects that:
(1) On the basis of traditional multi-feeder harmonic quantitative evaluation, the invention 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 invention 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 invention.
Drawings
FIG. 1 is a schematic diagram of a large-scale photovoltaic parallel-grid multi-feeder structure in an embodiment of the present invention;
fig. 2 is a flowchart of a method for quantitatively evaluating distributed photovoltaic harmonics in an embodiment of the present invention.
Detailed Description
It should be noted that the following detailed description is illustrative and is intended to provide further explanation of the present 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 example embodiments in accordance with 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
Figure SMS_4
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 impedance
Figure SMS_5
Andhsubequivalent harmonic current source->
Figure SMS_6
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 comprisinghInferior gradeEffective harmonic impedance->
Figure SMS_7
Andhsubequivalent harmonic current source->
Figure SMS_8
Is formed by connecting the two parts in parallel,irepresenting the first of the usersiAnd (5) a feeder line.
Feeder current matrix of this embodiment
Figure SMS_9
The method comprises the following steps:
Figure SMS_10
wherein,
Figure SMS_11
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:
Figure SMS_12
wherein,
Figure SMS_13
representing feederiAnd feeder linejA kind of electronic devicehMutual information value of the secondary current.
Figure SMS_14
Representing feederiA kind of electronic devicehThe entropy of the secondary current data is calculated by the following formula:
Figure SMS_15
Figure SMS_16
representing feederiA kind of electronic devicehProbability value of occurrence of a certain number of secondary current data.
Figure SMS_17
Representing feederiAnd feeder linejA kind of electronic devicehThe calculation formula of the joint entropy value of the secondary current data is as follows:
Figure SMS_18
Figure SMS_19
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 matrix
Figure SMS_20
The method comprises the following steps:
Figure SMS_21
wherein (1)>
Figure SMS_22
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:
Figure SMS_23
in the method, in the process of the invention,
Figure SMS_24
for observing signals, +.>
Figure SMS_25
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:
Figure SMS_26
in the method, in the process of the invention,
Figure SMS_27
is the mean signal of the observed signals.
Whitening the data, removing correlation between the metrology signals:
Figure SMS_28
in the method, in the process of the invention,
Figure SMS_30
for whitening matrix +.>
Figure SMS_32
,/>
Figure SMS_36
Is->
Figure SMS_31
Eigenvalues of matrix +.>
Figure SMS_34
Is that
Figure SMS_35
Characteristic value matrix of matrix, ">
Figure SMS_37
Matrix after data whitening, ++>
Figure SMS_29
And->
Figure SMS_33
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
Figure SMS_38
:
Figure SMS_39
In the method, in the process of the invention,Wfor an iteratively solved parameter matrix, wherein
Figure SMS_40
Is 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;
Figure SMS_41
wherein,
Figure SMS_42
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, obtainhThe harmonic quantitative evaluation result of subharmonic is:
Figure SMS_43
in the method, in the process of the invention,
Figure SMS_44
is->
Figure SMS_45
Is the first of (2)lColumn (1)kColumn harmonic quantitative evaluation results->
Figure SMS_46
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 be
Figure SMS_47
In the feeder line with mutual information value greater than the set value, the firstmFeeder current and +.>
Figure SMS_48
Correlation coefficient of->
Figure SMS_49
The method specifically comprises the following steps:
Figure SMS_50
wherein (1)>
Figure SMS_51
In the feeder line with mutual information value larger than the set value, the firstmFirst of the feeder linesi-individual current data; />
Figure SMS_52
In the feeder line with mutual information value larger than the set value, the firstmThe average value of the currents of the feeder lines; />
Figure SMS_53
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:
Figure SMS_54
wherein,
Figure SMS_55
is->
Figure SMS_56
Is a cumulative value of (a); />
Figure SMS_57
Is the current fluctuation coefficient; />
Figure SMS_58
In the feeder line with mutual information value larger than the set value, the firstmThe average value of the currents of the feeder lines; />
Figure SMS_59
In the feeder line with mutual information value larger than the set value, the firstmThe maximum current value of the strip feed line. />
Figure SMS_60
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 invention has been presented in conjunction with the drawings, it should be understood that it is not intended to limit the scope of the invention, but rather, it is intended to cover all modifications or variations within the scope of the invention as defined by the claims of the present invention.

Claims (10)

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.
2. A distributed photovoltaic harmonic quantification assessment method according to claim 1, wherein the method compriseshIn the harmonic quantitative evaluation result of subharmonic, the firstkColumn harmonic quantitative evaluation results
Figure QLYQS_1
The 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 that
Figure QLYQS_2
In the feeder line with mutual information value greater than the set value, calculating the firstmFeeder current and +.>
Figure QLYQS_3
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.
3. A distributed photovoltaic harmonic quantification assessment method according to claim 2, wherein the first calculation is performedmFeeder current of strip feeder
Figure QLYQS_4
Correlation coefficient of->
Figure QLYQS_5
The method specifically comprises the following steps:
Figure QLYQS_6
wherein,
Figure QLYQS_7
in the feeder line with mutual information value larger than the set value, the firstmFirst of the feeder linesi-individual current data; />
Figure QLYQS_8
In the feeder line with mutual information value larger than the set value, the firstmThe average value of the currents of the feeder lines; />
Figure QLYQS_9
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.
4. 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:
Figure QLYQS_10
wherein,
Figure QLYQS_11
is->
Figure QLYQS_12
Is a cumulative value of (a); />
Figure QLYQS_13
Is the current fluctuation coefficient; />
Figure QLYQS_14
In the feeder line with mutual information value larger than the set value, the firstmThe average value of the currents of the feeder lines; />
Figure QLYQS_15
In the feeder line with mutual information value larger than the set value, the firstmA current maximum value of the strip feeder; />
Figure QLYQS_16
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.
5. 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:
Figure QLYQS_17
wherein,
Figure QLYQS_18
representing feederiA kind of electronic devicehEntropy of subharmonic current; />
Figure QLYQS_19
Representing feederjA kind of electronic devicehEntropy of subharmonic current>
Figure QLYQS_20
Representing feederiAnd feeder linejA kind of electronic devicehA joint entropy value of the subharmonic current; />
Figure QLYQS_21
Representing feederiAnd feeder linejA kind of electronic devicehMutual information value of subharmonic current.
6. 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:
Figure QLYQS_22
wherein,
Figure QLYQS_23
is the firstlHarmonic quantitative evaluation results of the strip feeder lines; />
Figure QLYQS_24
Is a feeder linelRegression coefficients of (a); />
Figure QLYQS_25
Is that
Figure QLYQS_26
Is the first of (2)lColumn (S)/(S)>
Figure QLYQS_27
Is a separation matrix; />
Figure QLYQS_28
Represents n 1hA secondary common coupling point harmonic voltage matrix.
7. A distributed photovoltaic harmonic quantitative evaluation method as claimed in claim 6 wherein the regression coefficients for each feed line are solved using:
Figure QLYQS_29
wherein, Z 1 、Z 2 、…、Z n the feed lines 1,2 …,nregression coefficient of E is n* An identity matrix of 1,U s is a harmonic voltage on the system side.
8. 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.
9. 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-7.
10. 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-7.
CN202310658384.5A 2023-06-06 2023-06-06 Distributed photovoltaic harmonic quantitative evaluation method, system, terminal equipment and medium Active CN116388190B (en)

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 true CN116388190A (en) 2023-07-04
CN116388190B 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 (10)

* Cited by examiner, † Cited by third party
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
US20180100912A1 (en) * 2016-10-10 2018-04-12 Utilidata, Inc. Systems and methods for system measurements integrity determination
CN108089059A (en) * 2017-12-19 2018-05-29 西南交通大学 A kind of harmonic impedance evaluation method based on sparse independent component analysis
US20210223295A1 (en) * 2020-01-20 2021-07-22 Tongji University Method and system for detecting harmonic current in synchronous motors
CN114660362A (en) * 2022-05-19 2022-06-24 四川大学 System side harmonic impedance estimation method and device based on mutual information data optimization

Patent Citations (10)

* Cited by examiner, † Cited by third party
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
US20180100912A1 (en) * 2016-10-10 2018-04-12 Utilidata, Inc. Systems and methods for system measurements integrity determination
CN108089059A (en) * 2017-12-19 2018-05-29 西南交通大学 A kind of harmonic impedance evaluation method based on sparse independent component analysis
US20210223295A1 (en) * 2020-01-20 2021-07-22 Tongji University Method and system for detecting harmonic current in synchronous motors
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 (2)

* Cited by examiner, † Cited by third party
Title
PILIMERA K: ":Independent component analysis for harmonic source estimation from piecewise Constant parameter mixed measurments", 《SYSTEM THEORY》, pages 81 - 86 *
徐方维等: "基于互信息数据优选的系统侧谐波阻抗估计", 《电力工程技术》, pages 67 - 74 *

Also Published As

Publication number Publication date
CN116388190B (en) 2023-08-18

Similar Documents

Publication Publication Date Title
Bücher et al. A note on bootstrap approximations for the empirical copula process
Anderson et al. Nonparametric estimation and inference about the overlap of two distributions
Zang et al. Adaptive method for harmonic contribution assessment based on hierarchical K‐means clustering and Bayesian partial least squares regression
Dahale et al. Sparsity based approaches for distribution grid state estimation-a comparative study
CN112907062B (en) Power grid electric quantity prediction method, device, medium and terminal integrating temperature characteristics
Ma et al. Topology identification of distribution networks using a split-EM based data-driven approach
Morzfeld et al. Gaussian approximations in filters and smoothers for data assimilation
Warin Variations on branching methods for non linear PDEs
CN116388190B (en) Distributed photovoltaic harmonic quantitative evaluation method, system, terminal equipment and medium
Afrasiabi et al. Wide-area composite load parameter identification based on multi-residual deep neural network
Parpaei et al. Rational random walk‐based optimal placement of phasor measurement units to enhance the initializing and guiding the optimization processes
JP2008152619A (en) Data processor and data processing program
CN110703038B (en) Harmonic impedance estimation method suitable for fan access power distribution network
CN116879634A (en) Method, apparatus, device, medium and program product for monitoring line resistance of ground electrode
Yao et al. Efficient importance sampling for high-sigma yield analysis with adaptive online surrogate modeling
He et al. SOH estimation for lithium-ion batteries: An improved GPR optimization method based on the developed feature extraction
CN108198173B (en) A kind of online test method, device and the terminal device in distress in concrete region
CN106897328A (en) A kind of image search method and device
Acurio et al. Design and Implementation of a Machine Learning State Estimation Model for Unobservable Microgrids
CN110502715B (en) Click probability prediction method and device
Lee et al. Maternity length of stay modelling by gamma mixture regression with random effects
CN112381456A (en) Equipment development phase testability evaluation and system based on conflict data fusion
Alvarez-Esteban et al. Assessing when a sample is mostly normal
CN116796926B (en) Low-voltage side new energy resource aggregation capacity assessment method and device, medium and equipment
Lim et al. Bad data detection and estimation in high dimensional measurement data

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