CN105675635B - Tight rock component relative amount and brittleness index determine method and apparatus - Google Patents
Tight rock component relative amount and brittleness index determine method and apparatus Download PDFInfo
- Publication number
- CN105675635B CN105675635B CN201511031926.8A CN201511031926A CN105675635B CN 105675635 B CN105675635 B CN 105675635B CN 201511031926 A CN201511031926 A CN 201511031926A CN 105675635 B CN105675635 B CN 105675635B
- Authority
- CN
- China
- Prior art keywords
- relative amount
- indicatrix
- brittleness index
- rock
- mineral constituent
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/20—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/05—Investigating materials by wave or particle radiation by diffraction, scatter or reflection
- G01N2223/056—Investigating materials by wave or particle radiation by diffraction, scatter or reflection diffraction
Landscapes
- Chemical & Material Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides a kind of tight rock component relative amounts and brittleness index to determine method and apparatus, wherein, this method includes:It is analyzed based on X diffraction total rock, establishes each mineral constituent relative amount indicatrix in tight rock;Reconstruct each mineral constituent relative amount indicatrix and brittleness index indicatrix in tight rock;Each mineral constituent relative amount indicatrix and brittleness index indicatrix carry out seismic waveform instruction inverting, obtain each mineral constituent relative amount and the inverting data volume of brittleness index in tight rock as constraints using in the tight rock of reconstruct.The present invention solves each mineral constituent relative amount and brittleness index relation unobvious in compact reservoir in the prior art and causes be difficult on seismic profile the technical issues of, has reached the simple technique effect for accurately determining relation between each mineral constituent relative amount and brittleness index.
Description
Technical field
The present invention relates to technical field of geological exploration, more particularly to a kind of tight rock component relative amount and brittleness index
Determine method and apparatus.
Background technology
With the increasing of conventional oil and gas resource exploration development difficulty, unconventional petroleum resources will play important role.
In recent years, not only fine and close petroleum resources is found that in large-scale terrestrial lake basin, also have breakthrough in succession in small-sized fault depressed lacustrine basin.In view of
This, under the background of compact reservoir, the problem of prediction of " dessert " reservoir becomes extensive concern.Especially compact reservoir rock ore deposit
Object constituent content and brittleness index are both one of object that unconventionaloil pool exploration and development is paid close attention to and unconventionaloil pool
Explore the important parameter of " dessert " area selection.
The content of brittle mineral determines the effect of later stage fracturing reform in compact reservoir rock forming mineral component, and directly
Affect oil and gas production.Rock brittleness with high quartz content or high-carbon Carbonate Rocks content is relatively high, in fracturing process
Failure by shear can be generated and advantageously form complicated netted seam, by volume fracturing technology, well yield can be improved.In depth
Hollow band or sloped region have the rock of high potassium feldspar content, secondary pore growing, rich accumulation of oil and gas.
Therefore, how effectively to predict that brittle mineral constituent content and brittleness index spatial distribution are directly closed in compact reservoir
It is the exploration deployment work to next step.
At present, less to the research of rock forming mineral component and brittleness index, rock forming mineral component is mainly captured using element
Well logging (ECS) calculates mineral content, expensive, and the content of each mineral constituent can only can be just calculated in specific software.It is crisp
Sex index can indirectly be predicted using 2 elastic parameters of Young's modulus and Poisson's ratio, but cannot be quantified and be asked in rock respectively
Mineral constituent content can not disclose rock brittleness Producing reason.
How rock forming mineral component and brittleness index to be studied as a whole, not yet propose effective solution party at present
Case.
The content of the invention
An embodiment of the present invention provides a kind of tight rock component relative amounts and brittleness index to determine method, this method bag
It includes:
It is analyzed based on X diffraction total rock, establishes each mineral constituent relative amount indicatrix in tight rock;The fine and close rock of reconstruct
Each mineral constituent relative amount indicatrix and brittleness index indicatrix in stone;With each mineral constituent in the tight rock of reconstruct
Relative amount indicatrix and brittleness index indicatrix carry out seismic waveform instruction inverting, obtain densification as constraints
Each mineral constituent relative amount and the inverting data volume of brittleness index in rock.
In one embodiment, each mineral constituent relative amount indicatrix in tight rock is reconstructed, including:Pass through survey
Using multiple linear regression analysis method, it is bent to reconstruct each mineral constituent relative amount feature in tight rock for well curve sensitivity analysis
Line.
In one embodiment, each mineral constituent relative amount indicatrix is expressed as in the tight rock of reconstruct:
Y=a × AC+b × CN+c × DEN+d × GR+e × RD+f × RS+g
Wherein, AC represents interval transit time, and CN represents neutron, and DEN represents density, and GR represents natural gamma, and RD represents deep side
To RS represents shallow lateral, and a, b, c, d, e, f, g are the parameter after reconstruct, is constant.
In one embodiment, the brittleness index indicatrix in tight rock is reconstructed, including:According to the densification of reconstruct
Each mineral constituent relative amount indicatrix in rock reconstructs the brittleness index indicatrix in tight rock, wherein, the cause
The brittleness index indicatrix of close rock is expressed as:
Wherein, BI represents brittleness index, and Q represents the relative amount of quartz, and F represents the relative amount of feldspar, CARBONATE
Represent the relative amount of carbonate rock, CLAY represents the relative amount of clay.
In one embodiment, referred to each mineral constituent relative amount indicatrix and brittleness in the tight rock of reconstruct
Number indicatrix carries out seismic waveform instruction inverting, obtains each mineral constituent relative amount in tight rock as constraints
And the inverting data volume of brittleness index, including:With each mineral constituent relative amount indicatrix in the tight rock of reconstruct and crisp
Sex index indicatrix is as constraints, by analyzing prospect pit, determine the sample number of model foundation and high frequency into
Point;According to the sample number and radio-frequency component determined, initial model is established;Based on the initial model, pass through seismic wave
Shape indicates inverting, is finally inversed by the spatial distribution and brittleness index of the relative amount of each mineral constituent on well profile and even well profile
Spatial distribution;The spatial distribution for crossing the relative amount of each mineral constituent on well profile and even well profile that inverting is obtained
It is bent with each mineral constituent relative amount indicatrix and brittleness index feature of the reconstruct with the spatial distribution of brittleness index
Line is compared;If difference is less than predetermined error threshold, it is anti-that the initial model is applied to entire seismic data cube
Drill to obtain each mineral constituent relative amount and the inverting data volume of brittleness index in rock.
The embodiment of the present invention additionally provides a kind of tight rock component relative amount and brittleness index determining device, the device
Including:
Module is established, for being based on the analysis of X diffraction total rock, it is bent to establish each mineral constituent relative amount feature in tight rock
Line;Reconstructed module, for reconstructing each mineral constituent relative amount indicatrix and brittleness index indicatrix in tight rock;Instead
Drill module, for each mineral constituent relative amount indicatrix in the tight rock using reconstruct and brittleness index indicatrix as
Constraints carries out seismic waveform instruction inverting, obtains the anti-of each mineral constituent relative amount and brittleness index in tight rock
Drill data volume.
In one embodiment, the reconstructed module is specifically used for passing through logging sensitivity analy-sis, and utilization is polynary
Linear regression method reconstructs each mineral constituent relative amount indicatrix in tight rock.
In one embodiment, each mineral constituent relative amount feature is bent in the tight rock of the reconstructed module reconstruct
Line is expressed as:
Y=a × AC+b × CN+c × DEN+d × GR+e × RD+f × RS+g
Wherein, AC represents interval transit time, and CN represents neutron, and DEN represents density, and GR represents natural gamma, and RD represents deep side
To RS represents shallow lateral, and a, b, c, d, e, f, g are the parameter after reconstruct, is constant.
In one embodiment, the reconstructed module is specifically used for each mineral constituent phase in the tight rock according to reconstruct
To content characteristics curve, the brittleness index indicatrix in tight rock is reconstructed, wherein, the brittleness index of the tight rock is special
Sign curve is expressed as:
Wherein, BI represents brittleness index, and Q represents the relative amount of quartz, and F represents the relative amount of feldspar, CARBONATE
Represent the relative amount of carbonate rock, CLAY represents the relative amount of clay.
In one embodiment, the inverting module includes:Determination unit, for each ore deposit in the tight rock with reconstruct
Object component relative amount indicatrix and brittleness index indicatrix, by analyzing prospect pit, are determined as constraints
Go out the sample number and radio-frequency component of model foundation;Unit is established, for according to the sample number and radio-frequency component determined, establishing just
Beginning model;First inverting unit, for based on the initial model, indicating inverting by seismic waveform, being finally inversed by well
The spatial distribution of relative amount of each mineral constituent and the spatial distribution of brittleness index on section and even well profile;Comparison unit,
For the spatial distribution and brittleness for crossing the relative amount of each mineral constituent on well profile and even well profile for obtaining inverting
The spatial distribution of index, each mineral constituent relative amount indicatrix and brittleness index indicatrix with the reconstruct carry out
Comparison;Second inverting unit, in the case of being less than predetermined error threshold in definite difference, by the initial model application
Each mineral constituent relative amount and the inverting data volume of brittleness index in rock are obtained in entire seismic data cube inverting.
In embodiments of the present invention, it is special that each mineral constituent relative amount indicatrix and brittleness index in tight rock are reconstructed
Levy curve, and using in the tight rock of reconstruct each mineral constituent relative amount indicatrix and brittleness index indicatrix as about
Beam condition carries out seismic waveform instruction inverting, obtains the inverting of each mineral constituent relative amount and brittleness index in tight rock
Data volume solves each mineral constituent relative amount and brittleness index relation unobvious in compact reservoir in the prior art and causes
Be difficult on seismic profile the technical issues of, reach simple and accurately determined each mineral constituent relative amount and brittleness index
Between relation technique effect.
Description of the drawings
Attached drawing described herein is used for providing a further understanding of the present invention, forms the part of the application, not
Form limitation of the invention.In the accompanying drawings:
Fig. 1 is the method flow that tight rock constituent content according to embodiments of the present invention and brittleness index determine method
Figure;
Fig. 2 is each mineral group in utilization multiple linear regression analysis method reconstruct L1 well tight rocks according to embodiments of the present invention
Divide relative amount and brittleness index indicatrix schematic diagram;
Fig. 3 is that multiple linear regression analysis method according to embodiments of the present invention predicts mineral constituent content and measured value correlation
Analyze schematic diagram;
Fig. 4 is that quartz content (Q) inverting according to embodiments of the present invention connects well profile schematic diagram;
Fig. 5 is that feldspar content (F) inverting according to embodiments of the present invention connects well profile schematic diagram;
Fig. 6 is that brittleness index (BI) inverting according to embodiments of the present invention connects well profile schematic diagram;
Fig. 7 is the plane isopleth distribution schematic diagram of quartz content (Q) according to embodiments of the present invention;
Fig. 8 is the plane isopleth distribution schematic diagram of feldspar content (F) according to embodiments of the present invention;
Fig. 9 is the plane isopleth distribution schematic diagram of brittleness index (BI) according to embodiments of the present invention;
Figure 10 is tight rock constituent content according to embodiments of the present invention and the structure diagram of brittleness index determining device.
Specific embodiment
Understand to make the object, technical solutions and advantages of the present invention clearer, it is right with reference to embodiment and attached drawing
The present invention is described in further details.Here, the exemplary embodiment and its explanation of the present invention be for explaining the present invention, but simultaneously
It is not as a limitation of the invention.
In embodiments of the present invention, provide a kind of tight rock constituent content and brittleness index determines method, such as Fig. 1 institutes
Show, comprise the following steps:
Step 101:It is analyzed based on X diffraction total rock, establishes each mineral constituent relative amount indicatrix in tight rock;
Step 102:Reconstruct each mineral constituent relative amount indicatrix and brittleness index indicatrix in tight rock;
Specifically, using multiple linear regression analysis method, can be reconstructed by logging sensitivity analy-sis in tight rock
Each mineral constituent relative amount indicatrix, each mineral constituent relative amount indicatrix can represent in the tight rock of reconstruct
For:
Y=a × AC+b × CN+c × DEN+d × GR+e × RD+f × RS+g
Wherein, AC represents interval transit time, and CN represents neutron, and DEN represents density, and GR represents natural gamma, and RD represents deep side
To RS represents shallow lateral, and a, b, c, d, e, f, g are the parameter after reconstruct, is constant.
Correspondingly, can each mineral constituent relative amount indicatrix in the tight rock according to reconstruct, reconstruct fine and close rock
Brittleness index indicatrix in stone, wherein, the brittleness index indicatrix of tight rock can be expressed as:
Wherein, BI represents brittleness index, and Q represents the relative amount of quartz, and F represents the relative amount of feldspar, CARBONATE
Represent the relative amount of carbonate rock, CLAY represents the relative amount of clay.
When specific implementation, it can be analyzed by logging sensitivity analy-sis with reference to different methods, with
Reconstruct each mineral constituent relative amount indicatrix in tight rock.
Step 103:It is bent with each mineral constituent relative amount indicatrix in the tight rock of reconstruct and brittleness index feature
Line carries out seismic waveform instruction inverting, obtains each mineral constituent relative amount and brittleness in tight rock and refer to as constraints
Several inverting data volumes.
I.e., it is possible to each mineral constituent relative amount indicatrix and brittleness index indicatrix in the tight rock of reconstruct
As constraints, first to crossing well profile and even the instruction inverting of well profile progress seismic wave, each mineral group in tight rock is obtained
Divide relative amount and the inverting data volume of brittleness index.
Specifically, seismic waveform can be carried out by following steps and indicates inverting:
S1:Using in the tight rock of reconstruct each mineral constituent relative amount indicatrix and brittleness index indicatrix as
Constraints by analyzing prospect pit, determines the sample number and radio-frequency component of model foundation;
S2:According to the sample number and radio-frequency component determined, initial model is established;
S3:Based on the initial model, inverting is indicated by seismic waveform, was finally inversed by well profile and even well profile
The upper spatial distribution of relative amount of each mineral constituent and the spatial distribution of brittleness index;
S4:The spatial distribution for crossing the relative amount of each mineral constituent on well profile and even well profile that inverting is obtained
It is bent with each mineral constituent relative amount indicatrix and brittleness index feature of the reconstruct with the spatial distribution of brittleness index
Line is compared;
S5:If difference is less than predetermined error threshold, it is anti-that the initial model is applied to entire seismic data cube
Drill to obtain each mineral constituent relative amount and the inverting data volume of brittleness index in rock.
Each component content in above-mentioned tight rock and brittleness index are determined with reference to a specific embodiment method into
Row explanation, however, it is noted that, which is not formed to the present invention not merely to the present invention is better described
Work as restriction.
In this example primarily to solving how " dessert " storage is effectively predicted under compact reservoir background in the prior art
The problem of layer rock forming mineral component and brittleness index spatial distribution, is proposed one kind and is combined based on experiment, well logging, seismic facies therefore
Compact reservoir in rock forming mineral component and brittleness index spatial prediction techniques.With experimental analysis (analysis of X diffraction total rock) for base
Plinth, it is effective using seismic waveform instruction inversion technique different from conventional wave impedance inversion on Logging Curves analysis foundation
Rock forming mineral constituent content and brittleness index distribution and development degree in compact reservoir are predicted, has promoted unconventional reservoir
The development of middle rock forming mineral constituent content and brittleness index Predicting Technique.
Specifically, each component content and brittleness index determine that method mainly includes the following aspects in the tight rock:
1) brittleness index is calculated using each mineral constituent relative amount of experimental analysis in rock, can included:
S1:Rock core to typical prospect pit or typical outcrop samples systematic sampling in sampling process, can be set between sampling
0.5-1m is divided into, preferably has sample spot distribution in the different depth section of interval of interest, it so could utility well-log information
The indicatrix precision of forward modeling is higher;
S2:X diffraction total rock experimental analyses are carried out to the sample of acquisition;
S3:X diffraction total rock analysis results are classified, you can be divided into quartzy (Q), feldspar (F), carbonate rock
(CARBONATE), clay (CLAY) and other (OTHERS) five major class mineral, if conditions permit, can also segment feldspar
For potassium feldspar (KF) and plagioclase (XF) etc., the measured value for dividing these mineral and brittleness index (BI) is finely playbacked in lithology
On the premise of, it demarcates on well, for carrying out logging sensitivity analy-sis.
Specifically, the relative amount of each mineral constituent in tight rock can be utilized to obtain the expression of brittleness index (BI)
Formula:
Wherein, BI represents rock brittleness index, and V represents mineral volume, if V is well logging calculated value, then V is represented just
It is a volumn concentration, that is, represents the relative amount of each ore deposit component.
2) each mineral constituent relative amount and brittleness index in well logging quantitative judge tight rock, can include:
S1:Each mineral constituent relative amount and brittleness index log parameter sensitivity analysis in tight rock:
Each mineral constituent of rock core or outcrop samples in tight rock at different depth can be determined by above-mentioned analysis
Then relative amount and brittleness index size, can sum up different minerals component log response according to various conventional log
Feature.For example, it can sum up:Quartzy (Q) relative amount electrical property feature is more apparent, and resistivity value is higher, natural gamma value,
Middle subvalue and density value are relatively low.More than one high three low logging response character be quartzy relative amount main log response it is special
Sign;Clay mineral (CLAY) relative amount is more sensitive to natural gamma, density, interval transit time and resistivity, shows as nature
Gamma, density and interval transit time are higher, and resistivity is relatively low.The above are quartz and clay mineral relative amount logging response character,
Other mineral constituent relative amounts and brittleness index logging response character are different, the log type sensitive to its
Difference, so-called sensitive reference is exactly height of the log to the influence degree of each mineral constituent relative amount.
S2:Using multiple linear regression analysis method, each mineral constituent relative amount and brittleness index in quantitative forecast rock.
Since single factor test logging method multi-solution is larger, log parameter is multifactor concentrated expression.Therefore, can combine
Log parameter sensitivity analysis by analyzing distinct methods, preferably goes out each in multiple linear regression analysis method reconstruct tight rock
Mineral constituent relative amount and brittleness index indicatrix.For example, reconstruct quartzy (Q), feldspar (F), carbonate rock
(CARBONATE), clay (CLAY) and other (OTHERS) mineralogical composition relative amounts and the indicatrix of brittleness index can be with
It is expressed as:
Y=a × AC+b × CN+c × DEN+d × GR+e × RD+f × RS+g
Wherein, AC represents interval transit time, and CN represents neutron, and DEN represents density, and GR represents natural gamma, and RD represents deep side
To RS represents shallow lateral, and Q represents the relative amount of quartz, and F represents the relative amount of feldspar, and CARBONATE represents carbonate rock
The relative amount of (main including mineral such as calcite and dolomites), CLAY represent the relative amount of clay, a, b, c, d, e, f, g
It is constant for the parameter that multiple linear regression analysis method is fitted.
It will be compared by the result of each mineral constituent relative amount and experiment test in the tight rock for calculating of logging well,
And ensureing its correlation more than 0.7, correlation is higher, illustrates that fitting effect is preferable.
3) each mineral constituent relative amount and brittleness index in seismic inversion identification tight rock.
Using each mineral constituent relative amount in above-mentioned reconstruct tight rock and the indicatrix of brittleness index as constraints,
By analyzing prospect pit, preferably go out rational sample number and suitable radio-frequency component, rational initial model is established, with introductory die
Based on type, inverting is indicated with seismic waveform, first the company's of being finally inversed by well profile different minerals component relative amount and brittleness index
Spatial distribution, by inversion result and reconstruct indicatrix compare, if inversion result and match result on well are preferable,
Indicate that inverting finally obtains each mineral constituent relative amount and the inverting data volume of brittleness index in rock by seismic waveform.
Above-mentioned seismic waveform instruction inverting takes full advantage of the intensive seismic waveform information of spatial distribution in the horizontal, vertical
High-resolution well-log information is taken full advantage of upwards.Seismic waveform is densely distributed space structure data, and it is heavy to reflect
The spatial variations of product environment and lithology combination, seismic waveform instruction inverting utilize the preferably related well sample sheet of seismic waveform similitude,
Sample for reference spatial distribution distance and curve distribution feature establish initial model, instead of variogram analysis space variant structure,
Unbiased optimal estimation is carried out to radio-frequency component.
Specifically, which can include following three steps:
S1:Loading experiment analysis of data, well-log information and seismic data, and by the experimental analysis data in seismic data
It is demarcated on well-log information;
S2:It is preferred that related well sample sheet, sample for reference spatial distribution distance and curve distribution feature establish initial model;
S3:Preferably go out near well point with the most similar seismic waveform of indicatrix, using global optimization approach, promote
To no wellblock, inverting certainty greatly enhances, from completely random to progressively determining.By company's well profile of inverting and indicatrix value
Comparison is repeated, until establishing suitable geological model, geological model is applied to entire seismic data cube, it is anti-so as to reach
Drill the purpose of each mineral constituent relative amount and brittleness index data volume in tight rock.
Method is determined to each component content in above-mentioned tight rock and brittleness index with a specific experimental implementation process
It illustrates, comprises the following steps:
Step 1:Each mineral constituent relative amount and brittleness index in well logging quantitative judge tight rock.
First, log parameter sensitivity analysis is done to each mineral constituent relative amount and brittleness index in tight rock,
It is secondary by analyzing distinct methods, preferably go out multiple linear regression analysis method reconstruct tight rock in each mineral constituent relative amount and
Brittleness index indicatrix.As shown in Fig. 2, for using multiple linear regression analysis method reconstruct L1 well tight rock Mineralss and
Brittleness index indicatrix can sum up each mineral constituent relative amount logging response character in tight rock by the figure,
That is, quartzy (Q) relative amount electrical property feature is more apparent, and resistivity value is higher, natural gamma value, middle subvalue and density value compared with
It is low.More than one high three low logging response character be quartzy relative amount main logging response character;Clay mineral (CLAY)
Relative amount is more sensitive to natural gamma value, density value, interval transit time and resistivity, shows as natural gamma value, density value
Higher with interval transit time, resistivity value is relatively low.The above are quartz and clay mineral relative amount logging response character, other ore deposits
Object component relative amount and brittleness index logging response character are different.
Using the method for multiple linear regression, each mineral constituent relative amount in tight rock can be reconstructed, for example, stone
English relative amount (Q), feldspar relative amount (F), clay mineral relative amount (CLAY) and brittleness index (BI) and main well logging
The relation of response characteristic curve reconstructs indicatrix, and basis is provided for subsequent inverting work.
Quartz content (Q), feldspar content (F), clay mineral content (CLAY) and brittleness index (BI) feature of reconstruct are bent
Line is as follows:
Q=0.31*AC-1.23*CN-66.52*DEN+0.02*GR-0.04*RD+0.14*RS+136. 61;
F=-0.06*AC+1.44*CN-22.28*DEN-0.34*GR+0.002*RD+0.06*RS+97 .55;
CLAY=0.10*AC+3.28*CN+40.33*DEN-0.26*GR-0.08*RD+0.234*RS- 141.63;
It is analyzed by cross plot as shown in Figure 3, it can be with comparing calculation value (the fine and close rock calculated using log
Each mineral constituent relative amount in stone) with measured value (passing through X diffraction total rocks analysis result) correlation, correlation is higher, shows
Precision of prediction and effect are better.As shown in figure 3, quartz, feldspar, clay mineral and other mineral facies to content correlation all compared with
Height, substantially more than 0.75, in this way, reconstruct quartz and feldspar relative amount indicatrix calculated value are related to measured value
Property is more preferable, and related coefficient is up to more than 0.9.
Step 2:Each mineral constituent relative amount and brittleness index in seismic inversion identification rock:
Analyzed with reference to the logging character of step 1, with each mineral constituent relative amount in the tight rock that is reconstructed in step 1 and
The indicatrixes such as brittleness index are constraints, preferably go out rational sample number and suitable radio-frequency component, are established rational initial
Model crosses each mineral constituent relative amount and brittleness index on well profile and even well profile with seismic waveform instruction inversion prediction
Spatial distribution, wherein, spatial distribution refers to the geology distribution rule of each mineral constituent relative amount and brittleness index in the plane
Rule.If efficiency of inverse process is preferable with coincideing on well, indicate that inverting finally obtains each mineral group in tight rock by seismic waveform
Divide relative amount and the inverting data volume of brittleness index, wherein, inverting data volume is to have passed through seismic wave to original earthquake data
The data volume that shape inverting obtains, this data volume can have geologic rule that can also not have geologic rule.Divide as shown in Figs. 4-6
Not Wei in the inverting data volume quartzy relative amount (Q), feldspar relative amount (F), brittleness index (BI) it is anti-by seismic waveform
The company's well profile drilled, the quartzy relative amount (Q) of the reconstruct such as the log in Fig. 4 to 6 respectively Jing Guo squared processing,
The indicatrix of feldspar relative amount (F), brittleness index (BI).It can be seen that by company's well profile in Fig. 4 to 6 and pass through seismic wave
Quartzy relative amount (Q), feldspar relative amount (F) and brittleness index (BI) obtained by shape inverting and statistical result (well column on well
The signature logging curve on sub- side) it is consistent substantially, this explanation is with each in multiple linear regression analysis method reconstruct tight rock
Mineral constituent relative amount and brittleness index indicatrix are indicated inverting by seismic waveform, fine and close rock can be gone out with Accurate Prediction
Each mineral constituent relative amount and brittleness index distribution characteristics and development degree in stone.
However, due to each mineral constituent relative amount in rock and brittleness index and wave impedance relation unobvious, in addition
It is difficult on seismic profile, it is opposite that each mineral constituent in rock cannot be preferably predicted on conventional wave impedance inversion data volume
Content and brittleness index distribution characteristics and development degree, and use each mineral group in multiple linear regression analysis method reconstruct tight rock
Divide relative amount and brittleness index indicatrix, inverting is indicated by waveform, each mineral in tight rock can be gone out with Accurate Prediction
Component relative amount and brittleness index distribution characteristics and development degree, this is also the advantage place of the application.
According to each mineral constituent relative amount and brittleness index in the rock as shown in Fig. 7 to 9 on inversion result
Plane distribution situation and development degree carry out " dessert " reservoir prediction under compact reservoir background system evaluation, and brittleness are produced
Raw reason has carried out rational geologic interpretation, and the exploration portion of fine and close oil gas and shale oil gas is instructed by these achievements of acquirement
Administration, wherein, Fig. 7 show the plane isopleth distribution schematic diagram of quartz content (Q), and Fig. 8 show the plane of feldspar content (F)
Distribution of contours schematic diagram, Fig. 9 show the plane isopleth distribution schematic diagram of brittleness index (BI).Because by quartz content or
Based on carbonate rock content, the Brittleness of generation can generate failure by shear in fracturing process and advantageously form complexity
Netted seam by volume fracturing technology, can improve well yield.It is special by feldspar, particularly brittleness caused by potassium feldspar content
Sign, in deep hollow band or sloped region, secondary pore growing, secondary pore growing zone is lower step exploration targets.
In upper example, experimental analysis technology (analysis of X diffraction total rock), characteristic curve construction technology and earthquake mainly make use of
Waveform inversion technology, based on X diffraction total rock analysiss of data, using Logging Curves, with the method weight of multiple linear regression
The indicatrixes such as each mineral constituent content and brittleness index in structure tight rock, and using the indicatrix of reconstruct as constraint item
Part predicts rock Mineralss relative amount and brittleness index spatial distribution characteristic with seismic waveform instruction inversion technique
And development degree.
Mainly it is considered that brittle mineral relative amount and brittleness index and wave impedance relation unobvious in compact reservoir,
Be difficult on seismic profile, although wave impedance inversion and more attribution inversion part well point prediction results are preferable, it is phased not
Substantially, no wellblock precision can not meet exploration needs.In this example, with seismic waveform inversion technology, sample for reference space point
It is more demanding to the selection of sample number and the selection of radio-frequency component when cloth distance and curve distribution feature establish initial model.Tool
Body, it is analyzed based on X diffraction total rock, by log analysis, chooses the well logging to each mineral constituent content sensitivity in tight rock
Curve reconstructs indicatrix with the method for multiple linear regression;Each mineral constituent relative amount and brittleness index in rock
Prediction needs to carry out sensitivity curves analysis on the basis of well log interpretation, logging explanation and formation testing conclusion, and then selection is suitable
Sample number and radio-frequency component establish rational initial model, with each in seismic waveform instruction inversion technique prediction compact reservoir
Mineral constituent relative amount and brittleness index spatial feature and development degree.
It is unknown that above-described embodiment solves brittle mineral relative amount and brittleness index and wave impedance relation in compact reservoir
Aobvious, the difficult point being difficult on seismic profile is combined by experiment-well logging-seismic inversion, Accurate Prediction brittle mineral
Relative amount and brittleness index distribution and development degree, the skill being combined breakthroughly with well log interpretation with seismic inversion
Art effectively predicts brittle mineral relative amount and brittleness index distribution characteristics in compact reservoir.The Predicting Technique has obtained well
Application effect, play production actual effect well, pass through the brittle mineral relative amount and brittleness index of the technological prediction point
Cloth feature and the goodness of fit on well are higher, and poor oil reservoir, oil reservoir feldspar content are high, and brittleness index is high.Especially in commercial oil Jing Chu,
This feature becomes apparent, and illustrates that the main reason for brittleness generates in the example is feldspar content.
Based on same inventive concept, a kind of tight rock component relative amount and brittleness are additionally provided in the embodiment of the present invention
Index determining device, as described in the following examples.Due to tight rock component relative amount and brittleness index determining device solution
Certainly the principle of problem determines that method is similar to tight rock component relative amount and brittleness index, therefore tight rock component is opposite
The implementation of content and brittleness index determining device may refer to tight rock component relative amount and brittleness index determines method
Implement, overlaps will not be repeated.Used below, term " unit " or " module " can realize the software of predetermined function
And/or the combination of hardware.Although following embodiment described device is preferably realized with software, hardware or soft
The realization of the combination of part and hardware is also what may and be contemplated.Figure 10, which is that the tight rock component of the embodiment of the present invention is opposite, to be contained
A kind of structure diagram of amount and brittleness index determining device, as shown in Figure 10, including:Establish module 1001, reconstructed module 1002
With inverting module 1003, the structure is illustrated below.
Module 1001 is established, for being based on the analysis of X diffraction total rock, it is special to establish each mineral constituent relative amount in tight rock
Levy curve;
Reconstructed module 1002, it is special for reconstructing each mineral constituent relative amount indicatrix and brittleness index in tight rock
Levy curve;
Inverting module 1003 refers to for each mineral constituent relative amount indicatrix and brittleness in the tight rock with reconstruct
Number indicatrix carries out seismic waveform instruction inverting, obtains each mineral constituent relative amount in tight rock as constraints
And the inverting data volume of brittleness index.
In one embodiment, reconstructed module 1002 is specific can be used for, by logging sensitivity analy-sis, utilizing
Multiple linear regression analysis method reconstructs each mineral constituent relative amount indicatrix in tight rock.
In one embodiment, each mineral constituent relative amount feature is bent in the tight rock that reconstructed module 1002 reconstructs
Line can be expressed as:
Y=a × AC+b × CN+c × DEN+d × GR+e × RD+f × RS+g
Wherein, AC represents interval transit time, and CN represents neutron, and DEN represents density, and GR represents natural gamma, and RD represents deep side
To RS represents shallow lateral, and a, b, c, d, e, f, g are the parameter after reconstruct, is constant.
In one embodiment, reconstructed module 1002 specifically can be used for each mineral group in the tight rock according to reconstruct
Split-phase reconstructs the brittleness index indicatrix in tight rock to content characteristics curve, wherein, the brittleness of the tight rock refers to
Number indicatrix is expressed as:
Wherein, BI represents brittleness index, and Q represents the relative amount of quartz, and F represents the relative amount of feldspar, CARBONATE
Represent the relative amount of carbonate rock, CLAY represents the relative amount of clay.
In one embodiment, inverting module 1003 can include:Determination unit, in the tight rock with reconstruct
Each mineral constituent relative amount indicatrix and brittleness index indicatrix are as constraints, by analyzing prospect pit,
Determine the sample number and radio-frequency component of model foundation;Unit is established, for according to the sample number and radio-frequency component determined, building
Vertical initial model;First inverting unit, for based on the initial model, indicating inverting by seismic waveform, being finally inversed by
Cross the spatial distribution of relative amount of each mineral constituent and the spatial distribution of brittleness index on well profile and even well profile;Comparison is single
Member described crosses well profile and the even spatial distribution of the relative amount of each mineral constituent and crisp on well profile for obtain inverting
The spatial distribution of sex index, each mineral constituent relative amount indicatrix and brittleness index indicatrix with the reconstruct, into
Row comparison;Second inverting unit, should by the initial model in the case of being less than predetermined error threshold in definite difference
Each mineral constituent relative amount and the inverting data volume of brittleness index in rock are obtained for entire seismic data cube inverting.
It can be seen from the above description that the embodiment of the present invention realizes following technique effect:It reconstructs in tight rock
Each mineral constituent relative amount indicatrix and brittleness index indicatrix, and with each mineral constituent phase in the tight rock of reconstruct
To content characteristics curve and brittleness index indicatrix as constraints, seismic waveform instruction inverting is carried out, obtains fine and close rock
Each mineral constituent relative amount and the inverting data volume of brittleness index, solve each mineral in compact reservoir in the prior art in stone
Component relative amount and brittleness index relation unobvious and cause be difficult on seismic profile the technical issues of, to have reached letter
Single technique effect for accurately determining relation between each mineral constituent relative amount and brittleness index.
Obviously, those skilled in the art should be understood that each module of the above-mentioned embodiment of the present invention or each step can be with
It is realized with general computing device, they can concentrate on single computing device or be distributed in multiple computing devices
On the network formed, optionally, they can be realized with the program code that computing device can perform, it is thus possible to by it
Store and performed in the storage device by computing device, and in some cases, can be held with the order being different from herein
They are either fabricated to each integrated circuit modules or will be multiple in them by the shown or described step of row respectively
Module or step are fabricated to single integrated circuit module to realize.In this way, the embodiment of the present invention be not restricted to it is any specific hard
Part and software combine.
The foregoing is only a preferred embodiment of the present invention, is not intended to limit the invention, for the skill of this field
For art personnel, the embodiment of the present invention can have various modifications and variations.Within the spirit and principles of the invention, made
Any modification, equivalent substitution, improvement and etc. should all be included in the protection scope of the present invention.
Claims (6)
1. a kind of tight rock component relative amount and brittleness index determine method, which is characterized in that including:
It is analyzed based on X diffraction total rock, establishes each mineral constituent relative amount indicatrix in tight rock;
Reconstruct each mineral constituent relative amount indicatrix and brittleness index indicatrix in tight rock;The densification of the reconstruct
Each mineral constituent relative amount indicatrix is expressed as in rock:
Y=a × AC+b × CN+c × DEN+d × GR+e × RD+f × RS+g
Wherein, AC represents interval transit time, and CN represents neutron, and DEN represents density, and GR represents natural gamma, and RD represents deep lateral, RS
Represent shallow lateral, a, b, c, d, e, f, g are the parameter after reconstruct, are constant;
Each mineral constituent relative amount indicatrix and brittleness index indicatrix are as constraint item using in the tight rock of reconstruct
Part carries out seismic waveform instruction inverting, obtains the inverting data of each mineral constituent relative amount and brittleness index in tight rock
Body, including:Each mineral constituent relative amount indicatrix and brittleness index indicatrix are as about using in the tight rock of reconstruct
Beam condition by analyzing prospect pit, determines the sample number and radio-frequency component of model foundation;According to the sample number determined
And radio-frequency component, establish initial model;Based on the initial model, inverting is indicated by seismic waveform, was finally inversed by well
The spatial distribution of relative amount of each mineral constituent and the spatial distribution of brittleness index on section and even well profile;Inverting is obtained
It is described cross well profile and the even spatial distribution of relative amount of each mineral constituent and the spatial distribution of brittleness index on well profile,
With each mineral constituent relative amount indicatrix and brittleness index indicatrix of the reconstruct, compared;If difference is small
In predetermined error threshold, then the initial model is obtained into each mineral constituent in rock applied to entire seismic data cube inverting
The inverting data volume of relative amount and brittleness index.
2. the method as described in claim 1, which is characterized in that each mineral constituent relative amount feature is bent in reconstruct tight rock
Line, including:
By logging sensitivity analy-sis, using multiple linear regression analysis method, it is opposite to reconstruct each mineral constituent in tight rock
Content characteristics curve.
3. the method as described in claim 1, which is characterized in that the brittleness index indicatrix in reconstruct tight rock, including:
According to each mineral constituent relative amount indicatrix in the tight rock of reconstruct, the brittleness index reconstructed in tight rock is special
Curve is levied, wherein, the brittleness index indicatrix of the tight rock is expressed as:
<mrow>
<mi>B</mi>
<mi>I</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>Q</mi>
<mo>+</mo>
<mi>F</mi>
<mo>+</mo>
<mi>C</mi>
<mi>A</mi>
<mi>R</mi>
<mi>B</mi>
<mi>O</mi>
<mi>N</mi>
<mi>A</mi>
<mi>T</mi>
<mi>E</mi>
</mrow>
<mrow>
<mi>Q</mi>
<mo>+</mo>
<mi>F</mi>
<mo>+</mo>
<mi>C</mi>
<mi>A</mi>
<mi>R</mi>
<mi>B</mi>
<mi>O</mi>
<mi>N</mi>
<mi>A</mi>
<mi>T</mi>
<mi>E</mi>
<mo>+</mo>
<mi>C</mi>
<mi>L</mi>
<mi>A</mi>
<mi>Y</mi>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
</mrow>
Wherein, BI represents brittleness index, and Q represents the relative amount of quartz, and F represents the relative amount of feldspar, and CARBONATE is represented
The relative amount of carbonate rock, CLAY represent the relative amount of clay.
4. a kind of tight rock component relative amount and brittleness index determining device, which is characterized in that including:
Module is established, for being based on the analysis of X diffraction total rock, establishes each mineral constituent relative amount indicatrix in tight rock;
Reconstructed module, for reconstructing each mineral constituent relative amount indicatrix and brittleness index indicatrix in tight rock;
Each mineral constituent relative amount indicatrix is expressed as in the tight rock of the reconstructed module reconstruct:
Y=a × AC+b × CN+c × DEN+d × GR+e × RD+f × RS+g
Wherein, AC represents interval transit time, and CN represents neutron, and DEN represents density, and GR represents natural gamma, and RD represents deep lateral, RS
Represent shallow lateral, a, b, c, d, e, f, g are the parameter after reconstruct, are constant;
Inverting module, it is bent for each mineral constituent relative amount indicatrix in the tight rock with reconstruct and brittleness index feature
Line carries out seismic waveform instruction inverting, obtains each mineral constituent relative amount and brittleness in tight rock and refer to as constraints
Several inverting data volumes;Wherein, the inverting module includes:
Determination unit, it is bent for each mineral constituent relative amount indicatrix in the tight rock with reconstruct and brittleness index feature
Line, by analyzing prospect pit, determines the sample number and radio-frequency component of model foundation as constraints;
Unit is established, for according to the sample number and radio-frequency component determined, establishing initial model;
First inverting unit, for based on the initial model, indicating inverting by seismic waveform, being finally inversed by well profile
And connect the spatial distribution of relative amount of each mineral constituent and the spatial distribution of brittleness index on well profile;
Comparison unit, for the sky for crossing the relative amount of each mineral constituent on well profile and even well profile for obtaining inverting
Between the spatial distribution of distribution and brittleness index, it is special with each mineral constituent relative amount indicatrix and brittleness index of the reconstruct
Curve is levied, is compared;
Second inverting unit, in the case of being less than predetermined error threshold in definite difference, by the initial model application
Each mineral constituent relative amount and the inverting data volume of brittleness index in rock are obtained in entire seismic data cube inverting.
5. device as claimed in claim 4, which is characterized in that the reconstructed module is specifically used for through log sensibility
Analysis using multiple linear regression analysis method, reconstructs each mineral constituent relative amount indicatrix in tight rock.
6. device as claimed in claim 4, which is characterized in that the reconstructed module is specifically used for the tight rock according to reconstruct
In each mineral constituent relative amount indicatrix, reconstruct tight rock in brittleness index indicatrix, wherein, it is described densification rock
The brittleness index indicatrix of stone is expressed as:
<mrow>
<mi>B</mi>
<mi>I</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>Q</mi>
<mo>+</mo>
<mi>F</mi>
<mo>+</mo>
<mi>C</mi>
<mi>A</mi>
<mi>R</mi>
<mi>B</mi>
<mi>O</mi>
<mi>N</mi>
<mi>A</mi>
<mi>T</mi>
<mi>E</mi>
</mrow>
<mrow>
<mi>Q</mi>
<mo>+</mo>
<mi>F</mi>
<mo>+</mo>
<mi>C</mi>
<mi>A</mi>
<mi>R</mi>
<mi>B</mi>
<mi>O</mi>
<mi>N</mi>
<mi>A</mi>
<mi>T</mi>
<mi>E</mi>
<mo>+</mo>
<mi>C</mi>
<mi>L</mi>
<mi>A</mi>
<mi>Y</mi>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
</mrow>
Wherein, BI represents brittleness index, and Q represents the relative amount of quartz, and F represents the relative amount of feldspar, and CARBONATE is represented
The relative amount of carbonate rock, CLAY represent the relative amount of clay.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511031926.8A CN105675635B (en) | 2015-12-31 | 2015-12-31 | Tight rock component relative amount and brittleness index determine method and apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511031926.8A CN105675635B (en) | 2015-12-31 | 2015-12-31 | Tight rock component relative amount and brittleness index determine method and apparatus |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105675635A CN105675635A (en) | 2016-06-15 |
CN105675635B true CN105675635B (en) | 2018-06-01 |
Family
ID=56298489
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511031926.8A Active CN105675635B (en) | 2015-12-31 | 2015-12-31 | Tight rock component relative amount and brittleness index determine method and apparatus |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105675635B (en) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106248494B (en) * | 2016-08-29 | 2019-03-01 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | A method of for shale gas well reservoir brittleness overall merit |
CN106547034B (en) * | 2016-11-09 | 2018-11-02 | 西南石油大学 | A kind of method of determining compact reservoir rock brittleness index |
CN108519281B (en) * | 2018-02-27 | 2020-08-11 | 中国石油天然气股份有限公司 | Method, device and system for determining brittleness index of rock |
CN110320571B (en) * | 2018-03-29 | 2021-06-15 | 中国石油化工股份有限公司 | Compact sandstone reservoir rock brittleness logging evaluation method |
CN108825223B (en) * | 2018-06-22 | 2021-10-22 | 西南石油大学 | Method for extracting bedding characteristics of shale stratum |
CN108827774B (en) * | 2018-06-23 | 2019-06-07 | 东北石油大学 | Coal petrography brittleness evaluation method |
CN110094198A (en) * | 2019-06-05 | 2019-08-06 | 中国地质大学(北京) | Oil well grey matter content assaying method based on SPSS software |
CN110513098B (en) * | 2019-07-15 | 2022-07-01 | 中国海洋石油集团有限公司 | Method for improving longitudinal resolution of element logging curve |
CN113125413A (en) * | 2019-12-31 | 2021-07-16 | 中石化石油工程技术服务有限公司 | Shale brittleness online evaluation method based on laser-induced breakdown spectroscopy |
CN111429012A (en) * | 2020-03-27 | 2020-07-17 | 中国石油天然气股份有限公司大港油田分公司 | Shale brittle dessert evaluation method |
CN112505154B (en) * | 2020-11-11 | 2021-06-29 | 中国地质大学(北京) | Shale reservoir mineral component content analysis and lithofacies identification characterization method |
CN114791482B (en) * | 2021-01-25 | 2024-03-08 | 中国石油化工股份有限公司 | Method and device for determining rock mineral content coefficient and storage medium |
CN113030440B (en) * | 2021-03-22 | 2022-02-15 | 西南石油大学 | Method for predicting fracture brittleness index of shale in hydraulic fracturing process of shale |
-
2015
- 2015-12-31 CN CN201511031926.8A patent/CN105675635B/en active Active
Non-Patent Citations (3)
Title |
---|
南方古生界3套富有机质页岩矿物组成与脆性分析;陈吉,肖贤明;《煤炭学报》;20130531;第38卷(第5期);第824页右栏第2-6段 * |
基于岩石物理实验的致密油储层脆性指数预测—以柴西南A井区为例;张平 等;《中国石油学会2015年物探技术研讨会论文集》;20150513;第772页右栏第2段-第775页左栏第2段 * |
致密地层岩石脆性指数的测井优化建模;范卓颖 等;《石油学报》;20151130;第36卷(第11期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105675635A (en) | 2016-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105675635B (en) | Tight rock component relative amount and brittleness index determine method and apparatus | |
CN107817535B (en) | The determination method and apparatus of short lap | |
Grana | Joint facies and reservoir properties inversion | |
CN101930082B (en) | Method for distinguishing reservoir fluid type by adopting resistivity data | |
Verma et al. | Estimation of total organic carbon and brittleness volume | |
CN105510993A (en) | Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method | |
CN102012526A (en) | Method for discriminating type of reservoir fluid by using resistivity data | |
Wu et al. | Seismic sedimentology, facies analyses, and high-quality reservoir predictions in fan deltas: a case study of the Triassic Baikouquan Formation on the western slope of the Mahu Sag in China's Junggar Basin | |
CN105116449B (en) | A kind of recognition methods of weak reflection reservoir | |
CN105301647B (en) | The method for distinguishing grey matter mud stone and sandstone | |
CN106249300B (en) | It shakes to combine based on well and determines source rock TOC contents and anisotropism method and apparatus | |
CN102096098A (en) | Seismic attribute extracting method | |
Qian et al. | Prediction and modeling of petrophysical parameters of deep-buried, low permeability glutenite reservoirs in Yubei area, Turpan-Hami Basin, China | |
Leary et al. | Power-law scaling of spatially correlated porosity and log (permeability) sequences from north-central North Sea Brae oilfield well core | |
Richardson | Well correlation and Petrophysical analysis, a case study of “Rickie” field onshore Niger Delta | |
Ryberg et al. | Tomographic V p and V s structure of the California Central Coast Ranges, in the vicinity of SAFOD, from controlled-source seismic data | |
Khan et al. | Three-dimensional Structural Modeling (3D SM) and Joint Geophysical Characterization (JGC) of Hydrocarbon Reservoir: A Case Study of the Kadanwari field in Middle Indus Basin (MIB), Southeastern Pakistan | |
Gao et al. | Seismic structure and texture analyses for fractured reservoir characterization: An integrated workflow | |
Wu et al. | Combination of seismic attributes using clustering and neural networks to identify environments with sandstone-type uranium mineralization | |
CN112147676A (en) | Method for predicting thickness of coal bed and gangue | |
Hu et al. | Seismic interpretation of sandstone-type uranium deposits in the Songliao Basin, Northeast China | |
CN110174701A (en) | A kind of prediction technique and device of Igneous rock invasion range | |
Yuan et al. | Interpretation of fractures and joint inversion using multicomponent seismic data—Marcellus Shale example | |
Sun et al. | 3D spatial characterization of sand body for uranium reservoirs based on geostatistical resistivity inversion | |
Dong et al. | An introduction—Carbonate geophysics |
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 |