CN102081168B  Method for improving surface wave attenuation performance in seismic data processing  Google Patents
Method for improving surface wave attenuation performance in seismic data processing Download PDFInfo
 Publication number
 CN102081168B CN102081168B CN2009102413953A CN200910241395A CN102081168B CN 102081168 B CN102081168 B CN 102081168B CN 2009102413953 A CN2009102413953 A CN 2009102413953A CN 200910241395 A CN200910241395 A CN 200910241395A CN 102081168 B CN102081168 B CN 102081168B
 Authority
 CN
 China
 Prior art keywords
 ground roll
 surface wave
 apparent velocity
 seismic trace
 point
 Prior art date
Links
 230000002238 attenuated Effects 0.000 claims abstract description 8
 230000000875 corresponding Effects 0.000 claims abstract description 4
 238000005070 sampling Methods 0.000 claims description 11
 238000000034 methods Methods 0.000 claims description 7
 238000001914 filtration Methods 0.000 claims description 6
 239000000203 mixtures Substances 0.000 description 1
Abstract
Description
Technical field
The present invention relates to geophysical exploration technology, is to improve the method that geological data is handled signal to noise ratio (S/N ratio) and resolution during a kind of physical prospecting geological data is handled.
Background technology
In the physical prospecting seismic data processing technology; Surface wave attenuation is the important step during geological data is handled; Also be the important measures that improve geological data signal to noise ratio (S/N ratio) and resolution, the surface wave attenuation method of using at present mainly is to adopt the method for highpass filtering that the ground roll that local distributes is decayed.The cardinal principle of this highpass filtering surface wave attenuation method is according to the lower characteristics of earthquake data surface wave frequency; General earthquake data surface wave frequency all concentrates on about 10 hertz; And usable reflection geological data frequency is generally about 8 to 150 hertz; So just can design the Hipass filter more than 10 hertz, the geological data below 10 hertz is carried out attenuation filter according to the difference decay ground roll of ground roll and significant wave frequency distribution; Keep geological data composition more than 10 hertz, the ground roll energy just has been attenuated like this.This damped system is that single seismic trace ground roll is handled; Do not distinguish useful signal and ground roll during decay; In the decay ground roll, will be attenuated inevitably less than 10 hertz useful signal like this, this has limited the raising of signal to noise ratio (S/N ratio) and resolution during geological data is handled.
Summary of the invention
The object of the invention is to provide a kind of surface wave attenuation method, and this method can be avoided the decay to useful signal in surface wave attenuation, thereby realizes the raising of geological data signal to noise ratio (S/N ratio) and resolution.
Process of the present invention comprises following treatment step:
1) earthquakewaveexciting and gather the single track geological data forms the big gun collection record in the geological data;
2) confirm the apparent velocity of ground roll according to single shot record;
Step 2) apparent velocity of described definite ground roll is to choose one section ground roll lineups; Be the hypotenuse of rightangle triangle with the ground roll lineups simultaneously; Geophone offset with last seismic trace place of this section ground roll deducts the propagation distance of the geophone offset at first seismic trace place as ground roll then; The same time with last seismic trace place of this section ground roll deducts the traveltime of the time at first seismic trace place as ground roll, and ground roll passes distance obtains ground roll divided by the time apparent velocity.
3) single big gun seismologic record data are carried out little time window and spatial window division;
It is that time window is got 50 milliseconds that described time window of step 3 and spatial window are divided, and selects 510 seismic trace on the space.
4) carry out the linear energy The dip sweep to having divided spatial window, select the inclination angle of linear energy the most by force;
The The dip sweep of step 4) is that the maximum time sampling point with first seismic trace of spatial window is the point of rotation; Be made as the O point, any one sampling point B on last seismic trace is linked to be straight line OB with the O point in the spatial window; B forms a projection C point at the bottom of time window; This subpoint and O point form straight line OC, and the angle of establishing OB and OC is Q, establish again B to time of C apart from representing be T; Calculating the surface wave speed that ground roll propagates on the OB straight line is: V=OB*COS (Q)/T; Calculate amplitude quadratic sum then through all sampling points of OB this straight line, calculate according to the method described above last seismic trace of spatial window have a few and O point line the amplitude quadratic sum of sample value of process, search out the maximal value of amplitude quadratic sum then.
5) when the size of the corresponding speed V of the peaked straight line of amplitude quadratic sum in step 2) distribution between the maximum apparent velocity of determined ground roll and the minimum apparent velocity, affirmation is a ground roll, need be attenuated, otherwise be useful signal, does not need to decay;
6) according to the result of step 5), carry out surface wave attenuation with the method for highpass filtering, can avoid decay simultaneously, thereby realize the raising of geological data signal to noise ratio (S/N ratio) and resolution useful signal along geological data on the pairing straight line of the maximal value of amplitude quadratic sum.
The present invention is through the linear energy method for scanning; Well separate ground roll with useful signal, then ground roll is decayed, finally realized protection useful signal; Having reached geological data simultaneously improves noise in handling and than the purpose of resolution, has improved the surface wave attenuation performance.
Description of drawings
Fig. 1 is the original earthquake data single shot record;
Fig. 2 the present invention has carried out the single shot record after the surface wave attenuation to Fig. 1.
Embodiment
The method of surface wave attenuation performance during raising geological data of the present invention is handled through openair earthquake data acquisition, forms the big gun collection record in the geological data; According to this big gun collection record, our roughly eye estimate apparent velocity scope is then according to the visual observation ground roll apparent velocity scope that accurate Calculation is concrete again; The general apparent velocity scope that requires to be slightly larger than true ground roll of this velocity range is carried out very little time window to single big gun seismologic record data then and spatial window is divided, if ground roll is more weak; And disperse severe, so time window choose can select a little bit smaller, as 50 milliseconds; Select also corresponding a little bit smallerly on the space, can select 5 seismic traces, if the ground roll energy is more intense; And energy is more concentrated, and time window is chosen and can be selected more so, as 100 milliseconds; Select also greatly correspondingly on the space, can select 9 seismic traces.Spatial window to having divided carries out the linear energy The dip sweep; Select the inclination angle of the strongest linear energy,, then think ground roll if the inclination angle of strong linear energy is just in time in the distribution range between maximum apparent velocity and minimum apparent velocity; So extract the earthquake sampling point data of ground roll lineups; With highpass filtering the sampling point data that extract are carried out surface wave attenuation, so just reached the purpose of decay ground roll, improved the signal to noise ratio (S/N ratio) and the resolution of geological data simultaneously.
The concrete implementation procedure of the present invention comprises following treatment step:
1) gathers the single track geological data, form the big gun collection record in the geological data, like Fig. 1;
2) the roughly apparent velocity scope according to single shot record Fig. 1 calculating and acquisition ground roll is that this velocity range of 800 metre per second (m/s)s2200 metre per second (m/s) is generally bigger than ground roll apparent velocity scope;
Ground roll is constantly decay in communication process; The ground roll apparent velocity is unstable; The estimated apparent velocity of step 2 is the apparent velocity of strong ground roll just, and the present invention need be that the center provides a float value and confirms ground roll maximum apparent velocity and minimum apparent velocity at the surface wave speed that estimates with step 2.
As: step 2 estimates that the surface wave speed that comes out is 1000 a meters/per second; At this moment need provide 300 meters/per second of a float value; The present invention can confirm that the maximum apparent velocity of ground roll is 1000 meters/per second+300 meter/per second=1300 meter/per second so, and minimum apparent velocity is 1000 meters/per second300 meter/per second=700 meter/per second.The distribution range that can conveniently confirm the ground roll apparent velocity like this is 700 meters/per second to 1300 meter/per second.
3) then single big gun seismologic record data are carried out very little time window and spatial window division, Fig. 1 ground roll is strong, not too disperses, and time window can be located 100 milliseconds, can select 9 seismic traces on the space;
4) carrying out the linear energy The dip sweep to having divided spatial window, is the point of rotation with the maximum time sampling point of first seismic trace of spatial window, is made as the O point; Any one sampling point B on last seismic trace in the spatial window; Be linked to be straight line OB with the O point, B forms a projection C point at the bottom of time window, and this subpoint and O point form straight line OC; If the angle of OB and OC is Q; Establishing B again is T to time of the distance expression of C, and can calculate the surface wave speed that ground roll propagates so be on this straight line of OB: V=OB*COS (Q)/T calculates the amplitude quadratic sum through all sampling points of this straight line of OB then; By same method calculate last seismic trace of spatial window have a few and O point line the amplitude quadratic sum of sample value of process, search out the maximal value of amplitude quadratic sum then;
5) if in the distribution range of size between maximum apparent velocity of the determined ground roll of step 3) and minimum apparent velocity of the speed V of the peaked straight line correspondence of this amplitude quadratic sum; Then think ground roll; Need be attenuated, otherwise be useful signal, not need decay;
6) owing to ground roll is separated with useful signal well; At this moment carrying out surface wave attenuation with the method for highpass filtering, so just reached and not only decayed ground roll simultaneously but also protected the purpose of useful signal along geological data on the pairing straight line of the maximal value of amplitude quadratic sum;
Embodiment 1 is the process that the method for surface wave attenuation performance applied to actual earthquake data sheet big gun record during raising geological data of the present invention was handled, and Fig. 1 is that original earthquake data single shot record: Fig. 2 has carried out the single shot record after the surface wave attenuation to Fig. 1.Is 800 meter per seconds2200 meter per seconds to Fig. 1 single shot record through calculating the surface wave speed approximate range, and the ground roll energy is stronger relatively, and more concentrated, so time window has been selected 100 milliseconds, the space seismic trace is selected 9 seismic traces.
Can see that from Fig. 1 and Fig. 2 comparative analysis the ground roll among Fig. 1 is had letter to imitate decay, and the lineups of useful signal are very fuzzy in Fig. 1, but but very clear among the Fig. 2 after ground roll is attenuated.The present invention carries out the method for surface wave attenuation along the ground roll lineups huge superiority is to improve the signal to noise ratio (S/N ratio) and the resolution of geological data, but also can avoid useful signal to be attenuated.
Claims (2)
 One kind improve geological data handle in the method for surface wave attenuation performance, it is characterized in that: comprise following treatment step:1) earthquakewaveexciting and gather the single track geological data forms the big gun collection record in the geological data;2) confirm the apparent velocity of ground roll according to single shot record;The apparent velocity of described definite ground roll is to choose one section ground roll lineups; Be the hypotenuse of rightangle triangle with the ground roll lineups simultaneously; Geophone offset with last seismic trace place of this section ground roll deducts the propagation distance of the geophone offset at first seismic trace place as ground roll then; The same time with last seismic trace place of this section ground roll deducts the traveltime of the time at first seismic trace place as ground roll, and the surface wave propagation distance obtains the apparent velocity of ground roll divided by the time;3) single big gun seismologic record data are carried out little time window and spatial window division;4) readyportioned spatial window is carried out the linear energy The dip sweep, select the inclination angle of linear energy the most by force;The dip sweep is that the maximum time sampling point with first seismic trace of spatial window is the point of rotation; Be made as the O point, any one sampling point B on last seismic trace is linked to be straight line OB with the O point in the spatial window; B forms a projection C point at the bottom of time window; This subpoint and O point form straight line OC, and the angle of establishing OB and OC is Q, establish again B to time of C apart from representing be T; Calculating the surface wave speed that ground roll propagates on the OB straight line is: V=OB*COS (Q)/T; Calculate amplitude quadratic sum then through all sampling points of OB this straight line, calculate according to the method described above last seismic trace of spatial window have a few and O point line the amplitude quadratic sum of sample value of process, search out the maximal value of amplitude quadratic sum then;5) when the size of the corresponding speed V of the peaked straight line of amplitude quadratic sum in step 2) distribution between the maximum apparent velocity of determined ground roll and the minimum apparent velocity, affirmation is a ground roll, need be attenuated, otherwise be useful signal, does not need to decay;6) according to the result of step 5), carry out surface wave attenuation with the method for highpass filtering, can avoid decay simultaneously, thereby realize the raising of geological data signal to noise ratio (S/N ratio) and resolution useful signal along geological data on the pairing straight line of the maximal value of amplitude quadratic sum.
 2. method according to claim 1 is characterized in that it is that time window is got 50 milliseconds that described time window of step 3) and spatial window are divided, and selects 510 seismic trace on the space.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN2009102413953A CN102081168B (en)  20091201  20091201  Method for improving surface wave attenuation performance in seismic data processing 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN2009102413953A CN102081168B (en)  20091201  20091201  Method for improving surface wave attenuation performance in seismic data processing 
Publications (2)
Publication Number  Publication Date 

CN102081168A CN102081168A (en)  20110601 
CN102081168B true CN102081168B (en)  20120718 
Family
ID=44087243
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN2009102413953A CN102081168B (en)  20091201  20091201  Method for improving surface wave attenuation performance in seismic data processing 
Country Status (1)
Country  Link 

CN (1)  CN102081168B (en) 
Families Citing this family (6)
Publication number  Priority date  Publication date  Assignee  Title 

CN103091714B (en) *  20111028  20151118  中国石油化工股份有限公司  A kind of selfadaptation surface wave attenuation method 
CN103217709B (en) *  20120118  20150923  中国石油天然气集团公司  A kind of surface wave attenuation method improving geological data signal to noise ratio (S/N ratio) and resolution 
CN104570115B (en) *  20131016  20170308  中国石油天然气集团公司  A kind of surface wave attenuation method and device 
CN104330826A (en) *  20140731  20150204  中国石油化工股份有限公司  A method for removing various noises under the condition of complex surface 
CN104199104B (en) *  20140815  20170308  中国石油天然气集团公司  3D seismic data surface wave attenuation method and device 
CN105242309B (en) *  20151019  20171110  中国石油天然气集团公司  The drawing method and device of rule interference in a kind of geological data 
Citations (2)
Publication number  Priority date  Publication date  Assignee  Title 

US5521881A (en) *  19940902  19960528  Exxon Production Research Company  Method of processing seismic data having multiple reflection noise 
CN1215424C (en) *  20020326  20050817  周熙襄  Handling method for eliminating strong correlated interference waves in seismic exploration recordings 

2009
 20091201 CN CN2009102413953A patent/CN102081168B/en active IP Right Grant
Patent Citations (2)
Publication number  Priority date  Publication date  Assignee  Title 

US5521881A (en) *  19940902  19960528  Exxon Production Research Company  Method of processing seismic data having multiple reflection noise 
CN1215424C (en) *  20020326  20050817  周熙襄  Handling method for eliminating strong correlated interference waves in seismic exploration recordings 
NonPatent Citations (1)
Title 

高少武等.自适应相干噪声衰减技术.《石油地球物理勘探》.2003,第38卷(第3期),242246. * 
Also Published As
Publication number  Publication date 

CN102081168A (en)  20110601 
Similar Documents
Publication  Publication Date  Title 

Gazdag  Wave equation migration with the phaseshift method  
US7941273B2 (en)  System and method for using timedistance characteristics in acquisition, processing, and imaging of TCSEM data  
Cassidy et al.  Ground penetrating radar data processing, modelling and analysis  
CN101910870B (en)  Seismic sensor devices  
US7821869B2 (en)  Methods of enhancing separation of primary reflection signals and noise in seismic data using radon transformations  
USRE41656E1 (en)  Method and system for reducing effects of sea surface ghost contamination in seismic data  
Jokat  Seismic investigations along the western sector of Alpha Ridge, Central Arctic Ocean  
US7239578B2 (en)  Removal of noise from seismic data using radon transformations  
CA2600657C (en)  Radon transformations for removal of noise from seismic data  
EP2189818B1 (en)  Method of wavefield extrapolation for singlestation, dualsensor towed streamer signals  
ChávezGarcía et al.  Lateral propagation effects observed at Parkway, New Zealand. A case history to compare 1D versus 2D site effects  
Barr et al.  Attenuation of watercolumn reverberations using pressure and velocity detectors in a waterbottom cable  
US4353121A (en)  High resolution, marine seismic stratigraphic system  
CN104570125B (en)  A kind of method utilizing well data to improve image taking speed model accuracy  
US7366054B1 (en)  TauP filters for removal of noise from seismic data  
CN102305941B (en)  Method for determining stratum stack quality factor by direct scanning of prestack time migration  
CN104635269B (en)  A kind of method that fractured igneousrock type reservoir is predicted based on prestack orientation angle gathers  
AU2013245550B2 (en)  Systems and methods for highresolution imaging using separated wavefields  
EP2778718B1 (en)  Systems and methods for frequencydomain filtering and spacetime domain discrimination of seismic data  
CN1952691A (en)  Seismic streamer system and method  
Quinn et al.  Optimal processing of marine highresolution seismic reflection (Chirp) data  
CN101014881A (en)  System for attenuation of water bottom multiples in seismic data recorded by pressure sensors and particle motion sensors  
EA026344B1 (en)  System and method for acquisition and processing of elastic wavefield seismic data  
CN101598809A (en)  A kind of selfadaptation is eliminated the method for linear programming noise and multiple reflection interference  
CN101334483B (en)  Method for attenuating rayleigh wave scattered noise in earthquake datahandling 
Legal Events
Date  Code  Title  Description 

PB01  Publication  
C06  Publication  
SE01  Entry into force of request for substantive examination  
C10  Entry into substantive examination  
GR01  Patent grant  
C14  Grant of patent or utility model 