NO318870B1 - Fremgangsmate for lokalisering og identifisering av seismiske anomalier til et medium - Google Patents

Fremgangsmate for lokalisering og identifisering av seismiske anomalier til et medium Download PDF

Info

Publication number
NO318870B1
NO318870B1 NO19975868A NO975868A NO318870B1 NO 318870 B1 NO318870 B1 NO 318870B1 NO 19975868 A NO19975868 A NO 19975868A NO 975868 A NO975868 A NO 975868A NO 318870 B1 NO318870 B1 NO 318870B1
Authority
NO
Norway
Prior art keywords
correlation
seismic
anomalies
maximum correlation
time
Prior art date
Application number
NO19975868A
Other languages
English (en)
Other versions
NO975868L (no
NO975868D0 (no
Inventor
Naamen Keskes
Andre Bidegaray
Original Assignee
Elf Exploration Prod
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Elf Exploration Prod filed Critical Elf Exploration Prod
Publication of NO975868D0 publication Critical patent/NO975868D0/no
Publication of NO975868L publication Critical patent/NO975868L/no
Publication of NO318870B1 publication Critical patent/NO318870B1/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/282Application of seismic models, synthetic seismograms

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Light Guides In General And Applications Therefor (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Description

Foreliggende oppfinnelse angår en fremgangsmåte for å lokalisere og identifisere anomalier til et medium, av den type som de som er funnet i løpet av seismiske feltundersø-kelser.
Bildet av et seismisk areal blir generelt vist som ett eller mange todimensjonale seismiske snitt, referert til som seismisk 2D, definert av aksene x og t, eller ved hjelp av tredi-mensjonale seismiske snitt, hhv. seksjoner, referert til som seismiske 3D,definert av aksen x, y og t eller z, hvor t er tiden og z er dybden.
I en seismisk blokk blir en seismisk hendelse er funnet delvis ved ett eller flere skudd-punkter og fra mottakere forbundet med skuddpunktet definert ved deres koordinater i x og y aksene, og delvis av tiden t det tar fra skuddpunktet til den korresponderende mottaker eller dybdens z hvor den er lokalisert.
En anomali observeres av en tolker som en seismisk hendelse. Studiet av anomalier i et medium gir en bedre forståelse av mediet, siden noen av anomaliene kan være tegn på tilstedeværelse eller fravær av hydrokarboner (vann og olje) i det gitte mediet.
Detektering av anomalier på et seismisk snitt innbefatter et stort antall seismiske traser tidligere satt sammen med innenfor gitte kriterier, som f.eks. ved felles midtpunkter (CMP), ved en felles mottaker, etc. gjøres manuelt av en tolker. Deretter avhenger val-get eller avvisningen av en anomali følgelig kun på bedømmelse foretatt av tolkeren og hans evne/dyktighet til å tolke seismiske snitt riktig. Fig. 1 viser et seismisk snitt (x, t) på hvilke anomaliene, Al og Bl f.eks. er påført av tolkeren. De opptegnede anomaliene kan bli betraktet som enten forskjellige, identiske eller av samme art. Som det fremgår av fig. 1, til venstre og i den andre nedre halvdelen, kan mange anomalier Al som overlapper hverandre bli funnet, noe som gjør det vanskelig å bestemme deres art eller natur. Dersom det er en forkastning i et medium, står man ofte overfor to anomalier som ikke er beslektet. Disse to anomaliene er på forskjellige nivåer og kan bli betraktet som samme fra en ende til en annen av forkastningen mens tolkeren observerer hver som delt og begrenset av forkastningen uten å være i stand til å oppklare tvetydigheten.
I ethvert tilfelle kan manuell detektering av en anomali ikke ta med i betraktningen alle anomalier, spesielt de som er knapt synlige eller de som er skjult av andre anomalier. Foreliggende oppfinnelse er rettet mot en fremgangsmåte for å lokalisere og identifisere hver anomali som forefinnes mellom to forutbestemte nivåer som enten er virkelige eller fiktive horisonter av mediet som skal bli undersøkt.
Oppfinnelsen vedrører således en fremgangsmåte av den innledningsvis nevnte art som er kjennetegnet ved trekkene angitt i karakteristikken til det selvstendige patentkrav 1.
Fordelaktige utførelser av oppfinnelsen er angitt i de uselvstendige patentkravene.
En fordel med oppfinnelsen er at den automatisk sveiper en hel seismisk seksjon ved å markere suksessive intervaller tilliggende den ikke ovenfornevnte seismiske seksjon.
En annen fordel er å detektere hver anomali ved deres konnekse komponenter.
En fordel ligger i det faktum at hver anomali er tatt med i betraktningen og at relevante anomalier av forutbestemte kriterier kan velges. I dette tilfelle blir hver anomali som er av ingen eller liten relevans avvist og betraktet som uviktig.
Ifølge et annet trekk er hver homogen sone betraktet individuelt av en forplanter på en måte for å dekke ovenfornevnte sone i enhver retning mens det foretas kontroll av korrelasjonen med seismiske nabotraser.
Det er således mulig, som følge av foreliggende oppfinnelse, å detektere anomalier som er dekket og vanskelig å detektere manuelt mens man begrenser deres ytre kontur.
Ifølge et annet karakteristisk trekk er homogene soner og deres utvidelser danner anomale soner organisert i et antall lag på en slik måte at innenfor hvert lag er to anomale soner ikke dekket av en annen, idet lagene f.eks. er i en avtagende rekkefølge fra maksimal korrelasjon.
En fordel som ligger i dette trekket skyldes det faktum at anomaliene kan bli kategori-sert i forhold til hverandre, f.eks. ved å følge en reduksjon i maksimal korrelasjon og forskjellsnivåene til anomaliene anbringes i lageret.
Ifølge et annet karakteristisk trekk blir hver anomali validert ved å styre tidskorrelasjonsforholdet på en måte slik at kun anomaliene innenfor en maksimal korrelasjon over en minimumstid, blir holdt.
Dette tillater blant annet å validere hver anomali og danne markeringskort for hver anomali. Hvert kort kan innbefatte amplitude, origo, rommessige koordinater, overflate, etc.
I det påfølgende skal oppfinnelsen beskrives nærmere ved hjelp av en foretrukket utfø-relse og med henvisning til tegningene, hvor:
Fig. 1 viser et 2D seismisk snitt.
Fig. 2 viser en skjematisk fremstilling av en seismisk blokk (x,y,t).
Fig. 3 viser en skjematisk fremstilling av en seismisk trase i et gitt intervall og resultat
av korrelasjonen med en anomali modell.
Fig. 4 viser er forstørret riss av en del av en seismisk seksjon (x, t) inneholdende anomalier. Fig. 5a og 5b viser skjematisk en delvis fremstilling av et kart over maksimal korrelasjon og et isokronalt kart.
Fig. 6 viser en skjematisk fremstilling av beslektede komponenter.
Fig. 7 viser en skjematisk fremstilling av en ordning eller sortering av anomalier.
Fig. 8 viser en skjematisk fremstilling av en validert anomali.
Fig. 9 og 10 viser anomalier sortert på forskjellige nivåer og oppstående fra sorteringen
i skjemaet vist på fig. 7.
Ifølge oppfinnelsen er en 3D seismisk blokk utført som representerer mediet (fig. 2) og som inneholder et stort antall seismiske traser som et resultat, f.eks. av en samling av traser i fellesmidtpunktet (CMP). I blokken 1 defineres et avgrenset temporært intervall med et øvre nivå 2 og et nedre nivå 3, idet de øvre nivåene 2 og de nedre nivåene 3 korresponderer med virkelige eller fiktive horisonter, men at for alle praktiske formål korresponderer for betraktende intervall med en gitt minimumstid tmin og en maksimumstid tmaks-
På fig. 3 er vist en del av den seismiske trasen 4 innenfor nivåene 2 og 3 som korresponderer med hhv. tiden tmjn og W.,. En modell 5, som representerer en anomali er representert som et signal.
Ved første trinn blir modellen 5 korrelert med en del av trasen 4 på en måte for å få korrelert signalet 6 hvor maksimum korrelasjon er TM. Tiden eller indeks to av maksi-mumskorrelasjonen Tm tas. Så blir dette trinnet utført for hver del av den seismiske trasen innenfor nivåene 2 og 3 på en måte for å tilveiebringe verdier for maksimal korrelasjon Tm og merkene tj. Dette tillater dannelse av to kart 7 og 8, hvor ett 7 korresponderer med maksimal korrelasjon f (fig. 5a) og den andre 8 med indeks tj (fig. 5b), sistnevnte er kalt isokronalt kart. De to kartene 7 og 8 er av lik størrelse med den rommessige stør-relsen på blokken 1 og lokalisert i samme aksiale system, x, y for eksempel.
Ved et andre trinn blir anomaliene sortert. En terskel S på korrelasjonen bestemmes og kun de anomaliene med en maksimal korrelasjon rM over terskelen S blir holdt, så blir hver anomali utvidet for å finne konnekskomponentene til anomalien, idet ovenfornevnte utvidelse utføres på anomalier hvor maksimal korrelasjon er større enn S. Søkingen for konnekskomponentene gjøres på et isokronalt kart (fig. 6) på hvilket f.eks. er rappor-tert fire punkter Pl til P4 for koordinatene (ti, xi), fø, X2), fø, X3), og fø, X4) og korresponderer med fire etter hverandre følgende seismiske traser. To punkter, P, Q, tilhører en konneks komponent dersom en bane dannet av punktene til konnekskomponentene forbinder P med Q. De to nabopunktene Pj (xj, tj) og Pj (xj, tj) tilhører samme konnekse komponenter dersom | tj-t;| < At hvor At er en verdi for den forutbestemte terskelen. På fig. 6 tilhører således punktene Pl til P3 samme konnekse komponent på grunn av |t2-t]|
< At og |t3-tz| < At. På den andre siden tilhører punktene P4(x4, t4) ikke disse konnekse komponentene siden |t4-t3| = At' > At.
Et annet kriterievalg eller sortering kan bestå av avvisning av alle anomalier med en maksimal korrelasjon høyere enn S, men hvor størrelsen er mindre enn en gitt størrelse.
Konnekskomponentene Ci, C2, C3 og C4... blir så referert til eller nummerert slik at den isokronale avbildningen har nummererte konnekskomponenter.
Ved et tredje trinn blir hver anomali utvidet med en utbreder for å løse problemet med skjulte anomalier.
På fig. 4 fremgår det at anomalien A2 er unik og har ikke blitt overlappet av andre nær-liggende anomalier. Anomaliene A3 og A4 overlapper hverandre delvis, men på grunn av den dannede utvidelsen av forplanteren er de skilt fra hverandre med klare grenser.
Utvidelsen av hver anomali utføres i hver retning og er styrt nøye ved deres korrelasjon med de seismiske traser til grensene eller konturen til ovenfornevnte anomali med nabo-liggende seismiske traser. Utvidelsen av anomaliene stoppes når ovenfornevnte styrte korrelasjon blir en maksimumskorrelasjon under terskelen S. Forplanteren finner således komplementet av anomalien delvis skjult av en eller flere andre anomalier.
Ved et fjerde trinn blir anomaliene, som er strukket av forplanteren, klassifisert i form av multippelkart hvorved hver inneholder anomalier som ikke dekker hverandre. Kartet til disse anomaliene (fig. 7) er fortrinnsvis klassifisert av avtagende rekkefølge av deres maksimale korrelasjon. Det øvre kartet 10 på fig. 7 korresponderer med høyeste maksimale korrelasjon, mens det nedre kartet 11 korresponderer med svakeste maksimale korrelasjon, idet de to andre kartene 12 og 13 korresponderer med mellomliggende maksimale korrelasjon.
Ved et femte trinn blir anomaliene validert ved å peke ut maksimal korrelasjon, som korresponderer med minimumstider (fig. 8).
Det er mulig å etablere markeringskort for hver anomali, idet hvert markeringskort inneholder informasjon angående ovenfornevnte anomali som f.eks. amplitude, størrelse, overflate, opprinnelse etc.

Claims (8)

1. Fremgangsmåte for å lokalisere og identifisere seismiske anomalier til et medium, karakterisert ved bruk av en gitt seismisk blokk (1) sammensatt av seismiske traser lokalisert fra deres rommessige koordinater, avgrensning i den seismiske blokken ved i det minste ett temporært intervall mellom et øvre nivå (2) og et nedre nivå (3), valg av en temporær modell (5) av en anomali, korrelering av modellen (5) med hver av trasene (4) innenfor det temporære in-tervallet, bestemmelse for hver trase den maksimale korrelasjonen (rM ) og tiden (tj) korresponderende med den maksimale korrelasjon, utførelse av et kart (7) av de maksimale korrelasjonene som er lik den romlige dimensjonen til den seismiske blokken (1) og et kart (8) av tiden som korresponderer med den maksimale korrelasjonen, idet tidskartet har samme dimensjonene og referanse i det aksiale systemet (x, y) som kartet til den maksimale korrelasjonen.
2. Fremgangsmåte ifølge krav 1, karakterisert ved at tidskartet er segmentert i flere konnekse og homogene soner (Ci til G»), idet hver konnekse homogene sone er slik at et punkt i sonen innbefatter i det minste et nabopunkt som ikke inneholder den diskontinuitet i tiden høyere enn en gitt terskel (At).
3. Fremgangsmåte ifølge krav 2, karakterisert ved at alle homogene soner velges hvor det er i det minste ett punkt med maksimal korrelasjon over den forutbestemte terskelen med korrelasjon (S).
4. Fremgangsmåte ifølge krav 2, karakterisert ved at hver homogene sone er en overflate større enn en gitt verdi.
5. Fremgangsmåte ifølge krav 2 til 4, karakterisert ved at hver homogene sone er tatt hånd om individuelt av en forplanter for å spre ovenfornevnte sone i alle retninger, mens korrelasjonen styres med naboseismiske spor.
6. Fremgangsmåte ifølge krav 1 til 5, karakterisert ved at homogene soner og deres utvidelse danner anomale soner som er organisert i et antall lag (10 til 13) på en slik måte at i hvert lag, to anomale soner, hva som helst, ikke dekker hverandre.
7. Fremgangsmåte ifølge krav 6, karakterisert ved at lagene er i en avtagende rekkefølge i forhold til deres maksimale korrelasjon.
8. Fremgangsmåte ifølge ett av kravene 1 til 7, karakterisert v e d at hver anomali valideres ved å styre tidskorrelasjonsforholdet slik at kun anomalier med maksimal korrelasjon for en minimumstid holdes.
NO19975868A 1996-04-15 1997-12-12 Fremgangsmate for lokalisering og identifisering av seismiske anomalier til et medium NO318870B1 (no)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR9604681A FR2747476B1 (fr) 1996-04-15 1996-04-15 Procede de localisation et d'identification des anomalies d'un milieu
PCT/FR1997/000628 WO1997039366A1 (fr) 1996-04-15 1997-04-09 Procede de localisation et d'identification des anomalies d'un milieu

Publications (3)

Publication Number Publication Date
NO975868D0 NO975868D0 (no) 1997-12-12
NO975868L NO975868L (no) 1998-02-06
NO318870B1 true NO318870B1 (no) 2005-05-18

Family

ID=9491202

Family Applications (1)

Application Number Title Priority Date Filing Date
NO19975868A NO318870B1 (no) 1996-04-15 1997-12-12 Fremgangsmate for lokalisering og identifisering av seismiske anomalier til et medium

Country Status (7)

Country Link
US (1) US6032103A (no)
EP (1) EP0839328B1 (no)
DE (1) DE69710410D1 (no)
FR (1) FR2747476B1 (no)
NO (1) NO318870B1 (no)
OA (1) OA10644A (no)
WO (1) WO1997039366A1 (no)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7554883B2 (en) * 2004-10-11 2009-06-30 Landmark Graphics Corporation Fault filter for seismic discontinuity data
US20100153126A1 (en) * 2008-12-17 2010-06-17 Gulfsafe Llc Method and system for verifying an operation in a body of water

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4747054A (en) * 1984-11-28 1988-05-24 Conoco Inc. Method for non-linear signal matching
US4953142A (en) * 1989-01-06 1990-08-28 Marathon Oil Company Model-based depth processing of seismic data
US4964103A (en) * 1989-07-13 1990-10-16 Conoco Inc. Three dimensional before stack depth migration of two dimensional or three dimensional seismic data
GB9206676D0 (en) * 1992-03-27 1992-05-13 Geco As Method of locating hydrocarbon reserves
US5706194A (en) * 1995-06-01 1998-01-06 Phillips Petroleum Company Non-unique seismic lithologic inversion for subterranean modeling

Also Published As

Publication number Publication date
FR2747476B1 (fr) 1998-06-19
NO975868L (no) 1998-02-06
WO1997039366A1 (fr) 1997-10-23
FR2747476A1 (fr) 1997-10-17
EP0839328A1 (fr) 1998-05-06
NO975868D0 (no) 1997-12-12
EP0839328B1 (fr) 2002-02-13
DE69710410D1 (de) 2002-03-21
OA10644A (en) 2002-09-19
US6032103A (en) 2000-02-29
MX9710140A (es) 1998-06-28

Similar Documents

Publication Publication Date Title
CN104951817B (zh) 边界检测装置和边界检测方法
CN101512601B (zh) 从图像中确定深度图的方法以及确定深度图的设备
Rensink et al. Early completion of occluded objects
US7596056B2 (en) Method and program for fracture plane reconstruction
JP2927955B2 (ja) 実時間三次元感知装置
Berger Satellite hydrocarbon exploration: interpretation and integration techniques
EP0489135B1 (en) Method for attribute tracking in seismic data
CN104091324B (zh) 一种基于连通域分割的棋盘格图像快速特征匹配算法
CN104614766B (zh) 地震层位骨架化
NO302544B1 (no) Fremgangsmåte for automatisk plotting og hjelp ved tolkning av seismiske tverrsnitt, spesielt ved anvendelse av bildeanalyse-teknikker
Crosta et al. Geological mapping using Landsat thematic mapper imagery in Almeria Province, South-east Spain
Jazayeri et al. Interest operators for feature‐based matching in close range photogrammetry
NO326418B1 (no) Fremgangsmate for seismisk signalprosessering og utforskning
NO319052B1 (no) Anordning og fremgangsmate ved seismisk signalprosessering og undersokelse
CN105844616B (zh) 激光散射斑点辅助下的双目立体匹配算法与装置
NO332156B1 (no) Fremgangsmate og anordning for a detektere seismiske hendelser, og for a detektere og korrigere geometri- og statikk-feil i seismiske data
NO300608B1 (no) Fremgangsmåte ved plotting av flater i et 3D-volum
CN105954797A (zh) 地震数据的断裂识别方法和装置
Treash et al. Automatic road detection in grayscale aerial images
CN106558143A (zh) 一种100元人民币拼接纸币的识别方法及装置
NO318870B1 (no) Fremgangsmate for lokalisering og identifisering av seismiske anomalier til et medium
CN106772358B (zh) 一种基于cplex的多传感器分配方法
NO334943B1 (no) Fremgangsmåte og apparat for generering av kryssplott i attributtrom fra flere sett med attributt-data og å generere et klassedatasett fra kryssplottet
CA2126700C (en) Error detection in seismic data
CA2225511C (en) Method for locating and identifying site anomalies

Legal Events

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