NO335854B1 - Automatisert system for modellering av flerverdihorisonter med forkastninger - Google Patents
Automatisert system for modellering av flerverdihorisonter med forkastningerInfo
- Publication number
- NO335854B1 NO335854B1 NO20035446A NO20035446A NO335854B1 NO 335854 B1 NO335854 B1 NO 335854B1 NO 20035446 A NO20035446 A NO 20035446A NO 20035446 A NO20035446 A NO 20035446A NO 335854 B1 NO335854 B1 NO 335854B1
- Authority
- NO
- Norway
- Prior art keywords
- fault
- horizon
- data
- model
- reservoir
- Prior art date
Links
- 238000000034 method Methods 0.000 claims abstract description 61
- 238000004590 computer program Methods 0.000 claims description 15
- 238000001914 filtration Methods 0.000 claims description 14
- 238000000638 solvent extraction Methods 0.000 claims description 12
- 230000015654 memory Effects 0.000 claims description 10
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 description 10
- 238000005755 formation reaction Methods 0.000 description 10
- 238000004422 calculation algorithm Methods 0.000 description 9
- 229930195733 hydrocarbon Natural products 0.000 description 9
- 150000002430 hydrocarbons Chemical class 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 239000004215 Carbon black (E152) Substances 0.000 description 7
- 238000012545 processing Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005192 partition Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000005094 computer simulation Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 239000002689 soil Substances 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000013499 data model Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000013213 extrapolation Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- RXZBMPWDPOLZGW-HEWSMUCTSA-N (Z)-roxithromycin Chemical compound O([C@@H]1[C@@H](C)C(=O)O[C@@H]([C@@]([C@H](O)[C@@H](C)C(=N\OCOCCOC)/[C@H](C)C[C@@](C)(O)[C@H](O[C@H]2[C@@H]([C@H](C[C@@H](C)O2)N(C)C)O)[C@H]1C)(C)O)CC)[C@H]1C[C@@](C)(OC)[C@@H](O)[C@H](C)O1 RXZBMPWDPOLZGW-HEWSMUCTSA-N 0.000 description 1
- 108091000069 Cystinyl Aminopeptidase Proteins 0.000 description 1
- 102100020872 Leucyl-cystinyl aminopeptidase Human genes 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000006227 byproduct Substances 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010237 hybrid technique Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 238000004181 pedogenesis Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012800 visualization Methods 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Description
OPPFINNELSENS BAKGRUNN
Oppfinnelsens felt
Den foreliggende oppfinnelsen vedrører datamaskinmodellering av et reservoar med geologiske forkastningshorisonter, og nærmere bestemt en fremgangsmåte, et apparat og et datamaskinprogramprodukt for å modellere et slikt reservoar med i det minste én oppskyvning (eng: reverse fault) og som konsekvens av dette har i det minste én flerverdi-horisont.
Relatert teknikk
Energiindustrien er kontinuerlig involvert i undersøkelse av undergrunnsforekomst av hydrokarboner, slik som olje og naturgass. For å kunne lokalisere og estimere mengden hydrokarbon som er tilstede i slik forekomst på en kostnadseffektiv måte, og også for å tilveiebringe måter for å gjenvinne disse på effektiv måte med den eksisterende teknologien, gjør energiindustrien utstrakt bruk av datamaskinmodelleringsteknikker.
Slike modelleringsteknikker er beskrevet i ytterligere detalj i de følgende US patenter, som herved innlemmes ved referanse i denne beskrivelsen: 1) US patent nr. 5 982 707, med tittel "Method And Apparatus For Determining Geologic Relationships For Intersecting Faults", meddelt 9. november 1999, William E. Abbott. 2) US patent nr. 6 014 343, med tittel "Automatic Non-Artifically Extended Fault Surface Based Horizon Modeling System", meddelt 11. januar 2000, Graf o.a. 3) US patent nr. 6 138 076, med tittel "Automatic Non-Artifically Extended Fault Surface Based Horizon Modeling System", meddelt 24. oktober 2000, Graf o.a..
Jordformasjoner omfatter typisk horisonter og forkastninger. Det er kjent at hydrokarbonforekomster er tilbøyelige til å plasseres nærliggende skjæringspunktene mellom disse horisonter og forkastningene. Tidligere kjente modelleringsteknikker var typisk i stand til å håndtere normale forkastninger men ikke oppskyvninger. F.eks., i en tidligere teknikk var en oppskyvning behandlet som om det var en vertikal forkastning, og derved eliminerte man overlappingshorisontoverflater. Imidlertid reduserer en slik forenkling trofastheten ved modellen og øker feilen.
I en annen tidligere kjent løsning, er hydrokarbonjordreservoaret delt i felter, henvist til som blokker, med anvendelse av alle forkastningene og ved å kunstig utvide dem for å danne en lukket avdeling i reservoaret. Imidlertid, fører slike teknikker til ikke-realistiske og kompliserte datamaskinmodeller og derved setter man i fare muligheten til å forstå jordformasjonen og å avlede nøyaktig estimater for hydrokarbonreserver. Dynamic Graphics "Earth Vision" anvender denne modelleringstilnærming, som et eksempel.
I enda en annen tidligere kjent løsning, krever modelleringsteknikken at brukeren lager reservoarmodellen for hånd. F.eks. Technoguides "Petrell"-system eller Roxars "IRAP RMS"-system anvender denne modelleringsteknikken som et eksempel. Imidlertid er slike teknikker langtekkelige og tidskrevende. Denne teknikken er tilbøyelig til å gjøre mennesket til et ekstrautstyr for datamaskinen ved å kreve at brukerne på interaktiv måte bygger geometriske modeller, heller enn å beregne modellene ut fra dataene. I tillegg er disse teknikkene ikke i stand til effektivt å modellere reservoarer som har sammensatte forkastninger.
Chambers, K.T. et al.: Geologic modeling, upscaling and simulation of faulted reservoirs using faulted stratigraphic grids. Proceedings of the SPE reservoir simulation symposium; Houston, TX, USA. 1999.02.14 - 1999.02.17. Nr. 51889; Sider 1 - 12. XP00221591, Wheeler M.F. et al.: Parallel multiblock/multidomain approach for reservoir simulation, Proceedings of the SPE reservoir simulation symposium; Houston, TX, USA. 1999.02.14 - 1999.02.17. Nr. 51884; Sider 1-11. XP002215913, US 5740342 og US 6014343 omhandler alle automatisk modellering av forkastninger ved reservoarer i undergrunnen omfattende å inndele reservoaret i 3D-gitter-celler.
Ut fra det ovennevnte bør man se at det finnes et økende behov for å ha en modelleringsteknikk som kan håndtere flere typer forkastninger, omfatter oppskyvninger, og som er i stand til å produsere realistiske og nøyaktige modeller. I tillegg er det ønskelig at modelleringsteknikken støtter automatisering og derved krever veldig lite brukerpåvirkning.
Kort beskrivelse av oppfinnelsen
Det foregående behov er møtt ved den foreliggende oppfinnelsen. Ifølge en form av oppfinnelsen, i en fremgangsmåte for å modellere et reservoar som har en horisontflate med en oppskyvning, er et volum i reservoaret oppdelt i forkastningsblokker. En grenseflate i hver av forkastningsblokkene inkluderer en del av oppskyvningen. De geologiske horisontene på innsiden av de således definerte forkastningsblokkene er representert med overfiatemodeller med enkeltverdier. De rå horisontdataene er oppdelt i flatebiter (eng: patches), hvor hver flatebit tilsvarer en forkastningsblokk. Horisontoverflateflatebitmodellen er da avledet fra dataflatebiten inne i hver forkastningsblokk. Horisontflatebitmodellen omfatter en oppover eller nedover kastet forkastningstrase som definerer skjæringen av overflateflatebiten med den grensende oppskyvningen. Kontinuitetsbegrensninger er avledet automatisk for å sikre at horisontoverflateflatebitmodellene på motsatte sider av en oppskyvning samsvarer langs disse partiene i forkastningen som ikke representerer den virkelige geologiske forkastning, men som er en utvidelse av den virkelige forkastningsmodellen som er anvendt for å dele reservoarets volum i veldefinerte, lukkede blokker.
Ifølge et annet aspekt ved oppfinnelsen, er det modellert et reservoar med horisontflater og forkastninger, både oppskyvninger og normale forkastninger. Ved ett trekk ved denne fremgangsmåten, er forkastningsblokkene oppdelt med anvendelse av en undergruppe av alle forkastningene i reservoaret, imidlertid, omfatter undergruppen alle oppskyvninger. Som konsekvens av dette er hver oppskyvning inkludert i en grenseflate for i det minste én av forkastningsblokkene, og alle de oppdelte forkastningsblokkene har bare enkeltverdihorisontflater.
Andre utførelser, og også hensikter og fordeler ved oppfinnelsen vil komme frem ved å lese den følgende detaljerte beskrivelse og med henvisning til de vedlagte tegningene.
Kort beskrivelse av tegningene
De nye trekk som er antatt å være karakteristiske for oppfinnelsen er satt frem i de vedlagte kravene. Oppfinnelsen i seg selv, imidlertid, og også en foretrukket anvendelsesmåte, ytterligere hensikter og fordeler derav, vil forstås bedre med henvisning til den følgende detaljerte beskrivelse av en illustrerende utførelse når den leses sammen med de vedlagte tegningene, hvor: Fig. 1 illustrerer en tredimensjonal representasjon av et parti av en undergrunnsformasjon. Fig. 2A, 2B, 2C, 2D og 2E illustrerer typer forkastninger som typisk oppstår i en undergrunnsformasjon. Fig. 3A, 3B og 3C illustrerer enkeltverdi- og flerverdi-horisontoverflater, omfattende aspekter av en utførelse av oppfinnelsen. Fig. 4A og 4B illustrerer karakteristikker til en skjæringskurve mellom en forkastning og en horisont, ifølge en utførelse av oppfinnelsen. Fig. 5 illustrerer en modelleringsteknikk for å utvikle en horisontmodell i et tredimensjonalt reservoar med i det minste en oppskyvning, ifølge en utførelse av oppfinnelsen. Fig. 6A, 6B, 6C, 6D, 6E, 6F og 6G illustrerer ytterligere aspekter ved deling av reservoaret i forkastningsblokker, ifølge en utførelse av oppfinnelsen. Fig. 7A og 7B illustrerer dannelse av flatebitene, ifølge en utførelse av oppfinnelsen.
Fig. 8 illustrerer filtrering av data, ifølge en utførelse av oppfinnelsen.
Fig. 9A og 9B illustrerer utvikling av kontinuitetsgrensene, ifølge en utførelse av oppfinnelsen. Fig. 10 illustrerer en teknikk for å utvikle de flerverdi-horisontalmodeller ifølge en utførelse av oppfinnelsen. Fig. 11 illustrerer et blokkdiagram av et datamaskinsystem for å implementere fremgangsmåte eller apparataspektet ved den foreliggende oppfinnelsen ifølge en utførelse av oppfinnelsen. Fig. 12 illustrerer et blokkskjema av et klientserverdatamaskinsystem, ifølge en utførelse av den foreliggende oppfinnelsen. Fig. 13 illustrerer et blokkskjema av et parallellprosesseringsdatamaskinsystem, ifølge en utførelse av den foreliggende oppfinnelsen.
Detaljert beskrivelse av den foretrukkede utførelsen
I den følgende detaljerte beskrivelsen av de foretrukkede utførelsene, viser man til de vedlagte tegningene som illustrerer utførelser hvor oppfinnelsen kan utøves. Man bør imidlertid forstå at tegningene og den detaljerte beskrivelsen ikke er ment å begrense oppfinnelsen til den spesielle formen som beskrives, men derimot, er meningen å dekke alle endringer, ekvivalenter og alternativer som faller innenfor ånden og rammen for den foreliggende oppfinnelsen som definert av de vedlagte kravene.
Med henvisning til fig. 1, er en tredimensjonal (3D) representasjon av et parti av en undergrunnsformasjon, dvs. et reservoar, illustrert. Partiet av undergrunnsformasjonen omfatter et antall horisonter 110, 120, 130, 140, og forkastninger 160 og 170, som skjærer horisontene. Hydrokarbonforekomster som er inkludert i reservoaret 115 er tilbøyelige til å plasseres nærliggende skjæringskurvene mellom horisonter og forkastninger, f.eks. ved skjæring 125 av horisont 110 med forkastning 160. Den endelige hensikten ved horisontmodelleringssystemet ifølge den foreliggende oppfinnelsen er å hjelpe en bruker, f.eks. en geofysiker eller en reservoaringeniør, i oppgaven å tolke innmatingsdata, f.eks. brønnlogg og seismiske data, for å forutsi optimale lokaliseringer av undergrunnsforekomster av hydrokarbonenergi som er lagret innenfor jordformasjonen. Modellen er typisk anvendt for å vurdere mengden hydrokarboner i reservoaret og også effektive måter for å gjenvinne disse. Denne analysen er da anvendt for å ta forretningsavgjørelser vedrørende å støtte utvikling og håndtering av reservoaret.
Det tredimensjonale partiet av reservoaret 115 er en lukket blokk hvor grensen er definert av forkastningen 160 og partier av grensen til hele reservoaret. I det illustrerte eksemplet definerer grensen et firkantet fast legeme med et volum 135 med størrelse 80,46km (50 miles) lang, 32,19 km (20 miles) bred og 8,04km (5 miles) dypt.
En geologisk horisont er grensesnittet mellom to sedimentasjonsjordlag. En forkastning er et brudd i jordlaget og de nærliggende horisontflater, og over dette finnes det betraktelig forskyvning.
Med henvisning til fig. 2A, 2B, 2C, 2D og 2E, er flere typer forkastninger illustrert. Siden det er sannsynlig at hydrokarbonforekomst er lokalisert nærliggende skjæringen mellom horisont og forkastninger, er fastsettelse av lokaliseringen og type forkastning, og også en nøyaktig løsning av geometrien til alle de relaterte geologiske trekk, viktig. Avhengig av den relative retningen til forskyvningen mellom horisontflatene, er en forkastning beskrevet som en normal forkastning 210 vist i fig. 2A, en oppskyvning 220 vist i fig. 2B, en horisontal forkastning 230 vist i fig. 2C eller en blandet forkastning (ikke vist). For tverrsnittrisset av den normale forkastningen 210, overlapper ikke de horisontale flatene 21 la og b på respektive sider av forkastningsoverflaten 212 hverandre vertikalt som vist i fig. 2D. For tverrsnittrisset av oppskyvningen 220, overlapper horisontflatene 221a og b hverandre vertikalt som vist i fig. 2E. Overlappingsområdet er henvist til som flerverdiområdet 215. For den horisontale forkastning 230, er horisontflatene 231a og b forskjøvet lateralt men ikke vertikalt. En forkastning er blandet hvis, i noen områder langs forkastningsoverflaten, er den normal og i andre er den en oppskyvning eller en horisontal forkastning.
Med henvisning til fig 3A, 3B og 3C, illustreres skjæringskurven med et x-y-plan for enkeltverdi og flerverdi horisontflater 211 og 221 som resulterer fra forkastninger 210 henholdsvis 220. Fig. 3A illustrerer at en normal forkastning 210 resulterer i enkeltverdioverfiaten 211 i x-y-planet siden det ikke er noen vertikal overlapping mellom de tilsvarende horisontflater. Det vil si at i x-y-planet som skjærer horisont 211 finnes det én og bare én y-koordinatverdi. F.eks. i x-y-planet for en enkelt x-koordinat xl, finnes det bare et horisontpunkt eller en gitternode 310 for horisont 211a.
Fig. 3B illustrerer en flerverdig horisontflate 221 som resulterer fra oppskyvningen 220. Et område i et reservoar med flerverdi-horisontflater er kjent som flerverdiområdet 215.1 horisont 221 som skjærer x-y-planet, finnes det flere horisontpunkter eller gitternoder, hvor hver gitternode er tilknyttet en tilsvarende horisontflate. F.eks., for enkeltkoordinat x2, finnes det to tilsvarende gitternoder 325 og 345 i horisont 221. Gitternode 325 tilsvarer en oppoverrettet forkastningstrase for horisontflate 221a og gitternode 345 tilsvarer en nedoverrettet forkastningstrase for horisontflate 221b. Begge disse trasene vil beskrives ytterligere i fig. 4A, 4B og 4C nedenfor.
Trasen til en horisontflate som skjærer en forkastningsoverflate er en éndimensjonal kurve i det tredimensjonale Euclidsk rom. F.eks., der hvor en flat kontakt som representerer en forkastningsoverflate skjærer en flat horisontflate, er trasen (kartmønster) til kontakten en rett linje. Hvis kontakten ikke er plan er trasen en generell kurve. En forkastningstrase er en kurve som representerer skjæringskurven for horisonten og forkastningen. De oppoverrettede og nedoverrettede forkastningstrasene omfatter par av oppoverrettede og nedoverrettede skjæringslinjer mellom horisonten og hver forkastningsflate.
Fig. 3C illustrerer en flerverdi-horisontflate 351 som resulterer fra to oppskyvninger som er i nærheten av hverandre, og som danner en sammensatt oppskyvning. I tilfellet ved denne sammensatte oppskyvningen 380, for en enkelt koordinat x3 i x-y-planet som vises, finnes det tre tilsvarende gitternoder 350, 360 og 370 som skjærer horisontflaten 351.
Med henvisning til fig. 4A, 4B og 4C, illustreres karakteristikker for forkastningstraser for en skjæring mellom forkastning 160 og horisont 110 i fig. 1. Et forkastningsområde som omfatter en forkastningstrase som går oppover 475 og en forkastningstrase som går nedover 485 er tilveiebragt av skjæringskurven 125 mellom forkastningen 160 og horisonten 110. Forkastningsområdet er en åpning som er dannet i horisonten 110 når forkastningen 160 passerer gjennom horisonten 110. Den venstre siden av forkastningsområdet er vist løftet oppover 420, og den høyre side er vist senket nedover 430. Fig. 4B illustrerer et toppriss av horisonten 110 som viser forkastningsområdet 410. Fig. 4C illustrerer siderisset av horisonten 110 og en kryssforkastning 160 i fig. 4A tatt langs snittlinjer 11-11 i fig. 4A. Merk at den venstre siden av horisonten 110 er hevet sammenlignet med den høyre side av horisonten 110.
Med henvisning til fig. 5, er en modelleringsteknikk for å utvikle en horisontmodell for 3D reservoaret 115 som har i det minste én oppskyvning 220 illustrert, ifølge en utførelse av den foreliggende oppfinnelsen. Den tredimensjonale horisontmodell i den foreliggende oppfinnelsen er tilpasset for å motta horisontoverflatedata, forkastningsoverflate og forholdsdata som innmatninger. Kjente teknikker for å generere horisontinnmatingsdata, også henvist til som horisontdata, er beskrevet i sin helhet i US patentreferansene som var nevnt tidligere. Horisontdata omfatter data som er oppnådd for aktiviteter slik som boring og også oppsamling av dataprøver og utførelse av seismiske undersøkelser for å samle seismiske data. I tillegg til horisontdata, omfatter innmatingskravene forkastningsoverflater og en beskrivelse av forkastningsforhold, f.eks. navn av forkastninger som krysser hverandre i deres største/minste forhold. Kjente teknikker for å generere forkastningsoverflater og forholdsdata, f.eks. informasjon som beskriver forkastningsoverflater og forkastningsforhold er også beskrevet i sin helhet i US patentreferansene som er nevnt tidligere.
Den nye modelleringsteknikken omfatter i trinn 510, å dele et reservoar som har flerverdihorisontflater i forkastningsblokker. Reservoaret kan ha flere forkastninger, hvor bare noen av dem er oppskyvninger. Alle oppskyvningene er inkludert blant oppdelingene som er valgt som grenser for forkastningsblokker, slik at alle de oppdelte forkastningsblokkene har bare enkeltverdihorisontflater, selv om de kan være forkastet med normale forkastninger. De delte forkastningsblokkene med enkeltverdihorisontflatene er da utviklet videre med anvendelse av kjente teknikker som er beskrevet i sin helhet i US patent referansene som nevnt tidligere.
I trinn 510, er et volum i reservoaret delt i et flertall av forkastningsblokker. I én utførelse, omfatter alle avdelingene oppskyvninger, slik at i det minste én grenseflate for hver av flertallet i forkastningsblokkene omfatter et parti av et respektivt en av oppskyvningene. Oppdelingsprosessen er gjentatt for hver av de ytterligere oppskyvningene som kan være tilstede i reservoaret 115. For å gjenta, ved enden av delingsprosessen, på innsiden av hver av flertallet av forkastningsblokkene har de geologiske horisontene enkelt verdi. Avdelingsprosessen i trinn 510 for reservoaret med en oppskyvning er beskrevet i ytterligere detalj i fig. 6A og 6B, ifølge en utførelse av den foreliggende oppfinnelsen. Avdelingsprosessen i trinn 510 for reservoaret med flere oppskyvninger er beskrevet i ytterligere detalj i fig. 6C, 6D, 6E, 6F og 6G, ifølge en utførelse av den foreliggende oppfinnelsen.
I trinn 520, er horisontdataene som beskriver reservoaret 115 avdelt i datagrupper, hvor hver datagruppe tilsvarer forkastningsblokkene generert i trinn 510. For en gitt horisontflate, omfatter hver forkastningsblokk bare datapunkter som definerer enkeltverdi horisontmodell på innsiden av forkastningsblokken, som beskrevet i trinn 510.
En horisontflate innenfor en forkastningsblokk er kjent som en horisontflatebit, eller simpelthen en flatebit. Det vil si at for en gitt horisont, finnes i det meste en horisontflatebit, muligens med forkastninger, for hver forkastningsblokk, men noen av forkastningsblokkene kan mangle flatebiter. Antall datagrupper er mindre enn eller lik antall forkastningsblokker generert ved avdeling av reservoaret 115. Hver datagruppe tilsvarer en horisontflatebit. Gruppen med alle horisontflatebitene danner den fullstendige horisontmodellen i reservoaret. Dannelse av flatebitene for hvert trinn 520 er beskrevet i nærmere detalj i fig. 7.
I trinn 530, kan en valgfri filtrering av data utføres. Filtreringen kan være nødvendig for å fjerne forskjeller mellom innmatingshorisontdata og geometrien for forkastningsblokkene definert ved avdelingsforkastningene. Filtrering av data for trinn 530 er beskrevet i ytterligere detalj i fig. 8.
I trinn 540, er kontinuitetsbegrensninger utviklet langs kunstig utvidede forkastningsblokkgrenser for å beholde integriteten til horisontflaten og dens gradient. Nærmere bestemt, sikrer kontinuitetsbegrensningene at man opprettholder integriteten til en kontinuerlig horisontoverflatemodell og gradienten langs kunstig utvidede forkastingsblokkgrenser som anvendes for å bygge lukkede forkastningsblokkavdelinger. I en utførelse, resulterer kontinuitetsbegrensningene i samsvarende endepunkter for den oppoverrettede og den nedoverrettede forkastningstrasen til oppskyvningen. Utvikling av kontinuitetsbegrensningene for trinn 540 er beskrevet i nærmere detalj i fig. 9, ifølge en utførelse.
For å gjenta, det foregående har beskrevet å velge oppskyvninger for å avdele reservoaret slik at alle forkastningsblokkene i reservoaret har enkeltverdihorisonter. Med henvisning til fig. 6A, 6B, 6C, 6D, 6E, 6F og 6G er avdeling av et reservoar i forkastningsblokker, ved trinn 510 illustrert. Forkastningsblokkdelingsteknikken innlemmer den klassiske matematiske enhetsdelingskonsept.
Gitt en ramme for overflater og forkastninger som representerer brudd i reservoaret, og gruppen med innmatingshorisontdata (typisk en sky med punkter i det tredimensjonale Euclidske rom som inneholder reservoarets volum) som beskrevet i tidligere nevnte US patenter, er avdeling av volumet til reservoaret i lukkede undervolum er utviklet. Undervolumene er vanligvis kjent som forkastningsblokker. Bristene i reservoaret er ikke nødvendigvis koblet til hverandre, slik at de tilsvarende forkastnings overflatene definerer ikke en lukket forkastningsblokkpartisjon. For å overvinne dette, er den følgende prosedyren anvendt for å ekstrapolere, eller utvide, forkastningsoverflaten.
Et kartesisk koordinatsystem (x,y,z) for reservoaret og parameterisering ZF(x,y) for en forkastningsoverflate F er definert, hvor Z er en kontinuerlig deriverbar funksjon i definisjonsdomenetZ)F(x,y) for forkastningen. Reservoaret kan ses som et lukket volum hvor overflaten er sammenslutningen av overflater Si for i = 1, ..., n, for et naturlig tall n. Hver overflate Si tillater en parameterisering Si i et definisjonsdomene Dsi(x,y). Den opprinnelige forkastningsoverflate er utvidet gjennom reservoaret med anvendelse av en velkjent matematisk prosedyre for overflateekstrapolering. Ekstrapoleringen, eller utvidelsen, av forkastningsoverflaten kan uttrykkes som en ny parameterisering Z<A>F, definert over hele reservoaret, f.eks. sammensetning av alle Dsi(x,y) og slik at Z<A>F(x,y) = ZF(x,y) for hvert punkt (x,y) i DF(x,y). Naturligvis er D F(x,y) en undergruppe av sammensetningen av alle Dsi(x,y).
Med henvisning til fig. 6A er et flytskjema for avdeling av reservoaret med en oppskyvning illustrert, ifølge en utførelse. I trinn 610, er forkastningsoverflaten som omfatter oppskyvningen 220 kunstig utvidet for å krysse reservoargrensene. I trinn 620, er den kunstig utvidede forkastningsoverflate innført i volumet 135 til reservoaret og derved splittes eller avdeles volumet til reservoaret i i det minste to forkastningsblokker.
Med henvisning til fig. 6B, hvis man tar den første ekstrapolerte forkastning ZAF, f.eks., utvidet forkastningsoverflate 680, og innfører den i volum 135 til reservoaret vil det splitte det opprinnelige volum 135 i to eller flere lukkede undervolum eller forkastningsblokker 650 og 660 som illustrert, ifølge en utførelse. I det minste en grenseoverflate for hver av forkastningsblokkene 650 og 660 omfatter et parti av oppskyvningen 220 med endepunkter 610 og 620. Den viste oppoverrettede forkastningstrase 475 nedoverrettede forkastningstrase 485 omfatter endepunktene 610 og 620.1 en utførelse, er forkastningsblokkene formet ved kunstig å utvide forkastningsoverflater i oppskyvningen 220 til en skjæringskurve med en annen horisont med en annen forkastning.
Ytterligere forkastningsblokker er generert for ytterligere oppskyvninger, når de er til stede. Med henvisning til fig. 6C, er det illustrert et flytskjema for å avdele eller partisjonere reservoaret med en andre oppskyvning, ifølge en utførelse. I trinn 630, er forkastningsoverflaten som inkluderer den andre oppskyvningen kunstig utvidet for å skjære reservoargrensene. I trinn 640, er den kunstig utvidede forkastningsoverflate i den andre oppskyvningen innført i forkastningsblokken som omfatter den andre oppskyvningen. Innføringen resulterer i splitting eller partisjonering eller avdeling av volumet i forkastningsblokken i i det minste to nye forkastningsblokker. I det valgfrie trinnet 650, hvis flere enn to oppskyvninger er inkludert i reservoaret vil da trinn 630 og 640 gjentas for hver av de ytterligere oppskyvningene.
Fig. 6D, 6E, 6F og 6G illustrerer partisjonering av reservoaret med to oppskyvninger inkludert i forkastningsoverflatene 645 og 655. Ved å innføre den andre oppskyvningen i en hvilken som helst av disse nye blokkene vil den splitte den i to eller flere nye forkastningsblokker. Hvis man fortsetter denne prosessen til den siste ekstrapolerte forkastning er innført definerer man en gruppe med lukkede forkastningsblokker 655, 625 og 635 hvor unionen er det opprinnelige
reservoarvolumet 135, som illustrert i fig. 6D ifølge en utførelse. Fig. 6E illustrerer en lukket forkastningsblokk 615 og forkastningsoverflate 645 i oppskyvningen. Fig. 6F illustrerer lukket forkastningsblokk 625 og forkastningsoverflate 655 i den andre oppskyvningen. Fig. 6G illustrerer lukket forkastningsblokk 635. En ekstrapolert
forkastning splitter bare forkastningsblokken hvor den er innført. En ekstrapolert forkastningsoverflate kan innføres i mer enn én forkastningsblokk om ønskelig.
Et viktig element i forkastningsblokkpartisjoneringsprosessen, ved trinn 510, og som beskrevet ovenfor er friheten til å danne lukkede forkastningsblokker med anvendelse av bare en undergruppe av alle forkastningene i reservoaret. F.eks. kan en undergruppe anvendes som inkluderer bare oppskyvninger, men som inkluderer alle oppskyvningene i reservoaret. Dette vil med fordel tillate at oppskyvninger, som er viktige, modelleres mer nøyaktig enn en modell med forenklede oppskyvninger ved å behandle dem som vertikale forkastninger. Som et alternativ vil undergruppen av forkastningene som anvendes for å partisjonere reservoaret i forkastningsblokker ikke nødvendigvis være begrenset til oppskyvninger. Et hvilket som helst antall eksisterende normale forkastninger kan anvendes i partisjoneringsprosessen beskrevet i trinn 510.
Med henvisning til fig. 7 A, et flytskjema for å partisjonere dannelse av flatebiter er illustrert, ifølge en utførelse. I trinn 710, for en gitt horisont, er forkastningsgrensene kunstig utvidet. I trinn 720, er horisonten partisjonert eller avdelt i to flatebiter ved kunstig å utvide forkastningsgrensene. De to fiatebitene er på motsatte sider av de utvidede forkastningsgrensene.
Med henvisning nå til fig. 7B, er dannelse av fiatebitene illustrert, ifølge en utførelse. Det første trinn i å utvikle horisontflatebitene i horisontdataene er å partisjonere horisontinnmatingsdataene i forhold til et sett med forkastningsblokker, ved trinn 520, hvor unionen er hele reservoaret. Resultatet av å partisjonere horisontdataene er en samling av datagrupper. Antall datagrupper i denne samlingen er mindre enn eller likt antall forkastningsblokker som er generert ved å splitte eller avdele reservoaret. For en gitt horisont 110, vil flatebiter 710 og 720 dannes ved kunstig å utvide forkastningsblokkgrensene 680 som beskrevet i fig. 6B. De partisjonerte datagrupper med hver datagruppe tilsvarende en horisontflatebit er da anvendt i de senere trinnene i modelleringsalgoritmen.
Horisontmodellen ved innsiden av hver forkastningsblokk er representert som en enkelt verdi overflate. På denne måten, vil forkastningshorisontmodellering være en hybrid teknikk, som tillater både volumsplitting og ikke-splittingsforkastninger som valgt av brukeren.
Med henvisning til fig. 8, illustreres filtrering av data, ifølge en utførelse. Innmatingsdatapunktene 810 kan om ønskelig filtreres i forhold til grenseoverflatene til forkastningsblokken. Som beskrevet tidligere, kan filtreringen være nødvendig for å fjerne avvik mellom innmatingshorisontdataene og geometrien til forkastningsblokken som er definert ved splittingsforkastningene beskrevet i trinn 510 på fig. 5. Hovedkilden for innmatingshorisontdataene er typisk oppnådd fra seismisk tolkning. Det er kjent at tolkingen vil være tilbøyelig til å være unøyaktig i nærheten av forkastninger eller sprekker, f.eks. innenfor et smalt filtreringsområde 830. Ofte vil horisontdatapunkter tolkes over forkastningen uten å ta i betraktning den vertikale forskyvningen som en horisont 815 underkastes ved lokaliseringen av bristen. Hvis den beholdes ubehandlet, vil denne feilaktige tolkningen kunne føre til feil i horisontmodellen som er beregnet og deretter i modellen for hele reservoaret. En applikasjon som er i stand til å fjerne innmatingsdatapunkter basert på deres avstand fra forkastningsoverflaten kan anvendes for å utføre det valgfrie filtreringstrinnet. Den euklidske avstanden fra punktet til den opprinnelige (ikke ekstrapolerte) forkastningsoverflate er anvendt. Denne applikasjonen kan konfigureres for å fjerne punkter, slik som punkter 820, selektivt på én eller på begge sider av en bestemt forkastning 805.
Med henvisning til fig. 9A, er et blokkskjema for dannelse av flatebiter illustrert, ifølge en utførelse. I trinn 910, som er samme som trinn 710, for en gitt horisont, er forkastningsgrensene kunstig utvidet. I trinn 920, er kontinuitetsbegrensninger i form av funksjonsverdier og gradienter utviklet på hver side av de kunstig utvidede forkastningsgrensene.
Med henvisning til fig. 9B, er utvikling av kontinuitetsbegrensningene 905 illustrert, ved trinn 540 ifølge en utførelse. Kontinuitetsbegrensningene 905 må generelt sikre at en horisontflate og dens gradient er kontinuerlige over kunstig utvidede forkastningsblokkbegrensninger. Uten innføring av
kontinuitetsbegrensninger 905, vil horisontflatebitene 710 og 720 fra begge sider av den kunstig utvidede blokkgrensen 680 være diskontinuerlige. En slik diskontinuitet vil også føre til uoverenstemmelse ved endepunktene 610 og 620 i sprekkåpningen til horisonten forårsaket av det virkelige partiet 920 til oppskyvningen 220. Dette vil resultere i en lavkvalitets horisontmodell og føre til unøyaktige modelleringsresultater i de senere trinnene i arbeidsflyten, f.eks., ved eiendomsmodellering og strømningssimulering.
For å imøtekomme slike problemer, er kontinuitetsbegrensninger 905 i form av funksjonsverdier og gradient beregnet og pålagt langs den kunstige blokkgrensen
680. Hvis (xt,yt, zt) er parameteriseringen av krysset mellom den ekstrapolerte delen av en forkastning F og en horisontflatebit A og (x<A>t,y<A>t,z<A>t) er parameteriseringen av skjæringen mellom den samme delen F med horisontflatebit B på den andre siden av F, vil systemet sikre at (xt,yt, zt) = (x<A>t,y<A>t,z<A>t) for alle verdier av det reelle
parameter t.
I tillegg la HA(x,y) og HB(x,y) være parameteriseringene til horisontflatebitene A og B, f.eks. flatebiter 520 og 530, i nærheten av krysset med forkastning F. Horisontflatebiten på innsiden av en gitt forkastningsblokk kan ikke antas å tillate en mangfoldig (ikke forkastet) representasjon på grunn av antagelsen av at innsiden kan inkludere normale forkastninger. Imidlertid vil horisontoverflateområdene mellom sprekker kunne parameteriseres ved kontinuerlig dreriverbare kart H(x,y) hvor (x,y) hører til en undergruppe av unionen av alle Dsi(x,y). Som en konstruksjon, HA(x,y) = (xt,yt,zt) for alle verdiene til det reelle paramter t.
Da vil systemet sikre at grad HA(xt,yt) = grad HB(xt,yt) for alle verdiene av det reelle parameter t.
Skjæringskurvene til flatebitoverflaten ved grensene i den tilsvarende forkastningsblokk er kjent som flatebittrimmekurver (eng: patch trim loops). De er et viktig element i flerflatebithorisontmodellen. Kontinuitetsbegrensningene 905 sikrer at trimmekurvene for to nabohorisontflatebiter 710 og 720 sammenfaller over de ekstrapolerte delene av de felles forkastnings overflatene. Imidlertid vil trimmekurvene divergere i området for den virkelige delen 920 av oppskyvningen 220 og således definere forkastningsområdet 410, også kjent som sprekkåpningen i horisonten forårsaket av forkastningen.
Med henvisning nå til fig. 10, er en teknikk for å utvikle flerverdihorisontmodeller illustrert i et flytskjema, ifølge en utførelse av oppfinnelsen. Teknikken beskrevet nedenfor innlemmer trinn 520, 530 og 540 i fig. 5 som trinn 1010, 1020 og henholdsvis 1030. De etterkommende trinn 1040 til 1080 i teknikken følger logikken til algoritmen for enkeltverdihorisont beskrevet i US patentene som er vist til tidligere men anvendt bare på innsiden av forkastningsblokken. En database 1090 som omfatter en representasjon av en jordmodell er om ønskelig konfigurert for å levere innmating til og/eller motta utmating fra hvert av trinnene 1010 til 1080.1 en utførelse, er innmatingene til og/eller utmatingene fra hvert av trinnene 1010 til 1080 lagret i minnet. Overflateoperasj onene som er utført for å utvikle forkastningsflatebitmodellen er begrenset til innsiden av den tilsvarende forkastningsblokk, på grunn av det faktum at horisontflatebiten ikke er definert utenfor den tilsvarende forkastningsblokken. Dette eliminerer behovet for å klassifisere de resterende forkastningene fra rammen i forhold til forkastningsblokkpartisjonering av reservoaret.
I trinn 1010, partisjoneres innmatingsdataene lagret i databasen 1090 i forhold til forkastningsblokksplitting av reservoaret og defineres en gruppe med flatebiter for horisontdata, med hver flatebit tilsvarende en forkastningsblokk. De partisjonerte horisontdataene er lagret i databasen 1090.
I valgfritt trinn 1020, filtreres innmatingsdataene i forhold til grenseoverflatene i forkastningsblokkene for å mate ut filtrerte, partisjonerte horisontdata som er lagret i databasen 1090. Filtreringstrinnet kan utføres på de partisjonerte flatebitdataene for å eliminere feil sidepunkter eller punkter 820 som er plassert innenfor et smalt filtrerings område 830.
I trinn 1030, utvikles kontinuitetsbegrensninger for horisontflatebiter langs kunstig utvidede forkastningsblokkgrenser for å sikre kontinuiteten til overflaten og dens gradient over disse grenser. Kontinuitetsbegrensningene er lagret i databasen 1090.
I trinn 1040, anvendes flatebitdataene og de beregnede kontinuitetsbegrensningene for å utvikle en første, ikke forkastet horisontmodell på innsiden av hver blokk. Således er en ikke forkastet flatebitmodell av innsiden av hver forkastet blokk bygget med anvendelse av det tilsvarende parti av innmatingsdata og kontinuitetsbegrensninger. Det beregnes også en overflateflatebitgrense, f.eks. trimmekurver for en flatebit, ved å skjære den ikke forkastede flatebitmodell med grenseoverflatene til forkastningsblokken. De ikke-forkastede flatebitmodeller og den opprinnelige trimmekurven er lagret i databasen 1090.
I trinn 1050, anvendes flatebitdataene, kontinuitetsbegrensningsdataene, og flatebitgrensen for å beregne opprinnelig forkastningslokaliseringer for innsiden av hver flatebit med anvendelse av algoritmer beskrevet i US patentene som det er vist til tidligere. Således vil de opprinnelige forkastningsbeliggenhetene oppnås ved å skjære den ikke-forkastede flatebitmodell i trinn 1040 for hver forkastningsblokk med forkastningsoverflatene. Deretter filtreres og filtreringsavstanden økes adaptivt for hver forkastningslokalisering. Det gjentas helt til den maksimale distorsjons vinkel for hver forkastningslokalisering ligger under terskelen. Trinn 1050 er fullført ved å filtrere dataene i forhold til de opprinnelige forkastningslokaliseringene for å oppdatere den ikke-forkastede flatebitmodellen. Oppdater overflateflatebitgrensen, f.eks. flatebittrimmekurver ved å beregne skjæringen mellom den ikke-forkastede innsidehorisontmodell som resulterer fra den opprinnelige forkastningslokaliseringsberegning med grenseoverflatene til blokken. De opprinnelige forkastningslokaliseringene lagres, sammen med filtrerte flatebitdata, den oppdaterte ikke-forkastede modell og den oppdaterte flatebittrimmekurve i databasen 1090.
I trinn 1060, anvend de filtrerte flatebitdata, kontinuitetsbegrensningsdataene, overflateflatebitgrensen og de opprinnelige forkastningslokaliseringene fra trinn 1050 for å beregne en opprinnelig forkastet horisontmodell for innsiden av hver flatebit med anvendelse av algoritmen beskrevet i de tidligere refererte US patentene. Oppdater flatebittrimmekurven. Lagre den opprinnelige forkastede horisontmodell og oppdater flatebittrinnkurve i databasen 1090.
I trinn 1070, anvend den opprinnelige forkastede horisontmodell for innsiden av hver flatebit fra trinn 1060 for å beregne endelige forkastningslokaliseringer for innsideforkastningene med anvendelse av algoritmen beskrevet i US patenter som det er vist til tidligere. Lagre endelig innsideflatebitforkastningslokaliseringer for hver forkastningsblokk i databasen 1090.
I trinn 1080, anvend den endelige innsideforkastningslokalisering fra trinn 1070, flatebitdataene og kontinuitetsbegrensningene for å beregne en endelig forkastet horisontmodell for innsiden av hver flatebit med anvendelse av algoritmen beskrevet i US patenter som det er vist til tidligere. Oppdater overflateflatebitgrensen ved å beregne skjæringen av horisontmodellene som således er beregnet med grenseflatene for blokken. Lagre endelig forkastet flatebitmodell og oppdatert flatebittrinnkurve i databasen 1090.
Det tredimensjonale horisontmodelleringssystemet beskrevet ovenfor vil også støtte konformhorisontmodellering som beskrevet i de tidligere henviste US patentene. Det nye elementet her er at utformingsdataene er avledet i forbindelse med forkastningsblokkpartisjoneringen av reservoaret, ved trinn 510. De to hovedversjonene av konformmodellering er støttet: 1) med anvendelse av en begrenset gruppe med innmatingsdatapunkter som tilveiebringer den generelle lokaliseringen av modellen og en referansehorisontlappmodell som tilveiebringer utformingen, og 2) med anvendelse av en tykkelsesmodell og en referanselappmodell. Tykkelsesmodellen kan være enten en global, enkelt verdimodell for hele horisonten eller en gruppe av fiatelapptykkelsesmodeller. En enkelt modell er anvendelig når de geologiske deformasjonene ikke fører til forskjellig bergkompresjon på begge sider av en gitt oppskyvning. Flatebittykkelsesmodellen er anvendelig når de geologiske deformasjonene forårsaker forskjellig bergkompresjon over en oppskyvning. Det er opplagt at dette fenomen ikke kan representeres nøyaktig med en enkelt verditykkelsesmodell. Således vil en flerlappsmodell være vel egnet for slike anvendelser.
Kombinasjonen av tykkelse og referansemodeller er anvendelig ved å beregne en horisontmodell med spesifikke ønskede egenskaper. Imidlertid kan definisjonene av en god tykkelsesmodell være vanskelig. Systemet tillater at brukeren starter med tilnærming 1) ovenfra. Tykkelsesmodeliene er beregnet som et biprodukt av modelleringsprosessen. Når dette er fullført, kan systemets bruker endre tykkelsesmodellene som ønsket og anvende tilnærming 2) ovenfra.
Med henvisning nå til fig. 11, er et datamaskinsystem 1110 vist som er generelt anvendelig for de forskjellige utførelsene beskrevet. Systemet 1110 omfatter en prosessor 1115, et flyktig minne 1120, f.eks. RAM, et tastatur 1125, en pekeanordning 1130, f.eks. en mus, et ikke-flyktig minne 1135, f.eks. ROM, harddisk, floppy disk, CD-ROM og DVD, og fremvisningsanordning 1105 med en fremvisningsskjerm. Minnet 1120 og 1135 er for å lagre programinstruksjoner, som kan utføres av prosessor 1115 for å implementere forskjellige utførelser av en fremgangsmåte ifølge den foreliggende oppfinnelsen. Minnene 1120 og 1135 kan anvendes for å lagre partier av databasen 1090. Komponenter inkludert i systemet 1110 er sammenkoblet ved bussen 1140. En kommunikasjonsanordning (ikke vist) kan også kobles til bussen 1140 for å tillate informasjonsutveksling mellom system 1110 og andre anordninger slik som andre datamaskinsystemer via et nettverk slik som internett.
I flere utførelser inntar system 1110 flere fasonger, omfattende et personlig datamaskinsystem, klientserversystem, hovedramme-datamaskinsystem, parallellprosesseringsdatamaskinsystem, arbeidsstasjon, internettapplikasjon, PDA, innskutt prosessor med minne, osv. Det vil si at det bør forstås at begrepet datamaskinsystem er ment til å omfatte en hvilken som helst anordning med en prosessor som utfører instruksjoner fra et minnemedium.
Minnemediet vil fortrinnsvis lagre instruksjoner (også kjent som softwareprogram) for å implementere flere utførelser av en fremgangsmåte ifølge den foreliggende oppfinnelsen. I flere utførelser vil det ene eller flere softwareprogrammer være implementert på forskjellige måter, omfattende prosedyrebaserte teknikker, komponentbaserte teknikker, og/eller objektorienterte teknikker, blant annet. Spesifikke eksempler omfatter XML, C, C++, Java og Microsoft Foundation Classes
(MFC).
Med henvisning til fig. 12, er et klient-serverdatamaskinsystem 1230 anvendt for å implementere en utførelse av den foreliggende oppfinnelsen. Systemet modellerer flerverdihorisontmodeller i sammenhengen for en jordmodell, hvor kompleksitetene ved de forskjellige modellelementene i oppskyvningsstrukturer er skjult under enkelte føringer til komponenter i jordmodellen som er lagret på databasen 1090. Et softwaredesignmønster som kan påføres effektivt her er klient-servermodellen. Serveren 1210 er ansvarlig for å opprettholde og vedlikeholde databasens 1090 lagring og betjene anmodninger for datatilgang og lagring. Serveren 1210 er porten til det felles oppbevaringsstedet. De forskjellige komponenter (f.eks. forkastninger, horisonter osv.) i reservoaret 115 er representert i databasen 1090 og dataene for hver av dem er lagret med anvendelse av en spesiell datamodell. Serveren 1210 vet og håndhever forholdet mellom de forskjellige datamodellkomponenter. Jordmodelldataene er tilgjengelige for én eller flere klient 1220 datamaskinarbeidsstasjoner som kan samvirke med den via et nettverk 1250. Klient 1220 arbeidsstasjonen utfører et datamaskinprogram som virker med en spesiell del av jordmodellen og kan reagere overfor endringer i modellen utført av en annen klient gjennom databaseserveren 1210. Denne beskrivelsen er basert på angivelsen av US patentsøknad nr. 08/758,833, av Shamin Ahmed og Serge Dacic, og er innlemmet herved ved referanse. Datamaskinteknikken som beskrives tilveiebringer eller gir anledning for modulær softwaredesign og applikasjonsintegrering. F.eks. kan horisontmodelleringsfunksjonaliteten som beskrives i denne beskrivelsen utvikles til en applikasjon som er i stand til å samvirke med jordmodellen gjennom klient-serversystemet 1230. Den tar dataene for en horisont, forkastningsramme og overflater, og forkastningsblokkavdeling, ved trinn 510, og beregner og anbringer modellen på jordmodellen. En komplementær visualiseringsanvendelse kan utvikles på tilsvarende måte. Den kan, f.eks. vise frem på fremvisnings anordningen 1105 den beregnede horisontmodell og også en hvilken som helst av innmatingsmodellene. Begge applikasjonene vil kommunisere gjennom den felles jordmodell som er bragt på databasen 1090 og uttrykker interessebegivenheter som er sendt fra en applikasjon til en annen via serveren 1210.
Med henvisning til fig. 13, er et parallelt prosesseringsdatamaskinsystem 1390 som omfatter arbeidsstasjoner 1310, 1320 og 1330 anvendt for å implementere en utførelse av den foreliggende oppfinnelsen. Teknikken beskrevet i fig. 10 er vel egnet for parallell prosessering. Oppgaven å beregne horisontmodellkomponentene for hver forkastningsblokk, f.eks. 1340, 1350 og 1360 er delegert til en annen arbeidsstasjon. Trinn 1010 fra denne algoritme kan beregnes på en enkelt arbeidsstasjon. Etter dette, vil hver arbeidsstasjon utføre trinn 1020-1080 fra disse algoritmer og lagre resultatene i databasen 1090. På denne måten, vil mengden tid som er nødvendig for å beregne den fullstendige modellen kunne reduseres med fordel. I tillegg, kan mye større modeller beregnes med denne teknikken. Behovet for å synkronisere mellom beregningsarbeidsstasjoner som utfører beregningstrinnet 1020 er redusert, siden trinnene 1020-1080 er uavhengige av hverandre og kan utføres parallelt. Hver arbeidsstasjon er i stand til å kommunisere med en hvilken som helst annen arbeidsstasjon, og databasen 1090 kobles til serveren 1210 via et kommunikasjonssystem 1380. Serveren 1210 som er en felles ressurs for datamaskinarbeidsstasjonene 1310, 1320 og 1330, kan oppleve økt belastning. Imidlertid er forholdet mellom I/O og beregningssyklusene for algoritmene fra fig. 10 i favør for beregningstiden. Således er den felles lenke med databasen 1090 ikke en seriøs flaskehals for ytelsen til parallellprosesseringssystemet.
Beskrivelsen av den foreliggende utførelsen har vært presentert med den hensyn å illustrere og er ikke ment å være fullstendig eller å begrense oppfinnelsen til utførelsen som er beskrevet. Flere endringer og variasjoner vil være opplagt for fagmannen på området. For å si det på nytt, var utførelsene valgt og beskrevet for å forklare på best måte prinsippene bak oppfinnelsen, den praktiske applikasjon, og for å gjøre andre fagmenn på området i stand til å forstå oppfinnelsen. Flere andre utførelser med flere modifikasjoner kan være egnet for en spesiell anvendelse som man tenker på, men kan være innenfor rammen for den foreliggende oppfinnelsen. Fagmannen på området vil forstå at hardware og metodene som illustreres her kan variere avhengig av implementeringen. I tillegg er det viktig å merke at mens den foreliggende oppfinnelsen har vært beskrevet i forbindelse med et datamaskinsystem med en prosessor og et minne, vil fagmannen på området kunne skjønne at prosessene ved den foreliggende oppfinnelsen er i stand til å fordeles som datamaskinlesbart medium for instruksjoner på mange forskjellige måter og at den foreliggende oppfinnelsen vil gjelde på samme måte uansett den spesielle type signal der mediet som er anvendt i virkeligheten for å utføre fordelingen. Eksempler på datamaskinlesbare media omfatter registrerbare media f.eks. floppy disk, harddiskdrive, RAM, CD-ROM, CD-RW, DVD og overføringstypemedium slik som digitale og analoge kommunikasjonslenker.
For å gjenta, er flere ytterligere aspekter, modifikasjoner og endringer også tenkt og ment til å omfattes innenfor rammen for de følgende krav. I tillegg bør man forstå at i de følgende krav, er ikke handlingene nødvendigvis utført i den spesielle rekkefølgen som de er beskrevet i.
Claims (25)
1. Fremgangsmåte for modellering av et reservoar (115) med en horisontflate med en oppskyvning (220),
karakterisert vedat fremgangsmåten omfatter: a) å avdele et volum (135) av reservoaret (115) i forkastningsblokker (615, 625 og 635), hvor en grenseflate (645, 655) for hver av forkastningsblokkene (615, 625 og 635) omfatter et parti av oppskyvningen, hvor en innside av hver av forkastningsblokkene (615, 625, 635) omfatter bare enkeltverdihorisontflater; b) å avdele horisontdata som beskriver reservoaret (115) i datagrupper, hvor datagruppene definerer flatebiter (710 og 720), hvor hver flatebit tilsvarer en forkastningsblokk, hvor fiatebitene (710 og 720) omfatter en oppoverrettet (475) og nedoverrettet (485) forkastningstrase på motsatte sider av oppskyvningen (220); og c) å utvikle kontinuitetsbegrensninger (905) langs kunstig utvidede forkastningsblokkgrenser (680), hvor kontinuitetsbegrensningene (905) tilsvarer endepunkter (610 og 620) i den oppoverrettede (475) og nedoverrettede (485) til forkastningstrase.
2. Fremgangsmåte ifølge krav 1, a) omfattende kunstig å utvide forkastningsoverflaten (680) i oppskyvninger (220) til en grense i volum (135).
3. Fremgangsmåte ifølge krav 1, hvor forkastningsblokkene (615, 625 og 635) er formet ved kunstig å utvide forkastningsoverflater (680) i oppskyvningen (220) til en oppskjæring med en annen horisont med en annen forkastning.
4. Fremgangsmåte ifølge krav 1, hvor hvert datapunkt inkludert i horisontdataene er tildelt til bare én forkastningsblokk inkludert i forkastningsblokkene (615, 625 og 635).
5. Fremgangsmåte ifølge krav 1, a) omfattende å velge en undergruppe av alle forkastningene i reservoaret (115) hvor undergruppen omfatter oppskyvningen (220).
6. Fremgangsmåte ifølge krav 1, omfattende: å beregne en opprinnelig ikke-forkastet horisontmodell på innsiden av hver forkastningsblokk som svar til mottagelse av flatebitdata, og data som beskriver kontinuitetsbegrensningene (905); og å beregne en overflateflatebitgrense ved å skjære den opprinnelige ikke-forkastede horisontmodell med grenseflatene (645 og 655) i forkastningsblokken.
7. Fremgangsmåte ifølge krav 6, omfattende: å beregne opprinnelig forkastningslokaliseringer for innsiden av hver flatebit som svar til mottagelse av flatebitdata, kontinuitetsbegrensningen og overflateflatebitgrensen; og å oppdatere overflateflatebitgrensen ved å skjære den opprinnelige ikke-forkastede horisontmodell med grenseflater (645 og 655) i forkastningsblokken.
8. Fremgangsmåte ifølge krav 7, omfattende å beregne en opprinnelig forkastet modell for innsiden av hver forkastningsflatebit som svarer til mottagelse av flatebitdata, kontinuitetsbegrensningsdata, overflateflatebitgrense og de opprinnelige forkastningslokaliseringer.
9. Fremgangsmåte ifølge krav 8, omfattende å beregne endelige forkastningslokaliseringer for innvendige forkastninger som svar til den opprinnelige forkastede modell for innsiden av hver forkastningsflatebit.
10. Fremgangsmåte ifølge krav 9, omfattende: å beregne en endelig forkastet modell for innsiden av hver forkastningsflatebit som svar til mottagelse av de endelige forkastingslokaliseringer, flatebitdataene, og kontinuitetsbegrensningsdataene; og å oppdatere overflateflatebitgrensen ved å skjære den endelige forkastede horisontmodell med grenseflatene (645 og 655) i den endelige blokken.
11. Fremgangsmåte ifølge krav 1, omfattende å filtrere horisontdataene for å fjerne datapunkter (820) som er på feil side, hvor filtreringen er om ønskelig utført før utvikling av kontinuitetsbegrensningene (905).
12. Fremgangsmåte ifølge krav 1, hvor modelleringen definerer en modell for reservoaret (115) hvor modellen er tilpasset for å motta horisontdata og forkastningsoverflatene og forholdsdata.
13. Datamaskinprogramprodukt for et modelleringssystem som kan opereres for å tilveiebringe en representasjon av et reservoar (115) med en horisontflate med en oppskyvning (220),
karakterisert vedat datamaskinprogramproduktet omfatter: instruksjoner for å avdele et volum (135) i reservoaret (115) til et flertall av forkastningsblokker (615, 625 og 635), hvor en grenseflate (645, 655) for hvert flertall av forkastningsblokkene (615, 625 og 635) omfatter et parti av forkastningsblokken (220), hvor hver av flertallet av forkastningsblokker (615, 625 og 635) har en innside som er beskrevet av en tilsvarende enkel horisontverdimodell; instruksjoner for å avdele horisontdata som beskriver reservoaret (115) i datagrupper, hvor datagruppene definerer flatebiter (710 og 720), hvor hver flatebit tilsvarer en forkastningsblokk; og instruksjoner for å utvikle kontinuitetsbegrensninger (905) langs kunstig utvidede forkastningsblokkgrenser (680), hvor kontinuitetsbegrensningene (905) resulterer i tilsvarende endepunkter (610 og 620) for oppoverrettet (475) og nedoverrettet (485) traser for oppskyvningen (220).
14. Datamaskinprodukt ifølge krav 13, a) omfattende kunstig å utvide forkastningsoverflater (680) til oppskyvningen (220) til en grense i volum (135).
15. Datamaskinprogramprodukt ifølge krav 13, hvor forkastningsblokkene (615, 625 og 635) er dannet ved kunstig å utvide forkastningsoverflater (680) i oppskyvningen (220) til en skjæring av en annen horisont med en annen forkastning.
16. Datamaskinprogramprodukt ifølge krav 13, hvor hvert datapunkt inkludert i horisontdataene er tildelt til bare én forkastningsblokk inkludert i forkastningsblokkene (615, 625 og 635).
17. Datamaskinprogramprodukt ifølge krav 13, a) omfattende å velge en undergruppe av alle forkastninger i reservoaret (115), hvor undergruppen omfatter oppskyvningen (220).
18. Datamaskinprogramprodukt ifølge krav 13, omfattende: å beregne en opprinnelig ikke-forkastet horisontmodell på innsiden av hver forkastningsblokk som svar til mottagelse av flatebitdata, og data som beskriver kontinuitetsbegrensninger (905); og å beregne en overflatebitgrense ved å skjære den opprinnelige ikke-forkastede horisontmodell med grenseflatene (645 og 655) for forkastningsblokken.
19. Datamaskinprogramprodukt ifølge krav 18, omfattende: å beregne opprinnelige forkastningslokaliseringer for innsiden av hver flatebit som svar til mottagelse av flatebitdataene, kontinuitetsbegrensningsdataene og overflatebitgrensen; og å oppdatere overflateflatebitgrensen ved å skjære den opprinnelige ikke-forkastede horisontmodell med grenseflatene (645 og 655) i forkastningsblokken.
20. Datamaskinprogramprodukt ifølge krav 19, omfattende å beregne en opprinnelig forkastningsmodell for innsiden av hver forkastningsflatebit som svar til mottagelse av flatebitdata, kontinuitetsbegrensningsdata, overflateflatebitgrense og de opprinnelige forkastningslokaliseringer.
21. Datamaskinprogramprodukt ifølge krav 20, omfattende å beregne endelige forkastningslokaliseringer for innvendige forkastninger som svar til den opprinnelige forkastningsmodell for innsiden av hver forkastningsflatebit.
22. Datamaskinprogramprodukt ifølge krav 21, omfattende: å beregne en endelig forkastningsmodell for innsiden av hver forkastningsflatebit som svar til mottagelse av de endelige forkastningslokaliseringer, flatebitdata og kontinuitetsbegrensningsdata; og å oppdatere overflateflatebitgrensen ved å skjære den endelige forkastede horisontmodell med grenseflatene (645 og 655) for forkastningsblokken.
23. Datamaskinprogramprodukt ifølge krav 13, omfattende å filtrere horisontdataene for å fjerne datapunkter (820) som er på feil side, hvor filtreringer er om ønskelig utført før utvikling av kontinuitetsbegrensningene (905).
24. Datamaskinprogramprodukt ifølge krav 13, hvor modelleringen definerer en modell for reservoaret (115), hvor modellen er tilpasset til å motta horisontdata og forkastningsoverflater og forholdsdata.
25. Datamaskinprogramprodukt ifølge ett av kravene 13-24, i et modelleringssystem som kan opereres for å tilveiebringe representasjon av et reservoar (115) med en horisontflate med en oppskyvning (220),karakterisert vedat systemet omfatter: en prosessor (115); og et minne (1120, 1135) som lagrer instruksjoner som kan utføres med prosessoren (1115), hvor instruksjonene i henhold til krav 25 er utført,
hvor innsiden av hver forkastningsblokk (615, 625 og 635) omfatter bare enkeltverdihorisontflater og
hvor fiatebitene (710 og 720) omfatter en oppoverrettet (475) og nedoverrettet (485) forkastningstrase på motsatte sider av oppskyvningen (220).
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US29638501P | 2001-06-06 | 2001-06-06 | |
PCT/US2002/018095 WO2002099464A1 (en) | 2001-06-06 | 2002-06-05 | Automated system for modeling faulted multi-valued horizons |
Publications (2)
Publication Number | Publication Date |
---|---|
NO20035446D0 NO20035446D0 (no) | 2003-12-05 |
NO335854B1 true NO335854B1 (no) | 2015-03-09 |
Family
ID=23141785
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
NO20035446A NO335854B1 (no) | 2001-06-06 | 2003-12-05 | Automatisert system for modellering av flerverdihorisonter med forkastninger |
Country Status (7)
Country | Link |
---|---|
US (1) | US7424415B2 (no) |
EP (1) | EP1397703B1 (no) |
AT (1) | ATE387635T1 (no) |
CA (1) | CA2449286C (no) |
DE (1) | DE60225277D1 (no) |
NO (1) | NO335854B1 (no) |
WO (1) | WO2002099464A1 (no) |
Families Citing this family (49)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6853922B2 (en) * | 2001-07-20 | 2005-02-08 | Tracy Joseph Stark | System for information extraction from geologic time volumes |
US7899657B2 (en) * | 2003-01-24 | 2011-03-01 | Rockwell Automoation Technologies, Inc. | Modeling in-situ reservoirs with derivative constraints |
US7283911B2 (en) * | 2004-10-22 | 2007-10-16 | Landmark Graphics Corporation | System and method for interpreting reverse faults and multiple z-valued seismic horizons |
DE602004015297D1 (de) * | 2004-10-26 | 2008-09-04 | Total Sa | Verfahren und Computerprogramm zur Fehleroberflächenkonstruktion |
US8776895B2 (en) * | 2006-03-02 | 2014-07-15 | Exxonmobil Upstream Research Company | Method for quantifying reservoir connectivity using fluid travel times |
US9074454B2 (en) | 2008-01-15 | 2015-07-07 | Schlumberger Technology Corporation | Dynamic reservoir engineering |
US9372943B2 (en) * | 2008-05-05 | 2016-06-21 | Exxonmobil Upstream Research Company | Modeling dynamic systems by visualizing and narrowing a parameter space |
US7916576B2 (en) | 2008-07-16 | 2011-03-29 | Westerngeco L.L.C. | Optimizing a seismic survey for source separation |
AU2009314461B2 (en) | 2008-11-14 | 2015-11-05 | Exxonmobil Upstream Research Company | Forming a model of a subsurface region |
WO2010104535A1 (en) | 2009-03-13 | 2010-09-16 | Exxonmobil Upstream Research Company | Method for predicting fluid flow |
US8395966B2 (en) * | 2009-04-24 | 2013-03-12 | Westerngeco L.L.C. | Separating seismic signals produced by interfering seismic sources |
US8600708B1 (en) | 2009-06-01 | 2013-12-03 | Paradigm Sciences Ltd. | Systems and processes for building multiple equiprobable coherent geometrical models of the subsurface |
US8275593B2 (en) * | 2009-07-16 | 2012-09-25 | University Of Regina | Reservoir modeling method |
US8743115B1 (en) | 2009-10-23 | 2014-06-03 | Paradigm Sciences Ltd. | Systems and methods for coordinated editing of seismic data in dual model |
EP2317348B1 (en) * | 2009-10-30 | 2014-05-21 | Services Pétroliers Schlumberger | Method for building a depositional space corresponding to a geological domain |
WO2011066021A1 (en) | 2009-11-30 | 2011-06-03 | Exxonmobil Upstream Research Company | Adaptive newton's method for reservoir simulation |
IN2012DN05167A (no) | 2010-02-12 | 2015-10-23 | Exxonmobil Upstream Res Co | |
WO2011136861A1 (en) | 2010-04-30 | 2011-11-03 | Exxonmobil Upstream Research Company | Method and system for finite volume simulation of flow |
US20110310101A1 (en) * | 2010-06-22 | 2011-12-22 | Schlumberger Technology Corporation | Pillar grid conversion |
FR2962835B1 (fr) * | 2010-07-16 | 2013-07-12 | IFP Energies Nouvelles | Methode pour generer un maillage hexa-dominant d'un bassin geometriquement complexe |
AU2011283190A1 (en) | 2010-07-29 | 2013-02-07 | Exxonmobil Upstream Research Company | Methods and systems for machine-learning based simulation of flow |
AU2011283196B2 (en) | 2010-07-29 | 2014-07-31 | Exxonmobil Upstream Research Company | Method and system for reservoir modeling |
AU2011283193B2 (en) | 2010-07-29 | 2014-07-17 | Exxonmobil Upstream Research Company | Methods and systems for machine-learning based simulation of flow |
GB2502432B (en) | 2010-09-20 | 2018-08-01 | Exxonmobil Upstream Res Co | Flexible and adaptive formulations for complex reservoir simulations |
AU2011332274B2 (en) | 2010-11-23 | 2017-02-23 | Exxonmobil Upstream Research Company | Variable discretization method for flow simulation on complex geological models |
US9534485B2 (en) | 2011-04-18 | 2017-01-03 | Halliburton Energy Services, Inc. | Method for real-time downhole processing and detection of bed boundary for geosteering application |
US8688426B2 (en) | 2011-08-02 | 2014-04-01 | Saudi Arabian Oil Company | Methods for performing a fully automated workflow for well performance model creation and calibration |
US8731892B2 (en) | 2011-08-02 | 2014-05-20 | Saudi Arabian Oil Company | Systems and program product for performing a fully automated workflow for well performance model creation and calibration |
WO2013039606A1 (en) | 2011-09-15 | 2013-03-21 | Exxonmobil Upstream Research Company | Optimized matrix and vector operations in instruction limited algorithms that perform eos calculations |
US9759826B2 (en) * | 2012-04-03 | 2017-09-12 | Paradigm Sciences Ltd. | System and method for generating an implicit model of geological horizons |
CA2883169C (en) | 2012-09-28 | 2021-06-15 | Exxonmobil Upstream Research Company | Fault removal in geological models |
DE112013007652T5 (de) * | 2013-11-26 | 2016-08-11 | Landmark Graphics Corporation | Sequenz zum Beseitigen von Verwerfungen beim Aufbauen eines globalen Rasters für komplexe Verwerfungsnetztopologien |
WO2016018723A1 (en) | 2014-07-30 | 2016-02-04 | Exxonmobil Upstream Research Company | Method for volumetric grid generation in a domain with heterogeneous material properties |
US10359523B2 (en) | 2014-08-05 | 2019-07-23 | Exxonmobil Upstream Research Company | Exploration and extraction method and system for hydrocarbons |
EP3213127A1 (en) | 2014-10-31 | 2017-09-06 | Exxonmobil Upstream Research Company | Managing discontinuities in geologic models |
US10803534B2 (en) | 2014-10-31 | 2020-10-13 | Exxonmobil Upstream Research Company | Handling domain discontinuity with the help of grid optimization techniques |
AU2015339883B2 (en) | 2014-10-31 | 2018-03-29 | Exxonmobil Upstream Research Company | Methods to handle discontinuity in constructing design space for faulted subsurface model using moving least squares |
US10215886B2 (en) | 2016-02-19 | 2019-02-26 | Baker Hughes, A Ge Company, Llc | Plane-surface intersection algorithm with consistent boundary support |
US10466388B2 (en) | 2016-09-07 | 2019-11-05 | Emerson Paradigm Holding Llc | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles |
FR3058798A1 (fr) * | 2016-11-17 | 2018-05-18 | Landmark Graphics Corporation | Interpretation et visualisation d'horizon a valeurs z multiples a partir de donnees sysmiques |
AU2017360852A1 (en) * | 2016-11-17 | 2019-03-28 | Landmark Graphics Corporation | Multi-Z horizon interpretation and editing within seismic data |
HUE064459T2 (hu) | 2016-12-23 | 2024-03-28 | Exxonmobil Technology & Engineering Company | Eljárás és rendszer stabil és hatékony tározó szimulációhoz stabilitási proxyk alkalmazásával |
US10209380B2 (en) * | 2016-12-30 | 2019-02-19 | Schlumberger Technology Corporation | Methods and systems for juxtaposition across geological discontinuities |
CN109519169A (zh) * | 2018-10-25 | 2019-03-26 | 广东石油化工学院 | 复杂断块油藏分类评价方法 |
US10520644B1 (en) | 2019-01-10 | 2019-12-31 | Emerson Paradigm Holding Llc | Imaging a subsurface geological model at a past intermediate restoration time |
US11156744B2 (en) | 2019-01-10 | 2021-10-26 | Emerson Paradigm Holding Llc | Imaging a subsurface geological model at a past intermediate restoration time |
CN112505758B (zh) * | 2020-11-17 | 2021-10-26 | 中国石油集团工程咨询有限责任公司 | 一种基于断层分叉结构模型的复杂地质构造地震断层图像处理方法 |
CN114519784B (zh) * | 2022-02-14 | 2023-11-24 | 煤炭科学研究总院有限公司 | 三维断层的建模方法及其装置 |
WO2024063797A1 (en) * | 2022-09-19 | 2024-03-28 | Schlumberger Technology Corporation | Automated machine learning fault modeling with grouping |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5740342A (en) | 1995-04-05 | 1998-04-14 | Western Atlas International, Inc. | Method for generating a three-dimensional, locally-unstructured hybrid grid for sloping faults |
US6014343A (en) * | 1996-10-31 | 2000-01-11 | Geoquest | Automatic non-artificially extended fault surface based horizon modeling system |
-
2002
- 2002-06-05 US US10/480,099 patent/US7424415B2/en not_active Expired - Fee Related
- 2002-06-05 AT AT02741905T patent/ATE387635T1/de not_active IP Right Cessation
- 2002-06-05 DE DE60225277T patent/DE60225277D1/de not_active Expired - Lifetime
- 2002-06-05 CA CA002449286A patent/CA2449286C/en not_active Expired - Fee Related
- 2002-06-05 WO PCT/US2002/018095 patent/WO2002099464A1/en active IP Right Grant
- 2002-06-05 EP EP02741905A patent/EP1397703B1/en not_active Expired - Lifetime
-
2003
- 2003-12-05 NO NO20035446A patent/NO335854B1/no not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
CA2449286A1 (en) | 2002-12-12 |
US7424415B2 (en) | 2008-09-09 |
NO20035446D0 (no) | 2003-12-05 |
WO2002099464A1 (en) | 2002-12-12 |
EP1397703A1 (en) | 2004-03-17 |
ATE387635T1 (de) | 2008-03-15 |
EP1397703B1 (en) | 2008-02-27 |
DE60225277D1 (de) | 2008-04-10 |
US20040193960A1 (en) | 2004-09-30 |
CA2449286C (en) | 2007-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
NO335854B1 (no) | Automatisert system for modellering av flerverdihorisonter med forkastninger | |
US10995592B2 (en) | Method and system for analyzing the uncertainty of subsurface model | |
EP3329307B1 (en) | Assignment of systems tracts | |
EP2653893B1 (en) | Faulted geological structures containing unconformities | |
Cherpeau et al. | Stochastic structural modelling in sparse data situations | |
US10670753B2 (en) | History matching of time-lapse crosswell data using ensemble kalman filtering | |
RU2600944C1 (ru) | Формирование моделей распознаваемых геологических структур на основании набора узловых точек | |
NO335459B1 (no) | Modellering av geologiske objekter i undergrunnssoner med flere forkastninger | |
NO328482B1 (no) | Fremgangsmate, produkt og system for automatisert, ikke-kunstig, tredimensjonal forkastningsoverflate-basert horisontmodellering av undergrunnen | |
CA2724113A1 (en) | Virtual petroleum system | |
CA2371410C (en) | Method and apparatus for mapping uncertainty and generating a cube based on conditional simulation of random variables | |
EP3929630B1 (en) | Computer implemented method for correcting a reservoir model of a reservoir geological formation based on seismic images | |
US11163079B2 (en) | Method for producing a geological vector model | |
Milliotte et al. | From seismic interpretation to reservoir model: an integrated study accounting for the structural complexity of the Vienna Basin using an unstructured reservoir grid | |
CN112649893B (zh) | 一种面向薄储层的多资料多参数融合建模方法及系统 | |
Horna et al. | Isometric Unfolding of Stratigraphic Grid Units for Accurate Property Populating–Mathematical Concepts | |
Crane | Approach and application of industry software to structural investigations in the subsurface of Mercury's thrust fault-related landforms | |
Becerra et al. | Petroleum reservoir uncertainty mitigation through the integration with production history matching | |
Hinsch et al. | Uncertainties of balanced sections and the role of basement in the southern Fars area, Iran | |
CN106291748A (zh) | 一种获取烃源岩有机碳含量的方法 | |
CN113051697A (zh) | 古气藏地层孔隙压力演化史曲线生成方法及装置 | |
CN115469361B (zh) | 一种碎屑岩地层三维地质建模方法 | |
Suter | Novel Principles for Effective Earth Model Grid Management while Geosteering | |
Zhuang et al. | Estimation of the volume of a rock layer using kriging method for energy efficient excavation | |
Lacaze et al. | Fault analysis based on innovative watertight model method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM1K | Lapsed by not paying the annual fees |