地震数值计算中台阶状倾斜地层效应的消除

为了消除采用矩形网格的地震波正演计算中,台阶状倾斜界面引发的绕射波和网状干扰波,Zeng提出网格细分加权平滑方法。本文对该方法中加权函数的性质和参数的求取进行了讨论,给出了加权函数的图形和参数选取的办法。规格化因子C可以通过非规格化的加权函数计算。根据需要,定义了表征加权函数平滑度的参数s。由确定的网格大小与s计算衰减因子σ,就可以确定加权函数。对于含有一个倾斜界面的模型进行了数值模拟。一般情形下,参数s必须大于或等于0.2。通过对不同倾角倾斜界面的数值计算,显示了相同参数选择对于地层倾角变化的适应性。