CN112012726B - Lithology recognition method - Google Patents
Lithology recognition method Download PDFInfo
- Publication number
- CN112012726B CN112012726B CN201910463463.4A CN201910463463A CN112012726B CN 112012726 B CN112012726 B CN 112012726B CN 201910463463 A CN201910463463 A CN 201910463463A CN 112012726 B CN112012726 B CN 112012726B
- Authority
- CN
- China
- Prior art keywords
- value
- curve
- limestone
- representing
- range
- 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
- 238000000034 method Methods 0.000 title claims abstract description 62
- 235000019738 Limestone Nutrition 0.000 claims abstract description 69
- 239000006028 limestone Substances 0.000 claims abstract description 69
- 238000013507 mapping Methods 0.000 claims abstract description 46
- 238000003384 imaging method Methods 0.000 claims abstract description 4
- 101100409194 Rattus norvegicus Ppargc1b gene Proteins 0.000 claims description 44
- 208000019888 Circadian rhythm sleep disease Diseases 0.000 claims description 13
- 208000001456 Jet Lag Syndrome Diseases 0.000 claims description 13
- 208000033915 jet lag type circadian rhythm sleep disease Diseases 0.000 claims description 13
- 238000012937 correction Methods 0.000 claims description 12
- 239000012530 fluid Substances 0.000 claims description 6
- 239000011148 porous material Substances 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 3
- 230000005477 standard model Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 238000012545 processing Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 6
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
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
-
- 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
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Geochemistry & Mineralogy (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
A lithology recognition method, comprising: step one, respectively determining corresponding imaging numerical ranges according to the obtained different logging parameter values of the limestone skeleton of the work area to be analyzed; step two, respectively carrying out mapping operation on each logging parameter curve of the work area to be analyzed according to the mapping data value range; and step three, judging whether an overlapping area exists in the curve obtained through the mapping operation, wherein if the overlapping area exists, the lithology of the overlapping area is judged to be limestone. The method can simply, conveniently, rapidly and accurately realize the identification of the lithology of the limestone, and is suitable for field application. Compared with the existing method, the method has the advantages that the data processing operation is simpler and more convenient, the limestone lithology identification efficiency can be effectively improved, the limestone lithology can be rapidly identified, and therefore an effective auxiliary means is provided for actual production explanation.
Description
Technical Field
The invention relates to the technical field of oil and gas exploration and development, in particular to a lithology recognition method.
Background
For the identification of various lithologies, the most direct and reliable methods are based on sidewall coring and imaging logging, but the implementation cost of the methods is too high and the time consumption is long, and the methods cannot be applied to each well. So how to quickly identify lithology using conventional log curves is a common challenge in the industry.
Common methods for rapid identification of lithology using conventional log curves include a plate method, a two-dimensional and three-dimensional intersection map method, and some automatic identification algorithms (e.g., an artificial neural network method, a fuzzy clustering method, etc.). However, the automatic recognition algorithm requires a long-time data preprocessing process, is relatively complex to operate, is prone to over-fitting, and has a problem of large errors in certain work areas. The parameters often considered by the plate method and the two-dimensional intersection image method are too few, and the three-dimensional intersection image method is not favorable for rapid judgment of explanatory staff visually.
Therefore, a method capable of rapidly and accurately identifying lithology is needed.
Disclosure of Invention
In order to solve the problems, the invention provides a lithology recognition method, which comprises the following steps:
step one, respectively determining corresponding imaging numerical ranges according to the obtained different logging parameter values of the limestone skeleton of the work area to be analyzed;
step two, mapping operation is carried out on each logging parameter curve of the work area to be analyzed according to the mapping data value range;
and step three, judging whether an overlapping area exists in the curve obtained through the mapping operation, wherein if the overlapping area exists, the lithology of the overlapping area is judged to be limestone.
According to one embodiment of the invention, in the first step, different logging parameter values of the limestone framework are determined by using a lithology standard model in combination with the pure limestone core test measurement results of the work area to be analyzed.
According to one embodiment of the invention, the logging parameters include at least two of the following:
acoustic time difference, compensating neutrons and density.
According to one embodiment of the invention, in said step two,
based on a preset initial compensation neutron mapping range, respectively determining the ratio of the compensation neutron numerical value of the limestone framework to the whole mapping range from the upper mapping boundary distance to the lower mapping boundary distance in the mapping numerical value range according to the compensation neutron numerical value, and determining the compensation neutron mapping numerical value range according to the ratio.
According to one embodiment of the invention, the ratio of the compensated neutron value of the limestone skeleton to the overall range from the upper boundary distance of the plot under the range of the plotted values is determined according to the following expression:
wherein perc high Representing the ratio of the compensated neutron value of the limestone skeleton to the integral range of the upper boundary distance of the diagramming under the diagramming value range, CNL high And CNL low Respectively representing the upper limit value and the lower limit value of the mapping range of the preset initial compensation neutrons, CNL lim And the compensated neutron value of the limestone skeleton of the work area to be analyzed is represented.
According to one embodiment of the invention, the ratio of the compensated neutron value of the limestone skeleton to the overall range from the mapped lower boundary distance under the mapped value range is determined according to the following expression:
perc low ＝1-perc high
wherein perc low Representing the ratio of the compensated neutron value of the limestone skeleton to the overall range of the mapped lower boundary distance under the mapped numerical range, perc high The compensated neutron value of the limestone skeleton is represented by the ratio of the distance from the upper boundary of the diagramming to the whole range in the diagramming value range.
According to one embodiment of the invention, in said step two,
determining a lower boundary of a density curve graph according to the upper limit of the porosity of the work area to be analyzed and the density value of the limestone skeleton;
and determining the upper boundary of the density curve graph according to the lower boundary of the density curve graph and the numerical range of the compensation neutron curve graph.
According to one embodiment of the invention, the lower graphically boundary of the density curve is determined according to the following expression:
DEN low ＝α[DEN lim -φ high (DEN lim -DEN f )]
wherein DEN low Representing the lower boundary of the density curve, alpha representing the first correction factor, DEN lim Representing the density value of the limestone skeleton, phi high Represents the upper limit of porosity, DEN f Representing the density value of the pore fluid.
According to one embodiment of the invention, the upper graphically boundary of the density curve is determined according to the following expression:
DEN high ＝[(DEN lim -DEN low )perc high +DEN lim perc low ]/perc low
wherein DEN high Representing the density curve as the upper boundary of the graph, DEN lim Representing the density value of the limestone skeleton, DEN low Representing the lower boundary of the graph of the density curve, perc high And perc low Representing the upper and lower limits of the graphical numerical range of the compensated neutron curve.
According to one embodiment of the invention, in said step two,
determining an acoustic wave time difference curve to form an upper boundary of a graph according to the upper limit of the porosity of the work area to be analyzed and the acoustic wave time difference value of the limestone framework;
and determining the lower boundary of the sound wave time difference curve according to the upper boundary of the sound wave time difference curve and the numerical range of the compensation neutron curve.
According to one embodiment of the invention, the acoustic moveout curve is determined as an upper boundary of the graph according to the following expression:
AC high ＝β[φ high (AC f -AC lim )+AC lim ]
wherein AC high Representing the upper boundary of the acoustic time difference curve, beta representing the second correction factor, AC lim Representing the acoustic time difference value phi of the limestone skeleton high Represents the upper limit of porosity, AC f Representing the sonic jet lag value of the pore fluid.
According to one embodiment of the invention, the lower graphically boundary of the sonic jet lag curve is determined according to the following expression:
AC low ＝[AC lim perc high -(AC high -AC lim )perc low ]/perc high
wherein AC low Representing the lower boundary of the graph of the acoustic time difference curve, AC high Representing the upper boundary of the graph of the acoustic time difference curve, AC lim Sonic jet lag value, perc, representing limestone skeleton high And perc low Representing the upper and lower limits of the graphical numerical range of the compensated neutron curve.
The lithology recognition method provided by the invention can simply, conveniently, rapidly and accurately realize the recognition of the lithology of the limestone, and is suitable for field application. Compared with the existing method, the method has the advantages that the data processing operation is simpler and more convenient, the limestone lithology identification efficiency can be effectively improved, the limestone lithology can be rapidly identified, and therefore an effective auxiliary means is provided for actual production explanation.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following description will briefly explain the drawings required in the embodiments or the description of the prior art:
fig. 1 is a schematic flow chart of an implementation of a lithology recognition method according to an embodiment of the present invention.
FIG. 2 is a schematic diagram of an implementation flow of determining a density curve into a graphical numerical range according to one embodiment of the invention;
FIG. 3 is a schematic diagram of an implementation flow of determining a graphical range of acoustic moveout curves in accordance with one embodiment of the present invention;
fig. 4 is a schematic diagram of lithology recognition results according to one embodiment of the present invention.
Detailed Description
The following will describe embodiments of the present invention in detail with reference to the drawings and examples, thereby solving the technical problems by applying technical means to the present invention, and realizing the technical effects can be fully understood and implemented accordingly. It should be noted that, as long as no conflict is formed, each embodiment of the present invention and each feature of each embodiment may be combined with each other, and the formed technical solutions are all within the protection scope of the present invention.
In the following description, meanwhile, for the purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of the embodiments of the invention. It will be apparent, however, to one skilled in the art that the present invention may be practiced without some of these specific details or in the specific manner described herein.
Additionally, the steps illustrated in the flowcharts of the figures may be performed in a computer system such as a set of computer executable instructions, and although a logical order is illustrated in the flowcharts, in some cases the steps illustrated or described may be performed in an order other than that herein.
Aiming at the problems in the prior art, the invention provides a novel lithology recognition method which can realize rapid recognition of the lithology of limestone.
Fig. 1 shows a schematic implementation flow chart of the lithology recognition method provided in this embodiment.
As shown in fig. 1, in the lithology recognition method provided in this embodiment, a corresponding mapping numerical range is preferably determined in step S101 according to the obtained different logging parameter values of the limestone skeleton of the work area to be analyzed.
Specifically, in this embodiment, the logging parameters preferably include: acoustic time difference, compensating neutrons and density. That is, the logging parameter values acquired in step S101 of the method include the sonic jet lag value AC lim Compensated neutron numerical value CNL lim And density value DEN lim 。
The method preferably utilizes a lithology standard model and combines the measurement results of the pure limestone core test of the work area to be analyzed to obtain different logging parameter values of the limestone framework. Of course, in other embodiments of the present invention, the method may also obtain different logging parameter values of the limestone skeleton of the work area to be analyzed in other reasonable manners, which is not particularly limited by the present invention.
In this embodiment, in step S101, the method preferably uses a compensated neutron numerical value CNL of the limestone skeleton in the area to be analyzed based on a preset initial compensated neutron mapping range lim And respectively determining the ratio of the distance from the upper boundary of the diagrammatical image to the distance from the lower boundary of the diagrammatical image to the whole range of the compensated neutron numerical value of the limestone skeleton, and then determining the diagrammatical numerical value range of the compensated neutron according to the ratio. The upper limit value of the compensating neutron mapping numerical range is the ratio of the compensating neutron numerical value of the limestone framework to the whole range from the upper mapping boundary distance under the mapping numerical range, and the lower limit value of the compensating neutron mapping numerical range is the ratio of the compensating neutron numerical value of the limestone framework to the whole range from the lower mapping boundary distance under the mapping numerical range.
For example, the method may preferably determine in step S101 the ratio of the compensated neutron numerical value of the limestone skeleton to the overall range from the upper boundary distance of the plot under the range of the plotted numerical values according to the following expression:
wherein perc high Representing the ratio of the compensated neutron value of the limestone skeleton to the integral range of the upper boundary distance of the diagramming under the diagramming value range, CNL high And CNL low Respectively representing the upper limit value and the lower limit value of the mapping range of the preset initial compensation neutrons, CNL lim And the compensated neutron value of the limestone skeleton of the work area to be analyzed is represented.
While the compensated neutron value obtained from the limestone skeleton is at a distance perc from the upper boundary of the map high The method may then preferably determine the ratio of the compensated neutron value of the limestone skeleton to the overall range from the mapped lower boundary distance over the mapped value range according to the following expression:
perc low ＝1-perc high (2)
wherein perc low And the ratio of the distance from the lower boundary of the diagramming to the whole range is represented by the compensated neutron value of the limestone framework under the diagramming value range.
Of course, in other embodiments of the present invention, the method may also employ other reasonable ways to obtain the compensated neutron mapping value range according to the actual needs.
Fig. 2 is a schematic diagram of an implementation flow of determining a mapped numerical range of a determined density curve in the present embodiment.
In the present embodiment, as shown in fig. 2, the method preferably first determines a lower boundary of a density curve according to an upper limit of porosity of a work area to be analyzed and a density value of a limestone skeleton in step S201.
Specifically, in the present embodiment, the method preferably determines the lower boundary of the density curve graph according to the following expression:
DEN low ＝α[DEN lim -φ high (DEN lim -DEN f )] (3)
wherein DEN low Representing the lower boundary of the density curve, alpha representing the first correction factor, DEN lim Representing the density value of the limestone skeleton, phi high Represents the upper limit of porosity, DEN f Representing the density value of the pore fluid.
At the lower boundary DEN of the obtained density curve graph low Then, the method will map the lower boundary DEN according to the density curve in step S202 low And compensating the mapped numerical range of the neutron curve to determine the mapped upper boundary of the density curve, thereby obtaining the mapped numerical range of the density curve.
Specifically, in this embodiment, the method may preferably determine the upper boundary of the density curve in step S202 according to the following expression:
DEN high ＝[(DEN lim -DEN low )perc high +DEN lim perc low ]/perc low (4)
wherein DEN high Representing the density curve as the upper boundary of the graph, DEN lim Representing the density value of the limestone skeleton, DEN low Representing the lower boundary of the graph of the density curve, perc high And perc low Representing the upper and lower limits of the graphical numerical range of the compensated neutron curve.
Of course, in other embodiments of the invention, the method may also determine the mapped numerical range of the density curve in other reasonable manners, and the invention is not limited thereto.
Fig. 3 is a schematic diagram of an implementation flow of determining a graphical numerical range of a sound wave time difference curve in the present embodiment.
As shown in FIG. 3, when determining the range of the acoustic time difference curve, the method provided in this embodiment first determines in step S301 the upper limit phi of the porosity of the work area to be analyzed high Sonic time difference value AC of limestone skeleton lim And determining the acoustic wave time difference curve to form an upper boundary of the graph.
For example, the method may preferably determine the upper boundary of the sonic jet lag curve in step S301 according to the following expression:
AC high ＝β[φ high (AC f -AC lim )+AC lim ] (5)
wherein AC high Representing the upper boundary of the acoustic time difference curve, beta representing the second correction factor, AC lim Representing the acoustic time difference value phi of the limestone skeleton high Represents the upper limit of porosity, AC f Representing the sonic jet lag value of the pore fluid.
At the upper boundary AC of the obtained sonic jet lag curve high The method then generates an upper boundary AC in step S302 according to the above-mentioned sonic jet lag curve high And compensating the mapped numerical range of the neutron curve to determine a mapped lower boundary of the acoustic moveout curve.
For example, in this embodiment, the method may preferably determine the lower boundary of the acoustic moveout curve in step S302 according to the following expression:
AC low ＝[AC lim perc high -(AC high -AC lim )perc low ]/perc high (6)
wherein AC low Representing the lower boundary of the graph of the acoustic time difference curve, AC high Representing the upper boundary of the graph of the acoustic time difference curve, AC lim Sonic jet lag value, perc, representing limestone skeleton high And perc low Representing the upper and lower limits of the graphical numerical range of the compensated neutron curve.
Of course, in other embodiments of the present invention, the method may also determine the range of values of the sonic moveout curve in other reasonable manners, and the present invention is not limited thereto.
Meanwhile, it should be noted that, in different embodiments of the present invention, specific values of the first correction factor and the second correction factor may be configured to different reasonable values according to actual needs, and the present invention is not limited to specific values of the first correction factor and the second correction factor.
For example, the value AC of the acoustic time difference of the rock skeleton of the work area to be analyzed lim At 47.5us/ft, the compensated neutron value CNL lim A density value DEN of 0 lim 2.71g/cm 3 . Presetting an upper limit value CNL of a neutron mapping range of initial compensation high And a lower limit value CNL low 45 and-15, respectively. Then according to the expression (1), the ratio (i.e. the upper limit value of the compensated neutron mapping numerical range) perc of the distance from the upper boundary of the mapping to the whole range of the compensated neutron mapping numerical range of the limestone skeleton can be obtained high The value of (2) is 0.75. According to the expression (2), the ratio (i.e. the lower limit value of the compensated neutron mapping numerical range) perc of the distance from the lower boundary of the mapping to the whole range of the compensated neutron mapping numerical range of the limestone skeleton can be obtained low The value of (2) is 0.25.
For the range of values of the density curve, it is assumed that the upper limit of the porosity phi of the work area to be analyzed high 45% of the upper limit phi of the porosity high Corresponding density value DEN f Is 1g/cm 3 If the value of the first correction factor alpha is 1, the lower boundary DEN of the density curve graph can be determined according to the expression (3) low The value of (C) is 1.96g/cm 3 The upper boundary DEN of the density curve can be further determined according to the expression (4) high The value of (C) is 2.96g/cm 3 。
For the range of values of the acoustic moveout curve, the upper porosity limit φ is assumed high Corresponding sonic time difference value AC f And the second correction factor beta is 1.2, then the upper boundary AC of the acoustic wave time difference curve graph can be determined according to the expression (5) high The value of (2) is 120us/ft, and the lower boundary AC of the sound wave time difference curve graph can be further determined according to the expression (6) low Is 20us/ft.
Referring to fig. 1 again, in this embodiment, after determining the mapping value range corresponding to the logging parameter value, the method performs mapping operation on each logging parameter curve of the work area to be analyzed according to the mapping data value range obtained in step S101 in step S102.
Subsequently, the method determines in step S103 whether or not the curves obtained through the above-described mapping operation have an overlapping region. Wherein if there is an overlap region, the method may also determine the lithology of this overlap region as limestone in step S104. Fig. 4 is a diagram showing a lithology recognition result in the present embodiment.
It should be noted that, in other embodiments of the present invention, the method may also use only at least two of the above logging parameters for lithology recognition according to actual needs, which is not specifically limited by the present invention.
From the above description, it can be seen that the lithology recognition method provided by the invention can simply, rapidly and accurately recognize lithology of limestone, and is suitable for field application. Compared with the existing method, the method has the advantages that the data processing operation is simpler and more convenient, the limestone lithology identification efficiency can be effectively improved, the limestone lithology can be rapidly identified, and therefore an effective auxiliary means is provided for actual production explanation.
It is to be understood that the disclosed embodiments are not limited to the specific structures or process steps disclosed herein, but are intended to extend to equivalents of these features as would be understood by one of ordinary skill in the relevant arts. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting.
Reference in the specification to "one embodiment" or "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment of the invention. Thus, the appearances of the phrase "one embodiment" or "an embodiment" in various places throughout this specification are not necessarily all referring to the same embodiment.
While the above examples are intended to illustrate the principles of the invention in one or more applications, it will be apparent to those skilled in the art that various modifications in form, use and details of implementation may be made without departing from the principles and concepts of the invention. Accordingly, the invention is defined by the appended claims.
Claims (8)
1. A lithology recognition method, the method comprising:
step one, respectively determining corresponding imaging numerical ranges according to obtained different logging parameter values of a limestone skeleton of a work area to be analyzed, wherein the logging parameters comprise at least two of the following items: acoustic time difference, compensating neutrons and density;
step two, mapping operation is carried out on each logging parameter curve of the work area to be analyzed according to the mapping numerical range;
judging whether an overlapping area exists in the curve obtained through the mapping operation, wherein if so, judging the lithology of the overlapping area as limestone;
in the second step, determining a lower boundary of a density curve graph according to the upper limit of the porosity of the work area to be analyzed and the density value of the limestone skeleton; determining a density curve mapping upper boundary according to the density curve mapping lower boundary and the compensation neutron curve mapping numerical range;
determining the lower graphically boundary of the density curve according to the following expression:
DEN low ＝α[DEN lim -φ high (DEN lim -DEN f )]
wherein DEN low Representing the lower boundary of the density curve, alpha representing the first correction factor, DEN lim Representing the density value of the limestone skeleton, phi high Represents the upper limit of porosity, DEN f The density value representing the pore fluid determines the upper boundary of the graph of the density curve according to the following expression:
DEN high ＝[(DEN lim -DEN low )perc high +DEN lim perc low ]/perc low
wherein DEN high Representing the density curve as the upper boundary of the graph, DEN lim Representing the density value of the limestone skeleton, DEN low Representing the lower boundary of the graph of the density curve, perc high And perc low Representing the upper and lower limits of the graphical numerical range of the compensated neutron curve.
2. The method of claim 1, wherein in step one, different logging parameter values of the limestone framework are determined using a lithology standard model in combination with pure limestone core test measurements of the work area to be analyzed.
3. The method according to claim 1, wherein in the second step,
based on a preset initial compensation neutron mapping range, respectively determining the ratio of the compensation neutron numerical value of the limestone framework to the whole range of the mapping upper boundary distance and the mapping lower boundary distance under the mapping numerical value range according to the compensation neutron numerical value, and determining the compensation neutron mapping numerical value range according to the ratio.
4. A method according to claim 3, wherein the ratio of the compensated neutron value of the limestone skeleton to the total range from the upper graphically boundary distance in the graphically numeric range is determined according to the following expression:
wherein perc high Representing the ratio of the compensated neutron value of the limestone skeleton to the integral range of the upper boundary distance of the diagramming under the diagramming value range, CNL high And CNL low Respectively representing the upper limit value and the lower limit value of the mapping range of the preset initial compensation neutrons, CNL lim And the compensated neutron value of the limestone skeleton of the work area to be analyzed is represented.
5. The method of claim 3 or 4, wherein the ratio of the compensated neutron numerical value of the limestone skeleton to the total range from the mapped lower boundary distance over the mapped numerical range is determined according to the following expression:
perc low ＝1-perc high
wherein perc low Representing the ratio of the compensated neutron value of the limestone skeleton to the overall range of the mapped lower boundary distance under the mapped numerical range, perc high The compensated neutron value of the limestone skeleton is represented by the ratio of the distance from the upper boundary of the diagramming to the whole range in the diagramming value range.
6. The method according to claim 1, wherein in the second step,
determining an acoustic wave time difference curve to form an upper boundary of a graph according to the upper limit of the porosity of the work area to be analyzed and the acoustic wave time difference value of the limestone framework;
and determining the lower boundary of the sound wave time difference curve according to the upper boundary of the sound wave time difference curve and the numerical range of the compensation neutron curve.
7. The method of claim 6, wherein the sonic jet lag curve is graphically bounded according to the following expression:
AC high ＝β[φ high (AC f -AC lim )+AC lim ]
wherein AC high Representing the upper boundary of the acoustic time difference curve, beta representing the second correction factor, AC lim Representing the acoustic time difference value phi of the limestone skeleton high Represents the upper limit of porosity, AC f Representing the sonic jet lag value of the pore fluid.
8. The method of claim 6 or 7, wherein the sonic jet lag curve is determined as a lower graphically boundary according to the following expression:
AC low ＝[AC lim perc high -(AC high -AC lim )perc low ]/perc high
wherein AC low Representing the lower boundary of the graph of the acoustic time difference curve, AC high Representing the upper boundary of the graph of the acoustic time difference curve, AC lim Sonic jet lag value, perc, representing limestone skeleton high And perc low Representing the upper and lower limits of the graphical numerical range of the compensated neutron curve.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910463463.4A CN112012726B (en) | 2019-05-30 | 2019-05-30 | Lithology recognition method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910463463.4A CN112012726B (en) | 2019-05-30 | 2019-05-30 | Lithology recognition method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112012726A CN112012726A (en) | 2020-12-01 |
| CN112012726B true CN112012726B (en) | 2023-12-12 |
Family
ID=73500890
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910463463.4A Active CN112012726B (en) | 2019-05-30 | 2019-05-30 | Lithology recognition method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112012726B (en) |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101906963A (en) * | 2010-07-23 | 2010-12-08 | 中国石油化工集团公司 | Method for determining saturation by using C/O and formation macro capture cross-section intersecting technique |
| CN102071928A (en) * | 2009-11-25 | 2011-05-25 | 中国石油天然气股份有限公司 | Three-dimensional space volcanic lithology identification method |
| CN102606150A (en) * | 2012-03-30 | 2012-07-25 | 中国海洋石油总公司 | Method and system for identifying fragmental rock lithological characters |
| CN103850679A (en) * | 2014-04-01 | 2014-06-11 | 北京师范大学 | Method for reconstructing interval transit time curve by virtue of multiple logging curves |
| CN104213899A (en) * | 2013-06-04 | 2014-12-17 | 中国石油化工股份有限公司 | Logging identifying method of stratum rock framework |
| CN104514553A (en) * | 2013-09-29 | 2015-04-15 | 中国石油化工股份有限公司 | Recognition method of gas reservoir in lithic sandstone layer and application thereof |
| CN104989392A (en) * | 2015-07-10 | 2015-10-21 | 中国石油天然气股份有限公司 | Lithology identification method |
| CN105114067A (en) * | 2015-08-26 | 2015-12-02 | 中国石油天然气股份有限公司 | Lithology logging facies method |
| CN105275456A (en) * | 2014-06-23 | 2016-01-27 | 中国石油化工股份有限公司 | Method for identifying high-quality shale by using logging information |
| CN106285623A (en) * | 2015-06-08 | 2017-01-04 | 中国石油化工股份有限公司 | Determine the method and system of total content of organic carbon |
| CN106988737A (en) * | 2017-04-28 | 2017-07-28 | 中国石油大港油田勘探开发研究院 | A kind of method that utilization lithology combination recognizes sedimentary facies |
| CN108957587A (en) * | 2018-07-10 | 2018-12-07 | 中国石油天然气股份有限公司 | Igneous rock lithology identification method and device |
| CN109387869A (en) * | 2017-08-04 | 2019-02-26 | 中国石油化工股份有限公司 | Compact limestone Lithofacies Identification method and system |
-
2019
- 2019-05-30 CN CN201910463463.4A patent/CN112012726B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102071928A (en) * | 2009-11-25 | 2011-05-25 | 中国石油天然气股份有限公司 | Three-dimensional space volcanic lithology identification method |
| CN101906963A (en) * | 2010-07-23 | 2010-12-08 | 中国石油化工集团公司 | Method for determining saturation by using C/O and formation macro capture cross-section intersecting technique |
| CN102606150A (en) * | 2012-03-30 | 2012-07-25 | 中国海洋石油总公司 | Method and system for identifying fragmental rock lithological characters |
| CN104213899A (en) * | 2013-06-04 | 2014-12-17 | 中国石油化工股份有限公司 | Logging identifying method of stratum rock framework |
| CN104514553A (en) * | 2013-09-29 | 2015-04-15 | 中国石油化工股份有限公司 | Recognition method of gas reservoir in lithic sandstone layer and application thereof |
| CN103850679A (en) * | 2014-04-01 | 2014-06-11 | 北京师范大学 | Method for reconstructing interval transit time curve by virtue of multiple logging curves |
| CN105275456A (en) * | 2014-06-23 | 2016-01-27 | 中国石油化工股份有限公司 | Method for identifying high-quality shale by using logging information |
| CN106285623A (en) * | 2015-06-08 | 2017-01-04 | 中国石油化工股份有限公司 | Determine the method and system of total content of organic carbon |
| CN104989392A (en) * | 2015-07-10 | 2015-10-21 | 中国石油天然气股份有限公司 | Lithology identification method |
| CN105114067A (en) * | 2015-08-26 | 2015-12-02 | 中国石油天然气股份有限公司 | Lithology logging facies method |
| CN106988737A (en) * | 2017-04-28 | 2017-07-28 | 中国石油大港油田勘探开发研究院 | A kind of method that utilization lithology combination recognizes sedimentary facies |
| CN109387869A (en) * | 2017-08-04 | 2019-02-26 | 中国石油化工股份有限公司 | Compact limestone Lithofacies Identification method and system |
| CN108957587A (en) * | 2018-07-10 | 2018-12-07 | 中国石油天然气股份有限公司 | Igneous rock lithology identification method and device |
Non-Patent Citations (1)
| Title |
|---|
| 大宁—吉县上古生界主要砂体储层物性及含气特征;董良等;《能源与环保》;第第40卷卷(第第4期期);第87-94页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112012726A (en) | 2020-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11016214B2 (en) | Dolomite reservoir prediction method and system based on well and seismic combination, and storage medium | |
| CN103775057B (en) | Method and device for identifying effective reservoir of tight oil and gas reservoir | |
| CN106372402B (en) | The parallel method of fuzzy region convolutional neural networks under a kind of big data environment | |
| WO2017114443A1 (en) | Method and device for determining dissolution degree of karst reservoir | |
| CN109100793B (en) | Method for quantitatively analyzing influence of fracture factors on reservoir | |
| CN106526675B (en) | Tomography spatial parameter extraction method | |
| US11873709B2 (en) | Log based diagenetic rock typing and sweet spot identification for tight gas sandstone reservoirs | |
| CN113672853A (en) | Automatic standardized processing method and system for logging curve | |
| Lu et al. | A novel unbalanced weighted KNN based on SVM method for pipeline defect detection using eddy current measurements | |
| CN106446514A (en) | Fuzzy theory and neural network-based well-log facies recognition method | |
| CN110634187A (en) | House point cloud model generation method and device based on house type graph | |
| CN112012726B (en) | Lithology recognition method | |
| CN106778794B (en) | Computer imaging control method based on image recognition | |
| CN105700016A (en) | Organic bank reservoir earthquake prediction method | |
| CN104121009B (en) | Method and system for generating temperature and pressure curves | |
| CN109669211B (en) | Method and device for predicting seabed lithologic thin sand body reservoir, storage medium and terminal | |
| CN110673215A (en) | Automatic complex lithology distinguishing method based on intersection graph and Fisher distinguishing | |
| CN106568918B (en) | Shale organic carbon content TOC prediction method | |
| CN106250458A (en) | A kind of Unit Cell point geometry information extracting method | |
| CN104314558B (en) | Method for judging reservoir effectiveness by Stoneley wave energy loss degree | |
| CN110895704A (en) | Microorganism dune complex reservoir type identification method and device and storage medium | |
| CN116930023A (en) | Fine interpretation method and device for dense sandstone phase-control classified porosity logging | |
| CN115422830A (en) | Coal-bearing stratum coal seam thickness prediction method and device based on improved Krigin interpolation | |
| CN109116440B (en) | Dense limestone reservoir fracture identification method | |
| CN102096102B (en) | Digital modeling method for seismic exploration |
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 |