In the Kerch-Tamagne areas, mud volcanoes and associated "pressed synclines" are evolving. The problem of genesis of mud volcanoes is debated; researchers of mud volcanoes do not have a common opinion about the mechanism of their formation. The necessity of studying the mud volcanoes structures is also connected with the fact that some mud volcanoes in the Tamagne mud volcanoes provinces represent a real danger for the population at moments of intensive eruptions because of the immediate proximity from settlements.

It is not excluded that processes associated with general decontamination of the Earth underlies the mechanism of formation of mud volcanoes. It is known that millions of cubic meters of methane per day are derived from the Black sea bottom.

Earlier, the proposed vibroseismic method of monitoring the magmatic structures with a controllable seismic source will make possible to receive a new knowledge about the structure of volcanoes, their origin and dynamics of the behavior of dilatancy structures of living volcanoes with vibroseismic monitoring of these structures.

The solution to the fundamental problem of studying the structures of mud volcanoes is related to the development of theoretical, methodical and experimental bases of sounding the dilatancy zones of volcanoes with application of powerful vibrators. One of the key problems of studying volcanoes of the given type is the construction an adequate mathematical models for given structures describing the propagation of elastic waves generated by a powerful vibrator.

It is known that mathematical modeling is an effective tool to study the processes of propagation of elastic waves in various models of complex subsurface geometries. Currently, a wide spectrum of numerical methods as applied to the full wave fields simulation in inhomogeneous elastic media is available. The difference method and the finite element method are the most flexible from all known methods of numerical modeling of elastic waves propagation in the case of complicated three-dimensional

inhomogeneous elastic media. However, their use demands high computer costs even with application of the cluster super computers.

The developed program is intended for the numerical modeling of elastic waves propagation in three-dimensional inhomogeneous models of elastic media, using a finite difference method.

The results of the first vibroseismic experiments of research into the mud volcano Shugo are presented. It was found out that the wave field of seismic waves has a complex structure in the volcano vicinity. Further direction of experimental research will be defined by the results of numerical modeling to be used when interpreting the data of vibroseismic sounding of mud volcanoes in the Tamagne province.

We propose the further investigation of volcanic pipes using vibroseismic methods with application of seismic vibrators and areal recording systems. It is assumed that the choice of the scheme of observation and interpretation of results of field experiments will be made on the basis of results of mathematical modeling.