切换到宽版
  • 5274阅读
  • 7回复

[分享]数值模拟不收敛问题 [复制链接]

上一主题 下一主题
 
发帖
2422
财富
6580
威望
1501
交易币
0
只看楼主 倒序阅读 使用道具 0楼 发表于: 2008-01-02 | 石油求职招聘就上: 阿果石油英才网
模型不收敛的原因很多,网格参数,属性参数,流体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输出报告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):检查相渗曲线的导数,导数要光滑。
   (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.


作者: gulfmoon

评价一下你浏览此帖子的感受

精彩

感动

搞笑

开心

愤怒

无聊

灌水
思想无罪,准我一生以沉默相许!
离线niyajiala
发帖
1570
财富
1711
威望
12
交易币
0
只看该作者 1楼 发表于: 2008-03-14 | 石油求职招聘就上: 阿果石油英才网
— (起舞弄清影) 帖子下面注明了是“gulfmoon”海湾月光的文章了哦! (2008-07-30 16:01) —
这不是精准石油上的吗?
All is well
离线carol_chenyu
发帖
9
财富
43
威望
3
交易币
0
只看该作者 2楼 发表于: 2008-04-07 | 石油求职招聘就上: 阿果石油英才网
谢谢咯~~~~~~~~~
离线huoshanxq
发帖
27
财富
329
威望
1
交易币
0
只看该作者 3楼 发表于: 2008-05-19 | 石油求职招聘就上: 阿果石油英才网
thank you
离线spearshield
发帖
1178
财富
1517
威望
176
交易币
0
只看该作者 4楼 发表于: 2009-03-24 | 石油求职招聘就上: 阿果石油英才网
— (起舞弄清影) 帖子下面注明了是“gulfmoon”海湾月光的文章了哦! (2009-03-30 19:27) —
呵呵,好支持顶
[ 此贴被spearshield在2009-03-30 19:29重新编辑 ]
离线chensenswpu
发帖
1116
财富
292
威望
15
交易币
0
只看该作者 5楼 发表于: 2009-03-30 | 石油求职招聘就上: 阿果石油英才网
何必那么在乎出处,大家能收益就行了
我爱石油
发帖
2422
财富
6580
威望
1501
交易币
0
只看该作者 6楼 发表于: 2009-05-15 | 石油求职招聘就上: 阿果石油英才网
一、笔算开平方的步骤
1.将被开方数的整数部分从个位起向左每隔两位划为一段,用撇号分开,分成几段,表示所求平方根是几位数;
2.根据左边第一段里的数,求得平方根的最高位上的数;
3.从第一段的数减去最高位上数的平方,在它们的差的右边写上第二段数组成第一个余数;
4.把求得的最高位数乘以20去试除第一个余数,所得的最大整数作为试商
5.用商的最高位数的20倍加上这个试商再乘以试商.如果所得的积小于或等于余数,试商就是平方根的第二位数;如果所得的积大于余数,就把试商减小再试
6.用同样的方法,继续求平方根的其他各位上的数.
二、通过例子掌握开平方的方法
如何求1024的平方根呢?按照上面的6步走一下。
1、把1024每隔两位分开,10’24。
2、左边第一段10最高位是3。
3、10-3*3=1和24组成余数124。
4、124/(3*20)=2,2作为试商
5、(3*20+2)*2=124,正好等于余数。
6、1024的平方根就是32。

笔算开平方运算较繁,在实际中直接应用较少,但用这个方法可求出一个数的平方根的具有任意精确度的近似值.
作者:tshfang
来源: 泥胚文章写作 http://www.nipei.com 原文地址: http://www.nipei.com/article/2134
思想无罪,准我一生以沉默相许!
离线鱼小贝
发帖
77
财富
6
威望
0
交易币
0
只看该作者 7楼 发表于: 2012-01-05 | 石油求职招聘就上: 阿果石油英才网
很有帮助,谢谢楼主

网站事务咨询:QQ:1392013 | 26189883
阿果石油网为免费个人网站,为石油人提供免费的在线即时技术交流场所,拒绝任何人以任何形式在本论坛发表与中华人民共和国法律相抵触的言论和行为!
如有言论或会员共享的资料涉及到您的权益,请立即通知网站管理员,本站将在第一时间给予配合处理,谢谢!