【技术实现步骤摘要】
本专利技术涉及地球物理勘探，具体涉及储层属性反演方法、装置、设备及存储介质。
技术介绍
1、在油田开发中，地震数据的反演对于理解和解决储层连接性、流体流动障碍和隔板等问题至关重要。地震数据的反演通过量化储层压力、水饱和度和气饱和度的变化，为地下储层的管理提供关键信息；还可以通过对多个监测点的反演结果进行定量分析，以更精确地表征隔板的渗透性和位置。因此，地震数据的反演在油田开发中发挥着重要作用，为优化储层管理和生产策略提供了关键的技术支持。但是，受地震噪声、采集可重复性、地震成像精度等问题的影响，仅仅依靠地震数据反演存在着很大的不确定性。特别是在近海地震数据采集过程中，受涌浪等的影响，地震数据噪声对地震数据精度的影响更加明显。中国专利cn117471540a公开了一种基于叠后地震数据预测储层物性参数的方法、装置和设备，该方法根据叠后地震资料，利用提出的鲁棒谱分解方法得到不同频率的分频地震数据；根据不同频率的分频地震数据确定层间品质因子数据；对最终层间品质因子数据进行非线性同步反演得到孔隙度、渗透率、含气饱和度参数数据，与常规单项参数反演相比，更能够利用孔隙度、渗透率、含气饱和度三参数间的物理联系，提高反演结果的准确性。但是，该方法并未进行先验信息的约束，虽然通过对叠后地震数据精细处理提高了反演精度，但并不能克服地震资料本身精度的局限性。cn117388944a 公开了一种地质模型约束的多物性参数反演方法，首先依据井资料或先验地质知识构建基准岩石物理模型，使用地震正演方法计算基准岩石物理模型的地震响应，求取不同的实际观测数据相对于基
技术实现思路
1、针对现有技术存在的问题，本专利技术提供一种储层属性反演方法、装置、设备及存储介质，根据工区已有地质资料和测井数据建立先验模型，并将先验模型的数据同化到地震数据的储层属性反演中，对地震数据进行约束，提高储层属性反演的精度。
2、第一方面，本专利技术提供了一种储层属性反演方法，所述方法包括：
3、步骤s101：根据先验模型的储层属性估计，构建未知向量n的先验概率；
4、步骤s102：根据实际采集角道集地震数据的负振幅映射属性和先验模型模拟的角道集地震数据的负振幅映射属性，构建似然函数；
5、步骤s103：根据所述先验概率和所述似然函数，构建后验概率；
6、步骤s104：求解未知向量n，获取储层属性反演结果。
7、在一种可能的实现方式中，所述先验概率为：
8、，
9、其中，为先验概率，为未知向量，为先验模型的储层属性估计，为自身的协方差矩阵，i为单位矩阵，为概率因子；所述先验模型为根据工区已有地质资料和测井数据所建立的先验模型，所述储层属性包括压力变化值，水饱和度变化值和油气饱和度变化值；
10、所述未知向量为：
11、；
12、所述先验模型的储层属性估计为：
13、，
14、其中，为先验模型的压力变化值，为先验模型的水饱和度变化值，为先验模型的油气饱和度变化值；
15、为：
16、，
17、其中，为压力变化值的加权参数，为水饱和度变化值的加权参数，为油气饱和度变化值的加权参数，为压力变化值和水饱和度变化值的相关系数，为压力变化值和油气饱和度变化值的相关系数，为水饱和度变化值和油气饱和度变化值的相关系数；与为正交对角化的正交矩阵，为协方差矩阵的特征值构成的对角阵。
18、在一种可能的实现方式中，所述似然函数为：
19、，
20、其中，为似然函数，为实际采集角道集地震数据的负振幅映射属性，函数为先验模型模拟的角道集地震数据的负振幅映射属性，i为单位矩阵；
21、实际采集角道集地震数据的负振幅映射属性为：
22、，
23、其中，为实际采集的近偏移距角道集地震数据的负振幅映射属性，为实际采集的中偏移距角道集地震数据的负振幅映射属性，为实际采集的远偏移距角道集地震数据的负振幅映射属性；
24、函数为：
25、，
26、其中，为模拟的近偏移距角道集地震数据的负振幅映射属性，为模拟的中偏移距角道集地震数据的负振幅映射属性，为模拟的远偏移距角道集地震数据的负振幅映射属性；
27、，
28、其中，为近偏移距角道集地震数据的负振幅映射属性的加权参数，为中偏移距角道集地震数据的负振幅映射属性的加权参数，为远偏移距角道集地震数据的负振幅映射属性的加权参数，为实际采集的近偏移距角道集地震数据的负振幅映射属性和实际采集的中偏移距角道集地震数据的负振幅映射属性的相关系数，为实际采集的近偏移距角道集地震数据的负振幅映射属性和实际采集的远偏移距角道集地震数据的负振幅映射属性的相关系数，为实际采集的中偏移距角道集地震数据的负振幅映射属性和远偏移距角道集地震数据的负振幅映射属性的相关系数；与为正交对角化的正交矩阵，为协方差矩阵的特征值构成的对角阵。
29、在一种可能的实现方式中，所述后验概率为：
30、，
31、其中，为后验概率，为似然函数，为先验概率，d为实际采集角道集地震数据的负振幅映射属性，n为未知向量。
32、在一种可能的实现方式中，所述步骤s104包括：
33、步骤1：随机定义未知向量n的当前假设，计算所述当前假设对应的先验概率和似然函数，并根据所述当前假设对应的先验概率和似然函数计算所述当前假设对应的后验概率；对所述当前假设中的每个元素添加一个随机扰动，生成未知向量n的一个新假设，计算所述新假设对应的先验概率和似然函数，并根据所述新假设对应的先验概率和似然函数计算所述新假设对应的后验概率；如果新假设对应的后验概率高于当前假设对应的后验概率，则接受新假设,并将新假设作为当前假设，重复循环步骤1；
34、如果新假设对应的后验概率低于当前假设对应的后验概率，则计算比率r，r为新假设对应的后验概率与当前假设对应的后验概率之比；从0到1的均匀分布中抽取一个随机数u，如果u＜r，则接受新假设，将新假设作为当前假设，并重复循环步骤1；如果u≥r，保持当前假设不变，并重复循环步骤1；
35、步骤2：循环迭代一定次数后，检查连续迭代之间的后验概率的变化情况，如果连续迭代之间的后验概率的变化已经接近稳定，表示连续迭代之间的后验概率已经收敛，则终止循环。
36、步骤3：输出已经收敛的后验概率对应的未知向量n的一系列假设，对所述已经收敛的后验概率对应的未知向量n的一本文档来自技高网...
【技术保护点】
1.一种储层属性反演方法，其特征在于，所述方法包括：
2.根据权利要求1所述的一种储层属性反演方法，其特征在于，所述先验概率为：
3.根据权利要求2所述的一种储层属性反演方法，其特征在于，所述似然函数为：
4.根据权利要求3所述的一种储层属性反演方法，其特征在于，所述后验概率为：
5.根据权利要求4所述的一种储层属性反演方法，其特征在于，所述步骤S104包括：
6.根据权利要求2所述的一种储层属性反演方法，其特征在于，计算所述先验模型的压力变化值的标准差，作为压力变化值的加权参数；计算所述先验模型的水饱和度变化值的标准差，作为水饱和度变化值的加权参数；计算所述先验模型的油气饱和度变化值的标准差，作为油气饱和度变化值的加权参数；计算所述先验模型的压力变化值和水饱和度变化值的相关系数，作为压力变化值和水饱和度变化值的相关系数，计算所述先验模型的压力变化值和油气饱和度变化值的相关系数，作为压力变化值和油气饱和度变化值的相关系数，计算所述先验模型的水饱和度变化值和油气饱和度变化值的相关系数，作为水饱和度变化值和油气饱和度变化值的
7.根据权利要求3所述的一种储层属性反演方法，其特征在于，所述近偏移距角道集地震数据的负振幅映射属性的加权参数，所述中偏移距角道集地震数据的负振幅映射属性的加权参数和所述远偏移距角道集地震数据的负振幅映射属性的加权参数的计算公式为：
8.一种储层属性反演装置，其特征在于，所述装置包括：
9.一种电子设备，其特征在于，包括：
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质包括存储的程序，其中，在所述程序运行时控制所述计算机可读存储介质所在设备执行权利要求1至7中任意一项所述的方法。
...【技术特征摘要】
1.一种储层属性反演方法，其特征在于，所述方法包括：
2.根据权利要求1所述的一种储层属性反演方法，其特征在于，所述先验概率为：
3.根据权利要求2所述的一种储层属性反演方法，其特征在于，所述似然函数为：
4.根据权利要求3所述的一种储层属性反演方法，其特征在于，所述后验概率为：
5.根据权利要求4所述的一种储层属性反演方法，其特征在于，所述步骤s104包括：
6.根据权利要求2所述的一种储层属性反演方法，其特征在于，计算所述先验模型的压力变化值的标准差，作为压力变化值的加权参数；计算所述先验模型的水饱和度变化值的标准差，作为水饱和度变化值的加权参数；计算所述先验模型的油气饱和度变化值的标准差，作为油气饱和度变化值的加权参数；计算所述先验模型的压力变化值和水饱和度变化值的相关系数，作为压力变化值和水饱和度变化...
【专利技术属性】
技术研发人员：许艳，王奎峰，张太平，张红军，张龚，王圣雯，
申请(专利权)人：山东省地质科学研究院，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。