search

OPM / ResInsight

3D viewer and post processing of reservoir models
http://resinsight.org/
GNU General Public License v3.0
174 stars 92 forks source link

Flow Diagnostics result in "Inf" for non-isothermal single-phase flow #4037

Open skahroba opened 5 years ago

skahroba commented 5 years ago

Using flow diagnostics to calculate TOF for the water phase in a non-isothermal single-phase flow results in inf at all time steps. capture

magnesj commented 5 years ago

Hi @skahroba, thanks for your feedback. @atgeirr , do you have any input on this issue?

magnesj commented 5 years ago

@skahroba : If possible, please upload a test simulation case reproducing your problems.

@bska : Do you have any input on this?

bska commented 5 years ago

@bska : Do you have any input on this?

The most obvious cause of infinite time-of-flight is dividing by zero—which could either be the result of summing zero rates (no inflow, or no flow at all). Did this simulation output phase fluxes or would interface fluxes be calculated after the fact from pressure/potential differentials and derived densities/mobilities?

If possible, please upload a test simulation case reproducing your problems.

Yes, that would definitely be useful.

skahroba commented 5 years ago

If possible, please upload a test simulation case reproducing your problems.

Thanks for the follow up. Attached is a simple simulation test case. nonisothermal_test_case.zip

bska commented 5 years ago

Thanks for the follow up. Attached is a simple simulation test case.

Much appreciated. I've had a look at the information and I think this may be a simulator (OPM Flow) problem rather than a ResInsight problem. In particular, all of the well data looks like this

 'XWEL    '         260 'DOUB'                                                  
   0.00000000000000D+00   0.38399999922420D+04   0.00000000000000D+00           
   0.38399999922420D+04                nanD-2147483648   0.00000000000000D+00   
   0.25689386621320D+02   0.10000000000000D+01   0.00000000000000D+00

and the fifth item, the PROD well's reservoir voidage rate, is somewhat suspicious (NaN). There's also the connection flow rate data which looks a lot like

'XCON    '         116 'DOUB'                                                  
   0.00000000000000D+00   0.38399999881190D+04   0.00000000000000D+00           
   0.00000000000000D+00   0.00000000000000D+00   0.00000000000000D+00

in which everything except the connection's water flow rate at surface conditions is zero.

That, in turn, means that there's no source terms for the wells and that might well explain why the flow diagnostics calculator is unable to derive reasonable time-of-flight values.

Are you able to run ECLIPSE on the same case, and if so, do you get other/better results then?

atgeirr commented 5 years ago

Just to be clear: OPM Flow is currently not writing fluxes to restart files, so ResInsight is not able to run flow diagnostics on Flow output (this will be fixed eventually).

skahroba commented 5 years ago

Are you able to run ECLIPSE on the same case, and if so, do you get other/better results then?

Many thanks for your comprehensive reply. I have tried the same case with E100 and E300 and in both cases ResInsight crashes. Note that in all three cases, heat front development as well as water rate in the producer (surface and reservoir conditions) match with each other.