gdal.Translate(self.output, raster_dataset, format='GPKG', creationOptions=['APPEND_SUBDATASET=YES'])
It works.