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.     

Pressure effects on the joint elastic‐electrical properties of reservoir sandstones

We conducted a laboratory study of the joint elastic‐electrical properties of sixty‐three brine‐saturated sandstone samples to assess the likely impact of differential pressure (confining minus pore fluid pressures) in the range 8–60 MPa on the joint interpretation of marine seismic and controlled‐source electromagnetic survey data. The samples showed a wide range of petrophysical properties representative of most sandstone reservoirs. We found that a regression equation comprising both a constant and an exponential part gave a good fit to the pressure dependence of all five measured geophysical parameters (ultrasonic P‐ and S‐wave velocity, attenuation and electrical resistivity). Electrical resistivity was more pressure‐sensitive in clay‐rich sandstones with higher concentrations of low aspect ratio pores and micropores than in clean sandstones. Attenuation was more pressure‐sensitive in clean sandstones with large open pores (macropores) than in clay‐rich sandstones. Pore shape did not show any influence on the pressure sensitivity of elastic velocity. As differential pressure increases, the effect of the low aspect ratio pores and micropores on electrical resistivity becomes stronger than the effect of the macropores on attenuation. Further analysis of correlations among the five parameters as a function of pressure revealed potentially diagnostic relationships for geopressure prediction in reservoir sandstones.     

Joint PP‐ and PSV‐wave amplitudes versus offset and azimuth inversion for fracture compliances in horizontal transversely isotropic media

We propose a robust approach for the joint inversion of PP‐ and PSV‐wave angle gathers along different azimuths for the elastic properties of the homogeneous isotropic host rock and excess compliances due to the presence of fractures. Motivated by the expression of fluid content indicator in fractured reservoirs and the sensitivity of Lamé impedances to fluid type, we derive PP‐ and PSV‐wave reflection coefficients in terms of Lamé impedances, density, and fracture compliances for an interface separating two horizontal transversely isotropic media. Following a Bayesian framework, we construct an objective function that includes initial models. We employ the iteratively reweighted least‐squares algorithm to solve the inversion problem to estimate unknown parameters (i.e., Lamé impedances, density, and fracture compliances) from PP‐ and PSV‐wave angle gathers along different azimuths. Synthetic tests reveal that the unknown parameters estimated using the joint inversion approach match true values better than those estimated using a PP‐wave amplitude inversion only. A real data test indicates that reasonable results for subsurface fracture detection are obtained from the joint inversion approach.     

Joint inversion of Rayleigh‐wave dispersion data and vertical electric sounding data: synthetic tests on characteristic sub‐surface models

In the traditional inversion of the Rayleigh dispersion curve, layer thickness, which is the second most sensitive parameter of modelling the Rayleigh dispersion curve, is usually assumed as correct and is used as fixed a priori information. Because the knowledge of the layer thickness is typically not precise, the use of such a priori information may result in the traditional Rayleigh dispersion curve inversions getting trapped in some local minima and may show results that are far from the real solution. In this study, we try to avoid this issue by using a joint inversion of the Rayleigh dispersion curve data with vertical electric sounding data, where we use the common‐layer thickness to couple the two methods. The key idea of the proposed joint inversion scheme is to combine methods in one joint Jacobian matrix and to invert for layer S‐wave velocity, resistivity, and layer thickness as an additional parameter, in contrast with a traditional Rayleigh dispersion curve inversion. The proposed joint inversion approach is tested with noise‐free and Gaussian noise data on six characteristic, synthetic sub‐surface models: a model with a typical dispersion; a low‐velocity, half‐space model; a model with particularly stiff and soft layers, respectively; and a model reproduced from the stiff and soft layers for different layer‐resistivity propagation. In the joint inversion process, the non‐linear damped least squares method is used together with the singular value decomposition approach to find a proper damping value for each iteration. The proposed joint inversion scheme tests many damping values, and it chooses the one that best approximates the observed data in the current iteration. The quality of the joint inversion is checked with the relative distance measure. In addition, a sensitivity analysis is performed for the typical dispersive sub‐surface model to illustrate the benefits of the proposed joint scheme. The results of synthetic models revealed that the combination of the Rayleigh dispersion curve and vertical electric sounding methods in a joint scheme allows to provide reliable sub‐surface models even in complex and challenging situations and without using any a priori information.     

Numerical models and ductile ultimate deformation response of post‐tensioned self‐centering moment connections

New steel moment‐resisting connections that incorporate post‐tensioning elements to provide a self‐centering capacity and devices to dissipate seismic input energy have recently been proposed and experimentally validated. Experimental studies have confirmed that these connections are capable of undergoing large lateral deformations with negligible residual drifts. To facilitate their implementation, accurate modeling of the behavior of systems incorporating post‐tensioned connections must be readily available to designers and researchers. A number of simplified models have been suggested in the literature by researchers trying to capture experimental results at the beam–column connections and thereby to predict the global response of structures incorporating such connections. To provide a clear set of guidelines for the modeling of post‐tensioned steel frames, for practicing engineers as well as researchers, in this paper three types of numerical models of increasing complexity are presented: (i) a sectional analysis procedure, (ii) a lumped plasticity spring frame leveled approach and (iii) a non‐linear solid finite element analysis to predict the response at ultimate deformation levels. The analytical results obtained from the numerical models predict well the structural behavior of these connections when compared with available experimental data. Even at the ultimate deformation level, analytical results are in good agreement with test results. Furthermore, detailing requirements are proposed to assure that flexural hinges form in the beams in order to improve the cyclic response of steel self‐centering connections when drifts exceeding the design drifts are imposed to the system. Experimental and analytical studies demonstrate that steel post‐tensioned self‐centering connections incorporating the proposed detailing in the beams develop an increased deformation capacity and thereby exhibit a ductile response while avoiding a sudden loss of their strength and stiffness. Copyright © 2008 John Wiley & Sons, Ltd.     

On the winter evolution of snow thermophysical properties over land‐fast first‐year sea ice

The geophysical, thermodynamic and dielectric properties of snow are important state variables that are known to be sensitive to Arctic climate variability and change. Given recent observations of changes in the Arctic physical system (Arctic Climate Impact Assessment, 2004), it is important to focus on the processes that give rise to variability in the horizontal, vertical and temporal dimensions of the life‐history of snow on sea ice. The objectives in this study are to present these ‘state’ variables and to investigate the processes that govern variability in the vertical, horizontal and temporal dimension by using a case study over land‐fast first‐year sea ice for the period December 2003 to June 2004. Results from two sampling areas (thin and thick snowpacks) show that differences in snowpack thickness can substantially change the vertical and temporal evolution of snow properties. During the late fall and early winter (cooling period) we measured no significant changes in the physical properties, except for thin snow‐cover salinity, which decreased throughout the period. Fall‐snow desalination was only observed under thin snowpacks with a rate of ?0·12 ppt day^?1. Significant changes occurred in the late winter and early spring (warming period), especially for snow grain size. Snow grain kinetic growth of 0·25–0·48 mm·day^?1 was measured coincidently with increasing salinity and wetness for both thin and thick snowpacks. Copyright © 2006 John Wiley & Sons, Ltd.     

Assessment of rainfall‐runoff models based upon wavelet analysis

A basic hypothesis is proposed: given that wavelet‐based analysis has been used to interpret runoff time‐series, it may be extended to evaluation of rainfall‐runoff model results. Conventional objective functions make certain assumptions about the data series to which they are applied (e.g. uncorrelated error, homoscedasticity). The difficulty that objective functions have in distinguishing between different realizations of the same model, or different models of the same system, is that they may have contributed in part to the occurrence of model equifinality. Of particular concern is the fact that the error present in a rainfall‐runoff model may be time dependent, requiring some form of time localization in both identification of error and derivation of global objective functions. We explore the use of a complex Gaussian (order 2) wavelet to describe: (1) a measured hydrograph; (2) the same hydrograph with different simulated errors introduced; and (3) model predictions of the same hydrograph based upon a modified form of TOPMODEL. The analysis of results was based upon: (a) differences in wavelet power (the wavelet power error) between the measured hydrograph and both the simulated error and modelled hydrographs; and (b) the wavelet phase. Power difference and wavelet phase were used to develop two objective functions, RMSE(power) and RMS(phase), which were shown to distinguish between simulated errors and model predictions with similar values of the commonly adopted Nash‐Sutcliffe efficiency index. These objective functions suffer because they do not retain time, frequency or time‐frequency localization. Consideration of wavelet power spectra and time‐ and frequency‐integrated power spectra shows that the impacts of different types of simulated error can be seen through retention of some localization, especially in relation to when and the scale over which error was manifest. Theoretical objections to the use of wavelet analysis for this type of application are noted, especially in relation to the dependence of findings upon the wavelet chosen. However, it is argued that the benefits of localization and the qualitatively low sensitivity of wavelet power and phase to wavelet choice are sufficient to warrant further exploration of wavelet‐based approaches to rainfall‐runoff model evaluation. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical simulation of reservoir sediment and effects on hydro‐dynamic response of arch dams

Based on the dynamic theory for saturated porous media by Biot (Journal of the Acoustical Society of America 1956; 28 : 168–178), a numerical model is presented to analyse the reflection behaviours of reservoir sediment and compared with those from the visco‐elastic model. It is concluded that the two models give very similar results of reflection coefficient α within the frequency range of interest. Then, using the two models, the change of the reflection coefficients α with various sedimentation parameters and excitation frequencies are studied in detail. The results are further used in the analysis of response functions of hydro‐dynamic pressures on, and structural displacements of the Xiang Hong Dian arch dam, for which some results from a field vibration test are available. It appears that effects of water compressibility with sediment reflection on hydro‐dynamic pressures and structural response are not significant for this specific case. Copyright © 2001 John Wiley & Sons, Ltd.     

Assessing the effect of vegetation‐related bank strength on channel morphology and stability in gravel‐bed streams using numerical models

Bank strength due to vegetation dominates the geometry of small stream channels, but has virtually no effect on the geometry of larger ones. The dependence of bank strength on channel scale affects the form of downstream hydraulic geometry relations and the meandering‐braiding threshold. It is also associated with a lateral migration threshold discharge, below which channels do not migrate appreciably across their floodplains. A rational regime model is used to explore these scale effects: it parameterizes vegetation‐related bank strength using a dimensionless effective cohesion, C_r*. The scale effects are explored primarily using an alluvial state space defined by the dimensionless formative discharge, Q*, and channel slope, S, which is analogous to the Q–S diagrams originally used to explore meandering‐braiding thresholds. The analyses show that the effect of vegetation on both downstream hydraulic geometry and the meandering‐braiding threshold is strongest for the smallest streams in a watershed, but that the effect disappears for Q* > 10⁶. The analysis of the migration threshold suggests that the critical discharge ranges from about 5 m³/s to 50 m³/s, depending on the characteristic rooting depth for the vegetation. The analysis also suggests that, where fires frequently affect riparian forests, channels may alternate between laterally stable gravel plane‐bed channels and laterally active riffle‐pool channels. These channels likely do not exhibit the classic dynamic equilibrium associated with alluvial streams, but instead exhibit a cyclical morphologic evolution, oscillating between laterally stable and laterally unstable end‐members with a frequency determined by the forest fire recurrence interval. Copyright © 2008 John Wiley & Sons, Ltd.     

Near‐fault ground motions,and the response of elastic and inelastic single‐degree‐of‐freedom (SDOF) systems

In order to investigate the response of structures to near‐fault seismic excitations, the ground motion input should be properly characterized and parameterized in terms of simple, yet accurate and reliable, mathematical models whose input parameters have a clear physical interpretation and scale, to the extent possible, with earthquake magnitude. Such a mathematical model for the representation of the coherent (long‐period) ground motion components has been proposed by the authors in a previous study and is being exploited in this article for the investigation of the elastic and inelastic response of the single‐degree‐of‐freedom (SDOF) system to near‐fault seismic excitations. A parametric analysis of the dynamic response of the SDOF system as a function of the input parameters of the mathematical model is performed to gain insight regarding the near‐fault ground motion characteristics that significantly affect the elastic and inelastic structural performance. A parameter of the mathematical representation of near‐fault motions, referred to as ‘pulse duration’ (T_P), emerges as a key parameter of the problem under investigation. Specifically, T_P is employed to normalize the elastic and inelastic response spectra of actual near‐fault strong ground motion records. Such normalization makes feasible the specification of design spectra and reduction factors appropriate for near‐fault ground motions. The ‘pulse duration’ (T_P) is related to an important parameter of the rupture process referred to as ‘rise time’ (τ) which is controlled by the dimension of the sub‐events that compose the mainshock. Copyright © 2004 John Wiley & Sons, Ltd.     

Response evaluation of interconnected electrical substation equipment using real‐time hybrid simulation on multiple shaking tables

In an attempt to quantify the conductor cable effect on substation electrical equipment, real‐time hybrid simulation (RTHS) is conducted on interconnected equipment using two shaking tables. For this purpose, the existing RTHS system with advanced control capabilities at the Pacific Earthquake Engineering Research Center structural laboratory is enhanced to accommodate the simultaneous use of two shaking tables. An experimental parametric study is conducted to investigate the conductor cable effect using this system with a two‐table RTHS setup. Post insulators of disconnect switches, important components of substations that are usually tested with conventional methods for evaluating their seismic performance, are utilized as experimental substructures for realistic representation of the electrical equipment. Various global and local response parameters, including accelerations, forces, displacements, and strains, are considered to evaluate the effect of the tested conductor cable configuration for a wide range of support structure configurations, which are modeled in the computer as analytical substructures. The experimental parametric study results indicate that the conductor cable has a significant effect on the response of the interconnected equipment over the whole range of investigated support structures and needs to be explicitly considered for seismic testing of electrical equipment. Copyright © 2016 John Wiley & Sons, Ltd.     

Evaluation of predictors of non‐linear seismic demands using ‘fishbone’ models of SMRF buildings

Predictors (or estimates) of seismic structural demands that are less computationally time‐consuming than non‐linear dynamic analysis can be useful for structural performance assessment and for design. In this paper, we evaluate the bias and precision of predictors that make use of, at most, (i) elastic modal vibration properties of the given structure, (ii) the results of a non‐linear static pushover analysis of the structure, and (iii) elastic and inelastic single‐degree‐of‐freedom time‐history analyses for the specified ground motion record. The main predictor of interest is an extension of first‐mode elastic spectral acceleration that additionally takes into account both the second‐mode contribution to (elastic) structural response and the effects of inelasticity. This predictor is evaluated with respect to non‐linear dynamic analysis results for ‘fishbone’ models of steel moment‐resisting frame (SMRF) buildings. The relatively small number of degrees of freedom for each fishbone model allows us to consider several short‐to‐long period buildings and numerous near‐ and far‐field earthquake ground motions of interest in both Japan and the U.S. Before doing so, though, we verify that estimates of the bias and precision of the predictor obtained using fishbone models are effectively equivalent to those based on typical ‘full‐frame’ models of the same buildings. Copyright © 2003 John Wiley & Sons, Ltd.     

Multi‐depth monitoring of electrical conductivity and temperature of groundwater at a multilayered coastal aquifer: Jeju Island,Korea

To supplement conventional geophysical log data, this study presents temporal variations in electrical conductivity (EC) and temperature with depth in a multilayered coastal aquifer, on the eastern part of Jeju Island, Korea. One‐month time‐series data obtained at eight points from a multi‐depth monitoring system showed that semidiurnal and semimonthly tidal variations induced dynamic fluctuations in EC and temperature. At some depths, EC ranged from 1483 to 26 822 µS cm^?1, while some points showed no significant variations. The results of EC log and time‐series data revealed that a sharp fresh‐saltwater interface occurred at low tide, but the diffusion zone broadened to 20 m at high tide. EC, temperature, and tide level data were used for the cross‐correlation analysis. The response time of EC and temperature to tide appears to range from less than 30 min to 11 h. Using end‐member mixing analysis (EMMA), the fraction of variations of chloride concentration in the multilayered aquifer was explained, and a conceptual model was developed which subdivided the coastal aquifer into four vertical zones. The percentage of water derived from seawater varied from 2 to 48 at specific depth, owing to tidal fluctuations. Continuous observations of EC and temperature at multiple depths are powerful tools for quantifying the transport of saline water by tidal variations in multilayered coastal aquifers. Copyright © 2008 John Wiley & Sons, Ltd.     

Complex resistivity of mineral rocks in the context of the generalised effective‐medium theory of the induced polarisation effect

This paper develops the generalised effective‐medium theory of induced polarisation for rock models with elliptical grains and applies this theory to studying the complex resistivity of typical mineral rocks. We first demonstrate that the developed generalised effective‐medium theory of induced polarisation model can correctly represent the induced polarisation phenomenon in multiphase artificial rock samples manufactured using pyrite and magnetite particles. We have also collected representative rock samples from the Cu–Au deposit in Mongolia and subjected them to mineralogical analysis using Quantitative Evaluation of Minerals by Scanning Electron Microscopy technology. The electrical properties of the same samples were determined using laboratory complex resistivity measurements. As a result, we have established relationships between the mineral composition of the rocks, determined using Quantitative Evaluation of Minerals by Scanning Electron Microscopy analysis, and the parameters of the generalised effective‐medium theory of induced polarisation model defined from the laboratory measurements of the electrical properties of the rocks. These relationships open the possibility for remote estimation of types of mineralisation and for mineral discrimination using spectral induced polarization data.     

IM‐based and EDP‐based decision models for the verification of the seismic collapse safety of buildings

Decision models for the verification of seismic collapse safety of buildings are introduced. The derivations are based on the concept of the acceptable (target) annual probability of collapse, whereas the decision making involves comparisons between seismic demand and capacity, which is familiar to engineering practitioners. Seismic demand, which corresponds to the design seismic action associated with a selected return period, can be expressed either in terms of an intensity measure (IM) or an engineering demand parameter (EDP). Seismic capacity, on the other hand, is defined by dividing the near‐collapse limit‐state IM or EDP by an appropriate risk‐targeted safety factor (γ _im or γ _edp ), which is the only safety factor used in the proposed decision model. Consequently, the seismic performance assessment of a building should be based on the best possible estimate. For a case study, it is shown that if the target collapse risk is set to 10^?4 (0.5% over a period of 50 years), and if the seismic demand corresponds to a return period of 475 years (10% over a period of 50 years), then it can be demonstrated that γ _im is approximately equal to 2.5 for very stiff buildings, whereas for buildings with long periods the value of γ _im can increase up to a value of approximately 5. The model using γ _edp is equal to that using γ _im only if it can be assumed that displacements, with consideration of nonlinear behavior, are equal to displacements from linear elastic analysis.     

Bayesian inversion of time‐lapse seismic data for the estimation of static reservoir properties and dynamic property changes

Seismic conditioning of static reservoir model properties such as porosity and lithology has traditionally been faced as a solution of an inverse problem. Dynamic reservoir model properties have been constrained by time‐lapse seismic data. Here, we propose a methodology to jointly estimate rock properties (such as porosity) and dynamic property changes (such as pressure and saturation changes) from time‐lapse seismic data. The methodology is based on a full Bayesian approach to seismic inversion and can be divided into two steps. First we estimate the conditional probability of elastic properties and their relative changes; then we estimate the posterior probability of rock properties and dynamic property changes. We apply the proposed methodology to a synthetic reservoir study where we have created a synthetic seismic survey for a real dynamic reservoir model including pre‐production and production scenarios. The final result is a set of point‐wise probability distributions that allow us to predict the most probable reservoir models at each time step and to evaluate the associated uncertainty. Finally we also show an application to real field data from the Norwegian Sea, where we estimate changes in gas saturation and pressure from time‐lapse seismic amplitude differences. The inverted results show the hydrocarbon displacement at the times of two repeated seismic surveys.     

Critical damping of structures with elastically supported visco‐elastic dampers

This paper presents a new formulation for critical damping of structures with elastically supported visco‐elastic dampers.Owing to the great dependence of damper performance on the support stiffness, this model is inevitable for reliable modelling of structures with visco‐elastic dampers. It is shown that the governing equation of free vibration of this model is reduced to a third‐order differential equation and the conventional method for defining the critical damping for second‐order differential equations cannot be applied to the present model. It is demonstrated that the region of overdamped vibration is finite in contrast to that (semi‐infinite) for second‐order differential equations and multiple critical damping coefficients exist. However, it turns out that the smaller one is practically meaningful. Copyright © 2001 John Wiley & Sons, Ltd.     

Comparative analysis of relationships between NLDAS‐2 forcings and model outputs

This study reports results from an analysis of the relationship between atmospheric forcing and model‐simulated water and energy fluxes for the North American Land Data Assimilation System Project Phase 2 (NLDAS‐2). The relationships between mean monthly precipitation and total runoff are stronger in the Sacramento (SAC) and variable infiltration capacity (VIC) models, which grew out of the hydrological community, than in the Noah and Mosaic models, which grew out of the soil‐vegetation‐atmosphere transfer (SVAT) community. The reverse is true for the relationship between mean monthly precipitation and evapotranspiration. In addition, surface energy fluxes in VIC are less sensitive to model forcing (except for air temperature) than those in the Noah and Mosaic model. Notwithstanding these general conclusions, the relationships between forcings and model‐simulated water and energy fluxes for all models vary for different seasons, variables, and regions. These findings will ultimately inspire a combination of SVAT‐type model energy components with hydrological model water components to develop a SVAT‐hydrology model to improve both evapotranspiration and total runoff simulations. Copyright © 2011 John Wiley & Sons, Ltd.