Synthetic seismograms in heterogeneous media by one-return approximation

When reverberations between heterogeneities or resonance scattering can be neglected but accumulated effects of forward scattering are strong, the Born approximation is not valid but the De Wolf approximation can be applied in such cases. In this paper, renormalized MFSB (multiple-forescattering single-backscattering) equations and the dual-domain expression for scalar, acoustic and elastic waves are derived by a unified approach. Two versions of the one-return method (using MFSB approximation) are given: One is the wide-angle dual-domain formulation (thin-slab approximation); the other is the screen approximation. In the screen approximation, which involves a small-angle approximation for the wave-medium interaction, it can be seen clearly that the forward scattered, or transmitted waves are mainly controlled by velocity perturbations; while the backscattered or reflected waves, by impedance perturbations. The validity of the method and the wide-angle capability of the dual-domain implementation are demonstrated by numerical examples. Reflection coefficients of a plane interface derived from numerical simulations by the wide-angle method match the theoretical curves well up to critical angles. For the reflections of a low-velocity slab, the agreement between theory and synthetics only starts to deteriorate for angles greater than 70°. The accuracy of the wide-angle version of the method could be further improved by optimizing the wave-number filtering for the forward propagation and shrinking the step length along the propagation direction.     

A Comparison of Two Canopy Radiative Models in Land Surface Processes

This paper compares the predictions by two radiative transfer models-the two-stream approximation model and the generalized layered model （developed by the authors） in land surface processes -for different canopies under direct or diffuse radiation conditions. The comparison indicates that there are significant differences between the two models, especially in the near infrared （NIR） band. Results of canopy reflectance from the two-stream model are larger than those from the generalized model. However, results of canopy absorptance from the two-stream model are larger in some cases and smaller in others compared to those from the generalized model, depending on the cases involved. In the visible （VIS） band, canopy reflectance is smaller and canopy absorptance larger from the two-stream model compared to the generalized model when the Leaf Area Index （LAI） is low and soil reflectance is high. In cases of canopies with vertical leaf angles, the differences of reflectance and absorptance in the VIS and NIR bands between the two models are especially large. Two commonly occurring cases, with which the two-stream model cannot deal accurately, are also investigated. One is for a canopy with different adaxial and abaxial leaf optical properties; and the other is for incident sky diffuse radiation with a non-uniform distribution. Comparison of the generalized model within the same canopy for both uniform and non-uniform incident diffuse radiation inputs shows smaller differences in general. However, there is a measurable difference between these radiation inputs for a canopy with high leaf angle. This indicates that the application of the two-stream model to a canopy with different adaxial and abaxial leaf optical properties will introduce non-negligible errors.     

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.     

Spatial derivatives and perturbation derivatives of amplitude in isotropic and anisotropic media

Explicit equations for the spatial derivatives and perturbation derivatives of amplitude in both isotropic and anisotropic media are derived. The spatial and perturbation derivatives of the logarithm of amplitude can be calculated by numerical quadratures along the rays. The spatial derivatives of amplitude may be useful in calculating the higher-order terms in the ray series, in calculating the higher-order amplitude coefficients of Gaussian beams, in estimating the accuracy of zero-order approximations of both the ray method and Gaussian beams, in estimating the accuracy of the paraxial approximation of individual Gaussian beams, or in estimating the accuracy of the asymptotic summation of paraxial Gaussian beams. The perturbation derivatives of amplitude may be useful in perturbation expansions from elastic to viscoelastic media and in estimating the accuracy of the common-ray approximations of the amplitude in the coupling ray theory.     

最优可分表示法三维波场延拓

适于复杂介质的高精度波场延拓算子是叠前深度偏移研究的重要内容。本文采用最优可分表示方法,运用正反傅立叶变换构造了三维单程波场延拓算子,算子实现了波数域变量与空间(速度)域变量分离。波数域内进行相移计算,在空间域对因介质横向变速引起的时移作修正。脉冲响应显示在区域内各速度的脉冲计算值与理论值基本一致,说明最优可分表示法叠前深度偏移可适用于强变速条件下复杂介质的成像需求。SEG／EAGE模型和实测数据的成像结果验证了本文方法对复杂构造的成像能力。     

Surface-wave polarization data and global anisotropic structure

In the past few years, seismic tomography has begun to provide detailed images of seismic velocity in the Earth's interior which, for the first time, give direct observational constraints on the mechanisms of heat and mass transfer. The study of surface waves has led to quite detailed maps of upper-mantle structure, and the current global models agree reasonably well down to wavelengths of approximately 2000 km. Usually, the models contain only elastic isotropic structure, which provides an excellent fit to the data in most cases. For example, the variance reduction for minor and major arc phase data in the frequency range 7–15 mHz is typically 65–92 per cent and the data are fit to within 1–2 standard deviations. The fit to great-circle phase data, which are not subject to bias from unknown source or instrument effects, is even better. However, there is clear evidence for seismic anisotropy in various places on the globe. This study demonstrates how much (or little) the fit to the data is improved by including anisotropy in the modelling process. It also illuminates some of the trade-offs between isotropic and anisotropic structure and gives an estimate of how much bias is introduced by neglecting anisotropy. Finally, we show that the addition of polarization data has the potential for improving recovery of anisotropic structure by diminishing the trade-offs between isotropic and anisotropic effects.     

The 3-D Truncation Filtering Methodology Defined for Planar and Spherical Models: Interpreting Gravity Data Generated by Point Masses

This paper focuses on one particular way of linear filtering the gravity data to facilitate gravity inversion or interpretation. With the use of integral transforms the gravity anomalies are transformed into new quantities that allow an easier interpretation with the help of pattern recognition. As the integral transforms are in fact filters, and as the regions of integration are caps with a variable radius, which can be systematically changed as a free parameter, we refer to such methodology as the truncation filtering. Such filters may be understood as weighted spherical windows moving over the surface, on which the gravity anomaly is defined, the kernel of the transform being the weight function. The objective of this paper is to define and deploy the truncation filtering for a planar model, i.e. for a homogenous horizontally infinite layer with embedded anomalous masses, and for a spherical model, i.e., for a homogenous massive sphere with embedded anomalous masses. Instead of the original gravity anomaly, the quantities resulting from the truncation filtering are interpreted/inverted. As we shall see, this approach has certain benefits. The fundamental concept of the truncation filtering methodology is demonstrated here in terms of the model consisting of one point mass anomaly.The relationship between the depth of the point mass and the instant of the onset of the dimple pattern observed in sequences produced by truncation filtering the synthetic gravity data generated by point masses is, for both the planar and spherical models, compiled by computer simulations, as well as derived analytically. It is shown, that the dimple pattern is a consequence of truncating the domain of the filter and is free of the choice of the kernel of the filter. It is shown, that in terms of the mean earth and depths of point masses no greater than some 100 km the spherical model may be replaced by a planar model from the perspective of the truncation filtering methodology. It is also shown, that from the viewpoint of the truncation filtering methodology the rigorous gravity anomaly may be approximated by the vertical component of the gravity disturbance. The relationship between the instant of the dimple onset and the depth of the point mass thus becomes linear and independent of the magnitude (mass) of the point mass.     

基于Biot-Squirt方程的波场模拟

Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制,对地震波和声波的传播均产生重要影响. Dvorkin和Nur提出了同时包含Biot流动和喷射流动力学机制的统一的BISQ(Biot-Squirt)模型,基于这一模型,尽管有关弹性波在多孔隙介质中的衰减和频散问题已被广泛研究,然而,基于BISQ波传播方程的波场数值模拟至今仍未见报道. 本文从同时包含两种力学机制的孔隙弹性波方程出发,利用FCT有限差分法对含流体孔隙各向同性介质中的地震波和声波进行了数值模拟,并与基于Biot流动的Biot理论之模拟结果进行比较. 数值模拟结果表明:同时包含Biot流动和喷射流动影响的地震波和声波速度比仅包含Biot流动作用的地震波和声波速度慢,慢P波的衰减比根据Biot理论模拟的慢P波衰减更强.