An intelligent demand response (DR) program is developed for multi-energy industrial micro-grid consisting of manufacturing facilities, photovoltaic (PV) panels, and battery energy storage system (BESS). The proposed DR program tackles the practical challenges of components in the micro-grid including industrial process represented by a discrete manufacturing production model, uncertainty of PV generation, and operational cost of the BESS. The proposed DR program optimizes day-ahead production scheduling for manufacturing facilities and operation regime for the BESS in response to time of use electricity price and probabilistic forecasting of PV power. To capture the uncertainty of PV power, a data-driven PV power probabilistic forecasting model is developed and a copula-based approach is deployed for the sampling of temporally correlated scenarios of PV power over the scheduling horizon from probabilistic forecasts. The multi-energy management optimization problem is formulated as a scenario-based stochastic nonconvex mixed integer nonlinear programming (MINLP). A hybrid optimization method integrating the evolutionary algorithm and the branch-and-bound algorithm for mixed integer liner programming is proposed to solve the nonconvex MINLP. Simulation studies illustrate that the proposed DR program efficiently reduces the operational cost for manufacturing production and releases the stress of the main grid by making full use of flexibility of all the components in the micro-grid.