这是 gempy.core.model.Project 类的一个方法,用于设置地层模型的底部关系。
set_bottom_relation(series: Union[int, str], relation: Union[str, np.ndarray])
series:int 或 str,代表地层序列的名称或序号。relation:str 或 numpy.ndarray,代表地层序列顶部和底部之间的关系。关系可选项包括 'Erosion'、'Onlap'、'Unconformity'、'Nonconformity'、'Conformable',或者一个与 gempy.core.data_modules.grid_modules.Grid 类相同形状和顺序的布尔型数组。无返回值。该方法直接对 gempy.core.data_modules.grid_modules.Grid 对象中对应的 input_data 数组进行修改。
import gempy as gp
# 创建地层模型
geo_model = gp.create_model('ExampleModel')
# 设置地层系列
gp.init_data(geo_model, [0, 100, 200], [0, 10, 20], [1, 0], [10, 10])
# 设置地层面的关系
geo_model.set_bottom_relation(series=0, relation='Onlap')
geo_model.set_bottom_relation(series=1, relation=np.array([[False, True], [False, False]]))
# 显示该地层模型
gp.plot_2d(geo_model)
在上述示例中,我们首先创建了一个 ExampleModel 的地层模型,并使用 init_data 方法设置了两个不同的地层系列。然后,我们使用 set_bottom_relation 方法设置了地层系列 0 的底部关系为 'Onlap',并将地层系列 1 的底部关系设置为一个 2x2 布尔型数组。最后,我们调用了 plot_2d 方法将结果可视化。
TypeError:系列名称或编号不正确;地层序列顶部和底部之间的关系类型不正确或与 Grid 对象中的形状不匹配。ValueError:所提供的顶部和底部关系不符合地质学上的逻辑或不允许推断出正确的层序。