DEVICE AND METHOD FOR TOTAL CROSS-SECTION MEASUREMENT OF MASS FLOW RATE OF GAS, LIQUID AND SOLID IN MULTIPHASE FLOW

Abstract

A device for total cross-section measurement of a mass flow rate of gas, liquid and solid in a multiphase flow includes a gamma-ray source, a gamma-ray detector, and a differential pressure type flowmeter. The differential pressure type flowmeter includes a throat section, and the gamma-ray source and the gamma-ray detector are respectively disposed at opposite positions on both sides of the throat section. The gamma-ray detector is an array including a plurality of detection units, and the gamma-ray source is configured to emit gamma rays covering the measurement cross-section of the throat section. The gamma-ray detector is configured to receive the gamma rays passing through the measurement cross-section of the throat section.

Claims

1. A device, comprising a gamma-ray source, a gamma-ray detector, and a differential pressure type flowmeter; wherein the differential pressure type flowmeter comprises a throat section, and the gamma-ray source and the gamma-ray detector are respectively disposed at opposite positions on both sides of the throat section; the gamma-ray detector is in the form of an array comprising a plurality of detection units, and the gamma-ray source is configured to emit gamma rays covering a measurement cross-section of the throat section, and the gamma-ray detector is configured to receive the gamma rays passing through the measurement cross-section of the throat section.

2. The device of claim 1, wherein the gamma-ray source is configured to emit gamma rays of at least three energy levels.

3. The device of claim 1, wherein the gamma-ray source and the gamma-ray detector cooperate with each other to form a detection pair, and at least one detection pair is disposed along an axial direction of the differential pressure type flowmeter.

4. The device of claim 3, wherein when a plurality of detection pairs is disposed along the axial direction of the differential pressure type flowmeter, a distance between every two adjacent detection pairs is 1 to 10 cm.

5. The device of claim 4, wherein the plurality of detection pairs is staggered with each other along the axial direction of the differential pressure type flowmeter and distributed around the axial direction.

6. The device of claim 1, wherein a scintillation crystal is connected to one end of each detection unit receiving the gamma rays, and the detection unit is a silicon photomultiplier (SiPM).

7. The device of claim 1, wherein the scintillation crystal and the detection unit are fixed to each other by a couplant.

8. The device of claim 1, further comprising a temperature and pressure sensor for measuring a temperature and pressure of a multiphase fluid, and a differential pressure sensor for measuring a pressure difference between an inlet of the differential pressure type flowmeter and the throat section.

9. A method for total cross-section measurement of a mass flow rate of gas, liquid and solid in a multiphase flow using the device of claim 1, comprising: a) measuring a temperature T of a multiphase liquid by a temperature sensor, measuring a pressure difference ΔP between an inlet of the differential pressure type flowmeter and the throat section; and measuring transmission intensities N.sub.x1, N.sub.x2, N.sub.x3 of three gamma rays by the gamma-ray detector; b) calculating a total mass flow rate of the multiphase fluid and respective mass flow rates of gas, liquid and solid phases of the multiphase fluid according to the following formula: total mass flow rate: $Q_m=\frac{C}{\sqrt{1-{\beta }^4}}\mathrm{.Math.}\frac{\pi }{4}{D}^2\sqrt{\Delta P{\rho }_{mix}};$
Q.sub.m,g=Q.sub.m*GMF;  gas mass flow rate:
Q.sub.m,l=Q.sub.m*LMF;  liquid mass flow rate:
Q.sub.m,s=Q.sub.m*SMF;  solid mass flow rate: where, gas mass fraction, $GMF=\frac{Q_g}{Q_g+Q_l+Q_s},$ liquid mass fraction, $LMF=\frac{Q_l}{Q_g+Q_l+Q_s},$ solid mass fraction, $SMF=\frac{Q_s}{Q_g+Q_l+Q_s},$ where, Q.sub.g, Q.sub.l, Q.sub.s are respectively linear masses of the gas, liquid, and solid phases, and $Q_g=\theta -\frac{\left(k_2-1\right)\left(d_1-d_2\right)+\left(k_4-k_2\right)\left(d_1-\theta \right)}{\left(k_2-1\right)\left(k_1-k_3\right)+\left(k_4-k_2\right)\left(k_1-1\right)}-\frac{\left(k_3-1\right)\left(d_1-d_2\right)+\left(k_3-k_1\right)\left(d_1-\theta \right)}{\left(k_3-1\right)\left(k_2-k_4\right)+\left(k_3-k_1\right)\left(k_4-1\right)};$ $Q_l=\frac{\left(k_2-1\right)\left(d_1-d_2\right)+\left(k_4-k_2\right)\left(d_1-\theta \right)}{\left(k_2-1\right)\left(k_1-k_3\right)+\left(k_4-k_2\right)\left(k_1-1\right)};$ $Q_s=\frac{\left(k_3-1\right)\left(d_1-d_2\right)+\left(k_3-k_1\right)\left(d_1-\theta \right)}{\left(k_3-1\right)\left(k_2-k_4\right)+\left(k_3-k_1\right)\left(k_4-1\right)};$ $\mathrm{where},\theta =\frac{\pi }{4}{\left(\frac{c}{\sqrt{1-{\beta }^4}}\mathrm{.Math.}D\right)}^2\Delta P\frac{t}{l};$ $k_1=\frac{a_{l,1}-a_{l,2}}{a_{g,1}-a_{g,2}};$ $k_2=\frac{a_{s,1}-a_{s,2}}{a_{g,1}-a_{g,2}};$ $d_1=\frac{1}{a_{g,1}-a_{g,2}}\ln \left(\frac{f_1}{N_{x,2}}\right)-\ln \left(N_{x,1}\right);$ $k_3=\frac{a_{l,1}-a_{l,3}}{a_{g,1}-a_{g,3}};$ $k_4=\frac{a_{s,1}-a_{s,3}}{a_{g,1}-a_{g,3}};$ $d_2=\frac{1}{a_{g,1}-a_{g,3}}\ln \left(\frac{f_2}{N_{x,3}}\right)-\ln \left(N_{x,1}\right);$ letters in each formula are defined as follows: C: a discharge coefficient of a differential pressure type flowmeter; ε: a compression correction factor of the multiphase fluid; β: a diameter ratio of the differential pressure type flowmeter; D: a thickness measured by gamma rays, also a diameter of a fluid pipe; ΔP: the pressure difference; f1, f2: initial intensity ratios of a second gamma ray and a third gamma ray to a first gamma ray; ρ.sub.mix: an average density of the measured fluid; ${\rho }_{mix}=\left(Q_g+Q_l+Q_s\right)/S\frac{t}{l};$ S is an area of the measured cross-section; $S=\frac{\pi }{4}{D}^2,$ l is a unit length, t is a unit time; a is a mass absorption coefficient of the multiphase fluid for gamma rays, Q is the mass flow rate of the multiphase fluid to be measured, subscripts 1, 2 and 3 respectively represent gamma rays of different energy levels; all calculation results are weighted and averaged to obtain a final and accurate linear mass value of each phase according to the following formulas:
Q.sub.g=Σ(Xi/D*Q.sub.gi)/Σ(Xi/D)
Q.sub.l=Σ(Xi/D*Q.sub.li)/Σ(Xi/D)
Q.sub.s=Σ(Xi/D*Q.sub.si)/Σ(Xi/D) Q.sub.g, Q.sub.l, Q.sub.s are final linear masses of gas, liquid and solid of the multiphase flow; D is a diameter of the throat section of the differential pressure type flowmeter, X is a distance passed by each gamma ray in the throat section, Q.sub.gi, Q.sub.li, Q.sub.si are linear mass values measured after each gamma ray passes through the throat section.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0053] FIG. 1 is a schematic view of a device for total cross-section measurement of a mass flow rate of gas, liquid and solid in a multiphase flow of Example 1;

[0054] FIG. 2 is a schematic diagram of a path of the gamma rays in Example 1;

[0055] FIG. 3 is a schematic view of a gamma-ray detector in Example 1;

[0056] FIG. 4 is a schematic diagram of a path of the gamma rays in Example 2;

[0057] FIG. 5 is a schematic view of a device for total cross-section measurement of a mass flow rate of gas, liquid and solid in a multiphase flow of Example 3;

[0058] FIG. 6 is a three-dimensional diagram of a device for total cross-section measurement of a mass flow rate of gas, liquid and solid in a multiphase flow of Example 4;

[0059] FIG. 7 is an exploded view of a device for total cross-section measurement of a mass flow rate of gas, liquid and solid in a multiphase flow of Example 4; and

[0060] FIG. 8 is a schematic diagram of a path of the gamma rays in Example 4.

DETAILED DESCRIPTION

[0061] To further illustrate, embodiments detailing a device and a method for total cross-section measurement of a mass flow rate of gas, liquid and solid in a multiphase flow are described below. It should be noted that the following embodiments are intended to describe and not to limit the disclosure.

Example 1

[0062] As shown in FIG. 1, a device for total cross-section measurement of mass flow rate of gas, liquid and solid in a multiphase flow comprises a gamma-ray source 1, a gamma-ray detector 2, and a differential pressure type flowmeter 3. The differential pressure type flowmeter 3 comprises a throat section 31, and the gamma-ray source 1 and the gamma-ray detector 2 are respectively disposed at opposite positions on both sides of the throat section 31; the device further comprises a temperature and pressure sensor for measuring the temperature and pressure of the multiphase fluid and a differential pressure sensor for measuring the pressure difference between the inlet of the differential pressure type flowmeter 3 and the throat section 31.

[0063] The basic principle of the differential pressure type flowmeter 3 is as follows. In a round pipe filled with fluid, throttling devices such as a venturi, an orifice plate or a nozzle are disposed, and the place with the smallest diameter is called a throat. When the multiphase fluid flows through the throttling device, a static pressure difference will be generated between its upstream and the throat. There is a fixed functional relationship between the static pressure difference and the flow rate. As long as the static pressure difference is measured, the flow rate can be obtained from the flow formula.

[0064] As shown in FIG. 3, the gamma-ray detector 2 is in the form of an array comprising a plurality of semiconductor detectors 21. In this example, the semiconductor detector 21 is SIPM. The gamma rays emitted by the gamma-ray source 1 cover the total cross-section where the throat section 31 is disposed, and all gamma rays passing through the cross-section where the throat section 31 is disposed can be received by the gamma-ray detector 2.

[0065] As shown in FIG. 2, the gamma-ray source 1 and gamma-ray detector 2 are arranged during mounting, such that the rays emitted from one gamma-ray source 1 can fully cover the cross-section of the throat section 31 of the differential pressure type flowmeter 3, and the gamma rays can fully pass through the multiphase fluid flowing through the cross-section and are received by the gamma-ray detector 2 on the other side of the differential pressure type flowmeter 3, and the total cross-section of the multiphase fluid is detected spatially. The spatial inhomogeneity of the fluid is solved properly to make the measured data more accurate than those measured when only part of the gamma rays passes through the cross-section. In this example, only one radioactive source emits gamma rays, which will not generate interference and make the measured data more accurate.

[0066] To receive all the gamma rays passing through the multiphase fluid cross-section, the range of the gamma-ray detector 2 must be large. Therefore, in this example, an array comprising a plurality of semiconductor detectors 21 is adopted. As shown in FIG. 3, the semiconductor detector 21 forms a 4*4 square matrix, and nine 4*4 SIPM square matrices form the entire gamma-ray detector 2. The gamma-ray detector 2 is set to have a receiving range as large as possible, and can receive each ray without omission, making measurement more accurate; the semiconductor detector 21 is arranged in a square matrix to make it easier to calculate the phase fraction.

[0067] As shown in FIG. 2, in this example, a scintillation crystal 22 is connected to one end of each semiconductor detector 21 that receives gamma rays.

[0068] Scintillation crystals 22, as a material often used in ray detection technology, are capable of converting high-energy gamma rays into low-energy visible light fluorescence, which are then detected by a silicon photomultiplier 21 to convert into electrical signals. At present, commonly used scintillation crystals 22 include sodium (thallium) iodide NaI (Tl), lutetium oxyorthosilicate (LSO), etc. A silicon photomultiplier (SiPM) is a novel detector. After photons are absorbed, a current is generated in the SiPM and multiplied, which can output a larger current signal and can be received by the modular circuit. Its detection efficiency of gamma rays is higher and its volume is smaller.

[0069] The scintillation crystal 22 and the semiconductor detector 21 are fixed to each other by a couplant.

[0070] The scintillation crystal 22 is a high-density crystal. The surface of the silicon photomultiplier 21 is provided with a layer of epoxy resin. When the light is emitted from the scintillation crystal 22 to the silicon photomultiplier 21, it is emitted from an optically denser medium to an optically thinner medium. If air exists between them, total reflection easily occurs, causing light loss. Optical couplants, especially optical couplers, are transparent media with large refractive index. When couplants are disposed between the scintillation crystal 22 and the silicon photomultiplier 21, the air can be effectively eliminated and the light loss caused by total reflection can be significantly reduced. The couplant may be a silica gel to bond the scintillation crystal 22 to the silicon photomultiplier 21, effectively reducing the loss of light from the scintillation crystal 22 to the silicon photomultiplier 21 and improving the photoelectric conversion efficiency.

[0071] In this example, the gamma-ray source 1 is configured to emit gamma rays of at least three energy levels.

[0072] A multi-energy radiation source that can naturally emit gamma rays with more than three types of energy is used. Because the ratio of intensity of the naturally emitted gamma rays with three types of energy is inherent and constant that cannot be altered by humans and is not affected by any change in the external temperature and pressure, it brings great convenience and simplification to the solution of the flow formula of the disclosure. For example, using .sup.133Ba, the gamma rays emitted by the radiation source have three main energy levels, 31 keV, 81 keV, and 356 keV respectively; or using .sup.176Lu, the gamma rays emitted by the radiation source have at least three energy levels, namely, 307 keV, 202 keV, and 88 keV.

[0073] The working principle of the device is as follows. The multiphase fluid flows through the differential pressure type flowmeter 3, and the gamma rays emitted by the gamma-ray source 1 pass through the total cross-section of the multiphase fluid and are received by the gamma-ray detector 2 for photoelectric conversion. After the conversion, the data imaging, analysis and calculation are performed, to obtain the mass phase fraction of the multiphase fluid.

Example 2

[0074] As shown in FIG. 4, the difference between this example and Example 1 is that the device in Example 1 comprises only one gamma-ray source 1 and one gamma-ray detector 2. The one gamma-ray source 1 and the one gamma-ray detector 2 are defined as a detection pair. In this example, two detection pairs are disposed around the throat section 31 of the differential pressure type flowmeter 3, and the two detection pairs are arranged along the axis of the differential pressure type flowmeter 3.

[0075] Each detection pair can perform total cross-section measurement of the phase fractions of the multiphase fluid, and the data measured by each detection pair can be averaged, to make the measurement result to be more accurate.

[0076] In this example, the distance between any two adjacent detection pairs is 10 cm.

[0077] Theoretically, the state of the multiphase fluid changes little in a short distance and a short period of time. Therefore, the closer the distance between the two detection pairs, the smaller the state change of the multiphase fluid, and the more accurate the measured data. Since the velocity of the multiphase fluid is generally 10 m/s and the distance is set to 10 cm, the state of the multiphase fluid is almost unchanged within this distance, so that the data measured by each detection pair are considered to be from the same cross-section.

[0078] In this example, the two detection pairs can be staggered by 90° from each other, that is, the radiation angle of the gamma-ray source 1 is staggered by 90°. Thus, the multiphase fluid cross-section is measured from different angles, to make the data more accurate.

Example 3

[0079] As shown in FIG. 5, the difference between this example and Example 1 is that the device in Example 1 comprises only one gamma-ray source 1 and one gamma-ray detector 2. The one gamma-ray source 1 and the one gamma-ray detector 2 are defined as a detection pair. In this example, three detection pairs are disposed around the throat section 31 of the differential pressure type flowmeter 3, and the three detection pairs are arranged along the axis of the differential pressure type flowmeter 3.

[0080] Each detection pair can perform total cross-section measurement of the phase fractions of the multiphase fluid, and the data measured by each detection pair can be averaged, to make the measurement result to be more accurate.

[0081] In this example, the distance between any two adjacent detection pairs is 10 cm.

[0082] Theoretically, the state of the multiphase fluid changes little in a short distance and a short period of time. Therefore, the closer the distance between the two detection pairs, the smaller the state change of the multiphase fluid, and the more accurate the measured data. Since the velocity of the multiphase fluid is generally 10 m/s and the distance is set to 10 cm, the state of the multiphase fluid is almost unchanged within this distance, so that the data measured by each detection pair are considered to be from the same cross-section.

Example 4

[0083] As shown in FIG. 6 to FIG. 8, the difference between this example and Example 3 is that the detection pairs are staggered along the axis of the differential pressure type flowmeter 3.

[0084] Since both the gamma-ray source 1 and the gamma-ray detector 2 have a certain volume, if all the detection pairs are disposed in the same position and the distance therebetween is too small, they will easily touch each other, and if the distance therebetween is too large, the measurement accuracy will decrease. Therefore, the detection pairs are staggered to make full use of the space for the layout of the gamma-ray source 1 and the gamma-ray detector 2, which can ensure a short distance between the two detection pairs without interfering with each other and occupying space; and after staggering, adjacent gamma-ray sources 1 emit gamma rays in different directions, which will not easily interfere with each other and will improve measurement accuracy.

[0085] In addition, in this example, the detection pairs that are staggered from each other can be evenly distributed and disposed around the axis, which is more reasonable in layout and more accurate in detection data.

Example 5

[0086] This embodiment is a method for total cross-section measurement of mass flow rate of gas, liquid and solid in a multiphase flow in the above embodiment.

[0087] To facilitate the understanding of this disclosure, some terms in the field of multiphase fluid measurement are described as follows.

[0088] “Mass flow rate” refers to the mass of fluid flowing in a unit of time. In the SI unit system, its dimension can be kg/s.

[0089] “Volume flow rate” refers to the volume of fluid flowing in a unit of time. In the SI unit system, its dimension can be m.sup.3/s.

[0090] According to the nature of the penetrated fluid, there are three mass flow rates Q.sub.g, Q.sub.l, Q.sub.s, respectively, namely gas mass flow rate, liquid mass flow rate and solid mass flow rate. Based on the mass flow rates of gas, liquid, and solid, the relationship between the total mass flow rate and the diameter of the pipe is as follows:

[00008]$\begin{array}{cc}{Q}_g+Q_l+Q_s=\frac{\pi }{4}{\left(\frac{c}{\sqrt{1-{\beta }^4}}\mathrm{.Math.}D\right)}^2\Delta P\frac{t}{l}.& \left(8\right)\end{array}$

[0091] “Radial direction” refers to the direction along the diameter of the cross-sectional circle of the fluid pipe.

[0092] In the following text, the measurement method of multiphase fluid mass flow rate of this embodiment is described in detail.

[0093] In this example, a conventional differential pressure type flowmeter 3, such as a Venturi flowmeter is used to measure the pressure difference, and then the total mass flow rate of the multiphase fluid is calculated by the following formula:

[00009]$\begin{array}{cc}{Q}_m=\frac{c}{\sqrt{1-{\beta }^4}}\mathrm{.Math.}\frac{\pi }{4}{D}^2\sqrt{\Delta P{\rho }_{mix}};& \left(9\right)\end{array}$

[0094] Where, C is the discharge coefficient of the differential pressure type flowmeter, ε is the multiphase fluid compression correction factor, β is the diameter ratio of differential pressure type flowmeter, ΔP is the pressure difference, ρ.sub.mix is the multiphase fluid density (mixed density in terms of the multiphase fluid), and Dis the diameter of the fluid pipe.

[0095] Next, by using the gamma-ray detector 2 of the multi-energy radiation source, the mass flow rate of each of the three phases (gas, liquid and solid) in the multiphase fluid is measured.

[0096] Firstly, based on the gamma ray absorption equations:

[0097] The absorption equation of gamma ray 1:

[00010]$\begin{array}{cc}\ln \left(\frac{N_{0,1}}{N_{x,1}}\right)=a_{g,1}{Q}_g+a_{l,1}{Q}_l+a_{s,1}{Q}_s;& \left(10\right)\end{array}$

[0098] The absorption equation of gamma ray 2:

[00011]$\begin{array}{cc}\ln \left(\frac{N_{0,2}}{N_{x,2}}\right)=a_{g,2}{Q}_g+a_{l,2}{Q}_l+a_{s,2}{Q}_s;& \left(11\right)\end{array}$

[0099] The absorption equation of gamma ray 3:

[00012]$\begin{array}{cc}\ln \left(\frac{N_{0,3}}{N_{x,3}}\right)=a_{g,3}{Q}_g+a_{l,3}{Q}_l+a_{s,3}{Q}_s;& \left(12\right)\end{array}$

[0100] Secondly, based on the relationship between the mass flow rate measured by Venturi and the linear mass, and the following equation:

[00013]$\begin{array}{cc}{Q}_g+Q_l+Q_s=\frac{\pi }{4}{\left(\frac{c}{\sqrt{1-{\beta }^4}}\mathrm{.Math.}D\right)}^2\Delta P\frac{t}{l};& \left(13\right)\end{array}$

[0101] Where, Q.sub.g, Q.sub.l, Q.sub.s are the respective mass flow rates of the gas, liquid, and solid phases.

[0102] According to the characteristics of the radioactive source, there is a proportional relationship among N.sub.o,1, N.sub.o,2 and N.sub.o,3:

[0103] N.sub.0,2=f.sub.1N.sub.0,1, N.sub.0,3=f.sub.2N.sub.0,1, where, f.sub.1 and f.sub.2 are the known proportional coefficients, which are naturally constant coefficients and do not change with any measurement conditions. Due to the existence of the proportional coefficients, the three unknown quantities N.sub.0,2, N.sub.0,3,N.sub.0,1 can actually be regarded as only one unknown quantity N.sub.0,1.

[0104] In this way, the four unknown quantities (N.sub.0,1, Q.sub.g, Q.sub.l, Q.sub.s) can be solved directly and accurately through the four equations (10) to (13) above, thereby eliminating the need for measurement or calibration of N.sub.0,1. Since it is not required to calibrate N.sub.0,1 (i.e. the count value of empty pipe), the effect of temperature drift in the gamma ray receiver on the measurement is fundamentally avoided, and there is no need to dispose a constant temperature device in the gamma ray receiver.

[0105] In the equation set, a.sub.g,1, a.sub.g,2, a.sub.g,3, a.sub.l,1, a.sub.l,2, a.sub.l,3 and a.sub.s,1, a.sub.s,2, a.sub.s,3 are the linear mass absorption coefficients of gas, liquid, and solid to the gamma ray 1, gamma ray 2 and gamma ray 3 under working conditions, and, f.sub.1, f.sub.2 are fixed values, which can be obtained by calibration. N.sub.x,1, N.sub.x,2, N.sub.x,3, ΔP are the measured values, so the mass flow rates (Q.sub.g, Q.sub.l, Q.sub.s) can be solved directly according to the following equations:

[00014]$\begin{array}{cc}{Q}_g=\theta -\frac{\left(k_2-1\right)\left(d_1-d_2\right)+\left(k_4-k_2\right)\left(d_1-\theta \right)}{\left(k_2-1\right)\left(k_1-k_3\right)+\left(k_4-k_2\right)\left(k_1-1\right)}-\frac{\left(k_3-1\right)\left(d_1-d_2\right)+\left(k_3-k_1\right)\left(d_1-\theta \right)}{\left(k_3-1\right)\left(k_2-k_4\right)+\left(k_3-k_1\right)\left(k_4-1\right)};& \left(14\right)\\ Q_l=\frac{\left(k_2-1\right)\left(d_1-d_2\right)+\left(k_4-k_2\right)\left(d_1-\theta \right)}{\left(k_2-1\right)\left(k_1-k_3\right)+\left(k_4-k_2\right)\left(k_1-1\right)};& \left(15\right)\\ Q_s=\frac{\left(k_3-1\right)\left(d_1-d_2\right)+\left(k_3-k_1\right)\left(d_1-\theta \right)}{\left(k_3-1\right)\left(k_2-k_4\right)+\left(k_3-k_1\right)\left(k_4-1\right)};& \left(16\right)\end{array}$

[0106] Then, according to the equation of the mass flow rate calculated by Venturi:

[00015]$Q_t=\frac{C}{\sqrt{1-{\beta }^4}}\mathrm{.Math.}\frac{\pi }{4}{D}^2\sqrt{\Delta P{\rho }_{mix}}$

and the definition of mass phase fraction, the mass flow rates of gas, liquid, and solid phases and the total mass flow rate are calculated according to the following equations:

Q.sub.m,g=Q.sub.m*GMF  (17);

Q.sub.m,l=Q.sub.m*LMF  (18);

Q.sub.m,s=Q.sub.m*SMF  (19);

[0107] In the above equations,

[00016]$\begin{array}{c}\theta =\frac{\pi }{4}{\left(\frac{c}{\sqrt{1-{\beta }^4}}\mathrm{.Math.}D\right)}^2\Delta P\frac{t}{l};\\ k_1=\frac{a_{l,1}-a_{l,2}}{a_{g,1}-a_{g,2}};\\ k_2=\frac{a_{s,1}-a_{s,2}}{a_{g,1}-a_{g,2}};\\ d_1=\frac{1}{a_{g,1}-a_{g,2}}\ln \left(\frac{f_1}{N_{x,2}}\right)-\ln \left(N_{x,1}\right);\\ k_3=\frac{a_{l,1}-a_{l,3}}{a_{g,1}-a_{g,3}};\\ k_4=\frac{a_{s,1}-a_{s,3}}{a_{g,1}-a_{g,3}}\\ d_2=\frac{1}{a_{g,1}-a_{g,3}}\ln \left(\frac{f_2}{N_{x,3}}\right)-\ln \left(N_{x,1}\right);\end{array}$

[0108] C: The discharge coefficient of the differential pressure type flowmeter;

[0109] ε: The multiphase fluid compression correction factor;

[0110] β: The diameter ratio of differential pressure type flowmeter;

[0111] D: The thickness measured by gamma rays, namely, the diameter of the fluid pipe;

[0112] ΔP: The pressure difference;

[0113] ρ.sub.mix: The average density of the measured fluid;

[00017]${\rho }_{mix}=\left(Q_g+Q_l+Q_s\right)/S\frac{t}{l};$

[0114] S is the area of the measured cross section,

[00018]$S=\frac{\pi }{4}{D}^2,$

is the unit length, t is the unit time;

[0115] Gas mass fraction,

[00019]$GMF=\frac{Q_g}{Q_g+Q_l+Q_s};$

[0116] Liquid mass fraction,

[00020]$LMF=\frac{Q_l}{Q_g+Q_l+Q_s};$

[0117] Solid mass fraction,

[00021]$SMF=\frac{Q_s}{Q_g+Q_l+Q_s};$

[0118] Q.sub.g, Q.sub.l, Q.sub.s are the mass flow rates of gas, liquid and solid to be solved respectively;

[0119] a is the mass absorption coefficient of the multiphase fluid for gamma rays, Qis the mass flow rate of the multiphase fluid to be measured, the subscripts 1, 2 and 3 respectively represent gamma rays of different energy levels;

[0120] Since the path of each gamma ray passing through the multiphase fluid is different, the calculation of the phase fractions will also be different. Finally, all calculation results are weighted and averaged to obtain the final and accurate linear mass value of each phase according to the following formulas:

Q.sub.g=Σ(Xi/D*Q.sub.gi)/Σ(Xi/D)

Q.sub.l=Σ(Xi/D*Q.sub.li)/Σ(Xi/D)

Q.sub.s=Σ(Xi/D*Q.sub.si)/Σ(Xi/D)

[0121] Q.sub.g, Q.sub.l, Q.sub.s are the final linear masses of gas, liquid and solid to be solved. D is the diameter of the throat section (31) of the differential pressure type flowmeter (3), X is the distance passed by each gamma ray in the throat section (31), Q.sub.gi, Q.sub.li, Q.sub.si are the linear mass values measured after each gamma ray passes through the throat section (31); conduct statistics on all gamma rays and their linear mass data of each phase according to the formula, to calculate the final accurate values.

[0122] The measurement device and measurement method described in this example are described for measuring and calculating the mass flow rates of three phases (gas, liquid, and solid) in a multiphase fluid. The device and method are also applicable to the measurement of two-phase flow and the calculation of the respective mass flow rates of the gas phase and the liquid phase. Accordingly, the principle and method of calculating the mass flow rates according to the two energy levels of the gamma-ray radiation sources can be analogized based on the above content.

[0123] It will be obvious to those skilled in the art that changes and modifications may be made, and therefore, the aim in the appended claims is to cover all such changes and modifications.