A Rijke tube with a distributed heat source is investigated. Driven by the widely existing thermoacoustic instability in lean premixed gas turbine combustors, this work aims to explore the physicochemical underpinning and assist in the elucidation and analysis of this problem. The heat release model consists of a row of distributed heat sources with individual heat release rates. The integrated heat release rate is then coupled with the acoustic perturbation for thermoacoustic analysis. A continuation approach is employed to conduct the bifurcation analysis and capture the nonlinear behaviour inherent in the system. Unlike the conventional approach by the Galerkin method, the acoustic equations are originally discretized using the Method of Lines (MOL) to build up a dynamic system. The model is first validated and shown to yield good predictions with available experimental data. Influences of multiple heat sources, time delay, and heat release distribution are then studied to reveal the extensive nonlinear characteristics involved in the case of a distributed heat source. It is found that distributed heat source plays an important role in determining the stability of a thermoacoustic system.