ex21.data

来自「油田化学驱模拟的经典 fortran源代码」· DATA 代码 · 共 654 行 · 第 1/2 页

10468.63   7726.30   8274.36  11501.81   6263.75   6924.69
3840.92   3547.84  10738.49   6498.75
 5324.77   5155.97   4699.81   4438.36   4975.74   4556.38
 2086.02   3623.64   6495.75   6324.55
 3179.49   5744.87   5386.04   6983.89   4997.27   6072.15
 11416.56   8367.33   7140.67  
 3304.28   3515.29   3945.95   6141.03   8669.01  13914.33
 13287.43  10351.76   5533.01   5026.14
 6280.35   6949.74   6172.62   5630.68   2815.24   2416.38
 7370.20  11925.61  15745.65   5190.22
 3132.79   9830.07   7819.78   4433.85  15645.95  25319.40
 12283.62   7389.76   9344.47  10480.85
12122.90   4043.20  10651.00  12279.79   5600.20   8500.72
24024.40  21278.79   6963.28   8900.99
 9839.15   9497.49  12868.06  11922.57  18643.98  10208.87
 4746.14   3670.54   4512.34  
11163.06  16398.55  11032.31  11640.58   5657.36   7362.08
6029.21   6049.42   7783.38  16846.66
 9067.25   6548.27   5295.09   5912.19   8035.48  12795.84
 9439.49   3562.96   2433.96   2700.37
 6230.24   7036.53   7400.94   5436.56   5659.47   5561.34
 7692.90  11101.21  14211.09  14673.87
 4563.25   9016.72  18031.22   4392.95   2905.64   8920.62
 5390.33   2202.28   4075.59   4843.10
19359.18  13485.60   7230.65   2639.01   4853.58   7444.95
9193.50   5537.60   8441.97  
11663.72  11237.70   5500.50   6200.01   6368.65   6642.54
6568.81   7283.87   4665.76   6060.08
 5278.61   4280.67   6301.39   4787.36   6416.86   5519.09
 3735.57  16566.71  13609.30   9149.26
 9719.80  24263.75   7810.49   5463.35  11374.72   5757.10
 7712.98   9819.11   9824.67  12937.24
 9507.59   5979.93   4215.48   6181.56   3677.49   4768.31
 4501.42   8396.36  10921.76  10122.28
19815.61  12731.68   8954.00   6449.72   7234.27  12143.87
14348.91   9929.63  12481.54  
CC 3.3.8
CC CONSTANT Y-PERMEABILITY (MILIDARCY) FOR LAYER K = 1,NZ
*----PERMY(1)             
      8300.0     
CC
CC KV = .5 KX FOR LAYER K = 1,NZ
*----FACTZ(1)             
     0.5
CC
CC FLAG FOR CONSTANT OR VARIABLE DEPTH, PRESSURE, WATER SATURATION
*----IDEPTH  IPRESS  ISWI  ICWI
      0        1       0    -1
CC
CC CONSTANT DEPTH (FT)
*----D111
      0.0 
CC
CC INITIAL PRESSURE (PSIA)
*----PINIT   DEPTH  
      14.7    0.0      
CC
CC INITIAL WATER SATURATION
*----SWI
 1.00 
CC 
CC CONSTANT CHLORIDE AND CALCIUM CONCENTRATIONS (MEQ/ML)
*----C50       C60
    0.009       0.009      
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.0     .0001    0
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    .1006   .191   0.0267    .363    0.1006
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
     0.0198 0.0231   0.     0.
CC
CC THE CSE SLOPE PARAMETER FOR CALCIUM AND ALCOHOL 1 AND ALCOHOL 2
*----BETA6  BETA7  BETA8
     1.0     -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 (FOR TWO ALCOLHOL SYSTEM)
*----NALMAX   EPSALC
     20        0.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
     0
CC
CC INTERFACIAL TENSION PARAMETERS
*----G11  G12     G13   G21      G22      G23
     13.  -14.8   .007  13.0  -14.5  0.010
CC
CC LOG10 OF OIL/WATER INTERFACIAL TENSION 
*----XIFTW
     1.6500
CC
CC FLAG FOR OIL SOLUBILITY IN WATER
*---- IMASS  ICOR
       1      0
CC
CC OIL PART. IN WATER WITH NO SURF., MASS TRANSFER COEF.
*---- WSOL       CNEM2   ISOL
      0.00015    0.0      0
CC
CC CAPILLARY DESATURATION PARAMETERS FOR PHASE 1, 2, AND 3
*----ITRAP  T11    T22     T33
     0       0.      0.      0.
CC
CC FLAG FOR DRAINAGE OR IMBIBITION REL. PERM. CURVES
*----IPERM 
     1 
CC
CC CONSTANT RES. SATURATION OF PHASES 1,2,AND 3 AT LOW CAPILLARY NO.
*----ISRW IPRW IEW                                          
     0    0     0
CC
CC CONSTANT RES. SATURATION OF PHASES 1,2,AND 3 AT LOW CAPILLARY NO.
*----S1RWC  S2RWC  S3RWC                                           
     .240    .000  .240 
CC
CC ENDPOINT REL. PERM. OF PHASES 1,2,AND 3 AT LOW CAPILLARY NO.
*----P1RW  P2RW  P3RW
     1.000  0.650  1.000    
CC
CC REL. PERM. EXPONENT OF PHASES 1,2,AND 3 AT LOW CAPILLARY NO.
*----E1W    E2W   E3W
     2.85  2.70  2.85
CC
CC WATER AND OIL VISCOSITY , RESERVOIR TEMPERATURE
*----VIS1   VIS2  TEMPV
     1.0    0.89     0.
CC
CC VISCOSITY PARAMETERS
*----ALPHA1 ALPHA2  ALPHA3  ALPHA4  ALPHA5
     3.4     1.0     3.0     1.     1.
CC
CC PARAMETERS TO CALCULATE POLYMER VISCOSITY AT ZERO SHEAR RATE
*----AP1     AP2     AP3
     52.     2430.   40000.
CC
CC PARAMETER TO COMPUTE CSEP,MIN. CSEP, AND SLOPE OF LOG VIS. VS. LOG CSEP 
*----BETAP CSE1  SSLOPE (FOR POLYMER PROCESSES)
     2.     .01  .175
CC
CC PARAMETER FOR SHEAR RATE DEPENDENCE OF POLYMER VISCOSITY
*----GAMMAC  GAMHF  POWN
     4.      20.    1.1    
CC
CC FLAG FOR POLYMER PARTITIONING, PERM. REDUCTION PARAMETERS
*----IPOLYM EPHI3 EPHI4 BRK    CRK
     1      1.    1.    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  
   .433    .7036  0.7036  .4980   .346  0.   2 
CC
CC FT3 OR RES. BBLS
*---- 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.0       0.0           0.        0.        0.
CC
CC CONSTANT OR VARIABLE CAPILLARY PRESSURE CURVES, FLAG FOR OIL OR WATER WET
*---- ICPC   IEPC   IOW
      0       0      0
CC
CC CAPILLARY PRESSURE PARAMETERS,CPC
*----CPC 
     1.4
CC
CC CAPILLARY PRESSURE PARAMETERS,EPC
*---- EPC 
      -0.52
CC
CC MOLECULAR DIFFUSIVITY OF KCTH COMPONENT IN PHASE 1 (D(KC),KC=1,N)
*D(1) D(2) D(3) D(4) D(5) D(6) D(7) D(8) D(9) D(10) D(11)                   
 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) D(4) D(5) D(6) D(7) D(8) D(9) D(10) D(11)                    
 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) D(4) D(5) D(6) D(7) D(8) D(9) D(10) D(11)                    
 0.   0.   0.   0.   0.   0.                 
CC
CC LONGITUDINAL AND TRANSVERSE DISPERSIVITY OF PHASE 1 (FEET)
*----ALPHAL(1)     ALPHAT(1)
      0.10          0.03
CC
CC LONGITUDINAL AND TRANSVERSE DISPERSIVITY OF PHASE 2
*----ALPHAL(2)     ALPHAT(2)
      0.10          0.03
CC
CC LONGITUDINAL AND TRANSVERSE DISPERSIVITY OF PHASE 3
*----ALPHAL(3)     ALPHAT(3)
     0.10           0.03
CC
CC
*----IADSO
      0
CC
CC SURFACTANT AND POLYMER ADSORPTION PARAMETERS
*----AD31  AD32  B3D    AD41   AD42   B4D   IADK  IADS1  FADS REFK
     1.0    0.5  1000.  .7      0.    100.    0    0      0    100
CC
CC PARAMETERS FOR CATION EXCHANGE OF CLAY AND SURFACTANT
*----QV     XKC   XKS  EQW
     0.0    0.0   0.0   1.  
CC
CC**********************************************************************
CC                                                                     *
CC VII    WELL DATA                                                    *
CC                                                                     *
CC**********************************************************************
CC
CC
CC  FLAG FOR BOUNDARIES
*---- IBOUND  IZONE
        1      0
CC
CC  LEFT BOUNDARY, RIGHT BOUNDARY FLAGS
*--- IBL   IBR  IBT  WTABL  WTABR
     1     1    0     0      0
CC
CC PRESSURE, CONCENTRATION AT THE LEFT BOUNDARY
*---- PBL, C1BL, C5BL, C6BL
     15.028  1.0  0.009 0.009
CC
CC PRESSURE, CONCENTRATION AT THE RIGHT BOUNDARY
*---- PBR, C1BR, C5BR, C6BR
     14.7  1.0  0.009 0.009
CC
CC TOTAL NO. OF WELLS, WELL RADIUS MODEL FLAG
*----NWELL  IRO  ITSTEP  NWREL
      3     2     1       3
CC REMEDIATION FLUID (WATER OR SURF.) INJECTOR, RATE CONTROLLED:
CC WELL LOCATIONS, FLAG FOR SPECIFYING WELL TYPE, WELL RADIUS, SKIN
*----IDW  IW   JW   IFLAG    RW     SWELL  IDIR  KFIRST  KLAST IPRF
     1    16    1     1     1.000     0.   3      1      24      0
CC
CC NAME OF THE WELL
*---- WELNAM
W1 REMED.INJ.
CC
CC MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHECK  PWFMIN   PWFMAX  QTMIN   QTMAX
      0      0.0     5000.   0.0     10000.
CC INJECTOR TO MODEL OIL SPILLAGE, RATE CONTROLLED :
CC WELL LOCATIONS, FLAG FOR SPECIFYING WELL TYPE, WELL RADIUS, SKIN
*----IDW  IW   JW   IFLAG    RW     SWELL  IDIR  KFIRST  KLAST IPRF
     2    25    1     1      1.000     0.   3      1     1       0
CC
CC NAME OF THE WELL
*---- WELNAM
W2 PCE SPILLER
CC
CC MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHECK  PWFMIN   PWFMAX  QTMIN   QTMAX
      0      0.0     5000.   0.0     10000.
CC  REMEDIATION PRODUCER, RATE CONTROLLED:
CC WELL LOCATIONS, FLAG FOR SPECIFYING WELL TYPE, WELL RADIUS, SKIN
*----IDW  IW   JW   IFLAG    RW     SWELL  IDIR  KFIRST  KLAST  IPRF
     3    41    1     4      1.000     0.   3      1     24       0
CC
CC NAME OF THE WELL
*---- WELNAM
W3 REMED.PROD.
CC
CC MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHECK PWFMIN   PWFMAX  QTMIN   QTMAX
      0     0.0     5000.   0.0     10000.
CC  W1 REMED.INJ. IFLAG=1       
CC INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----QI(M,L)  C(M,KC,L),                            
 1   000.     1.0     0. 0.      0. 0.009 0.009 
 1   0.       0.      0. 0.      0. 0.    0. 
 1   0.       0.      0. 0.      0. 0.    0.
CC  W2 PCE SPILLER IFLAG=1
CC INJ. RATE AND INJ. COMP. FOR RATE CONS. WELLS FOR EACH PHASE (L=1,3)
*----QI(M,L)  C(M,KC,L), M = 2  (OIL SPILLER)
 2    0.00      0.  0. 0. 0. 0. 0.              
 2    33.333    0.  1. 0. 0. 0. 0.            
 2    0.        0.  0. 0. 0. 0. 0.            
CC  W3 REMED.PROD.IFLAG=2 
CC MAXIMUM RATE OF PRODUCTION (NEG. FOR A PRODUCER)
*----ID  Q(M)                          
     3   0.00    
CC
CC CUM. IN.PRO:MAPS. PRO:WELL   .SUM      .HIS1     .CONC   RESTART      
*----TINJ    CUMPR1   CUMHI1     WRHPV     WRPRF   RSTC
     30.0    15.0     15.0       5.0       5.0     15.0 
CC
CC FOR IMES=3 ,THE INI. TIME STEP,CONC. TOLERANCE,MIN. AND MAX. TIME STEP
*----DT    DCLIM                              CNMAX  CNMIN
     0.03  0.003 0.003 0.003 0.003 0.003 0.003  0.2   0.04