CN112329353A - 窜流量确定方法、装置、计算机设备及存储介质 - Google Patents

窜流量确定方法、装置、计算机设备及存储介质 Download PDF

Info

Publication number
CN112329353A
CN112329353A CN202011025535.6A CN202011025535A CN112329353A CN 112329353 A CN112329353 A CN 112329353A CN 202011025535 A CN202011025535 A CN 202011025535A CN 112329353 A CN112329353 A CN 112329353A
Authority
CN
China
Prior art keywords
model
fluid
phase particles
target
determining
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.)
Granted
Application number
CN202011025535.6A
Other languages
English (en)
Other versions
CN112329353B (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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202011025535.6A priority Critical patent/CN112329353B/zh
Publication of CN112329353A publication Critical patent/CN112329353A/zh
Application granted granted Critical
Publication of CN112329353B publication Critical patent/CN112329353B/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
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • 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
    • G06F30/25—Design optimisation, verification or simulation using particle-based methods
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本申请实施例公开了一种窜流量确定方法、装置、计算机设备及存储介质，属于油气勘探开发领域。该方法包括：构建第一模型；模拟第一模型中的流体渗流，直至第一模型中的流体处于稳定状态，确定第一模型的出口端裂缝的流量，作为第一流量；构建第二模型，第二模型为裂缝模型，且第一模型与第二模型中的裂缝相同；模拟第二模型中的流体渗流，直至第二模型中的流体处于稳定状态，确定第二模型的出口端裂缝的流量，作为第二流量；将第一流量与第二流量之间的差值，确定为裂缝‑孔隙双重介质的窜流量。该方法提高了效率，而且避免了人工操作中主观因素的影响，提供了准确的窜流量数据，有助于操作人员客观准确地认识油气井的产能情况。

Description

窜流量确定方法、装置、计算机设备及存储介质
技术领域
本申请实施例涉及油气勘探开发领域，特别涉及一种窜流量确定方法、装置、计算机设备及存储介质。
背景技术
在油气勘探开发领域，窜流量反映了基质流体向裂缝窜流能力的大小，是研究碳酸盐岩油气藏生产能力的重要参数之一。目前通常是由操作人员对试井曲线进行分析，得到储层平均窜流能力，但是对操作人员的要求较高，分析效率低，会耗费大量人力成本，并且不同操作人员得出的结果可能存在差异，准确性低。
发明内容
本申请实施例提供了一种窜流量确定方法、装置、计算机设备及存储介质，提供了一种准确确定裂缝-孔隙双重介质的窜流量的方法，有助于操作人员客观准确地认识油气井的产能情况，为后续工作提供了数据支持。所述技术方案如下：
一方面，提供了一种窜流量确定方法，所述方法包括：
构建第一模型，所述第一模型为裂缝-孔隙双重介质模型；
模拟所述第一模型中的流体渗流，直至所述第一模型中的流体处于稳定状态，确定所述第一模型的出口端裂缝的流量，作为第一流量；
构建第二模型，所述第二模型为裂缝模型，且所述第一模型与所述第二模型中的裂缝相同；
模拟所述第二模型中的流体渗流，直至所述第二模型中的流体处于稳定状态，确定所述第二模型的出口端裂缝的流量，作为第二流量；
将所述第一流量与所述第二流量之间的差值，确定为所述裂缝-孔隙双重介质的窜流量。
在一种可能实现方式中，所述模拟所述第一模型中的流体渗流，直至所述第一模型中的流体处于稳定状态，确定所述第一模型的出口端裂缝的流量，作为第一流量，包括：
对所述第一模型中的流体施加预设压力梯度，以使所述流体按照格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述第一模型的出口端裂缝的流量，作为所述第一流量。
在另一种可能实现方式中，所述模拟所述第二模型中的流体渗流，直至所述第二模型中的流体处于稳定状态，确定所述第二模型的出口端裂缝的流量，作为第二流量，包括：
对所述第二模型中的流体施加所述预设压力梯度，以使所述流体按照所述格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述第二模型的出口端裂缝的流量，作为所述第二流量。
在另一种可能实现方式中，所述第一模型包括固相颗粒、固相颗粒之间的孔隙和裂缝，且所述第一模型由网格构成；所述构建第二模型，包括：
为所述第一模型中的孔隙对应的网格添加第一标记，得到所述第二模型，所述第一标记表示对应的网格为固相颗粒。
在另一种可能实现方式中，所述构建所述第一模型，包括：
在目标三维区域中添加所述固相颗粒，以使添加的固相颗粒之间构成孔隙，得到第三模型，所述第三模型为孔隙模型；
对所述第三模型中的流体施加预设压力梯度，以使所述流体按照格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态；
在所述第三模型中添加裂缝，得到所述第一模型。
在另一种可能实现方式中，所述在目标三维区域中添加固相颗粒，以使固相颗粒之间构成孔隙，得到第三模型，包括：
获取目标参数，所述目标参数包括固相颗粒的半径及孔隙度，所述孔隙度指示所述第三模型中的孔隙空间体积之和与所述第三模型体积的比值；
根据所述半径，在所述目标三维区域中添加具有所述半径的球体固相颗粒，直至所述目标三维区域中的孔隙空间体积之和与所述目标三维区域体积的比值等于所述孔隙度时，将所述目标三维区域作为所述第三模型。
在另一种可能实现方式中，所述固相颗粒的半径为多个，所述目标参数还包括每个半径对应的体积占比，所述体积占比为同一半径对应的固相颗粒的体积与固相颗粒总体积的比值，所述固相颗粒总体积根据所述孔隙度确定；
所述根据所述半径，在所述目标三维区域中添加具有所述半径的球体固相颗粒，直至所述目标三维区域中的孔隙空间体积之和与所述目标三维区域体积的比值等于所述孔隙度时，将所述目标三维区域作为所述第三模型，包括：
按照多个半径的排列顺序，遍历所述多个半径；
根据第一个半径，在所述目标三维区域中添加具有所述第一个半径的球体固相颗粒，直至所述目标三维区域中所述第一个半径对应的固相颗粒的体积与所述固相颗粒总体积的比值，等于所述第一个半径对应的体积占比；
根据下一个半径，在所述目标三维区域中添加具有所述下一个半径的球体固相颗粒，直至所述目标三维区域中所述下一个半径对应的固相颗粒的体积与所述固相颗粒总体积的比值，等于所述下一个半径对应的体积占比。
在另一种可能实现方式中，所述根据所述半径，在所述目标三维区域中添加具有所述半径的球体固相颗粒之前，所述方法还包括：
为所述目标三维区域中的每个网格添加第二标记，所述第二标记表示对应的网格为孔隙；
所述根据所述半径，在所述目标三维区域中添加具有所述半径的球体固相颗粒，包括：
在所述目标三维区域中随机指定一个坐标点；
为以所述坐标点为球心，具有所述半径的球体区域内的网格添加所述第一标记。
在另一种可能实现方式中，所述在所述第三模型中添加裂缝，得到所述第一模型，包括：
获取裂缝开度；
根据所述裂缝开度，在所述第三模型中随机确定裂缝对应的网格；
为所述裂缝对应的网格添加所述第三标记，所述第三标记表示对应的网格为裂缝。
在另一种可能实现方式中，所述对所述第一模型中的流体施加预设压力梯度，以使所述流体按照所述格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述第一模型的出口端裂缝的流量，作为所述第一流量，包括：
对所述第一模型中的流体施加预设压力梯度，以使所述流体按照所述格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述格子-玻尔兹曼演化关系中的流体粒子分布函数的函数值；
根据所述函数值，确定所述第一模型中的流体速度场；
从所述流体速度场中获取裂缝出口端的流体速度，将所述裂缝出口端的流体速度之和，作为所述第一流量。
在另一种可能实现方式中，所述对所述第二模型中的流体施加预设压力梯度，以使所述流体按照所述格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述第二模型的出口端裂缝的流量，作为所述第二流量，包括：
对所述第二模型中的流体施加所述预设压力梯度，以使所述流体按照所述格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述格子-玻尔兹曼演化关系中的流体粒子分布函数的函数值；
根据所述函数值，确定所述第二模型中的流体速度场；
从所述流体速度场中获取裂缝出口端的流体速度，将所述裂缝出口端的流体速度之和，作为所述第二流量。
另一方面，提供了一种窜流量确定装置，所述装置包括：
第一构建模块，用于构建第一模型，所述第一模型为裂缝-孔隙双重介质模型；
第一模拟模块，用于模拟所述第一模型中的流体渗流，直至所述第一模型中的流体处于稳定状态，确定所述第一模型的出口端裂缝的流量，作为第一流量；
第二构建模块，用于构建第二模型，所述第二模型为裂缝模型，且所述第一模型与所述第二模型中的裂缝相同；
第二模拟模块，用于模拟所述第二模型中的流体渗流，直至所述第二模型中的流体处于稳定状态，确定所述第二模型的出口端裂缝的流量，作为第二流量；
确定模块，用于将所述第一流量与所述第二流量之间的差值，确定为所述裂缝-孔隙双重介质的窜流量。
在一种可能实现方式中，所述第一模拟模块，用于：
对所述第一模型中的流体施加预设压力梯度，以使所述流体按照格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述第一模型的出口端裂缝的流量，作为所述第一流量。
在另一种可能实现方式中，所述第二模拟模块，用于：
对所述第二模型中的流体施加所述预设压力梯度，以使所述流体按照所述格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述第二模型的出口端裂缝的流量，作为所述第二流量。
在另一种可能实现方式中，所述第一模型包括固相颗粒、固相颗粒之间的孔隙和裂缝，且所述第一模型由网格构成；所述第二构建模块，包括：
第一标记添加单元，用于为所述第一模型中的孔隙对应的网格添加第一标记，得到所述第二模型，所述第一标记表示对应的网格为固相颗粒。
在另一种可能实现方式中，所述第一构建模块，包括：
颗粒添加单元，用于在目标三维区域中添加所述固相颗粒，以使添加的固相颗粒之间构成孔隙，得到第三模型，所述第三模型为孔隙模型；
流体演化单元，用于对所述第三模型中的流体施加预设压力梯度，以使所述流体按照格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态；
裂缝添加单元，用于在所述第三模型中添加裂缝，得到所述第一模型。
在另一种可能实现方式中，所述颗粒添加单元，用于：
获取目标参数，所述目标参数包括固相颗粒的半径及孔隙度，所述孔隙度指示所述第三模型中的孔隙空间体积之和与所述第三模型体积的比值；
根据所述半径，在所述目标三维区域中添加具有所述半径的球体固相颗粒，直至所述目标三维区域中的孔隙空间体积之和与所述目标三维区域体积的比值等于所述孔隙度时，将所述目标三维区域作为所述第三模型。
在另一种可能实现方式中，所述固相颗粒的半径为多个，所述目标参数还包括每个半径对应的体积占比，所述体积占比为同一半径对应的固相颗粒的体积与固相颗粒总体积的比值，所述固相颗粒总体积根据所述孔隙度确定；
所述颗粒添加单元，用于：
按照多个半径的排列顺序，遍历所述多个半径；
根据第一个半径，在所述目标三维区域中添加具有所述第一个半径的球体固相颗粒，直至所述目标三维区域中所述第一个半径对应的固相颗粒的体积与所述固相颗粒总体积的比值，等于所述第一个半径对应的体积占比；
根据下一个半径，在所述目标三维区域中添加具有所述下一个半径的球体固相颗粒，直至所述目标三维区域中所述下一个半径对应的固相颗粒的体积与所述固相颗粒总体积的比值，等于所述下一个半径对应的体积占比。
在另一种可能实现方式中，所述第一构建模块，还包括：
第二标记添加单元，用于为所述目标三维区域中的每个网格添加第二标记，所述第二标记表示对应的网格为孔隙；
所述颗粒添加单元，用于在所述目标三维区域中随机指定一个坐标点；为以所述坐标点为球心，具有所述半径的球体区域内的网格添加所述第一标记。
在另一种可能实现方式中，所述裂缝添加单元，用于获取裂缝开度；根据所述裂缝开度，在所述第三模型中随机确定裂缝对应的网格，为所述裂缝对应的网格添加所述第三标记，所述第三标记表示对应的网格为裂缝。
在另一种可能实现方式中，所述第一模拟模块，用于：
对所述第一模型中的流体施加预设压力梯度，以使所述流体按照所述格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述格子-玻尔兹曼演化关系中的流体粒子分布函数的函数值；
根据所述函数值，确定所述第一模型中的流体速度场；
从所述流体速度场中获取裂缝出口端的流体速度，将所述裂缝出口端的流体速度之和，作为所述第一流量。
在另一种可能实现方式中，所述第二模拟模块，用于：
对所述第二模型中的流体施加所述预设压力梯度，以使所述流体按照所述格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述格子-玻尔兹曼演化关系中的流体粒子分布函数的函数值；
根据所述函数值，确定所述第二模型中的流体速度场；
从所述流体速度场中获取裂缝出口端的流体速度，将所述裂缝出口端的流体速度之和，作为所述第二流量。
另一方面，提供了一种计算机设备，所述计算机设备包括处理器和存储器，所述存储器中存储有至少一条程序代码，所述至少一条程序代码由所述处理器加载并执行以实现如上述方面所述的窜流量确定方法中所执行的操作。
另一方面，提供了一种计算机可读存储介质，所述计算机可读存储介质中存储有至少一条程序代码，所述至少一条程序代码由处理器加载并执行以实现如上述方面所述的窜流量确定方法中所执行的操作。
再一方面，提供了一种计算机程序产品或计算机程序，所述计算机程序产品或计算机程序包括计算机程序代码，所述计算机程序代码存储在计算机可读存储介质中。计算机设备的处理器从计算机可读存储介质读取该计算机程序代码，处理器执行所述计算机程序代码，使得所述计算机设备实现如上述方面所述的窜流量确定方法中所执行的操作。
本申请实施例提供的技术方案带来的有益效果至少包括：
本申请实施例提供的方法、装置、计算机设备及存储介质，通过分别模拟第一模型和第二模型中的流体渗流，确定第一模型和第二模型出口端裂缝的流量，由于第一模型为裂缝-孔隙双重介质模型，第二模型为裂缝模型，且第一模型与第二模型中的裂缝相同，因此通过将第一模型和第二模型出口端裂缝的流量进行对比，能够确定裂缝-孔隙双重介质的窜流量，以该窜流量来表示孔隙窜流对裂缝渗流能力的影响，提供了一种准确确定裂缝-孔隙双重介质的窜流量的方法，无需人工对试井曲线进行分析，节省了人力成本，提高了效率，而且避免了人工操作中主观因素的影响，提供了准确的窜流量数据，有助于操作人员客观准确地认识油气井的产能情况，为后续工作提供了数据支持。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请实施例的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本申请实施例提供的一种窜流量确定方法的流程图；
图2是本申请实施例提供的一种窜流量确定方法的流程图；
图3是本申请实施例提供的一种孔隙模型的示意图；
图4是本申请实施例提供的一种固相颗粒的多个半径对应的体积占比的示意图；
图5是本申请实施例提供的一种裂缝-孔隙双重介质模型的示意图；
图6是本申请实施例提供的一种裂缝出口端的流体速度的示意图；
图7是本申请实施例提供的一种窜流量确定装置的结构示意图；
图8是本申请实施例提供的一种终端的结构示意图；
图9是本申请实施例提供的一种服务器的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合附图对本申请实施方式作进一步地详细描述。
本申请所使用的术语“第一”、“第二”、“第三”等可在本文中用于描述各种概念，但除非特别说明，这些概念不受这些术语限制。这些术语仅用于将一个概念与另一个概念区分。举例来说，在不脱离本申请的范围的情况下，可以将第一模型称为第二模型，且类似地，可将第二模型称为第一模型。
本申请所使用的术语“多个”、“每个”，多个包括两个或两个以上，而每个是指对应的多个中的每一个。举例来说，多个半径包括7个半径，而每个是指这7个半径中的每一个半径。
图1是本申请实施例提供的一种窜流量确定方法的流程图，如图1所示，该方法包括：
101、构建第一模型。
其中，第一模型为裂缝-孔隙双重介质模型。
102、模拟第一模型中的流体渗流，直至第一模型中的流体处于稳定状态，确定第一模型的出口端裂缝的流量，作为第一流量。
103、构建第二模型，第二模型为裂缝模型，且第一模型与第二模型中的裂缝相同。
104、模拟第二模型中的流体渗流，直至第二模型中的流体处于稳定状态，确定第二模型的出口端裂缝的流量，作为第二流量。
105、将第一流量与第二流量之间的差值，确定为裂缝-孔隙双重介质的窜流量。
本申请实施例提供的方法，通过分别模拟第一模型和第二模型中的流体渗流，确定第一模型和第二模型出口端裂缝的流量，由于第一模型为裂缝-孔隙双重介质模型，第二模型为裂缝模型，且第一模型与第二模型中的裂缝相同，因此通过将第一模型和第二模型出口端裂缝的流量进行对比，能够确定裂缝-孔隙双重介质的窜流量，以该窜流量来表示孔隙窜流对裂缝渗流能力的影响，提供了一种准确确定裂缝-孔隙双重介质的窜流量的方法，该方法无需人工对试井曲线进行分析，节省了人力成本，提高了效率，而且避免了人工操作中主观因素的影响，提供了准确的窜流量数据，有助于操作人员客观准确地认识油气井的产能情况，为后续工作提供了数据支持。
图2是本申请实施例提供的另一种窜流量确定方法的流程图，应用于计算机设备，该计算机设备为终端或服务器，如图2所示，该方法包括：
201、在目标三维区域中添加固相颗粒，以使添加的固相颗粒之间构成孔隙，得到第三模型，第三模型为孔隙模型。
在本申请实施例中，为了确定窜流量，需要先构建第一模型和第二模型，第一模型为裂缝-孔隙双重介质模型，第二模型为裂缝模型，且第一模型与第二模型中的裂缝相同，以便后续将第一模型与第二模型中的裂缝出口端的流量进行对比，从而确定窜流量，以该窜流量来表示孔隙窜流对裂缝渗流能力的影响。
首先对第一模型的构建过程进行说明，在构建第一模型的过程中，会先创建目标三维区域，该目标三维区域是用来构造模型的三维构造区域，该目标三维区域包括多个网格，由多个网格拼接而成。
例如，该目标三维区域具有X、Y和Z三个方向，在X方向上的网格数为A，在Y方向上的网格数为B，在Z方向上的网格数为C。通过确定三个方向上的网格数，能够将相应的网格排列起来，得到该目标三维区域。
之后即可在该目标三维区域中添加固相颗粒，从而形成第三模型，并且添加的固相颗粒之间构成孔隙，因此该第三模型为孔隙模型。
可选地，该步骤201包括以下步骤301-302：
301、获取目标参数，目标参数包括固相颗粒的半径及孔隙度。
在本申请实施例中，所添加的固相颗粒为球体颗粒，因此可以设置固相颗粒的半径，根据该半径来添加具有该半径的球体固相颗粒。孔隙度指示第三模型中的孔隙空间体积之和与第三模型体积的比值，因此能够指示第三模型中的孔隙空间所占的比例。
可选地，考虑到一般的碳酸盐岩中会包括不同半径的固相颗粒，因此所获取的目标参数中可以包括多个半径，以便根据该多个半径构建出的第三模型中能够包括不同半径的固相颗粒，增强了准确性。
可选地，目标参数还包括每个半径对应的体积占比，体积占比为同一半径对应的固相颗粒的体积与固相颗粒总体积的比值，固相颗粒总体积根据孔隙度确定。例如，根据第三模型的总体积和孔隙度，确定第三模型中的孔隙空间对应的总体积以及固相颗粒的总体积。之后根据每个半径对应的体积占比和该固相颗粒的总体积，确定每个半径对应的固相颗粒的总体积。
可选地，固相颗粒的半径及孔隙度可以由操作人员设置，也可以采用其他方式设置。
可选地，固相颗粒的半径及孔隙度可以根据一般的碳酸盐岩中的固相颗粒的半径及孔隙度确定，从而以碳酸盐岩作为依据，使第三模型更接近于真实的碳酸盐岩，增强了第三模型的准确性，进而增强了后续确定的窜流量的准确性。
302、根据该半径，在目标三维区域中添加具有该半径的球体固相颗粒，直至目标三维区域中的孔隙空间体积之和与目标三维区域体积的比值等于孔隙度时，将目标三维区域作为第三模型。
针对于一个半径，每次可以在目标三维区域中添加具有该半径的球体固相颗粒，每次所添加的球体固相颗粒的数量可以任意设置，可以为1或者更多。并且，所添加的球体固相颗粒可以位于目标三维区域的任意位置，例如可以采用随机算法来随机确定球体固相颗粒的位置，从而在目标三维区域中随机添加固相颗粒。
通过随机添加固相颗粒，使固相颗粒和孔隙的分布更接近碳酸盐岩，增强了第三模型的准确性，进而增强了后续确定的窜流量的准确性。
另外，该目标三维区域中，除球体固相颗粒之外的其他区域可以作为孔隙，则在添加球体固相颗粒的过程中，可以统计目标三维区域中的孔隙空间体积之和，即为孔隙的总体积，确定该孔隙的总体积与该目标三维区域的体积的比值，判断该比值是否等于该孔隙度，如果该比值大于该孔隙度，表示目前的目标三维区域中的孔隙过多，则继续在目标三维区域中添加球体固相颗粒，直至该比值等于该孔隙度时，停止添加，将当前的目标三维区域作为第三模型。
如图3所示，正方体区域为目标三维区域，通过在目标三维区域中添加固相颗粒301，以使固相颗粒301之间构成孔隙302，从而得到第三模型。
可选地，在固相颗粒的半径为多个，目标参数还包括每个半径对应的体积占比的情况下，该步骤303可以包括：
按照多个半径的排列顺序，遍历多个半径；根据第一个半径，在目标三维区域中添加具有第一个半径的球体固相颗粒，直至目标三维区域中第一个半径对应的固相颗粒的体积与固相颗粒总体积的比值，等于第一个半径对应的体积占比；根据下一个半径，在目标三维区域中添加具有下一个半径的球体固相颗粒，直至目标三维区域中下一个半径对应的固相颗粒的体积与固相颗粒总体积的比值，等于下一个半径对应的体积占比，以此类推，直至将该多个半径对应的固相颗粒添加完成为止，将当前的目标三维区域作为第三模型。
可选地，按照多个半径的排列顺序遍历多个半径，可以是按照半径从小到大的排列顺序遍历多个半径，也可以是其他排列顺序。例如，第三模型每个方向的网格数都为80，孔隙度εm为0.07，固相颗粒的多个半径及每个半径对应的体积占比如图4所示。
根据如图4所示的多个固相颗粒的半径和每个半径对应的体积占比，在目标三维区域中添加具有半径的球体固相颗粒，直至目标三维区域中的孔隙空间体积之和与目标三维区域体积的比值等于孔隙度时，将目标三维区域作为第三模型。之后，按照3μm(微米)、5μm、7μm、9μm、11μm、13μm、15μm的排列顺序，遍历七个半径。首先，在目标三维区域中添加具有3μm半径的球体固相颗粒，直至目标三维区域中3μm半径对应的固相颗粒的体积与固相颗粒总体积的比值，等于10％；之后在目标三维区域中添加具有5μm半径的球体固相颗粒，直至目标三维区域中5μm半径对应的固相颗粒的体积与固相颗粒总体积的比值，等于14％。
可选地，在步骤303之前，该方法还包括：为目标三维区域中的每个网格添加第二标记，第二标记表示对应的网格为孔隙。通过这种添加第二标记的方式，能够先将目标三维区域中的网格设置为孔隙。
相应地，该步骤303包括：在目标三维区域中随机指定一个坐标点；为以该坐标点为球心，具有半径的球体区域内的网格添加第一标记，该第一标记表示对应的网格为固相颗粒，也即是将该球体区域设置为了固相颗粒，从而实现了在目标三维区域中添加固相颗粒的目的。
例如，在目标三维区域中指定一个坐标点(x，y，z)，在((i－x)2+(j－y)2+(k－z)2)0.5≤r的范围内的网格添加第一标记。其中，(x，y，z)表示目标三维区域中任一个点的三个方向上的坐标，(i，j，k)表示添加第一标记的网格的坐标，r表示固相颗粒的半径。
其中，第一标记与第二标记为不同的标记，例如第一标记为1，第二标记为0。例如，在目标三维区域中为每个网格确定对应的像素值，按照每个网格对应的像素值来显示目标三维区域，在显示目标三维区域时，由于像素值的不同，会导致不同的网格显示出不同的颜色。因此将像素值看做是标记，将添加固相颗粒的区域的像素值设置为1，将孔隙的区域的像素值设置为0，或者，将添加固相颗粒的区域的像素值设置为0，将孔隙的区域的像素值设置为1，不仅起到了区分孔隙和固相颗粒的目的，而且还能通过显示出的颜色直观地确定哪部分区域是孔隙，哪部分区域是固相颗粒。
另外，在固相颗粒的半径包括多个的情况下，在为目标三维区域中的每个网格添加第二标记之后，对于每个半径，在目标三维区域中随机指定一个坐标点；为以该坐标点为球心，具有该半径的球体区域内的网格添加第一标记，直至目标三维区域中第一个半径对应的固相颗粒的体积与固相颗粒总体积的比值，等于第一个半径对应的体积占比；根据下一个半径，在目标三维区域中随机指定一个坐标点，为以该坐标点为球心，具有该下一个半径的球体区域内的网格添加第一标记，直至目标三维区域中下一个半径对应的固相颗粒的体积与固相颗粒总体积的比值，等于下一个半径对应的体积占比，以此类推，直至将该多个半径对应的固相颗粒添加完成为止，将当前的目标三维区域作为第三模型。
如图4所示，按照半径从小到大的顺序将固相颗粒随机分布于构造区域，首先按照最小的半径3μm，在目标三维区域范围内运用随机函数随机指定一个球形孔隙的中心坐标，将以该中心坐标为球心、半径为3μm的范围内的网格添加像素值为1的标记，表示该范围为固相颗粒，重复该过程直至半径为3μm固相颗粒所占体积比重达到15％，随后开始5μm、7μm、9μm、11μm、13μm、15μm固相颗粒的生成。
202、对第三模型中的流体施加预设压力梯度，以使流体按照格子-玻尔兹曼演化关系进行演化，直至流体处于稳定状态。
通过采用格子-玻尔兹曼方法模拟第三模型中的流体渗流，能够获得流体处于稳定状态的第三模型，避免了流体渗流对后续处理过程的影响，提高了准确性。
其中，模拟第三模型中的流体渗流的具体方式，与后续的第一模型和第二模型类似，在此暂不做说明。
203、在第三模型中添加裂缝，得到第一模型。
可选地，该步骤203包括：获取裂缝开度，根据裂缝开度，在第三模型中随机确定裂缝对应的网格，为裂缝对应的网格添加第三标记，第三标记表示对应的网格为裂缝，从而实现了在第三模型中添加裂缝。例如裂缝开度为20μm的，在第三模型设置开度为20μm的水平裂缝，对裂缝对应的网格添加像素值为2的标记，得到第一模型。
可选地，裂缝开度可以由操作人员设置，也可以采用其他方式设置。另外，裂缝的走向也可以由操作人员设置，也可以采用其他方式设置。
可选地，第三标记可以指代网格的像素值，第三标记与第一标记和第二标记均不同，例如第三标记可以为2。
由于第三模型为孔隙模型，在第三模型中添加裂缝，得到的第一模型即为裂缝-孔隙双重介质模型。该第一模型可以如图5所示，该第一模型不仅包括固相颗粒301和孔隙302，还包括裂缝303。
204、模拟第一模型中的流体渗流，直至第一模型中的流体处于稳定状态，确定第一模型的出口端裂缝的流量，作为第一流量。
可选地，对第一模型中的流体施加预设压力梯度，以使流体按照格子-玻尔兹曼演化关系进行演化，直至流体处于稳定状态，确定格子-玻尔兹曼演化关系中的流体粒子分布函数的函数值，根据该函数值，确定第一模型中的流体速度场，从流体速度场中获取裂缝出口端的流体速度，将裂缝出口端的流体速度之和，作为第一流量。
例如，格子-玻尔兹曼演化关系如下：
fi(x+eiδt,t+δt)-fi(x,t)＝-Ω[fi(x,t)-fi eq(x,t)]；
其中，fi为粒子分布函数，单位为kg/m3(千克/平方米)；x表示第一模型中的网格的坐标，单位为m(米)，ei为格子速度，单位为m/s(米/秒)；t为时间，单位为s(秒)；δt为时间差，单位为s；Ω为碰撞矩阵，Ω＝M-1SM，M为粒子分布函数fi的变换矩阵，S表示对角函数；fi eq为稳定状态下的粒子分布函数。
对第一模型中的流体施加预设压力梯度
Figure BDA0002702016180000145
以使流体按照上述格子-玻尔兹曼演化关系进行演化，直至流体处于稳定状态，确定格子-玻尔兹曼演化关系中的流体粒子分布函数的函数值，根据该函数值，采用以下公式，确定第一模型的流体速度场：
Figure BDA0002702016180000141
其中，fi eq为稳定状态下的粒子分布函数值；ρ为流体密度，单位为kg/m3；wi为权系数；ei为格子速度，单位为m/s；u为流体速度，单位为m/s；cs为声速，单位为m/s。
另外，还可以采用以下公式，确定第一模型的渗透率，渗透率能够表示流体在孔隙中流动的难易程度：
Figure BDA0002702016180000142
其中，km为第一模型的渗透率；μ为流体动力粘度，单位为Pa·s(帕斯卡·秒)；
Figure BDA0002702016180000143
为第一模型水平方向上的流体速度ux的平均值，单位为m/s；
Figure BDA0002702016180000144
为预设压力梯度。例如，使用水作为模拟流体，流体密度ρ为983.1kg/m3，流体动力粘度μ为0.0004699Pa·s，计算得到模型的渗透率km为2.20×10-15m2。
通过确定模型的渗透率，体现了油气通过上述模型模拟的地层岩石流向井底的能力，为确定地层岩石的油气生产能力提供了数据支持。
可选地，第一模型中的流体速度场可以包括第一模型中各个位置处的流体速度。
如图6所示，图6示出了在裂缝出口端的流体速度，从图6可以看出在宽度为40-65、高度为22-80的区间内流体速度突然增大，表示在该宽度区间内的裂缝存在孔隙，受到孔隙流体窜流的影响而产生了窜流量，导致流体速度突然增大。
205、为第一模型中的孔隙对应的网格添加第一标记，得到第二模型，第一标记表示对应的网格为固相颗粒。
其中，第一模型包括固相颗粒、固相颗粒之间的孔隙和裂缝，且第一模型由网格构成。通过为第一模型中的孔隙对应的网格添加第一标记，替换掉该网格中原有的第二标记，以便将这些网格从孔隙更改为固相颗粒，从而得到了包含固相颗粒和裂缝，而不包含孔隙的第二模型，因此该第二模型为裂缝模型。
206、模拟第二模型中的流体渗流，直至第二模型中的流体处于稳定状态，确定第二模型的出口端裂缝的流量，作为第二流量。
可选地，对第二模型中的流体施加预设压力梯度，以使流体按照格子-玻尔兹曼演化关系进行演化，直至流体处于稳定状态，确定格子-玻尔兹曼演化关系中的流体粒子分布函数的函数值，根据该函数值，确定第二模型中的流体速度场；从流体速度场中获取裂缝出口端的流体速度，将裂缝出口端的流体速度之和，作为第二流量。
其中，获取第二模型对应的第二流量的过程，与上述获取第一流量的过程类似，在此不再赘述。另外本申请实施例也能够获取第二模型的渗透率，获取过程与上述获取第一模型的渗透率的过程类似，在此也不再赘述。
207、将第一流量与第二流量之间的差值，确定为裂缝-孔隙双重介质的窜流量。
通过确定裂缝-孔隙双重介质的窜流量，反映了第一模型的孔隙中的流体向裂缝窜流能力的大小，进而反映了第一模型模拟的碳酸盐岩的油气藏生产能力。
例如，采用格子-玻尔兹曼方法模拟裂缝-孔隙双重介质模型中的流体渗流，处于稳定状态后，确定出口端裂缝的第一流量q1为28.9107m/s。之后采用类似方法，确定裂缝模型的出口端裂缝的第二流量q2为28.6675m/s，则窜流量qex为0.2432m/s。
本申请实施例提供的方法，通过分别模拟第一模型和第二模型中的流体渗流，确定第一模型和第二模型出口端裂缝的流量，由于第一模型为裂缝-孔隙双重介质模型，第二模型为裂缝模型，且第一模型与第二模型中的裂缝相同，因此通过将第一模型和第二模型出口端裂缝的流量进行对比，能够确定裂缝-孔隙双重介质的窜流量，以该窜流量来表示孔隙窜流对裂缝渗流能力的影响，提供了一种准确确定裂缝-孔隙双重介质的窜流量的方法，该方法无需人工对试井曲线进行分析，节省了人力成本，提高了效率，而且避免了人工操作中主观因素的影响，提供了准确的窜流量数据，有助于操作人员客观准确地认识油气井的产能情况，为后续工作提供了数据支持。
图7是本申请实施例提供的一种窜流量确定装置的结构示意图，如图7所示，该装置包括：
第一构建模块701，用于构建第一模型，第一模型为裂缝-孔隙双重介质模型；
第一模拟模块702，用于模拟第一模型中的流体渗流，直至第一模型中的流体处于稳定状态，确定第一模型的出口端裂缝的流量，作为第一流量；
第二构建模块703，用于构建第二模型，第二模型为裂缝模型，且第一模型与第二模型中的裂缝相同；
第二模拟模块704，用于模拟第二模型中的流体渗流，直至第二模型中的流体处于稳定状态，确定第二模型的出口端裂缝的流量，作为第二流量；
确定模块705，用于将第一流量与第二流量之间的差值，确定为裂缝-孔隙双重介质的窜流量。
本申请实施例提供的装置，通过分别模拟第一模型和第二模型中的流体渗流，确定第一模型和第二模型出口端裂缝的流量，由于第一模型为裂缝-孔隙双重介质模型，第二模型为裂缝模型，且第一模型与第二模型中的裂缝相同，因此通过将第一模型和第二模型出口端裂缝的流量进行对比，能够确定裂缝-孔隙双重介质的窜流量，以该窜流量来表示孔隙窜流对裂缝渗流能力的影响，提供了一种准确确定裂缝-孔隙双重介质的窜流量的方法，该方法无需人工对试井曲线进行分析，节省了人力成本，提高了效率，而且避免了人工操作中主观因素的影响，提供了准确的窜流量数据，有助于操作人员客观准确地认识油气井的产能情况，为后续工作提供了数据支持。
在一种可能实现方式中，第一模拟模块702，用于：
对第一模型中的流体施加预设压力梯度，以使流体按照格子-玻尔兹曼演化关系进行演化，直至流体处于稳定状态，确定第一模型的出口端裂缝的流量，作为第一流量。
在另一种可能实现方式中，第二模拟模块704，用于：
对第二模型中的流体施加预设压力梯度，以使流体按照格子-玻尔兹曼演化关系进行演化，直至流体处于稳定状态，确定第二模型的出口端裂缝的流量，作为第二流量。
在另一种可能实现方式中，第一模型包括固相颗粒、固相颗粒之间的孔隙和裂缝，且第一模型由网格构成；第二构建模块703，包括：
第一标记添加单元，用于为第一模型中的孔隙对应的网格添加第一标记，得到第二模型，第一标记表示对应的网格为固相颗粒。
在另一种可能实现方式中，第一构建模块701，包括：
颗粒添加单元，用于在目标三维区域中添加固相颗粒，以使添加的固相颗粒之间构成孔隙，得到第三模型，第三模型为孔隙模型；
流体演化单元，用于对第三模型中的流体施加预设压力梯度，以使流体按照格子-玻尔兹曼演化关系进行演化，直至流体处于稳定状态；
裂缝添加单元，用于在第三模型中添加裂缝，得到第一模型。
在另一种可能实现方式中，颗粒添加单元，用于：
获取目标参数，目标参数包括固相颗粒的半径及孔隙度，孔隙度指示第三模型中的孔隙空间体积之和与第三模型体积的比值；
根据半径，在目标三维区域中添加具有半径的球体固相颗粒，直至目标三维区域中的孔隙空间体积之和与目标三维区域体积的比值等于孔隙度时，将目标三维区域作为第三模型。
在另一种可能实现方式中，固相颗粒的半径为多个，目标参数还包括每个半径对应的体积占比，体积占比为同一半径对应的固相颗粒的体积与固相颗粒总体积的比值，固相颗粒总体积根据孔隙度确定；
颗粒添加单元，用于：
按照多个半径的排列顺序，遍历多个半径；
根据第一个半径，在目标三维区域中添加具有第一个半径的球体固相颗粒，直至目标三维区域中第一个半径对应的固相颗粒的体积与固相颗粒总体积的比值，等于第一个半径对应的体积占比；
根据下一个半径，在目标三维区域中添加具有下一个半径的球体固相颗粒，直至目标三维区域中下一个半径对应的固相颗粒的体积与固相颗粒总体积的比值，等于下一个半径对应的体积占比。
在另一种可能实现方式中，第一构建模块701，还包括：
第二标记添加单元，用于为目标三维区域中的每个网格添加第二标记，第二标记表示对应的网格为孔隙；
颗粒添加单元，用于在目标三维区域中随机指定一个坐标点；为以坐标点为球心，具有半径的球体区域内的网格添加第一标记。
在另一种可能实现方式中，裂缝添加单元，用于获取裂缝开度；根据裂缝开度，在第三模型中随机确定裂缝对应的网格，为裂缝对应的网格添加第三标记，第三标记表示对应的网格为裂缝。
在另一种可能实现方式中，第一模拟模块702，用于：
对第一模型中的流体施加预设压力梯度，以使流体按照格子-玻尔兹曼演化关系进行演化，直至流体处于稳定状态，确定格子-玻尔兹曼演化关系中的流体粒子分布函数的函数值；
根据函数值，确定第一模型中的流体速度场；
从流体速度场中获取裂缝出口端的流体速度，将裂缝出口端的流体速度之和，作为第一流量。
在另一种可能实现方式中，第二模拟模块704，用于：
对第二模型中的流体施加预设压力梯度，以使流体按照格子-玻尔兹曼演化关系进行演化，直至流体处于稳定状态，确定格子-玻尔兹曼演化关系中的流体粒子分布函数的函数值；
根据函数值，确定第二模型中的流体速度场；
从流体速度场中获取裂缝出口端的流体速度，将裂缝出口端的流体速度之和，作为第二流量。
图8示出了本申请一个示例性实施例提供的终端800的结构框图。该终端800可以是便携式移动终端，比如：智能手机、平板电脑、MP3播放器(Moving Picture ExpertsGroup Audio Layer III，动态影像专家压缩标准音频层面3)、MP4(Moving PictureExperts Group Audio Layer IV，动态影像专家压缩标准音频层面4)播放器、笔记本电脑或台式电脑。终端800还可能被称为用户设备、便携式终端、膝上型终端、台式终端等其他名称。
终端800包括有：处理器801和存储器802。
处理器801可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器801可以采用DSP(Digital Signal Processing，数字信号处理)、FPGA(Field-Programmable Gate Array，现场可编程门阵列)、PLA(Programmable Logic Array，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器801也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称CPU(Central ProcessingUnit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器801可以集成有GPU(Graphics Processing Unit，图像处理器)，GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器801还可以包括AI(Artificial Intelligence，人工智能)处理器，该AI处理器用于处理有关机器学习的计算操作。
存储器802可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器802还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器802中的非暂态的计算机可读存储介质用于存储至少一个程序代码，该至少一个程序代码用于被处理器801所执行以实现本申请中方法实施例提供的窜流量确定方法。
在一些实施例中，终端800还可选包括有：外围设备接口803和至少一个外围设备。处理器801、存储器802和外围设备接口803之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口803相连。具体地，外围设备包括：射频电路804、显示屏805、摄像头组件806、音频电路807、定位组件808和电源809中的至少一种。
外围设备接口803可被用于将I/O(Input/Output，输入/输出)相关的至少一个外围设备连接到处理器801和存储器802。在一些实施例中，处理器801、存储器802和外围设备接口803被集成在同一芯片或电路板上；在一些其他实施例中，处理器801、存储器802和外围设备接口803中的任意一个或两个可以在单独的芯片或电路板上实现，本实施例对此不加以限定。
射频电路804用于接收和发射RF(Radio Frequency，射频)信号，也称电磁信号。射频电路804通过电磁信号与通信网络以及其他通信设备进行通信。射频电路804将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转换为电信号。可选地，射频电路804包括：天线系统、RF收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路804可以通过至少一种无线通信协议来与其它终端进行通信。该无线通信协议包括但不限于：万维网、城域网、内联网、各代移动通信网络(2G、3G、4G及5G)、无线局域网和/或WiFi(Wireless Fidelity，无线保真)网络。在一些实施例中，射频电路804还可以包括NFC(Near Field Communication，近距离无线通信)有关的电路，本申请对此不加以限定。
显示屏805用于显示UI(User Interface，用户界面)。该UI可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏805是触摸显示屏时，显示屏805还具有采集在显示屏805的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器801进行处理。此时，显示屏805还可以用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中，显示屏805可以为一个，设置在终端800的前面板；在另一些实施例中，显示屏805可以为至少两个，分别设置在终端800的不同表面或呈折叠设计；在另一些实施例中，显示屏805可以是柔性显示屏，设置在终端800的弯曲表面上或折叠面上。甚至，显示屏805还可以设置成非矩形的不规则图形，也即异形屏。显示屏805可以采用LCD(Liquid Crystal Display，液晶显示屏)、OLED(Organic Light-Emitting Diode,有机发光二极管)等材质制备。
摄像头组件806用于采集图像或视频。可选地，摄像头组件806包括前置摄像头和后置摄像头。前置摄像头设置在终端的前面板，后置摄像头设置在终端的背面。在一些实施例中，后置摄像头为至少两个，分别为主摄像头、景深摄像头、广角摄像头、长焦摄像头中的任意一种，以实现主摄像头和景深摄像头融合实现背景虚化功能、主摄像头和广角摄像头融合实现全景拍摄以及VR(Virtual Reality，虚拟现实)拍摄功能或者其它融合拍摄功能。在一些实施例中，摄像头组件806还可以包括闪光灯。闪光灯可以是单色温闪光灯，也可以是双色温闪光灯。双色温闪光灯是指暖光闪光灯和冷光闪光灯的组合，可以用于不同色温下的光线补偿。
音频电路807可以包括麦克风和扬声器。麦克风用于采集用户及环境的声波，并将声波转换为电信号输入至处理器801进行处理，或者输入至射频电路804以实现语音通信。出于立体声采集或降噪的目的，麦克风可以为多个，分别设置在终端800的不同部位。麦克风还可以是阵列麦克风或全向采集型麦克风。扬声器则用于将来自处理器801或射频电路804的电信号转换为声波。扬声器可以是传统的薄膜扬声器，也可以是压电陶瓷扬声器。当扬声器是压电陶瓷扬声器时，不仅可以将电信号转换为人类可听见的声波，也可以将电信号转换为人类听不见的声波以进行测距等用途。在一些实施例中，音频电路807还可以包括耳机插孔。
定位组件808用于定位终端800的当前地理位置，以实现导航或LBS(LocationBased Service，基于位置的服务)。定位组件808可以是基于美国的GPS(GlobalPositioning System，全球定位系统)、中国的北斗系统或俄罗斯的伽利略系统的定位组件。
电源809用于为终端800中的各个组件进行供电。电源809可以是交流电、直流电、一次性电池或可充电电池。当电源809包括可充电电池时，该可充电电池可以是有线充电电池或无线充电电池。有线充电电池是通过有线线路充电的电池，无线充电电池是通过无线线圈充电的电池。该可充电电池还可以用于支持快充技术。
在一些实施例中，终端800还包括有一个或多个传感器810。该一个或多个传感器810包括但不限于：加速度传感器811、陀螺仪传感器812、压力传感器88、指纹传感器814、光学传感器815以及接近传感器816。
加速度传感器811可以检测以终端800建立的坐标系的三个坐标轴上的加速度大小。比如，加速度传感器811可以用于检测重力加速度在三个坐标轴上的分量。处理器801可以根据加速度传感器811采集的重力加速度信号，控制显示屏805以横向视图或纵向视图进行用户界面的显示。加速度传感器811还可以用于游戏或者用户的运动数据的采集。
陀螺仪传感器812可以检测终端800的机体方向及转动角度，陀螺仪传感器812可以与加速度传感器811协同采集用户对终端800的3D动作。处理器801根据陀螺仪传感器812采集的数据，可以实现如下功能：动作感应(比如根据用户的倾斜操作来改变UI)、拍摄时的图像稳定、游戏控制以及惯性导航。
压力传感器88可以设置在终端800的侧边框和/或显示屏805的下层。当压力传感器88设置在终端800的侧边框时，可以检测用户对终端800的握持信号，由处理器801根据压力传感器88采集的握持信号进行左右手识别或快捷操作。当压力传感器88设置在显示屏805的下层时，由处理器801根据用户对显示屏805的压力操作，实现对UI界面上的可操作性控件进行控制。可操作性控件包括按钮控件、滚动条控件、图标控件、菜单控件中的至少一种。
指纹传感器814用于采集用户的指纹，由处理器801根据指纹传感器814采集到的指纹识别用户的身份，或者，由指纹传感器814根据采集到的指纹识别用户的身份。在识别出用户的身份为可信身份时，由处理器801授权该用户执行相关的敏感操作，该敏感操作包括解锁屏幕、查看加密信息、下载软件、支付及更改设置等。指纹传感器814可以被设置在终端800的正面、背面或侧面。当终端800上设置有物理按键或厂商Logo时，指纹传感器814可以与物理按键或厂商Logo集成在一起。
光学传感器815用于采集环境光强度。在一个实施例中，处理器801可以根据光学传感器815采集的环境光强度，控制显示屏805的显示亮度。具体地，当环境光强度较高时，调高显示屏805的显示亮度；当环境光强度较低时，调低显示屏805的显示亮度。在另一个实施例中，处理器801还可以根据光学传感器815采集的环境光强度，动态调整摄像头组件806的拍摄参数。
接近传感器816，也称距离传感器，设置在终端800的前面板。接近传感器816用于采集用户与终端800的正面之间的距离。在一个实施例中，当接近传感器816检测到用户与终端800的正面之间的距离逐渐变小时，由处理器801控制显示屏805从亮屏状态切换为息屏状态；当接近传感器816检测到用户与终端800的正面之间的距离逐渐变大时，由处理器801控制显示屏805从息屏状态切换为亮屏状态。
本领域技术人员可以理解，图8中示出的结构并不构成对终端800的限定，可以包括比图示更多或更少的组件，或者组合某些组件，或者采用不同的组件布置。
图9是本申请实施例提供的一种服务器的结构示意图，该服务器900可因配置或性能不同而产生比较大的差异，可以包括一个或一个以上处理器(Central ProcessingUnits，CPU)901和一个或一个以上的存储器902，其中，存储器902中存储有至少一条程序代码，至少一条程序代码由处理器901加载并执行以实现上述各个方法实施例提供的方法。当然，该服务器还可以具有有线或无线网络接口、键盘及输入输出接口等部件，以便进行输入输出，该服务器还可以包括其他用于实现设备功能的部件，在此不做赘述。
服务器900可以用于执行上述窜流量确定方法中计算机设备所执行的步骤。
本申请实施例还提供了一种计算机设备，该计算机设备包括处理器和存储器，该存储器中存储有至少一条程序代码，该至少一条程序代码由处理器加载并执行以实现上述实施例的窜流量确定方法中所执行的操作。
本申请实施例还提供了一种计算机可读存储介质，该计算机可读存储介质中存储有至少一条程序代码，该至少一条程序代码由处理器加载并执行以实现上述实施例的窜流量确定方法中所执行的操作。
本申请实施例还提供了一种计算机程序产品或计算机程序，该计算机程序产品或计算机程序包括计算机程序代码，该计算机程序代码存储在计算机可读存储介质中。计算机设备的处理器从计算机可读存储介质读取该计算机程序代码，处理器执行该计算机程序代码，使得该计算机设备实现如上述实施例的窜流量确定方法中所执行的操作。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上所述仅为本申请实施例的可选实施例，并不用以限制本申请实施例，凡在本申请实施例的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。

Claims (14)

1.一种窜流量确定方法，其特征在于，所述方法包括：
构建第一模型，所述第一模型为裂缝-孔隙双重介质模型；
模拟所述第一模型中的流体渗流，直至所述第一模型中的流体处于稳定状态，确定所述第一模型的出口端裂缝的流量，作为第一流量；
构建第二模型，所述第二模型为裂缝模型，且所述第一模型与所述第二模型中的裂缝相同；
模拟所述第二模型中的流体渗流，直至所述第二模型中的流体处于稳定状态，确定所述第二模型的出口端裂缝的流量，作为第二流量；
将所述第一流量与所述第二流量之间的差值，确定为所述裂缝-孔隙双重介质的窜流量。
2.根据权利要求1所述的方法，其特征在于，所述模拟所述第一模型中的流体渗流，直至所述第一模型中的流体处于稳定状态，确定所述第一模型的出口端裂缝的流量，作为第一流量，包括：
对所述第一模型中的流体施加预设压力梯度，以使所述流体按照格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述第一模型的出口端裂缝的流量，作为所述第一流量。
3.根据权利要求2所述的方法，其特征在于，所述模拟所述第二模型中的流体渗流，直至所述第二模型中的流体处于稳定状态，确定所述第二模型的出口端裂缝的流量，作为第二流量，包括：
对所述第二模型中的流体施加所述预设压力梯度，以使所述流体按照所述格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述第二模型的出口端裂缝的流量，作为所述第二流量。
4.根据权利要求1所述的方法，其特征在于，所述第一模型包括固相颗粒、固相颗粒之间的孔隙和裂缝，且所述第一模型由网格构成；所述构建第二模型，包括：
为所述第一模型中的孔隙对应的网格添加第一标记，得到所述第二模型，所述第一标记表示对应的网格为固相颗粒。
5.根据权利要求1所述的方法，其特征在于，所述构建所述第一模型，包括：
在目标三维区域中添加所述固相颗粒，以使添加的固相颗粒之间构成孔隙，得到第三模型，所述第三模型为孔隙模型；
对所述第三模型中的流体施加预设压力梯度，以使所述流体按照格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态；
在所述第三模型中添加裂缝，得到所述第一模型。
6.根据权利要求5所述的方法，其特征在于，所述在目标三维区域中添加固相颗粒，以使固相颗粒之间构成孔隙，得到第三模型，包括：
获取目标参数，所述目标参数包括固相颗粒的半径及孔隙度，所述孔隙度指示所述第三模型中的孔隙空间体积之和与所述第三模型体积的比值；
根据所述半径，在所述目标三维区域中添加具有所述半径的球体固相颗粒，直至所述目标三维区域中的孔隙空间体积之和与所述目标三维区域体积的比值等于所述孔隙度时，将所述目标三维区域作为所述第三模型。
7.根据权利要求6所述的方法，其特征在于，所述固相颗粒的半径为多个，所述目标参数还包括每个半径对应的体积占比，所述体积占比为同一半径对应的固相颗粒的体积与固相颗粒总体积的比值，所述固相颗粒总体积根据所述孔隙度确定；
所述根据所述半径，在所述目标三维区域中添加具有所述半径的球体固相颗粒，直至所述目标三维区域中的孔隙空间体积之和与所述目标三维区域体积的比值等于所述孔隙度时，将所述目标三维区域作为所述第三模型，包括：
按照多个半径的排列顺序，遍历所述多个半径；
根据第一个半径，在所述目标三维区域中添加具有所述第一个半径的球体固相颗粒，直至所述目标三维区域中所述第一个半径对应的固相颗粒的体积与所述固相颗粒总体积的比值，等于所述第一个半径对应的体积占比；
根据下一个半径，在所述目标三维区域中添加具有所述下一个半径的球体固相颗粒，直至所述目标三维区域中所述下一个半径对应的固相颗粒的体积与所述固相颗粒总体积的比值，等于所述下一个半径对应的体积占比。
8.根据权利要求6所述的方法，其特征在于，所述根据所述半径，在所述目标三维区域中添加具有所述半径的球体固相颗粒之前，所述方法还包括：
为所述目标三维区域中的每个网格添加第二标记，所述第二标记表示对应的网格为孔隙；
所述根据所述半径，在所述目标三维区域中添加具有所述半径的球体固相颗粒，包括：
在所述目标三维区域中随机指定一个坐标点；
为以所述坐标点为球心，具有所述半径的球体区域内的网格添加所述第一标记。
9.根据权利要求5所述的方法，其特征在于，所述在所述第三模型中添加裂缝，得到所述第一模型，包括：
获取裂缝开度；
根据所述裂缝开度，在所述第三模型中随机确定裂缝对应的网格；
为所述裂缝对应的网格添加所述第三标记，所述第三标记表示对应的网格为裂缝。
10.根据权利要求2所述的方法，其特征在于，所述对所述第一模型中的流体施加预设压力梯度，以使所述流体按照所述格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述第一模型的出口端裂缝的流量，作为所述第一流量，包括：
对所述第一模型中的流体施加预设压力梯度，以使所述流体按照所述格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述格子-玻尔兹曼演化关系中的流体粒子分布函数的函数值；
根据所述函数值，确定所述第一模型中的流体速度场；
从所述流体速度场中获取裂缝出口端的流体速度，将所述裂缝出口端的流体速度之和，作为所述第一流量。
11.根据权利要求3所述的方法，其特征在于，所述对所述第二模型中的流体施加预设压力梯度，以使所述流体按照所述格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述第二模型的出口端裂缝的流量，作为所述第二流量，包括：
对所述第二模型中的流体施加所述预设压力梯度，以使所述流体按照所述格子-玻尔兹曼演化关系进行演化，直至所述流体处于稳定状态，确定所述格子-玻尔兹曼演化关系中的流体粒子分布函数的函数值；
根据所述函数值，确定所述第二模型中的流体速度场；
从所述流体速度场中获取裂缝出口端的流体速度，将所述裂缝出口端的流体速度之和，作为所述第二流量。
12.一种窜流量确定装置，其特征在于，所述装置包括：
第一构建模块，用于构建第一模型，所述第一模型为裂缝-孔隙双重介质模型；
第一模拟模块，用于模拟所述第一模型中的流体渗流，直至所述第一模型中的流体处于稳定状态，确定所述第一模型的出口端裂缝的流量，作为第一流量；
第二构建模块，用于构建第二模型，所述第二模型为裂缝模型，且所述第一模型与所述第二模型中的裂缝相同；
第二模拟模块，用于模拟所述第二模型中的流体渗流，直至所述第二模型中的流体处于稳定状态，确定所述第二模型的出口端裂缝的流量，作为第二流量；
确定模块，用于将所述第一流量与所述第二流量之间的差值，确定为所述裂缝-孔隙双重介质的窜流量。
13.一种计算机设备，其特征在于，所述计算机设备包括处理器和存储器，所述存储器中存储有至少一条程序代码，所述至少一条程序代码由所述处理器加载并执行以实现如权利要求1-11任一项所述的窜流量确定方法中所执行的操作。
14.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有至少一条程序代码，所述至少一条程序代码由处理器加载并执行，以实现如权利要求1-11任一项所述的窜流量确定方法中所执行的操作。
CN202011025535.6A 2020-09-25 2020-09-25 窜流量确定方法、装置、计算机设备及存储介质 Active CN112329353B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011025535.6A CN112329353B (zh) 2020-09-25 2020-09-25 窜流量确定方法、装置、计算机设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011025535.6A CN112329353B (zh) 2020-09-25 2020-09-25 窜流量确定方法、装置、计算机设备及存储介质

Publications (2)

Publication Number Publication Date
CN112329353A true CN112329353A (zh) 2021-02-05
CN112329353B CN112329353B (zh) 2022-11-01

Family

ID=74304054

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011025535.6A Active CN112329353B (zh) 2020-09-25 2020-09-25 窜流量确定方法、装置、计算机设备及存储介质

Country Status (1)

Country Link
CN (1) CN112329353B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112949050A (zh) * 2021-02-09 2021-06-11 中国石油天然气股份有限公司 酸蚀参数的确定方法、装置、计算机设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106446396A (zh) * 2016-09-20 2017-02-22 中国石油大学(华东) 用于确定裂缝溶洞分布对储层渗透性影响的方法和装置
CN106837315A (zh) * 2015-12-03 2017-06-13 中国石油化工股份有限公司 裂缝性碳酸盐岩基质与裂缝耦合作用表征方法
CN109887083A (zh) * 2019-01-29 2019-06-14 中国石油集团测井有限公司西南分公司 一种裂缝-孔隙双重介质耦合渗透率模型的建立方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106837315A (zh) * 2015-12-03 2017-06-13 中国石油化工股份有限公司 裂缝性碳酸盐岩基质与裂缝耦合作用表征方法
CN106446396A (zh) * 2016-09-20 2017-02-22 中国石油大学(华东) 用于确定裂缝溶洞分布对储层渗透性影响的方法和装置
CN109887083A (zh) * 2019-01-29 2019-06-14 中国石油集团测井有限公司西南分公司 一种裂缝-孔隙双重介质耦合渗透率模型的建立方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李勇等: "碳酸盐岩油藏三重介质油藏数值模拟研究", 《石油天然气学报》 *
班凡生等: "裂缝-孔隙型双重介质油藏流体窜流规律", 《辽宁工程技术大学学报(自然科学版)》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112949050A (zh) * 2021-02-09 2021-06-11 中国石油天然气股份有限公司 酸蚀参数的确定方法、装置、计算机设备及存储介质
CN112949050B (zh) * 2021-02-09 2022-11-04 中国石油天然气股份有限公司 酸蚀参数的确定方法、装置、计算机设备及存储介质

Also Published As

Publication number Publication date
CN112329353B (zh) 2022-11-01

Similar Documents

Publication Publication Date Title
CN110064200B (zh) 基于虚拟环境的物体构建方法、装置及可读存储介质
CN111005715A (zh) 确定气井产量的方法、装置及存储介质
CN108764530B (zh) 油井抽油机的工作参数配置方法及装置
CN112308103A (zh) 生成训练样本的方法和装置
CN112329353B (zh) 窜流量确定方法、装置、计算机设备及存储介质
CN111255434A (zh) 气井的试井方法、装置和计算机存储介质
CN109833623B (zh) 基于虚拟环境的物体构建方法、装置及可读存储介质
CN111287718B (zh) 确定射孔参数的方法、装置及存储介质
CN111754564B (zh) 视频展示方法、装置、设备及存储介质
CN108733902B (zh) 油藏渗透率的确定方法、装置及存储介质
CN109281648B (zh) 确定油藏的合理井网密度的方法和装置
CN109281664B (zh) 水平井产油量的预测方法和装置
CN116071184A (zh) 油气层有效厚度下限确定方法、装置、设备及存储介质
CN113018865B (zh) 攀爬线生成方法、装置、计算机设备及存储介质
CN114280686B (zh) 碳酸盐岩储层的岩心物性分析方法及设备
CN110399688B (zh) 自动驾驶的环境工况确定方法、装置及存储介质
CN111274662B (zh) 对流连通老腔的筛选方法及装置
CN113240784A (zh) 图像处理方法、装置、终端及存储介质
CN113418848B (zh) 岩心的渗透参数的确定方法、装置、设备及存储介质
CN112182986B (zh) 确定面积波及信息的方法、装置、设备及存储介质
CN113689484B (zh) 深度信息的确定方法、装置、终端及存储介质
CN111375204B (zh) 显示数据获取方法、装置、设备及存储介质
CN110058306B (zh) 确定三维速度体的方法、装置及计算机存储介质
CN115478836A (zh) 孔隙参数的阈值信息确定方法、装置、设备及存储介质
CN115223665A (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
GR01 Patent grant
GR01 Patent grant