Isostructural zirconium-based metal–organic frameworks (Zr-MOFs) have attracted the attention of researchers because of their remarkable stability at high temperatures and high pressures and their chemical stabilities against acids and bases. Due to this stability, Zr-MOFs can be utilized in adsorption research, and the adsorption performance of a Zr-MOF depends on the pore size and the surroundings of the MOF. In this study, as the dimensions changed and the adsorption was carried out, the Zr-MOF material remained stable, and the adsorption of the best state was achieved at 235 mg/g. Through the simulation of theoretical kinetic models of Zr-MOFs, we initially postulated that the adsorption capacity is proportional to the pore size and that acid orange 7 (AO7) was adsorbed by the MOFs. Afterwards, we verified our hypotheses through a series of Brunauer–Emmett–Teller (BET) data analysis; non-local density function theory (NLDFT) was mainly used to analyze the data. Moreover, we determined that physical adsorption occurs on the surface of the MOFs during the adsorption process, while chemisorption occurs in the form of dye molecules combining with active sites. Ultimately, we concluded that the larger the pore size, the stronger the adsorption capacity, and this contribution casts a new light on the issue of wastewater treatment.