ex21.data

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

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