CN107704685B - 一种网格划分方法及装置 - Google Patents

一种网格划分方法及装置 Download PDF

Info

Publication number
CN107704685B
CN107704685B CN201710934317.6A CN201710934317A CN107704685B CN 107704685 B CN107704685 B CN 107704685B CN 201710934317 A CN201710934317 A CN 201710934317A CN 107704685 B CN107704685 B CN 107704685B
Authority
CN
China
Prior art keywords
grid
fracture
matrix
crack
grids
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
Application number
CN201710934317.6A
Other languages
English (en)
Other versions
CN107704685A (zh
Inventor
姜汉桥
糜利栋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201710934317.6A priority Critical patent/CN107704685B/zh
Publication of CN107704685A publication Critical patent/CN107704685A/zh
Application granted granted Critical
Publication of CN107704685B publication Critical patent/CN107704685B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Mining
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling

Abstract

本申请实施例公开了一种网格划分方法及装置。所述方法提供有目的油气藏的裂缝分布图像，其中，裂缝分布图像用于表征目的油气藏在水平二维面上的裂缝和基质的分布位置，所述方法包括：从裂缝分布图像中确定至少两个裂缝网格；其中，裂缝网格包括裂缝分布图像中相邻两个裂缝节点之间的裂缝部分，裂缝节点包括：裂缝分布图像中的裂缝交点与裂缝端点；确定裂缝分布图像中基质点到裂缝网格的距离，并将与基质点距离最近的裂缝网格作为基质点对应的裂缝网格；将对应同一个裂缝网格的基质点划分至同一个初始基质网格中；对初始基质网格进行网格划分处理，得到多个目标基质网格。本申请实施例提供的技术方案，可以降低网格数量，进一步提高计算效率。

Description

一种网格划分方法及装置
技术领域
本申请涉及油气藏数值模拟技术领域，特别涉及一种网格划分方法及装置。
背景技术
网格划分对油气藏数值模拟具有重要作用，所划分的网格数量与分布直接影响油气藏数值模拟的计算效率、计算精度和计算时间，而且网格对复杂求解区域的适应能力和逼近程度以及网格生产的自动化和智能化程度等，直接关系到油气藏数值模拟和数值结果的好坏。
在油气藏开发过程中，往往会通过人工压裂的方式在油气藏的基质天然裂缝的基础上增加人工裂缝，并将油气藏中的人工裂缝和基质天然裂缝，作为油气藏中油气的流动通道。针对这种裂缝性油气藏，目前通常采用直角平分(Perpendicular Bisection,PEBI)网格划分方法进行裂缝网格划分，再对裂缝网格划分后的油气藏进行数值模拟并分析油气藏中油气的流动情况。然而，当油气藏中的裂缝非均匀分布时，在裂缝交点和裂缝附近会产生大量尺寸较小的网格，导致网格数量较大，进而导致计算效率较低，计算时间较长。
发明内容
本申请实施例的目的是提供一种网格划分方法及装置，以降低网格数量，进一步提高计算效率，缩短计算时间。
为解决上述技术问题，本申请实施例提供一种网格划分方法及装置是这样实现的：
一种网格划分方法，提供有目的油气藏的裂缝分布图像，其中，所述裂缝分布图像用于表征所述目的油气藏在水平二维面上的裂缝和基质的分布位置，所述方法包括：
从所述裂缝分布图像中确定至少两个裂缝网格；其中，所述裂缝网格包括所述裂缝分布图像中相邻两个裂缝节点之间的裂缝部分，所述裂缝节点包括：所述裂缝分布图像中的裂缝交点与裂缝端点；
确定所述裂缝分布图像中基质点到所述裂缝网格的距离，并将与所述基质点距离最近的裂缝网格作为所述基质点对应的裂缝网格；
将对应同一个裂缝网格的基质点划分至同一个初始基质网格中；
对所述初始基质网格进行网格划分处理，得到多个目标基质网格。
优选方案中，所述基质点到所述裂缝网格的距离表示所述基质点与所述裂缝网格内指定像素点之间的距离，所述指定像素点表示所述裂缝网格内的与所述基质点的距离最短的像素点。
优选方案中，采用下述表达式表征所述裂缝网格的数量：
Figure BDA0001429459360000021
其中，Nfracture-grid表示所述裂缝网格的数量，Nnode-i表示所述裂缝分布图像中第i条裂缝对应的裂缝交点的数量，fracture-n表示所述裂缝分布图像中裂缝的数量。
优选方案中，所述对所述初始基质网格进行网格划分处理，得到多个目标基质网格，包括：
将所述初始基质网格转换为矩形初始基质网格；其中，所述矩形初始基质网格与所述初始基质网格一一对应，所述矩形初始基质网格的面积与所述初始基质网格的面积相同，且所述矩形初始基质网格的指定边长的长度与所述初始基质网格对应的裂缝网格的裂缝长度相同；
对所述矩形初始基质网格进行对数加密处理，以使得将所述矩形初始基质网格划分为多个目标基质网格；其中，所述目标基质网格的一边与所述指定边长平行，且长度相同；所述多个目标基质网格中相邻两个目标基质网格中第一基质网格的面积为第二基质网格的面积的指定倍数；所述第一基质网格表示所述相邻两个目标基质网格中距离与所述初始基质网格对应的裂缝网格较远的目标基质网格；所述第二基质网格表示所述相邻两个目标基质网格中距离与所述初始基质网格对应的裂缝网格较近的目标基质网格。
优选方案中，所述指定倍数的取值范围为2～5。
优选方案中，当一个所述初始基质网格划分为多个目标基质网格时，所述多个目标基质网格中目标基质网格的数量的取值范围为小于或等于10、且大于或等于2的整数。
优选方案中，在从所述裂缝分布图像中确定至少两个裂缝网格之前，对所述裂缝分布图像进行二值化处理，得到二值化图像；
相应地，所述从所述裂缝分布图像中确定至少两个裂缝网格，包括：从所述二值化图像中确定至少两个裂缝网格；其中，所述二值化图像包括：裂缝和基质。
一种网格划分装置，所述装置提供目的油气藏的裂缝分布图像，其中，所述裂缝分布图像用于表征所述目的油气藏在水平二维面上的裂缝和基质的分布位置，所述装置包括：裂缝网格确定模块、基质点对应的裂缝网格确定模块、初始基质网格确定模块和目标基质网格确定模块；其中，
所述裂缝网格确定模块，用于从所述裂缝分布图像中确定至少两个裂缝网格；其中，所述裂缝网格包括所述裂缝分布图像中相邻两个裂缝节点之间的裂缝部分，所述裂缝节点包括：所述裂缝分布图像中的裂缝交点与裂缝端点；
所述基质点对应的裂缝网格确定模块，用于确定所述裂缝分布图像中基质点到所述裂缝网格的距离，并将与所述基质点距离最近的裂缝网格作为所述基质点对应的裂缝网格；
所述初始基质网格确定模块，用于将对应同一个裂缝网格的基质点划分至同一个初始基质网格中；
所述目标基质网格确定模块，用于对所述初始基质网格进行网格划分处理，得到多个目标基质网格。
优选方案中，所述基质点到所述裂缝网格的距离表示所述基质点与所述裂缝网格内指定像素点之间的距离，所述指定像素点表示所述裂缝网格内的与所述基质点的距离最短的像素点。
优选方案中，所述目标基质网格确定模块包括：基质网格转换模块和对数加密处理模块；其中，
所述基质网格转换模块，用于将所述初始基质网格转换为矩形初始基质网格；其中，所述矩形初始基质网格与所述初始基质网格一一对应，所述矩形初始基质网格的面积与所述初始基质网格的面积相同，且所述矩形初始基质网格的指定边长的长度与所述初始基质网格对应的裂缝网格的裂缝长度相同；
所述对数加密处理模块，用于对所述矩形初始基质网格进行对数加密处理，以使得将所述矩形初始基质网格划分为多个目标基质网格；其中，所述目标基质网格的一边与所述指定边长平行，且长度相同；所述多个目标基质网格中相邻两个目标基质网格中第一基质网格的面积为第二基质网格的面积的指定倍数；所述第一基质网格表示所述相邻两个目标基质网格中距离与所述初始基质网格对应的裂缝网格较远的目标基质网格；所述第二基质网格表示所述相邻两个目标基质网格中距离与所述初始基质网格对应的裂缝网格较近的目标基质网格。
本申请实施例提供了一种网格划分方法及装置，可以从所述裂缝分布图像中确定至少两个裂缝网格；其中，所述裂缝网格包括所述裂缝分布图像中相邻两个裂缝节点之间的裂缝部分，所述裂缝节点包括：所述裂缝分布图像中的裂缝交点与裂缝端点；可以确定所述裂缝分布图像中基质点到所述裂缝网格的距离，并将与所述基质点距离最近的裂缝网格作为所述基质点对应的裂缝网格；可以将对应同一个裂缝网格的基质点划分至同一个初始基质网格中；可以对所述初始基质网格进行网格划分处理，得到多个目标基质网格。本申请方法通过所述裂缝节点对所述裂缝分布图像中的裂缝进行离散划分，所得到的裂缝网格的数量较少，从而可以减少最终划分的目标基质网格的数量，进一步提高计算效率。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本申请一种网格划分方法实施例的流程图；
图2是本申请实施例中确定裂缝网格的示意图；
图3是本申请实施例中基质像素点至裂缝网格之间距离的示意图；
图4是本申请实施例中网格划分处理的示意图；
图5是本申请实施例中裂缝分布图像和网格划分结果的示意图；
图6是本申请实施例中CMG软件和本申请的网格划分方法分别对应的模型生产井的日产量动态特征曲线；
图7本申请实施例中CMG软件和本申请的网格划分方法分别对应的模型生产井的累产量动态特征曲线；
图8是本申请网格划分装置实施例的组成结构图；
图9是本申请网格划分装置实施例中目标基质网格确定模块的组成结构图。
具体实施方式
本申请实施例提供一种网格划分方法及装置。
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
本申请实施例提供了一种网格划分方法。所述网格划分方法提供有目的油气藏的裂缝分布图像。
在本实施方式中，所述目的油气藏可以是页岩气藏。
在本实施方式中，所述裂缝分布图像用于表征所述目的油气藏在水平二维面上的裂缝和基质的分布位置。所述裂缝分布图像中可以包括裂缝和基质。
图1是本申请一种网格划分方法实施例的流程图。如图1所示，所述网格划分方法，包括以下步骤。
步骤S101：从所述裂缝分布图像中确定至少两个裂缝网格；其中，所述裂缝网格包括所述裂缝分布图像中相邻两个裂缝节点之间的裂缝部分，所述裂缝节点包括：所述裂缝分布图像中的裂缝交点与裂缝端点。
在本实施方式中，可以采用下述表达式表征所述裂缝网格的数量：
Figure BDA0001429459360000051
其中，Nfracture-grid表示所述裂缝网格的数量，Nnode-i表示所述裂缝分布图像中第i条裂缝对应的裂缝交点的数量，fracture-n表示所述裂缝分布图像中裂缝的数量。
例如，图2是本申请实施例中确定裂缝网格的示意图。图2中(a)为裂缝的示意图，其中，彼此相交的黑色线条表示裂缝。图2中(b)为从图2中(a)中确定裂缝网络的示意图，如图2所示，可以将图2中(b)中虚线箭头所指的相邻的节点(即所述裂缝节点)之间的裂缝部分作为所述裂缝网格。例如，可以相邻的裂缝端点和裂缝交点，或相邻的裂缝交点和裂缝交点之间的裂缝部分作为所述裂缝网格，即图2中(b)中虚线箭头所指的裂缝网格。
步骤S102：确定所述裂缝分布图像中基质点到所述裂缝网格的距离，并将与所述基质点距离最近的裂缝网格作为所述基质点对应的裂缝网格。
在本实施方式中，可以确定所述裂缝分布图像中基质点到所述裂缝网格的距离，并将与所述基质点距离最近的裂缝网格作为所述基质点对应的裂缝网格。其中，所述基质点到所述裂缝网格的距离表示所述基质点与所述裂缝网格内指定像素点之间的距离，所述指定像素点表示所述裂缝网格内的与所述基质点的距离最短的像素点。所述基质点表示所述裂缝分布图像中的基质对应的像素点。
例如，图3是本申请实施例中基质像素点至裂缝网格之间距离的示意图。图3中(a)为所述裂缝分布图像中基质对应的像素点到裂缝网格的距离的示意图，其中，像素点j表示所述裂缝分布图像中基质对应的像素点，裂缝网格i表示所述裂缝网格。图3中(b)为所述裂缝分布图像中基质对应的像素点到多个裂缝网格的距离的示意图，其中，黑色线条1、2、3和4表示所述裂缝网格。
步骤S103：将对应同一个裂缝网格的基质点划分至同一个初始基质网格中。
在本实施方式中，可以将对应同一个裂缝网格的基质点划分至同一个初始基质网格中。
步骤S104：对所述初始基质网格进行网格划分处理，得到多个目标基质网格。
在本实施方式中，对所述初始基质网格进行网格划分处理，得到多个目标基质网格，具体可以包括，可以将所述初始基质网格转换为矩形初始基质网格。其中，所述矩形初始基质网格与所述初始基质网格一一对应，所述矩形初始基质网格的面积与所述初始基质网格的面积相同，且所述矩形初始基质网格的指定边长的长度与所述初始基质网格对应的裂缝网格的裂缝长度相同。可以对所述矩形初始基质网格进行对数加密处理，以使得将所述矩形初始基质网格划分为多个目标基质网格；其中，所述目标基质网格的一边与所述指定边长平行，且长度相同；所述多个目标基质网格中相邻两个目标基质网格中第一基质网格的面积为第二基质网格的面积的指定倍数；所述第一基质网格表示所述相邻两个目标基质网格中距离与所述初始基质网格对应的裂缝网格较远的目标基质网格；所述第二基质网格表示所述相邻两个目标基质网格中距离与所述初始基质网格对应的裂缝网格较近的目标基质网格。
例如，图4是本申请实施例中网格划分处理的示意图。图4中(a)、(b)和(c)分别为所述初始基质网格、矩形初始基质网格和对数加密处理后的基质网格的示意图。图4中(a)中的虚线框内部分表示与所述初始基质网格对应的矩形初始基质网格，即图4中(b)所示的矩形初始基质网格。图4中(c)中的虚线表示将所述矩形初始基质网格划分为多个目标基质网格时的划分线。
在本实施方式中，所述指定倍数的取值范围可以为2～5。
在本实施方式中，当一个所述初始基质网格划分为多个目标基质网格时，所述多个目标基质网格中目标基质网格的数量的取值范围可以为小于或等于10、且大于或等于2的整数。
在另一种实施方式中，在从所述裂缝分布图像中确定至少两个裂缝网格之前，可以对所述裂缝分布图像进行二值化处理，可以得到二值化图像。
相应地，所述从所述裂缝分布图像中确定至少两个裂缝网格，具体可以包括，从所述二值化图像中确定至少两个裂缝网格。其中，所述二值化图像可以包括：裂缝和基质。
所述确定所述裂缝分布图像中基质点到所述裂缝网格的距离，并将与所述基质点距离最近的裂缝网格作为所述基质点对应的裂缝网格，具体可以包括，可以确定所述二值化图像中基质点到所述裂缝网格的距离，并将与所述基质点距离最近的裂缝网格作为所述基质点对应的裂缝网格。
例如，图5是本申请实施例中裂缝分布图像和网格划分结果的示意图。图5中(a)为二值化处理后的裂缝分布图像，即所述二值化图像。图5中(a)中的黑色粗线表示水平井，黑色细线表示裂缝。图5中(b)为所述初始基质网格的网格划分结果。图5中(b)中不同颜色区块表示不同的初始基质网格的编号。
在一个具体实施场景中，图6是本申请实施例中CMG软件和本申请的网格划分方法分别对应的模型生产井的日产量动态特征曲线。图7本申请实施例中CMG软件和本申请的网格划分方法分别对应的模型生产井的累产量动态特征曲线。可以分别基于油气藏数值模拟CMG软件和本申请的网格划分方法所得的网格划分结果对某页岩气藏进行数值模拟，可以得到CMG软件和本申请的网格划分方法分别对应的模型生产井动态特征曲线。如图6和图7所示，基于本申请的网格划分方法对应的模型生产井动态特征曲线与CMG软件对应的模型生产井动态特征曲线基本一致，表明本申请的网格划分方法比较合理，准确度较高。同时，在相同准确度的情况下，相比CMG软件中的网格划分方法，本申请的网格划分方法所划分的裂缝网格和目标基质网格数量较少，可以进一步提高油气藏数值模拟的计算效率，缩短计算时间。
所述网格划分方法实施例，可以从所述裂缝分布图像中确定至少两个裂缝网格；其中，所述裂缝网格包括所述裂缝分布图像中相邻两个裂缝节点之间的裂缝部分，所述裂缝节点包括：所述裂缝分布图像中的裂缝交点与裂缝端点；可以确定所述裂缝分布图像中基质点到所述裂缝网格的距离，并将与所述基质点距离最近的裂缝网格作为所述基质点对应的裂缝网格；可以将对应同一个裂缝网格的基质点划分至同一个初始基质网格中；可以对所述初始基质网格进行网格划分处理，得到多个目标基质网格。本申请方法通过所述裂缝节点对所述裂缝分布图像中的裂缝进行离散划分，所得到的裂缝网格的数量较少，从而可以减少最终划分的目标基质网格的数量，进一步提高计算效率。
图8是本申请网格划分装置实施例的组成结构图。所述网格划分装置提供目的油气藏的裂缝分布图像，其中，所述裂缝分布图像用于表征所述目的油气藏在水平二维面上的裂缝和基质的分布位置。如图8所示，所述网格划分装置可以包括：裂缝网格确定模块100、基质点对应的裂缝网格确定模块200、初始基质网格确定模块300和目标基质网格确定模块400。
所述裂缝网格确定模块100，可以用于从所述裂缝分布图像中确定至少两个裂缝网格；其中，所述裂缝网格包括所述裂缝分布图像中相邻两个裂缝节点之间的裂缝部分，所述裂缝节点包括：所述裂缝分布图像中的裂缝交点与裂缝端点。
所述基质点对应的裂缝网格确定模块200，可以用于确定所述裂缝分布图像中基质点到所述裂缝网格的距离，并将与所述基质点距离最近的裂缝网格作为所述基质点对应的裂缝网格。其中，所述基质点到所述裂缝网格的距离表示所述基质点与所述裂缝网格内指定像素点之间的距离，所述指定像素点表示所述裂缝网格内的与所述基质点的距离最短的像素点。
所述初始基质网格确定模块300，可以用于将对应同一个裂缝网格的基质点划分至同一个初始基质网格中。
所述目标基质网格确定模块400，可以用于对所述初始基质网格进行网格划分处理，得到多个目标基质网格。
图9是本申请网格划分装置实施例中目标基质网格确定模块的组成结构图。如图9所示，图8中的目标基质网格确定模块400可以包括：基质网格转换模块410和对数加密处理模块420。
所述基质网格转换模块410，可以用于将所述初始基质网格转换为矩形初始基质网格；其中，所述矩形初始基质网格与所述初始基质网格一一对应，所述矩形初始基质网格的面积与所述初始基质网格的面积相同，且所述矩形初始基质网格的指定边长的长度与所述初始基质网格对应的裂缝网格的裂缝长度相同。
所述对数加密处理模块420，可以用于对所述矩形初始基质网格进行对数加密处理，以使得将所述矩形初始基质网格划分为多个目标基质网格；其中，所述目标基质网格的一边与所述指定边长平行，且长度相同；所述多个目标基质网格中相邻两个目标基质网格中第一基质网格的面积为第二基质网格的面积的指定倍数；所述第一基质网格表示所述相邻两个目标基质网格中距离与所述初始基质网格对应的裂缝网格较远的目标基质网格；所述第二基质网格表示所述相邻两个目标基质网格中距离与所述初始基质网格对应的裂缝网格较近的目标基质网格。
所述网格划分装置实施例与所述网格划分方法实施例相对应，可以实现网格划分方法实施例的技术方案，并取得方法实施例的技术效果。
在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(Hardware Description Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等，目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的装置、模块，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中，内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括短暂电脑可读媒体(transitory media)，如调制的数据信号和载波。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于装置实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (8)

1.一种网格划分方法，其特征在于，提供有目的油气藏的裂缝分布图像，其中，所述裂缝分布图像用于表征所述目的油气藏在水平二维面上的裂缝和基质的分布位置，所述方法包括：
从所述裂缝分布图像中确定至少两个裂缝网格；其中，所述裂缝网格包括所述裂缝分布图像中相邻两个裂缝节点之间的裂缝部分，所述裂缝节点包括：所述裂缝分布图像中的裂缝交点与裂缝端点；
确定所述裂缝分布图像中基质点到所述裂缝网格的距离，并将与所述基质点距离最近的裂缝网格作为所述基质点对应的裂缝网格；
将对应同一个裂缝网格的基质点划分至同一个初始基质网格中；
对所述初始基质网格进行网格划分处理，得到多个目标基质网格，其中，包括：将所述初始基质网格转换为矩形初始基质网格；其中，所述矩形初始基质网格与所述初始基质网格一一对应，所述矩形初始基质网格的面积与所述初始基质网格的面积相同，且所述矩形初始基质网格的指定边长的长度与所述初始基质网格对应的裂缝网格的裂缝长度相同；对所述矩形初始基质网格进行对数加密处理，以使得将所述矩形初始基质网格划分为多个目标基质网格；其中，所述目标基质网格的一边与所述指定边长平行，且长度相同；所述多个目标基质网格中相邻两个目标基质网格中第一基质网格的面积为第二基质网格的面积的指定倍数；所述第一基质网格表示所述相邻两个目标基质网格中距离与所述初始基质网格对应的裂缝网格较远的目标基质网格；所述第二基质网格表示所述相邻两个目标基质网格中距离与所述初始基质网格对应的裂缝网格较近的目标基质网格。
2.根据权利要求1所述的一种网格划分方法，其特征在于，所述基质点到所述裂缝网格的距离表示所述基质点与所述裂缝网格内指定像素点之间的距离，所述指定像素点表示所述裂缝网格内的与所述基质点的距离最短的像素点。
3.根据权利要求1所述的一种网格划分方法，其特征在于，采用下述表达式表征所述裂缝网格的数量：
Figure FDA0002250163800000011
其中，Nfracture-grid表示所述裂缝网格的数量，Nnode-i表示所述裂缝分布图像中第i条裂缝对应的裂缝交点的数量，fracture-n表示所述裂缝分布图像中裂缝的数量。
4.根据权利要求1所述的一种网格划分方法，其特征在于，所述指定倍数的取值范围为2～5。
5.根据权利要求1所述的一种网格划分方法，其特征在于，当一个所述初始基质网格划分为多个目标基质网格时，所述多个目标基质网格中目标基质网格的数量的取值范围为小于或等于10、且大于或等于2的整数。
6.根据权利要求1所述的一种网格划分方法，其特征在于，在从所述裂缝分布图像中确定至少两个裂缝网格之前，对所述裂缝分布图像进行二值化处理，得到二值化图像；
相应地，所述从所述裂缝分布图像中确定至少两个裂缝网格，包括：从所述二值化图像中确定至少两个裂缝网格；其中，所述二值化图像包括：裂缝和基质。
7.一种网格划分装置，其特征在于，所述装置提供目的油气藏的裂缝分布图像，其中，所述裂缝分布图像用于表征所述目的油气藏在水平二维面上的裂缝和基质的分布位置，所述装置包括：裂缝网格确定模块、基质点对应的裂缝网格确定模块、初始基质网格确定模块和目标基质网格确定模块；其中，
所述裂缝网格确定模块，用于从所述裂缝分布图像中确定至少两个裂缝网格；其中，所述裂缝网格包括所述裂缝分布图像中相邻两个裂缝节点之间的裂缝部分，所述裂缝节点包括：所述裂缝分布图像中的裂缝交点与裂缝端点；
所述基质点对应的裂缝网格确定模块，用于确定所述裂缝分布图像中基质点到所述裂缝网格的距离，并将与所述基质点距离最近的裂缝网格作为所述基质点对应的裂缝网格；
所述初始基质网格确定模块，用于将对应同一个裂缝网格的基质点划分至同一个初始基质网格中；
所述目标基质网格确定模块，用于对所述初始基质网格进行网格划分处理，得到多个目标基质网格，其中，包括：将所述初始基质网格转换为矩形初始基质网格；其中，所述矩形初始基质网格与所述初始基质网格一一对应，所述矩形初始基质网格的面积与所述初始基质网格的面积相同，且所述矩形初始基质网格的指定边长的长度与所述初始基质网格对应的裂缝网格的裂缝长度相同；对所述矩形初始基质网格进行对数加密处理，以使得将所述矩形初始基质网格划分为多个目标基质网格；其中，所述目标基质网格的一边与所述指定边长平行，且长度相同；所述多个目标基质网格中相邻两个目标基质网格中第一基质网格的面积为第二基质网格的面积的指定倍数；所述第一基质网格表示所述相邻两个目标基质网格中距离与所述初始基质网格对应的裂缝网格较远的目标基质网格；所述第二基质网格表示所述相邻两个目标基质网格中距离与所述初始基质网格对应的裂缝网格较近的目标基质网格。
8.根据权利要求7所述的一种网格划分装置，其特征在于，所述基质点到所述裂缝网格的距离表示所述基质点与所述裂缝网格内指定像素点之间的距离，所述指定像素点表示所述裂缝网格内的与所述基质点的距离最短的像素点。
CN201710934317.6A 2017-10-10 2017-10-10 一种网格划分方法及装置 Active CN107704685B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710934317.6A CN107704685B (zh) 2017-10-10 2017-10-10 一种网格划分方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710934317.6A CN107704685B (zh) 2017-10-10 2017-10-10 一种网格划分方法及装置

Publications (2)

Publication Number Publication Date
CN107704685A CN107704685A (zh) 2018-02-16
CN107704685B true CN107704685B (zh) 2020-04-17

Family

ID=61184759

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710934317.6A Active CN107704685B (zh) 2017-10-10 2017-10-10 一种网格划分方法及装置

Country Status (1)

Country Link
CN (1) CN107704685B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110473255B (zh) * 2019-07-17 2020-12-01 南京畅淼科技有限责任公司 一种基于多重网格划分的船舶系船柱定位方法
CN116596921B (zh) * 2023-07-14 2023-10-20 济宁市质量计量检验检测研究院(济宁半导体及显示产品质量监督检验中心、济宁市纤维质量监测中心) 一种焚烧炉渣分选方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101598804A (zh) * 2008-06-04 2009-12-09 中国石油天然气集团公司 确定地下油气藏构造的三维方法
CN102110183A (zh) * 2010-12-30 2011-06-29 中国石油化工股份有限公司胜利油田分公司地质科学研究院 一种反映流体沿储层大裂缝窜流的数值模拟方法
CN102243680A (zh) * 2011-07-21 2011-11-16 中国科学技术大学 一种网格划分方法及系统
CN104392109A (zh) * 2014-11-04 2015-03-04 中国石油天然气股份有限公司 一种油藏的动态离散裂缝模拟方法和装置
WO2016011064A2 (en) * 2014-07-15 2016-01-21 Petroleum Fractured Reservoir Solutions, Llc Discrete irregular cellular models for simulating the development of fractured reservoirs
CN105913494A (zh) * 2016-03-30 2016-08-31 北京大学 多尺度裂缝精细地质建模及数值模拟方法和装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101598804A (zh) * 2008-06-04 2009-12-09 中国石油天然气集团公司 确定地下油气藏构造的三维方法
CN102110183A (zh) * 2010-12-30 2011-06-29 中国石油化工股份有限公司胜利油田分公司地质科学研究院 一种反映流体沿储层大裂缝窜流的数值模拟方法
CN102243680A (zh) * 2011-07-21 2011-11-16 中国科学技术大学 一种网格划分方法及系统
WO2016011064A2 (en) * 2014-07-15 2016-01-21 Petroleum Fractured Reservoir Solutions, Llc Discrete irregular cellular models for simulating the development of fractured reservoirs
CN104392109A (zh) * 2014-11-04 2015-03-04 中国石油天然气股份有限公司 一种油藏的动态离散裂缝模拟方法和装置
CN105913494A (zh) * 2016-03-30 2016-08-31 北京大学 多尺度裂缝精细地质建模及数值模拟方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
二维PEBI网格的生成;蔡强;《工程图学学报》;20051231(第2期);第69-72页 *

Also Published As

Publication number Publication date
CN107704685A (zh) 2018-02-16

Similar Documents

Publication Publication Date Title
JP6998964B2 (ja) ジオフェンスのインデックスグリッドを判断するための方法及び装置
Zhang et al. Optimizing high-resolution Community Earth System Model on a heterogeneous many-core supercomputing platform
CN108106964B (zh) 一种确定页岩含气量的方法及装置
CN109992640B (zh) 位置网格的确定方法及装置、设备及存储介质
CN107704685B (zh) 一种网格划分方法及装置
JP2022508028A (ja) ３次元画像処理におけるデータの読み書き方法とシステム、記憶媒体及び端末
CN108133005B (zh) 一种基于内存数据库的环境模型模拟方法、终端设备及存储介质
CN107292002B (zh) 一种数字岩心重构的方法及装置
Geyer et al. Limits of reproducibility and hydrodynamic noise in atmospheric regional modelling
US10318676B2 (en) Techniques for statistical frequency enhancement of statically timed designs
US9576085B2 (en) Selective importance sampling
Wang et al. The GPU version of LASG/IAP Climate System Ocean Model version 3 (LICOM3) under the heterogeneous-compute interface for portability (HIP) framework and its large-scale application
CN103713953A (zh) 一种内存数据的搬移装置及方法
CN105653785A (zh) 一种雷达杂波模拟方法、装置及系统
Oberhuber et al. TNL: Numerical library for modern parallel architectures
CN115374237B (zh) 一种基于北斗网格码的矢量空间数据存储与查询方法
CN107844637B (zh) 一种油藏数值确定方法及装置
US8713485B2 (en) Categorization of design rule errors
CN107563617B (zh) 一种确定油气运移路径上的油气量的方法及装置
Anuradha et al. Efficient workload characterization technique for heterogeneous processors
CN107256281B (zh) 一种基于裁剪法的fpga可重构资源非矩形布局方法
US10339107B2 (en) Multi-level colocation and processing of spatial data on MapReduce
Feinberg et al. Modeling analog tile-based accelerators using sst
Téllez‐Velázquez et al. A Spark image processing toolkit
Stamoulias et al. Hardware accelerators for financial applications in HDL and High Level Synthesis

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