Nanochannels are a powerful technique for detecting a wide range of biomolecules without labeling. The ion transport phenomena in nanochannel arrays differ from those in single nanochannels and are caused by interchannel communication. This study uses a fully coupled Poisson-Nernst-Planck (PNP) and Navier-Stokes model to investigate ion transport in nanochannel arrays. Instead of being set at a constant value, the surface charge density used in this study is established by the protonation and deprotonation of the silanol groups that are present on the walls of the silicon-based nanochannels. The surface charge density of the nanochannel walls varies with the number of nanochannels, the channel lateral distance, and the background solution properties, which consequently influence the ionic concentration distribution, flow velocity, and electric field strength. For example, in different numbers of nanochannel systems, the ion concentration in nanochannels is not much different, but it is different in reservoirs, especially near the openings of nanochannels. The number of nanochannels and the distance between nanochannels can also affect the formation of electro-convective vortex zones under certain conditions. These findings can aid in optimizing the nanochannel array design by regulating the number and distance of nanochannels and facilitating the construction of solid-state nanochannel arrays with any desired nanochannel dimensions.