A model‐based data‐driven dictionary learning for seismic data representation

Planar waves events recorded in a seismic array can be represented as lines in the Fourier domain. However, in the real world, seismic events usually have curvature or amplitude variability, which means that their Fourier transforms are no longer strictly linear but rather occupy conic regions of the Fourier domain that are narrow at low frequencies but broaden at high frequencies where the effect of curvature becomes more pronounced. One can consider these regions as localised “signal cones”. In this work, we consider a space–time variable signal cone to model the seismic data. The variability of the signal cone is obtained through scaling, slanting, and translation of the kernel for cone‐limited (C‐limited) functions (functions whose Fourier transform lives within a cone) or C‐Gaussian function (a multivariate function whose Fourier transform decays exponentially with respect to slowness and frequency), which constitutes our dictionary. We find a discrete number of scaling, slanting, and translation parameters from a continuum by optimally matching the data. This is a non‐linear optimisation problem, which we address by a fixed‐point method that utilises a variable projection method with ?₁ constraints on the linear parameters and bound constraints on the non‐linear parameters. We observe that slow decay and oscillatory behaviour of the kernel for C‐limited functions constitute bottlenecks for the optimisation problem, which we partially overcome by the C‐Gaussian function. We demonstrate our method through an interpolation example. We present the interpolation result using the estimated parameters obtained from the proposed method and compare it with those obtained using sparsity‐promoting curvelet decomposition, matching pursuit Fourier interpolation, and sparsity‐promoting plane‐wave decomposition methods.     

Generalized non‐hyperbolic approximation for qP‐wave relative geometrical spreading in a layered transversely isotropic medium with a vertical symmetry axis

Compensation for geometrical spreading along the ray‐path is important in amplitude variation with offset analysis especially for not strongly attenuative media since it contributes to the seismic amplitude preservation. The P‐wave geometrical spreading factor is described by a non‐hyperbolic moveout approximation using the traveltime parameters that can be estimated from the velocity analysis. We extend the P‐wave relative geometrical spreading approximation from the rational form to the generalized non‐hyperbolic form in a transversely isotropic medium with a vertical symmetry axis. The acoustic approximation is used to reduce the number of parameters. The proposed generalized non‐hyperbolic approximation is developed with parameters defined by two rays: vertical and a reference rays. For numerical examples, we consider two choices for parameter selection by using two specific orientations for reference ray. We observe from the numerical tests that the proposed generalized non‐hyperbolic approximation gives more accurate results in both homogeneous and multi‐layered models than the rational counterpart.     

Piecewise‐continuous velocity inversion

We suggest a new method to determine the piecewise‐continuous vertical distribution of instantaneous velocities within sediment layers, using different order time‐domain effective velocities on their top and bottom points. We demonstrate our method using a synthetic model that consists of different compacted sediment layers characterized by monotonously increasing velocity, combined with hard rock layers, such as salt or basalt, characterized by constant fast velocities, and low velocity layers, such as gas pockets. We first show that, by using only the root‐mean‐square velocities and the corresponding vertical travel times (computed from the original instantaneous velocity in depth) as input for a Dix‐type inversion, many different vertical distributions of the instantaneous velocities can be obtained (inverted). Some geological constraints, such as limiting the values of the inverted vertical velocity gradients, should be applied in order to obtain more geologically plausible velocity profiles. In order to limit the non‐uniqueness of the inverted velocities, additional information should be added. We have derived three different inversion solutions that yield the correct instantaneous velocity, avoiding any a priori geological constraints. The additional data at the interface points contain either the average velocities (or depths) or the fourth‐order average velocities, or both. Practically, average velocities can be obtained from nearby wells, whereas the fourth‐order average velocity can be estimated from the quartic moveout term during velocity analysis. Along with the three different types of input, we consider two types of vertical velocity models within each interval: distribution with a constant velocity gradient and an exponential asymptotically bounded velocity model, which is in particular important for modelling thick layers. It has been shown that, in the case of thin intervals, both models lead to similar results. The method allows us to establish the instantaneous velocities at the top and bottom interfaces, where the velocity profile inside the intervals is given by either the linear or the exponential asymptotically bounded velocity models. Since the velocity parameters of each interval are independently inverted, discontinuities of the instantaneous velocity at the interfaces occur naturally. The improved accuracy of the inverted instantaneous velocities is particularly important for accurate time‐to‐depth conversion.     

Improved control algorithm for real‐time substructure testing

Real‐time substructure testing is a novel method of testing structures under dynamic loading. The complete structure is separated into two substructures, one of which is tested physically at large scale and in real time, so that time‐dependent non‐linear behaviour of the substructure is realistically represented. The second substructure represents the surrounding structure, which is modelled numerically. In the current formulation this numerical substructure is assumed to remain linear. The two substructures interact in real‐time so that the response of the complete structure, incorporating the non‐linear behaviour of the physical substructure, is accurately represented. This paper presents several improvements to the linear numerical modelling of substructures for use in explicit time‐stepping routines for real‐time substructure testing. An extrapolation of a first‐order‐hold discretization is used which increases the accuracy of the numerical model over more direct explicit methods. Additionally, an integral form of the equation of motion is used in order to reduce the effects of noise and to take into account variations of the input over a time‐step. In order to take advantage of this integral form, interpolation of the model output is performed in order to smooth the output. The improvements are demonstrated using a series of substructure tests on a simple portal frame. While the testing approach is suitable for cases in which the physical substructure behaves non‐linearly, the results presented here are for fully linear systems. This enables comparisons to be made with analytical solutions, as well as with the results of tests based on the central difference method. Copyright © 2001 John Wiley & Sons, Ltd.     

基于改进的接收点插值算法的频率域海洋可控源电磁法2.5维正演

在海洋可控源电磁法勘探中,接收站常置于海底.在进行海洋电磁场模拟时,由于海水和海底介质存在显著电性差异,这给海底接收点处场值的求取带来困难.本文提出一种新的接收点插值算法,该算法考虑到海底电场法向分量不连续性问题,用法向电流分量进行插值以准确求取海底任意接收点处电磁场值.本文利用交错网格有限差分法实现了二维介质中频率域海洋可控源法(CSEM)正演.对构造走向做傅里叶变换,将三维电磁模拟问题转换为波数域2.5维问题,即三维场源激励下针对二维地电模型的电磁模拟问题.使用交错网格有限差分法,基于一次场/二次场分离方法导出波数域二次电场离散形式,并进一步求得波数域电磁场.采用本文提出的改进的插值算法可求得海底任意接收点处波数域电磁场,采用傅里叶逆变换对波数域电磁场进行积分可得到接收点处空间域电磁场.模型算例表明,与常规的线性插值和严格插值算法相比,本文提出的改进的插值算法具有更高的精度.     

A unified formulation of the piecewise exact method for inelastic seismic demand analysis including the P‐delta effect

The non‐linear analysis of single‐degree‐of‐freedom (SDOF) systems provides the essential background information for both strength‐based design and displacement‐based evaluation/design methodologies through the development of the inelastic response spectra. The recursive solution procedure called the piecewise exact method, which is efficiently used for the response analysis of linear SDOF systems, is re‐formulated in this paper in a unified format to analyse the non‐linear SDOF systems with multi‐linear hysteresis models. The unified formulation is also capable of handling the P‐delta effect, which generally involves the negative post‐yield stiffness of the hysteresis loops. The attractiveness of the method lies in the fact that it provides the exact solution when the loading time history is composed of piecewise linear segments, a condition that is perfectly satisfied for the earthquake excitation. Based on simple recursive relationships given for positive, negative and zero effective stiffnesses, the unified form of the piecewise exact method proves to be an extremely powerful and probably the best tool for the SDOF inelastic time‐history and response spectrum analysis including the P‐delta effect. A number of examples are presented to demonstrate the implementation of the method. Copyright © 2003 John Wiley & Sons, Ltd.     

Insight into the marine controlled‐source electromagnetic signal propagation

In marine controlled‐source electromagnetic (CSEM) surveys the subsurface is explored by emitting low‐frequency signals from an electric dipole source close to the sea‐bed. The main goal is often to detect and describe possible thin resistive layers beneath the sea‐bed. To gain insight into how CSEM signals propagate, it is informative to study a stratified model. The electromagnetic field is then given in terms of integrals over TE‐ and TM‐polarized plane‐wave constituents. An asymptotic evaluation of the field integrals for large propagation distances results in explicit spatial expressions for the field components and the derived expressions can be used to analyse how the CSEM signals propagate. There are two major signal pathways in a standard CSEM model. One of these pathways is via the thin resistive layer and the resulting response is accounted for by a pole in the reflection response for the TM mode. The signal is propagating nearly vertically down to the resistor from the source, then guided while attenuated along the reservoir, before propagating nearly vertically up to the receiver. The response is slightly altered by the sea‐bed interface and further modified in shallow water due to multiple reflections between the sea‐surface and sea‐bed at both the source and receiver sides. The other major signal pathway is via the resistive air half‐space, the so‐called airwave. The airwave is generated by the TE mode and interacts with the subsurface via vertically propagating signals reflected between the sea‐surface and subsurface at both the source and receiver sides.     

A 27‐point scheme for a 3D frequency‐domain scalar wave equation based on an average‐derivative method

Based on an average‐derivative method and optimization techniques, a 27‐point scheme for a 3D frequency‐domain scalar wave equation is developed. Compared to the rotated‐coordinate approach, the average‐derivative optimal method is not only concise but also applies to equal and unequal directional sampling intervals. The resulting 27‐point scheme uses a 27‐point operator to approximate spatial derivatives and the mass acceleration term. The coefficients are determined by minimizing phase velocity dispersion errors and the resultant optimal coefficients depend on ratios of directional sampling intervals. Compared to the classical 7‐point scheme, the number of grid points per shortest wavelength is reduced from approximately 13 to approximately 4 by this 27‐point optimal scheme for equal directional sampling intervals and unequal directional sampling intervals as well. Two numerical examples are presented to demonstrate the theoretical analysis. The average‐derivative algorithm is also extended to a 3D frequency‐domain viscous scalar wave equation.     

Regularized sparse‐grid geometric sampling for uncertainty analysis in non‐linear inverse problems

This paper introduces an efficiency improvement to the sparse‐grid geometric sampling methodology for assessing uncertainty in non‐linear geophysical inverse problems. Traditional sparse‐grid geometric sampling works by sampling in a reduced‐dimension parameter space bounded by a feasible polytope, e.g., a generalization of a polygon to dimension above two. The feasible polytope is approximated by a hypercube. When the polytope is very irregular, the hypercube can be a poor approximation leading to computational inefficiency in sampling. We show how the polytope can be regularized using a rotation and scaling based on principal component analysis. This simple regularization helps to increase the efficiency of the sampling and by extension the computational complexity of the uncertainty solution. We demonstrate this on two synthetic 1D examples related to controlled‐source electromagnetic and amplitude versus offset inversion. The results show an improvement of about 50% in the performance of the proposed methodology when compared with the traditional one. However, as the amplitude versus offset example shows, the differences in the efficiency of the proposed methodology are very likely to be dependent on the shape and complexity of the original polytope. However, it is necessary to pursue further investigations on the regularization of the original polytope in order to fully understand when a simple regularization step based on rotation and scaling is enough.     

Non‐iterative time integration scheme for non‐linear dynamic FEM analysis

This paper proposes a non‐iterative time integration (NITI) scheme for non‐linear dynamic FEM analysis. The NITI scheme is constructed by combining explicit and implicit schemes, taking advantage of their merits, and enables stable computation without an iteration process for convergence even when used for non‐linear dynamic problems. Formulation of the NITI scheme is presented and its stability is studied. Although the NITI scheme is not unconditionally stable when applied to non‐linear problems, it is stable in most cases unless stiffness hardening occurs or the problem has a large velocity‐dependent term. The NITI scheme is applied to dynamic analysis of the non‐linear soil–structure system and computation results are compared with those by the central difference method (CDM). Comparison shows that the stability of the NITI scheme is superior to that of the CDM. Accuracy of the NITI scheme is verified because its results are identical with those by the CDM in which the time step is set as 1/10 of that for the NITI scheme. The application of the NITI scheme to the mesh‐partitioned FEM is also proposed. It is applied to dynamic analysis of the linear soil–structure system. It yields the same results as a conventional single‐domain FEM analysis using the Newmark β method. This result verifies the usability of mesh‐partitioned FEM analysis using the NITI scheme. Copyright © 2003 John Wiley& Sons, Ltd.     

Neuro‐control of seismically excited steel structure through sensitivity evaluation scheme

The neuro‐controller training algorithm based on cost function is applied to a multi‐degree‐of‐freedom system; and a sensitivity evaluation algorithm replacing the emulator neural network is proposed. In conventional methods, the emulator neural network is used to evaluate the sensitivity of structural response to the control signal. To use the emulator, it should be trained to predict the dynamic response of the structure. Much of the time is usually spent on training of the emulator. In the proposed algorithm, however, it takes only one sampling time to obtain the sensitivity. Therefore, training time for the emulator is eliminated. As a result, only one neural network is used for the neuro‐control system. In the numerical example, the three‐storey building structure with linear and non‐linear stiffness is controlled by the trained neural network. The actuator dynamics and control time delay are considered in the simulation. Numerical examples show that the proposed control algorithm is valid in structural control. Copyright © 2001 John Wiley & Sons, Ltd.     

基于拟态有限体积法的频率域可控源三维正演计算

大规模地球物理电磁数据的定量解释需要发展高效、稳定的三维正反演算法.本文通过求解离散化的三维电场矢量Helmholtz方程,实现了基于有限体积法的频率域可控源电磁(CSEM)三维正演算法.为模拟具有强电性差异的三维电性介质,该算法采用拟态有限体积法(MFV)对Maxwell方程组进行离散化;另外,为获得稳定、高精度的正演数值结果,采用直接矩阵分解技术来求解离散所得到的大型稀疏线性方程组.对于具有多个发射源的CSEM测量来说,一次矩阵分解结果能够用于同频率下所有场源的正演计算.为降低场源奇异性及边界条件对数值精度的影响,采用虚拟场源校正技术,避免了散射场公式中在构建场源项时所需的大量时间.对于具有多个频率的CSEM的模拟计算,采用分频并行策略来加快三维正演计算.最后,通过与一维层状模型及三维模型的数值结果的对比验证了本文所开发的正演算法对频率域CSEM模拟计算的准确性及有效性,表明该正演算法能够有效应用于三维介质的数值计算.另外,对于多频率CSEM的并行测试结果表明基于分频并行策略的并行计算能够显著地降低正演计算时间.     

Seismic data interpolation using generalised velocity‐dependent seislet transform

Data interpolation is an important step for seismic data analysis because many processing tasks, such as multiple attenuation and migration, are based on regularly sampled seismic data. Failed interpolations may introduce artifacts and eventually lead to inaccurate final processing results. In this paper, we generalised seismic data interpolation as a basis pursuit problem and proposed an iteration framework for recovering missing data. The method is based on non‐linear iteration and sparse transform. A modified Bregman iteration is used for solving the constrained minimisation problem based on compressed sensing. The new iterative strategy guarantees fast convergence by using a fixed threshold value. We also propose a generalised velocity‐dependent formulation of the seislet transform as an effective sparse transform, in which the non‐hyperbolic normal moveout equation serves as a bridge between local slope patterns and moveout parametres in the common‐midpoint domain. It can also be reduced to the traditional velocity‐dependent seislet if special heterogeneity parametre is selected. The generalised velocity‐dependent seislet transform predicts prestack reflection data in offset coordinates, which provides a high compression of reflection events. The method was applied to synthetic and field data examples, and the results show that the generalised velocity‐dependent seislet transform can reconstruct missing data with the help of the modified Bregman iteration even for non‐hyperbolic reflections under complex conditions, such as vertical transverse isotropic (VTI) media or aliasing.     

Erratic noise suppression using iterative structure‐oriented space‐varying median filtering with sparsity constraint

Erratic noise often has high amplitudes and a non‐Gaussian distribution. Least‐squares–based approaches therefore are not optimal. This can be handled better with non–least‐squares approaches, for example based on Huber norm which is computationally expensive. An alternative method has been published which involves transforming the data with erratic noise to pseudodata that have Gaussian distributed noise. It can then be attenuated using traditional least‐squares approaches. This alternative method has previously been used in combination with a curvelet transform in an iterative scheme. In this paper, we introduce a median‐filtering step in this iterative scheme. The median filter is applied following the slope direction of the seismic data to maximally preserve the energy of useful signals. The new method can suppress stronger erratic noise compared with the previous iterative method, and can better deal with random noise compared with the single‐step implementation of the median filter. We apply the proposed robust denoising algorithm to a synthetic dataset and two field data examples and demonstrate its advantages over three different noise attenuation algorithms.     

Earthquake response of elastic SDF systems with non‐linear fluid viscous dampers

The steady‐forced and earthquake responses of SDF systems with a non‐linear fluid viscous damper (FVD) are investigated. The energy dissipation capacity of the FVD is characterized by the supplemental damping ratio ζsd and its non‐linearity by a parameter designated α. It is found that the structural response is most effectively investigated in terms of ζ_sd and α because (1) these two parameters are dimensionless and independent, and (2) the structural response varies linearly with the excitation intensity. Damper non‐linearity has essentially no influence on the peak response of systems in the velocity‐sensitive spectral region, but differences up to 14% were observed in the other spectral regions. The structural deformation is reduced by up to 25% when ζ_sd= 5%; and by up to 60% when ζ_sd= 30%. Non‐linear FVDs are advantageous because they achieve essentially the same reduction in system responses but with a significantly reduced damper force. For practical applications, a procedure is presented to estimate the design values of structural deformation and forces for a system with non‐linear FVD directly from the design spectrum. It is demonstrated that the earthquake‐induced force in a non‐linear FVD can be estimated from the damper force in a corresponding system with linear FVD, its peak deformation, and peak relative velocity; however, the relative velocity should not be approximated by the pseudo‐velocity as this approximation introduces a large error in the damper force. Finally, a procedure is presented to determine the non‐linear damper properties necessary to limit the structural deformation to some design value or the structural capacity for a given design spectrum. Copyright © 2002 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.     

Non‐linear‐Maxwell‐element‐type hysteretic control force

A simple non‐linear control law is proposed for reducing structural responses against seismic excitations. This law defines control force dynamics by one differential equation involving a non‐linear term that restrains the control force amplitude. If non‐linearity is neglected, the control force becomes the force in a Maxwell element, so it is called the non‐linear‐Maxwell‐element‐type (NMW) control force. The NMW control force vs. deformation relation plots hysteretic curves. The basic performance of an SDOF model with the NMW control force is examined for various conditions by numerical analyses. Furthermore, the control law is extended to fit an MDOF structural model, and an application example is shown. The computational results show that the NMW control force efficiently reduces structural responses. Copyright © 2000 John Wiley & Sons, Ltd.     

Adsorption Kinetic Modeling of Safranin onto Rice Husk Biomatrix Using Pseudo‐first‐ and Pseudo‐second‐order Kinetic Models: Comparison of Linear and Non‐linear Methods

Batch kinetic studies were carried out for the removal of safranin from aqueous solution using a biomatrix prepared from rice husk. The adsorption kinetic data were modeled using the pseudo‐first‐order and pseudo‐second‐order kinetic equations. The linear and non‐linear forms of these two widely used kinetic models were compared in this study. In order to determine the best‐fitting equation, the coefficient of determination (r²), the sum of the squares of the errors (SSE), sum of the absolute errors (SAE), average relative error (ARE), hybrid fractional error function (HYBRID), Marquardt's percent standard deviation (MPSD), and the Chi‐squared test (χ²) were used as error analysis methods. Results showed that the non‐linear forms of pseudo‐first‐order and pseudo‐second‐order models were more suitable than the linear forms for fitting the experimental data. Non‐linear method is thus more appropriate for estimating the kinetic parameters and should primarily be used to describe adsorption kinetics.     

Automatic stacking‐velocity estimation using similarity‐weighted clustering

Local seismic event slopes contain subsurface velocity information and can be used to estimate seismic stacking velocity. In this paper, we propose a novel approach to estimate the stacking velocity automatically from seismic reflection data using similarity‐weighted k‐means clustering, in which the weights are local similarity between each trace in common midpoint gather and a reference trace. Local similarity reflects the local signal‐to‐noise ratio in common midpoint gather. We select the data points with high signal‐to‐noise ratio to be used in the velocity estimation with large weights in mapped traveltime and velocity domain by similarity‐weighted k‐means clustering with thresholding. By using weighted k‐means clustering, we make clustering centroids closer to those data points with large weights, which are more reliable and have higher signal‐to‐noise ratio. The interpolation is used to obtain the whole velocity volume after we have got velocity points calculated by weighted k‐means clustering. Using the proposed method, one obtains a more accurate estimate of the stacking velocity because the similarity‐based weighting in clustering takes into account the signal‐to‐noise ratio and reliability of different data points in mapped traveltime and velocity domain. In order to demonstrate that, we apply the proposed method to synthetic and field data examples, and the resulting images are of higher quality when compared with the ones obtained using existing methods.     

Shaking table 2‐D models of a concrete gravity dam

One of the most famous and studied cases of dams subjected to earthquake loading is the Koyna Dam in India. In this study, a two‐dimensional model of Koyna Dam at 1/50 scale was used on a shake table to simulate effects and serve as data for non‐linear computer model calibration. A new concrete mix was designed for the non‐linear similitude modelling. This new mix provided the correct kinematic failure of concrete at scale. Two models were tested to failure: one with an initial shrinkage crack and one monolith. Reservoir effects were not modelled. The results of both models are discussed and compared. The ability to model non‐linear effects is discussed. Published in 2000 by John Wiley & Sons, Ltd.