数值模拟入门指南

《数值模拟入门指南》由会员分享，可在线阅读，更多相关《数值模拟入门指南（23页珍藏版）》请在金锄头文库上搜索。

1、黑油 第一：从掌握一套商业软件入手。 我给所有预从事油藏数值模拟领域工作的人员第一个建议是先从学一套商业数值模拟软件开始。起点越高 越好，也就是说软件功能越强越庞大越好。现在在市场上流通的ECLIPSE,VIP和CMG都可以。如果先学 小软件容易走弯路。有时候掌握一套小软件后再学商业软件会有心里障碍。 对于软件的学习，当然如果能参加软件培训最好。如果没有机会参加培训，这时候你就需要从软件安装时 附带的练习做起。油藏数值模拟软件通常分为主模型，数模前处理和数模后处理。主模型是数模的模拟器， 即计算部分。这部分是最重要的部分也是最难掌握的部分。它可以细分为黑油模拟器，组分模拟气，热采 模拟器，流线

2、法模拟器等。数模前处理是一些为主模拟器做数据准备的模块。比如准备油田的构造模型， 属性模型，流体的PVT参数，岩石的相渗曲线和毛管压力参数，油田的生产数据等。数模后处理是显示模 拟计算结果以及进行结果分析。以 ECLIPSE 软件为例，ECLIPSE100,ECLIPSE300 和 FrontSim 是主模拟器。ECLISPE100 是对黑油 模型进行计算，ECLISPE300是对组分模型和热采模拟进行计算，FrontSim是流线法模拟器。前处理模 块有Flogrid,PVTi,SCAL,Schedule,VFPi等。Flogrid用于为数值模拟建立模拟模型，包括油田构 造模型和属性模型；PV

3、Ti用于为模拟准备流体的PVT参数,对于黑油模型，主要是流体的属性随地层压力 的变化关系表，对于组分模型是状态方程； SCAL 为模型准备岩石的相渗曲线和毛管压力输入参数； Schedule处理油田的生产数据，输出ECLIPSE需要的数据格式（关键字）；VFPi是生成井的垂直管流 曲线表，用于模拟井筒管流。ECLIPSE OFFICE和FLOVIZ是后处理模块，进行计算曲线和三维场数 据显示和分析， ECLIPSE OFFICE 同时也是 ECLIPSE 的集成平台。对于初学者，不但要学主模型，也需要学前后处理。对于ECLISPE的初学者，应该先从ECLISPEOFFICE 学起，把ECLIS

4、PE OFFICE的安装练习做完。然后再去学Flogrid, Schedule和SCAL。PVTi主要 用于组分模型，做黑油模型可以不用。第二：做油藏数值模拟都需要准备什么参数 在照着软件提供的安装例子做练习时经常遇到的问题是：虽然一步一步按照手册的说明做，但做的时候不 明白每一步在做什么，为什么要这么做。这时候的重点在于你要知道你一开始做的工作都是为数值模拟计 算提供满足软件格式要求的基础参数。有了这些基础参数你才能开始进行模拟计算。这些基础参数包括以 下几个部分：1。模拟工作的基本信息：设定是进行黑油模拟，还是热采或组分模拟；模拟采用的单位制（米制或英制） 模拟模型大小（你的模型在X,Y,

5、Z三方向的网格数）；模拟模型网格类型（角点网格，矩形网格，径向网格 或非结构性网格）；模拟油藏的流体信息（是油，气，水三相还是油水或气水两相，还可以是油或气或水单 相，有没有溶解气和挥发油等） ;模拟油田投入开发的时间；模拟有没有应用到一些特殊功能（局部网格加 密，三次采油，端点标定，多段井等）；模拟计算的解法（全隐式，隐压显饱或自适应）。2。油藏模型：模型在 X,Y,Z 三方向的网格尺寸大小，每个网格的顶面深度，厚度，孔隙度，渗透率，净 厚度（或净毛比）。网格是死网格还是活网格。断层走向和断层传导率。3。流体PVT属性：油，气，水的地面密度或重度；油，气的地层体积系数，粘度随压力变化表；溶解

6、油 气比随压力的变化表；水的粘度，体积系数，压缩系数；岩石压缩系数。如果是组分模型，需要提供状态 方程。4。岩石属性：相对渗透率曲线和毛管压力曲线。如果是油，气，水三相，需要提供油水，油气相对渗透率 曲线和毛管压力曲线（软件会自动计算三相流动时的相对渗透率曲线）；如果是油，水两相或气，水两相， 只需要提供油水或气水两相相对渗透率曲线和毛管压力曲线。5。油藏分区参数：如果所模拟的油田横向或纵向流体属性，岩性变化比较大，或者存在不同的油水界面， 这时需要对模型进行PVT分区（不同区域用不同的PVT流体参数表），岩石分区（不同区域用不同的相对 渗透率曲线和毛管压力曲线）或者平衡分区（不同平衡区用不同

7、的油水界面）。另外如果想掌握油藏不同断 块的储量或采收率，可以对模型进行储量分区（不同储量区可以输出不同的储量，产量，采收率，剩余储 量等）。6。初始化计算参数：油藏模型初始化即计算油藏模型初始饱和度，压力和油气比的分布，从而得到油藏模 型的初始储量。这部分需要输入模型参考深度，参考深度处对应的初始压力，油水界面以及气水界面；油 气比或饱和压力随深度的变化；如果是组分模型，需要输入组分随深度的变化。7。输出控制参数：即要求软件在计算时输出哪些结果参数。比如要求输出模型计算油田的油，气，水产量 变化曲线；油田压力变化曲线；单井油，气，水产量变化曲线；单井井底压力变化曲线；单井含水，油气 比变化曲

8、线等。8。生产参数：对于已开发油田，这部分的数据量非常大。包括油田每口井的井位，井轨迹，井的射孔位置， 井的生产或注入历史（油，气，水产量，注入量，井底压力，井口压力等），井的作业历史等。第三：如何准备各部分参数 上面介绍了做油藏数值模拟所需要的参数，那么这些参数是如何得来的？又应该如何输入到数模模型中 哪？下面具体介绍各部分数据的来源即处理办法。1。油藏模型： 大多数油公司现在都在采用一体化工作流程，数模工程师不再需要自己去建立油藏模型。地质学家对油田 进行详细的油藏描述工作，基于地震解释数据，测井解释数据，岩芯数据以及结合地质家对油田的认识建 立三维地质模型。根据数模工程师的要求，地质家对

9、三维地质模型进行粗化处理，可以直接为数模工程师 输出符合数模工程师需要的油藏模拟模型。如果你不幸工作在一个没有很好地采用一体化工作流程的油公司或学校，那你就需要自己建立油藏模拟模 型了。这时你需要用到数模软件的前处理建模模块（比如Flog rid）.如果你从地质家处拿到的是三维地质模 型，你的工作量还不是很大，主要是采用前处理建模模块对地质模型进行粗化，生成数值模拟模型。如果 你拿到的是两维数据体，那你的工作量就要大的多。首先你要确定你必须要得到以下几方面的两维数据体：* 每层的顶面深度* 每层的厚度* 每层的孔隙度分布* 每层的渗透率分布* 每层的净厚度或净毛比分布* 断层数据有了这些数据体

10、，你就可以在前处理建模模块中建立三维数值模拟模型。还有更不幸的情况，你有可能工作在软硬件环境都很落后的油田（比如国内的四川油田），你根本连两维数 据体也拿不到，你能拿到的只是一张一张的等值线图，或者只是单井井点数据，这时候你首先需要做的是 将这些等值线图件用数值化仪数值化为两位数据体，或者根据井点数据生成等值线，然后再建三维模型。你在开始建三维数模模型时首先需要确定你的模型应该采用什么样的网格类型。目前数值模拟软件常用的 网格有正交网格，角点网格，径向网格,非结构化网格（比如PEBI）和动态网格等。你需要知道这些网格的不 同之处及其优缺点:A: 正交网格正交网格是最常见网格,也是最早用来描述油

11、藏的网格类型，目前仍然被广泛应用 .由于其计算速度快的特 点,一些大型油气田经常采用此网格类型.有研究认为正交网格计算结果比其他网格精确.正交网格的数学描述也比较简单。以ECLISPE为例，TOPS描述油藏顶部深度，DZ描述油藏每层厚度， DX描述每个网格X方向长度，DY描述每个网格Y方向长度。B: 角点网格 角点网格的特点是网格的走向可以延着断层线，边界线或尖灭线，也就是说网格可以是扭曲的。这样角点 网格克服了正交网格的不灵活性,可以用来方便地模拟断层,边界,尖灭.但由于角点网格网格之间不正交,这 种不正交一方面给传导率计算带来难度,增加模拟计算时间,另一方面也会对结果的精度有影响.角点网格

12、的数学模型很复杂，必须由前处理软件来生成。以 ECLISPE 为例， COORD 用来描述模型网格 的顶底坐标线(X,YZ), ZCORN描述每个网格八个角点的深度。C: 径向网格 径向网格比较简单，主要用于单井模拟。径向网格可以更合理的描述井附近流体的径向流动。D：非结构网格(PEBI网格)PEBI网格源于1908年就产生的Voronoi网格起主要特点是灵活而且正交.PEBI网格体系提供了方便的 方法来建立混合网格,比如模型整体采用正交网格 ,而对断层,井,边界等采用径向 ,六边型或其他网格 .网格 间的传达率可以自动计算.PEBI网格的灵活性对模拟直井或水平井的锥进问题非常有用另外PEBI

13、网格可 以用来精确模拟试井问题还有PEBI网格降低了网格走向对结果的影响.PEBI网格的缺点是矩阵比其他网 格要复杂的多,需要更加有效的解法.E: 动态网格动态网格是指网格可以随时间而改变.通常用于动态网格加密或动态粗化 .比如说在井生产时采用局部加密 而当井关闭时则采用正常网格.有时候建立全油田整体模型后,对于压力及饱和度变化快的区域,常常需要进行局部网格加密.采用局部加密 可以准确的描述井附近流体的细微变化。网格局部网格可以是正交网格，或是径向网格.Aziz认为(JPT 1993 年)在正交网格中进行正交网格局部加密,有时并不会对结果有改善.他建议采用混合网格,及在正交 网格内采用径向网格

14、加密,这样可以精确地模拟含水和气油比的变化规律.知道了这几种网格类型，那么在你建立模型时应该选择什么网格类型哪？你在建立网格时又应该注意些什 么哪？A: 在条件许可情况下尽量采用正交网格,而且尽量使网格保持均匀.尽量避免大网格直接连接小网格,这样 会带来严重的收敛问题.如果你的模型很大，最好采用正交网格。B: 角点网格已经非常成熟,但在建立角点网格时不要过分扭曲网格，要尽量保持模型的正交性。如果你的 边界与你的主断层相对平行，那么建立的网格系统正交性会比较好。你在建立网格后可以用前处理软件计 算模型正交性。C: 目前 PEBI 网格在解法上还不成熟,应避免使用.但相信在五年内它会成为主导网格.

15、D: 使用局部网格加密要小心 ,最小的网格不能小于井半径.而且局部网格加密部分要覆盖饱和度变化大的 网格.如果是水平井，局部网格加密要覆盖水平段。E: 网格越多模拟结果就越精确的概念是不对的.可以建立单井模型研究多大网格尺寸足够描述地质上的非 均质性.F: 网格走向会影响计算结果. 在天然裂缝油气藏,要使网格走向与主裂缝方向一直.G：DX/DY 应接近于1，不要大于 3。H: 井之间应有不少于三个网格。I: 模型的属性分布也很重要。尤其是数模人员自己插值，属性分布经常不合理。我见过许多模型，井只是 穿过单个的网格，井穿过网格的附近却全是死网格，这时模型计算根本不能收敛。数模人员如果需要自己 建

16、模，最好要找地质人员参谋，或者自己有很好的地质背景。初学者还有几个基本概念需要掌握：A：网格的I,J,K在做模拟时你需要知道每个网格的位置。模型每个网格的位置用网格的I,J,K表示。在ECLISPE软件中， 左上角的第一个网格为（1, 1, 1）,在X方向的第二个网格为（2, 1, 1）,在Y方向的第二个网格为（1, 2, 1）,在Z方向的第二个网格为（1, 1, 2）。其他类推。B：死网格 死网格是模型中不参予流动计算的网格。通常把模型中的泥岩设为死网格。模型自动会把孔隙度和净厚度 为 0 的网格设为死网格。如果你的油藏水体很大,你也不需要把水体都建在模型中。你可以把水体部分设 为死网格,然后用解析水体模拟油藏水体的贡献。C: 在数值模拟模型中如何正确应用有效孔隙度,净毛比？ 要明确区分定义.总孔隙度： 孔隙体积占总体积的百分数。有效孔隙度： 连通孔隙体积占总体积