Effect of cumulative seismic damage and corrosion on the life‐cycle cost of reinforced concrete bridges

Bridge design should take into account not only safety and functionality, but also the cost effectiveness of investments throughout a bridge life‐cycle. This paper presents a probabilistic approach to compute the life‐cycle cost (LCC) of corroding reinforced concrete (RC) bridges in earthquake‐prone regions. The approach is developed by combining cumulative seismic damage and damage associated with corrosion due to environmental conditions. Cumulative seismic damage is obtained from a low‐cycle fatigue analysis. Chloride‐induced corrosion of steel reinforcement is computed based on Fick's second law of diffusion. The proposed methodology accounts for the uncertainties in the ground motion parameters, the distance from the source, the seismic demand on the bridge, and the corrosion initiation time. The statistics of the accumulated damage and the cost of repairs throughout the bridge life‐cycle are obtained by Monte‐Carlo simulation. As an illustration of the proposed approach, the effects of design parameters on the LCC of an example RC bridge are studied. The results are valuable in better estimating the condition of existing bridges and, therefore, can help to schedule inspection and maintenance programs. In addition, by taking into consideration the two deterioration processes over a bridge life‐cycle, it is possible to estimate the optimal design parameters by minimizing, for example, the expected cost throughout the life of the structure. A comparison between the effects of the two deterioration processes shows that, in seismic regions, the cumulative seismic damage affects the reliability of bridges over time more than the corrosion even for corrosive environments. Copyright © 2008 John Wiley & Sons, Ltd.     

Accidental eccentricity in symmetric buildings due to wave passage effects arising from near‐fault pulse‐like ground motions

This article investigates the characteristics of the accidental eccentricity in symmetric buildings due to torsional response arising from wave passage effects in the near‐fault region. The soil–foundation–structure system is modeled as a symmetric cylinder placed on a rigid circular foundation supported on an elastic halfspace and subjected to obliquely incident plane SH waves simulating the action of near‐fault pulse‐like ground motions. The translational response is computed assuming that the superstructure behaves as a shear beam under the action of translational and rocking base excitations, whereas the torsional response is calculated using the mathematical formulation proposed in a previous study. A broad range of properties of the soil–foundation–structure system and ground motion input are considered in the analysis, thus facilitating a detailed parametric investigation of the structural response. It is demonstrated that the normalized accidental eccentricity is most sensitive to the pulse period (T_P) of the near‐fault ground motions and to the uncoupled torsional‐to‐translational fundamental frequency ratio (Ω) of the structure. Furthermore, the normalized accidental eccentricities due to simplified pulse‐like and broadband ground motions in the near‐fault region are computed and compared against each other. The results show that the normalized accidental eccentricity due to the broadband ground motion is well approximated by the simplified pulse for longer period buildings, while it is underestimated for shorter period buildings. For symmetric buildings with values of Ω commonly used in design practice, the normalized accidental eccentricity due to wave passage effects is less than the typical code‐prescribed value of 5%, except for buildings with very large foundation radius. Copyright © 2017 John Wiley & Sons, Ltd.     

2012年江苏高邮、宝应交界MS4.9地震发震断层参数测定

基于一维单侧有限移动震源模式,根据地震波传播过程中的多普勒效应,分别利用P波和S波拐角频率的方位变化,反演2012年7月20日江苏高邮、宝应交界MS4.9地震的发震断层面参数。P波和S波拐角频率的反演结果一致显示:本次地震的断层面破裂方向为232°左右,破裂面呈NE-SW向;地震马赫数v/c为0.2左右,平均破裂速度小于S波速度,破裂长度较短,为0.2~0.3km左右。破裂面方位与震源机制解、宏观烈度调查和余震精定位的研究结果具有一致性,结合震区周边的地质构造背景,分析认为滁河断裂很可能是高邮、宝应交界MS4.9地震的发震构造。     

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.     

Rotation and strain seismology

Asymmetric continuum theory points to the equal roles of the rotation motions and those related directly to shear or confining strains. The strain motions could be quite independent or mutually related with the eventual phase shifts, while the displacements have only a mathematical sense; a real displacement may appear along the fracture slip only. Formally, any deformation could be presented as related to the displacements; however, its origin in a fracture source should be considered either as belonging to an individual process or to complex correlated events; in these cases, the confining, shear and rotation strains can be related mathematically to the different displacement fields. Some of these related deformations could be emitted from a source with a phase shift, while the observed displacements (deduced from records) result as a sum of these independent displacements. An important influence on a source process and on the premonitory micro-events has the material defects, their distribution, and mobility. The defect arrays lead to a concentration of stresses and their local reorganization. Thus, in this paper, we consider the induced stresses and strains related to defect content and to its modification and redistribution. Moreover, an important role in understanding the complex correlated events in a source plays the release?Crebound mechanism. The release?Crebound mechanism in an earthquake source processes leads to a possible direct or phase-shifted correlations between the emitted motions; in this aspect, a propagation of the coupled strain and rotation waves is discussed. In particular, we consider the point fracture events as related either to a confining load or/and to the shear and rotation processes; we discuss the related effects and their meaning when discussing the fault plane mechanism and emitted waves. In some important seismic regions, we have the recording system which permits to record the strains and rotations. However, we should point that the wordwide seismological network is not adequate to record the complex strain deformations released in the fracture processes and remains quite insufficient to understand the global stress changes and related strain waves of a very low period. Consideration on a recording mechanism of the long displacement waves indicates the insufficiencies of the present global recording system and points that recording of the global strain and rotation waves is an important and urgent task.     

Input and system identification of the Hualien soil–structure interaction system using earthquake response data

This paper presents an input and system identification technique for a soil–structure interaction system using earthquake response data. Identification is carried out on the Hualien large‐scale seismic test structure, which was built in Taiwan for international joint research. The identified quantities are the input ground acceleration as well as the shear wave velocities of the near‐field soil regions and Young's moduli of the shell sections of the structure. The earthquake response analysis on the soil–structure interaction system is carried out using the finite element method incorporating the infinite element formulation for the unbounded layered soil medium and the substructured wave input technique. The criterion function for the parameter estimation is constructed using the frequency response amplitude ratios of the earthquake responses measured at several points of the structure, so that the information on the input motion may be excluded. The constrained steepest descent method is employed to obtain the revised parameters. The simulated earthquake responses using the identified parameters and input ground motion show excellent agreement with the measured responses. Copyright © 2003 John Wiley & Sons, Ltd.     

Bivariate Gaussian mixture–based equivalent linearization method for stochastic seismic analysis of nonlinear structures

To address challenges in stochastic seismic analysis of nonlinear structures, this paper further develops a recently proposed Gaussian mixture–based equivalent linearization method (GM‐ELM). The GM‐ELM uses a Gaussian mixture distribution model to approximate the probabilistic distribution of a nonlinear system response. Using properties of the Gaussian mixture model, GM‐ELM can decompose the non‐Gaussian response of a nonlinear system into multiple Gaussian responses of linear single–degree of freedom oscillators. With the set of the equivalent linear systems identified by GM‐ELM, response statistics as crossing rate and first‐passage probability can be computed conveniently using theories of linear random vibration analysis. However, the original version of GM‐ELM may lead to an inaccurate estimate because of the heuristic parameters of the linear system introduced to supplement insufficient information. To overcome this limitation and define unique equivalent linear systems, this paper proposes a further developed version of GM‐ELM, which uses a mixture of bivariate Gaussian densities instead of univariate models. Moreover, to facilitate the use of elastic response spectra for estimating the mean peak responses of a nonlinear structure, a new response spectrum combination rule is proposed for GM‐ELM. Two numerical examples of hysteretic structural systems are presented in this paper to illustrate the application of the bivariate GM‐ELM to nonlinear stochastic seismic analysis. The analysis results obtained by the bivariate GM‐ELM are compared with those obtained by the univariate GM‐ELM, the conventional equivalent linearization method, the tail equivalent linearization method, and Monte Carlo simulation. The supporting source code and data are available for download at https://github.com/yisangri/GitHub‐bGM‐ELM‐code.git     

Joint elastic‐electrical properties of reservoir sandstones and their relationships with petrophysical parameters

We measured in the laboratory ultrasonic compressional and shear‐wave velocity and attenuation (0.7–1.0 MHz) and low‐frequency (2 Hz) electrical resistivity on 63 sandstone samples with a wide range of petrophysical properties to study the influence of reservoir porosity, permeability and clay content on the joint elastic‐electrical properties of reservoir sandstones. P‐ and S‐wave velocities were found to be linearly correlated with apparent electrical formation factor on a semi‐logarithmic scale for both clean and clay‐rich sandstones; P‐ and S‐wave attenuations showed a bell‐shaped correlation (partial for S‐waves) with apparent electrical formation factor. The joint elastic‐electrical properties provide a way to discriminate between sandstones with similar porosities but with different clay contents. The laboratory results can be used to estimate sandstone reservoir permeability from seismic velocity and apparent formation factor obtained from co‐located seismic and controlled source electromagnetic surveys.     

Basement structural architecture and hydrocarbon conduit potential of polygonal faults in the Cooper‐Eromanga Basin,Australia

A thorough and complete understanding of the structural geology and evolution of the Cooper‐Eromanga Basin has been hampered by low‐resolution seismic data that becomes particularly difficult to interpret below the thick Permian coal measures. As a result, researchers are tentative to interpret the basement fault architecture within the basin, which is largely undefined. To provide a better understanding of the basement fault geometry, all available two‐dimensional seismic lines together with 12 three‐dimensional seismic surveys were structurally interpreted with assistance from seismic attribute analysis. The Upper Cretaceous Cadna‐owie Formation and top Permian reflectors were analysed using a common seismic attribute technique (incoherency) that was used to infer the presence of faults that may have otherwise been overlooked. Detailed basement fault maps for each seismic survey were constructed and used in conjunction with two‐dimensional seismic data interpretation to produce a regional basement fault map. Large north‐northeast–south‐southwest‐striking sinistral strike–slip faults were identified within the Patchawarra Trough appearing to splay from the main northeast–southwest‐striking ridge. These sinistral north‐northeast–south‐southwest‐striking faults, together with field‐scale southeast–northwest‐striking dextral strike–slip faults, are optimally oriented to have potentially developed as a conjugated fault set under a south‐southeast–north‐northwest‐oriented strike–slip stress regime. Geomechanical modelling for a regionally extensive system of Cretaceous polygonal faults was performed to calculate the Leakage Factor and Dilation Tendency of individual faults. Faults that extend into Lower Cretaceous oil‐rich reservoirs with strikes of between 060°N and 140°N and a high to near‐vertical dip angle were identified to most likely be acting as conduits for the tertiary migration of hydrocarbons from known Lower Cretaceous hydrocarbon reservoirs into shallow Cretaceous sediments. This research provides valuable information on the regional basement fault architecture and a more detailed exploration target for the Cooper‐Eromanga Basin, which were previously not available in literature.     

Response of nonlinear soil‐MDOF structure systems subjected to distinct frequency‐content components of near‐fault ground motions

This paper is devoted to investigate the effects of near‐fault ground motions on the seismic responses of nonlinear MDOF structures considering soil‐structure interaction (SSI). Attempts are made to take into account the effects of different frequency‐content components of near‐fault records including pulse‐type (PT) and high‐frequency (HF) components via adopting an ensemble of 54 near‐fault ground motions. A deep sensitivity analysis is implemented based on the main parameters of the soil‐structure system. The soil is simulated based on the Cone model concept, and the superstructure is idealized as a nonlinear shear building. The results elucidate that SSI has approximately increasing and mitigating effects on structural responses to the PT and HF components, respectively. Also, a threshold period exists above which the HF component governs the structural responses. As the fundamental period of the structure becomes shorter and structural target ductility reduces, the contribution of the HF component to the structural responses increases, elaborately. Soil flexibility makes the threshold period increase, and the effect of the PT component becomes more significant than the HF one. In the case of soil‐structure system, slenderizing the structure also increases this threshold period and causes the PT component to be dominant. Copyright © 2013 John Wiley & Sons, Ltd.     

Risk assessment of an interacting structure–soil system due to liquefaction

The present study deals with the non‐linear stochastic dynamic analysis of a soil–structure interacting system. The ultimate objective is to determine the risk of damage to the system due to liquefaction under a wide range of earthquake intensities. A Monte Carlo simulation approach is followed in conjunction with advanced finite element procedures. The stochastic spatial variability of soil properties and the randomness of the seismic excitation are taken into account in order to estimate the statistics of the response, measured in terms of uniform foundation settlement and tilting. Specifically, soil properties are modelled as non‐Gaussian random fields and seismic excitations as non‐stationary random processes. The probabilistic characteristics of the stochastic field modelling soil properties are established from in situ tests. The risk of damage to the soil–structure system due to liquefaction is assessed by establishing fragility curves, which are of paramount importance for risk assessment and management studies of such systems. Fragility curves express the probability of exceeding various thresholds in the response. The relative effect of the variability of various soil parameters on the variability of the response is also examined. Copyright © 2002 John Wiley & Sons, Ltd.     

Seismic intensity distribution of shallow earthquakes due to rupture velocities and faulting modes

Isoseismals of seismic intensity distributions are represented by earthquake source size, faulting mode, and rupture velocity of fault propagation. Unilateral faulting forms egg—shaped isoseimals, while bilateral faulting forms elliptical ones. It is found that the ratio of major to minor axes of isoseismals is sensitive to rupture velocity. Rupture velocity, faulting mode, and fault trend have been determined from the seismic intensity maps of the 1964 Niigata and the 1983 Japan Sea earthquakes in China by matching theoretical isoseismals. Rupture velocities thus estimated are mostly 70% to 90% of shear wave velocity. The difference would be considered as follows: short—waves which determine the seismic intensity are strongly dependent on the rupture of small—scale fault heterogeneities and on the jerky—onsets and abrupt terminations of local rupture propagations. On the other hand, rupture velocity from long—waves represents an average rupture propagation along the whole fault length. Faulting mode and fault trend estimated from seismic intensity maps match with each earthquake faulting independently determined. This suggests that the present method would be applicable to some historical earthquakes with known seismic intensity distribution to obtain detailed information on the faulting process.     

Global-scale seismic interferometry: theory and numerical examples

Progress in the imaging of the mantle and core is partially limited by the sparse distribution of natural sources; the earthquake hypocenters are mainly along the active lithospheric plate boundaries. This problem can be approached with seismic interferometry. In recent years, there has been considerable progress in the development of seismic interferometric techniques. The term seismic interferometry refers to the principle of generating new seismic responses by cross‐correlating seismic observations at different receiver locations. The application of interferometric techniques on a global scale could create sources at locations where no earthquakes occur. In this way, yet unknown responses would become available for the application of travel‐time tomography and surface‐wave dispersion studies. The retrieval of a dense‐enough sampling of source gathers would largely benefit the application of reflection imaging. We derive new elastodynamic representation integrals for global‐scale seismic interferometry. The relations are different from other seismic interferometry relations for transient sources, in the sense that they are suited for a rotating closed system like the Earth. We use a correlation of an observed response with a response to which free‐surface multiple elimination has been applied to account for the closed system. Despite the fact that the rotation of the Earth breaks source‐receiver reciprocity, the seismic interferometry relations are shown to be valid. The Coriolis force is included without the need to evaluate an extra term. We synthesize global‐scale earthquake responses and use them to illustrate the acoustic versions of the new interferometric relations. When the sampling of real source locations is dense enough, then both the responses with and without free‐surface multiples are retrieved. When we do not take into account the responses from the sources in the direct neighborhood of the seismic interferometry‐constructed source location, the response with free‐surface multiples can still be retrieved. Even when only responses from sources at a certain range of epicentral distances are available, some events in the Green's function between two receiver locations can still be retrieved. The retrieved responses are not perfect, but the artefacts can largely be ascribed to numerical errors. The reconstruction of internal events – the response as if there was a source and a receiver on (major) contrasts within the model – could possibly be of use for imaging. With modelling it is possible to discover in which region of the correlation panel stationary phases occur that contribute to the retrieval of events. This knowledge opens up a new way of filtering out undesired events and of discovering whether specific events could be retrieved with a given source‐receiver configuration.     

Integrated foundation–structure seismic assessment through non‐linear dynamic analyses

This paper aims at clarifying the role of dynamic soil–structure interaction in the seismic assessment of structure and foundation, when the non‐linear coupling of both subsystems is accounted for. For this purpose, the seismic assessment of an ideal set of bridge piers on shallow foundations is considered. After an initial standard assessment, based on capacity design principles, the evaluation of the seismic response of the piers is carried out by dynamic simulations, where both the non‐linear responses of the superstructure and of the foundation are accounted for, in the latter case through the macro‐element modeling of the soil–foundation system. The results of the dynamic simulations point out the beneficial effects of the non‐linear response of the foundation, which provides a substantial contribution to the overall energy dissipation during seismic excitation, thus allowing the structural ductility demand to decrease significantly with respect to a standard fixed‐base or linear‐elastic base assessment. Permanent deformations at the foundation level, such as rotation and settlement, turn out to be of limited amount. Therefore, an advanced assessment approach of the integrated non‐linear system, consisting of the interacting foundation and superstructure, is expected to provide more rationale and economic results than the standard uncoupled approach, which, neglecting any energy dissipation at the foundation level, generally overestimates the ductility demand on the superstructure. Copyright © 2016 John Wiley & Sons, Ltd.     

Retrieval of local surface wave velocities from traffic noise – an example from the La Barge basin (Wyoming)

In regions where active source seismic exploration is constrained by limitations of energy penetration and recovery, cost and logistical concerns, or regulatory restrictions, analysis of natural source seismic data may provide an alternative. In this study, we investigate the feasibility of using locally‐generated seismic noise in the 2–6 Hz band to obtain a subsurface model via interferometric analysis. We apply this technique to three‐component data recorded during the La Barge Passive Seismic Experiment, a local deployment in south‐western Wyoming that recorded continuous seismic data between November 2008 and June 2009. We find traffic noise from a nearby state road to be the dominant source of surface waves recorded on the array and observe surface wave arrivals associated with this source up to distances of 5 kms. The orientation of the road with respect to the deployment ensures a large number of stationary points, leading to clear observations on both in‐line and cross‐line virtual source‐receiver pairs. This results in a large number of usable interferograms, which in turn enables the application of standard active source processing methods like signal processing, common offset stacking and traveltime inversion. We investigate the dependency of the interferograms on the amount of data, on a range of processing parameters and on the choice of the interferometry algorithm. The obtained interferograms exhibit a high signal‐to‐noise ratio on all three components. Rotation of the horizontal components to the radial/transverse direction facilitates the separation of Rayleigh and Love waves. Though the narrow frequency spectrum of the surface waves prevents the inversion for depth‐dependent shear‐wave velocities, we are able to map the arrival times of the surface waves to laterally varying group and phase velocities for both Rayleigh and Love waves. Our results correlate well with the known geological structure. We outline a scheme for obtaining localized surface wave velocities from local noise sources and show how the processing of passive data benefits from a combination with well‐established exploration seismology methods. We highlight the differences with interferometry applied to crustal scale data and conclude with recommendations for similar deployments.     

Simulations of non‐stationary frequency content and its importance to seismic assessment of structures

To verify the importance of the non‐stationary frequency characteristic of seismic ground motion, a joint time–frequency analysis technique of time signals, called chirplet‐based signal approximation, is developed to extract the non‐stationary frequency information from the recorded data. The chirplet‐based signal approximation is clear in concept, similar to Fourier Transform in mathematical expressions but with different base functions. Case studies show that the chirplet‐based signal approximation can represent the joint time–frequency variation of seismic ground motion quite well. Both the random models of uniform modulating process and evolutionary process are employed to generate artificial seismic waves. The joint time–frequency modulating function in the random model of evolutionary process is determined by chirplet‐based signal approximation. Finally, non‐linear response analysis of a SODF system and a frame structure is performed based on the generated artificial seismic waves. The results show that the non‐stationary frequency characteristic of seismic ground motion can significantly change the non‐linear response characteristics of structures, particularly when a structure goes into collapse phase under seismic action. It is concluded that non‐stationary frequency characteristic of seismic ground motion should be considered for the assessment of seismic capacity of structures. Copyright © 2002 John Wiley & Sons, Ltd.     

Retrieving virtual reflection responses at drill‐bit positions using seismic interferometry with drill‐bit noise

In the field of seismic interferometry, researchers have retrieved surface waves and body waves by cross‐correlating recordings of uncorrelated noise sources to extract useful subsurface information. The retrieved wavefields in most applications are between receivers. When the positions of the noise sources are known, inter‐source interferometry can be applied to retrieve the wavefields between sources, thus turning sources into virtual receivers. Previous applications of this form of interferometry assume impulsive point sources or transient sources with similar signatures. We investigate the requirements of applying inter‐source seismic interferometry using non‐transient noise sources with known positions to retrieve reflection responses at those positions and show the results using synthetic drilling noise as source. We show that, if pilot signals (estimates of the drill‐bit signals) are not available, it is required that the drill‐bit signals are the same and that the phases of the virtual reflections at drill‐bit positions can be retrieved by deconvolution interferometry or by cross‐coherence interferometry. Further, for this case, classic interferometry by cross‐correlation can be used if the source power spectrum can be estimated. If pilot signals are available, virtual reflection responses can be obtained by first using standard seismic‐while‐drilling processing techniques such as pilot cross‐correlation and pilot deconvolution to remove the drill‐bit signatures in the data and then applying cross‐correlation interferometry. Therefore, provided that pilot signals are reliable, drill‐bit data can be redatumed from surface to borehole depths using this inter‐source interferometry approach without any velocity information of the medium, and we show that a well‐positioned image below the borehole can be obtained using interferometrically redatumed reflection responses with just a simple velocity model. We discuss some of the practical hurdles that restrict the application of the proposed method offshore.     

The effects of reservoir geometry on the seismic response of gravity dams

Conventional seismic analysis of gravity dams assumes that the behaviour of the dam–water–soil system can be represented using a 2‐D model since dam vertical contraction joints between blocks allow them to vibrate independently from each other. The 2‐D model assumes the reservoir to be infinite and of constant width, which is not true for certain types of reservoirs. In this paper, a boundary element method (BEM) model in the frequency domain is used to investigate the influence of the reservoir geometry on the hydrodynamic dam response. Important conceptual conclusions about the dam–reservoir system behaviour are obtained using this model. The results show that the reservoir shape influences the seismic response of the dam, making it necessary to account for 3‐D effects in order to obtain accurate results. In particular, the 3‐D pressure and displacement responses can be substantially larger than those computed with the 2‐D model. Copyright © 2007 John Wiley & Sons, Ltd.     

Seeing coal‐seam top ahead of the drill bit through seismic‐while‐drilling

Blast damage to the tops of coal seams due to incorrect blast standoff distances is a serious issue, costing the industry in Australia about one open‐cut mine for every ten operating mines. The current approach for mapping coal‐seam tops is through drilling and pierce‐point logging. To provide appropriate depth control with accuracy of ±0.2 m for blast hole drilling, it is typically necessary to drill deep reconnaissance boreholes on a 50 m x 50 m grid well in advance of overburden removal. Pierce‐point mapping is expensive and can be inaccurate, particularly when the seam is disturbed by rolls, faults, and other obstacles.Numerical modelling and prototype‐field testing are used in this paper to demonstrate the feasibility of two seismic‐while‐drilling‐based approaches for predicting the approach to the top of coal during blast hole drilling: (i) reverse “walk‐away” vertical seismic profiling recording, in which the drill bit vibration provides the source signal and the geophones are planted on the surface near the drill rig, and (ii) in‐seam seismic recording, in which channel waves, driven by the coupling to the coal of the seismic signal emitted by the approaching drill bit, are guided by the seam to geophones located within the seam in nearby or remote boreholes.