Orkiszewski算法的C语言实现.pdf

Orkiszewski 算法的 C 语言实现 Orkiszewski 的C 语言实现 中国石油大学（华东） 石油工程学院 Orkiszewski 方法是计算气液两相垂直管流的方法，由于 中间步骤涉及大量计算及迭代，使用计算机编程是实现计算的 简便方法。笔者使用 C 语言编程。整个代码分为三大块，主 序、综合函数、功能函数。 编程思路：多相管流计算所需的具体参数已经在程序中给 出。采用从井底向上计算，按深度增量迭代。读者可以根据自 己的需要修改参数和迭代方法。 由于在一口井中一般很难同时出现四种流态。对于低含气 井只会出现泡流或段塞流。下面所给程序中，没有给出过渡流 的计算代码，请读者根据自身需要加以完善。 编程中的主要难点是参数过多，所需计算的物性参数很 多，导致定义变量时容易造成混乱。由于计算物性参数的程序 和 Orkiszewski 程序的编写是分开进行的，讲两者很好的结合 又是一个很大的问题。读者可以采用代码复用或文件包含的方 法将二者结合。限于笔者的水平，以下所给的程序是将两者放 在一个程序里，导致变量很多且混乱。请读者自行完善！如有 好的建议，欢迎来信！ 参考资料： 1. 张琪.采油工程原理与设计.—东营：中国石油 大学出版社 第 1 页 Orkiszewski 算法的 C 语言实现 2. 步玉环、谷建伟、薛建泉.石油工程综合设 计.—东营：中国石油大学出版社 zxyzgsydx@ //奥奇斯泽斯基算法多相管流计算。质量含水率 0.44，初 始产液量设定为 50t/d. //假设给定温度、压力、液体流速、流体密度、流体表面 张力、管子直径、溶解气油比、生产气油比、压缩因子、原油 体积系数 //vs 为滑脱速度：0.244m/s #include stdio.h #include math.h double tempc(double Q,double H,double h,double tr,double t0,double Fw); double Rsc(double gamao,double gamag,double ta,double pa,double Rp); double Boc(double Rs,double gamag,double gamao,double ta); double dopc(double gamao,double gamag,double Rs,double Bo); 第 2 页 Orkiszewski 算法的 C 语言实现 double Zc(double gamag,double ta,double pa); double vlc(double gamaAPI,double ta,double Rs,double fw); double vgc(double gamag,double pa,double Z,double ta); double sigmac(double ta,double pa,double gamaAPI,double fw); double zhhs(double p1,int j);//定义综合函数。 double dz=100,H=2400;//定义两个全局变量，计算段长度 和油层中深 int main() { char c; int j=0; double p1; for(j=1;; j++) { if(j==1) { printf(Please enter