Multiple-factor analysis and optimization play a critical role in the the ability to maximizethe stimulated reservoir volume (SRV) and the success of economic shale gas production. In this paper, taking the typical continental naturally fractured silty laminae shale in China as anexample, response surface methodology (RSM) was employed to optimize multiple hydraulic fracturing parameters to maximize the stimulated area in combination with numerical modeling based on the coupled flow-stress-damage (FSD) approach. This paper demonstrates hydraulic fracturing effectiveness by defining two indicesnamelythe stimulated reservoir area (SRA) and stimulated silty laminae area (SLA). Seven uncertain parameters, such as laminae thickness, spacing, dip angle, cohesion, internal friction angle (IFA), in situ stress difference (SD), and an operational parameter-injection rate (IR) with a reasonable range based on silty Laminae Shale, Southeastern Ordos Basin, are used to fit a response of SRA and SLA as the objective function, and finally identity the optimum design under the parameters based on simultaneously maximizingSRA and SLA. In addition, asensitivity analysis of the influential factors is conducted for SRA and SLA. The aim of the study is to improve the artificial ability to control the fracturing network by means of multi-parameteroptimization. This work promises to provide insights into the effective exploitation of unconventional shale gas reservoirs via optimization of the fracturing design for continental shale, Southeastern Ordos Basin, China.