有限单元法模拟地下水流的水头反常现象——对流量计算的作用

地下水流数值模拟是定量研究地下水水量和水质的重要手段。为揭示有限单元法水头反常现象形成的原因,以三维非均质模型为例,刻画了水头反常现象。储量集中法(Lumped Mass方法)是有限单元法中储释水量计算的替代算法,可以改善水头反常现象。对结点控制区域的均衡性质进行了研究,分析了使用Lumped Mass方法时的水头反常现象,结果表明:结点控制区域的对流量计算受到该结点所在层及其相邻层水力坡度的影响;在抽水初始阶段,抽水井结点水头快速下降,导致相邻层中相邻结点控制区域的侧向流入量为正值,控制区域为储水过程,结点水头上升,表现出水头反常现象。因此,有限单元法中对流量的计算是引起水头反常现象的原因。     

The use of natural analogues to assess radionuclide transport

Predicting the long-term safety and performance of a repository for intermediate- to high-level radioactive wastes requires the evaluation of various scenarios which may influence the integrity of the repository system. One such scenario, radionuclide transport, is described, and some of the geochemical processes which serve to enhance or retard transport are illustrated with examples selected from natural analogue or natural system studies. These studies, as distinct from laboratory simulations, help considerably in understanding how near- and far-field radionuclide transport mechanisms work over geological timescales (thousands to hundreds of thousands of years), more in line with the anticipated lifespan of a repository.

Processes addressed include: stability of the spent fuel UO₂ matrix; bentonite backfill as a barrier to radionuclide diffusion; retention of radionuclides by absorption on fracture minerals; transport and diffusion of radionuclides controlled by interconnected bedrock porosity; geochemical influence and behaviour of repository construction materials on radionuclide transport.     

穿透曲线反演管壁溶质通量的实验研究

泉口穿透曲线提供了一个便捷有效的方法来研究岩溶含水层中污染物的运移。在旱季,滞留在石灰岩中的污染物通 过岩壁释放入管道,与管道水相混合并运输到下游泉口,因此,泉口穿透曲线能够反演出污染物释放时的管壁溶质通量。 笔者利用一个基于运移稀释模型的解析解,由40条实验室测量的穿透曲线来反演出相应的管壁溶质通量。结果表明,其基 本形状与穿透曲线相似,管壁透水的非均匀性在穿透曲线上明显表现出来,而在管壁溶质通量上表现为强振荡。因此,管 壁溶质通量比泉口穿透曲线对管壁透水的非均匀性有更高的分辨率。该通量是一个连接泉口穿透曲线和污染物在岩石中早 期分布的中间量,本文提供了一种研究管壁污染物释放过程的新方法。     

河网非恒定水沙数学模型研究

建立了能够适用于多种情况的河网非恒定流水沙数学模型。推导提出了河网泥沙方程组的分级解法,提出了适用于河网计算的汊点分沙模式,分析和对比了多种分沙模式的适用性,对河网区域中可蓄水汊点也提出了相应的处理办法。采用汉江杜家台分蓄洪区的资料作了模型的验证和应用。     

Hydrogeological modeling of radionuclide transport in low permeability media: a comparison between Boom Clay and Ypresian Clay

Deep low-permeability clay layers are considered as suitable environments for disposal of high-level radioactive waste. In Belgium, the Boom Clay is the reference host formation and the Ypresian Clay an alternative host formation for research and safety and feasibility assessment of deep disposal of nuclear waste. In this study, two hydrogeological models are built to calculate the radionuclide fluxes that would migrate from a potential repository through these two clay formations. Transport parameter heterogeneity is incorporated in the models using geostatistical co-simulations of hydraulic conductivity, diffusion coefficient and diffusion accessible porosity. The calculated radionuclide fluxes in the two clay formations are compared. The results show that in the Ypresian Clay larger differences between the fluxes through the lower and the upper clay boundary occur, larger total output radionuclide amounts are calculated and a larger effect of parameter heterogeneity on the calculated fluxes is observed, compared to the Boom Clay.     

油气田上地气物质迁移机制研究

将地气测量技术应用于油气田勘查及评价工作中,地气测量反映出４０００多ｍ的深部构造形成的弱异常,地气异常间接批示了该地区的含气性。地气物质迁移机理研究结果表明,采样器上沉积的物质是以纳米级微粒及其集合体形式存在的。这些纳米级微粒由于具有一系列物理和化学特性,使得它可以与地壳中上升气流的气体相吸咐,随着垂向迁移的气体而上升。不断产生的纳米级微粒被上升气流带到地表,从而在地表形成地气异常。     

Stochastic analysis of the effect of spatial variability of diffusion parameters on radionuclide transport in a low permeability clay layer

Most studies that incorporate subsurface heterogeneity in groundwater flow and transport models only analyze and simulate the spatial variability of hydraulic conductivity. Heterogeneity of the other flow and transport parameters are usually neglected. This approach is often justified, but there are, however, cases in which disregarding the heterogeneity of the other flow and transport parameters can be questionable. In low permeability media, for instance, diffusion is often the dominant transport mechanism. It therefore seems logical to incorporate the spatial variability of the diffusion parameters in the transport model. This study therefore analyses and simulates the spatial variability of the effective diffusion coefficient and the diffusion accessible porosity with geostatistical techniques and incorporates their heterogeneity in the transport model of a low permeability formation. The formation studied was Boom clay (Belgium), a candidate host rock for the deep geological disposal of high-level radioactive waste. The calculated output radionuclide fluxes of this model are compared with the fluxes calculated with a homogeneous model and a model with a heterogeneous hydraulic conductivity distribution. This analysis shows that the heterogeneity of the diffusion parameters has a much larger effect on the calculated output radionuclide fluxes than the heterogeneity of hydraulic conductivity in the low permeability medium under study.     

TEM observation of geogas-carried particles from the Changkeng concealed gold deposit,Guangdong Province,South China

Geogas prospecting is a relatively new technique in exploration for deep-seated and concealed mineral deposits. Now, the technique is to reveal the presence of concealed ore bodies by analyzing the element concentrations of substances that are transported by geogas from concealed ore bodies to the Earth's surface. However, the element concentrations provide us with limited information about the concealed ore bodies. In our work, to investigate the category, size, shape, chemical component of geogas particles, copper grids for transmission electron microscopy measurement were directly used to capture the geogas-carried solid particles in Quaternary sediments overlying the Changkeng concealed gold deposit, Guangdong Province, South China. In the samples, we found Au, PbSO₄, WO₃, Fe₂O₃, SiO₂, TiO₂, and Al₂O₃ particles, etc. using a transmission electron microscopy. Of which, gold particles, and Hg-, Zn-, Pb-, W-bearing particles may originate from concealed ore bodies and can provide more direct information than those provided by element concentrations.     

中国西北地区水汽的平流输送和辐合输送

分析了中国西北地区上空的水汽汇聚过程———平流输送与辐合输送,为西北地区降水过程数值模拟的参数化提供依据.分析表明:中国西北大部分地区为负的水汽平流输送,从上风方向得到的水汽输入小于水汽输出,是水汽的净流失区;中国西北地区中部为大范围负的水汽辐合输送,当地的水汽由风场向外扩散,加剧了当地的水分流失,是水汽扩散区;在水汽净输送的年际变化中,中国西北地区小于华北地区,平流输送呈上升趋势,辐合输送下降,平流的相对贡献在增加;在与大气环流系统的关系中,中国西北地区水汽的平流输送主要受西风带影响,辐合输送与南亚夏季风有显著相关.变化趋势表明西风的作用在加强.     

Stochastic analysis of the effect of heterogeneity and fractures on radionuclide transport in a low-permeability clay layer

Deep low-permeability clay layers are considered as safe environments for disposal of high-level radioactive waste. In Belgium, the Boom Clay is a candidate host rock for deep geological disposal. In this study, we analyze the effects of fractures and spatially variable hydraulic conductivity on radionuclide migration through the clay. Fracture geometry and properties are simulated with Monte Carlo simulation. The heterogeneity of hydraulic conductivity is simulated by direct sequential co-simulation using measurements of hydraulic conductivity and four types of secondary variables. The hydraulic conductivity and fracture simulations are used as input for a transport model. Radionuclide fluxes computed with this heterogeneous model are compared with fluxes obtained with a homogeneous model. The output fluxes of the heterogeneous model differ at most 8% from the homogeneous model. The main safety function of the Boom Clay is thus not affected by the fractures and the spatial variability of hydraulic conductivity.     

我国北方型裂隙岩溶水流及污染物运移数值模拟研究进展

简要回顾了近几十年来我国北方型裂隙岩溶水流及污染物运移数值模拟方面的研究成果,评述了裂隙岩溶水资源评价和污染物运移模拟的计算方法和模拟技术,其中对一些与数值模拟相关的重要问题、近年来研究的热点问题,如模型参数识别及其尺度效应问题、数值弥散和数值振荡问题、Non-Darcy流和Non-Fickian运移问题以及随机方法研究等进行了重点评述,对已有的解决办法进行了比较和对照,并探讨了今后的研究趋势.     

流体-岩石相互作用过程中氧同位素交换地球化学动力学评述

氧同位素研究对于示踪流体—岩石相互作用过程中流体的时间累积流量(或流体／岩石比)、流动方向和组成具有重要意义。基于质量平衡原理可以建立“封闭”体系和开放体系的氧同位素交换模型。“封闭”体系又可分为封闭、批式挥发和瑞利挥发体系。瑞利挥发较批式挥发造成岩石更大的^18O亏损，但在地质过程中两者差异并不显著。开放体系连续模型中氧同位素迁移的机制包括扩散称散和平流。由流体流动速率和扩散称散系数定义的Peclet数决定了上述两种机制在一定尺度上对氧同位素迁移的相对贡献。流体—岩石交换受表面动力学控制，当交换速率快于流体流动时，可认为流体和岩石达到了氧同位素分馆平衡，反之则没有达到平衡。由流体流动速率和流体—岩石反应速率常数定义的Damkoehler数决定了反应接近平衡的程度。如果采用多种矿物相监控，则矿物内部分馆可有效地区分这两种反应模型。对流体—岩石相互作用过程中氧同位素变化的地球化学动力学进行了系统评述，其原理和模型也可扩展到对其他元素的研究。     

Laboratory simulation of diffusion in contaminated marine sediments

Wastes from offshore oil drilling activities are often discharged to the marine environment. Solid wastes that settle onto the bottom sediment may pose a health threat to marine organisms and eventually to man through the food chain. We need to understand their fate in order to predict the chemical concentration levels and life-times in the sediment and adjoining aquatic boundary layer. A laboratory simulation of selected in-bed processes that addresses contaminant leaching from the sediment is proposed. The process chosen for simulation in this study is the coupled desorption-diffusion of contaminants from the bed to the water column. A simple mathematical model of the process is also proposed. Preliminary results using organic chemicals for both the simulation and the model are presented. The results suggest that the experimental procedure represents a good way of estimating the diffusive leaching rates of hydrophobic compounds from contaminated sediments.     

三门峡库区一维非恒定非均匀泥沙输移数学模型

黄河是高含沙河流,含沙量季节差异显著。因此研究黄河干流的非恒定水沙输移规律以及库区泥沙的淤积问题极为重要。应用圣维南方程组以及非恒定泥沙连续方程建立了非恒定非均匀泥沙含沙量计算公式,并根据沙量平衡方程推求出三门峡库区河底高程的变化规律。使用已有的资料对模型进行的验证表明:模型计算与实测资料符合良好,该模型具有较好的应用前景。     

地气测量在叶县—邓县—南漳地学剖面研究中的应用

本剖面横跨华北地台南缘,秦岭造山带和杨子地台北缘。在剖面上系统地采集了壤中气、土壤和气溶胶样品,分别测定Ｃｒ、Ｃｏ、Ｎｉ、Ｃｕ、Ｐｈ、Ｚｎ等微量元素以及Ｈｇ气、Ｈｇ土、Ｒｎ、ＣＨ４、ＣＯ２等气体组分,研究了它们的空间分布特征,并结合剖面沿线的地质和地震资料,讨论了异常与岩性分布、地质构造和找矿靶区之间的关系。研究结果表明,宽坪群和二郎坪群是找寻多金属矿的有利层位;龙王集至龙门集一带是找油气的有利地段;潦河附近的Ｆ（１４）断裂带是一条规模较大的含金破碎带。这次试验的结果还表明,地气测量用于研究厚层的松散层覆盖下的地质构造具有良好的前景。     

北山花岗岩热破裂室内模拟试验研究

岩石热破裂是高放废物地质处置工程中需深入研究的课题。对我国高放废物重点预选场址甘肃北山的花岗岩开展室内热破裂模拟试验研究,采用多通道温度测试仪、声发射、波速层析成像和数码显微镜等手段研究了该花岗岩热破裂过程。试验表明,（1）热破裂从试件端部开始产生,逐步向内缓慢扩展,表现出分段性和独立性;（2）根据声发射撞击率可将热破裂可分为稳定热损伤、宏观裂纹形成、宏观裂纹扩展、裂纹冷却闭合4个阶段,声发射定位的时空演化规律清楚地揭示了裂纹从试件上端部向内部扩展的规律;（3）波速层析成像指示了宏观裂纹位置及高温对岩石造成显著损伤的区域,热应力产生的损伤集中在试件边界,范围小,损伤严重,高温造成的损伤集中在钻孔附近高温区,范围较大,损伤略轻微;（4）监测多通道温度,获得了试件内的温度场并为数值模拟参数选取提供验证,采用有限元程序进行了热力耦合数值模拟,从机制上初步解释了热破裂现象,研究认为综合声发射实时监测热破裂过程和波速层析成像能实现对热损伤的量化的特性可实现岩石热破裂的动态监测和损伤量化,为今后地下实验室相关试验的开展和认识高放废物处置长期稳定性做了有意义的探索。     

波浪作用下沙坝剖面形成过程的数值模拟

建立了波浪、底部离岸流、泥沙运动和沙滩剖面演变耦合数学模型,模拟了实验室中波浪作用下沙坝剖面的形成过程,讨论了各个物理参数对剖面上的水动力和地形变化的影响。通过参数率定,较好地复演了不规则波作用下沙滩地形由均匀斜坡向沙坝剖面的演变,波高、底部离岸流、输沙率和剖面演变的计算结果与实测数据吻合良好。采用数值实验和误差分析方法,讨论了破碎波水滚倾角、泥沙扩散系数和床面休止角等物理参数对计算结果的影响。研究表明,所建立的模型能较好地描述波浪作用下沙坝剖面形成过程中的重要物理机制。     

土石坝溃决模拟及水流计算研究进展

对土石坝溃决模拟技术的最新研究成果进行了总结和评价,在此基础上对该研究领域今后的研究工作提出了若干建议,包括:应研究提出不同坝型溃决可能性分析评价方法;针对不同坝型开展溃坝离心模型试验,揭示在不同致灾因子作用下,土石坝溃决机理和溃口发展过程;开展土石坝初始管涌形成以后发展过程的试验研究,揭示孔流转变为堰流的控制条件;针对溃坝水流的流线曲率较大、溃坝水流大多是非恒定超临界流以及筑坝材料粒径级配范围变化大等特点,开展大型溃坝水工水力学模型试验,揭示不同坝型的溃口流量过程、泥沙输移及下游河道洪水演进规律。     

Application and Current Status of Geogas Prospecting in Concealed Uranium Deposits Exploration

The research achievements of geogas prospecting through the mechanisms of geogas prospecting and its applied research in the area of concealed uranium deposits prospecting were analyzed. Geogas prospecting is a new means of deeply concealed deposits exploration that started to develop gradually from the 1980s and is at an initial stage of research in concealed uranium deposits exploration. Furthermore, comments and suggestions were provided for following research work by summarizing the shortcomings and problems of current studies in this paper. A general upward micro-flow of geogas could carry the metal nanoparticles associated with ore body underground to the surface of earth to come into being anomaly halos of elements contents above the ore body, and with these anomaly halos, the relevant information of deeply concealed uranium deposits (or other classes of metal ore bodies) would be reflected. The depth that geogas measurement could achieve for prospecting concealed deposits under the covering layers is up to a thousand meter. Geogas prospecting can be used for investigating in small scales, surveying generally and in detail in big scales and locating the ore bodies. At the present stage, there are three primary under-mentioned problems in the studies of geogas prospecting: ①On the one hand, the researches of how mediums (covering layer) would affect the migration process of geogas are far from enough. On the other hand, appropriate mathematical theory models have not been set up yet. ②In China, there are currently few accurate and practical physical models for geogas prospecting to carry out controllable and repetitive series experiments as well as theoretical studies. It is essential that an appropriate physical model should be established timely in order to provide an experimental platform and a kind of technical assistance for the mechanism study of geogas prospecting. ③To achieve the aim of popularizing the geogas prospecting in the uranium exploration area, the research of trapping solution and sampling and measuring methods should be further developed and the methods of rapid geogas in-site measuring which are suitable for various geologic structures of uranium deposits should be made reasonably to cater to the need of rapid and accurate field measurement.