map_stack_to_surfaces 方法用于根据地层序列将每个网格的层位映射到相应的地层表面上。
| 参数名 | 数据类型 | 描述 |
|---|---|---|
| stack_series | numpy.ndarray | 地层序列 |
| order_surface | numpy.ndarray | 排序后的地层表面 |
| id_formations | numpy.ndarray | 不同地层名的编号 |
| 参数名 | 数据类型 | 描述 |
|---|---|---|
| points | numpy.ndarray | 空间点坐标 |
| formation | numpy.ndarray | 地层名 |
| series | numpy.ndarray | 地层编号 |
| id_points | numpy.ndarray | 网格点编号 |
无返回值。
import gempy as gp
geo_model = gp.create_model('my_model')
# 指定地质模型的图层结构
gp.init_data(geo_model, [0, 20, 40, 60, 80], [0, 1, 1, 2, 3], [3, 3, 3, 4, 4],
[2, 2, 2, 1, 1], [2, 2, 2, 1, 1], [1, 1, 1, 0, 0],
njit=True)
# 对地质模型进行堆栈建模
gp.compute_model(geo_model)
# 将每个网格的层位映射到相应的地层表面上
gp.map_stack_to_surfaces(geo_model._grid.topography.values, geo_model.series.array_series,
geo_model.surfaces.df.index.array)
stack_series、order_surface 和 id_formations 输入参数应当与地质模型 geo_model 中的相应参数保持一致。J. R. Mathieu, M. Robbi, and L. M. Heagy, An overview of the Gempy software for automatic 3D geological modeling and inversion, in SEG Technical Program Expanded Abstracts 2019, 2019, pp. 4847–4851. doi:10.1190/segam2019-3210219.1.