Next Article in Journal
A Power Hardware-in-the-Loop Based Method for FAPR Compliance Testing of the Wind Turbine Converters Control
Next Article in Special Issue
Viscosity Models for Drilling Fluids—Herschel-Bulkley Parameters and Their Use
Previous Article in Journal
Methane Hydrate Formation and Dissociation in Sand Media: Effect of Water Saturation, Gas Flowrate and Particle Size
Previous Article in Special Issue
Experimental Investigation of Inhibitive Drilling Fluids Performance: Case Studies from United States Shale Basins
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Methods for Investigation of Drilling Fluid Displacement in Irregular Annuli

1
SINTEF Industry, S.P. Andersens vei 15 B, N-7031 Trondheim, Norway
2
Department of Energy and Petroleum Engineering, University of Stavanger, N-4036 Stavanger, Norway
*
Author to whom correspondence should be addressed.
Energies 2020, 13(19), 5201; https://doi.org/10.3390/en13195201
Submission received: 31 August 2020 / Revised: 29 September 2020 / Accepted: 1 October 2020 / Published: 6 October 2020
(This article belongs to the Special Issue Advances in Drilling Fluid Technology)

Abstract

:
Experimental methods are still indispensable for fluid mechanics research, despite advancements in the modelling and computer simulation field. Experimental data are vital for validating simulations of complex flow systems. However, measuring the flow in industrially relevant systems can be difficult for several reasons. Here we address flow measurement challenges related to cementing of oil wells, where main experimental issues are related to opacity of the fluids and the sheer size of the system. The main objective is to track the propagation of a fluid-fluid interface during a two-fluid displacement process, and thereby to characterize the efficiency of the displacement process. We describe the implementation and use of an array of electrical conductivity probes, and demonstrate with examples how the signals can be used to recover relevant information about the displacement process. To our knowledge this is the most extensive use of this measurement method for studying displacement in a large-scale laboratory setup. Optical measurements and visual observations are challenging and/or costly in such large-scale systems, but can still provide qualitative information as shown in this article. Using electrical conductivity probes is a robust and fairly low-cost experimental method for characterizing fluid-fluid displacement in large-scale systems.

Graphical Abstract

1. Introduction

Physical experiments remain an indispensable tool for many fluid mechanical problems due to the complexity or size of the problem. This is also the case for the fluid-fluid displacement which is part of a primary cementing job where drilling fluid is displaced by well cement, only separated by prewashes and spacer fluids. The process of displacement with cement is properly described by Daccord et al. [1]. The sheer size of the physical domain of interest, typically a section of wellbore annulus of several hundred meters, makes full three-dimensional CFD simulations far too computer-intensive. Simplified models have therefore been developed to increase simulation speed, while capturing the main physics [2,3,4]. Such models must, however, be validated against experiments which, in turn, represent downhole conditions as realistically as possible. Furthermore, to avoid excessive waste handling issues, it is necessary to use laboratory fluids (non-hardening model fluids) instead of reactive fluids like cement slurries or other cementitious fluids.
Since the 1970s many experimental studies have been published on the displacement process for primary cementing applications [5,6,7,8,9,10,11,12,13]. Some of these were more academic, using model fluids in small scale setup, while others have been more industry oriented, using curing cement in full scale flow rigs, while yet other researchers have studied displacement of drilling fluids with real cement slurries inside metal casings [14,15].
The experimental equipment presented here is designed for displacement measurements with non-hardening laboratory fluids, which is an applicable way of evaluating displacements experimentally [2,4,16,17,18]. The two main approaches for monitoring the motion of the interface between the displaced and displacing fluids have been optical/visual and electrical resistance. The first approach is very useful for small scale experiments, providing nearly complete information on the interface position as function of time. However, it requires good optical contrast between the fluids and carefully designed lighting to minimize reflections. Also, differences in refractive index of materials should be avoided or accounted for. This becomes expensive for larger scale experiments and is not possible when using field fluids or other opaque fluids. Electrical resistance monitoring has proved to be a fairly simple and reliable method of detecting/discriminating between fluids when one of them can be made to be more conductive by addition of salts. This method was used in the works by Deawwanich [8], Tehrani et al. [6] and by Lockyear et al. [7].

2. Experimental Design

2.1. General Fluid Design

The test setup is constructed for evaluating different displacement processes in well operations. To be able to run a large number of tests with repetitions, it was necessary to use non-hardening fluids with viscous properties similar to drilling fluids, spacers and cements. The design requirements for the selected fluids to be used in the displacement experiments were:
(1)
viscous properties reflecting drilling fluid, spacer and cement slurry viscosities
(2)
viscosity contrast reflecting those of field applications
(3)
density contrast to simulate density differences in field operations
(4)
electrical conductivity contrast to be able to measure the displacement efficiency
(5)
optical transparency to be able to measure the displacement efficiency
(6)
low toxicity to ensure no harm for involved personnel and reduce cost of laboratory waste handling
(7)
low cost
The displaced and the displacing fluids were designed to have viscous properties similar to those expected in a real cementing operation, i.e., shear-thinning also with a yield stress. In addition, the displacing fluid was designed to be heavier and more viscous than the displaced fluid. To facilitate visual observation, the displaced fluid was designed to be as transparent as possible. Two different fluid systems were used, as shown in Table 1.
Barite was used as weight material of the displacing fluid in fluid system 1. However, table sugar (sucrose) was found to suit the test system better and was used to increase the density in fluid system 2. Sucrose dissolves in water to give the increased density without the addition of solid particles. Also, as a minor point, sucrose will increase the index of refraction in the fluids, which is also advantageous for visualization. Colouring agents (uranine and methylene blue) were added to the displacing fluid in fluid system 2 to improve visual contrast between the fluids. Finally, NaCl was added to the displacing fluid to provide electrical conductivity contrast.
A study of the fluid mechanical basis for conducting laboratory experiments with the scope of presenting displacement results that are relevant for industrial application concluded that it can be possible to use laboratory fluids in laboratory test equipment and produce results that can be scaled up to field applications [19]. Occasionally, the casings are rotated in well cementing operations to improve the displacement. The rotation rates in the laboratory will in general be different than those in the field as a result of using different viscous properties along with operating in a different geometry.

2.2. Test Rig

2.2.1. Overview

The experiments were conducted at ambient conditions using a purpose-built flow rig representing a wellbore. The test section was ca 10 m long and mounted on an inclinable frame, as shown in Figure 1. It is constructed as an annulus with a constant diameter inner pipe, representing the casing or liner. A process and instrumentation diagram is shown in Figure 2.
The second to the last fifth length of the outer pipe was replaced by a larger diameter pipe, simulating a wash out section as indicated by the colored section in Figure 2.
Transparent acrylic pipe sections of 2 m length were used both to represent the wellbore and the casing inside the wellbore. The diameter size of the gauge annular section was 5” by 6.5” (0.127 m by 0.1651 m). The wellbore of the washout section had a larger diameter and was concentric with respect to the rest of the wellbore. Figure 3 shows the details of the washout section. Note that the diameter at both ends of the washout section was tapered with an angle of 50° relative to wellbore axis from 11” (0.2794 m) to 6.5” (0.1651 m). Key data for the test section is listed in Table 2.
The inner pipe was supported by bearing joints (see Figure 4) every 2 m at the wellbore section connections to prevent inner pipe sagging. The bearing joint consists of an outer flange supporting the outer pipe, an inner pipe bearing set to support the inner pipe and to transmit torque between sections of inner pipe, and bolts to fix the inner pipe bearing laterally and to allow adjustment of eccentricity. The bolts each have a diameter of 10 mm. This represents about 4.4% of the cross-section area of the annulus and is thus expected to have a minimal impact on the flow and displacement.
The inner pipe was also mounted to a torque motor at the downstream end, allowing rotation of the inner pipe at a controlled rate. The eccentricity ε of the inner pipe, defined as the offset of the inner pipe axis relative to the wellbore axis divided by the concentric annular gap, could be adjusted manually using screws at the joints. Experiments were conducted with eccentricities of ε = 0 and ε = −0.42, with negative values for inner pipe axis below the wellbore axis.
The displaced and displacing fluid were stored in separate tanks prior to the experiments (see Figure 2). Within each tank there was a thermostat-controlled heating element to keep the temperature of the fluids at a temperature slightly above ambient temperature. Two valves at the flow loop inlet, see Figure 5, were used to select the source of flow through the test section being either displaced or displacing fluid.
In Figure 5 it is seen that the fluid makes a 90 degree turn at the inlet. This does not provide optimal inlet flow conditions, in particular for fluids with elasticity. We did not measure the elastic properties of the fluids used in this setup. However, we have not noticed any clear indications that the inlet bend has affected the results.

2.2.2. Instrumentation

The flow rate was measured using am Optiflux 4000 electromagnetic flow meter (Krohne, Duisburg, Germany) that was mounted downstream of the main pump Q2 as shown in Figure 2. To measure the pressure losses along the annulus, three differential pressure (DP) cells of type Fuji were used, see Table 3. These measured the differential pressure across three sequential 2 m sections of the test section. The pressure cell DP1 was used for measuring the gauge hole annular pressure losses. The other two DP cells had one or two pressure outlets mounted near the washout section and the corresponding pressure measurements were affected by the diameter changes at the ends of the washout section. It turned out that DP3 was faulty and no data are available from this sensor.
The location of the pressure ports were azimuthally between probes 6 and 7, and axially at the same position as the conductivity probes, with two of the ports shared between neighboring DP cells.
The propagation of the displacing fluid was monitored using conductivity probes, each consisting of a pair of electrodes. For each probe the conductivity between the electrodes was measured. 24 conductivity probes were located at four different axial stations along the test section (stations A, B, C and D in Figure 2). Six of eight slots for conductivity probes were used at each station; slots 2 and 4 were installed but not used in measurements (Figure 6). The number of actual probes was limited by the number of available channels of the amplifier. Note that no measurement, except visual, could be conducted in the wash out section. Here the measurement probes were positioned immediately upstream and downstream this section.
Each conductivity probe consisted of a pair of 3 mm diameter steel screws which were directed radially into the annulus and positioned 8 mm apart (center-center distance). Thus, the gap between the electrodes was 5 mm. In order to reduce any effects of the flow disturbance caused by the electrodes on the measurements, the electrodes in each pair were mounted at the same axial position, i.e., separated in the azimuthal direction. The penetration depth was 5 mm, while the width of the concentric annulus was ca 19 mm. Thus, the electrodes covered about one quarter of the radial distance in concentric annulus configuration.
The conductivity probes were connected to a signal amplifier which was specially designed for use with conductivity probes for measuring water wave heights (Wave Amplifier type 108, DHI, Hørsholm, Denmark). This amplifier gives an electrical AC signal around 2750 Hz. The resistance λ per unit length between a pair of infinitely long parallel cylindrical electrodes of common diameter δ, is given as shown in Equation (1) [20]:
λ = ψ π cosh 1 ( Δ δ )
where ψ is the resistivity of the fluid around the electrodes and Δ is their center-center distance. The geometrical factor of Equation (1) thus changes by 8% with an increase in distance Δ from 7.5 mm to 8.5 mm. Since the electrodes do not span the entire annular gap, a key question is to what extent they will detect conducting fluid which is not in direct contact with the electrodes. The measurement principle is illustrated in Figure 7. As the displacement process progresses, the contact with the high conductivity fluid increases and a general increase in conductivity is measured. In this case, which is representative for the displacement experiments with concentric annulus, the displaced fluid is stratified above the displacing fluid due to lower density. The interface is here strongly inclined with respect to the cross-section plane, and almost parallel to the wellbore axis. Thus the length of electrode contacting the conducting fluid increases with time as the thickness of the displaced fluid layer diminished. This gives a fairly slow transition from minimum to maximum conductivity. When the interface is nearly normal to the wellbore axis, and thus parallel to the electrode axes, the transition from minimum to maximum conductivity will be much more abrupt, unless there is strong mixing between the fluids.
Except for an offset value, the measured voltage signal, Vo, is proportional to the average conductance G between the electrodes. Assuming that boundary effects can be ignored, the conductance G is proportional to the length Lp of the electrodes of a probe, i.e.:
G = L p λ = π σ L p cosh 1 ( Δ δ )
where σ ≡ 1/ψ is the fluid conductivity. The conductivity increases with increasing length of probe being immersed in the high conductive fluid.
The wave amplifier ensures that the voltage across the electrodes is constant; being independent of the conductance between the electrodes. Instead, a change in the conductance will result in a change in the amplification in the wave amplifier, and thus a change in the output voltage. Thus, accounting for any voltage offset we may write the measured voltage signal as given in Equation (3):
V o = V o f f s e t + K   G
where K is a gain factor of the measurement system.
Both the displaced and the displacing fluids have some conductivity. However, the conductivity of the displaced fluid is an order of magnitude smaller than that of the displacing fluid, and we may approximate it as being non-conducting. The conductivity probe signal response will partly be due to changes in the local conductivity, σ, due to fluid mixing, and partly due to the position of the interface between low and high conductivity fluid relative to the probes.
The measured voltage for each pair of electrodes will vary between a minimum value Vo,min corresponding to the displaced fluid with conductivity Gmin and a maximum value Vo,max corresponding to the displacing fluid with conductivity Gmax. Assuming linearity between conductivity and volume fraction of displacing (high-conductivity) fluid we have:
G = ϕ G max + ( 1 ϕ ) G min
Using Equation (3) we then obtain:
ϕ = G G min G max G min = V o V o , min V o , max V o , min V r
Equation (5) holds true for a homogeneous mixture of the two fluids, and is also expected to hold true in some average sense for a heterogeneous system where there is a macroscopic interface between the two fluids.
Although a given value of ϕ can be obtained from this expression both from a homogeneous mixture of two fluids and from a given configuration of a sharp interface between the same fluids, we may distinguish between these two cases by considering the time derivative of the signal Vo(t). A section with mixing will yield a lesser pronounced change of the signal with time than a sharp displacement front.
Let S be the total fluid domain around a pair of electrodes. Let further ϕ(r) be the local volume fraction of displacing fluid and define an interface ∂Si between displaced and displacing fluid by ϕ(r) = 0.5. Although the fluids in the experiments were miscible (both were water-based), we assume that the interface was fairly well defined through each displacement experiment, since the duration of each experiment (on the order of one minute) is short relative to the time for molecular diffusion over a typical length scale (annular gap width). The fluid-fluid interface may become convoluted due to convective processes driven by gravity or hydrodynamic instabilities, but microscopic mixing of the fluids before the washout is expected to be primarily due to turbulence in the entrance region to the test section.
Since Maxwell’s equations are linear, we postulate that the normalized signal response Vr will not depend on the absolute values of the conductivities. We also postulate that it will depend only weakly on the shape of the wall interface ∂Sw of ∂S = ∂Sw + ∂Si, where ∂S is the total boundary of the fluid domain. We here ignore the inlets and outlets of the test section, which are both far removed from the conductivity probe electrodes. Thus, we conclude that Vr almost exclusively depends on the distribution ϕ(r) and, assuming there is little mixing between the fluids, mainly on the interface ∂Si. Note that this interface is not necessarily simply connected, a well-defined fluid front. For example, there may be isolated pockets of bypassed displaced fluid. If, however, the interface ∂Si is a well-defined fluid front, separating the displaced and displacing fluids axially, we expect that the signal response Vr will be a monotonous function of the axial position of this displacement front.
An additional conductivity probe, type CTI-500 (Jumo, Fulda, Germany, hereafter called C0), was mounted on a spool section in the fluid outlet downstream of the test section. The spool was connected to the test section with a 2.5” flexible hose of length approximately 9.5 m. The probe is a combined inductive conductivity and temperature transmitter and measured the electrical conductivity and the temperature of the fluid which had exited the test section. From the measured conductivity the volume fraction of the conducting fluid and thus the overall efficiency of the displacement process can be determined. There was no special mixing element upstream of the probe C0. However, it is assumed that the two fluids were well mixed in the transversal plane at the location of the C0, both because of the length of the connecting flexible hose, the higher fluid velocity in the hose (1.38 m/s), and because of the change in flow direction at the outlet of the test section.
The objectives of the measurements with the conductivity probes were three-fold:
(1)
detect the arrival of the displacing fluid at a particular probe location at the corresponding time where the signal Vr changes from low to high
(2)
assess qualitatively the degree of mixing of the two fluids at the probe location by measuring the duration of the transition from low to high signal level
(3)
measure and quantify the concentration of displacing fluid at the probe location
For the two first objectives it is only necessary to assume that the measured signal is a monotonously increasing function of the relative fraction ϕ of displacing fluid and that the two fluids have sufficiently large conductivity contrast. These were the main objectives of the distributed conductivity probes. The third objective requires a known relationship between the measured signal and the volume concentration of the displacing fluid and that the concentration is uniform in the sampling volume of the probe. In the present experiments the concentration of electrolyte, NaCl, in the displacing fluid was so low that a linear relationship can be assumed, as discussed below.
Since the scope of the measurements is to determine the displacement fluid volume fraction ϕ rather than the effective conductance G between the electrodes, we will in the following work only with the normalized probe response Vr. Since the absolute value V0 of each probe was not critical, the probes were not calibrated to give the same value for the same fluid (displaced and displacing). The offset and gain factors for each channel of the signal amplifier were adjusted prior to conducting the experiments to avoid saturation of the signal; meaning being out of range. Still this has happened in a few cases.
The linearity of ϕ(G) (see Equation (5)) was checked before the experiments by measuring the conductivity of a mixture of a non-conductive and a conductive water-based fluid. Thus, two prototype fluids A0 (water with 1.4% by weight of the (synthetic) clay laponite) and B0 (water with 0.8% by weight xanthan biopolymer and 1.0% by weight NaCl) were mixed in different ratios and the conductivity was measured using a handheld conductivity probe (Hanna Instruments HI933000). The results are shown in Figure 8. The measurement points indicate a linear relationship between conductivity and volume fraction, as shown by the stippled line. For such a mixture the conductivity is proportional to the volume concentration of the electrolytic ions when the concentration is sufficiently low. Note that the fluids A0 and B0 were not used in the actual displacement experiments, but represent fluid system 1, except for the omission of barite, soda ash and biocide in the displacing fluid. This omission is not expected to have any impact on the linearity of the conductivity versus composition.
While conducting the experiments, the following data were logged 50 times per second:
  • Flow rate (FM)
  • Differential pressure (DP1, DP2, DP3)
  • Inner pipe rotation speed (RPM)
  • Mixed fluid conductivity downstream test section (probe C0)
  • Mixed fluid temperature downstream test section (probe T0)
  • Tank temperature of displacing fluid (probe T1)
  • Local conductivity at stations A, B, C and D (probes Aj, Bj, Cj, Dj with j = 1, 3, 5, 6, 7, 8)
The displacement process was also set up to be monitored visually using six low-cost web-cameras mounted along the test section as indicated in Table 4.
Fluids were also characterized in terms of electrical conductivity using a HI933000 handheld conductivity probe (Hanna Instruments, Woonsocket, RI, USA). In addition, the pH was frequently measured. To keep the volume of test fluids at a minimum, as much as possible was reused. Therefore, there would always be a little degree of contamination of the fluids with the other fluid. As seen in Figure 9, the conductivity of the displaced fluid increased over time, while the conductivity of the displacing fluid was more stable. Still the contrast is sufficiently large to be able to measure the displacement properly.

2.3. Test Procedures

A detailed test procedure was established for running the experiments. The physical items in the test procedure are illustrated in Figure 3. The test procedure is as follows:
(1)
Prepare data logging program:
(a)
check/set appropriate limit values for conductivity probe C0 (see Figure 8). These values control the outlet valves.
(b)
check/set correct flow rate value (0.00436 m3/s, corresponding to 0.5 m/s in the regular size annulus)).
(2)
Start data logging
(3)
Start circulation of fluids with displaced fluid circulating in the test section and displacing fluid circulating in the bypass line.
(4)
Take sample of displaced fluid.
(5)
Flush DP cells with displaced fluid.
(6)
Stop circulation of displaced fluid for zero DP measurements.
(7)
Restart circulation of displaced fluid.
(8)
Turn on video lights and start video recording
(9)
Start displacement experiment by switching to displacing fluid using computer-controlled valve
(10)
Switch to automatic outlet valve control mode in logging program
(11)
Stop video recordings when outlet valves have switched to position 3 (displacing fluid)
(12)
Continue circulation with displacing fluid for at least one minute after stopping experiment.
(13)
Take sample of displacing fluid
(14)
Stop circulation

3. Data evaluation and Processing

3.1. Pressure Drop Measurements

The pressure gradient values were logged as indicated in Figure 2. The viscosity flow curves were measured using a manual six speed Fann35 viscometer and the fluid viscosity is described reasonably well using the Herschel-Bulkley model. The model parameter values for the fluid are shown in Table 5 together with calculated values for Reynolds number and pressure gradient for concentric and eccentric casing. In the table the surplus stress τs is a Herschel-Bulkley parameter following Saasen & Ytrehus [21], and is defined as:
τ s = τ ( γ ˙ s ) τ y
where γ ˙ s is a shear rate which is characteristic for the flow. τ ( γ ˙ s ) is interpolated from the Fann measurements. The characteristic shear rate, γ ˙ s , is defined as
γ ˙ s = 12 U D o D i = 158   s 1
where U is the average fluid velocity in the regular section. The surplus stress is presented to visualize the viscosity difference at the shear rates relevant for the fluid flow.

3.2. Pump Rate Variability

The delivered power to the pump was adjusted automatically in order to maintain constant flowrate as the fluid source is switched from displaced fluid to displacing fluid. Initially there were unacceptably large fluctuations in the flow rate adjustments. Before finishing the development of the set-up, the flow rate control was improved, leading to much smaller fluctuations. Figure 10 shows the normalized flow rates versus times for the experiments after the improved algorithm was implemented.

3.3. Conductivity Probe Signals

In order to measure the concentration of the displacing fluid and, hence, measure the displacement efficiency it is necessary to first normalise the probe signals. Thereafter, the normalised signals are converted to concentration values of the displacing fluid. In this way the local measurements of the displacing fluid volume fraction are obtained. Sample time traces are shown in Figure 11.
By comparing the measurements of neighbouring measurement points, both axially and around the perimeter, the shape and kinematics of the interface can be reconstructed. As an example we show measurements from two experiments, named N and R using same naming convention as in [22]. Both experiments are with concentric annulus at inclination 60 degrees from vertical, where experiment R is a repetition of experiment N to verify repeatability. These measurements are shown in Figure 11 for data recorded at the symmetric (with respect to a vertical plane containing the wellbore axis) probes A1 and A8 for two experiments. Probes A1 and A8 are at the axial position of station A, being the first measurement station after the test section inlet. The circumferential positions P1 and P8 are illustrated in Figure 6. Since there is no rotation of the inner pipe in these experiments, the signals of A1 and A8 should be equal by symmetry, and a comparison of these also provides a quality check of the data. It is noticed that the symmetry is quite well satisfied, even at later times. Reproducibility is also quite good. However, at later times there is a discrepancy as the displacement in the original experiment seems to be somewhat better than the repetition results. This is anticipated to be caused by the unintentional mixing of the fluids as the fluids are reused.
As discussed above, the conductivity probe signals are assumed to be proportional to the local fraction of displacing fluid at the location of the probe. Although these are only point measurements, they provide an extensive amount of information on the displacement process. We illustrate this by the experimental results shown in Figure 14 where two experiments with horizontal, concentric annulus, without inner pipe rotation are compared, as they represent the simplest experimental conditions. Also, the comparison is interesting, because the viscous properties of the fluids used in the two experiments (A and O) were fairly similar, both in terms of observed pressure drop values and in terms of viscosity measurements. Thus, we expect the results to be similar. In Figure 12 the local fraction of displacing fluid at station A, 2 m from the inlet, is plotted as function of the normalized pump volume (see Equation (13)), thus correcting for any variations in flow rate.
The following information can be obtained from measurement plots like the one shown in Figure 11:
(1)
The time at which signal value starts increasing indicates when displacing fluid arrives at the different locations (arrival time).
(2)
The difference in arrival time of signals at different axial positions (stations) but at the same azimuth position can be used to estimate the axial flow velocity of the interface.
(3)
The steepness of the curves indicates the degree of mixing of displaced and displacing fluid at the location of the probe. A steep transition indicates a well-defined interface with little mixing. A very gradual interface could be due to mixing of the two fluids, but it could also be caused by an interface which is slanted with respect to the axial flow direction.
(4)
Similarity in shape of probe signals at subsequent axial stations (A, B, C) in the regular section indicates that the displacement front propagates without deformation like a travelling wave. In practice there will always be some distortion, which can be inferred from the differences in shape among the probes.
It is noticed for Experiment O in Figure 12 that in the lower half of the annulus (probes A3, A5 and A6) the displacing fluid arrives sooner than a piston displacement, most likely as an effect of the larger density of the displacing fluid. The vertical line labelled piston displacement corresponds to the time when the total volume of displacing fluid in the test section is equal to the test section volume upstream of station A. For the measurement at the top of the annulus near the top of the annulus (probes A1 and A8) the signals are close to synchronous At the piston displacement time the lower part of the annulus contains about 70% displacing fluid in the O experiment, but only 10–20% displacing fluid in the A experiment. Still, in both experiments the displacing fluid front is seen to be inclined with respect to the pipe axis, since displacing fluid arrives earlier near the bottom of the annulus than near the top.
Video cameras covered the area between the measurement probe arrays and provided valuable supplement to the conductance measurements. Simple videos are more challenging to use for detailed analysis since the interface becomes difficult to track. Particles designed to remain on the fluid-fluid interface could enhance the value of videos. Analysis of the videos alone is difficult to use because of increasing opacity of the fluids, positioning of illumination and reflection of the illuminated light on the pipe wall. Still this method was applicable as additional information to the conductance probe data.

3.4. Calculation of Displacement Efficiency

We define the displacement efficiency following Deawwanich [8]:
χ ( t ) = Q ( t t 0 ) V ;   t t b χ ( t ) = Q ( t b t 0 ) + t b t d t [ c A ( t ) ] V ;   t > t b
where t0 is the time when the displacing fluid enters the test section, V is the annular volume of the test section, and tb is the time when the first displacing fluid leaves the test section. After the breakthrough time tb the concentration cA(t) of displaced fluid at the outlet decreases and the displacement efficiency increases at a slower rate. As time goes to infinity the displacement efficiency theoretically goes to unity provided no immobilized displaced fluid is left behind. An implicit assumption is that the displaced and displacing fluid move with the same velocity, so that the flow rate of the displaced fluid at any location can be written as:
Q A ( t ) = Q ( t ) c A ( t )
In our experiments we did not log the starting time t0 directly. However, we define the starting time as the time when the inlet valves were switched from pumping displaced fluid to pumping displacing fluid. This time t0 was inferred from the measurements of flow rate and differential pressure in the test section as described in the following. When the inlet valves are switched, a pressure pulse travels through the test section and there was a transient in the flow rate through the test section. This transient was due to capacity limitations of the pump for the displaced fluid. Before the switching of the inlet valves the flow in the test section was driven by the pump for the displaced fluid at full capacity, which corresponded to a rate somewhat lower than the scheduled rate of the experiment. Thus, at the start of the experiment, corresponding to the switching of the inlet valves the fluid in the test section was accelerated to the target flow rate. In Figure 10 this is seen as a transient increase in the flowrate as the control of the pump power became over-compensated. Since the flow rate was not constant, the expression for the displacement efficiency was calculated as:
χ ( t ) = t b t d t Q ( t ) [ c A ( t ) ] V T S
for times t > t0, where VTS is the fluid volume between the inlet point and the conductivity probe C0 at the outlet and cA is the concentration of displaced fluid at the outlet. We assume, based on the results shown in Figure 8, that the volumetric concentration cA(t) of the displaced fluid at the outlet is a linear function of the global conductivity σ0 similarly to the predictions given by Equation (11).
c A ( t ) = σ 0 * ( t ) σ 0 B σ 0 A σ 0 B
where the conductivities of the displaced and displacing fluids are σ0A and σ0B, and we limit the concentration to be a number between zero and one:
σ 0 * ( t ) = min ( σ 0 B , max ( σ 0 ( t ) , σ 0 A ) )
These conductivities σ0A and σ0B are calculated by time averages over suitable time intervals before and after the displacement process.
To be applicable for comparison with field operation data it is necessary to plot the displacement efficiency versus normalized pumped volume rather than versus time. Thus, we define the normalized pumped volume v(t) as independent variable:
v ( t ) = t 0 t d t Q ( t ) V T S
We also calculate the breakthrough time, i.e., the time of arrival of first displacing fluid at the outlet. Note, however, that the outlet is defined to be at the location of the conductivity probe, C0, which is located on a flexible hose about 9.5 m downstream of the annular test section outlet. This distance corresponds to a volume of ca 30 L with a transit time dependent on the nominal flow rate. Thus, to calculate the breakthrough time at the outlet of the annular section, a flow dependent time value should be subtracted from the measured breakthrough time. A corresponding adjustment should then be done for the relative breakthrough volume. For some experiments this would give negative values which could be caused by noise in the signal from the global conductivity probe. However, we are here primarily interested in a relative comparison among the experiments, and not in the absolute values of the breakthrough times.

3.5. Reconstruction of Displacement Front Velocity

The axial velocities of the displacement front are reconstructed based on the signals from conductivity probes along the test section and this was subsequently used to reconstruct the total fluid flux. Since the latter is also measured, the flux was also used to renormalize the calculated axial velocities. The calculations are based on a cross correlation of signals from probes at neighbouring measurement stations (A, B, C, D) but at the same circumferential position. This approach differs from the one used by Renteria et al. [2], who based the calculation on the time t when the signal level reached 20% of the maximum signal level. The motivation for the present approach was to avoid the use of an ad hoc signal threshold level.
Using the calculated local velocities, we estimate the total liquid flow rate by numerical integration over the annular cross section. Since the total liquid flow rate is measured using a Coriolis flow meter, the accuracy of the calculated velocities can be calculated. This accuracy is generally very good.
The cumulative flux of displaced fluid at the different measurement stations are estimated by assuming that the local axial velocities at each station are time-independent and by assuming that the average local concentration of displacing fluid is proportional to the normalized signal level. For most experiments and most measurement locations, this calculation gives displacement efficiencies being far in excess of unity. This indicates that there is a significant velocity difference between the displaced and the displacing fluid. Such a difference can be explained by assuming that the displacing fluid advances more rapidly in the centre of the annular channel than closer to the wall. Since the probes only protrude 5 mm from the outer wall into the annulus, the signal will be biased towards the fluid present at the wall. Although not directly confirmed by experimental observation, the assumption mentioned above is supported theoretically by the no-slip velocity condition at the wall.
First, we estimate the local axial velocity vi,j at circumferential position i and station j (j = 1, 2, 3, 4 for stations A, B, C, D respectively). Thus:
v i , j = Δ L j Δ T i , j
where ΔLj has been measured to be 1962 mm while the time delay ΔTi,j is calculated from the cross-correlation of normalized signals Ci,j and Ci,j-1.
We define a cross-correlation function S as
S i , j ( k ) = t min t max [ C i , j ( t ) C i , j 1 ( t k Δ t ) ] 2 d t
and:
Δ T i , j = k * Δ t ;   S i , j ( k * ) S i , j ( k )   k k *
Then the total axial volume flux Φj of fluid at any axial station j is estimated for the cases without inner pipe rotation as:
Φ j = 2 r i = 5 , 6 , 7 , 8 v i , j A i , j = Q
using symmetry (the factor 2), where Ai,j is the annular flow area associated with probe i at station j. The factor r is a renormalization factor which is determined by setting the flux equal to the measured flow rate. Hence a renormalized velocity is defined as
v i , j r = r v i , j = Q 2 k = 5 , 6 , 7 , 8 v k , j A k , j v i , j
In some cases, Equation (14) gives an unphysically large velocity, larger than a set value vmax. Then the value is discarded, and instead an area weighted average over the other probe positions is used.
We assume that the axial liquid velocity at any circumferential position is time-independent, and estimate the axial flux of displaced fluid at any station as
Φ j , d ( t ) = 2 i = 5 , 6 , 7 , 8 v i , j A i , j [ 1 C i , j ( t ) ]
The areas Ai,j are all equal to Aa/8 for a concentric annulus. For an eccentric annulus, the areas depend on the circumferential position i, but not on the station j.
In order to calculate the areas Ai,j for an eccentric annulus we introduce the local gap H(θ) depending on the azimuthal angle θ. As shown in Figure 13, θ is the angle between two radius lines from the centre of the inner cylinder (casing) where one line connects the centre with a point P on the perimeter of the outer cylinder wellbore wall) and the other line is vertical upwards from the centre. By trigonometry we have
( R i + H ( θ ) ) 2 + δ 2 2 ( R i + H ( θ ) ) δ cos θ = R o 2
where:
δ = ε ( R o R i )
Thus:
H ( θ ) = R o 2 δ 2 sin 2 θ R i + δ cos θ
and:
A i , j = θ j Δ θ / 2 θ j + Δ θ / 2 R i R i + H ( θ ) r d r d θ
where Δθ = π/4.
For example:
A i , 1 = A i , 8 = 0 π / 4 R i R i + H ( θ ) r d r d θ
It is convenient to normalize by the cross-sectional area of the casing, as shown in Equation (25). Thus:
a i , j A i , j π R i 2 = 1 π θ j Δ θ / 2 θ j + Δ θ / 2 1 1 + h ( θ ) x d x d θ
where:
h ( θ ) = r o 2 ε 2 ( r o 1 ) 2 sin 2 θ 1 + ε ( r o 1 ) cos θ
where ro = Ro/Ri.
Let Vj be the fluid volume between inlet valves and station j (j = A, B, C, D), and define displacement efficiency χj for the respective stations as:
χ j ( t ) = 1 V j 0 t Φ j , d ( t ) d t
χ j ( t ) = 2 k = 5 , 6 , 7 , 8 v k , j A k , j 0 t [ 1 C k , j ( t ) ] d t V j
We also estimate the displacement efficiency for the sections between stations j and j+1 and, in particular, the displacement efficiency for the washout section between stations C and D as:
χ j , j + 1 ( t ) = 0 t [ Φ j + 1 , d ( t ) Φ j , d ( t ) ] d t V j + 1 V j = 0 t [ V j + 1 χ j + 1 ( t ) V j χ j ( t ) ] d t V j + 1 V j

4. Sample Measurements

Figure 14 shows time traces from experiment in concentric annulus with 60 degrees well inclination. Average pressure gradients between 2 m and 4 m from inlet (DP1) and between 4 m and 6 m from inlet (DP2) are shown in black and red, respectively. The displacement process starts at time t = 0. Each experiment includes a sequence without flow in the test section before starting the displacement (here between-180 and-150 s), to obtain a zero reference for the pressure measurement.
Note that before the start of the displacement process the displacing fluid is recirculated with the target flow rate in a bypass loop, while the displaced fluid is circulated through the test section (see Figure 2). The flow through the test section is briefly stopped before the start of the displacement (here between −180 and −150 s) to obtain a zero reference for the DP cells. The flow meter (see symbol FM in Figure 14) is, however, mounted so that it always measures the flow rate of the displacing fluid, and the logged flow rate FM is therefore constant throughout the experiment. The flow rate of the displaced fluid before the start of the displacement, which corresponds to the pressure measurement before t = 0, was not logged, but was somewhat smaller than the flow rate FM of the displacement experiment. A small transient in FM is therefore observed as the displacement process starts due to the inertia of the fluid. In this experiment the flow rate briefly increases at t = 0 due to a small overcompensation of the pump control.
At t = 0 the measured pressure gradient increases mainly because the flow rate increases, but also because of the change in viscous properties of the fluid as the displacement proceeds.
The yellow curve (C0) shows absolute fluid conductivity at the sensor at 9.5 m downstream of the test section outlet, and the blue curve shows the corresponding temperature T0. The small change in temperature is well correlated with the change in conductivity, representing the arrival of the displacing fluid at this location. The pink (C0_corr) and green (TempC0_corr) curves show the same data corrected for the time offset due to the transit time from the test section to the sensor.
Figure 15 shows pressure gradients DP1 for the experiment shown in Figure 14 and the pressure gradient of a repetition test. The repetition test was conducted to verify the reproducibility of the experiments. In the repetition experiment a zero flow reference was not conducted, so the pressure reading before start of the displacement was selected as reference, It is seen that the reproducibility in terms of the pressure gradient during the displacement process is very good.
From the differential pressure measurements it is possible to calculate pressure gradients in the steady state situation (i.e., latter parts of the time traces such as shown in Figure 15). The frictional part of the pressure gradient has been obtained in the usual way by accounting for the difference between the process liquid density and the measurement liquid density (i.e., of the liquid in the measurement lines), according to:
| d p f d x | = | d p d x | Δ ρ g L cos θ
The measured pressure gradient dp/dx and the corresponding frictional pressure gradient dpf/dx from sensor DP1, as calculated from Equation (30) have been plotted in Figure 16, assuming a density difference of Δρ = 100 kg/m3.
The frictional pressure gradient should not depend on inclination, assuming that the annulus between the pressure ports is filled with displacing fluid only during the measurements. However, we notice that both for concentric and eccentric annulus the correction for the hydrostatic head gives an overcompensation. A similar over-compensation is seen for the sensor DP2. This indicates that the over-compensation could, at least partly, be due to a deviation in the assumed density difference. The change in measured pressure gradient with inclination corresponds to approximately Δρ = 50 kg/m3 for both eccentric and concentric annulus. Sensor DP2 consistently gave lower pressure gradients than sensor DP1. This can be explained by the proximity of the downstream port of sensor DP2 to the washout section, disturbing that measurement. The measured pressure gradient for sensor DP1 for horizontal annulus, agreed well with calculations both for concentric and eccentric annulus, as shown in Figure 16.
As mentioned in Section 2, the test section was equipped with an electrical motor which allowed the inner pipe to be rotated about its own axis, simulating cementing with rotating casing. Sample results are shown in Figure 17, demonstrating the significant impact of inner pipe rotation on displacement in an eccentric annulus. Without rotation (solid lines) the displacement process is severely delayed in the narrow annulus near the bottom (red, solid curve). With rotation (dashed lines), the fluid in the narrow part of the annulus (red, dashed curve) is displaced at the same time as or even before displacement occurs in the upper (wide) part of the annulus (green and black curves).
Figure 17 also shows how the rotation breaks the symmetry in the response of probes 1 and 8 (these probes are located symmetrically with respect to the vertical diameter). Video cameras mounted along the test section provided additional qualitative information on the nature of the displacement. To illustrate this, we have included in Figure 18 two images from video recordings of the washout section, illustrating the effect of eccentricity on the displacement in the washout section. This is of particular interest since there were no conductivity probes mounted inside the washout section. The images are taken from the side at the center part of the washout section (between 6 and 8 m from the inlet), covering the entire wellbore diameter, with flow from right to left. The left and right images show displacement with concentric and eccentric horizontal annulus, respectively. In the image to the left the interface is almost flat, and the corresponding video shows how this interface moves vertically upwards as the washout section is being filled with displacing fluid. In the image of the eccentric annulus the interface is strongly inclined, and the corresponding video shows how the displacing fluid fills the washout like a flooding wave which gradually becomes more horizontal. The volumetric flow rate of the experiments is 4.37 l/s while the volume of the washout section is approximately 97 L, corresponding to a minimum filling (or displacement) time of 22 s for the washout section. The videos indicate that the actual filling time is longer than this, which is also confirmed by analysis of the conductivity probe signals, see Figure 19. The actual filling time can be estimated as the delay between signal D1 or D8 relative to signal C5, which is about 70 s (comparing the time traces of probes C5 and D1).

5. Conclusions

A measurement system for liquid-liquid displacement experiments in concentric and eccentric annuli has been presented. Electrical conductance probes provide quantitative point measurements of the displacing fluid fraction at different axial and azimuthal positions along the wellbore, while video recordings provided supplementary qualitative information. The probes worked well, as demonstrated by overall consistency and repeatability of the signals. Albeit their individual response from point locations in the same vicinity from the wall, the combination of conductance probes gives an excellent measurement of the local average displacement process including mixing of miscible fluids. This electrical conductance technique has also been used successfully by other researchers [6,7,8] to generate displacement data and it should be accepted as a reliable method for measurement of displacement efficiency.
The video cameras used here covered the area between the arrays. Even with simple web-cameras as used in these experiments, video recordings can provide valuable supplement to and corroborate the conductance measurements. In the present experiments the analysis of the videos was difficult due to increasing opacity of the fluids, positioning of illumination and reflection of the illuminated light on the pipe wall. The recordings still provided qualitatively useful information. That was illustrated above by comparing the displacement pattern in the washout section for concentric and eccentric annuli, respectively.
The measured pressure gradient for sensor DP1 for horizontal annulus, agreed well with calculations both for concentric and eccentric annulus. A main motivation of the pressure measurements was to have a fluid quality check. As the calculations were based on viscometric measurements, the results demonstrated that the pressure measurements were consistent with the viscosimeter measurements.
Finally, the suitability of the device that was applied in the displacement studies presented by Renteria et al. [2], Tardy [4], Ytrehus et al. [16], Lund et al. [17,22] and Skadsem et al. [18] has been verified. Note that some of these studies have been used in development of displacement software.

Author Contributions

Conceptualization, J.D.Y. and A.S.; methodology, B.L. and A.T.; formal analysis, B.L.; investigation, A.T.; writing—original draft preparation, B.L.; writing—review and editing, B.L., J.D.Y., A.T. and A.S.; visualization, B.L. and A.T.; project administration, J.D.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded through the Research Centre DrillWell, supported by the Research Council of Norway, Aker BP, ConocoPhillips, Equinor and Wintershall. This publication was also financed by the research grant 294688, Optimized Wellbore Hydraulics.

Acknowledgments

The authors would like to thank the Research Council of Norway, Aker BP, ConocoPhillips, Statoil (now Equinor) and Wintershall, and for funding this work through the research centre DrillWell. The authors would also like to thank MI-SWACO, Schlumberger for providing the chemicals used for the fluid design and the technical support by Knud Richard Gyland. The authors also thanks Gunnar Lia Giil, Espen Johansen and Birgitte Ruud Kosberg for their assistance with the fluid design and control during their master program at NTNU, and PhD student Shreyansh Divyankar for assisting with displacement flow experiments. Last, but not least, the support with instrumentation and test preparation provided by our colleagues Terje Øyangen and Arne Erik Rekkebo is highly appreciated.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

SymbolDescription(Default) value
DiInner diameter of annulus (outer diameter of inner pipe)0.127 m (5”)
DoInner diameter of gauge outer pipe0.1651 m (6.5”)
DwoInner diameter of washout section0.2794 m (11”)
GElectrical conductance between electrodes of a probe
GmaxConductance in high-conductivity fluid
GminConductance in low-conductivity fluid
KGain factor of wave amplifier
LpLength of probe electrodes (length into annulus)5 mm
LsLength of each pipe section1.92 m
LtotTotal length of test section (ca)10 m
Lwo1Full bore length of washout section1.66m
Lwo2Length of increased diameter washout section1.805 m
QVolumetric flow rate
UAverage flow velocity in regular section of annulus0.5 m/s
VAnnular volume of the test section
VoMeasured voltage signal of a conductivity probe
Vo,maxMeasured voltage signal of a conductivity probe in high-conductivity fluid
Vo,minMeasured voltage signal of a conductivity probe in low-conductivity fluid
VoffsetOffset voltage signal of a conductivity probe
VrRelative voltage response
γ ˙ Fluid shear rate
γ ˙ s Fluid characteristic shear rate
ΔCenter-center distance of electrodes8 mm
εInner pipe eccentricity0 and −0.42
χDisplacement efficiency
δDiameter of electrodes3 mm
λResistance per unit length
θWellbore angle from vertical
σElectrical conductivity of fluid
τFluid shear stress
τsFluid surplus stress
τyFluid yield stress
ψElectrical resistivity of fluid

References

  1. Daccord, G.; Guillot, D.; Nilsson, F. Mud removal. In Well Cementing, 2nd ed.; Nelson, E.B., Guillot, D., Eds.; Schlumberger: Sugarland, TX, USA, 2006; pp. 143–189. [Google Scholar]
  2. Renteria, A.; Maleki, A.; Frigaard, I.A.; Lund, B.; Taghipour, A.; Ytrehus, J.D. Effects of irregularity on displacement flows in primary cementing of highly deviated wells. J. Pet. Sci. Eng. 2019, 172, 662–680. [Google Scholar] [CrossRef]
  3. Maleki, A.; Frigaard, I.A. Primary cementing of oil and gas wells in turbulent and mixed regimes. J. Eng. Math. 2017, 107, 201–230. [Google Scholar] [CrossRef]
  4. Tardy, P.M.J. A 3D model for annular displacements of wellbore completion fluids with casing movement. J. Pet. Sci. Eng. 2018, 162, 114–136. [Google Scholar] [CrossRef]
  5. Wilson, M.A.; Sabins, F.L. A laboratory investigation of cementing horizontal wells (SPE16928). SPE Drill. Eng. 1988, 3, 275–280. [Google Scholar] [CrossRef]
  6. Tehrani, A.; Ferguson, J.; Bittleston, S.H. Laminar Displacement in Annuli: A Combined Experimental and Theoretical Study. In Proceedings of the 67th Annual Technical Conference and Exhibition, Washington, DC, USA, 4–7 October 1992; Society of Petroleum Engineers: Richardson, TX, USA, 1992. SPE 24569. [Google Scholar]
  7. Lockyear, C.F.; Ryan, D.F.; Gunningham, M.M. Cement channeling: How to predict and prevent (SPE 19865). SPE Drill. Eng. 1990, 5, 201–208. [Google Scholar] [CrossRef]
  8. Deawwanich, T. Flow and Displacement of Viscoplastic Fluids in Eccentric Annuli. Ph.D. Thesis, The University of Adelaide, Adelaide, SA, Australia, March 2013. [Google Scholar]
  9. Jakobsen, J.; Sterri, N.; Saasen, A.; Aas, B.; Kjøsnes, I.; Vigen, A. Displacements in Eccentric Annuli During Primary Cementing in Deviated Wells. In Proceedings of the Production Operations Symposium, Oklahoma City, OK, USA, 7–9 April 1991; Society of Petroleum Engineers: Richardson, TX, USA, 1991. SPE 21686. [Google Scholar]
  10. Vefring, E.H.; Bjørkevoll, K.S.; Hansen, S.A.; Sterri, N.; Sævareid, O.; Aas, B. Optimization of Displacement Efficiency During Primary Cementing. In Proceedings of the Fifth Latin American and Caribbean Petroleum Engineering Conference and Exhibition, Rio de Janeiro, Brazil, 30 August–3 September 1997; Society of Petroleum Engineers: Richardson, TX, USA, 1997. SPE 39009. [Google Scholar]
  11. Kimura, K.; Takase, K.; Griffith, J.E.; Gibson, R.A.; Porter, D.S.; Becker, T.E. Custom-Blending Foamed Cement for Multiple Challenges. In Proceedings of the SPE/IADC Middle East Drilling Technology Conference, Abu Dhabi, UAE, 8–10 November 1999; Society of Petroleum Engineers: Richardson, TX, USA, 1999. SPE 57585. [Google Scholar]
  12. Clark, C.R.; Carter, L.G. Mud displacement with cement slurries (SPE 4090). SPE J. 1973, 25, 775–783. [Google Scholar]
  13. Zuiderwijk, J.J.M. Mud displacement in primary cementation. In Proceedings of the SPE-European Spring Meeting, Amsterdam, The Netherlands, 29–30 May 1974; Society of Petroleum Engineers: Richardson, TX, USA, 1974. SPE 4830. [Google Scholar]
  14. Crook, R.J.; Keller, S.R.; Wilson, M.A. Deviated Wellbore Cementing: Part 2 Solutions. J. Pet. Tech. 1987, 39, 961–966. [Google Scholar] [CrossRef]
  15. Smith, D.K. Chapter 7. In Cementing; SPE Monograph Series; Society of Petroleum Engineers: Richardson, TX, USA, 1990; Volume 4. [Google Scholar]
  16. Ytrehus, J.D.; Lund, B.; Taghipour, A.; Divyankar, S.; Saasen, S. Experimental investigation of wellbore fluid displacement in concentric and eccentric annulus. In Proceedings of the 36th International Conference on Ocean, Offshore and Arctic Engineering, Trondheim, Norway, 25–30 June 2017; American Society of Mechanical Engineers: New York, NY, USA, 2017. OMAE-62028. [Google Scholar]
  17. Lund, B.; Ytrehus, J.D.; Taghipour, A.; Divyankar, S.; Saasen, S. Fluid-fluid displacement for primary cementing in deviated washout sections. In Proceedings of the 37th International Conference on Ocean, Offshore and Arctic Engineering, Madrid, Spain, 17–22 June 2018; American Society of Mechanical Engineers: New York, NY, USA, 2018. OMAE-78707. [Google Scholar]
  18. Skadsem, H.J.; Kragset, S.; Lund, B.; Ytrehus, J.D.; Taghipour, A. Annular displacement in a highly inclined irregular wellbore: Experimental and three-dimensional numerical simulations. J. Pet. Sci. Eng. 2019, 172, 998–1013. [Google Scholar] [CrossRef]
  19. Saasen, A.; Lund, B.; Ytrehus, J.D. Theoretical Basis for Prediction of Drilling Fluid Removal in Annuli. In Proceedings of the 36th International Conference on Ocean, Offshore and Arctic Engineering, Trondheim, Norway, 25–30 June 2017; American Society of Mechanical Engineers: New York, NY, USA, 2017. OMAE-61030. [Google Scholar]
  20. Kennelly, A.E. The Linear Resistance between Parallel Conducting Cylinders in a Medium of Uniform Conductivity. Proc. Am. Phil. Soc. 1909, 48, 142–165. [Google Scholar]
  21. Saasen, A.; Ytrehus, J.D. Rheological Properties of Drilling Fluids–Use of Dimensionless Shear Rates in Herschel-Bulkley Models and Power-Law Models. Appl. Rheol. 2018, 28, 201854515. [Google Scholar]
  22. Lund, B.; Ytrehus, J.D.; Taghipour, A.; Saasen, A. Displacement Efficiency In Eccentric Annuli. In Proceedings of the 39th International Conference on Ocean, Offshore and Arctic Engineering, Virtual Conference, Fort Lauderdale, FL, USA, 3–7 August 2020; American Society of Mechanical Engineers: New York, NY, USA, 2020. OMAE-18707. [Google Scholar]
Figure 1. Test section with inlet to the right and outlet to the left. A torque motor for the rotating inner pipe is mounted at the outlet.
Figure 1. Test section with inlet to the right and outlet to the left. A torque motor for the rotating inner pipe is mounted at the outlet.
Energies 13 05201 g001
Figure 2. Process and instrumentation diagram for the flow rig. Note that probe C0 and valve V3 were positioned about 9.5 m downstream of the test section exit point.
Figure 2. Process and instrumentation diagram for the flow rig. Note that probe C0 and valve V3 were positioned about 9.5 m downstream of the test section exit point.
Energies 13 05201 g002
Figure 3. Details of the irregular section.
Figure 3. Details of the irregular section.
Energies 13 05201 g003
Figure 4. Bearing joint connecting two sections of inner pipe.
Figure 4. Bearing joint connecting two sections of inner pipe.
Energies 13 05201 g004
Figure 5. Piping and valves at inlet to test section. The valves V1 and V2 of Figure 2 were controlled from the PC using pneumatic actuators. A drain valve (blue handle) with corresponding hose (not shown in Figure 2) was used to drain the test section after an experiment.
Figure 5. Piping and valves at inlet to test section. The valves V1 and V2 of Figure 2 were controlled from the PC using pneumatic actuators. A drain valve (blue handle) with corresponding hose (not shown in Figure 2) was used to drain the test section after an experiment.
Energies 13 05201 g005
Figure 6. Schematic cross section of test section showing position of conductivity probes and inner pipe in eccentric position. Pressure ports (not shown) were located between probes 6 and 7.
Figure 6. Schematic cross section of test section showing position of conductivity probes and inner pipe in eccentric position. Pressure ports (not shown) were located between probes 6 and 7.
Energies 13 05201 g006
Figure 7. Principle of operation for conductivity measurement. Lower density displaced fluid (yellow) stratified above higher density displacing fluid (blue).
Figure 7. Principle of operation for conductivity measurement. Lower density displaced fluid (yellow) stratified above higher density displacing fluid (blue).
Energies 13 05201 g007
Figure 8. Measured conductivity (points) versus volume fraction of fluid B0. Stippled line shows the linear regression of the measurement points.
Figure 8. Measured conductivity (points) versus volume fraction of fluid B0. Stippled line shows the linear regression of the measurement points.
Energies 13 05201 g008
Figure 9. Fluid conductivity measured during 10 consecutive displacement experiments using the conductivity probe C0.
Figure 9. Fluid conductivity measured during 10 consecutive displacement experiments using the conductivity probe C0.
Energies 13 05201 g009
Figure 10. Normalized fluid flow rate versus time from start of displacement for the experiments with the automatic pump control.
Figure 10. Normalized fluid flow rate versus time from start of displacement for the experiments with the automatic pump control.
Energies 13 05201 g010
Figure 11. Local fraction of displacing fluid near the inlet (station A) and near the top (positions 1 and 8) for an experiment with inclined (60 degrees) concentric annulus (experiment N), and a repetition experiment (R), both without rotation of inner pipe.
Figure 11. Local fraction of displacing fluid near the inlet (station A) and near the top (positions 1 and 8) for an experiment with inclined (60 degrees) concentric annulus (experiment N), and a repetition experiment (R), both without rotation of inner pipe.
Energies 13 05201 g011
Figure 12. Displacement versus normalized pump volume for equivalent experiments in Experiments A and O. All measurements were conducted at Station A being 2 m from inlet.
Figure 12. Displacement versus normalized pump volume for equivalent experiments in Experiments A and O. All measurements were conducted at Station A being 2 m from inlet.
Energies 13 05201 g012
Figure 13. Definition of azimuthal angle θ and angle-dependent gap H(θ) for eccentric annulus.
Figure 13. Definition of azimuthal angle θ and angle-dependent gap H(θ) for eccentric annulus.
Energies 13 05201 g013
Figure 14. Experiment in concentric annulus with 60 degrees well inclination. Time traces of pressure gradients (DP1: black, DP2: red), flow rate (FM: green), conductivity of mixed fluid downstream of test section (C0: yellow), and temperature of mixed fluid downstream of test section (T0: blue).
Figure 14. Experiment in concentric annulus with 60 degrees well inclination. Time traces of pressure gradients (DP1: black, DP2: red), flow rate (FM: green), conductivity of mixed fluid downstream of test section (C0: yellow), and temperature of mixed fluid downstream of test section (T0: blue).
Energies 13 05201 g014
Figure 15. DP1 pressure gradients for the experiments shown in Figure 14 (black) and a repetition test (red).
Figure 15. DP1 pressure gradients for the experiments shown in Figure 14 (black) and a repetition test (red).
Energies 13 05201 g015
Figure 16. Measured pressure gradients DP1 for horizontal and inclined wellbore from experiments without inner pipe rotation, and corresponding frictional pressure gradients DP1_f, assuming a fluid density difference of 100 kg/m3. Crosses show corresponding calculations using the values from Table 5.
Figure 16. Measured pressure gradients DP1 for horizontal and inclined wellbore from experiments without inner pipe rotation, and corresponding frictional pressure gradients DP1_f, assuming a fluid density difference of 100 kg/m3. Crosses show corresponding calculations using the values from Table 5.
Energies 13 05201 g016
Figure 17. Effect of inner pipe rotation on displacement in inclined wellbore (60 degrees from vertical) with eccentric annulus (e = −0.42). Curves showing fraction of displacing fluid at station C (near inlet to washout section) without inner pipe rotation (experiment E1) and with inner pipe rotation (experiment F1). Azimuthal positions of probes (1, 5 and 8) are as shown in Figure 6.
Figure 17. Effect of inner pipe rotation on displacement in inclined wellbore (60 degrees from vertical) with eccentric annulus (e = −0.42). Curves showing fraction of displacing fluid at station C (near inlet to washout section) without inner pipe rotation (experiment E1) and with inner pipe rotation (experiment F1). Azimuthal positions of probes (1, 5 and 8) are as shown in Figure 6.
Energies 13 05201 g017
Figure 18. Displacement in washout section for horizontal flow in concentric annulus (a) and for eccentric annulus (b). Arrows indicate position of interface between displaced fluid and displacing fluid, where displacing fluid is below arrows (a) or below and to the right of arrows (b).
Figure 18. Displacement in washout section for horizontal flow in concentric annulus (a) and for eccentric annulus (b). Arrows indicate position of interface between displaced fluid and displacing fluid, where displacing fluid is below arrows (a) or below and to the right of arrows (b).
Energies 13 05201 g018
Figure 19. Measurement of displacing fluid fraction immediately upstream and downstream of the washout section for an experiment in horizontal wellbore with concentric inner pipe. Probes C1/D1 and C8/D8 measure fluid fraction near top upstream/downstream of washout section, while probes C5/D5 measure corresponding fluid fractions near the bottom.
Figure 19. Measurement of displacing fluid fraction immediately upstream and downstream of the washout section for an experiment in horizontal wellbore with concentric inner pipe. Probes C1/D1 and C8/D8 measure fluid fraction near top upstream/downstream of washout section, while probes C5/D5 measure corresponding fluid fractions near the bottom.
Energies 13 05201 g019
Table 1. Fluid composition by mass percentage for the fluids used.
Table 1. Fluid composition by mass percentage for the fluids used.
Fluid System 1Fluid System 2
ComponentDisplacedDisplacingDisplacedDisplacing
Tap water98.68598 73
Sucrose000 24.5
Laponite1.401.5 0
Xanthan Gum00.670.1 0
CMC0000.78
NaCl00.801.2
Soda Ash00.3700.42
Biocide00.090.080.1
Barite012.500
Table 2. Flow rig key parameters (see Figure 2).
Table 2. Flow rig key parameters (see Figure 2).
ParameterSymbolValue
Total length of test section (ca)Ltot10 m
Inner diameter of annulus (outer diameter of inner pipe)Di0.127 m (5”)
Inner diameter of gauge outer pipeDo0.1651 m (6.5”)
Inner diameter of washout sectionDwo0.2794 m (11”)
Length of each pipe sectionLs1.920 m
Full bore length of washout sectionLwo11.660 m
Length of increased diameter washout sectionLwo21.805 m
Eccentricityε0 and −0.42
Table 3. Differential pressure transmitters.
Table 3. Differential pressure transmitters.
NameTypePosition
DP1 (DP10)Fuji FCX, Type FKCW35V2AXAYYAE, Range-130 to 130 kPaBetween 2 m and 4 m from inlet
DP2 (DP4)Fuji FCX-AII, Type FKCW35V4AXCYYAE, Range-130 to 130 kPaBetween 4 m and 6 m from inlet
DP3 (DP1)Fuji FCX, Type FKCW43V2-AXCYY-AF, Range-32 kPa to +32 kPa.Between 6 m and 8 m from inlet (across irregular section)
Table 4. Video camera positions.
Table 4. Video camera positions.
CameraPosition
1Section 2 (between stations A and B), vertical view from bottom
2Section 2 (between stations A and B), side view
3Section 3 (between stations B and C), side view
4Section 4 (washout section, between stations C and D), vertical view from bottom
5Section 4 (washout section, between stations C and D), side view
6Section 5 (after station D), side view
Table 5. Transport properties for the barite free two-fluid systems.
Table 5. Transport properties for the barite free two-fluid systems.
PropertiesDisplacedDisplacing
Fluid density ρ (kg/m3)10001100
Yield stress τy (Pa)4.91.8
Consistency index K (Pa.sn)0.470.87
Fluid flow index n (-)0.550.57
Surplus stress τs (Pa) @ γ ˙ s = 158 s-17.215.3
Reynolds number @ U = 0.5 m/s (-)199164
Pressure gradient @ U = 0.5 m/s (Pa/m) (eccentric/concentric)1293/15801729/2108

Share and Cite

MDPI and ACS Style

Lund, B.; Taghipour, A.; Ytrehus, J.D.; Saasen, A. Experimental Methods for Investigation of Drilling Fluid Displacement in Irregular Annuli. Energies 2020, 13, 5201. https://doi.org/10.3390/en13195201

AMA Style

Lund B, Taghipour A, Ytrehus JD, Saasen A. Experimental Methods for Investigation of Drilling Fluid Displacement in Irregular Annuli. Energies. 2020; 13(19):5201. https://doi.org/10.3390/en13195201

Chicago/Turabian Style

Lund, Bjørnar, Ali Taghipour, Jan David Ytrehus, and Arild Saasen. 2020. "Experimental Methods for Investigation of Drilling Fluid Displacement in Irregular Annuli" Energies 13, no. 19: 5201. https://doi.org/10.3390/en13195201

APA Style

Lund, B., Taghipour, A., Ytrehus, J. D., & Saasen, A. (2020). Experimental Methods for Investigation of Drilling Fluid Displacement in Irregular Annuli. Energies, 13(19), 5201. https://doi.org/10.3390/en13195201

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop