Reservoarsimulering

Reservoarsimulering er bruk av datamaskiner til å løse et sett av matematiske ligninger som beskriver flyt av væsker (typisk gass, olje og vann) i et petroleumsreservoar.  
Slik simulering blir brukt i utstrakt grad av oljeselskap som ledd i arbeid med å forstå oppførselen til et reservoar og 
til å forutsi og planlegge framtidig produksjon fra reservoaret. Programvaren som benyttes kalles en reservoarsimulator.
Typiske resultat fra reservoarsimulering er produksjonsprofiler for reservoaret, det vil si simulerte verdier av produsert volum fra brønner og felt som funksjon av tid.  En reservoarsimulator vil også beskrive
utvikling av trykk og væskefordeling i reservoaret.
Reservoarsimulering undervises i høyere utdanning som fag både i matematikk og i petroleumsteknologi. I matematikk er fokus på utledning av ligninger og på effektive numeriske løsningsmetoder.
I petroleumsteknologi er fokus på bruk av reservoarsimulatorer i praktisk arbeid.
Programvare for reservoarsimulering er utviklet både som åpen kildekode til fri distribusjon, som kommersielt tilgjengelige applikasjoner og som intern programvare i oljeselskap.


== Bruk av reservoarsimulering ==
Før et felt kan bygges ut, stiller norske myndigheter krav om at utbygger leverer en plan for utbygging og drift.  Denne planen skal inneholde data for forventet produksjon fra reservoaret, og 
disse produksjonsestimatene er nesten alltid generert ved hjelp av reservoarsimulering.  Også for selskapsrapportering til børs stilles krav om informasjon om forventet produksjon. Både operatørselskap
og partnere vil derfor gjennomføre reservoarsimuleringer for så å si alle petroleumsforekomster, enten de er i drift eller de vurderes for framtidig produksjon.
Reservoarsimuleringer er også et viktig grunnlag for å vurdere alternative utbygningsløsninger og for å optimalisere produksjonen fra feltet.
I Norge vil også Oljedirektoratet gjennomføre egne reservoarsimuleringer: «Ved å foreta eget faglig arbeid med blant annet geologiske modeller og reservoarsimulering kan OD utfordre selskapenes vurderinger av ressurspotensialet».Bruk av reservoarsimulering på et gitt gass- eller oljefelt kan omfatte flere steg:  Modellbygging, historietilpasning og prediksjon.


=== Modellbygging ===
Før en simulering kan gjennomføres, må grunnleggende data for reservoaret defineres som inngangsdata til reservoarsimulatoren.  Dette omfatter mange typer data:

En beskrivelse av reservoargeometri, typisk basert på innsamlet data fra seismikk og brønnlogger.  Reservoaret beskrives som et lukket område, begrenset av horisonter og forkastninger.
Petrofysiske reservoaregenskaper, inkludert fordeling av porøsitet og permeabilitet.
Væske-egenskaper, som tetthet, kompressibilitet og viskositet.
Opprinnelig væskefordeling i reservoaret, definert ved dyp til væskekontakter (gass/olje og olje/vann) og ved likevekt mellom kapillærkrefter og tyngdekrefter.
Posisjon og konfigurasjon av brønner.Summen av all inngangsdata kalles for en simuleringsmodell for feltet, og en slik modell vil typisk være basert på en geologisk reservoarmodell.  Simuleringsmodeller blir vedlikeholdt
ved regelmessige oppdateringer, etter hvert som en samler erfaring med feltet og bruken av modellene.
En modell som omfatter hele petroleumsforekomsten kalles en fullfeltsmodell.  I motsetning til dette vil en sektormodell bare inkludere en del av reservoaret.  En brønnmodell kan lages for
å studere produksjonen fra en enkelt brønn eller et fåtall av brønner.


=== Historietilpasning ===
For felt som har vært i produksjon i en periode vil gyldigheten til en simuleringsmodell kunne prøves ut ved å se hvor godt den er i stand til å reprodusere observert feltoppførsel.  
Avvik mellom observasjoner og simulert resultat vil så å si alltid 
forekomme, fordi det er så stor usikkerhet knyttet til alle deler av reservoarbeskrivelsen. Historietilpasning er en arbeidsprosess der en forsøker å endre reservoarbeskrivelsen i modellen for å oppnå bedre samsvar
mellom observert og simulert resultat. En vellykket historietilpasning vil øke troverdigheten til modellen og redusere usikkerhet når modellen brukes til å forutsi framtidig reservoaroppførsel.
I historietilpasning er det viktig å ha en god beskrivelse av usikkerhet i modellen, det vil si en beskrivelse av hvilke deler av modellen en stoler mest og minst på.
Det eksisterer en rik litteratur på teknikker for historietilpasning, og også spesielle programvareverktøy for computer-assistert historietilpasning.


=== Prediksjon ===
Reservoarsimulering vil typisk bli brukt til å forutsi framtidig reservoaroppførsel, altså til å utføre en prediksjon.  
På grunn av usikkerhet i reservoarbeskrivelsen, vil arbeidet også omfatte forsøker på å estimere usikkerheten i prediksjonen.


== Matematisk modellering ==
Reservoarsimulering har mange fellestrekk med numerisk væskedynamikk, både i ligningsformulering og i løsningsmetoder.


=== Modelligninger ===
De matematiske ligningene som løses er sammensatt av tre typer ligninger:

Kontinuitetsligninger.  Dette er ligninger som uttrykker at masse er bevart, i reservoarsimulering massen av en komponent.  En komponent er en (tilnærmet) kjemisk uavhengig bestanddel i systemet.
Bevegelsesligninger.  I et porøst medium er væskebevegelsen beskrevet av Darcys lov.
Konstitutive relasjoner.  Dette er relasjoner mellom ulike material- eller væske-egenskaper, inkludert termodynamiske relasjoner.Kombinert vil disse ligningene lede til et sett av partielle differensialligninger som beskriver transporten av hver væskekomponent i reservoaret.  De ukjente variablene i ligningssettet er typisk
størrelser som beskriver væskefordeling, samt trykk, som funksjon av rom og tid.
Ofte antar en at reservoarprosessen er isotermal, det vil si at prosessen skjer uten temperaturendringer.  Dersom endringer i temperatur er viktig, så må ligningssettet utvides med en kontinuitetsligning 
som beskriver konservering av energi.


=== Løsningsmetoder ===
Ligningssettet blir typisk løst med en endelig differensmetode, selv om også alternative metoder er i bruk.  Reservoaret blir delt i et stort antall celler eller blokker, og hver celle har variable som 
beskriver væsketilstanden i cella. Inndelingen i celler gjøres som en del av modellbyggingen.
Arbeidet med å løse ligningssettet vil avhenge sterkt av hvor mange celler en deler reservoaret inn i. De første reservoarsimulatorene kunne håndtere bare et hundretall celler, men med utvikling av bedre matematiske
metoder og kraftigere datamaskiner, har modellene blitt større og større. En vanlig modell i dag vil ha et antall celler i størrelsesorden hundre-tusen til millioner, men det er gjennomført simuleringer med over 
en milliard celler.


== Simulatortyper ==
Reservoarsimulatorer kan klassifiseres som svartolje-simulatorer eller komposisjonelle simulatorer, basert
på beskrivelsen av væskene i reservoaret.  En kan også klassifisere simulatorer etter matematiske egenskaper.  Eldre simulatorer er gjerne spesialiserte, mens mange 
nyere simulatorer kombinerer mange typer inn i en og samme applikasjon.


=== Svartolje-simulatorer ===
I en svartolje-simulator (black oil-simulator) modelleres væskene i reservoaret som bestående av tre komponenter gass, olje og vann.  Gasskomponenten består av alle molekyler som er gass ved overflatevilkår, 
ved standard trykk og temperatur, og tilsvarende for de to andre komponentene, olje og vann.
Videre antar en at i reservoaret har tre væskefaser, som også kalles gass, olje og vann.  I reservoaret vil endringer i trykk føre til faseoverganger, for eksempel at olje fordamper til gass.  Dette håndteres i 
den matematiske formuleringen ved at gasskomponenten kan opptre både som gassfase og oljefase. I en faseovergang vil gasskomponenten da bevege seg fra gassfasen til oljefasen.
I en tradisjonell svartolje-formulering vil oljekomponenten alltid opptre i oljefasen og vannkomponenten alltid i vannfasen.  Faseoverganger påvirker altså kun fordeling av gasskomponenten i de to petroleumsfasene.
I en utvidet svartolje-formulering tillater en at oljekomponenten kan opptre både i oljefasen og i gassfasen.  En slik formulering er nødvendig for reservoar der trykkendringer fører til kondensering av gass til væske.
Svartolje-simuleringer er den vanligste formen for reservoarsimulering.


=== Komposisjonelle simulatorer ===
I en komposisjonell simulator brukes en beskrivelse der de involverte væskene deles inn i et antall uavhengige komponenter, typisk fem til 10 komponenter.  Hver komponent blir definert til å inneholde en gruppe av 
kjemiske komponenter.
Komposisjonelle simulatorer brukes for felt der en mer nøyaktig væskebeskrivelse er nødvendig, for eksempel i forbindelse med økt oljeutvinning (EOR) der spesielle kjemikalier er i bruk.


=== Termiske simulatorer ===
Termiske simulatorer brukes for reservoarprosesser der temperaturendringer er viktige, for eksempel til modellering av dampinjeksjon i tungolje-reservoar.


=== Strømlinjesimulatorer ===
Strømlinjesimulering er en alternativ teknikk til tradisjonell reservoarsimulering.  En strømlinje er en linje som er overalt tangent til væskehastighetsfeltet, og ved å utnytte disse kan en lage simulatorer
der modelligningene lar seg løse svært raskt. For å kunne oppnå dette må en imidlertid gjøre en del forenklinger i modelleringen av reservoarfysikken, slik at simuleringsresultatet kan være gyldig
under begrensede forutsetninger.


== Reservoarsimulatorer ==
Reservoarsimulatorer med åpen kildekode:

BOAST – Black Oil Applied Simulation Tool (Boast), fri svartolje-simulator fra U.S. Department of Energy. I hovedsak tiltenkt for utdanningsformål.
MRST – MATLAB Reservoir Simulation Toolbox (MRST), utviklet av SINTEF som en samling av verktøy i MATLAB. Verktøykasse for prototyping av matematiske metoder for reservoarsimulering.
OPM – Open Porous Media (OPM).  Svartolje-simulator utviklet i et samarbeid mellom NORCE, SINTEF Digital, Equinor, Ceetron Solutions og HPC-Simulation-Software.Kommersielle reservoarsimulatorer:

ECHELON - Svartolje-simulator fra Stone Ridge Technology.
ECLIPSE - Svartolje-simulator (ECLIPSE 100 og 200) og komposisjonell simulator (ECLIPSE 300) fra Schlumberger.  Opprinnelig utviklet av det engelske firmaet ECL (Exploration Consultants Limited), og simulatornavnet var et akronym for «ECL´s Implicit Program for Simulation Engineering».
GEM - Komposisjonell simulator fra Computer Modelling Group (CMG).
IMEX - Svartolje-simulator fra Computer Modelling Group (CMG).
INTERSECT - Kombinert svartolje, komposisjonell og termisk simulator fra Schlumberger.
Nexus - Kombinert svartolje, komposisjonell og termisk simulator fra Landmark, eid av Halliburton.
Sensor - Kombinert svartolje og komposisjonell simulator fra Coats Engineering Inc.
STARS - Termal simulator fra Computer Modelling Group (CMG), designet spesielt for avanserte kjemiske prosesser.
Tempest MORE - Kombinert svartolje, komposisjonell og termisk simulator fra Roxar, et selskap eid av Emerson Electric Company.
tNavigator - Kombinert svartolje, komposisjonell og termisk simulator fra Rock Fluid Dynamics.
3DSL - Strømlinjesimulator utviklet av STREAMsim.Reservoarsimulatorer utviklet av oljeselskap til internt bruk:

EMPower - Utviklet av Exxon Mobil.
GigaPOWERS - Utviklet av Saudi Aramco.
MoReS - Utviklet av Royal Dutch Shell.  Shell samarbeider også med Computer Modelling Group om utvikling av en ny simulator med navn CoFlow.
PSim - Utviklet av ConocoPhillips.I tillegg eksisterer det flere reservoarsimulatorer utviklet og brukt i akademiske miljø, til forskning og undervisningsformål.


== Sammenligningsstudier ==
Society of Petroleum Engineers har arrangert flere studier der resultat fra ulike reservoarsimulatorer har vært sammenlignet og vurdert, i såkalte «Compartive Solution Projects».

SPE 1 (1981):  Tofase svartolje-simulering i tre romdimensjoner.  Sju reservoarsimulatorer deltok, fra Amoco, CMG, Exxon, Intercomp, Mobil, Shell og Scientific Software.
SPE 2 (1986):  Trefase simulering av koningseffekter.  Elleve selskap deltok:  Arco, Chevron, D&S Research and Development, Franlab, Gulf, Harwell, Intercomp, McCord-Lewis Energy Services, J.S.Nolen and Associates, Scientific Software og Shell.
SPE 3 (1987):  Komposisjonell simulering av gass-sirkulering i et reservoar med gass-kondensat.  Ni simulatorer deltok, fra Arco, Chevron, Core Lab, CMG, Elf Aquitaine, Intercomp, Marathon,  McCord-Lewis Energy Services og Petek.
SPE 4 (1987):  Simulatorer for dampinjeksjon.  To-dimensjonal radial strøm. Seks simulatorer deltok.
SPE 5 (1987):  Komposisjonell simulering av blandbare væsker.  Sju simulatorer deltok, fra Arco, British Petroleum, CMG, Chevron, Energy Resource Consultant, Reservoir Simulation Research Corporation, samt Todd, Dietrich and Chase,
SPE 6 (1990):  Simulatorer for oppsprukne reservoar. Ti simulatorer deltok, fra Chevron, CMG, Dancomp, ECL, Franlab, Japan Oil, Marathon, Phillips Petroleum, Simulation and Modelling Consultency og Scientific Software-Intercomp,


== Historie ==
Med framveksten av de eksperimentelle datamaskinene på 1940- og 1950-tallet økte interessen for å ta slike maskiner i bruk også i petroleumsindustrien.    
Fram til ca 1960 ble imidlertid storparten av analyser for petroleumsreservoar utført ved håndregning eller ved hjelp av bord-kalkulatorer..  
De første forsøkene på bruk av datamaskiner var preget av sterke modell-forenklinger, som for eksempel bruk av modeller for to-fase-strøm av olje og vann i to romdimensjoner.  Problemene som ble løst var knyttet
til simulering av trykksenking (deplesjon) og trykkvedlikehold.Det første forsøket på å lage en modell i tre romdimensjoner ble presentert i 1967.  Modellen som ble simulert hadde 10 x 10 x 5 celler.  Programmet var skrevet i Fortran og ble kjørt på en 
IBM 7044 med en lagerkapasitet på 32 000 ord.
Utvikling av raskere og billigere datamaskiner førte til at stadig mer av beregningene ble utført på datamaskiner.  Utvikling av algoritmer og teknologi for reservoarsimulering har fulgt 
parallelt med utviklingen av datamaskinene, fra de første forsøk på å konstruere maskiner, til stormaskiner, til PC-er, til vektor-maskiner og parallell-arkitektur.
Det første kommersielle selskapet for reservoarsimulering var Intercomp Resource Development and Engineering, Inc, grunnlagt i 1968 av Keith Coats.  Selskapet ble i 1982 omdøpt til 
Scientific Software-Intercomp.
Arbeid med å lage de første reservoarsimulatorene kunne bygge på mye etablert teori for løsning av en enkelt partiell differensialligning.  Teorien måtte imidlertid generaliseres til det system av ligninger som 
ble brukt til å beskrive væskeflyt i et reservoar.  Mye tidlig arbeid var derfor fokusert på å finne ligningsformuleringer og numeriske metoder som kombinerte stabilitet og effektivitet.  
En IMPES-formulering, med en
implisitt trykkløsning og en eksplisitt metningsløsning, var lenge dominerende.  Etter hvert som datamaskinene ble mer effektive, ble en fullt implisitt formulering vanlig brukt.
I 1968 arrangerte Society of Petroleum Engineers det første symposiet innen numerisk reservoarsimulering (Symposium on Numerical Simulation of Reservoir Performance), i april i Dallas.
I det første tiåret hadde alle modellene en enkel bokslignende eller radiell geometri. På starten av 70-tallet kom de første forsøkene på å modellere en mer realistisk geometri, med bruk av kurvelineære koordinater 
og irregulære gitter.  Hjørnepunktsgeometri, som er standard i bruk i dag, ble introdusert først i 1992.Prisstigning på olje på 1970-tallet førte til at oljeselskapene tok i bruk mange nye utvinningsteknikker, og dette stilte nye krav til at reservoarsimulering måtte kunne modellere effekter knyttet til overflatespenning, 
komplekse fase-vekselvirkninger, CO2-injeksjon, damp-injeksjon og in-situ-forbrenning.  Dette førte til at en fikk utviklet mange nye simulatorer for spesielle formål.
Etterhvert som modellene økte i størrelse, ble det behov for bedre algoritmer for å løse det lineære ligningssystemet som utgjør kjernen i den numeriske håndteringen av modelligningene. 
Tidlige simulatorer brukte typisk en direkte
metode som Gauss-eliminasjon for løsningen, men ganske raskt overtok iterative metoder, som oftest basert på en form for konjugerte gradienter. Orthomin-metoden ble introdusert i 1976 og 
nøstet faktorisering i 1983.
På 1980-tallet og 1990-tallet utviklet brønnteknologien seg sterkt, med lengre og lengre brønner, vertikale, skrå-boret og horisontale. Enkle brønnmodeller basert på vertilkal hydrostatisk likevekt var ikke lenger tilstrekkelig, 
og en rekke nye brønnmodeller ble tatt i bruk for å kunne håndtere kryss-strøm, friksjonseffekter, avansert brønnstyring, grenbrønner og «smarte» brønner.  Simuleringsmodellen ble også i 
økende grad utvidet til å omfatte væskestrøm opp til overflaten og i prosesseringsutstyr på overflaten.


== Referanser ==


== Litteratur ==
K.Aziz, A.Settari (2002). Petroleum reservoir simulation. Calgary: Blitzprint Ltd. ISBN 0-9730614-0-5. J.R.Gilman, C.Ozgen (2013). Reservoir simulation: History matching and forecasting. Richardson, TX: Society of Petroleum Engineers. ISBN 978-1613992920. A.Datta-Gupta, M.J.King (2007). Streamline simulation: Theory and practice. Richardson, TX: Society of Petroleum Engineers. ISBN 978-1-55563-111-6.