Effect of wellbore storage and finite thickness skin on flow to a partially penetrating well in a phreatic aquifer

An analytical model is presented for the analysis of constant flux tests conducted in a phreatic aquifer having a partially penetrating well with a finite thickness skin. The solution is derived in the Laplace transform domain for the drawdown in the pumping well, skin and formation regions. The time-domain solution in terms of the aquifer drawdown is then obtained from the numerical inversion of the Laplace transform and presented as dimensionless drawdown–time curves. The derived solution is used to investigate the effects of the hydraulic conductivity contrast between the skin and formation, in addition to wellbore storage, skin thickness, delayed yield, partial penetration and distance to the observation well. The results of the developed solution were compared with those from an existing solution for the case of an infinitesimally thin skin. The latter solution can never approximate that for the developed finite skin. Dimensionless drawdown–time curves were compared with the other published results for a confined aquifer. Positive skin effects are reflected in the early time and disappear in the intermediate and late time aquifer responses. But in the case of negative skin this is reversed and the negative skin also tends to disguise the wellbore storage effect. A thick negative skin lowers the overall drawdown in the aquifer and leads to more persistent delayed drainage. Partial penetration increases the drawdown in the case of a positive skin; however its effect is masked by the negative skin. The influence of a negative skin is pronounced over a broad range of radial distances. At distant observation points the influence of a positive skin is too small to be reflected in early and intermediate time pumping test data and consequently the type curve takes its asymptotic form.     

Transient analysis for fluid injection into a dome reservoir

A dome-shaped layer can be selected as a storage site for fluid injection. In this study, we develop a mathematical model for simulating transient head distribution in a heterogeneous and anisotropic dome-shaped layer due to a constant-head injection in a fully penetrating well. In the model, a form of step change is adopted to approximate the upper and lower boundaries of the dome and then the layer is split into two regions. The Laplace-domain solution of the model is developed using the Laplace transform and method of separation of variables. The transient injection rate at wellbore can then be obtained based on Darcy’s law and Bromwich integral method. The predicted head contours from the head solution show significant vertical flow components near the location of step change in the dome reservoir. The results of sensitivity analysis indicate that the hydraulic conductivity is the most sensitive parameter and the specific storage is the least sensitive one to the injection rate after a short period of injection time. In addition, the injection rate for a dome reservoir is also very sensitive to the change of the height for the reservoir near the injection well (first region) at a very early injection time. In contrast, the injection rate is more sensitive to the change of the height of the second region than that of the first region at late time. This analytical solution may be used as a primary tool to assess the capacity of fluid injection to various dome reservoirs.     

Analytical models of contaminant transport in coastal aquifers

The Henry formulation, which couples subsurface flow and salt transport via a variable-density flow formulation, can be used to evaluate the extent of sea water intrusion into coastal aquifers. The coupling gives rise to nontrivial flow patterns that are very different from those observed in inland aquifers. We investigate the influence of these flow patterns on the transport of conservative contaminants in a coastal aquifer. The flow is characterized by two dimensionless parameters: the Péclet number, which compares the relative effects of advective and dispersive transport mechanisms, and a coupling parameter, which describes the importance of the salt water boundary on the flow. We focus our attention on two regimes – low and intermediate Péclet number flows. Two transport scenarios are solved analytically by means of a perturbation analysis. The first, a natural attenuation scenario, describes the flushing of a contaminant from a coastal aquifer by clean fresh water, while the second, a contaminant spill scenario, considers an isolated point source.     

考虑水平力作用的弹性地基梁位移及内力解析解

考虑水平力作用下弹性地基梁受力特点,针对Winkler地基模型不能考虑梁底面间与地基切向相互作用的不足,在地基与梁体接触面上引入一系列各自独立的水平弹簧。假定弹簧反力与梁底和地基间的切向相对位移成正比,根据弹性力学基本原理,结合小变形条件下梁体几何方程和线弹性本构关系,建立力学分析模型,导出基本微分方程组并利用微分算子级数法进行求解,获得了考虑水平力作用时地基梁位移及内力解析解。最后通过实例分析,验证了该解析解的正确性和合理性,并对切向地基反力系数等影响因素进行了探讨。其结果可为工程设计提供借鉴参考。     

面波的质点运动轨迹在小爆破识别与判据研究中的应用

本项工作主要研究宽频带数字地震记录中小当量爆破(即矿山爆破)与小震级地震事件的识别问题。以前我们利用北京遥测地震台网的短周期模拟记录资料和DAPS(台网中心数字化)系统(郑秀芬等,2000;王凤霞,2000)产出的资料作过一些研究(郑秀芬等,2006),取得了一定的进展。在此基础上,利用宽频带数字地震记录资料,进行了小当量爆破与小震级地震事件在近场源短周期面波质点运动轨迹方面的识别研究,并给出了一定的判据。     

摩擦与粘性流体耗能器串联体系动力分析

采用粘性屈服模型来模拟摩擦耗能器的力－速度关系，建立了摩擦与粘性流体耗能器串联耗能体系的动力分析方法．为考虑框架杆件和支撑材料的几何非线性，采用增量型Ｒｏｓｅｎｂｒｏｃｋ二级三阶半隐式Ｒｕｎｇｅ－Ｋｕｔｔａ法求解动力方程．比较了仅有摩擦耗能器体系与串联组合耗能体系的减振效果，分析了粘性流体耗能器参数对组合耗能体系减振效果的影响。     

Experimental and analytical study on the performance of particle tuned mass dampers under seismic excitation

A particle tuned mass damper (PTMD), which is a creative integration of a traditional tuned mass damper and an efficient particle damper in the vibration control area, is proposed. This paper presents a comprehensive study that involves experimental, analytical, and computational approaches. The vibration control effects of a PTMD that is attached to a five‐story steel frame under seismic input are investigated by a series of shaking table tests. The influence of some parameters (auxiliary mass ratio, gap clearance, mass ratio of particles to the total auxiliary mass, frequency characteristics, and amplitude level of the input) is explored, and the performance of the PTMD with/without buffered material is compared. The experimental results show that the PTMD can achieve significant damping effects under seismic excitations, and the bandwidth of the suppression frequency is expanded, showing the device's robustness and control efficiency. In addition, an approximately analytical solution that is based on the concept of an equivalent single‐particle damper is presented, and the method to determine the corresponding system parameters is introduced. A comparative study between experimental and numerical results is conducted to verify the feasibility and accuracy of this analytical model. Copyright © 2016 John Wiley & Sons, Ltd.     

粘粒含量对南京粉细砂液化影响的试验研究

南京粉细砂是一种以片状颗粒成分为主的粉细砂，与通常的圆形颗粒石英砂有一定区别，片状颗粒成分使得南京粉细砂具有各向异性的性质。通过对不同粘粒含量的南京粉细砂进行液化试验，分析其试验结果发现：对于粘粒含量ρc＝5％、10％和15％三种情况，在均压固结情况下，当南京粉细砂达到初始液化时，采用变形标准和孔压标准是一致的；而偏压固结时，只能采用变形标准作为南京粉细砂达到初始液化的标准，此时振动孔压仅达到围压的50％-70％。粘粒含量对南京粉细砂的抗液化强度影响很大，随粘粒含量增加，在不同固结比时，其抗液化强度并不是单调变化，而是在粘粒含量ρc＝10％左右处于一个低谷，其抗液化强度最低。     

Determination of the suspension competence of debris flows based on particle size analysis

The determination of the critical particle size between solid and fluid phases, i.e., the suspension competence, is fundamental for debris flow. A method for determining suspension competence based on particle size analysis is presented in this paper. Suspension competence of static experimental water-debris mixtures prepared with the sediment of Jiangjia Gully is -0.025 mm if the bulk density is less than 1,800 kg m-3 and it increases with bulk density of more concentrated mixtures. Suspension competence of natural debris flows in Jiangjia Gully increases exponentially with the bulk density. These two data sets are compared in order to understand the suspension mechanism. It is concluded that turbulence may play a leading role in particle suspension in non-viscous and sub-viscous debris flows, while in viscous debris flows both matrix strength and excess pore water pressure play important roles.     

Application of the qda–md method of environmental discrimination to particle size analyses of fine sediments by pipette and sedigraph methods: A comparative study

The results of particle size analyses by the low-technology pipette method and by the modern SediGraph, of a suite of fine lacustrine bottom sediment samples, are compared. Good correlations between the two data sets are observed, although the SediGraph consistently indicates finer grain sizes and better sorting of distributions. The QDa-Md method of environmental discrimination, originally established on the basis of traditional methods of size analysis, is found to be sufficiently versatile to accommodate the SediGraph data into the appropriate trend envelopes, despite the considerable differences in derived grain size parameters.     

基于粒子群优化的多目标体中梯电阻率异常反演

在采用中梯装置的电阻率剖面法应用中,特别是在环境、水文和工程等领域,经常遇到需要对多个异常目标体进行快速定位以便及时进行相应处理的情况.利用倾斜椭球体来近似模拟这些电阻率异常目标体,使得正演计算可以采用解析表达式来实现,提高正演计算的时效性.相应的地球物理模型也可简化为由椭球体个数、中心点位置、倾角、轴径以及电阻率等参数构成的粒子,多个粒子组成的粒子群在粒子群优化算法的控制下在给定的搜索空间中并行地搜索最优模型.通过对粒子群优化算法参数的合理设计,利用其良好的全局与局部均衡的搜索能力实现对多个异常目标体的同时反演.数值实验结果表明该反演方法能有效实现对多个目标体的同时反演,计算速度快、反演拟合精度较高、同时具有一定的抗噪音能力.快速的多目标体反演,可以实时准确的定量解译中梯剖面法圈定的异常目标体,较好地满足工程等领域的高时效性要求.     

Application of multivariate analysis to suspended matter particle size distribution in a karst aquifer

Numerous quantitative and qualitative variables control suspended sediment dynamics in karst systems. The objective of this study was to identify the hydrodynamic variables controlling the transport properties of particles in a karst aquifer (western Paris basin). The particle size distribution of suspended sediment infiltrating via a swallow hole was compared to that in discharge from a spring, allowing identification of the particle transport properties of the karst system. Hill and Smith analysis, a type of multivariate analysis that allows joint examination of quantitative and qualitative variables, was used to identify the hydrodynamic parameters controlling the transport properties of the suspended matter. The results demonstrate that the particle size distribution discharging at the karst spring is controlled by spring discharge and the hydraulic gradient of the system. The hydraulic gradient is defined by the piezometric level and the stage of the Seine River, which is in turn controlled by the tide. This study illustrates the use of Hill and Smith analysis to identify those variables which control suspended sediment transport. It also illustrates the application of this analysis to identify boundary conditions and evaluate variables which control the behaviour of the hydrologic system. Copyright © 2007 John Wiley & Sons, Ltd.     

液体黏滞阻尼器在大跨度桥梁中对伸缩缝保护作用的探讨

为了限制动荷载作用下梁端的位移,经常要在梁端安装阻尼器,与伸缩缝在结构上形成并联装置。当桥面突然发生大规模的刹车,或在突加的地震、风荷载等荷载的冲击作用下,阻尼器是否会对伸缩缝起到保护作用就成为一个令人关心的问题。由一般多自由度体系动力方程出发,通过矩阵变换得到结构局部关于阻尼器单自由度方程。基于此单自由度方程,对阻尼器在结构中具有局部减震效果提出了理论上的论证,在此基础上就阻尼器在各种冲击荷载作用下对伸缩缝的保护作用进行了探讨。     

Analytical and numerical analysis of pressure drawdown in a poroelastic reservoir with complete overburden effect considered

Starting from an analytical reservoir model that incorporates full interaction with an elastic overburden, a new hybrid mathematical approach is developed by combining two numerical discretization methods. A tabular reservoir (petroleum reservoir or an aquifer) in an infinite or semi-infinite domain is viewed as a macroscopic displacement discontinuity, allowing use of the efficient displacement discontinuity mathematical method to calculate stresses and displacements that arise because of pressure changes. A 3-D finite element method using a poroelastic formulation is used to discretize the reservoir itself. By coupling the displacement discontinuity and finite element methods, a 3-D large-scale poroelastic reservoir can be simulated within an infinite or semi-infinite domain. The numerical model has been verified through comparison to known solutions, and some time-dependent pressure drawdown problems are analyzed. Results indicate that including the complete overburden (reservoir surroundings) response has a significant effect on pressure drawdown in a poroelastic reservoir during pumping, and should be incorporated in appropriate applications such as well test equations and subsidence analyses.     

Effects of spatial wind inhomogeneity and turbulence anisotropy on circulation in an elongated basin: A simplified analytical solution

The response of a closed basin to wind forcing has been studied extensively, but the role of turbulence anisotropy coupled with spatial wind inhomogeneity has never been considered explicitly. An analytical solution is presented for steady state hydrodynamics, considering the central part of a constant-depth, elongated basin with a laterally varying wind, and a homogeneous, yet anisotropic, eddy viscosity tensor. The solution is derived both for a single-layer and for a two-layer basin, which is representative of a stratified lake with a well-developed thermocline. Since the focus is on the short-term barotropic reaction to wind forcing, which determines the type of lake circulation, baroclinic effects are neglected as a first approximation. In this case, the development of planimetric (depth-averaged) circulation superimposed on circulation in the vertical plane can be determined as a function of wind lateral variation and a few dimensionless parameters. The relevance of such an analytical solution is twofold. Firstly, the knowledge of the prevailing circulation can help in the choice of the best type of numerical model (three-dimensional vs. two-dimensional, depth- or lateral-averaged). Secondly, it shows the importance of correct estimates of both vertical and horizontal eddy viscosity, whereas the latter is not usually considered as an important parameter in lake hydrodynamics modeling.     

A study on characteristics of tectonic block motion and tectonic setting of strong earthquakes in northern part of the Shanxi fault depression zone

Introduction Shanxi fault depression zone (SFDZ) is one of important Cenozoic fault basin zones and strong earthquake belts in Chinese mainland. Its northern part has aroused wide research interests due to the complicated tectonics and high activity of strong earthquakes there. Early researches on this depression zone were carried out since 60s of last century (DENG, et al, 1973; DENG, YOU, 1985; LU, DING, 1985; XU, 1990; XU, et al, 1996, 2002). In 90s of last century, the geologica…     

保定市设计地震动参数的研究

本文根据保定市及周围地区的地震地质环境，在地震危险性分析的基础上，采用等效线性一维波动方程进行土层的地震反应分析。给出50年超越概率63％、10％、2％基岩和地面的水平向峰值加速度、反应谱（场址基本烈度Ⅶ度）和地震影响系数最大值。该结果为抗震设计提供了可靠依据，具有应用价值。     

Effect of particle density and inflow concentration of suspended sediment on bedload transport in rill flow

Laboratory flume experiments were carried out to evaluate the effect of particle density on bedload transport of sand‐sized particles and the effect of a suspended load of clay particles (kaolinite) on bedload transport of sand‐sized particles in rill flow conditions. Three materials in the range 400–600 µm were selected to simulate bedload transport of primary particles and aggregates: sand (2650 kg/m³), crushed brick (2450 kg/m³) and anthracite (1300–1700 kg/m³). In the two first experiments, two different methods were applied to determine bedload transport capacity of coarse particles for various conditions of flow discharge (from 2 to 15 L/min) and slope (2.2, 3 and 4%). In the third experiment, clear water was replaced with kaolinite–water mixture and bedload transport capacity of crushed brick particles was determined for a 4% slope and different concentrations of kaolinite (0, 7, 41 and 84 g/L). The results showed that bedload transport increased significantly with the decrease in particle density but the effect of particle density on transport rates was much less important than flow discharge. Velocity measurements of clear flow, flow mixed with coarse particles and coarse particles confirmed the existence of a differentiation between suspended load and bedload. In these experimental conditions, suspended load of kaolinite did not affect bedload rates of crushed brick particles. Three transport capacity formulae were tested against observed bedload rates. A calibration of the Foster formula revealed that the shear stress exponent should be greater than 1.5. The Low and the Govers unit stream power (USP) equations were then evaluated. The Low equation was preferred for the prediction of bedload rates of primary particles but it was not recommended in the case of aggregates of low density because of the limited experimental conditions applied to derive this equation. Copyright © 2008 John Wiley & Sons, Ltd.