Here is a syntax which works:

from scipy import optimize
from math import log

bounds = ((0, 1), (0, 1), (0, 1), (0, 1))

def some_function(x):
    return log(1+x[0])/(1+x[1]) + log(1+x[2])/(1+x[3])

results = dict()
print(optimize.shgo(some_function, bounds))