A New Method for Research on Unsteady Pressure Dynamics and Productivity of Ultralow-Permeability Reservoirs

In the numerous low-permeability reservoirs, knowing the real productivity of the reservoir became one of the most important steps in its exploitation. However, the value of permeability interpreted by a conventional well-test method is far lower than logging, which further leads to an inaccurate skin factor. This skin factor cannot match the real production situation and will mislead engineer to do an inappropriate development strategy of the oilfield. In order to solve this problem, key parameters affecting the skin factor need to be found. Based on the real core experiment and digital core experiment results, stress sensitivity and threshold pressure gradient are verified to be the most influential factors in the production of low-permeability reservoirs. On that basis, instead of a constant skin factor, a well-test interpretation mathematical model is established by defining and using a time-varying skin factor. The time-varying skin factor changes with the change of stress sensitivity and threshold pressure gradient. In this model, the Laplace transform is used to solve the Laplace space solution, and the Stehfest numerical inversion is used to calculate the real space solution. Then, the double logarithmic chart of dimensionless borehole wall pressure and pressure derivative changing with dimensionless time is drawn. The influences of parameters in expressions including stress sensitivity, threshold pressure, and variable skin factor on pressure and pressure derivative and productivity are analyzed, respectively. At last, the method is applied to the well-test interpretation of low-permeability oil fields in the eastern South China Sea. The interpretation results turn out to be reasonable and can truly reflect the situation of low-permeability reservoirs, which can give guidance to the rational development of low-permeability reservoirs.


Introduction
As the exploration of oil and gas gets deeper, the scale of proven reserves in offshore low-permeability oilfields gets larger and larger. Many topics worth study show up in offshore low-permeability reservoirs [1]. One of the most tricky issues is the precise evaluation of the physical parameters, which are variant and hard to predict in the lowpermeability reservoirs. The accurate evaluation model needs to be established and can give guidance to the future development of offshore reservoirs.
The seepage mechanism of low-permeability reservoirs is very complex. Due to the stress sensitivity effect and threshold pressure gradient in low-permeability reservoirs [1][2][3], the current research on the theory and application of seepage flow in offshore low-permeability oil reservoirs has limitation, including two common features: (i) as the two variables, stress sensitivity effect and threshold pressure gradient, both have big influences on seepage, most models only consider one of them as factor; (ii) in order to prevent the blowout, using of the mud with high density in offshore drilling and mud pollution caused by this method is not taken into consideration in most models.
After years of exploration, researchers all over the world have done a lot of work figuring out the seepage mechanism of low-permeability reservoirs. Stress-sensitive effect and threshold pressure gradient are accepted as the most influential factors. At present, there are two main ways to consider the influence of stress sensitivity: one is to redefine a pseudopressure, and the other is to derive from permeability variation function as the productivity of oil and gas well is usually affected by formation permeability [4]. Dou et al. introduced a new program to quantitatively analyze the stress sensitivity of permeability by using the traditional linear analysis method and determined the permeability modulus and other formation parameters based on the concepts of modified pseudoparameters and progressive approach [5]. Many scholars have also done a lot of research work in the stress-sensitive unstable seepage theory and the application of multiple fractured horizontal wells in unconventional oil reservoirs [2,[6][7][8][9]. When dealing with the nonlinear term of seepage differential equation, Pedrosa's variable substitution and perturbation technique are commonly used [10][11][12]. Zhang and Ambastha [13] further studied the influence of stress-sensitive permeability on unstable pressure dynamics by using the subsection permeability model. Wu and Pruess have applied the integral method [13,14]. Zeng and Zhao established an approximate equation to analyze the unstable pressure performance of reservoirs [15]. Zhang et al. put forward a well-test interpretation model of natural fractured reservoir considering stress sensitivity and threshold pressure gradient, but the model only got numerical solution, not analytical solution [16]. Qanbari and Clarkson provided a method for analyzing the transient linear flow in stress-sensitive reservoirs by using production data, but this method could not stably display all flow stages of pressure dynamic state (such as quasisteady flow stage) [17].
Several phenomena of non-Darcy flow in lowpermeability reservoirs (threshold pressure gradient, lowvelocity non-Darcy) can be explained theoretically from the perspective of force [18]. When the pressure is particularly high, the adsorption of liquid molecules on the rock surface is strong, and the liquid molecules adsorbed on the rock surface cannot move. The higher the pressure, the lower the liquid permeability, which is the characteristic of threshold pressure gradient. Malhasin discussed the principle of threshold pressure gradient generation from the perspective of micromolecular structure and force [18]. Miller and Low discovered the phenomenon of threshold pressure gradient when discussing the seepage of water in clay [19]. Pascal et al. considered the influence of the differential pressure threshold of low-speed flow and solved the consolidation problem in geotechnical engineering using the finite difference method in the numerical method [20]. Prada and Civan proposed that Darcy's formula should be corrected for the threshold pressure gradient [21]. Zeng et al. measured the small threshold pressure gradient value of ultralowpermeability reservoirs and proposed a method to determine the value using the pseudothreshold pressure gradient and nonlinear coefficients [3]. Other scholars have also done a lot of work in experimental methods and new algorithms [22][23][24][25][26], which provides a strong theoretical basis for the study of seepage laws in low-permeability reservoirs.
The cost and risk of offshore drilling are higher than land drilling. To prevent blowouts, the pressure of the mud column must be greater than the reservoir pressure; thus, mud particles and filtrate can enter the reservoir. When the filtrate is not compatible with the minerals of the reservoir, especially the clay mineral fluid, it will cause the clay minerals to swell, disperse, precipitate, etc. If blockage of the flow channel happens, it will result in the decrease of permeability and increase of skin surface and even can lower the productivity of the reservoir. For wells in low-permeability reservoirs, when mud invasion happens, the particles in the formation gradually accumulate in the near-well zone. Thus, the skin factor will slowly increase obviously as time passes. Therefore, it is inappropriate to use constant skin factor in the models, which will come up with wrong well-test interpretation results [27]. However, in the current mathematical model of well-test interpretation, skin factor is usually regarded as a constant [28].
Therefore, we consider seepage in offshore lowpermeability oilfield as a comprehensive result of threshold pressure gradient, stress sensitivity, and mud pollution. In this paper, the author established the mathematical model of unsteady seepage in offshore low-permeability oilfield by using the law of conservation of mass and Laplace transform. It is developed not only considering the influences of threshold pressure gradient and stress sensitivity but also considering the influence of mud pollution in the actual seepage process of offshore low permeability. By drawing the typical curves with different influencing factors by numerical inversion and programming, their influences are analyzed. In addition, from the result of the interpretation in actual lowpermeability oilfield, this model is more reasonable than conventional methods and can give guidance to the rational development of low-permeability oilfield.

Establishment and Solution of Unsteady
Seepage Mathematical Model for Low-Permeability Oilfield 2.1. Influencing Factors. In order to further analyze the influence of the three factors of threshold pressure gradient, stress sensitivity, and mud pollution and to study reasonably mathematical representation methods, this paper first uses the offshore low-permeability core experiment to study the expressions of stress sensitivity and threshold pressure gradient. In addition, the mud intrusion process is modeled by a numerical core simulation method, and the reasonable expression of three factors is determined.

Geofluids
It can be seen from the above figure (Figures 1-3) the exponent type has a wide range of changes in the early and middle stages. The change rates of power formula are larger in the early stage and smaller in the later stage. The overall change of the S S method is relatively smooth and uniform, the rate of change in the previous stage is smaller than that in the former two, and there is no smooth section in the two later stage's curves. Through the above analysis, it is concluded that the exponential form fitting is the best, so the exponential form method is used to fit the experimental data. Figure 4 shows the change of permeability of 4 cores during pressure increasing and pressure decreasing stages and proves that permeability is sensitive to effective stress. It can be seen from the figure that the lower the permeability, the greater the permeability decline rate when the effective stress increases: sample No.6 has the largest permeability decline rate, while sample No.10 has the smallest permeability decline rate. Therefore, the hysteresis effect is obvious, and the lower the permeability is, the more obvious the nonlinear effect is. The confining pressure range 38-60 MPa is measured, and 60 MPa is used as reference point to fit the experimental data.
The exponential model is the most common mathematical model of permeability variation. Assuming that the deformation of reservoir fluid viscosity and permeability is changing due to stress sensitivity accords with the Hooke elastic rheological law, similar to the definition of compressibility, Pedrosa defines the permeability modulus, which is used to describe the relationship between reservoir permeability and pressure and can be written as homogeneous medium.
The results of stress sensitivity test got the relationship between stress sensitivity coefficient and reference point permeability. As shown in the figures, the stress sensitivity coefficient and the reference point permeability satisfy the exponential relationship, and the correlation is good. The weaker the nonlinear effect, the apparent permeability changes linearly with the pressure. The smaller the permeability is, the stronger the nonlinear effect is, and the more obvious the influence of the pore compression on the flow. The nonlinear effect in the boosting process is much greater than that in the depressurization process.
Stress-sensitive coefficient method    Through the derivation process, it can be seen that c 1 reflects the existence of the yield stress of the fluid and the influence of the boundary layer on seepage, and c 2 mainly reflects the influence of the boundary layer on seepage. Because of the microscale flow effect, the influence of surface force on seepage cannot be ignored. The greater the surface force, the greater the binding force of boundary layer on fluid, which makes the fluid show stronger non-Newtonian and finally leads to the more deviation of seepage law from Darcy's law. At the same time, the greater the binding force, the greater the offset displacement force, making the fluid more easily adsorbed on the boundary layer, making the boundary layer thicker and further aggravating the microscale flow effect, which is a mutual coupling effect.
13 cores of WC oil group in LF were taken, and 7 of them were used in the threshold pressure gradient experiment of oil phase. In the flow experiment of low-permeability cores, the "pressure difference-flow method" is used. In the Cartesian coordinate system of flow and pressure gradient, not only a straight line is presented but also an upturned curve and a straight line are formed, so as determining the threshold pressure gradient. As shown in Figure 5, the experimental results of each core are fitted based on the nonlinear seepage model equation, the slope of the straight line segment is the permeability, and the intercept is the threshold pressure gradient.
The permeability obtained by experimental fitting (slope of linear segment) is different from that measured by gas. It is generally smaller because of the use of crude oil as medium. The permeability measured by liquid is generally lower than that measured by gas (absolute permeability). The lower the absolute permeability is, the greater the decrease of liquid permeability will be. This reflects that the denser the reservoir, the greater the additional resistance to liquid. Therefore, for low-permeability tight reservoirs, the relationship between fluid seepage velocity and pressure gradient is no longer a straight line passing through the origin. The velocity equation of low-speed non-Darcy seepage considering the 4 Geofluids threshold pressure gradient can be described by the following pseudothreshold pressure gradient model: 2.1.3. Study on Pollution Impact Caused by Offshore Mud Intrusion. Figure 6 shows the core numerical simulation. In order to clarify the law of core pollution by drilling fluid, a core numerical model is established, the experimental results of core pollution by drilling fluid are fitted, and the correla-tion between core damage degree and drilling fluid flux is obtained.
The basic experimental information is as follows: inject polymer drilling fluid into the core and measure the core permeability at the same time. After injecting the polymer drilling fluid, the core permeability decreased from 13.0 mD to 5.4 mD, and the permeability loss was 58.3%. The coredamaged model obtained by fitting is as follows: The radial grid model of single well is established to simulate the invasion of drilling filtrate into formation, and the

Geofluids
formation mechanism of pollution zone and its influence on productivity are studied. Combined with the basic geological characteristics of the target area, the theoretical model is established. Figure 7 shows the recovery under different conditions. The following two situations are compared: polluting reservoir and nonpolluting reservoir after mud filtrate immersion in reservoir.
(1) If the pollution of drilling fluid around the well is not considered, the invasion of drilling fluid will be little further, but this shows very tiny difference. (2) In all reservoirs with different original permeability value, the recovery degree will decrease when pollution happens. Figure 8 shows the influence of soaking time on pollution. Through simulation, it can be found that the higher the permeability, the smaller the skin caused by stress sensitivity; the larger the skin caused by pollution, the larger the pollution radius; the more mud invasion, and the higher the recovery loss.
As shown in Figures 9 and 10, reservoir pollution caused by mud invasion during drilling is usually the result of physical and chemical effects. If the average permeability of polluted area is regarded as a function of time, it will change obeying the rule as follows: (i) in the early stage, because the fine particles in the mud fluid can easily accumulate near the well zone, the average permeability in the polluted area will change rapidly; (ii) as time passes by, the number of this kind of particles is getting smaller, so the change of permeability will get slower; and (iii) the average permeability of the polluted area will become a constant. According to the results of numerical core fitting and the shape of pollution invasion curve, there is a certain exponential relationship between soaking time and permeability damage; the law is summarized as follows:   6 Geofluids is C. Non-Darcy radial seepage occurs in low-permeability reservoirs.

Establishment and Solution of Mathematical Model.
Since the threshold pressure gradient and stress sensitivity of low-permeability fluid cannot be ignored, and porous media and fluid are both compressible, it is necessary to consider them as porous media and elastic fluid. The basic differential equation of unstable seepage in low-permeability reservoir considering two factors at the same time can be obtained as follows (Appendix A): The above formula is a partial differential equation with strong nonlinearity. Since the simplified partial differential equations are still relatively nonlinear, in order to facilitate the production of the plates and form a unified pressure and pressure derivative plate, a dimensionless conversion is required. The following dimensionless definition is introduced (Table 2).
After dimensionless transformation, the simplified differential equation of dimensionless seepage flow considering both the threshold pressure gradient and stress sensitivity is obtained:

Geofluids
There is a gradient square term on the left side of the above equation, which is a nonlinear differential equation and cannot be solved directly. Make the following transformation: According to the canonical perturbation theory, considering that the dimensionless permeability modulus is usually very small (δ D ≪ 1), the zero-derivative perturbation solution can meet the engineering accuracy requirements, so it can be simplified: After conversion, the corresponding definite conditions under the conditions of skin factor and well storage are External interface condition : lim The influence of mud contamination is not considered in the current low-permeability mathematical model, and the corresponding initial conditions are ξ wD ðT D Þ = ξðR D , T D Þ − Sð∂ξ/∂R D Þj RD=1 . When the mud invasion is not considered, the skin s is usually regarded as a fixed value. In this paper, considering the influence of mud pollution, the difference between the initial time of mud threshold and the pressure recovery test is T 0 , and the dimensionless immersion time is introduced: In the mud invasion process, the dimensionless permeability initial rate of change is u iD = u i ϕ i μðC i + C m ÞR 2 w /K i . Above the partial differential equation is a nonhomogeneous second-order linear ordinary differential equation, where −ðλ BD /sr D Þ is a nonhomogeneous term, and the general solution of the inhomogeneous equation is the sum of the general solution and the particular solution of the inhomogeneous equation. In this paper, the constant variation method is used to find the particular solution (Appendix A).
In order to describe the permeability recovery characteristics of polluted areas, the skin factor is defined as Assuming that the radius of the polluted area is constant, that is, R s is a fixed value, defining the skin factor S i at the initial time (t = 0) as a constant, and S is a function of time (because K s is a function of time), the relationship between S and S i can be established by using the relationship between K s and K si ; then, there is a relational expression: According to the relationship between formula (10) and formula (11), the following relationship can be obtained: Thus, the relationship between the skin and time can be obtained: A partial differential equation can be obtained. Based on the partial differential equation, Laplace transform is applied to obtain the partial differential equations and boundary conditions: Based on the initial condition and formula (15), we can get:

Dimensionless radius
La is the symbol of laplace transformation, and the dimensionless bottom-hole flowing pressure in the Laplace space can be obtained, taking into account the threshold pressure gradi-ent, stress sensitivity, mud pollution (variable skin), and the influence of well storage: Among them are In order to further illustrate the scalability of the model, the model is simplified by using an asymptotic solution. When u iD and G D approach 0, only stress-sensitive models are considered. The simplified result is the same as the textbook [29]; u iD and α approach what value? When it is 0, it is a model that only considers the threshold pressure gradient. The simplified result is in consistent with the textbook [29], which confirms the scalability of the paper.

Analysis and Application of Well-
Test/Productivity Curve Considering Stress Sensitivity, Threshold Pressure, and Reservoir Damage 3.1. Typical Curve Analysis. By using the Stehfest numerical inversion method, the dimensionless bottom-hole pseudopressure in the Laplace space is inverted, and the double logarithmic typical curve of unsteady seepage bottom-hole pseudopressure in low-permeability reservoir considering stress sensitivity, threshold pressure, and reservoir damage can be obtained by programming (Figures 11-16). Moreover, the threshold pressure gradient, the threshold pressure gradient, stress sensitivity, the initial gradient rate of permeability, and the influence of a on typical curves are considered, respectively. It can be seen from the productivity curve in Figure 11 that the threshold pressure gradient has a great influence on productivity, and the greater the threshold pressure 9 Geofluids gradient, the faster the productivity decreases. Figure 12 shows that the pressure and pressure derivative curves increase with the increase of the threshold pressure gradient (G), and the pressure derivative curves is close to 0.5 horizontal line when not considering the influence of the threshold pressure gradient. The existence of the threshold pressure gradient makes the pressure derivative curves of lowpermeability reservoirs no longer tend to the 0.5 horizontal line. When the threshold pressure gradient gets larger, the upward warping becomes bigger. Therefore, the influence of threshold pressure gradient in low-permeability reservoirs cannot be ignored.
From the productivity curve in Figure 13, the stress sensitivity mainly affects productivity in the middle and late periods. The greater the stress sensitivity, the earlier the impact and the faster the productivity decreases. Figure 14 shows that in low-permeability reservoirs, the pressure and pressure derivative curves rise with the increase of threshold pressure gradient (G) and stress sensitivity (δ). The larger threshold pressure gradient (G) and stress sensitivity (δ) indicate a more difficult condition of the fluid flow. As the curve upward warping gets higher, the bigger influence of threshold pressure gradient and stress sensitivity in lowpermeability reservoirs should be considered.
It can be seen from the productivity curve in Figure 15 that the mud pollution coefficient has great influence on productivity. Mud invasion has an impact on productivity from the beginning, and the more serious the invasion, the faster the productivity declines. Figure 16 shows that in lowpermeability reservoirs, the greater the mud pollution (α), the greater its influence on permeability will be. This result indicates that more serious mud pollution and larger skin

Geofluids
show similar effect as the skin mainly affects the early hump stage of the pressure derivative curve. All three factors have certain influence on productivity and should not be ignored in well-test interpretation. Moreover, since the influence stages of each factor on productivity and pressure are different, chart analysis can have certain guiding significance for low-permeability oilfield development.

Comparison of Interpretation Results of Examples.
In order to illustrate the universality of the model, the model is applied in the conventional oil reservoir. Taking the offshore HZ reservoir as an example, the permeability of the commercial well-test interpretation software is 395.8 mD. The result explained in this paper model is 400 mD and the explained skin is -1.4. The explained results are very close to reality. This demonstrates the rationality of the model derived in the article to be interpreted in conventional oil reservoirs.
Through the analysis of multiple low-permeability wells offshore, the upturn at the end of the pressure derivative curve is typical in low-permeability reservoirs. However, when fitting the curve, the most conventional well-test interpretation software adopts the radial composite model without considering the influence of the special threshold pressure gradient, stress sensitivity, and mud pollution in the low-permeability formation. It results in a serious discrepancy between the well-test interpretation values and the actual situation, which is mainly manifested in two aspects: as shown in Table 3: (1) when there is no stimulation measures, the skin interpreted of some wells are negative; (2) the permeability of well-test interpretation is far lower than that of logging.
To solve the problems, the offshore low-permeability wells are reinterpreted using our model. Taking an offshore low-permeability reservoir as an example, the basic parameters of DST well-test interpretation are shown in Table 4.
According to the basic parameter (Table 4), the parameter fitting range is determined, and the history fitting is finished. Figure 17 shows the good fitting results. Table 5 is the basic formation parameter information of lowpermeability reservoir after fitting. Table 6 is a comparison of the interpretation results of the method in this paper and the conventional method. First of all, the permeability interpreted by this method is closer to that interpreted by logging than the conventional method. In addition, the skin interpreted by the conventional method is -1.01, which is inconsistent with the actual situation because the well has no stimulation measures and mud pollution. The method in this paper is more in line with the percolation theory of low-permeability reservoir. The skin after pollution is explained to be 5.5, which is more reasonable than the commercial well-test software. Permeability and skin affect oil reservoir production allocation, which is a key part in the development and production process. Skin factor in near-well zone is an important parameter to evaluate oil well productivity and completion efficiency. Through skin factor, formation conditions can be understood more deeply.

Skin Decomposition.
Taking the XJ Oilfield in this paper as an example, considering the influence of threshold pressure gradient, stress sensitivity, and mud invasion, the interpreted skin factor is bigger than the conventional method, which can provide a more accurate basis for reservoir protection and stimulation measures. It makes the calculation of pollution range and degree, and the prediction of stimulation rate and output more accurate. By applying good reservoir protection measures, the output of single well and oilfield can be increased, which will be helpful to the reasonable development evaluation of low-permeability oilfield and promote the development of low-permeability oilfield. Figure 18 shows the relationship between skin and invasion time during mud invasion process. The ultimate goal of formation testing is to obtain reservoir seepage parameters, formation pressure, and production information. Engineers can use skin factor to see the damage degree of reservoir. Based on that, the effect of increase measures can be evaluated by predicting the increase oil rate and giving evidence of the implementation. The skin factor obtained by formation test is usually the sum of the real skin damage factor and various pseudo-skin factors caused by various factors during drilling and completion. If the sum of the skins is used in the evaluation of reservoir damage degree, it will inevitably lead to deviation to the real value and poor effect of the stimulation measures. The total skin factor of reservoirs changes with different physical properties. Only through the fine decomposition of skin factor can the skin effect produced by each link be clearer, and the main reasons and links that affect the total skin factor can be found.
Based on the method provided in this paper, the relationship between skin and mud invasion time can be obtained by curve fitting. The initial skin is 5.5, and the total skin after the final mud invasion is 6.87. According to the relationship curve, the mud invasion skin (s) is equal to 1.37, and the ratio of mud invasion skin is 20.0%. Figure 19 shows the factors affecting the XJX-2Sa well. In order to further decompose the threshold pressure skin and stress-sensitive skin, so as to evaluate the influence of various factors, it is necessary to decompose the initial skin 5.5, which can be realized by the following steps: (1) Based on the above fitting parameters, when the skin is 5.5, the total pressure drop ΔP considering the threshold pressure and stress sensitivity is obtained (2) When the threshold pressure skin is S1 (stress sensitivity is 0), the pressure drop only considering the threshold pressure gradient (ΔP 11 ) is calculated. At the same time, when the stress-sensitive skin is 5.5-S1 (the threshold pressure gradient is 0), the pressure drop only considering stress sensitivity (ΔP 12 ) is obtained (3) When the threshold pressure skin is S2 (stress sensitivity is 0), the pressure drop only considering the threshold pressure gradient (ΔP 21 ) is calculated. At 11 Geofluids the same time, when the stress-sensitive skin is 5.5-S2 (the threshold pressure gradient is 0), the pressure drop only considering stress sensitivity (ΔP 22 ) is obtained (4) In this way, different ΔP n1 and ΔP n2 can be obtained through a large number of different S n . By calculating ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ½ðΔP 2 − ðΔP n1 + ΔP n2 Þ 2 Þ q , when the standard deviation (S n ) is the smallest, S n is the S * that we need to get, so that we can get the decomposed threshold pressure skin S * and stress-sensitive skin (5:5 − S * ). In this paper, the threshold pressure gradient skin is 4, and the stress-sensitive skin is 1.5. The threshold pressure skin accounts for 58.2%, and the stresssensitive skin accounts for 21.8%.    As shown in Figure 20, in order to further confirm the proportion of the various factors and verify the rationality of the method, this paper is verified by the numerical simulation method. On the basis of a low-permeability geological model of an oil reservoir in the South China Sea, the stress sensitivity and threshold pressure gradient laws of the target area are added to the numerical simulator, and the process of drilling filtrate invading the formation is simulated for 11 times, with permeability values ranging from 1 mD to 40 mD. Combined with the pressure curve under the condition of pollution, the skin factor is calculated and decomposed, and the pollution is quantitatively evaluated, so that the skin maps corresponding to reservoirs with different permeability can be obtained. Then, the skin factor obtained under different permeability is normalized by the total skin, and the contribution graph of different factors to the skin can be obtained.
As shown in Figure 21, we can see that the permeability 5.33 mD calculated by this method is about 19%, the stresssensitive skin is about 26%, and the threshold pressure is about 55%. When the permeability is 5.33 mD, the threshold pressure skin accounts for the largest proportion, the stresssensitive skin is the second, and the mud pollution skin is the smallest. The results of simulation match with the result using mathematic model in this paper, which proves that the method in this paper has wide applicability.
The analysis of the application effect of skin factor decomposition shows that the causes of various skin effects and the proportion of injuries caused by various factors can be clarified through the decomposition of skin factor. Taking the XJX-2 well in this paper as an example, we can see the proportion of skin-influencing factors in this 4low-      13 Geofluids permeability well. In order to ensure the full release of production capacity of low-permeability well, it is suggested that the threshold pressure gradient should be appropriately reduced by water injection and fracturing. Therefore, the skin factor decomposition technology has practical value in reservoir engineering.

Conclusion
(1) Since the effects of threshold pressure gradient, stress sensitivity, and reservoir pollution on seepage flow in low-permeability reservoirs cannot be ignored, to more precisely explain and predict its flow mechanism, in this paper, a mathematical model of unsteady pressure seepage flow is derived and established (2) In this paper, the numerical model is solved by the Laplace transform, and the double logarithmic chart of dimensionless borehole wall pressure and pressure derivative changing with dimensionless time is drawn by programming. After studying the sensitivity of the factors, it is found that threshold pressure gradient, stress sensitivity, and reservoir pollution in low-permeability reservoirs have different effects on different stages of pressure performance and should be considered, respectively (3) If not considering the threshold pressure gradient, stress sensitivity reservoir pollution, and the permeability, the skin factor calculated by the model will be lower than reality. Considering these factors, the interpretation results will be more reasonable. The results are applied to the well-test interpretation of the low-permeability oilfield in the eastern South China Sea, which successfully reflect the actual situation of low-permeability reservoirs. It has a certain guiding significance for the rational development of low-permeability reservoirs (4) The fine decomposition of skin factor can give more accurate evidence for applying reservoir protection and stimulation measures. It can also help the prediction of the effect of the measures, including the calculation of pollution range and degree and the prediction of stimulation rate and output. By separately evaluating the influence of each factor, the real effect and level each factor contributes to the change of the total skin factor can be clearly found.
Because porous medium and fluid are compressible, the state equation of porous medium and elastic fluid is needed. ðA:4Þ Consider that the compression coefficient is much smaller than the stress sensitivity coefficient (c m ≪ γ). To further simplify the seepage differential equation, we can get Þe δ Pi−P ð Þ ∂P ∂T : ðA:5Þ According to the Laplace formula (14), initial condition (15), boundary condition formula (16), and formula (17), the partial differential equation and the condition of definite solution in the Laplace space can be obtained.
The upper partial differential equation is a nonhomogeneous two-order linear ordinary differential equation.−G D /zR D . For nonhomogeneous terms, the general solution of the nonhomogeneous equation is the sum of the general solution of the corresponding homogeneous equation and the special solution of the nonhomogeneous equation. The corresponding homogeneous linear ordinary differential equation of order two is