CAJ | 학술논문

汪荣江提出一个简单的正交归一化技术来克服经典的Thomson-Haskell传播矩阵方法中存在的数值不稳定问题.为了进一步提高计算效率,给出该方法的2种改进.一种改进方法是将传播矩阵中与频率无关的部分分离出来,对于某一固定的水平慢度,这些矩阵只需计算一次;另一个改进是利用Langer块对角化的技术,将传播矩阵分解为几个稀疏矩阵的乘积.我们将改进之后的算法应用于计算水平分层模型中的广义反射系数.较之原有方案,提出的改进能节省一半计算时间.
왕영강제출일개간단적정교귀일화기술래극복경전적Thomson-Haskell전파구진방법중존재적수치불은정문제.위료진일보제고계산효솔,급출해방법적2충개진.일충개진방법시장전파구진중여빈솔무관적부분분리출래,대우모일고정적수평만도,저사구진지수계산일차;령일개개진시이용Langer괴대각화적기술,장전파구진분해위궤개희소구진적승적.아문장개진지후적산법응용우계산수평분층모형중적엄의반사계수.교지원유방안,제출적개진능절성일반계산시간.
Wang proposed a simple orthonormalization technique to overcome the numerical instability in the original Thomson-Haskell propagator algorithm. In this paper, two kinds of improvements to the method for computational efficiency are presented. One is to separate the frequency independent matrices from the propagator, these terms can be computed once for a fixed slowness; the other is to decompose the propagator by the Langer block-diagonalization into a product of several sparse matrices. The new developments are applied to calculate the generalized reflection coefficients for a multi-layered model. It leads nearly half reduction in computation time compared with the original scheme.