[发明专利]一种地震震源反演不确定性分析方法、存储介质及服务器

说明书

本发明公开一种地震震源反演不确定性分析方法、存储介质及服务器，所述方法包括步骤：采用贝叶斯统计方法对震源反演参数进行解算获取相应震源反演参数的先验概率密度函数；根据所述震源反演参数的先验概率密度函数获取震源信息量；根据所述震源信息量得出随机不确定性函数以及认知不确定性函数；结合所述随机不确定性函数以及认知不确定性函数对震源反演初步结果进行评估判定，得出准确的震源反演解译结果。本发明通过引入信息熵理论，将这两类不确定性统一在信息量的基准上进行衡量，从而实现了两类不确定性定量化评估以及统一，为InSAR地震震源反演结果正确解译提供科学依据。

技术领域

本发明涉及震源反演领域，尤其涉及一种地震震源反演不确定性分析方法、存储介质及服务器。

背景技术

InSAR是近几十年发展的先进空间对地测量技术，是传统的SAR遥感技术与射电天文干涉技术相结合的产物，该技术可获取近距离和高密度地震位移场的观测数据。基于InSAR获取的震源空间信息更加精细准确，在现有的地震工程相关应用也起到很积极的作用。

虽然InSAR获取的地震信息在空间分辨率上具备很大优势，但是必须认识到InSAR断层滑动反演仍然存在多解性，基于不同方法和数据反演得到的滑动模型差别非常大。在地震危险性分析中，基于不同InSAR反演结果得出的结论差异很大。发表在顶级期刊NatureGeoscience和Natrue上的两篇关于2010年智利8.8级大地震的文章，对于该地震是否完全释放了震间积累的能量给出了截然不同的结论。

因此，对InSAR地震震源不确定性进行分析研究是InSAR震源目标标准化构建的基础，对保障地震危险性评估等后续应用可靠性有着重要意义。

发明内容

鉴于上述现有技术的不足，本发明的目的在于提供一种地震震源反演不确定性分析方法、存储介质及服务器，旨在解决现有技术在使用InSAR 地震震源反演结果进行地震解译时，未同时考虑反演结果的随机不确定性和认知不确定性，导致解译结果不准确，阻碍对地震认知的问题。

本发明的技术方案如下：

一种地震震源反演不确定性分析方法，其中，包括步骤：

采用贝叶斯统计方法对震源反演参数进行解算获取相应震源反演参数的先验概率密度函数；

根据所述震源反演参数的先验概率密度函数获取震源信息量；

根据所述震源信息量得出随机不确定性函数以及认知不确定性函数；

结合所述随机不确定性函数以及认知不确定性函数对震源反演初步结果进行评估判定，得出准确的震源反演解译结果。

所述的地震震源反演不确定性分析方法，其中，所述震源反演参数包括：断层块滑动分布向量s、断层几何参数m、InSAR数据集权重σ²以及平滑系数α。

所述的地震震源反演不确定性分析方法，其中，所述震源信息量包括：震源反演参数需要的信息量、InSAR数据提供的信息量以及震源反演中利用的信息量。

所述的地震震源反演不确定性分析方法，其中，所述震源反演参数需要的信息量为震源反演参数先验概率密度函数的信息熵；所述InSAR数据提供的信息量为震源反演参数与InSAR观测数据之间的互信息；所述震源反演中利用的信息量为InSAR反演参数与震源反演参数之间的互信息。

所述的地震震源反演不确定性分析方法，其中，所述随机不确定性函数为震源反演参数需要的信息量减去InSAR数据提供的信息量。

所述的地震震源反演不确定性分析方法，其中，所述认知不确定性函数为InSAR数据提供的信息量减去震源反演中利用的信息量。