新安江右坝裂隙岩体渗流同位素示踪研究

应用同位素示踪方法对新安江右坝的裂隙岩体渗流场进行了研究。在单孔中应用同位素示踪技术可得到流速、流向、垂向流等水力学参数,结合相关的地下水井流理论,就可计算出裂隙等效水力隙宽、渗透系数等反映裂隙渗透性的参数,对其进行分析,就确定了渗流场的分布以及集中渗漏通道的位置。     

基于流量衰减曲线的岩溶含水层水文地质参数推求方法

分析流量衰减曲线是研究岩溶含水层水文地质参数的重要手段。在阐明流量衰减分析理论的基础上,通过连续监测一年寨底地下河出口的小时流量过程,共获取8 784组流量数据。首先选取典型的2013年9月25日3:00暴雨后至2013年10月13日17:00的446组流量数据作为研究对象,并利用MATLAB软件对该时间段的流量过程进行曲线拟合,从而确定了不同阶段流量衰减系数分别为0.642 5、0.025 8、0.001 5;然后通过流量衰减系数计算寨底不同含水介质导水系数T与给水度S的比值;最后利用裂隙率近似给水度的方法确定了寨底不同含水介质的水文地质参数。结果表明,寨底岩溶管道、裂隙、基岩介质导水系数分别为:258.87、3 865.98、40.64 m²/d。说明辅以必要的水文地质调查,确定管道分布和径流距离,流量衰减曲线可用于计算岩溶含水层不同介质的导水性能和给水度。     

Reflection and refraction of SH‐waves at a corrugated interface between two monoclinic elastic half‐spaces

A problem of reflection and transmission of shear wave incident upon a corrugated interface between two monoclinic solid half‐spaces have been investigated. Rayleigh's method of approximation is used to investigate the reflection and transmission coefficients for first and second approximation of the corrugation. For a special interface, closed‐form expressions of these coefficients for the first order approximation of the corrugation are obtained. It is found that these coefficients are functions of amplitude of corrugation, elastic parameters of the media, frequency of the incident wave and angle of incidence. The numerical computations are performed for a specific model and the results obtained are presented graphically. It is found that the reflection and transmission coefficients are strongly influenced by the corrugation and the elastic properties of the media. Results of some earlier workers in this field have been reduced as particular case from the present formulation. Copyright © 2004 John Wiley & Sons, Ltd.     

Implementation of the finite element method in the three‐dimensional discontinuous deformation analysis (3D‐DDA)

A modified three‐dimensional discontinuous deformation analysis (3D‐DDA) method is derived using four‐noded tetrahedral elements to improve the accuracy of current 3D‐DDA algorithm in practical applications. The analysis program for the modified 3D‐DDA method is developed in a C++ environment and its accuracy is illustrated through comparisons with several analytical solutions that are available for selected problems. The predicted solutions for these problems using the modified 3D‐DDA approach all show satisfactory agreement with the corresponding analytical results. Results presented in this paper demonstrate that the modified 3D‐DDA method with discontinuous modeling capabilities offers a useful computational tool to determine stresses and deformations in practical problems involving fissured elastic media with reasonable accuracy. Copyright © 2008 John Wiley & Sons, Ltd.     

裂隙性粘土边坡渐进性破坏的FLAC模拟

裂隙性粘土是一种具有特殊性质的粘土,这种粘土边坡的破坏是渐进性的,传统的极限平衡法无法描述分析它的破坏过程和机理,也无法精确地判断它的长期稳定性。用FLAC程序对裂隙性粘土边坡的渐进性破坏过程进行了描述,算例表明在一定程度上用FLAC程序来模拟边坡渐进破坏和进行稳定性判断是可行的。     

隧道开挖的双重介质水力耦合模型有限元分析

基于混合物理论孔隙-裂隙岩体的双重孔隙介质水力耦合计算的微分方程,利用伽辽金有限元法提出了相应的有限元公式,并基于岩体分类指标（RQD,RMR）提出了与岩体应变状态相关的渗透系数计算公式。编制了裂隙岩体双重介质流固耦合的2-D有限元程序,给出的验证算例表明,该程序是合理和实用的。同时将该程序用于隧道开挖的模拟计算,探讨渗流效应对开挖隧道围岩变形与渗流场的影响。计算结果表明,在隧道设计中不考虑渗流的影响是偏于不安全的。     

改进的遗传算法及其在渗流参数反演中的应用

利用水头实测资料,以裂隙组的渗透系数比例因子为待反演的参数向量,在采用基本遗传算法进行参数反演研究的基础上,针对裂隙岩体无压渗流参数反问题计算量过大,目标参数众多以及参数可能变化范围大等特点,提出了一种混合遗传算法求解此类问题,力求克服简单遗传算法在解决此类问题时存在的局部搜索能力弱、易出现早熟收敛及计算量大等缺陷,并通过典型岩坡渗流算例进行验证,同时给出了基本遗传算法、传统单纯形算法的反演成果。计算结果表明,该方法保持了基本遗传算法优点,并有效地提高了算法的运行效率,从而为求解裂隙岩体无压渗流参数反问题等计算量大的系列问题提供了新的途径。     

应力腐蚀和压力溶解引起的裂隙开度-刚度变化对THMM耦合过程影响的数值模拟

在使用Yasuhara等建立的裂隙开度的应力腐蚀和压力溶解模型的基础上,考虑裂隙闭合量对裂隙刚度的影响,针对一个假设的位于非饱和双重孔隙-裂隙岩体中且有核素泄漏的高放废物地质处置模型,拟定两种计算工况：（1）裂隙的刚度系数是裂隙闭合量的线性函数;（2）裂隙刚度为常数,进行热-水-应力-迁移耦合的二维有限元数值模拟,考察了岩体中的温度、裂隙开度的闭合速率、闭合量、孔（裂）隙水压力、地下水流速、核素浓度、裂隙刚度和正应力的变化、分布等情况。结果表明,两种工况岩体中的温度场、孔（裂）隙水压力、地下水流速、核素浓度无明显差别;裂隙闭合基本由应力腐蚀产生;在相同计算时间内两种工况的裂隙闭合量较为接近,工况1略大;工况1中离玻璃固化体越近,裂隙刚度值越高,并且在玻璃固化体附近的应力值较大,且集中程度较高。     

裂隙岩体流固耦合双重介质模型的有限元计算

基于孔隙-裂隙岩体的双重孔隙介质流固耦合计算的微分方程,利用伽辽金有限元法提出的相应有限元公式,并基于岩体分类指标（RQD,RMR）提出了与岩体应力状态相关的渗透系数计算公式。编制了相应的有限元程序并给出了应用算例,将计算结果与相关文献作了比较,得出相关结论。     

Numerical modelling of transient contaminant migration problems in infinite porous fractured media using finite/infinite element technique. Part I: Theory

The leakage effect in porous fissured media has been considered in a general sense by introducing a new expression of the leakage term in this paper. The double porosity concept is employed and the related expressions are formulated using the upwind finite element approach. Considering the infinite extension of the problem domain, a mapped transient infinite element has been presented to simulate the far field of the infinite medium. Since the mass transfer function of the present mapped transient infinite element is dependent on both space and time variables, the mechanism of transient contaminant migration problems in infinite porous fractured media can be rigorously simulated because the property matrices of the element are evaluated at any time of interest. By comparing the current numerical results with the analytical ones, the accuracy, correctness and effectiveness of the present method have been established. Three different time discretization schemes were examined and it was found that either the central difference or the backward difference approximation is suitable for the upwind finite element simulation of transient contaminant migration problems.     

Spreading of a NAPL lens in a double-porosity medium

Spreading of a non-aqueous phase liquid (NAPL) denser than water (DNAPL) lens (mound) in the unsaturated zone of double-porosity aquifer above an impervious plane boundary is investigated. The double-porosity aquifer is conceptualized as a fracture network surrounding pervious blocks. Vertical gravity equilibrium is assumed to prevail in each one of the two media, fractures and blocks. Through vertical integration, two coupled partial differential equations for the DNAPL content in each medium are obtained. The mass exchange rate between high- and low-permeability media is considered as a function of NAPL content. The dominant effect is gravity, whereas capillary forces are negligible. Analytical solutions for one-dimensional and axisymmetric problems of mound spreading are obtained.     

A naturally efficient numerical technique for porous convection stability with non‐trivial boundary conditions

A highly efficient numerical technique is presented for solving eigenvalue problems which arise in complicated convection — instability studies in porous media. The differential equations are written as a system of ‘natural’ variables which are suggested by the way the boundary conditions arise. The method easily gives high resolution in boundary layers, yields all the eigenvalues and eigenfunctions, deals with complex coefficients, and can handle spatially dependent coefficients in a very efficient manner. The numerical technique is motivated by the practical problem of salinization in porous sands in arid zones as is beautifully modelled by Gilman and Bear. Since the salinization study of Gilman and Bear is a prototype for the field of convective motion in unsaturated porous soils and this field is one which is increasingly occupying geotechnical attention, we believe the numerical method presented here has much potential. Copyright © 2000 John Wiley & Sons, Ltd.     

An approach for the geophysical assessment of fissuring of estuary and river flood embankments: validation against two case studies in England and Scotland

This paper proposes a two-stage geophysical approach to map the vertical cracking and the structural integrity of flood embankments made up of clay geomaterials susceptible to fissuring. The first stage is based on a ‘coarse-resolution’ investigation using conventional electrical resistivity tomography (ERT) equipment to identify the fissured zones in the embankment. This step is complemented by an additional geophysical technique, electromagnetic, to verify the ERT measurements. The second stage is based on a ‘high-resolution’ investigation using a miniature ERT system previously developed at the laboratory scale for detailed mapping of the fissure patterns. The ‘coarse-resolution’ stage is the major focus of this paper and was validated against two case studies in England and Scotland. Longitudinal ERT survey provided a tomographic picture of the upper desiccated zones of the embankments and fissured areas in 2-D, validating the range of resistivity results obtained previously on a fissured clay model in the laboratory. A transversal embankment resistivity tomography was also completed to show the positions of fissured zones in detail in the field. The electromagnetic technique as a fast screening tool allowed cross checking the ERT results and was also efficient in detecting high and low conductivity zones, indicating areas of potential weakness during flash floods and heavy rain. The southern embankment in England showed more fluctuations in the conductivity and resistivity than the north embankment in Scotland, likely to be due to the differences in climate, vegetation and location characteristics between the two sites. Conclusions were also drawn on the potential weaknesses for both embankments and the effect of vegetation on conductivity measurements.     

受振动荷载扰动裂隙性黄土单轴压缩力学行为研究

通过振动及单轴压缩试验,研究了受振动荷载扰动裂隙性黄土的单轴压缩力学行为。结果表明：裂隙性黄土的单轴压缩破坏模式表现为压裂破坏、滑移破坏、滑移-压裂复合破坏以及压剪破坏4种类型;振动扰动对单轴压缩条件下裂隙性黄土的破坏模式无显著影响,其破坏模式主要由初始裂隙的倾角控制。振动幅值和频率对裂隙性黄土应力-应变曲线的类型及特征无显著影响,不同振动参数条件下试样的应力-应变曲线均表现为应变软化型,且45°倾角试样的应力-应变曲线呈现出第二峰值强度高于第一峰值强度的“双峰”变化特征。单轴抗压强度随振动幅值和频率的增大均呈现出近似线性减小的变化规律;不同振动参数条件下试样单轴抗压强度随裂隙倾角增大近似呈现出“双V”变化特征。构建了受振动荷载扰动裂隙性黄土的二元介质本构模型,可较好预测其单轴压缩过程的应力-应变关系及单轴抗压强度。     

FEM validation of a double porosity elastic model for consolidation of structurally complex clayey soils

Laboratory consolidation of structured clayey soils is analysed in this paper. The research is carried out by two different methods. The first one treats the soil as an isotropic homogeneous equivalent Double Porosity (DP) medium. The second method rests on the extensive application of the Finite Element Method (FEM) to combinations of different soils, composing 2D or fully 3D ordered structured media that schematically discretize the complex material. Two reference problems, representing typical situations of 1D laboratory consolidation of structured soils, are considered. For each problem, solution is obtained through integration of the equations governing the consolidation of the DP medium as well as via FEM applied to the ordered schemes composed of different materials. The presence of conventional experimental devices to ensure the drainage of the sample is taken into account through appropriate boundary conditions. Comparison of FEM results with theoretical results clearly points out the ability of the DP model to represent consolidation processes of structurally complex soils. Limits of applicability of the DP model may arise when the rate of fluid exchange between the two porous systems is represented through oversimplified relations. Results of computations, obtained having assigned reasonable values to the meso‐structural and to the experimental apparatus parameters, point out that a partially efficient drainage apparatus strongly influences the distribution along the sample and the time evolution of the interstitial water pressure acting in both systems of pores. Data of consolidation tests in a Rowe's cell on samples of artificially fissured clays reported in the literature are compared with the analytical and numerical results showing a significant agreement. Copyright © 2000 John Wiley & Sons, Ltd.     

复杂介质地震波场研究进展

从传统的层状介质研究发展到对非均匀各向异性介质乃至随机介质研究,是解决当今高难度油气勘探和高标准油气开发问题的需要。面对这种发展趋势,近年来国内外有关学者进行了大量的复杂介质地震波场研究。本文对其在裂缝介质、双相介质及随机介质等方面的研究所取得的成果作一综述,以期为今后在该领域的研究积累较为系统的基础资料。     

多尺度有限单元法求解非均质多孔介质中的三维地下水流问题

应用多尺度有限单元法模拟非均质多孔介质中的三维地下水流问题。与传统有限单元法相比，多尺度有限单元法的基函数具有能反映单元内参数变化的优点，所以这种方法能在大尺度上抓住解的小尺度特征获得较精确的解。在介绍多尺度有限单元法求解非均质多孔介质中三维地下水流问题的基本原理之后，对参数水平方向渐变垂直方向突变的非均质多孔介质中的三维地下水流和Borden实验场的三维地下水流分别用多尺度有限单元法和传统等参有限单元法进行了计算，结果表明在模拟高度非均质多孔介质中的三维地下水流问题时，多尺度有限单元法比传统有限单元法有效，既节省计算量又有较高的精度；在模拟非均质性弱的多孔介质中的三维地下水流问题时，多尺度有限单元法虽然也能在大尺度上获得较为精确的解，但效果不明显。     

B‐bar based algorithm applied to meshfree numerical schemes to solve unconfined seepage problems through porous media

The object of this work is to establish a meshfree framework for solving coupled, steady and transient problems for unconfined seepage through porous media. The Biot's equations are formulated in displacements (or u − w) assuming an elastic solid skeleton. The free surface location and its evolution in time are obtained by interpolation of pore water pressures throughout the domain. Shape functions based on the principle of local maximum entropy are chosen for the meshfree approximation schemes. In order to avoid the locking involved in the fluid phase of the porous media, a B‐bar based algorithm is devised to compute the average volumetric strain in a patch composed of various integration points. The efficiency of such an implementation for one phase problems is shown through the Benchmark problem, Cook's membrane loaded by a distributive shear load. The proposed methodology is firstly applied to various classical examples in unconfined steady seepage problems through earth dams, then to the dynamic consolidation of a soil column. The results obtained for both problems are quite satisfactory and demonstrate the feasibility of the proposed method in solving coupled problems in porous media. Copyright © 2015 John Wiley & Sons, Ltd.     

A fully coupled element‐free Galerkin model for hydro‐mechanical analysis of advancement of fluid‐driven fractures in porous media

Hydraulic fracturing (HF) of underground formations has widely been used in different fields of engineering. Despite the technological advances in techniques of in situ HF, the industry uses semi‐analytical tools to design HF treatment. This is due to the complex interaction among various mechanisms involved in this process, so that for thorough simulations of HF operations a fully coupled numerical model is required. In this study, using element‐free Galerkin (EFG) mesh‐less method, a new formulation for numerical modeling of hydraulic fracture propagation in porous media is developed. This numerical approach, which is based on the simultaneous solution of equilibrium and continuity equations, considers the hydro‐mechanical coupling between the crack and its surrounding porous medium. Therefore, the developed EFG model is capable of simulating fluid leak‐off and fluid lag phenomena. To create the discrete equation system, the Galerkin technique is applied, and the essential boundary conditions are imposed via penalty method. Then, the resultant constrained integral equations are discretized in space using EFG shape functions. For temporal discretization, a fully implicit scheme is employed. The final set of algebraic equations that forms a non‐linear equation system is solved using the direct iterative procedure. Modeling of cracks is performed on the basis of linear elastic fracture mechanics, and for this purpose, the so‐called diffraction method is employed. For verification of the model, a number of problems are solved. According to the obtained results, the developed EFG computer program can successfully be applied for simulating the complex process of hydraulic fracture propagation in porous media. Copyright © 2016 John Wiley & Sons, Ltd.     

断裂多孔岩石中的复杂反应流体动力学

:断裂多孔介质中的流体流动和组分输运过程对于理解金属成矿过程、油气运移、污染物质的输运和分布十分重要 ,近 2 0年来得到了广泛的研究。文中综述了这一领域的进展 ,涉及断裂多孔介质中的流体流动和组分输运的连续介质模型、多孔介质和断裂多孔介质的渗透率关系 ,特别是重点介绍了渗滤理论及其在断裂多孔介质中的流体流动和组分输运中的应用