CN105806780A - Multi-vision-based system and method for measuring gas-liquid two-phase flow phase fraction of small channel - Google Patents

Multi-vision-based system and method for measuring gas-liquid two-phase flow phase fraction of small channel Download PDF

Info

Publication number
CN105806780A
CN105806780A CN201610136061.XA CN201610136061A CN105806780A CN 105806780 A CN105806780 A CN 105806780A CN 201610136061 A CN201610136061 A CN 201610136061A CN 105806780 A CN105806780 A CN 105806780A
Authority
CN
China
Prior art keywords
flow
image
phase
phase content
bubble
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
CN201610136061.XA
Other languages
Chinese (zh)
Other versions
CN105806780B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201610136061.XA priority Critical patent/CN105806780B/en
Publication of CN105806780A publication Critical patent/CN105806780A/en
Application granted granted Critical
Publication of CN105806780B publication Critical patent/CN105806780B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light

Landscapes

  • Measuring Volume Flow (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)

Abstract

The invention discloses a multi-vision-based system and method for measuring gas-liquid two-phase flow phase fraction of a small channel. The system comprises a high-speed camera, two light sources, a right-angle prism, a small pipeline, a glycerinum groove and a computer. The method comprises the following steps: firstly, simultaneously obtaining gas-liquid two-phase flow flowing image sequences of the small channel from two vertical angles; secondly, carrying out image pretreatment and optical distortion correction and building phase fraction measurement models of four typical flow patterns, namely bubble flow, slug flow, annular flow and wavy flow; thirdly, extracting feature vectors of images and carrying out flow pattern identification; and finally, selecting the phase fraction measurement model of the corresponding flow pattern and calculating the phase fraction. By a multi-vision technology, the phase fraction measurement model is built for each typical flow pattern, the flow pattern information is introduced into the phase fraction measurement model, and thus the measurement accuracy is improved. The corresponding device has the advantages of simplicity in structure, non-contact measurement, visualization and the like.

Description

Passage aisle gas-liquid two-phase flow containing rate based on many visions measures system and method
Technical field
The present invention relates to multiphase flow parameter measurement field, particularly relate to a kind of passage aisle gas-liquid two-phase flow containing rate based on many visions and measure system and method.
Background technology
In recent years, along with the emergence of micro-processing technology and new material technology, industrial equipment presents the development trend of microminiaturization, miniaturization.Owing to microminiaturized, compact apparatus have efficient, quick, the easily-controllable and advantage such as economize in raw materials, it is widely used in fields such as biology, medical treatment, chemical industry, pharmacy.Therefore, passage aisle gas-liquid two-phase flow parameter measurement becomes a current study hotspot and difficulties.
Compared to regular channels, passage aisle gas-liquid two-phase flow parameter measurement has certain difficulty.Under miniature scale, due to the reduction of hydraulic diameter so that the impact of surface tension and viscosity effect becomes big, and the impact of action of gravity weakens relatively, and theoretical model and empirical equation in regular channels are no longer applicable.
At present, existing passage aisle gas-liquid two-phase flow containing rate measuring method mainly has fast valve method, optical method, electric method and high speed video process etc..By contrast, high speed video process has intuitive, the advantage such as untouchable, is widely used in passage aisle two phase flow phase content is measured.But, the existing phase content measuring method majority based on high-speed camera is to obtain two phase flow image from an angle, and the shape information of gas-liquid two-phase spatial information and different angles is obtained deficiency.Many vision techniques utilize multiple camera or single camera and optics to constitute system and obtain the image information of object of study from multiple angles, efficiently solve the problem that single vision acquisition of information is not enough, provide probability for solving complex parameters measurement problem.
Summary of the invention
Not enough for the existing passage aisle gas-liquid two-phase flow containing rate measuring method acquisition of information based on high-speed photography, the problem that order of accuarcy is not high enough, the present invention is by repeatedly studying, for passage aisle biphase gas and liquid flow, propose a kind of passage aisle gas-liquid two-phase flow containing rate based on many visions and measure system and method, adopt many visual systemes to obtain two images being mutually perpendicular to angle of two phase flow simultaneously, carry out Image semantic classification and optic aberrance revising, set up the phase content measurement model of four kinds of typical flow patterns, extract characteristic vector and carry out meteor trail echoes, select the phase content measurement model of corresponding flow pattern, calculate phase content.
Passage aisle gas-liquid two-phase flow containing rate measurement system based on many visions includes passage aisle pipe-line system, many visions high-speed image sampling unit and microcomputer.Passage aisle pipe-line system is made up of the passage aisle being immersed in glycerol groove.Many visions high-speed image sampling unit is made up of two LED light sources, corner cube prism, high-speed cameras.Corner cube prism is placed near glycerol groove, and one right-angle side is perpendicular to the imaging plane of high-speed camera.Two light sources provide the backlight of two vertical direction, and the high-speed camera being positioned over glycerol groove front photographs two two phase flow images being mutually perpendicular to angle simultaneously, and the image transmitting collected calculates for image procossing and phase content in microcomputer.
Passage aisle gas-liquid two-phase flow containing rate measuring method step based on many visions is as follows:
The first step, obtains image: utilize many visions high-speed image sampling unit, obtains two passage aisle biphase gas and liquid flow flow image sequences being mutually perpendicular to angle simultaneously;
Second step, sets up phase content measurement model: carry out Image semantic classification and optic aberrance revising, and four kinds of typical flow patterns are modeled respectively;
3rd step, meteor trail echoes: extract the characteristic vector of image, carry out meteor trail echoes;
4th step, calculates phase content: select the phase content measurement model of corresponding flow pattern, calculates phase content.
Wherein, specifically comprising the following steps that of phase content measurement model is set up
1) Image semantic classification.Filled by difference shadow method, medium filtering, rim detection and binaryzation and carry out pretreatment, obtain bianry image.
2) optic aberrance revising.Correction coefficient is obtained, it is achieved the correction of biphase gas and liquid flow flow image by optical path analysis.For circular glass pipeline, its tube wall refraction rate and glycerol approximately equal.Pipeline is immersed in glycerol groove and can reduce the impact that captured image is caused by tube wall.As it is shown on figure 3, the foursquare lucite cuboid glycerol groove adopting cross section to be 25mm × 25mm, being immersed in by pipeline in glycerol groove, the center of pipeline center position and cross-sectional square shape coincides.For object point P in pipeline, (x, y), this point there is first time refraction at liquid phase and tube wall interface in the light sent, and in glycerol groove cell wall and Air Interface, second time refraction occurs.Two adjacent ray sent by P intersection point P'(x', the y' of reverse extending line after second time refraction) be a P to it should be observed that picture point, may thereby determine that the correction coefficient alpha=y/y' at a P place.After obtaining correction coefficient, bianry image is corrected.
3) four kinds of typical flow pattern modelings.Use for reference the thought of " segmentation, approximate summation ", set up bubble flow and slug flow phase content measurement model, annular flow phase content measurement model, wave flow phase content measurement model respectively.
Bubble flow and slug flow phase content measurement model: owing to the hydraulic diameter of passage aisle is less, the same cross section of pipeline only exists a bubble simultaneously, do not have same cross section in conventional pipeline and have the situation of multiple bubble, it is possible to determine the size of bubble in both direction by the image of two angles.Consider the situation of bubble flow and slug flow, at the same section of pipeline, be oval by region equivalent shared by gas phase.Within the scope of the axial length in pixels of pipeline, think that the area of section shared by gas phase is equal, therefore gas phase in this segment limit can be equivalent to an Elliptic Cylinder, and the major and minor axis of this Elliptic Cylinder bottom surface ellipse is determined by two bianry images being mutually perpendicular to angle respectively.Thus, it is possible to each bubble or slug to be equivalent to the superposition of several Elliptic Cylinders, each bubble or the volume of gas phase corresponding to vent plug can be tried to achieve.
As shown in fig. 5, it is assumed that captured single set of image exists n bubble/vent plug.Number of pixels shared by i-th bubble/vent plug length is li, number of pixels respectively a shared by the major and minor axis of the jth pixel place bubble section ellipse that i-th bubble/vent plug is correspondingijWith bij, then the volume V of gas phase in captured image correspondence pipeline section can be tried to achievegas1For:
V g a s 1 = π 4 μ 3 Σ i = 1 n Σ j = 1 l i α i j β i j a i j b i j
Wherein, μ is the physical length that each square pixel is corresponding, αijFor aijCorresponding correction coefficient, βijFor bijCorresponding correction coefficient.The then phase content ε of bubble flow/slug flow in corresponding pipeline section1For:
ϵ 1 = V g a s 1 V p i p e = Σ i = 1 n Σ j = 1 l i α i j β i j a i j b i j LD 2
Wherein, VpipeFor the conduit volume that image taking arrives, L is the duct length that image taking arrives, and D is the pipe diameter that image taking arrives, and the physical length μ that phase content is corresponding with pixel is unrelated, and only to each dimensional parameters of two phase flow image shared by the number of pixel relevant.
In actual photographed process, owing to length of pipe section taken in viewing field of camera is shorter, the calculated phase content of single set of image can not accurately reflect the phase content under current flow pattern.Therefore, choose the image construction image sequence being continuously shot by N group, take the meansigma methods of image sequence phase contentAs the phase content measured value of bubble flow/slug flow under corresponding experiment condition, it may be assumed that
ϵ ‾ b u b b l e / s l u g = Σ k = 1 N ϵ k N
Annular flow phase content measurement model: in horizontal pipe, by the impact of action of gravity, duct bottom liquid film is thicker;And liquid film is distributed around even tube wall in vertical pipeline.Similar with slug flow and bubble flow, under annular flow flow pattern, gas phase in pipeline is approximately the superposition of several Elliptic Cylinders, gaseous phase volume can be tried to achieve accordingly, it is hereby achieved that phase content.As shown in Figure 6, in acquired annular flow bianry image, number of pixels respectively a shared by the major and minor axis that gas phase cross section, ith pixel place is ovaliWith bi, then single set of image correspondence pipeline section inner annular flow containing rate ε2For:
ϵ 2 = V g a s 2 V p i p e = Σ i = 1 L α i β i a i b i LD 2
Wherein, αiFor aiCorresponding correction coefficient, βiFor biCorresponding correction coefficient.
Take the N group image construction image sequence being continuously shot, calculate the phase content that often group image is corresponding, then averageAs the measured value of annular flow phase content under corresponding experiment condition, it may be assumed that
ϵ ‾ a n n u l a r = Σ k = 1 N ϵ k N
Wave flow phase content measurement model: under wave flow flow pattern, in pipeline, the volume of gas phase can be obtained by the volume superposition of several cylinders.Different from other three kinds of typical flow patterns, owing to can uniquely determine the area of section shared by gas phase and liquid phase after determining height shared by gas phase and liquid phase, therefore in horizontal pipe, the image of the single angle that the phase content of wave flow can be shot by pipeline side is determined.As it is shown in fig. 7, in acquired wave flow image, the height in gas phase cross section, ith pixel place is hi, then the volume V of gas phase in single image correspondence pipeline sectiongas3For:
V g a s 3 = μ 3 Σ i = 1 L [ π 4 D 2 - D 2 4 arccos ( 2 h i D - 1 ) + ( h i - D 2 ) h i D - h i 2 ]
The then phase content ε of wave flow in single image correspondence pipeline section3For:
ϵ 3 = V g a s 3 V p i p e = 1 L Σ i = 1 L [ 1 - 1 π arccos ( 2 h i D - 1 ) + 2 π D ( 2 h i - D ) h i D - h i 2 ]
Take the N width image construction image sequence being continuously shot, calculate the phase content that each image is corresponding, then averageAs the measured value of annular flow phase content under corresponding experiment condition, it may be assumed that
ϵ ‾ w a v y = Σ k = 1 N ϵ k N .
Accompanying drawing explanation
Accompanying drawing 1 is the present invention passage aisle gas-liquid two-phase flow containing rate measurement system schematic based on multi-visual information integration technology, in Fig. 1, and 1-passage aisle pipe-line system, 2-many visions high-speed image sampling unit, 3-microcomputer.
Accompanying drawing 2 is many visions high-speed image sampling unit and passage aisle pipe-line system structural representation, in Fig. 2, and 4-the first LED light source, 5-the second LED light source, 6-glycerol groove, 7-passage aisle, 8-corner cube prism, 9-high-speed camera.
Accompanying drawing 3 is the index path that light is propagated in passage aisle and glycerol groove.
Accompanying drawing 4 is based on the passage aisle gas-liquid two-phase flow containing rate measurement procedure figure of many visions.
Accompanying drawing 5 is bubble flow and slug flow phase content measurement model schematic diagram.
Accompanying drawing 6 is annular flow phase content measurement model schematic diagram.
Accompanying drawing 7 is wave flow phase content measurement model schematic diagram.
Accompanying drawing 8 is the phase content measurement result figure of the horizontal pipe of 4.0mm internal diameter and vertical pipeline.
Detailed description of the invention
With reference to accompanying drawing 1, the present invention utilizes the passage aisle gas-liquid two-phase flow containing rate based on many visions to measure system, it is achieved passage aisle gas-liquid two-phase flow containing rate is measured.Passage aisle gas-liquid two-phase flow containing rate measurement system based on many visions includes passage aisle pipe-line system, many visions high-speed image sampling unit and microcomputer three part.
As shown in Figure 2, passage aisle pipe-line system 1 is made up of the passage aisle 7 being immersed in glycerol groove 6, many visions high-speed image sampling unit 2 is made up of high-speed camera 9, corner cube prism the 8, first LED light source the 4, second LED light source 5, many visions high-speed image sampling unit is responsible for two biphase gas and liquid flow flow images being mutually perpendicular to angle of collection simultaneously, many visions high-speed image sampling unit 2 image transmitting collected carries out image procossing to computer 3 and phase content calculates.
The first described LED light source 4 is perpendicular to passage aisle 6 and places, second LED light source 5 is both perpendicular to the first LED light source 5 and passage aisle 6 and in being generally aligned in the same plane with the first LED light source 5, corner cube prism 8 is positioned at passage aisle 6 other side relative to the first LED light source 4, high-speed camera 9 is positioned at passage aisle 6 other side relative to the second LED light source 5, and a right-angle side of corner cube prism 8 is perpendicular to the imaging plane of high-speed camera 6.
With reference to accompanying drawing 4, carry out phase content calculating according to the passage aisle gas-liquid two-phase flow containing rate measurement procedure figure based on many visions, concretely comprise the following steps:
The first step, obtains image: utilize many visions high-speed image sampling unit, obtains two passage aisle biphase gas and liquid flow flow image sequences being mutually perpendicular to angle simultaneously;
Second step, sets up phase content measurement model: carry out Image semantic classification and optic aberrance revising, and four kinds of typical flow patterns are modeled respectively;
3rd step, meteor trail echoes: extract the characteristic vector of image, carry out meteor trail echoes;
4th step, calculates phase content: select the phase content measurement model of corresponding flow pattern, calculates phase content.
Wherein, second step sets up concretely comprising the following steps of phase content measurement model:
1) Image semantic classification: filled by difference shadow method, medium filtering, rim detection and binaryzation and carry out pretreatment, obtain bianry image.
2) optic aberrance revising: obtain correction coefficient by optical path analysis, it is achieved the correction to bianry image.Adopting cross section is foursquare lucite cuboid glycerol groove, is immersed in by pipeline in glycerol groove, and the center of pipeline center position and cross-sectional square shape coincides.For object point P in pipeline, (x, y), this point there is first time refraction at liquid phase and tube wall interface in the light sent, and in glycerol groove cell wall and Air Interface, second time refraction occurs.Two adjacent ray sent by P intersection point P'(x', the y' of reverse extending line after second time refraction) be a P to it should be observed that picture point, may thereby determine that the correction coefficient alpha=y/y' at a P place.After obtaining correction coefficient, to correct image.
3) four kinds of typical flow pattern modelings: use for reference the thought of " segmentation, approximate summation ", set up bubble flow and slug flow phase content measurement model, annular flow phase content measurement model, wave flow phase content measurement model respectively.
Bubble flow and slug flow phase content measurement model: owing to the hydraulic diameter of passage aisle is less, the same cross section of pipeline only exists a bubble simultaneously, do not have same cross section in conventional pipeline and have the situation of multiple bubble, it is possible to determine the size of bubble in both direction by the image of two angles.Consider the situation of bubble flow and slug flow, at the same section of pipeline, be oval by region equivalent shared by gas phase.Within the scope of the axial length in pixels of pipeline, think that the area of section shared by gas phase is equal, therefore gas phase in this segment limit can be equivalent to an Elliptic Cylinder, and the major and minor axis of this Elliptic Cylinder bottom surface ellipse is determined by two bianry images being mutually perpendicular to angle respectively.Thus, it is possible to each bubble or slug to be equivalent to the superposition of several Elliptic Cylinders, each bubble or the volume of gas phase corresponding to vent plug can be tried to achieve.
As shown in fig. 5, it is assumed that captured single set of image exists n bubble/vent plug.Number of pixels shared by i-th bubble/vent plug length is li, number of pixels respectively a shared by the major and minor axis of the jth pixel place bubble section ellipse that i-th bubble/vent plug is correspondingijWith bij, then the volume V of gas phase in captured image correspondence pipeline section can be tried to achievegas1For:
V g a s 1 = π 4 μ 3 Σ i = 1 n Σ j = 1 l i α i j β i j a i j b i j
Wherein, μ is the physical length that each square pixel is corresponding, αijFor aijCorresponding correction coefficient, βijFor bijCorresponding correction coefficient.The then phase content ε of bubble flow/slug flow in corresponding pipeline section1For:
ϵ 1 = V g a s 1 V p i p e = Σ i = 1 n Σ j = 1 l i α i j β i j a i j b i j LD 2
Wherein, VpipeFor the conduit volume that image taking arrives, L is the duct length that image taking arrives, and D is the pipe diameter that image taking arrives, and the physical length μ that phase content is corresponding with pixel is unrelated, and only to each dimensional parameters of two phase flow image shared by the number of pixel relevant.
In actual photographed process, owing to length of pipe section taken in viewing field of camera is shorter, the calculated phase content of single set of image can not accurately reflect the phase content under current flow pattern.Therefore, choose the image construction image sequence being continuously shot by N group, take the meansigma methods of image sequence phase contentAs the phase content measured value of bubble flow/slug flow under corresponding experiment condition, it may be assumed that
ϵ ‾ b u b b l e / s l u g = Σ k = 1 N ϵ k N
Annular flow phase content measurement model: in horizontal pipe, by the impact of action of gravity, duct bottom liquid film is thicker;And liquid film is distributed around even tube wall in vertical pipeline.Similar with slug flow and bubble flow, under annular flow flow pattern, gas phase in pipeline is approximately the superposition of several Elliptic Cylinders, gaseous phase volume can be tried to achieve accordingly, it is hereby achieved that phase content.As shown in Figure 6, in acquired annular flow bianry image, number of pixels respectively a shared by the major and minor axis that gas phase cross section, ith pixel place is ovaliWith bi, then single set of image correspondence pipeline section inner annular flow containing rate ε2For:
ϵ 2 = V g a s 2 V p i p e = Σ i = 1 L α i β i a i b i LD 2
Wherein, αiFor aiCorresponding correction coefficient, βiFor biCorresponding correction coefficient.
Take the N group image construction image sequence being continuously shot, calculate the phase content that often group image is corresponding, then averageAs the measured value of annular flow phase content under corresponding experiment condition, it may be assumed that
ϵ ‾ a n n u l a r = Σ k = 1 N ϵ k N
Wave flow phase content measurement model: under wave flow flow pattern, in pipeline, the volume of gas phase can be obtained by the volume superposition of several cylinders.Different from other three kinds of typical flow patterns, owing to can uniquely determine the area of section shared by gas phase and liquid phase after determining height shared by gas phase and liquid phase, therefore in horizontal pipe, the image of the single angle that the phase content of wave flow can be shot by pipeline side is determined.As it is shown in fig. 7, in acquired wave flow image, the height in gas phase cross section, ith pixel place is hi, then the volume V of gas phase in single image correspondence pipeline sectiongas3For:
V g a s 3 = μ 3 Σ i = 1 L [ π 4 D 2 - D 2 4 arccos ( 2 h i D - 1 ) + ( h i - D 2 ) h i D - h i 2 ]
The then phase content ε of wave flow in single image correspondence pipeline section3For:
ϵ 3 = V g a s 3 V p i p e = 1 L Σ i = 1 L [ 1 - 1 π arccos ( 2 h i D - 1 ) + 2 π D ( 2 h i - D ) h i D - h i 2 ]
Take the N width image construction image sequence being continuously shot, calculate the phase content that each image is corresponding, then averageAs the measured value of annular flow phase content under corresponding experiment condition, it may be assumed that
ϵ ‾ w a v y = Σ k = 1 N ϵ k N .
Utilize in the present invention propose system and method internal diameter be 4.0mm horizontally and vertically transparent passage aisle pipe-line system in carried out phase content experiments of measuring.
In experiment, the excursion of gas phase flow rate is 6ml/min-41700ml/min, liquid phase changes in flow rate ranges for 6ml/min-4170ml/min, having carried out the phase content under bubble flow, slug flow, annular flow and four kinds of typical flow patterns of wave flow in horizontal pipe to measure, the phase content having carried out other three kinds of typical flow patterns except wave flow in vertical pipeline is measured.Adopt the phase content that fast valve method records as calibration value.
As shown in Figure 8, a () figure is the phase content experiments of measuring result in 4.0mm horizontal pipe, b () figure is the phase content experiments of measuring result in 4.0mm vertical pipeline, it can be seen that horizontally and vertically the phase content in small pipeline is measured absolute error and is respectively less than 6%.Test result indicate that, it is feasible, effective that the passage aisle gas-liquid two-phase flow containing rate based on many visions proposed by the invention measures system and method.

Claims (5)

1. the passage aisle gas-liquid two-phase flow containing rate based on many visions measures system, it is characterized in that: this system is by passage aisle pipe-line system (1), many visions high-speed image sampling unit (2) and microcomputer (3) composition, passage aisle pipe-line system (1) is made up of the passage aisle (7) being immersed in glycerol groove (6), many visions high-speed image sampling unit (2) are by high-speed camera (9), corner cube prism (8), first LED light source (4), second LED light source (5) forms, many visions high-speed image sampling unit is responsible for two biphase gas and liquid flow flow images being mutually perpendicular to angle of collection simultaneously, the image transmitting collected by many visions high-speed image sampling unit (2) carries out image procossing to computer (3) and phase content calculates;Described the first LED light source (4) is perpendicular to passage aisle (6) and places, second LED light source (5) is both perpendicular to the first LED light source (5) and passage aisle (6) and in being generally aligned in the same plane with the first LED light source (5), corner cube prism (8) is positioned at passage aisle (6) other side relative to the first LED light source (4), high-speed camera (9) is positioned at passage aisle (6) other side relative to the second LED light source (5), and a right-angle side of corner cube prism (8) is perpendicular to the imaging plane of high-speed camera (6).
2. the passage aisle gas-liquid two-phase flow containing rate measuring method based on many visions based on system described in claim 1, it is characterised in that the method comprises the steps:
1) obtain image: utilize many visions high-speed image sampling unit, obtain two passage aisle biphase gas and liquid flow flow image sequences being mutually perpendicular to angle simultaneously;
2) set up phase content measurement model: carry out Image semantic classification and optic aberrance revising, four kinds of typical flow patterns are modeled respectively;
3) meteor trail echoes: extract the characteristic vector of image, carry out meteor trail echoes;
4) calculate phase content: select the phase content measurement model of corresponding flow pattern, calculate phase content.
3. phase content measuring method according to claim 2, it is characterised in that described step 2) in the method setting up phase content measurement model particularly as follows:
1) Image semantic classification: filled by image segmentation, image denoising, rim detection and binaryzation and carry out pretreatment, obtain bianry image;
2) optic aberrance revising: obtain correction coefficient by optical path analysis, it is achieved the correction to bianry image;
3) four kinds of typical flow pattern modelings: set up bubble flow and slug flow phase content measurement model, annular flow phase content measurement model, wave flow phase content measurement model respectively.
4. according to Claims 2 or 3, set up phase content measurement model method, it is characterized in that described step 2) in optic aberrance revising method particularly as follows: adopt cross section be foursquare lucite cuboid glycerol groove, pipeline is immersed in glycerol groove, the center making pipeline center position and cross-sectional square shape coincides, for object point P (x in pipeline, y), first time refraction is there is in the light sent by this point at liquid phase and tube wall interface, second time refraction is there is in glycerol groove cell wall and Air Interface, analysis light path is known, the intersection point P'(x' of two adjacent ray sent by P reverse extending line after second time refraction, y') picture point corresponding to a P it is, so that it is determined that the correction coefficient alpha=y/y' at some P place, after obtaining correction coefficient, bianry image is corrected.
5. according to Claims 2 or 3, set up phase content measurement model method, it is characterised in that described step 3) in four kinds of typical flow patterns modelings method particularly as follows:
1) bubble flow and slug flow phase content measurement model: the hydraulic diameter in passage aisle is less, the same cross section of pipeline only exists a bubble simultaneously, the size of bubble in both direction can be determined by the image of two angles, consider the situation of bubble flow and slug flow, at the same section of pipeline, it is oval by region equivalent shared by gas phase, within the scope of the axial length in pixels of pipeline, think that the area of section shared by gas phase is equal, therefore gas phase in this segment limit can be equivalent to an Elliptic Cylinder, the major and minor axis of this Elliptic Cylinder bottom surface ellipse is determined by two bianry images being mutually perpendicular to angle respectively, thus, each bubble or slug can be equivalent to the superposition of several Elliptic Cylinders, each bubble or the volume of gas phase corresponding to vent plug can be tried to achieve,
Assuming there is n bubble/vent plug in captured single set of image, number of pixels shared by i-th bubble/vent plug length is li, number of pixels respectively a shared by the major and minor axis of the jth pixel place bubble section ellipse that i-th bubble/vent plug is correspondingijWith bij, then the volume V of gas phase in captured image correspondence pipeline section can be tried to achievegas1For:
V g a s 1 = π 4 μ 3 Σ i = 1 n Σ j = 1 l i α i j β i j a i j b i j
Wherein, μ is the physical length that each square pixel is corresponding, αijFor aijCorresponding correction coefficient, βijFor bijCorresponding correction coefficient, then the phase content ε of bubble flow/slug flow in corresponding pipeline section1For:
ϵ 1 = V g a s 1 V p i p e = Σ i = 1 n Σ j = 1 l i α i j β i j a i j b i j LD 2
Wherein, VpipeFor the conduit volume that image taking arrives, L is the duct length that image taking arrives, and D is the pipe diameter that image taking arrives, and the physical length μ that phase content is corresponding with pixel is unrelated, and only to each dimensional parameters of two phase flow image shared by the number of pixel relevant;
In actual photographed process, owing to length of pipe section taken in viewing field of camera is shorter, the calculated phase content of single set of image can not accurately reflect the phase content under current flow pattern, therefore, choose the image construction image sequence being continuously shot by N group, take the meansigma methods of image sequence phase contentAs the phase content measured value of bubble flow/slug flow under corresponding experiment condition, it may be assumed that
ϵ ‾ b u b b l e / s l u g = Σ k = 1 N ϵ k N
2) annular flow phase content measurement model: in horizontal pipe, by the impact of action of gravity, duct bottom liquid film is thicker;And liquid film is distributed around even tube wall in vertical pipeline, similar with slug flow and bubble flow, under annular flow flow pattern, gas phase in pipeline is approximately the superposition of several Elliptic Cylinders, gaseous phase volume can be tried to achieve accordingly, it is hereby achieved that phase content, in acquired annular flow bianry image, number of pixels respectively a shared by the major and minor axis that gas phase cross section, ith pixel place is ovaliWith bi, then single set of image correspondence pipeline section inner annular flow containing rate ε2For:
ϵ 2 = V g a s 2 V p i p e = Σ i = 1 L α i β i a i b i LD 2
Wherein, αiFor aiCorresponding correction coefficient, βiFor biCorresponding correction coefficient;
Take the N group image construction image sequence being continuously shot, calculate the phase content that often group image is corresponding, then averageAs the measured value of annular flow phase content under corresponding experiment condition, it may be assumed that
ϵ ‾ a n n u l a r = Σ k = 1 N ϵ k N
3) wave flow phase content measurement model: under wave flow flow pattern, in pipeline, the volume of gas phase can be obtained by the volume superposition of several cylinders, different from other three kinds of typical flow patterns, owing to can uniquely determine the area of section shared by gas phase and liquid phase after determining height shared by gas phase and liquid phase, therefore in horizontal pipe, the image of the single angle that the phase content of wave flow can be shot by pipeline side is determined, in acquired wave flow image, the height in gas phase cross section, ith pixel place is hi, then the volume V of gas phase in single image correspondence pipeline sectiongas3For:
V g a s 3 = μ 3 Σ i = 1 L [ π 4 D 2 - D 2 4 arccos ( 2 h i D - 1 ) + ( h i - D 2 ) h i D - h i 2 ]
The then phase content ε of wave flow in single image correspondence pipeline section3For:
ϵ 3 = V g a s 3 V p i p e = 1 L Σ i = 1 L [ 1 - 1 π arccos ( 2 h i D - 1 ) + 2 π D ( 2 h i - D ) h i D - h i 2 ]
Take the N width image construction image sequence being continuously shot, calculate the phase content that each image is corresponding, then averageAs the measured value of annular flow phase content under corresponding experiment condition, it may be assumed that
ϵ ‾ w a v y = Σ k = 1 N ϵ k N .
CN201610136061.XA 2016-03-10 2016-03-10 The phase content measurement method of passage aisle gas-liquid two-phase flow containing rate measuring system based on more visions Active CN105806780B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610136061.XA CN105806780B (en) 2016-03-10 2016-03-10 The phase content measurement method of passage aisle gas-liquid two-phase flow containing rate measuring system based on more visions

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610136061.XA CN105806780B (en) 2016-03-10 2016-03-10 The phase content measurement method of passage aisle gas-liquid two-phase flow containing rate measuring system based on more visions

Publications (2)

Publication Number Publication Date
CN105806780A true CN105806780A (en) 2016-07-27
CN105806780B CN105806780B (en) 2018-09-07

Family

ID=56468123

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610136061.XA Active CN105806780B (en) 2016-03-10 2016-03-10 The phase content measurement method of passage aisle gas-liquid two-phase flow containing rate measuring system based on more visions

Country Status (1)

Country Link
CN (1) CN105806780B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106482784A (en) * 2016-10-11 2017-03-08 天津大学 Two phase flow spatial complex network visualization analysis method based on grid sensor
CN107121287A (en) * 2017-05-25 2017-09-01 哈尔滨电机厂有限责任公司 A kind of many seat in the plane principal and subordinate triggered synchro imaging methods of model turbine fluidised form observation
CN108844959A (en) * 2018-04-12 2018-11-20 西安交通大学 The measurement of gas-liquid two-phase ring-type flow section phase content and modification method in a kind of round tube
CN111307809A (en) * 2020-02-24 2020-06-19 浙江大学 Small pipeline gas-liquid two-phase flow phase distribution optical detection system and method

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0875744A2 (en) * 1997-04-30 1998-11-04 British Aerospace Public Limited Company Method and apparatus for inspecting a test piece by optical metrology
CN101231300A (en) * 2008-02-04 2008-07-30 浙江大学 Apparatus and method for measuring microtubule gas-liquid diphasic flow rate
CN102590030A (en) * 2012-01-18 2012-07-18 浙江大学 Small-passage gas-liquid phase flow pattern identification device and method based on photovoltaic array sensor
CN102680203A (en) * 2012-05-09 2012-09-19 浙江大学 Micro-channel gas-liquid two-phase flow voidage measuring device and method
CN104457703A (en) * 2014-11-28 2015-03-25 浙江大学 Multiple-vision-based System and method for measuring parameters of gas-liquid two-phase flow in small channel based on
CN105426830A (en) * 2015-11-11 2016-03-23 浙江大学 Small-channel gas liquid two-phase flow pattern identification system and method based on multi-vision information fusion technology

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0875744A2 (en) * 1997-04-30 1998-11-04 British Aerospace Public Limited Company Method and apparatus for inspecting a test piece by optical metrology
CN101231300A (en) * 2008-02-04 2008-07-30 浙江大学 Apparatus and method for measuring microtubule gas-liquid diphasic flow rate
CN102590030A (en) * 2012-01-18 2012-07-18 浙江大学 Small-passage gas-liquid phase flow pattern identification device and method based on photovoltaic array sensor
CN102680203A (en) * 2012-05-09 2012-09-19 浙江大学 Micro-channel gas-liquid two-phase flow voidage measuring device and method
CN104457703A (en) * 2014-11-28 2015-03-25 浙江大学 Multiple-vision-based System and method for measuring parameters of gas-liquid two-phase flow in small channel based on
CN105426830A (en) * 2015-11-11 2016-03-23 浙江大学 Small-channel gas liquid two-phase flow pattern identification system and method based on multi-vision information fusion technology

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张立伟等: "小通道小流量气液两相流热式测量方法的研究", 《仪器仪表学报》 *
林翀等: "气液两相流相含率测量新方法研究", 《工程热物理学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106482784A (en) * 2016-10-11 2017-03-08 天津大学 Two phase flow spatial complex network visualization analysis method based on grid sensor
CN106482784B (en) * 2016-10-11 2018-11-16 天津大学 Two phase flow spatial complex network visualization analysis method based on grid sensor
CN107121287A (en) * 2017-05-25 2017-09-01 哈尔滨电机厂有限责任公司 A kind of many seat in the plane principal and subordinate triggered synchro imaging methods of model turbine fluidised form observation
CN108844959A (en) * 2018-04-12 2018-11-20 西安交通大学 The measurement of gas-liquid two-phase ring-type flow section phase content and modification method in a kind of round tube
CN108844959B (en) * 2018-04-12 2020-05-22 西安交通大学 Method for measuring and correcting gas-liquid two-phase annular flow section phase content in circular tube
CN111307809A (en) * 2020-02-24 2020-06-19 浙江大学 Small pipeline gas-liquid two-phase flow phase distribution optical detection system and method
CN111307809B (en) * 2020-02-24 2021-01-12 浙江大学 Small pipeline gas-liquid two-phase flow phase distribution optical detection system and method

Also Published As

Publication number Publication date
CN105806780B (en) 2018-09-07

Similar Documents

Publication Publication Date Title
CN109544679B (en) Three-dimensional reconstruction method for inner wall of pipeline
CN109489620B (en) Monocular vision distance measuring method
CN100557634C (en) A kind of camera marking method based on double 1-dimension drone
CN105806780A (en) Multi-vision-based system and method for measuring gas-liquid two-phase flow phase fraction of small channel
Wang et al. Liquid-level measurement using a single digital camera
CN109751964B (en) High-precision non-contact pipe diameter measuring method and device
CN107358628B (en) Linear array image processing method based on target
CN104596898A (en) Micro-bubble dimension online measurement apparatus and method thereof
CN104457703A (en) Multiple-vision-based System and method for measuring parameters of gas-liquid two-phase flow in small channel based on
Li et al. Factors influencing the accuracy of the photogrammetry-based deformation measurement method
WO2024001467A1 (en) Fatigue crack propagation rate test method and device based on deep learning
Hall et al. Uncertainty characterization of particle location from refocused plenoptic images
CN109974618A (en) The overall calibration method of multisensor vision measurement system
CN109490251A (en) Underwater refractive index self-calibrating method based on light field multilayer refraction model
Xue et al. Error analysis and liquid film thickness measurement in gas–liquid annular flow
CN114722549B (en) Dynamic parameter inversion method, device and medium for leakage source of underwater gas pipeline
CN105426830A (en) Small-channel gas liquid two-phase flow pattern identification system and method based on multi-vision information fusion technology
Zhao et al. Metric calibration of unfocused plenoptic cameras for three-dimensional shape measurement
CN110378884A (en) The method that Visual retrieval video carries out eccentric correction and expansion transformation in managing
CN110702015B (en) Method and device for measuring icing thickness of power transmission line
CN107256563A (en) Underwater 3 D reconstructing system and its method based on difference liquid level image sequence
KR100438212B1 (en) Extraction method for 3-dimensional spacial data with electron microscope and apparatus thereof
Du et al. Study on the reconstruction method of stereo vision in glass flume
CN109297554B (en) Method for measuring liquid phase flow in T-shaped pipe
Maas Contributions of digital photogrammetry to 3-D PTV

Legal Events

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