This paper introduces a mathematical model designed to optimise the operation of natural gas distribution networks, considering the injection of hydrogen in multiple nodes. This proposal is devised to optimise the quantity of hydrogen injected to maintain pressure, gas flows, and gas quality indexes (Wobbe Index (WI) and the Higher Heating Value (HHV)) within admissible limits. The model has been applied to a gas network case study with three distinct scenarios and implemented using Python. The findings from the case study show the maximum permissible volume of hydrogen in the network, quantify the total savings in natural gas, and estimate the reduction in carbon dioxide emissions.