Most analyses of geological heterogeneity are based on the assumption of repetitive deposition and stationary transition probabilities. However

the effects of tectonism and weathering are ignored

which leads to the absence of consistency verification of geological heterogeneity features（GHFs）. A test method for the consistency of GHFs is herein proposed

as well as a partition method to detect stationary Markovian zones. The correlation coefficients between local transition probability matrices are adopted as quantitative criteria of consistency. The stratigraphic distribution is simulated by the coupled Markov chain model

and the deformation analysis of surrounding rock uncertainty is accomplished through Monte-Carlo simulations. Taking advantage of the borehole data in Tsingtao

the effect on stratigraphic distribution and uncertainty analyses induced by the consistency of GHFs is discussed

and the effectiveness of the proposed method is illustrated. The results indicate that the accuracy of uncertainty analyses cannot be improved with increased number of boreholes when there are conflicting GHFs. After the consistency verification of GHFs and the partition of stationary Markovian zones

the absence of local features can be avoided to improve the accuracy of uncertainty analyses of surrounding rock deformation.