天然气液化中气液相平衡计算研究精要.pptVIP

天然气液化中气液相平衡计算研究 论文主要内容 1.文献综述（天然气液化，相平衡研究） 2.相平衡模型建立与方程选择(立方型模型，SRK方程） 3.算法建立（SRK方程，逸度方程） 4.计算机编程（C++编程） 5.工程实例分析（简单，复杂对比） 6.结论 论文研究思路和步骤 1.文献综述： ①对天然气液化总述（意义，优势） ②天然气液化方式（三种：阶式制冷工艺、混合制冷工艺和膨胀制冷工艺。使用最多：丙烷预冷混合制冷剂液化流程。） ③我国研究现状 相平衡分类 一些常见相平衡方程： ①范德华方程 ②活度系数方程 ③二阶维里型方程 ④RK方程 相平衡分类 相平衡比较 范德华单流体混合法则 综合选择 ①状态方程选择： SRK方程：形式简单，改善了气液相逸度计算效果，能解析求根，所需基础数据少 ②混合法则选择： 范德华单流体混合法则：烃类适用，形式简单，经验丰富 算法建立 重要公式： 1.SRK方程 2.压缩因子方程 3.逸度方程形式 4.相平衡常数方程： 编程思路与方法 1.方程简化： 2.内嵌算法：压缩因子方程（试算法） 3.循环节控制：迭代精度，初值选取 1.方程简化 double S[m_size],X1[m_size],AA[m_size]; for (int i=0;im_size;++i) double XK=0.0; S[i]=(bi[i]/b)*(m_Zl-1)-log(m_Zl-B); for(int j=0;jm_size;++j) XK+=pow(ai[i]*ai[j],0.5)*Zi[j]; X1[i]=(A/B*(2/a*XK-bi[i]/b))*log(1+B/m_Zl); AA[i]=S[i]-X1[i]; 2.内嵌算法 double a1=1,b1=-1,c1,d1,x=0,v; a1=1; b1=-1; c1=A-B-pow(B,2); d1=-(A*B); double Z[m_size]; double m_Zl=0,m_Zv=0; int X,f=0; for (X=0;X=1000000;x=(++X)/1000000.0) v=a1*x*x*x+b1*x*x+c1*x+d1; if (v-0.000001 v0.000001) Z[f]=(double)x; //printf(%f\n,x); ++f; 工程实例分析 方法：对比研究，误差分析 1.简单组分 2.复杂组分 误差分析：最大，最小，平均误差 简单组分实验数据 简单成分计算结果 复杂成分基本数据 复杂成分计算结果 误差分析 结论 （1）通过绪论和各种方程，模型的对比研究，我们可以看出，SRK方程在不失立方型方程形式简单的情况下，大大的改善了计算气液相平衡的效果，并且对方程所需的基础数据少，能够很好的获得纯组分及混合物的气液相组成摩尔分数，气化率等。 （2）根据实验数据与计算值的对比，我们可以看出，SRK方程对简单成分的计算精度较高，实验一中，各数据平均误差控制在7.00%以内，对复杂成分计算精度略为下降，实验二中，最高平均误差上升至9.83%。 结论 （3）SRK计算液相组分摩尔分数误差小于气相组成摩尔分数。由数据可以看出，两个实验中，气相组分摩尔分数误差分别为：5.58%，9.84%，4.67%，液相组分摩尔分数误差分别为：2.52%，3.74%，4.12%。 （4）SRK方程对轻烃的预测相对重烃要准确。在较轻组分的计算中呈现负偏差，而在较重烃中出现正偏差居多。这意味着运用在工程实际现场，会导致轻质烃成分偏大，重质烃成分偏小。。 本文不足和缺陷 1.算法编程中，并未对天然气状态进行判断 （是否处于气液两相区，如何判断，在算法中如何实现收敛） 2.SRK方程用于复杂组分计算结果还不太令人满意，如何继续提高精度，有待进一步研究 谢谢！ 1.简单成分 2.复杂成分 * * * * * WE ARE SEEING THE WORLD IN COLOR * 指导老师： 相平衡研究发展 1.相平衡定义与判定（定义，逸度判定法） 特征： 2.相平衡研究进程（三大阶段） 范德华状态方程 RK等方程，对状态方程一些列修正和改进 具有统计力学基础的状态方程 相平衡发展三大阶段 1 2 3 相平衡计算 状态方程法 活度系数法 立方型 范德华 维里型 SRK RK 活度系数模型 VanLaar模型 PT 虽然对各状态方程的评价都有实验数据支撑，但是做出绝对公正的判断是很困难的。一方面，状态方程的准确性很大程度上依赖于状态方程参数的优选；另一方面状态方