Seismic response of arch dams considering infinite radiation damping and joint opening effects

Effects of two important factors on earthquake response of high arch dams are considered and combined into one program. These factors are: effects of radiation damping of the infinite canyon and local non-linearity of the contraction joint opening between the dam monoliths. For modeling of rock canyon, the discrete parameters are obtained based on a curve fitting, thus allowing the nonlinear dam system to be solved in the time domain. The earthquake uniform free-field input at the dam-canyon interface is used. An engineering example is given to demonstrate the significant effects of the radiation damping on the structure response.     

Steady flow to a partially penetrating blind‐wall well in a confined aquifer

Wells in aquifers of loose collapsible sediment are cased so that they have a blind wall and gain water only from the bottom. The hydraulic gradient established at the bottom of these wells during pumping brings the aquifer materials in a quicksand state, which may cause abrasion of pipes and pumps and even the destruction of well structure. To examine the quicksand occurrence, an analytical solution for the steady flow to a partially penetrating blind‐wall well in a confined aquifer is developed. The validity of the proposed solution is evaluated numerically. The sensitivity of maximum vertical gradient along the well bottom in response to aquifer and well parameters is examined. The solution is presented in the form of dimensionless‐type curves and equations that can be easily used to design the safe pumping rate and optimum well geometry to protect the well against sand production. The solution incorporates the anisotropy of aquifer materials and can also be used to determine the hydraulic conductivity of the aquifer. Copyright © 2012 John Wiley & Sons, Ltd.     

Vertical earthquake response of megawatt‐sized wind turbine with soil‐structure interaction effects

The dynamic response of a wind turbine on monopile is studied under horizontal and vertical earthquake excitations. The analyses are carried out using the finite element program SAP2000. The finite element model of the structure is verified against the results of shake table tests, and the earthquake response of the soil model is verified against analytical solutions of the steady‐state response of homogeneous strata. The focus of the analyses in this paper is the vertical earthquake response of wind turbines including the soil‐structure interaction effects. The analyses are carried out for both a non‐homogeneous stratum and a deep soil using the three‐step method. In addition, a procedure is implemented which allows one to perform coupled soil‐structure interaction analyses by properly tuning the damping in the tower structure. The analyses show amplification of the ground surface acceleration to the top of the tower by a factor of two. These accelerations are capable of causing damage in the turbine and the tower structure, or malfunctioning of the turbine after the earthquake; therefore, vertical earthquake excitation is considered a potential critical loading in design of wind turbines even in low‐to‐moderate seismic areas. Copyright © 2015 John Wiley & Sons, Ltd.     

Three‐dimensional models of reservoir sediment and effects on the seismic response of arch dams

The important effects of bottom sediments on the seismic response of arch dams are studied in this paper. To do so, a three‐dimensional boundary element model is used. It includes the water reservoir as a compressible fluid, the dam and unbounded foundation rock as viscoelastic solids, and the bottom sediment as a two‐phase poroelastic domain with dynamic behaviour described by Biot's equations. Dynamic interaction among all those regions, local topography and travelling wave effects are taken into account. The results obtained show the important influence of sediment compressibility and permeability on the seismic response. The former is associated with a general change of the system response whereas the permeability has a significant influence on damping at resonance peaks. The analysis is carried out in the frequency domain considering time harmonic excitation due to P and S plane waves. The time‐domain results obtained by using the Fourier transform for a given earthquake accelerogram are also shown. The possibility of using simplified models to represent the bottom sediment effects is discussed in the paper. Two alternative models for porous sediment are tested. Simplified models are shown to be able to reproduce the effects of porous sediments except for very high permeability values. Copyright © 2004 John Wiley & Sons, Ltd.     

Analytical and statistical approach for evaluating the effects of a river barrage on river–aquifer interactions

The construction of a river barrage can increase groundwater levels upstream of the barrage during the rainy season. Analytical and statistical approaches were applied to evaluate the relationship between groundwater and river water at the Changnyeong–Haman river barrage in Korea using time series data of water level and electrical conductivity from June 2011 to September 2014. An artificial neural network based time series model was designed to filter out the effect of rainfall from the groundwater level data in the study area. Aquifer diffusivity and river resistance were estimated from the analytical solution of a one‐dimensional unit step response function by using the filtered groundwater level data. River resistance increased in response to groundwater level fluctuations. Cross‐correlation analyses between the groundwater and the river water showed that the lag time increased during the observation period for both the water level and the electrical conductivity while the cross‐correlation function declined for the same period. The results indicated that a constant river stage maintained at the river barrage can weaken the hydrologic stress and reduce the exchange of material between the river and the adjacent aquifer because of the deposition of fine sediment on the river bottom and walls. Copyright © 2016 John Wiley & Sons, Ltd.     

锦屏一级水电站库区中小地震震源机制解与局部应力场特征

利用四川及邻区数字地震台网2013—2018年的波形资料, 基于CAP和HASH方法反演得到锦屏一级水电站(以下简称锦屏水库)库区及周边M_L≥2.5地震的震源机制解, 根据DRSSI方法计算得到库区局部应力场信息, 结合该区域的地质资料, 从震源机制参数和局部应力场变化等方面讨论了水库加卸载对该区域地震活动的影响.取得的主要认识如下: (1)库区地震活动类型以走滑型为主, 不同加卸载阶段震源机制类型存在一定差异, 库区蓄水前期震源机制类型呈现多样性, 蓄水后期与区域背景场趋于一致; (2)库区应力场最大主应力方向以NNW-SSE向为主, 分时段结果显示, 蓄水后期局部应力场最大主应力方向由NNW向逐渐发生偏转至NW向; (3)库区周边的地质构造和岩性条件有利于库水的渗透, 蓄水后重力载荷和流体渗透的共同作用导致岩体断层面库仑应力的变化, 从而呈现出库区蓄水前后中小地震活跃程度的显著差异特征.

     

Groundwater response to tidal fluctuations in a leaky confined coastal aquifer with a finite length

The study on the hydraulic properties of coastal aquifers has significant implications both in hydrological sciences and environmental engineering. Although many analytical solutions are available, most of them are based on the same basic assumption that assumes aquifers extend landward semi‐infinitely, which does not necessarily reflect the reality. In this study, the general solutions for a leaky confined coastal aquifer have been developed that consider both finitely landward constant‐head and no‐flow boundaries. The newly developed solutions were then used to examine theoretically the joint effects of leakage and aquifer length on hydraulic head fluctuations within the leaky confined aquifer, and the validity of using the simplified solution, which assumes the aquifer is semi‐infinite. The results illustrated that the use of the simplified solution may cause significant errors, depending on joint effects of leakage and aquifer length. A dimensionless characteristic parameter was then proposed as an index for judging the applicability of the simplified solution. In addition, practical application of the general solution for the constant‐head inland boundary was used to characterize the hydraulic properties of a leaky confined aquifer using the data collected from a field site at the Seine River estuary, France, and the versatility of the general solution was further justified.     

Discussion on ‘an investigation into the effects of damping and nonlinear geometry models in earthquake analysis’ by Andrew Hardyniec and Finley Charney

The discussed paper explores how various assumptions for damping and geometric nonlinearity affect the seismic collapse behavior of steel‐frame buildings. The recommended damping scheme is questioned, and an alternative is suggested. Additional explanation is sought to justify the differences observed between the P‐delta and co‐rotational models for geometric nonlinearity. Copyright © 2016 John Wiley & Sons, Ltd.     

On the interaction between a radiatively damped planetary wave and the zonally averaged circulation in the middle atmosphere

The interaction between a planetary wave damped by cooling to space and the zonally averaged circulation in the middle atmosphere is examined for a steady-state situation in middle latitudes. Quasi-geostrophic scaling of Type 2 is assumed (i.e. the space scales are planetary and the superrotation is small).A set of mean equations is derived for this scaling which is complementary to the set of perturbation equations previously studied. The mean equations show that a planetary wave induces a mean meridional circulation which is balanced by an eddy momentum forcing function and a mean diabatic heating which is balanced by an eddy heat flux forcing function. The vertical eddy fluxes enter the forcing at the same order as the horizontal eddy fluxes.An analytical wave solution is found for the case of an atmosphere in constant superrotation. The eddy fluxes and forcing functions are evaluated for this special case. It is found that they are very sensitive to the values of the radiative damping coefficient and the superrotation. Since the damping coefficient depends on the ozone concentration and the intensity of the solar ultraviolet flux, the results suggest that changes in these quantities can strongly modify the wave-mean flow interaction in the middle atmosphere. Possible implications for climate change are discussed.     

高烈度艰险山区跨断层隧道减震层减震技术研究

为进一步提高高烈度艰险山区跨断层隧道的减震性能,依托丽香铁路蒙古哨隧道工程,利用有限差分数值软件FLAC~(3D)对跨断层隧道施设减震层不同模式的减震效果进行研究。研究结果表明:减震层施设于初支与围岩间,其横向位移减小4.17%、竖向位移减小14.32%、最大主应力减小47.89%、最小主应力减小25.93%、最大剪切应力减小27.74%、最小安全系数提高34.62%～59.40%;减震层施设于初支与二衬间,其横向位移减小3.72%、竖向位移减小7.73%、最大主应力减小16.11%、最小主应力减小20.21%、最大剪切应力减小20.73%、最小安全系数提高0.04%～28.38%;减震层施设于初支与围岩间的减震效果优于施设于初支与二衬间。研究成果对于跨断层隧道的减震设计及减震技术有着重要的意义。     

Analytical solution of transient flow in a sloping soil layer with recharge

Abstract

An analytical solution of planar flow in a sloping soil layer described by the linearized extended Boussinesq equation is presented. The solution consists of the sum of steady-state and transient-series solutions, the latter in a separation-of-variables form, and can satisfy an arbitrary initial condition via collocation; this feature reduces the number of series terms, making the solution efficient. Key parameter is the dimensionless linearization depth η _o (R), R being the dimensionless recharge. The variable η _o (R), not the slope, characterizes the flow as kinematic or diffusive, and R ≈ 0.2 demarcates the two regimes. The transient series converges rapidly for large η _o (large R, near-diffusive flow) and slowly as η _o → 0 (kinematic flow). The quasi-steady (QS) state method of Verhoest & Troch is also analysed and it is shown that the QS depth profiles approximate the transient ones well, only if Δt exceeds a system-dependent transition time between flow states (possibly >>1 day). In an application example for a 30-day recharge series, the QS solution differs from the transient one by as much as 20% (RMSE = 15%), does not track recharge changes as well and fails to conserve mass.     

Analysis of solute transport in a divergent flow tracer test with scale‐dependent dispersion

It has been known for many years that dispersivities increase with solute displacement distance in a subsurface. The increase of dispersivities with solute travel distance results from significant variation in hydraulic properties of porous media and was identified in the literature as scale‐dependent dispersion. In this study, Laplace‐transformed analytical solutions to advection‐dispersion equations in cylindrical coordinates are derived for interpreting a divergent flow tracer test with a constant dispersivity and with a linear scale‐dependent dispersivity. Breakthrough curves obtained using the scale‐dependent dispersivity model are compared to breakthrough curves obtained from the constant dispersivity model to illustrate the salient features of scale‐dependent dispersion in a divergent flow tracer test. The analytical results reveal that the breakthrough curves at the specific location for the constant dispersivity model can produce the same shape as those from the scale‐dependent dispersivity model. This correspondence in curve shape between these two models occurs when the local dispersivity at an observation well in the scale‐dependent dispersivity model is 1·3 times greater than the constant dispersivity in the constant dispersivity model. To confirm this finding, a set of previously reported data is interpreted using both the scale‐dependent dispersivity model and the constant dispersivity model to distinguish the differences in scale dependence of estimated dispersivity from these two models. The analytical result reveals that previously reported dispersivity/distance ratios from the constant dispersivity model should be revised by multiplying these values by a factor of 1·3 for the scale‐dependent dispersion model if the dispersion process is more accurately characterized by scale‐dependent dispersion. Copyright © 2007 John Wiley & Sons, Ltd.     

An investigation into the effects of damping and nonlinear geometry models in earthquake engineering analysis

The growing emphasis of considering the behavior of structures at extreme performance states, such as collapse, has necessitated the characterization of the effects of varying attributes of the structural model. One source of variability that has not previously been considered is variability in the mathematical model. This study investigated the effects of changing the geometric nonlinearity approach and damping model on a four‐story buckling restrained braced frame, a four‐story steel moment resisting frame, and an eight‐story steel moment resisting frame. The variations in behavior are quantified using the maximum interstory drift ratio as the performance metric and qualified by comparing the relative displaced shapes and component response histories at the collapse performance state. Copyright © 2015 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.     

稀疏台网的中小地震震源机制解反演研究——以2017年1月4日西藏仲巴4.7级地震为例

本文应用ISOLA近震全波形方法,以2017年1月4日西藏仲巴4.7级地震为例,反演稀疏台网记录的中小地震震源机制解。该地震反演所得最佳双偶机制参数为：节面Ⅰ的走向109°/倾角85°/滑动角-177°,节面Ⅱ的走向19°/倾角88°/滑动角-4°,最佳矩心位置为30.590°N、83.784°E,最佳质心深度为6km,矩震级M_W4.6。震源机制反演结果表明此次地震是一次走滑型为主的事件,其与震源区域附近历史地震震源机制解具有相同性质。本文还应用CSPS初动扫描法,利用P波初动资料和近震波形联合约束反演此次地震的震源机制,并与ISOLA近震全波形反演结果进行比较,结果表明,联合少量台站的的三分量波形数据,能够定量地判断最佳震源机制解,降低了P波初动反演结果的非唯一性,同时也约束了由于少量台站参与全波形反演引起的解的不稳定性。本文研究为中小地震震源参数测定提供了一种简单有效的方法,具有较高的稳定性和可靠性。     

Spatio‐temporal changes in bog pool bottom topography – temperature effect and its influence on pool development: an example from a raised bog in Estonia

Increases in pool water and peat temperature in summer accelerate peat decomposition and production of biogenic gases, which can be trapped in peat pores and cause oscillation of peatland surfaces and the rise of peat from the bottom of bog pools. Associated changes in peat water conductivity, holding capacity and transpiration also affect bog hydrology. Our multi‐year study is the first to show in detail the extent and dynamics of changes in bog pool depth and bottom topography associated with changes in temperature, peat type and other factors. The true seasonal rise of peat from the pool bottom begins once the water temperature at the pool bottom exceeds 13–14 °C, although the speed and extent of the rise depends on peat properties, making the rise more erratic than its subsequent descent. The more rapid descent occurs after the first large drop in the temperature of the pool's surface water at the end of summer, resulting from the combination of reduced methane production and increased gas solubility with less influence by peat properties. Much higher dissolved organic carbon concentrations (216 ± 26 mg l^?1) in the pore water of peat risen from the bottom to the pool surface compared with that in the same type of peat at the pool bottom (62 ± 20 mg l^?1) indicate an acceleration of peat decomposition at the warmer pool surface. We show the extent and character of changes in pool depth and bottom topography and how annual differences relate to temperature. Only a few degrees' increase in pool water temperature could induce the pool bottom to rise faster and more extensively for a longer period and enhance decomposition in the peat at the pool surface. This should be evaluated in greater detail to assess the effects of temperature increase on the carbon budget and hydrology of peatlands. Copyright © 2012 John Wiley & Sons, Ltd.     

Tide‐induced groundwater head fluctuation in a coastal aquifer system with a submarine outcrop covered by a thin silt layer

We present an analytical solution of groundwater head response to tidal fluctuation in a coastal multilayered aquifer system consisting of an unconfined aquifer, a leaky confined aquifer and a semi‐permeable layer between them. The submarine outcrop of the confined aquifer is covered by a thin silt layer. A mathematical model and the analytical solution of this model are given. The silt layer reduces the amplitude of the hydraulic head fluctuation by a constant factor, and shifts the phase by a positive constant (time lag), both of which depend on the leakances of the silt layer and the semi‐permeable layer. The time lag is less than 1·5 h and 3·0 h for semi‐diurnal and diurnal sea tides respectively. When the leakance of the semi‐permeable layer or the silt layer assumes certain special values, the solution becomes the existing solutions derived by previous researchers. The amplitude of the hydraulic head fluctuation in the confined aquifer increases with the leakance of the silt layer and decreases with the leakance of the semi‐permeable layer, whereas the phase shift of the fluctuation decreases with both of them. A hypothetical example shows that neglecting the silt layer may result in significant parameter estimation discrepancy between the amplitude attenuation and the time‐lag fittings. Copyright © 2007 John Wiley & Sons, Ltd.     

Out-of-plane (SH) soil-structure interaction: a shear wall with rigid and flexible ring foundation

Soil-structure interaction (SSI) of a building and shear wall above a foundation in an elastic half-space has long been an important research subject for earthquake engineers and strong-motion seismologists. Numerous papers have been published since the early 1970s; however, very few of these papers have analytic closed-form solutions available. The soil-structure interaction problem is one of the most classic problems connecting the two disciplines of earthquake engineering and civil engineering. The interaction effect represents the mechanism of energy transfer and dissipation among the elements of the dynamic system, namely the soil subgrade, foundation, and superstructure. This interaction effect is important across many structure, foundation, and subgrade types but is most pronounced when a rigid superstructure is founded on a relatively soft lower foundation and subgrade. This effect may only be ignored when the subgrade is much harder than a flexible superstructure: for instance a flexible moment frame superstructure founded on a thin compacted soil layer on top of very stiff bedrock below. This paper will study the interaction effect of the subgrade and the superstructure. The analytical solution of the interaction of a shear wall, flexible-rigid foundation, and an elastic half-space is derived for incident SH waves with various angles of incidence. It found that the flexible ring (soft layer) cannot be used as an isolation mechanism to decouple a superstructure from its substructure resting on a shaking half-space.     

饱和土半空间中地下圆形衬砌洞室对平面SV波的散射

在Biot饱和多孔介质动力学理论的基础上,利用Fourier-Bessel级数展开法,通过对舍有衬砌洞室的局部场地进行波场分析,得到饱和土半空间中圆柱形衬砌洞室对平面SV波的散射问题的解析解。经验证,本文得到的解可以退化为半空间单相介质的情况。通过与已有的单相弹性介质半空间中圆柱形衬砌洞室对平面SV波散射问题的解析解的对比,验证了此解的正确性。在解析解的基础上,数值计算给出洞口动应力集中放大系数,分析了入射频率和孔洞埋深对柱面上的应力集中因子的影响。