【摘要】 频率空间域地震波数值模拟具有独特的优势:可以同时模拟多源的波传播、每个频率之间独立并行地计算、计算频带选择灵活、不存在累计误差、容易模拟粘弹性介质中地震波传播。但是该方法的最大瓶颈是对于计算机内存的巨大需求。我们使用压缩存储系数矩阵的方法,极大地减少了计算机内存的需求量。同时为了减少短差分算子的数值频散,引用了频率空间域25点弹性波波动方程的差分格式,并使用了最小二乘意义下求出的优化差分系数。为了克服边界反射,采用了最佳匹配层吸收边界条件。数值模拟试验证明:用压缩存储系数矩阵及优化差分系数的频率空间域25点差分格式进行弹性波正演模拟,可以减少数值频散,提高计算精度。使用较大的网格间距,降低计算机内存需求,并保持较高的计算效率。该正演方法为后续弹性波偏移和弹性参数反演提供较好的基础。更多还原

【Abstract】 Numerical simulation in the frequency-space domain has inherent advantages,such as:it is possible to simulate wave propagation from multiple sources simultaneously;there are no cumulative errors;only the interesting frequencies can be selected;and it is more suitable for wave propagation in viscoelastic media.The only obstacle to using the method is the requirement of huge computer storage.We extend the compressed format for storing the coefficient matrix.It can reduce the required computer storage dramatically.We get the optimal coefficients by least-squares method to suppress the numerical dispersion and adopt the perfectly matched layer(PML) boundary conditions to eliminate the artificial boundary reflections.Using larger grid intervals decreases computer storage requirements and provides high computational efficiency.Numerical experiments demonstrate that these means are economic and effective,providing a good basis for elastic wave imaging and inversion.更多还原