A method of predicting model parameters (R1, R2, R3, . . . ) of a geological formation under investigation, wherein said geological formation is distinguished by reservoir parameters including observable data parameters and the model parameters (R1, R2, R3, . . . ) to be predicted, comprises the steps of calculating at least one model constraint (M1, M2, M3, . . . ) of the model parameters (R1, R2, R3, . . . ) by applying at least one rock physics model (f1, f2, f3, . . . ) on the model parameters (R1, R2, R3, . . . ), said at least one model constraint (M1, M2, M3, . . . ) including modelled data of at least one of the data parameters, and applying an inverse model solver process on observed input data (d1, d2, d3, . . . ) of at least one of the data parameters, including calculating predicted model parameters, which comprise values of the model parameters (R1, R2, R3, . . . ) which give a mutual matching of the input data and the modelled data, wherein the modelled data are provided with probability distribution functions, the inverse model solver process is conducted based on the probability distribution functions, wherein multiple predicted values of the model parameters are obtained comprising values of the model parameters (R1, R2, R3, . . . ) which give the mutual matching of the input data and the modelled data, and model probabilities of the predicted model parameters are calculated in dependency on the probability distribution functions.