Production Optimization from Multi-Fractured Horizontal Shale Gas Wells with a Dynamic Two-Phase Duel-Porosity Model