In　this　paper,　a　block-centered　finite　difference　method　is　derived　for　the　Darcy-Forchheimer　compressible　wormhole　propagation　in　porous　media　by　introducing　an　auxiliary　flux　variable　to　guarantee　full　mass　conservation.　Error　estimates　for　the　pressure,　velocity,　porosity,　concentration,　and　auxiliary　flux　in　different　discrete　norms　are　established　rigorously　and　carefully　on　nonuniform　grids.　Finally,　some　numerical　experiments　are　demonstrated　to　verify　the　theoretical　analysis　and　effectiveness　of　the　given　scheme.