CN102116755B - Method for measuring multiphase flow based on multi-section impedance type long-waist inner core and related speed measurement - Google Patents

Abstract

The invention belongs to the technical field of fluid measurement, and relates to a method for measuring a multiphase flow based on a multi-section impedance type long-waist inner core and related speed measurement. An adopted sensor comprises a cone throttling member, a shrinkage pressure drop differential pressure transmitter, a permanent pressure loss differential pressure transmitter and an electrical sensitive sensor. The method comprises the following steps of: acquiring two pairs of differential pressure signals generated by the cone throttling member; acquiring an electrical signal by using the electrical sensitive sensor; calculating volume flow of mixed fluid according to a correlation method; measuring data according to the differential pressure signals and the electrical signal; classifying, identifying and prejudging a flow pattern by using a support vector machine; measuring dielectric constant of the mixed fluid by using a capacitive model; measuring conductivity of the mixed fluid by using a conductivity model; and calculating a procedure parameter of an oil-gas-water flow. The method is not influenced in spite of conductivity of a continuous phase of the multiphase flow, and has the characteristics of high measuring accuracy, high reliability, more acquirable measuring information and the like.

Description

Multiphase flow measurement method based on long waist inner cone of multi-section impedance type and correlation speed measurement

Technical field

The invention belongs to the fluid measurement technical field, be specifically related to a kind of measuring method that merges differential pressure measurement principle and electrical measurement principle, be used for confirming that each phase component of polyphasic flow in external phase is from conductive phase to non-conductive phase scope contains parameters such as rate, flow velocity and flow.The present invention is measured as description object with oil/gas/water multiphase, but is not limited in this application, and measuring method of the present invention still is suitable in other industrial processs and chemical reaction.

Technical background

Polyphasic flow is often referred to and has flowing of two or more material simultaneously, comprises biphase gas and liquid flow, Dual-Phrase Distribution of Gas olid, liquid fixed double phase flow, liquid-liquid diphasic flow and gas-liquid liquid and gas-liquid-solid multiphase flow etc.Polyphasic flow often appears in modern project field such as power, nuclear energy, chemical industry, oil, refrigeration, metallurgy, pipeline transportation, aerospace, medicine, food and the equipment, as not realizing just can't guaranteeing the accurate measurement of its flow process parameter the reliable design and the safe operation of relevant devices in the pipeline.The popularity that polyphasic flow exists in nature and industrial process and the importance of application make it in national economy and daily life, occupy crucial status.Along with power industry and petrochemical industry develop rapidly, and, impel the work of polyphasic flow area research to develop rapidly, become at present and give the front subject greatly paid close attention to both at home and abroad the pay attention to day by day of environmental protection.

Because each alternate interfacial effect and relative velocity of existing of polyphasic flow, phase interface all is random variation on time and space, causes the flow characteristics of polyphasic flow complicated more than single-phase flow, and characteristic parameter is also many than single-phase flow.The oil gas water three phase flow measuring method can be divided into two aspects at present: 1, regard polyphasic flow as measurement that biphase gas and liquid flow carries out flow.2, each phase component in the polyphasic flow is contained the measurement of rate.The method of wherein regarding polyphasic flow as gas-liquid two-phase flow measurement flow mainly contains: correlation method, volumetric method, flow limit method, turbine type detection method, laser Doppler method, PIV particle imaging velocimetry, hot line hot film anemometer, process tomographic imaging technology, nuclear magnetic resonance method and Coriolis method etc.; This several method has wide range of applications, but needs other subsidiary means (like densitometer etc.) to measure void fraction.The measuring method of phase content has two kinds, and the one, indirect method promptly returns through experiment or theoretical derivation draws the phase content calculating formula that physical parameter and technological parameter with the gas-liquid two-phase medium are function; The 2nd, the measurement mechanism direct method of measurement is like Electromagnetic Wave Detection method, conductance method, capacitance method, calorifics method, microwave method, γ component table method and nuclear magnetic resonance method etc.In actual measurement, generally adopt the measurement mechanism direct method of measurement to confirm phase content.

Using that traditional single-phase flow measurement instrument combines with two-phase flow auxiliary parameter measurement means is a kind of main thought that solves polyphasic flow procedure parameter problems of measurement in the pipeline.Differential pressure flowmeter like orifice plate, Venturi tube, nozzle and inner cone flow meter etc., is to use maximum online flow measurement instruments of a kind of industry at present, and is also commonplace in the research of polyphasic flow procedure parameter field of measurement.Wherein inner cone flow meter notion that the fluid throttling is retracted to pipeline center's near axis fundamentally change into utilize the coaxial cone that is installed in pipeline center with fluid gradually throttling be retracted to inner-walls of duct, confirm fluid flow through differential pressure before and after measuring cone.Interior cone flow meter can be under shorter straight length condition; With wideer range ratio clean or dirty fluid is realized more accurate, more effective flow measurement; But also there are some shortcomings in himself design feature: there is pulsation in the pressure signal to exist the whirlpool district to cause here at the cone afterbody; Measuring accuracy and repeatability are reduced, and whole crushing is bigger, like the flow performance according to polyphasic flow inner cone flow meter is done malformation; Like inner venturi tube, spindle flowmeter and Shuttle Type Flow Meter etc., can effectively improve its precision and stability in multiphase flow measurement.Based on the multiphase flow measurement method of the responsive principle of electricity existing theory and application foundation preferably; Electrical parameter (electrical impedance, specific inductive capacity etc.) through the measuring channel inner fluid changes the measurement that each phase component of realization polyphasic flow contains rate, and this method is divided into mode such as condenser type, conductance type according to the difference of measurement parameter.If adopt multi-modal sensor and multi-operation mode (excitation measurement strategies) combined method, utilize the measurement advantage of each mode, can realize the comprehensive extraction of polyphasic flow procedure parameter through Combinatorial Optimization, and the acquisition all-sidedly and accurately of flowing information.

Patent CN99257650.4, CN200420093532.6, CN200420061026.9 have announced a kind of special-shaped Venturi tube, spindle flowmeter, Shuttle Type Flow Meter that is used for restriction device respectively, be overcome in the cone flow meter some shortcomings design.The interior relatively cone flow meter of special-shaped Venturi tube is at waist; Increase the equal diameter section; Form annular free air space, help to improve the stability of efflux coefficient, but it is short and angle is big that kinetic energy is converted into the conical diffuser of outlet of potential energy; Fluid produces whirlpool through separating after the cone, and the interior relatively cone flow meter of permanent crushing does not have obvious improvement.Spindle flowmeter and Shuttle Type Flow Meter are very approaching; Only difference to some extent on appearance profile: the throttling element of spindle flowmeter is fairshaped spindle; And the fusiformis body that the throttling element of Shuttle Type Flow Meter is made up of precentrum, middle equal diameter segment body and posterior pyramids; Throttling element is installed along conduit axis, and the centre has one section sufficiently long equal diameter section, and forms uniform circular passage between the inner-walls of duct.These two kinds of deformable bodys are in efflux coefficient stability; Anti-interference fluidity and low pressure loss property all are superior to interior cone flow meter; But throttling element long (4 times to internal diameter of the pipeline); When caliber is big, seeming, comparison is heavy, and it is single relatively to measure the polyphasic flow procedure parameter, needs to cooperate the measurement that other flowmeters could the complete process parameters.

Patent CN86207384, CN200810151346.6, CN03234526.7 have announced " mixer+double-venturi tube " moisture measuring device that Venturi tube restriction device and throttle orifice plate apparatus combination, inner cone restriction device and the combination of Venturi tube restriction device, the combination of two inner cone restriction device and Britain Solartron company propose respectively; Adopt two differential pressure methods to measure gas-liquid two phase contents and flow, i.e. two kinds of similar or inhomogeneous restriction device multiple measurement biphase gas and liquid flow procedure parameters.Utilize the homogeneous flow model; The differential pressure signal that obtains on the restriction device to the different flow coefficient is handled; Obtain the gas phase quality and contain rate, obtain gas-liquid separate phase flow rate in the fluid-mixing by the total mass flow calculation of fluid-mixing then, simultaneously the gas-liquid flow is carried out the temperature and pressure compensation.This combined method only is applicable to measure mixes two-phase fluid uniformly; And two restriction devices fit together the consistance that can not guarantee on the space; Therefore two-phase flow flow pattern and procedure parameter in flow process is difficult to be consistent; Run counter to easily that combined type measures precondition---promptly there is original reason error in the constant principle of its basic parameter when two-phase fluid flows through each ingredient.

Patent US 5485743, US 5130661, US 6857323, US 6940286, CN99209502.6, CN200680055193.4 etc. have announced employing process tomographic imaging technology measurement polyphasic flow procedure parameter.Process tomographic imaging is made up of the sensor space array, data acquisition and the processing unit that obtain the measured object field information, image reconstruction and thing field parameter extraction unit; The sensor array is set up responsive successively in measurement space under certain mode of operation; Through data acquisition and processing unit the signal that reflects detected material field dielectric distribution in the different sensitizing ranges is gathered and handled; Use image algorithm to reconstruct the image of reflection dielectric distribution then, adopt certain information processing method to extract required parameter on this basis.Yet chromatography imaging technique is unmet industrial application requirements still at aspects such as real-time, measuring accuracy and picture qualities; And the quick variation characteristic of polyphasic flow conflicts with the complicacy of tomography image reconstruction algorithm, may cause result of calculation to exist than mistake.

Summary of the invention

The objective of the invention is deficiency, provide a kind of and can measure the multiphase flow procedure parameter, and realize that based on each sensor own characteristic the signal between sensor is complementary and redundant, realize measuring optimized multiphase flow measurement method to above-mentioned prior art existence.The present invention can realize the polyphasic flow procedure parameter in external phase is from conductive phase to non-conductive phase scope; Accurate measurement like a minute phase content (cross section phase content, quality phase content), phase-splitting flow velocity, volumetric flow rate and mass rate; Carry out flow pattern identification through the analysis to measure signal, to optimize measurement result.Technical scheme of the present invention is following:

A kind of multiphase flow measurement method based on long waist inner cone of multi-section impedance type and correlation speed measurement; The sensor that adopts comprises segment of cone stream part, shrinks the pressure drop differential pressure transmitter; Electricity sensitive sensor, described cone throttling element are a kind of long waist inner cone restriction device, by the leading portion cone, constitute for three sections as stage casing column, the back segment cone of cone waist; Through the coaxial pipeline center of being fixed in of support, form annular free air space with the tube wall of pipeline; Two pressure ports of described contraction pressure drop differential pressure transmitter, one is arranged on cone throttling element upstream end, and another is arranged on annular free air space place or cone afterbody; Described electricity sensitive sensor comprises four groups of electrodes that are distributed on the annular free air space different cross section; There is the structure of two arrays of electrodes all identical in these four groups of electrodes with the spatial disposition mode; Be used to detect the conductivity of polyphasic flow, it is right to be called resistance electrode, and the structure of two arrays of electrodes is all identical with the spatial disposition mode in addition; Be used to detect the specific inductive capacity of polyphasic flow, it is right to be called capacitance electrode; Every group of electrode is made up of electrode in one or more and one or more external electrode, and interior electrode is fixed on cone waist surface, and external electrode is fixed on inner-walls of duct; Gating is positioned at electrode and its corresponding external electrode on the same cross section respectively as excitation/potential electrode and ground-electrode during measurement, and this measuring method comprises the following steps:

1) gathers two pairs of differential pressure signals that produce by the cone throttling element: shrink pressure drop differential pressure signal Δ p and permanent crushing differential pressure signal δ p;

2) utilize the electricity sensitive sensor to gather electric signal;

3) calculate the fluid-mixing volume flow Q according to cross-correlation method _v:

A, according to formula

Calculate the relevant speed V of liquid phase ₁Relevant speed V with gas phase ₂, in the formula, L ₁, L ₂Be respectively resistance electrode to the spacing between the two arrays of electrodes and capacitance electrode to the spacing between the two arrays of electrodes, τ ₀₁, τ ₀₂Be respectively resistance electrode to signal delay time and capacitance electrode to signal delay time.

B, according to formula V=K _lV ₁+ V _dAnd Q _v=A ₀V calculates the fluid-mixing volume flow Q _v, in the formula, V is the fluid-mixing mean flow rate, K ₁Be speed correction factor, V _dBe alternate sliding velocity, A ₀Circulation area for annular free air space;

4) according to differential pressure signal and electric signal measurement data, carry out classification, the identification and judgement of flow pattern, judge that external phase is conductive phase or non-conductive phase in the fluid;

5) whether be conductive phase according to external phase, adopt Different Calculating Models to calculate the procedure parameter of oil-gas-water three-phase flow,

If external phase is non-conductive phase, the computation model of employing is following:

A. according to the voltage signal of capacitance electrode, calculate the DIELECTRIC CONSTANTS of fluid-mixing to gathering _m

B. in the polyphasic flow with cross section oil content α _o, cross section water oil rate α _w, void fraction α _gEquivalence contains rate, known oil phase DIELECTRIC CONSTANTS for the volume of Oil, Water, Gas respectively _o, the water DIELECTRIC CONSTANTS _wAnd gas phase DIELECTRIC CONSTANTS _g, according to formula

And φ _o+ φ _w+ φ _g=1 calculates the oil phase volume contains rate φ _o, the water volume contains rate φ _wContain rate φ with gaseous phase volume _g, ε in the formula _mBe the specific inductive capacity of fluid-mixing, ρ _o, ρ _w, ρ _gBe oil phase density, water density, density of gas phase, C ₁For shrinking pressure drop Δ p is the efflux coefficient that differential pressure signal calibrates, and Δ p is for shrinking differential pressure;

C. according to formula

Calculate the oil phase volume flow Q _Vo, the water volume flow Q _VwWith gas phase volume flow rate Q _Vg

D. according to formula

Calculate oil phase flow velocity V _o, aqueous phase flow rate V _wWith gas phase flow velocity V _g, in the formula, A ₀Be the circulation area of annular free air space, α _oBe cross section oil content, α _wBe cross section water percentage, α _gBe void fraction;

If external phase is conductive phase, the computation model of employing is following:

A. according to the voltage signal of resistance electrode, calculate the conductivity of fluid-mixing to gathering _m

B. known water conductivity _w, according to formula

With formula α _l+ α _g=1 calculates the oil phase cross section contains rate α _o, the water cross section contains rate α _wContain rate α with the gas phase cross section _g, in the formula, σ _mBe profit gas mixed conductivity;

C. according to formula Calculate the oil phase volume flow Q _Vo, the water volume flow Q _VwWith gas phase volume flow rate Q _Vg

D. according to formula

Calculate oil phase flow velocity V _o, aqueous phase flow rate V _wWith gas phase flow velocity V _g

As preferred implementation, constitute the length range of each right electrode of capacitance electrode in the pipeline axial direction, promptly width range is 20mm-50mm, utilizes the cross correlation measurement method to measure the relevant speed V of gas phase ₁Constitute the length range of each right electrode of resistance electrode in the pipeline axial direction, promptly width range is 1mm-10mm, utilizes the relevant speed V of cross correlation measurement method measurement of liquid ₂

Beneficial effect of the present invention and advantage are following:

1, the measurement mechanism of the present invention's employing; The annular free air space that cone throttling element and duct wall form carries out rectification to the multiphase flow fluid-mixing; Flow in the circular passage of formation standard; Channel inner pressure is stablized and pressure is linear decline along axis, makes the repeatability of measurement and accuracy raising become possibility; Reliable mounting structure flexibly can be changed the cone throttling element of different-diameter ratio based on the different measuring object, with the scope of application of expansion measurement.

2, the present invention is through introducing flow pattern anticipation device or directly adopting electric sensor and the simultaneously-measured advantage of differential pressure (flow) sensor; Utilization contains rate and predicts the outcome and carry out classification, identification and the anticipation of flow pattern with the volume forecasting result; When adopting arc-shaped electrode like electric sensor; Also can utilize the differentiation that the rate distributed intelligence realizes instantaneous flow pattern that contains of its measurement, adopt the accurate measurement of Different Calculating Models implementation procedure parameter according to variations in flow patterns.

3, the polyphasic flow fluid-mixing flows through the cone throttling element and in annular free air space, produces and shrink pressure drop Δ p; And in the cone downstream pressure recovery place of certain distance produces permanent crushing δ p; Adopt contraction pressure drop Δ p to make differential pressure signal, comprised flowing information more fully with permanent crushing δ p; Utilize the simple crosscorrelation speed-measuring method to cooperate flow model, can obtain the two-phase flow speed in the annular circulating area more accurately based on polyphasic flow mechanism.

Description of drawings

Following accompanying drawing has been described embodiment chosen of the present invention, is exemplary drawings and non exhaustive or restricted, wherein:

The long waist inner cone of the multi-section impedance type multiphase flow measurement device front view that Fig. 1 measuring method of the present invention adopts;

The long waist inner cone of the multi-section impedance type multiphase flow measurement device lateral cross section view that this measuring method of Fig. 2 adopts;

Fig. 3 (a) is vertical C-C cross section view of the cone throttling element of the present invention's employing; (b) be the vertical view of the embodiment of the invention 1 cone; (c) be the horizontal B-B cross section view of cone;

The vertical view of second kind of sensor that Fig. 4 (a) measuring method of the present invention adopts (b) is its horizontal B-B cross section view; (c) be its vertical A-A cross section view; (d) be its horizontal C-C cross section view, wherein, interior electrode is four cross section ring electrodes, and dispatch from foreign news agency is four cross sections, 8 arc-shaped electrodes very;

The vertical view of the third sensor that Fig. 5 (a) measuring method of the present invention adopts (b) is its horizontal B-B cross section view; (c) be its vertical A-A cross section view; (d) be its horizontal C-C cross section view, wherein, interior electrode is four cross sections, 8 arc-shaped electrodes, and dispatch from foreign news agency is four cross sections, 8 arc-shaped electrodes very;

The vertical view of the 4th kind of sensor that Fig. 6 (a) measuring method of the present invention adopts (b) is its horizontal B-B cross section view; (c) be its vertical A-A cross section view; (d) be its horizontal C-C cross section view;

The vertical view of the 5th kind of sensor that Fig. 7 (a) measuring method of the present invention adopts (b) is its horizontal B-B cross section view; (c) be its vertical A-A cross section view; (d) be its horizontal C-C cross section view, wherein, interior electrode is single wide electrode, and dispatch from foreign news agency is four cross section same widths ring electrodes very;

The computing unit structural drawing that Fig. 8 adopts for the present invention;

Signal in Fig. 9 computing unit takes place and the mode selection module structural drawing;

Signal condition modular structure figure in Figure 10 computing unit;

AD in Figure 11 computing unit gathers and the signal pre-processing module structural drawing;

Figure 12 apparatus of the present invention measuring method process flow diagram;

Label declaration among the figure: the different cross section electrode is used a respectively, b, and c, d distinguishes, a among the label 16-21, b, c, d represent the 16-21 of different cross section electrode.

Embodiment

Below describing manufacturing in detail and operate step of the present invention, be intended to describe as embodiments of the invention, is not to be should be comprised within the scope of the invention other embodiment that can realize identical function by unique form of making or utilizing yet.

Specify the preferred embodiments of the present invention below in conjunction with Figure of description.

With reference to accompanying drawing 1-4; Two ends of measurement mechanism that the present invention adopts have the pipeline 1 of flange 2; 7, three pressure ports of a cone throttling element: upper reaches pressure port 3, middle pressure port 4,5, two transmitters of downstream pressure port: shrink pressure drop differential pressure transmitter 24, permanent crushing differential pressure transmitter 25; Be installed in the electricity sensitive sensor on annular free air space 12 corresponding pipelines 1 and the cone waist 9: external electrode 17, interior electrode 16, computing unit 27 and other accessory structures.Cone throttling element 7 be on the basis of traditional inner cone with throat downstream equal diameter lengthening form cone waist 9 and get; Cone throttling element 7 is through support 11 and the support fixture splice 6 coaxial pipe centerline places that are installed in; And the cone 7 that can change different throttling ratios according to measured matter and operating mode is with the expansion measurement range, and cone throttling element 7 is included in constant cone waist 9 and the diameter of entrance 8, diameter that the flow direction diameter increases gradually and is decreased to zero outlet section 10 gradually, when polyphasic flow fluid-mixing 13 flows through this cone throttling element 7; One of circulation area experience is dwindled gradually; And the process that enlarges gradually after in certain distance, stablizing, during polyphasic flow fluid-mixing 13 process cone entrances 8, circulation area reduces reposefully; Speed is increased regularly; Static pressure is corresponding to be reduced, and along with circulation area reduces, the rate of change of speed is big more; When fluid flows through cone waist 9 places, owing to inertia effect continues to quicken, can the fluid that continue to quicken be adjusted because of the rectified action of cone throttling element 7 again, flowed in its circular passage that forms standard at annular free air space 12 places.Cone throttling element 7 adopts resistant material (high performance metal material, industrial ceramics or other corrosion-resistant materials) in case the fluid corrosion.At cone throttling element 7 upper reaches; Annular free air space 12 and downstream pressure recover corresponding tube wall pressure port are set; Because cone throttling element 7 is coaxially installed in the pipeline 1;, fluid-mixing 13 can produce contraction pressure drop Δ p with middle pressure port 4 by pressure port 3 when flowing through at the upper reaches of cone throttling element 7; Pressure port 3 produces permanent crushing δ p with pressure recovery place downstream pressure port 5 at the upper reaches, shrinks pressure drop transmitter 24 and gathers these two pairs of differential pressure signals with permanent crushing transmitter 25 through pressure guiding pipe 23, and the electrical quantity signal that produces is transferred to computing unit 27 through transducer signal line 26.

Be embedded in cone waist 9 in the device that the present invention adopts and can in annular free air space 12, form symmetrical electric field with electricity sensitive sensor on the pipeline 1, through measurement fluid electrical parameter (as electricity lead, electric capacity, specific inductive capacity etc.) confirm that each phase (component) contains rate.The electricity sensitive sensor comprises four groups of electrodes that are distributed on the annular free air space different cross section; There is the structure of two arrays of electrodes all identical in these four groups of electrodes with the spatial disposition mode; Be used to detect the conductivity of polyphasic flow; The structure of two arrays of electrodes is all identical with the spatial disposition mode in addition, is used to detect the specific inductive capacity of polyphasic flow; Every group of electrode constitutes by interior electrode and external electrode, and interior electrode is cone throttling element itself or is fixed on cone waist surface that external electrode is fixed on inner-walls of duct.The electrode of two groups of same structures can be realized measurement of correlation.Among the figure, interior electrode 16 is embedded in cone waist 9 top layers, and external electrode 17 is embedded in pipeline 1 inwall top layer, and electrode and the corresponding top layer that is embedded into body maintain an equal level, and promptly electrode can not produce disturbance to polyphasic flow fluid-mixing 13.Electrode 16 is the metal ring electrode in four cross sections; The interior electrode resistance 16b, the 16d that comprise interior electrode 16a, 16c and the measurement fluid-mixing conductivity of measuring the fluid-mixing specific inductive capacity are linked to each other through separately circuit in the inner signal wire passage 22 of cone throttling element 7 and support 11 and computing unit 27 by terminals 18 separately, signal wire 20; Four cross section external electrodes 17 are the circular metal electrode with interior electrode 16 same widths and thickness; The external electrode 17b, the 17d that comprise external electrode 17a, 17c and the measurement fluid-mixing conductivity of measuring the fluid-mixing specific inductive capacity are linked to each other by separately circuit in terminals 19 separately, signal wire 21 and the computing unit 27.External electrode 17 is as exciting electrode in the present embodiment; Insulation and ground connection between interior electrode 16 and the cone throttling element 7 are measured the voltage signal between internal and external electrode, and measuring-signal is through signal wire 20,21; Be transferred to computing unit 27; Differential pressure signal simultaneous with above-mentioned obtains the polyphasic flow procedure parameter according to computation model, like a minute phase content (cross section phase content, quality phase content), phase-splitting flow velocity, volumetric flow rate and mass rate; To export signal 28 then and show perhaps teletransmission to the Long-distance Control interface, the measurement of completion total system.Measurement mechanism of the present invention is installed in the actual pipeline through flange, and the diameter of concrete pipeline and length is border situation and deciding factually.

Multi-section electrode form and array mode that the present invention adopts can have diversified characteristics, are 0 ° to 180 ° like electrode center of circle angular region; Electrode width scope 0mm is to L/4 (cone throttling element waist length L); Internal and external electrode can be that identical shaped electrode can be a difformity electrode etc. also, so the electrode structure of this device has various ways, and Fig. 2, Fig. 4-7 has described wherein five kinds of typical array modes, and other array configurations can make up on this basis.

Fig. 2 has described first kind of sensor construction that the present invention adopts: internal and external electrode is the ring electrode cross section view (preferred electrode structure) of same widths; One in the internal and external electrode as another ground connection of exciting electrode, and energisation mode is current excitation or voltage drive.

Fig. 4 has described second kind of sensor construction that the present invention adopts: dispatch from foreign news agency very is distributed in 8 arc-shaped electrodes on four cross sections; The electrode 17b, the 17d that comprise electrode 17a, 17c and the measurement fluid-mixing conductivity of measuring the fluid-mixing specific inductive capacity; Interior electrode is that same widths is distributed in four ring electrodes on the cross section; The electrode 16b, the vertical view (a) of 16d, transverse cross sectional view (b, d), the longitudinal profile view (c) that comprise electrode 16a, 16c and the measurement fluid-mixing conductivity of measuring the fluid-mixing specific inductive capacity; Wherein electrode 17a, 17c, 16a, 16c are for measuring the electrode of fluid-mixing specific inductive capacity; Width is 30mm, and electrode 17b, 17d, 16b, 16d are electrode, and width is 5mm.The structure of electrode, fixed form, signal wire etc. are consistent with interior electrode shown in the embodiment 1 in this array configuration; Each circular arc external electrode 17 all is connected computing unit 27 through terminals 19 separately with signal wire 21.The preferred central angle of external electrode is 11 °, can know according to experiment conclusion, is distributed in 4 °-40 ° of 8 arc-shaped electrode central angle ranges of choice on the cross section, and arc-shaped electrode number variable (central angle changes thereupon).One in the internal and external electrode as another ground connection of exciting electrode.Energisation mode is current excitation or voltage drive.

Fig. 5 dispatch from foreign news agency very is distributed in 8 arc-shaped electrodes on four cross sections; The electrode 17b, the 17d that comprise electrode 17a, 17c and the measurement fluid-mixing conductivity of measuring the fluid-mixing specific inductive capacity; Interior electrode is that same widths is distributed in 8 arc-shaped electrodes on four cross sections; The electrode 16b, the vertical view (a) of 16d, transverse cross sectional view (b, d), the longitudinal profile view (c) that comprise electrode 16a, 16c and the measurement fluid-mixing conductivity of measuring the fluid-mixing specific inductive capacity; Wherein electrode 17a, 17c, 16a, 16c are for measuring the electrode of fluid-mixing specific inductive capacity; Width is 30mm, and electrode 17b, 17d, 16b, 16d are for measuring the electrode of fluid-mixing conductivity, and width is 5mm.The internal and external electrode number is consistent in this array configuration; And the position is corresponding one by one; Electrode 16a, 16b, 16c, 16d in each circular arc are by terminals 18, signal wire 20 link to each other with computing unit 27 through the signal wire passages 22 in cone throttling element 7 and support 11 inside separately; Each circular arc external electrode 17a, 17b, 17c, 17d all are connected computing unit 27 through terminals 19 separately with signal wire 21.The preferred central angle of inside and outside circle arc-shaped electrode is 11 °, can know according to experiment conclusion, is distributed in 4 °-40 ° of 8 arc-shaped electrode central angle ranges of choice on the cross section, and arc-shaped electrode number variable (central angle changes thereupon).One in the internal and external electrode as another ground connection of exciting electrode.Energisation mode is current excitation or voltage drive.

Fig. 6 has described the 4th kind of sensor construction that the present invention adopts; Constitute two groups of right internal and external electrodes of capacitance electrode and be ring electrode; Constitute vertical view (a), transverse cross sectional view (b, d), longitudinal profile view (c) that two groups of right internal and external electrodes of resistance electrode are 8 arc-shaped electrodes; Wherein, the external electrode 17a, 17c and interior electrode 16a, the 16c that measure the fluid-mixing specific inductive capacity are ring electrode, and width is 30mm; The external electrode 17b, 17d and interior electrode 16b, the 16d that measure the fluid-mixing conductivity are 8 arc-shaped electrodes, and width is 5mm.The internal and external electrode number is consistent in this array configuration; And the position is corresponding one by one; The interior electrode 16b of electrode 16a, 16c and each circular arc, 16d in the annular are by terminals 18, signal wire 20 link to each other with computing unit 27 through the signal wire passages 22 in cone throttling element 7 and support 11 inside separately; Annular external electrode 17a, 17c and each circular arc external electrode 17b, 17d all are connected computing unit 27 through terminals 19 separately with signal wire 21.The preferred central angle of inside and outside circle arc-shaped electrode is 11 °, can know according to experiment conclusion, can know according to experiment conclusion, is distributed in 4 °-40 ° of 8 arc-shaped electrode central angle ranges of choice on the cross section, and arc-shaped electrode number variable (central angle changes thereupon).One in the internal and external electrode as another ground connection of exciting electrode.Energisation mode is current excitation or voltage drive.

Fig. 7 has described the 5th kind of sensor construction that the present invention adopts: dispatch from foreign news agency very is distributed in four same widths ring electrodes on the cross section; Comprise resistive electrodes 17b, the 17d of capacitive electrode 17a, 17c and the measurement fluid-mixing conductivity of measuring the fluid-mixing specific inductive capacity, interior electrode is that width is vertical view (a), transverse cross sectional view (b, d), the longitudinal profile view (c) of the single wide electrode 16 of L.In this array configuration, electrode 16 in the wide annular is linked to each other with computing unit 27 through the signal wire passages 22 in cone throttling element 7 and support 11 inside by terminals 18, signal wire 20; Annular external electrode 17a, 17b, 17c, 17d all are connected computing unit 27 through terminals 19 separately with signal wire 21.Dispatch from foreign news agency is exciting electrode very, and interior electrode grounding, energisation mode are current excitation or voltage drive.

Computing unit is shown in accompanying drawing 8～11: mainly forms by signal generation and mode selection module, signal condition module, AD collection and signal pre-processing module, computer bus and CPU or computing machine, and as shown in Figure 8.Wherein:

Signal generation and mode selection module concrete structure mainly are made up of circuit such as computer bus, logic control circuit, signal generation apparatus, constant voltage or constant current source and electrode channel selections shown in accompanying drawing 9.Computer bus is mainly accomplished communicating by letter and data transmission of data acquisition system (DAS) and CPU or computing machine; The logic control circuit major function is the working forms that circuit takes place to select with electrode channel control signal, realizes the selection of system incentive strategy and mode of operation, and and signal condition module and AD collection and signal pre-processing module between co-ordination; It is the pumping signal of the controlled various ways of realization frequency that signal generating circuit selects major function, for example single, double polarity pulse signal and sinusoidal signal etc.; Constant voltage or constant-current source circuit adopt voltage control voltage source (VCVS) or voltage control current source (VCCS) scheme, and the pumping signal that mainly signal generating circuit is produced converts amplitude and adjustable constant voltage or the current signal of phase place into; Electrode channel selects circuit to adopt the digital simulation switch chip, transmits passage according to working method and policy logic control circuit foundation excitation with measuring-signal.

Signal condition module concrete structure mainly is made up of electrode drive circuit, logic control circuit and signal amplification circuit shown in accompanying drawing 10.Electrode drive circuit adopts the analog drive circuit technology, improves the load capacity of system and reduces the influence of spurious signals such as contact resistance, electric capacity and inductance to measured signal; The section selection module is selected the digital simulation switch chip for use, realizes the multi-section alternation through logic control circuit; Amplifying circuit is made up of high-precision meter operational amplifier or programmable amplifier chip, but the induced signal that electrode (adjacent electrode) is upward small in the pair cross-section and the transmission signal of differential pressure transmitter amplify the raising signal to noise ratio (S/N ratio).The setting of the logical address of logic control circuit realization amplifying circuit and the selection output of enlargement factor, and the logic communication between responsible and AD collection and signal pre-processing module and signal generation and the mode selection module.

AD collection and signal pre-processing module structure are shown in accompanying drawing 11.Mainly form by AD converter array, digital signal processing and logic control circuit, data storage computer bus etc.The AD converter array circuit adopts the AD conversion chip to constitute, and can convert tested simulating signal and the differential pressure signal passed back through the signal condition module to digital signal.Digital signal processing and logic control circuit realize digital signal filtering, separate and be in harmonious proportion function such as storage control, and the data after will handling are transferred to CPU or computer module through computer bus.

The data acquisition system (DAS) course of work that the present invention adopts: at first CPU or computer module carry out initialization operation through computer bus to signal generation and mode selection module; Comprise the exciting current selection; Excitation frequency is selected and excitation/measurement strategies selection etc., and is realized by logic control circuit.Pumping signal is alternately set up responsive respectively through the cross section handover module in each cross section, circular passage; Change the differential pressure signal that electrode voltage changes and differential pressure transmitter is passed back that causes by dielectric distribution in the field and pass through the cable parallel transfer to the signal condition module; Electric signal after the conditioning gets into AD and gathers and signal pre-processing module; Accomplish analog to digital conversion and the digital signal after the conversion is carried out processing such as filtering and demodulation by digital signal processing and logic control circuit control AD chip; Effective information after the processing is conveyed into CPU or computing machine through the computer bus interface, combines to carry out the analysis and the processing of measurement data with measurement model.

Describe in the face of multiphase flow measurement method of the present invention down:

(1). gather two pairs of differential pressure signals that produce by the cone throttling element: shrink pressure drop differential pressure signal Δ p and permanent crushing differential pressure signal δ p;

(2). utilize the electricity sensitive sensor to gather electric signal;

(3). calculate fluid-mixing total volumetric flow rate Q according to cross-correlation method _v

Electric field in the measure annular free air space 12; The time varying signal that 16d/17d, 16b/17b and capacitance electrode is obtained 16c/17c, 16a/17a from resistance electrode has reflected that each phase component of fluid-mixing 13 contains rate and changes, and can disclose some hydrodynamic characteristics relevant with flow pattern.The flow pattern of supposing fluid 13 remained unchanged in short time period; Promptly satisfy " solidifying " flow model; When certain disturbance such as bubble in the biphase gas and liquid flow or water droplet; Pass through resistance electrode successively and during to 16c/17c, 16a/17a, can produce time varying signal x 16d/17d, 16b/17b and capacitance electrode ₁(t), y ₁(t) and x ₂(t), y ₂(t), utilize cross-correlation method to measure the relevant speed of discrete phase in the fluid 13, wherein cross correlation function R _Xy(τ) be:

$R_{\mathrm{xy}1}\left(\mathrm{\τ}\right)=\underset{T\→\∞}{\lim }\frac{1}{T}{\∫}_0^T{y}_1\left(t\right)x_1(t-\mathrm{\τ})\mathrm{dt},R_{\mathrm{xy}2}\left(\mathrm{\τ}\right)=\underset{T\→\∞}{\lim }\frac{1}{T}{\∫}_0^T{y}_2\left(t\right)x_2(t-\mathrm{\τ})\mathrm{dt}---\left(1\right)$

In the formula:

x ₁(t), x ₂(t) the expression resistance electrode is to 16d/17d and the capacitance electrode time varying signal to the 16c/17c generation;

y ₁(t), y ₂(t) the expression resistance electrode is to 16b/17b and the capacitance electrode time varying signal to the 16a/17a generation;

τ representes the independent variable of cross correlation function;

Therefore through coherent signal x ₁(t), y ₁(t) in the Fluid Computation 13 from upstream electrode to 16d/17d to the transit time of downstream electrode to 16b/17b, and from upstream electrode to 16c/17c to the transit time of downstream electrode to 16a/17a, also be related function R _Xy1(τ), R _Xy2(τ) the pairing time τ in the peak value of figure (maximal value) position ₀₁, τ ₀₂It is exactly the delay time T between two signals.The big hydrodynamic noise of the capacitance electrode of general large-size to being used in the test fluid 13; Like big bubble; The relevant speed that this electrode pair is measured can be thought the relevant speed of gas phase; The right three-dimensional resolving power of the resistance electrode of reduced size is sensitiveer, can measure in the fluid 13 the small flow noise that the random variation by the droplet spatial distribution state causes, and the relevant speed that this electrode pair is measured can be thought the relevant speed of liquid phase.So relevant speed V of liquid phase in the fluid-mixing 13 ₁Relevant speed V with gas phase ₂For:

$V_1=\frac{L_1}{{\mathrm{\τ}}_{01}}$ Or $V_2=\frac{L_2}{{\mathrm{\τ}}_{02}}---\left(2\right)$

In the formula, L ₁, L ₂Be respectively spacing and capacitive electrode between the conductance type electrode pair two arrays of electrodes to the spacing between the two arrays of electrodes, τ ₀₁, τ ₀₂Be respectively signal delay time and capacitive electrode between the conductance type electrode pair two arrays of electrodes to the signal delay time between the two arrays of electrodes.Because flowing of fluid 13 can not meet " solidifying " flow model hypothesis fully, always there are certain difference in relevant speed and real fluid cross section average velocity, therefore introduce speed correction factor K ₁, K ₂,

V=K _lV ₁+ V _dOr V=K ₂V ₂+ V _d(3)

Wherein:

K ₁, K ₂Be speed correction factor, V _dBe alternate sliding velocity, can unite definite by experimental study and experimental formula;

V is a fluid-mixing average velocity;

Then the fluid-mixing total volumetric flow rate is:

Q _v＝A ₀V (4)

In the formula:

Q _vBe the fluid-mixing total volumetric flow rate;

A ₀Circulation area for annular free air space.

(4). utilize Throttle Principle to calculate and divide phase content;

The differential pressure signal (measurement mechanism that the present invention adopts through differential pressure flowmeter; Have and shrink pressure drop differential pressure transmitter and two of permanent crushing differential pressure transmitters, what the present invention selected for use is the measuring-signal that shrinks the pressure drop differential pressure transmitter) with pipeline in the total volumetric flow rate formula put out the expression formula of the total volumetric flow rate of fluid in order:

$Q_v=\frac{\mathrm{\ϵ}\·C_1}{\sqrt{1-{\mathrm{\β}}^4}}{A}_0\sqrt{\frac{2\mathrm{\Δp}}{{\mathrm{\ρ}}_m}}---\left(5\right)$

For oil-gas-water three-phase flow, fluid hybrid density ρ _mCan be expressed as:

ρ _m＝α _oρ _o+α _wρ _w+α _gρ _g (6)

In the formula: C ₁For shrinking pressure drop Δ p is the efflux coefficient that differential pressure signal calibrates,

ρ _o, ρ _w, ρ _gBe oil phase density, water density, density of gas phase;

α _o, α _w, α _gBe cross section oil content, cross section water oil rate, void fraction;

Simultaneous formula (4), (5), (6) obtain hybrid density ρ _mAnd shrink the expression formula between the differential pressure Δ p:

${\mathrm{\α}}_o{\mathrm{\ρ}}_o+{\mathrm{\α}}_w{\mathrm{\ρ}}_w+{\mathrm{\α}}_g{\mathrm{\ρ}}_g=\frac{2{(\mathrm{\ϵ}\·C_1)}^2\mathrm{\Δp}}{V(1-{\mathrm{\β}}^4)}---\left(7\right)$

For the gas of oil product, water and known static pressure and the temperature confirmed, can think the oil phase density p _o, the water density p _wAnd density of gas phase ρ _gFor known, and known

α _o+α _w+α _g＝1 (8)

Therefore wherein contain rate and can draw fluid hybrid density ρ in the cross section of a certain phase if record oil gas water three phase _m

(5). analyze differential pressure signal and electric signal; Carry out classification, identification and the anticipation of flow pattern; Judging object fluid external phase is conductive phase (water) or non-conductive phase (gas phase or oil phase); Right according to the corresponding potential electrode of the different choice of external phase then, if external phase is gas phase or oil phase, utilize capacitance electrode that the voltage signal of 16a/17a and 16c/17c is calculated; If external phase is water then utilizes resistance electrode that the voltage signal of 16b/17b and 16d/17d is calculated.

The general step of flow pattern identification is: at first from measurement data, extract the eigenwert that can reflect variations in flow patterns, utilize sorting algorithm that eigenwert is sorted out identification then, such as SVMs, artificial neural network etc.The present invention is that example is carried out flow pattern identification with the support vector machine method, and concrete steps are following:

1, extracts image data eigenwert makeup time sequence in 5 seconds, from time series, extract three kinds of statistical nature composition characteristic vectors

$T_1=\left[\begin{array}{ccc}\stackrel{\‾}{x}& S& \mathrm{SK}\end{array}\right],$

Time series mean value

has been represented the average energy that electric signal or differential pressure signal comprised of electricity sensitive sensor collection:

$\stackrel{\‾}{x}=\frac{1}{n}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i---\left(9\right)$

Wherein n is the eigenwert quantity in the time series, x _iBe single V _mValue.

Seasonal effect in time series standard deviation (S), reacted the degree of scatter of measuring-signal:

$S=\sqrt{\frac{1}{n-1}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{(x_i-\stackrel{\‾}{x})}^2}---\left(10\right)$

Seasonal effect in time series gradient (SK), reacted the seasonal effect in time series inclined degree:

$\mathrm{SK}=\frac{n}{(n-1)(n-2)}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\left(\frac{x_i-\stackrel{\‾}{x}}{S}\right)}^3---\left(11\right)$

2, SVMs (SVM) is realized a between class distance maximization and a type interior distance minimization through setting up a lineoid, can be used as sorter and realizes flow pattern identification.The optimal classification function does

$f\left(x\right)=\mathrm{sgn}\{({\mathrm{\ω}}^{*}\·x)+b^{*}\}=\mathrm{sgn}\{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{a}_i^{*}{y}_{i}K(x_i\·x)+b^{*}\}---\left(12\right)$

In the formula, sgn () is-symbol function, ω ^*Be the normal vector of lineoid, b ^*Be the threshold value of classification,

Be optimum solution, y _iBe constraint condition, K (x _iX) be kernel function; Select kernel function according to actual needs.

3, with the input of proper vector, the data that collect are classified as SVM.The target of this patent flow pattern identification is water quality type flow pattern, makings type flow pattern and oily type flow pattern three major types, to SVMs of every kind of flow pattern structure, output y _i∈ 1,1}, (i=1,2,3), then the output of 3 SVMs is one 3 dimensional vectors, whether each representation in components sample to should flow pattern, (1 ,-1 ,-1) expression water quality type flow pattern, (1,1 ,-1) expression makings type flow pattern, the oily type flow pattern of (1 ,-1,1) expression.Select for use radially basic kernel function

wherein ζ get 0.05, the SVMs that trains can be realized the identification of convection.

(6). utilize electric method to calculate each phase content.

Introduce electrical parameter measuring method measurement cross section and contain rate, comprise two kinds of electrical methods among the present invention, the one, through capacitance electrode the voltage signal between 16a/17a or the 16c/17c is calculated the specific inductive capacity of fluid-mixing and calculates the cross section phase content; The 2nd, through resistance electrode the voltage signal between 16b/17b or the 16d/17d is calculated the conductivity of fluid-mixing and calculates the cross section phase content.Gather the electrical signal of the electricity sensitive sensor in each cross section simultaneously; According to the flow pattern of step (five) differentiate result's (external phase whether conduct electricity) selection wherein a kind of measured value of method carry out subsequent treatment, so whether the present invention is not conducted electricity by continuous phase of multiphase flow to influence.Concrete principle is following: the electrical quantity of fluid 13 (as electricity lead, electric capacity, specific inductive capacity etc.) depend on each mutually shared ratio of fluid-mixing; And can set up each phase content and the relation of mixing electrical quantity through some computation models, as the electrical quantity in the Boyle model can through in the fluid-mixing each minute the phase fluid specific inductive capacity function representation.Discrete phase 15 is to be distributed in the liquid phase of being made up of profit two phase fluid-mixings in the continuous gas phase 14, and then the specific inductive capacity of fluid-mixing 13 can use H-B (H-Bnon) formula to calculate:

$\sqrt{{\mathrm{\ϵ}}_m}={\mathrm{\φ}}_g\sqrt{{\mathrm{\ϵ}}_g}+(1-{\mathrm{\φ}}_g)({\mathrm{\φ}}_w\sqrt{{\mathrm{\ϵ}}_w}+{\mathrm{\φ}}_o\sqrt{{\mathrm{\ϵ}}_o})---\left(13\right)$

In the formula:

ε _o, ε _w, ε _g, ε _mBe respectively the specific inductive capacity of oil phase, water, gas phase, fluid-mixing;

φ _o, φ _w, φ _gThe volume that is respectively oil phase, water, gas phase contains rate;

The polyphasic flow middle section contains rate α and is difficult to direct measurement, and the average external volume of general available measuring section contains rate φ and replaces.Wherein the specific inductive capacity of oil phase and water can obtain through tabling look-up or demarcating; Because the fluid-mixing specific inductive capacity is proportional to the electric capacity between internal and external electrode; And the electric capacity between internal and external electrode is proportional to the voltage magnitude between internal and external electrode; Therefore can calculate the specific inductive capacity of fluid-mixing through measuring capacitance electrode to the voltage signal between 16a/17a or the 16c/17c, through simultaneous formula (7), (8), (13), and then the cross section that obtains oil gas water three phase contains rate α _o, α _w, α _gAnd total volumetric flow rate Q _v

When the water percentage in the fluid-mixing 13 was too high, the use of capacitance method received certain limitation, but electricity consumption this moment is led as tested electrical quantity.Because the gas phase conductivity is approximately zero, and liquid phase is a conducting medium, so the fluid mixed conductivity depends primarily on the phase content of liquid phase and distribution mutually.According to Maxwell resistivity approximation theory, then the conductivity of fluid-mixing 13 is:

${\mathrm{\σ}}_m=\frac{2{\mathrm{\α}}_l}{3-{\mathrm{\α}}_l}\frac{2{\mathrm{\α}}_w}{3-{\mathrm{\α}}_w}{\mathrm{\σ}}_w---\left(14\right)$

In the formula:

α _l, α _wBe liquid holdup and cross section water percentage;

σ _m, σ _wBe respectively oil gas water mixed conductivity and water conductivity;

Because the voltage magnitude between mixed conductivity and internal and external electrode is inversely proportional to, therefore, can calculates the conductivity of fluid-mixing, and then obtain cross section water percentage α through the voltage signal between measuring resistance electrode pair 16b/17b or the 16d/17d _wThrough simultaneous formula (7), (8), (14), and then the cross section that obtains oil gas water three phase contains rate α _o, α _w, α _gAnd total volumetric flow rate Q _v

The definition that contains rate according to volume can calculate each phase volume flow of oil gas water:

Q _vo＝φ _oQ _v

Q _vw＝φ _wQ _v (15)

Q _vg＝φ _gQ _v

In the formula:

φ _o, φ _w, φ _gFor oil phase, water, gaseous phase volume contain rate;

Q _Vo, Q _Vw, Q _Vg, Q _vBe oil phase, water, gas phase volume flow rate and total volumetric flow rate;

Each phase content is according to the phase velocity definition:

$V_o=\frac{Q_{\mathrm{vo}}}{{\mathrm{\α}}_o{A}_0}$

$V_w=\frac{Q_{\mathrm{vw}}}{{\mathrm{\α}}_w{A}_0}---\left(16\right)$

$V_g=\frac{Q_{\mathrm{vg}}}{{\mathrm{\α}}_g{A}_0}$

In the formula:

V _o, V _w, V _gBe oil phase flow velocity, aqueous phase flow rate, gas phase flow velocity;

α _o, α _w, α _gBe cross section oil content, cross section water percentage, void fraction.

Claims (2)

1. multiphase flow measurement method based on long waist inner cone of multi-section impedance type and correlation speed measurement; The sensor that adopts comprises segment of cone stream part, shrinks the pressure drop differential pressure transmitter; Electricity sensitive sensor, described cone throttling element are a kind of long waist inner cone restriction device, by the leading portion cone, constitute for three sections as stage casing column, the back segment cone of cone waist; Through the coaxial pipeline center of being fixed in of support, form annular free air space with the tube wall of pipeline; Two pressure ports of described contraction pressure drop differential pressure transmitter, one is arranged on cone throttling element upstream end, and another is arranged on annular free air space place or cone afterbody; Described electricity sensitive sensor comprises four groups of electrodes that are distributed on the annular free air space different cross section; There is the structure of two arrays of electrodes all identical in these four groups of electrodes with the spatial disposition mode; Be used to detect the conductivity of polyphasic flow, it is right to be called resistance electrode, and the structure of two arrays of electrodes is all identical with the spatial disposition mode in addition; Be used to detect the specific inductive capacity of polyphasic flow, it is right to be called capacitance electrode; Every group of electrode is made up of electrode in one or more and one or more external electrode, and interior electrode is fixed on cone waist surface, and external electrode is fixed on inner-walls of duct; Gating is positioned at electrode and its corresponding external electrode on the same cross section respectively as excitation/potential electrode and ground-electrode during measurement, and this measuring method comprises the following steps:

1) gathers two pairs of differential pressure signals that produce by the cone throttling element: shrink pressure drop differential pressure signal Δ p and permanent crushing differential pressure signal δ p;

2) utilize the electricity sensitive sensor to gather electric signal;

3) calculate the fluid-mixing volume flow Q according to cross-correlation method _v:

A, according to formula

Calculate the relevant speed V of liquid phase ₁Relevant speed V with gas phase ₂, in the formula, L ₁, L ₂Be respectively resistance electrode to the spacing between the two arrays of electrodes and capacitance electrode to the spacing between the two arrays of electrodes, τ ₀₁, τ ₀₂Be respectively resistance electrode to signal delay time and capacitance electrode to signal delay time.

B, according to formula V=K1 _V1+ _VD and Q _v=A ₀V calculates the fluid-mixing volume flow Q _v, in the formula, V is the fluid-mixing mean flow rate, K ₁Be speed correction factor, V _dBe alternate sliding velocity, A ₀Circulation area for annular free air space;

4) according to differential pressure signal and electric signal measurement data, carry out classification, the identification and judgement of flow pattern, judge that external phase is conductive phase or non-conductive phase in the fluid;

5) whether be conductive phase according to external phase, adopt Different Calculating Models to calculate the procedure parameter of oil-gas-water three-phase flow,

If external phase is non-conductive phase, the computation model of employing is following:

A. according to the voltage signal of capacitance electrode, calculate the DIELECTRIC CONSTANTS of fluid-mixing to gathering _m

B. in the polyphasic flow with cross section oil content α _o, cross section water oil rate α _w, void fraction α _gEquivalence contains rate for the volume of Oil, Water, Gas respectively, has added the oil phase DIELECTRIC CONSTANTS _o, the water DIELECTRIC CONSTANTS _wAnd gas phase DIELECTRIC CONSTANTS _g, according to formula

φ _o+ φ _w+ φ _g=1 calculates the oil phase volume contains rate φ _o, the water volume contains rate φ _wContain rate φ with gaseous phase volume _g, ε in the formula _mBe the specific inductive capacity of fluid-mixing, ρ _o, ρ _w, ρ _gBe oil phase density, water density, density of gas phase, C ₁For being the efflux coefficient that differential pressure signal calibrates to shrink pressure drop Δ p;

C. according to formula $\begin{array}{c}{Q}_{\mathrm{Vo}}={\mathrm{\φ}}_o{Q}_v\\ Q_{\mathrm{Vw}}={\mathrm{\φ}}_w{Q}_v\\ Q_{\mathrm{Vg}}={\mathrm{\φ}}_g{Q}_v\end{array}$ Calculate the oil phase volume flow Q _Vo, the water volume flow Q _VwWith gas phase volume flow rate Q _Vg

D. according to formula $\begin{array}{c}{V}_o=\frac{Q_{\mathrm{Vo}}}{{\mathrm{\α}}_o{A}_0}\\ V_w=\frac{Q_{\mathrm{Vw}}}{{\mathrm{\α}}_w{A}_0}\\ V_g=\frac{Q_{\mathrm{Vg}}}{{\mathrm{\α}}_g{A}_0}\end{array}$ Calculate oil phase flow velocity V _o, aqueous phase flow rate V _wWith gas phase flow velocity V _g, in the formula, A ₀Be the circulation area of annular free air space, α _oBe cross section oil content, α _wBe cross section water percentage, α _gBe void fraction;

If external phase is conductive phase, the computation model of employing is following:

A. according to the voltage signal of resistance electrode, calculate the conductivity of fluid-mixing to gathering _m

B. known water conductivity _w, according to formula ${\mathrm{\α}}_o{\mathrm{\ρ}}_o+{\mathrm{\α}}_w{\mathrm{\ρ}}_w+{\mathrm{\α}}_g{\mathrm{\ρ}}_g=\frac{{2(\mathrm{\ϵ}\·C_1)}^2\mathrm{\Δ\; p}}{V(1-{\mathrm{\β}}^4)},$ ${\mathrm{\σ}}_m=\frac{{2\mathrm{\α}}_l}{3-{\mathrm{\α}}_l}\frac{{2\mathrm{\α}}_w}{3-{\mathrm{\α}}_w}{\mathrm{\σ}}_w$ With formula α _l+ α _g=1 calculates the oil phase cross section contains rate α _o, the water cross section contains rate α _wContain rate α with the gas phase cross section _g, in the formula, σ _mBe profit gas mixed conductivity;

C. according to formula $\begin{array}{c}{Q}_{\mathrm{Vo}}={\mathrm{\φ}}_o{Q}_v\\ Q_{\mathrm{Vw}}={\mathrm{\φ}}_w{Q}_v\\ Q_{\mathrm{Vg}}={\mathrm{\φ}}_g{Q}_v\end{array}$ Calculate the oil phase volume flow Q _Vo, the water volume flow Q _VwWith gas phase volume flow rate Q _Vg

D. according to formula $\begin{array}{c}{V}_o=\frac{Q_{\mathrm{Vo}}}{{\mathrm{\α}}_o{A}_0}\\ V_w=\frac{Q_{\mathrm{Vw}}}{{\mathrm{\α}}_w{A}_0}\\ V_g=\frac{Q_{\mathrm{Vg}}}{{\mathrm{\α}}_g{A}_0}\end{array}$ Calculate oil phase flow velocity V _o, aqueous phase flow rate V _wWith gas phase flow velocity V _g

2. measuring method according to claim 1 is characterized in that, constituting each right electrode of capacitance electrode is 20mm-50mm at the length range of pipeline axial direction, utilizes the relevant speed V of cross correlation measurement method measurement of liquid ₁Constituting each right electrode of resistance electrode is 1mm-10mm at the length range of pipeline axial direction, utilizes the cross correlation measurement method to measure the relevant speed V of gas phase ₂

Images