In this paper, an approach for neutral Higgs bosons search is described based on 2HDM type-I at electron-positron linear colliders operating at √ s = 1 TeV. The beam is assumed to be unpolarized and fast detector simulation is included. The signal process produces a fully hadronic final state through e + e − → AH → b bb b where both CP-even and CP-odd Higgs bosons (H and A) are assumed to decay to a pair of b-jets. Several benchmark scenarios are introduced as the baseline for the analysis taking m H/A in the range 150-300 GeV. In order to avoid Higgs boson conversion A → ZH, Higgs boson masses are chosen with mA − mH < mZ . It is shown that with a proper kinematic correction applied on final state b-jet four momenta, true combinations of b-jets can be found for simultaneous reconstruction of both Higgs bosons through b b invariant mass calculation. Results show that at integrate luminosity of 500 f b −1 , observable signals can be achieved with statistical significance exceeding 5σ.