Created an incorrect geological formation

[wdk11]

I have many geological profile maps, and I extracted the geological segmentation line data from the geological profile maps. I used gempy for 3D modeling, and now I am verifying the correctness of gempy modeling. For example, I used the 3D model modeled by gempy to restore my geological profile lines. After comparing the modeling, I found that the Qgl geological layer was interpolated and appeared in the layers that should not have appeared, as shown in the red box below. The Qgl layer should not have appeared in the red box.


So how can I solve this problem and make the Qgl formation appear in the right place? Thank you.
This is my code and data:

import numpy as np
import gempy as gp
import gempy_viewer as gpv

data_path = 'C:/Users/Administrator/Desktop'

geo_model: gp.data.GeoModel = gp.create_geomodel(
project_name='Tutorial_ch1_1_Basics',
extent=[562298, 563540, 3032848, 3033348, 1800, 2700],
resolution=[50, 50, 50], # * Here we define the resolution of the voxels
# refinement=4, # * Here we define the number of octree levels. If octree levels are defined, the resolution is ignored.
importer_helper=gp.data.ImporterHelper(
path_to_orientations=data_path + "/orientations1.csv",
path_to_surface_points=data_path + "/points.csv",
)
)
gp.map_stack_to_surfaces(
gempy_model=geo_model,
mapping_object= # TODO: This mapping I do not like it too much. We should be able to do it passing the data objects directly
{
"Strat_Series": ('dmx', 'Qdl', 'Qgl')
# "Strat_Series": ('dmx', 'Qdl')
}
)
gp.set_section_grid(
grid=geo_model.grid,
section_dict={
'section1':([563915.18,3033408.73],[562473.68,3031254.95],[500,500]),
'section2':([563668.52,3033498.51],[562094.92,3032038.43],[500,500])
}
)
sol = gp.compute_model(geo_model,
engine_config=gp.data.GemPyEngineConfig(
backend=gp.data.AvailableBackends.PYTORCH,
use_gpu=True,
dtype="float64"
))

gpv.plot_2d(geo_model,section_names=['section1'])

orientations1.csv
points.csv

[javoha]

Hi,
how familar are you with the gempy infrastrucuture. In this case the problem appears to be that you grouped all your elements into a single strucutural group. This enforces a somewhat parallel behavior so that the element can not just end.
I suggest aking a look at the new gempy tutprial video series that explains most of the basic concepts in detail: First video tutorial

EDIT: I also converted this to a discussion as I am pretty sure it is not a bug or a similar issue.

[wdk11]

What actions should I take to achieve the effect in my original picture?
Do we need to set up a fault to forcibly terminate the geological line, or can other operations be done, such as setting up two series of structural groups? Thank you

[javoha]

As you want the red layer (in your gemy result) to fade out you probably need an uncomformity and yes this requires separating theelements into different strucutural groups. But again it would really help to watch the videos to understand these basic concepts of how gempy interpolates. The strucutural frame is very important.

[wdk11]

Thanks for your help.
I have successfully built the model I want