CN104564042B - A kind of brittle evaluation method of shale reservoir - Google Patents
A kind of brittle evaluation method of shale reservoir Download PDFInfo
- Publication number
- CN104564042B CN104564042B CN201410790723.6A CN201410790723A CN104564042B CN 104564042 B CN104564042 B CN 104564042B CN 201410790723 A CN201410790723 A CN 201410790723A CN 104564042 B CN104564042 B CN 104564042B
- Authority
- CN
- China
- Prior art keywords
- modulus
- density
- young
- msup
- msub
- 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
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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
Abstract
The invention provides a kind of brittle evaluation method of shale reservoir, belongs to unconventionaloil pool reservoir assessment technology field, and this method includes the elastic parameter and brittle mineral parameter for obtaining shale reservoir to be evaluated, and the elastic parameter includes Young's modulus * density;Young's modulus * density and other elastic parameters and/or brittle mineral parameter are subjected to cross analysis, and the codomain of Young's modulus * density is determined according to the result of cross analysis;Young's modulus * density attributes bodies are built according to the wave impedance in length and breadth obtained by prestack inversion;Young's modulus * density attributes bodies are evaluated the fragility of shale reservoir to be evaluated in codomain.The present invention improves the brittle evaluation precision of shale reservoir, realizes to whole research area's fragility planar distribution prediction, is had broad application prospects in shale gas exploration and development field.
Description
Technical field
The present invention relates to a kind of brittle evaluation method of shale reservoir, belongs to unconventionaloil pool reservoir assessment technology field.
Background technology
Shale gas is a kind of unconventional gas resource, because shale belongs to ultralow hole low permeability reservoir, causes its Cheng Zangmo
There is very big difference in formula, gas accumulation mode and development scheme, therefore the exploitation of shale gas must with conventional gas and oil reservoir
Reservoir fracturing improvement must be implemented.Reservoir fracturing improvement is to improve shale gas reservoir permeability, ensure shale gas is smoothly exploited one
Item key measure.There are some researches show in the big region of reservoir fragility, because rock is more easy to be pressed off and be more likely formed seam net ditch
It is logical, so increase reservoir reconstruction volume can improve shale gas recovery ratio.Therefore, it is to influence shale gas to find the high brittle zone of shale
One key factor of yield.
Existing shale reservoir brittleness evaluation method mainly wraps two kinds, and the first is the quartz in rock forming mineral component
The fragility of shale reservoir containing amount instruction is for second based on three dimensional seismic data and utilizes the result of prestack elastic parameter inversion
Characterize shale reservoir fragility.First method is based on well logging or core test, although precision is of a relatively high, can only be reflected
Single wellhole fragility change, it can not meet to whole research area's fragility planar distribution prediction.Second method utilizes prestack elasticity
The inversion result structure brittleness index of parameter, because the existing acquisition to elastic parameter is more difficult, therefore corresponding elasticity ginseng
Several precision is relatively low, it is difficult to realizes the exact evaluation to shale reservoir elasticity.
The content of the invention
The present invention is to solve existing for existing shale reservoir brittleness evaluation technology not meeting to whole research area's fragility
The problem of planar distribution is predicted and evaluation precision is relatively low, and then a kind of brittle evaluation method of shale reservoir is proposed, specifically
Including following technical scheme:
A kind of brittle evaluation method of shale reservoir, including:
The elastic parameter and brittle mineral parameter of shale reservoir to be evaluated are obtained, the elastic parameter includes Young's modulus *
Density;
The Young's modulus * density and other elastic parameters and/or the brittle mineral parameter cross point
Analyse, and the codomain of the Young's modulus * density is determined according to the result of the cross analysis;
Young's modulus * density attributes bodies are built according to the wave impedance in length and breadth obtained by prestack inversion;
The Young's modulus * density attributes body is commented the fragility of the shale reservoir to be evaluated in the codomain
Valency.
In the brittle evaluation method of shale reservoir of the present invention, by the Young's modulus * density with it is other described
Elastic parameter and/or the brittle mineral parameter, which carry out cross analysis, to be included:
Any at least two elastic parameters including the Young's modulus * density are subjected to cross analysis or will be included
Elastic parameter described in any two of Young's modulus * density carries out cross analysis with any one of brittle mineral parameter, and
The codomain of the Young's modulus * density is determined according to the result of the analysis.
In the brittle evaluation method of shale reservoir of the present invention, the elastic parameter also include in length and breadth wave impedance,
Extra large constant * density and modulus of shearing * density is drawn, the brittle mineral parameter contains including quartz content, shale content and calcite
Amount.
In the brittle evaluation method of shale reservoir of the present invention, the Young's modulus is calculated by below equation and obtained
:
E=3k (1- σ)
Wherein, E represents Young's modulus, and k represents bulk modulus;σ represents Poisson's ratio.
In the brittle evaluation method of shale reservoir of the present invention, the Poisson's ratio is calculated by below equation and obtained
:
Wherein, VpRepresent velocity of longitudinal wave, VsRepresent shear wave velocity.
In the brittle evaluation method of shale reservoir of the present invention, the Young's modulus * density attributes body by with
Lower formula, which calculates, to be obtained:
Wherein, ρ represents density, IpRepresent p-wave impedance, IsRepresent S-wave impedance.
In the brittle evaluation method of shale reservoir of the present invention, the brittle mineral parameter contains according to the quartz
Quantitative relationship between amount, shale content and calcite content and the Young's modulus * density obtains.
In the brittle evaluation method of shale reservoir of the present invention, the elastic parameter of the shale reservoir to be evaluated and
Brittle mineral parameter is obtained by the well-log information of the shale reservoir to be evaluated.
The beneficial effects of the invention are as follows:Shale reservoir fragility is evaluated by the product of Young's modulus and density, due to Young
The product of modulus and density can be characterized by the combined result of wave impedance in length and breadth, and wave impedance phase in prestack inversion in length and breadth
It is more accurately, to thereby eliminate caused by density item accuracy is poor to shale reservoir brittleness evaluation precision to density
Influence, so as to improve brittleness evaluation precision, and to whole research area's fragility planar distribution prediction, in shale gas exploration and development
Field has broad application prospects.
Brief description of the drawings
Fig. 1 is the flow chart for showing the brittle evaluation method of shale reservoir in a schematic manner.
Fig. 2 is the figure that crosses of modulus of shearing * density and Young's modulus * density that present embodiment provides.
Fig. 3 is that the Young's modulus * density that present embodiment provides characterizes shale reservoir fragility to be evaluated in shale storage
Response diagram on layer section.
Fig. 4 is the fragility planar distribution figure for the shale reservoir to be measured that present embodiment provides.
Embodiment
In the state of the art, because the extraction in the inversion result to density information needs wide-angle information
(i.e. long geophone offset), and the arrangement of the geophone offset of long length can not undertake in actual production, and wide-angle information earthquake
Data quality is also difficult to ensure that.Therefore, density ask for be the brittle evaluation method of shale reservoir key, while be also ability
The technical barrier of domain for a long time.The technology hair shaft that present embodiment proposes is directed to existing shale reservoir brittleness evaluation side
Method proposes a kind of brittle evaluation method of shale reservoir, passes through poplar because of the not high influence brought to evaluation result of density item precision
The product of family name's modulus and density evaluates fragility, because the product of attribute Young's modulus and density can be directly by the group of wave impedance in length and breadth
Result is closed to be characterized, and wave impedance relative density in prestack inversion is that more accurately, thus eliminating the need because close in length and breadth
Degree item is forbidden the influence brought, so as to improve brittleness evaluation precision.
Shown in brittle evaluation method combination Fig. 1 of shale reservoir that present embodiment proposes, including:
Step 11, the elastic parameter and brittle mineral parameter of shale reservoir to be evaluated are obtained, the elastic parameter includes poplar
Family name's modulus * density.
The elastic parameter can be according to well-log information, if sound wave, density, dipole logging are (if can be rule of thumb without actual measurement
Formula fitting) etc. obtain, corresponding elastic parameter can include:Wave impedance (I in length and breadthp, Is), draw extra large constant * density (λ ρ), cut
Shear modulu * density (μ ρ) and Young's modulus * density (E ρ).
The brittle mineral parameter includes quartz content (Vquartz), shale content (Vclay) and calcite content
(Vcalcite), it can calculate quartz content, shale content and calcite content by existing log, and according to by by institute
The quantitative relationship stated between quartz content, shale content and calcite content and the Young's modulus * density obtains the quartz
Percentage composition, shale content percentage and calcite content percentage.
Step 12, the Young's modulus * density and the elastic parameter and/or the brittle mineral parameter are crossed
Analyze, and the codomain of the Young's modulus * density is determined according to the result of the cross analysis.
The codomain of the Young's modulus * density is can determine by cross analysis, can will preferably include Young's modulus *
Any two of density or three elastic parameters carry out cross analysis, or any two bullet that will include Young's modulus * density
Property parameter and any one brittle mineral parameter carry out cross analysis, it is and true according to the figure that crosses obtained by cross analysis
Determine the codomain where Young's modulus * density.
Step 13, Young's modulus * density attributes bodies are built according to the wave impedance in length and breadth obtained by prestack inversion.
Due to the acquisition that wave impedance can be more accurate during prestack inversion in length and breadth, following side can be preferably used
Method:First based on sound wave, shear wave, density log curve, different angular-trace gathers is then extracted from prestack trace gather, typically
Three different angular-trace gathers are at least wanted, wave resistance in length and breadth can be solved by growing approximate equation (such as Sheuy, 1985) according to Zuo Buni
Anti-, density item.So the data computational accuracy of the brittleness evaluation process of shale reservoir can be improved.Young's modulus therein is table
The rigid physical quantity of rock is levied, can be represented with bulk modulus and Poisson's ratio:
E=3k (1- σ) (1)
Wherein, E represents Young's modulus, and k represents bulk modulus;σ represents Poisson's ratio.The Poisson's ratio is represented by:
Wherein, wherein, VpRepresent velocity of longitudinal wave, VsRepresent shear wave velocity.Formula (2) is substituted into formula (1), then obtained:
Wherein, ρ represents density.ρ, and I are multiplied by formula (3) both sides simultaneouslyp=ρ Vp,Is=ρ Vs, then formula (3) can be write as:
Wherein, IpRepresent p-wave impedance, IsRepresent S-wave impedance.
Formula (4) is then the Young's modulus * density attributes bodies built by wave impedance in length and breadth, due to not including in the attribute volume
Density item ρ, and the wave impedance in length and breadth obtained during prestack inversion has higher precision, therefore pass through the Young mould
Amount * density attributes body, which is evaluated the fragility of shale reservoir to be measured, can also obtain higher rating precision.
Step 14, by the Young's modulus * density attributes body in the codomain to the crisp of the shale reservoir to be evaluated
Property is evaluated.
The Young's modulus * density attributes bodies built according to formula (4), and cross analysis institute is utilized in convolution (2) and formula (3)
Codomain where the Young's modulus * density of determination, is evaluated the fragility of shale reservoir to be measured, brittle so as to obtain shale
Planar distribution.
Cross analysis is carried out as embodiment to the shale using Young's modulus * density, Poisson's ratio and quartz content below
The brittle evaluation method of reservoir is described in detail.
Shown in Fig. 2 is the figure that crosses of Young's modulus * density and Poisson's ratio, according to well logging information by Young's modulus * density
Attribute (E ρ), Poisson's ratio (Poisson Ratio) and quartz content (Vquartz) three's progress cross analysis, it may be determined that quartz contains
Amount is more than 60% (Vquartz>60%) boxed area in Fig. 2 is in region, Young's modulus * density corresponding to the boxed area
Codomain scope is more than 120GPa*g/cc.Value of the Young's modulus * density more than 120GPa*g/cc is projected to anti-by prestack
Drill on the Young's modulus * density attributes bodies of structure, be that Young's modulus * density characterizes shale reservoir fragility to be evaluated in shale storage
Response on layer section.Understood with reference to shown in Fig. 3, the higher region of shale reservoir fragility to be evaluated and log coincide compared with
It is good.Similarly, value of the Young's modulus * density more than 120GPa*g/cc is projected on Young's modulus * density attributes bodies, you can
Obtain the shale reservoir fragility planar distribution (quartz content to be evaluated shown in Fig. 4>60%).In Fig. 4, colour code represents Young mould
The property value of * density is measured, it is E* ρ that red arrives white to yellow again>120GPa*g/cc, corresponding fragility is more than 60% region, red
Color arrives white brittle again to yellow gradually to be strengthened, and white brittle is maximum.
The brittle evaluation method of shale reservoir proposed using present embodiment, pass through multiplying for Young's modulus and density
Product evaluation shale reservoir fragility, because the product of Young's modulus and density can be by the combined result of wave impedance in length and breadth come table
Sign, and wave impedance relative density in prestack inversion is more accurately, to thereby eliminate because density item accuracy is poor in length and breadth
The caused influence to shale reservoir brittleness evaluation precision, so as to improve brittleness evaluation precision, and it is crisp to whole research area
Mild-natured face prediction of spread, had broad application prospects in shale gas exploration and development field.
Present embodiment is that technical scheme is clearly and completely described, and embodiment therein is only
Only it is the part of the embodiment of the present invention, and is not whole embodiments.Based on the embodiment in the present invention, art technology
All other embodiment that personnel are obtained on the premise of not by creative work belongs to the protection model of the present invention
Enclose.
Claims (7)
- A kind of 1. brittle evaluation method of shale reservoir, it is characterised in that including:The elastic parameter and brittle mineral parameter of shale reservoir to be evaluated are obtained, the elastic parameter includes Young's modulus * density;Directly only the Young's modulus * density and other elastic parameters and the brittle mineral parameter cross point Analyse, and the codomain of the Young's modulus * density is determined according to the result of the cross analysis;Young's modulus * density attributes bodies are built according to the wave impedance in length and breadth obtained by prestack inversion;The Young's modulus * density attributes body is evaluated the fragility of the shale reservoir to be evaluated in the codomain;The Young's modulus * density is carried out into cross analysis with other elastic parameters and the brittle mineral parameter includes:Elastic parameter described in any two including Young's modulus * density and any one of brittle mineral parameter are handed over Converge and analyze, and the codomain of the Young's modulus * density is determined according to the result of the analysis.
- 2. the method as described in claim 1, it is characterised in that the elastic parameter also includes wave impedance, the extra large constant * of drawing in length and breadth Density and modulus of shearing * density, the brittle mineral parameter include quartz content, shale content and calcite content.
- 3. the method as described in claim 1, it is characterised in that the Young's modulus is calculated by below equation and obtained:E=3k (1- σ)Wherein, E represents Young's modulus, and k represents bulk modulus;σ represents Poisson's ratio.
- 4. method as claimed in claim 3, it is characterised in that the Poisson's ratio is calculated by below equation and obtained:<mrow> <mi>&sigma;</mi> <mo>=</mo> <mfrac> <mrow> <msup> <msub> <mi>V</mi> <mi>p</mi> </msub> <mn>2</mn> </msup> <mo>-</mo> <mn>2</mn> <msup> <msub> <mi>V</mi> <mi>s</mi> </msub> <mn>2</mn> </msup> </mrow> <mrow> <mn>2</mn> <msup> <msub> <mi>V</mi> <mi>p</mi> </msub> <mn>2</mn> </msup> <mo>-</mo> <mn>2</mn> <msup> <msub> <mi>V</mi> <mi>s</mi> </msub> <mn>2</mn> </msup> </mrow> </mfrac> </mrow>Wherein, VpRepresent velocity of longitudinal wave, VsRepresent shear wave velocity.
- 5. method as claimed in claim 4, it is characterised in that the Young's modulus * density attributes body passes through below equation meter Calculate and obtain:<mrow> <mi>E</mi> <mi>&rho;</mi> <mo>=</mo> <msup> <msub> <mi>I</mi> <mi>s</mi> </msub> <mn>2</mn> </msup> <mfrac> <mrow> <mn>3</mn> <msup> <msub> <mi>I</mi> <mi>p</mi> </msub> <mn>2</mn> </msup> <mo>-</mo> <mn>4</mn> <msup> <msub> <mi>I</mi> <mi>s</mi> </msub> <mn>2</mn> </msup> </mrow> <mrow> <msup> <msub> <mi>I</mi> <mi>p</mi> </msub> <mn>2</mn> </msup> <mo>-</mo> <msup> <msub> <mi>I</mi> <mi>s</mi> </msub> <mn>2</mn> </msup> </mrow> </mfrac> </mrow>Wherein, ρ represents density, IpRepresent p-wave impedance, IsRepresent S-wave impedance.
- 6. method as claimed in claim 2, it is characterised in that the brittle mineral parameter is according to the quartz content, shale Quantitative relationship between content and calcite content and the Young's modulus * density obtains.
- 7. the method as described in claim 1, it is characterised in that the elastic parameter and brittle mineral of the shale reservoir to be evaluated Parameter is obtained by the well-log information of the shale reservoir to be evaluated.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410790723.6A CN104564042B (en) | 2014-12-18 | 2014-12-18 | A kind of brittle evaluation method of shale reservoir |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410790723.6A CN104564042B (en) | 2014-12-18 | 2014-12-18 | A kind of brittle evaluation method of shale reservoir |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104564042A CN104564042A (en) | 2015-04-29 |
CN104564042B true CN104564042B (en) | 2017-12-01 |
Family
ID=53081027
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410790723.6A Active CN104564042B (en) | 2014-12-18 | 2014-12-18 | A kind of brittle evaluation method of shale reservoir |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104564042B (en) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104865124B (en) * | 2015-05-30 | 2017-10-27 | 重庆地质矿产研究院 | Shale brittleness index determination method based on rock stress-strain curve and ultrasonic longitudinal wave velocity |
CN105426650A (en) * | 2015-06-08 | 2016-03-23 | 中国石油集团川庆钻探工程有限公司 | Method for selecting out shale gas fracturing modification well section by using logging data |
CN105089661B (en) * | 2015-07-16 | 2017-12-05 | 中国石油天然气股份有限公司 | A kind of continuous method for judging speed stratum |
CN106569266B (en) * | 2015-10-12 | 2018-11-20 | 中国石油化工股份有限公司 | Shale gas reservoir brittle mineral content prediction method |
CN106568919A (en) * | 2015-10-13 | 2017-04-19 | 中国石油化工股份有限公司 | Shale brittleness prediction method based on rock physical analysis |
CN105370268B (en) * | 2015-10-23 | 2018-10-16 | 中国石油天然气集团公司 | A kind of method and apparatus of optimum level well staged fracturing parameter |
CN105257287A (en) * | 2015-11-05 | 2016-01-20 | 成都理工大学 | Method for evaluating parameters of oil reservoir of oil shale |
CN105865955A (en) * | 2016-03-23 | 2016-08-17 | 成都创源油气技术开发有限公司 | Logging evaluation method for brittleness of shale |
CN110552690A (en) * | 2018-05-30 | 2019-12-10 | 中国石油化工股份有限公司 | Shale reservoir brittleness evaluation method |
CN110186755A (en) * | 2019-04-23 | 2019-08-30 | 陕西国防工业职业技术学院 | One kind is with brill shale brittleness evaluation method |
CN110568150B (en) * | 2019-04-28 | 2022-03-01 | 中国石油天然气股份有限公司 | Oil shale identification method and device |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102967883A (en) * | 2012-11-20 | 2013-03-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for predicting rock brittleness probability through pre-stack elasticity parameter inversion of shale gas |
CN103122762A (en) * | 2013-01-30 | 2013-05-29 | 中国石油天然气股份有限公司 | Detection method and device for effective fractured intervals in unconventional shale oil and gas reservoir |
KR101415196B1 (en) * | 2013-10-15 | 2014-07-04 | 한국지질자원연구원 | Method for selecting fracking intervals of horizontal drilling zone in sweet spot range using resistivity and density logging data in shale play |
CN103982178A (en) * | 2014-04-16 | 2014-08-13 | 孙赞东 | Mineral content-based brittleness evaluation method for shale gas reservoir |
-
2014
- 2014-12-18 CN CN201410790723.6A patent/CN104564042B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102967883A (en) * | 2012-11-20 | 2013-03-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for predicting rock brittleness probability through pre-stack elasticity parameter inversion of shale gas |
CN103122762A (en) * | 2013-01-30 | 2013-05-29 | 中国石油天然气股份有限公司 | Detection method and device for effective fractured intervals in unconventional shale oil and gas reservoir |
KR101415196B1 (en) * | 2013-10-15 | 2014-07-04 | 한국지질자원연구원 | Method for selecting fracking intervals of horizontal drilling zone in sweet spot range using resistivity and density logging data in shale play |
CN103982178A (en) * | 2014-04-16 | 2014-08-13 | 孙赞东 | Mineral content-based brittleness evaluation method for shale gas reservoir |
Non-Patent Citations (1)
Title |
---|
泥页岩脆性地球物理预测技术;董宁等;《石油地球物理勘探》;20131231;第48卷(第S1期);第69-74页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104564042A (en) | 2015-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104564042B (en) | A kind of brittle evaluation method of shale reservoir | |
CN105182424B (en) | A kind of method and apparatus based on patchy saturation quantitative forecast reservoir porosity | |
CN105257286B (en) | A kind of method and device obtaining formation rock constituent content | |
CN104516017B (en) | A kind of carbonate rock petrophysical parameter seismic inversion method | |
CN104345345B (en) | A kind of shale reservoir total organic carbon TOC content prediction methods | |
CN104635269B (en) | A kind of method that fractured igneous-rock type reservoir is predicted based on prestack orientation angle gathers | |
CN104278991B (en) | Saline Lake Facies hydrocarbon source rock organic carbon and the polynary well logging computational methods of hydrocarbon potential | |
CN105675635B (en) | Tight rock component relative amount and brittleness index determine method and apparatus | |
CN104278989B (en) | Method for obtaining saturability index of low porosity and low permeability reservoir | |
CN103590827B (en) | Based on the compact clastic rock natural gas well PRODUCTION FORECASTING METHODS of Reservoir Classification | |
CN105986815A (en) | Method for recognizing geological sweet spots in shale stratum | |
CN102455436B (en) | Method for detecting anisotropic fracture of longitudinal noise attenuation prestack wave at limited azimuth angles | |
CN104005760B (en) | Based on the crack detection method of azimuthal anisotropy elastic impedance | |
CN109388817A (en) | A kind of Reservoir Fracture three-dimensional modeling method | |
CN104252007A (en) | Compatibility rock physics modeling method | |
CN104863574B (en) | A kind of Fluid Identification Method suitable for tight sandstone reservoir | |
CN103867198B (en) | Method for distinguishing formation density of carbonatite natural gas reservoir | |
CN108897066A (en) | Carbonate fracture density quantitative forecasting technique and device | |
CN105093313A (en) | Predicting method and apparatus for production capacity of single well in Karst oil-gas reservoir | |
CN110058323A (en) | A kind of tight sand formation brittleness index calculation method | |
CN105445441A (en) | Method and device for determining reservoir stratum gas saturation | |
CN106597543A (en) | Stratigraphic sedimentary facies division method | |
CN109425900A (en) | A kind of Seismic Reservoir Prediction method | |
CN104714252A (en) | Method for analyzing fluid factor sensibility | |
CN104727813A (en) | Porosity measuring method for resource storage commensalism type stratums |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |