EI、Scopus 收录
中文核心期刊

波动分量耦合效应引起的人工边界大反射问题

LARGE REFLECTION ON AN ARTIFICIAL BOUNDARY CAUSED BY THE COUPLING EFFECT OF DIFFERENT WAVE COMPONENTS

  • 摘要: 无限域空间有限化的人工边界条件(artificial boundary condition, ABC)对于确保地震波动数值模拟的精度和计算效率非常关键. 廖氏透射边界和Higdon单程波方程等外推型人工边界已被广泛使用, 但是弹性P波-S波耦合作用下该类ABC的基本性能尚未被充分揭示. 通过将外推型ABC的弹性波反射系数拆解为与纵横波速比cp/cs相关的物理因子和与边界计算波速caj (j = 1, ···, N)相关的N个(N为边界阶次)控制因子的乘积, 详细分析了参数取值处于cp/cs = 1.5 ~ 4, N = 1, 2, 3, caj = cs + q × (cpcs), q = 0, 0.1, 0.2, ···, 1.0范围内边界反射系数的变化规律. 研究发现当所有计算波速caj取值都接近cp时, S-P反射模式常常会出现一种幅值超过1(对应于100%反射) 的巨大反射, 导致反射误差随着边界阶次升高不降反增. 这种由于不同波动分量之间耦合效应所引起的人工边界大反射问题打破了廖氏透射边界和Higdon边界已有认识的增加边界阶次总能稳步提高精度的预期, 在理论和应用中应予以重视. 以反射系数的85%分位值作为精度指标, 分别对可能导致耦合大反射和能够导出高精度ABC的边界参数取值范围进行了界定. 数值试验证实了不同计算参数取值对外推型ABC精度的影响规律, 并与黏弹性ABC进行了比较.

     

    Abstract: Artificial boundary conditions (ABCs) for elastic waves play a key role in improving the accuracy and computational efficiency of numerical simulations of seismic wave propagation. Extrapolation-type artificial boundary conditions, including Liao’s multi-transmitting formula, Higdon’s one-way wave equations and etc., have been widely used in this field. However, the fundamental characteristics of these ABCs affected by the coupling effect of elastic P- and S-waves have, up to now, not been completely revealed. This work factorizes each of the elastic-wave reflection coefficients of extrapolation-type ABCs into the product of a physical part and a set of control parts. The physical part is a function of the ratio of compressional and shear waves cp/cs, which is the physical property of the elastic wave problem. The control parts are functions of the ABC’s computational velocities caj (j = 1, ···, N) and boundary order N, which are the user-defined control parameters of the ABC. Using this factorization, the characteristics of the boundary reflection coefficients are thoroughly analyzed across the following parameter ranges: cp/cs = 1.5 ~ 4, N = 1, 2, 3, caj = cs + q × (cpcs), with j = 1,···, N and q = 0, 0.1, 0.2, ···, 1.0. The results indicate a key finding that when all computational velocities caj’s approach the value of cp, the S-P reflection mode often exhibits a significant spurious reflection with an amplitude exceeding 1 (corresponding to 100% reflection), and thus causes the boundary reflection error to increase with the boundary order N, rather than decease. This large spurious reflection, caused by the coupling effect of different wave components, goes against the existing understanding in Liao’s and Higdon’s boundary conditions that increasing the boundary order always makes the solution more accurate. Consequently, this issue must be paid special attention in both theory and application. Using the 85% percentile value of reflection coefficients as an accuracy metric, the applicable ranges of boundary parameters that may cause the above-mentioned large spurious reflection and those can lead to high-precision ABCs are delineated, respectively. Numerical experiments confirm the previously revealed influence of different computational parameters on the accuracy of extrapolation-type ABCs, and present a comparison with the prevalent viscous-spring boundary (a typical stress-type ABC).

     

/

返回文章
返回