The differential tomographic reconstruction method implements, for a modeled object, a method for estimating the amount of energy received in a THz detection system after having emitted by a source and after having passed through said object, method of estimation which comprises: - a step of simulating a THz beam emitted by the source composed of a plurality of rays (30 to 33) whose orientation and eccentricity with respect to the optical axis are chosen randomly according to a Gaussian model based on Monte Carlo type random methods corresponding to the opening of a detection system, - a step of simulation of propagation of each THz ray in the materials composing the object according to Snell-Descartes laws taking into account account of the losses by reflection and of the deviation of the ray at each interface and - a simulation step of attenuation of the intensity of each ray by considering a model of attenuation n exponential according to Beer-Lambert law. Le procédé de reconstruction tomographique différentielle, met en œuvre, pour un objet modélisé, un procédé d'estimation de la quantité d'énergie reçue dans un système de détection THz après avoir émise par une source et après avoir traversé ledit objet, procédé d'estimation qui comporte : - une étape de simulation d'un faisceau THz émis par la source composé d'une pluralité de rayons (30 à 33) dont l'orientation et l'excentrement par rapport à l'axe optique sont choisis aléatoirement selon un modèle Gaussien basé sur les méthodes aléatoires de type Monte Carlo correspondant à l'ouverture d'un système de détection, - une étape de simulation de propagation de chaque rayon THz dans les matériaux composant l'objet selon les lois de Snell-Descartes prenant en compte des pertes par réflexion et de la déviation du rayon à chaque interface et - une étape de simulation d'atténuation de l'intensité de chaque rayon en considérant un modèle d'atténuation exponentielle selon la loi de Beer-Lambert.
