In this study, the thermoacoustic instabilities of partially premixed hydrogen flames in a lean direct injection (LDI) multi-cluster combustor are characterised using dynamical systems theory. The combustor is operated at a range of bulk velocities (30 - 90 m/s) and equivalence ratios(0.2-0.6), and time-resolved pressure oscillations and integrated OH* chemiluminescence measurements were taken. The thermoacoustic system reveals a variety of dynamical states in pressure such as period-1 Limit Cycle Oscillation (LCO) with a single characteristic frequency, period-2 LCO with two characteristic frequencies, intermittent, quasi-periodic and chaotic states as either bulk velocity or equivalence ratio is varied. At a bulk velocity of 30 m/s, as the equivalence ratio is gradually decreased from 0.6 to 0.2, the dynamical behaviour follows a sequence from an intermittent state to a period-1 LCO, then to a quasi-periodic state, and eventually reaches a chaotic state. As the equivalence ratio is decreased for a bulk velocity of 60 m/s, the pressure oscillations evolve from a period-2 LCO to quasi-periodic state before flame blows off. The emergence of period-2 and quasi-periodic states indicates the presence of strong non-linear interactions among the cavity acoustic modes. The emergence of period-2 and quasi-periodic states indicates strong non-linear interactions among the cavity acoustic modes. These modes and their spatial behaviour are investigated using a Helmholtz solver, showing that hydrogen flames can excite a wide range of cavity acoustic modes