In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid colu...In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns (radius -27 mm). Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole (the difference is 〉30 times), the stability and efficiency of the perfectly matched layer (PML) scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain (FDTD) in the ALWD simulation, including field-splitting PML (SPML), multiaxial PML(M- PML), non-splitting PML (NPML), and complex frequency-shifted PML (CFS-PML). The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor (a0) and the scaling factor (β0) depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.展开更多
Using staggered-grid finite difference method to solve seismic wave equation,large spatial grid and high dominant frequency of source cause numerical dispersion,staggeredgrid finite difference method,which can reduce ...Using staggered-grid finite difference method to solve seismic wave equation,large spatial grid and high dominant frequency of source cause numerical dispersion,staggeredgrid finite difference method,which can reduce the step spatial size and increase the order of difference,will multiply the calculation amount and reduce the efficiency of solving wave equation.The optimal nearly analytic discrete(ONAD)method can accurately solve the wave equation by using the combination of displacement and gradient of spatial nodes to approach the spatial partial derivative under rough grid and high-frequency condition.In this study,the ONAD method is introduced into the field of reverse-time migration(RTM)for performing forward-and reverse-time extrapolation of a two-dimensional acoustic equation,and the RTM based on ONAD method is realized via normalized cross-correlation imaging condition,effectively suppressed the numerical dispersion and improved the imaging accuracy.Using ONAD method to image the groove model and SEG/EAGE salt dome model by RTM,and comparing with the migration sections obtained by staggered-grid finite difference method with the same time order 2 nd and space order 4 th,results show that the RTM based on ONAD method can effectively suppress numerical dispersion caused by the high frequency components in source and shot records,and archive accurate imaging of complex geological structures especially the fine structure,and the migration sections of the measured data show that ONAD method has practical application value.展开更多
The phenomenon of the resonant activation (RA) of a particle over a fluctuating potential barrier with a four-value noise is investigated. It is shown that the mean first passage time (MFPT) displays six minima as...The phenomenon of the resonant activation (RA) of a particle over a fluctuating potential barrier with a four-value noise is investigated. It is shown that the mean first passage time (MFPT) displays six minima as the function of the transition rates γ1, γ2, γ3, γ4, γ5, and 76 of the four-value noise, respectively. In addition, the effect of other parameters of the system, such as the noise strength D of the additive Gaussian white noise and the parameter value a, b, c, and d of the four-value noise, on the RAs is also investigated.展开更多
Time-frequency peak filtering (TFPF) is highly efficient in suppressing random noise in seismic data. Although the hypothesis of stationary Gaussian white noise cannot be fulfilled in practical seismic data, TFPF can ...Time-frequency peak filtering (TFPF) is highly efficient in suppressing random noise in seismic data. Although the hypothesis of stationary Gaussian white noise cannot be fulfilled in practical seismic data, TFPF can effectively suppress white and colored random noise with different intensities, as can be theoretically demonstrated by detecting such noise in synthetic seismic data. However, a "zero-drift" effect is observed in the filtered signal and is independent of the average power and variance of the random noise, but related to its mean value. Furthermore, we consider the situation where the local linearization of the seismic data cannot be satisfied absolutely and study the "distortion" characteristics of the filtered signal using TFPF on a triangular wave. We found that over-compensation is possible in the frequency band for the triangular wave. In addition, it is nonsymmetrical and has a relationship to the time-varying curvature of the seismic wavelet. The results also present an improved approach for TFPF.展开更多
Full annulus simulations of the flow which develops in a transonic centrifugal compressor are performed at two stable operating points (peak efficiency and near surge) and during the path to surge. At stable conditi...Full annulus simulations of the flow which develops in a transonic centrifugal compressor are performed at two stable operating points (peak efficiency and near surge) and during the path to surge. At stable conditions, the flow field properties are analyzed by comparisons with experimental data and numerical simulations using a phase lagged approach previously carried out. Regarding the stage overall performance, an excellent agreement is obtained between the numerical results (both with time lagged approach and full-annulus calculation) and the ex- periments. From the full-annnlus simulations, the change in flow pattern from peak efficiency to surge is found to be perfectly similar to that obtained from the simulations using the time lagged approach. In particular, pro- vided that the operating point is stable, the flow proves to be chorochronic. The full-annulus simulations were continued after a unique small change in the throttle law applied at the exit of the numerical domain. The mass flow, pressure ratio and efficiency then significantly drop all the more the time progresses. The simulation becomes unstable and the surge inception well underway. The path to surge is found to be due to the enlargement of the boundary layer separation on the suction side of the diffuser vanes in accordance with the conclusions drawn from the chorochronic simulations and experiments. But as the time progresses, the flow loses its chorochronic character. Stall cells rotating at around 7% of the rotor speed develop and lead to surge in around 5 revolutions.展开更多
We report two models of the lateral displacement of acoustic-wave scattering on a fluid-solid interface that reveal an acoustic analog of the Goos-Hainchen effect in optics. This acoustic analog is called the acoustic...We report two models of the lateral displacement of acoustic-wave scattering on a fluid-solid interface that reveal an acoustic analog of the Goos-Hainchen effect in optics. This acoustic analog is called the acoustic Goos-Hainchen effect. Using newly proposed models, we made numerical calculations for the system ofa water-Perspex interface. Specifically, in the post-critical-angle region, we observed a lateral displacement (and transition time) of the reflected P-wave with respect to the incident P-wave. The first arrival of the acoustic signal from the interface is found to be a reflected P-wave rather than the sliding-refraction P-wave usually described in traditional acoustic-logging sliding P-wave theory. For both proposed models, the effective propagation speed of the reflected P-wave along the interface depends on not only the physical properties of the interracial media but also the incident angle. These observations are intriguing and warrant further investigation.展开更多
基金supported by NSFC(No.41174118)one of the major state S&T special projects(No.2008ZX05020-004)+1 种基金a Postdoctoral Fellowship of China(No.2013M530106)China Scholarship Council(No.2010644006)
文摘In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns (radius -27 mm). Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole (the difference is 〉30 times), the stability and efficiency of the perfectly matched layer (PML) scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain (FDTD) in the ALWD simulation, including field-splitting PML (SPML), multiaxial PML(M- PML), non-splitting PML (NPML), and complex frequency-shifted PML (CFS-PML). The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor (a0) and the scaling factor (β0) depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.
基金financially supported by the National Key R&D Program of China(No.2018YFC1405900)the National Natural Science Foundation of China(No.41674118)+1 种基金the Fundamental Research Funds for the Central Universities(No.201822011)the National Science and Technology Major Project(No.2016ZX05027-002)。
文摘Using staggered-grid finite difference method to solve seismic wave equation,large spatial grid and high dominant frequency of source cause numerical dispersion,staggeredgrid finite difference method,which can reduce the step spatial size and increase the order of difference,will multiply the calculation amount and reduce the efficiency of solving wave equation.The optimal nearly analytic discrete(ONAD)method can accurately solve the wave equation by using the combination of displacement and gradient of spatial nodes to approach the spatial partial derivative under rough grid and high-frequency condition.In this study,the ONAD method is introduced into the field of reverse-time migration(RTM)for performing forward-and reverse-time extrapolation of a two-dimensional acoustic equation,and the RTM based on ONAD method is realized via normalized cross-correlation imaging condition,effectively suppressed the numerical dispersion and improved the imaging accuracy.Using ONAD method to image the groove model and SEG/EAGE salt dome model by RTM,and comparing with the migration sections obtained by staggered-grid finite difference method with the same time order 2 nd and space order 4 th,results show that the RTM based on ONAD method can effectively suppress numerical dispersion caused by the high frequency components in source and shot records,and archive accurate imaging of complex geological structures especially the fine structure,and the migration sections of the measured data show that ONAD method has practical application value.
基金National Natural Science Foundation of China under Grant No.10375009K.C.Wong Magna Fund in Ningbo University of ChinaNational Natural Science Foundation of China under Grant No.10647134
文摘The phenomenon of the resonant activation (RA) of a particle over a fluctuating potential barrier with a four-value noise is investigated. It is shown that the mean first passage time (MFPT) displays six minima as the function of the transition rates γ1, γ2, γ3, γ4, γ5, and 76 of the four-value noise, respectively. In addition, the effect of other parameters of the system, such as the noise strength D of the additive Gaussian white noise and the parameter value a, b, c, and d of the four-value noise, on the RAs is also investigated.
基金supported by the National Natural Science Foundation of China (Grant Nos.40574051,41130421,40930418 & 40974064)the National Special Project of Science and Technology of China (Grant No.Sinoprobe-03)
文摘Time-frequency peak filtering (TFPF) is highly efficient in suppressing random noise in seismic data. Although the hypothesis of stationary Gaussian white noise cannot be fulfilled in practical seismic data, TFPF can effectively suppress white and colored random noise with different intensities, as can be theoretically demonstrated by detecting such noise in synthetic seismic data. However, a "zero-drift" effect is observed in the filtered signal and is independent of the average power and variance of the random noise, but related to its mean value. Furthermore, we consider the situation where the local linearization of the seismic data cannot be satisfied absolutely and study the "distortion" characteristics of the filtered signal using TFPF on a triangular wave. We found that over-compensation is possible in the frequency band for the triangular wave. In addition, it is nonsymmetrical and has a relationship to the time-varying curvature of the seismic wavelet. The results also present an improved approach for TFPF.
基金the HPC resources of CINES under the allocation 2012- 2a6356 and 2013-2a6356
文摘Full annulus simulations of the flow which develops in a transonic centrifugal compressor are performed at two stable operating points (peak efficiency and near surge) and during the path to surge. At stable conditions, the flow field properties are analyzed by comparisons with experimental data and numerical simulations using a phase lagged approach previously carried out. Regarding the stage overall performance, an excellent agreement is obtained between the numerical results (both with time lagged approach and full-annulus calculation) and the ex- periments. From the full-annnlus simulations, the change in flow pattern from peak efficiency to surge is found to be perfectly similar to that obtained from the simulations using the time lagged approach. In particular, pro- vided that the operating point is stable, the flow proves to be chorochronic. The full-annulus simulations were continued after a unique small change in the throttle law applied at the exit of the numerical domain. The mass flow, pressure ratio and efficiency then significantly drop all the more the time progresses. The simulation becomes unstable and the surge inception well underway. The path to surge is found to be due to the enlargement of the boundary layer separation on the suction side of the diffuser vanes in accordance with the conclusions drawn from the chorochronic simulations and experiments. But as the time progresses, the flow loses its chorochronic character. Stall cells rotating at around 7% of the rotor speed develop and lead to surge in around 5 revolutions.
基金the Xi’an University of Posts and Telecommunicationsthe Physical Sciences Division at the University of Chicagothe Scientific Research Program(Grant No.15JK1685)of the Shaanxi Provincial Education Department
文摘We report two models of the lateral displacement of acoustic-wave scattering on a fluid-solid interface that reveal an acoustic analog of the Goos-Hainchen effect in optics. This acoustic analog is called the acoustic Goos-Hainchen effect. Using newly proposed models, we made numerical calculations for the system ofa water-Perspex interface. Specifically, in the post-critical-angle region, we observed a lateral displacement (and transition time) of the reflected P-wave with respect to the incident P-wave. The first arrival of the acoustic signal from the interface is found to be a reflected P-wave rather than the sliding-refraction P-wave usually described in traditional acoustic-logging sliding P-wave theory. For both proposed models, the effective propagation speed of the reflected P-wave along the interface depends on not only the physical properties of the interracial media but also the incident angle. These observations are intriguing and warrant further investigation.