Generalized Born inversion of seismic reflection data

Summary . Born inverse methods give accurate and stable results when the source wavelet is impulsive. However, in many practical applications (reflection seismology) an impulsive source cannot be realized and the inversion needs to be generalized to include an arbitrary source function. In this paper, we present a Born solution to the seismic inverse problem which can accommodate an arbitrary source function and give accurate and stable results. It is shown that the form of the generalized inversion algorithm reduces to a Wiener shaping ***filter, which is solved efficiently using a Levinson recursion algorithm. Numerical examples of synthetic and real field data illustrate the validity of our method.     

Time-lapse traveltime change of singly scattered acoustic waves

We present a technique based on the single-scattering approximation that relates time-lapse localized changes in the propagation velocity to changes in the traveltime of singly scattered waves. We describe wave propagation in a random medium with homogeneous statistical properties as a single-scattering process where the fluctuations of the velocity with respect to the background velocity are assumed to be weak. This corresponds to one of two end-member regimes of wave propagation in a random medium, the first being single scattering, and the second multiple scattering. We present a formulation that relates the change in the traveltime of the scattered waves to a localized change in the propagation velocity by means of the Born approximation for the scattered wavefield. We validate the methodology with synthetic seismograms calculated with finite differences for 2-D acoustic waves. Potential applications of this technique include non-destructive evaluation of heterogeneous materials and time-lapse monitoring of heterogeneous reservoirs.     

Backscattering from spherical elastic inclusions and accuracy of the Kirchhoff approximation for curved interfaces

The Kirchhoff (or tangent plane) approximation, derived from the theoretically complete Kirchhoff–Helmholtz integral representation for the seismic wavefield, has been used extensively for the analysis of seismic-wave scattering from irregular interfaces; however, the accuracy of this method for curved interfaces has not been rigorously established. This paper describes an efficient Kirchhoff algorithm to simulate scattered waves from an arbitrarily curved interface in an elastic medium. Synthetic seismograms computed using this algorithm are compared with exact synthetics computed using analytical formulae for scattering of plane P waves by a spherical elastic inclusion. A windowing technique is used to remove strong internal reverberations from the analytical solution. Although the Kirchhoff method tends to underestimate the total scattering intensity, the accuracy of the approximation improves with increasing value of the wavenumber-radius product, kR . The arrival times and pulse shapes of primary reflections from the sphere are well approximated using the Kirchhoff approach regardless of curvature of the scattering surface, but the amplitudes are significantly underestimated for kR ≤ 5. The results of this work provide some new guidelines to assess the accuracy of Kirchhoff-synthetic seismograms for curved interfaces.     

Construction and computation of geometries on the ellipsoid

This article describes how to construct a wide range of geometry objects (called GeographicGeometry objects) in the coordinate system of an ellipsoid such as the Geographic coordinate system. Each construction process is formulated analytically and algorithmically using a combination of a set of fairly well-known mathematical methods such as ellipsoid geodesic construction functions, spherical trigonometry and iterative refinement methods. Each such geometry object may efficiently be converted to a corresponding Cartesian geometry object in any map projection coordinate system using an approximation algorithm. This property makes them particularly useful as a coordinate-system-independent geometry representation. A geographic geometry object is normally topologically equivalent to its Cartesian geometry counterpart except for some discontinuity and singularity cases.     

Improved version of BTOPMC model and its application in event-based hydrologic simulations

In this paper, a grid-based distributed hydrological model BTOPMC (Block-wise use of TOPMODEL) is introduced, which was developed from the original TOPMODEL. In order to broaden the model’s application to arid regions, improvement methodology is also implemented. The canopy interception and soil infiltration processes were incorporated into the original BTOPMC to model event-based runoff simulation in large arid regions. One de-signed infiltration model with application of time compression approximation method is emphasized and validated for improving model’s performance for event hydrological simulations with a case study of Lushi River basin.     

基于GIS的任意发生元Voronoi图逼近方法

许多地理问题的空间分析中需要采用Voronoi图,但是目前我们尚缺乏一些简单的易于实现的构建任意发生元Voronoi图的方法,也缺乏一个能直接生成任意发生元Voronoi图的软件,为此我们提出了一种基于GIS的构建任意发生元的未加权Voronoi图的逼近方法。首先用有限点来逼近原始发生元,然后构建这些点发生元Voronoi图,最后消除那些属于同一发生元的顶点和边,即得到原始发生元的逼近的Voronoi图。在该算法的具体实现过程中,充分利用了现有GIS软件可以生成点发生元Voronoi图的特性和处理空间数据的能力。试验结果表明,这种方法可以生成未加权的任意形状发生元的逼近Voronoi图,能满足地理问题空间分析的需要。如地理客体可以是点状地理客体(城市、县城、交通枢纽、商业中心和金融中心等)、线状地理客体(交通运输线、经济地带和河系等)、面状地理客体(经济区、公园和绿地等)或者它们的组合,它们的空间影响范围或空间服务范围都可以采用Voronoi图来界定。     

沙尘暴中球形沙粒局部带电对电磁波的交叉去极化效应

 基于小粒子散射的Rayleigh近似,给出局部带电小圆球粒子在不同极化方式入射波作用下的前向散射振幅函数,以及差分衰减、差分相移表达式,同时计算了线极化毫米波在带电沙尘暴中传播的交叉去极化分辨率。结果表明,圆球粒子局部带电时也会使得入射极化波产生去极化效应,且入射波的去极化效应随着电荷分布角的增加先增强后减弱,在90°角达到最大值;随着荷质比、入射波频率的增加不断增强,随着能见度的增加不断减小。结果说明,在沙尘暴中的无线电接收装置设计时必须考虑到沙粒带电对通信信号的去极化影响。     

改进的BTOPMC模型及其在水文模拟中的应用

In this paper, a grid-based distributed hydrological model BTOPMC （Block-wise use of TOPMODEL） is introduced, which was developed from the original TOPMODEL. In order to broaden the model＇s application to arid regions, improvement methodology is also implemented. The canopy interception and soil infiltration processes were incorporated into the original BTOPMC to model event-based runoff simulation in large arid regions. One designed infiltration model with application of time compression approximation method is emphasized and validated for improving model＇s performance for event hydrological simulations with a case study of Lushi River basin.     

Split-step complex Padé-Fourier depth migration

We present a split-step complex Padé-Fourier migration method based on the one-way wave equation. The downward-continuation operator is split into two downward-continuation operators: one operator is a phase-shift operator and the other operator is a finite-difference operator. A complex treatment of the propagation operator is applied to mitigate inaccuracies and instabilities due to evanescent waves. It produces high-quality images of complex structures with fewer numerical artefacts than those obtained using a real approximation of a square-root operator in the one-way wave equation. Tests on zero-offset data from the SEG/EAGE salt data show that the method improves the image quality at the cost of an additional 10 per cent computational time compared to the conventional Fourier finite-difference method.