Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.     

Fracture analysis in the south-western Corinth rift (Greece) and implications on fault hydraulic behavior

This paper reviews the data concerning the fracture network and the hydraulic characteristics of faults in an active zone of the Gulf of Corinth. Pressure gap measured through fault planes shows that in this area the active normal faults (Aigion, Helike) act, at least temporarily and locally, as transversal seal. The analysis of the carbonate cements in the fractures on both the hangingwall and the footwall of the faults also suggests that they have acted as local seals during the whole fault zone evolution. However, the pressure and the characteristics of the water samples measured in the wells indicate that meteoric water circulates from the highest part of the relief to the coast, which means it goes through the fault zones. Field quantitative analysis and core studies from the AIG-10 well have been performed to define both regional and fault-related fracture networks. Then laboratory thin section observations have been done to recognize the different fault rocks characterizing the fault zone components. These two kinds of approach give information on the permeability characteristics of the fault zone. To synthesize the data, a schematic conceptual 3D fluid flow modeling has been performed taking into account fault zone permeability architecture, sedimentation, fluid flow, fault vertical offset and meteoric water influx, as well as compaction water flow. This modeling allows us to fit all the data with a model where the fault segments act as a seal whereas the relays between these segments allow for the regional flow from the Peloponnese topographic highs to the coast.     

Peak discharge of a Pleistocene lava-dam outburst flood in Grand Canyon, Arizona, USA

The failure of a lava dam 165,000 yr ago produced the largest known flood on the Colorado River in Grand Canyon. The Hyaloclastite Dam was up to 366 m high, and geochemical evidence linked this structure to outburst-flood deposits that occurred for 32 km downstream. Using the Hyaloclastite outburst-flood deposits as paleostage indicators, we used dam-failure and unsteady flow modeling to estimate a peak discharge and flow hydrograph. Failure of the Hyaloclastite Dam released a maximum 11 × 10⁹ m³ of water in 31 h. Peak discharges, estimated from uncertainty in channel geometry, dam height, and hydraulic characteristics, ranged from 2.3 to 5.3 × 10⁵ m³ s⁻¹ for the Hyaloclastite outburst flood. This discharge is an order of magnitude greater than the largest known discharge on the Colorado River (1.4 × 10⁴ m³ s⁻¹) and the largest peak discharge resulting from failure of a constructed dam in the USA (6.5 × 10⁴ m³ s⁻¹). Moreover, the Hyaloclastite outburst flood is the oldest documented Quaternary flood and one of the largest to have occurred in the continental USA. The peak discharge for this flood ranks in the top 30 floods (>10⁵ m³ s⁻¹) known worldwide and in the top ten largest floods in North America.     

Attenuation of aqueous metal transport at two natural acid rock drainage occurrences in the Yukon Territory, Canada

Naturally acidic drainage associated with pyritic black shale has been observed in many locations in the Yukon Territory. While not necessarily linked to known mineral deposits, most of these natural acid rock drainage occurrences show elevated dissolved concentrations of trace elements, especially zinc, nickel, copper, cadmium and arsenic. Based on field observations, microbial investigation, chemical analyses and geochemical modeling, the fate and transport of potentially deleterious elements at two natural acid drainage occurrences with slightly different settings are examined. The Macintosh Creek is a small, acidic stream （pH 2.98-3.40）, 2 km long, located in the Macmillan Pass area of east-central Yukon amidst known sedimentary exhalative massive sulfide mineralization but remains undisturbed by exploration activities. Its trace metal content is apparently derived from groundwater discharges, which gave as much as 5.0, 2.5, 0.7, 0.13 and 0.03 mg/L ofZn, Ni, Cu and As, respectively. Interaction and sorption reactions with algal mats, biofilms and iron oxyhydroxides appear to be the dominant mechanisms attenuating aqueous contaminant transport along the stream. Cryogenic precipitation further consolidates the ferricrete formation and reduces the mobility of the sorbed metals. The tributaries of the Engineering Creek along the Dempster Highway in northern Yukon drain through a series of dolomite, phyllite, argillite, limestone, black shale, sandstone and conglomerate with no known concentration of mineralization. In this area, the water chemistry fully reflects the local geology with acidic streams invariably associated with black shale occurrences. Groundwater seeps in the headwaters area of the km-180 Creek completely enclosed in black shale gave pH 3.0 and as much as 148, 39, 2.9 and 9.1 mg/L of Zn, Ni, Cu and As, respectively. Sorption with iron oxyhydroxide and organic matter appear to dominate the attenuation of contaminant transport along the stream. However, once entered into carbonate-dominated terrains, secondary carbonate minerals exercise additional geochemical control on the local water chemistry as a result of neutralization.     

A multi-method approach to study the stability of natural slopes and landslide susceptibility mapping

In this paper, a multi-method approach for the assessment of the stability of natural slopes and landslide hazard mapping applied to the Dakar coastal region is presented. This approach is based on the effective combination of geotechnical field and laboratory works, of GIS, and of mechanical (deterministic and numerical) stability analysis. By using this approach, valuable results were gained regarding instability factors, landslide kinematics, simulation of slope failure and coastal erosion. This led to a thorough assessment and strong reduction in the subjectivity of the slope stability and hazard assessment and to the development of an objective landslide danger map of the SW coast of Dakar. Analysis of the results shows that the slides were influenced by the geotechnical properties of the soil, the weathering, the hydrogeological situation, and the erosion by waves. The landslide susceptibility assessment based on this methodological approach has allowed for an appropriate and adequate consideration of the multiple factors affecting the stability and the optimization of planning and investment for land development in the city.     

吉林台爆破料应力应变关系的三轴试验研究

通过室内大型三轴实验，研究了吉林台水库爆破料在不同级配下的应力应变关系，得到了在一定击实功下爆破料的最大干密度随细料含量变化的规律，分析了爆破料在不同级配和不同围压下应力与应变的变化规律、轴向应变与体积应变的关系及抗剪强度变化特性．从微观的角度说明了变化规律产生的原因，得出爆破料的抗剪强度随级配的变化而变化的规律．     

Effects of detailed soil spatial information on watershed modeling across different model scales

Hydro-ecological modelers often use spatial variation of soil information derived from conventional soil surveys in simulation of hydro-ecological processes over watersheds at mesoscale (10–100 km²). Conventional soil surveys are not designed to provide the same level of spatial detail as terrain and vegetation inputs derived from digital terrain analysis and remote sensing techniques. Soil property layers derived from conventional soil surveys are often incompatible with detailed terrain and remotely sensed data due to their difference in scales. The objective of this research is to examine the effect of scale incompatibility between soil information and the detailed digital terrain data and remotely sensed information by comparing simulations of watershed processes based on the conventional soil map and those simulations based on detailed soil information across different simulation scales. The detailed soil spatial information was derived using a GIS (geographical information system), expert knowledge, and fuzzy logic based predictive mapping approach (Soil Land Inference Model, SoLIM). The Regional Hydro-Ecological Simulation System (RHESSys) is used to simulate two watershed processes: net photosynthesis and stream flow. The difference between simulation based on the conventional soil map and that based on the detailed predictive soil map at a given simulation scale is perceived to be the effect of scale incompatibility between conventional soil data and the rest of the (more detailed) data layers at that scale. Two modeling approaches were taken in this study: the lumped parameter approach and the distributed parameter approach. The results over two small watersheds indicate that the effect does not necessarily always increase or decrease as the simulation scale becomes finer or coarser. For a given watershed there seems to be a fixed scale at which the effect is consistently low for the simulated processes with both the lumped parameter approach and the distributed parameter approach.     

Groundwater-flow modeling in the Yucatan karstic aquifer, Mexico

The current conceptual model of the unconfined karstic aquifer in the Yucatan Peninsula, Mexico, is that a fresh-water lens floats above denser saline water that penetrates more than 40 km inland. The transmissivity of the aquifer is very high so the hydraulic gradient is very low, ranging from 7–10 mm/km through most of the northern part of the peninsula. The computer modeling program AQUIFER was used to investigate the regional groundwater flow in the aquifer. The karstified zone was modeled using the assumption that it acts hydraulically similar to a granular, porous medium. As part of the calibration, the following hypotheses were tested: (1) karstic features play an important role in the groundwater-flow system; (2) a ring or belt of sinkholes in the area is a manifestation of a zone of high transmissivity that facilitates the channeling of groundwater toward the Gulf of Mexico; and (3) the geologic features in the southern part of Yucatan influence the groundwater-flow system. The model shows that the Sierrita de Ticul fault, in the southwestern part of the study area, acts as a flow barrier and head values decline toward the northeast. The modeling also shows that the regional flow-system dynamics have not been altered despite the large number of pumping wells because the volume of water pumped is small compared with the volume of recharge, and the well-developed karst system of the region has a very high hydraulic conductivity. Electronic Publication     

Computational modeling of cyclic mobility and post-liquefaction site response

Under seismic excitation, liquefied clean medium to dense cohesionless soils may regain a high level of shear resistance at large shear strain excursions. This pattern of response, known as a form of cyclic mobility, has been documented by a large body of laboratory sample tests and centrifuge experiments. A plasticity-based constitutive model is developed with emphasis on simulating the cyclic mobility response mechanism and associated pattern of shear strain accumulation. This constitutive model is incorporated into a two-phase (solid–fluid), fully coupled finite element code. Calibration of the constitutive model is described, based on a unique set of laboratory triaxial tests (monotonic and cyclic) and dynamic centrifuge experiments. In this experimental series, Nevada sand at a relative density of about 40% is employed. The calibration effort focused on reproducing the salient characteristics of dynamic site response as dictated by the cyclic mobility mechanism. Finally, using the calibrated model, a numerical simulation is conducted to highlight the effect of excitation frequency content on post-liquefaction ground deformations.     

ARTIFICIAL NEURAL NETWORK CORRECTION OF ROUGH ERRORS OF OBSERVATIONS*

In the context of tower measured radiation datasets.following the correction principle meeting a diagnostic equation in data quality control and in terms of a technique for model construction on data and ANN (artificial neural network) retrieval for BP correction of radiation measurements with rough errors available,a BP model is presented.Evidence suggests that the developed model works well and is superior to a convenient multivariate linear regression model,indicating its wide applications.