The purpose of this work to reproduces the flux of lines emitted by the ionized gas, resulting in a good agreement between observed and predicted line fluxes, temperatures and electron density with the chemical abundances of ionic species with carbon and oxygen rich planetary. The model has been calculated using clouds of dust grains in order to calculate the electron temperature, density, and emission fluxes. The dust grain with polycyclic aromatic hydrocarbons (PAH) physics combined with thick shell geometry solves the problem of the heating and allows the fitting of fluxes lines, the electron temperatures and the ionization structure of the nebula. Among other lines, we determined line ratio of CIIλ4267A, [NII] λ5755A, 6584A, 6548A, [SII] λ4078A, 4070A, 6716A, and 6731A and compared the results with the observation data by changing the abundances of C, N, O and S, which are enhanced by a factor of 2. The analysis of the ionization structure of the gas, as well as this calculation of the flux of ions, is presented in this work. Therefore, our result revealed that flux lines like [NII] λ6584A, [OII] λ3727A+29A and [OIII] λ5007A greater than the observed fluxes. The nebulae are found to have low electron temperature than found by the previous author, but the electron density is high.