METHOD FOR MONITORING A CORIOLIS MASS FLOW METER

20230184573 · 2023-06-15

    Inventors

    Cpc classification

    International classification

    Abstract

    The present disclosure relates to a method used to monitor a Coriolis mass flow meter, which has an oscillator with at least one measurement tube, the method including: exciting the oscillator so as to cause flexural vibrations of a first antisymmetric vibration mode by an excitation signal at a resonance frequency of the first antisymmetric vibration mode; sensing a vibration amplitude of the first antisymmetric vibration mode at the resonance frequency of the first antisymmetric vibration mode; sensing a time constant of the decaying free vibrations of the first antisymmetric vibration mode; and determining a modal elastic property of the oscillator with respect to the first antisymmetric vibration mode on the basis of the vibration amplitude of the first antisymmetric vibration mode, the excitation signal, and the time constant.

    Claims

    1-12. (canceled)

    13. A method for monitoring a Coriolis mass flow meter, the flow meter comprising an oscillator including at least one measurement tube, an exciter and two oscillation sensors, the method comprising: exciting the oscillator using the exciter as to cause flexural vibrations of a first antisymmetric vibration mode in response to an excitation signal applied to the oscillator at a resonance frequency of the first antisymmetric vibration mode wherein the exciter is offset from a center of the at least one measurement tube in a longitudinal direction, and wherein the exciter is offset not more than 5% of a measurement tube length in the longitudinal direction relative to the center of the at least one measurement tube such that symmetry breaking caused by an asymmetric mounting of that exciter is limited, which symmetry breaking causes a phase difference between vibrations of the two oscillation sensors upon excitation of the oscillator with an eigenfrequency of a symmetric drive mode; sensing a vibration amplitude of the first antisymmetric vibration mode at the resonance frequency of the first antisymmetric vibration mode; sensing a time constant of the decaying free vibrations of the first antisymmetric vibration mode; and determining a modal elastic property of the oscillator with respect to the first antisymmetric vibration mode based on the vibration amplitude of the first antisymmetric vibration mode, the excitation signal and the time constant.

    14. The method of claim 13, wherein the exciter is offset by at least 0.5% of the measurement tube length.

    15. The method of claim 13, wherein the exciter is offset by at least 1% of the measurement tube length.

    16. The method of claim 13, further comprising determining a modal quality of the oscillator based on the time constant, wherein the determining of the modal elastic property of the oscillator then is performed based on the vibration amplitude, the excitation signal and the modal quality.

    17. The method of claim 16, wherein the determining of the modal quality of the oscillator is further based on the resonance frequency of the oscillator.

    18. The method of claim 13, further comprising determining a variation in the modal elastic property of the oscillator by comparison with at least one reference value of the modal elastic property.

    19. The method of claim 13, wherein the mass flow meter is characterized by a calibration factor, which serves to determine a mass flow measurement value which is proportional to the calibration factor and to a time difference between in-phase points of the signals of two vibration sensors of the meter, and wherein the method further comprises adapting the calibration factor depending on a variation in the modal elastic property of the oscillator.

    20. The method of claim 19, further comprising: sensing a series of values of the modal elastic property; and determining a trend for the modal elastic property or a trend for the calibration factor.

    21. The method of claim 20, further comprising: determining a time period in which the modal elastic property or the calibration factor continues to lie within a permissible value range; and signaling the time period, and/or outputting an alarm signal, if the time period falls below a limit value, wherein the time period is not less than one week.

    22. The method of claim 21, wherein the time period is not less than eight weeks.

    23. The method of claim 13, further comprising: determining the modal elastic property of at least one further vibration mode; and calculating a current relationship between the modal elastic property of the first antisymmetric vibration mode and the modal elastic property of the at least one further vibration mode.

    24. The method of claim 23, further comprising evaluating the current relationship between the modal elastic property of the first antisymmetric vibration mode and the modal elastic property of the at least one further vibration mode.

    25. The method of claims 23, further comprising determining an extent and a type of wear of the at least one measurement tube of the oscillator based on the relationship between the modal elastic property of the first antisymmetric vibration mode and the modal elastic property of the at least one further vibration mode.

    26. The method of claim 23, wherein the at least one further vibration mode comprises the first symmetric flexural vibration mode and/or the second symmetric flexural vibration mode.

    27. The method of claim 13, wherein the modal elastic property comprises a modal flexural stiffness, a modal flexibility or a calibration factor.

    Description

    [0023] The invention is now explained in more detail on the basis of the exemplary embodiments shown in Figures.

    [0024] The following are shown:

    [0025] FIG. 1a: a schematic representation of an exemplary embodiment of a Coriolis mass flow meter for implementing the method according to the invention;

    [0026] FIG. 1b: a schematic representation of electromechanical transducers of the exemplary embodiment of the Coriolis mass flow meter from FIG. 1a;

    [0027] FIG. 2: a diagram of the vibration modes of a Coriolis mass flow meter;

    [0028] FIG. 3: a flow chart of a first exemplary embodiment of the method according to the invention;

    [0029] FIG. 4: a diagram of the relationship between variations of the calibration factor and the modal flexibility of the drive mode; and

    [0030] FIG. 5: a flow chart of a second exemplary embodiment of the method according to the invention.

    [0031] The Coriolis mass flow meter 1 shown in FIG. 1a comprises an oscillator 8 with curved measurement tubes 10 traveling substantially in parallel, as well as an exciter 11 which acts between the measurement tubes 10 in the direction of flow in order to excite them so as to excite these to flexural vibrations counter to one another. The exciter 11 is arranged offset, in the longitudinal direction of the measurement tubes, by approximately 2.5% of the length L of the measurement tubes relative to the center of the measurement tubes. Upon excitation of the oscillator with the exciter 11, a sufficient asymmetric force therefore acts in order to excite the first antisymmetric vibration mode, what is known as the f2 mode or first Coriolis mode, to resonant vibrations if the excitation of the oscillator takes place with a resonance frequency of the first antisymmetric vibration mode. Furthermore, the Coriolis mass flow meter 1 has two vibration sensors 12.1, 12.2, which are arranged symmetrically in the longitudinal direction relative to the center of the measurement tubes 10 in order to detect the relative movement of the measurement tubes 10 vibrating counter to one another. The measurement tubes 10 extend between two flow dividers 16, which fluidically combine the measurement tubes 10 and are respectively connected to a flange 18, which serves for the installation of the Coriolis mass flow meter 1 in a pipeline. A rigid support tube 60 which connects the flow dividers to one another extends between said flow dividers 16 in order to suppress vibrations of the flow dividers 16 counter to one another in the frequency range of the flexural vibration modes of the measurement tubes 10 counter to one another. The support tube can furthermore carry an electronics housing 80 in which a measuring and operating circuit 77 is contained which is configured to operate the meter and to implement the method according to the invention.

    [0032] As shown in FIG. 1b, the exciter 11 and the vibration sensors 12 are especially designed as electrodynamic transducers that respectively have an excitation magnet 14 or sensor magnet 14.1, respectively, and an excitation coil 13 or sensor coil 13.1, respectively, which are mechanically connected opposite one another to one of the measurement tubes 10. The excitation coil 13 is configured to be supplied by the operating circuit 77 with an alternating current whose frequency corresponds to the instantaneous eigenfrequency of a flexural vibration mode to be excited. The resulting magnetic field alternately effects an attractive and repulsive force on the excitation magnet 14, whereby the measurement tubes 10 are set into vibration counter to one another. Accordingly, the relative movements of the sensor magnets 14.1 vibrating with the measurement tubes 10 relative to the sensor coils 13.1 induce a voltage in the sensor coils 13.1, which depends especially on the relative velocity of the measurement tubes relative to each other. The measuring and operating circuit 77 is configured to sense and evaluate the induced voltages in order to determine therefrom the relative velocities or the deflection of the vibration sensors 12.1, 12.2 or of the measurement tubes 10, a modal deflection of the measurement tubes 10 for different vibration modes.

    [0033] The mode-dependent deflection of a measurement tube is shown schematically in FIG. 2. The curve a.sub.A hereby shows the bending line of a measurement tube for the first symmetric vibration mode, which is also called the drive mode or f1 mode. The curve a.sub.C1 shows the bending line of the measurement tube for the first Coriolis mode or the first antisymmetric vibration mode, in which the measurement tube is deflected by the Coriolis forces if a mass flow flows through the measurement tube vibrating with the first symmetric vibration mode. The first antisymmetric vibration mode has a vibration node in the tube center at z=0 in the longitudinal direction of the measurement tube. An exciter at this position would not be able to excite a vibration of the first antisymmetric vibration mode. Therefore, here the exciter is offset from the center by approximately 2.5% of the measurement tube length, i.e., approximately 5% of half the measurement tube length. The measurement tube length is hereby the length of a measurement tube center line, following the curved course of a measurement tube, between the inlet-side and outlet-side flow dividers 16 in which the measurement tubes 10 are fixed by their ends. In the offset position, the exciter can excite the first antisymmetric vibration mode if it impresses an excitation force F.sub.E at the resonance frequency of the first antisymmetric vibration mode.

    [0034] The positions of the vibration sensors are selected symmetrically in the longitudinal direction with respect to the measurement tube center of the measurement tubes, such that the deflections X.sub.S1, X.sub.S2 of the vibration sensors produce a sufficient measurement signal given both vibrations in the drive mode and the first antisymmetric vibration mode. Furthermore, shown in FIG. 2 is the bending line a.sub.C2 for the second antisymmetric vibration mode or the second Coriolis mode in which the measurement tube is deflected if the measurement tube through which a mass flow passes vibrates in the second symmetric drive mode (not shown here), the f3 mode. Similarly, the second Coriolis mode can be excited if the exciter impresses an excitation force F.sub.E at the resonance frequency of the second Coriolis mode.

    [0035] Due to the high modal quality Q of between 1000 and 10000, for example, the amplitudes of the vibration modes of the oscillator or of its measurement tubes exhibit a strong resonance exaggeration. In order to be able to infer the modal stiffness or flexibility of the individual modes, the modal quality Q is also to be determined in addition to the vibration amplitudes at the respective resonance frequencies. For this purpose, especially a decay curve of the respective vibration mode can be sensed after the excitation force has been switched off. The vibration amplitude normalized with the quality Q and the excitation force F.sub.E is a measure of the modal flexibility.

    [0036] The modal flexibility of the first antisymmetric vibration mode impresses a calibration factor calf which, in a first approximation, is inversely proportional to this modal flexibility, and which relates the mass flow rate dm/dt to a time delay Δt between zero crossings of the two vibration sensors, i.e.:

    [00001]dm/dt=calfΔt

    [0037] A monitoring of the modal flexibility of the first antisymmetric vibration mode with the method according to the invention thus directly enables monitoring and correction of the calibration factor calf, or a validation of the mass flow measurement value dm/dt.

    [0038] The method steps according to a first exemplary embodiment 100 of the method according to the invention are explained using FIG. 3. The method 100 can, for example, be implemented continuously, periodically, or in an event-controlled manner, wherein a triggering event can be, for example, a user request or the determination of a change in another monitoring variable of the measuring device.

    [0039] The method 100 begins with the excitation 110 of the oscillator to flexural vibrations of a first antisymmetric vibration mode with a modal excitation signal F.sub.c1 at a resonance frequency ω.sub.c1 of the first antisymmetric vibration mode. This first antisymmetric vibration mode is the first Coriolis mode or f2 mode, as explained in conjunction with FIG. 2.

    [0040] In the steady state of this first antisymmetric vibration mode, the sensing 120 of the vibration amplitude Xc.sub.1 of the first antisymmetric vibration mode takes place at its resonance frequency. For this purpose, the velocity-proportional induction voltage of the electrodynamic vibration sensors is evaluated at the resonance frequency of the first antisymmetric vibration mode.

    [0041] This is followed by the sensing 130 of a time constant τ.sub.c1 of the decaying free vibrations of the first antisymmetric vibration mode, for which purpose the excitation signal at the resonance frequency of the first antisymmetric vibration mode is partially or completely switched off, and the decaying induction voltage amplitudes of the vibration sensors are sensed at the resonance frequency.

    [0042] Finally, the determination 140 of a modal elastic property of the oscillator with respect to the first antisymmetric vibration mode takes place on the basis of its vibration amplitude, the excitation signal, and the time constants. For this purpose, for example, a modal quality Q.sub.c1 can initially be determined on the basis of the time constants. The modal quality Q.sub.c1 can, for example, be determined as follows:

    [00002]Qc1=τc1ωc12,

    where ω.sub.c1 is the resonance frequency of the considered vibration mode.

    [0043] The determination of the modal elastic property of the oscillator then takes place on the basis of the vibration amplitude, the excitation signal, and the modal quality.

    [0044] The modal elastic property can be, for example, the modal flexibility N.sub.c1, which is proportional to the modal vibration amplitude X.sub.c1 divided by the modal quality Q.sub.c1 and the amplitude of the modal excitation signal F.sub.c1, i.e.: [0045] N.sub.c1 = K.sub.c1 .Math. X.sub.c1 /(F.sub.c1 .Math. Q.sub.c1), where K.sub.c1 is a mode-specific constant.

    [0046] By comparison 150 of the modal flexibility N.sub.c1 determined in this way with a reference value N.sub.c1-0, a variation in the modal elastic property of the oscillator can be determined, wherein the reference value represents, for example, the state upon startup of the mass flow meter.

    [0047] As mentioned above, the calibration factor calf of the mass flow meter is substantially inversely proportional to the modal flexibility N.sub.c1. In this respect, the calibration factor calf is also available as an elastic property to be monitored of the first antisymmetric vibration mode, wherein the calibration factor calf can be determined as follows: calf = K.sub.calf / N.sub.c1, where K.sub.calf is a device-specific proportionality factor.

    [0048] The adaptation 160 of the calibration factor calf depending on a variation in the modal flexibility N.sub.c1 furthermore enables precise mass flow measurements even given wear of the measurement tubes. After repeated adaptation of the calibration factor given modified modal flexibilities N.sub.c1, a trend analysis 170 of the calibration factor calf can furthermore take place, and a prediction of remaining service life 180 can be provided, relating to the point in time up to which the meter can still be operated, assuming the same media properties. Details in this regard are explained further below in conjunction with FIG. 4.

    [0049] In addition, the method according to the invention can advantageously be combined with the method for monitoring the state of measurement tubes according to the international publication WO 2012 062551 A1, which teaches the monitoring of the modal flexibility N.sub.a of the first symmetric drive mode. This can especially be determined independently of quality via excitation outside of resonance. The relative deviation ΔN.sub.a of a current modal flexibility N.sub.a of the first symmetric drive mode from a reference state N.sub.a-0, for example in the brand-new state or after a certification, is likewise an indicator of a variation in the measurement tube. FIG. 4 relates the relative deviation ΔN.sub.a of the modal flexibility of the first symmetric drive mode to the relative deviation Δcalf of the calibration factor. Test series have yielded that, depending on the cause of the variation in the measurement tubes, two separate regimes occur for the relationship of the two monitoring variables Δcalf and ΔN.sub.a. Both regimes have the common starting point of a measurement tube in the reference state at (0,0). In the event of corrosion, the wall of a measurement tube is substantially uniformly attacked, such that the moments of inertia of all measurement tube cross sections change uniformly. Accordingly, the modal stiffnesses of the relevant vibration modes are affected uniformly so that, given corrosion, a very good correlation is to be observed between Δcalf and ΔN.sub.a. This corrosion regime is shown cross-hatched below the line b in FIG. 4. By contrast, given abrasion another picture emerges. Abrasion usually arises in heterogeneous media that comprise a liquid with a solid load. Depending on the Reynolds number, concentration, density distribution, and size distribution, different spatial distributions of the abrasion can occur, wherein, during its genesis, linear trajectories within the region shown in dotted lines in FIG. 4 above the line a in FIG. 4 were observed in a first approximation for the monitoring variables Δcalf and ΔN.sub.a. Since the two regimes are distinctly separate from one another, a plant operator is thus provided with a means to recognize abrasion and corrosion processes in the measuring device at an early stage, and to identify the type of material erosion, on the basis of the relationship of Δcalf and ΔN.sub.a. Thus, if a pair of values (Δcalf, ΔN.sub.a) is above the line a, abrasion is to be assumed, whereas if it is below the line b, this indicates corrosion. FIG. 5 shows a flow chart of a second exemplary embodiment 200 of a method according to the invention which realizes this aspect of the invention.

    [0050] The first method steps proceed analogous to the first exemplary embodiment, up to the determination 240 of the calibration factor calf on the basis of the amplitude Xc.sub.1 of the first antisymmetric vibration mode, its decay time τ.sub.c1, and the associated excitation signal F.sub.c1 at the resonance frequency ω.sub.c1. The determination 250 of the modal flexibility N.sub.a of the first symmetric vibration mode takes place in parallel with this. This can take place analogously to the determination of the modal flexibility in the first exemplary embodiment in resonance, or independently of quality with excitation outside of resonance, as described in WO 2012 062551 A1. In fact, the measurements for determining the calibration factor calf and the modal flexibility N.sub.a of the first symmetric vibration mode can take place simultaneously, since the vibrations can be excited in a superposed manner. If current values for the calibration factor calf and the modal flexibility N.sub.a are obtained, the formation 260 of a relationship between the relative deviations Δcalf and ΔN.sub.a takes place from their respective reference values, wherein the relative deviations of a variable x are determined in accordance with Δx = (x - x.sub.ref) / x.sub.ref, where x is the calibration factor calf or the modal flexibility N.sub.a, and where x.sub.ref refers to the state of the respective variables upon startup of the meter. The evaluation 270 of the relationship then takes place in the form of a classification, wherein it is determined whether the relationship indicates a corrosion or abrasion. A classification can only be reliably implemented when the determined wear has already reached a certain extent, for example if the value pairs in the illustration of FIG. 4 lie outside the inner elliptical arc i. Using the duration which requires a trajectory of the value pairs from the inner elliptical arc i to the middle elliptical arc ii, it can then be extrapolated when the value pairs reach a critical wear limit, which is represented by, for example, the outer elliptical arc iii. This point in time can be provided as a notice to plan a maintenance measure. Moreover, an alarm can be generated if the time until reaching the critical wear limit falls below a limit value of, for example, a quarter and/or a month.

    [0051] In the exemplary embodiments, the modal elastic property of the oscillator was described as the modal flexibility of the oscillator or the measurement tubes. Of course, the modal flexural stiffness or the calibration factor calf can similarly be used to model or describe the wear.

    [0052] In addition to the first symmetric flexural vibration mode, the second symmetric flexural vibration mode, which is also referred to as an f3 mode, can also be used as a further mode.