破碎、较破碎岩体渗透系数与埋深关系的研究——以某石灰石矿山为例

岩体的渗透性与其埋深关系紧密,一般表现为岩体的渗透性随埋深增加而逐渐减小的规律,本文开展破碎、较破碎岩体随单一埋深关系的渗透系数研究,采用负指数相关性模型和幂次函数模型对丽江古城某石灰石矿山在地表浅部55 m深度范围进行的5次压水水样数据进行拟合相关模型计算参数,发现幂次函数模型的相关性(R=0.916 5)比负指数模型相关性(R=0.884 2)优越。并根据拟合参数与CHEN Y F(2018)通过我国西南地区12个大坝上万组注水试验数据得出的拟合参数进行对比,得出基于常数头注水试验在地表以下较浅深度段的渗透系数可以代替地表渗透系数,采用多段不同埋深渗透系数拟合K₀时宜采用负指数相关模型,确定地表渗透系数后或已知地表渗透系数的情况下宜采用幂次函数模型预测不同深度渗透系数。     

一种基于地质指标的裂隙岩体渗透系数估算模型

裂隙岩体的渗透特性受控于裂隙的发育特征、连通特性和充填情况,并与岩体的地应力水平具有显著的相关性。基于裂隙岩体渗透性的影响因素,并考虑现有渗透系数估算模型的不足,利用钻孔压水试验和钻孔电视图像资料,建立考虑埋深(Z)、岩石质量指标(RQD)以及充填物指标(FSD)等3个指标的渗透系数估算ZRF模型,并应用到牙根二级水电站及其他工程区的渗透系数估算中。结果表明,与已有的渗透系数估算模型相比,ZRF模型较好地反映了岩体渗透性的影响因素,且模型参数物理意义明确,便于获取,对分析裂隙岩体渗透性具有一定的工程参考价值。     

利用电导率测井与压水试验联合评价岩体渗透性的方法

进行地下水封洞库的水封条件分析评价必须获得准确可靠的建库岩体渗透性参数,而获取岩体渗透系数常用的传统水文地质试验方法存在明显的不足。为了试验数据的精确性,文章基于广义径向流（GRF）理论,依托烟台某地下水封洞工程,以丙烷洞库交通巷道钻孔为例,开展压水试验并采用非稳定流理论的GRF模型优化解析试验数据,结合电导率测井试验确定导水裂隙位置并求出裂隙范围内的渗透系数。试验结果表明:GRF模型比稳定流模型解析结果大1~2倍。原因在于裂隙岩体进行分段压水时,各段水流维度不一,传统稳定流理论假设水流维度只有二维流,而GRF模型为空间n维裂隙流用压水过程全部数据进行拟合,不同时段在相应维度下进行计算,因此其求算的渗透系数K更接近于试验段真值,具有更好的兼容性和实用性。同时利用电导率测井试验计算长度（导水裂隙范围）远远小于压水试验段计算长度的特点,可将GRF模型解析得出的分段渗透系数做进一步细化平均以提高压水试验解析结果的精度,为水封洞库效果评价、洞库涌水量预测提供更加科学可靠的数据基础。     

磁共振探测估算含水层渗透系数的原位试验研究

通过一个完整的磁共振探测（MRS）、钻探、抽水试验过程,开展了MRS估算含水层渗透系数的原位试验研究。在分析MRS估算渗透系数准确性的基础上,系统剖析了引起准确性差异的主要因素,并对进一步提高估算效果提出了建议。结果显示：（1）采用现阶段常用的参数设置,与抽水试验计算值相比,MRS估算渗透系数的差值为抽水试验计算值的17.59%。（2）MRS推断的含水层顶板、底板埋深与钻孔揭示的信息相比,差值分别为4.11m、1.03m,表明其对估算渗透系数准确性的影响较小。（3）CP值是影响估算准确性的重要因素,其值为10-9数量级符合大多数地层的特点。另外,通过在已知渗透系数的钻孔附近进行MRS,从而获取CP参考值,应用该值估算的渗透系数准确性高于常用的参数设置。（4）指数a、b设置方面,应用Seevers公式（a=1、b=2）的估算效果优于Kenyon公式（a=4、b=2）。本成果有助于提高MRS估算渗透系数方法在野外条件下的适用性。     

一种面向对象的最优分割尺度计算模型

作为信息提取和分类的前提,面向对象的影像分割尺度参数的设置直接影响到提取和分类的精度。本文以GF-2影像数据为例,在已有分割理论和方法的基础上提出一种基于最优分割尺度的计算模型（OS模型）。该模型以主成分分析所得的主成分以及新建的归一化植被指数（normalized vegetation index,NDVI）特征层作为分割参考层,综合考虑均质因子的影响,构建加权尺度评价指数,插值拟合最优分割尺度。构建误差系数（E_c）对模型进行评价,结果表明：OS模型误差系数（E_c=1.15%）小于传统模型（E_c=3.28%）,且分割对象更均匀、与实际地物更接近。     

基于SUFI-2算法和SWAT模型的妫水河流域水文模拟及参数不确定性分析

敏感性分析和不确定性分析是分布式水文模型参数校准和模型应用的基础。本文以妫水河流域为例,基于SWAT模型和SUFI-2算法进行模型参数敏感性分析,结合手动调参和自动率定,通过SUFI-2的P因子和R因子进行模型不确定性分析,构建妫水河流域分布式水文模型。本次妫水河流域月尺度模拟中:率定期,确定系数R2=0.59,效率系数NSE=0.56;验证期,确定系数R2=0.82,效率系数NSE=0.80;P因子均大于0.5,R因子均小于0.5。结果表明,妫水河流域SWAT模型水文模拟效果较好。     

新场岩体定压力压水试验数据解译方法的对比

以高放废物地质处置北山预选区新场岩体三个深钻孔的定压力压水试验为例,分别从压水阶段数据求解、压力恢复阶段数据求解和特殊试验段参数求解3个方面,对稳定流公式法和非稳定流参数拟合法进行了对比,得到了如下结论：非稳定流参数拟合法较稳定流公式法适用的渗流类型更多,使用数据量更大,结果更为可靠;在压水阶段和压力恢复阶段数据完整时,优先选择压力恢复阶段非稳定流参数拟合结果作为试验段的渗透系数.依据上述筛选原则,得到了新场岩体BS17～BS19三个深钻孔渗透系数的概率分布图,可初步判断新场岩体为低渗透性岩体,适宜于高放废物地质处置.     

土壤水分特征曲线模型模拟性能评价

土壤水分特征曲线模型作为实验测定土壤水分特征数据的一种替代方法,因其具有计算方便快捷和便于嵌入数值模拟程序的优点,开始受到越来越广泛的关注。虽然文献中存在众多的土壤水分特征曲线模型,但是这些模型的适用范围及拟合性能尚不明确。为了获得更加准确适用的土壤水分特征曲线,在实际应用中通常需要花费大量时间和精力去测试各种模型。为了解决上述问题,在国内外研究成果的基础上收集整理了12种典型的土壤水分特征曲线模型,并利用包含不同质地、有机质含量及容重的8种土壤的实测土壤水分特征数据来评估比较这些模型的模拟性能。模型性能通过均方根误差（RMSE）、平均偏差（AD）、AIC准则（Akaike Information Criterion）和纳什效率系数（NSE）4个指标评估。研究结果表明：大部分的模型能够提供比较接近于实际的拟合结果,评价指标值也比较相近。其中,KCGS2006（包含3个参数）和K1999模型（包含2个参数）拟合效果最好,而Gregson1987（包含1个参数）的拟合效果最差。该研究可以深入了解各种土壤水分特征曲线模型的适用性与局限性,更好地为生态环境建设和农业可持续发展研究中土壤水力参数的选取提供依据和参考。     

利用RQD估算岩体不同深度的平均渗透系数和平均变形模量

渗透系数是进行裂隙岩体渗流模拟的必备参数,变形模量是工程岩体数值模拟的必备参数。大量研究中忽视了渗透系数和变形模量随深度变化这一重要规律,从而影响模拟结果的可靠性。为此,探讨了利用极易获取的RQD（岩石质量指标）资料估算不同深度的渗透系数和变形模量的可行性。通过以某花岗岩体为例,研究发现,RQD均值随深度增大,渗透系数均值随深度减小,其相关性很好。因此,利用RQD估算不同深度的平均渗透系数是可行的。根据RQD随深度的变化,利用经验公式估算了不同深度的变形模量均值和变化范围。估算得到的变形模量与实测结果基本一致。     

地下水污染强度及渗透系数场的反演识别研究

在制定地下水污染修复方案时,污染源参数和渗透系数场是最重要的地下水数值模型参数,但前人研究多集中于单一类型参数的识别。文章中采用地下水污染物运移模型（MT3DMS）和数据同化方法（迭代局部更新集合平滑器,ILUES）构成地下水污染源识别的求解框架,并利用Karhunen-Loève展开技术实现渗透系数场的参数降维,最后通过同化水头与浓度数据实现地下水污染源强和渗透系数场的联合反演。结果表明：（1）ILUES算法能精确识别污染源参数和渗透系数场,并且具有很高的普适性;（2）精确表征渗透系数在空间上呈现出的非均质性,是预测污染物迁移路径、反演污染强度的关键;（3）ILUES算法参数影响着反演效果,综合考虑计算效率和计算精度等,可以得到算例的最佳样本集合大小（Ne=4 000）和ILUES算法最佳参数组合（局部临近样本集合占比α=0.4,相对权重b=4）。但在实际工程案例中,如果对精度的要求不是过高,经验组合（α=0.1,b=1）更值得推荐。研究结果对于区域地下水资源调查、评价和管理等工作具有较强的实践意义,并可为后期地下水污染预测及地下水监测井网优化提供技术支撑。     

三峡库区黄土坡滑坡非饱和水力参数研究

非饱和水力参数在计算滑坡降雨入渗过程与稳定性时是至关重要的材料参数。在三峡库区黄土坡滑坡上进行双环渗透试验,获取黄土坡滑坡表土层的饱和渗透系数。对黄土坡滑坡表土层的含水率和基质吸力进行实时监测,采集了黄土坡滑坡表土层中含水率和基质吸力随时间的变化数据,采用van Genuchten土-水特征曲线模型拟合了4个实时监测剖面的土-水特征曲线及其拟合参数。将饱和渗透系数与土-水特征曲线拟合参数代入van Genuchten渗透系数函数模型,求出了黄土坡滑坡表土层在非饱和条件下的渗透系数函数,为黄土坡滑坡在降雨作用下的稳定性计算提供了可靠的水力参数     

Numerical simulations for effects of anisotropy and heterogeneity on well responses to slug tests

Because of their simplicity, assumptions of isotropy and homogeneity are often applied to slug test analyses, determining hydraulic conductivity of heterogeneous and anisotropic formations, in spite of their inherent unfitness. In this study, it was examined how anisotropy and heterogeneity of the formation affect well recovery curves and therefore estimates of hydraulic conductivity in slug tests. Using a finite element method (FEM) slug test model, several hypothetical slug tests in anisotropic and heterogeneous formations with known hydraulic parameters and actual formations of unknown parameters were analyzed. Error factors as a result of the assumptions of isotropy and homogeneity were quantified through sensitivity analyses of well recovery curves. Results of slug tests analyzed in this study indicate that: (1) the well recovery curves in slug tests do not show unique shape or signs reflecting the anisotropy and heterogeneity of the tested formation; (2) the assumption of isotropy does not deviate estimates of hydraulic conductivity significantly when the radial hydraulic conductivity is larger than the vertical one in the formation; (3) the assumption of homogeneity in a layered heterogeneous system has a high potential for errors in hydraulic conductivity estimation, however it can be avoided by placing the screened section apart from the adjacent layers; (4) a low permeability well skin skews the hydraulic conductivity of the original formation, depending on the size and hydraulic conductivity of the skin.     

An empirical model for estimating hydraulic conductivity of highly disturbed clastic sedimentary rocks in Taiwan

This paper presents the establishment of an empirical HC model for estimating rock mass hydraulic conductivity of highly disturbed clastic sedimentary rocks in Taiwan using high-resolution borehole acoustic televiewer and double packer hydraulic tests. Four geological parameters including rock quality designation (RQD), depth index (DI), gouge content designation (GCD), and lithology permeability index (LPI) were adopted for establishing the empirical HC model. To verify rationality of the proposed HC model, 22 in-situ hydraulic tests were carried out to measure the hydraulic conductivity of the highly disturbed clastic sedimentary rocks in three boreholes at two different locations in Taiwan. Besides, the model verification using another borehole data with four additional in-situ hydraulic tests from similar clastic sedimentary rocks was also conducted to further verify the feasibility of the proposed empirical HC model. The field results indicated that the rock mass in the study area has a conductivity between the order of 10^− 10 m/s and 10^− 6 m/s at the depth between 34 m and 275 m below ground surface. Results demonstrate that the empirical HC model may provide a useful tool to predict hydraulic conductivity of the highly disturbed clastic sedimentary rocks in Taiwan based on measured HC-values.     

Study on the selection of unsaturated flow model for the different types of soil and soft rock

Precise estimation of unsaturated hydraulic properties of porous media is indispensable in various study areas, such as analyzing the moisture flow, the drying process occurring from the surface, and the pollutant migration beneath the ground surface. Although many empirical/theoretical models describing the unsaturated hydraulic properties have been proposed by several previous researchers, the best model for the different types of soil/rock may not be identical. Thus, the model selection process and the estimation technique of the parameters included in the models should be developed. In the present study, the inverse technique based on the transient evaporation change was investigated to select the model and estimate the model parameters. The experimental work was based on a relatively low permeable soft rock and a relatively high permeable sandy soil (Toyoura standard sand). Experimental equipment was developed to precisely measure the evaporation rate for the high permeable sandy soil. The Genetic Algorithm (GA) was adopted in the inverse technique as an optimization tool. In order to simplify the problem, only the drying process from the saturated condition was considered. It was established that the information concerning the transient evaporation change could be used for the model selection and parameter estimation. Further, the saturation distribution could be used for the selection of the models. The present study provides important information for the development of the model selection process.     

基于遗传算法-支持向量机的滑坡渗透系数反演

求解库岸边坡岩土体的渗透系数是研究滑坡渗流场及多场演化的基础,一般通过原位试验和室内试验求得,但试验成本较高且试验位置具有一定的随机性。本文以三峡库区马家沟滑坡为例,提出一种利用地下水位动态观测资料反演滑坡岩土层渗透系数的方法。具体步骤为:（1）依据滑坡的勘察资料和水位观测数据,构建滑坡数值模型;（2）利用SPSS生成不同渗透系数正交试验组合,并将渗透系数代入数值模型中计算监测井的水位,得到不同渗透系数及其对应的模拟水位数据;（3）应用遗传算法优化的支持向量机构建坡体模拟水位与渗透系数的非线性映射关系,再通过代入实际动态监测水位值求得滑坡岩土层的渗透系数;（4）将求得的渗透系数代入数值模型,用计算的模拟水位与实际观测水位进行对比验证。研究结果表明:遗传算法优化的支持向量机具有良好的学习预测效果,能准确预测渗透系数与水位的关系。该反演方法具有高效、准确的优点,反演结果的精度满足实际应用需要。     

Classifying zones of suitability for manual drilling using textural and hydraulic parameters of shallow aquifers: a case study in northwestern Senegal

A method is proposed that uses analysis of borehole stratigraphic logs for the characterization of shallow aquifers and for the assessment of areas suitable for manual drilling. The model is based on available borehole-log parameters: depth to hard rock, depth to water, thickness of laterite and hydraulic transmissivity of the shallow aquifer. The model is applied to a study area in northwestern Senegal. A dataset of boreholes logs has been processed using a software package (TANGAFRIC) developed during the research. After a manual procedure to assign a standard category describing the lithological characteristics, the next step is the automated extraction of different textural parameters and the estimation of hydraulic conductivity using reference values available in the literature. The hydraulic conductivity values estimated from stratigraphic data have been partially validated, by comparing them with measured values from a series of pumping tests carried out in large-diameter wells. The results show that this method is able to produce a reliable interpretation of the shallow hydrogeological context using information generally available in the region. The research contributes to improving the identification of areas where conditions are suitable for manual drilling. This is achieved by applying the described method, based on a structured and semi-quantitative approach, to classify the zones of suitability for given manual drilling techniques using data available in most African countries. Ultimately, this work will support proposed international programs aimed at promoting low-cost water supply in Africa and enhancing access to safe drinking water for the population.     

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.     

Application of statistical approaches to analyze geological,geotechnical and hydrogeological data at a fractured-rock mine site in Northern Canada

Mine site characterization often results in the acquisition of geological, geotechnical and hydrogeological data sets that are used in the mine design process but are rarely co-evaluated. For a study site in northern Canada, bivariate and multivariate (hierarchical) statistical techniques are used to evaluate empirical hydraulic conductivity estimation methods based on traditional rock mass characterisation schemes, as well as to assess the regional hydrogeological conceptual model. Bivariate techniques demonstrate that standard geotechnical measures of fracturing are poor indicators of the hydraulic potential of a rock mass at the study site. Additionally, rock-mass-permeability schemes which rely on these measures are shown to be poor predictors of hydraulic conductivity in untested areas. Multivariate techniques employing hierarchical cluster analysis of both geotechnical and geological data sets are able to identify general trends in the data. Specifically, the geological cluster analysis demonstrated spatial relationship between intrusive contacts and increased hydraulic conductivity. This suggests promise in the use of clustering methods in identifying new trends during the early stages of hydrogeological characterization.