As a key resource in the oil and gas sector, shale gas development profoundly influences the advancement of the global energy industry. Deep shale gas reservoirs, typically found at depths exceeding 3500 m, represent a significant portion of total shale gas reserves. Currently, pore network models primarily simulate middle and shallow shale gas, insufficiently addressing the unique challenges posed by high-temperature and high-pressure conditions in deep shale gas formations. There is an urgent need to explore the microscale flow dynamics of deep shale gas. In this work, we use the pore network model to simulate the flow dynamics of deep shale gas, particularly under nanoconfinement conditions. The simulation integrates adsorption, slip flow, Knudsen diffusion, and bulk flow phenomena. Utilizing the dual-site Langmuir adsorption model tailored for high-temperature and high-pressure conditions enhances the accuracy of gas conductivity calculations for the adsorbed phase, thereby reflecting the specific characteristics of deep shale gas reservoirs. Moreover, this research investigates the flow characteristics of deep shale gas under varying sensitivity parameters, such as water saturation, temperature, and pressure. It explores methane flow dynamics, focusing on effective diffusion coefficients and permeability. The modified approach to calculating adsorption phase conductivity using the dual-site Langmuir adsorption model accurately captures the adsorption behaviors and characteristics of deep shale gas reservoirs.