Directional drilling 3D visualization modeling about gempy

[xiaohua00921]

Unnamed: 0,0,x,y,name,num,z,year,7,8,9,altitude,base,formation,top,_top_abs,md
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,1035.5970765000002,topo,0.0,108.7132874,0.0
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,1652.2456364,etchegoin,1035.5970765000002,1144.3103639,1035.5970765000002
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,2109.9936834,macoma,1652.2456364,1760.9589238,1652.2456364
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,2799.0279844,fruitvale,2109.9848934,2218.6981808,2109.9848934
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,2110.0036834,mclure,2109.9936834,2218.7069708,2109.9936834
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,3062.3277894,round_mountain,2799.0279844,2907.7412718,2799.0279844
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,3162.6180724,olcese,3062.3277894,3171.0410768,3062.3277894
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,3323.9220274,vedder,3162.6180724,3271.3313598,3162.6180724
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,4352.0370174,,3323.9220274,3432.6353148,3323.9220274
1,2906949,312385,3917847,KCL_B45,45,-3107,XX,14,0,0.275,121.8172989,833.7528457999999,topo,0.0,121.8172989,0.0
This is a sample data from the gempy website.
Listen to the people on gempy's github explaining
This example data is probably provided by you.
So I'm here to ask a few questions
I would appreciate your reply.
First of all, this example data
I would like to ask that the xyz value is not even changed.
But the z-values of the points shown in the model have changed.
I'm guessing that this is due to the fact that the source code changes the z-values according to the depth.
I'm not sure if that's correct.
Secondly, what I'm really trying to do is
I want to do a directional drill on the data
That is, it's not that the xy values don't change.
It's that it's all changing.
I see other examples on gempy that require inclination and azimuth.
I think I have all this data, including a direct xyz value.
And then based on the lithology at each point of the trajectory.
Interpolate the lithology in space
But I don't know how to give this data in the format I have.
Can you give me some help?
Or better yet, can you provide me with a sample data that meets my needs?

[xiaohua00921]

Also, I'd like to know about this interpolation algorithm that interpolates the lithologies in space based on the different lithologies, what is this interpolation algorithm?
Where can I see the source code of the algorithm?
I will show you the code that you use in this, but I don't know what the algorithm is.

调用 add_orientations 函数

gp.add_orientations(
x=x,
y=y,
z=z,
elements_names=elements_names,
pole_vector=pole_vector,
geo_model=geo_model
)

geo_model.interpolation_options

%%

gp.compute_model(geo_model)

g3d = gpv.plot_3d(geo_model, show_lith=False, show=False)
g3d.p.add_mesh(well_mesh)
g3d.p.show()

[AlexanderJuestel]

@xiaohua00921 literature and the source code are mentioned in gempy-project/gempy#973 (comment), so please avoid duplicate questions.

Also, I'd like to know about this interpolation algorithm that interpolates the lithologies in space based on the different lithologies, what is this interpolation algorithm? Where can I see the source code of the algorithm? I will show you the code that you use in this, but I don't know what the algorithm is.

调用 add_orientations 函数

gp.add_orientations( x=x, y=y, z=z, elements_names=elements_names, pole_vector=pole_vector, geo_model=geo_model )

geo_model.interpolation_options

%%

gp.compute_model(geo_model)

g3d = gpv.plot_3d(geo_model, show_lith=False, show=False) g3d.p.add_mesh(well_mesh) g3d.p.show()

[AlexanderJuestel]

@xiaohua00921

Could you please provide a link to the example you are referring to?

What do you mean by the xyz values are not changed? These would be the input data for your modelling. Why would they change?

You have to provide two csv files. One for the interfaces points with x,y,z, and formation column and one for the orientations with x,y,z, formation, dip, azimuth, polarity columns.

But I don't know how to give this data in the format I have. Can you give me some help? Or better yet, can you provide me with a sample data that meets my needs?

[xiaohua00921]

@xiaohua00921

您能否提供您所引用的示例的链接？

xyz 值不改变是什么意思？这些将是建模的输入数据。他们为什么要改变？

您必须提供两个 csv 文件。一个用于 x，y，z 和形成列的界面点，另一个用于 x，y，z、形成、倾角、方位角、极性列的方向。

但我不知道如何以我拥有的格式提供这些数据。您能帮我一些忙吗？或者更好的是，您能否提供满足我需求的示例数据？
Hello, I will explain my question to you in detail
First of all this problem all comes from this example of yours
The name of the example is “Unknown Model: Importing Borehole Data and Building a 3D Geological Model with GemPy”.
The data in this example is the same data I provided you with above.
I'm going to provide it again for you here

Unnamed: 0,0,x,y,name,num,z,year,7,8,9,altitude,base,formation,top,_top_abs,md
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,1035.5970765000002,topo,0.0,108.7132874,0.0
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,1652.2456364,etchegoin,1035.5970765000002,1144.3103639,1035.5970765000002
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,2109.9936834,macoma,1652.2456364,1760.9589238,1652.2456364
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,2799.0279844,fruitvale,2109.9848934,2218.6981808,2109.9848934
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,2110.0036834,mclure,2109.9936834,2218.7069708,2109.9936834
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,3062.3277894,round_mountain,2799.0279844,2907.7412718,2799.0279844
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,3162.6180724,olcese,3062.3277894,3171.0410768,3062.3277894
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,3323.9220274,vedder,3162.6180724,3271.3313598,3162.6180724
0,2904240,303412,3913997,1-1,12,-3894,XX,14,0,0.409,108.7132874,4352.0370174,,3323.9220274,3432.6353148,3323.9220274
1,2906949,312385,3917847,KCL_B45,45,-3107,XX,14,0,0.275,121.8172989,833.7528457999999,topo,0.0,121.8172989,0.0
The 3D model generated by running this example source code looks like this
I'll give you a screenshot of it below
I'm thinking that the xyz does not change from one set of points to the next.
It's only the depth of each rock type that changes.
In this xyz-axis 3D space.
Shouldn't it be just one point?
Because the z-value doesn't change.
What's going on here?

[xiaohua00921]

@xiaohua00921
您能否提供您所引用的示例的链接？
xyz 值不改变是什么意思？这些将是建模的输入数据。他们为什么要改变？
您必须提供两个 csv 文件。一个用于 x，y，z 和形成列的界面点，另一个用于 x，y，z、形成、倾角、方位角、极性列的方向。

但我不知道如何以我拥有的格式提供这些数据。您能帮我一些忙吗？或者更好的是，您能否提供满足我需求的示例数据？您好，我将向您详细解释我的问题 首先，这个问题都来自您的这个示例，示例的名称是“未知模型：导入钻孔数据并使用 GemPy 构建 3D 地质模型”。此示例中的数据与我在上面提供的数据相同。我将在这里再次为你提供它

未命名： 0,0，x，y，名称，数量，z，年，7,8,9，海拔，基地，地层，顶部，_top_abs，md 0,2904240,303412,3913997,1-1,12，-3894，XX，14,0,0.409,108.7132874,1035.5970765000002，topo，0.0,108.7132874,0.0 0,2904240,303412,3913997,1-1,12，-3894，XX，14,0,0.409,108.7132874,1652.2456364，etchegoin，1035.5970765000002,1144.3103639,1035.5970765000002 0,2904240,303412,3913997,1-1,12，-3894，XX，14,0,0.409,108.7132874,2109.9936834，macoma，1652.2456364,1760.9589238,1652.2456364 0,2904240,303412,3913997,1-1,12，-3894，XX，14,0,0.409,108.7132874,2799.0279844，Fruitvale，2109.9848934,2218.6981808,2109.9848934 0,2904240,303412,3913997,1-1,12，-3894，XX，14,0,0.409,108.7132874,2110.0036834，mclure，2109.9936834,2218.7069708,2109.9936834 0,2904240,303412,3913997,1-1,12，-3894，XX，14,0,0.409,108.7132874,3062.3277894，round_mountain，2799.0279844,2907.7412718,2799.0279844 0,2904240,303412,3913997,1-1,12，-3894，XX，14,0,0.409,108.7132874,3162.6180724，olcese，3062.3277894,3171.0410768,3062.3277894 0,2904240,303412,3913997,1-1,12，-3894，XX，14,0,0.409,108.7132874,3323.9220274，vedder，3162.6180724,3271.3313598,3162.6180724 0,2904240,303412,3913997,1-1,12，-3894，XX，14,0,0.409,108.7132874,4352.0370174，，3323.9220274,3432.6353148,3323.92202741,2906949,312385,3917847，KCL_B45,45，-3107，XX，14,0,0.275,121.8172989,833.7528457999999，topo，0.0,121.8172989,0.0 运行此示例源码生成的 3D 模型如下所示，我将在下面给大家一个屏幕截图，我认为 xyz 不会从一组点更改为下一组点。只是每种岩石类型的深度会发生变化。在此 xyz 轴 3D 空间中。难道不应该只是一个点吗？因为 z 值不会改变。 这是怎么回事？

Don't I think the xyz data is at the point where the hole is drilled vertically?
According to the 3D model that was run
Shouldn't his z-value vary?
I don't know if you're clear about my question.

[xiaohua00921]

I will provide you with links to examples.
https://docs.gempy.org/examples/real/mik.html#sphx-glr-examples-real-mik-py
This is the sample website link.