Natural　gas　hydrate　(NGH)　is　considered　as　a　type　of　clean　energy　to　replace　coal　and　oil.　During　exploitation,　permeability　is　one　of　the　key　parameters　controlling　production　efficiency,　reservoir　stability,　and　greenhouse　gas　sequestration.　Limited　by　experimental　and　numerical　simulation　tools,　in　current　research,　the　directionality　of　permeability　is　usually　ignored.　In　this　work,　a　DEM-CFD　coupling　simulation　method　is　developed　to　compute　the　anisotropic　permeability.　The　sedimentary　process　of　reservoir　sediments　is　reconstructed,　enabling　the　acquisition　of　numerical　models　that　possess　pore　structures　consistent　with　the　actual　fabric　characteristics.　The　fluid　transport　process　in　various　directions　can　be　simulated　with　a　finite　element　method.　Taking　the　natural　gas　hydrate　reservoir　in　the　Shenhu　area　of　the　South　China　Sea　as　an　example,　the　proposed　method　is　validated　and　applied　to　explore　the　effect　of　compaction　stress　on　permeability　anisotropy.　With　the　increase　in　compaction　stress,　the　permeability　anisotropy　exhibited　a　rapid　initial　increase,　followed　by　a　sustained　stabilization.　The　primary　cause　is　the　rearrangement　of　sediment　particles.　The　non-spherical　particles　are　driven　to　align　in　a　predominantly　horizontal　orientation,　thereby　enhancing　anisotropy.　The　proposed　method　provides　a　tool　for　the　efficient　exploitation　of　hydrate　resources.