Integration of high-resolution geophysical data for the characterization of alluvial aquifers : new approaches and implications for hydrological predictions

Abstract Accurate characterization of the spatial distribution of hydrological properties in heterogeneous aquifers at a range of scales is a key prerequisite for reliable modeling of subsurface contaminant transport, and is essential for designing effective and cost-efficient groundwater management and remediation strategies. To this end, high-resolution geophysical methods have shown significant potential to bridge a critical gap in subsurface resolution and coverage between traditional hydrological measurement techniques such as borehole log/core analyses and tracer or pumping tests. An important and still largely unresolved issue, however, is how to best quantitatively integrate geophysical data into a characterization study in order to estimate the spatial distribution of one or more pertinent hydrological parameters, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first develop a strategy for the assimilation of several types of hydrogeophysical data having varying degrees of resolution, subsurface coverage, and sensitivity to the hydrologic parameter of interest. In this regard a novel simulated annealing (SA)-based conditional simulation approach was developed and then tested in its ability to generate realizations of porosity given crosshole ground-penetrating radar (GPR) and neutron porosity log data. This was done successfully for both synthetic and field data sets. A subsequent issue that needed to be addressed involved assessing the potential benefits and implications of the resulting porosity realizations in terms of groundwater flow and contaminant transport. This was investigated synthetically assuming first that the relationship between porosity and hydraulic conductivity was well-defined. Then, the relationship was itself investigated in the context of a calibration procedure using hypothetical tracer test data. Essentially, the relationship best predicting the observed tracer test measurements was determined given the geophysically derived porosity structure. Both of these investigations showed that the SA-based approach, in general, allows much more reliable hydrological predictions than other more elementary techniques considered. Further, the developed calibration procedure was seen to be very effective, even at the scale of tomographic resolution, for predictions of transport. This also held true at locations within the aquifer where only geophysical data were available. This is significant because the acquisition of hydrological tracer test measurements is clearly more complicated and expensive than the acquisition of geophysical measurements. Although the above methodologies were tested using porosity logs and GPR data, the findings are expected to remain valid for a large number of pertinent combinations of geophysical and borehole log data of comparable resolution and sensitivity to the hydrological target parameter. Moreover, the obtained results allow us to have confidence for future developments in integration methodologies for geophysical and hydrological data to improve the 3-D estimation of hydrological properties.

Geophysical and geological analysis of the Collingwood Member of the Trenton Formation

The Collingwood Member is a mid to late Ordovician self-sourced reservoir deposited across the northern Michigan Basin and parts of Ontario, Canada. Although it had been previously studied in Canada, there has been relatively little data available from the Michigan subsurface. Recent commercial interest in the Collingwood has resulted in the drilling and production of several wells in the state of Michigan. An analysis of core samples, measured laboratory data, and petrophysical logs has yielded both a quantitative and qualitative understanding of the formation in the Michigan Basin. The Collingwood is a low permeability and low porosity carbonate package that is very high in organic content. It is composed primarily of a uniformly fine grained carbonate matrix with lesser amounts of kerogen, silica, and clays. The kerogen content of the Collingwood is finely dispersed in the clay and carbonate mineral phases. Geochemical and production data show that both oil and gas phases are present based on regional thermal maturity. The deposit is richest in the north-central part of the basin with thickest deposition and highest organic content. The Collingwood is a fairly thin deposit and vertical fractures may very easily extend into the surrounding formations. Completion and treatment techniques should be designed around these parameters to enhance production.

Geophysical logging on sediment core and hole CRP-3 from the Ross Sea, Antarctica

Cape Roberts drillhole CRP-3 in the northern part of McMurdo Sound (Ross Sea, Antarctica) targeted the western margin of the Victoria Land basin to investigate Neogene to Palaeogene climatic and tectonic history by obtaining continuous core and downhole logs (Cape Roberts Science Team, 2000). The CRP-3 drillhole extended to 939.42 mbsf (meters below seafloor) at a water depth of 297 m. The first downhole measurements after drilling were the temperature and salinity logs. Both were measured at the beginning and at the end of each of the three logging phases. Although an equilibrium temperature state may not have been fully reached after drilling, the temperature and salinity profiles seem to be scarcely disturbed. The average overall temperature gradient calculated from all temperature measurements is 28.5 K/km; remarkably lower than the temperature gradients found in other boreholes in the western Ross See and the Transantarctic Mountains. Anomalies in the salinity profiles at the beginning of each logging phase were no longer present at the end of the corresponding logging phase. This pattern indicates that drilling mud invaded the formation during drilling operations and flowed back into the borehole after drilling ceased. Thus, zones of temperature and salinity anomalies identify permeable zones in the formation and may be pathways for fluid flow. Radiogenic heat production, calculated from the radionuclide contents, is relatively low, with average values between 0.5 and 1.0 pW/m3. The highest values (up to 2 µW/m3) were obtained for the lower part of the Beacon Sandstone below 855 mbsf. The heat flow component due to radiogenic heat production integrated over the entire borehole is 0.7 mW/m2. Thermal conductivities range from 1.3 to 3 W/mK with an average value of 2.1 W/mK over the Tertiary section. Together with the average temperature gradient of 28.5 K/km this yields an average heat flow value of 60 mW/m2.

Core derived downhole logs for holes CRP-2 and CRP-2A

In the northern McMurdo Sound (Ross Sea, Antarctica), the CRP-2/2A drillhole targeted the western margin of the Victoria Land Basin to investigate Neogene to Palaeogene climatic and tectonic history by obtaining continuous core and downhole logs. Well logging of CRP-2/2A has provided a complete and comprehensive dataset of in situ geophysical measurements. This paper describes the evaluation and interpretation of the downhole logging data using multivariate statistical methods. Two major types of multivariate statistical methods were each yielding a different perspective: (1) Factor analysis was used as an objective tool for classification of the drilled sequence based on physical and chemical properties. The factor logs are mirroring the basic geological controls (i.e., grain size, porosity, clay mineralogy) behind the measured geophysical properties, thereby making them easier to interpret geologically. (2) Cluster analysis of the logs groups similar downhole geophysical properties into one cluster, delineating individual logging or sedimentological units. These objectively and independently defined units, or statistical electrofacies, are helpful in differentiating lithological and sedimentological characterisations (e.g. grain size, provenance). The multivariate statistical methods of factor and cluster analysis proved to be powerful tools for fast, reliable, and objective characterisation of downhole geophysical properties at CRP-2/2A, resulting in interpretations which are consistent with sedimentological findings.

Determination of in situ uniaxial compressive strength of coal seams based on geophysical data

The central issue for pillar design in underground coal mining is the in situ uniaxial compressive strength (sigma (cm)). The paper proposes a new method for estimating in situ uniaxial compressive strength in coal seams based on laboratory strength and P wave propagation velocity. It describes the collection of samples in the Bonito coal seam, Fontanella Mine, southern Brazil, the techniques used for the structural mapping of the coal seam and determination of seismic wave propagation velocity as well as the laboratory procedures used to determine the strength and ultrasonic wave velocity. The results obtained using the new methodology are compared with those from seven other techniques for estimating in situ rock mass uniaxial compressive strength.

Study of the Saldanha Massif (MAR, 36 degrees 34 ' N): Constrains from rock magnetic and geophysical data

We present a study of the magnetic properties of a group of basalt samples from the Saldanha Massif (Mid-Atlantic Ridge - MAR - 36degrees 33' 54" N, 33degrees 26' W), and we set out to interpret these properties in the tectono-magmatic framework of this sector of the MAR. Most samples have low magnetic anisotropy and magnetic minerals of single domain grain size, typical of rapid cooling. The thermomagnetic study mostly shows two different susceptibility peaks. The high temperature peak is related to mineralogical alteration due to heating. The low temperature peak shows a distinction between three different stages of low temperature oxidation: the presence of titanomagnetite, titanomagnetite and titanomaghemite, and exclusively of titanomaghemite. Based on established empirical relationships between Curie temperature and degree of oxidation, the latter is tentatively deduced for all samples. Finally, swath bathymetry and sidescan sonar data combined with dive observations show that the Saldanha Massif is located over an exposed section of upper mantle rocks interpreted to be the result of detachment tectonics. Basalt samples inside the detachment zone often have higher than expected oxidation rates; this effect can be explained by the higher permeability caused by the detachment fault activity.

Descoberta de conhecimento em Logs de tentativas de intrusão. Um estudo de caso em Instituições de Ensino Superior

Perante a evolução constante da Internet, a sua utilização é quase obrigatória. Através da web, é possível conferir extractos bancários, fazer compras em países longínquos, pagar serviços sem sair de casa, entre muitos outros. Há inúmeras alternativas de utilização desta rede. Ao se tornar tão útil e próxima das pessoas, estas começaram também a ganhar mais conhecimentos informáticos. Na Internet, estão também publicados vários guias para intrusão ilícita em sistemas, assim como manuais para outras práticas criminosas. Este tipo de informação, aliado à crescente capacidade informática do utilizador, teve como resultado uma alteração nos paradigmas de segurança informática actual. Actualmente, em segurança informática a preocupação com o hardware é menor, sendo o principal objectivo a salvaguarda dos dados e continuidade dos serviços. Isto deve-se fundamentalmente à dependência das organizações nos seus dados digitais e, cada vez mais, dos serviços que disponibilizam online. Dada a mudança dos perigos e do que se pretende proteger, também os mecanismos de segurança devem ser alterados. Torna-se necessário conhecer o atacante, podendo prever o que o motiva e o que pretende atacar. Neste contexto, propôs-se a implementação de sistemas de registo de tentativas de acesso ilícitas em cinco instituições de ensino superior e posterior análise da informação recolhida com auxílio de técnicas de data mining (mineração de dados). Esta solução é pouco utilizada com este intuito em investigação, pelo que foi necessário procurar analogias com outras áreas de aplicação para recolher documentação relevante para a sua implementação. A solução resultante revelou-se eficaz, tendo levado ao desenvolvimento de uma aplicação de fusão de logs das aplicações Honeyd e Snort (responsável também pelo seu tratamento, preparação e disponibilização num ficheiro Comma Separated Values (CSV), acrescentando conhecimento sobre o que se pode obter estatisticamente e revelando características úteis e previamente desconhecidas dos atacantes. Este conhecimento pode ser utilizado por um administrador de sistemas para melhorar o desempenho dos seus mecanismos de segurança, tais como firewalls e Intrusion Detection Systems (IDS).

Evolução tectónica mesozóica da Bacia Lusitaniana

The evolution of the Lusitanian Basin, localized on the western Iberian margin, is closely associated with the first opening phases of the North Atlantic. It persisted from the Late Triassic to the Early Cretaceous, more precisely until the end of the Early Aptian, and its evolution was conditioned by inherited structures from the variscan basement. The part played by the faults that establish its boundaries, as regards the geometric and kinematic evolution and the organization of the sedimentary bodies, is discussed here, as well as with respect to important faults transversal to the Basin. A basin evolution model is proposed consisting of four rifting episodes which show: i) periods of symmetrical (horst and graben organization) and asymmetrical (half graben organization) geometric evolution; ii) diachronous fracturing; iii) rotation of the main extensional direction; iv) rooting in the variscan basement of the main faults of the basin (predominantly thick skinned style). The analysis and regional comparison, particularly with the Algarve Basin, of the time intervals represented by important basin scale hiatuses near to the renovation of the rifting episodes, have led to assume the occurrence of early tectonic inversions (Callovian–Oxfordian and Tithonian–Berriasian). The latter, however, had a subsequent evolution distinct from the first: there is no subsidence renovation, which is discussed here, and it is related to a magmatic event. Although the Lusitanian Basin is located on a rift margin which is considered non-volcanic, the three magmatic cycles as defined by many authors, particularly the second (approx. 130 to 110 My ?), performed a fundamental part in the mobilization of the Hettangian evaporites, resulting in the main diapiric events of the Lusitanian Basin. The manner and time in which the basin definitely ends its evolution (Early Aptian) is discussed here. Comparisons are established with other west Iberian margin basins and with Newfoundland basins. A model of oceanization of this area of the North Atlantic is also presented, consisting of two events separated by approximately 10 My, and of distinct areas separated by the Nazaré fault. The elaboration of this synthesis was based on: - information contained in previously published papers (1990 – 2000); - field-work carried out over the last years, the results of which have not yet been published; - information gathered from the reinterpretation of geological mapping and geophysical (seismic and well logs) elements, and from generic literature concerning the Mesozoic of the west iberian margin.

Learning to manage in knowledge-based organizations: the role of learning logs or portfolios

Paper presented at the 8th European Conference on Knowledge Management, Barcelona, 6-7 Sep. 2008 URL: http://www.academic-conferences.org/eckm/eckm2007/eckm07-home.htm

Mapping the Pasture Steppe in Bayankhongor, Mongolia: comparison of classification methods, using Landsat-8 and geophysical data

Grasslands in semi-arid regions, like Mongolian steppes, are facing desertification and degradation processes, due to climate change. Mongolia’s main economic activity consists on an extensive livestock production and, therefore, it is a concerning matter for the decision makers. Remote sensing and Geographic Information Systems provide the tools for advanced ecosystem management and have been widely used for monitoring and management of pasture resources. This study investigates which is the higher thematic detail that is possible to achieve through remote sensing, to map the steppe vegetation, using medium resolution earth observation imagery in three districts (soums) of Mongolia: Dzag, Buutsagaan and Khureemaral. After considering different thematic levels of detail for classifying the steppe vegetation, the existent pasture types within the steppe were chosen to be mapped. In order to investigate which combination of data sets yields the best results and which classification algorithm is more suitable for incorporating these data sets, a comparison between different classification methods were tested for the study area. Sixteen classifications were performed using different combinations of estimators, Landsat-8 (spectral bands and Landsat-8 NDVI-derived) and geophysical data (elevation, mean annual precipitation and mean annual temperature) using two classification algorithms, maximum likelihood and decision tree. Results showed that the best performing model was the one that incorporated Landsat-8 bands with mean annual precipitation and mean annual temperature (Model 13), using the decision tree. For maximum likelihood, the model that incorporated Landsat-8 bands with mean annual precipitation (Model 5) and the one that incorporated Landsat-8 bands with mean annual precipitation and mean annual temperature (Model 13), achieved the higher accuracies for this algorithm. The decision tree models consistently outperformed the maximum likelihood ones.

Permeability estimates from poro-elastic dispersion analyses of sonic logs

Modern sonic logging tools designed for shallow environmental and engineering applications allow for P-wave phase velocity measurements over a wide frequency band. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and source frequencies ranging from approximately 1 to 30 kHz, the observable poro-elastic P-wave velocity dispersion is sufficiently pronounced to allow for reliable first-order estimations of the underlying permeability structure. These predictions have been tested on and verified for a surficial alluvial aquifer. Our results indicate that, even without any further calibration, the thus obtained permeability estimates as well as their variabilities within the pertinent lithological units are remarkably close to those expected based on the corresponding granulometric characteristics.

Control of landslide retrogression by discontinuities: Evidence by the integration of airbone- and ground-based geophysical information

The objective of this work is to present a multitechnique approach to define the geometry, the kinematics, and the failure mechanism of a retrogressive large landslide (upper part of the La Valette landslide, South French Alps) by the combination of airborne and terrestrial laser scanning data and ground-based seismic tomography data. The advantage of combining different methods is to constrain the geometrical and failure mechanism models by integrating different sources of information. Because of an important point density at the ground surface (4. 1 points m?2), a small laser footprint (0.09 m) and an accurate three-dimensional positioning (0.07 m), airborne laser scanning data are adapted as a source of information to analyze morphological structures at the surface. Seismic tomography surveys (P-wave and S-wave velocities) may highlight the presence of low-seismic-velocity zones that characterize the presence of dense fracture networks at the subsurface. The surface displacements measured from the terrestrial laser scanning data over a period of 2 years (May 2008?May 2010) allow one to quantify the landslide activity at the direct vicinity of the identified discontinuities. An important subsidence of the crown area with an average subsidence rate of 3.07 m?year?1 is determined. The displacement directions indicate that the retrogression is controlled structurally by the preexisting discontinuities. A conceptual structural model is proposed to explain the failure mechanism and the retrogressive evolution of the main scarp. Uphill, the crown area is affected by planar sliding included in a deeper wedge failure system constrained by two preexisting fractures. Downhill, the landslide body acts as a buttress for the upper part. Consequently, the progression of the landslide body downhill allows the development of dip-slope failures, and coherent blocks start sliding along planar discontinuities. The volume of the failed mass in the crown area is estimated at 500,000 m3 with the sloping local base level method.

Combining geophysical and laboratory measurements for civil Engineering.

The study of wave propagation at sonic frequency in soil leads to elasticity parameter determination. These parameters are compatible to those measured simultaneously by static loading. The acquisition of in situ elasticity parameter combined with laboratory description of the elastoplastic behaviour can lead to in situ elastoplastic curves. - L'étude de la propagation des ondes acoustiques permet la détermination des paramètres d'élasticité dans les sols. Ces paramètres sont cohérents avec des mesures statiques simultanées. L'acquisition des paramètres d'élasticité in situ associée à une description du comportement élasto-plastique mesuré en laboratoire permet d'obtenir des courbes d'élastoplasticité in situ.