The extensive depletion of the development target has triggered the demand for infill drilling in the upside target of multilayer unconventional reservoirs. To optimize the hydraulic fracturing design of newly drilled wells, we need to investigate the stress changes in the upside target induced by parent-well production. In this work, an integrated parent-child workflow is presented to model the spatial-temporal stress evolution and propose the optimal development strategy for the upside target using a data set from the Permian Basin. The stress dependence of matrix permeability and fracture conductivity is determined based on available experimental data and incorporated in our reservoir simulation with the aid of an embedded discrete fracture model (EDFM). With calibrated reservoir properties from history matching of an actual well in the development target (i.e., 3rd BS Sand), we run the finite element method (FEM) based geomechanics simulator to predict the 3D spatial-temporal evolution of the local principal stresses. A displacement discontinuity method (DDM) hydraulic fracture model is then applied to simulate the multi-cluster fracture propagation in the upside target (i.e., L2BSSh) with the updated heterogeneous stress field. Numerical results indicate that stress field redistribution associated with parent-well production not only occurs within the development target but also vertically propagates to the upside target. A smaller parent-child horizontal offset induces a severer deviation of child-fractures towards the parent wellbore, resulting in more substantial well interference and less desirable oil and gas production. The parent-child fracture overlapping ratio in our study is in 0.6 ~ 0.8 for the 400 ft horizontal offset and 0.2 ~ 0.5 for the 600 ft horizontal offset. Varying the parent-child vertical offset gives the same trend as we change the horizontal offset. But with a delayed infill time, placing child-well in different layers causes more significant variation in the ultimate recovery. Moreover, infill operations at an earlier time are less affected by parent-well depletion because of the more homogeneous stress state. The candidate locations to implement infill-wells are suggested in the end for different infill timing by co-simulation of the parent-child production. With the reservoir-geomechanics-fracture model, this work provides a general workflow to optimize the child-well completion in multilayer unconventional reservoirs. The conclusions drawn from this study are of guiding significance to the subsequent development in the Permian Basin.