Modelling of In-Situ Scale Deposition: The Impact of Reservoir and Well Geometries and Kinetic Reaction Rates

Abstract

Abstract Previous work has demonstrated how and where the mixingof incompatible brines occurs in waterflooded reservoirs, and what the impact would be on scale prevention strategies in terms of timing and placement of squeeze treatments. This paper extends this work, by modelling the resulting in-situ depositionprocess. The location of maximum scale deposition and the resulting brine compositions at the production well are calculated for a range of sensitivities, including reservoir geometry (1D, 2D areal, 2D vertical, 3D), well geometry (location and orientation within field and with respect to other wells and the aquifer), and reaction rate (ranging from no precipitation to equilibrium). In conventional systems with no aquifer, it is demonstrated that maximum scale deposition occurs in the immediate vicinity of the production wellbore, and therefore low produced cation concentrations indicate inadequate squeeze treatments. In systems where water injection is into the aquifer, low cation concentrations may also result from deposition deeper within the reservoir. Maximum scale dropout still occurs as the fluids approach the production well, but sufficiently far from the wellbore to be unaffected by squeeze treatments, or to have any major impact on productivity. The reaction rate is critical in determining the amount of scale deposition, but even under equilibrium conditions, sufficient concentrations of scaling ions are delivered to the production well to necessitate squeezing the well, although using lower volumes of inhibitor. Once cation concentrations have been reduced, it is predicted that they will never pick up again. This paper also discusses some of the limitations of modelling such systems, which include the determination of the kinetic reaction rates, size of the mixing zone, and the impact on permeability. Although the thermodynamics are fairly well understood, the kinetics are much more difficult. The size of the mixing zone is affected by numerical dispersion, and computationally intensive techniques are required to overcome this problem. Previous experience shows that formation damage factors are very difficult to extrapolate from coreflood data because there is a great difference between the dimensions of the mixing zone in the reservoir and the core plug. Introduction Previously presented work has shown the effect of brine mixing under various flow conditions, both idealised1,2(1D, 2D vertical, 2D areal, 3D) and actual reservoir conditions3,4(Alba Reservoir, North Sea). For the majority of these calculations a conventional reservoir simulator has been used. The advantages of using a conventional simulator are that a large proportion of waterflooded reservoirs have a field model dataset already available, the addition of tracer tracking to model the propagation of the mixing zone is relatively straightforward to implement, and the results are easy to visualise. This technique is quite powerful for demonstrating the movement of the water front and also the mixing zone relative to the production wells. Even within a given field the behaviour may vary quite markedly, depending on the configuration of neighbouring wells and the reservoir geometry, as was shown in the Alba case3.