Wet gas flow measuring method and apparatus

Abstract

A wet gas flow measuring method, wherein measuring total flow differential pressure value P of wet gas in a pipeline by a differential pressure flow measuring device (201), measuring section gas contents of the wet gas in the pipeline by at least two phase fraction meters respectively (202), obtaining optimized section gas content value GVFopt by a flow calculating module based on the section gas contents respectively measured by the at least two phase fraction meters (203); and calculating gas volume flow rate Qg and liquid volume flow rate Ql by the flow calculating module based on the total flow differential pressure value P of the wet gas and the optimized section gas content value GVFopt (204). As the section gas content of the wet gas in the pipeline is detected by the redundant phase fraction meters, the gas volume flow rate Qg and the liquid volume flow rate Ql can be measured accurately, which meets the requirements on production measurements of oil and gas field and facilitates management improvement and production optimization of oil-gas reservoir. A wet gas flow measuring apparatus is also provided.

Claims

1. A method for measuring wet gas flow, characterized in that the method comprises the steps of: (a) measuring the total flow differential pressure value P of wet gas in a pipeline via a differential pressure type flow measuring device disposed on the pipeline; (b) measuring gas void fraction of the wet gas in the pipeline via at least two phase fraction meters disposed on the pipeline downstream from the differential pressure type flow measuring device, respectively; (c) obtaining optimized gas void fraction GVF.sub.opt by a flow calculating module based on the gas void fractions respectively measured by the at least two phase fraction meters; and (d) calculating the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l by the flow calculating module based on the total flow differential pressure value P and the optimized gas void fraction GVF.sub.opt of the wet gas; characterized in that: the step (c) comprises the following steps: (c1) calculating the short-term local reliability function F( ) of the gas void fraction Xi measured by each phase fraction meter: F(X.sub.i)=D(X.sub.i)=E([X.sub.iE(X.sub.i)].sup.2), wherein Xi denotes the time sequence of the gas void fraction measured by the phase fraction meter i, 1iN; N denotes the number of phase fraction meters; D( ) denotes a variance function; and E( ) denotes an expectation function; (c2) constructing a long-term whole reliability function LF( ) for the gas void fraction Xi measured by each phase fraction meter:
S.sub.i[0]=X.sub.i[0],
S.sub.i[m]=S.sub.i[m1]+(1)X.sub.i[m],
LF(S.sub.i)=D(S.sub.i), wherein 01, and 0mM1, where M denotes the number of the gas void fractions included in the time sequence, and X.sub.i[m] denotes the m.sup.th gas void fraction in the time sequence of the gas void fraction measured by the phase fraction meter i; (c3) constructing a comprehensive reliability function NICE( ) for the gas void fraction Xi measured by each phase fraction meter by utilizing the short-term local reliability function F( ) and the long-term whole reliability function LF( ): $\mathrm{NICE}\left(X_i\right)=\frac{1}{F\left(X_i\right)+LF\left(S_i\right)};$ and (c4) obtaining the optimized gas void fraction GVF.sub.opt by utilizing the comprehensive reliability function NICE( ) $GV{F}_{\mathrm{opt}}\left[m\right]=\frac{\underset{i=1}{\overset{N}{.Math.}}{X}_i\left[m\right]\mathrm{NICE}\left(X_i\right)}{\underset{i=1}{\overset{N}{.Math.}}\mathrm{NICE}\left(X_i\right)},$ wherein GVF.sub.opt[m] denotes the m.sup.th optimized value in the time sequence of the optimized gas void fraction GVF.sub.opt.

2. The method according to claim 1, characterized in that: the axial lines of said at least two phase fraction meter are respectively intersected with and perpendicular to the axial direction of the pipeline.

3. The method according to claim 1, characterized in that: the step (d) comprises the following steps: (d1) calculating the mixed density .sub.mix of the wet gas according to the optimized gas void fraction GVF.sub.opt; (d2) calculating the total volume flow rate Q of the wet gas according to the total flow differential pressure value P and the mixed density .sub.mix; and (d3) calculating the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l according to the total volume flow rate Q and the optimized gas void fraction GVF.sub.opt.

4. The method according to claim 3, characterized in that: the step (d1) comprises utilizing the equation .sub.mix=.sub.gasGVF.sub.opt+.sub.liquid (1GVF.sub.opt)to calculate the mixed density .sub.mix of the wet gas, wherein .sub.gas denotes the gas density, and .sub.liquid denotes the liquid density.

5. The method according to claim 3, characterized in that the step (d2) comprises utilizing the equation $Q=K\sqrt{\frac{P}{{}_{\min }}}$ to calculate the flow volume flow rate Q of the wet gas, wherein the parameter K is a systematic parameter.

6. The method according to claim 3, characterized in the step (d3) comprises the steps of: calculating the gas volume flow rate Q.sub.g by utilizing the equation Q.sub.g=QGVF.sub.opt; and calculating the liquid volume flow rate Q.sub.l by utilizing the equation Q.sub.g=Q(1GVF.sub.opt).

7. A wet gas flow measuring apparatus, comprising a pipeline, a differential pressure type flow measuring device, at least two phase fraction meters, and a flow calculating module, wherein the differential pressure type flow measuring device is disposed on the pipeline and said at least two phase fraction meters are disposed on the pipeline downstream from the differential pressure type flow measuring device and wherein: the pipeline is used to convey the wet gas; the differential pressure type flow measuring device is used to measure the total flow differential pressure value P of wet gas in the pipeline; the at least two phase fraction meters are respectively used to measure the gas void fractions of wet gas in the pipeline; and the flow calculating module is used to calculate the optimized gas void fraction GVF.sub.opt based on the gas void fractions respectively measured by the at least two phase fraction meters; and the flow calculating module is used to calculate the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l based on the total flow differential pressure value P and the optimized gas void fraction GVF.sub.opt; characterized in that the flow calculating module calculates the short-term local reliability function F( )of the gas void fraction Xi measured by each phase fraction meter: F(X.sub.i)=D(X.sub.i)=E([X.sub.iE(X.sub.i)].sup.2), wherein Xi denotes the time sequence of the gas void fraction measured by the phase fraction meter i, 1iN; N denotes the number of phase fraction meters; D( ) denotes a variance function; and E( )denotes an expectation function; the flow calculating module constructs a long-term whole reliability function LF( )for the gas void fraction Xi measured by each phase fraction meter:
S.sub.i[0]=X.sub.i[0],
S.sub.i[m]=S.sub.i[m1]+(1)X.sub.i[m],
LF(S.sub.i)=D(S.sub.i), which 01, and 0mM1, where M denotes the number of the gas void fractions included in the time sequence, and X.sub.i[m] denotes the m.sup.th gas void fraction in the time sequence of the gas void fraction measured by the phase fraction meter i; the flow calculating module constructs a comprehensive reliability function NICE( )of the gas void fraction Xi measured by each phase fraction meter by utilizing the short-term local reliability function F( )and the long-term whole reliability function LF( ) wherein $\mathrm{NICE}\left(X_i\right)=\frac{1}{F\left(X_i\right)+LF\left(S_i\right)};$ and the flow calculating module utilizes the comprehensive reliability function NICE( )to obtain the optimized gas void fraction GVF.sub.opt: $GV{F}_{\mathrm{opt}}\left[m\right]=\frac{\underset{i=1}{\overset{N}{.Math.}}{X}_i\left[m\right]\mathrm{NICE}\left(X_i\right)}{\underset{i=1}{\overset{N}{.Math.}}\mathrm{NICE}\left(X_i\right)},$ wherein GVF.sub.opt[m] denotes the m.sup.th optimized value in the time sequence of the optimized gas void fraction GVF.sub.opt.

8. The wet gas flow measuring apparatus according to claim 7, characterized in that the axial lines of said at least two phase fraction meters are intersected and perpendicular to the axial direction of the pipeline.

9. The wet gas flow measuring apparatus according to claim 7, characterized in that the flow calculating module calculates the mixed density mix of wet gas based on the optimized gas void fraction GVF.sub.opt; and, the flow calculating module calculates the total volume flow rate Q of wet gas based on the total flow differential pressure value P and the mixed density .sub.mix of wet gas; and, the flow calculating module calculates the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l based on the total volume flow rate Q and the optimized gas void fraction GVF.sub.opt.

10. The wet gas flow measuring apparatus according to claim 9, characterized in that the flow calculating module utilizes the equation .sub.mix=.sub.gasGVF.sub.opt+.sub.liquid (1GVF.sub.opt) to calculate the mixed density .sub.mix of the wet gas, wherein .sub.gas denotes the gas density, and .sub.liquid denotes the liquid density.

11. The wet gas flow measuring apparatus according to claim 9, characterized in that the flow calculating module utilizes the equation $Q=K\sqrt{\frac{P}{{}_{\min }}}$ to calculate the total volume flow rate Q, wherein the parameter K is a systematic parameter.

12. The wet gas flow measuring apparatus according to claim 9, characterized in that the flow calculating module utilizes the equation Qg=QGVFopt to calculate the gas volume flow rate Qg and utilizes the equation Ql=Q(1GVFopt) to calculate the gas volume flow rate Ql.

13. The wet gas flow measuring apparatus according to claim 7, characterized in that said at least two phase fraction meters are installed on the same one section of the pipeline.

14. The wet gas flow measuring apparatus according to claim 7, characterized in that said at least two phase fraction meters are installed on different sections of the pipeline.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a diagram which describes a scheme for measuring wet gas flow in the prior art.

(2) FIG. 2 is a diagram which describes an example for the wet gas flow measuring method of the invention.

(3) FIG. 3 is a diagram which describes an example for calculating the gas volume flow rate and the liquid volume flow rate of the invention.

(4) FIG. 4 is a diagram which describes an example for calculating the optimized gas void fraction of the invention.

(5) FIG. 5 is a diagram which describes another example for the wet gas flow measuring method of the invention.

(6) FIG. 6 is a diagram which describes an example for the wet gas flow measuring apparatus of the invention.

(7) FIG. 7 is a diagram which describes another example for the wet gas flow measuring apparatus of the invention.

(8) FIG. 8 is a diagram which describes an example for horizontal arrangement of the wet gas flow measuring apparatus of the invention.

(9) FIG. 9 is a diagram which describes an example for vertical arrangement of the wet gas flow measuring apparatus of the invention.

SUMMARY OF THE INVENTION

(10) The technical problem to be solved by the invention is to provide a method for measuring wet gas flow and apparatus thereof, i.e., as the gas void fraction of the wet gas in the pipeline is detected by a redundant phase fraction meter, the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l can be measured accurately, which meets the requirements on production measurements of oil and gas field and facilitates management improvement and production optimization of oil-gas reservoir.

(11) An aspect of the invention is to provide a method for measuring a wet gas flow, comprising the following steps: (a) measuring the total flow differential pressure value P of wet gas in a pipeline via a differential pressure type flow measuring device; (b) measuring gas void fraction of the wet gas in the pipeline via at least two phase fraction meters, respectively; (c) obtaining optimized gas void fraction GVF.sub.opt by a flow calculating module based on the gas void fractions respectively measured by the at least two phase fraction meters; and (d) calculating the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l by the flow calculating module based on the total flow differential pressure value P and the optimized gas void fraction GVF.sub.opt of the wet gas.

(12) Another aspect of the invention is to provide a wet gas flow measuring apparatus, comprising a pipeline, a differential pressure type flow measuring device, at least two phase fraction meters, and a flow calculating module, wherein the differential pressure type flow measuring device and said at least two phase fraction meters are respectively installed on the pipeline, and wherein: the pipeline is used to convey the wet gas; the differential pressure type flow measuring device is used to measure the total flow differential pressure value P of wet gas in the pipeline; the at least two phase fraction meters are respectively used to measure the gas void fractions of wet gas in the pipeline; and the flow calculating module is used to calculate the optimized gas void fraction GVF.sub.opt based on the gas void fractions respectively measured by the at least two phase fraction meters; and the flow calculating module is used to calculate the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l based on the total flow differential pressure value P and the optimized gas void fraction GVF.sub.opt.

(13) In accordance with the invention, a differential pressure type flow measuring device is used to measure the total flow differential pressure value of wet gas in a pipeline; at least two phase fraction meters are respectively used to measure the gas void fractions of wet gas in the pipeline; a flow calculating module, based on the gas void fractions respectively measured by the at least two phase fraction meters, can be used to calculate the optimized gas void fraction GVFopt, and the flow calculating module, based on the total flow differential pressure value P and the optimized gas void fraction GVFopt, can be used to calculate the gas volume flow rate Qg and the liquid volume flow rate Ql. As the gas void fraction of the wet gas in the pipeline is detected by a redundant phase fraction meter, the gas volume flow rate Qg and the liquid volume flow rate Ql can be measured accurately, which meets the requirements on production measurements of oil and gas field and facilitates management improvement and production optimization of oil-gas reservoir.

DETAILED DESCRIPTION OF THE INVENTION

(14) By referring to the drawings, more all-around descriptions to the invention are made as follow, where exemplary examples of the invention are set forth.

(15) FIG. 2 is a diagram which describes one example for the wet gas flow measuring method of the invention, and as shown in FIG. 2, the wet gas flow measuring method in the embodiment may comprise the following steps:

(16) Step 201, measuring the total flow differential pressure value P of wet gas in a pipeline via a differential pressure type flow measuring device;

(17) Step 202, measuring the gas void fraction of the wet gas in the pipeline via at least two phase fraction meters, respectively;

(18) Step 203, obtaining optimized gas void fraction GVF.sub.opt by a flow calculating module based on the gas void fractions respectively measured by the at least two phase fraction meters; and;

(19) Step 204, calculating the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l by the flow calculating module based on the total flow differential pressure value P and the optimized gas void fraction GVF.sub.opt of the wet gas.

(20) In accordance with the wet gas flow measuring method as provided in the above example of the invention, a differential pressure type flow measuring device is used to measure the total flow differential pressure value P of wet gas in a pipeline; at least two phase fraction meters are respectively used to measure the gas void fractions of wet gas in the pipeline; a flow calculating module, based on the gas void fractions respectively measured by the at least two phase fraction meters, can be used to calculate the optimized gas void fraction GVF.sub.opt, and the flow calculating module, based on the total flow differential pressure value P and the optimized gas void fraction GVF.sub.opt, can be used to calculate the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l. As the gas void fraction of the wet gas in the pipeline is detected by a redundant phase fraction meter, the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l can be measured accurately, which meets the requirements on production measurements of oil and gas field and facilitates management improvement and production optimization of oil-gas reservoir.

(21) It is preferred that the axial lines of the at least two phase fraction meters are respectively intersected with the axial direction of the pipeline and perpendicular thereto. Since the axial line of each phase fraction meter is respectively intersected with the axial direction of the pipeline and perpendicular thereto, the precision of measuring result can be further assured.

(22) It is preferred that in above step 204, the method as shown in FIG. 3 may be used to calculate the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l. A person skilled in the art will know that the method as shown in FIG. 3 merely is illustrative, and thus a person skilled in the art can use other alternative methods to take the calculation. As shown in FIG. 3, the method according to the invention for calculating the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l comprises the following steps:

(23) Step 301, calculating the mixed density .sub.mix of wet gas according to the optimized gas void fraction GVF.sub.opt;

(24) Step 302, calculating the total volume flow rate Q of wet gas according to the total flow differential pressure value P and mixed density .sub.mix of wet gas; and

(25) Step 303, calculating the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l according to the total volume flow rate Q and the optimized gas void fraction GVF.sub.opt.

(26) Preferably, in the step 301, the following equation can be used to calculate the mixed density .sub.mix of wet gas: .sub.mix=.sub.gasGVF.sub.opt+.sub.liquid (1GVF.sub.opt), wherein .sub.gas denotes the gas density, and .sub.liquid denotes the liquid density.

(27) Preferably, in the step 302, the following equation can be used to calculate the total volume flow rate Q of wet gas:

(28) $Q=K\sqrt{\frac{P}{{}_{\mathrm{mix}}}},$
wherein the parameter K is a systematic parameter.

(29) Preferably, in the step 303, the equation Q.sub.g=QGVF.sub.opt can be used to calculate the gas volume flow rate Q.sub.g; and the equation Q.sub.l=Q(1GVF.sub.opt) can be used to calculate the liquid volume flow rate Q.sub.l.

(30) A person skilled in the art will know that the specific equations used in the above steps 301-303 are merely illustrative, and thus a person skilled in the art can use other alternative equations to take the calculations.

(31) In the above step 203, the flow calculating module, with a method for acquiring a mean value, can calculate the optimized gas void fraction GVF.sub.opt by acquiring the mean value of the above gas void fractions which are respectively measured by at least two phase fraction meters. For example, the calculation can be carried out with a method for calculating an arithmetic mean value, a geometric mean volume, or a root mean square mean value.

(32) Preferably, in step 203, the flow calculating module can calculate the optimized gas void fraction GVF.sub.opt by a method for calculating a weighted mean value as shown in FIG. 4.

(33) FIG. 4 is a diagram which describes an example of the invention for calculating the optimized gas void fraction by utilizing a weighed means value. As shown in FIG. 4, the method for calculating the optimized gas void fraction may comprises the following steps:

(34) Step 401, calculating the short-term local reliability function F( ) of the gas void fraction Xi measured by each phase fraction meter: F (X.sub.i)=D (X.sub.i)=E([X.sub.iE(X.sub.i)].sup.2), wherein Xi denotes the time sequence of the gas void fractions measured by the phase fraction meter i, 1iN; N denotes the number of phase fraction meters; D( ) denotes a variance function; and E( ) denotes an expectation function;

(35) Step 402, constructing a long-term whole reliability function LEO for the gas void fraction Xi measured by each phase fraction meter:
S.sub.i[0]=X.sub.i[0],
S.sub.i[m]=S.sub.i[m1]+(1)X.sub.i[m],
LF(S.sub.i)=D(S.sub.i),
wherein 01, where this parameter a refers to a coefficient which denotes a weight; different weights may be set according to specific experimental conditions, and an optimized weight value can be found out by trial and error method; and 0mM1, where M denotes the number of the gas void fractions included in the time sequence, and X.sub.i[m] denotes the m.sup.th gas void fraction in the time sequence of the gas void fraction measured by the phase fraction meter i;

(36) Step 403, constructing a comprehensive reliability function NICE( ) for the gas void fraction Xi measured by each phase fraction meter by utilizing the short-term local reliability function F( ) and the long-term whole reliability function LF( ):

(37) $\mathrm{NICE}\left(X_i\right)=\frac{1}{F\left(X_i\right)+LF\left(S_i\right)}.$

(38) Step 404, obtaining the optimized gas void fraction GVF.sub.opt by utilizing the comprehensive reliability function NICE( ):

(39) $GV{F}_{\mathrm{opt}}\left[m\right]=\frac{\underset{i=1}{\overset{N}{.Math.}}{X}_i\left[m\right]\mathrm{NICE}\left(X_i\right)}{\underset{i=1}{\overset{N}{.Math.}}\mathrm{NICE}\left(X_i\right)},$
wherein GVF.sub.opt[.sub.m] denotes the m.sup.th optimized value in the time sequence of the optimized gas void fraction GVF.sub.opt.

(40) A person skilled in the art can know that in a specific example, the comprehensive reliability function NICE( ) is deemed as the weight.

(41) A person skilled in the art can know that the method as shown in FIG. 4 merely is illustrative, and the person skilled in the art can use other alternative methods to take the calculation. For example, when calculating a weight mean value, other methods can be used to obtain the weight value, and the forms of the weight value are not limited to it. Meantime, a weight value to which a signal having a high quality corresponds is relatively high.

(42) FIG. 5 is a diagram which describes another example for the wet gas flow measuring method of the invention. As compared with the example as shown in FIG. 2, the example as shown in FIG. 5 further comprises the step 501 of the slip correction for the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l after they are obtained in the step 204.

(43) It is preferred to use a calculating method for the precise slip analytic solution of a gas and liquid annular mist flow to conduct the slip correction to the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l. The method is described specifically as follows:

(44) In step 501, a slip correction module is used to correct the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l, so as to acquire a corrected gas volume flow rate Q.sub.g and a corrected liquid volume rate Q.sub.l:

(45) $Q_1^{}=\frac{Q_1}{\left(1-GV{F}_{\mathrm{opt}}+\mathrm{Slip}*GV{F}_{\mathrm{opt}}\right)},Q_g^{}=\frac{Q_g*\mathrm{Slip}}{\left(1-GV{F}_{\mathrm{opt}}+\mathrm{Slip}*GV{F}_{\mathrm{opt}}\right)},$
wherein the equation for calculating the slip factor is shown as follows:

(46) $\mathrm{Slip}=\frac{\left(1-d\right)\left(\frac{1}{{}_k^{}}\frac{d}{\mathrm{.Math.}}-2\left(\frac{d}{\mathrm{.Math.}}-1\right)\right)}{{\left(1-\frac{d}{\mathrm{.Math.}}\right)}^2+\left(1-\mathrm{.Math.}\right)\left(\frac{1}{{}_k^{}}{\left(\frac{d}{\mathrm{.Math.}}\right)}^2-2\left({\left(\frac{d}{\mathrm{.Math.}}\right)}^2-\frac{d}{\mathrm{.Math.}}\right)\right)},$
where [d, 1], d=(r.sup.2/R.sup.2), wherein r denotes the radius of the gas phase column; R denotes the radius of the pipeline; .sub.k=.sub.k+1, wherein .sub.k denotes the viscosity ratio of gas phase to liquid phase; denotes the gas volume fraction in the gas core of the annular mist flow; d denotes the gas volume fraction at the section of the pipeline. In the equation, when tends to d, the situation shows that the flow pattern of the wet gas is a pure mist flow without any liquid film, and at this time, there is no slip between gas phase and liquid phase and S=1. When tends to 1, the situation shows that the flow pattern of the wet gas is a liquid and gas annular flow. When the flow pattern of wet gas is found between the pure mist flow and the gas and liquid annular flow model, takes a certain value between d and 1, and the value can be experientially selected by a person skilled in the art.

(47) Since a difference between the velocities of gas phase and liquid phase (i.e., the slip) can result in some errors, a solution which is usually used in the art for solving the technical problem is to fit relevant experimental data so as to correct the gas flow rate and the liquid flow rate. However, the disadvantage of the correction resides in strong dependences of experiential models on experimental data and measuring conditions, so that the method cannot achieve the balance of the universality and the precision. The defect in the prior art can be overcome by utilizing the slip correction method provided by the present invention, i.e., the precise analytic solution for the slip of gas and liquid annular mist flow, to correct the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l, thereby to further increase measuring precisions.

(48) Preferably, in the example as shown in FIG. 5, the step 203 also can be carried out according to the example as show in FIG. 4, and the step 204 can be carried out according to the example as show in FIG. 3.

(49) FIG. 6 is a diagram which describes an example for the wet gas flow measuring apparatus of the invention. As show in FIG. 6, the wet gas flow measuring apparatus in the embodiment comprises a pipeline 601, a differential pressure type flow measuring device 602, at least two phase fraction meters 603, and a flow calculating module 604, wherein the differential pressure type flow measuring device 602 and said at least two phase fraction meters 603 are respectively installed on the pipeline, and wherein the pipeline 601 is used to convey the wet gas; the differential pressure type flow measuring device 602 is used to measure the total flow differential pressure value of wet gas in the pipeline; at least two phase fraction meters 603 are respectively used to measure the gas void fractions of wet gas in the pipeline; a flow calculating module 604, according to the gas void fractions respectively measured by the at least two phase fraction meters 603, can be used to calculate the optimized gas void fraction GVF.sub.opt, and the flow calculating module 604, according to the total flow differential pressure value P and the optimized gas void fraction GVF.sub.opt, can be used to calculate the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l.

(50) By utilizing the wet gas flow measuring apparatus provided in the above example in accordance with the invention, a differential pressure type flow measuring device is used to measure the total flow differential pressure value of wet gas in a pipeline; at least two phase fraction meters are respectively used to measure the gas void fractions of wet gas in the pipeline; a flow calculating module, according to the gas void fractions respectively measured by the at least two phase fraction meters, can be used to calculate the optimized gas void fraction GVF.sub.opt, and the flow calculating module, according to the total flow differential pressure value P and the optimized gas void fraction GVF.sub.opt, can be used to calculate the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l. As the gas void fraction of the wet gas in the pipeline is detected by a redundant phase fraction meter, the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l can be measured accurately, which meets the requirements on production measurements of oil and gas field and facilitates management improvement and production optimization of oil-gas reservoir.

(51) It is preferred that the axial lines of the at least two phase fraction meters are respectively intersected with the axial direction of the pipeline and perpendicular thereto. Since the axial line of each phase fraction meter is respectively intersected with the axial direction of the pipeline and perpendicular thereto, the precision of measuring result can be further assured.

(52) Preferably, said at least two phase fraction meters may either be installed on the same one section, or be installed on different sections of the pipeline.

(53) Preferably, the flow calculating module 604, in particular, according to the optimized gas void fraction GVF.sub.opt, can calculate the mixed density .sub.mix of wet gas; it, according to the total flow differential pressure value P and the mixed density .sub.mix of wet gas, can calculate the total volume flow rate Q of wet gas; and it, according to the total volume flow rate Q and the optimized gas void fraction GVF.sub.opt, can calculate the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l.

(54) Preferably, the flow calculating module 604 specifically uses the following equation .sub.mix=.sub.gasGVF.sub.opt+.sub.liquid (1GVF.sub.opt) to calculate the mixed density .sub.mix of wet gas, wherein .sub.gas denotes the gas density, and .sub.liquid denotes the liquid density.

(55) Preferably, the flow calculating module 604 specifically uses the equation

(56) $Q=K\sqrt{\frac{P}{{}_{\mathrm{mix}}}}$
to calculate the total volume flow rate Q of wet gas, wherein the parameter K is a systematic parameter.

(57) Preferably, the flow calculating module 604 uses the equation Q.sub.g=QGVF.sub.opt to calculate the gas volume flow rate Q.sub.g, and uses the equation Q.sub.g=Q(1GVF.sub.opt) to calculate the liquid volume flow rate Q.sub.l.

(58) The flow calculating module 604, with a method for acquiring a mean value, can calculate the optimized gas void fraction GVF.sub.opt by calculating the mean value of the above gas void fractions which are respectively measured by at least two phase fraction meters. For example, the calculation can be carried out with a method for calculating an arithmetic mean value, a geometric mean value, or a root mean square mean value, etc.

(59) Preferably, the flow calculating module 604 further can utilize a method for calculating a weighed mean value to calculate the optimized gas void fraction GVF.sub.opt. The flow calculating module 604 specifically calculates the short-term local reliability function F( ) of the gas void fraction Xi measured by each phase fraction meter: F (X.sub.i)=D (X.sub.i)=E ([X.sub.iE(X.sub.i)].sup.2), wherein X.sub.i denotes the time sequence of the gas void fraction measured by the phase fraction meter i, 1iN; N denotes the number of the phase fraction meters; D( ) denotes a variance function; and E( ) denotes an expectation function.

(60) The flow calculating module 604 specifically constructs a long-term whole reliability function LEO for the gas void fraction X.sub.i measured by each phase fraction meter:
S.sub.i[0]=X.sub.i[0]
S.sub.i[m]=S.sub.i[m1]+(1)X.sub.i[m]
LF (S.sub.i)=D(S.sub.i), which 01,
where this parameter refers to a coefficient which denotes a weight; different weights may be set according to specific experimental conditions, and an optimized weight value can be found out by trial and error method; 0mM1, wherein M denotes the number of the gas void fractions included in the time sequence, and X.sub.i[m] denotes the m.sup.th gas void fraction in the time sequence of the gas void fraction measured by the phase fraction meter i.

(61) The flow calculating module 604 constructs a comprehensive reliability function NICE( ) of the gas void fraction X.sub.i measured by each phase fraction meter by utilizing the short-term local reliability function F( ) and the long-term whole reliability function LE( ):

(62) $\mathrm{NICE}\left(X_i\right)=\frac{1}{F\left(X_i\right)+LF\left(S_i\right)}.$

(63) The flow calculating module 604 specifically utilizes the comprehensive reliability function NICE( ) to obtain the optimized gas void fraction GVF.sub.opt:

(64) $GV{F}_{\mathrm{opt}}\left[m\right]=\frac{\underset{i=1}{\overset{N}{.Math.}}{X}_i\left[m\right]\mathrm{NICE}\left(X_i\right)}{\underset{i=1}{\overset{N}{.Math.}}\mathrm{NICE}\left(X_i\right)},$
wherein GVF.sub.opt[m] denotes the m.sup.th optimized value in the time sequence of the optimized gas void fraction GVF.sub.opt.

(65) FIG. 7 is a diagram which describes another example for the wet gas flow measuring apparatus of the invention. As compared with the example as shown in FIG. 6, the wet gas flow measuring apparatus further comprises a slip correction module 701 for the correction of the gas volume flow rate Q.sub.g and the liquid volume flow rate Q.sub.l so as to acquire a corrected gas volume flow rate Q.sub.g and a corrected liquid volume rate Q.sub.l:

(66) $Q_1^{}=\frac{Q_1}{\left(1-GV{F}_{\mathrm{opt}}+\mathrm{Slip}*GV{F}_{\mathrm{opt}}\right)},Q_g^{}=\frac{Q_g*\mathrm{Slip}}{\left(1-GV{F}_{\mathrm{opt}}+\mathrm{Slip}*GV{F}_{\mathrm{opt}}\right)},$
wherein the equation for calculating the slip factor is shown as follows:

(67) 0$\mathrm{Slip}=\frac{\left(1-d\right)\left(\frac{1}{{}_k^{}}\frac{d}{\mathrm{.Math.}}-2\left(\frac{d}{\mathrm{.Math.}}-1\right)\right)}{{\left(1-\frac{d}{\mathrm{.Math.}}\right)}^2+\left(1-\mathrm{.Math.}\right)\left(\frac{1}{{}_k^{}}{\left(\frac{d}{\mathrm{.Math.}}\right)}^2-2\left({\left(\frac{d}{\mathrm{.Math.}}\right)}^2-\frac{d}{\mathrm{.Math.}}\right)\right)}.$

(68) Where [d, 1], d=(r.sup.2/R.sup.2), wherein r denotes the radius of the gas phase column; R denotes the radius of the pipeline; .sub.k=.sub.k+1, wherein .sub.k denotes the viscosity ratio of gas phase to liquid phase; denotes the gas volume fraction in the gas core of the annular mist flow; d is the gas volume fraction at the section of the pipeline. In the equation, when tends to d, the situation shows that the flow pattern of the wet gas is a pure mist flow without any liquid film, and at this time, there is no slip between gas phase and liquid phase and S=1. When tends to 1, the situation shows that the flow pattern of the wet gas is a liquid and gas annular flow. When the flow pattern of the wet gas is found between the pure mist flow and the gas and liquid annular flow model, takes a certain value between d and 1, and the value can be experientially selected by a person skilled in the art.

(69) It is preferred that the wet gas flow measuring apparatus may be horizontally or vertically installed. FIG. 8 is a diagram which describes an example for horizontal arrangement of the wet gas flow measuring apparatus of the invention. In FIG. 8, the arrow represents the flowing direction of wet gas; sign 801 denotes the pipeline; sign 802 denotes a differential pressure type flow measuring device; sign 803 denotes the phase fraction meter (for the aim of simplifying the figure, merely two phase fraction meters are illustrated). FIG. 9 is a diagram which describes an example for vertical arrangement of the wet gas flow measuring apparatus of the invention. In FIG. 9, the arrow represents the flowing direction of the wet gas; the sign 901 denotes the pipeline; sign 902 denotes a differential pressure type flow measuring device; sign 903 denotes the phase fraction meter (for the aim of simplifying the figure, merely two phase fraction meters are illustrated).

(70) Preferably, on the pipeline, a pressure transmitter and a temperature transmitter are further installed for measuring pressure and temperature in the pipeline respectively, thereby to facilitate management improvement and production optimization of oil-gas reservoir.

(71) The descriptions to the invention are made aimed for exemplifying and describing the invention, but not for exclusive or limiting the invention to the disclosures. Thus, many modification and variations are obvious for a person skilled in the art. To select and describe the examples is aimed to better set forth the theory and actual applications of the invention, and so, a person skilled in the art can understand the invention so that various examples with various modifications can be designed to be adapted for certain uses.