Summary of the invention
After accessing DG for power distribution network, the decline of conventional failure locating scheme positional accuracy, low efficiency and power distribution network three
The asymmetric problem of phase parameter, the present invention propose a kind of based on the active power distribution network fault location for improving Karrenbauer transformation
New method.It the characteristics of feature and fault current of foundation active power distribution network parameter matrix, searches out and is suitable for the distribution that do not replace
The phase-model transformation method and fault location criterion of net three-phase network high-precision decoupling, realize the quick essence to the active power distribution network containing DG
True fault location.Technical solution is as follows:
A kind of active power distribution network Fault Locating Method based on improvement Karrenbauer transformation, including the following steps:
(1) control centre obtains shown in the form such as formula (1) of active power distribution network impedance matrix Z, Z, first uses standard
Karrenbauer transformation matrix does diagonalization processing to impedance matrix, obtains preliminary transformation matrix Y:
Y=K-1ZK (2)
To preliminary transformation matrix Y amendment Karrenbauer matrix K2Diagonalization processing is done, the diagonal of degree of precision is obtained
Change matrix λ:
λ=(K2)-1YK2=(KK2)-1Z(KK2) (3)
Wherein:
Enable S=KK2, S is to improve Karrenbauer transformation matrix;
(2) active power distribution network is divided into several sections, and the node of each section and its both ends is marked respectively
Number, the three-phase current at each two end node of section of protection device synchronous acquisition and progress at each two end node of section
Filtering processing recycles improvement Karrenbauer transformation matrix to convert to obtain α, β, 0 mould electric current;
(3) local data processing unit calculates its failure extra current after taking section both ends α mould electric currentWith quick
Fourier transformation obtains section both ends α mould failure extra current phase angle theta1(i)、θ1(j), calculation of sector α mould failure extra current phase
Angular difference:
Δθ1(i, j)=θ1(i)-θ1(j)
Gained Δ θ1(i, j) data are sent to control centre by communication, and control centre receives and stores Δ θ1(i, j), it is raw
At real-time system status information matrix;
(4) the break down threshold value δ of differentiation of control centre's section takes 7.5 °~8 °, judges first in system status information matrix
Whether element value is greater than threshold value, if certain element is lower than this threshold value, illustrates that the section does not break down;If certain element is greater than the threshold
Value, then differentiate that the section breaks down;
(5) control centre is detecting and is carrying out troubleshooting after positioning failure immediately, the protection to fault section both ends
Device sending action signal, isolated fault.
Compared with prior art, this programme proposes a kind of base on the basis of fully considering distributed generation resource access
In improve Karrenbauer transformation active power distribution network fault location new method, reach it is following the utility model has the advantages that
(1) the characteristics of being directed to distribution network line parameter unbalance, based on the phase-model transformation side for improving Karrenbauer transformation
The high-precision decoupling of dissymmetric network has may be implemented in method.The transform method calculating process is simple, calculation amount is smaller, decoupling precision
It is higher.
(2) on the basis of dissymmetric network high-precision decouples, active power distribution network section fault orientation criterion is derived.This
Criterion by transition Resistance Influence, does not meet distribution network parameters actual features, and active power distribution network may be implemented and position rapidly, without adding
Directional element is filled, there is higher reliability.
Specific embodiment
For the asymmetric actual conditions of power distribution network triphase parameter containing distributed generation resource, the invention proposes based on improvement
The active power distribution network fault location new method of Karrenbauer transformation.Power distribution network is first depending on because of line parameter circuit value caused by not replacing
It is asymmetric but have the characteristics that parameter impedance matrix central symmetry, to distribution network line impedance matrix progress standard
Karrenbauer is converted to increase the sparsity of matrix, provides premise for accurate decomposition;Further preliminary transformation matrix is done
Diagonalization is handled to obtain Karrenbauer correction matrix, is decoupled to obtain the achievable power distribution network three-phase high-precision that do not replace
Phase mode transformation matrix, i.e., improvement Karrenbauer matrix.On the basis of three-phase network high-precision decouples, pass through theory point
α mould failure extra current component phase angle difference after active power distribution network three-phase decouples is analysed, deriving not is had by transition Resistance Influence
The criterion of source power distribution network section fault positioning, is realized to the accurate fault location containing distributed power distribution network.The scheme packet
Include following steps:
1, the improvement Karrenbauer transformation of dissymmetric network accurate decomposition is realized
Actual low and medium voltage distribution network three-phase conducting wire is often in horizontally or vertically to arrange, the no longer complete phase of distance between three-phase
Deng this also results in the asymmetry of system line parameter.Analysis finds that the spread pattern of low and medium voltage distribution network conducting wire makes route
Impedance matrix is centrosymmetric the form of battle array:
Because of line parameter circuit value asymmetry caused by not replacing but have the characteristics that parameter impedance matrix central symmetry for power distribution network,
Impedance matrix is subjected to standard K arrenbauer transformation to increase the sparsity of matrix, provides premise for accurate decomposition;Herein
On the basis of, further diagonalization processing is done to obtain correction matrix to preliminary transformation matrix, and then obtain can be achieved not replace to match
The phase mode transformation matrix of power grid three-phase high-precision decoupling, i.e. improvement Karrenbauer matrix.
(1) standard K arrenbauer is converted
Known Karrenbauer transformation matrix is as follows:
For mathematically, the Uncoupled procedure of impedance matrix, that is, seek phase mode transformation matrix S and realize Z pairs of impedance matrix
The process of angling is S with formulae express-1ZS=λ.Impedance matrix Z is carried out using standard K arrenbauer transformation matrix preliminary
Transformation can obtain:
As it can be seen that for the low and medium voltage distribution network line impedance matrix as shown in formula (1), passing through for form
Karrenbauer transformation, the degree of coupling between α mold component and β mold component obtains to weaken to a certain extent.This degree of coupling
Reduction will further decrease the complexity of Z decoupling, be advantageously implemented the high-precision decoupling to impedance matrix Z.
(2) Karrenbauer transformation is improved
By formula (3) it is found that the decoupling to finally realize matrix Z, need to only realize matrix Y and decouple, that is to say, that warp
It crosses after Karrenbauer tentatively converts, the decoupling problem of Z can be further converted to the decoupling problem of Y: seek phase mode transformation matrix
K2Realize the diagonalization of matrix Y.Mathematical expression is as follows:
(K2)-1YK2=λ (5)
To modulus transformation matrix K2, start with first from solution eigenvalue λ:
Characteristic value can be acquired by formula (6):
In this way, the phase mode transformation matrix of Y are as follows:
In summary transformation twice, available:
Therefore, impedance matrix Z has just obtained the diagonalization of degree of precision by above-mentioned transformation.New phase mode transformation matrix are
S=KK2, K2It can be regarded as the amendment to original Karrenbauer transformation matrix, referred to as amendment Karrenbauer matrix.
New phase mode transformation matrix S=KK2Referred to as improve Karrenbauer transformation matrix.
2, based on the power distribution network section fault orientation criterion for improving Karrenbauer transformation
(1) when internal fault α mould network failure extra current phase angle difference
By principle of stacking it is found that the equivalent sequence diagrams of α mould can be decomposed into normal condition and failure additivity two parts, point
Not as shown in Fig. 1 (a), (b), wherein the reference direction and normal current of regulation failure extra currentDirection it is identical.
It can be with calculation of sector both ends failure extra current according to Fig. 1 (b) are as follows:
In formula:
Because the parameter of route is equally distributed, lumped parameter along routeIt can be expressed as Form, Z indicate the impedance of route unit length.Substituting into formula (11) final finishing can obtain:
It can deduce the phase angle difference of failure extra current are as follows:
According to formula (13) it is found that for α lay wire network, when troubles inside the sample space occurs, the phase angle of failure extra current
Difference is 180 °.
(2) when external fault α mould network failure extra current phase angle difference
When breaking down outside section, the section first, last both ends are identical by electric current, so α mould failure is additional at this time
The phase angle difference of electric current is 0 °.
Comprehensive front analysis is it is found that the phase angle difference of α mould failure extra current is 180 ° when intra-segment breaks down;
When breaking down outside section, the phase angle difference of α mould failure extra current is 0 °.
In view of measurement error that may be present and the calculating error of phase-model transformation in practice, to be set for failure criterion
One threshold value δ just assert that the section breaks down when fault component current phase angle difference is greater than the threshold value.Meter and Current Mutual Inductance
The angle error of device, general load are selected by 10% error curve, are thought of as 7 °;Calculating error can consider by 0.5 °~1 °.Institute
7.5 °~8 ° are generally taken with the threshold value of criterion.
Therefore, shown in active power distribution network fault location criterion such as formula (14):
3, influence of the transition resistance to fault location
The metallicity failure that power distribution network occurs is seldom, and short dot has transition resistance mostly.For the network after phase-model transformation
For, what transition resistance influenced is additional impedance.So transition resistance for α mould, can only change the value of additional impedance.By
Formula (13) is it is found that Δ θ1(i, j)=arg (- k2/k1), it is unrelated with transition resistance.So being based on α mould failure extra current phase angle
The failure criterion of difference is not influenced by transition resistance.
In order to verify the validity of proposed fault location strategy herein, asymmetrical three-phase active power distribution network is chosen herein and is calculated
Example, as shown in Figure 3.Corresponding simulation model has been built in Matlab/Simulink.
Line impedance matrix are as follows:
(1) impedance matrix decouples
Impedance matrix is decoupled using symmetrical component method, as a result as follows:
Impedance matrix is decoupled using Karrenbauer transformation is improved, as a result as follows:
Compare above two decoupling as a result, it has been found that: for parameter unbalance matrix, symmetrical component method weakens coupling journey
The ability of degree is smaller, and mentions improve the impedance matrix decoupling that Karrenbauer transformation is then able to achieve degree of precision herein.
(2) fault location emulates
Table 1 is mentions the Fault Locating Method based on improvement Karrenbauer transformation according to this paper, in different sections, no
With the simulation scenarios in the case of transition resistance.Comparison whether exists it is found that when different types of faults occur for different sections
Transition resistance, mentioned Fault Locating Method can realize accurate fault section location herein.
Table 1 is based on the fault location result for improving Karrenbauer transformation
When transition resistance value difference, when section 3 breaks down, the section is based on the α for improving Karrenbauer transformation
The value of modal transformation based on fault component current phase angle difference is as shown in table 2.As seen from table, when transition resistance difference, fault component phase angle difference
It is worth identical, the presence of transition resistance and value have no influence to mentioned Fault Locating Method herein, with front theoretical analysis result
Unanimously.
Fault component current phase angle difference (unit/0) when 2 difference transition resistance of table