CN105242309A - Earthquake data regular interference suppressing method and device - Google Patents
Earthquake data regular interference suppressing method and device Download PDFInfo
- Publication number
- CN105242309A CN105242309A CN201510677052.7A CN201510677052A CN105242309A CN 105242309 A CN105242309 A CN 105242309A CN 201510677052 A CN201510677052 A CN 201510677052A CN 105242309 A CN105242309 A CN 105242309A
- Authority
- CN
- China
- Prior art keywords
- interference
- geological data
- record
- data
- rule
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides an earthquake data regular interference suppressing method and device. The method comprises the steps that: seismic traces are arranged according to the geographical position ordering relation to form an earthquake data shot gather record to be processed; a regular noise spectrum analysis for an earthquake data area covered by the regular interference is performed to acquire the estimation basic frequency of the regular interference; the regular interference is separated to form a regular interference data record which is converted into an FK area earthquake data record; the front K inclination angles with the most powerful energy in earthquake data are determined by using energy scanning. In the regular interference data record, if N positive amplitude values of the inclination angles are continuously recorded, the determination is made that regular interference occurs in the inclination angles, and then filtering is performed in the inclination angle direction to form a new interference data record; a calculation is performed by subtracting the new interference data record from the earthquake data shot gather record to be processed. The earthquake data regular interference suppressing method and device can effectively reduce the valid signal damage in denoising of earthquake data regular interference so that the fidelity of the valid signals is improved after the noise suppression treatment.
Description
Technical field
The present invention relates to geophysical survey data processing technique, is reach by suppressing rule interference the method improving geological data signal to noise ratio (S/N ratio) and resolution in a kind of seismic data process.
Background technology
During geological data survey acquisition, the seismic signal that the wave detector that field is arranged receives, usually containing multiple perturbation wave, comprising can be used for the significant wave solving geological tasks, also comprises the interference wave affecting geological data imaging.
General, interference wave can be divided into regular interference and random disturbance ripple by the distribution characteristics according to interference wave.Regular interference can be called as coherent noise again, and this ripple usually has certain regularity on geological data spatio-temporal distribution.And random disturbance ripple is caused by random noise usually, the distribution on geological data on Time and place is without any rule.This two classes interference wave can affect geological data image quality usually in seismic data process, in order to improve geological data signal to noise ratio (S/N ratio) and resolution, often need to suppress this two classes interference wave, and pressing result also will directly affect geological data image quality.Therefore for the compacting of seismic data with low signal-to-noise ratio process noise, be especially an important step in seismic data process to the compacting of coherent noise.In field of seismic exploration, conventional relevant denoising drawing method comprises such as FK territory organized noise compact technique, KL transformation rule noise compact technique, t-p territory organized noise pressure technology, T-X territory and F-X territory organized noise compact technique etc., and these regular disturbance suppression technology can weaken useful signal energy to a certain extent, especially larger damage can be caused to low frequency useful signal, destroy useful signal even sometimes, cause to the geological data of post-processed not accurate enough, even bring the loss that cannot estimate, be not suitable with the current low frequency exploration engineering developed.
The compact technique of organized noise of the prior art especially can cause larger damage to low frequency effective information to the useful signal in geological data usually; cause to the geological data of post-processed not accurate enough; even bring the loss that cannot estimate, be not suitable with the current low frequency exploration engineering developed.
Summary of the invention
The object of the invention is to provide drawing method and the device of rule interference in a kind of geological data, the useful signal damage during interference denoising of geological data rule can be effectively reduced, while the interference of compacting rule, retain the useful signal in data as much as possible, improve the fidelity of the useful signal after noise compression process.
Drawing method and the device of rule interference in a kind of geological data provided by the invention are achieved in that
A drawing method for rule interference in geological data, described method comprises:
Seismic trace in the geological data of collection is arranged according to geographic position ordinal relation, forms pending geological data big gun collection record;
To being judged as in described pending geological data big gun collection record that organized noise spectrum analysis is carried out in the geological data region that there is rule interference covering, obtain described regular interference assessment dominant frequency;
From described pending geological data big gun collection record, the rule interference that frequency is described estimation dominant frequency is isolated, formation rule interfering data record by the wave filter of design; FK conversion is carried out to described regular interfering data record, obtains FK territory seismic data recording;
Within the scope of the rule interference inclination angle preset, energy scan is carried out to described FK territory geological data, determine front K the inclination angle that in geological data, energy is the strongest, K >=1;
In described regular interfering data record, respectively geological data amplitude statistics is carried out to a described K inclination angle; If count on N number of positive amplitude in described inclination angle continuously, then judge the geological data existence rule interference that described inclination angle is corresponding, N >=2;
In described pending geological data big gun collection record, judge that the inclination direction that there is rule interference carries out filtering along described, extract the regular interference on the described inclination direction that there is rule interference, form new interfering data record;
Described new interfering data record is deducted, the geological data after formation rule disturbance suppression with described pending geological data big gun collection record.
In preferred enforcement, described method also comprises: utilize the bandpass filter presetting frequency range to carry out filtering to described new interfering data record, interfering data record after acquisition filtering; Described default frequency range comprises the band limits that there is useful signal of thinking of user's setting of acquisition;
Accordingly, describedly deduct described new interfering data record with described pending geological data big gun collection record and comprise: deduct interfering data record after described filtering with described pending geological data big gun collection record.
In preferred enforcement, the span of described N is: 5≤N≤32.
In preferred enforcement, the value of described N is: 16.
In preferred enforcement, the span stating K is: 1≤K≤100.
In preferred enforcement, the value of described K is: 20.
A pressure setting for rule interference in geological data, described device comprises:
Sequence processing module, for being arranged according to geographic position ordinal relation by the seismic trace in the geological data of collection, forms pending geological data big gun collection record;
Interference dominant frequency estimation block, for being judged as in described pending geological data big gun collection record that organized noise spectrum analysis is carried out in the geological data region that there is rule interference covering, obtains described regular interference assessment dominant frequency;
Interference separation and conversion module, for being isolated the rule interference that frequency is described estimation dominant frequency from described pending geological data big gun collection record by the wave filter of design, formation rule interfering data record, and FK conversion is carried out to described regular interfering data record, obtain FK territory seismic data recording;
Scan module, for obtaining the K value of setting, and, within the scope of the rule interference inclination angle preset, energy scan is carried out to described FK territory geological data, determine front K the inclination angle that in geological data, energy is the strongest, K >=1;
Amplitude statistical module, for obtaining the N value of setting, and, in described regular interfering data record, respectively geological data amplitude statistics is carried out to a described K inclination angle; If count on N number of positive amplitude in described inclination angle continuously, then judge the geological data existence rule interference that described inclination angle is corresponding, N >=2;
Along described, dip filtering module, in described pending geological data big gun collection record, judges that the inclination direction that there is rule interference carries out filtering, extracts the regular interference on the described inclination direction that there is rule interference, form new interfering data record;
Interference processing module, for deducting described new interfering data record, the geological data after formation rule disturbance suppression with described pending geological data big gun collection record.
In preferred enforcement, described device also comprises: regular interfering data filtration module, for utilizing the bandpass filter of default frequency range to carry out filtering to described new interfering data record, and interfering data record after acquisition filtering; Described default frequency range comprises the band limits that there is useful signal of thinking of user's setting of acquisition;
Accordingly, deduct described new interfering data record with described pending geological data big gun collection record described in described interference processing module to comprise: deduct interfering data record after described filtering with described pending geological data big gun collection record.
In preferred enforcement, the span of described N is: 5≤N≤32.
In preferred enforcement, the value of described N is: 16.
In preferred enforcement, the span of described K is: 1≤K≤100.
In preferred enforcement, the value of described K is: 20.
The drawing method of rule interference in a kind of geological data provided by the invention and device, can by the geographic position relation rearrangement of the seismic trace in geological data according to seismic trace during image data, the linear distribution rule of rule interference can be kept like this, increase the coherence of rule interference, make it possible to effectively, accurately identify interfere information.Then, the present invention according to some differences of the aspects such as the inclination angle of regular interference noise and useful signal and energy, can judge the inclination angle of rule interference on big gun record in FK territory.Further can according to the noise inclination angle judged in FK territory, then the feature of rule interference lineups crest continuous distribution in time domain is adopted, again identify in time domain and finally determine regular interference, and obtaining relatively accurate reliable regular interfering data record by filtering.Then the method deducting regular interfering data record from original big gun record can be adopted to reach the object of compacting noise.The inventive method can improve geological data signal to noise ratio (S/N ratio) and resolution, and compared with FK territory linear voice denoising method simple in prior art, after denoising, useful signal fidelity is higher, can establish solid foundation for later stage geological data fidelity process.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, the accompanying drawing that the following describes is only some embodiments recorded in the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is the schematic flow sheet of a kind of embodiment of drawing method of rule interference in a kind of geological data of the present invention;
Fig. 2 is the schematic flow sheet of the another kind of embodiment of drawing method of rule interference in a kind of geological data of the present invention;
Fig. 3 is original single big gun seismogram before regular disturbance suppression in a kind of application scenarios of the present invention;
Fig. 4 is the regular interference big gun record utilizing small echo frequency division to leach from Fig. 1 in a kind of application scenarios of the present invention;
Fig. 5 is interfering data record after the filtering judged data shown in Fig. 4 by the inventive method;
Fig. 6 is the big gun record result that single primitively big gun seismologic record big gun record of Fig. 3 deducts after the filtering of Fig. 5 after interfering data record formation rule disturbance suppression;
Fig. 7 uses conventional FK filtering to Fig. 1 data rule disturbance suppression result;
Fig. 8 is the noise files that described conventional FK filter rule disturbance suppression exports;
Fig. 9 is the spectral contrast figure of Fig. 1, Fig. 6 and Fig. 7 data record;
Figure 10 is the modular structure schematic diagram of a kind of embodiment of pressure setting of rule interference in geological data of the present invention;
Figure 11 is the modular structure schematic diagram of the another kind of embodiment of pressure setting of rule interference in geological data of the present invention.
Embodiment
Technical scheme in the present invention is understood better in order to make those skilled in the art person, below in conjunction with the accompanying drawing in the embodiment of the present invention, technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, should belong to the scope of protection of the invention.
Be described in detail below in conjunction with the drawing method of accompanying drawing to rule interference in geological data of the present invention.Fig. 1 is the method flow diagram of a kind of embodiment of drawing method of rule interference in a kind of geological data of the present invention.Although the invention provides as following embodiment or method operation steps shown in the drawings, based on conventional or more or less operation steps can be comprised in the process without the need to performing creative labour.Logicality does not exist in necessary causal step, the execution sequence of these steps is not limited to the execution sequence that the embodiment of the present invention provides.In the data migration process in practice of described method or when device performs, can perform or executed in parallel (environment of such as parallel processor or multiple threads) according to embodiment or method shown in the drawings order.
In concrete provided by the invention a kind of geological data, as shown in Figure 1, described method can comprise a kind of embodiment of the drawing method of rule interference:
S1: arranged according to geographic position ordinal relation by the seismic trace in the geological data of collection, forms pending geological data big gun collection record.
The geological data of field acquisition can comprise the data according to certain way sequences such as road sequence or sequential usually.After obtaining acquiring seismic data in the present invention, can rearrange according to the geographic position ordinal relation of seismic trace in geological data, form pending geological data big gun collection record.Concrete, such as can arrange etc. according to seismic trace geophone offset order from small to large in survey line.Certainly, those skilled in the art also can take other mode to arrange according to seismic trace geographic position order, to realize the coherence strengthening regular interference noise, keep the linear distribution rule of rule interference, are convenient to identification and the prediction of regular interference noise.
S2: to being judged as in described pending geological data big gun collection record that organized noise spectrum analysis is carried out in the geological data region that there is rule interference covering, obtain described regular interference assessment dominant frequency.
Can judge to find the geological data region that there is rule interference and cover in the present invention in described pending geological data.Described determination methods can comprise the predetermined algorithm of employing or data processing method confirms to draw to after described pending geological data calculation process, such as, based on to the apparent velocity of earthquake data estimation or the amplitude of geological data or coherence calculation etc.In one embodiment of the present invention, the rule interference that also rule of thumb can be identified on described pending geological data big gun collection record by operating personnel and the existence rule interference distributing and then determine.
Then, organized noise spectrum analysis is carried out in the geological data region that can cover in described rule interference, obtains described regular interference assessment dominant frequency.The concrete the dominant frequency such as rule estimated can disturbed is as described regular interference assessment dominant frequency.
S3: isolate the rule interference that frequency is described estimation dominant frequency by the wave filter of design from described pending geological data big gun collection record, formation rule interfering data record, FK conversion is carried out to described regular interfering data record, obtains FK territory seismic data recording.
Can design a wavelet transform filter according to the estimation dominant frequency of the noise estimated in the present embodiment, the dominant frequency of filtering can be identical with described estimation dominant frequency, or in the predetermined threshold range centered by described estimation dominant frequency.The present embodiment can carry out frequency-division filter process by small echo frequency division to geological data, by frequency division, rule interference is separated from big gun record, forms one based on the regular interfering data record of rule interference.Then can do FK conversion to this regular interfering data record, form FK territory seismic data recording.
S4: within the scope of the rule interference inclination angle preset, energy scan is carried out to described FK territory geological data, determine front K the inclination angle that in geological data, energy is the strongest, K >=1.
A rule interference inclination angle scope can be set in the present invention.The apparent velocity scope that can comprise the regular interferencing propagation determined of the rule interference inclination angle half described in the application.Such as, described apparent velocity can be set in a kind of embodiment of the present invention select 500m/s-2000m/s, within the scope of 500m/s-2000m/s, inclination angle energy scan be carried out to the FK territory geological data formed, determine front K the inclination angle that in geological data, energy is the strongest.In a kind of embodiment of the present invention, the span of described inclination angle number K can comprise 1-00.Concrete can be arranged according to primary treatment demand.In a kind of embodiment of the present invention, the value of described K can be 20, namely can scan 20 inclination angles that energy is the strongest.Then the size at these inclination angles can be recorded.
S5: respectively geological data amplitude statistics is carried out to a described K inclination angle in described regular interfering data record; If count on N number of positive amplitude in described inclination angle continuously, then judge the geological data existence rule interference that described inclination angle is corresponding, N >=2.
A kind of embodiment of the present invention, can in the regular interfering data record of small echo frequency division filtering, geological data amplitude statistics is carried out at 20 inclination angles the strongest along energy, carries out noise compacting, be just left intact on the contrary to the inclination angle counting on the individual positive amplitude of K continuously.
The value of the N described in the present invention can carry out arranging, determining according to the concrete distribution situation of rule interference.General, the span of the N described in the present invention can comprise: 5≤N≤32.N described in the present embodiment can value be 16, namely carries out noise compacting to the inclination angle counting on 16 positive amplitudes continuously.
Along described, S6: in described pending geological data big gun collection record, judges that the inclination direction that there is rule interference carries out filtering, extracts the regular interference on the described inclination direction that there is rule interference, form new interfering data record.
Noise compacting can be carried out to judging the data existed on the inclination direction of rule interference.In an embodiment of the present invention, concrete drawing method can comprise: according to the rule interference inclination angle judged and the noise dominant frequency estimated, design a wave filter, filtering is carried out along the rule interference inclination direction judged in noise big gun record, extract all regular interferences existed on rule interference inclination angle, form a new interfering data record.
S7: deduct described new interfering data record, the geological data after formation rule disturbance suppression with described pending geological data big gun collection record.
After utilizing processing mode described above to obtain described new interfering data record, the new interfering data record that the embodiment of the present invention obtains can be deducted with the described pending geological data big gun record before denoising, the noise compacting of data record can be completed.
In a preferred embodiment of the present invention, after the described new interfering data record of acquisition, can also further screen in described new interfering data record, to filter the useful signal that may exist, and then the infringement reduced useful signal, the fidelity of useful signal after raising data processing.Fig. 2 is the schematic flow sheet of the another kind of embodiment of drawing method of rule interference in a kind of geological data of the present invention, and as shown in Figure 2, in another kind of embodiment of the present invention, described method can also comprise:
S701: utilize the bandpass filter presetting frequency range to carry out filtering to described new interfering data record, interfering data record after acquisition filtering; Described default frequency range comprises the band limits that there is useful signal of thinking of user's setting of acquisition;
Accordingly, describedly deduct described new interfering data record with described pending geological data big gun collection record and comprise: deduct interfering data record after described filtering with described pending geological data big gun collection record.
In concrete implementation, can to the new regular interfering data record of above-mentioned middle acquisition, bamboo product bandpass filter, user can determine the scope of described bandpass filtering with reference to described regular interference assessment dominant frequency.The useful signal that also may exist in described new interfering data record thought by this wave filter mainly filtering user, such as, can be 1Hz, 5Hz, 20Hz or 30Hz.
The drawing method of rule interference in the geological data that the embodiment of the present invention provides, can by the geographic position relation rearrangement of the seismic trace in geological data according to seismic trace during image data, the linear distribution rule of rule interference can be kept like this, increase the coherence of rule interference, make it possible to effectively, accurately identify interfere information.The inventive method not only can improve geological data signal to noise ratio (S/N ratio) and resolution, and compared with FK territory linear voice denoising method simple in prior art, after denoising, useful signal fidelity is higher, can establish solid foundation for later stage geological data fidelity process.
Here is the embodiment that the method for the invention applies to that actual seismic data single shot record carries out an application scenarios of linear voice decay.
Fig. 3 is original single big gun seismogram before regular disturbance suppression in application scenarios.Original single big gun seismologic record shown in Fig. 3 sorts by geographic position, the display of Fig. 3 horizontal ordinate direction be seismic trace Taoist monastic name, the display of ordinate direction be sampling point time of each seismic trace.Regular interference noise on estimation big gun record and distributed areas, then carry out organized noise spectrum analysis in the geological data region that regular interference noise covers, in Fig. 3, roughly the dominant frequency of estimated noise is 10Hz.Design a wavelet transform filter, the dominant frequency of filtering is 10Hz, by small echo frequency division, rule interference is separated from big gun record, form one based on the new noise big gun record of rule interference, the regular interference big gun record utilizing small echo frequency division to leach as shown in Figure 4.As can be seen from Figure 4 also there is useful signal in the regular interference big gun record that small echo frequency division is formed.Therefore, FK conversion can be done to these noise big gun record data in this application scene, form FK territory geological data.A given rule interference inclination angle scope, be also generally the apparent velocity scope of regular interferencing propagation, the present embodiment can select 500m/s-2000m/s, within the scope of 500m/s-2000m/s, carry out inclination angle energy scan to the FK data formed.The present embodiment has scanned 20 the strongest inclination angles of energy, and the size at these inclination angles is recorded, and in the rule interference big gun record of small echo frequency division filtering, geological data amplitude statistics is carried out at 20 inclination angles the strongest along energy.Noise compacting is carried out to the inclination angle counting on 16 positive amplitudes continuously, is just left intact on the contrary.Concrete noise drawing method can comprise: according to the rule interference inclination angle judged and the noise dominant frequency estimated, design a wave filter, filtering is carried out along the rule interference inclination direction judged in noise big gun record, extract all regular interferences existed on rule interference inclination angle, form a new regular interference big gun record.To the new regular interference big gun record obtained, bamboo product free transmission range is (1Hz, 5Hz, 20Hz, bandpass filter 30Hz), carries out filtering process by this bandpass filter to new regular interference big gun record, obtains interfering data record after filtering.Fig. 5 is interfering data record after the filtering judged data shown in Fig. 4 by the inventive method.Finally can deduct with the original geological data big gun record before denoising and invent interfering data record after the filtering more accurately that obtains, complete the noise compacting of single big gun seismic data recording.Fig. 6 is the big gun collection record result that single primitively big gun seismologic record big gun record of Fig. 3 deducts after the filtering of Fig. 5 after interfering data record formation rule disturbance suppression.
The drawing method of rule interference in geological data provided by the invention, the useful signal damage during interference denoising of geological data rule can be effectively reduced, while the interference of compacting rule, retain the useful signal in data as much as possible, improve the fidelity of the useful signal after noise compression process.Fig. 7 uses conventional FK filtering to Fig. 1 data rule disturbance suppression result, Fig. 8 is the noise files that described conventional FK filter rule disturbance suppression exports, comparison diagram 8 and Fig. 5 obviously can find out in Fig. 8 also there is part useful signal energy, and almost can't see useful signal in Fig. 5.Simultaneously; Fig. 9 is the spectral contrast figure of Fig. 1, Fig. 6 and Fig. 7 data record; from Fig. 9, (line_1 is Fig. 1 big gun record frequency spectrum before denoising; line_2 is Fig. 6 big gun record frequency spectrum after the regular disturbance suppression of utilization invention; line_3 be conventional FK territory rule disturbance suppression after Fig. 7 big gun record frequency spectrum) contrast can find; at about 10Hz; the inventive method more can protect useful signal energy than conventional FK denoising; also almost completely the same at the frequency spectrum of HFS and original earthquake data after the present invention's rule disturbance suppression, and there is the difference of any in conventional FK filtering.This can illustrate that the drawing method of rule interference in geological data provided by the invention can improve regular interference prediction precision, improves signal to noise ratio (S/N ratio) and the resolution of geological data to a great extent.Simultaneously carry out noise attentuation owing to adopting in time domain, compared with suppressing with conventional FK territory organized noise, more can keep the waveform of geological data after regular disturbance suppression, after avoiding denoising, produce alias.
Based on the drawing method of rule interference in geological data of the present invention, the invention provides the pressure setting of rule interference in a kind of geological data.Figure 10 is the modular structure schematic diagram of a kind of embodiment of pressure setting of rule interference in geological data of the present invention, and as shown in Figure 10, described device can comprise:
Sequence processing module 101, may be used for the seismic trace in the geological data of collection to arrange according to geographic position ordinal relation, forms pending geological data big gun collection record;
Interference dominant frequency estimation block 102, may be used for, to being judged as in described pending geological data big gun collection record that organized noise spectrum analysis is carried out in the geological data region that there is rule interference covering, obtaining described regular interference assessment dominant frequency;
Interference separation and conversion module 103, may be used for from described pending geological data big gun collection record, isolating by the wave filter of design the rule interference that frequency is described estimation dominant frequency, formation rule interfering data record, and FK conversion is carried out to described regular interfering data record, obtain FK territory seismic data recording;
Scan module 104, may be used for obtaining the K value arranged, and, within the scope of the rule interference inclination angle preset, energy scan is carried out to described FK territory geological data, determine front K the inclination angle that in geological data, energy is the strongest, K >=1;
Amplitude statistical module 105, may be used for obtaining the N value arranged, and, in described regular interfering data record, respectively geological data amplitude statistics is carried out to a described K inclination angle; If count on N number of positive amplitude in described inclination angle continuously, then judge the geological data existence rule interference that described inclination angle is corresponding, N >=2;
Dip filtering module 106, may be used in described pending geological data big gun collection record, judge that the inclination direction that there is rule interference carries out filtering along described, extract the regular interference on the described inclination direction that there is rule interference, form new interfering data record;
Interference processing module 107, may be used for deducting described new interfering data record, the geological data after formation rule disturbance suppression with described pending geological data big gun collection record.
Figure 11 is the modular structure schematic diagram of the another kind of embodiment of pressure setting of rule interference in geological data of the present invention, and as shown in figure 11, described device can also comprise:
Rule interfering data filtration module 108, may be used for utilizing the bandpass filter of default frequency range to carry out filtering to described new interfering data record, interfering data record after acquisition filtering; Described default frequency range comprises the band limits that there is useful signal of thinking of user's setting of acquisition;
Accordingly, deduct described new interfering data record with described pending geological data big gun collection record described in described interference processing module 107 to comprise: deduct interfering data record after described filtering with described pending geological data big gun collection record.
The specific descriptions that in device of the present invention, modules relates to can with reference to the description of method described above, and such as, the span of described N can comprise: 5≤N≤32, and the value of N described in preferred embodiment can be: 16.Or in other embodiment, the span of described K can comprise: 1≤K≤100, in preferred embodiment, the value of described K can be: 20 etc., between each embodiment can mutually with reference to, do not repeat at this.Described device can be understood as the implementation of a kind of product realizing the method for the invention.
The invention provides drawing method and the device of rule interference in a kind of geological data, small echo frequency division first burbling noise and useful signal can be passed through, and maintain the waveform character of regular interference noise to a certain extent, converted by FK, the organized noise energy at identical inclination angle is stacked up in FK territory, enhance the energy of rule interference, be convenient to judge rule interference inclination angle and identify in FK territory, then utilize these inclination angles to carry out noise in time domain again judge and identify, by carrying out bandpass filtering to noise, better protect the energy of useful signal.Compared with conventional FK rule interference filter, the waveform character of geological data after denoising can be kept better.
Although mention the data transformation of data filtering, FK conversion, energy scan, amplitude statistics etc., the description of process in content of the present invention, it must be the data transformation of complete standard or mentioned mode, the situation of process that the present invention is not limited to.Foregoing description in the present invention involved by each embodiment is only the application in some embodiments in the present invention, and on the basis of some standard, method, amended slightly disposal route also can carry out the scheme of each embodiment of the invention described above.Certainly, meet the process method step described in the various embodiments described above of the present invention other without creationary distortion, still can realize identical application, not repeat them here.
Although the invention provides the method operation steps as described in embodiment or process flow diagram, based on conventional or more or less operation steps can be comprised without performing creative labour.The sequence of steps enumerated in embodiment is only a kind of mode in numerous step execution sequence, does not represent unique execution sequence.When device in practice or client production perform, can perform or executed in parallel (environment of such as parallel processor or multiple threads) according to embodiment or method shown in the drawings order.
The device that above-described embodiment is illustrated or module, specifically can be realized by computer chip or entity, or be realized by the product with certain function.For convenience of description, various module is divided into describe respectively with function when describing above device.Certainly, the function of each module can be realized in same or multiple software and/or hardware when implementing of the present invention.Certainly, also the module realizing certain function can be combined realization by multiple submodule or subelement.
Method described in the present invention, device or module can realize by any suitable mode with the controller of computer readable program code, such as, controller can be taked such as microprocessor or processor and store the computer-readable medium of the computer readable program code (such as software or firmware) that can be performed by this (micro-) processor, logic gate, switch, special IC (ApplicationSpecificIntegratedCircuit, ASIC), the form of programmable logic controller (PLC) and embedding microcontroller, the example of controller includes but not limited to following microcontroller: ARC625D, AtmelAT91SAM, MicrochipPIC18F26K20 and SiliconeLabsC8051F320, Memory Controller can also be implemented as a part for the steering logic of storer.Those skilled in the art also know, except realizing except controller in pure computer readable program code mode, controller can be made to realize identical function with the form of logic gate, switch, special IC, programmable logic controller (PLC) and embedding microcontroller etc. by method step being carried out programming in logic completely.Therefore this controller can be considered to a kind of hardware component, and to the structure that also can be considered as the device realizing various function in hardware component that its inside comprises.Or even, the device being used for realizing various function can be considered as not only can be implementation method software module but also can be structure in hardware component.
Part of module in device of the present invention can describe in the general context of computer executable instructions, such as program module.Usually, program module comprises the routine, program, object, assembly, data structure, class etc. that perform particular task or realize particular abstract data type.Also can put into practice the present invention in a distributed computing environment, in these distributed computing environment, be executed the task by the remote processing devices be connected by communication network.In a distributed computing environment, program module can be arranged in the local and remote computer-readable storage medium comprising memory device.
As seen through the above description of the embodiments, those skilled in the art can be well understood to the mode that the present invention can add required hardware by software and realizes.Based on such understanding, technical scheme of the present invention can embody with the form of software product the part that prior art contributes in essence in other words, also can by embodying in the implementation process of Data Migration.This computer software product can be stored in storage medium, as ROM/RAM, magnetic disc, CD etc., comprising some instructions in order to make a computer equipment (can be personal computer, mobile terminal, server, or the network equipment etc.) perform the method described in some part of each embodiment of the present invention or embodiment.
Each embodiment in this instructions adopts the mode of going forward one by one to describe, between each embodiment same or analogous part mutually see, what each embodiment stressed is the difference with other embodiments.Of the present invention can be used in whole or in part numerous general or special purpose computing system environments or configuration in.Such as: personal computer, server computer, handheld device or portable set, laptop device, mobile communication terminal, multicomputer system, the system based on microprocessor, programmable electronic equipment, network PC, small-size computer, mainframe computer, the distributed computing environment comprising above any system or equipment etc.
Although depict the present invention by embodiment, those of ordinary skill in the art know, the present invention has many distortion and change and do not depart from spirit of the present invention, and the claim appended by wishing comprises these distortion and change and do not depart from spirit of the present invention.
Claims (12)
1. a drawing method for rule interference in geological data, it is characterized in that, described method comprises:
Seismic trace in the geological data of collection is arranged according to geographic position ordinal relation, forms pending geological data big gun collection record;
To being judged as in described pending geological data big gun collection record that organized noise spectrum analysis is carried out in the geological data region that there is rule interference covering, obtain described regular interference assessment dominant frequency;
From described pending geological data big gun collection record, the rule interference that frequency is described estimation dominant frequency is isolated, formation rule interfering data record by the wave filter of design; FK conversion is carried out to described regular interfering data record, obtains FK territory seismic data recording;
Within the scope of the rule interference inclination angle preset, energy scan is carried out to described FK territory geological data, determine front K the inclination angle that in geological data, energy is the strongest, K >=1;
In described regular interfering data record, respectively geological data amplitude statistics is carried out to a described K inclination angle; If count on N number of positive amplitude in described inclination angle continuously, then judge the geological data existence rule interference that described inclination angle is corresponding, N >=2;
In described pending geological data big gun collection record, judge that the inclination direction that there is rule interference carries out filtering along described, extract the regular interference on the described inclination direction that there is rule interference, form new interfering data record;
Described new interfering data record is deducted, the geological data after formation rule disturbance suppression with described pending geological data big gun collection record.
2. the drawing method of rule interference in a kind of geological data as claimed in claim 1, it is characterized in that, described method also comprises:
The bandpass filter presetting frequency range is utilized to carry out filtering to described new interfering data record, interfering data record after acquisition filtering; Described default frequency range comprises the band limits that there is useful signal of thinking of user's setting of acquisition;
Accordingly, describedly deduct described new interfering data record with described pending geological data big gun collection record and comprise: deduct interfering data record after described filtering with described pending geological data big gun collection record.
3. the drawing method of rule interference in a kind of geological data as claimed in claim 1 or 2, it is characterized in that, the span of described N is: 5≤N≤32.
4. the drawing method of rule interference in a kind of geological data as claimed in claim 3, it is characterized in that, the value of described N is: 16.
5. the drawing method of rule interference in a kind of geological data as claimed in claim 1 or 2, it is characterized in that, the span of described K is: 1≤K≤100.
6. the drawing method of rule interference in a kind of geological data as claimed in claim 5, it is characterized in that, the value of described K is: 20.
7. a pressure setting for rule interference in geological data, it is characterized in that, described device comprises:
Sequence processing module, for being arranged according to geographic position ordinal relation by the seismic trace in the geological data of collection, forms pending geological data big gun collection record;
Interference dominant frequency estimation block, for being judged as in described pending geological data big gun collection record that organized noise spectrum analysis is carried out in the geological data region that there is rule interference covering, obtains described regular interference assessment dominant frequency;
Interference separation and conversion module, for being isolated the rule interference that frequency is described estimation dominant frequency from described pending geological data big gun collection record by the wave filter of design, formation rule interfering data record, and FK conversion is carried out to described regular interfering data record, obtain FK territory seismic data recording;
Scan module, for obtaining the K value of setting, and, within the scope of the rule interference inclination angle preset, energy scan is carried out to described FK territory geological data, determine front K the inclination angle that in geological data, energy is the strongest, K >=1;
Amplitude statistical module, for obtaining the N value of setting, and, in described regular interfering data record, respectively geological data amplitude statistics is carried out to a described K inclination angle; If count on N number of positive amplitude in described inclination angle continuously, then judge the geological data existence rule interference that described inclination angle is corresponding, N >=2;
Along described, dip filtering module, in described pending geological data big gun collection record, judges that the inclination direction that there is rule interference carries out filtering, extracts the regular interference on the described inclination direction that there is rule interference, form new interfering data record;
Interference processing module, for deducting described new interfering data record, the geological data after formation rule disturbance suppression with described pending geological data big gun collection record.
8. the pressure setting of rule interference in a kind of geological data as claimed in claim 7, it is characterized in that, described device also comprises:
Rule interfering data filtration module, for utilizing the bandpass filter of default frequency range to carry out filtering to described new interfering data record, interfering data record after acquisition filtering; Described default frequency range comprises the band limits that there is useful signal of thinking of user's setting of acquisition;
Accordingly, deduct described new interfering data record with described pending geological data big gun collection record described in described interference processing module to comprise: deduct interfering data record after described filtering with described pending geological data big gun collection record.
9. the pressure setting of rule interference in a kind of geological data as described in claim 7 or 8, it is characterized in that, the span of described N is: 5≤N≤32.
10. the pressure setting of rule interference in a kind of geological data as claimed in claim 9, it is characterized in that, the value of described N is: 16.
The pressure setting of rule interference in 11. a kind of geological datas as claimed in claim 7 or 8, it is characterized in that, the span of described K is: 1≤K≤100.
The pressure setting of rule interference in 12. a kind of geological datas as claimed in claim 11, it is characterized in that, the value of described K is: 20.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510677052.7A CN105242309B (en) | 2015-10-19 | 2015-10-19 | The drawing method and device of rule interference in a kind of geological data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510677052.7A CN105242309B (en) | 2015-10-19 | 2015-10-19 | The drawing method and device of rule interference in a kind of geological data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105242309A true CN105242309A (en) | 2016-01-13 |
| CN105242309B CN105242309B (en) | 2017-11-10 |
Family
ID=55040009
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510677052.7A Active CN105242309B (en) | 2015-10-19 | 2015-10-19 | The drawing method and device of rule interference in a kind of geological data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105242309B (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107422382A (en) * | 2016-05-24 | 2017-12-01 | 中国石油化工股份有限公司 | The filtering method and device of self-adaption two-dimensional fitting of a polynomial |
| CN108363103A (en) * | 2018-02-01 | 2018-08-03 | 中国石油天然气集团有限公司 | Regular disturbance suppression method and device |
| CN111722286A (en) * | 2019-03-18 | 2020-09-29 | 中石化石油工程技术服务有限公司 | Random function median filtering technique |
| CN111781646A (en) * | 2020-06-29 | 2020-10-16 | 广州海洋地质调查局 | Amplitude-preserving and fidelity-preserving linear coherent noise suppression method and processing terminal |
| CN112083490A (en) * | 2019-06-12 | 2020-12-15 | 中国石油天然气集团有限公司 | Seismic data noise attenuation method and device |
| CN112799132A (en) * | 2019-11-13 | 2021-05-14 | 中国石油天然气股份有限公司 | Micro local linear noise suppression method and device |
| CN112882102A (en) * | 2019-11-29 | 2021-06-01 | 中国石油天然气集团有限公司 | Seismic noise attenuation method and device |
| CN113126164A (en) * | 2020-01-10 | 2021-07-16 | 中国石油天然气集团有限公司 | Seismic data denoising method and device |
| CN114740530A (en) * | 2021-01-07 | 2022-07-12 | 中国石油天然气股份有限公司 | Medium-high frequency quasi-linear noise suppression method and device based on hyperbolic time window constraint |
| CN115685320A (en) * | 2022-11-02 | 2023-02-03 | 中国矿业大学(北京) | Seismic shot gather data denoising method and device |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101930081A (en) * | 2009-06-26 | 2010-12-29 | 中国石油集团东方地球物理勘探有限责任公司 | Automatic identification and suppression method of single-frequency interference in seismic records |
| CN102081168A (en) * | 2009-12-01 | 2011-06-01 | 中国石油天然气集团公司 | Method for improving surface wave attenuation performance in seismic data processing |
| US20110276274A1 (en) * | 2010-05-06 | 2011-11-10 | Margaret Chen Yu | Seismic Interference Noise Elimination |
| AU2014201718A1 (en) * | 2013-03-22 | 2014-10-09 | Cgg Services Sa | Vector-dip filtering of seismic data in the time-frequency domain |
| CN104181600A (en) * | 2014-08-25 | 2014-12-03 | 中国石油天然气集团公司 | Seismic data linear noise attenuation method and device |
| CN104849758A (en) * | 2015-05-05 | 2015-08-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Regular disturbance suppression method for seismic data |
-
2015
- 2015-10-19 CN CN201510677052.7A patent/CN105242309B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101930081A (en) * | 2009-06-26 | 2010-12-29 | 中国石油集团东方地球物理勘探有限责任公司 | Automatic identification and suppression method of single-frequency interference in seismic records |
| CN102081168A (en) * | 2009-12-01 | 2011-06-01 | 中国石油天然气集团公司 | Method for improving surface wave attenuation performance in seismic data processing |
| US20110276274A1 (en) * | 2010-05-06 | 2011-11-10 | Margaret Chen Yu | Seismic Interference Noise Elimination |
| AU2014201718A1 (en) * | 2013-03-22 | 2014-10-09 | Cgg Services Sa | Vector-dip filtering of seismic data in the time-frequency domain |
| CN104181600A (en) * | 2014-08-25 | 2014-12-03 | 中国石油天然气集团公司 | Seismic data linear noise attenuation method and device |
| CN104849758A (en) * | 2015-05-05 | 2015-08-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Regular disturbance suppression method for seismic data |
Non-Patent Citations (3)
| Title |
|---|
| QIAN ZHONGPING,ET AL.: "A regular noise elimination method for prestack 3D seismic data", 《SEG TECHNICAL PROGRAM EXPANDED ABSTRACTS 2003》 * |
| 夏洪瑞等: "均值加权消除相干干扰", 《石油物探》 * |
| 甘其刚等: "叠前时空域线性干扰的衰减及应用", 《石油物探》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107422382A (en) * | 2016-05-24 | 2017-12-01 | 中国石油化工股份有限公司 | The filtering method and device of self-adaption two-dimensional fitting of a polynomial |
| CN108363103A (en) * | 2018-02-01 | 2018-08-03 | 中国石油天然气集团有限公司 | Regular disturbance suppression method and device |
| CN108363103B (en) * | 2018-02-01 | 2019-10-11 | 中国石油天然气集团有限公司 | Regular disturbance suppression method and device |
| CN111722286A (en) * | 2019-03-18 | 2020-09-29 | 中石化石油工程技术服务有限公司 | Random function median filtering technique |
| CN111722286B (en) * | 2019-03-18 | 2022-12-30 | 中石化石油工程技术服务有限公司 | Random function median filtering technique |
| CN112083490A (en) * | 2019-06-12 | 2020-12-15 | 中国石油天然气集团有限公司 | Seismic data noise attenuation method and device |
| CN112799132B (en) * | 2019-11-13 | 2023-08-22 | 中国石油天然气股份有限公司 | Micro-local linear noise suppression method and device |
| CN112799132A (en) * | 2019-11-13 | 2021-05-14 | 中国石油天然气股份有限公司 | Micro local linear noise suppression method and device |
| CN112882102A (en) * | 2019-11-29 | 2021-06-01 | 中国石油天然气集团有限公司 | Seismic noise attenuation method and device |
| CN112882102B (en) * | 2019-11-29 | 2024-05-28 | 中国石油天然气集团有限公司 | Method and device for attenuating seismic noise |
| CN113126164A (en) * | 2020-01-10 | 2021-07-16 | 中国石油天然气集团有限公司 | Seismic data denoising method and device |
| CN113126164B (en) * | 2020-01-10 | 2024-04-30 | 中国石油天然气集团有限公司 | Seismic data denoising method and device |
| CN111781646A (en) * | 2020-06-29 | 2020-10-16 | 广州海洋地质调查局 | Amplitude-preserving and fidelity-preserving linear coherent noise suppression method and processing terminal |
| CN111781646B (en) * | 2020-06-29 | 2022-04-15 | 广州海洋地质调查局 | Amplitude-preserving and fidelity-preserving linear coherent noise suppression method and processing terminal |
| CN114740530A (en) * | 2021-01-07 | 2022-07-12 | 中国石油天然气股份有限公司 | Medium-high frequency quasi-linear noise suppression method and device based on hyperbolic time window constraint |
| CN115685320A (en) * | 2022-11-02 | 2023-02-03 | 中国矿业大学(北京) | Seismic shot gather data denoising method and device |
| CN115685320B (en) * | 2022-11-02 | 2024-01-26 | 中国矿业大学(北京) | Method and device for denoising seismic shot gather data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105242309B (en) | 2017-11-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105242309A (en) | Earthquake data regular interference suppressing method and device | |
| Han et al. | Microseismic and seismic denoising via ensemble empirical mode decomposition and adaptive thresholding | |
| GB2500302A (en) | Removing non-orthogonal coherent noise from seismic data | |
| CN104007469A (en) | Weak seismic signal reconstruction method based on curvelet transform | |
| CN105116442A (en) | Lithologic oil-gas reservoir weak-reflection seismic signal reconstruction method | |
| CN104932017A (en) | Anomalous amplitude suppression method and apparatus | |
| Zhang et al. | Signal preserving and seismic random noise attenuation by Hurst exponent based time–frequency peak filtering | |
| CN102944895B (en) | Suppression method and suppression device for surface wave interference for seismic record | |
| Honório et al. | Progress on empirical mode decomposition-based techniques and its impacts on seismic attribute analysis | |
| CN107561588B (en) | Seismic data noise suppression method and device | |
| CN105093282A (en) | Energy replacement surface wave suppressing method based on frequency constraint | |
| Wu et al. | Multi-scale progressive fusion attention network based on small sample training for DAS noise suppression | |
| CN107255833B (en) | The determination method and apparatus of Seismic Stacked Section | |
| CN108693558B (en) | Seismic data processing method and device | |
| CN105301656A (en) | Method and device for suppressing ghost reflection signals | |
| Liao et al. | 3-D Butterworth filtering for 3-D high-density onshore seismic field data | |
| CN104570118A (en) | Method for automatically identifying and removing industrial interface based on dual factors | |
| CN112782763A (en) | Method, device and equipment for estimating seismic quality factor and storage medium | |
| CN105425297A (en) | Method and device of suppressing ghost signals | |
| CN113050167B (en) | Noise suppression method and device for seismic data | |
| CN109116425A (en) | Utilize the method for the frequency spectrum design filter removal noise of back wave | |
| CN112083491B (en) | Wellbore wave pressing method and device based on signal attribute characteristics | |
| CN113640882B (en) | Acquisition footprint noise removal method, electronic device and computer readable storage medium | |
| CN112230286B (en) | Method and system for reducing noise of tunnel wireless seismograph | |
| CN112799132B (en) | Micro-local linear noise suppression method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |