CN106886045A - A kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving - Google Patents
A kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving Download PDFInfo
- Publication number
- CN106886045A CN106886045A CN201510935169.0A CN201510935169A CN106886045A CN 106886045 A CN106886045 A CN 106886045A CN 201510935169 A CN201510935169 A CN 201510935169A CN 106886045 A CN106886045 A CN 106886045A
- Authority
- CN
- China
- Prior art keywords
- carried out
- distribution
- deformation coal
- coal
- out based
- 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.)
- Pending
Links
- 239000003245 coal Substances 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000004088 simulation Methods 0.000 claims abstract description 7
- 239000011435 rock Substances 0.000 claims abstract description 6
- 239000012530 fluid Substances 0.000 claims description 4
- 238000001615 p wave Methods 0.000 claims description 3
- 208000035126 Facies Diseases 0.000 claims description 2
- 235000009754 Vitis X bourquina Nutrition 0.000 claims 1
- 235000012333 Vitis X labruscana Nutrition 0.000 claims 1
- 240000006365 Vitis vinifera Species 0.000 claims 1
- 235000014787 Vitis vinifera Nutrition 0.000 claims 1
- 238000013213 extrapolation Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000003292 diminished effect Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6242—Elastic parameters, e.g. Young, Lamé or Poisson
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides a kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving, including：Crossed the step of obtaining deformation coal Rock Elastic Parameters and change using forward simulation；Offset distance domain CRP road collections are transformed into incident angle domain and the step of branch is superimposed；Extract wavelet respectively from each partial stack data volume and the step of interpolation goes out corresponding wavelet at any incident angle；The step of building initial model；The step of back analysis being carried out at known well point and asks for inverted parameters；Prestack Simultaneous Retrieving obtains Rock Elastic Parameters and predicts the step of deformation coal is distributed with this.
Description
Technical field
The invention belongs to coal seismic prospecting field, and in particular to a kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving.
Background technology
With the continuous intensification of coal mining depth, requirement to coalfield 3-d seismic exploration technology also more and more higher.In terms of the prediction of deformation coal distributed area, the wave impedance information that the inverting of poststack P-wave impedance is provided can not well distinguish primary coal and deformation coal, can not meet the requirement of deformation coal prediction.The technological synthesis of prestack Simultaneous Retrieving carry out Simultaneous Inversion using prestack incidence angle geological data, have obtained velocity of longitudinal wave, shear wave velocity and density this 3 basic elastic parameters.From this 3 basic parameters, more elastic parameters can be obtained in turn（Poisson's ratio, Lame Coefficient, Young's modulus, modulus of shearing, volume-based model etc.）, comprehensively utilize these elastic parameters carry out stratum lithology prediction and reservoir containing fluid properties detection, not only increase precision of prediction, and reduce the multi-solution of inversion result, enhance the ability for portraying reservoir properties.
The content of the invention
The present invention, to the not enough present situation of deformation coal distributed area forecasting research, according to the difference of deformation coal elastic parameter, proposes a kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving for geophysics means.
To reach above-mentioned purpose, technical scheme includes：
Drawn using forward simulation cross analysis, deformation coal Rock Elastic Parameters show as λ ρ and the μ equal step-downs of ρ values；
Drawn using forward simulation cross analysis, diminished than the Zp values of deformation coal with primary coal facies, Vp/Vs values become big；
CRP road collections are transformed into by incident angle domain by offset distance domain using stack velocity；
Angle gathers are divided into two or three parts by incidence angle distance to be overlapped, partial stack data volume is obtained；
Wavelet is extracted respectively from each partial stack data volume, and corresponding wavelet at any incident angle is gone out according to the wavelet interpolation for obtaining；
Build initial model, including P-wave And S surge impedance model and density model；
Back analysis is carried out at known well point, is that inverted parameters solution is i.e. right k、kc 、mWithmc Solution；
Velocity of longitudinal wave, shear wave velocity and the density obtained by prestack Simultaneous Retrieving, some other elastic parameter can be calculated, more especially parameter more sensitive to lithology and more sensitive to fluid, the parameter such as such as p-wave impedance, in length and breadth ripple ratio, λ ρ and μ ρ, the distribution of deformation coal is predicted according to the change of these elastic parameters.
The beneficial effects of the present invention are the rock physicses elastic parameter drawn there is provided a kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving, comprehensive utilization prestack Simultaneous Retrieving（Poisson's ratio, Lame Coefficient, Young's modulus, modulus of shearing, volume-based model etc.）, it is capable of the distribution of accurate prediction deformation coal.
Brief description of the drawings
Fig. 1 show a kind of flow chart of the embodiment of the method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving of the present invention.
Fig. 2 show forward simulation λ ρ-μ ρ and crosses figure.
Fig. 3 show forward simulation Vp/Vs-Zp- and crosses figure.
Specific embodiment
Specific implementation process of the present invention is described in detail below in conjunction with specific accompanying drawing.
Fig. 1 show a kind of flow chart of the embodiment of the method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving of the present invention.
Below in conjunction with a kind of embodiment of the method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving of the present invention shown in Fig. 1, is discussed in detail, it includes：
Step 101：Primary coal and deformation coal forward model are built, being crossed by λ ρ-μ ρ is drawn and cause λ ρ and μ ρ values less and less by the change of primary coal to deformation coal, therefore with the increase of degree of crushing, λ ρ and the μ equal step-downs of ρ values；
Step 102：Primary coal and deformation coal forward model are built, is intersected by Vp/Vs-Zp and is shown that deformation coal has smaller Zp and bigger Vp/Vs values；
Step 103：CRP road collections are gone into incident angle domain by offset distance domain using stacking velocity field；
Step 104：Angle gathers data by incidence angle distance are divided into two or three ranges of incidence angles it are overlapped and obtains partial stack body；
Step 105：Wavelet is extracted respectively from each partial stack body, and corresponding wavelet at any incident angle is then gone out according to the wavelet interpolation for obtaining
Step 106：Initial model, including P-wave And S surge impedance model and density model are built by well-log information and layer position data；
Step 107：At known point（At well point）Back analysis is carried out, suitable inverted parameters are determined k、kc 、mWithmc After carry out whole district's inverting；
Step 108：Velocity of longitudinal wave, shear wave velocity and density are obtained by prestack Simultaneous Retrieving, and then some other elastic parameter can also be obtained by mathematical operation, more especially parameter more sensitive to lithology and more sensitive to fluid, such as ripple ratio, λ ρ and μ ρ parameters in length and breadth, the distribution of deformation coal is predicted according to the change of these elastic parameters.
Claims (8)
1. a kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving, it is characterised in that including：
The step of discrimination model of deformation coal being obtained using forward simulation cross analysis；
The step of CRP road collections in offset distance domain are transformed into incident angle domain；
The step of angle gathers data be divided into two or three scopes being overlapped；
Extract wavelet respectively from each partial stack body and the step of interpolation goes out corresponding wavelet at any incident angle；
The step of building initial model；
The step of back analysis obtains inverted parameters；
The step of whole district's inverting obtains rock physicses elastic parameter and predicts deformation coal distribution with this.
2. a kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving according to claim 1, it is characterized in that, the discrimination model that the AVO forward simulations cross analysis are obtained with deformation coal crosses and Vp/Vs-Zp crosses including λ ρ-μ ρ, cause λ ρ and μ ρ values less and less by the change of primary coal to deformation coal, diminish than the Zp values of deformation coal and Vp/Vs values become big with primary coal facies.
3. a kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving according to claim 1, it is characterised in that the CRP road collections in offset distance domain are transformed into incident angle domain use stack velocity to described.
4. a kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving according to claim 1, it is characterized in that, angle gathers data are divided into two or three parts to described, be by the distance of incidence angle be divided into remote, nearly two angle gathers or it is remote, in, nearly three angle gathers parts, then the data of each part be overlapped obtain superposition of data body respectively.
5. a kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving according to claim 1, it is characterized in that, wavelet is extracted respectively from each partial stack body to described, corresponding wavelet is obtained by institute's grape ripple interpolation at any incident angle degree.
6. a kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving according to claim 1, it is characterized in that, the step of to the structure initial model, includes P-wave And S impedance model and density model, based on well-log information and well data, with reference to seismic horizon information, interpolation and extrapolation are carried out to log data using Kriging regression method, smooth, closure a physical model is produced（Such as surge impedance model and density model）.
7. a kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving according to claim 1, it is characterized in that, to it is described inverted parameters are obtained by back analysis the step of be that back analysis experiment is carried out at known well point, determine suitable inverted parameters i.e. rightk、kc 、mWithmc Solution, it is after taking natural logrithm to p-wave impedance, S-wave impedance and density respectively, then to carry out intersecting what is obtained respectively that it asks for process.
8. a kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving according to claim 1, it is characterized in that, the rock physicses elastic parameter obtained to whole district's inverting includes velocity of longitudinal wave, shear wave velocity and density, further some other elastic parameter can also be gone out by mathematical computations, more especially parameter more sensitive to lithology and more sensitive to fluid, the parameter such as such as p-wave impedance, in length and breadth ripple ratio, λ ρ and μ ρ, the distribution of deformation coal is predicted according to the change of these parameters.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510935169.0A CN106886045A (en) | 2015-12-15 | 2015-12-15 | A kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510935169.0A CN106886045A (en) | 2015-12-15 | 2015-12-15 | A kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106886045A true CN106886045A (en) | 2017-06-23 |
Family
ID=59174502
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510935169.0A Pending CN106886045A (en) | 2015-12-15 | 2015-12-15 | A kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106886045A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108761531A (en) * | 2018-04-12 | 2018-11-06 | 中国矿业大学 | A kind of Coal Pore Structure intelligent Forecasting based on SVM theories |
| CN109239777A (en) * | 2018-08-22 | 2019-01-18 | 中国地质大学（北京） | A method of it is developed using joint inversion method detection deformation coal |
| CN110632663A (en) * | 2019-09-27 | 2019-12-31 | 岭南师范学院 | Longitudinal and transverse wave velocity ratio inversion method without well control |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103163553A (en) * | 2013-02-01 | 2013-06-19 | 中国石油天然气股份有限公司 | Earthquake hydrocarbon detection method and detection device based on multiple pore medium model |
| CN104850897A (en) * | 2015-02-25 | 2015-08-19 | 中国矿业大学 | Prediction method for coal and gas outburst based on seismic information |
-
2015
- 2015-12-15 CN CN201510935169.0A patent/CN106886045A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103163553A (en) * | 2013-02-01 | 2013-06-19 | 中国石油天然气股份有限公司 | Earthquake hydrocarbon detection method and detection device based on multiple pore medium model |
| CN104850897A (en) * | 2015-02-25 | 2015-08-19 | 中国矿业大学 | Prediction method for coal and gas outburst based on seismic information |
Non-Patent Citations (3)
| Title |
|---|
| 孙学凯 等: "联合弹性波阻抗反演与同步反演确定构造煤的分布", 《煤炭学报》 * |
| 李娟娟: "煤与瓦斯突出预测的岩性地震反演方法研究", 《中国博士学位论文全文数据库 工程科技I辑》 * |
| 王远 等: "利用地震反演信息划分煤体结构", 《煤田地质与勘探》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108761531A (en) * | 2018-04-12 | 2018-11-06 | 中国矿业大学 | A kind of Coal Pore Structure intelligent Forecasting based on SVM theories |
| CN108761531B (en) * | 2018-04-12 | 2019-06-21 | 中国矿业大学 | A kind of Coal Pore Structure intelligent Forecasting based on SVM theory |
| CN109239777A (en) * | 2018-08-22 | 2019-01-18 | 中国地质大学（北京） | A method of it is developed using joint inversion method detection deformation coal |
| CN110632663A (en) * | 2019-09-27 | 2019-12-31 | 岭南师范学院 | Longitudinal and transverse wave velocity ratio inversion method without well control |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dix | Seismic velocities from surface measurements | |
| Patacca et al. | Structural architecture of the central Apennines: Interpretation of the CROP 11 seismic profile from the Adriatic coast to the orographic divide | |
| McGilvery et al. | The influence of local gradients on accommodation space and linked depositional elements across a stepped slope profile, offshore Brunei | |
| Collot et al. | Origin of a crustal splay fault and its relation to the seismogenic zone and underplating at the erosional north Ecuador–south Colombia oceanic margin | |
| Deveugle et al. | Characterization of stratigraphic architecture and its impact on fluid flow in a fluvial-dominated deltaic reservoir analog: Upper Cretaceous Ferron Sandstone Member, Utah | |
| US20140129479A1 (en) | Method to aid in the exploration, mine design, evaluation and/or extraction of metalliferous mineral and/or diamond deposits | |
| Osinowo et al. | Reservoir description and characterization of Eni field Offshore Niger Delta, southern Nigeria | |
| Trippanera et al. | Fault and graben growth along active magmatic divergent plate boundaries in Iceland and Ethiopia | |
| Alves | Scale-relationships and geometry of normal faults reactivated during gravitational gliding of Albian rafts (Espírito Santo Basin, SE Brazil) | |
| Forte et al. | Structural geometries and magnitude of shortening in the eastern Kura fold‐thrust belt, Azerbaijan: Implications for the development of the Greater Caucasus Mountains | |
| Grimaldi et al. | Fault framework and kinematic evolution of inversion structures: Natural examples from the Neuquén Basin, Argentina | |
| Monaldi et al. | Preserved extensional structures in an inverted Cretaceous rift basin, northwestern Argentina: Outcrop examples and implications for fault reactivation | |
| Sarhan | Wrench tectonics of Abu Gharadig Basin, Western Desert, Egypt: a structural analysis for hydrocarbon prospects | |
| CN102121995A (en) | Seismic inversion reservoir stratum forecasting method for overthrust-containing complex structure | |
| Johnson et al. | Influence of fault trend, bends, and convergence on shallow structure and geomorphology of the Hosgri strike-slip fault, offshore central California | |
| CN108333624A (en) | A kind of virtual log construction method based on geophysical information | |
| CN110609319B (en) | Method for establishing time domain sequence stratigraphic section by using forward modeling of sedimentation | |
| Shuck et al. | Strike‐slip enables subduction initiation beneath a failed rift: New seismic constraints from Puysegur Margin, New Zealand | |
| CN105242307A (en) | Complex carbonate stratum earthquake porosity obtaining method and apparatus | |
| CN104975851B (en) | For amplitude with the reservoir model optimization method of geophone offset variation road set analysis | |
| Anderson et al. | Facies, architecture, and origin of a reservoir-scale sand-rich succession within submarine canyon fill: insights from Wagon Caves Rock (Paleocene), Santa Lucia Range, California, USA | |
| CN106886045A (en) | A kind of method that deformation coal forecast of distribution is carried out based on prestack Simultaneous Retrieving | |
| Molezzi et al. | Mesoarchaean-Palaeoproterozoic crustal-scale tectonics of the central Witwatersrand basin-Interpretation from 2D seismic data and 3D geological modelling | |
| Uski et al. | The 2003 earthquake swarm in Anjalankoski, south-eastern Finland | |
| CN104422955A (en) | Method for using travel-time variation quantity to extract anisotropism parameters |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170623 |
|
| WD01 | Invention patent application deemed withdrawn after publication |