Study on impurity hole plasmas by global neoclassical simulation