Background Z score utility is emphasized in classifying coronary artery lesions in Kawasaki disease patients. The present study is the largest such multicenter Chinese pediatric study about coronary artery diameter reference values and Z score regression equation to date. It is useful in Chinese pediatric echocardiography. Methods A multicenter cohort was assembled, which consisted of 852 healthy children between 1 month and 17 years of age, ten children were excluded because their ultrasound images were not clear, or lost in following up. Diameters of the right coronary artery, left coronary artery, and left anterior descending coronary artery were assessed using echocardiography. Data were body surface area (BSA)-corrected using BSA calculated via either the Stevenson BSA formula or the Haycock BSA formula. Coronary artery diameter reference values and Z score regression equations were established for use in the Chinese pediatric population. Results No difference was observed between coronary artery diameter data corrected using BSAste or BSAhay. Of the five assessed regression models, the exponential model exhibited the best fit and was therefore selected as the basis for derivation of the SZ method. When comparing Z scores, those produced by the SZ method conformed to the standard normal distribution, while those produced by the D method did not. In addition, there was a statistically significant difference between Z scores produced by the SZ and D methods (P < 0.05). Conclusions Coronary artery diameter reference values for echocardiography were successfully established for use in the Chinese pediatric population, and a Z score regression equation more suitable for clinical use in this population was successfully developed.