The main goal of the work was to further improve the efficiency of a numerical algorithm designed to correct artifacts and distortions that arise during image restoration in optical-acoustic tomography problems. The proposed iterative scheme for improving biomedical optoacoustic images is based on Banach's fixed point theorem. The challenge was to determine the fastest and most efficient acceleration method applicable to the developed algorithm. In a numerical optical-acoustic experiment, a biological environment with a reconstructed object built into it was simulated. The simulations were run for two-dimensional and three-dimensional numerical models to study the effectiveness of the proposed algorithms. The quality of the reconstruction was determined by both quantitative and visual assessments of the results obtained. The convergence of different acceleration methods and the quality of the output images were analyzed in terms of PSNR (peak signal-to-noise ratio) and SSIM (structure similarity index).