gempy-project / gempy

GemPy is an open-source, Python-based 3-D structural geological modeling software, which allows the implicit (i.e. automatic) creation of complex geological models from interface and orientation data. It also offers support for stochastic modeling to address parameter and model uncertainties.
https://gempy.org
European Union Public License 1.2
941 stars 232 forks source link

Problems exporting flac3D .f3grid file format #837

Closed nearfuture625 closed 2 months ago

nearfuture625 commented 9 months ago

Hi all,

Thank you for developing such a useful software as Gempy, which allowed me to model a four-layer geological formation with two faults. I noticed that f3grid mesh files that can be used in FLAC3D can be exported via export.export_flac3D_input (geo_model, path='', filename='') in gempy.However, after exporting to f3grid file, I found that the mesh file does not contain a fault structure, please tell me how to export the mesh file related to the fault

54ab4b681d4e2a61efcccf0354bb1c5 Looking forward to your answer

MoiseRousseau commented 9 months ago

@nearfuture625: Let's talk here about this problem and the other you send me and related to exporting to FLAC3D 6.00 (https://github.com/MoiseRousseau/Variational-Tet-Mesh/issues/1).

  1. For instance, only the block mesh are exported, and there is not current code (from my side at least) to export the fault zone. Doing so would requires to adapt the mesh to the fault and export the face groups, or to create kind of a new material and export it as a fault zone.
  2. It seems you successfully export a mesh and reimport in FLAC3D on this issue, but you fail in the other issue (https://github.com/MoiseRousseau/Variational-Tet-Mesh/issues/1). What is different between both issues ? Can you send your GemPy input file ?
nearfuture625 commented 9 months ago

@MoiseRousseau :Thanks for your reply, last time I proposed that the f3grid file exported using the export_flac3D_input() function provided in gempy could not be displayed in flac3D, I have solved it, the zone part is added before B8, the zone sequence number is added after B8, and the node column has 9 columns after the sequence number, I adjusted the order of the node columns, new_column_order = [1, 2, 4, 5, 3, 8, 6, 7] to output an F3Grid file that can be displayed in FLAC. I'm going to attach my modified export_flac3D_input() function file, but the problem I'm facing now is that there is a lack of a way to export a mesh file containing faults to analyze the stresses near the faults using FLAC3D. 1017p.csv 1017o.csv input_flac.txt Looking forward to your insights on the method of tomographic derivation

ottaviani2 commented 4 months ago

Thank you both for the amazing opportunity to export file .f3grid from gempy. I need to update the code, there is an asterisk that prevents the creation of zones.

The line out.write('\n*ZONES') should be out.write('\nZONES')

I attach the updated file input_flac.txt

javoha commented 2 months ago

Resolved? Please feel free to reopen if not. I am also not sure how much of this is a gempy issue. But definitely a useful discussion.