Геоакустическая эмиссия — это процесс генерации упругих волн горными породами в результате динамической перестройки их структуры. Результаты наблюдений показывают, что на динамику геоакустической эмиссии влияют механические процессы, протекающие в очаге готовящегося землетрясения. Ранее с целью обоснования связи между вариациями геоакустической эмиссии и процессом подготовки землетрясений было проведено моделирование зон геоакустической эмиссии — областей поверхности земной коры с деформациями порядка 10^(−8)–10^(−5). Результаты проведенного ранее моделирования показывают, что уровень расчетных деформаций в пунктах наблюдений превышает приливные, но на порядок отличается от зарегистрированных. Это может быть связано с тем, что земная кора рассматривалась в виде однородной среды. В действительности земная кора состоит из слоев горных пород, часть из которых находятся в закритическом состоянии и проявляют пластические и квазипластические свойства. Настоящая статья посвящена моделированию влияния неоднородностей земной коры на пространственное распределение зон геоакустической эмиссии. Неоднородности описываются системой простых сил, распределенной по поверхности сферического включения. Интенсивность действия сил полагается постоянной. Получены решения краевой задачи линейной теории упругости в виде свертки функций Грина для однородного изотропного упругого полупространства. Проведены вычислительные эксперименты, и построены линии уровней компонентов поля вектора смещений поверхности земной коры. Показано, что сферические включения оказывают влияние на поле вектора смещений поверхности земной коры. Характер влияния зависит от количества неоднородных включений и их расположения относительно очага готовящегося землетрясения. Geoacoustic emission is the process of elastic wave generation by rocks as the result of dynamic reconstruction of their structure. Observation results show that mechanic processes, occurring in the source of a preparing earthquake, affect the geoacoustic emission dynamics. Modeling of geoacoustic emission zones, the regions of the earth crust surface with deformations of the order 10^(−8)–10^(−5), has been earlier carried out to prove the relation between geoacoustic emission variations and the process of earthquake preparation. Results of the modeling, which was performed earlier, show that the level of calculated deformations at observation sites exceeds the tidal ones but differs by one order from the recorded deformations. This may be associated with the fact that the earth crust was considered as a homogeneous environment. In reality, the earth crust consists of rock layers, some part of which has supercritical state and manifests plastic and quasi-plastic properties. The present paper is devoted to the modeling of the earth crust inhomogeneities impact on spatial distribution of geoacoustic emission zones. Inhomogeneities are described by simple force system distributed over spherical inclusion surface. Intensity of the force action was assumed to be constant. Solutions for the boundary problem of elasticity linear theory were obtained in the form of Green’s functions convolution for homogeneous isotropic elastic half-space. Computational experiments were carried out, and lines of the field component levels of the displacement vectors of the earth crust surface were constructed. It was shown that spherical inclusions affect displacement vector field of the earth crust surface. The impact character depends on the number of inhomogeneous inclusions and their locations relative to the source of a preparing earthquake.