Recent enhancements in computational capacity provide an opportunity for harnessing the enormous amount of reservoir data already acquired and extracting useful information for hydrocarbon exploration, development, and production. This article reports a three-step clustering technique to determine well groups based on subsurface geological heterogeneity using feature extraction, hierarchical ensemble clustering, and spatial mapping. The first step of the presented methodology is to group the wells into different clusters based on the formation rock composition and property features extracted from well logs using the expectation maximization algorithm. The one-dimensional (1D) stacking pattern of each well log curve is expressed through a two-dimensional (2D) transformation scheme. Thus, the clustering can capture the vertical stacking patterns of well logs, which is essential for reservoir heterogeneity characterization. This base clustering process generated a feature matrix which is further grouped through the hierarchical ensemble clustering in a latent space of well logs in the second step. Through the ensemble clustering, different clustering proposals obtained from the base clustering are integrated corroboratively to reflect a comprehensive feature of all studied logs. In the third step, the spatial clustering is performed based on the ensemble results, considering the spatial distances between well locations in the target area. The results of the 2D spatial map may provide insights into the sedimentary depositional environment in terms of the lateral geological heterogeneity features. Therefore, the proposed clustering technique can present a fast geological modeling method to integrate geological heterogeneity features presented in multiple well logs, which is not yet fully utilized in traditional geomodeling approaches. The results can also support further reservoir studies, such as petrophysical modeling, reservoir modeling, and fluid flow simulation studies.