Numerical Simulation Characteristics of Logging Response in Water Injection Well by Reproducing Kernel Method

Reproducing kernel Hilbert space method (RKHSM) is an effective method. This paper, for the first time, uses the traditional RKHSM for solving the temperature field in two phase flows of multilayer water injection well. According to 2D oil-water temperature field mathematical model of two phase flows in cylindrical coordinates, selecting the properly initial and boundary conditions, by the process of Gram-Schmidt orthogonalization, the analytical solution was given by reproducing kernel functions in a series expansion form, and the approximate solution was expressed by n-term summation.The satisfied numerical results were carried out by Mathematica 7.0, showing that the larger the difference between injected water temperature and initial borehole temperature or water injection conditions, the more obvious the indication of water accepting zones. The numerical examples evidence the feasibility and effectiveness of the proposed method of the two phase flows in engineering.


Introduction
As we known, during the middle and late stage of water flooding oil field it is difficult for temperature log to reflect actual water injection profile in normal water injection conditions; at the same time, with the deepening of oil field development and the extension of exploitation time, a phenomenon of the viscosity of oil in the well became stickier and the production is reduced, leading to a large amount of stored oil in underground which can not be mined.Gradually people found that water injection can save this problem.By injecting water, it can keep formation energy, enhance crude oil recovery ratio, and insure long-term high and stable production; at the same time, pressure in oil layer is improved, which can guarantee the high oil field production.
In recent years, many researchers in order to improve temperature log to determine water injection profile, lots of mathematical methods were applied in temperature log such as the ADI method [1], the full implicit form of finite difference method [2], the finite difference method [3], the implicit finite difference method [4], the alternating direction implicit method (ADI), and speedup method [5].
In this paper, RKHSM method is used to simulate the temperature model.It is well known that RKHSM is a simple and accurate method.In recent years, there has been a growing interest in using RKHSM to solve mathematical problems [6][7][8][9][10][11][12][13][14], but there is no one applying this method into characteristics of logging response in two phase flows problems.This paper tries to study the temperature field in water injection with RKHSM.By analyzing the simulated results we can conclude that temperature logging can distinguish water entry zones from non-water entry zones through hot-water injection for short periods of time, enlarging water injection velocity, selecting reasonable shut-in time.It is economical and reasonable way to use short, fast, hot-water injection in order to change the condition that temperature log can not delineate water injection profile.(

Mathematical Modeling
is the temperature,   is the thermal conductivity of water, and   is water saturation. is the depth of the oil wellbore, and  is the radial distance.  is the density of water, and   are the specific heat of water.]  is the velocity of fluid inside the oil wellbore.The initial temperature of the formation and wellbore is at the geothermal condition which can be written as | =0 =  0 =  + ;  is the formation temperature when  = 0,  is the geothermal gradient, the boundary of adiabatic condition is ((  )/)| =0 = 0, the upper boundary of adiabatic condition is ((  )/)| = = 0, and the under boundary of adiabatic condition is ((  )/)| =0 = 0.
In order to solve (1), we introduce the reproducing kernel spaces.

Reproducing Kernel Hilbert Spaces Method (RKHSM).
There are several reproducing kernel spaces (RKHS) to introduce: An inner product is The reproducing kernel is (2 An inner product is The reproducing kernel is ( An inner product is The reproducing kernel is (4) The reproducing kernel space is where  = [0, 1] × [0, 1] and Ṽ() and model problem (1) changes to the following problem: () (, , ) =  (, , ) .

Numerical Experiment.
In order to show RKHSM is effective, a nonlinear equation is used as a test to prove the high accuracy of the method.Select the above conditions and boundary conditions and let   = 1 and  = (0, 1) × (0, 1); the equation has an exact solution of the form (, , ) =  − ⋅ sin 10 ⋅ cos 10.By Mathematica 7.0, using our method, we choose the number of nodes as 23 and obtain the approximate solution  23 ; the comparison of the results with the exact solution is shown in Table 1, and they are found to be in good agreement with each other.Figure 1 is the exact solution (, , 1), and Figure 2 is the error ( −  23 ).

Numerical Simulation Results and Discussion
Using the common body structure of vertical wells [5], assuming that the fluid is Newton fluid, the radius of oil well is 0.1 m, and the geothermal gradient is 0.025 ∘ C/m; ignoring the influence of casing and cement, the initial formation temperature is 15 ∘ C. The interlayer between the two injections is 20 m.The thermal conductivity of water and oil is 4182 J/kg⋅K and 2234.2J/kg⋅K, the density of water and oil is 1000 kg/m 3 and 800 kg/m 3 , and the specific heat of water and oil is 0.618 W/m⋅K and 0.156 W/m⋅K.The boundary of adiabatic condition is | =0 = 15 ∘ C, ((  )/)| =0 = 0, the upper boundary of adiabatic condition is ((  )/)| =1 = 0, and the under boundary of adiabatic condition is ((  )/)| =0 = 0.

Effect of Different Injection
Rate.The injection temperature is 80 ∘ C; temperatures of numerical simulations are given in Figure 1 at  = 1 h when the injection rate is 0.0059 m/s, 0.095 m/s, and 0.0368 m/s.From Figure 3, it can be seen that the temperature of injection layer is more changed than others when ]  = 0.0059 m/s.The temperature in injection layer is obviously changed when the injection rate is small, because the fluid and borehole wall exchange the heat transfer sufficiently in the wellbore.

Effect of Different Injection
Temperature.When injection rate ]  = 0.0059 m/s, the temperature of numerical simulation when injection temperature is 20 ∘ C, 60 ∘ C, and 80 ∘ C is given in Figure 4 at  = 1 h.
Figure 4 shows the temperature profile is less obvious when the injection temperature is close to the original formation temperature (42.5 ∘ C), and the temperature profile is more obviously changed when the injection temperature is further deviated from the original formation temperature (20 ∘ C, 80 ∘ C).

Effect of Different Injection Time. When injection rate
]  = 0.0059 m/s, the injection temperature is 80 ∘ C, and the temperature of numerical simulation when injection time is 0.1 d, 10 d, and 100 d is given in Figure 5. Figure 5 presents that the temperature of injection layer is less changed compared with others when  = 0.1 d; this is because the fluid and borehole wall do not exchange the heat transfer sufficiently.As the temperature of fluid and borehole wall becomes consistent when injection time is long, temperature profile can not be distinguished.Therefore, reasonable injection time is good for temperature field of the recovery.

Effect of Different Oil Layer Thickness. When injection rate
]  = 0.0059 m/s, the injection temperature is 80 ∘ C, and injection time is  = 1 h; the temperature of numerical simulation when oil layer thickness is 1 m, 6 m, and 10 m is given in Figure 6.
In Figure 6, when oil layer thickness is 1 m, the temperature of injection layer is less changed compared with others.When the oil layer is thick, the temperature is high, because Mathematical Problems in Engineering the fluids are fully blended in the entrance of the injection layer.

Conclusions
In this paper, RKHSM is used to solve a class of linear or nonlinear partial differential equations which can describe the temperature field in water injection well; it can be seen that the method is feasible.According to the numerical simulation results in Section 3, the following conclusions are given.
During the water injection wellbore production, different injection-production conditions make different temperature profile.By analyzing the simulated results we can conclude that temperature logging can distinguish water entry zones from non-water entry zones through hot-water injection for short periods of time, enlarging water injection velocity, selecting reasonable shut-in time, so we can use the short, fast, hot-water injection in order to change the condition that temperature log can not delineate water injection profile.