Progress of Seepage Law and Development Technologies for Shale Condensate Gas Reservoirs

: With the continuous development of conventional oil and gas resources, the strategic transformation of energy structure is imminent. Shale condensate gas reservoir has high development value because of its abundant reserves. However, due to the multi-scale ﬂow of shale gas, adsorption and desorption, the strong stress sensitivity of matrix and fractures, the abnormal condensation phase transition mechanism, high-speed non-Darcy seepage in artiﬁcial fractures, and heterogeneity of reservoir and multiphase ﬂows, the multi-scale nonlinear seepage mechanisms are extremely complicated in shale condensate gas reservoirs. A certain theoretical basis for the engineering development can be provided by mastering the percolation law of shale condensate gas reservoirs, such as improvement of productivity prediction and recovery efﬁciency. The productivity evaluation method of shale condensate gas wells based on empirical method is simple in calculation but poor in reliability. The characteristic curve analysis method has strong reliability but a great dependence on the selection of the seepage model. The artiﬁcial intelligence method can deal with complex data and has a high prediction accuracy. Establishing an efﬁcient shale condensate gas reservoir development simulation technology and accurately predicting the production performance of production wells will help to rationally formulate a stable and high-yield mining scheme, so as to obtain better economic beneﬁts.


Introduction
With the reduction of conventional oil and gas resources and the strategic transformation of energy structure, the sustainable development of shale gas reservoir has the potential to fully meet the growing energy demand.Shale condensate gas reservoir plays an increasingly important role in the development of unconventional gas resources in the world.This resource target is relatively new and has great development potential.Its exploitation technology is challenging, which is mainly characterized by incorporating complicated multi-scale seepage law involved in the development of shale condensate gas reservoirs.Currently, great research progress on the multi-scale seepage law and the exploitation technology has been made in recent years.However, some difficult research problems remain to be solved in the future, such as the limitation of experimental equipment and nonupgraded research methods.In this paper, novel theories of seepage law emerging, experimental methods, and exploitation techniques of shale condensate gas reservoirs are reviewed and discussed.The paper is not only very helpful to provide ideas and guidance in time for upcoming difficulties in research and technology development, but also provides technical support for exploitation strategy innovation, so as to improve shale gas recovery efficiency and maximize economic benefits of shale condensate gas reservoirs.
In the process of development, the better expansibility of shale condensate gas reservoir can alleviate the reduction of reservoir pressure.Moreover, some condensate gas reservoirs exist in the liquid phase and have strong fluidity.Hence, condensate gas reservoirs usually have a high output.That is, in actual engineering, the transition area of condensate and gas have the highest output, the best economic benefits, and better development prospects [1].However, when the pressure drops below the dew point pressure, the condensate will be blocked.The flow resistance of gas phase and liquid phase will increase, and the productivity of producing wells will decrease.Moreover, the matrix porosity and permeability of shale condensate gas reservoir are very low, with micro/nanopores.The diameter and distribution of pores will affect the change of the fluid phase state, and also increase the critical condensate oil saturation.That is, the minimum saturation of condensate oil flow increases, and the recovery rate of condensate gas reservoir decreases.Natural gas production and liquid recovery will be significantly reduced with the reduction of oil-gas two-phase flow or even the three-phase flow of oil, gas, and water under water-bearing conditions [2][3][4].
The multi-scale nonlinear seepage law of shale condensate gas reservoir is complicated.First, in micro/nanopores, gas flow has Knudsen diffusion, molecular diffusion, and other multi-scale flow characteristics.Second, there is an adsorption-desorption mechanism.Furthermore, shale matrix and fractures have strong stress sensitivity.Finally, with the decrease in reservoir pressure, the formation pressure of shale condensate gas reservoir is gradually lower than the dew point pressure, and the condensate begins to precipitate.When the local reservoir pressure is lower than the critical flow pressure, the precipitated condensate flows near the production wellbore, forming the oil-gas two-phase flow or even the oil-gas-water three-phase flow under the condition of water.This special phase behavior may have an important impact on reservoir development.Shale condensate gas reservoirs are very dense, with an extremely low permeability and difficult fluid flow, and they are usually developed by large-scale hydraulic fracturing, which leads to the existence of multi-scale pore space composed of a tight matrix, microfractures, and hydraulic fractures [1].Therefore, it is particularly important to master the multi-scale complex seepage law of shale condensate gas reservoir development based on multiphysical field coupling, so as to establish an efficient simulation technology for shale condensate gas reservoir.
The complete production process of production wells in shale condensate gas reservoirs includes two parts: early fracturing flowback stage and long-term production stage, and their production rules are different.In the long-term production stage, fracturing fractures may be closed, the adsorbed gas may be desorbed, and condensate oil is released, all of which will affect the productivity of production wells [2][3][4].The above factors should be taken into account in the accurate production evaluation method of shale condensate gas reservoir.The establishment of an efficient simulation technology for shale condensate gas reservoir development can accurately predict the production performance of production wells, which is helpful to make reasonable and stable production plans with a high yield, so as to obtain better economic benefits.

Reservoir Characteristics, Resource Distribution, and Exploitation Methods for Shale Condensate Gas Reservoirs
Shale gas usually exists in the adsorbed state and free state.Adsorbed shale gas usually exists in organic matter.About 20%-82% of shale gas exists in the adsorbed state, and different shales have different adsorbed gas ratios, with an average of about 50% [5].Free shale gas usually exists in pores.Rock porosity is an important parameter for determining free gas content.Shale is very dense, with a small pore size and a large specific surface area, which is the main space for gas occurrence [6].The pore types are diverse and the structures are complex.The pore structure of shale condensate gas reservoir is quite different from that of dry gas reservoir.In shale condensate gas reservoir, there are mainly organic matter (OM) pores, clay mineral (CM) pores, intragranular dissolution (ID) pores, and microfractures [7].Micropores and mesoporous pores are mainly provided by OM and CM pores, while macropores are mainly provided by CM pores and microfractures [7].
The shale gas development of China has been led by Sinopec and PetroChina.The Jurassic shale in Sichuan Basin is characterized by abundant resources and great exploitation prospects.Hu et al. found that the first shale condensate gas reservoir in China developed under an abnormally high pressure and low geothermal system in the Jurassic strata of the Fuxing area in eastern Sichuan [8].The shale condensate reservoir is buried at a moderate depth of 2200-2950 m [9].The average permeability is greater than 0.1 mD, and the average porosity is more than 4%.The pore types are diverse, and the proportion of macropores and mesoporous pores is large, while the proportion of micropores is small.Eighty percent of shale gas in the reservoir is stored in inorganic pores in the free state [9].Macrofractures do not develop, and there are a few microfractures and nanoscale organic shrinkage fractures [9].
The dew-point pressure of Fuye 10 HF well is 46.07 MPa, the P-T diagram clearly shows the characteristics of shale condensate gas reservoir, and the gas-oil ratio of blowout is 1243-3142 m 3 /m 3 [9].The "dense-fracture volume stimulation" has a good effect, but the production law of oil and gas in high-pressure shale condensate gas reservoir is still unclear, so it is urgent to develop reasonable and efficient production methods [9].
Condensate gas reservoirs are mainly found in Eagle Ford, Duvernay, Marcellus, and other shale formations.The Duvernay shale in the Simonette block is rich in condensate gas and volatile oil, with an effective thickness of 30-45 m, a high TOC of 2-6%, an effective porosity of 3-6%, and a permeability of 0.0001-0.0003mD.According to the gas-oil ratio, it can be divided into a volatile oil zone, an ultra-high condensate oil zone, and a high condensate oil zone.In the process of development, the ultra-high condensate oil zone and the high condensate oil zone with a high recovery rate should be used first, and then the volatile oil zone should be used to increase the total production [10][11][12].The Eagle Ford shale has strong anisotropy and is weaker when failure along the bedding is possible [13].
Large-scale hydraulic fracturing and horizontal wells are mainly adopted for shale condensate gas reservoirs.The combination of these two technologies has greatly improved the production of shale gas reservoirs.Fracturing has been used in oil and gas production since the 1950s, and horizontal wells have been used in oil production since the early 1980s [14].In the 1980s and 1990s, Mitchell Energy and Development conducted experiments in the Barnett shale in north Central Texas, resulting in commercialization of deep shale gas production and the emergence of large-scale shale gas production [14].The comprehensive application of horizontal wells and multistage fracturing makes the production volume of shale gas reservoir expand from the radial flow state in a limited radius centered on vertical wells to the ellipsoidal flow state centered on horizontal wells, which greatly improves the productivity [15].The production comparison between a vertical well and a horizontal well is shown in Figure 1.The network horizontal well group, which has become the mainstream model in the world, makes the production space of shale condensate gas reservoir more three-dimensional, and the production efficiency is higher and the cost is significantly reduced.
The production process of shale condensate gas reservoirs at home and abroad is mainly divided into four steps: well drilling, fracturing and completion, production, and external centralized transport.The basic process is shown in Figure 2.
Depletion production is often used in condensate gas reservoirs in China, and, thus, in the process of production, gas reservoir pressure drops, condensate oil is precipitated, and the interfacial tension of condensate oil is constantly changing [16].Moreover, the reservoir stress sensitivity of shale condensate gas reservoir is relatively strong.If the pressure is released, the hydraulic fracture is rapidly released during the mining process, the permeability of the reservoir in the fracture area and its adjacent area decreases, and the peripheral gas is difficult to enter the main fracture system.The cumulative production of individual wells will be reduced [17].Therefore, the shale condensate gas reservoir should be controlled pressure production.
--Figure 1.Production comparison between a vertical well and a horizontal well (modified from [15]).
Relevant studies have shown that controlled pressure production can significantly delay the rate of production decline and improve the Estimated Ultimate Recovery (EUR) of a single well [17][18][19][20].In practical engineering, cyclic CO 2 injection is usually used to significantly increase the output of shale condensate gas reservoirs, and at the same time, CO 2 can be stored underground to achieve carbon sequestration [21][22][23][24][25][26].In addition, flue gas injection is also a development method of oil and gas recovery enhancement with higher economic benefits; the main components of flue gas are N 2 and CO 2 , as well as a small amount of SO 2 , CO, and NO x [27,28].Carbon capture, utilization, and storage (CCUS) can help mitigate global warming, reach the goal of carbon peak and carbon neutrality, and promote the low-carbon transformation of socio-economic systems ultimately.At present, CCUS is still at the early stage of technology development in the development of shale gas reservoirs, and relevant technical theories are not mature.Furthermore, the field pilot test of CCUS is rare in shale condensate gas reservoirs.Therefore, there is a broad research prospect in CCUS technology and involved complicated coupled mechanics in shale condensate gas reservoirs.
The time-lapse seismic methods are very sensitive to gas changes and may be used to monitor production in shale condensate reservoirs.Seismic time lapse studies are efficient to monitor the production of gas in any phase, due to the sensitivity of some seismic attributes to fluid contents.The elastic waveform inversion method can be used in time-lapse studies, where the velocity can be estimated from seismic data acquired at different periods [29].Furthermore, seismic wave attenuation, well known by its high sensitivity to fluids (gas and liquid) content, can be used for such a purpose [30].This can be explained by scattering and intrinsic attenuation mechanisms, which are due to heterogeneities and saturation of fluids [31].However, seismic wave attenuation remains difficult to be accurately estimated due to its high sensitivity to seismic noise, and thus it requires a robust methodology.
Shale condensate gas reservoirs are characterized by tightness, low porosity, and low permeability, and must be fractured before commercial exploitation.Fracturing fractures are quite different from natural fractures in terms of morphology and conductivity [32].Therefore, it is very important to establish realistic flow models for a shale condensate gas reservoir to accurately predict productivity.

Research Status of Seepage Law of Shale Condensate Gas Reservoir
Shale condensate gas reservoir has nanoscale pores and fluid migrates mainly in microfractures and fracturing fractures in the reservoir.The migration process of shale condensate gas reservoir has multi-scale complex nonlinear seepage law.The main influencing factors include the multi-scale flow characteristics of shale gas, the adsorption-desorption mechanism, the strong stress sensitivity of the shale matrix and fractures, the abnormal condensate phenomenon around the gas well caused by the drop in reservoir pressure during the development of the shale gas reservoir, the high speed non-Darcy flow in the artificial fractures of shale, the heterogeneity of the reservoir and the multiphase flow.Mastering the seepage law of shale condensate gas reservoir can provide a theoretical basis for the engineering technology development of the shale gas reservoir.

Multiscale Migration Mechanism of Shale Gas
Shale gas can exist in four states: fracture free gas, pore free gas, pore wall adsorbed gas, and kerogen dissolved gas [33].At the nanoscale, according to the intensity of the interaction between gas molecules and pore walls, the transport mechanism can be divided into continuous flow, slip flow, transitional flow, Knudsen diffusion, single-layer adsorbed gas surface diffusion, multi-layer adsorbed gas surface diffusion, and configuration diffusion [34].
Based on the physical adsorption characteristics of shale gas in nanoscale pores, the International Union of Pure and Applied Chemistry (IUPAC) classified pores into macropores (pore diameter > 50 nm), mesoporous pores (2-50 nm), and micropores (<2 nm) [35].The free path of gas molecules is much smaller than the pore diameter of macropores, the collision mainly occurs between gas molecules, and the collision frequency between gas molecules and the pore wall is very low.At this time, it is mainly continuous flow [36,37].With the decrease in pore diameter or gas pressure, the free path of gas molecules will increase.When the free path of gas molecules and pore diameter are of the same scale, the collision frequency between gas molecules and the pore wall increases to a degree that cannot be ignored.In this case, it is slip-off flow, transition flow, or Knudsen diffusion [37,38].When the pore diameter and the gas molecule diameter are at the same scale, only a single gas molecule is allowed to pass through the pore, and the collision frequency between the gas molecule and the pore wall is large, resulting in configuration diffusion [37,39,40].
When gas molecules are adsorbed by the pore wall, surface diffusion occurs.Surface diffusion is a very complex physical phenomenon, which is the interaction between diffused gas molecules and nanoscale pore walls, and it is affected by the interaction force between diffused gas molecules (such as van der Waals force, electrostatic force, etc.), so factors affecting surface diffusion mainly include pressure, temperature, properties of the nanopore wall, properties of the gas molecules, interactions between the gas molecules and the nanopore wall, etc. [40].Surface diffusion, an activation process of adsorbed gas molecules, is a continuous process in which diffused particles jump randomly between adsorption sites, and each jump requires a minimum activation energy, which is usually proportional to the adsorption energy of the adsorbed gas molecules and nanoscale pore walls, and passes through an activation transition state [41].
In nanoscale pores, besides free gas transport, there is also adsorbed gas transport on the pore wall.When the free gas is adsorbed by the wall surface, the adsorbed gas migrates under the action of the adsorption potential field.The gas surface transport quantity is the gas quantity of adsorbed gas passing through the unit cross-sectional area within a unit time, as follows [42]: where C s is the adsorption gas concentration and v s is the migration rate.
In the process of shale gas transmission, when the free path of gas molecules is of the same scale as the nanoscale pore, the collision frequency between the gas molecules and the pore wall is large, the gas is not in thermodynamic equilibrium, and gas slip-off flow occurs.K n [43] was defined in order to determine whether the hydrodynamic equilibrium was established: where λ is the ratio of the gas mean-free-path, that is, the distance between two collisions of gas molecules, m. d n is the characteristic length of nanopore, m.
As is shown in Figure 3, based on the hydrodynamic balance and gas flow, the gas flow region and the gas molecular flow are divided as follows, according to K n [44][45][46][47][48]: (1) When K n < 0.001, gas intermolecular collision is dominant, gas flow meets the continuity condition, and the N-S equation and Fourier's law of thermal conductivity apply.This section is a continuum region where the gas molecular migration is continuous flow, which can be expressed by the Hagen-Poiseuille equation, as follows [15]: where J v is the flow rate of gas continuous flow, mol/(m 2 •s); r is the pore radius, m; µ g is gas viscosity, Pa•s; R is the gas constant, J/(mol•K); T is the formation temperature, K; P is reservoir pressure, Pa; l is the distance of the gas transmission direction, m; ξ mb is the correction coefficient of gas flow in the porous medium.In this case, the gas migration capacity is proportional to the square of the nanoscale pore radius.That is, the gas migration capacity is very sensitive to the pore structure; gas migration capacity is proportional to pressure; since gas viscosity is a function of T 1/2 , gas migration capacity is proportional to T −3/2 , that is, temperature has a significant effect on the gas migration capacity.
(2) When 0.001 < K n < 0.1, the collision between gas molecules and shale pore wall gradually increases, but the collision frequency between gas molecules is still higher than that between gas molecules and pore wall, and gas molecules have velocity on nanopore wall.
In the calculation, the N-S equation and Fourier's law of thermal conductivity are usually used in the vast gas area, but the boundary should consider the velocity slip and temperature jump at the junction of the gas and solid.This section is a slip region where the gas molecular migration is slip flow, which can be expressed as [47]: where J vs is the rate of slip flow, mol/(m 2 •s); b is the gas slippage constant, dimensionless; α is the rarefied effect coefficient and has no dimensionality.With the increase in K n , the slippage effect becomes more obvious and the gas transmission volume increases.
(3) When 0.1 < K n < 10, the collision frequency between gas molecules is equivalent to the collision frequency between gas molecules and the pore wall surface.This section is a transition region where the gas molecular migration is transitional flow.
(4) When K n > 10, the collision between gas molecules and pore wall is very frequent so that the collision between gas molecules can be ignored.This section is a free region where gas molecular migration is Knudsen diffusion, which can be expressed as [33]: where J K is the amount of gas in Knudsen diffusion, mol/(m 2 •s); m is the molar mass of the gas, kg/mol.Knudsen diffusion amount of gas increases as the molecular weight of the gas decreases; Knudsen diffusivity is independent of gas pressure and viscosity; Knudsen diffusion capacity is proportional to T −1/2 , that is, Knudsen diffusion capacity is relatively insensitive to temperature.Free gas in shale nanoscale pores is transported under a pressure gradient or concentration gradient, and the main influencing factors are the pore diameter and gas phase pressure which affects the free path of gas molecules.In a typical shale development process, the pore diameter is 2-1000 nm and the pressure is 1-50 MPa, so K n is around 0.0002-6 [48].K n can be used to determine the transport mechanism of gas molecules in shale nanoscale pores.

Adsorption-Desorption Mechanism
Adsorbed gas and free gas coexist in shale condensate gas reservoir, and the adsorptiondesorption of gas molecules is an important factor affecting its nonlinear seepage mechanism, as shown in Figure 4. Therefore, understanding the adsorption and desorption mechanism of shale condensate gas reservoir is helpful to make better production plans, so as to improve productivity.The desorption and adsorption of shale gas is reversible [49].When the adsorption and desorption of shale condensate gas reservoir are in dynamic equilibrium, the adsorption velocity and desorption velocity are equal.In the process of production, free gas in fractures is constantly extracted, reservoir pressure decreases, and the desorption rate of gas molecules is greater than the adsorption rate, that is, desorption occurs [50].
The specific area of nanoscale pores in shale is large and most of them are oil-wet, which has a strong adsorption effect on gas molecules [39].The adsorption of shale gas can be divided into physical adsorption and chemical adsorption, and physical adsorption is the main one.The reversible physical adsorption is caused by the van der Waals force between the adsorbed gas and the nanoscale pore wall and the adsorbed gas molecules lose one degree of freedom, so the kinetic energy is converted into the heat related to adsorption and the heat consumed by physical adsorption is less.Chemisorption is irreversible and is adsorbed by ionic bonds, which requires a large energy to open ionic bonds during desorption [15].
The main factors affecting the adsorption capacity of shale are organic matter content, reservoir temperature, formation pressure, shale water content, gas reservoir composition, etc.
Lu and Li [51] discovered that the adsorption capacity of shale samples was approximately linearly related to the total organic matter content, that is, gas dissolution in organic matter was the main mechanism of adsorption gas storage; the adsorbed gas volume decreased with increasing temperature.Yu and Sepchroori [52] showed that the adsorption capacity of CH 4 in Barnett shale increases with the increase in pressure.Yuan et al. [53] found that in shale with high water content, the adsorption capacity and diffusion coefficient of gas decreased significantly, that is, the gas storage capacity and gas transport rate of shale were reduced by the presence of water.
Typical gas adsorption models, namely the relationship between adsorption capacity and pressure under isothermal conditions, mainly include the following models [54,55]: (1) Henry adsorption theory [56,57] The Henry adsorption formula is as follows: where M is the mass of adsorbed gas on solid per unit mass, p is the equilibrium pressure of gas, and k is a constant.When the pressure is low, the isotherm is close to a straight line, which is in line with Henry adsorption theory.
(2) Freundlich adsorption theory [58] The Freundlich adsorption formula is as follows: where n is a constant, which depends on the adsorbent, the type of adsorbent, and the adsorption temperature.
(3) Langmuir adsorption theory [59] The Langmuir monolayer adsorption formula: The above formula is more suitable for single-molecule adsorption [60].As is shown in Figure 5, Molecular layer adsorption model can be divided into single-molecule adsorption and poly molecule adsorption.Through experiments, Lu and Li [61] proved that Langmuir model could accurately characterize adsorption at a single temperature, but was not applicable to multiple temperatures, while the bi-Langmuir model could accurately characterize the influence of pressure and temperature.Ambrose et al. [62,63] calculated the volume of pore space occupied by free gas and the adsorbed layer based on the Langmuir adsorption model, improved the traditional method, proposed a new algorithm to calculate shale gas geological reserves, and studied the adsorption characteristics of methane in organic matter.The results showed that it was very important to consider the adsorption of multiple components in shale gas reservoir calculation.Multi-component adsorption with the extended Langmuir model was simulated by Yang et al. and the role of multi-component adsorption in the development of Eagle Ford shale condensate reservoir was studied [64].The adsorption law of pores with different shapes still conformed to the Langmuir single-layer adsorption law was proved by Song et al. [65].

= +
(4) BET adsorption theory: Brunauer et al. [66] studied the multi-molecular layer adsorption theory, namely the BET adsorption theory.The multi-molecular layer adsorption formula is as follows: Generally, when p/p 0 = 0.05-0.35, the above formula is valid.Yu et al. [67] conducted methane adsorption experiments on Marcellus shale core samples, and found for the first time that gas adsorption in shale gas reservoirs was characterized by multilayer adsorption, and the experimental results deviated from Langmuir isotherm and matched the BET isotherm.By comparing historical data, it is found that gas adsorption obeying the BET isotherm contributes more to the total recovery than gas adsorption obeying the Langmuir isotherm.
Related studies found that adsorption significantly increased the volume of condensate [68], resulting in fluid accumulation and blockage in the reservoir, and reduction in production output.Therefore, the technology of determining desorption pressure and re-injecting gas into the reservoir to slow down the desorption rate, that is, recycling gas injection to improve oil recovery, has broad research prospects.

Strong Stress Sensitivity Characteristics of Shale Matrix and Fracture
The matrix pores of shale gas reservoirs are at the nanometer scale.There are a certain number of micron-sized pores and a large number of microfractures in shale reservoirs.In the process of shale gas production, with the continuous exploitation of shale gas, the pore pressure decreases, and the natural fractures may gradually close, which leads to the decrease in reservoir permeability.It shows strong stress sensitivity characteristics.Not only is there stress sensitivity in the shale matrix, the stress sensitivity in the stimulated reservoir volume (SRV) area is more serious.Stress sensitivity is an important characteristic of shale reservoir, which has a great influence on permeability.The diagram of strong stress sensitivity characteristics of fractures is shown in Figure 6.At present, some scholars have conducted a series of studies on the strong stress sensitivity of reservoirs.From 2019 to 2022, based on core displacement experiments of fractured shale [69], pressure pulse attenuation experiments [69,70], and triaxial core compression experiments [71], the damage mechanism of stress sensitivity of the reservoir was explored by Xue et al.In addition, in 2019, the effect of stress sensitivity on material permeability under different pressures and gas compositions was explored by Singh et al. [72,73].From 2019 to 2022, the influence of stress sensitivity coefficients of matrix [74,75], microfractures [74], and hydraulic fractures [74,75] on gas production was analyzed by Liu et al.Based on the above research, stress sensitivity is an important characteristic of shale reservoir, which has a great influence on the permeability field.

Special Phase Transition Mechanism of Shale Gas Anomalous Condensation
Compared with conventional reservoirs, the fluid flow in shale is not only complex, but also affected by various transport mechanisms, such as confining phase behavior.Nanopores in shale matrix may lead to significant changes in capillary pressure and critical properties (critical pressure and critical temperature).The phase equilibrium behavior in shale reservoirs may also differ significantly from the bulk fluid phase equilibrium behavior in PVT units [74,75].
Many scholars have begun to study the mechanism of phase transition in nanopores.In 2011, Sapmanee [76] et al. found that when the pores reached the micrometer and nanometer scales, the pore confinement effect would change the fluid properties.Due to the pore confinement effect, the space limitation in the tiny pores would affect the phase behavior, fluid properties and flow behavior of the fluid.In 2012, Du et al. [77] obtained the effect of capillarity on fluid phase balance by studying PVT properties of confined fluid in nanopores.In 2017, Bui and Akkutlu [78] et al. found that the fluid in nanopores changes significantly with the pore size, and the change of capillary pressure and critical properties (critical pressure and critical temperature) in nanopores may have a significant impact on the fluid phase behavior.In 2022, Tian [79] et al. established a tight reservoir relative permeability model in order to study the effect of nano-confinement on relative permeability and reservoir flow, and the results showed that with the decrease of pore size, the effect of nano-confinement on relative permeability became increasingly significant.
Due to the existence of the pore confinement effect, the phase equilibrium of the fluid will also be different.Different from the conventional phase state calculation process, many scholars have proposed phase state calculation methods considering the pore confinement effect.Currently, two parts of the methodology exist to account for nanoconfinement effects, namely, critical property shift and consideration of capillary pressure.
(1) Critical property shift The critical temperature displacement equation is based on the state equation [80]: where T c and P c are the bulk critical temperature and pressure, T cp and P cp are the pore critical temperature and pressure, respectively, the critical pressure and temperature are corrected by the above equations.
(2) Capillary pressure Due to the large capillary pressure in nanopores and the phase pressure difference when oil and gas phases coexist, capillary pressure should be considered for calculation.
The well-known Yang Laplace equation is usually used for capillary pressure in porous media [81]: where σ is interfacial tension, which is the surface tension between solid and fluid; θ is the contact angle; and r is the radius of the capillary in the porous medium.The interfacial tension [82] can be calculated by the following equation: where ρ v and ρ L are the densities of the gas and liquid phases, respectively, and P ch,i are the parachute model parameters of the liquid or gas phases, which can be obtained by empirical formulas.x i and y i are the molar fractions of gas phase and liquid phase components, which can be obtained by the phase equilibrium calculation method proposed by Yuan Zhang [83] et al.The specific steps are shown in Figure 7. Through the above calculation process, the phase state diagram can be drawn for different pore diameters, and the phase state diagram considering the pore confinement effect is roughly shown in Figure 8 [84].Many scholars have undertaken research on phase shift and pore confinement and achieved certain results.In 2011, a simple, fast, and reliable alternative to a similar method in flash calculations was proposed by Yinghui [85], which realized the rapid convergence of phase calculation.In 2017, Brian [86] et al. found that the smaller the pore size, the higher the ratio of the specific surface area of the pore to the volume of the fluid, and the smaller pore can inhibit the gas bubble point.Based on the above characteristics, Brian carried out simulations through the relative permeability model to explain the capillary pressure and the interaction between the fluid and the matrix.In 2017, Arash [87] et al. carried out thermodynamic modeling of fluid reservoirs with different PVT properties considering the effect of pore size, and proposed detailed steps of finding the bubble point line and dew point line.The results showed that the established thermodynamic model could better simulate the PVT properties of shale condensates.In 2018, Liu [74] et al. established a mixed model of condensate gas reservoir by considering the influence of various transport mechanisms, such as the behavior of the pore-confined phase in the shale matrix.They simulated the phase transition under nano-confinement by using the quasi-critical properties in nanopores and introduced the specific calculation process.In 2019, An improved artificial neural network model was proposed by Wang [88] to achieve high-precision phase stability prediction, which not only reduces the computational cost, but also naturally and accurately predicts the saturation pressure, so as to accurately describe the phase behavior of the fluid.In 2020, Bai [89] et al. used the component model considering the limiting phase behavior to simulate the multiphase flow in the reservoir, introduced the capillary pressure into the workflow of the phase stability test and flash calculation, and evaluated the influence of the limiting phase behavior on the production dynamics.The results showed that the limiting phase behavior effect inhibited the gas bubble point pressure, thereby improving the production.In 2021, Chen [90]  The neural network was used for simulation and prediction, which had good accuracy and prediction effect, and could also significantly capture the influence of capillary action.It is a powerful means to predict phase behavior.In 2022, a multi-component hydrocarbon transport model based on a pore network was proposed by Song [92], and elaborated the thermodynamic mechanism of multi-component hydrocarbon transport.On the pore network model, considering the change in the oil-gas capillary pressure with oil saturation, a coupling model of thermodynamic phase equilibrium and fluid occurrence was further established to characterize the phase behavior.In 2022, Kong [84] et al. conducted a large number of studies on the behavior of fluid phase constrained by nanopores in tight shale reservoirs, optimized the equation of state (PR EOS), calculated the behavior of phase constrained by nanopores by changing the key properties of different components, and proposed a characterization method for the behavior of fluid phase in nanopores.In 2022, a calculation process of micro-nano confined space phase was proposed by Wei Bing [93], which took into account the changes in capillary pressure, critical temperature, and critical pressure.Based on the improved EOS equation, the convergence rate of multi-component mixtures was greatly accelerated, and the traditional calculation method was simplified.In 2022, Du [94] et al. considered the influence of molecular adsorption and confinement effects in nanoscale shale pores on fluid phase behavior, and proposed an equation of state considering adsorption to study fluid properties and phase behavior changes during gas injection in condensate gas reservoirs.In 2022, Zhang [95] et al. established a comprehensive thermodynamic model to simulate the shale condensates gas reservoir, which is suitable for bulk fluids and confined fluids.Through experimental simulation, it was shown that capillary interaction and adsorption will lead to an obvious deviation of the phase envelope.In addition, pore geometry also has a significant impact on phase behavior.This model provides useful insights into the phase behavior of nanopores.
In recent years, some scholars have found that CO 2 injection will also have an impact on the behavior of the fluid phase.In 2021, a new approach proposed by Zheng [96] is to study the phase behavior of hydrocarbon systems considering nanopore confinement and modify the phase equilibrium calculation.The results show that the decrease in nanopore size leads to a significant difference in the phase envelope, and the confinement of nanopore reduces the bubble point pressure.With the decrease in pore size, the interfacial tension decreases, and the capillary pressure increases significantly.Zhang [97] et al. also found the effect of CO 2 injection on the fluid phase balance in nanopore confinement.The results showed that the nanopore confinement effect could not be ignored, and the smaller the pore size, the greater the nanopore confinement effect, and the confinement helped to suppress the bubble point pressure.In 2022, Ruan Hongjiang [98] analyzed the phase characteristics of condensate gas based on the gas-liquid two-phase equilibrium theory through the phase state experiment of CO 2 injection.The experimental results showed that with the increase of the proportion of CO 2 injection, the phase state of condensate gas reservoir would shift to different degrees, and the parameter yield of condensate gas reservoir could also be effectively improved.According to the research of the above scholars, with the injection of CO 2 , the phase behavior will shift to different degrees, and unlike the pore limitation mentioned above, the injection of CO 2 helps to suppress the bubble point pressure, as shown in Figure 9 [97].CO 2 injection can not only change the phase balance and have a certain impact on the phase behavior of the fluid, but also improve the efficiency of oil and gas flow, evaporation, and condensation flooding.This makes the technology of CO 2 injection have great potential to store CO 2 to cope with global warming and effectively reduce greenhouse gas emissions, which will have great potential in future oil and gas extraction.

High-Speed Non-Darcy Flow in Artificial Fractures of Shale Reservoirs
However, in fact, high velocity flow cannot be accurately described by Darcy's law.When the flow reaches a high velocity or a high Reynolds number, the inertial effect will also become significant.Sometimes an inertial term is added to the Darcy equation, called the Forchheimer term, to explain the nonlinear behavior of pressure difference and velocity data [99]: where β is non-Darcy flow coefficient, and the parameters depend on reservoir characteristics, including permeability, porosity, channel complexity, or curvature; P is reservoir pressure, MPa; ρ is the density of flowing fluid, kg/m 3 ; k is the permeability of porous media, mD; µ is gas viscosity, Pa•s; v is the gas flow rate, m•s −1 .The microfractures of shale reservoir show Darcy seepage characteristics.In the artificial fracture area of shale, due to the multi-scale pore medium structure and multireservoir characteristics of shale reservoir, the shale migration mechanism is very complex, and the flow state of fluid may have high-speed non-Darcy flow.The Darcy equation that only considers viscous flow cannot be described.In order to determine whether the flow of shale gas in artificial fractures belongs to high-speed non-Darcy flow, the Reynolds number obtained by the following formula is usually used as the basis [15]: where Re is the Reynolds number; ϕ is the porosity of porous media.Through experiments, Zhu et al. [99] found that the flow in artificial fractures outside the critical Reynolds number Re of 0.2-0.3belongs to high-speed non-Darcy flow.Some experimental [100] or theoretical [101][102][103] research on high-speed non-Darcy flow has been carried out.In 2018, using a new proposed experimental measurement method of permeability and inertia coefficient of porous media under high-speed non-Darcy flow, the parameters of high-speed non-Darcy flow of polyethylene porous materials were tested by Zhang et al. [100].From 2019 to 2022, transport phenomena such as high-speed non-Darcy flow during shale gas production were characterized by Chen et al. [101][102][103].
The above research provides a certain reference for the study of high-speed non-Darcy flow of shale gas reservoirs.

Reservoir Heterogeneity
Compared with conventional gas reservoirs, shale gas reservoirs have completely different reservoir geological, physical, and mechanical properties.The pore structure of shale reservoirs is complex and multi-scale, which is mainly composed of four different types of porous media.Due to the variety of shale gas storage methods and different organic matter content, there are great differences in reservoir properties.
Shale condensate gas reservoirs have the characteristics of a small porosity and low permeability.Reservoir heterogeneity is strong, and natural fractures are generally developed.Due to the differences in the formation environment of different types of organic-rich shale reservoirs, shale reservoirs have multiple vertical heterogeneity, and the lateral permeability of shale reservoirs is also quite different from the vertical permeability.Due to the relatively complex matrix composition of shale reservoirs, the matrix itself has certain heterogeneity.Because bedding and fractures control the shape and spatial distribution of fracturing fractures, with the increase in depth, shale reservoirs have strong heterogeneity and anisotropy in the horizontal direction, which results in the vertical heterogeneity of shale oil and gas reservoirs being significantly stronger than the horizontal heterogeneity.Parameters such as organic matter abundance, porosity, permeability, reservoir fracability, in situ stress, gas content, and absorption ratio change greatly in the vertical direction.
After hydraulic fracturing, the SRV area around main fracture will produce irregular complex network fractures, and the density of microfractures gradually decreases to the surrounding area, which will lead to strong heterogeneity of permeability and porosity in the SRV area around main fracture.The fractal characteristics of the reservoir in the SRV area are more obvious.The diagram of reservoir heterogeneity at the SRV area is shown in Figure 10.The expression of non-uniform porosity at the SRV area considering stress sensitivity is as follows [104,105]: The expression of non-uniform permeability at the SRV area considering stress sensitivity is as follows [104,106]: where y is the distance to the main fracture, m; ϕ wi is the initial porosity at y = w F /2; D is the fractal dimension; w F is the width of main fracture, m; k wi is the initial permeability at y = w F /2, mD; θ 1 is the fractal index; p i is the initial pressure, MPa; C k is the permeability compression coefficient, MPa −1 ; C ϕ is the porosity compression coefficient, MPa −1 .At present, a series of experimental [107,108] and theoretical [107,[109][110][111][112] studies on reservoir heterogeneity characteristics have been carried out.In 2021, using high-precision high-pressure mercury intrusion (HPMI) experimental technology, the micro/nanopore structure characteristics of shale oil were studied, and based on fractal theory, its heterogeneity was characterized by Xia et al. [107].In 2022, using computed tomography (CT), the pore distribution of shale at a full scale was characterized by Zhan et al., and the characteristics showed horizontal uniformity and vertical heterogeneity at the macroscale [108].In addition, in the theoretical research on reservoir heterogeneity characteristics, from 2019 to 2022, based on fractal theory [109,110], the Gaussian distribution function [109], and the lognormal probability density function [111], the distribution of reservoir permeability and seepage law were characterized by Qi et al.In 2019, the expression of the apparent permeability of nanopore gas was also deduced by Liu et al. [15].The above research provides some guidance for the characterization of reservoir heterogeneity.

Multiphase Flow
During the production of shale condensate gas reservoirs, the formation pressure keeps dropping.When the bottom hole pressure is lower than the dew point pressure, the condensate will begin to precipitate, and when the local layer pressure is lower than the critical flow pressure, the precipitated condensate will form a continuous phase and begin to flow near the wellbore, and eventually form the seepage state of oil, gas, and water two-phase flow, or even the three-phase flow of oil, gas, and water under water conditions.
Relative permeability data are important basic data in the process of oil and gas reservoir development, which can be widely used in many aspects, such as numerical simulation of oil and gas reservoirs, dynamic analysis of production data, and EUR, so as to evaluate and reasonably adjust production measures and planning [112,113].The inaccuracy of relative permeability estimation will adversely affect reservoir resource evaluation and development [114].
During reservoir engineering development, oil, water, and gas phases may coexist, and three-phase flow usually occurs when water saturation is higher than residual water saturation and oil and gas are mobile phases [113].Three-phase relative permeability is required to calculate reservoir dynamic data, under production conditions such as carbon dioxide, burning reservoir, water-gas flooding, steam flooding, micellar injection, nitrogen injection, etc. [113].
The methods of estimating and measuring relative permeability curve directly in laboratory are mainly divided into the steady-state method and the unsteady-state method.The steady-state method is used to measure relative permeability by injecting fluid into porous media with a constant oil-water ratio until saturation and pressure are in equilibrium, and then the saturation, pressure, and flow under this equilibrium state are obtained, and the relative permeability is calculated by Darcy formula [115].Each measurement point requires the fluid saturation distribution in the rock to reach a stable state, which takes a long time and is not applicable to dense shale [116].The measurement of relative permeability by the unsteady-state method assumes that the distribution of oil and water saturation in rock samples is a function of time and distance in the process of water flooding, so the relative permeability of two phases can be calculated by measuring the flow rate (or pressure drop) in the process of water flooding with a constant pressure drop (or constant velocity) [115].It is difficult to measure three-phase flow experimentally.Therefore, the relative permeability curves of oil-water and oil-gas are usually obtained experimentally at home and abroad, and then calculated using the mathematical model proposed by Stone [117,118] in 1970 and 1973 and modified by later generations.As far as we know, it is difficult to directly measure the relative permeability of shale condensate gas reservoir with the existing instruments and equipment, so it is urgent to develop more suitable experimental equipment for shale condensate gas reservoir.Shale condensate gas reservoirs have a small porosity, a low permeability, and a large capillary force [119], which will have a certain impact on the relative permeability.In 2004, the multiphase permeability under different water saturation and net overburden stress was studied by Shanley et al. and capillary pressure and net overburden stress were analyzed to better understand the relationship between rock structure and gas production [120].In 2017, the obtained capillary pressure curve was integrated by Brian Stimpson and Maria Barrufet to obtain the relative permeability curve [86].In 2022, Su et al. [121] improved the traditional JBN method and established a processing method for experimental data of oil-water relative permeability.The results showed that with the increase in capillary force, the relative permeability of the oil phase increased, but the relative permeability of the water phase did not change significantly.In 2022, considering the change in oil and gas capillary pressure with oil saturation, a multi-component hydrocarbon transport model and a coupled thermodynamic phase equilibrium fluid occurrence model in the pore network were established by Song et al. [91] to study the thermodynamic mechanism of multi-component hydrocarbon transport.The results showed that the initial pressure gradient was affected by yield stress and pore size.
Relative permeability is affected by reservoir pressure performance, water saturation, pore size, pore distribution, and pore connectivity.In 2011, Orangi et al. [122] found that reservoir fracturing performance was very sensitive to fracture permeability and matrix relative permeability.In 2017, Ojha et al. [115] found that the relative permeability of water phase decreased sharply with the decrease in water saturation, while the relative permeability of hydrocarbons decreased at a relatively slow and constant rate with the increase in water saturation.However, this method did not consider the effects of relative permeability hysteresis and mixed wettability.In 2020, Li et al. [123] found that the relative permeability curve was mainly affected by pore diameter distribution.The smaller the pore radius was, the stronger the gas migration capacity was, and the larger the relative permeability of the gas phase was.If pore connectivity was known, the flow capacity of shale gas under different humidity and the flow behavior of shale gas in pores could be obtained.In 2022, Tian et al. [79] discovered that the influence of the abnormal viscosity effect on the relative permeability of water and oil phases was related to the average pore radius.With the increase in the slip length of the oil phase, the relative permeability of the water phase decreased, as did the relative permeability of the oil phase.With a more significant contact angle effect, the relative permeability of the water phase increased, but the relative permeability of the oil phase had no obvious change.However, the effect of pore connectivity and fluid distribution was ignored.
It is difficult to measure the phase behavior of multi-hydrocarbon fluids in micro/nanoscale channels by conventional experimental methods.Advanced microfluidic experimental technology can be used to study the characteristics of fluid phase in the nanoscale limited space, which can provide certain data support for unconventional shale gas development engineering technology.In 2021, Song et al. [124] designed a microfluidic experiment and fed the processed experimental data into a deep neural network (DNN) model, which was a deep learning algorithm.The results showed that phase saturation, wetting sequence and pore topology were the key factors affecting the relative permeability of the two phases.
However, as far as we know, with insufficient instruments and equipment, microfluidic technology only stays in the phase visualization stage, because the nanochannel is too narrow to measure some flow parameters, that is, to simulate the dynamic process of fluid flow in the nanopores.In addition, the analysis chip is mostly disposable, and the experimental cost is generally high.At present, the application of microfluidic technology in the phase state research of shale condensate gas reservoirs is just being developed, and a large room for development still exists in the future for the production of more accurate microfluidic chips, equipment upgrading to measure fluid dynamic flow characteristics, and other aspects.

Theory, Methods, and Technology of Shale Condensate Gas Production Well Dynamic Data Analysis
The production well dynamic data analysis of shale condensate gas reservoir is that according to the production dynamic characteristics of production wells, reservoir engineering methods are used comprehensively to analyze the physical properties of reservoirs and the production decline law of production wells, so as to accurately predict the production capacity of production wells and provide theoretical basis and technical support for rational production allocation and efficient development [125,126].At present, shale condensate gas production well dynamic data analysis methods are mainly divided into the empirical method, the characteristic curve analysis method, and the artificial intelligence method.

Empirical Method
The common empirical methods for the analysis of production well dynamic data of shale condensate gas reservoir include the Arps decline model, PLE decline model, Duong decline model, GEPD decline model, SEPD decline model, LGM decline model, etc.This kind of method is simple to calculate, but the characteristics of reservoir seepage in shale condensate gas development are not fully considered.It is greatly affected by geography, production mode and other factors, and its reliability is poor.
(1) Arps decline model Arps decline model does not require the parameters of the reservoir or well test, and can predict and calculate production based on the empirical relationship between production and time.This method is suitable for the analysis of production decline in shale gas condensate wells with constant pressure production and flow reaching the boundary control flow [127][128][129][130].The Arps decline model can be further divided into exponential decline, hyperbolic decline and harmonic decline.Specific models are as follows: where q i : initial daily output, 10 4 m 3 /d; q t : daily output at time t, 10 4 m 3 /d; D i : initial decline rate, d −1 ; t: decline time, d; a: decreasing index, dimensionless.When a = 0, it is exponential decline; when 0 < a < 1, it is hyperbolic decline; when a = 1, it is harmonic decline.The cumulative output is where Q t is the cumulative output at time t, 10 4 m 3 /d.The hyperbolic decline model is the most widely used in shale gas production decline prediction.In practical application, the decline index a often exceeds 1, so it is usually called the generalized hyperbolic decline model.
(2) PLE decline model Ilk et al. [131] modified the traditional exponential decline model and proposed the PLE decline model considering that the decline rate of shale gas well production would decrease in the later period: where a is the time index, dimensionless; D ∞ is the decline rate at infinite time, d −1 ; D 1 is the first-day decline rate, d −1 .This method is suitable for shale gas condensate production wells with a constant pressure production (or close to constant pressure production), and the flow is in the linear flow, bilinear flow, or transition flow before boundary control flow.
(3) Duong decline model Based on the characteristics of linear fracture flow in shale condensate gas production wells during long-term production, Duong [132] proposed a new decline model for production decline after fracturing of tight gas wells: After separating variables and integrating, the cumulative output can be written as where q 1 is the first-day tight gas well production; d is the fracture flow index, for linear flow, d = 1/2, for bilinear flow, d = 1/4.This method is suitable for shale gas condensate production wells with constant pressure production (or close to constant pressure production), and the flow is in the linear flow or bilinear flow.The Duong decline model is more suitable for the early production decline of production wells [133] and it is easily affected by the fluctuation of production data and prone to multiple solutions [134].
(3) GEPD decline model In 2021, the GEPD decline model proposed by Chen is [135] q = q i exp(−t δ /c) (24) where δ is the universal exponent, 0 < δ < 1; when t = 0, q = q i , and q is the initial theoretical production of the well.The well controlled production can be derived as follows: In 2021, Chen [135] used two production decline models to evaluate shale gas well production, and found that the GEPD decline model was closer to the reality than the hyperbolic decline model.
Most typical production decline models of shale condensate gas reservoirs are fitted to production data, and different parameters have different effects on production data, so parameter determination is particularly critical.In the actual development process, the flow of shale condensate gas reservoir has different stages, but the empirical method only uses a single function to predict the output of the whole production process, which has large errors and lacks rigor.

Characteristic Curve Analysis Method
The characteristic curve analysis method is to use the known production data to carry out inversion and fitting, so as to predict the productivity more accurately.However, the characteristic curve analysis method is related to the production mode and stage of gas reservoir, which has higher requirements for the selection of seepage model and should consider complex effects.Therefore, the selection of the seepage mathematical model is particularly important.
In 2017, considering the characteristics of fluid phase transition and fracture network during development, using a trilinear flow model to characterize SRV and a Langmuir isothermal adsorption model to characterize the adsorption-desorption of gas, a production well productivity prediction model for shale condensate gas reservoirs was established by Wu et al. [136].In 2019, considering the half-length of hydraulic fractures, the number of fractures, fracture spacing, and the adsorption-desorption mechanism of gas molecules, the semi-analytical solution of the production decline curve under variable pressure was obtained by Lu et al. [137].In 2020, Zhang and Meybodi [138] clarified the fracture properties of a multi-fractured horizontal well (MFHW) by iteratively calculating fracture volume, permeability, fracture compressibility, and permeability modulus, and hydraulic fracture dynamics was determined by analyzing flowback data and long-term production data.In 2020, based on the pressure transient analysis (PTA) and rate transient analysis (RTA) of production data of shale gas condensate reservoirs, a semi-analytical mathematical model for multi-stage hydraulically fractured wells (MHFW) of shale condensate gas reservoirs was established by Dahim et al. [139].According to the rich gas production of MHFW, the type curves of RTA and PTA were drawn.The effects of the horizontal well length, controls of linear and elliptical flows, the number of multi-fractures, and other factors on pseudo pressure and yield data were studied.In 2021, Luo et al. [140] used the source function method to discretized fractures, superimposed the hydraulic fracture model, and coupled the two models to obtain the unsteady seepage flow and productivity model of MFHW in shale gas reservoirs.They then analyzed the influences of fracture conductivity, fracture half-length, fracture spacing, skin factor, reservoir ratio, and leakage coefficient on productivity.In 2021, a new RTA was developed by Ren and Lau [141] to analyze shale gas production data by establishing trilinear flow models and simplifying equations to characterize bilinear and linear flow states.In 2021, considering the distribution of kerogen and the complex migration mechanism of gas, a production analysis model for shale gas reservoirs was established by Zeng et al. [142] to evaluate production.In 2023, a production prediction method for shale condensate gas reservoir based on the convolution equation was proposed by Wang et al. [143], which was applied to predict the dynamic production of horizontal wells in Duvernay shale condensate gas reservoir.Compared with the production decline analysis method in HIS Harmony RTA, this method has a better fitting effect and a higher accuracy in predicting EUR.
The analysis of dynamic production data of oil and gas reservoirs by the characteristic curve analysis method is usually based on some linear assumptions, but the production principle of shale condensate gas reservoirs is characterized by complexity and nonlinearity.Therefore, nonlinear problems should be considered more carefully when applying this method to dynamic data analysis, and the linearization solution of nonlinear problems should be further studied.Therefore, the applicability of this method in the field of production dynamic data analysis of shale condensate gas reservoir needs to be further improved.

Artificial Intelligence Method
In practical engineering, artificial intelligence technology plays an important role in analyzing complex dynamic production data and accurately predicting output.
In 2017, according to the characteristics of reservoir reconstruction data, an adaptive threshold denoising production prediction model was proposed based on a BP neural network nonlinear prediction algorithm by Zhu et al. [144].In 2017, a principal component analysis method was proposed to analyze Eagle Ford production data to predict production by Khanal et al. [145].In 2021, a machine learning model was established to predict the five-year cumulative production curve of MHFW by Chaikine and Gates [146].In 2021, a machine learning method was used to analyze the production data of Duvernay shale condensate gas reservoir and Kong et al. [147] found that the important factors affecting cumulative production in the first 12 months were percentage, condensate/gas ratio (CGR), completion length, fracturing fluid volume/length, and reservoir pressure.Furthermore, CGR is also the most important factor affecting long-term production.In 2021, based on the deep learning algorithm, the EUR evaluation method for shale gas wells was designed by Liu et al. [148].In 2022, based on the machine learning model, a production prediction model was established by Bhattacharyya and Vyas [149].The Eagle Ford production data was analyzed so that the production decline curve and EUR were predicted.Comparing the results with the actual data, it was found that the machine learning model could simulate the reservoir well.In 2022, based on machine learning methods, Hui et al. [150] estimated that the main factors affecting unconventional shale gas production were production metrics, formation pressure, effective porosity, total organic carbon, gas saturation, and shale thickness.
The production data analysis method based on artificial intelligence can process complex data with large amount of information and quantify the influence of multiple factors on productivity.In application, a certain amount of basic data is required as learning samples.However, the existent production dynamic data of shale condensate gas reservoirs are relatively limited in comparison with conventional reservoirs, so physical constraints for the development of shale condensate gas reservoirs should be added to improve the accuracy of prediction results when applying the artificial intelligence method to productivity prediction.

Theory, Method, and Technology of Numerical Simulation of the Shale Condensate Gas Reservoir
Numerical simulation of shale condensate gas reservoir is where an established model or commercial software is used to simulate the complex phase change and flow state of multi-components in the production process of shale condensate gas reservoir, to understand the characteristics of fluid flow and the main control factors affecting production well output, so as to provide a theoretical basis and technical scheme for production practice.
When modeling is used for simulation, the fluid-structure coupling model is usually established based on the characteristics of reservoir containing multi-scale pores and fractures and the characteristics of fluid multi-component and multi-phase state.Typical fluid-structure coupling models mainly include the equivalent continuum model, dual medium fluid-structure coupling model, discrete fracture model, etc. [151].The equivalent continuum model is to characterize fractures based on equivalent mechanics, treat the reservoir as a continuum, and then analyze it further based on continuum theory.This model is relatively simple and has large errors; the material exchange between fractures and matrix rocks is considered, but the anisotropy of matrix is ignored, so it is more suitable for the spatial scale across small fractured matrix; the development of discrete fracture model is relatively mature [151].The embedded discrete fracture model is a branch of the discrete fracture model, which can simulate the complex fracture matrix and does not require mesh encryption near the fracture as the discrete fracture model does, thus running faster [152].In addition, some scholars [153] proposed to establish a hybrid model, that is, to establish a single porous medium model, a dual medium model, and an embedded discrete fracture model, respectively, to simulate the unmodified fracturing area, the reformed fracturing area, and the hydraulic fracture, and to establish a seepage model to simulate fluid flow.
Many scholars have carried out a series of research on the numerical simulation of shale condensate gas reservoir.In 2014, Haghshenas et al. [68] established a heterogeneous reservoir model and used commercial software CMG to study the effects of heavy hydrocarbon component desorption, pore structure and connectivity, fluid composition, and flowing bottom pressure on production.In 2015, the method of numerical analysis was used by Wan et al. to analyze experimental data of cyclic nitrogen injection, and the results showed that molecular diffusion played an important role in the process of oil extraction, and cyclic injection technology could significantly improve production [154].In 2016, considering multicomponent apparent permeability, adsorption, and molecular diffusion, the shale gas condensate reservoir model was established, an equation of the state model was used to simulate the complex phase changes of the fluid, and the effects of capillary pressure and multi-component mechanism on the production were studied by Jiang and Younis [4].It was also found that CO 2 huff-n-puff could significantly improve the hydrocarbon recovery.In 2018, Wan and Mu [3] numerically simulated the effect of cyclic CO 2 injection under nanopore effect on production well recovery in an Eagle Ford shale condensate gas reservoir.In 2018, Pankaj et al. [155] used a complex hydraulic fracture model to simulate the whole wellbore and found that cyclic CO 2 injection could significantly improve the oil recovery by numerical simulation.In 2019, Liu et al. [74] developed a hybrid model of the composition and migration mechanism of fluid flow in condensate gas reservoirs, and microfractures and hydraulic fractures are conducted by the multiple interacting continua model and the embedded discretized fracture model.Based on this model, the effect of gas migration mechanism on the production of multistage fractured horizontal wells was analyzed.The results showed that Knudsen diffusion mainly occurred when simulated reservoir volume was absent and matrix permeability was very low.The adsorption effect was mainly related to the reservoir pressure level.In 2021, Jing et al. [156] established a reservoir model of MFHW, used an equation of state model to simulate the flow state of Duvernay shale condensate, and quantified the effects of the nanoconfinement effect, adsorption, and molecular diffusion on production.It was found that recovery could be improved by adjusting production-injection ratios.In 2022, an RTA-assisted numerical history-matching technique, which could improve the accuracy of numerical simulation of shale condensate gas reservoirs, was proposed by Hamdi et al. [157].The calibrated model could be used to optimize oil production technology.
When numerical simulation is used to simulate the oil-gas-water phase characteristics of shale condensate gas reservoir, much more formation and fluid parameters of shale condensate gas reservoir (such as relative permeability, multi-scale flow characteristics, heterogeneity of shale, etc.) are involved and are also difficult to estimate.As a result, there will be great uncertainties in the history fitting process for the numerical simulation.It is suggested that the experimental ability to test the formation and fluid parameters of shale condensate gas reservoir should be emphasized in future research, which can largely help to more accurately describe the characteristics of complicated flow progress and provide strong support for numerical simulation technology.

Summary and Prospect
In this paper, the reservoir characteristics, resource distribution, and development mode of shale condensate gas reservoir are summarized, as well as the research progress of seepage law and development technology of shale condensate gas reservoir at home and abroad in recent years.The following conclusions are drawn: (1) The permeability and porosity of shale condensate gas reservoir are very low, and the recoverable reserves are very rich, which has a considerable development prospect.
The combination of large-scale hydraulic fracturing and horizontal wells, and controlled pressure production through the injection of gas (such as CO 2 , N 2 , flue gas, etc.) during production, have greatly increased the production of shale condensate gas reservoirs while achieving carbon sequestration.Theories related to CCUS are still at the early stage of development in the development of shale gas reservoirs.(2) Compared with common reservoirs, the seepage law of shale condensate gas reservoir is more complicated: the migration mechanism of shale gas is related to pore diameter, temperature, pressure, gas molecular properties, and other factors.There is adsorption and desorption in the process of migration of shale condensate gas reservoir; during the production process, pore pressure decreases, natural fractures gradually close, and reservoir permeability decreases, showing strong stress sensitivity; the decrease in reservoir pressure leads to abnormal condensate around the gas well during development; in the artificial fracture area of shale, the migration mechanism of shale is very complicated due to the characteristics of multi-scale pore medium structure and multi-reservoir mode, and the fluid flow state may be high-speed non-Darcy flow; shale reservoir has complex heterogeneous characteristics; in the production process of shale condensate gas reservoir, when the formation pressure is lower than the dew point pressure, the condensate will start to precipitate, and when the formation pressure is lower than the critical flow pressure, the condensate will form continuous phases and start to flow near the wellbore, and eventually form a multiphase flow.(3) The dynamic data analysis method of production wells in shale condensate gas reservoirs predicts productivity by analyzing the physical properties of oil and gas reservoirs and the production decline law of production wells.Among them, the empirical method is simple to calculate, but it is greatly affected by geography, mining mode and other factors, and its reliability is poor.The characteristic curve analysis method is more widely used in engineering and has strong reliability, but it depends on the selection of seepage model.In addition, the linear solution of nonlinear problems should be studied to improve its applicability when it is applied to shale condensate gas reservoirs involving very complicated seepage law.Artificial intelligence method can process complex data and determine the influence degree of different influencing factors, and it has high prediction accuracy and needs certain basic data as support.
However, the existing production data of shale condensate gas reservoirs are relatively limited in comparison with conventional reservoirs, so it is necessary to add physical constraints to improve the accuracy of production prediction.(4) Mastering the complex seepage law involved in shale condensate gas reservoir development by designing new experiments so as to obtain more accurate data, establishing efficient shale gas reservoir development simulation technology (productivity prediction based on the analysis of dynamic production data, numerical simulation of oil, and gas reservoirs, the injection of CO 2 , N 2 , or flue gas to enhance oil and gas recovery, etc.), understanding the main controlling factors that affect fluid flow and production, and accurately predicting the production performance of production wells are conducive to an efficient development and improvement of economic benefits.First-day decline rate q 1 First-day tight gas well production d Fracture flow index δ Universal exponent

Figure 6 .
Figure 6.Diagram of strong stress sensitivity characteristics of fractures.
et al. improved the NPT-GEMC method through molecular simulation technology to study the constraint effect of phase behavior of different component reservoir fluids in multi-scale (nanoscale + macroscale) media.The results show that a large number of recombinant fractions are trapped in the nanopores, resulting in a much lower condensate yield than expected.In 2021, Zhang Tao [91] et al. described the dynamic process of phase equilibrium based on the P-R equation of state, and developed a set of adaptive deep learning algorithms on this basis.

Figure 10 .
Figure 10.Diagram of reservoir heterogeneity at the SRV area.
Mole fraction of component i in the total hydrocarbon mixture nAmount of total substance in the hydrocarbon mixture n L Amount of total substance in the liquid phase n V Amount of total substance in the gas phase K ij