Numerical local upscaling of elastic geomechanical properties for heterogeneous continua