相态测试理论与技术2.ppt

4 重馏分特征化处理技术 4.4 基于相态实验数据的最优化拟合方法 问题的提出 应用状态方程方法进行油气藏流体相态模拟计算， 必须能明确给出模拟所需的各烃类组分的计算参 数如分子量，临界温度，临界压力，偏心因子，以 及状态方程的可调参数 ，二元交互作用系 数kij等。根据前面的分析，一般对于低碳原子数组 分和非烃类气体，都有准确数据可查并可存储在计 算程序中供计算选择；但对于高碳数（如Cn+重馏 分）组分，虽经过重馏分特征化处理后可得到初始 数据，但直接使用初始数据计算可能会产生较大 的 误差。因此在进行完整的相平衡计算之前，需通过 对相态实验数据的拟合来调整重馏分的热力学参数， 以便使计算结果更好地符合实验数据而达最佳匹配。 4 重馏分特征化处理技术 4.4 基于相态实验数据的最优化拟合方法 拟合计算方法的选择 拟合计算一般采用两种方法。一是自动拟合法， 即以实测相态数据为目标函数，采用非线性回 归方法如单纯形最优化拟合等方法进行热力学 参数的匹配；二是采用经验法，即凭人工对各 种热力学参数对计算过程的敏感性的认识，确 定参数的调整范围，通过多次试凑计算，以获 得与实验数据相匹配的结果。两种方法各有其 优劣，实用时常先采用人工法进行快速试凑以 逼近收敛域，然后通过非线性回归计算得到精 确值。 4 重馏分特征化处理技术 4.4 基于相态实验数据的最优化拟合方法 非线性回归目标函数的建立 式中：F(X1,X2,…,Xn) -- 目标函数； X1,X2,…,Xn -- 所选择的拟合变量(Tcn+， Pcn+, 等）； Wi -- 第i个拟合变量的权重； Yci,Yei -- 分别为计算和实测的 相态参数(如露点压力， 泡点压力，反凝析油 饱和度，采收率等)。 4 重馏分特征化处理技术 4.4 基于相态实验数据的最优化拟合方法 拟合变量的选择 油气体系中重馏分的临界温度，临界压力，偏心 因子和重馏分与C1H4的二元交互作用系数。 状态方程中可调的 等参数 拟合目标变量（实测相态数据）的选择 实测的油气体系露点压力，泡点压力； 等组成膨胀，定容衰竭，地面分离器闪蒸等数据。 5 相平衡计算数值计算方法 问题的提出 相平衡计算数学模型 在数学上属于二次正定型的高度非线性化的方程组 5 相平衡计算数值计算方法 相平衡计算非线性方程组数值算法选择 牛顿-拉夫森迭代算法=》求解雅可比矩阵=》解析式 拟牛顿迭代算法=》差分矩阵=》适于隐式 逐步迭代算法=》非矩阵直接算法=》适于隐式 两类相平衡计算模型及数值计算方法 5.1 适合于采用牛顿-拉夫森算法的相平衡计算模型 5.1.1 露点计算模型 牛顿-拉夫森相平衡算法的特点 同时对计算变量P（或T或V）和xi 等n+1个变量进行迭代 （1）首先把物料平衡方程转化为压力的显式 露点特征：V 1，L 0 同时对计算变量P（或T）和xi 等n+1个变量进行迭代 5 相平衡计算数值计算方法 5.1.1 露点计算模型 （2）建立露点计算相平衡计算模型 5.1.2 泡点计算模型 （1）首先把物料平衡方程转化为压力的显式 同时对计算变量P（或T）和yi 等n+1个变量进行迭代 泡点特征：V 0，L 1 5 相平衡计算数值计算方法 5.1.2 泡点计算模型 （2）建立泡点计算相平衡计算模型 5.1.3 闪蒸计算模型 对计算变量L(或V)和x1,…xi,…,xn 等n+1个变量进行迭代 闪蒸计算特征： 需把n个逸度相等方程和第 n+1个物料平衡方程联立求解 5.1.4 牛顿-拉夫森算法结构分析 （1）露点计算 同时对计算变量P（或T）和xi 等n+1个变量进行迭代 相平衡方程组的矩阵结构 5 相平衡计算数值计算方法 5.1.4 牛顿-拉夫森算法结构分析 （1）露点计算 迭代变量的矩阵结构 牛顿-拉夫森算法迭代格式 5 相平衡计算数值计算方法 5.1.4 牛顿-拉夫森算法结构分析 （1）露点计算 Jacobi 矩阵的结构 5 相平衡计算数值计算