The reverse time migration(RTM)of ground penetrating radar(GPR)is usually implemented in its two-dimensional(2D)form,due to huge computational cost.However,2D RTM algorithm is difficult to focus the scattering signal ...The reverse time migration(RTM)of ground penetrating radar(GPR)is usually implemented in its two-dimensional(2D)form,due to huge computational cost.However,2D RTM algorithm is difficult to focus the scattering signal and produce a high precision subsurface image when the object is buried in a complicated subsurface environment.To better handle the multi-off set GPR data,we propose a three-dimensional(3D)prestack RTM algorithm.The high-order fi nite diff erence time domian(FDTD)method,with the accuracy of eighth-order in space and second-order in time,is applied to simulate the forward and backward extrapolation electromagnetic fi elds.In addition,we use the normalized correlation imaging condition to obtain pre-stack RTM result and the Laplace fi lter to suppress the low frequency noise generated during the correlation process.The numerical test of 3D simulated GPR data demonstrated that 3D RTM image shows excellent coincidence with the true model.Compared with 2D RTM image,the 3D RTM image can more clearly and accurately refl ect the 3D spatial distribution of the target,and the resolution of the imaging results is far better.Furthermore,the application of observed GPR data further validates the eff ectiveness of the proposed 3D GPR RTM algorithm,and its fi nal image can more reliably guide the subsequent interpretation.展开更多
Supercavitating flow around a slender symmetric wedge moving at variable velocity in static fluid has been studied. Singular integral equation for the flow has been founded through distributing the sources and sinks o...Supercavitating flow around a slender symmetric wedge moving at variable velocity in static fluid has been studied. Singular integral equation for the flow has been founded through distributing the sources and sinks on the symmetrical axis. The supereavity length at each moment is determined by solving the singular integral equation with finite difference method. The supercavity shape at each moment is obtained by solving the partial differential equation with variable coefficient. For the case that the wedge takes the impulse and uniformly variable motion, numerical results of time history of the supercavity length and shape are presented. The calculated results indicate that the shape and the length of the supercavity vary in a similar way to the case that the wedge takes variable motion, and there is a time lag in unsteady supercavitating flow induced by the variation of wedge velocity.展开更多
基金This work is supported by the National Natural Science Foundation of China(No.41604039,41604102,41764005,41574078)Guangxi Natural Science Foundation project(No.2020GXNSFAA159121,2016GXNSFBA380215).
文摘The reverse time migration(RTM)of ground penetrating radar(GPR)is usually implemented in its two-dimensional(2D)form,due to huge computational cost.However,2D RTM algorithm is difficult to focus the scattering signal and produce a high precision subsurface image when the object is buried in a complicated subsurface environment.To better handle the multi-off set GPR data,we propose a three-dimensional(3D)prestack RTM algorithm.The high-order fi nite diff erence time domian(FDTD)method,with the accuracy of eighth-order in space and second-order in time,is applied to simulate the forward and backward extrapolation electromagnetic fi elds.In addition,we use the normalized correlation imaging condition to obtain pre-stack RTM result and the Laplace fi lter to suppress the low frequency noise generated during the correlation process.The numerical test of 3D simulated GPR data demonstrated that 3D RTM image shows excellent coincidence with the true model.Compared with 2D RTM image,the 3D RTM image can more clearly and accurately refl ect the 3D spatial distribution of the target,and the resolution of the imaging results is far better.Furthermore,the application of observed GPR data further validates the eff ectiveness of the proposed 3D GPR RTM algorithm,and its fi nal image can more reliably guide the subsequent interpretation.
基金Sponsored by the National Natural Science Foundation of China(Grant No.10832007)
文摘Supercavitating flow around a slender symmetric wedge moving at variable velocity in static fluid has been studied. Singular integral equation for the flow has been founded through distributing the sources and sinks on the symmetrical axis. The supereavity length at each moment is determined by solving the singular integral equation with finite difference method. The supercavity shape at each moment is obtained by solving the partial differential equation with variable coefficient. For the case that the wedge takes the impulse and uniformly variable motion, numerical results of time history of the supercavity length and shape are presented. The calculated results indicate that the shape and the length of the supercavity vary in a similar way to the case that the wedge takes variable motion, and there is a time lag in unsteady supercavitating flow induced by the variation of wedge velocity.
