Gravity anomaly interpretation of 2D fault morphologies by means of nonplanar fault planes and exponential density contrast model: a space domain technique

Boundary faults associated with thick sedimentary basins are more often curved in cross section rather than planar. We develop a space domain-based automatic gravity inversion technique to quantify such listric fault sources from a set of observed gravity anomalies. The density contrast within the hanging wall of fault morphology is presumed to be known according to a prescribed exponential law. Furthermore, the fault plane is described by a polynomial function of arbitrary but specific degree, whose coefficients become the unknown parameters to be estimated from a set of observed gravity anomalies in addition to the thickness of the fault structure. Using a set of characteristic anomalies, the present inversion identifies approximate parameters pertaining to the origin of fault plane and depth to decollement horizon. Based on the errors between the observed and model gravity anomalies of the structure, the algorithm constructs and solves a system of normal equations to estimate the improvements in depth and coefficients of the polynomial in an iterative approach until one of the specified convergence criteria is fulfilled. The efficacy of the algorithm is shown with the analysis of gravity anomalies attributable to a synthetic model of a listric fault source in the presence of pseudorandom noise. Application of the proposed inversion technique on the observed gravity anomalies of the Ahri-Cherla master fault of the Godavari subbasin in India using the derived exponential density contrast model has yielded an interpretation that is consistent with the available/reported information.