CN113126179A - Volcanic pore identification method - Google Patents
Volcanic pore identification method Download PDFInfo
- Publication number
- CN113126179A CN113126179A CN201911420004.4A CN201911420004A CN113126179A CN 113126179 A CN113126179 A CN 113126179A CN 201911420004 A CN201911420004 A CN 201911420004A CN 113126179 A CN113126179 A CN 113126179A
- Authority
- CN
- China
- Prior art keywords
- volcanic
- porosity
- pore
- data
- acquiring
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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
Abstract
The invention provides a volcanic air hole identification method, which comprises the following steps: acquiring well drilling detection data X and Y, wherein X is acoustic time difference logging data, and Y is volume density logging data; obtaining acoustic time difference porosity from X
Obtaining bulk density porosity from Y
Will be provided with
And
substituting formula 1: establishing a pore identification parameter curve respectively taking the well depth value and the pore size value as coordinate axes according to the formula 1; and acquiring visible atmosphere data according to the atmosphere identification parameter curve. Technical method provided by the applicationThe method can solve the problem that the method in the prior art cannot distinguish visible big air holes.
Description
Technical Field
The invention relates to the technical field of volcanic rock measurement, in particular to a volcanic rock air hole identification method.
Background
With the development of oil and gas exploration technology, volcanic oil and gas reservoirs are attracting the attention of a plurality of scholars as an important resource exploration field. The volcanic oil and gas reservoir has various reservoir space types and complex diagenesis. The filling function is one of important volcanic diagenesis, and is mainly characterized in that volcanic air holes and cracks are filled with secondary minerals, namely, air hole and crack fillers are formed. Volcanic pore and crack fillers are the products of ancient fluid activities, and are formed by precipitation of supersaturated ions in a fluid under appropriate physicochemical conditions.
Volcanic rock is a product formed by volcanic eruption, the connectivity of an original volcanic rock storage space is poor, the seepage capability of a reservoir is poor, and an effective reservoir does not develop; the volcanic eruption intermittent period or later lifting is subjected to weathering erosion, volcanic rock is altered, a storage space is developed, the connectivity is good, primary pores, erosion pores and cracks coexist, a pore-crack dual medium storage layer is taken as a main component, and the heterogeneity of the storage layer is strong. The complete volcanic weathering crust formed after the volcanic rock is altered has a six-layer structure, and a weathering clay layer, a hydrolysis zone, a leaching zone, a disintegration zone I, a disintegration zone II and an original volcanic rock are arranged from the top of the weathering crust to the bottom in sequence. The erosion degree is weakened from the top of the weathering crust to the bottom in sequence, the weathering clay layer and the hydrolysis zone are mainly clay, the pores are mainly non-communicated ineffective pores, an effective reservoir cannot be formed, and the effective reservoir mainly develops in the erosion zone and the disintegration zone I. The leaching zone mainly develops the double-medium reservoir of the erosion holes and the cracks, and the disintegration zone I mainly develops the double-medium reservoir of the cracks and the erosion holes. The characteristics of lithology, physical properties, fluid seepage and the like of volcanic rock reservoirs with different alteration degrees are greatly different, so that the logging response characteristics are greatly different. With the enhancement of the volcanic rock alteration degree, logging values such as acoustic wave time difference, compensation neutrons and natural gamma become larger, logging values such as resistivity and lithologic density become smaller, and the physical property of the reservoir layer reflected by logging response becomes better with the enhancement of the alteration degree.
Volcanic rock pores are of various types and can be divided into visible large pores and invisible pores according to size. Most of volcanic visible air holes are distributed in volcanic rock in an isolated and disconnected mode, in the field of well logging, volcanic visible air holes are difficult to identify, at present, research is mostly focused on analysis and identification of the total porosity of volcanic rock, the identification degree of visible air holes is low, and the degree of fine distinguishing is not achieved. Therefore, a method for identifying the visible pores of the volcanic rock, which can be quickly identified, is lacked in the prior art.
Disclosure of Invention
The invention provides a volcanic air hole identification method, which aims to solve the problem that the method in the prior art cannot distinguish visible large air holes.
The invention provides a volcanic air hole identification method, which comprises the following steps:
acquiring well drilling detection data X and Y, wherein X is acoustic time difference logging data, and Y is volume density logging data;
obtaining acoustic time difference porosity from X
Obtaining bulk density porosity from Y
Will be provided with
And
substituting formula 1:
establishing a pore identification parameter curve respectively taking the well depth value and the pore size value as coordinate axes according to the formula 1;
and acquiring visible atmosphere data according to the atmosphere identification parameter curve.
Further, acquiring the time difference porosity of the sound wave according to X
The method specifically comprises the following steps:
substituting X into equation 2:
wherein a1 and b1 are coefficients.
Further, bulk density porosity is obtained from Y
The method specifically comprises the following steps:
substituting Y into equation 3:
wherein a2 and b2 are coefficients.
Further, the acoustic moveout porosity is obtained from X before substituting X into equation 2
Further comprising:
acquiring a reserve value of a well;
correcting a1 and b1 according to the reserve value;
the corrected a1 and b1 are substituted into equation 2.
Further, before substituting Y into equation 3, the bulk density porosity is obtained from Y
Also comprises：
Acquiring a reserve value of a well;
correcting a2 and b2 according to the reserve value;
the corrected a2 and b2 are substituted into equation 3.
By applying the technical scheme of the invention, the drilling detection data X is used for obtaining the acoustic time difference porosity, the volume density porosity is obtained by Y, the data are substituted into the formula 1 provided by the application, the formula 1 is used for establishing a pore identification parameter curve, and the visible big pore data can be obtained according to the curve. Through the identification method provided by the application, the data of the visible big air holes can be further obtained on the basis that only total pore data can be obtained at present, so that the direction and the size of the visible big air holes can be determined, the comprehensiveness of detection results is improved, and the method has important significance for finding the distribution rule of oil and gas reservoirs to the maximum extent.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this application, illustrate embodiments of the invention and, together with the description, serve to explain the invention and not to limit the invention. In the drawings:
fig. 1 shows a schematic view of a visible macropore distribution curve provided in accordance with an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. The following description of at least one exemplary embodiment is merely illustrative in nature and is in no way intended to limit the invention, its application, or uses. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The embodiment of the invention provides a volcanic air hole identification method, which comprises the following steps:
step 1, acquiring well drilling detection data X and Y, wherein X is acoustic time difference logging data, and Y is volume density logging data;
step 2, acquiring the time difference porosity of sound waves according to X
Step 3, obtaining the volume density porosity according to Y
Step 4, mixing
And
substituting formula 1:
wherein QKZS in formula 1 represents the air hole identification parameter.
And 5, establishing a pore identification parameter curve respectively taking the well depth value and the pore size value as coordinate axes according to the formula 1.
Specifically, the value of the well depth is taken as the ordinate of the curve, and the value of the pore size is taken as the abscissa of the curve.
And 6, acquiring visible atmosphere data according to the atmosphere identification parameter curve. Specifically, the visible air holes having a hole size value larger than 0 are defined, that is, the region on the right side of the ordinate axis is the visible air hole distribution region.
By applying the technical scheme of the invention, the drilling detection data X is used for obtaining the acoustic time difference porosity, the volume density porosity is obtained by Y, the data are substituted into the formula 1 provided by the application, the formula 1 is used for establishing a pore identification parameter curve, and the visible big pore data can be obtained according to the curve. Through the identification method provided by the application, the data of the visible big air holes can be further obtained on the basis that only total pore data can be obtained at present, so that the direction and the size of the visible big air holes can be determined, the comprehensiveness of detection results is improved, and the method has important significance for finding the distribution rule of oil and gas reservoirs to the maximum extent.
In the embodiment, when step 2 is executed, the time difference porosity of the sound wave is obtained according to X
The method specifically comprises the following steps:
substituting X into equation 2:
wherein a1 and b1 are coefficients.
In the present embodiment, in performing step 3, the bulk density porosity is obtained from Y
The method specifically comprises the following steps:
substituting Y into equation 3:
wherein a2 and b2 are coefficients. a1, b1, a2 and b2 can be determined by bulk density porosity and differential porosity in sound wave.
Specifically, the acoustic moveout porosity is obtained from X before substituting X into formula 2
Further comprising:
acquiring a reserve value of a well;
correcting a1 and b1 according to the reserve value;
the corrected a1 and b1 are substituted into equation 2.
By correcting a1 and b1 using the reserves of the drilled well, the accuracy of equation 2 can be improved, further improving the accuracy of the visible atmospheric data.
Specifically, bulk density porosity is obtained from Y before substituting Y into equation 3
Further comprising:
acquiring a reserve value of a well;
correcting a2 and b2 according to the reserve value;
the corrected a2 and b2 are substituted into equation 3.
By correcting a2 and b2 using the drilling reserves, the accuracy of equation 3 can be improved, further improving the accuracy of the visible atmospheric data.
For the purpose of facilitating an understanding of the present application, the following description is made in conjunction with experiments:
1. by taking visible atmosphere hole identification of a rock-carbon volcanic lava reservoir in an oil field in Xinjiang as an example, acoustic time difference logging data and volume density logging data of a volcanic lava section are selected by using J001 well drilling logging data, and the depth range of the reservoir section is 1700-1745 m.
2. From the above data, a 1-0.535, b 1-26.56 and a 2-60.015, b 2-164.2 were determined, and the porosity was calculated
And
and then, obtaining the visible atmospheric hole identification index QKZS of the volcanic rock, wherein the specific calculation formula is as follows:
ΦΔt0.535X-26.56 formula 2;
Φρ-60.015 x Y +164.2 formula 3;
QKZS＝2*(Φρ-ΦΔt)/(Φρ+ΦΔt) Formula 1;
the data obtained are shown in table 1:
TABLE 1 volcanic rock visible atmosphere recognition index QKZS calculation data table
3. Establishing a volcanic lava visible pore recognition 'QKZS' curve according to a series of QKZS values calculated by the formula 1;
4. as shown in fig. 1, the data obtained by cutting out a part of the depth in fig. 1 is obtained by defining a 0 line of a longitudinal scale in a graph (a right graph in fig. 1), and a curve part which is larger than 0 is filled, namely a distribution area of visible macropores of the volcanic lava.
After the data are obtained, imaging logging is measured on the J001 well site, the data obtained by imaging logging is compared with the data obtained by the method, the visible macropore identification is accurate, and the accuracy can reach 90%. The comparison result shows that the identification method provided by the application can realize effective identification of visible macropores of the volcanic lava on site without depending on imaging logging, is low in use cost, and can be popularized in a large scale in an area.
The identification method provided by the invention is an important expansion of a conventional logging means in the volcanic reservoir geological field, and has the characteristics of low cost, obvious achievement and definite effect directivity. The method provides help for the capability of innovatively developing and further developing the volcanic lava pore recognition logging instrument.
In the process of volcanic reservoir research, the conventional logging information is utilized to effectively identify the volcanic lava pores, starting from the reconstruction of conventional logging, and considering the reaction mechanism of the target stratum pores in logging, a model for identifying the visible macropores of the volcanic lava is established, and whether the volcanic lava develops the visible macropores or not is judged according to the corresponding judgment curve of the model, so that a basis is provided for volcanic exploration.
The method can be applied to the field of volcanic reservoir research, can be used for rapidly distinguishing the volcanic lava macropores, can be used for researching the precision, and can also provide a reference for comprehensively researching the volcanic. The method can effectively improve the comprehensive research precision of the volcanic reservoir and has important significance for finding the distribution rule of the oil and gas reservoir to the maximum.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments according to the present application. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof, unless the context clearly indicates otherwise.
The relative arrangement of the components and steps, the numerical expressions and numerical values set forth in these embodiments do not limit the scope of the present invention unless specifically stated otherwise. Meanwhile, it should be understood that the sizes of the respective portions shown in the drawings are not drawn in an actual proportional relationship for the convenience of description. Techniques, methods, and apparatus known to those of ordinary skill in the relevant art may not be discussed in detail but are intended to be part of the specification where appropriate. In all examples shown and discussed herein, any particular value should be construed as merely illustrative, and not limiting. Thus, other examples of the exemplary embodiments may have different values. It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, further discussion thereof is not required in subsequent figures.
In the description of the present invention, it is to be understood that the orientation or positional relationship indicated by the orientation words such as "front, rear, upper, lower, left, right", "lateral, vertical, horizontal" and "top, bottom", etc. are usually based on the orientation or positional relationship shown in the drawings, and are only for convenience of description and simplicity of description, and in the case of not making a reverse description, these orientation words do not indicate and imply that the device or element being referred to must have a specific orientation or be constructed and operated in a specific orientation, and therefore, should not be considered as limiting the scope of the present invention; the terms "inner and outer" refer to the inner and outer relative to the profile of the respective component itself.
Spatially relative terms, such as "above … …," "above … …," "above … …," "above," and the like, may be used herein for ease of description to describe one device or feature's spatial relationship to another device or feature as illustrated in the figures. It will be understood that the spatially relative terms are intended to encompass different orientations of the device in use or operation in addition to the orientation depicted in the figures. For example, if a device in the figures is turned over, devices described as "above" or "on" other devices or configurations would then be oriented "below" or "under" the other devices or configurations. Thus, the exemplary term "above … …" can include both an orientation of "above … …" and "below … …". The device may be otherwise variously oriented (rotated 90 degrees or at other orientations) and the spatially relative descriptors used herein interpreted accordingly.
It should be noted that the terms "first", "second", and the like are used to define the components, and are only used for convenience of distinguishing the corresponding components, and the terms have no special meanings unless otherwise stated, and therefore, the scope of the present invention should not be construed as being limited.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (5)
1. A volcanic pore identification method is characterized by comprising the following steps:
acquiring well drilling detection data X and Y, wherein X is acoustic time difference logging data, and Y is volume density logging data;
obtaining acoustic time difference porosity from X
Obtaining bulk density porosity from Y
Will be provided with
And
substituting formula 1:
establishing a pore identification parameter curve respectively taking the well depth value and the pore size value as coordinate axes according to the formula 1;
and acquiring visible atmosphere data according to the atmosphere identification parameter curve.
2. The volcanic pore identification method as claimed in claim 1, wherein the acoustic moveout porosity is obtained from X
The method specifically comprises the following steps:
substituting X into equation 2:
wherein a1 and b1 are coefficients.
3. The volcanic pore identification method as claimed in claim 1, wherein the bulk density porosity is obtained from Y
The method specifically comprises the following steps:
substituting Y into equation 3:
wherein a2 and b2 are coefficients.
4. The method of claim 2The volcanic pore identification method is characterized in that before X is substituted into formula 2, the sound wave time difference porosity is obtained according to X
Further comprising:
acquiring a reserve value of a well;
correcting a1 and b1 according to the reserve value;
the corrected a1 and b1 are substituted into equation 2.
5. The volcanic pore identification method of claim 3, wherein before substituting Y into equation 3, the bulk density porosity is obtained from Y
Further comprising:
acquiring a reserve value of a well;
correcting a2 and b2 according to the reserve value;
the corrected a2 and b2 are substituted into equation 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911420004.4A CN113126179B (en) | 2019-12-31 | 2019-12-31 | Volcanic pore recognition method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911420004.4A CN113126179B (en) | 2019-12-31 | 2019-12-31 | Volcanic pore recognition method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113126179A true CN113126179A (en) | 2021-07-16 |
| CN113126179B CN113126179B (en) | 2023-08-22 |
Family
ID=76769471
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911420004.4A Active CN113126179B (en) | 2019-12-31 | 2019-12-31 | Volcanic pore recognition method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113126179B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6411902B1 (en) * | 1999-04-19 | 2002-06-25 | Michael John Wiltshire | Shale compaction and sonic logs |
| CN105240005A (en) * | 2015-10-28 | 2016-01-13 | 中国石油化工股份有限公司 | Low-hole extra-low-permeability unconventional reservoir dissolution hole identification method |
| US20180163535A1 (en) * | 2014-01-14 | 2018-06-14 | Halliburton Energy Services, Inc. | Effective porosity determination for tight gas formations |
| CN108802192A (en) * | 2017-05-03 | 2018-11-13 | 中国石油化工股份有限公司 | A kind of calcarenaceous sandstone reservoir pore space kind identification method |
-
2019
- 2019-12-31 CN CN201911420004.4A patent/CN113126179B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6411902B1 (en) * | 1999-04-19 | 2002-06-25 | Michael John Wiltshire | Shale compaction and sonic logs |
| US20180163535A1 (en) * | 2014-01-14 | 2018-06-14 | Halliburton Energy Services, Inc. | Effective porosity determination for tight gas formations |
| CN105240005A (en) * | 2015-10-28 | 2016-01-13 | 中国石油化工股份有限公司 | Low-hole extra-low-permeability unconventional reservoir dissolution hole identification method |
| CN108802192A (en) * | 2017-05-03 | 2018-11-13 | 中国石油化工股份有限公司 | A kind of calcarenaceous sandstone reservoir pore space kind identification method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113126179B (en) | 2023-08-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU2020102001A4 (en) | Method for Identifying and Describing Deep Carbonate Karst Structures | |
| CN110412661B (en) | Evaluation method and device for dominant segment cluster of fine-grained rock oil and gas reservoir dessert segment | |
| Dutton et al. | Calcite cement distribution and its effect on fluid flow in a deltaic sandstone, Frontier Formation, Wyoming | |
| Ya et al. | Geologic features and genesis of the barriers and intercalations in carbonates: A case study of the Cretaceous Mishrif Formation, West Qurna oil field, Iraq | |
| CN103670383B (en) | A kind of method and apparatus identifying shale oil deposit Effective Reservoirs | |
| CN111880234B (en) | Fine identification method for tight reservoir fractures based on conventional well logging | |
| CN109061763B (en) | Comprehensive evaluation method for carbonate rock solution reservoir cave logging | |
| CN109116428B (en) | Fracture-cavity carbonate reservoir uncertainty modeling method and device | |
| CN104213899B (en) | A kind of Logging Identification Method of formation rock skeleton | |
| CN107829731B (en) | Clay alteration volcanic porosity correction method | |
| CN109143399B (en) | Method for identifying carbonate rock sequence interface | |
| CN106054254B (en) | A kind of log curve correction method and device | |
| CN107861162B (en) | Microelectrode logging data-based natural crack identification method and system | |
| Liu et al. | Natural fractures in tight gas volcanic reservoirs and their influences on production in the Xujiaweizi depression, Songliao Basin, China | |
| CN103792575B (en) | Densified sandstone air layer identification evaluation method free from influence of frame | |
| Al-Fahmi et al. | New insights on fracture roughness and wall mismatch in carbonate reservoir rocks | |
| Jiang et al. | Analyzing the origin of low resistivity in gas-bearing tight sandstone reservoir | |
| Smith et al. | Discrimination of effective from ineffective porosity in heterogeneous Cretaceous carbonates, Al Ghubar field, Oman | |
| CN115390155A (en) | Well logging interpretation method, device, electronic equipment and medium for heterogeneous reservoir | |
| CN110208874A (en) | A kind of tight sandstone reservoir effective reservoir recognition methods | |
| CN108802192B (en) | Calcium debris sandstone reservoir pore type identification method | |
| CN111537417B (en) | Rock sample pore development condition evaluation method | |
| CN113126179A (en) | Volcanic pore identification method | |
| Zhou et al. | Fluid flow through porous sandstone with overprinting and intersecting geological structures of various types | |
| Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China |
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 |