The diverse structural, electronic, and magnetic properties of silicon (Si)-substituted armchair and zigzag graphene nanoribbons (AGNRs and ZGNRs) were investigated using spin-polarized density functional theory (DFT) calculations. Pristine AGNRs belong to a nonmagnetic semiconductor with a direct bandgap of 1.63/1.92 eV determined by PBE/HSE06 functionals. Under various Si substitutions, nonmagnetic bandgaps were tuned at 1.49/1.87, 1.06/1.84, 0.81/1.45, 1.04/1.71, 0.89/1.05, and 2.38/3.0 eV (PBE/HSE06) in the single Si edge-, single Si non-edge-, double Si ortho-, double Si meta-, double Si para-, and 100% Si-substituted AGNR configurations, respectively. Meanwhile, pristine ZGNRs displayed antiferromagnetic semiconducting behavior with a spin degenerate bandgap of 0.52/0.81 eV (PBE/HSE06) and becomes a ferromagnetic semimetal in the single Si configurations or an unusual ferromagnetic semiconductor in the 100% Si configuration. Under the developed first-principles theoretical framework, the formation of quasi π (C-2pz and Si-3pz) and quasi σ (C-2s, -2pxy and Si-3s and -3pxy) bands was identified in the Si-substituted configurations. These quasi π and quasi σ bands showed weak separation, resulting in weak quasi sp2 hybridization in Si-C bonds, in which the identified hybridization mechanism was a strong evidence for the formation of stable planar 1D structures in the Si-substituted configurations. Our complete revelation of the essential properties of Si-substituted GNRs can provide a complete understanding of their chemically doped 1D materials for various practical applications.