Bessel beam array has been widely used in femtosecond laser processing, particle capture, optical microscopy, optical communication, and other fields. Especially in the field of industrial processing, the Bessel beam plays an important role in the process of pore structures with the ratio of height to depth due to its long focal depth characteristics. For the preparation of large-area periodic pore micro-nano structures, the parallel processing method of the Bessel beam array can significantly improve the machining efficiency. The machining quality of materials is closely related to the quality of the light field of the Bessel beam array, so it is significant to study the generation method of Bessel beam array with high quality. The traditional Bessel beam array generation methods include: multi-axicon phase serial superposition method, Dammann grating, and axicon phase superposition method, which can generate parallel and divergent Bessel beam array, respectively. However, the generated Bessel beam array has problems of poor uniformity and low diffraction efficiency. Therefore, two computational holography methods are proposed in this paper, which can generate high-quality parallel and divergent Bessel beam arrays respectively. Firstly, the computational hologram model of the proposed method is established, and the multi-axicon phase parallel splicing method is proposed, which effectively reduce the background noise of the optical field by improving the“aperture utilization ratio”of the window; The multi-lens and axicon phase superposition method is proposed, the multi-lens phase superposition method is used to generate multiple focus distributions on the observation plane, and then the beams of each focus are modulated into Bessel beams by superimposing axicon phase, thus forming Bessel beam array,the key of this method is the multi-lens phase superposition to generate multi-focus distributions with controllable position. Secondly, holograms of a 3×3 Bessel beam array are generated by the proposed method and the traditional method, and then simulated respectively, the transverse optical field distribution and diffraction pattern of the Bessel beam array in free space are obtained, the uniformity and diffraction efficiency of the Bessel beam array generated by the proposed method and the traditional method are compared and analyzed. The simulation results show that the uniformity and diffraction efficiency of the parallel Bessel beam array generated by the proposed method are 98.94% and 78.12%, respectively; the uniformity and diffraction efficiency of the diverging Bessel beam array generated by the proposed method are 97.95% and 79.23%, respectively. Finally, the images of 120 mm, 130 mm and 140 mm along the transmission direction of Bessel beam array are collected through experiments, which are highly consistent with the simulation results. Compared with traditional methods, the uniformity of parallel and divergent Bessel beam arrays produced by the proposed method is increased by 2.97% and 4.70%, respectively, and the diffraction efficiency is increased by 48.22% and 54.75%. The method proposed in this paper provides a technical approach to the generation of high quality Bessel beam arrays and has certain engineering application value.