selecionar features cuja soma da tabela de atributos ordenada é >=x

select_sorted_sum.py

# pega o layer
layer = qgis.utils.iface.activeLayer()

# acessa a tabela de atributos
features = layer.getFeatures()

# ordena a tabela ded atributos
# aqui um teste
for feat in features:
    print(sorted(feat["producao"]))
    
# valor que será o limiar de seleção
x = 2000

# seleciona todas as features cuja soma ordenada pelo campo producao" é <=2000

Segue

field_name = 'producao'
def field_to_order(f):
    """Essa funcao retorna o valor do campo a ser ordenado"""
    return f[field_name]

features = sorted(layer.getFeatures(), key=field_to_order)

for feat in features:
    print(feature[field_name]