NO337238B1

NO337238B1 - Seismic interference removal method

Info

Publication number: NO337238B1
Application number: NO20130090A
Authority: NO
Inventors: Thomas Elboth; Sindre Jansen
Original assignee: Fugro Geoteam As
Priority date: 2013-01-15
Filing date: 2013-01-15
Publication date: 2016-02-22
Also published as: NO20130090A1

Description

Oppfinnelsens tekniske område Technical field of the invention

Foreliggende oppfinnelse vedrører en seismisk ...... ifølge ingressen til det medfølgende patentkrav 1. The present invention relates to a seismic ... according to the preamble to the accompanying patent claim 1.

Bakgrunn for oppfinnelsen Background for the invention

Innledning Introduction

Et ofte opptredende problem med seismiske data er tilstedeværelsen av ulike typer koherent støy og særlig marin seismisk interferens (Sl). Seismisk interferens oppleves når forskjellige seismiske fartøy opererer samtidig i tett nærhet. Om amplituden og/eller moveout av den seismiske interferensen overstiger visse grenser, så vil de opererende fartøyene normalt måtte begynne med tidsdeling. Uheldigvis er dette kostbart og kan også lede til signifikant forsinkelse i fullførelsen av søket. De opererende fartøyene har vært kjent for å begynne tidsdelingen ved avstander opp til 100 km. Som en grov retningslinje er imidlertid seismisk interferens ofte sett på som problematisk når fartøyene er nærmere enn 40 km, noe som ofte er tilfelle i travle sommersesonger offshore i Nord-Europa og i den meksikanske gulfen. A frequently occurring problem with seismic data is the presence of various types of coherent noise and especially marine seismic interference (Sl). Seismic interference is experienced when different seismic vessels operate simultaneously in close proximity. If the amplitude and/or moveout of the seismic interference exceeds certain limits, then the operating vessels will normally have to start time sharing. Unfortunately, this is costly and can also lead to a significant delay in the completion of the search. The operating vessels have been known to begin time sharing at distances of up to 100 km. However, as a rough guideline, seismic interference is often seen as problematic when the vessels are closer than 40 km, which is often the case during busy summer seasons offshore in northern Europe and in the Gulf of Mexico.

Generelt kan seismisk interferens og algoritmer for koherent støyfjerning bli klassifisert i to grupper. Den første gruppen er basert på erkjennelsen av at koherent energi i skudd-domenet ofte opptrer som vilkårlig støy i andre domener som den felles offset, eller det felles midtpunktområdet (common midpoint (CMP)) (Larner et al., 1983). Verktøyene for vilkårlig støyfjerning, slik som FX-prediksjonsfiltrering (Canales 1984, Soubaras 1994) eller terskelmetoder (Hawkins 1998, Houck2000, Elboth et al., 2010 eller Schonwille et al., 2008), blir deretter anvendt på dataen, før dataen blir sortert tilbake til skudddomenet. Denne Sl-dempingstilnærming er blitt anvendt av Akbulut et al. (1984), Huaien et al. (1989) og mer nylig av Gulunay In general, seismic interference and coherent noise removal algorithms can be classified into two groups. The first group is based on the recognition that coherent energy in the shot domain often appears as arbitrary noise in other domains such as the common offset, or the common midpoint area (CMP) (Larner et al., 1983). The arbitrary denoising tools, such as FX prediction filtering (Canales 1984, Soubaras 1994) or thresholding methods (Hawkins 1998, Houck2000, Elboth et al., 2010 or Schonwille et al., 2008), are then applied to the data, before the data are sorted back to the shooting domain. This Sl damping approach has been used by Akbulut et al. (1984), Huaien et al. (1989) and more recently by Gulunay

(2008). I tilfellet der skuddpunktinter-vallene til interfererende fartøy imidlertid er synkronisert, opptrer seismisk interferens på de samme ankomsttidene i påfølgende skuddsamlinger. Seismisk interferens kan derfor ikke bli randomisert når de er sortert til andre domener og hele tilnærmingen bryter sammen. (2008). However, in the case where the shot point intervals of interfering vessels are synchronized, seismic interference occurs at the same arrival times in successive shot collections. Seismic interference cannot therefore be randomized when sorted into other domains and the whole approach breaks down.

Andre begrensninger vedrører fremgangsmåter som involverer ampli-tudeterskling. De krever at Sl-amplitudene må være større enn refleksjonsampli-tudene, hvilket ikke alltid er anvendbart. Other limitations relate to methods involving amplitude thresholding. They require that the Sl amplitudes must be greater than the reflection amplitudes, which is not always applicable.

Det er kjent å estimere et vektorfelt for refleksjonsmoveout ved å anvende kunnskapen om undersøkelsesgeometrien i kombinasjon med estimerte grunnformasjonshastigheter. It is known to estimate a vector field for reflection moveout by applying the knowledge of the survey geometry in combination with estimated basic formation velocities.

Den andre gruppen med Sl-fjerningsverktøy er basert på støymodellering og subtraksjon. Som tidlig eksempel er Kirlin og Done (1990), som anvender singulær verdidekomposisjon for å identifisere koherente hendelser i dataen for deretter å subtrahere disse. Mer nylige tilnærminger estimerer til slutt kildeposisjonen og/eller avfyringstidspunktet til den seismiske interferensen. Den seismiske interferensen blir deretter modellert og subtrahert, slik som i Warner et al. (2004); Brittan et al. (2008) og Kommedal et al. (2006). Suksessen med disse fremgangsmåtene avhenger sterkt av evnen til å bygge opp en nøyaktig modell av den seismiske interferensen. The second group of Sl removal tools is based on noise modeling and subtraction. An early example is Kirlin and Done (1990), who use singular value decomposition to identify coherent events in the data and then subtract these. More recent approaches ultimately estimate the source position and/or firing time of the seismic interference. The seismic interference is then modeled and subtracted, as in Warner et al. (2004); Brittan et al. (2008) and Kommedal et al. (2006). The success of these methods strongly depends on the ability to build up an accurate model of the seismic interference.

Fra den kjente teknikk skal det vises til US 2011/0276274 vedrørende en fremgangsmåte for prosessering av seismiske data, omfattende å tilveiebringe en samling av seismiske traser med seismiske forstyrrelser, dele den seismiske trasen i flerhet av vinduer, og bestemme seismiske forstyrelsesutbredelsesestimat. From the prior art, reference should be made to US 2011/0276274 regarding a method for processing seismic data, comprising providing a collection of seismic traces with seismic disturbances, dividing the seismic trace into a plurality of windows, and determining seismic disturbance propagation estimates.

Denne oppfinnelsen presenterer en fremgangsmåte for en automatisk estimering av Sl-movouten, hvilken deretter blir anvendt som et innledende trinn i to teknikker. Den første teknikken består av den automatiske genereringen av en tau-p demping, mens den andre teknikken anvender en fremgangsmåte for linjeintegralkonvolusjon og et estimat av vektorfelt for å separere seismisk interferens fra den seismiske refleksjonsdataen. This invention presents a method for an automatic estimation of the Sl moveout, which is then used as an initial step in two techniques. The first technique consists of the automatic generation of a tau-p damping, while the second technique uses a line integral convolution method and a vector field estimate to separate seismic interference from the seismic reflection data.

Oppsummering av oppfinnelsen Summary of the invention

Foreliggende oppfinnelse tilveiebringer en fremgangsmåte for prosessering av seismiske undersøkelsesdata, omfattende å tilveiebringe en original innsamling av seismiske spor med seismisk interferens, og å dele den originale innsamlingen opp i et flertall romvinduer; å bestemme for hvert av nevnte romvinduer et vektorfelt som er indikativt for moveouten av den seismiske interferens og et estimat av seismisk interferensmoveout. The present invention provides a method for processing seismic survey data, comprising providing an original collection of seismic traces with seismic interference, and dividing the original collection into a plurality of spatial windows; determining for each of said room windows a vector field indicative of the moveout of the seismic interference and an estimate of the seismic interference moveout.

I fremgangsmåten for prosesseringen av den innhentede dataen kan nevnte fastlegging inkludere et estimat av seismisk interferensmoveout inkludere sett med lokale moveouts av respektive enkelte av nevnte spor i hvert av nevnte romvinduer og fastlegging av et sett med lokale moveouts som skiller seg fra sett med lokale moveouts assosiert med seismiske undersøkelsesreflek-sjoner. In the method of processing the acquired data, said determination may include an estimate of seismic interference moveout, include sets of local moveouts of respective individual of said tracks in each of said spatial windows, and determination of a set of local moveouts that differ from sets of local moveouts associated with seismic survey reflections.

I fremgangsmåten for prosessering av seismiske undersøkelsesdata, kan nevnte fastlegging av et estimat av en seismisk interferensmoveout av moveout-vektorer i hvert av nevnte romvinduer, og fastlegging av et sett med lokale moveouts som skiller seg fra sett med lokale moveouts assosiert med seismiske undersøkelsesrefleksjoner i nevnte romvinduer. In the method of processing seismic survey data, said determining an estimate of a seismic interference moveout of moveout vectors in each of said spatial windows, and determining a set of local moveouts that differ from sets of local moveouts associated with seismic survey reflections in said room windows.

I fremgangsmåten for prosessering av seismiske undersøkelsesdata, kan nevnte fastlegging inkludere å identifisere på tvers av nevnte romvindu et sett med lokale moveoutkomponenter av moveoutvektorer som i det vesentlige har konstante moveoutverdier i rom og tid. In the method for processing seismic survey data, said determination may include identifying across said spatial window a set of local moveout components of moveout vectors having substantially constant moveout values in space and time.

I fremgangsmåten for prosessering av seismiske undersøkelsesdata kan et estimat av en seismisk interferensemoveout inkludere identifisering av den laveste relative statistiske spredning over nevnte romvindu. In the method of processing seismic survey data, an estimate of a seismic interference moveout may include identification of the lowest relative statistical spread over said spatial window.

I fremgangsmåten for prosessering av seismisk undersøkelsesdata kan nevnte relative statistiske spredning bli et normalisert standardavvik over nevnte romvinduer. In the method for processing seismic survey data, said relative statistical spread can become a normalized standard deviation over said spatial windows.

I fremgangsmåten for prosessering av seismiske undersøkelsesdata, kan nevnte fastlegging av settene med lokale moveouts for hvert av nevnte romvinduer inkludere krysskorrelerte par av nevnte spor i et glidende vindu for hvert av nevnte romvinduer, og glidning av nevnte glidende vindu i tid og rom over hvert av de respektive en av nevnte romvinduer. In the method of processing seismic survey data, said determination of the sets of local moveouts for each of said spatial windows may include cross-correlated pairs of said tracks in a sliding window for each of said spatial windows, and sliding said sliding window in time and space over each of the respective one of said room windows.

Foreliggende oppfinnelse tilveiebringer en første fremgangsmåte for å fjerne seismisk interferens i en original samling av seismisk data med seismisk interferens, omfattende The present invention provides a first method for removing seismic interference in an original collection of seismic data with seismic interference, comprising

- å fremskaffe et estimat ( p) av seismisk interferens moveouts generert på basis av et vektorfelt for seismisk interferensmoveout som er generert fra nevnte originale samling, - å anvende en tau-p transformasjon på nevnte originale samling for å generere en tau-p transformert samling, - å bestemme en tau-p demping i overensstemmelse med nevnte estimat av seismisk interferensmoveout - å anvende nevnte demping av tau-p på nevnte tau-p-transformerte samling for å generere en delvis dempet tau-p-transformert samling, og å anvende en invers tau-p-transformasjon på nevnte delvis dempede tau-p-transformerte samling for å generere en delvis dempet samling. - to obtain an estimate (p) of seismic interference moveouts generated on the basis of a seismic interference moveout vector field generated from said original collection, - to apply a tau-p transformation to said original collection to generate a tau-p transformed collection , - to determine a tau-p damping in accordance with said estimate of seismic interference moveout - to apply said damping of tau-p to said tau-p-transformed assembly to generate a partially damped tau-p-transformed assembly, and to apply an inverse tau-p transform on said partially damped tau-p transformed collection to generate a partially damped collection.

I nevnte første fremgangsmåte for å fjerne seismisk interferens i en original samling av seismiske spor med seismisk interferens, inkluderer fast-legging av nevnte tau-p demping å bestemme et første område i tau-p-domenet som inkluderer et estimat ( p) og nevnte anvendelse av nevnte tau-p-demping på nevnte tau-p-demping på nevnte tau-p-transformerte samling inkluderer å dempe ut data i tau-p-transformerte samling som er på utsiden av nevnte første område for å generere nevnte delvis dempede tau-p-transformerte samling, og å subtrahere nevnte delvis dempede samling fra nevnte originale samling, for å generere en samling med redusert seismisk interferens. In said first method of removing seismic interference in an original collection of seismic traces with seismic interference, determining said tau-p attenuation includes determining a first region in the tau-p domain that includes an estimate ( p ) and said applying said tau-p damping to said tau-p damping to said tau-p-transformed array includes damping out data in tau-p-transformed array that is outside of said first region to generate said partially damped tau -p-transformed assembly, and subtracting said partially damped assembly from said original assembly, to generate an assembly with reduced seismic interference.

I nevnte første fremgangsmåte for å fjerne seismisk interferens i en original samling av seismiske spor med seismisk interferens, inkluderer determinering av nevne tau-p-demping å determinere et første område i tau-p-domenet som inkluderer estimatet ( p ) av interferensmoveouten og nevnte anvendelse av nevnte tau-p-demping på nevnte tau-p-transformerte samling inkluderer å dempe ut dataen i den tau-p-transformerte samling som er på innsiden i det første området for å generere nevnte delvis dempede tau-p-transformerte samling. In said first method of removing seismic interference in an original collection of seismic traces with seismic interference, determining said tau-p attenuation includes determining a first region in the tau-p domain that includes the estimate ( p ) of the interference moveout and said applying said tau-p damping to said tau-p-transformed collection includes damping out the data in the tau-p-transformed collection that is inside the first region to generate said partially damped tau-p-transformed collection.

Nevnte første fremgangsmåte for fjerning av seismisk interferens i en original samling av seismiske spor med seismisk interferens, kan omfatte raffinering av tau-p-demping ved estimering av den sentrale ankomsttiden til den seismiske interferensen i skudd-domenet, for derigjennom å skaffe tilveie den sentrale ankomsttiden til den seismiske interferensen i tau-p-domenet. Said first method for removing seismic interference in an original collection of seismic traces with seismic interference may comprise refining tau-p damping by estimating the central arrival time of the seismic interference in the shot domain, thereby providing the central the arrival time of the seismic interference in the tau-p domain.

Foreliggende oppfinnelse fremskaffer en andre fremgangsmåte for å fjerne seismisk interferens i en original samling av seismiske spor med seismisk interferens, omfattende - å skaffe tilveie minst ett vektorfeltestimat av seismisk interferensmoveout på basis av nevnte originale samling - å anvende en avbildningsteknikk av typen linjeintegralkonvolusjon (Line integrated convolution (LIC)), derfor hvert sample av inputsamling (t, x) blir en lokal streamline, som starter ved denne sentersample, kalkulert i fremover og bakover retninger for 2L+1 sampler som følger vektorfeltet til nevnte seismiske interferensmoveout, hvorved seismisk interferens blir antatt å være koherent langs linjeintegralet, det adderer opp konstruktivt, og omvendt blir refleksjonshyperbolene som ikke er koherent over det samme linjeintegralet stakket ut ved å bruke LICen for effektiv filtrering av nevnte originale samling langs lokale strømningslinjer definert av nevnte seismiske vektorfelt for den seismiske interferensmoveout, for å generere et estimat av seismisk interferens med større mengder detaljer, og å subtrahere de nevnte utviklede estimat av seismiske interferens fra nevnte originale samling for å produsere en demping av den seismiske interferensen, der den seismiske interferensen i det vesentlige er blitt fjernet. The present invention provides a second method for removing seismic interference in an original collection of seismic traces with seismic interference, comprising - providing at least one vector field estimate of seismic interference moveout on the basis of said original collection - using an imaging technique of the type line integral convolution (Line integrated convolution (LIC)), therefore each sample of input collection (t, x) becomes a local streamline, starting at this center sample, calculated in forward and backward directions for 2L+1 samples following the vector field of said seismic interference moveout, whereby seismic interference becomes assumed to be coherent along the line integral, it adds up constructively, and conversely, the reflection hyperbolas which are not coherent over the same line integral are extracted using the LIC for efficient filtering of said original collection along local flow lines defined by said seismic vector field of the seismic interferer smoveout, to generate an estimate of seismic interference with a greater amount of detail, and to subtract said developed estimate of seismic interference from said original collection to produce a damping of the seismic interference, wherein the seismic interference has been substantially removed.

Foreliggende oppfinnelse skaffer til veie en tredje fremgangsmåte for fjerning av seismisk interferens i en original samling av seismiske spor med seismisk interferens, omfattende: The present invention provides a third method for removing seismic interference in an original collection of seismic traces with seismic interference, comprising:

- å skaffe tilveie minst ett vektorfeltestimat for en refleksjonshyperbel-moveout på basis av nevnte originale samling, - å bruke en avbildningsteknikk av typen linjeintegralkonvolusjon (LIC), der en lokal strømlinje som starter ved denne sentersample for hver sample av inputsamling (t,x) blir kalkulert i retning forover og bakover for 2L+1-sampler, etterfulgt av nevnte vektorfelt for refleksjonshyperbelmoveout, hvorved refleksjonshyperbolene blir antatt å være koherent langs linjeintegralet, der de summer seg opp konstruktivt, og motsatt at den seismiske interferensen som ikke er koherent over det samme linjeintegral blir stakket ut, - å benytte LICen for effektiv filtrering av nevnte originale samling langs lokale strømlinjer definert av nevnte vektorfelt for refleksjonshyper-belmoveout, å generere en Sl-dempet samling med store mengde detaljer, der den seismiske interferensen er blitt i det vesentlig fjernet. - to provide at least one vector field estimate for a reflection hyperbola moveout based on said original collection, - to use a line integral convolution (LIC) type imaging technique, where a local streamline starting at this center sample for each sample of input collection (t,x) is calculated in the forward and backward directions for 2L+1 samples, followed by the aforementioned vector field for reflection hyperbola moveout, whereby the reflection hyperbolas are assumed to be coherent along the line integral, where they add up constructively, and opposite that the seismic interference which is not coherent above it the same line integral is extracted, - to use the LIC for efficient filtering of said original collection along local streamlines defined by said vector field for reflection hyper-belmoveout, to generate an Sl-attenuated collection with a large amount of detail, where the seismic interference has become essentially removed.

Foreliggende oppfinnelse tilveiebringer et softwareprodukt for computere omfattende computerinstruksjoner beregnet på å gjennomføre en hvilken som helst av ovennevnte beskrevne fremgangsmåter. The present invention provides a software product for computers comprising computer instructions intended to carry out any of the above described methods.

Foreliggende oppfinnelse tilveiebringer et softwareprodukt for computere som omfatter computerinstruksjoner tilpasset for å utføre en hvilken som helst utførelsesform av ovennevnte beskrevne fremgangsmåter. The present invention provides a software product for computers comprising computer instructions adapted to perform any embodiment of the above described methods.

Foreliggende oppfinnelse tilveiebringer et apparat som omfatter funksjonelle elementer tilpasset for å operere ifølge et hvilket som helst av ovennevnte beskrevne fremgangsmåter. The present invention provides an apparatus comprising functional elements adapted to operate according to any of the above described methods.

Foreliggende oppfinnelse tilveiebringer et apparat som omfatter funksjonelle elementer tilpasset for å operere ifølge en hvilken som helst utførelsesform av ovennevnte beskrevne fremgangsmåter. The present invention provides an apparatus comprising functional elements adapted to operate according to any embodiment of the above described methods.

Detaljert beskrivelse av oppfinnelsen Detailed description of the invention

I det følgende vil foreliggende oppfinnelse bli beskrevet i form av et eksempel og med henvisning til de medfølgende tegninger der: figur 1 er en skjematisk tegning som viser et ideelt eksempel på tre syntetiske spor med en overlagret refleksjonshendelse. Ved å kryss-korrelere nabosporene blir et estimat for hvor mye de blir forsinket i forhold til hverandre i tid (vt) produsert. I en aktuell undersøkelse blir sporavstanden (vx) i den horisontale retningen gitt av sensorlayouten og er typisk konstant. Som illustrert med vektorene i figuren kan moveouten av en lokal refleksjonshendelse derfor bli representert av en vektor v = In the following, the present invention will be described in the form of an example and with reference to the accompanying drawings where: figure 1 is a schematic drawing showing an ideal example of three synthetic tracks with a superimposed reflection event. By cross-correlating the neighboring traces, an estimate of how much they are delayed relative to each other in time (vt) is produced. In a current survey, the track distance (vx) in the horizontal direction is given by the sensor layout and is typically constant. As illustrated with the vectors in the figure, the moveout of a local reflection event can therefore be represented by a vector v =

(vxvt). Ved å beregne slike vektorer for hver sample i den seismiske samling, er et (vxvt). By calculating such vectors for each sample in the seismic collection, et

vektorfelt som nøyaktig representerer tids- og romvarierende moveout av data oppnådd; vector fields that accurately represent the time- and space-varying moveout of data obtained;

figur 2a er en første graf som illustrerer distribusjonen av vt-vektorkom-ponentene beregnet i åtte horisontal like størrelsesromvinduer i det seismiske skuddsamlingseksempel illustrert i figur 3a. I fjernoffsetten representert av de siste få horisontale romvinduene, er det åpenbart fra plottet at to forskjellige distribusjoner er representert, av hvilke én kommer fra refleksjonsdataen, og den andre fra den seismiske interferensen;. figure 2a is a first graph that illustrates the distribution of the vt vector components calculated in eight horizontal windows of equal size in the seismic shot collection example illustrated in figure 3a. In the far offset represented by the last few horizontal space windows, it is obvious from the plot that two different distributions are represented, one of which comes from the reflection data, and the other from the seismic interference;.

Figur 2b er en andre graf som illustrerer det relative standardavviketOjfor hver moveoutverdi beregnet ifølge eksempelet i figur 2a. I dette eksempelet blir minimum0\oppnådd ved en moveout på -1800, som fremskaffer et estimat av moveouten som kommer fra seismisk interferens, illustrert i figur 1; Figur 3a er en eksemplifisert representasjon av seismisk dataplot av et første eksempel på en faktisk skuddsamling som er forurenset av seismisk interferens. I den horisontale retningen er det~600 spor som representerer registreringer fra en 8 km lang seismisk kabel, mens den vertikale retningen representerer tid. I det øvre høyre hjørnet kan en se klart noen seismiske refleksjonsdata, mens de nesten horisontale stripene som ses i den nedre delen av avbildningen er seismisk interferens. Det er denne interferensen som vi ønsker å fjerne uten å skade den underliggende svake refleksjonsdataen. De vertikale linjene er bare trukket på toppen av avbildningen for å illustrere åtte horisontal romvinduene som er benyttet. I det seismiske dataplottet er signalstyrken representert ved nivået på mørkheten. Følgelig representerer lyse farger svake signaler, mens det mørke (svarte) representerer sterke signaler. Figur 3b er et plot av seismisk dataavstand forberedt på basis av en eksemplifisert representasjon av seismisk dataplott illustrert i figur 3a, som illustrerer hvordan ankomsttiden til den seismiske interferensen i en skuddsamling blir estimert. Denne ankomsttiden blir benyttet for å sette en foredlet demping i tau-p-domenet. X-en beskriver et minimum som tilnærmet representerer den sentrale tiden til seismisk interferens, f, som vist i figur 3a. Figur 3c er et seismisk dataplott av tau-p-domene, som illustrerer den seismiske skuddsamlingen i figur 3a etter at dataen har blitt transformert til tau-p-domenet. På den horisontale aksen har vi nå p-verdier som representerer moveout i millisekunder, ms, mens tau er et mål på tid langs den vertikale aksen. I dette plottet kan en observere den store flekken i midten. Denne flekken er der den seismiske interferensen er avbildet. Legg også merke til «x»-en som igjen representerer moveout-estimatet av den seismiske interferensen langs den horisontale aksen og den sentrale tiden til seismisk interferens langs den vertikale aksen. Boksen som representerer den lokaliserte Sl-dempingen er også vist. På en ikke-restriktiv måte viser figur 3c en rektangulær boks. Formen på boksen kan imidlertid fordelaktig bli tilpasset Sl-avbildningen i tau-p-domenet. Boksen kan følgelig være trapesformet eller asymmetrisk rundt x-merket i figuren. Formen blir definert ved å benytte moveoutestimatet av den seismiske interferensen og sentraltiden til det seismiske interferensestimatet. I det ikke-restriktive eksemplet til den rektangulære boksen, som vist i figur 3, kan p^ =p- Aplog/?max = p + A2, der Api og Ap2for eksempel kan være satt av brukeren basert på iterative tester for å finne optimale verdier. På samme måte kan de øvre og nedre grenser langs tidsaksen for eksempel være satt av brukeren basert på iterative tester for å finne optimale verdier. Figur 3d er et seismisk dataplott som illustrerer resultatet av anvendelse av fjerning av seismisk interferens ifølge foreliggende oppfinnelse fra det første eksemplet til en faktisk skuddsamling som er forurenset av seismisk interferens som er blitt illustrert i figur 3a. Figur 3e er en seismisk data som illustrerer forskjellen mellom det første eksempelet til en faktisk skuddsamling som er forurenset av seismisk interferens, som er blitt illustrert i figur 3a og resultatet av anvendelse av seismisk interferensfjerning ifølge foreliggende oppfinnelse, som er illustrert i figur 3d. I forskjellsplottet i figur 3e, kan bare seismisk interferens klart bli sett, mens ingen åpenbar seismisk refleksjonsdata er synlig. Slike forskjellsplott blir ofte benyttet som en essensiell del av en kvalitetskontroll som er gjennomført under den seismiske prosesseringen; Figurene 4a, 4b og 4c er seismiske dataplott assosiert med en ytterligere, andre faktisk eksemplifisert skuddsamling som er forurenset av seismisk interferens fra foran, i hvilket figur 4a er et plott av den originale skuddsamlingen i dette eksempelet; figur 4b er et plot av skuddsamlingen i dette eksempelet etter at seismisk interferens er blitt fjernet ved anvendelse av foreliggende oppfinnelse; og figur 4c er et forskjellsplott oppnådd fra data i plottet i figur 4a. I figur 4a er også en strømlinje som følger moveouten til den seismiske interferensen vist med stiplede linjer med et pilhode. Figurene 5a, 5b og 5c er seismiske dataplott assosiert med en ytterligere, andre faktiske eksemplifiserte skuddsamling som er forurenset av seismisk interferens fra akter, i hvilket figur 5a er et plot av den originale skuddsamlingen i dette eksemplet; figur 5b er et plott av skuddsamlingen i dette eksempelet etter at seismisk interferens er blitt fjernet ved anvendelse av foreliggende oppfinnelse; og figur 5c er et forskjellsplott oppnådd ved å plotte resultatene ved å subtrahere dataen i plottet i figur 5b fra dataen i plottet i figur 5a. Figur 6 er et flytskjema som illustrerer en fremgangsmåte ifølge foreliggende oppfinnelse for seismisk interferens, moveout estimering. Figur 7 er et flytskjema som illustrerer en første fremgangsmåte ifølge foreliggende oppfinnelse for å fjerne seismisk interferens fra en seismisk skuddsamling. Figur 8 er et flytskjema som illustrerer en andre fremgangsmåte ifølge foreliggende oppfinnelse for å fjerne seismisk interferens fra en seismisk skuddsamling. Figure 2b is a second graph illustrating the relative standard deviation 0j for each moveout value calculated according to the example in Figure 2a. In this example, the minimum 0\ is reached by a moveout of -1800, which provides an estimate of the moveout coming from seismic interference, illustrated in Figure 1; Figure 3a is an exemplified seismic data plot representation of a first example of an actual shot assemblage contaminated by seismic interference. In the horizontal direction there are ~600 traces representing records from an 8 km long seismic cable, while the vertical direction represents time. In the upper right corner one can clearly see some seismic reflection data, while the almost horizontal stripes seen in the lower part of the image are seismic interference. It is this interference that we want to remove without damaging the underlying weak reflection data. The vertical lines are only drawn at the top of the image to illustrate the eight horizontal room windows that are used. In the seismic data plot, the signal strength is represented by the level of darkness. Consequently, light colors represent weak signals, while the dark (black) represent strong signals. Figure 3b is a plot of seismic data distance prepared on the basis of an exemplified representation of the seismic data plot illustrated in Figure 3a, which illustrates how the arrival time of the seismic interference in a shot collection is estimated. This arrival time is used to set a refined damping in the tau-p domain. The X describes a minimum that approximately represents the central time of seismic interference, f, as shown in Figure 3a. Figure 3c is a tau-p domain seismic data plot, illustrating the seismic shot collection in Figure 3a after the data has been transformed to the tau-p domain. On the horizontal axis we now have p-values representing moveout in milliseconds, ms, while tau is a measure of time along the vertical axis. In this plot you can observe the large spot in the middle. This spot is where the seismic interference is imaged. Also note the "x" which again represents the moveout estimate of the seismic interference along the horizontal axis and the central time of seismic interference along the vertical axis. The box representing the localized Sl attenuation is also shown. In a non-restrictive manner, Figure 3c shows a rectangular box. However, the shape of the box can advantageously be adapted to the S1 image in the tau-p domain. The box can therefore be trapezoidal or asymmetrical around the x mark in the figure. The shape is defined by using the moveout estimate of the seismic interference and the central time of the seismic interference estimate. In the non-restrictive example of the rectangular box, as shown in Figure 3, p^ =p- Aplog/?max = p + A2, where Api and Ap2 can for example be set by the user based on iterative tests to find optimal values. Similarly, the upper and lower limits along the time axis can, for example, be set by the user based on iterative tests to find optimal values. Figure 3d is a seismic data plot illustrating the result of applying the seismic interference removal of the present invention from the first example to an actual shot collection contaminated by seismic interference illustrated in Figure 3a. Figure 3e is a seismic data illustrating the difference between the first example of an actual shot assemblage contaminated by seismic interference, which has been illustrated in Figure 3a and the result of applying seismic interference removal according to the present invention, which is illustrated in Figure 3d. In the difference plot in Figure 3e, only seismic interference can be clearly seen, while no obvious seismic reflection data is visible. Such difference plots are often used as an essential part of a quality control carried out during the seismic processing; Figures 4a, 4b and 4c are seismic data plots associated with a further, second actually exemplified shot assemblage contaminated by frontal seismic interference, in which Figure 4a is a plot of the original shot assemblage in this example; Figure 4b is a plot of the shot collection in this example after seismic interference has been removed using the present invention; and Figure 4c is a difference plot obtained from data in the plot of Figure 4a. In Figure 4a, a streamline that follows the moveout of the seismic interference is also shown by dashed lines with an arrowhead. Figures 5a, 5b and 5c are seismic data plots associated with an additional second actual exemplified shot assemblage contaminated by seismic interference from aft, in which Figure 5a is a plot of the original shot assemblage in this example; Figure 5b is a plot of the shot collection in this example after seismic interference has been removed using the present invention; and Figure 5c is a difference plot obtained by plotting the results of subtracting the data in the plot of Figure 5b from the data in the plot of Figure 5a. Figure 6 is a flowchart illustrating a method according to the present invention for seismic interference, moveout estimation. Figure 7 is a flowchart illustrating a first method according to the present invention for removing seismic interference from a seismic shot collection. Figure 8 is a flowchart illustrating a second method according to the present invention for removing seismic interference from a seismic shot collection.

Seismisk interferens Sl, moveout estimering Seismic interference Sl, moveout estimation

I en kommersiell setting der en prosessor må gå gjennom gigabytes av data trenger brukerinteraksjonen å bli minimert. In a commercial setting where a processor must go through gigabytes of data, user interaction needs to be minimized.

En fremgangsmåte blir presentert her, som er i stand til å fungere for individuelle skuddsamlinger og som er i stand til å estimere Sl-moveouten uten krav om brukerinteraksjon. A method is presented here which is able to work for individual shot collections and which is able to estimate the Sl moveout without requiring user interaction.

Denne fremgangsmåten for Sl-moveout kan bli benyttet videre ved teknikker for Sl-fjerning, noe som vil bli forklart i senere seksjoner i denne beskrivelsen. This procedure for Sl moveout can be used further with techniques for Sl removal, which will be explained in later sections of this description.

I den følgende skal estimering av Sl-moveout ifølge oppfinnelsen bli forklart under henvisning til flytskjemaet i den medfølgende figur 6, og illustrasjonene i figurene 1, 2a, 2b og 3a. Først blir skuddsamlingen som inneholder M ganger N sampler delt opp i et antall J romvinduer der hvert romvindu inneholder alle samplene, og fortrinnsvis det samme antallet spor. Tallet J til romvinduer blir valgt på basis av kravet til oppløsning og prosesseffektivitet, og kan bli også valgt avhengig av antallet spor i skuddsamlingen. Det er observert at J valgt for å inkludere minst noen få titalls spor per vindu har gitt gode resultater. Hvert sporvindu blir scannet over skuddsamlingen ved å anvende et lite glidende vindu som strekker seg i tid og rom, der liten betyr at det glidende vinduet i retningen til spornummerdimensjonen blir dimensjonert for å dekke to spor som er lokalisert umiddelbart nær hver andre, eller to spor som er plassert nær nok hverandre til å inneholde seismisk interferens av interesse. For at det glidende vinduet skal dekke to spor som er lokalisert nær hverandre, men ikke ligge umiddelbart nær hverandre, vil følgelig det glidende vinduet bli deloppdelt i to delvinduer som tilsvarende er i avstand fra hverandre. I retningen til tidsret-ningen (sampler), er det glidende vinduet fortrinnsvis dimensjonert for å dekke nok sampler til å bestemme den normale moveouten til de faktiske undersøkelsesdataen, som i en faktisk undersøkelse typisk kunne være rundt ti sampler, og for å unngå senere ankomsttrekk av spor å bli inkludert ved å bestemme moveoutestimatet. I hver glidende vindusposisjon blir to spor, fortrinnsvis nabospor, krysskorrelert for å oppnå en vektor: v= (vxvt)T som representerer den lokale moveouten i dataen som er dekket av det glidende vinduet. vxer lik avstanden til sporene innenfor det glidende vinduet. I tilfellet der sporene er plassert umiddelbart ved siden av hverandre, er følgelig sporav-standen vxlik 1. Komponenten til tidsforsinkelsen vttil den lokale moveout-vektoren, typisk gitt i millisekunder (ms), blir estimert ved krysskorrelasjon til sporene som er inne I det glidende vinduet, eller inne I det glidende delvinduet, om ikke-tilstøtende sporer blir tatt I betraktning. De resulterende lokale moveoutvektorene beskriver et vektorfelt der vi en skuddsamling uten seismisk interferens gir den lokale moveout til refleksjonshyperblene til den innhentede dataen, men om seismisk interferens er til stede kan v også representere moveouten til seismisk interferens. Ved å scanne hele skuddsamlingen, blir en total av K-vektorer v deretter estimert for hvert romvindu, som vist i figur 1. Antallet med K-vektorer er gitt av størrelsen til det glidende vinduet og trinnet anvendt i scanningen for å bevege vinduet i sporvinduretningen og i tidsretningen (sample). Det skal anføres at om korrelasjonen av sporene inne i det glidende vinduet i en spesifikk lokasjon i skuddsamlingen ikke gir et brukbart resultat, blir ingen korresponderende vektor v bestemt for denne spesielle lokasjonen. In the following, estimation of Sl-moveout according to the invention will be explained with reference to the flowchart in the accompanying Figure 6, and the illustrations in Figures 1, 2a, 2b and 3a. First, the shot collection containing M by N samples is divided into a number of J spatial windows where each spatial window contains all the samples, and preferably the same number of tracks. The number J for room windows is chosen on the basis of the requirement for resolution and process efficiency, and can also be chosen depending on the number of tracks in the shot collection. It has been observed that J chosen to include at least a few tens of traces per window has given good results. Each track window is scanned over the shot collection using a small sliding window that spans time and space, where small means that the sliding window in the direction of the track number dimension is sized to cover two tracks that are located immediately close to each other, or two tracks which are placed close enough to each other to contain seismic interference of interest. In order for the sliding window to cover two tracks which are located close to each other, but not immediately close to each other, the sliding window will consequently be divided into two partial windows which are correspondingly at a distance from each other. In the direction of the time direction (samples), the sliding window is preferably sized to cover enough samples to determine the normal moveout of the actual survey data, which in an actual survey could typically be around ten samples, and to avoid later arrival moves of tracks to be included in determining the moveout estimate. At each sliding window position, two tracks, preferably neighboring tracks, are cross-correlated to obtain a vector: v= (vxvt)T representing the local moveout in the data covered by the sliding window. vxer equal to the distance of the tracks within the sliding window. In the case where the tracks are located immediately next to each other, the track spacing vx is therefore equal to 1. The component of the time delay vt to the local moveout vector, typically given in milliseconds (ms), is estimated by cross-correlation to the tracks contained in the sliding window, or inside the sliding partial window, if non-adjacent tracks are taken into account. The resulting local moveout vectors describe a vector field where we give a shot collection without seismic interference the local moveout to the reflection hyperbola of the acquired data, but if seismic interference is present v can also represent the moveout to seismic interference. By scanning the entire shot collection, a total of K-vectors v is then estimated for each spatial window, as shown in Figure 1. The number of K-vectors is given by the size of the sliding window and the step used in the scan to move the window in the track window direction and in the time direction (sample). It should be noted that if the correlation of the traces within the sliding window at a specific location in the shot collection does not yield a usable result, no corresponding vector v is determined for that particular location.

To forskjellige vektorfelt blir nå identifisert fra K-vektorene v. Én av to forskjellige vektorfelt indikerer moveouten til seismisk interferens, mens den andre indikerer moveouten til refleksjonshyperbelen til den seismiske undersøkelsen. Ifølge oppfinnelsen blir moveouten til Sl-vektorfeltet identifisert. Two different vector fields are now identified from the K vectors v. One of two different vector fields indicates the moveout to seismic interference, while the other indicates the moveout to the reflection hyperbola of the seismic survey. According to the invention, the moveout to the Sl vector field is identified.

Den seismiske interferensen ankommer typisk rimelig rettlinjet gjennom en skuddsamling, i kontrast til hva som er tilfelle med refleksjonshyperbelen. Ved å anvende denne observasjonen kan en vurdere distribusjonen av komponentene vtK-moveouten til vektorene vfor hvert romvindu, mens en tar hensyn til antagelsen at vektorene som er indikative for seismisk interferens, og som kommer fra en signifikant avstand, har en heller konstant moveout i rom og tid, mens vektorer som er indikative for refleksjonshyperbler, har større variasjoner. Særlig er vt-komponenten til den lokale moveoutvektor typisk mindre ved næroffset og øker med offset. Sl-moveouten er følgelig tatt i betraktning å være den med minimum statistisk spredning, fordelaktig det minimale relative standardavvik. I den seismiske interferensen, moveoutestimatet ifølge foreliggende oppfinnelse, blir en statistisk distribusjon vt,iav komponentene vttil vektoren vinne i hvert romvindu, beregnet. For den i<th->moveoutverdi til distribusjonen vt,i, blir standardavviket ov,itil antallet hendelser rijjf0r j=i, ■ ■ ■, J, kalkulert over den tilgjengelige J-representasjonen. Det relative standardavviket blir deretter beregnet for å oppnå det normaliserte standardavviket ou gitt her på en ikke-restriktiv måte i prosent av: The seismic interference typically arrives in a reasonably straight line through a shot collection, in contrast to what is the case with the reflection hyperbola. By applying this observation, one can assess the distribution of the components vtK-moveout of the vectors vfor each space window, while taking into account the assumption that the vectors which are indicative of seismic interference, and which come from a significant distance, have a rather constant moveout in space and time, while vectors indicative of reflection hyperbolas have greater variations. In particular, the vt component of the local moveout vector is typically smaller at near offset and increases with offset. The Sl moveout is therefore considered to be the one with the minimum statistical spread, preferably the minimum relative standard deviation. In the seismic interference, moveout estimate according to the present invention, a statistical distribution vt,iav of the components vttil the vector win in each spatial window is calculated. For the i<th->moveout value of the distribution vt,i, the standard deviation ov,itil the number of events rijjf0r j=i, ■ ■ ■, J, is calculated over the available J representation. The relative standard deviation is then calculated to obtain the normalized standard deviation ou given here in a non-restrictive manner as a percentage of:

der max («,,) er det maksimale antallet hendelser ny ved den i<th>-moveoutverdi. where max («,,) is the maximum number of events new at the i<th> moveout value.

Resultatet av beregningene av distribusjonene vt,ifor hver av et antall romvinduer er blitt illustrert i form av eksempel i figur 2a, som viser åtte grafer i respektive en av åtte romvinduer, av hvilke hver graf representerer distribu-sjonen av forekomsten av verdiene til komponenten vttil vektorene v som blir beregnet inne i dens respektive romvindu. The result of the calculations of the distributions vt,i for each of a number of room windows has been illustrated in the form of an example in Figure 2a, which shows eight graphs in each one of eight room windows, of which each graph represents the distribution of the occurrence of the values of the component vttil the vectors v which are calculated inside its respective spatial window.

For ytterligere å illustrere estimat av seismisk interferensmoveout ifølge foreliggende oppfinnelse skal nå virkelige dataeksempler vist i figur 3a bli vurdert, der seismisk interferens opptrer på grunn av et enkelt forstyrrende fartøy. Skuddsamlingen til 648 spor i figur 3a er blitt delt opp i J=8 romvinduer som hver inneholder 81 spor. Distribusjonen av komponentene vttil vektorene v for hvert romvindu er vist i figur 2a. Figur 2b viser det relative standardavviket0\for hver moveoutverdi. Minimumet o, er oppnådd ved en moveoutverdi som er gitt av: To further illustrate the estimate of seismic interference moveout according to the present invention, real data examples shown in figure 3a will now be considered, where seismic interference occurs due to a single disturbing vessel. The shot collection of 648 tracks in Figure 3a has been divided into J=8 room windows, each containing 81 tracks. The distribution of the components vt to the vectors v for each room window is shown in Figure 2a. Figure 2b shows the relative standard deviation0\for each moveout value. The minimum o, is achieved by a moveout value given by:

som i dette spesielle tilfellet (skuddsamlingen i figur 3a) er tilnærmet lik 1800 ms. which in this particular case (the collection of shots in Figure 3a) is approximately equal to 1800 ms.

For å validere dette estimatet blir dette datasettet transformer til tau-p-domenet, og moveouten på omlag 850 ms blir lokalisert (se figur 3c). Denne verdien ligger nært opp til estimatet p til Sl-moveout, som er vurdert å være tilstrekkelig nøyaktig. To validate this estimate, this data set is transformed into the tau-p domain, and the moveout of approximately 850 ms is located (see Figure 3c). This value is close to the estimate p for Sl-moveout, which is considered to be sufficiently accurate.

Til slutt lager vi et Sl-vektorfelt i skuddsamlingen v0<SI>ved å beholde bare vektorene med estimert moveout innen [ pmin /<?>max]-intervaller, dermin= p- Apl ogmax<=>p<+>A2. I figur 3c er A1-<Ap>2- 300 ms. Siden fremgangsmåte 2 krever en tett representasjon av vektorfeltet, er de manglende vektorene opp-nådd ved interpolering og blir gjennomsnittet over de eksisterende verdier, med fortrinnsvis én vektor per samplelokasjon. Denne tette representasjonen av Sl-vektorfeltet er benevnt som v<sl>. Finally, we create an Sl vector field in the shot collection v0<SI>by keeping only the vectors with estimated moveout within [ pmin /<?>max] intervals, wheremin= p- Apl andmax<= >p<+>A2. In figure 3c, A1-<Ap>2- is 300 ms. Since method 2 requires a dense representation of the vector field, the missing vectors are obtained by interpolation and are averaged over the existing values, preferably with one vector per sample location. This dense representation of the Sl vector field is named v<sl>.

Tilsvarende lager vi en refleksjonshyperbelvektorfelt i skuddsamlingen v0<R>ved å beholde bare vektorene som ikke tilhører til Sl-vektorfeltet v0<SI>. Siden fremgangsmåte 2 krever en tett representasjon av vektorfeltet, blir de mang-lende vektorene fremskaffet ved interpolasjon og blir gjennomsnittet over de eksiterende vektorene for å fremskaffe fortrinnsvis én vektor per sampleloka-sjon, Denne tette representasjonen av refleksjonshyperblenes vektorfelt er betegnet med V<R>. Similarly, we create a reflection hyperbola vector field in the shot collection v0<R> by keeping only the vectors that do not belong to the Sl vector field v0<SI>. Since method 2 requires a dense representation of the vector field, the missing vectors are provided by interpolation and are averaged over the exciting vectors to provide preferably one vector per sample location. This dense representation of the vector field of the reflection hyperbola is denoted by V<R> .

Fremgangsmåte 1 for Sl-fjerning: Automatisk tau-p demping. Method 1 for Sl removal: Automatic rope-p damping.

I det følgende skal den første fremgangsmåten for fjerning av seismisk interferens ifølge oppfinnelsen bli forklart. Med automatisk menes at fremgangsmåten kan bli gjennomført av en maskin, slik som for eksempel en computer gjennomfører et egnet skriftlig computerprogram uten behov for interaksjon av en menneskelig operatør. Den første fremgangsmåten for å fjerne seismisk interferens inkluderer generering av en tau-p-demping, basert på estimatet p av Sl-moveout og fortrinnsvis den sentrale tiden f. Etter å ha generert tau-p-dempingen som avgrenser området i tau-p-domenet der seismisk interferens har kartlagt, blir dataen i tau-p-domenet dempet ut tilsvarende for effektivt å fjerne eller redusere seismisk interferens i samlingen. In the following, the first method for removing seismic interference according to the invention will be explained. By automatic is meant that the procedure can be carried out by a machine, such as for example a computer carrying out a suitable written computer program without the need for interaction by a human operator. The first procedure to remove seismic interference includes the generation of a tau-p attenuation, based on the estimate p of Sl moveout and preferably the central time f. After generating the tau-p attenuation that delimits the area in the tau-p- the domain where seismic interference has mapped, the data in the tau-p domain is attenuated accordingly to effectively remove or reduce seismic interference in the collection.

En første utførelsesform av den første fremgangsmåten er forklart i det følgende under henvisning til flytskjemaet i den medfølgende figur 7. I denne første utførelsesform dempes dataen i tau-p-transformasjonsamlingen som ikke ligger i den genererte tau-p-dempingen. Deretter blir så en invers tau-p-transformasjon gjennomført på data som er igjen etter at tau-p-dempingen er blitt anvendt, for å generere en delvis dempet samling. Den delvis dempede samling blir deretter subtrahert eller adaptivt subtrahert fra den originale skuddsamling, for slik å fremskaffe en Sl-undertrykket skuddsamling (se eksempel illustrert i figur 3d). Subtraksjonen fremskaffer god preservering av refleksjonsdata, som vist i det illustrerende plottet i figur 3e. A first embodiment of the first method is explained below with reference to the flowchart in the accompanying Figure 7. In this first embodiment, the data in the tau-p transform collection that is not in the generated tau-p damping is damped. An inverse tau-p transformation is then performed on data remaining after the tau-p damping has been applied to generate a partially damped collection. The partially attenuated collection is then subtracted or adaptively subtracted from the original shot collection, so as to produce a Sl-suppressed shot collection (see example illustrated in Figure 3d). The subtraction provides good preservation of reflectance data, as shown in the illustrative plot in Figure 3e.

En andre utførelsesform av den første fremgangsmåten er forklart i det følgende. Denne andre utførelsesformen involverer fremover-invers tau-p-transformasjon. Den fremskaffer et resultat som kan avvike fra resultatene oppnådd med den første utførelsesformen, siden fremover-invers tau-p-transformasjon ikke er ansett som signalpreserverende. En tau-p-demping basert på estimatet av Sl-moveout, og fordelaktig, den sentrale tiden f. I denne andre utførelsesformen blir all data som ligger inne i den genererte tau-p-dempingen, dempet ut i tau-p-domenet. Deretter blir en invers tau-p-transformasjon gjennomført på dataen som blir igjen etter tau-p-dempingen anvendt for å generere en delvis dempet samling som fremskaffer en Sl-dempet skuddsamling. A second embodiment of the first method is explained in the following. This second embodiment involves the forward-inverse tau-p transform. It provides a result that may differ from the results obtained with the first embodiment, since the forward-inverse tau-p transform is not considered signal-preserving. A tau-p damping based on the estimate of Sl moveout, and advantageously, the central time f. In this second embodiment, all data contained in the generated tau-p damping is damped out in the tau-p domain. Then, an inverse tau-p transform is performed on the data remaining after the tau-p damping is applied to generate a partially damped collection which provides an S1 damped shot collection.

For å foredle tau-p-dempingen, blir også et estimat av den sentrale ankomsttiden til seismisk interferens i skudd-domenet f generert, for derigjennom å fremskaffe f i tau-p-domenet, for slik å muliggjøre en mer nøyaktig demping. For å gjøre dette, er kalkulasjonen for hvert sample av det første sporet (f 1 )T, gjort for den totale avstanden (gitt for eksempel i sampler) for alle vektorene (f,- x/)<T>av Sl-vektorfeltet v0<SI>, av; In order to refine the tau-p damping, an estimate of the central arrival time of seismic interference in the shot domain f is also generated, thereby providing f in the tau-p domain, thus enabling a more accurate damping. To do this, the calculation for each sample of the first trace (f 1 )T is made for the total distance (given for example in samples) of all the vectors (f,- x/)<T> of the Sl vector field v0 <SI>, of;

Deretter blir et estimat av f oppnådd ved minimum avstand dtot, ved Then an estimate of f is obtained at the minimum distance dtot, at

f sargmin^ f sargmin^

Fastleggingen av f for det reelle datasetteksempel er vist i figur 3b. Estimatet av f og p definerer en rimelig nøyaktig tau-p-demping vist som et rektangulært vindu i figur 3c. Tau-p-dempingen kunne bli designet som nevnt ovenfor, med andre former enn den eksemplifiserte rektangulære formen illustrert i figur 3c. Grensene eller begrensningene av arealet til dempingen kunne fordelaktig bli fastlagt på basis av erfaring oppnådd ved bruken av fremgangsmåten ifølge oppfinnelsen, eller ved å benytte prosesseringsteknikker av typen plottavbildning eller statistisk analyse av den tau-p-transformerte skuddsamling. The determination of f for the real data set example is shown in Figure 3b. The estimate of f and p defines a reasonably accurate tau-p damping shown as a rectangular window in Figure 3c. The Tau-p damping could be designed as mentioned above, with other shapes than the exemplified rectangular shape illustrated in Figure 3c. The limits or limitations of the area of the damping could advantageously be determined on the basis of experience gained in the use of the method according to the invention, or by using processing techniques such as plotting or statistical analysis of the tau-p-transformed shot collection.

Fremgangsmåte 2 for Sl-fjerning: Linjeintegralkonvolusjon Method 2 for Sl removal: Line integral convolution

I det følgende blir den andre fremgangsmåten for seismisk interferens, Sl-fjerning ifølge oppfinnelsen forklart med referanse til flytskjemaet i den medfølgende figur 8. Tilsvarende som for den første fremgangsmåten kan denne andre fremgangsmåten bli gjennomført av en maskin, slik som for eksempel en computer som gjennomfører et egnet skrevet computerprogram uten krav til interaksjon fra en operatør. Den andre fremgangsmåten benytter en fremgangsmåte referert til som linjeintegralkonvolusjon (Cabral og Leedom 1993). In the following, the second method for seismic interference, Sl removal according to the invention is explained with reference to the flowchart in the accompanying figure 8. Similar to the first method, this second method can be carried out by a machine, such as for example a computer which executes a suitable written computer program without requiring interaction from an operator. The second method uses a method referred to as line integral convolution (Cabral and Leedom 1993).

I en første utførelsesform av den andre fremgangsmåten ifølge oppfinnelsen er teknikken benevnt linjeintegralkonvolusjon LIC (= Line Integral Convolution) blitt tilpasset for Sl-fjerning, noe som tar fordel av det faktum at seismisk interferens kan forventes lokalt å være heller koherent. I denne første utførelsesformen, som i det følgende skal beskrives nærmere for hver sample av inputsamling (t, x), blir en lokal strømlinje som starter ved dette sentersample kalkulert i retning fremover og bakover for2L+1 sampler, etterfølgende inputvektorfeltet V<SI>som er den tette representasjonen av Sl-vektorfeltet. Outputverdien i (t, x) til kalkulasjonen er medianen eller gjennomsnittsverdien til alle amplitudene langs denne strømningslinjen. Siden seismisk interferens er antatt å være koherent langs linjeintegralet blir de addert opp konstruktivt. Motsatt er refleksjonshyperblene fra den seismiske dataen ikke koherent over det samme linjeintegral, og typisk bare stakket ut. Ifølge oppfinnelsen filtrerer LICen inputskuddsamlingen effektivt langs lokal strømlinjer definert av V<R>for å generere et estimat av den seismiske interferensen med store mengder detaljer. Ved å ha estimatet av det av den seismiske interferensen generert på denne måten, blir den subtrahert fra input skuddsamling å produsere en Sl-dempet samling, i hvilken den seismiske interferensen er blitt i alt vesentlig fjernet. In a first embodiment of the second method according to the invention, the technique called line integral convolution LIC (= Line Integral Convolution) has been adapted for Sl removal, which takes advantage of the fact that seismic interference can be expected locally to be rather coherent. In this first embodiment, which will be described in more detail in the following for each sample of the input collection (t, x), a local streamline starting at this center sample is calculated in the forward and backward direction for 2L+1 samples, following the input vector field V<SI> as is the dense representation of the Sl vector field. The output value in (t, x) of the calculation is the median or average value of all the amplitudes along this streamline. Since seismic interference is assumed to be coherent along the line integral, they are added up constructively. Conversely, the reflection hyperbola from the seismic data is not coherent over the same line integral, and typically only protrudes. According to the invention, the LIC effectively filters the input shot collection along local streamlines defined by V<R> to generate an estimate of the seismic interference with a large amount of detail. By having the estimate of that of the seismic interference generated in this way, it is subtracted from the input shot collection to produce an Sl damped collection, in which the seismic interference has been substantially removed.

Ifølge oppfinnelsen, også i en andre utførelsesform av den andre fremgangsmåten som er forklart i det følgende, blir fremgangsmåten referert til som linjeintegralkonvolusjon. Ifølge oppfinnelsen har teknikken med linjeintegralkonvolusjon, LIC, blitt tilpasset for Sl-fjerning, der en tar fordel av det faktum at seismisk interferens kan forventes lokalt å bli heller koherent. For hvert sample av inputsamlingen (t, x), blir en lokal strømlinje som starter ved denne sentersample i retningen framover og bakover for 2L+1 sampler, som følger inputvektorfelt v<R>som er den tette representasjonen av vektorfeltets refleksjonshyperbel. Outputverdien i (t, x) av kalkulasjonen er median- eller gjennomsnittsverdien for alle amplitudene langs denne strømlinjen. Siden refleksjonshyperblene er antatt å være koherente langs integrallinje, summerer de opp konstruktivt. Motsatt er ikke den seismiske interferensen fra seismisk data koherent over det samme linjeintegral, og typisk bare stakker opp. Ifølge oppfinnelsen filtrerer LICen effektivt inputskuddsamlingen langs lokale strømlinjer definert av V<R>| for å generere et estimat av Sl-dempede samling, der den seismiske interferensen er blitt i det vesentlige fjernet. According to the invention, also in a second embodiment of the second method which is explained in the following, the method is referred to as line integral convolution. According to the invention, the technique of line integral convolution, LIC, has been adapted for Sl removal, where one takes advantage of the fact that seismic interference can be expected locally to become rather coherent. For each sample of the input collection (t, x), a local streamline starting at this center sample in the forward and backward directions for 2L+1 samples follows the input vector field v<R> which is the dense representation of the vector field's reflection hyperbola. The output value in (t, x) of the calculation is the median or average value of all the amplitudes along this streamline. Since the reflection hyperbolas are assumed to be coherent along the integral line, they sum up constructively. Conversely, the seismic interference from seismic data is not coherent over the same line integral, and typically just stacks up. According to the invention, the LIC effectively filters the input shot collection along local streamlines defined by V<R>| to generate an estimate of Sl damped assemblage, where the seismic interference has been substantially removed.

Figurene 4a, b og c og figur 5a, b og c er relatert til hver sin av to reelle skuddsamlinger, kontaminert av seismisk interferens fra et enkelt forstyrrende fartøy: figurene 4a, b og c er relatert til en skuddsamling fra front, og figurene 5a, b og c er relatert til en skuddsamling med seismisk interferens fra akter. Begge skuddsamlingene, illustrert henholdsvis i figurene 4a og 5a, er prosessert med fremgangsmåte 2 for L=30 sampler. Outputsamlinger (henholdsvis figurene 4b og 5b) er rimelig rent, med seismisk interferens suksessfullt dempet, og forskjellsplottene (henholdsvis figurene 4c og 5c) av resultatet fra subtraksjonen viser at refleksjonshyperblene til den underliggende refleksjonsdataen er vel preservert. Figures 4a, b and c and figures 5a, b and c relate to each of two real shot collections, contaminated by seismic interference from a single disturbing vessel: figures 4a, b and c relate to a shot collection from the front, and figures 5a , b and c are related to a shot assemblage with seismic interference from aft. Both shot collections, illustrated respectively in figures 4a and 5a, are processed with method 2 for L=30 samples. Output assemblies (Figures 4b and 5b, respectively) are reasonably clean, with seismic interference successfully attenuated, and the difference plots (Figures 4c and 5c, respectively) of the result from the subtraction show that the reflection hyperbolas of the underlying reflection data are well preserved.

Fordeler Benefits

De to effektive fremgangsmåtene for Sl-fjerning basert på estimering av automatisk Sl-moveout er begge fullt ut datadrevet, og trenger bare anvendelse av én skuddsamling. Prosesseringstiden er kort og fremgangsmåten har derfor potensial for å kunne anvendes i sanntid etter hvert som undersøkelsen blir gjennomført. Siden ankomsttiden og retningen til seismisk interferens kan endres fra én skuddsamling til en annen skuddsamling, genererer den automatiske deteksjonen av attributter automatisk nøyaktige dempinger med liten interaksjon fra prosessoren. The two effective Sl removal methods based on automatic Sl moveout estimation are both fully data driven, requiring only the application of one shot collection. The processing time is short and the method therefore has the potential to be used in real time as the survey is carried out. Since the arrival time and direction of seismic interference can change from one shot collection to another shot collection, the automatic detection of attributes automatically generates accurate attenuations with little interaction from the processor.

Claims

1. Method for processing seismic survey data, comprising: obtaining an original collection of seismic traces with seismic interference, dividing the original collection into a plurality of spatial windows, determining for each of said spatial windows, a vector field indicative of moveout of seismic interference and to determine an estimate of seismic interference moveout.

2. Method according to claim 1, wherein said determination of an estimate for a seismic interference moveout includes: - determining sets of local moveouts for each respective one of said tracks in each of said space windows, and - detecting a set of local moveouts that deviate from set of local moveouts associated with the seismic survey reflections.

3. Method according to claim 2, wherein said determination of an estimate of a seismic interference moveout includes determining a distribution of local moveout components of moveout vectors in each of said spatial windows and detecting a set of local moveouts that deviate from the sets of local moveouts associated with seismic survey reflection in said room windows.

4. Method according to claim 3, where said detection includes identifying across said space windows a set of local moveout components of moveout vectors which have essentially constant moveout values in space and time.

5. Method according to claim 3 or 4, wherein said determination of an estimate of a seismic interference moveout includes identifying the lower relative statistical spread over said spatial windows.

6. Method according to claim 5, where said relative statistical spread is a normalized standard deviation over said room windows.

7. Method according to any one of claims 2-6, wherein said fixing set of local moveouts for each of said space windows includes cross-correlating pairs of said tracks in a sliding window for each of said space windows, and sliding said sliding window in time and room above its respective one of said room windows.

8. An apparatus comprising functional elements adapted to operate according to any one of claims 1 to 7.