We characterized the dusty circumstellar nebula and central star of the C-rich bipolar planetary nebula (PN) NGC 6781 using our own Herschel data augmented with the archival data from UV to radio and constructed one of the most comprehensive photoionization PN models ever produced consisting of the ionized, atomic and molecular gas components as well as the dust component. We reproduced the observed spectral energy distribution (SED), constrained by 136 observational data points. The total nebula mass was estimated to be 0.41 M⊙, with a significant fraction (about 70 %) of it existing in the photo-dissociation region (PDR) surrounding the ionized nebula. This finding demonstrates the critical importance of the PDR in PNe, which are typically recognized as the hallmark of ionized/H+ region. It is therefore essential to characterize the PDR of the circumstellar nebula to understand material recycling in the Milky Way and other galaxies.