首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
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.  相似文献   

2.
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.  相似文献   

3.
We investigate the dependence of quasi P-wave phase velocity propagating in orthotropic media on particular elasticity parameters. Specifically, due to mathematical facilitation, we consider the squared-velocity difference, , resulted from propagation in two mutually perpendicular symmetry planes. In the context of the effective medium theory, may be viewed as a parameter evaluating the influence of cracks – embedded in the background medium – parallel to one or both aforementioned planes. Our investigation is both theoretical and numerical. Based on Christoffel's equations, we propose two accurate approximations of . Due to them, we interpret the aforementioned squared-velocity difference as being twice more dependent on , than on . To describe the magnitude of the dependence, we consider the proportions between the partial derivatives of . Further, it occurs that is influenced by the ratio of vertically propagating quasi P-wave to vertically propagating quasi S-wave. Anomalously high might be caused by the low P/S ratio, which in turn can be an indicator of the presence of gas in natural fractures or aligned porosity. Also, we carry out numerical sensitivity study, according to which is approximately twice more dependent on than on , twice more sensitive to than to , and equally dependent on as on . The dependence on and can be neglected, especially for small phase angles. We verify the approximations and perform the sensitivity study, using eight examples of the elasticity tensors.  相似文献   

4.
Full waveform inversion in transversely isotropic media with a vertical symmetry axis provides an opportunity to better match the data at the near and far offsets. However, multi-parameter full waveform inversion, in general, suffers from serious cycle-skipping and trade-off problems. Reflection waveform inversion can help us recover a background model by projecting the residuals of the reflected wavefield along the reflection wavepath. Thus, we extend reflection waveform inversion to acoustic transversely isotropic media with a vertical symmetry axis utilizing the proper parameterization for reduced parameter trade-off. From a radiation patterns analysis, an acoustic transversely isotropic media with a vertical symmetry axis is better described by a combination of the normal-moveout velocity and the anisotropic parameters η and δ for reflection waveform inversion applications. We design a three-stage inversion strategy to construct the optimal resulting model. In the first stage, we only invert for the background by matching the simulated reflected wavefield from the perturbations of and δ with the observed reflected wavefield. In the second stage, the background and η are optimized simultaneously and the far-offset reflected wavefield mainly contribute to their updates. We perform Born modelling to compute the reflected wavefield for the two stages of reflection waveform inversion. In the third stage, we perform full waveform inversion for the acoustic transversely isotropic media with a vertical symmetry axis to delineate the high-wavenumber structures. For this stage, the medium is described by a combination of the horizontal velocity , η and ε instead of , η and δ. The acoustic multi-parameter full waveform inversion utilizes the diving waves to improve the background as well as utilizes reflection for high-resolution information. Finally, we test our inversion algorithm on the modified Sigsbee 2A model (a salt free part) and a two-dimensional line from a three-dimensional ocean bottom cable dataset. The results demonstrate that the proposed reflection waveform inversion approach can recover the background model for acoustic transversely isotropic media with a vertical symmetry axis starting from an isotropic model. This recovered background model can mitigate the cycle skipping of full waveform inversion and help the inversion recover higher resolution structures.  相似文献   

5.
Ebb-tidal deltas are highly dynamic environments affected by both waves and currents that approach the coast under various angles. Among other bedforms of various scales, these hydrodynamics create small-scale bedforms (ripples), which increase the bed roughness and will therefore affect hydrodynamics and sediment transport. In morphodynamic models, sediment transport predictions depend on the roughness height, but the accuracy of these predictors has not been tested for field conditions with strongly mixed (wave–current dominated) forcing. In this study, small-scale bedforms were observed in the field with a 3D Profiling Sonar at five locations on the Ameland ebb-tidal delta, the Netherlands. Hydrodynamic conditions ranged from wave dominated to current dominated, but were mixed most of the time. Small-scale ripples were found on all studied parts of the delta, superimposed on megaripples. Even though a large range of hydrodynamic conditions was encountered, the spatio-temporal variations in small-scale ripple dimensions were relatively small (height 0.015 m, length 0.11 m). Also, the ripples were always highly three-dimensional. These small dimensions are probably caused by the fact that the bed consists of relatively fine sediment. Five bedform height predictors were tested, but they all overestimated the ripple heights, partly because they were not created for small grain sizes. Furthermore, the predictors all have a strong dependence on wave- and current-related velocities, whereas the ripple heights measured here were only related to the near-bed orbital velocity. Therefore, ripple heights and lengths in wave–current-dominated, fine-grained coastal areas ( mm) may be best estimated by constant values rather than values dependent on the hydrodynamics. In the case of the Ameland ebb-tidal delta, these values were found to be m and m. ©2019 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.  相似文献   

6.
Laboratory experiments were conducted to investigate the formation of river bedforms under sediment supply-limited conditions, i.e. when a motionless substratum is bared by the dynamics of the mobile sediments. Three series of experiments were organized in a laboratory flume by fixing all the hydrodynamic and morphodynamic parameters but varying the thickness of the initial layer of mobile sediments which covers the rigid bottom of the flume. At the end of all the experiments, which lasted for the same amount of time, the formation of transverse sand dunes was observed. For decreasing , the rigid bottom of the flume was bared progressively earlier during the experiment and the measurements showed a clear tendency of the bedforms to lengthen, i.e. to increase their crest-to-crest distance. Moreover, under strong supply limitation, the two-dimensional transverse dunes turned into three-dimensional barchanoid forms and into isolated barchan dunes characterized by an abrupt reduction in bedform heights. A two-dimensional Fourier analysis of the bottom profile was performed, providing the amplitude of the main streamwise and spanwise harmonic components of the bottom morphology as a function of . © 2019 John Wiley & Sons, Ltd.  相似文献   

7.
To more accurately predict the migration behavior of pollutants in porous media, we conduct laboratory scale experiments and model simulation. Aniline (AN) is used in one-dimensional soil column experiments designed under various media and hydrodynamic conditions. The advection-dispersion equation (ADE) and the continuous-time random walk (CTRW) were used to simulate the breakthrough curves (BTCs) of the solute transport. The results show that the media and hydrodynamic conditions are two important factors affecting solute transport and are related to the degree of non-Fickian transport. The simulation results show that CTRW can more effectively describe the non-Fickian phenomenon in the solute transport process than ADE. The sensitive parameter in the CTRW simulation process is , which can reflect the degree of non-Fickian diffusion in the solute transport. Understanding the relationship of with velocity and media particle size is conducive to improving the reactive solute transport model. The results of this study provide a theoretical basis for better prediction of pollutant transport in groundwater.  相似文献   

8.
Logistic regression, neural networks and support vector machines are tested for their effectiveness in isolating surface waves in seismic shot records. To distinguish surface waves from other arrivals, we train the algorithms on three distinguishing features of surface-wave dispersion curves in the domain: spectrum coherency of the trace's magnitude spectrum, local dip and the frequency range for a fixed wavenumber k in the spectrum. Numerical tests on synthetic data show that the kernel-based support vector machines algorithm gives the highest accuracy in predicting the surface-wave window in the domain compared to neural networks and logistic regression. This window is also used to automatically pick the fundamental dispersion curve. The other two methods correctly pick the low-frequency part of the dispersion curve but fail at higher frequencies where there is interference with higher-order modes.  相似文献   

9.
Probabilistic seismic analysis of structures involves the construction of seismic demand models, often stated as probabilistic models of structural response conditioned on a seismic intensity measure. The uncertainty introduced by the model is often a result of the chosen intensity measure. This paper introduces the concept of using fractional order intensity measures (IMs) in probabilistic seismic demand analysis and uses a single frame integral concrete box‐girder bridge class and a seismically designed multispan continuous steel girder bridge class as case studies. The fractional order IMs considered include peak ground response and spectral accelerations at 0.2 and 1.0 s considering a single degree of freedom system with fractional damping, , as well as a linear single degree of freedom system with fractional response, . The study reveals the advantage of fractional order IMs relative to conventional IMs such as peak ground acceleration, peak ground velocity, or spectral acceleration at 0.2 and 1.0 s. Metrics such as efficiency, sufficiency, practicality, and proficiency are measured to assess the optimal nature of fractional order IMs. The results indicate that the proposed fractional order IMs produce significant improvements in efficiency and proficiency, whereas maintaining practicality and sufficiency, and thus providing superior demand models that can be used in probabilistic seismic demand analysis. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

10.
The proportion of water younger than 2–3 months (young water fraction, Fyw) has become increasingly investigated in catchment hydrology. Fyw is typically estimated by comparing seasonal tracer cycles in precipitation and streamflow, through water sampling. However, some open research questions remain, such as: (i) whether part of the summer precipitation should be discarded because the high evapotranspiration demand, (ii) how well Fyw serves as a metric to compare catchments, and (iii) how sampling frequency affects Fyw estimates. To address these questions, we investigated Fyw in soil-, ground- and stream waters for the small Mediterranean Can Vila catchment. Rainfall was sampled at 5-mm intervals. Mobile soil water and groundwater were sampled fortnightly. Stream water was sampled depending on flow at variable time intervals (30 min to 1 week). Over 58 months, this sampling provided 1,529 δ18O determinations. Isotopic analyses results led us to include summer precipitation in the input signal. We found the highest Fyw in mobile soil waters (34%), while this was almost zero for groundwater except during wet periods. For stream waters, Fyw depended on the discharge variations, so that the flow-weighted young water fraction () was 22.6%, whereas the time-weighted Fyw was just 6.2%. Both and its discharge sensitivity (Sd) varied when different 12-month sampling periods were investigated. The young water fraction that would be obtained from a virtual thorough sampling () was estimated from the Sd and the observed stream flow. This showed an underestimation of by 25% for the frequent dynamic sampling and 66% for weekly sampling, due to missing high flows. Our results confirm that Fyw and its discharge sensitivity are metrics very sensitive to meteorological forcing during the analysed period. Thus, comparisons between catchments need long-term mean annual values and their variability. Our findings also support the dependence of Fyw estimates on the sampling rate and show the advantages of flow-weighted analysis. Finally, catchment water turnover investigations should be accompanied by the analysis of flow duration curves.  相似文献   

11.
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.  相似文献   

12.
Organic carbon (OC) in valley bottom downed wood and soil that cycles over short to moderate timescales (101 to 105 years) represents a large, dynamic, and poorly quantified pool of carbon whose distribution and residence time affects global climate. We sought to quantify this potentially important OC pool at the watershed scale to estimate its magnitude and age, as well as determine the controls on its variability within watersheds. To do this, we compared four disparate mountain river basins to show that mountain river valley bottoms store substantial estimated OC stocks in floodplain soil and downed wood (median OC of MgC/ha, n = 178). Although soil OC is generally young (exhibiting a median radiocarbon fraction modern value of , n = 121), geomorphic processes regulate soil burial and processes that limit microbial respiration, preserving aged OC in especially deep, unconfined, wet, and/or high-elevation floodplain soils. We statistically modeled OC stocks to show that valley bottom morphology and hydrology regulate variability in floodplain soil retention and resulting variability in OC stock and age in floodplain soil throughout river networks. Comparing the distribution of OC stocks between wood and soil, we find that where floodplain soils are present, their OC stocks are generally greater than OC stocks stored in wood. Our results suggest that although mountain rivers may accumulate large OC stocks relatively rapidly, those stocks are highly sensitive to alterations in soil and wood retention, implying that human alterations to either disturb or restore floodplain wood and soil storage may have substantial impacts on OC storage in river corridors. © 2020 John Wiley & Sons, Ltd.  相似文献   

13.
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.  相似文献   

14.
Stream channel morphology forms the template upon which hydraulic aspects of aquatic habitat are created, yet spatial and temporal variability in habitat imposed by changing morphology is not well understood. This paper presents a conceptual model linking sediment supply patterns to spatial and temporal variability in channel form and aquatic habitat. To evaluate this model, change over time in three habitat variables is quantified using a 2D hydrodynamic modeling approach. A 45-year record of topographic data from Carnation Creek, a catchment in coastal British Columbia, is used for the flow modeling. Using the Nays2DH modeling platform, water depths and velocities are simulated in eight channel segments located at different positions relative to locations of historical colluvial input using seven flow levels ranging from 3% to 400% of mean annual discharge (0.02 to 3.31 m s ). Results indicate that habitat availability changes through time as a result of sediment supply-driven changes to channel morphology and wood loads, but patterns in habitat vary as a function of dominant channel segment morphology. Spatial and temporal variability in morphology also influences the relationship between habitat availability and river discharge, leading to non-stationary habitat-discharge rating curves. When habitat areas are predicted by applying these curves to daily flow series spanning annual dry seasons, over 50% of the variance in cumulative seasonal habitat area can be explained by year-to-year changes in channel morphology and wood loading, indicating that changing morphology is an important factor for driving temporal habitat variability. This variance is related to the morphological variability of a channel segment, which in turn is associated with the segment position relative to zones of colluvial input. Collectively, these results suggest that variability in habitat is impacted by channel morphology, and can be evaluated partly on the basis of a channel's sediment supply regime. © 2019 John Wiley & Sons, Ltd.  相似文献   

15.
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.  相似文献   

16.
Contamination of water due to bromate is a severe health hazard. The aim of the present study was to remove bromate from water using a crosslinked polystyrene based strongly basic anion exchange resin De‐Acidite FF‐IP. Batch experiments were performed to study the influence of various experimental parameters such as effect of pH, contact time, temperature, and effect of competing anions on bromate removal by De‐Acidite FF‐IP resin. At optimum parameters, the removal rate of bromate was very fast and 90% removal took place in 5 min and equilibrium was established within 10 min. The presence of competitive anions reduced the bromate adsorption in the order of Cl? > F? > CO > SO > NO > PO. The practical utility of this resin has been demonstrated by removing bromate in some of the commercial bottled water from Saudi Arabia. The level of bromate was determined using a very sensitive, precise and rapid method based on ultra‐performance liquid chromatography‐tendem mass spectrometry (UPLC‐MS/MS).  相似文献   

17.
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.  相似文献   

18.
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, TS, (b) the period of a bedrock site of equal thickness, Tb, (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 TS 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.  相似文献   

19.
  相似文献   

20.
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.  相似文献   

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

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