【技术实现步骤摘要】
本专利技术涉及人工智能，特别地涉及一种基于离散余弦变换的三维油藏历史拟合方法及系统。
技术介绍
1、在油藏数值模拟领域，准确的历史拟合对油藏动态预测有着至关重要的意义；然而，现有的常规商业油藏数值模拟软件在进行历史拟合时，主要采用整体或分区域调整渗透率的方式；这种方法存在诸多弊端，一方面，其调整过程主观性较强，且缺乏明确的分区依据；另一方面，多数情况下调整后的渗透率与孔隙度匹配不佳，无法正确反映原始地层参数间的关系，进而导致预测数据出现偏差。
2、在油藏历史拟合方面，现有的相关研究大多局限于二维场参数化和河道相，在处理三维复杂非均值多层系油藏时存在严重不足；即便最新的研究也仍停留在二维场参数化和河道相，实际应用价值有限，且三维复杂的常规砂岩、碳酸盐岩真实油藏的历史拟合。
技术实现思路
1、有鉴于此，本专利技术提出一种基于离散余弦变换的三维油藏历史拟合方法及系统，能够科学转换渗透率场、精准筛选关键频率域，实现数据压缩，大幅提升计算效率，有效调整参数，提升油藏历史拟合精度。
2、为实现上述目的，本专利技术提供如下技术方案：
3、一种基于离散余弦变换的三维油藏历史拟合方法，包括：
4、通过离散余弦变换将渗透率场从空间域转换为m1个频率域；
5、利用蒙特卡洛方法对所述m1个频率域进行扰动，生成n个先验地质模型；
6、对生成的n个先验地质模型进行数值模拟计算，得到每个先验地质模型对应的计算结果；
7、通过pa
8、对先验地质模型的模拟数据与实际数据进行分析，得到模拟数据的误差数据，判断所述误差数据是否超出预设误差阈值；
9、若超出，基于所述误差数据从所述m2个频率域中确定待调整参数，通过es - mda同化方法对所述待调整参数进行调整，直至生成的后验地质模型的误差数据未超出预设误差阈值；
10、若未超出，则判定三维油藏历史拟合已完成。
11、在上述技术方案的基础上，本专利技术还可以做如下改进：
12、进一步地，所述通过离散余弦变换将渗透率场从空间域转换为m1个频率域，包括：
13、通过公式（1）计算离散余弦变换后的系数；
14、公式（1）；
15、式中，为离散余弦变换后的系数，、、为归一化系数，、、分别为在x、y、z方向上的网格点数，为原始三维空间场中的函数值，i、j、k分别是在x、y、z方向上的空间索引。
16、进一步地，所述利用蒙特卡洛方法对所述m1个频率域进行扰动，生成n个先验地质模型，包括：
17、将每个频率域乘以对应的离散余弦变换后的系数然后加权求和得到n个先验地质模型，其中，每个频率域保留不同的渗透率场特征信息。
18、进一步地，所述通过pareto相关性分析方法对所述计算结果进行分析，得到与目标函数相关性最大的前m2个频率域，包括：
19、通过公式（2）计算两个频率域x和y之间的相关系数；
20、公式（2）；
21、式中，为两个频率域x和y之间的相关系数，为频率域x的第i个数据点，为频率域y的第i个数据点，为频率域x的平均值，为频率域y的平均值，为数据点的数量。
22、进一步地，所述基于所述误差数据从所述m2个频率域中确定待调整参数，包括：
23、通过公式（3）计算后验地质模型参数值；
24、公式（3）；
25、式中，为后验地质模型参数值，为先验地质模型参数值，为区域相关系数，为先验数据之间的协方差系数，为实际数据的方差系数，为在第j个位置的实际数据，为先验地质模型在第j个位置对应的模拟数据。
26、进一步地，所述基于所述误差数据从所述m2个频率域中确定待调整参数，通过es- mda同化方法对所述待调整参数进行调整，包括：
27、对误差数据进行多维度分解，从所述m2个频率域中确定m3个待调整频率域，将m3个待调整频率域对应的参数确定为待调整参数；
28、确定待调整参数的调整方向和参数调整幅度范围；
29、基于所述待调整参数的调整方向和参数调整幅度范围对所述待调整参数进行调整。
30、进一步地，所述通过es - mda同化方法对所述待调整参数进行调整，直至生成的后验地质模型的误差数据未超出预设误差阈值，包括：
31、通过es - mda同化方法对所述待调整参数进行调整后，对后验地质模型的模拟数据与实际数据进行分析，得到模拟数据对应的误差数据，判断所述误差数据是否超出预设误差阈值，不断迭代，直至所述误差数据小于预设误差阈值。
32、一种基于离散余弦变换的三维油藏历史拟合系统，包括：
33、离散余弦变换模块，用于通过离散余弦变换将渗透率场从空间域转换为m1个频率域；
34、先验地质模型生成模块，用于利用蒙特卡洛方法对所述m1个频率域进行扰动，生成n个先验地质模型；
35、数值模拟计算模块，用于对生成的n个先验地质模型进行数值模拟计算，得到每个先验地质模型对应的计算结果；
36、相关性分析模块，用于通过pareto相关性分析方法对所述计算结果进行分析，得到与目标函数相关性最大的前m2个频率域；
37、误差分析模块，用于对先验地质模型的模拟数据与实际数据进行分析，得到模拟数据的误差数据，判断所述误差数据是否超出预设误差阈值；
38、若超出，基于所述误差数据从所述m2个频率域中确定待调整参数，通过es - mda同化方法对所述待调整参数进行调整，直至生成的后验地质模型的误差数据未超出预设误差阈值；
39、若未超出，则判定三维油藏历史拟合已完成。
40、一种电子设备，包括存储器、处理器及存储在存储器上并在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如所述方法的步骤。
41、一种非暂态计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现所述方法的步骤。
42、本专利技术具有如下优点：
43、本专利技术中基于离散余弦变换的三维油藏历史拟合方法，通过离散余弦变换将渗透率场转换为频率域，能挖掘数据深层特征，如同对数据进行了条理化梳理，为后续精准分析提供良好基础。
44、本专利技术中基于离散余弦变换的三维油藏历史拟合方法，利用蒙特卡洛方法扰动频率域生成n个先验地质模型，充分考虑多种可能性，使先验模型更具多样性，更全面地涵盖油藏可能的地质状况，避免模型单一导致的拟合偏差。
45、本专利技术中基于离散余弦变换的三维油藏历史拟合方法，pareto相关性分析可筛选出与目标函数相关性最大的前m2个频率域，能精准定位对油藏关键指标（如产油量、产水量等）影响显著的因素，减少不必要的计算量，提高分析效率。
46本文档来自技高网...
【技术保护点】
1.一种基于离散余弦变换的三维油藏历史拟合方法，其特征在于，包括：
2.根据权利要求1所述的基于离散余弦变换的三维油藏历史拟合方法，其特征在于，所述通过离散余弦变换将渗透率场从空间域转换为M1个频率域，包括：
3.根据权利要求1所述的基于离散余弦变换的三维油藏历史拟合方法，其特征在于，所述利用蒙特卡洛方法对所述M1个频率域进行扰动，生成N个先验地质模型，包括：
4.根据权利要求1所述的基于离散余弦变换的三维油藏历史拟合方法，其特征在于，所述通过Pareto相关性分析方法对所述计算结果进行分析，得到与目标函数相关性最大的前M2个频率域，包括：
5.根据权利要求1所述的基于离散余弦变换的三维油藏历史拟合方法，其特征在于，所述基于所述误差数据从所述M2个频率域中确定待调整参数，包括：
6.根据权利要求1所述的基于离散余弦变换的三维油藏历史拟合方法，其特征在于，所述基于所述误差数据从所述M2个频率域中确定待调整参数，通过ES - MDA同化方法对所述待调整参数进行调整，包括：
7.根据权利要求6所述的基于离散余弦变换
8.一种基于离散余弦变换的三维油藏历史拟合的系统，其特征在于，包括：
9.一种电子设备，包括存储器、处理器及存储在存储器上并在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至7中的任一项所述的方法的步骤。
10.一种非暂态计算机可读介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至7中的任一项所述的方法的步骤。
...【技术特征摘要】
1.一种基于离散余弦变换的三维油藏历史拟合方法，其特征在于，包括：
2.根据权利要求1所述的基于离散余弦变换的三维油藏历史拟合方法，其特征在于，所述通过离散余弦变换将渗透率场从空间域转换为m1个频率域，包括：
3.根据权利要求1所述的基于离散余弦变换的三维油藏历史拟合方法，其特征在于，所述利用蒙特卡洛方法对所述m1个频率域进行扰动，生成n个先验地质模型，包括：
4.根据权利要求1所述的基于离散余弦变换的三维油藏历史拟合方法，其特征在于，所述通过pareto相关性分析方法对所述计算结果进行分析，得到与目标函数相关性最大的前m2个频率域，包括：
5.根据权利要求1所述的基于离散余弦变换的三维油藏历史拟合方法，其特征在于，所述基于所述误差数据从所述m2个频率域中确定待调整参数，包括：
6.根据权利要求1所述的基于离散余弦变...
【专利技术属性】
技术研发人员：杨迪，杨涛，龚斌，刘玄，苏瑾，
申请(专利权)人：特雷西能源科技股份有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。