CN111562629B - Saturation determination method and device based on equivalent pore section index - Google Patents
Saturation determination method and device based on equivalent pore section index Download PDFInfo
- Publication number
- CN111562629B CN111562629B CN201910112284.6A CN201910112284A CN111562629B CN 111562629 B CN111562629 B CN 111562629B CN 201910112284 A CN201910112284 A CN 201910112284A CN 111562629 B CN111562629 B CN 111562629B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- resistivity
- index
- equivalent
- 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.)
- Active
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 277
- 238000000034 method Methods 0.000 title claims abstract description 93
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 42
- 239000008398 formation water Substances 0.000 claims description 107
- 239000011435 rock Substances 0.000 claims description 64
- 238000004364 calculation method Methods 0.000 claims description 43
- 238000004457 water analysis Methods 0.000 claims description 22
- 238000011010 flushing procedure Methods 0.000 claims description 21
- 238000004590 computer program Methods 0.000 claims description 17
- 229920006395 saturated elastomer Polymers 0.000 claims description 17
- 230000015572 biosynthetic process Effects 0.000 claims description 15
- 238000003860 storage Methods 0.000 claims description 14
- 239000000047 product Substances 0.000 description 22
- 238000005755 formation reaction Methods 0.000 description 12
- 238000002474 experimental method Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 10
- 230000008859 change Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 9
- 238000004891 communication Methods 0.000 description 8
- 238000011158 quantitative evaluation Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 5
- 230000007423 decrease Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 210000004884 grey matter Anatomy 0.000 description 3
- 230000035515 penetration Effects 0.000 description 3
- 230000035699 permeability Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a saturation determination method and device based on equivalent pore section indexes, wherein the method comprises the following steps: obtaining an equivalent pore section index; calculating a pore structure index using the equivalent pore section index; and determining the water saturation of the reservoir by using the porosity structural index. The invention can establish a method and a device for accurately calculating the saturation by simultaneously considering the porosity and the pore structure in a complex reservoir.
Description
Technical Field
The invention relates to the field of petroleum exploration, in particular to a geophysical exploration technology, and particularly relates to a saturation calculation method and device based on equivalent pore section indexes.
Background
With the increasing complexity of oil and gas exploration and development objects, complex reservoirs are increasing, such as: carbonate, mudstone fracture reservoirs, volcanic reservoirs and the like, the difficulty in well logging evaluation of reservoir saturation is continuously increased, so that the well logging interpretation coincidence rate is low, the oil gas discovery rate is low, and the reason is the problem of selecting a saturation model and parameters thereof. Currently, a plurality of known reservoir saturation evaluation models are provided, and a large number of Archie models and derivative forms thereof are formed in the reservoir saturation evaluation process by using a main Alqi formula. The most critical of these is how to accurately determine the petroelectricity parameters in the alchi model: lithology coefficient a, saturation coefficient b, pore structure index m and cementation index n.
The previous experimental research results and a large number of exploration and development practices show that for most reservoirs, including complex sand shale reservoirs, the values of a and b are more than 1, and the values of n are more than 2; the most studied and the biggest controversy is the determination of the pore structure index m value. The most common method for determining the parameter is to use a rock core to carry out a rock electric experiment, but the method has the problem that the scales of the rock core and the well logging are inconsistent, is limited by the number of cores and the representative problem of the rock core, and the obtained rock electric parameter is mostly fixed, and the actual reservoir pore structure is complex and changeable, so that the method is difficult to objectively and continuously reflect the real situation of the reservoir. Later, well logging students at home and abroad put forward a method for changing the m value of the pore structure, the method considers the influence of the porosity of the reservoir on the m value, so that the m values of the pore structures of reservoirs with different porosities are different, the purpose of changing the m value of the pore structure is achieved, and a good application effect is achieved, but the method only considers the influence of the porosity on the m value, and does not consider the influence of the pore structure of the reservoir on the m value, but the influence is the most important, and the result that the interpretation result of the saturation of the reservoir of complex sand shale is inconsistent with the test conclusion is brought.
Therefore, how to provide a method and a device for accurately calculating the oil saturation in a complex reservoir by considering the porosity and the pore structure at the same time is a problem to be solved.
Disclosure of Invention
Aiming at the problems in the prior art, the invention can establish the method and the device for accurately calculating the oil saturation by considering the porosity and the pore structure in the complex reservoir. The saturation calculation precision and interpretation coincidence rate are greatly improved, so that the purpose of calculating the complex reservoir saturation with high precision of the conventional logging series is achieved, and the problem of quantitative evaluation of the reservoir saturation is effectively solved.
In order to solve the technical problems, the invention provides the following technical scheme:
in a first aspect, the present invention provides a saturation determination method based on equivalent pore section index, comprising:
obtaining an equivalent pore section index;
calculating a pore structure index by using the equivalent pore section index;
reservoir water saturation is determined using the porosity structural index.
In one embodiment, obtaining the equivalent pore cross-section index comprises:
creating an equivalent rock volume model;
and obtaining an equivalent pore section index by utilizing an equivalent rock volume model according to the parallel resistance principle.
In one embodiment, the obtaining the equivalent pore section index using the equivalent rock volume model according to the principle of parallel resistance includes:
Acquiring logging curve data;
acquiring formation water resistivity, reservoir total porosity and reservoir flushing zone resistivity according to logging curve data;
and calculating the equivalent pore section index according to the product of the formation water resistivity and the total porosity of the reservoir and the resistivity of the reservoir flushing zone.
In one embodiment, the calculation formula for calculating the equivalent pore section index from the product of the formation water resistivity, the total porosity of the reservoir, and the reservoir flushing zone resistivity is:
wherein, IPS is equivalent pore section index, dimensionless; r is R w Is the formation water resistivity, Ω·m; phi is the total porosity of the reservoir,%; r is R xo The reservoir is flushed with resistivity, Ω·m.
In one embodiment, the obtaining the equivalent pore section index using the equivalent rock volume model according to the principle of parallel resistance includes:
acquiring rock core experimental data, formation water analysis data and logging curve data;
acquiring formation water resistivity, total porosity of the reservoir and reservoir resistivity when the reservoir is saturated with formation water by 100% according to rock core experimental data, formation water analysis data and logging curve data;
and calculating the equivalent pore section index according to the product of the formation water resistivity and the total porosity of the reservoir and the reservoir resistivity when the reservoir is saturated with formation water by 100 percent.
In one embodiment, the calculation formula for calculating the equivalent pore section index from the product of the formation water resistivity, the total porosity of the reservoir, and the reservoir resistivity at 100% saturation of the reservoir with formation water is:
wherein, IPS is equivalent pore section index, dimensionless; r is R w Is the formation water resistivity, Ω·m; phi is the total porosity of the reservoir,%; r is R 0 Reservoir when the reservoir is 100% saturated with formation waterResistivity, Ω·m.
In one embodiment, calculating the pore structure index using the equivalent pore cross section index comprises: the pore structure index formula calculated by using the equivalent pore section index is as follows:
m＝-logIPS
wherein m is a pore structure index, and is dimensionless.
In one embodiment, the pore structure index m has a value in the range of 0 to 3.
In one embodiment, determining reservoir water saturation using the porosity structural index includes:
and calculating the water saturation of the reservoir according to the product of the formation water resistivity, the lithology coefficient and the saturation coefficient, the reservoir undisturbed formation resistivity, the cementation index, the reservoir total porosity and the pore structure index.
In one embodiment, the calculation formula for calculating the reservoir water saturation according to the formation water resistivity, the product of the lithology coefficient and the saturation coefficient, the reservoir undisturbed formation resistivity, the cementation index, the reservoir total porosity and the pore structure index is as follows:
Wherein, a: lithology coefficient, dimensionless; b: saturation coefficient, dimensionless; n: cementing index, dimensionless; sw: reservoir water saturation,%; r is R t : the resistivity of the undisturbed stratum of the reservoir is omega-m.
In a second aspect, the present invention provides a saturation determination apparatus based on equivalent pore section index, the apparatus comprising:
the equivalent pore section index obtaining unit is used for obtaining the equivalent pore section index;
a pore structure index calculation unit for calculating a pore structure index using the equivalent pore section index;
and the water saturation determining unit is used for determining the water saturation of the reservoir by using the porosity structural index.
In one embodiment, the equivalent pore section index obtaining unit includes:
the model creation module is used for creating an equivalent rock volume model;
and the equivalent pore section index acquisition module is used for acquiring the equivalent pore section index by utilizing an equivalent rock volume model according to the parallel resistance principle.
In one embodiment, the equivalent pore section index obtaining module includes:
logging curve data acquisition module: the method comprises the steps of acquiring log curve data;
a logging curve data calculation module: the method comprises the steps of obtaining formation water resistivity, reservoir total porosity and reservoir flushing zone resistivity according to logging curve data;
Equivalent pore section index first calculation module: and the method is used for calculating the equivalent pore section index according to the product of the formation water resistivity and the total porosity of the reservoir and the resistivity of the reservoir flushing zone.
In one embodiment, the equivalent pore section index obtaining module includes:
rock core experimental data and formation water analysis data acquisition module: the method comprises the steps of acquiring rock core experimental data, stratum water analysis data and logging curve data;
rock core experimental data and formation water analysis data calculation module: the method comprises the steps of obtaining formation water resistivity, total porosity of a reservoir and reservoir resistivity when the reservoir is saturated with formation water by 100% according to rock core experimental data, formation water analysis data and logging curve data;
the equivalent pore section index calculation module is a second calculation module: and the method is used for calculating the equivalent pore section index according to the product of the formation water resistivity and the total porosity of the reservoir and the reservoir resistivity when the reservoir is saturated with formation water by 100 percent.
In an embodiment, the water saturation determination unit is specifically configured to: and calculating the water saturation of the reservoir according to the product of the formation water resistivity, the lithology coefficient and the saturation coefficient, the reservoir undisturbed formation resistivity, the cementation index, the reservoir total porosity and the pore structure index.
In a third aspect, the invention provides an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the steps of a saturation determination method based on equivalent pore cross-section index when the program is executed by the processor.
In a fourth aspect, the present invention provides a computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of a saturation determination method based on equivalent pore cross section index.
As can be seen from the above description, the method and the device for determining saturation based on equivalent pore section indexes are provided in the present invention, firstly, on the basis of an equivalent rock volume model, the equivalent pore section indexes are calculated by using two methods of a core experiment and a conventional logging curve, the equivalent pore section indexes can reflect the core pore structure, then, the equivalent pore section indexes are used to define a pore structure index, the pore structure index simultaneously considers the porosity of the reservoir and the pore structure of the reservoir, and finally, the saturation is calculated by using the pore structure index and the alchi formula. The saturation calculation precision and interpretation coincidence rate are greatly improved, so that the purpose of calculating the complex reservoir saturation with high precision of the conventional logging series is achieved, and the problem of quantitative evaluation of the reservoir saturation is effectively solved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a saturation determination method based on equivalent pore cross-section index in an embodiment of the present invention;
FIG. 2 is a flow chart of step 100 in an embodiment of the invention;
FIG. 3 is a flow chart of step 102 in an embodiment of the invention;
FIG. 4 is a flowchart of step 102 in another embodiment of the present invention;
FIG. 5 is a flow chart of a saturation determination method based on equivalent pore section index in an embodiment of the present invention;
FIGS. 6a and 6b are schematic views of equivalent rock volume models based on the equivalent pore cross-section index saturation determination method in an embodiment of the present invention;
FIG. 7 is a flow chart of a saturation determination method based on equivalent pore section index in an embodiment of the present invention;
FIG. 8 is a schematic diagram of the result of explaining a well H using a saturation determination method based on equivalent pore section index in an embodiment of the present invention;
FIG. 9 is a schematic structural diagram of a saturation determination device based on equivalent pore cross-section index in an embodiment of the present invention;
fig. 10 is a schematic structural diagram of an electronic device in an embodiment of the invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
An embodiment of the present invention provides a specific implementation manner of a saturation determining method based on an equivalent pore section index, referring to fig. 1, the method specifically includes the following contents:
step 100: and obtaining the equivalent pore section index.
It is understood that the equivalent pore section index may reflect the core pore structure.
Step 200: and calculating the pore structure index by using the equivalent pore section index.
It will be appreciated that step 200 is to calculate a porosity structural index using the porosity section index calculated in step 100, and to give a variation range of the porosity structural index in combination with actual data.
Step 300: reservoir water saturation is determined using the porosity structural index.
Step 300 may be: obtaining a lithology coefficient, a saturation coefficient and a cementation index through a rock core experiment; and obtaining the resistivity of the undisturbed stratum of the reservoir and the total porosity of the reservoir through logging data, and calculating the water saturation of the core by using the Archie formula and the equivalent pore section index calculated in the step 200, so that the oil saturation of the core can be calculated.
As can be seen from the above description, the present invention provides a saturation determining method based on an equivalent pore section index, by obtaining an equivalent pore section index capable of reflecting a pore structure of a core, and defining a pore structure index by the equivalent pore section index, the pore structure index considers both the porosity of a reservoir and the pore structure of the reservoir, and finally calculates the saturation according to the pore structure index and the alchi formula. The saturation calculation precision and interpretation coincidence rate are greatly improved, so that the purpose of calculating the complex reservoir saturation with high precision of the conventional logging series is achieved, and the problem of quantitative evaluation of the reservoir saturation is effectively solved.
In one embodiment, referring to fig. 2, step 100 comprises:
step 101: creating an equivalent rock volume model;
the reservoir rock is divided into two parts of stratum water in a core skeleton and a core pore in volume.
Step 102: and obtaining an equivalent pore section index by utilizing an equivalent rock volume model according to the parallel resistance principle.
Based on the step 101, the reservoir rock Dan Dianzu is calculated by using the rock skeleton and the stratum water resistance through the resistance parallel principle, and the equivalent pore section index is obtained.
In one embodiment, referring to fig. 3, step 102 comprises:
step 1021: acquiring logging curve data;
in one embodiment, the log data is conventional log data comprising: three lithology log data, three porosity log data, and three electrical log data, in one embodiment, the log data may be: gamma curve, borehole diameter curve and natural potential curve (three lithology log data); neutron porosity curve, sonic jet time log, density log (three-porosity log data); deep lateral resistivity, shallow lateral resistivity curve, and microsphere-type focused resistivity (three-electrical log data).
Step 1022: acquiring formation water resistivity, reservoir total porosity and reservoir flushing zone resistivity according to logging curve data;
in one embodiment, the formation water resistivity may be calculated using a natural potential curve, or may be calculated using a visual formation water resistivity method and conventional log data, although the invention is not limited thereto. Reservoir total porosity may be calculated using one or more combinations of three porosity curves and reservoir flushing zone resistivity may be calculated using a shallow lateral resistivity curve.
Step 1023: and calculating the equivalent pore section index according to the product of the formation water resistivity and the total porosity of the reservoir and the resistivity of the reservoir flushing zone.
In one embodiment, in step 1023, the calculation formula for calculating the equivalent pore section index according to the product of the formation water resistivity and the total porosity of the reservoir and the resistivity of the reservoir flushing zone may be:
wherein IPS is equivalent pore section index,dimensionless; r is R w Is the formation water resistivity, Ω·m; phi is the total porosity of the reservoir,%; r is R xo The reservoir is flushed with resistivity, Ω·m.
In another embodiment, referring to fig. 4, step 102 comprises:
It will be appreciated that the core experimental data in step 102a includes: core resistivity, core skeleton resistivity.
In one embodiment, the step 102a core experimental data further includes: the lithology coefficient, saturation coefficient and cementation index may be the following steps: the determination of reservoir water saturation provides a more accurate parameter.
It will be appreciated that formation water resistivity may be obtained from formation water analysis data, reservoir resistivity at 100% saturation of the reservoir with formation water may be obtained from core experimental data, and reservoir total porosity may be calculated from log data in a manner similar to the reservoir total porosity calculation method of step 1022.
In one embodiment, the calculation formula for calculating the equivalent pore section index in step 102c based on the product of the formation water resistivity and the total porosity of the reservoir and the reservoir resistivity when the reservoir is 100% saturated with formation water may be:
Wherein, IPS is equivalent pore section index, dimensionless; r is R w Is the formation water resistivity, Ω·m; phi isReservoir total porosity,%; r is R 0 Is the reservoir resistivity at 100% saturation of the reservoir with formation water, Ω·m.
In one embodiment, the calculation formula for calculating the pore structure index using the equivalent pore section index in step 200 may be:
m＝-logIPS
wherein m is a pore structure index, and is dimensionless.
It will be appreciated that the pore structure index m calculated in step 200 considers both reservoir porosity and reservoir pore structure.
In one embodiment, step 200 calculates the pore structure index m to have a value ranging from 0 to 3.
It will be appreciated that the value of m itself should have certain constraints, the upper limit of which should generally not be greater than 3. In argillite sandstone and siltstone, due to the narrow pore roar and low permeability, at least in the medium-high pore formations, argillite sandstone and siltstone have higher m values relative to the pure formation. When sandstone contains grey matter, pore penetration decreases and its m tends to increase.
In one embodiment, step 300 may be: and calculating the water saturation of the reservoir according to the product of the formation water resistivity, the lithology coefficient and the saturation coefficient, the reservoir undisturbed formation resistivity, the cementation index, the reservoir total porosity and the pore structure index.
In one embodiment, the calculation formula for calculating the reservoir water saturation in step 300 according to the formation water resistivity, the product of the lithology coefficient and the saturation coefficient, the reservoir undisturbed formation resistivity, the cementation index, the reservoir total porosity and the pore structure index may be:
wherein, a: lithology coefficient, dimensionless; b: saturation coefficient, dimensionless; n: cementing index, dimensionless; sw: reservoir water saturation,%; r is R t : the resistivity of the undisturbed stratum of the reservoir is omega-m.
In one embodiment, the present invention also provides an embodiment of the saturation determination method based on equivalent pore cross-section index, see fig. 5.
Step M01: an equivalent rock volume model is created.
The method comprises the following steps: the equivalent rock volume model consists of two parts of stratum water in a core skeleton and a core pore.
Step M02: and obtaining an equivalent pore section index by utilizing an equivalent rock volume model according to the parallel resistance principle.
It can be understood that, in this step, as shown in fig. 6a and 6b, the resistance of the whole core is regarded as the parallel connection of the core skeleton and the resistance of the pore fluid in the core by using an equivalent rock volume model, which is specifically as follows:
in formula (1): r is (r) 0 、r ma And r w And respectively representing the resistances of the rock core, the framework and formation water in the pores of the rock core, and omega. It should be noted that the equivalent volume model may be based on a core at 100% saturation with formation water.
Since the core skeleton resistance tends to infinity, equation (1) is written as:
the electrical resistance is defined by a value of the electrical resistance,
wherein R is 0 Is the resistivity of the core, omega.m; l is the length of the core and is meter; a is the cross-sectional area of the core and square meters.
The same principle is as follows:
wherein R is w Is the formation water resistivity, Ω·m; l (L) w The length of the core is meter; a is that w Is the cross-sectional area of pores in the core (i.e., the cross-sectional area of formation water in the core), square meters; bringing equations (3) and (4) into equation (2) has:
in the formula (5) of the present invention,
a defined formula for porosity; so equation (5) can be changed to:
in the formula (6), phi t Is the total porosity of the core;
equivalent pore section index is defined as:
in the formula (7), IPS epsilon (0, 1), thereby reflecting the change rule of the pore structure index m value.
From the above description, it can be seen that, in this embodiment, steps M01 and M02 can calculate, based on the equivalent rock volume model, an equivalent pore section index by using the formation water resistivity, the core total porosity and the core resistivity when the core is saturated with formation water by 100%, where the equivalent pore section index can reflect the core pore structure.
M03: and calculating the pore structure index by using the equivalent pore section index.
In the step M03, from the form of the formula (7) of the equivalent pore section index, when the porosity of the core is larger and the resistivity is lower, the equivalent pore section index IPS is larger; when the pore structure of the core is poor, the porosity is reduced, and the higher the resistivity is, the smaller the equivalent pore section index IPS is.
Taking the logarithm of the two sides of the formula (7), the change range of IPS is (- ≡0), and then letting:
m＝-log IPS (8)
according to the rock electricity experiment, the m value of the eastern part (Daqing, dagang, victory, liaohe and Jiang Han) of China is in the range of 1.5-3. According to the experimental data of the mexico core, the m value of the sandstone is 0.5-2.6; in order to make the equivalent pore section index accurately reflect the change of m value, according to the numerical simulation result of domestic researchers, the numerical change range of the formula (8) should be controlled between 0 and 6. The value of m should itself have certain constraints, and its upper limit should generally not be greater than 3.0. In argillite sandstones and siltstone sandstones, the permeability is low due to the narrow pore roar, so that at least in medium-high pore reservoirs, argillite sandstones and siltstone sandstones have a higher m value than in pure reservoirs. When sandstone contains grey matter, pore penetration decreases and its m tends to increase.
Step M04: and calculating the water saturation of the reservoir according to the product of the formation water resistivity, the lithology coefficient and the saturation coefficient, the reservoir undisturbed formation resistivity, the cementation index, the reservoir total porosity and the pore structure index.
The specific method of step M04 may be: the value of the pore section index M calculated according to the step M03, and R obtained from water analysis data w The total water saturation of the reservoir Swt may be calculated in combination with the porosity log and the deep resistivity log in the log. The following Alqi formula may be used to calculate:
in the formula (9), swt is the water saturation of the core,%; r is R t The true resistivity of the stratum is shown, omega-m, n are cementation indexes, and dimensionless; the core oil and gas saturation Sog was then calculated according to Sog =1-Swt.
From the above description, the invention provides a saturation determining method based on equivalent pore section indexes, firstly, on the basis of an equivalent rock volume model, equivalent pore section indexes are calculated by using two methods of a core experiment and a conventional logging curve, the equivalent pore section indexes can reflect a core pore structure, then, a pore structure index is defined by the equivalent pore section indexes, the pore structure index simultaneously considers the porosity of a reservoir and the pore structure of the reservoir, and finally, the saturation is calculated by using the pore structure index and an Archie formula. The saturation calculation precision and interpretation coincidence rate are greatly improved, so that the purpose of calculating the complex reservoir saturation with high precision of the conventional logging series is achieved, and the problem of quantitative evaluation of the reservoir saturation is effectively solved.
To further illustrate the present solution, the present invention provides a specific application example of the saturation determining method based on the equivalent pore section index, taking a certain oilfield well H as an example, where the specific application example specifically includes the following, see fig. 7.
And (I) obtaining an equivalent pore section index.
S0: and obtaining an equivalent rock volume model.
Referring to fig. 6a, the equivalent rock volume model consists of the core skeleton and formation water in the core pores. It should be noted that the equivalent rock volume model is based on a core at 100% saturation with formation water.
S1: logging curve data is acquired.
It will be appreciated that the log data is conventional log data, including: three lithology log data, three porosity log data, and three electrical log data.
S2: and calculating the formation water resistivity by using a apparent formation water resistivity method according to the resistivity logging curve data.
It will be appreciated that S2 may also use natural potential curve data or other methods to calculate formation water resistivity.
S3: and calculating the total porosity of the reservoir according to the acoustic time difference logging curve data.
It will be appreciated that S3 may also be calculated using one or more combinations of three porosity curves.
S4: and calculating the resistivity of the reservoir flushing zone according to the shallow lateral resistivity curve data.
S5: and calculating the equivalent pore section index according to the product of the formation water resistivity and the total porosity of the reservoir and the resistivity of the reservoir flushing zone.
The calculation formula in S5 may be:
wherein, IPS is equivalent pore section index, dimensionless; r is R w Is the formation water resistivity, Ω·m; phi is the total porosity of the reservoir,%; r is R xo The reservoir is flushed with resistivity, Ω·m.
From the previous steps 102a to 102c, steps S1 to S5 may be replaced with:
s1': acquiring core resistivity, core skeleton resistivity, lithology coefficient, saturation coefficient, cementing index, formation water analysis data and logging curve data;
it is understood that step S1' refers to the core resistivity at 100% saturation of the core with formation water, i.e. the reservoir resistivity at 100% saturation of the reservoir with formation water.
S2': and calculating the formation water resistivity according to the formation water analysis data.
S3': and calculating the total porosity of the reservoir according to neutron porosity curve data.
It will be appreciated that the principle of step S3' is similar to step S3.
S4': and calculating the equivalent pore section index according to the product of the formation water resistivity and the total porosity of the reservoir and the reservoir resistivity when the reservoir is saturated with formation water by 100 percent.
The calculation formula in S4' may be:
wherein IPS is equivalent pore section indexDimensionless; r is R w Is the formation water resistivity, Ω·m; phi is the total porosity of the reservoir,%; r is R 0 Is the reservoir resistivity at 100% saturation of the reservoir with formation water, Ω·m.
It can be understood that the steps S1 'to S4' are to use an equivalent rock volume model to consider the resistance of the whole core as the parallel connection of the core skeleton and the resistance of the pore fluid in the core, specifically as shown in the formula (1). Since the core skeleton resistance tends to infinity, formula (1) is written as formula (2), and in addition, the definition of the resistance is defined by formulas (3) and (4), formulas (3) and (4) are brought into formula (2) and are combined with the definition of the porosity to obtain a formula for calculating the equivalent pore section index, as shown in formula (7).
It will be appreciated that: the equivalent pore section indexes IPS calculated in the steps S1 to S5 are more generalized and more convenient to operate, when the core experimental data are not included, the method has a wide application range, but the calculation accuracy is relatively slightly poor, the equivalent pore section indexes IPS calculated in the steps S1 'to S4' are calculated based on the rock experimental data, the result is more accurate, certain requirements on the core data are met in terms of quantity and quality, certain limitation is met, and different calculation methods can be selected by an oil field according to the self situation.
And (II) calculating the pore structure index by using the equivalent pore section index.
S6: and calculating the pore structure index by using the equivalent pore section index.
In step S6, from the form of equation (7) of the equivalent pore section index, when the porosity of the core is larger and the resistivity is lower, the equivalent pore section index IPS is larger; when the pore structure of the core is poor, the porosity is reduced, and the higher the resistivity is, the smaller the equivalent pore section index IPS is.
Taking the logarithm of the two sides of the formula (7), the variation range of IPS is (- ≡0), and the formula for calculating the pore structure index can be obtained as shown in the formula (8).
According to the rock electricity experiment, the m value of the eastern part (Daqing, dagang, victory, liaohe and Jiang Han) of China is in the range of 1.5-3. According to the experimental data of the mexico core, the m value of the sandstone is 0.5-2.6; in order to make the equivalent pore section index accurately reflect the change of m value, according to the numerical simulation result of domestic researchers, the numerical change range of the formula (8) should be controlled between 0 and 6. The value of m should itself have certain constraints, and its upper limit should generally not be greater than 3.0. In argillite sandstone and siltstone, due to the narrow pore roar and low permeability, at least in the medium-high pore formations, argillite sandstone and siltstone have higher m values relative to the pure formation. When sandstone contains grey matter, pore penetration decreases and its m tends to increase.
And thirdly, determining the water saturation of the reservoir by using the porosity structural index.
S7: and determining the water saturation of the reservoir by using a porosity structural index and an Archie formula.
From the value of the equivalent pore section index m calculated in step S6, and R calculated in S w The core total water saturation Swt can be calculated in combination with the porosity log and the deep resistivity log in the log. The core water saturation may be calculated using the alchi equation (equation 9), and the core oil and gas saturation Sog may be calculated according to Sog =1-Swt.
The interpretation results of well H are shown in fig. 8, in interpretation layer No. 1 in the first wire frame and interpretation layer No. 2 in the second wire frame, when the tight layer sandwiched in the gas layer causes the porosity to decrease and the resistivity to increase significantly, i.e. the sandstone section sandwiched by siltstone at 1987-1988 m in layer No. 1, the corresponding depth section sandwiched by siltstone in interpretation layer No. 2, such as: 2001-2005 m, 2005.9-2006.8 m and 2015.7-2016.4 m, the saturation of the hydrocarbon gas calculated by the present method (trace 5, SW2 curve) is correspondingly low, while the saturation calculated by the conventional method (trace 5, water saturation curve) reflects little. The equivalent pore section index comprehensively reflects the change of the porosity and the resistivity, and can better characterize the influence of the change of the pore structure in the reservoir on the saturation of the oil-gas.
It will be appreciated that the saturation calculated by this method, in combination with conventional saturation calculations, can also intuitively indicate the gas-bearing thickness and distribution characteristics of the effective reservoir, and can be used for rapid compartmentalization of the effective reservoir.
From the above description, the invention provides a saturation determining method based on equivalent pore section indexes, firstly, on the basis of an equivalent rock volume model, equivalent pore section indexes are calculated by using two methods of a core experiment and a conventional logging curve, the equivalent pore section indexes can reflect a core pore structure, then, a pore structure index is defined by the equivalent pore section indexes, the pore structure index simultaneously considers the porosity of a reservoir and the pore structure of the reservoir, and finally, the saturation is calculated by using the pore structure index and an Archie formula. The saturation calculation precision and interpretation coincidence rate are greatly improved, so that the purpose of calculating the complex reservoir saturation with high precision of the conventional logging series is achieved, and the problem of quantitative evaluation of the reservoir saturation is effectively solved.
Based on the same inventive concept, the embodiments of the present application also provide a saturation determination device based on an equivalent pore cross-section index, which can be used to implement the method described in the above embodiments, as described in the following embodiments. Since the principle of solving the problem by the saturation determining device based on the equivalent pore section index is similar to that of the saturation determining method based on the equivalent pore section index, the implementation of the saturation determining device based on the equivalent pore section index can be implemented by referring to the saturation determining method based on the equivalent pore section index, and the repetition is omitted. As used below, the term "unit" or "module" may be a combination of software and/or hardware that implements the intended function. While the system described in the following embodiments is preferably implemented in software, implementation in hardware, or a combination of software and hardware, is also possible and contemplated.
An embodiment of the present invention provides a specific implementation manner of a saturation determining device based on an equivalent pore section index, which can implement a saturation determining method based on an equivalent pore section index, and referring to fig. 9, the saturation determining device based on an equivalent pore section index specifically includes the following contents:
An equivalent pore section index determination unit 10 for obtaining an equivalent pore section index;
a pore structure index calculating unit 20 for calculating a pore structure index using the equivalent pore section index;
a water saturation determination unit 30 for determining reservoir water saturation using the porosity structural index.
In one embodiment, the equivalent pore section index obtaining unit includes:
the model creation module is used for creating an equivalent rock volume model;
and the equivalent pore section index acquisition module is used for acquiring the equivalent pore section index by utilizing an equivalent rock volume model according to the parallel resistance principle.
In one embodiment, the equivalent pore section index obtaining module includes:
logging curve data acquisition module: the method comprises the steps of acquiring log curve data;
a logging curve data calculation module: the method comprises the steps of obtaining formation water resistivity, reservoir total porosity and reservoir flushing zone resistivity according to logging curve data;
equivalent pore section index first calculation module: and the method is used for calculating the equivalent pore section index according to the product of the formation water resistivity and the total porosity of the reservoir and the resistivity of the reservoir flushing zone.
In one embodiment, the equivalent pore section index obtaining module includes:
Rock core experimental data and formation water analysis data acquisition module: the method comprises the steps of acquiring rock core experimental data, stratum water analysis data and logging curve data;
rock core experimental data and formation water analysis data calculation module: the method comprises the steps of obtaining formation water resistivity, total porosity of a reservoir and reservoir resistivity when the reservoir is saturated with formation water by 100% according to rock core experimental data, formation water analysis data and logging curve data;
the equivalent pore section index calculation module is a second calculation module: and the method is used for calculating the equivalent pore section index according to the product of the formation water resistivity and the total porosity of the reservoir and the reservoir resistivity when the reservoir is saturated with formation water by 100 percent.
In an embodiment, the water saturation determination unit is specifically configured to: and calculating the water saturation of the reservoir according to the product of the formation water resistivity, the lithology coefficient and the saturation coefficient, the reservoir undisturbed formation resistivity, the cementation index, the reservoir total porosity and the pore structure index.
From the above description, the invention provides a saturation determining device based on equivalent pore section indexes, firstly, on the basis of an equivalent rock volume model, equivalent pore section indexes are calculated by using two methods of a core experiment and a conventional logging curve, the equivalent pore section indexes can reflect a core pore structure, then, a pore structure index is defined by the equivalent pore section indexes, the pore structure index simultaneously considers the porosity of a reservoir and the pore structure of the reservoir, and finally, the saturation is calculated by using the pore structure index and an Archie formula. The saturation calculation precision and interpretation coincidence rate are greatly improved, so that the purpose of calculating the complex reservoir saturation with high precision of the conventional logging series is achieved, and the problem of quantitative evaluation of the reservoir saturation is effectively solved.
The embodiment of the present application further provides a specific implementation manner of an electronic device capable of implementing all the steps in the saturation determining method based on the equivalent pore section index in the foregoing embodiment, and referring to fig. 10, the electronic device specifically includes the following contents:
a processor 1201, a memory 1202, a communication interface (Communications Interface) 1203, and a bus 1204;
wherein the processor 1201, the memory 1202 and the communication interface 1203 perform communication with each other through the bus 1204; the communication interface 1203 is configured to implement information transmission between related devices such as a server device, a detection device, and a user device.
The processor 1201 is configured to invoke a computer program in the memory 1202, and when the processor executes the computer program, the processor implements all the steps in the saturation determining method based on equivalent pore section index in the above embodiment, for example, when the processor executes the computer program, the processor implements the following steps:
step 100: and obtaining the equivalent pore section index.
Step 200: and calculating the pore structure index by using the equivalent pore section index.
Step 300: reservoir water saturation is determined using the porosity structural index.
From the above description, it can be known that, in the electronic device in the embodiment of the present application, firstly, on the basis of an equivalent rock volume model, an equivalent pore section index is calculated by using two methods of a core experiment and a conventional logging curve, the equivalent pore section index can reflect a core pore structure, then, a pore structure index is defined by the equivalent pore section index, the pore structure index considers both the porosity of a reservoir and the pore structure of the reservoir, and finally, the saturation is calculated by using the pore structure index and an alchi formula. The saturation calculation precision and interpretation coincidence rate are greatly improved, so that the purpose of calculating the complex reservoir saturation with high precision of the conventional logging series is achieved, and the problem of quantitative evaluation of the reservoir saturation is effectively solved.
The embodiments of the present application also provide a computer-readable storage medium capable of implementing all the steps in the saturation determination method based on equivalent pore section index in the above embodiments, and a computer program stored on the computer-readable storage medium, which when executed by a processor implements all the steps in the saturation determination method based on equivalent pore section index in the above embodiments, for example, the following steps are implemented when the processor executes the computer program:
step 100: and obtaining the equivalent pore section index.
Step 200: and calculating the pore structure index by using the equivalent pore section index.
Step 300: reservoir water saturation is determined using the porosity structural index.
As can be seen from the above description, the computer readable storage medium in the embodiments of the present application firstly calculates an equivalent pore section index based on an equivalent rock volume model by using two methods of core experiments and conventional logging curves, the equivalent pore section index can reflect the core pore structure, then defines a pore structure index by the equivalent pore section index, the pore structure index simultaneously considers the reservoir porosity and the reservoir pore structure, and finally calculates the saturation by the pore structure index and the alchi formula. The saturation calculation precision and interpretation coincidence rate are greatly improved, so that the purpose of calculating the complex reservoir saturation with high precision of the conventional logging series is achieved, and the problem of quantitative evaluation of the reservoir saturation is effectively solved.
In this specification, each embodiment is described in a progressive manner, and identical and similar parts of each embodiment are all referred to each other, and each embodiment mainly describes differences from other embodiments. In particular, for a hardware+program class embodiment, the description is relatively simple, as it is substantially similar to the method embodiment, as relevant see the partial description of the method embodiment.
The foregoing describes specific embodiments of the present disclosure. Other embodiments are within the scope of the following claims. In some cases, the actions or steps recited in the claims can be performed in a different order than in the embodiments and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.
Although the present application provides method operational steps as described in the examples or flowcharts, more or fewer operational steps may be included based on conventional or non-inventive labor. The order of steps recited in the embodiments is merely one way of performing the order of steps and does not represent a unique order of execution. When implemented by an actual device or client product, the instructions may be executed sequentially or in parallel (e.g., in a parallel processor or multi-threaded processing environment) as shown in the embodiments or figures.
Although the present description provides method operational steps as described in the examples or flowcharts, more or fewer operational steps may be included based on conventional or non-inventive means. The order of steps recited in the embodiments is merely one way of performing the order of steps and does not represent a unique order of execution. When implemented in an actual device or end product, the instructions may be executed sequentially or in parallel (e.g., in a parallel processor or multi-threaded processing environment, or even in a distributed data processing environment) as illustrated by the embodiments or by the figures. The terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, it is not excluded that additional identical or equivalent elements may be present in a process, method, article, or apparatus that comprises a described element.
For convenience of description, the above devices are described as being functionally divided into various modules, respectively. Of course, when implementing the embodiments of the present disclosure, the functions of each module may be implemented in the same or multiple pieces of software and/or hardware, or a module that implements the same function may be implemented by multiple sub-modules or a combination of sub-units, or the like. The above-described apparatus embodiments are merely illustrative, for example, the division of the units is merely a logical function division, and there may be additional divisions when actually implemented, for example, multiple units or components may be combined or integrated into another system, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or units, which may be in electrical, mechanical or other form.
Those skilled in the art will also appreciate that, in addition to implementing the controller in a pure computer readable program code, it is well possible to implement the same functionality by logically programming the method steps such that the controller is in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers, etc. Such a controller can be regarded as a hardware component, and means for implementing various functions included therein can also be regarded as a structure within the hardware component. Or even means for achieving the various functions may be regarded as either software modules implementing the methods or structures within hardware components.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In one typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include volatile memory in a computer-readable medium, random Access Memory (RAM) and/or nonvolatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of computer-readable media.
Computer readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of storage media for a computer include, but are not limited to, phase change memory (PRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape disk storage or other magnetic storage devices, or any other non-transmission medium, which can be used to store information that can be accessed by a computing device. Computer-readable media, as defined herein, does not include transitory computer-readable media (transmission media), such as modulated data signals and carrier waves.
It will be appreciated by those skilled in the art that embodiments of the present description may be provided as a method, system, or computer program product. Accordingly, the present specification embodiments may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present description embodiments may take the form of a computer program product on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein.
The present embodiments may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The embodiments of the specification may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
In this specification, each embodiment is described in a progressive manner, and identical and similar parts of each embodiment are all referred to each other, and each embodiment mainly describes differences from other embodiments. In particular, for system embodiments, since they are substantially similar to method embodiments, the description is relatively simple, as relevant to see a section of the description of method embodiments. In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the embodiments of the present specification. In this specification, schematic representations of the above terms are not necessarily directed to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, the different embodiments or examples described in this specification and the features of the different embodiments or examples may be combined and combined by those skilled in the art without contradiction.
The foregoing is merely an example of an embodiment of the present disclosure and is not intended to limit the embodiment of the present disclosure. Various modifications and variations of the illustrative embodiments will be apparent to those skilled in the art. Any modification, equivalent replacement, improvement, or the like, which is within the spirit and principles of the embodiments of the present specification, should be included in the scope of the claims of the embodiments of the present specification.
Claims (9)
1. A saturation determination method based on equivalent pore cross-section index, comprising:
obtaining an equivalent pore section index;
calculating a pore structure index using the equivalent pore section index;
determining reservoir water saturation using the pore structure index;
obtaining the equivalent pore cross-section index includes:
creating an equivalent rock volume model;
according to the parallel resistance principle, the equivalent pore section index is obtained by utilizing the equivalent rock volume model;
the obtaining the equivalent pore section index by using the equivalent rock volume model according to the resistance parallel principle comprises the following steps:
acquiring logging curve data;
acquiring formation water resistivity, reservoir total porosity and reservoir flushing zone resistivity according to the logging curve data;
Calculating the equivalent pore section index according to the product of the formation water resistivity and the reservoir total porosity and the reservoir flushing zone resistivity;
the calculation formula for calculating the equivalent pore section index according to the product of the formation water resistivity and the reservoir total porosity and the reservoir flushing zone resistivity is as follows:
wherein, IPS is the equivalent pore section index, dimensionless; r is R w Is the formation water resistivity, Ω·m; phi is the total porosity of the reservoir,%; r is R xo Flushing the reservoir with resistivity, Ω·m;
the calculating the pore structure index using the equivalent pore cross section index comprises:
m＝-logIPS
wherein m is a pore structure index, and is dimensionless.
2. The saturation determination method of claim 1, wherein the obtaining the equivalent pore section index using the equivalent rock volume model according to a resistive parallel principle comprises:
acquiring rock core experimental data, formation water analysis data and logging curve data;
acquiring formation water resistivity, total porosity of the reservoir and reservoir resistivity when the reservoir is saturated with formation water by 100% according to the rock core experimental data, the formation water analysis data and the logging curve data;
calculating the equivalent pore section index according to the product of the formation water resistivity and the total porosity of the reservoir and the reservoir resistivity when the reservoir is saturated with formation water by 100 percent;
The calculation formula for calculating the equivalent pore section index according to the product of the formation water resistivity and the total porosity of the reservoir and the reservoir resistivity when the reservoir is saturated with formation water by 100 percent is as follows:
wherein, IPS is the equivalent pore section index, dimensionless; r is R w Is the formation water resistivity, Ω·m; phi is the total porosity of the reservoir,%; r is R 0 Is the reservoir resistivity at 100% saturation of the reservoir with formation water, Ω·m.
3. The saturation determination method of claim 1, wherein the pore structure index m value ranges between 0-3.
4. The saturation determination method of claim 1, wherein the determining reservoir water saturation using the pore structure index comprises:
and calculating the water saturation of the reservoir according to the product of the formation water resistivity, the lithology coefficient and the saturation coefficient, the reservoir undisturbed formation resistivity, the cementation index, the reservoir total porosity and the pore structure index.
5. The saturation determination method of claim 4, wherein the calculation formula for calculating the reservoir water saturation from the formation water resistivity, the product of the lithology coefficient and the saturation coefficient, the reservoir undisturbed formation resistivity, the cementing index, the reservoir total porosity and the pore structure index is:
Wherein, a: lithology coefficient, dimensionless; b: saturation coefficient, dimensionless; n: cementing index, dimensionless; sw: reservoir water saturation,%; r is R t : the resistivity of the undisturbed stratum of the reservoir is omega-m.
6. A saturation determination device based on equivalent pore cross-section index, comprising:
the equivalent pore section index obtaining unit is used for obtaining the equivalent pore section index;
a pore structure index calculation unit for calculating a pore structure index using the equivalent pore section index;
a water saturation determination unit for determining reservoir water saturation using the pore structure index;
the equivalent pore section index obtaining unit includes:
the model creation module is used for creating an equivalent rock volume model;
the equivalent pore section index acquisition module is used for acquiring the equivalent pore section index by utilizing the equivalent rock volume model according to the resistance parallel connection principle;
the equivalent pore section index obtaining module comprises:
logging curve data acquisition module: the method comprises the steps of acquiring log curve data;
a logging curve data calculation module: the method comprises the steps of obtaining formation water resistivity, reservoir total porosity and reservoir flushing zone resistivity according to logging curve data;
Equivalent pore section index first calculation module: the equivalent pore section index is calculated according to the product of the formation water resistivity and the total porosity of the reservoir and the resistivity of a reservoir flushing zone;
the equivalent pore section index obtaining module comprises:
rock core experimental data and formation water analysis data acquisition module: the method comprises the steps of acquiring rock core experimental data, stratum water analysis data and logging curve data;
rock core experimental data and formation water analysis data calculation module: the method comprises the steps of obtaining formation water resistivity, total porosity of a reservoir and reservoir resistivity when the reservoir is 100% saturated with formation water according to rock core experimental data, formation water analysis data and logging curve data;
the equivalent pore section index calculation module is a second calculation module: the equivalent pore section index is calculated according to the product of the formation water resistivity and the total porosity of the reservoir and the reservoir resistivity when the reservoir is 100% saturated with formation water;
the calculation formula for calculating the equivalent pore section index according to the product of the formation water resistivity and the reservoir total porosity and the reservoir flushing zone resistivity is as follows:
wherein, IPS is the equivalent pore section index, dimensionless; r is R w Is the formation water resistivity, Ω·m; phi is the total porosity of the reservoir,%; r is R xo Flushing the reservoir with resistivity, Ω·m;
the calculating the pore structure index using the equivalent pore cross section index comprises:
m＝-logIPS
wherein m is a pore structure index, and is dimensionless.
7. The saturation determination apparatus of claim 6, wherein the water saturation determination unit is specifically configured to: and calculating the water saturation of the reservoir according to the product of the formation water resistivity, the lithology coefficient and the saturation coefficient, the reservoir undisturbed formation resistivity, the cementation index, the reservoir total porosity and the pore structure index.
8. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the steps of the saturation determination method according to any one of claims 1 to 5 when the program is executed by the processor.
9. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of the saturation determination method of any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910112284.6A CN111562629B (en) | 2019-02-13 | 2019-02-13 | Saturation determination method and device based on equivalent pore section index |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910112284.6A CN111562629B (en) | 2019-02-13 | 2019-02-13 | Saturation determination method and device based on equivalent pore section index |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111562629A CN111562629A (en) | 2020-08-21 |
| CN111562629B true CN111562629B (en) | 2023-04-25 |
Family
ID=72069485
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910112284.6A Active CN111562629B (en) | 2019-02-13 | 2019-02-13 | Saturation determination method and device based on equivalent pore section index |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111562629B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114458305A (en) * | 2020-11-03 | 2022-05-10 | 中国石油天然气集团有限公司 | Method and device for determining pore structure coefficient |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6470274B1 (en) * | 1998-12-30 | 2002-10-22 | Baker Hughes Incorporated | Water saturation and sand fraction determination from borehole resistivity imaging tool, transverse induction logging and a tensorial dual water saturation model |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5349528A (en) * | 1990-03-12 | 1994-09-20 | Halliburton Logging Services, Inc. | Method apparatus for determination of porosity lithological composition |
| CN103422853B (en) * | 2012-05-24 | 2016-10-26 | 中国石油天然气集团公司 | A kind of method and device determining stratum water saturation |
| CN105114064B (en) * | 2015-08-04 | 2018-03-13 | 中国石油天然气股份有限公司 | Method for Determining Saturation of Tight Sandstone Reservoir |
| US10393641B2 (en) * | 2015-10-30 | 2019-08-27 | Schlumberger Technology Corporation | Methods of determining cementation exponent and saturation exponent in porous media from dielectric dispersion data |
-
2019
- 2019-02-13 CN CN201910112284.6A patent/CN111562629B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6470274B1 (en) * | 1998-12-30 | 2002-10-22 | Baker Hughes Incorporated | Water saturation and sand fraction determination from borehole resistivity imaging tool, transverse induction logging and a tensorial dual water saturation model |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111562629A (en) | 2020-08-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7933750B2 (en) | Method for defining regions in reservoir simulation | |
| Cacace et al. | Modelling of fractured carbonate reservoirs: outline of a novel technique via a case study from the Molasse Basin, southern Bavaria, Germany | |
| CA2692425C (en) | Method, program and computer system for scaling hydrocarbon reservoir model data | |
| Hosseini et al. | Geostatistical modeling and spatial distribution analysis of porosity and permeability in the Shurijeh-B reservoir of Khangiran gas field in Iran | |
| Al-Fatlawi et al. | A new practical method for predicting equivalent drainage area of well in tight gas reservoirs | |
| Arekhov et al. | The role of diffusion on reservoir performance in underground hydrogen storage | |
| CN111562629B (en) | Saturation determination method and device based on equivalent pore section index | |
| Abdulredah et al. | Building 3D geological model using non-uniform gridding for Mishrif reservoir in Garraf oilfield | |
| Baker et al. | Reservoir characterizations and reservoir performance of Mishrif Formation in Amara Oil Field | |
| Zhang et al. | Using digital rock modeling to estimate permeability and capillary pressure from NMR and geochemical logs | |
| Ogunyomi et al. | History matching and rate forecasting in unconventional oil reservoirs using an approximate analytical solution to the double porosity model | |
| Al-Thuwaini et al. | Innovative approach to assist history matching using artificial intelligence | |
| Watson et al. | Flow diagnostics for model ensembles | |
| Rastegarnia et al. | Estimation of flow zone indicator distribution by using seismic data: a case study from a central Iranian oilfield | |
| Mohamed et al. | A new approach to improve reservoir modeling via machine learning | |
| Junker et al. | Modern Approach to Estimation of Uncertainty of Predictions With Dynamic Reservoir Simulation—A Case Study of a German Rotliegend Gas Field | |
| Worthington | Maximizing the effectiveness of integrated reservoir studies: practical approaches to improving the process and results | |
| Sahni et al. | History match case study: Use of assisted history match tools on single-well models in conjunction with a full-field history match | |
| Lun et al. | A procedure for integrating geologic concepts into history matching | |
| Irving et al. | An uncertainty modelling workflow for structurally compartmentalized reservoirs | |
| Dong et al. | An introduction—Carbonate geophysics | |
| Reedy et al. | Analysis of finely laminated deep marine turbidites: Integration of core and log data yields a novel interpretation model | |
| Damsleth et al. | Scale consistency from cores to geologic description | |
| Fassihi et al. | Aquifer Influx Versus Water Injection in GoM | |
| Berentsen et al. | An efficient method for explicit hydraulic fracture representation in full field reservoir simulation |
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 |