Method and system for calculating stability margin domain of control parameters of hydraulic turbine regulating system

Abstract

The disclosure provides a method and a system for calculating a stability margin domain of the control parameters of a hydraulic turbine regulating system, which belongs to the technical field of stability margin analysis of the hydraulic turbine regulating system. The method includes: determining a high-order state space model of the hydraulic turbine regulating system; determining a dominant eigenvalue and an interference eigenvalue of a state matrix of the high-order state space model; determining an associated state variable of the interference eigenvalue; determining the interference eigenvalues under different target control parameter conditions; determining whether a target control parameter meets a stability margin requirement; and determining the stability margin domain of the control parameter. The present disclosure can accurately calculate the stability margin domain of the complex hydraulic turbine system and help operators quickly identify the control parameters adjustment range meeting adjustment and attenuation requirements of the system.

Claims

1. A method for calculating a stability margin domain of the control parameters of a hydraulic turbine regulating system, comprising following steps: step 1: constructing a high-order state space model of the hydraulic turbine regulating system; wherein, the hydraulic turbine regulating system comprises: a PID governor, a water diversion module comprising a pressure regulating well, a hydraulic turbine, and a generator module; step 2: randomly selecting a set of the control parameter within a PI stability domain, finding a conjugate pole without any zero within a preset distance and closest to an imaginary axis in a zero-pole distribution diagram after converting the high-order state space model as a dominant eigenvalue, and determining an eigenvalue closer to the imaginary axis than the dominant eigenvalue as an interference eigenvalue; step 3: calculating a participation factor of each state quantity in the hydraulic turbine regulating system to the interference eigenvalue, and determining the state variable corresponding to the participation factor exceeding a preset value as an associated state variable of the interference eigenvalue; step 4: inputting a target control parameter into the high-order state space model to calculate the eigenvalue, and judging a certain eigenvalue as the interference eigenvalue if a ratio of the participation factor of the associated state variable to the certain eigenvalue is more than 50% of the participation factor of all the state variables to the certain eigenvalue; step 5: judging whether a real part of the eigenvalue after removing the interference eigenvalue is all less than-m, if yes, the target control parameter meets a stability margin requirement; wherein m is the stability margin requirement; step 6: moving the control parameter (0,0) according to a step size, updating the target control parameter by adopting the control parameter obtained after each movement, and proceeding the step 4 until all target control parameter sets meeting the stability margin requirement are obtained to form a parameter stability margin domain; the step 4 specifically comprises following steps: step 4.1: inputting the target control parameter into the high-order state space model to solve for the eigenvalue of a state matrix; step 4.2: calculating the participation factor of all the state variable to all the eigenvalue using the participation factor; step 4.3: calculating a ratio ${}_{g,k}/\underset{i=1}{\overset{n}{.Math.}}{}_{i,k}$ of the participation factor of a current state variable to a current eigenvalue in the participation factor of all the state variable to the current eigenvalue, and considering the current eigenvalue to be the interference eigenvalue if ${}_{g,k}/\underset{i=1}{\overset{n}{.Math.}}{}_{i,k}>0.5;$ wherein, .sub.g,k is the participation factor of a current g-th state variable to a current k-th eigenvalue; n is a total number of the eigenvalue; step 4.4: updating the current eigenvalue to a next eigenvalue, and executing step 4.3 until all the interference eigenvalue is found.

2. The method for calculating the stability margin domain of the control parameters of the hydraulic turbine regulating system according to claim 1, wherein the step 2 specifically comprises following steps: step 2.1: randomly selecting the set of the control parameter within the PI stability domain, and substituting the control parameter and a system parameter into the high-order state space model to solve for the eigenvalue; and then converting the high-order state space model into a zero-pole model based on the solved eigenvalue to construct the zero-pole distribution diagram; step 2.2: finding the conjugate pole closest to the imaginary axis and without any zero nearby as the dominant eigenvalue according to a pole and a zero in the zero-pole distribution diagram, and taking the pole or the zero with the real part being greater than the real part of the dominant eigenvalue as the interference eigenvalue.

3. The method for calculating the stability margin domain of the control parameters of the hydraulic turbine regulating system according to claim 1, wherein a calculation method of the participation factor is: $\begin{array}{cc}{}_{i,k}=.Math.\frac{v_{i,k}{u}_{i,k}}{V_k^T{U}_k}.Math.& V_k^{T}A={}_k{V}_k,A{U}_k={}_k{U}_k\end{array}$ wherein .sub.i,k is the participation factor of an i-th state variable of the system to a k-th eigenvalue, V.sub.k and U.sub.k are respectively a left eigenvector and a right eigenvector corresponding to an eigenvalue .sub.k; v.sub.i,k and u.sub.i,k are respectively i-th row elements of V.sub.k and U.sub.k; and A is a system matrix.

4. The method for calculating the stability margin domain of the control parameters of the hydraulic turbine regulating system according to claim 1, wherein the step 6 specifically comprises following steps: step 6.1: taking a PI stability domain plane coordinate origin (0,0) as a reference point, keeping a value K.sub.p constant, giving a step size dk.sub.i in a direction K.sub.i, finding a next target control parameter (0, K.sub.i=K.sub.i+dk.sub.i) in the direction K.sub.i to be substituted into the high-order state space model, judging whether the target control parameter meets the stability margin requirement according to the step 4 and the step 5, and continuously updating K.sub.i according to dk.sub.i until a critical stability value of the maximum K.sub.i under the current K.sub.p is updated; step 6.2: giving the step size dk.sub.p in the direction K.sub.p, updating K.sub.p=K.sub.p+dk.sub.p, returning to the step 6.1, calculating all stability value K.sub.i under the current K.sub.p until K.sub.p reaches the maximum stability value K.sub.p; and constructing the set of all the control parameter meeting the stability margin requirement as the parameter stability margin domain, wherein K.sub.p is a proportional coefficient; K.sub.i is an integral coefficient.

5. A system for calculating a stability domain of control parameters of a hydraulic turbine regulating system, comprising: a construction module of a high-order state space model, used to construct the high-order state space model of the hydraulic turbine regulating system, wherein the hydraulic turbine regulating system comprises a PID governor, a water diversion module comprising a pressure regulating well, a hydraulic turbine, and a generator module; an eigenvalue screening module, used to randomly select a set of the control parameter within a PI stability domain, find a conjugate pole without any zero within a preset distance and closest to an imaginary axis in a zero-pole distribution diagram after converting the high-order state space model as a dominant eigenvalue, and determine an eigenvalue closer to the imaginary axis than the dominant eigenvalue as an interference eigenvalue; an associated state variable obtaining module, used to calculate a participation factor of each state quantity in the hydraulic turbine regulating system to the interference eigenvalue, and determine the state variable corresponding to the participation factor exceeding a preset value as an associated state variable of the interference eigenvalue; an interference eigenvalue inverse solution module, used to input a target control parameter into the high-order state space model to calculate the eigenvalue, and judge a certain eigenvalue as the interference eigenvalue if a ratio of the participation factor of the associated state variable to the certain eigenvalue is more than 50% of the participation factor of all the state variables to the certain eigenvalue; a stability margin judgement module, used to judge whether a real part of the eigenvalue after removing the interference eigenvalue is all less than m, if yes, the target control parameter meets a stability margin requirement; wherein m is the stability margin requirement; a target control parameter updating module, used to move the control parameter (0,0) according to a step size, and update the target control parameter by adopting the control parameter obtained after each movement, wherein the set of all the control parameters meeting the stability margin requirement forms the parameter stability margin domain; the interference eigenvalue inverse solution module comprises: an eigenvalue solution module, used to input the target control parameter into the high-order state space model to solve for the eigenvalue of a state matrix; a participation factor calculation unit, used to calculate the participation factors of all the state variables to all the eigenvalue using the participation factors; an interference feature confirmation unit, used to calculate a ratio ${}_{g,k}/\underset{i=1}{\overset{n}{.Math.}}{}_{i,k}$ of the participation factor of a current state variable to a current eigenvalue in the participation factor of all the state variable to the current eigenvalue, and consider the current eigenvalue to be the interference eigenvalue if ${}_{g,k}/\underset{i=1}{\overset{n}{.Math.}}{}_{i,k}>0.5;\mathrm{wherein}{}_{g,k}$ is a participation factor of a current g-th state variable to a current k-th eigenvalue; n is a total number of the eigenvalue; an iterative driving unit, used to update the current eigenvalue to a next eigenvalue, and drive the interference feature confirmation unit to find all the interference eigenvalue.

6. The system for calculating the stability domain of the control parameters of the hydraulic turbine regulating system according to claim 5, wherein the eigenvalue screening module comprises: a zero-pole distribution diagram construction unit, used to randomly select the set of the control parameter within the PI stability domain, and substitute the control parameters and a system parameter into the high-order state space model to solve for the eigenvalue; and then convert the high-order state space model into a zero-pole model based on the solved eigenvalue to construct the zero-pole distribution diagram; an eigenvalue confirmation unit, used to find the conjugate pole closest to the imaginary axis and without any zero nearby as the dominant eigenvalue according to a pole and a zero in the zero-pole distribution diagram, and take the pole or the zero with the real part being greater than the real part of the dominant eigenvalue as the interference eigenvalue.

7. The system for calculating the stability domain of the control parameters of the hydraulic turbine regulating system according to claim 5, wherein a calculation method of the participation factor in the associated state variable obtaining module is: $\begin{array}{cc}{}_{i,k}=.Math.\frac{v_{i,k}{u}_{i,k}}{V_k^T{U}_k}.Math.& V_k^{T}A={}_k{V}_k,A{U}_k={}_k{U}_k\end{array}$ wherein .sub.i,k is the participation factor of an i-th state variable of the system to a k-th eigenvalue, V.sub.k and U.sub.k are respectively a left eigenvector and a right eigenvector corresponding to an eigenvalue .sub.k; v.sub.i,k and u.sub.i,k are respectively i-th row elements of V.sub.k and U.sub.k; and A is a system matrix.

8. The system for calculating the stability domain of the control parameters of the hydraulic turbine regulating system according to claim 5, wherein the target control parameter updating module comprises: a K.sub.i updating unit, used to take a PI stability domain plane coordinate origin (0,0) as a reference point, keep a value K.sub.p constant, give a step size dk.sub.i in a direction K.sub.i, and find a next target control parameter (0, K.sub.i=K.sub.i+dk.sub.i) in the direction K.sub.i; a K.sub.p updating unit, used to give the step size dk.sub.p in the direction K.sub.p, and update K.sub.p=K.sub.p+dk.sub.p, wherein K.sub.p is a proportional coefficient; K.sub.i is an integral coefficient.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a flowchart of a method for calculating a stability margin domain of the control parameters of a hydraulic turbine regulating system provided by an embodiment of the disclosure.

(2) FIG. 2 is a model block diagram of a hydraulic turbine regulating system provided by an embodiment of the disclosure.

(3) FIG. 3 is a zero-pole distribution diagram of a system provided by an embodiment of the disclosure.

(4) FIG. 4 is a preliminary determination diagram of a system dominant eigenvalue-interference eigenvalue diagram provided by an embodiment of the disclosure.

(5) FIG. 5 is a flowchart of judging whether the control parameters meets a stability margin requirement provided by an embodiment of the disclosure.

(6) FIG. 6 is a flowchart of solving a stability margin domain of the control parameters provided by an embodiment of the disclosure.

(7) FIG. 7 is a result diagram of solving a stability margin domain of the control parameters provided by an embodiment of the disclosure.

(8) FIG. 8 is a schematic diagram of a system for calculating a stability margin domain of the control parameters of a hydraulic turbine regulating system provided by an embodiment of the disclosure.

DESCRIPTION OF THE EMBODIMENTS

(9) In order for the objectives, technical solutions, and advantages of the disclosure to be clearer, the disclosure will be further described in detail below with reference to the drawings and the embodiments. It should be understood that the specific embodiments described here are only used to explain the disclosure and are not intended to limit the disclosure.

(10) The term and/or in the disclosure describes an association relationship between associated objects, indicating that there may be three relationships. For example, A and/or B may indicate the three cases where A exists alone, A and B exist simultaneously, and B exists alone. The symbol / in the disclosure indicates that the associated objects are in an or relationship. For example, A/B indicates A or B.

(11) In the embodiments of the disclosure, words such as exemplary or for example are used to represent examples, illustrations, or explanations. Any embodiment or design described as exemplary or for example in the embodiments of the disclosure should not be construed as preferred or advantageous over other embodiments or designs. Rather, the use of the words such as exemplary or for example is intended to present concepts in a specific manner.

(12) In the description of the embodiments of the disclosure, unless otherwise specified, the meaning of multiple refers to two or more than two.

Embodiment 1

(13) The embodiment takes a dual-machine tailwater hydropower system of a hydropower station as a research object. Parameters of the dual-machine tailwater hydropower system of the hydropower station are specifically shown in Table 1.

(14) TABLE-US-00001 TABLE 1 Unit number Component Parameter 1# Pressure pipe T.sub.w1 = 1.9559 s, r.sub.1 = 0.0296, L.sub.1 = 559.6 m, A.sub.1 = 83.4 m.sup.2 Draft tube T.sub.w2 = 0.4473 s, r.sub.2 = 0.0077, L.sub.2 = 228.5 m, A.sub.2 = 175.4 m.sup.2 T.sub.w3 = 2.0105 s, r.sub.3 = 0.0029, L.sub.3 = 618.7 m, A.sub.3 = 83.6 m.sup.2 2# Pressure pipe T.sub.w4 = 0.4013 s, r.sub.4 = 0.0007, L.sub.4 = 252.6 m, A.sub.4 = 171.2 m.sup.2 Draft tube T.sub.c = 272.2 s, T.sub.q = 21.5 s, F = 1452.2 m.sup.2, .sub.c = 45.7 m.sup.2 Public Pressure regulating well T.sub.wS = 1.9395 s, r.sub.5 = 0.017, pipe Tailwater tunnel L.sub.5 = 1031.8 m, A.sub.5 = 289.3 m.sup.2 Hydraulic Rated working H.sub.r = 202 m, Q.sub.r = 538.8 m.sup.3/s, turbine condition P.sub.r = 1035MW Governor / Six parameters of hydraulic turbine Dynamo [e.sub.y, e.sub.h, e.sub.x, e.sub.qy, e.sub.qh, e.sub.qx] = [1, 1.5, 1, 1, 0.5, 0] K.sub.D = 0.1, T.sub.d = 0.01, b.sub.p = 0.02 T.sub.a = 12.07, e.sub.g = 0.176, T.sub.ae = 12.07, e.sub.ge = 0.176, k.sub.1 = k.sub.2 = 0.5

(15) As shown in FIG. 1, the disclosure provides a method for calculating a stability margin domain of the control parameters of a hydraulic turbine regulating system, which includes the following steps.

(16) In step 1, as shown in FIG. 2, a state matrix of a high-order state space model of the hydraulic turbine regulating system is determined. The hydraulic turbine regulating system includes four modules, a PID governor, a water diversion module including a pressure regulating well, a hydraulic turbine, and a generator. The high-order state space model is used to express a relationship between a state variable of a multi-machine co-pipeline hydraulic turbine regulating system and a system matrix composed of each parameter module through a state equation, and the state matrix of the high-order state space model of the hydraulic turbine regulating system is determined, which is specifically introduced below.

(17) Taking a certain dual-machine co-tailwater hydraulic turbine regulating system as an example, a governor sub-model is expressed as:

(18) $\{\begin{array}{c}\frac{d{x}_D}{dt}=-\frac{K_d}{T_D}\frac{dx}{dt}-\frac{1}{T_D}{x}_D\\ x_P=K_P\left(c_1-x-b_p{y}_1\right)\\ \frac{d{x}_I}{dt}=K_i{c}_1-K_{i}x-K_i{b}_p{y}_1\\ \frac{d{y}_1}{dt}=\frac{1}{T_{y1}}\left(x_P+x_I+x_D-y_1\right)\\ \frac{dy}{dt}=\frac{1}{T_y}\left(y_1-y\right)\end{array}$ where K.sub.p is a proportional coefficient, K.sub.i is an integral coefficient, K.sub.d is a differential coefficient, T.sub.D is a differential time constant, b.sub.p is a permanent slip coefficient, x.sub.p, x.sub.1, and x.sub.D are respectively a proportional signal, an integral signal, and a differential signal, x is a rotational speed deviation relative value of the unit, c.sub.1 is a set rotational speed relative value, y.sub.1 and y are respectively an auxiliary relay output signal and a guide vane opening deviation relative value, T.sub.y1 is an auxiliary relay reaction time constant, and T.sub.y is a main relay reaction time constant.

(19) A water diversion system sub-model including the pressure regulating well may be expressed as:

(20) $\{\begin{array}{c}\frac{d{q}_{t1}}{dt}=\frac{1}{T_{w1}+T_{w2}}\left[-h_{t1}-z_s-\left(r_1+r_2\right)q_{t1}\right]\\ \frac{d{q}_{t2}}{dt}=\frac{1}{T_{w3}+T_{w4}}\left[-h_{t2}-z_s-\left(r_3+r_4\right)q_{t2}\right]\\ \frac{d{q}_5}{dt}=\frac{1}{T_{w5}}\left(-r_5{q}_5+z_s\right)\\ \frac{d{z}_s}{dt}=\frac{1}{T_c}\left[q_{t1}+q_{t2}-q_5+\frac{T_{q}d\left(q_{t1}+q_{t2}-q_5\right)}{dt}\right]\end{array}$ where T.sub.w1 and T.sub.w2 are respectively a water flow inertia time constant of a pressure pipe and a tailwater pipe of Unit No. 1, T.sub.w3 and T.sub.w4 are respectively a water flow inertia time constant of a pressure pipe and a tailwater pipe of Unit No. 2, q.sub.i1 and q.sub.i2 are respectively a deviation relative value of flow rates of No. 1 and No. 2 hydraulic turbines, h.sub.i1 and h.sub.i2 are respectively a deviation relative value of operating heads of No. 1 and No. 2 hydraulic turbines, r.sub.1 and r.sub.2 are respectively a head loss coefficient of a pressure pipe and a tailwater pipe of Unit 1, r.sub.3 and r.sub.4 are respectively a head loss coefficient of a pressure pipe and a tailwater pipe of Unit 2, T.sub.w5, r.sub.5, and q.sub.5 are respectively a water flow inertia time constant of a tailwater tunnel, a pipe head loss coefficient, and a flow rate deviation relative value, z.sub.s is a head deviation relative value of the pressure regulating well, T.sub.c is a time constant of the pressure regulating well, and T.sub.q is a flow time constant of the pressure regulating well.

(21) A hydraulic turbine sub-model may be expressed as:

(22) $\{\begin{array}{c}{q}_t=e_{qx}x+e_{\mathrm{qy}}y+e_{qh}{h}_t\\ m_t=e_{x}x+e_{y}y+e_h{h}_t\end{array}$ where m.sub.t is a torque deviation relative value of the hydraulic turbine, q.sub.t is a flow deviation relative value of the hydraulic turbine, x is a rotational speed deviation relative value, h.sub.t is a head deviation relative value of the hydraulic turbine, y is a guide vane opening deviation relative value, e.sub.y is a transmission coefficient of hydraulic turbine torque to guide vane opening, e.sub.x is a transmission coefficient of hydraulic turbine torque to rotational speed, e.sub.h is a transmission coefficient of hydraulic turbine torque to working head, e.sub.qy is a transmission coefficient of hydraulic turbine flow to guide vane opening, e.sub.qx is a transmission coefficient of hydraulic turbine flow to rotational speed, and e.sub.qh is a transmission coefficient of hydraulic turbine flow to working head.

(23) A generator and load sub-model may be expressed as:

(24) 0$T_{ae}\frac{dx}{dt}=m_{te}-m_{g0}-e_{ge}x$$\{\begin{array}{c}{T}_{ae}=\underset{i}{\overset{n}{.Math.}}{k}_i{T}_{ai}\\ m_{\mathrm{te}}=\underset{i}{\overset{n}{.Math.}}{k}_i{m}_{\mathrm{ti}}\\ e_{ge}=\underset{i=1}{\overset{n}{.Math.}}{k}_i{e}_{gi}\end{array}$ where T.sub.ae is an equivalent unit inertia time constant, m.sub.te is a deviation relative value of an equivalent active torque of the unit, m.sub.g0 is a deviation relative value of an equivalent load torque, e.sub.ge is a self-adjustment coefficient of an equivalent generator load, k.sub.i is a rated capacity ratio of an i-th unit, T.sub.ai is an inertia time constant of the i-th unit, e.sub.gi is a self-adjustment coefficient of an i-th generator load, m.sub.ti is a deviation relative value of an active torque of the i-th unit, n represents a total of n units, and there are two units here.

(25) Combining the four sub-models above, the state space model of the dual-machine co-tailwater hydraulic turbine regulating system may be expressed as:

(26) $\{\begin{array}{c}\stackrel{}{X}=A_{13,13}{X}_{13,1}+B_{13,3}{U}_{3,1}\\ Y=C_{1,13}{X}_{13,1}+D_{1,3}{U}_{3,1}\end{array}$ where X.sub.13,1=[x,y,y,y.sub.1,y.sub.1,x.sub.1,x.sub.1,x.sub.D,x.sub.D,q.sub.t,q.sub.t,q.sub.5,z.sub.s].sup.T, U=[c,c,m.sub.g].sup.T,

(27) $A_{13,13}=\left[\begin{array}{ccc}{a}_{11}& .Math.& a_{1,13}\\ .Math.& & .Math.\\ a_{13,1}& .Math.& a_{13,13}\end{array}\right],$$B_{13,3}=\left[\begin{array}{ccc}{b}_{1,1}& b_{2l}& b_{31}\\ .Math.& .Math.& .Math.\\ b_{13l}& b_{132}& b_{13,3}\end{array}\right],$$C_{1,13}=\left[c_{1,1}.Math.c_{1,13}\right],$$D_{1,3}=0.$ where A.sub.13,13 is a system state matrix, X.sub.13,1 is a system state vector, B.sub.13,3 is a control matrix, U.sub.3,1 is a control vector, {dot over (X)} is a system state differential vector, Y is a system output, C.sub.1,13 is a state output matrix, and D.sub.1,3 is a control output matrix. where state variables of Unit No. 1 and Unit No. 2 are distinguished by superscript , and the variables with the superscript represent the variables of Unit No. 2.
a.sub.1,1=(e.sub.gek.sub.1e.sub.xk.sub.2e.sub.x+k.sub.1e.sub.he.sub.qx/e.sub.qh+k.sub.2e.sub.he.sub.qx/e.sub.qh)/T.sub.a, a.sub.1,2=k.sub.1(e.sub.ye.sub.he.sub.qy/e.sub.qh)/T.sub.ae
a.sub.1,3=k.sub.2(e.sub.ye.sub.he.sub.qy/e.sub.qh)/T.sub.ae, a.sub.1,10=k.sub.1e.sub.h/(e.sub.qhT.sub.ae), a.sub.1,11=k.sub.2e.sub.h/(e.sub.qhT.sub.ae), b.sub.1,3=1/T.sub.ae
a.sub.2,2=1/T.sub.y, a.sub.2,4=1/T.sub.y, a.sub.3,3=1/T.sub.y, a.sub.3,5=1/T.sub.y
a.sub.4,1=K.sub.P/T.sub.y1, a.sub.4,4=(1+K.sub.Pb.sub.p)/T.sub.y1, a.sub.4,6=1/T.sub.y1, a.sub.4,8=1/T.sub.y1, b.sub.4,1=K.sub.P/T.sub.y1
a.sub.5,1=K.sub.P/T.sub.y1, a.sub.5,5=(1+K.sub.Pb.sub.p)/T.sub.y1, a.sub.5,7=1/T.sub.y1, a.sub.5,9=1/T.sub.y1, b.sub.5,2=K.sub.P/T.sub.y1
a.sub.6,1=K.sub.1, a.sub.6,4=K.sub.1b.sub.p, b.sub.6,1=K.sub.1, a.sub.7,1=K.sub.1, a.sub.7,5=K.sub.1b.sub.p, b.sub.7,2=K.sub.1
a.sub.8,1=K.sub.D/T.sub.d*a.sub.1,1, a.sub.8,2=K.sub.D/T.sub.d*a.sub.1,2, a.sub.8,3=K.sub.D/T.sub.d*a.sub.1,3, a.sub.8,8=1/T.sub.d
a.sub.8,10=K.sub.D/T.sub.d*a.sub.1,10, a.sub.8,11K.sub.D/T.sub.d*a.sub.1,11, b.sub.8,3=K.sub.D/T.sub.d*b.sub.1,3
a.sub.9,1=K.sub.D/T.sub.d*a.sub.1,1, a.sub.9,2=K.sub.D/T.sub.d*a.sub.1,2, a.sub.9,3=K.sub.D/T.sub.d*a.sub.1,3, a.sub.9,9=1/T.sub.d
a.sub.9,10=K.sub.D/T.sub.d*a.sub.1,10, a.sub.9,11=K.sub.D/T.sub.d*a.sub.1,11, b.sub.9,3=K.sub.D/T.sub.d*b.sub.1,3
a.sub.10,1=e.sub.qx/e.sub.qh/(T.sub.w1+T.sub.w2), a.sub.10,2=e.sub.qy/e.sub.qh/(T.sub.w1+T.sub.w2), a.sub.10,10=(r.sub.1+r.sub.2+1/e.sub.qh)/(T.sub.w1+T.sub.w2), a.sub.10,13=1/(T.sub.w1+T.sub.w2)
a.sub.11,1=e.sub.qx/e.sub.qh/(T.sub.w3+T.sub.w4), a.sub.11,3=e.sub.qy/e.sub.qh/(T.sub.w3+T.sub.w4), a.sub.11,11=(r.sub.3+r.sub.4+1/e.sub.qh)/(T.sub.w3+T.sub.w4), a.sub.11,13=1/(T.sub.w3+T.sub.w4)
a.sub.12,12=r.sub.5/T.sub.w5, a.sub.12,13=1/T.sub.w5, a.sub.13,1=T.sub.q*(a.sub.10,1+a.sub.11,1)/T.sub.c, a.sub.13,2=T.sub.q*a.sub.10,2/T.sub.c, a.sub.13,3=T.sub.q*a.sub.11,3/T.sub.c
a.sub.13,10=(1+T.sub.q*a.sub.10,10)/T.sub.c, a.sub.13,11=(1+T.sub.q*a.sub.11,11)/T.sub.c, a.sub.13,12=(T.sub.q*a.sub.12,12+1)/T.sub.c, a.sub.13,13=T.sub.q*(a.sub.10,13+a.sub.11,13a.sub.12,13)/T.sub.c where except for the elements above, the remaining elements in A.sub.13,13 and B.sub.13,3 are all 0.

(28) In step 2, a dominant eigenvalue and an interference eigenvalue of the state matrix of the high-order state space model are preliminarily determined. All conjugate eigenvalues closest to an imaginary axis are found, a zero-pole distribution diagram of the system is analyzed, a conjugate pole closest to the imaginary axis and without any zero nearby is found as a dominant pole, and an eigenvalue closer to the imaginary axis than the dominant eigenvalue is determined as the interference eigenvalue. A more specific introduction is as follows.

(29) In step 2.1, control parameters (K.sub.p=2, K.sub.i=0.15) within a PI stability domain are randomly selected, the control parameters and a system parameter (Table 1) are substituted into the high-order state space model to solve for the eigenvalue, and the high-order state space model is then converted into a zero-pole model. The conversion relationship is:

(30) $\begin{array}{c}\mathrm{sys}\left(A,B,C,D\right)=\mathrm{sys}\left(P,Z\right)\\ \{\begin{array}{cc}{}_i={}_i+j{}_i& i=1,2,.Math.,13\\ p_i={}_i& i=1,2,.Math.,13\\ z_k={}_k+j{}_k& k=1,2,.Math.,K\left(K13\right)\end{array}\end{array}$ where A, B, C, and D are matrices in the state space model, P and Z are respectively a pole and a zero in the zero-pole model, .sub.i is an i-th eigenvalue of the system, p.sub.i is an i-th pole, Z.sub.k is a k-th zero, .sub.i is a real part of the eigenvalue .sub.i, .sub.i is an imaginary part of .sub.i, .sub.k is a real part of the pole p.sub.i, and .sub.k is an imaginary part of p.sub.i. The eigenvalues are as shown in Table 2, and a zero-pole distribution of the system is as shown in FIG. 3.

(31) In Step 2.2, a conjugate eigenvalue .sub.8,9 closest to the imaginary axis and without any zero nearby is found according to P and Z in the zero-pole model, and .sub.8,9 is used a dominant eigenvalue .sub.m=.sub.m.sub.mj. The necessary and sufficient condition for the dominant eigenvalue is the conjugate pole closest to the imaginary axis and without any zero nearby. The sufficient condition for the interference eigenvalue is .sub.i>.sub.m, and an eigenvalue with a real part being greater than the real part of the dominant eigenvalue is an interference eigenvalue .sub.g=.sub.g.sub.gj. In the embodiment, according to the sufficient condition for the interference eigenvalue further being .sub.i>.sub.8,9, eigenvalues .sub.10,11 and .sub.13 with real parts being greater than the real part of the dominant eigenvalue are the interference eigenvalues, as shown in FIG. 4.

(32) TABLE-US-00002 TABLE 2 Eigenvalue .sub.i (.sub.i = .sub.i + .sub.ij, .sub.i and .sub.i are real part and imaginary part of eigenvalue) .sub.1 99.9878 .sub.2 25.0353 .sub.3 100 .sub.4 24.9971 .sub.5 2.3938 .sub.6 0.6100 .sub.7 0.8432 .sub.8,9 0.0809 0.0791j .sub.10,11 0.0220 0.0406j .sub.12 2.1768 .sub.13 0.0029

(33) In step 3, an associated state variable x.sub.g of the interference eigenvalue is determined. The participation factors of each state variable of the system to the interference eigenvalues (.sub.10,11 and .sub.13) are calculated. The state variable with a greater participation factor is the main associated state variable of the interference eigenvalue.

(34) It should be pointed out that in step 3, the associated state variable x.sub.g of the interference eigenvalue is judged using the participation factor. The participation factors of all the state variables to the interference eigenvalue are calculated. The state variable with a greater participation factor is the associated state variable x.sub.g. A calculation method of the participation factor is:

(35) ${}_{i,k}=.Math.\frac{v_{i,k}{u}_{i,k}}{V_k^T{U}_k}.Math.$$V_k^{T}A={}_k{V}_k,A{U}_k={}_k{U}_k$ where .sub.i,k is a participation factor of an i-th state variable of the system to a k-th eigenvalue, V.sub.k and U.sub.k are respectively a left eigenvector and a right eigenvector corresponding to an eigenvalue .sub.k, and v.sub.i,k and u.sub.i,k are respectively i-th row elements of V.sub.k and U.sub.k. The greater the .sub.i,k, the stronger the correlation between the i-th state variable and the k-th eigenvalue. Calculation results of the participation factors are shown in Table 3. According to the size of the participation factor, strongly correlated state variables q.sub.5 and z.sub.s of the interference eigenvalue .sub.10,11 may be determined, and strongly correlated state variables of the interference eigenvalue .sub.13 are x.sub.l and x.sub.l.

(36) TABLE-US-00003 TABLE 3 1 2 3 4 5 6 7 8 9 10 11 12 13 x y y y.sub.1 y.sub.1 x.sub.1 x.sub.1 x.sub.D x.sub.D q.sub.t1 q.sub.t2 q.sub.t5 z.sub.t6 State .sub.10.11 0.044 0 0 0 0 0.046 0.046 0 0 0.078 0.078 0.540 0.644 variable Participation factor .sub.13 0 0 0 0 0 0.5 0.5 0 0 0 0 0 0

(37) Step 1 to step 3 are mainly to obtain the strongly correlated state variables of the interference eigenvalues. Under different target control parameter conditions, the number of dominant eigenvalues and interference eigenvalues and the corresponding interference eigenvalue positions may change, but the interference eigenvalues are always strongly correlated with the fixed state variables of the system. Therefore, the previously obtained interference eigenvalues may be used as the strongly correlated state variables to calculate the participation factors of the strongly correlated variables in each eigenvalue of the system, and inversely solve for the interference eigenvalues of the system under different control parameters according to the size of the participation factor, thereby preventing the interference eigenvalue positions from changing as the control parameters change.

(38) In step 4, as shown in FIG. 5, the interference eigenvalues under different target control parameter conditions are determined. An eigenvalue of a state matrix under the current PI parameter condition is calculated, wherein if the ratio of the participation factor of the state variable x.sub.g to the eigenvalue is more than 50% of the participation factors of all the state quantities to the eigenvalue, the eigenvalue is the interference eigenvalue.

(39) A more detailed explanation is as follows.

(40) In step 4, the interference eigenvalue of the system under the current PI control parameter condition is determined using the ratio of the main associated state variable of the interference eigenvalue to the participation factor of the eigenvalue. The sufficient condition for the eigenvalue of the system is

(41) ${}_{g,k}/\underset{i=1}{\overset{n}{.Math.}}{}_{i,k}>0.5,$
that is, if a ratio of a participation factor of a strongly correlated state quantity x.sub.g of the interference eigenvalue to a certain eigenvalue is more than 50% of participation factors of all state quantities to the eigenvalue, the eigenvalue is considered to be the interference eigenvalue of the system under the current control parameters.

(42) The specific operations are as follows.

(43) In step 4.1, a current PI parameter is substituted into the state space model, and a system matrix A is updated to solve for an eigenvalue .sub.i of the matrix A.

(44) In step 4.2, participation factors .sub.i,k of all the state variables to the eigenvalue is calculated using the calculation method of the participation factor in step 3, where i=1 . . . 13, k=1 . . . 13.

(45) In step 4.3, a ratio

(46) ${}_{g,k}/\underset{i=1}{\overset{n}{.Math.}}{}_{i,k}$
of a participation factor of x.sub.s to a certain eigenvalue .sub.k in the participation factors of all the state variables to the eigenvalue is calculated, and .sub.k is considered to be the interference eigenvalue if

(47) ${}_{g,k}/\underset{i=1}{\overset{n}{.Math.}}{}_{i,k}>0.5.$

(48) In the embodiment, the interference eigenvalue among the eigenvalues is excluded. From a 1-st eigenvalue to a 13-th eigenvalue, if a ratio of participation factors of q.sub.5 and z.sub.s to a certain eigenvalue is more than 50% of all the state quantities to the eigenvalue, that is,

(49) $\underset{g=11}{\overset{12}{.Math.}}{}_{g,k}/\underset{i=1}{\overset{13}{.Math.}}{}_{i,k}>0.5,$
or if a ratio of participation factors of x.sub.l and x.sub.l to a certain eigenvalue is more than 50% of all the state quantities to the eigenvalue, the eigenvalue is excluded.

(50) In step 4.4, the judgment of step 4.3 is performed on each eigenvalue .sub.i to find all interference eigenvalues .sub.g.

(51) In practical applications, the difference in m represents the speed of the oscillation attenuation rate of the hydraulic turbine. In the prior art, it is impossible to effectively distinguish the interference eigenvalues, and it is difficult to accurately calculate the control parameter range meeting the attenuation speed requirement of the hydraulic turbine regulating system. For some complex systems, such as a multi-frequency oscillation coupled hydraulic turbine regulating system, the prior art cannot solve directly. Therefore, in step 5 and step 6, the positioning method of the interference eigenvalue in step 4 is introduced and the interference eigenvalue is eliminated to overcome the shortcomings of low calculation accuracy and inability to adapt to complex systems in the prior art. A control stability margin domain corresponding to different attenuation requirements m is established. The operators can quickly adjust the control parameters according to the stability margin domain, so that the system meets the attenuation requirements.

(52) In step 5, whether the target control parameter meets the stability margin requirement is determined. Under a target PI parameter condition, whether real parts of all eigenvalues after the interference eigenvalue is removed are less than-m (m represents the stability margin requirement) is determined. If the real part is less than-m, the PI parameter meets the stability margin requirement, and vice versa.

(53) A more specific introduction is as follows.

(54) In step 5, whether the system meets the stability margin requirement is determined using the real parts of the eigenvalues other than the interference eigenvalue. The sufficient condition for the system to meet the stability margin requirement is .sub.i.Math..sub.g,.sub.i<m (i=1 . . . n). If .sub.i.Math..sub.g,.sub.i<m (i=1 . . . n), the system does not meet the stability margin requirement, and m is the specific stability margin requirement.

(55) In step 6, as shown in FIG. 6, the stability margin domain of the control parameter is determined. A PI parameter is enumerated within a PI parameter stability domain, whether the PI parameter meets the stability margin requirement is judged according to the above judgment criteria. All PI parameters meeting the stability margin requirement are retained. A set of the parameters is the stability margin domain of the PI parameters.

(56) A more specific introduction is as follows.

(57) In step 6.1, a PI stability domain plane coordinate origin (0,0) is taken as a reference point, a step size dk.sub.i in a direction K.sub.i is given, a next target control parameter (0, K.sub.i=K.sub.i+dk.sub.i) in the direction K.sub.i is found and substituted into the high-order state space model, whether the target control parameter meets the stability margin requirement is judged according to step 4 (shown in FIG. 5), and K.sub.i is continuously updated according to dk.sub.i until a critical stability value of the maximum K.sub.i under the current K.sub.p is updated.

(58) In step 6.2, a value K.sub.p is updated, a step size dk.sub.p=0.1 in the direction K.sub.p is given, K.sub.p=K.sub.p+dk.sub.p is updated, step 6.1 is repeated, and all the stability values K.sub.i under the current K.sub.p are calculated until K.sub.p reaches the maximum stability value K.sub.p. The set of all the control parameters meeting the stability margin requirement is the parameter stability margin domain. Calculation results of the stability margin domain are as shown in FIG. 7.

Embodiment 2

(59) An embodiment of the disclosure provides a system for calculating a stability domain of the control parameters of a hydraulic turbine regulating system, which includes the following.

(60) A construction module 101 of a high-order state space model is used to construct a high-order state space model of the hydraulic turbine regulating system. The hydraulic turbine regulating system includes a PID governor, a water diversion module including a pressure regulating well, a hydraulic turbine, and a generator module.

(61) An eigenvalue screening module 201 is used to randomly select a set of control parameter within a PI stability domain, find a conjugate pole without any zero within a preset distance and closest to an imaginary axis in a zero-pole distribution diagram after converting the high-order state space model as a dominant eigenvalue, and determine an eigenvalue closer to the imaginary axis than the dominant eigenvalue as an interference eigenvalue.

(62) An associated state variable obtaining module 301 is used to calculate a participation factor of each state quantity in the hydraulic turbine regulating system to the interference eigenvalue, and determine the state variable corresponding to the participation factor exceeding a preset value as an associated state variable of the interference eigenvalue.

(63) An interference eigenvalue inverse solution module 401 is used to input a target control parameter into the high-order state space model to calculate an eigenvalue of a state matrix, and judge a certain eigenvalue of the state matrix as the interference eigenvalue if a ratio of a participation factor of the associated state variable to the certain eigenvalue of the state matrix is more than 50% of participation factors of all state variables to the certain eigenvalue of the state matrix.

(64) A stability margin judgement module 501 is used to judge whether a real part of the eigenvalue of the state matrix after removing the interference eigenvalue is all less than-m. If yes, the target control parameter meets a stability margin requirement, where m is the stability margin requirement.

(65) A target control parameter updating module 601 is used to move the control parameters (0,0) according to a step size, and update the target control parameter by adopting the control parameters obtained after each movement.

(66) The set of all the control parameters meeting the stability margin requirement forms the parameter stability margin domain.

(67) Further preferably, the eigenvalue screening module 201 includes the following.

(68) A zero-pole distribution diagram construction unit is used to randomly select the set of the control parameter within the PI stability domain, and substitute the control parameters and a system parameter into the high-order state space model to solve for the eigenvalue; and then convert the high-order state space model into a zero-pole model based on the solved eigenvalue to construct the zero-pole distribution diagram.

(69) An eigenvalue confirmation unit is used to find the conjugate pole closest to the imaginary axis and without any zero nearby as the dominant eigenvalue according to a pole and a zero in the zero-pole distribution diagram, and take the pole or the zero with a real part being greater than a real part of the dominant eigenvalue as the interference eigenvalue.

(70) Further preferably, a calculation method of the participation factor in the associated state variable obtaining module 301 is:

(71) $\begin{array}{cc}{}_{i,k}=.Math.\frac{v_{i,k}{u}_{i,k}}{V_k^T{U}_k}.Math.& V_k^{T}A={}_k{V}_k,A{U}_k={}_k{U}_k\end{array}$

(72) where .sub.i,k is a participation factor of an i-th state variable of the system to a k-th eigenvalue, V.sub.k and U.sub.k are respectively a left eigenvector and a right eigenvector corresponding to an eigenvalue .sub.k, and v.sub.i,k and u.sub.i,k are respectively i-th row elements of V.sub.k and U.sub.k.

(73) Further preferably, the interference eigenvalue inverse solution module 401 includes the following.

(74) An eigenvalue solution module is used to input the target control parameter into the high-order state space model to solve for the eigenvalue of the state matrix.

(75) A participation factor calculation unit is used to calculate the participation factors of all the state variables to all the eigenvalues using the participation factors.

(76) An interference feature confirmation unit is used to calculate a ratio

(77) 0${}_{g,k}/\underset{i=1}{\overset{n}{.Math.}}{}_{i,k}$
of a participation factor of a current state variable to a current eigenvalue in the participation factors of all the state variables to the current eigenvalue, and consider the current eigenvalue to be the interference eigenvalue if

(78) ${}_{g,k}/\underset{i=1}{\overset{n}{.Math.}}{}_{i,k}>0.5,$
where .sub.g,k is a participation factor of a current g-th state variable to a current k-th eigenvalue and n is a total number of eigenvalues.

(79) An iterative driving unit is used to update the current eigenvalue to a next eigenvalue, and drive the interference feature confirmation unit to find all the interference eigenvalues.

(80) Further preferably, the target control parameter updating module 601 includes the following.

(81) A K.sub.i updating unit is used to take a PI stability domain plane coordinate origin (0,0) as a reference point, keep a value K.sub.p constant, give a step size dk.sub.i in a direction K.sub.i, and find a next target control parameter (0, K.sub.i=K.sub.i+dk.sub.i) in the direction K.sub.i.

(82) A K.sub.p updating unit is used to give a step size dk.sub.p in the direction K.sub.p, and update K.sub.p=K.sub.p+dk.sub.p.

(83) K.sub.p is a proportional coefficient and K.sub.i is an integral coefficient.

(84) In practical applications, appropriate parameters (Kp and Ki) of a controller need to be set in a governor of a turbine to meet a speed regulating performance requirement (stability margin requirement) of the turbine. According to different stability margin requirements, any corresponding set of controller parameters (Kp and Ki) in the control parameters stability margin domain obtained by the control parameters stability domain calculating system of a turbine regulating system is found as the parameters of the controller, which can ensure that the speed regulation of the PID governor meets actual requirements. More specifically, when regulating the governor, a set of control parameters may be manually selected in the control parameter stability margin domain or the control parameter stability margin domain may be read through a processor in a computer, thereby randomly selecting a set of data to control the PID governor to ensure that the turbine regulating system is within a safe adjustment range of the control parameters, so that the oscillation of the turbine regulating system can quickly attenuate under the premise of ensuring long-term safe and stable operation of the unit.

(85) In summary, compared with the prior art, the disclosure has the following advantages.

(86) In the method and the system for calculating the stability margin domain of the control parameter of the hydraulic turbine regulating system provided by the disclosure, the high-order state space model and the zero-pole model of the hydraulic turbine regulating system are established. Through zero-pole complex plane analysis and participation factor calculation, the interference eigenvalue and the strongly correlated state variable thereof can be accurately found to facilitate locating a source of a sub-wave oscillation of the regulating system and understanding features of the sub-wave oscillation.

(87) In the method and the system for calculating the stability margin domain of the control parameter of the hydraulic turbine regulating system provided by the disclosure, the interference eigenvalue of the system is judged for the ratio of the participation factor of the system feature under different control parameter conditions being more than 50% according to the strongly correlated state variable of the interference eigenvalue. The accurate positioning of the interference eigenvalue of the system under different control parameters can be implemented. A key judgment indicator is provided for extracting the sub-wave oscillation under different parameters of the complex hydraulic turbine regulating system to reduce the difficulty for the operators to comprehensively analyze the system oscillation characteristics.

(88) In the method and the system for calculating the stability margin domain of the control parameter of the hydraulic turbine regulating system provided by the disclosure, all the control parameters are traversed according to the fixed step size within a governor PI control parameter stability domain. After eliminating the corresponding interference eigenvalue under each parameter condition, whether the real parts of the remaining eigenvalues are less than the stability margin requirement is judged, and the control parameter set meeting the stability margin requirement is determined to accurately calculate the stability margin domain of the complex hydraulic turbine system, thereby helping the operators intuitively and quickly adjust the control parameters, so that the system meets a control attenuation requirement, and the technical threshold for adjusting the control parameters of the hydraulic turbine regulating system is lowered.

(89) It should be understood that the above system is used to execute the method in the above embodiment. The corresponding program modules in the system have implementation principles and technical effects similar to those described in the above method. For the working process of the system, reference may be made to the corresponding process in the above method, which will not be described again here.

(90) Based on the method in the above embodiment, an embodiment of the disclosure provides an electronic apparatus. The apparatus may include at least one memory for storing a program and at least one processor for executing the program stored in the memory. When the program stored in the memory is executed, the processor is used to execute the method described in the above embodiment.

(91) Based on the method in the above embodiment, an embodiment of the disclosure provides a computer-readable storage medium. The computer-readable storage medium stores a computer program. When the computer program is run on a processor, the processor executes the method in the above embodiment.

(92) Based on the method in the above embodiment, an embodiment of the disclosure provides a computer program product. When the computer program product is run on a processor, the processor executes the method in the above embodiment.

(93) It can be understood that the processor in the embodiment of the disclosure may be a central processing unit (CPU), other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), field programmable gate arrays (FPGA), other programmable logic devices, transistor logic devices, hardware components, or any combination thereof. The general-purpose processor may be a microprocessor or any conventional processor.

(94) The steps of the method in the embodiment of the disclosure may be implemented through hardware or by a processor executing a software command. The software command may be composed of a corresponding software module. The software module may be stored in a random access memory (RAM), a flash memory, a read-only memory (ROM), a programmable ROM (PROM), an erasable PROM (EPROM), an electrically EPROM (EEPROM), a register, a hard disk, a mobile hard disk, a CD-ROM, or any other form of storage medium well-known in the art. An exemplary storage medium is coupled to the processor, so that the processor may read information from the storage medium and write information to the storage medium. Of course, the storage medium may also be an integral part of the processor. The processor and the storage medium may be located in the ASIC.

(95) The above embodiments may be implemented in whole or in part by software, hardware, firmware, or any combination thereof. When implemented using software, the embodiments may be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer commands. When the computer program commands are loaded and executed on a computer, processes or functions described according to the embodiments of the disclosure are generated in whole or in part. The computer may be a general-purpose computer, a specific-purpose computer, a computer network, or other programmable systems. The computer command may be stored in a computer-readable storage medium or transmitted through the computer-readable storage medium. The computer command may be transmitted from one website, computer, server, or data center to another website, computer, server, or data center through wired (for example, coaxial cable, optical fiber, and digital subscriber line (DSL)) or wireless (for example, infrared, wireless, microwave, etc.) means. The computer-readable storage medium may be any available medium that can be accessed by a computer or a data storage apparatus such as a server or a data center including one or more available media. The available media may be magnetic media (for example, floppy disk, hard disk, and magnetic tape), optical media (for example, DVD), semiconductor media (for example, solid state disk (SSD)), etc.

(96) It can be understood that various reference numerals involved in the embodiments of the disclosure are only for convenience of description and are not used to limit the scope of the embodiments of the disclosure.

(97) It is easy for persons skilled in the art to understand that the above descriptions are only preferred embodiments of the disclosure and are not intended to limit the disclosure. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the disclosure should be included in the protection scope of the disclosure.