CN115614029A - Method for identifying dolomite reservoir and computer-readable storage medium - Google Patents
Method for identifying dolomite reservoir and computer-readable storage medium Download PDFInfo
- Publication number
- CN115614029A CN115614029A CN202110795597.3A CN202110795597A CN115614029A CN 115614029 A CN115614029 A CN 115614029A CN 202110795597 A CN202110795597 A CN 202110795597A CN 115614029 A CN115614029 A CN 115614029A
- Authority
- CN
- China
- Prior art keywords
- curve
- reservoir
- indication
- indication curve
- natural gamma
- 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.)
- Pending
Links
- 229910000514 dolomite Inorganic materials 0.000 title claims abstract description 65
- 239000010459 dolomite Substances 0.000 title claims abstract description 65
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000011161 development Methods 0.000 claims description 43
- 229930195733 hydrocarbon Natural products 0.000 claims description 14
- 150000002430 hydrocarbons Chemical class 0.000 claims description 14
- 238000005259 measurement Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 abstract description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 7
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 230000035699 permeability Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 244000309464 bull Species 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000009395 breeding Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a method for identifying a dolomite reservoir and a computer readable storage medium, wherein the method for identifying the dolomite reservoir comprises the following steps: 1) Collecting compensated neutron, natural gamma and deep lateral resistivity data; 2) Reconstructing the compensated neutrons, the natural gamma rays and the deep lateral resistivity to obtain an indication curve; 3) Intercepting the indication curve according to the threshold value of the indication curve of the area; 4) And the position corresponding to the intercepted indication curve is the dolomite reservoir. A computer readable storage medium storing a program implementing the method for identifying a dolomite reservoir, implementing the steps of data taking and processing in the method for identifying a dolomite reservoir.
Description
Technical Field
The invention relates to a method for identifying a dolomite reservoir, and belongs to the technical field of geological exploration of a dolomite reservoir hydrocarbon reservoir.
Background
The dolomite is a sedimentary carbonate rock, and carbonate rock reservoirs are wide in distribution range and rich in reserves in China, and have great potential in the aspect of increasing storage and yield. In the exploration and development process of the Marfan five-section multi-breeding dolomite reservoir, an important work in the exploration of reservoir identification oil and gas reservoirs is effectively and accurately carried out. When identifying a dolomite reservoir, the conventional method mostly utilizes natural gamma, acoustic time difference, gas measurement of total hydrocarbon, compensation of neutrons, deep lateral resistivity of 5 logging curves and porosity, permeability and gas saturation obtained by calculation according to the natural gamma, the acoustic time difference and the deep lateral resistivity to judge the development position of the reservoir. However, this method requires comprehensive utilization of multiple logging curve data and establishment of calculation formulas for porosity, permeability and gas saturation, and requires too many logging curves and a complicated calculation process.
The patent document with publication number CN111751905A discloses a carbonate reservoir effectiveness identification method based on a rock three-porosity model. The method firstly utilizes the time difference of the compensated density, the compensated neutron and the sound wave to calculate the porosity phi N Porosity phi of the porous material D 、φ b In the use of phi N And phi D Calculating the total porosity phi T . Then according to the total porosity phi T And acoustic porosity phi b Obtaining three porosities of the carbonate rock, namely the porosity phi of the communicated fracture-cavity respectively 2 Porosity phi of non-connected slot nc And the porosity of the matrix phi m . According to the total porosity phi T And porosity phi of connected slot 2 Porosity phi of non-connected slot n c and the porosity of the matrix phi m Obtaining a formation cementation index m; according to the porosity phi of the connected seam hole 2 And porosity of non-connected slot nc Obtaining the sum of the secondary porosity and the total porosity, quantitatively judging the effectiveness of the reservoir through a standard diagram of a ratio BZ of the secondary porosity to the total porosity and a cementation index m, and judging the effectiveness of the reservoir when the ratio BZ of the secondary porosity to the total porosity is more than 0.02 and the cementation index m is between 1 and 2.18Belonging to the stratum of origin. Although the reservoir stratum can be determined, the method needs a large amount of calculation processes and is complex to implement.
Disclosure of Invention
The invention aims to provide a method for identifying a dolomite reservoir and a computer readable storage medium, which are used for solving the problems of excessive requirements on logging curve data and complex calculation process in the process of identifying the dolomite reservoir.
In order to achieve the above object, the scheme of the invention comprises:
the invention discloses a method for identifying a dolomite reservoir, which comprises the following steps:
1) Collecting compensated neutron, natural gamma and deep lateral resistivity data;
2) Substituting the compensated neutrons, the natural gamma rays and the deep lateral resistivity into a reconstruction formula to obtain an indication curve;
3) Intercepting the indication curve according to the threshold value of the indication curve of the area;
4) The depth corresponding to the intercepted indication curve is the dolomite reservoir;
the method for obtaining the threshold value of the indication curve comprises the following steps:
a) Collecting logging data of the region, wherein the logging data comprise: deep lateral resistivity, acoustic time difference, compensated neutrons, gas-measured total hydrocarbons, natural gamma values;
b) Substituting the compensated neutron curve, the natural gamma curve and the deep lateral resistivity data into a reconstruction formula to obtain an indication curve;
c) Determining the development position of a dolomite reservoir according to natural gamma, deep lateral resistivity, compensation density, gas measurement of total hydrocarbons and a sonic time difference curve;
d) And comparing the development position of the dolomite reservoir with the indication curve, and reserving the part of the indication curve corresponding to the development position of the reservoir, wherein the minimum value of the part of the indication curve is the threshold value of the indication curve.
The beneficial effects of doing so are: and superposing the indication curve and a parameter curve of the conventional judgment reservoir development position to determine a threshold value, and determining the development position of the whole five-section reservoir better through the threshold value. The method solves the problems of multiple logging curves and complex calculation process required by identification of the Mawu-section dolomite reservoir, simply and quickly judges the development position of the Mawu-section dolomite reservoir, and provides effective support for exploration and development of the dolomite reservoir hydrocarbon reservoir.
Further, in step 2), the reconstruction formula is:
wherein A is an indication curve, CNL is a compensation neutron parameter, GR is a natural gamma curve parameter, and LLD is a deep lateral resistivity parameter.
Further, in the step d), the method for comparing the development position of the dolomite reservoir with the indication curve is superposition display.
The beneficial effects of doing so are: the indication curve threshold is obtained more straightly and clearly by superposition display.
Further, a computer readable storage medium for identifying a dolomite reservoir, wherein the storage medium has stored therein instructions for implementing the method of:
1) Collecting compensated neutron, natural gamma and deep lateral resistivity data;
2) Substituting the compensated neutrons, the natural gamma rays and the deep lateral resistivity into a reconstruction formula to obtain an indication curve;
3) Intercepting the indication curve according to the threshold value of the indication curve of the area;
4) The depth corresponding to the intercepted indication curve is the dolomite reservoir;
the method for obtaining the threshold value of the indication curve comprises the following steps:
a) Collecting logging data of an area, wherein the logging data comprise: deep lateral resistivity, acoustic time difference, compensated neutrons, gas measured total hydrocarbons, natural gamma values;
b) Substituting the compensated neutron curve, the natural gamma curve and the deep lateral resistivity data into a reconstruction formula to obtain an indication curve;
c) Determining the development position of a dolomite reservoir according to natural gamma, deep lateral resistivity, compensation density, gas measurement of all hydrocarbons and acoustic time difference curves;
d) And comparing the development position of the dolomite reservoir with the indication curve, and reserving the part of the indication curve corresponding to the development position of the reservoir, wherein the minimum value of the part of the indication curve is the threshold value of the indication curve.
Further, in step 2), the reconstruction formula is:
wherein A is an indication curve, CNL is a compensation neutron parameter, GR is a natural gamma curve parameter, and LLD is a deep lateral resistivity parameter.
Further, in the step d), the method for comparing the development position of the dolomite reservoir with the indication curve is superposition display.
Drawings
FIG. 1 is a flow chart of a method of identifying a dolomitic reservoir of the present invention;
FIG. 2 is a graph of the effect of identifying development sites of a dolomite reservoir according to the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments.
The method comprises the following steps:
as shown in fig. 1, a method for identifying a dolomite reservoir according to the present invention is implemented by the following steps:
1) Collecting compensation neutron data, natural gamma data and deep lateral resistivity data;
2) Substituting the collected compensation neutron data, natural gamma data and deep lateral resistivity data into a reconstruction formula to obtain an indication curve;
3) Intercepting an indication curve according to the threshold value of the indication curve of the area where the dolomite is located;
4) The depth corresponding to the intercepted indication curve is the dolomite reservoir;
the method for obtaining the threshold value of the indication curve comprises the following steps:
a) The method comprises the following steps of collecting logging data of a region where the dolomite is located, wherein the logging data comprise: deep lateral resistivity data, sound wave time difference data, compensation neutron data, gas measurement all hydrocarbon data and natural gamma value data;
b) Substituting the collected compensation neutron data, natural gamma data and deep lateral resistivity data into a reconstruction formula to obtain an indication curve;
c) According to a method in the prior art (such as the method described in CN 111751905A), determining a development position of a dolomite reservoir based on the collected natural gamma data, deep lateral resistivity data, compensated density data, gas measurement all hydrocarbon data and acoustic wave time difference data;
d) And superposing and displaying the development position of the dolomite reservoir obtained according to the prior art and the indication curve, and determining the threshold value of the indication curve.
The conception of the invention is as follows: according to logging data acquired in the prior art, natural gamma data in the logging data can reflect lithology under a certain condition, compensation subdata can reflect a reservoir development condition under a certain condition, and deep lateral resistivity data can reflect a reservoir gas condition under a certain condition. Therefore, the method acquires three groups of data including natural gamma, compensated neutrons and deep lateral resistivity, reconstructs the three groups of data to form an indication curve, compares the dolomite reservoir position curve obtained by the prior art method with the indication curve, intercepts the indication curve after comparison, reserves the part of the indication curve corresponding to the reservoir position, and takes the value corresponding to the intercepted position as the threshold value.
For the identification of the positions of the dolomite reservoir at other positions in the area, only the compensation neutron data, the natural gamma data and the deep lateral resistivity data are collected and reconstructed to obtain an indication curve, the indication curve is intercepted according to the obtained threshold value of the position of the area, and the position (depth) corresponding to the remaining part after interception is judged as the position of the dolomite reservoir.
According to the method, after the threshold value is obtained for one area, the dolomite reservoir can be identified only by acquiring the compensated neutron data, the natural gamma data and the deep lateral resistivity data for other positions and performing new reconstruction, so that the acquisition amount and the calculation amount of logging data are greatly reduced, and the identification complexity of the dolomite reservoir is reduced.
Taking five-section reservoir of Ordos Danigen gas field horse as an example, firstly collecting five-section logging data of Ordos Danigen field horse: natural gamma data, deep lateral resistivity data, gas logging total hydrocarbon data, compensation neutron data, and acoustic time difference data.
In order to better reflect the development position of the reservoir, the acquired Compensation Neutrons (CNL), natural gamma curve (GR) and deep lateral resistivity (LLD) data are reconstructed and processed:
the indicating curve a is obtained.
And superposing and displaying the reservoir development position obtained by the prior art and the indication curve to determine the threshold value of the indication curve of the developing dolomite reservoir.
The method for judging the development position of the reservoir in the prior art comprises the following steps: and 4 logging curves of a gas-logging all-hydrocarbon logging curve (QT), a natural gamma curve (GR), an acoustic time difference (AC) and deep lateral resistivity (LLD) are utilized, and meanwhile, 3 curves of the natural gamma curve (GR), the acoustic time difference (AC) and the deep lateral resistivity (LLD) are utilized to calculate porosity (phi), permeability (K), shale content (SH) and gas saturation (Sg), so that the development position of the reservoir stratum is judged.
The lower limit of the parameters for development of the horse five-section reservoir is generally: AC is more than or equal to 157 MuS/M, QT is more than or equal to 2 percent, LLD is more than or equal to 45 ohm M and less than or equal to 1300 ohm M, phi is more than or equal to 2.5 percent, K is more than or equal to 0.03mD, SH is less than or equal to 13.0 percent, and Sg is more than or equal to 50.0 percent.
The method for determining the threshold value of the indicating curve comprises the following steps: and (3) displaying the development position of the reservoir layer and the indicating curve in an overlapping mode, and reserving the overlapping part of the indicating curve and the development position of the reservoir layer, wherein the overlapping part is not a continuous curve due to discontinuous reservoir layers, the overlapping position of each residual section of curve corresponds to a fixed numerical value in theory, and the fixed numerical value is an indicating curve threshold value.
If the actual value fluctuates within a small range, those skilled in the art should understand that the threshold curve value is calculated by a statistical method, such as a data processing method of averaging the data in the fluctuation region.
And (5) adjusting (intercepting) the indication curve according to a threshold value, and only reserving a part which is overlapped with the development position of the reservoir, wherein the part is the development position of the whole reservoir of the Mawu section of dolomite.
As shown in fig. 2, the indication curve of five sections of dolomite in the mare in the bulls-land gas field is superimposed with the parameter curve of the reservoir development position to obtain an indication curve with a threshold value of 0.8, and only the part of the indication curve greater than 0.8 is reserved by adjusting the indication curve according to the threshold value, and the depth corresponding to the part is the reservoir development position. When the development positions of the five sections of dolomite in the bulls land gas field are identified later, the conventional method is not needed to be used again for calculating and determining the development positions of the dolomite, and only the indication curve is needed to be adjusted according to 0.8, namely the development positions of the reservoirs of the five sections of dolomite in the bulls land gas field are obtained. If the development positions of the dolomite need to be identified in other areas of the place, only the threshold value of the indication curve of the place needs to be obtained according to the method in the first identification, the development positions of the reservoirs of the whole section of the dolomite need to be determined according to the threshold value of the indication curve, and the development positions of the reservoirs of the whole section of the dolomite need to be determined only by combining the first obtained threshold value of the indication curve and the indication curve in the later identification.
Computer-readable storage medium embodiments:
a computer readable storage medium storing a program implementing the method for identifying a dolomite reservoir, implementing the steps of data taking and processing in the method for identifying a dolomite reservoir.
The media referred to in this embodiment is programmable data processing apparatus having stored thereon computer program instructions. The media described above may also be one or more computer-readable storage media having computer-usable program code embodied therewith.
Claims (6)
1. A method of identifying a dolomite reservoir, the method of identifying a dolomite reservoir comprising:
1) Collecting compensated neutron, natural gamma and deep lateral resistivity data;
2) Substituting the compensated neutrons, the natural gamma rays and the deep lateral resistivity into a reconstruction formula to obtain an indication curve;
3) Intercepting the indication curve according to the threshold value of the indication curve of the area;
4) The position corresponding to the intercepted indication curve is the dolomite reservoir;
the method for obtaining the threshold value of the indication curve comprises the following steps:
a) Collecting logging data of an area, wherein the logging data comprise: deep lateral resistivity, acoustic time difference, compensated neutrons, gas-measured total hydrocarbons, natural gamma values;
b) Substituting the compensated neutron curve, the natural gamma curve and the deep lateral resistivity data into a reconstruction formula to obtain an indication curve;
c) Determining the development position of a dolomite reservoir according to natural gamma, deep lateral resistivity, compensation density, gas measurement of all hydrocarbons and acoustic time difference curves;
d) And comparing the development position of the dolomite reservoir with the indication curve, and reserving a part corresponding to the development position of the reservoir on the indication curve, wherein the minimum value of the part is the threshold value of the indication curve.
2. A method for identifying a dolomite reservoir as claimed in claim 1, wherein in step 2), the reconstruction formula is:
wherein A is an indication curve, CNL is a compensation neutron parameter, GR is a natural gamma curve parameter, and LLD is a deep lateral resistivity parameter.
3. A method of identifying a dolomite reservoir as claimed in claim 1 wherein in step d) the method of comparing the development location of the dolomite reservoir to the indicative curve is a superimposed display.
4. A computer readable storage medium for identifying a dolomite reservoir, the storage medium having stored therein instructions for performing a method comprising:
1) Collecting compensated neutron, natural gamma and deep lateral resistivity data;
2) Substituting a reconstruction formula into the compensation neutrons, the natural gamma rays and the deep lateral resistivity to obtain an indication curve;
3) Intercepting the indication curve according to the threshold value of the indication curve of the area;
4) The depth corresponding to the intercepted indication curve is the dolomite reservoir;
the method for obtaining the threshold value of the indication curve comprises the following steps:
a) Collecting logging data of an area, wherein the logging data comprise: deep lateral resistivity, acoustic time difference, compensated neutrons, gas-measured total hydrocarbons, natural gamma values;
b) Substituting the compensated neutron curve, the natural gamma curve and the deep lateral resistivity data into a reconstruction formula to obtain an indication curve;
c) Determining the development position of a dolomite reservoir according to natural gamma, deep lateral resistivity, compensation density, gas measurement of total hydrocarbons and a sonic time difference curve;
d) And comparing the development position of the dolomite reservoir with the indication curve, and reserving the part of the indication curve corresponding to the development position of the reservoir, wherein the minimum value of the part is the threshold value of the indication curve.
5. The computer-readable storage medium for identifying a dolomite reservoir as set forth in claim 4, wherein in the step 2), the reconstruction formula is:
wherein A is an indication curve, CNL is a compensation neutron parameter, GR is a natural gamma curve parameter, and LLD is a deep lateral resistivity parameter.
6. The computer readable storage medium for identifying a dolomite reservoir as set forth in claim 4, wherein the method for comparing the development position of the dolomite reservoir with the indication curve in step d) is a folding display.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110795597.3A CN115614029A (en) | 2021-07-14 | 2021-07-14 | Method for identifying dolomite reservoir and computer-readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110795597.3A CN115614029A (en) | 2021-07-14 | 2021-07-14 | Method for identifying dolomite reservoir and computer-readable storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115614029A true CN115614029A (en) | 2023-01-17 |
Family
ID=84856007
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110795597.3A Pending CN115614029A (en) | 2021-07-14 | 2021-07-14 | Method for identifying dolomite reservoir and computer-readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115614029A (en) |
-
2021
- 2021-07-14 CN CN202110795597.3A patent/CN115614029A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103670388B (en) | Method for evaluating organic carbon content of shale | |
| CN101930082B (en) | Method for discriminating reservoir fluid type by using resistivity data | |
| CN104278991B (en) | Saline Lake Facies hydrocarbon source rock organic carbon and the polynary well logging computational methods of hydrocarbon potential | |
| CN107449707A (en) | Three-dimensional characterization determination method and device for quantification of pores with different scales in shale reservoir | |
| CN103883322B (en) | Shale gas reservoirs exploration method and device | |
| CN107795320B (en) | Calculation method for horizontal well carbonate reservoir parameters | |
| CN110955982B (en) | Metamorphic rock reservoir permeability calculation method and device and computer storage medium | |
| CN103698811A (en) | Carbonate rock structure component well logging quantitative identification method and application thereof | |
| CA2867583C (en) | Fracking method for fracking intervals of a horizontal drilling zone in a sweet spot range based on measurements of resistivity and neutron logging data in the horizontal drillingzone | |
| CN105447762A (en) | Calculation method for low permeability reservoir flooding information of fluid replacement | |
| CN112177605B (en) | Method for determining main control factors of favorable gathering area of sea natural gas hydrate | |
| US20230160304A1 (en) | Method and system for predicting relative permeability curve based on machine learning | |
| CN111008482B (en) | Metamorphic rock reservoir water saturation calculation method and device | |
| CN104834007B (en) | Method for calculating filling degree of carbonate fracture-cave reservoir in seismic inversion process | |
| CN112861890A (en) | Reservoir evaluation model construction method and reservoir identification method | |
| CN112185469B (en) | Method for predicting favorable gathering area of sea natural gas hydrate | |
| CN106156456A (en) | Sand-conglomerate body invalid reservoir subtraction method | |
| CN113775326A (en) | Method and device for evaluating movable water saturation, electronic device, and medium | |
| CN115614029A (en) | Method for identifying dolomite reservoir and computer-readable storage medium | |
| CN116930023A (en) | Fine interpretation method and device for dense sandstone phase-control classified porosity logging | |
| Roberts-Ashby et al. | An evaluation of porosity and potential use for carbon dioxide storage in the Upper Cretaceous Lawson Formation and Paleocene Cedar Keys Formation of south-central and southern Florida | |
| CN103353618A (en) | Calculating method of gas resource amount | |
| Shen et al. | Hierarchical approach to modeling karst and fractures in carbonate karst reservoirs in the Tarim Basin | |
| CN105528732B (en) | Method for predicting productivity of gas testing well | |
| CN110633904A (en) | Metamorphic rock reservoir classification method and device |
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 |