东北石油大学石油工程课程设计抽油设计C语言编程.docx

#include stdio.h #include stdlib.h #include math.h #define PI 3.1415926 #define g 9.8 float H=1545; //H:油层深度float Do=0.0889; //Do:油管内径float Dt=0.19; //Dt:套管直径float O=0.0323; //O:地温梯度float t0=79.8; //t0:井底温度float Ph=10.75; //Ph:地层压力float Pb=9.98; //Pb:饱和压力float K1=2.69; //K1:传热系数float Q=29.1; //Q:试井产液float Pwf=4.68;//Pwf:试井流压float fw=0.145; //fw:体积含水率float Zo=903.87; //Zo:原油密度float Zw=1000; //Zw:水的密度float Co=2305.61; //Co:原油比热float Cw=4318.41; //Cw:地层水比热float Hs=201; //Hs:设计沉没度float Qs=27.6; //Qs:设计排量float s=3; //s:抽油机冲程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,q1=0,W,fz; fz=1*Zw*fw/(1*Zw*fw+(1-fw)*Zo); W=(fz*Cw+(1-fz)*Co)*((Q*fz*Zw+(1-fz)*Q*Zo)/(24*60*60)); T=(W*O+q1)/K1*(1-exp(-K1/W*L))+t0-O*L; return (T); } //原油粘温关系floatniandu(float T) { float u; //原油动力粘度u=pow(10,a)/pow(T,b); //粘度公式return(u); } //井底流压floatliuya() { floatqmax,Pwfo; qmax=Q/(1-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-1)),0.5))/2/0.8*pow(Ph,2); printf(设计排量下的井底流压：%fMPa\n,Pwfo); return(Pwfo); } //沉没段油、水混合液的平均密度floatpingjunmidu() { floatZp; Zp=Zo*(1-fw)+Zw*fw; return(Zp); } //泵吸入口压力floatbengya(float Zp) { float Ps; Ps=Zp*g*Hs/pow(10,6); printf(泵吸入口压力:%fMPa\n,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\n,Dp); return (Dp); } //悬点载荷计算及抽油杆强度校核floatchouyougan(floa