CN113126179B - Volcanic pore recognition method - Google Patents

Volcanic pore recognition method Download PDF

Info

Publication number
CN113126179B
CN113126179B CN201911420004.4A CN201911420004A CN113126179B CN 113126179 B CN113126179 B CN 113126179B CN 201911420004 A CN201911420004 A CN 201911420004A CN 113126179 B CN113126179 B CN 113126179B
Authority
CN
China
Prior art keywords
formula
porosity
volcanic
substituting
data
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
Application number
CN201911420004.4A
Other languages
Chinese (zh)
Other versions
CN113126179A (en
Inventor
陈宏飞
杜雪彪
黄伟强
郑爱萍
高军
巨世昌
李宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201911420004.4A priority Critical patent/CN113126179B/en
Publication of CN113126179A publication Critical patent/CN113126179A/en
Application granted granted Critical
Publication of CN113126179B publication Critical patent/CN113126179B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The application provides a volcanic air hole identification method, which comprises the following steps: acquiring drilling detection data X and Y, wherein X is acoustic time difference logging data, and Y isVolumetric density log data; obtaining acoustic time difference porosity from XObtaining bulk Density porosity from YWill beAndsubstituting formula 1: establishing an air hole identification parameter curve which takes a well depth value and an air hole size value as coordinate axes according to a formula 1; and obtaining visible atmospheric hole data according to the pore identification parameter curve. The technical scheme provided by the application can solve the problem that the visible atmospheric holes cannot be distinguished by the method in the prior art.

Description

Volcanic pore recognition method
Technical Field
The application relates to the technical field of volcanic rock measurement, in particular to a volcanic rock pore identification method.
Background
With the development of oil and gas exploration technology, volcanic oil and gas reservoirs are regarded as important resource exploration fields, and attention of a plurality of students is drawn. Volcanic hydrocarbon reservoirs are various in reservoir space types and complex in diagenetic effect. The filling effect is one of important volcanic rock diagenetic effects, and is mainly represented by filling volcanic rock pores and cracks with secondary minerals, namely pore and crack filling materials are formed. Volcanic pore and fracture fillers are the products of paleo-fluid activity, formed by precipitation of supersaturated ions in the fluid under suitable physicochemical conditions.
Volcanic is a product formed by volcanic eruption, the original volcanic has poor reservoir space connectivity, the reservoir seepage capability is poor, and an effective reservoir does not develop; volcanic eruption intermittent or later lifting suffers from weathering and erosion, volcanic is changed, reservoir space is developed, connectivity is good, primary holes, erosion holes and cracks coexist, a hole-seam dual-medium reservoir is mainly used, and reservoir heterogeneity is strong. The complete volcanic rock weathering crust formed after volcanic rock alteration 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 undisturbed volcanic rock are arranged from the top of the weathering crust downwards in sequence. The alteration degree is weakened from the top of the weathered shell downwards in sequence, the weathered clay layer and the hydrolysis zone are mainly made of clay, the pores are mainly non-communicated invalid pores, an effective reservoir cannot be formed, and the effective reservoir mainly develops in the leaching zone and the disintegration zone I. Wherein, the etching zone mainly develops a double-medium reservoir layer with dissolving holes and cracks, and the disintegrating zone I mainly develops a crack and a dissolving Kong Shuangchong medium reservoir layer. The characteristics of lithology, physical properties, fluid seepage and the like of volcanic reservoirs with different alteration degrees are greatly different, so that the logging response characteristics are greatly different. As the volcanic rock alteration degree is enhanced, logging values such as acoustic time difference, compensated neutrons, natural gamma and the like become larger, logging values such as resistivity, lithologic density and the like become smaller, and reservoir physical properties reflected by logging response become better as the alteration degree is enhanced.
Volcanic pores are of many kinds and can be classified into visible atmospheric pores and invisible pores according to size. The visible pores of the volcanic are distributed in the volcanic in an isolated and non-communicated form, so that the identification of the visible pores of the volcanic is difficult in the field of well logging, the research is concentrated on the analysis and identification of the overall porosity of the volcanic, the identification degree of the visible atmospheric pores is low, and the degree of fine distinction is not achieved. Therefore, the prior art lacks a method for identifying the visible pores of volcanic rock capable of being rapidly identified.
Disclosure of Invention
The application provides a volcanic air hole identification method, which aims to solve the problem that visible air holes cannot be distinguished by the method in the prior art.
The application provides a volcanic air hole identification method, which comprises the following steps:
acquiring 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 beAnd->Substituting formula 1:
establishing an air hole identification parameter curve which takes a well depth value and an air hole size value as coordinate axes according to a formula 1;
and obtaining visible atmospheric hole data according to the pore identification parameter curve.
Further, obtaining acoustic time difference porosity from XThe method specifically comprises the following steps:
substituting X into equation 2:
where a1 and b1 are coefficients.
Further, bulk density porosity is obtained from YThe method specifically comprises the following steps:
substituting Y into equation 3:
where a2 and b2 are coefficients.
Further, the acoustic time difference porosity is obtained from X before X is substituted into formula 2Further comprises:
acquiring a reserve value of the well;
correcting a1 and b1 according to the reserve value;
substituting the corrected a1 and b1 into formula 2.
Further, the bulk density porosity is obtained from Y before substituting Y into equation 3Further comprises:
acquiring a reserve value of the well;
correcting a2 and b2 according to the reserve value;
substituting the corrected a2 and b2 into formula 3.
By applying the technical scheme of the application, through the identification method, the acoustic time difference porosity is acquired by using the drilling detection data X, the volume density porosity is acquired by using Y, the data are substituted into the formula 1 provided by the application, the air hole identification parameter curve is established by using the formula 1, and the visible atmospheric hole data can be acquired according to the curve. The identification method provided by the application can further obtain the data of the visible atmospheric holes on the basis that only the total pore data can be obtained in the prior art, so that the azimuth and the size of the visible atmospheric holes can be determined, the comprehensiveness of the detection result is improved, and the identification method has important significance in maximally finding the distribution rule of the oil and gas reservoirs.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the application. In the drawings:
fig. 1 shows a schematic diagram of a visible atmospheric pore distribution curve provided according to an embodiment of the present application.
Detailed Description
The following description of the embodiments of the present application will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present application, but not all embodiments. The following description of at least one exemplary embodiment is merely exemplary in nature and is in no way intended to limit the application, its application, or uses. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
The embodiment of the application provides a volcanic air hole identification method, which comprises the following steps:
step 1, acquiring drilling detection data X and Y, wherein X is acoustic time difference logging data and Y is volume density logging data;
step 2, obtaining the acoustic time difference porosity according to X
Step 3, obtaining the volume density porosity according to Y
Step 4, willAnd->Substituting formula 1:
wherein QKZS in formula 1 represents the pore recognition parameter.
And 5, establishing an air hole identification parameter curve which takes the well depth value and the air hole size value as coordinate axes according to the formula 1.
Specifically, the well depth value is taken as the ordinate of the curve, and the pore size value is taken as the abscissa of the curve.
And step 6, obtaining visible atmospheric hole data according to the pore identification parameter curve. Specifically, the region on the right side of the ordinate axis, which is the visible atmospheric holes, having the pore size value larger than 0 is the distribution region of the visible atmospheric holes.
By applying the technical scheme of the application, through the identification method, the acoustic time difference porosity is acquired by using the drilling detection data X, the volume density porosity is acquired by using Y, the data are substituted into the formula 1 provided by the application, the air hole identification parameter curve is established by using the formula 1, and the visible atmospheric hole data can be acquired according to the curve. The identification method provided by the application can further obtain the data of the visible atmospheric holes on the basis that only the total pore data can be obtained in the prior art, so that the azimuth and the size of the visible atmospheric holes can be determined, the comprehensiveness of the detection result is improved, and the identification method has important significance in maximally finding the distribution rule of the oil and gas reservoirs.
In the present embodiment, when step 2 is performed, the acoustic time difference porosity is obtained from XThe method specifically comprises the following steps:
substituting X into equation 2:
where a1 and b1 are coefficients.
In this embodiment, the bulk density void is obtained from Y while step 3 is performedDegree ofThe method specifically comprises the following steps:
substituting Y into equation 3:
where a2 and b2 are coefficients. a1, b1, a2, and b2 can be determined by bulk density porosity and acoustic moveout porosity.
Specifically, the acoustic time difference porosity is obtained from X before X is substituted into formula 2Further comprises:
acquiring a reserve value of the well;
correcting a1 and b1 according to the reserve value;
substituting the corrected a1 and b1 into formula 2.
By correcting a1 and b1 with the reserve value of the well, the accuracy of equation 2 can be improved, further improving the accuracy of the visible atmospheric data.
Specifically, the bulk density porosity is obtained from Y before substituting Y into equation 3Further comprises:
acquiring a reserve value of the well;
correcting a2 and b2 according to the reserve value;
substituting the corrected a2 and b2 into formula 3.
By correcting a2 and b2 with the reserve value of the well, the accuracy of equation 3 can be improved, further improving the accuracy of the visible atmospheric data.
In order to facilitate the understanding of the present application, the following description is made in connection with experiments:
1. taking visible atmospheric pore identification of a carbon volcanic lava reservoir of a certain oil field in Xinjiang as an example, using J001 well drilling logging data to select sonic time difference logging data and volume density logging data of a volcanic lava section, wherein the reservoir section is in a depth range of 1700-1745 m.
2. A1=0.535, b1= -26.56 and a2= -60.015, b2=164.2 were determined from the above data, and the porosity was calculatedAnd->Further, the volcanic visible atmospheric pore identification index QKZS is obtained, and the specific calculation formula is as follows:
Φ Δt =0.535 x-26.56 equation 2;
Φ ρ -60.015 x y+164.2 formula 3;
QKZS=2*(Φ ρΔt )/(Φ ρΔt ) Equation 1;
the data obtained are shown in table 1:
TABLE 1 volcanic visible atmospheric identification index QKZS calculation data sheet
3. Establishing a volcanic lava visible air hole identification 'QKZS' curve according to a series of QKZS values calculated in the formula 1;
4. as shown in fig. 1, the data obtained by cutting a part of depth in fig. 1 is marked with 0 line of a longitudinal scale in a graph (right graph in fig. 1), and a curve part larger than 0 is filled, namely a distribution area of visible atmospheric holes of the lava.
After the data are obtained, imaging logging is measured on the J001 well site, and the data obtained by the imaging logging are compared with the data obtained by the application, so that the identification of the atmospheric holes is accurate, and the accuracy can reach 90%. The comparison result shows that on the premise of not depending on imaging logging, the effective identification of the visible atmospheric holes of the volcanic lava can be realized by only using the identification method provided by the application, the use cost is low, and the method can be widely popularized in the area.
The identification method provided by the application is an important expansion of conventional logging means in the volcanic reservoir geological field, and has the characteristics of low cost, obvious results and definite effect directivity. And provides help for innovatively developing and further developing the volcanic lava pore identification logging instrument.
In the process of volcanic reservoir research, volcanic lava pores can be effectively identified by utilizing conventional logging data, and an identification model of volcanic lava visible atmospheric pores is established from the aspect of reaction mechanism of target layer pores in logging, so that a corresponding discrimination curve of the model is utilized to judge whether volcanic lava develops visible atmospheric pores or not, and a basis is provided for volcanic rock exploration.
The method can be applied to the volcanic reservoir research field, can be used for rapidly distinguishing the volcanic lava air holes, can be used for researching the precision, and can also provide a reference for comprehensively researching the volcanic rocks. Drilling without imaging logging is difficult in judging visible atmospheric pores of the volcanic reservoir, and the method can effectively improve comprehensive research precision of the volcanic reservoir and has important significance for maximally finding the distribution rule of the oil and gas reservoirs.
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 exemplary embodiments according to the present application. As used herein, the singular is also intended to include the plural unless the context clearly indicates otherwise, and furthermore, it is to be understood that the terms "comprises" and/or "comprising" when used in this specification are taken to specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof.
The relative arrangement of the components and steps, numerical expressions and numerical values set forth in these embodiments do not limit the scope of the present application unless it is specifically stated otherwise. Meanwhile, it should be understood that the sizes of the respective parts shown in the drawings are not drawn in actual scale for convenience of description. Techniques, methods, and apparatus known to one of ordinary skill in the relevant art may not be discussed in detail, but should be considered part of the specification where appropriate. In all examples shown and discussed herein, any specific values should be construed as merely illustrative, and not a limitation. Thus, other examples of the exemplary embodiments may have different values. It should be noted that: like reference numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, no further discussion thereof is necessary in subsequent figures.
In the description of the present application, it should be understood that the azimuth or positional relationships indicated by the azimuth terms such as "front, rear, upper, lower, left, right", "lateral, vertical, horizontal", and "top, bottom", etc., are generally based on the azimuth or positional relationships shown in the drawings, merely to facilitate description of the present application and simplify the description, and these azimuth terms do not indicate and imply that the apparatus or elements referred to must have a specific azimuth or be constructed and operated in a specific azimuth, and thus should not be construed as limiting the scope of protection of the present application; the orientation word "inner and outer" refers to inner and outer relative to the contour of the respective component itself.
Spatially relative terms, such as "above … …," "above … …," "upper surface at … …," "above," and the like, may be used herein for ease of description to describe one device or feature's spatial location relative 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 in use or operation in addition to the orientation depicted in the figures. For example, if the device in the figures is turned over, elements described as "above" or "over" other devices or structures would then be oriented "below" or "beneath" the other devices or structures. Thus, the exemplary term "above … …" may include both orientations of "above … …" and "below … …". The device may also be positioned in other different ways (rotated 90 degrees or at other orientations) and the spatially relative descriptors used herein interpreted accordingly.
In addition, the terms "first", "second", etc. are used to define the components, and are only for convenience of distinguishing the corresponding components, and the terms have no special meaning unless otherwise stated, and therefore should not be construed as limiting the scope of the present application.
The above description is only of the preferred embodiments of the present application and is not intended to limit the present application, but various modifications and variations can be made to the present application by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the protection scope of the present application.

Claims (1)

1. A volcanic pore identification method, characterized in that the identification method comprises:
acquiring drilling detection data X and Y, wherein X is acoustic time difference logging data and Y is volume density logging data;
acquiring acoustic time difference porosity phi from X △t
Obtaining bulk density porosity phi from Y ρ
Will phi △t And phi ρ Substituting the air hole identification parameter QKZS into the formula 1 to obtain the air hole identification parameter QKZS:
QKZS=2*(φ ρΔt )/(φ ρΔt ) Equation 1;
establishing a pore identification parameter curve which takes a well depth value as an ordinate and takes a Kong Shibie parameter QKZS as an abscissa according to a formula 1;
obtaining visible atmospheric hole data according to the air hole identification parameter curve;
acquiring acoustic time difference porosity phi from X △t The method specifically comprises the following steps:
substituting X into equation 2:
φ △t =a1 x+b1 formula 2;
wherein a1 and b1 are coefficients;
obtaining bulk density porosity phi from Y ρ The method specifically comprises the following steps:
substituting Y into equation 3:
φ ρ =a2×y+b2 formula 3;
wherein a2 and b2 are coefficients;
acquiring acoustic time difference porosity phi from X before substituting X into formula 2 △t Further comprises:
acquiring a core analysis porosity value of a well;
correcting a1 and b1 according to the core analysis porosity value;
substituting the corrected a1 and b1 into formula 2;
obtaining the bulk density porosity phi from Y before substituting Y into equation 3 ρ Further comprises:
acquiring a core analysis porosity value of a well;
correcting a2 and b2 according to the core analysis porosity value;
substituting the corrected a2 and b2 into formula 3.
CN201911420004.4A 2019-12-31 2019-12-31 Volcanic pore recognition method Active CN113126179B (en)

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 CN113126179A (en) 2021-07-16
CN113126179B true 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 (3)

* Cited by examiner, † Cited by third party
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
CN108802192A (en) * 2017-05-03 2018-11-13 中国石油化工股份有限公司 A kind of calcarenaceous sandstone reservoir pore space kind identification method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
MX2016007733A (en) * 2014-01-14 2016-09-07 Halliburton Energy Services Inc Effective porosity determination for tight gas formations.

Patent Citations (3)

* Cited by examiner, † Cited by third party
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
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
CN113126179A (en) 2021-07-16

Similar Documents

Publication Publication Date Title
AU2020102001A4 (en) Method for Identifying and Describing Deep Carbonate Karst Structures
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
US20210404331A1 (en) Fine identification method of tight reservoir fracture based on conventional logging data
Pechnig et al. Integrated log interpretation in the German Continental Deep Drilling Program: Lithology, porosity, and fracture zones
CN102128025B (en) Method for determining density of salt-gypsum layer deflecting drilling fluid
CN104213899B (en) A kind of Logging Identification Method of formation rock skeleton
MX2012010692A (en) Multi-step borehole correction scheme for multi-component induction tools.
CN105986815A (en) Method for recognizing geological sweet spots in shale stratum
CN107941670B (en) Rock debris porosity determination method
Swarbrick Challenges of porosity based pore pressure prediction
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
Liu et al. Natural fractures in tight gas volcanic reservoirs and their influences on production in the Xujiaweizi depression, Songliao Basin, China
Rizzo et al. Detecting the onset of strain localization using two‐dimensional wavelet analysis on sandstone deformed at different effective pressures
CN106897531A (en) A kind of method for quantitatively evaluating of hyposmosis limestone reservoir permeability
EP1569163A1 (en) Method of determining planar events from borehole or core images
CN113126179B (en) Volcanic pore recognition method
Liu et al. The critical control of arkosic sandstone porosity on deformation band formation: Insights from the Shulu across-fault borehole in the Bohai Bay Basin, China
Duan et al. Architectural characterization of Ordovician fault-controlled paleokarst carbonate reservoirs, Tahe oilfield, China
Zhou et al. Fluid flow through porous sandstone with overprinting and intersecting geological structures of various types
Ning et al. Theory on logging interpretation of igneous rocks and its application
CN105589098B (en) Reservoir inversion method and system for removing transverse influence of sedimentary argillaceous rocks by carbonate rocks
CN106646668B (en) A kind of method for building up of radar logging standard well model
CN112878999B (en) Method and device for calculating water saturation of anisotropic stratum
CN116446861A (en) Method for evaluating underground mountain stratum logging based on elemental analysis

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