In KAERI, a site descriptive model for stress field estimation had already been constructed by using integrated field data within KURT site scale. A sub-divided rock block domain containing major fracture zones has spatial rock mass and fault properties. The properties were decided based on the rock classification results of several borehole investigations. Modeling for maximum and minimum horizontal stress field estimation was performed and compared with the in-situ data. As a result, a depth-dependent stress ratio was adopted to obtain numerical results closer to actual in-situ data. Although the results were suitable at a relatively low depth (~500 m), there is still some deviation trend at a deep depth. This study aims to improve these modeling results by incorporating not only depth-dependent stress ratio but also changes in rock mass properties along the depth. The deep borehole of DB2 in the KURT site indicated fracture distribution corresponding to the property changes. Natural fractures are typically randomly oriented, and the fracture frequency decreases with increasing depth. The increase in P-wave velocity log data accompanies these features. A discrete fracture network (DFN) model can be used to simulate fractured rock explicitly, but DFN modeling is not feasible for site scale analysis because of its numerical efficiency. Therefore, as a preliminary model in this study, the effect of fracture distribution was considered by substituting the influence for the depth-dependent property. The properties were estimated from the fracture frequency and P-wave velocity log data. The influence of elastic modulus and density on the site stress field was dominant, with decreasing the deviation trend between modeling and in-situ data at a deep depth. Considering that the depth of the repository construction is within about 500 m, it may not be necessary to consider the change of rock properties with depth. However, it was determined that the rock property effect might need to be considered when the loading conditions change due to subsidence in the long-term evolution scenario. Continuously, this site descriptive modeling will be interdependently conducted with a representative DFN block model for deriving equivalent properties in fractured rock.