In this work, a free boundary problem is presented for the attachment process in the initial phase of multispecies biofilm formation. The free boundary is represented by the biofilm thickness and it is assumed to be initially zero. The growth of attached species is governed by nonlinear hyperbolic PDEs. The free boundary evolution is governed by a first-order differential equation depending on the attachment, detachment, biomass velocity and substrates. The quasi-static diffusion of substrates is modelled by a system of semi-linear elliptic PDEs. The qualitative analysis of solutions leads to prove existence, uniqueness and some properties of solutions. We highlight that the free boundary velocity is greater than the characteristic velocity during the first instants of biofilm formation and the free boundary is a space-like line. It is proved that the attachment function depends linearly on the concentrations of all the attaching species. The first phase of biofilm growth is shown to be completely determined by environmental conditions and characterized by a specific mathematical inequality. The opposite inequality describes the further phase where the bulk liquid stops to directly affect the biofilm life. The mentioned inequalities could be assumed as rigorous definitions of non-mature and mature biofilms, respectively.