Seismic observations have revealed two large low shear velocity provinces (LLSVPs) in the lowermost mantle beneath Pacific and Africa. One hypothesis for the origin of LLSVPs is that they are caused by accumulation of subducted oceanic crust on the core-mantle boundary (CMB). Here, I perform high resolution geodynamical calculations to test this hypothesis. The result shows that it is difficult for a thin (~ 6 km) subducted oceanic crust to accumulate on the CMB, and the major part of it is viscously stirred into the surrounding mantle. Another hypothesis for the origin of LLSVPs is that they are caused by thermochemical piles of more-primitive material which is remnant of Earth's early differentiation. In such case, a significant part of the subducted oceanic crust would enter the more-primitive reservoir, while other parts are either directly entrained into mantle plumes forming on top of the more-primitive reservoir or stirred into the background mantle. As a result, mantle plumes entrain a variable combination of compositional components including more-primitive material, old oceanic crust which first enters the more-primitive reservoir and is later entrained into mantle plumes with the more-primitive material, young oceanic crust which is directly entrained into mantle plumes without contacting the more-primitive reservoir, and depleted background mantle material. The result reconciles geochemical observation of multiple compositional components and varying ages of oceanic crust in the source of ocean-island basalts. Seismic studies have detected ultra-low velocity zones (ULVZs) in some localized regions on the CMB. Here, I present 3D thermochemical calculations to show that the distribution of ULVZs provides important information about their origin. ULVZs with a distinct composition tend to be located at the edges of LLSVPs, while ULVZs solely caused by partial melting tend to be located inboard from the edges of LLSVPs. This indicates that ULVZs at the edges of LLSVPs are best explained by distinct compositional heterogeneity, while ULVZs located insider of LLSVPs are better explained by partial melting. The results provide additional constraints for the origin of ULVZs.