NO812801L

NO812801L - F-K GEOPHYSICAL OPERATIONS INCLUDING SEISMIC RECORDING FILTERING

Info

Publication number: NO812801L
Application number: NO812801A
Authority: NO
Inventors: William Harold Ruehle
Original assignee: Mobil Oil Corp
Priority date: 1980-08-20
Filing date: 1981-08-19
Publication date: 1982-02-22
Also published as: GB2082771A

Abstract

Seismiske signaler blir transformert til en f-k rekke med digitale sampler som viser amplituden som en funksjon av frekvensen og bølgetallet. f-k rekken blir så filtrert ved vekttallgivning til samplene i lukkede frekvensområder med vekttallgivningen avhengig av signal og støy. Seismisk opptegning som har en øket frem-stilling av jordens formasjon blir så frembragt fra den filtrerte rekken.Seismic signals are transformed into an f-k series of digital samples showing the amplitude as a function of frequency and wave count. The f-k series is then filtered by weighting the samples in closed frequency ranges with the weighting depending on signal and noise. Seismic plots that have an increased representation of the earth's formation are then generated from the filtered array.

Description

Foreliggende oppfinnelse angår geofysisk undersøkelse ogThe present invention relates to geophysical investigation and

mer bestemt geofysiske operasjoner utført på f-k transformering av seismiske opptegninger. more specifically geophysical operations performed on f-k transformation of seismic records.

Ved seismiske undersøkelser blir seismogrammer tilveiebragtDuring seismic surveys, seismograms are provided

i felten behandlet for å gjøre de lettere å tyde. Feltseismogrammer innbefatter tilfeldig, eller usammenhengende støy. som er bevirket av de benyttede instrumentene eller de geofysiske formasjonenes spesielle uklare egenskaper, og sammenhengende støy som innbefatter ekko-bilder, flerrefleks joner, og overflatebølger. Geofysiske operasjoner som blir utført på seismogrammer for å gjøre de mer synlige innbefatter konvolvering, dekonvulvering, signaturvurdering, designering, reduksjon av flerrefleks joner, og korrelasjon. Noen geofysiske operasjoner som reduserer støy kan bli utført i felten. F.eks. kan geofoner bli anordnet i et spesielt mønster for å dekonvolvere, eller filtrere en bestemt støy-type. Den mer fremherskende praksisen er imidlertid å ut-føre disse geofysiske operasjoner på seismogrammer etter at de er blitt opptegnet og transportert til et sentralt behandlingssenter. in the field treated to make them easier to interpret. Field seismograms include random, or incoherent, noise. which is caused by the instruments used or the special unclear properties of the geophysical formations, and coherent noise which includes echo images, multiple reflections and surface waves. Geophysical operations performed on seismograms to make them more visible include convolution, deconvulsion, signature estimation, designating, reduction of multiple reflections, and correlation. Some geophysical operations that reduce noise can be carried out in the field. E.g. geophones can be arranged in a special pattern to deconvolve, or filter a particular noise type. The more prevalent practice, however, is to perform these geophysical operations on seismograms after they have been recorded and transported to a central processing centre.

I løpet av tidligere utvikling av seismiske undersøkelserDuring the previous development of seismic surveys

ble seismogrammene opptegnet som analoge spenninger på magnetbånd. Seismogrammer blir fremdeles vanligvis opptegnet på magnetbånd, men de blir opptegnet som digitale sampler som fremstiller amplituden til et mottatt seismisk signal.som the seismograms were recorded as analogue voltages on magnetic tape. Seismograms are still usually recorded on magnetic tape, but they are recorded as digital samples that represent the amplitude of a received seismic signal.

en funksjon av tiden. Siden seismogrammene vanligvis blir tilveiebragt langs en undersøkelselinje på jordoverflaten kan digitale sampler bli formet til x-t rekker med hver sampel i rekken representererende amplituden til det seismiske signalet som en funksjon av den horisontale avstanden av tiden.. Når slike rekker blir visuelt reprodusert ved plotting eller lignende blir seismiske seksjoner frembragt. En seismisk seksjon viser underjordiske lag til en a function of time. Since the seismograms are usually provided along a survey line on the Earth's surface, digital samples can be formed into x-t arrays with each sample in the array representing the amplitude of the seismic signal as a function of horizontal distance in time. When such arrays are visually reproduced by plotting or the like seismic sections are generated. A seismic section shows underground layers to a

jordseksjon. Det er det prinsippielle verktøyet som geo-fysikerne studerer for å bestemme egenskapen til jordens underjordiske formasjoner. Før ;.en. rekke med seismiske sampler kan bli omformet til en seismisk seksjon som kan bli tydet av geofysikeren må rekken bli utstrakt behandlet av en eller flere geofysiske, behandlingsoperasjoner nevnt tidligere. En av disse geofysiske behandlingsoperasjonene er filtrering (konvolvering), som er den mest vanlige. soil section. It is the principle tool that geo-physicists study to determine the nature of the Earth's underground formations. For one. array of seismic samples can be transformed into a seismic section that can be interpreted by the geophysicist, the array must be extensively processed by one or more geophysical processing operations mentioned earlier. One of these geophysical processing operations is filtering (convolution), which is the most common.

Feltarbeiderne har omformet x-t rekkene til rekker som representerer amplituden som en funksjon av frekvensen The field workers have transformed the x-t series into series representing the amplitude as a function of frequency

og bølgenummeret. Dette er .vanligvis omtalt som en f-k transformasjon. Kutrona, et al., beskriver i "Optical and the wave number. This is usually referred to as an f-k transformation. Kutrona, et al., describe in "Optical

Data Processing" I.E.E.E. Proce.edings (1959) f-k filtreringData Processing" I.E.E.E. Proce.edings (1959) f-k filtering

i optisk behandling, av data opptegnet som variable tetthets fotografier.. Ved denne behandlingstypen blir plant mono-kromatisk lys ført gjennom en fotografisk film med variabel tetthet, så gjennom en kuleformet linse for å frembringe en todimensjonal Fourier transformasjon av gjenstanden på filmen. Jackson, P.L., beskriver i "Analysis of Variable Density Seismograms By Optical Data Processor", i Geophysics, No. 1, (1965) filtrering av f-k seismiske seksjoner ved optisk behandling. US patent 3 370 268 og Dobrin, M.V., in optical processing, of data recorded as variable density photographs.. In this type of processing, plane monochromatic light is passed through a photographic film of variable density, then through a spherical lens to produce a two-dimensional Fourier transform of the object on the film. Jackson, P.L., describes in "Analysis of Variable Density Seismograms By Optical Data Processor", in Geophysics, No. 1, (1965) filtering of f-k seismic sections by optical processing. US Patent 3,370,268 and Dobrin, M.V.,

i artikkelen "Short Note on Optical Processing and Filtering" i Geophysics, no. 6 (1965) beskriver også optisk filtrering av f-k seismiske seksjoner. Artikkelen av Dobrin beskriver begrensningene ved optisk behandling. Pga. disse begrensningene blir største delen av geofysisk behandling utført på digial data behandlere, og ikke optiske behandlere. in the article "Short Note on Optical Processing and Filtering" in Geophysics, no. 6 (1965) also describes optical filtering of f-k seismic sections. The article by Dobrin describes the limitations of optical processing. Because of. these limitations mean that the largest part of geophysical processing is carried out on digital data processors, and not optical processors.

Ved digial behandling har f-k transformasjonen blitt benyttes som et instrument for å studere seismisk filtrering. US patentene nr. 3 284 763 og 3 274 541 beskriver hastighets-filtrering som blir utført i tidsdomener på digitale seismogrammer. Disse patentene beskriver en type hastighets-filtrering som vanligvis er henvist til som "kakestykke" . pga. det filtrerte områdets form i f-k transformasjonen til den seismiske seksjon. De to ovenfornevnte US patentene utfører ikke en filtreringsoperasjon på f-k seksjonen, men all filtreringen er i tidsdomenen på seksjonene i de normale x-t rekkene. US patent nr. 3 550 073 er en forbedring av "kakestykket" filtreringen. In digital processing, the f-k transformation has been used as an instrument to study seismic filtering. US Patent Nos. 3,284,763 and 3,274,541 describe velocity filtering that is performed in time domains on digital seismograms. These patents describe a type of rate filtering commonly referred to as "piece of cake". because of. the shape of the filtered area in the f-k transformation of the seismic section. The two above-mentioned US patents do not perform a filtering operation on the f-k section, but all the filtering is in the time domain on the sections in the normal x-t series. US Patent No. 3,550,073 is an improvement on the "piece of cake" filtration.

Britisk patent nr. 1 099 483 beskriver en geofysisk behandlings-operasjon noen ganger henvist til som en horisontal stabling. f-k transformasjonen er nyttig for å forestå horisontalstablingen, imiderltid er alle operasjonene i horisontalstablingen utført i tidsdomenen, ikke på f-k transformerte seksjoner. British Patent No. 1,099,483 describes a geophysical processing operation sometimes referred to as a horizontal stack. The f-k transformation is useful for performing the horizontal stacking, however, all operations in the horizontal stacking are performed in the time domain, not on f-k transformed sections.

Alle ovenfornevnte typer med romfilter har ulempen av at deAll the above-mentioned types with room filters have the disadvantage that they

har svak steilhet mellom passeringsbåndet og sprrebåndet.has a slight steepness between the passing band and the scattering band.

Disse filtrene mangler således selektivtitet ved diskrimi-neringen mellom ønskelige og uønskelige tilfeller som er nær hverandre i f-k transformasjonen. These filters thus lack selectivity when discriminating between desirable and undesirable cases that are close to each other in the f-k transformation.

Det er et formål med foreliggende oppfinnelse å tilveiebringe et filter som diskriminerer mellom tilfeller med hastigheter eller fall som er svært tett sammen i f-k transformasjon. It is an object of the present invention to provide a filter that discriminates between cases with speeds or drops that are very close together in f-k transformation.

Grunnen for at geofysisk behandlingsoperasjoner ikke tidligere er blitt utført på f-k transformerte seismiske seksjoner er at man har trodd at. en slik operasjon ville bryte et funda-mentalt sampel teorem. Man har imidlertid funnet at visse geofysiske behandlingsoperasjoner, spesielt filtrering, kan på nyttig måte bli utført på f-k transformasjoner til seismiske seksjoner. The reason why geophysical processing operations have not previously been carried out on f-k transformed seismic sections is that it has been believed that. such an operation would break a fundamental sample theorem. However, it has been found that certain geophysical processing operations, particularly filtering, can be usefully performed on f-k transforms to seismic sections.

Det er følgelig et formål med foreliggende oppfinnelse-It is therefore an object of the present invention-

å utføre geofysikalske behandlingsoperasjoner på seis-to carry out geophysical treatment operations on seismic

miske seksjoner (rekker), som representerer amplitude og fase som en funksjon av frekvens og bølgetall. misc sections (rows), which represent amplitude and phase as a function of frequency and wavenumber.

Ifølge foreliggende oppfinnelse representerer seismiske signaler amplituden til seismiske refleksjoner som en funksjon av tiden og avstanden' blir transformert til en seksjon som representerer amplituden som en funksjon av frekvensen og bølgetallet. Geofysiske behandlingeoperasjon-er blir utført på disse f-k seksjonene. I samsvar med oppfinnelsen blir f-k seksjonen filtrert ved vekttallgivelse til alle samplene i foreskrevede tette områder av frekvens. og- bølgetallet. According to the present invention, seismic signals representing the amplitude of seismic reflections as a function of time and distance are transformed into a section representing the amplitude as a function of frequency and wave number. Geophysical treatment operations are carried out on these f-k sections. In accordance with the invention, the f-k section is filtered by weighting all the samples in prescribed tight ranges of frequency. and- the wave number.

I samsvar med et ytterligere•trekk ved oppfinnelsen blir alle samplene vraket med unntak av de som representerer refleksjonene fra geofysiske former som har en foreskrevet form. Ved riktig definering av grensene til områdene av de vrakede samplene i f-k transformasjonen kan refleksjonene fra formasjonene med en bestemt form bli filtrert eller . sluppet igjennom.. Når f-k transformasjonen er konvertert tilbake til en normal seismisk seksjon vil de ønskede formasjonene bli fremvist. Det er f.eks. mulig å filtrere ref-, leksjoner fra bratte hellende formasjoner, eller fra formasjoner som har intet fall. In accordance with a further feature of the invention, all the samples are scrapped with the exception of those representing the reflections from geophysical forms having a prescribed shape. By correctly defining the boundaries of the areas of the wrecked samples in the f-k transformation, the reflections from the formations with a certain shape can be filtered or . let through.. When the f-k transformation is converted back to a normal seismic section the desired formations will be displayed. It is e.g. possible to filter reflections from steeply sloping formations, or from formations that have no dip.

I samsvar med foreliggende oppfinnelse har dessuten f-k filteret svært steile sider s.lik at fallselektiviteten er forbedret i forhold til den til optiske behandlere og hastighetsfiltere til tidligere kjente innretninger. In accordance with the present invention, the f-k filter also has very steep sides such that the drop selectivity is improved in relation to that of optical processors and velocity filters of previously known devices.

Ovenfornevnte og andre formål, trekk og fordeler ved foreliggende oppfinnelse skal nærmere beskrives med henvisning til de medfølgende tegninene, hvor:. The above and other purposes, features and advantages of the present invention shall be described in more detail with reference to the accompanying drawings, where:

Fig. 1 viser en seismisk seksjon i- dens normale amplitude-form som en funksjon av avstand og tid,' og f-k transformasjonen til denne seksjonen, Fig. 1 shows a seismic section i- its normal amplitude form as a function of distance and time,' and the f-k transformation of this section,

fig. 2 viser en seismisk seksjon og dens f-k transformasjon, fig. 3 viser et f-k filter av foreliggende oppfinnelse, fig. 2 shows a seismic section and its f-k transformation, fig. 3 shows an f-k filter of the present invention,

fig. 4 viser en seismisk seksjon, filtrering av f-k transformasjonen ocj den seismiske seksjonen etter at den er blitt filtrert, fig. 4 shows a seismic section, filtering the f-k transform and the seismic section after it has been filtered,

fig. 5 viser et flytediagram over filtreringen ifølge oppfinnelsen, fig. 5 shows a flow diagram of the filtration according to the invention,

fig. 6 viser impulsreaksjonen. til et tidsdomenefilter for fig. 6 shows the impulse response. to a time domain filter for

fjerning av vannbundede reflekterte brytninger,removal of water-bound reflected refractions,

fig. 7 viser en seismisk seksjonsmodell,fig. 7 shows a seismic section model,

fig. 8 viser f-k transformasjoner av seksjonen på fig. 7, fig. 9 viser et f-k filter som vraker alle sampler unntatt de som representerer refleksjonen av grensesnittet fig. 8 shows f-k transformations of the section in fig. 7, fig. 9 shows an f-k filter which discards all samples except those representing the reflection of the interface

som har fall på 0,032 sek. pr. spor,which has a drop of 0.032 sec. per traces,

fig.10 viser f-k transformasjonen til fig.. 7, filtrert med fig.10 shows the f-k transformation of fig.. 7, filtered with

filteret på. fig. 9,filter on. fig. 9,

fig.11 viser den seismiske seksjonen som følge av at fig. 10 fig.11 shows the seismic section as a result of fig. 10

er omformet til en x-t rekke,is transformed into an x-t array,

fig.12 viser et f-k filter for å vrake alle unntatt 0-falls fig.12 shows an f-k filter to discard all but 0 cases

refleksjonene,the reflections,

fig.13 viser et f-k filter for å vrake 0-fall refleksjonen, fig.14 viser en seismisk seksjonsmodell, fig.13 shows an f-k filter to wreck the 0-fall reflection, fig.14 shows a seismic section model,

fig.15 viser f-k transformasjoner av seksjonen på fig. 14, fig.16 viser f-k seksjonen på fig. 15 filtrert med filteret fig.15 shows f-k transformations of the section in fig. 14, fig. 16 shows the f-k section of fig. 15 filtered with the filter

på fig. 12.,on fig. 12th,

fig.17 viser x-t transformasjonen på fig. 16,fig.17 shows the x-t transformation in fig. 16,

fig.18 viser f-k transformasjonen på fig. 15, filtrertfig.18 shows the f-k transformation in fig. 15, filtered

med filteret på fig. 13,with the filter in fig. 13,

fig.19 viser den seismiske seksjon etter filtreringen vist fig.19 shows the seismic section after the filtering shown

på fig.. 18, fig.20 viser en seismisk seksjonsmodell, on fig. 18, fig. 20 shows a seismic section model,

fig.21 viser den seismiske seksjonen på fig. 20 etter at den er blitt filtrert i samsvar med foreliggende fig.21 shows the seismic section of fig. 20 after it has been filtered in accordance with the present

oppfinnelse,invention,

fig.22 viser den seismiske' seksjonen på fig. 20 etter' at den er blitt filtrert med et tidligere kjent kake-stykkefilter, fig.22 shows the seismic section of fig. 20 after it has been filtered with a previously known pie-piece filter,

fig. 23 viser f-k transformasjonen på fig. 20,fig. 23 shows the f-k transformation in fig. 20,

fig. 24bg 25 er feltseismogrammer,fig. 24bg 25 are field seismograms,

fig. 26 viser et f-k filter for signal/støyøkning,fig. 26 shows an f-k filter for signal/noise gain,

fig. 26a viser f-k transformasjonen på fig. 24,fig. 26a shows the f-k transformation in fig. 24,

fig. 27 viser feltopptegningen til fig. 24 etter filtrering med foreliggende oppfinnelse,- fig. 27 shows the field diagram of fig. 24 after filtering with the present invention,-

fig. 28 viser feltopptegningen til fig. 25. etter at denfig. 28 shows the field drawing of fig. 25. after that the

er blitt filtrert i samsvar méd foreliggende oppfinnelse, og has been filtered in accordance with the present invention, and

fig. 29 viser et eksempel på en databehandlingsinnretningfig. 29 shows an example of a data processing device

på hvilken oppfinnelsen kan bli utført.on which the invention may be carried out.

Før beskrivelsen av foreliggende oppfinnelse er det nyttigBefore the description of the present invention, it is useful

å betrakte den seismiske seksjonen på fig. 1 og dens f-k transformasjon. En seismisk seksjon viser amplituden til seismiske signaler som en funksjon av avstanden x langs en undersøkelses linje og tiden t etter skuddet som frembringer refleksjonene. En seismisk seksjon viser ideelt underjordiske lag som frembringer de seismiske refleksjonene, idet de idealiserte grenseflatene 11, 12 og 13 er vist på fig. 1. En seismisk seksjon eller opptegning' er vist med rekker med digitale sampler.Digitalsamplene i kolonnen 46 utgjør, amplituden for de seismiske refleksjonene som en refleksjon av tiden ved stedet 46 langs undersøkelses-linjen. Amplituden er betegnet A(x,t). For å være nøyakt-ig :er selve digitalsamplene henvist til som en rekke, mens derimot utgangsplotten til disse samplene er en seismisk seksjon. Det er imidlertid vanlig innenfor geofysikken å be-nytte ordet "s.eksjon" for å beskrive både plotten og rekken med digitale sampler og denne praksisen vil også to consider the seismic section of fig. 1 and its f-k transformation. A seismic section shows the amplitude of seismic signals as a function of the distance x along a survey line and the time t after the shot that produces the reflections. A seismic section shows idealized subsurface layers which produce the seismic reflections, the idealized interfaces 11, 12 and 13 being shown in fig. 1. A seismic section or plot' is shown with rows of digital samples. The digital samples in column 46 constitute the amplitude of the seismic reflections as a reflection of time at location 46 along the survey line. The amplitude is denoted A(x,t). To be precise: the digital samples themselves are referred to as a series, while the output plot of these samples is a seismic section. However, it is common in geophysics to use the word "section" to describe both the plot and the series of digital samples, and this practice will also

bli fulgt her. Med henvisning til en geofysisk behandlings-operasjon på enseismisk seksjon skal det bemerkes at det her er tale om digitale sampler som skal bli behandlet. be followed here. With reference to a geophysical processing operation on an seismic section, it should be noted that these are digital samples to be processed.

Fig. 1 viser også en f-k transformasjon av en vanlig seismisk seksjon. Fourier transformasjoner som omformer rekken med amplituder som en funksjon av tiden og avstanden til rekken av amplituden som en funksjon av frekvensen og bølgetallet er velkjent. Den høyre delen av fig. 1 viser rekken av digitale sampler for hver frekvens og bølgetall-verdi i rekken. De digitale samplene ér betegnet A(f,k). Fig. 1 also shows an f-k transformation of a normal seismic section. Fourier transforms that transform the series of amplitudes as a function of time and distance into the series of amplitudes as a function of frequency and wavenumber are well known. The right part of fig. 1 shows the series of digital samples for each frequency and wavenumber value in the series. The digital samples are denoted A(f,k).

Fig. 2 viser noen av egenskapene til f-k transformasjonene som er nyttige ved séismisk behandling. Linjen 14 på den vanlige seismiske seksjonen transformeres til linjen 15 Fig. 2 shows some of the properties of the f-k transformations that are useful in seismic processing. Line 14 on the normal seismic section is transformed into line 15

på f.-k tranformasjonen. Linjene 16, som er på forskjellige steder,, men har samme vinkelen eller fall, transformeres til en enkel, linje 17. Linjen 18 transformeres til linjen 19. Sinusbølgene 20 transformeres til punktet 21. on the f.-k transformation. The lines 16, which are in different places, but have the same angle or dip, are transformed into a single line 17. The line 18 is transformed into the line 19. The sine waves 20 are transformed into the point 21.

Man har funnet at filtreringen av f-k transformeringene til seismiske seksjoner kan bli tilveiebragt ved å vrake alle samplene innenfor et foreskrevet lukketområde med frekvens og bølgetall. Et slikt filter er vist på fig. 3.. Fig. 3 viser en rekke med digital sampler som representerer amplituden som en funksjon av frekvensen f og bølgetallet k. It has been found that the filtering of the f-k transformations of seismic sections can be provided by scraping all the samples within a prescribed closed range of frequency and wavenumber. Such a filter is shown in fig. 3.. Fig. 3 shows a series of digital samples representing the amplitude as a function of the frequency f and the wavenumber k.

En vertikal dimensjon har blitt tillagt til standard, f-k plottet på fig. 3 for å vise amplituden. Amplituden W(i,j) A vertical dimension has been added to the standard, f-k plotted in fig. 3 to show the amplitude. The amplitude W(i,j)

i dette eksempel er enten 1 eller 0 som er vanlig. I samsvar med foreliggende oppfinnelse blir alle samplene innenfor in this example, either 1 or 0 is common. In accordance with the present invention, all the samples are within

et foreskrevet lukket område med frekvens og bølgetall som har grensene 22-29 vraket. Ved den enkleste utførelsen er denne vrakningen på enkel måte tilveiebragt ved å multiplisere alle samplene i dette lukkede område med 0. På fig. 3 er dette vist ved at alle samplene i det lukkede området som har grensene 22-29 er 0, mens alle andre sampler er 1. a prescribed closed area of frequency and wavenumber which has the limits 22-29 wreck. In the simplest embodiment, this scraping is simply provided by multiplying all the samples in this closed area by 0. In fig. 3, this is shown by the fact that all the samples in the closed area that have the boundaries 22-29 are 0, while all other samples are 1.

Nytten av et slikt filter ved seismisk behandling er vistThe usefulness of such a filter in seismic processing is shown

på fig. 4, hvor en idealisert seismisk seksjon 30. er vist.. Seksjonen viser horisontale reflekterende grenseflater 31 on fig. 4, where an idealized seismic section 30. is shown. The section shows horizontal reflective boundaries 31

og 32 og fallgrenseflåtene 33, 34 og 35 som alle har samme fall. Den seismiske seksjonen viser også vannbundne reflekterte brytninger, som fremvist ved hjelp av linjene 36-39. and 32 and the drop limit rafts 33, 34 and 35 which all have the same drop. The seismic section also shows water bound reflected refractions, as shown by lines 36-39.

Dette er støy som skulle bli fjernet av den seismiske seksjon ved behandlingen. This is noise that should be removed by the seismic section during the processing.

Ved f-k transformasjonen 40 til den seismiske seksjonen blir de horisontale linjene 31 og 32 transformert til enkle vertikale linjer 41 og fallgrenseflåtene .33-35 er blitt transformert til linjen 42. De vannbundne reflekterte brytningene er blitt transformert til linjene 43-46. I samsvar med foreliggende oppfinnelse blir disse brytningene filtrert ut av den seismiske seksjonen ved å vrake alle samplene i det lukkede område med frekvens og bølgetall antydet med 4?. Pga. samplingsintervallet er det nødvendig at området 47 omhyller et større område enn i hvilket linjene 43-46 ligger. Filtreringen blir tilveiebragt ved å multiplisere alle samplene i dette området med 0; Kun linjene 48 og 49 forblir ved denne filtrerte f-k transformasjonen av den seismiske seksjonen. Når denne blir konvertert tilbake til en vanlig seismisk seksjon forblir grenseflatene 31-35, In the f-k transformation 40 of the seismic section, the horizontal lines 31 and 32 are transformed into simple vertical lines 41 and the dip boundary rafts .33-35 have been transformed into line 42. The water bound reflected refractions have been transformed into lines 43-46. In accordance with the present invention, these refractions are filtered out of the seismic section by wrecking all the samples in the closed area with frequency and wavenumber indicated by 4?. Because of. sampling interval, it is necessary that the area 47 encloses a larger area than in which the lines 43-46 lie. The filtering is provided by multiplying all the samples in this range by 0; Only lines 48 and 49 remain with this filtered f-k transformation of the seismic section. When this is converted back to a normal seismic section, the boundary surfaces 31-35 remain,

men de vannbundne reflekterte brytningene er blitt filtrert ut av seksjonen. Nytten av foreliggende oppfinnelse kan bli demonstrert ved å betrakte vanskeligheten med å fjerne disse vannbundne reflekterte brytningene ved vanlig tid-domenefiltrering. Fig. 6 viser•impulsreaksjonen til et tidsdomenefilter som ville fjerne vannbundne reflekterte brytninger fra en seismisk seksjon. Fagmannen på området vil se at et slikt filter er nesten umulig å realisere på et tidsromfilter. Det er ytterst komplisert og det vil ta så lang tid å utføre det at dets utførelse er i virkelig-heten umulig. but the water bound reflected refractions have been filtered out of the section. The utility of the present invention can be demonstrated by considering the difficulty in removing these water-bound reflected refractions by conventional time-domain filtering. Fig. 6 shows•the impulse response of a time-domain filter that would remove water-bound reflected refractions from a seismic section. The expert in the field will see that such a filter is almost impossible to realize on a time domain filter. It is extremely complicated and it will take so long to carry it out that its execution is in reality impossible.

Som et annet eksempel skal betraktes f-k filtrering avAnother example should be considered f-k filtering of

den idealiserte seismiske seksjonen på fig. 7. På fig. 7the idealized seismic section of fig. 7. In fig. 7

er refleksjonene antydet ved syv forskjellige opptegnings-tider betegnet 50-56 og det er fem refleksjoner ved hver linje. Dvs. refleksjonene 50a-50e har en TQ på 0,1 sek., refleksjonene 51a-51e har en Tq på 0,2 sekunder osv. the reflections are indicated at seven different recording times denoted 50-56 and there are five reflections at each line. That is reflections 50a-50e have a TQ of 0.1 sec., reflections 51a-51e have a Tq of 0.2 seconds, etc.

v v

Ved hvert tidspunkt' er det refleksjoner med fem forskjellige fall mellom 6t = 0 og 6t = 0,032 sek/spor. Disse trettifem refleksjonene transformeres til fem linjer i f-k transformasjonen på fig. 8. Refleksjonene 50å, 51a, 52a... 56a på fiy. 7 transformeres til de vertikale linjene 57 på fig. 8, refleksjonen 50d. 56b transformeres til linjen 58 på fig. 8, At each time' there are reflections with five different dips between 6t = 0 and 6t = 0.032 sec/trace. These thirty-five reflections are transformed into five lines in the f-k transformation in fig. 8. The reflections 50å, 51a, 52a... 56a on fiy. 7 is transformed into the vertical lines 57 in fig. 8, the reflection 50d. 56b is transformed into line 58 in fig. 8,

og refleksjonen 50c... 56c transformeres til linjen 59. Refleksjonene 50d...56d transformeres til linjeseymentene 60a py 60b på fig. 8. Det skal bemerkes at disse segmentene er en enkel linje som vikles rundt k grensen som et resultat av uriktiy rommessig sampling av den seismiske seksjon. Dette er et vanliy tilfelle. De bratteste refleksjonene 50e... and the reflection 50c...56c is transformed into the line 59. The reflections 50d...56d are transformed into the line elements 60a py 60b in fig. 8. It should be noted that these segments are a simple line that wraps around the k boundary as a result of uriktiy spatial sampling of the seismic section. This is a vanliy case. The steepest reflections 50e...

56e på fig. 7 transformeres likeledes inn i en linje som har segmentene 61a, 61b og 61c på fig. 8. 56e in fig. 7 is likewise transformed into a line having segments 61a, 61b and 61c in fig. 8.

Fig. 9 viser et f-k filter som vraker alle refleksjonene unntatt de som nar de steileste fallene på 0,0 32 sek. pr. Fig. 9 shows an f-k filter which discards all the reflections except for those with the steepest drops of 0.0 32 sec. per

spor, dvs. refleksjonene 50e...56e på fig. 7. Det lukkede område i hvilket samplene ikke er vraket har grensene 6a, traces, i.e. the reflections 50e...56e in fig. 7. The closed area in which the samples are not wrecked has the boundaries 6a,

63a, 63b, 63c, 64a, 64b, 64a og 65 som er lineære funksjoner av frekvens og bølgetall. Ved.vraking av alle samplene på utsiden av dette området blir alle refleksjonene som.ikke har et fall på 0,032 sek. pr. spor filtrert fra den seismiske seksjonen. Fig. 10 viser f-k transformasjonen av den seismiske seksjon som er blitt filtrert ved vraking av alle samplene på utsiden av området definert av lineæryrensene 62-65 (fig. 9). Kune linjen med segmentene 61a, 61b oy 61c forblir på fig. 10.. dår denne blir transformert tilbake til en vanlig x-t rekke, fåes den seismiske seksjonen på fig. 11. Kun de steileste refleksjonene 50e....56e forblir etter 63a, 63b, 63c, 64a, 64b, 64a and 65 which are linear functions of frequency and wave number. By scrapping all the samples on the outside of this area, all the reflections that do not have a drop of 0.032 sec. per traces filtered from the seismic section. Fig. 10 shows the f-k transformation of the seismic section which has been filtered by scraping all the samples on the outside of the area defined by the linear filters 62-65 (fig. 9). Kune line with segments 61a, 61b and 61c remains in fig. 10.. when this is transformed back into a regular x-t series, the seismic section in fig. 11. Only the steepest reflections 50e....56e remain behind

f-k filtreringen. Skarpheten er mye bedre enn den som er tilveiebragt med hastighet eller kakestykkfilter. På fig. 11 er den sinusformede støyen mellom refleksjonene et resultat av øket intensitet i f-k rommet langs det ønskede fallet frembragt ved interferens med vrakede tilfeller. f-k filtering. The sharpness is much better than that provided by speed or pie filter. In fig. 11, the sinusoidal noise between the reflections is a result of increased intensity in the f-k space along the desired drop produced by interference with wrecked cases.

Fig. 12 viser et f-k filter som vil vrake alle unntatt 0-fall . tilfellene, dvs. den vil vrake alle refleksjonene untatt refleksjonene 56a...56e på fig. 7. Alle samplene utenfor området spm nar de lineære grensene 66-69 blir vraket. Grensene 67 Fig. 12 shows an f-k filter that will reject all but 0 cases. the cases, i.e. it will destroy all the reflections except the reflections 56a...56e in fig. 7. All the samples outside the range spm when the linear boundaries 66-69 are destroyed. The borders 67

og 69 er konstante bølgetallgrenser. Grensen 67 er et bølgetall som er en liten endedel av det høyeste, bølgetallet +k. Det høyeste bølgetallet står i forhold til det horisontale sample-intervallet som er den norisontale avstanden langs undersøk-elseslinjen mellom stillingene ved hvilken CDP innstilles eller sporene på en feltopptegning eller med en CDP samler er opp-nådd. Dvs. det høyeste bølgenummeret er k = l/(2Ax), hvor Ax er det norisontale sampelihtervallet mellom CDP innstill-ingene. Som vist på fig. 12 er den konstante bølgetallgrens-en 67 tilnærmet 0,2 av det høyeste bølgetallet. Det er klart at denne grensen kan være en hver liten endedel av det høy-este bølgetallet. Grensen 67 er + a/2Ax og grensen 69 er - a/2Ax, hvor a er den lille endedelen. Grensene 66 og 68 er konstante frekvensgrenser, spesielt den laveste og høyeste and 69 are constant wavenumber limits. The limit 67 is a wave number which is a small end part of the highest, the wave number +k. The highest wavenumber is in relation to the horizontal sample interval which is the horizontal distance along the survey line between the positions at which the CDP is set or the tracks on a field plot or with a CDP collector are reached. That is the highest wavenumber is k = l/(2Ax), where Ax is the norizontal sample interval between the CDP settings. As shown in fig. 12, the constant wavenumber limit 67 is approximately 0.2 of the highest wavenumber. It is clear that this limit can be every small end part of the highest wave number. The limit 67 is + a/2Ax and the limit 69 is - a/2Ax, where a is the small end part. Limits 66 and 68 are constant frequency limits, especially the lowest and highest

tilstedeværende frekvensen. De laveste og høyeste frekvens-ene står i forhold til tidssamplingsintervallet. Dvs., f = l/(2A'r). present frequency. The lowest and highest frequencies are in relation to the time sampling interval. That is, f = l/(2A's).

Fig. 14-17 viser et eksempel på driften av filteret påFig. 14-17 shows an example of the operation of the filter on

fig. 12. Fig. 14 viser en seismisk seksjon med tre reflek-sjonsfamilier: 74, 75 og 76. Hver av de har tre par med refleksjoner med forskjellige fall. f-k transformasjonen på fig. 14 er vist på fig. 15. Dette viser resultatet av rommessig aliasering på transformasjonen av steilest fall, fig. 12. Fig. 14 shows a seismic section with three reflection families: 74, 75 and 76. Each of them has three pairs of reflections with different dips. the f-k transformation in fig. 14 is shown in fig. 15. This shows the result of spatial aliasing on the steepest-fall transformation,

en folding av transformasjonen og Nyquist bølgetallet. Når dette er filtrert ved vraking av alle samplene unntatt de som representerer flate, 0-fall, refleksjonene, tilveiebringes f-k transformasjonen på fig. 16. Det skal bemerkes at denne er analog for å behandle normale utgående bevegelser med primær refleksjonshastighet. Når f-k transformasjonen på fig. 16- blir konvertert tilbake til en normal seismisk seksjon a convolution of the transform and the Nyquist wavenumber. When this is filtered by decrementing all the samples except those representing the flat, 0-drop, reflections, the f-k transformation of fig. 16. It should be noted that this is analogous to treating normal outgoing movements with primary reflection velocity. When the f-k transformation in fig. 16- is converted back to a normal seismic section

tilveiebringes, fig. 17. Kun de flate refleksjonene blir tilbake. Dette eksempel demonstrerer følsomheten til f-k 0-ings filteret i tilfeller som er tett sammen ved fall. Dette er svært vanskelig å tilveiebringe ved tidsdomene-filtrering. is provided, fig. 17. Only the flat reflections remain. This example demonstrates the sensitivity of the f-k 0-ings filter in cases close together by drop. This is very difficult to provide with time domain filtering.

Fig. 13 viser et filter som vil vrake flate refleksjoner. Filteret 13 vil vrake alle sampler i området som har de lineære grensene 70-73. Dette har den virkning at de flate refleksjonene 50a....56a på fig. 7 blir filtrert ut. Som kjent blir refleksjoner korrigert ved flerhastighet ha en 0-ut bevegelse pr. spor^Foreliggende oppfinnelse frembringer derfor en effektiv teknikk for skarp diskriminering mot flere refleksjoner. Grensene 71 er igjen en liten endedel av det høyeste bølgetallet, +a/2Ax og grensen 73 er -a/2Ax. Fig. 13 shows a filter that will destroy flat reflections. The filter 13 will discard all samples in the range that have the linear boundaries 70-73. This has the effect that the flat reflections 50a....56a in fig. 7 is filtered out. As is known, reflections corrected at multiple speeds have a 0-out movement per trace^The present invention therefore provides an effective technique for sharp discrimination against multiple reflections. The limits 71 are again a small end part of the highest wave number, +a/2Ax and the limit 73 is -a/2Ax.

Når f-k transformasjonen på fig. 15 blir filtrert med filteret vist på fig. 13 fremkommer seksjonen på fig. 18. Alle samplene i området definert av grensene 70-72 har blitt vraket. Når denne rekken blir.transformert tilbake til en x-t seksjon fremkommer fig. 19. Denne viser vrakningen av de flate, refleksjonene, som kan være flere refleksjoner. Dette, viser igjen følsomheten til f-k 0-ings filteret i tilfeller som er tett sammen i fallet. When the f-k transformation in fig. 15 is filtered with the filter shown in fig. 13 the section on fig. 18. All samples in the area defined by boundaries 70-72 have been scrapped. When this series is transformed back into an x-t section, fig. 19. This shows the wreckage of the flat, reflections, which can be several reflections. This again shows the sensitivity of the f-k 0-ings filter in cases that are close together in the case.

Fig. 20-2 3 er nyttige for demonstrering av forbedringen av filteret i følge foreliggende oppfinnelse i forhold til hastighetsfilteret av den typen vanligvis henvist til som "kakastykke". Fig. 20 viser en seismisk seksjonsmodell og fig. 2.3 viser dens f-k transformasjon. De bratteste fall-refleksjonene blir filtrert ut ved å vrake alle samplene i det lukkede området som har grensene 74a-74b, 75a-75b og 76. Når en slik seksjon er filtrert på denne måten, fremkommer den seismiske seksjonen vist på fig. 21. Figs. 20-2 3 are useful for demonstrating the improvement of the filter of the present invention over the velocity filter of the type commonly referred to as "piece of cake". Fig. 20 shows a seismic section model and fig. 2.3 shows its f-k transformation. The steepest dip reflections are filtered out by discarding all the samples in the closed area having the boundaries 74a-74b, 75a-75b and 76. When such a section is filtered in this way, the seismic section shown in fig. 21.

Det skal bemerkes at de bratteste fall-refleksjonene erIt should be noted that the steepest dip reflections are

ikke tilstede i fig. 21.not present in fig. 21.

Når en slik seksjon blir filtrert med tidligere kjent domene "kakestykke" filter er virkningen at alle samplene i det kake-formede området definert av grensene 74a og 75a på fig. 23 blir vraket. Resultatet av en slik kakestykke-filtrering er vist på fig. 22. Lav-frekvens komponenter av de bratteste fall^-refleksjonene er fjernet, men refleksjonene fremkommer fremdeles i seksjonen på fig. 22 pga. høye frekvenskomponenter. I samsvar med foreliggende oppfinnelse blir refleksjonene fullstendig fjernet ved vraking av alle samplene som representerer den refleksjonen innbefattet i den aliaeserte delen innenfor grensen 74b, 75b og 76 på fig. 23. When such a section is filtered with the previously known domain "pie piece" filter, the effect is that all the samples in the pie-shaped area defined by the boundaries 74a and 75a in FIG. 23 is wrecked. The result of such pie-piece filtering is shown in fig. 22. Low-frequency components of the steepest fall^ reflections have been removed, but the reflections still appear in the section of fig. 22 due to high frequency components. In accordance with the present invention, the reflections are completely removed by scrapping all the samples representing that reflection included in the aliased portion within the boundaries 74b, 75b and 76 of FIG. 23.

For å forhindre aliasering slik som vist på fig. 23 må det horisontale samplirigsintervallet være mindre enn: To prevent aliasing as shown in fig. 23, the horizontal sampling interval must be less than:

hvor.t er den korteste perioden ved pulsen og 6t er utgangs-bevegelsestetthéten AT/AX til refleksjonen. Dersom dette mellomrommet ikke blir tilveiebragt må man filtrere seismo-grammet slik at alias frekvenskomponentene blir fjernet. where.t is the shortest period of the pulse and 6t is the output motion density AT/AX of the reflection. If this space is not provided, the seismogram must be filtered so that the alias frequency components are removed.

Alle frekvenser må bli vraket over f = (2Ax6t)~<l>hvor AxAll frequencies must be broken over f = (2Ax6t)~<l>where Ax

er det virkelige spormeilomrommet. f-k O-filtéret ifølge oppfinnelsen tillater vraking av hele refleksjonen ved inn-befatting av den aliaserte komponenten i et lukket område is the real track space. The f-k O filter according to the invention allows the wrecking of the entire reflection by including the aliased component in a closed area

av vrakede sampler.of wrecked samples.

Som et eksempel for signal/støyøkning vises til feltopp-tegningene vist på fig. 24 og 25. Feltseksjonen på fig. 24 As an example of signal/noise increase, reference is made to the field plan drawings shown in fig. 24 and 25. The field section in fig. 24

ble omformet til f-k seksjonen vist på fig. 26a og filtrert med filteret vist på fig. 26. Dette filteret vraker alle sampler i det lukkede området definert av grensene 77-93. was reshaped into the f-k section shown in fig. 26a and filtered with the filter shown in fig. 26. This filter discards all samples in the closed range defined by limits 77-93.

Det skal bemerkes at den virkelige gjennomføringen avIt should be noted that the actual implementation of

grensene slik som 78, 84, 86, 88, 90 og 91 er av segmenter, idet hver av de har en konstant sekvens eller et konstant bølge-tall. På denne måten er alle grensene til filteret i dets the boundaries such as 78, 84, 86, 88, 90 and 91 are of segments, each having a constant sequence or a constant wavenumber. In this way, all the limits of the filter are in its

virkelige Utførelse konstante bølgetall eller konstante frekvenser. real Execution constant wave numbers or constant frequencies.

Når feltopptegningen på fig. 24 er blitt filtrert med filteret på fig. 26 og konvertert tilbake til x-t seksjon, blir fig., 27 frembragt. Av fig. 27 fremgår refleksjonsgrenseflaten 94 klart. Av feltopptegningen på fig. 24 fremgår tilsvarende refleksjoner 94a klart synlig for en trenet geofysiker. Dette er et godt eksempel på den effektive filtreringevnen til foreliggende oppfinnelse. When the field plot on fig. 24 has been filtered with the filter in fig. 26 and converted back to x-t section, Fig. 27 is produced. From fig. 27, the reflection boundary surface 94 is clearly visible. From the field drawing on fig. 24 corresponding reflections 94a are clearly visible to a trained geophysicist. This is a good example of the effective filtering ability of the present invention.

Når feltseksjonen 25 er filtrert ved hjelp av et f-k filter ifølge foreliggende oppfinnelse, blir den seismiske seksjonen vist på fig. 28 frembragt. Filtreringsoperasjonen øker be-traktelig igjen tydeligheten av jordlagene. Fig. 24-28 viser et godt eksempel på signal/støyforholdøkning som kan tilveiebringes ved filtrering i samsvar med foreliggende oppfinnelse. When the field section 25 is filtered by means of an f-k filter according to the present invention, the seismic section is shown in fig. 28 produced. The filtering operation considerably increases again the clarity of the soil layers. Fig. 24-28 show a good example of signal/noise ratio increase which can be provided by filtering in accordance with the present invention.

Fig. 5 viser et flytdiagram over en datamaskinutført behandling for utførelse av foreliggende oppfinnelse. Denne behandlingen gjør bruk av en kjent hurtig Fourier transformasjon (FFT). Denne hurtige Fourier transformasjonen vil operere i begge retninger, dvs. fra x-t til f-k og fra f-k til x-t. Både den reelle og imaginære seksjon må derfor bli tilført som en inngang. Som vist på fig. 5 representerer rekken med digitale sampler et x-t felt seismogram som blir tilført som antydet med henvisningstallet 100. Den imaginære inngangen til en rekke med 0'er er vist med henvisningstållet 101. FFT operasjonen er antydet med henvisningtallet 102. Denne omformet den. seismiske seksjonen til en rekke med digitale sampler som representerer amplituden som en funksjon av frekvensen og horisontalavstanden. Den reelle del av denne frekvens transformasjon er vist med henvisningstållet 103 og den imaginære delen er vist med henvisningstållet 104. Disse seksjonene blir multiplekset ved 105 for å omordne samplene. .Dvs. alle de første samplene blir samplet sammen og alle de andre samplene blir samlet sammen og bortetter. Dette frem-bringerden reelle f transformasjonen 106 og den imaginære f transformasjonen 107. Igjen blir FFT operasjonen utført som vist med henvisningstållet 108 for å transformere rekkene til en reell f-k transformasjon 109 og en imaginær- f-k transformasjon 110. Fig. 5 shows a flow diagram of a computer-executed process for carrying out the present invention. This processing makes use of a well-known fast Fourier transform (FFT). This fast Fourier transformation will operate in both directions, ie from x-t to f-k and from f-k to x-t. Both the real and imaginary sections must therefore be added as an input. As shown in fig. 5, the row of digital samples represents an x-t field seismogram which is fed as indicated by the reference number 100. The imaginary input to a row of 0's is shown by the reference number 101. The FFT operation is indicated by the reference number 102. This reshaped it. the seismic section into a series of digital samples representing the amplitude as a function of frequency and horizontal distance. The real part of this frequency transformation is shown by the reference number 103 and the imaginary part is shown by the reference number 104. These sections are multiplexed at 105 to rearrange the samples. .Ie all the first samples are sampled together and all the other samples are collected together and deleted. This produces the real f-transformation 106 and the imaginary f-transformation 107. Again, the FFT operation is performed as shown by reference number 108 to transform the rows into a real f-k transformation 109 and an imaginary f-k transformation 110.

Filtreringen blir utført ved å multiplisere hver sampel i rekken med en vektfunksjon W(i,j) og W(i,j). Den reelle delen av denne operasjonen er betegnet med 111 og den imaginære delen er betegnet med tallet 112. W(i,j), W(i,j) kan være 0 i alle delene av rekken innenfor det lukkede området til frekvens og bølgetall. F.eks. kan den være 0 for alle samplene 1 rekken innenfor grensene 22-29 på fig. 3. Ellers kan rekken W(i,j), W(i,j) være 1. De filtrerte seksjoner blir igjen til-ført den hurtige Fourier transformasjonen ved 113 for å transformere seksjonen tilbake til en., x-t seismisk seksjon vist The filtering is performed by multiplying each sample in the row with a weight function W(i,j) and W(i,j). The real part of this operation is denoted by 111 and the imaginary part is denoted by the number 112. W(i,j), W(i,j) can be 0 in all parts of the series within the closed range of frequency and wave number. E.g. can it be 0 for all the samples 1 row within the limits 22-29 in fig. 3. Otherwise, the row W(i,j), W(i,j) may be 1. The filtered sections are again fed to the fast Fourier transform at 113 to transform the section back to one., x-t seismic section shown

ved 114.at 114.

Behandlingen av foreliggende oppfinnelse innbefatter den fundamentale operasjonen av utførelsen av en hurtig Fourier transformasjon og en f-k 0-filtreringsoperasjon. Hver av disse operasjoner vil bli bedre forstått fra den følgende detaljerte beskrivelsen. The processing of the present invention includes the fundamental operation of performing a fast Fourier transform and an f-k 0 filtering operation. Each of these operations will be better understood from the following detailed description.

Den hurtige Fourier transformasjon modulen benyttet vedThe fast Fourier transformation module used by

én utførelsesform av oppfinnelsen er av Cooley-Tukey typen beskrevet i Cooley, J.W.; Tukey, J.W. "An Algorithm for Machine Calculation of Complex Fourier Series", Mathematical Computation, Vol. 19, 1965 sidene 297-301. Denne, transformasjon er et eksempel på en teknikk som kan benyttes ved seismisk databehandling. Følgende definisjoner for uttrykk-ene tillater denne bruken. one embodiment of the invention is of the Cooley-Tukey type described in Cooley, J.W.; Tukey, J.W. "An Algorithm for Machine Calculation of Complex Fourier Series", Mathematical Computation, Vol. 19, 1965 pages 297-301. This transformation is an example of a technique that can be used in seismic data processing. The following definitions for the expressions allow this use.

Fourier serien med en to-dimensjonal funksjon f(t,x( er gitt ved: hvor The Fourier series with a two-dimensional function f(t,x() is given by: where

n m n.m

to dimensjonserien F^,7- er kompleks i formen: the two dimensional series F^,7- is complex in the form:

1 Li 1 Li

hvor where

som kan bliuttrykt som et par funksjoner, amplituder og fase for de variable -,- . which can be expressed as a pair of functions, amplitudes and phase for the variables -,- .

TL TL

For seismisk databehandling er det tale om variablene x og t For seismic data processing, these are the variables x and t

og tilsvarende transformasjonsvariable k = og f = -and corresponding transformation variables k = and f = -

L T L. T

hvor L er perioden for x variable og T er perioden for den tidsvariable. where L is the period for the x variable and T is the period for the time variable.

t,x mellomrom . f-k mellomrom t,x space . f-k space

variable samplingintervall variable sampling interval

Nyquistpunkt. Nyquist point.

frekvens frequency

utgangsbevegélsestetthet hastighet Nøkkelparameter Utgangsbevege.lsestettheten kan være positiv eller negativ som følge av en positiv eller negativ hastighet i f-k mellomrommet, f-k transformasjonene blir derfor fremstilt over området output movement density velocity Key parameter The output movement density can be positive or negative as a result of a positive or negative velocity in the f-k space, the f-k transformations are therefore produced over the area

Ved geofysisk undersøkelse har en vilkårlig definisjon a<y>fallet eller AT/Ax blitt tatt i bruk. Denne vilkårligheten stemmer fra det faktum at seismiske seksjoner blir fremvist under en konvensjon som plasserer nord og østenden på linjene til høyre på den seismiske seksjonsfremvisningen uavhengig av retningen i hvilken linen ble skutt. f-k transformasjonene beskrevet i spesifikasjonene avhenger av ordenen som de seismiske sporene tilsynekommer på båndet og på hvilket de blir lagret, idet første sporet definerer "0" tall til x-aksen og påfølgende spor tilsynekommer i økende positiv x retning. Med denne konvensjonen tilfeller som faller nedover i positiv x blir sagt å ha et positivt fall. Fall oppover i positiv In geophysical surveying, an arbitrary definition a<y>fall or AT/Ax has been adopted. This arbitrariness stems from the fact that seismic sections are displayed under a convention that places the north and east ends of the lines to the right of the seismic section display regardless of the direction in which the line was shot. The f-k transformations described in the specifications depend on the order in which the seismic traces appear on the tape and on which they are stored, with the first trace defining "0" numbers to the x-axis and subsequent traces appearing in the increasing positive x direction. By this convention cases that fall downward in positive x are said to have a positive fall. Fall upwards in positive

x har negativt fall.x has a negative slope.

Filtreringen av f-k transformasjonene, angitt med trinnene 111 og 112 på fig. 5 blir utført ved en operasjon som vraker eller nuller ut alle samplene til en rekke innenfor et lukket område. Hovedlisten for et datamaskinprogramfor utførelse av denne operasjonen er gitt nedenfor. The filtering of the f-k transformations, indicated by steps 111 and 112 in fig. 5 is performed by an operation which discards or zeros out all the samples of a range within a closed range. The main listing for a computer program to perform this operation is given below.

PARAMETER LIST, X, X, NT, TO, TL, XIPARAMETER LIST, X, X, NT, TO, TL, XI

SEISDATA (FIRSY, LAST, Tl, T2, FILE1, FILE2)SEISDATA (FIRSY, LAST, Tl, T2, FILE1, FILE2)

PAD TRACES (STARTNO, VALUÉ1, ENDNO, VALUE2)PAD TRACES (STARTNO, VALUE1, ENDNO, VALUE2)

PAD Time (TSTARTN, TVALUE1, TENDN, TVALUE2)PAD Time (TSTARTN, TVALUE1, TENDN, TVALUE2)

TRANSFORM TIME (FILE2, FILE3, FILE4, START, STOP)TRANSFORM TIME (FILE2, FILE3, FILE4, START, STOP)

MULTIPLEX (FILE3, FILE4, FIRST, LAST, REAL)MULTIPLEX (FILE3, FILE4, FIRST, LAST, REAL)

MULTIPLEX (FILE4, FILE5, FIRST, LAST IM)MULTIPLEX (FILE4, FILE5, FIRST, LAST IM)

TRANSFORM SPACE (FILE5, FILE6, FILE7, FILE8, START, STOP) PLOT FK (X, T, f, K, LEVELS, KEY, FILE7, FILE8) TRANSFORM SPACE (FILE5, FILE6, FILE7, FILE8, START, STOP) PLOT FK (X, T, f, K, LEVELS, KEY, FILE7, FILE8)

FILTER (FR, Fl, ROW, COLUMN, KEY)FILTER (FR, Fl, ROW, COLUMN, KEY)

FILTER FK (FR, Fl, FILE7, FILE8, FILE9, FILElO, KEY)FILTER FK (FR, Fl, FILE7, FILE8, FILE9, FILElO, KEY)

PLOT FK (X, T, f, K, LEVELS, KEP, FILE9, FILElO)PLOT FK (X, T, f, K, LEVELS, KEP, FILE9, FILElO)

TRANSFORM K (FILE9, FILElO, FILE11, FILE12, STARTK, STOPK) MULTIPLEX (FILE11, FILE12, FIRST, LAST, REAL) TRANSFORM K (FILE9, FILElO, FILE11, FILE12, STARTK, STOPK) MULTIPLEX (FILE11, FILE12, FIRST, LAST, REAL)

MULTIPLEX (FILE12, FILE 14, FIRST, LAST, IM)MULTIPLEX (FILE12, FILE 14, FIRST, LAST, IM)

TRANSFORM F (FILE13, FILE14, FILE15, START, STOP)TRANSFORM F (FILE13, FILE14, FILE15, START, STOP)

SEISDATASEIS DATA

Tar data mellom spor FIRSTR og LASTR fra FILEl og'lasterTakes data between tracks FIRSTR and LASTR from FILEl and'loads

i FILE2.in FILE2.

PAD TRACESPAD TRACES

Utfyllingen ENDNO og STARTNO til sporene ved enden av FILE1 til VALUE1 og VALUE2 er vanligvis 0. The padding ENDNO and STARTNO to the slots at the end of FILE1 to VALUE1 and VALUE2 is usually 0.

PAD TIMEPAD HOUR

Utfyllingstiden før og etter TSTARTN og TENDN med verdiene TVALUE1 og TVALUE2 henholdsvis, er vanligvis 0. The fill time before and after TSTARTN and TENDN with the values TVALUE1 and TVALUE2 respectively, is usually 0.

TRANSFORM TIMETRANSFORM TIME

Utfører FOURIER transformasjon på findata fra FILE2 ogPerforms FOURIER transformation on fine data from FILE2 and

putter reelle resultater inn i FILE3 og imaginære i FILE4. puts real results into FILE3 and imaginary into FILE4.

MULTIPLEXMULTIPLEX

Gjenkjenningsfrekvenssekvensmessige serier for FILE3 ogRecognition frequency sequential series for FILE3 and

FILE4 inn i romsekvensmessige serier i FILE5 og FILE6. FILE4 into spatially sequential series in FILE5 and FILE6.

TRANSFORM SPACETRANSFORM SPACE

Utfører Fourier transform på mellomromsekvens serier medPerforms Fourier transform on space sequence series with

FILE 5 og 6 med reelle resultater i FILE7, imaginære i FILE8.. FILE 5 and 6 with real results in FILE7, imaginary in FILE8..

PLOTFKPLOT FK

Plotter f-k transdormasjoner.Plot f-k transformations.

FILTERFILTER

Konstruerer vektfunksjonen FR og Fl for reelle og imaginære serier i f-k. mellomrommet eller ved å innstille nøkkelen vil et filter bli konstruert for å operere kun størrelsen til f-k transformasjonen. Et andre valg er 1/ER+FI. Constructs the weight function FR and Fl for real and imaginary series in f-k. the space or by setting the key a filter will be constructed to operate only on the magnitude of the f-k transform. Another choice is 1/ER+FI.

Andre valg er mulige.Other choices are possible.

FILTER FKFILTER FK

Utfører produktet for FR, Fl elementene ved transformasjon på FILE 7 og. 8 med resultatene on FILE 9 og 10. Tasten blir benyttet for å signalisere valget som i filter. Performs the product for the FR, Fl elements by transformation on FILE 7 and. 8 with the results on FILE 9 and 10. The key is used to signal the selection as in filter.

PLOT FKPLOT FK

Plotter resulterende FILTER FK transformasjon.Plot resulting FILTER FK transformation.

TRANSFORM KTRANSFORM K

Utfører Fourier transformasjonen på K sekvensmessig serie med FILES 9 og 10 som anbringer resultatene i FILTER 11 og 12. Performs the Fourier transform on K sequential series with FILES 9 and 10 placing the results in FILTERS 11 and 12.

MULTIPLEXMULTIPLEX

Gjenkjenner mellomromsekvensseriene på FILTERS 13 og 14. Recognizes the space sequence series on FILTERS 13 and 14.

TRANSFORM FTRANSFORM F

Utfører Fourier transformasjon på frekvenssekvensmessig serier på FILTER 13 og 14 med påfølgende tidsekvensmessig serier på FILE 15. Performs Fourier transformation on frequency sequence-wise series on FILTER 13 and 14 with subsequent time sequence-wise series on FILE 15.

Resultatet av operasjonen er et tidsmellomrom eller frekvens-bølgetall-filteroperasjon på inngangsdataseksjonen eller opptegningen. The result of the operation is a time slot or frequency-wavenumber filter operation on the input data section or plot.

Ovenfornevnte er en Fortran liste av en serie med subrutiner, funksjoner eller moduler. Hver utfører en rekke med operasjoner i modulær form for-behandlingen. The above is a Fortran listing of a series of subroutines, functions or modules. Each performs a series of operations in modular form for the processing.

Fig. 29 viser ved en utførelsesform av oppfinnelsen filter-operasjonen som ble utført på et "Control. Data Corporation" datamaskinsystem, modell nr. 175/760, som har en 6600 CPU Fig. 29 shows, in one embodiment of the invention, the filter operation performed on a "Control. Data Corporation" computer system, Model No. 175/760, having a 6600 CPU

og en CYPER 175 CPU med følgende perifere utstyr, CDC 6250and a CYPER 175 CPU with the following peripherals, the CDC 6250

9 Track Tape Transporter, CDC MAP III Array Processor og CDC Extended Core Storage. Spesifikasjonene for systemet er: 9 Track Tape Transporter, CDC MAP III Array Processor and CDC Extended Core Storage. The specifications for the system are:

Sentral prosessor: 18 og 60 bit registerCentral processor: 18 and 60 bit register

Periferisk prosessor: 12 bit registerPeripheral processor: 12 bit register

Utvidet hukommelse (ESC): 100K 60 bit ordExtended memory (ESC): 100K 60 bit words

Sentral hukommelse: 131K 60 bit ord (6600)Central memory: 131K 60 bit words (6600)

262K 60 bit ord (C175)262K 60 bit words (C175)

Rekke prosessor: CDC MAP IIIProcessor range: CDC MAP III

Masselager: 16 model 844-41Bulk stock: 16 model 844-41

348 million ord348 million words

Bånddrev: 2-7 spor 556/800 BPITape drive: 2-7 tracks 556/800 BPI

4-9 spor 800/1600 BPI 30 - 9 spor 1600/6250 BPI 4-9 slots 800/1600 BPI 30 - 9 slots 1600/6250 BPI

2 -21 spor 356 BPI2 -21 tracks 356 BPI

Opptegnings-drawing

utstyr: 2 - modél 405 kortleserequipment: 2 - modél 405 card reader

1 - model 415 korthuller1 - model 415 card holes

4 - linjeskrivere4 - line printers

Programmeringen som er nødvendig for å utføre oppfinnelsen vil. fremgå for fagmann på området av det foregående og fra håndboken for den spesielle anvendte datamaskin. The programming necessary to carry out the invention will. apparent to a person skilled in the art from the foregoing and from the manual for the particular computer used.

Mens spesielle utførelsesformer har blitt vist og beskrevetWhile particular embodiments have been shown and described

er forskjellige modifikasjoner innenfor rammen av oppfinnelsen mulig. Medfølgende krav er derfor ment å dekke alle slike modifikasjoner. different modifications within the scope of the invention are possible. The accompanying requirements are therefore intended to cover all such modifications.

Claims

Method for seismic survey of earth formations, characterized by generating seismic signals representing the amplitude of seismic reflections as a function of time and distance along a survey line, transforming the signals into one f-k sequence of digital samples representing the amplitude as a function of frequency and the wavenumber, operation on the f-k series to perform geophysical processing operations,

and. production from the series that has been operated on, a seismic recording that has an increased presentation of the formations.

2. Method according to claim 1, characterized by operation on the f-k range including filtering by discarding all the samples defined by prescribed closed ranges of frequency and wavenumber..

3. Method according to claim 2, characterized in that the filtering is carried out by weighting all the samples, the weighting being dependent on signal and noise in the seismic signals.

4. Method according to claim 1, characterized in that the wrecked samples are outside the closed region having boundaries that are linear with functions of the frequency and the wavenumber on both sides of a line having a prescribed angle to wreck all the samples except those representing reflections from the formations that have a prescribed dip.

5. Method according to claim i, characterized in that the closed area has boundaries, where one is a constant wave number above 0 and the other is a constant wave number below 0 and that the areas of wrecked samples include all samples with approximately 0 wave numbers representing several reflections.

6. Method according to claim 5, characterized in that the constant wave numbers are approximately a/2Ax where a is a small part of the highest wave number and Ax is the horizontal distance along the survey line between the positions where the seismic signals are provided.

7. Method according to claim 1, characterized in that one of the limits is a constant wave number greater than 0 and another of the limits is a constant wave number less than 0 and that the area of wrecked samples includes all samples that have a wave number greater than a limit and less than another limit whereby only reflections that have 0 dip are represented in the seismic record.

8. Method according to claim 1, characterized in that the area which has linear boundaries which are segments ..

with constant frequency and constant wave number.

9. Method of release according to claim 1, characterized .

in that the filtering discards samples in an area, where the samples represent noise whereby the signal/noise ratio of the filtered section is improved.

10. Method according to claim 1, characterized by back-transforming the filtered series to display the reflection time as a reflection of the distance along a survey line, and in that the production of the seismic record includes plotting a seismic section of the reflection times as a function of the distance. •

11. Method according to claim 1, characterized by back-transformation of the filtered series to display the reflection times as a function of the distance along the survey line, and that the production of a seismic record includes storing a series of digital samples representing a seismic section with reflection times as a function of distance.