NO330675B1 - Fremgangsmate for prosessering av et seismisk datasett - Google Patents

Fremgangsmate for prosessering av et seismisk datasett Download PDF

Info

Publication number
NO330675B1
NO330675B1 NO20041809A NO20041809A NO330675B1 NO 330675 B1 NO330675 B1 NO 330675B1 NO 20041809 A NO20041809 A NO 20041809A NO 20041809 A NO20041809 A NO 20041809A NO 330675 B1 NO330675 B1 NO 330675B1
Authority
NO
Norway
Prior art keywords
seismic data
coordinate
transformation
data set
inverse
Prior art date
Application number
NO20041809A
Other languages
English (en)
Other versions
NO20041809L (no
Inventor
Michel Albert Schonewille
Original Assignee
Pgs Exploration Uk Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Pgs Exploration Uk Ltd filed Critical Pgs Exploration Uk Ltd
Publication of NO20041809L publication Critical patent/NO20041809L/no
Publication of NO330675B1 publication Critical patent/NO330675B1/no

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362Effecting static or dynamic corrections; Stacking

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Control Of Vending Devices And Auxiliary Devices For Vending Devices (AREA)

Description

Bakgrunn for oppfinnelsen
1. Teknisk område
Foreliggende oppfinnelse angår generelt området geo-fysiske undersøkelser. Mer spesielt angår oppfinnelsen behandling av seismiske data. Spesielt er oppfinnelsen en fremgangsmåte for å korrigere for tidsforskyvninger i seismiske data som et resultat av azimutal variasjon.
2. Beskrivelse av beslektet teknikk
Nåværende seismiske datainnsamlinger fører til at seismiske data blir samplet uregelmessig langs romkoordinatene. For konvensjonelle seismiske datasett er disse koordinatene vanligvis linjemidtpunktet, tverrlinjemidtpunktet, offset og azimut. Denne uregelmessige samplingen kan generere problemer for tidsmessig seismikk og avbilding, innbefattende forstakk-avbildning. Den uregelmessige samplingen i midt-punktene og offset kan regulariseres ved å benytte konvensjonelle regulariserings- eller rekonstruksjons-teknikker, slik som Fourier-regularisering. Disse regulariserings-teknikkene kalkulerer eller anslår nye verdier for linjemidtpunktet, tverrlinjemidtpunktet og offset-variable slik at disse variable blir regelmessig samplet. Disse teknikkene som innbefatter Fourier-regularisering, kompenserer imidlertid ikke for variasjon i den azimutvariable. Likevel kan azimutale variasjoner ha stor innvirkning på behandlingen av seismiske data.
Hvis et enkelt fall-lag i en homogen undergrunn blir tatt i betraktning og to traser blir sammenlignet med det samme midtpunkt og absolutt offset, men med forskjellige azimut, så vil refleksjonshendelsen forskyves i tid. Tidsforskyvningen er avhengig av fallvinkelen og retningen til laget, hastigheten i undergrunnen, azimut og offset. Disse tidsforskyvningene som skyldes azimutvariasjoner, er vanligvis begrenset i størrelsesorden noen få millisekunder. For generelle seismiske dataavbildningsbehov vil neglisjering av disse tidsforskyvningene ha begrensede konsekvenser. For tidsforløpsdata hvor grunnlinje- og overvåknings-undersøkelser blir holdt fra hverandre, kan selv en tidsforskyvning på 4 ms føre til feil av samme størrelse som differansen i de målte signaler. For sterkt fallende lag, spesielt med skjeve fallretninger, er neglisjering av azimutvariasjonen ikke en effektiv løsning for tidsforløpsdata. Her er repeterbarhet essensiell.
Etter Fourier-regularisering blir linje- og krysslinje-midtpunktene og de absolutte offset i de regulariserte traser regelmessig samplet. For flere ytterligere behandlingsmetoder, slik som fallutflyttingskorreksjon og forstakkmigrering, er kilde- og mottaker-koordinatene nødvendige. For å utlede disse posisjonene fra midtpunkt- og offset-posisjonene, er det nødvendig med en azimut. En løsning er å anta at azimut er null i forhold til linjeretningen. Dette er seilingsretningen i en marin seismisk undersøkelse. Antagelse av null azimut i forhold til seilingsretningen er imidlertid ikke ideell fordi det seismiske signalet vil være avhengig av azimut. I prinsippet, hvis azimut blir endret, skal data korrigeres for denne spesielle endring i azimut. En annen løsning er å anta at hver ny regularisert trase har nesten den samme azimut som den nærmeste inngangstrasen. Estimering av azimutverdier ved hjelp av den nærmeste inngangstrasen er fysisk mer korrekt enn å anta en azimutverdi lik null i forhold til
seilingsretningen. Denne løsningen på azimutestimering kan imidlertid føre til problemer ved ytterligere
behandlingsalgoritmer. Det er for eksempel ønskelig for
fallutflyttingskorreksjon og forstakkmigrering å ha data som er regelmessig samplet i midtpunkt, offset og azimut. Dette er diskutert av Canning, A. og Gardner, G.H.F., 1996, "Another look at the question of azimuth:" The Leading Edge, 15, nr. 07, 821-823.
Duijindam. A.J.W. m/fl., 1999: "A general reconstruction scheme for dominant azimuth 3D seismic data", 69. Ann. Internat. Mtg: Soc. of Expl.Geophys., Expanded Abstracts, beskriver en fremgangsmåte for rekonstruksjon (regularisering) av uregelmessig samplede seismiske data under anvendelse av Fourier- eller Radon-regularisering. Deres rekonstruksjonsmåte omfatter å ompostere dataene langs mottakerlinjene til nøyaktige tverrlinjeposisjoner, fulgt av en minste kvadraters rekonstruksjon i midtpunkt-offset-domene langs tverrlinjer og etter NMO-korreksjon. De forutsetter tilfelle med seismisk datainnsamling med en fremherskede azimut for lange offset. De ignorerer imidlertid azimutvariasjon.
Duijndam, A.J.W. m/fl. påpeker at en effektiv regulariseringsmetode vil ha et antall gunstige anvendelser ved seismisk databehandling. Det kan forbedre genereringen av pseudo-nulloffsetdata for konvensjonelle bingestakkings-teknikker. Det kan regularisere og generere manglende data for forstakkbehandling som krever tett og regelmessig sampling, slik som tredimensjonal forstakkavbildning og tredimensjonal, overflaterelatert multipldempning. Det kan forbedre tilpasningen av seismiske tidsforløpsdata å forbedre amplitude/vinkel-analyse (AVA). Det kan forbedre koherent støydempning på forstakkdata, noe som muliggjør migrering med høy oppløsning av enkle, felles offset-datavolumer.
Det finnes derfor et behov for en fremgangsmåte til regularisering av uregelmessig samplede seismiske data som tilveiebringer korreksjoner for de tidsforskyvninger som skyldes azimutvariasjoner. Dette vil forbedre repeterbarheten til seismisk tidsforløpsdata-sampling og -behandling.
KORT OPPSUMMERING AV OPPFINNELSEN
Oppfinnelsen er en fremgangsmåte for å korrigere for tidsforskyvninger i seismiske data som et resultat av azimutal variasjon. Et seismisk datasett blir behandlet ved å anvende en transformasjon på det seismiske datasettet som korrigerer for tidsforskyvninger i det seismiske datasettet som et resultat av azimutal variasjon.
Alternativt blir et seismisk datasett sortert til felles offset-samlinger. De følgende trinn blir så anvendt på hver felles offset-samling: En transformasjon blir anvendt på den felles offset-samling som korrigerer for tidsforskyvninger i den felles offset-samling som er et resultat av azimutal variasjon. Den transformerte felles offset-samling blir inverstransformert.
KORT BESKRIVELSE AV TEGNINGENE
Oppfinnelsen og dens fordeler vil bli lettere forstått under henvisning til den følgende detaljerte beskrivelse og de vedføyde tegninger, hvor: Fig. 1 er et flytskjema som illustrerer behandlingstrinnene i en konvensjonell fremgangsmåte til Fourier-regularisering av seismiske data med minst to uregelmessig samlede romkoordinater og en tids-koordinat; Fig. 2 er et flytskjema som illustrerer behandlingstrinnene i en utførelsesform av fremgangsmåten ifølge oppfinnelsen for behandling av seismiske data med minst to romkoordinater og en tidskoordinat, for å korrigere for tidsforskyvninger i de seismiske data som et resultat av azimutal variasjon; Fig. 3 er et flytskjema som illustrerer behandlingstrinnene i en konvensjonell fremgangsmåte for regularisering av seismiske data med minst tre uregelmessig samplede romkoordinater og en tidskoordinat; Fig. 4 er et flytskjema som illustrerer behandlingstrinnene i en utførelsesform av fremgangsmåten ifølge oppfinnelsen for behandling av seismiske data med minst tre romkoordinater og en tidskoordinat, for å koordinere for tidsforskyvninger i de seismiske data som er et resultat av azimutal variasjon; Fig. 5 er et planriss av innsamlingsgeometrien til basisundersøkelsen og overvåkningsundersøkelsene i eksempelet; Fig. 6 er en plotting av azimutvariasjonene som funksjon av tverrlinjekoordinaten i grunnlinje- og overvåknings-undersøkelsene for den innsamlingsgeometri som er vist på fig. 5; Fig. 7a er et tverrsnitt i tverrlinjeretningen for den felles offset-seksjon på 2000 meter for grunnlinje-undersøkelsen etter regularisering uten azimut-korreksjon; Fig. 7b er et tverrsnitt i tverrlinjeretningen for den felles offset-seksjonen på 2000 meter for over-våkningsundersøkelsen etter regularisering uten az imutkorreksj on; Fig. 8a er en plotting av ankomsttidene ved maksimal amplitude som funksjon av tverrlinjekoordinaten for grunnlinje- og overvåknings-undersøkelsene, etter regularisering uten azimutkorreksjon; Fig. 8b er en plotting av tidsforskyvningene mellom de to undersøkelsene sammenlignet med azimutdifferansen mellom de to undersøkelsene, som funksjon av tverr linjekoordinaten, etter regularisering uten azimut-korreksjon; Fig. 9a er et tverrsnitt i tverrlinjeretningen for den felles offset-seksjon på 2000 meter som illustrerer differansen mellom grunnlinje- og overvåknings-undersøkelsen, og en plotting av NRMS-differansen som funksjon av tverrlinjekoordinaten, etter regularisering uten azimutkorreksjon; Fig. 9b er en plotting av den normaliserte rotmiddelkvadrat-differansen sammenlignet med azimutdifferansen mellom grunnlinje- og overvåknings-undersøkelsene som funksjon av tverrlinjekoordinaten, etter regularisering uten azimutkorreksjon; Fig. 10a er et tverrsnitt i tverrlinjeretningen for den felles offset-seksjon på 2000 meter for grunnlinje-undersøkelsen etter regularisering med azimutkorreks jon; Fig. 10b er et tverrsnitt i tverrlinjeretningen for den felles offset-seksjon på 2000 meter for overvåk-ningsundersøkelsen etter regularisering med azimutkorreks jon; Fig. Ila er en plotting av ankomsttidene ved maksimal amplitude som en funksjon av tverrslinjekoordinaten for basis- og overvåkningsundersøkelsene etter regularisering med azimutkorreksjon; Fig. 11b er en plotting av tidsforskyvningene mellom de to undersøkelsene sammenlignet med azimutdifferansene mellom de to undersøkelsene, som funksjon av tverrlinjekoordinaten, etter regularisering med azimutkorreks j on; Fig. 12a er et tverrsnitt i tverrlinjeretningen for den felles offset-seksjon på 2000 meter som illustrerer differansen mellom grunnlinje- og overvåkings- undersøkelsen, og en plotting av NRMS-differansen som funksjon av tverrlinjekoordinaten, etter regularisering med azimutkorreksjon,• Fig. 12b er en plotting av den normaliserte rotmiddelkvadrat-differansen sammenlignet med azimutdifferansen mellom grunnlinje- og overvåknings-undersøkelsene som funksjon av tverrlinjekoordinaten, etter regularisering med og uten azimutkorreksjon,• Fig. 13a er et tverrsnitt i tverrlinjeretningen som illustrerer tverrlinjene til den felles offset-seks jon på 2000 meter for grunnlinjeundersøkelsen; Fig. 13b er et tverrsnitt i tverrlinjeretningen som illustrerer differansen mellom grunnlinje- og overvåknings-undersøkelsene etter regularisering uten azimutkorreksjon; og Fig. 13c er et tverrsnitt i tverrlinjeretningen som illustrerer differansen mellom grunnlinje- og overvåkningsundersøkelsene etter regularisering med azimutkorreksj on.
Selv om oppfinnelsen vil bli beskrevet i forbindelse med dens foretrukne utførelsesformer, vil man forstå at oppfinnelsen ikke er begrenset til disse. Tvert imot er oppfinnelsen ment å dekke alle alternativer, modifikasjoner og ekvivalenter som kan innbefattes innenfor rammen av oppfinnelsen slik den er definert i de vedføyde patentkrav.
DETALJERT BESKRIVELSE AV OPPFINNELSEN
Oppfinnelsen er en fremgangsmåte for behandling av seismiske data for å korrigere for tidsforskyvninger i de seismiske data som skyldes azimutal variasjon. I spesielle utførelsesformer er oppfinnelsen en fremgangsmåte for regularisering av seismiske data ved å korrigere for tids forskyvninger i de seismiske data som er et resultat av azimutal variasjon.
Seismiske data kommer i forskjellige former avhengig av hvordan de er innsamlet. De behandlede data i et vanlig to-dimensjonalt seismisk datasett inneholder to koordinater som kan representere et bilde av en to-dimensjonal skive av jorden, men bare etter å ha blitt behandlet. Rådataene for et to-dimensjonalt, seismisk datasett blir imidlertid vanligvis innsamlet med tre koordinater. Disse omfatter to romkoordinater og en tidskoordinat. Den ene tidskoordinaten er vanligvis en toveis forplantningstid for et seismisk bølgefelt fra en seismisk kildeposisjon til en seismisk mottakerposisjon lagt ut på en en-dimensjonal linje. De to romkoordinatene er vanligvis de to tilsvarende posisjonene, hver definert ved hjelp av en koordinat, for den seismiske kilden og den seismiske mottakeren på den en-dimensjonale linjen langs hvilken dataene blir innsamlet. Denne linjen kan for eksempel være en slepekabellinje ved marine, seismiske dataunder-søkelser. Denne linjen blir så vanligvis tatt som en en-dimensjonal koordinatakse, og romkoordinatene blir så definert langs denne horisontale koordinataksen for enkelthets skyld.
Ved vanlig behandling av et to-dimensjonal, seismisk datasett blir de to romkoordinatene i rådataene omformet til en romkoordinat. Denne ene romkoordinaten er vanligvis den ene koordinaten til midtpunktet mellom kilde- og mottaker-posis jonene langs den definerende, en-dimensjonale koordinataksen. Den ene tidskoordinaten kan forbli som en tidskoordinat eller konverteres til dybde, en annen romkoordinat, men i vertikalretningen. Bare etter behandling blir derfor de seismiske rådata med tre koordinater omformet til et to-dimensjonalt, seismisk datasett med to koordinater.
Likeledes inneholder de behandlede data i et standard, tre-dimensjonalt seismisk datasett tre koordinater, men igjen bare etter behandling. Rådataene for et tre-dimensjonalt, seismisk datasett blir vanligvis innsamlet med fem koordinater som omfatter fire romkoordinater og en tidskoordinat. Den ene tidskoordinaten er igjen typisk en toveis forplantningstid for et seismisk bølgefelt fra en seismisk kildeposisjon til en seismisk mottakerposisjon, idet denne tiden er lokalisert i et to-dimensjonalt område. De fire romkoordinatene er typisk de tilsvarende to-koordinatposisjonene til den seismiske kilden og den seismiske mottakeren på et todimensjonalt overflateområde over hvilket dataene blir innsamlet. Dette overflateområde kan for eksempel være det område som effektivt undersøkes ved hjelp av en rekke slepekabler ved marine, seismiske dataundersøkelser. Dette området blir så vanligvis brukt til å definere et to-dimensjonalt system av to koordinatakser og romkoordinatene blir så definert langs disse to horisontale koordinataksene for enkelhets skyld.
Dette koordinatsystemet er typisk et kartesisk koordinatsystem hvor innsamlingen av rådataene for det seismiske datasettet lett kan beskrives. Dette koordinatsystemet vil typisk bli definert ved hjelp av linje- og tverrlinje-retningene til den seismiske undersøkelsen hvor de seismiske rådata ble innsamlet. Den første rommessige koordinataksen vil derfor typisk være i linjeretningen, og den andre rommessige koordinatkasen vil så være i tverrlinjeretningen. Denne koordinatakserepresentasjonen er imidlertid ikke en begrensning av fremgangsmåten ifølge foreliggende oppfinnelse.
Ved vanlig seismisk databehandling for et tre-dimensjonalt seismisk datasett, blir de fire romkoordinatene omformet til to romkoordinater. Disse to romkoordinatene er typisk de to koordinatene for midtpunktet til kilde- og mottaker-posisjonene. Igjen kan den ene tidskoordinaten forbli som en tidskoordinat eller omformes til dybde, en annen romkoordinat i vertikalretningen.
Alternativt kan rådataene i et tre-dimensjonalt seismisk datasett representeres av fem forskjellige koordinater. Det er igjen den ene tidskoordinaten, men fire forskjellige romkoordinater. Disse fire romkoordinatene er de to koordinatene for et midtpunkt, en offset-koordinat og en azimut-koordinat. De første og andre midtpunktkoordinatene er punktkoordinater, offset-koordinaten er en lengdekoordinat og azimut-koordinaten er en vinkelkoordinat. Dette er derfor nå en kombinasjon av et kartesisk koordinatsystem for de første og andre midtpunktkoordinatene og et polarkoordinatsystem for offset-koordinaten og azimut-koordinaten.
Generelt regulariserer en en-dimensjonal regulariseringsmetode en uregelmessig samplet romkoordinat i et seismisk datasett med minst en romkoordinat og en tidskoordinat. Det kan være ytterligere romkoordinater tilstede som ikke blir regularisert. En to-dimensjonal regulariseringsmetode regulariserer likeledes to uregelmessige samplede romkoordinater i et seismisk datasett med minst to romkoordinater og en tidskoordinat, mens en tre-dimensjonal regulariseringsmetode regulariserer tre uregelmessig samplede romkoordinater i et seismisk datasett med minst tre romkoordinater og en tidskoordinat. Fremgangsmåten ifølge oppfinnelsen vil bli beskrevet ved hjelp av utførelsesformer som innbefatter en to-dimensjonal regulariseringsmetode og en tre-dimensjonal regulariseringsmetode.
En konvensjonell fremgangsmåte for regularisering av seismiske data er Fourier-regularisering. Fourier-regularisering korrigerer imidlertid ikke for tidsforskyvninger i de seismiske data som er et resultat av azimutal variasjon. Fremgangsmåten ifølge oppfinnelsen korrigerer imidlertid for tidsforskyvninger i de seismiske data som er et resultat av azimutal variasjon. Som en hjelp til å forstå foreliggende oppfinnelse og for å skjelne foreliggende opp finnelse fra teknikkens stand, vil en konvensjonell to-dimensjonal Fourier-regulariseringsmetode først bli beskrevet under henvisning til fig. 1. Så vil en utførelsesform av foreliggende oppfinnelse bli beskrevet under henvisning til fig. 2.
Fig. 1 viser et flytskjema som illustrerer de typiske behandlingstrinnene i en konvensjonell fremgangsmåte for regularisering av seismiske data med minst to uregelmessig samplede romkoordinater og en tidskoordinat. Den følgende diskusjon av fig. 1 vil gi en kort oversikt over denne konvensjonelle fremgangsmåten. Den konvensjonelle fremgangsmåten vil bli beskrevet ved hjelp av Fourier-transformasjoner og inverse av disse. For nærmere diskusjon kan det vises til Duijndam m.fl. (1999), som er nevnt ovenfor, og Schonewille, M.A., Ph.D. avhandling, "Fourier reconstruction of irregularly sampled seismic data", Delft University of Technology, 2000.
Ved trinn 101 blir et seismisk datasett valgt for regularisering. Det seismiske datasettet blir valgt med minst to romkoordinater og en tidskoordinat. De to romkoordinatene vil bli kalt den første og den andre romkoordinaten. Den første og den andre romkooridnaten vil bli representert av henholdsvis de variable x og y.Tidskoordinaten kan representeres av forplantningstiden t. Tidskoordinaten kan alternativt omformes til dybde z ved å benytte kjennskap til eller estimater av den akustiske hastigheten i de lokale media. Her vil imidlertid tidskoordinaten være representert av tidsfrekvensen co. Denne representasjonen vil lette de ligninger som brukes til å beskrive regulariseringsmetoden. Det seismiske datasettet blir derfor definert i (x,y,co)-domenet. (x, y, co)-domenet vil her bli kalt romdomenet.
Det seismiske datasettet er vanligvis en samling av registrerte seismiske traser. Samlingen av registrerte seismiske traser i det seismiske datasettet kan her være uregelmessig samplet i de første og andre romkoordinatene. Dette betyr at posisjonen til x- og y-koordinatene til de seismiske trasene kan være uregelmessig adskilt. Fordi det seismiske datasettet kan være uregelmessig samplet, blir en spesiell foroverrettet Fourier-transformasjon beregnet for å håndtere de uregelmessig samplede dataene.
Ved trinn 102 blir en invers Fourier-transformasjon for uregelmessig samplede data utledet for det seismiske datasett som er valgt i trinn 101. Denne inverse Fourier-transformasjonen transformerer data fra (kx, ky, co)-domenet til uregelmessig samplede seismiske data i romdomenet. (kx, ky, co)-domenet blir her kalt Fourier-domenet. Her er kxog ky de bølgenumre som svarer til de første og andre romkoordinatene x og y. Utledningen av denne inverse Fourier-transformasjonen begynner med den vanlige diskrete, inverse Fourier-transformasjonen, som er velkjent på området. En utførelsesform er gitt ved:
Her er P(x,y,co) en funksjon som representerer de seismiske data i romdomenet og P(kx,ky,co) er en funksjon som representerer de transformerte, seismiske data i Fourier-domenet. De variable Akxog Aky er regelmessige samplingsintervaller i kx- og ky-koordinatretningene i Fourier-domenet. Uttrykket exp[x] er en alternativ (enklere å lese) formulering av eksponensialuttrykket ex.
Den inverse Fourier-transformasjonen som er gitt ved ligning (1) er definert for regelmessig samplede seismiske data i Fourier-domenet, men er gyldig for en hvilken som helst romposisjon definert av de første og andre romkoordinatene x og y i romdomenet. Denne inverse transformasjonen er derfor i stand til å transformere det uregelmessig samplede, seismiske datasettet som er valgt i trinn 101.
I formuleringen av den inverse Fourier-transformasjonen som er gitt i ligning (1), er videre alle kombinasjoner av m og 1 spesifisert. Dette betyr at et rektangulært område i Fourier-domenet blir brukt. Dette området kalles et rektangulært"understøttelsesområde". Denne begrensningen til et rektangulært område er ikke nødvendig. Valg av et enkelt understøttelsesområde kan imidlertid gi for mange parametre selv om det fører til en formulering med matriser som ikke kan inverteres. I en minste kvadraters invertering kan det derfor være fordelaktig å redusere antallet Fourier-koeffisienter som beregnes.
En invers Fourier-transformasjon for de uregelmessig samplede seismiske data som er valgt i trinn 101, blir derfor utledet fra den inverse Fourier-transformasjonen for regelmessig samplede data som definert i ligning (1). Dette kan utføres ved å omskrive ligning (1) som:
Her er kx,rog ky,rde bølgetall som svarer til de 2R+1 samplede punktene (kx,r/ky(r) for r= -R,...,R, i Fourier-domenet.
Det skal bemerkes at ligning (2) kan uttrykkes mer generelt som:
Potensen n kan innta en av verdiene n=l eller n=2 i ligning (3).Potensen n=l representerer enten en lineær Fourier-transformasjon eller en lineær Radon-transformasjon over offset-aksen. Potensen n=2 representerer en Fourier-transformasjon med en kvadratisk utstrekking av offset-aksen eller en parabolsk Radon-transformasjon. For den lineære Radon-trans formasjonen er kx=copx, hvor pxer langsomhets-parameteren i x-retningen. For den kvadratiske Fourier-transformasjonen eller den parabolske Radon-transformasjonen er kx=coqx, hvor ck er krumningen i x-retningen. Lignende definisjoner gjelder for ky. Se Duijndam m.fl. (1999), som diskutert ovenfor, for ytterligere detaljer. For tydelighets skyld vil imidlertid den gjenværende diskusjon av den konvensjonelle regulariseringsmetoden på fig. 1 bli illustrert ved hjelp av Fourier-transformasjoner.
Ved trinn 103 blir en foroverrettet Fourier-transformasjon beregnet fra den inverse Fourier-transformasjonen som er utledet i trinn 102. Denne foroverrettede Fourier-transformasjonen transformerer seismiske data i romdomenet til Fourier-domenet. Side den inverse Fourier-transformasjonen fra trinn 102 ble utledet fra uregelmessig samplede seismiske data i romdomenet, gjelder den foroverrettede Fourier-transformasjonen også uregelmessig samplede seismiske data i romdomenet.
Denne beregningen av den foroverrettede Fourier-transformasjonen blir vanligvis utført ved hjelp av minste kvadraters invertering. Denne minste kvadraters invertering er velkjent på området og vil bare kort bli beskrevet her.
Ligning (2) kan omskrives på matrise/vektor-form for K uregelmessig samplede samplingsposisjoner (xk,yk) for k=0,..., K-1, i romdomenet, som:
Her er vektorelementene for p gitt ved for k=0,...,K-l; idet vekt or element ene for p er gitt ved r=-R,...,R; og matriseelementet for A er gitt ved
For båndbegrensede, seismiske data innenfor under-støttelsesområdet, vil ligning (7) være nøyaktig. I praksis er imidlertid de seismiske dataene vanligvis ikke nøyaktig båndbegrenset. De seismiske datakomponentene utenfor den spesifiserte båndbredden, utgjør støyen i forovermodellen og bør inkorporeres i et støyuttrykk n. Ligning (7) blir da:
Ligning (8) representerer en forovermodell for de seismiske dataene, og er en to-dimensjonal transformasjon fra Fourier-domenet til rom domenet. For å oppnå den ønskede transformasjon av dataene fra (det uregelmessige samplede) romdomenet til Fourier-domenet, blir en minste kvadraters invertering ifølge ligning (8) brukt:
Hvor p er den minste kvadraters estimering av Fourier-spektret til de uregelmessige samplede dataene p.
Hvis understøttelsesområdet er rektangulært, så har matrisen:
en Block-Toeplitz-Toeplitz-Block-struktur (BTTB-struktur) for hvilken hurtige inverteringsmetoder er velkjente på området. Hvis understøttelsesområdet er et delsett av det rektangulære understøttelsesområdet, så kan BTTB-strukturen fremdeles brukes i en konjugert gradientmetode ved å maskere de ekstra Fourier-koeffisientene, som velkjent på området. Den ikke-uniforme, diskrete Fourier-transformasjonen kan beregnes ved å bruke en ikke-uniform hurtig Fourier-transformasjon. Den minste kvadraters formulering i ligning (9) blir utvidet til å innbefatte veie- og diagonalstabiliserings-teknikker.
Det vises igjen til fig. 1 hvor det seismiske datasett som er valgt i trinn 101, ved trinn 104 blir transformert fra romdomenet til Fourier-domenet. Denne transformasjonen blir utført ved å bruke den foroverrettede Fourier-transformasjonen som er beregnet i trinn 103. Denne transformasjonen genererer et Fourier-transformert datasett.
Ved trinn 105 blir det Fourier-transformerte datasettet fra trinn 104 inverst transformert fra Fourier-domenet tilbake til romdomenet. Denne inverse transformasjonen blir vanligvis utført ved hjelp av en invers Fourier-transformasjon. Den inverse Fourier-transformasjonen er typisk en invers, diskret Fourier-transformasjon, eller alternativt en invers, hurtig-transformasjon. Begge disse inverse Fourier-transformasjonene er velkjente på området. Denne inverse transformasjonen genererer et inverst Fourier-transformert datasett. De transformerte seismiske trasene i det inverst Fourier-transf ormerte seismiske datasettet vil nå være regelmessig samplet.
Den konvensjonelle fremgangsmåten med to-dimensjonal Fourier-regularisering som er beskrevet under henvisning til fig. 1, korrigerer ikke for tidsforskyvninger i de seismiske data som er et resultat av azimutal variasjon. En utfør-elsesform av fremgangsmåten ifølge oppfinnelsen som korrigerer for tidsforskyvninger i et seismisk datasett, som er et resultat av azimutal variasjon, vil imidlertid nå bli beskrevet under henvisning til fig. 2. Denne utførelsesformen av fremgangsmåten i henhold til oppfinnelsen kan også regularisere uregelmessig samplede seismiske data selv om regulariseringen ikke er et krav i henhold til fremgangsmåten ifølge foreliggende oppfinnelse.
Fig. 2 viser et flytskjema som illustrerer behandlingstrinnene i en utførelsesform av fremgangsmåten ifølge oppfinnelsen for behandling av seismiske data med minst to romkoordinater og en tidskoordinat, ved å korrigere for tidsforskyvninger i de seismiske data som er et resultat av azimutal variasjon.
Ved trinn 201 blir et seismisk datasett valgt for regularisering. Det seismiske datasettet blir valgt som et seismisk rådatasett med minst to romkoordinater og en tids-koordinat. De to romkoordinatene vil bli kalt den første og den andre romkoordinaten. Den første og den andre romkoordinaten vil henholdsvis bli representert av de variable x og y. De første og andre romkoordinatene vil vanligvis være orientert i linje- og tverrlinje-retningene til den seismiske undersøkelsen hvor det seismiske datasettet blir samlet inn, men dette er ikke en begrensning av oppfinnelsen. Tidskoordinaten kan representeres av forplantningstiden t eller kan alternativt omformes til dybde z ved å benytte kjennskap til eller estimater av akustisk hastighet i de lokale media.Tidskoordinaten vil imidlertid her bli representert av tidsfrekvensen co. Denne representasjonen vil forenkle de ligninger som brukes til å beskrive regulariseringsmetoden. Det seismiske datasettet er derfor definert i (x, y, co)-domenet.
(x, y, co)-domenet vil her bli kalt romdomenet.
Det seismiske datasettet er typisk en samling av registrerte seismiske traser. Samlingen av registrerte seismiske traser i det seismiske datasettet kan her være uregelmessig samplet i en eller begge av de første og andre romkoordinatene. Dette er imidlertid ikke noen begrensning av oppfinnelsen. Fremgangsmåten ifølge oppfinnelsen kan anvendes på seismiske datasett hvor en eller begge de første og andre romkoordinater er regelmessig samplet. Fremgangsmåten ifølge oppfinnelsen tilveiebringer likevel azimutal tidsforskjøvet korreksjon uansett om romkoordinatene er regelmessig eller uregelmessig samplet. Fremgangsmåten ifølge oppfinnelsen vil bli illustrert ved hjelp av tilfelle med uregelmessig samplede data, men dette er kun for illustrasjonsformål og ikke ment som en begrensning av oppfinnelsen.
Ved trinn 202 blir en invers transformasjon som korrigerer for tidsforskyvninger i seismiske data som et resultat av azimutal variasjon og som er i stand til å transformere uregelmessig samplede seismiske data, utledet for det seismiske datasettet som er valgt i trinn 201. Den inverse transformasjonen er fortrinnsvis en invers Radon-transformasjon, men dette er ikke en begrensning av oppfinnelsen. Denne inverse transformasjonen transformerer eventuelle, uregelmessige samplede seismiske data i et falldomene til regelmessig samplede seismiske data i romdomenet. I det foretrukne tilfelle med en invers Radon-transformasjon, er dette en transformasjon fra (px,py,co)-domenet til romdomenet. (px,py, co)-domenet blir ofte kalt Radon-domenet. Her er pxog py langsomhetsparameterne som representerer fall, svarende til de første og andre romkoordinatene x og y. Typisk representerer pxog py fall i henholdsvis linje- og tverrlinje-retningene, men dette er ikke en begrensning av fremgangsmåten ifølge oppfinnelsen.
Utledningen av denne inverse transformasjonen vil bli illustrert ved hjelp av et tilfelle med en invers Radon-transformasjon, men dette er kun for illustrasjonsformål og er ikke ment å være en begrensning av oppfinnelsen. En hvilken som helst invers transformasjon som transformerer seismiske data fra et falldomene til romdomenet, er ment å være innbefattet i fremgangsmåten ifølge foreliggende oppfinnelse. Utledningen av en invers Radon-transformasjon begynner med den diskrete, standard, inverse Radon-transformasjonen for regelmessig samplede seismiske data, som er velkjent på området. En utførelsesform er gitt ved:
Her er P en funksjon som representerer de seismiske dataene i romdomenet, og P er en funksjon som representerer de transformerte seismiske data i Radon-domenet. De variable Apxog Apy er regelmessig samplingsintervaller i henholdsvis px- og py-koordinatretningene i Radon-domenet. pxog py-koordinatretningene er typisk henholdsvis i linje- og tverrlinje-retningene, men dette er ikke en begrensning av fremgangsmåten ifølge oppfinnelsen.
En invers Radon-transformasjon som er i stand til å transformere uregelmessig samplede seismiske data blir deretter utledet fra den inverse Radon-transformasjonen for regelmessig samplede data som er gitt ved ligning (11). Denne inverse Radon-transformasjonen transformerer regelmessig samplede seismikkdata i Radondomenet til samplede, muligens uregelmessig, data i romdomenet. Utledningen av denne inverse Radon-transformasjonen kan utføres ved å omskrive ligning (11) til:
hvor px,rog py,rer bølgetallene som svarer til til 2R+1 samplede punktene (kx,r/ky,r) for r=-R,...,R, i Radon-domenet.
Deretter blir en invers Radon-transformasjon som korrigerer for tidsforskyvninger i seismiske data og som er et resultat av azimutal variasjon, og som er i stand til å transformere eventuelle uregelmessig samplede seismiske data, utledet fra den inverse Radon-transformasjonen for de uregelmessig samplede data gitt ved ligning (12). En tidsforskyvning At i romdomenet transformeres til en lineær fasedreining i frekvensdomenet, som svarer til en multiplikasjon med e"-^. En f oroverrettet Radon-transf ormas jon for samplede seismiske data med en tidsforskyvningskorreksjon blir derved gitt ved å omskrive ligning (12) til:
Legg merke til at tidsforskyvningen kan være en funksjon av en hvilken som helst "kjent" parameter. Som et illustrerende eksempel som ikke er begrensende for fremgangsmåten ifølge oppfinnelsen, tidsforskyvningen kan være en funksjon av de variable som beskriver traseposisjonen (for eksempel skuddkoordinater x,y, mottakerkoordinater x,y, og offset-koordinat h) og langsomhetsparameterne i begge retninger (px,py) . Med andre ord er
Deretter blir en invers Radon-transformasjon for eventuelle uregelmessig samplede seismiske data med tidsforskyvningskorreksjon for et fallende lag i en homogen undergrunn utledet fra den inverse Radon-transformasjonen for eventuelt uregelmessig samplede data med en tidsforskyvnings-korreks jon, slik som gitt ved ligning (13). Betrakt tidsforskyvninger i de to forskjellige ankomsttidene, ti og t2, for en refleksjonshendelse på grunn av azimutvariasjoner for et enkelt fallende lag i en homogen undergrunn. Chemingui, N., og Baumstein, A.I., "Handling azimuth variations in multi-streamer marine surveys", 2000, 70th Ann.Internat. Mtg: Soc. of Expl. Geophys., Expanded Abstracts, side 1-4, gir følgende formel for tidsforskyvningen i denne situasjonen:
hvor h er halve of f set-avstanden, v er hastigheten, cp er fallvinkelen, ij/i og ^ 2 er de første og andre azimut i henholdsvis de første og andre koordinatretninger, og 6 er fallretningen. og \|/2 er vanligvis azimutverdier i henholdsvis linje- og tverrlinje-retningene, men dette er ikke noen begrensning av fremgangsmåten ifølge oppfinnelsen.
Fra ligning (15) kan det observeres at de tidsforskyvninger som skyldes azimutvariasjon blir større for større offset-verdier, for lavere hastighet, for mindre ti, for brattere fall, for større azimutdifferanser (opp til1^2-^1=^/2, så blir tidsforskyvningene mindre igjen), og for større differanser mellom fallretningen 0 og den midlere undersøkelsesazimut (^1—^2) /2, opp til ji/4 (skjevt fall). For større differanser blir tidsforskyvningene igjen mindre.
Hvis differansen mellom den midlere undersøkelsesazimut og fallretningen 0 videre er null (linjefall) eller hvis den er ji/4 (tverrlinjefall) , så blir tidsforskyvningene meget små. Nøyaktige tverrlinje- og linje-fall er derfor ikke representert ved undersøkelse av effekter hvor azimut-vairasjoner er viktige. Tidsforskyvningene synes i sterk grad avhengig av offset, imidlertid skal det bemerkes at for store offset-verdier, blir azimutvariasjonen ofte mindre (spesielt for marine undersøkelser), noe som motvirker den totale effekten. Tidsforskyvningene synes også sterkt avhengige av hastighet, men legg merke til at for et lag ved en viss dybde, vil ti bli mindre hvis hastigheten øker, igjen i motsetning til den totale effekten. Tidsforskyvningen er sterkt avhengig av fallvinkelen.
Hastigheten, fallvinkelen og fallretningen er ikke kjent, og derfor må denne formelen i ligning (15) tilpasses slik at den blir avhengig av de variable hvis målte verdier er kjent, slik som xs,ys,xr,yr,h,pxog py. De følgende trinn blir fortrinnsvis brukt til å tilpasse ligning (15).Langsomhetsparameterne kan transformeres til polarkoordinater gitt av den absolutte langsomhet
og langsomhetsretningen hvor sign(py)=l for py ^ 0 og sign(py) = -1 for py < 0. For et eneste fallende lag i en homogen undergrunn, er fallretningen lik langsomhetsregningen hvis x og y relaterer til midtpunktet for en konstant of f set-seksjon. For null offset blir 2sin(cp)/v lik den absolutte langsomheten |p|, eller
Den første azimut kan derfor utledes fra skudd- og mottaker-posisjonene og den annen azimut blir valgt å være lik null i forhold til den første koordinatretningen, vanligvis linjeretningen. Alternativt kan den første azimut velges å være null i forhold til den annen koordinatretning, typisk tverrlinjeretningen. Innretting av de første og andre azimut i linje- og tverrlinje-retningene er kun for illustrasjonsformål og er ikke ment som en begrensning av oppfinnelsen. Azimutkorreksjoner kan utføres i en hvilken som helst ønsket retning.
Ved å bruke disse erstatningene som er gitt ved ligningene (16)-(18), blir formelen for tidsforskyvning i ligning (15) lik:
Den inverse Radon-transformasjonen med tidsforskyvningskorreksjon kan n gis ved å omskrive ligning (11) som:
Legg merke til at tidsforskyvningene er avhengige av ankomsttiden ti. Siden beregningene er utført med frekvenser i Radon-domenet, kan en glidende tidsvindu-løsning med overlappende tidsvinduer foretrekkes. Bruken av en løsning med glidende tidsvindu med overlappende rommessige vinduer blir likeledes foretrukket.
Det vises igjen til fig. 2 hvor en forovertransformasjon ved trinn 203 blir beregnet fra den inverse transformasjon som er utledet i trinn 202. Forovertransformasjonen er fortrinnsvis en foroverrettet Radon-transformasjon for den foretrukne inverse Radon-transformasjon som er diskutert ovenfor. Denne foroverrettede transformasjonen transformerer seismiske data i romdomenet til falldomenet. Siden den inverse transformasjonen fra trinn 202 ble utledet for å korrigere for tidsforskyvninger i de seismiske data som skyldes azimutale variasjoner, korrigerer forovertransformasjonen også for tidsforskyvninger i de seismiske data som er et resultat av azimutale variasjoner. I det tilfelle at en eller begge av de første og andre romkoordinater i det seismiske datasettet eventuelt var uregelmessig samplet, så ble den inverse transformasjon fra trinn 202 utledet for eventuelle uregelmessig samplede seismiske data i romdomenet. I dette - tilfellet gjelder forovertransformasjonen også eventuelle uregelmessig samplede seismiske data i romdomenet.
Beregningen av forovertransformasjonen vil bli illustrert ved hjelp av beregningen av en foroverrettet Radon-transf ormas jon fra den inverse Radon-transformasjon som er diskutert ovenfor, men dette er kun for illustrasjonsformål og er ikke ment som en begrensning av oppfinnelsen. Denne beregningen av den foroverrettede Radon-transformasjonen blir fortrinnsvis utført ved hjelp av et minste kvadraters invertering. Ligning (2) kan omskrives på matrise-vektor-form for K uregelmessig samplede samplingssteder (xk,yK) , for k=0,...,K-l, i romdomenet, som
Her er vektorelementene til p gitt ved for k=0,K-l; idet vektorelementene til p er gitt ved: for r=-R,...,R; og matriseelemetene i A er gitt ved
for k=0,...,K-l og r=-R,...,R.
Den lineære Radon-transformasjonen representerer de seismiske data som en sum av lineære fallhendelser. I praksis kan imidlertid ikke alle seismiske data representeres som en sum av lineære fallhendelser. I det tilfelle vil ligning (21), som definert av ligningene (22)-(24), ikke være eksakte. De datakomponenter som ikke kan representeres av den lineære
Radon-transformasjonen, bør derfor inkorporeres i et vektorstøyuttrykk n. Ligning (21) blir da
Ligning (25) representerer en forovermodell for de seismiske dataene, og er en to-dimensjonal transformasjon fra Radon-domenet til romdomenet. For å oppnå den ønskede transformasjon av dataene fra det eventuelt uregelmessig samplede romdomenet til Radon-domenet, blir det brukt en minste kvadraters invertering av ligning (25):
hvor p er det minste kvadraters estimat av Radon-spekteret for de uregelmessig samplede data p.
Den minste kvadraters formulering i ligning (26) kan utvides til å innbefatte en hvilken som helst stabiliser-ingsteknikk som er kjent på området. Dette kan innbefatte, men er ikke begrenset til, diagonale stabiliseringsteknikker for eksempel, som er velkjente på området. Anvendelse av diagonal stabilisering vil endre ligning (26) til det følgende:
hvor k er en stabiliseringskonstant og I er identitets-matrisen. Fremgangsmåter for stabilisering med høy oppløsning som for eksempel beskrevet i Sacchi, M.D. ogUlrych, T.J., 1995,"High resolution velocity gathers and offset space reconstruction", Geophysics, vol. 60, 1169-1177, kan også anvendes.
Det vises igjen til fig. 2 hvor det seismiske datasett som er valgt i trinn 201, ved 204 blir transformert fra romdomenet til falldomenet. Falldomenet er fortrinnsvis Radon-domenet. Denne transformasjonen blir fortrinnsvis utført ved å anvende den foroverrettede transformasjonen som er beregnet i trinn 203. Denne transformasjonen genererer et transformert datasett.
Ved trinn 205 blir det transformerte datasettet fra trinn 204 invers transformert fra falldomenet tilbake til romdomenet. Denne inverse transformasjonen blir fortrinnsvis gjennomført ved hjelp av en invers Radon-transformasjon. Den foretrukne inverse Radon-transformasjonen er fortrinnsvis en invers, diskret Radon-transformasjon som er velkjent på området. Denne inverse transformasjonen genererer et inverst transformert datasett. De transformerte seismiske traser i det inverst transformerte, seismiske datasettet vil nå bli korrigert for tidsforskyvninger i de seismiske data som er et resultat av azimutal variasjon. I tillegg vil nå eventuelt uregelmessig samplede seismiske data bli regelmessig samplet.
Den utførelsesform av oppfinnelsen som er beskrevet under henvisning til fig. 2, anvender et seismisk datasett med minst to romkoordinater, som eventuelt er uregelmessig samplet, og en tidskoordinat. Hvis Radon-transformasjoner brukes, kan den utførelsesform av oppfinnelsen som er beskrevet under henvisning til fig. 2, kalles to-dimensjonal Radon-regularisering med azimutal tidsforskyvningskorreksjon. Fremgangsmåten ifølge oppfinnelsen har en utførelsesform som anvender seismiske datasett med minst tre romkoordinater, som eventuelt er uregelmessig samplet, og en tidskoordinat som vil bli beskrevet under henvisning til fig. 4. Hvis Radon-transf ormas joner blir benyttet, kan den utførelsesform av oppfinnelsen som er beskrevet under henvisning til fig. 4, kalles en tre-dimensjonal Radon-regularisering med azimutal tids forskyvningskorreks j on.
Som en hjelp til å forstå foreliggende oppfinnelse og for å skjelne foreliggende oppfinnelse fra teknikkens stand, vil først en konvensjonell tre-dimensjonal Fourier-regularisering bli beskrevet under henvisning til fig. 3. Så blir en utførelsesform av foreliggende oppfinnelse beskrevet under henvisning til fig. 4.
Fig. 3 viser et flytskjema som illustrerer de typiske behandlingstrinnene i henhold til en konvensjonell fremgangsmåte for regularisering av seismiske data med minst tre uregelmessig samplede romkoordinater og en tidskoordinat.
Ved trinn 301 blir det valgt et seismisk datasett for regularisering. Det seismiske datasettet blir typisk valgt med minst tre uregelmessig samplede romkoordinater og en tidskoordinat. Det seismiske datasettet blir vanligvis valgt for å inneholde to koordinater for et midtpunkt, en offset-koordinat, og en azimut-koordinat som den fjerde romkoordinat. De første tre romkoordinatene, det vil si de to midtpunkt-koordinatene og offset-koordinatene, er de tre romkoordinatene som regulariseres. Variasjoner i den fjerde romkoordinaten, azimut-koordinaten, blir ikke kompensert for.
I trinn 302 blir det seismiske datasettet som er valgt i trinn 301, sortert til felles skudd-samlinger.
I trinn 303 blir en-dimensjonal regularisering anvendt på hver fellesskudd-samling fra trinn 302 for å regularisere den første midtpunkt-koordinaten. Vanligvis er den første midtpunkt-koordinaten 1injemidtpunkt-koordinaten. Den en-dimensjonale regularisering som anvendes, er typisk en en-dimensjonal Fourier-regularisering. Følgelig blir hver linjemidtpunkt-koordinat omposisjonert hovedsakelig i midten av bingene i linjeretningen.
I trinn 304 blir det regulariserte datasettet fra trinn 303 sortert til felles første midtpunktkoordinat-samling. Disse samlingene vil være volumer som omfatter den annen midtpunkt-koordinat, offset og tid. Det regulariserte datasettet blir vanligvis sortert til "tverrlinjer" eller felles linjemidtpunkt-koordinatsamlinger. Disse samlingene er volumer som omfatter tverrlinjemidtpunkt-koordinat, offset og tid.
I trinn 305 blir to-dimensjonal regularisering påført hver felles, første midtpunktkoordinat-samling fra trinn 304. Den to-dimensjonale regulariseringen er en typisk Fourier-regularisering. De to romdimensjonene som regulariseres, er den annen midtpunkt-koordinaten og offset. Vanligvis er den andre midtpunkt-koordinaten tverrlinjemidtpunkt-koordinaten, slik at de to romdimensjonene som regulariseres, er tverrlinjemidtpunkt-koordinaten og offset.
Den konvensjonelle fremgangsmåten for tre-dimensjonal Fourier-regularisering som er beskrevet under henvisning til fig. 3, korrigerer ikke for tidsforskyvninger i de seismiske data som et resultat av azimutal variasjon. En utførelsesform av fremgangsmåten ifølge oppfinnelsen som korrigerer for tidsforskyvninger i et seismisk datasett som et resultat av azimutal variasjon, vil imidlertid nå bli beskrevet under henvisning til fig. 4. Denne utførelsesformen av fremgangsmåten ifølge oppfinnelsen kan også regularisere uregelmessig samplede seismiske data, selv om regularisering ikke er noe krav i henhold til fremgangsmåten ifølge oppfinnelsen.
Fig. 4 viser et flytskjema som illustrerer behandlingstrinnene ifølge en utførelsesform av fremgangsmåten i henhold til oppfinnelsen for behandling av seismiske data med minst tre romkoordinater og en tidskoordinat, som korrigerer for tidsforskyvninger i de seismiske data, som er et resultat av azimutal variasjon.
Ved trinn 401 blir det valgt et seismisk datasett for behandling. Det seismiske datasettet er valgt som et seismisk rådatasett med minst tre romkoordinater og en tidskoordinat.
I trinn 402 blir det seismiske datasettet som er valgt i trinn 401, omformet til en datarepresentasjon hvor de fire romkoordinatene vil være to midtpunkt-koordinater, offset og azimut. De første tre romkoordinatene, det vil si de første og andre midtpunkt-koordinatene og offset-koordinaten, vil være de tre romkoordinatene som kan regulariseres hvis en av de første tre romkoordinatene er uregelmessig samplet.Regularisering er ikke et krav i henhold til oppfinnelsen.Tidsforskyvning som et resultat av variasjon i den fjerde romkoordinaten, azimut-koordinaten, vil bli korrigert ifølge denne fremgangsmåten i henhold til oppfinnelsen. De første og andre midtpunkt-koordinatene vil typisk være henholdsvis linje- og tverrlinje-midtpunktkoordinatene, og denne konvensjonen kan anvendes her for i1lustrasjonsformål. Dette er imidlertid ikke en begrensning av oppfinnelsen. Spesielt kan koordinatretningene defineres slik at azimut-korreksjonene er i en hvilket som helst ønsket retning.
I trinn 403 blir det seismiske datasett som er omformet i trinn 402, sortert til felles skudd-samlinger.
I trinn 404 blir en-dimensjonale regulariseringen fortrinnsvis anvendt på hver felles skudd-samling fra trinn 403 hvis regularisering ar ønsket. Den eventuelle en-dimensjonale regularisering som anvendes, er fortrinnsvis en en-dimensjonal Fourier-regularisering. Den eventuelle en-dimensjonale regularisering regulariserer fortrinnsvis offset-koordinaten.
Det to-dimensjonale regulariseringstrinn, trinn 406 nedenfor, i denne tre-dimensjonale utførelsesformen av fremgangsmåten i henhold til oppfinnelsen, vil fortrinnsvis bli utført på et datasett som omfatter de første og andre midtpunkt-koordinatene og tidskoordinaten. Fallparameterne i et falldomene, slik som spesielt langsomhetsparameterne pxog py i Radon-domenet, bør svare til fallene i lagene i undergrunnen som representeres av det seismiske datasett. Det en-dimensjonale regulariseringstrinnet, hvis det utføres, blir fortrinnsvis utført på offset-koordinaten.
I trinn 405 blir datasettet fra trinn 404, eventuelt regularisert, sortert til felles offset-samlinger. Disse samlingene vil være volumer som omfatter første og andre midtpunkt-koordinater, azimut og tid. De første og andre midtpunkt-koordinatene vil typisk være henholdsvis linje- og tverr1inje-midtpunktkoordinatene. Dette er imidlertid kun for illustrasjonsformål og ikke en begrensning av oppfinnelsen.
I trinn 406 blir den to-dimensjonale utførelsesform av fremgangsmåten ifølge oppfinnelsen, som beskrevet under henvisning til fig. 2 ovenfor, anvendt på hver felles offset-samling fra trinn 405. Den to-dimensjonale utførelsesformen av fremgangsmåten ifølge oppfinnelsen er fortrinnsvis en to-dimensjonal Radon-regularisering med korreksjon av azimut-tidsforskyvning. Hver felles offset-samling blir korrigert for tidsforskyvninger i de seismiske data som er et resultat av azimutal variasjon. Azimutverdien blir følgelig hovedsakelig konvertert til null azimut. Azimutverdien er nå hovedsakelig reorientert inn i den første midtpunkt-koordinatretningen, typisk linjeretningen.
Alternativt kan azimutverdier hovedsakelig reorienteres til den andre midtpunkt-koordinatretningen, vanligvis tverrlinjeretningen. Innretting av azimutverdiene i linje-eller tverrlinje-retningene er kun for illustrasjonsformål og er ikke ment som en begrensning av oppfinnelsen. Azimut korreksjoner kan utføres i en hvilken som helst ønsket retning.
De to romdimensjonene som kan regulariseres om ønsket, er de første og andre midtpunkt-koordinatene, vanligvis linje- og tverr1inj e-midtpunktkoordinatene.
Den utførelsesform av fremgangsmåten ifølge oppfinnelsen som er beskrevet under henvisning til fig. 4, korrigerer seismiske data med minst tre romkoordinater og en tids-koordinat for tidsforskyvninger i det seismiske datasettet som er et resultat av azimutal variasjon. Denne utførelsesformen av fremgangsmåten i henhold til oppfinnelsen kan også regularisere uregelmessig samplede seismiske data i de minst tre romkoordinatene, selv om regularisering ikke er en nødvendighet ifølge fremgangsmåte i henhold til oppfinnelsen.
Eksempler
Fremgangsmåten ifølge oppfinnelsen vil for eksempel bli anvendt på en modell med et eneste fallplan i en homogen undergrunn. Dataene representerer en tidsforløpsundersøkelse. De parametre som beskriver innsamlingsgeometrien er bistabil (flip-flop) avfyring med et 25 m skuddintervall (flip til flop), 16 slepekabler, avstand på 75 m mellom slepekablene og ingen vinkelinnstilling. To undersøkelser er blitt modellert som representerer en tidsforløpsmåling. Grunnlinje- og overvåknings-undersøkelsen er skutt i samme retning, men seilretningene er ikke repetert. Tverrlinjeavstanden mellom de to undersøkelsene er 100 m. Fig. 5 viser et planriss av innsamlingsgeometrien til grunnlinjeundersøkelsen 501 og overvåkningsundersøkelsen 502. Den sorte boksen 503 viser den delen som blir brukt til regulariseringen. En felles offset-seksjon blir brukt. Fig. 6 vier en plotting av azimutvariasjonene som funksjon av tverrlinje-koordinaten i tilfelle med overvåk- ningsundersøkelser for den innsamlingsgeometrien som er vist på fig. 5. Fig. 6 viser azimutvariasjonen for grunnlinje-undersøkelsen 601 for en tverrlinje, azimutvariasjonen for overvåkningsundersøkelsen 602 og differansen mellom de to azimutvariasjonene 603.
I eksemplet blir fremgangsmåten ifølge oppfinnelsen anvendt på en homogen undergrunn med et eneste fallende lag. Et fallende lag med 7 graders fall og en 45 graders fall-retning blir brukt her. En felles offset-seksjon på 2000 m blir valgt. Hastigheten er 1480 m/s, og dybden er omkring 1000 m. Fig. 7a og 7b, 8a og 8b, og 9a og 9b viser resultatene av regularisering uten korreksjon for azimut-tidsforskyvning. Fig. 10a og 10b, lia og 11b og 12a og 12b viser de tilsvarende resultater av regularisering med korreksjon for azimut-tidsforskyvning. Fig. 7a og 7b viser tverrsnitt i tverrlinje-retningen for felles offset-seksjonen på 2000 m etter konvensjonell regularisering uten korreksjon for azimut-tidsforskyvning. Fig. 7a viser tverrlinjene for grunnlinjeundersøkelsen og fig. 7b viser tverrlinjene for overvåkningsundersøkelsen. Tidene for maksimal amplitude blir plukket ut etter delsampling, slik at nøyaktigheten er bedre enn samplingsintervallet på 4 ms. Linjene 701 på fig. 7a og 702 på fig. 7b viser tiden for den maksimale amplituden. Fig. 8a og 8b viser flere resultater av regularisering uten korreksjon for azimut-tidsforskyvning. Fig. 8a viser en plotting av ankomsttidene ved maksimal amplitude som funksjon av tverrlinje-koordinaten for grunnlinjen 801 og overvåknings-undersøkelsene 802. Fig. 8b viser en plotting av tidsforskyvningene 803 mellom de to undersøkelsene sammenlignet med azimut-differansen 804 mellom de to undersøkelsene, som funksjon av tverrlinje-koordinaten. Som ventet korrelerer tidsforskyvningene mellom grunnlinje- og overvåknings- undersøkelsen meget sterkt med azimut-differansene. Fig. 8b viser at bare små tidsforskyvninger finnes mellom de to undersøkelsene. Den maksimale tidsforskyvning mellom undersøkelsene er omkring 2,5 ms, mindre enn det vanlige tidssamplingsmellomrommet på 4 ms. Fig. 9a og 9b viser flere resultater av regularisering uten azimutkorreksjon. Fig. 9a viser et tverrsnitt av den felles offset-seksjon på 2000 m som illustrerer differansen 901 mellom grunnlinje- og overvåkningsundersøkelsen i tverrslinje-retningen, og en plotting av den normaliserte rotmiddelkvadrat-differansen (NRMS-differansen) 902 som funksjon av tverrlinje-koordinaten. Fig. 9b viser en plotting av NRMS-differansen 903 sammenlignet med azimut-differansen 904 mellom grunnlinje- og overvåknings-undersøkelsene som funksjon av tverrlinje-koordinaten. Fig. 9a viser at hvis grunnlinje- og overvåknings-undersøkelsen blir subtrahert, så betyr disse tidsforskyvningene at differansen ikke er lik null. Fig. 9b viser igjen at differansen mellom undersøkelsen korrelerer sterkt med azimut-differansen, og dermed tidsforskyvningene. Den totale NRMS-differansen uten azimutkorreks jon er 25%, og den maksimale NRMS-differansen er 40%. Fig. 10a og 10b viser tverrsnitt av tverrlinjer som funksjon av tverrlinje-koordinaten for den 2000 m felles offset-seksjonen etter regularisering med azimutkorreksjon. Fig. 10a viser tverrlinjene for grunnlinjeundersøkelsen, og fig. 10b viser tverrlinjene for overvåkningsundersøkelsen. Linjene 1001 på fig. 10a og 1002 på fig. 10b viser tids-punktene for den maksimale amplitude. Fig. 10a og 10b kan sammenlignes med fig. 7a og 7b. Fig. Ila og 11b viser flere resultater av regularisering med azimut-korreksjon. Fig. Ila viser en plotting av ankomsttidene ved maksimal amplitude som funksjon av tverrlinje-koordinaten for grunnlegningen 1101 og overvåkningsunder- søkelsene 1102. Fig. 11b viser en plotting av tidsforskyvningene 1103 mellom de to undersøkelsene sammenlignet med azimut-differansen 1104 mellom de to undersøkelsene som funksjon av tverrlinje-koordinaten. Fig. Ila viser at ankomsttidene ved maksimal amplitude for grunnlinje- og overvåknings-undersøkelsen nå er nesten lik. Fig. 11b viser at tidsforskyvningene er i størrelsesorden noen titalls millisekunder og nå ikke oppviser korrelasjon med azimut-dif f eransene . Fig. Ila og 11b kan sammenlignes med fig. 8a og 8b.
Fig. 12a og 12b viser flere resultater av regularisering med azimut-korreksjon. Fig. 12a viser et tverrsnitt i tverrlinjeretningen til felles offset-seksjonen på 2000 m, som illustrerer differansen 1201 mellom grunnlinje- og overvåknings-undersøkelsen, og en plotting av NRMS-differansen 1202 som funksjon av tverrlinje-koordinaten. Sammenligning av fig. 12a med fig. 9a viser at azimut-korreksjonen fører til meget bedre repeterbarhet. Fig. 12b viser en plotting av NRMS-dif f eransen 1203 mellom grunnlinje- og
overvåkningsundersøkelsene sammenlignet med azimut-differansen 1204 mellom grunnlinje og overvåknings-undersøkelsene, som funksjon av tverrlinje-koordinaten. Til sammenligning erNRMS-dif f eransen 1205 mellom grunnlinje og overvåknings-undersøkelsene uten azimut-korreksjon fra fig. 6b, vist.
Til slutt viser fig. 13a, 13b og 13c tverrsnitt i tverrlinjeretningen, alle i samme skala for sammenligningens skyld. Fig. 13a viser tverrlinjene for den felles offset-seksjon på 2000 m i grunnlinjeundersøkelsen. Fig. 13b viser differansen mellom grunnlinje- og overvåkningsundersøkelsene etter regularisering uten azimutkorreksjon. Fig. 13c viser differansen mellom grunnlinje- og overvåknings-undersøkelsene etter regularisering med azimut-korreksjon- Den totale NRMS-dif f eransen med azimut-korreksjon er 3,1% og den maksimale NRMS-differansen er 6,0%. Uten azimut-korreksjon var disse tallene 25% og 40%.
Oppfinnelsen er en fremgangsmåte for å korrigere for tidsforskyvninger i det seismiske datasettet som er et resultat av azimutal variasjon. Oppfinnelsen kan i tillegg brukes til regularisering av seismiske datasett. Korreksjonen for tidsforskyvningene er blitt eksplisitt beskrevet i både to-dimensjonale og tre-dimensjonale Radon-regulariserings-utførelsesformer.
Korreksjon av tidsforskyvninger er illustrert med tilfelle av et enkelt fallende lag i en homogen undergrunn, men oppfinnelsen er ikke begrenset av dette. Azimutkorreks jonen er en fall- og azimut-avhengig tidsforskyvning som illustrert i forovermodellen i henhold til den to-dimensjonale, minste kvadraters Radon-transformasjon.
Det skal bemerkes at det foregående kun er en detaljert beskrivelse av spesielle utførelsesformer av oppfinnelsen, og at mange endringer, modifikasjoner og alternativer i forhold til de beskrevne utførelsesformer kan foretas i samsvar med den gitte beskrivelse uten å avvike fra oppfinnelsens ramme. Den foregående beskrivelse er derfor ikke ment å begrense omfanget av oppfinnelsen. Omfanget av oppfinnelsen skal i stedet kun bestemmes av de vedføyde patentkrav og deres ekvivalenter.

Claims (36)

1. Fremgangsmåte for behandling av et seismisk datasett,karakterisert ved: å anvende en transformasjon på det seismiske datasettet som korrigerer for tidsforskyvninger i det seismiske datasett som er et resultat av azimutal variasjon.
2. Fremgangsmåte ifølge krav 1, hvor transformasjonen er beregnet fra en invers transformasjon som korrigerer for tidsforskyvninger i det seismiske datasett som er et resultat av azimutal variasjon.
3. Fremgangsmåte ifølge krav 2, videre omfattende: å beregne en forovertransformasjon fra den inverse trans formas j onen.
4. Fremgangsmåte ifølge krav 3, videre omfattende: å omforme det seismiske datasett ved å benytte den beregnede forovertransformasjonen; og å inverstransformere det transformerte, seismiske datasett.
5. Fremgangsmåte ifølge krav 2, hvor inverstransformasjonen er definert fra et falldomene til et romdomene.
6. Fremgangsmåte ifølge krav 5, hvor falldomenet er et Radon-domene.
7. Fremgangsmåte ifølge krav 2, hvor den inverse transformasjon er en invers Radon-transformasjon.
8. Fremgangsmåte ifølge krav 7, hvor den inverse Radon- transformasjonen er definert ved:
hvor P er en funksjon som representerer de seismiske data i romdomenet, P er en funksjon som representerer de transformerte seismiske data i Radon-domenet, x og y er romkoordinater, co er tidsfrekvens, pxog py er langsomhetsparametre som svarer til henholdsvis x og y, henholdsvis Apxog Apy er regulære samplingsintervaller i px- og py-koordinatretningene, h er en halv offset, ti er en første ankomsttid for en refleksjonshendelse, ijfier en første azimut, \ p\ ~ yjpl + Py 2 er en absolutt langsomhet og pe=sign (py) cos"<1>(Px/|p|) er en langsomhetsretning.
9. Fremgangsmåte ifølge krav 1, videre omfattende et innledende trinn som omfatter: å omforme det seismiske datasettet til en representasjon med en azimut-koordinat, minst en annen romkoordinat og en tidskoordinat.
10. Fremgangsmåte ifølge krav 9, hvor minst en av den minst ene andre romkoordinat er uregelmessig samplet.
11. Fremgangsmåte ifølge krav 10, hvor den transformasjonen som innbefatter en korreksjon for tidsforskyvninger i det seismiske datasettet som et resultat av azimutal variasjon, blir benyttet til å regularisere minst en av den minst ene av de andre romkoordinatene.
12. Fremgangsmåte ifølge krav 4, hvor trinnet med invers transformering omfatter: å anvende en invers Radon-transformasjon på det transformerte datasettet.
13. Fremgangsmåte ifølge krav 4, hvor trinnet med å beregne en foroverrettet transformasjon videre omfatter: å anvende en minste kvadraters invertering på den inverse transformasjon.
14. Fremgangsmåte ifølge krav 13, hvor trinnet med å anvende en minste kvadraters invertering videre omfatter. å anvende en stabiliseringsmetode på den minste kvadraters invertering.
15. Fremgangsmåte ifølge krav 14, hvor stabiliseringsmetoden er en diagonal stabiliseringsmetode.
16. Fremgangsmåte ifølge krav 14, hvor stabiliseringsmetoden er en stabilisering med høy oppløsning.
17. Fremgangsmåte ifølge krav 1, videre omfattende: å sortere det seismiske datasettet til felles offset-samlinger; å anvende transformasjonen på den felles offset-samler og inverstransformere den transformerte, felles offset-samling.
18. Fremgangsmåte ifølge krav 17, hvor trinnet med å anvende transformasjonen på den felles offset-samling videre omfatter følgende innledende trinn: å utlede en invers transformasjon som innbefatter en korreksjon for tidsforskyvninger i den felles offset-samling som er et resultat av azimutal variasjon; og å beregne den anvendte transformasjonen ut fra den utledede, inverse transformasjon.
19.Fremgangsmåte ifølge krav 18, hvor den utledede inverse transformasjon blir definert fra et falldomene til et romdomene.
20. Fremgangsmåte ifølge krav 19, hvor falldomenet er et Radon-domene.
21. Fremgangsmåte ifølge krav 18, hvor den utledede, inverse transformasjon er en invers Radon-transformasjon.
22. Fremgangsmåte ifølge krav 21, hvor den inverse Radon-transf ormas jon er definert ved:
hvor P er en funksjon som representerer de seismiske data i romdomenet, P er en funksjon som representerer de transformerte, seismiske data i Radon-domenet, x og y er romkoordinater, co er tidsfrekvens, pxog pY er langsomhetsparametre som svarer til henholdsvis x og y, Apxog Apy er regulære samplingsintervaller i px- og py- koordinatretningene, h er halve offset-verdien, ti er en første ankomsttid for en refleksjonshendeIse, tyier en første azimut, \ p\ = yjp*- + Pl er en absolutt langsomhet og pe=sign (py) cos"<1>(Px/|p|) er en langsomhetsretning.
23. Fremgangsmåte ifølge krav 17, videre omfattende et innledende trinn med: å omforme det seismiske datasettet til en representasjon med en første midtpunkt-koordinat, en annen midtpunkt-koordinat, en offset-koordinat, en azimut-koordinat og en tids-koordinat.
24.Fremgangsmåte ifølge krav 23, hvor den første midtpunkt-koordinat er en linjemidtpunkt-koordinat og en annen midtpunkt-koordinat er en tverrlinjemidtpunkt-koordinat.
25.Fremgangsmåte ifølge krav 23, hvor minst en av den første midtpunkt-, annen midtpunkt- og offset-koordinat er uregelmessig samplet.
26. Fremgangsmåte ifølge krav 17, videre omfattende de innledende trinn: å sortere det seismiske datasett til felles skudd-samlinger; å anvende en-dimensjonal regularisering på hver felles skudd-samling; å sortere de regulariserte felles skudd-samlingene til felles offset-samlinger.
27.Fremgangsmåte ifølge krav 26, hvor den en-dimensjonale regularisering er en en-dimensjonal Fourier-regularisering.
28. Fremgangsmåte ifølge krav 26, hvor den en-dimensjonale regularisering blir anvendt for å regularisere offset-verdier.
29. Fremgangsmåte ifølge krav 25, hvor transformasjonen som innbefatter en korreksjon for tidsforskyvninger i det seismiske datasett på grunn av azimutal variasjon, blir utledet for å regularisere minst en av de første og andre midtpunkt-koordinater.
30. Fremgangsmåte ifølge krav 18, hvor trinnet med å beregne den anvendte transformasjon videre omfatter: å anvende en minste kvadraters invertering på den utledede inverse Radon-transformasjon.
31. Fremgangsmåte ifølge krav 30, hvor trinnet med å påføre en minste kvadraters invertering videre omfatter: å anvende en stabiliseringsmetode på den minste kvadraters invertering.
32. Fremgangsmåte ifølge krav 31, hvor stabiliseringsmetoden er en diagonal stabiliseringsmetode.
33. Fremgangsmåte ifølge krav 31, hvor stabiliseringsmetoden er en stabiliseringsmetode med høy oppløsning.
34. Fremgangsmåte ifølge krav 17, hvor trinnet med invers transformering videre omfatter: å anvende en invers Radon-transformasjon på det transformerte datasettet.
35. Fremgangsmåte ifølge krav 18, hvor trinnet med å utlede en invers transformasjon benytter en glidende tidsvindu-løsning med overlappende tidsvinduer.
36. Fremgangsmåte ifølge krav 18, hvor trinnet med å utlede en invers transformasjon benytter en glidende romvindu-løsning med overlappende rommessige vinduer.
NO20041809A 2003-05-30 2004-05-04 Fremgangsmate for prosessering av et seismisk datasett NO330675B1 (no)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US10/448,809 US6889142B2 (en) 2003-05-30 2003-05-30 Method of correcting for time shifts in seismic data resulting from azimuthal variation

Publications (2)

Publication Number Publication Date
NO20041809L NO20041809L (no) 2004-11-30
NO330675B1 true NO330675B1 (no) 2011-06-06

Family

ID=32682491

Family Applications (1)

Application Number Title Priority Date Filing Date
NO20041809A NO330675B1 (no) 2003-05-30 2004-05-04 Fremgangsmate for prosessering av et seismisk datasett

Country Status (5)

Country Link
US (1) US6889142B2 (no)
AU (1) AU2004202317B2 (no)
BR (1) BRPI0401885B1 (no)
GB (1) GB2402217B (no)
NO (1) NO330675B1 (no)

Families Citing this family (47)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7626887B2 (en) * 2006-04-19 2009-12-01 Westerngeco L.L.C. Displacement field calculation
US7830747B2 (en) * 2006-11-01 2010-11-09 Pgs Geophysical As Method for multi-azimuth prestack time migration for general heterogeneous, anisotropic media
US7969818B2 (en) * 2007-12-19 2011-06-28 Pgs Geophysical As Method for regularizing offset distribution in towed seismic streamer data
US8139440B2 (en) * 2007-12-27 2012-03-20 Westerngeco L.L.C. Spectral conditioning for surface seismic data
EP2376944A4 (en) * 2008-12-17 2017-02-22 Exxonmobil Upstream Research Company Method for imaging of targeted reflectors
US8724429B2 (en) 2008-12-17 2014-05-13 Exxonmobil Upstream Research Company System and method for performing time-lapse monitor surverying using sparse monitor data
AU2009333603B2 (en) 2008-12-17 2014-07-24 Exxonmobil Upstream Research Company System and method for reconstruction of time-lapse data
US8199611B2 (en) * 2009-02-05 2012-06-12 Westerngeco L.L.C. Deriving tilt-corrected seismic data in a multi-axis seismic sensor module
US9536022B1 (en) 2009-06-01 2017-01-03 Paradigm Sciences Ltd. Systems and methods for modeling faults in the subsurface
US8600708B1 (en) 2009-06-01 2013-12-03 Paradigm Sciences Ltd. Systems and processes for building multiple equiprobable coherent geometrical models of the subsurface
US8711140B1 (en) 2009-06-01 2014-04-29 Paradigm Sciences Ltd. Systems and methods for building axes, co-axes and paleo-geographic coordinates related to a stratified geological volume
US9418182B2 (en) 2009-06-01 2016-08-16 Paradigm Sciences Ltd. Systems and methods for building axes, co-axes and paleo-geographic coordinates related to a stratified geological volume
US8743115B1 (en) 2009-10-23 2014-06-03 Paradigm Sciences Ltd. Systems and methods for coordinated editing of seismic data in dual model
US9013956B2 (en) * 2009-10-27 2015-04-21 Chevron U.S.A Inc. Method and system for seismic imaging and earth modeling using beam tomography
US8126652B2 (en) * 2010-01-20 2012-02-28 Pgs Geophysical As Azimuth correction for data reconstruction in three-dimensional surface-related multiple prediction
US8478531B2 (en) * 2010-01-20 2013-07-02 Pgs Geophysical As Dip-based corrections for data reconstruction in three-dimensional surface-related multiple prediction
US10054703B2 (en) * 2010-03-23 2018-08-21 Pgs Geophysical As Method of imaging the subsurface using stacked seismic data from azimuthally varying velocity and amplitude information
US8301385B2 (en) * 2010-06-29 2012-10-30 Chevron U.S.A. Inc. Shot gather data beamer and debeamer
US20120215453A1 (en) * 2011-02-22 2012-08-23 Cggveritas Services Sa Device and method for multi-dimensional coherency driven denoising data
US9075162B2 (en) 2011-11-10 2015-07-07 Pgs Geophysical As Method and system for separating seismic sources in marine simultaneous shooting acquisition
US9217803B2 (en) 2012-01-12 2015-12-22 Cggveritas Services Sa Device and method for estimating time-shifts
US10114134B2 (en) 2012-03-02 2018-10-30 Emerson Paradigm Holding Llc Systems and methods for generating a geological model honoring horizons and faults
US9759826B2 (en) 2012-04-03 2017-09-12 Paradigm Sciences Ltd. System and method for generating an implicit model of geological horizons
US20140088879A1 (en) * 2012-09-27 2014-03-27 Chevron U.S.A. Inc. System and method for noise attenuation in seismic data
US9477010B2 (en) 2013-03-15 2016-10-25 Paradigm Sciences Ltd. Systems and methods to build sedimentary attributes
US9575194B2 (en) 2013-05-01 2017-02-21 Cgg Services Sas Method apparatus and system for migration noise attenuation and image enhancement
US9651695B2 (en) 2013-09-19 2017-05-16 Pgs Geophysical As Construction and application of angle gathers from three-dimensional imaging of multiples wavefields
US10795053B2 (en) 2013-10-29 2020-10-06 Emerson Paradigm Holding Llc Systems and methods of multi-scale meshing for geologic time modeling
CN104749620A (zh) * 2013-12-26 2015-07-01 中国石油化工股份有限公司 一种方位角抽取道集一致性处理方法
US20150198729A1 (en) * 2014-01-13 2015-07-16 Cgg Services Sa Regularization of spatially aliased seismic data
US10670757B2 (en) * 2014-02-26 2020-06-02 Pgs Geophysical As Methods and systems for quantifying coherency and constraining coherency-based separation in simultaneous shooting acquisition
US10422923B2 (en) 2014-03-28 2019-09-24 Emerson Paradigm Holding Llc Systems and methods for modeling fracture networks in reservoir volumes from microseismic events
US9690002B2 (en) 2015-06-18 2017-06-27 Paradigm Sciences Ltd. Device, system and method for geological-time refinement
CN105929447B (zh) * 2016-05-31 2017-11-21 中国石油大学(华东) 考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法
US10466388B2 (en) 2016-09-07 2019-11-05 Emerson Paradigm Holding Llc System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles
US10732311B2 (en) 2016-10-19 2020-08-04 Pgs Geophysical As Method for adjusting attributes of marine seismic survey data
US11016212B2 (en) 2017-04-11 2021-05-25 Saudi Arabian Oil Company Compressing seismic wavefields in three-dimensional reverse time migration
CN108008447B (zh) * 2017-10-23 2019-09-10 中国石油天然气集团公司 一种炮集记录处理方法、装置及计算机存储介质
US10684382B2 (en) * 2018-01-23 2020-06-16 Saudi Arabian Oil Company Generating target-oriented acquisition-imprint-free prestack angle gathers using common focus point operators
US10969466B2 (en) * 2018-11-13 2021-04-06 GM Global Technology Operations LLC Sliding window integration scheme for object detection in a radar system
US10520644B1 (en) 2019-01-10 2019-12-31 Emerson Paradigm Holding Llc Imaging a subsurface geological model at a past intermediate restoration time
US11047999B2 (en) 2019-01-10 2021-06-29 Chevron U.S.A. Inc. System and method for seismic imaging
US11156744B2 (en) 2019-01-10 2021-10-26 Emerson Paradigm Holding Llc Imaging a subsurface geological model at a past intermediate restoration time
CN111722274B (zh) * 2019-03-21 2023-02-07 中国石油天然气集团有限公司 拾取地震波时距曲线的方法及装置
US11353609B2 (en) 2019-12-20 2022-06-07 Saudi Arabian Oil Company Identifying geologic features in a subterranean formation using angle domain gathers sampled in a spiral coordinate space
US11656378B2 (en) 2020-06-08 2023-05-23 Saudi Arabian Oil Company Seismic imaging by visco-acoustic reverse time migration
CN116626750A (zh) * 2023-03-23 2023-08-22 华北科技学院(中国煤矿安全技术培训中心) 一种基于Radon变换的矿山微震层析成像定位算法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4797861A (en) 1985-11-18 1989-01-10 Western Atlas International, Inc. Method of processing seismic data
US5719822A (en) * 1996-10-04 1998-02-17 Vector Seismic Data Processing, Inc. Seismic data radon dip moveout method
US6272435B1 (en) * 1999-09-23 2001-08-07 Pgs Tensor, Inc. Migration with dipping acquisition plane system and method
US6292754B1 (en) * 1999-11-11 2001-09-18 Bp Corporation North America Inc. Vector recomposition of seismic 3-D converted-wave data

Also Published As

Publication number Publication date
AU2004202317A1 (en) 2004-12-16
GB2402217A (en) 2004-12-01
GB0411899D0 (en) 2004-06-30
AU2004202317B2 (en) 2008-10-30
US20040243312A1 (en) 2004-12-02
BRPI0401885B1 (pt) 2017-06-06
US6889142B2 (en) 2005-05-03
NO20041809L (no) 2004-11-30
GB2402217B (en) 2006-11-15
BRPI0401885A (pt) 2005-07-26

Similar Documents

Publication Publication Date Title
NO330675B1 (no) Fremgangsmate for prosessering av et seismisk datasett
AU2013213704B2 (en) Device and method for directional designature of seismic data
EP2350693B1 (en) Continuous adaptive surface wave analysis for three-dimensional seismic data
CN108139498B (zh) 具有振幅保持的fwi模型域角度叠加
US9250341B2 (en) Device and method for extrapolating specular energy of reverse time migration three dimensional angle gathers
NO20121031A1 (no) Prosess for a karakterisere utviklingen av er reservoar
NO320103B1 (no) Seismisk prosessering med generelle ikke-hyperbolske gangtidskorreksjoner
EP3094993B1 (en) Device and method for deghosting seismic data using sparse tau-p inversion
NO20121524A1 (no) Innretning og fremgangsmate for a fjerne spokelser fra data registrert av streamer med variabel dybde
Van Dalen et al. Retrieving surface waves from ambient seismic noise using seismic interferometry by multidimensional deconvolution
EP3529640A1 (en) Generating subterranean imaging data based on vertical seismic profile data and ocean bottom sensor data
Perez et al. Improving lateral and vertical resolution of seismic images by correcting for wavelet stretch in common-angle migration
Biondi et al. Transformation of 3-D prestack data by azimuth moveout (AMO)
EP3052967B1 (en) System and method for seismic adaptive imaging
Lehujeur et al. Eikonal Tomography Using Coherent Surface Waves Extracted From Ambient Noise by Iterative Matched Filtering—Application to the Large‐N Maupasacq Array
CN108549102B (zh) 联合梯度结构张量和多窗分析的地层结构曲率估计方法
CN108957545A (zh) 气枪阵列子波方向性反褶积方法及系统
CN115598704A (zh) 一种基于最小二乘逆时偏移生成保幅角道集的方法、设备及可读存储介质
NO20130064A1 (no) Innretning og fremgangsmate for a estimere tidsforskyvning.
CHENG et al. Prestack Depth migration with the Finite‐Difference Method in Frequencyspace Domain
US9014985B2 (en) System and method for compensating time and offset varying near-surface effects in seismic data background
Fam et al. Multi-Focusing stacking using the Very Fast Simulated Annealing global optimization algorithm
Palmer Can full waveform refraction imaging overcome the century-long impasse at the first breaks?
Jang et al. Research Article Application of Reverse Time Migration to Faults Imaging in Rakhine Basin, Myanmar
Van Leeuwen The reconstruction of missing marine seismic data