|
83 | 83 | # %%
|
84 | 84 | geo_model.orientations_copy
|
85 | 85 |
|
86 |
| -# %% |
87 |
| -# Declaring the Sequential Order of Geological Formations |
88 |
| -# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
89 |
| -# In our model, we want the geological units to appear in the correct chronological order. |
90 |
| -# Such order could be determined by a sequence of stratigraphic deposition, unconformities |
91 |
| -# due to erosion, or other lithological genesis events like igneous intrusions. A similar |
92 |
| -# age-related order is declared for faults in our model. In GemPy, we use the function |
93 |
| -# `gempy.map_stack_to_surfaces` to assign formations or faults to different sequential series |
94 |
| -# by declaring them in a Python dictionary. |
95 |
| -# |
96 |
| -# The correct ordering of series is crucial for model construction! It's possible to assign |
97 |
| -# several surfaces to one series. The order of units within a series only affects the color |
98 |
| -# code, so we recommend maintaining consistency. The order can be defined by simply changing |
99 |
| -# the order of the lists within `gempy.core.data.StructuralFrame.structural_groups` and |
100 |
| -# `gempy.core.data.StructuralGroups.elements` attributes. |
101 |
| -# |
102 |
| -# Faults are treated as independent groups and must be younger than the groups they affect. |
103 |
| -# The relative order between different faults defines their tectonic relationship |
104 |
| -# (the first entry is the youngest). |
105 |
| -# |
106 |
| -# For a model with simple sequential stratigraphy, all layer formations can be assigned to |
107 |
| -# one series without an issue. All unit boundaries and their order would then be determined |
108 |
| -# by interface points. However, to model more complex lithostratigraphical relations and |
109 |
| -# interactions, separate series definition becomes important. For example, modeling an |
110 |
| -# unconformity or an intrusion that disrupts older stratigraphy would require declaring a |
111 |
| -# "newer" series. |
112 |
| -# |
113 |
| -# By default, we create a simple sequence inferred from the data: |
114 |
| -# |
115 |
| -# |
116 |
| -# Setting Up the Structural Framework |
117 |
| -# """"""""""""""""""""""""""""""""""" |
118 |
| -# |
119 |
| -# Each data entry includes an ID that corresponds to a formation name, which GemPy uses to create a base structural framework. |
120 |
| -# However, the sequence of these formations is still arbitrary, as they have all automatically been assigned to the default |
121 |
| -# structural group in our `geo_model` object. We will fix this in the subsequent steps. |
122 | 86 |
|
123 | 87 | # %%
|
124 |
| -geo_model.structural_frame |
125 |
| - |
126 |
| -# %% |
127 |
| -# Our example model comprises four main layers (plus an underlying |
128 |
| -# basement that is automatically generated by GemPy) and one main normal |
129 |
| -# fault displacing those layers. Assuming a simple stratigraphy where each |
130 |
| -# younger unit was deposited onto the underlying older one, we can assign |
131 |
| -# these layer formations to one series called "Strat\_Series". For the |
132 |
| -# fault, we declare a respective "Fault\_Series" as the first key entry in |
133 |
| -# the mapping dictionary. We could give any other names to these |
134 |
| -# series, the formations however have to be referred to as named in the |
135 |
| -# input data. |
136 |
| -# |
137 |
| -# |
138 | 88 | # Declaring the Sequential Order of Structural Elements (Geological Formations)
|
139 | 89 | # """""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""
|
140 | 90 | # In our model, we want the geological units to appear in the correct chronological order. This order could be determined by a sequence
|
141 | 91 | # of stratigraphic deposition, unconformities due to erosion, or other lithological genesis events like igneous intrusions. A similar
|
142 | 92 | # age-related order is declared for faults in our model. In GemPy, we use the function `gempy.map_stack_to_surfaces` to assign formations
|
143 | 93 | # or faults to different sequential series by declaring them in a Python dictionary.
|
144 |
| - |
| 94 | +# |
145 | 95 | #
|
146 | 96 | # The correct ordering of series is crucial for model construction! It's possible to assign several surfaces to one series. The order of
|
147 | 97 | # units within a series only affects the color code, so we recommend maintaining consistency. The order can be defined by simply changing
|
|
179 | 129 | # Note how the structural frame still indicates the "Fault Series" group to have a relation type "erode".
|
180 | 130 | # We still need to tell GemPy that we want this group to be a fault. We do this using the function `set_is_fault`.
|
181 | 131 |
|
182 |
| -geo_model.structural_frame # Display the resulting structural frame |
183 |
| - |
184 | 132 | # %%
|
185 | 133 | gp.set_is_fault(
|
186 | 134 | frame=geo_model.structural_frame,
|
|
0 commit comments