This paper focuses on the study of the limit cycle bifurcation of a general cubic Z2-equivariant Hamiltonian differential system with two nilpotent singular points. Initially, the necessary and sufficient conditions for the occurrence of two nilpotent singularities are provided in the unperturbed system, along with a complete description of all possible phase portraits on the plane. Subsequently, using a combination of the Melnikov function method and various analytical techniques, we derive the approximate expansion of the first order Melnikov function at the interval endpoint, as well as formula for its coefficients. These results can be utilized to analyze double homoclinic loop bifurcation of the perturbed system, finding the lower bound of limit cycles.