CN105204066B - A kind of direct indicating means in coal seam Igneous rock invasion position based on spectral factorization - Google Patents
A kind of direct indicating means in coal seam Igneous rock invasion position based on spectral factorization Download PDFInfo
- Publication number
- CN105204066B CN105204066B CN201510700647.XA CN201510700647A CN105204066B CN 105204066 B CN105204066 B CN 105204066B CN 201510700647 A CN201510700647 A CN 201510700647A CN 105204066 B CN105204066 B CN 105204066B
- Authority
- CN
- China
- Prior art keywords
- coal seam
- frequency
- igneous rock
- cross plot
- intercept
- 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.)
- Expired - Fee Related
Links
Abstract
A kind of direct indicating means in coal seam Igneous rock invasion position based on spectral factorization, belongs to coal seam igneous rock method for indicating position.Method is:1) by man-machine interaction, coal seam reflection line-ups is picked up on conventional three-dimensional seismic data cube;2) coal seam reflected waveform data is extracted along coal seam reflection line-ups, generates subdata body;3) utilize the seismic channel in spectral factorization algorithm, subdata body to carry out spectral factorization, obtain corresponding time-frequency spectrum;4) for each time-frequency spectrum, the maximum response of its each frequency curve is searched for, maximum response corresponding time and frequency are depicted as cross plot;5) scatterplot in linear regression analysis method, fitting cross plot is utilized, intercept and Grad is obtained;6) intercept and Grad of all seismic channels are depicted as cross plot, indicate the intrusion position of igneous rock.Prediction is true, reliable, Forecasting Methodology science, simple and direct, predicts the outcome quantitative, directly perceived, high-precision, meets the requirement of the highly efficient and productive back production in mine district coal seam.
Description
Technical field
The present invention relates to a kind of coal seam igneous rock method for indicating position, particularly a kind of coal seam igneous rock based on spectral factorization
Invade the direct indicating means in position.
Background technology
In many coal mining enterprises of China, there is more serious Igneous rock invasion in its main mining coal seam.Except part coal
Beyond the igneous rock that layer is broken into is substituted completely, the igneous rock that most of coal seam is only partly broken into is substituted.When intrusion fire into
When petrosa point substitutes coal seam, its position in coal seam mainly has three kinds, i.e., at the top of, bottom and middle part.
Due to igneous rock hardness be much larger than coal seam hardness, during down-hole coal excavation, the plant equipment such as cutting coal machine be easy to because
Igneous rock is cut and accelerated wear test, while coal cutting efficiency is also reduced, so as to cause larger economic loss.In coal seam back production
Before, if accurate location of the igneous rock that can be invaded with Accurate Prediction in coal seam.In the coal seam back production design phase, it is possible to pre-
First consider the influence of coal seam Igneous rock invasion factor.Designed by optimizing coal seam back production, it is to avoid occur coal-cutting machinery because of cutting fire
Diagenesis and the mechanical accelerated wear test caused and production efficiency reduction.
For the prediction of coal seam Igneous rock invasion position, current most popular method is borehole standardization method.I.e.
By the coal seam Igneous rock invasion position disclosed by existing drilling in exploiting field, the coal seam Igneous rock invasion position in exploiting field is predicted.
Because coal seam is for its roof and floor country rock and igneous rock, than relatively soft.Igneous rock is in invaded coal layer, and it invades position
May because of locus difference, occur more violent change.Therefore, the igneous rock position of borehole standardization invaded coal layer is passed through
Put, precision is relatively low, reliability is poor.
On the other hand, igneous rock can typically be invaded in invaded coal layer along coal seam.Therefore, in invaded coal layer fire into
Rock in layered distribution, forms sill as coal seam.Carry out coalfield 3-d seismic exploration when, though coal seam whether by fire into
Rock is invaded, and typically can all have preferable coal seam back wave.It is different from normal coal seam back wave, when the fire containing intrusion in coal seam
During diagenesis, corresponding change can all occur for the kinematics and dynamic characteristic of coal seam back wave.Moreover, kinematics and dynamics are special
The change levied, can have certain difference because of the difference of Igneous rock invasion position.
The content of the invention
The invention aims to provide it is a kind of have prediction true, reliable, it is Forecasting Methodology science, simple and direct based on spectrum point
The direct indicating means in coal seam Igneous rock invasion position of solution.
The object of the present invention is achieved like this:The direct indicating means in igneous rock position:1) by man-machine interaction, in routine
Coal seam reflection line-ups is picked up on 3-d seismic data set;2) according to the feature of coal seam back wave, along coal seam reflected wave in phase
Axle extracts coal seam reflected waveform data, generates subdata body;3) each seismic channel in spectral factorization algorithm, subdata body is utilized
Spectral factorization is carried out, corresponding time-frequency spectrum is obtained;4) for each time-frequency spectrum, the peak response of its each frequency curve is searched for
Value, cross plot is depicted as by maximum response corresponding time and frequency;5) linear regression analysis method is utilized, cross plot is fitted
In scatterplot, obtain intercept and Grad;6) intercept and Grad of all seismic channels are depicted as cross plot, according to each
Position of the seismic channel in cross plot, directly indicates the intrusion position of igneous rock.
The direct indicating means in igneous rock position, is comprised the following steps that:
1) by man-machine interaction, coal seam reflection line-ups is picked up on conventional three-dimensional seismic data cube;
Specific method is as follows:1. existing well-log information in exploiting field is collected, synthetic seismogram is made;2. earthquake will be synthesized
Record is contrasted with well lie, recognizes position and the phase of coal output layer back wave;3. the coal seam back wave position that will identify that
With phase as control point, coal seam reflection line-ups is picked up in the whole district using man-machine interaction;
2) according to the feature of coal seam back wave, coal seam reflected waveform data is extracted along coal seam reflection line-ups, subnumber is generated
According to body;
Specific method is as follows:1. according to the feature of coal seam back wave, determine data extract when the time window upper bound, time window
Lower bound and time window length;2. the time window function of extraction is defined as bell-shaped Gaussian window;3. reflected along the coal seam of interaction pickup
Ripple lineups, extract the coal seam reflected waveform data of each seismic channel successively, generate subdata body;
3) utilize each seismic channel in spectral factorization algorithm, subdata body to carry out spectral factorization, obtain corresponding time-frequency
Spectrum;
Specific method is as follows:1. extract in subdata body per coal seam reflected waveform data together;2. according to coal seam back wave
Feature, the initial frequency of reasonable selection time-frequency spectrum, terminate the parameter such as frequency and step-length;3. suitable spectral factorization algorithm is chosen,
For wavelet transformation or S-transformation etc.;4. the time-frequency spectrum of each seismic channel in data volume is calculated;
4) for each time-frequency spectrum, the maximum response of its each frequency curve is searched for, maximum response is corresponding
Time and frequency are depicted as cross plot;
Specific method is as follows:1. the maximum response of each frequency curve in time-frequency spectrum is searched for;2. each frequency is read bent
The line maximum response corresponding time;3. by abscissa of frequency, the maximum response corresponding time be ordinate, draw intersection
Figure;
5) scatterplot in linear regression analysis method, fitting cross plot is utilized, intercept and Grad is calculated;
Specific method is as follows:1. according to the characteristics of scatterplot is distributed in cross plot, the data area of linear regression, starting are chosen
Frequency elects 20Hz as, terminates frequency and elects 100Hz as;2. according to principle of least square method, line is entered to the scatterplot in the range of data
Property return;3. the linear equation obtained according to linear regression, the intercept and Grad of digital simulation straight line;
6) intercept and Grad of all seismic channels are depicted as cross plot, according to each seismic channel in cross plot
Position, directly indicates the intrusion position of igneous rock;
Specific method is as follows:1. the intercept and gradient corresponding to all seismic channels are sorted successively;2. intercept and ladder are drawn
The cross plot of degree;3. according to the characteristic distributions of scatterplot in cross plot, threshold line is defined, it is preliminary to indicate Igneous rock invasion position;④
It will tentatively indicate that result is contrasted with known drilling exposure situation, and verify whether the threshold line of definition will be reasonable;If 5. rationally,
The then preliminary direct instruction for indicating that igneous rock position is coal seam Igneous rock invasion position;If unreasonable, threshold value is redefined
Line, repeat step is 4. untill obtaining satisfactory result.
Beneficial effect, as a result of such scheme, by extracting the kinematics and dynamic characteristic of coal seam back wave, in advance
The intrusion position of coal seam igneous rock is surveyed, spectral factorization is carried out to coal seam back wave, the time-frequency spectrum of coal seam back wave is obtained;Recycle line
Property regression analysis, quantitatively obtain variation relation of the coal seam back wave response maximum corresponding time with wavelet frequency;Finally, lead to
Setting threshold line is crossed, the intrusion position of igneous rock in coal seam is directly indicated;The 3D seismic data space lattice of mine district is close
Degree improves the precision of prediction of coal seam Igneous rock invasion position much larger than the space lattice density of drilling.
1) prediction achievement is more directly perceived, and the intrusion position with coal seam igneous rock is directly corresponding.The prediction achievement of the present invention is direct
The Igneous rock invasion position in correspondence coal seam, that is, it is from top, bottom or middle part invaded coal layer to correspond to igneous rock in coal seam.Therefore,
It is used directly for instructing the design of coal seam stoping scheme, advantageously reduces the accelerated wear test of the plant equipment such as cutting coal machine, also have
Beneficial to the raising of coal seam back production efficiency.
2) precision of detection of coal seam Igneous rock invasion position is higher.Because the present invention uses 3-D seismics in prediction
Data volume, the space sampling densities of data are significantly larger than drill hole density.Therefore, ratio is portrayed for coal seam Igneous rock invasion position
It is more accurate, hence it is evident that to reduce the horizontal swing for predicting coal seam Igneous rock invasion position, predict coal seam Igneous rock invasion position
Horizontal swing is less than 20m.
3) reliability of detection achievement is higher.Due to initial data of the present invention be 3D seismic data, without
It is drilling data, the space sampling densities of data reach 5m × 5m, even more high.Therefore, with higher reliability and prediction
Accuracy rate.
4) Forecasting Methodology is simple, and production efficiency is higher.
Advantage:The precision of prediction of coal seam Igneous rock invasion position is improved, makes mine district coal seam Igneous rock invasion position
Quantitative forecast be possibly realized, it is true, reliable with prediction, it is Forecasting Methodology science, simple and direct, predict the outcome quantitative, directly perceived, high-precision
The features such as spending, fully meets the requirement of the highly efficient and productive back production in mine district coal seam.
Brief description of the drawings:
Fig. 1 picks up schematic diagram for the coal seam back wave layer position interaction of the present invention.
Fig. 2 is coal seam back wave seismic channel correspondence time-frequency spectrum schematic diagram of the invention.
The frequency-time cross plot schematic diagram that Fig. 3 is extracted in the time-frequency spectrum for the present invention.
Fig. 4 is linear regression rear cut-off distance-gradient cross plot schematic diagram of the invention.
Fig. 5 directly indicates coal seam Igneous rock invasion position flow chart for the intercept-gradient cross plot of the present invention.
In figure:1st, seismic channel;2nd, the layer position of interaction pickup;3rd, the time window upper bound;4th, time window lower bound;5th, time-frequency spectrum is responded
Curve;6th, the maximum of points of time-frequency spectrum response curve;7th, the scatterplot in frequency-time cross plot;8th, linear regression fit straight line;
Scatterplot when the 9th, being invaded at the top of igneous rock in intercept-gradient cross plot;Intercept-gradient cross plot when the 10th, being invaded in the middle part of igneous rock
In scatterplot;11st, scatterplot when igneous rock bottom is invaded in intercept-gradient cross plot;12nd, intrusion and middle part are invaded at the top of igneous rock
The threshold line entered;13rd, the threshold line of the intrusion of igneous rock bottom and middle part intrusion.
Embodiment
The invention will be further described below in conjunction with the accompanying drawings
Embodiment 1:A kind of direct indicating means in coal seam Igneous rock invasion position based on spectral factorization comprises the following steps:
1) by man-machine interaction, coal seam reflection line-ups is picked up on conventional three-dimensional seismic data cube;2) according to coal seam
The feature of back wave, coal seam reflected waveform data is extracted along coal seam reflection line-ups, generates subdata body;3) calculated using spectral factorization
Each seismic channel in method, subdata body carries out spectral factorization, obtains corresponding time-frequency spectrum;4) for each time-frequency spectrum,
The maximum response of its each frequency curve is searched for, maximum response corresponding time and frequency are depicted as cross plot;5) it is sharp
The scatterplot in linear regression analysis method, fitting cross plot is used, intercept and Grad is obtained;6) by the intercept of all seismic channels and
Grad is depicted as cross plot, according to position of each seismic channel in cross plot, directly indicates the intrusion position of igneous rock.
The direct indicating means in igneous rock position, is comprised the following steps that:
1) by man-machine interaction, coal seam reflection line-ups is picked up on conventional three-dimensional seismic data cube.Specific method is such as
Under:
1. existing well-log information in exploiting field is collected, synthetic seismogram is made;
2. synthetic seismogram is contrasted with well lie, recognizes position and the phase of coal output layer back wave;
3. the coal seam back wave position that will identify that and phase are as control point, using man-machine interaction in whole district's pickup as schemed
Coal seam reflection line-ups shown in 1.
2) according to the feature of coal seam back wave, coal seam reflected waveform data is extracted along coal seam reflection line-ups, subnumber is generated
According to body.Specific method is as follows:
1. according to the feature of coal seam back wave, the time window upper bound, time window lower bound and the time window when data are extracted are determined
The lower bound of time window and the difference in the upper bound are time window length in length (be usually 11ms, 21ms, 31ms or 41ms etc.), such as Fig. 1;
2. the time window function of extraction is defined as bell-shaped Gaussian window, its function expression is h (t);
3. along the coal seam reflection line-ups of interaction pickup, the coal seam reflected waveform data of each seismic channel is extracted successively, is pressed
Calculated according to formula one, generate subdata body.
subdata(t-t0)=data (t) * h (t-t0) formula one
Wherein, t0For the time window upper bound as shown in Figure 1, t is the sampling point corresponding time in seismic channel, when data (t) is t
The initial data volume data at quarter, subdata (t) is t subdata volume data.
3) utilize each seismic channel in spectral factorization algorithm, subdata body to carry out spectral factorization, obtain corresponding time-frequency
Spectrum.Specific method is as follows:
1. extract in subdata body per coal seam reflected waveform data subdata (t) together;
2. according to the feature of coal seam back wave, the initial frequency f of reasonable selection time-frequency spectrum0, terminate frequency fendAnd step delta
The parameters such as f;
3. suitable spectral factorization algorithm, such as wavelet transformation or S-transformation are chosen;
4. formula two is utilized, the time-frequency spectrum S (t, f) of each seismic channel in data volume is calculated.
Formula two
Wherein the t times, f is frequency, Ω (τ, generating function when f) for spectral factorization.
4) for each time-frequency spectrum, the maximum response of its each frequency curve is searched for, maximum response is corresponding
Time and frequency are depicted as cross plot.Specific method is as follows:
1. the maximum response of each frequency curve in time-frequency spectrum is searched for, as shown in Figure 2;
2. each frequency curve maximum response corresponding time is read;
3. by abscissa of frequency, the maximum response corresponding time be ordinate, be depicted as intersection as shown in Figure 3
Figure.
5) scatterplot in linear regression analysis method, fitting cross plot is utilized, intercept and Grad is calculated.Specific method is such as
Under:
1. according to the characteristics of scatterplot is distributed in cross plot, the data area of linear regression is chosen, general initial frequency is elected as
20Hz, terminates frequency and elects 100Hz as;
2. according to principle of least square method, linear regression is carried out to the scatterplot in the range of data;
3. the linear equation obtained according to linear regression, the intercept and Grad of digital simulation straight line, such as the institute of formula three
Show.
T=Gf+P formula three
Wherein, t is the time, and f is frequency, and G is gradient, and P is intercept.
6) intercept and Grad of all seismic channels are depicted as cross plot, according to each seismic channel in cross plot
Position, directly indicates the intrusion position of igneous rock, and specific method is as follows:
1. the intercept P and gradient G corresponding to all seismic channels are sorted successively;
2. it is depicted as intercept-gradient cross plot as shown in Figure 4;
3. according to the characteristic distributions of scatterplot in cross plot, threshold line as shown in Figure 4 is defined, Igneous rock invasion is tentatively indicated
Position;
4. it will tentatively indicate that result is contrasted with known drilling exposure situation, and verify whether the threshold line of definition will be reasonable;
If 5. rationally, the preliminary direct instruction for indicating that igneous rock position is coal seam Igneous rock invasion position;If no
Rationally, then threshold line is redefined, repeat step is 4. untill obtaining satisfactory result.Implement process as shown in Figure 5.
Claims (1)
1. a kind of direct indicating means in coal seam Igneous rock invasion position based on spectral factorization, it is characterized in that:
1)By man-machine interaction, coal seam reflection line-ups is picked up on conventional three-dimensional seismic data cube;
Specific method is as follows:1. existing well-log information in exploiting field is collected, synthetic seismogram is made;2. by synthetic seismogram
Contrasted with well lie, recognize position and the phase of coal output layer back wave;3. the coal seam back wave position that will identify that and phase
Coal seam reflection line-ups is picked up using man-machine interaction in the whole district in position as control point;
2)According to the feature of coal seam back wave, coal seam reflected waveform data is extracted along coal seam reflection line-ups, subdata body is generated;
Specific method is as follows:1. according to the feature of coal seam back wave, determine data extract when the time window upper bound, time window lower bound
And time window length;2. the time window function of extraction is defined as bell-shaped Gaussian window;3. the coal seam back wave along interaction pickup is same
Phase axle, extracts the coal seam reflected waveform data of each seismic channel successively, generates subdata body;
3)Using spectral factorization algorithm, each seismic channel in subdata body carries out spectral factorization, obtains corresponding time-frequency spectrum;
Specific method is as follows:1. extract in subdata body per coal seam reflected waveform data together;2. according to the spy of coal seam back wave
Levy, initial frequency, termination frequency and the step parameter of reasonable selection time-frequency spectrum;3. suitable spectral factorization algorithm is chosen, is small echo
Conversion or S-transformation;4. the time-frequency spectrum of each seismic channel in data volume is calculated;
4)For each time-frequency spectrum, the maximum response of its each frequency curve is searched for, by the maximum response corresponding time
Cross plot is depicted as with frequency;
Specific method is as follows:1. the maximum response of each frequency curve in time-frequency spectrum is searched for;2. each frequency curve is read most
The big response corresponding time;3. by abscissa of frequency, the maximum response corresponding time be ordinate, draw cross plot;
5)Using linear regression analysis method, the scatterplot in cross plot is fitted, intercept and Grad is calculated;
Specific method is as follows:1. according to the characteristics of scatterplot is distributed in cross plot, the data area of linear regression, initial frequency are chosen
Elect 20Hz as, terminate frequency and elect 100Hz as;2. according to principle of least square method, the scatterplot in the range of data is linearly returned
Return;3. the linear equation obtained according to linear regression, the intercept and Grad of digital simulation straight line;
6)The intercept and Grad of all seismic channels are depicted as cross plot, according to position of each seismic channel in cross plot
Put, directly indicate the intrusion position of igneous rock;
Specific method is as follows:1. the intercept and gradient corresponding to all seismic channels are sorted successively;2. intercept and gradient are drawn
Cross plot;3. according to the characteristic distributions of scatterplot in cross plot, threshold line is defined, it is preliminary to indicate Igneous rock invasion position;4. will be just
Step indicates that result is contrasted with known drilling exposure situation, verifies whether the threshold line of definition is reasonable;If 5. rationally, just
Step indicates the direct instruction that igneous rock position is coal seam Igneous rock invasion position;If unreasonable, threshold line is redefined, weight
Multiple step is 4. untill obtaining satisfactory result.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510700647.XA CN105204066B (en) | 2015-10-26 | 2015-10-26 | A kind of direct indicating means in coal seam Igneous rock invasion position based on spectral factorization |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510700647.XA CN105204066B (en) | 2015-10-26 | 2015-10-26 | A kind of direct indicating means in coal seam Igneous rock invasion position based on spectral factorization |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105204066A CN105204066A (en) | 2015-12-30 |
CN105204066B true CN105204066B (en) | 2017-07-21 |
Family
ID=54951841
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510700647.XA Expired - Fee Related CN105204066B (en) | 2015-10-26 | 2015-10-26 | A kind of direct indicating means in coal seam Igneous rock invasion position based on spectral factorization |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105204066B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110208860B (en) * | 2019-07-02 | 2021-03-26 | 中国煤炭地质总局地球物理勘探研究院 | Method and device for predicting intrusion range of igneous rock |
CN113050193B (en) * | 2019-12-27 | 2023-08-22 | 中国石油天然气股份有限公司 | Distribution identification method, device and system for basement igneous rock |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5671136A (en) * | 1995-12-11 | 1997-09-23 | Willhoit, Jr.; Louis E. | Process for seismic imaging measurement and evaluation of three-dimensional subterranean common-impedance objects |
US5831935A (en) * | 1996-03-05 | 1998-11-03 | Chevron U.S.A. Inc. | Method for geophysical processing and interpretation using seismic trace difference for analysis and display |
CN102707317A (en) * | 2010-10-27 | 2012-10-03 | 中国石油化工股份有限公司 | Method of using absorption and attenuation characteristics of seismic wave for reservoir analysis |
CN103412332A (en) * | 2013-01-22 | 2013-11-27 | 中国地质大学(北京) | Method for quantitative calculation of thickness of thin reservoir layer |
CN103527184A (en) * | 2013-10-28 | 2014-01-22 | 北京大学 | Method and system for predicting dolomite reservoir |
-
2015
- 2015-10-26 CN CN201510700647.XA patent/CN105204066B/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5671136A (en) * | 1995-12-11 | 1997-09-23 | Willhoit, Jr.; Louis E. | Process for seismic imaging measurement and evaluation of three-dimensional subterranean common-impedance objects |
US5831935A (en) * | 1996-03-05 | 1998-11-03 | Chevron U.S.A. Inc. | Method for geophysical processing and interpretation using seismic trace difference for analysis and display |
US6317384B1 (en) * | 1996-03-05 | 2001-11-13 | Chevron U.S.A., Inc. | Method for geophysical processing and interpretation using seismic trace difference for analysis and display |
CN102707317A (en) * | 2010-10-27 | 2012-10-03 | 中国石油化工股份有限公司 | Method of using absorption and attenuation characteristics of seismic wave for reservoir analysis |
CN103412332A (en) * | 2013-01-22 | 2013-11-27 | 中国地质大学(北京) | Method for quantitative calculation of thickness of thin reservoir layer |
CN103527184A (en) * | 2013-10-28 | 2014-01-22 | 北京大学 | Method and system for predicting dolomite reservoir |
Non-Patent Citations (3)
Title |
---|
AVO属性分析方法的改进及其在煤系地层中的应用;舒梦珵 等;《地球物理学进展》;20100830;第25卷(第4期);第1323-1331页 * |
煤层岩浆岩侵入区的交会图定量预测技术—以卧龙湖煤矿为例;陈同俊 等;《煤炭学报》;20121231;第37卷(第12期);第2070-2075页 * |
谱分解技术在预测煤层岩浆岩侵入区的应用;陈同俊 等;《中国煤炭地质》;20100830;第22卷(第8期);第16-17、57页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105204066A (en) | 2015-12-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103454686B (en) | Small scale sedimentary facies based on strata slicing carries out the method and system of reservoir prediction | |
CN111257926B (en) | Method for predicting ancient valley uranium reservoir by using old seismic data | |
CN104749617B (en) | A kind of multi-scale facture reservoir forward method for establishing model | |
CN103899285B (en) | Multilamellar sandstone oil reservoir near limit Water-cut Period is rotated water drive method | |
CN105527653B (en) | A kind of virtual log construction method based on geological information | |
CN104280770B (en) | Prediction method of compact transition rock reservoir stratum | |
CN103382838A (en) | Reservoir stratum analysis method and device based on pressing-ability of fracturing geological body | |
CN104502969A (en) | Channel sandstone reservoir identification method | |
CN106597543B (en) | A kind of Sedimentary Facies division methods | |
RU2610454C2 (en) | Method of formation fracturing intensification location points determination using mineralogical composition, as well as corresponding system and software | |
CN104747144B (en) | Geological stratum selection method based on natural gas multi-stratum series partial pressure combined solution | |
CN111090709A (en) | Big data geological analysis method for sandstone-type uranium ore mineralization prediction | |
CN108442882B (en) | Shale gas large-displacement horizontal well while-drilling geosteering method based on two-dimensional seismic data | |
KR101415196B1 (en) | Method for selecting fracking intervals of horizontal drilling zone in sweet spot range using resistivity and density logging data in shale play | |
CN107132573A (en) | A kind of method that application wavelet decomposition reconfiguration technique recognizes the lower lithological pool of strong impedance shielding | |
CN103065051A (en) | Method for performing grading and sectionalizing on rock mass automatically | |
CN105242307A (en) | Complex carbonate stratum earthquake porosity obtaining method and apparatus | |
KR101415197B1 (en) | Method for selecting fracking intervals of horizontal drilling zone in sweet spot range using resistivity and neutron logging data in shale play | |
CN111161403A (en) | Method for acquiring top and bottom surfaces of sedimentary stratum based on spatial interpolation of drilling data | |
CN105204066B (en) | A kind of direct indicating means in coal seam Igneous rock invasion position based on spectral factorization | |
CN108663714A (en) | A kind of sedimentary micro depicting method | |
CN105301647A (en) | Method for distinguishing calcareous mudstone from sandstone | |
CN106353806B (en) | It is a kind of to detect the connective method of reef reservoir | |
CN106285642A (en) | A kind of formation pore pressure Forecasting Methodology based on seismic data | |
CN106597545A (en) | Seismic prestack inversion method and device for horizontal cracks |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170721 Termination date: 20201026 |