rms数值模拟软件教程

1、RMS 初 级 培 训 教 程北京万格迪信息技术有限公司Beijing Vangand IT Ltd. Co.,前 言RMS是一个功能强大的油藏随机模拟软件，其有灵活多样地数据输入方式和多种模拟方法；该手册旨在对RMS用户进行初级培训。手册中的数据输入和模拟方式以简单的地质情况为例，有关RMS详细的使用方法，请参考用户手册或随机帮助。北京万格迪信息技术有限公司Beijing Vangand IT Ltd. Co.,Tel:+86-10-62321436Fax:+86-10-623471951关于RMS11 RMS功能简介：RMS 是一个较先进的油藏模拟软件。通过精确的储层建模和油藏描述，确定储

2、层的空间分布结构、沉积相空间展布状况，建立储层物性三维数据体、流体分布情况。为提新井位提供直观而可信的参考资料，为老井开采方案的确定提供依据。同时，RMS还能针对有利目标、帮助钻井工程完成精确而科学的钻井方案设计。12 RMS工作流程及模块： 模块 相应功能RMSgeoformmmRMSgeomodRMSgeoplexRMSsimgridRMSstreamRMS wellplanStructural from workFacies modelingPetrophysical modelingsimulation grid design and upsceellingRanking and mo

3、del analysisWell design 13 RMS工作平台工程面板：数据的输入、输出、流程计算、模块运行Menu bar: 菜单条Data folder: 数据文件夹Jobs folder: 作业列表文件夹Workflous : 作业流程设计区域Properfies: 属性窗口多窗口显示面板：图形显示、编辑、图形文件的输出、图形的打印Menu bar: 菜单条Tool bar: 工具条Graphical display area: 显示区域1.4 RMS启动UNIX：%>RMSWindow: 双击RMS 或 Start programs Roxar RMS 15 环境变量设置1

4、51 坐标系统设置 Optiona Project Coordinate system 缺省选项为ANY,即本地区（X，Y，Z）坐标系统 Create/edit coordinate system可设置或编辑新的地理或大地投影坐标系统。152 单位设置 OptionsUnit set 选取 “Project” , 工程单位设置 选取 “Export/import”, 输入输出单位设置 选取 “Both”, 前二者使用相同单位时的设置Add输入新的单位设置名称，允许自定义参数的单位。153 个性化设置 Options Preferences进入User preferences窗口2 础数据的输入

5、21 地震解释层位的输入211 层位名定义 Step1: 产生层位名horizonshorizons admin 层类型： 地震解释层面 由井点或其它数据计算的层面 二界面间的地层 点取“Add/remove/split”按钮，并点取“interpreted horizon”或 “calculated horizon”，创建解释层面和计算层面 Step2：自动插入二层面间的地层 点取“Auto insert” 点取“insert isochores between horizon”Insert Step3: 修改层面或地层名称 直接在name栏内修改即可 点取isochore右边的按钮，选择该

6、层位的顶界面applyOk212 地震解释层面的输入 地震解释结果以（x.y.z）离散点数据文件方式保存，每层一个文件（还有其它格式，输入方法不同）。 以输入层为top的层为例： 文件名：top.dat X1 Y1 Z1 : : : Xn Yn Zn Xn+1 Yn+1 Zn+1 在top下：points（MB3） import ASCII IRAP classic选择相应层文件top.datOK 22 井文件输入 井文件的输入方式灵活，允许多种格式的井位坐标、井径（井斜）数据、测井数据、物性数据、分层数据、沉积相等数据的输入。 以“ASCII RMS WELL”格式为例。 ASCII RMS

7、 WELL文件f36.dat：1.0waterf36 670863.0 653516.86ZONE DISC 0 top 1 f221 2 f221m 3 f222 4 f222m 5 f223 6 f223m 7 f224 8 f224m 9 f2311FACIES DISC 0 background 1 channel 2 mouthbar 3 far_bar 4 sheet_sand 5 overbank DLT unit scaleGR unit scaleSP unit scale670863.04 653516.79 1300.124 0.0000 1 324 12 6 67086

8、3.04 653516.79 1310.183 0.0000 1 240 45 24670863.04 653516.79 1320.241 0.0000 1 245 9 -3 . Well(MB3) import ASCII RMS WELL 选择需要输入的井文件f36.datok 注意：井文件中的层位名要与地层层位名定义一致23 工区边界的输入工区边界文件boundary.dat,即为一多边形X1 Y1 Z1: : :Xn Yn Zn (Z为0)999.000 999.000 999.000 （表示数据结束）Cultural dataintersect lines Project boun

9、dary(MB3) import ASCII IRAP Classic选择相应的边界多边形文件boundary.datok3断层模型的建立31 断层多边形数据的输入例如：有top层断层多边形文件top fault.txt X1 Y1 Z1 : : : Xn Yn Zn -999 -999 -999 horizonsTopMid lines depth(MB3) Import ASCII IRAP Classic 选择top_fault.txt ok 3.2 创建断层模型 Faults（MB3） Create new mode输入一个新的断层模型名：new_mode0133创建断层 faults

10、Input data(MB3) Create new faults输入所要产生的断层名f1,f234 产生断层框架 new_mode01(MB3) generate Network ok3. 5 断层赋值 faults Input data(MB3) Aaaign input data to faults A:选择f1和Fault polygons, 点击Copy selection, 在3D View 窗口中点击相应的断层多边形，add data 再在3Dview中点击bottom层的f2 add data B：再选择f2，重复上述过程Close36 编辑断层框架 将net work设置为编

11、辑状态 net work（MB3） sent to Edit (或直接拖入3D View) A：异常点的编辑 在（选取）状态下，点击要编辑的点，拖到恰当的位置。 B：断层截止点的编辑 点击次级断层的终点MB3extend from node，将鼠标移近主断层。这样，在主层上产生一个次级断层的截止点。 用同样的操作方法，可产生主断层与边界的交点 C：定义断层类型 Faults(MB3) Faut grouping Define groups: 定义所有断层类型 Assign : 指定断层的类型 ok3.7 产生断层模型 new_mode(MB3) Fault modeling general:

12、选取 Fault surfaces 选取 Faul lines 选取 Replace existing surfaces/lines 选择所需断层 fault surface top surface: estimate from horizon data top base surface: estimate from horizon data bottom 分别点击每一个断层、拍定数据类型、算法等参数 Fault lines Sletet horizons: 选择需要产生断层线的层 Import data: 拍定输入层位数据类型：Depth Surface Setup : 点击每个断层，设置外

13、延距离和影响范围 Apply Execute ok4地层构造模型41 计算井与层面的交点horizons（MB3） well intersection选取: Create intersection points using well logs选择所有需要计算井交叉点的层界面和层Append/Replace apply选取: Estimate isochore thicknesses Append/Replace apply ok4.2 层位网格化 horizons（MB3） horizons gridding select horizons: 选择需要网格化的层面 main input: 选择

14、主要输入数据,选取 points Extra input: 选择额外输入数据well points（此参数可选，如果有另外的数据时，可选择产生输入） Faults: 选择断层约束类型,Faults polygons (如果有断层模型) Smooth： 设置平滑参数 Well correction 设置井点数据校正参数 点击“Panels”下左边的按钮，设置网格大小及边界条件 applyexeculteok4.3 地层构造建模 horizons(MB3)Stratigraphic mod select parts of the horizon sequence 中选择: ALL Horizon

15、building 中选取 build from top Isochore cowection 中选取 correct proportion Well cowection 中选取 use applyok 地层模型的建立便完成了！44 将层位模型调整到断层模型中 Horizons(MB3) Horizons administration Adjust4.5 利用断层模型约束的地层构造建模 Horizons strafigraphic mod 其它的选择与4.3内容一致，在option input 中, 选取 use fault model 确定断层模型 ApplyExeculeok5地层厚度数据分

16、析目的：了解各地层的最大厚度、平均厚度、最大厚度出现的概率，为确定地层网格纵向率的大小作准备。51创建分析变量 Data analysis（MB3） Create Univariate 将“well”拖曳到 “Data to analyse”框中. 在select log中选择ZONELOG 选取 Zone log,并选择“ZONELOG” 在Investigate中选择inteval thicknessok便产生一个分析变量A_ZONELOG5. 2 创建柱状图 A_ZONELOG(MB3) Create Histogram 便产生了一个柱状图Histol5. 3 显示柱状图 在3D Mul

17、tiview中，点击柱状图按钮，点击MB1并拖拉鼠标，产生一个柱状图图片框，在数据区点击“Histol”便能显示产生的柱状图。54 设置滤波分层统计地层厚度 Filter(MB3) Create 点击Zone file 并选择相应地层apply。 在3D multi_view中便显示该地层的厚度信息。6建地层网格化数据体三维网格是所有建模的框架61 创建一个网格数据体Zone(MB3) Create Zone选取:Depth Surface例如需要创建三个zone的三维网格数据体：层D到层C，层C到层B，层B到层A。 在地层选择框中点击D、C、B、A OK便产生一个Zone1数据体 Zone1

18、MB3) Inforation可修改 Zone的名称：3DZone 62 产生一个三维地质建模网格 在3DzoneGrid(MB3) Create Modelling Grid Grid Type: Corner point Columns: Increment :25 Rows : Increment :25 选取 Use fault 选择所需断层模型 选取 Auto_Calculate best fit range 在Vertical layout框内，对每一个Zone,分别给出网格厚度或网格化层的数量 apply Execulte ok 6.3 显示网格信息 A：显示网格线 在3Dzon

19、eGrid(MB3) Visual settings 将“lines”选中 B：显示网格信息 在3Dzone下，Grid(mb3) Information 可观察网格大小，最大、最小X、Y、Z等64 创建区域索引参数 A：创建子网格参数 Grid(MB3) Create Region index parameter输入 subgrid 点击“subgrid”execute ok B: 创建断块索引参数 Grid(MB3) Create Region index parameter输入 fault_segment 点击“Grid fault splits”executeok C: 检查索引参数（

20、以fault_segment为例） 点击fault_segment,让其在3Dview中显示 Grid(MB3) Filte 点击“value fileter” 选取fault_segment，上、下二个参数都设置为：3 Apply 可观察3Dview中只具体相应的fault_segment7. 井数据粗化71 产生粗化数据体 Grid(MB3) Create Blocked wells输入“BW”为数据体名Data Selection： wells：选择需要进行粗化处理的井，一般选择建模范围内所有井 将ZONELOG拍定为“Zone log” 将FACIES、poro、 perm、 sw等曲

21、线指定为“scale up logs” parameters： 选择ZONELOG： 选择“shift and scale logs to match surgrids” 设置“cale layer aweraging”为enable 选择FACIES： 设置“call layer averaging”为enable选择poro、 perm、 sw等： 选择“FACIES”作为Bias log apply execute ok7.2 观察井数据粗化结果BW-ZONELOG（MB3） well log editor/caculator选取“show raw logs”,比较粗化前后的数据8确定性

22、储层物性建模 目的：在随机沉积相模拟和随机储层物性模拟之前，简单而快速地确定性模拟有助于提高随机建模精度。 内容：三维内插利用“函数关系”对含水饱和度进行模拟 三维内插: 找到内插椭圆上所有数据通过该点与所有数据点的距离计算加权值利用加权平均法对所有点赋值某些现有地质构造或沉积相数据可用来约束内插运算，一般运用“沉积相”和“断层模型”来约束三维内插 质量控制： 比较粗化井数据与输出数据的最大、最小、平均值。比较值的分布情况（柱状图、趋势）比较在井点处的值的分布情况3Dzone Grid(MB3)CreatIntepulated parametersFramework： 选择poro曲线进行内插

23、 选择“seperate interpolation in each subgrid”,并在subgrid框内并选择“All”Algorithms： 输入内插半径： X Y Z Orientation Subgrid1 4000 2000 20 130 Subgrid2 4000 2700 20 0 Executeok9. 沉积相模拟 在储层建模中，有两大基本内容：沉积相模拟和储层韧性模拟。RMS提供了四种沉积相模拟技术： Elementary: RMSgeomod moduleChannels: RMSgeoplex module Composite: RMSgeoplex module B

24、elts: RMSgeoplex module 模拟质量监控方法：比较输入数据（粗化井数据及趋势数据）和模拟输出数据的最大、最小和平均值；利用柱状图、趋势分析等方法比较输入数据和输出数据的分布状况；比较输出数据体中井点的属性值与井数据；（以Channel modelling为例）91 统计河流沉积相所占比例在所创建的Zone中，粗化井数据体为Bwzone Bwzone (MB3) statistics在selectlog中选择FACIES，分别选择不同的层和井。在“blocked data”中观察河流相的百分比。92 观察河流相的垂直和水平分布趋势Bwzone(MB3) editor/calc

25、ulator在3Dview观察河流相的水平分布趋势93 产生河流相厚度分布柱状图（与5.2步骤相同）94 运行“channel facies modeling”Grid(MB3)stochastic modeling Facies:channels将“the name for the modeling out put”设置为channel打开“Create body parameteter”选择粗化井数据及所有井将FACIES赋给facies log打开“condition on well data”打开“Separate Subgrid modes”并选择“Surgrid 2” （因为 Su

26、rgrid 1无河流相）95 对Subgrid2进行设置打开“sand body mode”将“crevasses”设置为辅助相（auxiliary facies）将channal、 crevasse、 background 分布放入chnnal 、crevasse 和background框中Volum frac. ：Channel: volume fractio:0.4 Tolerauce: 0.02 Crevasse: valume:0.13 Tolerance:0.02在Geometry和 Crevasse/barrier中，分别定义沉积体的厚度和宽度参数。excuteok10储层随机模

27、拟RMS储层韧性模拟有多种实现方式： 确定性三维内插 随机储层韧性模拟 关系函数模拟实现步骤： 数据转换 （data transformation） 将数据转换成满足高斯分布。即期望为0，标准方差为1，没有趋势。变量图分析 （Experimental Variogram） 变量分析表示随着距离的增加，属性值相似程度的变化情况。变量图模拟（Variogram Modelling） 变量图分析结果能不能直接运用到储层物性模拟之中，变量图模拟产生一个数学模型，将这个模型运用到随机模拟之中。随机模拟EXERCISE1: Data Analysis of Petrophysical DataCreate

28、 univariate analysis objects for the blocked porosity and permeability log:_ MB3 Æ Data analysis folder Æ Create Univariate_ MB2 Æ drag and drop the BiasedBW from zone 3Dgeomodel into the Create analysis object panel._ Select Poro from the list of logs and toggle ON Subgrid index_ Tog

29、gle ON Facies log and select FACIES from the list of logs_ Investigate “Scaled up data”_ Press Apply, and a new data analysis object will be created_ Select Perm from the list of logs, and press OK.Open the Data analysis folder (if it is not already open)._ Under Visual settings for both new data an

30、alysis objects toggle on Colour table and change the Variable value from Selection to Facies as in the figure on the next page.Create a histogram for both the Porosity and the Permeability analysis objects:_ MB3 Æ Data analysis Æ BiasedBW Poro ÆCreate Æ Histogram_ Insert a data a

31、nalysis view into the Multiviewer._ Display the histograms in a Data analysis View.Filter out log samples from Subgrid 1:_ MB3 Æ Data analysis Æ BiasedBW Poro ÆFilter ÆCreate Æ Toggle Zone filter and select Subgrid2_ Ok.Create a multivariate analysis of porosity vs. permeabi

32、lity:_ MB3 Æ Data analysis in the Project data list Æ Create MultivariateA new and empty multivariate analysis object will be created. This must be linked to two univariate analysis objects:_ MB3 Æ New Multivariate analysis object Æ Connect to Univariates_ Drop the porosity and p

33、ermeability univariates into the drop site.These univariate analysis objects are now linked to the multivariate analysis object._ Open the Visual settings panel for the multivariate analysis object, and toggle on Colour table. Change the Variable value from Selection to Facies to map a colour table

34、to the facies distribution. Then create a scatter plot of blocked porosity vs. permeability:_ MB3 Æ Multivariate analysis Æ Create Æ ScatterplotSelect porosity and permeability._ Create a Data Analysis View and display the scatter plot in it._ Plot the regression line by clicking on t

35、he Regression line icon in the Multiviewer panel._ Save the project.Optional Step_ Filter on each sand facies type one by one to see whether the correlation between porosity and permeability is facies dependent.EXERCISE2: Transformation of Petrophysical DataCreate the following scatter plots for the

36、 porosity and permeability data analysis objects made in a previous exercise:Porosity vs. X (simbox), call it SimX.Porosity vs. Z (simbox), call it SimZ.Porosity vs. Z (intrabody in simbox), call it SimZbody._ MB3 Æ Data analysis Æ BiasedBW Poro Æ Create Æ Scatterplot_ Once creat

37、ed give the objects relevant names._ MB3 Æ Data analysis Æ BiasedBW Poro Æ Scatter 1 Æ Information.Optional:Do the same for permeability:Permeability vs. Z (simbox)Permeability vs. X (simbox)Permeability vs. Z (intrabody in simbox)In this step we will transform the blocked porosi

38、ty data for each facies in Subgrid 2 so that they are normally distributed. Filter on the subgrid and facies as indicated below. The suggested transformations steps are given below, actual values may differ slightly when the data analysis is carried out on live data. For all facies, apply a Truncate

39、 realizations transformation to limit the porosity to an acceptable and reasonable data range.Subgrid 2; channelFilter on the channel facies type in Subgrid 2:_ MB3 Æ Data analysis Æ BiasedBW Poro Æ Filter Æ Edit. Æ Toggle Facies filter and select channel Æ OkOpen the t

40、ransformation panel:_ MB3 Æ Data analysis Æ BiasedBW Poro Æ Transf Æ Create._ Select Truncate realizations in the transformation list._ Estimate the range (Press Estimate)._ Toggle on Truncate smaller values & Truncate larger values, and enter 0.0 and 0.35 respectively._ Add

41、this transformation; Press Add button._ Keep the transformations panel open throughout the exercise to view the various transformation steps as they are added to the transformation sequence for each facies._ Display the scatter plots per facies in turn in a Data Analysis View.An intrabody trend is a

42、pparent within the channel facies, however well_F contains an incorrect interpretation of either the petrophysical or facies logs. _ Apply a well filter in the Data analysis to disregard this data point._ Keep the transformations panel open throughout the exercise to view the various transformation

43、steps as they are added to the transformation sequence for each facies._ Display the scatter plots per facies in turn in a Data Analysis View.An intrabody trend is apparent within the channel facies, however well_F contains an incorrect interpretation of either the petrophysical or facies logs. _ Ap

44、ply a well filter in the Data analysis to disregard this data point.Subgrid 2; crevasse_ Filter on the crevasse facies type in Subgrid 2:_ MB3 Æ Data analysis Æ BiasedBW Poro Æ Filter Æ Edit._ Use the Truncate realizations transformation to limit the simulation to an acceptable a

45、nd reasonable data range._ Toggle on Truncate smaller values & Truncate larger values, and enter 0.0 and 0.16 respectively._ Check for trends in the scatter plots. No trends are apparent.Display the histogram in the other Data analysis view._ Use the Truncate observations transformation to remov

46、e high outlier data._ Estimate the range (Press Estimate)._ Toggle on Truncate larger values and enter 0.17._ Add this transformation; Press Add button._ Apply a Shift and scale transformation in the histogram view.Check that the mean value is 0.0 and the standard deviation close to 1.0.Subgrid 2; b

47、ackground_ Filter on the background facies type in Subgrid 2:_ MB3 Æ Data analysis Æ BiasedBW Poro Æ Filter Æ Edit._ Check for trends in the scatter plots. No trends are apparent._ Use Truncate realizations transformation to limit the simulation to an acceptable and reasonable da

48、ta range._ Estimate the range (Press Estimate)._ Toggle on Truncate smaller values & Truncate larger values, and enter 0.0 and 0.15 respectively._ Add this transformation; Press Add button._ Display the histogram in the other Data Analysis View and apply a Normal score transformation.Check that

49、the mean value is 0.0 and the standard deviation close to 1.0.Once completed, remove all filters and look at a summary of all the transformations:_ Open the transformations panel, if it is not already opened_ MB3 Æ Transf Æ Edit _Press the Transformation summary button_ Close the panel.Optional stepTransformations also need to be applied to permeability for each facies type: _ Repeat the steps above using the data analysis to guide the choice of transformations. (The transformation sequences will be similar.)The permeability distributions are normally significan