Summary: A new method of modeling Gaussian homogeneous isotropic random fields on the basis of \(\text{Sub}(\Omega)\) theory is presented. The approximation of the model with given accuracy and reliability in the space \(L_p,p\geq1\) is considered.