In this article, the principles of the gray model and BP neural network model are analyzed, and the characteristics of land-use change and spatial and temporal distribution are studied in-depth, and at the same time, to explore the influence of land-use change on ESV, the relationship between the two is analyzed using gray correlation degree, and a mathematical model is constructed to maximize the benefits of the regional system, coupling economic and ecological benefits, combined with Geo SOS-FLUS model to achieve the optimization of land use. This article constructs a combined prediction model of a gray neural network. The gray differential equation parameters correspond to the weights and thresholds of the neural network, and the optimized parameters are determined by training the neural network to make it stable. Then the training results of the BP neural network are fitted with the results obtained from the gray GM (1.1) model. Finally, the prediction results of the three models, gray GM (1.1), BP God Meridian, and gray neural network model, are compared and analyzed. The global spatial autocorrelation and local spatial aggregation patterns of regional soil erosion and its erosion factors are analyzed using the Exploratory Spatial Data Analysis (ESDA) method in spatial measurement theory.