聚驱油田有杆泵抽油系统效率计算方法研究 第13页
、聚驱井的系统效率计算程序
#include "stdio.h"
#include "math.h"
float qc[10]={36,13,73,13,17,14,102};
float H [10] ;
float Nd[10]={11.16,6.69,11.31,6.37,6.49,6.27,8.39} ;
float pt[10]={0.38,0.51,0.36,0.35,0.4,0.47,0.36} ;
float pc[10]={0.9,0.57,0.65,0,0.63,1.22,0.46} ;
float r[10]={0};
float Hd[10]={930.6,719.4,896.2,729,943,686,302.2} ;
float Hp[10]={981.35,924.57,1042.56,923.92,950.29,1000.85,1002};
float Hz [10];
float Hs [10];
float r0=0.86;
float Ng[10];
float L[10]={981,898,1042,900,950,1000,1002};
float d[10]={22,22,22,22,22,22,22} ;
float D[10]={62,62,62,62,62,62,62};
float A[10]={349.62,167.917,1041.61,219.594,662.26,779.337,517.781};
float sd[10]={0.01,0.01,0.02,0.01,0.02,0.02,0.02};
float fd[10]={1.75,1.75,2.01,1.75,1.75,1.75,1.75,1.75,1.75,1.75};
float cn[10]={6,6,3.5,5,5,4,7};
float H1[10];
float fw[10]={0.985,0.929,0.977,0.875,0.713,0.903,0.938};
/* qc-油井实际产液量 */
/* H-有效举液高度 */
/* Nd-电动机功率 */
/* pt-油井油压 */
/* pc-油井套压 */
/* r-油井内混合液体密度 */
/* Hd-油井动液面高度 */
/* Hp-下泵深度 */
/* Hz-折算有效举液高度 */
/* Hs-折算沉没度 */
/* r0-原油密度 */
/* Ng-光杆效率 */
/* L -油管长度 */
/* d -抽油杆直径mm */
/* D -油管直径 mm */
/* cn -冲程 m */
/* p- 沉没压力 */
/* K--气油比 */
/* fw-含水率 */
/* A-示功图面积 */
/* sd-示功图减程比 */
/* fd-示功图力比 */
#define N 10 /*N-井数*/
void changguifa(float n[N]) /* 常规法 */
{int i;
for (i=0;i<N;i++)
{r[i] =1.0*fw[i]+r0*(1-fw[i]);
H[i] =Hd[i]+102*(pt[i]-pc[i])/r[i];
n[i]=qc[i]*H[i]/8812.8/Nd[i];
}
}
void gaijinfa(float n1[N]) /* 改进算法 */
{int i;
for(i=0;i<N;i++)
for(i=0;i<N;i++)
{ r[i]=1.0*fw[i]+r0*(1-fw[i]);
Ng[i]=A[i]*sd[i]*fd[i]*cn[i]/60;
Hz[i]=Hp[i]+100*pt[i]/r[i];
Hs[i]=(Hp[i]-Hd[i])*(r0/r[i])+100*pc[i]/r[i]/10000;
n1[i]=Ng[i]/Nd[i]*Hz[i]*qc[i*86400)/(Ng[i]+Hs[i]*qc[i]*1000/102/86400);
}
}751com.cn
void nizipenjingfa(float n2[N]) /* 拟自喷井法 */
{int i;
float f[N],D1[N]; /*f需要计算*/
float Hn[N]; /* Hn-拟自喷井下有效的举液高度 */
for(i=0;i<N;i++)
{r[i]=fw[i]*1.0+(1-fw[i])*r0 ;
D1[i]=sqrt(D[i]*D[i]-d[i]*d[i]) / 100 ;
f[i]=(qc[i]/(10390.8*D1[i]*D1[i]*D1[i]*D1[i])) *r[i];
Hn[i]=L[i]-(L[i]-Hd[i])*(r0/(r[i]+f[i]))+102*(pt[i]-pc[i])/(r[i]+f[i]) ;
n2[i]=9.8*qc[i]*Hn[i]/Nd[i]/86400;
}
}
void nengliangfenxifa(float n3[N]) /*能量分析法*/
{int i;
float E[N]={0};
for(i=0;i<N;i++)
{r[i]=1.0*fw[i]+0.86*(1-fw[i]);
751com.cn
void ceshifa(float n4[N]) /*测试法*/
{int i;
for(i=0;i<N;i++)
{ r[i]=fw[i]*1.0+(1-fw[i])*r0 ;
H1[i]=Hd[i]+1000*(pt[i]-pc[i])/r[i]/9.8;
n4[i]=H1[i]*qc[i]*r[i]*9.8/Nd[i]/86400;
}
}
main()
{int i;
FILE *fp;
float n[N];
float n1[N];
float n2[N];
float n3[N];
float n4[N] ;
changguifa(n) ;
gaijinfa(n1) ;
nizipenjingfa(n2) ;
nengliangfenxifa(n3) ;
ceshifa(n4);
fp=fopen("e:\\shuiqu.txt","wt");
if(fp==NULL) printf("open error");
for(i=0;i<N;i++)
{fprintf(fp,"n[%d]=%f ",i,n[i]);
fprintf(fp,"n1[%d]=%f ",i,n1[i]);
fprintf(fp,"n2[%d]=%f ",i,n2[i]);
fprintf(fp,"n3[%d]=%f ",i,n3[i]);
fprintf(fp,"n4[%d]=%f ",i,n4[i]);
fprintf(fp,"\n"); }
fclose(fp);
}