Canada has the second largest proven oil reserves next to Saudi Arabia which is mostly located in Alberta and Saskatchewan but is unconventional heavy oil and bitumen. The tar sands are found at the surface and are mined, yet 80% of the 173 billion barrels of heavy oil and bitumen exist in-situ according to the Canadian Association of Petroleum Producers (CAPP). Two factors inhibit the economic extraction and processing of Canadian heavy oil; its enormous viscosity ranging from 1000 to over 1 million mPa.s and the asphaltene content (high molecular weight molecules containing heavy metals and sulphur). Heavy oil and bitumen were only included in the reserves estimates through the efforts of Canadian enhanced oil recovery (EOR) research. Viscosity reduction is the one common element of in-situ methods of heavy oil recovery with the exception of cold production. Currently, steam assisted gravity drainage (SAGD) and cyclic steam stimulation (CSS) are being used commercially in the field where the oil’s viscosity is reduced by injecting steam. Thermal methods are energy intensive requiring vast volumes of water such that any improvement would be beneficial. Solvent extraction is one alternative requiring no water, the solvent is recoverable and reusable, and depending on the mode of operation the heavy oil is upgraded in-situ. Vapour Extraction (VAPEX) and enhanced solvent extraction (N-SolvTM) are two such methods. VAPEX and N-Solv reduce the bitumen’s viscosity via mass transfer and a combination of mass and heat transfer, respectively. A light hydrocarbon solvent (instead of steam) is injected into an upper horizontal well where the solvent mixes with the heavy oil, reduces its viscosity and allows the oil to drain under gravity to a bottom production well. The idea of using solvents for heavy oil extraction has been around since the 1970s and both VAPEX and N-Solv are patented processes. However, there is still much to be learned about how these processes physically work. Research to date has focused on varying system parameters (including model dimensions, permeability, heavy oil viscosity, solvent type and injection rate, etc.) to observe the effect on oil production from laboratory scale models. Based on an early mass balance model by Butler and Mokrys (1989) and an improvement by Das (1995), molecular diffusion alone cannot account for the produced oil rates observed from laboratory models. Until recently, very little progress had been made towards qualifying and quantifying the mass transfer mechanisms with the exception of the diffusivity of light hydrocarbons in heavy oil. Mass transfer can only be by diffusion and convection. Differentiating and quantifying the contribution of each is complex due to the nature and viscosity of the oil. The goal of this thesis is to investigate the mass transfer mechanisms during the solvent recovery of heavy oil. Quantifying the diffusion of light hydrocarbon solvents has been an active topic of research with limited success since the mid 1990’s. The experimental approach presented here focused on capturing the rate of solvent mass transfer into the bitumen by measuring the bitumen swelling and the butane uptake independently. Unlike early pressure decay methods, the pressure is held constant to not violate the assumed equilibrium solvent concentration at the interfacial boundary condition. The high solubility of solvent in heavy oil complicates the physical modeling because simplifying assumptions of a constant diffusion coefficient, constant density and a quiescent liquid should not be used. The model was developed from first principles to predict the bitumen swelling. The form of the concentration dependent diffusivity was assumed and the diffusivity coefficients initially guessed. The swelling (moving boundary) was fixed by defining a new dimensionless space coordinate and the set of partial differential equations solved using the method of lines. Using the non-linear regression (lsqnonlin) function in MATLAB®, optimising for the difference in predicted and experimentally found bitumen heights and independently validating the result using the solvent uptake, the diffusivity of butane in heavy oil (at 25oC) was found to be Dsb = 4.78 x 10-6ωs + 4.91 x 10-6 cm2/s where ωs is the solvent mass fraction. Diffusion alone has proven inadequate in predicting oil recovery rates from laboratory scale models. It is logical to assume that convective mass transfer plays a role at mixing the solvent and bitumen while draining via gravity through the reservoir porous matrix. Solvent extraction experiments were conducted in etched glass micromodels to observe the pore scale phenomena. The pore scale mechanisms were found to differ depending on how the solvent extraction was operated, with non-condensing (VAPEX) or condensing (N-SolvTM) solvent. Observations show increased convective mixing and an increased rate of interface advancement when the solvent condenses on the bitumen surface. Evidence of trapped butane vapour being mobilised with the draining live oil and a technique of observing solvent extraction using UV light confirm that the draining live oil is on average one pore deep. While the interface appears from a distance to be uniform, at the pore scale it is not. Live oil can drain from one to two pores via capillary displacement mechanisms in one section of the interface and via film flow only in another area (James and Chatzis 2004; James et al. 2008). This work also shows the detrimental impact of having a non-condensable gas present during solvent extraction (James and Chatzis 2008). In summary, this work emphasises the mass transfer and drainage displacement mechanisms of non-condensing (VAPEX) and condensing (N-Solv) solvent extraction methods of heavy oil recovery.

Canada has the second largest proven oil reserves next to Saudi Arabia which is mostly located in Alberta and Saskatchewan but is unconventional heavy oil and bitumen. The tar sands are found at the surface and are mined, yet 80% of the 173 billion barrels of heavy oil and bitumen exist in-situ according to the Canadian Association of Petroleum Producers (CAPP). Two factors inhibit the economic extraction and processing of Canadian heavy oil; its enormous viscosity ranging from 1000 to over 1 million mPa.s and the asphaltene content (high molecular weight molecules containing heavy metals and sulphur). Heavy oil and bitumen were only included in the reserves estimates through the efforts of Canadian enhanced oil recovery (EOR) research. Viscosity reduction is the one common element of in-situ methods of heavy oil recovery with the exception of cold production. Currently, steam assisted gravity drainage (SAGD) and cyclic steam stimulation (CSS) are being used commercially in the field where the oil’s viscosity is reduced by injecting steam. Thermal methods are energy intensive requiring vast volumes of water such that any improvement would be beneficial. Solvent extraction is one alternative requiring no water, the solvent is recoverable and reusable, and depending on the mode of operation the heavy oil is upgraded in-situ. Vapour Extraction (VAPEX) and enhanced solvent extraction (N-SolvTM) are two such methods. VAPEX and N-Solv reduce the bitumen’s viscosity via mass transfer and a combination of mass and heat transfer, respectively. A light hydrocarbon solvent (instead of steam) is injected into an upper horizontal well where the solvent mixes with the heavy oil, reduces its viscosity and allows the oil to drain under gravity to a bottom production well. The idea of using solvents for heavy oil extraction has been around since the 1970s and both VAPEX and N-Solv are patented processes. However, there is still much to be learned about how these processes physically work. Research to date has focused on varying system parameters (including model dimensions, permeability, heavy oil viscosity, solvent type and injection rate, etc.) to observe the effect on oil production from laboratory scale models. Based on an early mass balance model by Butler and Mokrys (1989) and an improvement by Das (1995), molecular diffusion alone cannot account for the produced oil rates observed from laboratory models. Until recently, very little progress had been made towards qualifying and quantifying the mass transfer mechanisms with the exception of the diffusivity of light hydrocarbons in heavy oil. Mass transfer can only be by diffusion and convection. Differentiating and quantifying the contribution of each is complex due to the nature and viscosity of the oil. The goal of this thesis is to investigate the mass transfer mechanisms during the solvent recovery of heavy oil. Quantifying the diffusion of light hydrocarbon solvents has been an active topic of research with limited success since the mid 1990’s. The experimental approach presented here focused on capturing the rate of solvent mass transfer into the bitumen by measuring the bitumen swelling and the butane uptake independently. Unlike early pressure decay methods, the pressure is held constant to not violate the assumed equilibrium solvent concentration at the interfacial boundary condition. The high solubility of solvent in heavy oil complicates the physical modeling because simplifying assumptions of a constant diffusion coefficient, constant density and a quiescent liquid should not be used. The model was developed from first principles to predict the bitumen swelling. The form of the concentration dependent diffusivity was assumed and the diffusivity coefficients initially guessed. The swelling (moving boundary) was fixed by defining a new dimensionless space coordinate and the set of partial differential equations solved using the method of lines. Using the non-linear regression (lsqnonlin) function in MATLAB®, optimising for the difference in predicted and experimentally found bitumen heights and independently validating the result using the solvent uptake, the diffusivity of butane in heavy oil (at 25oC) was found to be Dsb = 4.78 x 10-6ωs + 4.91 x 10-6 cm2/s where ωs is the solvent mass fraction. Diffusion alone has proven inadequate in predicting oil recovery rates from laboratory scale models. It is logical to assume that convective mass transfer plays a role at mixing the solvent and bitumen while draining via gravity through the reservoir porous matrix. Solvent extraction experiments were conducted in etched glass micromodels to observe the pore scale phenomena. The pore scale mechanisms were found to differ depending on how the solvent extraction was operated, with non-condensing (VAPEX) or condensing (N-SolvTM) solvent. Observations show increased convective mixing and an increased rate of interface advancement when the solvent condenses on the bitumen surface. Evidence of trapped butane vapour being mobilised with the draining live oil and a technique of observing solvent extraction using UV light confirm that the draining live oil is on average one pore deep. While the interface appears from a distance to be uniform, at the pore scale it is not. Live oil can drain from one to two pores via capillary displacement mechanisms in one section of the interface and via film flow only in another area (James and Chatzis 2004; James et al. 2008). This work also shows the detrimental impact of having a non-condensable gas present during solvent extraction (James and Chatzis 2008). In summary, this work emphasises the mass transfer and drainage displacement mechanisms of non-condensing (VAPEX) and condensing (N-Solv) solvent extraction methods of heavy oil recovery.

Abstract The subject of the thesis was to analyze by an engineering approach the global CO₂ balance and CO₂ utilisation. The aim was to apply methods and knowledge used in engineering sciences to describe the global CO₂ balance and the role of CO₂ in anthropogenic utilisation applications. Moreover barriers restricting commercialisation of new applications are discussed. These subjects were studied by literature reviews and calculations based on thermodynamics models. Engineering methods have shown to be applicable to describe the global balance of CO₂ and to define by a numerical way the Earth’s system carrying capacity. Direct and indirect actions, which mitigate the overload situation, were derived from the results. To screen out the attractive CO₂ properties in utilisation applications a mapping analysis was carried out. Properties, which enhance mass and heat transfer, are one of the most meaningful characteristics from the chemical engineering point of view. Attractive properties are often achieved at the supercritical state. Engineering thermodynamic methods were used in fluid phase determination of the case studies. Even simple methods are sufficient to advice experimental research work. The thermodynamic knowledge is the basement in creation of industrial scale chemical processes. If detailed information on system properties is needed, a model development due to the special requirements of high pressure systems and CO₂ features is required. This knowledge covers property information from all the components involved in chemical reactions. In addition to engineering knowledge successful technology transfer requires positive social structure as well. Finally, if the humankind is willing to mimic Nature and use light of the Sun as an energy source in engineering systems, development of thermodynamic methods is required also in this area. Especially the work terms, originally defined in classical mechanical thermodynamics and afterwards formulised also in other parts of the engineering fields, play a key role. If this development work is successful, we may see the shift from thermodynamics approach to ‘photodynamics’. Mitigation of global warming is a problem, which needs several kinds of activities. As a result of this study, there are listed a few engineering actions, which have a possibility to contribute to the work towards the carbon neutral society. / Tiivistelmä Väitöskirjatyössä sovelletaan insinööritieteissä käytettyjä metodeja ja tietämystä määriteltäessä ilmakehän CO₂-tase sekä antropogeenisten hyötykäyttökohteiden merkitys teollisissa prosesseissa ja globaaleissa CO₂-virroissa. Lisäksi pohditaan uusien CO₂-hyötykäyttösovellusten kaupallistamiseen liittyviä rajoitteita. Näitä aiheita on tutkittu käymällä läpi tieteellistä kirjallisuutta ja tekemällä laskelmia. Insinööritieteistä tutun taselaskennan avulla tarkastellaan ilmakehän CO₂-virtoja. Sen pohjalta määritetään numeerisesti maapallon CO₂-kantokyky. Tuloksista johdetaan suoria ja epäsuoria toimenpide-ehdotuksia, joiden avulla voidaan lieventää ilmakehän CO₂-ylikuormaa. Kartoitusmenetelmän avulla selvitetään hyötykäytön kannalta edulliset CO₂:n aineominaisuudet. Kemiantekniikan näkökulmasta ominaisuudet, jotka parantavat aineen- ja lämmönsiirtoa, ovat kiinnostavimpia. Nämä ominaisuudet tulevat esille silloin, kun fluidi on ylikriittisessä olomuodossa. Termodynaamisia laskentamenetelmiä sovelletaan esimerkkiseosten olomuodon eli faasin määrityksessä. Tulokset osoittavat, että jopa verraten yksinkertaiset menetelmät antavat tietoja, jotka auttavat ymmärtämään laboratoriokokeiden faasikäyttäytymistä. Teollisen mittakaavan kemiallisten prosessien kehityksessä ja suunnittelussa termodynamiikan hallitseminen on keskeinen edellytys. Jos CO₂:n kiinnostavia ominaisuuksia toivotaan hyödynnettävän teollisesti, korkeapaineisten systeemien termodynaamisen teorian hallinta sekä aineominaisuuksien määrittäminen kaikille systeemiin osallistuville komponenteille ja niiden seoksille nousee merkittävään asemaan. Läpikotainen teorian ja teknisten perusteiden hallitseminen ei vielä takaa menestyksellistä teknologiansiirtoa pienestä suureen mittakaavaan. Lisäksi tarvitaan myönteinen ja kannustava yhteiskuntajärjestelmä. Mikäli tavoitellaan vielä rohkeampaa kehitysnäkymää, tilannetta, jossa luonnon tavoin CO₂-prosessien energianlähteenä käytettäisiin auringonvaloa, havaitaan, että tämäkin askel edellyttäisi termodynaamista menetelmäkehitystä. Keskeinen termodynaaminen konsepti on työ. Työ siirtää energiaa ympäristön ja systeemin välillä. Tämä on määritelty jo klassisessa mekaniikassa; kappaleen siirto tietystä paikasta toiseen. Kemiantekniikassa työlle on kehitetty käyttökelpoisia kaavoja paine–tilavuus–lämpötila-systeemeihin. Mikäli työn elementit kyettäisiin määrittelemään auringonvalon fotoenergialle, avaisi se uusia näkymiä reaktiokemiaan. Silloin termodynamiikan sijaan voitaisiin ehkä mieluummin puhua 'photodynamiikasta'. Ilmaston lämpeneminen on ongelma, jonka lieventämiseen tarvitaan useanlaisia toimia. Etsittäessä tietä kohti hiilineutraalia yhteiskuntaa, insinöörit voivat avustaa suunnan löytämisessä hyödyntämällä tieteenalallaan käytettyjä metodeja ja teorioita sekä tarpeen vaatiessa kehittää niitä edelleen uusille alueille.