The combination of preoperative chemotherapy and surgical treatment has been shown to significantly enhance the prognosis of colorectal cancer with liver metastases (CRLM) patients. Nevertheless, as a result of variations in clinicopathological parameters, the prognosis of this particular group of patients differs considerably. This study aimed to develop and evaluate Cox proportional risk regression model and competing risk regression model using two patient cohorts. The goal was to provide a more precise and personalized prognostic evaluation system. We collected information on individuals who had a pathological diagnosis of colorectal cancer between 2000 and 2019 from the Surveillance, Epidemiology, and End Results (SEER) Database. We obtained data from patients who underwent pathological diagnosis of colorectal cancer and got comprehensive therapy at the hospital between January 1, 2010, and June 1, 2022. The SEER data collected after screening according to the inclusion and exclusion criteria were separated into two cohorts: a training cohort (training cohort) and an internal validation cohort (internal validation cohort), using a random 1:1 split. Subgroup Kaplan-Meier (K-M) survival analyses were conducted on each of the three groups. The data that received following screening from the hospital were designated as the external validation cohort. The subsequent variables were chosen for additional examination: age, gender, marital status, race, tumor site, pretreatment carcinoembryonic antigen level, tumor size, T stage, N stage, pathological grade, number of tumor deposits, perineural invasion, number of regional lymph nodes examined, and number of positive regional lymph nodes. The primary endpoint was median overall survival (mOS). In the training cohort, we conducted univariate Cox regression analysis and utilized a stepwise regression approach, employing the Akaike information criterion (AIC) to select variables and create Cox proportional risk regression models. We evaluated the accuracy of the model using calibration curve, receiver operating characteristic curve (ROC), and area under curve (AUC). The effectiveness of the models was assessed using decision curve analysis (DCA). To evaluate the non-cancer-related outcomes, we analyzed variables that had significant impacts using subgroup cumulative incidence function (CIF) and Gray's test. These analyses were used to create competing risk regression models. Nomograms of the two models were constructed separately and prognostic predictions were made for the same patients in SEER database. This study comprised a total of 735 individuals. The mOS of the training cohort, internal validation cohort, and QDU cohort was 55.00months (95%CI 46.97-63.03), 48.00months (95%CI 40.65-55.35), and 68.00months (95%CI 54.91-81.08), respectively. The multivariate Cox regression analysis revealed that age, N stage, presence of perineural infiltration, number of tumor deposits and number of positive regional lymph nodes were identified as independent prognostic risk variables (p < 0.05). In comparison to the conventional TNM staging model, the Cox proportional risk regression model exhibited a higher C-index. After controlling for competing risk events, age, N stage, presence of perineural infiltration, number of tumor deposits, number of regional lymph nodes examined, and number of positive regional lymph nodes were independent predictors of the risk of cancer-specific mortality (p < 0.05). We have developed a prognostic model to predict the survival of patients with synchronous CRLM who undergo preoperative chemotherapy and surgery. This model has been tested internally and externally, confirming its accuracy and reliability.