In this work, we are interested in the mathematical model of reaction diffusion systems. The originality of our study is to work with concentrations appearing in reactors together with measure initial data. To validate this model, we prove the existence of global weak solutions. The "j" technique introduced by Pierre and Martin [18] is suitable for this type of solutions. However, its adaptation has some new technical difficulties that we have to overcome.