Geostatistical seismic inversion for non‐stationary patterns using direct sequential simulation and co‐simulation

Geostatistical seismic inversion methods are routinely used in reservoir characterisation studies because of their potential to infer the spatial distribution of the petro‐elastic properties of interest (e.g., density, elastic, and acoustic impedance) along with the associated spatial uncertainty. Within the geostatistical seismic inversion framework, the retrieved inverse elastic models are conditioned by a global probability distribution function and a global spatial continuity model as estimated from the available well‐log data for the entire inversion grid. However, the spatial distribution of the real subsurface elastic properties is complex, heterogeneous, and, in many cases, non‐stationary since they directly depend on the subsurface geology, i.e., the spatial distribution of the facies of interest. In these complex geological settings, the application of a single distribution function and a spatial continuity model is not enough to properly model the natural variability of the elastic properties of interest. In this study, we propose a three‐dimensional geostatistical inversion technique that is able to incorporate the reservoir's heterogeneities. This method uses a traditional geostatistical seismic inversion conditioned by local multi‐distribution functions and spatial continuity models under non‐stationary conditions. The procedure of the proposed methodology is based on a zonation criterion along the vertical direction of the reservoir grid. Each zone can be defined by conventional seismic interpretation, with the identification of the main seismic units and significant variations of seismic amplitudes. The proposed method was applied to a highly non‐stationary synthetic seismic dataset with different levels of noise. The results of this work clearly show the advantages of the proposed method against conventional geostatistical seismic inversion procedures. It is important to highlight the impact of this technique in terms of higher convergence between real and inverted reflection seismic data and the more realistic approximation towards the real subsurface geology comparing with traditional techniques.     

Non‐stationary seismic response of tanks with soil interaction by wavelets

A wavelet‐based random vibration theory has been developed for the non‐stationary seismic response of liquid storage tanks including soil interaction. The ground motion process has been characterized via estimates of statistical functionals of wavelet coefficients obtained from a single time history of ground accelerations. The tank–liquid–soil system has been modelled as a two‐degree‐of‐freedom (2‐DOF) system. The wavelet domain equations have been formulated and the wavelet coefficients of the required response state are obtained by solving two linear simultaneous algebraic equations. The explicit expression for the instantaneous power spectral density function (PSDF) in terms of the functionals of the input wavelet coefficients has been obtained. The moments of this PSDF are used to estimate the expected pseudo‐spectral acceleration (PSA) response of the tank. Parametric variations are carried out to study the effects of tank height, foundation natural frequency, shear wave velocity of soil and ratio of the mass of tank (including liquid) to the mass of foundation on the PSA responses of tanks. Copyright © 2001 John Wiley & Sons, Ltd.     

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.     

Wavelet‐based non‐stationary seismic rocking response of flexibly supported tanks

The non‐stationary rocking response of liquid storage tanks under seismic base excitations including soil interaction has been developed based on the wavelet domain random vibration theory. The ground motion has been characterized through statistical functionals of wavelet coefficients of the ground acceleration history. The tank–liquid–foundation system is modelled as a multi‐degree‐of‐freedom (MDOF) system with both lateral and rocking motions of vibration of the foundation. The impulsive and convective modes of vibration of the liquid in the tank have been considered. The wavelet domain coupled dynamic equations are formulated and then solved to get the expressions of instantaneous power spectral density function (PSDF) in terms of functionals of input wavelet coefficients. The moments of the instantaneous PSDF are used to obtain the stochastic responses of the tank in the form of coefficients of hydrodynamic pressure, base shear and overturning base moment for the largest expected peak responses. Parametric variations are carried out to study the effects of various governing parameters like height of liquid in the tank, height–radius ratio of the tank, ratio of total liquid mass to mass of foundation, and shear wave velocity in the soil medium, on the responses of the tank. Copyright © 2003 John Wiley & Sons, Ltd.     

Wavelet‐based non‐stationary response analysis of a friction base‐isolated structure

A wavelet‐based stochastic formulation has been presented in this paper for the seismic analysis of a base‐isolated structural system which is modelled as a two‐degree‐of‐freedom (2‐DOF) system. The ground motion has been modelled as a non‐stationary process (both in amplitude and frequency) by using modified Littlewood–Paley basis wavelets. The proposed formulation is based on replacing the non‐linear system by an equivalent linear system with time‐dependent damping properties. The expressions of the instantaneous damping and the power spectral density function (PSDF) of the superstructure response have been obtained in terms of the functionals of input wavelet coefficients. The proposed formulation has been validated by simulating a ground motion process. The effect of the frequency non‐stationarity on the non‐linear response has also been studied in detail, and it has been clearly shown how ignoring the frequency non‐stationarity in the ground motion leads to inaccurate non‐linear response calculations. Copyright © 2000 John Wiley & Sons, Ltd.     

The use of low frequencies in a full‐waveform inversion and impedance inversion land seismic case study

Velocity model building and impedance inversion generally suffer from a lack of intermediate wavenumber content in seismic data. Intermediate wavenumbers may be retrieved directly from seismic data sets if enough low frequencies are recorded. Over the past years, improvements in acquisition have allowed us to obtain seismic data with a broader frequency spectrum. To illustrate the benefits of broadband acquisition, notably the recording of low frequencies, we discuss the inversion of land seismic data acquired in Inner Mongolia, China. This data set contains frequencies from 1.5–80 Hz. We show that the velocity estimate based on an acoustic full‐waveform inversion approach is superior to one obtained from reflection traveltime inversion because after full‐waveform inversion the background velocity conforms to geology. We also illustrate the added value of low frequencies in an impedance estimate.     

The non‐stationary cross‐correlation coefficients in the seismic analysis of MDOF structural systems

A method concerning the evaluation, in a very compact form, of the non‐stationary modal cross‐correlation coefficients of MDOF structural systems subjected to seismic excitations is presented. It is available both in the case when the excitation is considered as a white‐noise process and when it is considered as a filtered process. The evaluation of these coefficients is required when a transient seismic analysis is performed by the use of the modal response spectrum approach. This is necessary when the strong‐motion phase of the earthquake is significantly short with respect to the fundamental period of the structure. Copyright © 2003 John Wiley & Sons, Ltd.     

Seismic response of a cable‐stayed bridge deck under multi‐component non‐stationary random ground motion

A Markov method of analysis is presented for obtaining the seismic response of cable‐stayed bridges to non‐stationary random ground motion. A uniformly modulated non‐stationary model of the random ground motion is assumed which is specified by the evolutionary r.m.s. ground acceleration. Both vertical and horizontal components of the motion are considered to act simultaneously at the bridge supports. The analysis duly takes into account the angle of incidence of the earthquake, the spatial correlation of ground motion and the quasi‐static excitation. A cable‐stayed bridge is analysed under a set of parametric variations in order to study the non‐stationary response of the bridge. The results of the numerical study indicate that (i) frequency domain spectral analysis with peak r.m.s. acceleration as input could provide more r.m.s. response than the peak r.m.s. response obtained by the non‐stationary analysis; (ii) the longitudinal component of the ground motion significantly influences the vertical vibration of the bridge; and (iii) the angle of incidence of the earthquake has considerable influence on the deck response. Copyright © 2003 John Wiley & Sons, Ltd.     

Data‐driven time–frequency analysis of seismic data using non‐stationary Prony method

Empirical mode decomposition aims to decompose the input signal into a small number of components named intrinsic mode functions with slowly varying amplitudes and frequencies. In spite of its simplicity and usefulness, however, empirical mode decomposition lacks solid mathematical foundation. In this paper, we describe a method to extract the intrinsic mode functions of the input signal using non‐stationary Prony method. The proposed method captures the philosophy of the empirical mode decomposition but uses a different method to compute the intrinsic mode functions. Having the intrinsic mode functions obtained, we then compute the spectrum of the input signal using Hilbert transform. Synthetic and field data validate that the proposed method can correctly compute the spectrum of the input signal and could be used in seismic data analysis to facilitate interpretation.     

Alternative non‐linear demand estimation methods for probability‐based seismic assessments

Alternative non‐linear dynamic analysis procedures, using real ground motion records, can be used to make probability‐based seismic assessments. These procedures can be used both to obtain parameter estimates for specific probabilistic assessment criteria such as demand and capacity factored design and also to make direct probabilistic performance assessments using numerical methods. Multiple‐stripe analysis is a non‐linear dynamic analysis method that can be used for performance‐based assessments for a wide range of ground motion intensities and multiple performance objectives from onset of damage through global collapse. Alternatively, the amount of analysis effort needed in the performance assessments can be reduced by performing the structural analyses and estimating the main parameters in the region of ground motion intensity levels of interest. In particular, single‐stripe and double‐stripe analysis can provide local probabilistic demand assessments using minimal number of structural analyses (around 20 to 40). As a case study, the displacement‐based seismic performance of an older reinforced concrete frame structure, which is known to have suffered shear failure in its columns during the 1994 Northridge Earthquake, is evaluated. Copyright © 2008 John Wiley & Sons, Ltd.     

A method of calculating the non‐linear seismic response of a building braced with viscoelastic dampers

Buildings are continually subject to dynamic loads, such as wind load, seismic ground motion, and even the load from internal utility machines. The recent trend of constructing more flexible high‐rise buildings underscores the importance of including viscoelastic dampers in building designs. Viscoelastic dampers are used to control the dynamic response of a building. If the seismic design is based only on the linear response spectrum, considerable error may occur when calculating the seismic response of a building; rubber viscoelastic dampers show non‐linear hysteretic damping that is quite different from viscous damping. This study generated a non‐linear response spectrum using a non‐linear oscillator model to simulate a building with viscoelastic dampers installed. The parameters used in the non‐linear damper model were obtained experimentally from dynamic loading tests. The results show that viscoelastic dampers effectively reduce the seismic displacement response of a structure, but transmit more seismic force to the structure, which essentially increases its seismic acceleration response. Copyright © 2003 John Wiley & Sons, Ltd.     

Modelling exterior unreinforced beam‐column joints in seismic analysis of non‐ductile RC frames

The seismic response of non‐ductile reinforced concrete (RC) buildings can be affected by the behaviour of beam‐column joints involved in the failure mechanism, especially in typical existing buildings. Conventional modelling approaches consider only beam and column flexibility, although joints can provide a significant contribution also to the overall frame deformability. In this study, the attention is focused on exterior joints without transverse reinforcement, and a possible approach to their modelling in nonlinear seismic analysis of RC frames is proposed. First, experimental tests performed by the authors are briefly presented, and their results are discussed. Second, these tests, together with other tests with similar features from literature, are employed to calibrate the joint panel deformability contribution in order to reproduce numerically the experimental joint shear stress–strain behaviour under cyclic loading. After a validation phase of this proposal, a numerical investigation of the influence of joints on the seismic behaviour of a case study RC frame – designed for gravity loads only – is performed. The preliminary failure mode classification of the joints within the analysed frame is carried out. Structural models that (i) explicitly include nonlinear behaviour of beam‐column joints exhibiting shear or anchorage failure or (ii) model joints as elements with infinite strength and stiffness are built and their seismic performance are assessed and compared. A probabilistic assessment based on nonlinear dynamic simulations is performed by means of a scaling approach to evaluate the seismic response at different damage states accounting for uncertainties in ground‐motion records. Copyright © 2016 John Wiley & Sons, Ltd.     

Dynamic analysis of train–bridge system subjected to non‐uniform seismic excitations

Based on the theory of dynamic wheel–rail interactions, a dynamic model of coupled train–bridge system subjected to earthquakes is established, in which the non‐uniform characteristics of the seismic wave input from different foundations are considered. The bridge model is based on the modal comprehension analysis technique. Each vehicle is modelled with 31 degrees of freedom. The seismic loads are imposed on the bridge by using the influence matrix and exerted on the vehicles through the dynamic wheel–rail interaction relationships. The normal wheel–rail interaction is tackled by using the Hertzian contact theory, and the tangent wheel–rail interaction by the Kalker linear theory and the Shen–Hedrick–Elkins theory. A computer code is developed. A case study is performed to a continuous bridge on the planned Beijing–Shanghai high‐speed railway in China. Through input of typical seismic waves with different propagation velocities to the train–bridge system, the histories of the train running through the bridge are simulated and the dynamic responses of the bridge and the vehicles are calculated. The influences of train speed and seismic wave propagation velocity on the dynamic responses of the bridge–vehicle system are studied. The critical train speeds are proposed for running safety on high‐speed railway bridges under earthquakes of various intensities. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

Earthquake‐induced interaction between adjacent reinforced concrete structures with non‐equal heights

The influence of the structural pounding on the ductility requirements and the seismic behaviour of reinforced concrete structures designed to EC2 and EC8 with non‐equal heights is investigated. Special purpose elements of distributed plasticity are employed for the study of the columns. Two distinct types of the problem are identified: Type A, where collisions may occur only between storey masses; and Type B, where the slabs of the first structure hit the columns of the other (72 Type A and 36 Type B pounding cases are examined). Type A cases yielded critical ductility requirements for the columns in the pounding area mainly for the cases where the structures were in contact from the beginning of the excitation. In both pounding types the ductility requirements of the columns of the taller building are substantially increased for the floors above the highest contact storey level probably due to a whiplash behaviour. The most important issue in the pounding type B is the local response of the column of the tall structure that suffers the hit of the upper floor slab of the adjacent shorter structure. In all the examined cases this column was in a critical condition due to shear action and in the cases where the structures were in contact from the beginning of the excitation, this column was also critical due to high ductility demands. It can be summarized that in situations of potential pounding, neglecting its possible effects leads to non‐conservative building design or evaluation that may become critical in some cases. Copyright © 2004 John Wiley & Sons, Ltd.     

A new method for non‐stationary random critical excitation

Since earthquake ground motions are very uncertain even with the present knowledge, it is desirable to develop a robust structural design method taking into account these uncertainties. Critical excitation approaches are promising and a new non‐stationary random critical excitation method is proposed. In contrast to the conventional critical excitation methods, a stochastic response index is treated as the objective function to be maximized. The power (area of power spectral density (PSD) function) and the intensity (magnitude of PSD function) are fixed and the critical excitation is found under these restrictions. It is shown that the original idea for stationary random inputs can be utilized effectively in the procedure for finding a critical excitation for non‐stationary random inputs. The key for finding the new non‐stationary random critical excitation is the exchange of the order of the double maximization procedures with respect to time and to the power spectral density function. Copyright © 2001 John Wiley & Sons, Ltd.     

A macro‐element model for non‐linear soil–shallow foundation–structure interaction under seismic loads: theoretical development and experimental validation on large scale tests

In this paper, different formulations of a macro‐element model for non‐linear dynamic soil‐structure interaction analyses of structures lying on shallow foundations are first reviewed, and secondly, a novel formulation is introduced, which combines some of the characteristics of previous approaches with several additional features. This macro‐element allows one to model soil‐footing geometric (uplift) and material (soil plasticity) non‐linearities that are coupled through a stiffness degradation model. Footing uplift is introduced by a simple non‐linear elastic model based on the concept of effective foundation width, whereas soil plasticity is treated by means of a bounding surface approach in which a vertical load mapping rule is implemented. This mapping is particularly suited for the seismic loading case for which the proposed model has been conceived. The new macro‐element is subsequently validated using cyclic and dynamic large‐scale laboratory tests of shallow foundations on dense sand, namely: the TRISEE cyclic tests, the Public Works Research Institute and CAMUS IV shaking table tests. Based on this comprehensive validation process against a set of independent experimental results, a unique set of macro‐element parameters for shallow foundations on dense sand is proposed, which can be used to perform predictive analyses by means of the present model. Copyright © 2011 John Wiley & Sons, Ltd.     

A new passive rolling‐pendulum vibration absorber using a non‐axial‐symmetrical guide to achieve bidirectional tuning

A new, passive, vibroprotective device of the rolling‐pendulum tuned mass damper type is presented that, relying on a proper three‐dimensional guiding surface, can simultaneously control the response of the supporting structure in two mutually orthogonal horizontal directions. Unlike existing examples of ball vibration absorbers, mounted on spherical recesses and effective for axial‐symmetrical structures, the new device is bidirectionally tuneable, by virtue of the optimum shape of the rolling cavity, to both fundamental structural modes, even when the corresponding natural frequencies are different, in such a case recurring to an innovative non‐axial‐symmetrical rolling guide. Based on Appell's non‐holonomic mechanics, a non‐linear dynamic model is first derived for the bidirectional absorber mounted on a 1‐storey 3‐degrees‐of‐freedom linear structure translating under the effect of both imparted base motion and applied dynamic forces. A laboratory‐scaled prototype of the device is then tested to experimentally demonstrate the bidirectional tuning capability and to validate the mathematical model. The design procedure and the seismic performance of the absorber are finally exemplified through numerical simulation. Copyright © 2009 John Wiley & Sons, Ltd.     

Modeling of a large‐scale magneto‐rheological damper for seismic hazard mitigation. Part I: Passive mode

Magneto‐rheological (MR) dampers are a promising device for seismic hazard mitigation because their damping characteristics can be varied adaptively using an appropriate control law. During the last few decades researchers have investigated the behavior of MR dampers and semi‐active control laws associated with these types of dampers for earthquake hazard mitigation. A majority of this research has involved small‐scale MR dampers. To investigate the dynamic behavior of a large‐scale MR damper, characterization tests were conducted at the Lehigh Network for Earthquake Engineering Simulation equipment site on large‐scale MR dampers. A new MR damper model, called the Maxwell Nonlinear Slider (MNS) model, is developed based on the characterization tests and is reported in this paper. The MNS model can independently describe the pre‐yield and post‐yield behavior of an MR damper, which makes it easy to identify the model parameters. The MNS model utilizes Hershel–Bulkley visco‐plasticity to describe the post‐yield non‐Newtonian fluid behavior, that is, shear thinning and thickening behavior, of the MR fluid that occurs in the dampers. The predicted response of a large‐scale damper from the MNS model along with that from existing Bouc–Wen and hyperbolic tangent models, are compared with measured response from various experiments. The comparisons show that the MNS model achieves better accuracy than the existing models in predicting damper response under cyclic loading. Copyright © 2012 John Wiley & Sons, Ltd.     

Preconditioned non‐linear conjugate gradient method for frequency domain full‐waveform seismic inversion

We present preconditioned non‐linear conjugate gradient algorithms as alternatives to the Gauss‐Newton method for frequency domain full‐waveform seismic inversion. We designed two preconditioning operators. For the first preconditioner, we introduce the inverse of an approximate sparse Hessian matrix. The approximate Hessian matrix, which is highly sparse, is constructed by judiciously truncating the Gauss‐Newton Hessian matrix based on examining the auto‐correlation and cross‐correlation of the Jacobian matrix. As the second preconditioner, we employ the approximation of the inverse of the Gauss‐Newton Hessian matrix. This preconditioner is constructed by terminating the iteration process of the conjugate gradient least‐squares method, which is used for inverting the Hessian matrix before it converges. In our preconditioned non‐linear conjugate gradient algorithms, the step‐length along the search direction, which is a crucial factor for the convergence, is carefully chosen to maximize the reduction of the cost function after each iteration. The numerical simulation results show that by including a very limited number of non‐zero elements in the approximate Hessian, the first preconditioned non‐linear conjugate gradient algorithm is able to yield comparable inversion results to the Gauss‐Newton method while maintaining the efficiency of the un‐preconditioned non‐linear conjugate gradient method. The only extra cost is the computation of the inverse of the approximate sparse Hessian matrix, which is less expensive than the computation of a forward simulation of one source at one frequency of operation. The second preconditioned non‐linear conjugate gradient algorithm also significantly saves the computational expense in comparison with the Gauss‐Newton method while maintaining the Gauss‐Newton reconstruction quality. However, this second preconditioned non‐linear conjugate gradient algorithm is more expensive than the first one.