Nowadays, the chemical etching technique had become one of the most important processes in the advanced semiconductor manufacturing and microelectromechanical systems. In this study, a theoretical model based on the phase-field method is developed to investigate the formation and evolution of surface morphologies during chemical etching processes. In the formation mechanisms of etched surface morphologies, the anisotropic diffusion of atoms and the anisotropic etching kinetics are both taken into consideration. To enhance the applications of this theoretical model, various etching parameters such as the crystalline structure of substrate and the anisotropic etching rate, corresponding to the different etching conditions are also examined in the numerical calculations. Accompanying with the quantitative analyses of the simulations, the influence of the etching parameters on the formation of surface morphologies during the etching process could be distinctly demonstrated. With the systematic investigations, this study is expected to improve the technique of the chemical etching and the development of the new manufacturing parameters for the applications in the realistic fabrication.