Prestack nonstationary deconvolution based on variable-step sampling in the radial trace domain

The conventional nonstationary convolutional model assumes that the seismic signal is recorded at normal incidence. Raw shot gathers are far from this assumption because of the effects of offsets. Because of such problems, we propose a novel prestack nonstationary deconvolution approach. We introduce the radial trace （RT） transform to the nonstationary deconvolution, we estimate the nonstationary deconvolution factor with hyperbolic smoothing based on variable-step sampling （VSS） in the RT domain, and we obtain the high-resolution prestack nonstationary deconvolution data. The RT transform maps the shot record from the offset and traveltime coordinates to those of apparent velocity and traveltime. The ray paths of the traces in the RT better satisfy the assumptions of the convolutional model. The proposed method combines the advantages of stationary deconvolution and inverse Q filtering, without prior information for Q. The nonstationary deconvolution in the RT domain is more suitable than that in the space-time （XT） domain for prestack data because it is the generalized extension of normal incidence. Tests with synthetic and real data demonstrate that the proposed method is more effective in compensating for large-offset and deep data.

DISCRETE ENERGY ANALYSIS OF THE THIRD-ORDER VARIABLE-STEP BDF TIME-STEPPING FOR DIFFUSION EQUATIONS

This is one of our series works on discrete energy analysis of the variable-step BDF schemes.In this part,we present stability and convergence analysis of the third-order BDF(BDF3)schemes with variable steps for linear diffusion equations,see,e.g.,[SIAM J.Numer.Anal.,58:2294-2314]and[Math.Comp.,90:1207-1226]for our previous works on the BDF2 scheme.To this aim,we first build up a discrete gradient structure of the variable-step BDF3 formula under the condition that the adjacent step ratios are less than 1.4877,by which we can establish a discrete energy dissipation law.Mesh-robust stability and convergence analysis in the L^(2) norm are then obtained.Here the mesh robustness means that the solution errors are well controlled by the maximum time-step size but independent of the adjacent time-step ratios.We also present numerical tests to support our theoretical results.

Variable-step integration method for milling chatter stability prediction with multiple delays

This paper analyzes the stability of milling with variable pitch cutter and tool runout cases characterized by multiple delays,and proposes a new variable-step numerical integration method for efficient and accurate stability prediction. The variable-step technique is emphasized here to expand the numerical integration method,especially for the low radial immersion cases with multiple delays. First,the calculation accuracy of the numerical integration method is discussed and the variable-step algorithm is developed for milling stability prediction for single-delay and multiple-delay cases,respectively. The milling stability with variable pitch cutter is analyzed and the result is compared with those predicted with the frequency domain method and the improved full-discretization method. The influence of the runout effect on the stability boundary is investigated by the presented method. The numerical simulation shows that the cutter runout effect increases the stability boundary,and the increasing stability limit is verified by the milling chatter experimental results in the previous research. The numerical and experiment results verify the validity of the proposed method.

Analysis of the second-order BDF scheme with variable steps for the molecular beam epitaxial model without slope selection

In this work,we are concerned with the stability and convergence analysis of the second-order backward difference formula(BDF2)with variable steps for the molecular beam epitaxial model without slope selection.We first show that the variable-step BDF2 scheme is convex and uniquely solvable under a weak time-step constraint.Then we show that it preserves an energy dissipation law if the adjacent time-step ratios satisfy r_(k):=τ_(k)/τ_(k-1)<3.561.Moreover,with a novel discrete orthogonal convolution kernels argument and some new estimates on the corresponding positive definite quadratic forms,the L^(2)norm stability and rigorous error estimates are established,under the same step-ratio constraint that ensures the energy stability,i.e.,0<r_(k)<3.561.This is known to be the best result in the literature.We finally adopt an adaptive time-stepping strategy to accelerate the computations of the steady state solution and confirm our theoretical findings by numerical examples.