SIMULTANEOUS DETERMINATION OF MULTIPHASE FLOWRATES AND CONCENTRATIONS

Abstract: The present invention relates to a method and a device for measuring volume flow rates of liquid phase components and gas and determinating their volume concentrations in a multiphase mixture along a pipeline. Measurements are executed with an ultrasonic system which includes a set of local acoustic transducers arranged in the interior of the pipeline. Each pair of an emitter and a receiver of the transducer forms a sampling volume of a medium being under control. Volume concentrations of mixture components are determined by timing of passage of acoustic pulses through the sampling volume of the medium. Volume flow rates of the mixture components are calculated by measuring phase velocities and volume concentrations in two pipeline divisions with different cross-section areas located in series at a distance one from the other in flow direction.

Description: Title: Simultaneous determination of multiphase flowrates and concentrations . FIELD OF THE INVENTION The invention relates to a method and a device for determination volume flow rates and volume concentrations of liquid phase components and gas in multiphase liquid/gas mixtures such as oil well fluids, composed of oil, water and gas. DESCRIPTION OF THE RELATED ART An effluent flowing along a pipeline from an oil well is a multiphase mixture of oil, water and gas. Accurate and simultaneous measurement of flow rates and volume concentrations of mixture components is important for control of the well operation. Methods and devices for measuring these flow characteristics accepted up to the present demand preliminary separation of gas in special separators which are installed in measuring devices on the territory of oil fields. This fact leads to considerable capital expenses by execution of such measurements. Techniques for measuring of the said multiphase flow characteristics without a prior separation of gas phase are also known. These methods and devices are based on various physical principles: difference of density and electromagnetic characteristics of the components, interaction with gamma-rays and ultrasonic waves, and others. US patent 5.287.752 describes a device for determination flow rates of multiphase fluids by means of a set of capacitors placed on two parallel plates which are arranged inside a horizontal or inclined pipeline parallel to a flow direction. In order to determine a water/oil volume proportion and a part of pipeline cross-section occupied by liquid phase impedances of a medium located at the moment in measuring cells of all elementary capacitors are measured. A velocity of the liquid phase is determined by measuring and cross-correlation impedances of elementary capacitors located in a matrix row situated in a part of cross-section occupied by the liquid phase. A velocit_ of gas is determined by measuring passage time of structural deformations of the flow in the upper part of the pipeline. Volume flow rates of the phases are determined taking into account the parts of the pipeline cross-section occupied by the liquid and gas phases of the flow. The proposed method has limited sphere of application since it can be effectively used only by an intermittent flow regime. Besides a type of emulsion and dispersion of components aren' t taken into account in this method. US patent 5.367.911 describes an apparatus for sensing fluid behaviour in a conduit that can be used as a flow meter. The measuring device includes at least two sensors arranged within a pipeline, one downstream the other. The sensors may include acoustic transducers or electrical conductivity (or resistivity) sensors. Each sensor provides an output data signal, indicative of the measured physical property of a medium flowing in respective sampling volumes. Output signals are processed in a circuit and are cross-correlated. Since a distance between the sensors is known a calculation of the flow speed is made. However the authors of the patent don' t take into account that a gas phase moves relative to a liquid phase in multiphase flows. European patent A 0684458 relates to a method and a device for measuring flow rate of multiphase fluids. The device comprises two sections situated at a distance one from the other along a duct. The sections are implemented in the form of Venturis. Each section includes a passage with different cross-section areas provided with means for inducing a change of speed therein and a variation of dynamic pressure correspondingly. Besides the device includes appropriate means for measuring the resulting pressure differences. Two pressure difference signals obtained in the said respective sections are suitable for cross-correlation to produce a third signal representative of a total volume flow rate. In order to determine the flow rates of phases another pressure difference is measured in a venturi type flow section and the signal that is a function of the total mass flow rate Q and density p of a mixture is obtained. One more pressure difference is measured in the section of a "gradiomanometer" type. This section is disposed in a portion of a vertical duct that has constant cross-section area. In a conventional manner the latter pressure difference represents a static pressure difference that is proportional, to a first approximation, to the density p of the mixture. In addition, the device includes a module situated in the duct, that provides one or more indications of composition of the multiphase fluid, in other words, determines volume or mass concentrations of components constituting the fluid. A processor calculates the mass flow rate Q via parallel paths depending on gas content. When the gas content is moderate (< 65%) a first processing way is used on the basis of a signal proportional to the expression Q 2 /p and a signal representative of the density p of the mixture. A second processing way based on a signal representative of the total volume flow rate and a signal representative of the density p of the mixture is used when the gas content is high (> 65%) . It is necessary to note that devices with venturi type sections have a small dynamic range and hence may be used in a limited range of flow rate measurements. Besides, a flow rate factor of such devices significantly depends on gas content that influences on the measurement accuracy. In the described device the gas content is calculated but not measured that also decreases the measurement accuracy. Using the device for measuring the flow rate of oilwell effluents containing crude oil may lead to clogging of its pressure takeoffs . SUMMARY OF THE INVENTION The present invention provides a method and a device for determination volume flow rates of multiphase mixture components along a portion of a pipeline without prior separation of gas. The present invention provides measuring volume concentrations of the multiphase mixture components along the portion of the pipeline. The invention also provides a method and a device for measuring the above-mentioned characteristics of a multiphase medium with different types of flow. The invention ensures obtaining reliable data of measuring characteristics of effluents with different sizes of gas inclusions. Besides the invention secures compactness of the device and its simple portability. The above-mentioned features are achieved through the method for determination volume flow rates of liquid components and gas in the multiphase mixture flowing along a pipeline according to which a measuring flow cell installed into a pipeline includes two pipeline divisions located in series in the flow direction and having different flow section areas: F 2 = k F : (diameter D = Di k) k ≠ 1. When k « 0,5 changing of the flow section area induces significant changing of a liquid phase velocity and correspondingly of a real gas phase velocity in the measuring pipeline divisions ( w g ,ι < w qr 2 ) while changing of a relative velocity of gas inclusions and of a real volume gas concentration φ in the mixture is insignificant. Calculation analysis of a mixture flow model has allowed to derive a formula for determination a liquid phase volume flow rate for the multiphase flow passing through the calibrated pipeline divisions: Q, = (k / 1 - k) Fi [ w q , 2 (1 - φ 2 ) - w g ,ι (1 - ] . A gas volume flow rate is determined by the following formula: Q g = Fi Wg, i • φi = F 2 W g , • φ2 • The real velocities of gas phase w q , the volume concentrations of gas φ, the volume concentrations of liquid components, such as water W and oil (1-W) in the calibrated pipeline divisions are determined by means of ultrasonic sensing sampling volumes of the multiphase flow with a set of transducers arranged within the measuring pipeline divisions along a radius of a flow section. These transducers serve as emitters and receivers of acoustic signals in sampling volumes. Obtained values of local characteristics of the multiphase flow are then averaged over cross-section areas of the measuring pipeline divisions. Measuring of the real gas velocity is executed by cross- correlation of sensor signals or by Doppler' s method. Measuring of the volume concentration of gas is executed through indication of acoustic conductivity of the sampling volumes of the medium. Ultrasonic measuring of the volume concentration of liquid phase components is based on determination of time of acoustic pulses pass through the sampling volume since it was found that in a fluid such as water/oil mixture the time of signal pass depends practically linearly on a proportion of the volume concentrations of these components regardless emulsion type. The mentioned features are also provided by the apparatus for determination volume flow rates and volume concentrations of liquid components and gas of a liquid-gas multiphase mixture flowing along a pipeline including a measuring flow cell being installed in a pipeline. The measuring flow cell comprises two pipeline divisions disposed in series in a flow direction and having different flow section areas: F 2 = k Fi (diameter D 2 = D k) k ≠ 1. Changing of a flow section area (when k « 0,5) induces significant changing of a liquid phase velocity and a real gas velocity in the measuring pipeline divisions ( w q < w g , 2 ) while changing of a relative velocity of gas inclusions and of a real volume gas concentration φ in a mixture is insignificant. A volume flow rate of a liquid phase is determined by a difference of products of the real gas phase velocity w q by a part of a pipeline section occupied by a liquid phase (1 - φ) in the first and the second measuring pipeline divisions: Q t = (k / 1 - k) Fi [ w q , 2 (1 - φ 2 ) - w q , ! (1 - ] . Gas volume flow rate is determined by the following formula: Qg = i Wq, 1 • φl = F 2 W q , 2 ■ φ 2 . The real velocities of gas phase w q , the volume concentrations of gas φ, the volume concentrations of liquid components, such as water W and oil (1-W) in the calibrated pipeline divisions of the pipeline are determined by means of ultrasonic sensing local volumes of the multiphase flow with a set of transducers arranged within the measuring pipeline divisions along a radius A a flow section. The principle of operation of a local gas velocity meter is based on determination a cross-correlation function of an amplitude of a signal of an acoustic conductivity transducer. The two transducers are placed at a fixed distance one upstream the other in the flow direction. The acoustic transducer consists of an emitter and a receiver of ultrasonic pulses providing acoustic illumination of the sampling volume. The transducer may be used as an emitter and a receiver of reflected signals in the "emission-reception" mode. An electro-acoustic channel of the meter operates in the following way: voltage pulses from a pulse generator come to the emitter where they are converted into ultrasonic pulses. After passing through the sampling volume they are received by the receiver, converted into the voltage pulses, amplified and are sent to a peak detector input which is controlled by strobe pulses. The strobe pulses determine a time interval during which a signal reception is expected. A voltage at the input of the peak detector is proportional to an amplitude of a received signal and is determined by acoustic energy losses in a sensor sampling volume. Output signals of peak detectors come to a calculator which determines a cross-correlation or an autocorrelation function (in case of one transducer) and calculates the real local velocity of gas phase or liquid phase without gas. Besides the expounded principle Doppler' s method for measuring the gas phase local velocity through sensing a medium by ultrasonic pulses directed upstream the flow may be used. In this variant emitter and receivers are also disposed inside the measuring pipeline divisions. The principle of operation of a volume gas concentration meter is based on indication of acoustic conductivity of a sampling volume. A signal from a voltage pulse generator is sent to an emitter consisting of a transmitter and a waveguide. After conversion acoustic pulses reach the sampling volume through the waveguide, pass through the volume and through a receiving waveguide come to the transmitter where they are transformed into a voltage signal that after amplification comes to a peak detector. A strobe pulse former opens the peak detector for a time while the signal coming is expected. From the peak detector an output signal proportional to an amplitude of the received signal comes to a comparator which compares the output signal of the peak detector with a discrimination level set by a former of discrimination level. An output signal of the comparator comes to a calculator which determines the volume gas content in a medium as a ratio of time of gas phase presence in the sampling volume to the full time of the measurement. The principle of operation of ultrasonic meter of volume concentrations of liquid components is based on determination time of passage of ultrasonic pulses through a sampling volume of a multiphase flow since it was found that in liquid phase, such as water/oil mixture the time of passage of ultrasonic signal practically linearly depends on relationship of volume concentrations of liquid components regardless emulsion type. A distance between emitter and receiver is chosen so that penetration of large gas inclusions with sizes more than 1 mm is prevented. Voltage pulses from a generator are sent to the ultrasonic emitter which forms acoustic pulses. The acoustic pulses pass through a sampling volume, are received by the receiver and are transformed into a voltage signal which is amplified and then sent to a comparator being strobed. The comparator opens for the time while the reception of the signal is being expected by means of the strobe pulse former that ensures high interference immunity of the scheme. Simultaneously with forming the emitting pulses a scheme forming a pulse duration is activated. This scheme is stopped by the signal coming from the comparator output. So the duration of the output signal is equal to the time of the ultrasonic signal pass from the emitter to the receiver. Then the pulse is transformed into an amplitude signal and comes into calculator which determines volume concentration of liquid phase components. A processor functioning according to set programs controls operation of the meters of local flow parameters w q , w 9ι 2 , Φi, φ 2> W, averages these parameters over the cross-sections of measuring pipeline divisions and calculates volume flow rates of liquid phase components and gas . BRIEF DESCRIPTION OF FIGURES AND DRAWINGS Figure la, lb. Configuration of ti multiphase flow rate device proposed within the framework of the present invention. Figure 2. Block diagram of the gas phase local velocity meter for the variant of sequential location of transducers. Figure 3. Voltage diagram of the signal processed in the block diagram presented in Figure 2. Figure 4. Typical form of a cross-correlation function for output signals of peak detectors. Figure 5. Block diagram of the local gas phase velocity meter for the variant of transducers located in series ("emission- reception" mode) . Figure 6. Voltage diagram of the signal processed in the block diagram presented in Figure 5. Figure 7. Block diagram of the local gas phase velocity meter while working of a pair of transducers in the mode of autocorrelation function forming. Figure 8. Typical form of an autocorrelation function for output signals of peak detectors. Figure 9. Block diagram of the local gas phase velocity meter while working of one transducer in the "emission-reception" mode. Figure 10. Voltage diagram of the signal processed in the block diagram presented in the Figure 9. Figure 11. Autocorrelation function form of the of the output signal of the peak detector in the variant of using one transducer. Figure 12. Block diagram of the ultrasonic Doppler' s meter of the gas phase local velocity. Figure 13. Typical form of signals in the block diagram of ultrasonic Doppler' s meter of velocity presented in the Figure 12. Figure 14. Second variant of the block diagram of the ultrasonic Doppler' s meter for measuring of the gas phase local velocity. Figure 15. Voltage diagram of a signal processed in the block diagram presented in Figure 14. Figure 16. Block diagram of the gas content meter in a multiphase mixture is presented. Figure 17. Diagram of a signal processed in the block diagram of the gas content meter. Figure 18. Diagram of an output signal of the peak detector. Figure 19. Block diagram of the ultrasonic meter of volume concentrations of liquid components. Figure 20. Voltage diagram of a signal processed in the block diagram for measuring liquid phase volume concentrations. Figure 21. Second variant of the block diagram of the ultrasonic meter of volume concentrations of liquid components. Figure 22. Voltage diagram of a signal for the second variant of the block diagram of the volume concentration meter of liquid components. VARIANTS OF THE INVENTION REALISATION Configuration of a measuring flow cell of the apparatus for determination volume flow rates of liquid components and gas in a multiphase mixture is represented in Figure la and Figure lb. The measuring flow cell is installed into a pipeline by means of flange connections. The flow cell comprises two measuring pipeline divisions 1 and 2 disposed in series in a flow direction and having different flow section areas: F 2 = k F λ (diameter D 2 = ϋ k) . For Figure 1 k < 1 holds. Changing of a flow section area induces significant changing of a liquid phase velocity and a real gas phase velocity in the measuring pipeline divisions with the cross-section areas F x and F 2 . In order to ensure minimum hydrodynamic disturbance of the flow a transition from a first division to a second division and back to the initial cross-section area Fi of the pipeline is realized through transition pipeline divisions 3 and 4. Each sensor 5 and 6 of a real velocity meter and a gas content meter includes a set of transducers placed inside the measuring pipeline divisions along radii of the sections. A sensor 7 of the volume concentration meter of liquid components comprises a set of transducers located in a cavity of the first pipeline divisions. To accelerate the process of viscous medium change within the transducer's volumes and for removal of paraffin deposits the sensors are equipped with mechanical cleaning arrangements or electric heaters . The sensors are installed in such a way that they can be removed from the measuring pipeline divisions, for example, for execution of technical maintenance or substitution. Let us consider separately the meters of multiphase flow parameters being a part of the dev ce and a calculation model of the multiphase mixture flow used for determination of volume flow rates of mixture components. A calculation model of a liquid/gas mixture where gas inclusions of various size represent gas phase is used to determine flow rates of phases. Averaged physical values are used in formulas. The real volume concentration of the gas in the i-th cross- section of a flow is: φ, =F g , 1 / F x , (1) where F x = π/4 D is a cross-section area of the i-th pipeline divisions, F g(1 = Φi • i is a cross-section area occupied by gas. Since F v = F g>1 + F (>1 , where F ((1 is a cross-sectional area occupied by liquid, we can write instead of formula (1): t g,! + /,,! * Wg,_ / w, t x where w r q> 1 = Q gιl / F x is a reduced velocity of the gas phase in the i-th pipeline divisions, where Q, gιl is a volume flow rate of the gas phase in the i-th pipeline divisions; w 1 ! , ! = Q»,ι / F L is a reduced velocity of the liquid phase in the i-th pipeline divisions, where Q, #1 is a volume flow rate of the liquid phase in the i-th pipeline divisions; W g , ! = Q g , ι / Fg (1 is a real velocity of the gas phase in the i-th pipeline divisions, where Q gιl is a volume flow rate of the gas phase in the i-th pipeline divisions; w t f = Q t , ι I F (ιl is a real velocity of the liquid phase in the i-th pipeline divisions, where Q (rl is a volume flow rate of the liquid phase in the i-th pipeline divisions; and F, (1 = (1- φ Fi is a section area of the i-th pipeline divisions occupied by liquid. Since besides w (/1 = w r tιl / (1 - φ , and w qιl = w lι∑ + w g,1 rel , where w q, f el is a relative velocity of the gas phase in the i-th pipeline divisions, so we have: W'n w r π + t ,, x + φ q rel According to experimental data relative velocity of the flow of gas bubbles Wg #1 rel (group velocity of floating-up) is connected with the real volume concentration φ x by the following relationship: w q ,f el = w g ,. / (1 - q>α ) , (4) where w g , ∞ is average velocity of the individual bubble floating-up in the infinite liquid medium. Real velocities w q and w qι2 in measuring pipeline divisions are connected with relative velocities in the following way: w g ,ι = w t ,ι + Wg, ! rel and w g , 2 = w, ι2 + w q< X el . (5) Subtracting the first equality (5) from the second one (5) we get the following equality: W g , 2 - W g(1 = ΔW g = (W,, 2 - W tΛ ) - (W g , 2 rel - W q , ! " 1 ) , (6) which can be written in the form: Δ w g = vf lι2 /{l - φ 2 ) - w r tιl /{l - φ : ) + w g,∞ [1/(1 - φ 2 ) - 1/(1 Supposing that we have the following relationships F 2 = kF x , where k ≠ 1, and taking into account that w r ltl = Q, #1 /F we get: Δ w g = [ 1] + [ 1] (8) Fi (1 - φi) k (1 - φ 2 ) 1 - Φi 1 - φ2 since Q ( i = Q ( . It follows from the relationships (3) and (4) that Φi = w r q , 1 I X q , x + tf {ι + w g , „) . (9) After execution of appropriate transformations and taking into account that Q gιl = Q g we get: 1/ φ, = 1 + Q, /Q g + F 1 w qι∞ /Q g . (10) Substituting the value Q g = F x φ x w qιX , where w gιl and φ x - are values being measured we obtain: 1 1 Q ( q> m = 1 + ( + ) , (11) Q / Wg, «, where from φ = 1 - - . (12) F w q , x w qtX Consequently Q, = F x [ w qι (1 - φ x ) - w q , . ] . (13) It should be noted that in the case of stationary liquid (Q ( = 0) it follows from the formula (13) the following relationship: w q , x = w g , ∞ /(l - φ , what coincides with definition of the relative velocity (4) so that in this case t^ g/1 = q ,f el . It follows from the formula (13) that O f = X [ W g , l (1 - φi) - W g , . ] (14) and Q( = F 2 [ w g , 2 (1 - φ 2 ) - w q , ∞ ] . (15) Having equated relationships (14) and (15) and taking into account that F 2 = kFi, where k ≠ 1, we get: Q« = Fi [ W g , 2 (1 - φ 2 ) - w q (1 - φj) ] k/(l - k) . (16) So the volume flow rate of liquid phase in the calibrated pipeline divisions is determined according to the relationship (16) from the measured real velocities and volume concentrations of gas phase in the first and the second measuring pipeline divisions. If F z = 0,5Fι the expression (16) becomes: Q, = Fi [ w q , 2 (1 - φ 2 ) - W g , ! (1 - φ x ) ] . (17) Besides it should be noted that if φi = φ 2 = 0 velocity meters shall fix acoustic inhomogeneity of liquid phase and correspondingly velocities w, #1 and w t ι 2 . So the relationship (16) transforms into the formula Q, = F λ ■ w l , and if φ x = φ = 1 this relationship takes the form Q, = 0. Volume flow rates of liquid phase components are determined by the formulas: Q o u = Q," (1 - ) and Q w = Q, • W , (18) where W is the volume concentration of water in the emulsion. Volume flow rate of the gas phase is determined by the following relationship: Qg = w a Φi w a φ 2 (19) Figure 2 shows a block diagram of an ultrasonic meter of local gas velocity w q of a multiphase mixture. A circuit of the meter includes: a generator of voltage pulses 8, a first transducer 9 connected in series with the generator 8 and comprising an emitter 10 and a receiver 11 (a gap between them forms a first sampling volume 12) , a first amplifier 13 and a first peak detector 14 being strobed. The following elements are connected in series to the generator 8: a second transducer 15 comprising an emitter 16 and a receiver 17 (a gap between them forms the second sampling volume 18) , the second amplifier 19 and also the second peak detector 20 being strobed. Besides, former 21 of delayed strobe pulses and the first and the second peak detectors 14 and 20 are connected to generator 8. The latter peak detectors are connected to a calculator 24 and a display 25 respectively through a first and a second analogue-to-digital converter (ADC) 22 and 23. Sensors 9 and 15 are placed inside the pipeline in such a way that a flow at first pas. es through the one sampling volume, for example, volume 18, and then through the other, for example, volume 12. The sizes of transducers are chosen in such a way that they induce minimum disturbances of a flow (diameter of transducers < 3 mm) . Distance δ between emitter and receiver makes up about 2 mm and distance 1 between the lower and the upper pairs of transducers is equal to 3 ÷ 5 mm. Waveguides of the first and the second pairs of transducers in plan view of the sensor are located perpendicularly to one another that also improves hydrodynamics of a flow. An ultrasonic local velocity meter operates in the following way. Voltage pulses from the generator 8 are transmitted to the emitters 3 and 9, transformed into ultrasonic pulses and pass through the sampling volumes 12 and 18, then they are received by the receivers 11 and 17, transformed into a voltage signal, amplified by the amplifiers 13 and 19 and transmitted to the peak detectors 14 and 20 being strobed. Simultaneously with transmitting of the ultrasonic pulses the passage time of which is determined by a distance between the emitter and the receiver by the fixed pulse frequency the strobe pulses come to strobe inputs of the peak detectors 14 and 20. The strobe pulses switches the peak detectors into an active state. As a result voltage levels proportional to amplitudes of the received acoustic signals are formed at outputs of the peak detectors (see the voltage diagram shown in Figure 3) . After analogue-to-digital conversion in the ADC 22 and the ADC 23 the voltage signals are transmitted to the calculator 24 that calculates a cross-correlation function (CCF) for the received acoustic signals and presents it on the display 25. Because of discrete structure the multiphase mixture is an acoustically inhomogeneous medium. Hence an amplitude of received signals will fluctuate. Acoustic diffusers (the major part of them are gas inclusions, making the main contribution to diffusion of ultrasonic pulses) induce fluctuation at first when they pass through the second sampling volume. As a result an amplitude of an output signal at the second peak detector 20 changes and then with some delay equal to time of an acoustic diffuser pass from the second sampling volume to the first sampling volume τ, an amplitude of an output signal at the second 14 also changes. Statistical data accumulation of output signals of peak detectors provides a formation of a CCF maximum, its coordinate along a time axis is equal to τ. So a local gas velocity is determined by the expression: w q = 1 / τ , where 1 is the distance between the first and the second sampling volumes . A typical form of the cross-correlation function is shown in Figure 4. Another variant of acoustic sensing multiphase mixture by measuring of a local gas velocity is also possible. In this case two located in series acoustic transducers which operate in the mode "emission-reception" are used. The illustration of such a decision is presented in Figure 5. In this variant the velocity meter consists of two identical electro-acoustic channels, each one contains the following connected in series elements: an acoustic sensor 26, an amplifier 13, a peak detector 14 being strobed, an analogue-to-digital converter (ADC) 22 and also an electric pulse generator 8 connected to the sensor 26 through a sampling volume resistor 27 and a former of delayed strobe pulses 21. The former 21 is connected to a strobe input of the peak detector 14. Outputs of the channels are connected to a calculator 24 and then a display 25. Acoustic sensors are located inside a pipeline so that a flow 28 successively passes at first through a sampling volume of the first channel and then through a sampling volume of the second channel. The meter operates in the following way. Electric pulses from the generator 8 are transmitted to the acoustic sensor 26 where they are transformed into ultrasonic signals and emitted into the flow 28. Then a part of acoustic energy reflects from medium diffusers and comes back to the sensor 26, is amplified by the amplifier 13 and transmitted to the peak detector 14 being strobed. Simultaneously the delayed strobe pulse from the former 21 is transmitted to the strobe input of the peak detector 14 (see voltage diagram in Figure 6) . The resistor 27 executes uncoupling of an output of the generator 8 and an input of the amplifier 13. At the output of the peak detector 14 the voltage level proportional to an amplitude a of received signal is formed. Time of the strobe pulse delay relative to th pulse of the generator 8 t 0 (see Figure 6) is set taking into account passage time of the ultrasonic signal from the sensor to the sampling volume and back. A signal amplitude at the output of the peak detector fluctuates according to emergence of acoustic diffusers in a sampling volume. Since the diffusers at first pass through the sampling volume of the first sensor and then through the sampling volume of the second sensor the maximum on their CCF is formed. Co- ordinate τ of this maximum along the time axis is determined by passage time of the diffusers from the first sensor to the second one. A velocity of the diffusers contained in a medium is determined by the following formula: w g = I / τ , where 1 is the distance between the first and the second sensors. For calculation the CCF signals from the outputs of the peak detectors of the first and the second channels come through the ADC to the calculator 24. Results of calculations are shown on the display 25. Besides the above-described variant the local gas velocity meter can be implemented using one transducer with a pair of emitter and receiver of acoustic signals and also located inside the pipeline. The emitter and receiver are placed opposite to each other and form a sampling volume. The distance between them is chosen so that a mixture freely flows through the sampling volume. By an acoustic diffuser passage through a gap an ultrasonic signal is damped for a time equal to the time of the diffuser passage through the sampling volume. On the basis of these events an autocorrelation function of output signals is formed and the time of the diffuser passage through the sampling volume is determined. Illustration of this variant of the local gas velocity meter is shown in Figure 7. In this case the circuit contains elements connected in series: electric pulses generator 8, emitter 10 acoustically connected with receiver 11, amplifier 13, peak detector 14 being strobed, analogue- to-digital converter (ADC) 22, calculator 24 and display 25. The generator 8 is also connected through a former of delayed strobe pulses 21 with a strobe input of the peak detector. The space between the emitter (10) and receiver (11) represents a sampling volume 12. The velocity meter operates in the following way. Electric pulses from the generator 8 come to the emitter 10, are transformed into ultrasonic signals and through the sampling volume 12 come to the receiver 11, then to the amplifier 13 and to the peak detector 14. Simultaneously strobe pulses from the former 21 delayed for the time of the signal propagation from the emitter to the receiver are sent to the strobe input of the peak detector. Voltage from the input of the peak detector 14 proportional to the amplitude of a received signal is transferred to the ADC 22, then to the calculator 24 and the display 25. When diffusers of acoustic signal with sizes of particles less than the sampling volume are present in the flow every diffuser penetrating into the sampling volume will induce amplitude fluctuation of received signal. To the first approximation the time of amplitude fluctuation is equal to the time of diffuser passage through the sampling volume. Autocorrelation function determines the average time for statistical sampling of data. A typical form of the autocorrelation function is shown in Figure 8. So the local gas velocity can be calculated by the formula: where d is a linear size of a piezotransmitter plate in the flow direction, Xj is a main lobe width of the autocorrelation function (Figure 8) . One more variant of the local gas velocity meter is also possible. Its technical realization is shown in the Figure 9. In this variant a circuit of the velocity meter includes the following elements connected n series: an acoustic transducer 26, an amplifier 13, a peak detector 14 being strobed, an analog-to-digital transmitter 22, a calculator 24 and a display 25 and also a generator 8 connected through a resistor 27 with the transducer 26 and connected through a former of delayed strobe pulses 21 with a strobe input of the peak detector 14. The transducer 26 is situated inside a pipeline so that multiphase flow 28 crosses an acoustic field of the tr nsducer 26 perpendicularly to the flow direction. The meter operates in the following way. Voltage pulses from the generator 8 are transmitted through the resistor 27 to the transducer 26 where they are transformed into acoustic signals and emitted into the flow 28 perpendicularly to its direction. A portion of acoustic energy is reflected from acoustic diffusers of the multiphase medium (the main part of them are gas inclusions) and returns to the transducer 26 where it is transformed into electric signals which through the amplifier 13 come to the peak detector 14. Simultaneously a delayed strobe pulse from the former 21 is transmitted to the strobe input of the peak detector 14 (see the voltage diagram shown in Figure 10) . A resistor uncouples the output of the generator 8 and the input of the amplifier 13. The voltage amplitude at the output of the peak detector 14 is proportional to an amplitude of the received signal . The time of delay of a strobe pulse t 0 (see Figure 10) relative to a generator 8 pulse is set according to the time of ultrasonic signal passage from the transducer 26 to the sampling volume and back. A signal amplitude at the output of the peak detector fluctuates according to emergence of acoustic diffusers in the sampling volume. To the first approximation the time of fluctuation ' is equal to the time of diffuser pass through the sampling volume. Under condition that the sizes of diffusers are much smaller than the size of the sampling volume the local gas velocity can be determine by autocorrelation of the signals by the formula: W g = d / τ : , where d is a linear size of a piezotransmitter plate in the flow direction, τ λ is a main lobe width of the autocorrelation function (Figure 11) . Besides the above described variants another implementation of the ultrasonic local gas velocity meter using Doppler' s method of velocity determination is also possible. In this case the emitter and receiver with linear sizes up to 3 mm are set inside calibrated pipeline divisions under a fixed angle relative to each other. The circuit of the meter is shown in Figure 12. The meter contains an electric pulses generator 8 connected to emitter 10. Receiver 11 is connected through an amplifier 13 to a phase detector-multiplier 29. The following elements are connected in series to the output of the detector 29: a low-pass filter 30, the second amplifier 31, a signal spectrum calculator 32 and a display 25. A signal is processed in the measurement circuit in the following way. After the reflection of emitted ultrasonic oscillations from acoustic diffusers of a flow the acoustic signals come to the receiver 11, are transformed into voltage signals, are transmitted through the amplifier 13 to a first input of the phase detector 29. A voltage signal from the generator 8 is sent to a second input of the detector 29. From an output of the detector 29 low frequency signals are sent through a filter 30 and an amplifier 31 to the calculator 32 where a Doppler' s frequency proportional to a velocity of acoustic diffusers approach to a transducer is determined and then a local gas velocity is calculated. Results of processing are sent to the display 9. Signal processing in the circuit is shown in Figure 13. Another variant of technical implementation of the ultrasonic Doppler' s local gas velocity meter is demonstrated in Figure 14. Emitter and receiver with linear sizes up to 3 mm are also disposed inside calibrated pipeline divisions at a fixed angle relative to each other. Measuring circuit of the meter contains a voltage pulses generator 8 connected to an emitter 10. A receiver 11 is joined through an amplifier 13 to a phase detector-multiplier 29, its output is connected to a "sampling-storage" block 30. The second input of the phase detector 29 is linked to the generator 8. An input of the "sampling-storage" block 30 is connected to generator 8 through the former of delayed strobe pulses 21. An output of the block 30 is connected to the calculator 32 and then the display 25. The meter operates in the following way. Voltage pulses from the generator 8 are transmitted to the emitter 10 and induce acoustic pulses propagating opposite flow direction. Pulses reflected from acoustic diffusers mainly from gas bubbles come to the receiver 11 and through the amplifier 13 are transmitted to the first input of the phase detector-multiplier 29. A high-frequency signal from the generator 8 is sent to the second input of the detector 2 A A low-frequency signal from the detector 29 is sent to the "sampling-storage" block 30 which registers the signal at its input at the moments of time determined by the time position of the delayed strobe pulse from the former 21. Spectral processing of the signal from the "sampling-storage" block 30 is executed in the calculator 32 where the Doppler' s frequency proportional to the approach velocity of acoustic diffusers to the transmitter is isolated and a local gas velocity is calculated. Results of processing are shown on the display 25. Signal processing in the circuit is shown in Figure 15. Ultrasonic gas content meter (see Figure 16) contains a voltage pulses generator 8 connected in series to an emitter 10, which is acoustically linked to a receiver 11, an amplifier 33 and a peak detector 34 being strobed. The generator 8 is also connected to a strobe input of the peak detector 34 through a former of delayed strobe pulses 35. Output of the peak detector is linked to a direct input of a first comparator 36, to an inverse input of a second comparator 37 and to a calculator 24. Outputs of comparators 36 and 37 are also connected to the calculator 24 and then to a display 25. An inverse input of the first peak detector and a direct input of the second peak detector are connected to a first voltage setting device 38 and a second voltage setting device 39 respectively. The emitter and the receiver 10 and 11 are fixed one against the other forming a sampling volume 40. The meter operates in the following way. Rectangular voltage pulses produced by the generator 8 are transformed by the emitter 10 into ultrasonic pulses which are emitted into the sampling volume 40, reach the receiver 11, are transformed into voltage pulses and transmitted trough the amplifier 33 to the peak detector 34. A diagram of the signal processing in the elements of a meter measuring circuit is presented in the Figure 17. At the output of the peak detector 34 a level proportional to the amplitude of a signal being received at the moment of coming of a delayed strobe pulse is formed. The amplitude of the signal being received is determined by the gas volume concentration in the sampling volume 40. When the sampling volume is filled with liquid without gas inclusions the amplitude of the signal being received is at maximum and the voltage level at the input of peak detector 34 is higher than the voltage of the setting device 38 (Ul). That induces actuation of the comparator 36 and forming of an individual logic signal at its output. The logic signal is sent to the calculator 24 and is considered by the calculator 24 as a situation with gas volume concentration φ = 0 (see Figure 18) . Sizes of gas inclusions in a real multiphase flow are different and can be both less and greater than a size of the sampling volume 40. When the sizes of bubbles or gas plugs exceed the size of the sampling volume the propagation of the ultrasonic pulses is totally blocked, the amplitude of the signal being received diminishes to a minimum determined by a level of noises, and a voltage level at the output of the peak detector 34 is also minimum and is below a voltage of the setting device 39 (U2) . In this case the comparator 37 actuates and produces an individual logic signal considered by the calculator 24 as a situation with a gas volume concentration φ = 1. When the sizes of the gas bubbles are smaller than the size of the sampling volume 40 an output signal amplitude of the of peak detector 34 ranges from Ul to U2 (see Figure 18) and is described by the following relationship: U = U max exp (- k • n b • d b ) , (20) where U max is an amplitude of the signal when the liquid phase fills the controlled volume, k is a proportionality factor determined by geometric sizes of the sensor, by the ultrasonic frequency and so on, n b is concentration of gas bubbles, d b is diameter of gas bubbles. Taking into account that the concentration of bubbles changes continuously in the sampling volume owing to the mixture flowing the signal amplitude also fluctuates. Number of bubbles in the sampling volume is determined by Puasson' s law. So through the measurement of an average value of the received signal and of its dispersion values n b and d b are calculated using a known mathematical model by the calculator 24. The volume gas content is determined according to the formula : π d 3 1 φ 3 = N • • — , (21) 6 V where V is the sampling volume, N = n b • V is a number of the bubbles in the sampling volume. The gas phase concentration in case of variable composition of the gas inclusions in the flow is determined by the relationship: t ' 1+t * φ φ = , (22) T T = ti + t 2 + t 3 is time of averaging, where ti is a period of time when the gas inclusions are absent in the sampling volume, t 2 is a period of time when the gas inclusions presented by the bubbles of big diameter and also the gas plugs are present in the sampling volume, t 3 is a period of time when small bubbles are present in the sampling volume. The size of the sampling volume is chosen according to conditions of either technical implementation or application of the sensor, as a rule the size is smaller than 1 mm 3 . Block diagram of the ultrasonic meter of volume concentrations of liquid components is shown in the Figure 19. A measuring circuit of the meter comprises a generator of voltage pulses 8 and the following elements connected to it in series: an emitter 10 acoustically linked to a receiver 11, an amplifier 41, a first comparator 42, a first element 2& 43, a first RS-trigger 44, a second element 2& 45, a second RS-trigger 46 and a "duration- amplitude" converter 47. The generator 8 is also connected to a former of delayed strobe pulses 48 and to the second inputs of the RS-triggers 44 and 46. The second input of the first comparator 42 is linked to a voltage setting device 49. An output of the amplifier 41 is connected to a second comparator 50, its output is linked to the second input of the second element 2& 45. An output of the former of delayed strobe pulses 48 is connected to the second input of the first element 2& 43. The emitter and receiver are mounted in a tool body 51 one opposite the other thus forming a sampling volume 52. The tool body 51 is equipped by a heater 53 and an element 54 for mechanical cleaning of the sampling volume 52. The ultrasonic meter of volume concentration operates in the following way. Rectangular voltage pulses produced by the generator 8 are transformed into ultrasonic pulses by the emitter 10. After passing through the sampling volume 52 they reach the receiver 11 and are transformed into electric pulses. Then the signal through the amplifier 41 comes to the direct input of the first comparator 42. Simultaneously with sending the voltage pulse the first RS- trigger 44 is switched into a state "zero" and the second RS-trigger 46 - into a state "one". Since an inverse input of the comparator 42 is connected to the voltage setting device 49, actuation of the comparator 42 occurs when the amplitude of the received signal excesses a set voltage. Pulses from an output of the comparator 42 are transmitted to an S- input of the first RS-trigger 44 through the first element 2& 43 being strobed by pulses from the output of the former of delayed strobe pulses 48 and switch it into the state "one" (see the diagram of voltages in Figure 20) . A time lag is determined the by time of the ultrasonic pulses propagation from the emitter 10 to the receiver 11. Using a delay element excludes false actuations of the meter stipulated by electric and acoustic noises. Since one of the inputs of the second comparator 50 is connected to the ground wire it produces voltage pulses every time when the amplitude of the received signal crosses a "zero" mark so fixing even a weak signal (see Figure 20) . An output signal of the comparator doesn't depend on the amplitude of the received signal. A signal from an output of the first RS-trigger 44 transmitting to one of the inputs of the second element 2& 45 permits a passage through it of the signal from the second comparator 50 which indicate that the received signal crossed a "zero" mark. The first "zero" mark crossing induces actuation of the second RS-trigger 46 thus switching it into "zero" state. The voltage pulses formed so have durations proportional to passage time of the ultrasonic pulses from the emitter 10 to the receiver 11 and don't depend on ultrasonic pulses amplitudes. Then these pulses are transformed in the converter 47 into an amplitude signal proportional to their durations which is transmitted to the calculator and the monitor. In the second variant of the ultrasonic meter of volume concentrations of liquid components (see Figure 21) a voltage setting device is implemented as a peak detector 55 being strobed (see Figure 21) . Its input is connected to the output of the amplifier 41, a strobe input is linked to the output of the former of the delayed strobe pulses 48 and an output of the peak detector 55 is connected to the second input of the first comparator 42 through a voltage divider 56. The voltage setting device operates in the following way. A voltage signal from the amplifier 41 is transmitted to the peak detector 55. Simultaneously with a time lag determined by the passage time of the ultrasonic pulses from the emitter 10 to the receiver 11 the signal from the former of the delayed strobe pulses 48 comes to its strobe input (see Figure 22) . As a result a voltage potential equal to a maximum value of the signal amplitude is formed at the output of the peak detector 55. The voltage signal passed through the divider 56 is damped so that secure actuation of the first comparator 42 is guaranteed at the selected half-wave of the signal by its variations due to changes of the controlled medium properties and temperature and owing to aging of measuring circuit elements an so on. The use of such voltage setting device permits to support automatically the level of the comparator actuation by significant (10-fold) changes of signal attenuation in the medium stipulated, for example, by emergence of gas bubbles within the sampling volume, by change of dispersion of components and by other reasons. The operation of the meters of local values w g ,ι, w g , 2 , φi, Φ2 and W is controlled according to a set program by processor. Time averaging and cross-section averaging of the calibrated pipeline divisions of the above indicated values are also executed by means of the processor. The volume flow rates of the multiphase flow components such as liquid, oil, water and gas Q,, Qou, Q W Qg are moreover determined according to the formulas (16, 18, 19) by the processor . Though the invention described in particular for use with an oil, water and gas mixture, it must be understood that the principle of the invention, as set forth in the appended claims, is applicable for other mixtures also. Further, though the example of Figure 1 shows an order of pipeline divisions with decreasing cross-section areas seen in the flow direction an opposite order of pipeline line divisions, that is with increasing cross-section areas seen in the flow direction, can be used.

Claims: C L A I M S 1. Method for determination of characteristics of a flow of a multiphase liquid/gas mixture along a pipeline, comprising the steps of: a. changing a flow velocity of the mixture in a first pipeline division (1) and a second pipeline division (2) which is in series with the first pipeline division (1); b. measuring local flow characteristics of the flow in the first and the second pipeline divisions (1, 2); c. determining volume flow rates of gas phase (Q g ) , a first component of liquid phase (Q x ) and a second component of liquid phase (Q 2 ) of the mixture from the measured local flow properties, characterized in that d. for each pipeline division (1, 2) a real velocity (w) of at least one phase of the mixture is measured; e. in each pipeline division (1, 2) acoustic conductivity of the mixture is measured; f. a volume concentration of gas phase (φ) of the mixture is determined for each pipeline division (1,2) on the basis of the measured acoustic conductivity values; g. the volume concentration of liquid components (W) of the mixture is determined for pipeline division (1, 2) on the basis of measured acoustic conductivity values; h. volume flow rate values are determined on the basis of the velocity values and concentration values of the phases in each pipeline division (1, 2) . 2. Method according to claim 1 characterized in that the area (Fi) of the cross-section of the first pipeline division (1) differs from the area (F 2 ) of the cross-section of pipeline division (2)- F 2 = k Fi, where k φ 1. 3. Method according to claim 2, characterized in that the volume flow rate value of liquid phase is determined by the formula: Q,= k /(k - 1) Fi [w 2 (1 - φ 2 ) - ! (1 - φ ] , where i, w is an average real velocity of gas phase in the first pipeline division (1) and the second pipeline division (2) respectively, φi, φ 2 is an average real volume gas concentration in the mixture in the first pipeline division (1) and the second pipeline division (2) respectively; the volume flow rate of gas phase is determined by: Q g =F 1 w J φ 1 or Q g =F 2 w 2 φ 2 , the volume flow rate of the first component of liquid phase is determined by: Qι= W Q, , and the volume flow rate of the second component of liquid phase is determined by: Q 2 = (1 - W) Q, . 4. Method according to one of the preceding claims, characterized in that the velocity (w) of gas phase is measured at different radial locations in each of said cross-sections of the first and second pipeline divisions (1, 2), and measured local velocity values for each cross-section are averaged to provide a value for use as velocity value in calculations. 5. Method according to one of the preceding claims, characterized in that the concentration of gas phase (φ) is measured at different radial locations in each of said cross-sections of the first and the second pipeline divisions (1, 2), and measured volume concentration values for each cross-section are averaged to provide a value for use as concentration value in calculations. 6. Method according to one of the preceding claims, characterized in that measurements are carried out by the use of ultrasonic transducers . 7. Method according to a preceding claim, characterized in that volume concentrations of liquid phase components (W) of the mixture are determined by using ultrasonic transducers in at least one cross-section in at least one of the pipeline divisions (1, 2) and by measuring time of passage of ultrasonic pulses through the mixture from the transducers. 8. Method according to a preceding claim, characterized in that volume concentrations of gas phase components (φ) of the mixture are determined by using ultrasonic transducers in at least one cross- section of the first and the second pipeline divisions (1, 2) and by measuring amplitude of ultrasonic pulses passed through the mixture from the transducers. 9. Method according to a preceding claim, characterized in that velocities (w) of the mixture phases are determined by using ultrasonic transducers in at least one cross-section of the first and the second pipeline divisions (1, 2) and by cross-correlation or autocorrelation methods. 10. Method according to a preceding claim, characterized in that velocities (w) of the mixture phases are determined by using ultrasonic transducers in at least one cross-section of the first and the second pipeline divisions (1, 2) and by measuring Doppler frequency of ultrasonic pulses from the transducers. 11. Method according to claims 1, 2, 3, 4, 5, characterized in that measurements are carried out by the use of electrical conductivity transducers. 12. Method according to claims 1, 2, 3, 4, 5, characterized in that measurements are carried out by the use of electrical capacity transducers . 13. Method according to claims 1, 2, 3, 4, 5, characterized in that measurements are carried out by the use of optical sensors. 14. Method according to claim 8, characterized in that the liquid phase components of the mixture are water and oil. 15. Device for determination characteristics of a flow of a multiphase liquid/gas mixture along a pipeline, comprising: i. a flow cell for changing a flow velocity of the mixture in a first pipeline division (1) and a second pipeline division (2) which is in series with the first division (1) while the both divisions are disposed in the flow cell; j . means for measuring local flow characteristics of the flow in the first and the second pipeline divisions (1, 2) ; k. means for determining volume flow rates of gas phase (Q g ) , a first component of liquid phase (Q and a second component of liquid phase (Q 2 ) of the mixture from the measured local flow properties. 16. Device according to claim 15, in which for each pipeline division (1, 2) the following means for measuring local flow characteristics of at least one phase of the mixture are used: 1. an ultrasonic gas velocity meter for measuring a real gas velocity (w) of the mixture based on correlation or Doppler' s methods; m. an ultrasonic volume gas concentration meter; n. an ultrasonic meter of volume concentrations of liquid components. 17. Device according to claim 15, in which meters with electrical capacity or electrical conductivity transducers for measuring local flow properties of at least one phase of the mixture are used. 18. Device according to claim 15, in which for each pipeline division (1, 2 ) a gamma-meter for determining the volume gas concentration is used. 19. Device according to claim 15, in which for each pipeline division (1, 2) disposed vertically an apparatus for determining the volume gas concentration by means of measuring a static pressure difference is used. 20. Device according to claims 15, 16, 17, in which the gas velocity meter measures the velocity (w) of gas phase at different radial locations in each of said cross-sections of the first and the second pipeline divisions (1, 2) , and measured local velocity values for each cross-section are averaged to provide a value for use as velocity value in calculations. 21. Device according to claims 15, 16, 17, in which the volume gas concentration meter measures the concentration of gas phase (φ) at different radial locations in each of said cross-sections of the first and the second pipeline divisions (1, 2) , and measured volume concentration values for each cross-section are averaged to provide a value for use as concentration value in calculations. 22. Device according to claims 15, 16, 17, 18, 19, 20, 21, in which the area (F x ) of the cross-section of the first pipeline division (1) differs from the area (F 2 ) of the cross-section of pipeline division (2) F 2 = k Fi, where k Φ 1. 23. Device according to claims 15, 16, 17, 18, 19, 20, 21, including a processor calculating the volume flow rate value of liquid phase by the formula: Q ( = k /(k - 1) Fx [w 2 (1 - φ 2 ) - wx (1 - Φl ) ] , where w L , w 2 is an average real velocity of gas phase in the first pipeline division (1) and the second pipeline division (2) respectively, ψι , φ 2 is an average real volume gas concentration in the mixture in the first pipeline division (1) and the second pipeline division (2) respectively; calculating the volume flow rate of gas phase by the formulas: Q g =FiW^ x or Q g =F 2 w 2 φ 2 , calculating the volume flow rate of the first component of liquid phase by the formula: Qι= W Q ( , and calculating the volume flow rate of the second component of liquid phase by the formula: Q 2 = (1 - W) Q, . 24. Device according to claims 15, 16, 17, 18, 19, 20, 22, 23, in which the sensed liquid phase components of the multiphase flow are water and oil.

Inventors: DROBKOV VLADIMIR, MELNIKOV VLADIMIR, SHUSTOV ANDREY, DROBKOV, VLADIMIR, MELNIKOV, VLADIMIR, SHUSTOV, ANDREY
Applicants: DROBKOV VLADIMIR, MELNIKOV VLADIMIR, NEST INTERNAT N V, SHUSTOV ANDREY, DROBKOV, VLADIMIR, MELNIKOV, VLADIMIR, NEST INTERNATIONAL N.V., SHUSTOV, ANDREY
Application year: 2001
Main classification: G01F1
IPC codes: G01F1, G01N29