A METHOD OF INCREASING THE NUMERICAL STABILITY OF THE WAVE FIELD CALCULATION BASED ON THE MATRIX METHOD OF THOMSON?HASKELL

Purpose. It is known that the area of application of the Thomson?Haskell calculation scheme is in principle limited during the direct computer realization. The purpose of this work is the development of approaches to fix this problem for the case of arbitrary dipole effective-point source which is located in any layer of horizontal-layered half-space. Design/methodology/approach. The suggested approach is based on the method of transition from image of the wave field using characteristic matrices of the 4th order (method of Thomson?Haskell) to image, using characteristic matrices of the 6th order. Findings. The result of this work is the wave field presentation which does not include the characteristic matrices of the 4th order. Each minor of these matrices is replaced with the corresponding element of characteristic matrices of the 6th order. There is a doubling of columns and rows in the resulting matrices of the 6th order. Because of this, the order of characteristic matrices is lowered to five by reducing them to a quasi diagonal form. This will save computer time during calculation of interference fields. Practical value/implications. The developed approach increases the numerical stability during the calculation of wave fields. As a result, the computation of synthetic seismograms can be done for very thick layers of medium and higher wave frequencies.