At present, the inversion of controlled source audio-frequency magnetotelluric data is mainly based on the isotropic media. Because of the rock fracture or the shape of mineral, the resistivity is different in each direction in the complicated geoelectric body. To fit actual geological conditions better, this paper mainly probes into the inversion of controlled source audio-frequency magnetotelluric data in layered media with azimuthal anisotropy induced by ground orthogonal horizontal electric dipoles. We develop an Occam inversion algorithm to invert controlled source audio-frequency magnetotelluric data for azimuthal anisotropy in layered media. Designing the theoretical geoeletric models, the synthetic examples demonstrate the validity and stability of the Occam inversion algorithm. At the same time, the inversion results of different data suggest that in order to get better inversion results of azimuthally anisotropic layered media the apparent resistivities ρ_sxy,ρ_syx, and phases φ_xy,φ_yx, should be used.