On account of the indeterminacy and subjectivity of decision makers (DMs) in complexity decision-making environments, the evaluation information over alternatives presented by DMs is usually fuzzy and ambiguous. As the generalization of intuitionistic fuzzy sets (IFSs), the Pythagorean fuzzy set (PFS) is more useful in expressing fuzzy and ambiguous information. Meanwhile, in order to consider human hesitance, dual hesitant Pythagorean fuzzy sets (DHPFSs) are presented, which can be more valid for handling real multiple attribute decision-making (MADM) problems. To fuse the information in DHPFSs more effectively, in this article, some dual hesitant Pythagorean fuzzy Heronian mean operators, which can consider the relationships between arguments being fused, are defined and studied. Evidently, the new proposed operators can obtain more exact results than other existing methods. In addition, some important properties of these Heronian mean (HM) operators are discussed. Subsequently, the defined aggregation operators are used in MADM with dual hesitant Pythagorean fuzzy numbers (DHPFNs), and the MADM model is developed. In accordance with the defined operators and the built model, the dual hesitant Pythagorean fuzzy generalized weighted Heronian mean (DHPFGWHM) operator and dual hesitant Pythagorean fuzzy generalized geometric weighted Heronian mean (DHPFGGWHM) operator are applied to deal with the green supplier selection in supply chain management, and the availability and superiority of the proposed operators are analyzed by comparing them with some existing approaches. The method presented in this paper can effectively solve the MADM problems in which the decision-making information is expressed by DHPFNs and the attributes are interactive.