岩体裂隙饱和/非饱和渗流机理初步研究

基于裂隙局部开度为平行板模型的假定,按一定分布抽样生成裂隙面上开度随机数,并分别对单裂隙饱和及非饱和渗流进行了数值模拟.裂隙饱和及非饱和渗流在渗流机理上都是由于水势差引起的,只不过在非饱和渗流情形下,湿润相水压力小于非湿润相空气的压力.饱和渗流时,水流主要集中在大开度,即“沟槽”中;而非饱和渗流时,水流运动多发生在小开度,即“狭缝”中.     

含水层层状非均质对地下水流系统的影响

区域尺度上含水层非均质具有复杂的结构性和随机性,难以准确刻画,造成非均质对区域地下水流系统的影响机制研究不够深入。本文以鄂尔多斯盆地白垩系地下水流系统为研究实例,选择典型剖面,采用剖面二维随机数值模拟方法,通过对比不同非均质刻画方法下地下水流场的变化,探讨含水层层状非均质对地下水流系统的影响机制。结果显示,均质条件下模型各向异性（含水层水平和垂向渗透系数比值Kh/Kv）取值为1000时,地下水流场与实际条件较为接近;非均质条件下,渗透系数方差取值0.91,水平相关长度取值5000 m,Kh/Kv取值150时,接近实际条件。研究表明,在大尺度地下水流模拟研究中,采用水平相关长度、渗透系数方差和各向异性值三个变量生成的随机场能很好地刻画含水层的层状非均质特征及其对水流系统的影响控制作用。由于含水层不同尺度层状非均质的叠加效应,采用均质各向异性介质等效概化含水层层状非均质性会造成等效各向异性值偏大失真的效应。     

降雨入渗条件下边坡岩体饱和非饱和渗流计算

简要分析了连续介质饱和非饱和渗流数学模型,并讨论了渗流有限元计算中的有关问题,同时研究了降雨入渗机理及模拟方法,在此基础上编写了非饱和渗流程序SUSC。运用该程序模拟了某边坡降雨过程渗流场的变化情况。计算结果表明,在降雨入渗作用下边坡顶部迅速被饱和,随后表面雨水逐渐向边坡深部下渗,形成从坡顶往深部压力水头等值线由高(零)→低→高的封闭现象。随着降雨的继续,边坡顶部负压区进一步缩小,且负压绝对值减小。降雨结束后,由于上部地下水的继续下渗,在边坡的一定范围、一定时间内压力水头继续升高。根据计算结果可知,局部地方的压力水头最高值出现在降雨结束2 d左右,往后整个边坡非饱和区地下水压力水头全面下降,逐渐恢复原状。模拟结果总体可靠,可作为边坡稳定性分析评价及边坡排水加固的参考依据。同时表明,运用上述饱和非饱和渗流模型及降雨模拟方法,计算降雨条件下边坡岩体的渗流场是可行的。     

GeoStudio软件在土坡饱和—非饱和渗流分析中的应用

在对饱和-非饱和渗流数学模型简要分析的基础上,对影响渗流有限元计算分析结果的一些参数选取方法进行讨论,如非饱和土的水土特征曲线、渗透系数及毛细水头等。利用GeoStud io软件对降雨条件下的土质边坡内部渗流场进行模拟,分析地下水位线、坡体饱和度、压力水头随降雨时间的变化规律。模拟结果与实际情况比较吻合,可为土质边坡稳定性的分析提供参考依据。     

地下水位在非饱和水流数据同化中的应用

为理解地下水位观测信息在非饱和水流数据同化中的数据价值,建立了基于地下水位动态观测信息的一维饱和-非饱和水流集合卡尔曼滤波,通过虚拟数值实验检验了地下水位观测信息在非饱和水力参数估计和水分校正中的潜在价值。研究结果表明:在以地下水位为唯一观测数据时,同时更新参数和水头比仅更新水头能更好地校正土壤剖面的水头分布;当多层单个水力参数未知时,地下水位观测可以为参数估计提供有效信息;当多层多个参数未知时,地下水位与多层多个参数之间的复杂关系导致观测信息难以估计出最优的(唯一的)参数值;地下水位可作为辅助信息,与含水量观测等信息联合运用改善参数估计和含水量预测精度。     

多孔介质随机参数的相关性对非饱和水流的影响

基于van Genuchten-Mualem非饱和水分特征模型,建立了非饱和流运动的随机数值模型。将饱和水力传导度和孔隙大小分布参数视为服从对数正态分布的随机场,用Karhunen-Loeve展开分解,水头表示为混沌多项式展开。通过摄动方法得到一系列关于水头展开的偏微分方程,并用有限差分法进行求解。应用本文的模型分析了两随机场在统计不相关和完全相关模式下对水流随机分析的影响,结果表明两种模式下的水头均值相同,完全相关模式下的水头标准差较不相关模式下的明显偏小。     

垃圾土渗透性和持水性的试验研究

通过室内常水头试验测定了垃圾土的饱和渗透系数,通过室内压力板仪试验测得垃圾土的土-水特征曲线。基于土-水特征曲线预测了垃圾土的非饱和渗透系数,并通过室内入渗试验进行了初步验证。常水头试验得到深层、中层和浅层垃圾土的饱和渗透系数分别为3.56×10-4、3.50×10-3、4.81×10-2 cm/s。土-水特征曲线试验表明,垃圾土饱和含水率和残余含水率较高,进气值很小,土-水特征曲线在低基质吸力时存在陡降段,其中浅层垃圾土的陡于深层垃圾。验证试验表明,预测得到的非饱和渗透系数与实测结果接近,基于土-水特征曲线预测垃圾土非饱和渗透系数的方法基本可行。     

降雨入渗条件下非饱和土朗肯土压力分析

传统的土压力分析仅考虑了土体饱和强度对土压力产生的贡献,忽略了基质吸力及其变化对土压力的影响。运用非饱和土有效应力原理和饱和土朗肯土压力公式推导了非饱和土朗肯土压力公式,结合Iverson降雨入渗解析解,推导出降雨入渗条件下非饱和土压力公式。该公式将降雨入渗时的非饱和土压力表示为时间和深度的函数,更符合实际情况。研究结果表明：采用该方法计算得到的土压力值相对于传统计算结果偏大,作用点偏高;此外,随着降雨的发生、入渗和停止,主动土压力呈现“减小-增大-减小-稳定”趋势,被动土压力呈现“增大-减小-增大-稳定”趋势,该现象由降雨过程中基质吸力改变所致。由该公式获得的土压力分布及变化规律可用于挡土工程结构的设计。     

黄河大堤非饱和土土-水特性试验研究

模拟工程中非饱和土体的实际受力状态,利用改进的非饱和土三轴仪对黄河大堤非饱和土在不同围压条件下的土-水特性进行了试验研究。根据试验结果,分析了黄河大堤非饱和土体在外荷作用下土-水特征曲线及其变化规律。并对相应的土-水特征曲线模型进行了拟合,得出其拟合函数。根据土-水特征曲线可知,非饱和土的基质吸力随着土体所受周围压力的增大而减小,随着含水量的增大而减小;对于某一类土而言,在高含水量(含水量＞20 %) 时,土的基质吸力随含水量的变化较小,而在低含水量 (含水量≤20 %)时,土的基质吸力随含水量的变化较大。随着含水量的增加,非饱和土的工程性质向弱性变化。     

水位波动下非饱和心墙土坝体的渗流和稳定性

基于Richards方程,对坝体的饱和-非饱和渗流场进行了模拟,再根据饱和-非饱和渗流场和非饱和土抗剪强度公式,对坝体的稳定性进行了分析。结果表明,水位骤降过程中坝体的安全系数通常呈现先缓慢增加后迅速减小的变化过程,分析坝体失稳时塑性区和位移场发现,水位下降的初期,坝体左侧坡体的安全系数要低于坝体右侧坡体,但水位下降到一定程度,右侧坡体的安全系数迅速减小,并先于左侧坡体失稳;采用有限元强度折减法用于多坡面边坡稳定分析时,只能获得最小安全系数的包络线;心墙具有隔水防渗的作用,对水位变化渗流具有阻尼作用。     

Understanding the effect of single-fracture heterogeneity from single-well injection-withdrawal (SWIW) tests

The single-well injection-withdrawal (SWIW) tracer test is a method used to estimate the tracer retardation properties of a fracture or fracture zone. The effects of single-fracture aperture heterogeneity on SWIW-test tracer breakthrough curves are examined by numerical modelling. The effects of the matrix diffusion and sorption are accounted for by using a particle tracking method through the addition of a time delay added to the advective transport time. For a given diffusion and sorption property (P _m) value and for a heterogeneous fracture, the peak concentration is larger compared to a homogeneous fracture. The cumulative breakthrough curve for a heterogeneous fracture is similar to that for a homogeneous fracture and a less sorptive/diffusive tracer. It is demonstrated that the fracture area that meets the flowing water, the specific flow-wetted surface (sFWS) of the fracture, can be determined by matching the observed breakthrough curve for a heterogeneous fracture to that for a homogeneous fracture with an equivalent property parameter. SWIW tests are also simulated with a regional pressure gradient present. The results point to the possibility of distinguishing the effect of the regional pressure gradient from that of diffusion through the use of multiple tracers with different P _m values.     

A practical method for identifying parameters in the seismic design of embankments

The effect of heterogeneity of soils on a seismic design procedure based on residual displacement according to safety factors is examined. Although soils are heterogeneous in nature, at present the seismic designs of embankments implement the assumption that soils are homogeneous. In this study, a seismic reliability analysis using a Monte Carlo simulation and considering the heterogeneity of soil strength was applied to the Newmark sliding block method. The data of the analysis in this study were based on in-situ soil test results of actual airport embankments. This paper indicates a rational and practical approach to modifying the actual heterogeneous soil conditions into homogeneous design parameters for the Newmark sliding block method in homogeneous ground conditions.     

Computer simulation of Pb and Sr isotope evolution of the Earth's crust and upper mantle

Isotope evolution in a differentiated (crust and upper mantle) and chemically heterogeneous Earth has been computed for a model with isotopic exchange between crust and mantle at an exponentially decreasing rate. To simulate the effects of subduction, cells of crust and mantle are selected at random, their contents mixed and then redistributed into the cells from which they came—to give new chemically heterogeneous but momentarily isotopically homogeneous systems. Daughter isotopes in each cell grow according to equations appropriate for closed chemical systems. Rock age distributions and isotopic data created by the computer calculation mimic nature. Pb isotope data changes through geologic time are illustrated to demonstrate that two-stage interpretations applied to Pb data for rocks with complex histories may be misleading. The intercept of the Pb ore growth curve and regressions fit to Pb data gives minimum values only for the duration of heterogeneous $\text{U}\text{Pb}$ systems, not the time when heterogeneous distribution first occurred. An intercept derived time of about 3 b.y., in the Pb system is shown to be quantitatively compatible with an average $\text{Rb}\text{Sr}$ age of crustal rocks of 1.5 b.y. and with a constant degree of chemical heterogeneity for all of Earth history.     

Seasonal changes in basal conditions at Briksdalsbreen, Norway: the winter–spring transition

The winter–spring transition is a dynamic time within the glacier system, because it marks a period of instability as the glacier undergoes a change in state from winter to summer. This period is normally associated with sudden pressure fluctuations resulting in hydrological instabilities within the subglacial drainage system. New data are presented from wireless multi-sensor subglacial probes incorporated within the till at Briksdalsbreen, Norway. Water pressure readings recorded a two-phase winter–spring transition. Event 1 occurred early in the year (December–January) and marked the start of activity within the subglacial environment following the winter. However, this did not result in any permanent changes in subglacial activity and was followed by a period of quiescence. Event 2 occurred later in the year in accordance with changing external weather conditions and the retreat of the snow pack. It was characterized by high-magnitude pressure peaks and diurnal oscillations in connected regions. The variations in sensor trends that followed this event suggested that a transition in the morphology of the subglacial drainage system had occurred in response to these pressure fluctuations. Event 2 also showed some similarities with spring events recorded at valley glaciers in the Alps. A conceptual model is presented associating the form of the winter–spring transition with respect to the location of the probes within connected and unconnected regions of the subglacial drainage system. These data provide further evidence for temporal and spatial heterogeneous subglacial drainage systems and processes. The identification and analysis of subglacial activity during the winter–spring transition can contribute to the interpretation of hydro-mechanical processes occurring within the subglacial environment and their effect on glacier dynamics.     

条形荷载作用下梯度饱和土的动力响应分析

根据Biot波动理论,研究了条形均布荷载作用下非均匀饱和土地基的动力响应问题。利用Fourier积分变换,通过Helmholtz矢量分解原理,建立了饱和土层在动荷载作用下的回传射线矩阵法计算列式,考虑饱和土的物理力学性质沿深度方向按幂函数连续变化,采用数值Fourier逆变换获得了饱和土地基的位移、应力和孔隙压力等物理量的数值解。分析讨论了材料非均匀性对饱和土介质动力特性的影响。结果表明,非均匀饱和土的动力行为与均匀饱和土有着明显的不同,当土体的非均匀程度越高,条形荷载中点下流体压力和应力幅值越大,而位移、流体压力以及应力等物理量在水平方向的振动频率均随着土体非均匀变化程度的增强而增大。     

Deposition of heterogeneous radiocolloid from groundwater on enclosing rocks

A model of precipitation of radiocolloid particles on enclosing rocks at colloid-facilitated transport of radionuclides by groundwater is considered. The model proposed is based on probabilistic analysis of deposition of heterogeneous colloidal particles on fixed sites of groundwater filtration channels, which are formed by connected systems of pore_ and fracture voids in enclosing rocks. The model takes into account heterogeneity of the system of colloidal particles and competition of different particles for the same potential deposition sites. Theoretical analysis is based on the solution of the system of Kolmogorov-Feller equations. The limit values of the obtained solutions correspond to the equilibrium distribution of radiocolloid particles of different sizes between the groundwater and the enclosing rocks. In the case of a homogeneous system of particles, the obtained distribution is reduced to the Langmuir isotherm.     

含水层非均质性不同刻画方法对地下水流和溶质运移预测的影响

为探讨含水层非均质性不同刻画方法对地下水流和溶质运移预测的影响,基于非均质含水层砂箱实验,分别用传统等效均质模型、克立金插值和水力层析刻画含水层渗透系数场,并探讨了先验信息对水力层析结果的影响.将不同方法估算的渗透系数场用以预测地下水流和溶质运移过程,以此判断不同方法估算结果的优劣,分析含水层非均质性对地下水流和溶质运移的影响.结果表明:与克立金插值法相比,水力层析法可以更好地刻画含水层非均质性,较准确地预测地下水流和溶质运移过程;钻孔岩心渗透系数样本值作为先验信息可以提高水力层析法估算结果的精度;传统等效均质模型无法准确预测地下水流和溶质运移过程.含水层非均质性的增强将导致溶质污染羽分布形态和运移路径的空间变异性增强,并且优势通道直接决定溶质的分布及运移路径.      

溶蚀作用下古岩溶盆地系统中介质场演化模拟

为了研究古岩溶盆地在溶蚀作用下碳酸盐岩油气储层介质场的演化规律,基于地下水渗流理论和碳酸盐溶蚀动力理论,用数值模型模拟了均质盆地和非均质盆地2种岩溶盆地含水系统发育过程。岩溶盆地含水系统演化过程中总伴随着裂隙差异溶蚀,不管均质盆地还是非均质盆地,都越来越非均质化,强溶蚀带集中在潜水位及优势裂隙附近,在侵蚀基准面和构造裂隙及层理等大裂隙处会形成高孔隙率、渗透率的良好储层。受介质场非均质化反作用,岩溶高地区潜水位随系统演化不断下降,潜水面处的强烈溶蚀随水位不断下切而使高地区改造成竖直裂隙发育的厚储层;坡地区水平径流活跃,易于形成水平裂隙发育的储层;平原区水位相对稳定,溶蚀作用主要发生大裂隙和侵蚀基准面处,但在大裂隙网内经溶蚀可形成局部高孔隙率储层。对比两子模型发现,在岩溶盆地含水系统中,主导渗流场都要经历由局部流场向二级流场、二级流场向全局流场的转变,前者发生在坡地区内部,后者发生在坡地区与平原区之间;流畅的渗流场转换更有利于介质场发育,受大裂隙网导水作用,非均质盆地比均质盆地更快捷、更顺利地实现主导流场升级,溶蚀作用更强烈,3 ka后非均质盆地比均质盆地的孔隙率增幅大34%。