In this project, we conducted a comprehensive study on nuclear graphite properties with terahertz (THz) imaging. Graphite samples from Idaho National Laboratory were carefully imaged by continuous wave (CW) THz. The CW THz imaging of graphite shows that the samples from different billet with different fabricating conditions have different pore size and structure. Based on this result, we then used a phase sensitive THz system to study the graphite properties. In this exploration, various graphite were studied. By imaging nuclear graphite samples in reflection mode at nine different incident polarization angles using THz time-domain spectroscopy, we find that different domain distributions and levels of porosity will introduce polarization dependence in THz reflectivity. Sample with higher density is less porous and has a smaller average domain distribution. As a consequence, it is less polarization-dependent and the polarization-dependent frequency is higher. The results also show that samples oxidized at higher temperatures tend to be more polarization dependent. The graphite from the external billet is more polarization dependent compared to that from the center billet. In addition, we performed laser-based ultrasonic measurements on these graphite samples. The denser, unoxidized samples allow surface acoustic waves to propagate more rapidly than in the samples that had already undergone oxidation. Therefore, for the oxidized samples, the denser samples show less polarization-dependence, higher polarization-dependent frequency, and allow the surface acoustic waves propagate faster.