This article takes the measured Total Electron Content (TEC) from the GPS points of the China Regional Crust Observation Network as the starting point to establish a regional ionospheric empirical model. The model’s performance is enhanced by considering solar flux and geomagnetic activity data. The refinement function model of the ionospheric TEC diurnal variation component, seasonal variation component, and geomagnetic component is studied. Using the nonlinear least squares method to fit undetermined coefficients, MEFM-ITCR (Multi-parameter Empirical Fusion Model–Ionospheric TEC China Regional Model) is proposed to forecast the regional ionosphere TEC in China. The results show that the standard deviation of MEFM-ITCR residuals is 3.74TECU, and MEFM-ITCR fits the modeling dataset well. Analyses of geographic location variation, seasonal variation, and geomagnetic disturbance were carried out for MEFM-ITCR performance. The results indicate that in the Chinese region, MEFM-ITCR outperforms IRI2020 and NeQuick2 models in terms of forecast accuracy, linear correlation, and model precision for TEC measured using GPS points under different latitudes and longitudes, different seasons, and different geomagnetic disturbances. The empirical TEC model built for the Chinese region in this paper provides a new ionospheric delay correction method for GNSS single frequency users and is of great significance for establishing other new and improving existing ionospheric empirical models.