We analyze Bouguer anomalies in/around the focal region of the 2016 Kumamoto earthquake to examine features, such as continuity, segmentation and faulting type, of the active fault zones related to the earthquake. Several derivatives and structural parameters calculated from a gravity gradient tensor are applied to highlight the features. First horizontal and vertical derivatives, as well as a normalized total horizontal derivative, characterize well the continuous subsurface fault structure along the Futagawa fault zone. On the other hand, the Hinagu fault zone is not clearly detected by these derivatives, especially in the case of the Takano-Shirahata segment, suggesting a difference of cumulative vertical displacement between the two fault zones. The normalized total horizontal derivative and the dimensionality index indicate a discontinuity of the subsurface structure of the Hinagu fault zone, that is, a segment boundary between the Takano-Shirahata and the Hinagu segments. The aftershock distribution does not extend beyond this segment boundary. In other words, this segment boundary controls the southern end of the rupture area of the foreshock. We also recognize normal fault structures dipping to the northwest in some areas of the fault zones from estimations of dip angles.