ex45.data

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

CC*******************************************************************
CC                                                                  *
CC    BRIEF DESCRIPTION OF DATA SET : UTCHEM (VERSION 9.0)          *
CC                                                                  *
CC*******************************************************************
CC                                                                  *
CC  GEOCHEMISTRY 1D TEST RUN                                        *
CC                                                                  *
CC  LENGTH (FT) : 100               PROCESS : SINGLE PHASE WATER    * 
CC  THICKNESS (FT) : 10.            INJ. RATE (FT3/DAY) : 10        *
CC  WIDTH (FT) : 10.                COORDINATES : CARTISIAN         *
CC  POROSITY : 0.25                                                 *
CC  GRID BLOCKS : 100x1x1                                           *
CC  DATE : 04/24/2000                                                *
CC                                                                  *
CC*******************************************************************
CC
CC*******************************************************************
CC                                                                  *
CC    RESERVOIR DESCRIPTION                                         *
CC                                                                  *
CC*******************************************************************
CC  
CC
*----RUNNO
 ex45
CC  
CC
*----HEADER
run ex45 using utchem v9.0 
GEOCHEMISTRY - PB AND CR
SINGLE PHASE FLOW
CC
CC SIMULATION FLAGS
*---- IMODE IMES IDISPC ICWM ICAP IREACT IBIO ICOORD ITREAC ITC IGAS ieng
        1    2    3      0    0     2      0   1     0      1    0   0
CC
CC NUMBER OF GRID BLOCKS AND FLAG SPECIFIES CONSTANT OR VARIABLE GRID SIZE
*----NX   NY  NZ  IDXYZ  IUNIT
    100    1    1  0       0          
CC
CC  CONSTANT GRID BLOCK SIZE IN X, Y, AND Z
*----DX1       DY1      DZ1
     1.        10.       10.
CC
CC TOTAL NO. OF COMPONENTS, NO. OF TRACERS, NO. OF GEL COMPONENTS
*---- N   NO NTW   NTA   NGC   NG  NOTH
      19  0  0      0    11    0   0
cc
cc
*--- species name
water
oil
surf.
polymer
anion
calcium
comp.7
comp.8
gc-1
gc-2
gc-3
gc-4
gc-5
gc-6
gc-7
gc-8
gc-9
gc-10
gc-11
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  0  0  1  1  1  1  1  1  1 1 1 1 1 1 1 1 1 
CC
CC*******************************************************************
CC                                                                  *
CC    OUTPUT OPTIONS                                                *
CC                                                                  *
CC*******************************************************************
CC  
CC
CC  FLAG FOR PV OR DAYS for output and 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  1  1  0  0  1  1  1  1  1  1  1 1 1 1 1 1 1 1 1 1 1 
CC
CC FLAG FOR PRES,SAT.,TOTAL CONC.,TRACER CONC.,CAP.,GEL, ALKALINE PROFILES
*----IPPRES IPSAT IPCTOT IPBIO IPCAP IPGEL IPALK itemp ipobs 
      0      0      1      0      0     0    1      0   0  
CC
CC FLAG FOR WRITING SEVERAL PROPERTIES 
*----ICKL IVIS IPER ICNM ICSE  ihystp   ifoamp  inoneq
      0     0    0    0    0    0         0       0 
CC
CC FLAG FOR WRITING SEVERAL PROPERTIES TO UNIT 4 (PROF)
*---- IADS  IVEL IRKF IPHSE
      0     0    0    0 
CC
CC*******************************************************************
CC                                                                  *
CC    RESERVOIR PROPERTIES                                          *
CC                                                                  *
CC*******************************************************************
CC  
CC
CC MAX. SIMULATION TIME ( days)
*---- TMAX
      150.
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      0     0      0    0
CC
CC CONSTANT POROSITY
*----PORC1
     .25              
CC
CC CONSTANT X-PERMEABILITY (MILIDARCY) FOR LAYER K = 1,NZ
*----PERMX(1)
     1000.
CC
CC CONSTANT Y-PERMEABILITY (MILIDARCY) FOR LAYER K = 1,NZ
*----PERMY(1)
     1000.
CC
CC CONSTANT Z-PERMEABILITY
*----PERMZC (MILIDARCY)
     100.
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
     145.04
CC
CC CONSTANT INITIAL WATER SATURATION
*----SWI
     1.0
CC
CC CONSTANT CHLORIDE AND CALCIUM CONCENTRATIONS (MEQ/ML)
*----C50       C60
    0.05892   .00271
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
*--- 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.     .016    0.   .016     0.0   .016
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
     .65   .9   0.     0.
CC
CC THE CSE SLOPE PARAMETER FOR CALCIUM AND ALCOHOL 1 AND ALCOHOL 2
*----BETA6  BETA7  BETA8
     0.0    0.    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
     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.477
cc
cc flag to allow solubility of oil in water
*---- imass  icor
        0    0
CC
CC CAPILLARY DESATURATION PARAMETERS FOR PHASE 1, 2, AND 3
*----ITRAP   T11        T22        T33
     0       1865.      28665.46      364.2 
cc
cc flag for direction of rel. perm. and pc curves
*---- 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
     .147    .4    .147
CC
CC CONSTANT ENDPOINT REL. PERM. OF PHASES 1,2,AND 3 AT LOW CAPILLARY NO.
*----P1RW  P2RW     P3RW
     .13771   .91480   .13771
CC
CC CONSTANT REL. PERM. EXPONENT OF PHASES 1,2,AND 3 AT LOW CAPILLARY NO.
*----E1W     E2W  E3W
     2.1817   1.40475   2.1817
CC
CC WATER AND OIL VISCOSITY , RESERVOIR TEMPERATURE
*----VIS1   VIS2  TEMPV
     0.42   2.5   0.0
CC
CC VISCOSITY PARAMETERS
*----ALPHA1 ALPHA2  ALPHA3  ALPHA4  ALPHA5
     2.25811963     0.99671125      0.0    .9      .7
CC
CC PARAMETERS TO CALCULATE POLYMER VISCOSITY AT ZERO SHEAR RATE
*----AP1     AP2     AP3
     74.3     1106.   6376.
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
     0.      10.    1.8
CC
CC FLAG FOR POLYMER PARTITIONING, PERM. REDUCTION PARAMETERS
*----IPOLYM EPHI3 EPHI4 BRK    CRK
     1      1.    0.9    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 
     .4368  .3462 0.3462 .4333  .346  0.  1  
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 PARAMETERS, CPC
*----CPC 
     0.
CC
CC CAPILLARY PRESSURE PARAMETERS, EPC
*---- EPC
      2.
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) D(12)
     0.   0.   0.   0.   0.   0.   0.   0.   0.   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) D(4) D(5) D(6) D(7) D(8) D(9)  D(10)  D(11) D(12)
     0.   0.   0.   0.   0.   0.   0.   0.   0.    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) D(4) D(5) D(6) D(7) D(8) D(9) D(10)  D(11) D(12)
     0.   0.   0.   0.   0.   0.   0.   0.   0.   0.     0.     0.
     0.   0.   0.   0.   0.   0.   0.
CC
CC LONGITUDINAL AND TRANSVERSE DISPERSIVITY OF PHASE 1
*----ALPHAL(1)     ALPHAT(1)
     0.02          0.0
CC
CC LONGITUDINAL AND TRANSVERSE DISPERSIVITY OF PHASE 2