No description available.

A short-duration production test was undertaken at the Mallik site in Canada’s Mackenzie Delta in April 2007 as part of the JOGMEC/NRCan/Aurora Mallik 2007 Gas Hydrate Production Research Well Program. Reservoir stimulation was achieved by depressurization of a concentrated gas hydrate interval between 1093 and 1105m (RKB). Geologic and porous media conditions of the production interval have been quantified by geophysical studies undertaken in 2007 and geophysical and core studies undertaken by previous international partnerships in 1998 and 2002. These investigations have documented that the production interval consists of a sand-dominated succession with occasional silty sand interbeds. Gas hydrate occurs mainly within the sediment pore spaces, with concentrations ranging between 50-90%. Laboratory experiments conducted on reconstituted core samples have quantified the effects of pore water salinity and porous media conditions on pressure-temperature stability, suggesting that the partition between gas hydrate stability and instability should be considered as a phase boundary envelope or zone, rather than a discrete threshold. Strength testing on natural core samples has documented the dramatic changes in physical properties following gas hydrate dissociation, with sediments containing no hydrate behaving as unconsolidated sands. While operational problems limited the duration of the production test, a vigorous reservoir response to pressure draw down was observed with increasing gas flow during the testing period. We interpret that pressure temperature (P-T) conditions within the test zone were close to the gas hydrate phase equilibrium threshold, with dissociation initiated at 10 MPa bottomhole pressure (BHP), approximately 1 MPa below in situ conditions. The observation of an increase in production rates at approximately 8.2 MPa BHP may be consistent with the notion of an indistinct gas hydrate stability threshold, with rates increasing as P-T conditions traverse the phase boundary envelope. Significant sand inflow to the well during the test is interpreted to result from the loss of sediment strength during gas hydrate dissociation, with the sediment behaving as a gasified slurry. The increase in gas production rates during the final hours of the test may result from non-uniform gas hydrate dissociation and be affected by accelerated dissociation along water filled natural fractures or fine-scale geologic heterogeneities. These may initiate worm hole or high permeability conduits in association with sand production.

In the numerical simulation of fluid flow and solute transport in porous media, finite element method (FEM) has long been utilized and has been proven to be efficient. In this work, an alternative approach called random walk particle tracking (RWPT) method is proposed. In this method, a finite number of particles represent the distribution of a solute mass. Each particle carries a certain fraction of the total mass and moves in the porous media according to the velocity field. The proposed RWPT model is established on a scientific software platform OpenGeoSys (OGS), which is an open source initiative for numerical simulation of thermo-hydro-mechanical-chemical (THMC) processes in porous media. The flow equation is solved using finite element method in OGS. The obtained hydraulic heads are numerically differentiated to obtain the velocity field. The particle tracking method does not solve the transport equation directly but deals with it in a physically stochastic manner by using the velocity field. Parallel computing concept is included in the model implementation to promote computational efficiency. Several benchmarks are developed for the particle tracking method in OGS to simulate solute transport in porous media and pore space. The simulation results are compared to analytical solutions and other numerical methods to test the presented method. The particle tracking method can accommodate Darcy flow as it is the main consideration in groundwater flow. Furthermore, other flow processes such as Forchheimer flow or Richards flow can be combined with as well. Two applications indicate the capability of the method to handle theoretical real-world problems. This method can be applied as a tool to elicit and discern the detailed structure of evolving contaminant plumes. / Bei der numerischen Simulation von Strömung und Stofftransport in porösen Medien hat die Nutzung der Finite-Elemente-Methode (FEM) eine lange Tradition und wird sich als effizient erweisen. In dieser Arbeit wird ein alternativer Ansatz, die random walk particle tracking (RWPT) Methode vorgeschlagen. Bei diesem Verfahren stellt eine endliche Anzahl von Partikeln die Verteilung eines gelösten Stoffes dar. Jedes Teilchen trägt einen bestimmten Bruchteil der Gesamtmasse und bewegt sich in den porösen Medien gemäß des Geschwindigkeitsfeldes. Das vorgeschlagene RWPT Modell basiert auf der wissenschaftlichen Softwareplattform OpenGeoSys (OGS), die eine Open-Source-Initiative für die numerische Simulation thermo-hydro-mechanisch-chemischen (THMC) in porösen Medien darstellt. Die Strömungsgleichung wird in OGS mit der Finite-Elemente-Methode gelöst. Der Grundwasserstand wird numerisch berechnet, um das Geschwindigkeitsfeld zu erhalten. Die Partikel-Tracking-Methode löst die Transportgleichung nicht direkt, sondern befasst sich mit ihr in einer physikalisch stochastische Weise unter Nutzung des Geschwindigkeitsfeldes. Zur Berücksichtigung der Recheneffizienz ist ein Parallel Computing-Konzept in der Modell-Implementierung enthalten. Zur Simulation des Stofftransports in porösen Medien und im Porenraum wurden mehrere Benchmarks für die Partikel-Tracking-Methode in OGS entwickelt. Die Simulationsergebnisse werden mit analytischen Lösungen und andere numerische Methoden verglichen, um die Aussagefähigkeit des vorgestellten Verfahrens zu bestätigen. Mit der Partikel-Tracking-Methode kann die Darcy-Strömung gelöst werden, die das wichtigste Kriterium in der Grundwasserströmung ist. Außerdem bewältigt die Methode auch andere Strömungsprozesse, wie die Forchheimer-Strömung und die Richards-Strömung. Zwei Anwendungen zeigen die Leistungsfähigkeit der Methode bei der prinzipiellen Handhabung von Problemen der realen Welt. Die Methode kann als ein Instrument zur Aufdeckung Erkennung der detaillierte Struktur von sich entwickelnden Schadstofffahnenangewendet werden.

A procedure is presented to analyze select geometric and effective properties of the porous transport layer (PTL) of the polymer electrolyte membrane fuel cell (PEMFC) in com- pressed and uncompressed states using micro-computed X-ray tomography (Micro CT). A method of compression using a novel device design was employed to mimic the non-homogeneous compression conditions found in functioning fuel cells. The process also features open source image processing and CFD analysis through the use of software packages Fiji and OpenFOAM (proprietary software is also used such as Matlab). Tomographic images of a PTL sample in different compressive states are first analyzed by measuring local porosity values in the through-plane and both in- plane directions. The objective of this study was to develop a method for imaging the PTL structure to show directionality within its properties using relatively inexpensive and non-destructional means. Three different PTL types were tested, one without any additives, one with Polytetrafluoroethylene (PTFE) and one with PTFE and a microporous layer (MPL). Non-homogeneous porosity was shown to exist with the highest and least variable porosity values obtained from the in-plane direction that was in-line with the direction of fibres. Porosity values compared well with values obtained from the literature. The profile of the PTL with MPL added was unattainable using this procedure as the resolution of the Micro CT was too low to resolve its pore space. The next stage involved the effective properties analysis which included effective electronic conductivity and effective diffusivity. It was found that the through-plane values for the effective electronic conductivity study were higher than expected. The ratio between through-plane and in-plane was found to be much higher than expected from literature. Lack of sufficient resolution of fibre contacts has been shown to play a role in this discrepancy. These contact problems were shown not too affect measurements of diffusivity in the pore phase. The in-plane direction parallel to the direction of fibres was found to have the highest values of effective transport properties. Effective diffusivity ratios of between 0.1 and 0.37 were found to be reasonable with the limited experimental evidence found in literature. The it was found that the Bruggeman relation for calculating diffusivity and percolation theory by Tomadakis and Sotirchos over predicted the values for diffusion within the PTL and it is suggested that these theories are not suitable for predicting diffusivity for this material. / Thesis (Master, Mechanical and Materials Engineering) -- Queen's University, 2012-02-02 15:46:29.395

No description available.

This thesis is motivated by problems studied as part of the Diavik Waste Rock Pile Project. Located at the Diavik Diamond Mine in the Northwest Territories, with academic support from the University of Waterloo, the University of Alberta, and the University of British Columbia, this project focuses on constructing mine waste rock piles and studying their physical and chemical properties and the transport processes within them. One of the main reasons for this investigation is to determine the effect of environmental factors on acid mine drainage (AMD) due to sulfide oxidation and the potential environmental impact of AMD. This research is concerned with modeling pressure driven flow through waste rock piles. Unfortunately, because of the irregular shape of the piles, very little data for fluid flow about such an obstacle exists, and the numerical techniques available to work with this domain are limited. Since this restricts the study of the mathematics behind the flow, this thesis focuses on a cylindrical domain, since flow past a solid cylinder has been subjected to many years of study. The cylindrical domain also facilitates the implementation of a pseudo-spectral method. This thesis examines a pressure driven flow through a cylinder of variable permeability subject to turbulent forcing. An equation for the steady flow of an incompressible fluid through a variable permeability porous medium is derived based on Darcy's law, and a pseudo-spectral model is designed to solve the problem. An unsteady time-dependent model for a slightly compressible fluid is then presented, and the unsteady flow through a constant permeability cylinder is examined. The steady results are compared with a finite element model on a trapezoidal domain, which provides a better depiction of a waste rock pile cross section.

In this investigation, the flow friction associated with laminar pulsating flows through porous media was numerically studied. The problem is of interest for understanding the regenerators of Stirling and pulse tube cryocoolers. Two-dimensional flow in a system composed of a number of unit cells of generic porous structures was simulated using a CFD tool, with sinusoidal variations of flow with time. Detailed numerical data representing the oscillating velocity and pressure variations for five different generic porous structure geometries in the porosity range of 0.64 to 0.84, with flow pulsation frequency of 40 Hz were obtained, and special attention was paid to the phase shift characteristics between the velocity and pressure waves. Based on these detailed numerical data, the standard unsteady volume-averaged momentum conservation equation for porous media was then applied in order to obtain the instantaneous as well as cycle-averaged permeability and Forchheimer coefficients. It was found that the cycle-averaged permeability coefficients were nearly the same as those for steady flow, but the cycle-averaged Forchheimer coefficients were about two times larger than those for steady flow. Significant phase lags were observed with respect to the volume-averaged velocity and pressure waves. The parametric trends representing the dependence of these phase lags on porosity and flow Reynolds number were discussed. The phase difference between pressure and velocity waves, which is important for pulse tube cryocooling, depended strongly on porosity and flow Reynolds number.

Experimental and numerical studies simulating the gas-cooled divertor plate design concept have been carried out. While thermo-fluid and thermo-mechanical analyses have been previously performed to show the feasibility of the divertor plate design and its ability to accommodate a maximum heat flux of up to 10 MW/m2, no experimental data have heretofore been published to support or validate such analyses. To that end, this investigation has been undertaken. A test module with prototypical cross-sectional geometry has been designed, constructed, and instrumented. Experiments spanning the prototypical Reynolds numbers of the helium-cooled divertor have been conducted using pressurized air as the coolant. A second test module where the planar jet exiting the inlet manifold is replaced by a two-dimensional hexagonal array of circular jets over the entire top surface of the inlet manifold has also been tested. The thermal performance of both test modules with and without a porous metallic foam layer in the gap between the outer surface of the inlet manifold and the cooled surfaces of the pressure boundary were directly compared. For a given mass flow rate, the slot design with the metallic foam insert showed the highest heat transfer coefficient, with a pressure drop lower than that of the array of circular jets without foam. Additionally, numerical simulations matching the experimental operating conditions for the two cases without foam were performed using the computational fluid dynamics software package, FLUENT® v6.2. Comparisons of the experimental and numerical pressure drop, temperature, and heat transfer coefficient were made.

Το κίνητρο της παρούσας εργασίας είναι η μελέτη των μηχανισμών καταβύθισης μέσω της in-situ ανάμειξης και καταβύθισης δυσδιάλυτων αλάτων( ανθρακικό ασβέστιο) πάνω στις επιφάνειες αλλά και κατά μήκος των πορώδων υλικών. Οι επικαθίσεις δυσδιάλυτων αλάτων σε πορώδη σχηματισμούς, αποτελεί ένα σημαντικό πρόβλημα σε βιομηχανικές εφαρμογές. Για παράδειγμα, στην παγκόσμια βιομηχανία πετρελαίου και φυσικού αερίου πολλές περιοχές άντλησης πετρελαίου αντιμετωπίζουν τεράστιο πρόβλημα επικαθίσεων αλάτων, με αποτέλεσμα την δραματική μείωση της παραγωγής πετρελαίου. Επιπλέον, προβλήματα εναπόθεσης αλάτων μειώνουν την εκμετάλλευση γεωθερμικής ενέργειας, τη διάρκεια ζωής μεμβρανών αφαλάτωσης, καθώς και την απομόνωση του CO2 σε υπόγεια πηγάδια, κτλ. Οι επικαθίσεις αλάτων σε πηγάδια άντλησης πετρελαίου προέρχονται είτε από απευθείας καταβύθιση των διαλυμένων αλάτων που βρίσκονται στο νερό που υπάρχει σε υπόγειες κοιλότητες και καταβυθίζονται εξαιτίας της διαφοροποίησης των τοπικών συνθηκών(κυρίως πίεσης και θερμοκρασίας), είτε από τις μεγάλες ποσότητες νερού που χρησιμοποιείται για την εκτόπιση του πετρελαίου στα πηγάδια παραγωγής. Ένα από τα σημαντικότερα προβλήματα που αντιμετωπίζουν οι εταιρίες εξόρυξης πετρελαίου είναι οι μειωμένοι ρυθμοί άντλησης λόγω της μείωσης του τοπικού πορώδους και της διαπερατότητας από τη δημιουργία και την ανάπτυξης επικαθίσεων στα πηγάδια εξόρυξης. Οι επικαθίσεις αλάτων μπορούν επιπλέον να φράξουν τις σωληνώσεις ή να σχηματίσουν ένα λεπτό στρώμα στα τοιχώματα των σωληνώσεων που χρησιμοποιούνται στην παραγωγή. Ο σχηματισμός επικαθίσεων έχει σαν αποτέλεσμα την δημιουργία λειτουργικών προβλημάτων και δυσκολιών, που οδηγούν σε επιπλέον κόστος και σε εξαιρετικές περιπτώσεις στην εγκατάλειψη του πηγαδιού άντλησης πετρελαίου. Αρχικά έγινε μελέτη της αυθόρμητης καταβύθισης του ανθρακικού ασβεστίου σε αντιδραστήρες διαλείποντος έργου. Μελετήθηκε επίσης και η επίδραση της ίδιας ποσότητας άμμου σε διαφορετικούς αρχικούς υπερκορεσμούς κατά την αυθόρμητη καταβύθιση του ανθρακικού ασβεστίου. Καθ’ όλη τη διάρκεια των πειραμάτων συλλέγονταν δείγματα από τον αντιδραστήρα και ανιχνευόταν η μείωση της συγκέντρωσης των ιόντων ασβεστίων με τη μέθοδο της ατομικής απορρόφησης. Από τη μείωση της συγκέντρωσης των ιόντων ασβεστίου υπολογίστηκε ο ρυθμός καταβύθισης και στη συνέχεια η φαινόμενη τάξη της αντίδρασης και η σταθερά ταχύτητας. Σε όλα τα πειράματα ανιχνεύτηκε με περίθλαση ακτίνων Χ και ηλεκτρονική μικροσκοπία σάρωσης, ο ασβεστίτης ως η μοναδική καταβυθιζόμενη φάση. Κατά την αυθόρμητη καταβύθιση του ανθρακικού ασβεστίου απουσία άμμου, ο κύριος μηχανισμός καταβύθισης βρέθηκε ότι είναι η επιφανειακή διάχυση των δομικών μονάδων σε κατάλληλες θέσεις στην επιφάνεια των κρυστάλλων, ενώ στα πειράματα παρουσία ποσότητας άμμου οι κρύσταλλοι αναπτύσσονται με βάση το πολυπυρηνικό πρότυπο. Στην παρούσα εργασία μελετάται επίσης η διεργασία καταβύθισης του δυσδιάλυτου άλατος του ανθρακικού ασβεστίου (CaCO3) σε δισδιάστατα πορώδη δοκίμια Plexiglas®. Μέσω της οπτικής παρατήρησης της ανάπτυξης των επικαθίσεων, επιχειρήθηκε η εξαγωγή ποιοτικών και ποσοτικών πληροφοριών που αφορούν τους μηχανισμούς και τον τρόπο με τον οποίο αναπτύσσονται οι κρύσταλλοι ανθρακικού ασβεστίου στα τοιχώματα του προσομοιωτή πορώδους μέσου κάτω από ρεαλιστικές συνθήκες ροής, θερμοκρασίας και συγκέντρωσης. Δύο ευδιάλυτα διαλύματα χλωριούχου ασβεστίου, CaCl2, και όξινου ανθρακικού ασβεστίου, NaHCO3, αναμιγνύονταν λίγο πριν την είσοδό τους στο πορώδες μέσο με τη βοήθεια εμβολοφόρων αντλιών. Η καταβύθιση του ανθρακικού ασβεστίου λάμβανε μέρος μέσα στο πορώδες γραμμικό μέσο. Χρησιμοποιώντας ένα οπτικό μικροσκόπιο σε συνδυασμό με βιντεοκάμερα ήταν δυνατή η παρακολούθηση και η καταγραφή των κρυστάλλων που εμφανίζονταν καθώς και η ανάπτυξη τους κατά μήκος του καναλιού, στις θέσεις 2 cm, 6 cm, 8 cm και 10 cm από την είσοδο του. Σε κάθε θέση μετρήθηκε η μεταβολή του pH και της συγκέντρωσης του ασβεστίου με το χρόνο και επιπλέον καταγράφηκε ο σχηματισμός των κρυστάλλων του ανθρακικού ασβεστίου με τη βοήθεια της κάμερας του μικροσκοπίου. Για τη μέτρηση της συγκέντρωσης των ιόντων ασβεστίων χρησιμοποιήθηκε η φασματοφωτομετρία ατομικής απορρόφησης ενώ ο υπολογισμός του μεγέθους των κρυστάλλων έγινε με τη βοήθεια του προγράμματος Adobe Photoshop CS5.1. Κατά την διάρκεια του πειράματος ο υπερκορεσμός μειωνόταν κατά μήκος του καναλιού, ως αποτέλεσμα της in situ ανάμειξης των διαλυμάτων του χλωριούχου ασβεστίου και ανθρακικού νατρίου. Τέλος, κοντά στην είσοδο του καναλιού ο αριθμός και το μέγεθος των κρυστάλλων που καταβυθιστήκαν ήταν μεγαλύτερο, συγκριτικά με αυτούς που καταβυθιστήκαν κοντά στην έξοδο. / The deposition of sparingly soluble salts in porous formations is a severe problem encountered in industrial applications. For example, in international oil and gas industry many oil fields in operation face the major problem of scale which reduces dramatically oil production. Moreover, scaling problems limit the utilization of geothermal energy, the operational life of desalination membranes, in CO2 sequestration in subsoil wells, etc. Scale deposits encountered in oil fields are formed either by direct precipitation of the dissolved salts present in the formation water and precipitate because of the change of the local conditions (mainly pressure and temperature) or because of the large quantities of the injection water used for the displacement of the oil in the production wells. In the course of scale formation and development at wellbores, both local porosity and permeability may be reduced resulting in lower oil recovery rates. Scale deposits may also block flow by clogging perforations or by forming a thick layer on the wall of the production tubing. The formation of scales may result in a series of operational problems and difficulties, and in turn, result to additional capital cost, substantial operating cost and in extreme cases in the abandonment of the oil wells. The injection water (sea water) used in oil production is a stream most often incompatible with the well formation water. Mixing the incompatible streams interact chemically resulting in the precipitation of minerals. The high concentrations of SO42- in sea water and the low concentrations of Ba2+/Sr2+ in the formation waters constitute solutions where the concentrations of SO42- are low while the concentrations of Ca2+, Ba2+ and Sr2+ are high. Mixing of these streams, there is a high probability of CaSO4, BaSO4 and/or SrSO4 precipitation. Some of the more common ions frequently encountered in oil field waters and cause precipitation in incompatible waters are Ca2+, Ba2+, Sr2+, Fe2+, HCO3- and SO42-. In most cases, the low solubility and tenaciously adhering calcium carbonate and barium sulfate scale deposits are the most commonly salts found in well scale deposits. Knowledge of the mechanisms underlying the formation of scale deposits and the concomitant effects on the changes of the local porosity and permeability is of paramount importance for the efficient design of scale control processes. On the other hand there are beneficial applications involving the precipitation of insoluble salts. The controlled precipitation of sparingly soluble salts may be used for the effective consolidation of unconsolidated or poorly consolidated formations, like the sandy soils or granular materials (applications: protection of cultivated soils from erosion by raining water, protection of sandy reservoirs from the intrusion of sand grains within the oil wells, waterproofing of underground constructions, e.g. tunnels or other concrete constructions, etc). The main motivation for the present work is the investigation of the precipitation mechanisms through in-situ mixing and the precipitation of sparingly soluble salts (calcium carbonate) in porous materials. The proposed method may be applicable for the prevention of sand entrainment into well bores during oil production in poorly consolidated oil reservoirs. In the present study two series of experiments were done in order to investigate the kinetics of calcium carbonate precipitation in batch reactors using calcium chloride and sodium bicarbonate as the initial reagents in the absence and presence of sand. In all supersaturated solutions used in the present work, the stoichiometric of total calcium ions to total bicarbonate ions was equal to 1:1. During the experiments samples were taken, filtered through membrane filters and analyzed for total calcium by atomic absorption spectroscopy. The analysis of the kinetic data showed that in spontaneous precipitation experiments the crystals grow with spiral mechanism (n=2) and in the presence of sand kinetics data were best fitted by the polynuclear model (n>2). At the end of each experiment, the fluid suspension was filtered with a 0.22 μm filter and the solids were collected, dried at room temperature, and characterized by powder X-ray diffraction (XRD) and scanning electron microscopy (SEM). Finally, also a series of experiments were done using two dimensional porous (Plexiglas®) medium models in which insoluble salts were nucleated from supersaturated solutions. Through the direct observation of the evolution of the deposits it was attempted to obtain qualitative and quantitative information concerning the mechanisms and pattern of deposition and growth of calcium carbonate scale on pore walls under realistic geometrical, flow, temperature and concentration conditions. Two soluble salts of calcium chloride, CaCl2, and sodium bicarbonate, NaHCO3, were mixed just before the injection point using syringe pumps. Precipitation of calcium carbonate took place within the porous medium. Using an optical microscope in combination with a video- camera it was possible to monitor crystal appearance and growth at a steady state. The growth of the calcium carbonate crystals in different positions within the etched linear porous medium as a function of time was observed. The visual observation of the crystal growth took place in different positions along the channel (2 cm, 6 cm, 8 cm and 10 cm from the inlet of the channel). In every position the changes of the pH and calcium concentration were measured and the formation of crystals was recorded through image capture. During the flow of the mixture through the channel, supersaturation decreased as a result of the in situ mixing of calcium and bicarbonate ions. As it was expected precipitation was favoured near the inlet of the cell and few crystals were precipitated near the outlet. Furthermore, the size of the crystal was bigger near the inlet of the channel in comparison to the crystals that have precipitate near the outlet of the channel.

Since 1996, super-critical CO$_2$ has been injected at a rate of $\sim$0.85~Mt~yr$^{-1}$ into a pristine, saline aquifer at the Sleipner carbon capture and storage project. A suite of time-lapse, three-dimensional seismic reflection surveys have been acquired over the injection site. This suite includes a pre-injection survey acquired in 1994 and seven post-injection surveys acquired between 1999 and 2010. Nine consistently bright reflections within the reservoir, mapped on all post-injection surveys, are interpreted to be thin layers of CO$_2$ trapped beneath mudstone horizons. The areal extents of these CO$_2$ layers are observed to either increase or remain constant with time. However, volume flux of CO$_2$ into these layers has proven difficult to measure accurately. In addition, the complex planform of the shallowest layer, Layer 9, has proven challenging to explain using reservoir simulations. In this dissertation, the spatial distribution of CO$_2$ in Layer~9 is measured in three dimensions using a combination of seismic reflection amplitudes and changes in two-way travel time between time-lapse seismic reflection surveys. The CO$_2$ volume in this layer is shown to be growing at an increasing rate through time. To investigate CO$_2$ flow within Layer~9, a numerical gravity current model that accounts for topographic gradients is developed. This vertically-integrated model is computationally efficient, allowing it to be inverted to find reservoir properties that minimise differences between measured and modelled CO$_2$ distributions. The best-fitting reservoir permeability agrees with measured values from nearby wells. Rapid northward migration of CO$_2$ in Layer~9 is explained by a high permeability channel, inferred from spectral decomposition of the seismic reflection surveys. This numerical model is found to be capable of forecasting CO$_2$ flow by comparing models calibrated on early seismic reflection surveys to observed CO$_2$ distributions from later surveys. Numerical and analytical models are then used to assess the effect of the proximity of an impermeable base on the flow of a buoyant fluid, motivated by the variable thickness of the uppermost reservoir. Spatial gradients in the confinement of the reservoir are found to direct the flow of CO$_2$ when the current is of comparable thickness to the reservoir. Finally, CO$_2$ volume in the second shallowest layer, Layer~8, is measured using structural analysis and numerical modelling. CO$_2$ in Layer~8 is estimated to have reached the spill point of its structural trap by 2010. CO$_2$ flux into the upper two layers is now $\sim$40\% of total CO$_2$ flux injected at the base of the reservoir, and is increasing with time. This estimate is supported by observations of decreasing areal growth rate of the lower layers. The uppermost layers are therefore expected to contribute significantly to the total reservoir storage capacity in the future. CO$_2$ flow within Layer~9 beyond 2010 is forecast to be predominantly directed towards a topographic dome located $\sim$3~km north of the injection point. This dissertation shows that advances in determining the spatial distribution and flow of CO$_2$ in the sub-surface can be made by a combination of careful seismic interpretation and numerical flow modelling.