CN109444190A - The method for improving complex hydrocarbon layer water saturation evaluation precision is analyzed using digital cores - Google Patents
The method for improving complex hydrocarbon layer water saturation evaluation precision is analyzed using digital cores Download PDFInfo
- Publication number
- CN109444190A CN109444190A CN201811296212.3A CN201811296212A CN109444190A CN 109444190 A CN109444190 A CN 109444190A CN 201811296212 A CN201811296212 A CN 201811296212A CN 109444190 A CN109444190 A CN 109444190A
- Authority
- CN
- China
- Prior art keywords
- water saturation
- layer
- complex hydrocarbon
- hydrocarbon layer
- saturation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/20—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
- G01N23/203—Measuring back scattering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/22—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/05—Investigating materials by wave or particle radiation by diffraction, scatter or reflection
- G01N2223/053—Investigating materials by wave or particle radiation by diffraction, scatter or reflection back scatter
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/07—Investigating materials by wave or particle radiation secondary emission
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/616—Specific applications or type of materials earth materials
Abstract
The invention discloses a kind of methods analyzed using digital cores and improve complex hydrocarbon layer water saturation evaluation precision, belong to the technical field of petroleum exploration domain, comprising the following steps: (1) prepare test sample；(2) test sample carries out backscattered electron and secondary electron imaging；(3) secondary electron is imaged and carries out pore volume segmentation and statistics, to obtain pore morphology parameter；(4) automatic ensaying is carried out to Precise imaging, carries out automatic ensaying to acquire its characteristic atomic spectrogram to characteristic atomic spectrogram, to obtain additional conductive mineral content；(5) saturation litho-electric parameters determine, and establish the transformational relation of pore morphology parameter and m and the transformational relation of additional conductive mineral and n；(6) water saturation of complex hydrocarbon layer is calculated；To reach the distribution by extracting the smaller aperture for influencing saturation degree evaluation and the content of crucial additional electric conductivity mineral, and then obtain the purpose of accurate quantification water saturation numerical value.
Description
Technical field
The invention belongs to the technical fields of petroleum exploration domain, are mentioned in particular to a kind of using digital cores analysis
The method of high complex hydrocarbon layer water saturation evaluation precision.
Background technique
In oil and gas exploration and development process, water saturation seek be most critical oil and gas zone evaluation parameter,
Currently, industrial application is sought with Archie formula, Archie formula is specific as follows:
The formula shows stratum water saturation size and formation water resistivity Rw is almost directly proportional and formation resistivity
Rt is almost inversely proportional, meanwhile, by rock porosity sizePore structure factor m, saturation degree factor of n.
In archie formula, a regulation coefficient related with shale ingredient in rock, content and distribution form, generally 1
Left and right；B is as the additional saturation degree micro Distribution coefficient of uncertainty of wetability, generally 1 or so；M be pore structure because
Son is the measurement of hole tortuosity caused by skeleton and pore network, and classical archie sandstone m is near 2；N be saturation degree because
Son, classical archie sandstone n is near 2；
The porosity of rock can be obtained by neutron, density and acoustic logging；
Rw: the series of strata formation water resistivity measured value；
Rt: stratum actual resistivity measured value；
This method is coincide very well in the high middle and high infiltration reservoir in hole with real data in routine, is used widely, usually logical
Core experiment analysis is crossed, the porosity size of reservoir is obtainedThe litho-electric parameters such as a, b, m, n, since boring sample is few, rock core is real
It tests that analysis cost is high, the period is long, leads to not to each layer of progresss sample test, a usual area is in several pieces of samples of test
Afterwards, a regional experience value is counted, is used for all layers of reservoir evaluation later.A set of stratum of the parameters of formation water Rw in somewhere
In generally also think constant, the factor of influence saturation degree so just only leaves porosity φ and resistivity Rt, in many areas,
The porosity difference on certain set stratum is little, and such resistivity just parameter the most key at identification of hydrocarbon, resistivity is bigger, ground
The oiliness of layer is better.But due to each stratum sedimentary difference, each practical litho-electric parameters of layer differ greatly, and cause resistivity is high to have
When be water layer, and the sometimes oil reservoir that resistivity is low, and practical hydrocarbon characteristic is widely different, directly affects and utilizes Archie public
Formula carries out the precision that water saturation is sought.
To overcome drawbacks described above, people have carried out multiple improvement on the basis of Archie formula, as Indonesia's formula, Su Men are answered
Cured formula, dual water model etc., but effect be not it is ideal, can not objective quantitative ground-to-ground layer water saturation is asked
Take, especially in recent years, oil-gas exploration and development from conventional oil gas reservoir be transitioned into complicated reservoirs (low porosity and low permeability reservoir, compact reservoir,
Shale gas, low-resistivity reservoir etc.), quantitative saturation degree is sought becoming the bottleneck for restricting complicated reservoirs evaluation.
Oil field total rock heart coring is at high cost, and time-consuming, can not layer-by-layer coring, but with the raising of technology, currently, oily
Field scene can be carried out with relatively low cost, carry out sidewall coring in each emphasis layer, while there are also a large amount of for Oil Field
Landwaste data, these data are the important materials of later period oil reservoir saturation degree fine evaluation, but due to during usual coring by
Mud flush, mechanical impact etc. lead to not directly carry out seeking for original water saturation using side-wall coring etc..
Summary of the invention
In consideration of it, in order to solve the above problems existing in the present technology, utilizing digital cores point the present invention provides a kind of
The method of analysis raising complex hydrocarbon layer water saturation evaluation precision can utilize live side-wall coring or landwaste data to reach,
Stratum Analysis of micro pore structure and mineral composition analysis are carried out, point of the fine pore and venturi that influence saturation degree evaluation is extracted
The content of cloth and crucial additional electric conductivity mineral, and then obtain the purpose of accurate quantification water saturation numerical value.
The technical scheme adopted by the invention is as follows: it is a kind of to improve complex hydrocarbon layer water saturation using digital cores analysis
The method of evaluation precision, comprising the following steps:
(1) test sample is prepared；
(2) backscattered electron is carried out to the test sample of step (1) and secondary electron is imaged, to obtain the essence of rock surface
It is fine into picture；
(3) pore volume segmentation and statistics are carried out to Precise imaging, to obtain pore morphology parameter；
(4) automatic ensaying is carried out to test sample, and acquires its characteristic atomic spectrogram, obtain the distribution of various elements
Spectral line, and then obtain additional conductive mineral content；
(5) saturation degree litho-electric parameters determine: establishing the transformational relation of pore morphology parameter and m and obtain each layer of specific aim
Porosity structure factor m*, establish the transformational relation of additional conductive mineral content and n and obtain each layer of saturation exponent
n*；Wherein, m is the pore structure factor in conventional Archie formula, and n is the saturation degree factor in conventional Archie formula；
(6) step (5) resulting porosity structure factor m* and saturation exponent * is used, containing for complex hydrocarbon layer is calculated
Water saturation Sw。
Further, as follows to the preparation process of test sample in the step (1):
1) pass through drilling shaft lining coring or take chip sample；
2) regular shape is condensed into the sample of step 1) by the easy gelling material of plasticity, and cutting and polishing is carried out to it
Observation surface is formed, to show the rock deposition state of prime stratum.
Further, specific step is as follows for backscattered electron imaging and secondary electron imaging in the step (2):
Observation surface in step 2) is put into electron microscope, adjusts field range and focusing range, selection has representative
Property sample area and carry out the scanning of point-by-point secondary electron and backscattered electron scanning, finally obtain rock surface it is fine at
Picture.
Further, to the solid mineral grain and hole sky in the Precise imaging of step 2) acquisition in the step (3)
Between separated, and interstitial space is divided into larger hole and smaller aperture, smaller gap includes venturi, counts smaller respectively
Every morphological feature of hole is to obtain pore morphology parameter.
Further, specific step is as follows in the step (4):
A. point-by-point automatic ensaying is carried out to the representative region on sample using the energy disperse spectroscopy on electron microscope,
Characteristic atomic spectrogram is acquired and generated, the distribution spectral line of various elements is obtained；
B. using the component of the common various rocks of earth's surface and content data library, element spectrogram inverting is obtained into various mines
Object content；
C. additional conductive mineral content is counted.
Further, specific step is as follows for the step (5):
I, combines typical oil reservoir and typical water layer laboratory core-analysis number by the pore morphology parameter in step (3)
According to, it is established that the transforming function transformation function f between pore morphology parameter and single layer litho-electric parameters m；
II, combines typical oil reservoir and typical water layer experiment by the resulting additional conductive mineral content of statistics in step (4)
Room core analysis data, it is established that the transforming function transformation function F between additional conductive mineral content and single layer litho-electric parameters n；
III, establishes the transformational relation between pore morphology parameter and m: m*=m X f；Establish additional conductive mineral and n it
Between transformational relation: n*=n X F, finally obtaining each layer, targetedly porosity structure factor m* and each layer be targetedly
Saturation exponent *.
The invention has the benefit that
1. improving the evaluation of complex hydrocarbon layer water saturation using a kind of analyze using digital cores provided by the present invention
The method of precision carries out stratum Analysis of micro pore structure and mine by the data using live side-wall coring or chip sample
Object composition analysis extracts distribution and the pass of the fine pore and venturi (being referred to as smaller gap) that influence water saturation evaluation
Key adds the content of electric conductivity mineral, and then obtains accurate quantification water saturation numerical value, overcomes due to during usual coring
By mud flush, mechanical impact etc., lead to not directly carry out seeking for original water saturation using sidewall coring etc., it should
Method is verified by practical oil field data, is significantly improved reservoir saturation degree Explanation Accuracy, can be fast implemented saturation degree and seek, have
Effect directs oil field development production.
2. in step (3) by Precise imaging solid mineral grain and interstitial space separate, be conducive to
Statistics obtains the various features of smaller aperture, pore morphology parameter is ultimately formed, since the various features of smaller aperture are to rock
The content of middle irreducible water and the height of resistivity have a very big impact, therefore, for pore morphology parameter determination to the later period
The correction of porosity structure factor reliable ensure is provided.
3. element spectrogram inverting can be obtained various mineral contents in step (4), due to adding conductive mineral content
Larger to resistivity effects, the correction for being counted additional conductive mineral content to the later period for parameter n provides reliably
It ensures.
Detailed description of the invention
Fig. 1 is the side provided by the invention for being analyzed using digital cores and improving complex hydrocarbon layer water saturation evaluation precision
The flow diagram of method；
Fig. 2 is the side provided by the invention for being analyzed using digital cores and improving complex hydrocarbon layer water saturation evaluation precision
The Precise imaging schematic diagram of the rock surface of step (2) in method；
Fig. 3 is the side provided by the invention for being analyzed using digital cores and improving complex hydrocarbon layer water saturation evaluation precision
The scan image that sample pore volume is divided in the step of method (3)；
Fig. 4 is the side provided by the invention for being analyzed using digital cores and improving complex hydrocarbon layer water saturation evaluation precision
The scan image of the automatic ensaying of sample in the step of method (4).
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.
It should also be noted that similar label and letter indicate similar terms in following attached drawing, therefore, once a certain Xiang Yi
It is defined in a attached drawing, does not then need that it is further defined and explained in subsequent attached drawing.
In the Formula Series using resistivity logging data evaluation reservoir saturation degree, Archie formula is still using most
One of extensive formula.Litho-electric parameters a, b, m and n are the key that whether Archie formula can accurately calculate reservoir saturation degree ginseng
Number, the especially reservoir of the Complicated Pore Structures such as carbonate rock, volcanic rock.According to research shows that n is to calculating if m can be controlled
The influence of reservoir water saturation is conducive to the saturation degree that Complicated Pore Structures reservoir is precisely calculated using Archie formula,
In Archie formulaThe porosity of rock can be obtained by neutron, density and acoustic logging；Rw: the series of strata water flooding
Resistivity measurements；Rt: stratum actual resistivity measured value.
The present invention provides a kind of sides that raising complex hydrocarbon layer water saturation evaluation precision is analyzed using digital cores
Method, comprising the following steps:
(1) test sample is prepared；It is as follows to the preparation process of test sample in the step:
1) pass through drilling shaft lining coring or the sample for taking landwaste the most basic；
2) regular shape is condensed into the sample of step 1) by the easy gelling material of plasticity, passes through cutting machine and polishing machine
Cutting and polishing is carried out to the sample for condensing into regular shape to form observation surface, to show prime stratum on observation surface
Rock deposition state.Preferably, epoxy resin can be selected in the easy gelling material of plasticity.It is carried out using sample of this method to basis
Processing can overcome in existing method due to, by mud flush, mechanical impact etc., leading to not directly during usual coring
The defect that original water saturation is sought is carried out using side-wall coring etc..
(2) backscattered electron is carried out to the test sample of step (1) and secondary electron is imaged；Backscattered electron in the step
It is as follows with the specific processing step of secondary electron imaging:
Observation surface in step 2) is put into electron microscope, adjusts field range and focusing range, selection has representative
Property sample area carry out the scanning of point-by-point secondary electron and backscattered electron scanning, respectively obtain the Precise imaging of rock surface.
Wherein, the secondary electron imaging mainly shape characteristic of reflection sample surfaces 10nm or so, as contrast be pattern
The formation of contrast, contrast is mainly taken at inclination angle of the sample surfaces relative to incident beam.Secondary electron image high resolution, without bright
Aobvious shadow effect, field depth is big, three-dimensional sense is strong, is the main imaging mode of scanning electron microscope, especially suitable for rough samples surface
Morphology observation has a wide range of applications in the fields such as material and life science.
The composition characteristics of the main reflection sample surfaces of backscattered electron imaging, i.e. sample average atomic number Z big position
Stronger backscattered electron signal is generated, forms brighter region on fluorescent screen.
(3) secondary electron is imaged and carries out pore volume segmentation and statistics, to obtain pore morphology parameter；It is right in the step
Solid mineral grain and interstitial space in the Precise imaging of step 2) acquisition are separated, and interstitial space are divided into larger
Hole and smaller aperture, smaller gap include venturi, count the various features of smaller aperture respectively to obtain pore morphology parameter.
Wherein, count smaller aperture every morphological feature mainly include the accounting of smaller aperture, average shape factor, length-width ratio, partially
Degree etc..The various features of smaller aperture have a very big impact the content of irreducible water in rock and the height of resistivity.
(4) to the automatic ensaying of representational sample area progress is chosen in step (2), to acquire its characteristic atomic
Spectrogram；Utilize the energy disperse spectroscopy on electron microscope to representational sample area (with the sample area in step (2) in the step
Domain range is consistent) point-by-point automatic ensaying is carried out, characteristic atomic spectrogram is acquired and generated, the distribution profile of various elements is obtained
Line, to obtain additional conductive mineral content；Specific step is as follows for the step:
A. point-by-point automatic ensaying is carried out to the representative region on sample using the energy disperse spectroscopy on electron microscope,
Characteristic atomic spectrogram is acquired and generated, the distribution spectral line of various elements is obtained；
B. the component of the common various rocks of earth's surface and content data library, the composition member of the common various rocks of earth's surface are utilized
Element and content data library belong to the record data of existing literature, can directly carry out using inverting when calculating various mineral contents,
Element spectrogram inverting is obtained into various mineral contents；
C. additional conductive mineral content is counted, additional conductive mineral those of larger to resistivity effects is especially paid close attention to and contains
Amount, such as total shale content, smectite content.
(5) saturation litho-electric parameters determine: establishing the transformational relation of pore morphology parameter and m and obtain each layer targetedly
Porosity structure factor m* establishes the transformational relation of additional conductive mineral content and n and obtains each layer of saturation exponent *；
Wherein, m is the pore structure factor in conventional Archie formula, and n is the saturation degree factor in conventional Archie formula；The step
Specific step is as follows:
I, combines typical oil reservoir and typical water layer laboratory core-analysis number by the pore morphology parameter in step (3)
According to, it is established that the transforming function transformation function f between pore morphology parameter and single layer litho-electric parameters m；Wherein, m is that skeleton draws with pore network
The measurement of the hole tortuosity risen, hole tortuosity is higher, and m value is bigger, otherwise when rock matrix is the independent part being separated from each other
Timesharing, that is, when hole does not have tortuosity, m=1.The pore morphology factor, venturi of every block of rock can be depicted by step (3)
The tortuosities such as intermediate value characterising parameter (i.e. pore morphology parameter), in conjunction with oil reservoir typical on the well for sampling region and typical water layer
Laboratory core-analysis data, it is established that the transforming function transformation function f between pore morphology parameter and single layer litho-electric parameters m, the transformation letter
Number will be different in different basins, DIFFERENT DEPOSITIONAL ENVIRONMENTS, and then can quickly obtain each layer targetedly porosity factor m*.
II, combines typical oil reservoir and typical water layer experiment by the resulting additional conductive mineral content of statistics in step (4)
Room core analysis data, it is established that the transforming function transformation function F between additional conductive mineral content and single layer litho-electric parameters n；
Wherein, n is saturation exponent, is the correction to saturation degree micro Distribution inhomogeneities, generally by shale content and
The influence of shale patterns affect, especially smectite content, for complicated reservoirs, single value can not then make up and down more set reservoirs with
And the explanation of plane multiple wells is reasonable simultaneously.By the automatic ensaying of step (4), obtains shale in each reservoir samples and contain
The content of amount and montmorillonite, in conjunction with oil reservoir typical on the well for sampling region and typical water layer laboratory core-analysis number
According to, it is established that the transforming function transformation function between shale type and shale content (i.e. additional conductive mineral content) and single layer litho-electric parameters n
F, the transforming function transformation function will be different in different basins, DIFFERENT DEPOSITIONAL ENVIRONMENTS, and then the saturation degree that can quickly obtain each layer refers to
Number n*.
III, establishes the transformational relation between pore morphology parameter and m: m*=m X f；Establish additional conductive mineral and n it
Between transformational relation: n*=n X F, finally obtaining each layer, targetedly porosity structure factor m* and each layer be targetedly
Saturation exponent *.
The determination for why carrying out saturation litho-electric parameters, due to for conventional sandstone, using single litho-electric parameters method base
This is applicable in, but for complicated reservoirs, and single litho-electric parameters can not then make the explanation of more set reservoirs and plane multiple wells up and down same
Shi Heli.To improve saturation degree Explanation Accuracy, need to each layer of reservoir using targetedly litho-electric parameters m* and n*, it is necessary to
Consideration is changed with blowhole tortuosity and shale.
(6) water saturation of complex hydrocarbon layer is calculated, using the litho-electric parameters on stratum are obtained in committed step (5), i.e.,
Each layer targetedly porosity factor m*, each layer of saturation exponent *, and the porosity of rockPass through neutron, density
It is obtained with acoustic logging, a regulation coefficient related with shale ingredient in rock, content and distribution form, generally 1 or so；B is
As the additional saturation degree micro Distribution coefficient of uncertainty of wetability, generally 1 or so and the stratum water power on this set stratum
The resistance rate Rw and resistivity Rt on this set stratum, the water saturation on stratum can be acquired using Archie formula.
The present invention is not limited to above-mentioned optional embodiment, anyone can show that other are various under the inspiration of the present invention
The product of form, however, make any variation in its shape or structure, it is all to fall into the claims in the present invention confining spectrum
Technical solution, be within the scope of the present invention.
Claims (6)
1. a kind of analyze the method for improving complex hydrocarbon layer water saturation evaluation precision using digital cores, which is characterized in that
The following steps are included:
(1) test sample is prepared；
(2) backscattered electron is carried out to the test sample of step (1) and secondary electron is imaged, with obtain rock surface it is fine at
Picture；
(3) pore volume segmentation and statistics are carried out to Precise imaging, to obtain pore morphology parameter；
(4) automatic ensaying is carried out to test sample, and acquires its characteristic atomic spectrogram, obtain the distribution profile of various elements
Line, and then obtain additional conductive mineral content；
(5) saturation degree litho-electric parameters determine: establishing the transformational relation of pore morphology parameter and m and obtain each layer targetedly hole
Porosity structure factor m* establishes the transformational relation of additional conductive mineral content and n and obtains each layer of saturation exponent *；Its
In, m is the pore structure factor in conventional Archie formula, and n is the saturation degree factor in conventional Archie formula；
(6) step (5) resulting porosity structure factor m* and saturation exponent * is used, the aqueous full of complex hydrocarbon layer is calculated
With degree Sw。
2. according to claim 1 analyze the side for improving complex hydrocarbon layer water saturation evaluation precision using digital cores
Method, which is characterized in that as follows to the preparation process of test sample in the step (1):
1) pass through drilling shaft lining coring or take chip sample；
2) regular shape is condensed into the sample of step 1) by the easy gelling material of plasticity, and cutting and polishing is carried out to it and is formed
Surface is observed, to show the rock deposition state of prime stratum.
3. according to claim 1 analyze the side for improving complex hydrocarbon layer water saturation evaluation precision using digital cores
Method, which is characterized in that specific step is as follows for backscattered electron imaging and secondary electron imaging in the step (2):
Observation surface in step 2) is put into electron microscope, adjusts field range and focusing range, is chosen representational
Sample area simultaneously carries out point-by-point secondary electron scanning and backscattered electron scanning, finally obtains the Precise imaging of rock surface.
4. according to claim 1 analyze the side for improving complex hydrocarbon layer water saturation evaluation precision using digital cores
Method, which is characterized in that in the step (3) to step 2) acquisition Precise imaging in solid mineral grain and interstitial space into
Row separation, and interstitial space is divided into larger hole and smaller aperture, smaller gap includes venturi, counts smaller aperture respectively
Every morphological feature to obtain pore morphology parameter.
5. according to claim 1 analyze the side for improving complex hydrocarbon layer water saturation evaluation precision using digital cores
Method, which is characterized in that specific step is as follows in the step (4):
A. point-by-point automatic ensaying, acquisition are carried out to the representative region on sample using the energy disperse spectroscopy on electron microscope
And characteristic atomic spectrogram is generated, obtain the distribution spectral line of various elements；
B. using the component of the common various rocks of earth's surface and content data library, element spectrogram inverting is obtained into various mineral and is contained
Amount；
C. additional conductive mineral content is counted.
6. according to claim 1 analyze the side for improving complex hydrocarbon layer water saturation evaluation precision using digital cores
Method, which is characterized in that specific step is as follows for the step (5):
I, combines typical oil reservoir and typical water layer laboratory core-analysis data by the pore morphology parameter in step (3), builds
Erect the transforming function transformation function f between pore morphology parameter and single layer litho-electric parameters m；
II, combines typical oil reservoir and typical water layer laboratory rock by the resulting additional conductive mineral content of statistics in step (4)
The heart analyzes data, it is established that the transforming function transformation function F between additional conductive mineral content and single layer litho-electric parameters n；
III, establishes the transformational relation between pore morphology parameter and m: m*=m X f；It establishes between additional conductive mineral and n
Transformational relation: n*=n X F, finally obtaining each layer, targetedly porosity structure factor m* and each layer are targetedly saturated
Spend index n*.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811296212.3A CN109444190B (en) | 2018-11-01 | 2018-11-01 | Method for improving evaluation precision of water saturation of complex oil and gas reservoir by utilizing digital core analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811296212.3A CN109444190B (en) | 2018-11-01 | 2018-11-01 | Method for improving evaluation precision of water saturation of complex oil and gas reservoir by utilizing digital core analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109444190A true CN109444190A (en) | 2019-03-08 |
| CN109444190B CN109444190B (en) | 2021-11-09 |
Family
ID=65550626
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811296212.3A Active CN109444190B (en) | 2018-11-01 | 2018-11-01 | Method for improving evaluation precision of water saturation of complex oil and gas reservoir by utilizing digital core analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109444190B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110108739A (en) * | 2019-05-07 | 2019-08-09 | 中国石油化工股份有限公司 | The method for characterizing rough surface rock sample Crude Oil occurrence status |
| CN115234216A (en) * | 2022-05-25 | 2022-10-25 | 东营市永昇能源科技有限责任公司 | Method for establishing logging identification layout and application thereof |
| CN117054460A (en) * | 2023-10-12 | 2023-11-14 | 中国科学院地质与地球物理研究所 | Method for calculating saturation of different fluids based on refrigeration electron microscope |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1719875B1 (en) * | 2000-05-23 | 2010-08-11 | Saudi Arabian Oil Company | Pyrolytic oil-productivity index method for predicting reservoir rock and oil characteristics |
| CN101930082A (en) * | 2009-06-24 | 2010-12-29 | 中国石油集团川庆钻探工程有限公司 | Method for distinguishing reservoir fluid type by adopting resistivity data |
| US20120065951A1 (en) * | 2010-09-02 | 2012-03-15 | Schlumberger Technology Corporation | Method to predict dense hydrocarbon saturations for high pressure high temperature |
| CN103375166A (en) * | 2012-04-26 | 2013-10-30 | 中国石油集团长城钻探工程有限公司 | Method for determining content of a plurality of mineral constituents in stratum |
| CN103422853A (en) * | 2012-05-24 | 2013-12-04 | 中国石油天然气集团公司 | Method and device for determining water saturation of stratum |
| CN103543474A (en) * | 2012-07-16 | 2014-01-29 | 中国石油化工股份有限公司 | Hydrocarbon saturation evaluation method based on non-conductive pore model |
| CN103615230A (en) * | 2013-10-26 | 2014-03-05 | 中国石油化工集团公司 | Method for establishing double-shale-indicative-factor water saturation model |
| CN104712329A (en) * | 2015-01-29 | 2015-06-17 | 中国石油大学(华东) | Computing model of mud shale oil gas saturability |
| CN105334149A (en) * | 2015-11-24 | 2016-02-17 | 中国石油大学(北京) | Micro-pore structure evaluation and reservoir classification method for tight reservoirs |
| CN105464654A (en) * | 2015-12-15 | 2016-04-06 | 中国石油天然气股份有限公司 | Method and device for determining saturation exponent of tight sandstone |
| CN105954496A (en) * | 2016-06-13 | 2016-09-21 | 中国石油天然气股份有限公司 | Reconstruction method for pore throat and fluid distribution of sandstone reservoir |
| CN107389396A (en) * | 2017-06-23 | 2017-11-24 | 东北石油大学 | Realize the preparation method and its experimental method of heterogeneous core in the layer that dispensing point is adopted |
| CN107780923A (en) * | 2017-11-01 | 2018-03-09 | 中石化石油工程技术服务有限公司 | A kind of foundation of the water-saturation model based on Shale Correction, emulation mode |
-
2018
- 2018-11-01 CN CN201811296212.3A patent/CN109444190B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1719875B1 (en) * | 2000-05-23 | 2010-08-11 | Saudi Arabian Oil Company | Pyrolytic oil-productivity index method for predicting reservoir rock and oil characteristics |
| CN101930082A (en) * | 2009-06-24 | 2010-12-29 | 中国石油集团川庆钻探工程有限公司 | Method for distinguishing reservoir fluid type by adopting resistivity data |
| US20120065951A1 (en) * | 2010-09-02 | 2012-03-15 | Schlumberger Technology Corporation | Method to predict dense hydrocarbon saturations for high pressure high temperature |
| CN103375166A (en) * | 2012-04-26 | 2013-10-30 | 中国石油集团长城钻探工程有限公司 | Method for determining content of a plurality of mineral constituents in stratum |
| CN103422853A (en) * | 2012-05-24 | 2013-12-04 | 中国石油天然气集团公司 | Method and device for determining water saturation of stratum |
| CN103543474A (en) * | 2012-07-16 | 2014-01-29 | 中国石油化工股份有限公司 | Hydrocarbon saturation evaluation method based on non-conductive pore model |
| CN103615230A (en) * | 2013-10-26 | 2014-03-05 | 中国石油化工集团公司 | Method for establishing double-shale-indicative-factor water saturation model |
| CN104712329A (en) * | 2015-01-29 | 2015-06-17 | 中国石油大学(华东) | Computing model of mud shale oil gas saturability |
| CN105334149A (en) * | 2015-11-24 | 2016-02-17 | 中国石油大学(北京) | Micro-pore structure evaluation and reservoir classification method for tight reservoirs |
| CN105464654A (en) * | 2015-12-15 | 2016-04-06 | 中国石油天然气股份有限公司 | Method and device for determining saturation exponent of tight sandstone |
| CN105954496A (en) * | 2016-06-13 | 2016-09-21 | 中国石油天然气股份有限公司 | Reconstruction method for pore throat and fluid distribution of sandstone reservoir |
| CN107389396A (en) * | 2017-06-23 | 2017-11-24 | 东北石油大学 | Realize the preparation method and its experimental method of heterogeneous core in the layer that dispensing point is adopted |
| CN107780923A (en) * | 2017-11-01 | 2018-03-09 | 中石化石油工程技术服务有限公司 | A kind of foundation of the water-saturation model based on Shale Correction, emulation mode |
Non-Patent Citations (6)
| Title |
|---|
| 刘伟 等: "数字岩心技术在致密砂岩储层含油饱和度评价中的应用", 《断块油气田》 * |
| 孙游雪 等: "水淹层评价中Ｒｗ、ｍ、ｎ 值的确定", 《山东理工大学学报（自然科学版）》 * |
| 李先鹏: "泥质地层中饱和度指数的确定", 《岩性油气藏》 * |
| 杜园青: "复杂断块油藏流体识别主控因素分析及解释方法研究", 《长江大学学报（自科版）》 * |
| 殷树军: "凝灰质储层复杂油水层测井识别方法", 《石油天然气学报(江汉石油学院学报)》 * |
| 胡英杰 等: "吉林油田低阻油气层岩心实验研究与解释", 《河南石油》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110108739A (en) * | 2019-05-07 | 2019-08-09 | 中国石油化工股份有限公司 | The method for characterizing rough surface rock sample Crude Oil occurrence status |
| CN115234216A (en) * | 2022-05-25 | 2022-10-25 | 东营市永昇能源科技有限责任公司 | Method for establishing logging identification layout and application thereof |
| CN115234216B (en) * | 2022-05-25 | 2023-09-19 | 东营市永昇能源科技有限责任公司 | Method for establishing logging identification layout and application thereof |
| CN117054460A (en) * | 2023-10-12 | 2023-11-14 | 中国科学院地质与地球物理研究所 | Method for calculating saturation of different fluids based on refrigeration electron microscope |
| CN117054460B (en) * | 2023-10-12 | 2023-12-08 | 中国科学院地质与地球物理研究所 | Method for calculating saturation of different fluids based on refrigeration electron microscope |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109444190B (en) | 2021-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107228816B (en) | Different type pore evolution evaluation method in a kind of mud shale | |
| Anselmetti et al. | Quantitative characterization of carbonate pore systems by digital image analysis | |
| US10466386B2 (en) | Methods for determining a saturation-height function in oil and gas reservoirs | |
| CN109444190A (en) | The method for improving complex hydrocarbon layer water saturation evaluation precision is analyzed using digital cores | |
| CN104698500A (en) | Method for predicting reservoir lithogenous phase through geology and logging information | |
| CN109030311A (en) | Based on nuclear magnetic resonance T2Compose the pore structure classification and recognition methods of sensitive parameter | |
| CN104989386B (en) | A kind of generation method of fine and close oily property identification method and plate | |
| CN108956424A (en) | A kind of method of hole quantitatively characterizing in shale | |
| Yamada et al. | Revisiting porosity analysis from electrical borehole images: integration of advanced texture and porosity analysis | |
| CN107003246A (en) | The fluid behaviour of porous material LIBS | |
| CN106501870B (en) | A kind of lacustrine facies densification shelly limestone comparative good-quality reservoir stratum identification method | |
| CN114897767A (en) | Multi-scale representation and reservoir classification method for reservoir space of compact mixed-reservoir rock reservoir | |
| CN113075731B (en) | Deep reservoir continuity wellbore digital modeling method and device | |
| CN110470584A (en) | A method of evaluation imbibition and water lock comprehensive effect | |
| CN110019594A (en) | A kind of geologic datebase method for building up based on numerical simulation | |
| CN103499594A (en) | Method for measuring lower limit of feature size of crude oil movable pore throat | |
| Zhang et al. | Characterization of the pore size distribution with SEM images processing for the tight rock | |
| CN109580678B (en) | Method for rapidly identifying and evaluating low-resistance oil-gas reservoir by using digital core technology | |
| Kazak et al. | Integration of large-area Scanning-Electron-Microscopy imaging and automated mineralogy/petrography data for selection of nanoscale pore-space characterization sites | |
| Gerard et al. | Petrographic image analysis: an alternate method for determining petrophysical properties | |
| CN117233064B (en) | Method for calculating shale nuclear magnetic resonance logging fracture cut-off value | |
| CN117129509B (en) | Method for calculating shale fracture nuclear magnetic resonance logging T2 cut-off value based on KI-CT | |
| Tawfeeq et al. | Digital rock analysis: an alternative method to predict petrophysical properties, case study from Mishrif Formation | |
| CN113137225A (en) | Gravel reservoir evaluation method based on gravel component analysis | |
| CN110399649A (en) | A kind of metamorphic rock quantization recognition methods based on diagenesis indicator element |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |