3.3.7. Model (Sensitivity) Weights FileΒΆ

The weighting matrix is generally used to implement a sensitivity weighting which may be needed when there are few source locations or few frequencies; that is, when the data do not adequately constrain the distribution of conductivity with depth. In these circumstances, the model recovered without any weighting may place structures at the surface. The weighting matrix can be used to find solutions with alternative distributions of conductivity with depth. The methodology is similar to that employed when inverting surface gravity or magnetic data, so the weights in the matrix can be inversely proportional to some power of the depth of the associated cell.

The structure and ordering of the weights file is identical to that of a model file, i.e.:


\(w_{1,1,1}\)
\(w_{1,1,2}\)
\(\;\vdots\)
\(w_{1,1,nz}\)
\(w_{1,2,1}\)
\(\;\vdots\)
\(w_{1,nx,nz}\)
\(w_{2,1,1}\)
\(\;\vdots\)
\(w_{i,j,k}\)
\(\;\vdots\)
\(w_{ny,nx,nz}\)


where \(N=nx \times ny \times nz\) is the total number of cells in the corresponding tensor mesh. Cells in the x and y directions are organized E to W and S to N, respectively, while cells are organized top to bottom in the vertical. Thus \(w_{1,1,1}\) is the weight for the top southwesternmost cell and \(w_{ny,nx,nz}\) is the weight for the bottom northeastermost cell. For sensitivity weighting, only relative values matter. Regions where cells have smaller weights will have more structure. For example, assigning weights of surface cells to a value of unity, and deeper cells to values between 0 and 1 will force structure to occur at greater depths.

Note

Weights for inactive cells, including the air, do not play a role in the solution and so they can be assigned any number.