求cmg大神帮忙改错！SAGD数值模拟

求大神帮忙，本人是做加拿大Christina Lake油藏的历史拟合，但是model一直不能拟合上，日产量特别小，求大神帮忙
以下是dat文件

INUNIT SI
WSRF WELL 1
WSRF GRID TIME
WSRF SECTOR TIME
OUTSRF GRID CCHLOSS CCHLOSSCUM FLUXRC FPOROS KRG KRO KRW MASDENO PRES SG SO
            SW TEMP VISO X Y
OUTSRF WELL LAYER ALL
OUTSRF WELL COMPONENT ALL
OUTSRF WELL DOWNHOLE
WPRN GRID 0
OUTPRN GRID NONE
OUTPRN RES NONE
**  Distance units: m
RESULTS XOFFSET           0.0000
RESULTS YOFFSET           0.0000
RESULTS ROTATION           0.0000  **  (DEGREES)
RESULTS AXES-DIRECTIONS 1.0 -1.0 1.0
** ***************************************************************************
** Definition of fundamental cartesian grid
** ***************************************************************************
GRID VARI 51 1 30
KDIR DOWN
DI IVAR
51*1.5
DJ JVAR
750
DK ALL
1530*1
DTOP
51*359
PERMI CON         5000
**  0 = null block, 1 = active block
NULL CON            1
POR CON         0.34
PERMK EQUALSI * 0.8
**  0 = pinched block, 1 = active block
PINCHOUTARRAY CON            1
PERMJ EQUALSI
END-GRID
ROCKTYPE 1
PRPOR 1300
CPOR 1e-6
ROCKCP 2.01e6 0
THCONR 2.74e5
THCONW 5.35e4
THCONO 1.15e4
THCONG 2.5e3
THCONMIX COMPLEX
HLOSSPROP OVERBUR 2.35e6 1.5e5
         UNDERBUR 2.35e6 1.5e5
** Model and number of components
** Model and number of components
MODEL 3 3 3 1
COMPNAME 'Water' 'Dead_Oil' 'Soln_Gas'
CMM
0 0.272455 0.0188305
PCRIT
0 0 4585.25
TCRIT
0 0 -75.0825
KV1
0 0 525180
KV2
0 0 -0.0940592
KV3
0 0 210.162
KV4
0 0 -879.991
KV5
0 0 -265.99
PRSR 101
TEMR 13
PSURF 101
TSURF 16.85
MASSDEN
1009.41 932.499 161.049
CP
4.88997e-007 1e-005 1e-005
CT1
-1.02097e-005 0.00118773 0.00118773
AVG
0 0 3.69196e-005
BVG
0 0 1
VISCTABLE
**      temp                                 
            5         0  14627862.15       115  ** Live oil visc (P=1300) = 262.21
           12         0  4420687.104      98.1  ** Live oil visc (P=1300) = 128.328
           20         0  1291721.594      68.4  ** Live oil visc (P=1300) = 58.7954
           30         0  331656.0328      54.1  ** Live oil visc (P=1300) = 6.05785
           40         0  100864.5923      43.4  ** Live oil visc (P=1300) = 1.6005
           50         0   35419.2117      35.2  ** Live oil visc (P=1300) = 0.701469
           60         0   14062.1725      28.9  ** Live oil visc (P=1300) = 0.410779
           70         0    6203.0386        24  ** Live oil visc (P=1300) = 0.286455
           80         0    2996.2711      20.1  ** Live oil visc (P=1300) = 0.22295
           90         0    1565.6454        17  ** Live oil visc (P=1300) = 0.186406
          100         0     875.9733      14.6  ** Live oil visc (P=1300) = 0.16354
          110         0     520.2501      12.6  ** Live oil visc (P=1300) = 0.14833
          120         0     325.5854      10.9  ** Live oil visc (P=1300) = 0.137732
          130         0     213.3676      9.56  ** Live oil visc (P=1300) = 0.130078
          140         0     145.6378      8.43  ** Live oil visc (P=1300) = 0.124389
          150         0     103.0622      7.49  ** Live oil visc (P=1300) = 0.120059
          160         0      75.3141       6.7  ** Live oil visc (P=1300) = 0.116699
          170         0       56.737      6.03  ** Live oil visc (P=1300) = 0.114048
          180         0      43.6994      5.46  ** Live oil visc (P=1300) = 0.111927
          190         0      34.5021      4.97  ** Live oil visc (P=1300) = 0.110208
          200         0      27.8104      4.55  ** Live oil visc (P=1300) = 0.108801
          210         0      22.8386      4.18  ** Live oil visc (P=1300) = 0.107637
          220         0      19.0741      3.86  ** Live oil visc (P=1300) = 0.106666
          230         0      16.1745      3.58
          240         0      13.9058      3.34
          250         0      12.1056      3.12
          260         0      10.6583      2.93
          270         0       9.4811      2.76
          280         0       8.5129      2.61
          300         0       7.0349      2.36
VSMIXCOMP 'Soln_Gas'
VSMIXENDP 0.00662251 0.058
VSMIXFUNC 0.00662251 0.0119493 0.0172194 0.0224344 0.0275941 0.0327005 0.0377557 0.0427622 0.0477216 0.0526689 0.0576162
ROCKFLUID
RPT 1 WATWET
**        Sw          krw        krow
SWT
         0.15            0       0.834
     0.190625  0.000390625    0.733008
     0.271875   0.00351562     0.55057
       0.3125      0.00625    0.469125
     0.353125   0.00976563    0.394195
      0.39375    0.0140625    0.325781
     0.434375    0.0191406    0.263883
        0.475        0.025      0.2085
     0.515625    0.0316406    0.159633
      0.55625    0.0390625    0.117281
     0.596875    0.0472656   0.0814453
       0.6375      0.05625    0.052125
     0.678125    0.0660156   0.0293203
      0.71875    0.0765625   0.0130313
     0.759375    0.0878906  0.00325781
          0.8          0.1           0
**        Sl       krg        krog
SLT
        0.155     0.992           0
     0.207813  0.871875  0.00325781
     0.260625    0.7595   0.0130312
     0.313438  0.654875   0.0293203
      0.36625     0.558    0.052125
     0.419063  0.468875   0.0814453
     0.471875    0.3875    0.117281
     0.524687  0.313875    0.159633
       0.5775     0.248      0.2085
     0.630313  0.189875    0.263883
     0.683125    0.1395    0.325781
     0.735938  0.096875    0.394195
      0.78875     0.062    0.469125
     0.841562  0.034875     0.55057
     0.894375    0.0155    0.638531
     0.947187  0.003875    0.733008
            1         0       0.834
KRTEMTAB SORG SORW
**      TEMP      SORG      SORW
           13      0.22       0.5
          200      0.06       0.2
INITIAL
VERTICAL DEPTH_AVE

INITREGION 1
REFPRES 1500
REFDEPTH 359
DWOC 389
DGOC 359
TEMP CON           13
SO CON         0.85
MFRAC_OIL 'Soln_Gas' CON    0.0476449
MFRAC_OIL 'Dead_Oil' CON     0.952355
NUMERICAL
UPSTREAM KLEVEL
RUN
DATE 2007 2 1

DTWELL 0.001
**
WELL  'Injector'
INJECTOR MOBWEIGHT IMPLICIT 'Injector'
INCOMP  WATER  1.0  0.0  0.0
TINJW  217.0
QUAL  0.8
OPERATE  MAX  BHP  2100.0  CONT REPEAT
**          rad  geofac  wfrac  skin
GEOMETRY  J  0.086  0.249  1.0  0.0
      PERF      GEOA  'Injector'
** UBA              ff          Status  Connection  
    1 1 24         1.0  OPEN    FLOW-FROM  'SURFACE'
**
WELL  'Prod'
PRODUCER 'Prod'
OPERATE  MIN  BHP  1400.0  CONT REPEAT
OPERATE  MAX  STEAM  1.0  CONT REPEAT
**          rad  geofac  wfrac  skin
GEOMETRY  J  0.086  0.249  1.0  0.0
      PERF      GEOA  'Prod'
** UBA              ff          Status  Connection  
    1 1 29         1.0  OPEN    FLOW-TO  'SURFACE'
HTWELL 'Injector'
   HTWTEMP 240
HTWELL 'Prod'
   HTWTEMP 240
DATE 2007 3  1.00000
DATE 2007 4  1.00000
DATE 2007 5  1.00000
HTWELL 'Injector' OFF
HTWELL 'Prod' OFF
DATE 2007 6  1.00000
DATE 2007 7  1.00000
DATE 2007 8  1.00000
DATE 2007 9  1.00000
DATE 2007 10  1.00000
DATE 2007 11  1.00000
DATE 2007 12  1.00000
DATE 2008 1  1.00000
DATE 2008 2  1.00000
DATE 2008 3  1.00000
DATE 2008 4  1.00000
DATE 2008 5  1.00000
DATE 2008 6  1.00000
DATE 2008 7  1.00000
DATE 2008 8  1.00000
DATE 2008 9  1.00000
DATE 2008 10  1.00000
DATE 2008 11  1.00000
DATE 2008 12  1.00000
DATE 2009 1  1.00000
DATE 2009 2  1.00000
DATE 2009 3  1.00000
DATE 2009 4  1.00000
DATE 2009 5  1.00000
DATE 2009 6  1.00000
DATE 2009 7  1.00000
DATE 2009 8  1.00000
DATE 2009 9  1.00000
DATE 2009 10  1.00000
DATE 2009 11  1.00000
DATE 2009 12  1.00000
DATE 2010 1  1.00000
DATE 2010 2  1.00000
DATE 2010 3  1.00000
DATE 2010 4  1.00000
DATE 2010 5  1.00000
DATE 2010 6  1.00000
DATE 2010 7  1.00000
DATE 2010 8  1.00000
DATE 2010 9  1.00000
DATE 2010 10  1.00000
DATE 2010 11  1.00000
DATE 2010 12  1.00000
DATE 2011 1  1.00000
DATE 2011 2  1.00000
DATE 2011 3  1.00000
DATE 2011 4  1.00000
DATE 2011 5  1.00000
DATE 2011 6  1.00000
DATE 2011 7  1.00000
DATE 2011 8  1.00000
DATE 2011 9  1.00000
DATE 2011 10  1.00000
DATE 2011 11  1.00000
DATE 2011 12  1.00000
DATE 2012 1  1.00000
DATE 2012 2  1.00000
DATE 2012 3  1.00000
DATE 2012 4  1.00000
DATE 2012 5  1.00000
DATE 2012 6  1.00000
DATE 2012 7  1.00000
DATE 2012 8  1.00000
DATE 2012 9  1.00000
DATE 2012 10  1.00000
DATE 2012 11  1.00000
DATE 2012 12  1.00000
DATE 2013 1  1.00000
DATE 2013 2  1.00000
DATE 2013 3  1.00000
DATE 2013 4  1.00000
DATE 2013 5  1.00000
DATE 2013 6  1.00000
DATE 2013 7  1.00000
DATE 2013 8  1.00000
DATE 2013 9  1.00000
DATE 2013 10  1.00000
DATE 2013 11  1.00000
DATE 2013 12  1.00000
DATE 2014 1  1.00000
DATE 2014 2  1.00000
DATE 2014 3  1.00000
DATE 2014 4  1.00000
DATE 2014 5  1.00000
DATE 2014 6  1.00000
DATE 2014 7  1.00000
DATE 2014 8  1.00000
DATE 2014 9  1.00000
DATE 2014 10  1.00000
DATE 2014 11  1.00000
DATE 2014 12  1.00000
DATE 2015 1  1.00000
DATE 2015 2  1.00000
DATE 2015 3  1.00000
DATE 2015 4  1.00000
DATE 2015 5  1.00000
DATE 2015 6  1.00000
DATE 2015 7  1.00000
DATE 2015 8  1.00000
DATE 2015 9  1.00000
DATE 2015 10  1.00000
DATE 2015 11  1.00000
RESULTS PVTIMEX VISCREGION 1
RESULTS PVTIMEX PVTREGION 1 FALSE
RESULTS PVTIMEX TABLECOLS P RS BO BG VISO VISG DENOIL DENGAS
RESULTS PVTIMEX TABLE 101.325 0.541318 0.996153 0.984998 172.931 0.0104224 930.952 0.806529
RESULTS PVTIMEX TABLE 181.237 0.734647 0.996461 0.549432 169.893 0.010433 930.819 1.44591
RESULTS PVTIMEX TABLE 261.149 0.937119 0.996784 0.380432 166.798 0.0104447 930.679 2.08823
RESULTS PVTIMEX TABLE 341.06 1.14736 0.997121 0.290627 163.676 0.0104573 930.531 2.73351
RESULTS PVTIMEX TABLE 420.972 1.36439 0.997469 0.234915 160.548 0.0104706 930.38 3.38177
RESULTS PVTIMEX TABLE 500.884 1.58748 0.997828 0.196981 157.428 0.0104846 930.223 4.03304
RESULTS PVTIMEX TABLE 580.795 1.81604 0.998197 0.169484 154.328 0.0104991 930.061 4.68734
RESULTS PVTIMEX TABLE 660.707 2.04961 0.998575 0.148639 151.258 0.0105142 929.894 5.34471
RESULTS PVTIMEX TABLE 740.618 2.2878 0.998962 0.132291 148.223 0.0105298 929.724 6.00516
RESULTS PVTIMEX TABLE 820.53 2.53028 0.999357 0.119128 145.23 0.0105458 929.549 6.66872
RESULTS PVTIMEX TABLE 900.442 2.77679 0.999759 0.108301 142.283 0.0105623 929.371 7.33541
RESULTS PVTIMEX TABLE 980.353 3.02708 1.00017 0.0992388 139.386 0.0105792 929.188 8.00523
RESULTS PVTIMEX TABLE 1060.27 3.28094 1.00059 0.0915424 136.539 0.0105965 928.999 8.67827
RESULTS PVTIMEX TABLE 1140.18 3.53819 1.00101 0.0849242 133.747 0.0106143 928.814 9.35457
RESULTS PVTIMEX TABLE 1220.09 3.79866 1.00144 0.0791734 131.009 0.0106324 928.622 10.034
RESULTS PVTIMEX TABLE 1300 4.0622 1.00188 0.0741295 128.328 0.0106509 928.423 10.7168
RESULTS PVTIMEX TABLE 2040 6.62872 0.994521 0.0461902 128.328 0.0108395 935.293 17.1991
RESULTS PVTIMEX TABLE 2780 9.38023 0.987268 0.0331219 128.328 0.0110581 942.163 23.985
RESULTS PVTIMEX TABLE 3520 12.2779 0.980121 0.0255475 128.328 0.0113071 949.034 31.0962
RESULTS PVTIMEX TABLE 4260 15.2977 0.973077 0.020606 128.328 0.011588 955.904 38.5533
RESULTS PVTIMEX TABLE 5000 18.4229 0.966133 0.0171307 128.328 0.0119033 962.774 46.3746
RESULTS PVTIMEX TRES 13
RESULTS PVTIMEX BPP 15
RESULTS PVTIMEX BWI 0.999453
RESULTS PVTIMEX DENSITYWATER 1010
RESULTS PVTIMEX VISCOSITYWATER 1.29792
RESULTS PVTIMEX WATERCVW 0
RESULTS PVTIMEX DENSITYOIL 926.941
RESULTS PVTIMEX GASGRAVITY 0.65
RESULTS PVTIMEX WATERCOMP 4.88997e-007
RESULTS PVTIMEX REFPW 1300
RESULTS PVTIMEX CO 1e-005
RESULTS PVTIMEX CVO 0
RESULTS PVTIMEX RATIODEADPVT 1
RESULTS PVTIMEX VISCPRESSURE 101.3
RESULTS PVTIMEX COMPOSITION 2 0.952355 0.0476449
RESULTS PVTIMEX KVALUETEMP FALSE 400 -99999 0 0.264
RESULTS PVTIMEX END

RESULTS SPEC 'Permeability I'  
RESULTS SPEC SPECNOTCALCVAL -99999      
RESULTS SPEC REGION 'All Layers (Whole Grid)'
RESULTS SPEC REGIONTYPE 'REGION_WHOLEGRID'
RESULTS SPEC LAYERNUMB 0
RESULTS SPEC PORTYPE 1
RESULTS SPEC CON 5000        
RESULTS SPEC SPECKEEPMOD 'YES'
RESULTS SPEC STOP


RESULTS SPEC 'Permeability J'  
RESULTS SPEC SPECNOTCALCVAL -99999      
RESULTS SPEC REGION 'All Layers (Whole Grid)'
RESULTS SPEC REGIONTYPE 'REGION_WHOLEGRID'
RESULTS SPEC LAYERNUMB 0
RESULTS SPEC PORTYPE 1
RESULTS SPEC EQUALSI 0 1           
RESULTS SPEC SPECKEEPMOD 'YES'
RESULTS SPEC STOP


RESULTS SPEC 'Permeability K'  
RESULTS SPEC SPECNOTCALCVAL -99999      
RESULTS SPEC REGION 'All Layers (Whole Grid)'
RESULTS SPEC REGIONTYPE 'REGION_WHOLEGRID'
RESULTS SPEC LAYERNUMB 0
RESULTS SPEC PORTYPE 1
RESULTS SPEC EQUALSI 1 0.8         
RESULTS SPEC SPECKEEPMOD 'YES'
RESULTS SPEC STOP


RESULTS SPEC 'Porosity'  
RESULTS SPEC SPECNOTCALCVAL -99999      
RESULTS SPEC REGION 'All Layers (Whole Grid)'
RESULTS SPEC REGIONTYPE 'REGION_WHOLEGRID'
RESULTS SPEC LAYERNUMB 0
RESULTS SPEC PORTYPE 1
RESULTS SPEC CON 0.34        
RESULTS SPEC SPECKEEPMOD 'YES'
RESULTS SPEC STOP


RESULTS SPEC 'Oil Saturation'  
RESULTS SPEC SPECNOTCALCVAL -99999      
RESULTS SPEC REGION 'All Layers (Whole Grid)'
RESULTS SPEC REGIONTYPE 'REGION_WHOLEGRID'
RESULTS SPEC LAYERNUMB 0
RESULTS SPEC PORTYPE 1
RESULTS SPEC CON 0.85        
RESULTS SPEC SPECKEEPMOD 'YES'
RESULTS SPEC STOP


RESULTS SPEC 'Temperature'  
RESULTS SPEC SPECNOTCALCVAL -99999      
RESULTS SPEC REGION 'All Layers (Whole Grid)'
RESULTS SPEC REGIONTYPE 'REGION_WHOLEGRID'
RESULTS SPEC LAYERNUMB 0
RESULTS SPEC PORTYPE 1
RESULTS SPEC CON 13         
RESULTS SPEC SPECKEEPMOD 'YES'
RESULTS SPEC STOP


RESULTS SPEC 'Grid Thickness'  
RESULTS SPEC SPECNOTCALCVAL -99999      
RESULTS SPEC REGION 'All Layers (Whole Grid)'
RESULTS SPEC REGIONTYPE 'REGION_WHOLEGRID'
RESULTS SPEC LAYERNUMB 0
RESULTS SPEC PORTYPE 1
RESULTS SPEC CON 1           
RESULTS SPEC SPECKEEPMOD 'YES'
RESULTS SPEC STOP


RESULTS SPEC 'Grid Top'  
RESULTS SPEC SPECNOTCALCVAL -99999      
RESULTS SPEC REGION 'Layer 1 - Whole layer'
RESULTS SPEC REGIONTYPE 'REGION_LAYER'
RESULTS SPEC LAYERNUMB 1
RESULTS SPEC PORTYPE 1
RESULTS SPEC CON 359         
RESULTS SPEC SPECKEEPMOD 'YES'
RESULTS SPEC STOP

978264313

我大概看了一下，也不知道对不对。模拟完成后，只有很小的一部分区域温度达到100度以上，这样的话，油的粘度还是很小，日产量当然就低了。又看了一下注入量，注热水，受井底压力的限制，除了刚开始注入量大一点，后面的注入量大概每天43m3。可能是注热量达不到。
当然了，数模还是需要从整体考虑。希望对你有用。

数模不是物模

不知道·，绑定