In experiments the two-dimensional systems are realized mainly on solid substrates which introduce quenched disorder due to some inherent defects. The defects of substrates influence the melting scenario of the systems and have to be taken into account in the interpretation of the experimental results. We present the results of the molecular dynamics simulations of the two dimensional system with the core-softened potential in which a small fraction of the particles is pinned, inducing quenched disorder.The potentials of this type are widely used for the qualitative description of the systems with the water-like anomalies. In our previous publications it was shown that the system demonstrates an anomalous melting scenario: at low densities the system melts through two continuous transition in accordance with the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) theory with the intermediate hexatic phase, while at high densities the conventional first order melting transition takes place. We find that the well-known disorder-induced widening of the hexatic phase occurs at low densities, while at high density part of the phase diagram random pinning transforms the first-order melting into two transitions: the continuous KTHNY-like solid-hexatic transition and first-order hexatic-isotropic liquid transition. Despite almost forty years of investigations, the controversy about the microscopic nature of melting in two dimensions (2D) still lasts. In a crystal in contrast to an isotropic liquid two symmetries are broken: translational and rotational. These two symmetries are not independent, since a rotation of one part of an ideal crystal with respect to another part disrupts not only the orientational order but also the translational order. However, it is possible to imagine the state of matter with orientational order, but without the translational one. As it was shown by Mermin [1] in two dimensions the long-range translational order can not exist because of the thermal fluctuations and transforms to the quasi-long-range one. On the other hand, the real long range orientational order does exist in this case.These ideas were used in the widely accepted KTHNY theory [2][3][4][5] where it was proposed that, in contrast to the 3D case where melting is the first-order transition, the 2D melting can occur through two continuous transitions. In the course of the first transition the bound dislocation pairs dissociate at some temperature T m transforming the quasi-long range translational order into the short-range one, and long-range orientational order into the quasi-long range order. At this transition the decay of the translational correlation function will change from algebraic to exponential, and the orientational correlation function will obey the algebraic decrease. The new intermediate phase with the quasi-long range orientational order is called the hexatic phase. After consequent dissociation of the disclination pairs at some temperature T i the hexatic phase transforms into the isotropic liquid with the exponen...