Eclipse数值模拟软件问答

更新时间:2023-08-17 18:22:01 阅读量: 资格考试认证 文档下载

说明:文章内容仅供预览,部分内容可能不全。下载后的文档,内容与下面显示的完全一致。下载之前请确认下面内容是否您想要的,是否完整无缺。

数值模拟应用

Eclipse数值模拟软件问答(初级)

1. ECLIPSE输出结果文件是哪些?

.GRID或.FGRID: 网格文件

.EGRID: 网格文件,与GRID格式不同,文件要小的多。(用关键字GRIDFILE来控制输出类型) .INIT或.FINIT: 属性文件。(用关键字INIT来控制输出)

.PRT: 报告输出。 文件很大,模型处理及计算结果详细报告。(RPTGRID,RPTPROP,RPTSOL,RPTSCHED控制输出)

.LOG: 后台作业时的输出报告,文件比PRT要小很多。可用于错误检查。

.DBG: Debug文件,一般不用。可用于检查ECLIPSE如何处理输入参数。

.SAVE: 用于快速重启。(用关键字SAVE来控制输出)

.RFT:RFT计算结果。(用关键字WRFTPLT来控制输出)

.FLUX: 流动边界。 (用关键字DUMPFLUX来控制输出)

.Snnnn或.UNSMRY: 图形文件输出(在SUMMARY部分定义)

.Xnnnn或.UNRST: 重启文件输出(用RPTRST,RPTSOL或RPTSCHED来控制输出)

2. ECLIPSE输出文件都有什么格式?

格式化输出: 可读文件,文件大。(用关键字FMTOUT来控制)

非格式化输出:不可读文件,文件小。

多输出文件: 每一时间步一个输出文件。

单文件输出: 所有时间步输出到一个文件。(用关键字UNIFOUT来控制)

ECLIPSE缺省输出: 非格式化,多文件输出。

3. ECLIPSE数据文件分几部分,各部分定义什么数据类型?

ECLIPSE数据类型分八部分,各部分内的关键字除几个个别的外不能混用。

RUNSPEC: 定义模型维数以及模型基本类型,包括模型网格维数,最大井数,井组数,流体类型,输出类型控制等。

 GRID: 定义模型网格和属性,包括顶部深度,厚度,孔隙度,渗透率,净毛比,一般由前处理软件Flogrid或Petrel输出。

 EDIT: 编辑孔隙体积,传导率。

 PROPS: 流体PVT及岩石数据,包括油、气体积系数,粘度随压力变化,水的体积系数,粘度;油,气,水地面密度等。岩石数据是相渗曲线和毛管压力。

 REGIONS: 分区数据,包括流体分区,岩石分区,储量区,平衡区等。

 SOLUTION: 平衡区数据,包括油水界面、油气界面、参考压力、参考深度、水体参数。

 SUMMARY: 计算结果输出,包括油田,井组,单井的油、气、水产量,压力输出,网格的压力,饱

数值模拟应用

和度输出等。

 SCHEDULE: 动态数据部分,包括定义井位,射孔,产量,压力,历史拟合,预测等。

4. 在RUNSPEC中定义什么内容?

模型大小:DIMENS

模型流体:OIL,GAS,WATER,DISGAS,VAPOIL,根据油藏实际情况选用不同组合。

模型单位:METRIC, FIELD, LAB

模型维数:WELLDIMS(井,射孔和井组维数),VFPIDIMS(注水井VFP表),VFPPDIMES(生产井VFP表),EQLDIMS(平衡区),TABDIMS(流体及岩石维数)

模型输出格式控制:UNIFOUT,FMTOUT,UNIFIN,FMTIN

特殊功能:API,BRINE,DUALPERM,DUALPORO,MISCIBLE,TRACES

标题: TITLES

模拟开始时间:START

5.在GRID部分需要定义什么参数?

模型构造: COORD,ZCORN或DX,DY,DZ,TOPS;

模型孔隙度:PORO;

模型渗透率:PERMX,PERMY,PERMZ;

模型净毛比或净厚度:NTG,DZNET;

断层: FAULTS。

可以采用正交网格或角点网格,还可以采用PEBI网格。网格部分数据最好由前处理软件提供。

6. 如何修改GRID部分参数?

COPY: 拷贝

 ADD: 增加

 MULTIPLY: 乘

7. 如何修改断层传导率?

用MULTFLT修改每条断层的传导率。

8. 如何定义流体?

流体数据来自于试验报告,描述油、气体积系数,粘度随压力的变化,油、气、水地面密度等。  油、气、水、溶解气和挥发油: PVTO(或PVCO), PVTG,PVTW;

 油、气、水和溶解气: PVTO(或PVCO),PVDG(或PVZG),PVTW;

 油、水: PVDO(或PVCDO),PVTW;

 气、水: PVDG,PVTW。

数值模拟应用

9. 如何定义不同类型流体?

提供多个流体定义表,在REGIONS部分定义PVT分区,每一个网格将对应一套PVT属性表。

10. 如何定义岩石?

用关键字ROCK定义。需要定义在参考压力下岩石的压缩系数。如果岩石随地层压力变化,孔隙度和渗透率会发生变化,可以用ROCKTAB来定义岩石压缩表。这样当断层压力降低造成上浮压力上升,ECLIPSE会通过ROCKTAB来模拟孔隙度和渗透率的变化。

11. 如何定义相对渗透率和毛管压力曲线?

 来自于试验报告,两相相对渗透率和毛管压力曲线。

 油水: SWOF或SWFN,SOF2;

 油气: SGOF或SGFN,SOF3。

12. 如何定义分区?

实际油藏非常复杂,在不同部位可能流体属性不同,同时由于非均质造成不同的岩石属性,需要应用不同的相对渗透率和毛管压力曲线,或者油藏有多个油水界面,这是都需要对模型进行分区。可以进行以下分区:

 储量分区: FIPNUM 不同分区可以输出不同的储量;

 PVT分区: PVTNUM 对应不同的流体类型;

 岩石分区: SATNUM 对应不同的相对渗透率和毛管压力曲线;

 平衡分区: EQUNUM 对应不同的油气或油水界面。

 通常需要用前处理软件FLOGRID或PETREL来帮助进行分区,比如根据孔隙度的不同范围进行岩石分区,使孔隙度不同范围的网格应用不同的相渗曲线。

13. 在SOLUTION部分如何进行模型初始化?

需要提供参考深度、参考深度对应的原始压力、油水界面、油水界面处的毛管压力、气水界面、气水界面处的毛管压力和饱和压力或溶解油气比表。

如果油水界面为自由水面,则油水界面处的毛管压力为0。

 ECLIPSE将计算原始条件下压力分布,油,气,水饱和度分别。计算原始油,气,水储量及各分区储量。

14. 如何开始历史拟合?

需要定义井口位置(WELSPECS);

井轨迹和井射孔位置(COMPDAT);

井的类型(WELSPECS),包括生产井和注水,注气井;

井的实际产油量、产气量、产水量、井口压力和井底流压(WCONHIST)。

 这部分数据量非常大,需要用ECLIPSE前处理软件SCHEDULE来完成。

数值模拟应用

15. 历史拟合阶段控制什么生产?

可以定油量、气量、水量、液量和定油藏产量,井口压力和井底压力此时只用于与计算结果拟合,不会对生产进行控制。

16. 如何拟合井口压力?

首先需要为每口井提供其垂直管流表(VFP),此表由前处理软件VFPi生成,描述不同井口压力在不同产量、含水、油气比条件下对应的井底压力。然后在WCONHIST里为每口井设置其所用的VFP表,这样ECLIPSE就会在每一时间步输出井口压力,对比此压力与实际压力来完成历史拟合。

压力拟合可以有两种方式,一种是拟合在一定油、气、水产量和井口压力时不同井筒深度对应的压力,也就是压力剖面。另一种是拟合对应不同产量和井口压力时的不同井底压力。

对于第一种方式,输入THP, tubing pressure是对应不同深度的压力剖面。

对于第二种方式,输入THP, tubing pressure就是BHP。

17. 如何定义注水和注气井?

对于注水和注气井,历史拟合和预测阶段用不同的关键字。在历史拟合阶段用WCONINJH,此时ECLIPSE控制地面注入量,而BHP和THP用于拟合,不用于控制。在预测阶段用WCONINJE,此时注入井可以控制地面注入量,油藏注入量,BHP或THP控制以及井组控制。

18. 如何定义产量预测?

用WCONPROD来定义产量预测。可以定地面油、气、水、液量或油藏产量生产,也可定BHP、THP或井组产量生产。要注意此时所有的数据项都起控制作用,比如井在定油量生产时,其产水量、产气量受WCONPROD上限控制,BHP和THP受WCONPROD下限控制。

19. 如何控制经济极限?

可以用CECON、WECON或\GECON控制产层、井和井组最小产油量、产气量;最大含水、油气比。在违背经济极限时可以采用一系列作业,比如关井、封层(WPLUG)、修井(WORKTHP)、气举(WLIFT)、换油管(WLIFT)、钻新井(QDRILL)等。

20. 如何定义井组生产?

用GCONPROD定义井组生产。井组生产可以定地面油、气、水、液量或油藏产量。该井组可以受上一级井组产量控制或独立于上一级井组。

21. 如何定义井组注入?

用GCONINJE定义井组注入。井组可以定地面注入量、油藏注入量、地面注采比和油藏亏空。

22. 如何控制ECLIPSE计算时间步?

 通常情况下ECLIPSE缺省计算参数设置是最好的,所以一般情况下不需要修改ECLIPSE计算参数设置。在模型出现收敛性问题时可以用TUNING来修改ECLIPSE计算参数。比如最大时间步、最小时间步、非线形收敛误差、物质平衡误差和最大时间截断误差等。

数值模拟应用

建议需要时只修改最大时间步,不要修改其他参数。

23. 什么是重启计算?

 历史拟合结束后需要进行产量预测,在进行产量预测计算时,不需要再从历史拟合开始时进行计算,可以直接从历史拟合结束的时间接着往下算。这种应用上一次计算的输出作为下一次计算的初始输入计算就叫重启计算。

 要进行重启计算,首先要定义重启时间步的输出。可以用RPTRST来定义输出每时间步,每月,每年或每隔几月几年重启时间步文件。如果采用多文件格式输出,则文件后缀为:.X0000, .X0001等,如果是单文件输出,则输出文件为.UNRST.

 重启文件记录了每时间步模型压力分布,饱和度分布,溶解油气比分布,同时也记录所有井的井位,射孔位置,产量控制。不过重启文件没有记录垂直管流表(VFP表),所以在应用垂直管流表时要记住重启时需加上垂直管流表。

24. 如何进行重启计算?

 ECLIPSE有两种重启计算方法,快速重启和完全重启。

 快速重启不需要重新处理RUNSPEC,GRID,EDIT,PROPS和REGIONS部分,如果在历史拟合计算时设了SAVE关键字,这些部分将保存在输出的SAVE文件中,这样在重启计算时不用再计算传导率。

 完全重启需要重新处理RUNSPEC,GRID,EDIT,PROPS和REGIONS部分,要重新计算传导率。 完全重启步骤:

在历史拟合部分用RPTRST要求输出重启文件。

在PRT文件中检查重启时间对应的重启文件步。

将历史拟合文件拷贝为重启文件。

 删掉SOLUTION部分中的EQUIL和水体部分,用RESTART关键字设重启。

在SCHEDULE部分用SKIPREST或删掉所有历史拟合时间步。

如果有VFP表,要保留VFP表。

增加新时间步进行预测计算。

快速重启步骤:

在历史拟合部分用SAVE和RPTRST要求输出SAVE文件和重启文件。

将历史拟合文件拷贝为重启文件。

删掉所有SUMMARY以前部分。

用LOAD关键字装载SAVE文件。

用RESTART设重启时间。

在SCHEDULE部分用SKIPREST或删掉所有历史拟合时间步。

如果有VFP表,要保留VFP表。

数值模拟应用

增加新时间步进行预测计算。

快速重启和完全重启的优劣:

 快速重启不用重新计算传导率。对于大模型,可以节省时间。

完全重启可以修改历史拟合部分参数。

25. 如何将CMG模型转换为ECLIPSE模型?

(1)CMG软件使用"**$"作为注释行,ECLIPSE软件使用“--”作为注释行。

 ECLIPSE软件许多关键字结束后需要用“/”来表示结束,而CMG不需要。

 在采用公制时,CMG压力单位KPA,而ECLIPSE为BARS.

(2)网格部分的转换:

 CMG Gridbuilder可以输出ECLIPSE格式网格数据,所以只需要将CMG模型输入到CMG Gridbuilder,然后再输出ECLIPSE格式就可以了。

 属性参数只需要改关键字。

 CMG用MOD修改参数,而ECLIPSE用ADD,MULTIPLY。

(3)PVT数据与相渗曲线转换:

 CMG与VIP类似,主要是关键字的转换,数据格式不一致。不过在EXCEL的帮助下都可以得到转换。ECLIPSE定义PVT数据时将油气、油水分开用不同的关键字定义,而CMG同时定义。所以在转换时需要分开数据。相渗曲线没有区别。但注意毛管压力单位,ECLISPE米制用BARS,CMG用KPA。

 对于组分模型,区别也主要是关键字不同。EOS是一样的。

(3)分区数据:

 只需要转换关键字。

(4)平衡数据:

数据量很小,转换简单。要注意压力单位,用公制的话,ECLIPSE参考面压力单位是BARS,而CMG是KPA。

 水体定义不同,需要做拟合。

(5) 生产数据:

这是最复杂的部分,建议不要转换,有原始数据的话,重新从头做。VFP表也不同,最好重做。  如果必须转换,需要编程。

ECLIPSE对于历史拟合和产量预测用不同的控制关键字,历史拟合时用WCONHIST,此时定义的BHP和THP不起控制作用,只用于拟合,在预测时用WCONPROD,此时BHP和THP起控制作用。CMG对于历史拟合和产量预测没有区别。

CMG会将所有井排序,在定义产量控制时,直接用井的序列号。ECLIPSE用井名或井组定义。 (6)检查

数值模拟应用

检查三维模型,是否发生坐标倒转,属性与网格是否对应。

检查储量是否拟合。

检查井位是否正确。

产量是否拟合。

26. 如何将VIP模型转换为ECLIPSE模型?

VIP软件提供转换器,可以自动将ECLIPSE模型转换为VIP模型。但ECLIPSE并不提供此类转换器,需要手工转换。下面介绍详细步骤:

(1)VIP软件使用"C"作为注释行,ECLIPSE软件使用“--”作为注释行。

 ECLIPSE软件许多关键字结束后需要用“/”来表示结束,而VIP不需要。

 在采用公制时,VIP压力单位KPA,而ECLIPSE为BARS。

 VIP模型初始化和生产分两部分,而ECLIPSE在一个文件。

 ECLIPSE模型数据体分为RUNSPEC(维数定义),GRID(网格定义),EDIT(修改),PROPS(流体及岩石参数),SOLUTION(油水、气水界面和初始压力),REGIONS(分区),SUMMARY(输出参数),SCHEDULE(生产数据)。

 VIP模型数据体先定义输出参数,然后定义油水、气水界面,初始压力,流体及岩石参数和定义网格。在另一个文件里定义生产数据。

 ECLIPSE通常缺省计算方法是最好的,一般不用设置控制计算参数,除非在不收敛情况下。VIP往往要设置计算控制参数。

(2) 网格部分的转换:

 VIP用一个关键字COPR来定义角点网格,ECLIPSE使用两个关键字COORD和ZCORN来定义网格。COORD定义从顶到底的坐标,ZCORN定义每个角点的深度。可以用ECLIPSE软件中的GRID来帮助将VIP网格转化为ECLIPSE网格。步骤为:

"Edt model"

Enter the name of the model

or RETURN for menu options 'dummy'

Model DUMMY does not exist.

Create new model ? (Y/n) 'Y'

AVAILABLE MAP PROJECTIONS



1. Universal Transverse Mercator

2. Oblique Mercator (Hotine)

3. Transverse Mercator

数值模拟应用

4. Mercator

5. Lambert Conformal Conic

6. Albers Equal Area

7. Universal Polar Stereographic

8. Azimuthal Equidistant

9. 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 METRES 'METRES'

Enter 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 & columns

as NROW & NCOL or XINC & YINC ? (R/i) 'R'



Enter default 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 'xxx'

Grid XXX does not exist. Create new grid ? (Y/n) 'Y'

Enter name of map ( H for Help ) ' '

数值模拟应用

====================================================

====================================================

"Input grid"

Enter depth units for grid or RETURN for FEET 'FEET'

Enter origin as MAP coords, or G for grid default or D to digitize...&

Enter MAP coords or option (g/D) 'D'

'0 0'

Enter a point along the X axis '100 0'

Does the Y axis point upwards? (y/N) 'Y'

Okay to continue with these axes ? (Y/n) 'Y'



Input Eclipse EGRID or GRID file, data file or VIP file ? (e/g/D/v) V'

Enter name of file to be read '

Enter units of input grid or RETURN for METRES 'feet'

Is the grid in map coordinates ? (y/N) 'y'

Output ECLIPSE Grid

也可以用GRIDSIM来转还网格数据。

属性数据将关键字转换后可以直接用,但一定要先检查好坐标方向,确定两模型方向一致,否则将会上下颠倒。

VIP数据修改直接在所需修改数据下用MOD,而ECLIPSE用ADD,MULTIPLY来修改参数。

(2)PVT数据与相渗曲线转换:

 主要是关键字的转换,数据格式不一致。不过在EXCEL的帮助下都可以得到转换。ECLIPSE定义PVT数据时将油气,油水分开用不同的关键字定义,而VIP同时定义。所以在转换时需要分开数据。相渗曲线没有区别。但注意毛管压力单位,ECLISPE米制用BARS,VIP用KPA。

对于组分模型,区别也主要是关键字不同。EOS是一样的。

(3)分区数据:

 只需要转换关键字。不过一定要检查坐标方向,否则是不对的。

(4) 平衡数据:

 数据量很小,转换简单。VIP直接在IEQUIL里定义饱和压力,ECLIPSE用PDVD或RSVD来定义随深度变化。要注意压力单位,用公制的话,ECLIPSE参考面压力单位是BARS,而VIP是KPA.

 水体定义不同,需要做拟合。

(5) 生产数据:

数值模拟应用

 这是最复杂的部分,建议不要转换,有原始数据的话,重新从头做。VFP表也不同,最好重做。  如果必须转换,需要编程。

ECLIPSE对于历史拟合和产量预测用不同的控制关键字,历史拟合时用WCONHIST,此时定义的BHP和THP不起控制作用,只用于拟合,在预测时用WCONPROD,此时BHP和THP起控制作用。VIP对于历史拟合和产量预测没有区别。

 VIP会将所有井排序,在定义产量控制时,直接用井的序列号。ECLIPSE用井名或井组定义。  (6)检查

检查三维模型,是否发生坐标倒转,属性与网格是否对应。

检查储量是否拟合。

检查井位是否正确。

产量是否拟合。

如何加快数模计算以及如何解决数模计算的收敛性问题?本文将以ECLIPSE软件为例从两方面介绍收敛性问题。第一方面介绍数值模拟计算与收敛有关的一些概念。第二部分介绍如何通过修改模型数据来加速计算,解决收敛性问题。

一:数模计算的收敛性: 在了解收敛性之前,应该首先了解几个基本概念:

1。报告步: 一个数模作业包括多个报告步,报告步是用户设置要求多长时间输出运行报告,比如可以每个月,每季度或每年输出运行报告,运行报告包括产量报告和动态场(重启)报告。在ECLIPSE软件中,报告步是通过DATES和TSTEP关键字来设置的。

2。时间步: 一个报告步包括多个时间步,时间步是软件自动设置(VIP需要用户设置)即通过多个时间步的计算来达到下一个报告步,以ECLIPSE为例,假如报告步为一个月,在缺省条件下,ECLISPE第一个时间步取一天,然后以三倍增加,即第二个时间步取三天,然后取九天,下一个时间步是17天来达到30天的报告步,然后会以每30天的时间步来计算。时间步可以通过TUNING关键字来修改。

3。非线形迭代:一个时间步包括多次非线形迭代。在缺省情况下,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-2008)

 “STEP 10” : 说明这是第10个时间步。

 “TIME= 100.00 DAYS”: 说明现在模拟到第100天。

 “+10.0 DAYS”: 说明这个时间步是10天。

数值模拟应用

 “REPT" : 说明为什么选10天做为时间步,REPT是指由于到了下一个报告步。

 “5 ITS": 说明此10天时间步需要5次非线形迭代。

 ”(1-FEB-2008)“: 现在的模拟时间。

 模拟计算的时间取决于时间步的大小,如果模型没有发生时间步的截断而且能保持长的时间步,那表明该模型没有收敛性问题,反之如果经常发生时间步截断,那模型计算将很慢,收敛性差。时间步的大小主要取决于非线形迭代次数。如果模型只用一次非线形迭代计算就可以收敛,那表明模型很容易收敛,如果需要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) (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: 最大时间步

 HALF: 接近报告步时的时间步取半

 DIFF: 时间步截断CHOP之后的增长

 (在ECLIPSE技术手册的第125页还会找到更多)

 如果模型中有很多CHOP,DIFF,MINS,那模型有严重的收敛性问题。

二:如何解决数模计算的收敛性问题:

 模型不收敛的原因很多,网格参数,属性参数,流体PVT参数,岩石相渗曲线,毛管压力曲线,相渗曲线端点标定,初始化,井轨迹,垂直管流表都会造成模型不收敛,下面分别介绍各部分如何造成不收敛及如何解决。

1。网格部分:

网格正交性差和网格尺寸相差太大是导致不收敛的主要原因之一。正交性差会给矩阵求解带来困难,而网格尺寸相差大会导致孔隙体积相差很大,大孔隙体积流到小孔隙体积常会造成不收敛。

 解决办法:

 网格正交性差通常是在建角点网格时为描述断层或裂缝的走向而造成的。在此情况下,最好能使边

数值模拟应用

界与主断层或裂缝走向平行,这样一方面网格可以很好地描述断层或裂缝,另一方面正交性也很好。  在平面上最好让网格大小能够较均匀,在没有井的地方网格可以很大,但最好能够从大到小均匀过渡。纵向上有的层厚,有的薄,最好把厚层能再细分。在检查模型时应该每层每层都在三维显示中检查。径向局部网格加密时里面最小的网格不要太小。

在ECLIPSE里用MINPV关键字可以把小于设定孔隙体积的网格设为死网格,这样通常会有用。

2。属性参数:

 不合理的插值计算会导致属性分布很差,如果是从地质模型粗化为数模模型,通常问题不大,只是有时候数模人员自己插值时会有问题。

 解决办法:

 有可能尽量用地质模型的数据,自己插值时可以加一些控制点使属性合理分布。X,Y方向的渗透率最好相等或级差不大。在井连通网格的Z方向渗透率不要设为0,如果想控制垂向流动,可给一个很小的值。

3。流体PVT参数

 流体PVT参数会有两种可能的问题,一是数据不合理导致了负总压缩系数,二是压力或气油比范围给的不够导致模型对PVT参数进行了外插。

 解决办法:

 检查PRT文件中的WARNING信息,如果在油藏压力范围内有负总压缩系数的警告,应该修改PVT参数,否则的化会有收敛性问题。如果负总压缩系数是在油藏压力范围之外,可以忽略该警告。此部分的修改主要可以小规模修改油和气的FVF和RS。关于ECLIPSE如何计算负总压缩系数,参照ECLIPSE技术手册第936页。

 在ECLIPSE中加EXTRAPMS关键字可以要求输出如果发生PVT插值后的警告信息。在提供PVT表时,压力应该覆盖所有范围,包括注水后的压力上升。RS值也应该考虑到气在油中的重新溶解。

4。岩石相渗曲线和毛管压力曲线。

 ECLIPSE不会对输入模型的相渗曲线和毛管压力曲线进行光滑,将会应用每一个输入饱和度和相渗值,所以要保证输入的参数是合理的。通常的问题有:

 (1): 饱和度和相对渗透率的数据位数过多。

 (2): 饱和度值太接近,导致相渗曲线的倾角变化很大。

 (3): 饱和度有很小变化但相对渗透率发生了很大变化。

 解决办法:

 (1):饱和度和相对渗透率最多给两位小数就够了。

 (2):检查相渗曲线的导数?梢杂τ肊CLIPSE中的SCAL模块),导数要光滑。

 ??3):将临界饱和度和束缚饱和度设为不同的值。

数值模拟应用

5。端点标定:

 在应用端点标定时,有时标定完后的相渗曲线倾角很大,标定后的毛管了很大。

 解决办法:

 在三维显示中检查标定完的PCW,可以给PCW一个最大值来控制毛管压力。

 输出每个网格标定后的相渗曲线进行检查。

6。初始化:

 初始化最容易发生的问题是在初始时模型不稳定,流体在初始条件下就会发生流动,这也会导致模型不收敛。造成模型初始不稳定的主要有:

 (1) 手工赋网格饱和度和压力值。

 (2)拟合初始含水饱和度。

 解决办法:

 (1)尽量不要直接为网格赋压力和饱和度值,尽量由模型通过油水界面及参考压力来进行初始化计算。

 (2)要想拟合地质提供的初始含水饱和度分布,应该进行毛管压力的端点标定,这样毛管压力会稳住每个网格的水,在初始条件下不会流动。

 (3)可以通过让模型在没有任何井的情况下计算十年来检查初始条件下模型是否稳定,如果10年的计算模型压力和饱和度度没有变化,说明模型初始是稳定的。

7。井轨迹:

 在进行井处理时井可能以之字型在网格中窜过,有可能发生井的实际窜过方向与模型关键字定义的方向不符,这也会导致不收敛。

 解决办法:

 (1) 在三维显示中检查井轨迹。

 (2) 如果井已经关掉,在模拟时不要给零产量,要用关键字把井关掉。

 (3) 检查井射孔,井不要射在孤立的网格上。

8:垂直管流曲线:

 有了垂直管流曲线很容易导致模型不收敛,这有两种可能:

 (1)曲线有交叉。

 (2)曲线发生了外插。

 解决办法:

 (1) 用前处理软件(ECLIPSE中的VFPi)检查曲线。

 (2)在ECLIPSE中加EXTRAPMS关键字可以要求输出如果发生VFP插值后的警告信息。

 (3)曲线应该覆盖所有井口压力,含水,油气比及产量。

数值模拟应用

 (4)在ECLIPSE用WVFPEXP。

9:其他解决办法:

如果模型数据没有问题,可以调整模拟器的收敛计算参数,对于ECLIPSE,可以做以下调整:

(1)调整TUNING中的最大时间步。如果模型每计算到30天就会截断时间步,可以将最大时间步调整为20天,这样计算会快很多。

(2)调整TUNING中的最大线形迭代次数到70次。

(3)降低TUNING中的线形收敛误差标准

(4)对于组分模型用FREEZEPC和DPCDT.

数值模拟应用

 因为师兄给了我一些地质图件和比如原始地质储量的某些参数,让我完成地质模型,但我不知道从何处着手。你可不可以详细地介绍一下有效模型的建立步骤(就以没有任何可引用的文件为例,也就是说参数都要自己输入),我还投了邮件在你邮箱里,请多关照!谢谢!

解答:

 你的问题在国内比较普遍。问题主要在于各部门之间缺乏合作。很多油田没有项目组的概念,勘探部门做完的项目把结果交给开发部门任务就算完成,而且结果经常以图形和报告形式,而不是开发部门更需要的数据文件,这样导致开发部门无法有效共享勘探部门的研究。这里面有软件的一体化问题,最重要的是管理问题。

 现在说你的问题。我不知道你拿到的是井点参数还是等值图。假如你拿到的是井点参数,而且你用的软件是Flogrid,以下是步骤:

1。 进入Flogrid的TOOLS-2D MAPPING CARVAS.

2。 在2D MAPPING CARVAS窗口,进入CREATE-SCATTER SET。

3。 选择属性。对于数值模拟,你应该完成每个储层孔隙度,渗透率,厚度,净厚度,深度等值图。 4。设置你图形的范围,最大最小值。

5。 对应每一井点,输入其数值。

6。 进入CREATE MESH MAP

7。选择属性,选择SCATTER SET, 生成MESH MAP.

8. 将所有图件全部完成后, 参照FLOGRID手册中的练习十: Generating structure and properties from 2D mapping.

FloGrid 网格生成器

FloGid是一个一体化的产品,它支持三维油藏描述标准格式RESCUE(POSC)输入。它也可以直接读入油藏描述软件生成的二维地质图,在其基础上自动(或手工)形成油藏模拟网格系统。

 生成可由Eclipse直接使用的各种油藏模拟网格系统(正交、径向、角点、PEBI)

 支持三维RESCUE(POSC标准)格式输入

 直接从大多数绘图软件包和油藏描述系统读入数据

 支持多相upscaling(物性粗化)

 支持三维非规则网格的局部加密

 在处理断层交叉面时, 可采用分段坐标系

 采用多点流动方式来反映网格间流动的非正交效应

 在建立油藏模拟网格系统时能自动识别和考虑地质特

 全三维数据(井、各种图、地质网格、断层和模拟网格)可视化

 提供了一组把地质模型或随机模型的细网格粗化成油藏模拟网格的工具

数值模拟应用

Eclipse Office数模一体化管理模块

EclipseOffice可以打开和管理Eclipse数模家族的任意软件,允许你在数模运行中随时查看计算结果;可以编辑和评价数模计算结果,并且可生成结果报告;可以快捷的建一个数据研究模型并进行计算。EclipseOffice是一个一整化的桌面环境,.提供了五个特色模块,给用户控制管理 数模流程提供了极大的方便。

 项目管理-在Office环境下管理特定的模拟研究项目

 数据管理-以数据面板的形式列入所有ECLIPSE的标字号,并在特定的油田条件下,列出 有关标字号的含义、数值和单位

 运行管理-启动及监测模拟运行

 结果显示-显示曲线和二维、三维结果

 报告输出

Schedule 生产动态数据和完井数据准备

 精确地、高效地为ECLIPSE准备井的生产数据和完井数据

 以交互方式工作

 快速地定义井、 井组的逻辑结构

 灵活地定义时间步

 自动计算完井数据

 自动生成ECLIPSE所需的数据卡片

PVTi 流体PVT数据分析

PVTi是根据相应的状态方程、样品测试生成和分析PVT相态特性

 是一个通用的、基于状态方程的、为Eclipse准备输入数据的PVT数据分析软件包

 具有众多的状态方程模型和分析方法,可对流体样品进行分组或拟组分化

 可以交互方式或批处理方式工作

 四个可用的状态方程:

Redlick–Kwong

Sovae–Kwong

Peng–Robinson

Zudkevitch–Joffe

 两种粘度计算方法:

Pedersenet.al

Lohrenz–Bray–Clark

 三种流体定义方法:

数值模拟应用

 标准(Library)组分只需输入相应的组分即可

 特征组分

 拟组分

 物质平衡检查

 试验过程的模拟

 实验数据用状态方程进行回归

 生成Eclipse所需的PVT数据

SCAL 准备相渗透率和毛细管压力数据

 读入实验室数据,完成诸如曲线光滑、分组等工作

 把实验室数据转换成ECLIPSE所需的格式

 为用户提供了一个开放的界面和编程语言

 数据的规一化、平均或内插

 生成模拟所需的关键字

 3D显示模拟网格及其属性

 端点值的调整

VFPi 井筒垂直管流计算



 为Eclipse计算出所需的复杂的举升曲线

 六种多相流计算方法用于计算压力剖面

 可处理水平井、 地面流动管线和垂直井

 可考虑气举、井底泵、气体压缩机及地面油嘴的影响

 黑油及组分流体计算能力

SimOpt自适应式历史拟合

Simopt是一个全新的历史拟合软件。它通过计算数模模型与实际历史数据间的梯度差,借助RMS均方差的调整来达到拟合。

 一次可同时对多个参数(>500)进行梯度计算

 运用均方差、汉森(Hessian)矩阵、相关性矩阵、协方差矩阵对结果进行分析

 引入梯度带分析技术对地质模型进行优化

 可充分利用地质统计学的成果

 应用置信度限制,通过线性预测分析,进行趋势历史拟合

FloViz三维可视化

 FloViz支持各种网格形式,包括:块中心网格、角点网格、PEBI网格、笛卡尔和径向加密网格及

数值模拟应用

粗化网格。

 交互式地对图形进行旋转、平移、缩放、光照、多重照相、透视和透明功能。

 各种(井、井组及全油田)的曲线图的显示

 诸如断层、剖面等各种网格图的显示

 饱和度、压力等计算结果的等值图

 监测模拟的运行

 提供了一组宏命令来重复经常性工作

 可进行水驱和气体锥进的动态显示、可进行颜色填充的场图显示、箭头绘图,说明流动方向和主要的传导率显示功能。

 显示任一时间步计算结果

 动画显示运算结果

 观察油藏的内部切片

 利用门限值显示任意部分网格单元

 放大显示油藏模型

 创建和修正彩色场图

 显示油气水三元相图

 任意参数的区间显示

 完井及水平井显示

 用户可按自己的爱好来定义或修改颜色和注释

问题00001:

gulf moon,您好! 我是新手,刚安装了微机版,不会用,是先用FLOGRID建立地质模型 吗? 可是点击FLOGRID,出现目录,不知道如何选 ?里面有例子吗?我把出现的界面 附后!

回答:

你最好先学ECLIPSE Office, 你会知道ECLIPSE基本流程。

如果是微机版,你应该有随机手册。你可以做Office练习二。

你可以参照我的Blog里关于ECLIPSE初学者常见问题。

问题00002:

你的数值模拟工作者天地我几乎天天都上,看看有没有什么新信息。您辛苦了!向您致敬!

我是河南油田搞数模和建模的,由于搞的时间不长,所以问题较多。

我有个问题想请教你一下,就是Eclipse在初始化模拟以后怎样把含油饱和度、孔隙度、渗透率、 NTG等物性模型输出,为导入Petrel所用呢?或者是各个时间步长的物性导出模型为Petrel所用。 因为我建模时没有键饱和度模型。希望能与你常交流!祝您工作顺利!

数值模拟应用

回答:

ECLIPSE计算输出文件为:

.GRID: 网格文件

.INIT: 属性文件,包括孔隙度,渗透率,NTG等

.UNRST或.X0000: 动态结果,包括饱和度,压力等。

要想将ECLIPSE结果输入到Petrel, 你需要先运行ECLIPSE生成以上文件。然后就可以直接输入到Petrel了。

有问题可以再联系我。

问题00003:

我按照你给我的指导,已经成功的把.GRID: 网格文件 .INIT: 属性文件加载进了Petrel中了!

但是我用的是Petrel2002SE版本,好像不能加载动态结果。最新的版本应该能够加载,我到时再试一下。 另外我还要请教你的是,我是用Eclipse做数模的,历史拟合老是出现计算出来的油水没有实际的多, 我用的是定油量拟合,问题大部分出现在前期拟合的不好,后期拟合的较好,我修改NTG 和相渗曲线, 但是没什么大的变化, 咨询了别人,他们说什么从水体的能量和渗透率等方面考虑,我不知道怎样才能把水体的能量加强呢?Case Definition 中我选择的是解析水体,另外好像就没有水体相关方面的设置了呀?所以还得请教您呀!我们的驱块注过聚合物,目前为聚后水驱。您的网站和回复使我学习到了很多知识呀!非常感谢!

回答:

1。你的模型是定油量生产,如果模型计算油量达不到你设定的油量,那说明模型此时已经受井底压力的控制了。

这可能有几方面的原因:

a. 参考面的压力太低。(EQUIL)

b. 渗透率太低。 (PERMX,PERMY,PERMZ)

c.井污染系数太大。 (COMPDAT)

d.油粘度太大。 (PVTO)

2。当模型能够满足你的油量后,如果产水量偏低,可以考虑以下方面:

a. 相对渗透率的形状,考虑增加前期水相相对渗透率。(SWFN或SWOF)

b. 水体。你选完水体后在Initialization部分还需要定义水体,包括水体大小,

属性以及水体与模型的连接。 (参照关键字:AQUCT,AQUANCON)

c. 模型渗透率及油水粘度比。

你最好先定油藏产量(RESV)拟合油田和单井压力水平,然后定油量拟合含水。

另外每次参数调整最好有所依据,一定要在合理范围之内。

数值模拟应用

问题00004:

高手您好!现在有个关于模型粗化方面的疑问想请教一下,我以前只是做到建立三维储层模型就结束,没有跟后面的数模结合过,所以现在想请教一下,对三维模型进行粗化后的结果是不是仍是三维网格,只是网格节点数目减少了很多?

特别是垂向网格的粗化,比如说,两层之间的垂向多个网格是不是在粗化后垂向上仍有多个网格而非单个网格?

只是网格数目减少了那?谢谢!!

我是地质大学(北京)的一名在读学生,很感谢您的博客网站,希望能多多交流!

回答:

首先谢谢你访问我的网站!

你说的对,网格粗化的目的就是为了减小网格数。数值模拟模型要进行模拟计算,但由于计算机内存,计算速度的限制或由于项目时间的要求,不可能直接应用地质模型。通常地质模型网格的数量级在百万或千万,但数模模型的数量级在十万或百万。将地质模型合并为数模模型的过程就是粗化。 粗化分两步,先是网格的合并,然后是属性的合并。网格合并可以是平面,也可以是垂向。 在垂向合并时不要将隔层与储层合并在一起。属性合并包括孔隙度,渗透率,净毛比等。

问题00005:

高手,你好!

谢谢回复!

我想问一下,如果将垂向网格粗化至一个,即变换为平面图以后,这个平面图与按三维网格提取的有效厚度平面图有什么区别呢?或者说,是不是根据三维网格提取的有效厚度平面图只是一个单纯的垂向有效网格节点的累加,而粗化后的平面图程是根据不同属性对应的不同算法得到的呢?也就是说,粗化的平面图,其计算过程要比提取有效厚度平面图的过程复杂?

如果是这样,那么,这两个平面图有什么区别吗?你认识哪个能更好的反映储层有效厚度的平面展布规律呢?

回答:

这里面概念不同。比如你把100×100×10的地质模型粗化为100×100×1的数模模型,此时每个数模网格都具有孔隙度,渗透率,厚度,深度和有效厚度属性。对于每个属性,粗化的方法不同,孔隙度的粗化采用体积加权法,而渗透率粗化则要采用流动计算法,有效厚度的粗化要考虑到孔隙度的截至值,应保证孔隙度,渗透率和有效厚度粗化后的一致。你所说的有效厚度平面图通常针对某小层,而且通常是等值线图,往往是将井点处某小层有效厚度相加然后平面插值生成。与粗化是不同概念。简单的说,粗化就是为了数模计算节省计算时间和内存占用而采用的一系列数学计算方法。

问题00006:

数值模拟应用

海湾月光:

你好,我想学习数模软件eclipse,能不能告诉我从什么开始学习比较好

回答:

我建议你先从ECLIPSE OFFICE开始。你需要知道ECLIPSE计算所需要的几部分参数,包括模型,流体,相渗,生产数据等。

然后学习如何通过三维饱和度,压力场或计算曲线分析计算结果。你可以先完成ECLISPE OFFICE里面的练习一和二。

完成此两个练习后再与我联系。

问题00007:

老师你好:

我要模拟一个井组的话,是不是应该把周围的一些井也建到模型中去?然后在裁剪?如果是的话具体应该怎么做?谢谢!

回答:

在你模拟一个井组时,此井组肯定要受到附近井组的生产或注入的影响。如果你有全油田模型,你可以通过流动边界的方法来模拟一个井组。所谓流动边界是指在你进行全油田模拟计算时,ECLIPSE将记录每时间步流体流入或流出你设定边界处的流量。生成流动边界后你就可以对单一井组进行模拟而不需要模拟全油田。所以在你建立全油田模型后,先运行DUMPFLUX生成FLUX文件,然后运行USEFLUX模拟井组。你可以参照ECLIPSE技术手册中关于流动边界(FLUX BOUNDARY)的介绍。

问题00008:

你能给我把PEBI网格作一下介绍吗?有ECLIPSE里面怎么用这个东西.在FLOGRID里这们的例子有没有?

回答:

PEBI是新一代网格体系,最早的矩形网格是为了满足计算的需要,但不能很好地描述油藏真实构造, 角点网格的推出弥补了此缺陷,但角点网格由于网格不正交,给计算速度和精度带来影响。PEBI完善了网格描述,不但能模型油藏真实构造,而且保证正交。 现在PEBI技术还没有完全成熟,不过你在FLOGRID中输出PEBI网格后,可以INCLUDE到你的GRID部分,就可以计算了。

关于FLOGRID如何建立PEBI网格,可参照FLOGRID联系4:Building Unstructured simulation model,以及ECLIPSE OFFICE练习3:Constructing a PEBI simulation grid.

问题00009:

我是ECL的初学者,想向你请教FLOGRID的问题,我有测井的相关数据、井位坐标、井斜数据,如何建立地质模型?你能给我讲解一下具体的操作步骤和提供一些FLOGRID的学习资料吗?谢谢!

本文来源:https://www.bwwdw.com/article/oi8j.html

Top