When acidification radius is greater than or equal to pollution zone radius ( r a ≥ r d ), pollution zone is improved by acidification, the acidification zone become a new “pollution zone”. The radius of acidification zone becomes new radius of “pollution zone”, the permeability of acidification zone becomes a new permeability of “pollution zone”.

Therefore, according to the hypothesis in reference [13], the skin factor in section x of horizontal well is expressed as follows:

S a ( x ) = ( K / K a ) ln [ 1 β + 1 ( r ah ( x ) r w + r ah 2 ( x ) r w 2 + β 2 + 1 ) ] (1)

S a ( x ) is local skin factor, f; K is reservoir permeability, mD; K a is acidification zone permeability, mD; r ah ( x ) is the distance of horizontal well heel in acidification zone, m; β is reservoir anisotropy coefficient, f; r w is wellbore radius, m.

When acidification radius is less than pollution zone radius ( r a < r d ), acids does not penetrate the whole pollution zone. In the original polluted zone, formed a composite zone consisting of acidification zone and polluted zone.

On the yz plane, the horizontal well can be regarded as the vertical well with radius r_w in a reservoir with thickness L, considering reservoir anisotropy, the steady seepage equation on yz plane is as follows:

K z ∂ 2 P ∂ z 2 + K y ∂ 2 P ∂ y 2 = 0 (2)

The internal boundary condition is:

z 2 + y 2 = r w ; P 2 = P w f 2 (3)

On the yz plane, the isobaric lines around the wellbore were concentric ellipses, uv coordinate system was introduced, the isobars around the wellbore become confocal concentric ellipses, as shown in Figure 2.

Introduced uv coordinate system, and conversed coordinate, there are:

u = ( K y / K z ) 1 / 4 z ; v = ( K z / K y ) 1 / 4 y (4)

Formula (2) can be changed to:

∂ 2 P ∂ u 2 + ∂ 2 P ∂ v 2 = 0 (5)

Formula (3) can be changed to:

( K z / K y ) 1 / 4 u 2 + ( K y / K z ) v 2 = r w 2 ; P 2 = P w f 2 (6)

In order to solve the above Laplace equation with elliptic inner boundary conditions, conformal transformation is introduced to transform the elliptic coordinate system into the ρθ coordinate system, there are:

u = b cosh ρ cos θ ; v = b sinh ρ sin θ (7)

When ρ = ρ w at the wellbore, derived from Equation (7):

u 2 b 2 cosh 2 ρ w + v 2 b 2 sinh 2 ρ w = 1 (8)

Comparison formula (5) and formula (7), there are:

tanh 2 ρ w = K z / K y ;     b 2 = r w 2 ( K y − K z ) / ( K z K y ) 1 / 2 (9)

The average radius of the elliptical isobars can be defined as the average of the long and the short axes, there are:

r ¯ = a cosh ρ + a sinh ρ 2 = b 2 e ρ (10)

After coordinate transformation, as long as the average radius of the wellbore and the average radius of the reservoir damage zone are calculated separately, it can be solved by using the formula of skin factor calculation.

Comparison formula (9) and formula (10), the average radius of wellbore is:

r w ¯ = b 2 e ρ w = r w ( K y / K z + 1 ) / 2 ( K y / K z ) 1 / 4 (11)

The anisotropic problem is transformed into isotropic problem by coordinate transformation. Where the anisotropic coefficient β = K y / K z , then formula (11) becomes:

r w ¯ = r w ( β + 1 ) / 2 β (12)

Formula (4) can be changed to:

u = β z ; v = y / β (13)

The long and short axes of isobaric lines in equivalent isotropic formations are obtained from Equation (7) respectively:

u o = b cosh ρ ; v o = b sinh ρ (14)

For the elliptical polluted zone, combination formula (13) and formula (14), when ρ = ρ d , as follows:

ρ d = sinh − 1 ( v d / b ) = sinh − 1 ( r dh / b β ) (15)

Combination formula (10) and formula (15), average radius of elliptical polluted zone is:

r d , eq ¯ = b 2 e ρ d = r w 2 β − 1 / β [ ( r dh r w β 2 − 1 ) 2 + 1 ] (16)

Formation permeability K and composite zone permeability are defined. The expression of anisotropic formation skin factor derived from reference [14] is used:

Combination formula (12) formula (16) and formula (17), skin factor of horizontal well is:

of composite zone is obtained by average of radial flow permeability in reference [15] :

In document [16], the radius of polluted zone and acidification zone distributes heterogeneous along the horizontal well section, and the permeability of polluted zone and acidification zone are constant, as shown in Figure 3. The radius of polluted zone and acidification zone are functions of the distance form horizontal well heel x, ,.

Formula (19) can be changed to:

Formula (18) can be changed to:

is reservoir skin factor, f; K is reservoir permeability, mD; is permeability in composite zone, mD; is radius of pollution zone, m; is reservoir anisotropy coefficient, f; is wellbore radius, m.

In reference [13], while the horizontal well was polluted, the heel of horizontal well is the greatest, the toe is the smallest. From the heel to the toe, the radius of the pollution zone varies linearly. Based on previous studies, the author

assumed that the acidification zone and polluted zone radius d is tribute parabolically along the wellbore. Because the contact time between the toe of horizontal well and acid fluid is relatively short, it can be assumed that the formation is not polluted. When the radius of acidification zone reaches the maximum, the distribution functions of the radius of pollution zone and acidification zone along the wellbore direction are as follows:

is maximum depth pollution invasion, m; is maximum depth acid invasion, m; L is horizontal well length, m.

Combination formula (20) formula (22) formula (23) and formula (21), the distribution pattern of local epidermal factors after acidification of horizontal well can be obtained.

Assumed that the upper and lower boundaries of the reservoir were closed, z_w is the distance between the horizontal well and the lower boundaries of the reservoir, the horizontal section is divided into N micro-segments along the length direction, and the micro-lines on the micro-segments converge radially and uniformly, while the radial flow of micro-segments is different.

According to the mirror reflection principle, the micro-element of the production section in section i() is transformed into the production well row in the infinite formation. From the potential superposition principle, the potential of the micro-element section in section i in the boundary reservoir at any point in space can be obtained as follows:

Among, is potential of section i at any point in the reservoir, is flow of in section i, is potential at the boundary.

Among:;

;

.

The pressure of boundary is, the flow pressure at the heel of horizontal well is, the pressure at the midpoint of section i is, and the potential at the midpoint of section j () is, the relationship between the radial flow rate and the flow pressure. The horizontal well seepage model as follows:

Among, is the pressure of boundary, MPa; is flow of in section i, m³/s; is the pressure at the midpoint of section i, MPa; K is reservoir permeability, mD; is oil viscosity, mPa·s; is oil density, kg/m³; is the distance between horizontal well and lower boundary of reservoir, m; is the distance between upper and lower boundary of the reservoir

Assumed the flow in horizontal wellbore was single-phase incompressible fluid without gravity loss. The pressure loss of the i section in the production section were, includes friction loss caused by wellbore wall, acceleration loss caused by radial inflow of fluid, additional pressure drop caused by skin factor:

The friction loss is:

Among, is friction loss of i section, MPa; D is wellbore diameter, m; is upstream flow of the i section, m³/s; is radial inflow of the i section, m³/s; is oil density, kg/m³; is correction coefficient considering the effect of radial inflow, f; is wellbore friction coefficient, when the

wellbore flow is laminar flow, when, ,

Among, is Re of wellbore inflow, dimension; When the wellbore flow is

turbulent, , ,

Among, is absolute roughness of wellbore, m.

The acceleration loss is:

Selected the segment from distance the heel of horizontal well, the additional pressure drop caused by skin factor is, the flow rate of radial flow into wellbore is, so the formula for calculating the additional pressure drop caused by skin factor is defined as follows:

Among, is the additional pressure drop caused by skin factor, MPa; is the flow rate of reservoir fluid flowing radially into the micro-element section, m³/d; K is the formation permeability, mD.

Combination formula (26) formula (27) formula (29), the pressure loss analysis model of the micro-element section considering the change of skin factor can be obtained:

For production section of horizontal well, there are the following relationships between flow rate and flow pressure:

; (31)

Among,;, is flow pressure at heel of production section,.

Formula (25) formula (30) and formula (31) constitute a closed system of equations containing 2N equations and 2N unknown variables (, ()). It is horizontal well productivity model considering skin factor in the horizontal section after acidification.

The basic parameters of a horizontal well were as follows: horizontal length, wellbore radius, reservoir thickness, formation permeability, pollution zone permeability, acidification zone permeability, formation anisotropy coefficient, oil viscosity, oil density, oil volume coefficient, the pressure at boundary, bottom flow pressure, oil release radius, maximum mud invasion distance, maximum acid action distance.

Used the formula (21) of the skin factor obtained in this paper, combined with the actual data above, the fig of the distribution of the skin factor can be drawn (Figure 4).

It can be seen from Figure 4 that the l skin factor of horizontal well decreases gradually from heel (x = 0 m) to toe (x = 200 m) after acidification, and the skin factor is less after acidification than the skin factor before acidification. It shows that formation damage near wellbore is removed by acidification method, and the skin factor was reduced, achieved the goal of increase the well productivity.

Yield multiplier ratio was usually used to evaluate the acidification effect. The yield increase ratio is the ratio of the yield after acidification to that before acidification. According to the calculation results of the previous productivity model, analyzed the factors affecting the acidification effect of horizontal well.

1) Effect of horizontal well length on yield multiplier ratio of horizontal well

Considered the horizontal well length L = 100 m, 200 m, 300 m, 400 m, 500 m, other parameters in the above examples. The productivity model proposed in this paper is used to calculate yield multiplier ratio of horizontal well. Figure 5 is the relationship curve of the influence of different horizontal section length on the productivity increase ratio after acidizing. As can be seen from Fig, during the increase of horizontal well length, yield multiplier ratio decreases gradually, and the decreasing range tends to be stable. This shows that the shorter the length of horizontal, the better the acidification effect.

2) Effect of reservoir anisotropy coefficient on yield multiplier ratio of horizontal well

Considered the reservoir anisotropy coefficient β = 1, 3, 5, 7, 9 and other parameters in the above examples, used the productivity model proposed in this paper to yield increase ratio of horizontal well. Figure 6 is the relationship curve of the influence of reservoir anisotropy coefficient on yield multiplier ratio of

horizontal well. As can be seen from Figure 6, when other parameters are certain, during the increase of reservoir anisotropy ratio, the yield multiplier ratio is also increasing. This shows that the more serious the formation anisotropy, the better the acidification effect of horizontal well.

3) Effect of pollution depth on yield multiplier ratio of horizontal well after acidification

Considered pollution depth (maximum pollution depth of mud) R_m,max = 1 m, 1.6 m, 2 m, 2.4 m, 2.8 m, and other parameters as the parameters of the above examples. The effect of maximum pollution depth of mud on the yield increase ratio after acidification is shown in Figure 7. As can be seen from Figure 7, when other parameters are certain, the yield increase ratio of acidification

increases with the increase of pollution depth. This shows that the deeper the mud invasion near the wellbore, the better the effect of acidification on increasing production.

4) Effect of pollution degree on yield multiplier ratio of horizontal well after acidification

Considered the damage degree (the formation permeability/polluted area permeability) as 5, 10, 15, 20, 25, and other parameters as the parameters of the above examples, Figure 8 shows the effect of different damage degree on the yield increase ratio of horizontal well after acidification. As can be seen from Figure 8, when other parameters are given, the higher the degree of pollution, the greater the ratio of yield increase of acidification. It shows that the acidification effect of horizontal wells is better when the pollution of horizontal wells is more serious. The more serious the pollution around the wellbore, the better the effect of acidification stimulation is, which is also the reason why acidification stimulation is used in the development of oil and gas fields.