Determination of the geometry of bedrock, by nonlinear inverse modeling of gravity data, is the aim of this paper. In this method, reliable geological structures can be obtained by minimum geology priori information. The usual practice of inverting gravity anomalies of two-dimensional bodies replaced by n-sides polygon for determining location of the vertical that best explain the observed anomalies. In this method, the geometry of the bedrock is replaced by a series of juxtaposing prisms. Finally the length of each prism is the depth of the bedrock at that point.
The algorithm uses a nonlinear iterative procedure for simulation of bedrock geometry. At the first step, the nonlinear problem changes to a linear problem by a proper approximation and standard method. The second step is the parameterization of the model. Finally, an initial model is suggested on the basis of geological and geophysical assumption and using the numerical analysis, Jacobean matrix is calculated. In each iteration the inversion will improve the initial model, considering the differences between observed and calculated gravity anomalies, based on Levenberg-Marquardt's method.
The practical effectiveness of this method is demonstrated by inversion of synthetic (free noise and noise contaminated data) and real examples. The real data is acquired over the Moghan area and the results compared with the geological information.