CN106646615B - A kind of data processing method and device of surface wave frequency dispersion curve - Google Patents
A kind of data processing method and device of surface wave frequency dispersion curve Download PDFInfo
- Publication number
- CN106646615B CN106646615B CN201611244008.8A CN201611244008A CN106646615B CN 106646615 B CN106646615 B CN 106646615B CN 201611244008 A CN201611244008 A CN 201611244008A CN 106646615 B CN106646615 B CN 106646615B
- Authority
- CN
- China
- Prior art keywords
- seismic
- surface wave
- frequency
- phase difference
- phase
- 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
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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
Abstract
The application provides the data processing method and device of a kind of surface wave frequency dispersion curve.The described method includes: the seismic target earthquakes road in selected frequency domain seismic data to be processed；Determine the adjacent seismic channel quantity threshold and phase difference threshold in the seismic target earthquakes road；M adjacent seismic channels are chosen altogether in the seismic target earthquakes road；The M adjacent seismic channels and the seismic target earthquakes road are calculated separately in the phase difference of each reference frequency point；Determine that phase difference is greater than the road number of the seismic channel of the phase difference threshold on the reference frequency point；Determine corresponding reference frequency point whether be surface wave frequency dispersion curve cumulative frequency point.Using embodiment each in the application, the noise section in frequency-phase velocity frequency dispersion spectrum can effectively eliminate, to effectively improve the precision of the surface wave frequency dispersion curve picked up from the frequency-phase velocity frequency dispersion spectrum.
Description
Technical field
This application involves seismic data processing technology field, in particular to a kind of data processing method of surface wave frequency dispersion curve
And device.
Background technique
In seismic data processing technology, carrying out near surface structure inverting using surface wave is to carry out earth formation analysis to commonly use
Technology.Carrying out near surface inverting using surface wave mainly includes two processing stages, and first processing stage is to carry out surface wave dispersion
Analysis generates surface wave frequency dispersion curve；Second processing stage is anti-using the surface wave frequency dispersion curve progress near surface structure generated
It drills.Therefore, the precision of this inversion procedure mode depends on the precision of surface wave dispersion analysis.Obtain the surface wave dispersion of higher precision
Curve is the key that improve inversion result precision.
In the prior art, surface wave frequency dispersion curve in order to obtain, usually adds up to seismic data, and line frequency of going forward side by side dissipates place
Reason obtains frequency-phase velocity frequency dispersion spectrum, picks up surface wave frequency dispersion curve from the frequency-phase velocity frequency dispersion spectrum.But tired
During adding, due to other than linear surface wave, further including other nonlinear seismic waves in the seismic data, these ground
The phase of seismic wave and the phase difference of surface wave are larger, will result in the interference of opposite wave accumulation result in this way, cause frequency-phase velocity
Spending in frequency dispersion spectrum includes noise section, and then reduces the precision for picking up obtained surface wave frequency dispersion curve.
At least there are the following problems in the prior art: in the prior art, when adding up to seismic data, in seismic data
Contain nonlinear seismic wave, the nonlinear seismic wave and surface wave phase difference cause surface wave accumulation result away from larger
Interference causes in frequency-phase velocity frequency dispersion spectrum comprising noise section, and then causes to pick up from the frequency-phase velocity frequency dispersion spectrum
The surface wave frequency dispersion curve precision obtained is lower.
Summary of the invention
The purpose of the embodiment of the present application is to provide the data processing method and device of a kind of surface wave frequency dispersion curve, to eliminate frequency
Noise section in rate-phase velocity frequency dispersion spectrum, obtains the surface wave frequency dispersion curve of higher precision to pick up.
The embodiment of the present application provides a kind of data processing method of surface wave frequency dispersion curve and device is achieved in that
A kind of data processing method of surface wave frequency dispersion curve, which comprises
In the corresponding seismic channel of seismic data to be processed, selected target seismic channel；
Determine the adjacent seismic channel quantity threshold and phase difference threshold in the seismic target earthquakes road, the adjacent seismic channel number
Threshold value is equal to N, N >=2；
According to the phase difference threshold, determine the M adjacent seismic channels in seismic target earthquakes road and the seismic target earthquakes road every
The phase equalization of a Frequency point, N≤M≤2N；
According to determining as a result, determining in each Frequency point, with seismic target earthquakes road phase for the phase equalization
The road number of inconsistent adjacent seismic channel；
According to described in each Frequency point, the road number of inconsistent adjacent seismic channel with seismic target earthquakes road phase, really
It makes and may participate in cumulative Frequency point.
In preferred embodiment, the method also includes:
The corresponding seismic data of cumulative Frequency point will be may participate in and carry out the cumulative scattered processing of line frequency of going forward side by side, obtain surface wave dispersion
Curve.
It is described according to the phase difference threshold in preferred embodiment, determine the M in seismic target earthquakes road Yu the seismic target earthquakes road
The mode of phase equalization of a adjacent seismic channel in each Frequency point, comprising:
At each Frequency point, according to the phase difference threshold, each of described M neighboring track is determined respectively
The phase equalization of seismic channel and the seismic target earthquakes road.
It is described according to the phase difference threshold in preferred embodiment, with determining each of described M neighboring track respectively
Shake the mode of the phase equalization in road and the seismic target earthquakes road, comprising:
Calculate separately out the phase difference value of each seismic channel Yu the seismic target earthquakes road；
According to the phase difference value and phase difference threshold, determine each seismic channel in the M neighboring track with it is described
The phase equalization in seismic target earthquakes road.
It is described according to the phase difference value and phase difference threshold in preferred embodiment, it determines in the M neighboring track
The mode of the phase equalization of each seismic channel and the seismic target earthquakes road, comprising:
If the phase difference value is less than the phase difference threshold, it is determined that the result of phase equalization is that phase is consistent；
If the phase difference value is not less than the phase difference threshold, it is determined that the result of phase equalization is that phase is different
It causes.
In preferred embodiment, it is described according in each Frequency point, the inconsistent phase with seismic target earthquakes road phase
The road number of adjacent seismic channel determines the mode that may participate in cumulative Frequency point, comprising:
Determine the inconsistent road number threshold value of phase；
According to the road number and institute in each Frequency point and the inconsistent adjacent seismic channel of seismic target earthquakes road phase
The inconsistent road number threshold value of phase is stated, it is cumulative to determine whether each Frequency point may participate in, according to definitive result, counting can join
With cumulative Frequency point.
In preferred embodiment, the basis is in each Frequency point and the inconsistent adjacent earthquake of seismic target earthquakes road phase
The road number in road and the inconsistent road number threshold value of the phase, determine whether each Frequency point may participate in cumulative mode, packet
It includes:
It is different no more than the phase if the described and inconsistent adjacent seismic channel of seismic target earthquakes road phase road number
Cause road number threshold value, it is determined that corresponding Frequency point may participate in cumulative；
If it is inconsistent to be greater than the phase for the described and inconsistent adjacent seismic channel of seismic target earthquakes road phase road number
Road number threshold value, it is determined that corresponding Frequency point can not be accumulated.
In preferred embodiment, the adjacent seismic channel quantity threshold in the determination seismic target earthquakes road and phase difference threshold
Mode includes:
Surface wave dispersion processing is carried out to the seismic data to be processed, obtains frequency-phase velocity frequency dispersion spectrum；
According to the frequency-phase velocity frequency dispersion spectrum, determine the seismic target earthquakes road adjacent seismic channel quantity threshold and
Phase difference threshold.
A kind of data processing equipment of surface wave frequency dispersion curve, described device include:
Seismic channel chosen module, in the corresponding seismic channel of seismic data to be processed, selected target seismic channel；
Threshold determination module, for determining the adjacent seismic channel quantity threshold and phase difference threshold in the seismic target earthquakes road,
It is also used to determine the inconsistent road number threshold value of phase；
Phase equalization determining module, for determining seismic target earthquakes road with the target according to the phase difference threshold
Shake the phase equalization of the M adjacent seismic channels in each Frequency point in road；
Road number determining module, for determining as a result, determining in each Frequency point, with institute according to the phase equalization
State the road number of the inconsistent adjacent seismic channel of seismic target earthquakes road phase；
Frequency point determining module is used in each Frequency point according to, inconsistent with seismic target earthquakes road phase
The road number of adjacent seismic channel is determined to may participate in cumulative Frequency point.
In preferred embodiment, described device further include:
Surface wave frequency dispersion curve generation module add up simultaneously for will may participate in the corresponding seismic data of cumulative Frequency point
Frequency dispersion processing is carried out, surface wave frequency dispersion curve is obtained.
In preferred embodiment, the phase equalization determining module:
Phase difference value computing module, for calculating separately out the phase difference of each seismic channel Yu the seismic target earthquakes road
Value；
Phase equalization discrimination module, for determining that the M is a adjacent according to the phase difference value and phase difference threshold
The phase equalization of each seismic channel and the seismic target earthquakes road in road；
Recognition result output module is used for output phase consistency definitive result.
In preferred embodiment, the Frequency point determining module includes:
The inconsistent road number threshold determination module of phase, for determining the inconsistent road number threshold value of phase；
Can additive discrimination module, for according to described inconsistent in each Frequency point and seismic target earthquakes road phase
The road number of adjacent seismic channel and the inconsistent road number threshold value of the phase, it is cumulative to determine whether each Frequency point may participate in；
Can cumulative frequency point statistical module, for whether may participate in cumulative determination according to each Frequency point as a result,
It counts and may participate in cumulative Frequency point.
It, can be by determining earthquake using a kind of data processing method of surface wave frequency dispersion curve provided by the embodiments of the present application
The corresponding seismic channel of data determines the phase of adjacent seismic channel in seismic data not in the phase equalization threshold value of each Frequency point
The unanimously Frequency point of (phase difference is larger).In seismic data cumulative process, the consistent (phase of the phase for the adjacent seismic channel that only adds up
Potential difference is different small) the corresponding seismic data of Frequency point, can effectively eliminate the noise section in frequency-phase velocity frequency dispersion spectrum, from
And effectively improve the precision of the surface wave frequency dispersion curve picked up from the frequency-phase velocity frequency dispersion spectrum.Utilize the application reality
The data processing equipment for applying a kind of surface wave frequency dispersion curve of example offer, can be performed automatically the method, does not need implementation personnel
Excessive participation, so that it may the surface wave frequency dispersion curve for obtaining degree of precision effectively increases user experience.
Detailed description of the invention
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
The some embodiments recorded in application, for those of ordinary skill in the art, in the premise of not making the creative labor property
Under, it is also possible to obtain other drawings based on these drawings.
Fig. 1 is a kind of method flow of the data processing method of the surface wave frequency dispersion curve provided in the application one embodiment
Schematic diagram；
Fig. 2 is a kind of method stream of the data processing method of the surface wave frequency dispersion curve provided in another embodiment of the application
Journey schematic diagram；
Fig. 3 is a kind of modular structure of the data processing equipment of the surface wave frequency dispersion curve provided in the application one embodiment
Schematic diagram；
Fig. 4 is a kind of module knot of the data processing equipment of the surface wave frequency dispersion curve provided in another embodiment of the application
Structure schematic diagram.
Fig. 5 is the seismic data to be processed provided in one example of the application；
Fig. 6, Fig. 7, Fig. 8 are to carry out surface wave dispersion analysis to the seismic data to be processed in Fig. 5 in one example of the application
Obtained surface wave frequency-phase velocity spectrum.
Specific embodiment
The embodiment of the present application provides the data processing method and device of a kind of surface wave frequency dispersion curve.
In order to make those skilled in the art better understand the technical solutions in the application, below in conjunction with the application reality
The attached drawing in example is applied, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described implementation
Example is merely a part but not all of the embodiments of the present application.Based on the embodiment in the application, this field is common
The application protection all should belong in technical staff's every other embodiment obtained without creative efforts
Range.
Fig. 1 is a kind of a kind of method flow signal of embodiment of data processing method of herein described surface wave frequency dispersion curve
Figure.Although this application provides as the following examples or method operating procedure shown in the drawings or apparatus structure, based on conventional
It either may include more or less operating procedure or module list in the method or device without creative labor
Member.In the step of there is no necessary causalities in logicality or structure, the execution sequence of these steps or the module of device
Structure is not limited to the embodiment of the present application or execution shown in the drawings sequence or modular structure.The method or modular structure
Device or end product in practice according to embodiment or method shown in the drawings or modular structure in application, can carry out suitable
Sequence execute or it is parallel execute (such as parallel processor or multiple threads environment, even include distributed treatment reality
Apply environment).
A kind of method flow of the data processing method of the surface wave frequency dispersion curve provided in Fig. 1 the application one embodiment shows
Be intended to, specifically, as shown in Figure 1, the method may include:
S1: the seismic target earthquakes road in selected frequency domain seismic data to be processed.
The seismic data to be processed, refers to the seismic data for carrying out surface wave dispersion processing.
In the application one embodiment, the seismic data to be processed is frequency domain.In another embodiment of the application
In, the seismic data to be processed is time-domain, need to be converted to the seismic data to be processed of time-domain frequency domain to
Handle seismic data.
Selecting for the seismic target earthquakes road, mainly according to the seismic data to be processed, is selected in the seismic data
Each reference frequency point where seismic channel.
S2: the adjacent seismic channel quantity threshold and phase difference threshold in the seismic target earthquakes road, the adjacent seismic channel are determined
Quantity threshold is equal to 2N, N >=1.
The adjacent seismic channel quantity threshold and phase difference threshold, experience can be combined according to the actual situation certainly by implementing personnel
Row determines.Wherein, the value range of the phase difference threshold is to be less than or equal to 360 degree more than or equal to 0 degree.
In the application one embodiment, the adjacent seismic channel quantity threshold and phase difference in the determination seismic target earthquakes road
The mode of threshold value may include:
Surface wave dispersion processing is carried out to the seismic data to be processed, obtains frequency-phase velocity frequency dispersion spectrum.
According to the frequency-phase velocity frequency dispersion spectrum, determine the seismic target earthquakes road adjacent seismic channel quantity threshold and
Phase difference threshold.
Typically, for the weaker situation of the frequency-phase velocity frequency dispersion spectrum noise, the adjacent seismic channel number threshold
Value and phase difference threshold can choose smaller value.Conversely, for the stronger situation of noise of the frequency-phase velocity frequency dispersion spectrum,
The adjacent seismic channel quantity threshold and phase difference threshold can choose the larger value.
Certainly, in the application other embodiments, the determination of specific adjacent seismic channel quantity threshold and phase difference threshold
Mode and determining standard, can voluntarily be determined with operation experience according to the actual situation by implementation personnel.
S3: according to the adjacent seismic channel quantity threshold, choosing M adjacent seismic channels in the seismic target earthquakes road altogether, and N≤
M≤2N。
The road number of the adjacent seismic channel, the adjacent seismic channel including seismic target earthquakes road, the adjacent seismic channel is little
In the adjacent seismic channel number threshold value N, adjacent seismic channel number, which is added, is equal to the M.
If the seismic channel number in the seismic target earthquakes road is not less than N, respectively there is N number of adjacent seismic channel on both sides, at this point, M
=2N.
If the seismic channel number of seismic target earthquakes road side is 0, and the seismic channel of the seismic target earthquakes road other side
Number is not less than N, then M=N.
The value of the M is determined according to the seismic channel number in actual seismic target earthquakes road, as long as the seismic channel number of any side is not
Less than N, then the adjacent seismic channel number of this side is with regard to value N, if the seismic channel number of side is less than N, this side is adjacent to
Shake road number is equal to the seismic channel number of this side.
S4: the M adjacent seismic channels and the seismic target earthquakes road are calculated separately in the phase difference of each reference frequency point.
Each Frequency point refers to the Frequency point for including in the corresponding seismic data in the seismic target earthquakes road.
S5: the phase difference is compared with the phase difference threshold, determines the phase on the reference frequency point
Difference is greater than the road number of the seismic channel of the phase difference threshold.
The reference frequency point includes the Frequency point for including in the corresponding seismic data in seismic target earthquakes road.
S6: whether the road number is greater than preset phase equalization threshold value, determines corresponding ginseng according to comparison result
Than Frequency point whether be surface wave frequency dispersion curve cumulative frequency point.
The phase equalization threshold value, which is referred to, is greater than the phase difference threshold with the phase difference in the seismic target earthquakes road
The threshold value of seismic channel road number.
In the application one embodiment, whether the road number is greater than preset phase equalization threshold value, according to comparing
As a result determine corresponding reference frequency point whether be surface wave frequency dispersion curve cumulative frequency point mode, may include:
If the road number is greater than preset phase equalization threshold value, corresponding reference frequency point is not surface wave frequency dispersion curve
Cumulative frequency point.
If the road number is less than or equal to preset phase equalization threshold value, corresponding reference frequency point is surface wave dispersion
The cumulative frequency point of curve.
It, can be with using a kind of embodiment of the data processing method of the surface wave frequency dispersion curve provided in the various embodiments described above
By determining that the corresponding seismic channel of seismic data in the phase equalization of each Frequency point, determines adjacent seismic channel in seismic data
Phase inconsistent (phase difference is larger) Frequency point.In seismic data cumulative process, the phase for the adjacent seismic channel that only adds up
The corresponding seismic data of Frequency point of position consistent (phase difference is small) can effectively eliminate making an uproar in frequency-phase velocity frequency dispersion spectrum
Part point, to effectively improve the precision of the surface wave frequency dispersion curve picked up from the frequency-phase velocity frequency dispersion spectrum.
Fig. 2 is a kind of method stream of the data processing method of the surface wave frequency dispersion curve provided in another embodiment of the application
Journey schematic diagram, specifically, may include:
S1: the seismic target earthquakes road in selected frequency domain seismic data to be processed.
S2: the adjacent seismic channel quantity threshold and phase difference threshold in the seismic target earthquakes road, the adjacent seismic channel are determined
Quantity threshold is equal to 2N, N >=1.
S3: according to the adjacent seismic channel quantity threshold, choosing M adjacent seismic channels in the seismic target earthquakes road altogether, and N≤
M≤2N。
S4: the M adjacent seismic channels and the seismic target earthquakes road are calculated separately in the phase difference of each reference frequency point.
S5: the phase difference is compared with the phase difference threshold, determines that phase difference is greater than institute on the reference frequency point
State the road number of the seismic channel of phase difference threshold.
S6: whether the road number is greater than preset phase equalization threshold value, determines corresponding ginseng according to comparison result
Than Frequency point whether be surface wave frequency dispersion curve cumulative frequency point.
S7: the corresponding seismic data of the cumulative frequency point is added up, the seismic data after being added up, to described
Seismic data after cumulative carries out frequency dispersion processing, obtains surface wave frequency-phase velocity spectrum, from the surface wave frequency-phase velocity spectrum
Pick up surface wave frequency dispersion curve.
Using embodiment provided by the above embodiment, the available surface wave frequency for eliminating noise section-phase velocity spectrum,
The higher surface wave frequency dispersion curve of precision that can be picked up from the surface wave frequency-phase velocity spectrum.
A kind of embodiment of the data processing method of surface wave frequency dispersion curve of the various embodiments described above offer, Ke Yitong are provided
It crosses and determines that the corresponding seismic channel of seismic data in the phase equalization of each Frequency point, determines adjacent seismic channel in seismic data
The Frequency point of phase inconsistent (phase difference is larger).In seismic data cumulative process, the phase for the adjacent seismic channel that only adds up
The unanimously corresponding seismic data of the Frequency point of (phase difference is small) can effectively eliminate the noise in frequency-phase velocity frequency dispersion spectrum
Part, to effectively improve the precision of the surface wave frequency dispersion curve picked up from the frequency-phase velocity frequency dispersion spectrum.
In one example of the application, using herein described method, frequency dispersion point is carried out to seismic data as shown in Figure 5
Analysis.
In Fig. 5, what abscissa indicated is seismic channel number, and what ordinate indicated is the time.
Fig. 6 is that surface wave frequency-phase velocity obtained without use the method in this example as a comparison is composed.
Fig. 7 is to use herein described method in this example, carries out dispersion analysis to the seismic data in Fig. 5, obtains
A kind of surface wave frequency-phase velocity spectrum.
During the corresponding dispersion analysis of Fig. 7, for 3, phase equalization threshold value is 1 for the adjacent seismic channel number threshold value that uses.
Fig. 8 is to use herein described method in this example, carries out dispersion analysis to the seismic data in Fig. 5, obtains
A kind of surface wave frequency-phase velocity spectrum.
During the corresponding dispersion analysis of Fig. 8, for 6, phase equalization threshold value is 1 for the adjacent seismic channel number threshold value that uses.
In above-mentioned Fig. 6, Fig. 7, Fig. 8, the number on abscissa represents the numerical value of phase velocity, the number equal generation on ordinate
The numerical value of table frequency.
By the way that tri- width figure of Fig. 6, Fig. 7, Fig. 8 is compared it can be found that surface wave frequency-phase velocity that Fig. 7 and Fig. 8 are obtained
Distribution of color ratio Fig. 6 in spectrogram more evenly illustrates that noise has obtained preferable compacting.In addition, occurring the two of white in Fig. 8
Point illustrates that wave energy enhances, and signal-to-noise ratio improves, more preferable to pressing result ratio Fig. 7 of noise.
Pass through data processed result obtained in examples detailed above, it may be said that it is bright, frequency dispersion is carried out using herein described method
Analysis can effectively improve surface wave frequency-phase velocity spectrum signal-to-noise ratio, so as to pick up the higher face of accuracy from spectrum
Wave dispersion curve.
Based on the data processing method of surface wave frequency dispersion curve described herein, the application provides a kind of surface wave frequency dispersion curve
Data processing equipment, the module can be used in the functional unit of near surface structure inverting.Fig. 3 is the application one implementation
The modular structure schematic diagram of the data processing equipment of a kind of surface wave frequency dispersion curve provided in example, specifically, as shown in figure 3, institute
Stating device may include:
Seismic channel chosen module 101 can be used for the seismic target earthquakes road in selected frequency domain seismic data to be processed.
Threshold determination module 102, the threshold value being determined in seismic data to be processed.
Phase difference calculating module 103, can be used for calculating separately the M adjacent seismic channels and the seismic target earthquakes road exists
The phase difference of each reference frequency point.
Road number determining module 104 can be used for for the phase difference being compared with the phase difference threshold, determine
Phase difference is greater than the road number of the seismic channel of the phase difference threshold on the reference frequency point.
Frequency point determining module 105, can be used for whether road number described in comparison is greater than preset phase equalization threshold value, root
According to comparison result determine corresponding reference frequency point whether be surface wave frequency dispersion curve cumulative frequency point.
It, can be automatic using a kind of embodiment of the data processing equipment of surface wave frequency dispersion curve provided by the above embodiment
The method is executed, the excessive participation of implementation personnel is not needed, effectively increases user experience.
In another embodiment of the application, the threshold determination module may include:
First threshold determining module is determined for the adjacent seismic channel quantity threshold and phase in the seismic target earthquakes road
Poor threshold value.
Second threshold determining module is determined for phase equalization threshold value.
In another embodiment of the application, the first threshold determining module may include:
Frequency dispersion processing module can be used for carrying out surface wave dispersion processing to the seismic data to be processed, obtain frequency-
Phase velocity frequency dispersion spectrum；
Threshold value chooses module, can be used for determining the seismic target earthquakes road according to the frequency-phase velocity frequency dispersion spectrum
Adjacent seismic channel quantity threshold and phase difference threshold.
Fig. 4 is a kind of module knot of the data processing equipment of the surface wave frequency dispersion curve provided in another embodiment of the application
Structure schematic diagram, specifically, the apparatus may include:
Seismic channel chosen module 101 can be used for the seismic target earthquakes road in selected frequency domain seismic data to be processed.
Threshold determination module 102, the threshold value being determined in seismic data to be processed.
Phase difference calculating module 103, can be used for calculating separately the M adjacent seismic channels and the seismic target earthquakes road exists
The phase difference of each reference frequency point.
Road number determining module 104 can be used for for the phase difference being compared with the phase difference threshold, determine
Phase difference is greater than the road number of the seismic channel of the phase difference threshold on the reference frequency point.
Frequency point determining module 105, can be used for whether road number described in comparison is greater than preset phase equalization threshold value, root
According to comparison result determine corresponding reference frequency point whether be surface wave frequency dispersion curve cumulative frequency point.
Surface wave frequency dispersion curve generation module 106, for generating surface wave frequency dispersion curve.
In another embodiment of the application, the surface wave frequency dispersion curve generation module may include:
Seismic data accumulator module, for the corresponding seismic data of cumulative frequency point to add up, after being added up
Seismic data.
Frequency dispersion processing module obtains surface wave frequency-phase velocity for carrying out frequency dispersion processing to the seismic data after described add up
Degree spectrum.
Dispersion curve pickup model, for picking up surface wave frequency dispersion curve from the surface wave frequency-phase velocity spectrum.
Using a kind of data processing equipment for surface wave frequency dispersion curve that the various embodiments described above provide, can be performed automatically described
Method does not need the excessive participation of implementation personnel, so that it may which the surface wave frequency dispersion curve for obtaining degree of precision effectively increases user
Experience.
In the data processing equipment of the surface wave frequency dispersion curve, the selected target seismic channel determines the seismic target earthquakes
The adjacent seismic channel quantity threshold and phase difference threshold in road calculate phase difference, determine the phase difference on the reference frequency point
It is counted to greater than the road of the seismic channel of the phase difference threshold and determines whether corresponding reference frequency point is the tired of surface wave frequency dispersion curve
The extension of the embodiment of Frequency point is added to be referred to the associated description of preceding method.
Although mentioning the data processing method of different surface wave frequency dispersion curves in teachings herein, from selected target earthquake
Road, the adjacent seismic channel quantity threshold for determining the seismic target earthquakes road and phase difference threshold calculate phase difference, determine described
The road that phase difference is greater than the seismic channel of the phase difference threshold on reference frequency point counts to whether determine corresponding reference frequency point
For the description of the various sequential manners of the cumulative frequency point of surface wave frequency dispersion curve, data acquisition/processing/way of output etc., still,
The application is not limited to be situation described in professional standard or embodiment etc., and certain professional standards or use are made by oneself
Right way of conduct formula or embodiment description practice processes on embodiment modified slightly also may be implemented above-described embodiment it is identical, etc.
The implementation result being anticipated that after same or close or deformation.Using these modifications or deformed data acquisition, processing, exports, sentences
The embodiment of disconnected mode etc., still 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 may include more or less operating procedure.The step of enumerating in embodiment sequence is only numerous steps
One of execution sequence mode, does not represent and unique executes sequence.It, can when device or client production in practice executes
To execute or parallel execute (such as at parallel processor or multithreading according to embodiment or method shown in the drawings sequence
The environment of reason, even distributed data processing environment).The terms "include", "comprise" or its any other variant are intended to contain
Lid non-exclusive inclusion, so that process, method, product or equipment including a series of elements are not only wanted including those
Element, but also including other elements that are not explicitly listed, or further include for this process, method, product or equipment
Intrinsic element.In the absence of more restrictions, be not precluded include the process, method of the element, product or
There is also other identical or equivalent elements in person's equipment.
Device that above-described embodiment illustrates or module etc. can specifically realize 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 when description apparatus above with function to retouch respectively
It states.Certainly, the function of each module can be realized in the same or multiple software and or hardware when implementing the application,
The module for realizing same function can be realized by the combination of multiple submodule etc..Installation practice described above is only
Schematically, for example, the division of the module, only a kind of logical function partition, can there is other draw in actual implementation
The mode of dividing, such as multiple module or components can be combined or can be integrated into another system, or some features can be ignored,
Or it does not execute.
It is also known in the art that other than 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 so that controller with logic gate, switch, specific integrated circuit, programmable
Logic controller realizes identical function with the form for being embedded in microcontroller etc..Therefore this controller is considered one kind
Hardware component, and the structure that the device for realizing various functions that its inside includes can also be considered as in hardware component.Or
Person even, can will be considered as realizing the device of various functions either the software module of implementation method can be hardware again
Structure in component.
The application can describe in the general context of computer-executable instructions executed by a computer, such as program
Module.Generally, program module includes routines performing specific tasks or implementing specific abstract data types, programs, objects, group
Part, data structure, class etc..The application can also be practiced in a distributed computing environment, in these distributed computing environments,
By executing task by the connected remote processing devices of communication network.In a distributed computing environment, program module can
To be located in the local and remote computer storage media including storage equipment.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
It realizes by means of software and necessary general hardware platform.Based on this understanding, the technical solution essence of the application
On in other words the part that contributes to existing technology can be embodied in the form of software products, the computer software product
It can store in storage medium, such as ROM/RAM, magnetic disk, CD, including some instructions are used so that a computer equipment
(can be personal computer, mobile terminal, server or the network equipment etc.) executes each embodiment of the application or implementation
Method described in certain parts of example.
Each embodiment in this specification is described in a progressive manner, the same or similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.The application can be used for crowd
In mostly general or special purpose computing system environments or configuration.Such as: personal computer, server computer, handheld device or
Portable device, laptop device, multicomputer system, microprocessor-based system, set top box, programmable electronics set
Standby, network PC, minicomputer, mainframe computer, distributed computing environment including any of the above system or equipment etc..
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application there are many deformation and
Variation is without departing from spirit herein, it is desirable to which the attached claims include these deformations and change without departing from the application's
Spirit.
Claims (9)
1. a kind of data processing method of surface wave frequency dispersion curve, which is characterized in that the described method includes:
Seismic target earthquakes road in selected frequency domain seismic data to be processed；
Determine the adjacent seismic channel quantity threshold and phase difference threshold in the seismic target earthquakes road, the adjacent seismic channel quantity threshold
Equal to 2N, N >=1；
According to the adjacent seismic channel quantity threshold, M adjacent seismic channels, N≤M≤2N are chosen altogether in the seismic target earthquakes road；
The M adjacent seismic channels and the seismic target earthquakes road are calculated separately in the phase difference of each reference frequency point；
The phase difference is compared with the phase difference threshold, determines that phase difference is greater than institute on the reference frequency point
State the road number of the seismic channel of phase difference threshold；
Compare whether the road number is greater than preset phase equalization threshold value, corresponding reference frequency point is determined according to comparison result
Whether be surface wave frequency dispersion curve cumulative frequency point.
2. a kind of data processing method of surface wave frequency dispersion curve as described in claim 1, which is characterized in that the method is also wrapped
It includes:
The corresponding seismic data of the cumulative frequency point is added up, the seismic data after being added up, after described add up
Seismic data carry out frequency dispersion processing, obtain surface wave frequency-phase velocity spectrum, pick up face from the surface wave frequency-phase velocity spectrum
Wave dispersion curve.
3. a kind of data processing method of surface wave frequency dispersion curve as described in claim 1, which is characterized in that described in the comparison
Whether road number is greater than preset phase equalization threshold value, determines whether corresponding reference frequency point is surface wave frequency according to comparison result
The mode of the cumulative frequency point of non-dramatic song line, comprising:
If the road number is greater than preset phase equalization threshold value, corresponding reference frequency point is not the tired of surface wave frequency dispersion curve
Add Frequency point；
If the road number is less than or equal to preset phase equalization threshold value, corresponding reference frequency point is surface wave frequency dispersion curve
Cumulative frequency point.
4. a kind of data processing method of surface wave frequency dispersion curve as described in any one of claim 1 or 3, feature exist
In the phase equalization threshold value includes:
It is greater than the threshold value of the seismic channel road number of the phase difference threshold with the phase difference in the seismic target earthquakes road.
5. a kind of data processing method of surface wave frequency dispersion curve as described in claim 1, which is characterized in that described in the determination
The adjacent seismic channel quantity threshold in seismic target earthquakes road and the mode of phase difference threshold include:
Surface wave dispersion processing is carried out to the seismic data to be processed, obtains frequency-phase velocity frequency dispersion spectrum；
According to the frequency-phase velocity frequency dispersion spectrum, the adjacent seismic channel quantity threshold and phase in the seismic target earthquakes road are determined
Poor threshold value.
6. a kind of data processing equipment of surface wave frequency dispersion curve, which is characterized in that described device includes:
Seismic channel chosen module, for selecting the seismic target earthquakes road in frequency domain seismic data to be processed；
Threshold determination module, for determining the threshold value in seismic data to be processed；
Wherein, the threshold determination module includes:
First threshold determining module, for determining the adjacent seismic channel quantity threshold and phase difference threshold in the seismic target earthquakes road；
Second threshold determining module, for determining phase equalization threshold value；
Phase difference calculating module, for calculating separately M adjacent seismic channels and the seismic target earthquakes road in each reference frequency point
Phase difference；
Road number determining module is determined for the phase difference to be compared with the phase difference threshold in the reference frequency
Phase difference is greater than the road number of the seismic channel of the phase difference threshold on rate point；
Whether Frequency point determining module is greater than preset phase equalization threshold value for the road number, according to comparison result
Determine corresponding reference frequency point whether be surface wave frequency dispersion curve cumulative frequency point.
7. a kind of data processing equipment of surface wave frequency dispersion curve as claimed in claim 6, which is characterized in that the first threshold
Determining module includes:
Frequency dispersion processing module obtains frequency-phase velocity frequency for carrying out surface wave dispersion processing to the seismic data to be processed
Dissipate spectrum；
Threshold value chooses module, for composing according to the frequency-phase velocity frequency dispersion, determines the adjacent earthquake in the seismic target earthquakes road
Road quantity threshold and phase difference threshold.
8. a kind of data processing equipment of surface wave frequency dispersion curve as claimed in claim 6, which is characterized in that described device is also wrapped
It includes:
Surface wave frequency dispersion curve generation module, for generating surface wave frequency dispersion curve.
9. a kind of data processing equipment of surface wave frequency dispersion curve as claimed in claim 8, which is characterized in that the surface wave dispersion
Curve generation module includes:
Seismic data accumulator module, for the corresponding seismic data of cumulative frequency point to be added up, the earthquake after being added up
Data；
Frequency dispersion processing module obtains surface wave frequency-phase velocity for carrying out frequency dispersion processing to the seismic data after described add up
Spectrum；
Dispersion curve pickup model, for picking up surface wave frequency dispersion curve from the surface wave frequency-phase velocity spectrum.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611244008.8A CN106646615B (en) | 2016-12-29 | 2016-12-29 | A kind of data processing method and device of surface wave frequency dispersion curve |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611244008.8A CN106646615B (en) | 2016-12-29 | 2016-12-29 | A kind of data processing method and device of surface wave frequency dispersion curve |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106646615A CN106646615A (en) | 2017-05-10 |
| CN106646615B true CN106646615B (en) | 2018-12-25 |
Family
ID=58835586
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611244008.8A Active CN106646615B (en) | 2016-12-29 | 2016-12-29 | A kind of data processing method and device of surface wave frequency dispersion curve |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106646615B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6945895B2 (en) * | 2017-10-12 | 2021-10-06 | サウザン・ユニバーシティ・オブ・サイエンス・アンド・テクノロジー | Surface wave exploration methods and terminal devices |
| CN108562938B (en) | 2018-03-23 | 2019-09-06 | 中国石油天然气股份有限公司 | A kind of method, apparatus and system for eliminating frequency dispersion effect |
| CN109799530A (en) * | 2018-12-25 | 2019-05-24 | 核工业北京地质研究院 | Rayleigh waves dispersion curve inversion method for seismic surface wave exploration |
| CN112748459A (en) * | 2019-10-29 | 2021-05-04 | 中国石油天然气集团有限公司 | Method and device for automatically picking up frequency dispersion curve |
| CN112987090A (en) * | 2019-12-02 | 2021-06-18 | 中国石油天然气集团有限公司 | Method and device for picking up surface wave frequency dispersion curve |
| CN112987089A (en) * | 2019-12-02 | 2021-06-18 | 中国石油天然气集团有限公司 | Method and device for enhancing surface wave frequency dispersion spectrum signal |
| CN111290017B (en) * | 2020-03-04 | 2020-12-15 | 南方科技大学 | Surface wave exploration method for jointly extracting Rayleigh wave frequency dispersion characteristics through seismic electric wave field |
| US11754744B2 (en) | 2020-03-04 | 2023-09-12 | Southern University Of Science And Technology | Surface wave prospecting method for jointly extracting Rayleigh wave frequency dispersion characteristics by seismoelectric field |
| CN111766625B (en) * | 2020-07-06 | 2021-12-14 | 中国科学技术大学 | Seismic background noise dispersion curve extraction method based on deep learning |
| CN112051610B (en) * | 2020-10-21 | 2021-06-11 | 中国地质大学（北京） | Vector field multi-mode surface wave frequency dispersion calculation method and system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102478670A (en) * | 2010-11-24 | 2012-05-30 | 中国石油天然气集团公司 | Method for forecasting reservoir fluid property through earthquake attenuation attribute |
| CN103048684A (en) * | 2011-10-11 | 2013-04-17 | 中国石油化工股份有限公司 | Multi-component seismic data surface wave pressing method |
| CN103238158A (en) * | 2010-12-01 | 2013-08-07 | 埃克森美孚上游研究公司 | Simultaneous source inversion for marine streamer data with cross-orrelation objective function |
| CN104614769A (en) * | 2015-02-05 | 2015-05-13 | 中铁隧道集团有限公司 | Beam-forming filtering method for suppressing seismic surface waves |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10338242B2 (en) * | 2015-05-20 | 2019-07-02 | Conocophillips Company | Surface wave tomography using sparse data acquisition |
-
2016
- 2016-12-29 CN CN201611244008.8A patent/CN106646615B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102478670A (en) * | 2010-11-24 | 2012-05-30 | 中国石油天然气集团公司 | Method for forecasting reservoir fluid property through earthquake attenuation attribute |
| CN103238158A (en) * | 2010-12-01 | 2013-08-07 | 埃克森美孚上游研究公司 | Simultaneous source inversion for marine streamer data with cross-orrelation objective function |
| CN103048684A (en) * | 2011-10-11 | 2013-04-17 | 中国石油化工股份有限公司 | Multi-component seismic data surface wave pressing method |
| CN104614769A (en) * | 2015-02-05 | 2015-05-13 | 中铁隧道集团有限公司 | Beam-forming filtering method for suppressing seismic surface waves |
Non-Patent Citations (1)
| Title |
|---|
| "逆频散方法压制频散面波";周腾飞 等;《石油地球物理勘探》;20121231;第47卷(第6期);第882-886页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106646615A (en) | 2017-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106646615B (en) | A kind of data processing method and device of surface wave frequency dispersion curve | |
| Callegari et al. | Spectral properties of chaos-based FM signals: Theory and simulation results | |
| Kardoutsou et al. | A new pulse indicator for the classification of ground motions | |
| CN107728210B (en) | Multiple instruments acquire the determination method and apparatus that road is lacked in data | |
| Rezaeian et al. | Validation of simulated earthquake ground motions based on evolution of intensity and frequency content | |
| CN109061724A (en) | A kind of seismic data noise-reduction method based on adaptive variation mode decomposition | |
| CN104599681B (en) | The method and apparatus of audio frequency process | |
| CN110457524A (en) | Model generating method, video classification methods and device | |
| CN105244041A (en) | Song audition evaluation method and device | |
| CN104484232B (en) | A kind of method and device for improving image reconstruction speed | |
| CN107229075B (en) | The determination method and apparatus of Depth Domain seismic wavelet | |
| CA2462325A1 (en) | Reverberation sound generating apparatus | |
| CN104364845A (en) | Processing apparatus, processing method, program, computer readable information recording medium and processing system | |
| CN105911590B (en) | A kind of compression process method and device of harmonic noise | |
| CN104407278A (en) | Extraction method and device of ultrasonic wave characteristic parameters | |
| CN106338768B (en) | It is a kind of to generate the processing method of reservoir prediction attribute data, apparatus and system | |
| CN109901222A (en) | Frequency method, storage medium and equipment are opened up in earthquake based on dynamic dictionary match tracing | |
| CN108469633A (en) | The computational methods and device of a kind of ground interval quality factors | |
| CN104111477B (en) | Method and system for acquiring sliding time for harmonic suppression | |
| CN107688197B (en) | Thin bed prediction method and apparatus | |
| CN106197466B (en) | Data processing method and device and terminal | |
| CN108062959A (en) | A kind of noise reduction method and device | |
| Shin et al. | Transformation approach to model online gaming traffic | |
| CN108226636A (en) | Automatic filter method and apparatus | |
| Alimuradov | Research of frequency-selective properties of empirical mode decomposition methods for speech signals' pitch frequency estimation |
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 |