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 PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims description 12
- 230000002547 anomalous effect Effects 0.000 claims description 4
- 230000002123 temporal effect Effects 0.000 claims description 4
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 230000000875 corresponding effect Effects 0.000 description 3
- 230000002596 correlated effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application 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.
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)
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)
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 |
-
1996
- 1996-04-15 FR FR9604681A patent/FR2747476B1/fr not_active Expired - Fee Related
-
1997
- 1997-04-09 US US08/973,771 patent/US6032103A/en not_active Expired - Fee Related
- 1997-04-09 WO PCT/FR1997/000628 patent/WO1997039366A1/fr active IP Right Grant
- 1997-04-09 DE DE69710410T patent/DE69710410D1/de not_active Expired - Lifetime
- 1997-04-09 EP EP97918189A patent/EP0839328B1/fr not_active Expired - Lifetime
- 1997-12-12 NO NO19975868A patent/NO318870B1/no not_active IP Right Cessation
- 1997-12-15 OA OA70163A patent/OA10644A/en unknown
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 |