Holographic Data Storage (HDS) is one of the next generation storage technologies that can actualize high data capacity and high data transfer rate. Since information is recorded 3-dimensionally in a thick medium, data capacity of the HDS is not constrained by diffraction limit. However, behavior of wavefront in an inhomogeneous thick medium is highly complex, and it is hard to handle propagation of wavefront in the medium analytically. Therefore, we establish a numerical technique for analysis of volume holograms. The proposed technique is based on the scalar diffraction theory, which is described as the volume integral equation. By applying Born approximation and angular spectrum method to the volume integral equation, the technique can be applicable for various problems. We analyze characteristics of the volume hologram with spherical reference wave, and confirm effectiveness of the proposed technique. Compared to conventional techniques such as coupled wave analysis, beam propagation method, and finite-difference time domain method, the proposed technique has application potentiality for various problems, and it is easy to implement. In this study, we show effectiveness of the proposed technique by applying to analysis of the volume hologram with spherical reference wave. It can be expected that the proposed technique may become a tool for design of HDS systems.