Hi again,
it is really hard to evaluate if gempy is suited for this with this little information. I will try to give a general answer:

Gempy does prefer layer stacks and does not have any specific method to model intrusion bodies. Depending on many factors, mainly the shape of your intrusions, their relation to the surrounding strucuture and also the amount of these bodies you want to include in your model it might still be very possible to model with Gempy. The ebst way to evaluate this would be to understand the modeling principles underlying gempy and the structure that is used to define the interaction between structural groups and elements. I would recommend ourvideo tutorial series to get a grasp if you have not already done so.

If a pure gempy solution is not suitable there is also the option to use a gempy model as a basis and either add objects to i or perfom some kind of geostatistical operation within certain units as a post processing step to add intrusion bodies.

I hope this helps, feel free to ask follow-ups if anything is unclear.

Cheers,
Jan

Hi,
again this might be an unsatisfying answer - but I can't give a clear yes or no if this is generally possible. Some general points. Just working with a map (no depth data) is always tricky, You would also need orientation data for you units which Gempy rewuires. As with every complex strucutral model you might need to somplify some areas or make some manual adjustments. Unit thickness can not be set as input in Gempy.

Inserting intrusions in Gempy is not possible. You could either model them as a new unit on top of you othere units and then just cut it with the topography (this would be fully within GemPy but I cant guarantee you would be satisfied with the result), or just take the gempy solution and add your intrusions as objects (e.g by changing the resulting lithology block or inserting new meshes) afterwards (outside of gempy).

Best regards,
Jan

Let us know if you have any problems/questions along the way!