To address the issue of high computational load and long processing time in the imaging calculation process of lightning radiation sources using time-reversal algorithms in the frequency domain, a series of optimized acceleration methods for time-reversal algorithms are proposed. Firstly, the introduction of Graphics Processing Unit (GPU) parallel computing is further expanded to dual-GPU cluster computing, while adopting multi-threading concurrent programming methods of Central Processing Unit (CPU) and GPU to cover the time consumption caused by mutual waiting; secondly, the maximum energy search algorithm and phase difference filtering algorithm and other programs are written to the GPU to achieve parallel search and calculation of phase difference values; finally, the data acquisition and transmission module, GPU processing module, and positioning imaging module are integrated to achieve the integrated function from data collection to positioning imaging of lightning data. For a 500 ms segment of Very High Frequency (VHF) lightning data, the actual measurement only takes 18 minutes, which is about 703 times faster than the previous CPU single-line processing process, and the number of effective positioning points is only reduced by 9.3%.