从地震序列的模拟看数值地震预报的前景

汶川大地震的伤痛还未愈合，又发生了玉树（Ms=7．1）、智利（Mw=8．8）、海地（Ms=7．5）及日本（Mw=9．0）等特大地震。地震难道就不能预报？用什么方法可以准确预报地震？这些问题再一次成为科学界与社会关注的热点。本文以发生于2008年汶川地震的龙门山断裂带上的强震序列的模拟为出发点，通过基于严格数学、力学原理的有限元方法计算给出龙门山断裂带上大地震的时间序列、同震位移场及同震应力场。通过研究发现，计算得到的龙门山断裂带上的大地震平均复发间隔为～3300年，这与古地震调查及其他方法给出的结果有一致性；模拟的同震应力场形态与汶川余震的深度分布之间存在着关联，同震位移与观测结果符合得较好。值得特别指出的是，文中的模型能够模拟断层由闭锁到解锁产生地震的全过程。模拟结果进一步显示，完全有可能利用科学计算方法给出未来强震发生的时、空、强三要素。类似于气象预报，本文将这种预测地震的方法称为数值地震预报。最后，文中对数值地震预报的可能性及需要解决的关键问题等进行了深入分析和研究，认为数值地震预报在科学上是可能的、技术上是可行的，数值地震预报是成功预报地震的最有效途径。