dk =Matrix([[op],[0],[-o]])
dk*2

>>> op,o = symbols('op o')
>>> dk =Matrix([[2*op],[o+op**2],[-op/o]])
>>> print([[cell.as_coeff_exponent(op) for cell in row] for row in dk.tolist()])
[[(2, 1)], [(o + op**2, 0)], [(-1/o, 1)]]