该方法可以创建一个正则网格,并将正则网格传递给模型来进行建模和仿真。
extent
(list):网格的空间范围,包括模型中的所有地表和地下区域的最小和最大边界 [xmin, xmax, ymin, ymax, zmin, zmax]
。
resolution
(list):网格的分辨率 ,即该网格中的每个单元的体积大小,它是一个具有三个正整数的列表 [dx, dy, dz]
。
origin
(list):网格原点的位置,它是一个具有三个正整数的列表 [x,y,z]
。
grid
(gempy.core.grid.RegularGrid
):正则网格对象。import gempy as gp
# 创建一个空的gempy项目
project = gp.create_project('my_project')
# 创建一个正则网格
extent = [-60, 60, -60, 60, -150, 0]
resolution = [25, 25, 25]
origin = [0, 0, 0]
regular_grid = project.set_regular_grid(extent=extent, resolution=resolution, origin=origin)
# 在正则网格上添加GeoDataFrame
geo_data = gp.create_data('my_data', extent=extent)
geo_data.set_coordinates_from_grid(regular_grid)
project.add_data(geo_data)
# 定义地层系列并添加到gempy项目中
series = project.add_series('my_series')
series.add_surfaces(['surface1', 'surface2', 'surface3'])
series.set_is_fault(['surface3'], is_fault=True)
# 构建gempy模型
project.add_orientations()
project.update_to_interpolator()