CN103206208B - A kind of method of microcosmic different existence state remaining oil macro-level quantitative - Google Patents
A kind of method of microcosmic different existence state remaining oil macro-level quantitative Download PDFInfo
- Publication number
- CN103206208B CN103206208B CN201310085454.9A CN201310085454A CN103206208B CN 103206208 B CN103206208 B CN 103206208B CN 201310085454 A CN201310085454 A CN 201310085454A CN 103206208 B CN103206208 B CN 103206208B
- Authority
- CN
- China
- Prior art keywords
- remaining oil
- oil
- remaining
- different
- existence state
- 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.)
- Active
Links
Abstract
The invention discloses a kind of method of microcosmic different existence state remaining oil macro-level quantitative, by macroscopical numerical simulation, obtain remaining oil reserve distribution field and water saturation distribution field; Determine the underground water saturation scope under different washing degree condition according to the washing degree criteria for classifying and original water saturation, and ask for different washing degree macroscopic view remaining oil reserves in conjunction with remaining oil reserve distribution field and water saturation distribution field; Experiment is corrected by getting core well water saturation, determine core washing degree under current state, by getting core well oil-containing sheet microscopic observation experiment, determine microcosmic different existence state remaining oil relative amount, and statistics draws the relative amount of different existence state microscopic remaining oil under different washing degree; According to the relative amount of different existence state microscopic remaining oil under different washing degree macroscopic view remaining oil reserves and different washing degree, application microcosmic different existence state Fuel Oil Remaining calculation expression carries out macro-level quantitative to microscopic remaining oil.
Description
Technical field
The invention belongs to technical field of petroleum extraction, particularly relate to a kind of method of microcosmic different existence state remaining oil macro-level quantitative.
Background technology
From the eighties in 20th century, research remaining oil distribution has caused the common concern of each Oil Producing Countries of the world, remaining oil number and occurrence status be the basic condition improving recovery ratio further.Current remaining oil research technology is mainly divided into two aspects: one be between reflection well, macroscopical remaining oil of interlayer remaining oil distribution; Two is the microscopic remaining oils reflecting remaining oil distribution in rock pore structure.Macroscopic view remaining oil can use numerical simulation technology to accomplish quantification, and microscopic remaining oil can obtain by getting the analysis of core well oil-containing sheet, but its analysis result is just based on the relative quantity of oil-containing sheet, the percentage namely shared by different existence state remaining oil.As everyone knows, macroscopical remaining oil is made up of the microscopic remaining oil of different existence state, and when macroscopical Fuel Oil Remaining one timing, wherein comprised different existence state microscopic remaining oil absolute magnitude will be the key determining oilfield well network adjustment, chemical flooding agent design.Because current remaining oil research technology method cannot know the absolute magnitude of different existence state remaining oil, therefore, usually use the method for laboratory core displacement test when driving agent design, be similar to " black box " theoretical, result is counter by experiment pushes away the chemical system mated with oil reservoir.
The difference of the research method of current microscopic remaining oil experimentally equipment can be divided into two aspects: the first kind is the relative amount obtaining different existence state remaining oil based on two-dimentional core wafer analysis, and research method comprises casting body flake, ESEM, fluorescence microscope, glass etching oil displacement experiment etc.; Equations of The Second Kind is the relative quantity and the absolute magnitude that are obtained different existence state remaining oil after successively scanning based on the little rock core of three-dimensional by image processing techniques, and research method comprises nanometer CT technology, nuclear magnetic resonance.
Based on the microscopic remaining oil research that two-dimentional core wafer is analyzed, be merely able to the relative amount obtaining different existence state remaining oil, and research precision is not high; Microscopic remaining oil based on the little core analysis of three-dimensional is studied, although can at three-dimensional perspective study of rocks mineral and organic matter, and pass through relative quantity and the absolute magnitude that the image processing techniques after successively scanning obtains different existence state remaining oil, but, the core sample that can scan due to CT or nuclear-magnetism and small, volume only has grade to arrive Centimeter Level size, and laboratory apparatus is expensive, therefore, the method as theoretical research, can only not be suitable for produced on-site practice at present.
Microscopic remaining oil is the remaining oil be trapped in different existence state (membranaceous, tufted, corner shape etc.) in rock pore structure, because current research method is all from microcosmic angle, mutually isolate with macroscopical remaining oil, do not set up dependency relation with it, therefore, all can only obtain the absolute magnitude in the relative quantity of different existence state microscopic remaining oil or tiny sampler, also just cannot characterize different existence state microscopic remaining oil absolute magnitude size in actual oil reservoir, and this absolute magnitude is the key determining oilfield well network adjustment, chemical flooding agent design.
Summary of the invention
The invention provides a kind of method of microcosmic different existence state remaining oil macro-level quantitative, be intended to solve the absolute magnitude in the relative quantity or tiny sampler that can only obtain different existence state microscopic remaining oil at present, the problem of different existence state microscopic remaining oil absolute magnitude size in actual oil reservoir cannot be characterized.
The object of the present invention is to provide a kind of method of microcosmic different existence state remaining oil macro-level quantitative, the method comprises the following steps:
Step one, by macroscopical numerical simulation, obtains remaining oil reserve distribution field and water saturation distribution field;
Step 2, determines the underground water saturation scope under different washing degree condition according to the washing degree criteria for classifying and original water saturation, and asks for different washing degree macroscopic view remaining oil reserves in conjunction with remaining oil reserve distribution field and water saturation distribution field;
Step 3, experiment is corrected by getting core well water saturation, determine core washing degree under current state, by getting core well oil-containing sheet microscopic observation experiment, determine microcosmic different existence state remaining oil relative amount, and statistics draws the relative amount of different existence state microscopic remaining oil under different washing degree;
Step 4, the microcosmic different existence state Fuel Oil Remaining calculation expression that to determine with different existence state microscopic remaining oil relative amount under different washing degree macroscopic view remaining oil reserves and different washing degree be function;
Step 5, the different washing degrees macroscopic view remaining oil reserves of trying to achieve according to step 2 and step 3 add up the relative amount of different existence state microscopic remaining oil under the different washing degrees drawn, the microcosmic different existence state Fuel Oil Remaining calculation expression of applying step four carries out macro-level quantitative to microscopic remaining oil.
Further, in step, by macroscopical numerical simulation, when obtaining remaining oil reserve distribution field and water saturation distribution field, on Reservoir Description basis, application CMG reservoir numerical simulation software carries out semi-automatic history matching for production to target reservoir, and then obtains macroscopical remaining oil reserve distribution field and macroscopical water saturation distribution field data volume, and concrete semi-automatic processing procedure is as follows:
The first step, determines adjustable reservoir parameter, comprising: permeability saturation curve (RP), the permeability Kx in X, Y, Z tri-directions, Ky, Kz, the conductivity Tx in X, Y, Z tri-directions, Ty, Tz, totally 7 parameters;
Second step, to each parameter weight, weight span is [0-100], by weight normalization, weighted value after normalization is the adjusting range of this parameter, parameter after adjustment is applied in field data calculating, dopes the production fluid of current oil reservoir, oil productive capacity and pressure maintenance level;
3rd step, keep level three parameters as the foundation of weight allocation next time using production fluid, oil productive capacity and pressure, if prediction liquid-producing capacity and pressure keep level lower than actual conditions, then permeability saturation curve (RP) and Kx, Ky distribute larger weight, if oil productive capacity is lower than actual conditions, permeability saturation curve (RP) distributes larger weight, repetition second step calculates, if three indexs of prediction conform to substantially with actual conditions, then terminate weight allocation, stop calculating.
Further, during adjustment permeability saturation curve (RP), if oil productive capacity is lower than actual conditions, then the adjustment of oil relative permeability forward increases, and aqueous phase relative permeability oppositely adjusts reduction; If oil productive capacity is higher than actual conditions, then adjustment direction is just the opposite;
The adjusting range of permeability is limited, after permeability adjusting range is more than 50%, auto-steering is adjusted conductivity;
Weight is arranged cannot revise permeability saturation curve form, needs finger prosthesis to determine permeability saturation curve form.
Further, the method imports water saturation and the remaining oil reserves parameter field of any time by residual oil quantitative analytical system, is grade classification index, calculates the remaining oil reserves in different demarcation interval with water saturation.
Further, described residual oil quantitative analytical system comprises:
Data acquisition module, for receiving input data volume file, adopting sequential organization algorithm for design to realize, judging wrong position and reason in file input process, the relevant entry in manual modification input data volume file, and again reads;
Residual oil quantitative analysis module, is connected with described data acquisition module, for receiving the data of described data collecting module collected, calculating the remaining oil reserves in each interval respectively, and form standard output form with given criteria for classifying interval;
Visualized data output module, is connected with described residual oil quantitative analysis module, exports with form and EXECL standard pie chart form for the standard output form remaining oil reserves data volume that will generate.
Further, described data acquisition module and residual oil quantitative analysis module all adopt sequential organization algorithm for design to realize.
Further, microscopic remaining oil occurrence status adopts the experiment of sealed coring well oil-containing sheet laser confocal scanning to obtain, and detailed process is:
First, adopt freezing tabletting technology, carry out cutting sample at low temperature environment, when guaranteeing abrasive disc, original form of hole inner fluid is not destroyed;
Secondly, body of casting automatic image-analysis microscope, petrographic microscope and fluorescence microscope are combined, body of casting automatic image-analysis microscope is utilized to extract pore throat characteristic parameter, rock particles feature, the character of petrographic microscope identification mineral, fluorescence microscopy chooses Ultraluminescence filter, distinguish oils water boundary, the extraction of remaining oil occurrence status information is completed by residual oil quantitative analytical system, by observing fluoroscopic image, the remaining oil distribution state in hole being differentiated, characterizing the feature of remaining oil micro Distribution with this result.
The method of microcosmic different existence state remaining oil macro-level quantitative provided by the invention, by macroscopical numerical simulation, obtains remaining oil reserve distribution field and water saturation distribution field; Determine the underground water saturation scope under different washing degree condition according to the washing degree criteria for classifying and original water saturation, and ask for different washing degree macroscopic view remaining oil reserves in conjunction with remaining oil reserve distribution field and water saturation distribution field; Experiment is corrected by getting core well water saturation, determine core washing degree under current state, by getting core well oil-containing sheet microscopic observation experiment, determine microcosmic different existence state remaining oil relative amount, and statistics draws the relative amount of different existence state microscopic remaining oil under different washing degree; According to the relative amount of different existence state microscopic remaining oil under the different washing degrees that the different washing degree macroscopic view remaining oil reserves of trying to achieve and statistics draw, application microcosmic different existence state Fuel Oil Remaining calculation expression carries out macro-level quantitative to microscopic remaining oil; Microscopic remaining oil oil-containing sheet analysis result based on laser confocal scanning organically combines with numerical simulation macroscopic view remaining oil research technology result by the present invention, achieve the macro-level quantitative of different existence state microscopic remaining oil, precision is high, cost is low, meet the actual needs of oil field development, improve the research level of remaining oil, from the progressive accuracy research for quantification of conceptual research of non-quantitation, there is stronger propagation and employment and be worth.
Accompanying drawing explanation
Fig. 1 is the realization flow figure of the method for the microcosmic different existence state remaining oil macro-level quantitative that the embodiment of the present invention provides;
Fig. 2 is the distribution characteristics schematic diagram of the remaining oil reserve distribution field that the embodiment of the present invention provides;
Fig. 3 is the distribution characteristics schematic diagram of the water saturation distribution field that the embodiment of the present invention provides;
Fig. 4 is the main interface of the residual oil quantitative analytical system that the embodiment of the present invention provides;
Fig. 5 is the realization flow figure of data acquisition module in the residual oil quantitative analytical system that provides of the embodiment of the present invention;
Fig. 6 is the realization flow figure of residual oil quantitative analysis module in the residual oil quantitative analytical system that provides of the embodiment of the present invention;
Fig. 7 is the realization flow figure of visualized data output module in the residual oil quantitative analytical system that provides of the embodiment of the present invention;
Fig. 8 is the target reservoir different layers position different washing degree macroscopic view remaining oil reserves proportion distribution map that the embodiment of the present invention provides;
Fig. 9 is the operating principle flow chart of the residual oil quantitative analytical system that the embodiment of the present invention provides;
Figure 10 is the remaining oil microscopic fluorescence original image figure that the embodiment of the present invention provides;
Figure 11 is the rear image graph of employing threshold method segmentation that the embodiment of the present invention provides;
Figure 12 is the remaining oil type identification result figure that the embodiment of the present invention provides, and wherein, green portion is water, pink colour part is dissimilar remaining oil;
Figure 13 is all types of remaining oil microscopic pattern schematic diagram after two rivers oilfield water flooding of providing of the embodiment of the present invention;
Figure 14 is the different existence state microscopic remaining oil reserve distribution figure that the embodiment of the present invention provides.
Detailed description of the invention
In order to make object of the present invention, technical scheme and advantage clearly understand, below in conjunction with drawings and Examples, the present invention is described in further detail.Should be appreciated that specific embodiment described herein only in order to explain the present invention, and be not used in restriction invention.
Fig. 1 shows the realization flow of the method for the microcosmic different existence state remaining oil macro-level quantitative that the embodiment of the present invention provides.
The method comprises the following steps:
Step one, by macroscopical numerical simulation, obtains remaining oil reserve distribution field and water saturation distribution field;
Step 2, determines the underground water saturation scope under different washing degree condition according to the washing degree criteria for classifying and original water saturation, and asks for different washing degree macroscopic view remaining oil reserves in conjunction with remaining oil reserve distribution field and water saturation distribution field;
Step 3, experiment is corrected by getting core well water saturation, determine core washing degree under current state, by getting core well oil-containing sheet microscopic observation experiment, determine microcosmic different existence state remaining oil relative amount, and statistics draws the relative amount of different existence state microscopic remaining oil under different washing degree;
Step 4, the microcosmic different existence state Fuel Oil Remaining calculation expression that to determine with different existence state microscopic remaining oil relative amount under different washing degree macroscopic view remaining oil reserves and different washing degree be function;
Step 5, the different washing degrees macroscopic view remaining oil reserves of trying to achieve according to step 2 and step 3 add up the relative amount of different existence state microscopic remaining oil under the different washing degrees drawn, the microcosmic different existence state Fuel Oil Remaining calculation expression of applying step four carries out macro-level quantitative to microscopic remaining oil.
In embodiments of the present invention, in step, by macroscopical numerical simulation, when obtaining remaining oil reserve distribution field and water saturation distribution field, on Reservoir Description basis, application CMG reservoir numerical simulation software carries out semi-automatic history matching for production to target reservoir, and then obtains macroscopical remaining oil reserve distribution field and macroscopical water saturation distribution field data volume, and concrete semi-automatic processing procedure is as follows:
The first step, determines adjustable reservoir parameter, comprising: permeability saturation curve (RP), the permeability Kx in X, Y, Z tri-directions, Ky, Kz, the conductivity Tx in X, Y, Z tri-directions, Ty, Tz, totally 7 parameters;
Second step, to each parameter weight, weight span is [0-100], by weight normalization, weighted value after normalization is the adjusting range of this parameter, parameter after adjustment is applied in field data calculating, dopes the production fluid of current oil reservoir, oil productive capacity and pressure maintenance level;
3rd step, keep level three parameters as the foundation of weight allocation next time using production fluid, oil productive capacity and pressure, if prediction liquid-producing capacity and pressure keep level lower than actual conditions, then permeability saturation curve (RP) and Kx, Ky distribute larger weight, if oil productive capacity is lower than actual conditions, permeability saturation curve (RP) distributes larger weight, repetition second step calculates, if three indexs of prediction conform to substantially with actual conditions, then terminate weight allocation, stop calculating.
In embodiments of the present invention, during adjustment permeability saturation curve (RP), if oil productive capacity is lower than actual conditions, then the adjustment of oil relative permeability forward increases, and aqueous phase relative permeability oppositely adjusts reduction; If oil productive capacity is higher than actual conditions, then adjustment direction is just the opposite;
The adjusting range of permeability is limited, after permeability adjusting range is more than 50%, auto-steering is adjusted conductivity;
Weight is arranged cannot revise permeability saturation curve form, needs finger prosthesis to determine permeability saturation curve form.
In embodiments of the present invention, the method imports water saturation and the remaining oil reserves parameter field of any time by residual oil quantitative analytical system, is grade classification index, calculates the remaining oil reserves in different demarcation interval with water saturation.
In embodiments of the present invention, residual oil quantitative analytical system comprises:
Data acquisition module, for receiving input data volume file, adopting sequential organization algorithm for design to realize, judging wrong position and reason in file input process, the relevant entry in manual modification input data volume file, and again reads;
Residual oil quantitative analysis module, is connected with data acquisition module, for receiving the data of data collecting module collected, calculating the remaining oil reserves in each interval respectively, and form standard output form with given criteria for classifying interval;
Visualized data output module, is connected with residual oil quantitative analysis module, exports with form and EXECL standard pie chart form for the standard output form remaining oil reserves data volume that will generate.
In embodiments of the present invention, data acquisition module and residual oil quantitative analysis module all adopt sequential organization algorithm for design to realize.
In embodiments of the present invention, microscopic remaining oil occurrence status adopts the experiment of sealed coring well oil-containing sheet laser confocal scanning to obtain, and detailed process is:
First, adopt freezing tabletting technology, carry out cutting sample at low temperature environment, when guaranteeing abrasive disc, original form of hole inner fluid is not destroyed;
Secondly, body of casting automatic image-analysis microscope, petrographic microscope and fluorescence microscope are combined, body of casting automatic image-analysis microscope is utilized to extract pore throat characteristic parameter, rock particles feature, the character of petrographic microscope identification mineral, fluorescence microscopy chooses Ultraluminescence filter, distinguish oils water boundary, the extraction of remaining oil occurrence status information is completed by residual oil quantitative analytical system, by observing fluoroscopic image, the remaining oil distribution state in hole being differentiated, characterizing the feature of remaining oil micro Distribution with this result.
Below in conjunction with drawings and the specific embodiments, application principle of the present invention is further described.
The technic relization scheme of microcosmic different existence state remaining oil macro-level quantitative is as follows:
(1) first carry out macroscopical numerical simulation study, obtain remaining oil reserve distribution field and water saturation distribution field;
(2) determine the underground water saturation scope under different washing degree condition according to the washing degree criteria for classifying and original water saturation, and ask for different washing degree macroscopic view remaining oil reserves in conjunction with remaining oil reserve distribution field and water saturation distribution field;
(3) carry out get core well water saturation correct experiment, determine core washing degree under current state, carry out and get core well oil-containing sheet microscopic observation experiment, determine microcosmic different existence state remaining oil relative amount, and statistics draws the relative amount of different existence state microscopic remaining oil under different washing degree;
(4) the microcosmic different existence state Fuel Oil Remaining calculation expression that to propose with different existence state microscopic remaining oil relative amount under different washing degree macroscopic view remaining oil reserves and different washing degree be function;
(5) the different washing degrees macroscopic view remaining oil reserves of trying to achieve according to (2) step and (3) step add ups the relative amount of different existence state microscopic remaining oil under the different washing degrees drawn, the microcosmic different existence state Fuel Oil Remaining calculation expression applying (4) step quantizes microscopic remaining oil.
Fig. 1 is the macro-level quantitative method realization flow figure of microcosmic different existence state remaining oil.Flow chart shows, the present invention is take washing degree as the microcosmic different existence state residual oil calculation formula of function is bridge, macroscopical remaining oil research technology result and microscopic remaining oil result of study is combined, obtains microscopic remaining oil quantized result.Below for Henan Oil Field two river VI oil group, elaborate implementation of the present invention.
The acquisition of macroscopic field data:
On Reservoir Description basis, application CMG reservoir numerical simulation software carries out history matching for production to target reservoir, and then obtains macroscopical remaining oil reserves field and macroscopical water saturation field data body.The acquisition of field data of the present invention is different from traditional artificial treatment mode, but the processing mode of semi-automation, processing procedure is as follows:
(1) determine adjustable reservoir parameter, comprise, permeability saturation curve (RP), the permeability Kx in X, Y, Z tri-directions, Ky, Kz, the conductivity Tx in X, Y, Z tri-directions, Ty, Tz, totally 7 parameters;
(2) to each parameter weight, weight span is [0-100], by weight normalization, weighted value after normalization is the adjusting range of this parameter, in the field data of the parameter application after adjustment being calculated, the production fluid of current oil reservoir, oil productive capacity and pressure maintenance level etc. can be doped;
(3) keep level three parameters as the foundation of weight allocation next time using production fluid, oil productive capacity and pressure, if prediction liquid-producing capacity and pressure keep level lower than actual conditions, then RP and Kx, Ky distribute larger weight, if oil productive capacity is lower than actual conditions, RP distributes larger weight, repeat (2) step to calculate, if three indexs of prediction conform to substantially with actual conditions, terminate weight allocation, stop calculating.
Present treatment process need illustrates 3 points, one is the adjustment more complicated of permeability saturation curve, if when oil productive capacity is lower than actual conditions, then oil relative permeability forward adjustment (increase), aqueous phase relative permeability oppositely adjusts (reduction); If oil productive capacity is higher than actual conditions, then adjustment direction is just the opposite.Two are adjusting ranges of permeability is limited, therefore, after permeability adjusting range is more than 50%, auto-steering is adjusted conductivity.Three be weight arrange be to revise permeability saturation curve form, therefore need manually to specify its form, therefore the method is referred to as semi-automatic processing mode.
Fig. 2, Fig. 3 are respectively the reserves and water saturation field distribution feature that adopt above-mentioned processing mode to obtain.
Different washing degree macroscopic view remaining oil reserves: because CMG numerical simulation software cannot be classified quantitative analysis, for this reason, have developed residual oil quantitative analytical system software NewStoneOILANSv1.0.Software can realize importing the water saturation of any time and remaining oil reserves parameter field, and be then grade classification index with water saturation, calculate the remaining oil reserves in different demarcation interval, Fig. 4 is the main interface of residual oil quantitative analytical system software.
Residual oil quantitative analytical system is divided into three modules: data acquisition module, residual oil quantitative analysis module, visualized data output module.
Data acquisition module: this module is mainly used to receive input data volume file, sequential organization algorithm for design is adopted to realize, program realization flow figure is shown in Fig. 5, judging in file input process after wrong position and reason can relevant entry in amendment input data volume file manually, and again reads.
Residual oil quantitative analysis module: this module receives the data that data collecting module collected arrives, interval with the given criteria for classifying, calculate the remaining oil reserves in each interval respectively, and form standard output form.Adopt sequential organization algorithm for design to realize, program realization flow figure is shown in Fig. 6.
Visualized data output module: the standard output form remaining oil reserves data volume generated exports with form and EXECL standard pie chart form by this module, and program realization flow figure is shown in Fig. 7.
According to the oil reservoir water saturation value (0.52,0.68) of the weak washing determined, strong washing correspondence, application OILANS software calculates different washing degree macroscopic view remaining oil reserves form (table 1) in target reservoir different layers position, and reserves proportion distribution map (Fig. 8).
Table 1 different washing degree macroscopic view remaining oil reserves
Microscopic remaining oil occurrence status:
The microcosmic occurrence status of remaining oil adopts the experiment of sealed coring well oil-containing sheet laser confocal scanning to obtain.The general thought of experimental technique is: first, adopts freezing tabletting technology, carries out cutting sample at low temperature environment, and when guaranteeing abrasive disc, original form of hole inner fluid can not be destroyed; Secondly, that body of casting automatic image-analysis microscope, petrographic microscope and fluorescence microscope are combined, the advantage of comprehensive three, body of casting automatic image-analysis microscope is utilized to extract pore throat characteristic parameter, rock particles feature, the character of petrographic microscope identification mineral, fluorescence microscopy chooses Ultraluminescence filter, distinguish oils water boundary, the extraction of remaining oil occurrence status information is completed with remaining oil analysis software, by observing fluoroscopic image, the remaining oil distribution state in hole being differentiated, characterizing the feature of remaining oil micro Distribution with this result.
Analysis process:
Image capture: select high-sensitive colourful CCD video camera to gather R, G, B data image signal.
Color transformed: due to R, G, B tri-the colour space be linear, namely the equivalent R of zones of different, G, B are that brightness, colourity, saturation ratio are all different to the visual experience of people, be unfavorable for Iamge Segmentation, therefore set up R, G, B and H, Mathematical Modeling between V, C (brightness, colourity, saturation ratio), realize R, G, B to H, V, C conversion.
Color enhancement: adopt involution enhancing and histogram equalization.Select the valley point of R, G, B curve as coefficient, carry out space and do involution enhancing.
Color classification: utilize the value of R, G, B and the peak value of H curve to carry out color classification.
Iamge Segmentation: by color of image by Objective extraction out.
Feature extraction: select suitable gray scale thresholding target to be separated from background by grey level histogram, and calculate oil area, containing parameters such as water area, remaining oil type ratios.
Iamge Segmentation: for remaining oil microscopic fluorescence image (as shown in Figure 10), filemot part is crude oil, blue part seawater, the task of Iamge Segmentation will extract crude oil (yellowish-brown) and water (blueness) exactly from the background of rock (dark color), but carries out follow-up analytical work again.
Threshold method is a kind of traditional image partition method, image (as shown in figure 11) after segmentation, and purple is crude oil, and green is water, and the algorithm that the present invention adopts is global threshold partitioning algorithm and auto-thresholding algorithm.
Individual identification: after above-mentioned Iamge Segmentation, can identify fuel-displaced and water from microscopic fluorescence image, detects the oil of all single types further from the foreground pixel of image.By adding up the size of remaining oil, shape, type and area, calculate the content of oil and water.But the shape of remaining oil and complex distribution in actual conditions, as shown in figure 12, the overlap of remaining oil continuous distributed or connection, need the edge line detecting all single remaining oils from image.Therefore, when carrying out remaining oil type distributes and analyzing, the higher algorithm of some automaticities is proposed to carry out the rim detection of single remaining oil.At present, Weak link labeling algorithm, corrosion propagation and principal curve cluster is had to combine algorithm, spot detection algorithm etc., but these algorithms the analysis of inapplicable remaining oil.Have employed the restraining growth algorithm having very much progress relative to algorithm above in automaticity and rim detection effect.
Feature extraction: after individual identification is carried out to microscopic remaining oil, independently remaining oil can be extracted.In the bianry image obtained, by statistical computation object pixel, obtain the parameters such as the size and dimension of remaining oil.Introduce the computational methods (by undefined computational process all in units of pixel) of Common Parameters in the analysis of remaining oil size and dimension below.
(1) oil area: the quantity of namely adding up the corresponding foreground pixel of this remaining oil area;
(2) containing water area: the quantity of namely adding up the corresponding foreground pixel of this water;
(3) girth: the quantity of namely adding up remaining oil marginal foreground pixel;
(4) remaining oil type.Namely the quantity of the corresponding foreground pixel of the type remaining oil is added up;
(5) oil-water ratio: oil area and the ratio containing the water surface;
(6) remaining oil type ratio: the ratio of the overall remaining oil in different remaining oil type place.
Reservoir microscopic remaining oil type after the oilfield water flooding of two rivers: the laser confocal scanning experiment having carried out 46 pieces of samples altogether, because reservoir micropore structure, exploitation water out behavior etc. there are differences, after waterflooding development, there is some difference for reservoir remaining oil distribution.Shuan He oil field VI oil group remaining oil type mainly comprises following several form (Figure 13) according to the observation: hole table film-form, granular absorption shape, tufted, intergranular adsorption-like, corner shape.
Different washing degree core microscopic remaining oil distribution feature:
Get core well saturation correction result to show, strong water logging during each layer position of oil reservoir belongs to substantially, macroscopically there is some difference for remaining oil distribution, and on microcosmic, also there is some difference for different washing degree microscopic remaining oil.The relative amount (table 2) of different washing degree different existence state microscopic remaining oil in 46 pieces of core samples has been added up according to above-mentioned remaining oil treatment technology, global regularity shows as: along with washing degree increases, free state and bound state remaining oil are in regular change, along with washing degree increases, the hole table film-form of bound state increases gradually, granular absorption shape is also substantially in increasing trend, and the tufted remaining oil of free state reduces gradually.
Table 2 different washing degree different existence state microscopic remaining oil relative amount
The macro-level quantitative result of microcosmic different existence state remaining oil
Obtained macroscopical absolute magnitude and the microcosmic relative quantity of remaining oil above by the method for method for numerical simulation and microcosmos experiment, in order to the two organically being combined, be defined as follows:
G
i; Microscopic remaining oil reserves, unit: 10
4t;
G
j; Macroscopic view remaining oil reserves, unit: 10
4t;
A
ij; Microscopic remaining oil relative amount, unit, %;
I: remaining oil microcosmic occurrence status;
J: washing degree type;
Therefore: for the oil reservoir having j kind washing degree, the reserves of microscopic remaining oil i in oil reservoir can be represented by the formula:
Two Variables is there is, wherein microscopic remaining oil relative amount a on the right side of above formula
ijcan be obtained after normalization by table 2, table 3 is a after normalization
ijbinary interaction factor table.
Table 3a
ijbinary interaction factor table
Macroscopic view remaining oil reserves G
jprovided (table 1) by first segment macroscopic view remaining oil research technology result, so far, just can ask for target reservoir different existence state microscopic remaining oil reserves (table 4).Contrast different existence state microscopic remaining oil reserves find (Figure 14), and because this oil reservoir washing degree is high, the remaining oil in oil reservoir shows film-form, corner shape, intergranular adsorption-like and granular absorption shape based on hole, only has a small amount of tufted remaining oil.Because the washing degree of different layers position is different, therefore, different existence state microscopic remaining oil proportion is also not quite similar.
Table 4 different existence state microscopic remaining oil reserves
The method of the microcosmic different existence state remaining oil macro-level quantitative that the embodiment of the present invention provides, by macroscopical numerical simulation, obtains remaining oil reserve distribution field and water saturation distribution field; Determine the underground water saturation scope under different washing degree condition according to the washing degree criteria for classifying and original water saturation, and ask for different washing degree macroscopic view remaining oil reserves in conjunction with remaining oil reserve distribution field and water saturation distribution field; Experiment is corrected by getting core well water saturation, determine core washing degree under current state, by getting core well oil-containing sheet microscopic observation experiment, determine microcosmic different existence state remaining oil relative amount, and statistics draws the relative amount of different existence state microscopic remaining oil under different washing degree; According to the relative amount of different existence state microscopic remaining oil under the different washing degrees that the different washing degree macroscopic view remaining oil reserves of trying to achieve and statistics draw, application microcosmic different existence state Fuel Oil Remaining calculation expression carries out macro-level quantitative to microscopic remaining oil; Microscopic remaining oil oil-containing sheet analysis result based on laser confocal scanning organically combines with numerical simulation macroscopic view remaining oil research technology result by the present invention, achieve the macro-level quantitative of different existence state microscopic remaining oil, precision is high, cost is low, meet the actual needs of oil field development, improve the research level of remaining oil, from the progressive accuracy research for quantification of conceptual research of non-quantitation, there is stronger propagation and employment and be worth.
These are only preferred embodiment of the present invention, not in order to limit the present invention, all any amendments done within the spirit and principles in the present invention, equivalent replacement and improvement etc., all should be included within protection scope of the present invention.
Claims (7)
1. a method for microcosmic different existence state remaining oil macro-level quantitative, it is characterized in that, the method comprises the following steps:
Step one, by macroscopical numerical simulation, obtains remaining oil reserve distribution field and water saturation distribution field;
Step 2, determines the underground water saturation scope under different washing degree condition according to the washing degree criteria for classifying and original water saturation, and asks for different washing degree macroscopic view remaining oil reserves in conjunction with remaining oil reserve distribution field and water saturation distribution field;
Step 3, experiment is corrected by getting core well water saturation, determine core washing degree under current state, by getting core well oil-containing sheet microscopic observation experiment, determine microcosmic different existence state remaining oil relative amount, and statistics draws the relative amount of different existence state microscopic remaining oil under different washing degree;
Step 4, the microcosmic different existence state Fuel Oil Remaining calculation expression that to determine with different existence state microscopic remaining oil relative amount under different washing degree macroscopic view remaining oil reserves and different washing degree be function;
Step 5, the different washing degrees macroscopic view remaining oil reserves of trying to achieve according to step 2 and step 3 add up the relative amount of different existence state microscopic remaining oil under the different washing degrees drawn, the microcosmic different existence state Fuel Oil Remaining calculation expression of applying step four carries out macro-level quantitative to microscopic remaining oil.
2. the method for claim 1, it is characterized in that, in step, by macroscopical numerical simulation, when obtaining remaining oil reserve distribution field and water saturation distribution field, be on Reservoir Description basis, application CMG reservoir numerical simulation software carries out semi-automatic history matching for production to target reservoir, and then obtain macroscopical remaining oil reserve distribution field and macroscopical water saturation distribution field data volume, concrete semi-automatic processing procedure is as follows:
The first step, determines adjustable reservoir parameter, comprising: relative permeability RP curve, the permeability Kx in X, Y, Z tri-directions, Ky, Kz, the conductivity Tx in X, Y, Z tri-directions, Ty, Tz, totally 7 parameters;
Second step, to each parameter weight, weight span is [0-100], by weight normalization, weighted value after normalization is the adjusting range of this parameter, parameter after adjustment is applied in field data calculating, dopes the production fluid of current oil reservoir, oil productive capacity and pressure maintenance level;
3rd step, keep level three parameters as the foundation of weight allocation next time using production fluid, oil productive capacity and pressure, if prediction liquid-producing capacity and pressure keep level lower than actual conditions, then relative permeability RP curve and Kx, Ky distribute larger weight, if oil productive capacity is lower than actual conditions, relative permeability RP curve distributes larger weight, repetition second step calculates, if three indexs of prediction conform to substantially with actual conditions, then terminate weight allocation, stop calculating.
3. method as claimed in claim 2, is characterized in that, during adjustment relative permeability RP curve, if oil productive capacity is lower than actual conditions, then the adjustment of oil relative permeability forward increases, and aqueous phase relative permeability oppositely adjusts reduction; If oil productive capacity is higher than actual conditions, then adjustment direction is just the opposite;
The adjusting range of permeability is limited, after permeability adjusting range is more than 50%, auto-steering is adjusted conductivity;
Weight is arranged cannot revise permeability saturation curve form, needs finger prosthesis to determine permeability saturation curve form.
4. the method for claim 1, it is characterized in that, the method imports water saturation and the remaining oil reserves parameter field of any time by residual oil quantitative analytical system, is grade classification index, calculates the remaining oil reserves in different demarcation interval with water saturation.
5. method as claimed in claim 4, it is characterized in that, described residual oil quantitative analytical system comprises further:
Data acquisition module, for receiving input data volume file, adopting sequential organization algorithm for design to realize, judging wrong position and reason in file input process, the relevant entry in manual modification input data volume file, and again reads;
Residual oil quantitative analysis module, is connected with described data acquisition module, for receiving the data of described data collecting module collected, calculating the remaining oil reserves in each interval respectively, and form standard output form with given criteria for classifying interval;
Visualized data output module, is connected with described residual oil quantitative analysis module, exports with form and EXECL standard pie chart form for the standard output form remaining oil reserves data volume that will generate.
6. method as claimed in claim 5, is characterized in that, described residual oil quantitative analysis module adopts sequential organization algorithm for design to realize.
7. the method for claim 1, is characterized in that, microscopic remaining oil occurrence status adopts the experiment of sealed coring well oil-containing sheet laser confocal scanning to obtain, and detailed process is:
First, adopt freezing tabletting technology, carry out cutting sample at low temperature environment, when guaranteeing abrasive disc, original form of hole inner fluid is not destroyed;
Secondly, body of casting automatic image-analysis microscope, petrographic microscope and fluorescence microscope are combined, body of casting automatic image-analysis microscope is utilized to extract pore throat characteristic parameter, rock particles feature, the character of petrographic microscope identification mineral, fluorescence microscopy chooses Ultraluminescence filter, distinguish oils water boundary, the extraction of remaining oil occurrence status information is completed by residual oil quantitative analytical system, by observing fluoroscopic image, the remaining oil distribution state in hole being differentiated, characterizing the feature of remaining oil micro Distribution with this result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310085454.9A CN103206208B (en) | 2013-03-08 | 2013-03-08 | A kind of method of microcosmic different existence state remaining oil macro-level quantitative |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310085454.9A CN103206208B (en) | 2013-03-08 | 2013-03-08 | A kind of method of microcosmic different existence state remaining oil macro-level quantitative |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103206208A CN103206208A (en) | 2013-07-17 |
| CN103206208B true CN103206208B (en) | 2016-01-20 |
Family
ID=48753594
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310085454.9A Active CN103206208B (en) | 2013-03-08 | 2013-03-08 | A kind of method of microcosmic different existence state remaining oil macro-level quantitative |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103206208B (en) |
Families Citing this family (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103439302B (en) * | 2013-09-17 | 2015-09-16 | 黑龙江八一农垦大学 | The method of organic composite distribution and content in the burnt three-dimensional values mud stone of a kind of two-photon laser copolymerization and tight sand |
| CN104100263B (en) * | 2013-10-29 | 2017-01-25 | 中国石油化工股份有限公司 | Method for determining residual oil saturation of non-uniform water flooded layer |
| CN104747154B (en) * | 2015-03-05 | 2017-01-25 | 中国石油天然气股份有限公司 | Method for using oil displacement efficiency ratios for improving steam drive remaining oil research accuracy |
| CN106150454B (en) * | 2015-03-27 | 2018-11-30 | 中国石油化工股份有限公司 | Offshore oilfield vector quantization production and injection proration method |
| CN104833789B (en) * | 2015-04-09 | 2016-08-24 | 中国石油大学(北京) | Utilize the method that tight sand micropore structure determines water-gas relation |
| CN106337680A (en) * | 2015-07-13 | 2017-01-18 | 中国石油化工股份有限公司 | Method of determining recoverable reserves in oil reservoir production wells |
| CN106770377B (en) * | 2016-11-23 | 2019-05-07 | 东北石油大学 | The device and method of remaining oil distribution are monitored during carbon dioxide displacement core |
| CN108931497B (en) * | 2018-04-20 | 2021-06-01 | 中国石油天然气股份有限公司 | Method and device for evaluating hydrocarbon source rock through laser confocal three-dimensional reconstruction |
| CN109025924B (en) * | 2018-05-08 | 2024-04-05 | 中国海洋石油集团有限公司 | Oil saturation dynamic monitoring platform based on microcosmic rock slice |
| CN108468536B (en) * | 2018-05-08 | 2023-04-25 | 中国海洋石油集团有限公司 | Oil saturation monitoring method based on ultraviolet spectrum microcosmic oil displacement simulation experiment |
| CN109209312B (en) * | 2018-09-07 | 2020-11-13 | 中国石油化工股份有限公司 | Resource development method suitable for polymer flooding |
| CN111075439B (en) * | 2018-10-19 | 2023-02-10 | 中国石油天然气股份有限公司 | Method, device and storage medium for determining effective thickness of conglomerate reservoir |
| CN111175207B (en) * | 2018-11-12 | 2023-04-18 | 中石化石油工程技术服务有限公司 | Method for detecting effective porosity of rock |
| CN110439518B (en) * | 2019-08-14 | 2020-05-12 | 中国石油大学(华东) | Quantitative evaluation method for ternary composite displacement effect |
| CN110644979B (en) * | 2019-09-03 | 2021-08-27 | 中国石油大学(北京) | Method and device for acquiring initial occurrence state of pore fluid |
| CN110984971B (en) * | 2019-10-14 | 2023-06-02 | 中国石油天然气股份有限公司 | Microscopic residual oil classification method based on formation mechanism analysis |
| CN111236934B (en) * | 2020-02-25 | 2021-10-08 | 中国石油大学(北京) | Method and device for determining flooding level |
| CN112150424B (en) * | 2020-09-16 | 2023-02-24 | 中国石油大学(华东) | Microcosmic residual oil quantitative analysis method based on fluorescent thin sheet |
| CN112282748B (en) * | 2020-11-12 | 2022-04-15 | 中国石油大学(北京) | Residual oil occurrence form discrimination method based on micro-fluidic and machine learning |
| CN112576238B (en) * | 2020-12-02 | 2022-10-28 | 中国石油大学(华东) | System, method and application for determining position and content of residual oil in low-permeability reservoir |
| CN113959998B (en) * | 2021-09-30 | 2022-06-28 | 岩基电(苏州)精密仪器有限公司 | Method for determining crude oil components and occurrence states thereof based on fluorescence lifetime |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101074600A (en) * | 2007-06-18 | 2007-11-21 | 大庆石油学院 | Core-dissecting micro-analysis |
| CN201269098Y (en) * | 2008-10-23 | 2009-07-08 | 大庆油田有限责任公司 | Chart board used for seeking flooding degree |
| CN101942995A (en) * | 2010-09-16 | 2011-01-12 | 中国石油天然气股份有限公司 | Method for evaluating watered-out degree of ultralow-permeability water-drive reservoir sandstone oil layer |
| CN102926723A (en) * | 2012-11-20 | 2013-02-13 | 中国科学院地质与地球物理研究所 | Remaining oil cave extracting method of exhausted oil field |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2455016B (en) * | 2006-09-08 | 2011-10-19 | Bp Exploration Operating | Hydrocarbon recovery process |
-
2013
- 2013-03-08 CN CN201310085454.9A patent/CN103206208B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101074600A (en) * | 2007-06-18 | 2007-11-21 | 大庆石油学院 | Core-dissecting micro-analysis |
| CN201269098Y (en) * | 2008-10-23 | 2009-07-08 | 大庆油田有限责任公司 | Chart board used for seeking flooding degree |
| CN101942995A (en) * | 2010-09-16 | 2011-01-12 | 中国石油天然气股份有限公司 | Method for evaluating watered-out degree of ultralow-permeability water-drive reservoir sandstone oil layer |
| CN102926723A (en) * | 2012-11-20 | 2013-02-13 | 中国科学院地质与地球物理研究所 | Remaining oil cave extracting method of exhausted oil field |
Non-Patent Citations (3)
| Title |
|---|
| 微CT扫描重建低渗气藏微观孔隙结构_以新场气田上沙溪庙组储层为例;苏娜等;《石油与天然气地质》;20111031;第32卷(第54期);792-795 * |
| 微观水驱油实验及剩余油形成机理研究;苏娜等;《断块油气田》;20071130;第14卷(第6期);50-51 * |
| 高含水期油田剩余油分布研究;陆建林等;《石油学报》;20010930;第22卷(第5期);48-52 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103206208A (en) | 2013-07-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103206208B (en) | A kind of method of microcosmic different existence state remaining oil macro-level quantitative | |
| Ewaid et al. | Groundwater hydrochemistry assessment of North Dhi-Qar Province, South of Iraq using multivariate statistical techniques | |
| CN108416686A (en) | A kind of Eco-Geo-Environment Type division method based on Coal Resource Development | |
| CN106650812B (en) | A kind of urban water-body extracting method of satellite remote-sensing image | |
| CN112308290B (en) | Method for investigating and distributing points in polluted site | |
| Trevisani et al. | Surface texture analysis of a high-resolution DTM: Interpreting an alpine basin | |
| Torkashvand et al. | New hybrid evolutionary algorithm for optimizing index-based groundwater vulnerability assessment method | |
| CN103345566A (en) | Geochemical exploration anomaly recognizing and evaluating method based on contents contained in geology | |
| Gárfias‐Soliz et al. | Time series and stochastic analyses to study the hydrodynamic characteristics of karstic aquifers | |
| CN112459776A (en) | Quantitative evaluation method and device for deep and ultra-deep carbonate rock ancient oil reservoirs | |
| CN113720639A (en) | Method for investigation, supplement and sampling of polluted site | |
| Haider et al. | Identification of peneplains by multi‐parameter assessment of digital elevation models | |
| Jamshidi Gohari et al. | Improving permeability estimation of carbonate rocks using extracted pore network parameters: a gas field case study | |
| CA2921107C (en) | Method of quantitative calculation and assessment of bituman content in ancient carbonate reservoirs | |
| Dawson et al. | A method for the automated extraction of environmental variables to help the classification of rivers in Britain | |
| Dupas et al. | Landscape spatial configuration influences phosphorus but not nitrate concentrations in agricultural headwater catchments | |
| Wang et al. | Data-driven analysis of soil consolidation with prefabricated vertical drains considering stratigraphic variation | |
| CN116070147A (en) | Method for analyzing shale lithofacies and structure based on gray phase | |
| CN116109813A (en) | Anchor hole drilling identification method, system, electronic equipment and medium | |
| CN111691878B (en) | Lithofacies quality grading evaluation method and device for shale | |
| Gu et al. | Prediction of shell content from thin sections using hybrid image process techniques | |
| Bi et al. | Evaluating natural background levels of heavy metals in shallow groundwater of the Pearl River Delta via removal of contaminated groundwaters: Comparison of three preselection related methods | |
| CN114463621A (en) | Remote sensing quantitative analysis method and device for modern sediment source-sink system | |
| CN106468653B (en) | Method for judging clastic rock deposition microphase type | |
| Tanh | Quantifying and modeling the surficial karst of the Nullarbor Plain (southern Australia) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| ASS | Succession or assignment of patent right |
Owner name: INSTITUTE OF PETROLEUM EXPLORATION DEVELOPMENT OF Free format text: FORMER OWNER: PETROLEUM EXPLORATION DEVELOPMENT RESEARCH INSTITUTE, SINOPEC HENAN OILFIELD BRANCH Effective date: 20150526 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20150526 Address after: Henan City, Henan province 473132 Wancheng District of Nanyang oilfield Institute Applicant after: Institute of Petroleum Exploration Development of Henan Oilfield Branch, Sinopec Address before: Henan City, Henan province 473132 Wancheng District of Nanyang oilfield Institute Applicant before: Petroleum Exploration Development Research Institute, Sinopec Henan Oilfield Branch |
|
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |