Convective water-tank is used to model the development of the convective boundary layer on thermal heterogeneous surface. Mosaic-like cover is put on the bottom of the water-tank to achieve thermal heterogeneity. In the convective boundary layer development, temperature profiles were measured on different surface region and with which heat flux profiles are calculated. The results show that normalized heat flux profiles behave the non-linear characteristics, that is, the convex and concave features. Convex and concave features oscillate over time. This is obviously different from the case of homogeneous surface. From heat flux contour map, the oscillation of the overall average heat flux profiles is smaller than that of the uncovered and the covered surface regions. And the oscillation on the uncovered surface region has a higher frequency and greater amplitude. The results show that some turbulence structures of convective boundary layer on thermal heterogeneous surface are quite different from that on the homogeneous surface. The law of the parameters for homogeneous surface can not be applied alike for heterogeneous surface. The experimental results can provide useful reference for numerical simulation and field observations.