Examples of computer-implemented method for geostatistical reservoir modeling include: obtaining a prior probability distribution function using primary data; obtaining a likelihood probability distribution function, via a computer processor, using secondary data, wherein the likelihood probability distribution function is obtained using a Gaussian mixture model that models non-linear relationship between the primary data and secondary data; combining the prior probability distribution function with the likelihood probability distribution function to generate a posterior probability distribution function; and outputting a reservoir model based on the posterior probability distribution function.