CN107895076B - 一种基于目标的非贯穿型河道建模方法及系统 - Google Patents

一种基于目标的非贯穿型河道建模方法及系统 Download PDF

Info

Publication number
CN107895076B
CN107895076B CN201711105289.3A CN201711105289A CN107895076B CN 107895076 B CN107895076 B CN 107895076B CN 201711105289 A CN201711105289 A CN 201711105289A CN 107895076 B CN107895076 B CN 107895076B
Authority
CN
China
Prior art keywords
river channel
model
river
dimensional
boundary
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
CN201711105289.3A
Other languages
English (en)
Other versions
CN107895076A (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.)
Yangtze University
Original Assignee
Yangtze 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 Yangtze University filed Critical Yangtze University
Priority to CN201711105289.3A priority Critical patent/CN107895076B/zh
Publication of CN107895076A publication Critical patent/CN107895076A/zh
Application granted granted Critical
Publication of CN107895076B publication Critical patent/CN107895076B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation

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)
  • Processing Or Creating Images (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于目标的非贯穿型河道建模方法及系统,其中,生成河道中线,建立河道模型,构建不规则的边界线,通过样条函数定量描述边界线,将截断河道模型、河道边界线赋值到网格模型中实现河道的截断,无论是多套地层或一套地层内存在多组河道的情况下,做到一组河道对应一条边界,多组河道分别被多条边界影响,模拟结果比较好地刻画了河道与边界的关系,实现了多组河道在研究区内部终止,实现复杂边界条件的河道建模。建立符合实际需要的地质模型,提高油气采收率。

Description

一种基于目标的非贯穿型河道建模方法及系统
技术领域
本发明涉及油藏描述技术领域,尤其是涉及一种基于目标的非贯穿型河道建模方法及系统。
背景技术
储层地质建模作为现代油藏描述研究的核心内容,国内外学者开展了大量的研究,其中C.V.Deutsch设计的基于目标的建模方法(Fluvsim)能够较好地表征不同储层结构单元的成因关系,是一种重要的河流相储层建模方法。国内在利用该算法描述复杂河道形态及预测河道储层分布方面有不少成功案例。在Fluvsim算法中,河道的模拟是沿主流线方向一直贯穿整个研究区,无法人为控制河道的终止点。在某些情况下,河道在研究区内部终止,这时候传统的基于目标的河道建模方法无法刻画这种非贯穿整个研究区的河道。对于只有一组河道的情况,还可以采用网格分区的方式解决,即:根据这组河道的边界对网格进行分区,然后在模拟河道的时候限制在某个区块内部。当存在多套地层或一套地层内存在多组河道的情况下,不同组河道的方向可能不同,且不同组河道的边界也不一样,这时候利用网格分区的方式就无法解决多组非贯穿型河道的建模。鉴于现有建模技术存在的上述问题,亟待研究一种新的河道建模方法,以解决当需在研究区内建立复杂边界条件来刻画非贯穿型河道的问题。
发明内容
本发明的目的在于克服上述技术不足,提出一种基于目标的非贯穿型河道建模方法及系统,解决现有技术中的上述技术问题。
为达到上述技术目的,本发明的技术方案提供一种基于目标的非贯穿型河道建模方法,包括:
S1、生成若干河道中线节点,连接若干河道中线节点形成河道中线,河道中线贯穿预设三维研究区域;
S2、以河道中线为中心在三维研究区域内建立三维河道模型;
S3、在三维研究区域内确定河道边界线,并获取河道边界线对应的曲线函数,河道边界线在三维研究区域中截断所述三维河道模型;
S4、三维研究区域网格化后为第一模型,河道边界线在第一模型中穿过的网格为边界网格,识别出截断河道模型,截断河道模型为三维河道模型中从河道源头到边界网格的部分;建立网格化的第二模型,第二模型的大小和第一模型相同,将第一模型中截断河道模型、河道边界线映射到第二模型中,截断河道模型、河道边界线在第一模型、第二模型中的位置相同;
S5、循环执行步骤S1-S4,生成不同的三维河道模型和确定该三维河道模型的河道边界线,在第二模型中模拟出各截断河道模型及对应的河道边界线。
本发明还提供一种基于目标的非贯穿型河道建模系统,包括:
河道中线生成模块:生成若干河道中线节点,连接若干河道中线节点形成河道中线,河道中线贯穿预设三维研究区域;
河道模型生成模块:以河道中线为中心在三维研究区域内建立三维河道模型;
河道边界生成模块:在三维研究区域内确定河道边界线,并获取河道边界线对应的曲线函数,河道边界线在三维研究区域中截断所述三维河道模型;
河道截断模块:三维研究区域网格化后为第一模型,河道边界线在第一模型中穿过的网格为边界网格,识别出截断河道模型,截断河道模型为三维河道模型中从河道源头到边界网格的部分;建立网格化的第二模型,第二模型的大小和第一模型相同,将第一模型中截断河道模型、河道边界线映射到第二模型中,截断河道模型、河道边界线在第一模型、第二模型中的位置相同;
循环模块:循环执行河道中线生成模块至河道截断模块,生成不同的三维河道模型和确定该三维河道模型的河道边界线,在第二模型中模拟出各截断河道模型及对应的河道边界线。
与现有技术相比,本发明的有益效果包括:本发明通过建立河道模型,构建不规则的边界线,通过样条函数定量描述边界线,将截断河道模型、河道边界线赋值到网格模型中实现河道的截断,无论是多套地层或一套地层内存在多组河道的情况下,做到一组河道对应一条边界,多组河道分别被多条边界影响,模拟结果比较好地刻画了河道与边界的关系,实现了多组河道在研究区内部终止,实现复杂边界条件的河道建模。建立符合实际需要的地质模型,提高油气采收率。
附图说明
图1是本发明提供的一种基于目标的非贯穿型河道建模方法流程图;
图2是本发明提供的一种基于目标的非贯穿型河道建模系统结构框图;
图3是河道方向与河道中线示意图;
图4是多组河道边界示意图。
附图中:1、基于目标的非贯穿型河道建模系统,11、河道中线生成模块,12、河道模型生成模块,13、河道边界生成模块,14、河道截断模块,15、循环模块。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
如图1,本发明提供了一种基于目标的非贯穿型河道建模方法,包括:
S1、生成若干河道中线节点,连接若干河道中线节点形成河道中线,河道中线贯穿预设三维研究区域;
S2、以河道中线为中心在三维研究区域内建立三维河道模型;
S3、在三维研究区域内确定河道边界线,并获取河道边界线对应的曲线函数,河道边界线在三维研究区域中截断所述三维河道模型;
S4、三维研究区域网格化后为第一模型,河道边界线在第一模型中穿过的网格为边界网格,识别出截断河道模型,截断河道模型为三维河道模型中从河道源头到边界网格的部分;建立网格化的第二模型,第二模型的大小和第一模型相同,将第一模型中截断河道模型、河道边界线映射到第二模型中,截断河道模型、河道边界线在第一模型、第二模型中的位置相同;
S5、循环执行步骤S1-S4,生成不同的三维河道模型和确定该三维河道模型的河道边界线,在第二模型中模拟出各截断河道模型及对应的河道边界线。
本发明所述的基于目标的非贯穿型河道建模方法,步骤S1中:
如图3所示,A为河道主方向,D为河道中线节点,B为河道中线,C为河道中线节点相对于河道主方向的偏移量;
若干河道中线节点相对于河道主方向的各偏移量呈高斯分布,河道主方向根据预设河道方位角确定,河道主方向可为一条直线;河道中线节点相对于河道主方向的偏移量为河道中线节点到河道主方向直线的垂直距离。
本发明所述的基于目标的非贯穿型河道建模方法,步骤S2中:
三维河道模型中各河道中线节点处的河道宽度呈高斯分布,三维河道模型中各河道中线节点处的河道厚度呈高斯分布。
本发明所述的基于目标的非贯穿型河道建模方法,步骤S3中:
确定河道边界线后,如果河道边界线为不规则曲线,则将河道边界线离散成若干密集的河道边界点,利用三次样条插值的方法获取连接若干河道边界点的样条曲线的函数即为河道边界线对应的曲线函数;
具体的,河道边界线的形态可以为规则的,也可以为不规则的,因此定义边界线也可以有很多种方式,为了方法的适用性,满足不同情况下建模的需要,本发明确定边界线可采用手工绘图的方式确定,地质工作者根据已有的地质资料和建模时的要求,确定河道边界线的位置。手工绘图的河道边界线往往为不规则边界线,难以用现成的公式去刻画,也就是难以定量化描述整个边界线,这个时候将不规则边界线离散成若干密集的河道边界点,离散出的河道边界点越密集反映的河道边界线越准确,利用样条曲线函数表示河道边界线,从而实现河道边界线定量化描述。河道边界线在第一模型中穿过的网格为边界网格,如果河道边界线没有定量化,无法准确知晓哪些是边界网格,河道边界线定量化后,则可以计算出边界网格。
本发明所述的基于目标的非贯穿型河道建模方法,步骤S4中识别截断河道模型的步骤为:
遍历第一模型中三维河道模型内的每一个网格,当一网格和河道源头所处的网格位于边界网格的同一侧时,该网格为截断河道网格,所有截断河道网格组成截断河道模型。
上述技术方案中,预设三维研究区域大小为500×500×30米,第一模型和第二模型的每个网格大小均为1.25×1.25×0.05米。
采用本发明的技术方案后,如图4中,两条河道分别被两条河道边界线截断,从图中可以看出河道1-2被边界1-1正确截断,而不受边界2-1的影响,河道2-2被边界2-1正确截断,而不受边界1-1影响。做到了多组河道分别被多组边界影响,模拟结果比较好地刻画了河道与边界的关系,实现了河道在研究区内部终止。
如图2,本发明还提供一种基于目标的非贯穿型河道建模系统1,包括:
河道中线生成模块11:生成若干河道中线节点,连接若干河道中线节点形成河道中线,河道中线贯穿预设三维研究区域;
河道模型生成模块12:以河道中线为中心在三维研究区域内建立三维河道模型;
河道边界生成模块13:在三维研究区域内确定河道边界线,并获取河道边界线对应的曲线函数,河道边界线在三维研究区域中截断所述三维河道模型;
河道截断模块14:三维研究区域网格化后为第一模型,河道边界线在第一模型中穿过的网格为边界网格,识别出截断河道模型,截断河道模型为三维河道模型中从河道源头到边界网格的部分;建立网格化的第二模型,第二模型的大小和第一模型相同,将第一模型中截断河道模型、河道边界线映射到第二模型中,截断河道模型、河道边界线在第一模型、第二模型中的位置相同;
循环模块15:循环执行河道中线生成模块至河道截断模块,生成不同的三维河道模型和确定该三维河道模型的河道边界线,在第二模型中模拟出各截断河道模型及对应的河道边界线。
本发明所述的基于目标的非贯穿型河道建模系统1,河道中线生成模块11中:
若干河道中线节点相对于河道主方向的各偏移量呈高斯分布,河道主方向根据预设河道方位角确定。
本发明所述的基于目标的非贯穿型河道建模系统1,河道模型生成模块12中:
三维河道模型中各河道中线节点处的河道宽度呈高斯分布,三维河道模型中各河道中线节点处的河道厚度呈高斯分布。
本发明所述的基于目标的非贯穿型河道建模系统1,河道边界生成模块13中:
确定河道边界线后,如果河道边界线为不规则曲线,则将河道边界线离散成若干密集的河道边界点,利用三次样条插值的方法获取连接若干河道边界点的样条曲线的函数即为河道边界线对应的曲线函数。
本发明所述的基于目标的非贯穿型河道建模系统1,河道截断模块14中:
遍历第一模型中三维河道模型内的每一个网格,当一网格和河道源头所处的网格位于边界网格的同一侧时,该网格为截断河道网格,所有截断河道网格组成截断河道模型。
与现有技术相比,本发明的有益效果包括:本发明通过建立河道模型,构建不规则的边界线,通过样条函数定量描述边界线,将截断河道模型、河道边界线赋值到网格模型中实现河道的截断,无论是多套地层或一套地层内存在多组河道的情况下,做到一组河道对应一条边界,多组河道分别被多条边界影响,模拟结果比较好地刻画了河道与边界的关系,实现了多组河道在研究区内部终止,实现复杂边界条件的河道建模。建立符合实际需要的地质模型,提高油气采收率。
以上所述本发明的具体实施方式,并不构成对本发明保护范围的限定。任何根据本发明的技术构思所做出的各种其他相应的改变与变形,均应包含在本发明权利要求的保护范围内。

Claims (8)

1.一种基于目标的非贯穿型河道建模方法,其特征在于,包括:
S1、生成若干河道中线节点,连接若干河道中线节点形成河道中线,所述河道中线贯穿预设三维研究区域;
S2、以所述河道中线为中心在所述三维研究区域内建立三维河道模型;
S3、确定河道边界线后,如果所述河道边界线为不规则曲线,则将所述河道边界线离散成若干密集的河道边界点,利用三次样条插值的方法获取连接若干河道边界点的样条曲线的函数即为河道边界线对应的曲线函数,河道边界线在三维研究区域中截断所述三维河道模型;
S4、三维研究区域网格化后为第一模型,所述河道边界线在所述第一模型中穿过的网格为边界网格,识别出截断河道模型,所述截断河道模型为三维河道模型中从河道源头到边界网格的部分;建立网格化的第二模型,所述第二模型的大小和第一模型相同,将第一模型中截断河道模型、河道边界线映射到所述第二模型中,截断河道模型、河道边界线在所述第一模型、所述第二模型中的位置相同;
S5、循环执行步骤S1-S4,生成不同的三维河道模型和确定该三维河道模型的河道边界线,在所述第二模型中模拟出各截断河道模型及对应的河道边界线。
2.如权利要求1所述的基于目标的非贯穿型河道建模方法,其特征在于,步骤S1中:
若干河道中线节点相对于河道主方向的各偏移量呈高斯分布,所述河道主方向根据预设河道方位角确定。
3.如权利要求1所述的基于目标的非贯穿型河道建模方法,其特征在于,步骤S2中:
三维河道模型中各河道中线节点处的河道宽度呈高斯分布,所述三维河道模型中各河道中线节点处的河道厚度呈高斯分布。
4.如权利要求1所述的基于目标的非贯穿型河道建模方法,其特征在于,步骤S4中识别截断河道模型的步骤为:
遍历所述第一模型中三维河道模型内的每一个网格,当一网格和河道源头所处的网格位于边界网格的同一侧时,该网格为截断河道网格,所有所述截断河道网格组成截断河道模型。
5.一种基于目标的非贯穿型河道建模系统,其特征在于,包括:
河道中线生成模块:生成若干河道中线节点,连接若干河道中线节点形成河道中线,所述河道中线贯穿预设三维研究区域;
河道模型生成模块:以所述河道中线为中心在所述三维研究区域内建立三维河道模型;
河道边界生成模块:确定河道边界线后,如果所述河道边界线为不规则曲线,则将所述河道边界线离散成若干密集的河道边界点,利用三次样条插值的方法获取连接若干河道边界点的样条曲线的函数即为河道边界线对应的曲线函数,河道边界线在三维研究区域中截断所述三维河道模型;
河道截断模块:所述三维研究区域网格化后为第一模型,所述河道边界线在所述第一模型中穿过的网格为边界网格,识别出截断河道模型,所述截断河道模型为三维河道模型中从河道源头到边界网格的部分;建立网格化的第二模型,所述第二模型的大小和第一模型相同,将第一模型中截断河道模型、河道边界线映射到所述第二模型中,截断河道模型、河道边界线在所述第一模型、所述第二模型中的位置相同;
循环模块:循环执行河道中线生成模块至河道截断模块,生成不同的三维河道模型和确定该三维河道模型的河道边界线,在所述第二模型中模拟出各截断河道模型及对应的河道边界线。
6.如权利要求5所述的基于目标的非贯穿型河道建模系统,其特征在于,河道中线生成模块中:
若干河道中线节点相对于河道主方向的各偏移量呈高斯分布,所述河道主方向根据预设河道方位角确定。
7.如权利要求5所述的基于目标的非贯穿型河道建模系统,其特征在于,河道模型生成模块中:
三维河道模型中各河道中线节点处的河道宽度呈高斯分布,所述三维河道模型中各河道中线节点处的河道厚度呈高斯分布。
8.如权利要求5所述的基于目标的非贯穿型河道建模系统,其特征在于,河道截断模块中:
遍历所述第一模型中三维河道模型内的每一个网格,当一网格和河道源头所处的网格位于边界网格的同一侧时,该网格为截断河道网格,所有所述截断河道网格组成截断河道模型。
CN201711105289.3A 2017-11-10 2017-11-10 一种基于目标的非贯穿型河道建模方法及系统 Active CN107895076B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711105289.3A CN107895076B (zh) 2017-11-10 2017-11-10 一种基于目标的非贯穿型河道建模方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711105289.3A CN107895076B (zh) 2017-11-10 2017-11-10 一种基于目标的非贯穿型河道建模方法及系统

Publications (2)

Publication Number Publication Date
CN107895076A CN107895076A (zh) 2018-04-10
CN107895076B true CN107895076B (zh) 2021-04-27

Family

ID=61804911

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711105289.3A Active CN107895076B (zh) 2017-11-10 2017-11-10 一种基于目标的非贯穿型河道建模方法及系统

Country Status (1)

Country Link
CN (1) CN107895076B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108614950A (zh) * 2018-05-18 2018-10-02 中国电建集团华东勘测设计研究院有限公司 一种基于bim的河道类线性工程设计建模方法
CN111997604B (zh) * 2020-10-16 2023-12-08 陕西延长石油(集团)有限责任公司 一种单期河道边界的识别方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105022856A (zh) * 2015-04-17 2015-11-04 长江大学 预测高弯度曲流河储层内部结构的储层建模方法
CN105741045A (zh) * 2016-02-05 2016-07-06 中国水利水电科学研究院 一种洪水风险动态分析与展示系统及方法
CN105844709A (zh) * 2016-03-25 2016-08-10 中国水利水电科学研究院 复杂河道地形流域洪水演进虚拟仿真的淹没线追踪方法
CN106934860A (zh) * 2017-01-12 2017-07-07 天津大学 一种基于t样条的三维地质建模方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105022856A (zh) * 2015-04-17 2015-11-04 长江大学 预测高弯度曲流河储层内部结构的储层建模方法
CN105741045A (zh) * 2016-02-05 2016-07-06 中国水利水电科学研究院 一种洪水风险动态分析与展示系统及方法
CN105844709A (zh) * 2016-03-25 2016-08-10 中国水利水电科学研究院 复杂河道地形流域洪水演进虚拟仿真的淹没线追踪方法
CN106934860A (zh) * 2017-01-12 2017-07-07 天津大学 一种基于t样条的三维地质建模方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
FLUVSIM: a program for object-based stochastic modeling of fluvial depositional systems;C.V. Deutsch等;《Computers & Geosciences》;20021231;第525-535页 *
河流相建模算法Fluvsim的一种改进方法;陈迪等;《测绘与空间地理信息》;20141231;第37卷(第12期);第56-59页 *

Also Published As

Publication number Publication date
CN107895076A (zh) 2018-04-10

Similar Documents

Publication Publication Date Title
CN105334536B (zh) 致密砂岩储层网状裂缝系统有效性评价方法
CN105651676B (zh) 一种水平井规则开发井网下的储层非均质性表征方法
CN104616350A (zh) 缝洞型碳酸盐油藏三维物理模型建立方法
CN105626023A (zh) 低渗透油藏垂直压裂裂缝方位试井确定方法
CN105279790A (zh) 裂缝网络三维数字岩心建模方法
CN105184864A (zh) 一种用于天然地基换填工程量计算的场地地层三维地质结构模型生成方法
CN107895076B (zh) 一种基于目标的非贯穿型河道建模方法及系统
CN103728671B (zh) 确定碳酸盐岩潜山含油气性的方法和装置
CN105631078B (zh) 天然裂缝性油藏自适应介质的数值模拟方法
CN110955982B (zh) 一种变质岩储层渗透率计算方法、装置及计算机存储介质
CN104933208B (zh) 一种不同尺度裂缝油藏数值模拟综合处理方法
MX2010013221A (es) Distribucion de propiedades en un modelo volumetrico tridimensional utilizando un campo de continuidad maxima.
CN103942842A (zh) 嵌入式曲流河砂体建模方法
CN108053477A (zh) 一种管道内变形的数值处理方法
CN103941288A (zh) 一种基于平面导航的三维目标地质体解释方法
CN106227943B (zh) 一种二维不同表面粗糙度的随机孔隙生成方法
CN112394404A (zh) 一种渐进式储层精细表征方法
CN105095540A (zh) 一种井间相对位置关系的法面扫描解析方法
CN110362726B (zh) 一种城市地下排水管网数据补充完善方法及系统
CN105653746B (zh) 基于嵌入式离散裂缝模型的压裂井建模及模拟方法
CN106780740A (zh) 基于gocad软件的坝基岩体优势渗透通道识别方法
CN105068141A (zh) 确定断缝系统的岩溶储层分布的方法和装置
CN106652033B (zh) 一种地质剖面自然网格的剖分方法
CN108595782A (zh) 一种离散裂缝中基质与裂缝间的传质计算方法
CN108242069B (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