双重介质含水层中井水位潮汐相位滞后现象的理论分析

简略地介绍了双重介质含水层模型；推导出引潮力作用下双重介质含水层中渗流运动的偏微分方程；将此泛定的方程与潜水的Boulton方程进行数学类比，从而得出一个新的迟后给水项，为解释固体潮致井水位振荡相位滞后现象奠定了物理基础；通过对该泛定方程的分析，还发现该迟后给水项与深部含水层中存在的一种潮汐液流振荡现象有关．      

地下水潮汐现象的物理机制和统一数学方程

用新的分层承压含水层模式 ,不但考虑含水层的力学压缩性质 ,而且考虑含水层的渗流特性 ,并结合扰动信息源的频率特性 ,分别研究扰动源地球固体潮、大气潮和地表负荷潮对承压井水位和流量的影响机理 ,给出相应的偏微分方程。从方程的解释或数值解讨论扰动源与承压井含水层的力学压缩参数、渗流特性参数及与频率特性频数的关系 ,进而给出承压井水位和流量对地球固体潮、大气潮和地表负荷潮汐响应的统一数学方程及其潮汐响应函数 ,并揭示了上述几类潮汐扰动信息源对承压井水位和流量影响的物理机理     

中国地下水潮汐的观测研究和分析

本文简单介绍了中国地下水潮汐的观测研究情况,介绍了对承压水位潮汐机理研究的理论结果,给出了固体潮体应变、气压和海洋潮汐对承压水头影响的偏微分方程,从这些方程的简化解讨论了它们与含水层参数间的关系.并从全国选取了33口有代表性的井孔,对这些井孔的水位和气压观测资料采取反复调和分析的方法,求出了它们的水位固体潮体应变各分波的潮汐系数和位相滞后,也求出了水位气压各分波的系数和位相滞后.      

深井水位的固体潮效应

本文从体应变固体潮对深井水位影响的偏微分方程出发,考虑到含水层和井孔之间相互渗流的边界条件,用叠加原理、冲量定理和分离变量法等方法得出了方程的解.通过对这个解中水井含水层参数给予一些可能的值进行数值计算,讨论了水井固体潮系数和位相滞后与水井含水层参数间的关系,较好地解释了井水位对固体潮响应的位相滞后现象.计算表明,井孔的半径、含水层的孔隙度及固体骨架的体压缩系数愈大,含水的导水系数愈小,则水井水位的固体潮系数愈小,而水位对固体潮响应的位相滞后愈大.井水对长周期的潮汐响应比对短周期的更好.      

Beitrag zur analytischen Beschreibung von Transportprozessen im Grundwasserleiter bei diffusem Stoffeintrag

A model is described for the relation between water quantity and water quality in the aquifer with the diffuse substance input by the gravitation water being taken into account. The model is based on a partial differential equation of the second order for a flow pipe. For constant transport parameters, a sink of the kinetics of the first order and a time-dependent source term an exact analytical solution is presented and explained by a model. At a diffuse substance input, the dispersion has only a slight influence on the transport of substances contained in water. Potential applications of these solutions to different problems are mentioned.     

线弹性含水层井水位、孔压对引潮位响应的研究及其应用

应用线弹性介质力均衡方程,研究不排水条件下饱水岩体在潮汐力作用下的体应变,提出了承压含水层孔压对引潮高的线性响应方程,并给出了该方程响应系数(E)的物理意义.结合Hsieh 等(1987)提出的井水位对孔压响应的振幅比(A)和位相差(α1)公式,进一步推导出了井水位-引潮高振幅比M=EA和位相差α=α1+α2公式.M和...     

A double-porosity model for a fractured aquifer with non-Darcian flow in fractures

Abstract

Non-Darcian flow in a finite fractured confined aquifer is studied. A stream bounds the aquifer at one side and an impervious stratum at the other. The aquifer consists of fractures capable of transmitting water rapidly, and porous blocks which mainly store water. Unsteady flow in the aquifer due to a sudden rise in the stream level is analysed by the double-porosity conceptual model. Governing equations for the flow in fractures and blocks are developed using the continuity equation. The fluid velocity in fractures is often too high for the linear Darcian flow so that the governing equation for fracture flow is modified by Forcheimer's equation, which incorporates a nonlinear term. Governing equations are coupled by an interaction term that controls the quasi-steady-state fracture—block interflow. Governing equations are solved numerically by the Crank-Nicolson implicit scheme. The numerical results are compared to the analytical results for the same problem which assumes Darcian flow in both fractures and blocks. Numerical and analytical solutions give the same results when the Reynolds number is less than 0.1. The effect of nonlinearity on the flow appears when the Reynolds number is greater than 0.1. The higher the rate of flow from the stream to the aquifer, the higher the degree of nonlinearity. The effect of aquifer parameters on the flow is also investigated. The proposed model and its numerical solution provide a useful application of nonlinear flow models to fractured aquifers. It is possible to extend the model to different types of aquifer, as well as boundary conditions at the stream side. Time-dependent flow rates in the analysis of recession hydrographs could also be evaluated by this model.     

深井水位对固体潮和气压的响应

本文以体应变固体潮对深井水位影响的偏微分方程为基础，考虑到含水层与井孔之间相互渗流的边界条件，用叠加原理、冲量定理的分离变量等方法得出了方程的解。把水井含水层的参量给予一些可能的值，通过数值计算讨论了水井固体潮系数、位相滞后和含水层参数之间的关系，解释了井水位对固体潮响应的位相滞后现象。结果表明，井孔的半径、含水层的孔隙度及固体骨架的体压缩系数愈大，含水层的导水系数愈小，则水井水位的固体潮系数愈小     

Effect of earth tide on deep well water level

On the basis of the partial differential equations of the water level in a deep well that is affected by the bulk strain tide, the boundary conditions of permeating each other between aquifer and well are considered. The solutions of the equations have been obtained by the superposition principle, inpulse theorem and separation of the variations. Some suitable values are given to the parameters of well aquifer in the solutions. By means of the numerical calculation the relations among the well tide coefficients and phase lag phenomenon of response of the well level to earth tides is well interpreted. The parametric studies and the obtained graphs show that the more the radius of casing in a well, the porosity of an aquifer and the bulk compressibility for solid matrix of an aquifer, the less the transmissibility of an aquifer, then the less the tidal coefficient of well level, while the more the phase lag for well level to respond to the tide is. The response of the well level to the tide of longer period is better than the shorter period events. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 66–75, 1991.     

A quasi three-dimensional model for flow and transport in unsaturated and saturated zones: 1. Implementation of the quasi two-dimensional case

A quasi three-dimensional (QUASI 3-D) model is presented for simulating the subsurface water flow and solute transport in the unsaturated and in the saturated zones of soil. The model is based on the assumptions of vertical flow in the unsaturated zone and essentially horizontal groundwater flow. The 1-D Richards equation for the unsaturated zone is coupled at the phreatic surface with the 2-D flow equation for the saturated zone. The latter was obtained by averaging 3-D flow equation in the saturated zone over the aquifer thickness. Unlike the Boussinesq equation for a leaky-phreatic aquifer, the developed model does not contain a storage term with specific yield and a source term for natural replenishment. Instead it includes a water flux term at the phreatic surface through which the Richards equation is linked with the groundwater flow equation. The vertical water flux in the saturated zone is evaluated on the basis of the fluid mass balance equation while the horizontal fluxes, in that equation, are prescribed by Darcy law. A 3-D transport equation is used to simulate the solute migration. A numerical algorithm to solve the problem for the general quasi 3-D case was developed. The developed methodology was exemplified for the quasi 2-D cross-sectional case (QUASI2D). Simulations for three synthetic problems demonstrate good agreement between the results obtained by QUASI2D and two fully 2-D flow and transport codes (SUTRA and 2DSOIL). Yet, simulations with the QUASI2D code were several times faster than those by the SUTRA and the 2DSOIL codes.     

Hydrodynamic characterization of a Sahelian coastal aquifer using the ocean tide effect (Dridrate Aquifer,Morocco)

Abstract

In wells tapping coastal aquifers, piezometric fluctuations can be observed in response to the ocean tide. Simultaneous recordings of the ocean tide and of the piezometric variations may provide a basis for characterizing the hydrodynamics of the aquifer. This approach was attempted to characterize the Dridrate aquifer, located on the Atlantic coast of Morocco. This aquifer accounts for most of the regional drinking water resources. However, its hydrodynamic characteristics are very poorly known. The study compares observed and simulated piezometric fluctuations, under various assumptions (confined, semi-confined aquifer). The model, which best explains the hydrodynamic behaviour of this aquifer is a semi-confined and strongly heterogeneous aquifer model (calculated hydraulic diffusivity values vary over several orders of magnitude). This result is new and rather surprising, since to date this aquifer was considered confined in view of its geological setting. Consequently, new questions are raised regarding the protection and management of the groundwater resources of this aquifer.     

承压井水位对气压动态过程的响应

本文从气压影响承压井水位的偏微分方程出发,考虑到井孔和含水层之间相互渗流的边界条件,得出了方程的解。给出了承压井水位对气庄动态过程响应的气压效率和位相滞后。讨论了气压效率和位相滞后与水井含水层系统各参数的关系。并指出气压效率越大的承压井,其含水层对应力应变的敏感程度越差。     

冷水源井水温固体潮现象及其机理研究以黄村井为例

该文对黄村冷水源井的水温固体潮现象进行了观测研究,初步解释了冷水源井的水温固体潮汐的形成机理.通过对黄村观测井进行水温梯度的详细观测及不同深度水温的对比观测研究,得到了黄村井水温潮汐现象的观测结果:① 黄村冷水源井的水温固体潮的相位与水位的相位是相反的;② 黄村井水温梯度曲线呈下凹型,特别是在含水层及受含水层进出水影响较大的附近区域下凹程度大,且随着与含水层底板距离逐渐变大,下凹程度也变小;③ 水温传感器的观测值与含水层观测距离存在一定的规律性:距离含水层越远,水温潮汐差越小, 直至潮汐变化消失.这说明冷水源井与热水井的水温潮汐现象是不同的,前者是吸热过程, 后者是放热过程,由此造成二者水-热动力学特征的不同.      

日本秋田地震(1983年7.7级)引起我国应力场的震时调整变化

本文通过水平分层承压含水层模型,推导出含水层的水头和含水层应力变化之间的定量关系.利用我国在日本秋田地震时有水位阶变的52口井孔的观测资料,通过水井水位的体应变固体潮系数,估算了震时引起我国应力场的调整变化.讨论了这种应力调整时应力增加区和降低区交替出现的可能原因.      

A new theoretical solution of the effect of atmospheric pressure on water level

Based on the partial differential equation governing the effect of atmospheric pressure on water level of confined well, deriving the boundary condition and considering the seepage water between well and aquifer, the author obtained the analytical solution of water level change in time domain under the action of an atmospheric pressure history with the Laplace transform method. This solution is composed of two terms:stable and retarded terms. The stable term is the multiplication of barometric efficiency and simultaneous atmospheric pressure, and it implies the value of water level after infinite time when the atmospheric pressure is a constant from the time in question. The retarded term is the transient process due to the time lag of water exchange between well and aquifer. From the solution, it is obtained that the interference of atmospheric pressure on water level is the integral superimposition of the contribution of all atmospheric pressure changes before the time in question. So that, we further found out the response function of pulsive atmospheric pressure history. Calculation shows: (1) The pulsive response function starts from zero and tends to a steady value, which is proportional to the barometric efficiency, when the time tends to infinity; (2) The retarded time depends on the mechanical property of aquifer and the radius of well. The larger the seepage coefficient, the smaller the radius of well and the thicker the aquifer, then the shorter the retarded time gets. This solution can be used as the theoretical basis for further analysis of the atmospheric effect and practical correcting method in the future.     

自流井水位响应能力与水柱高度关系的分析

分析了万全等４口自流井水位阶变、固体潮幅度与水柱高度的关系。结果表明，当含水层所受应力状态不变，水位观测条件相同时，水位阶变、固体潮幅度度与水柱高度有正相关系。     

华蓥山地区4口井水位潮汐动态特征研究

本文选择沿华蓥山断裂带分布的荣昌等4口观测井,利用Baytap-G潮汐分析方法,计算各井水位和气压及理论固体潮的潮汐振幅谱,比较其潮汐频谱差异,通过对主要潮汐分波振幅的回归计算定量分析各井水位受气压潮和固体潮影响的大小。基于对井水位正常动态的认识,选择各井水位潮汐的主要分波,对井水位长时序数据进行分析计算,提取水位潮汐响应特征参数(振幅比和相位差),进而探讨特征参数动态变化特征。最后对井水位受气压潮和固体潮影响的差异原因进行了初步探讨。结果表明,荣昌井水位主要受气压作用的影响,北碚、大足、南溪三口井水位受固体潮-气压潮综合作用的影响,而荣昌井水位只受气压潮影响可能与该井所处含水层裂隙发育且该井未下设止水套管有关;荣昌井P_1S_1K_1波和南溪井M_2波振幅比和相位差在几次大震后没有明显变化,说明地震波没有使井孔与含水层之间的水流交换发生显著变化,而北碚井和大足井M_2波振幅比和相位差分别在汶川和芦山地震时发生变化,反映了地震波的疏通影响。     

井水位的固体潮效应和气压效应与含水层参数间的定量关系

本文分析了水井水位的固体潮效应和气压效应与承压含水层参数之间的定量关系。结果表明,含水层的固体骨架体积压缩系数和孔隙度越小,水井水位的固体潮系数则越大;含水层的固体骨架的体积压缩系数越小,孔隙度越大,則水井水位的气压效率越大。本文还把理论结果和固体潮及气压的实际观测结果进行了对比,两者符合得较好。     

Biased monitoring of fresh water-salt water mixing zone in coastal aquifers

In coastal aquifers, significant vertical hydraulic gradients are formed where fresh water and underlying salt water discharge together upward to the seafloor. Monitoring boreholes may act as "short circuits" along these vertical gradients, connecting between the higher and the lower hydraulic head zones. When a sea tide is introduced, the fluctuations of both the water table and the depth of the mixing zone are also biased due to this effect. This problem is intensified in places of long-screen monitoring boreholes, which are common in many places in the world. For example, all approximately 500 boreholes of the fresh water-salt water mixing zone in the coastal aquifer of Israel are installed with 10 to 50 m long screens. We present field measurements of these fluctuations, along with a three-dimensional numerical model. We find that the in-well fluctuation magnitude of the mixing zone is an order of magnitude larger than that in the porous media of the actual aquifer. The primary parameters that affect the magnitude of this bias are the anisotropy of the aquifer conductivity and the borehole hydraulic parameters. With no sea tide, borehole interference is higher for the anisotropic case because the vertical hydraulic gradients are high. When tides are introduced, the amplitude of the mixing zone fluctuation is higher for the isotropic case because the overall effective hydraulic conductivity is greater than the conductivity in the anisotropic case. In the aquifer, the fresh water-salt water mixing zone fluctuations are dampened, and tens of meters inland from the shoreline, the fluctuations are on the order of few centimeters.