NO20131614L - Fremgangsmåte for marine seismiske undersøkelser - Google Patents

Fremgangsmåte for marine seismiske undersøkelser

Info

Publication number
NO20131614L
NO20131614L NO20131614A NO20131614A NO20131614L NO 20131614 L NO20131614 L NO 20131614L NO 20131614 A NO20131614 A NO 20131614A NO 20131614 A NO20131614 A NO 20131614A NO 20131614 L NO20131614 L NO 20131614L
Authority
NO
Norway
Prior art keywords
noise
hydrophones
seismic
output signal
signal
Prior art date
Application number
NO20131614A
Other languages
English (en)
Other versions
NO335663B1 (no
Inventor
Ali Ozbek
James Martin
Nils Lunde
Simon Hastings Bittleston
Original Assignee
Schlumberger Holdings
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
Publication of NO20131614L publication Critical patent/NO20131614L/no
Application filed by Schlumberger Holdings filed Critical Schlumberger Holdings
Publication of NO335663B1 publication Critical patent/NO335663B1/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/364Seismic filtering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/16Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
    • G01V1/20Arrangements of receiving elements, e.g. geophone pattern
    • G01V1/201Constructional details of seismic cables, e.g. streamers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/30Noise handling
    • G01V2210/32Noise reduction
    • G01V2210/324Filtering
    • G01V2210/3246Coherent noise, e.g. spatially coherent or predictable

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (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)
  • Geophysics And Detection Of Objects (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

En fremgangsmåte og et system for å utføre en marin seismisk undersøkelse er beskrevet, inkludert tauing av minst en seismisk streamer som omfatter et flertall av hydrofoner distribuert med intervaller på i gjennomsnitt ikke mer enn 500 cm langs denne, i vannet over området som skal undersøkes; å rette akustiske signaler ned gjennom vannet og inn i den underliggende grunnen; å motta, med hydrofonene, seismiske signaler reflektert fra lag i grunnen undervannet; å digitalisere utgangssignalet fra hver hydrofon separat; og å filtrere dette utgangssignalet for å redusere støyen som er tilstede i utgangssignalet og å generere et signal med et redusert støyinnhold, idet filtreringsprosessen benytter, som et ytterligere inngangssignal, det digitale utgangssignalet til i det minste én nærliggende hydrofon. Filtreringen anvendes ved enkeltsensorregistreringer før gruppedannelse og der dermed i stand til å detektere og redusere innholdet av støy med en koherenslengde på 20 meter eller mindre. Den reduserer støy slik som streamer- eller utbulingsstøy.

Description

Foreliggende oppfinnelse angår marine seismiske undersøkelser, og er mer spesielt knyttet til støyreduksjonsmetoder og anordninger for bruk ved marine seismiske undersøkelser.
For å utføre en 3D marin seismisk undersøkelse taues et flertall seismiske streamere (hydrofonkabler), hver typisk flere tusen meter lange og med kjeder av hydrofoner og tilhørende elektronisk utstyr distribuert langs sin lengde, med omkring 5 knop bak et seismisk undersøkelsesfartøy, som også tauer en eller flere seismiske kilder, vanligvis luftkanoner. Akustiske signaler produsert av de seismiske kildene rettes ned gjennom vannet inn i grunnen under, hvor de reflekteres fra de forskjellige lag (strata). De reflekterte signalene mottas av hydrofonene i streamerne, digitaliseres og sendes så til det seismiske undersøkelsesfartøy et, hvor de registreres og i det minste gjennomgår en delvis behandling med det endelige mål å bygge opp en representasjon av lagene i området som undersøkes.
Typisk taues det opp til 12 streamere, hvor hver streamer er flere kilometer lang. Streamerne utgjøres av seksjoner som typisk kan være 100-200 meter lange; hver seksjon består av hydrofoner inne i en ytre kappe (skin) som kan være fylt med olje, skum, eller en mer solid substans. Belastningsvaiere og avstandsstykker utgjør det indre skjelettet til streameren.
Signalene som mottas av hydrofonene i streamerne er alltid forurenset av støy fra et flertall kilder. Fartøyets krenging, spesielt i høy sjø, forårsaker vibrasjoner i belastningsvaierne, som virker inn på sammenkoblingene og den oljefylte kappen, noe som genererer utbulingsbølger (bulge waves) (eller pustende bølger (breathing waves)) som forplanter seg ned langs streamerne. Trykkvariasj onene detekteres av hydrofonene, noe som tilføyer støy og forvrenger de registrerte seismiske signalene. Når streameren beveger seg gjennom vannet, forårsaker grensesjiktturbulens trykkfluktuasjoner ved ytterveggen til kappen, som igjen kobles til hydrofonene.
Utbulingsbølger kan også forårsakes av virvelløsriving ved elliptisk vannbevegelse rundt streameren forårsaket av bølger.
For tiden omfatter en av de viktigste teknikkene for å redusere denne støyen, fast sammenkobling (hard-wiring) av tilgrensende hydrofoner for å summere deres respektive analoge utgangssignaler: typisk inneholder en gruppe åtte uniformt fordelte hydrofoner, og sentrene til hydrofonene ligger typisk med intervaller på 6,25 meter. Et slikt arrangement er beskrevet i vårt US-patent 5,351,218, som også beskriver hvordan par av tilgrensende grupper kan forbindes med hverandre selektivt for å danne grupper på seksten tilgrensende hydrofoner hvis gruppesentre ligger med 12,5 meters intervaller.
Siden de individuelle hydrofonene i hver gruppe ligger relativt tett, typisk bare 90 cm fra hverandre, antas det at alle hydrofonene i en gitt gruppe mottar i hovedsak det samme seismiske signalet. Det seismiske signalet forsterkes derfor ved summering av de analoge utgangssignalene fra hydrofonene i gruppen, mens støyen som påvirker hver hydrofon, dersom den er tilfeldig ukorrelert, vil tendere mot å kanselleres av summeringsprosessen. Gruppene av åtte eller seksten hydrofoner kan dermed anses ekvivalent med enkelte hydrofoner med en 6,25 eller 12,5 meter avstand, en forbedring på åtte eller seksten i forhold til en individuell hydrofon innenfor en gruppe, og gi en svært god avvisning av tilfeldig støy.
Imidlertid er en vesentlig kilde til støy som påvirker hydrofonene, bevegelsen av overflaten av vannet i området for undersøkelsen, på grunn av bølger og dønninger. Støy på grunn av bølger eller dønninger, som heretter vil refereres til simpelthen som "dønningsstøy", er ikke virkelig tilfeldig i forhold til gruppene av åtte eller seksten hydrofoner ifølge kjent teknikk, slik at summeringen av de analoge utgangssignalene til hydrofonene i hver gruppe ikke er særlig effektive i å redusere den. Videre øker betydningen av dønningsstøy dramatisk når høyden til bølgene eller dønningene øker, i den grad at når været forårsaker høyden til å overstige et visst nivå, typisk 2 til 4 meter, reduseres signal-/støy forholdet så mye at undersøkelsen må avbrytes til været bedres. Denne "weather downtime" kan øke de totale kostnadene ved undersøkelsen betydelig.
En metode for å anvende adaptiv signalbehandling ved dempning av utbulingsbølger er beskrevet i US-patent 4,821,241. Der er det foreslått å anordne belastningssensorer sammen med hydrofonene i streameren. Belastningssensorene er følsomme for mekaniske belastninger som påføres kabelen, men er i hovedsak ikke følsomme for akustiske bølger som forplanter seg i flytende medier. Utgangssignalene fra belastningssensorene kombineres med utgangssignalene fra de tilsvarende hydrofonene som er plassert på samme sted, for å kansellere falske signaler på grunn av utbulingsbølger.
En annen metode for å anvende adaptiv signalbehandling for å dempe utbulingsbølgene er beskrevet i US-patent 5,251,183.1 dette patentet er det foreslått å benytte et akselerometer festet mellom innføringsdelen av streameren og hydrofonen. Intra-shot- og inter-shot-akselerometer og hydrofonsignaler registreres. Metoden benytter inter-shot- og intra-shot adaptive behandlingssløyfer. Den inter-shot-adaptive behandlingssløyfen utleder inter-shot komplekse vekttall fra inter-shot akselerometersignaler og inter-shot hydrofonsignaler. Den intra-shot adaptive behandlingssløyfen modellerer utbulingsbølgestøy i intra-shot hydrofonsignalene ved å kombinere de inter-shot komplekse vekttallene med intra-shot akselerometersignalene. Utbulingsbølgestøydempingen oppnås ved å subtrahere intra-shotbløgestøymodellen fra intra-shot seismisk detektorsignalene.
Andre typer støy, slik som tverrstrømsstøy (crossflow noise) generert når streameren utsettes for tverrgående strømninger (cross-currents), har karakteristikker som er like med karakteristikkene til utbulingsbølgestøyen. Tverrstrømsstøy er et spesielt viktig problem når det seismiske undersøkelsesfartøyet og streamerne snus i løpet av en seismisk undersøkelse. Nivået av tverrstrømsstøy er typisk så overveldende ved snuing at registreringen av seismiske data simpelthen stoppes mens fartøyet snus. Kjente signalbehandlingsmetoder for seismiske data har ikke lykkes å oppnå en adekvat demping av disse typene støy heller.
Det er derfor et formål med oppfinnelsen å fremskaffe fremgangsmåter og anordninger for å redusere effektene av støy, slik som dønningsstøy, ved marine seismiske undersøkelser.
Ifølge en side ved foreliggende oppfinnelse er det fremskaffet et seismisk registreringssystem som omfatter: en streamer for mottak av seismiske signaler reflektert fra lag under en vannmengde, hvor streameren omfatter et flertall av hydrofoner distribuert ved gjennomsnittlige intervaller på ikke mer enn 500cm langs denne; analog-til-digitalomformerutstyr som mottar og digitaliserer analoge signaler fra hydrofonene for å generere, for hver hydrofon, et separat digitalt signal; og et filter som mottar som inngangssignal, det digitaliserte signalet sammen med det digitaliserte signalet fra i det minste én nærliggende hydrofon og som genererer et signal med redusert støyinnhold.
Ifølge en annen side ved oppfinnelsen er det fremskaffet en fremgangsmåte for å utføre en marin seismisk undersøkelse, hvor fremgangsmåten omfatter: - å taue minst en seismisk streamer med et flertall av hydrofoner fordelt ved gjennomsnittlige intervaller på ikke mer enn 500cm langs denne, i vannet over det området som skal undersøkes; - å rette akustiske signaler ned gjennom vannet og inn i grunnen under; - å motta, med hydrofonene, seismiske signaler reflektert fra lag i grunnen under vannet;
- å digitalisere et utgangssignal fra hver hydrofon separat; og
- å filtrere nevnte utgangssignal for å redusere støyen som er tilstede i nevnte utgangssignal og å generere et signal med et redusert støyinnhold, idet nevnte filtreringstrinn benytter, som et ytterligere inngangssignal, det digitale utgangssignalet til i det minste én nærliggende hydrofon.
Oppfinnelsen er basert på den oppdagelse at de dominerende støykomponentene typisk har en koherenslengde på mindre enn 20 meter i det lave frekvensområdet. Ved å sample bølgefeltet ved intervaller som er vesentlig under den avstanden for å unngå - eller i det minste redusere - aliasing, og å anvende en egnet støyfiltreringsteknikk, er det mulig å redusere mengden støy i registrerte streamerdata mer effektivt enn ved kjente gruppedannelsesmetoder. Koherenslengden til støyen bestemmes av forskjellige parametre, inkludert streamerdesign, konstruksjon, tauing og værforhold.
Spesielt foreslår oppfinnelsen å sample bølgefeltet ved en gjennomsnittlig samplingsavstand på mindre enn 500cm. Denne samplingstettheten kan ikke oppnås ved konvensjonell "gruppedannede" innsamlingsdata, hvor utgangene til tilgrensende hydrofoner er koblet sammen. Denne kjente teknikken midler målingene over koherenslengden til støyen, og gjør den dermed uegnet for konvensjonelle koherent støyfiltreringsmetoder.
Oppfinnelsen fremskaffer et system som kan konvertere det seismiske signalet slik det er målt av en enkelt hydrofon, til en støyredusert versjon av samme signal. Dette støyreduserte signalet kan være tilgjengelig for etterfølgende behandlingstrinn, enten i det såkalte "pre-stack" domene eller signalet kan kombineres med de støyreduserte signalene til andre hydrofoner i en prosess som normalt refereres til som "stacking".
Betegnelsen "tilgrensende hydrofoner" er ment å inkludere direkte tilstøtende hydrofoner, men også tilfeller hvor det benyttes inngangssignaler fra hydrofoner som ligger nærme (men ikke direkte ved siden av).
Fortrinnsvis ligger avstanden mellom hydrofonene i området 200cm til 330cm, og kan for eksempel være omkring 205cm til 210cm eller omkring 305cm til 315cm.
Romlig-tidsmessig (spatio-temporal) filtrering anvendes fortrinnsvis for å dempe støyen fra de mottatte signalene. Vanligvis er det ved seismiske undersøkelser mulig å skille generelt mellom støy og seismiske signaler i form av frekvens, forplantningsretning og (tilsynelatende) hastighet.
I overensstemmelse med en annen side ved foreliggende oppfinnelse virker filteret som en stråleformer (beamformer), dvs at den har en respons som er følsom for det romlige og/eller tidsmessige spektrale innholdet av inngangssignalene.
Fortrinnsvis er filteret et adaptivt filter, enda mer foretrukket er et filter som tilpasser sine filtervekttall eller koeffisienter under et forhåndsdefinert sett av begrensninger.
Ved en spesielt foretrukket utførelse omfatter filteret M romlig og/eller tidsmessig lokale flerkanals adaptive filtre med K kanaler, hver med en lengde L. For de fleste anvendelser er tallene M, K og L lik eller større enn to.
Bruken av en filterbank for støydemping av seismiske signaler er beskrevet i internasjonal patentsøknad WO 97/25632. Imidlertid krever ikke foreliggende oppfinnelse at det defineres en referansekanal for å beregne de adaptive filterbankkoeffisientene. Med andre ord inngår det ingen støyestimater i adapteringsprosessen. Derfor kan foreliggende fremgangsmåte anvendes ved støyforurensede seismiske signaler hvor det ikke er tilgjengelig noen uavhengig måling eller estimat av støyen.
Ifølge en side ved oppfinnelsen er koeffisientene til filterbanken begrenset slik at dets respons svarer til det for en stråleformer med en spesifisert synsretning.
Fremgangsmåten kan utføres på lagrede data eller på rå seismiske data etter behov. Slik kan rå seismiske data ifølge fremgangsmåten filtreres på datainnsamlingsstedet. Dette sikrer at et "renset" signal er tilgjengelig fra datainnsamlingsstedet og kan lastes ned direkte fra stedet i denne form. Dette reduserer mengden data som må sendes for analyse "off-site" og reduserer kostnadene og lagringsproblemene knyttet til akkumulering av tilstrekkelige mengder støyfylte data for analyse off-site. Fremgangsmåten kan anvendes for enkeltsensorregistreringer, dvs registreringer som går forut for dannelse av grupper som kombinerer signalene fra to eller flere seismiske sensorer.
Disse og andre trekk ved oppfinnelsen, foretrukne utførelsesformer og varianter av disse, mulige anvendelser og fordeler, vil oppfattes og forstås av fagmenn på området ut fra den følgende detaljerte beskrivelse og tegningene. Fig. IA er en generalisert illustrasjon av et seismisk undersøkelsesfartøy som tauer en sammensetning av seismiske streamere i overensstemmelse med en første utførelse av oppfinnelsen; Fig. IB er et skjematisk tverrsnitt som illustrerer en del av en streamerkabel som
benytter en utførelse av foreliggende oppfinnelse; og
Fig. 2 viser et generelt blokkdiagram av en adaptiv stråleformer i samsvar med
foreliggende oppfinnelse.
Fig. IA viser et seismisk fartøy som tauer en akustisk kilde og en streamer gjennom en vannmengde. Det seismiske fartøyet 10 trekker minst en seismisk kilde 11 og minst en seismisk streamer 12. Streameren 12 er festet til fartøyet 10 ved hjelp av en innføringskabel 13 som er festet til en kableoppbevaringstrommel 14 plassert om bord i fartøyet. En halebøye 15 er festet til den distale enden av kabelen ved hjelp av et langt tau av lignende materiale. Halebøyen kan alternativt være utrustet med en akustisk, elektromagnetisk eller visuell anordning til bruk for å lokalisere enden av streamerkabelen.
Som vist i fig. IB inneholder hver seksjon av streamerkabelen 12 inneholder et flertall hydrofoner 21 slik det er kjent fra det seismiske faget. Hydrofonene er sammenkoblet med en transmisjonslinje (ikke vist) til en fjerntliggende registreringsanordning plassert om bord i fartøyet. Tilgrensende hydrofoner er ikke koblet sammen slik at de danner grupper som genererer et enkelt utgangssignal. Isteden er hver hydrofon utformet for å avgi et separat utgangssignal gi, som deretter filtreres i en prosess som beskrevet under.
I tillegg inneholder streamerkabelen belastningsdeler 22 som er utformet for å absorbere strekkbelastninger som kabelen utsettes for når den taues. Registrerings-og styrkekomponentene er omgitt av en plasthylse 23 i form av et langstrakt rør. Hylsen er fortrinnsvis fylt med et lettvekts ballastfluid for å gjøre seksjonen nøytral eller med en viss oppdrift. Når den er full av fluid er det indre av hylsen i hovedsak ved atmosfærisk trykk. Den sylindriske formen til hylsen opprettholdes av et flertall av skillevegger som ikke er vist.
Når streamerkabelen taues gjennom vannmengdene, avfyres luftkanoner 11 og den akustiske energien som genereres på denne måten, forplanter seg gjennom vannlaget og formasjonene under havbunnen. Ved forskjellige refleksjonspunkter eller plan, reflekteres en del av den akustiske energien. Hydrofonene 21 kavelen mottar det direkte bølgefeltet og alle reflekterte og avbøyde bølgefelter som traverserer streameren. Det mottatte bølgefeltet er i de fleste tilfeller kraftig forurenset av støy fra forskjellige kilder.
For å dempe den uønskede støyen i det mottatte signalet kan hydrofonene 21 være adskilt med intervaller på 3,125 meter. Selv om eksperimentelle data tyder på at bedre støydempning kan oppnås med samplingsintervaller på 2,25 meter, er optimal avstand avhengig av en rekke begrensninger, slik som tilgjengelig båndbredde for datatransmisjon og registrering eller produksjonskostnader. De ovennevnte verdier for hydrofonavstand er utledet fra en våt (oljefylt - kerosene filled) streamer kabel og modifikasjon av samplingsintervallet kan være nødvendig for andre typer streamere, slik som faste og halv-faste streamere.
I kombinasjon med egnede støyfiltreringsmetoder, som for eksempel diskutert under, resulterer den samplingsavstanden som er beskrevet over, i reduksjon av uønsket støy, spesielt koherent støy, slik som utbulingsbølgestøy, dønningstøy og tverrstrømsstøy.
Det vises nå til figur 2, hvor det er vist et generelt blokkdiagram for en adaptiv bølgeformer brukt som filter for å redusere støyen til de enkelte sensorregistreringer. Det antas at det foreligger K sensorer plassert ved rk, med k = 1,...,K. Hver sensor k registrerer signal gk(n) med n = 1,...,N, ved bruk av en A/D-omformer. bokstaven 'n' benyttes som en indeks for diskrete tidssampler. Samplingsintervallet er At. Signalene gk(n) er strålestyrte ved bruk av forsinkelser Tk mot en generell "signalretning". Dette er den generelle retningen fra hvilken de seismiske signalene antas å ankomme. De strålestyrte kanaleneXk(n) behandles av lokale multikanals adaptive filtre for å produsere utgangssignalet:
hvorWikv(t)er de justerbare koeffisientene til det adaptive filteret, h;(n) er vinduene som anvendes ved utgangsenden, M er antallet lokale multikanals adaptive filtre (eller antallet utgangsvinduer), og L = Li + L2+1 er antallet koeffisienter per kanal. Her og under angir en understrekning under en bokstav en vektor (liten bokstav) eller en matrise (stor bokstav).
Ligning [1] kan omskrives som en (vindus-) sum over et skalarprodukt ved bruk av en tapp-inngangsvektor x(n) ved tiden t, definert som:
og en tappe-vekt vektorWjdefinert som Ved å benytte definisjonene [2] og [3], blir ligning [1] til
Ligningene [1] og [4] beskriver hvordan man finner utgangssignalet fra stråleformeren eller filterbanken så snart de M tappevektene Wi har blitt angitt. Disse vektorene beregnes som løsningen på et optimeringsproblem som er beskrevet under.
Optimeringsproblemet er definert som
med begrensningene
KL er det totale antall filterkoeffi si enter, og ||.|| angir L2normen, denne kostnadsfunksjonen er en lineær kombinasjon av utgangseffekten til stråleformeren (det første leddet i ligning [5]), og den såkalte "hvit-støy forsterkningen" til strål eformeren veid med inngangseffekten (det andre leddet i ligning [5]). Den relative vekten til de to leddene justeres av 8<2>termen. Inkluderingen av "hvit-støy forsterkningen" til strål eformeren i kostnadsfunksjonen har som hensikt å øke stråleformerens robusthet ved usikkerheter i modellering av signalet (av og til referert til som perturbasjoner) og numerisk korrelasjon mellom signalet og støyen.
Ligning [6] beskriver Q lineære begrensninger på de tillatte løsningene på optimeringsproblemet. Her er KLxQ-matrisen C begrensningsmatrisen, og Q-vektoren f er responsvektoren. Den faktiske utformingen av de lineære begrensningene diskuteres under.
En mulig løsning på optimeringsproblemet avhenger av å påføre de følgende to begrensninger på vindusfunksjonen h;(n):
for j <> i-l, i, i+1. Den første begrensningen sikrer at filterbanken er ekvivalent med enkeltfiltertilfellet dersom alle de lokale filtrene (wj) er identiske. Den andre begrensningen sikrer at vinduene har kompakt støtte. Optimeringsproblemet kan i stor grad løsrives ved bruk av den andre betingelsen (ligning [10]), og approksimasjonen
Approksimasjonen i ligning [11] krever at filtre som ligger ved siden av hverandre produserer lignende resultater når de benyttes på samme inngangsdata i tidsområder hvor tilstøtende vinduer overlapper hverandre, isteden for å kreve at tilstøtende filtre er like på en punkt-til-punkt basis. På denne måten er approksimasjonen tilsvarende det å kreve at integralet av to funksjoner ligger nær hverandre heller enn selve funksjonene.
Med denne approksimasjonen blir det første leddet i kostnadsfunksjonen, Ji
Det andre leddet i kostnadsfunksjonen kan omskrives som:
hvor "tr" angir trasen til matrisen.
Ved å kombinere ligningene (5), (12), (14) og reorganisering av leddene, kan den totale kostnadsfunksjonen skrives som:
hvor I angir KLxKL-enhetsmatrisen. Det avkoblede optimeringsproblemet kan løses for hver av de M tidsvinduene under begrensningene [6]. Ved å benytte metoden med Lagrangemultiplikatorer, er de optimale tappevektene i hvert vindu gitt av med
Det andre leddet i den modifiserte lokale korrelasjonsmatrisen Oj~kan anses som et regulariseringsledd med 8<2>som den regulariserende parameter. I kjede-signalbehandlingslitteraturen har regularisering av korrelasjonsmatriser med addisjon av en skalert enhetsmatrise blitt foreslått for å øke robustheten ved tilstedeværelse av perturbasjoner, i sammenheng med smalbånds-stråleforming. Her omfatter kostnadsfunksjonen [5] regulariseringsleddet fra begynnelsen, noe som fører til en generalisering for bredbånds adaptiv stråleforming. Dermed endres filterresponsen som en funksjon av frekvensen til signalet.
Når inngangsdataene til stråleformeren karakteriseres av romlig og tidsmessig ukorrelert (eller hvit) støy, blir både korrelasjonsmatrisen Oj og den modifiserte korrelasjonsmatrisen Oj~proporsjonal med enhetsmatrisen. I dette tilfelle blir den optimale vekt-vektoren
Vekt-vektoren wq kalles den passive (quiescent) løsningen til det optimale stråleformer-problemet, og den tilhørende responsen er kjent som den passive responsen. Legg merke til at den passive løsningen avhenger fullstendig av begrensningsmatrisen C og responsvektoren f.
Den optimale vekt-vektoren w, går mot den passive vekt-vektorenWq selv for generelle støyfelt ettersom regulariseringsparameteren 82 øker. I dette tilfelle går den modifiserte korrelasjonsmatrisen Oj~mot enhetsmatrisen (jf. [17]). Regulariseringsparameteren 8<2>vekter derfor den optimale løsningen mellom en løsning som er fullstendig avhengig av de mottatte data og en løsning som er uavhengig av dataene. For 8<2>= 1 har begge løsningene samme vekt i den forstand at deres korresponderende korrelasjonsmatriser har samme traseverdi. I situasjoner hvor perturbasjonen er høyere, dvs at antagelsene om den seismiske innsamlingsgeometrien ikke er nøyaktige, kan det å finne en stråleformerrespons med en høyere grad av regularisering gi mer robuste resultater.
En annen side ved oppfinnelsen knytter seg til utformingen av lineære begrensninger (ligning [6]) som skal påtrykkes stråleformeren.
En type lineære begrensninger som kan påtrykkes stråleformeren er de som er utformet for å bevare seismiske signaler som er innfallende fra en målretning, mens interferens som faller inn fra en annen retning undertrykkes. Styringsforsinkelser tk, slik som de vist i figur 2, definerer en enkelt "synsretning". Signaler som faller inn i denne retningen er i fase, og for disse signalene kan systemet anses å være et enkelt FIR (finite impulse response) filter. Verdiene til koeffisientene for denne ekvivalente prosessoren er lik med summen av de korresponderende koeffisientene i den adaptive prosessoren. Hver lokal stråleformer w, består av de adaptive filtreneWji,Wi2, ..,WiK, som behandler data fra hver kanal, og en summeringsenhet. Summen av de individuelle filtrene wu,Wi2,WiKer begrenset til å giWgq, som er den ønskede responsen for signaler som faller inn i synsretningen, f.eks. en enhetspuls i synsretningen:
for i = 1, ..., M ogWikoppdelt ifølge De lineære begrensningsligningene [6] kan omskrives som matriseligningen hvor KLxL-matrisen
er begrensningsmatrisen, og I er LxL-enhetsmatrisen.
For en forvrengningsfri respons i synsretningen, kan w^q velges som en enhetsimpuls, f.eks.
Den passive responsen blir da den for en fast-vektet stråleformer med enkle like vekter for alle elementer. I frekvens-bølgetall-domenet tilsvarer dette en sync-funksjon som er konstant i f-retningen. Derfor bevarer stråleformeren, for økende verdier av regulariseringsparameteren 8<2>, signaler som er innfallende ikke bare fra synsretningen, men også fra tilgrensende retninger.
Selv om det ikke på noen måte må tolkes begrensende, kan eksempler på behandlingsparametre inkludere bruk av syv inngangstraser med 3,125 meter mellomrom, en filterlengde på 11 punkter, en 83 punkters vinduslengde, 16 millisekunder resampling ved bruk av kvadraturspeil-filterteknikk, standard FK dekomposisjon for høyere delbånd, og en 0.2 regulariseringsfaktor.
Fremgangsmåten og systemet ifølge oppfinnelsen kan implementeres ved hjelp av en lang rekke alternative utførelsesformer. Hydrofonene kan være distribuert langs streameren med regelmessige eller uregelmessige intervaller, og streamerene kan være våte, faste eller halv-faste. Analog-til-digitalomformermidlene er fortrinnsvis plassert nær hydrofonene, men de kan alternativt plasseres fjernere fra disse, eller de kan til og med være plassert i det seismiske undersøkelsesfartøy et. Filtrerte utgangssignaler kan genereres og registreres for hver av hydrofonene. Dette vil typisk gjøres ved filtrering av de digitaliserte signalene som er produsert av hver hydrofon sammen med de digitale signalene produsert av de nærmeste N hydrofonene på hver side av hydrofonen (hvor N er et heltall større enn eller lik 1, slik som 3). Alternativt kan filtreringsprosessen benyttes i forbindelse med digital gruppedannelse. Ved denne typen utførelse kan for eksempel hydrofonavstanden være 3,125 meter og den romlige separasjonen av utgangstrasene kan være 12,5 meter. Ved denne typen utførelse vil filtreringen typisk benytte det digitaliserte signalet fra hydrofonen plassert ved posisjonen til den ønskede utgangstrase med de digitaliserte signalene fra mellom 1 og 4 av de nærmeste hydrofonene på hver side av denne hydrofonen. Mens filtreringen typisk vil utnytte et odde antall digitaliserte signaler (tre eller flere) når hydrofonavstanden er regulær (dvs produserer en utgangstrase for posisjonen til senterhydrofonen), kan filtreringsprosessen også benyttes for å generere en "syntetisk" trase plassert ved senter av et like antall digitaliserte signaler.
Som diskutert i den siste seksjonen, er det mulig, ved bruk av synsretnings-begrensninger og regulari sering, å bevare signaler som er innfallende fra retninger nær synsretningen.

Claims (3)

1. Fremgangsmåte for å behandle marine seismiske data innsamlet ved tauing av i det minste en seismisk streamer med et flertall av i hovedsak uniformt fordelte hydrofoner distribuert ved intervaller på ikke mer enn 500cm langs denne, i vannet over området som skal undersøkes, å rette akustiske signaler ned gjennom vannet og inn i grunnen under vannet, og å motta, med hydrofonene, seismiske signaler reflektert fra lag i grunnen under vannet, karakterisert ved å digitalisere et utgangssignal fra hver hydrofon separat, og å filtrere dette utgangssignalet for å redusere støyen som er tilstede i nevnte utgangssignal og å generere et signal med et redusert støyinnhold, idet nevnte filtreringstrinn benytter, som et ytterligere inngangssignal, det digitale utgangssignalet til i det minste én nærliggende hydrofon.
2. Fremgangsmåte ifølge krav 1, karakterisert ved at nevnte filtrering vesentlig demper dønningsstøy som er til stede i nevnte utgangssignal.
3. Fremgangsmåte ifølge krav 1, karakterisert ved at nevnte filtrering vesentlig demper tverrstrømsstøy som er til stede i nevnte utgangssignal.
NO20131614A 1998-05-20 2013-12-05 Fremgangsmåte for å redusere støy i marine seismiske registreringer NO335663B1 (no)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
GBGB9810706.3A GB9810706D0 (en) 1998-05-20 1998-05-20 Marine seismic acquisition system and method
PCT/GB1999/001544 WO1999060421A1 (en) 1998-05-20 1999-05-14 Marine seismic acquisition system and method

Publications (2)

Publication Number Publication Date
NO20131614L true NO20131614L (no) 2001-01-19
NO335663B1 NO335663B1 (no) 2015-01-19

Family

ID=10832298

Family Applications (2)

Application Number Title Priority Date Filing Date
NO20005865A NO335023B1 (no) 1998-05-20 2000-11-20 System og fremgangsmåte for å redusere støy i marine seismiske registreringer
NO20131614A NO335663B1 (no) 1998-05-20 2013-12-05 Fremgangsmåte for å redusere støy i marine seismiske registreringer

Family Applications Before (1)

Application Number Title Priority Date Filing Date
NO20005865A NO335023B1 (no) 1998-05-20 2000-11-20 System og fremgangsmåte for å redusere støy i marine seismiske registreringer

Country Status (14)

Country Link
US (1) US6684160B1 (no)
EP (1) EP1084426B1 (no)
CN (1) CN1203323C (no)
AU (1) AU758743B2 (no)
BR (1) BR9911036A (no)
CA (1) CA2327365A1 (no)
DE (1) DE69926631D1 (no)
EA (1) EA002612B1 (no)
GB (2) GB9810706D0 (no)
ID (1) ID27045A (no)
NO (2) NO335023B1 (no)
NZ (1) NZ507839A (no)
OA (1) OA12191A (no)
WO (1) WO1999060421A1 (no)

Families Citing this family (64)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9906456D0 (en) * 1999-03-22 1999-05-12 Geco Prakla Uk Ltd Method and system for reducing effects of sea surface ghost contamination in seismic data
GB0003593D0 (en) * 2000-02-17 2000-04-05 Geco As Marine seismic surveying
GB2396013B (en) 2002-12-04 2006-03-08 Westerngeco Seismic Holdings Processing seismic data
GB2410551B (en) * 2004-01-30 2006-06-14 Westerngeco Ltd Marine seismic acquisition system
US7359283B2 (en) * 2004-03-03 2008-04-15 Pgs Americas, Inc. System for combining signals of pressure sensors and particle motion sensors in marine seismic streamers
US20050194201A1 (en) * 2004-03-03 2005-09-08 Tenghamn Stig R.L. Particle motion sensor for marine seismic sensor streamers
AU2004317795B2 (en) 2004-03-17 2008-08-28 Westerngeco Seismic Holdings Ltd. Marine seismic survey method and system
AU2004319619B2 (en) 2004-05-04 2009-10-01 Westerngeco Seismic Holdings Limited Enhancing the acquisition and processing of low frequencies for sub-salt imaging
US7573781B2 (en) 2004-07-30 2009-08-11 Teledyne Technologies Incorporation Streamer cable with enhanced properties
CN100339723C (zh) * 2004-12-30 2007-09-26 中国科学院地质与地球物理研究所 高频海底数字地震仪
US8477561B2 (en) 2005-04-26 2013-07-02 Westerngeco L.L.C. Seismic streamer system and method
GB2428089B (en) * 2005-07-05 2008-11-05 Schlumberger Holdings Borehole seismic acquisition system using pressure gradient sensors
US7379391B2 (en) * 2005-11-18 2008-05-27 Westerngeco L.L.C. Marine seismic air gun timing
US7400552B2 (en) * 2006-01-19 2008-07-15 Westerngeco L.L.C. Methods and systems for efficiently acquiring towed streamer seismic surveys
US7426439B2 (en) * 2006-05-11 2008-09-16 Westerngeco L.L.C. Method and apparatus for marine seismic data acquisition
US7466625B2 (en) * 2006-06-23 2008-12-16 Westerngeco L.L.C. Noise estimation in a vector sensing streamer
US7379386B2 (en) * 2006-07-12 2008-05-27 Westerngeco L.L.C. Workflow for processing streamer seismic data
RU2321023C1 (ru) * 2006-07-17 2008-03-27 Калининградский Пограничный Институт Фсб России Сейсмический компенсатор помех блока увч рыбопоискового эхолота
US7492665B2 (en) 2006-07-24 2009-02-17 Westerngeco L.L.C. Seismic data acquisition and source-side derivatives generation and application
US7298672B1 (en) * 2006-08-22 2007-11-20 Pgs Geophysical Marine seismic streamer having acoustic isolation between strength members and sensor mounting
US8593907B2 (en) * 2007-03-08 2013-11-26 Westerngeco L.L.C. Technique and system to cancel noise from measurements obtained from a multi-component streamer
US8559265B2 (en) 2007-05-17 2013-10-15 Westerngeco L.L.C. Methods for efficiently acquiring wide-azimuth towed streamer seismic data
US8488409B2 (en) * 2007-05-17 2013-07-16 Westerngeco L.L.C. Acquiring azimuth rich seismic data in the marine environment using a regular sparse pattern of continuously curved sail lines
US7616522B2 (en) * 2007-05-18 2009-11-10 Input/Output, Inc. Seismic streamer with irregularly spaced hydrophones
US8136383B2 (en) * 2007-08-28 2012-03-20 Westerngeco L.L.C. Calibrating an accelerometer
US9158015B2 (en) * 2007-10-04 2015-10-13 Westerngeco L.L.C. Seismic streamer platform
US8553490B2 (en) * 2007-11-09 2013-10-08 Pgs Geophysical As Array grouping of seismic sensors in a marine streamer for optimum noise attenuation
US9001618B2 (en) * 2007-12-05 2015-04-07 Pgs Geophysical As Method of attenuating noise in marine seismic streamers utilizing varied sensor spacing and position-dependent band-pass filters
US8730761B2 (en) * 2007-12-17 2014-05-20 Westerngeco L.L.C. Attenuating noise in seismic data
US8781749B2 (en) 2007-12-28 2014-07-15 Westerngeco L.L.C. Attenuating noise in seismic data
US8681580B2 (en) 2008-05-15 2014-03-25 Westerngeco L.L.C. Multi-vessel coil shooting acquisition
US9857491B2 (en) 2008-05-15 2018-01-02 Westerngeco L.L.C. Multi-vessel coil shooting acquisition
US8724426B2 (en) 2008-06-03 2014-05-13 Westerngeco L.L.C. Marine seismic streamer system configurations, systems, and methods for non-linear seismic survey navigation
CN101604028B (zh) * 2008-06-11 2011-09-07 中国石油天然气集团公司 一种利用面炮合成记录确定海底高程的方法
US9594181B2 (en) 2008-06-13 2017-03-14 Westerngeco L.L.C. Filtering and presentation of heading observations for coil shooting
US9052411B2 (en) 2008-06-13 2015-06-09 Westerngeco L.L.C. Method to determine the deviation of seismic equipment from a planned curved path
US8687461B2 (en) * 2008-09-02 2014-04-01 Westerngeco L.L.C. Marine seismic source handling system
US7974151B2 (en) * 2008-09-17 2011-07-05 Westerngeco L.L.C. Cetacean protection system
US8483008B2 (en) * 2008-11-08 2013-07-09 Westerngeco L.L.C. Coil shooting mode
US20110158048A1 (en) * 2009-12-28 2011-06-30 Guigne Jacques Y Spiral sensor configuration for seismic beamforming and focusing
US8711654B2 (en) * 2009-12-30 2014-04-29 Westerngeco L.L.C. Random sampling for geophysical acquisitions
US8681581B2 (en) * 2009-12-30 2014-03-25 Westerngeco L.L.C. Randomization of data acquisition in marine seismic and electromagnetic acquisition
US8588025B2 (en) * 2009-12-30 2013-11-19 Westerngeco L.L.C. Method and apparatus for acquiring wide-azimuth marine data using simultaneous shooting
US8737163B2 (en) * 2010-02-17 2014-05-27 Westerngeco L.L.C. Wide seismic source systems
US8543342B1 (en) * 2010-09-30 2013-09-24 The United States Of America As Represented By The Secretary Of The Navy Towed array flow noise test apparatus
US8634270B2 (en) 2010-10-01 2014-01-21 Westerngeco L.L.C. Determining sea conditions in marine seismic spreads
US8947973B2 (en) 2010-11-17 2015-02-03 WesternGeco L.L.P. Active detection of marine mammals during seismic surveying
US8614580B2 (en) 2010-12-13 2013-12-24 Westerngeco L.L.C. Dynamically activating different subsets of a plurality of electrodes
US8963549B2 (en) 2010-12-13 2015-02-24 Westerngeco L.L.C. Electromagnetic measurements using a group of two or more electromagnetic receivers
US8949030B2 (en) 2011-07-29 2015-02-03 Westerngeco L.L.C. Attenuating sea-surface ghost wave effects in seismic data
US10139505B2 (en) 2011-08-09 2018-11-27 Pgs Geophysical As Digital sensor streamers and applications thereof
US8650963B2 (en) 2011-08-15 2014-02-18 Pgs Geophysical As Electrostatically coupled pressure sensor
US8717845B2 (en) 2011-08-24 2014-05-06 Pgs Geophysical As Quality-based steering methods and systems for 4D geophysical surveys
US9103942B2 (en) 2011-10-28 2015-08-11 Westerngeco L.L.C. Methods and systems for survey designs
US9274239B2 (en) 2012-01-13 2016-03-01 Westerngeco L.L.C. Wavefield deghosting
US20130265849A1 (en) 2012-04-04 2013-10-10 Westerngeco L.L.C. Methods and devices for enhanced survey data collection
US20140078860A1 (en) * 2012-09-19 2014-03-20 Cgg Services Sa Interference noise attenuation method and apparatus
CN102955172B (zh) * 2012-10-11 2015-12-16 中国水电顾问集团贵阳勘测设计研究院 水上走航式地震勘探方法及装置
US9423522B2 (en) 2012-12-11 2016-08-23 Westerngeco L.L.C. Communication systems for water vehicles
US20140241119A1 (en) * 2013-02-22 2014-08-28 Cgg Services Sa Far-field detection device, system and method
US9423519B2 (en) 2013-03-14 2016-08-23 Pgs Geophysical As Automated lateral control of seismic streamers
CN109001794A (zh) * 2018-06-15 2018-12-14 山东大学 适用于城市环境的车载移动式地震探测系统和方法
US11035969B2 (en) 2019-06-19 2021-06-15 Magseis Ff Llc Marine diffraction survey for small object detection
CN111694057B (zh) * 2020-06-03 2021-03-23 西安交通大学 一种压制地震资料涌浪噪声的方法、存储介质及设备

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4310904A (en) * 1974-09-30 1982-01-12 Sanders Associates, Inc. Ambient sea noise elimination method and apparatus
US4319347A (en) * 1976-03-08 1982-03-09 Western Geophysical Co. Of America Seismic method and system of improved resolution and discrimination
US4254480A (en) * 1978-09-11 1981-03-03 Standard Oil Company (Indiana) Frequency independent directionally sensitive array in seismic surveying
US4744064A (en) * 1984-07-20 1988-05-10 Exxon Production Research Company Method and apparatus for seismic exploration
GB2180341B (en) * 1985-09-13 1989-11-15 Amoco Corp Method of acquiring and processing seismic energy signals and marine seismic array
US4821241A (en) * 1988-05-23 1989-04-11 Teledyne Exploration Co. Noise-cancelling streamer cable
NO303305B1 (no) * 1991-02-01 1998-06-22 Geco As FremgangsmÕte ved omkobling og sammenkobling av sensorgrupper
US5251183A (en) * 1992-07-08 1993-10-05 Mcconnell Joseph R Apparatus and method for marine seismic surveying utilizing adaptive signal processing
US5448531A (en) * 1994-05-05 1995-09-05 Western Atlas International Method for attenuating coherent noise in marine seismic data
GB2309082B (en) * 1996-01-09 1999-12-01 Geco As Noise filtering method
GB2337591B (en) * 1998-05-20 2000-07-12 Geco As Adaptive seismic noise and interference attenuation method

Also Published As

Publication number Publication date
AU758743B2 (en) 2003-03-27
GB9911155D0 (en) 1999-07-14
EA200001211A1 (ru) 2001-06-25
ID27045A (id) 2001-02-22
NO20005865L (no) 2001-01-19
WO1999060421A1 (en) 1999-11-25
EP1084426A1 (en) 2001-03-21
CN1301350A (zh) 2001-06-27
EA002612B1 (ru) 2002-06-27
GB2337592B (en) 2000-05-03
GB2337592A (en) 1999-11-24
NO335663B1 (no) 2015-01-19
BR9911036A (pt) 2001-02-06
NO20005865D0 (no) 2000-11-20
CN1203323C (zh) 2005-05-25
GB9810706D0 (en) 1998-07-15
US6684160B1 (en) 2004-01-27
OA12191A (en) 2006-05-09
NZ507839A (en) 2003-11-28
NO335023B1 (no) 2014-08-25
DE69926631D1 (de) 2005-09-15
EP1084426B1 (en) 2005-08-10
AU3942399A (en) 1999-12-06
CA2327365A1 (en) 1999-11-25

Similar Documents

Publication Publication Date Title
NO20131614L (no) Fremgangsmåte for marine seismiske undersøkelser
AU2005208086B2 (en) Marine seismic acquisition system
AU2008249205B2 (en) Method and apparatus for marine seismic data acquisition
EP2060934B1 (en) Array grouping of seismic sensors in a marine streamer for optimum noise attenuation
MXPA05002448A (es) Sistema para combinar senales de sensores de presion y sensores del movimiento de particulas en streamers sismicos marinos.
NO334323B1 (no) Apparat og fremgangsmåte for innsamling av multikomponente geofysiske data
MXPA06004595A (es) Metodo y sistema de cable de registro sismico.
NO339003B1 (no) Fremgangsmåte for demping av støy i marinseismiske lyttekabler
US7167413B1 (en) Towed streamer deghosting
NO341505B1 (no) Demping av støy i seismiske streamere ved variert sensoravstand og posisjonsavhengige båndpassfiltre
EP3112907B1 (en) Separation of up-going and down-going wavefields including the direct arrival
GB2461418A (en) Seismic streamer including a multicomponent section
MX2012005137A (es) Sistema y tecnica para aumentar la separacion de los sensores del movimiento de particulas en su cable sismico.
MXPA06008525A (en) Marine seismic acquisition system
NO20151668L (no) Seismisk hydrofonkabelsystem

Legal Events

Date Code Title Description
MM1K Lapsed by not paying the annual fees