CN113126179B - Volcanic pore recognition method - Google Patents
Volcanic pore recognition method Download PDFInfo
- 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
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000001514 detection method Methods 0.000 claims abstract description 8
- 238000005553 drilling Methods 0.000 claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims description 5
- 239000011435 rock Substances 0.000 description 11
- 230000004075 alteration Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 239000004927 clay Substances 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- 230000007062 hydrolysis Effects 0.000 description 2
- 238000006460 hydrolysis reaction Methods 0.000 description 2
- 238000002386 leaching Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000945 filler Substances 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- 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
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.
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)
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)
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. |
-
2019
- 2019-12-31 CN CN201911420004.4A patent/CN113126179B/en active Active
Patent Citations (3)
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 |