ex51.data

油田化学驱模拟的经典 fortran源代码

CC 
CC PBL=PRESSURE AT THE CENTER OF THE TOP LAYER AT THE LEFT BOUNDARY
*---- PBL(PSI)   C1BL(C1)    C5BL     C6BL 
      16.107     1.00       0.0017  0.0017
CC
CC PBR=PRESSURE AT THE CENTER OF THE TOP LAYER AT THE RIGHT BOUNDARY
*---- PBR(PSI)   C1BR(C1)    C5BR     C6BR 
      16.107     1.00       0.0017  0.0017
CC  
CC TOTAL NUMBER OF WELLS, WELL RADIUS FLAG, TIME STEP CONTROL FLAG
*----NWELL   IRO   ITIME   NWREL
     8      2     1          8 
CC  IFLAG =1:RI,2:PP,3:PI,4:RP 
CC WELL ID,LOCATIONS,AND FLAG FOR SPECIFYING WELL TYPE, WELL RADIUS, SKIN
*----IDW   IW    JW    IFLAG    RW     SWELL  IDIR   IFIRST  ILAST  IPRF
      1    7     13      4     0.170    0.      3      2        2     0  
CC
CC WELL NAME
*---- WELNAM
 U21(EXT-1)  
CC
CC ICHEK  MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHEK   PWFMIN   PWFMAX  QTMIN   QTMAX
      0      0.0     50000.   0.0    10000.
CC  IFLAG =1:RI,2:PP,3:PI,4:RP 
CC WELL ID,LOCATIONS,AND FLAG FOR SPECIFYING WELL TYPE, WELL RADIUS, SKIN
*----IDW   IW    JW    IFLAG    RW     SWELL  IDIR   IFIRST  ILAST  IPRF
      2     7     9      4     0.170    0.      3      5        5     0  
CC
CC WELL NAME
*---- WELNAM
 SB1(EXT-2)   
CC
CC ICHEK  MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHEK   PWFMIN   PWFMAX  QTMIN   QTMAX
      0      0.0     50000.   0.0    10000.
CC  IFLAG =1:RI,2:PP,3:PI,4:RP 
CC WELL ID,LOCATIONS,AND FLAG FOR SPECIFYING WELL TYPE, WELL RADIUS, SKIN
*----IDW   IW    JW    IFLAG    RW     SWELL  IDIR   IFIRST  ILAST  IPRF
      3     7     5      4     0.170    0.      3      4        4     0  
CC
CC WELL NAME
*---- WELNAM
 SB5(EXT-3)    
CC
CC ICHEK  MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHEK   PWFMIN   PWFMAX  QTMIN   QTMAX
      0      0.0     50000.   0.0    10000.
CC  IFLAG =1:RI,2:PP,3:PI,4:RP 
CC WELL ID,LOCATIONS,AND FLAG FOR SPECIFYING WELL TYPE, WELL RADIUS, SKIN
*----IDW   IW    JW    IFLAG    RW     SWELL  IDIR   IFIRST  ILAST  IPRF
      4    14     13    1      0.170    0.      3      2        2   0
CC
CC WELL NAME
*---- WELNAM
 SB3(INJ-1)    
CC
CC ICHEK  MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHEK   PWFMIN   PWFMAX  QTMIN   QTMAX
      0      0.0     50000.   0.0    10000.
CC  IFLAG =1:RI,2:PP,3:PI,4:RP 
CC WELL ID,LOCATIONS,AND FLAG FOR SPECIFYING WELL TYPE, WELL RADIUS, SKIN
*----IDW   IW    JW    IFLAG    RW     SWELL  IDIR   IFIRST  ILAST  IPRF
      5    14     9     1      0.170    0.      3      6        6   0
CC
CC WELL NAME
*---- WELNAM
 SB2(INJ-2)   
CC
CC ICHEK  MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHEK   PWFMIN   PWFMAX  QTMIN   QTMAX
      0      0.0     50000.   0.0    10000.
CC  IFLAG =1:RI,2:PP,3:PI,4:RP 
CC WELL ID,LOCATIONS,AND FLAG FOR SPECIFYING WELL TYPE, WELL RADIUS, SKIN
*----IDW   IW    JW    IFLAG    RW     SWELL  IDIR   IFIRST  ILAST  IPRF
      6     14    5     1      0.170    0.      3      4        4   0
CC
CC WELL NAME
*---- WELNAM
 SB4(INJ-3)
CC
CC ICHEK  MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHEK   PWFMIN   PWFMAX  QTMIN   QTMAX
      0      0.0     50000.   0.0    10000.
CC  IFLAG =1:RI,2:PP,3:PI,4:RP 
CC WELL ID,LOCATIONS,AND FLAG FOR SPECIFYING WELL TYPE, WELL RADIUS, SKIN
*----IDW   IW   JW IFLAG    RW     SWELL  IDIR   ZFIRST  ZLAST  IPRF 
      7    11   9    4     0.170    0.      3      6        6     0  
CC
CC WELL NAME
*---- WELNAM
 SB6(M1)      
CC
CC ICHEK  MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHEK   PWFMIN   PWFMAX  QTMIN   QTMAX
      0      0.0     50000.   0.0    10000.
CC  IFLAG =1:RI,2:PP,3:PI,4:RP 
CC WELL ID,LOCATIONS,AND FLAG FOR SPECIFYING WELL TYPE, WELL RADIUS, SKIN
*----IDW   IW   JW IFLAG    RW     SWELL  IDIR   ZFIRST  ZLAST  IPRF 
      8    18   9    1     0.170    0.      3      4        4     0  
CC
CC WELL NAME
*---- WELNAM
 SB8(HC)
CC
CC ICHEK  MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHEK   PWFMIN   PWFMAX  QTMIN   QTMAX
      0      0.0     50000.   0.0    10000.
CC
CC ID,INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----ID,QI(M,L)  
     1    -712.  
CC
CC ID,INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----ID,QI(M,L)  
     2    -712.  
CC
CC ID,INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----ID,QI(M,L) 
     3    -712.  
CC
CC ID,INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----ID,QI(M,L)  C(M,KC,L)  
    4   481.   1.00 0. 0.00 0.00 0.1249 0.0051 
    4   0.0    0.    0.0 0.0  0.0  0.0    0.0
    4   0.0    0.    0.0 0.0  0.0  0.0    0.0
CC
CC ID,INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----ID,QI(M,L)  C(M,KC,L)  
    5   481.   1.00 0. 0.00 0.00 0.1249 0.0051
    5   0.0    0.    0.0 0.0  0.0  0.0    0.0
    5   0.0    0.    0.0 0.0  0.0  0.0    0.0
CC
CC ID,INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----ID,QI(M,L)  C(M,KC,L)  
     6   481.   1.00 0. 0.00 0.00 0.1249 0.0051
     6   0.0    0.    0.0 0.0  0.0  0.0    0.0
     6   0.0    0.    0.0 0.0  0.0  0.0    0.0
CC
CC ID,INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----ID,QI(M,L)  C(M,KC,L)  
     7  -0.00 
CC
CC ID,INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----ID,QI(M,L)  C(M,KC,L)  
     8   1348.   1.00 0. 0.00 0.00 0.0051 0.0051 0. 0. 5*0.0 
     8   0.0    0.    0.0 0.0  0.0  0.0    0.0
     8   0.0    0.    0.0 0.0  0.0  0.0    0.0
CC  DAYS    MAPS     PROD.DATA  ELL.DATA  MAPS    RESTART
CC          PROFILE  PROFILE    HISTORIES  CONC   
*----TINJ   CUMPR1   CUMHI1     WRHPV      WRPRF    RSTC
      1     1.0      1          0.1         0.999   200.  
CC
CC FOR IMES=3 ,THE INI. TIME STEP,CONC. TOLERANCE,MIN. AND MAX. TIME STEP
*----DT      DCLIM      CNMAX   CNMIN 
     1E-8   0.003       0.100  0.005
CC  ********************** surfactant flush *********************  
CC
*--- IBMOD
     0
CC
CC WELL RADIUS FLAG, TIME STEP CONTROL FLAG, FLAG FOR WELLS
*----IRO   ITIME   IFLAG(M),M=1,NWELL
      2      1      4  4  4  1  1  1  4  1
CC
CC  NO. OF WELLS WITH CHANGES IN LOCATION ETC
*--- NWEL1
     0
CC
CC  NO. OF WELLS WITH CHANGES IN RATE OR CONC.
*---- NWEL2  ID
      3      4  5  6  
CC
CC ID,INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----ID,QI(M,L)  C(M,KC,L)   0.1249 SALINITY ADJUSTED FOR 0.9245 WATER FRACT)
     4   481.0  0.9245    0.  0.0755 0.00 0.1155 0.0051 
     4   0.0       0.     0.0  0.0   0.0   0.    0.0
     4   0.0       0.     0.0  0.0   0.0   0.    0.0
CC
CC ID,INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----ID,QI(M,L)  C(M,KC,L)  
     5   481.0  0.9245 0. 0.0755 0.00 0.1155 0.0051
     5   0.0       0.     0.0  0.0   0.0   0.    0.0
     5   0.0       0.     0.0  0.0   0.0   0.    0.0
CC
CC ID,INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----ID,QI(M,L)  C(M,KC,L)  
     6   481.0  0.9245 0. 0.0755 0.00 0.1155 0.0051 
     6   0.0       0.     0.0  0.0   0.0   0.    0.0
     6   0.0       0.     0.0  0.0   0.0   0.    0.0
CC  DAYS    MAPS     PROD.DATA WELL.DATA  MAPS    RESTART
CC          PROFILE  PROFILE   HISTORIES  CONC   
*----TINJ   CUMPR1   CUMHI1    WRHPV      WRPRF   RSTC
      4.4    1.7      5       0.050       1.699    17  
CC
CC FOR IMES=3 ,THE INI. TIME STEP,CONC. TOLERANCE,MIN. AND MAX. TIME STEP
*----DT      DCLIM      CNMAX   CNMIN 
     1E-7   0.003       0.050     0.001
CC  ******************PII POST SURF.WATERFLUSH***************  
CC
*--- IBMOD
     0
cc
CC WELL RADIUS FLAG, TIME STEP CONTROL FLAG, FLAG FOR WELLS
*----IRO   ITIME   IFLAG(M),M=1,NWELL
      2      1      4  4  4  1  1  1  4  1 
CC
CC  NO. OF WELLS WITH CHANGES IN LOCATION ETC
*--- NWEL1
     0
CC
CC  NO. OF WELLS WITH CHANGES IN RATE OR CONC.
*---- NWEL2  ID
      3      4  5  6
CC
CC ID,INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----ID,QI(M,L)  C(M,KC,L)  
     4   481.0  1.00     0.   0.0   0.0  0.0051 0.0051 
     4   0.0    0.       0.0  0.0   0.0   0.    0.0
     4   0.0    0.       0.0  0.0   0.0   0.    0.0
CC
CC ID,INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----ID,QI(M,L)  C(M,KC,L)  
     5   481.0  1.00   0.   0.0   0.0 0.0051 0.0051
     5   0.0    0.     0.0  0.0   0.0   0.    0.0
     5   0.0    0.     0.0  0.0   0.0   0.    0.0
CC
CC ID,INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----ID,QI(M,L)  C(M,KC,L)  
     6   481.0  1.00      0.  0.0    0.0   0.0051 0.0051 
     6   0.0       0.     0.0  0.0   0.0   0.     0.0
     6   0.0       0.     0.0  0.0   0.0   0.     0.0
CC  DAYS    MAPS     PROD.DATA WELL.DATA  MAPS    RESTART
CC          PROFILE  PROFILE    HISTORIES  CONC   
*----TINJ   CUMPR1   CUMHI1     WRHPV      WRPRF   RSTC
      22.4    9.0     9.0       0.10        8.999    22  
CC
CC FOR IMES=3 ,THE INI. TIME STEP,CONC. TOLERANCE,MIN. AND MAX. TIME STEP
*----DT      DCLIM      CNMAX   CNMIN 
     1E-7   0.003       0.1     0.005