The modification of Gaussian elimination is described here as numerical method to determine the physically adequate solution of the linear equations system that is affected by complexity of the data error (measurements, parameterization et. al.). Other modifications and algorithms are mainly established to efficiently overcome calculating issues. In order to avoid influence of errors in the data the proposed technique selects a set of invertible sub-matrixes. Developing the Tikhonov regularization theory the regularized subsequence is determined for each unknown under consideration. Application of the method to local seismicity data derived high resolution areas and revealed that the hydrothermal field of high concentration CH4 and CO2 is surrounded by rock having high seismic velocity. Based on this result we conclude that nature protects the gas flow due to special structure, which is likely presented by dense basalt. Decrypt of satellite images shows that the zone of lowest elevation does not intersect this structure. Particularities of background seismicity are examined for the gas area being studied. Factors learned from nature experience permit us to formulate conditions to select a place for underground gas storage.