Modeling of flow in porous and fractured media is a very important problem in reservoir engineering. As for numerical simulations conventional Navier-Stokes codes are applied to flow in both porous and fractured media. But they have long computation times, poor convergence and problems of numerical instabilities. Therefore, it is desired to develop another computational method that is more efficient and use simple rules to represent the flow in fractured media rather than partial differential equations. In this thesis Lattice Boltzmann Automaton Model will be used to represent the single phase fluid flow in two dimensional synthetic fractures and the simulation results obtained from this model are used to train Artificial Neural Networks. It has been found that as the mean aperture-fractal dimension ratio increases permeability increases. Moreover as the anisotropy factor increases permeability decreases with a second order polynomial relationship.

Abstract The Late Devonian-Early Carboniferous Drummond Basin of central Queensland, Australia is one of the oldest units of the New England Fold Belt. Samples of authigenic clay minerals have been collected from across the Drummond Basin and analysed using mineralogic, isotopic and geochronologic techniques. X-Ray diffraction analysis on 36 samples selected from six exploration wells was carried out on volcaniclastic and argillaceous sandstones and mudstones of the Drummond Basin sequence in order to ascertain authigenic illite content. Sixteen samples were found to be suitable for K-Ar dating purposes, containing abundant illite and illite/smectite mixed layers (I/S). K-Ar age dating was carried out on these selected samples yielding three age groups at 300, 250 and 200 Ma. Temperature of formation was estimated for the authigenic illite and illite/smectite minerals, based on illite crystallinity and illite content of illite/smectite mixed layers. Two groups of temperatures have been delineated, with a high temperature assemblage of between 275¨¬C-320¨¬C and a low temperature assemblage of approximately 100¨¬C. When compared with the age groups, it is seen that the older event was typically hotter. Stable isotope compositions for the authigenic clay minerals were analysed with ¥äD and ¥ä18O values of between -68¢¶ to -117¢¶ and +7.5¢¶ to +14.4¢¶ (SMOW) respectively. Isotopic compositions for the fluid in equilibrium with the clay minerals were calculated using palaeotemperatures estimated from the clay mineralogy and range between -90¢¶ and -41¢¶ and -0.6¢¶ and +8.7¢¶ (SMOW) for hydrogen and oxygen respectively. These compositions are more depleted in deuterium and more enriched in 18O than those previously reported for basins in the region. The data also produced three groups of fluid compositions, which correlate with the age distributions delineated by the K-Ar geochronology. The first fluid composition, which correlated with the Permo-Carboniferous age (~300 Ma) is interpreted to result from a mixing of Permo-Carboniferous meteoric water with highly evolved formation water trapped within the Drummond Basin since deposition of the basin. The regional extensional event at the Permo-Carboniferous boundary, which initiated the opening of the Bowen Basin to the east of the Drummond Basin allowed meteoric water to percolate deep into the basin, enabling the growth of illite rich clay minerals at this time. The second fluid group is considered to represent the evolution of Late Triassic meteoric water migrating through the Drummond Basin during another extensional period in the region at approximately 200 Ma. This event precipitated a new period of clay mineral growth in the basin (reflected in the K-Ar dates), raised the geothermal gradient in the region and extended the crust, forming new sedimentary basins within the New England Fold Belt. The third assemblage of fluid compositions were found to be very similar to the Late Triassic group, with ¥äD and ¥ä18O values plotting very close to the Late Triassic meteoric water evolution trend. This is interpreted as representing the same fluid as that which produced the Late Triassic compositions and therefore the 250 Ma ages may represent partial resetting of older ages by the Late Triassic thermal episode. However, if it is assumed that these K-Ar ages are not old ages partially reset by a thermal event, the stable isotope data could represent a mixing of Mid-Late Permian meteoric water with formation water trapped in the basin. This event may represent the regional compressive Hunter-Bowen event and may be responsible for dewatering the basin before the Late Triassic extensional episode.

Geodynamic modelling, via computer simulations, offers an easily controllable method for investigating the behaviour of an Earth system and providing feedback to conceptual models of geological evolution. However, most available computer codes have been developed for engineering or hydrological applications, where strains are small and post-failure deformation is not studied. Such codes cannot simultaneously model large deformation and porous fluid flow. To remedy this situation in the face of tectonic modelling, a numerical approach was developed to incorporate porous fluid flow into an existing high-deformation code called Ellipsis. The resulting software, with these twin capabilities, simulates the evolution of highly deformed tectonic regimes where fluid flow is important, such as in mineral provinces. A realistic description of deformation depends on the accurate characterisation of material properties and the laws governing material behaviour. Aside from the development of appropriate physics, it can be a difficult task to find a set of model parameters, including material properties and initial geometries, that can reproduce some conceptual target. In this context, an interactive system for the rapid exploration of model parameter space, and for the evaluation of all model results, replaces the traditional but time-consuming approach of finding a result via trial and error. The visualisation of all solutions in such a search of parameter space, through simple graphical tools, adds a new degree of understanding to the effects of variations in the parameters, the importance of each parameter in controlling a solution, and the degree of coverage of the parameter space. Two final applications of the software code and interactive parameter search illustrate the power of numerical modelling within the feedback loop to field observations. In the first example, vertical rheological contrasts between the upper and lower crust, most easily related to thermal profiles and mineralogy, exert a greater control over the mode of crustal extension than any other parameters. A weak lower crust promotes large fault spacing with high displacements, often overriding initial close fault spacing, to lead eventually to metamorphic core complex formation. In the second case, specifically tied to the history of compressional orogenies in northern Nevada, exploration of model parameters shows that the natural reactivation of early normal faults in the Proterozoic basement, regardless of basement topography or rheological contrasts, would explain the subsequent elevation and gravitationally-induced thrusting of sedimentary layers over the Carlin gold trend, providing pathways and ponding sites for mineral-bearing fluids.

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.

Made available in DSpace on 2014-10-09T12:41:29Z (GMT). No. of bitstreams: 0 / Made available in DSpace on 2014-10-09T14:03:20Z (GMT). No. of bitstreams: 0 / Tese (Doutoramento) / IPEN/T / Instituto de Pesquisas Energeticas e Nucleares - IPEN-CNEN/SP

Made available in DSpace on 2014-10-09T12:38:39Z (GMT). No. of bitstreams: 0 / Made available in DSpace on 2014-10-09T14:05:15Z (GMT). No. of bitstreams: 1 06042.pdf: 9369626 bytes, checksum: 37950a2f878d6535f671de4a025da71c (MD5) / Tese (Doutoramento) / IPEN/T / Instituto de Pesquisas Energeticas e Nucleares - IPEN/CNEN-SP

Made available in DSpace on 2014-10-09T12:55:46Z (GMT). No. of bitstreams: 0 / Made available in DSpace on 2014-10-09T14:05:12Z (GMT). No. of bitstreams: 0 / Tese (Doutoramento) / IPEN/T / Instituto de Pesquisas Energeticas e Nucleares - IPEN-CNEN/SP

Made available in DSpace on 2014-10-09T12:54:00Z (GMT). No. of bitstreams: 0 / Made available in DSpace on 2014-10-09T14:09:26Z (GMT). No. of bitstreams: 1 12431.pdf: 23287148 bytes, checksum: ca5f3b291fb81099ce731e9a131ce645 (MD5) / Tese (Doutoramento) / IPEN/T / Escola Politecnica, Universidade de Sao Paulo - POLI/USP

Made available in DSpace on 2014-10-09T12:41:29Z (GMT). No. of bitstreams: 0 / Made available in DSpace on 2014-10-09T14:03:20Z (GMT). No. of bitstreams: 0 / Circuitos de convecção natural ou sistemas de circulação natural são empregados em diversas áreas da engenharia. Reatores nucleares refrigerados a água utilizam circuitos de circulação natural como método passivo de seguranca. Em situações críticas, sem qualquer controle externo, o sistema permanece em segurança por suas próprias características de funcionamento (intrinsecamente seguro). O trabalho proposto consiste em estudar numericamente o circuito de circulação natural de água, localizado no Instituto de Pesquisas Energéticas e Nucleares / Comissão Nacional de Energia Nuclear em São Paulo, por meio do uso de modelos matemáticos, objetivando determinar o padrão do escoamento em condições sem mudança de fase líquido-vapor. A comparação dos resultados de temperatura obtidos por cada um dos modelos de turbulência aos pontos instrumentados no circuito experimental, na condição transitória, revelou desvios significativos nas respostas do modelo de zero equação. Desvios intermediário foram observados nos modelos de transporte da viscosidade turbulenta (EVTE), k - ω, SST e SSG e resultados melhores foram vericados nos modelos k - ε e DES (com significativa superioridade do primeiro modelo). / Tese (Doutoramento) / IPEN/T / Instituto de Pesquisas Energeticas e Nucleares - IPEN-CNEN/SP

Made available in DSpace on 2014-10-09T12:38:39Z (GMT). No. of bitstreams: 0 / Made available in DSpace on 2014-10-09T14:05:15Z (GMT). No. of bitstreams: 1 06042.pdf: 9369626 bytes, checksum: 37950a2f878d6535f671de4a025da71c (MD5) / Tese (Doutoramento) / IPEN/T / Instituto de Pesquisas Energeticas e Nucleares - IPEN/CNEN-SP