山西漫水井水位微动态特征及干扰因素分析

选取2008-2018年漫水井水位数据进行统计分析发现,漫水井水位受多种干扰因素影响,包括降雨、雷电、抽水干扰,传感器故障,路由器和数采死机。对其微动态特征研究发现漫水井水位长期趋势特征受抽水干扰及降雨共同影响,月动态、日动态特征可以看到明显的潮汐效应,同震响应能力较好。了解其微动态特征并准确识别相关干扰因素有助于识别地震前兆信息。     

天津井水位同震效应分析

通过整理天津张道口、静海、宝坻等多口观测井的数字化水位的观测资料,并针对全球近年来发生的8级以上地震分析了井水位同震变化的特征。结果表明,天津地区这几口观测井对远、大地震水位同震效应显著,地震中水震波变化为振荡型与震荡阶变型。     

天津井水位同震效应分析

通过整理天津张道口、静海、宝坻等多口观测井的数字化水位的观测资料,并针对全球近年来发生的8级以上地震分析了井水位同震变化的特征。结果表明,天津地区这几口观测井对远、大地震水位同震效应显著,地震中水震波变化为振荡型与震荡阶变型。     

天津两井水位同震效应分析

从不同的方位对天津地区两口地震观测专用井记录到的国外7级、国内6级以上的地震进行分析。两口井所处的地理位置、地质结构、井深及其含水层岩性不同,对同一地震的表现形式也不尽相同。通过对比分析,得到井水位对地震的响应特点及几点认识：观测井孔一般对远场大震记录比较灵敏,而对近震及地方震反映不灵敏;井孔水位仪记录的水位波幅与震级、震中距有关;井孔同震效应的能力和井孔对固体潮的反映能力不一致。     

Geoelectrical exploration for groundwater in Al Maha Forest,Ain Jouhra,Morocco

An area of about 30 km² located in Ain Jouhra, south of Rabat, Morocco, was the subject of a geoelectric resistivity investigation. The main goal of the investigation was the assessment of the groundwater potential of the uppermost aquifer. The aquifer conditions such as depth, thickness and boundaries were also investigated. The obtained apparent resistivity curves were first analysed qualitatively and classified using simple curve shapes. Thereafter, the data were converted to resistivity and thickness pairs semi‐quantitatively by means of master curves and then quantitatively by computer modelling using ATO and Winsev software (Zohdy, 1989; Zohdy and Bisdrof, 1989). Lithological control from the available single well with a stratigraphic log aided in the correlation of the resistivity values to different rock units. Three different AB‐spacing iso‐resistivity maps, an isopach map of the main groundwater‐bearing horizon, the depth to the aquifer substratum map and five geoelectric cross‐sections were constructed. The interpretation of these soundings indicates the presence of an unconfined to semi‐confined sandy aquifer with relatively important extent and varying thickness. The maximal thickness of the aquifer is recorded in the central part of the investigated area and is thinning southwards to pinch out farther to the south. Geophysical as well as field data indicate a hydraulic connection between the upper and deeper aquifers. Indeed, the two aquifers are separated from each other by a marly substratum that is indicated throughout the area by the lowest values of the interpreted true resistivity. The value of this resistivity varies laterally, most likely due to the lateral variation in the shale‐to‐sand ratio. The altitude of the substratum decreases towards the north, and increases southwards. Regarding the availability of the groundwater in the study area, zones with high potential are theoretically expected to occur in the central part where the transversal resistance is greatest. However, sufficient water supply and high flow rates from wells intended to produce restrictively from the most upper aquifer are not likely to exist. This conclusion, which seems to be very pessimistic, is evidenced from two real field and experimental observations. The first is the rapid fall of the level of Gharnoug lake, despite the ongoing feeding by three wells. Hence, the amount of water level drop cannot be accounted for by the evaporation alone. That means that the deeper aquifer is continuously draining the upper aquifer at a high flow rate. Very low rates are recorded in all the wells that penetrated only the upper aquifer, the exception being the well that reached deeper into the lower aquifer. The flow rate in this lower aquifer measured 18 litre s^?1. Copyright © 2007 John Wiley & Sons, Ltd.     

Obtaining the steady‐state drawdown solutions of constant‐head and constant‐flux tests

The solutions of constant‐head and constant‐flux tests are commonly used to predict the temporal or spatial drawdown distribution or to determine aquifer parameters. Theis and Thiem equations, for instance, are well‐known transient and steady‐state drawdown solutions, respectively, of the constant‐flux test. It is known that the Theis equation is not applicable to the case where the aquifer has a finite boundary or the pumping time tends to infinity. On the other hand, the Thiem equation does not apply to the case where the aquifer boundary is infinite. However, the issue of obtaining the Thiem equation from the transient drawdown solution has not previously been addressed. In this paper, the drawdown solutions for constant‐head and constant‐flux tests conducted in finite or infinite confined aquifers with or without consideration of the effect of the well radius are examined comprehensively. Mathematical verification and physical interpretation of the solutions to these two tests converging or not converging to the Thiem equation are presented. The result shows that there are some finite‐domain solutions for these two tests that can converge to the Thiem equation when the time becomes infinitely large. In addition, the time criteria to give a good approximation to the finite‐domain solution by the infinite‐domain solution and the Thiem equation are investigated and presented. Copyright © 2008 John Wiley & Sons, Ltd.     

Modelling groundwater levels in an urban coastal aquifer using artificial neural networks

The prediction of groundwater levels in a basin is of immense importance for the management of groundwater resources, especially in coastal regions where the water table fluctuations are to be limited to avoid seawater intrusion. In this paper, an Artificial Neural Network (ANN) methodology is presented to predict groundwater levels in individual wells with one month lead. Groundwater levels were also predicted in neighboring wells using model parameters from the best network of a well. This methodology is applied to an urban coastal aquifer in Andhra Pradesh state, India. The results suggest that the feed forward neural network with Levenberg Marquardt (LM) algorithm is a good choice for predicting groundwater levels in individual wells. Bayesian Regularization (BR) model parameters of Balaji Nagar well are also used successfully to predict groundwater levels in the study area. It was observed that the ANN‐based algorithms were a better choice for the prediction of groundwater levels with limited hydrological parameters. Copyright © 2007 John Wiley & Sons, Ltd.     

Optimal pumping locations of skimming wells

Abstract

A real-life problem involving pumping of groundwater from a series of existing wells along a river flood plain underlain with geologically saline water is examined within a conceptual framework. Unplanned pumping results in upconing of saline water. Therefore, it is necessary to determine optimal locations of fixed capacity pumping wells in space and time from a set of pre-selected candidate wells that minimize total salinity concentration in space and time. The nonlinear, non-convex, combinatorial problem involving zero—one decision variables is solved in a simulation—optimization (S/O) framework. Optimization is accomplished by using simulated annealing (SA)—a search algorithm. The computational burden is primarily managed by replacing the numerical model with a surrogate simulator—artificial neural network (ANN). The computational burden is further reduced through intuitive algorithmic guidance. The model results suggest that the skimming wells must be operated from optimal locations such that they are staggered in space and time to obtain least saline water.     

Using the extended finite element method for simulation of transient well leakage in multilayer aquifers

The extended finite element (XFEM) is applied to the problem of transient leakage from abandoned or free-flowing artesian wells in perforated aquifer-aquitard systems. To more accurately capture the singularities in potentiometric head at the wells, the standard linear finite element basis is locally augmented with asymptotic analytical solutions which enable more accurate calculations of leakage rates between aquifers. Highly accurate flux estimates are obtained without the need for higher mesh resolution near wells. Simulations are carried out to test both the accuracy and convergence properties of the XFEM implementation, and the XFEM results are compared to those of a high-resolution standard finite element model. It is seen that for the type of singularity-driven problem posed here, the standard FEM is unable to resolve leakage rates without very fine discretization, but that the XFEM performs robustly with fewer degrees of freedom. The impact of aquifer geometric heterogeneity on leakage rates is assessed and seen to be an important factor in determining total leakage. It is demonstrated that the XFEM may be a valuable tool in many water resources applications where small-scale effects can impact global system behavior.     

开发地震构造建模技术的研究方法

某油田已注水开发40年,高井网密度条件下仍对某些微构造、断层认识不足。运用三维地震解释资料与开发井资料解释两项技术结合,优势互补,相互验证,形成了一套开发地震可视化构造建模技术,应用该项技术显著提高了断层的首尾部延伸、复杂断裂带断层的组合关系、无井断点控制的小断层识别精度,对油田特高含水期提高钻井成功率、油田精细调整挖潜具有重要意义。     

基于黏弹滑移界面理论揭示套管中模式波的传播特征

目前固井质量评价主要是利用套管中模式波的幅度或衰减变化判断套管与水泥间的胶结状况.在超声测量频段（70~500 kHz）,现有声波测井仪器工作时在套管中激发的模式波类似薄板中的对称兰姆波（拉伸波）或反对称兰姆波（弯曲型兰姆波）.本文将黏弹滑移界面理论应用到套管和水泥的耦合界面,探讨了套管中的准兰姆波（包括拉伸波和弯曲型兰姆波）与套管后物质的耦合方式,拉伸波与套后物质主要以剪切刚度耦合,即便是0.01 mm的微环也会使得拉伸波的衰减明显降低,接近套后是流体时的响应特征;而弯曲型兰姆波与套管后物质主要以法向刚度耦合,这使得在套后胶结轻质水泥时弯曲型兰姆波对微环不敏感,在套后耦合常规水泥时微环的存在又会使得弯曲型兰姆波的衰减明显增强,衰减的增加说明其与套后物质的声耦合更紧密,也表明微环与常规水泥的等效声阻抗适当降低,更有利于弯曲型兰姆波传播时与套后介质实现法向刚度的匹配.拉伸波和弯曲型兰姆波对同一胶结状况的不同响应特征也为固井质量综合评价时区分不同的水泥胶结状况提供了可能.     

1994年2月16日共和5.8级地震短临预报实况

本文真实地再了１９９４年２月１６日青海共和５．８级地震前前兆异常的发展过程及对该地震的短临预报过程，同时系统地分析了前兆异常在时间，空间及种数上的分布特征。指出，在有一定前兆监测能力的地区，对一些较大地震作出一定的或成功的短临预报是有希望的。     

潜水泵变频稳流抽水观测技术系统

潜水泵变频稳流抽水观测技术系统适用于不同观测条件,可以保证原有水文地球化学各测项的数字化和模拟观测,还可实施数字化水位与水温的同井观测,可以记录清晰的潮汐效应和同震效应,不低干天然自流状态下观测资料质量水平,是值得推广的新技术.     

井孔水位记震能力变化与近大地震关系研究

通过对云南省及邻近地区7级大震,震前群井水位记震能力变化的研究及异常、预报指标检验,表明震前记震能力减弱或消失与近大地震有关。其异常特征表现为形态相似、时间同步、幅度明显等特点。井孔水位记震能力变化可能是地下水中期预报地震的一种新方法。     

关于地震地下水观测井建设的若干问题

根据有关的国家与地震行业技术标准,结合建井的实践经验,本文论述了地震地下水观测井的定位、观测井结构设计与钻井施工的技术要求。     

UNSTEADY MOVEMENT OF FRESH WATER IN THICK UNCONFINED SALINE AQUIFERS

Abstract

Fresh-water lenses are formed in unconfined saline aquifers in response to deep percolation from rainfall, artificial recharge, and seepage from irrigation waters and/or in response to injecting fresh water through vertical or horizontal wells. An approximate differential equation is derived in terms of the depth of the fresh-salt water interface below the initial position of the saline-water table. This equation is analogous to that of the ground-water motion in two dimensions. The wealth of knowledge available from solving the latter equation is used to obtain approximate expressions for the movement of the fresh-salt water interface in several flow systems wherein this interface does not reach the bottom of the aquifer. These approximate solutions as well as others for related quantities of interest may afford useful tools for rationally planning the extraction of usable waters from such flow systems.     

唐山地震前不同区域油水井异常变化及其机理分析

对唐山地震前胜利、大港、辽河3个油田的油井和渤中2井的资料以及北京附近部分水井的资料进行整理和分析,描述了油水井的短期和临震异常特征;分析了唐山地震前同一地区和不同地区间的油水井前兆异常之间的关系,对油水井异常变化的机理进行了研究和分析,认为不同的油水井前兆异常是由于唐山地震孕震过程中造成的不同地区的应力状态不同所致,油水井前兆异常可以用来监测和预报地震。     

Improvement on the estimation of constant-rate drawdown in large-diameter wells

Abstract

Papadopulos and Cooper's (PC) solution can be used to describe the drawdown resulting from pumping with a constant rate at a large-diameter well. However, this solution is too complicated to be accurately evaluated due to the oscillatory nature of the Bessel functions. The PC approach resulted in tabulated values of dimensionless drawdown at the well with an accuracy of four or fewer digits for selected values of dimensionless storage coefficient versus dimensionless time. Some researchers have fitted the tabulated values with interpolation formulas that are easy to use in engineering applications. Those formulas may be more accurate if the tabulated values are computed with greater accuracy. In this study, we propose an efficient numerical procedure, including a root search scheme, to find the roots of the integrand, Gaussian quadrature for numerical integration, and Shanks transform to accelerate convergence of infinite series. The proposed procedure can evaluate the dimensionless drawdown with greater accuracy and is useful in practice if there is a need for high accuracy for the observation either at the well or in the aquifer at some distance from the pumping well

Editor D. Koutsoyiannis

Citation Chang, Y.C., Yeh, H.D., and Wang, C.T., 2013. Improvement on the estimation of constant-rate drawdown in large-diameter wells. Hydrological Sciences Journal, 58 (3), 716–727.     

Horizontal rough fracture type curves for groundwater movement

Abstract

Type curves are derived analytically for radial flow in rough horizontal fractures toward a well. The basic assumptions are that there is no turbulent flow near the borehole and the well storage is ignored. The basis of the methodology is to write explicit expressions for the continuity and cubic law flow equations, which are combined using a Boltzmann transformation leading to a simple ordinary differential equation for groundwater movement. Solutions are presented as a set of type curves for different fracture apertures. It is observed that the solutions provide a method of uniquely identifying fracture hydraulic parameters when the fracture is smooth, but pose ambiguity for rough fracture parameter estimations. However, large time portions of these type curves appear as straight lines on semi-logarithmic paper, which provides a unique way for rough fracture parameter determination. Identification of the fracture parameters, namely, the aperture and relative roughness, is possible in a unique manner with the use of these lines and the dimensionless time drawdown concept. The cubic law is the asymptotic behaviour, either for large times or large fracture apertures. Prior to this asymptotic part, there is a non-cubic portion which gives rise to systematic deviations from the cubic law. The technique presented is useful, especially for evaluating pumping tests from a single major fracture isolated by packers.     

利用深井地下水位动态研究大华北地区现今构造应力场状态

在水平层状含水层(一维)模式下, 根据含水层应力与井水位变化之间的定量关系表达式, 利用华北地区40多口深井水位动态变化资料, 探索研究华北地区现今构造应力场状态。 即水位上升代表构造应力场压性增强张性减弱; 水位下降代表构造应力场压性减弱张性增强。 初步研究结果表明, 1995~1998年华北北部地区张家口-渤海构造带和山西构造带北部水位普遍上升, M_S5.0以上地震活动频繁; 1999年以后华北地区水位普遍处于下降状态, 反映了华北地区1999年以来构造应力场状态以张性为主, 整个华北地区M_S5.0以上地震活动处于相对平静状态。 该结果可为华北地区地震活动性预测研究提供应力场背景信息。