大型三维矩形、U形渡槽地震响应对比分析

本文应用任意拉格朗日-欧拉(ALE)有限元方法对比分析大型三维排架支撑矩形、U形渡槽结构地震响应,以获得渡槽中动水压力资料。研究表明,渡槽中水体的晃动幅度十分显著,水体的大幅晃动显著地改变两种截面形式渡槽中水压力、槽体应力值及U形渡槽动水压力和应力分布,对矩形渡槽动水压力和应力分布规律影响不大。动应力分布的改变直接影响到大型预应力U形截面渡槽结构中预应力钢筋的有效性。在渡槽的抗震设计中,应充分考虑U形截面渡槽动水压力和动应力分布变化问题,动水压力应该作为一项主要的外荷载予以考虑。     

Quantitative detection of fluid distribution using time-lapse seismic

Although previous seismic monitoring studies have revealed several relationships between seismic responses and changes in reservoir rock properties, the quantitative evaluation of time‐lapse seismic data remains a challenge. In most cases of time‐lapse seismic analysis, fluid and/or pressure changes are detected qualitatively by changes in amplitude strength, traveltime and/or Poisson's ratio. We present the steps for time‐lapse seismic analysis, considering the pressure effect and the saturation scale of fluids. We then demonstrate a deterministic workflow for computing the fluid saturation in a reservoir in order to evaluate time‐lapse seismic data. In this approach, we derive the physical properties of the water‐saturated sandstone reservoir, based on the following inputs: V_P, V_S, ρ and the shale volume from seismic analysis, the average properties of sand grains, and formation‐water properties. Next, by comparing the in‐situ fluid‐saturated properties with the 100% formation‐water‐saturated reservoir properties, we determine the bulk modulus and density of the in‐situ fluid. Solving three simultaneous equations (relating the saturations of water, oil and gas in terms of the bulk modulus, density and the total saturation), we compute the saturation of each fluid. We use a real time‐lapse seismic data set from an oilfield in the North Sea for a case study.     

Evaluation of vertical seismic response for a large‐scale beam‐supported aqueduct

By the theories of potential flow and structural vibration, the formulae for evaluating the ‘wet’ (with water) frequencies and mode shapes of the beam‐supported aqueduct are derived through a simplified fluid‐structure interaction analysis. The time‐history formulae of structural responses to the vertical seismic excitation are obtained. Applying the response‐spectrum principle, the equivalent vertical earthquake load exerted on the beam and the corresponding effects are also derived. Several illustrative examples are conducted. The analytical results show that: (i) The ‘wet’ frequencies of the structure are lower than the corresponding ‘dry’ (without water) frequencies due to the participating water mass, but the ‘wet’ mode shapes are identical to the corresponding ‘dry’ ones. (ii) The water mass plays an important role in the vertical seismic response, which varies with the different geological sites. For the different seismic inputs, the deeper the water is, the greater are the structural responses. (iii) The vertical seismic effects on the beam are generally not too small to be neglected and should be considered in the structural designs of a beam‐supported aqueduct. Copyright © 2002 John Wiley & Sons, Ltd.     

A coupled model for simulating surface water and groundwater interactions in coastal wetlands

Coastal wetlands are characterized by strong, dynamic interactions between surface water and groundwater. This paper presents a coupled model that simulates interacting surface water and groundwater flow and solute transport processes in these wetlands. The coupled model is based on two existing (sub) models for surface water and groundwater, respectively: ELCIRC (a three‐dimensional (3‐D) finite‐volume/finite‐difference model for simulating shallow water flow and solute transport in rivers, estuaries and coastal seas) and SUTRA (a 3‐D finite‐element/finite‐difference model for simulating variably saturated, variable‐density fluid flow and solute transport in porous media). Both submodels, using compatible unstructured meshes, are coupled spatially at the common interface between the surface water and groundwater bodies. The surface water level and solute concentrations computed by the ELCIRC model are used to determine the boundary conditions of the SUTRA‐based groundwater model at the interface. In turn, the groundwater model provides water and solute fluxes as inputs for the continuity equations of surface water flow and solute transport to account for the mass exchange across the interface. Additionally, flux from the seepage face was routed instantaneously to the nearest surface water cell according to the local sediment surface slope. With an external coupling approach, these two submodels run in parallel using time steps of different sizes. The time step (Δt_g) for the groundwater model is set to be larger than that (Δt_s) used by the surface water model for computational efficiency: Δt_g = M × Δt_s where M is an integer greater than 1. Data exchange takes place between the two submodels through a common database at synchronized times (e.g. end of each Δt_g). The coupled model was validated against two previously reported experiments on surface water and groundwater interactions in coastal lagoons. The results suggest that the model represents well the interacting surface water and groundwater flow and solute transport processes in the lagoons. Copyright © 2011 John Wiley & Sons, Ltd.     

Empirical model for estimating Fourier amplitude spectra of ground acceleration in Taiwan region

A collection of ground‐motion recordings (1070 acceleration records) of moderate (5.1⩽M_L⩽6.5) earthquakes obtained during the execution of the Taiwan Strong Motion Instrumentation Program (TSMIP) since 1991 was used to study source scaling model and attenuation relations for a wide range of earthquake magnitudes and distances and to verify the models developed recently for the Taiwan region. The results of the analysis reveal that the acceleration spectra of the most significant part of the records, starting from S‐wave arrival, can be modelled accurately using the Brune's ω‐squared source model with magnitude‐dependent stress parameter Δσ, that should be determined using the recently proposed regional relationships between magnitude (M_L) and seismic moment (M₀) and between M₀ and Δσ. The anelastic attenuation Q of spectral amplitudes with distance may be described as Q=225 ƒ^1.1 both for deep (depth more than 35 km) and shallow earthquakes. The source scaling and attenuation models allow a satisfactory prediction of the peak ground acceleration for magnitudes 5.1⩽M⩽6.5 and distances up to about 200 km in the Taiwan region, and may be useful for seismic hazard assessment. Copyright © 2000 John Wiley & Sons, Ltd.     

大型矩形渡槽-水耦合体系的动力性能分析

本文把任意拉格朗日-欧拉(ALE)方法用于分析大型矩形渡槽结构-水耦合体动力性能,并针对不同的水位、不同截面深宽比,研究渡槽结构在谐波激励、地震激励下的振动反应。研究表明,渡槽结构中水体的泼溅具有明显的非线性特征;渡槽中水位对耦合体动力特性影响极具规律性,但是随截面深宽比减小,水体振荡的非线性明显增强;渡槽结构截面深宽比选择合理,水体泼溅作用产生的倾覆力矩不会使槽体发生倾覆失稳,但是水体的大幅晃动产生的动水压力对渡槽槽身应力有显著的影响,在渡槽抗震设计中应给予足够重视。     

Seismic risk of circum-pacific earthquakes I. Strain energy release

Commonly used earthquake “whole process” frequency - magnitude and strain energy - magnitude laws are merged to obtain an analytic expression for an upper bound magnitude to regional earthquake occurrenceM ₃, which is expressed primarily in terms of the annual maximum magnitudeM ₁ and the magnitude equivalent of the annual average total strain energy releaseM ₂. Values ofM ₃ are also estimated graphically from cumulative strain energy release diagrams. Both methods are illustrated by application to the high seismicity of the circum-Pacific belt, using Duda’s (1965) data and regionalisation. Values ofM ₃ obtained analytically, with their uncertainties, are in agreement with those obtained graphically. Empirical relations are then obtained betweenM ₁,M ₂, andM ₃, which could be of general assistance in regional seismic risk considerations if they are found to be of a universal nature. For instance.M ₃ andM ₂ differ by one magnitude unit in subregions of the circum-Pacific.     

Multi‐variable relations for soil effects on seismic ground motion

Soil effects on peak ground acceleration, velocity and elastic response spectra (5% damping) are expressed by simple approximate relations in terms of five key parameters: (a) the fundamental vibration period of the non‐linear soil, T_S, (b) the period of a bedrock site of equal thickness, T_b, (c) the predominant excitation period, Te, (d) the peak seismic acceleration at outcropping bedrock, a, and (e) the number of significant excitation cycles, n. Furthermore, another relation is proposed for the estimation of T_S in terms of the soil thickness H, the average shear wave velocity of the soil V?_S,o and a. The aforementioned parameters were first identified through a simplified analytical simulation of the site excitation. The multivariable approximate relations were then formulated via a statistical analysis of relevant data from more than 700 one‐dimensional equivalent‐linear seismic ground response analyses, for actual seismic excitations and natural soil conditions. Use of these relations to back‐calculate the numerical results in the database gives an estimate of their error margin, which is found to be relatively small and unbiased. The proposed relations are also independently verified through a detailed comparison with strong motion recordings from seven well‐documented case studies: (a) two sites in the San Fernando valley during the Northridge earthquake, and (b) five different seismic events recorded at the SMART‐1 accelerometer array in Taiwan. It is deduced that the accuracy of the relations is comparable to that of the equivalent‐linear method. Hence, they can be readily used as a quick alternative for routine applications, as well as for spreadsheet computations (e.g. GIS‐aided seismic microzonation studies) where numerical methods are cumbersome to implement. Copyright © 2003 John Wiley & Sons, Ltd.     

Hysteresis‐based analysis of overland metal transport

Introducing a concept of equivalent mass depth of flow, this study describes the phenomenon of non‐point source pollutant (metal) transport for pavement (or overland) flow in analogy with wave propagation in wide open channels. Hysteretic and normal mass rating curves are developed for runoff rate and mass of 12 dissolved and particulate‐bound metal elements (pollutants) using the rainfall‐runoff and water quality data of the 15 × 20 m² instrumented pavement in Cincinnati, USA. Normal mass rating curves developed for easy computation of pollutant load are found to be of a form similar to Manning's, which is valid for open channel flows. Based on the hysteresis analysis, wave types for dissolution and mixing of particulate‐bound metals are identified. The analysis finds that the second‐order partial‐differential equation normally used for metal transport does not have the efficacy to describe fully the strong non‐linear phenomena such as is described for various metal elements by dynamic waves. In addition, the proportionality concept of the popular SCS‐CN concept is extended for determining the potential maximum metal mass M_p of all the 12 elements transported by a rain storm and related to the antecedent dry period (ADP). For the primary metal zinc element, M_p is found to increase with the ADP and vice versa. Copyright © 2003 John Wiley & Sons, Ltd.     

An insightful parametrization for the flatlander's interpretation of time‐lapsed seismic data

An approximation is developed that allows mapped 4D seismic amplitudes and time‐shifts to be related directly to the weighted linear sum of pore pressure and saturation changes. The weights in this relation are identified as key groups of parameters from a petroelastic model and include the reservoir porosity. This dependence on groups of parameters explains the inherent non‐uniqueness of this problem experienced by previous researchers. The proposed relation is of use in 4D seismic data feasibility studies and inversion and interpretation of the 4D seismic response in terms of pore pressure and water saturation changes. A further result is drawn from analysis of data from the North Sea and West Africa, which reveals that the relative interplay between the effects of pore pressure and saturation changes on the seismic data can be simplified to the control of a single, spatially variant parameter C_S/C_P. Combining these results with those from published literature, we find that C_S/C_P = 8 appears to be a generality across a range of clastic reservoirs with a similar mean porosity. Using this C_S/C_P value, an in situ seismic‐scale constraint for the rock stress sensitivity component of the petroelastic model is constructed considering this component carries the largest uncertainty.     

Spectral magnitude and seismic radiated energy computed from CDSN records

We have employed 10 digital records and computed the spectral magnitude and the seismic radiated energy for 18 large earthquakes (M _s≥6) occurred in Eur-asian belt during 1986–1989. The nine digital stations (CD-SN) distribute all over China and one in Germany. The spectral magnitudes of various period have different stability among stations. The stability is better for maximum spectral magnitudemi and seismic radiated energyE, their differences among stations are smaller, especially for the stations where the ray path main penetrates the low mantle. But the stability of corner period is usually not good. The relation between seismic radiated energy and seismic moment magnitudeM _w is lg (E)=1.5Mw+c, wherec is a constant. The maximum spectral magnitudemi=M _w+0.1, it is consistant with theoretical prediction. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 418–426, 1993. This work supported by the Deutsche Forschungsgemeinschaft, Bonn, F. R. Germany. The support is grateful acknowledge.     

大型渡槽有限条法动力建模研究

根据渡槽结构特点，采用有限条划分渡槽槽身，推导出渡槽槽身结构的单元刚度矩阵、相容质量矩阵的显式表达式；渡槽支架采用空间梁单元模拟；联接槽身和支架的盆式橡胶支座采用弹性元件单元模拟。编程程序具体计算了某大型渡槽的模态，并用ANSYS对大型渡槽进行模态分析，计算结果表明：该模型计算的渡槽结构的固有频率和ANSYS计算结果接近。模型可用于大型渡槽的动力分析。     

Characteristics of coseismic water level changes at Tangshan well for the Wenchuan MS8.0 earthquake and its larger aftershocks

Coseismic water level changes which may have been induced by the Wenchuan M_S8.0 earthquake and its 15 larger aftershocks (M_S≥?5.4) have been observed at Tangshan well. We analyze the correlation between coseismic parameters (maximum amplitude, duration, coseismic step and the time when the coseismic reach its maximum amplitude) and earthquake parameters (magnitude, well-epicenter distance and depth), and then compare the time when the coseismic oscillation reaches its maximum amplitude with the seismogram from Douhe seismic station which is about 16.3 km away from Tangshan well. The analysis indicates that magnitude is the main factor influencing the induced coseismic water level changes, and that the well-epicenter distance and depth have less influence. M_S magnitude has the strongest correlation with the coseismic water level changes comparing to M_W and M_L magnitudes. There exists strong correlation between the maximum amplitude, step size and the oscillation duration. The water level oscillation and step are both caused by dynamic strain sourcing from seismic waves. Most of the times when the oscillations reach their maximum amplitudes are between S and Rayleigh waves. The coseismic water level changes are due to the co-effect of seismic waves and hydro-geological environments.     

核电结构PCS水箱液动压力分析等效模型

中国核电厂抗震设计规范推荐采用的Housner模型不适用于复杂形状核电储液结构的流固耦合分析。对于AP1000和CAP1400核电站屏蔽厂房顶部非能动安全壳冷却系统重力水箱(简称PCS水箱),基于圆柱形水箱的Housner等效质量-弹簧模型,通过引入水箱体积修正参数,提出PCS水箱的三维等效质量-弹簧模型。采用有限元软件ADINA建立水箱结构流固耦合整体有限元模型以进行模态分析,计算分析PCS水箱和对应环形水箱在不同尺寸和液体深度条件下的液体晃动自振特性。对比整体有限元模型与三维等效质量-弹簧模型计算结果发现,提出的PCS水箱三维等效质量-弹簧模型能给出其内液体晃动各阶振型的液动压力合理估计值,适用于具有复杂形状的PCS水箱液动压力分析。本文的等效模型方法可推广应用于其他复杂形状水箱的液动压力分析。     

Seismic risk of circum-Pacific earthquakes: II. Extreme values using Gumbel's third distribution and the relationship with strain energy release

In a previous paper (Makropoulos andBurton, 1983) the seismic risk of the circum-Pacific belt was examined using a whole process technique reduced to three representative parameters related to the physical release of strain energy, these are:M ₁, the annual modal magnitude determined using the Gutenberg-Richter relationship;M ₂, the magnitude equivalent to the total strain energy release rate per annum, andM ₃, the upper bound magnitude equivalent to the maximum strain energy release in a region.The risk analysis is extended here using the part process statistical model of Gumbel's IIIrd asymptotic distribution of extreme values. The circum-Pacific is chosen being a complete earthquake data set, and the stability postulate on which asymptotic distributions of extremes are deduced to give similar results to those obtained from whole process or exact distributions of extremes is successfully checked. Additionally, when Gumbel III asymptotic distribution curve fitting is compared with Gumbel I using reduced chi-squared it is seen to be preferable in all cases and it also allows extensions to an upper-bounded range of magnitude occurrences. Examining the regional seismicity generates several seismic risk results, for example, the annual mode for all regions is greater thanm(1)=7.0, with the maximum being in the Japan, Kurile, Kamchatka region atm(1)=7.6. Overall, the most hazardous areas are situated in this northwestern region and also diagonally opposite in the southeastern circum-Pacific. Relationships are established between the Gumbel III parameters and quantitiesm ₁(1),X ₂ and , quantities notionally similar toM ₁,M ₂ andM ₃ although is shown to be systematically larger thanM; thereby giving a physical link through strain energy release to seismic risk statistics. Inall regions of the circum-Pacific similar results are obtained forM ₁,M ₂ andM ₃ and the notionally corresponding statistical quantitiesm ₁(1),X ₂ and , demonstrating that the relationships obtained are valid over a wide range of seismotectonic enviroments.     

Spectral responses of water level in tidal river and groundwater

The water level of five river stages and seven groundwater wells in the Taipei Basin were analysed by spectral analysis in the frequency domain. The diurnal, semi‐diurnal and quarter‐diurnal tidal components of the Tanshui River appear to relate closely to astronomical tides as K₁, M₂ and M₄, respectively. It is also found that the diurnal component reveals a reversed phase angle in the middle section of the Tanshui River; the phase of the quarter‐diurnal component is also found to be reversed at stations upstream in the Tanshui River and Hsintien Stream. It is believed that these phenomena could be caused by local variation in the river channel topography. The autospectrum and cross‐spectrum between groundwater elevation and nearby river stage were observed to correlate highly with the frequency of the astronomical tides K₁, M₂ and M₄. From the study of the phase shift and time lag of water level fluctuations at river stages and groundwater wells, it was found that the tidal effects of diurnal, semi‐diurnal, and quarter‐diurnal components were significantly different. The relationships between phase and the fluctuated range of atmospheric pressure and water level imply that change in atmospheric pressure does not affect water level fluctuation in the river stage and groundwater well. Copyright © 1999 John Wiley & Sons, Ltd.     

Investigation of seismic behavior of fluid–rectangular tank–soil/foundation systems in frequency domain

Main purpose of this study is to evaluate the dynamic behavior of fluid–rectangular tank–soil/foundation system with a simple and fast seismic analysis procedure. In this procedure, interaction effects are presented by Housner's two mass approximations for fluid and the cone model for soil/foundation system. This approach can determine; displacement at the height of the impulsive mass, the sloshing displacement and base forces for the soil/foundation system conditions including embedment and incompressible soil cases. Models and equations for proposed method were briefly explained for different tank–soil/foundation system combinations. By means of changing soil/foundation conditions, some comparisons are made on base forces and sloshing responses for the cases of embedment and no embedment. The results showed that the displacements and base shear forces generally decreased, with decreasing soil stiffness. However, embedment, wall flexibility, and soil–structure interaction (SSI) did not considerably affect the sloshing displacement.     

The regional characteristics of seismic activity in China

ＴｈｅｒｅｇｉｏｎａｌｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｓｅｉｓｍｉｃａｃｔｉｖｉｔｙｉｎＣｈｉｎａＺｈｅｎ－ＬｉａｎｇＳＨＩ，ＪｉａｎＷＡＮＧａｎｄＸｉａｏ－ＤｏｎｇＺＨＡＮＧ（时振梁，王健，张晓东）（ＩｎｓｔｉｔｕｔｅｏｆＧｅｏｐｈｙｓｉｃｓ，Ｓ...     

大型渡槽结构模态分析

大型渡槽在南水北调水利工程中起着重要作用。本文根据大型渡槽的结构特点,考虑渡槽横向弯扭耦合振动、约束扭转变形、渡槽支架变形和盆式橡胶支座等胶地支座等对渡槽动力作用的影响,应用梁段有限元方法,对某大型渡槽的模态进行分析,研究了该渡槽设计水位、渡槽支架高度、渡槽支架截面尺寸、渡槽盆式橡胶支座刚度变化等对大型渡槽模态的影响,探讨了该大型渡槽的上述４种工况下的模态变化范围,所得结果可为大型渡槽的抗震设计提     

Effect of boundary element details on the seismic deformation capacity of structural walls

The objective of this study is to investigate the effect of boundary element details of structural walls on their deformation capacities. Structural walls considered in this study have different sectional shapes and/or transverse reinforcement content at the boundaries of the walls (called boundary element details hereafter). Four full‐scale wall specimens (3000mm (h_w)×1500mm (l_w)×200mm (T)) were fabricated and tested. Three specimens are rectangular in section and the other specimen has a barbell‐shaped cross‐section (a wall with boundary columns). The rectangular wall specimens are reinforced according to the common practice used for reinforced concrete residence buildings in Korea and Chile. In this study, the primary variable for these rectangular specimens is the content of transverse reinforcement to confine the boundary elements of a wall. The barbell‐shaped specimen was designed in compliance with ACI 318‐95. The response of the barbell‐shaped specimen is compared with those of other rectangular specimens. The effective aspect ratio of the specimens is set to two in this study. Based on the experimental results, it is found that the deformation capacities of walls, which are represented by displacement ductility, drift ratio and energy dissipation capacities, are affected by the boundary element details. Copyright © 2002 John Wiley & Sons, Ltd.