CN102124185B - 使用基于单元的流动模拟结果计算流线轨迹的系统和方法 - Google Patents

Abstract

本发明公开了一种使用有限差分模拟结果计算流线轨迹的系统和方法，可以使用其与其它流线技术以分析结果。

Description

使用基于单元的流动模拟结果计算流线轨迹的系统和方法

相关申请的交叉引用

本申请要求于2007年12月7日提交的美国临时专利申请第61/005,786号的优先权，将其通过引用并入于此。

关于联邦政府资助研究声明

不适用。

技术领域

本发明涉及用于计算流线轨迹(也称作流线)的系统和方法，可以用于分析基于单元的流动模拟结果。更具体地，本发明使用有限差分流动模拟器的结果计算流线轨迹，可以与其它流线技术一起使用以分析结果。

背景技术

在石油和天然气生产领域，当向油藏(reservoir)中钻了一个或更多井眼并且通过井眼从油藏中抽出液态烃(例如，石油和天然气)时，油气藏被开采了。在这种情况下，井眼一般被称作出油井眼或产油井。流体，通常是水，在远离产油井的一个或更多点被引入油藏，从而替代液态或气态烃，这引起其通过产油井从油藏排出，将水引入油藏的点通常称作注入孔(injector)。

在产油井喷出的液态和/或气态烃和在注入孔注入的水建立了油藏的液体流动模式，其可能受到重力的影响。油藏内的流体流动可以模拟为时间的函数，从而预测来自特定的产油井的液态和/或气态烃的开采在油藏的使用期内是如何变化的。

可以使用本领域的已知技术构建给定油藏中的流体流动的数学模型。然而，这些技术不准确并且可以使用其它数值技术实施，从而达到估计的解，例如，流线技术和有限差分技术。美国专利第7,164,990中更特别地描述了这些技术及其缺点，通过引用并入本申请。

‘990专利一般涉及使用流线模拟器和有限差分模拟器确定涉及油气藏的流体流动并处理流体流动数据的方法。‘990专利描述了确定包含两种或更多流体组分的一组流体的流动，其包括确定一个或更多来自压力场的流线和求解沿着每一条流线的流体成分。该方法还包括求解沿着每一条流线的压力并建议了使用有限差分技术以求解沿着每一条流线的流体成分和/或压力。‘990专利包括例示了传统流线技术的流程图(图2)并计算了来自压力场的流线。然而，‘990专利没有教导或建议使用有限差分油藏模拟结果计算流线的系统或方法。此外，‘990专利没有教授或建议包括使用流线轨迹的流线技术，以分析有限差分流体模拟结果。

因此，需要使用流线计算而不需要传统流线模拟执行的流线技术。还需要包括使用流线轨迹以分析有限差分流动模拟结果的流线技术。最后，需要使用有限差分流动模拟结果以计算流线。

发明内容

因此，本发明通过提供使用基于单元的流动模拟结果以计算流线轨迹的系统和方法，满足了上述需要并克服了现有技术中的一个或更多的缺点。

在一个实施例中，本发明包括计算流线的方法，包括：(i)初始化模拟网格单元内的种子点，(ii)计算种子点的速度矢量，(iii)在速度矢量方向上选择增量距离，(iv)根据增量距离和速度矢量的方向计算新种子点的位置坐标；v)计算种子点的另一个速度矢量；vi)在另一个速度矢量的方向上选择增量距离；vii)根据增量距离和另一个速度矢量的方向计算另一个新种子点的位置坐标；viii)根据每一个各自的速度矢量和各自的增量距离的速度或每一个各自的另一个速度矢量和各自的增量距离的速度，计算每一个增量距离的增量飞行时间；ix)绘制每一个各自的增量距离的每一个增量飞行时间的颜色；x)根据每一个增量飞行时间的总和计算累计飞行时间；xi)在距注入井最近的流线的末端，将累计飞行时间调整至零；以及xii)显示每一种颜色，其代表距注入井最近的流线的末端和距开采井或滞点区域最近的流线的另一个末端之间的流线上的累计飞行时间。

在另一个实施例中，本发明包括具有计算机可执行指令以计算流线的程序载体装置。指令是可执行的，以实现(i)初始化模拟网格单元内的种子点；(ii)计算种子点的速度矢量；(iii))在速度矢量方向上选择增量距离；(iv)根据增量距离和速度矢量的方向计算新种子点的位置坐标；v)计算种子点的另一个速度矢量；vi)在另一个速度矢量的方向上选择增量距离；vii)根据增量距离和另一个速度矢量的方向计算另一个新种子点的位置坐标；viii)根据每一个各自的速度矢量和各自的增量距离的速度或每一个各自的另一个速度矢量和各自的增量距离的速度，计算每一个增量距离的增量飞行时间；ix)绘制每一个各自的增量距离的每一个增量飞行时间的颜色；x)根据每一个增量飞行时间的总和计算累计飞行时间；xi)在距注入井最近的流线的末端，将累计飞行时间调整至零；以及xii)显示每一种颜色，其代表距注入井最近的流线的末端和距开采井或滞点区域最近的流线的另一个末端之间的流线上的累计飞行时间。

从下面的各种实施方式和相关附图的描述，本发明的其它方面、优点和实施方式对本领域的技术人员将会变得显而易见。

附图说明

下面将参考附图描述本发明，其中，相同的部件使用相同的参考标号，其中：

图1A是例示了并入本发明的工作流程的流程图。

图1B是例示了实施本发明的方法的流程图。

图2A例示了根据本发明计算的流线轨迹的示例性显示图(display)。

图2B例示了根据本发明的流线轨迹的示例性计算。

图3A例示了本发明的流线分析技术产生的流动场的示例性显示图。

图3B例示了根据相对于注入端的飞行时间过滤流线追踪的分析技术的示例性显示。

图3C例示了根据相对于开采端的飞行时间过滤流线追踪的分析技术的示例性显示。

图4例示了本发明的流线分析技术产生的井分布因素的示例性显示图。

图5例示了根据本发明计算的流线示例性显示。

图6是图5中例示的流线的部分二维显示图，其具有飞行时间信息。

图7是图5中例示的流线的三维显示图，其包括网格块。

图8是例示了实施本发明的计算机系统的方框图。

具体实施方式

下面具体地描述了本发明的主题，但描述本身并不限制本发明的范围。因而，结合其它现有或未来的技术，本发明的主题也可以以其它方式实施，以包括不同的步骤或类似于在此描述的步骤的结合。此外，虽然在此可以使用术语“步骤”以描述采用的方法的不同组成部分，但不应该将术语解释为暗示在此公开的各种步骤之间的任何特定的顺序，除非通过特定顺序的描述明确限制。

方法描述

现在参考图1A，例示了实施本发明的工作流程的流程图。在步骤101中，通过用户界面将标准有限差分流动模拟数据输入计算机系统。

在步骤102中，输入数据进入流动模拟的流动模拟器。可接受的流动模拟器可以包括，例如，兰德马克图形公司(Landmark GraphicsCorporation)销售的商业软件应用

和

流动模拟器通过由三维模拟(网格)单元的六个侧面(表面)计算相位通量而产生通量数据，其在步骤103作为每一个单元表面的模拟通量写入计算机系统中的数据库，例如，

结果数据库(“VDB”)。例如，归一化的速度被输出到数据库并可以被定义为：

速度×表面面积

总体积(Bulk Volume)

如下面的表1所描述的，对于每一相，绘制了六个速度数组(arrays)，每一个坐标有两个数组。

相	输出速度数组	描述
			油	VLNOX-VLNOX+VLNOY-VLNOY+VLNOZ-VLNOZ+	网格单元的表面的每一个方向的归一化油相速度：一个在+方向，一个在-方向。
气	VLNGX-VLNGX+VLNGY-VLNGY+VLNGZ-VLNGZ+	网格单元的表面的每一个方向的归一化气相速度：一个在+方向，一个在-方向。
			水	VLNWX-VLNWX+VLNWY-VLNWY+VLNWZ-VLNWZ+	网格单元的表面的每一个方向的归一化水相速度：一个在+方向，一个在-方向。

表1

在步骤104中，流线分析模块(“StreamCalc”)从数据库读取通量数据并计算流线轨迹和飞行时间(time-of-flight)推出的量。StreamCalc结果也可以包括示意图、平面线图和流线数据，其可以被写入用于读取通量数据的相同的数据库。数据库可以包括包含需要在步骤106中显示流线轨迹的流线和其它输出文件夹和文件的组。

图2A例示了在步骤104计算的流线轨迹的示例性显示图。沿着流线的色差是当粒子从流线的一端穿到另一端的飞行时间的量化。颜色编码流线替代色差处理流线可以优选为量化粒子的飞行时间。无论怎样，在此描述的颜色编码流线可以包括绘制相同的颜色作为不同的色调(即，色差)或绘制不同的颜色。通过在模拟网格单元200的中心开始，并且沿着在该位置的流动方向的速度矢量204a，增量地移动坐标为(x，y，z)的种子点202，从而计算流线206。速度矢量204a在新种子点位置被更新，并且重复该过程，直至新的种子点位置坐标到达产油井或零速度点(滞点)。该过程可以通过沿着在位置流动的反向的另一个速度矢量204b增量地移动坐标为(x，y，z)的种子点202而在相反的方向上重复。另一个速度矢量204b可以在新的种子点位置被更新，并且过程重复直至新的种子点位置坐标达到注入井或零速度点(滞点)。该过程可以在图2A中的每一个单元中重复，这使得将单元体积标记至注入和产油井或至滞点区域。对于附加的分辨率，可以在单元内使用多个种子点(开始点)，在这种情况下，与种子点的数量成比例地给每一条流线分配单元体积。

因此，本发明允许使用独特的基于模拟器的流线的流线分析技术，以分析有限差分油藏模拟结果。可以使用StreamCalc结果分析非对流条件，例如，初次采油、高度可压缩流体、重力控制的流动和对毛细管效应敏感的问题。可以使用StreamCalc结果确定波及系数，可视化成对的注入孔-产油井分组以及使用飞行时间方法理解其初次回采机制(primary recovery mechanisms)。从StreamCalc结果，模拟和非模拟用户可以提出建议，从而为有效的体积波及而修改水注入速率并最大化石油开采。结果，可以使用StreamCalc结果执行流线分析，而不需要传统的流线模拟。

使用StreamCalc结果的各种示例性流线分析技术包括，但不限于：1)使用点轨迹以显示流动场(图3A)；2)使用沿着轨迹的行程时间，也称作飞行时间(“TOF”)，作为单元的流体达到产油井或相反地作为注入的流体到达特定的单元所需时间的指示；3)识别由特定井排出或影响的油藏区域；4)计算井分布因素，即，穿过井的开采或注入百分比(图4)；5)计算油藏波及系数，即，油藏体积联系为采出或注入的流体体积的百分比；以及6)流体突破时间，表示为具有TOF最小值的流线，其连接成对的注入孔-产油井。在图3A中，流线300根据普通的成对的注入孔-产油井被着色，较浅的色差代表井1和井2之间的连接，而较深的色差代表井1和井3之间的连接。同样，颜色编码可以优选于色差处理流线300。在图4中，不同的色差代表补偿井对特定开采或注入井的百分比贡献。

在步骤105，以列表(分析表)的形式将步骤104的结果存在数据库和/或其它ASC II文件中。

可以将示例性输出文件以列出的形式描述在下面的表2中。

表2

在步骤106中，可以使用用于分析的兰德马克图形公司销售的另一种商业软件应用

显示流线轨迹和其它相关数据。

通过经流线显示和飞行时间显示量化并可视化扫描(sweeps)，提供了可以估算流动速率的工具。在图5中，例如，简单的流线500的示例性

显示图例示了各种注入井和产油井的连接。如图6中的示例性

显示图所例示的，也可以将二维飞行时间信息绘制为流线600。通过在任一流线中置入粒子，作为时间的函数，可以将沿着流线的距离绘制为速度场的流线。通过色差或颜色编码该信息，可以沿着流线产生飞行时间信息。如图3B和3C中例示的示例性

显示图，也可以使用该信息过滤飞行时间信息，并提供前缘的可视化表示，源自各个注入孔并在产油井结束。在图3B中，流线300在相对于流线注入端TOF＝7000天处被截，因而例示了相对于注入井1的前缘。绿色流线连接井1和井2，而蓝色流线连接井1和井3。在图3C中，流线300在相对于流线开采端TOF＝7000天处被截，因而例示了相对于产油井2和产油井3的前缘。绿色流线连接井1和井2，而蓝色流线连接井1和井3。

例如，可以产生如图5说明的相同类型信息的三维视图，并且如图7中示例性

显示图说明的，可以由流线702或网格块704编码。网格块704上的色差代表粒子随着时间的推移而流向的产油井，不同的色差代表属于不同井的不同排油区。作为替代，可以将色差或颜色编码绘制为流线，从而代表更传统的成对的注入井-产油井，如图3A-3C中的示例性

显示图所说明的。

现在参考图1B和2B，例示了用于计算图1A中的步骤104的流线轨迹的优选方法和示例性计算。在步骤104a中，将模拟网格单元200的中心的种子点202初始化(Xo，Yo，Zo)，并且将飞行时间设置为零。

在步骤104b中，计算种子点202的速度矢量204a。可以使用由Liu，Parker，及Camilleri发表在Society of Petroleum Engineers的文章“A New Particle Tracking Algorithm For Tracer Flow Simulation”(“SPE 51905”)(将其以引用的方式并入本申请)中描述的技术和算法计算速度矢量204a。使用SPE 51905中描述的方法和算法，可以计算在坐标空间的任何点(x，y，z)的速度矢量204a。

在步骤104c中，使用现有技术中的已知技术选择在速度矢量204a方向上的预定增量距离(DL)并计算种子点202的新坐标(Xi，Yi，Zi)。

在步骤104d中，使用速度矢量204a的增量距离(DL)和计算速度(VEL)计算增量飞行时间(DT)：DT＝DL/VEL。

在步骤104e中，计算累计的飞行时间：TOF＝TOF+DT。

在步骤104f中，使用新位置的坐标(Xi，Yi，Zi)重新初始化种子点202。

在步骤104g中，该方法可以从步骤104b开始重复。

如图2B所说明的，每一次流线穿入相邻的网格单元，可以使用SPE 51905中描述的技术计算小于DL的D1，以及相应的DL，使得将种子点202的新位置(Xi，Yi，Zi)计算为与网格单元边界恰好一致。然后，使用来自相邻的网格单元的速度使流线前进。

当在步骤104b中计算的速度矢量小于标示了滞点区域的某预定值，或当种子点202的位置移近产油井(阱)时，图1B中说明的方法可以终止。相反地，该方法可以从种子点202开始重复，但速度矢量乘以-1，以代表相反的流动方向。这时，当在步骤104b中计算的速度矢量小于标示了滞点区域的预定值，或当种子点202的位置移近注入井(源)时，可以将方法终止。因而描述的方法的重复的飞行时间结果可以被结合并调整，使得飞行时间在流线的源端呈现为零，并在朝阱端的方向增加。

系统描述

本发明可以通过计算机可执行程序指令实现，诸如程序模块，一般称作由计算机执行的软件应用或应用程序。软件可以包括，例如，执行特定任务或实现特定抽象数据类型的例行程序、程序、对象、组件和数据结构。软件形成接口使计算机根据输入源而反应。兰德马克图形公司销售的商业软件

可以作为接口应用，以实现本发明。软件结合接收数据的源，也可以与其它代码段合作，以开始响应于所接收的数据的多种任务。可以将软件存储和/或携载于任一种存储媒体上，诸如CD-ROM、磁盘、磁泡存储器和半导体存储器(例如，各种类型的RAM或ROM)上。而且，可以通过诸如光纤、金属线、自由空间的多种载体媒体和/或诸如互联网的多种网络中的任一种传输软件及其结果。

此外，本领域的技术人员将认识到本发明可以用多种计算机-系统配置实施，包括手持装置、多处理器系统、基于微处理器的或可编程消费类电子产品、微型计算机、大型计算机及类似。任意数量的计算机-系统和计算机网络都可以使用本发明。可以在分布式计算环境中实施本发明，在分布式计算环境中，任务可以由通过通信网络链接的远程处理装置执行。在分布式计算环境中，程序模块可以位于包括记忆存储装置的本地和远程计算机-存储媒体中。因此，可以连同各种硬件、软件或其结合在计算机系统或其它处理系统中实施本发明。

现在参考图8，说明了在计算机上实施本发明的系统的方框图。该系统包括计算单元，有时称作计算系统，包括存储器、应用程序、数据库、ASC II文件、用户接口和处理单元。计算单元只是一个适合的计算环境的示例并不限制本发明的用途或功能的范围。

存储器主要存储由计算单元执行以实现在此描述的工作流程和方法的应用程序，也可描述为包括计算机可执行指令的程序模块。因此，存储器包括实现图1B中描述的方法的StreamCalc模块和其它部分。

虽然所示计算单元具有广义存储器，该计算单元一般包括多种计算机可读媒体。通过举例的方式但不限制，计算机可读媒体可以包括计算机存储媒体和通信媒体。计算系统存储器可以包括诸如只读存储器(ROM)和随机存储器(RAM)的形式为易失性存储器和/或非易失性存储器的计算机存储媒体。基本输入/输出系统(BIOS)一般存储在ROM中，BIOS包括帮助在计算单元内的组成部分(element)之间传送信息的基本例程，诸如在启动过程中。RAM一般包括数据和/或程序模块，其可以由处理单元立即访问和/或正在操作。通过举例的方式但不限制，计算单元包括操作系统、应用程序、其它程序模块和程序数据。

存储器中示出的组件也可以包括在其它可擦除/不可擦除、易失性/非易失性计算机存储媒体中。仅例如，硬盘驱动器可以从不可擦除、非易失磁性媒体读取，或写入不可擦除、非易失磁性媒体，磁盘驱动器可以从可擦除、非易失磁盘中读取，或写入可擦除、非易失磁盘，以及光盘驱动可以从诸如CD ROM或其它光学媒体的可擦除、非易失光盘读取，或写入可擦除、非易失光盘。可以在示例性操作环境中使用的其它可擦除/不可擦除，易失性/非易失性计算机存储媒体可以包括但不限于盒式磁带、闪存卡、数字多用盘、数字视频磁带、固态RAM、固态ROM及类似。因此，上面讨论的驱动和其相关的计算机存储媒体存储和/或载有计算单元的计算机可读指令、数据结构、程序模块和其它数据。

用户可以通过用户接口将命令和信息输入计算单元，该用户接口可以是输入装置，诸如键盘和指点装置，其通常称作鼠标、轨迹球或触控板。输入装置可以包括麦克风、控制杆、卫星天线、扫描仪或类似。

通常这些和其它输入装置通过用户接口连接至处理单元，用户接口连接至系统总线，但是这些和其它输入装置也可以通过诸如并行端口或通用串行总线(USB)的其它接口和总线结构连接。监视器或其它类型的显示装置可以通过诸如视频接口的接口连接至系统。除了监视器以外，计算机也可以包括其它外围输出装置，诸如扬声器和打印机，其可以通过输出外围接口连接。

虽然没有示出计算单元的很多其它内部组件，本领域的普通技术人员应认识到这些组件和互连是众所周知的。

虽然已经结合优选实施方式描述了本发明，本领域的技术人员将理解本发明不限于这些实施方式。虽然本发明所说明的实施方式涉及在石油和天然气工业中使用速度通量数据，可以将本发明应用于其它领域和学科的任何适合的速度通量数据。因此，在不脱离所附权利要求及其等同的精神和范围的情况下，可以对公开的实施方式作出各种替代实施方式和变型。

Claims (14)

1.一种计算机实施的计算流线的方法，包括：

初始化模拟网格单元内的种子点；

计算所述种子点的速度矢量；

在所述速度矢量的方向上选择增量距离；以及

根据所述增量距离和所述速度矢量的方向计算新种子点的位置坐标；

计算所述种子点的另一个速度矢量；

在所述另一个速度矢量的方向上选择增量距离；

根据所述增量距离和所述另一个速度矢量的方向计算另一个新种子点的位置坐标；

根据每一个各自的速度矢量和各自的增量距离的速度或每一个各自的另一个速度矢量和各自的增量距离的速度，计算每一个增量距离的增量飞行时间；

绘制对每一个各自的增量距离的每一个增量飞行时间的颜色；

根据每一个增量飞行时间的总和计算累计飞行时间；

在距注入井最近的流线末端，将所述累计飞行时间调整至零；以及

显示每一种代表距所述注入井最近的流线末端和距开采井或滞点区域最近的流线另一个末端之间的所述流线上的累计飞行时间的颜色。

2.如权利要求1所述的方法，还包括：

a)使用所述新种子点的位置坐标重新初始化所述种子点；

b)使用所述另一个新种子点的位置坐标重新初始化所述种子点；

c)重复以下步骤：计算使用所述新种子点的位置坐标重新初始化的种子点的速度矢量，在所述速度矢量的方向上选择增量距离，并且根据所述增量距离和所述速度矢量的方向计算新种子点的位置坐标；

d)重复以下步骤：计算使用所述另一个新种子点的位置坐标重新初始化的种子点的另一个速度矢量，在所述另一个速度矢量的方向上选择增量距离，并且根据所述增量距离和所述另一个速度矢量的方向计算另一个新种子点的位置坐标；以及

e）重复步骤a-d，直至所述新种子点的位置坐标和所述另一个新种子点的位置坐标到达开采井、注入井或滞点区域。

3.如权利要求1所述的方法，还包括：对于每一个模拟网格单元重复权利要求1的所述步骤。

4.如权利要求1所述的方法，其中，所述种子点的速度矢量由所述网格单元的每一个单元表面的模拟通量计算。

5.如权利要求4所述的方法，其中，所述模拟通量由每一个单元表面的归一化油相速度、归一化气相速度和归一化水相速度表示。

6.如权利要求1所述的方法，其中，如果其包含在所述模拟网格单元中，则所选择的增量距离是预定的，并且如果所述流线穿入另一个模拟网格单元，则计算所选择的增量距离，使得所述新种子点的位置坐标与所述网格单元和所述另一个网格单元之间的边界一致，并且具有与所述边界一致的所述新种子点的位置坐标的所述种子点的速度矢量，由所述另一个网格单元的每一个单元表面的模拟通量计算。

7.如权利要求1所述的方法，还包括：

初始化所述模拟网格单元内的另一个种子点；

计算所述另一个种子点的速度矢量；

在所述速度矢量的方向上选择增量距离；以及

根据所述增量距离和所述速度矢量的方向计算新的另一个种子点的位置坐标。

8.一种计算机实施的计算流线的系统，其包括：

初始化模拟网格单元内的种子点的模块；

计算所述种子点的速度矢量的模块；

在所述速度矢量的方向上选择增量距离的模块；

根据所述增量距离和所述速度矢量的方向计算新的种子点的位置坐标的模块；

计算所述种子点的另一个速度矢量的模块；

在所述另一个速度矢量的方向上选择增量距离的模块；

根据所述增量距离和所述另一个速度矢量的方向计算另一个新种子点的位置坐标的模块；

根据每一个各自的速度矢量和各自的增量距离的速度或每一个各自的另一个速度矢量和各自的增量距离的速度来计算每一个增量距离的增量飞行时间的模块；

绘制对每一个各自的增量距离的每一个增量飞行时间的颜色的模块；

根据每一个增量飞行时间的总和计算累计飞行时间的模块；

在距注入井最近的流线末端将所述累计飞行时间调整至零的模块；以及

显示每一种代表距所述注入井最近的流线末端和距开采井或滞点区域最近的流线另一个末端之间的所述流线上的累计飞行时间的颜色的模块。

9.如权利要求8所述的系统，其还包括：

使用所述新种子点的位置坐标重新初始化所述种子点的模块a)；

使用所述另一个新种子点的位置坐标重新初始化所述种子点的模块b)；

重复计算使用所述新种子点的位置坐标重新初始化的种子点的速度矢量、在所述速度矢量的方向上选择增量距离并且根据所述增量距离和所述速度矢量的方向计算新种子点的位置坐标的步骤的模块c)；

重复计算使用所述另一个新种子点的位置坐标重新初始化的种子点的另一个速度矢量、在所述另一个速度矢量的方向上选择增量距离并且根据所述增量距离和所述另一个速度矢量的方向计算另一个新种子点的位置坐标的步骤的模块d)；以及

重复上述模块a)-d)的步骤直至所述新种子点的位置坐标和所述另一个新种子点的位置坐标到达开采井、注入井或滞点区域的模块e)。

10.如权利要求8所述的系统，其中，所述计算单元还包括：对每一个模拟网格重复权利要求8中各模块执行的步骤的模块。

11.如权利要求8所述的系统，其中，所述种子点的速度矢量由所述网格单元的每一个单元表面的模拟通量计算。

12.如权利要求11所述的系统，其中，所述模拟通量由每一个单元表面的归一化油相速度、归一化气相速度和归一化水相速度表示。

13.如权利要求8所述的系统，其中，如果其包含在所述模拟网格单元中，则所选择的增量距离是预定的，并且如果所述流线穿入另一个模拟网格单元，则计算所选择的增量距离，使得所述新种子点的位置坐标与所述网格单元和所述另一个网格单元之间的边界一致，并且具有与所述边界一致的所述新种子点的位置坐标的所述种子点的速度矢量，由所述另一个网格单元的每一个单元表面的模拟通量计算。

14.如权利要求8所述的系统，还包括：

初始化所述模拟网格单元内的另一个种子点的模块；

计算所述另一个种子点的速度矢量的模块；

在所述速度矢量的方向上选择增量距离的模块；以及

根据所述增量距离和所述速度矢量的方向计算新的另一个种子点的位置坐标的模块。

Images