Geophysics and Geochemistry - Structure Calculation from Gravity Data and Density Logs

Combination of gravity data and density logs makes it possible to determine the structural relief. Under a wide class of geologic conditions, the solution is unique, or limited to a single parametric system of solutions. As a useful byproduct, the method yields a new way to construct derivative maps, continue the field up or down, subtract regional anomalies, and perform smoothing. In the past, the problem of interpreting gravity anomalies in terms of structural geology was solved along two main lines: 1) A body of a certain shape and density contrast was assumed and the corresponding gravity anomaly calculated, analytically for simple bodies, graphically and/or numerically for more complex ones; then the computed anomaly was compared with the observed one, and the assumed shape corrected according to the differences between the two. 2) The anomaly is assumed to come from a warped interface between two media of different density; the relief of the warping is assumed small compared to the depth of burial, so that the structure can be approximated by a flat horizontal surface of variable density distribution. The first method is a standard technique, described in every text book on applied geophysics, and is well applicable, for bodies with predominant vertical dimension. The second one, probably first used by Tsuboi and Fuchida,' is mathematically more elegant. It applies to the fashionable field of integral transforms (Bullard and Cooper, 2 Tomoda and Aki, 3 Dean 4) and is useful mainly in cases when the vertical relief of the structures is comparatively small. 30th methods, tacitly or openly, assume a density discontinuity across the boundary of the two media, while within each medium the density is constant. However, the recent invention of the gamma-gamma