Methane hydrate bearing sediments(MHBS) are usually found undersea, which is generally categorized into different types, i.e. bond, pore-filling and so on. The mechanical properties of MHBS will be deteriorated as methane hydrate (MH) dissociates and many marine accidents will be caused; thus a lot of attention has been attracted in investigating MHBS properties changes during exploitation. A newly proposed 2D temperature-water pressure-mechanical bond model is implemented into PFC2D, the commercial software of the distinct element method (DEM), to simulate the dissociation of MH by thermal recovery and depressurization using soil in good exhaust and drainage conditions. By comparison with laboratory test results, it is validated that the newly proposed model is applicable to simulate the mechanical behavior of MHBS. At the same time, the micromechanical properties are also analyzed. With temperature increasing during thermal exploitation, the anisotropic degree of total contact distribution increases; the amount of bond contact, whose principal direction is horizontal all the time, decreases and the amount of unbonded contact, whose principal direction is vertical all the time, increases; the value of averaged pure rotation rate (APR) and its concentration degree increases all the time. With back (water) pressure decreasing during exploitation by depressurization, the total contact distribution changes from isotropic to anisotropic with a vertical principal direction, and the value of APR is small and uniformly distributed. After the recovery of back water pressure, the sample is damaged further, and the anisotropic degree of total contact distribution increases. Besides, the value of APR increases and the concentration degree of positive and negative value increases.