裂隙岩体渗透系数确定方法研究

裂隙岩体渗透系数以及渗透主方向的确定对研究岩体渗透性大小及各向异性具有重要意义。高放废物地质处置库介质岩体的渗透性能将直接影响其使用安全性。本文运用离散裂隙网络模拟的方法对我国高放废物处置库甘肃北山预选区3#钻孔附近裂隙岩体进行了渗透性质分析。通过对3#钻孔1715～1780m段压水试验数据的反演,标定了离散裂隙网络渗流模型中的裂隙渗透参数(导水系数T)。利用标定的离散裂隙网络模型对场区裂隙岩体进行了渗流模拟,确定了该区域裂隙岩体的渗流表征单元体(REV)的尺寸大小以及渗透主值和主渗透方向。运用离散裂隙网络模型计算得出的渗透主值的几何均值与现场压水试验计算结果较接近,证明了计算结果的有效性。     

金川Ⅲ矿区岩体分类及其力学特性参数预测

本文针对金川Ⅲ矿区岩体质量及参数,并在11个钻孔数据资料的基础上,采用数学统计和数值模型来分别分析矿区岩体质量和力学参数的特征。对11个钻孔的资料进行详细的统计分析,进行每一个钻孔不同深度的岩体节理裂隙统计和岩体分类评价,并由此获得了建立矿块力学参数模型的数据。再建立表征整个岩体特性的三维力学模型,并揭示了矿块力学特性的分布特性和变化规律。此研究为崩落法采矿设计和数值模拟提供可靠的依据。     

单轴压缩条件下溶蚀岩体力学特性数值试验研究

工程尺度的溶蚀岩体难以开展室内外力学试验,导致溶蚀岩体的力学参数获取存在困难,因此基于岩石矿物含量特征、室内力学试验和岩体结构面特征,采用三维颗粒流离散元法,通过建立等效孔隙型溶蚀岩体模型进行单轴压缩数值试验,进而分析其力学特性和变形破坏机制。研究表明：采用改进平行黏结模型,基于细观矿物特征,通过元胞自动机算法剔除颗粒可建立孔隙型溶蚀岩体模型;加载前期,岩体结构面首先快速破坏而产生以剪切为主的微裂纹;随加载进行,岩块内逐渐产生以拉裂纹为主的破坏,其微裂纹呈指数增加,而结构面微裂纹先激增后趋于稳定;相同轴向应变时,岩块内拉裂纹随溶蚀率增加而增加,而结构面剪切微裂纹减少;岩体的变形破坏分为结构面快速破坏、岩块弹性变形、岩块塑性变形和岩体完全破坏等4个阶段,其破坏形式随溶蚀率增加而从整体均匀性破坏转化为局部结构性破坏;溶蚀使岩体强度降低,溶蚀率与单轴抗压强度和变形模量分别呈反比和负指数函数关系。     

裂隙岩体局部化破坏分析的多重势面分岔模型

采用多重势面弹塑性分岔理论对裂隙岩体的局部化破坏进行计算模拟．基于多重势面弹塑性理论分析局部化问题，构造了适用于裂隙岩体破坏的多重势面不连续分岔模型，并使用数值方法求解局部化方阶在有限元方法的基础上，使用该模型计算裂隙岩体的局部化破坏条带．算例分析表明这一模型用于分析裂隙岩体的局部化破坏是有效的．     

裂隙岩体渗流模型研究现状与展望

裂隙岩体渗流对于边坡、地下工程及基础岩土体的承载能力有显著的制约作用。本文简要地介绍了多种裂隙岩体渗流模型研究现状,评述了几类比较有代表性的渗流模型特点以及存在的不足,为选取合理的数学模型用于求解具体的裂隙岩体渗流问题提供了参考依据,并在上述基础上提出了一些需要进一步研究的问题。指出就目前应用最为广泛的等效连续介质模型而言,裂隙岩体几何参数、有效孔隙度以及等效渗透张量的确定仍有待于更深入的研究,而从细观力学结构入手研究渗流耦合模型将具有重大的理论意义和实用价值。     

含单裂隙岩体蠕变破裂裂纹扩展演化规律研究

岩体内部存在大量不同尺寸裂隙，形成不同结构形式的裂隙岩体，裂隙的存在会降低岩体的力学性能，对于岩体工程的整体稳定性具有很大的影响。本研究基于滑动翼裂纹模型对恒定载荷下含单裂隙岩体蠕变变形进行了理论推导和分析，基于MAPLE软件直观给出恒定载荷下裂隙角度、长度对蠕变变形的影响，采用连续-非连续数值分析软件GDEM,对预制单裂隙在不同倾角、不同长度情况下岩体的蠕变破裂裂纹扩展进行数值模拟。通过理论分析和数值模拟可以看出，岩体蠕变速率随着裂隙长度的增大而增大；随着裂隙角度的增加出现先增大后减小的趋势，其中裂隙倾角为30°时蠕变速率最大。含单裂隙岩体蠕变破裂裂纹扩展规律的研究可为岩体工程蠕变破裂规律研究及蠕变治理提供一定的理论依据和参考。     

裂隙岩体中非饱和渗流与运移的概念模型及数值模拟

探讨了裂隙岩体中非饱和地下水渗流与溶质运移的几种概念模型的构造及数值模拟问题 ,如裂隙网络模型、连续体模型、等效连续体模型、双孔隙度 (单渗透率 )模型、双渗透率模型、多组份连续体模型等。在裂隙岩体中 ,非饱和地下水的渗流可能只局限于岩体中的岩石组份、或裂隙网络 ,也可能在裂隙和岩石中同时发生 ;对前一种情形只需考虑单一连续体中的流动 ,而后一种情况则需要包括地下水在岩石和裂隙之间的交换。岩体中的裂隙网络往往是溶质运移的主要通道 ;但当溶质在裂隙与岩石之间的渗透和扩散是重要的运移机制时 ,就需要考虑岩石与裂隙界面处的溶质交换。为了模拟岩石与裂隙之间地下水和溶质的交换 ,就需要了解岩石与裂隙之间相互作用的模式和范围 ,使得这类问题的概念模型较单一连续体模型多了一层不确定性、其数值模拟也变得更为困难。因为在实际问题中不易、甚至根本不能判别非饱和渗流的实际形态 ,具体采用哪种模型主要取决于分析的目的和对现场数据的掌握程度。不论哪种模型都会受到模型及参数不确定性的影响 ,因此必须考虑与其他辅助模型的比较.     

单裂隙岩体水力裂缝扩展机理的数值模拟

应用真实破裂过程分析系统软件RFPA2D2.0-Flow,采用细观方法表达岩石材料非均匀特性,建立了渗流-应力-损伤耦合模型;模拟分析了裂隙岩体在具有不同角度、不同长度、天然裂隙中点与井孔之间不同间距的单一天然裂隙条件下的水压致裂过程;通过对模型内拉伸裂缝萌生及扩展的描述,得出了天然裂隙角度、长度、天然裂隙中点与井孔的间距对裂隙岩体内水力裂纹的扩展及与天然裂隙相交模式的影响趋势。模拟分析结果表明:天然裂隙角度对模型稳定性影响较小,但对裂纹扩展模式的影响较大;天然裂隙与井孔间距的减小使模型的失稳水压值随着裂隙角度的增加而逐渐增大,但最大值与间距未减小时基本相同;随着天然裂隙长度的增加,模型的承载能力逐渐减弱,且减弱的幅度逐渐增大。     

裂隙岩体损伤局部化破坏分岔模型及其应用

采用概率统计方法分析节理裂隙岩体的几何特征，定义了反映裂隙岩体几何特征的组构张量．基于不可逆热力学理论，通过裂纹扩展的细观分析，得出了损伤的发展机理和演化方程，把损伤演化和裂隙的几何特征的变化联系起来，建立了弹塑性损伤本构关系．为分析含有节理裂隙岩体在发生局部化破坏时的特征，通过对发生局部化时的裂隙岩体的分析。构造了适用于节理裂隙岩体局部化分析的不连续分岔模型．利用非线性规划数值解法，可以得出局部化破坏的方向特征．在有限元方法中，根据该模型给出了节理裂隙岩体相关的算例，分析表明该模型用于分析裂隙岩体的局部化破坏是有效的．     

黄河黑山峡大柳树坝坝址岩体渗漏量及灌浆量数值模拟研究

大坝坝址岩体的渗漏研究对于大坝的安全运营至关重要。在进行渗漏治理的措施中,灌浆是常见的一种,但目前确定灌浆的各项参数,如配合比,孔距等,需通过灌浆试验花费大量的人财物力及时间才能确定。为快速查清每一种灌浆方案所灌浆量,本文针对模拟坝址水库渗漏灌浆效果的三维数值模拟开展了计算研究。在对大柳树坝坝址岩体地质情况及渗漏情况分析的基础上,将蓄水前的地下水运移规律模拟出,与实际通过钻孔平硐资料及分析得出的实际地下水运移规律做比较,验证对地下水运移规律的假设; 然后将蓄水后的地下水运移规律模拟出,确定地下水的运移规律,了解地下水运移的方向,验证与根据实际情况初步判断的地下水运移方向是否吻合; 最后通过数值模拟确定在一定灌浆浆液配合比的情况下所需灌浆量与实际的灌浆试验结果做比较。得出结果与实际灌浆试验结果吻合,证明通过数值模拟,建立浆液配合比与灌浆效果之间的关系,在类似情况快速地选取配合比是可行的。     

Stochastic Analysis of Macrodispersion in a Semi‐Confined Aquifer

In this paper, the macroscopic dispersion resulting from one and twodimensional flows through a semiconfined aquifer with spatially variable hydraulic conductivity K which is represented by a stationary (statistically homogeneous) random process is analyzed using the spectral representation technique. Stochastic fluctuation equations of the steady flow and solute transport are solved to construct the macroscopic dispersive flux and evaluate the resulting macrodispersivity tensor in terms of the leakage factor and input covariances describing the hydraulic conductivity in a semiconfined aquifer bounded by a leaky layer above and an impervious stratum below. The macrodispersivity tensor is studied using some convenient forms of the log hydraulic conductivity process. The sensitivity of the resulting macrodispersivity to the input covariances is discussed along with the influence of the leakage factor for both one and twodimensional flows. It is found that the longitudinal macrodispersivities are increased due to the presence of leakage, while the transverse macrodispersivities are reduced due to leakage.     

液压减振器随机热力学模型与油温变化研究

为了获得具有更多信息和更加接近工程实际的液压减振器油温热力学模型,将随机不确定性理论引入到液压减振器油温传统热力学模型中进行研究。将液压油密度、导热系数、比热容和运动粘度作为随机变量,运用求解函数数字特征的代数综合法建立减振器随机热力学模型,进而获得油液传热过程规律。将油温随机热力学模型研究结果和传统模型的计算结果与实验结果进行比较,证明随机不确定性理论的引入可行且随机热力学模型比传统模型更加优越。     

Retrieval of Spacewise Dependent Hydraulic Properties of Anisotropic Rocks from Transient Flow Experiments

An inverse finite-difference method (FDM) is developed to characterise the spatially dependent components of the hydraulic conductivity tensor together with the specific storage for anisotropic materials, using experimental data generated from the direct FDM. This simulated surface data serves as additional information to a genetic algorithm (GA) optimisation procedure, using a modified least squares functional, that minimises the difference between the experimental data and the FDM-predicted boundary pressure and/or average hydraulic flux measurements based on current hydraulic conductivity tensor and specific storage estimates.     

渗流问题单参数摄动随机反演方法

在正演随机模拟方法的基础上,结合Taylor展开式和随机变量的摄动方程,讨论了随机参数反演问题,提出了摄动随机反演方法,给出了一阶均值反演准则和二阶均值反演准则,提出了单随机变量的均值和方差的表示方法,计算均值时采用改进的遗传算法,计算方差时采用统计的方法。给出了一个Thies模型反演导水系数的例子,计算表明该方法简单实用,效果良好。     

THEORETICAL MODEL OF EFFECTIVE STRESS COEFFICIENT FOR ROCK/SOIL-LIKE POROUS MATERIALS

Physical mechanisms and influencing factors on the effective stress coefficient for rock/soil-like porous materials are investigated, based on which equivalent connectivity index is proposed. The equivalent connectivity index, relying on the meso-scale structure of porous material and the property of liquid, denotes the connectivity of pores in Representative Element Area （REA）. If the conductivity of the porous material is anisotropic, the equivalent connectivity index is a second order tensor. Based on the basic theories of continuous mechanics and tensor analysis, relationship between area porosity and volumetric porosity of porous materials is deduced. Then a generalized expression, describing the relation between effective stress coefficient tensor and equivalent connectivity tensor of pores, is proposed, and the expression can be applied to isotropic media and also to anisotropic materials. Furthermore, evolution of porosity and equivalent connectivity index of the pore are studied in the strain space, and the method to determine the corresponding functions in expressions above is proposed using genetic algorithm and genetic programming. Two applications show that the results obtained by the method in this paper perfectly agree with the test data. This paper provides an important theoretical support to the coupled hydro-mechanical research.     

Effect of 2D Image Resolution on 3D Stochastic Reconstruction and Developing Petrophysical Trend

Multi-resolution digital rock physics (DRP) makes it possible to up-scale petrophysical properties from micron size to core sample size using two-dimensional (2D) thin section images. Resolution of 3D images and sample size are challenging problems in DRP where high-resolution images are acquired from small samples using inefficient and expensive micro-CT facilities. Three-dimensional stochastic reconstruction is an alternative approach to overcome these challenges. In this paper, we use multi-resolution images and investigate effect of 2D image resolution on 3D stochastic reconstruction and development of petrophysical trends for our two sandstone and carbonate original representative volume elements (RVEs). The proposed method includes three steps. In the first step, the spatial resolution of our original RVEs is decreased synthetically. In the second step, stochastic RVEs are realized for each resolution using two perpendicular images, correlation functions, and phase recovery algorithm. In the reconstruction method, a full set of two-point correlation functions (TPCFs) is extracted from two perpendicular 2D images. Then TPCF vectors are decomposed and averaged to realize 3D stochastic RVEs. In the third step, petrophysical properties like relative and absolute permeability as well as porosity and formation factor are computed. The output is used to develop trends for petrophysical properties in different resolutions. Experimental results illustrate that the proposed method can be used to predict petrophysical properties and reconstruct 3D RVEs for resolutions unavailable in the acquired 2D or 3D data.     

Membrane Curvatures and Stress-strain Full Fields of Axisymmetric Bulge Tests from 3D-DIC Measurements. Theory and Validation on Virtual and Experimental results

The bulge test is mostly used to analyze equibiaxial tensile stress state at the pole of inflated isotropic membranes. Three-dimensional digital image correlation (3D-DIC) technique allows the determination of three-dimensional surface displacements and strain fields. In this paper, a method is proposed to determine also the membrane stress tensor fields for in-plane isotropic materials, independently of any constitutive equation. Stress-strain state is then known at any surface point which enriches greatly experimental data deduced from the axisymmetric bulge tests. Our method consists, first in calculating from the 3D-DIC experimental data the membrane curvature tensor at each surface point of the bulge specimen. Then, curvature tensor fields are used to investigate axisymmetry of the test. Finally in the axisymmetric case, membrane stress tensor fields are determined from meridional and circumferential curvatures combined with the measurement of the inflating pressure. Our method is first validated for virtual 3D-DIC data, obtained by numerical simulation of a bulge test using a hyperelastic material model. Afterward, the method is applied to an experimental bulge test performed using as material a silicone elastomer. The stress-strain fields which are obtained using the proposed method are compared with results of the finite element simulation of this overall bulge test using a neo-Hookean model fitted on uniaxial and equibiaxial tensile tests.     

Low-rank solution of an optimal control problem constrained by random Navier-Stokes equations

We develop a low-rank tensor decomposition algorithm for the numerical solution of a distributed optimal control problem constrained by two-dimensional time-dependent Navier-Stokes equations with a stochastic inflow. The goal of optimization is to minimize the flow vorticity. The inflow boundary condition is assumed to be an infinite-dimensional random field, which is parametrized using a finite- (but high-) dimensional Fourier expansion and discretized using the stochastic Galerkin finite element method. This leads to a prohibitively large number of degrees of freedom in the discrete solution. Moreover, the optimality conditions in a time-dependent problem require solving a coupled saddle-point system of nonlinear equations on all time steps at once. For the resulting discrete problem, we approximate the solution by the tensor-train (TT) decomposition and propose a numerically efficient algorithm to solve the optimality equations directly in the TT representation. This algorithm is based on the alternating linear scheme (ALS), but in contrast to the basic ALS method, the new algorithm exploits and preserves the block structure of the optimality equations. We prove that this structure preservation renders the proposed block ALS method well posed, in the sense that each step requires the solution of a nonsingular reduced linear system, which might not be the case for the basic ALS. Finally, we present numerical experiments based on two benchmark problems of simulation of a flow around a von Kármán vortex and a backward step, each of which has uncertain inflow. The experiments demonstrate a significant complexity reduction achieved using the TT representation and the block ALS algorithm. Specifically, we observe that the high-dimensional stochastic time-dependent problem can be solved with the asymptotic complexity of the corresponding deterministic problem.