Partial Differential Equations (PDEs) have been already widely used to simulate various complex phenomena in porous media. This paper is one of the first attempts to apply PDEs for simulating in real 3D structures. We apply this scheme to the specific case study of the microbial decomposition of organic matter in soil pore space. We got a 3D geometrical representation of the pore space relating to a network of volume primitives. A mesh of the pore space is then created by using the network. PDEs system is solved by free finite elements solver Freefem3d in the particular mesh. We validate our PDEs model to experimental data with 3D Computed Tomography (CT) images of soil samples. Regarding the current state of art on soil organic matter decay models, our approach allows taking into account precise 3D spatialization of the decomposition process by a pore space geometry description.