(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211422234.6 (22)申请日 2022.11.15 (71)申请人 中科数智能源科技 （深圳） 有限公司 地址 518000 广东省深圳市龙岗区龙城街 道黄阁坑社区腾飞路9号创投大厦 1903 (72)发明人 龚斌　王宇　侯壮　方军龄　 (74)专利代理 机构 北京中和立达知识产权代理 有限公司 1 1756 专利代理师 孟姣 (51)Int.Cl. G06F 30/20(2020.01) G06F 111/10(2020.01) G06F 119/02(2020.01) (54)发明名称 一种裂缝油藏裂缝 敏感性分析方法 (57)摘要 本发明提供了一种裂缝油藏裂缝敏感性分 析方法， 包括： S1： 确定出需要分析的裂缝属性； S2： 基于油藏基础裂缝模型的裂缝属性参数， 设 置每个裂缝属性的波动范围； S3： 基于蒙特卡洛 方法在每个裂缝属性的波动范围内进行全变量 随机采样， 形成每个裂缝属性的参数组合； S4： 将 每个裂缝属性的参数组合输入油藏数值模拟器， 计算出每个裂缝属性的参数组合对应的产量； S5： 基于每个裂缝属性的参数组合对应的产量进 行裂缝属性的全局敏感性分析， 获得裂缝属性敏 感性分析结果； 用以同时考虑 多种裂缝属性对油 藏产量的敏感性， 进而实现对油藏裂缝属性进行 全局敏感性分析， 可为裂缝性油藏的动态分析， 开发方案调整提供 可靠依据。 权利要求书3页 说明书10页 附图2页 CN 115470664 A 2022.12.13 CN 115470664 A 1.一种裂缝油藏裂缝 敏感性分析 方法， 其特 征在于， 包括： S1： 确定出需要分析的裂缝属性； S2： 基于油藏基础裂缝模型的裂缝属性 参数， 设置每 个裂缝属性的波动范围； S3： 基于蒙特卡洛方法在每个裂缝属性的波动范围内进行全变量随机采样， 形成每个 裂缝属性的参数组合； S4： 将每个裂缝属性的参数组合输入油藏数值模拟器， 计算出每个裂缝属性的参数组 合对应的产量； S5： 基于每个裂缝属性的参数组合对应的产量进行裂缝属性的全局敏感性分析， 获得 裂缝属性敏感性分析 结果。 2.根据权利要求1所述的一种裂缝油藏裂缝敏感性分析方法， 其特征在于， 所述裂缝属 性包括： 裂缝密度、 裂缝 方位、 裂缝长度以及裂缝开度。 3.根据权利要求1所述的一种裂缝油藏裂缝敏感性分析方法， 其特征在于， S2： 基于油 藏基础裂缝模型的裂缝属性 参数， 设置每 个裂缝属性的波动范围， 包括： 基于确定性建模方法搭建出目标油田的大尺度裂缝模型； 基于随机建模方法搭建出目标油田的小尺度裂缝模型； 基于对目标油田的油藏基础裂缝模型的裂缝属性参数的先验认识， 设置每个裂缝属性 的波动范围； 其中， 油藏基础裂缝模型包括： 大尺度裂缝模型和小尺度裂缝模型。 4.根据权利要求1所述的一种裂缝油藏裂缝敏感性分析方法， 其特征在于， S3： 基于蒙 特卡洛方法在每个裂缝属性的波动范围内进 行全变量随机采样， 形成每个裂缝属性的参数 组合， 包括： 基于蒙特卡洛方法， 在每个裂缝属性的波动范围内进行全变量随机采样， 获得对应裂 缝属性的 个采样数据； 确定出裂缝属性的总个数 ； 将每个裂缝属性的前 个采样数据作为采样矩阵的前 列数据， 将每个裂缝属性的后 个采样数据作为采样矩阵的后 列数据， 获得 个裂缝属性的采样矩阵， 其中， 采样矩阵的 大小为 ； 将对应采样矩阵的前 列当作矩阵 ， 将对应采用矩阵的后 列当作矩阵 ； 将矩阵 中的第 列替换为矩阵 中的第 列， 获得矩阵 ； 将矩阵 和矩阵 以及矩阵 当作对应 裂缝属性的参数组合。 5.根据权利要求1所述的一种裂缝油藏裂缝敏感性分析方法， 其特征在于， S5： 基于每 个裂缝属性的参数组合对应的产量进 行裂缝属性的全局敏感性分析， 获得裂缝属性敏感性 分析结果， 包括： 基于每个裂缝属性参数组合的产量进行裂缝属性的全局敏感性分析， 获得全局敏感性权　利　要　求　书 1/3 页 2 CN 115470664 A 2分析结果； 基于全局敏感性分析结果绘制出裂缝属性敏感性统计图， 基于裂缝属性敏感性统计图 获得裂缝属性敏感性分析 结果。 6.根据权利要求5所述的一种裂缝油藏裂缝敏感性分析方法， 其特征在于， 基于每个裂 缝属性参数组合的产量进行裂缝属性的全局敏感性分析， 获得全局敏感性分析 结果， 包括： 基于每个裂缝属性 参数组合的产量， 计算出每 个裂缝属性对产量方差的贡献； 将每个裂缝属性对产量方差的贡献当作全局敏感性分析 结果。 7.根据权利要求6所述的一种裂缝油藏裂缝敏感性分析方法， 其特征在于， 基于每个裂 缝属性参数组合的产量， 计算出每 个裂缝属性对产量方差的贡献， 包括： 基于全局敏感性指数计算公式和每个裂缝属性参数组合的产量， 计算出对应裂缝属性 对产量方差的全局敏感性指数； 将全局敏感性指数作为对应 裂缝属性对产量方差的贡献。 8.根据权利要求7所述的一种裂缝油藏裂缝敏感性分析方法， 其特征在于， 基于全局敏 感性指数计算公式和每个裂缝属性参数组合的产量， 计算出对应裂缝属性对产量方差的全 局敏感性指数， 包括： 其中： 其中， 为第 个裂缝属性， 且第i个裂缝属性作为矩阵 中的第 列元素以及作为矩阵 中的第 列元素， 为第 个裂缝属性对产量方差的全局敏感性指数， 为样本总数， 为N中 的第 个样本， 为油藏数值模拟器的函数表达式， 为裂缝属性参数组合 中第 个样 本的产量， 为裂缝属性参数组合 中第 个样本的产量， 为矩阵 的方差， 为裂缝属性参数组合 对应的产量， 为裂缝属性参数组合 对应的产量， 为产 量为 条件下的裂缝属性 的方差， 为 的期望。 9.根据权利要求7所述的一种裂缝油藏裂缝敏感性分析方法， 其特征在于， 基于全局敏 感性分析结果绘制出裂缝属性敏感性统计图， 基于裂缝属性敏感性统计图获得裂缝属性敏 感性分析 结果， 包括： 基于全局敏感性分析结果中包含的每个裂缝属性对产量方差的全局敏感性指数， 绘制 出裂缝属性敏感性统计图， 其中， 裂缝属性敏感性统计图为柱状图形式； 基于裂缝属性敏感性统计图， 确定出影响产量的最敏感因素； 将最敏感因素和裂缝属性敏感性统计图作为裂缝属性敏感性分析 结果。权　利　要　求　书 2/3 页 3 CN 115470664 A 3