As one of the follow-up studies of Herschel Planetary Nebula Survey (HerPlaNS; Ueta et al. 2014), we focus on a bipolar planetary nebula (PN) NGC6781 to characterize the dusty nebula and the central star based on our own Herschel data and the rich archival spectroscopic/photometric image data in the wavelengths from UV to far-IR. With CLOUDY, we constructed a comprehensive photoionization model of NGC6781 ever made including data from UV to radio. We succeeded to reproduce the observed spectral energy distribution (SED) and the atomic gas, H2, CO, and OH molecular line fluxes. We found that about 40% of the total dust mass would be from warm-cold dust components.