声波方程变网格有限差分正演模拟的虚假反射分析

摘要

变网格正演模拟方法因其低存储、高效率等特点,在含低速带或小尺度异常体的地下介质的数值模拟中发挥了重大作用.变网格模拟算法的精度与模拟时采用的频率、细网格长度、网格比、介质速度等参数密切相关.当模拟参数选取不合理时,在粗、细网格的分界处会产生强能量虚假反射,严重降低波场模拟的精度及分辨率.因此,选择合理的参数以压制虚假反射的产生在变网格正演模拟中显得尤为重要.本文基于一维假设和平面波理论,结合变系数变网格算法,推导出在变网格声波方程正演模拟中由于网格变化所引起的虚假反射的反射率函数公式,该公式为反射率函数与频率、细网格步长、网格比、介质速度的函数.基于该函数表达式,本文通过理论分析和模型数值算例,系统分析了频率、细网格步长、网格比、介质速度等对于变网格数值模拟中虚假反射的影响.理论分析和数值模拟实验均表明:在声波方程变网格正演模拟中,到达临界位置(即虚假反射率函数为1)之前,频率越高、细网格步长越大、网格比越大、介质速度越低,虚假反射率越大,即虚假反射现象越明显;反之,虚假反射率越小,虚假反射越弱.本文推导出的反射率函数可为实际变网格数值模拟中的各参数选取提供理论指导.