Technical Field
This invention relates generally to a control method and apparatus for antisurge control of turbocompressors having Surge Limit Lines with small slopes. More particularly, it relates to a method which determines a surge line and an operating point (both associated with a turbocompressor) by using specific combinations of invariant coordinates.
Background Art
A turbocompressor's Surge Limit Line, displayed in coordinates of volumetric flow rate (Q.sub.s) and polytropic head (H.sub.p), can be difficult to characterize if the slope of the line is small; that is, nearly horizontal. And it can be especially difficult to characterize if the surge line exhibits a local maximum or minimum, or both. This is often the case with axial compressors having adjustable inlet guide vanes, and for centrifugal compressors with both variable inlet guide vanes and diffuser vanes.
The present-day method of transforming a Surge Limit Line to common invariant spaces, such as reduced flow rate and pressure ratio (q.sub.r.sup.2, R.sub.c) or reduced flow rate and reduced head (q.sub.r.sup.2,h.sub.r), does not diminish the characterizing problem.
Surge Limit Lines exhibiting local maxima and minima are not functions of pressure ratio or of reduced head since the relationship is not one-to-one for either surge line. For this reason, it is impossible to accurately describe them (even for constant equivalent speed) using the standard approach to construct an antisurge parameter, S.sub.s =f(R.sub.c)/q.sub.r.sup.2.
Summary of the Invention
The purpose of this invention is to provide antisurge control for turbocompressors having Surge Limit Lines with small slopes (nearly horizontal). This proposed control method will easily and accurately determine a surge limit and an operating point (both associated with a turbocompressor) by using combinations of invariant coordinates that differ from those revealed in the prior art. The emphasis of this new technique is directed to axial compressors having adjustable inlet guide vanes, and to centrifugal compressors with both variable inlet guide vanes and diffuser vanes; although the method has application with many types of turbocompressors.
Turbocompressor antisurge control algorithms should compensate for variations in suction conditions; this is accomplished by calculating the operating point and the Surge Limit Line, utilizing specific coordinates referred to as invariant coordinates which are derived by using the notations of similitude or dimensional analysis. The result is that the surge limit is invariant (stationary) to suction conditions. Subsequently, the key to this invention is that any combination (linear or nonlinear) of invariant coordinates is also invariant; and the invariant coordinates of interest are: ##EQU1## where: p.sub.d =absolute pressure at discharge
p.sub.s =absolute pressure in suction ##EQU2## .DELTA.p.sub.o =differential pressure flow measurement (in suction or discharge)
p=pressure
N=rotational speed
Z=compressibility
R=gas constant
T=temperature
The derivation of the above invariant coordinates is introduced on pages 3-8 of ASME technical publication (96-GT-240) by Batson entitled, "Invariant Coordinate Systems for Compressor Control," which is incorporated herein by reference. Furthermore, information pertaining to invariant coordinates, applicable to this invention, is described in U.S. Pat. No. 5,508,943 by Batson and Narayanan entitled, "Method and Apparatus for Measuring the Distance of a Turbocompressor's Operating Point to the Surge Limit Interface," most particularly relative to the specification therein, which patent is incorporated herein by reference.
Brief Description of the Drawings
FIG. 1 shows a compressor map in the invariant space (q.sub.r.sup.2,R.sub.c) with Surge Limit Lines exhibiting a local maximum and minimum.
FIG. 2 shows a compressor map in the transformed space (q.sub.r.sup.2,R.sub.c q.sub.r.sup.n).
FIG. 3 shows a compressor map in the transformed space (q.sub.r.sup.2,R.sub.c q.sub.r.sup.n,N.sub.e).
FIG. 4A shows a schematic diagram of a turbocompressor and its control scheme for calculating a quantity (invariant coordinate), R.sub.c.sup.m q.sub.r.sup.n, using p.sub.s for the calculation of q.sub.r.sup.2.
FIG. 4B shows a schematic diagram of a turbocompressor and its control scheme for calculating a quantity (invariant coordinate), R.sub.c.sup.m q.sub.r.sup.n, using p.sub.d for the calculation of q.sub.r.sup.2.
FIG. 4C shows a schematic diagram of a turbocompressor and its control scheme for calculating a quantity (invariant coordinate), R.sub.c.sup.m q.sub.r.sup.n, using T.sub.d for the calculation of N.sub.e.
FIG. 4D shows a schematic diagram of a turbocompressor and its control scheme for calculating a quantity (invariant coordinate), R.sub.c.sup.m q.sub.r.sup.n, using T.sub.d for the calculation of N.sub.e, and p.sub.d for the calculation of q.sub.r.sup.2.
FIG. 5A shows a schematic diagram of a turbocompressor and its control scheme for calculating a quantity (invariant coordinate), h.sub.r.sup..alpha. q.sub.r.sup..beta., using p.sub.s for the calculation of q.sub.r.sup.2.
FIG. 5B shows a schematic diagram of a turbocompressor and its control scheme for calculating a quantity (invariant coordinate), h.sub.r.sup..alpha. q.sub.r.sup..beta., using p.sub.d for the calculation of q.sub.r.sup.2.
FIG. 5C shows a schematic diagram of a turbocompressor and its control scheme for calculating a quantity (invariant coordinate), h.sub.r.sup..alpha. q.sub.r.sup..beta., using T.sub.d for the calculation of N.sub.e.
FIG. 5D shows a schematic diagram of a turbocompressor and its control scheme for calculating a quantity (invariant coordinate), h.sub.r.sup..alpha. q.sub.r.sup..beta., using T.sub.d for the calculation of N.sub.e, and p.sub.d for the calculation of q.sub.r.sup.2.
Best Mode for Carrying Out the Invention
The task of characterizing a Surge Limit Line, displayed in volumetric flow rate (Q.sub.s) and polytropic head (H.sub.p), can be difficult or even impossible if the slope of this line is nearly horizontal. Transforming this surge limit to the commonly used invariant spaces (q.sub.r.sup.2,R.sub.c) or (q.sub.r.sup.2,h.sub.r) does not reduce the problem. This situation is not unusual with variable geometry compressors, such as axial compressors having adjustable inlet guide vanes, and centrifugal compressors with both variable inlet guide vanes and diffuser vanes. The characterizing task is even more difficult if the Surge Limit Lines exhibit a local maximum or minimum, or both.
In the coordinates (q.sub.r.sup.2,R.sub.c) or (q.sub.r.sup.2,h.sub.r), the surge limit for variable geometry compressors is fixed for a constant equivalent speed (N.sub.e). In other words, the surge limit is a surface that is intersected with a plane of constant N.sub.e to reduce it to a single curve--valid only for that value of N.sub.e.
FIG. 1 is an example of a compressor map with Surge Limit Lines exhibiting a local maximum and minimum in the invariant space (q.sub.r.sup.2,R.sub.c). In the construction of this figure, speed and molecular weight were constant, only the temperature caused the equivalent speed to vary. These surge lines are not functions of R.sub.c since the relationship is not one-to-one. For this reason, it is impossible to accurately describe the lines (even for constant N.sub.e) by using the standard approach to construct an antisurge parameter: ##EQU3##
As mentioned in Disclosure of the Invention, any combination (linear or nonlinear) of invariant coordinates is also invariant. Therefore, the characterizing problem stated above is easily remedied by utilizing the combinations R.sub.c q.sub.r.sup.n and h.sub.r q.sub.r.sup.n where n is real (n.di-elect cons.R). Results of the transformation of the Surge Limit Lines of FIG. 1 are now depicted in FIG. 2 which shows the lines at different temperatures. Over most of their range, the slope of these new curves is significantly greater than those in (q.sub.r.sup.2,R.sub.c) space. Moreover, the slope is, everywhere, nonnegative. FIG. 2 also shows that the surge limit interface is a surface in the coordinates (q.sub.r.sup.2,q.sub.r.sup.n R.sub.c.sup.m,N.sub.e) where, in this figure, n=m=1. Therefore, a parameter indicating proximity to surge could be constructed as ##EQU4## where the function f.sub.1 (.multidot.) returns the value of q.sub.r.sup.2 at surge.
FIG. 3 shows how the function in the numerator of Eq. (1) can be separated into two, such as ##EQU5## where the product of the functions f(.multidot.) and g(.multidot.) returns the value of q.sub.r.sup.2 at surge.
The form of the proximity to the surge variable defined by Eq. (2) is easier to commission than that of Eq. (1). Each of the functions f(q.sub.r.sup.n R.sub.c.sup.m) and g(N.sub.e) can be determined separately. In both FIG. 2 and FIG. 3, similar results would be realized if reduced head (h.sub.r) were used rather than pressure ratio (R.sub.c).
More generally, functions defined as R.sub.c.sup.m q.sub.r.sup.n and h.sub.r.sup..alpha. q.sub.r.sup..beta. (where m, n, .alpha., and .beta. are real-valued exponents) are useful for these cases. Usually, m and n, or .alpha. and .beta. will not be unique for a given Surge Limit Line, but are chosen to (1) eliminate regions of negative slope, and (2) provide for simple and accurate characterization as a function of N.sub.e.
FIG. 4A shows a schematic diagram of a turbocompressor installation and its control scheme for determining a surge line and an operating point as functions of a quantity (invariant coordinate), R.sub.c.sup.m q.sub.r.sup.n. The installation incorporates transmitters for detecting and generating a rotational speed signal (N) 1 and a suction temperature signal (T.sub.s) 2. Also included are devices for detecting and generating process input signals, such as differential pressure (.DELTA.p.sub.o) 3 across a differential flow measurement device, suction pressure (p.sub.s) 4, and discharge pressure (p.sub.d) 5.
In the block-diagram portion of FIG. 4A, speed and suction temperature transmitter data (N and T.sub.s) 1, 2 are acted on by algebraic operations to produce an equivalent speed value (N.sub.e) 6 which is then characterized as a function, g(N.sub.e) 7. Concurrently, a module calculates a reduced flow rate (q.sub.r.sup.2) 8 by dividing the differential pressure signal by the suction pressure signal. Following that, the quotient (q.sub.r.sup.2) is taken to the n/2 power as a reduced flow parameter (q.sub.r.sup.n) 9. Another module calculates a pressure ratio (R.sub.c) 10 which is taken to the m power as a pressure ratio parameter (R.sub.c.sup.m) 1, then multiplied by q.sub.r.sup.n to yield the product R.sub.c.sup.m q.sub.r.sup.n 12, and characterized as a function, f(R.sub.c.sup.m q.sub.r.sup.n) 13. At this time, a divider calculates a ratio 14 ##EQU6## that is multiplied by g(N.sub.e), resulting in a modified antisurge parameter 15 ##EQU7## Note that this modified antisurge parameter equation differs from the standard approach version by the inclusion of the modification factor g(N.sub.e), and of the reduced flow parameter (q.sub.r.sup.n) in the dividend: f(R.sub.c.sup.m q.sub.r.sup.n) versus f(R.sub.c).
The modified antisurge parameter (S.sub.s) is summed with a safety margin (b) to construct a surge control line (S=S.sub.s +b) 16. The value S is then processed by a subtracter to define the distance between an operating point and a surge control line (.delta.=1-S) 17. Finally, the distance (.delta.) and a predetermined set point (SP) are both transmitted to a Proportional-Integral-Differential (PID) controller 18 and processed into a signal for modulating an end control element 19.
FIG. 4B shows an identical diagram layout to that of FIG. 4A, but with a discharge pressure signal (p.sub.d) as the divisor for calculating a reduced flow rate (q.sub.r.sup.2) 8.
FIG. 4C shows an identical diagram layout to that of FIG. 4A, but with discharge temperature data (T.sub.d) 2 used to produce an equivalent speed value (N.sub.e) 6.
FIG. 4D shows an identical diagram layout to that of FIG. 4C, but with a discharge pressure signal (p.sub.d) as the divisor for calculating a reduced flow rate (q.sub.r.sup.2) 8.
FIG. 5A shows a schematic diagram of a turbocompressor installation and its control scheme for determining a surge line and an operating point as functions of a quantity (invariant coordinate), h.sub.r.sup..alpha. q.sub.r.sup..beta.. The installation incorporates transmitters for detecting and generating a rotational speed signal (N) 20, a suction temperature signal (T.sub.s) 21, and a discharge temperature signal (T.sub.d) 22. Also included are devices for detecting and generating process input signals, such as differential pressure (.DELTA.p.sub.o) 23 across a differential flow measurement device, suction pressure (p.sub.s) 24, and discharge pressure (p.sub.d) 25.
In the block-diagram portion of FIG. 5A, speed and temperature transmitter data (N and T.sub.s) 20, 21 are acted on by algebraic operations to produce an equivalent speed value (N.sub.e) 26 which is then characterized as a function, g(N.sub.e) 27. Concurrently, a module calculates a reduced flow rate (q.sub.r.sup.2) 28 by dividing the differential pressure signal by the suction pressure signal. Following that, the quotient (q.sub.r.sup.2) is taken to the .beta./2 power to calculate a reduced flow parameter (q.sub.r.sup..beta.) 29. Another module calculates a pressure ratio (R.sub.c) 30, while yet another calculates a temperature ratio (R.sub.T) 31. Next, the logarithms of both the pressure ratio log(R.sub.c)! 32 and the temperature ratio log(R.sub.T)! 33 are computed and jointly acted on by a divider to calculate an exponent (.sigma.) 34.
The pressure ratio (R.sub.c) is taken to the .sigma. power 35, reduced by unity (R.sub.c.sup..sigma. -1) 36, and divided by the exponent (.sigma.) to calculate a reduced head (h.sub.r) 37. Following that, reduced head is taken to the .alpha. power as a reduced head parameter (h.sub.r.sup..alpha.) 38, then multiplied by q.sub.r.sup..beta. to yield the product h.sub.r.sup..alpha. q.sub.r.sup..beta. 39, and characterized as a function, f(h.sub.r.sup..alpha. q.sub.r.sup..beta.) 40. At this time, a divider calculates a ratio 41 ##EQU8## that is multiplied by g(N.sub.e), resulting in a modified antisurge parameter 42 ##EQU9## Note that this modified parameter equation differs from the standard approach version by the inclusion of the modification factor g(N.sub.e), and of the invariant coordinate (h.sub.r.sup..alpha. q.sub.r.sup..beta.) in the dividend: f(h.sub.r.sup..alpha. q.sub.r.sup..beta.) versus f(R.sub.c).
The modified surge parameter (S.sub.s) is summed with a safety margin (b) to construct a surge control line (S=S.sub.s +b) 43. The value S is then processed by a subtracter to define the distance between an operating point and a surge control line (.delta.=1-S) 44. Finally, the distance (.delta.) and a predetermined set point (SP) are both transmitted to a PID controller 45 and processed into a signal for modulating an end control element 46.
FIG. 5B shows an identical diagram layout to that of FIG. 5A, but with a discharge pressure signal (p.sub.d) as the divisor for calculating a reduced flow rate (q.sub.r.sup.2) 28.
FIG. 5C shows an identical diagram layout to that of FIG. 5A, but with discharge temperature data (T.sub.d) 22 used to produce an equivalent speed value (N.sub.e) 26.
FIG. 5D shows an identical diagram layout to that of FIG. 5C, but with a discharge pressure signal (p.sub.d) as the divisor for calculating a reduced flow rate (q.sub.r.sup.2) 28.
In all turbocompressor installation schematics (FIGS. 4A-D and 5A-D), the flow measurement device, FT, can be located in either the suction or the discharge of the compressor.
This new technique has been described as being applicable to axial compressors having adjustable inlet guide vanes, and to centrifugal compressors with both variable inlet guide vanes and diffuser vanes; however the method has application with many types of turbocompressors not of the aforementioned types.
Obviously, many modifications and variations of the present inventions are possible in light of the above teachings. It is, therefore, to be understood that within the scope of the appended claims, the invention may be practiced otherwise than as specifically described.