阳光石油网|石油技术交流|石油人论坛

 找回密码
 欢迎注册
查看: 628|回复: 1

[ECLIPSE] eclipse计算时大量显示Zero determinant problem in LAINV是什么情况?

[复制链接]
  • TA的每日心情
    奋斗
    2024-1-13 13:48
  • 签到天数: 240 天

    [LV.8]以坛为家I

    发表于 2019-1-8 15:12:38 | 显示全部楼层 |阅读模式

    马上注册,下载丰富资料,享用更多功能,让你轻松玩转阳光石油论坛。

    您需要 登录 才可以下载或查看,没有账号?欢迎注册

    x
    eclipse计算时大量显示Zero determinant problem in LAINV是什么情况?
  • TA的每日心情
    奋斗
    2016-2-26 13:48
  • 签到天数: 508 天

    [LV.9]以坛为家II

    发表于 2019-1-8 16:05:00 | 显示全部楼层
    本帖最后由 joylin 于 2019-1-10 16:54 编辑

    简单分析下这个问题:
    首先你应该提供更详细的信息,比如模拟显示的关于该问题的详细信息。

    LAINV,这个词大家可能不清楚,其实我也不清楚,但是可以猜测,应该跟解线性方程组时的LU分解有关。
    如图:

    就是最左边的矩阵。
    模拟器在计算时大致是:非线性方程组→求导雅克比矩阵→线性方程组→解线性方程组
    线性方程组解算器中,用到了LU分解,Ax=B的问题变为LxU=B,具体的我也没深入,有兴趣的话自己研究。
    这时,需要对L求逆矩阵,也就是LAINV了,但是这个逆矩阵的行列式为0,这就不开心了。
    至此,Zero determinant problem in LAINV的含义大致清楚了。


    在2012版之前,出现这个问题会当做一个error,模拟器停止,2012以后会截断时间步,当做局部线性方程组不收敛。
    如果这个问题出现的不多,可以不用理会。如何和你一样很多,那就有问题了,最直接的影响就是你的计算结果已经偏离的正确的解,
    也就是说,即使计算完了,你的结果将毫无意义。


    如何解决:
    当成线性方程不收敛的问题解决。
    模拟器:TUNING关键字的Record 3部分控制了非线性方程和线性方程的解算,适当修改。
    静态模型:如果是出现很多线性问题,很有可能是静态模型有问题,做一下质量控制。
    动态模型:如果是概念模型,看看设置的生产参数是否合理。

    个人经验:静态模型出问题的可能性大一些,比如网格畸形,渗透率差异过大,传导率不合理等等。可以先质控下静态模型,MINPV截断下活网格试试。如果你用的是2014及以后的版本,可以尝试CPR解算器试试,也许有改进。

    总之,这种问题,一个模型一个样,很难定性定量分析,个人理解,仅供参考。







    回复 支持 反对

    使用道具 举报

    您需要登录后才可以回帖 登录 | 欢迎注册

    本版积分规则

    QQ|Archiver|手机版|小黑屋|阳光石油网 ( 鲁ICP备2021003870号-1 )

    GMT+8, 2025-1-9 17:28 , Processed in 0.047542 second(s), 20 queries .

    Powered by Discuz! X3.4 Licensed

    Copyright © 2001-2021, Tencent Cloud.

    快速回复 返回顶部 返回列表