CN109633743A - A method of based on waveform separation seismic facies technological prediction coal seam thickness - Google Patents

A method of based on waveform separation seismic facies technological prediction coal seam thickness Download PDF

Info

Publication number
CN109633743A
CN109633743A CN201910043205.0A CN201910043205A CN109633743A CN 109633743 A CN109633743 A CN 109633743A CN 201910043205 A CN201910043205 A CN 201910043205A CN 109633743 A CN109633743 A CN 109633743A
Authority
CN
China
Prior art keywords
seismic
coal seam
waveform
seam thickness
seismic facies
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
Application number
CN201910043205.0A
Other languages
Chinese (zh)
Inventor
杨文强
董守华
金学良
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Huaibei Mining Group Co Ltd
Original Assignee
China University of Mining and Technology CUMT
Huaibei Mining Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT, Huaibei Mining Group Co Ltd filed Critical China University of Mining and Technology CUMT
Priority to CN201910043205.0A priority Critical patent/CN109633743A/en
Publication of CN109633743A publication Critical patent/CN109633743A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles

Abstract

The invention discloses a kind of method based on waveform separation seismic facies technological prediction coal seam thickness, prediction and explanation before not disclosed suitable for coal seam thickness;The present invention utilizes the waveform separation seismic facies technology under borehole restrained, seismic facies recognition training is carried out by extracting seismic trace near well waveform as sample, and establish the mapping relations of coal seam thickness and seismic facies, based on neural network recognization work area seismic facies, and coal seam thickness information is converted by seismic facies information, instruct the prediction work of coal seam thickness.By the above-mentioned means, can be realized the prediction in zonal identical coal seam thickness section and explain work.

Description

A method of based on waveform separation seismic facies technological prediction coal seam thickness
Technical field
The present invention relates to a kind of methods based on waveform separation seismic facies technological prediction coal seam thickness, and it is pre- to belong to coal seam thickness Survey method.
Background technique
This noun of seismic facies comes from oil seismic exploration technology, eighties of last century the seventies latter stage oil seismic exploration neck The appearance of seismic stratigraphy in domain, this concept of seismic facies are gradually realized and are widely applied.1962, Sloss L L handle " phase " is defined as the summation of System Domain and its material exhibits when certain rock stratum generates.Nineteen eighty-two, Sheriff point out earthquake Mutually that seismic signature is formed by by depositional environment (such as marine facies or terrestrial facies), seismic facies it can be appreciated that sedimentary facies in the earthquake The summation showed on section.The seismic facies understood in seismic sequence is the three-dimensional space as defined by specific seismic reflection parameter Between seismic reflection unit, it is the seismic response of specific sedimentary facies or geologic body, and the unit is anti-in the earthquake of three-dimensional space It is different to penetrate feature unit adjacent thereto.
The division of seismic facies is of great significance to oil-gas exploration, and seismic facies parameter is to divide the foundation of seismic facies, also referred to as For earthquake facies marker.It is manual operations initially dividing seismic facies, it is time-consuming, it is abnormal especially to work as reflecting attribute on seismic profile When not prominent, this work is even more difficulty.Traditional seismic facies analysis method is relative to growing up in recent years quantitatively For seismic phase analysis.Traditional seismic facies analysis method is, continuity quality, frequency height strong and weak using seismic reflection amplitude Low, the reflection attributes such as geometric shape and internal structure work out seismic facies map, carry out the conversion of sedimentary facies between the opposite well of earthquake.It is this The artificial property of method is larger, has certain uncertainty.
Necessarily contain coal thickness information in seismic data, and amplitude, continuity and frequency bandwidth for being reflected in waveform etc. is a variety of Among waveform attributes feature, and seismic channel waveform is the concentrated expression of seismic response parameter (amplitude, phase, frequency etc.), only It is not enough to completely show the information of target zone by single attribute feature.
It can be eliminated using seismic waveform feature progress seismic facies classification and analyze bring using single Seismic Attribute Parameters Limitation has a good application prospect in three-dimensional work area.The principle of waveform separation seismic facies technology is according between waveform Difference, then the sample training of neural network is carried out, achieve the purpose that cluster.Waveform separation seismic facies analysis using neural network and Mode identification technology, to the actual seismic data track in a certain interval by its wave character of trace comparison and seismic properties feature, carefully Its cross directional variations is portrayed in cause, to obtain seismic anomaly planar distribution, i.e. waveform separation seismic facies map.
Summary of the invention
Goal of the invention: in order to overcome the deficiencies in the prior art, the present invention provides a kind of waveform separation seismic facies skill The method that art predicts coal seam thickness establishes contacting for coal thickness and seismic facies by waveform, achievees the purpose that predict coal seam thickness.
Technical solution: to achieve the above object, the technical solution adopted by the present invention are as follows:
A method of based on waveform separation seismic facies technological prediction coal seam thickness, include the following steps:
(1) according to the borehole data of the whole district, coal seam thickness is divided into N number of continuous coal thickness section;
(2) to n-th of coal thickness section, the most apparent R of selected characteristicn(drilling of selection is more, divides more smart for a drilling Seismic channel beside carefully) extracts model trace waveform, using model trace waveform as training sample as model trace;
(3) all training samples of neural network learning are utilized, the seismic facies distributed intelligence of the whole district is obtained;It specifically includes as follows Step:
(31) all model traces of the whole district that step (2) are extracted are divided into more than two groupings, define coal seam thickness to be predicted Position is target track, calculates the Manhattan distance between all target tracks (non-model trace) and all model traces by road, that is, counts Calculate the summation of the sampling difference value between all target track waveforms and each model trace waveform sampling point:
Wherein: A indicates model trace, AiIndicate i-th of wavelet of model trace A;B indicates target track B, BiIndicate target track B's I-th of wavelet;Target track and model trace are seismic channel, and each seismic channel has G wavelet, and M indicates model trace A and target track B Between Manhattan distance;
(32) target track is divided into Manhattan where the smallest model trace by the value for comparing Manhattan distance Group in, each target track and each model trace can only divide in a group, and according to this principle, all target tracks are divided Into different model trace groups;
(33) target track in same group is clustered, forms seismic waveform flat distribution map, i.e. seismic facies map;
(4) seismic facies map of the whole district is construed to corresponding coal thickness section, achievees the purpose that predict coal seam thickness.
Preferably, the step (4), which has, includes the following steps:
(41) response relation between the attribute and coal seam thickness of seismic channel waveform is analyzed;The seismic channel of different coal thickness reflections Waveform is inevitable different, and this phenomenon can be characterized often through certain waveform attributes features, such as the tuning effect of echo amplitude It answers, echo amplitude is correlation within λ/4 0≤H < in coal seam thickness H, and λ is seismic wave wavelength;But single category Property feature is not enough to completely show target zone information;
(42) relationship between the shape and seismic facies of seismic channel waveform is analyzed;Seismic channel waveform is the base of seismic data Plinth is the general characteristic of earthquake information, can sufficiently be excavated using seismic channel waveform separation seismic facies technology and to be contained in waveform Institute's difinite quality and quantitative information, by waveform separation seismic facies point analysis, waveform of different shapes will necessarily classify in different Seismic facies, but seismic facies technology cannot direct forecasting coal thickness information;
(43) using seismic channel waveform as medium, the mapping relations between coal seam thickness and seismic facies map are established;Drilling is taken off The seismic trace near well waveform shown carries out neural network learning as training sample and weaves seismic facies map, judges identical seismic facies area Domain is identical coal seam thickness hierarchical region.
The utility model has the advantages that the method for waveform separation seismic facies technological prediction coal seam thickness provided by the invention, with the prior art It compares, has the advantage that 1, compared with single attribute forecasting coal thickness information, the present invention has comprehensively considered each attribute in waveform Feature improves the precision of prediction;2, the present invention discloses information to drill as constraint, and traditional seismic facies qualitative analysis is carefully drawn For quantitative analysis, seismic facies information indirect is creatively construed to coal thickness information.
Detailed description of the invention
Fig. 1 is implementation flow chart of the invention;
Fig. 2 is three-dimensional wedge-shaped coal geology illustraton of model;
Fig. 3 is the part seismic channel waveform extracted from three-dimensional wedge-shaped coal geology model;
Fig. 4 is three-dimensional wedge-shaped coal geology model coal thickness plan view;
Fig. 5 is the coal thickness plan view predicted based on the present invention.
Specific embodiment
The present invention will be further explained with reference to the accompanying drawing.
A method of based on waveform separation seismic facies technological prediction coal seam thickness, include the following steps:
(1) according to the borehole data of the whole district, coal seam thickness is divided into N number of continuous coal thickness section;
(2) to n-th of coal thickness section, the most apparent R of selected characteristicn(drilling of selection is more, divides more smart for a drilling Carefully), the seismic trace near well waveform of these drillings is extracted as model trace training sample;
(3) all model trace training samples of neural network learning are utilized, the seismic facies distributed intelligence of the whole district is obtained;Specific packet Include following steps:
(31) all model traces of the whole district that step (2) are extracted are divided into more than two groupings, define coal seam thickness to be predicted Position is target track, calculates the Manhattan distance between all target tracks and all model traces by road, that is, calculates all targets The summation of sampling difference value between road waveform and each model trace waveform sampling point:
Wherein: A indicates model trace, AiIndicate i-th of wavelet of model trace A;B indicates target track B, BiIndicate target track B's I-th of wavelet;Target track and model trace are seismic channel, and each seismic channel has G wavelet, and M indicates model trace A and target track B Between Manhattan distance;The Manhattan distance of identical two model trace training samples is 0;
(32) target track is divided into Manhattan where the smallest model trace by the value for comparing Manhattan distance Group in, each target track and each model trace can only divide in a group, and according to this principle, all target tracks are divided Into different model trace groups;Classify according to similarity to all target tracks;
(33) target track in same group is clustered, forms the seismic waveform plane distribution with obvious geological Significance Figure, i.e. seismic facies map;
(4) seismic facies map of the whole district is construed to corresponding coal thickness section, achievees the purpose that predict coal seam thickness.
Further explanation is made to the present invention below with reference to example.
Step 1: 100 Inline lines are arranged in the three-dimensional wedge-shaped coal geology total 1000m of model shown in Fig. 2;Through analyzing Learn that coal seam is located at 400m, thickness change is 0~5m, density 1.41g/cm3, velocity of longitudinal wave 1820m/s.It is determined based on coal seam When window, when window need comprising roof and bottom plate information.
Step 2: coal seam thickness 0m linear increment to 5m, be divided into 0-0.7,0.7-1.4,1.4-2.0,2.0-2.6, 24 samples are extracted as model trace, mark in the eight continuous coal thickness sections 2.6-3.2,3.2-3.8,3.8-4.4,4.4-5.0 altogether Remember respective coal thickness information and extracts its waveform.
Step 3: by the Manhattan distance value in road computation model road and whole district's target track, with minimal distance principle to wave Shape is divided into eight classes, obtains the plan view with eight class seismic facies, as shown in Figure 3.
Step 4: being eight coal thickness sections, last earthquake by these eighth types of seismic facies interpretations in conjunction with the classification standard of step 2 Phasor is converted into coal thickness information distributing plan, as shown in Figure 4, Figure 5.
The above is only a preferred embodiment of the present invention, it should be pointed out that: for the ordinary skill people of the art For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered It is considered as protection scope of the present invention.

Claims (3)

1. a kind of method based on waveform separation seismic facies technological prediction coal seam thickness, characterized by the following steps:
(1) according to the borehole data of the whole district, coal seam thickness is divided into N number of continuous coal thickness section;
(2) to n-th of coal thickness section, the most apparent R of selected characteristicnSeismic channel beside a drilling extracts model as model trace Road waveform, using model trace waveform as training sample;
(3) all training samples of neural network learning are utilized, the seismic facies distributed intelligence of the whole district, the i.e. seismic facies of the whole district are obtained Figure;
(4) seismic facies map of the whole district is construed to corresponding coal thickness section, achievees the purpose that predict coal seam thickness.
2. the method according to claim 1 based on waveform separation seismic facies technological prediction coal seam thickness, it is characterised in that: The step (3), which has, to be included the following steps:
(31) all model traces of the whole district that step (2) are extracted are divided into more than two groupings, define coal seam thickness position to be predicted For target track, the Manhattan distance between all target tracks and all model traces is calculated by road, that is, calculates all target track waves The summation of sampling difference value between shape and each model trace waveform sampling point:
Wherein: A indicates model trace, AiIndicate i-th of wavelet of model trace A;B indicates target track B, BiIndicate the i-th of target track B A wavelet;Target track and model trace are seismic channel, and each seismic channel has G wavelet, and M is indicated between model trace A and target track B Manhattan distance;
(32) target track is divided into Manhattan apart from the group where the smallest model trace by the value for comparing Manhattan distance Interior, each target track and each model trace can only divide in a group, and according to this principle, all target tracks are divided into not In same model trace group;
(33) target track in same group is clustered, forms seismic waveform flat distribution map, i.e. seismic facies map.
3. the method according to claim 1 based on waveform separation seismic facies technological prediction coal seam thickness, it is characterised in that: The step (4), which has, to be included the following steps:
(41) due to different coal thickness have different seismic channel wave characters, according to drilling disclose coal seam thickness information, The response relation between seismic trace near well waveform and coal seam thickness is established, and seismic waveform by all wells is grouped;
(42) based on seismic channel waveform grouping braiding seismic facies map, therefore every kind of seismic facies represents one group of seismic channel waveform;
(43) using seismic channel waveform as medium, the mapping relations between coal seam thickness and seismic facies are established, seismic facies information is turned Turn to coal seam thickness information.
CN201910043205.0A 2019-01-17 2019-01-17 A method of based on waveform separation seismic facies technological prediction coal seam thickness Pending CN109633743A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910043205.0A CN109633743A (en) 2019-01-17 2019-01-17 A method of based on waveform separation seismic facies technological prediction coal seam thickness

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910043205.0A CN109633743A (en) 2019-01-17 2019-01-17 A method of based on waveform separation seismic facies technological prediction coal seam thickness

Publications (1)

Publication Number Publication Date
CN109633743A true CN109633743A (en) 2019-04-16

Family

ID=66062076

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910043205.0A Pending CN109633743A (en) 2019-01-17 2019-01-17 A method of based on waveform separation seismic facies technological prediction coal seam thickness

Country Status (1)

Country Link
CN (1) CN109633743A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112114360A (en) * 2019-06-21 2020-12-22 中国石油天然气集团有限公司 Seismic waveform analysis method and device
CN112213781A (en) * 2020-07-30 2021-01-12 中国煤炭地质总局地球物理勘探研究院 Method and system for predicting coal seam thickness under big data
CN112379442A (en) * 2020-11-02 2021-02-19 中国石油天然气集团有限公司 Seismic waveform classification method and device
CN112684497A (en) * 2019-10-17 2021-04-20 中国石油天然气集团有限公司 Seismic waveform clustering method and device

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1975462A (en) * 2006-09-05 2007-06-06 孟召平 Coal seam thickness analyzing method based on earthquake attribute
CN102478668A (en) * 2010-11-30 2012-05-30 中国石油天然气集团公司 Method for applying seismic multiattribute parameters to predicting coal seam thickness
CN102650702A (en) * 2012-05-03 2012-08-29 中国石油天然气股份有限公司 Seismic waveform analysis and reservoir prediction method and device
CN103487832A (en) * 2013-09-12 2014-01-01 电子科技大学 Method for classifying supervised waveforms in three-dimensional seismic signal
CN105068118A (en) * 2015-06-10 2015-11-18 中国矿业大学 Method for interpreting range of magmatic rock intrusion into coal layer based on seismic information

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1975462A (en) * 2006-09-05 2007-06-06 孟召平 Coal seam thickness analyzing method based on earthquake attribute
CN102478668A (en) * 2010-11-30 2012-05-30 中国石油天然气集团公司 Method for applying seismic multiattribute parameters to predicting coal seam thickness
CN102650702A (en) * 2012-05-03 2012-08-29 中国石油天然气股份有限公司 Seismic waveform analysis and reservoir prediction method and device
CN103487832A (en) * 2013-09-12 2014-01-01 电子科技大学 Method for classifying supervised waveforms in three-dimensional seismic signal
CN105068118A (en) * 2015-06-10 2015-11-18 中国矿业大学 Method for interpreting range of magmatic rock intrusion into coal layer based on seismic information

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
姚维益: "地震属性技术研究及其在胜利富林洼陷储层预测中的应用", 《中国优秀硕士学位论文全文数据库•基础科学辑》 *
张许良等: "煤层裂隙探测的地震波形分类技术研究", 《河南理工大学学报(自然科学版)》 *
李治欣等: "地震多属性技术解释煤矿薄煤区的应用", 《科技创新与应用》 *
秦永军 等: "地震波形分类技术在煤层分叉解释中的应用", 《中国煤炭地质》 *
袁峰: "永陇矿区煤层厚度预测研究", 《中国地质大学硕士专业学位论文》 *
郭彦省等: "地震属性及其在煤层厚度预测中的应用", 《中国矿业大学学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112114360A (en) * 2019-06-21 2020-12-22 中国石油天然气集团有限公司 Seismic waveform analysis method and device
CN112684497A (en) * 2019-10-17 2021-04-20 中国石油天然气集团有限公司 Seismic waveform clustering method and device
CN112684497B (en) * 2019-10-17 2023-10-31 中国石油天然气集团有限公司 Seismic waveform clustering method and device
CN112213781A (en) * 2020-07-30 2021-01-12 中国煤炭地质总局地球物理勘探研究院 Method and system for predicting coal seam thickness under big data
CN112213781B (en) * 2020-07-30 2023-04-25 中国煤炭地质总局地球物理勘探研究院 Method and system for predicting thickness of coal seam under big data
CN112379442A (en) * 2020-11-02 2021-02-19 中国石油天然气集团有限公司 Seismic waveform classification method and device

Similar Documents

Publication Publication Date Title
CN109061765B (en) Trap evaluation method for heterogeneous thin sandstone interbed reservoir
US6735526B1 (en) Method of combining directional seismic attributes using a supervised learning approach
CN107356958B (en) A kind of fluvial depositional reservoir substep seismic facies prediction technique based on geological information constraint
CN109633743A (en) A method of based on waveform separation seismic facies technological prediction coal seam thickness
CN106842307B (en) Reservoir stratum fine prediction method based on waveform classification and retrieval under forward constraint
CN108802812A (en) A kind of formation lithology inversion method of well shake fusion
CN102650702A (en) Seismic waveform analysis and reservoir prediction method and device
US7519476B1 (en) Method of seismic interpretation
CN105445800A (en) Thick sand body top differentiation lithologic reservoir identification method
CN104698497A (en) Method for fine calibration of position of sand shale thin interbed
CN108415075A (en) A kind of method for predicting reservoir under condition of less wells
CN109725348A (en) A method of sedimentary facies is identified based on seismic data
CN111239824A (en) Volcanic lithofacies prediction method
CN111175819A (en) Gravel rock sector sedimentary facies belt fine dividing method based on well-seismic multi-stage constraint
CN112505754B (en) Method for collaborative partitioning sedimentary microfacies by well-seismic based on high-precision sequence grid model
CN116168224A (en) Machine learning lithology automatic identification method based on imaging gravel content
CN104977611A (en) Reef reservoir engraving method
CN110927817A (en) Method for optimally designing and controlling thin-layer superposed riverway horizontal well track
CN109283577B (en) Seismic horizon calibration method
CN104714248B (en) Formation overlap-denudation line recognition method based on iteration earthquake DNA detection
CN110795513B (en) Method for predicting distribution of river facies source storage ectopic type compact oil gas dessert area
CN115561814A (en) Exploration method of shallow stratum-lithologic reservoir in subsidence basin gentle slope zone
CN105259576A (en) Hydrocarbon reservoir identification method by means of seismic statistical characteristic
RU2158939C1 (en) Method of search for oil and gas deposits with use of neurocomputer system for processing of seismic prospecting data
CN112444868B (en) Seismic phase analysis method based on improved K-means algorithm

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190416