[发明专利]叠前弹性阻抗随机反演方法、装置、存储介质及处理器

说明书

本发明实施例提供一种叠前弹性阻抗随机反演方法、装置、处理器及存储介质。方法包括：确定统计特征参数，统计特征参数包括横向自相关长度、纵向自相关长度、自相关角度中的至少一者；根据统计特征参数确定对应的自相关函数；对自相关函数进行处理，加入统计特征参数的均值和标准差，以得到基于统计特征参数的先验信息模型；根据先验信息模型与预先获取的地震数据，建立对应的随机反演的目标函数；对先验信息模型进行更新，通过目标函数得到对应的弹性阻抗反演结果。

技术领域

本发明涉及地球物理反演领域，具体地涉及一种叠前弹性阻抗随机反演方法、装置、存储介质及处理器。

背景技术

地震随机反演结合地质统计先验信息，融合测井数据中的高频成分，提高反演结果的垂向分辨率，对地下薄储层的有效刻画具有一定的优势。地质统计先验信息是地震随机反演的重要组成部分，影响着先验模型构建的准确性及后续随机反演结果的可靠性。常规的地震随机反演以测井数据作为硬数据，地震数据作为软约束，充分利用测井数据中的高频信息，通过变差函数来随机模拟得到多个模拟实现，并与地震反演相结合得到多个反演解。

在实际的地震随机反演过程中，为了降低反演的复杂程度，通常假设地下介质是层状均匀介质，但是实际地下介质由大量复杂的非均匀介质组成。常规地震随机反演先验信息的构建主要利用测井数据获得的变差函数来表征储层的空间结构特性，对地震数据利用不足，而且变差函数不足以有效表征实际复杂地层的非均质特性。由于用于描述地下储层空间特征的参数有限，所构建的先验模型不能有效表征地下小尺度的变化，进而会影响后续地震随机反演的精度。

实际地下地层介质具有非均匀的特性，对于含有非均质特性的地下地层介质可以通过随机介质进行表征，随机介质可以理解为在大尺度均匀背景下分布着小尺度非均匀扰动的理论模型。

发明内容

本发明实施例的目的是提供一种叠前弹性阻抗随机反演方法、装置、存储介质及处理器。

为了实现上述目的，本发明第一方面提供一种叠前弹性阻抗随机反演方法，包括：

确定统计特征参数，统计特征参数包括横向自相关长度、纵向自相关长度、自相关角度中的至少一者；

根据统计特征参数确定对应的自相关函数；

对自相关函数进行处理，加入统计特征参数的均值和标准差，以得到基于统计特征参数的先验信息模型；

根据先验信息模型与预先获取的地震数据，建立对应的随机反演的目标函数；

对先验信息模型进行更新，通过目标函数得到对应的弹性阻抗反演结果。

在本发明的实施例中，确定统计特征参数包括：获取叠前地震数据；确定与叠前地震数据相对应的估计子波的功率谱；根据估计子波的功率谱确定使用弹性阻抗表征的非平稳随机介质的功率谱；对非平稳随机介质的功率谱进行傅里叶变换，得到随机介质的自相关函数；通过特征向量法和自相关函数确定统计特征参数。

在本发明的实施例中，叠前地震数据为至少三个部分角度叠加地震数据。

在本发明的实施例中，对自相关函数进行处理包括：对自相关函数进行傅里叶变换；通过维纳辛钦定理得到傅里叶变换后的自相关函数的振幅谱；在振幅谱中加入随机相位后，进行傅里叶反变换得到对应的空间随机扰动序列；对空间随机扰动序列进行归一化处理。

在本发明的实施例中，根据先验信息模型与预先获取的地震数据，建立对应的随机反演的目标函数包括：获取地震数据；将地震数据融入至似然函数；通过统计特征参数的正则化约束项减少反演解的误差，以此建议对应的随机反演的目标函数。

在本发明的实施例中，对先验信息模型进行更新包括：通过快速模拟退火算法对先验模型进行扰动更新。在本发明的实施例中，根据选用不同弹性参数表示的弹性阻抗方程以及弹性阻抗反演结果的数据确定出不同弹性参数的值。