Can we handle finite/small faults (that do not go across through geomodel extents) in GemPy?

[geomodeller]

Hi,

I wish to ask if there is a way to include finite faults in GemPy. For example, the following figure shows a top view of geomodel with three faults.

I know how to include faults like "fault #1" that go across the extent. However, "fault #2 and #3" does not go up to the geomodel extent, but when importing in GemPy, I've seen that it will create faults like:

In a way to make them finite faults, I was thinking of adding two imaginary faults (i.e., green) that can define the ends of individual faults, as shown below:

To make those green faults not affect the surfaces, we set the faults to have "isFinite == False". But it didn't work as we wished :(

Any help would be highly appreciated!

Thank you,
Best!

Honggeun Jo

[geomodeller]

Oh, we set green faults' "isFinite" == True, not False

[Japhiolite]

Hi @whghdrms
Yes, indeed you can set faults as finite with the "isFinite" method. Be cautious, though, that method needs some work, don't trust the results too much.
Also, the faults will still be plotted across the whole model domain, but that is just visualization. Their influence on the model is more or less confined to the extent of the input-points for the respective fault.

[geomodeller]

Thank you for your prompt response! After posting this question, I have been 'playing' with GemPy to have two perpendicular faults where one of the faults goes across the other one while not affecting strata. Here is what I got:

On top of setting "isFinite == True," we also need to (1) sequence of fault #2 (i.e., an imaginary fault that define the boundary) comes after fault #1 and (2) change the relationship of fault as shown below:

I also found some visual artifacts, but when transformed to NumPy.array, the result looks reasonable. Thank you again for your help!

Best,
Honggeun Jo

[Japhiolite]

Good looking model Honggeun, nicely presents fault relationships. The visual artifacts may occur yes, due to the main_fault stopping at main_fault_2 and being deflected along the surface of main_fault_2 using a marching cubes algorithm.