  We study glueball G production in gluonic penguin decay B -> X_s G, using
next-to-leading order b -> s g^* gluonic penguin interaction and an effective
coupling of a glueball to two gluons. The effective coupling allows us to study
the decay rate of a glueball to two pseudoscalars in the framework of chiral
perturbation theory. Identifying the f_0(1710) to be a scalar glueball, we then
determine the effective coupling strength with the branching ratio of f_0(1710)
-> K KBar. We find that the charm penguin to be important and obtain a sizable
branching ratio for Br(B -> X_s G) in the range of (0.7 - 1.7) \times 10^{-4}.
Rare hadronic B decay data from BaBar and Belle can provide important
information about glueballs.
