N23建模步骤.doc

Petrel地质建模解析步骤（东胜N23A项目）一、整理数据，格式为不带符号的空格间隔，*.txt格式，需要资料：小层平面图数据、构造图数据、井对比分层数据、井斜数据、孔渗饱数据，砂体等厚图数据。1)井口坐标wellname y x kb 顶深(MD) 底深(MD)s1025120924.0021674358.00 4.60 1000 1900s1025120924.0021674358.00 4.60 1000 19002)井斜数据MD X Y Z TVD DX DY AZIM INCL1499.878 456979.063 6782712.412-1499.878 1499.878 0.000 0.000 99.853 42.2771500.031 456979.164 6782712.395 -1499.991 1499.9910.100-0.017 99.852 42.2781500.183 456979.265 6782712.377 -1500.104 1500.104 0.201 -0.035 99.851 42.281 1500.335 456979.3666782712.359 -1500.217 1500.2170.302-0.05299.850 42.2833)测井曲线DEP(MD) RESIS AC SP GR1400.0000 5.1703 374.2136 35.5975 127.01400.1000 5.2997 374.2136 35.7233 127.01400.2000 5.0606 372.9888 35.8568 126.8如有测井综合解释的孔、渗、饱资料，按相同格式加上。可直接读取测井*.las格式文件。4)沉积相的划分,或有效厚度，孔隙度，饱和度等深度(MD)孔隙度 饱和度 渗透率 岩性代码 有效厚度1046.3 0.22 0.42 3000 1 0.7 1047.2 0.27 0.53 9000 2 0.21047.6 0.68 0.22 15000 3 1 1048.6 0.32 0.32 12000 1 0.3准备相或孔渗饱曲线时,一般没有以上格式的曲线,但我们可以利用现成的数据库通过编写程序来实现,准备以上曲线时要以每口井的名称为文件名来描述该井的相或孔渗饱信息.这样可以通过批量输入来节省数据的输入时间.地质数据1)分层数据XY Z(TVD) Wellpoint 层名井名s10221674358.00 5120924.00 -1296.20 horizon pzs10221674358.00 5120924.00 -1315.80 horizon pbs10221674358.00 5120924.00 -1315.80 horizon pb s102或采用以下的格式:Z(MD) Wellpoint 层名 井名1272.70 horizon pd s102 1296.20 horizon pz s102 1315.80 horizon pb s102 2)等厚图（点或面,也可以利用PETREL计算 Isochors Point）X Y thickness 21674358.00 5120924.00 23.70 21674358.00 5120924.00 33.20 21674358.00 5120924.00 44.80 3)属性平面图（点，面）（N/G, Porosity, Permeability, Saturation）X Y property 21674358.00 5120924.00 23.70 21674358.00 5120924.00 33.20 21674358.00 5120924.00 44.80 4)断点数据（点）X Y z(TVD) 21674358.00 5120924.00 1123.70 21674358.00 5120924.00 1133.20 21674358.00 5120924.00 1144.80 断点数据可以通过General Point /line的方式输入,输入后要二、数据输入1、构造、零线、歼灭线、断层、边界线、井分层数据输入设定输入井孔渗饱的深度选择输入数据的类型，构造选点，断层、零线、歼灭线选线型输入断层线，在setting里设置成断层性质断层赋值到构造面上（先做构造面的生成，参考下面建模/构造面的生成），点击蓝表输入相应构造面，点击Z=A，点击Apply，赋值成功然后在下窗口make polygons，点击“剪子功能”剪出上下盘点击相应断层，右键Split，劈开断层上下盘三、建模1、断层建模打开Structual Modeling/Define Model,键入模型名字，点Ok在3D窗口中打开劈开的断层线将窗口转到模型点击右边工具“箭头”，shirt选中断层线点击断层成型图标，生成断层剖面，点颜色填充用构造面来核对断层的精确性2、构造面的生成点击Utilities/MakeSurface，点蓝标输入构造线点，加入边界boundary在Geometry，选中Automatic（自动），Grid increment（区块大小设置）填写参数点击Algorithm，选择图形属性MethodConvergent Interpolation（计算方法：收敛），参数设置如下选择井控制，有效单井控制距离800成图3、地层的生成（1）、在Welltop里剪出地层，全选tops，先点击顶，后在底右键Convert to iso,转化为厚度点在生成文件中点开，选thick右键Convert to point所要砂层，选中输入歼灭线，显出歼灭线将尖灭线数字化，Make polygons，顺着歼灭线画点界，完毕双击，“yes”点击“箭头”选中刚才所画的线，点击圈闭，生成可赋值的歼灭线将先生成的歼灭线右键，convert to point点击生成的点，右键setting Calculations A=0.4（也可以不赋值）点击相应的厚度面，Setting /Operations/Common Operations/Append Points 点击蓝标选入方才的歼灭线点值 ，点击run生成地层Make surface，点蓝标，加入厚度点，边界选中成图属性，选Iso厚度面计算方法不用井控制生成砂体展布图3、制造网格 Pillar Gridding，填写I、J参数，点T生成沉积方向点击制作沉积趋势方向沉积方向如图建立以断层为界的边界线，先点中断层，选择右边“B”，点线点中断层后点击在断层的起点点击，然后顺着边界成线到断层的另一头的第一点。点中断层的始终点4、生成层面，Make Horizons 点蓝标一次输入构造面以及构造面的深度（well top）断层断距参数设置井控制范围参数生成层面5、生成纵向地层，Make zones，一次输入构造面，以及构造面的顶、底按大层骨架输入地层厚度面及相应顶/底6、生成层内渗流网格（地层分辨率）（按砂体厚度制定网格大小，干层、隔层为1，比如1m的砂层可以定2个网格）7、输入测井曲线选择线属性选择线对应的英文名成图测井曲线离散化 Scale up well logs 选择孔渗饱曲线成图8、数据分析数据分析的结果可以直接被相建模和属性建模的模块调用.数据分析分为两类:对离散数据的分析和连续数据的分析.离散数据的分析:（1）.打开分析窗口,界面如下:（2）.可以选择分析的对象,是经过离散化的井点,还是未离散化的测井曲线,或整个模型.（3）.是否使用滤波功能（4）.分别按以下的标签,进行相应的分析（5）.打开每个标签后,按键,刷新显示（6）.在进行probility分析,需要一个连续型变量的属性模型,其目的是分析该离散属性和连续属性的相关性并求出其对应关系.如下图下图的意义是:coarse sand(红色部分) 在属性(Attribution)为1的地方出现的概率为55%,而在其他属性的地方出现的概率为0;sand and conglomerate(兰色部分)在属性(Attribution)为0的地方出现的概率为90%,在属性为1的地方出现的概率为5%.（7）.变差函数的分析1)首先设置主方向的分析参数,包括带宽,搜索半径,步长,容忍度等,然后再设置次方向和垂向上的参数,这些参数的设置需要用户对本地区数据的大概了解的基础上,否则分析的结果的可信度大大降低.在该例中的分析参数和结果如下图.在分析变差之前,首先大概了解数据的分布情况,然后再设置这些分析参数,这样才能达到比较好的分析效果.连续数据的分析（1）.打开分析窗口,界面如下:（2）.可以选择分析的对象,是经过离散化的井点,还是未离散化的测井曲线,或整个模型.（3）.是否使用滤波功能;是否用相约束（4）.分别按以下的标签,进行相应的分析（5）.打开每个标签后,按键,刷新显示（6）.在进行transformation分析时,可以能够多种数据的处理,包括,输入截断,输出截断,对数变换,奇异值的消除,1D,2D,3D趋势分析,正态分布变换等（7）.例如要进行输入截断和正态分布转换,处理的对话框如下:（7）.变差函数的分析（同离散数据的分析相似）1)首先设置主方向的分析参数,包括带宽,搜索半径,步长,容忍度等,然后再设置次方向和垂向上的参数,这些参数的设置需要用户对本地区数据的大概了解的基础上,否则分析的结果的可信度大大降低.在该例中的分析参数和结果如下图.在分析变差之前,首先大概了解数据的分布情况,然后再设置这些分析参数,这样才能达到比较好的分析效果.孔渗饱饱和度场9、生成有效厚度层面，用零线，制造净毛比场，过程同生成砂体层面选择刚生成的零线作为边界，赋值整个面为1值同理生成零面两个层面相交生成平面相8、生成空间平面属性面场 Faces Modeling，选择层，选择计算方法Assign Values 有有效储层的面选择surface，点