Starting from the hamiltonian for the Heisenberg ferromagnet which comprise randomly distributed nonmagnetic ions as impurities in a Bravais lattice, we express the spin operators by means of the Dyson-Maleev transformation in terms of the Bose operators of the second quantization. Then by using methods of quantum statistical field theory, we derive the partition function and the free energy for the system. We adopt the Matsubara thermal perturbation method to a portion of the hamiltonian which describes the interaction between magnons and the stationary field of nonmagnetic ions. Upon averaging over all possible distributions of impurities, we express the free energy of the system as a function of the mean impurity concentration. Subsequently, we set up the double-time single particle Green function at temperature T in the momentum space in terms of magnon operators and derive the equation of motion for the Green function through the Heisenberg equation of motion and then solve the resulting equation. From this, we calculate the self-energy and then the spectral density function for the system. We apply the formalism to the case of the simple cubic lattice and compute the density of states, the spectral density function and the lifetime of the magnons as a function of energy for several values of the mean concentration of nonmagnetic ions in the ferromagnetic lattice. We calculate the magnon energy spectrum as a function of the average impurity concentration fraction c, which shows that for low lying states, the excitation energy increases continuously with c in the studied range 0.1 ≤ c ≤ 0.7. We also use the spectral density function to compute some thermodynamical quantities through the magnon occupation number. We have obtained closed form expressions for the configurationally averaged physical quantities of interest in a unified fashion as functions of the mean concentration of nonmagnetic impurities c to any order of c applicable below a critical percolation concentration c p . The quantities of interest comprise the thermodynamic potential (free energy), the spin-wave self-energy and the spectral density function from which other quantities can be derived.