set_surface_order_from_solution
是Gempy的ImplicitCoKriging
模块中的一个函数,用于从模型的解决方案中设置地层表面的顺序。
set_surface_order_from_solution(self, solution)
solution
: numpy数组,模型的解决方案。该函数使用模型的解决方案来为地层表面设置顺序。 下面给出一个使用set_surface_order_from_solution
的示例:
import gempy as gp
# 导入数据
geo_data = gp.create_data([0,2000,0,2000,-2000,0],[50,50,50],
path_o = '/home/user/gempy_tutorial/notebooks/data/input_data/simple_fault_model_orientations.csv',
path_i = '/home/user/gempy_tutorial/notebooks/data/input_data/simple_fault_model_points.csv')
# 实例化一个Gempy模型,并进行插值
model = gp.create_model('simple_fault_model')
model.set lithology(['0', '1'], colors=['red', 'blue'], plot=True)
model.set_regular_grid(0,2000,-2000,2000,0,-2000,cell_sizes=[20,20,20])
# Gempy模型编译
model.surfaces.topography.values = geo_data.surface_points.topography.values
model.add_surface('topography')
model.set_is_fault(['fault'], [0])
model.set_section_line((0, 1000), (2000, 1000))
model.set_interpolator(interp_data='default', compile_theano=True, theano_optimizer='fast_compile')
# 执行地层建模
sol = model.compute_model(verbose=False)
# 设置地层表面顺序
model.set_surface_order_from_solution(sol)
该函数没有明确的输出,但是它将更新模型实例中的地层表面顺序。