1、1Eclipse 数值模拟软件问答(初级)1. ECLIPSE输出结果文件是哪些?.GRID或.FGRID: 网格文件.EGRID: 网格文件,与 GRID格式不同,文件要小的多。 (用关键字 GRIDFILE来控制输出类型).INIT 或.FINIT: 属性文件。 (用关键字 INIT来控制输出).PRT: 报告输出。 文件很大,模型处理及计算结果详细报告。(RPTGRID,RPTPROP,RPTSOL,RPTSCHED 控制输出).LOG: 后台作业时的输出报告,文件比 PRT要小很多。可用于错误检查。.DBG: Debug文件,一般不用。可用于检查 ECLIPSE如何处理输入参数。.SA
2、VE: 用于快速重启。 (用关键字 SAVE来控制输出).RFT :RFT 计算结果。 (用关键字 WRFTPLT来控制输出).FLUX: 流动边界。 (用关键字 DUMPFLUX来控制输出).Snnnn 或.UNSMRY: 图形文件输出(在 SUMMARY部分定义).Xnnnn 或.UNRST: 重启文件输出(用 RPTRST,RPTSOL 或 RPTSCHED来控制输出)2. ECLIPSE输出文件都有什么格式?格式化输出: 可读文件,文件大。 (用关键字 FMTOUT来控制) 非格式化输出:不可读文件,文件小。 多输出文件: 每一时间步一个输出文件。 单文件输出: 所有时间步输出到一个文
3、件。 (用关键字 UNIFOUT来控制)ECLIPSE 缺省输出: 非格式化,多文件输出。3. ECLIPSE数据文件分几部分,各部分定义什么数据类型?ECLIPSE数据类型分八部分,各部分内的关键字除几个个别的外不能混用。RUNSPEC: 定义模型维数以及模型基本类型,包括模型网格维数,最大井数,井组数,流体类型,输出类型控制等。 GRID: 定义模型网格和属性,包括顶部深度,厚度,孔隙度,渗透率,净毛比,一般由前处理软件 Flogrid或 Petrel输出。 EDIT: 编辑孔隙体积,传导率。 PROPS: 流体 PVT及岩石数据,包括油、气体积系数,粘度随压力变化,水的体积系数,粘度;油
4、,气,水地面密度等。岩石数据是相渗曲线和毛管压力。 REGIONS: 分区数据,包括流体分区,岩石分区,储量区,平衡区等。 SOLUTION: 平衡区数据,包括油水界面、油气界面、参考压力、参考深度、水体参数。 SUMMARY: 计算结果输出,包括油田,井组,单井的油、气、水产量,压力输出,网格的压力,饱2和度输出等。 SCHEDULE: 动态数据部分,包括定义井位,射孔,产量,压力,历史拟合,预测等。4. 在 RUNSPEC中定义什么内容?模型大小:DIMENS 模型流体:OIL,GAS,WATER,DISGAS,VAPOIL,根据油藏实际情况选用不同组合。 模型单位:METRIC, FIE
5、LD, LAB 模型维数:WELLDIMS(井,射孔和井组维数) ,VFPIDIMS(注水井 VFP表) ,VFPPDIMES(生产井 VFP表),EQLDIMS(平衡区) ,TABDIMS(流体及岩石维数) 模型输出格式控制:UNIFOUT,FMTOUT,UNIFIN,FMTIN特殊功能:API,BRINE,DUALPERM,DUALPORO,MISCIBLE,TRACES标题: TITLES 模拟开始时间:START5在 GRID部分需要定义什么参数?模型构造: COORD,ZCORN 或 DX,DY,DZ,TOPS; 模型孔隙度:PORO; 模型渗透率:PERMX,PERMY,PERMZ
6、; 模型净毛比或净厚度:NTG,DZNET; 断层: FAULTS。 可以采用正交网格或角点网格,还可以采用 PEBI网格。网格部分数据最好由前处理软件提供。6. 如何修改 GRID部分参数?COPY: 拷贝 ADD: 增加 MULTIPLY: 乘7. 如何修改断层传导率?用 MULTFLT修改每条断层的传导率。8. 如何定义流体?流体数据来自于试验报告,描述油、气体积系数,粘度随压力的变化,油、气、水地面密度等。 油、气、水、溶解气和挥发油: PVTO(或 PVCO), PVTG,PVTW; 油、气、水和溶解气: PVTO(或 PVCO),PVDG(或 PVZG),PVTW; 油、水: PV
7、DO(或 PVCDO),PVTW; 气、水: PVDG,PVTW。39. 如何定义不同类型流体?提供多个流体定义表,在 REGIONS部分定义 PVT分区,每一个网格将对应一套 PVT属性表。10. 如何定义岩石?用关键字 ROCK定义。需要定义在参考压力下岩石的压缩系数。如果岩石随地层压力变化,孔隙度和渗透率会发生变化,可以用 ROCKTAB来定义岩石压缩表。这样当断层压力降低造成上浮压力上升,ECLIPSE会通过 ROCKTAB来模拟孔隙度和渗透率的变化。11. 如何定义相对渗透率和毛管压力曲线? 来自于试验报告,两相相对渗透率和毛管压力曲线。 油水: SWOF 或 SWFN,SOF2;
8、油气: SGOF 或 SGFN,SOF3。12. 如何定义分区?实际油藏非常复杂,在不同部位可能流体属性不同,同时由于非均质造成不同的岩石属性,需要应用不同的相对渗透率和毛管压力曲线,或者油藏有多个油水界面,这是都需要对模型进行分区。可以进行以下分区: 储量分区: FIPNUM 不同分区可以输出不同的储量; PVT分区: PVTNUM 对应不同的流体类型; 岩石分区: SATNUM 对应不同的相对渗透率和毛管压力曲线; 平衡分区: EQUNUM 对应不同的油气或油水界面。 通常需要用前处理软件 FLOGRID或 PETREL来帮助进行分区,比如根据孔隙度的不同范围进行岩石分区,使孔隙度不同范围
9、的网格应用不同的相渗曲线。13. 在 SOLUTION部分如何进行模型初始化?需要提供参考深度、参考深度对应的原始压力、油水界面、油水界面处的毛管压力、气水界面、气水界面处的毛管压力和饱和压力或溶解油气比表。 如果油水界面为自由水面,则油水界面处的毛管压力为 0。 ECLIPSE将计算原始条件下压力分布,油,气,水饱和度分别。计算原始油,气,水储量及各分区储量。14. 如何开始历史拟合?需要定义井口位置(WELSPECS);井轨迹和井射孔位置(COMPDAT);井的类型(WELSPECS),包括生产井和注水,注气井;井的实际产油量、产气量、产水量、井口压力和井底流压(WCONHIST)。 这部
10、分数据量非常大,需要用 ECLIPSE前处理软件 SCHEDULE来完成。415. 历史拟合阶段控制什么生产?可以定油量、气量、水量、液量和定油藏产量,井口压力和井底压力此时只用于与计算结果拟合,不会对生产进行控制。16. 如何拟合井口压力?首先需要为每口井提供其垂直管流表(VFP),此表由前处理软件 VFPi生成,描述不同井口压力在不同产量、含水、油气比条件下对应的井底压力。然后在 WCONHIST里为每口井设置其所用的 VFP表,这样 ECLIPSE就会在每一时间步输出井口压力,对比此压力与实际压力来完成历史拟合。压力拟合可以有两种方式,一种是拟合在一定油、气、水产量和井口压力时不同井筒深
11、度对应的压力,也就是压力剖面。另一种是拟合对应不同产量和井口压力时的不同井底压力。 对于第一种方式,输入 THP, tubing pressure是对应不同深度的压力剖面。 对于第二种方式,输入 THP, tubing pressure就是 BHP。17. 如何定义注水和注气井?对于注水和注气井,历史拟合和预测阶段用不同的关键字。在历史拟合阶段用 WCONINJH,此时ECLIPSE 控制地面注入量,而 BHP 和 THP 用于拟合,不用于控制。在预测阶段用 WCONINJE,此时注入井可以控制地面注入量,油藏注入量,BHP 或 THP 控制以及井组控制。18. 如何定义产量预测? 用 WCO
12、NPROD 来定义产量预测。可以定地面油、气、水、液量或油藏产量生产,也可定BHP、THP 或井组产量生产。要注意此时所有的数据项都起控制作用,比如井在定油量生产时,其产水量、产气量受 WCONPROD 上限控制,BHP 和 THP 受 WCONPROD 下限控制。19. 如何控制经济极限? 可以用 CECON、WECON 或GECON 控制产层、井和井组最小产油量、产气量;最大含水、油气比。在违背经济极限时可以采用一系列作业,比如关井、封层(WPLUG)、修井(WORKTHP)、气举(WLIFT)、换油管(WLIFT)、钻新井(QDRILL)等。20. 如何定义井组生产? 用 GCONPRO
13、D定义井组生产。井组生产可以定地面油、气、水、液量或油藏产量。该井组可以受上一级井组产量控制或独立于上一级井组。21. 如何定义井组注入? 用 GCONINJE定义井组注入。井组可以定地面注入量、油藏注入量、地面注采比和油藏亏空。22. 如何控制 ECLIPSE计算时间步? 通常情况下 ECLIPSE缺省计算参数设置是最好的,所以一般情况下不需要修改 ECLIPSE计算参数设置。在模型出现收敛性问题时可以用 TUNING来修改 ECLIPSE计算参数。比如最大时间步、最小时间步、非线形收敛误差、物质平衡误差和最大时间截断误 差等。5 建议需要时只修改最大时间步,不要修改其他参数。23. 什么是
14、重启计算? 历史拟合结束后需要进行产量预测,在进行产量预测计算时,不需要再从历史拟合开始时进行计算,可以直接从历史拟合结束的时间接着往下算。这种应用上一次计算的输出作为下一次计算的初始输入计算就叫重启计算。 要进行重启计算,首先要定义重启时间步的输出。可以用 RPTRST来定义输出每时间步,每月,每年或每隔几月几年重启时间步文件。如果采用多文件格式输出,则文件后缀为:.X0000, .X0001等,如果是单文件输出,则输出文件为.UNRST. 重启文件记录了每时间步模型压力分布,饱和度分布,溶解油气比分布,同时也记录所有井的井位,射孔位置,产量控制。不过重启文件没有记录垂直管流表(VFP 表)
15、 ,所以在应用垂直管流表时要记住重启时需加上垂直管流表。24. 如何进行重启计算? ECLIPSE有两种重启计算方法,快速重启和完全重启。 快速重启不需要重新处理 RUNSPEC,GRID,EDIT,PROPS和 REGIONS部分,如果在历史拟合计算时设了 SAVE关键字,这些部分将保存在输出的 SAVE文件中,这样在重启计算时不用再计算传导率。 完全重启需要重新处理 RUNSPEC,GRID,EDIT,PROPS和 REGIONS部分,要重新计算传导率。完全重启步骤: 在历史拟合部分用 RPTRST要求输出重启文件。 在 PRT文件中检查重启时间对应的重启文件步。 将历史拟合文件拷贝为重启
16、文件。 删掉 SOLUTION部分中的 EQUIL和水体部分,用 RESTART关键字设重启。 在 SCHEDULE部分用 SKIPREST或删掉所有历史拟合时间步。 如果有 VFP表,要保留 VFP表。 增加新时间步进行预测计算。快速重启步骤: 在历史拟合部分用 SAVE和 RPTRST要求输出 SAVE文件和重启文件。 将历史拟合文件拷贝为重启文件。 删掉所有 SUMMARY以前部分。 用 LOAD关键字装载 SAVE文件。 用 RESTART设重启时间。 在 SCHEDULE部分用 SKIPREST或删掉所有历史拟合时间步。 如果有 VFP表,要保留 VFP表。6 增加新时间步进行预测计
17、算。快速重启和完全重启的优劣: 快速重启不用重新计算传导率。对于大模型,可以节省时间。 完全重启可以修改历史拟合部分参数。25. 如何将 CMG模型转换为 ECLIPSE模型?(1)CMG 软件使用“*$“作为注释行,ECLIPSE 软件使用“”作为注释行。 ECLIPSE软件许多关键字结束后需要用“/”来表示结束,而 CMG不需要。 在采用公制时,CMG 压力单位 KPA,而 ECLIPSE为 BARS. (2)网格部分的转换: CMG Gridbuilder可以输出 ECLIPSE格式网格数据,所以只需要将 CMG模型输入到 CMG Gridbuilder,然后再输出 ECLIPSE格式就
18、可以了。 属性参数只需要改关键字。 CMG用 MOD修改参数,而 ECLIPSE用 ADD,MULTIPLY。 (3)PVT 数据与相渗曲线转换: CMG与 VIP类似,主要是关键字的转换,数据格式不一致。不过在 EXCEL的帮助下都可以得到转换。ECLIPSE定义 PVT数据时将油气、油水分开用不同的关键字定义,而 CMG同时定义。所以在转换时需要分开数据。相渗曲线没有区别。但注意毛管压力单位,ECLISPE 米制用 BARS,CMG 用 KPA。 对于组分模型,区别也主要是关键字不同。EOS 是一样的。 (3)分区数据: 只需要转换关键字。 (4)平衡数据: 数据量很小,转换简单。要注意压
19、力单位,用公制的话,ECLIPSE 参考面压力单位是 BARS,而 CMG是 KPA。 水体定义不同,需要做拟合。 (5) 生产数据: 这是最复杂的部分,建议不要转换,有原始数据的话,重新从头做。VFP 表也不同,最好重做。 如果必须转换,需要编程。ECLIPSE 对于历史拟合和产量预测用不同的控制关键字,历史拟合时用 WCONHIST,此时定义的 BHP和 THP不起控制作用,只用于拟合,在预测时用 WCONPROD,此时 BHP和 THP起控制作用。CMG 对于历史拟合和产量预测没有区别。CMG 会将所有井排序,在定义产量控制时,直接用井的序列号。ECLIPSE 用井名或井组定义。 (6)
20、检查7 检查三维模型,是否发生坐标倒转,属性与网格是否对应。 检查储量是否拟合。 检查井位是否正确。 产量是否拟合。26. 如何将 VIP模型转换为 ECLIPSE模型?VIP软件提供转换器,可以自动将 ECLIPSE模型转换为 VIP模型。但 ECLIPSE并不提供此类转换器,需要手工转换。下面介绍详细步骤: (1)VIP 软件使用“C“作为注释行,ECLIPSE 软件使用“”作为注释行。 ECLIPSE软件许多关键字结束后需要用“/”来表示结束,而 VIP不需要。 在采用公制时,VIP 压力单位 KPA,而 ECLIPSE为 BARS。 VIP模型初始化和生产分两部分,而 ECLIPSE在
21、一个文件。 ECLIPSE模型数据体分为 RUNSPEC(维数定义) ,GRID( 网格定义) ,EDIT(修改) ,PROPS(流体及岩石参数) ,SOLUTION(油水、气水界面和初始压力),REGIONS(分区) ,SUMMARY(输出参数) ,SCHEDULE(生产数据) 。 VIP模型数据体先定义输出参数,然后定义油水、气水界面,初始压力,流体及岩石参数和定义网格。在另一个文件里定义生产数据。 ECLIPSE通常缺省计算方法是最好的,一般不用设置控制计算参数,除非在不收敛情况下。VIP 往往要设置计算控制参数。 (2) 网格部分的转换: VIP用一个关键字 COPR来定义角点网格,E
22、CLIPSE 使用两个关键字 COORD和 ZCORN来定义网格。COORD定义从顶到底的坐标,ZCORN 定义每个角点的深度。可以用 ECLIPSE软件中的 GRID来帮助将 VIP网格转化为 ECLIPSE网格。步骤为:“Edt model“Enter the name of the modelor RETURN for menu options dummyModel DUMMY does not exist.Create new model ? (Y/n) YAVAILABLE MAP PROJECTIONS1. Universal Transverse Mercator2. Obliq
23、ue Mercator (Hotine)3. Transverse Mercator84. Mercator5. Lambert Conformal Conic6. Albers Equal Area7. Universal Polar Stereographic8. Azimuthal Equidistant9. None ( normal cartesian coordinates )Choose an option ( H for HELP )or RETURN for 9 Enter units for model area ( H for Help )or RETURN for ME
24、TRES METRESEnter minimum X and Y coordinates for model area 0 0 (you can use your axis value in the model )Enter maximum X and Y coordinates for model area 10000 10000 (you can use your axis value in the model)Define default number of mesh rows & columnsas NROW & NCOL or XINC & YINC ? (R/i) REnter d
25、efault number of mesh rows or RETURN for 50 Enter default number of mesh columns or RETURN for 50 Enter maximum number of model strata ( 1-100 )or RETURN for 1 “Return“=Create a grid=“Edt grid“Enter the name of the grid xxxGrid XXX does not exist. Create new grid ? (Y/n) YEnter name of map ( H for H
26、elp ) 9=“Input grid“Enter depth units for grid or RETURN for FEET FEETEnter origin as MAP coords, or G for grid default or D to digitize.&Enter MAP coords or option (g/D) D0 0Enter a point along the X axis 100 0Does the Y axis point upwards? (y/N) YOkay to continue with these axes ? (Y/n) YInput Ecl
27、ipse EGRID or GRID file, data file or VIP file ? (e/g/D/v) VEnter name of file to be read Enter units of input grid or RETURN for METRES feetIs the grid in map coordinates ? (y/N) yOutput ECLIPSE Grid 也可以用 GRIDSIM来转还网格数据。属性数据将关键字转换后可以直接用,但一定要先检查好坐标方向,确定两模型方向一致,否则将会上下颠倒。VIP数据修改直接在所需修改数据下用 MOD,而 ECLIP
28、SE用 ADD,MULTIPLY来修改参数。 (2)PVT 数据与相渗曲线转换: 主要是关键字的转换,数据格式不一致。不过在 EXCEL的帮助下都可以得到转换。ECLIPSE 定义PVT数据时将油气,油水分开用不同的关键字定义,而 VIP同时定义。所以在转换时需要分开数据。相渗曲线没有区别。但注意毛管压力单位,ECLISPE 米制用 BARS,VIP 用 KPA。 对于组分模型,区别也主要是关键字不同。EOS 是一样的。 (3)分区数据: 只需要转换关键字。不过一定要检查坐标方向,否则是不对的。 (4) 平衡数据: 数据量很小,转换简单。VIP 直接在 IEQUIL里定义饱和压力,ECLIPS
29、E 用 PDVD或 RSVD来定义随深度变化。要注意压力单位,用公制的话,ECLIPSE 参考面压力单位是 BARS,而 VIP是 KPA. 水体定义不同,需要做拟合。 (5) 生产数据:10 这是最复杂的部分,建议不要转换,有原始数据的话,重新从头做。VFP 表也不同,最好重做。 如果必须转换,需要编程。ECLIPSE 对于历史拟合和产量预测用不同的控制关键字,历史拟合时用 WCONHIST,此时定义的 BHP和 THP不起控制作用,只用于拟合,在预测时用 WCONPROD,此时 BHP和 THP起控制作用。VIP 对于历史拟合和产量预测没有区别。 VIP会将所有井排序,在定义产量控制时,直
30、接用井的序列号。ECLIPSE 用井名或井组定义。 (6)检查 检查三维模型,是否发生坐标倒转,属性与网格是否对应。 检查储量是否拟合。 检查井位是否正确。 产量是否拟合。如 何 加 快 数 模 计 算 以 及 如 何 解 决 数 模 计 算 的 收 敛 性 问 题 ? 本文将以 ECLIPSE软件为例从两方面介绍收敛性问题。第一方面介绍数值模拟计算与收敛有关的一些概念。第二部分介绍如何通过修改模型数据来加速计算,解决收敛性问题。一:数模计算的收敛性: 在了解收敛性之前,应该首先了解几个基本概念:1。报告步: 一个数模作业包括多个报告步,报告步是用户设置要求多长时间输出运行报告,比如可以每个月
31、,每季度或每年输出运行报告,运行报告包括产量报告和动态场(重启)报告。在 ECLIPSE软件中,报告步是通过 DATES和 TSTEP关键字来设置的。2。时间步: 一个报告步包括多个时间步,时间步是软件自动设置(VIP 需要用户设置)即通过多个时间步的计算来达到下一个报告步,以 ECLIPSE为例,假如报告步为一个月,在缺省条件下,ECLISPE 第一个时间步取一天,然后以三倍增加,即第二个时间步取三天,然后取九天,下一个时间步是 17天来达到 30天的报告步,然后会以每 30天的时间步来计算。时间步可以通过 TUNING关键字来修改。3。非线形迭代:一个时间步包括多次非线形迭代。在缺省情况下
32、,ECLIPSE 如果通过 12次的非线形迭代没有收敛,ECLIPSE 将对时间步减小 10倍。比如下一个时间步应该是 30天,如果通过 12次的迭代计算不能达到收敛,ECLIPSE 将把时间步缩短为 3天。下一个时间步将以 1.25倍增长,即 3.75天,4.68天, 。 。 。 。如果在计算过程中经常发生时间步的截断,计算将很慢。4。线形迭代: 一个非线形迭代包括多次线形迭代。线形迭代是解矩阵。 在 ECLIPSE输出报告 PRT文件中可以找到时间步,迭代次数的信息, STEP 10 TIME= 100.00 DAYS ( +10.0 DAYS REPT 5 ITS) (1-FEB-200
33、8) “STEP 10” : 说明这是第 10个时间步。 “TIME= 100.00 DAYS”: 说明现在模拟到第 100天。 “+10.0 DAYS”: 说明这个时间步是 10天。11 “REPT“ : 说明为什么选 10天做为时间步,REPT 是指由于到了下一个报告步。 “5 ITS“: 说明此 10天时间步需要 5次非线形迭代。 ”(1-FEB-2008)“: 现在的模拟时间。 模拟计算的时间取决于时间步的大小,如果模型没有发生时间步的截断而且能保持长的时间步,那表明该模型没有收敛性问题,反之如果经常发生时间步截断,那模型计算将很慢,收敛性差。时间步的大小主要取决于非线形迭代次数。如果
34、模型只用一次非线形迭代计算就可以收敛,那表明模型很容易收敛,如果需要 2到 3次,模型较易收敛,如果需要 4到 9次,那模型不易收敛,大于 10次的化模型可能有问题,如果大于 12次,时间步将截断。在 PRT文件中如果看到以下信息: PROBLEM: AT TIME 200 DAYS (1-FEB-2009): NON-LINEAR EQUATION CONVERGENCE FAILURE ITERATION LIMIT REACHED - TIME STEP CHOPPED FROM 10 STEP 20 TIME= 200.00 DAYS ( +1.0 DAYS CHOP 5 ITS) (
35、1-FEB-2009) 那表明时间步发生了截断。 (注: 另外如果你见到如下信息: WARNING AT TIME 0.0 DAYS (1-MAR-2004): LINEAR EQUATIONS NOT FULLY CONVERGED - RUN MAY GO FASTER IF YOU INCREASE LITMX (=25 - TUNING KEYWORD) 你可以不必管。这只是线形方程不收敛) 除了 REPT,CHOP外,在 RPT文件中还常见以下信息来表明为什么选择现在的时间步: INIT: 表明是初始时间步 TRNC: 为满足时间截断误差 MINS: 最小时间步 MAXS: 最大时间
36、步 HALF: 接近报告步时的时间步取半 DIFF: 时间步截断 CHOP之后的增长 (在 ECLIPSE技术手册的第 125页还会找到更多) 如果模型中有很多 CHOP,DIFF,MINS,那模型有严重的收敛性问题。 二:如何解决数模计算的收敛性问题: 模型不收敛的原因很多,网格参数,属性参数,流体 PVT参数,岩石相渗曲线,毛管压力曲线,相渗曲线端点标定,初始化,井轨迹,垂直管流表都会造成模型不收敛,下面分别介绍各部分如何造成不收敛及如何解决。1。网格部分:网格正交性差和网格尺寸相差太大是导致不收敛的主要原因之一。正交性差会给矩阵求解带来困难,而网格尺寸相差大会导致孔隙体积相差很大,大孔隙
37、体积流到小孔隙体积常会造成不收敛。 解决办法:12 网格正交性差通常是在建角点网格时为描述断层或裂缝的走向而造成的。在此情况下,最好能使边界与主断层或裂缝走向平行,这样一方面网格可以很好地描述断层或裂缝,另一方面正交性也很好。 在平面上最好让网格大小能够较均匀,在没有井的地方网格可以很大,但最好能够从大到小均匀过渡。纵向上有的层厚,有的薄,最好把厚层能再细分。在检查模型时应该每层每层都在三维显示中检查。径向局部网格加密时里面最小的网格不要太小。在 ECLIPSE里用 MINPV关键字可以把小于设定孔隙体积的网格设为死网格,这样通常会有用。2。属性参数: 不合理的插值计算会导致属性分布很差,如果
38、是从地质模型粗化为数模模型,通常问题不大,只是有时候数模人员自己插值时会有问题。 解决办法: 有可能尽量用地质模型的数据,自己插值时可以加一些控制点使属性合理分布。X,Y 方向的渗透率最好相等或级差不大。在井连通网格的 Z方向渗透率不要设为 0,如果想控制垂向流动,可给一个很小的值。3。流体 PVT参数 流体 PVT参数会有两种可能的问题,一是数据不合理导致了负总压缩系数,二是压力或气油比范围给的不够导致模型对 PVT参数进行了外插。 解决办法: 检查 PRT文件中的 WARNING信息,如果在油藏压力范围内有负总压缩系数的警告,应该修改 PVT参数,否则的化会有收敛性问题。如果负总压缩系数是
39、在油藏压力范围之外,可以忽略该警告。此部分的修改主要可以小规模修改油和气的 FVF和 RS。关于 ECLIPSE如何计算负总压缩系数,参照 ECLIPSE技术手册第 936页。 在 ECLIPSE中加 EXTRAPMS关键字可以要求输出如果发生 PVT插值后的警告信息。在提供 PVT表时,压力应该覆盖所有范围,包括注水后的压力上升。RS 值也应该考虑到气在油中的重新溶解。4。岩石相渗曲线和毛管压力曲线。 ECLIPSE不会对输入模型的相渗曲线和毛管压力曲线进行光滑,将会应用每一个输入饱和度和相渗值,所以要保证输入的参数是合理的。通常的问题有: (1): 饱和度和相对渗透率的数据位数过多。 (2
40、): 饱和度值太接近,导致相渗曲线的倾角变化很大。 (3): 饱和度有很小变化但相对渗透率发生了很大变化。 解决办法: (1):饱和度和相对渗透率最多给两位小数就够了。 (2):检查相渗曲线的导数?梢杂 肊 CLIPSE中的 SCAL模块) ,导数要光滑。13 ?3):将临界饱和度和束缚饱和度设为不同的值。5。端点标定: 在应用端点标定时,有时标定完后的相渗曲线倾角很大,标定后的毛管了很大。 解决办法: 在三维显示中检查标定完的 PCW,可以给 PCW一个最大值来控制毛管压力。 输出每个网格标定后的相渗曲线进行检查。6。初始化: 初始化最容易发生的问题是在初始时模型不稳定,流体在初始条件下就会
41、发生流动,这也会导致模型不收敛。造成模型初始不稳定的主要有: (1) 手工赋网格饱和度和压力值。 (2)拟合初始含水饱和度。 解决办法: (1)尽量不要直接为网格赋压力和饱和度值,尽量由模型通过油水界面及参考压力来进行初始化计算。 (2)要想拟合地质提供的初始含水饱和度分布,应该进行毛管压力的端点标定,这样毛管压力会稳住每个网格的水,在初始条件下不会流动。 (3)可以通过让模型在没有任何井的情况下计算十年来检查初始条件下模型是否稳定,如果 10年的计算模型压力和饱和度度没有变化,说明模型初始是稳定的。7。井轨迹: 在进行井处理时井可能以之字型在网格中窜过,有可能发生井的实际窜过方向与模型关键字
42、定义的方向不符,这也会导致不收敛。 解决办法: (1) 在三维显示中检查井轨迹。 (2) 如果井已经关掉,在模拟时不要给零产量,要用关键字把井关掉。 (3) 检查井射孔,井不要射在孤立的网格上。8:垂直管流曲线: 有了垂直管流曲线很容易导致模型不收敛,这有两种可能: (1)曲线有交叉。 (2)曲线发生了外插。 解决办法: (1) 用前处理软件(ECLIPSE 中的 VFPi)检查曲线。 (2)在 ECLIPSE中加 EXTRAPMS关键字可以要求输出如果发生 VFP插值后的警告信息。14 (3)曲线应该覆盖所有井口压力,含水,油气比及产量。 (4)在 ECLIPSE用 WVFPEXP。9:其他
43、解决办法: 如果模型数据没有问题,可以调整模拟器的收敛计算参数,对于 ECLIPSE,可以做以下调整: (1)调整 TUNING中的最大时间步。如果模型每计算到 30天就会截断时间步,可以将最大时间步调整为 20天,这样计算会快很多。 (2)调整 TUNING中的最大线形迭代次数到 70次。(3)降低 TUNING中的线形收敛误差标准(4)对于组分模型用 FREEZEPC和 DPCDT.15 因为师兄给了我一些地质图件和比如原始地质储量的某些参数,让我完成地质模型,但我不知道从何处着手。你可不可以详细地介绍一下有效模型的建立步骤(就以没有任何可引用的文件为例,也就是说参数都要自己输入) ,我还
44、投了邮件在你邮箱里,请多关照!谢谢! 解答: 你的问题在国内比较普遍。问题主要在于各部门之间缺乏合作。很多油田没有项目组的概念,勘探部门做完的项目把结果交给开发部门任务就算完成,而且结果经常以图形和报告形式,而不是开发部门更需要的数据文件,这样导致开发部门无法有效共享勘探部门的研究。这里面有软件的一体化问题,最重要的是管理问题。 现在说你的问题。我不知道你拿到的是井点参数还是等值图。假如你拿到的是井点参数,而且你用的软件是 Flogrid,以下是步骤:1 。 进入 Flogrid的 TOOLS-2D MAPPING CARVAS.2 。 在 2D MAPPING CARVAS窗口,进入 CRE
45、ATE-SCATTER SET。3 。 选择属性。对于数值模拟,你应该完成每个储层孔隙度,渗透率,厚度,净厚度,深度等值图。4 。设置你图形的范围,最大最小值。5 。 对应每一井点,输入其数值。6 。 进入 CREATE MESH MAP7 。选择属性,选择 SCATTER SET, 生成 MESH MAP.8. 将所有图件全部完成后, 参照 FLOGRID手册中的练习十: Generating structure and properties from 2D mapping.FloGrid 网格生成器 FloGid 是一个一体化的产品,它支持三维油藏描述标准格式 RESCUE(POSC)输入
46、。它也可以直接读入油藏描述软件生成的二维地质图,在其基础上自动(或手工)形成油藏模拟网格系统。 生成可由 Eclipse直接使用的各种油藏模拟网格系统(正交、径向、角点、PEBI) 支持三维 RESCUE(POSC标准)格式输入 直接从大多数绘图软件包和油藏描述系统读入数据 支持多相 upscaling(物性粗化) 支持三维非规则网格的局部加密 在处理断层交叉面时, 可采用分段坐标系 采用多点流动方式来反映网格间流动的非正交效应 在建立油藏模拟网格系统时能自动识别和考虑地质特 全三维数据(井、各种图、地质网格、断层和模拟网格)可视化16 提供了一组把地质模型或随机模型的细网格粗化成油藏模拟网格
47、的工具Eclipse Office数模一体化管理模块EclipseOffice 可以打开和管理 Eclipse数模家族的任意软件,允许你在数模运行中随时查看计算结果;可以编辑和评价数模计算结果,并且可生成结果报告;可以快捷的建一个数据研究模型并进行计算。EclipseOffice 是一个一整化的桌面环境,.提供了五个特色模块,给用户控制管理 数模流程提供了极大的方便。 项目管理在 Office环境下管理特定的模拟研究项目 数据管理以数据面板的形式列入所有 ECLIPSE的标字号,并在特定的油田条件下,列出 有关标字号的含义、数值和单位 运行管理启动及监测模拟运行 结果显示显示曲线和二维、三维结果 报告输出S