Seed size traits, including seed length (SL), seed width (SW), and seed thickness (ST), are crucial appearance parameters that determine soybean seed weight, yield, and ultimate utilization. However, there is still a large gap in the understanding of the genetic mechanism of these traits. Here, 281 soybeans were utilized to analyze the genetic architecture of seed size traits in different years through multiple (single-locus and multi-locus) genome-wide association study (GWAS) models, and candidate genes were predicted by integrating information on gene function and transcriptome sequencing data. As a result, two, seven, and three stable quantitative trait nucleotides (QTNs) controlling SL, SW, and ST were detected in multiple environments using the single-locus GWAS model, and concurrently detected by the results of the multi-locus GWAS models. These stable QTNs are located on 10 linkage disequilibrium blocks, with single genome regions ranging in size from 20 to 440 kb, and can serve as the major loci controlling soybean seed size. Furthermore, by combining gene functional annotation and transcriptome sequencing data of seeds at different developmental stages from two extreme soybean accessions, nine candidate genes, including Glyma.05G038000, Glyma.05G244100, Glyma.05G246900, Glyma.07G070200, Glyma.11G010000, Glyma.11G012400, Glyma.17G165500, Glyma.17G166500, and Glyma.20G012600 within the major loci that may regulate soybean seed size, were mined. Overall, these findings offer valuable insights for molecular improvement breeding as well as gene functional studies to unravel the mechanism of soybean seed size.