In order to simulate the propagation characteristics of disturbances in solar chromosphere, transition layer and corona, we discuss numerically the propagation process of a disturbance produced by local heating in vicinity of transition layer, by adopting a 2-dimensional MHD simulation under a self-consistent non-isothermal, non-uniform initial state. The results indicate that the disturbance propagates to ambient plasma with a speed of local fast MHD wave. We find the numerical results can fit in with that computed via local fast MHD wave, its characteristics can explain wave events observed by SOHO/EIT.