A search for anomalous Z-Z-gamma and Z-gamma-gamma couplings in the process e+e- -> Z gamma has been performed using data collected with the OPAL detector at LEP2 energies. The study is based on the analysis of the total rate and of the shape of distributions of kinematic variables for the qqbar-gamma and nu-nubar-gamma final states. Consistency with the Standard Model expectations has been observed, indicating no evidence for anomalous couplings, and 95% CL upper limits have been derived.