一种针对嵌入式离散裂缝模型局部网格加密方法及装置

本技术涉及油藏数值模拟，尤其涉及一种针对嵌入式离散裂缝模型局部网格加密方法及装置。

背景技术：

1、嵌入式离散裂缝模型（英文：embedded discrete fracture model，简称edfm）是一种应用于非常规裂缝性储层的流动模型，在该模型中模拟裂缝和基质之间的相互作用。该模型将裂缝网络直接嵌入基岩结构化网格系统中，在高精度保留了裂缝真实几何信息的同时避免了复杂的非结构化网格剖分过程。虽然仍需要计算裂缝与网格之间的几何信息，但相对于复杂的非结构化网格剖分过程，其计算复杂度会降低，从而能够提高计算效率。

2、目前的局部网格加密技术均是基于几何方法进行加密，即在裂缝附近选取部分区域，将原网格进行加密并重新计算裂缝和背景网格之间的传导率，裂缝影响区域和加密倍数等参数不好确定、容易形成“锯齿状”非连通区域且缺乏流动意义。

技术实现思路

1、在本技术实施例中，通过提供一种针对嵌入式离散裂缝模型局部网格加密方法，解决了目前的局部网格加密技术均是基于几何方法进行加密，即在裂缝附近选取部分区域，将原网格进行加密并重新计算裂缝和背景网格之间的传导率，裂缝影响区域和加密倍数等参数不好确定、容易形成“锯齿状”非连通区域且缺乏流动意义的问题。

2、第一方面，本技术实施例提供了一种针对嵌入式离散裂缝模型局部网格加密方法，该方法包括：对原始嵌入式离散裂缝模型进行粗略网格划分，获得粗略网格模型；获取所述粗略网格模型中的粗略网格的信息；在所述粗略网格模型中选取目标网格并对所述目标网格进行网格细化，获得局部细化网格模型；其中，所述目标网格为所述粗略网格模型中的裂缝所在的粗略网格；获取所述局部细化网格模型中的细化网格的信息；在所述局部细化网格模型中建立拟稳态压力场；在所述拟稳态压力场中根据划分公式和划分规则进行加密网格划分，获得局部加密网格模型；获取所述局部加密网格模型中加密网格的信息；将所述粗略网格模型中的所述目标网格的信息替换为所计算出的所述加密网格的信息，获得最终嵌入式离散裂缝模型。

3、结合第一方面，在一种可能的实现方式中，所述在所述局部细化网格模型中建立拟稳态压力场，包括：对所述局部细化网格模型中的边界和裂缝分别赋予不同的压力；根据公式对所述局部细化网格模型中每个细化网格列出平衡方程并联立求解，从而在所述局部细化网格模型中建立拟稳态压力场；其中，c为非零常量，为细化网格i的体积，为细化网格i的压力，为细化网格j的压力，为细化网格间的传导率。

4、结合第一方面，在一种可能的实现方式中，计算所述细化网格间的传导率的公式如下：；其中，为细化网格i的渗透率，为细化网格i间的传导因子，为细化网格j的渗透率，为细化网格j间的传导因子。

5、结合第一方面，在一种可能的实现方式中，所述划分规则包括：所述局部加密网格模型中的任一加密网格包含于所述局部细化网格模型，对于任意的两个加密网格的交集为空集，并且所有加密网格的并集为所述局部细化网格模型。

6、结合第一方面，在一种可能的实现方式中，所述划分公式为：；其中，为细化网格i对应的压力，为细化网格j对应的压力，为预设阈值。

7、结合第一方面，在一种可能的实现方式中，所述获取所述局部加密网格模型中加密网格的信息，包括：根据公式计算加密网格的体积；其中，为细化网格i的体积，m为细化网格i对应的加密网格；根据公式计算加密网格的深度；其中，为细化网格i的体积，为细化网格i对应的加密网格m的体积，为细化网格i的深度；根据公式计算加密网格的压力；其中，为细化网格i的体积，为细化网格i对应的加密网格m的体积，为细化网格i的压力；根据公式计算加密网格间的传导率；其中，为加密网格n到加密网格m的总流量，为加密网格n的平均压力，是加密网格m的平均压力；根据公式计算加密网格n到加密网格m的总流量；其中，为细化网格i到细化网格j的总流量，为细化网格间的传导率，为细化网格i的压力，为细化网格j的压力。

8、第二方面，本技术实施例提供了一种针对嵌入式离散裂缝模型局部网格加密装置，该装置包括：获得粗略网格模型模块，用于对原始嵌入式离散裂缝模型进行粗略网格划分，获得粗略网格模型；获取粗略网格的信息模块，用于获取所述粗略网格模型中的粗略网格的信息；获得局部细化网格模型模块，用于在所述粗略网格模型中选取目标网格并对所述目标网格进行网格细化，获得局部细化网格模型；其中，所述目标网格为所述粗略网格模型中的裂缝所在的粗略网格；获取细化网格的信息模块，用于获取所述局部细化网格模型中的细化网格的信息；建立模块，用于在所述局部细化网格模型中建立拟稳态压力场；获得局部加密网格模型模块，用于在所述拟稳态压力场中根据划分公式和划分规则进行加密网格划分，获得局部加密网格模型；获取加密网格的信息模块，用于获取所述局部加密网格模型中加密网格的信息；获得最终嵌入式离散裂缝模型模块，用于将所述粗略网格模型中的所述目标网格的信息替换为所计算出的所述加密网格的信息，获得最终嵌入式离散裂缝模型。

9、结合第二方面，在一种可能的实现方式中，所述在所述局部细化网格模型中建立拟稳态压力场，包括：对所述局部细化网格模型中的边界和裂缝分别赋予不同的压力；根据公式对所述局部细化网格模型中每个细化网格列出平衡方程并联立求解，从而在所述局部细化网格模型中建立拟稳态压力场；其中，c为非零常量，为细化网格i的体积，为细化网格i的压力，为细化网格j的压力，为细化网格间的传导率。

10、结合第二方面，在一种可能的实现方式中，计算所述细化网格间的传导率的公式如下：；其中，为细化网格i的渗透率，为细化网格i间的传导因子，为细化网格j的渗透率，为细化网格j间的传导因子。

11、结合第二方面，在一种可能的实现方式中，所述划分规则包括：所述局部加密网格模型中的任一加密网格包含于所述局部细化网格模型，对于任意的两个加密网格的交集为空集，并且所有加密网格的并集为所述局部细化网格模型。

12、结合第二方面，在一种可能的实现方式中，所述划分公式为：；其中，为细化网格i对应的压力，为细化网格j对应的压力，为预设阈值。

13、结合第二方面，在一种可能的实现方式中，所述获取所述局部加密网格模型中加密网格的信息，包括：根据公式计算加密网格的体积；其中，为细化网格i的体积，m为细化网格i对应的加密网格；根据公式计算加密网格的深度；其中，为细化网格i的体积，为细化网格i对应的加密网格m的体积，为细化网格i的深度；根据公式计算加密网格的压力；其中，为细化网格i的体积，为细化网格i对应的加密网格m的体积，为细化网格i的压力；根据公式计算加密网格间的传导率；其中，为加密网格n到加密网格m的总流量，为加密网格n的平均压力，是加密网格m的平均压力；根据公式计算加密网格n到加密网格m的总流量；其中，为细化网格i到细化网格j的总流量，为细化网格间的传导率，为细化网格i的压力，为细化网格j的压力。

14、第三方面，本技术实施例提供了一种针对嵌入式离散裂缝模型局部网格加密服务器，包括存储器和处理器；所述存储器用于存储计算机可执行指令；所述处理器用于执行所述计算机可执行指令，以实现第一方面或第一方面任一种可能的实现方式所述的方法。

15、第四方面，本技术实施例提供了一种计算机可读存储介质，所述计算机可读存储介质存储有可执行指令，计算机执行所述可执行指令时能够实现第一方面或第一方面任一种可能的实现方式所述的方法。

16、本技术实施例中提供的一个或多个技术方案，至少具有如下技术效果：

17、本技术实施例提供了一种针对嵌入式离散裂缝模型局部网格加密方法，该方法对原始嵌入式离散裂缝模型进行粗略网格划分，获得粗略网格模型，获取粗略网格模型中的粗略网格的信息，在粗略网格模型中选取目标网格并对目标网格进行网格细化，获得局部细化网格模型，其中，目标网格为粗略网格模型中的裂缝所在的粗略网格，获取局部细化网格模型中的细化网格的信息，在局部细化网格模型中建立拟稳态压力场，在拟稳态压力场中根据划分公式和划分规则进行加密网格划分，获得局部加密网格模型，获取局部加密网格模型中加密网格的信息，将粗略网格模型中的目标网格的信息替换为所计算出的加密网格的信息，获得最终嵌入式离散裂缝模型，解决了目前的局部网格加密技术均是基于几何方法进行加密，即在裂缝附近选取部分区域，将原网格进行加密并重新计算裂缝和背景网格之间的传导率，裂缝影响区域和加密倍数等参数不好确定、容易形成“锯齿状”非连通区域且缺乏流动意义的问题。