Numerical modelling of stratified flow: a spectral approach

A numerical formulation is developed to solve the three-dimensional hydrodynamic equations which describe flow in a stratified sea.Arbitrary continuous physically realistic variations of density and eddy viscosity can be included in the model, which is sufficiently flexible to be applicable to sea areas of any horizontal extent and depth. A continuous current profile from sea surface to sea bed, is computed with the model. A method for expanding computed current profiles in terms of vertical modes is proposed and the contribution of these modes to the current profiles is considered.The time variation of the wind-induced circulation of a stratified lake in response to a suddenly imposed and maintained wind stress is examined. Calculations show that the wind-driven surface current is modulated by the internal seiche motion of the lake.     

Numerical modelling of spectral sky radiance during a total solar eclipse

Ground-based observations of the spectral sky radiance within the totality region during a total solar eclipse are of significant interest because the contribution from direct and single scattered light from the solar disk is eliminated. In the present paper, we develop a numerical model of the spectral sky radiance during totality, evaluate the contribution of double scattered sunlight to the sky radiance at totality and compare it to solar corona emissions. The results show that the single scattered coronal light is the major contaminant in ground-based observations of the corona, opposite to what was believed before. For observation directions outside the corona the double scattered sun radiation dominates over the single scattered coronal light.     

Seismic wave modeling in viscoelastic VTI media using spectral element method

Spectral element method (SEM) for elastic media is well known for its great flexibility and high accuracy in solving problems with complex geometries. It is an advanced choice for wave simulations. Due to anelasticity of earth media, SEM for elastic media is no longer appropriate. On fundamental of the second-order elastic SEM, this work takes the viscoelastic wave equations and the vertical transversely isotropic (VTI) media into consideration, and establishes the second-order SEM for wave modeling in viscoelastic VTI media. The second-order perfectly matched layer for viscoelastic VTI media is also introduced. The problem of handling the overlapped absorbed corners is solved. A comparison with the analytical solution in a two-dimensional viscoelastic homogeneous medium shows that the method is accurate in the wave-field modeling. Furtherly, numerical validation also presents its great flexibility in solving wave propagation problems in complex heterogeneous media. This second-order SEM with perfectly matched layer for viscoelastic VTI media can be easily applied in wave modeling in a limited region.     

方位各向异性粘弹性介质波场数值模拟

当地震信号通过复杂地球介质时,地层除了表现为各向异性,还表现为内在的粘弹性特征.因此,为准确描述地震波在地球介质中的传播特征,理想的地球介质模型应该能够模拟岩石的各向异性特征和衰减特征.本文给出了各向异性粘弹性介质模型的波动方程及其差分格式,并利用有限差分法实现了地震波波场数值模拟.结果表明了该介质模型中地震波场特征与各向异性主轴方位和介质的粘滞性参数之间的关系.     

任意各向异性介质三维非结构谱元法直流电阻率正演模拟研究

各向异性是地电异常解释中不可忽视的因素,广泛存在于裂隙或层理发育的地质环境中.本文针对任意各向异性条件下直流电阻率法三维正演问题进行研究,结合非结构谱元法建立模拟算法,充分利用谱方法的指数收敛性以及非结构有限元对地形和复杂异常体刻画能力,提高计算精度和效率.通过灵活的四面体网格剖分和高阶谱插值,实现了复杂介质任意各向异性模型电阻率响应的高精度数值模拟.我们首先通过层状各向异性模型验证本文非结构谱元法的计算精度,进而我们以半空间中立方体模型为例分析各向异性对电阻率响应的影响特征,并通过计算针对不同各向异性参数的视电阻率极性图,探究地下介质各向异性特征识别方法.最后,我们针对典型的山脊模型计算和分析存在地形效应条件下各向异性电流场分布及视电阻率特征.模型计算结果表明基于四面体网格的谱元法模拟带有复杂地形和异常体的任意各向异性模型具有很高的计算精度.本文的研究成果将在推进电阻率方法用于解决裂隙及层理等环境和工程地质问题中发挥积极作用.

     

Acoustic viscoelastic modeling by frequency-domain boundary element method

Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element methods. Targeted at viscoelastic numerical modeling for multilayered media, the constant-Q acoustic wave equation is transformed into the corresponding wave integral representation with its Green’s function accounting for viscoelastic coefficients. An efficient alternative for full-waveform solution to the integral equation is proposed in this article by extending conventional frequency-domain boundary element methods to viscoelastic media. The viscoelastic boundary element method enjoys a distinct characteristic of the explicit use of boundary continuity conditions of displacement and traction, leading to a semi-analytical solution with sufficient accuracy for simulating the viscoelastic effect across irregular interfaces. Numerical experiments to study the viscoelastic absorption of different Q values demonstrate the accuracy and applicability of the method.     

Numerical modelling of time-harmonic seismic wave fields in simple structures by the gaussian beam method. part I

m¶rt;au n unm ¶rt; u ¶rt;uau uu n m nm ¶rt; ¶rt;. a uau ¶rt;m a ma m nau a¶rt;a, amu aa mu. au mu uu naam, m m¶rt;au n ¶rt;am ¶rt;mam m mam ¶rt;a u amu aa mu (. u. 9). aumu amu, uauau n ma nam ¶rt;a . ¶rt; amu m am ¶rt;m ¶rt;am mumm mam a naama am, an. aa uuau n, naama uauau n, um¶rt;.     

Numerical modelling of time-harmonic seismic wave fields in simple structures by the Gaussian beam method. Part II

ma ama m n¶rt;u am[1]. u ama, m¶rt;au n unm ¶rt; uuauu ma m nau a¶rt;a. a uau ¶rt;m ma a muna . au mu mamau n [1], m m¶rt;au n ¶rt;am ¶rt;mam m mam ¶rt;a u amu aa mu. aumu amu, uauau n ma nuu ¶rt;a . m am ¶rt;am mumm mam a naama am, aa uu L ₀ au n. aa, m am mam namuu aum m L ₀ . amu aa mu u aumu amu, a m mam n ¶rt; u L ₀ . u uu L ₀ , anum¶rt; ma amu aa mu nuam, nu aua nam m, m uu L ₀ auum. aumu amu, u L ₀ um m uau m a uuu anum¶rt; u, mmmu a. uu L ₀ , uuu nam. ma uam mumm mam ¶rt;a ¶rt;u naama am. uau ¶rt;m ma anum¶rt; u ma muna S, S u SS.     

2D finite-difference viscoelastic wave modelling including surface topography

We have pursued two-dimensional (2D) finite-difference (FD) modelling of seismic scattering from free-surface topography. Exact free-surface boundary conditions for the particle velocities have been derived for arbitrary 2D topographies. The boundary conditions are combined with a velocity–stress formulation of the full viscoelastic wave equations. A curved grid represents the physical medium and its upper boundary represents the free-surface topography. The wave equations are numerically discretized by an eighth-order FD method on a staggered grid in space, and a leap-frog technique and the Crank–Nicholson method in time.
In order to demonstrate the capabilities of the surface topography modelling technique, we simulate incident point sources with a sinusoidal topography in seismic media of increasing complexities. We present results using parameters typical of exploration surveys with topography and heterogeneous media. Topography on homogeneous media is shown to generate significant scattering. We show additional effects of layering in the medium, with and without randomization, using a von Kármán realization of apparent anisotropy. Synthetic snapshots and seismograms indicate that prominent surface topography can cause back-scattering, wave conversions and complex wave patterns which are usually discussed in terms of inter-crust heterogeneities.     

球坐标系下谱元法三维地震波场模拟

谱元法已成为区域性乃至大陆性尺度地震波场模拟的重要工具.对于区域或大陆尺度层析成像而言,地球曲率不可忽略,此时模拟地震波传播采用球坐标系更为合适.本文从球坐标系下弹性波动方程弱形式出发,基于球坐标系变分原理给出了球坐标系下求解三维地震波方程的谱元法.另一方面,计算Fréchet敏感核是进行全波形反演的关键,本文借助伴随原理,推导了全波走时层析成像三维Fréchet敏感核表达式.为了验证球坐标系下谱元法的精度,我们将数值模拟结果与normal mode方法得到的解析解在1-D PREM模型下进行了对比.同时,我们将此方法应用到华北克拉通区域,以期获得地球内部结构精确成像.基于3-D全球径向各向异性地幔模型S362ANI和3-D地壳模型Crust1.0,我们建立了华北克拉通初始3-D背景模型,并将数值模拟结果与实际观测台站记录波形资料进行对比分析,利用互相关方法提取走时残差,最后给出了Fréchet敏感核在3-D空间中的分布,这些工作为下一步开展球坐标系下三维大尺度全波形反演奠定了基础.

     

Numerical modelling of suspended sediments

The turbulent advection-diffusion mathematical model in three-dimensional space is solved by a mixed finite element finite difference method. Linear finite elements in the vertical direction and central finite differences in the horizontal directions are used coupled with the Galerkin error minimization procedure. The integration in time is performed in fractional steps (one explicit one implicit) by splitting the differential operator. The method is illustrated by application to the three-dimensional movement of suspended sediment. Its accuracy is checked by comparison to analytical solutions and its efficiency is gauged relative to finite elements and implicit finite difference solutions for two-dimensional suspended sediment transport over a dredged channel.     

GPR modelling by the Fourier method: improvement of the algorithm

We improve two aspects of the modelling scheme for the simulation of electromagnetic radio waves, based on the Fourier pseudospectral method.
When there are large contrasts in the material properties, use of the standard algorithm (regular grid) causes a series of artefacts, as, for instance, ringing and acausal events. These problems, due to the non-locality of the differential operator, are solved by using the staggered Fourier method (staggered grid).
Realistic radiation patterns can be obtained from simple combinations of magnetic and electric sources. If the directivity pattern of the antenna is known, from either a finite-difference simulation or an analytic evaluation or an experimental characterization, it can then be simulated by a composite-source concept. This effective source is implemented in the modelling algorithm by means of a perturbation technique, which first computes the intensity and directional spectra of the single electromagnetic sources. Their location is optimized to obtain the best fit with a minimum number of sources. The approach is, in principle, valid for the far-field radiation pattern of the antenna.     

黏弹性双相介质参数反演的同伦方法

本文以介质的动态时域响应（位移、速度、加速度）作为反演的依据,对黏弹性双相介质的材料参数进行了反分析研究.根据附加条件,并考虑到由于测量误差和测量数据的不充分造成的反演问题的不适定性,将黏弹性双相介质的参数反演问题转化为非线性泛函的极小值问题,然后应用大范围收敛的同伦方法求出使泛函极小的根作为反问题的解.最后,以二维半空间黏弹性双相介质模型为例进行了数值反演分析,数值结果表明本文方法在参数的反演过程中目标函数稳定收敛,且具有一定的抗“噪声”干扰性能;与利用速度和加速度信息反演的结果相比,利用位移信息反演的结果具有较高的精度和抗噪声能力.     

Attenuation from microseismic datasets by the peak frequency method benchmarked with the spectral ratio method

Recently proposed peak-frequency method is used to estimate the P- and S-wave quality factors from microseismic events. We use a downhole monitoring dataset of 10 high signal-to-noise ratio microseismic events to calculate P- and S-wave effective attenuation of a carbonate reservoir. We benchmark these results with the spectral ratio method and obtain mutually consistent results. Additionally we develop and test two techniques of peak frequency determination. We show that the peak frequency method can be successfully used in the estimation of the quality factor and it provides precise measurements of attenuation.     

Numerical Modeling of Transient Seismic Fields in Viscoelastic Media Based on the Laguerre Spectral Method

--The wave propagation in real media can be described within the theory of linear viscoelasticity. The presence of convolutional integral in Boltzmann's superposition principle poses the main difficulties in implementing the direct numerical methods in time domain. The paper presents a new algorithm, based on the application of the spectral Laguerre method for the approximation of temporal derivatives and convolution as applied to the problem of seismic wave propagation in the heterogeneous viscoelastic medium. Examples of the calculation of synthetic seismograms for different models of viscoelastic media are presented.     

Numerical modelling of drop load tests

Assessment of the attenuation of induced vibrations in the ground plays an important role in evaluating comfort and structural safety. Analytical and empirical wave attenuation relationships of increasing complexity and detail are presented in the paper, as well as a numerical model that accurately reproduces wave attenuation for a well-documented site, namely the one of the Tower of Pisa, Italy. A new source model is calibrated on near-field data and used as input for the dynamic coupled consolidation finite element analysis to achieve a satisfactory simulation. The accuracy of simpler analytical and empirical approaches is then comprehensively assessed through comparison with the validated numerical model and the field data obtained from geophones at various distances from the impact source.     

Numerical modelling of downstream migrating antidunes

A numerical model is presented that compute the geometrical dimensions and movement of downstream migrating antidunes. The model solves the Navier–Stokes equations together with the k‐epsilon turbulence model to find the water flow field over the bedforms. A two‐dimensional width‐averaged grid is used. The bed elevation changes are computed by solving the convection–diffusion equation for suspended sediments and bedload, together with the Engelund–Hansen sediment transport formula. The free surface is computed with an algorithm based on water continuity in the surface cells. Non‐orthogonal adaptive grids were used, moving vertically with the computed location of the bed and the free water surface. The numerical model was tested on data from a physical model study where regular downstream migrating antidunes had been observed. The numerical model started out with a flat bed and the trains of antidunes formed over time. Many of the physical processes observed in earlier studies were replicated by the numerical model. Four dune parameters were computed in the current tests: The antidune wavelength, height and celerity, together with the average water depth. The antidune wavelengths were best predicted with an accuracy of 3 to 8% compared with the measurements. The antidune heights were computed with a deviation of 11 to 25% compared with an empirical formula. The water depths over the antidunes were predicted with an accuracy of 3 to 9% related to the measured values. The average antidune celerity was the parameter with largest deviation: For the coarsest grid it was overpredicted with 37%. Copyright © 2017 John Wiley & Sons, Ltd.     

Numerical modelling of floating debris in the world's oceans

A global ocean circulation model is coupled to a Lagrangian particle tracking model to simulate 30 years of input, transport and accumulation of floating debris in the world ocean. Using both terrestrial and maritime inputs, the modelling results clearly show the formation of five accumulation zones in the subtropical latitudes of the major ocean basins. The relative size and concentration of each clearly illustrate the dominance of the accumulation zones in the northern hemisphere, while smaller seas surrounded by densely populated areas are also shown to have a high concentration of floating debris. We also determine the relative contribution of different source regions to the total amount of material in a particular accumulation zone. This study provides a framework for describing the transport, distribution and accumulation of floating marine debris and can be continuously updated and adapted to assess scenarios reflecting changes in the production and disposal of plastic worldwide.     

GPR spectral analysis for clay content evaluation by the frequency shift method

A recent approach relates the shift of the frequency peak of the Ground Penetrating Radar (GPR) spectrum with the increasing of the moisture content in the soil. The weakness characterizing this approach is represented by the needs of high resolution signals, whereas GPR spectra are affected by low resolution. The novelty introduced by this work is twofold. First, we evidence that clay content information is present in the location where the maximum amplitude of the GPR spectra occurs. Then, we propose three super resolution methods, namely parabolic, triangular, and sinc-based interpolators, to further refine the location of the frequency peak. In fact, it is really important to be able to find this location quite precisely, to obtain accurate estimates of clay content. We show that the peak location can be found best through sinc-interpolation in the frequency domain of the measured data. Our experimental results confirm the effectiveness of the proposed approach to resolve a frequency shift in the GPR spectrum, even for a small amount of clay.