A method for obtaining geological heterogeneity trends of a geological formation, including the steps: drilling wells (w1, w2, w3, w4, w5) that penetrate the formation, acquiring well logs (log 1, log 2) for each well (w1, w2, w3, w4, w5) as function of depth inter Vals (Dk) of the respective well, determining a third degree tensor (Tk, m, n), where a z-dimension denotes the depths, a x-dimension denotes the well logs, and a y-dimension denotes the wells, extracting matrices (L1k,m, L2k,m, . . . , LMk,n) from the tensor (Tk, m, n), clustering the matrices (L1k,n, L2k,n, . . . , LMk,n) based on the characteristics of the corresponding well logs (log 1, log 2) to a clustering result matrix, aggregating the clustering result matrix to a cluster ensemble (π1, π2, . . . , πM), and spatial partitioning the cluster ensemble (π1, π2, . . . , πM) to a map that shows the geological heterogeneity trends associated with cluster types of the wells (w1, w2, w3, w4, w5).