Evaluation of rock-mass permeability tensor and prediction of tunnel inflows by means of geostructural surveys and finite element seepage analysis   总被引：1，自引：0，他引：1  

N. Coli  G. Pranzini  A. Alfi  V. Boerio 《Engineering Geology》2008,101(3-4):174-184

  相似文献   

Sampling procedure for small-scale heterogeneities (crossbedding) for reservoir modelling     

D. Mikes   《Marine and Petroleum Geology》2006,23(9-10):961-977

Sampling permeability from individual laminae in crossbedded rocks is very well possible by making use of a so-called probe permeameter. The methodology and calculation of permeabilities from the output of this instrument are fairly well-established. However, what is hardly established is a sampling strategy that guarantees adequate representation of permeability contrast between laminae. Especially so, because finer foreset and bottomset are difficult to sample, yet their permeabilities are critical to fluid flow through a crossbed. Therefore we propose to systematically sample laminae. To this end we suggest a sampling strategy that is adapted to the heterogeneity under study, namely to the thickness of the finer foreset.  相似文献   

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations     

M. Gerritsen  J. V. Lambers 《Computational Geosciences》2008,12(2):193-208

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.  相似文献   

Accurate local upscaling with variable compact multipoint transmissibility calculations     

James V. Lambers  Margot G. Gerritsen  Bradley T. Mallison 《Computational Geosciences》2008,12(3):399-416

We propose a new single-phase local upscaling method that uses spatially varying multipoint transmissibility calculations. The method is demonstrated on two-dimensional Cartesian and adaptive Cartesian grids. For each cell face in the coarse upscaled grid, we create a local fine grid region surrounding the face on which we solve two generic local flow problems. The multipoint stencils used to calculate the fluxes across coarse grid cell faces involve the six neighboring pressure values. They are required to honor the two generic flow problems. The remaining degrees of freedom are used to maximize compactness and to ensure that the flux approximation is as close as possible to being two-point. The resulting multipoint flux approximations are spatially varying (a subset of the six neighbors is adaptively chosen) and reduce to two-point expressions in cases without full-tensor anisotropy. Numerical tests show that the method significantly improves upscaling accuracy as compared to commonly used local methods and also compares favorably with a local–global upscaling method.  相似文献   

Estimation of hydraulic permeability considering the micro morphology of rocks of the borehole YAXCOPOIL-1 (Impact crater Chicxulub,Mexico)     

S. I. Mayr  H. Burkhardt  Yu. Popov  A. Wittmann 《International Journal of Earth Sciences》2008,97(2):385-399

Internal surface, formation factor, Nuclear Magnetic Resonance (NMR)-T2 relaxation times and pore radius distributions were measured on representative core samples for the estimation of hydraulic permeability. Permeability is estimated using various versions of the classic Kozeny–Carman-equation (K–C) and a further development of K–C, the fractal PaRiS-model, taking into account the internal surface. In addition to grain and pore size distribution, directly connected to permeability, internal surface reflects the internal structure (“micro morphology”). Lithologies could be grouped with respect to differences in internal surface. Most melt rich impact breccia lithologies exhibit large internal surfaces, while Tertiary post-impact sediments and Cretaceous lithologies in displaced megablocks display smaller internal surfaces. Investigations with scanning electron microscopy confirm the correlation between internal surface and micro morphology. In addition to different versions of K–C, estimations by means of NMR, pore radius distributions and some gas permeability measurements serve for cross-checking and calibration. In general, the different estimations from the independent methods and the measurements are in satisfactory accordance. For Tertiary limestones and Suevites bulk with very high porosities (up to 35%) permeabilites between 10⁻¹⁴ and 10⁻¹⁶ m² are found, whereas in lower Suevite, Cretaceous anhydrites and dolomites, bulk permeabilites are between 10⁻¹⁵ and 10⁻²³ m².  相似文献   

Permeability development in vesiculating magmas: implications for fragmentation   总被引：2，自引：2，他引：0  

Caroline Klug  Katharine V. Cashman 《Bulletin of Volcanology》1996,58(2-3):87-100

 Fragmentation, or the "coming apart" of magma during a plinian eruption, remains one of the least understood processes in volcanology, although assumptions about the timing and mechanisms of fragmentation are key parameters in all existing eruption models. Despite evidence to the contrary, most models assume that fragmentation occurs at a critical vesicularity (volume percent vesicles) of 75–83%. We propose instead that the degree to which magma is fragmented is determined by factors controlling bubble coalescence: magma viscosity, temperature, bubble size distribution, bubble shapes, and time. Bubble coalescence in vesiculating magmas creates permeability which serves to connect the dispersed gas phase. When sufficiently developed, permeability allows subsequent exsolved and expanded gas to escape, thus preserving a sufficiently interconnected region of vesicular magma as a pumice clast, rather than fully fragmenting it to ash. For this reason pumice is likely to preserve information about (a) how permeability develops and (b) the critical permeability needed to insure clast preservation. We present measurements and calculations that constrain the conditions (vesicularity, bubble size distribution, time, pressure difference, viscosity) necessary for adequate permeability to develop. We suggest that magma fragments explosively to ash when and where, in a heterogeneously vesiculating magma, these conditions are not met. Both the development of permeability by bubble wall thinning and rupture and the loss of gas through a permeable network of bubbles require time, consistent with the observation that degree of fragmentation (i.e., amount of ash) increases with increasing eruption rate. Received: 5 July 1995 / Accepted: 27 December 1995  相似文献   

Fault zone geometry of a mature active normal fault: A potential high permeability channel (Pirgaki fault, Corinth rift, Greece)   总被引：2，自引：1，他引：2  

Yves Graud  Marc Diraison  Ney Orellana 《Tectonophysics》2006,426(1-2):61

We present results from petrophysical analysis of a normal fault zone with the aim of defining the flow pathways and their behavior during seismic and interseismic periods. Data are obtained on porosity geometry, strain structure and mineralogy of different domains of a normal fault zone in the Corinth rift. Data point out a close relationship between mineralogy of the clayey minerals, porosity network and strain structures and allow definition of a macroscopic anisotropy of the flow parameters with a strong control by microscopic ultracataclasite structures. The Pirgaki fault zone, developed within pelagic limestone, has a sharp asymmetric porosity profile, with a high porosity volume in the fault core and in the damage zone of the hanging wall. From porosity volumes and threshold measurements, a matrix permeability variation of 6 orders of magnitude could be expected between the protolith and the fault core. Modifications of this pathway during seismic and interseismic phases are depicted. Healing of cracks formed during seismic slip events occurred in the fault core zone and the porous network in the damage zone is sealed in a second step. The lens geometry of the fault core zone is associated with dissolution surfaces and open conduits where dissolved matter could move out of the fault core zone. These elementary processes are developed in particular along Riedel's structures and depend on the orientation of the strain surfaces relative to the local stress and depend also on the roughness of each surface type. P-surfaces are smooth and control shearing process. R-surfaces are rough and present two wavelengths of roughness. The long one controls localization of dissolution surfaces and conduits; the short one is characteristic of dissolution surfaces. The dissolved matter can precipitate in the open structures of the hanging wall damage zone, decreasing the connectivity of the macroscopic conduit developed within this part of the fault zone.  相似文献   

Radiomagnetotelluric mapping, groundwater numerical modelling and 18-Oxygen isotopic data as combined tools to determine the hydrogeological system of a landslide prone area     

P. Tullen  P. Turberg  A. Parriaux 《Engineering Geology》2006,87(3-4):195-204

Three methods were combined to determine the groundwater recharge and transfer processes of a landslide prone area. First, the radiomagnetotelluric method was used to investigate the distribution of electrical resistivity (ρ) of the subsurface and build a three-dimensional model of permeability (k), through an experimental relation between ρ and k. Second, this structural model of permeability and additional climatologic data were used to fix boundary and recharge conditions to perform a three-dimensional and transient numerical simulation of the groundwater flow. Finally 18-Oxygen time series observed at the main springs were used to validate the model. This association of methods led to an improved characterization of the groundwater flow system at local scale and a better understanding of the role of this system on the landslide phenomenon. This structured approach is thought to be useful to design specific remediation strategies to drain the unstable mass.  相似文献   

Attenuation of landfill leachate at two uncontrolled landfills     

Jin-Yong Lee  Jeong-Yong Cheon  Hyung-Pyo Kwon  Hee-Sung Yoon  Seong-Sun Lee  Jong-Ho Kim  Joung-Ku Park  Chang-Gyun Kim 《Environmental Geology》2006,51(4):581-593

Attenuation characteristics of landfill leachate were examined for two uncontrolled landfills in Korea. The two landfills containing municipal wastes without appropriate bottom liner and leachate treatment system have different landfill age, waste volume, and most importantly different hydrogeologic settings. One landfill (Cheonan landfill) is situated in an open flat area while the other (Wonju landfill) is located in a valley. Variations of various parameters including dissolved organic carbon (DOC), dissolved oxygen (DO), alkalinity, pH, electrical conductivity (EC), redox potential (ORP), ammonia (NH₃), nitrate (NO₃⁻), sulfate (SO₄²⁻), and chloride (Cl⁻) were examined along groundwater flow path. All these parameters were analyzed every month for a year. In the interior of the landfills, typical anaerobic conditions revealed by low DO and NO₃ concentrations, negative ORP values, high NH₃, alkalinity, and Cl⁻ concentrations were observed. Generally, higher levels of contaminants (DOC, NH₃, and Cl⁻) were detected in the dry season while they were greatly lowered in the wet season. Significantly, large decrease of Cl^- concentration in the wet season indicates that the dilution or mixing is one of dominant attenuation mechanisms of leachate. But detailed variation behaviors in the two landfills are different and they were largely dependent on permeability of surface and subsurface layers. The intermediately permeable surface of the landfills receives part of direct rainfall infiltration but most rainwater is lost to fast runoff. The practically impermeable surface of clayey silt (paddy field) at immediately adjacent to the Cheonan landfill boundary prevented direct rainwater infiltration and hence redox condition of the ground waters were largely affected by that of the upper landfill and the less permeable materials beneath the paddy fields prohibited dispersion of the landfill leachate into down gradient area. In the Wonju landfill, there are three different permeability divisions, the landfill region, the sandy open field and the paddy field. Roles of the landfill and paddy regions are very similar to those at the Cheonan. The very permeable sandy field receiving a large amount of rainwater infiltration plays a key role in controlling redox condition of the down gradient area and contaminant migration. This paper reports details of the attenuation and redox conditions of the landfill leachates at the two uncontrolled landfills.  相似文献   

Fluid permeability of sedimentary rocks in a complete stress–strain process     

J. -A. Wang  H. D. Park   《Engineering Geology》2002,63(3-4):291-300

The permeability of sedimentary rocks during triaxial compression tests was investigated to relate it to the complete strain–stress process. It was found that the permeability was not constant, but varied with the stress and strain states in the rocks. Prior to the peak strength, the permeability decreases with increasing load. A dramatic increase in permeability occurs during the strain softening period. In the present study, in situ measurements of fluid flow and pressure in floor strata was carried out in a double longwall mining face in the Yangzhuang colliery. These measurements show that both the strata pressure and the position with respect to the mining face influence the hydrogeologic properties. The permeability increased in the floor strata behind the mining face because those mining induced fractures opened as the strata pressure decreased. To better understand this change in hydraulic behavior around the mining faces, 3-D numerical modeling was carried out. The model provides the general picture of the stress distribution and failure zone both in the floor and roof strata. The field and model results demonstrate the importance of changes in the stress and strain states on the hydrogeology of a site.  相似文献