CN107807391A - Seismic data processing method, device, electronic equipment and computer storage medium - Google Patents
Seismic data processing method, device, electronic equipment and computer storage medium Download PDFInfo
- Publication number
- CN107807391A CN107807391A CN201710820589.3A CN201710820589A CN107807391A CN 107807391 A CN107807391 A CN 107807391A CN 201710820589 A CN201710820589 A CN 201710820589A CN 107807391 A CN107807391 A CN 107807391A
- Authority
- CN
- China
- Prior art keywords
- geological data
- data
- wave
- numerical value
- velocity
- 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.)
- Granted
Links
- 238000003860 storage Methods 0.000 title claims abstract description 27
- 238000003672 processing method Methods 0.000 title abstract description 3
- 238000012545 processing Methods 0.000 claims abstract description 110
- 238000000034 method Methods 0.000 claims abstract description 107
- 238000012937 correction Methods 0.000 claims abstract description 76
- 230000001629 suppression Effects 0.000 claims description 60
- 230000008569 process Effects 0.000 claims description 54
- 238000006243 chemical reaction Methods 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 7
- 241001269238 Data Species 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 238000013506 data mapping Methods 0.000 claims description 2
- 238000013341 scale-up Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 11
- 238000001914 filtration Methods 0.000 abstract description 4
- 229910052704 radon Inorganic materials 0.000 description 17
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 17
- 230000009471 action Effects 0.000 description 13
- 230000006870 function Effects 0.000 description 13
- 238000010586 diagram Methods 0.000 description 12
- 238000000926 separation method Methods 0.000 description 6
- 238000013508 migration Methods 0.000 description 5
- 230000005012 migration Effects 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 238000003825 pressing Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011017 operating method Methods 0.000 description 2
- 238000000844 transformation Methods 0.000 description 2
- XDVMCVGTDUKDHL-UHFFFAOYSA-N [amino(2-azaniumylethylsulfanyl)methylidene]azanium;dibromide Chemical compound Br.Br.NCCSC(N)=N XDVMCVGTDUKDHL-UHFFFAOYSA-N 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/52—Move-out correction
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 aims to provide a seismic data processing method, a seismic data processing device, an electronic device and a computer storage medium, wherein the method comprises the following steps: acquiring first seismic data, wherein the first seismic data comprise primary waves and multiple waves; performing dynamic correction processing on the seismic data according to the specified speed to obtain second seismic data; wherein the designated velocity comprises a numerical value of the primary velocity, a numerical value of the multiple velocity, or a numerical value between a numerical value of the primary velocity and a numerical value of the multiple velocity; performing preset processing on the second seismic data to obtain third seismic data; the presetting comprises reducing the offset distance in the second seismic data according to a specified rule, and/or increasing the time in the second seismic data according to a specified rule. The implementation mode of the invention can better distinguish the primary wave from the multiple wave, thereby obtaining better multiple wave filtering effect.
Description
Technical field
The present invention relates to field of seismic exploration, more particularly to a kind of seismic data processing technique, device, electronic equipment and meter
Calculation machine storage medium.
Background technology
Generally, more subwaves are considered as a kind of coherent noise of severe jamming primary wave in seismic data process.
More subwaves are superimposed with primary wave, it is difficult to separate, this not only lowers the signal to noise ratio of seismic data, and have a strong impact on
The authenticity and reliability of seismic imaging, adverse effect is brought to follow-up explanation work.
The method of multiple suppression is pressed typically according to the difference of more subwaves and primary wave using the method for filtering at present
System, such as Radon conversion, frequency-wavenumber (F-K frequency-wavenumber) conversion and predictive deconvolution etc..Due to multiple
The difference of the time difference of ripple and primary wave possible very little, especially interbed multiple and primary wave is very small, influences filter effect.
The content of the invention
The purpose of this specification embodiment is to provide a kind of seismic data processing technique, device, electronic equipment and calculating
Machine storage medium, being capable of effective multiple suppression.
To achieve the above object, this specification embodiment provides a kind of seismic data processing technique, and methods described includes：
The first geological data is obtained, wherein, first geological data includes primary wave and more subwaves；According to command speed to described
Shake data and move correction process, obtain the second geological data；Wherein, the numerical value of the command speed an including wave velocity,
The repeatedly numerical value between the numerical value of wave velocity or the numerical value of wave velocity and the numerical value of multiple wave velocity；Second ground
Shake data include multiple dimensions, and the multiple dimension includes geophone offset and time；Default place is carried out to second geological data
Reason, obtains the 3rd geological data；The default processing includes, according to reducing in second geological data specified rule
Geophone offset, and/or, the time in second geological data is increased according to specified rule.In second geological data
In, the primary wave understands with more subwaves to be separated.
To achieve the above object, this specification embodiment also provides a kind of seismic data process device, described device bag
Include：Acquiring unit, the acquiring unit are used to obtain the first geological data, wherein, first geological data includes primary wave
With more subwaves；First processing units, first processing units are used to move correction process to the geological data according to command speed, obtained
Obtain the second geological data；Wherein, the command speed includes the numerical value of a wave velocity, the numerical value of the repeatedly wave velocity
Or the numerical value between the numerical value of wave velocity and the numerical value of multiple wave velocity；Second geological data includes multiple dimensions,
The multiple dimension includes geophone offset and time；Second processing unit, second processing unit are used for second geological data
Default processing is carried out, obtains the 3rd geological data；The default processing includes, and reduces second earthquake number according to specified rule
The geophone offset in, and/or, the time in second geological data is increased according to specified rule.
To achieve the above object, this specification embodiment also provides a kind of electronic equipment, including：Input equipment, storage
Device, processor；The input equipment be used for obtain the first geological data, wherein, first geological data include primary wave with
More subwaves；The processor is used to move correction process to the geological data according to command speed, obtains the second geological data；Its
In, the command speed includes the numerical value of a wave velocity, the repeatedly numerical value of wave velocity or the number of a wave velocity
Numerical value between value and the repeatedly numerical value of wave velocity；Second geological data includes multiple dimensions, and the multiple dimension includes
Geophone offset and time；Default processing is carried out to second geological data, obtains the 3rd geological data；The default processing bag
Include, the geophone offset reduced according to specified rule in second geological data, and/or, according to increasing specified rule
The time in second geological data；The memory is used to access data, and the data include the input or are acquired
Geological data, the intermediate data and result data of processor processing.
To achieve the above object, this specification embodiment also provides a kind of computer-readable storage medium, and the computer is deposited
Storage media is stored with computer program instructions, is realized when the computer program instructions are performed：The first geological data is obtained,
Wherein, first geological data includes primary wave and more subwaves；Correction process is moved to the geological data according to command speed,
Obtain the second geological data；Wherein, the command speed includes the numerical value of a wave velocity, the number of the repeatedly wave velocity
Numerical value between the numerical value of the numerical value and multiple wave velocity of value or wave velocity；Second geological data includes multiple dimensions
Degree, the multiple dimension include geophone offset and time；Default processing is carried out to second geological data, obtains the 3rd earthquake number
According to；The default processing includes, the geophone offset reduced according to specified rule in second geological data, and/or, according to
Specified rule increases the time in second geological data.
To achieve the above object, this specification embodiment also provides a kind of seismic data processing technique, methods described bag
Include：The first geological data is obtained, wherein, first geological data includes primary wave and more subwaves；According to the primary wave
The velocity of wave of velocity of wave and more subwaves, determines command speed；The command speed between the primary wave velocity of wave with it is described more
Between the velocity of wave of subwave；Correction process is moved to the geological data according to the command speed, obtains the second geological data.Institute
State in the second geological data, primary wave and more subwave substantial separations.
To achieve the above object, this specification embodiment also provides a kind of seismic data process device, described device bag
Include：Acquiring unit, the acquiring unit are used to obtain the first geological data, wherein, first geological data includes primary wave
With more subwaves；First determining unit, first determining unit are used for velocity of wave and more subwaves according to the primary wave
Velocity of wave, determines command speed, and the command speed is between the velocity of wave of the primary wave and the velocity of wave of more subwaves；Second
Determining unit, first determining unit are used to move correction process to the geological data according to the command speed, obtain the
Two geological datas.
The technical scheme provided from above the application embodiment, the application embodiment is mainly by determining to specify
Speed, the command speed is between the primary wave velocity of wave and more subwave velocities of wave, to enter action correction, substantial separation
The primary wave and more subwaves, then by reducing the geophone offset in second geological data, and/or, described in increase
Time parameter in second geological data, the primary wave and more subwaves are further separated, be follow-up multiple suppression
Method provides condition.Above way make use of on the basis of routinely dynamic correction, using command speed to the primary wave overcorrect,
To more subwave undercorrections, the primary wave and more subwaves are separated, and pass through the change of geophone offset or time, further separation
More subwaves, it is possible to achieve to the more preferable pressing result of more subwaves.
Brief description of the drawings
, below will be to embodiment party in order to illustrate more clearly of this specification embodiment or technical scheme of the prior art
The required accompanying drawing used is briefly described in formula or description of the prior art, it should be apparent that, drawings in the following description are only
Only it is some embodiments described in this specification, for those of ordinary skill in the art, is not paying creative labor
On the premise of dynamic property, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 is a kind of flow chart for seismic data process that this specification provides；
Fig. 2 is the geological data schematic diagram containing two groups of primary waves and more subwaves that this specification provides；
Fig. 3 is the geological data schematic diagram entered to geological data after action correction that this specification provides；
Fig. 4 is the primary wave in specified coordinate system and more subwave schematic diagrames that this specification provides；
Fig. 5 is another primary wave in specified coordinate system and more subwave schematic diagrames that this specification provides；
Fig. 6 is the geological data schematic diagram after the reaction correction that this specification provides；
Fig. 7 is the geological data schematic diagram after another reaction correction that this specification provides；
Fig. 8 is CRP gather (CRP common reflection in the actual conditions that this specification provides
Point) phenogram；
Fig. 9 is the geological data sign entered in the actual conditions that this specification provides to geological data after action correction
Figure；
Figure 10 is the geological data phenogram after multiple suppression in the actual conditions that this specification provides；
Figure 11 is the phenogram of geological data after reaction correction in the actual conditions that this specification provides；
Figure 12 is the normal-moveout spectrum calculated data before multiple wave pressure system that this specification provides；
Figure 13 is the normal-moveout spectrum calculated data after multiple wave pressure system that this specification provides；
Figure 14 is the phenogram for directly carrying out migration before stack result that this specification provides；
Figure 15 is the phenogram of migration before stack result after the multiple suppression that this specification provides；
Figure 16 is a kind of electronic equipment schematic diagram that this specification provides.
Embodiment
Below in conjunction with the accompanying drawing in this specification embodiment, the technical scheme in this specification embodiment is carried out
Clearly and completely describing, it is clear that described embodiment is only a part of embodiment of this specification, rather than all
Embodiment.Based on the embodiment in this specification, those of ordinary skill in the art are not before creative work is made
The every other embodiment obtained is put, should all belong to the scope of the application protection.
Refer to Fig. 1.A kind of flow chart for seismic data processing technique that this specification provides.Methods described can include
Following steps.
In the present embodiment, the object for performing the seismic data process can be the electronics for having logical operation function
Equipment.The electronic equipment can be server and client side.The client can be desktop computer, tablet personal computer, notes
This computer, work station etc..Certainly, client is not limited to the above-mentioned electronic equipment with certain entity, and it can also be operation
Software in above-mentioned electronic equipment.A kind of program software formed by program development is can also be, the program software can be with
Run in above-mentioned electronic equipment.
Step S10：The first geological data is obtained, wherein, first geological data includes primary wave and more subwaves.
In the present embodiment, the geological data refers to the original earthquake gathered in seismic prospecting by various sensors
Obtained in the geological data that data or intermediate treatment are crossed, such as the data directly obtained from seismic prospecting.It can pass through
Geological data is handled and explained, may infer that the property and form of subterranean strata.The geological data includes primary wave
With more subwaves.The primary wave represents that, from focus, when reaching receiving point (geophone station), there occurs one between low interface
The ripple of secondary reflection, the repeatedly wave table shows from focus, when reaching receiving point (geophone station), between subsurface interface there occurs
Multiple reflections or the ripple of refraction.Wherein, the speed of the primary wave is more than the speed of more subwaves.
In the present embodiment, first geological data is obtained, can be set by the input of input equipment, the input
Standby can be keyboard, mouse, voice-input device etc.；It can also be set by the input of external storage equipment, the external storage
Standby can be USB flash disk, mechanical hard disk etc.；It can also be received by network, such as internet, LAN；Can also be by reading this
Ground data etc..
Step S12：Correction process is moved to the geological data according to command speed, obtains the second geological data；Wherein, institute
State the numerical value, the repeatedly numerical value of wave velocity or the numerical value of a wave velocity of command speed an including wave velocity with it is more
Numerical value between the numerical value of secondary wave velocity；Second geological data includes multiple dimensions, and the multiple dimension includes geophone offset
And the time.
In the present embodiment, the command speed can be the velocity of wave of the primary wave or the velocity of wave of more subwaves or
Speed between the two.Preferably, the command speed can be the primary wave velocity of wave and more subwave velocities of wave are averaged
Speed.When command speed is the average speed of the primary wave velocity of wave and more subwave velocities of wave, after moving correction process, more have
Beneficial to the separation of the primary wave and more subwaves.
In the present embodiment, referring to Fig. 2, first geological data is in first coordinate system, time-big gun
Inspection is away from coordinate system.First geological data includes, the reception time of the primary wave, the primary wave energy group, described one
The race inspection of subwave away from, more subwaves receive the time, the repeatedly wave energy group, more subwaves race inspection away from etc..It is described
It is dynamic correction be in the first geological data time carry out TEC time error correction, so second geological data still in the time-
Geophone offset coordinate system, second geological data include multiple dimensions, and the multiple dimension includes geophone offset and time.
In the present embodiment, it is described it is dynamic correction refer on the different each road of geophone offset come from same interface, same point
Back wave arrival time (receive time), the echo time (receiving the time) being corrected at CMP.I.e. conventional dynamic correction
The reflection interval of different geophone offsets is corrected to the process of the reflection interval of zero shot-geophone distance, it is described according to the command speed to institute
It is that the command speed is replaced into the speed in the dynamic correction of the routine to enter action correction to state geological data and move correction process.
Step S14：Default processing is carried out to second geological data, obtains the 3rd geological data；The default processing
Including, the geophone offset reduced according to specified rule in second geological data, and/or, institute is increased according to specified rule
State the time in the second geological data.
In the present embodiment, the default processing refers to the processing for amplifying the primary wave and more subwave differences.
The separating degree of the primary wave and more subwaves can be increased, good basis is provided for later stage multiple suppression.
In the present embodiment, the default processing includes the institute reduced according to specified rule in second geological data
Geophone offset is stated, and/or, the time in second geological data is increased according to specified rule.The specified rule can be with
Make zoom by a certain percentage.Such as the geophone offset in the second geological data is multiplied by 0.3 or 0.5 simultaneously, naturally it is also possible to
Geophone offset or time are reduced by exponential manner, still, what is put is too big, and curvature range is just very big, and one can increase amount of calculation,
As long as the two easily produces alias so be amplified to the degree can that can be separated, it is not to be the bigger the better, preferably presses
It is amplified and reduces according to certain proportion.It is described to reduce in second geological data in a specific embodiment
The geophone offset is the geophone offset reduced by a certain percentage in second geological data.In another specific embodiment party
In formula, it is described in second geological data to increase by a certain percentage to increase the time in second geological data
Time.In another particular embodiment of the invention, the geophone offset reduced according to specified rule in second geological data,
With the time in second geological data is increased according to specified rule and carried out simultaneously.Pass through the default processing, increase
The primary wave and the absolute value of the curvature of more subwaves, further separate the primary wave and more subwaves.Tie below
Close Fig. 3 to illustrate, in Fig. 3, a wave profile raises up, and curvature is negative, and the repeatedly wave profile is downward, and curvature is just.
Race by reducing corresponding points in Fig. 3 coordinate systems is examined away from then curve compresses to the origin of coordinates, and the curvature figure of the primary wave is special
The curvature figure characteristic parameter absolute value for levying parameter and more subwaves increases, further separate the primary wave with it is described repeatedly
Ripple.Similarly, by increasing the time in Fig. 3 coordinate systems, the curvature figure characteristic parameter of the primary wave and more subwaves
Curvature figure characteristic parameter absolute value can also increase, or the mode of two kinds of default processing is used in combination, and be attained by increasing
The effect of the curvature figure characteristic parameter and the absolute value of the curvature figure characteristic parameter of more subwaves of the big primary wave.
In a specific Sample Scenario, the seismic data processing technique of this specification embodiment offer, Ke Yitong
The software run in electronic equipment is crossed to realize.
In this Sample Scenario, the first geological data is obtained, wherein, first geological data is including primary wave and repeatedly
Ripple；It can also can also can be received by the input of input equipment by the input of external storage equipment by network,
The first geological data can be obtained by reading local data etc..Correction process is moved to the geological data according to command speed,
Obtain the second geological data；Wherein, the command speed includes the numerical value of a wave velocity, the number of the repeatedly wave velocity
Numerical value between the numerical value of the numerical value and multiple wave velocity of value or wave velocity.Dynamic correction is on the basis of routinely dynamic correction
The dynamic correction for replacing conventional NMO velocity to implement with command speed.Default processing is carried out to second geological data, is obtained
To the 3rd geological data；The default processing includes, and the big gun reduced according to specified rule in second geological data is examined
Away from, and/or, the time in second geological data is increased according to specified rule.For example, by second card data
The time that geophone offset is all multiplied by the 0.5, or second geological data is all multiplied by 2.The default processing can further separate described
The curvature of primary wave and more subwaves, is easy to filtering out for more subwaves.Above way make use of the primary wave with it is described
More subwaves have different speed, and make use of dynamic correction is original to carry out initial gross separation, further through default processing, i.e., using described
Primary wave is with more subwaves in the property of two dimensions of coordinate system, point of the processing amplification primary wave and more subwaves
From degree.
In one embodiment, correction process is moved to the geological data according to the command speed, obtains described the
Two geological datas, including：According to the command speed, obtain and move correction time difference data at different geophone offsets corresponding to geophone station；
The hourage that geophone station at different geophone offsets receives and the dynamic correction time difference data are corrected, obtain second ground
Shake data；Wherein, the hourage is to geophone station required time after seismic wave produces.
In the present embodiment, the hourage refers to back wave arrival time, and the reception time of wave detector.
In the present embodiment, the dynamic correction time difference that the dynamic correction of conventional levels obtains obtains according to equation below：
Wherein, the Δ t is the dynamic correction time difference, and x is geophone offset, and V is the spread speed of focus in media as well, i.e. statement pair
The primary wave velocity of wave or corresponding more subwave velocities of wave answered, t0When for the reception time of NMO corrected point being the reflection at zero shot-geophone distance
Between, t is the reception time at corresponding geophone offset.It is described that correction process is moved to the geological data according to the command speed and incited somebody to action
The command speed replaces the V.Dynamic correction in other cases in the same way or principle, is no longer repeated herein.
When the command speed is less than primary wave velocity of wave, during more than more subwave velocities of wave, the primary wave Δ t dynamic correction time differences can increase,
More subwave Δ t dynamic correction time differences can reduce, and primary wave is overcorrect, and more subwaves are undercorrections.With reference to Fig. 2
Illustrated with Fig. 3.Described Fig. 2 represents the model data containing two groups of primary waves Yu more subwaves, in each group, the primary wave
Approached with more subwaves, be unfavorable for multiple suppression.The data in Fig. 2 are entered at correction of taking action according to the command speed
Reason, output result is as shown in figure 3, the primary wave is overcorrect, and more subwaves are undercorrections.When command speed is once
During ripple velocity of wave, the primary wave correction is normal, more subwave undercorrections；When command speed is more subwave velocities of wave, described one
Subwave overcorrect, Method And Principle are no longer repeated herein with the dynamic correction process for using midrange speed.
In one embodiment, the geophone offset reduced according to the specified rule in second geological data,
And/or the time in specified rule increase second geological data includes：Scaled described second
The geophone offset in geological data, and/or, scale up the time in second geological data.
In the present embodiment, the geophone offset in scaled second geological data, and/or, press
Ratio increases the time in second geological data, such as the geophone offset in the second geological data is multiplied by into 0.3 or 0.5
Deng or the time is multiplied by the modes of 2 or 3 grades, be amplified curvature, can also be in scaled second geological data
The geophone offset and the time scaled up in second geological data while the mode of progress, are amplified song
Rate.This method is unlikely to put curvature figure characteristic parameter too greatly, avoids amount of calculation excessive and produces alias.
In one embodiment, the 3rd geological data at least conforms to fisrt feature, and the fisrt feature is described
The curvature figure characteristic parameter of the curvature figure characteristic parameter of primary wave and more subwaves is not at same positive number interval or negative
Number interval；After the 3rd geological data is obtained, methods described also includes：According to the fisrt feature, to the 3rd ground
Shake data and carry out multiple suppression processing, obtain the 3rd geological data after multiple suppression.
Refering to Fig. 3, in the present embodiment, when the command speed is used between the primary wave and more subwaves
Midrange speed enters timing of taking action, the primary wave overcorrect, more subwave undercorrections, and a wave profile raises up, institute
It is downward to state multiple wave profile.The primary wave curvature is negative, and the repeatedly curve rate is positive number.Similarly, when the specified speed
When degree uses the speed of the primary wave, more subwave undercorrections, the curvature figure characteristic parameter of more subwaves is in just
Number interval, the curvature figure characteristic parameter of the primary wave are not at positive number interval；When the command speed is using described multiple
During the speed of ripple, the primary wave overcorrect, the curvature figure characteristic parameter of the primary wave is in negative section, described multiple
The curvature figure characteristic parameter of ripple is not at negative section.Therefore the primary wave curvature figure characteristic parameter with it is described repeatedly
The curvature figure characteristic parameter of ripple is not at same positive number interval or negative section.
In the present embodiment, after the 3rd geological data is obtained, methods described also includes：It is special according to described first
Sign, multiple suppression processing is carried out to the 3rd geological data, obtains the 3rd geological data after multiple suppression.
In the present embodiment, it is not at by the curvature figure characteristic parameter of the primary wave and more subwaves same
Section, processing is filtered to the more subwaves wherein separated.In a specific embodiment, can be converted by Radon,
More subwaves are filtered in Radon domains.
In a Sample Scenario, obtain the 3rd geological data after multiple suppression and may comprise steps of.
Step S20：By the 3rd seismic data mapping to specified coordinate system, the 4th geological data is obtained；It is described to specify
Coordinate system includes curvature figure characteristic dimension and time dimension.In a specific embodiment, to the 3rd geological data
Parabolic Radon transform is carried out according to equation below.
M (τ, q)=Σ d (x, t=τ+qx2)
Wherein, m (τ, q) is the data in Radon domains (τ, q), and d (x, t) is the data of geophone offset-time-domain, is herein
For the 3rd geological data in first coordinate system, τ is the intercept time, and q is curvature parameters.Geophone offset-time-domain
Curvature of the lineups in Radon domains isWherein Δ t=tmax- τ, tmaxThe time for being lineups at maximum offset.
Illustrated with reference to Fig. 5, Fig. 5 is that the 3rd earthquake number of default processing acquisition is being carried out to second geological data
According to Radon conversion is carried out, as illustrated, the abscissa of the specified coordinate system represents curvature, ordinate represents the time.Shown one
Subwave curvature figure characteristic parameter is negative, and after Radon is converted, the energy group of the primary wave is shown more in negative section
Subwave curvature figure characteristic parameter is just, and after Radon is converted, the energy groups of more subwaves is easy to pair in positive number interval
The filtering process of more subwaves.
Step S22：Processing is filtered to the 4th geological data, obtains the 4th earthquake number after multiple suppression
According to.Referring to Fig. 5, in a specific embodiment, the primary wave energy group is in negative section, the repeatedly wave energy group
In positive number interval, the energy group for aligning number interval is filtered, and obtains more pure primary wave energy group.
Step S24：Processing is penetrated into the 4th geological data reflection after the multiple suppression, obtains the multiple suppression
The 3rd geological data afterwards.Specifically, anti-Radon conversion is carried out according to equation below to the 4th geological data.
D (x, t)=∑ m (q, τ=t-qx2)
Wherein, m (τ, q) is the data in Radon domains (τ, q), and d (x, t) is the data of geophone offset-time-domain, i.e., described
Second data volume in first coordinate system, τ are the intercept time, and q is curvature parameters.The lineups of geophone offset-time-domain exist
The curvature in Radon domains isWherein Δ t=tmax- τ, tmaxThe time for being lineups at maximum offset.
In one embodiment, after the 3rd geological data after obtaining multiple suppression, in addition to：Refer to according to second
Set pattern then, the geophone offset in amplifying after the 3rd geological data after the multiple suppression, and/or, specified according to second
Rule, the time in reducing after the 3rd geological data after the multiple suppression, obtain second after multiple suppression
Geological data；Wherein, second specified rule is the inverse transformation corresponding to the specified rule.For example, work as the default place
The mode of reason is that the geophone offset of first data volume is reduced into original 1/3rd, then now, the compacting is multiple
The geophone offset of the second data volume after ripple amplifies three times, recovers to geophone offset.Similarly, if the mode of the default processing is to put
The big time, then the specified rule is the time corresponding to diminution.By the processing of specified rule, second data volume can be made
In data recovery to actual conditions under geophone offset and the time.
In one embodiment, also include after the second geological data after obtaining the multiple suppression：To described
The second geological data after multiple suppression carries out reaction correction process, obtains the first geological data after multiple suppression.Its
In, the reaction correction process corresponds to the dynamic correction process, i.e., according to the command speed by the reflection of different geophone offsets
Time from the reflection interval of zero shot-geophone distance recovers to the process of different geophone offsets each self-corresponding reflection interval, to facilitate the later stage pair
The processing of geological data.
A specific Sample Scenario of this specification is illustrated with reference to Fig. 2 to Fig. 7.This specification embodiment
The seismic data processing technique of offer, it can be realized by running on the software in electronic equipment.The seismic data processing technique
So that primary wave is more preferably distinguished with more subwaves, so as to obtain more preferable more subwave filter effects, facilitate subsequent job.
In this implement scene, the first geological data is obtained, wherein, first geological data is including primary wave and repeatedly
Ripple.Refering to Fig. 2, the geological data schematic diagram of two groups of primary waves of acquisition and more subwaves.Obtain first geological data, Ke Yitong
The input of input equipment is crossed, the input equipment can be keyboard, mouse, voice-input device etc.；It can also be deposited by external
The input of equipment is stored up, the external storage equipment can be USB flash disk, mechanical hard disk etc.；It can also be received by network, such as because of spy
Net, LAN；Can also be by reading local data etc..
In this implement scene, correction process is moved to the geological data according to command speed, obtains the second geological data；
Wherein, the numerical value of the command speed an including wave velocity, the repeatedly numerical value of wave velocity or wave velocity
Numerical value between numerical value and the numerical value of multiple wave velocity；Second geological data includes multiple dimensions, the multiple dimension bag
Include geophone offset and time.The dynamic correction refers to the back wave from same interface, same point on the different each road of geophone offset
Arrival time (receives the time), the echo time (receiving the time) being corrected at CMP.It is i.e. conventional to move correction by different big guns
Examine away from reflection interval be corrected to zero shot-geophone distance reflection interval process, it is described according to the command speed to the earthquake number
It is that the command speed is replaced into the speed in the dynamic correction of the routine to enter action correction according to dynamic correction process.Refering to Fig. 3, pass through
The command speed determined enters action to geological data described in Fig. 2 and corrected, herein, the command speed for it is described once
The average speed of wave velocity and the repeatedly wave velocity, moves the primary wave overcorrect after correction, curvature is negative, more subwaves
Undercorrection, curvature is just.
In this implement scene, default processing is carried out to second geological data, obtains the 3rd geological data；It is described pre-
If processing includes, the geophone offset reduced according to specified rule in second geological data, and/or, according to specified rule
Increase the time in second geological data.Refering to Fig. 5, data in Fig. 3 are carried out with default processing, the default processing
To reduce geophone offset and/or increase time, and carry out high-resolution Radon conversion.Specifically, simultaneously by the second geological data
Geophone offset is multiplied by 0.3 or 0.5.By the default processing, increase the primary wave and the absolute value of the curvature of more subwaves,
Further separate the primary wave and more subwaves.Race by reducing corresponding points in Fig. 3 coordinate systems is examined away from then curve is to seat
Origin compression is marked, the curvature figure characteristic parameter absolute value of the curvature figure characteristic parameter of the primary wave and more subwaves increases
Greatly, the primary wave and more subwaves are further separated.Similarly, by increasing the time in Fig. 3 coordinate systems, the primary wave
The curvature figure characteristic parameter absolute value of curvature figure characteristic parameter and more subwaves can also increase, it is or pre- by two kinds
If the mode of processing is used in combination, the song of the curvature figure characteristic parameter Yu more subwaves that increase the primary wave is attained by
The effect of the absolute value of rate graphic feature parameter.Comparison diagram 5 and Fig. 4, the primary wave energy group and the repeatedly wave energy group
Separating effect is more preferable.Refering to Fig. 7, processing is filtered to more subwaves in Fig. 5 Radon domains, then carries out Radon inverse transformations,
According to the specified rule corresponding to the default processing, amplify geophone offset and/or reduce the time recover geophone offset and/or when
Between, reaction correction process has then been carried out, has obtained the geological data after multiple suppression.Comparison diagram 7 and Fig. 6, multiple wave pressure system
Effect is better than Fig. 6 multiple ripple pressing result.
With reference to Fig. 8 to Figure 14, treatment effect of the specific this specification embodiment to real data.Refering to Fig. 8, figure
8 represent specific CRP gather (CRP common reflection point), and certainly, present embodiment can also transport
For trace gathers such as common midpoint gathers (CMP common middle point).Enter action correction to the CRP trace gathers shown in Fig. 8,
Obtain the trace gather shown in Fig. 9, it can be seen that more subwave information of undercorrection be present.Default place is carried out to second geological data
Reason, obtains the 3rd geological data.Multiple suppression processing is carried out to the 3rd geological data, obtains the after multiple suppression
Three geological datas.The 3rd geological data after the multiple suppression is handled, obtains the second ground after multiple suppression
Data are shaken, obtain the phenogram after the multiple suppression shown in Figure 10.It is that reaction school is carried out to the data in Figure 10 refering to Figure 11
Just, obtain and suppressed geological data phenogram after more subwaves.Comparison diagram 11 and Fig. 8, it can be seen that more subwaves are fine
Compacting.Refering to Figure 12, Figure 13, Figure 12 corresponds to the normal-moveout spectrum shown in Fig. 8, and Figure 13 corresponds to Figure 11 normal-moveout spectrum, in figure
Dotted line represents the pickup leads of the primary wave velocity of wave, and pickup leads lower left represents multiple wave energy, and Figure 13 is in contrast to Figure 12 lower-lefts
The multiple wave energy of side is significantly filtered out.Refering to Figure 14, Figure 15, Figure 14 is to the number in the data before multiple wave pressure system i.e. Fig. 8
According to, directly carry out pre-stack time migration result, Figure 15 is to the data in the data after multiple wave pressure system i.e. Figure 11, is folded
The result of preceding time migration, out of ellipse and rectangle frame as can be seen that after multiple wave pressure system, the lineups in migrated section connect
Continuous property improves, and resolution ratio is also improved.For example, the lineups in rectangle frame, before multiple wave pressure system, two lineups have one
Determine the interference of degree；After multiple wave pressure system, the two lineups have obtained preferably separating.
A kind of seismic data process device is additionally provided in this specification embodiment, as described in following embodiment.
It is similar to a kind of seismic data processing technique to solve the principle of problem due to a kind of seismic data process device, therefore a kind of earthquake
The implementation of data processing equipment may refer to a kind of implementation of seismic data processing technique, repeats part and repeats no more.Following institute
Use, term " unit " or " module " can realize the combination of the software and/or hardware of predetermined function.Although implement below
Device described by example is preferably realized with software, but hardware, or the realization of the combination of software and hardware is also possible
And be contemplated.The device can specifically include：Acquiring unit, the first determining unit, first processing units, second processing list
Member.The structure is specifically described below.
Acquiring unit, the acquiring unit are used to obtain the first geological data, wherein, first geological data includes one
Subwave and more subwaves.
In the present embodiment, the geological data includes a wave velocity and multiple wave velocity.Wherein, the earthquake number
According to can be obtained from the data that seismic prospecting obtains；The primary wave is represented from focus, when reaching receiving point, low
There occurs the ripple of primary event between interface, the repeatedly wave table shows from focus, when reaching receiving point, subsurface interface it
Between there occurs multiple reflections or the ripple of refraction.
First processing units, first processing units are used to move correction process to the geological data according to command speed, obtained
Obtain the second geological data；Wherein, the command speed includes the numerical value of a wave velocity, the numerical value of the repeatedly wave velocity
Or the numerical value between the numerical value of wave velocity and the numerical value of multiple wave velocity；Second geological data includes multiple dimensions,
The multiple dimension includes geophone offset and time.
In the present embodiment, conventional dynamic correction is the reflection that the reflection interval of different geophone offsets is corrected to zero shot-geophone distance
The process of time, it is described according to the command speed to the geological data move correction process be by the command speed replace institute
State the speed in conventional dynamic correction.
Second processing unit, second processing unit are used to carry out second geological data default processing, obtain the 3rd
Geological data；The default processing includes, the geophone offset reduced according to specified rule in second geological data, and/
Or, the time in second geological data is increased according to specified rule.
In a present embodiment, after second processing unit, described device also includes：3rd processing unit, is used for
Obtain the 3rd geological data after multiple suppression.3rd geological data at least conforms to fisrt feature, the fisrt feature
Same positive number area is not at for the curvature figure characteristic parameter of the primary wave and the curvature figure characteristic parameter of more subwaves
Between or negative section；After the 3rd geological data is obtained, methods described also includes：According to the fisrt feature, to described
3rd geological data carries out multiple suppression processing, obtains the 3rd geological data after multiple suppression.
In a present embodiment, after the 3rd processing unit, in addition to：Fourth processing unit.The fourth process
Unit is according to the second specified rule, the geophone offset in amplifying after the 3rd geological data after the multiple suppression, and/
Or, according to the second specified rule, the time in reducing after the 3rd geological data after the multiple suppression, suppressed
The second geological data after more subwaves；Wherein, second specified rule is the inverse transformation corresponding to the specified rule.
In a present embodiment, after fourth processing unit, in addition to：5th processing unit；5th processing
Unit carries out reaction correction process to the second geological data after the multiple suppression, obtains the first ground after multiple suppression
Shake data.
Refer to Figure 16.A kind of electronic equipment is additionally provided in this specification embodiment, including：Input equipment, storage
Device, processor；The input equipment be used for obtain the first geological data, wherein, first geological data include primary wave with
More subwaves；The processor is used to move correction process to the geological data according to command speed, obtains the second geological data；Its
In, the command speed includes the numerical value of a wave velocity, the repeatedly numerical value of wave velocity or the number of a wave velocity
Numerical value between value and the repeatedly numerical value of wave velocity；Second geological data includes multiple dimensions, and the multiple dimension includes
Geophone offset and time；Default processing is carried out to second geological data, obtains the 3rd geological data；The default processing bag
Include, the geophone offset reduced according to specified rule in second geological data, and/or, according to increasing specified rule
The time in second geological data；The memory is used to access data, and the data include the input or are acquired
Geological data, the intermediate data and result data of processor processing.
The input equipment is one of main device for entering row information exchange between user and computer system.The input
Equipment includes keyboard, mouse, camera, scanner, light pen, writing input board, speech input device etc.；Input equipment is used for handle
Programs of these numbers of initial data and processing are input in computer.The input equipment can also obtain receive other modules,
The data that unit, equipment transmit.
The processor can be implemented in any suitable manner.For example, processor can take such as microprocessor or
Processor and storage can be by the calculating of the computer readable program code (such as software or firmware) of (micro-) computing device
Machine computer-readable recording medium, gate, switch, application specific integrated circuit (Application Specific Integrated Circuit,
ASIC), form of programmable logic controller (PLC) and embedded microcontroller etc..
It is used to protect stored memory device in the memory modern information technologies.The memory includes plurality of layers
It is secondary, in digital display circuit, as long as can preserve binary data can be memory；In integrated circuits, one does not have in fact
The circuit with store function of thing form is also memory, such as RAM, FIFO；In systems, there is the storage of physical form
Equipment is also memory, such as memory bar, TF card.
In the present embodiment, the function and effect of electronic equipment specific implementation, can be compareed with other embodiment
Explain, will not be repeated here.
A kind of computer-readable storage medium is additionally provided in this specification embodiment, the computer-readable storage medium is stored with
Computer program instructions, realized when the computer program instructions are performed：The first geological data is obtained, wherein, described
One geological data includes primary wave and more subwaves；Correction process is moved to the geological data according to command speed, obtains the second ground
Shake data；Wherein, the command speed includes the numerical value of a wave velocity, the numerical value or primary wave of the repeatedly wave velocity
Numerical value between the numerical value of speed and the numerical value of multiple wave velocity；Second geological data includes multiple dimensions, the multiple
Dimension includes geophone offset and time；Default processing is carried out to second geological data, obtains the 3rd geological data；It is described default
Processing includes, the geophone offset reduced according to specified rule in second geological data, and/or, increased according to specified rule
The time in big second geological data.
In the present embodiment, the memory includes but is not limited to random access memory (Random Access
Memory, RAM), read-only storage (Read-Only Memory, ROM), caching (Cache), hard disk (Hard Disk
Drive, HDD) or storage card (Memory Card).The memory can be used for storing computer program instructions.Network leads to
Believe that unit can be set according to standard as defined in communication protocol, for carrying out the interface of network connection communication.
In the present embodiment, the function and effect of the programmed instruction specific implementation of computer-readable storage medium storage, can
To compare explanation with other embodiment, will not be repeated here.
Another seismic data processing technique that this specification also provides, methods described include：The first geological data is obtained,
Wherein, first geological data includes primary wave and more subwaves；According to the velocity of wave of the primary wave and the ripple of more subwaves
Speed, determine command speed；The command speed is between the velocity of wave of the primary wave and the velocity of wave of more subwaves；According to institute
State command speed and correction process is moved to the geological data, obtain the second geological data.
In one embodiment, second geological data at least conforms to fisrt feature, and the fisrt feature is described
The curvature figure characteristic parameter of the curvature figure characteristic parameter of primary wave and more subwaves is not at same positive number interval or negative
Number interval；After second geological data is obtained, methods described also includes：According to the fisrt feature, to second ground
Shake data and carry out multiple suppression processing, obtain the second geological data after multiple suppression.
In one embodiment, also include after the second geological data after obtaining the multiple suppression：To described
The second geological data after multiple suppression carries out reaction correction process, obtains the first geological data after multiple suppression.
One embodiment of this specification is illustrated with reference to Fig. 2, Fig. 3, Fig. 4, Fig. 6.Refering to Fig. 2, two groups are obtained
The geological data schematic diagram of primary wave and more subwaves.Refering to Fig. 3, by the command speed of determination to earthquake described in Fig. 2
Data enter action correction, the primary wave overcorrect, and curvature is negative, more subwave undercorrections, and curvature is just.Refering to Fig. 4, figure
4 show the result that data in Fig. 3 are carried out with high-resolution Radon conversion, and as described in Figure, more subwaves are separated with primary wave.
Refering to Fig. 6, processing is filtered to more subwaves in Fig. 4 Radon domains, then carries out Radon inverse transformations, and carried out reaction school
Positive processing, obtains the geological data after multiple suppression, there it can be seen that more subwaves are significantly suppressed.
The first seismic data processing technique provided with reference to this specification, another geological data that this specification provides
Processing method needs not move through default processing and corresponding specified rule, and method is more succinct, but the pressing result of more subwaves
Not as the geological data after the multiple suppression by default processing acquisition, more subwaves of near migration range still have residual.
This specification embodiment propose seismic data processing technique it is most essential the characteristics of be, in one embodiment,
Enter action using the velocity of wave midrange speed of velocity of wave and the more subwaves of the primary wave to correct, separate the primary wave with repeatedly
Ripple.In another embodiment, by reducing the geophone offset parameter, and/or, increase the time parameter, increase is once
The separation degree of ripple and more subwaves.In a kind of preferred embodiment, at the same using the primary wave velocity of wave with it is described repeatedly
The velocity of wave midrange speed of ripple enters action correction and by reducing the geophone offset parameter, and/or, increase the time parameter, it is right
The primary wave is separated with more subwaves, and the compacting for more subwaves in geological data provides improved method.
This specification also provides a kind of seismic data process device, and described device includes：Acquiring unit, the acquiring unit
For obtaining the first geological data, wherein, first geological data includes primary wave and more subwaves；First determining unit, institute
The velocity of wave that the first determining unit is used for velocity of wave and more subwaves according to the primary wave is stated, determines command speed, the finger
Constant speed degree is between the velocity of wave of the primary wave and the velocity of wave of more subwaves；Second determining unit, described first determines list
Member is used to move correction process to the geological data according to the command speed, obtains the second geological data.
Although mention a kind of seismic data processing technique, device, electronic equipment or computer storage in teachings herein to be situated between
Matter, still, the application are not limited to be professional standard or the situation described by embodiment etc., some professional standards or
Embodiment amended slightly can also realize above-mentioned implementation on the practice processes described using self-defined mode or embodiment
Example it is identical, equivalent or close or deformation after it is anticipated that implementation result.Using the data acquisition after these modifications or deformation, place
The embodiment of reason, output, judgment mode etc., still it may belong within the scope of the optional embodiment of the application.
Although this application provides the method operating procedure as described in embodiment or flow chart, based on conventional or noninvasive
The means for the property made can include more or less operating procedures.The step of being enumerated in embodiment order is only numerous steps
A kind of mode in execution sequence, does not represent unique execution sequence., can when device or client production in practice performs
So that according to embodiment, either method shown in the drawings order is performed or parallel performed (such as at parallel processor or multithreading
The environment of reason, even distributed data processing environment).Term " comprising ", "comprising" or its any other variant are intended to contain
Lid nonexcludability includes, so that process, method, product or equipment including a series of elements not only will including those
Element, but also the other element including being not expressly set out, or it is this process, method, product or equipment also to include
Intrinsic key element.In the absence of more restrictions, be not precluded from the process including the key element, method, product or
Other identical or equivalent elements in person's equipment also be present.
Device that above-described embodiment illustrates or module etc., it can specifically be realized by computer chip or entity, or by having
There is the product of certain function to realize.For convenience of description, it is divided into various modules during description apparatus above with function to retouch respectively
State.Certainly, the function of each module can be realized in same or multiple softwares and/or hardware when implementing the application,
The module for realizing same function can be realized by the combination of multiple submodule etc..Device embodiment described above is only
Schematically, for example, the division of the module, only a kind of division of logic function, can there is other draw when actually realizing
The mode of dividing, such as multiple module or components can combine or be desirably integrated into another system, or some features can be ignored,
Or do not perform.
It is also known in the art that in addition to realizing controller in a manner of pure computer readable program code, it is complete
Entirely can by by method and step carry out programming in logic come controller with gate, switch, application specific integrated circuit, may be programmed
The form of logic controller and embedded microcontroller etc. realizes identical function.Therefore this controller is considered one kind
Hardware component, and what its inside was included is used to realize that the device of various functions can also to be considered as the structure in hardware component.Or
Person even, not only can be able to will be the software module of implementation method but also can be hardware for realizing that the device of various functions is considered as
Structure in part.
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure, class etc..The application can also be put into practice in a distributed computing environment, in these DCEs,
By performing task by communication network and connected remote processing devices.In a distributed computing environment, program module can
With in the local and remote computer-readable storage medium including storage device.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
Realized by the mode of software plus required general hardware platform.Based on such understanding, the technical scheme essence of the application
On the part that is contributed in other words to prior art can be embodied in the form of software product, the computer software product
It can be stored in storage medium, such as ROM/RAM, magnetic disc, CD, including some instructions are causing a computer equipment
(can be personal computer, mobile terminal, server, either network equipment etc.) perform each embodiment of the application or implementation
Method described in some parts of example.
Each embodiment in this specification is described by the way of progressive, same or analogous portion between each embodiment
Divide mutually referring to what each embodiment stressed is the difference with other embodiment.The application can be used for crowd
In more general or special purpose computing system environments or configuration.Such as：Personal computer, server computer, handheld device or
Portable set, laptop device, multicomputer system, the system based on microprocessor, set top box, programmable electronics are set
Standby, network PC, minicom, the DCE etc. of mainframe computer including any of the above system or equipment.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and
Change is without departing from spirit herein, it is desirable to which appended claim includes these deformations and changed without departing from the application.
Claims (14)
1. a kind of seismic data processing technique, it is characterised in that methods described includes：
The first geological data is obtained, wherein, first geological data includes primary wave and more subwaves；
Correction process is moved to the geological data according to command speed, obtains the second geological data；Wherein, the command speed bag
Include the numerical value, the repeatedly numerical value of wave velocity or the numerical value of a wave velocity and the number of multiple wave velocity of a wave velocity
Numerical value between value；Second geological data includes multiple dimensions, and the multiple dimension includes geophone offset and time；
Default processing is carried out to second geological data, obtains the 3rd geological data；The default processing includes, according to specified
The geophone offset that rule is reduced in second geological data, and/or, second geological data is increased according to specified rule
In the time.
2. the method as described in claim 1, it is characterised in that moved according to the command speed to the geological data at correction
Reason, second geological data is obtained, including：
According to the command speed, obtain and move correction time difference data at different geophone offsets corresponding to geophone station；
The hourage that geophone station at different geophone offsets receives and the dynamic correction time difference data are corrected, acquisition described the
Two geological datas；Wherein, the hourage is to geophone station required time after seismic wave produces.
3. the method as described in claim 1, it is characterised in that reduced according to the specified rule in second geological data
The geophone offset, and/or, increasing the time in second geological data according to the specified rule includes：By than
The geophone offset that example is reduced in second geological data, and/or, scale up described in second geological data
Time.
4. the method as described in claim 1, it is characterised in that the 3rd geological data at least conforms to fisrt feature, described
Fisrt feature is not at together for the curvature figure characteristic parameter of the primary wave and the curvature figure characteristic parameter of more subwaves
One positive number interval or negative section；After the 3rd geological data is obtained, methods described also includes：It is special according to described first
Sign, multiple suppression processing is carried out to the 3rd geological data, obtains the 3rd geological data after multiple suppression.
5. method as claimed in claim 4, it is characterised in that multiple suppression processing is carried out to the 3rd geological data,
Obtaining the 3rd geological data after the multiple suppression includes：
By the 3rd seismic data mapping to specified coordinate system, the 4th geological data is obtained；The specified coordinate system includes song
Rate graphic feature dimension and time dimension；
Processing is filtered to the 4th geological data, obtains the 4th geological data after multiple suppression；
Processing is penetrated into the 4th geological data reflection after the multiple suppression, obtains the 3rd earthquake after the multiple suppression
Data.
6. method as claimed in claim 5, it is characterised in that after the 3rd geological data after obtaining multiple suppression, also
Including：
According to the second specified rule, the geophone offset in amplifying after the 3rd geological data after the multiple suppression, and/
Or, according to the second specified rule, the time in reducing after the 3rd geological data after the multiple suppression, suppressed
The second geological data after more subwaves；Wherein, second specified rule is the inverse transformation corresponding to the specified rule.
7. method as claimed in claim 6, it is characterised in that after the second geological data after obtaining the multiple suppression
Also include：
Reaction correction process is carried out to the second geological data after the multiple suppression, obtains the first ground after multiple suppression
Shake data.
8. a kind of seismic data process device, it is characterised in that described device includes：
Acquiring unit, the acquiring unit are used to obtain the first geological data, wherein, first geological data includes primary wave
With more subwaves；
First processing units, first processing units are used to move correction process to the geological data according to command speed, obtain the
Two geological datas；Wherein, the command speed includes the numerical value of a wave velocity, the numerical value or one of the repeatedly wave velocity
Numerical value between the numerical value of secondary wave velocity and the numerical value of multiple wave velocity；Second geological data includes multiple dimensions, described
Multiple dimensions include geophone offset and time；
Second processing unit, second processing unit are used to carry out second geological data default processing, obtain the 3rd earthquake
Data；The default processing includes, the geophone offset reduced according to specified rule in second geological data, and/or, root
Increase the time in second geological data according to specified rule.
9. a kind of electronic equipment, it is characterised in that including：Input equipment, memory, processor；
The input equipment is used to obtain the first geological data, wherein, first geological data includes primary wave and more subwaves；
The processor is used to move correction process to the geological data according to command speed, obtains the second geological data；Wherein,
The numerical value of the command speed an including wave velocity, the repeatedly numerical value of wave velocity or the numerical value of a wave velocity with
Numerical value between the numerical value of multiple wave velocity；Second geological data includes multiple dimensions, and the multiple dimension is examined including big gun
Away from and the time；Default processing is carried out to second geological data, obtains the 3rd geological data；The default processing includes, root
The geophone offset reduced according to specified rule in second geological data, and/or, second ground is increased according to specified rule
Shake the time in data；
The memory is used to access data, and the data include the input or the geological data being acquired, the processor
The intermediate data and result data of processing.
A kind of 10. computer-readable storage medium, it is characterised in that the computer-readable storage medium is stored with computer program instructions,
Realized when the computer program instructions are performed：
The first geological data is obtained, wherein, first geological data includes primary wave and more subwaves；
Correction process is moved to the geological data according to command speed, obtains the second geological data；Wherein, the command speed bag
Include the numerical value, the repeatedly numerical value of wave velocity or the numerical value of a wave velocity and the number of multiple wave velocity of a wave velocity
Numerical value between value；Second geological data includes multiple dimensions, and the multiple dimension includes geophone offset and time；
Default processing is carried out to second geological data, obtains the 3rd geological data；The default processing includes, according to specified
The geophone offset that rule is reduced in second geological data, and/or, second geological data is increased according to specified rule
In the time.
11. a kind of seismic data processing technique, it is characterised in that methods described includes：
The first geological data is obtained, wherein, first geological data includes primary wave and more subwaves；
According to the velocity of wave of the primary wave and the velocity of wave of more subwaves, command speed is determined；The command speed is between described
Between the velocity of wave of the velocity of wave of primary wave and more subwaves；
Correction process is moved to the geological data according to the command speed, obtains the second geological data.
12. method as claimed in claim 11, it is characterised in that second geological data at least conforms to fisrt feature, institute
The curvature figure characteristic parameter for stating curvature figure characteristic parameter and more subwaves that fisrt feature is the primary wave is not at
Same positive number interval or negative section；After second geological data is obtained, methods described also includes：It is special according to described first
Sign, multiple suppression processing is carried out to second geological data, obtains the second geological data after multiple suppression.
13. method as claimed in claim 12, it is characterised in that the second geological data after the multiple suppression is obtained
Also include afterwards：
Reaction correction process is carried out to the second geological data after the multiple suppression, obtains the first ground after multiple suppression
Shake data.
14. a kind of seismic data process device, it is characterised in that described device includes：
Acquiring unit, the acquiring unit are used to obtain the first geological data, wherein, first geological data includes primary wave
With more subwaves；
First determining unit, first determining unit are used for the velocity of wave of the velocity of wave and more subwaves according to the primary wave,
Command speed is determined, the command speed is between the velocity of wave of the primary wave and the velocity of wave of more subwaves；
Second determining unit, first determining unit are used to move at correction the geological data according to the command speed
Reason, obtain the second geological data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710820589.3A CN107807391B (en) | 2017-09-13 | 2017-09-13 | Seismic data processing method, device, electronic equipment and computer storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710820589.3A CN107807391B (en) | 2017-09-13 | 2017-09-13 | Seismic data processing method, device, electronic equipment and computer storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107807391A true CN107807391A (en) | 2018-03-16 |
| CN107807391B CN107807391B (en) | 2020-02-14 |
Family
ID=61592140
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710820589.3A Active CN107807391B (en) | 2017-09-13 | 2017-09-13 | Seismic data processing method, device, electronic equipment and computer storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107807391B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108051859A (en) * | 2018-01-09 | 2018-05-18 | 中国石油天然气股份有限公司 | Method and device for processing road set record and storage medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102636811A (en) * | 2012-04-10 | 2012-08-15 | 恒泰艾普石油天然气技术服务股份有限公司 | Eliminating method of multiple waves in bidimensional seismic data on sea |
| CN103149592A (en) * | 2013-03-07 | 2013-06-12 | 天津城市建设学院 | Method for separating variable offset vertical seismic profile (VSP) wave fields |
| CN104280777A (en) * | 2013-07-12 | 2015-01-14 | 中国石油天然气集团公司 | Method for suppressing interference of seismic data multiples on land |
| CN104914466A (en) * | 2015-06-26 | 2015-09-16 | 中国石油大学(华东) | Method for improving seismic data resolution |
-
2017
- 2017-09-13 CN CN201710820589.3A patent/CN107807391B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102636811A (en) * | 2012-04-10 | 2012-08-15 | 恒泰艾普石油天然气技术服务股份有限公司 | Eliminating method of multiple waves in bidimensional seismic data on sea |
| CN103149592A (en) * | 2013-03-07 | 2013-06-12 | 天津城市建设学院 | Method for separating variable offset vertical seismic profile (VSP) wave fields |
| CN104280777A (en) * | 2013-07-12 | 2015-01-14 | 中国石油天然气集团公司 | Method for suppressing interference of seismic data multiples on land |
| CN104914466A (en) * | 2015-06-26 | 2015-09-16 | 中国石油大学(华东) | Method for improving seismic data resolution |
Non-Patent Citations (1)
| Title |
|---|
| 孙成禹 等: "一种无拉伸畸变的动校正方法", 《石油物探》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108051859A (en) * | 2018-01-09 | 2018-05-18 | 中国石油天然气股份有限公司 | Method and device for processing road set record and storage medium |
| CN108051859B (en) * | 2018-01-09 | 2020-01-07 | 中国石油天然气股份有限公司 | Method and device for processing road set record and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107807391B (en) | 2020-02-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| SA519410263B1 (en) | Generating common image gather using wave-field separation | |
| US20180089533A1 (en) | Automated methods and systems for locating document subimages in images to facilitate extraction of information from the located document subimages | |
| CN109001813B (en) | Method, device and system for suppressing multiple waves | |
| CN105116442A (en) | Lithologic oil-gas reservoir weak-reflection seismic signal reconstruction method | |
| CN108564550B (en) | Image processing method and device and terminal equipment | |
| AU2014345427B2 (en) | Method and device for processing seismic signals | |
| CN109308465A (en) | Table line detecting method, apparatus, equipment and computer-readable medium | |
| CN114255337A (en) | Method and device for correcting document image, electronic equipment and storage medium | |
| CN103955499A (en) | Visual experience enhancing method based on instant calculation and dynamic tracing | |
| CN104375186B (en) | A kind of method and system for removing abnormal first arrival | |
| CN106019374B (en) | The tomograph imaging method of similarity is dissipated based on reflection groove wave frequency | |
| Kim et al. | A novel hardware design for SIFT generation with reduced memory requirement | |
| CN107807391A (en) | Seismic data processing method, device, electronic equipment and computer storage medium | |
| CN104781698A (en) | System and method for isolating signal in seismic data | |
| US8694687B2 (en) | Computing-system identifier using software extraction of manufacturing variability | |
| CN109633744B (en) | Method, device and equipment for extracting seismic wavelets and storage medium | |
| CN108693558B (en) | Seismic data processing method and device | |
| CN109801300B (en) | Method, device and equipment for extracting coordinates of checkerboard corner points and computer-readable storage medium | |
| WO2002021159A2 (en) | Method and system for attenuation of acquisition geometry imprint in seismic data | |
| CN114740527A (en) | First arrival take-off time pickup method and device | |
| CN114969089A (en) | Data obtaining method and device, electronic equipment and storage medium | |
| CN112767337A (en) | Cigarette case detection method and device | |
| CN107356973A (en) | A kind of seismic channel set data processing method and device | |
| CN114170589A (en) | Rock lithology identification method based on NAS, terminal equipment and storage medium | |
| CN108717708B (en) | Method and device for calculating edge slope of regular image |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |