Stability of finite difference numerical simulations of acoustic logging-while-drilling with different perfectly matched layer schemes

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(MPML), 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 d0. The optimal parameter space for the maximum value of the linear frequency-shifted factor (α0) 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.

[1]  Keiiti Aki,et al.  Discrete wave-number representation of seismic-source wave fields , 1977, Bulletin of the Seismological Society of America.

[2]  D. Komatitsch,et al.  An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation , 2007 .

[3]  Jean Virieux,et al.  SH-wave propagation in heterogeneous media; velocity-stress finite-difference method , 1984 .

[4]  Weng Cho Chew,et al.  A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates , 1994 .

[5]  Weng Cho Chew,et al.  A three-dimensional finite difference code for the modeling of sonic logging tools , 1998 .

[6]  Kristel C. Meza-Fajardo,et al.  A Nonconvolutional, Split-Field, Perfectly Matched Layer for Wave Propagation in Isotropic and Anisotropic Elastic Media: Stability Analysis , 2008 .

[7]  Z. Alterman,et al.  Propagation of elastic waves in layered media by finite difference methods , 1968 .

[8]  Stephen D. Gedney,et al.  Convolution PML (CPML): An efficient FDTD implementation of the CFS–PML for arbitrary media , 2000 .

[9]  K. Marfurt.  Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations , 1984 .

[10]  Robert L. Higdon,et al.  Radiation boundary conditions for elastic wave propagation , 1990 .

[11]  Tsili Wang,et al.  Finite-difference modeling of elastic wave propagation: A nonsplitting perfectly matched layer approach , 2003 .

[13]  Wei Zhang,et al.  Unsplit complex frequency-shifted PML implementation using auxiliary differential equations for seismic wave modeling , 2010 .

[14]  N. Chen.  Borehole wave propagation in isotropic and anisotropic media : three-dimensional finite difference approach , 1994 .

[15]  Jean-Pierre Berenger,et al.  Improved PML for the FDTD solution of wave-structure interaction problems , 1997 .

[16]  Xiaoming Tang,et al.  Investigation Of Lwd Quadrupole Shear Measurement In Real Environments , 2003 .

[17]  Jianlin Zhu.  A transparent boundary technique for numerical modeling of elastic waves , 1999 .

[18]  Hiroshi Kawase,et al.  Time-domain response of a semi-circular canyon for incident SV, P, and Rayleigh waves calculated by the discrete wavenumber boundary element method , 1988 .

[19]  Warwick D. Smith.  A nonreflecting plane boundary for wave propagation problems , 1974 .

[20]  Ekkehart Tessmer,et al.  3-D elastic modeling with surface topography by a Chebychev spectral method , 1994 .

[21]  C. Tsogka,et al.  Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media , 2001 .

[22]  John Lysmer,et al.  A Finite Element Method for Seismology , 1972 .

[23]  J. Virieux.  P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method , 1986 .

[24]  Q.H. Liu,et al.  Using perfectly matched layers for elastodynamics , 1996, IEEE Antennas and Propagation Society International Symposium. 1996 Digest.

[25]  Maxim Dmitriev,et al.  Application of M-PML reflectionless boundary conditions to the numerical simulation of wave propagation in anisotropic media. Part I: Reflectivity , 2011 .

[26]  C. Torres‐Verdín,et al.  HP-adaptive multi-physics finite-element simulation of wireline borehole sonic waveforms , 2011 .

[27]  D. Schmitt,et al.  Full‐wave acoustic logging: synthetic microseismograms and frequency‐wavenumber analysis , 1985 .

[28]  E. A. Skelton,et al.  Guided elastic waves and perfectly matched layers , 2007 .

[29]  Roland Martin,et al.  A High-Order Time and Space Formulation of the Unsplit Perfectly Matched Layer for the Seismic Wave Equation Using Auxiliary Differential Equations (ADE-PML) , 2010 .

[30]  Raj Mittra,et al.  Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers , 1996 .

[31]  Guo Tao,et al.  Wavefield simulation and analysis with the finite-element method for acoustic logging while drilling in horizontal and deviated wells , 2013 .

[32]  Xiaojun Huang.  Effects of Tool Positions on Borehole Acoustic Measurements: a Stretched Grid Finite Difference Approach , 2003 .

[33]  Andrew L. Kurkjian,et al.  Acoustic multipole sources in fluid-filled boreholes , 1986 .

[34]  James S. Sochacki,et al.  Absorbing boundary conditions and surface waves , 1987 .

[35]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[36]  Thomas L. Marzetta,et al.  Semblance processing of borehole acoustic array data , 1984 .

[37]  M. Nafi Toksöz,et al.  Elastic wave propagation in a fluid-filled borehole and synthetic acoustic logs , 1981 .

[38]  J. Sochacki.  Absorbing boundary conditions for the elastic wave equations , 1988 .

[39]  J. Carcione.  The wave equation in generalized coordinates , 1994 .

[40]  R. Madariaga.  Dynamics of an expanding circular fault , 1976, Bulletin of the Seismological Society of America.

[41]  Moshe Reshef,et al.  A nonreflecting boundary condition for discrete acoustic and elastic wave equations , 1985 .

[42]  G. Tao,et al.  Wavefield simulation and data-acquisition-scheme analysis for LWD acoustic tools in very slow formations , 2011 .

[43]  Xiuming Wang,et al.  Finite difference modelling of dipole acoustic logs in a poroelastic formation with anisotropic permeability , 2013 .

[44]  X. He,et al.  Finite-difference modeling of the monopole acoustic logs in a horizontally stratified porous formation. , 2009, The Journal of the Acoustical Society of America.