A research approach of "dismantling and building" is proposed for high-precision ocean scene infrared simulation of sea surface with state level above 3. In this approach, the large-scale sea surface is described as a combination of rough sea surface, foam and breaking waves; then three kinds of featured multi-scale voxels, namely "rough sea surface", "rough sea surface covering foam", and "rough sea surface covering breaking waves", are abstracted from the complex sea surface. Finally, by means of sea surface gridding, featured voxels matching, rendering based on directional-spectral characteristics reconstruction and other methods, the large-scale sea surface radiation/scattering characteristics are reconstructed with the three kinds of featured multi-scale voxels and the calculation of large-scale sea surface "air-surface-body" coupling radiation/scattering characteristics is completed. This article conducted modeling research on multi-scale coupled radiation and scattering characteristics of three kinds of featured multi-scale voxels. A calculation model for the radiation/scattering directional-spectral characteristics of three kinds of featured multi-scale voxels was constructed and the factors influencing the radiation scattering direction spectral characteristics of the structures were analyzed. The calculation results indicate that with the increase of the sea surface wind speed, the thickness of foam and the concentration of bubbles in the typical structure on the sea surface gradually increase, which makes the scattering ability of the structure increase, thus increasing the bidirectional reflectance distribution function of the structure. With the increase of the detection wavelength, the absorption of seawater is significantly enhanced, which leads to a significant increase in the difference between the bidirectional reflectance distribution function of structures under different wind speeds. For different incident angles, the zenith angle corresponding to the maximum value of bidirectional reflectance distribution function of any structure gradually changes along with the incident zenith angle.