Cylindrical-Wave Approach for electromagnetic scattering by subsurface metallic targets in a lossy medium

An analytical solution is developed to the two-dimensional scattering problem of a plane-wave propagating in air, impinging on the interface with a dissipative soil, and interacting with a finite set of subsurface metallic targets. The Cylindrical Wave Approach is applied, the electromagnetic field scattered by the targets is expanded into cylindrical waves and use is made of the plane-wave spectrum to take into account the interaction of such waves with the planar interface between air and soil. The theoretical solution is implemented in a Fortran code. The numerical evaluation of the spectral integral relevant to reflected and transmitted cylindrical wave functions in the presence of lossy media is performed by means of Gaussian adaptive quadrature formulas. The method may return the field values in each point of the space, both in the near and far zones; moreover it may be applied for any polarization, and for arbitrary values of the cylinder sizes and positions.     

Dynamic displacements and stresses in the vicinity of a cylindrical cavity embedded in a half-space

The two-dimensional response of a viscoelastic half-space containing a buried, unlined, infinitely long cylindrical cavity of circular cross-section subjected to harmonic plane SH, P, SV and Rayleigh waves is obtained by an indirect boundary integral method based on the two-dimensional Green's functions for a viscoelastic half-space. An extensive critical review of the existing numerical results obtained by other techniques is presented together with some new numerical results describing the motion on the ground surface and the motion and stresses on the wall of the cavity for P, SV, SH and Rayleigh waves.     

On planar seismic wavefront modeling for estimating rotational ground motions: Case of 2-D SH line-source

At large hypocentral distances, it is convenient to approximate the curved transient seismic wavefronts as planar to estimate rotational ground motions from the single-station recordings of translational ground motions. In this paper, we investigate whether and when this approximation, referred to as the ‘plane-wave’ approximation, can be considered adequate close to the source. For this, we consider a simplistic source model comprising a two-dimensional, kinematic shear dislocation SH line-source buried in a homogenous, elastic half-space and assume this to be an equivalent representation of a finite-sized fault. The ‘plane-wave’ rotational motion is then synthesized from the exact translational motion solution to the assumed model and compared with the exact rotational motion solution for this model. The comparison between the two sets of rotational amplitudes in frequency domain suggests that the plane-wave approximation may be adequate, when the wavelengths of the seismic waves are much smaller than the source depth. When this is not true, the plane-wave approximation is seen to underestimate the Fourier amplitudes close to the source by several orders, particularly when the fault planes are vertically oriented. A similar comparison in the time domain indicates that a severe underestimation may also occur when the source rise time is longer than the shear-wave arrival time at the epicenter. Significant discrepancies are also observed between the waveforms of the exact and plane-wave rotational motions.     

饱和土中深埋圆柱形衬砌洞室对瞬态平面波的散射

基于Biot饱和多孔介质动力学理论,运用Laplace变换和波函数展开法,根据饱和土体与衬砌结构交界面的连续条件和衬砌结构内边界上的应力自由条件,得到饱和土中深埋圆柱形衬砌洞室对瞬态平面P波和SV波散射问题的解答,该解答可以退化成为饱和土中深埋圆柱形空穴或弹性夹塞物的情形,并很容易转换成为对稳态波散射的解.通过与已有的相关问题的解析解答进行对比,验证了该解答的正确性.同时利用Laplace逆变换的数值方法,给出了饱和土和衬砌中应力和位移场在时域内的数值解,通过算例,分析了衬砌厚度、刚度对衬砌内边界处应力集中因子的影响.     

Diffraction of obliquely incident waves by a cylindrical cavity embedded in a layered viscoelastic half-space

The three-dimensional harmonic response in the vicinity of an infinitely long, cylindrical cavity of circular cross-section buried in a layered, viscoelastic half-space is obtained when the half-space is subjected to homogeneous plane waves and surface waves impinging at an oblique angle with respect to the axis of the cavity. The solution is obtained by an indirect boundary integral method based on the use of moving Green's functions for the viscoelastic half-space. Numerical results describing the motion on the ground surface and the motion and stresses on the wall of the cavity are presented for obliquely incident P-, SV-, SH- and Rayleigh waves with different horizontal angles of incidence.     

A note on the surface motion of a semi-cylindrical canyon for incident cylindrical SH waves radiated by a finite fault

The plane-wave assumption for incident SH waves can be a good approximation for cylindrical and spherical waves radiated from finite sources, even when the source is as close as twice the size of the inhomogeneity, and when the source and the inhomogeneity are described within the same coordinate system. However, in a more general setting, and when the fault’s radiation pattern must be considered, the plane-wave approximation may not yield satisfactory answers for arbitrary orientation of the fault. Jalali et al. (2015) demonstrated this for a semi-cylindrical, sedimentary valley, and in this study we extend their results to a case in which the semi-circular, sedimentary valley is replaced by a canyon. We describe the effects of incident cylindrical waves on the amplitudes of surface motion in and near the semi-cylindrical canyon when the causative faults are at different distances and have different curvatures and orientations.     

Dynamic response of an eccentrically lined circular tunnel in poroelastic soil under seismic excitation

An exact analysis for two-dimensional dynamic interaction of monochromatic progressive plane compressional and shear seismic waves with a permeable circular tunnel lining of circumferentially varying wall thickness buried in a boundless porous elastic fluid-saturated formation is presented. The novel features of Biot dynamic theory of poroelasticity in conjunction with the translational addition theorems for cylindrical wave functions, along with the appropriate wave field expansions and the pertinent boundary conditions are employed to develop a closed-form solution in form of infinite series. The analytical results are illustrated with numerical examples in which an air-filled and water-saturated permeable tunnel lining of variable wall thickness, embedded within water-saturated surrounding formations of distinct frame properties (soft, stiff, and stiff viscoelastic soil), is insonified by fast compressional or shear waves at selected angles of incidence. The effects of liner eccentricity, interface permeability, formation material type, incident wave frequency, and angle of incidence on the hoop stress amplitude are evaluated and discussed for representative values of the parameters characterizing the system. Limiting cases are considered and good agreements with the solutions available in the literature are obtained.     

Detection of subsurface metallic utilities by means of a SAP technique: Comparing MUSIC- and SVM-based approaches

The identification of buried cables, pipes, conduits, and other cylindrical utilities is a very important task in civil engineering. In the last years, several methods have been proposed in the literature for tackling this problem. Most commonly employed approaches are based on the use of Ground Penetrating Radars, i.e., they extract the needed information about the unknown scenario starting from the electromagnetic field collected by a set of antennas. In the present paper, a statistical method, based on the use of smart antenna techniques, is used for the localization of a single buried object. In particular, two efficient algorithms for the estimation of the directions of arrival of the electromagnetic waves scattered by the targets, namely the MUltiple SIgnal Classification and the Support Vector Regression, are considered and their performances are compared.     

一种新的实时电磁逆散射方法

为解决介质圆柱体逆散射问题,提出一种新的在线逆散射方法,通过支持向量机将逆散射问题转化成一个回归估计问题. 该方法可应用于各种逆散射方面, 尤其是目标的几何与电磁参数重构和埋地目标探测. 文中首次将支持向量机方法应用到该领域,设置多个散射场的观测点,通过提取散射场的不同信息作为样本信息训练支持向量机, 建立了介质圆柱体的逆散射模型, 利用该模型重构了介质圆柱体的电磁参数,同时探测了埋地位置. 数值结果显示了该方法的有效性和准确性,为目标的实时逆散射研究提供了一种有效方法.     

Elastic wave propagation in an irregularly layered medium

The indirect boundary element method (IBEM) is used to simulate wave propagation in two-dimensional irregularly layered elastic media for internal line sources. The method is based on the integral representation for scattered elastic waves using single layer boundary sources. Fulfillment of the boundary conditions leads to a system of integral equations. Results are obtained in the frequency domain and seismograins are computed through Fourier synthesis. In order to test and validate the method we present various comparisons between our results and the time series obtained analytically for a buried line source in a half-space and by using the recently developed spectral element method (SEM).     

3-D seismic response analysis of long lined tunnels in half-space

The dynamic response of infinitely long lined tunnels with a uniform cross-section buried into an elastic or viscoelastic half-space to body and surface harmonic seismic waves is numerically determined by a special direct boundary element method in the frequency domain. The waves have an arbitrary direction of propagation with respect to the axis of the tunnel and this renders the problem three-dimensional. However, this problem is effectively reduced to a two-dimensional one by a coordinate transformation and appropriate integration of the full-space dynamic fundamental solution along the direction of the tunnel axis. Quadratic isoparametric boundary line elements and advanced numerical integration techniques for the treatment of singular integrals produce results of high accuracy. Numerical results are presented for the case of an infinitely long lined tunnel of circular cross-section and compared against those of a full three-dimensional boundary element analysis, as well as those of other methods. Thus the proposed method is illustrated and its merits demonstrated.     

圆弧形凹陷饱和土场地对平面P波散射问题的解析解

以Biot 饱和多孔介质动力学理论为基础，利用Fourier Bessel级数展开法, 得到饱和多孔介质半空间中圆弧形凹陷地形对平面P波的散射问题的解析解. 数值计算给出地表位移幅值，分析了入射波波长、入射角、圆弧高宽比对地表位移幅值的影响，并与现存的在单相介质情况下得到的结论进行对比.     

Application of the weighted residual method to diffraction by 2-D canyons of arbitrary shape: I. Incident SH waves

A solution for the two-dimensional scattering and diffraction of plane SH waves by canyons of arbitrary shape in an elastic half space is presented. The wave field for arbitrary geometry in this paper is computed numerically by the method of weighted residues (moment method). The wave displacement field computed by the present residual method for the case of a semi-circular canyon was shown to agree analytically and numerically with that computed by the exact closed form series solution. The same observations about ground amplifications, their dependence on frequencies and orientations of the incident waves, can be stated here for canyons of arbitrary shape as previously made for circular canyons.     

Diffraction of anti-plane SH waves by a semi-circular cylindrical hill with an inside concentric semi-circular tunnel

A closed-form analytic solution of two-dimensional scattering and diffraction of plane SH waves by a semicylindrical hill with a semi-cylindrical concentric tunnel inside an elastic half-space is presented using the cylindrical wave functions expansion method. The solution is reduced to solving a set of infinite linear algebraic equations. Fourier expansion theorem with the form of complex exponential function and cosine function is used. Numerical solutions are obtained by truncation of the infinite equations. The accuracy of the presented numerical results is carefully verified.     

Time-optimal path planning in dynamic flows using level set equations: theory and schemes

We develop an accurate partial differential equation-based methodology that predicts the time-optimal paths of autonomous vehicles navigating in any continuous, strong, and dynamic ocean currents, obviating the need for heuristics. The goal is to predict a sequence of steering directions so that vehicles can best utilize or avoid currents to minimize their travel time. Inspired by the level set method, we derive and demonstrate that a modified level set equation governs the time-optimal path in any continuous flow. We show that our algorithm is computationally efficient and apply it to a number of experiments. First, we validate our approach through a simple benchmark application in a Rankine vortex flow for which an analytical solution is available. Next, we apply our methodology to more complex, simulated flow fields such as unsteady double-gyre flows driven by wind stress and flows behind a circular island. These examples show that time-optimal paths for multiple vehicles can be planned even in the presence of complex flows in domains with obstacles. Finally, we present and support through illustrations several remarks that describe specific features of our methodology.     

三维偏移距平面波有限差分叠前时间偏移

本文提出了中点-半偏移距域内的三维偏移距平面波(offset plane-wave)方程,并给出了其有限差分解法.偏移距平面波可通过对CMP道集进行平面波分解(倾斜叠加或线性Radon变换)生成,然而这样做会产生严重的噪音干扰.本文提出了局部倾斜叠加方法(local slant-stacking)来消除离散线性Rado...     

圆形隧洞对地下地震动的影响:SH波入射

地下结构对地震波的散射改变了场地的动力特性,无论是地上建筑还是地下结构的安全评价和抗震设计中,目前均没有很好的考虑由于地下结构的存在对原场地地震动尤其是地下地震动的影响。基于弹性波动理论,运用波函数展开法和镜像原理,分析了弹性半空间中圆形隧洞对柱面SH波的散射问题,得到了含圆形隧洞的弹性半空间位移解析解。通过数值算例分析了圆形隧洞对原场地地震动的影响,重点考察了隧洞埋深、隧洞半径和围岩衬砌模量比等参数的影响规律。结果表明,地下结构对沿线场地的动力特性有着显著的影响,对其自身以及沿线工程结构的抗震设计提供一定参考价值。     

Scattering of plane SV waves by cylindrical canyons in saturated porous medium

On the basis of Biot dynamic theory, an analytic solution of two-dimensional scattering and diffraction of plane SV waves by circular cylindrical canyons in a half space of saturated porous media is presented in this paper for the first time. The solution is obtained by employing the Fourier–Bessel series expansion technique. Parametric studies had been carried out, which includes: the angle of incidence, the frequency of the incident SV wave, the porosity of saturated porous medium and the stiffness and Poisson's ratio of the solid-skeleton. All the outcomes are useful for the seismic analysis of the surface topography conditions.     

A strategy for the determination of the dielectric permittivity of a lossy soil exploiting GPR surface measurements and a cooperative target

In this paper we deal with an indirect measure of the dielectric permittivity of the soil starting from GPR surface data collected on a buried “cooperative” target, meant as an object buried on purpose and whose extent is known a-priori. This target is exploited in order to achieve, from its image obtained from a suitable GPR data processing, an indirect measure of the dielectric permittivity of the embedding soil. GPR data processing is based on a linear microwave tomographic approach funded on the Born Approximation. Using this Born approach on two-dimensional inversion tests, we investigate the effect of the soil's electrical conductivity and permittivity on this indirect measure and demonstrate that the electrical field scattered by a spot-like buried object permits an accurate estimation of the soil permittivity even when no information of the soil conductivity is available.     

Surface motion of a cylindrical hill of circular-arc cross-section for incident plane SH waves

A closed-form solution of two-dimensional scattering of plane SH waves by a cylindrical hill of circular-arc cross-section in a half-space is presented using the wave functions expansion method. The solution is reduced to solving a set of infinite linear algebraic equations, using the auxiliary functions and the exterior region form of the Graf's addition theorem. Numerical solutions are obtained by truncation of the infinite equations and their accuracies are demonstrated by convergence of the numerical results and by the extent to which the numerical results fit the exact boundary conditions with increasing the truncation order. The numerical results for some typical cases are then presented for checking accuracies of various numerical methods. The effects of the height-to-width ratio of the hill on surface ground motion are finally illustrated.