石油工程课程设计c语言,东北石油大学-石油工程抽油设计C语言编程.doc

float s=3;//s:抽油机冲程

#include #include #include #define PI 3.1415926 #define g 9.8float H=1545;//H:油层深度float D〇=0.0889;//Do:油管内径float Dt=0.19;//Dt:套管直径float 0=0.0323;//0:地温梯度float t0=79.8;//to:井底温度float Ph=10.75;//Ph:地层压力float Pb=9.98;//Pb:饱和压力float Kl=2.69;//K1:传热系数float Q=29.1;//Q:试井产液float Pwf=4.68;//Pwf:试井流压float fw=0.145;//fw:体积含水率float Z〇=903.87;//Zo:原油密度float Zw=1000;//Zw:水的密度float C〇=2305.61;//Co源油比热float Cw=4318.41;//Cw:地层水比热float Hs=201;//Hs:设计沉没度float Qs=27.6;//Qs:设计排量

float n=6;//n:抽油机冲次float Ab=0.7;//Ab:抽油机栗效float a=9-95,b=3.98; //excel 冋归得 a,b 的值floatwendu(float L);floatniandu(float T);floatliuya();floatpingjunmidu();floatbengya(float Zp);floatbengshen(float Pwfo,floatPs;floatZp);floatbengjing();floatchouyougan(float Lp,floatZp); floatjiaohe();main(void){ float Pwfo,Ps,Zp,Lp,Dp; //设计排量下的井底流压，亲口吸入压力，抽汲液平均密度，泵下入深度，泵的直径Pwfo=liuya();Zp=pingjunmidu();Ps=bengya(Zp);Lp=bengshen(Pwfo,Ps,Zp);Dp=bengjing();chouyougan(Lp,Zp);jiaohe();}//井温分布 floatwendu(float L){float T,q 1=0,W,fz;fz=l*Zw*fw/(l*Zw*fw+(l-fw)*Zo);W=(fz*Cw+(l-fz)*Co)*((Q*fz*Zw+(l-fz)*Q*Zo)/(24*60*60)); T 二(W*O+ql)/Kl*(l-exp(-Kl/W*L))+t0-O*L; return (T);//原油粘温关系 floatniandu(float T){float u;//原油动力粘度u=pow(10,a)/pow(T;b);//粘度公式return(u);{float Ps;

{float Ps;

//井底流压floatliuya(){floatqmax,Pwfo;qmax=Q/(l-0,2*Pwf/Ph-0.8*pow(Pwf/Ph,2)); printf("最大排量:％f m3/d\n",qmax);Pwfo=(-0.2/Ph+pow((pow(0.2/Ph/2)-(4*0.8)/pow(Ph,2)*(Qs/qmax-l))/0.5))/2/0.8*pow(Ph,2);printf(’^计排量下的井底流压：％fMPa\n",Pwfo); return(Pwfo);}//沉没段油、水混合液的平均密度floatpingjunmidu(){floatZp;Zp=Zo*(l-fw)+Zw*fw;return(Zp);}//泵吸入门压力 floatbengya(float Zp)Ps=Zp*g*Hs/pow(10,6);printf(n 泵吸入 U 压力:％f M Pa\nn,Ps);return(Ps);}//下泵深度floatbengshen(float Pwfo,floatPs,floatZp) { floatLp;Lp=H-(Pwfo-Ps)*pow(10,6)/g/Zp; printf("下荥深度:％f m\n",Lp); return (Lp);//确定录径floatbengjing(){floatDp;Dp=pow(Qs/(360*PI*s*n*Ab),0.5)； printf(，?泵径:％f m