首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 375 毫秒
1.
A simple relationship is proposed in this paper to construct damage‐based inelastic response spectra including the effect of ground motion duration that it can be used for damage control in seismic design of structures. This relation is established for three groups of ground motions with short‐duration, moderate‐duration, and long‐duration ranges. To develop the model, the duration effect is included in the cyclic ductility of structures by an energy‐based method, and then strength reduction factors are computed based on this modified ductility (named ). The strength reduction factors were calculated for 44 stiffness‐degrading oscillators having vibration periods between 0.05 and 4.0 s, four ultimate ductility capacities, and five damage levels subjected to 296 earthquake records. The results showed that ductility capacity, damage level, and ground motion duration are effective parameters in the energy dissipation of structures, which affect the spectra. The values of short‐period oscillators (e.g., low‐rise structures) under short‐duration records are generally greater than those under moderate‐duration and long‐duration records. Residual analysis has been made in terms of magnitude and distance to examine the validity of the proposed simple expression. Finally, the introduced spectra were compared with three previously published proposals. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

2.
Proposals are developed to update Tables 11.4‐1 and 11.4‐2 of Minimum Design Loads for Buildings and Other Structures published as American Society of Civil Engineers Structural Engineering Institute standard 7‐10 (ASCE/SEI 7–10). The updates are mean next generation attenuation (NGA) site coefficients inferred directly from the four NGA ground motion prediction equations used to derive the maximum considered earthquake response maps adopted in ASCE/SEI 7–10. Proposals include the recommendation to use straight‐line interpolation to infer site coefficients at intermediate values of (average shear velocity to 30‐m depth). The NGA coefficients are shown to agree well with adopted site coefficients at low levels of input motion (0.1 g) and those observed from the Loma Prieta earthquake. For higher levels of input motion, the majority of the adopted values are within the 95% epistemic‐uncertainty limits implied by the NGA estimates with the exceptions being the mid‐period site coefficient, Fv, for site class D and the short‐period coefficient, Fa, for site class C, both of which are slightly less than the corresponding 95% limit. The NGA data base shows that the median value of 913 m/s for site class B is more typical than 760 m/s as a value to characterize firm to hard rock sites as the uniform ground condition for future maximum considered earthquake response ground motion estimates. Future updates of NGA ground motion prediction equations can be incorporated easily into future adjustments of adopted site coefficients using procedures presented herein. Published 2014. This article is a U.S. Government work and is in the public domain in the USA. Earthquake Engineering & Structural Dynamics published by John Wiley & Sons Ltd.  相似文献   

3.
Correlation in spectral accelerations for earthquakes in Europe   总被引:1,自引:0,他引:1  
The shape of a uniform hazard spectrum has been criticized to be unrealistic for a site where the spectral ordinates of the uniform hazard spectrum at different periods are governed by different scenario events and conservative for long‐return‐period earthquake shaking. The conditional mean spectrum considering epsilon (CMS‐ε) takes into account the correlation of spectral demands (represented by values of ε) at different periods, to address these issues. This paper proposes new prediction models for the correlation coefficient of ε(T1) and ε(T2), a key component for developing a CMS, using Pan‐European earthquake records from a European ground motion database. Epsilon (ε) for each record is computed using the 2005 Ambraseys ground‐motion prediction equation. The model can be used to develop CMS for European sites, and it can be incorporated in the European seismic standards. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

4.
5.
We present deterministic ground motion simulations that account for the cyclic multiaxial response of sediments in the shallow crust. We use the Garner Valley in Southern California as a test case. The multiaxial constitutive model is based on the bounding surface plasticity theory in terms of total stress and is implemented in a high‐performance computing finite‐element parallel code. A major advantage of this model is the small number of free parameters that need to be calibrated given a shear modulus reduction curve and the ultimate soil strength. This, in turn, makes the model suitable for regional‐scale simulations, where geotechnical data in the shallow crust are scarce. In this paper, we first describe a series of numerical experiments designed to verify the model implementation. This is followed by a series of idealized large‐scale simulations in a 35 26 4.5 km domain that encompasses the Garner Valley downhole array site, which is an instrumented and well‐characterized site in Southern California. Material properties were extracted from the Southern California Earthquake Center Community velocity model, CVM‐S4.26, considering its optional geotechnical layer, while the modulus reduction curves and soil strength were selected empirically to constrain the nonlinear soil model parameters. Our nonlinear simulations suggest that peak ground displacements within the valley increase relative to the linear case, while peak ground accelerations can increase or decrease, depending on the frequency content of the excitation. The comparisons of our simulations against hybrid three‐dimensional–one‐dimensional site response analyses suggest the inadequacy of the latter to capture the complexity of fully three‐dimensional simulations.  相似文献   

6.
Scenario‐based earthquake simulations at regional scales hold the promise in advancing the state‐of‐the‐art in seismic risk assessment studies. In this study, a computational workflow is presented that combines (i) a broadband Green's function‐based fault‐rupture and ground motion simulation—herein carried out using the “UCSB (University of California at Santa Barbara) method”, (ii) a three‐dimensional physics‐based regional‐scale wave propagation simulation that is resolved at  Hz, and (iii) a local soil‐foundation‐structure finite element analysis model. These models are interfaced with each other using the domain reduction method. The innermost local model—implemented in ABAQUS—is additionally enveloped with perfectly matched layer boundaries that absorb outbound waves scattered by the structures contained within it. The intermediate wave propagation simulation is carried out using Hercules , which is an explicit time‐stepping finite element code that is developed and licensed by the CMU‐QUAKE group. The devised workflow is applied to a  km region on the European side of Istanbul, which was modeled using detailed soil stratigraphy data and realistic fault rupture properties, which are available from prior microzonation surveys and earthquake scenario studies. The innermost local model comprises a chevron‐braced steel frame building supported by a shallow foundation slab, which, in turn, rests atop a three‐dimensional soil domain. To demonstrate the utility of the workflow, results obtained using various simplified soil‐structure interaction analysis techniques are compared with those from the detailed direct model. While the aforementioned demonstration has a limited scope, the devised workflow can be used in a multitude of ways, for example, to examine the effects of shallow‐layer soil nonlinearities and surface topography, to devise site‐ and structure‐specific seismic fragilities, and for calibrating regional loss models, to name a few.  相似文献   

7.
The nonlinearity of the seismic amplitude‐variation‐with‐offset response is investigated with physical modelling data. Nonlinearity in amplitude‐variation‐with‐offset becomes important in the presence of large relative changes in acoustic and elastic medium properties. A procedure for pre‐processing physical modelling reflection data is enacted on the reflection from a water‐plexiglas boundary. The resulting picked and processed amplitudes are compared with the exact solutions of the plane‐wave Zoeppritz equations, as well as approximations that are first, second, and third order in , , and . In the low angle range of 0°–20°, the third‐order plane‐wave approximation is sufficient to capture the nonlinearity of the amplitude‐variation‐with‐offset response of a liquid‐solid boundary with , , and ρ contrasts of 1485–2745 m/s, 0–1380 m/s, and 1.00–1.19 gm/cc respectively, to an accuracy value of roughly 1%. This is in contrast to the linear Aki–Richards approximation, which is in error by as much as 25% in the same angle range. Even‐order nonlinear corrective terms are observed to be primarily involved in correcting the angle dependence of , whereas the odd‐order nonlinear terms are involved in determining the absolute amplitude‐variation‐with‐offset magnitudes.  相似文献   

8.
This study uses instrumented buildings and models of code‐based designed buildings to validate the results of previous studies that highlighted the need to revise the ASCE 7 Fp equation for designing nonstructural components (NSCs) through utilizing oversimplified linear and nonlinear models. The evaluation of floor response spectra of a large number of instrumented buildings illustrates that, unlike the ASCE 7 approach, the in‐structure and the component amplification factors are a function of the ratio of NSC period to the supporting building modal periods, the ground motion intensity, and the NSC location. It is also shown that the recorded ground motions at the base of instrumented buildings in most cases are significantly lower than design earthquake (DE) ground motions. Because ASCE 7 is meant to provide demands at a DE level, for a more reliable evaluation of the Fp equation, 2 representative archetype buildings are designed based on the ASCE 7‐16 seismic provisions and exposed to various ground motion intensity levels (including those consistent with the ones experienced by instrumented buildings and the DE). Simulation results of the archetype buildings, consistent with previous numerical studies, illustrate the tendency of the ASCE 7 in‐structure amplification factor, [1 + 2(z/h)] , to significantly overestimate demands at all floor levels and the ASCE 7 limit of to in many cases underestimate the calculated NSC amplification factors. Furthermore, the product of these 2 amplification factors (that represents the normalized peak NSC acceleration) in some cases exceeds the ASCE 7 equation by a factor up to 1.50.  相似文献   

9.
Pre‐stack seismic data are indicative of subsurface elastic properties within the amplitude versus offset characteristic and can be used to detect elastic rock property changes caused by injection. We perform time‐lapse pre‐stack 3‐D seismic data analysis for monitoring sequestration at Cranfield. The time‐lapse amplitude differences of Cranfield datasets are found entangled with time‐shifts. To disentangle these two characters, we apply a local‐correlation‐based warping method to register the time‐lapse pre‐stack datasets, which can effectively separate the time‐shift from the time‐lapse seismic amplitude difference without changing the original amplitudes. We demonstrate the effectiveness of our registration method by evaluating the inverted elastic properties. These inverted time‐lapse elastic properties can be reliably used for monitoring plumes.  相似文献   

10.
Seismic intensity measures (IMs) perform a pivotal role in probabilistic seismic demand modeling. Many studies investigated appropriate IMs for structures without considering soil liquefaction potential. In particular, optimal IMs for probabilistic seismic demand modeling of bridges in liquefied and laterally spreading ground are not comprehensively studied. In this paper, a coupled-bridge-soil-foundation model is adopted to perform an in-depth investigation of optimal IMs among 26 IMs found in the literature. Uncertainties in structural and geotechnical material properties and geometric parameters of bridges are considered in the model to produce comprehensive scenarios. Metrics such as efficiency, practicality, proficiency, sufficiency and hazard computability are assessed for different demand parameters. Moreover, an information theory based approach is adopted to evaluate the relative sufficiency among the studied IMs. Results indicate the superiority of velocity-related IMs compared to acceleration, displacement and time-related ones. In particular, Housner spectrum intensity (HI), spectral acceleration at 2.0 s (S a-20), peak ground velocity (PGV), cumulative absolute velocity (CAV) and its modified version (CAV 5) are the optimal IMs. Conversely, Arias intensity (I a ) and shaking intensity rate (SIR) which are measures often used in liquefaction evaluation or related structural demand assessment demonstrate very low correlations with the demand parameters. Besides, the geometric parameters do not evidently affect the choice of optimal IMs. In addition, the information theory based sufficiency ranking of IMs shows an identical result to that with the correlation measure based on coefficient of determination (R 2). This means that R 2 can be used to preliminarily assess the relative sufficiency of IMs.  相似文献   

11.
Evapotranspiration (ET) can cause diel fluctuations in the elevation of the water table and the stage in adjacent streams. The diel fluctuations of water levels change head gradients throughout the day, causing specific discharge through near‐stream sediment to fluctuate at the same time scale. In a previous study, we showed that specific discharge controls the residence time of groundwater in streambed sediment that, in turn, exerted the primary control on removal from groundwater passing through the streambed. In this study, we examine the magnitude of diel specific discharge patterns through the streambed driven by ET in the riparian zone with a transient numerical saturated–unsaturated groundwater flow model. On the basis of a first‐order kinetic model for removal, we predicted diel fluctuations in stream concentrations. Model results indicated that ET drove a diel pattern in specific discharge through the streambed and riparian zone (the removal zones). Because specific discharge is inversely proportional to groundwater travel time through the removal zones and travel time determines the extent of removal, diel changes in ET can result in a diel pattern in concentration in the stream. The model predictions generally matched observations made during summertime base‐flow conditions in a small coastal plain stream in Virginia. A more complicated pattern was observed following a seasonal drawdown period, where source components to the stream changed during the receding limb of the hydrograph and resulted in diel fluctuations being superimposed over a multi‐day trend in concentrations. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

12.
Release of nitrogen compounds into groundwater, particularly those compounds from excessive agricultural fertilization, is a major concern in an aquifer recharge. Among the nitrogen compounds, ammonium ( ) is a common one. In order to assess the risk of agricultural fertilizer contamination to an aquifer through infiltration, adsorption onto a loamy agricultural soil profile (0–0.60 m depth) was studied using a soil column experiment and modelling simulation. The soil used in the experiment was drawn from an agricultural field in Xinzhen, Fangshan district, Beijing, China, and reconstituted in laboratory soil columns. Column experiments were conducted using bromide (conservative tracer) and ‐bearing aqueous solutions. The ammonium concentrations in the soil water samples were measured, and their values were plotted as the breakthrough curves. The chemical's soil–water distribution coefficients (Kd) were calculated using breakthrough curves. Then the retardation factor (R) in saturated soil was calculated. For the ‐bearing aqueous solutions, the strongest adsorption occurred at the soil depth of 0.30–0.45 m. The convection–dispersion equation model and chemical non‐equilibrium model in Hydrus‐1D were used to simulate transport in the loamy soil. The two‐site chemical non‐equilibrium model in Hydrus‐1D was best to simulate transport through the soil column. Parameter sensitivity study was conducted to investigate the influences of solute transport by Kd, the fraction of exchange sites assuming to be in equilibrium with the solution phase (f), the longitudinal dispersivity (λ), and the first‐order rate coefficients (ω). The sensitivity analysis results indicate Kd is the most critical parameter.  相似文献   

13.
The conventional velocity scan can be computationally expensive for large‐scale seismic data sets, particularly when the presence of anisotropy requires multiparameter scanning. We introduce a fast algorithm for 3D azimuthally anisotropic velocity scan by generalizing the previously proposed 2D butterfly algorithm for hyperbolic Radon transforms. To compute semblance in a two‐parameter residual moveout domain, the numerical complexity of our algorithm is roughly as opposed to of the straightforward velocity scan, with N being the representative of the number of points in a particular dimension of either data space or parameter space. Synthetic and field data examples demonstrate the superior efficiency of the proposed algorithm.  相似文献   

14.
Although changes in rainfall characteristics have been noted across the world, few studies have reported those in mountainous areas. This study was undertaken to clarify spatial and temporal variations in rainfall characteristics such as annual rainfall amount (Pr), mean daily rainfall intensity (η), and ratio of rain days (λ) in mountainous and lowland areas in Taiwan. To this aim, we examined spatial and year‐to‐year variations and marginal long‐term trends in Pr, η, and λ, based on rainfall data from 120 stations during the period 1978–2008. The period mean rainfall () at the lowland stations had strong relationships with the period mean daily rainfall intensity () and the period mean ratio of rain days () during those 31 years. Meanwhile, was only strongly related to at mountainous stations, indicating that influences on spatial variations in were different between lowland and mountainous stations. Year‐to‐year variations in Pr at each station were primarily determined from the variation in η at most stations for both lowland and mountainous stations. Long‐term trend analysis showed that Pr and η increased significantly at 10% and 31% of the total 120 stations, respectively, and λ decreased significantly at 6% of the total. The increases in Pr were mostly accompanied by increases in η. Although stations with significant η increases were slightly biased toward the western lowland area, increases or decreases in Pr and λ were not common. These results contribute to understanding the impacts of possible climate changes on terrestrial hydrological cycles. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

15.
The reassignment method remaps the energy of each point in a time‐frequency spectrum to a new coordinate that is closer to the actual time‐frequency location. Two applications of the reassignment method are developed in this paper. We first describe time‐frequency reassignment as a tool for spectral decomposition. The reassignment method helps to generate more clear frequency slices of layers and therefore, it facilitates the interpretation of thin layers. The second application is to seismic data de‐noising. Through thresholding in the reassigned domain rather than in the Gabor domain, random noise is more easily attenuated since seismic events are more compactly represented with a relatively larger energy than the noise. A reconstruction process that permits the recovery of seismic data from a reassigned time‐frequency spectrum is developed. Two approaches of the reassignment method are used in this paper, one of which is referred to as the trace by trace time reassignment that is mainly used for seismic spectral decomposition and another that is the spatial reassignment that is mainly used for seismic de‐noising. Synthetic examples and two field data examples are used to test the proposed method. For comparison, the Gabor transform method, inversion‐based method and common deconvolution method are also used in the examples.  相似文献   

16.
While seismic reflection amplitudes are generally determined by real acoustical impedance contrasts, there has been recent interest in reflections due to contrasts in seismic‐Q. Herein we compare theoretical and modelled seismic reflection amplitudes for two different cases of material contrasts. In case A, we examine reflections from material interfaces that have a large contrast in real‐valued impedance () with virtually no contrast in seismic‐Q. In case B, we examine reflections from material interfaces that have virtually no contrast in but that have very large seismic‐Q contrasts. The complex‐valued reflection coefficient formula predicts non‐zero seismic reflection amplitudes for both cases. We choose physical materials that typify the physics of both case A and case B. Physical modelling experiments show significantly large reflections for both cases – with the reflections in the two cases being phase shifted with respect to each other, as predicted theoretically. While these modelling experiments show the existence of reflections that are predicted by theory, there are still intriguing questions regarding the size of the Q‐contrast reflections, the existence of large Q‐contrast reflections in reservoir rocks and the possible application of Q‐reflection analysis to viscosity estimation in heavy oilfields.  相似文献   

17.
We present an approach based on local‐slope estimation for the separation of scattered surface waves from reflected body waves. The direct and scattered surface waves contain a significant amount of seismic energy. They present great challenges in land seismic data acquisition and processing, particularly in arid regions with complex near‐surface heterogeneities (e.g., dry river beds, wadis/large escarpments, and karst features). The near‐surface scattered body‐to‐surface waves, which have comparable amplitudes to reflections, can mask the seismic reflections. These difficulties, added to large amplitude direct and back‐scattered surface (Rayleigh) waves, create a major reduction in signal‐to‐noise ratio and degrade the final sub‐surface image quality. Removal of these waves can be difficult using conventional filtering methods, such as an filter, without distorting the reflected signal. The filtering algorithm we present is based on predicting the spatially varying slope of the noise, using steerable filters, and separating the signal and noise components by applying a directional nonlinear filter oriented toward the noise direction to predict the noise and then subtract it from the data. The slope estimation step using steerable filters is very efficient. It requires only a linear combination of a set of basis filters at fixed orientation to synthesize an image filtered at an arbitrary orientation. We apply our filtering approach to simulated data as well as to seismic data recorded in the field to suppress the scattered surface waves from reflected body waves, and we demonstrate its superiority over conventional techniques in signal preservation and noise suppression.  相似文献   

18.
Our objective was to discover the effect of variations in fluid properties and fracture geometry on the velocity of seismic wave propagation in fluid‐saturated media with parallel planar fractures. We used numerical models calculated by analytical solutions to examine the behaviour of P‐wave phase velocity dispersion in the normal direction to layering, in non‐porous and porous media with planar fractures. We also examined the anisotropy of low frequency phase and group velocities of fast and slow P‐waves and angular‐dependent reflection coefficients in media with planar fractures, under conditions of saturation by fluids with varying bulk moduli, densities, and fracture apertures. We defined several parametre , , and characterising dispersion, characterising anisotropy, characterising the difference between fast and slow modes, and R0 and characterising reflection. Our results show that the behaviour of dispersion shows wider stopbands in the case of gas saturation. Concavity indicator of dispersion for gas saturation was greater than that for liquid saturation and is usually greater than one. Anisotropy is more sensitive to bulk modulus contrast than to density contrast between the solid and the fluid, and is more sensitive to density contrast than to bulk modulus contrast. The case of gas saturation usually had a greater negative R0 and a greater value of compared with those of brine and heavy and light oil saturations. Our results are helpful in distinguishing fluid types saturating geophysical fractures and estimating the aperture and spacing of planar fractures. In seismic exploration, bulk modulus and fluid density can provide useful information in distinguishing among brine, oil, and gas; fracture geometry is important to estimate the permeability of reservoirs.  相似文献   

19.
Under certain circumstances, seismic propagation within porous media may be associated to the conversion of mechanical energy to electromagnetic energy, which is known as a seismo‐electromagnetic phenomenon. The propagation of fast compressional P‐waves is more specifically associated to the manifestations of a seismoelectric field linked to the fluid flows within the pores. The analysis of seismoelectric phenomena, which requires the combination of the theory of electrokinetics and Biot's theory of poroelasticity, provides us with transfer function that links the coseismic seismoelectric field E to the seismic acceleration . To measure the transfer function, we have developed an experimental setup enabling seismoelectric laboratory observation in unconsolidated quartz sand within the kilohertz range. The investigation focused on the impact of fluid conductivity and water saturation over the coseismic seismoelectric field. During the experiment, special attention was given to the accuracy of electric field measurements. We concluded that, to obtain a reliable estimate of the electric field amplitude, the dipole from which the potential differences are measured should be of much smaller length than the wavelength of the propagating seismic field. Time‐lapse monitoring of the seismic velocities and seismoelectric transfer functions were performed during imbibition and drainage experiments. In all cases, the quantitative analysis of the seismoelectric transfer function was in good agreement with theoretical predictions. While investigating saturation variations from full to residual water saturation, we showed that the ratio undergoes a switch in polarity at a particular saturation , which also implies a sign change of the filtration, traducing a reversal of the relative fluid displacement with respect to the frame. This sign change at critical saturation stresses a particular behaviour of the poroelastic medium: the dropping of the coseismic electric field to zero traduces the absence of relative pore/fluid displacements representative of a Biot dynamically compatible medium. We concluded from our experimental study in loose sand that the measurements of the coseismic seismoelectric coupling may provide information on fluid distribution within the pores and that the reversal of the seismoelectric field may be used as an indicator of the dynamically compatible state of the medium.  相似文献   

20.
Compressed Sensing has recently proved itself as a successful tool to help address the challenges of acquisition and processing seismic data sets. Compressed sensing shows that the information contained in sparse signals can be recovered accurately from a small number of linear measurements using a sparsity‐promoting regularization. This paper investigates two aspects of compressed sensing in seismic exploration: (i) using a general non‐convex regularizer instead of the conventional one‐norm minimization for sparsity promotion and (ii) using a frequency mask to additionally subsample the acquired traces in the frequency‐space () domain. The proposed non‐convex regularizer has better sparse recovery performance compared with one‐norm minimization and the additional frequency mask allows us to incorporate a priori information about the events contained in the wavefields into the reconstruction. For example, (i) seismic data are band‐limited; therefore one can use only a partial set of frequency coefficients in the range of reflections band, where the signal‐to‐noise ratio is high and spatial aliasing is low, to reconstruct the original wavefield, and (ii) low‐frequency characteristics of the coherent ground rolls allow direct elimination of them during reconstruction by disregarding the corresponding frequency coefficients (usually bellow 10 Hz) via a frequency mask. The results of this paper show that some challenges of reconstruction and denoising in seismic exploration can be addressed under a unified formulation. It is illustrated numerically that the compressed sensing performance for seismic data interpolation is improved significantly when an additional coherent subsampling is performed in the domain compared with the domain case. Numerical experiments from both simulated and real field data are included to illustrate the effectiveness of the presented method.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号