Feature request: VTK import #15

Activity

1. acroucher commented on Mar 26, 2018

   

   acroucher

   on Mar 26, 2018

   Owner

   More actions

   You mean import from VTK to t2grid? That should be possible. Off the top of my head, I think it would involve:

   • constructing block names in some sensible way
   • calculating the block volumes and centroids, and connection areas and distances
   • assigning rock types? not sure if you have any rock data in your VTK- could just leave this out and assign a default rock type everywhere
   • assign default PARAM data

   Implementation-wise, one question would be whether to use the VTK Python bindings or an alternative like meshio, to read the VTK mesh into Python. VTK may have geometric functions included for computing block volumes etc.? On the other hand, meshio is often more portable and less problematic than VTK.
2. tirons commented on Apr 5, 2018

   

   tirons

   on Apr 5, 2018 via email

   Author

   More actions

   Thanks! Sorry to be getting back so late. I've been using PyVTK bindings, and currently can already convert SUTRA and PETREL grids into VTK, and convert from VTK into SUTRA and PETREL/ECLIPSE. Basically, like you said, I just need to work on the conversion to TOUGH grids, which are specified in quite a different fashion than the columnar data in PETREL. - The block volumes are easy to calculate, there is a built in function for that in VTK. - Geometrical centre is also easy to calculate, connection distances are also easy, I need more familiarity with the TOUGH specifications to understand "connection areas" better. I'll look into meshio, I have mainly been using vts (vtkStructuredGrid) which I don't think is supported by meshio. I'll keep you appraised of my progress, if I have something worth sharing, I'll submit a pull request for pyTOUGH. Thanks!

   …

   -Trevor

   On 26 March 2018 at 16:53, Adrian Croucher ***@***.***> wrote: You mean import from VTK to t2grid? That should be possible. Off the top of my head, I think it would involve: - constructing block names in some sensible way - calculating the block volumes and centroids, and connection areas and distances - assigning rock types? not sure if you have any rock data in your VTK- could just leave this out and assign a default rock type everywhere - assign default PARAM data Implementation-wise, one question would be whether to use the VTK Python bindings or an alternative like meshio, to read the VTK mesh into Python. VTK may have geometric functions included for computing block volumes etc.? On the other hand, meshio is often more portable and less problematic than VTK. — You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub <#15 (comment)>, or mute the thread <https://github.com/notifications/unsubscribe-auth/APbWFJ_C0Xj3M8NwBlSragXJd-LVdH1Xks5tiXFqgaJpZM4StBko> .