Analytical solution for enhanced recharge around a bedrock exposure caused by deep-aquifer dewatering through a variable thickness aquitard

In this study an analytical solution was developed to predict steady radially-symmetric percolation rates from an aquifer underlain by a variable thickness aquitard. The solutions consider an aquitard with constant thickness and with radial-symmetrically increasing thickness outward from the center. The solution was used to predict the percolation rate from a peat layer around a bedrock outcrop in the James Bay Lowland near the De Beers Victor diamond mine. In this case the marine sediment layer limited the direct connection between the peat layer and the bedrock as an aquitard. Our zero order solution with constant marine sediment thickness showed the best fit to the steady state water level data of June 2012. It was found that the enhanced recharge around bioherms (i.e., at rates greater than the regional average of 0.7 mm/day) will only occur in marine sediments less than 4.3 m thick, for extreme depressurization of 30 m.     

A hybrid inversion method of damped least squares with simulated annealing used for Rayleigh wave dispersion curve inversion

Surface wave methods are becoming increasingly popular in many geotechnical applications and in earthquake seismology due to their noninvasive characteristics.Inverse surface wave dispersion curves are a crucial step in most surface wave methods.Many inversion methods have been applied to surface wave dispersion curve inversion,including linearized inversion and nonlinearized inversion methods.In this study,a hybrid inversion method of Damped Least Squares(DLS) with Very Fast Simulated Annealing(VFSA) is developed for multi-mode Rayleigh wave dispersion curve inversion.Both synthetic and in situ fi eld data were used to verify the validity of the proposed method.The results show that the proposed method is superior to the conventional VFSA method in aiming at global minimum,especially when parameter searching space is adjacent to real values of the parameters.The advantage of the new method is that it retains both the merits of VFSA for global search and DLS for local search.At high temperatures,the global search dominates the runs,while at a low temperatures,the local search dominates the runs.Thus,at low temperatures,the proposed method can almost directly approach the actual model.     

Assessment of groundwater recharge processes through karst vadose zone by cave percolation monitoring

Recharge processes of karst aquifers are difficult to assess given their strong heterogeneity and the poorly known effect of vadose zone on infiltration. However, recharge assessment is crucial for the evaluation of groundwater resources. Moreover, the vulnerability of karst aquifers depends on vadose zone behaviour because it is the place where most contamination takes place. In this work, an in situ experimental approach was performed to identify and quantify flow and storage processes occurring in karst vadose zone. Cave percolation monitoring and dye tracing were used to investigate unsaturated zone hydrological processes. Two flow components (diffuse and quick) were identified and, respectively, account for 66% and 34% of the recharge. Quickflow was found to be the result of bypass phenomenon in vadose zone related to water saturation. We identify the role of epikarst as a shunting area, most of the storage in the vadose zone occurring via the diffuse flow component in low permeability zones. Relationship between rainfall intensity and transit velocity was demonstrated, with 5 times higher velocities for the quick recharge mode than the diffuse mode. Modelling approach with KarstMod software allowed to simulate the hybrid recharge through vadose zone and shows promising chances to properly assess the recharge processes in karst aquifer based on simple physical models.     

A general solution for tide‐induced groundwater fluctuation in an estuarine‐coastal confined and unconfined aquifer system

This paper presents an analytical solution to tide‐induced head fluctuations in a two‐dimensional estuarine‐coastal aquifer system that consists of an unconfined aquifer and a heterogeneous confined aquifer extending under a tidal river with a semipermeable layer between them. This study considers the joint effects of tidal‐river leakage, inland leakage, dimensionless transmissivity between the tidal‐river and inland confined aquifer, and transmissivity anisotropic ratios. The analytical solution for this model is obtained via the separation of variables method. Three existing solutions related to head fluctuation in one‐ or two‐dimensional leaky confined aquifers are considered as special cases in the present solution. This study shows that there is a threshold of tidal‐river confined aquifer length. When the tidal‐river length is greater than the threshold length, the inland head fluctuations remain sensitive to the leakage effect but become insensitive to the tidal‐river width and dimensionless transmissivity. Considering leakage and transmissivity anisotropy, this study also demonstrates that at a location farther from the river–inland boundary, head fluctuations increase with increasing leakage and transmissivity anisotropy; the maximum head fluctuation occurs when leakage and transmissivity anisotropy are both at their maximum values. The combined action of the 3 effects of loading, tidal‐river aquifer leakage, and inland aquifer leakage differs significantly according to various aquifer parameters. The analytical solution in this paper can be applied to demonstrate the behaviours of the head fluctuations of an estuarine‐coastal aquifer system, and the head fluctuations can be clearly described when the tidal and hydrogeological parameters are derived from field measurement data or hypothetical cases.     

Numerical modelling of stream–aquifer interaction: Quantifying the impact of transient streambed permeability and aquifer heterogeneity

Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.     

Rapid estimation of soil heavy metal nickel content based on optimized screening of near-infrared spectral bands

In near-infrared spectroscopy,the traditional feature band extraction method has certain limitations.Therefore,a band extraction method named the three-step extraction method was proposed.This method combines characteristic absorption bands and correlation coefficients to select characteristic bands corresponding to various spectral forms and then uses stepwise regression to eliminate meaningless variables.Partial least squares regression(PLSR)and extreme learning machine(ELM)models were used to verify the effect of the band extraction method.Results show that the differential transformation of the spectrum can effectively improve the correlation between the spectrum and nickel(Ni)content.Most correlation coefficients were above 0.7 and approximately 20%higher than those of other transformation methods.The model effect established by the feature variable selection method based on comprehensive spectral transformation is only slightly affected by the spectral transformation form.Infive types of spectral transformation,the RPD values of the proposed method were all within the same level.The RPD values of the PLSR model were concentrated between 1.6 and 1.8,and those of the ELM model were between 2.5 and2.9,indicating that this method is beneficial for extracting more complete spectral features.The combination of the three-step extraction method and ELM algorithm can effectively retain important bands associated with the Ni content of the soil.The model based on the spectral band selected by the three-step extraction method has better prediction ability than the other models.The ELM model of the first-order differential transformation has the best prediction accuracy(RP^2=0.923,RPD=3.634).The research results provide some technical support for monitoring heavy metal content spectrum in local soils.     

A simple LMS‐based approach to the structural health monitoring benchmark problem

A structure's health or level of damage can be monitored by identifying changes in structural or modal parameters. However, the fundamental modal frequencies can sometimes be less sensitive to (localized) damage in large civil structures, although there are developing algorithms that seek to reduce this difficulty. This research directly identifies changes in structural stiffness due to modeling error or damage using a structural health monitoring method based on adaptive least mean square (LMS) filtering theory. The focus is on computational simplicity to enable real‐time implementation. Several adaptive LMS filtering based approaches are used to analyze the data from the IASC–ASCE Structural Health Monitoring Task Group Benchmark problem. Results are compared with those from the task group and other published results. The proposed methods are shown to be very effective, accurately identifying damage to within 1%, with convergence times of 0.4–13.0 s for the twelve different 4 and 12 degree of freedom benchmark problems. The resulting modal parameters match to within 1% those from the benchmark problem definition. Finally, the methods developed require 1.4–14.0 Mcycles of computation and therefore could easily be implemented in real time. Copyright © 2004 John Wiley & Sons, Ltd.     

建立基于土地利用的交通发生量模型的研究

针对出行生成预测中变量众多以及多重相关性问题,以常州调查数据为样本,分别采用逐步回归分析方法和偏最小二乘回归（PLSR）方法,建立基于土地利用形态的交通发生量模型。并对两者建模方法进行分析比较．结果表明PLSR能有效地解决变量多重相关性问题,模型意义清晰,效果显著。     

一种迭代算法在结构物理参数识别中的应用

本文在利用复合反演方法进行剪切型结构物理参数识别过程中,引入一种迭代算法来求解关于参数的病态方程组,探讨了几种比例阻尼形式的比例系数识别的处理方法。以10层剪切型框架结构为例,利用该迭代算法研究了无阻尼、三种比例阻尼形式以及噪信比为5％、10％的高斯白噪声等情形下的参数识别问题。结果表明,该迭代算法能很好地求解参数的病态方程组,且对四种阻尼形式和两种噪信比的噪声有很好的适应性。     

重力异常分离的相关法

本文研究重力异常分离的相关法,它由测点坐标的多项式方程构建区域异常,由观测异常减去区域异常计算剩余异常,根据最大化已知测点的目标密度分界面深度与剩余异常的互相关原则来分离重力异常.互相关的最大化目标可写为最小二乘方程,并利用线性或非线性优化算法迭代求解出构建区域异常的多项式系数,进而可计算出全区的区域异常和目标密度分界面的剩余异常,实现重力异常分离.理论重力异常试验结果表明,相关法直接有效地利用了已知测点的目标密度分界面深度信息,分离效果良好、反映实际地质情况.