batch.data

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

CC*******************************************************************
CC                                                                  *
CC     BRIEF DESCRIPTION OF DATA SET : UTCHEM (VERSION 9.0  )       *
CC                                                                  *
CC*******************************************************************
CC                                                                  *
CC    batch calculation to obtain UTCHEM phase behavior             *
CC                                                                  *
CC  LENGTH (FT) :                PROCESS : SURFACTANT               *
CC  THICKNESS (FT) :              INJ. RATE (FT3/DAY) : 0.05        *
CC  WIDTH (FT) :                 COORDINATES : CARTESIAN            *
CC  POROSITY :                                                      *
CC  GRID BLOCKS :  5x 1x1                                           *
CC  DATE : 5/14/2000                                                *
CC                                                                  *
CC*******************************************************************
CC
CC*******************************************************************
CC                                                                  *
CC    RESERVOIR DESCRIPTION                                         *
CC                                                                  *
CC*******************************************************************
CC  
CC
*----RUNNO
 isofol/PCE
CC  
CC
*----HEADER
get the hand's parameter from batch samples (with pCE)
salinity of 700 ppm ca++ 
NO Correct fluid properties, only match of phase behavior parameters
CC
CC SIMULATION FLAGS
*---- IMODE IMES IDISPC ICWM ICAP IREACT IBIO ICOORD ITREAC ITC IGAS IENG 
        1    2    0      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
     5   1    1   0        0          
CC
CC  CONSTANT GRID BLOCK SIZE IN X, Y, AND Z
*----DX1       DY1      DZ1
      1       1        1          
CC
CC TOTAL NO. OF COMPONENTS, NO. OF TRACERS, NO. OF GEL COMPONENTS
*----N   no  NTW NTA  ngc NG  noth
     7   0  0   0   0    0    0
CC
CC  NAME OF SPECIES
*---- SPNAME(I) FOR I=1,N 
WATER
OIL
isofol SURF
POLYMER(no)
cl-
ca++
ipa
CC
CC FLAG INDICATING IF THE COMPONENT IS INCLUDED IN CALCULATIONS OR NOT
*----ICF(KC) FOR KC=1,N
     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  0  1  1  1 
CC
CC FLAG FOR PRES,SAT.,TOTAL CONC.,TRACER CONC.,CAP.,GEL, ALKALINE PROFILES
*----IPPRES IPSAT IPCTOT IPBIO IPCAP IPGEL IPALK IPTEMP IPOBS
      0      0      0      0      0     0    0     0       0 
CC
CC FLAG FOR WRITING SEVERAL PROPERTIES 
*----ICKL IVIS IPER ICNM  ICSE  IFOAMP  IHYSTP  INONEQ
      0     0    0    0    0    0        0       0  
CC
CC FLAG FOR WRITING SEVERAL PROPERTIES TO 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
       120
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
      1.            
CC
CC CONSTANT X-PERMEABILITY (MILIDARCY) FOR LAYER K = 1,NZ
*----PERMX(1)   
     1000000.  
CC
CC CONSTANT Y-PERMEABILITY (MILIDARCY) FOR LAYER K = 1,NZ
*----PERMY(1) 
     1000000      
CC
CC CONSTANT Z-PERMEABILITY
*----PERMZC (MILIDARCY)
    1000000. 
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
     1. 
CC
CC CONSTANT CHLORIDE AND CALCIUM CONCENTRATIONS (MEQ/ML)
*----C50 (700 mg/l Ca++)      C60
    0.035                  .0350
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.     .12     .0     .054   0.      0.171
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 lab. data including both calcium and alcohol
CC LOWER AND UPPER EFFECTIVE SALINITY FOR ALCOHOL 1 AND ALCOHOL 2
*----CSEL7  CSEU7(opt = 1000 ppm ca++)  CSEL8  CSEU8
     .04         0.06     0.   0.
CC
CC THE CSE SLOPE PARAMETER FOR CALCIUM AND ALCOHOL 1 AND ALCOHOL 2
*----BETA6(1)  BETA7 (-2.08) BETA8
     0.        0.    0.
CC
CC FLAG FOR ALCOHOL PART. MODEL AND PARTITION COEFFICIENTS
*----IALC  OPSK7O  OPSK7S  OPSK8O  OPSK8S
     0     0.04     2.8     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
     1 
CC
CC INTERFACIAL TENSION PARAMETERS
*--- chuh  ahuh
     0.3  13 
CC
CC LOG10 OF OIL/WATER INTERFACIAL TENSION 
*----XIFTW
     1.65
CC
CC MASS TRANSFER FLAG
*----IMASS  ICOR
     0       0
CC
CC CAPILLARY DESATURATION PARAMETERS FOR PHASE 1, 2, AND 3
*----ITRAP   T11        T22        T33
     0       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
    3*0 
CC
CC CONSTANT ENDPOINT REL. PERM. OF PHASES 1,2,AND 3 AT LOW CAPILLARY NO.
*----P1RW  P2RW  P3RW
    3*1. 
CC
CC CONSTANT REL. PERM. EXPONENT OF PHASES 1,2,AND 3 AT LOW CAPILLARY NO.
*----E1W  E2W  E3W
    3*1. 
CC
CC WATER AND OIL VISCOSITY , RESERVOIR TEMPERATURE
*----VIS1   VIS2  TEMPV
      1.    1.  0.
CC
CC VISCOSITY PARAMETERS
*----ALPHA1 ALPHA2  ALPHA3  ALPHA4  ALPHA5
     0.     0     0   0  0 
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  .433  0.433  0.433   0.433  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 
     0. 
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)
     0.         0 
CC
CC LONGITUDINAL AND TRANSVERSE DISPERSIVITY OF PHASE 2
*----ALPHAL(2)     ALPHAT(2)
     0.          0 
CC
CC LONGITUDINAL AND TRANSVERSE DISPERSIVITY OF PHASE 3
*----ALPHAL(3)     ALPHAT(3)
     0.          0 
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
     0.    .0    1000.  0.7    0.    100.  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 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   ITIME  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        1    0
CC
CC WELL NAME
*---- WELNAM
INJECTOR
CC
CC ICHEK, MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHEK PWFMIN   PWFMAX  QTMIN   QTMAX
      0    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    5   1     2       .5       0.     3     1         1       0
CC
CC WELL NAME
*---- WELNAM
PRODUCER
CC
CC ICHEK , MAX. AND MIN. ALLOWABLE BOTTOMHOLE PRESSURE AND RATE
*----ICHEK PWFMIN   PWFMAX  QTMIN   QTMAX
      0    0.0     5000.   0.0     50000.
CC adjust the opt. salinity based on the water conc
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    0.05    0.78  0. 0.02  0. 0.0273 0.0273   0.2 0.  0. 0.  0.        
     1    0.05    0.   1.  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   CUMHI1   WRHPV   WRPRF      RSTC
     500       5        5       0.5     5         50 
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