In this thesis a number of issues related to oceanographic problems have been dealt with on the basis of applying asymptotic methods.  The first study focused on the tidal generation of internal waves, a process which is quantifed by the conversion rates. These have traditionally been calculated by using the WKB approximation. However, the systematic imprecision of this theory for the lowest modes as well as turbulence at the seabed level affect the results. To handle these anomalies we introduced another asymptotic technique, homogenization theory, which led to signifcant improvements, especially for the lowest modes.  The second study dealt with the dynamical aspects of a nonlinear oscillator which can be interpreted as a variant of the classical two-box models used in oceanography. The system is constituted by two connected vessels containing a fluid characterised by a nonlinear equation of state and a large volume differences between the vessels is prescribed. It is recognised that the system, when performing relaxation oscillations, exhibits almost-discontinuous jumps between the two branches of the slow manifold of the problem. The lowest-order analysis yielded reasonable correspondence with the numerical results.  The third study is an extension of the lowest-order approximation of the relaxation oscillations undertaken in the previous paper. A Mandelstam condition is imposed on the system by assuming that the total heat content of the system is conserved during the discontinuous jumps.  In the fourth study an asymptotic analysis is carried out to examine the oscillatory behaviour of the thermal oscillator. It is found that the analytically determined corrections to the zeroth-order analysis yield overall satisfying results even for comparatively large values of the vessel-volume ratio. / At the time of the doctoral defense, the following papers were unpublished and had a status as follows: Paper 3: Manuscript. Paper 4: Manuscript.

The objective of this research effort is to provide a methodology for assessing the effectiveness of sampling techniques used to gather different types of geo-physical information by a robotic agent. We focus on assessing how well unique real-time sampling strategies acquire information that is, otherwise, too dangerous or costly to collect by human scientists. Traditional sampling strategies and informed search tech- niques provide the underlying structure for a navigating robotic surveyor whose goal is to collect samples that yield an accurate representation of the measured phenomena under realistic constraints. These sampling strategies are alternative improvements that provide greater information gain than current sampling technology allows. The contributions of this work include the following: 1) A method for estimating spa- tially distributed phenomena, using a partial sample set of information, that shows improvement over that of a more traditional estimation method. 2) A method for sampling this phenomena in the form of a navigation scheme for a mobile robotic survey system. 3) A method of ranking and comparing different navigation algorithms relative to one another based on performance (reconstruction error) and resource (distance) constraints. We introduce a specific class of navigation algorithms as example sampling strategies to demonstrate how our methodology allows different robot navigation options to be contrasted and the most practical strategy selected.

The dynamics of the world's oceans occur at a vast range of length scales. Although there are theories that aid in understanding the dynamics at planetary scales and microscales, the motions in between are still not yet well understood. This work discusses a numerical model to study barotropic wind-driven gyre flow that is capable of resolving dynamics at the synoptic, O(1000 km), mesoscale, O(100 km) and submesoscales O(10 km). The Quasi-Geostrophic (QG) model has been used predominantly to study ocean circulations but it is limited as it can only describe motions at synoptic scales and mesoscales. The Rotating Shallow Water (SW) model that can describe dynamics at a wider range of horizontal length scales and can better describe motions at the submesoscales. Numerical methods that are capable of high-resolution simulations are discussed for both QG and SW models and the numerical results are compared. To achieve high accuracy and resolve an optimal range of length scales, spectral methods are applied to solve the governing equations and a third-order Adams-Bashforth method is used for the temporal discretization. Several simulations of both models are computed by varying the strength of dissipation. The simulations either tend to a laminar steady state, or a turbulent flow with dynamics occurring at a wide range of length and time scales. The laminar results show similar behaviours in both models, thus QG and SW tend to agree when describing slow, large-scale flows. The turbulent simulations begin to differ as QG breaks down when faster and smaller scale motions occur. Essential differences in the underlying assumptions between the QG and SW models are highlighted using the results from the numerical simulations.

The extent and dynamics of permafrost are tightly linked to the distribution and movement of water in arctic landscapes. As the Arctic warms more rapidly than the global average, profound changes are expected in both permafrost and hydrology; however, much is still not known about the interactions between these two systems. The aim of this thesis is to provide new knowledge on the links between permafrost and hydrology under varying environmental conditions and across different scales. The objectives are to (i) determine how permafrost distributions and patterns in morphology are linked to hydrology, (ii) determine how groundwater flow influences ground temperature dynamics in permafrost landscapes, and (iii) explore the mechanisms that link permafrost to groundwater and streamflow dynamics. A range of methods have been applied within the four studies (papers I-IV) comprising the thesis: geophysical (ground penetrating radar and electrical resistivity tomography) and GIS techniques for mapping and analyzing permafrost distributions and related morphology; numerical modeling of coupled heat and water fluxes for mechanistic understanding permafrost-hydrological links; and statistical analyses for detecting trends in streamflow associated with permafrost thaw. Combining these various methods here allows for, and may be considered a prerequisite for, novel insights to processes. The thesis also presents statistical analyses of field observations of ground temperatures, ground- and surface water levels, as well as lake and shore morphological variables. Discontinuous permafrost peatlands are heterogeneous environments regarding permafrost distributions and thickness which is manifested in surface systems such as lake geometries. In these environments, lateral groundwater fluxes, which are not considered in most permafrost models, can significantly influence ground temperature dynamics, especially during high groundwater gradient conditions. River discharge data provide a potential for monitoring catchment-scale changes in permafrost, as the magnitude and seasonality of groundwater fluxes feeding into streams are affected by the distribution of permafrost. This thesis highlights the need to understand water and permafrost as an integrated system with potential internal feedback processes. For example, permafrost thaw can lead to increases in groundwater discharge which in turn can lead to increased heat transfer through the ground, resulting in further acceleration of permafrost thaw rates. / <p>At the time of the doctoral defense, the following paper was unpublished and had a status as follows: Paper 3: Manuscript.</p>

While geophysical reservoir characterization has been an area of research for the last three decades, geophysical reservoir monitoring, time-lapse studies, have recently become an important geophysical application. Generally speaking, the main target is to detect, estimate, and discriminate the changes in subsurface rock properties due to production. This research develops various sensitivity and feasibility analyses to investigate the effects of production-induced time-lapse changes on geophysical measurements including seismic and controlled-source electromagnetic (CSEM) data. For doing so, a realistic reservoir model is numerically simulated based on a prograding near-shore sandstone reservoir. To account for the spatial distribution of petrophysical properties, an effective porosity model is first simulated by Gaussian geostatistics. Dispersed clay and dual water models are then efficiently combined with other well-known theoretical and experimental petrophysical correlations to consistently simulate reservoir model parameters. Next, the constructed reservoir model is subjected to numerical simulation of multi-phase fluid flow to replicate a waterflooding scenario of a black oil reservoir and to predict the spatial distributions of fluid pressure and saturation. A modified Archie’s equation for shaly sandstones is utilized to simulate rock resistivity. Finally, a geologically consistent stress-sensitive rock physics model, combined with the modified Gassmann theory for shaly sandstones, is utilized to simulate seismic elastic parameters. As a result, the comprehensive petro-electro-elastic model developed in this dissertation can be efficiently utilized in sensitivity and feasibility analyses of seismic/CSEM data with respect to petrophysical properties and, ultimately, applied to reservoir characterization and monitoring research. Using the resistivity models, a base and two monitor time-lapse CSEM surveys are simulated via accurate numerical algorithms. 2.5D CSEM modeling demonstrates that a detectable time-lapse signal after 5 years and a strong time-lapse signal after 10 years of waterflooding are attainable with the careful application of currently available CSEM technology. To simulate seismic waves, I employ different seismic modeling algorithms, one-dimensional (1D) acoustic and elastic ray tracing, 1D full elastic reflectivity, 2D split-step Fourier plane-wave (SFPW), and 2D stagger grid explicit finite difference (FD). My analyses demonstrate that acoustic modeling of an elastic medium is a good approximation up to ray parameter (p) equal to 0.2 sec/km. However, at p=0.3 sec/km, differences between elastic and acoustic wave propagation is the more dominant effect compared to internal multiples. Here, converted waves are also generated with significant amplitudes compared to primaries and internal multiples. I also show that time-lapse modeling of the reservoir using SFPW approach is very fast compared to FD, 100 times faster for my case here. It is capable of handling higher frequencies than FD. It provides an accurate image of the waterflooding process comparable to FD. Consequently, it is a powerful alternative for time-lapse seismic modeling. I conclude that both seismic and CSEM data have adequate but different sensitivities to changes in reservoir properties and therefore have the potential to quantitatively map production-induced time-lapse changes. / text

The integrated use of Electromagnetic Ground Conductivity, Ground Electrical Resistivity and Seismic Refraction geophysics methods was made in Alia valley, Eritrea for the aim of developing a methodology best suited for ground water exploration in hard rock areas with limited professional and financial resources. A total area of 25 square kilometer was covered with the geophysical surveys. The effective use of information from wells and geological information from outcrops and existing literature helped to create a preliminary conceptual hydrogeologic model that in turn provided a general picture of the research site in terms of its groundwater potential and related factors. The optimum configuration of the geophysical methods and their interpretation was researched by the use of synthetic mathematical models and experimental surveys in areas of known hydrogeological characteristics like lithology and depth of water table. Two approximately parallel electromagnetic ground conductivity profile lines crossed the survey area in approximately north-south direction to assess the variation in the electrical conductivity of the major geological structures. Vertical electrical sounding (VES) surveys were conducted close to existing wells with lithological information and at some selected sites. Survey sites for the two dimensional electrical resistivity imaging and the seismic refraction surveys were selected based on the results from the electromagnetic ground conductivity surveys. The formal geophysical surveys results were correlated with nearby boreholes and with each other. The responses of the different geophysical methods for each lithological layer were recognized. In addition, the location and characteristics of the major geological structures of the area were studies based on the information from the geophysical surveys and other available information. At the end of the research, parts of the research area with better groundwater potential were identified. As a result of the research, an efficient and cost effective geophysical approach was developed for groundwater exploration in hard rock areas of Eritrea and similar places of the world. Recommendations were made on the wider and effective use of the approach with due consideration for the limited financial, human resource and material capacities available for ground water exploration in most countries of Africa and other countries of the world.

No description available.

Roof falls and accumulation of dangerous gasses are the most common hazards in any underground coal mine. Different mechanisms can jeopardize the stability of the roof in underground excavations and successful roof control can only be obtained if the failure mechanism is identified and understood properly. The presence of discontinuities, the inherent variability of the rock mass and discontinuity properties, and the uncertainties associated with directions and magnitudes of the in-situ stress makes the rock engineering problems challenging. The numerical modeling can assist the ground control engineers in designing and evaluating the stability of the underground excavations. If extensive geological and geotechnical data are available, then detailed predictions of deformation, stress and stability can be accomplished by performing numerical modeling. If not, still the numerical modeling can be used to perform parametric studies to gain insight into the possible ranges of responses of a system due to likely ranges of various parameters. The parametric studies can help to identify the key parameters and their impact on stability of underground excavations. The priorities of the material testing and site investigation can be set based on the selected key parameters from parametric studies. An underground coal mine in western Pennsylvania is selected as a case study mine to investigate the underlying causes of roof falls at this mine. The immediate roof at the case study mine consists of laminated silty shale, shale, or sandstone that changes from area to area, and the floor is shale or soft fireclay. This study was mainly focused in the stability analysis of the roofs with the laminated silty shale rock type, where the majority of roof falls had taken place in the roof with this type of roof material. Extensive laboratory tests were performed on the core samples obtained from the case study mine to estimate the intact rock and discontinuity properties of the materials that occur in large extent at the selected interest area of the case study mine. In this research, the three-dimensional distinct element method was used to investigate the stability of the roof in an underground room-and-pillar coal mine. The implemented technique was able to accurately capture the failure of the major discontinuities and rock masses which consist of intact rock and minor discontinuities. In order to accurately replicate the post failure behavior of the rock layers in the immediate roof area, the strain-softening material constitutive law was applied to this region. Extensive numerical parametric studies were conducted to investigate the effect of different parameters such as the variation of immediate roof rock mass strength properties, variation of discontinuity mechanical properties, orientations and magnitudes of the horizontal in-situ stresses, and the size of pillars and excavations on stability of the excavations. The distribution of post failure cohesion along with other measures such as accumulated plastic shear strain, distribution of Z-displacements at the roofline, failure state (joint slip and tensile failure) and displacement (normal and shear displacements) of discontinuities were used to accurately assess the roof stability in this case study. The research conducted in this dissertation showed that the bedding planes play an important role on the behavior of roof in underground excavations. Therefore, an appropriate numerical modeling technique which incorporates the effect of discontinuities should be employed to simulate the realistic behavior of the discontinuous rock masses such as the layered materials in roof strata of the underground coal mines. The three-dimensional distinct element method used in this research showed the clear superiority of this technique over the continuum based methods.

Mine to Mill optimization is considered as a key concept for metal mining recently. Targeting operational best practices on a highly varying environment is challenging. Impact of underperformed basic operations such as drilling and blasting will sustain this inefficiency in mineral processing. Data provided for each of these operations from software and hardware utilized on field reached a level where advanced data analytics becomes applicable. In order to represent the operations as close to reality, an integrated layer of data where transactional and process based data lives is crucial. Data warehousing and data mining are alternative tools that rely on a robust data structure. Data mining utilizes the integrated data layer for pattern discovery within the data itself. Relations that are unknown for now can be investigated by data mining algorithms that rely on vast amount of data. Empirical equations that are based on a limited set of data could be improved by using data mining algorithms. The main objective of optimizing the mine to mill value chain also challenges the concept of providing real-time feedback. This research proposes a data-driven mine-to-mill framework for modern mines.

No description available.