We study the dynamic orbit about Phobos of the Martian Moon eXploration (MMX) spacecraft and simulate gravity field estimation using Doppler, image landmarks, and LIght Detection And Ranging (LIDAR) data based on a mission plan of MMX and investigate whether the differences in the internal density structure of Phobos can be detected through this mission. Degree 2 Stokes coefficients of the Phobos gravity field, C¯20\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\overline{C} }_{20}$$\\end{document} and C¯22\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\overline{C} }_{22}$$\\end{document}, which are necessary to obtain the moment of inertia by combining with the libration amplitude, could be determined with an accuracy within the order of 0.1% using datasets from two-dimensional (2D) quasi-satellite orbits (QSOs). If observations from 3D-QSO could be realized at low altitude, coefficients up to degree and order 5 could be estimated, which could be used to detect regional density anomalies. Moreover, the observation data from the ascending trajectory after landing on the Phobos surface can be used to detect local density anomalies around the landing site, which could help interpret the origin of the samples and understand the internal structure of Phobos.Graphical