Solubility measurements from multicomponent gas systems are still relatively scarce. The few experimental datasets involving binary gas systems found and gathered in this work show that PHREEQC and Reaktoro quantify reasonably well the co-solubility of CO2, N2, and CH4 from either CO2-N2 or CH4-CO2 binary gas mixtures, providing a good description of the equilibrium compositions in the aqueous and gaseous phase. To test PHREEQC and Reaktoro under conditions that are more representative of natural systems, these were compared against the limited bubble point pressure measurements reported in the literature. In this case, the aqueous solutions were not simply synthetic brines prepared in a laboratory, but samples obtained from geothermal wells. Considering the compositional assortment of natural waters, the agreement between the predicted and measured bubble point pressures varied from sample to sample; in some cases, the consistency between the predicted and measured dataset is fairly good (e.g., , whereas in some other cases, a greater discrepancy is observed. There are several factors, stemming either from the experimental protocols or model limitations, that may explain the discrepancies. Firstly, the gas samples at many sites consist of a mixture of many different components. Yet, we showed that no database can handle all gas types with the same accuracy. Thus, the two databases tested may not perfectly describe the fluid mixtures considered. Secondly, the measured bubble point pressure is also impacted by how it was determined. Different sampling, fluid handling, and measurement techniques may lead to different results. Moreover, the definition of the bubble point pressure may also drive to contrasting modelled and measured values. Theoretically, the bubble point pressure is the pressure at which the first bubble of gas appears within a liquid. However, experimentally, depending on the technique used, the onset of degassing may not be recorded with the same precision. In many cases, the bubble point pressure is determined using a visual criterion, whereas within the models, we always reported the Pb as the pressure at which the volume of gas is above 1×10 -5 L/ L of solution. Lastly, a mismatch can also be introduced simply because the modelled system is not representative of the real experimental conditions. As mentioned in the results section, some of the experimental datasets either do not clearly state which conditions the GLR is given at or contain flaws in the reported solution and gas composition, i.e., charge and mass imbalance, respectively. To deal with these obstacles, in the models, we ensured electroneutrality by adding/removing chloride and by normalizing gas compositions. Yet, the fidelity of the models and their predictability can only be achieved by a proper chemical (experimental) characterization of the fluids.