An important goal of the High-Luminosity run of the LHC is the measurement of the Higgs boson pair production cross-section, thereby placing constraints on the Higgs boson self-couplings. The latest European Strategy for Particle Physics update indicates that, under fairly conservative assumptions, the HL-LHC is expected to discover this process, motivating further detailed study. Theoretically, this process is known to N3LO in the Heavy Top-quark Limit (HTL), with NLO QCD and EW corrections known including quark masses. One of the dominant remaining theoretical uncertainties comes from the dependence of the prediction on the top-quark mass scheme and scale. In this talk, we discuss a recent attempt to use the Method of Regions and the tools of Soft-Collinear Effective Theory (SCET) to understand the structure of the quark mass corrections in the high-energy limit, defined by s,|t|,|u| >> mT^2 >> mH^2. We argue that the mass corrections follow a predicable factorised pattern at leading power and we present initial work towards a next-to-leading power understanding of the high-energy behaviour.