CN106646615A - Data processing method and data processing device for surface wave dispersion curve - Google Patents
Data processing method and data processing device for surface wave dispersion curve Download PDFInfo
- Publication number
- CN106646615A CN106646615A CN201611244008.8A CN201611244008A CN106646615A CN 106646615 A CN106646615 A CN 106646615A CN 201611244008 A CN201611244008 A CN 201611244008A CN 106646615 A CN106646615 A CN 106646615A
- Authority
- CN
- China
- Prior art keywords
- threshold values
- frequency
- surface wave
- phase
- phase difference
- 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
- 239000006185 dispersion Substances 0.000 title claims abstract description 127
- 238000012545 processing Methods 0.000 title claims abstract description 25
- 238000003672 processing method Methods 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 claims abstract description 55
- 238000001228 spectrum Methods 0.000 claims abstract description 41
- 230000001186 cumulative effect Effects 0.000 claims description 35
- 230000000052 comparative effect Effects 0.000 claims description 7
- 238000010130 dispersion processing Methods 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 5
- 238000011017 operating method Methods 0.000 description 4
- 239000000047 product Substances 0.000 description 4
- 238000003860 storage Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000009826 distribution 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
- 238000003825 pressing Methods 0.000 description 1
- 230000000750 progressive effect Effects 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
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a data processing method and a data processing device for a surface wave dispersion curve. The method comprises the steps of selecting target seismic traces in to-be-processed seismic data in a frequency domain; determining an adjacent seismic trace number threshold and a phase difference threshold of the target seismic traces; totally selecting M adjacent seismic traces in the target seismic traces; respectively calculating phase differences between the M adjacent seismic traces and the target seismic traces at each reference frequency point; determining the number of the seismic traces in which the phase difference on the reference frequency point is larger than the phase difference threshold; and determining whether the corresponding reference frequency point is an added frequency point of the surface wave dispersion curve. By means of the data processing method and the data processing device, the noise part in a frequency-phase-speed frequency dispersion spectrum can be effectively eliminated, thereby effectively improving precision of the surface wave dispersion curve which is obtained from the frequency-phase-speed frequency dispersion spectrum.
Description
Technical field
The application is related to seismic data processing technology field, more particularly to a kind of data processing method of surface wave frequency dispersion curve
And device.
Background technology
In seismic data processing technology, it is to carry out earth formation analysis to commonly use to carry out near surface structure inverting using face ripple
Technology.Carrying out near surface inverting using face ripple 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 that to carry out near surface structure using the surface wave frequency dispersion curve for generating anti-
Drill.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 for improving inversion result precision.
In prior art, in order to obtain surface wave frequency dispersion curve, usually geological data is added up, line frequency of going forward side by side dissipates place
Reason, obtains frequency-phase velocity frequency dispersion spectrum, and from the frequency-phase velocity frequency dispersion spectrum surface wave frequency dispersion curve is picked up.But, tired
Plus during, due in the geological data in addition to linear face ripple, also including other nonlinear seismic waves, these ground
The phase place of seismic wave is larger with the phase difference of face ripple, so will result in the interference of opposite ripple accumulation result, causes frequency-phase velocity
Noise section is included in degree frequency dispersion spectrum, and then reduces the precision of the surface wave frequency dispersion curve that pickup is obtained.
At least there are the following problems in prior art：In prior art, when carrying out cumulative to geological data, in geological data
Nonlinear seismic wave is contained, the nonlinear seismic wave is larger with face wave phase gap, opposite ripple accumulation result is caused
Interference, causes in frequency-phase velocity frequency dispersion spectrum comprising noise section, and then causes from the frequency-phase velocity frequency dispersion spectrum to pick up
The surface wave frequency dispersion curve precision for obtaining is relatively low.
The content of the invention
The purpose of the embodiment of the present application is to provide a kind of data processing method and device of surface wave frequency dispersion curve, to eliminate frequency
Noise section in rate-phase velocity frequency dispersion spectrum, so as to pick up the surface wave frequency dispersion curve for obtaining higher precision.
The embodiment of the present application provides what a kind of data processing method and device of surface wave frequency dispersion curve was realized in：
A kind of data processing method of surface wave frequency dispersion curve, methods described includes：
In the pending corresponding seismic channel of geological data, selected target seismic channel；
Determine the adjacent seismic channel number threshold values and phase difference threshold values in the seismic target earthquakes road, the adjacent seismic channel number
Threshold values is equal to N, N >=2；
According to the phase difference threshold values, determine seismic target earthquakes road with the M adjacent seismic channel in the seismic target earthquakes road every
The phase equalization of individual Frequency point, N≤M≤2N；
According to the determination result of the phase equalization, determine in each Frequency point, with seismic target earthquakes road phase place
The road number of inconsistent adjacent seismic channel；
According to described in each Frequency point, the road number of the adjacent seismic channel inconsistent with seismic target earthquakes road phase place, really
Make and may participate in cumulative Frequency point.
In preferred embodiment, methods described also includes：
The corresponding geological data of cumulative Frequency point will be may participate in carries out the scattered process of cumulative line frequency of going forward side by side, and obtains surface wave dispersion
Curve.
It is described according to the phase difference threshold values in preferred embodiment, determine the M in seismic target earthquakes road and the seismic target earthquakes road
Individual adjacent seismic channel the phase equalization of each Frequency point mode, including：
At described each Frequency point, according to the phase difference threshold values, each in the M neighboring track is determined respectively
The phase equalization in seismic channel and the seismic target earthquakes road.
It is described according to the phase difference threshold values in preferred embodiment, each ground in the M neighboring track is determined respectively
The mode in road and the phase equalization in the seismic target earthquakes road is shaken, including：
The phase difference value of described each seismic channel and the seismic target earthquakes road is calculated respectively；
According to the phase difference value and phase difference threshold values, 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 values in preferred embodiment, in determining the M neighboring track
The mode of each seismic channel and the phase equalization in the seismic target earthquakes road, including：
If the phase difference value is less than the phase difference threshold values, it is determined that the result of phase equalization is that phase place is consistent；
If the phase difference value is not less than the phase difference threshold values, it is determined that the result of phase equalization differs for phase place
Cause.
In preferred embodiment, in each Frequency point described in the basis, the phase inconsistent with seismic target earthquakes road phase place
The road number of adjacent seismic channel, determines the mode of the Frequency point that may participate in cumulative, including：
Determine the inconsistent road number threshold values of phase place；
According to the road number of the adjacent seismic channel inconsistent with seismic target earthquakes road phase place in each Frequency point, and institute
State the inconsistent road number threshold values of phase place, it is determined that described each Frequency point whether may participate in it is cumulative, according to determine result, counting to join
With cumulative Frequency point.
In preferred embodiment, the basis is in the adjacent earthquake inconsistent with seismic target earthquakes road phase place of each Frequency point
The road number in road, and the inconsistent road number threshold values of the phase place, it is determined that whether described each Frequency point may participate in cumulative mode, bag
Include：
If the road number of the adjacent seismic channel inconsistent with seismic target earthquakes road phase place, no more than described phase place differs
Cause road number threshold values, it is determined that corresponding Frequency point may participate in cumulative；
If the road number of the adjacent seismic channel inconsistent with seismic target earthquakes road phase place, inconsistent more than the phase place
Road number threshold values, it is determined that corresponding Frequency point can not be accumulated.
In preferred embodiment, the adjacent seismic channel number threshold values for determining the seismic target earthquakes road and phase difference threshold values
Mode includes：
Surface wave dispersion process is carried out to the pending geological data, frequency-phase velocity frequency dispersion spectrum is obtained；
According to the frequency-phase velocity frequency dispersion spectrum, determine the seismic target earthquakes road adjacent seismic channel number threshold values and
Phase difference threshold values.
A kind of data processing equipment of surface wave frequency dispersion curve, described device includes：
Seismic channel chosen module, in the pending corresponding seismic channel of geological data, selected target seismic channel；
Threshold determination module, for determining the adjacent seismic channel number threshold values and phase difference threshold values in the seismic target earthquakes road,
It is additionally operable to determine the inconsistent road number threshold values of phase place；
Phase equalization determining module, for according to the phase difference threshold values, determining seismic target earthquakes road with target ground
Phase equalization of the M adjacent seismic channel in shake road in each Frequency point；
Road number determining module, for according to the determination result of the phase equalization, determining in each Frequency point, with institute
State the road number of the inconsistent adjacent seismic channel of seismic target earthquakes road phase place；
Frequency point determining module, for according to described in each Frequency point, it is inconsistent with seismic target earthquakes road phase place
The road number of adjacent seismic channel, determines and may participate in cumulative Frequency point.
In preferred embodiment, described device also includes：
Surface wave frequency dispersion curve generation module, carries out adding up simultaneously for will may participate in the corresponding geological data of cumulative Frequency point
Frequency dispersion process 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 the phase difference of described each seismic channel and the seismic target earthquakes road respectively
Value；
Phase equalization discrimination module, for according to the phase difference value and phase difference threshold values, determining that the M is individual adjacent
The phase equalization in each seismic channel in road and the seismic target earthquakes road；
Recognition result output module, for output phase uniformity result is determined.
In preferred embodiment, the Frequency point determining module includes：
The inconsistent road number threshold determination module of phase place, for determining the inconsistent road number threshold values of phase place；
Can additive discrimination module, for according to described inconsistent with seismic target earthquakes road phase place in each Frequency point
The road number of adjacent seismic channel, and the inconsistent road number threshold values of the phase place, it is determined that described each Frequency point whether may participate in it is cumulative；
Can cumulative frequency point statistical module, for whether may participate in cumulative determination result according to described each Frequency point,
Count and may participate in cumulative Frequency point.
A kind of data processing method of the surface wave frequency dispersion curve provided using the embodiment of the present application, can be by determining earthquake
The corresponding seismic channel of data determines the phase place of adjacent seismic channel in geological data not in the phase equalization threshold values of each Frequency point
The consistent Frequency point of (phase difference is larger).In geological data cumulative process, the consistent (phase of phase place of the adjacent seismic channel that only adds up
Potential difference is different little) the corresponding geological data of Frequency point, can effectively eliminate the noise section in frequency-phase velocity frequency dispersion spectrum, from
And effectively improve from the frequency-phase velocity frequency dispersion spectrum and pick up the precision of the surface wave frequency dispersion curve for obtaining.Using the application reality
A kind of data processing equipment of surface wave frequency dispersion curve of example offer is provided, methods described can be automatically performed, it is not necessary to implement personnel
Excessive participation, it is possible to obtain the surface wave frequency dispersion curve of degree of precision, effectively increase Consumer's Experience.
Description of the drawings
In order to be illustrated more clearly that the embodiment of the present application or technical scheme of the prior art, below will be to embodiment or existing
The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments described in application, for those of ordinary skill in the art, in the premise for not paying creative labor
Under, can be with according to these other accompanying drawings of accompanying drawings acquisition.
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 the application another embodiment
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 the application another embodiment
Structure schematic diagram.
Fig. 5 is the pending geological data provided in one example of the application；
Fig. 6, Fig. 7, Fig. 8 be in one example of the application to Fig. 5 in pending geological data carry out surface wave dispersion analysis
The face wave frequency rate, for obtaining-phase velocity spectrum.
Specific embodiment
The embodiment of the present application provides a kind of data processing method and device of surface wave frequency dispersion curve.
In order that those skilled in the art more fully understand the technical scheme in the application, below in conjunction with the application reality
The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described enforcement
Example is only some embodiments of the present application, rather than the embodiment of whole.Based on the embodiment in the application, this area is common
The every other embodiment that technical staff is obtained under the premise of creative work is not made, should all belong to the application protection
Scope.
Fig. 1 is that a kind of a kind of method flow of embodiment of data processing method of herein described surface wave frequency dispersion curve is illustrated
Figure.Although this application provides such as following embodiments or method operating procedure shown in the drawings or apparatus structure, based on routine
Either more or less operating procedures or module list can be included in methods described or device without the need for performing creative labour
Unit.In the step of there is no necessary causality 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 sequence shown in the drawings or modular structure.Described method or modular structure
When device in practice or end product application, can carry out according to embodiment or method shown in the drawings or modular structure suitable
Sequence perform either executed in parallel (environment of such as parallel processor or multiple threads, even including 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
It is intended to, specifically, as shown in figure 1, methods described can include：
S1：Seismic target earthquakes road in the pending geological data of selected frequency domain.
The pending geological data, refers to the geological data for carrying out surface wave dispersion process.
In the application one embodiment, the pending geological data is frequency domain.In another embodiment of the application
In, the pending geological data is time-domain, needs for the pending geological data of time-domain to be converted to treating for frequency domain
Process geological data.
The seismic target earthquakes road select, mainly according to the pending geological data, select the geological data in
The seismic channel that is located of each reference frequency point.
S2：Determine the adjacent seismic channel number threshold values and phase difference threshold values in the seismic target earthquakes road, the adjacent seismic channel
Number threshold values is equal to 2N, N >=1.
The adjacent seismic channel number threshold values and phase difference threshold values, enforcement personnel can incorporate experience into certainly according to actual conditions
Row determines.Wherein, the span of the phase difference threshold values is that 360 degree are less than or equal to more than or equal to 0 degree.
In the application one embodiment, the adjacent seismic channel number threshold values and phase difference for determining the seismic target earthquakes road
The mode of threshold values can include：
Surface wave dispersion process is carried out to the pending geological data, frequency-phase velocity frequency dispersion spectrum is obtained.
According to the frequency-phase velocity frequency dispersion spectrum, determine the seismic target earthquakes road adjacent seismic channel number threshold values and
Phase difference threshold values.
Typically, for the weaker situation of the noise of the frequency-phase velocity frequency dispersion spectrum, the adjacent seismic channel number valve
Value and phase difference threshold values can select smaller value.Conversely, for the stronger situation of the noise of the frequency-phase velocity frequency dispersion spectrum,
The adjacent seismic channel number threshold values and phase difference threshold values can select higher value.
Certainly, in the application other embodiment, the determination of specific adjacent seismic channel number threshold values and phase difference threshold values
Mode and determination standard, can voluntarily be determined by enforcement personnel according to actual conditions and operation experience.
S3：According to the adjacent seismic channel number threshold values, in the seismic target earthquakes road M adjacent seismic channel is chosen altogether, N≤
M≤2N。
The adjacent seismic channel, including the adjacent seismic channel in seismic target earthquakes road, the road number of described adjacent seismic channel is little
In the adjacent seismic channel number threshold values N, adjacent seismic channel number is added and is equal to the M.
If the seismic channel number in the seismic target earthquakes road is not less than N, respectively there are N number of adjacent seismic channel, now, M in both sides
=2N.
If the seismic channel number of seismic target earthquakes road side is 0, and the seismic channel of seismic target earthquakes road opposite side
Number is not less than N, then M=N.
The value of the M, determines according to the seismic channel number in actual seismic target earthquakes road, as long as arbitrarily the seismic channel number of 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 phase difference of the M adjacent seismic channels and the seismic target earthquakes road in each reference frequency point is calculated respectively.
Described each Frequency point, refers to the Frequency point included in the corresponding geological data in the seismic target earthquakes road.
S5：The phase difference is compared with the phase difference threshold values, the phase place on the reference frequency point is determined
Road number of the difference more than the seismic channel of the phase difference threshold values.
The reference frequency point includes the Frequency point included in the corresponding geological data in seismic target earthquakes road.
S6：Whether the comparison road number is more than default phase equalization threshold values, and according to comparative result corresponding ginseng is determined
Than the cumulative frequency point whether Frequency point is surface wave frequency dispersion curve.
The phase equalization threshold values is referred to and is more than the phase difference threshold values with the phase difference in the seismic target earthquakes road
The threshold values of seismic channel road number.
In the application one embodiment, whether the comparison road number is more than default phase equalization threshold values, according to comparing
As a result determine that whether corresponding reference frequency point is the mode of the cumulative frequency point of surface wave frequency dispersion curve, can include：
If the road number is more than default phase equalization threshold values, corresponding reference frequency point is not surface wave frequency dispersion curve
Cumulative frequency point.
If the road number is less than or equal to default phase equalization threshold values, corresponding reference frequency point is surface wave dispersion
The cumulative frequency point of curve.
Using a kind of embodiment of the data processing method of the surface wave frequency dispersion curve provided in the various embodiments described above, can be with
By determining that the corresponding seismic channel of geological data, in the phase equalization of each Frequency point, determines adjacent seismic channel in geological data
Phase place inconsistent (phase difference is larger) Frequency point.In geological data cumulative process, the phase of the adjacent seismic channel that only adds up
The corresponding geological data of Frequency point of position consistent (phase difference is little), can effectively eliminate making an uproar in frequency-phase velocity frequency dispersion spectrum
Part point, so as to effectively improve from the frequency-phase velocity frequency dispersion spectrum the precision for picking up the surface wave frequency dispersion curve for obtaining.
Fig. 2 is a kind of method stream of the data processing method of the surface wave frequency dispersion curve provided in the application another embodiment
Journey schematic diagram, specifically, can include：
S1：Seismic target earthquakes road in the pending geological data of selected frequency domain.
S2：Determine the adjacent seismic channel number threshold values and phase difference threshold values in the seismic target earthquakes road, the adjacent seismic channel
Number threshold values is equal to 2N, N >=1.
S3：According to the adjacent seismic channel number threshold values, in the seismic target earthquakes road M adjacent seismic channel is chosen altogether, N≤
M≤2N。
S4：The phase difference of the M adjacent seismic channels and the seismic target earthquakes road in each reference frequency point is calculated respectively.
S5：The phase difference is compared with the phase difference threshold values, determines that phase difference is more than institute on the reference frequency point
State the road number of the seismic channel of phase difference threshold values.
S6：Whether the comparison road number is more than default phase equalization threshold values, and according to comparative result corresponding ginseng is determined
Than the cumulative frequency point whether Frequency point is surface wave frequency dispersion curve.
S7：The corresponding geological data of the cumulative frequency point is added up, the geological data after being added up, to described
Geological data after cumulative carries out frequency dispersion process, obtains face wave frequency rate ,-phase velocity spectrum, from the face wave frequency rate ,-phase velocity spectrum
Pickup surface wave frequency dispersion curve.
The embodiment provided using above-described embodiment, the face wave frequency rate ,-phase velocity spectrum of the noise section that can be eliminated,
The higher surface wave frequency dispersion curve of the precision for obtaining can be picked up from the face wave frequency rate ,-phase velocity spectrum.
A kind of embodiment of the data processing method of the surface wave frequency dispersion curve provided using the various embodiments described above, Ke Yitong
The corresponding seismic channel of determination geological data is crossed in the phase equalization of each Frequency point, adjacent seismic channel in geological data is determined
The Frequency point of phase place inconsistent (phase difference is larger).In geological data cumulative process, the phase place of the adjacent seismic channel that only adds up
Unanimously the corresponding geological data of the Frequency point of (phase difference is little), can effectively eliminate the noise in frequency-phase velocity frequency dispersion spectrum
Part, so as to effectively improve from the frequency-phase velocity frequency dispersion spectrum the precision for picking up the surface wave frequency dispersion curve for obtaining.
In one example of the application, using herein described method, frequency dispersion point is carried out to geological data as shown in Figure 5
Analysis.
In Fig. 5, what abscissa was represented is seismic channel number, and what ordinate was represented is the time.
Fig. 6 is the face wave frequency rate ,-phase velocity spectrum not obtained using methods described in this example as a comparison.
Fig. 7 is in this example, to employ herein described method, and to the geological data in Fig. 5 dispersion analysis are carried out, and is obtained
A kind of face wave frequency rate ,-phase velocity spectrum.
During the corresponding dispersion analysis of Fig. 7, for 3, phase equalization threshold values is 1 to the adjacent seismic channel number threshold values for adopting.
Fig. 8 is in this example, to employ herein described method, and to the geological data in Fig. 5 dispersion analysis are carried out, and is obtained
A kind of face wave frequency rate ,-phase velocity spectrum.
During the corresponding dispersion analysis of Fig. 8, for 6, phase equalization threshold values is 1 to the adjacent seismic channel number threshold values for adopting.
In above-mentioned Fig. 6, Fig. 7, Fig. 8, the numeral on abscissa represents the numerical value of phase velocity, the generation of the numeral on ordinate
The numerical value of table frequency.
By the width figure of Fig. 6, Fig. 7, Fig. 8 tri- being contrasted it is found that face wave frequency rate ,-phase velocity for obtaining of Fig. 7 and Fig. 8
Distribution of color in spectrogram than Fig. 6 evenly, illustrates that noise has obtained preferable compacting.In addition, occurring the two of white in Fig. 8
Point, illustrates that wave energy strengthens, and signal to noise ratio is improved, more preferable than Fig. 7 to the pressing result of noise.
By the data processed result obtained in examples detailed above, it may be said that bright, using herein described method frequency dispersion is carried out
Analysis, can effectively improve the signal to noise ratio of face wave frequency rate ,-phase velocity spectrum, the face higher such that it is able to pick up the degree 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 one enforcement of the application
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 can include：
Seismic channel chosen module 101, can be used for selecting the seismic target earthquakes road in the pending geological data of frequency domain.
Threshold determination module 102, the threshold values being determined in pending geological data.
Phase difference calculating module 103, can be used for calculating respectively the M adjacent seismic channel and exists with the seismic target earthquakes road
The phase difference of each reference frequency point.
Road number determining module 104, can be used for being compared the phase difference with the phase difference threshold values, determine
Road number of the phase difference more than the seismic channel of the phase difference threshold values on the reference frequency point.
Whether Frequency point determining module 105, can be used for road number described in comparison more than default phase equalization threshold values, root
Determine that whether corresponding reference frequency point is the cumulative frequency point of surface wave frequency dispersion curve according to comparative result.
A kind of embodiment of the data processing equipment of the surface wave frequency dispersion curve provided using above-described embodiment, can be automatic
Perform methods described, it is not necessary to implement the excessive participation of personnel, effectively increase Consumer's Experience.
In the application another embodiment, the threshold determination module can include：
First threshold determination module, is determined for the adjacent seismic channel number threshold values and phase place in the seismic target earthquakes road
Difference limen value.
Second threshold determination module, is determined for phase equalization threshold values.
In the application another embodiment, first threshold determination module can include：
Frequency dispersion processing module, can be used for carrying out surface wave dispersion process to the pending geological data, obtain frequency-
Phase velocity frequency dispersion is composed；
Threshold values chooses module, can be used for, according to the frequency-phase velocity frequency dispersion spectrum, determining the seismic target earthquakes road
Adjacent seismic channel number threshold values and phase difference threshold values.
Fig. 4 is a kind of module knot of the data processing equipment of the surface wave frequency dispersion curve provided in the application another embodiment
Structure schematic diagram, specifically, described device can include：
Seismic channel chosen module 101, can be used for selecting the seismic target earthquakes road in the pending geological data of frequency domain.
Threshold determination module 102, the threshold values being determined in pending geological data.
Phase difference calculating module 103, can be used for calculating respectively the M adjacent seismic channel and exists with the seismic target earthquakes road
The phase difference of each reference frequency point.
Road number determining module 104, can be used for being compared the phase difference with the phase difference threshold values, determine
Road number of the phase difference more than the seismic channel of the phase difference threshold values on the reference frequency point.
Whether Frequency point determining module 105, can be used for road number described in comparison more than default phase equalization threshold values, root
Determine that whether corresponding reference frequency point is the cumulative frequency point of surface wave frequency dispersion curve according to comparative result.
Surface wave frequency dispersion curve generation module 106, for generating surface wave frequency dispersion curve.
In the application another embodiment, the surface wave frequency dispersion curve generation module can include：
Geological data accumulator module, for the corresponding geological data of cumulative frequency point to be added up, after being added up
Geological data.
Frequency dispersion processing module, for carrying out frequency dispersion process to the geological data after described adding up, obtains face wave frequency rate ,-phase velocity
Degree spectrum.
Dispersion curve pickup model, for picking up surface wave frequency dispersion curve from the face wave frequency rate ,-phase velocity spectrum.
A kind of data processing equipment of the surface wave frequency dispersion curve provided using the various embodiments described above, can be performed described automatically
Method, it is not necessary to implement the excessive participation of personnel, it is possible to obtain the surface wave frequency dispersion curve of degree of precision, effectively increase user
Experience.
In the data processing equipment of the surface wave frequency dispersion curve, the selected target seismic channel, the seismic target earthquakes are determined
The adjacent seismic channel number threshold values in road and phase difference threshold values, calculate phase difference, determine the phase difference on the reference frequency point
Count to more than the road of the seismic channel of the phase difference threshold values and determine whether corresponding reference frequency point is the tired of surface wave frequency dispersion curve
Plus the extension of the embodiment of Frequency point is 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 number threshold values and phase difference threshold values that determine the seismic target earthquakes road, calculate phase difference, determine described
Whether phase difference counts to the corresponding reference frequency point of determination more than the road of the seismic channel of the phase difference threshold values on reference frequency point
For the description of the various sequential manners of the cumulative frequency point of surface wave frequency dispersion curve, data acquisition/process/way of output etc., but,
The application is not limited to be professional standard or the situation described by embodiment etc., some professional standards or using making by oneself
Right way of conduct formula or embodiment description practice processes on embodiment amended slightly can also realize above-described embodiment it is identical, etc.
Foreseeable implementation result after same or close or deformation.Using these modifications or the data acquisition after deformation, process, export, sentencing
The embodiment of disconnected mode etc., within the scope of still may belong to 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 of the property made can include more or less operating procedures.The step of enumerating in embodiment order is only numerous steps
A kind of mode in execution sequence, does not represent unique execution sequence.When device or client production in practice is performed, can
So that according to embodiment, either method shown in the drawings order performs either executed in parallel (such as at parallel processor or multithreading
The environment of reason, even distributed data processing environment).Term " including ", "comprising" or its any other variant are intended to contain
Lid nonexcludability is included, so that a series of process, method, product or equipment including key elements not only will including those
Element, but also including other key elements being not expressly set out, or also include for this process, method, product or equipment
Intrinsic key element.In the absence of more restrictions, it is not excluded that including the process of the key element, method, product or
Also there are other identical or equivalent elements in person's equipment.
Device that above-described embodiment is illustrated or module etc., specifically can 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 with function when describing apparatus above to retouch respectively
State.Certainly, the function of each module can be realized in same or multiple softwares and/or hardware when the application is implemented,
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 have other drawing when actually realizing
Point mode, such as multiple module or components can with reference to 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 pure computer readable program code mode, it is complete
Entirely can by by method and step carry out programming in logic cause controller with gate, switch, special IC, may be programmed
The form of logic controller and embedded microcontroller etc. is realizing identical function.Therefore this controller is considered one kind
Hardware component, and the device for realizing various functions included to its inside can also be considered as the structure in hardware component.Or
Person even, can be used to realizing that the device of various functions be considered as not only being the software module of implementation method but also can be hardware
Structure in part.
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes execution particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure, class etc..The application can also be in a distributed computing environment put into practice, in these DCEs,
Task is performed by the remote processing devices connected by communication network.In a distributed computing environment, program module can
With positioned at including 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
Realize by the mode of software plus required general hardware platform.Based on such understanding, the technical scheme essence of the application
On prior art is contributed part in other words can be embodied in the form of software product, the computer software product
Can be stored in storage medium, such as ROM/RAM, magnetic disc, CD, including some instructions are used so that a computer equipment
(can be personal computer, mobile terminal, server, either network equipment etc.) perform each embodiment of the application or enforcement
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 was 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.For example：Personal computer, server computer, handheld device or
Portable set, laptop device, multicomputer system, set based on the system of microprocessor, set top box, programmable electronics
Standby, network PC, minicom, mainframe computer, including the DCE etc. of 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 deviating from spirit herein, it is desirable to which appended claim includes these deformations and changes without deviating from the application's
Spirit.
Claims (10)
1. a kind of data processing method of surface wave frequency dispersion curve, it is characterised in that methods described includes：
Seismic target earthquakes road in the pending geological data of selected frequency domain；
Determine the adjacent seismic channel number threshold values and phase difference threshold values in the seismic target earthquakes road, the adjacent seismic channel number threshold values
Equal to 2N, N >=1；
According to the adjacent seismic channel number threshold values, M adjacent seismic channel, N≤M≤2N are chosen altogether in the seismic target earthquakes road；
The phase difference of the M adjacent seismic channels and the seismic target earthquakes road in each reference frequency point is calculated respectively；
The phase difference is compared with the phase difference threshold values, determines that phase difference is more than institute on the reference frequency point
State the road number of the seismic channel of phase difference threshold values；
Whether the comparison road number is more than default phase equalization threshold values, and according to comparative result corresponding reference frequency point is determined
It is whether the cumulative frequency point of surface wave frequency dispersion curve.
2. a kind of data processing method of surface wave frequency dispersion curve as claimed in claim 1, it is characterised in that methods described is also wrapped
Include：
The corresponding geological data of the cumulative frequency point is added up, the geological data after being added up, after described adding up
Geological data carry out frequency dispersion process, obtain face wave frequency rate ,-phase velocity spectrum, pick up face from the face wave frequency rate ,-phase velocity spectrum
Wave dispersion curve.
3. a kind of data processing method of surface wave frequency dispersion curve as claimed in claim 1, it is characterised in that described in the comparison
Whether road number is more than default phase equalization threshold values, determines whether corresponding reference frequency point is face wave frequency according to comparative result
The mode of the cumulative frequency point of non-dramatic song line, including：
If the road number is more than default phase equalization threshold values, corresponding reference frequency point is not the tired of surface wave frequency dispersion curve
Plus Frequency point；
If the road number is less than or equal to default phase equalization threshold values, corresponding reference frequency point is surface wave frequency dispersion curve
Cumulative frequency point.
4. a kind of data processing method of the surface wave frequency dispersion curve as described in any one in claim 1 or 3, its feature exists
In the phase equalization threshold values includes：
With the threshold values of the seismic channel road number that the phase difference in the seismic target earthquakes road is more than the phase difference threshold values.
5. a kind of data processing method of surface wave frequency dispersion curve as claimed in claim 1, it is characterised in that described in the determination
The adjacent seismic channel number threshold values in seismic target earthquakes road and the mode of phase difference threshold values include：
Surface wave dispersion process is carried out to the pending geological data, frequency-phase velocity frequency dispersion spectrum is obtained；
According to the frequency-phase velocity frequency dispersion spectrum, the adjacent seismic channel number threshold values and phase place in the seismic target earthquakes road are determined
Difference limen value.
6. a kind of data processing equipment of surface wave frequency dispersion curve, it is characterised in that described device includes：
Seismic channel chosen module, for selecting the pending geological data of frequency domain in seismic target earthquakes road；
Threshold determination module, for determining pending geological data in threshold values；
Phase difference calculating module, for calculating the M adjacent seismic channel respectively with the seismic target earthquakes road in each reference frequency
The phase difference of rate point；
Road number determining module, for the phase difference to be compared with the phase difference threshold values, determines in reference frequency
Road number of the phase difference more than the seismic channel of the phase difference threshold values on rate point；
Whether Frequency point determining module, be more than default phase equalization threshold values, according to comparative result for relatively more described road number
Determine that whether corresponding reference frequency point is the cumulative frequency point of surface wave frequency dispersion curve.
7. a kind of data processing equipment of surface wave frequency dispersion curve as claimed in claim 6, it is characterised in that the threshold values determines
Module includes：
First threshold determination module, for determining the adjacent seismic channel number threshold values and phase difference threshold values in the seismic target earthquakes road；
Second threshold determination module, for determining phase equalization threshold values.
8. a kind of data processing equipment of surface wave frequency dispersion curve as claimed in claim 7, it is characterised in that first threshold values
Determining module includes：
Frequency dispersion processing module, for carrying out surface wave dispersion process to the pending geological data, obtains frequency-phase velocity frequency
Dissipate spectrum；
Threshold values chooses module, for according to the frequency-phase velocity frequency dispersion spectrum, determining the adjacent earthquake in the seismic target earthquakes road
Road number threshold values and phase difference threshold values.
9. a kind of data processing equipment of surface wave frequency dispersion curve as claimed in claim 6, it is characterised in that described device is also wrapped
Include：
Surface wave frequency dispersion curve generation module, for generating surface wave frequency dispersion curve.
10. a kind of data processing equipment of surface wave frequency dispersion curve as claimed in claim 9, it is characterised in that the face wave frequency
Non-dramatic song line generation module includes：
Geological data accumulator module, for the corresponding geological data of cumulative frequency point to be added up, the earthquake after being added up
Data；
Frequency dispersion processing module, for carrying out frequency dispersion process to the geological data after described adding up, obtains face wave frequency rate ,-phase velocity
Spectrum；
Dispersion curve pickup model, for picking up surface wave frequency dispersion curve from the face wave frequency rate ,-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 true CN106646615A (en) | 2017-05-10 |
| CN106646615B 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) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108562938A (en) * | 2018-03-23 | 2018-09-21 | 中国石油天然气股份有限公司 | Method, device and system for eliminating frequency dispersion effect |
| WO2019071515A1 (en) * | 2017-10-12 | 2019-04-18 | 南方科技大学 | Surface wave detection method and terminal device |
| CN109799530A (en) * | 2018-12-25 | 2019-05-24 | 核工业北京地质研究院 | Rayleigh waves dispersion curve inversion method for seismic surface wave exploration |
| CN111290017A (en) * | 2020-03-04 | 2020-06-16 | 南方科技大学 | Surface wave exploration method for jointly extracting Rayleigh wave frequency dispersion characteristics through seismic electric wave field |
| CN111766625A (en) * | 2020-07-06 | 2020-10-13 | 中国科学技术大学 | Seismic background noise dispersion curve extraction method based on deep learning |
| CN112051610A (en) * | 2020-10-21 | 2020-12-08 | 中国地质大学（北京） | Vector field multi-mode surface wave frequency dispersion calculation method and system |
| 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 |
| 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 |
Citations (5)
| 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 |
| US20160341839A1 (en) * | 2015-05-20 | 2016-11-24 | Conocophillips Company | Surface wave tomography using sparse data acquisition |
-
2016
- 2016-12-29 CN CN201611244008.8A patent/CN106646615B/en active Active
Patent Citations (5)
| 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 |
| US20160341839A1 (en) * | 2015-05-20 | 2016-11-24 | Conocophillips Company | Surface wave tomography using sparse data acquisition |
Non-Patent Citations (1)
| Title |
|---|
| 周腾飞 等: ""逆频散方法压制频散面波"", 《石油地球物理勘探》 * |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019071515A1 (en) * | 2017-10-12 | 2019-04-18 | 南方科技大学 | Surface wave detection method and terminal device |
| US10921477B2 (en) | 2018-03-23 | 2021-02-16 | Petrochina Company Limited | Method, apparatus, and system for eliminating frequency dispersion effect |
| CN108562938B (en) * | 2018-03-23 | 2019-09-06 | 中国石油天然气股份有限公司 | Method, device and system for eliminating frequency dispersion effect |
| CN108562938A (en) * | 2018-03-23 | 2018-09-21 | 中国石油天然气股份有限公司 | Method, device 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 |
| CN112748459B (en) * | 2019-10-29 | 2024-05-28 | 中国石油天然气集团有限公司 | Automatic frequency dispersion curve picking method and device |
| 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 |
| CN112987090B (en) * | 2019-12-02 | 2024-05-28 | 中国石油天然气集团有限公司 | Surface wave dispersion curve pickup method and device |
| 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 |
| CN111290017A (en) * | 2020-03-04 | 2020-06-16 | 南方科技大学 | Surface wave exploration method for jointly extracting Rayleigh wave frequency dispersion characteristics through seismic electric wave field |
| CN111766625A (en) * | 2020-07-06 | 2020-10-13 | 中国科学技术大学 | Seismic background noise dispersion curve extraction method based on deep learning |
| CN111766625B (en) * | 2020-07-06 | 2021-12-14 | 中国科学技术大学 | Seismic background noise dispersion curve extraction method based on deep learning |
| CN112051610A (en) * | 2020-10-21 | 2020-12-08 | 中国地质大学（北京） | Vector field multi-mode surface wave frequency dispersion calculation method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106646615B (en) | 2018-12-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106646615A (en) | Data processing method and data processing device for surface wave dispersion curve | |
| Bora et al. | On the relationship between Fourier and response spectra: Implications for the adjustment of empirical ground‐motion prediction equations (GMPEs) | |
| Callegari et al. | Spectral properties of chaos-based FM signals: Theory and simulation results | |
| CN104395779B (en) | System and method for estimating and weakening the noise in geological data | |
| CN110163827A (en) | Training method, image de-noising method, device and the medium of image denoising model | |
| CN105242309B (en) | The drawing method and device of rule interference in a kind of geological data | |
| CN107728210B (en) | Multiple instruments acquire the determination method and apparatus that road is lacked in data | |
| CN104932017B (en) | A kind of anomalous amplitude drawing method and device | |
| CN109061724A (en) | A kind of seismic data noise-reduction method based on adaptive variation mode decomposition | |
| CN107330128B (en) | Authentication abnormity judgment method and device | |
| CN104007315A (en) | Improved empirical mode decomposition processing method | |
| CN104599681B (en) | The method and apparatus of audio frequency process | |
| CN106680557A (en) | Electricity theft prevention method based on user behavior analysis | |
| CN109753985A (en) | Video classification methods and device | |
| CN109472758A (en) | A kind of seismic section image grain details Enhancement Method | |
| CN103763049A (en) | Cooperative spectrum sensing method based on FastICA algorithm | |
| Bodyanskiy et al. | On-line video segmentation using methods of fault detection in multidimensional time sequences | |
| CN109901222A (en) | Frequency method, storage medium and equipment are opened up in earthquake based on dynamic dictionary match tracing | |
| Salmasi et al. | Fractal dimension analysis for spike detection in low SNR extracellular signals | |
| CN108226636A (en) | Automatic filter method and apparatus | |
| Wang et al. | Time series sequences classification with inception and LSTM module | |
| CN102609469B (en) | Mining method for fuzzy rough monotonic data based on inclusion degree | |
| CN109633268A (en) | A kind of square wave fundamental frequency discrimination method based on B-spline and histogram | |
| CN106680879A (en) | Method and device for sedimentary cycle division | |
| Shin et al. | Transformation approach to model online gaming traffic |
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 |