CN106383078B - Method and device for determining water drive efficiency of rock - Google Patents
Method and device for determining water drive efficiency of rock Download PDFInfo
- Publication number
- CN106383078B CN106383078B CN201610835624.4A CN201610835624A CN106383078B CN 106383078 B CN106383078 B CN 106383078B CN 201610835624 A CN201610835624 A CN 201610835624A CN 106383078 B CN106383078 B CN 106383078B
- Authority
- CN
- China
- Prior art keywords
- value
- detected
- core sample
- porosity
- pixel
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 58
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 26
- 239000011148 porous material Substances 0.000 claims abstract description 66
- 238000006073 displacement reaction Methods 0.000 claims abstract description 12
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 18
- 229910052753 mercury Inorganic materials 0.000 claims description 18
- 229920006395 saturated elastomer Polymers 0.000 claims description 10
- 238000002347 injection Methods 0.000 claims description 9
- 239000007924 injection Substances 0.000 claims description 9
- 239000004744 fabric Substances 0.000 claims description 4
- 210000000867 larynx Anatomy 0.000 claims description 3
- 238000002591 computed tomography Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000035699 permeability Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000001035 drying Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The embodiment of the invention provides a method and a device for determining the water flooding efficiency of a rock, wherein the method comprises the following steps: calculating the porosity value of each pixel point based on the scanned image of the core sample to be detected, and calculating the distribution frequency of each porosity value; calculating the average major flow throat radius value of the core sample to be detected, wherein the major flow throat is a throat of which the size and connectivity in the core sample to be detected meet the oil displacement requirement; determining the pore distribution frequency of the core sample to be detected according to the average main flow throat radius value of the core sample to be detected, wherein the pore distribution frequency is the distribution frequency of all pores connected with the main flow throat with the radius larger than the average main flow throat radius value in the core sample to be detected; determining a porosity value with the distribution frequency matched with the pore distribution frequency as a porosity threshold value; and calculating the water flooding efficiency of the core sample to be detected according to the porosity value larger than the porosity threshold value. According to the scheme, the water flooding efficiency of the core sample to be detected can be rapidly and accurately calculated.
Description
Technical field
The present invention relates to oil-gas field development technical field, in particular to the determination method and dress of a kind of rock waterflood efficiency
It sets.
Background technique
Accurately and reliably core experiment result can provide number for establishment, the adjustment of Exploitation policy etc. of Oil and gas field development scheme
According to basis.Conventional core oil-displacement test is by experiment condition (for example, temperature, pressure), rock core property (for example, hypotonic, non-equal
Matter), human factor (for example, metering, calculate) etc. be affected, the time for completing experiment is longer, cost is also higher, it may be considered that
Other substitute technologies or method.
It is physico to have become the another progress after physical analogy for the computer modeling technique of flow event in pore media
Learn the important means of seepage research.Pore network model is on the basis of sufficiently obtaining rock core pore structure characteristic, with quantitatively pre-
Survey reservoir rock seepage characteristic be starting point, many aspects of physical chemistry seepage research achieve it is important into
Exhibition, and start to show its unique charm in many conventional methods almost helpless problem research.Now
Accurate porosity and permeability successfully can be calculated by the result of graphical analysis still, to imitate for water drive
The prediction of the dynamics displacement process such as rate is a difficult problem always.
Summary of the invention
The embodiment of the invention provides a kind of determination methods of rock waterflood efficiency, unpredictable in the prior art to solve
The technical issues of waterflood efficiency.This method comprises: calculating each pixel based on the scan image to core sample to be detected
Porosity value calculates the distribution frequency of each porosity value;Calculate the average mainstream throat radius of the core sample to be detected
Value, wherein mainstream venturi is the venturi that size and connectivity meet displacement of reservoir oil requirement in the core sample to be detected;According to
The average mainstream throat radius value of the core sample to be detected, determines the distribution of pores frequency of the core sample to be detected,
The distribution of pores frequency is the hole for the mainstream venturi connection that all radiuses are greater than the average mainstream throat radius value in institute
State the distribution frequency in core sample to be detected;Distribution frequency and the matched porosity value of distribution of pores frequency are determined as
Porosity threshold value;The waterflood efficiency of the core sample to be detected is calculated according to the porosity value for being greater than the porosity threshold value.
In one embodiment, based on the scan image to core sample to be detected, the porosity of each pixel is calculated
Value, comprising: the core sample to be detected under oven-dried condition is scanned using medical treatment CT, obtains the first CT of each pixel
Value;The core sample to be detected being completely soaked under saturated-water phase is scanned using medical CT, obtains each pixel
2nd CT value;According to the first CT value of each pixel and the 2nd CT value, the porosity value of each pixel is calculated.
In one embodiment, it is calculated every by following formula according to the first CT value of each pixel and the 2nd CT value
The porosity value of a pixel:Wherein, Φ is the porosity value of a pixel;CTdry
It is the first CT value of the next pixel of oven-dried condition;CTwetIt is the 2nd CT value for impregnating the next pixel of saturated-water phase;
CTwaterIt is the CT value of saturated water;CTairIt is the CT value of air.
In one embodiment, the average mainstream throat radius value of the core sample to be detected is calculated, comprising: to described
Core sample to be detected carries out constant-rate mercury injection, and the average mainstream larynx of the core sample to be detected is calculated according to intrusive mercury curve
Road radius value.
In one embodiment, according to the average mainstream throat radius value of the core sample to be detected, determine it is described to
Detect the distribution of pores frequency of core sample, comprising: be scanned, obtain to be detected to core sample to be detected using micron CT
The three-dimensional pore throat distributed image of core sample;In three-dimensional pore throat distributed image, according to being averaged for the core sample to be detected
Mainstream throat radius value determines the distribution of pores frequency of the core sample to be detected.
In one embodiment, by following formula according to the porosity value for being greater than the porosity threshold value calculate it is described to
Detect the waterflood efficiency of core sample, comprising:Wherein, EwtIt is waterflood efficiency;VvIt is one
The volume of a pixel;ΦmaxIt is maximum porosity value in all porosity values;ΦtIt is porosity threshold value, DΦBe it is described to
Detect the distribution of pores frequency of core sample;Φ is the porosity value of a pixel.
The embodiment of the invention also provides a kind of determining devices of rock waterflood efficiency, in the prior art can not be pre- with solution
The technical issues of surveying waterflood efficiency.The device includes: porosity value computing module, for based on sweeping to core sample to be detected
Trace designs picture, calculates the porosity value of each pixel, calculates the distribution frequency of each porosity value;Radius value computing module is used
In the average mainstream throat radius value for calculating the core sample to be detected, wherein mainstream venturi is the core sample to be detected
Size and connectivity meet the venturi of displacement of reservoir oil requirement in product;Distribution of pores frequency determining module, for according to described to be checked
The average mainstream throat radius value for surveying core sample, determines the distribution of pores frequency of the core sample to be detected, the hole
Distribution frequency is the hole for the mainstream venturi connection that all radiuses are greater than the average mainstream throat radius value described to be detected
Distribution frequency in core sample;Porosity threshold determination module, for matching distribution frequency with the distribution of pores frequency
Porosity value be determined as porosity threshold value;Waterflood efficiency computing module, for according to the hole for being greater than the porosity threshold value
Angle value calculates the waterflood efficiency of the core sample to be detected.
In one embodiment, further includes: medical CT equipment, for being carried out to the core sample to be detected under oven-dried condition
Scanning, obtains the first CT value of each pixel;The core sample to be detected being completely soaked under saturated-water phase is scanned,
Obtain the 2nd CT value of each pixel;The porosity value computing module, specifically for the first CT according to each pixel
Value and the 2nd CT value, calculate the porosity value of each pixel.
In one embodiment, the porosity value computing module, by following formula according to the first of each pixel
CT value and the 2nd CT value, calculate the porosity value of each pixel:Wherein, Φ is a picture
The porosity value of vegetarian refreshments;CTdryIt is the first CT value of the next pixel of oven-dried condition;CTwetIt is to impregnate one under saturated-water phase
2nd CT value of a pixel;CTwaterIt is the CT value of saturated water;CTairIt is the CT value of air.
In one embodiment, the radius value computing module, be specifically used for according to the core sample to be detected into
The intrusive mercury curve of row constant-rate mercury injection calculates the average mainstream throat radius value of the core sample to be detected.
In one embodiment, further includes: micron CT equipment is scanned core sample to be detected, obtains to be detected
The three-dimensional pore throat distributed image of core sample;The distribution of pores frequency determining module is specifically used in three-dimensional pore throat distribution map
As in, according to the average mainstream throat radius value of the core sample to be detected, the hole of the core sample to be detected is determined
Distribution frequency.
In one embodiment, waterflood efficiency computing module is specifically used for through following formula according to greater than the hole
The porosity value of degree threshold value calculates the waterflood efficiency of the core sample to be detected, comprising:Wherein, EwtIt is waterflood efficiency;VvIt is the volume of a pixel;ΦmaxIt is the largest hole
Porosity value;ΦtIt is porosity threshold value, DΦIt is the distribution of pores frequency of the core sample to be detected;Φ is a pixel
Porosity value.
In embodiments of the present invention, it by being scanned to core sample to be detected, can be calculated often based on scan image
The porosity value of a pixel calculates the distribution frequency of each porosity value;In the average mainstream for calculating core sample to be detected
After throat radius value, the distribution of pores frequency of core sample to be detected can be determined, and then matching hole point can be passed through
The mode of the distribution frequency of cloth frequency and porosity value, determines porosity threshold value, finally, according to the hole for being greater than porosity threshold value
Porosity value can calculate the water drive effect of core sample to be detected.Since distribution of pores frequency is to be greater than average mainstream based on radius
What the hole of the mainstream venturi connection of throat radius value calculated, enable the porosity threshold value determined to guarantee porosity value
Hole greater than the threshold value is substantially to connect with mainstream venturi, based on the porosity threshold value can quickly and accurately calculate to
Detect the waterflood efficiency of core sample.
Detailed description of the invention
The drawings described herein are used to provide a further understanding of the present invention, constitutes part of this application, not
Constitute limitation of the invention.In the accompanying drawings:
Fig. 1 is a kind of flow chart of the determination method of rock waterflood efficiency provided in an embodiment of the present invention;
Fig. 2 is the porosity value distribution frequency schematic diagram that a kind of medical CT scan provided in an embodiment of the present invention obtains;
Fig. 3 is a kind of constant speed intrusive mercury curve schematic diagram provided in an embodiment of the present invention;
Fig. 4 is the pore radius distribution frequency schematic diagram that a kind of micron of CT scan provided in an embodiment of the present invention obtains;
Fig. 5 is a kind of structural block diagram of the determining device of rock waterflood efficiency provided in an embodiment of the present invention.
Specific embodiment
To make the objectives, technical solutions, and advantages of the present invention clearer, right below with reference to embodiment and attached drawing
The present invention is described in further details.Here, exemplary embodiment and its explanation of the invention is used to explain the present invention, but simultaneously
It is not as a limitation of the invention.
Inventor has found that the mechanism of rock water drive oil is to increase pore pressure, by intrapore oil by supplementing energy
It is driven in venturi and flows out.There are many venturi quantity in rock, and throat radius is not of uniform size, but mainly pass through radius size greatly, even
The good venturi displacement of reservoir oil of the general character, i.e., by the displacement of reservoir oil of mainstream venturi, which refers to that size and connectivity meet the displacement of reservoir oil and want
The venturi asked, in order to accurately predict waterflood efficiency, inventor proposes a kind of determination method of rock waterflood efficiency, should
The determination method of rock waterflood efficiency determines the distribution of pores frequency of core sample to be detected according to average mainstream throat radius value,
The distribution frequency of distribution of pores frequency combination porosity value determines the porosity threshold value for calculating waterflood efficiency, and then according to big
The waterflood efficiency of core sample to be detected is calculated in the porosity value of porosity threshold value.
In embodiments of the present invention, a kind of determination method of rock waterflood efficiency is provided, as shown in Figure 1, this method packet
It includes:
Step 101: based on the scan image to core sample to be detected, calculating the porosity value of each pixel, calculate
The distribution frequency of each porosity value;
Step 102: calculating the average mainstream throat radius value of the core sample to be detected, wherein mainstream venturi is institute
State the venturi that size and connectivity in core sample to be detected meet displacement of reservoir oil requirement;
Step 103: according to the average mainstream throat radius value of the core sample to be detected, determining the rock core to be detected
The distribution of pores frequency of sample, the distribution of pores frequency are the mainstreams that all radiuses are greater than the average mainstream throat radius value
Distribution frequency of the hole of venturi connection in the core sample to be detected;
Step 104: distribution frequency and the matched porosity value of distribution of pores frequency are determined as porosity threshold value;
Step 105: the water drive of the core sample to be detected is calculated according to the porosity value for being greater than the porosity threshold value
Efficiency.
Process as shown in Figure 1 is it is found that in embodiments of the present invention, by being scanned to core sample to be detected, base
The porosity value of each pixel can be calculated in scan image, calculates the distribution frequency of each porosity value;It is to be checked calculating
It surveys after the average mainstream throat radius value of core sample, can determine the distribution of pores frequency of core sample to be detected, into
And porosity threshold value can be determined, finally, root by way of the distribution frequency of matching hole distribution frequency and porosity value
The water drive effect of core sample to be detected can be calculated according to the porosity value for being greater than porosity threshold value.Since distribution of pores frequency is
It is calculated based on the hole that radius is greater than the mainstream venturi connection of average mainstream throat radius value, so that the porosity determined
The hole that threshold value can guarantee that porosity value is greater than the threshold value is substantially to connect with mainstream venturi, can based on the porosity threshold value
Quickly and accurately to calculate the waterflood efficiency of core sample to be detected.
When it is implemented, in order to realize the porosity value for calculating each pixel, in the present embodiment, based on to be detected
The scan image of core sample calculates the porosity value of each pixel, comprising: uses medical treatment CT (Computed
Tomography, CT scan) core sample to be detected under oven-dried condition is scanned, obtain each picture
First CT value of vegetarian refreshments;The core sample to be detected being completely soaked under saturated-water phase is scanned using medical CT, is obtained
2nd CT value of each pixel;According to the first CT value of each pixel and the 2nd CT value, the hole of each pixel is calculated
Angle value.
Specifically, being scanned after the drying of step 1, core sample to be detected with medical CT, writing scan voltage, scanning
Position and the condition of scanning obtain the first CT value CT of each pixel of core sample to be detecteddry;
Step 2 evacuates core sample to be detected, after being completely soaked 100% saturated water, in scanning identical with step 1
Under voltage, the condition of scanning and scan position, the core sample to be detected being completely soaked under saturated-water phase is carried out with medical treatment CT
Scanning, obtains the 2nd CT value CT of each pixel of rock core of fully saturated waterwet;
Step 3, scanning obtain air and the saturation CT value CT of waterair、CTwater;
Step 4, by following formula according to the first CT value of each pixel and the 2nd CT value, calculate each pixel
Porosity value:
Wherein, Φ is the porosity value of a pixel;CTdryIt is the first CT value of the next pixel of oven-dried condition;
CTwetIt is the 2nd CT value for impregnating the next pixel of saturated-water phase;CTwaterIt is the CT value of saturated water;CTairIt is the CT of air
Value.
When it is implemented, after calculating the porosity value of each pixel, so that it may count point of each porosity value
Cloth frequency, i.e., each porosity value percentage shared in all porosity values are specific as shown in Fig. 2, abscissa is hole
Angle value, ordinate are frequency values.
When it is implemented, in order to which the average mainstream throat radius value of core sample to be detected can be calculated, in the present embodiment
In, calculate the average mainstream throat radius value of the core sample to be detected, comprising: carry out to the core sample to be detected permanent
Fast Mercury injection calculates the average mainstream throat radius value of the core sample to be detected according to intrusive mercury curve.Specifically, treating
It detects core sample and carries out constant-rate mercury injection, obtain intrusive mercury curve (as shown in figure 3, dashed curve indicates hole, line in Fig. 3
Curve indicates that venturi, the curve of lines combination dot indicate overall), according to the intrusive mercury curve calculating core sample to be detected
Average mainstream throat radius value, which can be permeability accumulation contribution up to middle and high infiltration 95%, low
Seep the weighted average of 80% pervious all throat radius values.
When it is implemented, being swept after obtaining the average mainstream throat radius value of core sample to be detected using micron CT equipment
Core sample to be detected is retouched, the three-dimensional pore throat distributed image of core sample to be detected, the root in three-dimensional pore throat distributed image are obtained
According to the average mainstream throat radius value of core sample to be detected, the distribution of pores frequency of core sample to be detected, the hole are determined
Distribution frequency refers to that radius value is greater than the radius of the hole of the mainstream venturi connection of average mainstream throat radius value in rock to be detected
Percentage in all holes of heart sample, it is specific as shown in Figure 4.
When it is implemented, after determining porosity threshold value, by following formula according to greater than the porosity threshold value
Porosity value calculates the waterflood efficiency of the core sample to be detected, comprising:
Wherein, EwtIt is waterflood efficiency;VvIt is the volume of a pixel;ΦmaxIt is maximum hole in all porosity values
Angle value;ΦtIt is porosity threshold value, DΦIt is the distribution of pores frequency of the core sample to be detected;Φ is the hole of a pixel
Porosity value.
The determination method of above-mentioned rock waterflood efficiency is described in detail below in conjunction with specific example, this method includes following step
It is rapid:
(1) core sample to be detected: two pieces of reservoir rocks, the porosity of rock core 1 are 27.2%, the permeability of air
1123mD;The porosity of rock core 2 is 11.4%, the permeability 1.54mD of air.
(2) instrument: GE LightSpeed8 medical CT scanners, MicroXCT-200 microns of CT scanners of Xradia,
ASPE-730 constant speed mercury injection apparatus.
(3) experimentation:
Two pieces of core samples to be detected are handled according to the following steps:
1. be scanned under 100kv scanning voltage after two pieces of core samples to be detected drying, writing scan position and
The condition of scanning obtains the first CT value CT of two pieces of dry each pixels of rock coredry;
2. two pieces of core samples to be detected are evacuated, after being completely soaked 100% saturated water, in scanning identical with step 1
Under voltage, the condition of scanning and scan position, two pieces of core samples to be detected are scanned, obtain two blocks of rocks of fully saturated water
2nd CT value CT of each pixel of the heartwet;
3. scanning obtains air and the saturation CT value CT of waterair、CTwater, calculate the hole of each pixel of every block of rock core
Porosity value, and count the distribution frequency of porosity value;
4. two blocks of rock cores (i.e. above-mentioned core sample to be detected) are dried again, then swept using micron CT scanner
It retouches, scanning voltage 40kV, 1.02 microns of resolution ratio;
5. pair two blocks of rock cores carry out constant-rate mercury injection, average mainstream throat radius value is calculated according to intrusive mercury curve, respectively
It is 1,13.5 microns of rock core, 2,3.1 microns of rock core;
6. image carries out three-dimensional reconstruction after two pieces of rock core micron CT scan, three-dimensional pore throat distributed image is obtained, according to average
Mainstream throat radius Data-Statistics radius value is greater than the distribution of the pore radius of the mainstream venturi connection of average mainstream throat radius value
Frequency, respectively rock core 1,15.7%, rock core 2,11.3%.By the progress of the distribution frequency of distribution of pores frequency and porosity value
Match, determines that the respective porosity threshold value of two blocks of rock cores is respectively 23.5%, 10.4%.
7. calculating theoretical waterflood efficiency using the porosity value for being greater than porosity threshold value, two blocks of rock cores are respectively rock core 1,
61.9%, rock core 2,53.5%.
8. carrying out rock core flood pot test simultaneously, the practical waterflood efficiency of rock core 1 is 60.3%, the practical water drive of rock core 2
Efficiency is 51.5%, and the error of prediction is respectively 2.7% and 3.9%.
Based on the same inventive concept, a kind of determining device of rock waterflood efficiency is additionally provided in the embodiment of the present invention, such as
Described in the following examples.The determination of the principle and rock waterflood efficiency that are solved the problems, such as due to the determining device of rock waterflood efficiency
Method is similar, therefore the implementation of the determining device of rock waterflood efficiency may refer to the reality of the determination method of rock waterflood efficiency
It applies, overlaps will not be repeated.Used below, the software of predetermined function may be implemented in term " unit " or " module "
And/or the combination of hardware.Although device described in following embodiment is preferably realized with software, hardware or soft
The realization of the combination of part and hardware is also that may and be contemplated.
Fig. 5 is a kind of structural block diagram of the determining device of the rock waterflood efficiency of the embodiment of the present invention, as shown in figure 5, packet
It includes:
Porosity value computing module 501, for calculating each pixel based on the scan image to core sample to be detected
Porosity value, calculate the distribution frequency of each porosity value;
Radius value computing module 502, for calculating the average mainstream throat radius value of the core sample to be detected,
In, mainstream venturi is the venturi that size and connectivity meet displacement of reservoir oil requirement in the core sample to be detected;
Distribution of pores frequency determining module 503, for the average mainstream throat radius according to the core sample to be detected
Value determines that the distribution of pores frequency of the core sample to be detected, the distribution of pores frequency are that all radiuses are greater than described put down
Distribution frequency of the hole of the mainstream venturi connection of equal mainstream throat radius value in the core sample to be detected;
Porosity threshold determination module 504 is used for distribution frequency and the matched porosity value of distribution of pores frequency
It is determined as porosity threshold value;
Waterflood efficiency computing module 505, it is described to be checked for being calculated according to the porosity value for being greater than the porosity threshold value
Survey the waterflood efficiency of core sample.
In one embodiment, further includes: medical CT equipment, for being carried out to the core sample to be detected under oven-dried condition
Scanning, obtains the first CT value of each pixel;The core sample to be detected being completely soaked under saturated-water phase is scanned,
Obtain the 2nd CT value of each pixel;The porosity value computing module, specifically for the first CT according to each pixel
Value and the 2nd CT value, calculate the porosity value of each pixel.
In one embodiment, the porosity value computing module, by following formula according to the first of each pixel
CT value and the 2nd CT value, calculate the porosity value of each pixel:Wherein, Φ is a picture
The porosity value of vegetarian refreshments;CTdryIt is the first CT value of the next pixel of oven-dried condition;CTwetIt is to impregnate one under saturated-water phase
2nd CT value of a pixel;CTwaterIt is the CT value of saturated water;CTairIt is the CT value of air.
In one embodiment, the radius value computing module, be specifically used for according to the core sample to be detected into
The intrusive mercury curve of row constant-rate mercury injection calculates the average mainstream throat radius value of the core sample to be detected.
In one embodiment, further includes: micron CT equipment is scanned core sample to be detected, obtains to be detected
The three-dimensional pore throat distributed image of core sample;The distribution of pores frequency determining module is specifically used in three-dimensional pore throat distribution map
As in, according to the average mainstream throat radius value of the core sample to be detected, the hole of the core sample to be detected is determined
Distribution frequency.
In one embodiment, waterflood efficiency computing module is specifically used for through following formula according to greater than the hole
The porosity value of degree threshold value calculates the waterflood efficiency of the core sample to be detected, comprising:Wherein, EwtIt is waterflood efficiency;VvIt is the volume of a pixel;ΦmaxIt is the largest hole
Porosity value;ΦtIt is porosity threshold value, DΦIt is the distribution of pores frequency of the core sample to be detected;Φ is a pixel
Porosity value.
In another embodiment, a kind of software is additionally provided, the software is for executing above-described embodiment and preferred reality
Apply technical solution described in mode.
In another embodiment, a kind of storage medium is additionally provided, above-mentioned software is stored in the storage medium, it should
Storage medium includes but is not limited to: CD, floppy disk, hard disk, scratch pad memory etc..
In embodiments of the present invention, it by being scanned to core sample to be detected, can be calculated often based on scan image
The porosity value of a pixel calculates the distribution frequency of each porosity value;In the average mainstream for calculating core sample to be detected
After throat radius value, the distribution of pores frequency of core sample to be detected can be determined, and then matching hole point can be passed through
The mode of the distribution frequency of cloth frequency and porosity value, determines porosity threshold value, finally, according to the hole for being greater than porosity threshold value
Porosity value can calculate the water drive effect of core sample to be detected.Since distribution of pores frequency is to be greater than average mainstream based on radius
What the hole of the mainstream venturi connection of throat radius value calculated, enable the porosity threshold value determined to guarantee porosity value
Hole greater than the threshold value is substantially to connect with mainstream venturi, based on the porosity threshold value can quickly and accurately calculate to
Detect the waterflood efficiency of core sample.
Obviously, those skilled in the art should be understood that each module of the above-mentioned embodiment of the present invention or each step can be with
It is realized with general computing device, they can be concentrated on a single computing device, or be distributed in multiple computing devices
On composed network, optionally, they can be realized with the program code that computing device can perform, it is thus possible to by it
Store and be performed by computing device in the storage device, and in some cases, can be held with the sequence for being different from herein
The shown or described step of row, perhaps they are fabricated to each integrated circuit modules or will be multiple in them
Module or step are fabricated to single integrated circuit module to realize.In this way, the embodiment of the present invention be not limited to it is any specific hard
Part and software combine.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field
For art personnel, the embodiment of the present invention can have various modifications and variations.All within the spirits and principles of the present invention, made
Any modification, equivalent substitution, improvement and etc. should all be included in the protection scope of the present invention.
Claims (10)
1. a kind of determination method of rock waterflood efficiency characterized by comprising
Based on the scan image to core sample to be detected, the porosity value of each pixel is calculated, calculates each porosity value
Distribution frequency;
Calculate the average mainstream throat radius value of the core sample to be detected, wherein mainstream venturi is the rock core to be detected
Size and connectivity meet the venturi of displacement of reservoir oil requirement in sample;
According to the average mainstream throat radius value of the core sample to be detected, the hole point of the core sample to be detected is determined
Cloth frequency, the distribution of pores frequency are the holes for the mainstream venturi connection that all radiuses are greater than the average mainstream throat radius value
Distribution frequency of the gap in the core sample to be detected;
Distribution frequency and the matched porosity value of distribution of pores frequency are determined as porosity threshold value;
The waterflood efficiency of the core sample to be detected is calculated according to the porosity value for being greater than the porosity threshold value;
The water drive of the core sample to be detected is calculated according to the porosity value for being greater than the porosity threshold value by following formula
Efficiency, comprising:
Wherein, EwtIt is waterflood efficiency;VvIt is the volume of a pixel;ΦmaxIt is maximum porosity in all porosity values
Value;ΦtIt is porosity threshold value, DΦIt is the distribution of pores frequency of the core sample to be detected;Φ is the hole of a pixel
Angle value.
2. the determination method of rock waterflood efficiency as described in claim 1, which is characterized in that based on to core sample to be detected
Scan image, calculate the porosity value of each pixel, comprising:
The core sample to be detected under oven-dried condition is scanned using medical CT, obtains the first CT value of each pixel;
The core sample to be detected being completely soaked under saturated-water phase is scanned using medical CT, obtains each pixel
2nd CT value;
According to the first CT value of each pixel and the 2nd CT value, the porosity value of each pixel is calculated.
3. the determination method of rock waterflood efficiency as claimed in claim 2, which is characterized in that by following formula according to each
First CT value of pixel and the 2nd CT value calculate the porosity value of each pixel:
Wherein, Φ is the porosity value of a pixel;CTdryIt is the first CT value of the next pixel of oven-dried condition;CTwetIt is
Impregnate the 2nd CT value of the next pixel of saturated-water phase;CTwaterIt is the CT value of saturated water;CTairIt is the CT value of air.
4. the determination method of rock waterflood efficiency as claimed any one in claims 1 to 3, which is characterized in that described in calculating
The average mainstream throat radius value of core sample to be detected, comprising:
Constant-rate mercury injection is carried out to the core sample to be detected, the core sample to be detected is calculated according to intrusive mercury curve
Average mainstream throat radius value.
5. the determination method of rock waterflood efficiency as claimed any one in claims 1 to 3, which is characterized in that according to described
The average mainstream throat radius value of core sample to be detected, determines the distribution of pores frequency of the core sample to be detected, comprising:
Core sample to be detected is scanned using micron CT, obtains the three-dimensional pore throat distributed image of core sample to be detected;
In three-dimensional pore throat distributed image, according to the average mainstream throat radius value of the core sample to be detected, determine described in
The distribution of pores frequency of core sample to be detected.
6. a kind of determining device of rock waterflood efficiency characterized by comprising
Porosity value computing module, for calculating the hole of each pixel based on the scan image to core sample to be detected
Angle value calculates the distribution frequency of each porosity value;
Radius value computing module, for calculating the average mainstream throat radius value of the core sample to be detected, wherein mainstream larynx
Road is the venturi that size and connectivity meet displacement of reservoir oil requirement in the core sample to be detected;
Distribution of pores frequency determining module is determined for the average mainstream throat radius value according to the core sample to be detected
The distribution of pores frequency of the core sample to be detected, the distribution of pores frequency are that all radiuses are greater than the average mainstream larynx
Distribution frequency of the hole of the mainstream venturi connection of road radius value in the core sample to be detected;
Porosity threshold determination module, for distribution frequency and the matched porosity value of distribution of pores frequency to be determined as hole
Porosity threshold value;
Waterflood efficiency computing module, for calculating the core sample to be detected according to the porosity value for being greater than the porosity threshold value
The waterflood efficiency of product;
Waterflood efficiency computing module, specifically for being calculated by following formula according to the porosity value for being greater than the porosity threshold value
The waterflood efficiency of the core sample to be detected, comprising:
Wherein, EwtIt is waterflood efficiency;VvIt is the volume of a pixel;ΦmaxIt is the largest porosity value;ΦtIt is porosity threshold
Value, DΦIt is the distribution of pores frequency of the core sample to be detected;Φ is the porosity value of a pixel.
7. the determining device of rock waterflood efficiency as claimed in claim 6, which is characterized in that further include:
Medical CT equipment obtains the first of each pixel for being scanned to the core sample to be detected under oven-dried condition
CT value;The core sample to be detected being completely soaked under saturated-water phase is scanned, the 2nd CT value of each pixel is obtained;
The porosity value computing module, specifically for calculating each according to the first CT value of each pixel and the 2nd CT value
The porosity value of pixel.
8. the determining device of rock waterflood efficiency as claimed in claim 7, which is characterized in that the porosity value calculates mould
Block calculates the porosity value of each pixel by following formula according to the first CT value of each pixel and the 2nd CT value:
Wherein, Φ is the porosity value of a pixel;CTdryIt is the first CT value of the next pixel of oven-dried condition;CTwetIt is
Impregnate the 2nd CT value of the next pixel of saturated-water phase;CTwaterIt is the CT value of saturated water;CTairIt is the CT value of air.
9. the determining device of the rock waterflood efficiency as described in any one of claim 6 to 8, which is characterized in that the radius
It is worth computing module, specifically for calculating institute according to the intrusive mercury curve for carrying out constant-rate mercury injection to the core sample to be detected
State the average mainstream throat radius value of core sample to be detected.
10. the determining device of the rock waterflood efficiency as described in any one of claim 6 to 8, which is characterized in that further include:
Micron CT equipment, is scanned core sample to be detected, obtains the three-dimensional pore throat distributed image of core sample to be detected;
The distribution of pores frequency determining module is specifically used in three-dimensional pore throat distributed image, according to the rock core to be detected
The average mainstream throat radius value of sample, determines the distribution of pores frequency of the core sample to be detected.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610835624.4A CN106383078B (en) | 2016-09-20 | 2016-09-20 | Method and device for determining water drive efficiency of rock |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610835624.4A CN106383078B (en) | 2016-09-20 | 2016-09-20 | Method and device for determining water drive efficiency of rock |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106383078A CN106383078A (en) | 2017-02-08 |
CN106383078B true CN106383078B (en) | 2019-05-07 |
Family
ID=57935808
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610835624.4A Active CN106383078B (en) | 2016-09-20 | 2016-09-20 | Method and device for determining water drive efficiency of rock |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106383078B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107525753A (en) * | 2017-07-25 | 2017-12-29 | 中国石油大学(北京) | Determine the method and device of hole tortuosity |
CN108507920B (en) * | 2018-03-02 | 2020-12-01 | 中国石油天然气股份有限公司 | Method for simulating relative permeability of oil phase and water phase |
CN108492331B (en) * | 2018-03-20 | 2019-04-23 | 中国科学院地质与地球物理研究所 | The method and apparatus of the pore-size distribution of rock are obtained based on simulation mercury injection method |
CN113124808B (en) * | 2021-04-12 | 2022-08-12 | 陕西延长石油(集团)有限责任公司 | Method for evaluating matching of initial radius of microsphere and pore throat structure of reservoir |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5359194A (en) * | 1992-05-01 | 1994-10-25 | Texaco Inc. | X-ray CT measurement of secondary (vugular) porosity in reservoir core material |
CN101344001A (en) * | 2008-08-05 | 2009-01-14 | 中国石化集团华北石油局 | Analytical method of X-ray fluorescence terrigenous clastic rock porosity in petroleum well drilling |
CN102094642A (en) * | 2010-12-17 | 2011-06-15 | 中国石油天然气股份有限公司 | In-layer heterogeneous model oil-water displacement efficiency evaluation system |
CN102519853A (en) * | 2011-12-05 | 2012-06-27 | 中国石油大学(华东) | Calculating method of porosity lower limit of recovery industry for effective reservoir under pore structure restraint |
CN103198363A (en) * | 2013-03-19 | 2013-07-10 | 中国石油天然气股份有限公司 | Reservoir gas production amount prediction method and device based on CT pore analysis |
CN103306671A (en) * | 2013-05-17 | 2013-09-18 | 中国石油天然气股份有限公司 | Four-quadrant reservoir type identification method and system |
CN103674802A (en) * | 2013-11-05 | 2014-03-26 | 中国石油天然气股份有限公司 | Rock closed porosity determination method |
CN103698803A (en) * | 2012-09-27 | 2014-04-02 | 中国石油天然气股份有限公司 | Rock pore structure characterization method and device |
CN104458525A (en) * | 2013-09-13 | 2015-03-25 | 中国石油天然气股份有限公司 | Acquisition and processing method for characterizing characteristics of micro-pore throat by conventional mercury intrusion curve |
CN104535475A (en) * | 2015-01-08 | 2015-04-22 | 中国石油大学(北京) | Carbonate rock microstructure determination method and device |
CN104748906A (en) * | 2015-04-10 | 2015-07-01 | 中国石油大学(北京) | Device and method for determining capillary pressure of carbonatite |
CN104794709A (en) * | 2015-04-10 | 2015-07-22 | 四川大学 | Three-dimensional core image pore and throat segmentation method |
CN104897545A (en) * | 2015-06-09 | 2015-09-09 | 中国石油天然气股份有限公司 | Rock core pore structure change detection and analysis method |
CN105891089A (en) * | 2016-05-18 | 2016-08-24 | 中国石油大学(北京) | Method and device for determining permeability of reservoir |
CA2871781C (en) * | 2012-05-18 | 2017-02-07 | Ingrain, Inc. | Method and system for estimating rock properties from rock samples using digital rock physics imaging |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8170799B2 (en) * | 2008-11-24 | 2012-05-01 | Ingrain, Inc. | Method for determining in-situ relationships between physical properties of a porous medium from a sample thereof |
US8583410B2 (en) * | 2010-05-28 | 2013-11-12 | Ingrain, Inc. | Method for obtaining consistent and integrated physical properties of porous media |
-
2016
- 2016-09-20 CN CN201610835624.4A patent/CN106383078B/en active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5359194A (en) * | 1992-05-01 | 1994-10-25 | Texaco Inc. | X-ray CT measurement of secondary (vugular) porosity in reservoir core material |
CN101344001A (en) * | 2008-08-05 | 2009-01-14 | 中国石化集团华北石油局 | Analytical method of X-ray fluorescence terrigenous clastic rock porosity in petroleum well drilling |
CN102094642A (en) * | 2010-12-17 | 2011-06-15 | 中国石油天然气股份有限公司 | In-layer heterogeneous model oil-water displacement efficiency evaluation system |
CN102519853A (en) * | 2011-12-05 | 2012-06-27 | 中国石油大学(华东) | Calculating method of porosity lower limit of recovery industry for effective reservoir under pore structure restraint |
CA2871781C (en) * | 2012-05-18 | 2017-02-07 | Ingrain, Inc. | Method and system for estimating rock properties from rock samples using digital rock physics imaging |
CN103698803A (en) * | 2012-09-27 | 2014-04-02 | 中国石油天然气股份有限公司 | Rock pore structure characterization method and device |
CN103198363A (en) * | 2013-03-19 | 2013-07-10 | 中国石油天然气股份有限公司 | Reservoir gas production amount prediction method and device based on CT pore analysis |
CN103306671A (en) * | 2013-05-17 | 2013-09-18 | 中国石油天然气股份有限公司 | Four-quadrant reservoir type identification method and system |
CN104458525A (en) * | 2013-09-13 | 2015-03-25 | 中国石油天然气股份有限公司 | Acquisition and processing method for characterizing characteristics of micro-pore throat by conventional mercury intrusion curve |
CN103674802A (en) * | 2013-11-05 | 2014-03-26 | 中国石油天然气股份有限公司 | Rock closed porosity determination method |
CN104535475A (en) * | 2015-01-08 | 2015-04-22 | 中国石油大学(北京) | Carbonate rock microstructure determination method and device |
CN104748906A (en) * | 2015-04-10 | 2015-07-01 | 中国石油大学(北京) | Device and method for determining capillary pressure of carbonatite |
CN104794709A (en) * | 2015-04-10 | 2015-07-22 | 四川大学 | Three-dimensional core image pore and throat segmentation method |
CN104897545A (en) * | 2015-06-09 | 2015-09-09 | 中国石油天然气股份有限公司 | Rock core pore structure change detection and analysis method |
CN105891089A (en) * | 2016-05-18 | 2016-08-24 | 中国石油大学(北京) | Method and device for determining permeability of reservoir |
Also Published As
Publication number | Publication date |
---|---|
CN106383078A (en) | 2017-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106383078B (en) | Method and device for determining water drive efficiency of rock | |
CN106323836B (en) | Borehole wall permeability calculation method and device | |
CN104902254B (en) | For the method and device that the disparity map of stereo-picture is estimated | |
CN109544442B (en) | Image local style migration method of double-countermeasure-based generation type countermeasure network | |
CN106056562A (en) | Face image processing method and device and electronic device | |
CN104751473B (en) | Method and device for determining characteristics of multi-scale holes of carbonatite | |
CN108802073A (en) | A kind of rock electrical parameters acquisition methods and device based on digital cores | |
CN105678248A (en) | Face key point alignment algorithm based on deep learning | |
CN109632604B (en) | Method for coarsening relative permeability of polymer flooding from pore size to core size | |
CN105844625A (en) | Movable profile image segmentation method fusing edge and area | |
CN105631872B (en) | Remote sensing image registration method based on multi-characteristic points | |
CN107993228A (en) | A kind of vulnerable plaque automatic testing method and device based on cardiovascular OCT images | |
WO2011131029A1 (en) | Method for detecting similar units based on outline belt graph | |
CN108830856A (en) | A kind of GA automatic division method based on time series SD-OCT retinal images | |
CN109142875A (en) | A method of tight sand rock electrology characteristic is obtained using digital cores | |
CN109663357A (en) | Game special implementation method and device | |
CN113899768B (en) | Fluid injection-based method for analyzing pore connectivity of compact sandstone core | |
CN106708723A (en) | Testing and verifying methods and apparatuses | |
Fessel et al. | Structuring precedes extension in percolating Physarum polycephalum networks | |
CN109934906B (en) | Crude oil saturation acquisition method and device, computer equipment and storage medium | |
CN107610216A (en) | Video camera based on the three-dimensional point cloud generation method of particle group optimizing various visual angles and application | |
CN108734099A (en) | Auth method and device, electronic equipment, computer program and storage medium | |
CN109389651A (en) | Magnetic resonance chemical shift coded imaging method and apparatus | |
CN111141768A (en) | Method and device for acquiring rock wettability change in crude oil displacement | |
CN107403448A (en) | Cost function generation method and cost function generating means |
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 |