pertrel地质建模-实操

pertrel地质建模-实操

pertrel地质建模-实操

[TOC]

HibisciDai/Petrel2009_testdata

数据输入

数据准备

wellheader(井头文件)

wellheader井头文件,包括

井名,Wellname;
x坐标,X-coord;
y坐标,Y-coord;
井底深度,KB;
补心海拔,Td(MD);
井别(油井,气井,注水井等用户可以自己定义),Symbol。

其中WellName(井名)是必须的且井名必须唯一,X、Y坐标也是必须的,其余为可选项。

Wellname X-coord Y-coord KB Td(MD) Symbol 备注
X1 239749.5300 4581430.2300 944.87 5754.00 2 干井
X2 245019.6500 4582221.3600 942.11 5370.00 3 油井
X3 246821.7400 4582047.7900 941.32 5750.00 4 低产油井
X4 244804.0900 4583215.8600 941.98 5580.00 3 油井
X5 243662.9300 4584608.3500 942.70 5602.00 3 油井
X6 243397.0300 4580738.9900 944.42 5633.65 3 油井
X7 241839.9400 4582251.2300 937.55 5620.00 15 注水井
X8 15241839.9400 4582251.2300 937.55 5826.00 19 报废井

每列的位置比较自由,没有固定的顺序。

加载井头文件实操

  • 井头源数据
1
2
3
/wellhead

jingtou.txt

加载井头文件实操-1

加载井头文件实操-2

加载井头文件实操-3

加载井头文件实操-4

加载井头文件实操-5

井斜(井轨迹)数据文件

井轨迹数据应该是ASCII码文件,按列排放。数据中包含几种可能的信息组合格式:

  • ①Measured depth(测深MD),inclination(井斜角),azimuth(方位角);
MD AZIM INCL
1499.88 99.85 42.28
1500.03 99.85 42.28
1500.18 99.85 42.28
1500.34 99.85 42.28
1500.49 99.85 42.29
1500.64 99.85 42.29
1500.79 99.85 42.29
  • ②True vertical depth(真垂深TVD), X-offset(X偏移量), Y-offset(Y偏移量) (MD可选) ;
MD(米) TVD(米) DX(米) DY(米)
0 0 0 0
0.125 0.125 -0.001 -0.004
0.25 0.25 -0.002 -0.007
0.375 0.375 -0.004 -0.011
0.5 0.5 -0.005 -0.014
  • ③True vertical depth(真垂深), X, Y (MD可选)
MD(米) TVD(米) X坐标(米) Y坐标(米)
0 0 15703611.07 4149596.08
0.125 0.125 15703611.07 4149596.083
0.25 0.25 15703611.07 4149596.086
0.375 0.375 15703611.07 4149596.089
0.5 0.5 15703611.07 4149596.092
0.625 0.625 15703611.08 4149596.094

加载井斜数据实操

  • 井斜源文件
1
2
3
4
5
6
7
8
9
10
/inclination

W1.txt
W2.txt
W3.txt
W4.txt
W5.txt
W6.txt
W7.txt
W8.txt

雷同加载井头文件实操-2导入

加载井斜文件实操-1

加载井斜文件实操-2

加载井斜文件实操-3

加载井斜文件实操-4

测井数据文件

Petrel软件对输入的测井曲线格式是灵活的,唯一的要求是测井数据文件按列排列,即每条测井曲线为列数据排列。
Petrel软件支持LAS格式文本格式的测井曲线数据。

  • 文本格式
DEPTH AC CNL DEN GR
5180 49.091 -0.2 2.745 33.432
5180.125 49.373 -0.2 2.741 33.375
5180.25 49.395 -0.2 2.74 32.041
5180.375 49.241 -0.1 2.743 32.195
5180.5 48.854 -0.1 2.746 31.385
5180.625 48.325 0 2.749 33.124
  • LAS格式

LAS格式测井数据

加载测井数据实操

  • 测井源文件
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
/logging

W1.txt
W2.txt
W3.txt
W4.txt
W5.txt
W6.txt

/logging_las

W1.LAS
W2.LAS
W3.LAS
W4.LAS
W5.LAS
W6.LAS
  • ascii转las格式
1
2
3
/tools/AscToLas

AscToLas.exe

ascii转las格式-1

ascii转las格式-1

雷同加载井头文件实操-2导入

加载测井文件实操-1

加载测井文件实操-2

加载测井文件实操-3

加载测井文件实操-4

地质分层数据

分层数据是指一个层位或者一个断层与井轨迹发生交叉,在井上的一系列标记。

MD Type Horizon Name Well Name
2068.47 HORIZON Base Cretaceous B4
1886.76 HORIZON Base Cretaceous B8
1836.97 HORIZON Base Cretaceous B9
2003.34 HORIZON Base Cretaceous C1
1998.97 HORIZON Base Cretaceous C2
2004.29 HORIZON Base Cretaceous C3
1980.43 HORIZON Base Cretaceous C4

地质分层数据实操

1
2
3
/stratification

stratification.txt

地质分层数据实操-1

地质分层数据实操-2

地质分层数据实操-3

加载工区边界

加载工区边界-1

需要对工区数据进行数值化。

加载工区边界实操

  • 数值化边界
1
2
3
4
5
6
7
8
9
/boundary

boundary.bmp | 工区边界图
boundary.prj | 数值化处理后项目
boundary.xyz | 工区边界数据

/tools/Numerical_boundary

R2V.exe

采用R2V数值化软件,点线数值化,打开bmp格式工区边界图片。

数值化边界-1

数值化边界-2

需要确定3-4个定位点作为工区基准,图中可以看到边界基准点的焦点。

数值化边界-3

雷同操作,进行四个基准点的添加。

数值化边界-4

通过四个基准点,可以确定改平面任一一点的位置,下边进行边界划分。

依次按照边界进行鼠标点击,连接成线,点位要按顺序点击。

数值化边界-5

数值化边界-6

选择.xyz格式,点击OK

边界数据是x,y,z格式,制作时候可以选择z值。

用记事本打开数据,要进行闭合,起始点和结束点必须相同

数值化边界-7

  • 添加工区边界源数据

加载工区边界实操-1

加载工区边界实操-2

加载工区边界实操-3

加载工区边界实操-4

加载工区边界实操-5

File-expert graphic 导出图片

加载工区边界实操结果图

创建工区

通过选择菜单栏的 File: Save Project ,Ctrl+S 或者单击Tools Bar上面的 Save Project 图标存储当前工区,软件会弹出“Save project as”窗口,提示输入工区名字,输入工区名称(如下图),OK。

创建工区

1
2
3
截止此步操作,进行保留备份工作空间

data\workspace_backup\importdate_finish

计算器

计算全局就点击well,计算单口井就选择单口井。

根据泥质含量去计算砂岩泥岩区分

计算器-1

计算器-2

if(SH<=40%,shayan,niyan)

计算器-3

点击enter,如果有井没有SH曲线则会提示错误,点击确定跳过。

计算器-4

打开一口井查看数据。

计算器-5

计算器-6

单个曲线运算

孔隙度*10。

计算器-7

多个曲线运算

新增曲线然后用曲线相乘,计算储量 = 孔隙度 * 含油饱和度

计算器-8

联井剖面

联井剖面-1

联井剖面-2

选择 1 2 3 井进行联井

联井剖面-3

勾选孔隙度和泥质含量

联井剖面-4

调整比例尺

联井剖面-5

调整曲线顺序,井深度曲线在中间,左边右边各一条曲线

联井剖面-6

调整曲线的刻度

联井剖面-7

联井剖面-8

添加分层数据到井段

联井剖面-9

通过层位拉平

联井剖面-10

适当调整比例尺

联井剖面-11

调整井曲线距离

联井剖面-12

调整曲线颜色,线形

联井剖面-13

输出图像

联井剖面-14

注意调整输出图像页面大小

联井剖面-15

联井剖面-16

建模

建模流程

新建模型

新建模型-1

断层模型

没有略

网格平面化

网格化之前确定边界(未边界化)

网格平面化-1

设置网格大小

网格平面化-2

查看3D划分

网格平面化-3

划分标准层

利用井数据进行构造形态划分。

做构造面有两种方式,地震数据或者井数据

利用生产边界建立工区

划分标准层-1

划分标准层-2

该层面数据来自于8口井插值法

划分标准层-3

插值以后改变了井点的原来数据,所以要采用矫正,使用A21层去矫正。

保证插值以后层面数据和井上数据是一致的。

划分标准层-4

插值后处理

划分标准层-5

  • Unchange 不改变
  • Truncated 截断(设置高度),%表示输入数据的百分比。
    10%,插值后如果超过输入数据10%,做截断操作。
    10,最大值不超过10,超过=10。
  • Eliminated 删除
    超过直接删除

设置显示格式

划分标准层-6

划分标准层-7

可以看出是一个单斜构造面

同理做其他三个层面

划分标准层-8

四个层面已经输入,需要导入至模型中

双击make horizons

划分标准层-9

划分标准层-10

纵向网格划分

纵向网格划分-1

纵向网格划分-2

属性模型

模型建立需要对测井曲线进行粗化

模型通过井的数据差值、数学算法得到井间数据

为了更接近需要用更多信息进行约束(数学运算)

相模型建立-1

相模型建立

属性模型建立-1

属性模型建立-2

主要依据井数据,通过数学方法得到数据建立模型。为了让模型更接近实际,需要用其他资料进行约束。

比如在建立碎屑岩地层,建立渗透率模型可以采用孔隙度来约束。孔隙度高的地方渗透率也会变高。通过井上的孔渗数据进行分析。

属性模型建立-3

属性模型建立-4

属性模型建立-5

属性模型建立-6

属性模型建立-7

属性模型建立-8

  • 样式改变

A21面向上平移50m

属性模型建立-9

通过计算器计算

属性模型建立-10

属性模型建立-11

同理A23下移50m

属性模型建立-12

井在目的层之上的段去掉(不可返回操作)

属性模型建立-13

调整井名称显示

属性模型建立-14

1
2
workspace_backup\k_model_finsih2
k_model_finsih2.pet

同理对砂泥岩曲线进行粗化

属性模型建立-15

对岩性曲线进行图像划分

属性模型建立-16

调整后的岩性曲线如下

属性模型建立-17

其他相模型构建

建立相模型前要进行测井曲线的粗化,保证一个网格一个数据。

岩性相

Property modeling -> Scale up well logs

选择岩性曲线

岩性相模型-1

粗化后的曲线有[U]标识,只是测井曲线上的粗化。

岩性相模型-2

井上边可以看到代表岩性的色块。

岩性相模型-3

双击Propetry modeling->Facies modeling,相模型

岩性相模型-4

采用截断高斯函数,对于岩相来说只有两种格式,泥岩和砂岩。

岩性相模型-5

需要先对岩性曲线进行修改

岩性相模型-6

岩性相模型-7

岩性相模型-8

返回相模型创建界面

岩性相模型-9

高斯截断对每一个相进行单独设置

设置变参函数

岩性相模型-10

由于没有约束条件,所以仅仅是简单的数学运算,没有体现地质规律。

一般用地震或者平面沉积相组合来进行约束。

岩性相模型-11

泥质相

粗化泥质曲线。

泥质相模型-1

设置岩石物理模型

双击Property modeling->Petrophysical modeling

泥质相模型-2

泥质相模型-3

孔隙度相

孔隙度相模型-1

孔隙度相模型-2

采用相控比较符合地层规律。

孔隙度相模型-3

渗透率相

粗化渗透率曲线,暂时选用算数平均法。

渗透率相模型-1

渗透率相模型-2

渗透率相模型-3

结果与孔隙度相模型有一定的相关度。

渗透率相模型-4

含油饱和度相

粗化含油饱和度曲线。

含油饱和度相模型-1

含油饱和度与孔隙度渗透率也有关系,原始含有饱和度跟储层物性有必然相关性。

采用渗透率去约束。

含油饱和度相模型-2

含油饱和度相模型-3

含油饱和度相模型-4

储量计算

基于三维网格计算,将有效网格储量累加起来,就是油藏的储量。

算储量需要参数

1.面积;
2.储层有效厚度;(网格高度是有效厚度)
3.孔隙度;
4.原始含油饱和度;
5.原油的密度;(需要输入)
6.原油的体积系数。(需要输入)

净毛比相模型

需要判断哪些网格有效。
即孔隙度、渗透率、含油和饱和度大于某个值,网格才有效。

净毛比,判断为储层的网格,值为0或者1。

生成净毛比模型,先复制一个相模型。

储量计算-1

计算器计算,注意单位类型。

储量计算-2

储量计算-3

油水边界确定

双击Structural modeling->Make contacts

储量计算-4

储量计算-5

计算体积

双击Utilities->Volume calculation

  • 设置油水界面

储量计算-6

  • 设置孔隙度

计算储量时候孔隙度的单位是小数,目前模型为百分数。

采用计算器对 孔隙度/100。

储量计算-8

  • 设置含油饱和度

采用计算器对 含油饱和度/100。

储量计算-9

  • 设置净毛比

如果没设置则会计算全部。

储量计算-7

  • 设置油藏属性

储量计算-10

计算结果

储量计算-11

结果未精确到小数点后,可以采用 % 方式计算最后结果 / 100即可

1
2
workspace_backup\petrophysical_modling_finsih
petrophysical_modling_finsih.pet

断层模型

1
fault_model

断层数据有两种类型
断层序Fault Stick和断层多边形Fault Polygon

1
2
fault_import_data\Fault_Polygon
fault_import_data\Fault_Stick

Fault Stick

Fault_Stick-1

Pertrel中井文件需要有专用文件夹,其他数据都可以自定义文件夹。

Fault_Stick-2

Fault_Stick-3

Fault_Stick-4

单条线代表断层切面。

新建Fault Stick模型。

Fault_Stick-5

Fault_Stick-6

Fault_Stick-7

Fault_Stick-8

Fault Polygon

数据来源于构造图,断层在构造图上的投影。

Fault_Polygon-1

Fault_Polygon-2

雷同stick,先新建模型。

Fault_Polygon-3

点击右侧按钮Create fault from fault polygons从断层多边形产生断层模型

Fault_Polygon-4

默认形成的断层形状,在两边有畸形体现。一般需要把两边的断层删掉。

Fault_Polygon-5

Fault_Polygon-6

Fault_Polygon-7

Fault_Polygon-8

利用外边界区域扩展范围。选中边界stick后点击右侧工具栏Add pillar to end

Fault_Polygon-9

Fault_Polygon-10

如果对断层stick间隔感到太宽,需要在两个柱子之间增加一个。多选两个柱子然后点击右侧工具栏Add pillar between

Fault_Polygon-11

Fault_Polygon-12

断层连接

勾选Fault Stick和Fault Polygon两个模型,调整显示格式。

断层连接-1

Polygon断层柱子有超过Stick模型,所以需要先删掉Polygon断层与Stick相切部分的边界柱子。

断层连接-2

观察相交趋势面,连接两个断层的stick。(不是同一个断层模型中的断层不可以链接)

断层连接-3

断层连接-4

断层连接-5

公用的stick变为纯色,两个断层链接成功。

断层连接-6

Polygon断层有穿过Stick断层,需要调整弯曲度。

对于当个柱子的调整需要选中整个柱子然后拖拽上下移动。

单选一个点可以进行上下移动。

断层连接-7

调整后连接的断层

断层连接-8

断层网格化

点击Processes->Pillar gridding

断层网格化-1

设置边界

设置断层边界,将断层作为边界的一部分。

选中断层边界,右边工具栏点击Set part of grid boundary

断层网格化-2

断层网格化-3

其他边界参考工区边界,手动画。

勾选显示工区边界。

断层网格化-4

点击右侧工具栏Create external grid boundary segment

断层网格化-5

从断层最末点开始点击,沿着工区边界划线,再到断层起始点结束。

断层网格化-6

划分网格

双击Processes->Pillar gridding

断层网格化-7

断层网格化-8

1
2
workspace_backup\fault_model_grid_finish
fault_model_grid_finish

后接创建层面操作。

裂缝模型

基础资料导入

1
2
workspace_backup\crack_model_start
crack_model_start

初始模型中已经将常规的测井资料导入

裂缝模型-1

导入裂缝数据

输入数据为成像测井解释裂缝数据,深度裂缝倾角裂缝倾向

1
2
data\crack_data
imaging_logging_fracture_parameters.prn

插入井点数据文件夹(输入裂缝资料属于点数据)。

裂缝模型-2

输入数据。

裂缝模型-3

裂缝模型-4

裂缝模型-5

连井剖面进行显示

裂缝模型-6

显示GR,裂缝参数(主要有玫瑰花图矢量图),

矢量图显示倾角信息

裂缝模型-7

裂缝模型-8

玫瑰花图显示倾向信息

裂缝模型-9

裂缝模型-10

倾角和倾向主要呈90度关系,可以对应验证裂缝信息。

裂缝离散网络模型

生成裂缝密度曲线

需要裂缝密度参数。需要根据成像测井数据来计算裂缝密度。

裂缝模型-11

根据输入裂缝点数据来生成裂缝密度。

裂缝模型-12

对于没有输入裂缝点的井,可以通过对已有裂缝资料的井建立样品,用人工神经网络方法进行预测。

人工神经网络预测

人工神经网络功能

  • 人工神经网络分类

例如对于成像测井,我们输入伽马,自动输出分类后的伽马曲线。

  • 人工神经网络预测

输入声波、密度、中子可以计算预测渗透率。

输入数据

输入数据,选择输入的数据,通过测井数据建立关系。

选择有裂缝密度的井。

  • 裂缝与那些测井曲线相关

GR(岩性)、声波、密度、电阻率、深浅电阻幅度差(设置新曲线幅度差、一般俩个电阻取对数再减)、井径

裂缝模型-13

生成新曲线RDRS,测井曲线深浅电阻幅度差。

设置参数

裂缝模型-14

裂缝模型-15

裂缝模型-16

误差满足或者迭代次数达到就退出训练。

预测完成后会生成非线性的关系式,软件中看不到。

根据预测结果生成

裂缝模型-17

裂缝模型-18

裂缝模型-19

裂缝密度三维曲线

导入蚂蚁体数据(地震数据)

直接导入

1
2
data\crack_data
ANT(dept)-2.zgy

选择地震数据

裂缝模型-20

裂缝模型-21

裂缝模型-22

此时输入的蚂蚁体数据并没有写入模型当中。

需要对地震数据进行重采样来对应模型网格大小。

地震数据重采样

裂缝模型-23

将地震数据放入模型。

裂缝模型-24

三维空间裂缝模型

裂缝密度曲线粗化

裂缝模型-25

裂缝密度三维空间

裂缝模型-26

裂缝模型-27

裂缝模型-28

裂缝离散网络模型

裂缝模型-29

裂缝模型-30

根据裂缝密度来进行确定裂缝平均长度,约小的长度,最终结果越密集。

裂缝模型-31

裂缝模型-32

最重要是确定裂缝与那些因素有关,用参数去尽可能的约束。

文章作者: HibisciDai
文章链接: http://hibiscidai.com/2020/11/05/pertrel地质建模-实操/
版权声明: 本博客所有文章除特别声明外,均采用 CC BY-NC-SA 4.0 许可协议。转载请注明来自 HibisciDai
支付宝打赏
微信打赏