Gaussian random field is usually used to model Gaussian spatial data. In practice, we may encounter non-Gaussian data that are skewed. One solution to model skew spatial data is to use a skew random field. Recently, many skew random fields have been proposed to model this type of data, some of which have problems such as complexity, non-identifiability, and non-stationarity. In this article, a flexible class of closed skew-normal distribution is introduced to construct valid stationary random fields, and some important properties of this class such as identifiability and closedness under marginalization and conditioning are examined. The reasons for developing valid spatial models based on these skew random fields are also explained. Additionally, the identifiability of the spatial correlation model based on empirical variogram is investigated in a simulation study with the stationary skew random field as a competing model. Furthermore, spatial predictions using a likelihood approach are presented on these skew random fields and a simulation study is performed to evaluate the likelihood estimation of their parameters.