Hydrologic testing during drilling: application of the flowing fluid electrical conductivity (FFEC) logging method to drilling of a deep borehole

Drilling of a deep borehole does not normally allow for hydrologic testing during the drilling period. It is only done when drilling experiences a large loss (or high return) of drilling fluid due to penetration of a large-transmissivity zone. The paper proposes the possibility of conducting flowing fluid electrical conductivity (FFEC) logging during the drilling period, with negligible impact on the drilling schedule, yet providing important information on depth locations of both high- and low-transmissivity zones and their hydraulic properties. The information can be used to guide downhole fluid sampling and post-drilling detailed testing of the borehole. The method has been applied to the drilling of a 2,500-m borehole at Åre, central Sweden, firstly when the drilling reached 1,600 m, and then when the drilling reached the target depth of 2,500 m. Results unveil eight hydraulically active zones from 300 m down to borehole bottom, with depths determined to within the order of a meter. Further, the first set of data allows the estimation of hydraulic transmissivity values of the six hydraulically conductive zones found from 300 to 1,600 m, which are very low and range over one order of magnitude.     

Paleoenvironmental reconstruction of Holocene coastal sediments adjacent to archaeological ruins in Aqaba,Jordan

Stratified cultural remains from the Early Roman/Nabataean to Byzantine periods in the coastal zone of Aqaba, Jordan, and analyses of thirteen sediment cores provide evidence for changes in the depositional environment during the Holocene. The overall trend in subsurface sediments is a basal marine transgressive layer overlain by a regressional sequence of embayment lagoonal sediments identified from microfossil analyses, and backshore pond, alluvium, and eolian deposits until the 1st century B.C., when mudbrick structures appear. Based on two radiocarbon dates, a brackish water coastal embayment formed prior to ca. 5900–5700 B.C. and was subsequently filled by siltation. Local tectonic subsidence along faults of the Dead Sea transform may have helped form the lagoon. Freshwater Candona sp. ostracods found in sand layers in the lagoon facies show signs of having been transported and redeposited. These data indicate that lakes or marshes were likely located inland of the study area. Supporting faunal and floral evidence for wetter climatic conditions, wetland habitats, and extensive water diversion and agricultural terraces has been excavated at late Chalcolithic (4000–3500 B.C.) sites in the Aqaba region (Khalil & Schmidt, 2009 ). By the 8th century B.C., the depositional environment along the coastal plain of Aqaba was dominated by distal alluvial fan and eolian sedimentation and the shoreline had prograded about 400m seaward. The migration of human settlements since the 8th century B.C. from the center of the valley toward the southeast may be driven by the changing course of Wadi Yutim and conditions along the coastal plain. © 2010 Wiley Periodicals, Inc.     

Integrated geophysical and geochemical study on AMD generation at the Haveri Au–Cu mine tailings,SW Finland

The Haveri tailings area contains 1.5 Mt of sulfide-bearing waste from the Au–Cu mine that operated during 1942–1961. Geophysical and geochemical methods were used to evaluate and characterize the generation of acid mine drainage (AMD). Correlations were examined among the electrical resistivity tomography (ERT) data, the total sulfide content and concentrations of sulfide-bound metals (Cu, Co, Fe, Mn, Ni, Pb and Zn) of tailings samples, and the resistivity and geochemistry of surface water. The resulting geophysical–geochemical model defines an area in the vadose tailings, where a low resistivity anomaly (<10 Ohm m) is correlated with the highest sulfide content, extensive sulfide oxidation and low pH (average 3.1). The physical and geochemical conditions, resulting from the oxidation of the sulfide minerals, suggest that the low resistivity anomaly is associated with acidic and metal-rich porewater (i.e., AMD). The lower resistivity values in the saturated zone of the central impoundment suggest the formation of a plume of AMD. The natural subsoil layer (silt and clay) and the bedrock surface below the tailings area were well mapped from the ERT data. The detected fracture zones of the bedrock that could work as leakage pathways for AMD were consistent with previous geological studies. The integrated methodology of the study offers a promising approach to fast and reliable monitoring of areas of potential AMD generation and its subsurface movement over large areas (ca. 9 ha). This methodology could be helpful in planning drill core sampling locations for geochemical and mineralogical analysis, groundwater sampling, and choosing and monitoring remedial programs.     

Active tectonic morphology and submarine deformation of the northern Gulf of Eilat/Aqaba from analyses of multibeam data

A high-resolution marine geophysical study was conducted during October-November 2006 in the northern Gulf of Aqaba/Eilat, providing the first multibeam imaging of the seafloor across the entire gulf head spanning both Israeli and Jordanian territorial waters. Analyses of the seafloor morphology show that the gulf head can be subdivided into the Eilat and Aqaba subbasins separated by the north-south-trending Ayla high. The Aqaba submarine basin appears starved of sediment supply, apparently causing erosion and a landward retreat of the shelf edge. Along the eastern border of this subbasin, the shelf is largely absent and its margin is influenced by the Aqaba Fault zone that forms a steep slope partially covered by sedimentary fan deltas from the adjacent ephemeral drainages. The Eilat subbasin, west of the Ayla high, receives a large amount of sediment derived from the extensive drainage basins of the Arava Valley (Wadi ’Arabah) and Yutim River to the north–northeast. These sediments and those entering from canyons on the south-western border of this subbasin are transported to the deep basin by turbidity currents and gravity slides, forming the Arava submarine fan. Large detached blocks and collapsed walls of submarine canyons and the western gulf margin indicate that mass wasting may be triggered by seismic activity. Seafloor lineaments defined by slope gradient analyses suggest that the Eilat Canyon and the boundaries of the Ayla high align along north- to northwest-striking fault systems—the Evrona Fault zone to the west and the Ayla Fault zone to the east. The shelf–slope break that lies along the 100 m isobath in the Eilat subbasin, and shallower (70–80 m isobaths) in the Aqaba subbasin, is offset by approx. 150 m along the eastern edge of the Ayla high. This offset might be the result of horizontal and vertical movements along what we call the Ayla Fault on the east side of the structure. Remnants of two marine terraces at 100 m and approx. 150 m water depths line the southwest margin of the gulf. These terraces are truncated by faulting along their northern end. Fossil coral reefs, which have a similar morphological appearance to the present-day, basin margin reefs, crop out along these deeper submarine terraces and along the shelf–slope break. One fossil reef is exposed on the shelf across the Ayla high at about 60–63 m water depth but is either covered or eroded in the adjacent subbasins. The offshore extension of the Evrona Fault offsets a fossil reef along the shelf and extends south of the canyon to linear fractures on the deep basin floor.     

Understanding the effect of single-fracture heterogeneity from single-well injection-withdrawal (SWIW) tests

The single-well injection-withdrawal (SWIW) tracer test is a method used to estimate the tracer retardation properties of a fracture or fracture zone. The effects of single-fracture aperture heterogeneity on SWIW-test tracer breakthrough curves are examined by numerical modelling. The effects of the matrix diffusion and sorption are accounted for by using a particle tracking method through the addition of a time delay added to the advective transport time. For a given diffusion and sorption property (P _m) value and for a heterogeneous fracture, the peak concentration is larger compared to a homogeneous fracture. The cumulative breakthrough curve for a heterogeneous fracture is similar to that for a homogeneous fracture and a less sorptive/diffusive tracer. It is demonstrated that the fracture area that meets the flowing water, the specific flow-wetted surface (sFWS) of the fracture, can be determined by matching the observed breakthrough curve for a heterogeneous fracture to that for a homogeneous fracture with an equivalent property parameter. SWIW tests are also simulated with a regional pressure gradient present. The results point to the possibility of distinguishing the effect of the regional pressure gradient from that of diffusion through the use of multiple tracers with different P _m values.     

Review: The state-of-art of sparse channel models and their applicability to performance assessment of radioactive waste repositories in fractured crystalline formations

Laboratory and field experiments done on fractured rock show that flow and solute transport often occur along flow channels. ‘Sparse channels’ refers to the case where these channels are characterised by flow in long flow paths separated from each other by large spacings relative to the size of flow domain. A literature study is presented that brings together information useful to assess whether a sparse-channel network concept is an appropriate representation of the flow system in tight fractured rock of low transmissivity, such as that around a nuclear waste repository in deep crystalline rocks. A number of observations are made in this review. First, conventional fracture network models may lead to inaccurate results for flow and solute transport in tight fractured rocks. Secondly, a flow dimension of 1, as determined by the analysis of pressure data in well testing, may be indicative of channelised flow, but such interpretation is not unique or definitive. Thirdly, in sparse channels, the percolation may be more influenced by the fracture shape than the fracture size and orientation but further studies are needed. Fourthly, the migration of radionuclides from a waste canister in a repository to the biosphere may be strongly influenced by the type of model used (e.g. discrete fracture network, channel model). Fifthly, the determination of appropriateness of representing an in situ flow system by a sparse-channel network model needs parameters usually neglected in site characterisation, such as the density of channels or fracture intersections.     

Flow dimension as an indicator of hydraulic behavior in site characterization of fractured rock

We examine the possibility of using the flow dimension identified from constant pressure injection tests as a tool for characterizing the hydraulic conditions of fractured media. The data comes from a low-conductivity crystalline rock site, from depths of up to 450 m, and is obtained with 2 m and 10 m measurement scales. In the analysis, the general solution for n-dimensional flow by Barker (1988) is applied. The results show that the most prominent characteristics of the medium can be identified; that is, linear and sublinear flow dimensions as distinguished from dimensions higher than two. In many cases, however, there is significant difficulty in distinguishing the dimensions n = 2, 2.5, and 3 from each other. This is usually because of the experimental difficulties in achieving the ideal conditions required by the theory during the early part of the experiment. In such cases, a full flow curve is not available for the type-curve fitting. In the nonunique cases the higher dimensions typically correspond to higher, sometimes unrealistically high, values of specific storage and to the less reliable and less representative early part of the experiment. Therefore, most of the dimensions in categories n = 3 can be excluded, thus leaving the majority observations in the categories of n = 2 and n = 2-2.5. The dominance of dimension n = 2 is more pronounced for data related to fracture zones in comparison to that related to "average" rock, in particular in the 2 m scale data. The proportion of low (n < 1.5) flow dimensions is small, but for the 10 m scale data it is relatively higher at greater depths and corresponds to lower conductivities. For the smaller 2 m scale data, the low dimensions are not linked to greater depths or systematically smaller conductivities, giving preliminary indication of different flow dimension behavior for the two different scales.     

Late Quaternary deformation and slip rates in the northern San Andreas fault zone at Olema Valley, Marin County, California

Quaternary sedimentary deposits along the structural depression of the San Andreas fault (SAF) zone north of San Francisco in Marin County provide an excellent record of rates and styles of neotectonic deformation in a location near where the greatest amount of horizontal offset was measured after the great 1906 San Francisco earthquake. A high-resolution gravity survey in the Olema Valley was used to determine the depth to bedrock and the thickness of sediment fill along and across the SAF valley. In the gravity profile across the SAF zone, Quaternary deposits are offset across the 1906 fault trace and truncated by the Western and Eastern Boundary faults, whose youthful activity was previously unknown. The gravity profile parallel to the fault valley shows a basement surface that slopes northward toward an area of present-day subsidence near the head of Tomales Bay. Surface and subsurface investigations of the late Pleistocene Olema Creek Formation (Qoc) indicate that this area of subsidence was located further south during deposition of the Qoc and that it has migrated northward since then. Localized subsidence has been replaced by localized contraction that has produced folding and uplift of the Qoc. This apparent alternation between transtension and transpression may be the result of a northward-diverging fault geometry of fault strands that includes the valley-bounding faults as well as the 1906 SAF trace. The Vedanta marsh is a smaller example of localized subsidence in the fault zone, between the 1906 SAF trace and the Western Boundary fault. Analyses of Holocene marsh sediments in cores and a paleoseismic trench indicate thickening, and probably tilting, toward the 1906 trace, consistent with coseismic deformation observed at the site following the 1906 earthquake.New age data and offset sedimentary and geomorphic features were used to calculate four late Quaternary slip rate estimates for the SAF at this latitude. Luminescence dates of 112–186 ka for the middle part of the Olema Creek Formation (Qoc), the oldest Quaternary deposit in this part of the valley, suggest a late Pleistocene slip rate of 17–35 mm/year, which replaces the unit to a position adjacent to its sediment source area. A younger alluvial fan deposit (Qqf; basal age 30 ka) is exposed in a quarry along the medial ridge of the fault valley. This fan deposit has been truncated on its western side by dextral SAF movement, and west-side-down vertical movement that has created the Vedanta marsh. Paleocurrent measurements, clast compositions, sediment facies distributions, and soil characteristics show that the Bear Valley Creek drainage, now located northwest of the site, supplied sediment to the fan, which is now being eroded. Restoration of the drainage to its previous location provides an estimated slip rate of 25 mm/year. Furthermore, the Bear Valley Creek drainage probably created a water gap located north of the Qqf deposit during the last glacial maximum 18 ka. The amount of offset between the drainage and the water gap yields an average slip rate of 21–30 mm/year. Finally, displacement of a 1000-year-old debris lobe approximately 20 m from its hillside hollow along the medial ridge indicates a minimum late Holocene slip rate of 21–25 mm/year. Similarity of the late Pleistocene rates to the Holocene slip rate, and to previous rates obtained in paleoseismic trenches in the area, indicates that the rates may not have changed over the past 30 ka, and perhaps the past 200–400 ka. Stratigraphic and structural observations also indicate that valley-bounding faults were active in the late Pleistocene and suggest the need for further study to evaluate their continued seismic potential.     

Genesis of Palaeoproterozoic iron skarns in the Misi region, northern Finland

Sodic alteration is widespread in Palaeoproterozoic greenstone and schist belts of the northern Fennoscandian shield. In the Misi region that forms the easternmost part of the Peräpohja schist belt, several small magnetite deposits show intimate spatial relationships with intensely albitised gabbros, raising the possibility that regional sodic alteration released iron, which was subsequently accumulated into deposits. Two of these magnetite deposits, Raajärvi and Puro display a typical paragenesis as follows (from oldest to youngest): (1) diopside, (2) actinolite/tremolite-magnetite ± chlorite, biotite, and (3) serpentine ± hematite, chlorite. Mass balance calculations suggest that significant amounts of Fe, Ca, Mg, K, Cu, V, and Ba were lost, and Na and Si gained during the albitisation of the gabbro, at near-constant Al, Ga, Ti, and Zr. Significant amounts of Si, Ca, Fe, and Na were enriched in the formation of skarn related to magnetite deposits. Fe and V leached from country rocks deposited during the skarn-alteration and formed the vanadium rich iron deposits while Cu passed through the system without significant precipitation due to low sulphur fugasity. Variations in Na, Ca, Mg, K, and Ba contents reflect the composition of the infiltrating fluid during alteration. Conventional heating-freezing measurements and proton-induced X-ray emission (PIXE) analyses of the fluid inclusions related to actinolite/tremolite-magnetite stage alteration indicate that the fluids that caused the alteration and the Fe-mineralisation were complex, oxidised, highly saline H₂O ± CO₂ fluids that contained high amounts of Na, Ca, K, Fe, and Ba as well as elevated concentrations of Cu, Zn, and Pb. The oxygen isotope thermometry suggest that temperature during the Fe-mineralisation stage was between 390 and 490°C. Calculated δ¹⁸O_fluid values of 6.1–9.8‰ SMOW and δ¹³C values of calcites in the ores and skarns were between ?7.7 and 10.9‰ PDB and most likely reflect admixture of ¹³C depleted, possibly magmatic fluids with the marble wall rocks that show δ¹³C_calcite values of 13‰ PDB. The SIMS U–Pb data on the zircons in the albitised gabbro next to the Raajärvi and Puro deposits suggest that intrusion of the gabbro took place at 2123±7 Ma and was accompanied by the formation of diopside skarn. The TIMS data on the metasomatic titanites related to sodic alteration yielded ages of 2062±3 and 2017±3 Ma. Iron was probably stripped from the mafic country rocks by sodic alteration between 2123 and 2017 Ma, driven by repeated brine influxes. Subsequently, the metal-rich brine was focused by a fault system and the iron was precipitated from this fluid by a combination of wall rock reaction, fluid mixing, and a drop in the temperature.