Compaction of heterogeneous earth foundation having creeping property has been studied in the given work. The creeping property of compactible water-saturated heterogeneous soil is described by Maslov-Arutyunyan theory of tensile-creeping body in the interpretation of Florin V.A.
For theoretical research the process of soil compaction, its heterogeneity, are considered in the work by deformation module of the compactible mass. Mathematically it is expressed by power or exponential function which is changed by the depth of compactible soil body. Besides it is accepted here that for strongly compressible water-saturated heavy soil in the initial instant of time a part of load, suddenly applied load q to the soil, equal in value of structural compressive strength рstr, is immediately taken by the soil skeleton, i.e. 
Mathematical statement of this problem comes to the determination in this soil body of distribution of pressure in the porous fluid р(х, z, t), tension in the soil skeleton σ(z, t), and vertical movement of points of the upper surface S(t) (setting). Pressure in the porous fluid is identified by solution of integral-differential equation at the initial and boundary conditions.
An outer crust of earth is usually characterized by a large-scale heterogeneity of earth and rocks composing it. It is due to the sufficiently complex geological-tectonic structure of rocks, where one or another construction project is built. Dropping of heterogeneity of geological structure of the outer crust of earth may lead in the future to damage of engineering constructions, due to the settlement, development in the subfoundation. As is known, nonrigid properties of heterogeneous soils, actually, are changed together with point positions, and assumption about their homogeneity represents idealization of real conditions. Certainly, providing the heterogeneity, the mathematical problem is incomparably complex and, therefore, in such cases we often apply for different simplification kinds of model, acceptable from one or another point of view. At that in one group of problems parameters characterizing properties of the material are sectionally constant. This means that the investigated body consists of several homogeneous bodies. In the other group of problems they represent continuous functions of point position. Moreover, the soil, which module of deformation continuously grows with a depth is called as a continuously heterogeneous. Such model of the soil was presented in the works of E.K. Klein [2] and Popov [3] for solution of some contact problems of the theory of elasticity. G.К. Klein at calculation of constructions laying on the solid base, for the module of deformation of the soil has following expression
(1)
where En is the module of deformation of the soil at a depth z = 1; the figure n in most cases is in the limit of 0 < n < 2. In the works of Popov at solution of similar problems the module of deformation of the soil is in the following form
(2)
where E0, α are experimental data.
Furthermore, intrinsically the soils also relate to the rheologic bodies and deformations connected with creeping phenomenon, and find their reflection in compaction of the earth masses with rates which do not exceed several centimeters a year. Moreover the process of compaction of a construction foundation can be in a slow permanent motion without any features of its final stabilization. Notwithstanding such insignificant intensity of deformation, in some instances they can appear to be inadmissible for constructions located on the deformed foundation.
Mechanics of such soils under load and in time are usually described by rheologic equations of the soils state which connect between themselves its deformation, tension and their derivatives in time, and they had a form of linear differential equation. Later, integral relations were used for solution of the earth medium compaction problems. One of the first who used them for description of soil skeleton state was V.A. Florin [4]. At this, relationship between soil porosity factor ε(t) and sum of the main tensions θ(t) in a skeleton of homogeneous soil can be presented as follows:
(3)
where ε(τ1) and ε(t) are relatively coefficients of porosity for initial and final time point; ξ is a coefficient of the lateral earth pressure; τ is a time of load application; t is a time for which creeping deformation is determined; n is a dimension of compaction which possesses values properly 1, 2, 3; δ(t, τ) is a function throwing back changes in the soil porosity to the time point t from the unit load applied in the time point τ. For heterogeneous soil we take it in the form:
(4)
Here а0 is a coefficient of instantaneous compaction; а1 and γ1 are parameters of creeping determined from tests. М is a point of compactible body; C(τ, t) is a measure of creeping for homogeneous soil subjected to the compaction.
Therefore, let compactible medium contain not only fluid filling soil pores, but also some entrapped air or other gas; a soil skeleton is heterogeneous on the structure, i.e. some features reflecting physical-mathematical properties are functions from space coordinates. Particularly, the soil deformation module can be changed on the coordinate determining depth of the considered point position; in the large the soil can possess by ageing property; flow of a fluid filling soil pores is subject to the generalized filtration law of Darcy-Gersevanov; a soil skeleton is tensile-creeping porous medium is subject to the hereditary creep theory of G.N. Maslov – N.Kh. Arutyunyan [4]; in the initial instant of time a part of load, suddenly applied load q to the soil, equal in value of structural compressive strength рstr, is immediately taken by the soil skeleton, i.e. 
At this, relationship between soil porosity factor and sum of the main tensions can be in the form (3); soil; according to impermeability soil is anisotropic. Therefore compaction of the earth masses taking into account (3), (4) relating to the pore pressure amounts to the study of following integral-differential equation with variable coefficients:
(5)
Initial condition of equation (5) has the following form:
(6)
where φ(t) is a function of ageing; θ*, p* are relatively sum of the main tensions and pore pressure for stabilized condition of compactible body; р is a pressure in porous fluid; 
is a factor considering volume compression.
In the given work the equation (5) at (6) for cases (1) and (2) are solved for one dimensional compaction of heterogeneous soil having creeping property. These solutions reflect distribution of pressure in the porous fluid, tension in the soil skeleton, and vertical movement of points of the upper surface of compactible foundation soils. They have been received in the form of Bessel’s functions combination, equations in the form (5) for heterogeneous boundary conditions have been solved in the works [5, 6].
The work is submitted to the International Scientific Conference «Modern high technologies», Israel (Tel Aviv), April 25 – May 2, 2014. came to the editorial office оn 27.01.2014.