CN116108617A - 一种页岩油藏局部网格加密方法、装置及存储介质 - Google Patents

一种页岩油藏局部网格加密方法、装置及存储介质 Download PDF

Info

Publication number
CN116108617A
CN116108617A CN202211444807.5A CN202211444807A CN116108617A CN 116108617 A CN116108617 A CN 116108617A CN 202211444807 A CN202211444807 A CN 202211444807A CN 116108617 A CN116108617 A CN 116108617A
Authority
CN
China
Prior art keywords
crack
grid
conductivity
grids
threshold value
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.)
Pending
Application number
CN202211444807.5A
Other languages
English (en)
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.)
Xian Shiyou University
Original Assignee
Xian Shiyou University
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 Xian Shiyou University filed Critical Xian Shiyou University
Priority to CN202211444807.5A priority Critical patent/CN116108617A/zh
Publication of CN116108617A publication Critical patent/CN116108617A/zh
Pending legal-status Critical Current

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
    • 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
    • G06T17/205—Re-meshing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • 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)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本申请公开了一种页岩油藏局部网格加密方法、装置及存储介质，涉及石油开采技术领域，解决了现有技术中不能针对局部网格进行加密的问题。该方法包括，将裂缝嵌入原始背景网格中，计算裂缝导流能力；根据裂缝导流能力划分裂缝控制区域，计算裂缝控制区域内的网格与裂缝之间的传导率，进一步根据传导率计算裂缝控制区域内的网格与裂缝之间的等效距离，设定门限阈值，判断等效距离与门限阈值的大小；对传导率大于等于门限阈值的网格，多次进行二倍网格加密，直至网格的传导率均小于门限阈值；并对网格的传导率重新进行计算。该方法实现了对局部网格进行加密，降低了非稳态渗流对流动计算精度的影响，提高了页岩油藏嵌入离散裂缝模型的数值模拟精度。

Description

一种页岩油藏局部网格加密方法、装置及存储介质
技术领域
本申请涉及石油开采技术领域，尤其涉及一种页岩油藏局部网格加密方法、装置及存储介质。
背景技术
页岩基质渗透率极低，仅为几百纳达西到微达西，是常规油藏储层的百万分之一，微裂缝、水平页理缝以及建设性成岩作用形成的次生孔隙是游离液态烃赋存的主要空间，常规开发方式难以有效动用。
目前，对于嵌入式离散裂缝模型体积压裂改造区内，需要更高精度的网格，以降低基质-裂缝非稳态渗透流对流动计算精度的影响。但是如果将全区采用高精度网格体系，则会大大提高模型计算节点数量，大幅度降低计算效率。如何针对局部网格进行加密是目前亟需解决的问题。
发明内容
本申请实施例通过提供一种油藏局部网格加密方法、装置及存储介质，解决了现有技术中不能针对局部网格进行加密的问题，实现了对局部网格进行加密，降低了非稳态渗流对流动计算精度的影响，提高了页岩油藏嵌入离散裂缝模型的数值模拟精度。
第一方面，本发明实施例提供了一种页岩油藏局部网格加密方法，该方法包括：
将裂缝嵌入原始背景网格中，计算裂缝导流能力；
根据所述裂缝导流能力划分裂缝控制区域，计算所述裂缝控制区域内的网格与所述裂缝之间的传导率，进一步根据所述传导率计算所述裂缝控制区域内的网格与所述裂缝之间的等效距离；
设定门限阈值，判断所述等效距离与所述门限阈值的大小；
对所述传导率大于等于所述门限阈值的所述网格，多次进行二倍网格加密，直至所述网格的传导率均小于所述门限阈值；
并对所述网格的传导率重新进行计算。
结合第一方面，在一种可能的实现方式中，所述将裂缝嵌入原始背景网格中，包括：根据所述裂缝与所述网格的相对位置，选择不同的嵌入方式将所述裂缝嵌入所述原始背景网格中。
结合第一方面，在一种可能的实现方式中，所述根据所述裂缝导流能力划分裂缝控制区域，包括：根据以下公式进行计算划定裂缝控制区域：
RD＝{P|D(f,P)＜c·tf}，
其中，c表示松弛系数；D(f,P)表示点P到所述裂缝的距离，其计算公式为：
Figure BDA0003949777690000021
tf表示裂缝导流能力，其计算公式为：tf＝kfdf，其中kf表示嵌入单位裂缝的渗透率；df表示裂缝宽度。
结合第一方面，在一种可能的实现方式中，所述计算所述裂缝控制区域内的网格与所述裂缝之间的传导率，包括：采用网格等效距离方法计算，计算公式为：
Figure BDA0003949777690000022
其中，Afm表示嵌入裂缝的面积；kf表示嵌入单位裂缝的渗透率；km表示背景单元的渗透率；＜d＞表示背景网格与嵌入式离散裂缝之间的等效距离；n表示嵌入裂缝的法向量；x表示背景网格内的一点到裂缝的距离向量，S表示背景网格的体积。
结合第一方面，在一种可能的实现方式中，所述进行二倍网格加密，包括：沿XYZ方向将所述网格分为两个网格，并对所述两个网格加表示密为8个网格或沿XY方向将所述网格分为两个网格，并对所述两个网格加密为4个网格。
结合第一方面，在一种可能的实现方式中，所述并对所述网格的传导率重新进行计算，包括：计算任意两个相邻网格之间的传导率，其计算公式为：
Figure BDA0003949777690000031
其中，Al表示所述任意两个相邻网格接触部分的面积；k表示网格渗透率；D1和D2分别表示所述任意两个相邻网格的中心点到公共面的距离；ccoor表示局部网格传导率校正系数。
第二方面，本发明实施例提供了一种页岩油藏局部网格加密装置，该装置包括：
导流计算模块，用于将裂缝嵌入原始背景网格中，计算裂缝导流能力；
等效距离计算模块，用于根据所述裂缝导流能力划分裂缝控制区域，计算所述裂缝控制区域内的网格与所述裂缝之间的传导率，进一步根据所述传导率计算所述裂缝控制区域内的网格与所述裂缝之间的等效距离；
判断模块，用于设定门限阈值，判断所述等效距离与所述门限阈值的大小；
网格加密模块，用于对所述传导率大于等于所述门限阈值的所述网格进行多次二倍网格加密，直至所述网格的传导率均小于所述门限阈值；
二次传导计算模块，用于并对所述网格的传导率重新进行计算。
结合第二方面，在一种可能的实现方式中，所述导流计算模块用于根据所述裂缝与所述网格的相对位置，选择不同的嵌入方式将所述裂缝嵌入所述原始背景网格中。
结合第二方面，在一种可能的实现方式中，所述等效距离计算模块用于根据以下公式进行计算划定裂缝控制区域：
RD＝{P|D(f,P)＜c·tf}，
其中，c表示松弛系数；D(f,P)表示点P到所述裂缝的距离，其计算公式为：
Figure BDA0003949777690000032
tf表示裂缝导流能力，其计算公式为：tf＝kfdf，其中kf表示嵌入单位裂缝的渗透率；df表示裂缝宽度。
结合第二方面，在一种可能的实现方式中，所述等效距离计算模块用于采用网格等效距离方法计算，计算公式为：
Figure BDA0003949777690000041
其中，Afm表示嵌入裂缝的面积；kf表示嵌入单位裂缝的渗透率；km表示背景单元的渗透率；＜d＞表示背景网格与嵌入式离散裂缝之间的等效距离；n表示嵌入裂缝的法向量；x表示背景网格内的一点到裂缝的距离向量，S表示背景网格的体积。
结合第二方面，在一种可能的实现方式中，所述网格加密模块用于沿XYZ方向将所述网格分为两个网格，并对所述两个网格加表示密为8个网格或沿XY方向将所述网格分为两个网格，并对所述两个网格加密为4个网格。
结合第二方面，在一种可能的实现方式中，所述二次传导计算模块用于计算任意两个相邻网格之间的传导率，其计算公式为：
Figure BDA0003949777690000042
其中，Al表示所述任意两个相邻网格接触部分的面积；k表示网格渗透率；D1和D2分别表示所述任意两个相邻网格的中心点到公共面的距离；ccoor表示局部网格传导率校正系数。
第三方面，本发明实施例提供了一种页岩油藏局部网格加密服务器，包括存储器和处理器；
所述存储器用于存储计算机可执行指令；
所述处理器用于执行所述计算机可执行指令，以实现第一方面及第一方面任一项所述的方法。
第四方面，本发明实施例提供了一种计算机可读存储介质，所述计算机可读存储介质存储有可执行指令，计算机执行所述可执行指令时能够实现第一方面及第一方面任一项所述的方法。
本发明实施例中提供的一个或多个技术方案，至少具有如下技术效果或优点：
本发明实施例通过采用了一种页岩油藏局部网格加密方法、装置及存储介质，该方法包括，将裂缝嵌入原始背景网格中，计算裂缝导流能力；根据裂缝导流能力划分裂缝控制区域，计算裂缝控制区域内的网格与裂缝之间的传导率，进一步根据传导率计算裂缝控制区域内的网格与裂缝之间的等效距离，设定门限阈值，判断等效距离与门限阈值的大小；对传导率大于等于门限阈值的网格，多次进行二倍网格加密，直至网格的传导率均小于门限阈值；并对网格的传导率重新进行计算；将裂缝嵌入原始背景网格中，使用模型对实际情况进行估算，门限阈值的大小是可以自由设置的，可以根据对计算结果的要求进行相应阈值大小的设置，筛选出需要进行加密的网格，进行加密，有效解决了现有技术中不能针对局部网格进行加密的问题，实现了对局部网格进行加密，降低了非稳态渗流对流动计算精度的影响，提高了页岩油藏嵌入离散裂缝模型的数值模拟精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对本发明实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的网格加密方法的流程图；
图2为本申请实施例提供的不同裂缝嵌入原始背景网格的示意图；
图3为本申请实施例提供的页岩油藏压裂水平模型示意图；
图4为本申请实施例提供的裂缝嵌入原始背景网格示意图；
图5为本申请实施例提供的划分裂缝控制区域示意图；
图6为本申请实施例提供的局部网格加密区域示意图；
图7为本申请实施例提供的局部网格第一次加密结果示意图；
图8为本申请实施例提供的局部网格第二次加密结果示意图；
图9本申请实施例提供的页岩油藏局部网格加密装置示意图；
图10本申请实施例提供的页岩油藏局部网格加密服务器示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
目前，我国页岩石油具有以下几个特点，(1)页岩基质渗透率极低，仅为几百纳达西到微达西，是常规油藏储层的百万分之一，微裂缝、水平页理缝以及建设性成岩作用形成的次生孔隙是游离液态烃赋存主要空间，常规开发方式难以有效动用；(2)由于热演化程度较高，处于液态烃与气态烃并存窗口，以密度较小原油为主，油质较轻、气油比较高、可流动性好，依靠水平井体积压裂技术可实现经济效益开发；(3)陆相混积岩、含碳酸盐质页岩地层脆性矿物含量高，利用人工压裂改造技术形成有效人造缝网系统，可实现较高单井产量和累计EUR(单井累计采出量)。
研究与实践表明，水平井+体积压裂是实现陆相页岩油藏规模化开发的有效方式：水力压裂形成的裂缝网络为流体流动提供了高渗通道，在提高地层渗透率的同时扩大了储层与井筒的接触面积，从而达到增产的效果。
为了进一步实现精细化与经济化开发，需要高效且相等的压裂水平井产能预测方法，为了精细刻画压裂裂缝网对流动的影响，人们开始基于离散裂缝模型对裂缝进行表征。离散裂缝模型(全称：Discrete_Fracture_Network，简称：DFM)是将裂缝按照其真实分布及形态，进行完整和显性的描述的模型，其发展来源于20世纪70年代的离散裂缝建模(全称：Discrete Fracture Network Model，简称：DFN)方法的提出，并随着非常规油气藏中人工压裂技术的广泛使用。早期的离散裂缝模型都是基于匹配型网格，即将裂缝作为内边界并以此为约束面来进行网格剖分。
由于裂缝几何形态的复杂性，需采用非结构化网格技术，其剖分过程非常复杂和繁琐，尤其当裂缝间距离或夹角很小时，常常因网格划分质量差导致计算困难。对此，Lee等人和Moinfar等人提出了嵌入式离散裂缝模型(全称：Embedded Discrete FractureModel，简称：EDFM)，该模型将裂缝网络直接嵌入基岩结构化网格系统中，避免了上述复杂的非结构化网格剖分过程。虽然仍需要计算裂缝与网格之间的几何信息，但相对于复杂的非结构化网格剖分过程，其计算复杂度大大降低，从而能够提高计算效率。对于嵌入式离散裂缝模型体积压裂改造区内，需要更高精度的(分辨率更高的)网格，以降低基质-裂缝非稳态渗流对流动计算精度的影响。但是如果将全区采用高精度网格体系，则将会大大提高模型计算节点数量，大幅降低计算效率。因此，最优方法是采用局部网格加密方法，仅在体积压裂改造区内进行网格加密。然而，目前缺少针对性的局部网格加密技术，导致仍然采用均一网格模型进行数值模拟计算。
本发明实施例提供了一种页岩油藏局部网格加密方法，如图1所示该方法包括以下步骤S101至S105。
S101，将裂缝嵌入原始背景网格中，计算裂缝导流能力。
S102，根据裂缝导流能力划分裂缝控制区域，计算裂缝控制区域内的网格与裂缝之间的传导率，进一步根据传导率计算裂缝控制区域内的网格与裂缝之间的等效距离。
S103，设定门限阈值，判断等效距离与门限阈值的大小。
S104，对传导率大于等于门限阈值的网格，多次进行二倍网格加密，直至网格的传导率均小于门限阈值。
S105，并对网格的传导率重新进行计算。
本申请提供的方法能够将针对需要加密的网格进行局部加密，进而减轻了计算的压力，并且能够提高计算效率。本申请提供的方法能够仅仅在体积压裂改造区内进行网格加密，这种方式能够降低基质-裂缝非稳态渗流对流动计算精度的影响，提高了页岩油藏嵌入式裂缝模型的数据模拟计算精度。
在步骤S101中，将裂缝嵌入原始背景网格中，包括：根据裂缝与网格的相对位置，选择不同的嵌入方式将裂缝嵌入原始背景网格中。将体积压裂封网进行嵌入离散裂缝进行建模，即：计算裂缝对变形与背景网格的几何相交关系，截取落在目标背景网格内的多边形部分，将裂缝多边形切分成很多小的多边形。根据裂缝与网格的相对位置，有不同的嵌入情形，如图2所示为不同裂缝嵌入原始背景网格的情形。
在步骤S102中，根据裂缝导流能力划分裂缝控制区域，包括：根据以下公式进行计算划定裂缝控制区域：
RD＝{P|D(f,P)＜c·tf}，
其中，c表示松弛系数；D(f,P)表示点P到裂缝的距离，其计算公式为：
Figure BDA0003949777690000081
tf表示裂缝导流能力，其计算公式为：tf＝kfdf，其中kf表示嵌入单位裂缝的渗透率；df表示裂缝宽度。
在步骤S102中，计算裂缝控制区域内的网格与裂缝之间的传导率，包括：采用网格等效距离方法计算，计算公式为：
Figure BDA0003949777690000082
其中，Afm表示嵌入裂缝的面积；kf表示嵌入单位裂缝的渗透率；km表示背景单元的渗透率；＜d＞表示背景网格与嵌入式离散裂缝之间的等效距离；n表示嵌入裂缝的法向量；x表示背景网格内的一点到裂缝的距离向量，S表示背景网格的体积。
步骤S102和步骤S103，根据裂缝导流能力划分裂缝控制区域，对于划分好的区域进行计算裂缝导流能力，并将其和设定的门限阈值进行比较，在这一过程中，由于门限阈值是设定的，即可以根据自身对结果的要求进行门限阈值的设定，需要计算结果更加精确，则门限阈值就设定的较小，若需要计算速度较快，则门限阈值就设定的较大。对于精度和计算速度，在使用的过程中根据实际情况不同进行取舍。
在步骤S104中，进行二倍网格加密，包括：沿XYZ方向将网格分为两个网格，并对两个网格加表示密为8个网格或沿XY方向将网格分为两个网格，并对两个网格加密为4个网格。在步骤S104中需要注意的是，为了提高模型的计算精度，网格不一定是裂缝实际相交的网格，即对于裂缝控制区域内的网格都需要进行判断。
对进行过一次加密后的网格继续进行判断，判断条件为步骤S103中条件，入不符合则在此进行二倍网格加密，当然此处提及的二倍网格知识一种加密网格，也可以为多倍网格加密。
由于模型采用了多重网格局部加密，常规的油藏数值模拟骤然间无法直接解析，需要重新计算相邻网格之间的传导率，作为油藏数值。在步骤S105中，并对网格的传导率重新进行计算，包括：计算任意两个相邻网格之间的传导率，其计算公式为：
Figure BDA0003949777690000091
其中，Al表示任意两个相邻网格接触部分的面积；k表示网格渗透率；D1和D2分别表示任意两个相邻网格的中心点到公共面的距离；ccoor表示局部网格传导率校正系数。ccoor的计算公式为：
Figure BDA0003949777690000092
其中DH12表示任意两个相邻网格中心连接线水平方向分量的平方，DV12表示任意两个相邻网格中心连接线垂直方向分量的平方。
Figure BDA0003949777690000101
其中，P1和P2分别为任意两个相邻网格中心坐标，
Figure BDA0003949777690000102
和
Figure BDA0003949777690000103
分别为x、y和z方向的单位向量。
本申请提供一个具体的实施例来对本方法进行说明。如图3所示为页岩油藏压裂水平井模型示意图。详细参数如表1所示。
表1机理模型参数
Figure BDA0003949777690000104
结合步骤S101和S102将裂缝嵌入原始背景网格，并划分裂缝控制区域，如图4中不规则网格区域。结合步骤S103计算等效距离并划分局部网格加密区，如图5所示。根据S104，对局部网格加密区进行网格加密，如图6所示，为局部网格加密区域示意图，图7和图8分别为一次加密和二次加密结果图。根据步骤S105计算传导率，将计算结果与原始模型进行对比，并与5m网格尺寸参考模型进行对比，如表2所示。
表2对比模型参数
Figure BDA0003949777690000105
本发明实施例提供了一种页岩油藏局部网格加密装置，如图9所示该装置包括：导流计算模块901、等效距离计算模块902、判断模块903、网格加密模块904以及二次传导计算模块905。
导流计算模块901，用于将裂缝嵌入原始背景网格中，计算裂缝导流能力；导流计算模块901用于根据裂缝与网格的相对位置，选择不同的嵌入方式将裂缝嵌入原始背景网格中。
等效距离计算模块902，用于根据裂缝导流能力划分裂缝控制区域，计算裂缝控制区域内的网格与裂缝之间的传导率，进一步根据传导率计算裂缝控制区域内的网格与裂缝之间的等效距离。等效距离计算模块902用于根据以下公式进行计算划定裂缝控制区域：RD＝{P|D(f,P)＜c·tf}，其中，c表示松弛系数；D(f,P)表示点P到裂缝的距离，其计算公式为：
Figure BDA0003949777690000111
tf表示裂缝导流能力，其计算公式为：tf＝kfdf，其中kf表示嵌入单位裂缝的渗透率；df表示裂缝宽度。等效距离计算模块902用于采用网格等效距离方法计算，计算公式为：
Figure BDA0003949777690000112
其中，Afm表示嵌入裂缝的面积；kf表示嵌入单位裂缝的渗透率；km表示背景单元的渗透率；＜d＞表示背景网格与嵌入式离散裂缝之间的等效距离；n表示嵌入裂缝的法向量；x表示背景网格内的一点到裂缝的距离向量，S表示背景网格的体积。
判断模块903，用于设定门限阈值，判断等效距离与门限阈值的大小。
网格加密模块904，用于对传导率大于等于门限阈值的网格进行多次二倍网格加密，直至网格的传导率均小于门限阈值；网格加密模块904用于沿XYZ方向将网格分为两个网格，并对两个网格加表示密为8个网格或沿XY方向将网格分为两个网格，并对两个网格加密为4个网格。
二次传导计算模块905，用于并对网格的传导率重新进行计算。二次传导计算模块905用于计算任意两个相邻网格之间的传导率，其计算公式为：
Figure BDA0003949777690000113
其中，Al表示任意两个相邻网格接触部分的面积；k表示网格渗透率；D1和D2分别表示任意两个相邻网格的中心点到公共面的距离；ccoor表示局部网格传导率校正系数。
上述实施例阐明的装置或模块，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。当然，也可以将实现某功能的模块由多个子模块或子单元组合实现。
本发明实施例提供了一种页岩油藏局部网格加密服务器，如图10所示包括存储器1001和处理器1002；存储器1001用于存储计算机可执行指令；处理器1002用于执行计算机可执行指令，以实现本申请提供的方法。
本发明实施例提供了一种计算机可读存储介质，计算机可读存储介质存储有可执行指令，计算机执行可执行指令时能够实现本申请提供的方法。
上述存储介质包括但不限于随机存取存储器(英文：Random Access Memory；简称：RAM)、只读存储器(英文：Read-Only Memory；简称：ROM)、缓存(英文：Cache)、硬盘(英文：Hard Disk Drive；简称：HDD)或者存储卡(英文：Memory Card)。所述存储器可以用于存储计算机程序指令。
本申请中所述的方法、装置或模块可以以计算机可读程序代码方式实现控制器按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(英文：Application Specific Integrated Circuit；简称：ASIC)、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C9051F320，存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的硬件的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，也可以通过数据迁移的实施过程中体现出来。该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，移动终端，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
以上实施例仅用以说明本申请的技术方案，而非对本申请限制；尽管参照前述实施例对本申请进行了详细的说明，本领域普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请技术方案的范围。

Claims (9)

1.一种页岩油藏局部网格加密方法，其特征在于，包括：
将裂缝嵌入原始背景网格中，计算裂缝导流能力；
根据所述裂缝导流能力划分裂缝控制区域，计算所述裂缝控制区域内的网格与所述裂缝之间的传导率，进一步根据所述传导率计算所述裂缝控制区域内的网格与所述裂缝之间的等效距离；
设定门限阈值，判断所述等效距离与所述门限阈值的大小；
对所述传导率大于等于所述门限阈值的所述网格，多次进行二倍网格加密，直至所述网格的传导率均小于所述门限阈值；
并对所述网格的传导率重新进行计算。
2.根据权利要求1所述的方法，其特征在于，所述将裂缝嵌入原始背景网格中，包括：根据所述裂缝与所述网格的相对位置，选择不同的嵌入方式将所述裂缝嵌入所述原始背景网格中。
3.根据权利要求1所述的方法，其特征在于，所述根据所述裂缝导流能力划分裂缝控制区域，包括：根据以下公式进行计算划定裂缝控制区域：
R D＝{P|D(f,P)＜c·tf}，
其中，c表示松弛系数；D(f,P)表示点P到所述裂缝的距离，其计算公式为：
Figure FDA0003949777680000011
tf表示裂缝导流能力，其计算公式为：tf＝kfdf，其中kf表示嵌入单位裂缝的渗透率；df表示裂缝宽度。
4.根据权利要求1所述的方法，其特征在于，所述计算所述裂缝控制区域内的网格与所述裂缝之间的传导率，包括：采用网格等效距离方法计算，计算公式为：
Figure FDA0003949777680000012
其中，Afm表示嵌入裂缝的面积；kf表示嵌入单位裂缝的渗透率；km表示背景单元的渗透率；D*表示背景网格与嵌入式离散裂缝之间的等效距离；n表示嵌入裂缝的法向量；x表示背景网格内的一点到裂缝的距离向量，S表示背景网格的体积。
5.根据权利要求1所述的方法，其特征在于，所述进行二倍网格加密，包括：沿XYZ方向将所述网格分为两个网格，并对所述两个网格加表示密为8个网格或沿XY方向将所述网格分为两个网格，并对所述两个网格加密为4个网格。
6.根据权利要求1所述的方法，其特征子在于，所述并对所述网格的传导率重新进行计算，包括：计算任意两个相邻网格之间的传导率，其计算公式为：
Figure FDA0003949777680000021
其中，Al表示所述任意两个相邻网格接触部分的面积；k表示网格渗透率；D1和D2分别表示所述任意两个相邻网格的中心点到公共面的距离；ccoor表示局部网格传导率校正系数。
7.一种页岩油藏局部网格加密装置，其特征在于，包括：
导流计算模块，用于将裂缝嵌入原始背景网格中，计算裂缝导流能力；
等效距离计算模块，用于根据所述裂缝导流能力划分裂缝控制区域，计算所述裂缝控制区域内的网格与所述裂缝之间的传导率，进一步根据所述传导率计算所述裂缝控制区域内的网格与所述裂缝之间的等效距离；
判断模块，用于设定门限阈值，判断所述等效距离与所述门限阈值的大小；
网格加密模块，用于对所述传导率大于等于所述门限阈值的所述网格进行多次二倍网格加密，直至所述网格的传导率均小于所述门限阈值；
二次传导计算模块，用于并对所述网格的传导率重新进行计算。
8.一种页岩油藏局部网格加密服务器，其特征在于，包括存储器和处理器；
所述存储器用于存储计算机可执行指令；
所述处理器用于执行所述计算机可执行指令，以实现权利要求1-6任一项所述的方法。
9.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有可执行指令，计算机执行所述可执行指令时能够实现如权利要求1-6任一项所述的方法。
CN202211444807.5A 2022-11-18 2022-11-18 一种页岩油藏局部网格加密方法、装置及存储介质 Pending CN116108617A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211444807.5A CN116108617A (zh) 2022-11-18 2022-11-18 一种页岩油藏局部网格加密方法、装置及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211444807.5A CN116108617A (zh) 2022-11-18 2022-11-18 一种页岩油藏局部网格加密方法、装置及存储介质

Publications (1)

Publication Number Publication Date
CN116108617A true CN116108617A (zh) 2023-05-12

Family

ID=86262728

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211444807.5A Pending CN116108617A (zh) 2022-11-18 2022-11-18 一种页岩油藏局部网格加密方法、装置及存储介质

Country Status (1)

Country Link
CN (1) CN116108617A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117390998A (zh) * 2023-12-08 2024-01-12 西安石油大学 一种针对嵌入式离散裂缝模型局部网格加密方法及装置
CN117610313A (zh) * 2024-01-18 2024-02-27 西安石油大学 页岩油藏人工裂缝的网格加密方法及装置

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117390998A (zh) * 2023-12-08 2024-01-12 西安石油大学 一种针对嵌入式离散裂缝模型局部网格加密方法及装置
CN117610313A (zh) * 2024-01-18 2024-02-27 西安石油大学 页岩油藏人工裂缝的网格加密方法及装置
CN117610313B (zh) * 2024-01-18 2024-04-26 西安石油大学 页岩油藏人工裂缝的网格加密方法及装置

Similar Documents

Publication Publication Date Title
CN116108617A (zh) 一种页岩油藏局部网格加密方法、装置及存储介质
CN103400020B (zh) 一种测算多条相交离散裂缝流动状况的油藏数值模拟方法
CN103399970B (zh) 用离散裂缝线处理进行数模测算油藏流动状况的方法
CN109829217B (zh) 压裂性裂缝油藏产能模拟方法及装置
CN111222271A (zh) 基于基质-裂缝非稳态窜流油藏裂缝数值模拟方法及系统
CN105160134B (zh) 致密储层多重介质中油气流动的混合介质模拟方法及装置
CN108399290B (zh) 基于裂隙网络连通性的裂隙流的modflow模拟方法
CN102968806A (zh) 一种不规则边界区域内自适应网格的生成方法
CN103279991A (zh) 一种利用离散裂缝端点变形数模提高油藏开发效果的方法
CN105260543A (zh) 基于双孔模型的多重介质油气流动模拟方法及装置
CN114154430A (zh) 一种压裂油藏co2驱油流动模拟方法
CN107066672B (zh) 一种代替气驱组分模型的数值模拟方法
CN105653746A (zh) 基于嵌入式离散裂缝模型的压裂井建模及模拟方法
CN116167302B (zh) 一种天然气水合物增产模拟中人工复杂裂缝的描述方法
Wang A hybrid dual-continuum discrete fracture modeling approach for numerical simulation of production from unconventional plays
CN109558614A (zh) 页岩气藏多尺度裂缝内气体流动的模拟方法及系统
Sun et al. A Powerful Thermal-EDFM for Modeling Enhanced Geothermal System With 3D Complex Fracture Networks
CN105512458A (zh) 缝洞型碳酸盐岩油藏衰竭开采特征的综合分析方法及装置
CN116971769A (zh) 一种砾岩水力裂缝扩展方向的预测方法
Wang et al. A non-uniform grid approach for high-resolution flood inundation simulation based on GPUs
CN114818551B (zh) 一种页岩油藏压裂缝网自适应半无限导流表征方法及装置
Jiang et al. Comparative analysis and risk assessment of dam-break floods: Taking Pingshuijiang Reservoir as an example
CN113360984A (zh) 一种考虑壁面捕获效应的支撑剂输送数值模拟方法
Yang et al. Advanced dual-porosity and dual-permeability model for tight rock integrated primary depletion and enhanced oil recovery simulation
CN117390998A (zh) 一种针对嵌入式离散裂缝模型局部网格加密方法及装置

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