CN105652317B - High-resolution seismic imaging method and device - Google Patents
High-resolution seismic imaging method and device Download PDFInfo
- Publication number
- CN105652317B CN105652317B CN201610074937.2A CN201610074937A CN105652317B CN 105652317 B CN105652317 B CN 105652317B CN 201610074937 A CN201610074937 A CN 201610074937A CN 105652317 B CN105652317 B CN 105652317B
- Authority
- CN
- China
- Prior art keywords
- imaging
- sampling point
- amplitude sampling
- point
- seismic
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 265
- 238000005070 sampling Methods 0.000 claims abstract description 143
- 238000000034 method Methods 0.000 claims abstract description 134
- 230000005012 migration Effects 0.000 claims abstract description 32
- 238000013508 migration Methods 0.000 claims abstract description 32
- 238000012545 processing Methods 0.000 claims abstract description 15
- 230000008569 process Effects 0.000 claims description 74
- 239000011159 matrix material Substances 0.000 claims description 36
- 230000008859 change Effects 0.000 claims description 12
- 239000013598 vector Substances 0.000 claims description 12
- 239000000203 mixture Substances 0.000 claims description 8
- 239000000284 extract Substances 0.000 claims description 6
- 238000000354 decomposition reaction Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 3
- 108010001267 Protein Subunits Proteins 0.000 claims description 2
- 239000003245 coal Substances 0.000 abstract description 7
- 238000010276 construction Methods 0.000 abstract description 6
- 238000005516 engineering process Methods 0.000 description 12
- 238000004458 analytical method Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 230000008901 benefit Effects 0.000 description 5
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000003325 tomography 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
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a high-resolution seismic imaging method and device and relates to the field of seismic exploration. The method comprises the following steps: obtaining the single-shot data and migration velocity data for imaging processing; calculating the seismic travel time information corresponding to the single-shot data according to the detector spatial position and the migration velocity data; selecting amplitude sampling points in the single-shot data to form a common image point gather according to the position point identifier consisting of the detector spatial position and the seismic travel time information; and calculating the common image point gather according to a multiple-signal classifying method to obtain a high-resolution seismic imaging result corresponding to the single-shot data. By calculating the obtained common image point gathered by the multiple-signal classifying method, the stable-phase point can be accurately judged, and the seismic waveform is selected accordingly; the high-resolution geological information is excavated from the collected data with low signal-to-noise ratio, and accurate seismic data is provided for the structure identification in coal exploration; and the method is of great significance to lowering the investment cost and reducing the construction risk.
Description
Technical field
The present invention relates to field of seismic exploration, in particular to a kind of high-resolution seismic exploration wave imaging method and device.
Background technology
Seismic prospecting is i.e. using underground medium elasticity and the difference of density, the ground produced by observation and analysis artificial earthquake
Propagation law of the seismic wave in underground, infers the geophysical exploration method of the property and form of subterranean strata；Wherein, above-mentioned analysis
Seismic wave mainly carries out imaging to the seismic wave that artificial earthquake is produced in the method for the propagation law of underground, then over the ground
The imaging results of seismic wave are analyzed.
At present, in field of seismic exploration, the Kirchhoff migration technologies based on Integral Solution of Wave Equation are still main flow
Seismic imaging technology (Leveille et al, 2011), the advantage of the technology be it is adaptable to complicated observation system,
And can realize that specific aim target area is imaged；But, the technology equally exists certain shortcoming, i.e., when the earthquake money of collection
When material signal to noise ratio is relatively low, seismic imaging is carried out using above-mentioned Kirchhoff migration technologies, often so that image quality is paid no attention to
Think.
In order to solve the above problems, multiple signal classification method is adopted at present, multiple signal classification method is in geological data
Process field has a well application, such as velocity analysis (Kirlin, 1992) and common reflection surface imaging (Gelius et al.,
2013), above-mentioned multiple signal classification method is mainly used in the effective band for estimating signal, removes noise jamming.Will be above-mentioned multiple
During Modulation recognition method is applied to seismic imaging, seismic imaging is carried out, it is possible to use its useful signal estimates advantage,
Compacting noise improves imaging resolution.Multiple signal classification method was once applied to seismic imaging by Gelius (2014), but was tied
Fruit shows that the method does not solve steady phase point and determines and seismic waveform On The Choice.
Inventor has found under study for action, above-mentioned Kirchhoff migration technologies and multiple signal classification method seismic wave into
Image space face is respectively present not good for the low image quality of seismic data of noise and cannot solve steady phase point determination and seismic waveform
The problem of selection, for the problems referred to above, not yet proposes effective settling mode at present.
The content of the invention
It is an object of the invention to provide a kind of high-resolution seismic exploration wave imaging method and device, can realize effective earthquake
The estimation of signal, improves seismic imaging resolution ratio.
In a first aspect, a kind of high-resolution seismic exploration wave imaging method is embodiments provided, including：
Obtain the corresponding migration velocity data of single-shot data and single-shot data of pending imaging；Wherein, single-shot number
According to the multiple seismic waves for reflecting including subterranean strata interface or reflect, seismic wave carries geological information；
According to the corresponding wave detector locus of single-shot data and migration velocity data, the shot position of single-shot data is calculated
Seimic travel time information of any one the imaging point position of Jing subsurface imagings space to geophone station position；Wherein, single-shot data pair
The wave detector answered is for multiple, and each wave detector is arranged on corresponding geophone station position；
According to the position point identification that wave detector locus and seimic travel time information are constituted, position is chosen from single-shot data
The corresponding amplitude sampling point of point identification is put, the corresponding common imaging gather of single-shot data is obtained；
Calculating process is carried out to common imaging gather according to multiple signal classification method, the corresponding high score of single-shot data is obtained
Resolution seismic imaging result.
With reference in a first aspect, embodiments provide the first possible embodiment of first aspect, wherein, root
Calculating process is carried out to common imaging gather according to multiple signal classification method, the corresponding high-resolution seismic exploration ripple of single-shot data is obtained
Imaging results include：
Along the horizontal direction of amplitude sampling point, to the corresponding amplitude sample in each imaging point position in common imaging gather
Point carries out scan process, obtains amplitude sampling point in the corresponding laterally main energy band scope of horizontal direction；Wherein, horizontal main energy
The corresponding amplitude sampling point sequence of band scope carries geological information；
Along the longitudinal direction of amplitude sampling point, earthquake is carried out to each the amplitude sampling point in the range of horizontal main energy band
The judgement of ripple size, obtains the longitudinal direction effectively seismic wave magnitude range corresponding to horizontal main energy band scope；
According to the imaging array model of the effective seismic wave magnitude range composition of horizontal main energy band scope and its corresponding longitudinal direction
Enclose, the corresponding imaging array of imaging array scope is extracted from common imaging gather；
Useful signal estimation process is carried out to imaging array, the corresponding high-resolution seismic exploration ripple imaging knot of single-shot data is obtained
Really.
With reference to the first possible embodiment of first aspect, second of first aspect is embodiments provided
Possible embodiment, wherein, place is scanned to the corresponding amplitude sampling point in each imaging point position in common imaging gather
Reason, obtaining amplitude sampling point includes in the corresponding laterally main energy band scope of horizontal direction：
Along the horizontal direction of amplitude sampling point, each corresponding amplitude in imaging point position in common imaging gather is determined
Transverse energy band scope of the sampling point in horizontal direction；
Summation process is carried out to all amplitude sampling points in the range of each transverse energy band, is obtained corresponding to each
Amplitude sampling point superposition value of the individual transverse energy with scope；
Relatively amplitude sampling point superposition value of each transverse energy with scope, chooses maximum amplitude sampling point superposition value
Corresponding transverse energy band scope is horizontal main energy band scope.
With reference to second possible embodiment of first aspect, embodiments provide first aspect the third
Possible embodiment, wherein, along the longitudinal direction of amplitude sampling point, to each amplitude in the range of horizontal main energy band
Sampling point carries out seismic wave size judgement, obtains the longitudinal direction effectively seismic wave magnitude range bag corresponding to horizontal main energy band scope
Include：
Along the longitudinal direction of amplitude sampling point, zero point is carried out to each the amplitude sampling point in the range of horizontal main energy band
Judge with change in polarity；
Detect the zero point in a longitudinal direction of each the amplitude sampling point in the range of horizontal main energy band and including
During the positive and negative extreme value of zero point, determine the corresponding longitudinal magnitude range of positive and negative extreme value be corresponding amplitude sampling point longitudinal direction effectively
Seismic wave magnitude range.
With reference to the third possible embodiment of first aspect, the 4th kind of first aspect is embodiments provided
Possible embodiment, wherein, useful signal estimation process is carried out to imaging array, the corresponding high-resolution of single-shot data is obtained
Seismic imaging result includes：
The method that the longitudinal direction amount of predetermined number is chosen according to the line direction in imaging array, is carried out at decomposition to imaging array
Reason, obtains multiple imaging submatrixs；Wherein, at least each two is continuously imaged submatrix and includes that identical is longitudinally measured；
The covariance matrix of each imaging submatrix is calculated respectively, obtains the sub- covariance corresponding to each imaging submatrix
Matrix；
The many sub- covariance matrixes to obtaining carry out average computation process, obtain average covariance matrices；
Characteristic vector R of average covariance matrices is calculated by below equationaV=λ v；Wherein, RaFor average covariance square
Battle array, characteristic vectors of the v for average covariance matrices, characteristic values of the λ for average covariance matrices；
Singular value decomposition process is carried out to average covariance matrices by below equation
To earthquake effective signal subspace matrix and noise subspace matrix；
According to the characteristic vector and noise subspace matrix of average covariance matrices, calculating includes seismic imaging result
The pseudo- spectrum of multiple signal classification.
With reference to the 4th kind of possible side of the first possible embodiment to first aspect of first aspect, first aspect
Any one possible embodiment in formula, embodiments provides the 5th kind of possible embodiment of first aspect,
Wherein, methods described also includes：
Obtain the high-resolution seismic exploration ripple imaging results of multiple single-shot data；
Process is overlapped to the high-resolution seismic exploration ripple imaging results of multiple single-shot data, shot gather data is obtained corresponding
High-resolution seismic exploration ripple imaging results；Wherein, shot gather data includes multiple single-shot data.
Second aspect, the embodiment of the present invention additionally provide a kind of high-resolution seismic exploration ripple imaging device, and device includes：
First acquisition module, for obtaining the single-shot data and the corresponding migration velocity of single-shot data of pending imaging
Data；Wherein, single-shot data include the reflection of subterranean strata interface or the multiple seismic waves for reflecting, and seismic wave carries geology letter
Breath；
Seimic travel time information computational module, for according to the corresponding wave detector locus of single-shot data and migration velocity
Data, calculate the earthquake of any one imaging point position of shot position Jing subsurface imagings space of single-shot data to geophone station position
Information when ripple is walked；Wherein, the corresponding wave detector of single-shot data is for multiple, and each wave detector is arranged at corresponding geophone station position
Put；
Common imaging gather sets up module, for the position constituted according to wave detector locus and seimic travel time information
Point identification, from the corresponding amplitude sampling point of chosen position point identification in single-shot data, obtains the corresponding imaging point altogether of single-shot data
Road collection；
Processing module is calculated, for carrying out calculating process to common imaging gather according to multiple signal classification method, is obtained
The corresponding high-resolution seismic exploration ripple imaging results of single-shot data.
With reference to second aspect, the first possible embodiment of second aspect is embodiments provided, wherein, meter
Calculating processing module includes：
Scan process unit, for the horizontal direction along amplitude sampling point, to each imaging point in common imaging gather
The corresponding amplitude sampling point in position is scanned process, obtains amplitude sampling point in the corresponding laterally main energy band model of horizontal direction
Enclose；Wherein, laterally the corresponding amplitude sampling point sequence of main energy band scope carries geological information；
Seismic wave size judging unit, for the longitudinal direction along amplitude sampling point, in the range of horizontal main energy band
Each amplitude sampling point carry out seismic wave size judgement, obtain the longitudinal direction effectively seismic wave corresponding to horizontal main energy band scope
Magnitude range；
Imaging array extraction unit, for according to horizontal main energy band scope and its effective seismic wave size in corresponding longitudinal direction
The imaging array scope of scope composition, extracts the corresponding imaging array of imaging array scope from common imaging gather；
Estimate processing unit, for useful signal estimation process is carried out to imaging array, obtain the corresponding height of single-shot data
Resolution seismic ripple imaging results.
With reference to the first possible embodiment of second aspect, second of second aspect is embodiments provided
Possible embodiment, wherein, scan process unit includes：
First determination subelement, for the horizontal direction along amplitude sampling point, determine in common imaging gather each into
Transverse energy band scope of the corresponding amplitude sampling point of image point position in horizontal direction；
Summation process subelement, for suing for peace to all amplitude sampling points in the range of each transverse energy band
Process, obtain the amplitude sampling point superposition value with scope corresponding to each transverse energy；
Comparing subunit, for comparing amplitude sampling point superposition value of each transverse energy with scope；
Subelement is chosen, the corresponding transverse energy band scope of the amplitude sampling point superposition value for choosing maximum is laterally main
Energy band scope.
With reference to second possible embodiment of second aspect, embodiments provide second aspect the third
Possible embodiment, wherein, seismic wave size judging unit includes：
Zero point and change in polarity judgment sub-unit, for the longitudinal direction along amplitude sampling point, to horizontal main energy band
In the range of each amplitude sampling point carry out zero point and change in polarity and judge；
Second determination subelement, for detecting each the amplitude sampling point in the range of horizontal main energy band in longitudinal direction
During zero point and the positive and negative extreme value including zero point on direction, determine that the corresponding longitudinal magnitude range of positive and negative extreme value is corresponding amplitude
The longitudinal direction of value sampling point effectively seismic wave magnitude range.
A kind of high-resolution seismic exploration wave imaging method provided in an embodiment of the present invention and device, including：Obtain it is pending into
As the single-shot data and the corresponding migration velocity data of single-shot data of process；According to the corresponding wave detector locus of single-shot data
With migration velocity data, the corresponding seimic travel time information of single-shot data is calculated；Then, according to wave detector locus and above-mentioned
The position point identification of seimic travel time information composition, from the corresponding amplitude sampling point of chosen position point identification in single-shot data, group
Into common imaging gather；Finally, calculating process is carried out to above-mentioned common imaging gather according to multiple signal classification method, obtains list
The corresponding high-resolution seismic exploration ripple imaging results of big gun data；
The low image quality of seismic data of noise is not good and nothing with existing in terms of seismic imaging in prior art
Method solves the problems, such as that steady phase point determines to choose with seismic waveform and compares which adopts multiple signal classification method inclined to Kirchhoff
The common imaging point road for obtaining for moving into picture carries out calculating process, can accurately judge steady phase point and thus choose seismic waveform, from
And the estimation of effective seismic signal can be realized, improve seismic imaging resolution ratio；In its gathered data from low signal-to-noise ratio,
High-resolution geological information is excavated, and more accurate seismic data is provided for the conformation identification in Coal Exploration exploitation,
Reducing input cost and reducing significant on construction risk.
To enable the above objects, features and advantages of the present invention to become apparent, preferred embodiment cited below particularly, and coordinate
Appended accompanying drawing, is described in detail below.
Description of the drawings
In order to be illustrated more clearly that the technical scheme of the embodiment of the present invention, below by to be used attached needed for embodiment
Figure is briefly described, it will be appreciated that the following drawings illustrate only certain embodiments of the present invention, thus be not construed as it is right
The restriction of scope, for those of ordinary skill in the art, on the premise of not paying creative work, can be with according to this
A little accompanying drawings obtain other related accompanying drawings.
The flow chart that Fig. 1 shows a kind of high-resolution seismic exploration wave imaging method that the embodiment of the present invention is provided；
The flow chart that Fig. 2 shows another kind of high-resolution seismic exploration wave imaging method provided by the embodiment of the present invention；
The flow chart that Fig. 3 shows another kind of high-resolution seismic exploration wave imaging method provided by the embodiment of the present invention；
The flow chart that Fig. 4 shows another kind of high-resolution seismic exploration wave imaging method provided by the embodiment of the present invention；
The flow chart that Fig. 5 shows another kind of high-resolution seismic exploration wave imaging method provided by the embodiment of the present invention；
Fig. 6 shows the schematic diagram of the imaging results of traditional Kirchhoff migration and imaging techniques；
Fig. 7 shows that a kind of imaging results of high-resolution seismic exploration wave imaging method that the embodiment of the present invention is provided are illustrated
Figure；
Fig. 8 shows a kind of structural representation of high-resolution seismic exploration ripple imaging device that the embodiment of the present invention is provided；
Fig. 9 calculates processing module in showing a kind of high-resolution seismic exploration ripple imaging device that the embodiment of the present invention is provided
With the structural representation of scan process unit；
Figure 10 shows seismic wave size in a kind of high-resolution seismic exploration ripple imaging device that the embodiment of the present invention is provided
The structural representation of judging unit and estimation processing unit.
Specific embodiment
Below in conjunction with accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Ground description, it is clear that described embodiment is only a part of embodiment of the invention, rather than the embodiment of whole.Generally exist
The component of the embodiment of the present invention described and illustrated in accompanying drawing can be arranged and be designed with a variety of configurations herein.Cause
This, the detailed description of the embodiments of the invention to providing in the accompanying drawings is not intended to limit claimed invention below
Scope, but it is merely representative of the selected embodiment of the present invention.Based on embodiments of the invention, those skilled in the art are not doing
The every other embodiment obtained on the premise of going out creative work, belongs to the scope of protection of the invention.
In prior art, Kirchhoff migration technologies and multiple signal classification method are deposited in terms of seismic imaging respectively
Not good for the low image quality of seismic data of noise and cannot solve the problems, such as that steady phase point determines and seismic waveform is chosen；Pin
To the problem, a kind of high-resolution seismic exploration wave imaging method and device are embodiments provided, Kirchhoff is imaged by which
Technology is combined with multiple signal classification method, has reached the purpose for improving seismic imaging resolution ratio；The technology of the present invention
Advantage is to be determined by steady phase point and the selection of seismic waveform self adaptation, realizes the earthquake wave height based on multiple signal classification
Resolution imaging, improves the detectivity of underground space geologic body.
With reference to Fig. 1, a kind of high-resolution seismic exploration wave imaging method is embodiments provided, methods described includes as follows
Step：
S101, the single-shot data and the corresponding migration velocity data of the single-shot data that obtain pending imaging；Its
In, the single-shot data include the reflection of subterranean strata interface or the multiple seismic waves for reflecting, and the seismic wave carries described
Geological information.
Specifically, focus and wave detector are laid in target area, from focus to underground Transmitted station, by the focus correspondence
Multiple wave detectors receive the seismic wave through subsurface reflective or reflect corresponding seismic wave as single-shot data, due to the list
The seismic wave that big gun data include carries geological information, therefore can carry out imaging to above-mentioned single-shot data, with according to imaging
As a result the work such as coal karst collapse col umn identification are carried out.
Above-mentioned migration velocity data carry out migration velocity analysis process and obtain to above-mentioned single-shot data；Wherein, it is above-mentioned
Migration velocity data, corresponding to formation velocity parameter, for calculating when seimic wave propagation is travelled (i.e. time), which is to single-shot number
According to carrying out migration velocity analysis acquisition；And above-mentioned migration velocity analysis is a kind of velocity modeling technology, the velocity modeling technology is led to
The seismic wave focusing crossed in analysis seismic acquisition data determines speed parameter, that is, obtain migration velocity data.
The embodiment of the present invention, the single-shot data and the single-shot data for obtaining above-mentioned pending imaging first are corresponding partially
Speed data is moved, to carry out subsequent treatment.
S102, according to the corresponding wave detector locus of the single-shot data and the migration velocity data, calculate described
Believe to the seimic travel time of geophone station position any one imaging point position of shot position Jing subsurface imagings space of single-shot data
Breath；Wherein, the corresponding wave detector of the single-shot data is for multiple, and each described wave detector is arranged at the corresponding detection
On point position.
Specifically, according to the corresponding wave detector locus of each single-shot data and the corresponding migration velocity of single-shot data
Data carry out seimic travel time information calculating, and seimic travel time information is tried to achieve by ray tracing eikonal equation, and ray tracing journey
Eikonal equation is the differential relationship of speed, locus and seimic travel time；Wherein, the locus for referring at this refers to single-shot number
According to corresponding wave detector locus.
Then according to the single-shot data and corresponding migration velocity data, the seismic wave shot point of the single-shot data is calculated first
Position (i.e. shot position) to any one imaging point position of subsurface imaging space it is corresponding walk when (in order to distinguish following walking
When, when becoming first herein and walking), as subsurface imaging space includes multiple imaging points, therefore the big gun of the single-shot data for calculating
Point position is also multiple when walking to any one imaging point position of subsurface imaging space corresponding first；
Then any one imaging point position of subsurface imaging space is calculated again to the seismic wave geophone station position of the single-shot data
Put (i.e. geophone station position) it is multiple walk when while walking (referred to herein as second), wherein, each second is walked with one first when walking
When it is corresponding.
The single-shot data pair are obtained by being added when corresponding second walks when walking with first when finally walk above-mentioned first
The seimic travel time answered；The seimic travel time for multiple, multiple seimic travel times correspond to respectively different subsurface imaging spaces into
Picture point.
S103, the position point identification constituted according to the wave detector locus and the seimic travel time information, from institute
The corresponding amplitude sampling point of the position point identification is chosen in stating single-shot data, the corresponding imaging point altogether of the single-shot data is obtained
Road collection.
Specifically, matrix of the single-shot data equivalent to one including multiple data values (i.e. amplitude samples point), and the single-shot
During the seimic travel time information that the wave detector locus of data and the wave detector are calculated respectively represents above-mentioned matrix respectively
(wherein, wave detector locus and corresponding seimic travel time information can be with abstract for the locus of data value and time location
For the line number and row number of above-mentioned matrix), and according to a locus and time location can with positional matrix in a data
Value (i.e. amplitude samples point), then finds the amplitude samples point from the corresponding matrix of single-shot data, and extracts the amplitude samples
Point；
Due to the corresponding wave detector of single-shot data and seimic travel time information be it is multiple, therefore according to the space bit of wave detector
Put and multiple data values (i.e. amplitude sampling point) can be positioned with corresponding seimic travel time, then from the corresponding matrix of single-shot data
In find and extract all corresponding amplitude samples points；By it is all amplitude samples points composition single-shot data extracted corresponding altogether into
Picture point road collection.Wherein, the different underground space imaging point of above-mentioned seimic travel time correspondence.
S104, calculating process is carried out to the common imaging gather according to multiple signal classification method, obtain the single-shot
The corresponding high-resolution seismic exploration ripple imaging results of data.
Specifically, by multiple signal classification method modification image function, (i.e. subsequent extracted is imaged square to the embodiment of the present invention
Battle array), then the image function (i.e. imaging array) to changing carries out useful signal estimation process, obtains the single-shot data correspondence
High-resolution seismic exploration ripple imaging results.
The function of geological information is mainly carried due to changing image function, therefore useful signal estimation is carried out to the function
Process, seismic imaging result is obtained for high-resolution seismic exploration ripple imaging results.
In seismic imaging in a kind of high-resolution seismic exploration wave imaging method provided in an embodiment of the present invention, with prior art
It is good and cannot solve that steady phase point determines and seismic waveform chooses that aspect does not exist for the low image quality of seismic data of noise
Problem is compared, and which adopts multiple signal classification method to calculate the common imaging point road for obtaining of Kirchhoff migration imagings
Process, can accurately judge steady phase point and thus choose seismic waveform such that it is able to realize the estimation of effective seismic signal, improve
Seismic imaging resolution ratio；In its gathered data from low signal-to-noise ratio, high-resolution geological information is excavated, be that coal is surveyed
The conformation identification visited in exploitation provides more accurate seismic data, has with reducing input cost is reduced on construction risk
Significance.
Used in the embodiment of the present invention, multiple signal classification method can excavate high score from low signal-to-noise ratio gathered data
Resolution geological information, provides more accurate seismic data for the conformation identification in Coal Exploration exploitation, is reducing input cost
It is significant on construction risk with reducing.Specifically, with reference to Fig. 2, in above-mentioned steps 104, according to multiple signal classification side
Method carries out calculating process to the common imaging gather, obtains the corresponding high-resolution seismic exploration ripple imaging results of the single-shot data
Concrete steps include：
S201, along the horizontal direction of the amplitude sampling point, to each imaging point position in the common imaging gather
Corresponding amplitude sampling point is scanned process, obtains the amplitude sampling point in the corresponding laterally main energy of the horizontal direction
Band scope；Wherein, the corresponding amplitude sampling point sequence of the laterally main energy band scope carries the geological information.
Specifically, common imaging gather is equally a matrix, and which is from the corresponding matrix of single-shot data to extract specific
The amplitude sampling point of position (i.e. the position of the locus of wave detector and seimic travel time message identification), then constitutes new square
Battle array；Above-mentioned scan process i.e. in the new matrix, along the horizontal direction of the amplitude sampling point, with each amplitude sample
The amplitude sampling point of the horizontal peripherad predetermined number of the sampling point centered on point, is selected, multigroup amplitude sampling point sequence is obtained, then
Summation process is carried out to every group of amplitude sampling point sequence, the maximum corresponding horizontal model of amplitude sampling point sequence of summed result is chosen
Enclose for horizontal main energy band scope.
Wherein, laterally the center of main energy band scope is steady phase point position, i.e. space corresponding to earthquake useful signal
Position.
S202, along the longitudinal direction of the amplitude sampling point, to each amplitude in the range of the laterally main energy band
Value sampling point carries out seismic wave size judgement, obtains the longitudinal direction effectively seismic wave size model corresponding to the laterally main energy band scope
Enclose.
Specifically, along the longitudinal direction of the corresponding amplitude sampling point sequence of horizontal main energy band scope for selecting, to this
Longitudinal sequence that each amplitude sampling point is located in sequence carries out zero point and change in polarity judges, detect zero point and including
During the positive and negative extreme value of the zero point, the scope of longitudinal sequence that the amplitude sampling point is located has as the longitudinal direction of the amplitude sampling point
Effect seismic wave magnitude range.
According to the method described above, for the amplitude sampling point in the range of each laterally main energy band, you can obtain corresponding to
Laterally the longitudinal direction of main energy band scope effectively seismic wave magnitude range, obtains shaking corresponding to the longitudinal direction of transverse-vibration amplitude sampling point sequence
Amplitude sampling point sequence.
S203, according to the effective seismic wave magnitude range composition of the laterally main energy band scope and its corresponding longitudinal direction
Imaging array scope, the corresponding imaging array of the imaging array scope is extracted from the common imaging gather.
Specifically, above-mentioned transverse-vibration amplitude sampling point sequence and pitch amplitude value sampling point sequence constitute an imaging array,
Then the imaging array is extracted from common imaging gather.Wherein, geological information is carried in the imaging array.
S204, useful signal estimation process is carried out to the imaging array, obtain the corresponding high-resolution of the single-shot data
Rate seismic imaging result.
Specifically, the imaging array to carrying geological information carries out imaging, you can obtain single-shot data corresponding
High-resolution seismic exploration ripple imaging results.
In the embodiment of the present invention, with reference to Fig. 3, the specific scanning process in above-mentioned steps 201, including following step
Suddenly：
S301, along the horizontal direction of the amplitude sampling point, determine each imaging point position in the common imaging gather
Put transverse energy band scope of the corresponding amplitude sampling point in the horizontal direction.
S302, summation process is carried out to all amplitude sampling points in the range of transverse energy band each described, obtained
The amplitude sampling point superposition value with scope corresponding to transverse energy each described.
S303, compare amplitude sampling point superposition value of each described transverse energy with scope, choose maximum amplitude
The corresponding transverse energy band scope of sampling point superposition value is the laterally main energy band scope.
Corresponding to above-mentioned steps 301-303, in the corresponding matrix of common imaging gather, along the amplitude sampling point
Horizontal direction, centered on each amplitude sampling point, selects the amplitude sampling point of the horizontal peripherad predetermined number of the sampling point, obtains
To the corresponding amplitude sampling point sequence of multigroup transverse energy band scope.
Summation process is carried out to above-mentioned every group of amplitude sampling point sequence, summed result (i.e. amplitude sampling point superposition value) is chosen
The maximum corresponding lateral extent of amplitude sampling point sequence is horizontal main energy band scope.
Above-mentioned predetermined number can be set as needed, and the embodiment of the present invention is not particularly limited to this.
In the embodiment of the present invention, the concrete determination method of the effectively seismic wave magnitude range of the longitudinal direction in above-mentioned steps 202, bag
Include following steps：
Along the longitudinal direction of the amplitude sampling point, to each the amplitude sampling point in the range of the laterally main energy band
Carry out zero point and change in polarity judges；
Detecting zero in the longitudinal direction of each the amplitude sampling point in the range of the laterally main energy band
When point and the positive and negative extreme value including the zero point, determine the corresponding longitudinal magnitude range of the positive and negative extreme value and shake described in corresponding
The longitudinal direction of amplitude sampling point effectively seismic wave magnitude range.
Specifically, along the longitudinal direction of the corresponding amplitude sampling point sequence of horizontal main energy band scope for selecting, to this
Longitudinal sequence that each amplitude sampling point is located in sequence carries out zero point and change in polarity judges, detect zero point and including
During the positive and negative extreme value of the zero point, the scope of longitudinal sequence that the amplitude sampling point is located has as the longitudinal direction of the amplitude sampling point
Effect seismic wave magnitude range.
According to the method described above, for the amplitude sampling point in the range of each laterally main energy band, you can obtain corresponding to
Laterally the longitudinal direction of main energy band scope effectively seismic wave magnitude range, obtains shaking corresponding to the longitudinal direction of transverse-vibration amplitude sampling point sequence
Amplitude sampling point sequence.
In the embodiment of the present invention, with reference to Fig. 4, in above-mentioned steps 204, the imaging array is carried out at useful signal estimation
The concrete grammar of reason comprises the steps：
S401, longitudinal method measured that predetermined number is chosen according to the line direction in the imaging array, to the imaging
Matrix carries out resolution process, obtains multiple imaging submatrixs；Wherein, at least continuous imaging submatrix of each two includes phase
Same longitudinal direction amount.
Specifically, the longitudinal direction amount of selection is preferably and has overlap, such as longitudinal direction amount 1,2,3,4,5,6,7, the imaging square
The line direction of battle array can choose 1234,2345,3456,4567 respectively as imaging submatrix；Wherein, the quantity of above-mentioned longitudinal direction amount
Can be configured as needed, the present invention is not particularly limited to this.
S402, the covariance matrix for calculating each imaging submatrix respectively, obtain sub corresponding to being imaged each described
The sub- covariance matrix of matrix.
S403, the multiple described sub- covariance matrix to obtaining carry out average computation process, obtain average covariance matrices.
S404, characteristic vector R that average covariance matrices are calculated by below equationaV=λ v；Wherein, RaFor average association side
Difference matrix, characteristic vectors of the v for average covariance matrices, characteristic values of the λ for average covariance matrices.
S405, below equation excessively carry out singular value decomposition process to the average covariance matricesObtain earthquake effective signal subspace matrix and noise subspace matrix.
S404, the characteristic vector according to the average covariance matrices and the noise subspace matrix, calculating include institute
State the pseudo- spectrum of multiple signal classification of seismic imaging result.
For field of seismic exploration, in practice, target area is to include multiple focus and multiple wave detectors, therefore for every
One target area, including the earthquake shot gather data (i.e. shot gather data) of multiple single-shot data carries complete geological information,
And each single-shot data typically just carries part geological information, with reference to Fig. 5, method provided in an embodiment of the present invention is also wrapped
Include：
S501, the high-resolution seismic exploration ripple imaging results for obtaining multiple single-shot data；
S502, the high-resolution seismic exploration ripple imaging results to multiple single-shot data are overlapped process, obtain big gun collection
The corresponding high-resolution seismic exploration ripple imaging results of data；Wherein, the shot gather data includes multiple single-shot data.
As the earthquake shot gather data for obtaining is only the data of the complete geological information for carrying reaction target area, and earthquake
Each single-shot data in shot gather data simply carry a part of data of the geological information of reaction target area；In Practical Calculation
During be individually to carry out above-mentioned imaging for each single-shot data, finally the imaging of each single-shot data is tied
Fruit is overlapped process, you can obtain the high-resolution seismic exploration ripple imaging results of earthquake shot gather data.
A kind of high-resolution seismic exploration wave imaging method provided in an embodiment of the present invention is carried out with reference to specific embodiment
Explanation：
By D seismic data processing, a kind of high-resolution seismic exploration wave imaging method based on multiple signal classification is illustrated
Application effect in infrastructure exploration.
(1) earthquake shot gather data and migration velocity file (i.e. above-mentioned migration velocity data) are read in；
(2) seimic travel time information is calculated by migration velocity file；
(3) tradition Kirchhoff imaging results, as shown in Figure 6；
(4) the high-resolution seismic exploration ripple imaging results of multiple signal classification are as shown in Figure 7；
Contrast Fig. 2 and Fig. 3, it should be apparent that at section vertical position 3.5 seconds, based on multiple signal classification
High-resolution seismic exploration ripple imaging section (Fig. 3) with more preferable seismic wave resolution capabilities, as shown in round frame in Fig. 6 and Fig. 7, and
And tomography is shown preferably, arrow pointed location in such as Fig. 6 and Fig. 7.In the imaging of deep seismic weak signal, tradition is compared
Kirchhoff imaging methods, the high-resolution seismic exploration wave imaging method based on multiple signal classification still have good advantage,
As shown in square frame position in Fig. 6 and Fig. 7.
In seismic imaging in a kind of high-resolution seismic exploration wave imaging method provided in an embodiment of the present invention, with prior art
It is good and cannot solve that steady phase point determines and seismic waveform chooses that aspect does not exist for the low image quality of seismic data of noise
Problem is compared, and which adopts multiple signal classification method to calculate the common imaging point road for obtaining of Kirchhoff migration imagings
Process, can accurately judge steady phase point and thus choose seismic waveform such that it is able to realize the estimation of effective seismic signal, improve
Seismic imaging resolution ratio；In its gathered data from low signal-to-noise ratio, high-resolution geological information is excavated, be that coal is surveyed
The conformation identification visited in exploitation provides more accurate seismic data, has with reducing input cost is reduced on construction risk
Significance.
The embodiment of the present invention additionally provides a kind of high-resolution seismic exploration ripple imaging device, and described device is used to perform above-mentioned one
High-resolution seismic exploration wave imaging method is planted, with reference to Fig. 8, described device includes：
First acquisition module 11, the corresponding skew speed of the single-shot data and single-shot data for obtaining pending imaging
Degrees of data；Wherein, single-shot data include the reflection of subterranean strata interface or the multiple seismic waves for reflecting, and seismic wave carries geology
Information；
Seimic travel time information computational module 12, for according to the corresponding wave detector locus of single-shot data and skew speed
Degrees of data, calculates the ground of any one imaging point position of shot position Jing subsurface imagings space of single-shot data to geophone station position
Information when seismic wave is walked；Wherein, the corresponding wave detector of single-shot data is for multiple, and each wave detector is arranged at corresponding geophone station
On position；
Common imaging gather sets up module 13, for the position constituted according to wave detector locus and seimic travel time information
Point identification is put, from the corresponding amplitude sampling point of chosen position point identification in single-shot data, the corresponding common imaging of single-shot data is obtained
Point road collection；
Processing module 14 is calculated, for carrying out calculating process to common imaging gather according to multiple signal classification method, is obtained
To the corresponding high-resolution seismic exploration ripple imaging results of single-shot data.
Further, with reference to Fig. 9, in above-mentioned high-resolution seismic exploration ripple imaging device, calculating processing module 14 includes：
Scan process unit 141, for the horizontal direction along amplitude sampling point, to each imaging in common imaging gather
The corresponding amplitude sampling point in point position is scanned process, obtains amplitude sampling point in the corresponding laterally main energy band of horizontal direction
Scope；Wherein, laterally the corresponding amplitude sampling point sequence of main energy band scope carries geological information；
Seismic wave size judging unit 142, for the longitudinal direction along amplitude sampling point, to horizontal main energy band scope
Each interior amplitude sampling point carries out seismic wave size judgement, obtains the longitudinal direction effectively earthquake corresponding to horizontal main energy band scope
Ripple magnitude range；
Imaging array extraction unit 143, for according to horizontal main energy band scope and its effective seismic wave in corresponding longitudinal direction
The imaging array scope of magnitude range composition, extracts the corresponding imaging array of imaging array scope from common imaging gather；
Estimate processing unit 144, for useful signal estimation process is carried out to imaging array, obtain single-shot data corresponding
High-resolution seismic exploration ripple imaging results.
Further, with reference to Fig. 9, in above-mentioned high-resolution seismic exploration ripple imaging device, scan process unit 141 includes：
First determination subelement 1411, for the horizontal direction along amplitude sampling point, determines every in common imaging gather
Transverse energy band scope of the corresponding amplitude sampling point in individual imaging point position in horizontal direction；
Summation process subelement 1412, for carrying out to all amplitude sampling points in the range of each transverse energy band
Summation process, obtains the amplitude sampling point superposition value with scope corresponding to each transverse energy；
Comparing subunit 1413, for comparing amplitude sampling point superposition value of each transverse energy with scope；
Subelement 1414 is chosen, the corresponding transverse energy band scope of the amplitude sampling point superposition value for choosing maximum is horizontal stroke
To main energy band scope.
Further, with reference to Figure 10, in above-mentioned high-resolution seismic exploration ripple imaging device, seismic wave size judging unit 142
Including：
Zero point and change in polarity judgment sub-unit 1421, for the longitudinal direction along amplitude sampling point, to horizontal main energy
Each amplitude sampling point in the range of amount band carries out zero point and change in polarity judges；
Second determination subelement 1422, for existing detecting each the amplitude sampling point in the range of horizontal main energy band
During zero point and the positive and negative extreme value including zero point on longitudinal direction, the corresponding longitudinal magnitude range of positive and negative extreme value is determined for corresponding
The longitudinal direction of amplitude sampling point effectively seismic wave magnitude range.
Further, with reference to Figure 10, in above-mentioned high-resolution seismic exploration ripple imaging device, estimate that processing unit 144 includes：
Resolution process subelement 1441, for the side of the longitudinal direction amount of predetermined number is chosen according to the line direction in imaging array
Method, carries out resolution process to imaging array, obtains multiple imaging submatrixs；Wherein, at least each two is continuously imaged submatrix
Longitudinally measure including identical；
First computation subunit 1442, for calculate respectively each imaging submatrix covariance matrix, obtain corresponding to
The sub- covariance matrix of each imaging submatrix；
Second computation subunit 1443, for carrying out average computation process to many sub- covariance matrixes for obtaining, obtains
Average covariance matrices；
3rd computation subunit 1444, for characteristic vector R of average covariance matrices is calculated by below equationaV=λ
v；Wherein, RaFor average covariance matrices, characteristic vectors of the v for average covariance matrices, features of the λ for average covariance matrices
Value；
Singular value decomposition processes subelement 1445, for average covariance matrices are carried out with singular value point by below equation
Solution processObtain earthquake effective signal subspace matrix and noise subspace matrix；
4th computation subunit 1446, for characteristic vector and noise subspace matrix according to average covariance matrices,
Calculating includes the pseudo- spectrum of the multiple signal classification of seismic imaging result.
Further, above-mentioned high-resolution seismic exploration ripple imaging device, also includes：
Second acquisition module, for obtaining the high-resolution seismic exploration ripple imaging results of multiple single-shot data；
The high-resolution seismic exploration ripple imaging results of multiple single-shot data are overlapped process, are obtained by overlap-add procedure module
The corresponding high-resolution seismic exploration ripple imaging results of shot gather data；Wherein, shot gather data includes multiple single-shot data.
In seismic imaging in a kind of high-resolution seismic exploration ripple imaging device provided in an embodiment of the present invention, with prior art
It is good and cannot solve that steady phase point determines and seismic waveform chooses that aspect does not exist for the low image quality of seismic data of noise
Problem is compared, and which adopts multiple signal classification method to calculate the common imaging point road for obtaining of Kirchhoff migration imagings
Process, can accurately judge steady phase point and thus choose seismic waveform such that it is able to realize the estimation of effective seismic signal, improve
Seismic imaging resolution ratio；In its gathered data from low signal-to-noise ratio, high-resolution geological information is excavated, be that coal is surveyed
The conformation identification visited in exploitation provides more accurate seismic data, has with reducing input cost is reduced on construction risk
Significance.
What the embodiment of the present invention was provided carried out a kind of computer program of high-resolution seismic exploration ripple imaging, including depositing
The computer-readable recording medium of program code is stored up, the instruction that described program code includes can be used to perform previous methods enforcement
Method described in example, implements and can be found in embodiment of the method, will not be described here.
Those skilled in the art can be understood that, for convenience and simplicity of description, the system of foregoing description,
The specific work process of device and unit, may be referred to the corresponding process in preceding method embodiment, will not be described here.
In several embodiments provided herein, it should be understood that disclosed system, apparatus and method, can be with
Realize by another way.Device embodiment described above is only schematically, for example division of the unit,
It is only a kind of division of logic function, there can be other dividing mode when actually realizing, and for example, multiple units or component can
To combine or be desirably integrated into another system, or some features can be ignored, or not perform.It is another, it is shown or beg for
By coupling each other direct-coupling or communication connection can be by the indirect of some communication interfaces, device or unit
Coupling is communicated to connect, and can be electrical, mechanical or other forms.
The unit as separating component explanation can be or may not be it is physically separate, it is aobvious as unit
The part for showing can be or may not be physical location, you can local to be located at one, or can also be distributed to multiple
On NE.Some or all of unit therein can be selected according to the actual needs to realize the mesh of this embodiment scheme
's.
In addition, each functional unit in each embodiment of the invention can be integrated in a processing unit, it is also possible to
It is that unit is individually physically present, it is also possible to which two or more units are integrated in a unit.
If the function is realized using in the form of SFU software functional unit and as independent production marketing or when using, can be with
It is stored in a computer read/write memory medium.Based on such understanding, technical scheme is substantially in other words
The part contributed to prior art or the part of the technical scheme can be embodied in the form of software product, the meter
Calculation machine software product is stored in a storage medium, is used including some instructions so that a computer equipment (can be individual
People's computer, server, or network equipment etc.) perform all or part of step of each embodiment methods described of the invention.
And aforesaid storage medium includes：USB flash disk, portable hard drive, read-only storage (ROM, Read-Only Memory), arbitrary access are deposited
Reservoir (RAM, Random Access Memory), magnetic disc or CD etc. are various can be with the medium of store program codes.
The above, the only specific embodiment of the present invention, but protection scope of the present invention is not limited thereto, any
Those familiar with the art the invention discloses technical scope in, change or replacement can be readily occurred in, should all be contained
Cover within protection scope of the present invention.Therefore, protection scope of the present invention should be defined by the scope of the claims.
Claims (10)
1. a kind of high-resolution seismic exploration wave imaging method, it is characterised in that include：
Obtain the single-shot data and the corresponding migration velocity data of the single-shot data of pending imaging；Wherein, the list
Big gun data include the reflection of subterranean strata interface or the multiple seismic waves for reflecting, and the seismic wave carries geological information；
According to the corresponding wave detector locus of the single-shot data and the migration velocity data, the single-shot data are calculated
Seimic travel time information of any one the imaging point position of shot position Jing subsurface imagings space to geophone station position；Wherein, institute
The corresponding wave detector of single-shot data is stated for multiple, and each described wave detector is arranged on the corresponding geophone station position；
According to the position point identification that the wave detector locus and the seimic travel time information are constituted, from the single-shot data
It is middle to choose the corresponding amplitude sampling point of the position point identification, obtain the corresponding common imaging gather of the single-shot data；
Calculating process is carried out to the common imaging gather according to multiple signal classification method, the single-shot data is obtained corresponding
High-resolution seismic exploration ripple imaging results；Wherein, the common imaging gather is carried out at calculating according to multiple signal classification method
Reason includes：Imaging array in the common imaging gather is extracted according to multiple signal classification method, the imaging array is entered
The estimation of row useful signal is processed；The imaging array carries the geological information.
2. method according to claim 1, it is characterised in that it is described according to multiple signal classification method to the common imaging
Point road collection carries out calculating process, and obtaining the corresponding high-resolution seismic exploration ripple imaging results of the single-shot data includes：
Along the horizontal direction of the amplitude sampling point, to each corresponding amplitude in imaging point position in the common imaging gather
Value sampling point is scanned process, obtains the amplitude sampling point in the corresponding laterally main energy band scope of the horizontal direction；Its
In, the corresponding amplitude sampling point sequence of the laterally main energy band scope carries the geological information；
Along the longitudinal direction of the amplitude sampling point, each the amplitude sampling point in the range of the laterally main energy band is carried out
The judgement of seismic wave size, obtains the longitudinal direction effectively seismic wave magnitude range corresponding to the laterally main energy band scope；
According to the imaging square of the effective seismic wave magnitude range composition of the laterally main energy band scope and its corresponding longitudinal direction
Battle array scope, extracts the corresponding imaging array of the imaging array scope from the common imaging gather；
Useful signal estimation process is carried out to the imaging array, obtain the corresponding high-resolution seismic exploration ripple of the single-shot data into
As result.
3. method according to claim 2, it is characterised in that the horizontal direction along the amplitude sampling point is right
In the common imaging gather, the corresponding amplitude sampling point in each imaging point position is scanned process, obtains the amplitude sample
Point includes in the corresponding laterally main energy band scope of the horizontal direction：
Along the horizontal direction of the amplitude sampling point, determine in the common imaging gather each imaging point position is corresponding and shake
Transverse energy band scope of the amplitude sampling point in the horizontal direction；
Summation process is carried out to all amplitude sampling points in the range of transverse energy band each described, is obtained corresponding to each
Amplitude sampling point superposition value of the individual transverse energy with scope；
Relatively amplitude sampling point superposition value of each described transverse energy with scope, chooses maximum amplitude sampling point superposition value
Corresponding transverse energy band scope is the laterally main energy band scope.
4. method according to claim 3, it is characterised in that the longitudinal direction along the amplitude sampling point is right
Each amplitude sampling point in the range of the laterally main energy band carries out seismic wave size judgement, obtains corresponding to the horizontal master
The longitudinal direction effectively seismic wave magnitude range of energy band scope includes：
Along the longitudinal direction of the amplitude sampling point, each the amplitude sampling point in the range of the laterally main energy band is carried out
Zero point and change in polarity judge；
In each the amplitude sampling point zero point in the longitudinal direction detected in the range of the laterally main energy band and
Including the zero point positive and negative extreme value when, determine the corresponding longitudinal magnitude range of the positive and negative extreme value be the corresponding amplitude
The longitudinal direction of sampling point effectively seismic wave magnitude range.
5. method according to claim 4, it is characterised in that described that the imaging array is carried out at useful signal estimation
Reason, obtaining the corresponding high-resolution seismic exploration ripple imaging results of the single-shot data includes：
The method that the longitudinal direction amount of predetermined number is chosen according to the line direction in the imaging array, is carried out to the imaging array point
Solution process, obtains multiple imaging submatrixs；Wherein, at least continuous imaging submatrix of each two includes identical longitudinal direction
Amount；
The covariance matrix of each imaging submatrix is calculated respectively, is obtained corresponding to the son association that submatrix is imaged each described
Variance matrix；
Multiple described sub- covariance matrix to obtaining carries out average computation process, obtains average covariance matrices；
Characteristic vector R of average covariance matrices is calculated by below equationaV=λ v；Wherein, RaFor average covariance matrices, v
For the characteristic vector of average covariance matrices, characteristic values of the λ for average covariance matrices；
Singular value decomposition process is carried out to the average covariance matrices by below equation
To earthquake effective signal subspace matrix and noise subspace matrix；
According to the characteristic vector and the noise subspace matrix of the average covariance matrices, calculating include the seismic wave into
As the pseudo- spectrum of the multiple signal classification of result.
6. the method according to claim 1-5 any one, it is characterised in that methods described also includes：
Obtain the high-resolution seismic exploration ripple imaging results of multiple single-shot data；
Process is overlapped to the high-resolution seismic exploration ripple imaging results of multiple single-shot data, shot gather data is obtained corresponding
High-resolution seismic exploration ripple imaging results；Wherein, the shot gather data includes multiple single-shot data.
7. a kind of high-resolution seismic exploration ripple imaging device, it is characterised in that described device includes：
First acquisition module, for obtaining the single-shot data and the corresponding migration velocity of the single-shot data of pending imaging
Data；Wherein, the single-shot data include the reflection of subterranean strata interface or the multiple seismic waves for reflecting, and the seismic wave is carried
There is geological information；
Seimic travel time information computational module, for according to the corresponding wave detector locus of the single-shot data and the skew
Speed data, calculates any one imaging point position of shot position Jing subsurface imagings space of the single-shot data to geophone station position
The seimic travel time information put；Wherein, the corresponding wave detector of the single-shot data is for multiple, and each described wave detector is respectively provided with
On the corresponding geophone station position；
Common imaging gather sets up module, for what is constituted according to the wave detector locus and the seimic travel time information
Position point identification, chooses the corresponding amplitude sampling point of the position point identification from the single-shot data, obtains the single-shot number
According to corresponding common imaging gather；
Processing module is calculated, for carrying out calculating process to the common imaging gather according to multiple signal classification method, is obtained
The corresponding high-resolution seismic exploration ripple imaging results of the single-shot data；Wherein, according to multiple signal classification method to it is described altogether into
Picture point road collection carries out calculating process to be included：Imaging square in the common imaging gather is extracted according to multiple signal classification method
Battle array, carries out useful signal estimation process to the imaging array；The imaging array carries the geological information.
8. device according to claim 7, it is characterised in that the calculating processing module includes：
Scan process unit, for the horizontal direction along the amplitude sampling point, in the common imaging gather each into
The corresponding amplitude sampling point of image point position is scanned process, obtains the amplitude sampling point in the corresponding horizontal stroke of the horizontal direction
To main energy band scope；Wherein, the corresponding amplitude sampling point sequence of the laterally main energy band scope carries the geology letter
Breath；
Seismic wave size judging unit, for the longitudinal direction along the amplitude sampling point, to the laterally main energy band model
Each amplitude sampling point in enclosing carries out seismic wave size judgement, obtains having corresponding to the longitudinal direction of the laterally main energy band scope
Effect seismic wave magnitude range；
Imaging array extraction unit, for according to laterally main energy band scope and its effective seismic wave in the corresponding longitudinal direction
The imaging array scope of magnitude range composition, extracts the corresponding imaging of the imaging array scope from the common imaging gather
Matrix；
Estimate processing unit, for useful signal estimation process is carried out to the imaging array, obtain the single-shot data correspondence
High-resolution seismic exploration ripple imaging results.
9. device according to claim 8, it is characterised in that the scan process unit includes：
First determination subelement, for the horizontal direction along the amplitude sampling point, determines every in the common imaging gather
Transverse energy band scope of the corresponding amplitude sampling point in individual imaging point position in the horizontal direction；
Summation process subelement, for suing for peace to all amplitude sampling points in the range of transverse energy band each described
Process, obtain the amplitude sampling point superposition value with scope corresponding to transverse energy each described；
Comparing subunit, for comparing amplitude sampling point superposition value of each described transverse energy with scope；
Subelement is chosen, the corresponding transverse energy band scope of the amplitude sampling point superposition value for choosing maximum is the horizontal master
Energy band scope.
10. device according to claim 9, it is characterised in that the seismic wave size judging unit includes：
Zero point and change in polarity judgment sub-unit, for the longitudinal direction along the amplitude sampling point, to the laterally main energy
Each amplitude sampling point in the range of amount band carries out zero point and change in polarity judges；
Second determination subelement, for detecting each the amplitude sampling point in the range of the laterally main energy band described
During zero point and the positive and negative extreme value including the zero point on longitudinal direction, the corresponding longitudinal magnitude range of the positive and negative extreme value is determined
For the longitudinal direction effectively seismic wave magnitude range of the corresponding amplitude sampling point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610074937.2A CN105652317B (en) | 2016-02-02 | 2016-02-02 | High-resolution seismic imaging method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610074937.2A CN105652317B (en) | 2016-02-02 | 2016-02-02 | High-resolution seismic imaging method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105652317A CN105652317A (en) | 2016-06-08 |
| CN105652317B true CN105652317B (en) | 2017-03-22 |
Family
ID=56488272
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610074937.2A Active CN105652317B (en) | 2016-02-02 | 2016-02-02 | High-resolution seismic imaging method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105652317B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107817523B (en) * | 2017-10-30 | 2018-10-09 | 中国矿业大学(北京) | The analysis method and device of diffracted wave migration velocity |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7388808B2 (en) * | 2004-09-23 | 2008-06-17 | Pgs Americas, Inc. | Method for depth migrating seismic data using pre-stack time migration, demigration, and post-stack depth migration |
| CN102162858A (en) * | 2010-12-06 | 2011-08-24 | 中国海洋石油总公司 | Method for analyzing dynamic correction velocity by using asymmetric travel time |
-
2016
- 2016-02-02 CN CN201610074937.2A patent/CN105652317B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105652317A (en) | 2016-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106772583B (en) | A kind of earthquake diffracted wave separation method and device | |
| CN104570125B (en) | Method for improving imaging speed model precision by using well data | |
| CN107505651B (en) | Seismic first break and back wave combine slope chromatography imaging method | |
| CN104237940B (en) | Diffracted wave imaging method and device based on dynamic characteristics | |
| CN102937721B (en) | Limited frequency tomography method for utilizing preliminary wave travel time | |
| US6574566B2 (en) | Automated feature identification in data displays | |
| CN106094032B (en) | A kind of method for building formation velocity model | |
| CN107817523B (en) | The analysis method and device of diffracted wave migration velocity | |
| CN103645503B (en) | A kind of three-dimensional time territory illumination analysis and vibration amplitude compensation method | |
| CN110308484A (en) | A kind of chromatography conversion method and system based on deep learning intelligent screening first arrival | |
| CN104268412B (en) | Angle gather ray chromatography migration velocity analysis method and device | |
| CN106443765A (en) | Urban engineering seismic detection comprehensive imaging method based on multi-component observation system | |
| CN1625699A (en) | A method for shallow water flow detection | |
| CN106772592B (en) | Diffracted wave focuses the analysis method and device of energy | |
| CN106772593B (en) | The imaging method and device of diffracted wave | |
| CN107218033A (en) | A kind of inversion method for the method and formation parameter for recognizing attitude of stratum | |
| CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
| CN106959467A (en) | Seimic wave velocity inversion method and device | |
| Marciniak et al. | Integrated geophysical imaging of a mountain landslide–a case study from the Outer Carpathians, Poland | |
| CN112305591A (en) | Tunnel advanced geological prediction method and computer readable storage medium | |
| CN102565852B (en) | Angle domain pre-stack offset data processing method aiming to detect oil-gas-bearing property of reservoir | |
| Huang et al. | Bayesian seismic refraction inversion for critical zone science and near‐surface applications | |
| CN104199088B (en) | Incident angle gather extraction method and system | |
| CN101609163B (en) | Multi-scale seismic data joint imaging method based on fluctuation theory | |
| Lellouch et al. | Microseismic analysis over a single horizontal distributed acoustic sensing fiber using guided waves |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |