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

float s=3; //s:抽油机冲程 #include stdio.h #include stdlib.h #include math-h #define PI 3.1415926 #define g 9.8 float 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\n,Dp); ret