1. A computer-implemented method for automated rock physics modeling, comprising: (a) pre-processing log data obtained from a suite of logging tools and corresponding auxiliary information over a range of depth of a subterranean formation, the pre-processing including; (i) dividing the log data and the corresponding auxiliary information into zones within the range of depth, (ii) coarsening the log data and the corresponding auxiliary information by determining an average value for the log data and the corresponding auxiliary information for each zone to produce a coarsened log data and a coarsened corresponding auxiliary information over the range of depth; (b) solving an inverse problem associated with the rock physics model with a method comprising, (i) identifying a list of parameters associated with the range of depth within the subterranean formation, (ii) dividing the list of parameters into a set of at least two subgroups of parameters, (iii) perform a local inversion for each of the at least two subgroups of parameters identified in step (b)(ii) to obtain optimum values for all parameters belonging to a given subgroup, (iv) collect the optimal values of parameters obtained in step (iii) to determine the optimal values for the list of parameters identified in step (b)(i) by comparing the local inversions from step (b)(iii) with the pre-processed log data from step (a) to identify the subgroup of parameters as the optimum set of parameters for use in a forward run of the rock physics model; and (c) computing a set of rock mechanical properties, comprising (i) a forward run of the rock physics model using the optimal parameters obtained in step (b)(iv) to calculate at least one of density log, a compressional velocity log and a shear velocity log, (ii) calculating an optimized set of rock mechanical properties for the range of depth use the at least one of the calculated density, compressional velocity and shear velocity logs of step (c)(i) in a rock physics model.