Vertical image waves in elliptically anisotropic media

By reparameterization of the kinematic expressions for remigration in elliptically anisotropic media using a new ellipticity parameter, we derive a new image-wave equation in elliptically anisotropic media describing the position of the reflector as a function of the medium ellipticity. This image wave equation, which is a kind of medium-dependent one-way wave equation, can be used for automatically stretching a time-migrated image in depth until wells are tied or other given geologic criteria are met. In this way, it provides a useful means to use a priori depth information for finding an estimate of the vertical velocity, which cannot be detected from time processing only. Simple numerical examples confirm the validity of the image-wave equation even for nonconstant velocity.     

A fast algorithm for 3D azimuthally anisotropic velocity scan

The conventional velocity scan can be computationally expensive for large‐scale seismic data sets, particularly when the presence of anisotropy requires multiparameter scanning. We introduce a fast algorithm for 3D azimuthally anisotropic velocity scan by generalizing the previously proposed 2D butterfly algorithm for hyperbolic Radon transforms. To compute semblance in a two‐parameter residual moveout domain, the numerical complexity of our algorithm is roughly as opposed to of the straightforward velocity scan, with N being the representative of the number of points in a particular dimension of either data space or parameter space. Synthetic and field data examples demonstrate the superior efficiency of the proposed algorithm.     

Variational Data Assimilation Technique in Mathematical Modeling of Ocean Dynamics

Problems of the variational data assimilation for the primitive equation ocean model constructed at the Institute of Numerical Mathematics, Russian Academy of Sciences are considered. The model has a flexible computational structure and consists of two parts: a forward prognostic model, and its adjoint analog. The numerical algorithm for the forward and adjoint models is constructed based on the method of multicomponent splitting. The method includes splitting with respect to physical processes and space coordinates. Numerical experiments are performed with the use of the Indian Ocean and the World Ocean as examples. These numerical examples support the theoretical conclusions and demonstrate the rationality of the approach using an ocean dynamics model with an observed data assimilation procedure.     

Nonlinear dynamic analysis of multi-base seismically isolated structures with uplift potential II： verification examples

The work presented in this paper serves as numerical verification of the analytical model developed in the companion paper for nonlinear dynamic analysis of multi-base seismically isolated structures. To this end, two numerical examples have been analyzed using the computational algorithm incorporated into program 3D-BASIS-ME-MB, developed on the basis of the newly-formulated analytical model. The first example concerns a seven-story model structure that was tested on the earthquake simulator at the University at Buffalo and was also used as a verification example for program SAP2000. The second example concerns a two-tower, multi-story structure with a split-level seismic-isolation system. For purposes of verification, key results produced by 3D-BASIS-ME-MB are compared to experimental results, or results obtained from other structural/finite element programs. In both examples, the analyzed structure is excited under conditions of bearing uplift, thus yielding a case of much interest in verifying the capabilities of the developed analysis tool.     

场地土层模型参数的地震动记录反演方法

考虑场地土层不同物理参数的综合反演,改进了水平与竖向谱比(Horizontal-to-Vertical Spectral Ratio,简称HVSR)混合全局优化反演方法,进一步通过土层的S、P波波速、厚度、泊松比、密度和S、P波品质因子等土层参数反演的敏感性分析,形成了可同时反演场地土层厚度及剪切波速的混合全局优化反演方法.以美国GVDA和日本IWTH27竖向强震动观测台阵为例,分别以理论HVSR及加速度观测记录获得的HVSR曲线为目标,反演获得了场地浅层速度结构,并与观测台阵场地钻孔揭示的土层模型进行比较,验证了发展的反演方法的合理性和适用性.本文研究表明,基于地震加速度记录的HVSR全局优化反演方法是获取场地土层速度结构的一种有效的途径.     

Pattern-guided dip estimation with plane-wave destruction filters

We propose to use pattern-guided dip estimation to mitigate aliasing problem that possibly exists in structure-oriented data processing. A straightforward and effective approach of generating pattern-guided dip is presented, which generally involves three rounds of standard dip estimation with plane-wave destruction filters. The first use of plane-wave destruction filter is for generating a mask operator distinguishing aliased and non-aliased data, based on measuring the uncertainty of the first dip estimation. The second plane-wave destruction filter uses the aliasing-free portions of the input data, and the dip in the aliasing-affected area is automatically padded with the ‘pattern’ dip by smoothing regularization. The result of the second plane-wave destruction filter is used as the initial dip for the inversion of the last-pass plane-wave destruction filter, which produces a pattern-guided dip. For some specific applications, the mask operator can be easily generated through other methods, and we can skip the first dip estimation. Two numerical examples, related to picking information using predictive painting and structure-oriented interpolation, respectively, demonstrate that our pattern-guided dip can effectively mitigate the aliasing problem in structure-oriented data processing.     

瑞雷波多阶模式频散曲线稀疏正则化反演方法研究

目前瑞雷波多阶模式频散曲线反演中仅考虑数据的拟合,缺乏对模型的约束,不能很好地刻画地层间断面的问题,针对此问题,研究了瑞雷波多阶模式频散曲线稀疏正则化反演方法.正演模拟基于广义反射-透射系数法,数值计算上采用一种快速求根方法,与二等分方法相比,能够在很短的时间内达到最优的收敛效果;反演建模时采用L1范数正则化方法对模型...     

NUMERCAL COMPUTATION OF VERTICAL TWO-DIMENSIONAL NON-EQUILIBRIUM SEDIMENT MOVEMENT

INON-EQUILIBRIUMMOVEMENTContinualexchangeremainsamongthebedmaterial,bedloadandsuspendedloadinariver.Thebedloadmovesonthebedsurfaceandjumpsinsteps.Itstransportrateperwidthvarieswithvaryingflowintensity;whilethesuspendedloadmoveswithalongstep,evenifthehydraulicfactorsbecomeweaker,itwillnotretUrntothebeduntilfinishingthefallingprocess.ThismeansthatahysteresisexistsbetWeenthechangeofthesuspendedsedimentmovementandflowvelocity.Foruniformandsteadyflow,thesedimentmovementkeepsinequilibriumi…     

无网格局部Petrov-Galerkin法大地电磁场二维正演模拟

有限差分法和有限单元法在大地电磁场数值模拟中已经得到了广泛的应用,但其数值结果的精度在很大程度上依赖于网格的离散程度.当模拟起伏地形、弯曲界面等复杂地电模型大地电磁场响应时,常常需要花费大量的时间以便得到较合理的离散网格.无网格局部Petrov-Galerkin法(MLPG)不同于有限差分法和有限元法,其形函数和权函数脱离了网格的束缚.本文详细推导了二维大地电磁场边值问题的弱式形式,并将其离散为局部积分域内的表达形式.通过模拟二维海洋地电模型大地电磁场响应,并与结构网格有限元结果进行对比,验证了本文算法和程序的正确性及精度.设计了一个含有弯曲界面的二维地电模型,讨论了不同离散网格对MLPG无网格法模拟结果的影响,并与结构有限元法结果进行了比较,结果表明MLPG无网格法模拟结果受离散网格影响较小.最后利用MLPG无网格法计算了两个海洋起伏地形模型的大地电磁响应,讨论了海底起伏地形对大地电磁响应的影响.     

接收函数确定Moho面速度和密度跃变的方法研究及应用

通过对单层模型反射和透射系数的推导,提出了利用接收函数一次转换波和多次波确定Moho面速度和密度跃变的速度-密度跃变(δβ-δρ)扫描叠加方法.利用反射率法计算了不同模型的远震理论地震图,按照与处理实际观测波形一致的方法和流程计算了理论接收函数;根据不同模型数值试验结果,深入分析了界面速度和密度跃变对接收函数震相幅度的影响.利用(δβ-δρ)扫描叠加方法,对理论接收函数进行了数值试验,结果证明了该方法的可行性.最后将该方法应用于位于青藏高原东北缘的高台(GTA)台和兰州(LZH)台,确定了两个台站下方Moho面的速度跃变分别约为(19±1)%和(20±1)%,密度跃变最小值为(4±2)%和(6±2)%.     

Subsurface offset behaviour in velocity analysis with extended reflectivity images

Migration velocity analysis with the constant‐density acoustic wave equation can be accomplished by the focusing of extended migration images, obtained by introducing a subsurface shift in the imaging condition. A reflector in a wrong velocity model will show up as a curve in the extended image. In the correct model, it should collapse to a point. The usual approach to obtain a focused image involves a cost functional that penalizes energy in the extended image at non‐zero shift. Its minimization by a gradient‐based method should then produce the correct velocity model. Here, asymptotic analysis and numerical examples show that this method may be too sensitive to amplitude peaks at large shifts at the wrong depth and to artefacts. A more robust alternative is proposed that can be interpreted as a generalization of stack power and maximizes the energy at zero‐subsurface shift. A real‐data example is included.     

Offset continuation (OCO) ray tracing using OCO trajectories

Offset continuation (OCO) is a seismic configuration transform designed to simulate a seismic section as if obtained with a certain source-receiver offset using the data measured with another offset. Since OCO is dependent on the velocity model used in the process, comparison of the simulated section to an acquired section allows for the extraction of velocity information. An algorithm for such a horizon-oriented velocity analysis is based on so-called OCO rays. These OCO rays describe the output point of an OCO as a function of the Root Mean Square (RMS) velocity. The intersection point of an OCO ray with the picked traveltime curve in the acquired data corresponding to the output half-offset defines the RMS velocity at that position. We theoretically relate the OCO rays to the kinematic properties of OCO image waves that describe the continuous transformation of the common-offset reflection event from one offset to another. By applying the method of characteristics to the OCO image-wave equation, we obtain a raytracing-like procedure that allows to construct OCO trajectories describing the position of the OCO output point under varying offset. The endpoints of these OCO trajectories for a single input point and different values of the RMS velocity form then the OCO rays. A numerical example demonstrates that the developed ray-tracing procedure leads to reliable OCO rays, which in turn provide high-quality RMS velocities. The proposed procedure can be carried out fully automatically, while conventional velocity analysis needs human intervention. Moreover, since velocities are extracted using offset sections, more redundancy is available or, alternatively, OCO velocities can be studied as a function of offset.     

Energy-consistent integration method and its application to hybrid testing

The time-integration algorithm is an indispensable element to determine response of the boundary of the numerical as well as physical parts in a hybrid test. Instability of the time-integration algorithm may directly lead to failure of the test, so stability of an integration algorithm is particularly important for hybrid testing. The explicit algorithms are very popular in hybrid testing, because iteration is not needed. Many unconditionally stable explicit-algorithms have been proposed for hybrid testing. However, the stability analysis approaches used in all these methods are valid only for linear systems. In this paper, a uniform formulation for energy-consistent time integrations, which are unconditionally stable, is proposed for nonlinear systems. The solvability and accuracy are analyzed for typical energy-consistent algorithms. Some numerical examples and the results of a hybrid test are provided to validate the effectiveness of energy-consistent algorithms.     

地震偏移分裂算法的稳定性

本文讨论用分裂算法解高阶偏移方程时所用的有限差分方程与原偏微分方程的相容性和算法的稳定性。根据Lax等价原理,这种相容性和稳定性可以保证数值方程收敛于原偏微分方程。我们证明了该分裂算法满足相容性和稳定性要求。这就不但以文献[1]中的实例说明了这种算法的实用性,而且从原理上论证了算法的正确性。     

过套管电阻率测井的地层电阻率快速反演算法

过套管电阻率测井是通过测量套管壁电势实现测量地层的视电阻率,基于传输线方程理论,针对层状地层,给出了套管壁电势、电流对地层横向电阻导数的微分方程(称Jacobi矩阵微分方程)及边界条件;利用Jacobi矩阵微分方程边值问题导出了过套管电阻率测井反演地层参数的Jacobi矩阵系数的解析表示,利用Marquardt方法实现了过套管测井的地层电阻率反演;通过计算对Jacobi矩阵的特性进行了探讨,并获得了较快的计算速度(因为Jacobi矩阵是用解析解表示的),反演结果与地层模型取得了较好的逼近.本文实现了过套管电阻率测井地层参数的Jacobi系数矩阵的快速计算及地层电阻率反演,为进一步开展电阻率测井数据处理提供了理论依据和快速反演算法.     

利用走时反演二维多孔介质渗透率

综合利用直接方法和优化方法研究二维多孔介质渗透率反演问题. 问题的提法与前人不同,是由区域内各点流体渗透的走时来反演渗透率. 该反问题的求解可分成两步进行,归结为两个相应的子问题进行研究. 首先由流体走时反演二维多孔介质速度场的数值算法,然后由二维多孔介质速度场反演二维多孔介质渗透率的数值算法,最后给出数值例子. 数值结果表明了所用数值方法的正确性和有效性.     

Optimal placement of dampers for passive response control

The effectiveness of viscous and viscoelastic dampers for seismic response reduction of structures is quite well known in the earthquake engineering community. This paper deals with the optimal utilization of these dampers in a structure to achieve a desired performance under earthquake‐induced ground excitations. Frequency‐dependent and ‐independent viscous dampers and viscoelastic dampers have been considered as the devices of choice. To determine the optimal size and location of these dampers in the structure, a genetic algorithm is used. The desired performance is defined in terms of several different forms of performance functions. The use of the genetic approach is not limited to any particular form of performance function as long as it can be calculated numerically. For illustration, numerical examples for different building structures are presented showing the distribution and size of different dampers required to achieve a desired level of reduction in the response or a performance index. Copyright © 2002 John Wiley & Sons, Ltd.     

Triplications for the converted wave in transversely isotropic media with a tilted symmetry axis

In seismic data processing, serious problems could be caused by the existence of triplication and need to be treated properly for tomography and other inversion methods. The triplication in transversely isotropic medium with a vertical symmetry axis has been well studied and concluded that the triplicated traveltime only occurs for S wave and there is no triplication for P and converted PS waves since the P wave convexity slowness always compensates the S wave slowness concavity. Compared with the vertical symmetry axis model, the research of the triplication in transversely isotropic medium with a tilted symmetry axis is still keeping blank. In order to analyse the triplication for the converted wave in the tilted symmetry axis model, we examine the traveltime of the triplication from the curvature of averaged P and S wave slowness. Three models are defined and tested in the numerical examples to illustrate the behaviour of the tilted symmetry axis model for the triplicated traveltime with the change of the rotation angle. Since the orientation of an interface is related to the orientation of the symmetry axis, the triplicated traveltime is encountered for the converted wave in the tilted symmetry axis model assuming interfaces to be planar and horizontal. The triplicated region is influenced by the place and level of the concave curvature of the P and S wave slowness.     

Implementation aspects of eigendecomposition‐based high‐resolution velocity spectra

In this paper, we discuss high‐resolution coherence functions for the estimation of the stacking parameters in seismic signal processing. We focus on the Multiple Signal Classification which uses the eigendecomposition of the seismic data to measure the coherence along stacking curves. This algorithm can outperform the traditional semblance in cases of close or interfering reflections, generating a sharper velocity spectrum. Our main contribution is to propose complexity‐reducing strategies for its implementation to make it a feasible alternative to semblance. First, we show how to compute the multiple signal classification spectrum based on the eigendecomposition of the temporal correlation matrix of the seismic data. This matrix has a lower order than the spatial correlation used by other methods, so computing its eigendecomposition is simpler. Then we show how to compute its coherence measure in terms of the signal subspace of seismic data. This further reduces the computational cost as we now have to compute fewer eigenvectors than those required by the noise subspace currently used in the literature. Furthermore, we show how these eigenvectors can be computed with the low‐complexity power method. As a result of these simplifications, we show that the complexity of computing the multiple signal classification velocity spectrum is only about three times greater than semblance. Also, we propose a new normalization function to deal with the high dynamic range of the velocity spectrum. Numerical examples with synthetic and real seismic data indicate that the proposed approach provides stacking parameters with better resolution than conventional semblance, at an affordable computational cost.     

Seismic reliability assessment of lifeline networks using clustering‐based multi‐scale approach

Seismic reliability assessment of lifeline networks gives rise to various technical challenges, which are mostly caused by a large number of network components, complex network topology, and statistical dependence between component failures. For effective risk assessment and probabilistic inference based on post‐hazard observations, various non‐simulation‐based algorithms have been developed, including the selective recursive decomposition algorithm (S‐RDA). To facilitate the application of such an algorithm to large networks, a new multi‐scale approach is developed in this paper. Using spectral clustering algorithms, a network is first divided into an adequate number of clusters such that the number of inter‐cluster links is minimized while the number of the nodes in each cluster remains reasonably large. The connectivity around the identified clusters is represented by super‐links. The reduced size of the simplified network enables the S‐RDA algorithm to perform the network risk assessment efficiently. When the simplified network is still large even after a clustering, additional levels of clustering can be introduced to have a hierarchical modeling structure. The efficiency and effectiveness of the proposed multi‐scale approach are demonstrated successfully by numerical examples of a hypothetical network, a gas transmission pipeline network, and a water transmission network. Copyright © 2014 John Wiley & Sons, Ltd.