Hydraulic stimulation technologies are currently widely applied across resource and power generation industries to increase productivity of oil/gas or hot water reservoirs. These technologies utilise a pressurised fluid, which is applied inside the well to initiate and drive fractures as well as to open a network of natural fractures. To prevent the opened fractures from complete closure during production stage, small particles (proppants) are normally injected with the pressurised fluid. These particles are subjected to confining stresses when the fluid pressure is removed, which lead to a partial closure of the opened fractures.The residual fracture openings are the main outcome of such hydraulic stimulations as these openings significantly affect the permeability of the reservoirs and, subsequently, well productivity. Past research was largely focused on the assessment of conditions and characteristics of fluid driven fractures as well as proppant placement techniques. Surprisingly, not much work was devoted to the assessment of the residual fracture profiles. In this work we develop a simplified non-linear model of residual closure of a crack filled with deformable particles under remote compressive stresses. It is demonstrated that the closure profile is significantly influenced by the distribution and compressibility of the particles, which are often ignored in the existing analytical or semi-analytical models.