[go: up one dir, main page]

CN109188016A - Oil-gas-water three-phase flow split-phase flow velocity acoustic-electric bimodal measurement method - Google Patents

Oil-gas-water three-phase flow split-phase flow velocity acoustic-electric bimodal measurement method Download PDF

Info

Publication number
CN109188016A
CN109188016A CN201810990581.6A CN201810990581A CN109188016A CN 109188016 A CN109188016 A CN 109188016A CN 201810990581 A CN201810990581 A CN 201810990581A CN 109188016 A CN109188016 A CN 109188016A
Authority
CN
China
Prior art keywords
water
phase
flow
signal
ultrasonic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810990581.6A
Other languages
Chinese (zh)
Other versions
CN109188016B (en
Inventor
谭超
史雪薇
董峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201810990581.6A priority Critical patent/CN109188016B/en
Publication of CN109188016A publication Critical patent/CN109188016A/en
Application granted granted Critical
Publication of CN109188016B publication Critical patent/CN109188016B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/24Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
    • G01P5/241Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by using reflection of acoustical waves, i.e. Doppler-effect
    • G01P5/242Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by using reflection of acoustical waves, i.e. Doppler-effect involving continuous, e.g. modulated or unmodulated, waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means

Landscapes

  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Electrochemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Measuring Volume Flow (AREA)

Abstract

本发明涉及一种油气水三相流分相流速声电双模态测量方法,包含如下步骤:针对油气水三相水基‑分散流,利用电学传感器采集水相含率信号,利用连续波超声传感器采集多普勒频移信号;根据去噪后的电学传感器的水相含率测量数据,结合特征提取与分类算法实现油气水三相流流型的识别;提取水相含率信号的局部波动特性;利用集合经验模态分解,对连续波超声多普勒频移信号进行分解并选取主分量;采用超声电学联合分析的方法,获得不同流型下水相含率和多普勒频移信号中不同波动尺度的之间的对应关系;对不同流型下、反映分相流速信息的相应本征模态函数进行快速傅里叶变换并计算平均多普勒频移,根据多普勒效应的基本原理,计算分相流速。

The invention relates to a method for measuring the split-phase flow velocity of oil, gas and water three-phase flow, which comprises the following steps: for the three-phase water-based-dispersed flow of oil, gas and water, an electrical sensor is used to collect a water phase holdup signal, and a continuous wave ultrasonic The sensor collects the Doppler frequency shift signal; according to the water phase holdup measurement data of the de-noised electrical sensor, combined with feature extraction and classification algorithms to realize the identification of the three-phase flow pattern of oil, gas and water; extract the local fluctuation of the water phase holdup signal characteristics; using ensemble empirical mode decomposition to decompose the continuous wave ultrasonic Doppler frequency shift signal and select the principal component; using the method of ultrasonic-electrical joint analysis to obtain the water phase holdup and Doppler frequency shift signal under different flow patterns Correspondence between different fluctuation scales; fast Fourier transform is performed on the corresponding eigenmode functions reflecting the information of the split-phase flow velocity under different flow patterns, and the average Doppler frequency shift is calculated. According to the basic Doppler effect The principle is to calculate the split-phase flow rate.

Description

Oil-gas-water three-phase flow split-phase flow velocity acoustic-electric bimodal measurement method
Technical field
The invention belongs to fluid measurement technical field, it is related to a kind of ultrasonic sensor and electric sensor combined type measurement side Method, the non-disturbance formula for three-phase flow split-phase flow velocity measure.
Technical background
Oil-gas-water three-phase flow is three groups of shuntings in gas/liquid/liquid three-phase fluid, it is widely present in the exploitation and fortune of petroleum During defeated.In recent years, oil field enters high water-cut stage, oil-gas-water three-phase flow mostly using water as continuous phase, oil droplet and bubble as from Dephasing flows together with continuous phase.Become the hot issue of multiphase flow area research to the measurement of its flow parameter, especially The measurement of phase content, flow velocity is of great significance for the accurate estimation and the guarantee of production safety of yield.However and single-phase flow It compares, complicated interaction and relative motion characteristic cause its fluidal texture sufficiently complex between three-phase, and flow pattern is complicated and changeable, right The accurate acquisition of flow parameter brings huge challenge.
It is varied for the measurement method of multiphase flow procedure parameter.In traditional measurement method, single-phase instrument is mostly used Table complete process parameter directly acquires, explicit physical meaning.For example than more typical example turbine flowmeter, differential pressure flowmeter, Vortex-shedding meter etc., but these methods belong to intrusive measurement more, can produce to flow process while obtaining flow parameter Raw disturbance, influences fluidal texture and reduces measurement accuracy.Therefore, in scientific research and industrial processes, using non-intruding The method of formula is completed to be of great significance to the accurate acquisition of multiphase flow procedure parameter.Some emerging measurement methods, such as electricity Method, ultrasonic method, microwave method, laser method etc. are gradually applied in multiphase flow measurement.Wherein electric method and ultrasonic method have non-intruding, Sensor structure is simple, easy for installation, low in cost and is concerned.
Electricity phase content measurement method by measurement fluid-mixing electrical parameter indirect gain fluid mutually distribution according to The difference of sensor structure, shape and energisation mode is divided into annular photoconductive array, conducting probe, electricity process tomographic imaging etc. Method.Wherein conductance method is based on annular electro impedance measuring sensor, is a kind of common phase content measurement method, has non-disturb Dynamic, the radiationless, advantages such as response is fast.It is the mixture using conducting medium as continuous phase that it, which requires measurand,.It is by swashing It encourages the electrode driving voltage constant to upper application and establishes electricity sensitivity field, it is sensitive when the phase content of detected fluid changes The impedance operator of field changes therewith, realizes the survey of detected fluid phase content to upper potential difference by obtaining measuring electrode Amount, is widely used in multiphase flow measurement.
Phenomena such as reflection, refraction, scattering occur when method for measuring ultrasonic flow speed is propagated in a fluid by ultrasound convection current The flow velocity of body measures, and has that non-disturbance, sensor structure be simple, the advantages such as low in cost.According to different measuring principles, It is broadly divided into correlation method, time difference method, Doppler method etc..Time difference method is flowed using the spread speed of ultrasonic wave in a fluid and fluid The measurement of the relational implementation flow of speed is chiefly used in the higher situation of discrete phase content.Correlation method using fluid noise signal from One sensor passes to the acquisition for the transition time used when another sensor realizing flow, it is based on cross correlation algorithm, fits For multiple fluid, but the physics practical significance of cross-correlation flow velocity is also and indefinite up to now.Doppler method utilizes ultrasound Doppler effect realize the measurement to the flow velocity of discrete phase, explicit physical meaning.Ultrasonic wave is emitted by ultrasonic transducer Fluid is received after the reflection or scattering of the discrete phase (drop or bubble) moved in heterogeneous fluid by ultrasonic transducer.Root According to Doppler effect, incident acoustic wave and the average flowing speed for receiving frequency difference and discrete phase scattering particles in measured zone between sound wave It spends directly proportional.It mainly include that how general pulsed wave ultrasound Doppler method and continuous wave ultrasound be when ultrasonic Doppler is measured for flow velocity Le method.Pulsed wave ultrasound wave Doppler method is intermittent in the form of a pulse to emit ultrasonic wave into fluid, and between each pulse Every being received in the time to acoustical signal, there is range gating ability, the position of scatterer can be obtained while obtaining flow velocity Confidence breath.But due to how the limitation of Gui this special sample frequency, farthest measure distance and maximum detection flow velocity be limited.Continuous wave Doppler technology continuously emits and receives ultrasonic wave, generally there is a fixed measured zone, can obtain scattering in measured zone The mean flow rate of particle, sensor structure is generally relatively simple, and cost is relatively low, and the limitation without maximum detection flow velocity.Continuous wave When ultrasonic Doppler technique is measured for oil-gas-water three-phase flow flow velocity, discrete phase includes drop and bubble, in flow process Can reflection and scattering process occur to ultrasound.And because the slippage between gas phase and liquid phase be can not ignore, drop is gentle Doppler frequency shift caused by steeping is different.Therefore doppler shifted signal is that have scatterer different in flow rate to super on different location The superposition of acoustical Doppler effect is a kind of sophisticated signal comprising the fluctuation of multi-frequency scale, the fluctuation of some of them dimensions in frequency Reflect out of phase flow rate information.So electric sensor and continuous wave ultrasound doppler sensor are applied in combination, obtain The water phase content and flow rate information of oil-gas-water three-phase flow.By the foundation of Time-Frequency Analysis Method and model, oil gas water three phase is realized The acquisition of flow point phase flow velocity.
Summary of the invention
The purpose of the present invention is utilizing multi-sensor information, a kind of accurate, non-turbulent oil gas water three phase flow point is provided Phase flow velocity extracting method.Sensor of the invention and measurement method are to combine ultrasonic Doppler sensor and electric sensor, It is respectively used to obtain the flow velocity and water phase content information of oil-gas-water three-phase flow.Using signal processing technology to multi-source information at Reason: to water phase content data carry out feature extraction and combining classification algorithm can realize oil-gas-water three-phase flow flow characteristics analysis and The correct identification of flow pattern is carried out set empirical mode decomposition to ultrasonic doppler measurements data on this basis and is joined using acoustic-electric The method for closing analysis can realize that the unperturbed dynamic formula of the split-phase flow velocity of oil gas water under different flow patterns accurately obtains.Technical side of the invention Case is as follows:
A kind of oil-gas-water three-phase flow split-phase flow velocity acoustic-electric bimodal measurement method, using a twin crystal ultrasonic transducer and electricity Learn sensor;The twin crystal ultrasonic transducer is for obtaining ultrasonic Doppler frequency shift signal, using the structure that transceiver is ipsilateral, Inside is respectively used to receive and emit ultrasonic wave, the twin crystal comprising receiving piezoelectric ceramic wafer and transmitting piezoelectric ceramic chip Ultrasonic transducer is installed in horizontal pipe bottom and guarantees that with the angle of horizontal direction be θ;The electric sensor is for obtaining Water intaking phase content signal, is installed in horizontal pipe simultaneously with ultrasonic transducer;Measurement method comprises the following steps:
1) for oil gas water three phase it is water base-dispersion stream, using electric sensor acquire water phase content signal H (t), t ∈ (0, T), doppler shifted signal f is acquired using continuous wave ultrasound sensord(t), (0, T) t ∈, wherein T is the sampling time, and is gone Except the high-frequency noise in two signals;
2) according to the water phase content measurement data of the electric sensor after denoising, binding characteristic is extracted to be realized with sorting algorithm The identification of oil-gas-water three-phase flow flow pattern;
3) it extracts the localised waving characteristic of water phase content signal: the water phase content signal under different flow patterns being carried out continuous small Wave conversion is obtained in the continuous wave wavelet conversion coefficient to dimensioning a and given time bWith the local wavelet energy coefficient of calculatingWherein ψ (t) is Mother wavelet function, using Daubechies4 small echo, out to out 500;When m- scale axis plot part wavelet energy system Number spectrogram obtains the localised waving characteristic of water phase content signal;
4) using set empirical mode decomposition, continuous wave ultrasound doppler shifted signal is decomposed and chooses main point Amount: using the method for set empirical mode decomposition, by DOPPLER ULTRASOUND SIGNAL fd(t), t ∈ (0, T) is decomposed into several step by step Levy mode function IMFj(t) (j=1,2 ..., n) and trend function rn(t) sum:Wherein n To decompose obtained intrinsic mode function number, preceding 3 grades of intrinsic mode function IMF are choseni(t) (i=1,2,3) is mentioned as flow velocity The principal component taken;
5) according to step 3), 4) in decomposition result obtained under different flow patterns using the method for ultrasonic electric Conjoint Analysis Corresponding relationship in water phase content and doppler shifted signal between different Fluctuation Scales, and the principal component obtained in step 4) The middle extraction chosen different component and be used for split-phase flow velocity, for water base-dispersion-bubble flow and water base-dispersion-slug flow, with Large scale, i.e. scale > 50 in the local wavelet energy coefficient spectrogram of water phase content signal, wave component there is the corresponding relationship to be The first order intrinsic mode function IMF that doppler shifted signal decomposes1(t), it is high fdrequency component and has in time Discontinuity reflects the velocity information for being interrupted the discrete phase bubble of appearance in flow process;And it is small with the part of water phase content signal It is that doppler shifted signal decomposes to obtain that the wave component of wave energy coefficient spectrogram Small and Medium Sized (scale < 50), which has corresponding relationship, The second level and third level intrinsic mode function IMF2(t)&IMF3(t), they are low frequency component and have in time continuous Property, reflect the velocity information of fully decentralized oil droplet in flow process;For water base-dispersion-plug-like stream, believe with water phase content Number local wavelet energy coefficient spectrogram in large scale, i.e. scale > 50, wave component to have corresponding relationship be Doppler frequency shift The first order and second level intrinsic mode function IMF that signal decomposition obtains1(t)&IMF2(t), they are high fdrequency components and in the time It is upper that there is discontinuity, reflect the velocity information for being interrupted the discrete phase vent plug of appearance in flow process;And with water phase content signal Local wavelet energy coefficient spectrogram Small and Medium Sized, i.e. scale < 50, wave component to have corresponding relationship be Doppler frequency shift letter Number decompose obtained third level intrinsic mode function IMF3(t), it is low frequency component and has continuity in time, reflects The velocity information of fully decentralized oil droplet in flow process.
6) under different flow patterns, reflection split-phase flow rate information corresponding intrinsic mode function carry out Fast Fourier Transform (FFT) simultaneously Average Doppler frequency shift is calculated, according to the basic principle of Doppler effect, calculates split-phase flow velocity: for water base-dispersion-bubble flow With water base-dispersion-slug flow, the average true flow velocity of gas phase are as follows:The average true flow velocity of drop isFor water base-dispersion-plug-like stream, the average true flow velocity of gas phase isDrop Average true flow velocity are as follows:Wherein f0For the driving frequency of ultrasonic wave transmitting probe, θ is ultrasonic wave acoustic beam direction With the angle of horizontal direction, c is the velocity of sound in fluid,To utilize the fast of different intrinsic mode functions The average Doppler frequency shift that fast Fourier transformation obtains, Sd*(fd*) it is corresponding to frequency displacement fd*Doppler spectral intensity, * is step 5) reflect the corresponding principal component of split-phase flow rate information under the different flow patterns determined in.
Substantive distinguishing features of the invention are: oil gas water three phase is water base-dispersion stream in, discrete phase includes oil droplet and bubble, and It will form different flow patterns according to the relationship between different gas-liquids.Under different flow patterns, different discrete phases are due to flowing over There is the difference of distribution and size, the wave that water phase content can be caused different in journey.Mutually containing based on electricity sensitivity principle Rate information detection technology can accurately obtain the volatility series of water phase content when the phase content of detected fluid is distributed and changes.Base In ultrasonic Doppler effect, can occur in oil droplet and bubble surface when continuous wave ultrasound signal is propagated in oil-gas-water three-phase flow Reflection and scattering phenomenon, and because the slippage between gas phase and liquid phase be can not ignore, so drop it is gentle steep oneself-meeting cause it is different more General Le frequency displacement.Therefore doppler shifted signal is a kind of sophisticated signal comprising the fluctuation of multi-frequency scale, some of them frequency ruler The fluctuation of degree reflects out of phase flow rate information.Therefore, by electricity phase content measurement sensor and continuous wave ultrasound Doppler Sensor combines-dispersion stream water base to oil gas water three phase and measures: the fluctuation of water phase content is obtained using electric sensor Time series obtains local wavelet energy coefficient spectrogram using continuous wavelet transform technology, obtains the part of water phase content signal Wave characteristic, and binding characteristic is extracted and the identification of classification and identification algorithm realization three-phase flow flow pattern;It is how general using continuous wave ultrasound It strangles sensor and obtains doppler shifted signal, and the signal is decomposed based on set empirical mode decomposition and is able to reflect The principal component of split-phase flow velocity;Electricity and ultrasonic signal are subjected to Conjoint Analysis, determining under different flow patterns can reflect in principal component The component of different discrete phase flowing informations, and the average how general of respective component is calculated by the Fast Fourier Transform (FFT) to principal component Strangle frequency displacement, the final acquisition for realizing split-phase flow velocity.Beneficial effects of the present invention and advantage are as follows:
1. this method is non-disturbance formula measurement means, any disturbance will not be generated to fluid;
2. this method can be obtained using the measurement method of electricity, ultrasonic double-mode state to flow process a more complete description;
3. this method uses ultrasonic electric Conjoint Analysis, can accurately be extracted from continuous wave ultrasound doppler shifted signal The split-phase flow velocity of oil-gas-water three-phase flow, it is simple and easy;
4. this method measurement is convenient, speed is fast, at low cost.
Detailed description of the invention
The following drawings describes the selected embodiment of the present invention, is exemplary diagram and non exhaustive or restricted, in which:
Ultrasound and electricity bimodal measuring system structural schematic diagram in measurement method Fig. 1 of the invention;
Ultrasonic Doppler sensor structure schematic diagram in measurement method Fig. 2 of the invention;
Conductance ring sensor structure schematic diagram in measurement method Fig. 3 of the invention;
Gather empirical mode decomposition flow chart in measurement method Fig. 4 of the invention;
The calculating step of oil-gas-water three-phase flow split-phase flow velocity in measurement method Fig. 5 of the invention;
Fig. 6 liquid phase apparent velocity verification result;
Fig. 7 gas phase apparent velocity verification result;
The total apparent velocity verification result of Fig. 8.
Specific embodiment
Below with reference to the instruction sheet embodiment that the present invention will be described in detail.
Fig. 1 is ultrasound and electricity bimodal measuring system structural schematic diagram in measurement method of the invention.The electricity passes Sensor is illustrated by taking conductance ring sensor as an example.The measuring system is by being installed on the one group of electric sensor tested on pipeline section 1 Array 2 and coupled electrical signal occur and acquisition unit 4, one group of ultrasonic sensor array 3 and coupled Ultrasonic signal occurs to constitute with acquisition unit 6.Electric sensor array 2 is embedded in the ring of inner wall of the pipe by four at certain intervals Shape metal electrode composition;Ultrasonic sensor array 3 is installed on the twin crystal ultrasonic wave transducer of the ipsilateral transceiver of pipeline by one group Device composition.(installation is suitable for installation at certain intervals on tested pipeline 1 with ultrasonic sensor array 3 for the electric sensor array 2 Sequence is unlimited) one dual-modality sensor array of composition, it works at the same time.When measured stream multiphase flow enters measured tube from direction of flow 0 When road, electric sensor array 2 is fluctuated by the water phase content that reflection detected fluid occurs to obtain with acquisition unit 4 for electrical signal The measurement data of information, by the data be sent into phase content data processing unit 5 carry out data processing complete analysis of flow characteristics with Flow Regime Ecognition;The flow velocity that reflection detected fluid is obtained with acquisition unit 6 occurs by ultrasonic signal for ultrasonic sensor array 3 simultaneously It is complete to be sent into the progress data processing of doppler shift data processing unit 7 by the measurement data of the Doppler frequency shift of information for the data At the decomposition of Fluctuation Scales different in data;Analysis of flow characteristics described above and Flow Regime Ecognition result and frequency-shift data are decomposed As a result it is sent into the acquisition that split-phase flow relocity calculation unit 8 completes the split-phase flow velocity of oil-gas-water three-phase flow.
Fig. 2 is the structural schematic diagram of ultrasonic Doppler sensor in measurement method of the invention.Surpass used in the invention patent Sound doppler transducer is the twin crystal ultrasonic transducer of ipsilateral transceiver, and internal includes a transmitting piezoelectric ceramic chip 3a With a reception piezoelectric ceramic wafer 3b, it is respectively attached on acoustic-electric coupling material 3c and 3d.Acoustic-electric coupling material 3c and 3d is straight It connects and is contacted with detected fluid 2, and with angle theta between test pipeline 10Installation makes ultrasonic beam access and detected fluid 2 It flows direction 0 and keeps angle theta.It joined sound insulating material 4 between acoustic-electric the coupling material 3c and 3d to prevent transmitting and reception sound Interference is generated between wave.The twin crystal ultrasonic transducer is installed on the bottom of test pipeline 1.Transmitting piezoelectric ceramic chip 3a hair Penetrate ultrasonic wave, sound wave is propagated in fluid 2, receive measurement space 5 in discrete phase (bubble 6, drop 7) reflection and scattering after by Piezoelectric ceramic wafer 3b is received to receive.It receives sound wave and emits the flow velocity of fluid in the frequency difference reflection measurement space 5 between sound wave Information.
Fig. 3 is the structural schematic diagram of conductance ring sensor in measurement method of the invention.The electric sensor is by four Circular metal electrode (1,2,3,4) is axially aligned and is embedded in inner wall of the pipe 0 and forms, and wherein electrode 1 and electrode 4 are excitation For electrode to for generating stable electricity sensitivity field therebetween, electrode 2 and electrode 3 are measuring electrode to for measuring electricity therebetween Potential difference.Using the working method of current excitation voltage measurement, by injecting constant square-wave excitation current into electrode 1 and will be electric Pole 4 is grounded, and stable electricity sensitivity field can be formed among electrode 1 and 4.The medium phase content when detected fluid flows through sensitivity field And distribution changes the potential difference between electrode 2 and electrode 3 can be caused accordingly to change, it can be real by measuring the potential difference The measurement of existing fluid water phase content.
Fig. 4 is set empirical mode decomposition flow chart used in measurement method of the invention.Gathering empirical mode decomposition is It is proposed to solve the problems, such as the modal overlap in Conventional wisdom mode decomposition, main thought is more to some classification in statistics The accuracy of measurement can be improved in the secondary average value sought that measures.Set empirical mode decomposition passes through in signal x (t) to be decomposed The different white noise ω of m group is repeatedly addedi(t), (i=1,2 ..., m) the different composite signal X of composition m groupi(t), (i=1, 2 ..., m), then empirical mode decomposition is carried out to every group of composite signal.Empirical mode decomposition is by the process of screening by the above institute The composite signal X statedi(t), (i=1,2 ..., m) present in the fluctuation of different scale separate step by step, generating a series of has The intrinsic mode functions IMF of different scaleij, (i=1,2 ..., m;J=1,2 ..., n), until remaining value sequence rinIt (t) is single Letter of transfer number stops decomposing, and wherein n is to decompose the intrinsic mode function number generated every time.Finally the whole gamut of experience mode decomposition is obtained The intrinsic mode function component of each group arrived is averaging the intrinsic mode function as set empirical mode decomposition
Fig. 5 is the calculating steps flow chart that present invention ultrasound obtains oil-gas-water three-phase flow split-phase flow velocity with electricity combination sensor Figure.For oil gas water three phase it is water base-dispersion stream, electricity of the invention, ultrasound data conjoint analysis method obtain split-phase flow velocity Steps are as follows for calculating:
Step 1: acquiring the water phase content signal H (t) of electric sensor, Doppler's frequency of t ∈ (0, T) and ultrasonic sensor Shifting signal fd(t), (0, T) t ∈, and the high-frequency noise (denoising method is unlimited) in two signals is removed, when wherein T is sampling Between;
Step 2: according to the electric sensor water phase content measurement data after denoising, binding characteristic is extracted real with sorting algorithm The identification of existing oil-gas-water three-phase flow flow pattern (water base-dispersion-bubble flow, water base-dispersion-plug-like stream, water base-dispersion-slug flow);
The general step of Flow Regime Ecognition are as follows: be able to reflect the spy of variations in flow patterns feature by extracting from measurement data first Value indicative, such as statistical method, Time-Frequency Analysis Method, nonlinear analysis method;Then support vector machines, artificial neural network are utilized Etc. sorting algorithms characteristic value is subjected to classification identification.By taking Continuous Wavelet Transform as an example, first to the water phase content after denoising Measurement data H (t) carries out continuous wavelet transform (morther wavelet is Daubechies4 small echo, out to out 500):
Wherein wa,bFor in the Continuous Wavelet Transform Coefficients for giving dimensioning a and given time b, ψ (t) is mother wavelet function.By Water phase content wave characteristic under different flow patterns is different, therefore m- small echo-energy spectrum parameter when calculating:
The parameter is able to reflect water phase content sequence along the energy variation of time shaft.The mean value and variance for calculating E (b) can Feature as Flow Regime Ecognition inputs trained artificial neural network, and achievable oil gas water three phase is water base-dispersion stream exemplary flow The identification of type.
Step 3: water phase content sequence being handled using Time-frequency Decomposition and the method for feature extraction, obtains different flow patterns Under water phase content signal in the localised waving characteristic as caused by different discrete phases;
Oil gas water three phase is water base-dispersion stream in, discrete phase includes oil droplet and bubble.Since different discrete phases are in flow process In discreteness and size difference: the size of discrete phase oil droplet is less than the size of discrete phase bubble, and oil droplet is thoroughly dispersed in It is flowed together in water with water phase, and bubble interruption occurs, this wave that water phase content can be caused different.Using when frequency division The method of solution decomposes the wave characteristic of the water phase content measurement data of the electric sensor after denoising, in combination with feature The method of parameter extraction can determine the different wave characteristics of the water phase content as caused by different discrete phases, to disclose flow process The flow behavior of middle difference discrete phase.By taking continuous wavelet transform as an example, water is calculated on the basis of continuous wavelet transform in step 2 Phase content signal is in the local wavelet energy coefficient for giving dimensioning a and given time b:
Local wavelet energy coefficient can be used for reflecting signal to dimensioning and when the fluctuating range inscribed, when m- ruler The localised waving characteristic of the available water phase content signal of degree axis plot part wavelet energy coefficient spectrogram: large scale (scale > 50) wave component mainly reflects the low frequency fluctuation of water phase content caused by discrete phase bubble, small scale (scale < 50) Wave component mainly reflect the high frequency fuctuation within a narrow range of water phase content caused by discrete phase oil droplet.
Step 4: using set empirical mode decomposition, continuous wave ultrasound doppler shifted signal being decomposed;
Based on the method for set empirical mode decomposition, by the process of screening by DOPPLER ULTRASOUND SIGNAL fd(t),t∈(0, T it) is decomposed into several intrinsic mode function IMF step by stepj(t) (j=1,2 ..., n) and trend function rn(t) sum:
Wherein n is to decompose obtained intrinsic mode function number.Under normal conditions, the component obtained by screening process In, it is most important component that the big high fdrequency component of energy, which has usually contained most significant, most important information in signal,.Therefore it selects Take preceding 3 grades of intrinsic mode function IMFi(t) principal component that (i=1,2,3) is extracted as flow velocity;
Step 5: according to Flow Regime Ecognition result and decomposition result, using ultrasonic electric Conjoint Analysis, dividing under different flow patterns Calculating of the different principal components for split-phase flow velocity is not chosen.
1) method for using ultrasonic electric Conjoint Analysis, obtains different flow patterns and is lauched in phase content and doppler shifted signal Corresponding relationship between different Fluctuation Scales, and different components is chosen for split-phase stream in the principal component that step 4 obtains The extraction of speed.Local wavelet energy for water base-dispersion-bubble flow and water base-dispersion-slug flow, with water phase content signal It is the first order eigen mode that doppler shifted signal decomposes that the wave component of large scale, which has corresponding relationship, in coefficient spectrogram State function IMF1(t), it is high fdrequency component and has discontinuity in time, reflects and is interrupted the discrete of appearance in flow process The velocity information of phase bubble;And have pair with the wave component of the local wavelet energy coefficient spectrogram Small and Medium Sized of water phase content signal Should be related to be the second level and third level intrinsic mode function IMF that doppler shifted signal decomposes2(t)&IMF3(t), it Be low frequency component and in time have continuity, reflect the velocity information of fully decentralized oil droplet in flow process.It is right The wave component of large scale has pair in water base-dispersion-plug-like stream, with the local wavelet energy coefficient spectrogram of water phase content signal Should be related to be the first order and second level intrinsic mode function IMF that doppler shifted signal decomposes1(t)&IMF2(t), it Be high fdrequency component and in time have discontinuity, reflect be interrupted in flow process the discrete phase vent plug of appearance speed letter Breath;And how general have corresponding relationship with the wave component of the local wavelet energy coefficient spectrogram Small and Medium Sized of water phase content signal is Strangle the third level intrinsic mode function IMF that frequency shift signal decomposes3(t), it is low frequency component and has in time continuous Property, reflect the velocity information of fully decentralized oil droplet in flow process.
2) under different flow patterns, reflection split-phase flow rate information corresponding intrinsic mode function carry out Fast Fourier Transform (FFT) simultaneously It calculates average Doppler frequency shift: the Doppler in Fast Fourier Transform (FFT) and sampling time is carried out to first order intrinsic mode function Average frequency shift:The adduction of the first order and second level intrinsic mode function component is carried out quick Fourier transformation simultaneously calculates Doppler's average frequency shift in the sampling time:To the second level Fast Fourier Transform (FFT) is carried out with the adduction of third level intrinsic mode function component and calculates the Doppler in liquid film time range Average frequency shiftFast Fourier Transform (FFT) is carried out to third level intrinsic mode function and is sampled Doppler's average frequency shift in time:Wherein fd1, fd2, fd3, fd4Respectively IMF1(t), IMF1(t)+IMF2(t), IMF2(t)+IMF3(t), IMF3(t) corresponding Doppler frequency shift, Sd1(fd1), Sd2(fd2), Sd3(fd3), Sd4(fd4) be respectively its corresponding Doppler frequency shift energy spectral intensity.
3) according to the basic principle of Doppler effect, split-phase flow velocity is calculated: for water base-dispersion-bubble flow and water base-point Dissipate-slug flow, the average true flow velocity of gas phase are as follows:The average true flow velocity of drop isFor water base-dispersion-plug-like stream, the average true flow velocity of gas phase isDrop Average true flow velocity are as follows:Wherein f0For the driving frequency of ultrasonic wave transmitting probe, θ is ultrasonic wave acoustic beam direction With the angle of horizontal direction, c is the velocity of sound in fluid,Be the 2 of step 5) in using different intrinsic mode functions quick Fu in The average Doppler frequency shift that leaf transformation obtains, * be step 5 1) in determine different flow patterns under reflect split-phase flow rate information phase The principal component answered.
Experimental verification is as follows:
Dynamic experiment is completed on oil-gas-water three-phase flow experimental provision.The fluid-transporting tubing of experimental provision is internal diameter 50mm Stainless steel tube, multiphase inflow entrance to outlet distance be about 16.6m.Electricity and ultrasonic combination sensor are installed on apart from multiphase flow The test pipeline section of entrance about 13m, so that flow pattern sufficiently develops.Meanwhile be equipped on test pipeline section pressure gauge and thermometer so as to Working condition is recorded, video camera is for observing and recording flow pattern.In experiment fluid used be respectively tap water (density 998kg/m3, 1.01 × 10-3Pas of dynamic viscosity), 15# technical white oil (density 790kg/m3,3.9 × 10-2Pas of dynamic viscosity) and dry Dry air (density 1.2kg/m3,1.81 × 10-5Pas of dynamic viscosity).It is each mutually real in inlet by a T-type mixer Now mixing is simultaneously measured by the single-phase flowmeter of standard to each phase flow rate is carried out before mixing.By changing out of phase stream Amount proportion, generates different nowed formings.Every group of experiment fixing oil phase and water phase flow gradually increase gas phase flow rate, are formed not Same flow pattern (water base-dispersion-bubble flow, water base-dispersion-plug-like stream, water base-dispersion-slug flow).In experiment, mark three under condition Mutually stream total flow range is in 6.66m3/ h to 14.97m3/ h, total apparent change in flow range is in 0.92m/s to 1.76m/s.It is each The time of test constantly is 10s under flow regime, completes the wave time sequence of water phase content and obtaining for doppler shifted signal It takes.
The wave time sequence and Doppler frequency shift to water phase content are completed according to the step in the specific embodiment The processing of signal, can obtain oil gas water three phase it is water base-split-phase true mean flow rate of the dispersion stream under different flow patterns.For verifying institute State the validity of method, it is assumed that slippage is not present between oil and water, the Real-time Water phase content measured using conductance ring can obtain meter Calculate real-time gas phase content and in real time total liquid phase content.
Wherein φwFor the Real-time Water phase content in liquid phase;βl, βg, βwFor real-time total liquid phase content, real-time gas phase content and Real-time Water phase content.βwIt is obtained by electricity phase content sensor, since slippage φ being not present between oil and waterwEntrance can be passed through It is obtained containing rate.Further calculate split-phase apparent velocity:
Wherein Jl, Jg, J is respectively liquid phase apparent velocity, gas phase apparent velocity and total apparent velocity.It is referred to entrance Apparent velocity compares verifying.As shown in figs 6-8, using root-mean-square error as evaluation index: liquid phase is apparent, gas phase for verification result Apparently, the root-mean-square error of summary table pipe is respectively 0.06m/s, 0.06m/s, 0.07m/s.Therefore, the method can be obtained accurately Take the split-phase flow velocity of oil-gas-water three-phase flow.

Claims (1)

1.一种油气水三相流分相流速声电双模态测量方法,采用一个双晶超声换能器和电学传感器;所述双晶超声换能器用于获取超声多普勒频移信号,采用收发一体同侧的结构,内部包含接收压电陶瓷晶片和发射压电陶瓷晶片,分别用于接收和发射超声波,所述双晶超声换能器被安装于水平管道底部并保证与水平方向的夹角为θ;所述电学传感器用于获取水相含率信号,与超声换能器同时安装于水平管道中;测量方法包含如下步骤:1. An acoustic-electric dual-modal measurement method for the split-phase flow velocity of oil-gas-water three-phase flow, using a twin-crystal ultrasonic transducer and an electrical sensor; the twin-crystal ultrasonic transducer is used to obtain an ultrasonic Doppler frequency shift signal, It adopts the structure of receiving and transmitting on the same side, and contains a receiving piezoelectric ceramic chip and a transmitting piezoelectric ceramic chip, which are used to receive and transmit ultrasonic waves respectively. The dual-crystal ultrasonic transducer is installed at the bottom of the horizontal pipe and ensures the The included angle is θ; the electrical sensor is used to obtain the water holdup signal, and is installed in the horizontal pipeline simultaneously with the ultrasonic transducer; the measurement method includes the following steps: 1)针对油气水三相水基-分散流,利用电学传感器采集水相含率信号H(t),t∈(0,T),利用连续波超声传感器采集多普勒频移信号fd(t),t∈(0,T),其中T为采样时间,并去除两个信号中的高频噪声;1) For the oil-gas-water three-phase water-based-dispersed flow, the water phase holdup signal H(t),t∈(0,T) is collected by the electrical sensor, and the Doppler frequency shift signal f d ( t), t∈(0,T), where T is the sampling time, and removes high-frequency noise in the two signals; 2)根据去噪后的电学传感器的水相含率测量数据,结合特征提取与分类算法实现油气水三相流流型的识别;2) According to the water phase holdup measurement data of the electrical sensor after denoising, combined with feature extraction and classification algorithms to realize the identification of the three-phase flow pattern of oil, gas and water; 3)提取水相含率信号的局部波动特性:对不同流型下的水相含率信号进行连续小波变换获得在给定尺度a和给定时刻b的连续波小波变换系数和计算局部小波能量系数其中ψ(t)为母小波函数,采用Daubechies4小波,总尺度为500;在时间-尺度轴上绘制局部小波能量系数谱图获得水相含率信号的局部波动特性;3) Extract the local fluctuation characteristics of the water phase holdup signal: perform continuous wavelet transform on the water phase holdup signal under different flow patterns to obtain the continuous wave wavelet transform coefficients at a given scale a and a given time b and calculate the local wavelet energy coefficients where ψ(t) is the mother wavelet function, Daubechies4 wavelet is used, and the total scale is 500; the local wavelet energy coefficient spectrum is drawn on the time-scale axis to obtain the local fluctuation characteristics of the water holdup signal; 4)利用集合经验模态分解,对连续波超声多普勒频移信号进行分解并选取主分量:利用集合经验模态分解的方法,将超声多普勒信号fd(t),t∈(0,T)逐级分解为若干本征模态函数IMFj(t)(j=1,2,...,n)和趋势函数rn(t)的和:其中n为分解得到的本征模态函数个数,选取前3级本征模态函数IMFi(t)(i=1,2,3)作为流速提取的主分量;4) Using the ensemble empirical mode decomposition, decompose the continuous wave ultrasonic Doppler frequency shift signal and select the principal component: using the method of ensemble empirical mode decomposition, the ultrasonic Doppler signal f d (t), t∈( 0,T) is decomposed into the sum of several eigenmode functions IMF j (t) (j=1,2,...,n) and trend function r n (t): where n is the number of eigenmode functions obtained by decomposition, and the first three eigenmode functions IMF i (t) (i=1, 2, 3) are selected as the principal components extracted from the flow velocity; 5)根据步骤3)、4)中的分解结果,采用超声电学联合分析的方法,获得不同流型下水相含率和多普勒频移信号中不同波动尺度的之间的对应关系,并在步骤4)得到的主分量中选取不同的分量用于分相流速的提取,对于水基-分散-泡状流和水基-分散-弹状流,与水相含率信号的局部小波能量系数谱图中大尺度,即尺度&gt;50,的波动分量有对应关系的为多普勒频移信号分解得到的第一级本征模态函数IMF1(t),它是高频分量且在时间上具有间断性,反映流动过程中间断出现的离散相气泡的速度信息;而与水相含率信号的局部小波能量系数谱图中小尺度(尺度&lt;50)的波动分量有对应关系的为多普勒频移信号分解得到的第二级和第三级本征模态函数IMF2(t)&IMF3(t),它们是低频分量且在时间上具有连续性,反映了流动过程中完全分散的油滴的速度信息;对于水基-分散-塞状流,与水相含率信号的局部小波能量系数谱图中大尺度,即尺度&gt;50,的波动分量有对应关系的为多普勒频移信号分解得到的第一级和第二级本征模态函数IMF1(t)&IMF2(t),它们是高频分量且在时间上具有间断性,反映了流动过程中间断出现的离散相气塞的速度信息;而与水相含率信号的局部小波能量系数谱图中小尺度,即尺度&lt;50,的波动分量有对应关系的为多普勒频移信号分解得到的第三级本征模态函数IMF3(t),它是低频分量且在时间上具有连续性,反映了流动过程中完全分散的油滴的速度信息。5) According to the decomposition results in steps 3) and 4), the method of ultrasonic-electrical joint analysis is used to obtain the water phase holdup under different flow patterns and the corresponding relationship between different fluctuation scales in the Doppler frequency shift signal, and in Step 4) Select different components from the principal components obtained for the extraction of phase-splitting flow velocity. For water-based-dispersion-bubble flow and water-based-dispersion-slug flow, the local wavelet energy coefficient of the water-phase holdup signal is The fluctuation component of the large scale in the spectrum, that is, the scale &gt; 50, has a corresponding relationship with the first-order eigenmode function IMF 1 (t) obtained by the decomposition of the Doppler frequency-shifted signal, which is a high-frequency component and is in the It is discontinuous in time, reflecting the velocity information of discrete phase bubbles that appear intermittently in the flow process; and the fluctuation component of the small scale (scale &lt; 50) in the local wavelet energy coefficient spectrum of the water phase holdup signal has a corresponding relationship as The second-order and third-order eigenmode functions IMF 2 (t) & IMF 3 (t) obtained from the decomposition of the Doppler frequency-shifted signal are low-frequency components and are continuous in time, reflecting the complete flow in the flow process. Velocity information of dispersed oil droplets; for water-based-dispersed-plug flow, there is a corresponding relationship with the large-scale wavelet component in the local wavelet energy coefficient spectrum of the water-phase holdup signal, that is, the scale &gt; 50. The first-order and second-order eigenmode functions IMF 1 (t) & IMF 2 (t) obtained from the decomposition of the Peller frequency-shifted signal, which are high-frequency components and are discontinuous in time, reflecting the discontinuity in the flow process The velocity information of the discrete phase gas plug that appears; and the fluctuation component corresponding to the small scale in the local wavelet energy coefficient spectrum of the water phase holdup signal, that is, the scale &lt; 50, is obtained by the decomposition of the Doppler frequency-shifted signal. The third-order eigenmode function IMF 3 (t), which is a low-frequency component and is continuous in time, reflects the velocity information of oil droplets that are completely dispersed during the flow. 6)对不同流型下、反映分相流速信息的相应本征模态函数进行快速傅里叶变换并计算平均多普勒频移,根据多普勒效应的基本原理,计算分相流速:对于水基-分散-泡状流和水基-分散-弹状流,气相的平均真实流速为:液滴的平均真实流速为对于水基-分散-塞状流,气相的平均真实流速为液滴的平均真实流速为:其中f0为超声波发射探头的激励频率,θ为超声波声束方向与水平方向的夹角,c为流体中声速,为利用不同本征模态函数的快速傅里叶变换得到的平均多普勒频移,Sd*(fd*)为对应于频移fd*的多普勒谱强度,*为步骤5)中确定的不同流型下反映分相流速信息的相应的主分量。6) Perform fast Fourier transform on the corresponding eigenmode functions reflecting the information of the split-phase flow velocity under different flow patterns and calculate the average Doppler frequency shift. According to the basic principle of the Doppler effect, calculate the split-phase flow velocity: for For water-based-dispersed-bubble flow and water-based-dispersed-slug flow, the average true flow velocity of the gas phase is: The average true flow rate of the droplets is For water-based-dispersion-plug flow, the average true flow rate of the gas phase is The average true flow rate of the droplets is: where f 0 is the excitation frequency of the ultrasonic transmitting probe, θ is the angle between the direction of the ultrasonic sound beam and the horizontal direction, c is the speed of sound in the fluid, is the average Doppler frequency shift obtained by using the fast Fourier transform of different eigenmode functions, S d* (f d* ) is the Doppler spectral intensity corresponding to the frequency shift f d* , * is step 5 The corresponding principal components that reflect the information of the split-phase flow velocity under different flow patterns determined in ).
CN201810990581.6A 2018-08-28 2018-08-28 Acoustic-electric bimodal measurement method for phase-splitting flow velocity of oil-gas-water three-phase flow Active CN109188016B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810990581.6A CN109188016B (en) 2018-08-28 2018-08-28 Acoustic-electric bimodal measurement method for phase-splitting flow velocity of oil-gas-water three-phase flow

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810990581.6A CN109188016B (en) 2018-08-28 2018-08-28 Acoustic-electric bimodal measurement method for phase-splitting flow velocity of oil-gas-water three-phase flow

Publications (2)

Publication Number Publication Date
CN109188016A true CN109188016A (en) 2019-01-11
CN109188016B CN109188016B (en) 2020-05-05

Family

ID=64916677

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810990581.6A Active CN109188016B (en) 2018-08-28 2018-08-28 Acoustic-electric bimodal measurement method for phase-splitting flow velocity of oil-gas-water three-phase flow

Country Status (1)

Country Link
CN (1) CN109188016B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114563470A (en) * 2022-01-28 2022-05-31 天津大学 Oil-water two-phase flow phase content measuring method based on nonlinear ultrasound
CN115106499A (en) * 2022-06-30 2022-09-27 北京科技大学 A method and system for judging abnormal fluctuation of crystallizer liquid level
DE102022105759A1 (en) 2022-03-11 2023-09-14 Krohne Ag Flow measuring device and method for operating the same
CN118655184A (en) * 2024-08-21 2024-09-17 南京尤力浦能源科技有限公司 A non-invasive multiphase flow content detection method and device based on continuous wave

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS58100704A (en) * 1981-12-11 1983-06-15 Nippon Kagaku Kogyo Kk Measuring apparatus of gas-liquid two phase flow
JPH02112757A (en) * 1988-10-21 1990-04-25 Babcock Hitachi Kk Instrument for measuring concentration of particulate matter in piping
EP0691527A1 (en) * 1994-07-05 1996-01-10 Institut Français du Pétrole Apparatus for method for measuring the velocity profile in a polyphase fluid
WO2001067050A1 (en) * 2000-03-09 2001-09-13 Nest International N.V. Simultaneous determination of multiphase flowrates and concentrations
CN101280680A (en) * 2008-05-23 2008-10-08 安东石油技术(集团)有限公司 Three-phase flow measuring device
CN103175986A (en) * 2013-03-08 2013-06-26 燕山大学 Method for measuring gas-phase flow rate and liquid-phase flow rates of oil-gas-water three-phase flow
CN104090020A (en) * 2014-07-10 2014-10-08 天津大学 Electric and ultrasonic-based bimodal multiphase flow measuring device
CN104089985A (en) * 2014-07-10 2014-10-08 天津大学 Visual multiphase flow test method based on electric and ultrasonic sensing principle
CN104155471A (en) * 2014-07-10 2014-11-19 天津大学 Multiphase flow testing method based on cross-correlation velocity measurement of ultrasonic and electric multiple sensors
CN105181996A (en) * 2015-08-12 2015-12-23 天津大学 Two-phase flow velocity acoustoelectric bimodal measuring method
CN106932606A (en) * 2017-02-27 2017-07-07 天津大学 Gas-liquid two-phase slug flow and slug flow flow velocity acoustic-electric bimodal measuring method
CN106996988A (en) * 2017-03-20 2017-08-01 天津大学 Oil gas water three phase plug-like disperses flow velocity measuring method
CN107024603A (en) * 2017-02-27 2017-08-08 天津大学 Gas-liquid two-phase bubbly flow flow velocity acoustic-electric bimodal measuring method
CN107153086A (en) * 2017-04-28 2017-09-12 天津大学 The acoustic-electric bimodal fusion measuring method of two phase flow procedure parameter

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS58100704A (en) * 1981-12-11 1983-06-15 Nippon Kagaku Kogyo Kk Measuring apparatus of gas-liquid two phase flow
JPH02112757A (en) * 1988-10-21 1990-04-25 Babcock Hitachi Kk Instrument for measuring concentration of particulate matter in piping
EP0691527A1 (en) * 1994-07-05 1996-01-10 Institut Français du Pétrole Apparatus for method for measuring the velocity profile in a polyphase fluid
WO2001067050A1 (en) * 2000-03-09 2001-09-13 Nest International N.V. Simultaneous determination of multiphase flowrates and concentrations
CN101280680A (en) * 2008-05-23 2008-10-08 安东石油技术(集团)有限公司 Three-phase flow measuring device
CN103175986A (en) * 2013-03-08 2013-06-26 燕山大学 Method for measuring gas-phase flow rate and liquid-phase flow rates of oil-gas-water three-phase flow
CN104090020A (en) * 2014-07-10 2014-10-08 天津大学 Electric and ultrasonic-based bimodal multiphase flow measuring device
CN104089985A (en) * 2014-07-10 2014-10-08 天津大学 Visual multiphase flow test method based on electric and ultrasonic sensing principle
CN104155471A (en) * 2014-07-10 2014-11-19 天津大学 Multiphase flow testing method based on cross-correlation velocity measurement of ultrasonic and electric multiple sensors
CN105181996A (en) * 2015-08-12 2015-12-23 天津大学 Two-phase flow velocity acoustoelectric bimodal measuring method
CN106932606A (en) * 2017-02-27 2017-07-07 天津大学 Gas-liquid two-phase slug flow and slug flow flow velocity acoustic-electric bimodal measuring method
CN107024603A (en) * 2017-02-27 2017-08-08 天津大学 Gas-liquid two-phase bubbly flow flow velocity acoustic-electric bimodal measuring method
CN106996988A (en) * 2017-03-20 2017-08-01 天津大学 Oil gas water three phase plug-like disperses flow velocity measuring method
CN107153086A (en) * 2017-04-28 2017-09-12 天津大学 The acoustic-electric bimodal fusion measuring method of two phase flow procedure parameter

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
XIAOXIAO DONG ET: "Measuring Oil-Water Two-Phase Flow Velocity With Continuous-Wave Ultrasound Doppler Sensor and Drift-Flux Modle", 《IEE TRANSACTIONS OV INSTRUMENTATION AND MEASUREMENT》 *
戴玮等: "水平管道油气水三相流含水率测量", 《传感器与微系统》 *
韩骏等: "小波神经网络在油水两相流软测量中的应用", 《天津大学学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114563470A (en) * 2022-01-28 2022-05-31 天津大学 Oil-water two-phase flow phase content measuring method based on nonlinear ultrasound
DE102022105759A1 (en) 2022-03-11 2023-09-14 Krohne Ag Flow measuring device and method for operating the same
CN115106499A (en) * 2022-06-30 2022-09-27 北京科技大学 A method and system for judging abnormal fluctuation of crystallizer liquid level
CN115106499B (en) * 2022-06-30 2024-02-20 北京科技大学 A method and system for identifying abnormal fluctuations in crystallizer liquid level
CN118655184A (en) * 2024-08-21 2024-09-17 南京尤力浦能源科技有限公司 A non-invasive multiphase flow content detection method and device based on continuous wave
CN118655184B (en) * 2024-08-21 2024-11-15 南京尤力浦能源科技有限公司 Continuous wave-based non-invasive multiphase flow content detection method and device

Also Published As

Publication number Publication date
CN109188016B (en) 2020-05-05

Similar Documents

Publication Publication Date Title
Shi et al. Flow rate measurement of oil-gas-water wavy flow through a combined electrical and ultrasonic sensor
CN104101687B (en) Based on the polyphasic flow visual testing method of ultrasonic Doppler and electricity multisensor
Garbini et al. Measurement of fluid turbulence based on pulsed ultrasound techniques. Part 2. Experimental investigation
CN109188016A (en) Oil-gas-water three-phase flow split-phase flow velocity acoustic-electric bimodal measurement method
JP3110048B2 (en) Method and apparatus for analyzing two-phase flow
CN104155471B (en) Multiphase flow testing method based on cross-correlation velocity measurement of ultrasonic and electric multiple sensors
Zhai et al. The ultrasonic measurement of high water volume fraction in dispersed oil-in-water flows
CN104965104B (en) Two phase flow split-phase flow velocity acoustic-electric bimodal measuring method
CN104155358A (en) Visual multiphase flow testing device with combination of ultrasonic Doppler and electrical sensor
JP2003526101A (en) Method and apparatus for simultaneously measuring the flow rate and concentration of a multiphase liquid / gas mixture
CN107632167B (en) Two-phase Flow Velocity Measurement Method Based on Ultrasonic Pulse Doppler and Electrical Multi-sensor
Meribout et al. A multisensor intelligent device for real-time multiphase flow metering in oil fields
CN105181793B (en) Method based on ultrasonic wave two-frequency signal measurement two phase flow void fraction
CN106813108A (en) A kind of leakage locating method based on speed difference
CN112067691A (en) Multi-frequency ultrasonic tomography method based on attenuation spectrum fusion of oil, gas and water three-phase medium
CN106290977B (en) Processing method for obtaining water flow velocity signal by using Doppler ultrasonic current meter
CN105181996B (en) Two-phase flow velocity acoustic-electric bimodal measuring method
Wang et al. A new method for processing ultrasonic gas flowmeter signal in wet gas
o'Sullivan et al. Ultrasonic measurement of gas flow using electrostatic transducers
CN107389974B (en) Gas-Liquid Slug Flow structure flow velocity acoustic-electric bimodal measurement method
CN112129832A (en) Sweep frequency ultrasonic attenuation measurement method for detecting content of dispersed phase
CN106443060B (en) Two-phase flow velocity measurement method based on continuous wave ultrasonic Doppler spectrum correction
Rychagov et al. Multipath flowrate measurements of symmetric and asymmetric flows
CN106932476A (en) A kind of device and method for measuring solid-liquid two-phase flow phase content
Abda et al. Ultrasonic device for real-time sewage velocity and suspended particles concentration measurements

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant