NO339499B1 - Generalisert 3D-overflate-multippel prediksjon - Google Patents

Generalisert 3D-overflate-multippel prediksjon Download PDF

Info

Publication number
NO339499B1
NO339499B1 NO20065127A NO20065127A NO339499B1 NO 339499 B1 NO339499 B1 NO 339499B1 NO 20065127 A NO20065127 A NO 20065127A NO 20065127 A NO20065127 A NO 20065127A NO 339499 B1 NO339499 B1 NO 339499B1
Authority
NO
Norway
Prior art keywords
registered
deviation
target
traces
pair
Prior art date
Application number
NO20065127A
Other languages
English (en)
Other versions
NO20065127L (no
Inventor
Jr William Henry Dragoset
Ian Moore
Richard Bisley
Original Assignee
Western Geco Seismic Holdings 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 Western Geco Seismic Holdings Ltd filed Critical Western Geco Seismic Holdings Ltd
Publication of NO20065127L publication Critical patent/NO20065127L/no
Publication of NO339499B1 publication Critical patent/NO339499B1/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
    • 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/56De-ghosting; Reverberation compensation

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Oceanography (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Description

Utførelser av den foreliggende oppfinnelse vedrører generelt marinseismiske undersøkelser og, mer særskilt, en fremgangsmåte for svekking av innvirkningen til overflatemultipler i seismiske signaler.
Seismiske undersøkelser er en fremgangsmåte for bestemmelse av strukturen til undergrunnsformasjoner. En seismisk undersøkelse benytter typisk seismikkenergikilder som genererer seismiske bølger, og seismikkmottakere som detekterer seismiske bølger. De seismiske bølger forplanter seg i formasjonene i grunnen, hvor en del av bølgene reflekteres fra grensesjiktene mellom undergrunnsformasjoner. Amplituden og polariteten til de reflekterte bølger bestemmes av forskjellene i akustisk impedans mellom berglagene som utgjør de undergrunnsformasjoner. Den akustiske impedansen til et berglag er produktet av den akustiske forplantningshastighet i laget og lagets tetthet. De seismiske mottakerne detekterer de reflekterte seismiske bølger og omdanner dem til representative elektriske signaler. Signalene blir typisk overført med elektriske, optiske, radio- eller andre midler til innretninger som registrerer signalene. Ved å analysere de registrerte signaler (eller traser), kan formen, posisjonen og sammensetningen av de undergrunnsformasjoner bestemmes.
En marinseismisk undersøkelse er en fremgangsmåte for bestemmelse av strukturen til undergrunnsformasjoner under vannmasser. En marinseismisk undersøkelse benytter typisk i vannet anordnede seismikkenergikilder og seismikkmottakere som enten slepes etter et fartøy eller plasseres på bunn i vannet ut fra et fartøy. Energikilden er typisk en eksplosiv innretning eller et trykkluftsystem som genererer seismisk energi, som så forplanter seg i vannmassen og inn i de underliggende grunnformasjoner. Når seismikkbølgene treffer grensesjikt mellom undergrunnsformasjoner, blir en del av seismikkbølgene reflektert tilbake gjennom jorden og vannet og til seismikkmottakerne, for detektering, overføring og registrering. De seismiske mottakere som typisk benyttes i marine seismikkundersøkelser, er trykksensorer, så som hydrofoner. I tillegg kan det benyttes bevegelsessensorer, så som akselerasjonsmålere. Både kilder og mottakere kan strategisk omplasseres for dekking av undersøkelsesområdet.
Seismikkbølger blir imidlertid reflektert fra andre grensesjikt enn de som foreligger mellom de ønskede undergrunnsformasjoner. Seismikkbølger blir også reflektert fra bunnen og fra vannflaten, og de resulterende bølger blir på sin side reflektert osv. Bølger som reflekteres flere ganger benevnes som "multipler". Bølger som reflekteres flere ganger i vannlag mellom vannflaten og bunnen benevnes som "vann-bunn-multipler". Vann-bunn-multipler har lenge vært betraktet som et problem ved marin seismikkbehandling og tolking, og det har derfor for håndtering av vann-bunn-multipler vært utviklet multippelsvekkingsmetoder basert på bølgeligninger. Imidlertid kan det defineres et større sett av multipler, innbefattende vann-bunn-multipler som et subsett. Det større sett inneholder multipler med oppadrettede refleksjoner fra grensesjikt mellom undergrunnsformasjoner i tillegg til oppadrettede refleksjoner fra bunnen. Multiplene i det større sett har til felles deres nedadrettede refleksjoner ved vannflaten og benevnes derfor som "overflatemultipler". Fig. 1, omtalt nedenfor, viser eksempler på ulike typer refleksjoner.
Fig. 1 viser rent skjematisk en marinseismisk undersøkelse. Prosedyren er generelt betegnet med henvisningstallet 100. Undergrunnsformasjoner som skal undersøkes, så som formasjonene 102 og 104, ligger under en vannmasse 106. Seismikkenergikilder 108 og seismikkmottakere 110 er plassert i vannmassen 106, typisk ved hjelp av ett eller flere seismikkfartøy (ikke vist). En seismikkilde, så som en luftkanon, tilveiebringer seismikkbølger i vannmassen 106, og en del av seismikkbølgene går ned gjennom vannet og mot de undergrunnsformasjoner 102 og 104 under vannmassen 106. Når seismikkbølgene når en seismikkreflektor vil en del av seismikkbølgene reflekteres oppover mens en del av seismikkbølgene fortsetter nedover. Seismikkreflektoren kan være vannbunnen 112 eller ett av grensesjiktene mellom to undergrunnsformasjoner, så som grensesjiktet 114 mellom formasjonene 102 og 104. Når de oppadrettede, reflekterte bølger når grensesjiktet mellom vann og luft ved vannflaten 116, vil en større andel av bølgene reflekteres nedover igjen. På denne måten kan seismikkbølgene reflekteres flere ganger mellom oppadreflektorer, så som vannbunnen 112 eller underliggende formasjonsgrensesjikt, og nedad-reflektoren ved vannflaten 116, slik det vil bli beskrevet nærmere nedenfor. Hver gang de reflekterte bølger forplanter seg forbi en seismikkmottaker 110, vil mottakeren 110 avføle de reflekterte bølger og generere representative signaler.
Primærrefleksjoner er de seismikkbølger som er reflektert bare én gang, fra vannbunnen 112 eller et grensesjikt mellom undergrunnsformasjoner, før de detekteres i en seismikkmottaker 110. Et eksempel på en primærrefleksjon er vist med strålebanene 120 og 122 i fig. 1. Primærrefleksjonene inneholder den ønskede informasjon vedrørende de undergrunnsformasjoner som er hensikten med den marinseismiske undersøkelse. Overflatemultipler er de bølger som er reflektert flere ganger mellom vannflaten 116 og enhver oppadreflektor, så som vannbunnen 112 eller formasjonsgrensesjiktene, før de avføles i en mottaker 110. Et eksempel på en overflatemultippel som er en vann-bunn-multippel, er vist med strålebanene 130, 132, 134 og 136. Det punkt i vannflaten 116 hvor bølgen reflekteres nedover, blir generelt betegnet som et nedadrefleksjonspunkt 133. Overflatemultiplen som starter i strålebanen 130, er en multippel av størrelsesorden én, fordi multiplen inneholder en refleksjon fra vannflaten 116. To eksempler på generelle overflatemultipler med oppoverrettede refleksjoner fra både vannbunnen 112 og formasjonsgrensesjikt er vist med strålebanene 140, 142, 144, 146, 148 og 150 og med strålebanene 160, 162, 164, 166, 168 og 170. Begge disse sistnevnte to eksempler på overflatemultipler er multipler av størrelsesorden to, fordi multiplene inneholder to refleksjoner fra vannflaten 116. Generelt vil en overflatemultippel være av størrelsesordenen i når den inneholder i refleksjoner fra vannflaten 116. Overflatemultipler er ekstern støy som forstyrrer det ønskede primære refleksjonssignal.
Overflatemultippelsvekking er en forstakkinversering av et registrert bølgefelt som fjerner alle størrelsesordner av samtlige overflatemultipler som forefinnes i det marinseismiske signal. Til forskjell fra noen bølgeligningsbaserte multippelsvekkingsalgoritmer, krever overflatemultippelsvekking ingen modellering av eller anslag vedrørende posisjoner, former og refleksjonskoeffisienter for de multippelbevirkende reflektorer. Isteden baserer overflatemultippelsvekking seg på den interne fysiske overensstemmelse mellom primær- og multippelhendelser som må foreligge i ethvert skikkelig registrert marint datasett. Den informasjon som er nødvendig for overflatemultippelsvekkingbehandlingen ligger allerede i seismikkdataene.
Ulike tidligere metoder har vært forsøkt for fjerning av overflatemultipler fra registrerte traser. Man har eksempelvis merket seg at gangtiden for en vann-bunn-multippel er en funksjon av "avviket", dvs. avstanden mellom kilden og mottaker, og av antall ganger multiplen reflekteres fra overflaten. Dersom eksempelvis multiplen reflekteres fra overflaten én gang før den mottas i mikrofonen, og avviket er null, så vil multiplens gangtid være nøyaktig to ganger den for de prinsipale bølger. Dette faktum har vært benyttet i mange ulike forsøk på å fjerne multipler.
Andre metoder innbefatter bruk av komplekse strålefølgemetoder hvor det genereres en syntetisk multippelbølge som trekkes fra den virkelige bølgen for dermed å oppnå en antatt multippelfri registrering. Disse metoder er imidlertid tungvinte fordi de krever signifikant kjennskap til den undersjøiske struktur så vel som til havbunnen før den syntetiske bølgen kan genereres. Lignende syntetiske multipler kan genereres med mer nøyaktige metoder som ikke direkte innbefatter strålefølging, eksempelvis feltforplantningsmetoder, men disse krever også detaljert kjennskap til i det minste havbunnen, så vel som kjennskap til de undergrunnsgr ens esj iktene, og er derfor ikke så praktiske som ønsket.
Dagens overflatemultippel-prediksjonsalgoritmer krever konvolusjoner av par av traser, hvor mottakerlokaliseringen for en trase i paret er sammenfallende med skuddlokaliseringen for den andre trasen. Fordi kilde- og mottakerlokaliseringer i registrerte data sjelden stemmer nøyaktig overens, er det vanlig praksis å regulere datasettene til en nominell geometri, slik at man oppnår overenstemmelsen mellom kilde- og mottakerlokaliseringer. Multiplene blir så prediktert for denne regulære geometri, og blir deretter deregulert til den opprinnelige geometri før subtraksjon. Uheldigvis vil slik regulering og (særlig) dereguleringen ofte være unøyaktig, hvilket ofte medfører signifikante feil i de predikterte multipler.
US-4 887 243 foreslår en overflate-multippel-prediksjonsalgoritme der par av traser konvolveres. I US-4 887 243 kreves at mottakerlokasjonen for en trase i paret sammenfaller med skuddlokasjonen for den andre trasen. Siden kilde- og mottakerlokasjonene i registrerte data sjelden sammenfaller presist, er det vanlig praksis å regularisere datasettene til en nominell geometri slik at dette sammenfallet mellom kilde- og mottakerlokasj oner oppnås. Multiplene blir så prediktert for denne regulære geometrien, og deretter de-regularisert til den originale geometrien før subtrahering. Dessverre er regulariserings- og (særlig) de-regulariseringsprosessene ofte unøyaktige, noe som ofte fører til signifikante feil i de predikterte multiplene.
Det foreligger derfor et behov for en bedret fremgangsmåte for fjerning av registrerte multippel-overflaterefleksjonshendelser fra seismiske registreringer for seismikkdatabehandlingsformål.
En eller flere utførelser av oppfinnelsen tar sikte på å minimere feil som skyldes regulering og deregulering, ved at bruken av en regulering minimeres mens dereguleringen helt unngås. I én utførelse benytter oppfinnelsen data fra mange undergrunnslinjer på én gang, for derved å minimere reguleringen. Som sådan kan oppfinnelsen innbefatte konvolusjoner av traser fra ulike undergrunnslinjer. I en annen utførelse blir algoritmen ifølge oppfinnelsen, for bibehold av den effektivitet som er tilordnet behandling som baserer seg på undergrunnslinjer, formulert for å arbeide med par av undergrunnslinjer samtidig.
I nok en utførelse er oppfinnelsen rettet mot en generalisert implementering av en tredimensjonal overflaterelatert multippelprediksjon-(SMP)-algoritme, som ikke foretar noen antagelser vedrørende regulariteten eller fordelingen av traser i det registrerte datasett, dvs. at det ikke forefinnes et nominelt geometrikonsept. Som sådan blir det registrerte datasett helt enkelt behandlet som en samling av traser definert av deres kilde- og mottakerlokaliseringer. På denne måten blir én eller flere utførelser av oppfinnelsen konfigurert for prediktering av multipler for et hvilket som helst sett av traser som defineres med deres kilde- og mottakerlokaliseringer. Særlig gjelder at oppfinnelsen kan prediktere multipler med riktig lokalisering, avvik og asimut.
En eller flere utførelser av oppfinnelsen kan også være utformet for å ta hensyn til samtlige irregulariteter i innhentingsgeometrien, hvilket vil kunne være særlig fordelaktig når man skal ta hensyn til kabelavviksvinkler. Nøyaktigheten av de predikterte multipler vil kunne være avhengig av fordelingen av trasene i det registrerte datasett.
Derfor er én eller flere utførelser av oppfinnelsen rettet mot en fremgangsmåte for prediktering av et antall overflatemultipler for et antall måltraser i en registrering av seismikkdata. I én utførelse innbefatter fremgangsmåten tilveiebringelsen av en fil som inneholder informasjon vedrørende et antall par av registrerte traser. Hvert par registrerte traser er i hovedsaken nærmest en ønsket skuddsidetrase og en ønsket mottakersidetrase. Fremgangsmåten innbefatter videre en konvolvering av parene av registrerte traser for derved å generere et antall konvolusjoner og stakke konvolusj onene for hver måltrase.
I en annen utførelse innbefatter fremgangsmåten (a) valg av en måltrase, (b) valg av et potensielt nedadrefleksjonspunkt for den valgte måltrase, (c) beregning av et ønsket skuddsidemidtpunkt, avvik og/eller asimut, og minst ett ønsket mottakersidemidtpunkt, avvik og/eller asimut, under utnyttelse av det valgte potensielle nedadrefleksjonspunkt og den valgte måltrasen. Det ønskede skuddsidemidtpunkt, avvik og asimut definerer en ønsket skuddsidetrase. Det ønskede mottakersidemidtpunkt, avvik og asimut definerer en ønsket mottakersidetrase. Fremgangsmåten innbefatter videre (d) bestemmelse av et par registrerte traser som i hovedsaken er nærmest den ønskede skuddsidetrase og den ønskede mottakersidetrase, og (e) konvolvering av paret av registrerte traser for generering av en konvolusj on.
Oppfinnelsen skal nå forklares nærmere under henvisning til utførelseseksempler, hvorav noen er vist på tegningen. Det skal imidlertid nevnes at tegningen bare viser typiske utførelser av oppfinnelsen og ikke er ment å være innsnevrende eller begrensende, idet oppfinnelsen også kan innbefatte andre like effektive utførelser.
På tegningen viser:
Fig. 1 er et skjematisk riss av en marinseismisk undersøkelse.
Fig. 2 viser et flytskjema for det første trinnet i en fremgangsmåte for gjennomføring av en tredimensjonal overflatemultippelprediksjon i samsvar med én eller flere utførelser av oppfinnelsen. Fig. 3A og 3B viser et flytskjema for det andre trinnet i en fremgangsmåte for gjennomføring av en tredimensjonal overflatemultippelprediksjon i samsvar med én eller flere utførelser av oppfinnelsen. Fig. 4 viser et grunnriss av en innhentingsgeometri i samsvar med én eller flere utførelser av oppfinnelsen. Fig. 5 viser et datamaskinnettverk hvor ulike utførelser av oppfinnelsen kan implementeres. Fig. 2 viser et flytskjema for det første trinnet 200 i en fremgangsmåte for gjennomføring av en tredimensjonal overflatemultippelprediksjon i samsvar med én eller flere utførelser av oppfinnelsen. I trinn 210 velges en måltrase. Et eksempel på en valgt måltrase er vist i fig. 4 som trasen (S, R). Måltraser definerer de lokaliseringer hvor multiplene skal predikteres. I trinnet 220 bestemmes eller defineres en apertur 305 for den valgte måltrasen. Aperturen 305 kan være et rektangulært område og være sentrert om et midtpunktsted M på måltrasen. Andre geometriske former av aperturen 305 vil kunne være mulige i utførelser av oppfinnelsen. Aperturen 305 defineres til å innbefatte i hovedsaken samtlige
potensielle nedadrefleksjonspunkter (DRP'er) til overflatemultiplene for måltrasen. Som eksempel er det i fig. 4 vist et potensielt nedadrefleksjonspunkt X. I trinnet 230 blir aperturen 305 gitret i et antall celler. I én utførelse er midtpunktet til måltrasen plassert i én av gitternodene (cellesenter). Gitteravstanden kan være vilkårlig. Gitternodene definerer de potensielle DRP'er for måltrasen.
I trinn 240 velges et potensielt DRP, så som et første DRP, for den valgte måltrasen. I trinnet 250 beregnes det ønskede skuddsidemidtpunkt Ms, avviket Xs og asimut 9s samt det ønskede mottakersidemidtpunkt MR, avviket XRog asimut 9R. Ms er midtpunktlokaliseringen mellom kilden og det valgte potensielle DRP. Avviket Xs er den horisontale avstanden mellom det valgte potensielle DRP og kilden S. Asimut 9s er vinkelen mellom den linjen som forbinder kilden S og det valgte potensielle DRP og en fast retning, som typisk vil være den såkalte in-line-retning. MR er midtpunktlokaliseringen mellom mottakeren R og det valgte potensielle DRP. Avviket XRer den horisontale avstand mellom det valgte potensielle DRP og mottakeren R. Asimut 9r er vinkelen mellom linjen som forbinder mottakeren R og det valgte potensielle DRP og en fast retning, som typisk er in-line-retningen. I én utførelse blir det ønskede skuddsidemidtpunkt Ms, avviket Xs og asimut 9s og det ønskede mottakersidemidtpunkt MR, avviket XRog asimut 9R beregnet basert på den valgte måltrase og det valgte potensielle DRP. Midtpunktene, avvikene og asimutene definerer sammen den ønskede skuddsidetrase (S, X) og den ønskede mottakersidetrase (X, R).
I trinn 260 blir den inngangstrasen som er nærmest den ønskede skuddsidetrasen og den inngangstrasen som er nærmest den ønskede mottakersidetrasen bestemt. I én utførelse bestemmes de nærmeste inngangstraser ved minimering av en objektiv funksjon, som definerer nærheten til to traser basert på deres midtpunkter, avvik og asimuter. Et eksempel på en objektivfunksjon er
D<2>= |Am|<2>+ wx|Ax|<2>+ we|A9|<2>
hvor D er et mål for nærheten mellom trasene, Am, Ax og A9 er forskjeller i henholdsvis midtpunkt, avvik og asimut, og wxog we er vekter som definerer den relative viktigheten av feil i avvik og asimut sammenlignet med feil i midtpunkter. wxer dimensjonsløs mens we har dimensjonen L2. I én utførelse settes we lik null som følge av dårlig asimutdekning i inngangsdatasettet. I en annen utførelse kan det foreligge en minimumsverdi for den minimerte objektivfunksjon, over hvilken minimumsverdi det anses at det ikke finnes en tilpasset trase.
I trinn 270 blir informasjon vedrørende de nærmeste inngangstraser lagret i en fil, som kan betegnes som en konvolusjon-indeksfil (CIF). Eksempelvis kan slik informasjon innbefatte identifikasjoner for de nærmeste inngangstraser, deres tilordnede undergrunns linjer, det valgte potensielle nedadrefleksjonspunkt X, det ønskede skuddsidemidtpunkt Ms, avvik Xs og asimut 9s, det ønskede mottakersidemidtpunkt MR, avvik XR og asimut 0R, og den valgte måltrase som skal predikteres. I trinnet 280 tas det en bestemmelse hvorvidt aperturen innbefatter et annet potensielt DRP for den valgte måltrasen. Dersom svaret er positivt går behandlingen tilbake til trinn 240, hvor det velges et annet potensielt DRP. Er svaret negativt så går behandlingen til trinn 285, hvor det foretas en bestemmelse hvorvidt det eksisterer en annen måltrase. Er svaret positivt, så går behandlingen tilbake til trinn 210, for valg av en annen måltrase. Er svaret negativt, så går behandlingen videre i trinn 290, hvor CIF deles i én eller flere subfiler i samsvar med par av undergrunnslinjer som inneholder nærmeste inngangstraser. På denne måten vil hver subfil inneholde informasjon vedrørende bare ett par undergrunnslinjer, idet hver undergrunnslinje inneholder en inngangstrase som er nærmest enten en ønsket skuddsidetrase eller en ønsket mottakersidetrase. Ordningen av undergrunnslinjer i paret er ikke kritisk.
Fig. 3A og 3B viser et flytskjema for det andre trinnet 300 i en fremgangsmåte for
gjennomføring av en tredimensjonal overflatemultippelprediksjon i samsvar med én eller flere utførelser av oppfinnelsen. I trinn 310 velges den første subfil. I trinn 320 leses informasjon fra den valgte subfil vedrørende et par av inngangstraser nærmest en ønsket skuddsidetrase og en ønsket mottakersidetrase for en valgt måltrase. I trinn 330 blir et par inngangstraser korresponderende med informasjonen vedrørende paret av nærmeste inngangstraser, uttrukket fra et sett av registrerte seismikkdata.
Det registrerte sett av seismikkdata kan lagres i enhver fil eller ethvert datalager som vil være kjent for fagpersonen. Settet med registrerte seismikkdata kan ekstrapoleres til null-avvik. Settet av registrerte seismikkdata kan være en samling av forstakktraser definert av midtpunkt, avvik og asimut. Hver trase i det registrerte seismikkdatasett kan ha en underjordisk linjeidentifikasjon og en unik traseidentifikasjon som kan benyttes for å identifisere inngangstrasen i det registrerte seismikkdatasett. Settet med registrerte seismikkdata kan være organisert i undergrunns linjer eller andre underoppdelinger, så som seilingslinjer.
I trinnet 340 benyttes en differensial-utflytingskorreksjon på paret av uttrukkede registrerte traser for korreksjon av avvikene til de uttrukkede registrerte traser med hensyn til det ønskede skuddsideavvik og det ønskede mottakersideavvik. I trinn 250 blir paret av korrigerte og uttrukkede registrerte traser konvolvert. I trinn 355 lagres konvolusj onen.
I trinnet 360 tas det en bestemmelse hvorvidt den valgte subfil inneholder et annet par inngangstraser som skal konvolveres. Er svaret positivt, så går behandlingen tilbake til trinn 320. Er svaret negativt, så fortsetter behandlingen i trinn 365, hvor konvolusj onene sorteres i samsvar med måltrasene. I trinnet 370 blir samtlige konvolusjoner for hver måltrase stakket for oppnåelse av én enkelt, stakket konvolusj on pr. måltrase for den valgte subfil.
I trinn 375 treffes det en avgjørelse hvorvidt det eksisterer en annen subfil av CIF. Er svaret positivt, så velges den subfilen (trinn 378) og behandlingen går tilbake til trinn 320. Er svaret negativt, så fortsetter behandlingen i trinn 380, hvor samtlige stakkede konvolusjoner fra enhver subfil i CIF sorteres i samsvar med måltrasene. I trinnet 385 blir de stakkede konvolusjoner for hver subfil i CIF stakket for hver måltrase for oppnåelse av én enkelt, stakket konvolusj on pr. måltrase fra alle subfiler.
I trinn 390 blir kildesignaturen dekonvolvert i samsvar med kjent teknikk. I trinn 395 benyttes et tredimensjonalt p-filter for korrigering av stakkings virkning en på småbølgen, i samsvar med teknikk som vil være velkjent for en fagperson.
Fig. 5 viser et datamaskinnettverk 500, hvor utførelser av oppfinnelsen kan implementeres. Datamaskinnettverket 500 innbefatter en systemdatamaskin 530, som kan være en hvilken som helst konvensjonell PC eller arbeidsstasjon, så som en UNIX-basert arbeidsstasjon. Systemdatamaskinen 530 kommuniserer med disklagringsinnretninger 529, 531 og 533, som kan være eksterne harddisk-lagringsinnretninger. Disklagringsinnretningene 529, 531 og 533 kan være konvensjonelle harddiskdrifter og kan som sådanne implementeres ved hjelp av et lokalt områdenettverk eller med fjernaksess. Selv om disklagringsinnretningene 529, 531 og 533 er vist som separate innretninger, kan det selvfølgelig benyttes én enkelt disklagringsinnretning for lagring av programinstruksjoner, måledata og resultater, alt etter behov.
I én utførelse lagres seismikkdata fra hydrofoner i disklagringsinnretningen 521. Systemdatamaskinen 530 kan motta egnede data fra disklagringsinnretningen 531 for gjennomføring av en 3-D-overflatemultippelprediksjon i samsvar med programinstruksjoner som svarer til de her beskrevne fremgangsmåter. Programinstruksjonene kan være skrevet i et datamaskinprogrammeringsspråk, så som C++, Java og lignende. Programinstruksj onene kan lagres i en datamaskinlesbar hukommelse, så som en programdisklagringsinnretning 533. Naturligvis kan det hukommelsesmedium som brukes for lagring av programinstruksj onene være av en hvilken som helst konvensjonell type som benyttes for lagring av datamaskinprogram, herunder harddiskdrifter, floppy disker, CD-ROM og andre optiske media, magnetbånd, og lignende.
Ifølge en foretrukket utførelse av oppfinnelsen har systemdatamaskinen 530 primært en grafikkvisning 527, eller alternativt en utgang via en printer 528. Systemdatamaskinen 530 kan lagre resultatene fra de foran beskrevne fremgangsmåter i disklageret 529, for senere bruk og videre analyse. Tastaturet 526 og pekeinnretningen (eksempelvis en mus, en kule eller lignende) 525 kan være tilordnet systemdatamaskinen 530 for derved å muliggjøre en interaktiv bruk.
Systemdatamaskinen 530 kan være plassert i et datasenter i en avstand fra undersøkelsesområdet. Systemdatamaskinen 530 har forbindelse med hydrofoner
(enten direkte eller via en registreringsenhet, ikke vist), for mottak av signaler som er indikative for den reflekterte seismikkenergi. Disse signaler, etter konvensjonell formatering og annen behandling, lagres av systemdatamaskinen 530 som digitale data i disklageret 531, for senere uttak og behandling som beskrevet foran. Mens fig. 5 viser et disklager 531 direkte forbundet med systemdatamaskinen 530, vil det også være mulig at disklagringsinnretningen 531 er tilgjengelig via et lokalt områdenettverk eller med fjernaksess. Videre, selv om disklagringsinnretningene 529, 531 er vist som separate innretninger for lagring av seismiske inngangsdata og analyseresultater, kan disklagringsinnretningene 529, 531 være implementert i én enkelt disk-drift (enten sammen med eller atskilt fra
programdisklagringsinnretningen 533), eller på annen konvensjonell måte som fagpersoner vil ha kjennskap til.

Claims (12)

1. Fremgangsmåte for prediktering av et antall overflatemultipler for et antall måltraser i en registrering av seismikkdata, innbefattende: (a) å velge en måltrase; (b) å velge et potensielt nedadrefleksjonspunkt for den valgte måltrase; (c) å beregne minst en av et ønsket skuddsidemidtpunkt, -avvik og -asimut, og minst en av et ønsket mottakersidemidtpunkt, -avvik og -asimut, basert på den valgte måltrase og det valgte potensielle nedadrefleksjonspunktet, hvor de ønskede skuddsidemidtpunkt, -avvik og -asimut definerer den ønskede skuddsidetrase, og de ønskede mottakssidemidtpunkt, -avvik og -asimut definerer den ønskede mottaksidetrase; (d) å bestemme et par av registrerte traser i hovedsak nærmest den ønskede skuddsidetrase og den ønskede mottakssidetrase; (e) å konvolvere paret av registrerte traser for å generere en konvolusj on; (f) å gjenta trinnene (a) til og med (e) for alle de potensielle nedadrefleksjonspunktene for hver måltrase; og (g) å stakke et flertall av konvolusjoner for hver måltrase.
2. Fremgangsmåte ifølge krav 1, videre innbefattende: bestemmelse av en apertur for den valgte måltrase.
3. Fremgangsmåte ifølge krav 2, videre innbefattende en gitring av aperturen i et antall celler.
4. Fremgangsmåte ifølge krav 2, hvor aperturen utformes for å inkludere i hovedsaken alle de potensielle nedadrefleksjonspunkter av overflatemultipler for den valgte måltrase.
5. Fremgangsmåte i samsvar med krav 1, hvor bestemmelsen av par av registrerte traser innbefatter minimering av en objektfunksjon.
6. Fremgangsmåte ifølge krav 1, videre omfattende å danne en fil som inneholder informasjon vedrørende paret av registrerte traser.
7. Fremgangsmåte ifølge krav 6, videre omfattende å inndele filen i én eller flere subfiler.
8. Fremgangsmåte ifølge krav 6, videre omfattende å inndele filen i én eller flere subfiler i samsvar med undergrunns-linjepar, idet hvert par av registrerte traser i hver subfil kommer fra et par undergrunns-linjer.
9. Fremgangsmåte ifølge krav 6, videre omfattende uttrekking av et antall registrerte traser som korresponderer med informasjonen vedrørende parene av registrerte traser.
10. Fremgangsmåte ifølge krav 9, videre omfattende korrigering av avvikene til de uttrukkede registrerte traser til et avvik for den ønskede skuddsidetrase og et avvik for den ønskede mottakersidetrase.
11. Fremgangsmåte ifølge krav 9, videre omfattende å benytte en differensial-utflyttingskorreksjon på de uttrukkede registrerte traser for korrigering av avvikene til de uttrukkede registrerte traser.
12. Fremgangsmåte ifølge krav 1, videre omfattende å sortere konvolusj onene i samsvar med måltrasene.
NO20065127A 2004-04-07 2006-11-07 Generalisert 3D-overflate-multippel prediksjon NO339499B1 (no)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US56012904P 2004-04-07 2004-04-07
PCT/US2004/023119 WO2005103764A1 (en) 2004-04-07 2004-07-16 Generalized 3d surface multiple prediction

Publications (2)

Publication Number Publication Date
NO20065127L NO20065127L (no) 2006-11-07
NO339499B1 true NO339499B1 (no) 2016-12-19

Family

ID=34958345

Family Applications (1)

Application Number Title Priority Date Filing Date
NO20065127A NO339499B1 (no) 2004-04-07 2006-11-07 Generalisert 3D-overflate-multippel prediksjon

Country Status (8)

Country Link
US (1) US7796467B2 (no)
EP (1) EP1735640B1 (no)
AT (1) ATE417288T1 (no)
AU (1) AU2004318849B2 (no)
MX (1) MXPA06011579A (no)
NO (1) NO339499B1 (no)
RU (1) RU2339056C2 (no)
WO (1) WO2005103764A1 (no)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7257492B2 (en) 2005-08-26 2007-08-14 Westerngeco L.L. Handling of static corrections in multiple prediction
US7987054B2 (en) 2008-05-23 2011-07-26 Exxonmobil Upstream Research Company Efficient multiple prediction in two and three dimensions
US8164978B2 (en) * 2008-09-19 2012-04-24 Westerngeco L.L.C. Interbed multiple prediction
US8300498B2 (en) * 2009-06-30 2012-10-30 Pgs Geophysical As Method for dynamic aperture determination for three-dimensional surface-related multiple elimination
US8478531B2 (en) * 2010-01-20 2013-07-02 Pgs Geophysical As Dip-based corrections for data reconstruction in three-dimensional surface-related multiple prediction
US8126652B2 (en) * 2010-01-20 2012-02-28 Pgs Geophysical As Azimuth correction for data reconstruction in three-dimensional surface-related multiple prediction
CN102236101B (zh) * 2010-05-07 2013-03-13 中国石油天然气股份有限公司 多层裂缝预测方法和装置
US8437218B2 (en) * 2010-06-29 2013-05-07 Westerngeco Llc Correcting geometry-related time and amplitude errors
US10295688B2 (en) 2010-08-10 2019-05-21 Westerngeco L.L.C. Attenuating internal multiples from seismic data
EP2831635B1 (en) 2012-03-29 2020-03-11 Geco Technology B.V. Seismic noise removal
US10054704B2 (en) 2013-07-01 2018-08-21 Westerngeco L.L.C. Predicting multiples in survey data
WO2015136379A2 (en) * 2014-03-14 2015-09-17 Cgg Services Sa Method and apparatus for estimating source signature in shallow water
US9928315B2 (en) 2014-07-30 2018-03-27 Chevron U.S.A. Inc. Re-ordered interpolation and convolution for faster staggered-grid processing
US10416326B2 (en) * 2015-12-18 2019-09-17 Schlumberger Technology Corporation Method for predicting multiples in survey data
US10295687B2 (en) 2016-08-30 2019-05-21 Schlumberger Technology Corporation Attenuation of multiple reflections
EP3649491B1 (en) 2017-07-05 2022-09-14 Services Pétroliers Schlumberger Reflection seismology internal multiple estimation
CN108828664B (zh) * 2018-06-07 2019-12-10 中国石油天然气股份有限公司 一种多次波识别方法及装置

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4887243A (en) * 1982-12-27 1989-12-12 Mobil Oil Corporation Removal of surface multiples

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3371310A (en) 1966-03-14 1968-02-27 Pan American Petroleum Corp Discriminating between primary and multiple seismic reflections
GB9612470D0 (en) * 1996-06-14 1996-08-14 Geco As Multiple attenuation method
US6094620A (en) * 1999-01-08 2000-07-25 Exxonmobil Upstream Research Company Method for identifying and removing multiples from seismic reflection data
US6169959B1 (en) 1999-05-27 2001-01-02 Baker Hughes Incorporated Method of predicting kinematics for surface multiples
EP1405105A2 (en) * 2001-05-25 2004-04-07 ExxonMobil Upstream Research Company Trapped water bottom multiple and peg-leg multiple suppression for ocean bottom seismic data
US6735527B1 (en) * 2003-02-26 2004-05-11 Landmark Graphics Corporation 3-D prestack/poststack multiple prediction
US7181347B2 (en) * 2003-09-23 2007-02-20 Westerngeco, L.L.C. Method for the 3-D prediction of free-surface multiples
EP1738201B1 (en) * 2004-04-07 2014-04-09 WesternGeco Seismic Holdings Limited Fast 3-d surface multiple prediction

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4887243A (en) * 1982-12-27 1989-12-12 Mobil Oil Corporation Removal of surface multiples

Also Published As

Publication number Publication date
AU2004318849B2 (en) 2008-11-13
EP1735640A1 (en) 2006-12-27
RU2006138605A (ru) 2008-05-20
ATE417288T1 (de) 2008-12-15
NO20065127L (no) 2006-11-07
EP1735640B1 (en) 2008-12-10
US20080043573A1 (en) 2008-02-21
MXPA06011579A (es) 2008-03-11
RU2339056C2 (ru) 2008-11-20
AU2004318849A1 (en) 2005-11-03
US7796467B2 (en) 2010-09-14
WO2005103764A1 (en) 2005-11-03

Similar Documents

Publication Publication Date Title
NO339499B1 (no) Generalisert 3D-overflate-multippel prediksjon
NO339072B1 (no) Hurtig 3-D-overflatemultippel-prediksjon
US8437218B2 (en) Correcting geometry-related time and amplitude errors
US9013956B2 (en) Method and system for seismic imaging and earth modeling using beam tomography
CA2801531C (en) Method and system for determining source signatures after source ghost removal
US10605937B2 (en) Device and method for smart picking surface waves dispersion curves
NO337875B1 (no) Fremgangsmåte og medium for forutsigelse av fri-overflate-multipler i 3D seismiske data
NO339243B1 (no) Fremgangsmåte for svekking av multipler fra vannlag i en samling seismiske datatraser
US20080130411A1 (en) Seismic imaging with natural green&#39;s functions derived from vsp data
US8300498B2 (en) Method for dynamic aperture determination for three-dimensional surface-related multiple elimination
NO343375B1 (no) Seismiske marine langtidsundersøkelser som benytter interpolerte multikomponents streamer-trykkdata
NO339050B1 (no) Fremgangsmåte for prosessering av seismikkdata
US9581711B2 (en) Interbed multiple prediction
US9733374B2 (en) Repeatability indicator based on shot illumination for seismic acquisition
GB2462195A (en) Method of summing dual-sensor towed streamer signals using seismic reflection velocities
WO2013093468A2 (en) Full waveform inversion quality control method
WO2013093467A1 (en) Method of, and apparatus for, full waveform inversion
CN105137479A (zh) 一种面元覆盖次数的计算方法及装置
US20130208566A1 (en) Method for processing dual-sensor streamer data with anti-alias protection
GB2503640A (en) Quality Assurance in a Full Waveform Inversion Process
EP0297852A2 (en) Method for real time display of marine seismic survey data coverage
Alsadi et al. Seismic Wave Transmission and Refraction