聚驱油田有杆泵抽油系统效率计算方法研究 第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);
}

 << 上一页  [11] [12] [13] 

聚驱油田有杆泵抽油系统效率计算方法研究 第13页下载如图片无法显示或论文不完整，请联系qq752018766