This study investigates terahertz (THz) wave scattering from a simulated lunar regolith surface, with a focus on the Brewster feature, backscattering, and bistatic scattering within the 325 to 500 GHz range. We employed a generalized power-law spectrum to characterize surface roughness and fabricated Gaussian correlated surfaces from Durable Resin V2 using 3D printing technology. The complex dielectric permittivity of these materials was determined through THz time-domain spectroscopy (THz-TDS). Our experimental setup comprised a vector network analyzer (VNA) equipped with dual waveguide frequency extenders for the WR-2.2 band, transmitter and receiver modules, polarizing components, and a scattering chamber. We systematically analyzed the effects of root-mean-square (RMS) height, correlation length, dielectric constant, frequency, polarization, and observation angle on THz scattering. The findings highlight the significant impact of surface roughness on the Brewster angle shift, backscattering, and bistatic scattering. These insights are crucial for refining theoretical models and developing algorithms to retrieve physical parameters for lunar and other celestial explorations.