To explore the clinical features, multimodal ultrasound features and multimodal ultrasound imaging features in predicting lymph node metastasis in the central cervical region of papillary thyroid carcinoma. A total of 129 patients with papillary thyroid carcinoma (PTC) confirmed by pathology were selected from our hospital from September 2020 to December 2022. According to the pathological results of cervical central lymph nodes, these patients were divided into metastatic group and non-metastatic group. Patients were randomly sampled and divided into training group (n = 90) and verification group (n = 39) according to the ratio of 7:3. The independent risk factors for central lymph node metastasis (CLNM) were determined by least absolute shrinkage and selection operator and multivariate logistic regression. Based on independent risk factors to build a prediction model, select the best diagnostic effectiveness of the prediction model sketch line chart, and finally, the line chart calibration and clinical benefits were evaluated. A total of 8, 11 and 17 features were selected from conventional ultrasound images, shear wave elastography (SWE) images and contrast-enhanced ultrasound (CEUS) images to construct the Radscore of conventional ultrasound, SWE and CEUS, respectively. After univariate and multivariate logistic regression analysis, male, multifocal, encapsulation, iso-high enhancement and multimodal ultrasound imaging score were independent risk factors for cervical CLNM in PTC patients (p < 0.05). Based on independent risk factors, a clinical combined with multimodal ultrasound feature model was constructed, and multimodal ultrasound Radscore were added to the clinical combined with multimodal ultrasound feature model to form a joint prediction model. In the training group, the diagnostic efficacy of combined model (AUC = 0.934) was better than that of clinical combined with multimodal ultrasound feature model (AUC = 0.841) and multimodal ultrasound radiomics model (AUC = 0.829). In training group and validation group, calibration curves show that the joint model has good predictive ability for cervical CLNM of PTC patients; The decision curve shows that most of the net benefits of the nematic chart are higher than those of clinical + multimodal ultrasound feature model and multimodal ultrasound radiomics model within a reasonable risk threshold range. Male, multifocal, capsular invasion and iso-high enhancement are independent risk factors of CLNM in PTC patients, and the clinical plus multimodal ultrasound model based on these four factors has good diagnostic efficiency. The joint prediction model after adding multimodal ultrasound Radscore to clinical and multimodal ultrasound features has the best diagnostic efficiency, high sensitivity and specificity, which is expected to provide objective basis for accurately formulating individualized treatment plans and evaluating prognosis.