In this paper the main criteria of the water curtain system for unlined rock caverns (URCs) is described. By the application of numerical modeling and genetic programming (GP), a method for water curtain system pre-design for Iranian crude oil storage URCs (common dimension worldwide) is presented. A comprehensive set of numerical simulations is performed using the finite element based commercial software (COMSOL 5.1) where the results are used as database for genetic programming. Describing equations of water inflow to the filled and empty caverns and water production rate of water curtain boreholes are generated using GP. By equating the proposed equations to each other, water curtain system can be pre-designed. Relative error of the generated GP equations shows their ability and accuracy. Applying a standard regression coefficient method, sensitivity analysis of parameters related to water curtain performance and water inflow to the caverns is performed as well. The results help the design of the water curtain system for crude oil storage caverns worldwide.