CN1270054C - Device and method for microscopic quantitatively processing dynamic color image for oil drive - Google Patents

Abstract

The present invention relates to a device and a method for quantitative processing of a microscopic oil-displacing dynamic colored image. The method can be used for carrying out various types of oil-displacing experimental processes on a microscopic network model (visible), accurately describing the continuously changed motion trajectories and the displacement laws of various displacing phases and rapidly and accurately calculating the quantity of residual oil in a calculation model to determine the corresponding relationship between the injection quantity and the crude oil recovery degree so that the analysis of microscopic seepage flow experimental data is converted from qualitative description into quantitative calculation. System hardware is mainly composed of a micro-quantity displacement pump, a microscopic simulation model, a microscope, a color video camera, a color image frame storing board, a color image decoder, a multifunctional card, a computer, a monitor, a color ink jet printer, etc. System software is composed of dynamic color image processing software, static color image processing software and universal color image processing software. Complete Chinese character operation is adopted in the process of image processing, and a chart automatically generated from obtained data and output can also be stored for secondary development and secondary editing in a database.

Description

Microcosmic oil drive dynamic color image quantizing method

Technical field

The invention belongs to oil field development experimental technique field, relate in particular to a kind of microcosmic oil drive dynamic color image quantification treatment device and processing method.

Background technology

Use the microcosmos network model and carry out the extensive use at home and abroad of indoor displacement of reservoir oil simulated experiment technology, but in the past to experimental data obtain and analyze only rest on videotape and means such as photograph on, can only make qualitative description.The analysis and processing method that obtains the quality and the data of improvement that how to improve experimental data is the problem that oil experimental work person is concerned about always and actively develops for many years, image analysis technology is the effective means that improves microscopic seepage experimental data analysis precision, by this technology, the analysis of experimental data is calculated to quantitative by qualitative description, make the variation of DATA PROCESSING generation matter.At present, the domestic effective ways that also are not used in the oil displacement experiment dynamic image analysis, part unit also can only carry out quantitative processing to static datas such as rock core casting body flakes from the instrument of external introduction.

Summary of the invention

In order to utilize computer quickly and accurately the microcosmic oil drive experimentation to be made quantitative Treatment, the invention provides a kind of microcosmic oil drive dynamic color image quantification treatment device and the processing method that is used for this device.

For achieving the above object, microcosmic oil drive dynamic color image quantification treatment device provided by the invention comprises:

The micro-displacement pump of displacement liquid;

The stratum simulation model of microscopic of experiment;

Amplify the microscope of micromodel;

The colored stylus of picked-up micromodel dynamic image can be installed on the microscope, and its signal can be imported video tape recorder and record, and also can import computer admission image;

Colour decoder is decomposed into R, G, B three-component with composite color video signal, so that frame is deposited three road A/D converter work on the plate;

The cromogram picture frame is deposited plate, can directly insert in the dead slot in the computer casing;

Colour picture monitor can be accepted the video input of multiple standard, deposits plate with frame and is connected;

Computer, and

Printer.

The microcosmic oil drive dynamic color image quantizing method that is used for above-mentioned treating apparatus provided by the invention comprises dynamic color image handling procedure and static chromatic image handling procedure two parts, is described below respectively:

1, dynamic color image handling procedure

(1) dynamic image admission mainly realizes the continuous admission of dynamic image, and record has related parameter when videotaping, as pattern number, experiment number, screen proportion, driving type, injection rate, multiplication factor etc.When videotaping on computers, the admission time of an images is 1/25 second, and the time interval of admission image is 1.5 seconds continuously, can manually enroll when enrolling image interval greater than 1.5 seconds, also can enroll automatically, and the time interval of admission can be set arbitrarily automatically.These picture intelligence losses are very little, and the resolution ratio height is the analysis and research displacement characteristics and measure the basis that related parameter is arranged.

(2) dynamic image shows, playback and the jigsaw puzzle function of dynamic image on monitor that main realization has been enrolled shows that simultaneously image has related parameter.This helps to choose the image that needs are handled, and measures required characteristic parameter.

(3) dynamic image is safeguarded, realizes the insertion and the deletion of dynamic image, to guarantee the continuity of displacement process changing features image.

(4) Tu Xiang compression and dump, the dynamic image of whole model is compressed to original 1/10～1/20 so that store on the floppy disk.

(5) the hole parameter is analyzed, and calculates the characteristic parameter of micromodel hole, venturi, and preserves and cut apart back skeleton template image GJ.JRP, uses for later analysis.

(6) recovery ratio tracing analysis utilizes the dynamic image of enrolling in the displacement process to measure the related parameter that has of displacing phase and displaced phase, calculates the recovery percent of reserves under different driving type, the different injection doubly, and draws displacement multiple and recovery percent of reserves relation curve.

(7) the water costa is analyzed, and extracts water costa and several continuous water leading edge images of describing and superpose from every width of cloth dynamic image, calculates the fltting speed and the water logging feature of water all directions in whole model.

(8) macroscopical trajectory analysis, the flow trace of display-object on whole model drawn the rate curve of this target on all directions.

(9) microcosmic trajectory analysis is selected suitable multiplication factor, follows the tracks of the admission target at intrapore dynamic image, and the target signature in each images is measured, and obtains the variation characteristic and the movement locus of target.

(10) demarcate multiplication factor, be determined at the length of the scale under certain multiplication factor, pixel length and area are converted to actual measured value.

2, static chromatic image handling procedure

This program can be carried out the automatic searching and the extraction of hole identification and venturi to the image of filling colour liquid or solid in the pore constriction, determines by operating personnel whether the venturi that is extracted is correct.Hole, venturi identification are being calculated the parameter and the distribution thereof of hole, venturi after judging, and exported corresponding chart.

The present invention has following advantage:

1, present technique dynamic color image process software can absorb image selectively and it is deposited in the hard disk of computer, so that call at any time during quantitative analysis displacement rule, has realized all admissions of image of displacement process.

2, realized the mutual conversion of R, G, B (red, green, blue) color space and H, V, C (lightness, tone, saturation ratio) color space, made chromatic image really realize the RGB processing, handled than black and white and pseudo-colours image and improved certainty of measurement.

3, can carry out the mensuration of 17 basic parameters to the target in the image.

4, realized the tracking and the measurement of microcosmic target in the hole.

5, are the important breakthroughs in this field to the different accurate quantifications that inject recovery percent of reserves under the multiple in the microcosmic oil drive experimentation.

6, in image processing procedure, adopt full Chinese character operation, and the data that obtain can be generated automatically chart output, also can deposit database in, be convenient to the secondary development editor and use.

Also in conjunction with the accompanying drawings the present invention is described in detail with embodiment below.

Description of drawings

Fig. 1 is apparatus of the present invention schematic diagram;

Fig. 2 is a dynamic color image handling procedure flow chart of the present invention;

Fig. 3 is the static chromatic image handling procedure of the present invention flow chart;

Fig. 4 is the pore-size distribution histogram;

Fig. 5 is the aperture area distribution histogram;

Fig. 6 is the venturi distribution histogram;

Fig. 7 is water drive-ternary composite driving recovery ratio curve;

Fig. 8 for polymer flooding for residual oil after the water drive;

Fig. 9 foam compound system displacement cecum residual oil;

Figure 10 is the different stacking charts of water costa constantly;

Figure 11 is that the water costa is at difference moment hodograph;

Figure 12 is microcosmic target (oil droplet) movement locus figure.

The specific embodiment

The device that present embodiment adopts is respectively:

Trace displacement pump 1, ISCO-100DX type, range of flow are 0.01 μ l/min～25ml/min, pressure limit 0～70Mpa, stability of flow can satisfy the requirement of various analog rates.

Simulation model of microscopic 2, use technology such as photographs, plate-making, the pore structure on the rock core casting body flake is etched on the flat glass plate after amplification, make the flat glass model through methods such as bonding, sintering, size is 4cm * 4cm, and its pore structure is similar to the reservoir rock pore structure.

Micro-border 3, object lens multiple 1 *, 4 * and 10 *, can make whole model is an images, can make an images only comprise several holes again, even only be the part of a hole, thereby make microcosmic and macroscopical target all enroll image and measured accurately.

Stylus 4, the WV-CP410/G of Panasonic type, videotaping sensitivity is 3LUX, plain 815 points of horizontal picture dot, resolution ratio is 600 row, can be installed on the micro-border, its signal can be imported video tape recorder and record, and also can import computer admission image.

Colour decoder 5, the VCB1 type is decomposed into R, G, B three-component with composite color video signal, so that frame is deposited three road A/D converter work on the plate.

The cromogram picture frame is deposited plate 6, the MVC32P type can directly insert in the dead slot in the computer casing, and four 512 * 512 * 8bit frames are deposited the true colour imagery capture card of body, image acquisition speed is every frame 1/25 second, but color picture signal can be simultaneously also timesharing be input to the RGB frame and deposit in the body.

Colour picture monitor 7, the PVM-1440QM of Sony can accept the video input of multiple standard.

Computer 8, association's Benyue 1100 types, 128MB internal memory, 12GB hard disk.

Printer 9, EPSON Stylus Color color inkjet printer.

More than device is connected to technique known.

The program that is used for said apparatus comprises dynamic color image handling procedure and static chromatic image handling procedure two partly, is described below respectively:

Dynamic color image handling procedure flow process:

1) picked-up is visual, absorbs whole images of required experimentation, is stored in the hard disk;

2) select analysis project, select the evaluating objects parameter;

The tracing analysis of A1 recovery ratio, the project of selecting analysis are the recovery ratios in the displacement test process;

A2 chooses the saturated oils image, the image when choosing model saturated oils before the not displacement;

A3 calculates the saturated oils area, calculates the area (model hole thickness is identical, and the amount of area can be considered as volume) of saturated oils in the initial mask;

A4 chooses the oil displacement process image, chooses the images after displacing fluid has entered the model hole, (by the priority sequencing of time);

A5 calculates the area of oil in the oil displacement process, calculates the area of oil in the selected image;

A6 calculates this process recovery ratio, and the area (volume) of selected image oil and the area (volume) of saturated oils are compared, and calculates by the fuel-displaced ratio of displacement;

A7 draws the recovery ratio curve, provides the recovery ratio form, provides the recovery ratio tables of data of whole each process of experiment, draws the recovery ratio curve of whole oil displacement process; (Fig. 7, table 4)

The B1 target trajectory is analyzed, and analyzes macroscopic view, movement locus and the kinematic parameter of microcosmic target in a part of hole of whole model or model;

B2 chooses initial aperture position image, chooses the image of the model hole at place in the target travel process;

B3 makes the pore structure template image, with in macroscopic view or the microcosmic target travel process the hole of process handle, make the negative of these holes, with reference as position and motion mode in the target travel process;

B4 chooses the target travel process image, chooses the images (by the priority sequencing of time) in the target travel process;

B5 carries out target location and parameter analysis, calculates target in whole model or the movement velocity in the part hole, shape and location parameter;

B6 render target movement locus provides the position and the kinematic parameter of target, and the resulting result of B5 is depicted as chart (Figure 10, Figure 11, Figure 12, table 5, table 6);

1) result preserves printout, and resulting chart can be deposited, also printable output.

Above-mentioned flow process as shown in Figure 2.

Static chromatic image handling procedure flow process:

1) picked-up image or read image, picked-up casting body flake or model hole image also can read casting body flake or the model hole image preserved;

2) defined analysis window is chosen the hole of the analysis part of wanting;

3) carry out image and cut apart, the difference according to hole and skeleton color splits hole;

4) venturi extracts automatically, is identified by big young pathbreaker's venturi of computer according to hole and venturi width;

5) the manual modification, the venturi that can't extract at computer can manually be clicked the mouse and be cut apart;

6) calculated hole diameters parameter, calculate pore diameter size, distribution and statistical parameter (table 1, Fig. 4), pore area size and distribute (table 2, Fig. 5) and throat size and distribution (table 3, Fig. 6);

7) next width of cloth that adds up can add up to several visual result of calculation;

8) result preserves printout, and result of calculation can be deposited, and is also exportable.

Above-mentioned flow process as shown in Figure 3.

The result of present embodiment is as follows:

1, hole, larynx parameter are analyzed

Pore throat parameter to practice ground layer emulation micromodel is carried out statistical analysis, can obtain:

Hole parameter and pore-size distribution-face rate, maximum coordination number, minimum ligancy, average coordination number, average pore diameter, sorting coefficient, the degree of bias, sharp degree, the coefficient of variation, average pore diameter ratio, pore-size distribution (table 1) and pore-size distribution histogram are seen Fig. 4.

Table 1 hole parameter and pore size distribution data table

Face rate (%)	Minimum ligancy	Maximum coordination number	Average coordination number	Average pore diameter (μ m)
					40.35	1	10	3.719	1623.32
Sorting coefficient	The degree of bias	The point degree	The coefficient of variation	Average pore throat diameter ratio
					0.58	1.82	5.62	0.95	6.99

Classification	Φ	The aperture	Number of apertures	Frequency (%)	Classification	Φ	Aperture (μ m)	Number of apertures	Frequency (%)
										1	＞3.32	＜100.00	0	0.00	11	-0.14	1100.00	3	5.26
2	2.32	200.00	0	0.00	12	-0.26	1200.00	0	0.00
										3	1.74	300.00	0	0.00	13	-0.38	1300.00	2	3.51
4	1.32	400.00	1	1.75	14	-0.49	1400.00	0	0.00
										5	1.00	500.00	0	0.00	15	-0.580	1500.00	3	5.26
6	0.74	600.00	2	3.51	16	-0.68	1600.00	3	5.26
										7	0.51	700.00	0	0.00	17	-0.77	1700.00	3	5.26
8	0.32	800.00	4	7.02	18	-0.85	1800.00	3	5.26
										9	0.15	900.00	1	1.75	19	-0.93	1900.00	0	0.00
10	0.00	1000.00	0	0.00	20	＜-0.93	＞1900.00	32	56.14

Simultaneously can provide pairing number of apertures of different aperture area (table 2) and corresponding histogram, see Fig. 5.

Table 2 pore area distributed data table

Classification	Classification area (μ m 2)	Hole area (μ m 2)	Number of apertures	Ratio (%)	Classification	Classification area (μ m 2)	Hole area (μ m 2)	Number of apertures	Ratio (%)
										1	＜7.85e+03	0.00e+00	0	0.00	11	1.58e+06	6.09e+06	4	5.14
2	1.65e+05	4.39e+05	3	0.37	12	1.74e+06	1.70e+06	1	1.43
										3	3.22e+05	7.46e+05	3	0.63	13	1.89e+06	9.16+06	5	7.73
4	4.79e+05	1.27e+06	3	1.07	14	2.05e+06	0.00e+00	0	0.00
										5	6.36e+05	4.90e+05	1	0.41	15	2.21e+06	6.34e+06	3	5.35
6	7.93e+05	4.21e+06	6	3.55	16	2.36e+06	2.25e+06	1	1.90

The venturi number (table 3) and the histogram of venturi parameter and venturi distribution-maximum throat diameter, minimum throat diameter, average throat diameter, different venturi size correspondences are seen Fig. 6.

Table 3 venturi distributed data table

Classification	Venturi width (μ m)	The venturi number	Frequency (%)	Classification	Venturi width (μ m)	The venturi number	Frequency (%)
								1	＜25.00	0	0.00	11	275.00	1	0.96
2	50.00	0	0.00	12	300.00	6	5.77
								3	75.00	1	0.96	13	325.00	1	0.96
4	100.00	14	13.46	14	350.00	5	4.81
								5	125.00	9	8.65	15	375.00	1	0.96
6	150.00	16	15.38	16	400.00	3	2.88
								7	175.00	7	6.73	17	425.00	1	0.96
8	200.00	12	11.54	18	450.00	1	0.96
								9	225.00	0	0.00	19	475.00	2	1.92
10	250.00	17	16.35	20	≥475.00	7	6.73

Maximum value=1137.68 μ m minimum value=58.59 μ m averages=232.19 μ m

2. inject multiple and recovery percent of reserves relation curve

Because it is the voids volume of micromodel has only tens microlitres, quite difficult to the metering of extraction liquid.For many years, seepage flow experimental work person can only adopt means such as microcosmic photograph and shooting that experimentation is made qualitative description, experimental result is quantized, and the needs that can not satisfy experiment (particularly tertiary oil recovery microcosmic oil drive experiment) research have hindered the raising of research level.This cover treatment system can be carried out quantitative analysis apace to the experiment dynamic image, determine the corresponding relation of injection rate and crude oil recovery percent of reserves, particularly can be effectively under the tertiary oil recovery different wetting condition, various displacement systems, different displacement program displacement of reservoir oil result compare (table 4), and then can make assay to oil displacement efficiency by the extraction curve, see 7.

Table 4 tertiary oil recovery microcosmic oil drive experimental result

Drive type	Breakthrough recovery (%)		Waterflood recovery efficiency factor (%)		Chemical flooding recovery ratio (%)
							Water is wet	Oil is wet	Water is wet	Oil is wet	Water is wet	Oil is wet
	Direct polymer flooding	56.3	10.0			56.5							48.1
Water drive-polymer flooding							30.5	11.0	31.4	37	51.0	48.7
	Direct ternary composite driving	77.0	63.0			81.9							76.3
Water drive-ternary is driven							25.2	10.8	29.4	38.5	97.3	75.4
	Water drive-polymer flooding-ternary is driven	30.5	11.0	31.4	37.0	97.5							72.0
Direct foam composite flooding							68.1	34.7			98.3	81.0
	Water drive-foam composite flooding	28.7	21.4	30.5	36.2	99.4							72.3
Water drive-polymer flooding-ternary is driven-foam flooding								11		37.0		76.4

3. microcosmic oil drive mechanism image

During the tertiary oil recovery different chemical agent displacement of reservoir oil, this quantization system can clearly be differentiated oil, water, chemical agent at any time to viewed ken admission image, and output apace, as Fig. 8, shown in Figure 9.

4. the displacing phase costa is analyzed

Utilize the abundant colour of computer to show, can provide the different costa stacking charts constantly of displacing phase, see Figure 10, and obtain seepage velocity value (table 5) and average seepage velocity value and the corresponding distribution map of displacing phase leading edge target, see Figure 11 in difference moment X, Y direction.

Table 5 water leading edge is at difference speed data table constantly

Sequence number	Time (time: divide: second)	Coordinate position	Speed (μ m/s)	Directions X speed (μ m/s)	Y direction speed (μ m/s)
						1	0:6:13	363，384	20.77	14.68	14.68
2	0:6:18	266，276	2409.38	1609.96	1792.53
						3	0:6:20	223，237	2408.78	1784.23	1618.26
4	0:6:22	158，146	4640.26	2697.10	3775.93
						5	0:6:27	82，63	1867.87	1261.41	1377.59
6	0:6:31	59，23	957.28	477.18	829.88
						7	0:6:32	10，9	4229.11	4066.39	1161.83

Average speed: 130.84 μ m/s

5. residual oil target analysis

In the displacement process, can carry out quantitative analysis to the residual oil target morphology, promptly difference constantly remaining oil droplet area, girth, maximum gauge, minimum diameter, equivalent circular radius, radius of curvature etc. (table 6) and can obtain participating in the movement locus of target, as shown in figure 12.

Table 6 Micro dynamic target signature survey data table

Sequence number	Time (every minute and second)	Area μ m 2)	Girth μ m)	Feret maximum gauge (μ m)	Feret minimum diameter (μ m)	Equivalent circular radius μ m)	Specific surface (μ m -1)	Form factor	Barycentric coodinates (X, Y)	Front end radius of curvature (μ m)	Hole is to oil droplet width (μ m)	By pearl width (μ m)	By pearl to pore width (μ m)
														1	00:00	51514	838	266	242	128	0.021	0.92	478，368	110	59	257	312
2	00:02	57417	869	281	260	135	0.019	0.96	402，324	120	120	260	160
														3	00:08	61367	900	289	266	140	0.019	0.95	347，241	120	314	272	203
4	00:09	55998	873	273	260	134	0.020	0.92	332，179	123	306	263	360
														5	00:10	58429	880	281	263	136	0.019	0.95	300，132	120	199	273	460
6	00:11	55126	869	273	255	132	0.020	0.92	268，116	133	166	278	454

Operation result shows that hardware capability of the present invention is complete, and running software is good, and parameters reaches the requirement of experimental precision, can satisfy the needs of microcosmic oil drive experiment.

Claims (1)

1, a kind of microcosmic oil drive dynamic color image quantizing method is made up of dynamic color image process software and static chromatic image process software two parts, and wherein, the dynamic color image process control method comprises:

1) picked-up is visual, absorbs whole images of required experimentation, is stored in the hard disk;

2) select analysis project, select the evaluating objects parameter;

The tracing analysis of A1 recovery ratio, the project of selecting analysis are the recovery ratios in the displacement test process;

A2 chooses the saturated oils image, the image when choosing model saturated oils before the not displacement;

A3 calculates the saturated oils area, calculates the area (model hole thickness is identical, and the amount of area can be considered as volume) of saturated oils in the initial mask;

A4 chooses the oil displacement process image, chooses a images after displacing fluid has entered the model hole by the sequencing of time;

A5 calculates the area of oil in the oil displacement process, calculates the area of oil in the selected image;

A6 calculates this process recovery ratio, and the area (volume) of selected image oil and the area (volume) of saturated oils are compared, and calculates by the fuel-displaced ratio of displacement;

A7 draws the recovery ratio curve, provides the recovery ratio form, provides the recovery ratio tables of data of whole each process of experiment, draws the recovery ratio curve of whole oil displacement process;

The B1 target trajectory is analyzed, and analyzes macroscopic view, movement locus and the kinematic parameter of microcosmic target in a part of hole of whole model or model;

B2 chooses initial aperture position image, chooses the image of the model hole at place in the target travel process;

B3 makes the pore structure template image, with in macroscopic view or the microcosmic target travel process the hole of process handle, make the negative of these holes, with reference as position and motion mode in the target travel process;

B4 chooses the target travel process image, and in chronological sequence order is chosen the images in the target travel process;

B5 carries out target location and parameter analysis, calculates target in whole model or the movement velocity in the part hole, shape and location parameter;

B6 render target movement locus provides the position and the kinematic parameter of target, and the resulting result of B5 is depicted as chart;

3) result preserves printout, and resulting chart can be deposited, also printable output;

Static chromatic image process software flow process:

1) picked-up image or read image, picked-up casting body flake or model hole image also can read casting body flake or the model hole image preserved;

2) defined analysis window is chosen the hole of the analysis part of wanting;

3) carry out image and cut apart, the difference according to hole and skeleton color splits hole;

4) venturi extracts automatically, is identified by big young pathbreaker's venturi of computer according to hole and venturi width;

5) the manual modification, the venturi that can't extract at computer can manually be clicked the mouse and be cut apart;

6) calculated hole diameters parameter calculates pore diameter size, distribution and statistical parameter, pore area size and distribution and throat size and distribution;

7) next width of cloth that adds up can add up to several visual result of calculation;

8) result preserves printout, and result of calculation can be deposited, and is also exportable.

Images