A high-resolution dispersion imaging method of seismic surface waves based on chirplet transform

The surface-wave analysis method is widely adopted to build a near-surface shear-wave velocity structure. Reliable dispersion imaging results form the basis for subsequent picking and inversion of dispersion curves. In this paper, we present a high-resolution dispersion imaging method (CSFK) of seismic surface waves based on chirplet transform (CT). CT introduces the concept of chirp rate, which could focus surface-wave dispersion energy well in time-frequency domain. First, each seismic trace in time-distance domain is transformed to time-frequency domain by CT. Thus, for each common frequency gather, we obtain a series of 2D complex-valued functions of time and distance, which are called pseudo-seismograms. Then, we scan a series of group velocities to obtain the slanting-phase function and perform a spatial Fourier transform on the slanting-phase function to get its amplitude. In addition, power operation is adopted to increase the amplitude difference between dispersion energy and noise. Finally, we generate the dispersion image by searching for the maximum amplitude of a slanting-phase function. Because the CSFK method considers the position of surface-wave energy in the time-frequency domain, this largely eliminates the noise interference from other time locations and improves the resolution and signal-to-noise ratio of the dispersion image. The results of synthetic test and field dataset processing demonstrate the effectiveness of the proposed method. In addition, we invert all 120 sets of dispersion curves extracted from reflected wave seismic data acquired for petroleum prospecting. The one-dimensional inversion shear-wave velocity models are interpolated into a two-dimensional profile of shear-wave velocity, which is in good agreement with the borehole data.

A Novel Fast Multiple-Scattering Approximate Model for Oceanographic Lidar

An effective lidar simulator is vital for its system design and processing algorithms. However, laser transmission is a complex process due to the effects of sea surface and various interactions in seawater such as absorption, scattering, and so on. It is sophisticated and difficult for multiple scattering to accurately simulate. In this study, a multiple-scattering lidar model based on multiple-forward-scattering-single-backscattering approximation for oceanic lidar was proposed. Compared with previous analytic models, this model can work without assuming a homogeneous water and fixed scattering phase function. Besides, it takes consideration of lidar system and environmental parameters including receiver field of view, different scattering phase functions, particulate sizes, stratified water, and rough sea surface. One should note that because the scattering phase function is difficult to determine accurately, the simulation accuracy may be reduced in a complex oceanic environment. The Cox–Munk model used in our method simulates capillarity waves but ignores gravity waves, and the pulse stretching is not included. The wide-angle scattering occurs in the dense subsurface phytoplankton, which sometimes makes it hard to use this model. In this study, we firstly derived this method based on an analytical solution by convolving Gaussians of the forward-scattering contribution of layer dr and the energy density at R in the small-angle-scattering approximation. Then, the effects of multiple scattering and water optical properties were analyzed using the model. Meanwhile, the validation with Monte Carlo model was implemented. Their coefficient of determination is beyond 0.9, the RMSE is within 0.02, the MAD is within 0.02, and the MAPD is within 8%, which indicates that our model is efficient for oceanographic lidar simulation. Finally, we studied the effects of FOV, SPF, rough sea surface, stratified water, and particle size. These results can provide reference for the design of the oceanic lidar system and contribute to the processing of lidar echo signals.

Phase Shift Analysis for Alpha-alpha Elastic Scattering using Phase Function Method for Gaussian Local Potential

The phase shifts for α- α scattering have been modeled using a two parameter Gaussian local potential. The time independent Schrodinger equation (TISE) has been solved iteratively using Monte-Carlo approach till the S and D bound states of the numerical solution match with the experimental binding energy data in a variational sense. The obtained potential with best fit parameters is taken as input for determining the phase-shifts for the S channel using the non-linear first order differential equation of the phase function method (PFM). It is numerically solved using 5th order Runge-Kutta (RK-5) technique. To determine the phase shifts for the ℓ=2 and 4 scattering state i.e. D and G-channel, the inversion potential parameters have been determined using variational Monte-Carlo (VMC) approach to minimize the realtive mean square error w.r.t. the experimental data.

Neutron-Proton Scattering Phase Shifts in S-Channel using Phase Function Method for Various Two Term Potentials

The scattering phase shifts for n-p scattering have been modeled using various two term exponential type potentials such as Malfliet-Tjon, Manning-Rosen and Morse to study the phase shifts in the S-channels. As a first step, the model arameters for each of the potentials are determined by obtaining binding energy of the deuteron using matrix methods vis-a-vis Variational Monte-Carlo (VMC) technique to minimize the percentage error w.r.t. the experimental value. Then, the first order ODE as given by phase function method (PFM), is numerically solved using 5th order Runge-Kutta (RK-5) technique, by substituting the obtained potentials for calculating phase shifts for the bound 3S1 channel. Finally, the potential parameters are varied in least squares sense using VMC technique to obtain the scattering phase-shifts for each of the potentials in the 1S0 channel. The numerically obtained values are seen to be matching with those obtained using other analytical techniques and a comparative analysis with the experimental values up to 300 MeV is presented.