Laplace-transform analytic element solution of transient flow in porous media

A Laplace-transform analytic element method (LT-AEM) is described for the solution of transient flow problems in porous media. Following Laplace transformation of the original flow problem, the analytic element method (AEM) is used to solve the resultant time-independent modified Helmholtz equation, and the solution is inverted numerically back into the time domain. The solution is entirely general, retaining the mathematical elegance and computational efficiency of the AEM while being amenable to parallel computation. It is especially well suited for problems in which a solution is required at a limited number of points in space–time, and for problems involving materials with sharply contrasting hydraulic properties. We illustrate the LT-AEM on transient flow through a uniform confined aquifer with a circular inclusion of contrasting hydraulic conductivity and specific storage. Our results compare well with published analytical solutions in the special case of radial flow.     

Flow and fracture maps for basaltic rock deformation at high temperatures

Understanding how the strength of basaltic rock varies with the extrinsic conditions of stress state, pressure and temperature, and the intrinsic rock physical properties is fundamental to understanding the dynamics of volcanic systems. In particular it is essential to understand how rock strength at high temperatures is limited by fracture. We have collated and analysed laboratory data for basaltic rocks from over 500 rock deformation experiments and plotted these on principal stress failure maps. We have fitted an empirical flow law (Norton’s law) and a theoretical fracture criterion to these data. The principal stress failure map is a graphical representation of ductile and brittle experimental data together with flow and fracture envelopes under varying strain rate, temperature and pressure. We have used these maps to re-interpret the ductile–brittle transition in basaltic rocks at high temperatures and show, conceptually, how these failure maps can be applied to volcanic systems, using lava flows as an example.     

Earthen barriers to control lava flows in the 2001 eruption of Mt. Etna

Preceded by four days of intense seismicity and marked ground deformation, a new eruption of Mt. Etna started on 17 July and lasted until 9 August 2001. It produced lava emission and strombolian and phreatomagmatic activity from four different main vents located on a complex fracture system extending from the southeast summit cone for about 4.5 km southwards, from 3000 to 2100 m elevation (a.s.l.). The lava emitted from the lowest vent cut up an important road on the volcano and destroyed other rural roads and a few isolated country houses. Its front descended southwards to about 4 km distance from the villages of Nicolosi and Belpasso. A plan of intervention, including diversion and retaining barriers and possibly lava flow interruption, was prepared but not activated because the flow front stopped as a consequence of a decrease in the effusion rate. Extensive interventions were carried out in order to protect some important tourist facilities of the Sapienza and Mts. Silvestri zones (1900 m elevation) from being destroyed by the lava emitted from vents located at 2700 m and 2550 m elevation. Thirteen earthen barriers (with a maximum length of 370 m, height of 10–12 m, base width of 15 m and volume of 25 000 m³) were built to divert the lava flow away from the facilities towards a path implying considerably less damage. Most of the barriers were oriented diagonally (110–135°) to the direction of the flow. They were made of loose material excavated nearby and worked very nicely, resisting the thrust of the lava without any difficulty. After the interventions carried out on Mt. Etna in 1983 and in 1991–1992, those of 2001 confirm that earthen barriers can be very effective in controlling lava flows.     

Fluid flow and mineralization of Youjiang Basin in the Yunnan-Guizhou-Guangxi area, China

Comprehensive studies, based on isotope geochemistry of C, H, O, S and Sr, chronology, common element and trace element geochemistry of fluid inclusions for the epithermal Au, As, Sb and Hg deposits in the Youjiang Basin and its peripheral areas, suggested that the ore fluid was the basin fluid with abundant metallic elements and the large-scale fluid flow of the same source in the late Yenshan stage was responsible for huge epithermal mineralization and silicification. The ore fluid flowed from the basin to the platform between the basin and the platform and migrated from the inter-platform basin to the isolated platform in the Youjiang Basin. The synsedimentary faults and paleokast surface acted respectively as main conduits for vertical and lateral fluid flow.     

Velocity structure of uppermost mantle beneath North China from Pn tomography and its implications

20301 Pn arrival time data are collected from the seismological bulletins of both national and regional seismic networks. Pn travel time residuals are tomographically inverted for the Pn velocity structure of uppermost mantle beneath North China. The result indicates that the average Pn velocity in North China is 7.92 km/s, and the velocity varies laterally from ?0.21 to +0.29 km/s around the average. The approximately NNE trending high and low velocity regions arrange alternatively west-eastward. From west to east we can see high velocity in the middle Ordos region, the Shanxi graben low, the Jizhong depression high, the west Shandong uplift and Bohai Sea low, and the high velocity region to the east of the Tanlu fault. In the southern boundary zone of the North China block, except for the high velocity in the Qingling Mountains region, the velocity is generally lower than the average. Obvious velocity anisotropy is seen in the Datong Cenozoic volcanic region, with the fast velocity direction in NNE-SSW. Notable velocity anisotropy is also seen around the Bay of Bohai Sea, and the fast velocity directions seem to show a rotation pattern, possibly indicating a flow-like deformation in the uppermost mantle there. The Pn velocity variations show a reversed correlation with the Earth's heat flow. The low Pn velocity regions generally show high heat flow, e.g., the Shanxi graben and Bohai Sea region. While the high Pn velocity regions usually manifest low heat flow, e.g., the region of Jizhong depression. This indicates that the Pn velocity variation in the study region is mainly aroused by the regional temperature difference in the uppermost mantle. Strong earthquakes in the crust tend to occur in the region with the abnormal low Pn velocity, or in the transition zone between high and low Pn velocity regions. The earthquakes in the low velocity region are shallower, while that in the transition zone are deeper.     

沉积盆地异常低压与低压油气藏成藏机理综述

地下异常低压主要有两种成因：抬升—剥蚀反弹和在介质孔隙度、渗透率非均质性条件下的区域地下水稳态流动，而化学渗透与流体“冷却”在低压形成中只起次要作用。根据圈闭类型、储盖组合及成藏过程，将低压油气藏分为三种类型：①常规地层型(除砂岩透镜体外)低压油气藏，低渗透岩石通常起遮挡作用，底水与边水不发育；②砂岩透镜体低压油气藏，通常分布于盆地中心的深部，具有不含水、充满油气的特点，油气的充注和水的排出与构造抬升之前压实作用、超压引起的水驱裂缝和毛细管力的作用有关，抬升—剥蚀引起的异常低压导致水由砂岩向页岩的流动有助于油气藏中水的排出；③深盆区低渗透储层低压气藏，通常分布在含水层的下倾方向(气水倒置)，异常低压是由于构造抬升致使超压向低压演化的结果。实例研究表明，构造抬升盆地中的低压系统是一个水动力相对封闭的体系，有利于油气的聚集与保存。     

Digital model of corrected accumulated flow for peak discharge data acquisition and drainage system design

The design of a drainage system for a roofing slate quarry was implemented by the enhancement of discharge peak estimation, and the uncertainty inevitably associated with the engineering model was reduced.

The development of a topographical, geological, and vegetation cover database developed from a Geographical Information System (GIS) allowed for the definition of the drainage network for a hydraulic system, along with the calculation of the runoff coefficient. This is applied to the digital model of accumulated flow (DMF) as a weight correction coefficient, using a matrix-based model at 5×5 m resolution. The new digital model of corrected accumulated flow (DMCF) is the result of combining the thematic maps with the map of slope <3%, which was previously created from the slope model. It is demonstrated that this new model allows to apply the “Rational Method” on cartographic units defined by the GIS.

The DMCF is compared with other traditional applications of the Rational Method based on the calculation of the discharge peak considering: (1) the drainage basin as a single watershed or (2) defining an average runoff coefficient in each sub-watershed. Both approaches have bigger discharge peaks than those obtained by the DMCF since the slope, lithology, and vegetation cover have average values, and the runoff coefficient is poorly defined, increasing the uncertainty in the discharge peak.     

Hazard Zoning for Landslides Connected to Torrential Floods in the Jerte Valley (Spain) by using GIS Techniques

The Jerte Valley is anortheast-southwest tending graben located in the mountainous region of west central Spain (Spanish Central System). Mass movements have been a predominant shaping process on the Valley slopes during the Quaternary. Present day activity is characterized as either `first-time failure' (shallow debris slides and debris flows) or `reactivations' of pre-existing landslides deposits.A delineation of landslide hazard zoningwithin the Valley has been carried out by using the detailed documentation of a particular event (a debris slide and a sequel torrential flood, which occurred on the Jubaguerra stream gorge), and GIS techniques. The procedure has had four stages, which are: (1) the elaboration of a susceptibility map (spatial prediction) of landslides; (2) the elaboration of a map of `restricted susceptibility' in the particular case of slopes that are connected to streams and torrents (gorges); (3) the elaboration of a digital model which relates the altitude to the occurrence probability of those particular precipitation conditions which characterized the Jubaguerra event and (4) the combination of the probability model with the `restricted susceptibility map', to establish `critical zones' or areas which are more prone to the occurrence of phenomena that have same typology as this one.     

Generalized Coarse Graining Procedures for Flow in Porous Media

This paper focuses on heterogeneous soil conductivities and on the impact their resolution has on a solution of the piezometric head equation: owing to spatial variations of the conductivity, the flow properties at larger scales differ from those found for experiments performed at smaller scales. The method of coarse graining is proposed in order to upscale the piezometric head equation on arbitrary intermediate scales. At intermediate scales large scale fluctuations of the conductivities are resolved, whereas small scale fluctuations are smoothed by a partialy spatial filtering procedure. The filtering procedure is performed in Fourier space with the aid of a low-frequency cut-off function. We derive the partially upscaled head equations. In these equations, the impact of the small scale variability is modeled by scale dependent effective conductivities which are determined by additional differential equations. Explicit results for the scale dependent conductivity values are presented in lowest order perturbation theory. The perturbation theory contributions are summed up with using a renormalisation group analysis yielding explicit results for the effective conductivity in isotropic media. Therefore, the results are also valid for highly heterogeneous media. The results are compared with numerical simulations performed by Dykaar and Kitanidis (1992). The method of coarse graining combined by a renormalisation group analysis offers a tool to derive exact and explicit expressions for resolution dependent conductivity values. It is, e.g., relevant for the interpretation of measurement data on different scales and for reduction of grid-block resolution in numerical modeling. This revised version was published online in July 2006 with corrections to the Cover Date.     

河北省黄壁庄水库坝基渗漏单井示踪试验研究

以河北省黄壁庄水库副坝坝基的渗漏塌坑为例，利用“九五”国家重点科技成果推广项目“智能化地下水动态参数测量仪”在天然流场下，测量出水库自然渗漏量地下水的渗透流速流向、垂直向流速流向，每米含水层的渗透系数等工程必须的水文地质参数。确定性地诊断和检测出黄壁水库坝基渗漏隐患出现的具体位置和高程，从而为大坝的渗漏处理和安全处理提供详细的现场测量数据。