In ultra-dense networks (UDN), multiple association can be regarded as a user-centric pattern in which a user can be served by multiple base stations (BSs). The data rate and quality of service can be improved. However, BSs in user-centric paradigm are required to serve more users due to this multiple association scheme. The improvement of system performance may be limited by the improving load of BSs. In this letter, we develope an analytical framework for the load distribution of BSs in heterogeneous user-centric UDN. Based on open loop power control (OLPC), a user-centric scheme is considered in which the clustered serving BSs can provide given signal to interference plus noise ratio (SINR) for any typical user. As for any BS in different tiers, by leveraging stochastic geometry, we derive the Probability Mass Function (PMF) of the number of the served users, the Cumulative Distribution Function (CDF) of total power consumption, and the CDF bounds of downlink sum data rate. The accuracy of the theoretical analysis is validated by numerical simulations, and the effect the system parameters on the load of BSs is also presented.