An approach to solving the optimization problem for a multistage hydraulic fracturing (MHF) design is proposed. Free optimization parameters are set as follows: the length of the horizontal well, the number of fractures and the amount of proppant loaded in each fracture. Optimization targets are the maximum of the cumulative well production, the maximum of the net present value (NPV), and the minimum of fracturing costs. As an optimization algorithm, we use genetic algorithm NGSA-II, which requires calculating three related values at each iteration step: the fracture geometry, the post-fracture well production, and economy indicators. The approach proposed is illustrated by the case of the lowpermeability oil reservoir under the following suggestions. It is assumed that the oil reservoir is rectangular, the horizontal well is positioned along the centerline of the reservoir, and hydraulic fractures are placed on equal distance perpendicular to the wellbore and symmetric about it. In addition, all fractures are identical to each other. The geometric characteristics of fractures (length and width) are determined by the amount of proppant injected and are calculated by empirical relationships. To obtain the value of post-fracture well production, approximate analytical formulas that take into account the final conductivity of fractures are applied. The main economic indicator that characterizes the economic profitability of the MHF is the NPV-based income. The case studies for different values of the average permeability of the reservoir are carried out. The numerical results show that the dependence of NPV on the well production is not always linear. The results show that after certain adjustment of the algorithm to parameters of the particular field, the model can be used as a tool for planning of the field development.