Inter-satellite optical communication systems are regarded as the foundation of future space-ground integrated network. However, the optical antennas exposed in outer space will be affected by complicated factors and result in the occurrence of diversified wavefront distortions. These concerneddistortions, inevitable and of irregular shapes, may cause the degradation of communication performance and even lead to the instbility of communication system. In this paper, we focus on modeling of diversified wavefront distortion which is impossible to be described by analytic functions. Considering such characteristics, we propose to model such distortions via wavelet method, in which optimal wavelet basisfunctions are selected based on both computation complexityandmodeling precision. Simulations illustrate that wavelet method has advantages over existing results in the above mentioned two aspects.