CN106330095B - Collect system fault locating method in a kind of large-sized photovoltaic power station - Google Patents
Collect system fault locating method in a kind of large-sized photovoltaic power station Download PDFInfo
- Publication number
- CN106330095B CN106330095B CN201610857465.8A CN201610857465A CN106330095B CN 106330095 B CN106330095 B CN 106330095B CN 201610857465 A CN201610857465 A CN 201610857465A CN 106330095 B CN106330095 B CN 106330095B
- Authority
- CN
- China
- Prior art keywords
- node
- photovoltaic power
- fault
- neg
- power station
- 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- 239000013598 vector Substances 0.000 claims abstract description 34
- 238000002347 injection Methods 0.000 claims abstract description 12
- 239000007924 injection Substances 0.000 claims abstract description 12
- 239000000243 solution Substances 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 27
- 238000005259 measurement Methods 0.000 claims description 23
- 238000009826 distribution Methods 0.000 claims description 19
- 238000010248 power generation Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- 238000005457 optimization Methods 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 2
- 238000007906 compression Methods 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 230000001629 suppression Effects 0.000 claims description 2
- 230000007704 transition Effects 0.000 abstract description 8
- 238000000819 phase cycle Methods 0.000 abstract 2
- 238000004891 communication Methods 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02S—GENERATION OF ELECTRIC POWER BY CONVERSION OF INFRARED RADIATION, VISIBLE LIGHT OR ULTRAVIOLET LIGHT, e.g. USING PHOTOVOLTAIC [PV] MODULES
- H02S50/00—Monitoring or testing of PV systems, e.g. load balancing or fault identification
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
Abstract
The invention belongs to new energy grid-connected power technical field, more particularly to system fault locating method is collected in a kind of large-sized photovoltaic power station, the control used based on inverter in actual photovoltaic plant, using photovoltaic generation unit do not export negative-sequence current under the conditions of asymmetrical voltage the characteristics of, node negative sequence voltage equation and compressive sensing theory corresponding to sparse measuring point are combined after photovoltaic plant collects the interior generation unbalanced fault of system, and fault location is carried out by improved Bayes's compressed sensing restructing algorithm solution node negative phase-sequence Injection Current vector；What greatest member corresponding node occurrence number was most in the long interior nodes negative phase-sequence Injection Current amplitude vector of data window is regarded as malfunctioning node.The present invention only needs the asynchronous metrical information of sparse measuring point, and is not influenceed by fault type and transition resistance, and anti-noise ability is strong, and precision is high, disclosure satisfy that practical application request in large-sized photovoltaic power station.
Description
Technical Field
The invention belongs to the technical field of new energy grid-connected power generation, and particularly relates to a fault positioning method in a collecting system of a large photovoltaic power station.
Background
The large-scale centralized photovoltaic power station is formed by collecting a large number of photovoltaic power generation units, and overhead wires and buried cables coexist in a collecting system. The collecting system is not provided with an automatic reclosing device generally, and faults need to be manually checked after the faults trip. The photovoltaic array area is large, the degree of the appearance area is small, line patrol is difficult, the recovery power generation time of a disconnected photovoltaic system is long, and the power generation efficiency is lost. Therefore, the method is beneficial to improving the power generation efficiency of the photovoltaic power station by accurately and quickly positioning the fault position in the collecting system.
At present, fault location methods in power systems mainly include an impedance method, a traveling wave method, an S signal injection method, and a wide area communication method. The impedance method is easily influenced by factors such as power supply parameters and load parameters, and when more branches in the network have complex structures, false fault points are easy to appear. The traveling wave method is difficult to realize the accurate synchronous acquisition of the multi-end traveling wave signals under the condition of numerous complicated branches of the network. The injection signal intensity of the S signal injection method is limited by PT capacity, and the positioning accuracy is greatly influenced by factors such as wire distributed capacitance, grounding resistance and the like. The positioning accuracy of the method based on the single-point measurement information is greatly influenced by a complex network structure and the feeding fault current of the photovoltaic system.
The wide area communication method utilizes fault information of a plurality of measuring points in the power distribution network, and determines a fault section according to the difference of fault characteristics of the fault section and a non-fault section. In the existing fault positioning method based on multipoint measurement information, the method based on information interaction among measurement points needs synchronous measurement information and has high communication requirements. Most methods based on voltage sag can not avoid the cycle traversal process of possible fault points of the whole network, and when the network scale is large, the operation process is very long and the positioning accuracy is poor. When the access of the distributed power supply is considered, the distributed power supply is simply equivalent to a constant current source or a controllable source under the system frequency by the conventional method, the distributed and centralized photovoltaic grid-connected power generation systems are obviously different in topology and control, so that the fault current characteristics are greatly different, the conventional equivalent method is difficult to apply, the photovoltaic power supply permeability in the collection system reaches 100%, and the control characteristic of a photovoltaic inverter cannot be ignored. Therefore, it is necessary to research a fault location method suitable for the control characteristics of the inverter in the large photovoltaic power station.
Disclosure of Invention
In order to solve the problems, the invention provides a method for positioning faults in a collecting system of a large photovoltaic power station, which is characterized by comprising the following steps of:
step 1, analyzing to obtain a fault current characteristic rule in a collection system of a photovoltaic power station by combining the photovoltaic power station controlled by an inverter;
step 2, offline line connection and parameters in a photovoltaic power station collection system form a node negative sequence impedance matrix, selecting a line where a sparse measurement node is located, and taking absolute values of all elements to form a perception matrix;
step 3, after asymmetric faults occur in the photovoltaic power station collection system, combining the characteristic that the fault current of the photovoltaic system in the photovoltaic power station collection system has no negative sequence component, calculating by using the voltage information of the measured nodes after the faults to obtain a measured point negative sequence voltage amplitude vector, forming an underdetermined equation set for solving the sparse node negative sequence injection current amplitude vector by using a sensing matrix, and solving by using an improved Bayesian compression sensing reconstruction algorithm to obtain a sparse solution;
and 4, counting the occurrence times of the node corresponding to the maximum element in the node negative sequence injection current amplitude vector in the data window length, and regarding the maximum as a fault node.
The specific process of the step 1 is as follows: the method is characterized in that the characteristic analysis of fault current in a photovoltaic power station collecting system is combined with the control of a photovoltaic power generation system in the photovoltaic power station collecting system, the suppression of negative sequence component in inverter grid-connected current is taken as a control target, and reference current under a double-synchronous rotating coordinate is expressed as follows:
in the formula (I), the compound is shown in the specification,andthe positive sequence components of the d and q axes of the grid-connected point voltage,andis a grid-connected reference current d and q axis positive sequence component,andnegative sequence components of d and q axes of grid-connected reference current are obtained;for a given active power reference value,a given reactive power reference value;
obtaining the positive sequence component amplitude of the grid-connected current of the photovoltaic inverter:wherein the content of the first and second substances,is the amplitude of the positive sequence component of the grid-connected current,for the magnitude of the positive sequence component of the grid-connected point voltage,
according to the apparent power actually provided by the photovoltaic power generation system during the fault, the amplitude of the grid-connected current positive sequence component can be written as follows:wherein gamma is a positive sequence voltage drop coefficient, E m The grid-connected point voltage amplitude before the fault; p 0 ' active Power direct Current component, Q, supplied to photovoltaic Power generating Unit during Fault 0 ' is a reactive power direct current component provided by the photovoltaic power generation unit during a fault;
obtaining a three-phase fault current expression:
in the formula I a 、I b And I c Three-phase fault currents respectively; intermediate variables
In the step 2, the perception matrix calculation method comprises the following steps:
node negative sequence impedance array Z formed according to line topology and parameters in photovoltaic power station collection system N×N ，
In the formula, N is the number of network nodes in the photovoltaic power station collection system,is the negative sequence mutual impedance between node k and node l, when k = lIs a negative sequence self-impedance.
Selecting a row corresponding to a measuring point from a node negative sequence impedance array according to the distribution position of the measuring points in the photovoltaic power station collection system, and taking absolute values of all elements to obtain a sensing matrix Z neg ，
In the formula, M is the installation number of the measuring points; i.e. i j The node number j =1,2, … and M corresponding to the j measuring point is shown;represents a measurement point i j And the negative-sequence trans-impedance modulus value with node l.
The underdetermined equation set in the step 3 is as follows:
the above formula is abbreviated as V neg ＝Z neg ·I neg ；V neg Is a negative sequence voltage amplitude vector of Mx 1 dimension measuring pointIn which the elementsRepresents a measurement point i j The node is a negative sequence voltage module value; z neg Is a M × N dimensional perceptual matrix, I neg Injecting a current amplitude vector for a node negative sequence to be solved in Nx 1 dimension;representing the negative sequence current injected into the negative sequence network by the fault point.
The improved Bayes compressed sensing reconstruction algorithm in the step 3 is based on a compressed sensing theory and a sparse Bayes learning SBL algorithm;
the compressed sensing theoretical model y = phi theta + e, y is observation data of M × 1 dimension, phi is a sensing matrix of M × N dimension, theta is a sparse vector to be reconstructed of N × 1 dimension, and e is a obedient N (0, sigma) of N × 1 dimension 2 ) White Gaussian noise of (1), wherein M<&N; the compressed sensing optimization reconstruction algorithm utilizes the sparsity of theta, and solves the underdetermined equation set to recover the sparse theta with high probability and high precision through a small amount of measuring point data; i is neg Corresponding to the sparse vector theta, Z to be reconstructed neg Corresponding to the sensing matrix phi, V neg Corresponds to the observed data y.
The SBL algorithm assumes that each element of θ obeys a gaussian distribution:α i is unknown hyper-parameter, is the inverse variance of the gaussian distribution, and obtains a compressed sensing gaussian likelihood model:in the SBL algorithm, a Relevance Vector Machine (RVM) layering prior is used as a sparse prior of theta, and the sparse prior and a Gaussian likelihood model are conjugated, so that the calculation is convenient. Defining independent variables of theta in the sparse vector and following Gaussian distribution of zero mean:go toStep, assuming that the hyper-parameter α obeys the gamma prior distribution:where a, b are hyper-parameters in the parametric model. Obtaining a prior distribution form of theta through marginalizing a hyper-parameter alpha:when theta is i As a result of the observation data,is a likelihood function, gamma distribution gamma (alpha) i | a, b) is the inverse variance α from the Gaussian distribution i Is conjugated. The theta prior distribution now satisfies the t-distribution. When the parameters a and b are reasonably selected, the t-distribution at theta can be obtained i The place where =0 is sharp, therefore, the prior distribution of theta ensures that most variables of theta are 0, namely, ensures sparsity in the reconstruction process. Likewise, a gamma prior p (β | c, d) = Γ (β | a, b) is introduced, where β = σ -2 For the inverse noise variance, c and d are the hyper-parameters in the parametric model. The RVM mechanism is to obtain the hyper-parameter alpha and the noise inverse variance sigma of the sparse vector from the learning iteration process -2 . A, b, c, d are all generally 0, in which case the hyperparameters α and σ -2 A uniform distribution is obeyed. According to the prior distribution and the observation data y, the following posterior distribution can be obtained by a Bayesian formula: p (theta, alpha, sigma) 2 |y)＝p(θ|y,α,σ 2 )p(α,σ 2 Y), which involves a number of matrix inversion operations. In order to avoid huge calculation amount, a rapid RVM algorithm provided by taping and Faul is adopted, point estimation of a sparse vector theta is obtained through posterior estimation, and theta sparsity can be automatically ensured in an optimization solving process without parameter control.
The specific process of the step 4 is
1) Firstly, a node negative sequence impedance matrix is generated in an off-line manner according to structural parameters in a photovoltaic power station collection system, and Z is formed according to measuring point distribution neg ；
2) Reading voltage data of a measuring point after a fault, and calculating V of a continuous window length neg ；
3) Will Z neg V of starting time in window length neg Sending the parameters into a Bayesian compressed sensing reconstruction algorithm for calculation to obtain I neg ；
4) Continuously reading V at the next moment neg Repeating the step 3) until the V of the whole window length is read neg ；
5) And counting the occurrence frequency of the node corresponding to the element with the maximum median of the reconstructed vectors at each moment, wherein the node with the maximum occurrence frequency is regarded as a fault node.
In the step 4, the number of non-zero elements in the reconstructed vector theoretically corresponds to the number of fault nodes. However, due to the influence of the actual performance of the algorithm, measurement, transmission noise and the like, it is difficult to ensure that the number of non-zero elements in the reconstructed vector corresponds to the number of fault nodes. In order to reduce the possibility of misjudgment, the method observes the reconstruction result of a section of continuous window length, and counts the occurrence frequency of the node corresponding to the element with the maximum median of the reconstruction vector at each moment, wherein the node with the maximum occurrence frequency is regarded as a fault node.
Advantageous effects
(1) The influence of photovoltaic output conditions is avoided;
(2) Measuring points are sparse, measuring information does not need to be synchronous, and data used for positioning is short;
(3) The device is not affected by fault types and transition resistance, and has strong noise resistance.
Drawings
FIG. 1 is a flow chart of a large photovoltaic power station collection system fault location method of the present invention
FIG. 2 is a schematic view of a topology of a collection system of a large photovoltaic power station
FIG. 3 is a schematic view of a negative sequence network of the convergence system
FIG. 4 is a flowchart of a fault location algorithm of the present invention
I in FIG. 5 when an ABG-10 Ω fault occurs at node 15 neg Reconstructing the results
I at BC-25 Ω failure at node 40 of FIG. 6 neg Reconstructing the results
I at node 8 of FIG. 7 in the event of a CA-10 Ω failure neg Reconstructing the results
Detailed Description
The present invention will be described in detail below with reference to the accompanying drawings. The invention provides a method for positioning faults in a collecting system of a large photovoltaic power station, and fig. 1 is a flow chart of the method for positioning the faults in the collecting system of the large photovoltaic power station, which comprises the following steps:
step 1, analyzing and obtaining a fault current characteristic rule in a collecting system by combining control adopted by an inverter in a photovoltaic power station;
step 2, collecting line connection and parameter offline in the photovoltaic power station to form a node negative sequence impedance matrix, selecting a line where a sparse measurement node is located, and taking absolute values of all elements to form a sensing matrix;
step 3, after asymmetric faults occur in the collection system, combining the characteristic that fault current of a photovoltaic system in the station has no negative sequence component, calculating by using voltage information of measured nodes after the faults to obtain a voltage amplitude vector of the negative sequence of the measured points, forming an underdetermined equation set for solving sparse node negative sequence injection current amplitude vectors by using a sensing matrix, and solving by using an improved Bayesian compressed sensing reconstruction algorithm to obtain a sparse solution;
and 4, counting the occurrence times of the node corresponding to the maximum element in the node negative sequence injection current amplitude vector in the data window length, and regarding the maximum as a fault node.
Fig. 2 is a topological structure diagram of access of a single collection station in a collection system of a large photovoltaic power station, wherein a three-winding transformer in the diagram represents a photovoltaic power generation unit with a rated capacity of 1MW, a plurality of photovoltaic power generation units are accessed to a collection cable, and a plurality of collection cables are accessed to the collection station in a photovoltaic array area and then accessed to a booster station through an overhead line for centralized grid connection.
Based on control adopted by an inverter in an actual photovoltaic power station, theoretical analysis is carried out on the characteristics of fault current in a collecting system, and a photovoltaic power generation unit for inhibiting negative sequence current control does not output negative sequence current under the condition of asymmetric voltage. When an asymmetric fault occurs at a node in the trunking system, only one negative-sequence power supply is available at the faulty node in the entire network, as shown in fig. 3. The fault point negative sequence power supply can be regarded as a negative sequence current source, and the negative sequence voltage of each node can be obtained by the product of the node negative sequence impedance array and a sparse node negative sequence injection current vector (containing only one non-zero element):
in the formula (I), the compound is shown in the specification,and the negative sequence current injected into the negative sequence network by the fault point is shown, and N is the number of nodes of the network of the in-station collection system.
Considering that voltage measuring devices are arranged on M (M < N) nodes in the collection system, extracting corresponding M rows in the above formula to obtain a node negative sequence voltage equation formed by rows with known node negative sequence voltages. In order to avoid phase calculation and synchronization requirements on measurement information, both sides of the equation are in the form of amplitude values:
wherein M is the number of measuring point installations. i.e. i k The node numbers (k =1,2, …, M) corresponding to the kth measurement point are shown. The above formula is abbreviated as V neg ＝Z neg ·I neg ，V neg Is a negative sequence voltage amplitude vector of Mx 1 dimension measuring point, Z neg Is a M × N dimensional perceptual matrix, I neg And injecting a current amplitude vector for the node negative sequence to be solved in the dimension of Nx 1.Representing the negative sequence current injected into the negative sequence network by the fault point.
The above formula is an underdetermined equation set, infinite solutions exist, and the traditional least square method cannot be used for obtaining the sparsityAnd the fault point is easy to be judged wrongly or difficultly. And the compressed sensing theory can effectively avoid the problems. Compressed sensing utilizes a small amount of measuring point data, and an original signal which is sparse enough can be recovered by solving an underdetermined equation set through an optimized reconstruction algorithm. The correlation between the array and the column of the sensing matrix obtained by the node impedance matrix is strong and the original signal I to be reconstructed neg The method has the characteristic of extremely high sparsity, and the Bayes compressed sensing reconstruction algorithm based on the sparse Bayes learning algorithm is adopted, so that the algorithm is more suitable for processing the condition of strong correlation between columns in a sensing matrix than other optimization reconstruction algorithms, and the most sparse solution and stronger anti-noise capability are more easily obtained.
Fig. 4 is a flow chart of the fault location provided by the present invention, and the location steps after the fault occurs in the convergence system are as follows:
1) Firstly, a node negative sequence impedance matrix is generated in an off-line manner according to the structural parameters of a collecting system in a photovoltaic power station, and Z is formed according to the distribution of measuring points neg ；
2) Reading voltage data of a measuring point after a fault, and calculating V of a continuous window length neg ；
3) Will Z neg V of starting time in window length neg Sending the parameters into a Bayesian compressed sensing reconstruction algorithm for calculation to obtain I neg ；
4) Continuously reading V at the next moment neg Repeating the step 3) until the V of the whole window length is read neg ；
5) Theoretically, the number of non-zero elements in the reconstructed vector corresponds to the number of fault nodes. However, due to the influence of the actual performance and measurement of the algorithm, transmission noise and the like, it is difficult to ensure that only one non-zero element (considering a single fault) is contained in the reconstructed vector. Therefore, in order to reduce the possibility of erroneous judgment, the reconstruction result of a section of continuous window length is observed, the frequency of occurrence of the node corresponding to the element with the largest mean value of the reconstruction vector at each moment is counted, and the node with the largest frequency of occurrence is regarded as a fault node.
Taking the collection system given in fig. 3 as an example model, the number of measured points is 12. Asymmetric faults of different fault types and transition resistances are applied to all nodes in the system, the situation that measured data contain high noise is considered, and the positioning result is classified according to the distance between the calculation result of the method and a real fault node, as shown in table 1. Table 1 shows the positioning results of traversing all the nodes by the method, and the results in table 1 show that when a fault occurs in the convergence system, the positions of most of the fault nodes can be found more accurately by using the sparse measurement point asynchronous amplitude information. Due to the problem of the system structure, part of branches can not circulate negative sequence current, so that errors exist in positioning, and faults are positioned at adjacent nodes. However, such errors do not have much impact on fault clearance considering the actual wiring distance within the station.
TABLE 1
FIG. 5 is I of the power supply when the AB two-phase short circuit occurs at the node 15 and the ground fault occurs through the 10 omega transition resistor neg Reconstructing the result, wherein the x axis in the graph is the node number, the y axis is the time window length, and the z axis represents I neg The magnitude of the middle element. Reconstruction result I at any time in window length neg The 15 th element in the data is the largest, and the rest elements are all near 0, so that the fault node can be visually judged through a reconstruction result and is consistent with the position of the real fault node.
FIG. 6 shows the BC two-phase 25 Ω transition resistance short circuit fault at node 40, and the voltage measurement signal with 3% noise neg And reconstructing a result. At this time, the reconstruction result I at each time point in the window length neg In which a plurality of non-zero elements with distinct amplitudes are present, but at each instant I neg The 40 th element is much larger than the rest elements, so that the fault node can still be visually judged through the reconstruction result and is consistent with the real fault node.
FIG. 7 shows I when two phases of CA occur short circuit fault at node 8 via 10 Ω transition resistance neg The maximum element corresponding nodes at different moments in the reconstruction result are 7-11, and the final positioning result is the node 7 which is a real fault node adjacent node.
In summary, the method utilizes short asynchronous measurement information after the fault of the sparse measurement point, and can accurately find most fault nodes under the conditions of different fault types, transition resistance and high measurement noise. Due to the structural problem of the system, negative sequence current cannot flow in part of the branches, when faults occur at the last measuring point and the downstream node on the collection branch, the node corresponding to the maximum element in the reconstruction result is distributed between the measuring point and the tail end node of the branch, the final positioning result cannot exceed the range, even in a smaller range, the error cannot bring great influence to actual fault clearing, and the actual application requirements can be met.
In view of the analysis basis and conditions of the method, the method is not influenced by the operation conditions (such as output conditions, power factors and the like) of the photovoltaic system, only short asynchronous measurement information after the fault of the sparse measurement point is needed, and good fault positioning precision can be ensured under the conditions of different fault types, transition resistance and high measurement noise.
It should be noted that the above-mentioned embodiments are only preferred embodiments of the present invention, and should not be construed as limiting the scope of the present invention, and any minor variations and modifications of the present invention are within the scope of the present invention without departing from the spirit of the present invention.
Claims (5)
1. A method for positioning faults in a collection system of a large photovoltaic power station is characterized by comprising the following steps:
step 1, analyzing to obtain a fault current characteristic rule in a collection system of the photovoltaic power station by combining an inverter control mode actually adopted in the photovoltaic power station;
step 2, offline line connection and parameters in a photovoltaic power station collection system form a node negative sequence impedance matrix, selecting lines where sparse measurement points are located, and taking absolute values of all elements to form a sensing matrix;
step 3, after asymmetric faults occur in the photovoltaic power station collection system, combining the characteristic that the fault current of the photovoltaic system in the photovoltaic power station collection system has no negative sequence component, calculating by utilizing voltage information of the measured points after the faults to obtain a measured point negative sequence voltage amplitude vector, forming an underdetermined equation set for solving the sparse node negative sequence injection current amplitude vector by utilizing a sensing matrix, and solving by utilizing an improved Bayesian compression sensing reconstruction algorithm to obtain a sparse solution;
step 4, counting the occurrence frequency of the node corresponding to the maximum element in the node negative sequence injection current amplitude vector in the data window length, regarding the most nodes as fault nodes,
the improved Bayes compressed sensing reconstruction algorithm in the step 3 is based on a compressed sensing theory and a sparse Bayes learning SBL algorithm;
the compressed sensing theoretical model y = phi theta + e, y is observation data of M × 1 dimension, phi is a sensing matrix of M × N dimension, theta is a sparse vector to be reconstructed of N × 1 dimension, and e is a obedient N (0, sigma) of N × 1 dimension 2 ) White Gaussian noise of (1), wherein M<<, N; the compressed sensing optimization reconstruction algorithm utilizes the sparsity of theta, and solves the underdetermined equation set to recover the sparse theta with high probability and high precision through a small amount of measuring point data; I.C. A neg Corresponding to the sparse vector theta, Z to be reconstructed neg Corresponding to the sensing matrix phi, V neg Corresponds to the observed data y.
2. The method for locating faults in a collection system of large photovoltaic power stations as claimed in claim 1, wherein the specific process of step 1 is as follows: the method is characterized in that the characteristic analysis of fault current in a photovoltaic power station collecting system is combined with the control of a photovoltaic power generation system in the photovoltaic power station collecting system, the suppression of negative sequence component in inverter grid-connected current is taken as a control target, and reference current under a double-synchronous rotating coordinate is expressed as follows:
in the formula (I), the compound is shown in the specification,andthe positive sequence components of the grid-connected point voltage d and q axes,andis a grid-connected reference current d and q axis positive sequence component,andnegative sequence components of d and q axes of grid-connected reference current are obtained; p is 0 * For a given active power reference value, Q 0 * A given reactive power reference value;
obtaining the positive sequence component amplitude of the grid-connected current of the photovoltaic inverter:wherein the content of the first and second substances,is the amplitude of the positive sequence component of the grid-connected current, for the positive sequence component amplitude of the grid-connected point voltage,
according to the apparent power actually provided by the photovoltaic power generation system during the fault, the amplitude of the positive sequence component of the grid-connected currentWherein gamma is a positive sequence voltage drop coefficient, E m The grid-connected point voltage amplitude before the fault; p is 0 ' active Power DC component, Q, supplied to photovoltaic Power Generation Unit during Fault 0 ' is a reactive power direct current component provided by the photovoltaic power generation unit during a fault;
obtaining a three-phase fault current expression:
in the formula I a 、I b And I c Three-phase fault currents respectively; intermediate variables
3. The method for locating faults in a collection system of a large photovoltaic power station as claimed in claim 1, wherein in the step 2, the sensing matrix calculation method comprises the following steps:
node negative sequence impedance array Z formed according to line topology and parameters in photovoltaic power station collection system N×N ，
In the formula, N is the number of network nodes in the photovoltaic power station collection system,is the negative sequence mutual impedance between node k and node l, when k = lIs a negative-sequence self-impedance,
selecting a row corresponding to a measuring point from a node negative sequence impedance array according to the distribution position of the measuring points in the photovoltaic power station collection system, and taking absolute values of all elements to obtain a sensing matrix Z neg ，
In the formula, M is the installation number of the measuring points; i.e. i j The node numbers j =1,2, … and M corresponding to the j measuring point are shown; represents the measurement point i j And the negative-sequence trans-impedance modulus value with node l.
4. The method for locating faults in a collection system of large photovoltaic power stations as claimed in claim 1, wherein the underdetermined equations in step 3 are:
the above formula is abbreviated as V neg ＝Z neg ·I neg ；V neg Is a negative sequence voltage amplitude vector of Mx 1 dimension measuring point, wherein the elementsRepresents a measurement point i j The node is a negative sequence voltage module value; z neg Is a M × N dimensional perceptual matrix, I neg Injecting a current amplitude vector for a node negative sequence to be solved in dimension Nx 1;representing the negative sequence current injected into the negative sequence network by the fault point.
5. The method for locating faults in a collection system of large photovoltaic power stations as claimed in claim 1, wherein the specific process of step 4 is
1)Firstly, a node negative sequence impedance matrix is generated in an off-line manner according to structural parameters in a photovoltaic power station collection system, and Z is formed according to measuring point distribution neg ；
2) Reading voltage data of a measuring point after a fault, and calculating V of a continuous window length neg ；
3) Will Z neg V of starting time in window length neg Sending the parameters into a Bayesian compressed sensing reconstruction algorithm for calculation to obtain I neg ；
4) Continuously reading V at the next moment neg Repeating the step 3) until the V of the whole window length is read neg ；
5) And counting the occurrence frequency of the node corresponding to the element with the maximum median of the reconstructed vectors at each moment, wherein the node with the maximum occurrence frequency is regarded as a fault node.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610857465.8A CN106330095B (en) | 2016-09-27 | 2016-09-27 | Collect system fault locating method in a kind of large-sized photovoltaic power station |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610857465.8A CN106330095B (en) | 2016-09-27 | 2016-09-27 | Collect system fault locating method in a kind of large-sized photovoltaic power station |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106330095A CN106330095A (en) | 2017-01-11 |
| CN106330095B true CN106330095B (en) | 2018-04-10 |
Family
ID=57820151
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610857465.8A Active CN106330095B (en) | 2016-09-27 | 2016-09-27 | Collect system fault locating method in a kind of large-sized photovoltaic power station |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106330095B (en) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108764997A (en) * | 2018-05-25 | 2018-11-06 | 北京京东金融科技控股有限公司 | Take in recognition methods, device and computer readable storage medium |
| CN109142986B (en) * | 2018-10-10 | 2020-09-29 | 华北电力大学 | Fault positioning method for alternating current distribution line |
| CN110045226A (en) * | 2019-03-18 | 2019-07-23 | 昆明理工大学 | A kind of compressed sensing based electrical power distribution network fault location method |
| CN110389282A (en) * | 2019-07-17 | 2019-10-29 | 华北电力大学 | Accurate fault distance-finding method suitable for more power distribution circuits |
| CN110930263B (en) * | 2019-11-15 | 2020-12-29 | 广东电网有限责任公司 | Medium-voltage distribution network short-circuit current calculation method containing photovoltaic power supply and induction motor based on black hole particle swarm algorithm |
| CN111337793B (en) * | 2020-04-22 | 2021-05-25 | 华北电力大学 | Sparse measurement point configuration method suitable for power distribution network compressed sensing fault positioning technology |
| CN113933644A (en) * | 2020-07-14 | 2022-01-14 | 上海华建电力设备股份有限公司 | Power transmission line asymmetric fault positioning method free of influence of line parameters |
| CN114039546A (en) * | 2021-11-18 | 2022-02-11 | 河海大学常州校区 | Photovoltaic array fault positioning method and device based on photovoltaic module I-V curve reconstruction |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6231900B2 (en) * | 2014-02-17 | 2017-11-15 | 株式会社ダイヘン | Photovoltaic power generation facility for grid connection |
| CN204168236U (en) * | 2014-11-25 | 2015-02-18 | 江苏蓝天光伏科技有限公司 | A kind of experimental rig for photovoltaic power generation quantity measuring and calculating |
-
2016
- 2016-09-27 CN CN201610857465.8A patent/CN106330095B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106330095A (en) | 2017-01-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106330095B (en) | Collect system fault locating method in a kind of large-sized photovoltaic power station | |
| Belagoune et al. | Deep learning through LSTM classification and regression for transmission line fault detection, diagnosis and location in large-scale multi-machine power systems | |
| Koley et al. | Detection and classification of faults on six phase transmission line using ANN | |
| CN106771868A (en) | A kind of flexible direct current distribution system fault locating method based on traveling wave | |
| Kumar et al. | Fault classification of phase to phase fault in six phase transmission line using Haar wavelet and ANN | |
| CN110348114B (en) | Non-precise fault identification method for power grid completeness state information reconstruction | |
| Xiong et al. | A data-driven approach for fault time determination and fault area location using random matrix theory | |
| CN107329045A (en) | Distribution Network Failure least square location algorithm | |
| CN107462810A (en) | A kind of fault section location method suitable for active power distribution network | |
| CN111654392A (en) | Low-voltage distribution network topology identification method and system based on mutual information | |
| CN106019043B (en) | Power grid fault diagnosis method based on fuzzy integral multi-source information fusion theory | |
| Zhang et al. | Detection of single-phase-to-ground faults in distribution networks based on Gramian Angular Field and Improved Convolutional Neural Networks | |
| El-Tawab et al. | Self-healing of active distribution networks by accurate fault detection, classification, and location | |
| CN113009279B (en) | Neo4 j-based power distribution network voltage sag fault positioning and visualization system | |
| Shu et al. | Fault phase selection and distance location based on ANN and S-transform for transmission line in triangle network | |
| Siddique et al. | Fault classification and location of a PMU-equipped active distribution network using deep convolution neural network (CNN) | |
| Ali et al. | The Transmission line fault detection and identification using Artificial Neural Network | |
| Chatterjee et al. | Cross-correlation based single-phase fault location scheme for transmission lines using one-end voltage | |
| Peng et al. | Fault Location Algorithm for Multi-Terminal Transmission Lines of Distribution System Based on Loop-Analysis | |
| Yu et al. | Research on Asymmetric Fault Location of Wind Farm Collection System Based on Compressed Sensing | |
| Liang | Single-phase Grounding Fault Location Method for Distribution Network Based on Convolutional Neural Network | |
| Datta et al. | Line to ground and line to line fault analysis in ieee standard 9 BUS system | |
| Luo et al. | Intelligent Location Method with Limited Measurement Information for Multi-Branch Distribution Networks | |
| Hossain et al. | Distribution fault location using wide-area voltage magnitude measurements | |
| Phukan et al. | Fault diagnosis of a high voltage transmission line using waveform matching approach |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |