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	报废井

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

加载井头文件实操

• 井头源数据

/wellhead

jingtou.txt

井斜（井轨迹）数据文件

井轨迹数据应该是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

加载井斜数据实操

• 井斜源文件

/inclination

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

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

测井数据文件

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格式

加载测井数据实操

• 测井源文件

/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格式

/tools/AscToLas

AscToLas.exe

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

地质分层数据

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

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

地质分层数据实操

/stratification

stratification.txt

加载工区边界

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

加载工区边界实操

• 数值化边界

/boundary

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

/tools/Numerical_boundary

R2V.exe

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

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

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

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

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

选择.xyz格式，点击OK。

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

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

• 添加工区边界源数据

File-expert graphic 导出图片

创建工区

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

截止此步操作，进行保留备份工作空间

data\workspace_backup\importdate_finish

计算器

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

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

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

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

打开一口井查看数据。

单个曲线运算

孔隙度*10。

多个曲线运算

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

联井剖面

选择 1 2 3 井进行联井

勾选孔隙度和泥质含量

调整比例尺

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

调整曲线的刻度

添加分层数据到井段

通过层位拉平

适当调整比例尺

调整井曲线距离

调整曲线颜色，线形

输出图像

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

建模

新建模型

断层模型

没有略

网格平面化

网格化之前确定边界（未边界化）

设置网格大小

查看3D划分

划分标准层

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

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

利用生产边界建立工区

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

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

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

插值后处理

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

设置显示格式

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

同理做其他三个层面

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

双击make horizons

纵向网格划分

属性模型

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

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

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

相模型建立

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

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

• 样式改变

A21面向上平移50m

通过计算器计算

同理A23下移50m

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

调整井名称显示

workspace_backup\k_model_finsih2
k_model_finsih2.pet

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

对岩性曲线进行图像划分

调整后的岩性曲线如下

其他相模型构建

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

岩性相

Property modeling -> Scale up well logs

选择岩性曲线

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

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

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

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

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

返回相模型创建界面

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

设置变参函数

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

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

泥质相

粗化泥质曲线。

设置岩石物理模型

双击Property modeling->Petrophysical modeling

孔隙度相

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

渗透率相

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

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

含油饱和度相

粗化含油饱和度曲线。

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

采用渗透率去约束。

储量计算

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

算储量需要参数

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

净毛比相模型

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

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

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

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

油水边界确定

双击Structural modeling->Make contacts

计算体积

双击Utilities->Volume calculation

• 设置油水界面

• 设置孔隙度

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

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

• 设置含油饱和度

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

• 设置净毛比

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

• 设置油藏属性

计算结果

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

workspace_backup\petrophysical_modling_finsih
petrophysical_modling_finsih.pet

断层模型

fault_model

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

fault_import_data\Fault_Polygon
fault_import_data\Fault_Stick

Fault Stick

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

单条线代表断层切面。

新建Fault Stick模型。

Fault Polygon

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

雷同stick，先新建模型。

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

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

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

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

断层连接

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

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

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

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

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

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

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

调整后连接的断层

断层网格化

点击Processes->Pillar gridding

设置边界

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

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

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

勾选显示工区边界。

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

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

划分网格

双击Processes->Pillar gridding

workspace_backup\fault_model_grid_finish
fault_model_grid_finish

后接创建层面操作。

裂缝模型

基础资料导入

workspace_backup\crack_model_start
crack_model_start

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

导入裂缝数据

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

data\crack_data
imaging_logging_fracture_parameters.prn

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

输入数据。

连井剖面进行显示

显示GR，裂缝参数（主要有玫瑰花图和矢量图），

矢量图显示倾角信息

玫瑰花图显示倾向信息

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

裂缝离散网络模型

生成裂缝密度曲线

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

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

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

人工神经网络预测

人工神经网络功能

• 人工神经网络分类

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

• 人工神经网络预测

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

输入数据

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

选择有裂缝密度的井。

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

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

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

设置参数

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

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

根据预测结果生成

裂缝密度三维曲线

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

直接导入

data\crack_data
ANT(dept)-2.zgy

选择地震数据

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

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

地震数据重采样

将地震数据放入模型。

三维空间裂缝模型

裂缝密度曲线粗化

裂缝密度三维空间

裂缝离散网络模型

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

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