This study investigates rock drillability characteristics and spatial distribution patterns to optimize drilling efficiency. Laboratory testing of core samples from six wells provided comprehensive mechanical data, including uniaxial compressive strength (147.81-306.82 MPa), triaxial compressive strength (310.56-785.00 MPa), and drillability indices (6.86-10.46). Multiple linear regression models correlating rock properties with well log parameters achieved high predictive accuracy (R² = 0.886-0.994). Spatial interpolation using kriging and radial basis function methods revealed that 75% of samples classify as difficult-to-drill formations. Kriging demonstrated superior performance with >90% prediction accuracy and significantly lower errors for triaxial strength prediction (5.6% vs. 19.2% mean absolute error). Strong correlations between rock strength and drillability (r > 0.8) confirm mechanical strength as the primary controlling factor. This methodology successfully characterizes three-dimensional drillability distribution, providing quantitative foundations for drilling optimization in complex continental shale systems.