CN117076956B - 缝洞油藏物理模型相似准数优化方法及装置 - Google Patents
缝洞油藏物理模型相似准数优化方法及装置 Download PDFInfo
- Publication number
- CN117076956B CN117076956B CN202311329733.5A CN202311329733A CN117076956B CN 117076956 B CN117076956 B CN 117076956B CN 202311329733 A CN202311329733 A CN 202311329733A CN 117076956 B CN117076956 B CN 117076956B
- Authority
- CN
- China
- Prior art keywords
- fracture
- oil reservoir
- determining
- physical model
- similarity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 77
- 238000005457 optimization Methods 0.000 title claims abstract description 44
- 238000011160 research Methods 0.000 claims abstract description 9
- 230000006870 function Effects 0.000 claims description 79
- 239000003921 oil Substances 0.000 claims description 62
- 239000010779 crude oil Substances 0.000 claims description 16
- 238000004519 manufacturing process Methods 0.000 claims description 15
- 238000002347 injection Methods 0.000 claims description 12
- 239000007924 injection Substances 0.000 claims description 12
- 230000001133 acceleration Effects 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 abstract description 10
- 238000006073 displacement reaction Methods 0.000 description 15
- 230000000694 effects Effects 0.000 description 15
- 230000008569 process Effects 0.000 description 12
- 239000000243 solution Substances 0.000 description 9
- 238000010586 diagram Methods 0.000 description 6
- 239000012530 fluid Substances 0.000 description 5
- 238000002360 preparation method Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011524 similarity measure Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000004141 dimensional analysis Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Hardware Design (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请公开了一种缝洞油藏物理模型相似准数优化方法及装置，该方法包括：根据缝洞油藏的特征确定物理模型的关键参数，并确定与缝洞油藏多相流动相关的主要相似准数；根据物理模型的研究要求确定主要相似准数的权重系数；根据主要相似准数与权重系数建立优化目标函数；利用优化目标函数确定关键参数的最终取值。解决了现有技术中的方程分析法虽然思路清晰，但对现象的数学描述要求比较严格，量纲分析法虽然不需要确定描述物理现象的数学方程，但相似准则的物理意义较差的问题。通过实施本申请中的方法能够根据所研究的问题较全面地考虑所有相似条件，进而能够较大程度地还原缝洞油藏的实际特征。
Description
技术领域
本申请涉及油藏物理模型与物模实验技术领域，尤其涉及一种缝洞油藏物理模型相似准数优化方法及装置。
背景技术
缝洞油藏储层具有空间随机分布、配置关系复杂、形状尺度变化多样等特点。在现有技术中采用物理模拟方法研究缝洞油藏中多相流的流动规律。物理模拟方法是根据相似原理设计物理模型（包括油藏和流体模型），能够在短时间内从一个观察尺度上测取所需的资料，再对资料进行分析归纳得到真实油藏尺度的流动规律，对油田生产具有重要的指导意义。
进行物理模拟实验时，为了使物理模型能够反映真实模型的特征，必须满足相似性条件。但同时考虑所有的相似性条件，将完整的油藏缝洞系统制作成一个较小的物理模型是无法实现的。因此，只能选取缝洞油藏的关键参数构建物理模型，现有技术中常用方程分析法和量纲分析法分析物理模型的关键参数。
方程分析法是在用准确的数学方程组和定解条件描述物理现象的情况下，根据物理方程的因次和谐性推导求得相似准数。方程分析法思路清晰，推导出的相似准则的物理意义明确。但对物理现象的数学描述要求较为严格，并且可能遗漏一些涉及到的相似准数。量纲分析法不需要描述物理现象的数学方程，只需要知道现象所涉及的物理量即可。用这种方法推导出的相似准则数群比较完整，但相似准则的物理意义没有方程分析法清楚。
发明内容
本申请实施例通过提供一种缝洞油藏物理模型相似准数优化方法，解决了现有技术中的方程分析法虽然思路清晰，但对现象的数学描述要求比较严格，量纲分析法虽然不需要确定描述物理现象的数学方程，但相似准则的物理意义较差的问题。实现了一种只需要知道现象所涉及的物理量，又能使相似度准则的物理意义清楚的相似准数优化方法。
第一方面，本申请实施例提供了一种缝洞油藏物理模型相似准数优化方法，包括：根据缝洞油藏的特征确定物理模型的关键参数，并确定与缝洞油藏多相流动相关的主要相似准数；根据所述物理模型的研究要求确定所述主要相似准数的权重系数；根据所述主要相似准数与所述权重系数建立优化目标函数；利用所述优化目标函数确定所述关键参数的最终取值。
结合第一方面，在第一种可能的实现方式中，所述关键参数包括溶洞直径、裂缝开度、原油粘度、原油密度、重力加速度、生产压差、产量与注入速度中的一个或多个。
结合第一方面，在第二种可能的实现方式中，所述优化目标函数如下：；式中，为所述优化目标函数，为第i个所述主要相似准数的所述
权重系数，为所述主要相似准数的偏差函数，为拟合参数向量，即所述关键参数，为观
测数据向量，即缝洞油藏的实际数据。
结合第一方面的第二种可能的实现方式，在第三种可能的实现方式中，所述偏差
函数，如下：；式中，为第i个所述主要相似准数的所述偏差函
数，为缝洞油藏的第i个所述主要相似准数，为所述物理模型的第i个所述主要
相似准数。
结合第一方面的第二种可能的实现方式，在第四种可能的实现方式中，所述利用所述优化目标函数确定所述关键参数的最终取值，包括：当所述优化目标函数满足预设条件时，确定所述关键参数的有效范围；根据所述有效范围确定所述关键参数的所述最终取值。
结合第一方面的第四种可能的实现方式，在第五种可能的实现方式中，所述确定
所述关键参数的所述最终取值的公式如下：；式中，为
所述优化目标函数的最小值，为第i个所述主要相似准数的所述权重系数，为所述主
要相似准数的所述偏差函数，为拟合参数向量，即所述关键参数，为观测数据向量，即缝
洞油藏的实际数据。
结合第一方面，在第六种可能的实现方式中，所述利用所述优化目标函数确定所述关键参数的最终取值后，还包括：根据所述关键参数的所述最终取值建立所述物理模型。
第二方面，本申请实施例提供了一种缝洞油藏物理模型相似准数优化装置，包括：确定模块：用于根据缝洞油藏的特征确定物理模型的关键参数，并确定与缝洞油藏多相流动相关的主要相似准数；权重系数模块，用于根据所述物理模型的研究要求确定所述主要相似准数的权重系数；建立模块，用于根据所述主要相似准数与所述权重系数建立优化目标函数；终值模块，用于利用所述优化目标函数确定所述关键参数的最终取值。
结合第二方面，在第一种可能的实现方式中，所述关键参数包括溶洞直径、裂缝开度、原油粘度、原油密度、重力加速度、生产压差、产量与注入速度中的一个或多个。
结合第二方面，在第二种可能的实现方式中，所述优化目标函数如下：；式中，为所述优化目标函数，为第i个所述主要相似准数的所述
权重系数，为所述主要相似准数的偏差函数，为拟合参数向量，即所述关键参数，为观
测数据向量，即缝洞油藏的实际数据。
结合第二方面的第二种可能的实现方式，在第三种可能的实现方式中，所述偏差
函数，如下：；式中，为第i个所述主要相似准数的所述偏差函
数，为缝洞油藏的第i个所述主要相似准数，为所述物理模型的第i个所述主要
相似准数。
结合第二方面的第二种可能的实现方式，在第四种可能的实现方式中，所述利用所述优化目标函数确定所述关键参数的最终取值，包括：当所述优化目标函数满足预设条件时，确定所述关键参数的有效范围；根据所述有效范围确定所述关键参数的所述最终取值。
结合第二方面的第四种可能的实现方式，在第五种可能的实现方式中，所述确定
所述关键参数的所述最终取值的公式如下：；式中，为
所述优化目标函数的最小值，为第i个所述主要相似准数的所述权重系数，为所述主
要相似准数的所述偏差函数，为拟合参数向量，即所述关键参数，为观测数据向量，即缝
洞油藏的实际数据。
结合第二方面，在第六种可能的实现方式中，所述利用所述优化目标函数确定所述关键参数的最终取值后，还包括：根据所述关键参数的所述最终取值建立所述物理模型。
第三方面，本申请实施例提供了一种设备，所述设备包括：处理器；用于存储处理器可执行指令的存储器；所述处理器执行所述可执行指令时，实现如第一方面或第一方面任一种可能实现的方式所述的方法。
第四方面，本申请实施例提供了一种非易失性计算机可读存储介质，所述非易失性计算机可读存储介质包括用于存储计算机程序或指令，当该计算机程序或指令被执行时，使如第一方面或第一方面任一种可能实现的方式所述的方法被实现。
本申请实施例中提供的一个或多个技术方案，至少具有如下技术效果或优点：
本申请实施例通过分析缝洞油藏的特征确定关键参数，以及确定与缝洞油藏多相流相关的主要相似准数，进而建立优化目标函数确定关键参数的最终取值，有效解决了现有技术中方程分析法对物理现象的数学描述要求比较严格，量纲分析法确定的相似准则的物理意义较差，进而实现了一种既不需要对物理现象进行数学描述，又能使相似准数的物理意义较为清楚的相似准数优化方法，能够根据所研究的问题较全面的考虑所有相似条件，进而能够较大程度地还原缝洞油藏的实际特征。
附图说明
为了更清楚地说明本申请实施例的技术方案，下面将对本申请实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的缝洞油藏物理模型相似准数优化方法的流程图；
图2为本申请实施例提供的缝洞油藏物理模型相似准数优化方法的装置结构图；
图3为本申请实施例提供的缝洞油藏的真实模型与物理模型的示意图；
图4为本申请实施例提供的不同物理模型与缝洞油藏的驱替-采出曲线图；
图5a至图5c为本申请实施例提供的缝洞油藏的实际驱替效果图与不同物理模型的驱替效果图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
以下对本申请实施例涉及的部分技术做出说明，以助于理解，应当将它们认为仅仅是示范性的。因此，本领域普通技术人员应当认识到，可以对这里描述的实施例做出各种改变和修改，而不会背离本申请的范围和精神。同样，为了清楚和简明，以下的描述中省略了部分对公知功能和结构的描述。
图1是本申请实施例提供的缝洞油藏物理模型相似准数优化方法的流程图，包括步骤101至步骤104。其中，图1仅为本申请实施例示出的一种执行顺序，并不代表缝洞油藏物理模型相似准数优化方法的唯一执行顺序，在可实现最终结果的情况下，图1所示出的步骤可以被并列或颠倒执行。
步骤101：根据缝洞油藏的特征确定物理模型的关键参数，并确定与缝洞油藏多相流动相关的主要相似准数。其中，关键参数包括溶洞直径、裂缝开度、原油粘度、原油密度、重力加速度、生产压差、产量与注入速度中的一个或多个。
在本申请实施例中，根据缝洞油藏开发的特征，确定建立物理模型涉及的关键参
数，包括溶洞直径、裂缝开度、原油粘度、原油密度、重力加速度、生产压差、产量与注入速度
八个关键参数。并根据定理方法，推导出了与缝洞油藏多相流的流动规律相关性最高的
七个主要相似准数，包括：1：缝洞导流能力比值，2：流度比（影响缝洞油藏的驱替效
果），3：生产压差与重力关系（影响缝洞油藏的驱替效果与物理模型的有效性），4：雷诺
数（惯性项与粘性项之比，影响流态的重要参数），5：注入量与流速之间的关系（影响缝洞
油藏的驱替效果与物理模型的有效性），6：邦德数（重力与流体表面张力的关系）与7：
弗劳德数（惯性力与重力之比）。
步骤102：根据物理模型的研究要求确定主要相似准数的权重系数。具体地，根据
物理模型的研究要求，可以修改不同主要相似准数的权重系数。在本申请实施例中，如果研
究缝洞油藏单元组合结构，则可以提高缝洞导流能力比值的权重系数。如果研究不同驱替
流体对缝洞油藏的影响，则可以提高流度比的权重系数。示例性地，在研究缝洞油藏的流动
规律时，分别确定主要相似准数的权重系数，将1的权重系数确定为0.2，2的权重系数
确定为0.2，3的权重系数确定为0.1，4的权重系数确定为0.3，5的权重系数确定为
0.1，6的权重系数确定为0，7的权重系数确定为0.1。
步骤103：根据主要相似准数与权重系数建立优化目标函数。具体地，根据主要相似准数，得到缝洞油藏与物理模型之间的方程组：
。为缝洞油藏的第i个主要相似准数，为物理模型的第
i个主要相似准数，并根据上述方程组定义偏差函数，如下：
。并且，，，式中，为第i个
主要相似准数的偏差函数，为缝洞油藏的第i个主要相似准数，为物理模型的
第i个主要相似准数，为拟合参数向量，即关键参数，为观测数据向量，即缝洞油藏的实
际数据。由于不可能同时满足所有主要相似准数，故偏差函数无法求得精确的解。为了减少
缝洞油藏与物理模型之间的相似准数的偏差，根据主要相似准数及其权重系数定义优化目
标函数，如下：
。式中，为优化目标函数，为第i个主要相似准
数的权重，为相似准数的偏差函数，为拟合参数向量，即关键参数，为观测数据向
量，即缝洞油藏的实际数据。进而将求取偏差函数转换为求取优化目标函数。
步骤104：利用优化目标函数确定关键参数的最终取值。具体地，当优化目标函数满足预设条件时，确定关键参数的有效范围。在本申请实施例中，将预设条件设置为当优化目标函数取最小值时。在有效范围中确定关键参数的最终取值。在本申请实施例中，受物理模型制备精度和/或实验器材等因素的限制，需要确定关键参数的有效范围。例如，受物理模型制备精度的限制，对于缝洞组合模型，裂缝的最小宽度为0.01m。又比如，受实验器材的材质限制，工况压力不能超过30MPa。故需要确定各个关键参数的有效范围。
在本申请实施例中，分别计算关键参数的具体取值的公式如下：
。式中，为优化目标函数的最小值，为第i个主要相似准数的权重，为第i个主要相似准数的偏差函数，为拟合参数向
量，即关键参数，为观测数据向量，即缝洞油藏的实际数据。
根据优化目标函数确定关键参数的最终取值后，根据关键参数的最终取值建立物理模型。如图3所示，是缝洞油藏的真实模型与物理模型的示意图。
本申请实施例通过确定关键参数与主要相似准数，能够避免描述物理现象的数学方程，并使相似准数的物理意义清楚。通过优化目标函数，能够降低缝洞油藏与物理模型的相似准数之间的偏差，进而使构建的物理模型能够较大程度地还原缝洞油藏的实际特征。如图4所示，为不同物理模型与缝洞油藏实际的驱替-采出曲线图。如图5a至图5c所示，图5a为缝洞油藏的实际驱替效果图，图5b为采用传统方法选取相似准数设计物理模型后的驱替效果图，图5c为采用本申请中的方法确定的相似准数设计物理模型后的驱替效果图。可以看出，使用本申请中的方法所建立的物理模型与缝洞油藏的实际效果较为接近。
虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。本实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照本实施例或者附图所示的方法顺序执行或者并行执行（例如并行处理器或者多线程处理的环境）。
如图2所示，本申请实施例还提供一种缝洞油藏物理模型相似准数优化装置200。该装置包括：确定模块201、权重系数模块202、建立模块203与终值模块204，具体如下。
确定模块201用于根据缝洞油藏的特征确定物理模型的关键参数，并确定与缝洞
油藏多相流动相关的主要相似准数。其中，关键参数包括溶洞直径、裂缝开度、原油粘度、原
油密度、重力加速度、生产压差、产量与注入速度中的一个或多个。确定模块201具体用于，
在本申请实施例中，根据缝洞油藏开发的特征，确定建立物理模型涉及的关键参数，包括溶
洞直径、裂缝开度、原油粘度、原油密度、重力加速度、生产压差、产量与注入速度八个关键
参数。并根据定理方法，推导出了与缝洞油藏多相流的流动规律相关性最高的七个主要
相似准数，包括：1：缝洞导流能力比值，2：流度比（影响缝洞油藏的驱替效果），3：生
产压差与重力关系（影响缝洞油藏的驱替效果与物理模型的有效性），4：雷诺数（惯性项
与粘性项之比，影响流态的重要参数），5：注入量与流速之间的关系（影响缝洞油藏的驱
替效果与物理模型的有效性），6：邦德数（重力与流体表面张力的关系）与7：弗劳德数
（惯性力与重力之比）。
权重系数模块202用于根据物理模型的研究要求确定主要相似准数的权重系数。
权重系数模块202具体用于，根据物理模型的研究要求，可以修改不同主要相似准数的权重
系数。在本申请实施例中，如果研究缝洞油藏单元组合结构，则可以提高缝洞导流能力比值
的权重系数。如果研究不同驱替流体对缝洞油藏的影响，则可以提高流度比的权重系数。示
例性地，在研究缝洞油藏的流动规律时，分别确定主要相似准数的权重系数，将1的权重
系数确定为0.2，2的权重系数确定为0.2，3的权重系数确定为0.1，4的权重系数确定
为0.3，5的权重系数确定为0.1，6的权重系数确定为0，7的权重系数确定为0.1。
建立模块203用于根据主要相似准数与权重系数建立优化目标函数。建立模块203具体用于，根据主要相似准数，得到缝洞油藏与物理模型之间的方程组：
。为缝洞油藏的第i个主要相似准数，为物理模型的第
i个主要相似准数，并根据上述方程组定义偏差函数，如下：
。并且，，，式中，为第i个
主要相似准数的偏差函数，为缝洞油藏的第i个主要相似准数，为物理模型的
第i个主要相似准数，为拟合参数向量，即关键参数，为观测数据向量，即缝洞油藏的实
际数据。由于不可能同时满足所有主要相似准数，故偏差函数无法求得精确的解。为了减少
缝洞油藏与物理模型之间的相似准数的偏差，根据主要相似准数及其权重系数定义优化目
标函数，如下：
。式中，为优化目标函数，为第i个主要相似准
数的权重，为相似准数的偏差函数，为拟合参数向量，即关键参数，为观测数据向
量，即缝洞油藏的实际数据。进而将求取偏差函数转换为求取优化目标函数。
终值模块204用于利用优化目标函数确定关键参数的最终取值。终值模块204具体用于，当优化目标函数满足预设条件时，确定关键参数的有效范围。在本申请实施例中，将预设条件设置为当优化目标函数取最小值时。在有效范围中确定关键参数的最终取值。在本申请实施例中，受物理模型制备精度和/或实验器材等因素的限制，需要确定关键参数的有效范围。例如，受物理模型制备精度的限制，对于缝洞组合模型，裂缝的最小宽度为0.01m。又比如，受实验器材的材质限制，工况压力不能超过30MPa。故需要确定各个关键参数的有效范围。
在本申请实施例中，分别计算关键参数的具体取值的公式如下：
。式中，为优化目标函数的最小值，为第i个主要相似准数的权重，为第i个主要相似准数的偏差函数，为拟合参数向
量，即关键参数，为观测数据向量，即缝洞油藏的实际数据。
根据优化目标函数确定关键参数的最终取值后，根据关键参数的最终取值建立物理模型。
本申请所述装置中的部分模块可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
上述申请实施例阐明的装置或模块，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述方便，描述以上装置时以功能分为各种模块分别描述。在实施本申请实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现。当然，也可以将实现某功能的模块由多个子模块或子单元组合实现。
本申请中所述的方法、装置或模块可以以计算机可读程序代码方式实现控制器按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该（微）处理器执行的计算机可读程序代码（例如软件或固件）的计算机可读介质、逻辑门、开关、专用集成电路（英文：Application Specific Integrated Circuit；简称：ASIC）、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320，存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本申请实施例还提供了一种设备，所述设备包括：处理器；用于存储处理器可执行指令的存储器；所述处理器执行所述可执行指令时，实现如本申请实施例所述的方法。
本申请实施例还提供了一种非易失性计算机可读存储介质，其上存储有计算机程序或指令，当该计算机程序或指令被执行时，使如本申请实施例中所述的方法被实现。
此外，在本发明的各个实施例中的各功能模块可以集成在一个处理模块中，也可以是各个模块单独存在，也可以两个或两个以上模块集成在一个模块中。
上述存储介质包括但不限于随机存取存储器（英文：Random Access Memory；简称：RAM）、只读存储器（英文：Read-Only Memory；简称：ROM）、缓存（英文：Cache）、硬盘（英文：Hard Disk Drive；简称：HDD）或者存储卡（英文：Memory Card）。所述存储器可以用于存储计算机程序指令。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的硬件的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，也可以通过数据迁移的实施过程中体现出来。该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备（可以是个人计算机，移动终端，服务器，或者网络设备等）执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施方式采用递进的方式描述，各个实施方式之间相同或相似的部分互相参见即可，每个实施方式重点说明的都是与其他实施方式的不同之处。本申请的全部或者部分可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、移动通信终端、多处理器系统、基于微处理器的系统、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
以上实施例仅用以说明本申请的技术方案，而非对本申请限制；尽管参照前述实施例对本申请进行了详细的说明，本领域普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请技术方案的范围。
Claims (8)
1.一种缝洞油藏物理模型相似准数优化方法，其特征在于，包括：
根据缝洞油藏的特征确定物理模型的关键参数，并确定与缝洞油藏多相流动相关的主要相似准数；其中，所述关键参数包括溶洞直径、裂缝开度、原油粘度、原油密度、重力加速度、生产压差、产量与注入速度中的一个或多个，所述主要相似准数，所述主要相似准数包括缝洞导流能力比值、流度比、生产压差与重力关系、雷诺数、注入量与流速之间的关系、邦德数以及弗劳德数；
根据所述物理模型的研究要求确定所述主要相似准数的权重系数；
根据所述主要相似准数与所述权重系数建立优化目标函数；其中，所述优化目标函数如下：
式中，f(x)为所述优化目标函数，wi为第i个所述主要相似准数的所述权重系数，dπi为所述主要相似准数的偏差函数，x为拟合参数向量，即所述关键参数，u为观测数据向量，即缝洞油藏的实际数据；
利用所述优化目标函数确定所述关键参数的最终取值。
2.根据权利要求1所述的方法，其特征在于，所述偏差函数，如下：
dπi＝(πi)实际-(πi)物模→0；式中，dπi为第i个所述主要相似准数的所述偏差函数，(πi)实际为缝洞油藏的第i个所述主要相似准数，(πi)物模为所述物理模型的第i个所述主要相似准数。
3.根据权利要求1所述的方法，其特征在于，所述利用所述优化目标函数确定所述关键参数的最终取值，包括：
当所述优化目标函数满足预设条件时，确定所述关键参数的有效范围；
根据所述有效范围确定所述关键参数的所述最终取值。
4.根据权利要求3所述的方法，其特征在于，所述确定所述关键参数的所述最终取值的公式如下：
式中，/>为所述优化目标函数的最小值，wi为第i个所述主要相似准数的所述权重系数，dπi为所述主要相似准数的所述偏差函数，x为拟合参数向量，即所述关键参数，u为观测数据向量，即缝洞油藏的实际数据。
5.根据权利要求1所述的方法，其特征在于，所述利用所述优化目标函数确定所述关键参数的最终取值，还包括：
根据所述关键参数的所述最终取值建立所述物理模型。
6.一种缝洞油藏物理模型相似准数优化装置，其特征在于，包括：
确定模块：用于根据缝洞油藏的特征确定物理模型的关键参数，并确定与缝洞油藏多相流动相关的主要相似准数；其中，所述关键参数包括溶洞直径、裂缝开度、原油粘度、原油密度、重力加速度、生产压差、产量与注入速度中的一个或多个，所述主要相似准数，所述主要相似准数包括缝洞导流能力比值、流度比、生产压差与重力关系、雷诺数、注入量与流速之间的关系、邦德数以及弗劳德数；
权重系数模块，用于根据所述物理模型的研究要求确定所述主要相似准数的权重系数；
建立模块，用于根据所述主要相似准数与所述权重系数建立优化目标函数；其中，所述优化目标函数如下：
式中，f(x)为所述优化目标函数，wi为第i个所述主要相似准数的所述权重系数，dπi为所述主要相似准数的偏差函数，x为拟合参数向量，即所述关键参数，u为观测数据向量，即缝洞油藏的实际数据；
终值模块，用于利用所述优化目标函数确定所述关键参数的最终取值。
7.一种用于执行缝洞油藏物理模型相似准数优化方法的设备，其特征在于，包括：
处理器；
用于存储处理器可执行指令的存储器；
所述处理器执行所述可执行指令时，实现如权利要求1至5中任意一项所述的方法。
8.一种非易失性计算机可读存储介质，其特征在于，包括用于存储计算机程序或指令，当该计算机程序或指令被执行时，使如权利要求1至5中任一项所述的方法被实现。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311329733.5A CN117076956B (zh) | 2023-10-16 | 2023-10-16 | 缝洞油藏物理模型相似准数优化方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311329733.5A CN117076956B (zh) | 2023-10-16 | 2023-10-16 | 缝洞油藏物理模型相似准数优化方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117076956A CN117076956A (zh) | 2023-11-17 |
| CN117076956B true CN117076956B (zh) | 2024-01-26 |
Family
ID=88702869
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311329733.5A Active CN117076956B (zh) | 2023-10-16 | 2023-10-16 | 缝洞油藏物理模型相似准数优化方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117076956B (zh) |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2066738C1 (ru) * | 1993-05-26 | 1996-09-20 | Устюжанин Александр Михайлович | Способ эксплуатации системы газлифтных скважин |
| CN101899972A (zh) * | 2010-06-30 | 2010-12-01 | 中国石油大学(北京) | 砂岩油藏注水开发可预测物理模型建立方法 |
| CN104616350A (zh) * | 2015-02-09 | 2015-05-13 | 西南石油大学 | 缝洞型碳酸盐油藏三维物理模型建立方法 |
| CN105117511A (zh) * | 2015-07-22 | 2015-12-02 | 西南石油大学 | 一种缝洞油藏井间连通通道及流动参数的表征方法 |
| CN106644362A (zh) * | 2016-11-24 | 2017-05-10 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种超音速空腔流动的马赫数敏感性分析方法 |
| CN108590639A (zh) * | 2018-03-30 | 2018-09-28 | 长江大学 | 三维油藏物理模型夹层设置方法 |
| CN110363325A (zh) * | 2019-05-06 | 2019-10-22 | 中国石油化工股份有限公司 | 复杂断块油藏多目标注采调控优化方法 |
| CN112861423A (zh) * | 2021-01-08 | 2021-05-28 | 中国石油大学(北京) | 数据驱动的注水油藏优化方法和系统 |
| CN113553746A (zh) * | 2021-08-18 | 2021-10-26 | 中国石油大学(北京) | 快速诊断缝洞油藏储层参数的方法和处理器 |
| CN113627640A (zh) * | 2020-05-08 | 2021-11-09 | 中国石油化工股份有限公司 | 一种针对缝洞型油藏油气井的产能试井预测方法及系统 |
| WO2021244521A1 (zh) * | 2020-06-04 | 2021-12-09 | 广州虎牙科技有限公司 | 模型训练、对象分类方法和装置、电子设备及存储介质 |
| CN115130274A (zh) * | 2022-04-02 | 2022-09-30 | 常州大学 | 一种稠油油藏相似准则研究方法 |
| CN115559714A (zh) * | 2022-08-31 | 2023-01-03 | 常州大学 | 一种多分支断溶体油藏缝洞参数的确定方法及系统 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2949586B1 (fr) * | 2009-08-26 | 2011-09-23 | Inst Francais Du Petrole | Methode pour exploiter un reservoir petrolier a partir d'un calage d'historique optimise |
| WO2016172827A1 (zh) * | 2015-04-27 | 2016-11-03 | 武汉武大卓越科技有限责任公司 | 一种逐步求精的路面裂缝检测方法 |
| EP3362640B1 (en) * | 2015-10-14 | 2023-07-26 | Landmark Graphics Corporation | History matching of hydrocarbon production from heterogenous reservoirs |
| US10983233B2 (en) * | 2019-03-12 | 2021-04-20 | Saudi Arabian Oil Company | Method for dynamic calibration and simultaneous closed-loop inversion of simulation models of fractured reservoirs |
| US20230267369A1 (en) * | 2022-01-31 | 2023-08-24 | Idaho State University | Generalized local adaptive fusion regression process based on physicochemical and physiochemical underlying hidden properties for quantitative analysis of molecular based spectroscopic data |
-
2023
- 2023-10-16 CN CN202311329733.5A patent/CN117076956B/zh active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2066738C1 (ru) * | 1993-05-26 | 1996-09-20 | Устюжанин Александр Михайлович | Способ эксплуатации системы газлифтных скважин |
| CN101899972A (zh) * | 2010-06-30 | 2010-12-01 | 中国石油大学(北京) | 砂岩油藏注水开发可预测物理模型建立方法 |
| CN104616350A (zh) * | 2015-02-09 | 2015-05-13 | 西南石油大学 | 缝洞型碳酸盐油藏三维物理模型建立方法 |
| CN105117511A (zh) * | 2015-07-22 | 2015-12-02 | 西南石油大学 | 一种缝洞油藏井间连通通道及流动参数的表征方法 |
| CN106644362A (zh) * | 2016-11-24 | 2017-05-10 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种超音速空腔流动的马赫数敏感性分析方法 |
| CN108590639A (zh) * | 2018-03-30 | 2018-09-28 | 长江大学 | 三维油藏物理模型夹层设置方法 |
| CN110363325A (zh) * | 2019-05-06 | 2019-10-22 | 中国石油化工股份有限公司 | 复杂断块油藏多目标注采调控优化方法 |
| CN113627640A (zh) * | 2020-05-08 | 2021-11-09 | 中国石油化工股份有限公司 | 一种针对缝洞型油藏油气井的产能试井预测方法及系统 |
| WO2021244521A1 (zh) * | 2020-06-04 | 2021-12-09 | 广州虎牙科技有限公司 | 模型训练、对象分类方法和装置、电子设备及存储介质 |
| CN112861423A (zh) * | 2021-01-08 | 2021-05-28 | 中国石油大学(北京) | 数据驱动的注水油藏优化方法和系统 |
| CN113553746A (zh) * | 2021-08-18 | 2021-10-26 | 中国石油大学(北京) | 快速诊断缝洞油藏储层参数的方法和处理器 |
| CN115130274A (zh) * | 2022-04-02 | 2022-09-30 | 常州大学 | 一种稠油油藏相似准则研究方法 |
| CN115559714A (zh) * | 2022-08-31 | 2023-01-03 | 常州大学 | 一种多分支断溶体油藏缝洞参数的确定方法及系统 |
Non-Patent Citations (3)
| Title |
|---|
| Bio-Inspired algorithms for the optimization of offshore oil production systems;Ian Nascimento Vieira 等;《Numerical Methods in Engineering》;第91卷(第10期);第1023-1044页 * |
| 油藏多相渗流体系物理模拟相似参数的数值分析方法;白玉湖 等;《2007年油藏地址建模与数值模拟技术应用研讨会》;第219-225页 * |
| 相似准则数值优化方法及其应用;周济福 等;《自然科学进展》;第15卷(第11期);第1297-1304页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117076956A (zh) | 2023-11-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU2017202259B2 (en) | Methods and systems for reservoir history matching for improved estimation of reservoir performance | |
| Ville et al. | Convected level set method for the numerical simulation of fluid buckling | |
| Khayyer et al. | Enhancement of stability and accuracy of the moving particle semi-implicit method | |
| Sodja | Turbulence models in CFD | |
| Tubbs et al. | GPU accelerated lattice Boltzmann model for shallow water flow and mass transport | |
| Liang | A structured but non‐uniform Cartesian grid‐based model for the shallow water equations | |
| Garzon et al. | Improvements of explicit crack surface representation and update within the generalized finite element method with application to three‐dimensional crack coalescence | |
| MX2013006155A (es) | Sistemas y metodos para reducir el tiempo de ejecucion de un modelo simulador de deposito. | |
| Azam et al. | A simple deterministic plastoelastohydrodynamic lubrication (PEHL) model in mixed lubrication | |
| Lai et al. | Numerical solution of the Saint-Venant equations by an efficient hybrid finite-volume/finite-difference method | |
| Kesserwani et al. | A conservative high‐order discontinuous Galerkin method for the shallow water equations with arbitrary topography | |
| Gu et al. | An effective way to control numerical instability of a nonordinary state‐based peridynamic elastic model | |
| CN117076956B (zh) | 缝洞油藏物理模型相似准数优化方法及装置 | |
| CN109635405B (zh) | 一种基于空间守恒的多相流动瞬态计算方法及系统 | |
| CN107451383B (zh) | 一种平面二维水沙数学模型初始床沙级配的率定方法 | |
| CN116734968A (zh) | 一种流量系数的校验方法、装置、设备及存储介质 | |
| Furnes | Flow in Pelton turbines | |
| CN112733415B (zh) | 一种薄壁弹性体边界的无网格处理方法、装置、终端设备及计算介质 | |
| Buachart et al. | A shallow water model by finite point method | |
| Yamazaki et al. | Comparison of uncertainty quantification approaches in a supersonic biplane airfoil problem | |
| JP2014241002A (ja) | 数値計算プログラム、数値計算方法及び情報処理装置 | |
| Augustyn et al. | Metrics-based auto scaling module for Amazon web services cloud platform | |
| Shobeyri | Simulating Free Surface Flows Using a New Incompressible SPH Model Improved by MPS Method | |
| Oertel | Numerical Modeling of free-surface flows in practical applications | |
| Hou | Development of a varying mesh scheme for finite element crack closure analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |