Reverse time migration with Gaussian beams using optimized ray tracing systems in transversely isotropic media

Prestack depth migration is a key technology for imaging complex reservoirs in media with strong lateral velocity variations. Prestack migrations are broadly separated into ray-based and wave-equation-based methods. Because of its efficiency and flexibility, ray-based Kirchhoff migration is popular in the industry. However, it has difficulties in dealing with the multi-arrivals, caustics and shadow zones. On the other hand, wave-equation-based methods produce images superior to that of the ray-based methods, but they are expensive numerically, especially methods based on two-way propagators in imaging large regions. Therefore, reverse time migration algorithms with Gaussian beams have recently been proposed to reduce the cost, as they combine the high computational efficiency of Gaussian beam migration and the high accuracy of reverse time migration. However, this method was based on the assumption that the subsurface is isotropic. As the acquired azimuth and maximum offsets increase, taking into account the influence of anisotropy on seismic migration is becoming more and more crucial. Using anisotropic ray tracing systems in terms of phase velocity, we proposed an anisotropic reverse time migration using the Gaussian beams method. We consider the influence of anisotropy on the propagation direction and calculate the amplitude of Gaussian beams with optimized correlation coefficients in dynamic ray tracing, which simplifies the calculations and improves the applicability of the proposed method. Numerical tests on anisotropic models demonstrate the efficiency and accuracy of the proposed method, which can be used to image complex structures in the presence of anisotropy in the overburden.     

星点亚像元细分定位精度分析

星点的精确定位为后续的导航解算提供数据支撑。基于高斯分布原理进行星点模拟,为星点定位提供图像数据。分别介绍了质心法和高斯曲面拟合法的亚像元细分定位原理,并研究系统误差、随机噪声、分割阈值、星等和星点坐标分布5个因素对星点定位精度的影响,分析两种方法的异同点。对模拟的星点进行定位分析发现,星点距离参考点越远定位误差越大,定位精度随噪声增大而降低。实验表明,相比高斯曲面拟合方法,质心法具有更好的实用性,其定位精度对星点的像元信息依赖性较小。     

Geostatistical analysis and conditional simulation for estimating the spatial variability of hydraulic conductivity in the Choushui River alluvial fan,Taiwan

This work evaluated the spatial variability and distribution of heterogeneous hydraulic conductivity (K) in the Choushui River alluvial fan in Taiwan, using ordinary kriging (OK) and mean and individual sequential Gaussian simulations (SGS). A baseline flow model constructed by upscaling parameters was inversely calibrated to determine the pumping and recharge rates. Simulated heads using different K realizations were then compared with historically measured heads. A global/local simulated error between simulated and measured heads was analysed to assess the different spatial variabilities of various estimated K distributions. The results of a MODFLOW simulation indicate that the OK realization had the smallest sum of absolute mean simulation errors (SAMSE) and the SGS realizations preserved the spatial variability of the measured K fields. Moreover, the SAMSE increases as the spatial variability of the K field increases. The OK realization yields small local simulation errors in the measured K field of moderate magnitude, whereas the SGS realizations have small local simulation errors in the measured K fields, with high and low values. The OK realization of K can be applied to perform a deterministic inverse calibration. The mean SGS method is suggested for constructing a K field when the application focuses on extreme values of estimated parameters and small calibration errors, such as in a simulation of contaminant transport in heterogeneous aquifers. The individual SGS realization is useful in stochastically assessing the spatial uncertainty of highly heterogeneous aquifers. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical evaluation of the dilogarithm of complex argument

A numerical algorithm to evaluate the dilogarithmic function of a complex argument is proposed. The use of the dilogarithm in celestial mechanics appears in the exact Delaunay normalization of some functions involving the equation of the centre.     

Modified precise time step integration method of structural dynamic analysis

The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.     

Discretization error for the superposition of Gaussian beams

Summary. A high-frequency asymptotic integral expansion of a time-harmonic wavefield into Gaussian beams was derived in a previous paper by Klimeš. The discretization error caused by replacing this integral superposition by a discrete summation of Gaussian beams is estimated in this paper.     

Testing the correctness of the sequential algorithm for simulating Gaussian random fields

The sequential algorithm is widely used to simulate Gaussian random fields. However, a rigorous application of this algorithm is impractical and some simplifications are required, in particular a moving neighborhood has to be defined. To examine the effect of such restriction on the quality of the realizations, a reference case is presented and several parameters are reviewed, mainly the histogram, variogram, indicator variograms, as well as the ergodic fluctuations in the first and second-order statistics. The study concludes that, even in a favorable case where the simulated domain is large with respect to the range of the model, the realizations may poorly reproduce the second-order statistics and be inconsistent with the stationarity and ergodicity assumptions. Practical tips such as the multiple-grid strategy do not overcome these impediments. Finally, extending the original algorithm by using an ordinary kriging should be avoided, unless an intrinsic random function model is sought after.     

区域集合预报系统2 m温度预报的校准技术

采用非齐次高斯回归 (NGR) 技术对国家气象中心区域集合预报系统的2 m温度预报结果开展了一阶偏差和二阶离散度的校准研究。对预报结果比较详尽的检验分析表明:校准后的2 m温度预报可靠性和预报技巧均显著提高,表现为校准后集合预报成员的均方根误差与离散度更为接近;原Talagrand直方图中的“L”形分布现象得到有效改善;Brier评分、最小连续分级概率评分 (CRPS) 明显减小,相对作用特征 (ROC) 面积增大,说明校准后的2 m温度预报表现出更好的预报技能。此外,NGR技术与自适应误差订正技术的对比试验表明,NGR在消除集合平均偏差和提高集合离散度两个方面均有优势。     

Determination of ammonia emission rates from a tunnel ventilated chicken house using passive samplers and a Gaussian dispersion model

Atmospheric deposition can provide a significant fraction of the nitrogen loading to coastal waters. The Delmarva Peninsula, on the eastern shore of the Chesapeake Bay, is a region with intense poultry production that may supply a significant source of atmospheric ammonia (NH₃). There is a need to quantify ammonia NH₃ from representative growing methods in this region in order to more accurately estimate agricultural NH₃ emissions and to develop best management practices. In this study, NH₃ emissions were determined at an 18,600-bird tunnel-ventilated chicken house using a modified sampling grid (a planar arrangement normal to the length of the house) with Ogawa passive samplers to characterize the emission plume downwind from the house. This improvement in the sampling strategy, compared to a previous study, simplified the inverse Gaussian plume analysis which improved the confidence in the emission factors. In this study, a total of four separate plume characterizations were conducted over the final 3 weeks of the 6-week broiler grow-out cycle. The mean emission factor observed at the tunnel-ventilated house, 0.13 g NH₃–N/bird-day, was an order of magnitude lower than that previously observed at a nearby side-wall ventilated house. Although not all growing variables were measured, the large difference in emission factors between the two ventilation regimes suggest that modern, tunnel-ventilated houses may result in a significant decrease in NH₃ emissions compared with traditional growing methods. This variability in emission factors underscores the need for characterizing chicken houses under various conditions and determining the factors that control these atmospheric emissions.     

Adaptive analysis of infinite beams dynamics problems using the periodic configuration update method in the time domain

This paper proposes a new approach for the assessment of the dynamic response of continuously supported infinite beams under high‐speed moving loads. A change in the representation of equations of motion in the dynamics of discrete structures is proposed to obtain an improved accuracy of the numerical integration in the time domain. The proposed numerical method called the “periodic configuration update” or “PCU method” is applied to solve the problem of a vertical moving harmonic load on an infinite classical Euler‐Bernoulli beam resting on a continuous viscoelastic foundation. This study shows the superiority of the proposed method in comparison with other methods presented in the literature that suffer from the material time derivative, i.e., convective terms, that arises from the Galilean transformation. To confront this numerical problem, the PCU method retains the principle of the spatial follow of loads while zeroing the relative velocity with the traversed beam via a step‐by‐step adaptive integration of the equation of structural dynamics. The dynamic load is modeled with high theoretical velocities that can reach the critical velocity of the studied beam with different angular frequencies belonging to moderate frequency range. A parametric study is carried out to analyze the influence of key parameters on the convergence. The obtained results show a high efficiency of the PCU method for solving these types of problems relative to the dynamics of high speed trains/tracks.