科吉思石油论坛

标题: 【FAQ】2.1 问题离散化及解决方案 [打印本页]

作者: 乔聪颖    时间: 2020-6-17 09:27
标题: 【FAQ】2.1 问题离散化及解决方案
本帖最后由 张鹏326 于 2020-7-14 11:34 编辑

2.1 Problem discretization and solution
2.1
问题离散化及解决方案

We are finding solution of differential equations system – 2.1.
我们将寻找微分方程组的解决方案 – 2.1。
Standard three-phase three-component isothermal filtration problem (black-oil) taking into account compositional in main variables – molar densities.
标准三相三组分等温渗流问题(黑油),考虑了主要变量(摩尔浓度)的组成。
Algorithm to find solution:  
查找解决方案的算法:
Time approximation (fully implicit) → half-discrete scheme
时间近似(全隐式)→半离散方案
Space approximation (finite volumes to fulfil conservation) → a system of non-linear equations on each step
空间近似(有限体积以实现守恒)→每步都求解一个非线性方程组
Newton method → a system of linear equations on each iteration of Newton
牛顿法→每次牛顿迭代的线性方程组
BCGStab or GMRES algorithms with preconditioners to solve system of linear equations with non-symmetric matrix → calculating preconditioner on each linear system solution iteration
带有预处理器的BCGStab或GMRES算法来求解具有非对称矩阵的线性方程组→在每次线性系统解决方案迭代中计算预处理器


Parameters which influence calculation
影响计算的参数

Linear problem solution
线性问题解决方案
The linear problem looks like
线性问题如下:
where A is matrix, b – the right-hand side. The simplest solution of this problem is given by minimal residual method, the ancestor of generalized minimal residual method (GMRES). The algorithm looks like
其中A是矩阵,b–右侧。这个问题的最简单解决方案是最小残差方法,即广义最小残差方法(GMRES)的前身。该算法看起来像
Where iterational parameter is chosen to minimize the length of residual
其中迭代参数用于最小化残差的长度
as

where ( , ) – vector scalar product.
其中 ( , ) –矢量标量积。(注:向量的点乘)
The common traits of all linear problem solution algorithms:
所有线性问题求解算法的共同特征:
new solution is a subtraction of previous solution and residual with multiplier :
So, the new solution is very sensible to the value of .
因此,新的解对的值非常敏感。
so the residual (and the error) goes to zero «in average», and the residuals in some blocks (usually near well perforations) may be higher than in other blocks. Different methods minimize different norms, so there may be different answers for wells even for the same accuracy.
因此,残差(和误差)“平均”应趋近于零,并且某些网格块(通常在井射孔附近)的残差可能高于其他网格块。不同的方法会最小化不同的规范,因此即使具有相同的精度,井的求解结果也会有所不同。
Work is distributed evenly between all cores
所有内核之间平均分配任务








欢迎光临 科吉思石油论坛 (http://bbs.colchis.cn/) Powered by Discuz! X3.2