keurfonluu / toughio

Pre- and post-processing Python library for TOUGH
Other
56 stars 8 forks source link

How can i define elative permeability curve by Tab data? #151

Open BourbonTiger opened 1 year ago

BourbonTiger commented 1 year ago

Hi, I found that when defining rock properties relative permeability curve and capillary pressure curve there are only examples of input using formulas. Is there a function in toughio to define relative permeability curve and capillary pressure curve using Tab data? The data format is as follows: ROCKS----1--------2--------3--------4--------5--------6--------7----*----8 Rock0 2 2600.e00 0.148 1.019E-14 1.019E-14 1.019E-14 2.51 920. 4.4e-10 15 .000 .00 28.5 20, 20 0.2100 0.0000 1.0000 0.000 0.2365 0.0009 0.9580 0.000 0.2630 0.0037 0.9001 0.000 0.2895 0.0083 0.8344 0.000 0.3160 0.0148 0.7449 0.000 0.3425 0.0231 0.6616 0.000 0.3690 0.0333 0.5844 0.000 0.3955 0.0454 0.5135 0.000 0.4220 0.0593 0.4486 0.000 0.4485 0.0750 0.3900 0.000 0.4750 0.0926 0.3375 0.000 0.5015 0.1120 0.2960 0.000 0.5280 0.1333 0.2511 0.000 0.5545 0.1565 0.2172 0.000 0.5810 0.1815 0.1794 0.000 0.6075 0.2083 0.1478 0.000 0.6340 0.2370 0.1223 0.000 0.6605 0.2676 0.0931 0.000 0.6870 0.3000 0.0000 0.000 1.0000 0.8000 0.0000 0.000 0.000 0.000 1.000 0.000 0.030 0.00002 0.916 0.000 0.072 0.00053 0.803 0.000 0.114 0.00249 0.697 0.000 0.157 0.00693 0.599 0.000 0.199 0.01444 0.508 0.000 0.241 0.02579 0.425 0.000 0.283 0.04157 0.349 0.000 0.326 0.06284 0.280 0.000 0.368 0.08902 0.219 0.000 0.410 0.12084 0.166 0.000 0.452 0.15853 0.120 0.000 0.494 0.20219 0.081 0.000 0.537 0.25309 0.050 0.000 0.579 0.30878 0.026 0.000 0.621 0.37021 0.010 0.000 0.663 0.43709 0.002 0.000 0.706 0.51083 0.000 0.000 0.748 0.58752 0.000 0.000 0.790 0.66826 0.000 0.000

keurfonluu commented 1 year ago

Hi @BourbonTiger,

toughio does not support such tabular data for capillarity and relative permeability. I don't pretend to know all the variations of TOUGH, but I have never seen these curves being defined in such fashion in block ROCKS. Is this input file working in the version of TOUGH you are using? If yes, what version exactly are you using?

BourbonTiger commented 1 year ago

thank you, My coworker and I confirmed that toughreact is being used, but we're not sure if the version we're using has been modified yet.

BourbonTiger commented 1 year ago

Relative permeability curve data in tab data format is useful in reservoir numerical simulation.

keurfonluu commented 1 year ago

I am not familiar with TOUGHREACT, but I've looked through the manual, and despite some modifications to the block ROCKS (which are already supported by toughio), there is no mention of such enhancement for relative permeability format.

Relative permeability curve data in tab data format is useful in reservoir numerical simulation.

I am not saying otherwise. However, toughio is designed for unmodified versions of TOUGH codes so that it can be used by most people, and is only supporting these official versions. If you are using a customized version of TOUGH3, feel free to fork toughio and support these features in your fork ;)

Usually, we fit these curves with a model that exists in TOUGH, instead of providing raw data points.

BourbonTiger commented 1 year ago

In addition, when I use toughio-export to export vtk file and visualize it with paraview, paraview reports error: Unsupported data type: vtktypeint32ERROR: In vtkDataReader.cxx, line 1977vtkUnstructuredGridReader (000001CA19EBA3A0): Unsupported data type: vtktypeint32, do you know how to fix it?

keurfonluu commented 11 months ago

In addition, when I use toughio-export to export vtk file and visualize it with paraview, paraview reports error: Unsupported data type: vtktypeint32ERROR: In vtkDataReader.cxx, line 1977vtkUnstructuredGridReader (000001CA19EBA3A0): Unsupported data type: vtktypeint32, do you know how to fix it?

Not sure why that happens. Please export to VTU format instead.