NO332880B1 - Fremgangsmate for a redusere effekter fra grov sjo pa seismiske data - Google Patents
Fremgangsmate for a redusere effekter fra grov sjo pa seismiske data Download PDFInfo
- Publication number
- NO332880B1 NO332880B1 NO20014595A NO20014595A NO332880B1 NO 332880 B1 NO332880 B1 NO 332880B1 NO 20014595 A NO20014595 A NO 20014595A NO 20014595 A NO20014595 A NO 20014595A NO 332880 B1 NO332880 B1 NO 332880B1
- Authority
- NO
- Norway
- Prior art keywords
- height
- sea surface
- seismic
- liquid medium
- sea
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000001603 reducing effect Effects 0.000 title claims description 4
- 230000000694 effects Effects 0.000 claims abstract description 19
- 238000012545 processing Methods 0.000 claims abstract description 14
- 230000036962 time dependent Effects 0.000 claims abstract description 14
- 239000007788 liquid Substances 0.000 claims description 15
- 238000005316 response function Methods 0.000 claims description 5
- 239000012530 fluid Substances 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 3
- 230000001902 propagating effect Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 25
- 238000012937 correction Methods 0.000 abstract description 10
- 238000004364 calculation method Methods 0.000 abstract description 4
- 238000001228 spectrum Methods 0.000 description 21
- 230000004044 response Effects 0.000 description 13
- 230000006870 function Effects 0.000 description 11
- 238000013213 extrapolation Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000010420 art technique Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 230000003334 potential effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004379 similarity theory Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3808—Seismic data acquisition, e.g. survey design
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/56—De-ghosting; Reverberation compensation
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Det er beskrevet en fremgangsmåte for å korrigere seismiske data for virkningene av grov sjø. Korreksjoner for virkningen av en grov havoverflate blir gjort ved å bestemme den tidsavhengige høyde av havoverflaten, enten ved direkte måling eller ved beregning fra de seismiske data. En dekonvolveringsoperator blir generert fra den målte eller beregnede høyde av havoverflaten og blir brukt til å redusere eller eliminere virkningene av den grove havoverflate. Når dataene er blitt korrigert for virkningen av den grove havoverflaten, kan de behandles på en hvilken som helst passende måte egnet for behandling av seismiske data frem- skaffet under forhold med rolig sjø eller plan havoverflate.
Description
Foreliggende oppfinnelse vedrører reduksjon av dobbeltrefleksjoner (spøk-elsesrefleksjoner) fra overflaten i seismiske data fremskaffet i et flytende medium. Oppfinnelsen angår spesielt en fremgangsmåte for å korrigere for virkningene av grov sjø på marine, seismiske data.
Fig. 1 er et skjematisk diagram som viser refleksjoner mellom en havoverflate (S), en havbunn (W) og en målreflektor (T). Forskjellige hendelser som vil bli registrert i seismogrammet, er vist og er merket i samsvar med den rekke med grenseflater som de er reflektert ved. Stjernene indikerer den seismiske kilde og pilhodene indikerer forplantningsretningen ved mottakeren. Hendelser som ender med "S" ble sist reflektert ved den grove havoverflate og blir kalt mottakerspøkel-seshendelser. Nedadgående spøkelsesrefleksjoner fra havoverflaten er en uøns-ket forurensningskilde som vanskeliggjør tolkningen av de oppadgående refleksjoner fra jordens undergrunn.
Fjerning av spøkelsesrefleksjonene fra seismiske data er for mange ekspe-rimentelle konfigurasjoner ekvivalent med opp/ned-bølgefeltseparasjon av de registrerte data. I slike konfigurasjoner representerer den oppadgående del av bølge-feltet spøkelsesrefleksjonene og det oppadgående bølgefelt representerer det ønskede signal.
Spøkelsesrefleksjoner fra havoverflaten vil inntreffe under alle sjøtilstander. Grov sjø er en ytterligere støykilde i seismiske data. Bortsett fra den ofte obser-verte dønningstøy blir ytterligere feil innført i refleksjonshendelsene ved spøkel-sesrefleksjon og spredning fra den grove havoverflaten. Spøkelseshendelser som er forstyrret av grov sjø innfører feil som er signifikante for tidsbestemte seismiske målinger og den pålitelige innsamling av repeterbare data for stratigrafisk in-versjon.
Virkningen av grov sjø er å forstyrre amplituden og ankomsttiden for spøkel-sesrefleksjoner fra havoverflaten og tilføye en spredningsavsluttende del, eller
hale, til spøkelsespulsen. Impulsresponsen kan beregnes ved endelige differanse-eller Kirchhoff-metoder (f.eks.) fra en spredende overflate som statistisk representerer typiske overflater for grov sjø. Retningsformen til Pierson-Moskowitz-spekteret som er beskrevet av Pierson, W. J. og Moskowitz, L, 1964 "A proposed Spec-tral Form for Fully Developed Wind Seas Based on the Similarity Theory of S. A. Kitaigorodskii" J. Geo. Res., 69, 24, 5181-5190, (heretter "Pierson og Moskowitz
(1964)"), and Hasselmann, D. E., Dunckel, M. og Ewing, J. A., 1980 "Directional Wave Spectra Observed During JONSWAP 1973", J. Phys. Oceanography, v10, 1264-1280, (heretter "Hasselmann et al, (1980)"). Både vindens hastighet og retning definerer spektrene. Den signifikante bølgehøyde (SWH, Significant Wave Height) er den subjektive bølgeamplitude fra bølgetopp til bølgebunn, og er lik omtrent 4 ganger RMS bølgehøyden. Forskjellige realiseringer blir oppnådd ved å multiplisere det todimensjonale overflatespektrum med gaussisk tilfeldige tall.
Fig. 2 viser et eksempel på impulser fra grov sjø langs en 400 meter lang todimensjonal linje (f.eks. en slepekabel) beregnet under en grov havoverflate med tredimensjonal SWH på 2 m. Slepekabelformen påvirker detaljene til impuls-ene, og i dette eksempelet er slepekabelen rett og horisontal. Fig. 2 viser, fra øverst til nederst: spøkelsesbølgens ankomsttid, spøkelsesbølgens maksimale amplitude, et snitt gjennom den grove sjøen over slepekabelen, og de beregnede impulser fra den grove sjøen. Den sorte toppen er den oppadgående bølge som er uforstyrret; den hvite bølgedalen er spøkelsesrefleksjonen som er reflektert fra den grove sjøen. Sistnevnte del av bølgen ved hver mottaker er spredningshalen fra stadig mer fjerntliggende deler av overflaten. Legg merke til at amplituden og ankomsttiden til spøkelsesforstyrrelser endres ganske langsomt med offset. Forstyrrelsene i ankomsttid blir styrt av de dominerende bølgelengder i havoverflaten, som er 100-200 m for 2-4 m SWH sjøer, og amplitudeforstyrrelsene blir styrt av krumningen til havoverflaten som har en RMS-radius på omkring 80 m og er ganske uavhengig av tilstanden i sjøen. Den diffraksjonshalen opptrer som kvasi-tilfeldig støy som følger spøkelsespulsen.
Forstyrrelsene fra den grove sjøen forårsaker en delvis fylling og en delvis
forskyvning av spøkelseshakket i frekvensdomenet. ("Spøkelseshakket" er et mini-mum i spekteret som forårsakes av destruktiv interferens mellom det direkte signal og spøkelsessignalet). De tilføyer også spekteret en liten rippel som beløper seg til 1-2 dB av feilen for vanlig havtilstander. I det etterfølgende stakkdomene omsettes dette til en feil i signalet som er omkring -20 dB for en 2 m SWH sjø.
Fig. 3 viser et eksempel på hvordan en slik feil kan være av betydning for tidsrelaterte målinger. Panelet øverst til venstre viser en etterstakket, tidsmigrert, syntetisk seismisk seksjon med endelig differanse. Det øvre midtre panel viser de samme data, men etter simulering av produksjon i oljereservoaret ved å forskyve olje/vann-kontakten med 6 m og innføre en 6 m delvis uttømmingssone over dette. Den lille differansen kan så vidt ses på den sorte delen av refleksjonen til høyre for forkastningen like under 2s toveis forplantningstid. Panelet øverst til høyre viser differansen mellom disse to seksjonene multiplisert med en faktor på 10. Dette er den ideelle seismiske respons fra tidsforløp-uregelmessigheten.
De venstre og midtre nedre paneler viser de samme seismiske seksjoner, men forstyrrelser fra grov sjø med 2 m SWH (som beskrevet ovenfor) er blitt tilføyd rådataene før behandling. Legg merke til at forskjellige grovsjø-effekter blir tilføyet hver modell for å representere de forskjellige sjøforhold ved innsamlingstids-punktet. Den differanse som er oppnådd mellom de to seksjonene, er vist på det nedre høyre panel (igjen multiplisert med en faktor lik 10). Feilene i reflektorampli-tuden og fasen (forårsaket av forstyrrelser fra den grove sjøen) innfører støy med amplitude tilsvarende den virkelige seismiske tidsforløp-respons. I stor grad blir den virkelige responsen maskert av disse forstyrrelsene fra den grove sjøen. En fremgangsmåte for å korrigere denne type feil er klart viktig i slike tilfeller, og med de økende behov for data med høy kvalitet og lavt støynivå, blir korreksjon for den grove sjøen nødvendig selv ved beskjedne sjøtilstander.
Nøyaktige filtre for opp/ned-atskillelse av bølgefeltmålinger med flere kom-ponenter i havbunnskabel-konfigurasjoner (OBC-konfigurasjoner) er blitt utledet av Amundsen og Ikelle, og er beskrevet i GB-patentsøknad nr. 9800741.2. En normal innfallstilnærmelse til spøkelsesfilmene for data innsamlet på havbunnen, ble beskrevet av Barr, F.J. i US-patent 4,979,150 utstedt 1990, med tittel "System for attenuating water-colum reflections" (heretter kalt "Barr (1990)"). For alle praktiske formål var dette tidligere beskrevet av White, J.E., i en artikkel fra 1965 med tittel "Seismic waves: radiation, transmission and attenuation", McGraw-Hill (heretter kalt "White (1965)"). Disse tidligere kjente teknikkene krever imidlertid måling av både hastighet og trykk. Dessuten korrigerer de ikke fullstendig de komplekse refleksjoner fra overflater med grov sjø, og en tilfredsstillende metode for å eliminere eller redusere virkningene av grov sjø på seismiske data er nødvendig.
Oppfinnelsen er definert i patentkravene.
Foreliggende oppfinnelse tilveiebringer en fremgangsmåte for å redusere virkningene på seismiske data av nedadgående reflektert og/eller spredt seismisk energi som forplanter seg i et flytende medium, idet fremgangsmåten omfatter:
a) å fremskaffe seismiske data ved å benytte en seismisk kilde og en seismisk mottaker anbrakt inne i det flytende medium; b) å bestemme høyden av minst en del av overflaten til det flytende medium som en funksjon av tid; og c) å behandle de seismiske data for å korrigere for variasjoner i høyden av det flytende medium på grunn av den grove overflaten, ved å bruke en prosess
som i det minste delvis er basert på et estimat av formen til overflaten ved posisjo-ner bort fra den minst ene del av overflaten, idet estimatet er basert på den be-stemte høyde av overflaten.
Foreliggende oppfinnelse tilveiebringer således en fremgangsmåte til å korrigere for virkningene av grov sjø på marine seismiske data. I motsetning til tidligere kjente teknikker som er nevnt ovenfor, krever oppfinnelsen dessuten bare en måling av høyden av havoverflaten.
Den del av havoverflaten hvis høyde skal måles, kan være lokalisert over den seismiske kilde eller over den seismiske mottaker. Høyden av havoverflaten kan måles ved bare et enkelt punkt over den seismiske kilde, eller høyden av et område av havoverflaten over kilden kan måles. Hvis mottakeren er en slepekabel, kan alternativt høyden av havoverflaten måles ved et antall punkter langs slepekabelen.
Høyden av havoverflaten kan måles direkte, f.eks. ved hjelp av egnede sensorer anordnet på den seismiske kilde eller den seismiske mottaker. Alternativt kan høyden av havoverflaten bestemmes fra de innsamlede seismiske data.
Det vises til de vedføyde tegninger, hvor:
fig. 1 viser eksempler på enkle seismiske strålebaner for hovedhendelser, og spøkelsesrefleksjoner som i det minste blir reflektert fra overflaten til den grove sjøen;
fig. 2 illustrerer den potensielle virkning av spøkelsesrefleksjon og spredning fra en grov sjøoverflate på konsistensen til den seismiske databølgeform;
fig. 3 illustrerer den potensielle virkning av den grove sjøen på kvaliteten av seismiske tidsforløp-data;
fig. 4 illustrerer skjematisk et eksempel på en dataprosessor som kan benyttes til å utføre foretrukne utførelsesformer av oppfinnelsen;
fig. 5 viser spredefunksjonen som fremskaffes ved å måle høyden av havoverflaten gjennom en sirkulær åpning;
fig. 6 viser restspredefunksjonen som oppnås fra spredefunksjonen på
fig. 5;
fig. 7 viser amplituden til spøkelsesspektrene før og etter korreksjon for den grove sjøoverflaten over kilden i henhold til en fremgangsmåte ifølge oppfinnelsen;
fig. 8 viser fasen til spøkelsesspektrene før og etter korreksjon for den grove havoverflate over kilden i henhold til en fremgangsmåte ifølge foreliggende oppfinnelse;
fig. 9 viser amplituden til spøkelsesspektrene før og etter korreksjon for den grove havoverflate over mottakeren i henhold til en fremgangsmåte ifølge foreliggende oppfinnelse; og
fig. 10 viser fasen til spøkelsesspektrene før og etter korreksjon for den grove havoverflate over mottakeren, i henhold til en fremgangsmåte ifølge foreliggende oppfinnelse.
I den foreliggende oppfinnelse blir virkningen av en grov havoverflate på de seismiske data korrigert for ved å bestemme variasjonene i høyden til en del av havoverflaten som en funksjon av tid. De tidsavhengige variasjoner i høyden av den del av havoverflaten kan måles direkte, eller de kan beregnes fra de seismiske data. De målte eller beregnede variasjoner i høyden av delen av havoverflaten, blir brukt under behandling av de seismiske data for å redusere eller eliminere virkningene av variasjonene i høyden av havoverflaten på de behandlede data. Mer formelt blir det utviklet en dekonvolveringsoperator fra den tidsvarierende høyde av vedkommende del av havoverflaten, og de seismiske data blir operert på av denne operatoren.
Foreliggende oppfinnelse gjør det mulig å fjerne virkningene av den grove sjøen, eller i det minste redusere dem, under behandlingen av de seismiske data. De data som fremskaffes under forhold med grov sjø blir derfor konvertert, med en rimelig tilnærmelse, til data som effektivt ble innsamlet under forhold med flat eller plan sjø. Den videre behandling av dataene etter korrigering for virkningene av den grove sjøen, kan innbefatte et trinn med fjerning av spøkelsesrefleksjoner fra dataene for å eliminere virkningene av refleksjoner fra havoverflaten, selv om dette ikke alltid er nødvendig.
Det er mulig å ekstrapolere fra de målte eller beregnede høyder av havoverflaten, for å anslå den tidsavhengige form av havoverflaten bort fra det område av havoverflaten hvor høydemålingene eller beregningene ble oppnådd. Det har vist seg at denne ekstrapoleringen tilveiebringer nyttige resultater opp til godt over 100 m fra det område hvor høydemålingene ble fremskaffet.
Når den tidsavhengige høyde av havoverflaten er blitt målt eller beregnet og resultatene er blitt ekstra polert for å estimere den tidsavhengige høyde av havoverflaten bort fra det område hvor målingene eller beregningene ble tatt, blir den tidsavhengige grove sjørespons-funksjon beregnet fra den tidsavhengige høyde av havoverflaten. Dette kan gjøres ved hjelp av enhver egnet fremgangsmåte slik som f.eks. Kirchoff-integrasjon. En dekonvolveringsoperator for å eliminere virkningene av den tidsavhengige høyde av havoverflaten, blir så beregnet, og blir på-ført de seismiske data for å korrigere for virkningene av den tidsavhengige høyde av havoverflaten.
I en utførelsesform hvor høyden av havoverflaten blir målt direkte, er den del av havoverflaten hvis høyde måles, fortrinnsvis lokalisert enten over den seismiske kilde eller over den seismiske mottaker. Dermed blir høyden av havoverflaten målt bare ved ett punkt eller over et lite område (hvis den målte del befinner seg over den seismiske kilde) eller hovedsakelig langs en linje med punkter (hvis den målte del befinner seg over en slepekabel). Vinkelfordelingen til den målte havoverflate vil følgelig ikke være korrekt. Det har imidlertid vist seg at selv om vinkelfordelingen til den målte havoverflate ikke er korrekt, er det mulig å utlede en responsfunksjon som er tilstrekkelig nøyaktig til å gjøre feilene som er igjen etter dekonvolvering, til mindre enn 0,5 dB i amplitude og 3,6° i fase.
Hvis den seismiske kilde er en slepekabel, krever fremgangsmåten at slepekabelen er forsynt med et antall dybdesensorer fordelt langs dens lengde for å gjøre det mulig å bestemme høyden til havoverflaten ved et antall punkter over slepekabelen over tid (en slepekabel forblir vanligvis ved hovedsakelig konstant høyde over havbunnen). En konvensjonell slepekabel er forsynt med én eller flere dybdesensorer for å overvåke slepekabelens dybde. Disse konvensjonelle dybdesensorene er imidlertid forholdsvis ufølsomme, og er bare ment å gi en indikasjon på at slepekabelen er omtrent vannrett og holdes i konstant høyde. De konvensjonelle sensorene er ikke i stand til å bestemme den tidsvarierende form av havoverflaten med den nøyaktighet som er nødvendig i forbindelse med foreliggende oppfinnelse, av to grunner.
For det første er konvensjonelle dybdesensorer på en slepekabel for mye atskilt fra hverandre til å gi en registrering av profilen til havoverflaten som er tilstrekkelig nøyaktig for fremgangsmåten ifølge denne utførelsesform av oppfinnelsen. Oppfinnelsen krever at høyden av havoverflaten blir målt ved laterale intervaller på omkring 5 m, for å fremskaffe en tilstrekkelig nøyaktig måling av høyden til havoverflaten (høyden av havoverflaten blir fortrinnsvis målt innenfor 5 til 10 m fra hver mottaker). I konvensjonelle slepekabler er imidlertid avstanden mellom sensorene meget større enn dette, vanligvis omkring 200 m.
For det annet er dybdesensorene på en konvensjonell slepekabel filtrert til for lav båndbredde til å være av nytte i foreliggende oppfinnelse. De konvensjonelle dybdesensorer som er anordnet på en slepekabel, er hydrostatiske dybdesensorer, og disse kan detektere frekvenser på opp til omkring 0,02 Hz. Det er da et dødbånd i de registrerte data mellom 0,02 Hz til omkring 3 Hz, som er den lav-este detekterte frekvens for hydrofoner. Bølger på havoverflaten opptar imidlertid frekvensbåndet fra omkring 0,05 Hz til 0,5 Hz, slik at de dybdesensorer som er anordnet på en konvensjonell slepekabel, er uegnede for å detektere profilen til havoverflaten med den nøyaktighet som kreves for å praktisere foreliggende oppfinnelse.
For å realisere en utførelsesform av foreliggende oppfinnelse hvor høyden til havoverflaten over en slepekabel blir målt, må slepekabelen være forsynt med
sensorer som er i stand til å måle høyden til havoverflaten. Hver sensor må være i stand til å måle høyden innenfor en nøyaktighet på omkring 5-10 cm, og må være følsom i bølgebåndområdet fra 0,05 Hz til 0,5 Hz. En egnet sensor er en trykksen-sor eller en akustisk ekkodybde-sensor. Sensorene er atskilt langs lengden av slepekabelen, fortrinnsvis med jevne mellomrom. Under drift måler hver sensor høy-den av havoverflaten over sensoren, enten kontinuerlig eller med diskrete intervaller, og resultatene av målingene som tas av sensorene, blir overført til en datama-skin eller en annen dataregistreringsanordning og blir lagret for bruk i databe-handlingen.
Alternativet til å foreta måling av høyden til havoverflaten ved et antall punkter langs en linje, f.eks. langs en linje over en slepekabel, er at målingene av høy-den til havoverflaten kan tas over et område over den seismiske kilde. Dette kan tenkes på som å måle høyden av havoverflaten gjennom en åpning, med oppnåd-de resultater for bare en del av havoverflaten som er "synlig" gjennom åpningen.
Hvis arealet til en seismisk kilde er lite, og en seismisk kilde kan ha et areal på omkring bare 15 kvadratmeter, kan det godt være tilstrekkelig til å måle høyden av havoverflaten ved bare ett enkelt punkt over den seismiske kilde. Dette kan tenkes på som å måle høyden av havoverflaten gjennom en punktåpning. Når høy-den av havoverflaten blir målt ved bare ett enkelt punkt, er det fremdeles mulig å utføre det ekstra poleringstrinn som er beskrevet ovenfor for å anslå den tidsavhengige høyde av havoverflaten bort fra det punkt hvor målingen ble tatt. Selv om estimatet har et dårlig utseende, har den dekonvolveringsoperatoren som utledes fra estimatet vist seg å være nøyaktig.
I én utførelsesform hvor høyden av et område av havoverflaten over den seismiske kilde blir målt, er området over hvilket høyden av havoverflaten blir målt, fortrinnsvis sirkulært, og er fortrinnsvis også sentrert over den seismiske kilde. Dette innfører symmetri og forenkler behandlingen. I prinsippet kan det imidlertid benyttes åpninger eller aperturer som har andre former. Det område over kilden hvor høydemålingene blir utført, har fortrinnsvis en dimensjon i størrelsesorden meter. En sirkel med en radius på omkring 7,5 meter kan f.eks. benyttes.
Når høyden av havoverflaten blir målt, ikke ved et enkelt punkt, men gjennom en åpning med endelige dimensjoner, blir høyden av havoverflaten i praksis ikke målt ved hvert punkt innenfor åpningen. En todimensjonal gruppe med dybdesensorer blir brukt til å måle høyden av havoverflaten ved et antall punkter innenfor åpningen, f.eks. ved 6 eller 8 punkter. Forutsatt at sensorene er tilstrekkelig nær hverandre, f.eks. med en avstand mellom tilstøtende sensor i området fra 5-10 m, er det ikke noe betydelig tap av nøyaktighet i den etterfølgende behandling av dataene. Som i ovennevnte eksempel måler hver sensor høyden av havoverflaten ved et punkt over den seismiske kilde, enten kontinuerlig eller ved diskrete intervaller, og resultatene av målingene tatt av sensorene, blir sendt til en data-maskin eller en annen dataregistreringsanordning og blir lagret for bruk ved data-behandlingen.
En seismisk kilde er vanligvis opphengt i et flyteorgan. Én hensiktsmessig måte å måle høyden til havoverflaten ved et punkt over en seismisk kilde, er å montere en sensor på flyteorganet for å registrere flyteorganets høyde. En egnet sensor for dette er en GPS-mottaker.
Hvis målingen av høyden til havoverflaten blir foretatt gjennom en sirkulær åpning sentrert over den seismiske kilde, så vil denne åpningen i x-y-domenet (det antas at havoverflaten strekker seg i x-retningen og y-retningen, mens bølgehøy-den strekker seg i z-retningen), bli en konvolvering i kx-ky med en punktsprede-funksjon (PSF) som er en sylindrisk Bessel-funksjon, som vist på fig. 5. PSF vil maskere mye av de fine detaljene i dette havoverflate-spekteret, og dette er ekvivalent med å si at åpningen som høydemålingene ble gjort gjennom, er utilstrekke-lig til å beskrive formen av havoverflaten fullstendig. Den første nullgjennomgang for PSF er ved omkring kB = 0,5 radianer/meter, hvor kB er den radiale variable av Bessel-funksjonen. Hver bølgespektrum-sampel vil bli spredt over et område med radius kB.
Høyden til havoverflaten blir målt bare for et diskret område av havoverflaten over den seismiske kilde, og som bemerket ovenfor, fører dette til et oppløs-ningstap i spekteret. Ifølge foreliggende oppfinnelse blir imidlertid høyden til havoverflaten målt over tid, og denne målingen over tid gir ekstra informasjon som kompenserer for oppløsningstapet forårsaket av den lille åpning som høydedata-ene innsamles gjennom.
Hvis de registrerte høydedata ble filtrert slik at bare én vinkelfrekvens, w, blir betraktet ved et bestemt tidspunkt, bør det virkelige spektrum ligge fullstendig på en sirkel i bølgenummerdomenet med en radius k gitt av:
I ligning (1) er c(w) hastigheten til overflatebølger som har en vinkelfrekvens w.
PSF vil spre hver sampel på denne sirkelen over et nærliggende område
med radius kB. Den nøyaktige form av vinkelspredningsfunksjonen er k-avhengig, og har form av en fullstendig sirkel med lave bølgetall, og har form av en bue ved høyere bølgetall. Dette svarer til at mangelen på retningsmessig informasjon hovedsakelig gjelder ved lange overflatebølgelengder.
Det er enkelt å dekonvolvere den radiale spredning av en enkelt frekvens-komponent i dataene. Verdiene blir ganske enkelt flyttet radialt bakover på den sirkel som har radius k = c(w) / w. Denne prosedyren kan gjentas for hver tidsmessig frekvens i spekteret for havoverflaten. Restspredningsfunksjonen for en sampel ved (kx, ky) er således skjæringsbuen mellom en sirkel med radius k sentrert på origo, med et sirkulært område med radius kB sentrert på (kx, ky). Et eksempel på den resulterende spredefunksjon er vist på fig. 6.
Ved bølgelengder som er meget korte sammenlignet med dimensjonene til den åpning som høydemålingene ble tatt gjennom, er åpningen stor nok til at retningen av bølgene blir anslått med rimelig nøyaktighet. Det er imidlertid meget lite energi ved disse korte bølgelengdene. Mesteparten av energien til overflatebølgen er ved bølgelengder som er lenger enn dimensjonene av den åpning som høyde-målingene foretas gjennom, og overflatebølger som har disse bølgelengdene, kan ikke måles direkte over en hel bølgelengde. I stedet blir disse bølgene estimert fra tidsavhengigheten til høyden av havoverflaten. Dette gir et godt estimat på amplituden til bølgen, men et dårlig estimat om dens retning. Når bølgene rekonstrue-res, betyr denne dårlige estimeringen av retningen til bølger med lang bølgelengde at den rekonstruerte bølgeoverflate vanligvis har sirkulær symmetri. Bare ved det egentlige sentrum for den rekonstruerte overflate, hvor målingene ble tatt, opptrer virkelige reelle detaljer.
Det kan innledningsvis synes at refleksjonsresponsen som er beregnet fra en rekonstruert bølgeoverflate som har slik dårlig vinkelmessig oppløsning, vil ha liten, om noen, relasjon til den virkelige respons. Oppfinnerne har imidlertid funnet at det i virkeligheten er mulig å utlede en god refleksjonsrespons fra den rekonstruerte bølgeformoverflate. Dette er fordi amplituden og fasen til bølgene i den estimerte bølgeformoverflate er korrekte, fordi de ble utledet fra høydemålinger tatt over tid. Den vertikale responsen ved en tid t medfører et integral over overflaten ved en radius t/c (hvor c er lydhastigheten i vann), slik at retningsfeilen i spekteret for overflatebølgene er mindre signifikant. Hvis responsen var rent lineær i høy-den, så ville responsen i virkeligheten være helt korrekt. Responsen som ligger utenfor vertikalen, vil være mer følsom for retningsfeilen i den rekonstruerte overflatebølgeform, fordi integrasjonsbølgebanene ikke er sirkler. Ved typiske strålingsvinkler som benyttes ved seismiske undersøkelser, er imidlertid feilen aksepterbar. Fig. 7 til 10 viser eksempler på hvor effektiv fremgangsmåten ifølge foreliggende oppfinnelse er når det gjelder å redusere 95% konfidens-feilsøylene på spektrene for mottaker- og kilde-spøkelsesrefleksjoner for seismiske data fremskaffet under forhold med grov sjø. Fig. 7 viser resultatene av korrigering for en overflate med grov sjø ved å ta høydemålinger over kilden. Trase (a) viser amplituden til spøkelsesspekteret som er fremskaffet for en grov havoverflate med en 2 m SWH, og trase (b) viser resultatet av å korrigere trasen (a) for den grove havoverflate ved å bruke den ovenfor beskrevne fremgangsmåte. Trase (c) svarer til trase (a), men ble tilveiebrakt for en overflate med grov sjø som har 4 m SWH og trase (d) viser resultatet av å korrigere trase (c) for den grove havoverflaten ved å bruke den ovenfor beskrevne fremgangsmåte. Fig. 8 svarer til fig. 7, men viser effekten av å korrigere for den grove havoverflate på fasen til spøkelsesspekteret. Trasene (a) og (b) ble fremskaffet for sjø med to meter SWH og trasene (c) og (d) ble fremskaffet for sjø med 4 m SWH. Fig. 9 og 10 viser resultatene av en korreksjon basert på å ta høydemålin-ger over mottakeren, for amplituden (fig. 9) og fasen (fig. 10) til spøkelsesspekte-ret. I hver av figurene 9 og 10 viser trase (a) data fremskaffet i sjø med 4 m SWH, og trase (b) viser resultatet av korreksjon for den grove havoverflate ved å bruke fremgangsmåten som er beskrevet ovenfor.
I de ovenfor beskrevne utførelsesformer har høyden til havoverflaten blitt målt direkte, ved hjelp av egnede dybdesensorer anordnet på kilden eller på mottakeren. Ifølge én modifikasjon av denne utførelsesformen av oppfinnelsen blir høyden av havoverflaten ikke målt direkte, men blir anslått fra de innsamlede seismiske data. Selv om spekteret til høyden av havoverflaten ligger i dødbåndet til det seismiske datainnsamlingssystem, som nevnt ovenfor, er det likevel innfanget indi-rekte. Dette er fordi visse egenskaper ved de seismiske data (slik som f.eks. spøk-elseshakk-frekvensen), som varierer med tid over trasen, har en lav frekvens-grense som er lavere enn den lave frekvensgrensen til de innsamlede data selv. Dette skjer fordi egenskapene ikke er lineære funksjoner av dataene. Hvis f.eks. en statistisk metode blir utviklet for fjerning av spøkelsesrefleksjoner på grunn av grov sjø, hvor tidsforskyvningene for direkte- og spøkelses-hendelser blir estimert fra dataene (ved å bruke tidligere informasjon om korrelasjoner utledet fra sjøtilstanden), kan disse data benyttes som inngang til overflaterekonstruksjons- og de-konvolverings-metoden som er beskrevet ovenfor.
I den utførelsesform som er beskrevet ovenfor, er trinnene med å måle eller beregne høyden av en del av havoverflaten, ekstrapolering for å estimere høy-den av andre deler av havoverflaten, beregning av den tidsavhengige grove sjø-respons og bestemmelsen av dekonvolveringsoperatoren, blitt beskrevet separat for klarhetens skyld. I prinsippet er det imidlertid mulig å kombinere ett eller flere av disse trinn til et enkelt behandlingstrinn. Trinnet med å beregne den tidsavhengige grove sjørespons kan f.eks. kombineres med ekstrapoleringstrinnet.
Fremgangsmåten for å korrigere data fremskaffet med en grov sjøoverflate ved å bruke måling av høyden til havoverflaten, kan utføres ved å bruke et hvilket som helst konvensjonelt behandlingssystem for seismiske data. Behandlingen blir fortrinnsvis utført på en dataprosessor konfigurert for å behandle store mengder med data. Fig. 4 illustrerer f.eks. en mulig konfigurasjon for en slik dataprosessor. Dataprosessoren består vanligvis av én eller flere sentralenheter 350, et hoved-minne 352, kommunikasjons- eller inn/ut-moduler 354, grafiske anordninger 356, en akselerator 358 med flytende punktum og masselagringsanordninger slik som bånd og plater 360.
Selv om foretrukne utførelsesformer av oppfinnelsen er blitt beskrevet, er beskrivelsene og figurene bare illustrerende og ikke ment å begrense foreliggende oppfinnelse. Den eneste seismiske mottaker som det er referert til i beskrivelsen av de foretrukne utførelsesformer av oppfinnelsen, er f.eks. en slepekabel. Oppfinnelsen kan imidlertid anvendes på andre mottakere, slik som f.eks. en vertikal mottakergruppe.
Claims (12)
1. Fremgangsmåte for å redusere virkningene i seismiske data av nedadgående, reflektert og/eller spredt seismisk energi som forplanter seg i et flytende medium,
karakterisert ved: a) å fremskaffe seismiske data ved å benytte en seismisk kilde og en seismisk mottaker anbrakt inne i det flytende medium; b) å bestemme høyden av minst en del av overflaten til det flytende medium som en funksjon av tid; og c) å behandle de seismiske data for å korrigere for variasjoner i høyden av det flytende medium på grunn av den grove overflaten, ved å bruke en prosess som i det minste delvis er basert på et estimat av formen til overflaten ved posisjo-ner bort fra den minst ene del av overflaten, idet estimatet er basert på den be-stemte høyde av overflaten.
2. Fremgangsmåte ifølge krav 1,
karakterisert vedat trinn (c) omfatter å beregne en dekonvolveringsoperator fra resultatene av bestemmelsen av høyden av overflaten til det flytende medium, og å anvende dekonvolveringsoperatoren på de seismiske data.
3. Fremgangsmåte ifølge krav 2,
karakterisert vedat trinn (c) omfatter å beregne havoverflatens responsfunksjon fra den tidsavhengige høyde av havoverflaten, og å beregne dekonvolveringsoperatoren fra havoverflatens responsfunksjon.
4. Fremgangsmåte ifølge krav 3,
karakterisert vedå estimere den tidsavhengige høyde av en ytterligere del av havoverflaten fra resultatene av trinn (b) før beregning av havoverflatens responsfunksjon.
5. Fremgangsmåte ifølge noen av de foregående krav,
karakterisert vedat trinn (b) omfatter å måle høyden av i det minste en del av overflaten til det flytende medium.
6. Fremgangsmåte ifølge krav 5,
karakterisert vedat trinn (b) omfatter å måle høyden av overflaten til det flytende medium ved ett eller flere punkter.
7. Fremgangsmåte ifølge krav 6,
karakterisert vedat den seismiske mottaker er en slepekabel, og at trinn (b) omfatter å måle høyden av overflaten til det flytende medium ved et antall punkter langs lengden av slepekabelen.
8. Fremgangsmåte ifølge krav 6,
karakterisert vedat trinn (b) omfatter å måle høyden av overflaten til det flytende medium ved et enkelt punkt lokalisert over den seismiske kilde.
9. Fremgangsmåte ifølge krav 5,
karakterisert vedat trinn (b) omfatter å måle høyden til overflaten av det flytende medium over et areal av overflaten til det flytende medium.
10. Fremgangsmåte ifølge krav 9,
karakterisert vedat trinn (b) omfatter å måle høyden av overflaten til det flytende medium over et areal som er hovedsakelig sirkulært.
11. Fremgangsmåte ifølge krav 10,
karakterisert vedat trinn (b) omfatter å måle høyden av overflaten til det flytende medium over et areal som er hovedsakelig sentrert over den seismiske kilde.
12. Fremgangsmåte ifølge noen av kravene 1 til 5,
karakterisert vedat trinn (b) omfatter å beregne høyden av overflaten til det flytende medium som en funksjon av tid fra de seismiske data.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GBGB9906456.0A GB9906456D0 (en) | 1999-03-22 | 1999-03-22 | Method and system for reducing effects of sea surface ghost contamination in seismic data |
PCT/GB2000/001066 WO2000057206A1 (en) | 1999-03-22 | 2000-03-21 | A method of reducing effects of a rough sea surface on seismic data |
Publications (3)
Publication Number | Publication Date |
---|---|
NO20014595D0 NO20014595D0 (no) | 2001-09-21 |
NO20014595L NO20014595L (no) | 2001-11-21 |
NO332880B1 true NO332880B1 (no) | 2013-01-28 |
Family
ID=10850039
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
NO20014596A NO333154B1 (no) | 1999-03-22 | 2001-09-21 | Fremgangsmate og utstyr for a redusere signalforurensende pavirkning pa seismiske data pa grunn av falske ekkosignaler fra sjooverflaten |
NO20014595A NO332880B1 (no) | 1999-03-22 | 2001-09-21 | Fremgangsmate for a redusere effekter fra grov sjo pa seismiske data |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
NO20014596A NO333154B1 (no) | 1999-03-22 | 2001-09-21 | Fremgangsmate og utstyr for a redusere signalforurensende pavirkning pa seismiske data pa grunn av falske ekkosignaler fra sjooverflaten |
Country Status (5)
Country | Link |
---|---|
US (4) | US6775618B1 (no) |
AU (2) | AU3311100A (no) |
GB (3) | GB9906456D0 (no) |
NO (2) | NO333154B1 (no) |
WO (2) | WO2000057207A1 (no) |
Families Citing this family (97)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6477470B2 (en) * | 2000-12-01 | 2002-11-05 | Pgs Americas, Inc. | Method and system for deghosting |
GB2379741B (en) * | 2001-09-18 | 2003-11-19 | Westerngeco Ltd | Method for reducing the effect of Sea-surface ghost reflections |
FR2831962B1 (fr) | 2001-11-08 | 2004-06-25 | Geophysique Cie Gle | Procede de traitement sismique, notamment pour la compensation de birefringence sur des traces sismiques |
GB2384068B (en) * | 2002-01-11 | 2005-04-13 | Westerngeco Ltd | A method of and apparatus for processing seismic data |
GB2389183B (en) | 2002-05-28 | 2006-07-26 | Westerngeco Ltd | Processing seismic data |
US7239577B2 (en) * | 2002-08-30 | 2007-07-03 | Pgs Americas, Inc. | Apparatus and methods for multicomponent marine geophysical data gathering |
GB2400662B (en) | 2003-04-15 | 2006-08-09 | Westerngeco Seismic Holdings | Active steering for marine seismic sources |
US7123543B2 (en) * | 2003-07-16 | 2006-10-17 | Pgs Americas, Inc. | Method for seismic exploration utilizing motion sensor and pressure sensor data |
GB2405473B (en) | 2003-08-23 | 2005-10-05 | Westerngeco Ltd | Multiple attenuation method |
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 |
GB2414299B (en) | 2004-05-21 | 2006-08-09 | Westerngeco Ltd | Interpolation and extrapolation method for seismic recordings |
US7800977B2 (en) * | 2004-06-01 | 2010-09-21 | Westerngeco L.L.C. | Pre-stack combining of over/under seismic data |
GB2415040B (en) * | 2004-06-12 | 2007-03-21 | Westerngeco Ltd | Three-dimensional deghosting |
US7100438B2 (en) * | 2004-07-06 | 2006-09-05 | General Electric Company | Method and apparatus for determining a site for an offshore wind turbine |
US7336561B2 (en) * | 2004-09-07 | 2008-02-26 | Pgs Americas, Inc. | System for attenuation of water bottom multiples in seismic data recorded by pressure sensors and particle motion sensors |
US20060133202A1 (en) * | 2004-12-22 | 2006-06-22 | Tenghamn Stig R L | Motion sensors in a marine seismic streamer |
US7499373B2 (en) | 2005-02-10 | 2009-03-03 | Westerngeco L.L.C. | Apparatus and methods for seismic streamer positioning |
US7319636B2 (en) * | 2005-03-14 | 2008-01-15 | Westerngeco, L.L.C. | Calibration of pressure gradient recordings |
US7433264B2 (en) * | 2005-03-18 | 2008-10-07 | Westerngeco L.L.C. | Methods and systems for determination of vertical correction of observed reflection seismic signals |
US7450467B2 (en) * | 2005-04-08 | 2008-11-11 | Westerngeco L.L.C. | Apparatus and methods for seismic streamer positioning |
US7372769B2 (en) * | 2005-04-08 | 2008-05-13 | Western Geco L.L.C. | Method and apparatus for adaptive over/under combination of seismic data |
US8477561B2 (en) * | 2005-04-26 | 2013-07-02 | Westerngeco L.L.C. | Seismic streamer system and method |
US20060256653A1 (en) * | 2005-05-05 | 2006-11-16 | Rune Toennessen | Forward looking systems and methods for positioning marine seismic equipment |
US7768869B2 (en) * | 2005-05-05 | 2010-08-03 | Pgs Americas, Inc. | Method for deghosting and water layer multiple reflection attenuation in marine seismic data |
GB2428089B (en) | 2005-07-05 | 2008-11-05 | Schlumberger Holdings | Borehole seismic acquisition system using pressure gradient sensors |
US7660191B2 (en) * | 2005-07-12 | 2010-02-09 | Westerngeco L.L.C. | Methods and apparatus for acquisition of marine seismic data |
US7804738B2 (en) * | 2006-03-21 | 2010-09-28 | Westerngeco L.L.C. | Active steering systems and methods for marine seismic sources |
US20080008038A1 (en) * | 2006-07-07 | 2008-01-10 | Johan Olof Anders Robertsson | Method and Apparatus for Estimating a Seismic Source Signature |
US7480204B2 (en) * | 2006-07-07 | 2009-01-20 | Westerngeco L.L.C. | Seismic data processing |
US7379386B2 (en) * | 2006-07-12 | 2008-05-27 | Westerngeco L.L.C. | Workflow for processing streamer seismic data |
US7523003B2 (en) | 2006-07-12 | 2009-04-21 | Westerngeco L.L.C. | Time lapse marine seismic surveying |
US7379385B2 (en) * | 2006-07-26 | 2008-05-27 | Westerngeco L.L.C. | Processing of seismic data acquired using over/under streamers and/or over/under sources |
US7773456B2 (en) * | 2006-10-02 | 2010-08-10 | Bp Corporation North America Inc. | System and method for seismic data acquisition |
GB2444953B (en) * | 2006-12-19 | 2009-07-22 | Westerngeco Seismic Holdings | Method for obtaining an image of a subsurface by regularizing irregularly sampled seismic data |
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 |
US7426438B1 (en) | 2007-03-16 | 2008-09-16 | Westerngeco L.L.C. | Technique to provide seismic data-based products |
US7679991B2 (en) * | 2007-03-16 | 2010-03-16 | Westerngeco L. L. C. | Processing of seismic data acquired using twin over/under streamers |
US8077543B2 (en) * | 2007-04-17 | 2011-12-13 | Dirk-Jan Van Manen | Mitigation of noise in marine multicomponent seismic data through the relationship between wavefield components at the free surface |
US7676327B2 (en) | 2007-04-26 | 2010-03-09 | Westerngeco L.L.C. | Method for optimal wave field separation |
GB2450122B (en) * | 2007-06-13 | 2009-08-05 | Westerngeco Seismic Holdings | Method of representing signals |
US7715988B2 (en) * | 2007-06-13 | 2010-05-11 | Westerngeco L.L.C. | Interpolating and deghosting multi-component seismic sensor data |
US8014228B2 (en) * | 2007-08-09 | 2011-09-06 | Westerngeco, L.L.C. | Marine seismic sources and methods of use |
US20090040872A1 (en) * | 2007-08-09 | 2009-02-12 | Jahir Pabon | Removing Vibration Noise from Multicomponent Streamer Measurements |
US8116166B2 (en) * | 2007-09-10 | 2012-02-14 | Westerngeco L.L.C. | 3D deghosting of multicomponent or over / under streamer recordings using cross-line wavenumber spectra of hydrophone data |
NO20083861L (no) * | 2007-09-14 | 2009-03-16 | Geco Technology Bv | Bruk av kildehoydemalinger for a fjerne sjoforstyrrelser |
GB2454747B (en) * | 2007-11-19 | 2010-08-04 | Westerngeco Seismic Holdings | Moveout correction of seismic data |
US8175765B2 (en) | 2007-12-13 | 2012-05-08 | Westerngeco L.L.C. | Controlling movement of a vessel traveling through water during a seismic survey operation |
US7646672B2 (en) | 2008-01-18 | 2010-01-12 | Pgs Geophysical As | Method for wavefield separation in 3D dual sensor towed streamer data with aliased energy in cross-streamer direction |
US8077544B2 (en) * | 2008-03-28 | 2011-12-13 | Westerngeco L.L.C. | Dual-wavefield multiple attenuation |
US7957906B2 (en) * | 2008-05-07 | 2011-06-07 | Pgs Geophysical As | Method for attenuating low frequency noise in a dual-sensor seismic streamer |
US8964501B2 (en) * | 2008-05-25 | 2015-02-24 | Westerngeco L.L.C. | System and technique to determine high order derivatives from seismic sensor data |
US7872942B2 (en) * | 2008-10-14 | 2011-01-18 | Pgs Geophysical As | Method for imaging a sea-surface reflector from towed dual-sensor streamer data |
US8085617B2 (en) | 2008-10-31 | 2011-12-27 | Sercel Inc. | System and method for reducing the effects of ghosts from the air-water interface in marine seismic exploration |
CA2744446C (en) * | 2008-11-21 | 2015-04-14 | Exxonmobil Upstream Research Company | Free charge carrier diffusion response transducer for sensing gradients |
US8164977B2 (en) * | 2008-12-10 | 2012-04-24 | Westerngeco L.L.C. | Simulating up-going pressure wavefield |
US8699297B2 (en) * | 2009-02-13 | 2014-04-15 | Westerngeco L.L.C. | Deghosting and reconstructing a seismic wavefield |
US8554484B2 (en) * | 2009-02-13 | 2013-10-08 | Westerngeco L.L.C. | Reconstructing seismic wavefields |
US20100211320A1 (en) * | 2009-02-13 | 2010-08-19 | Massimiliano Vassallo | Reconstructing a seismic wavefield |
US20100211322A1 (en) * | 2009-02-13 | 2010-08-19 | Massimiliano Vassallo | Interpolating a pressure wavefield along an undersampled direction |
US8050139B2 (en) * | 2009-03-27 | 2011-11-01 | Westerngeco L.L.C. | System and method for towing acoustic source arrays |
GB0905261D0 (en) | 2009-03-27 | 2009-05-13 | Geco Technology Bv | Processing seismic data |
US9103934B2 (en) * | 2009-08-05 | 2015-08-11 | Westerngeco L.L.C. | Method for reducing marine source volume while maintaining image quality |
US20110044127A1 (en) * | 2009-08-19 | 2011-02-24 | Clement Kostov | Removing free-surface effects from seismic data acquired in a towed survey |
US9285493B2 (en) * | 2009-08-27 | 2016-03-15 | Pgs Geophysical As | Sensor grouping for dual sensor marine seismic streamer and method for seismic surveying |
FR2955396B1 (fr) | 2010-01-15 | 2013-03-01 | Cggveritas Services Sa | Dispositif de traitement de donnees sismiques marines |
FR2955397B1 (fr) | 2010-01-15 | 2012-03-02 | Cggveritas Services Sa | Procede et dispositif d'acquisition de donnees sismiques marines |
US10545252B2 (en) | 2010-01-15 | 2020-01-28 | Westerngeco L.L.C. | Deghosting and interpolating seismic data |
US8730766B2 (en) * | 2010-01-22 | 2014-05-20 | Ion Geophysical Corporation | Seismic system with ghost and motion rejection |
US8467265B2 (en) * | 2010-01-26 | 2013-06-18 | Westerngeco L.L.C. | Interpolation and deghosting of seismic data acquired in the presence of a rough sea |
EP2372318B1 (de) * | 2010-03-26 | 2020-03-18 | VEGA Grieshaber KG | Störechospeicherung bei Behälterrauschen |
US8902699B2 (en) * | 2010-03-30 | 2014-12-02 | Pgs Geophysical As | Method for separating up and down propagating pressure and vertical velocity fields from pressure and three-axial motion sensors in towed streamers |
AU2014201114B2 (en) * | 2010-06-10 | 2014-07-31 | Cggveritas Services Sa | Method and apparatus for deghosting seismic data |
FR2961316A1 (fr) * | 2010-06-10 | 2011-12-16 | Cggveritas Services Sa | Procede de traitement de donnees sismiques marines |
DE112011102495T5 (de) | 2010-07-26 | 2013-06-06 | Exxonmobil Upstream Research Company | Seismisches Erfassungsverfahren zur Modentrennung |
CN103168254B (zh) * | 2010-09-02 | 2017-03-29 | 离子地球物理学公司 | 多分量声波传感器和方法 |
US8634270B2 (en) * | 2010-10-01 | 2014-01-21 | Westerngeco L.L.C. | Determining sea conditions in marine seismic spreads |
AU2011232767B2 (en) | 2010-10-14 | 2014-05-08 | Cggveritas Services Sa | Method and device to acquire seismic data |
US20130028048A1 (en) * | 2011-07-25 | 2013-01-31 | Soellner Walter | Methods and apparatus for seismic imaging which accounts for sea-surface variations |
US20130088939A1 (en) * | 2011-10-10 | 2013-04-11 | Pascal Edme | Wavefield separation using a gradient sensor |
US9103943B2 (en) | 2011-11-28 | 2015-08-11 | Fugro-Geoteam As | Acquisition and processing of multi-source broadband marine seismic data |
US9835750B2 (en) | 2012-01-20 | 2017-12-05 | Cgg Services Sas | Buoy based marine seismic survey system and method |
EP2831633B1 (en) * | 2012-03-27 | 2020-03-11 | WesternGeco Seismic Holdings Limited | Interferometry-based data redatuming and/or depth imaging |
US9442209B2 (en) * | 2012-07-10 | 2016-09-13 | Pgs Geophysical As | Methods and systems for reconstruction of low frequency particle velocity wavefields and deghosting of seismic streamer data |
US9625600B2 (en) * | 2012-12-04 | 2017-04-18 | Pgs Geophysical As | Systems and methods for removal of swell noise in marine electromagnetic surveys |
US9405028B2 (en) | 2013-02-22 | 2016-08-02 | Ion Geophysical Corporation | Method and apparatus for multi-component datuming |
US9274241B2 (en) * | 2013-03-14 | 2016-03-01 | Pgs Geophysical As | Method and system for suppressing swell-induced electromagnetic noise |
MX355746B (es) | 2013-03-14 | 2018-04-27 | Ion Geophysical Corp | Dispositivos, sistemas y metodos de sensor sismico que incluyen filtracion de ruido. |
US10520623B2 (en) | 2013-05-31 | 2019-12-31 | Westerngeco L.L.C. | Methods and systems for marine survey acquisition |
US10598807B2 (en) | 2014-02-18 | 2020-03-24 | Pgs Geophysical As | Correction of sea surface state |
GB2525072B (en) * | 2014-02-18 | 2020-08-19 | Pgs Geophysical As | Correction of sea surface state |
US10107929B2 (en) * | 2014-12-18 | 2018-10-23 | Pgs Geophysical As | Methods and systems to determine ghost operators from marine seismic data |
US10436922B2 (en) * | 2015-10-05 | 2019-10-08 | Cgg Services Sas | Device and method for constrained wave-field separation |
US10324210B2 (en) * | 2016-06-30 | 2019-06-18 | Schlumberger Technology Corporation | Method and apparatus for determining rough sea topography during a seismic survey |
US10732311B2 (en) | 2016-10-19 | 2020-08-04 | Pgs Geophysical As | Method for adjusting attributes of marine seismic survey data |
GB2560578B (en) * | 2017-03-17 | 2022-06-15 | Equinor Energy As | A method of deghosting seismic data |
CN109975868B (zh) * | 2019-03-26 | 2020-03-03 | 中国科学院地质与地球物理研究所 | 一种基于泰勒展开的波动方程鬼波压制方法 |
Family Cites Families (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2757356A (en) * | 1954-01-08 | 1956-07-31 | Texas Instruments Inc | Method and apparatus for canceling reverberations in water layers |
US3747055A (en) | 1971-06-16 | 1973-07-17 | Texaco Inc | High resolution shooting with space domain filtering |
CA1136255A (en) * | 1978-08-17 | 1982-11-23 | Samuel M. Theodoulou | Body motion compensation filter with pitch and roll correction |
US4353121A (en) | 1980-07-24 | 1982-10-05 | Fairfield Industries, Inc. | High resolution, marine seismic stratigraphic system |
US4380059A (en) | 1980-08-20 | 1983-04-12 | Mobil Oil Corporation | F-K Filtering of multiple reflections from a seismic section |
US4486865A (en) * | 1980-09-02 | 1984-12-04 | Mobil Oil Corporation | Pressure and velocity detectors for seismic exploration |
US4992992A (en) | 1988-10-21 | 1991-02-12 | Western Atlas International, Inc. | Processing for seismic data from slanted cable |
EG19158A (en) * | 1989-08-25 | 1996-02-29 | Halliburton Geophys Service | System for attenuation of water-column reverberation |
US5051961A (en) * | 1989-10-26 | 1991-09-24 | Atlantic Richfield Company | Method and apparatus for seismic survey including using vertical gradient estimation to separate downgoing seismic wavefield |
SE465643B (sv) * | 1990-02-22 | 1991-10-07 | Bertil Gateman | Elektrooptiskt sensorsystem foer insamling av marina seismiska data |
US4992993A (en) * | 1990-06-18 | 1991-02-12 | Western Atlas International, Inc. | Correction for variable water-column velocity in seismic data |
US5309360A (en) * | 1991-05-23 | 1994-05-03 | Halliburton Geophysical Services, Inc. | Method for attenuating undesirable data, such as multiples, using constrained cross-equalization |
US5293352A (en) * | 1993-01-07 | 1994-03-08 | Western Atlas International, Inc. | Method for removing noise due to near surface scatterers |
US5365492A (en) | 1993-08-04 | 1994-11-15 | Western Atlas International, Inc. | Method for reverberation suppression |
US5396472A (en) * | 1993-09-24 | 1995-03-07 | Western Atlas International | Method for deriving water bottom reflectivity in dual sensor seismic surveys |
US5581514A (en) * | 1993-11-10 | 1996-12-03 | Geco-Prakla, Inc. | Surface seismic profile system and method using vertical sensor |
US5521881A (en) * | 1994-09-02 | 1996-05-28 | Exxon Production Research Company | Method of processing seismic data having multiple reflection noise |
US5621699A (en) * | 1995-07-07 | 1997-04-15 | Pgs Ocean Bottom Seismic, Inc. | Apparatus and method of calibrating vertical particle velocity detector and pressure detector in a sea-floor cable with in-situ passive monitoring |
US5754492A (en) * | 1996-02-12 | 1998-05-19 | Pgs Tensor, Inc. | Method of reverberation removal from seismic data and removal of dual sensor coupling errors |
US5587965A (en) * | 1996-04-26 | 1996-12-24 | Western Atlas International, Inc. | Surface multiple attenuation via eigenvalue decomposition |
US5696734A (en) * | 1996-04-30 | 1997-12-09 | Atlantic Richfield Company | Method and system for eliminating ghost reflections from ocean bottom cable seismic survey signals |
US5850922A (en) | 1996-05-17 | 1998-12-22 | Ryobi North America, Inc. | Shipping and retail display pallet pack |
US5621700A (en) | 1996-05-20 | 1997-04-15 | Schlumberger Technology Corporation, Geco-Prakla Div. | Method for attenuation of reverberations using a pressure-velocity bottom cable |
US5850622A (en) | 1996-11-08 | 1998-12-15 | Amoco Corporation | Time-frequency processing and analysis of seismic data using very short-time fourier transforms |
GB9800741D0 (en) * | 1998-01-15 | 1998-03-11 | Geco As | Multiple attenuation of multi-component sea-bottom data |
GB9810706D0 (en) * | 1998-05-20 | 1998-07-15 | Geco As | Marine seismic acquisition system and method |
GB9906995D0 (en) | 1998-09-16 | 1999-05-19 | Geco Prakla Uk Ltd | Seismic detection apparatus and related method |
US6493636B1 (en) * | 1998-11-05 | 2002-12-10 | Shell Oil Company | Method of marine seismic exploration utilizing vertically and horizontally offset streamers |
US6477470B2 (en) | 2000-12-01 | 2002-11-05 | Pgs Americas, Inc. | Method and system for deghosting |
-
1999
- 1999-03-22 GB GBGB9906456.0A patent/GB9906456D0/en not_active Ceased
-
2000
- 2000-03-21 AU AU33111/00A patent/AU3311100A/en not_active Abandoned
- 2000-03-21 WO PCT/GB2000/001074 patent/WO2000057207A1/en active Application Filing
- 2000-03-21 US US09/936,863 patent/US6775618B1/en not_active Ceased
- 2000-03-21 GB GB0122579A patent/GB2363459B/en not_active Expired - Fee Related
- 2000-03-21 GB GB0122352A patent/GB2363457B/en not_active Expired - Fee Related
- 2000-03-21 US US11/501,195 patent/USRE41656E1/en not_active Expired - Lifetime
- 2000-03-21 US US12/264,784 patent/USRE43188E1/en not_active Expired - Lifetime
- 2000-03-21 US US09/936,745 patent/US6529445B1/en not_active Expired - Lifetime
- 2000-03-21 WO PCT/GB2000/001066 patent/WO2000057206A1/en active Application Filing
- 2000-03-21 AU AU33119/00A patent/AU3311900A/en not_active Abandoned
-
2001
- 2001-09-21 NO NO20014596A patent/NO333154B1/no not_active IP Right Cessation
- 2001-09-21 NO NO20014595A patent/NO332880B1/no not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
NO20014595L (no) | 2001-11-21 |
GB9906456D0 (en) | 1999-05-12 |
GB2363457A8 (en) | 2001-12-19 |
GB2363457B (en) | 2003-04-30 |
WO2000057206A1 (en) | 2000-09-28 |
NO333154B1 (no) | 2013-03-18 |
NO20014596L (no) | 2001-11-20 |
AU3311900A (en) | 2000-10-09 |
US6529445B1 (en) | 2003-03-04 |
GB0122352D0 (en) | 2001-11-07 |
NO20014596D0 (no) | 2001-09-21 |
WO2000057207A1 (en) | 2000-09-28 |
NO20014595D0 (no) | 2001-09-21 |
GB2363459A (en) | 2001-12-19 |
USRE41656E1 (en) | 2010-09-07 |
GB2363459B (en) | 2003-04-30 |
USRE43188E1 (en) | 2012-02-14 |
US6775618B1 (en) | 2004-08-10 |
GB2363457A (en) | 2001-12-19 |
GB0122579D0 (en) | 2001-11-07 |
AU3311100A (en) | 2000-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
NO332880B1 (no) | Fremgangsmate for a redusere effekter fra grov sjo pa seismiske data | |
AU2016200908B2 (en) | Method for acquiring and processing marine seismic data to extract and constructively use the up-going and down-going wave-fields emitted by the source(s) | |
NO20191413A1 (no) | Fremgangsmåte for å kombinere signaler fra trykksensorer og partikkelbevegelsessensorer i marine seismiske streamere | |
AU2002331953B2 (en) | Determination of the height of the surface of a fluid column | |
US7782708B2 (en) | Source signature deconvolution method | |
US9684086B2 (en) | Method for acquiring and processing marine seismic data to extract and constructively use the up-going and down-going wave-fields emitted by the source | |
NO337327B1 (no) | Fremgangsmåte for å kalibrere en seismisk mottaker | |
NO332303B1 (no) | Fremgangsmate for bestemmelse av lokale bolgehoyder og akustisk sensor i marine seismiske signaler | |
NO343375B1 (no) | Seismiske marine langtidsundersøkelser som benytter interpolerte multikomponents streamer-trykkdata | |
US9594180B2 (en) | Removing ghost reflections from marine seismic data | |
NO331334B1 (no) | Demping av stoy og multipler i seismikkdata ved hjelp av sensordata for trykk og vertikal partikkelbevegelse. | |
NO178987B (no) | Fremgangsmåte ved seismiske undersökelser | |
FR2955397A1 (fr) | Procede et dispositif d'acquisition de donnees sismiques marines | |
NO337140B1 (no) | Fremgangsmåte for bestemmelse av dekningsgrad ved marin seismisk kartlegging | |
AU2002331953A1 (en) | Determination of the height of the surface of a fluid column | |
CA2704141A1 (fr) | Methode pour imager une zone cible du sous-sol a partir de donnees de type walkaway | |
NO178648B (no) | Fremgangsmåte for korrigering av feilaktige tidspunkter ved seismiske observasjoner | |
GB2440276A (en) | Processing seismic data acquired using three or more over/under streamers and/or sources | |
AU2014203490A1 (en) | Survey techniques using streamers at different depths | |
US10338251B2 (en) | Method and apparatus for directional designature | |
NO339068B1 (no) | Prosessering av seismiske data som er representative for akselerasjonsbølgefeltet | |
Faggetter et al. | Time-lapse imaging using 3D ultra-high-frequency marine seismic reflection data | |
US20150003197A1 (en) | Inversion techniques using streamers at different depths | |
Hegna | Estimation of the acoustic wavefield generated by a seismic vessel from towed-streamer data | |
Day et al. | Dual-sensor Streamer Data–Calibration, Acquisition QC and Attenuation of Seismic Interferences and Other Noises |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM1K | Lapsed by not paying the annual fees |