An assessment of miniaturised electrical imaging equipment to monitor pollution plume evolution in scaled centrifuge modelling

This paper demonstrates the capabilities and limitations of a miniaturised electrical imaging technique (resistivity tomography) developed at Cardiff University to image contaminant plumes in scaled centrifuge models of the vadose zone. For this purpose a generic model of contaminant infiltration into unsaturated sand was designed. The imaging technique produces two-dimensional contoured plots of the resistivity distribution before and during contaminant infiltration experiments. During the experiments, dyed NaCl solution was released into the model and the change in resistivity associated with the contaminant plume evolution was imaged as a function of time and g-level. Capillary pressure was monitored constantly by matrix potential probes (tensiometers) in order to investigate the effect of capillary forces on plume evolution. Tests at 1g (static conditions) and 10g are described in this paper. Comparison of resulting two-dimensional tomography with observed plume geometry at the end of the 1g test showed this imaging technique to be highly effective.

Contaminant plume evolution in the unsaturated sand model was observed to be mainly gravity-driven, with plume migration and geometry being strongly affected by a tenfold increase in gravity in the centrifuge experiment. It is concluded that miniaturised electrical imaging can be a useful tool for monitoring pollution plume evolution during centrifuge tests, but when plume evolution is rapid, the time taken to interrogate each array restricts the effectiveness of the technique in monitoring changes in plume geometry. However, in such cases, resistivity tomography does provide valuable information on residual levels of contaminant fluid retained within the soil after passage of the plume.     

Evaluation of landfill disposal boundary by means of electrical resistivity imaging

This paper deals with an employment of electrical resistivity imaging (ERI) for survey of leachate content on the waste disposal site in Northern Israel. The research consisted of conducting ten ERI lines and drilling investigation wells. Data simulation used a 2D EarthImager inversion program. Analysis of 2D ERI interpretation results shows that determination of the boundary between the landfill body bottom intensively saturated with leachates and underlying layers of highly water saturated fat nonconsolidated clays presents a challenge. However, statistical analysis of ERI data indicates that standard deviation and confidence interval of a set of resistivity data measured in the landfill body are significantly larger than those in underlying clays. Moreover, maximum changes of these parameters are found on the boundary between landfill body and underlying soil, thus reflecting natural differences in scattering of resistivity data measured in these two objects.     

Application of physical modelling of debris flow triggering to field conditions: Limitations posed by boundary conditions

Debris flows are often triggered by Hortonian overland flow during high-intensity rainstorms. Data derived from debris flow trigger zones in the southern French Alps were fed into a physical model of debris flow triggering based on Takahashi. Using a Monte Carlo approach with 1000 runs, the results show a wide distribution of safety factor values, indicating that physical modelling based on actual field measurements may not always be practical.As all safety factor values obtained are well below 1 even though debris flows only occur during very high-intensity rainstorms, the model used must be inappropriate. Apparently, the composition of the overland flow plays an important role: during high-intensity rainstorms it usually has a very high sediment content and contains stones. This prevents it from flowing through the pores of coarse debris accumulations in the central gully of a trigger zone; it will rather run over the debris. This situation is more stable than with the fluid flowing through the pores. The behaviour switch of the fluid above a certain sediment and stone content thus drastically changes the triggering conditions for debris flows and it is concluded that debris flow triggering in the area requires the occurrence of both overland flow and landsliding.     

高温高压下岩石电导率不同测量方法的实验对比——以二辉橄榄岩为例

利用YJ-3000t紧装式六面顶压机,在1.0GPa、500～900℃和2.0GPa、500～950℃的实验条件下,同时采用交流阻抗谱法(AC0.05～10-6Hz)、直流法(DC)和单频交流法(0.1Hz)三种方法对采自河北大麻坪玄武岩中的尖晶石二辉橄榄岩进行了电导率的实验测量.实验结果表明:在实验的温度和压力范围内...     

高密度电法在柳州太阳村镇岩溶塌陷区调查中的应用研究

为了查明柳州太阳村镇岩溶塌陷形成的地质背景,本文采用高密度电法对该塌陷区进行调查研究。研究结果表明:高密度电法对于确定土层的厚度,区分土层和岩石的分界面具有较好的效果,通过电阻率的分布及变化趋势,可以较好的圈定地下岩溶裂隙带的发育位置,结合地质资料,对识别断裂破碎带也有较好的效果。三个塌陷区的致塌地质背景不尽相同:第Ⅰ,Ⅲ塌陷区塌陷的形成与人类活动密切相关,分别由机井抽水及输油管道铺设引发;第Ⅱ塌陷区塌陷的形成与自然因素相关,由断裂破碎带触发。通过高密度电法反演结果断面图,能较准确地查明塌陷的潜伏隐患区:第Ⅰ塌陷区的隐伏区为通过抽水井的近东西向的岩溶裂隙带,第Ⅱ塌陷区地质灾害的潜伏区为断层西侧的岩溶裂隙带,第Ⅲ塌陷区的隐患区为输油管道附近东侧的岩溶裂隙带,均位于人类活动频繁附近地段及构造发育带。     

Seismic imaging of lithospheric discontinuities and continental evolution

Discontinuities in physical properties within the continental lithosphere reflect a range of processes that have contributed to craton stabilization and evolution. A survey of recent seismological studies concerning lithospheric discontinuities is made in an attempt to document their essential characteristics. Results from long-period seismology are inconsistent with the presence of continuous, laterally invariant, isotropic boundaries within the upper mantle at the global scale. At regional scales, two well-defined interfaces termed H (60 km depth) and L (200 km depth) of continental affinity are identified, with the latter boundary generally exhibiting an anisotropic character. Long-range refraction profiles are frequently characterized by subcontinental mantle that exhibits a complex stratification within the top 200 km. The shallow layering of this package can behave as an imperfect waveguide giving rise to the so-called teleseismic Pn phase, while the L-discontinuity may define its lower base as the culmination of a low velocity zone. High-resolution, seismic reflection profiling provides sufficient detail in a number of cases to document the merging of mantle interfaces into lower continental crust below former collisional sutures and magmatic arcs, thus unambiguously identifying some lithospheric discontinuities with thrust faults and subducted oceanic lithosphere. Collectively, these and other seismic observations point to a continental lithosphere whose internal structure is dominated by a laterally variable, subhorizontal layering. This stratigraphy appears to be more pronounced at shallower lithospheric levels, includes dense, anisotropic layers of order 10 km in thickness, and exhibits horizontal correlation lengths comparable to the lateral dimensions of overlying crustal blocks. A model of craton evolution which relies on shallow subduction as a principal agent of craton stabilization is shown to be broadly compatible with these characteristics.     

EH-4电导率成像系统的特点及其在金属矿勘探中的应用

EH-4电导率成像系统是一种混场源频率域电磁观测系统,具有轻便、高效、勘探深度大、分辨率高的特点。通过在块状硫化物铜矿、斑岩铜矿已知矿区的试验研究,证实EH-4系统在这两种矿床类型的隐伏-半隐伏金属矿勘探中有较好的找矿效果。     

吐哈盆地地浸砂岩型铀矿成矿条件与盆地动力学演化

分析吐哈盆地气候、地貌景观、地下水动力学、大地构造及构造、层间氧化带发育条件、含矿目的层地质及地球化学特征等，认为吐哈盆地北部与南部(艾丁湖斜坡带、了南凹陷、南湖凹陷等)地浸砂岩型铀矿成矿条件存在明显差异，其根本原因是与盆地南北动力学演化的差异有关。依据盆地构造沉降、构造层序、受板块运动影响等特征，认为盆地构造演化分为5个阶段，其中第三演化阶段(C-T)挤压型前陆盆地沉积阶段奠定了盆地南北构造演化差异的基础，对中新生代地浸砂岩型铀矿的分布影响深远；第四阶段(J-K)弱伸展型断陷盆地阶段为南部地浸砂岩型铀矿的形成准备了物质基础；第五阶段新生代挤压型前陆盆地阶段(E-Q)为南部地浸砂岩型铀矿形成富集提供了充分的外部动力。     

Seismic imaging and potential field modelling to delineate structures hosting VHMS deposits in the Skellefte Ore District, northern Sweden

The Skellefte District in northern Sweden is a roughly 150 by 50 km² large early Proterozoic massive sulphide belt. Based on high-resolution reflection seismic data along two parallel seismic profiles, potential field modeling has been carried out and two geologic cross sections have been constructed that are consistent with the available geophysical data as well as surface geologic observations. The combined modelling suggests that the Kristineberg deposit occurs on the northern limb of a regional E–W striking syncline. The interpretations help to identify new prospective areas, both down-plunge from known ores, and on the ore-bearing horizon on the southern limb of the syncline. The new results suggest that the post-orogenic Revsund granites can be divided into two major types of intrusives, those which are intruded as domes/stokes with a maximum present day thickness of about 3–3.5 km and those which are intruded as thin sheets, with a maximum thickness of a few hundred meters. The margins of the intrusions are generally inclined inwards, suggesting that the current erosion level is near the middle, or toward the base, of the granites. The contact between the Skellefte volcanic rocks and the Bothnian Basin has been interpreted as a thrust fault. We also suggest that crustal thickening predates the Skellefte volcanism and that the interpreted Bothnian Basin rocks are either a structural basement or a separate terrane to the Skellefte volcanism. Diffraction patterns in the reflection seismic data can be interpreted as originating from either a mafic–ultramafic intrusion or a mineralization zone, similar to observations elsewhere in the world. The results obtained in this study have greatly improved our understanding of the tectonostratigraphic framework and architecture of the poly-deformed c. 1.9 Ga Skellefte VHMS belt and is a key step towards building a 3D geological model in the area.     

Electrical conductivity measurement of gneiss under mid- to lower crustal P–T conditions

We conducted electrical conductivity measurements perpendicular and parallel to mineral foliation in dry gneiss at up to 1000 K and a constant pressure of 1 GPa. The analyzed gneisses were obtained from the Higo metamorphic belt, Kyushu, Japan. As the metamorphic conditions of these rocks have been well determined by previous studies, we were able to select samples that were representative of the middle to lower crust. Prior to the conductivity measurements, the samples were maintained at the maximum temperature for a long period, until the electrical conductivity had stabilized. Our experiment results reveal linear and reproducible conductivity data between temperatures of 600 and 1000 K. Conductivity measured perpendicular and parallel to foliation differ by an order of magnitude over the same temperature window. A plausible explanation for this discrepancy in conductivity is the contrasting configuration of minerals in the two sample orientations, as observed by backscattered electron image (BEI) and electron probe microanalysis (EPMA). We evaluated the conductivity and computed activation energy for each of the samples and compared the results with those of previous studies; our results are consistent with the conductivity values reported for other types of rocks. We also compared the experiment results with data derived from electromagnetic (EM) soundings. Electrical conductivity measurements undertaken perpendicular to foliation can account for the subsurface conductivity structure beneath central Kyushu, Japan.     

就地应力条件下岩石静动态力学性质的研究

岩石力学性质参数可以通过模拟地层应力条件下的室内岩石力学实验得到,但它只能代表一些离散点,所以,在实际应用中的大多数情况下,石油工程师采用连续的测井曲线计算岩石力学参数,然而由于加载频率的不同,计算值与室内准静态条件下测得的参数值有一定的差异,只有校核后的测井曲线才能用于工程应用中。为了从室内实验结果中得到校核后的连续的测井曲线,对模拟就地应力条件下静动态岩石力学性质进行了初步的研究。     

非平衡吸附模型在研究渗滤液对土壤污染影响中的应用

分析了有机污染物在土壤中的迁移转化机理，建立了非平衡吸附作用下渗滤液中有机污染物在土壤中迁移转化的动力学模型，给出了模型的数值解法，模拟出渗滤液在非平衡吸附作用下的污染过程；同时还探讨了模型参数降雨量p，垃圾土土层的厚度h，含水率θ等对有机污染物运移的影响。结果表明：在总的污染源一定的情况下，降雨量的增大和污染土层厚度的增大会使得下层土壤中有机物的浓度降低，为定量研究有机污染物在土壤-水环境系统中分配与归宿提供理论依据，同时也可为监测、治理和恢复地下水污染提供一定的科学根据。     

高地应力卸荷条件下错动带塑性变形规律与硬化特征

为阐明高地应力地下工程开挖卸荷条件下错动带变形破坏的塑性规律与硬化特征,通过一系列不同围压、不同加卸荷应力路径下的三轴试验,着重研究了不同加卸荷应力路径下的错动带塑性变形规律;基于试验结果,在应力空间进一步探究了等效塑性功、等效塑性应变、塑性体应变等内变量作为硬化参数的应力路径相关性,并提出应力路径无关的硬化参量修正公式。研究结果表明:(1)高地应力下错动带的变形力学特性存在着明显的应力路径效应,在相同卸荷应力路径下,高围压会抑制错动带试样环向塑性变形的发展;在相同初始围压下,卸轴压卸围压应力路径下错动带试样的塑性体积残余应变(6%)明显大于其余卸荷应力路径(2%~4%);不同应力路径对错动带塑性体积变形的促进作用为卸轴压卸围压应力路径>恒轴压卸围压应力路径>增轴压卸围压应力路径;(2)错动带在不同卸荷路径变形破坏过程中,其塑性体应变、等效塑性应变、等效塑性功等内变量均存在一定的应力路径相关性,因而在错动带的弹塑性分析中,直接将以上任何一个状态参量作为其硬化参数并假定其与应力路径无关是不准确的。为此,提出等效塑性功修正方法,发现当修正等效塑性功公式中反映错动带材料特性的参数n_s=-0.4时,修正等效塑性功具有明显的应力路径无关性,更适宜作为高地应力卸荷条件下错动带的硬化参数来描述其卸荷塑性应变硬化特征。揭示的错动带塑性力学特性可进一步深化高地应力卸荷条件下错动带的变形破坏认识,为实际工程中错动带的破坏机制分析和支护控制提供理论依据。     

Size and shape evolution of pores in a viscoplastic matrix under compression

The frequent use of soils and earth materials for hydraulic capping and for geo‐environmental waste containment motivated our interest in detailed modelling of changes in size and shape of macro‐pores to establish links between soil mechanical behaviour and concurrent changes in hydraulic and transport properties. The objective of this study was to use finite element analysis (FEA) to test and extend previous analytical solutions proposed by the authors describing deformation of a single macro‐pore embedded in linear viscoplastic soil material subjected to anisotropic remote stress. The FEA enables to consider more complex pore geometries and provides a detailed picture of matrix yield behaviour to explain shortcomings of approximate analytical solutions. Finite element and analytical calculations agreed very well for linear viscous as well as for viscoplastic materials, only limited for the case of isotropic remote stress due to the simplifications of the analytical model related to patterns and onset of matrix‐yielding behaviour. FEA calculations were compared with experimental data obtained from a compaction experiment in which pore deformation within a uniform modelling clay sample was monitored using CAT scanning. FEA predictions based on independently measured material properties and initial pore geometry provided an excellent match with experimentally determined evolution of pore size and shape hence lending credence to the potential use of FEA for more complex pore geometries and eventually connect macro‐pore deformation with hydraulic properties. Copyright © 2006 John Wiley & Sons, Ltd.     

Application of unconstrained optimization and sensitivity analysis to calibration of a soil constitutive model

Large sets of soil experimental data (field and laboratory) are becoming increasingly available for calibration of soil constitutive models. A challenging task is to calibrate a potentially large number of model parameters to satisfactorily match many data sets simultaneously. This calibration effort can be facilitated by optimization techniques. The current study aims to explore systematic approaches for exercising optimization and sensitivity analysis in the area of soil constitutive modelling. Analytical, semi‐analytical and numerical optimization techniques are employed to calibrate a multi‐surface‐plasticity sand model. Calibration is based on results from a number of drained triaxial sample tests and a dynamic centrifuge liquefaction test. The analytical and semi‐analytical approaches and associated sensitivity analysis are applied to calibrate the model non‐linear shear stress–strain response. Thereafter, model parameters controlling shear–volume coupling effects (dilatancy) are calibrated using a solid–fluid fully coupled finite element program in conjunction with an advanced numerical optimization code. A related sensitivity study reveals the challenges often encountered in optimizing highly non‐linear functions. Overall, this study demonstrates applicability and limitations of optimization techniques for constitutive model calibration. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical modelling on water retention and permeability of compacted GMZ bentonite under free-swelling conditions

In this study, a solid–liquid–gas coupled equation was established to simulate water retention characteristics of highly compacted GMZ bentonite. Then, modelling results were compared with laboratory test results. Results indicate that GMZ bentonite has a strong moisture expansion (or a limit drying shrinkage) characteristic. The control equation can simulate the water absorption and deformation characteristics very well at high relative humidity (or low suction). Environmental scanning electron microscope (ESEM) observation reveals the course grain soil texture of the surface under low relative humidity (RH), while the surface of GMZ bentonite becomes smooth (more fine-grained soil texture) as RH increases. Differences were found between the porosities calculated by macroexperiment results and microscopic observations with ESEM method. This is because only the interaggregate pores can be observed by ESEM photographs. Additionally, we find that the simulated effective porosities are close to the results calculated by microscopic tests, while the effective porosity is considered as the main flow channel of flow. Further, the intrinsic permeability, the effective water and gas permeability are calculated based on the proposed model. The modelling results coincide well with the laboratory experimental results and support the reliability of the proposed model.     

Spatial dynamics of soil salinity under arid and semi-arid conditions: geological and environmental implications

Salinity is a highly important problem in arid and semi-arid regions. Consequently, the analysis of soil salinity plays a crucial role in environmental sciences. Environmental variables show spatial variability and statistical tools able to analyze and describe the spatial dependence and quantify the scale and intensity of the spatial variation are needed. Moreover, the spatial prediction has special importance for agricultural transformation of the land or for environmental restoration (selection of the most appropriate species adapted to soil salinity). In this paper we propose a methodological framework based on geostatistics to model the spatial dynamics of soil salinity measurements to further analyze the case of a region in southeastern Spain.     

加卸荷条件下灰岩能耗变化规律试验研究

对灰岩开展单轴、常规三轴与峰前卸荷试验,对比分析不同应力路径下灰岩变形过程的能量变化特征。试验分析表明,灰岩加荷破坏和卸荷破坏的能量变化规律存在明显差异。由于围压抑制了裂纹的扩展,与单轴压缩破坏相比,常规三轴压缩试验灰岩吸收的总能量 和积蓄的弹性能 更多。常规三轴试验灰岩积蓄的弹性能 在峰值强度前一直增加,而卸荷破坏的弹性能 在卸荷开始后基本不变,说明卸荷破坏释放的弹性能 主要是在加荷阶段累积完成的,岩体所处的初始地应力状态决定了其破坏的能量释放量。加荷和卸荷应力路径下,灰岩临近破坏点耗散能 都快速增加,但是卸荷破坏耗散能 增速远大于加荷破坏耗散能增速,耗散能的迅速增加表明岩石破坏的发生。随着围压升高,灰岩吸收的总能量 和储存的弹性能 逐渐升高。随着卸荷速度增加,灰岩吸收的总能量 、储存的弹性能 和耗散的应变能 逐渐减小。     

Improvement of a hypoplastic model to predict clay behaviour under undrained conditions

A shortcoming of the hypoplastic model for clays proposed by the first author is an incorrect prediction of the initial portion of the undrained stress path, particularly for tests on normally consolidated soils at isotropic stress states. A conceptually simple modification of this model, which overcomes this drawback, is proposed in the contribution. The modified model is applicable to both normally consolidated and overconsolidated soils and predicts the same swept-out-memory states (i.e., normal compression lines) as the original model. At anisotropic stress states and at higher overconsolidation ratios the modified model yields predictions similar to the original model.     

地球内部物质物性的原位高温高压研究:大体积压机与同步辐射源的结合

大体积压机(large-volume press,简称为LVP)与金刚石对顶砧(diamond-anvil cell,简称为DAC)是两种相辅相成的、在高压研究中不可或缺的重要技术。大体积压机在温压范围上较金刚石对顶砧小,但它具有更加均匀的温压环境以及高于金刚石对顶砧几千以至百万倍的试样体积。这些特色使大体积压机成为精确测量试样整体物性和系统研究复杂材料系列的理想工具。在过去的20年里,高压和同步辐射技术方面的进展使大体积高压研究取得了革命性的突破。文中选择近期地学界的一些重要科研成果和技术突破作为例子,力图对同步辐射和大体积高压技术在高温高压研究中的近期进展作综合评述,最后简略讨论了LVP的发展前景。