ImplicitCoKriging类的set_is_fault()
函数用于将模型中的断层标志设置为True,以便在模拟面波时对断层进行处理。
set_is_fault(fault_idx: int)
fault_idx
: 断层编号,从0开始。默认值为None。该方法没有返回值。
import gempy as gp
# 创建模型
model = gp.create_model("my_model")
# 添加不同层的区域
gp.init_data(model, [0, 50, 100], [0, 50, 100], [0, 50, 100])
# 定义深度结构
gp.set_interpolator(model,
output=['geology'],
compile_theano=True,
theano_optimizer='fast_compile',
dtype='float32')
# 将区域标记为断层
model.set_is_fault(0)
# 得到默认值的地质模型
gp.compute_model(model)
# 进行地质横切面绘图
gp.plot_2d(model,
cell_number=[25],
show_data=True,
show_results=True,
direction='y')