一种黏声介质有限差分地震波复走时计算方法

本发明公开了一种黏声介质有限差分地震波复走时计算方法，先建立黏声介质地震波复走时计算公式，接着对所需探测区域的速度模型进行初始化，实现了对震源点以及其余网格点分别进行初值设置，利用因式分解思想，将速度模型的非吸收衰减走时分成非吸收衰减走时固定项和非吸收衰减走时扰动项；慢度分解成慢度固定项和慢度扰动项，良好地解决了震源奇异性问题，降低了误差传播影响，提高了地震波走时的计算精度；后续采用Gauss‑Seidel扫描计算，得到各个网格点的非吸收衰减走时及耗散走时，最终得到各个位置的黏声介质地震波复走时；在该过程中采用了快速扫描方案，相对与传统算法节省了计算时间，有效提高了计算效率。