Electrical-hydraulic conductivity model for a weathered-fractured aquifer system of Olbanita,Lower Baringo Basin,Kenya Rift

Groundwater yield in the Kenya Rift is highly unsustainable owing to geological variability. In this study, field hydraulic characterization was performed by using geo-electric approaches. The relations between electrical-hydraulic (eh) conductivities were modeled hypothetically and calibrated empirically. Correlations were based on the stoch-astic models and field-scale hydraulic parameters were contingent on pore-level parameters. By considering variation in pore-size distributions over eh conduction interval, the relations were scaled-up for use at aquifer-level. Material-level electrical conductivities were determined by using Vertical Electrical Survey and hydraulic conductivities by analyzing aquifer tests of eight boreholes in the Olbanita aquifer located in Kenya rift. VES datasets were inverted by using the computer code IP2Win. The main result is that ln T=0.537 (ln Fa)+3.695; the positive gradient indicating eh conduction through poresurface networks and a proxy of weathered and clayey materials. An inverse (1/F-K) correlation is observed. Hydraulic parameters determined using such approaches may possibly contri-bute significantly towards sustainable yield management and planning of groundwater resources.     

咸淡水驱替过程中含水介质渗透性变化的试验研究

在野外水文地质调查的基础上,采集青岛市大沽河下游咸水入侵区砂样,首先对含水介质的粒度和矿物组成进行了分析,然后通过室内砂槽模拟试验,对咸淡水驱替过程中含水介质的渗透性变化规律进行了研究。结果表明,当咸淡水相互驱替时,含水介质的渗透性会发生显著变化,这种变化主要是由于驱替液盐浓度的变化,使得伊利石、高岭石和绿泥石等非膨胀性粘土矿物经释放、迁移、絮凝和沉积,从而重新分布所引起的。含水介质渗透性的变化具有明显的非均质性,即位于距入水处不同水平距离以及不同高度处的含水介质,其渗透性变化规律有所不同。     

The differential system method for parameter identification; unconfined aquifer case

We present the applicability of differential system (DS) method for identification of hydraulic conductivity and effective porosity in a phreatic aquifer. In the original setting, the first step of the DS system is to solve an overdetermined algebraic system in the least squares sense. A natural extension of the method is to pose a least squares problem in an appropriate functional space. We show an improvement of the identification by considering the least square problem in the space of square integrable functions in the time variable for a finite interval.     

Peralkaline syenite autoliths from Kilombe volcano, Kenya Rift Valley: Evidence for subvolcanic interaction with carbonatitic fluids

Mineral chemistry, textures and geochemistry of syenite autoliths from Kilombe volcano indicate that they crystallized in the upper parts of a magma chamber from peralkaline trachytic magmas that compositionally straddle the alkali feldspar join in the “residuum system” (ne = 0–1.03; qz = 0–0.77). Mineral reaction and/or overgrowth processes were responsible for the replacement of (i) Mg–hedenbergite by aegirine–augite, Ca–aegirine and/or aegirine, (ii) fayalite by amphiboles, and (iii) magnetite by aenigmatite. Ti–magnetite in silica-saturated syenites generally shows ilmenite exsolution, partly promoted by circulating fluids.

By contrast, the Fe–Ti oxides in the silica-undersaturated (sodalite-bearing) syenites show no signs of deuteric alteration. These syenites were ejected shortly after completion of crystallization. Ilmenite–magnetite equilibria indicate fO₂ between − 19.5 and − 23.1 log units (T 679–578 °C), slightly below the FMQ buffer. The subsequent crystallization of aenigmatite and Na-rich pyroxenes suggests an increase in the oxidation state of the late-magmatic liquids and implies the influence of post-magmatic fluids.

Irrespective of silica saturation, the syenites can be divided into (1) “normal” syenites, characterized by Ce/Ce ratios between 0.83 and 0.99 and (2) Ce-anomalous syenites, showing distinct negative Ce-anomalies (Ce/Ce 0.77–0.24). “Normal” silica-saturated syenites evolved towards pantelleritic trachyte. The Ce-anomalous syenites are relatively depleted in Zr, Hf, Th, Nb and Ta but, with the exception of Ce, are significantly enriched in REE.

The silica-saturated syenites contain REE–fluorcarbonates (synchysite-bastnaesite series) with negative Ce-anomalies (Ce/Ce 0.4–0.8, mean 0.6), corroded monazite group minerals with LREE-rich patches, and hydrated, Fe- and P-rich phyllosilicates. Each of these is inferred to be of non-magmatic origin. Fractures in feldspars and pyroxenes contain Pb-, REE- and Mn-rich cryptocrystalline or amorphous material. The monazite minerals are characterized by the most prominent negative Ce-anomalies (Ce/Ce_mean = 0.5), and in the most altered and Ca-rich areas (depleted in REE), Ce/Ce is less than 0.2.

It is inferred that carbonatitic fluids rich in F, Na and lanthanides but depleted in Ce by fractional crystallization of cerian pyrochlore, percolated into the subvolcanic system and interacted with the syenites at the thermal boundary layers of the magma chamber, during and shortly after their crystallization.

Chevkinite–(Ce), pyrochlore, monazite and synchysite-bastnaesite, occurring as accessory minerals, have been found for the first time at Kilombe together with eudialyte, nacareniobsite–(Ce) and thorite. These latter represent new mineral occurrences in Kenya.     

The hydraulic conductivity field and groundwater flow in the unconfined aquifer system of the Keta Strip,Ghana

This study investigates the hydraulic conductivity field and the groundwater flow pattern as predicted by a calibrated steady state groundwater flow model for the Keta Strip, southeastern Ghana. The hydraulic conductivity field is an important parameter in evaluating aquifer properties in space, and in general basin-wide groundwater resources evaluation and management. This study finds that the general hydraulic conductivity of the unconsolidated unconfined aquifer system of the Keta Strip ranges between 2 m/d and 20 m/d, with an average of 15 m/d. The spatial variation in horizontal hydraulic conductivity appears to take the trend in the variations in the nature of the material in space. Calibrated groundwater recharge suggests that 6.9–34% of annual precipitation recharges the shallow aquifer system. This amount of recharge is significant and suggests high fortunes in terms of groundwater resources development for agriculture and industrial activities in the area. A spatial distribution of groundwater recharge from precipitation is presented in this study. The spatial pattern appears to take the form of the distribution in horizontal hydraulic conductivity, and suggests that the vertical hydraulic conductivity takes the same pattern of spatial variation as the horizontal hydraulic conductivity. This is consistent with observations in other areas. The resulting groundwater flow is dominated by local flow systems as the unconfined system is quite shallow. A general northeast – southwest flow pattern has been observed in the study area.     

Pilot-scale field model tests for detecting landfill leachate intrusion into the subsurface using a grid-net electrical conductivity measurement system

The grid-net electrical conductivity measurement system for detecting exact locations of landfill leachate intrusion in the subsurface was developed in this study. Laboratory and pilot-scale field model tests were performed to evaluate the direct application of a grid-net electrical conductivity measurement system for the detection of landfill leachate. A significant increase in electrical conductivity of soil was observed by adding landfill leachate. This can be explained as an increase in electrical conductivity of pore fluid due to an increase in leachate constituents as charge carriers. In pilot-scale field model tests, leachate intrusion locations were accurately identified at the initial stage of landfill leachate release by the grid-net electrical conductivity measurement system. The electrical conductivity of the subsurface before leachate injection lay within a small range of 24.8–43.0 S/cm. The electrical conductivity values in detected points were approximately ten times more than the conductivity values of the subsurface without landfill leachate intrusion. The results in this study indicate that the grid-net electrical conductivity measurement method has a possible application for detecting locations of landfill leachate intrusion into the subsurface at the initial stage, and thus has great potential in monitoring leachate leakage at waste landfills.     

裂谷作用对层序地层充填样式的控制——以西非裂谷系Termit盆地下白垩统为例

Termit盆地位于尼日尔东南部,属于西非裂谷系的北延部分,是发育于前寒武系—侏罗系基底之上的中、新生代裂谷盆地。该盆地早白垩世—古近纪经历了"裂谷—坳陷—裂谷"的构造演化过程及"陆相—海相—陆相"的沉积演化过程,表现为晚白垩世大规模海侵、早白垩世和古近纪两期裂谷叠置的特点。基于构造作用影响裂谷盆地层序发育的观点,分析了Termit盆地下白垩统裂谷阶段内的层序地层充填样式。根据裂谷作用的强弱,将早白垩世裂谷阶段划分为裂谷初始期、裂谷深陷期及裂谷萎缩期3个阶段。裂谷初始期层序断裂活动弱,构造沉降小,长轴物源体系较为发育,陡坡带为加积至退积型河流或三角洲沉积,缓坡带发育加积型河流或三角洲体系。裂谷深陷期层序断裂活动强烈,构造沉降大,陡坡带形成退积型水下扇或滑塌扇沉积,缓坡带发育退积型三角洲体系,盆地中心为泥岩充填。裂谷萎缩期层序断裂活动减弱并趋于停止,陡坡带为进积型扇三角洲沉积,缓坡带发育进积型三角洲体系。研究表明:裂谷作用对层序地层充填样式具有明显的控制作用,以构造作用为主线的裂谷盆地层序地层分析方法,能有效预测沉积体系和储层分布。     

含水层渗透系数空间变异性对地下水数值模拟的影响

地下水数值模拟是目前定量研究地下水水量和水质的重要手段。使用基于随机理论的MonteCarlo方法来进行地下水数值模拟。这种方法能较好地考虑水文地质参数的空间变异性。主要将MonteCarlo方法和确定性模型模拟方法的模拟结果在渗透系数场、水头场、速度场和浓度场等方面进行了比较。结果表明:在模拟三维非均质含水层中的溶质运移问题时,充分考虑了含水层渗透系数空间变异性的MonteCarlo法比确定性方法更为有效,模拟精度提高了很多,且对模拟误差及误差来源有合理的数学解释。     

Conceptual model of a coastal aquifer system in northern Greece and assessment of saline vulnerability due to seawater intrusion conditions

This paper refers to the development of a conceptual model for the management of a coastal aquifer in northern Greece. The research presents the interpretation and analysis of the quantitative (groundwater level recordings and design of piezometric maps) regime and the formation of the upcone within the area of investigation. Additionally it provides the elaboration of the results of chemical analyses of groundwater samples (physicochemical parameters, major chemical constituents and heavy metals and trace elements) of the area which were taken in three successive irrigation periods (July–August 2003, July–August 2004 and July 2005), in order to identify areas of aquifer vulnerability. The study identifies the areas where ion exchange phenomena occur, as well as the parts of the aquifer where the qualitative degradation of the aquifer system is enhanced. The paper, finally, assesses the lack of any scientific groundwater resources management of the area by the local water authorities, as well as the current practices of the existing pumping conditions scheme as applied by groundwater users.An erratum to this article can be found at     

土壤非饱和导水率模型中参数的敏感性分析

针对邵明安根据再分布过程得出的非饱和导水率模型,对其中参数的敏感性进行分析,通过实测资料计算非饱和导水率对不同参数的敏感度,对比参数对非饱和导水率的影响程度。研究表明采用线性关系表示土壤湿润剖面的平均湿度和湿润锋处湿度之间的关系时,非饱和导水率对参数的敏感性比其他参数高得多。     

Application of a solute transport model under variable velocity conditions in a conduit flow aquifer: Olalde karst system, Basque Country, Spain

A model based on numerical solutions, which allows for solving the dispersion equation under variable recharge and velocity conditions, is developed to simulate solute transport in conduit flow aquifers during flow recession periods. As an example, the evolution of a tracer in the little known karst conduit that links the sinking stream of Oma valley to the Olalde spring is investigated in the karstic region of Santa Eufemia-Ereñozar (Basque Country, Spain). The model, with different hypothetical structures, allows for obtaining series of tracer breakthrough curves, which are fitted to experimental data using an optimization algorithm. These results, although they can be used to simulate the tracer evolution between the two points considered, do not allow for determining the internal structure and spatial disposition of contributions in the aquifer.     

利用滨海地下水潮汐效应估计含水层参数:以北海半岛为例

通过分析滨海含水层中观测孔的水位动态资料,确定含水层的水力性质和某些参数,或推测含水层在海底的延伸长度和露头处是否存在覆盖层,是水文地质学家一直很感兴趣的问题之一。文中充分分析了广西北海半岛滨海地区水文地质条件,将研究区第一个承压含水层系统概化为具有上覆弱透水层越流和含水层海底露头覆盖层的垂向三层结构的不稳定流模型。应用该模型的解析解对广西北海半岛滨海含水层系统的观测数据进行了分析,拟合效果很好,并利用最小平方罚函数法估算了研究区含水层的参数。     

Estimating hydraulic conductivity using grain-size analyses, aquifer tests, and numerical modeling in a riverside alluvial system in South Korea

Hydraulic conductivity (K) for an alluvial system in a riverbank filtration area in Changwon City, South Korea, has been studied using grain-size distribution, pumping and slug tests, and numerical modeling. The alluvial system is composed of layers: upper fine sand, medium sand, lower fine sand, and a highly conductive sand/gravel layer at the base. The geometric mean of K for the sand/gravel layer (9.89?×?10^?4 m s^?1), as determined by grain-size analyses, was 3.33 times greater than the geometric mean obtained from pumping tests (2.97?×?10^?4 m s^?1). The geometric mean of K estimates obtained from slug tests (3.08?×?10^?6 m s^?1) was one to two orders of magnitude lower than that from pumping tests and grain-size analyses. K estimates derived from a numerical model were compared to those derived from the grain-size methods, slug tests and pumping tests in order to determine the degree of deviation from the numerical model. It is considered that the K estimates determined by the slug tests resemble the uppermost part of the alluvial deposit, whereas the K estimates obtained by grain-size analyses and pumping tests are similar to those from the numerical model for the sand/gravel layer of the riverside alluvial system.     

侧向有限越流承压含水层中非完整井非稳定流模型及解析解

工程建设中当距离抽水井r=rb处水位基本没有变化或不受抽水影响时,或当此处存在止水帷幕时,含水层系统视为侧向有限延伸,rb为有限半径。为此,构建更加符合工程实际的侧向有限延伸的典型弱透水层-承压水层系统中非完整井非稳定流计算模型,同时考虑井径和井储效应的影响,应用Laplace变换和分离变量法得到了水位降深在拉氏空间下的解析解,并应用拉氏数值逆变换Stehfest法得到真实空间下的水位降深。新建立的解析解可以进一步退化为诸多已有解,并进一步将其与已知解和有限元数值解进行对比,验证了所得解的正确性和可靠性。基于新建解重点分析了侧向边界和井的完整性对承压水层水位降深的影响。结果表明:含水层系统的侧向有限边界仅对抽水后期的水位降深影响明显,含水层系统侧向无限延伸情况下的水位降深要大于情形1(在r=rb处为定水头边界)且明显小于情形2(在r=rb处为不透水边界)下的水位降深,rb越小,两者之间的误差越大;抽水井的完整性对整个抽水期间不同情形下的水位降深均有明显的影响,承压含水层顶板处的水位降深随着抽水井滤管的长度和埋深的增加而减小。     

裂隙岩溶含水系统溢流泉演化过程的数值模拟

本文利用非连续裂隙网络介质模型,耦合水动力条件和碳酸盐溶蚀动力条件,建立岩溶含水系统演化的数学模型,并用数值法进行模拟分析。在已建数值模型基础上,构造了两组实例模型(降水入渗补给的裂隙岩溶含水系统和河流补给的裂隙岩溶含水系统)来研究裂隙岩溶系统中溢流泉早期演化过程。模拟发现,裂隙岩溶含水系统演化过程中,裂隙不断被溶蚀拓宽,致使系统中潜水位不断下降,许多小泉被疏干;同时不同裂隙之间存在溶蚀差异,在潜水位及优势通道附近的裂隙溶蚀速度快,而其它位置则相对慢得多,差异溶蚀促使系统形成优势溶管,优势溶管逐渐袭夺水流,形成大泉。另外对比两模型发现,降水入渗为主要补给源的裂隙岩溶含水系统,泉管道从汇向源发育;而河流入渗为主要补给源的裂隙岩溶含水系统,泉管道则由源向汇发育。      

越流含水层系统中的溶质运移方程

越流含水层系统在我国分布十广泛,但国内有关越流含水层系统的地下水污染模型尚见报导,为建立越流含水层系统地下水污染的数学模型,本文导出了包括越流和井流项的溶质运移方程。     

Evaluation of the ability of an artificial neural network model to simulate the input-output responses of a large karstic aquifer: the La Rochefoucauld aquifer (Charente,France)

The ability of artificial neural networks (ANN) to model the rainfall-discharge relationships of karstic aquifers has been studied in the La Rochefoucauld karst system, south-west France, which supplies water to the city of Angoulême. A neural networks model was developed based on MLP (multi-layer perceptron) networks and the Levenberg-Marquardt optimization algorithm. Raw rainfall data were used without transformation into effective rainfall. This allowed for the elimination of certain non-verifiable simplifying assumptions and their subsequent introduction into the modeling procedure. The raw rainfall and discharge data were divided into three groups for the training, the validation and the prediction test of the ANN model. The training and validation phases led to a very satisfactory calibration of the ANN model. The attempt to predict discharges showed that the ANN model is able to simulate the karstic aquifer discharges. The shape of the simulated hydrographs was found to be similar to that of the actual hydrographs. These encouraging results make it possible to consider interesting and new prospects for the modeling of karstic aquifers, which are highly non-linear systems.

---

Resumen Se ha estudiado la capacidad de las redes artificiales neurales (ANN) para modelizar las relaciones de lluvia-descarga de acuíferos kársticos en el sistema kárstico La Rochefocauld, al suroeste de Francia, el cual abastece de agua a la ciudad de Angoulême. Se desarrolló un modelo de redes neurales en base a redes MLP (Perceptron Multi-Capas) y el algoritmo de optimización Levenberg-Marquardt. Se utilizaron datos de lluvia sin la transformación hacia lluvia efectiva. Esto permitió la eliminación de ciertos supuestos simplificadores no verificables y su subsiguiente introducción en el procedimiento de modelizado. Los datos brutos de descarga y lluvia se dividieron en 3 grupos para la preparación, validación y la prueba de predicción del modelo ANN. Las fases de preparación y validación llevaron a una calibración muy satisfactoria del modelo ANN. El intento por predecir descargas mostró que el modelo ANN es capaz de simular las descargas del acuífero kárstico. Se encontró que la forma de los hidrogramas sintéticos es similar a la de los hidrogramas reales. Estos resultados alentadores hacen posible considerar prospectos nuevos e interesantes para el modelizado de acuíferos kársticos los cuales son sistemas altamente no-lineares.

---

Résumé L’aptitude des réseaux de neurones artificiels (RNA) à modéliser les relations pluie-débit des aquifères karstiques a été évaluée sur le karst de La Rochefoucauld (Sud-Ouest de la France), qui fournit l’eau potable à la capitale régionale Angoulême. Un modèle RNA a été développé à cet effet, basé sur les réseaux PMC (Perceptron Multicouche) et l’algorithme d’optimisation de Levenberg-Marquardt. Les données de pluie utilisées concernent la pluie brute, sans transformation en pluie efficace, ce qui permet de s’affranchir de certaines hypothèses simplificatrices non vérifiables. Les données de pluie brute et de débit ont été réparties en 3 groupes pour l’apprentissage, la validation et le test de prédiction du RNA. Les phases d’apprentissage et de validation ont permis d’aboutir à une calibration très satisfaisante du modèle RNA. La tentative de prédiction a montré que le RNA est apte à simuler les débits de l’aquifère karstique à partir de la pluie brute. La forme des hydrogrammes simulés est semblable à celle des hydrogrammes réels. Les résultats obtenus sont très encourageants et permettent d’envisager des perspectives intéressantes et nouvelles de modélisation des aquifères karstiques, qui sont des systèmes hautement non-linéaires.

     

A comparison of electrical and electromagnetic methods for the detection of hydraulic pathways in a fractured rock aquifer,Clare Valley,South Australia

Within fractured rock, the irregular and often unpredictable distribution and geometry of hydraulically conductive fractures produces large spatial variations in bore yield and groundwater quality. As fractures act as conduits for flow of both groundwater and electrical charge, methods which can efficiently detect the distribution of electrical pathways can be used to infer characteristics of significant hydrological parameters. This study compares the capabilities and limitations of electrical data obtained from direct current (DC) and electromagnetic (EM) surface azimuthal measurements, and from DC borehole-to-surface and cross-borehole measurements, for the interpretation of major hydrological structures in Clare Valley, South Australia. Electrical and EM surface methods are limited by poor depth sensitivity and the presence of conductive overburden, but provide useful tools for determining directional variations in resistivity at sites lacking bedrock exposure and boreholes. Application of borehole-to-surface methods yielded a better-resolved interpretation of sub-vertical fracture strike and was useful in identifying lateral variations in bedrock heterogeneity. Improved flexibility and sensitivity to measurements at depth permitted cross-borehole electrical tomography data to be used in reconstructing the spatial distribution of sub-horizontal, laterally extensive, electrically conductive zones. While the technique is restricted to small-scale sites with multiple boreholes, inferences can be made on fluid connections over a much larger regional scale. It is important to note, however, that while electrical methods provide valuable information about in-situ hydraulic pathways, they do not provide a complete hydraulic characterisation. Such a task requires integration of surface and borehole geophysics, geologic mapping, sampling and pumping tests of wells with packed-off intervals.

---

Resumen La geometría, y con frecuencia impredecible, distribución irregular de fracturas hidráulicamente conductivas dentro de roca fracturada genera variaciones espaciales grandes en producción de pozos y calidad de agua subterránea. Debido a que las fracturas actúan como conductos de flujo de agua y carga eléctrica, los métodos que detectan eficientemente la distribución de trayectorias eléctricas pueden utilizarse para inferir las características de parámetros hidrológicos significativos. Este estudio compara las capacidades y limitaciones de datos eléctricos obtenidos de mediciones azimutales superficiales electromagnéticas (EM) y de corriente directa (DC), y de mediciones de DC realizadas en la superficie y pozos así como mediciones realizadas entre pozos para la interpretación de estructuras hidrológicas principales en el Valle Clare, sur de Australia. Los métodos superficiales eléctricos y EM están limitados por sensitividad de profundidad pobre y la presencia de cubierta conductiva, pero aportan herramientas útiles para determinar variaciones direccionales de resistividad en sitios que carecen de pozos y afloramientos rocosos. La aplicación de métodos superficiales y de pozos aportan una mejor interpretación del rumbo de fracturas sub-verticales y fue útil en identificar variaciones laterales en la heterogeneidad del macizo rocoso. El mejoramiento de flexibilidad y sensitividad en las mediciones profundas permitió que los datos de tomografía eléctrica de los pozos fuera utilizado en la reconstrucción de la distribución espacial de zonas eléctricas conductivas, sub-horizontales y lateralmente extensas. Aunque la técnica está restringida a sitios de pequeña escala con múltiples pozos, puede realizarse inferencias sobre relaciones entre fluidos en una escala regional mucho más grande. Sin embargo, es importante notar que aunque los métodos eléctricos aportan información valiosa acerca de las trayectorias hidráulicas in-situ, aún no proporcionan una caracterización hidráulica completa. Esta tarea requiere integrar geofísica superficial y de pozos, mapeo geológico, muestreo y pruebas de bombeo en pozos con intervalos sin empaque.

---

Résumé Au sein des roches fracturées, lirrégularité et limprédictible distribution et géométrie des fractures par lesquelles sécoulent leau souterraine produit de larges variations spatiales entre les débits des forages et les paramètres de la qualité de leau. Comme les fractures conduisent et leau et les charges électriques, les méthodes qui peuvent de manière efficiente détecter la distribution des courants électriques peuvent être utilisées pour analyser les caractéristiques des principaux paramètres hydrologiques. Cette étude compare les capacités et les limites des données produites par Courant Direct (DC en Anglais) et par mesure Electromagnétique azimutale de surface (EM en Anglais), et par DC entre surface et forage et entre plusieurs forages, pour linterprétation de la structure hydrologique de la Vallée de Clare, Australie du Sud. Les méthodes électriques et électromagnétiques sont limitées par la faible sensitivité à la profondeur et la présence dune couverture conductrice, mais procure des outils utiles pour déterminer les variations directionnelles en terme de résistivité là où la roche nest pas affleurante et dans les puits. Lapplication de la méthode dinvestigation «forage à surface» apporte une meilleure interprétation des fractures sub-verticales et des hétérogénéités latérales. Les tomographies électriques entre forage bénéficient de la flexibilité et de la sensibilité des mesures en profondeur, et permettent de délimiter lextension latérale des hétérogénéités sub-horizontales de zones conductrices. Alors que la technique est restreinte à de petits sites comprenant de nombreux forages, il est possible de reconstituer les connections hydrauliques à des échelles régionales. Il est important de noter que les méthodes fournissent des informations intéressantes mais pas des caractérisations hydrauliques complètes. Pour cela les données pourraient être complétées par des études plus poussées intégrant les différentes prospections géophysique, les données des cartes géologiques, des échantillonnages et des essais de pompage à différents intervalles de profondeur.

     

多孔介质非饱和导水率预测的分形模型

多孔介质非饱和导水率是地下水污染预测与评价的重要参数。根据分形几何的基本原理和方法,推导出了与Campbell经验公式在形式上完全一致的多孔介质非饱和导水率的预测公式。公式中的幂指数为介质孔隙分维和随机行走分维的函数,分别体现了多孔介质的静态性质与动态性质对其中水分运动的影响,但静态性质的影响是主要的,即导水率主要受多孔介质的结构控制。根据文献中报道的大量数据,利用笔者推导的预测公式计算得到的幂指数的统计值与试验测定的幂指数的统计值基本一致,说明推导的理论公式预测多孔介质非饱和导水率是较为可靠的。     

地下含水层储能两阶段热量运移数值模型研究

为了简单而准确地预测地下含水层储能情况,考虑了地下含水层储能过程时间跨度大、储能保温分阶段的实际特点,分析了地下含水层流动和换热模型物理参数的对比特点,以热平衡和热扩散原理为基础,建立分成两个连续阶段的地下含水层储能数值计算模型.模型求解中,采用控制容积法,以全隐格式进行热扩散方程的离散化,然后应用Jacobi方法迭代求解,模拟结果和实际观测数据吻合很好.该模型还分析了含水层储能循环采灌过程中抽出储能水的温度变化的一般特点.