A four-state one-dimension discrete spin model with the nearest neighbor and partial next nearest neighbor lattice interactions is proposed for describing electrical double layer formed by ionic liquid in cylindrical pore that can accommodate only one row of particles. The model, in which cation is modeled as a dimer formed by two hard spheres (HSs): one is positively charged with diameter , while the other is neutral with diameter , and anion is modeled by charged HS with diameter , is designed to investigate effects of pore diameter , ion transfer free energy, ion size and electrical valence on behaviors of specific differential capacitance and specific energy storage and its saturation value . Relevant Kramers–Wannier transfer matrix is , which gives system partition function by numerically solving maximum real root of the relevant characteristic equation. Main findings are briefly described below. (i) The ion valence has little effect on the value, but does reduce the threshold voltage strength , beyond which the energy storage reaches . However, the ion valence (whether co- or counter-ion) raises both the and values as the relevant ion size decreases enough. (ii) Decreasing of the value weakly reduces the value, but raises the peak height and the value obviously; correlation between the value and value becomes stronger with the ions getting more pore hostile. (iii) With the counter-ion or co-ion charged site size decreasing, both the value and value rise. However, the effect of the co-ion charged site size is weakened as the co-ion gets more pore hostile, and the cation neutral site size influences on the value less than the cation charged site size does. (iv) Pore friendly co-ion or pore hostile counter-ion always helps in raising the value; the two have significant synergistic effect with the co-ion pore friendliness playing more important role than the counter-ion pore hostility.