3.2.6. Cell and Face Weights¶
3.2.6.1. Cell Weights Input File¶
The parameters used to create cell weights are defined in the input file. The lines within the input file are as follows:
Line # |
Parameter |
Description |
|---|---|---|
1 |
path to octree mesh |
|
2 |
path to active cells model |
|
3 |
set number of cells below surface in which weighting horizontal smoothing is applied |
|
4 |
surface weights values |
|
5 |
name for output face weighting file |
3.2.6.1.1. Line Descriptions¶
OcTree Mesh: file path to the OcTree mesh file
Active Topography Cells: Here, the user can choose to specify the cells which lie below the surface topography. To do this, the user may supply the file path to an active cells model file or type “ALL_ACTIVE”. The active cells model has values 1 for cells lying below the surface topography and values 0 for cells lying above.
# Surface Layers: This represents the first parameter used to apply surface weighting. Here, the user specifies how many cell layers below the surface will have surface weighting in X and Y.
Horizontal Surface Weights: Here, the user specifies the weights on X and Y faces for every layer (from surface layer downwards). Essentially, we are invoking a smoothness along the X and Y directions that decreases with depth. The user must enter a set of decreasing values separated by spaces. The number of values that must be entered is equal to the integer value set on the previous line.
Output Name: File name for the output interface weights file.
3.2.6.2. Interface Weights Input File¶
The parameters used to create interface weights are defined in the input file. The lines within the input file are as follows:
Line # |
Parameter |
Description |
|---|---|---|
1 |
path to octree mesh |
|
2 |
path to active cells model |
|
3 |
path to a reference conductivity model |
|
4 |
set as log or linear mapping |
|
5 |
set threshold for largest gradients preserved in recovered model |
|
6 |
set number of cells below surface in which weighting horizontal smoothing is applied |
|
7 |
sets the weighting for horizontal smoothing for each layer |
|
8 |
name for output face weighting file |
Fig. 3.7 Example input file for creating interface weights (Download ).¶
3.2.6.2.1. Line Descriptions¶
OcTree Mesh: file path to the OcTree mesh file
Active Topography Cells: Here, the user can choose to specify the cells which lie below the surface topography. To do this, the user may supply the file path to an active cells model file or type “ALL_ACTIVE”. The active cells model has values 1 for cells lying below the surface topography and values 0 for cells lying above.
Conductivity Model: Here, the user provides reference model containing the conductivity structures they would like to preserve. The flag “NO_MODEL” may be used if only surface weighting is applied.
Mapping Type: Here, the user specifies whether the mapping between the model value and the physical property value on the mesh is linear or logarithmic. To specify, use the flags “LIN_MODEL” or “LOG_MODEL”. In the case of dcip, the inversion recovers the log-conductivity. As a result, the user should use “LOG_MODEL”.
Gradient Tolerance: Here, the user specifies the threshold for the largest gradients that can be preserved followed by a replacement value if the gradient is over the threshold (i.e. Val1 Val2). If we try to preserve very large gradients due to blocky models, we may fit these features at the expense of others. The tolerance threshold and replacement value should be determined from the expected physical property values and cell size. For log-conductivity and 100m cell widths, a good threshold might be: \(\sigma_{max}/dh\).
# Surface Layers: This represents the first parameter used to apply surface weighting. Here, the user specifies how many cell layers below the surface will have surface weighting in X and Y.
Horizontal Surface Weights: Here, the user specifies the weights on X and Y faces for every layer (from surface layer downwards). Essentially, we are invoking a smoothness along the X and Y directions that decreases with depth. The user must enter a set of decreasing values separated by spaces. The number of values that must be entered is equal to the integer value set on the previous line.
Output Name: File name for the output interface weights file.