Petrel2018-测井数据导入及基础绘图

[TOC]

Petrel2018-测井数据导入及基础绘图

新建项目

设置项目

设置项目坐标系统

新建井文件夹

导入井头信息

Petrel给出的井头文件示例

# Petrel well head
VERSION 1
BEGIN HEADER
Name
UWI
Well symbol
Surface X
Surface Y
Offset
TD (MD)
Cost
Spud date
Operator
END HEADER
"A 1" A1  3      452905.03 6781834.74 -0.00  2313.76 -999   "11/01/1991" "" 
"B 1" B1  4      456919.15 6787585.65 -0.00  1998.89 -999   "01/11/1995" ""

本次加载示例：

WELLNAME	UWI	KB	TMD	Y1	X1
3430	100120708324W600	686	2307.00 	6229837	00575118
3449	102121708325W600	810.2	4466.00 	6232960	00565650
3466	100022808325W602	778.3	4176.00 	6235099	00567926
3481	100161708325W600	777.5	2605.00 	6231633	00567843
3509	100040908422W600	637.3	1995.00 	6239205	00597766
3510	100062208422W603	0	1939	6242713	00599579

2维窗口显示井位

导入井斜数据

导入后会自动匹配井名

本次导入的井斜数据都是直井。

1 2	WELLNAME UWI KB TMD TVD Y1 X1 Y2 X2 20 100040605424W504 1271.9 3956.0 3910.3 5958323 861287 5958563 861177

如果是斜井、水平井、定向井，需要指定X、Y和TVD。

导入测井数据

给出示例数据：

This loader is completely general and can be used
with different well log formats.
 
Important: When importing multiple data files, make sure all files have the same format.
 
~Version Information Block
VERS .                 2.00:  CWLS LOG ASCII STANDARD - VERSION 2.00
use_wrap .                   NO:  One line per depth step
~Well information block
#MNEM.UNIT        Data Type   Information
#---------- ----------------  -----------
STRT .M             45.8724:  START DEPTH
STOP .M           3889.8576:  STOP DEPTH
STEP .M              0.1524:  STEP
NULL .          -999.250000:  NULL VALUE
COMP .                     :  COMPANY
WELL .                2_1-1:  WELL
FLD  .                     :  FIELD
LOC  .                     :  LOCATION
PROV .                     :  PROVINCE
SRVC .                     :  SERVICE COMPANY
DATE .                     :  LOG DATE
UWI  .                     :  UNIQUE WELL ID
~Curve Information Block
#MNEM.UNIT         API CODE   Curve Description
#---------- ----------------  -----------------
DEPTH.M                    :  
DT   .US/F                 :  DELTA-T
GR   .GAPI                 :  Gamma ray RAY
HPHI .LPU                  :  
ILD  .OHMM                 :  INDUCTION LOG DEEP
RHOB .G/C3                 :  BULK DENSITY
~Parameter Information Block
#MNEM.UNIT            Value   Description
#---------- ----------------  -----------
PROJECT.             ULA_GYDA:  
SET  .                 BEST:  
~A        DEPTH             DT             GR           HPHI            ILD           RHOB 
45.8724      -999.250000         0.2373      -999.250000      -999.250000      -999.250000 
46.0248      -999.250000         0.2083      -999.250000      -999.250000      -999.250000 
46.1772      -999.250000         0.1792      -999.250000      -999.250000      -999.250000 
...
...

只导入目标层目标曲线：

DEPTH       PERM         POR          SW           
954.0000     0.3740       9.7433       32.2979      
954.0508     0.2293       8.5035       38.6421      
954.1016     0.0976       6.3390       53.3102         
...
...

会提示对照井位，注意检查没有出现的井：

数据导入过程会提示有问题的井数据：

数据出错才会报错，需要检查输入数据。

导入层位数据

点选左边Well tops，右键导入数据

层位数据示例：

# Petrel well tops
# Lines starting with # are comments
VERSION 2
BEGIN HEADER
X
Y
Z
MD
TWT
# Valid types: HORIZON, FAULT, OTHER
Type
Surface name
Well name
Symbol
Interpreter
Dip angle
Dip azimuth
END HEADER
457920.31 6784738.68 -1783.84 -1701.53 HORIZON "Base Cretaceous" "B4" Unknown 2068.47 ""   ""   -999   -999   
457746.94 6787092.61 -1886.76 -999     HORIZON "Base Cretaceous" "B8" Unknown 1886.76 ""   ""   -999   -999   
456719.08 6785550.13 -1836.19 -1740.00 HORIZON "Base Cretaceous" "B9" Unknown 1836.97 ""   ""   -999   -999   
453078.12 6786788.35 -1996.92 -1894.14 HORIZON "Base Cretaceous" "C1" Unknown 2003.34 ""   ""   -999   -999   
454686.86 6787607.12 -1977.81 -1878.58 HORIZON "Base Cretaceous" "C2" Unknown 1998.97 ""   ""   -999   -999   
456306.69 6788724.75 -1989.68 -1883.87 HORIZON "Base Cretaceous" "C3" Unknown 2004.29 ""   ""   88.86  -54.66 
454693.38 6786210.63 -1965.23 -1859.14 HORIZON "Base Cretaceous" "C4" Unknown 1980.43 ""   ""   88.23  -70.43

导入层位格式：

Montney	3430	1995.25
Belloy	3430	2262
Montney	3449	2346.4
Belloy	3449	2484.7
Montney	3466	2272
Belloy	3466	2449.5
Montney	3481	2257
Belloy	3481	2536
Montney	3509	1654
Belloy	3509	1961.25
Montney	3510	1636.625
Belloy	3510	1939
...

移动层标识，将层顶底挪入stratigraphy中，自动切分地层。

同理，导入大层和小层，分两个层位文件夹：

工区边界导入

在左侧空白地方右键点击import file

示例工区边界：

FDASCI 0 2 "Computed" 0 99999 0
FDATTR 4 2 1 0
FDATNM 1 "Depth" "real"
!Coordinate System:
!  Projection System ID: 7740565      Initialized  Quick: 0
!  Secondary Horizontal: Units: m  Per/Meter: 1.000000000000
!        Geodetic Datum: ERP50-GB  Ellipsoid: INTNL
!     Projection Method: Transverse Mercator
!      Central Meridian: 3.00000000000
!       User Horizontal: Units: m  Per/Secondary: 1.000000000000
(3D25.16)
-> Cut_0
   0.3354165721485568D+06   0.6459140906319961D+07   0.1086545605468750D+05
   0.3355314799575806D+06   0.6459136397933483D+07   0.1059374804687500D+05
   0.3356463885192871D+06   0.6459131889515877D+07   0.1042337988281250D+05
   0.3357563007507324D+06   0.6459127577129364D+07   0.1012011328125000D+05
...

导入的工区边界：

FFASCI 0 1 "LINES" 0 1e+10
FFATTR 0 1
->1
556327	6212828	0
522316	6285298	0
472123	6488837	0
482725	6488435	0
499919	6487128	0
...

也可以选line-type为：generic boundary

导入后左侧可以看到工区边界

闭合工区边界：不闭合的话不可以进行平面图展布计算

顶部构造图-Z

Z值就是海拔

调整网格大小

选择算法

调整色标

先点击最大最小，然后适当调整色标

地层埋深图-MD

同顶部构造图类似，选择MD属性就可以

等厚图绘制-Thickness

是否使用过滤器，选择no

得到以下点集：

层位名称及点集规范：

创建surface：

孔隙度-渗透率平面图

1、创建孔隙度属性

Well tops → Attributes → 新建属性

选择连续参数

制作POR/SW连续参数

2、相同步骤，制作PERM的连续参数

3、制作层段孔渗SW数值

按enter，变粗后为当前所选

制作POR和PERM的层段数值，注意选择算法为Convergent interpolation

boundary必须闭合才可以进行计算，需要检查边界是否闭合。

提示数据，点击ok。

根据电脑速度运行

运行结束后左侧显示命名的surfacePOR，后自动成图。

4、调整参数-色标

显示刻度尺：
右键 Color Legend → Color Legend

5、出图

点击菜单栏 HOME window 创建 map window

在map菜单栏中可以设置指南针尺度等调整视域大小

中间工具条可以快速调整视域缩放

出图菜单栏map capture

图件导出

右上角capture→export graphic可以导出
选择emf格式最清楚(数据较大)
导出emf后在本地进行裁剪，删去右下角图标
再导出一个png格式，方便数据传输

插值方法说明

地层等厚度用Isochore interpolation，不会有负数
构造、埋深、孔隙度、渗透率、饱和度等用Convergent interpolation，反距离加权平均。

牛眼修改方法

wellpath用-Excel把每一行各导出为一个txt文本

excel→文件→选项→信任中心→信任中心设置→宏设置→启用所有宏

alt+f11调出vb编辑界面

VB代码：

Sub txt()

Dim i, j, arr(), brr(), myRow, myCol

arr = Sheet1.UsedRange

myRow = UBound(arr, 1)

myCol = UBound(arr, 2)

For i = 1 To myRow

Open ThisWorkbook.Path & "\" & arr(i, 1) & ".txt" For Output As #1

Print #1, Join(Application.Index(arr, 1), " ")

Print #1, Join(Application.Index(arr, i), " ")

Close #1

Next

End Sub

运行结果3