the changing relationship of hydrogeological parameters of dadianzi well-aquifer system

Hydrogeological parameter is an important index to characterize the hydrogeological properties of the aquifer, and has a clear physical basis and mechanism. Although the predecessors have made significant achievements in these areas, research is lacking on the changing law and relationship of the hydrogeological parameters of well-aquifer system. The digital water level and barometric pressure data of Dadianzi Well are used as the basis in this study. Based on the theories of elastic mechanics, rock mechanics and fluid mechanics, and using barometric pressure coefficient and tidal factor, the hydrogeological parameters in Dadianzi well-aquifer system in undrained conditions are studied. The corresponding water storage rate can also be obtained quantitatively. In addition, with the thickness of the aquifer, the pressure transmitting coefficient, the radius of the well and the frequency of the tidal wave, the permeability coefficient and transmissibility coefficient of well-aquifer system can be obtained, and the relationships between them are derived. The results show that: 1)There is an obvious power function relationship between porosity and solid skeleton volume compression coefficient, volume compression coefficient of water in aquifer, water storage rate, permeability coefficient and transmissibility coefficient. The volume compression coefficient of solid skeleton, water storage rate, permeability coefficient and transmissibility coefficient have a positive correlation with the porosity, the volume compression coefficient of water in aquifer decreases with increasing porosity. The volume compression coefficient of solid skeleton and water in aquifer can be well fitted to one of two quadratic polynomials. And the volume compression coefficient of water in aquifer is larger than the solid skeleton volume compression coefficient, water is more easily compressed. In addition, with the increase of water storage rate, the permeability coefficient and transmissibility coefficient also increase linearly; 2)Different from the traditional pumping test and indoor experiment, this paper uses the digital water level and other data, combined with the pressure coefficient and Venedikov tidal harmonic analysis results to access to the porosity, the volume compression coefficient of solid skeleton and water in aquifer medium, water storage rate, the permeability coefficient and the transmissibility coefficient. This method is simple and accurate.     

井-含水层系统地下水类型定量分析和判别

利用水位、 气压和理论固体潮数据, 采用卷积回归法中水位对气压的阶跃响应函数, 定量地分析和判定了华北北部板桥井、 大灰厂井、 黄骅井、 大甸子井、 丰镇井和三号地井的井-含水层系统的地下水类型, 并结合研究时段内各井的气压系数和M₂波潮汐因子的结果进行了对比分析． 结果表明: ① 各井的滞后时间与阶跃响应函数之间存在明显的以e为底的指数函数关系, 且底数e的系数的正负决定了井-含水层系统的地下水类型; ② 承压井的阶跃响应函数随滞后时间的增大而增大, 且最佳阶跃响应函数值越大, 相应的气压系数和M₂波潮汐因子也越大, 反之亦然; ③ 潜水井和半承压水井的阶跃响应函数随滞后时间的增大而减小, 其最佳阶跃响应函数与气压系数和M₂波潮汐因子间的关系不明显, 可能与含水层的水力特性、 井孔结构及固体潮汐波的频率有关．      

Estimation of coseismic deformation, poroelasticity, and fracturing of rocks from the data on water level in a borehole

Variations in the water level in boreholes emerge in response to tidal, baric, and tectonic forcing. We analyze the data on atmospheric pressure and water level recorded in the boreholes located in the mid-latitude Eurasia (45°?C55°N) from Western Europe (Belgium, Uccle), Siberia (coastal area of Lake Baikal, Talaya River) to Far Eastern Russia (the Bychikha borehole near Khabarovsk and the boreholes on the Kamchatka Peninsula and Kurils). The response of the water level in a borehole to periodic tidal and baric impacts is investigated. In this case, water level variations reflect areal (lateral) and vertical deformations, which allows estimating the elastic moduli and porosity of the confined aquifer in a static model. Measurements in the boreholes drilled in the fractured solid rock enable determining the geometric parameters of the fractures. The possibility to apply the method for evaluating tectonic deformations is discussed. Application of tidal coefficients of boreholes for determining the coseismic deformation is demonstrated by the example of the Kultuk earthquake (Lake Baikal, August 27, 2008, M = 6.3).     

Application of Wavelet Decomposition to Removing Barometric and Borehole Water Tidal Response in Level

Wavelet decomposition is used to analyze barometric fluctuation and earth tidal response in borehole water level changes. We apply wavelet analysis method to the decomposition of barometric fluctuation and earth tidal response into several temporal series in different frequency ranges. Barometric and tidal coefficients in different frequency ranges are computed with least squares method to remove barometric and tidal response. Comparing this method with general linear regression analysis method, we find wavelet analysis method can efficiently remove barometric and earth tidal response in borehole water level. Wavelet analysis method is based on wave theory and vibration theories. It not only considers the frequency characteristic of the observed data but also the temporal characteristic, and it can get barometric and tidal coefficients in different frequency ranges. This method has definite physical meaning.     

小波分析在井水位的气压和潮汐改正中的应用

本文将小波分析方法应用于地下水位观测资料的气压和潮汐改正中。该方法利用小波将气压观测和理论重力固体潮资料分解为不同频段的时间序列,用最小二乘法求出气压和固体潮不同频段的响应系数,以消除气压和固体潮对水位的影响。用该方法计算的结果与一般线形回归分析方法得到的结果进行对比,发现该方法能更有效地分析和消除气压和固体潮对地下水位的影响。小波分解在井水位气压、潮汐改正中的应用以波动理论和振动理论为基础,它不仅考虑了观测资料的频率域特性,而且还能直接在时间域内计算出不同频率范围的气压系数和潮汐系数,具有明确的物理意义。     

Application of Wavelet Decomposition to Removing Barometric and Tidal Response in Borehole Water Level

Wavelet decomposition is used to analyze barometric fluctuation and earth tidal response in borehole water level changes. We apply wavelet analysis method to the decomposition of barometric fluctuation and earth tidal response into several temporal series in different frequency ranges. Barometric and tidal coefficients in different frequency ranges are computed with least squares method to remove barometric and tidal response. Comparing this method with general linear regression analysis method, we find wavelet analysis method can efficiently remove barometric and earth tidal response in borehole water level. Wavelet analysis method is based on wave theory and vibration theories. It not only considers the frequency characteristic of the observed data but also the temporal characteristic, and it can get barometric and tidal coefficients in different frequency ranges. This method has definite physical meaning.     

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

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

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

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

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

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

Storage in confined aquifer: Spectral analysis of groundwater in responses to Earth tides and barometric effect

Dilatation of aquifer and associated water level fluctuation in groundwater well is known to be driven periodically from lunar, solar, or other tidal forces. Time‐dependent variables in groundwater system, such as water level, can be converted to power spectra in the frequency domain using Fourier transform to evaluate significant fluctuation. The major innovation of this research is to develop spectral representation in frequency domain for the groundwater system that the storage in confined aquifer can be determined considering dilatation affected by Earth tides and barometric effect. In order to verify applicability of the evolved method, time series of Earth tides and barograph are collected; aquifer storage is then determined inversely by selecting significant semidiurnal and diurnal components in spectra computation. It suggests that to discover groundwater storage using groundwater level with barograph and tidal potential of Earth in frequency domain becomes accessible and feasible.     

气压变化及其对地壳形变和深井水位的影响

根据北京塔院地区1984和1985年的微气压计记录采用频谱分析和调和分析两种方法得出了大气潮谱的详细结构;用负荷勒甫数代替潮汐勒甫数,仿照固体潮的有关公式,导出了大气潮引起的地壳应变和地倾斜公式;使用北京塔院井的水位观测资料结合当地的气压资料,分析了气压变化对深井水位的影响.理论计算和实测资料的分析结果基本相符,特别是两种结果都得出S₂大气潮引起的体膨胀约为S₂固体潮体膨胀的20％.     

井水位的“记忆”滞后效应

观测资料表明,井水位对信息响应存在的“记忆”滞后现象,它与一般的位相滞后不同,在鲁２９井现场试验也证明了井水位对井吕空气压力变化的响应存在的“记忆”滞后现象。利用水平层状承压含水层模式,从理论上解释了井水位对井口空气压力变化响应“记忆”滞后现象,认为这种现象与水井含水层的导水系数有关,含水层导不系数越小,这种现象越明显,用一般多元回归方法无法较好地扣除井水位中“记忆”滞后影响,作给出了一种可以扣     

井水位随气压变化之解析解

气压作用下, 井-含水层系统中地下水流是一类流体力学问题。 本文应用井壁水流通量边界条件和气压作用下井壁内外水(孔)压平衡条件, 提出了一个井水位随气压变化的解析公式。 解析式表明, 气压系数随时段长度增加而增大, 并趋于气压常数; 气压系数随时段长度的变化只依赖于导水系数与井半径平方的比值(T/r²_w), 而与气压变化过程无关; 气压常数只与含水层的一维荷载效率(B)有关, 而与导水系数和井半径无关。 解析解所反映的气压系数与时段长度的关系, 与南溪井实测序列数据分析结果具有很好的一致性。 根据气压系数随时段长度变化过程, 提出了一个参数估计方法, 应用于估计南溪井含水层气压常数和导水系数, 并对本文提出的参数估计方法进行了讨论。     

气压变化及其对地壳形变和深井水位的影响

根据北京塔院地区1984和1985年的微气压计记录采用频谱分析和调和分析两种方法得出了大气潮谱的详细结构;用负荷勒甫数代替潮汐勒甫数,仿照固体潮的有关公式,导出了大气潮引起的地壳应变和地倾斜公式;使用北京塔院井的水位观测资料结合当地的气压资料,分析了气压变化对深井水位的影响.理论计算和实测资料的分析结果基本相符,特别是两种结果都得出S_2大气潮引起的体膨胀约为S_2固体潮体膨胀的20％.     

Multi‐depth monitoring of electrical conductivity and temperature of groundwater at a multilayered coastal aquifer: Jeju Island,Korea

To supplement conventional geophysical log data, this study presents temporal variations in electrical conductivity (EC) and temperature with depth in a multilayered coastal aquifer, on the eastern part of Jeju Island, Korea. One‐month time‐series data obtained at eight points from a multi‐depth monitoring system showed that semidiurnal and semimonthly tidal variations induced dynamic fluctuations in EC and temperature. At some depths, EC ranged from 1483 to 26 822 µS cm^?1, while some points showed no significant variations. The results of EC log and time‐series data revealed that a sharp fresh‐saltwater interface occurred at low tide, but the diffusion zone broadened to 20 m at high tide. EC, temperature, and tide level data were used for the cross‐correlation analysis. The response time of EC and temperature to tide appears to range from less than 30 min to 11 h. Using end‐member mixing analysis (EMMA), the fraction of variations of chloride concentration in the multilayered aquifer was explained, and a conceptual model was developed which subdivided the coastal aquifer into four vertical zones. The percentage of water derived from seawater varied from 2 to 48 at specific depth, owing to tidal fluctuations. Continuous observations of EC and temperature at multiple depths are powerful tools for quantifying the transport of saline water by tidal variations in multilayered coastal aquifers. Copyright © 2008 John Wiley & Sons, Ltd.     

基于井潮、气压效应的苏皖地区突出水位异常分析

基于流体及岩石力学理论,利用气压和潮汐效应,研究了苏皖地区苏02井、苏03井、苏06井、苏18井和定远04井含水层介质在不排水状态下孔隙度和渗透系数的变化情况.结合井-含水层系统在受力和水均衡两种状态下的微动态变化模型,分析了各井水位自2014年下半年以来相继出现同步上升现象的水动力学成因机制.结果表明:① 5口井的气压系数与M₂波潮汐因子之间存在明显的线性关系,即气压系数越大,潮汐因子也越大;②各井的孔隙度和渗透系数出现一组同步增大现象,主要与补给量大于泄流量、补给增强有关.      

A semi‐analytical solution for groundwater responses to stream‐stage variations and tidal fluctuations in a coastal aquifer

A two‐dimensional semi‐analytical solution to analyse stream–aquifer interactions in a coastal aquifer where groundwater level responds to tidal effects is presented. The conceptual model considered is a two‐dimensional subsurface system with stream and coastline boundaries at right angles. The dimensional and non‐dimensional boundary value problems were solved for water level in the aquifer by successive application of Laplace and Fourier transform techniques, and the results were obtained by numerical inversion of the transformed solution. The solution was then verified by reducing the solutions to one‐dimensional known problems and comparing the results with those from previous studies. Hypothetical examples were used to examine the characteristics of water‐level variations due to the variations in stream stage and the fluctuations in tide level. Sensitivity analysis indicated that streambed leakance has no influence over the amplitude of groundwater fluctuations, but that the effect of stream stage increases with increasing leakance. Little difference was observed in the water level for different aquifer penetration ratios with narrow stream width. Increases in streambed leakance caused increases in the effect of aquifer penetration by the stream on the water level. An increased specific yield value resulted in decreased amplitude of water fluctuations and mean water level, and showed that water‐level variations due to stream and tidal boundaries are sensitive to specific yield. Copyright © 2006 John Wiley & Sons, Ltd.     

Teleconnection and climatic oscillation in aquifer water level in Kumamoto plain,Japan

Relative little is known about the interaction between climate change and groundwater. Analysis of aquifer response to climatic variability could improve the knowledge related to groundwater resource variations and therefore provides guidance on water resource management. In this work, seasonal and annual variations of groundwater levels in Kumamoto plain (Japan) and their possible interactions with climatic indices and El Niño Southern Oscillation (ENSO) were analyzed statistically. Results show the following: (1) The water level in the recharge area mainly fluctuates at 1‐ and 2‐year periods, whereas the significant periodicity for water level oscillation in the coastal aquifer is 0.5 year. (2) The aquifer water levels are possibly influenced by variability in precipitation, air temperature, barometric pressure, humidity variances and ENSO. Relative high correlations and large proportions of similarities in wavelet power patterns were found between these variables and water levels. (3) Aquifer response to climatic variances was evaluated using cross wavelet transform and wavelet coherence. In recharging aquifers, the ENSO‐induced annual variations in precipitation, air temperature, humidity and barometric pressure affect aquifer water levels. The precipitation, air temperature and humidity respond to ENSO with a 4‐, 6‐ and 8‐month time lag, respectively, whereas the ENSO imparts weak influence on the barometric pressure. Significant biennial variation of water levels during 1991–1995 is caused primarily by precipitation and humidity variations. In the coastal aquifer, the 0.5‐year variability in ENSO is transferred by precipitation, barometric pressure and humidity to aquifer water levels, and the precipitation/humidity influence is more significant comparing with the barometric pressure. Copyright © 2014 John Wiley & Sons, Ltd.     

深井水位对地震波、固体潮和气压的水力响应——以范县井为例

以河南范县井为例,利用不同的水力响应模型分析了井水位对地震波、固体潮和气压的响应特征,并基于相关水力响应模型反演估算了含水层的水力参数。结果显示:在高频加载作用过程中,井-含水层系统中的水流模式以水平向为主,而在低频加载作用过程中,则为水平向和垂直向共存的混合模式;利用周期为10—102 s的高频段的地震波响应模型估算的含水层导水系数值较大,为7.20×10-3 m2/s,利用周期为3.75×104 s的低频段的固体潮响应模型估算的含水层导水系数值较小,为2.02×10-6 m2/s,而利用周期为102—104 s的中等频率段的气压响应模型得到的估算值介于二者之间,为3.44×10-5 m2/s。由此分析认为,在周期性加载作用过程中,井-含水层系统内的水流模式与加载频率有关,基于不同水力响应模型反演估算的含水层水力参数存在尺度效应。本研究取得的认识,既可为井水位动态响应的机理解释提供理论基础,也可为目标含水层水力参数的原位测量提供技术支撑。     

Tidal and non-tidal signals in groundwater boreholes in the KTB area,Germany

Tides and barometric pressure variations cause pore pressure changes in the solid earth. In boreholes which are hydraulically connected to confined aquifers these pore pressure changes can be observed as water level variations. In case of confined aquifers boreholes can be regarded as volumetric strainmeters. From June 2004 until May 2005 a large scale injection experiment was realised in the pilot borehole (4000 m) at the German Continental Deep Drilling (KTB) site in Bavaria. During this test the injection induced deformation was observed by a tiltmeter array around KTB consisting of five tiltmeter stations at horizontal distances from 1.6 to 3 km from the injection borehole. At each station the local groundwater level was observed in depths between 30 and 50 m. The pore pressure was recorded in the boreholes. The time series were checked for tidal signals and injection influence. The injection was not verified presumably due to the fact that the penetration of the injected water was controlled by local geology and tectonic faults zones. Although the boreholes extending only to shallow depth tidal signals are clearly observable in at least two of five stations. On December the 26th 2004 the Sumatra–Andaman earthquake occurred which is clearly detectable in the pore pressure variations at all five stations.