Gempy 是一个开源 Python 库，用于生成完整的 3D 结构地质模型。该库是从界面、断层和层方向创建地质模型的完整开发，它还关联地质层序列以表示岩石侵入和断层顺序。
地质建模算法基于通用协同克里金插值法，并支持 Numpy、PyMC3 和 Theano 等高端 Python 数学库。
Gempy 创建了一个网格模型，可以使用 Matplotlib 将其可视化为 2D 截面，也可以将 3D 几何对象可视化为 VTK 对象，允许在 Paraview 上表示地质模型，以进行自定义切片、过滤、透明度和样式设置。
在这一部分中，我们导入教程所需的库。该脚本需要 Gempy 以及 Numpy 和 Matplotlib。我们在脚本单元（%matplotlib inline）之后为 Matplotlib 图形的交互式表示配置了一个 Jupyter 选项。请注意，警告只是用户在运行脚本时必须记住的消息，它们并不意味着代码失败。由于本教程是在 Windows 上，一些补充库无法安装，但地质建模代码的整体性能是完整的。
本教程在 2km x 2km x 2km 的扩展范围内创建一个 100 列 x 100 行 x 100 层的网格。更高的分辨率是可能的，但计算时间会更长。坐标系为局部坐标系，教程将使用 UTM 坐标评估 Gempy 的性能。方向和地质接触从 CSV 文件导入并转换为 Pandas 数据框。然后定义地质系列（断层/地层）以及地质构造序列。值得一提的是，必须独立插入故障，其中最近的排在第一个：
Gempy 有一些有用的特征来表示定义的地质序列和地层序列。
可以通过 Gempy 的“.get_”函数访问用于构建地质模型的不同数据集。
在这一部分中，进行了 2D 和 3D 表示以呈现接口和方向。
输入数据准备好后，我们可以使用 Gempy 库中的 InterpolatonData 方法定义插值的数据和参数。
地质模型在“compute_model”方法下计算。模型过程的结果是与 geo_data 具有相同阵列维度的岩性和断层。
岩性区块有两部分，第一部分包含有关岩性地层的信息，第二部分表示方向。在这一部分中，岩性分布和断层分离信息用直方图表示。
与任何其他 Numpy 数组一样，生成的岩性块可以在 Matplotlib 上表示。然而，Gempy 有特殊的横截面表示方法。通过使用 Jupyter 小部件，沿 Y 方向的地质横截面的交互式表示是使用把手沿行移动来执行的。