Based on the time domain Mueller matrix solution, a method for Very High Frequency(VHF, 30 MHz-300 MHz) radar pulse penetration on sea ice is developed. The polarimetric radar pulse echoes contain echoes from top and bottom interfaces, volumetric scattering of air-bubbles and interface-volume interactions. By changing model parameters such as the layer thickness, salinity, dielectric constant and interfacial roughness, the radar pulse echoes on sea ice are numerically simulated to verify the feasibility of the proposed method. The results demonstrate that the model can be used to explore sea ice, and the information such as depth and other structure properties can be inversed and estimated.