ex01.data

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

CC*******************************************************************
CC                                                                  *
CC     BRIEF DESCRIPTION OF DATA SET : UTCHEM (VERSION 9.0 )        *
CC                                                                  *
CC*******************************************************************
CC                                                                  *
CC  SURFACTANT/POLYMER FLOOD TEST (EX01), 11X11X2                   *
CC                                                                  *
CC  LENGTH (FT) : 250               PROCESS : SURFACTANT/POLYMER    *
CC  THICKNESS (FT) : 10             INJ. RATE (FT3/DAY) : 112.3     *
CC  WIDTH (FT) : 250                COORDINATES : CARTESIAN         *
CC  POROSITY : 0.20                                                 *
CC  GRID BLOCKS : 11X11X2                                           *
CC  DATE : 04/05/2000                                               *
CC                                                                  *
CC*******************************************************************
CC
CC*******************************************************************
CC                                                                  *
CC    RESERVOIR DESCRIPTION                                         *
CC                                                                  *
CC*******************************************************************
CC  
CC
*----RUNNO
EX01
CC  
CC
*----HEADER
EX01
TESTING  UTCHEM VERSION 9  ON  DEC 500au (MARATHON)
SURFACTANT/POLYMER, SALINITY GRADIENT RUN (.4 -> .3 -> .1). 
CC
CC SIMULATION FLAGS
*---- IMODE IMES IDISPC ICWM ICAP IREACT IBIO ICOORD ITREAC ITC IGAS IENG 
        1    2    3      0    0     0      0    1     0      0    0   0
CC
CC NO. OF GRIDBLOCKS,FLAG SPECIFIES CONSTANT OR VARIABLE GRID SIZE, UNIT
*----NX   NY  NZ  IDXYZ   IUNIT
     11   11  2   0        0          
CC
CC  CONSTANT GRID BLOCK SIZE IN X, Y, AND Z
*----DX1       DY1      DZ1
     22.727    22.727  5.0          
CC
CC TOTAL NO. OF COMPONENTS, NO. OF TRACERS, NO. OF GEL COMPONENTS
*----N   NO  NTW NTA  NGC NG  NOTH
     11   0  3   0   0    0    0
CC
CC  NAME OF SPECIES
*---- SPNAME(I) FOR I=1,N 
WATER
OIL
SURF.
POLYMER
ANION
CALCIUM
ALCOHOL
GAS(NO)
TRACER 1
TRACER 2
TRACER 3
CC
CC FLAG INDICATING IF THE COMPONENT IS INCLUDED IN CALCULATIONS OR NOT
*----ICF(KC) FOR KC=1,N
     1  1  1  1  1  1  1  0  1  1  1 
CC
CC*******************************************************************
CC                                                                  *
CC    OUTPUT OPTIONS                                                *
CC                                                                  *
CC*******************************************************************
CC  
CC
CC FLAG FOR PV OR DAYS TO PRINT OR TO STOP THE RUN
*----ICUMTM  ISTOP  IOUTGMS
       0         0    0
CC
CC FLAG INDICATING IF THE PROFILE OF KCTH COMPONENT SHOULD BE WRITTEN
*----IPRFLG(KC),KC=1,N
     1  1  1  1  0  0  0  0  0  0  0
CC
CC FLAG FOR individual map files
*----IPPRES IPSAT IPCTOT IPBIO IPCAP IPGEL IPALK IPTEMP IPOBS 
      1      1      1      0      0     0    0     0       0 
CC
CC FLAG for individual output map files
*----ICKL IVIS IPER ICNM ICSE ihystp  ifoamp  inoneq
      1    1    1    1    1   0        0       0 
CC
CC FLAG  for variables to PROF output file
*----Iads IVel Irkf Iphse 
      0    0    0    0 
CC
CC*******************************************************************
CC                                                                  *
CC    RESERVOIR PROPERTIES                                          *
CC                                                                  *
CC*******************************************************************
CC  
CC
CC MAX. SIMULATION TIME ( DAYS)
*---- TMAX
       1500
CC
CC ROCK COMPRESSIBILITY (1/PSI), STAND. PRESSURE(PSIA)
*----COMPR   PSTAND
      0.      0.
CC
CC FLAGS INDICATING CONSTANT OR VARIABLE POROSITY, X,Y,AND Z PERMEABILITY
*----IPOR1 IPERMX IPERMY IPERMZ  IMOD
       0      1     1      0      0
CC
CC CONSTANT POROSITY
*----PORC1
     .20             
CC
CC CONSTANT X-PERMEABILITY (MILIDARCY) FOR LAYER K = 1,NZ
*----PERMX(1)  PERMX(2)   
     500.      100.
CC
CC CONSTANT Y-PERMEABILITY (MILIDARCY) FOR LAYER K = 1,NZ
*----PERMY(1)  PERMY(2)
     500.      100.
CC
CC CONSTANT Z-PERMEABILITY
*----PERMZC (MILIDARCY)
     50.
CC
CC FLAG FOR CONSTANT OR VARIABLE DEPTH, PRESSURE, WATER SATURATION
*----IDEPTH  IPRESS  ISWI  ICWI
      0        0       0   -1 
CC
CC CONSTANT DEPTH (FT)
*----D111
      0.
CC
CC CONSTANT PRESSURE (PSIA)
*----PRESS1
     1.
CC
CC CONSTANT INITIAL WATER SATURATION
*----SWI
     .65
cc
cc
*--- c50    c60 
   0.4     0.003
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.     .0001    0
CC
CC FLAG INDICATING TYPE OF PHASE BEHAVIOR PARAMETERS
*---- 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.131   .1    .191   .026   .363   .028
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
     .177   .344   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
     1     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 MODEL FLAG
      0
CC
CC INTERFACIAL TENSION PARAMETERS
*----G11  G12     G13   G21   G22    G23
     13.  -14.8   .007  13.   -14.5  .010
CC
CC LOG10 OF OIL/WATER INTERFACIAL TENSION 
*----XIFTW
     1.3
CC
CC MASS TRANSFER FLAG
*----IMASS  icor
     0      0
CC
CC CAPILLARY DESATURATION PARAMETERS FOR PHASE 1, 2, AND 3
*----ITRAP   T11        T22        T33
     1       1865.      59074      364.2 
CC
CC RELATIVE PERM. FLAG (0:IMBIBITION COREY,1:FIRST DRAINAGE COREY)
*----IPERM
     0 
CC
CC FLAG FOR CONSTANT OR VARIABLE REL. PERM. PARAMETERS
*----ISRW  IPRW  IEW
     0      0    0
CC
CC CONSTANT RES. SATURATION OF PHASES 1,2,AND 3 AT LOW CAPILLARY NO.
*----S1RWC  S2RWC  S3RWC
     .37    .35    .37
CC
CC CONSTANT ENDPOINT REL. PERM. OF PHASES 1,2,AND 3 AT LOW CAPILLARY NO.
*----P1RW  P2RW  P3RW
     .11   .95   .11
CC
CC CONSTANT REL. PERM. EXPONENT OF PHASES 1,2,AND 3 AT LOW CAPILLARY NO.
*----E1W  E2W  E3W
     1.   2.16  1.
CC
CC  RES. SATURATION OF PHASES 1,2,AND 3 AT HIGH CAPILLARY NO.
*----S1RC  S2RC  S3RC
     .0    .0    .0
CC
CC ENDPOINT REL. PERM. OF PHASES 1,2,AND 3 AT HIGH CAPILLARY NO.
*----P1RC  P2RC  P3RC
     1.    1.    1.
CC
CC REL. PERM. EXPONENT OF PHASES 1,2,AND 3 AT HIGH CAPILLARY NO.
*----E13C  E23C  E31C
     1.     2.16   1.
CC
CC WATER AND OIL VISCOSITY , RESERVOIR TEMPERATURE
*----VIS1   VIS2  TEMPV
     0.86   4.   0.
CC
CC VISCOSITY PARAMETERS
*----ALPHA1 ALPHA2  ALPHA3  ALPHA4  ALPHA5
     2.5     2.3      10.    1.      1.
CC
CC PARAMETERS TO CALCULATE POLYMER VISCOSITY AT ZERO SHEAR RATE
*----AP1     AP2     AP3
     81.     2700.   2500.
CC
CC PARAMETER TO COMPUTE CSEP,MIN. CSEP, AND SLOPE OF LOG VIS. VS. LOG CSEP 
*----BETAP CSE1  SSLOPE
     10.    .01   .17
CC
CC PARAMETER FOR SHEAR RATE DEPENDENCE OF POLYMER VISCOSITY
*----GAMMAC  GAMHF  POWN
     20.      10.    1.8
CC
CC FLAG FOR POLYMER PARTITIONING, PERM. REDUCTION PARAMETERS
*----IPOLYM EPHI3 EPHI4 BRK    CRK
     1      1.    0.8    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  .368   0.368  .42  .346  0.  2 
CC
CC  FLAG FOR CHOICE OF UNITS ( 0:BOTTOMHOLE CONDITION , 1: STOCK TANK)
*-----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.
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 PARAMETER, CPC0 
*----CPC0 
     9. 
CC
CC CAPILLARY PRESSURE PARAMETER, EPC0 
*---- EPC0
      2.
CC
CC MOLECULAR DIFFUSION COEF. 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.   0.   0.   0.   0.    0.
CC
CC MOLECULAR DIFFUSION COEF. 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.   0.   0.   0.    0.     0.
CC
CC MOLECULAR DIFFUSION COEF. 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.   0.   0.   0.   0.     0.
CC
CC LONGITUDINAL AND TRANSVERSE DISPERSIVITY OF PHASE 1
*----ALPHAL(1)     ALPHAT(1)
     12.          .4
CC
CC LONGITUDINAL AND TRANSVERSE DISPERSIVITY OF PHASE 2
*----ALPHAL(2)     ALPHAT(2)
     12.           .4
CC
CC LONGITUDINAL AND TRANSVERSE DISPERSIVITY OF PHASE 3
*----ALPHAL(3)     ALPHAT(3)
     12.           .4
CC
CC FLAG TO SPECIFY ORGANIC ADSORPTION CALCULATION
*----IADSO
      0
CC
CC SURFACTANT AND POLYMER ADSORPTION PARAMETERS
*----AD31  AD32  B3D    AD41   AD42  B4D   IADK, IADS1, FADS REFK,md
     1.5    .5    1000.  0.7    0.    100.  0      0     0    500
CC
CC PARAMETERS FOR CATION EXCHANGE OF CLAY AND SURFACTANT
*----QV     XKC   XKS  EQW
     0.044  .25   .2   419.
CC
CC TRACER PARTITIONING COEFFICIENT (TK(IT),IT=1,NT)
*---- TK(1)  TK(2)  TK(3)
      .0     1.     0.
CC
CC SALINITY DEPENDENCE PART. COEFF., ref. salinity
*---- TKS(1)  TKS(2)  TKS(3)   C5INI 
      0.      0.       0.      0.4
CC
CC RADIACTIVE DECAY COEFFICIENT (RDC(IT),IT=1,NT)
*---- RDC(1)  RDC(2)   RDC(3)
      .0      .0       .0
CC
CC TRACER RETARDATION COEFFICIENT (RET(IT),IT=1,NT)
*---- RET(1)   RET(2)  RET(3)
      .0       0.      0.
CC
CC*******************************************************************
CC                                                                  *
CC    WELL DATA                                                     *
CC                                                                  *
CC*******************************************************************
CC
CC
CC FLAG FOR RIGHT AND LEFT BOUNDARY
*---- IBOUND izone
      0       0
CC  
CC TOTAL NUMBER OF WELLS, WELL RADIUS FLAG, FLAG FOR TIME OR COURANT NO.
*----NWELL   IRO   ITSTEP  NWREL
      2      2      1      2
CC
CC WELL ID,LOCATIONS,AND FLAG FOR SPECIFYING WELL TYPE, WELL RADIUS, SKIN
*----IDW   IW    JW    IFLAG    RW     SWELL  IDIR   IFIRST  ILAST  IPRF
      1    1     1      1       .5     0.      3      1        2    0
CC
CC WELL NAME
*---- WELNAM
INJECTOR
CC
CC ICHEK, MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHEK PWFMIN   PWFMAX  QTMIN   QTMAX
      2    0.0     5000.   0.0     1000.
CC
CC WELL ID, LOCATION, AND FLAG FOR SPECIFYING WELL TYPE, WELL RADIUS, SKIN
*----IDW  IW   JW   IFLAG    RW     SWELL  IDIR  IFIRST   ILAST    IPRF
     2    11   11   2       .5       0.     3     1         2       0
CC
CC WELL NAME
*---- WELNAM
PRODUCER
CC
CC ICHEK , MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHEK PWFMIN   PWFMAX  QTMIN   QTMAX
      2    0.0     5000.   0.0     50000.
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)  
     1    112.3     .94  0.  0.03  .05  0.3  0.001  0.03  0.  1. 1.  1.        
     1     0.        0.   0.  0.    0.   0.   0.     0.    0.  0. 0.  0. 
     1     0.        0.   0.  0.    0.   0.   0.     0.    0.  0. 0.  0.   
CC
CC ID, BOTTOM HOLE PRESSURE FOR PRESSURE CONSTRAINT WELL (IFLAG=2 OR 3)
*----ID   PWF
     2    1.
CC
CC CUM. INJ. TIME , AND INTERVALS (PV OR DAY) FOR WRITING TO OUTPUT FILES
*----TINJ    CUMPR1   cumh1  WRHPV   WRPRF      RSTC
     111.3   50        50     10     5      0.25
CC
CC FOR IMES=2 ,THE INI. TIME STEP,CONC. TOLERANCE,MAX.,MIN. COURANT NUMBERS
*----DT      DCLIM   CNMAX   CNMIN    
     0.05    0.003   0.2     0.01
cc
cc
*--- ibmod
      0
CC
CC IRO, ITIME, NEW FLAGS FOR ALL THE WELLS
*---- IRO ITSTEP IFLAG
       2   1     1  2
CC  
CC NUMBER OF WELLS CHANGES IN LOCATION OR SKIN OR PWF
*----NWEL1
      0
CC  
CC NUMBER OF WELLS WITH RATE CHANGES, ID
*----NWEL1   ID
     1        1
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)  
     1    112.3     1.  0.  0.  0.05  0.1  0.001  0.  0.  0. 0. 0.       
     1    0.        0.  0.  0.  0.    0.   0.     0.  0.  0. 0. 0.
     1    0.        0.  0.  0.  0.    0.   0.     0.  0.  0. 0. 0.
CC
CC CUM. INJ. TIME , AND INTERVALS (PV) FOR WRITING TO OUTPUT FILES
*----TINJ    CUMPR1 cumh   WRHPV     WRPRF   RSTC
     3300    360     360     30       360    730 
CC
CC FOR IMES=2 ,THE INI. TIME STEP,CONC. TOLERANCE,MAX.,MIN. TIME STEPS
*----DT      DCLIM   CNMAX   CNMIN    
     0.05    0.003   0.2     0.02