A non-reflecting boundary condition for the finite element modeling of infinite reservoir with layered sediment

The design of seismic resistant concrete gravity dam necessitates accurate determination of hydrodynamic pressure developed in the adjacent reservoir. The hydrodynamic pressure developed on structure is dependent on the physical characteristics of the boundaries surrounding the reservoir including reservoir bottom. The sedimentary material in the reservoir bottom absorbs energy at the bottom, which will affect the hydrodynamic pressure at the upstream face of the dam. The fundamental parameter characterizing the effect of absorption of hydrodynamic pressure waves at the reservoir bottom due to sediment is the reflection coefficient. The wave reflection coefficient is determined from parameters based on sediment layer thickness, its material properties and excitation frequencies. An analytical or a closed-form solution cannot account for the arbitrary geometry of the dam or reservoir bed profile. This problem can be efficiently tackled with finite element technique. The need for an accurate truncation boundary is felt to reduce the computational domain of the unbounded reservoir system. An efficient truncation boundary condition (TBC) which accounts for the reservoir bottom effect is proposed for the finite element analysis of infinite reservoir. The results show the efficiency of the proposed truncation boundary condition.     

Seismic pressures on rigid cantilever walls retaining elastic continuously non-homogeneous soil: An exact solution

The dynamic response of an elastic continuously nonhomogeneous soil layer over bedrock retained by a pair of rigid cantilever walls to a horizontal seismic motion and the associated seismic pressure acting on these walls are determined analytically–numerically. The soil non-homogeneity is described by a shear modulus increasing nonlinearly with depth. The problem is solved in the frequency domain under conditions of plane strain and its exact solution is obtained analytically. This is accomplished with the aid of Fourier series along the horizontal direction and solution of the resulting system of two ordinary differential equations with variable coefficients by the method of Frobenius in power series. Due to the complexity of the various analytical expressions, the final results are determined numerically. These results include seismic pressures, resultant horizontal forces and bending moments acting on the walls. The solution of the problem involving a single retaining wall can be obtained as a special case by assuming the distance between the two walls to be very large. Results are presented in terms of numerical values and graphs using suitable dimensionless quantities. The effect of soil non-homogeneity on the system response is assessed through comparisons for typical sets of the parameters involved.     

Boundary element reservoir model for seismic analysis of gravity dams

The boundary element method has been successfully applied in the past to the analysis of hydrodynamic forces in two- and three-dimensional finite water reservoirs subjected to seismic ground motions. In extending the method to an infinite reservoir, the loss of energy due to pressure waves moving away towards infinity must be taken into account. In addition, for both finite and infinite reservoirs, energy is lost owing to partial absorption of the waves incident on a flexible bottom consisting of alluvial deposits. This paper presents the results of more recent research on the application of the boundary element method to the analysis of 2D reservoir vibration. Two different formulations are used: a constant boundary element formulation and a linear boundary element formulation. Special boundary conditions to treat infinite radiation and foundation damping have been incorporated in both formulations. Numerical results have been obtained for each of the two alternative formulations and compared against each other as well as with classical solutions and results obtained by other researchers.     

Fluid–foundation interaction in the seismic response of gravity dams

The seismic response of a dam is strongly influenced by its interaction with the water reservoir and the foundation. The hydrodynamic forces in the reservoir are in turn affected by radiation of waves towards infinity, wave absorption at the reservoir bottom, and cross-coupling between the foundation below the dam and the reservoir bottom. The fluid–foundation interaction effect, i.e. the wave absorption along the reservoir bottom, can be accounted for by using either an approximate one-dimensional (1D) wave propagation model or a rigorous analysis of interaction between the flexible soil along the base and the water. The rigorous approach requires enormous computational effort because of (a) cross-coupling between the foundation of the dam and the soil below the reservoir and (b) frequency dependence of the boundary condition along the fluid-foundation interface. The analysis can be simplified by ignoring the cross-coupling and by using the approximate 1D wave propagation model. The effects of each of these two simplifications on the accuracy and computational efficiency of the procedure used for the seismic response analysis of a dam are examined. Analytical results are presented for the complex frequency-response functions as well as the time histories of the response of Pine Flat dam to Taft and E1 Centro ground motions.     

含峭壁V形峡谷对地震SH波散射的解析解

地表地形常引起地震动的局部放大,这是由于地震波传播至局部地形时产生了散射现象.本文利用波函数展开方法和区域匹配技术,提出了含峭壁V形峡谷对平面SH波散射问题的解析解,并进行了退化验证.通过频域内的参数分析,揭示了峭壁深度、入射波频率和角度等因素对峡谷场地地面运动的影响规律,发现上部峭壁会增强峡谷对地震动的地形放大效应....     

基于等效单自由度体系的结构地震输入能量的研究

通过对多自由度体系结构在地震作用下的能量分析,考虑地震动因素和结构自身特性,提出一种框架结构地震总输入能量的简化求解方法.实例计算结果表明:基于等效单自由度体系来估计多自由度体系的地震总输入能,是一种计算简便且较为准确的方法.     

古田地震水口水电站重力坝强震观测记录的谱结构特征

通过对水口水电站重力坝强震反应台站在古田地震中获取的强震反应观测资料进行信噪比、反应谱和功率谱分析,得到如下结论:①大坝在0.7～15 Hz频率段的振动特性较为可信;②坝基和自由场输入地震动富于高频,峰值加速度反应谱存在较大差异;③坝基输入地震动存在差异性,建议今后此类大坝抗震设计时考虑多点地震动输入;④单个卓越频率携带的能量对反应谱影响不大,反应谱是和输入地震动总能量相关的;⑤坝体刚度较大,此次地震中还处于线弹性状态。初步了解了强震记录的地震动特性和大坝结构的抗震性能,对认识水库地震近场地震动特性和重力坝地震反应有一定的参考意义。     

圆弧形凹陷在SH波入射下瞬态反应的解析解

针对现有地震波散射问题的解析解均在频域内给出,其适用的频带较窄,不能获取到可靠的瞬态反应解析解的问题,本文利用宽频带频域解析解进一步获取了不同空间点地震反应的宽频带频响函数,并基于傅里叶变换得出SH波入射下圆弧形凹陷地形瞬态动力响应的解析解.以上述为基础,研究了输入脉冲与实际地震动时凹陷的瞬态反应的基本特征,并重点分析...     

Seismic isolation effect of lined circular tunnels with damping treatments

The Havriliak-Negami model for dynamic viscoelastic material behavior and Biot＇s theory of poroelasticity are employed to develop an exact solution for three-dimensional scattering effect of harmonic plane P-SV waves from a circular cavity lined with a multilayered fluid-filled shell of infinite length containing viscoelastic damping materials and embedded within a fluid-saturated permeable surrounding soil medium. The analytical results are illustrated with numerical examples where the effects of liner/coating structural arrangement, viscoelastic material properties, liner-soil interface bonding condition, seismic excitation frequency, and angle of incidence on the induced dynamic stress concentrations are evaluated and discussed to obtain representative values of the parameters that characterize the system. It is demonstrated that incorporating viscoelastic damping materials with a low shear modulus in the constrained layer configuration is an efficient means of enhancing the overall seismic isolation performance, especially for near-normally incident seismic shear waves where the amplitudes of induced dynamic stresses may be reduced by up to one-third of those without isolation in a relatively wide frequency range. Some additional cases are considered and good agreements with solutions available in the literature are obtained.     

Effects of reservoir bottom absorption on earthquake response of concrete gravity dams

A procedure is presented to analyse the response of concrete gravity dams due to horizontal and vertical earthquake ground motion components considering dam-water interaction and partial absorption of hydrodynamic pressure waves at the reservoir bottom into the foundation medium. The effects of reservoir bottom absorption on the hydrodynamic force on a rigid dam are examined first. The harmonic response of an idealized dam cross-section is presented for a wide range of parameters characterizing the properties of the dam, the impounded water and the foundation medium. Based on these frequency response functions the effects of dam-water interaction and of reservoir bottom absorption in the response of dams due to horizontal and vertical components of ground motion are investigated.     

Research on seismic fluid identification driven by rock physics

Seismic fluid identification works as an effective approach to characterize the fluid feature and distribution of the reservoir underground with seismic data. Rock physics which builds bridge between the elastic parameters and reservoir parameters sets the foundation of seismic fluid identification, which is also a hot topic on the study of quantitative characterization of oil/gas reservoirs. Study on seismic fluid identification driven by rock physics has proved to be rewarding in recognizing the fluid feature and distributed regularity of the oil/gas reservoirs. This paper summarizes the key scientific problems immersed in seismic fluid identification, and emphatically reviews the main progress of seismic fluid identification driven by rock physics domestic and overseas, as well as discusses the opportunities, challenges and future research direction related to seismic fluid identification. Theoretical study and practical application indicate that we should incorporate rock physics, numerical simulation, seismic data processing and seismic inversion together to enhance the precision of seismic fluid identification.     

SH波入射下圆弧状沉积盆地地震反应瞬态解析解计算

地震波散射问题的解析解是研究局部场地、地形、盆地等不规则地层结构对地震动参数放大效应影响的重要理论工具.现有解析解大部分在频域内给出,无法直接用于研究不规则地层结构对地震动峰值、反应谱等参数的放大效应.本文基于平面SH波入射下圆弧状沉积盆地动力响应宽频带稳态解析解,通过Fourier变换,获取瞬态响应解析解.基于此,研...     

基于混合波场地震动输入技术的近海场地地震反应分析方法

基于人工边界子结构模型,提出一种利用混合波场实现近海场地中地震P波和SV波垂直输入的方法.该方法中用于波动输入的混合波场由计算模型两侧截断边界的自由波场和底面边界的入射波场构成,避免了不规则近海场地的自由波场求解.同时采用基于声流体单元的流固耦合算法模拟场地-海水动力相互作用,利用流体介质人工边界和黏弹性人工边界单元模...     

Effect of reservoir bottom on earthquake response of concrete dams

The absorption of hydrodynamic pressure waves at the reservoir bottom has dominant effects on the structural response of the dam when subjected to ground motion. In the present study, a model is proposed for the absorption effects of the reservoir bottom in the earthquake analysis of dams. The model utilizes the wave reflection coefficient approach and is based on the solution of the wave equation in a sediment layer of viscoelastic material with a constant thickness overlying an elastic, semi-infinite foundation. Numerical studies were conducted to evaluate the effect of the sediment layer thickness and material properties as well as the effect of reflection of waves from the underlying rock. It is shown that the current approach of assuming the wave reflection coefficient at the reservoir bottom based on the characteristics of the sediment material and excluding the effect of the reflected waves from the underlying rock, may significantly underestimate the seismic response of the dam.     

Seismic time–frequency analysis as a robust method to estimate the fluid saturation: A case study of carbonates reservoir,Iran

The propagation of seismic waves through a saturated reservoir compresses the fluid in the pore spaces. During this transition, parts of seismic energy would be attenuated because of intrinsic absorption. Rock physics models make the bridge between the seismic properties and petrophysical reality in the earth. Attenuation is one of the significant seismic attributes used to describe the fluid behaviour in the reservoirs. We examined the core samples using ultrasonic experiments at the reservoir conditions. Given the rock properties of the carbonate reservoir and experiment results, the patchy saturation mechanism was solved for substituted fluid using the theory of modulus frequency. The extracted relationship between the seismic attenuation and water saturation was used in time–frequency analysis. We performed the peak frequency method to estimate the Q factor in the Gabor domain and determined the water saturation based on the computed rock physics model. The results showed how the probable fault in the reservoir has stopped the fluid movement in the reservoir and caused touching the water‐bearing zone through drilling.     

Two-dimensional dynamic analysis of concrete gravity and embankment dams including hydrodynamic effects

An analysis procedure in the frequency domain is developed for determining the earthquake response of two-dimensional concrete gravity and embankment dams including hydrodynamic effects; responses of the elastic dams and compressible water are assumed linear. The dam and fluid domain are treated as substructures and modelled with finite elements. The only geometric restriction is that an infinite fluid domain must maintain a constant depth beyond some point in the upstream direction. For such an infinite uniform region, a finite element discretization over the depth is combined with a continuum representation in the upstream direction. The fluid domain model approximately accounts for interaction between the fluid and underlying foundation medium through a damping boundary condition applied along the reservoir bottom, while the dam foundation is assumed rigid. Several examples are presented to demonstrate the accuracy of the fluid domain model and to illustrate dam responses obtained from the analysis procedure.     

A model study on the effects of input motion on the seismic behaviour of tunnels

Tunnel behaviour under earthquake loading is affected by many factors such as shape, depth and stiffness of the tunnel lining and the nature of the input motion. However, current knowledge on the effects of these parameters on the seismic behaviour of tunnels is limited to lack of experimental or field data. Existing analytical methods are based on assumptions, the validity of which needs to be established using carefully conducted experimental studies and numerical analyses. This paper focuses on the effects of input motion characteristics on seismic behaviour of circular and square tunnels. Dynamic centrifuge tests were carried out on model tunnels using input motions of different amplitude and frequency. Accelerations and earth pressures around the tunnels were measured. Complementary Finite Element analyses were conducted with different types of input motions. Results show that magnitude of the maximum input acceleration plays a crucial role on the maximum and residual lining forces, which the tunnel experiences.     

地震低频伴影的数值模拟与应用

Strong low-frequency energy beneath a hydrocarbon reservoir is called a seismic low-frequency shadow and can be used as a hydrocarbon indicator （Tarter et al., 1979） bu the physical mechanism of the observed low-frequency shadow is still unclear. To stud） the mechanism, we performed seismic numerical simulation of geological models with a hydrocarbon-bearing zone using the 2-D diffusive-viscous wave equation which car effectively model the characteristics of velocity dispersion and transform the seismic dat~ centered in a target layer slice within a time window to the time-frequency domain by usinl time-frequency signal analysis and sort the frequency gathers to common frequency cubes. Then, we observe the characteristics of the seismic low-frequency shadow in the common frequency cubes. The numerical simulations reveal that the main mechanism of seismic lowfrequency shadows is attributed to high attenuation of the medium to high seismic frequency components caused by absorption in the hydrocarbon-filled reservoir. Results from a practical example of seismic low-frequency shadows show that it is possible to identify the reservoir by the low-frequency shadow with high S/N seismic data.     

Numerical simulation of earthquake excited dam-reservoirs with irregular geometries using an immersed boundary method

In this work, a ghost-cell immersed boundary method is proposed for the hydrodynamic response of earthquake excited dam-reservoirs. The numerical method employs a second order accurate two-step projection algorithm including compressibility effects in pressure field due to earthquake. The effects of reservoir bottom absorption are treated by introducing damping terms into the momentum equations. Hydrodynamic response of earthquake excited dam with a sloping face is simulated to demonstrate the accuracy of the present numerical method. Numerical results compared with previous numerical and analytical solutions show that the present immersed boundary method can accurately compute the hydrodynamic forces on inclined and curved dam faces including the effects of water compressibility and reservoir bottom absorption for the possibility of resonance. The proposed numerical method was shown to have significant advantages in computational time and memory usage for the hydrodynamic simulation of large dam-reservoirs with arbitrary geometries. Hydrodynamic forces on a double curvature arch dam subjected to real earthquake induced ground motion are also simulated to demonstrate the capability of the method.     

Correlations between Energy and Displacement Demands for Performance-Based Seismic Engineering

The development of a scientific framework for performance-based seismic engineering requires, among other steps, the evaluation of ground motion intensity measures at a site and the characterization of their relationship with suitable engineering demand parameters (EDPs) which describe the performance of a structure. In order to be able to predict the damage resulting from earthquake ground motions in a structural system, it is first necessary to properly identify ground motion parameters that are well correlated with structural response and, in turn, with damage. Since structural damage during an earthquake ground motion may be due to excessive deformation or to cumulative cyclic damage, reliable methods for estimating displacement demands on structures are needed. Even though the seismic performance is directly related to the global and local deformations of the structure, energy-based methodologies appear more helpful in concept, as they permit a rational assessment of the energy absorption and dissipation mechanisms that can be effectively accomplished to balance the energy imparted to the structure. Moreover, energy-based parameters are directly related to cycles of response of the structure and, therefore, they can implicitly capture the effect of ground motion duration, which is ignored by conventional spectral parameters. Therefore, the identification of reliable relationships between energy and displacement demands represents a fundamental issue in both the development of more reliable seismic code provisions and the evaluation of seismic vulnerability aimed at the upgrading of existing hazardous facilities. As these two aspects could become consistently integrated within a performance-based seismic design methodology, understanding how input and dissipated energy are correlated with displacement demands emerges as a decisive prerequisite. The aim of the present study is the establishment of functional relationships between input and dissipated energy (that can be considered as parameters representative of the amplitude, frequency content and duration of earthquake ground motions) and displacement-based response measures that are well correlated to structural and non-structural damage. For the purpose of quantifying the EDPs to be related to the energy measures, for comprehensive range of ground motion and structural characteristics, both simplified and more accurate numerical models will be used in this study for the estimation of local and global displacement and energy demands. Parametric linear and nonlinear time-history analyses will be performed on elastic and inelastic SDOF and MDOF systems, in order to assume information on the seismic response of a wide range of current structures. Hysteretic models typical of frame force/displacement behavior will be assumed for the local inelastic cyclic response of the systems. A wide range of vibration periods will be taken into account so as to define displacement, interstory drift and energy spectra for MDOF systems. Various scalar measures related to the deformation demand will be used in this research. These include the spectral displacements, the peak roof drift ratio, and the peak interstory drift ratio. A total of about 900 recorded ground motions covering a broad variety of condition in terms of frequency content, duration and amplitude will be used as input in the dynamic analyses. The records are obtained from 40 earthquakes and grouped as a function of magnitude of the event, source-to-site condition and site soil condition. In addition, in the data-set of records a considerable number of near-fault signals is included, in recognition of the particular significance of pulse-like time histories in causing large seismic demands to the structures.