CN105204090A - 复杂形态砂体的布尔模拟方法 - Google Patents
复杂形态砂体的布尔模拟方法 Download PDFInfo
- Publication number
- CN105204090A CN105204090A CN201510651446.5A CN201510651446A CN105204090A CN 105204090 A CN105204090 A CN 105204090A CN 201510651446 A CN201510651446 A CN 201510651446A CN 105204090 A CN105204090 A CN 105204090A
- Authority
- CN
- China
- Prior art keywords
- shape
- template
- nodes
- river channel
- sandbody
- 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
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种复杂形态砂体的布尔模拟方法,该方法先建立空白模板,依据河道剖面的宽度、深度绘制河道砂体横剖面的形状,得到形状模板,形状模板所有节点构成二维矩阵K,形状模板的节点在河道砂体横剖面的形状内部,矩阵元素km,n=1;当形状模板的节点在河道砂体横剖面的形状外部,矩阵元素km,n=0;定义维度为W×H的工区网格G,布尔算法序贯访问工区网格G的每个节点,把形状模板中与其对应的所有节点的值赋给它,得到复杂几何形态砂体的剖面模型。本发明不受椭圆形状参数来确定河道砂体横剖面形状的约束,建立复杂几何形态的砂体图像并进行布尔模拟,使模拟结果更加真实。
Description
技术领域
本发明属于油藏描述技术领域,具体地指一种复杂形态砂体的布尔模拟方法。
背景技术
储层三维地质建模技术和油藏数值模拟技术对于油田高效开发具有重要的支持作用,随着油田对于储层研究精细化的要求越来越高,对于这两项技术也提出了更高的要求。布尔方法是储层建模方法中最简单的一种算法,布尔方法的原理如下所示:
设U为坐标随机变量,Xk是表征第k类物体几何特征(形状、大小、方向)的参数随机变量;第k类几何物体中心点的分布构成一点过程U,它可以用形状随机过程Xk和表示第k类几何物体出现与否的指标随机过程Ik两者的联合分布“示性”,从而构成一示性点过程。其中
布尔方法就是依据一定的概率定律,按照空间物体分布统计规律产生这些物体中心点的空间分布,并通过2×k个随机函数Xk(u),Ik(u,k)(k=1,2,3,…,k,u∈定义域)的联合分布,确定中心点在此处的物体的几何形状、大小、属性。
国内外学者开展了大量的研究,其中斯坦福大学开发的GSLIB软件库中布尔模拟算法Ellipsim能简单表征砂体叠加分布特征,是一种以河道砂体为重建目标的储层建模方法。Ellipsim算法模拟椭圆形态河道砂体结构,砂体的空间分布从给定的均匀分布函数中通过随机抽样得到的。这样做能够满足常见的河流相储层砂体结构的建模需求,但是这种Ellipsim算法无法刻画几何形态复杂的砂体形态,例如具透镜状形态的河道砂体,其几何形态(宽度与厚度比)、侧向连续型等均有明显差异。
因此,急需研究一种新的河道砂体重建方法,以解决刻画复杂几何形态砂体的建模问题。
发明内容
本发明的目的是克服现有技术的不足,提供一种复杂形态砂体的布尔模拟方法。
为实现上述目的,本发明采用以下技术方案:
复杂形态砂体的布尔模拟方法,包括以下步骤:
a、确定单元尺寸,所述单元尺寸包括沿水平方向和垂直方向的尺寸;
b、建立至少一个空白模板,所述空白模板由多个步骤a所述的单元组成,每个单元构成一个节点,依据河道剖面的宽度、深度分别确定所述空白模板沿水平方向和垂直方向的单元数量M以及N;
c、在步骤b所述的空白模板中绘制河道砂体横剖面的形状,得到形状模板;
d、步骤c所述的形状模板的所有节点构成维度为M×N的二维矩阵K,当形状模板的节点在河道砂体横剖面的形状内部,矩阵元素km,n=1;当形状模板的节点在河道砂体横剖面的形状外部,矩阵元素km,n=0;
e、定义维度为W×H的工区网格G,所述工区网格G的所有节点构成二维矩阵Q,所有矩阵元素Qw,h均为空值;
f、布尔算法序贯访问所述工区网格G的每个节点,在模拟节点G(w,h)时,将步骤d所述形状模板中与其对应的所有节点的二维矩阵K元素km,n值赋予矩阵元素Qw,h,Qw,h=1表示砂岩相,Qw,h=0表示泥岩相,依次类推模拟工区网格G的所有节点,得到复杂几何形态砂体的剖面模型。
作为优选方案,在步骤c中,所述河道砂体横剖面的形状呈“V”型或“U”型。
进一步地,在步骤a中,垂直方向的单元尺寸不大于1米/格。在精细刻画垂向沉积特征,采用这样的网格尺寸,可得到更逼真的模拟结果。
本发明相对于现有技术具有以下有益效果:可根据研究区沉积特征,例如河流相砂体横剖面的透镜状,不受椭圆形状参数来确定河道砂体横剖面形状的约束,建立复杂几何形态的砂体图像并进行布尔模拟,使模拟结果更加真实,提高模型精度,提高油气采收率。
附图说明
图1是一种复杂形态砂体的布尔模拟方法的流程图。
图2是空白模板的示意图。
图3是河道砂体横剖面的形状。
图4是河道砂体横剖面的形状模板。
图5是布尔算法序贯访问工区网格G的流程图。
图6是复杂几何形态砂体的剖面模型。
具体实施方式
下面结合附图详细说明本发明的的实施情况,但它们并不构成对本发明的限定,仅作举例而已,同时通过说明本发明的优点将变得更加清楚和容易理解。
以四川中部上三叠统须家河组沉积储层为例,研究对象为该储层的河道砂体,河道具有透镜状几何形态。
为了刻画复杂形态河道砂体的变化和分布特征,给出河道形态的空白模板,该模板的每个单元在水平方向和垂直方向的尺寸分别为2米/单元、0.5米/单元,空白模板的河道横向宽度为75单元、垂向厚度为10单元,如图1所示。在空白模板中绘制河道砂体横剖面的形状,如图2所示,得到的形状模板。形状模板的所有节点构成维度为75×10的二维矩阵K,当形状模板的节点在河道砂体横剖面的形状内部,矩阵元素km,n=1;当形状模板的节点在河道砂体横剖面的形状外部,矩阵元素km,n=0,如图3所示。
定义工区网格G的水平宽度为600单元、垂向厚度为400单元,工区网格G的所有节点构成定义维度为600×400的二维矩阵Q,所有矩阵元素Qw,h均为空值。
如图5所示,布尔算法序贯访问工区网格G的每个节点,在模拟节点Q(w,h)时,把形状模板中与其对应的所有节点的二维矩阵K元素km,n值赋予矩阵元素Qw,h,Qw,h=1表示砂岩相,Qw,h=0表示泥岩相,依次类推模拟完工区网格G的所有节点,得到复杂几何形态砂体的剖面模型,如图6所示。
其它未详细说明的部分均属于现有技术。
Claims (3)
1.一种复杂复杂形态砂体的布尔模拟方法,其特征在于,包括以下步骤:
a、确定单元尺寸,所述单元尺寸包括沿水平方向和垂直方向的尺寸;
b、建立至少一个空白模板,所述空白模板由多个步骤a所述的单元组成,每个单元构成一个节点,依据河道剖面的宽度、深度分别确定所述空白模板沿水平方向和垂直方向的单元数量M以及N;
c、在步骤b所述的空白模板中绘制河道砂体横剖面的形状,得到形状模板;
d、步骤c所述的形状模板的所有节点构成维度为M×N的二维矩阵K,当形状模板的节点在河道砂体横剖面的形状内部,矩阵元素km,n=1;当形状模板的节点在河道砂体横剖面的形状外部,矩阵元素km,n=0;
e、定义维度为W×H的工区网格G,所述工区网格G的所有节点构成二维矩阵Q,所有矩阵元素Qw,h均为空值;
f、布尔算法序贯访问所述工区网格G的每个节点,在模拟节点G(w,h)时,将步骤d所述形状模板中与其对应的所有节点的二维矩阵K元素km,n值赋予矩阵元素Qw,h,Qw,h=1表示砂岩相,Qw,h=0表示泥岩相,依次类推模拟完工区网格G的所有节点,得到复杂几何形态砂体的剖面模型。
2.根据权利要求1所述的复杂形态砂体的布尔模拟方法,其特征在于步骤c中,所述河道砂体横剖面的形状呈“V”型或“U”型。
3.根据权利要求2所述的复杂形态砂体的布尔模拟方法,其特征在于,步骤a中,垂直方向的单元尺寸不大于1米/格。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510651446.5A CN105204090A (zh) | 2015-10-09 | 2015-10-09 | 复杂形态砂体的布尔模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510651446.5A CN105204090A (zh) | 2015-10-09 | 2015-10-09 | 复杂形态砂体的布尔模拟方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105204090A true CN105204090A (zh) | 2015-12-30 |
Family
ID=54951864
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510651446.5A Pending CN105204090A (zh) | 2015-10-09 | 2015-10-09 | 复杂形态砂体的布尔模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105204090A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109242952A (zh) * | 2018-08-14 | 2019-01-18 | 长江大学 | 一种基于双物源水槽沉积实验的沙体三维建模方法及系统 |
CN110428497A (zh) * | 2019-06-25 | 2019-11-08 | 长江大学 | 辫状河训练图像生成方法 |
CN111985124A (zh) * | 2020-07-15 | 2020-11-24 | 北京金阳普泰石油技术股份有限公司 | 一种单砂体模型构建方法、系统、设备及可读存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080086287A1 (en) * | 2004-12-06 | 2008-04-10 | Shiyu Xu | Integrated Anisotropic Rock Physics Model |
CN102147479A (zh) * | 2011-01-11 | 2011-08-10 | 中国海洋石油总公司 | 一种储层空间物性参数的建模方法 |
CN103679540A (zh) * | 2013-12-31 | 2014-03-26 | 中国石油大学(北京) | 一种基于储层构型要素界面的储层模型的粗化方法及系统 |
CN103699751A (zh) * | 2013-12-30 | 2014-04-02 | 中国石油大学(北京) | 一种基于空间矢量的砂体储层构型建模方法及系统 |
-
2015
- 2015-10-09 CN CN201510651446.5A patent/CN105204090A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080086287A1 (en) * | 2004-12-06 | 2008-04-10 | Shiyu Xu | Integrated Anisotropic Rock Physics Model |
CN102147479A (zh) * | 2011-01-11 | 2011-08-10 | 中国海洋石油总公司 | 一种储层空间物性参数的建模方法 |
CN103699751A (zh) * | 2013-12-30 | 2014-04-02 | 中国石油大学(北京) | 一种基于空间矢量的砂体储层构型建模方法及系统 |
CN103679540A (zh) * | 2013-12-31 | 2014-03-26 | 中国石油大学(北京) | 一种基于储层构型要素界面的储层模型的粗化方法及系统 |
Non-Patent Citations (3)
Title |
---|
乔勇 等: "基于改进布尔模拟的曲流河点坝建模", 《断块油气田》 * |
冯国庆 等: "井间砂体分布的布尔条件模拟", 《石油钻探技术》 * |
尹艳树 等: "储层随机建模研究进展", 《天然气地球科学》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109242952A (zh) * | 2018-08-14 | 2019-01-18 | 长江大学 | 一种基于双物源水槽沉积实验的沙体三维建模方法及系统 |
CN109242952B (zh) * | 2018-08-14 | 2024-04-16 | 长江大学 | 一种基于双物源水槽沉积实验的沙体三维建模方法及系统 |
CN110428497A (zh) * | 2019-06-25 | 2019-11-08 | 长江大学 | 辫状河训练图像生成方法 |
CN110428497B (zh) * | 2019-06-25 | 2021-08-06 | 长江大学 | 辫状河训练图像生成方法 |
CN111985124A (zh) * | 2020-07-15 | 2020-11-24 | 北京金阳普泰石油技术股份有限公司 | 一种单砂体模型构建方法、系统、设备及可读存储介质 |
CN111985124B (zh) * | 2020-07-15 | 2023-07-07 | 北京金阳普泰石油技术股份有限公司 | 一种单砂体模型构建方法、系统、设备及可读存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102609982B (zh) | 空间地质数据非结构化模式的拓扑发现方法 | |
CN108241777B (zh) | 基于非结构网格有限元法计算水合物沉积物中渗流速度场的方法 | |
CN104504754B (zh) | 一种油气储层多点统计建模的方法及装置 | |
WO2017007924A1 (en) | Improved geobody continuity in geological models based on multiple point statistics | |
CN104318618B (zh) | 一种广义三棱柱空间数据模型的三维剖切方法 | |
CN105738952B (zh) | 一种水平井区储层岩石相建模方法 | |
CN103606191A (zh) | 一种复杂地质模型的快速建模方法 | |
CN106227957A (zh) | 等效裂缝建模的方法 | |
CN105204090A (zh) | 复杂形态砂体的布尔模拟方法 | |
CN110599594A (zh) | 一种岩石物性结构三维建模的方法 | |
CN107481320A (zh) | 一种三维地质建模的无网格法 | |
CN104200528A (zh) | 基于矢量闭合的三维建模方法 | |
CN104574513A (zh) | 基于几何形态的准确刻画夹层三维分布的表征方法 | |
CN103207410A (zh) | 一种针对崎岖海底的混合网格模型建立方法 | |
CN105653746A (zh) | 基于嵌入式离散裂缝模型的压裂井建模及模拟方法 | |
Zu et al. | Overview of three-dimensional geological modeling technology | |
CN110765665A (zh) | 一种地学动态建模方法及系统 | |
CN110428497A (zh) | 辫状河训练图像生成方法 | |
CN110306968A (zh) | 不规则井网优化方法及其计算机可读存储介质 | |
CN112562080B (zh) | 一种基于钻孔数据的地质结构降维模型建模方法 | |
CN103236087A (zh) | 一种三棱柱地质模型的构建方法 | |
Li et al. | Alternative 3D modeling approaches based on complex multi-source geological data interpretation | |
Liu et al. | Study on a computing technique suitable for true 3D modeling of complex geologic bodies | |
CN104915990B (zh) | 一种精细刻画河流相储层空间分布的方法 | |
CN115469361B (zh) | 一种碎屑岩地层三维地质建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151230 |