An Extensive Study on Desorption Models Generated Based on Langmuir and Knudsen Diffusion

: Although gas desorption is a known phenomenon, modeling ﬂuid ﬂow in tight gas reservoirs often ignores the governing desorption effect, assuming that viscous transport is the predominant controller, resulting in an erroneous prediction of mass transport and ﬂuid ﬂow calculations. Thus, developing a new model accommodating all the major contributing forces in such a medium is essential. This work introduces a new comprehensive ﬂow model suitable for tight unconventional reservoirs, including viscous, inertia, diffusion, and sorption forces, to account for ﬂuid transport. Based on Langmuir law and Knudsen diffusion effect, three models were generated and compared with different known models using synthetic data. The model was solved and ana-lyzed for different scenario cases, and parametric studies were conducted to evaluate the desorption effect on different reservoir types using MATLAB. Results show that the contribution of the sorption mechanism to the ﬂow increases with the reducing permeability of the medium and lower viscosity of the ﬂowing ﬂuid and an additional pressure drop up to 10 psi was quantiﬁed.


Introduction
Unconventional reservoirs embody the upcoming revolution in the oil and gas industry due to the rapid growth of world energy demand anticipated by the United States Energy Information Association Outlook [1] and the continuous draining of the current conventional hydrocarbon resources. The increase in demand has led to an aggressive investigation of unconventional accumulations that were once shelved due to production complexities to be considered for exploring opportunities. Unconventional reservoirs are reservoirs with divergent geological characteristics [2], inconstant geochemical characteristics, intricate petrophysical properties, eccentricities in fluid phase behavior, and various governing flow mechanisms/forces [3]. Unconventional reservoirs include various types of reservoirs including gas hydrates, tight gas and oil, heavy oil, and shale oil and gas.
In their natural form, gas hydrates consist of water molecules (ice molecules) that act as hosts [4]. In contrast, guest molecules are methane, propane, isobutene, and ethane, among other gases that are chemically bonded with van der Waals forces in low degrees of temperature and, at the same time, under high pressure. Natural hydrate gas is primarily white and has the same appearance as ice. Natural hydrates are commonly known to have methane, which can burn, thus having the name "fire ice". Natural gas hydrates have long been deemed a nuisance for blocking transmission pipes, endangering the foundations of deep-water platforms and pipelines, and risking more disruption to the output of deepwater oil and gas [5], maybe a significant possible source of energy in the future. In the permafrost and beneath the ocean floor, huge deposits of natural gas hydrates are widely dispersed [6].
Geologically, shale is the source rock where oil was originated and trapped with part migration to the reservoir rock where migration path was available. Shale oil is mainly found in organic shale with nanoscale pore diameters [7,8]. The shale rock has two purposes: one is generating while the other is reserving shale oil. Shale oil exists in both adsorbed form and free state form. Generally, shale oil is lightweight with low viscosity levels as well. Due to the very low transmissibility of reservoir rock, the oil is hard to produce on a commercial scale. The term "shale oil" also applies to oil found in semi-previous shale where permeability is in micro-Darcies. In contrast to shale gas, shale oil was formed during the generative window of organic evolution. As the shale rock continues to produce, a large portion of shale oil is retained, thus making the reservoir saturated.
Coalbed methane (CBM), as the name implies, is primarily composed of methane generated in coal seams and stored in the adsorbed state. Its common name in coal mining is "marsh gas". Coalbed methane is more or less an unconventional gas that is mostly found in coal beds where coals are both the source rock as well as the reservoir. Ninety-five percent of coalbed methane is chiefly composed of methane and other dispensable amounts of hydrocarbons and nonhydrocarbons.
Tight oil corresponds to oil trapped in reservoirs where low or ultra-low rock permeability contributes to unfavorable reservoir flow efficiency. Rock matrix permeability is mostly of the range of 10 to 1 mD or less and has a porosity that ranges between 1% and 10%. For conventional growth, some sandstone and carbonate reservoirs may have even lower permeabilities [9]. Tight sandstone gas refers to gas originated in sandstone, which has a permeability matrix that is less than 10 −1 mD and is characterized by a porosity level below 10%. Sandstone interbedded with shale-mudstone, resulting in reservoirs immediately in contact with sources, no obvious traps or direct caprocks, but regional seals are well developed [10].
The forces found in unconventional reservoirs that control the fluid movement are diffusion, viscous, convection, desorption, inertial, capillary, sorption, and viscoelastic. These mechanisms are classified into two groups based on their controlling effect: trapability and/or displacement [11]. Trapability forces aid in trapping the wetting phase fluid, while displacement controls fluid mobility through the porous medium in nano, micro, and macro scales. Note here that these mechanisms can become of importance when discussing fluid flow and other relevant issues related to CO 2 and hydrogen storage as the model that will be created later on involves a term that represents the sorption phenomena [12].

Governing Forces
Several physical processes and mechanisms affect the flow in porous media. The effect of certain forces governs more than other forces depending on the reservoir types and conditions and the forces acting as trapability or displacement [11]. Thus, the currently used method-Darcy's equation-calculates pressure drop and flow concentration, but it lacks accuracy and validity as it ignores most of the governing forces. Diffusion is the movement of particles or materials (gas or liquid) in a substance from a high concentration area to a low concentration area. It is considered a displacement force [13]. Viscoelasticity describes the property of materials that can be either solids or liquid experiencing deformation forces, where both viscous and elastic characteristics forces can be observed [14,15]. The viscoelasticity force is composed of both displacement and trapability forces. Capillary pressure is a force that exists in saturated reservoirs with two or more immiscible fluids in reservoir rock's tiny pores [16]. The presence of capillary forces in a porous medium leads to hydrocarbon entrapment. Capillary forces are considered trapability forces due to their mechanism [17]. Inertial forces form due to the fluid momentum [18]. The fluid tends to resist velocity change by applying an equal or opposite force to the direction of the applied force. Inertial force is one of the forces affecting the displacement of fluid in unconventional reservoirs. Advection forces contribute to fluid and molecular movement through porous media, enhancing displacement and mass transfer [19]. Viscosity is the resistance of fluids to any alternations in shape or movement of nearby segments [18]. It could be described as internal friction within the fluid molecules, which resists velocity Energies 2021, 14, 6435 3 of 15 difference advancement. The viscosity of liquids decreases significantly as temperature increases, whereas the viscosity of gases increases as temperature increases. Sorption and desorption: by definition, sorption is the transfer of a substance from the gas phase or solution to the solid phase. Sorption has two types: adsorption and absorption [19]. Absorption occurs when molecules penetrate a solid medium. In contrast, adsorption occurs when gaseous molecules get attached to the surface of the solid and depend on the properties of the sorbate (gaseous molecules) and the sorbent (solid) [20]. Sorption is affected by many factors: pressure, temperature, porosity, permeability, the surface area of sorbent, and organic carbon content of the sorbate [20]. On the other hand, desorption occurs when reservoir pressure drops and the free gas that got sorbed by the sorbent gets released again [21][22][23].
Fluids can coexist in the same basin in the subsurface, forming reservoirs with two or more fluid types, such as oil and gas; oil and water; and water, oil, and gas reservoirs. Such cases led to further investigations to improve Darcy's equation to accommodate the multiphase flow phenomena. Consequently, a comprehensive and intensive study utilized different resources, theoretical research, experimental observation, data processing, simulation, and field applications. As a result, many publications have been produced to describe and shape the multiphase flow in reservoirs that enriched the petroleum engineering research outcomes [24,25]. The measurement of oil, gas, and water flow; output; and assessment of primary and secondary oil recovery and reservoir dynamics focus on multiphase flow physics in reservoirs [26]. Due to the nature of Darcy's equation and its limitations, it is replaced by the Forchheimer equation at higher flow velocities [27].

Forchheimer's Equation
As Darcy's equation has limited validity in the low-velocity domain, a better flow correlation to estimate higher velocity flow was crucially required. Forchheimer's equation was one of the first equations that addressed this issue and improved fluid flow prediction and calculations. Forchheimer's equation is an empirical equation that relates pressure drop across a given distance with fluid velocity and friction across porous media. Forchheimer's equation is a modified extension of Darcy's equation to accommodate, describe and calculate the effect of non-Darcy flow in porous media [28].
Forchheimer's equation is expressed as where β is the non-Darcy flow coefficient, which is measured using core samples in laboratories or determined from analysis of multi-rate pressure tests. The first term represents Darcy's equation, which denotes the viscous forces, while the second term denotes the inertial governing forces. Deviations in fluid flow behavior from Darcy's law have long been recognized in many applications [29,30]

Langmuir Adsorption Isotherm
In 1916, Irvin Langmuir developed the theory of Langmuir isotherm, which was used extensively to describe the sorption/desorption in the oil and gas industry. The theory is based on the following assumptions:

•
One molecule is adsorbed in each adsorption location.

•
No interaction between adsorbed gas molecules within the pores.

•
The energy at sorption/desorption location is equal.
Langmuir isotherm equation is presented as follows: where V is the volume of the adsorbed gas, V L is the Langmuir volume constant, P L is the Langmuir pressure constant, and P is the pressure at the reservoir as shown in Figure 1.
where is the non-Darcy flow coefficient, which is measured using core samples in laboratories or determined from analysis of multi-rate pressure tests. The first term represents Darcy's equation, which denotes the viscous forces, while the second term denotes the inertial governing forces. Deviations in fluid flow behavior from Darcy's law have long been recognized in many applications [29,30] 1. 3

. Langmuir Adsorption Isotherm
In 1916, Irvin Langmuir developed the theory of Langmuir isotherm, which was used extensively to describe the sorption/desorption in the oil and gas industry. The theory is based on the following assumptions: • One molecule is adsorbed in each adsorption location.

•
No interaction between adsorbed gas molecules within the pores.

•
The energy at sorption/desorption location is equal.
Langmuir isotherm equation is presented as follows: where is the volume of the adsorbed gas, is the Langmuir volume constant, is the Langmuir pressure constant, and is the pressure at the reservoir as shown in Figure  1.

Knudsen Diffusion
Knudsen diffusion, named after Martin Knudsen, occurs when the mean pore diameter is smaller than the mean free path of gas molecules, leading to collisions of gas molecules with the pore walls. Knudsen diffusion is considered a predominant transport mechanism.
The concentration flux can be calculated using the following equation: in which, is the Knudsen diffusivity and can be calculated using the following equation:

Knudsen Diffusion
Knudsen diffusion, named after Martin Knudsen, occurs when the mean pore diameter is smaller than the mean free path of gas molecules, leading to collisions of gas molecules with the pore walls. Knudsen diffusion is considered a predominant transport mechanism.
The concentration flux can be calculated using the following equation: (4) in which, D k is the Knudsen diffusivity and can be calculated using the following equation: where d p is the pore diameter, T is temperature, M is molar mass of gas, and k is the Boltzmann constant. The present work aims to establish a new model based on Langmuir and Knudsen diffusion to better estimate the fluid flow in tight porous media. Additionally, it also targets to couple diffusion model to capture the actual behavior of the fluid. Several investiga- tions were conducted using various case matrices and data from different unconventional reservoir types using MATLAB (developed by MathWorks in Massachusetts, USA) for simulation [22,31,32].

Model Creation
Several steps were taken to create the desired models, and a thorough comparison was conducted to determine the optimum model that quantitively represents viscous, inertial, diffusion, and desorption forces in a tight porous medium. The same concept that Forchheimer applied is implemented in this study to formulate a new equation that accounts for viscous, inertial, diffusion, and desorption forces. Note that Forchheimer's Model represents viscous and inertial forces only, and thus a new term that can represent both diffusion and desorption forces must be added.
The first step to formulate the desired model is to decide which Desorption/Diffusion law will be used to create the model and how it can be developed to fit the porous media characteristics. Three basic fundamental laws were used: Fick's law to represent diffusion forces, Knudsen law used in two primary forms: one to represent diffusion forces and the other to represent desorption forces, and Langmuir's law to represent desorption forces. Three models were developed in this study, and each will be discussed in detail in the following sections.

First Model
To create a credible model, a relationship between diffusion and desorption forces and the pressure drop should be established to add the desired term to Forchheimer's equation. The first term created is expected to capture the effect of both diffusion and desorption forces by coupling these two forces in a single term based on the understanding of Fick's and Langmuir's models. Note that Fick's first law was used to establish such a relationship as it states that the velocity is directly proportional to the pressure gradient. By understanding the merits of Fick's first law, it can be concluded that the concentration gradient controls diffusion. Furthermore, Langmuir's coefficient represented by the symbol θ is added to mathematically represent desorption forces as in Equation (8).
To conclude the addition of the following term developed the first model was created: We perform a dimensional analysis to check the term's validity: The second model was similarly created by adding a single term that represents both diffusion and desorption forces. It was generated using Fick's diffusion, Knudsen desorption, and Langmuir's models. The main difference between this model and the first model created is the added Knudsen desorption coefficient to enhance the model's accuracy and represent desorption forces better. It is believed that desorption forces cannot be represented by Langmuir's coefficient only. Similarly, the second model was developed by the addition of the following term:

Third Model
The addition of two individual terms-one representing diffusion forces and the second representing desorption forces-distinctively created the third model. The diffusion term, which was initially developed by AlDhuhoori [33], was represented by a velocity term based on the original concept of Fick's partial pressure law, a diffusivity coefficient which was derived from Fick's first law was included to increase the accuracy of the model and to account for the diffusion of a particular phase into another. Furthermore, a non-Darcy term that includes both porosity and permeability of a medium was added to emulate reservoir properties and characteristics.
As for the separate desorption term, generated by understanding the merits of Knudsen desorption law, it includes a viscosity term, a velocity term, a length term that represents the length of the reservoir (since it plays a significant role in estimating Knudsen desorption coefficient, which was also added to the new term). The new addition to Forchheimer's model can be represented mathematically by adding two terms represented in the following equation: ∂P ∂x

Mathematical Formulation
The same concept that Forchheimer applied is implemented in this study to formulate a new equation that accounts for viscous, inertial, diffusion, and desorption forces. The new modifications to the model address the fluid transport phenomena into three scales incorporating a diffusion-desorption term while accounting for the porous media properties. A detailed derivation of the second model is shown in this section to explain the mathematical expression and steps taken to reach the final form of the model and to achieve the main goal of the derivation: to be able to represent these forces in both time and space. The first step is to add the modified diffusion/desorption term to the equation, and thus we get Taking the derivative of both sides of the equation with respect to x: By rearranging we get This is the continuity equation in 1-D form for a single steady-state system, or By substituting the continuity equation from above in the previous equation, we get Using the chain rule, ∂ρ ∂t We get Moreover, taking ∂ρ ∂P = ρc We get Or on simplifying and rearranging which can be written as Note that the other two models were formulated using almost the same steps.

Numerical Modeling
A unique synthetic environment consisting of 1-D mass transport of a single gaseous phase has been created to study the validity of various models generated above.
The conceptual design of model has the following assumption: reservoir characteristics and features were adequately selected to ensure that parametric values properly represent the reservoir rock and fluid properties. Proper boundaries were selected, and the models were then discretized evenly to different blocks for numerical simulation. The central point of the block method was selected for better estimation of pressure and concentration solutions, using the implicit method of second-order partial derivatives, central and forward methods for the first order of partial derivatives. The following equations were generated and substituted in the generated model to solve via MATLAB Simulator: Finite difference formula for ∂P ∂t : Finite difference formula for ∂P ∂x : Finite difference formula for ∂ 2 P ∂x 2 : Finite difference formula for ∂ 2 C ∂x 2 : An iterative approach was employed, where at each time step, the model calculates concentration based on the pressure values. As for the pressure estimation, the model calculates the pressure at each time step based on the pressure estimated at the previous time step, considering the reservoir characteristics and fluid properties during the process. All the numerical models are set up to follow the same process and iterative approach.

Synthetic Case Study Validation
The synthetic case study shows a gaseous fluid with two slightly different concentrations and a reservoir with two main different boundary pressure conditions. The reservoir rock parameters are assumed to be constant, and thus this case is entirely homogenous. Additionally, fluid properties are assumed to be constant.
Furthermore, a grid system of 100 grid blocks is used to represent the 500 ft long reservoir. Table 1 summarizes the data used in the gas-gas phase medium containing mainly methane [32][33][34][35]. Several plots have been constructed in MATLAB to simulate the behavior of a tight unconventional reservoir for both pressure and concentration drop during a specific time and distance until pressure and concentration steady-state times are reached. The newly derived models were tested for validation by utilizing the data provided in Table 1, representing a tight unconventional reservoir gas-gas case. Figure 2 shows the results of the Desorption M3 model, created to represent the pressure distribution within a porous media, the concentration, and the rest of the model's pressure profiles results. Additionally, the pressure gradient and the concentration gradient at different time iterations in terms of hours until we reach a pressure steady-state time and a concentration steady-state time. Note that calculations showed that the time required for concentration to reach the steady-state condition for each model is less than the time required for the pressure to reach the steady-state phase for all the models.  Several plots have been constructed in MATLAB to simulate the behavior of a tight unconventional reservoir for both pressure and concentration drop during a specific time and distance until pressure and concentration steady-state times are reached. The newly derived models were tested for validation by utilizing the data provided in Table 1, representing a tight unconventional reservoir gas-gas case. Figure 2 shows the results of the Desorption M3 model, created to represent the pressure distribution within a porous media, the concentration, and the rest of the model's pressure profiles results. Additionally, the pressure gradient and the concentration gradient at different time iterations in terms of hours until we reach a pressure steady-state time and a concentration steady-state time. Note that calculations showed that the time required for concentration to reach the steady-state condition for each model is less than the time required for the pressure to reach the steady-state phase for all the models.

Comparison of the Newly Developed Models Pressure Gradients
By comparing the pressure gradients of the three newly developed models, it can be observed that the velocity and the pressure gradient have a directly proportional relationship ( Figure 3). As the velocity increase, so does the pressure gradient. Additionally, Desorption Model 1 has a much higher pressure gradient estimation compared to the other two models created, which showed almost the same result/trend with a pressure differential of almost ten psi. The initial presumption observed at this stage is that the first model is overestimating the pressure differential compared to the other two models despite being theoretically accurate. In the following section, a comparison will be carried out to investigate the credibility of all three models by comparing them to previously published desorption models.

Comparison of the Newly Developed Models Pressure Gradients
By comparing the pressure gradients of the three newly developed models, it can be observed that the velocity and the pressure gradient have a directly proportional relationship ( Figure 3). As the velocity increase, so does the pressure gradient. Additionally, Desorption Model 1 has a much higher pressure gradient estimation compared to the other two models created, which showed almost the same result/trend with a pressure differential of almost ten psi. The initial presumption observed at this stage is that the first model is overestimating the pressure differential compared to the other two models despite being theoretically accurate. In the following section, a comparison will be carried out to investigate the credibility of all three models by comparing them to previously published desorption models.

Comparison of Different Flow Models with the Newly Developed Models
A comparison between Darcy's model, Forchheimer's model, diffusion, and the newly developed model in the gas-gas case is performed to compare the pressure gradient drop using the three models. Figure 4 shows the predicted pressure drop across the porous medium at steady-state conditions using the models. Where the new models should accommodate the following forces: inertial loss, diffusion, and desorption.

Comparison of Different Flow Models with the Newly Developed Models
A comparison between Darcy's model, Forchheimer's model, diffusion, and the newly developed model in the gas-gas case is performed to compare the pressure gradient drop using the three models. Figure 4 shows the predicted pressure drop across the porous medium at steady-state conditions using the models. Where the new models should accommodate the following forces: inertial loss, diffusion, and desorption.   Figure 4 shows the newly developed model, which exhibited a semi-exponential behavior. Compared to the diffusion model created by Aldhuhoori [33], Forchheimer's model, and Darcy's model, Desorption M1 experiences a high-pressure gradient estima-  Figure 4 shows the newly developed model, which exhibited a semi-exponential behavior. Compared to the diffusion model created by Aldhuhoori [33], Forchheimer's model, and Darcy's model, Desorption M1 experiences a high-pressure gradient estimation of a difference of 66.402 psi with Diffusion model, a difference of 71.25 psi compared with Forchheimer's, and a difference of 142.650 psi with Darcy's. The first model overestimates the pressure drop, which indicates that this model failed to quantify the effect of both desorption and diffusion. On the other hand, Desorption M2 lays precisely on Forchheimer's Models as the difference of pressure between Desorption M2 and Forchheimer's is zero as shown in Table 2, leading to lower pressure drop estimation for desorption and diffusion effects and hence underestimating the pressure drop. Moreover, the Desorption M3 model represented viscous, inertia, diffusion, and desorption effects and resulted in a higher pressure gradient than the diffusion model. So far, the third desorption model is showing promising results. It can capture the effects of both diffusion and desorption without overestimating/underestimating compared to the other two models.

Comparison of Newly Developed Models with Well-Known Published Models
Furthermore, the newly developed model was compared with two different wellknown publication models to study the sorption effect only. The two models were developed by Shabro [36] and Boosari [37].
It can be observed in Figure 5 that compared with Shabro's and Boosari's models, Desorption M3 generates a more accurate pressure gradient compared with the other desorption models. In comparison, Desorption M1 has an overestimating pressure gradient, while Desorption M2 model is underestimating the pressure gradient and hence the first and second desorption models will be discarded, and further studies will be conducted to thoroughly investigate the effect of the third desorption model and its effectiveness in estimating the pressure drop.

Comparison of Newly Developed Models Using Different Permeabilities
Over the years, many definitions and methods to distinguish unconventional reservoirs from conventional reservoirs were developed by many researchers [38,39]. They mainly focused on a certain property or a specified characteristic to distinguish between those reservoirs. A couple of the most famous criteria used to characterize and distinguish

Comparison of Newly Developed Models Using Different Permeabilities
Over the years, many definitions and methods to distinguish unconventional reservoirs from conventional reservoirs were developed by many researchers [38,39]. They mainly focused on a certain property or a specified characteristic to distinguish between those reservoirs. A couple of the most famous criteria used to characterize and distinguish reservoirs are based on the permeability of the porous media in a specified reservoir. In this sensitivity study, the newly developed model will be tested against different permeabilities representing different types of reservoirs to examine the effectiveness of the desorption forces utilizing a criterion developed by Belhaj [40].
Belhaj classified various types of reservoirs based on permeability values ranging from 10 −4 mD up to 100 mD. It includes five main types of reservoirs: shale, tight reservoirs, limestone, CBM, and conventional reservoirs as summarized in Table 3. To investigate the effect of permeability on the newly developed model, Darcy's model, Forchheimer's, and AlDhuhoori's diffusion model are included in this study. Figure 6 shows the criteria developed by Belhaj to classify different types of reservoirs based on reservoir permeability.  The four models were compared at different permeabilities ranging from 10 −4 mD up to 10 −1 mD. As indicated from Figure 7 as the velocity increases, the difference between the pressure drop estimation using the newly developed model and Forchheimer's model becomes clear. The diffusion and desorption effects become much more noticeable. The four models were compared at different permeabilities ranging from 10 −4 mD up to 10 −1 mD. As indicated from Figure 7 as the velocity increases, the difference between the pressure drop estimation using the newly developed model and Forchheimer's model becomes clear. The diffusion and desorption effects become much more noticeable. It can be concluded that as the permeability increases, the difference between the pressure drop calculated using the newly developed model and Forchheimer's model decreases until both models generate the same result at a permeability of 10 −1 mD and thus indicating that the diffusion and desorption effects decrease as the permeability increases. Thus, it is apparent that the newly developed model can adequately capture the effect of diffusion and desorption forces in extremely tight media such as shale, tight reservoirs, and some limestones and fails to characterize that effect in higher permeability reservoirs such as CBM and conventional reservoirs in general. Note that at low permeabilities, the model exhibited the same trend as that of Forchheimer's, but at higher velocities, the model fails to capture those effects. Therefore, the model can only predict the pressure drop at low permeabilities.
It can also be observed that the desorption effects tend to have exponential behavior with low permeabilities, while with higher permeability, it tends to have a straight line behavior.

Conclusions
The following points summarize the work done on investigating the effect of diffusion forces in modeling fluid flow in tight unconventional reservoirs and highlights the main conclusions: • A new 1-D linear and radial mathematical models were derived, tested, and numerically simulated to investigate and address the effect of diffusion/desorption forces in tight unconventional reservoirs.

•
The newly developed model proved the importance of the contribution of diffusion/desorption forces to fluid flow in the porous media, as initially hypothesized. The significant effect of both forces is noted as the reservoir deviated towards unconventionality.

•
Like the establishment of steady-state time for pressure gradient, the concept of steady-state time for concentration gradient is introduced for the first time. This newly introduced concentration gradient steady-state time means that the medium It can be concluded that as the permeability increases, the difference between the pressure drop calculated using the newly developed model and Forchheimer's model decreases until both models generate the same result at a permeability of 10 −1 mD and thus indicating that the diffusion and desorption effects decrease as the permeability increases. Thus, it is apparent that the newly developed model can adequately capture the effect of diffusion and desorption forces in extremely tight media such as shale, tight reservoirs, and some limestones and fails to characterize that effect in higher permeability reservoirs such as CBM and conventional reservoirs in general. Note that at low permeabilities, the model exhibited the same trend as that of Forchheimer's, but at higher velocities, the model fails to capture those effects. Therefore, the model can only predict the pressure drop at low permeabilities.
It can also be observed that the desorption effects tend to have exponential behavior with low permeabilities, while with higher permeability, it tends to have a straight line behavior.

Conclusions
The following points summarize the work done on investigating the effect of diffusion forces in modeling fluid flow in tight unconventional reservoirs and highlights the main conclusions:

•
A new 1-D linear and radial mathematical models were derived, tested, and numerically simulated to investigate and address the effect of diffusion/desorption forces in tight unconventional reservoirs.

•
The newly developed model proved the importance of the contribution of diffusion/desorption forces to fluid flow in the porous media, as initially hypothesized. The significant effect of both forces is noted as the reservoir deviated towards unconventionality.

•
Like the establishment of steady-state time for pressure gradient, the concept of steady-state time for concentration gradient is introduced for the first time. This newly introduced concentration gradient steady-state time means that the medium reached a constant concentration gradient throughout the pores system (constant slope).

Patents
A new desorption model has been created that better estimates the desorption effect along with the diffusion effect. The new proposed model will be able to address the pressure drop in porous media by adhering to the diffusion and desorption effect.