The recognition of the geomechanical properties of methane hydrate-bearing soil (MHBS) is crucial to exploring energy resources. The paper presents the mechanical properties of a pore-filled MHBS at a critical state using the distinct element method (DEM). The pore-filled MHBS was simulated as cemented MH agglomerates to fill the soil pores at varying levels of methane hydration (MH) saturation. A group of triaxial compression (TC) tests were conducted, subjecting MHBS samples to varying effective confining pressures (ECPs). The mechanical behaviors of a pore-filled MHBS were analyzed, as it experienced significant strains leading to a critical state. The findings reveal that the proposed DEM successfully captures the qualitative geomechanical properties of MHBS. As MH saturation increases, the shear strength of MHBS generally rises. Moreover, higher ECPs result in increased shear strength and volumetric contraction. The peak shear strength of MHBS increases with rising MH saturation, while the residual deviator stress remains mainly unchanged at a critical state. There is a good correlation between fabric changes of the MHBS with variations in principal stresses and principal strains. With increasing axial strain, the coordination number (CN) and mechanical coordination number (MCN) increase to peak values as the values of MH saturation and ECPs increase, and reach a stable value at a larger axial strain.