This study considers a nonlinear grey Bernoulli forecasting model with conformable fractional-order accumulation, abbreviated as CFNGBM$(1, 1, \lambda)$, to study the gross regional product in the Cheng-Yu area. The new model contains three nonlinear parameters, the power exponent $\gamma$, the conformable fractional-order $\alpha$ and the background value $\lambda$, which increase the adjustability and flexibility of the CFNGBM$(1, 1, \lambda)$ model. Nonlinear parameters are determined by the moth flame optimization algorithm, which minimizes the mean absolute prediction percentage error. The CFNGBM$(1, 1, \lambda)$ model is applied to the gross regional product of 16 cities in the Cheng-Yu area, which are Chongqing, Chengdu, Mianyang, Leshan, Zigong, Deyang, Meishan, Luzhou, Suining, Neijiang, Nanchong, Guang'an, Yibin, Ya'an, Dazhou and Ziyang. With data from 2013 to 2021, several grey models are established and results show that the new model has higher accuracy in most cases.