Acoustic viscoelastic modeling by frequency-domain boundary element method

Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element methods. Targeted at viscoelastic numerical modeling for multilayered media, the constant-Q acoustic wave equation is transformed into the corresponding wave integral representation with its Green＇s function accounting for viscoelastic coefficients. An efficient alternative for full-waveform solution to the integral equation is proposed in this article by extending conventional frequency-domain boundary element methods to viscoelastic media. The viscoelastic boundary element method enjoys a distinct characteristic of the explicit use of boundary continuity conditions of displacement and traction, leading to a semi-analytical solution with sufficient accuracy for simulating the viscoelastic effect across irregular interfaces. Numerical experiments to study the viscoelastic absorption of different Q values demonstrate the accuracy and applicability of the method.

Dynamic interaction numerical models in the time domain based on the high performance scaled boundary finite element method

Consideration of structure-foundation-soil dynamic interaction is a basic requirement in the evaluation of the seismic safety of nuclear power facilities. An efficient and accurate dynamic interaction numerical model in the time domain has become an important topic of current research. In this study, the scaled boundary finite element method （SBFEM） is improved for use as an effective numerical approach with good application prospects. This method has several advantages, including dimensionality reduction, accuracy of the radial analytical solution, and unlike other boundary element methods, it does not require a fundamental solution. This study focuses on establishing a high performance scaled boundary finite element interaction analysis model in the time domain based on the acceleration unit-impulse response matrix, in which several new solution techniques, such as a dimensionless method to solve the interaction force, are applied to improve the numerical stability of the actual soil parameters and reduce the amount of calculation. Finally, the feasibility of the time domain methods are illustrated by the response of the nuclear power structure and the accuracy of the algorithms are dynamically verified by comparison with the refinement of a large-scale viscoelastic soil model.

Stability conditions of explicit integration algorithms when using 3D viscoelastic artificial boundaries

Viscoelastic artificial boundaries are widely adopted in numerical simulations of wave propagation problems.When explicit time-domain integration algorithms are used,the stability condition of the boundary domain is stricter than that of the internal region due to the influence of the damping and stiffness of an viscoelastic artificial boundary.The lack of a clear and practical stability criterion for this problem,however,affects the reasonable selection of an integral time step when using viscoelastic artificial boundaries.In this study,we investigate the stability conditions of explicit integration algorithms when using three-dimensional(3D)viscoelastic artificial boundaries through an analysis method based on a local subsystem.Several boundary subsystems that can represent localized characteristics of a complete numerical model are established,and their analytical stability conditions are derived from and further compared to one another.The stability of the complete model is controlled by the corner regions,and thus,the global stability criterion for the numerical model with viscoelastic artificial boundaries is obtained.Next,by analyzing the impact of different factors on stability conditions,we recommend a stability coefficient for practically estimating the maximum stable integral time step in the dynamic analysis when using 3D viscoelastic artificial boundaries.

The Analysis and Comparison of All Kinds of Buried Pipeline Model Based on Seismic Effect

The problem of seismic response of buried pipeline aimed at the interaction of soil around the pipeline and the complicated calculation model was considered, and the various simplified finite element model was calculated, and it was analyzed. Firstly, the ADINA finite element analysis software was used. The four nodes in shell unit were used by tube. The spring unit was used by soil spring. The analysis model of buried pipeline finite element numerical based on tension and compression spring was established. Seismic wave was input. The response to the simple boundary, viscoelastic boundary earthquake were calculated and analyzed by the finite element numerical simulation. The pipeline’s earthquake ground motion response was obtained, and was compared with the real soil model, and the most suitable simplified calculation model for numerical analysis of buried pipeline was found, which was the numerical analysis model of buried pipeline pressure spring tension finite element based on the viscoelastic boundary, and the theory basis for the seismic design of pipeline was provided.

An acoustic-wave theory for casing bond evaluation using viscoelastic slip boundary modeling

Slip boundary condition is commonly utilized to model elastic wave propagation through layered earth media. The same approach is used here to characterize acoustic wave propagation along a cased borehole with various cement bond conditions. By modeling the cement layer between casing and formation as a viscoelastic slip interface with complex coupling rigidity parameters, one can not only reduce the complexity in the classical elastic wave modeling of the problem, but also efficiently model various complicated wave phenomena that are difficult for the existing modeling. More specifically, the new theory can well describe the effect of the cement bond condition change and the location of the change(i.e., whether it is in the first interface between casing and cement, or the second interface between cement and formation) on the acoustic waves,demonstrating the good modeling capability and predicting power. Application of the theory to field data shows that the theory can correctly model the acoustic wave characteristics and interpret the cement bond condition, thus providing a useful fundament theory for casing bond evaluation using acoustic logging.

GENERAL DECAY RATE ESTIMATES FOR VISCOELASTIC WAVE EQUATION WITH VARIABLE COEFFICIENTS

The authors study decay properties of solutions for a viscoelastic wave equation with variable coefficients and a nonlinear boundary damping by the differential geometric approach.