In the article, universal methods of statistical modeling (Monte Carlo methods) of geophysical data using the Gaussian correlation function have been developed, which make it possible to solve the problems of generating adequate realizations of random fields on a grid in three-dimensional space of required regularity and detail. Since in geophysics, most of the results of object research are presented in digital form, the accuracy of which depends on various random influences, the problem of the condition of the maps arises in the case when the data cannot be obtained with the specified detail in some observation areas. It is proposed to apply statistical simulation of random fields methods, to solve the problems of conditional maps, supplement the required detail of research results with additional data, to achieve the required accuracy of observations, and other similar problems in geophysics. An algorithm for numerical modeling of realizations of homogeneous isotropic random fields in three-dimensional space with a Gaussian correlation function is formulated on the basis of the theorem on estimation of the mean-square approximation of such random fields by the partial sum of the "spectral decomposition" series. Using the example of data from aeromagnetic surveying in the area of the Ovruch depression, the proposed algorithm for statistical modeling of random fields is implemented in solving the problems of map fitness by supplementing the data with simulated adequate implementations to the required level of detail. When analyzing data by profiles, they are divided into deterministic (trend) and random components. The trend is proposed to approximate by cubic splines and the homogeneous isotropic random component is proposed to modeling on the basis "spectral decomposition" of random fields on 3-D space in the Ovruch depression. According to the algorithm, authors received random component implementations on the study area with twice detail for each profile. When checking their adequacy, authors made the conclusions that the relevant random components histogram has Gaussian distribution. The built variogram of these implementations has the best approximation by theoretical variogram which is connected to the Gaussian type correlation function. As a result of superimposing the simulated array of the random component on the spline approximation of the real data, a more detailed implementation was obtained for the data of geomagnetic observations in the selected area. A comparative analysis of the results of modeling realizations random fields with the Gaussian correlation function with other correlation functions is carried out. Therefore, the method of statistical modeling of realizations of random fields in three-dimensional space with the Gaussian correlation function makes it possible to supplement the results of measurements of the full magnetic field intensity vector with data with a given detail as much as possible.