1、参数初始化； 2、形成边界削波数据； 3、波场初始化； 4、Zmax层波场延拓到深度Zmax-1； 5、Zi+1层波场延拓到深度Zi； 6、重复5，从Iz Nz-1开始，直到Iz 1，得测线上的频率—空间域波场； 7、频率-空间域波场对频率做反傅里叶变换，得时间-空间波场； -1, parameter initialization 2, the formation of boundary clipping data 3, wave field init