电阻率层析成像技术在滑坡调查中的应用——土耳其研究实例

电阻率层析成像是一种近地表地球物理勘查方法,它广泛应用于滑坡在内的各种地质、环境和工程等地质勘查。本次研究中,对土耳其艾登市瑟凯地区的一个滑坡场地进行了电阻率层析成像（ERT）勘查。2003年,由于开挖作业,导致一所新建学校附近斜坡上的上第三系地层变得不稳定,并在一场大雨之后出现了移动。随之而来的滑坡覆盖了这所学校的部分建筑物。作者采用温纳排列形式的二维电阻率方法,沿滑坡体布设三条剖面,得到了滑坡体的几何结构和特征等有用信息。另外,为了解电阻率法对滑坡的响应曲线,在野外调查之前,利用电阻率方法对滑坡进行了二维合成电阻率模拟研究。在滑坡场地还打了八个钻孔。钻孔和电阻率层析成像勘查两样成果都指出该场地内有一断层。沿着滑坡轴的测线上电阻率数据反映了该断裂面。     

高密度电阻率成像法在水文地质领域中的应用

高密度电阻率成像法(ERT)作为一种高效、准确的地球物理探测手段,在水文地质领域得到了广泛的应用。本文通过文献追踪和分析,重点讨论了电阻率与含水率、渗透系数、含盐量等水文地质表征参数之间的相关性,并就该方法在地下水资源勘查、溶质运移动态监测以及水文地质参数反演等方面的应用进行了剖析,阐述了ERT探测过程中的主要影响因素及各因素之间的相互关系。最后,针对ERT方法在应用过程中存在的问题提出相应的改进建议。     

基于高密度电阻率法的土体干缩裂隙动态发育过程精细监测研究

土体在气候作用下发育干缩裂隙是一种常见的自然现象,裂隙的存在会极大弱化土体的工程性质,诱发许多岩土与地质工程问题。为了实时掌握黏性土中干缩裂隙网络的发展状态,提出了一套基于高密度电阻率层析成像技术（ERT）的土体干缩裂隙动态发展过程精细监测方法。分别开展模型试验及原位试验,利用自行研制的测定系统持续采集电流-电位差数据,随后利用自行开发的有限元法电阻率层析成像（FemERT）系统进行数据处理,获取了裂隙网络在不同发育阶段的空间分布特征。结果表明：（1）ERT可以实现土体裂隙发育过程的精细监测,具备监测三维裂隙网络几何形态的能力,裂隙宽度的识别精度达到毫米级,裂隙深度的识别精度达到厘米级;（2）ERT的感度分布特征解释了裂隙发育对于土体电阻率的影响规律,测定电阻值时程曲线因裂隙产生位置的不同而呈现不同的变化规律;（3）反演电阻率及其相对变化率（Rev）可以直观表征裂隙网络在不同阶段的空间几何形态,凸显裂隙动态发育过程对于土体导电性的影响。     

由电阻率计算水文地质参数的灰色模型

本文应用南京某水源地的地球物理数据(电阻率)与实际钻孔中含水层渗透系数之间存在的灰关系,建立了由电阻率计算渗透系数的灰色模型GM(1,1)。进而了解了该区含水层渗透系数的整个平面分布情况。经与勘探抽水试验成果对比表明,由此获得的资料是有实际应用价值的。     

三维电阻率视图在揭示煤矿采区水害中的应用

为满足煤田勘探、开发的需要,结合电磁法勘探的特点,利用三维空间数据场及计算机图形学技术,开发了具有灵活、方便等特点的三维电阻率可视化软件系统.应用和钻探证实,三维电阻率视图在揭示采区水害分布、确定积水区和导水通道等方面,取得了良好的地质效果.     

不同地质统计方法在确定渗透系数场中的对比研究

应用地质统计方法研究渗透系数场的空间变异性。利用MMR含水层场地实测数据,通过去类分析、特异值处理、正态变换,逐步逼近研究区渗透系数的稳健变差函数,得到三维渗透系数场的几何各向异性套合模型。在此基础上,采用普通克里格法和指示克里格法、高斯序列模拟法和指示序列模拟法分别对数据进行插值和条件模拟。最后结合具体的地质条件,对四种方法在渗透系数场生成中的应用进行对比分析和评价。     

高密度电阻率法探测DNAPLs污染的适宜性探讨

重非水相流体（DNAPLs）的电阻率一般较高,在地下介质中与周围的地下水形成明显的电性差别。利用这一特性,非侵入式的高密度电阻率法（ERT）在探测DNAPLs污染场址中污染范围展现出一定优势。然而在实际场址的调查应用中,ERT有时无法探测到显著的高阻异常,或探测的高阻异常是由于介质的非均质性产生,从而导致对DNAPLs污染分布的误判。针对该问题,文章首先基于有效介质电阻率模型,计算了不同DNAPLs饱和度情形下的地层总电阻率值,从理论上解释了当DNAPLs污染程度较低时不存在显著高阻异常的原因。进一步构建理想ERT正演模型,探讨ERT探测DNAPLs污染的适宜性。研究结果表明:介质的非均质性（例如高阻的砾石等）对ERT探测具有很强的干扰,可造成DNAPLs分布范围的误判,因此传统静态ERT探测DNAPLs污染仅适用于均质介质或者非均质介质中DNAPLs污染程度较高的条件。对于非均质介质DNAPLs污染程度较低时,需采用基于时间域的差分反演方法屏蔽介质非均质性的影响,从而准确判定DNAPLs分布。     

EH-4系统观测资料的非远区场校正研究

EH－4系统实测的卡尼亚视电阻率在非远区场时应作必要的校正。基于导出的EH－4系统的卡尼亚电阻率ρxy和ρyx的表达式，采用多项式分段逼近的方法，建立EH－4系统测深数据的全频域的视电阻率算法。根据建立的算法对均匀大地、二层大地和三层大地模型的EH－4系统测深数据进行非远区场校正的结果表明：均匀大地时，全频域的视电阻率与真电阻率吻合很好，D型、K型和Q型断面的视电阻率曲线与MT法的接近，G型、H型和A型的视电阻率曲线得到很大改善。     

基于不同地质统计方法的渗透系数场对污染物运移的影响

渗透系数场的空间变异性是影响污染物运移结果的决定因素,而地质统计方法是解决渗透系数空间变异性的主要技术手段。本文利用野外场地实测数据,采用普通克里格法和指示克里格法、顺序高斯模拟法和顺序指示模拟法四种地质统计方法,插值估测和模拟再现随机渗透系数场,进而对比研究四种渗透系数场对大尺度污染物运移的影响。研究结果表明,污染羽的质心位置（一阶矩）主要由渗透系数的平均值来决定;污染羽在空间上的展布范围（二阶矩）主要受渗透系数空间变异方差的影响;条件模拟克服了估计法的平滑效果,较好地再现真实曲线的波动性,渗透系数（ lnK）估计方差与污染羽空间二阶矩随着条件模拟次数的增加而减小,并且顺序指示模拟程度更加明显。     

电测深资料求解岩含水层渗透系数的可行性——以吴营水源地为例

本文从电流场与渗流场相似入手，从理论上进行探讨，用类比的方法剖析相关量的物理含义及其量纲的类型，把视电阻率与渗透系数作为相关量，讨论二者间的关系规律，提出用平均视电阻率来研究和计算渗透系数的方法，实践证明吴家营水源地灰岩含水层渗透系数（Ｋ）与平均视电阻率（ρｓｍ）之间具有一定的相关关系且符合理论模型，而且Ｋ随ρｓｍ的变化具有地段性特点，此乃水文地质条件分区不同富水性地段的特征显示，用回归方法模拟Ｋ     

渗透系数的空间变异性对污染物运移的影响研究

随机水文地质学方法,较传统的确定性方法而言,是解决非均质含水层中水流和溶质运移问题的一种更为合理的手段。据以往研究,假设渗透系数场遵循对数正态分布,利用直接傅立叶变换方法来生成渗透系数随机场。应用基于随机理论的蒙特卡罗方法,来研究渗透系数的空间变异性对污染物运移结果的影响。实例研究表明,污染物在含水层中运移过程中污染羽的展布范围(二阶矩)随着渗透系数空间变异方差的增大而扩大,而污染羽在空间上的质心位置(一阶矩)基本不受方差的影响,仅取决于渗透系数随机场的均值大小。另外还分析了污染羽在各点的浓度变化方差和变异系数分别随渗透系数变异方差的变化状况。     

Surface electrical resistivity tomography and hydrogeological characterization to constrain groundwater flow modeling in an agricultural field site near Ferrara (Italy)

A portion of an unconfined alluvial aquifer located in the Padana Plain (Italy) was characterized following an integrated hydro-geophysical approach. Initially an electrical resistivity tomography (ERT) survey was employed to localize the boundaries of a modest paleo-channel body and to design the installation of a groundwater monitoring network. Multilevel slug-tests were performed to estimate the aquifer’s saturated hydraulic conductivities. Determined permeability values together with electrical resistivity data were correlated. The correlation resulted in a site specific bi-logarithmic linear relationship. Based on this relationship, punctually determined hydraulic conductivities were spatially extended over the studied flow domain. Finally, continuously measured piezometric heads were used to calibrate a 3D flow model. Sensitivity analysis was performed to confirm the reliability of the reconstructed permeability field, as well as, to assess the minimum number of direct measurements needed to safely characterize the selected aquifer portion. The integration of the ERT survey results with the classical hydrogeological tests can be conveniently applied to constrain the permeability field during flow model calibration. Although the applicability of the determined relationship is site specific, the followed procedure is useful especially when there is a need to optimize the available resources and in case of small-scale pilot studies.     

Assessment of integrated electrical resistivity data on complex aquifer structures in NE Nuba Mountains – Sudan

Two geophysical techniques were integrated to map the groundwater aquifers on complex geological settings, in the crystalline basement terrain in northeast Nuba Mountains. The water flow is structurally controlled by the northwest–southeast extensional faults as one of several in-situ deformational patterns that are attributed to the collision of the Pan-African oceanic assemblage of the Nubian shield against the pre-Pan African continental crust to the west. The structural lineaments and drainage systems have been enhanced by the remote sensing technique. The geophysical techniques used are: vertical electrical soundings (VES) and electrical resistivity tomography (ERT), in addition to hydraulic conductivity measurements. These measurements were designed to overlap in order to improve the producibility of the geophysical data and to provide a better interpretation of the hydrogeological setting of the aquifer complex structure. Smooth and Block inversion schemes were attempted for the observed ERT data to study their reliability in mapping the different geometries in the complex subsurface. The VES data was conducted where ERT survey was not accessible, and inverted smoothly and merged with the ERT in the 3D resistivity grid. The hydraulic conductivity was measured for 42 water samples collected from the distributed dug wells in the study area; where extremely high saline zones were recorded and have been compared to the resistivity values in the 3D model.     

Application of Stochastic Fracture Network with Numerical Fluid Flow Simulations to Groundwater Flow Modeling in Fractured Rocks

ＩＮＴＲＯＤＵＣＴＩＯＮＧｒｏｕｎｄｗａｔｅｒｏｒｆｌｕｉｄｆｌｏｗｍｏｄｅｌｉｎｇｉｎｆｒａｃｔｕｒｅｄｒｏｃｋｓｉｓａｃｏｍｐｌｉｃａｔｅｄｔｈｅｏｒｅｔｉｃａｌａｎｄａｐｐｌｉｅｄｔｏｐｉｃ.Ｂｏｔｈｔｈｅｏｒｅｔｉｃａｌｌｙａｎｄｏｐｅｒａｔｉｏｎａｌｌｙ ,ｉｔｉｓｉｍｐｏｒｔａｎｔｉｎｍａｎｙｆｉｅｌｄｓｓｕｃｈａｓｇｅｏｌｏｇｉｃａｌａｎｄｈｙｄｒｏｇｅｏｌｏｇｉｃａｌｅｎｇｉｎｅｅｒｉｎｇ ,ｅｎｖｉｒｏｎｍｅｎｔａｌｅｎｇｉｎｅｅｒｉｎｇａｎｄｐｅｔｒｏｌｅｕｍｅｎｇｉｎｅｅｒｉｎ…     

基于EnKF融合地球物理数据刻画含水层非均质性

含水层非均质性的刻画是模拟地下水中污染物运移的关键。以渗透系数为研究对象,构建了综合集合卡尔曼滤波方法、有效电阻率模型与地下水运移模型的同化框架,通过融合地球物理观测数据与污染物浓度观测数据来推估渗透系数的空间分布。基于理想算例,验证了该同化框架刻画含水层非均质渗透系数场的有效性,并针对不同初始参数信息与观测类型对比了耦合与非耦合水文地球物理方法的适用性。研究结果表明:基于集合卡尔曼滤波方法同化多种类型的观测数据,可有效地推估非均质参数空间分布。当初始信息较准确时,耦合方法的参数推估精度更高;初始信息存在偏差时,非耦合方法有更好的同化效果。由于非耦合方法计算成本较低且对初始信息缺失时适用性更强,在实际应用中可先基于非耦合方法初步估计参数,再利用耦合方法进一步提高参数推估精度。融合多种类型观测数据可有效提高参数推估效果。     

含水介质各向异性对渗透系数空间变异性统计的影响

针对目前渗透系数实测样本假定满足各向同性的局限,以Borden含水层试验场实测数据为例,通过反证法进行相关数据分析,总结了含水介质各向异性对渗透系数空间变异性统计的影响,并指明渗透系数各向同性假设的适用条件以及不合理之处。同时相应地给出了研究尺度下渗透系数场能否采用平稳随机场描述的判定依据。最后为今后进一步开展渗透系数空间变异性研究提出几点建议。     

The hydraulic conductivity field and groundwater flow in the unconfined aquifer system of the Keta Strip,Ghana

This study investigates the hydraulic conductivity field and the groundwater flow pattern as predicted by a calibrated steady state groundwater flow model for the Keta Strip, southeastern Ghana. The hydraulic conductivity field is an important parameter in evaluating aquifer properties in space, and in general basin-wide groundwater resources evaluation and management. This study finds that the general hydraulic conductivity of the unconsolidated unconfined aquifer system of the Keta Strip ranges between 2 m/d and 20 m/d, with an average of 15 m/d. The spatial variation in horizontal hydraulic conductivity appears to take the trend in the variations in the nature of the material in space. Calibrated groundwater recharge suggests that 6.9–34% of annual precipitation recharges the shallow aquifer system. This amount of recharge is significant and suggests high fortunes in terms of groundwater resources development for agriculture and industrial activities in the area. A spatial distribution of groundwater recharge from precipitation is presented in this study. The spatial pattern appears to take the form of the distribution in horizontal hydraulic conductivity, and suggests that the vertical hydraulic conductivity takes the same pattern of spatial variation as the horizontal hydraulic conductivity. This is consistent with observations in other areas. The resulting groundwater flow is dominated by local flow systems as the unconfined system is quite shallow. A general northeast – southwest flow pattern has been observed in the study area.     

A coupled hydromechanical system defined through rock mass classification schemes

Changes in the hydraulic conductivity field, resulting from the redistribution of stresses in fractured rock masses, are difficult to characterize due to complex nature of the coupled hydromechanical processes. A methodology is developed to predict the distributed hydraulic conductivity field based on the original undisturbed parameters of hydraulic conductivity, Rock Mass Rating (RMR), Rock Quality Designation (RQD), and additionally the induced strains. The most obvious advantage of the methodology is that these required parameters are minimal and are readily available in practice. The incorporation of RMR and RQD, both of which have been applied to design in rock engineering for decades, enables the stress-dependent hydraulic conductivity field to be represented for a whole spectrum of rock masses. Knowledge of the RQD, together with the original hydraulic conductivity, is applied to determine the effective porosity for the fractured media. When RQD approaches zero, the rock mass is highly fractured, and fracture permeability will be relatively high. When RQD approaches 100, the degree of fracturing is minimal, and secondary porosity and secondary permeability will be low. These values bound the possible ranges in hydraulic behaviour of the secondary porosity within the system. RMR may also be applied to determine the scale effect of elastic modulus. As RMR approaches 100, the ‘softening’ effect of fractures is a minimum and results in the smallest strain-induced change in the hydraulic conductivity because the induced strain is uniformly distributed between fractures and matrix. When RMR approaches zero, the laboratory modulus must be reduced significantly in order to represent the rock mass. This results in the largest possible change in the hydraulic conductivity because the induced strain is applied entirely to the fracture system. These values of RMR bound the possible ranges in mechanical behaviour of the system. The mechanical system is coupled with the hydraulic system by two empirical parameters, RQD and RMR. The methodology has been applied to a circular underground excavation and to qualitatively explain the in situ experimental results of the macropermeability test in the drift at Stripa. Copyright © 1999 John Wiley & Sons, Ltd.     

Stochastic groundwater flow analysis in the presence of trends in heterogeneous hydraulic conductivity fields

Due to changes in lithostatic pressure, differential fracturing across bedding planes and irregularities in depositional environments, hydraulic conductivity exhibits heterogeneities and trends at various spatial scales. Using spectral theory, we have examined the effect of trends in hydraulic conductivity on (1) the solution of the mean equation for hydraulic head, (2) the covariance of hydraulic head, (3) the cross-covariances of hydraulic head and log-hydraulic conductivity perturbations and their gradients, and (4) the effective hydraulic conductivity. It is shown that the field of hydraulic head is sensitive to the presence of trends in ways that cannot be predicted by the classical analysis based on stationary hydraulic conductivity fields. The controlling variables for the second moments of hydraulic head are the mean hydraulic gradient, the correlation scale of log-hydraulic conductivity and its variance, and the slope of the trend in log-hydraulic conductivity. The mean hydraulic gradient introduces complications in the analysis since it is, in general, spatially variable. In this respect, our results are approximate, yet indicative of the true role of spatially variable patterns of log-hydraulic conductivity on groundwater flow systems.