The present disclosure relates to compressor control methods and systems. Embodiments disclosed herein specifically relate to wet compressors, in particular centrifugal wet compressors, which process gas that can contain a liquid phase, e.g. heavy hydrocarbons, water or the like.
Centrifugal compressors have been designed to process a so-called wet gas, i.e. gas that can contain a certain percentage of a liquid phase. Wet gas processing is often required in the oil and gas industry, where gas extracted from a well, such as a subsea well, can contain a liquid hydrocarbon phase, or water. For several reasons, it might be useful to know the liquid volume fraction (shortly LVF) of the gas processed by the compressor, i.e. the volume percentage of liquid in the fluid flow. Usually, the liquid volume fraction in the gas flow at the suction side of the compressor, however, is not known. Flowmeters capable of determining the liquid volume fraction are cumbersome and expensive and might not be suitable in certain applications in extreme environmental conditions.
A need therefore exists, for reliably and efficiently estimating the liquid volume fraction of a gas flowing through a compressor.
According to a first aspect, a method of determining a liquid volume fraction in a multi-phase gas processed by a compressor having a suction side and a delivery side is disclosed. The method can comprise the following steps:
a) measuring a first compressor operating parameter;
b) selecting a tentative liquid volume fraction of the gas processed by the compressor;
c) based on stored data representing a compressor operative curve for the tentative liquid volume fraction, determining an estimated value of a second compressor operating parameter as a function of the first compressor operating parameter;
d) measuring an actual value of the second compressor operating parameter;
e) comparing the actual value of the second compressor operating parameter to the estimated value of the second compressor operating parameter and determining an error therefrom;
f) based on the error, selecting a different tentative liquid volume fraction and repeating steps (c) to (e) until an error value equal to or lower than an error threshold is obtained.
The liquid volume fraction LVF contained in the gas processed by the compressor can thus be estimated without the need for direct measurement. The LVF determined by means of the above calculation can be used e.g. for adapting the anti-surge control of the compressor. An anti-surge control line can be selected based upon the liquid content in the wet gas, for optimal anti-surge operation.
The first compressor operating parameter can be the compression ratio or a parameter related thereto. In other embodiments, the first compressor operating parameter can be a parameter related to the compressor driving power, e.g. the corrected power. A definition of corrected power is given later on, reference being made to exemplary embodiments of the subject matter disclosed herein.
In some embodiments the second compressor operating parameter can be a parameter related to the compressor driving power, e.g. the corrected power. In other embodiments, the second compressor operating parameter can be the compression ratio or a parameter related thereto.
In some embodiments, the step of determining an estimated value of a second compressor operating parameter further comprises the step of:
based on stored data representing the compressor operative curve for the tentative liquid volume fraction, determining an estimated value of a third compressor operating parameter;
based on stored data representative of a further compressor operative curve for the tentative liquid volume fraction, and based on the estimated value of the third compressor operating parameter, determining the estimated value of the second compressor operating parameter.
According to a further aspect, disclosed herein is a system comprising:
a driver;
a compressor drivingly coupled to the driver and comprised of an anti-surge arrangement including an anti-surge line and an anti-surge control valve arranged there along;
a control unit functionally coupled to the anti-surge valve;
wherein the control unit is configured and controlled to perform a method as disclosed above.
According to another aspect, disclosed herein is a method of operating a wet-gas compressor, comprising the following steps:
running the compressor and processing a gas there through;
determining a liquid volume fraction of the gas at the suction side of the compressor;
selecting a surge control line as a function of the liquid volume fraction.
The method can further comprise the steps of:
providing sets of operating curves and surge control lines of the wet-gas compressor at different liquid volume fractions;
selecting the set of operating curves and respective surge control line corresponding to the determined liquid volume fraction.
According to some embodiments, the step of determining the liquid volume fraction at the suction side of the compressor can be performed repeatedly, e.g. at constant or variable time intervals, during operation of the compressor.
The step of determining the liquid volume fraction of the gas can comprise the step of detecting the amount of liquid in a multi-phase flow meter, or a step of estimating the amount of liquid, i.e. the liquid volume fraction, with an iterative method based upon operating parameters of the compressor.
According to a further embodiment, disclosed herein is a compressor system comprising:
a wet-gas compressor having a suction side and a delivery side;
an anti-surge arrangement comprising an anti-surge line fluidly coupling the delivery side and the suction side of the compressor and including an anti-surge control valve there along;
a control unit, functionally connected to the anti-surge control line, configured and arranged to: determine a liquid volume fraction of the gas at the suction side of the compressor; select a surge control line as a function of the liquid volume fraction; acting upon the anti-surge control valve to prevent the compressor from operating beyond the selected surge control line.
Features and embodiments are disclosed here below and are further set forth in the appended claims, which form an integral part of the present description. The above brief description sets forth features of the various embodiments of the present invention in order that the detailed description that follows may be better understood and in order that the present contributions to the art may be better appreciated. There are, of course, other features of embodiments of the invention that will be described hereinafter and which will be set forth in the appended claims. In this respect, before explaining several embodiments of the invention in details, it is understood that the various embodiments of the invention are not limited in their application to the details of the construction and to the arrangements of the components set forth in the following description or illustrated in the drawings. Embodiments of the invention are capable of other embodiments and of being practiced and carried out in various ways. Also, it is to be understood that the phraseology and terminology employed herein are for the purpose of description and should not be regarded as limiting.
As such, those skilled in the art will appreciate that the conception, upon which the disclosure is based, may readily be utilized as a basis for designing other structures, methods, and/or systems for carrying out the several purposes of embodiments of the present invention. It is important, therefore, that the claims be regarded as including such equivalent constructions insofar as they do not depart from the spirit and scope of embodiments of the present invention.
A more complete appreciation of the disclosed embodiments of the invention will be readily obtained as the same becomes better understood by reference to the following detailed description when considered in connection with the accompanying drawings, wherein:
The following detailed description of the exemplary embodiments refers to the accompanying drawings. The same reference numbers in different drawings identify the same or similar elements. Additionally, the drawings are not necessarily drawn to scale. Also, the following detailed description does not limit embodiments of the invention. Instead, the scope of embodiments of the invention is defined by the appended claims.
Reference throughout the specification to “one embodiment” or “an embodiment” or “some embodiments” means that the particular feature, structure or characteristic described in connection with an embodiment is included in at least one embodiment of the subject matter disclosed. Thus, the appearance of the phrase “in one embodiment” or “in an embodiment” or “in some embodiments” in various places throughout the specification is not necessarily referring to the same embodiment(s). Further, the particular features, structures or characteristics may be combined in any suitable manner in one or more embodiments.
In the following description of exemplary embodiments, methods and systems will be described wherein the liquid fraction volume (shortly LFV) is estimated and used to act upon an anti-surge control algorithm of a centrifugal compressor. More specifically, the LVF is used to optimize the surge control line used in the anti-surge algorithm. It shall however be understood that the disclosed methods and systems for LVF estimation can be used for other purposes, whenever a measure of the liquid volume fraction in a wet gas is desired or useful.
The driver 5 is drivingly coupled to the compressor 3 by means of a drive shaft 7. The compressor 3 can be a centrifugal, multi-stage compressor. The compressor 3 and the driver 5 can be integrated in a single casing, not shown, forming a motor-compressor unit.
The compressor 3 has a suction side 9 and a delivery side 11. The suction side 9 receives gas at a suction temperature Ts and at a suction pressure Ps. The pressure of the gas is boosted by the compressor 3 and gas at a delivery pressure Pd and delivery temperature Td is delivered at the compressor delivery side 11.
The compressor 3 can be provided with an anti-surge arrangement. In some embodiments the anti-surge arrangement comprises an anti-surge line with an anti-surge control valve arranged therealong, the anti-surge line fluidly connecting the delivery side 11 of the compressor 3 to the suction side 9 of the compressor 3. According to the schematic of
The anti-surge control valve 15 can be a bi-phase valve, i.e. a valve capable of handling a bi-phasic flow, containing gas and liquid.
The system 1 can be further comprised of a central control unit 17 and instrumentalities for measuring various operating parameters of the system 1. In some embodiments, a pressure transducer 21 and a temperature transducer 23 can be arranged and configured for measuring the suction pressure Ps and the suction temperature Ts. A pressure transducer 25 and a temperature transducer 27 can also be provided, to measure the delivery pressure Pd and the delivery temperature Td. In the exemplar embodiment of
The transducers 23-31 are functionally connected to the central control unit 17. This latter can be further provided with memory resources 33, wherein data representing operating curves, i.e. performance characteristics of the compressor 3 are stored. Possible operating curves useful to operate the methods of the present disclosure will be described here below. The data of the curves can be stored in the form of tables or matrices, for instance. In other embodiments, functions or algorithms can be stored to calculate the values of the operating curves.
The operating condition of the compressor 3 shall be carefully controlled to prevent surging phenomena. These occur when the compressor is operated in off-design conditions at low flowrate and high compression ratio. Surging affects the whole machine and is aerodynamically and mechanically undesirable. It can cause vibrations, lead to flow reversal and seriously damage the compressor and the compressor driver and can negatively affect the whole cycle operation. To prevent surging, the compressor is controlled such as to remain at a distance from a surge limit line defined in a compression ratio vs. corrected flowrate diagram. A surge control line, also known as surge avoidance line, is usually set at a distance of the surge limit line and the compressor is controlled such that the operating point thereof remains within an operating envelope delimited by the surge control line. When the operating point of the compressor approaches the surge control line, the anti-surge control valve 15 is opened and gas is returned from the compressor delivery side 11 to the compressor suction side 9. Thus, the compressor operating point in a compression ratio vs. flowrate diagram is moved away from the surge control line and back in a safety operation area.
Re-circulating gas through the anti-surge line 13 causes power losses, since part of the gas which has been compressed in a power-consuming compression process is returned to the suction side of the compressor at the suction pressure. The corresponding power which has been spent to compress the recirculated gas flow is wasted.
A careful setting of the surge control line and a careful control of the compressor are desirable in order to prevent surging phenomena but at the same time avoiding recirculation of unnecessarily large amounts of compressed gas.
The process gas entering the compressor 3 at the suction side 9 thereof can be in dry conditions, i.e. containing no liquid volume fraction (LVF=0). In some operating conditions, however, the process gas can contain a significant amount of liquid phase. The liquid volume fraction LVF can be e.g. from about 0% to about 3%, which can correspond to a liquid mass fraction (LMF) from about 0% to 30%. It shall be noted that the upper limit is given by way of example only and shall not be construed as a limiting value.
During compression, the gas temperature increases and the liquid volume fraction can drop or even become zero. In some operating conditions, however, liquid can be present also in the gas flow at the delivery side 11 of the compressor 3.
It has been noted that if wet gas is processed, the surge limit line in a compression ratio vs. flowrate diagram moves from the right to the left as the liquid volume fraction LVF increases.
It would therefore be useful to determine, with a reasonable degree of approximation, the amount of liquid present in the gas flow, i.e. the LVF, since the surge control line could be moved towards the vertical axis of the compression ratio vs. flowrate diagram based on the actual LVF value, such that gas recirculation can be reduced.
In some circumstances the amount of liquid volume fraction (LVF) contained in the gas flowing through the compressor inlet 9 can be difficult to measure and such measurement may require costly and complex instrumentalities. In some situations, direct measurement of LVF may be unfeasible or inappropriate. As an alternative to direct measurement of LVF, an iterative method can be used to provide a sufficiently precise estimation of the actual liquid volume fraction, starting from easily measurable parameters of operation of the compressor 3.
An embodiment of the method will now be described reference being made to
More specifically,
wherein:
{dot over (m)} is the actual mass flow
Tin and Pin are the temperature and the pressure, respectively, at the suction side of the compressor;
zin is the compressibility factor or compression factor;
R is the gas constant (also known as the molar, universal, or ideal gas constant).
The corrected rotation speed of the compressor can be expressed as
wherein n is the angular speed and the other parameters are defined above.
In
wherein W is the actual power and the remaining parameters are defined above.
In some embodiments, the above defined corrected values can be rendered dimensionless by referring the actual measured pressure and temperature values to respective pressure and temperature reference values.
The curve W(LVF0) applies for dry gas, i.e. for LVF=0%. Curves W(LVF1), W(LVFj), W(LVFj+1), W(LVFj+2) are corrected power operating curves at increasing liquid volume fractions plotted as a function of the flowrate related parameter, e.g. the corrected mass flowrate {dot over (m)}C. Once again, the curves of
The curves C(LVFj) and W(LVFj) can be represented in form of tables or matrices of numeric values, wherein to each corrected mass flowrate {dot over (m)}C a compression ratio value (PR=Pd/Ps) and a power value (W) are associated. As stated above the curves further depend upon the rotation speed of the compressor and the gas composition. The curves plotted in
The curve SCL(0) in
Referring to the flow chart of
LVF(j)=0% j=0 (4)
i.e. it is assumed that the inlet gas is in dry conditions.
The actual compression ratio PRA=Pd/Ps can be calculated by measuring the delivery pressure Pd and the suction pressure Ps of the compressor 3 using pressure transducers 21, 25. Once the actual pressure ratio or compression ratio PRA has been determined, an estimated flowrate related parameter, e.g. an estimated corrected mass flowrate {dot over (m)}CE can be calculated using curve C(LVF0) in
Based on the estimated corrected mass flowrate {dot over (m)}CE, an estimated corrected power WEj required to drive the compressor can be determined using the curve W(LVF0) of
E
W=(WA−WEj) (5)
If the gas processed by the compressor 3 is actually approximately dry (i.e. LVF=0%, approximately), the error EW is around zero. An error EW outside a given range of tolerance around zero, e.g. defined by an error threshold EW0, indicates the initially assumed value of LVF (in the present example LVF(j)=0, dry gas conditions) is incorrect and a new value for LVF(j+1) must be used at the next iterative step (j+1).
An appropriate increased value can be selected, e.g. each subsequent LVF(j) value can be increased by an amount ΔLVF=0.01% over the previous one, which means that at each jth iterative step the tentative LVF value LVF(j) is set as
LVF(j+1)=LVF(j)+ΔLVF (6)
The above described sequence of steps of the iterative loop is then repeated with the newly set tentative value LVF(j) of liquid volume fraction. The C(LVFj) curve for LVF=LVF(j) is selected in the diagram of
E
W
=W
A
−W
E(j) (7)
is calculated and compared with the threshold EW0.
The iterative process thus described ends when an error EW on the estimated power related parameter is achieved, which is equal to or lower than the error threshold EW0. The tentative value LVF(j) to which the iterative process has converged is the estimated liquid volume fraction at the current operating conditions (current speed compressor and gas composition).
The value of LVF(j) thus determined can be used to select the optimal SCL. According to other embodiments, the SCL can be selected at each iterative loop, rather than only once the error EW has been minimized.
In the above described iterative method, two sets of operating curves have been used, namely the curves representing the compression ratio (PR=Pd/Ps) as a function of the flowrate related parameter {dot over (m)}C (
Also in this case the method can start by setting a tentative liquid volume fraction value LVF=0% and choosing the PR(W) curve corresponding to the dry gas operating conditions. Based on the measured actual compression ratio PRA=(Pd/Ps), the estimated power related parameter WEj can be calculated using the PR(W) curve corresponding to LVF=0%. The estimated power related value WEj is then compared with the actual power related value WA measured using the power transducer 31. A power error EW is then calculated as
E
W
=W
A
−W
E
j (8)
The error EW is compared with a threshold value EW0 and, if the error is greater than the admissible error threshold EW0, a next iterative step is performed, by setting a new tentative liquid volume fraction value
LVF(j+1)=LVF(j)+ΔLVF (9)
At each generic jth iterative step a tentative value LVF(j) is used to select the operative curve PR(Wj) corresponding to the set tentative LVF(j) value and the above described calculations are repeated, until the iterative process converges to an error EW that is equal to or lower than the error threshold EW0. The corresponding tentatively LVF(j) value is assumed as the estimated LVF.
A different embodiment of the method summarized in
LVF(i)=LVF(i)+ΔLVF (10)
as shown in the flow chart of
In all embodiments disclosed so far, a first compressor operating parameter and a second compressor operating parameter are used. According to the embodiment of
In the embodiment of
In the above described embodiments the starting point of the iterative process is LVF=0, i.e. the first iterative loop is performed assuming that dry gas is processed. This is convenient, since if the calculated error is above an error threshold, there is only one way of implementing the next iterative step, namely by increasing the assumed LVF value. However, in embodiments of the method described herein, the starting point of the iterative process can be any value for LVF. A sort of perturb-and-observe method can then be implemented. If the calculated error is above the admitted threshold, the assumed LVF is either increased or decreased. If the error calculated at the next iterative step is higher than the previously calculated error, the subsequent iterative loop will start by modifying the LVF in the opposite direction: it will be decreased if the previous iterative loop was executed by increasing the LVF value; otherwise, it will be increased, if the previous iterative loop was executed by decreasing LVF value.
In some embodiments, the LVF of the gas being processed can be estimated on the basis of thermodynamic calculations. The estimated value can be used as the starting point for one of the iterative methods disclosed above. Since in this case the estimated LVF value is different than zero, a perturb-and-observe iterative process can be used. The estimation of the starting LVF value is determined e.g. based on the gas composition, and upon the following parameters: suction pressure (Ps), delivery pressure (Pd), suction temperature (Ts) and delivery temperature (Td) of the gas processed by the compressor 3.
Since the operating curves change as a function of the rotation speed of the compressor, the rotation speed or the corrected rotation speed as defined by equation (2) can be used as a further parameter to select the proper family of operating curves each time the iterative process is performed.
The same holds true for the mean molecular weight of the gas. Different operating curves apply for different chemical compositions of the gas processed through the compressor 3. The chemical composition, and thus the molecular weight, of the gas is usually a slow-changing parameter. For instance, in case of gas wells, the composition remains quasi-constant and an update of the gas composition can be performed e.g. once a day or even less frequently. The gas composition can be analyzed off-line, e.g. in a laboratory using gas samples. Based on the result of the analysis the proper operating curves can be selected manually, for instance. On-line gas composition analysis can also be performed, e.g. by means of a gas chromatograph. The proper operating curves can be selected automatically. The mean molecular weight of the gas can be calculated based on the chemical composition.
The above described calculation methods can be performed continuously, or at a given frequency to monitor the actual LVF of the gas at the suction side of the compressor 3. For instance, the above described calculations can be re-started at given time intervals.
However, in order to render the above calculations more efficient, and to reduce the computational load, in some embodiments measures can be met to reduce the number of iterative calculations performed, or else to reduce the frequency wherewith these calculations are performed.
For instance, since the liquid volume fraction depends upon the suction pressure Ps and the suction temperature Ts of the gas, the other parameters (e.g. compressor rotation speed and gas composition) being the same, once the iterative method used has converged towards an error below the error threshold, the iterative calculation can be stopped. A new calculation to estimate the LVF can be performed only upon detection of a pressure or temperature fluctuation at the suction side 9 of the compressor 3. In other embodiments, the iterative calculations can be repeated periodically, but with a frequency that can be made dependent upon the fluctuation of the pressure and/or temperature at the suction side of compressor 3, i.e. the larger the fluctuations the more frequent the repetition of the iterative calculation.
In order to further simplify the method and reduce computational load, measures can be taken in order to perform the above described iterative calculation only if a preliminary routine establishes that wet gas is present at the suction side 9 of compressor 3. If the preliminary routine determines that dry gas is present at the suction side 9 of compressor 3, no estimation of the LVF is performed, since the actual value of the liquid volume fraction is zero.
A possible embodiment of a preliminary routine will be described here below, reference being made to the flow chart of
The first step of the preliminary routine provides for measuring the volumetric flowrate QVD at the delivery side of the compressor 3, e.g. by means of flowmeter 29. Based upon the measured temperatures Ts and Td at the suction side and delivery side of the compressor 3, upon the measured pressures Ps and Pd at the suction side and delivery side, as well as on the basis of the gas composition and assuming that dry gas is present at the suction side 9 of the compressor 3, an estimated mass flow rate is calculated. The estimated corrected mass flowrate ({dot over (m)}CS)E at the suction side 9 of the compressor 3 can then be calculated using equation (1). Based on the estimated ({dot over (m)}CS)E value and using the C(LVF0) curve of
E
PR
=PR
A
−PR
E (11)
is then calculated and compared with an error threshold EPR0. If the error EPR is equal to or lower than the error threshold EPR0, the assumption that the gas is dry at both the delivery side and the suction side can be assumed to be correct. Otherwise, if the calculated error EPR is above the error threshold EPR0, the conclusion is drawn that wet gas conditions are present at least at the suction side 9 of the compressor 3. In the first case (dry gas), no routine will be started to determine the actual LVF. In the second case, one of the routines for estimating the actual LVF, as summarized in
E
W
=W
A
−W
E (12)
is then calculated and compared with an error threshold EW0. If the error EW is equal to or lower than the error threshold EW0, the assumption that the gas is dry at both the delivery side and the suction side can be assumed to be correct. Otherwise, if the calculated error EW is above the error threshold EW0, the conclusion is drawn that wet gas conditions are present at least at the suction side 9 of the compressor 3. In the first case (dry gas), no routine will be started to determine the actual LVF. In the second case, one of the routines for estimating the actual LVF, as summarized in
In both embodiments of
Based on the above described method, the estimated liquid volume fraction LVF of the gas processed by the compressor 3 is determined with a sufficient accuracy and the estimated value can be used to select the optimal surge control curve SCL(LVF=x %) (
The above described method of LVF estimation can be used also for purposes different than surge control, whenever the liquid volume fraction of a wet gas shall be calculated.
The above described embodiments use methods for calculating the liquid volume fraction LVF of the gas processed by the compressor, e.g. in order to select a proper surge control line, in order to adapt surge control to the actual content of liquid in a wet gas.
The calculation methods described so far allow the LVF to be determined avoiding measurement of the actual liquid content at the suction side of the compressor. However, in order to adapt the surge control to the potentially variable content of the liquid in the gas, measurement of the LVF, rather than estimation thereof based on the above iterative calculation methods is not excluded.
While the disclosed embodiments of the subject matter described herein have been shown in the drawings and fully described above with particularity and detail in connection with several exemplary embodiments, it will be apparent to those of ordinary skill in the art that many modifications, changes, and omissions are possible without materially departing from the novel teachings, the principles and concepts set forth herein, and advantages of the subject matter recited in the appended claims. Hence, the proper scope of the disclosed innovations should be determined only by the broadest interpretation of the appended claims so as to encompass all such modifications, changes, and omissions. In addition, the order or sequence of any process or method steps may be varied or re-sequenced according to alternative embodiments.
This written description uses examples to disclose the invention, including the preferred embodiments, and also to enable any person skilled in the art to practice the invention, including making and using any devices or systems and performing any incorporated methods. The patentable scope of the invention is defined by the claims, and may include other examples that occur to those skilled in the art. Such other examples are intended to be within the scope of the claims if they have structural elements that do not differ from the literal language of the claims, or if they include equivalent structural elements with insubstantial differences from the literal languages of the claims.
Number | Date | Country | Kind |
---|---|---|---|
102016000070842 | Jul 2016 | IT | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2017/066978 | 7/6/2017 | WO | 00 |