一种碳酸盐岩油藏酸压酸液滞留时间计算方法及装置制造方法及图纸

本申请公开了一种碳酸盐岩油藏酸压酸液滞留时间计算方法及装置，该方法包括：根据三维裂缝模型建立二维裂缝模型，将二维裂缝模型离散化划分成多个原网格；选择所需的多相流动模型和湍流模型；设置二维裂缝模型的边界条件和参数并进行仿真模拟，读取仿真模拟的结果；重新设定原网格的分辨率得到细网格，并建立粗化网格模型得到粗网格；将原网格中的酸液体积分数映射到粗网格中，并获取每个细网格位于粗网格中的酸液体积分数；获取粗网格中酸液总滞留时间

【技术实现步骤摘要】
一种碳酸盐岩油藏酸压酸液滞留时间计算方法及装置


[0001]本申请涉及石油开采
，尤其涉及一种碳酸盐岩油藏酸压酸液滞留时间计算方法及装置
。

技术介绍

[0002]碳酸盐岩蕴藏着全球
48%
的石油储量和
28%
的天然气储量
。
碳酸盐岩储层的高效开发对于满足国家能源需求具有重大意义
。
[0003]碳酸盐岩储层由于孔隙度和渗透率普遍较低，且非均质性强，开发效果普遍较差
。
近年来，以多级交替注入酸压技术为代表的深穿透酸压改造工艺是提高碳酸盐岩油藏开发效果的有效措施
。
其中，影响酸压效果的主要因素为酸蚀裂缝导流能力和酸蚀缝长，影响的情况与缝内酸液流动分布情况密切相关
。
[0004]但是，目前尚缺少可靠的交替酸压过程中酸液滞留时间的定量计算方法，制约了酸压效果评估与参数优化的水平，不能够定量表征酸液在裂缝中不均匀分布的强弱
。

技术实现思路

[0005]在本申请实施例中，通过提供一种碳酸盐岩油藏酸压酸液滞留时间计算方法及装置，解决了目前尚缺少可靠的交替酸压过程中酸液滞留时间的定量计算方法，制约了酸压效果评估与参数优化的水平，不能够定量表征酸液在裂缝中不均匀分布的强弱的问题
。
[0006]第一方面，本申请实施例提供了一种碳酸盐岩油藏酸压酸液滞留时间计算方法，该方法包括：根据三维裂缝模型建立二维裂缝模型，将所述二维裂缝模型离散化划分成多个原网格；选择所需的多相流动模型和湍流模型；设置所述二维裂缝模型的边界条件和参数并进行仿真模拟，读取所述仿真模拟的结果；重新设定所述原网格的分辨率得到细网格，并建立粗化网格模型得到粗网格；将所述原网格中的酸液体积分数映射到所述粗网格中，并获取每个所述细网格位于所述粗网格中的酸液体积分数；获取所述粗网格中酸液总滞留时间
。
[0007]结合第一方面，在一种可能的实现方式中，所述根据三维裂缝模型建立二维裂缝模型，包括：在所述三维裂缝模型中截取任意缝长与缝高剖面获得二维裂缝模型；在所述二维裂缝模型中设置酸液入口与酸液出口
。
[0008]结合第一方面，在一种可能的实现方式中，所述选择所需的多相流动模型和湍流模型，包括：选择流体体积模型为所述多相流动模型；根据公式和分别计算裂缝射孔数和有效射孔高度；其中，为裂缝高度，为射孔间距，为相位角度，为射孔孔径，射孔孔径小于或等于裂缝高度；根据公式计算酸液的注入速度；其中，为酸液注入排量，为裂缝宽度，为有效
射孔高度；根据公式计算雷诺数；其中，为流体密度，为粘性系数，为特征长度，为裂缝宽度，设定；设定临界雷诺数，判断所述雷诺数是否大于所述临界雷诺数；若判断结果为是，选择
k
‑
ε
模型为所述湍流模型；若判断结果为否，选择层流模型为所述湍流模型
。
[0009]结合第一方面，在一种可能的实现方式中，所述获取每个所述细网格位于所述粗网格中的酸液体积分数，包括：计算所述粗网格中所包含的每个所述细网格的体积；根据公式计算所述粗网格中所包含的每个所述细网格的权重；其中，为所计算的所述粗网格中所包含的每个所述细网格的体积；
G
为所述粗网格，
i
为所述细网格，为每个所述细网格中的酸液体积分数；根据公式计算每个所述细网格位于所述粗网格中的酸液体积分数
。
[0010]结合第一方面，在一种可能的实现方式中，所述获取所述粗网格中酸液总滞留时间，包括：根据公式计算所述粗网格中酸液总滞留时间；其中，为当前时刻
t
到下一时刻
t+1
的时间间距，为在
t
时刻所述粗网格中的酸液体积分数，为单位阶跃函数，为所设置的最低有效酸液体积分数阈值；根据公式判断时间段内，酸液是否在所述粗网格中滞留；若即判断结果为是，统计到所述粗网格的酸液总滞留时间中
。
[0011]第二方面，本申请实施例提供了一种碳酸盐岩油藏酸压酸液滞留时间计算装置，该装置包括：划分模块，用于根据三维裂缝模型建立二维裂缝模型，将所述二维裂缝模型离散化划分成多个原网格；选择模块，用于选择所需的多相流动模型和湍流模型；设置模块，用于设置所述二维裂缝模型的边界条件和参数并进行仿真模拟，读取所述仿真模拟的结果；重新设定模块，用于重新设定所述原网格的分辨率得到细网格，并建立粗化网格模型得到粗网格；映射模块，用于将所述原网格中的酸液体积分数映射到所述粗网格中，并获取每个所述细网格位于所述粗网格中的酸液体积分数；获取总滞留时间模块，用于获取所述粗网格中酸液总滞留时间
。
[0012]结合第二方面，在一种可能的实现方式中，所述划分模块具体用于：在所述三维裂缝模型中截取任意缝长与缝高剖面获得二维裂缝模型；在所述二维裂缝模型中设置酸液入口与酸液出口
。
[0013]结合第二方面，在一种可能的实现方式中，所述选择模块具体用于：选择流体体积模型为所述多相流动模型；根据公式和分别计算裂缝射孔数和有效射孔高度；其中，为裂缝高度，为射孔间距，为相位角度，为射孔孔径，射孔孔径小于或等于裂缝高度；根据公式计算酸液的注入速度；其中，为
酸液注入排量，为裂缝宽度，为有效射孔高度；根据公式计算雷诺数；其中，为流体密度，为粘性系数，为特征长度，为裂缝宽度，设定；设定临界雷诺数，判断所述雷诺数是否大于所述临界雷诺数；若判断结果为是，选择
k
‑
ε
模型为所述湍流模型；若判断结果为否，选择层流模型为所述湍流模型
。
[0014]结合第二方面，在一种可能的实现方式中，所述映射模块具体用于：计算所述粗网格中所包含的每个所述细网格的体积；根据公式计算所述粗网格中所包含的每个所述细网格的权重；其中，为所计算的所述粗网格中所包含的每个所述细网格的体积；
G
为所述粗网格，
i
为所述细网格，为每个所述细网格中的酸液体积分数；根据公式计算每个所述细网格位于所述粗网格中的酸液体积分数
。
[0015]结合第二方面，在一种可能的实现方式中，所述获取总滞留时间模块具体用于：根据公式计算所述粗网格中酸液总滞留时间；其中，为当前时刻
t
到下一时刻
t+1
的时间间距，为在
t
时刻所述粗网格中的酸液体积分数，为单位阶跃函数，为所设置的最低有效酸液体积分数阈值；根据公式判断时间段内，酸液是否在所述粗网格中滞留；若即判断结果为是，统计到所述粗网格的酸液总滞留时间中
。
[0016]第三方面，本申请实施例提供了一种碳酸盐岩油藏酸压酸液滞留时间计算服务器，包括存储器和处理器；所述存储器用于存储计算机可执行指令；所述处理器用于执行所述计算机可执行指令，以实现第一方面或第一方面任一种可能的实现方式所述的方法
。
[0017]第四方面，本申请实施例提本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种碳酸盐岩油藏酸压酸液滞留时间计算方法，其特征在于，包括：根据三维裂缝模型建立二维裂缝模型，将所述二维裂缝模型离散化划分成多个原网格；选择所需的多相流动模型和湍流模型；设置所述二维裂缝模型的边界条件和参数并进行仿真模拟，读取所述仿真模拟的结果；重新设定所述原网格的分辨率得到细网格，并建立粗化网格模型得到粗网格；将所述原网格中的酸液体积分数映射到所述粗网格中，并获取每个所述细网格位于所述粗网格中的酸液体积分数；获取所述粗网格中酸液总滞留时间
。2.
根据权利要求1所述的方法，其特征在于，所述根据三维裂缝模型建立二维裂缝模型，包括：在所述三维裂缝模型中截取任意缝长与缝高剖面获得二维裂缝模型；在所述二维裂缝模型中设置酸液入口与酸液出口
。3.
根据权利要求1所述的方法，其特征在于，所述选择所需的多相流动模型和湍流模型，包括：选择流体体积模型为所述多相流动模型；根据公式和分别计算裂缝射孔数和有效射孔高度；其中，为裂缝高度，为射孔间距，为相位角度，为射孔孔径，射孔孔径小于或等于裂缝高度；根据公式计算酸液的注入速度；其中，为酸液注入排量，为裂缝宽度，为有效射孔高度；根据公式计算雷诺数；其中，为流体密度，为粘性系数，为特征长度，为裂缝宽度，设定；设定临界雷诺数，判断所述雷诺数是否大于所述临界雷诺数；若判断结果为是，选择
k
‑
ε
模型为所述湍流模型；若判断结果为否，选择层流模型为所述湍流模型
。4.
根据权利要求1所述的方法，其特征在于，所述获取每个所述细网格位于所述粗网格中的酸液体积分数，包括：计算所述粗网格中所包含的每个所述细网格的体积；根据公式计算所述粗网格中所包含的每个所述细网格的权重；其中，为所计算的所述粗网格中所包含的每个所述细网格的体积；
G
为所...

【专利技术属性】
技术研发人员：李俊超，畅新鸽，
申请(专利权)人：西安石油大学，
类型：发明
国别省市：