Seismic inversion with generalized Radon transform based on local second-order approximation of scattered field in acoustic media

Sound velocity inversion problem based on scattering theory is formulated in terms of a nonlinear integral equation associated with scattered field. Because of its nonlinearity, in practice, linearization algorisms (Born/single scattering approximation) are widely used to obtain an approximate inversion solution. However, the linearized strategy is not congruent with seismic wave propagation mechanics in strong perturbation (heterogeneous) medium. In order to partially dispense with the weak perturbation assumption of the Born approximation, we present a new approach from the following two steps: firstly, to handle the forward scattering by taking into account the second-order Born approximation, which is related to generalized Radon transform (GRT) about quadratic scattering potential; then to derive a nonlinear quadratic inversion formula by resorting to inverse GRT. In our formulation, there is a significant quadratic term regarding scattering potential, and it can provide an amplitude correction for inversion results beyond standard linear inversion. The numerical experiments demonstrate that the linear single scattering inversion is only good in amplitude for relative velocity perturbation ( \( \delta_{c}/c_{0} \) ) of background media up to 10 %, and its inversion errors are unacceptable for the perturbation beyond 10 %. In contrast, the quadratic inversion can give more accurate amplitude-preserved recovery for the perturbation up to 40 %. Our inversion scheme is able to manage double scattering effects by estimating a transmission factor from an integral over a small area, and therefore, only a small portion of computational time is added to the original linear migration/inversion process.     

A Compact Explicit Difference Scheme of High Accuracy for Extended Boussinesq Equations

Presented here is a compact explicit difference scheme of high accuracy for solving the extended Boussinesq equations.For time discretization,a three-stage explicit Runge-Kutta method with TVD property is used at predicting stage,a cubic spline function is adopted at correcting stage,which made the time discretization accuracy up to fourth order;For spatial discretization,a three-point explicit compact difference scheme with arbitrary order accuracy is employed.The extended Boussinesq equations derived by Beji and Nadaoka are solved by the proposed scheme.The numerical results agree well with the experimental data.At the same time,the comparisons of the two numerical results between the present scheme and low accuracy difference method are made,which further show the necessity of using high accuracy scheme to solve the extended Boussinesq equations.As a valid sample,the wave propagation on the rectangular step is formulated by the present scheme,the modelled results are in better agreement with the experimental data than those of Kittitanasuan.     

单程波最优化可分近似偏移方法精度分析(英文)

单程波算子的可分近似形式是关于空间变量和波数变量的函数,这使得快速傅里叶变换的利用成为可能,从而极大提高了计算效率。从函数近似的角度来看,最优可分近似方法与其他几种可分近似方法一样都具有单程波算子可分表达的形式,但最优可分近似方法是其中唯一的函数整体近似的方法,这就造成了该方法在相位误差曲线、脉冲响应、模型偏移结果上表现出与其他可分近似方法不一样的特征:该方法精度较高,且随着阶数的增大,对速度变化的敏感度降低。     

一种模拟非均匀介质中弹性波传播新的k-space方法

传统的伪谱（PS）方法,采用傅里叶变换（FT）计算空间导数具有很高的精度,每个波长仅需要两个采样点,而时间导数采用有限差分（FD）近似因而精度较低.当采用大时间步长时,由于时空精度不平衡,PS法存在不稳定性问题.原始的k-space方法可以有效地克服这些问题但是却无法适用于非均匀介质.为了提高原始k-space方法模拟非均匀介质波动方程的精度,我们提出了一种新的k-space算子族.它是用非均匀介质的变速度代替原k-space算子中的常数补偿速度构造得到,引入低秩近似可以高效求解.我们将构造的新的k-space算子应用于耦合的二阶位移波动方程,而不是交错网格一阶速度应力波动方程,使模拟弹性波的计算存储量减少.我们从数学上证明了基于二阶波动方程的k-space方法与基于一阶波动方程的k-space方法是等价的.数值模拟实验表明,与传统的PS、交错网格PS和原始的k-space方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.

     

TTI介质qP波逆时偏移中伪横波噪声压制方法

在对地下复杂构造介质,特别是盐丘侧翼及岩下区域进行成像时,相对于传统的各向同性逆时偏移和VTI逆时偏移,具有倾斜对称轴的TTI逆时偏移成像效果最优.不仅反射同相轴更加的连续,而且能量得到了更好的聚焦.传统的各向异性介质全弹性波RTM的计算量大且计算效率低.由于目前仍以纵波勘探为主,因此TTI逆时偏移qP波波动方程的选取显得尤为重要.为了提高计算效率,采用将沿着对称轴方向的横波速度设为零的方法,简化得到qP波波动方程.然而,这样会引入一种严重影响成像效果的低速度、低振幅的qSV波人为干扰.本文建立了qP波方程的完全匹配层控制方程,而后借助于辅助波场采用一种高效的压制伪横波噪声传播的方法,通过模型测试验证了该方法的有效性.     

基于有效邻域波场近似的起伏地表保幅高斯束偏移

随着我国陆上地震勘探向复杂地表探区的转移,高精度、适应性强的地震成像方法在地震资料的处理、解释及后续属性分析、储层预测中具有重要意义.本文基于有效邻域波场近似理论发展了一种成像精度更高且适用于复杂起伏地表条件的叠前保幅高斯束偏移方法.在传统水平地表高斯束偏移的基础上,本文根据中心射线附近有效邻域内高斯束表征的近似波场,导出了起伏地表条件下具有相对振幅保持的高斯束偏移公式,并给出了一种精度更高的旁轴射线传播角度计算方法.同现有的高斯束偏移方法相比,本文方法不仅考虑了起伏地表对高斯束走时的线性影响,而且首次引入了由地表高程差异和近地表速度变化引起的二次时差校正项和振幅校正项,使得成像结果更加准确可靠.两个典型模型算例验证了本文方法的正确性和有效性.     

Ekman边界层动力学的理论研究

大气边界层及其与自由大气之间的相互作用具有明显的非线性特征,而这些特征是经典Ekman理论所不能描述的,因此,发展中等复杂程度(介于完全模式与经典Ekman模型之间)的大气边界层动力学模式,简称中间模式,对人们从理论上认识大气边界层动力学过程的非线性特征具有重要意义。本文对目前最具代表性的几个中间边界层模型:地转动量近似边界层模型、Ekman动量近似边界层模型以及弱非线性边界层模型进行了总结和分析,阐述了Ekman层主要动力学特征。通过分析上述各模型的理论框架,揭示了各模型的物理意义及其在描述Ekman边界层基本动力特征上的优点和局限性,并指出尽管在细节定量描述上有差异,但各中间模型对Ekman层动力学特征的定性描述具有很好的一致性。对于这些Ekman边界层近似理论模型的进一步应用问题,主要回顾和总结了利用上述模型探讨地形边界层结构、大气锋生过程、低层锋面结构和环流以及边界层日变化、低空急流形成等动力学问题的研究,并对这些研究所揭示的Ekman层动力学特征及其对自由大气低层运动的影响进行了分析,结果表明,这些Ekman边界层近似模型可以较好地揭示大气边界层动力学特征,在大气边界层动力学及其与自由大气相互作用的研究上具有重要价值。另外,还对目前Ekman边界层理论研究中存在的问题进行了一些分析,提出了有待进一步研究的科学问题。     

A two‐dimensional model of hydraulic performance of stormwater infiltration systems

Stormwater infiltration systems are a popular method for urban stormwater control. They are often designed using an assumption of one‐dimensional saturated outflow, although this is not very accurate for many typical designs where two‐dimensional (2D) flows into unsaturated soils occur. Available 2D variably saturated flow models are not commonly used for design because of their complexity and difficulties with the required boundary conditions. A purpose‐built stormwater infiltration system model was thus developed for the simulation of 2D flow from a porous storage. The model combines a soil moisture–based model for unsaturated soils with a ponded storage model and uses a wetting front‐tracking approach for saturated flows. The model represents the main physical processes while minimizing input data requirements. The model was calibrated and validated using data from laboratory 2D stormwater infiltration trench experiments. Calibrations were undertaken using five different combinations of calibration data to examine calibration data requirements. It was found that storage water levels could be satisfactorily predicted using parameters calibrated with either data from laboratory soils tests or observed water level data, whereas the prediction of soil moistures was improved through the addition of observed soil moisture data to the calibration data set. Copyright © 2012 John Wiley & Sons, Ltd.     

The sensitivity of simulated flow and water quality response to spatial heterogeneity on a hillslope in the Tarrawarra catchment,Australia

Models are widely used to simulate hydrological response and the generation and transport of constituents such as salt, phosphorus, and nitrogen from catchments to streams. Several models use a spatial representation with catchments divided into subcatchments. Variations in land use and other characteristics within subcatchments are represented by spatially lumped hydrologic response units (HRUs) or functional units instead of using fully distributed models. This approach disregards any spatial interaction between HRUs, including their connectivity to each other and to the stream and the influence of these interactions on water and constituent export. A spatially explicit hydrological model (Thales) was used to simulate a variety of theoretical catchments with soils dominated by combinations of infiltration excess, saturation excess, and subsurface stormflow processes and different soil constituent concentrations that were spatially interacting (i.e. located along a hillslope sequence). The modelling results show that the response of both runoff and concentration is sensitive to varying spatial arrangements due to interactions of runoff, infiltration, and chemical processes between the different soil types in many but not all situations. Results highlight the importance of considering connectivity of pathways when modelling hydrological response and constituents export. This is achieved by comparing pairs of simulations and the corresponding differences in the exported loads. Copyright © 2009 John Wiley & Sons, Ltd.