Hello geo_enth,
Currently there's no instrument for fast polygon selection using Python API based on the polygon size.
So the workaround would be only custom function creation that estimates the size of the mesh faces (as the triangle area based on the vertex coordinates) in the model, calculates the min/max values and then selects and removes the faces based on the input threshold. However, this approach would likely work slow for high-poly models.