Safe mine operations and optimal economical decision making in the context of lignite resources require an adequate level of knowledge about the spatial distribution of critical attributes in terms of geometry and quality in the deposit. Therefore, ore body models are generated using different approaches in geostatistics, depending on the problem to be solved. In this article the analysis of geostatistical methods used for deposits modeling has been presented. Based on exploration data concerning caloric value Q, models of one exemplary lignite deposit has been made. Two models of deposit were prepared using two different methods: ordinary kriging (OK) and sequential Gaussian conditional simulation (SGSIM). Different models of the same deposit were analyzed and compared with source data using criterion of fidelity to statistical attributes like: mean value, variance, statistical distribution. Models, which have been created based on exploration data, were compared with in-situ data gained from survey activities in the exploitation process. As a result of comparison correlation factor and measures of deviations were computed: average relative error, absolute relative error. Models were compared with in situ data, considering statistical features and local variability as well. In conclusion, the study gives valuable information into the benefits of using certain geostatistical approaches for variable tasks and problems in the lignite deposits design process. For the assessment of average values of deposit parameters ordinary kriging provides appropriate effects. Geostatisical simulation (e.g. sequential Gaussian simulation - SGSIM) provides much more relevant information for tasks connected to probability (or risk) of defined threshold exceedences than ordinary kriging. Models made with simulation method are characterized by high fidelity of spatial distribution in comparison to source data.

Bełchatów lignite deposit is located in the central partof Poland in the tectonic Kleszczów graben. It is dividedinto several parts, which are mining fields: Kamieńsk area (eastern part of the deposit), Bełchatów area (central partof the deposit) and Szczerców area (western part of the deposit). The subject of this study was the Belchatow area.The main issue of the investigations was the dependence of local, regional and global, horizontal variability of selected lignite qualitative parameters (moisture, ash content, calorific value and sulfur content in the as receivedstate) is a function of viewing direction. There was applied the geostatistical analysis of the lignite variability parameters with use of semivariograms.

The researches which were conducted at different scales of observation: in the locale scale – in small field size8 ́8 m called experimental area (local analysis), in larger homogeneous separated parts of the Belchatow area(regional analysis) and in the whole Bełchatów area scale (general analysis). The results proved the visible anisotropy of variability mine lignite parametres. Anisotropy structure observed in regional and global scale isconnected with tectonic structure of the Bełchatów Graben. The detailed studies show the variated level of anisotropy observed in different areas of Bełchatów field.

However, no dependence of the relative level of ash and total sulfur content anisotropy on the environment of sedimentation of the main coal deposit in different parts of the Belhchatów field has been observed. Both parameters characterize with strong or medium anisotropy level in examined fields. Moreover, anisotropy is alsovisible in the local scale. Conducted researches confirmed the thesis that zonal anisotropy is prevalent kind of anisotropy in the regional scale. In the range of the whole deposit the total sulfur content showed zonal anisotropy,whereas the ash content revealed geometric anisotropy