ex51.data

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

CC NUMBER OF REGIONS OF AQUIFER WITH MODIFIED Y PERMEABILITY
*---NMOD
        11    
CC
CC
*-----IMIN    IMAX    JMIN    JMAX    KMIN    KMAX    IFACT    FACTX
  7  7 13 13  1  2  1  34000.
  7  7  9  9  2  5  1  34000.
  7  7  5  5  2  4  1  34000.
 14 14 13 13  1  2  1  34000.
 14 14  9  9  2  6  1  34000.
 14 14  5  5  2  4  1  34000.
 11 11  9  9  2  6  1  34000.
 18 18  9  9  2  4  1  34000.
  8 13  3  6  1  6  2   0.70
 15 16  7 11  1  6  2   0.05
  8 13 12 14  2  2  2   0.70
CC
CC NUMBER OF REGIONS OF AQUIFER WITH MODIFIED Z PERMEABILITY
*---NMOD
        11   
CC
CC
*-----IMIN    IMAX    JMIN    JMAX    KMIN    KMAX    IFACT    FACTX
  7  7 13 13  1  2  1  32000.
  7  7  9  9  2  5  1  32000.
  7  7  5  5  2  4  1  32000.
 14 14 13 13  1  2  1  32000.
 14 14  9  9  2  6  1  32000.
 14 14  5  5  2  4  1  32000.
 11 11  9  9  2  6  1  32000.
 18 18  9  9  2  4  1  32000.
  8 13  3  6  1  6  2   0.70
 15 16  7 11  1  6  2   0.05
  8 13 12 14  2  2  2   0.70
CC    
CC CONSTANT CHLORIDE AND CALCIUM CONCENTRATIONS (MEQ/ML)
*----C50       C60
     0.0017  0.0017
CC
CC*******************************************************************
CC                                                                  *
CC    PHYSICAL PROPERTY DATA                                        *
CC                                                                  *
CC*******************************************************************
CC
CC
CC OIL CONC. AT PLAIT POINT FOR TYPE II(+)AND TYPE II(-), CMC
*---- C2PLC   C2PRC  EPSME  IHAND
      0.      1.     .002      1
CC
CC
*----IFGHBN
    0 
CC SLOPE AND INTERCEPT OF BINODAL CURVE AT ZERO, OPT., AND 2XOPT SALINITY
CC FOR ALCOHOL 1
*----HBNS70 HBNC70 HBNS71 HBNC71 HBNS72 HBNC72
     0.1   .13100   .191   0.020   .363   .13100
CC SLOPE AND INTERCEPT OF BINODAL CURVE AT ZERO, OPT., AND 2XOPT SALINITY
CC FOR ALCOHOL 2
*----HBNS80 HBNC80 HBNS81 HBNC81 HBNS82 HBNC82
     0.     0.     0.     0.     0.     0.
CC
CC LOWER AND UPPER EFFECTIVE SALINITY FOR ALCOHOL 1 AND ALCOHOL 2
*----CSEL7  CSEU7  CSEL8  CSEU8
     .1521   0.1521  0.   0.
CC
CC THE CSE SLOPE PARAMETER FOR CALCIUM AND ALCOHOL 1 AND ALCOHOL 2
*----BETA6  BETA7  BETA8
     .8     -2.    0.
CC
CC FLAG FOR ALCOHOL PART. MODEL AND PARTITION COEFFICIENTS
*----IALC  OPSK7O  OPSK7S  OPSK8O  OPSK8S
     0     0.      0.      0.      0.
CC
CC NO. OF ITERATIONS, AND TOLERANCE
*----NALMAX   EPSALC
     20       .0001
CC
CC ALCOHOL 1 PARTITIONING PARAMETERS IF IALC=1
*----AKWC7   AKWS7  AKM7  AK7     PT7
     4.671   1.79   48.   35.31   .222 
CC
CC ALCOHOL 2 PARTITIONING PARAMETERS IF IALC=1
*----AKWC8   AKWS8  AKM8  AK8     PT8
     0.      0.     0.    0.      0
CC
CC 
*-----IFT FLAG  (=0 FOR HEALY AND REED, =1 FOR CHUN HUH).
     1
CC  FOR CHUN HUN FIRST PARAMETER = CONSTANT, SECOND = A CONSTANT
CC INTERFACIAL TENSION PARAMETERS
*---- CHUH   AHUH
     0.22   9.  
CC
CC LOG10 OF OIL/WATER INTERFACIAL TENSION 
*----XIFTW  LOG IFT = 0.68 => IFT = 4.8  DYNES/CM
     0.6800
CC
CC OIL SOLUBILITY IN WATER FLAG
*---- IMASS  ICOR
      1       0
CC
CC TCE SOLUBILITY IN WATER
*----  WSOL(DNAPL SOL IN VOL/VOL.)
      0.00080 
CC
CC CAPILLARY DESATURATION PARAMETERS FOR PHASE 1, 2, AND 3
*----ITRAP   T11        T22        T33
     1       1865.     59074.    364.2 
CC
CC KRJ AND PC (0:IMBIBITION COREY, 1: DRAINAGE COREY), 2:VAN GENUCHTEN
*---- IPERM 
      0 
CC
CC FLAG FOR CONSTANT OR VARIABLE RESIDUAL SATURATIONS
*----ISRW  IPRW  IEW
     0     0      0
CC
CC CONSTANT RES. SATURATION  BY LAYER               
*----S1RWC  S2RWC   S3RWC  
     0.20   0.20   0.20
CC
CC ENDPOINT REL. PERM.                
*----P1RW  P2RW  P3RW 
     0.388  0.429 0.388
CC
CC REL. PERM. EXPONENT                  
*----E1W  E2W  E3W 
     3.3  2.8  3.3      
CC
CC  RES. SATURATION OF PHASES 1,2,AND 3 AT HIGH CAPILLARY NO.
*----S1RC  S2RC  S3RC
     0.0   0.0   0.0
CC
CC ENDPOINT REL. PERM. OF PHASES 1,2,AND 3 AT HIGH CAPILLARY NO.
*----P1RC  P2RC  P3RC
     1.0    1.0  1.0    
CC
CC REL. PERM. EXPONENT OF PHASES 1,2,AND 3 AT HIGH CAPILLARY NO.
*----E13C  E23C  E31C
     3.3   2.8   3.3
CC
CC WATER AND OIL VISCOSITY , RESERVOIR TEMPRATURE
*----VIS1   VIS2  TEMPV
     0.90    0.770    0.0
CC
CC VISCOSITY PARAMETERS
*----ALPHA1 ALPHA2  ALPHA3  ALPHA4  ALPHA5
     3.4    1.       3.      1.     1.  
CC
CC PARAMETERS TO CALCULATE POLYMER VISCOSITY AT ZERO SHEAR RATE
*----AP1     AP2     AP3
      73.    1006.   10809.
CC
CC PARAMETER TO COMPUTE CSEP,MIN. CSEP, AND SLOPE OF LOG VIS. VS. LOG CSEP 
*----BETAP CSE1  SSLOPE
     2.    .01   0     
CC
CC PARAMETER FOR SHEAR RATE DEPENDENCE OF POLYMER VISCOSITY
*----GAMMAC  GAMHF  POWN
     160.    188.   1.843
CC
CC FLAG FOR POLYMER PARTITIONING, PERM. REDUCTION PARAMETERS
*----IPOLYM EPHI3 EPHI4 BRK    CRK
     1      1.    1.0   1000.  0.0186
CC
CC SPECIFIC WEIGHT FOR COMPONENTS 1,2,3,7,AND 8 , AND GRAVITY FLAG
*----DEN1   DEN2  DEN23    DEN3 DEN7   DEN8   IDEN
    .4334  .5981  0.5981 0.4615 0.433 0.00000  2
CC
CC  FLAG FOR UNITS (0:RESERVOIR VOLUMES, 2: SURFACE VOLUMES)
*----ISTB
     0
CC
CC COMPRESSIBILITY FOR VOL. OCCUPYING COMPONENTS 1,2,3,7,AND 8 
*----COMPC(1)  COMPC(2)      COMPC(3)  COMPC(7)  COMPC(8) 
     0.000003  0.000003      0.000003     0.0         0.0
CC
CC CONSTANT OR VARIABLE PC PARAM., WATER-WET OR OIL-WET PC CURVE FLAG 
*----ICPC   IEPC  IOW 
     0      0      0
CC
CC CAPILLARY PRESSURE PARAMETERS, CPC
*-    (=2.9 FOR BORDEN IMBIBITION)
       2.900
CC
CC CAPILLARY PRESSURE PARAMETERS, EPC
*-
     2.0 
CC
CC MOLECULAR DIFFUSIVITY OF KCTH COMPONENT IN PHASE 1 (D(KC),KC=1,N)
*----D(1) D(2) D(3-6)
     0.   0.   0.0   0.0  0.0   0.0 
CC
CC MOLECULAR DIFFUSIVITY OF KCTH COMPONENT IN PHASE 2 (D(KC),KC=1,N)
*----D(1) D(2)  D(3-6)
     0.   0.    0.0   0.0  0.0  0.0 
CC
CC MOLECULAR DIFFUSIVITY OF KCTH COMPONENT IN PHASE 3 (D(KC),KC=1,N)
*----D(1) D(2)   D(3-6)
     0.   0.     0.0   0.0   0.0   0.0  
CC
CC LONGITUDINAL AND TRANSVERSE DISPERSIVITY OF PHASE 1
*----ALPHAL(1)     ALPHAT(1)
     0.010         0.00
CC
CC LONGITUDINAL AND TRANSVERSE DISPERSIVITY OF PHASE 2
*----ALPHAL(2)     ALPHAT(2)
     0.010         0.00
CC
CC LONGITUDINAL AND TRANSVERSE DISPERSIVITY OF PHASE 3
*----ALPHAL(3)     ALPHAT(3)
     0.010         0.00
CC
CC
*-- IADOIL
    0
CC
CC SURFACTANT AND POLYMER ADSORPTION PARAMETERS
*----AD31  AD32  B3D    AD41   AD42  B4D  IADK  IF F  REFK
     1.    .5    1000.  0.7    0.    100.  0    0  0   0.0
CC
CC PARAMETERS FOR CATION EXCHANGE OF CLAY AND SURFACTANT
*----QV     XKC   XKS  EQW
     0.0    .0     0.   419.
CC
CC*******************************************************************
CC                                                                  *
CC    WELL DATA                                                     *
CC                                                                  *
CC*******************************************************************
CC
CC 
CC FLAG FOR PRESSURE CONSTRAINED BOUNDARY
*---- IBOUND  IZONE
      1        0
CC
CC FLAG FOR PRESSURE CONSTRAINED BOUNDARY
*---- IBL     IBR   IBT   WTABL  WTABR                
      1         1    0    0      0