The main aim of the work was the development and research of a numerical algorithm designed to correct artifacts and distortions that occur during image reconstruction in optical-acoustic tomography tasks. The task was to develop an algorithm capable of compensating for the peculiarities of the reconstruction method. The proposed iterative scheme for enhancing optoacoustic images is based on Banach's fixed point theorem. The problem of eliminating distortions and artifacts in optoacoustic images has been modeled and tested for four iterative schemes. Two-dimensional and three-dimensional numerical models were used to investigate the efficiency of the proposed algorithm. The optical-acoustic experiment simulated a biological environment with an object to be reconstructed built into it. The linear (2D case) or planar (3D case) detecting acoustic arrays located on the surface of the studied samples were considered. The reconstruction quality was determined using both the quantitative and visual estimates of the obtained results. The efficiency of the iterative algorithm of enhancing the reconstruction quality was quantitatively estimated using the structure similarity index SSIM and relative reconstruction error E. It was shown that the developed algorithm makes it possible to considerably improve the image quality as compared with the traditional optical-acoustic tomography methods even with a small number of iterations. Algorithms of the k-Wave Matlab toolbox software package were used for comparing. The results obtained in the work can be important from the standpoint of the prospects of their further practical application in the tasks of biomedical visualization.