CN105808967B - Exciter saturation coefficient computational methods based on secant iterative method - Google Patents
Exciter saturation coefficient computational methods based on secant iterative method Download PDFInfo
- Publication number
- CN105808967B CN105808967B CN201610208149.8A CN201610208149A CN105808967B CN 105808967 B CN105808967 B CN 105808967B CN 201610208149 A CN201610208149 A CN 201610208149A CN 105808967 B CN105808967 B CN 105808967B
- Authority
- CN
- China
- Prior art keywords
- exciter
- value
- saturation coefficient
- secant
- saturation
- 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 15
- 238000000205 computational method Methods 0.000 title abstract 2
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 230000005284 excitation Effects 0.000 claims description 11
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 230000000750 progressive effect Effects 0.000 abstract 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
The present invention provides a kind of exciter saturation coefficient computational methods based on secant iterative method, by the progressive fitting exciter no-load curve of secant iterative method, to reach the purpose for calculating saturation coefficient.This method avoids error caused by manually recognizing, and it is suitable for the test data containing interference.The present invention has the result of calculation more more precisely consistent than artificial interpretation, the problem of avoiding different artificial sentence read result difference, suitable for the exciter saturation coefficient calculating based on exciter no-load characteristic.
Description
Technical Field
The invention relates to the technical field of generator excitation systems, in particular to a method for calculating saturation coefficients of an exciter based on a secant iteration method.
Background
Exciter saturation coefficient C in parameter calculation of excitation system of large-scale single-machine excitation generator set 1 And C 2 Is a key parameter, which is defined as follows:
in which I Ef After the no-load characteristic of the exciter is tested, the per unit value, U, of the exciting current of the exciter is measured fd The per unit value of the exciting voltage is measured after the no-load characteristic of the exciter is tested.
The parameter directly influences the dynamic characteristics of the generator excitation system and is one of the key parameters of the generator excitation system. The accurate calculation of the saturation coefficient by the exciter no-load curve is the key for accurately calculating the parameters of the generator excitation system.
The widely applied calculation method at present only uses an exciter air gap line and a no-load saturation curve, and calculates exciter exciting current values at 0.75 times and 1 time of maximum exciting voltage values by solving the following equation sets:
in which I fB 、I fJ Is the exciting current value at the position of the exciter air gap line corresponding to 1 time and 0.75 time of the maximum exciting voltage value respectively, I f0 、I fK Is the exciting current value, U, at the position corresponding to the maximum voltage value of 1 time and 0.75 time on the no-load curve of the exciter fM Is the maximum excitation voltage.
The saturation coefficient is determined by an exciter no-load characteristic test curve according to the formula, but the test result of the exciter no-load curve in the actual field test is influenced by various interference factors in the environment and often contains certain disturbance. And saturation factor pair I f0 、I fK The numerical value of (2) is high in requirement, large errors can occur when identification is carried out by means of manual experience, and the consistency of different interpretation results is difficult to guarantee.
Disclosure of Invention
The invention provides an exciter saturation coefficient calculation method based on a secant iteration method, aiming at a no-load curve obtained in a no-load characteristic test, the saturation coefficient can be calculated under the condition of containing interference data, the accurate interpretation of the test data is realized, and an accurate basis is provided for the parameter calculation of a generator excitation system.
An exciter saturation coefficient calculation method based on a secant iteration method comprises the following steps:
(1) Through the no-load characteristic test of the exciter, a group of no-load test data of the exciter are obtained:
X={I i ,U i },i=1,2,…,m;
wherein U is i For measured excitationMagnetic voltage value, I i M is a positive integer greater than or equal to 2 for the measured exciter exciting current value;
(2) The following set of iterative equations is established by the principle of least squares:
wherein C 1 And C 2 The exciter saturation coefficient can be obtained by the principle of least square method, and the saturation coefficient C satisfying the above formula 1 And C 2 The fitted curve has the minimum mean square error;
(3) From variable C 2 To a certain initial value C 2,1 Starting from the first two formulas of the step (2), solving the variable C 1 Two initial values of C 1,1 (1) And C 1,1 (2) :
(4) Two C to be solved 1 Substituting the initial values into the discriminant together:
F 1 =C 1,1 (1) -C 1,1 (2)
(5) If the result of the above-mentioned discriminant is greater than 0, C is decreased 2 Is such that it is C 2,2 (ii) a If the result is less than 0, increase C 2 Is such that it is C 2,2 :
C 2,2 =C 2,1 -0.1,if(C 1,1 (1) -C 1,1 (2) )>0,
C 2,2 =C 2,1 +0.1,if(C 1,1 (1) -C 1,1 (2) )<0,
(6) Using C obtained in step (5) 2 New value of C 2,2 Repeating the formula of steps (3) - (4) to obtain a new value F 2 :
(7) Calculating C by using a secant iteration formula 2 Correction value of (C) 2,3 :
(8) And (5) repeating the steps (6) to (7) and performing iterative solution:
(9) When the absolute value of the third equation in step (8) is less than 10 -9 Then, the saturation coefficient C satisfying the condition can be obtained 1 And C 2 In which C is 1 Take the last calculation C 1 (1) And C 1 (2) Average value of (a).
The invention has the beneficial effects that:
1. the method has more accurate and consistent calculation results than manual interpretation, and avoids the problem of different manual interpretation results;
2. fitting is carried out based on the principle of a least square method, and the mean square error of a fitting curve is minimum;
3. the used secant iteration method can obtain the calculation result in limited iteration times.
Drawings
FIG. 1 is a schematic representation of test data employed in the present invention;
FIG. 2 is a saturation coefficient C calculated using the present invention 2 Schematic diagram of the iteration.
Detailed Description
The technical solution of the present invention will be clearly and completely described below with reference to the accompanying drawings and specific examples.
The no-load characteristic curve obtained by measuring the no-load characteristic test of a certain synchronous generator exciter is shown in figure 1. Reading the test data to obtain a set of no-load characteristic data as follows:
where U is the measured excitation voltage value and I is the measured exciter excitation current value.
Given saturation coefficient C 2,1 Initial value of (C) 2,1 =1, substituting the first and second equations of the iterative equation set:
can be solved to obtain the saturation coefficient C 1 Two initial values of C 1,1 (1) And C 1,1 (2) :
And (3) calculating a discriminant:
F 1 =C 1,1 (1) -C 1,1 (2) =-1.97×10 -7
the discriminant is less than 0, and C is increased 2 The numerical value of (A):
C 2,2 =1+0.1=1.1
will C 2,2 Substituting the iteration equation set to continuously calculate C 1 :
And (3) calculating a discriminant:
F 2 =C 1,2 (1) -C 1,2 (2) =-4.51×10 -8
calculating C by secant formula 2 :
The iterative loop calculation is entered, and the results are shown in the following table:
when iterating to the fifth time, the absolute value of the discriminant is less than 10 -9 And ending iteration, wherein the obtained saturation coefficient is as follows:
saturation coefficient C 2 The iterative process of (a) is shown in figure 2.
Comparing the variance of the algorithm of the present invention with the conventional algorithm, the test data can be read:
solving the following equation set to obtain C 1 And C 2 :
Obtaining by solution:
the variance comparing the conventional algorithm to the inventive cut-line iteration is as follows:
therefore, the variance of the secant iteration method is far smaller than that of the conventional algorithm, and the fitting precision is higher.
The above description is only for the specific embodiments of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are also within the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (1)
1. An exciter saturation coefficient calculation method based on a secant iteration method is characterized by comprising the following steps of:
(1) Through the no-load characteristic test of the exciter, a group of no-load test data of the exciter are obtained:
X={I i ,U i },i=1,2,…,m;
wherein U is i For measured excitation voltage value, I i M is a positive integer greater than or equal to 2 for the measured exciter exciting current value;
(2) The following set of iterative equations is established by the principle of least squares:
wherein C 1 And C 2 For the exciter saturation systemNumber, as known from the principle of least squares, the saturation coefficient C satisfying the above formula 1 And C 2 The fitted curve has the minimum mean square error;
(3) From variable C 2 To a certain initial value C 2,1 Starting from the first two formulas of the step (2), solving the variable C 1 Two initial values of (C) 1,1 (1) And C 1,1 (2) :
(4) Two C to be solved 1 The initial values are substituted into the discriminant together:
F 1 =C 1,1 (1) -C 1,1 (2)
(5) If the result of the above-mentioned discriminant is greater than 0, C is decreased 2 Is such that it is C 2,2 (ii) a If the result is less than 0, increase C 2 Is such that it is C 2,2 :
C 2,2 =C 2,1 -0.1,if(C 1,1 (1) -C 1,1 (2) )>0,
C 2,2 =C 2,1 +0.1,if(C 1,1 (1) -C 1,1 (2) )<0,
(6) Using C obtained in step (5) 2 New value of (C) 2,2 Repeating the formula of steps (3) - (4) to obtain a new value F 2 :
(7) Calculating C by using a secant iteration formula 2 Correction value of (C) 2,3 :
(8) And (5) repeating the steps (6) to (7) and performing iterative solution:
(9) When the absolute value of the third equation in step (8) is less than 10 -9 Then, the saturation coefficient C satisfying the condition can be obtained 1 And C 2 In which C is 1 Take the last calculation C 1 (1) And C 1 (2) Average value of (a).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610208149.8A CN105808967B (en) | 2016-04-06 | 2016-04-06 | Exciter saturation coefficient computational methods based on secant iterative method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610208149.8A CN105808967B (en) | 2016-04-06 | 2016-04-06 | Exciter saturation coefficient computational methods based on secant iterative method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105808967A CN105808967A (en) | 2016-07-27 |
CN105808967B true CN105808967B (en) | 2018-01-02 |
Family
ID=56459597
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610208149.8A Active CN105808967B (en) | 2016-04-06 | 2016-04-06 | Exciter saturation coefficient computational methods based on secant iterative method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105808967B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101272116A (en) * | 2008-04-30 | 2008-09-24 | 山东电力研究院 | Modeling method of generator excitation system |
CN104135018A (en) * | 2014-08-04 | 2014-11-05 | 东南大学 | Distribution method for angular speed non-linear excitation controllers |
CN105403800A (en) * | 2015-12-30 | 2016-03-16 | 华北电力科学研究院有限责任公司 | Electric power system parameter determination method and device |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100794245B1 (en) * | 2006-08-22 | 2008-01-11 | 한국전기연구원 | An intelligent monitoring system of the reactive power limit of generator using machine model parameters and method the same |
-
2016
- 2016-04-06 CN CN201610208149.8A patent/CN105808967B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101272116A (en) * | 2008-04-30 | 2008-09-24 | 山东电力研究院 | Modeling method of generator excitation system |
CN104135018A (en) * | 2014-08-04 | 2014-11-05 | 东南大学 | Distribution method for angular speed non-linear excitation controllers |
CN105403800A (en) * | 2015-12-30 | 2016-03-16 | 华北电力科学研究院有限责任公司 | Electric power system parameter determination method and device |
Also Published As
Publication number | Publication date |
---|---|
CN105808967A (en) | 2016-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111046327B (en) | Prony analysis method suitable for low-frequency oscillation and subsynchronous oscillation identification | |
CN103809023B (en) | Synchronized harmonic phasor measuring method based on binary search | |
CN109142865B (en) | Frequency domain spectrum identification method considering polarization equivalent circuit parameters of oiled paper insulation interface | |
CN103399213B (en) | A kind of high-tension apparatus dielectric dissipation factor computing method based on fundamental frequency correction | |
CN113821938B (en) | Short-term prediction method and device for metering error state of mutual inductor | |
CN106383287A (en) | Method for improving power grid stable section measurement precision through improving data synchronism | |
CN103902451A (en) | Intelligent electricity meter software quality evaluation method | |
CN102323198B (en) | Method and system for correcting rock electricity experiment parameters m and n of core | |
CN109870635B (en) | Oiled paper insulation aging state evaluation method and system based on activation energy iterative correction | |
CN102904519A (en) | Robust excitation system parameter identification method based on exponential-type target function | |
CN105808967B (en) | Exciter saturation coefficient computational methods based on secant iterative method | |
US7577536B1 (en) | Determination of ohmic losses in electrical devices | |
CN104849530B (en) | A kind of measuring method of MOA resistive current first harmonics | |
JP6161783B2 (en) | Method for obtaining impedance of power transmission / distribution network by computer support, power generation apparatus and computer program for implementing the method | |
CN112511056A (en) | Robust generator dynamic state estimation method based on phasor measurement | |
CN105548882B (en) | Generator no-load characteristic linearity range discrimination method based on linear condensation degree | |
Serov et al. | Reseach of influency of mains frequency deviation on voltage spectrum measurement error by DFT method | |
CN116087698A (en) | Method and system for evaluating thermal ageing state of high-voltage sleeve | |
CN107358541A (en) | A kind of method for detecting abnormality of electric analog data | |
CN107917733A (en) | A kind of conductive structure thickness based on model and Eddy Conductivity detection method | |
CN113675883A (en) | Method for identifying control delay parameters of grid-connected inverter | |
CN102930077B (en) | Error-resistant excitation system parameter identification method based on improved target function | |
CN109975739B (en) | High-precision intelligent electric energy meter debugging and measuring method | |
JP2012037448A (en) | Thunder impulse withstand voltage test system, reference waveform calculation program, and thunder impulse withstand voltage test method | |
CN108038279B (en) | Method and device for evaluating electromechanical transient recovery process of fan under asymmetric power grid fault |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |