1. Introduction
Even as a general global goal is to consume energy from renewable energy sources, oil and gas are needed in the transition [
1]. Oil production is expected to increase during the next three decades globally, which entails the oil and gas industry to significantly impact the world’s energy consumption in the coming decades [
2]. The general trend towards more sustainable energy production also affects the offshore oil and gas industry as discharge legislation becomes stricter [
3,
4]. With stricter policies, produced water re-injection (PWRI) has gained growing attention to extend the reservoirs’ economic life and decrease produced water (PW) discharge to minimize environmental impacts [
5]. However, the re-injected PW and the injected seawater must have a continually high quality to prevent formation damage and unpredicted injectivity reduction [
6,
7,
8,
9,
10]. Currently, only ∼14.5% of the PW in the Danish sector of the North Sea is re-injected, which is a >50% reduction since 2009 due to reservoir challenges according to the Danish Environmental Protection Agency [
3,
11]. To increase the PWRI percentage, efficient management is essential to maintain a high PW and seawater quality, which involves proper treatment and accurate monitoring. Monitoring the PW also support protecting the receiving environment when discharging. Accurate information on the amount of oil and particles, sizes, and classification of particles in the injection water (IW) can essentially be used for decision support, reporting, or even advanced control to achieve better operation in the treatment process [
6,
12]. Therefore, the importance and awareness of accurate online measurements of oil and suspended particles have gained increasing attention within the oil and gas industry [
5,
6].
The water quality in offshore injection water treatment (IWT) processes is usually assessed in terms of particles’ plugging tendency, also addressed as suspended solids or total suspended solids (TSS) [
8,
13]. However, the IW quality should be economically viable to achieve and counterbalanced against other solutions, e.g., well-stimulation [
12,
13]. Currently, offshore IWT processes rely on off-line measurements of TSS concentration following the ISO 11923; thus, in case of a decrease in IW quality, the reaction time is long due to the use of onshore laboratory measurements [
14].
Even if at-line execution significantly reduces the reaction time of off-line TSS concentration measurements, it is difficult for operators to investigate where and what is the root cause of the decreased quality when there is no available information related to TSS. Especially as significant water quality decline often occur downstream [
8]. Online TSS and oil-in-water (OiW) concentration measurements have not become standardized in the oil and gas industry, despite the long history of measuring particle sizes online [
14]. The measurement technique is often not the main reason for the error source in the observed results [
15]. Different studies emphasize that the main source of errors often are due to improper installment and calibration [
15,
16,
17].
According to Latif [
18] only a few research studies of continuous measurement techniques for online concentration measurements have been published. Latif [
18] and Fjords Processing [
19] have done a comparison study on different techniques for measuring OiW concentration. However, their focus is on how different variations of parameters such as mixtures, pressures, chemicals, and temperatures can affect the measurements and not on the uncertainty related to the petroleum engineers’ ability to determine an acceptable calibration. Although PWRI is a solution for extending the economic life of oil production, it is beneficial to measure particle sizes to increase IWT processes’ performance. Several different methods for measuring particle sizes exist, each based on several design options from different manufacturers. Besides microscopy, which is based on direct observation, all other techniques are challenged by their property assumptions that only the equivalent diameter of a sphere is measured [
12]. The existing morphologies are not taken into account, which complicates the ability to classify particles [
12]. Another advantage of microscopy is the manual discrimination of particles captured to evaluate the results. This paper aims to examine two different online microscopy analyzers: Jorin ViPA and Canty InFlow.
Both microscopes are based on the same technique of utilizing a high-resolution video camera to capture images of the particles passing the view cell. Both monitors are based on bright field illumination techniques, presenting a dark image of particles passing with a bright background. The pixel scale of Jorin ViPA is fixed by the manufacturer, and the Canty InFlow is adjusted to a specific pixel scale. A fluorescence-based monitor (Turner TD-4100XDC) is used as a benchmark to evaluate the two online microscopes’ performance to measure OiW concentrations. The fluorescence-based monitor is selected based on a previous work by Hansen et al. [
6]. This paper addresses the problems occurring when petroleum engineers must subjectively decide what particles are considered in focus. Furthermore, the calibration procedure will be validated on known solid particle sizes and their ability to measure different OiW concentrations accurately and in real-time. Especially to determine the accurate volumetric concentration in each image is not trivial as the depth of field (DoF) would affect the estimate of particle concentration.
2. Materials and Methods
Two setups were constructed to execute the experiments presented in this paper. The two systems are shown in
Figure 1 and
Figure 2. The setup shown in
Figure 1 is used to calibrate both online microscopes using different known polystyrene particle sizes produced by BS-Partikel [
20].
The three different polystyrene particle sizes follow a statistically Gaussian distribution, and the measured means and standard deviations presented in parentheses, are given by the manufacturer:
Nominal diameter of µ: µ
Nominal diameter of µ: µ
Nominal diameter of µ: µ
Figure 2 illustrates a skid-mounted testing platform for validating the two online microscopes’ performance. The platform is equipped with a centrifugal pump (CP), pressure transmitters (P
x), flowmeters (Q
x), and pneumatic control valves (V
x). The system is configurable to direct the liquid through the sidestream and the mainstream by manipulating the control valves.
By taking advantage of on-line sampling, both microscopes: Jorin Visual Process Analysis (ViPA) B HiFlo (C2) and Canty InFlow VD4912-960 (C3) are installed on a sidestream, making them applicable in most installations regardless of the flow velocity and the pipeline’s dimensions. Although, on-line measurement complicates the sampling procedure as a maldistribution between the run and the branch can occur.
Both microscopes mainly consist of a camera, light source, flow cell, optical lens, and a computer. For C
2 a progressive CCD camera (Sony ICX445) with
captures images at maximum frame rate at full resolution of 31 frames per second (fps) with sensor pixel size of
µ
. A Comar 10 OA 25 lens is used to output a field of view of
µ
with a magnification of 10, a numerical aperture of 0.25, and a focal length of 16.9
. The manufacturer of C
3 retains most of the equipment specifications confidential, however some general specifications of both microscopes are listed in
Table 1.
A fluorescence-based monitor (C
1) that is sensitive to the aromatic content is installed to measure the OiW concentration. C
1 is installed as a benchmark to facilitate the investigation of measuring OiW concentrations by C
2 and C
3. The accuracy of C
1 has been extensively studied by Hansen et al. [
6], and will not be evaluated in this paper.
The difficulty of online microscopy analysis is related to the narrow focus area, where only a fraction of the entire flow is directed into the sidestream and through the view cell. A further complication occurs as only a narrow DoF of the passing flow is measured by the microscope; thus, not all TSS that passes the view cell will be observed and are therefore excluded from being measured as shown in
Figure 3.
To achieve statistical results, it is necessary to sample a sufficient amount of particles that is constricted by the small DoF. The quantity of captured particles strongly relates to the accuracy of the predicted particle size distribution (PSD). The accuracy of the PSD can be determined within a region of relative error,
, from a defined confidence level with its represented z-score value,
u. The number of sufficient particles,
, that are required is based on
and sample standard deviation,
s, has been proposed by Masuda and Iinoya [
21]. The approach is based upon the assumption that the PSD follows a log-normal distribution.
where
is the non-zero exponential constant that defines the particle distribution, for log-normal distribution
[
22].
c can be calculated as
where
is a basis number, for the count basis
. The number of particles required within a certain range of error changes depending on how the data is scattered. Thus, the number of particles increases as the deviation of the data increases. Simultaneously, the number of particles increases proportionately to the error of magnitude required. On the contrary, the relative error can be estimated based on the number of particles observed.
Experimental Design
This subsection describes the design and objective of all executed experiments on both setups from
Figure 1 and
Figure 2. The two microscopes were calibrated according to a known particle size from BS-Partikel, followed by a performance validation on known sizes, larger and smaller than the known particle size used for calibration. The two microscopes were installed in series together with the fluorescence-based monitor to observe their performance of measuring OiW concentrations. The following experiments are executed on the setup presented in
Figure 1:
Experiment I: calibration of both microscopes, by addition of known particle sizes with a mean and standard deviation: µ.
Experiment II: validation test of the calibration procedure executed in Experiment I, by addition of a known particle size larger and smaller than : µ and µ.
The pump speed was fixed at the same value in both Experiment I and II, which generated a flow rate of ∼2 . Most of the presented data in Experiment I are based on data from C2, although the exact same procedure was executed on C3. Both microscopes have three calibration parameters:
Threshold value (THV): an 8-bit integer, resulting in a greyscale image with pixel values in the range of 0–255, from black to white for both microscopes.
Edge strength/focus rejection value (ESV/FRV): an edge detection method, running a convolution kernel for estimating the gradient at each pixel on the image. The most common method is Sobel, which has been used on C3. Other methods are also available when using C3. C2’s edge detection method is confidential. The ESVs of C2 are in a range of 0–10, where ESV only includes particles with a large gradient. C3’s FRV ranges from 0–1000, although FRV did not include any particles in the calibration analysis.
Depth of field (DoF): defines the depth of the captured images. DoF normally defines the distance between the closest and farthest particle in the image that appears acceptably sharp. DoF is only necessary to adjust if the user is interested in the sample’s concentration as the DoF influence to field of view volume. To determine the DoF for these microscopes, the DoF is adjusted to present the known concentration based on the accepted measured particles fulfilling THV and ESV/FRV. The DoF is then estimated using
where
C is the known concentration,
k is the number of images,
n is the number of particles in each
jth image,
is the equivalent area diameter for each
ith particle in each
j-th image,
is the field of view, and
is the average DoF. By adjusting
in (
4) of the captured images within a defined time interval, it is possible to adjust the measured concentration from the microscopes to match the known concentration.
Only the first two calibration parameters were analyzed in Experiment I and validated in Experiment II by the addition of the particle sizes:
and
. The DoF value was tuned in Experiment III. The following experiment is executed on the setup presented in
Figure 2:
Experiment III: a performance evaluation of the microscopes, to measure OiW concentrations in steady-state and real-time measurements, benchmarked according to the known OiW concentration and the measurements obtained from C
1. Six nominal OiW concentrations were investigated
. The supply tank was filled with 163.59
of 1 µ
filtered tap water. The solution of oil added are listed in
Table 2 for each OiW concentration.
CP was kept at 100% pump speed to ensure that the droplet size distribution does not change over time due to the high shear occurring in CP. V2 was kept at a fixed opening degree 37.5%, and the flow rate through the sidestream was kept constant by manipulating V1.
3. Results
The results are divided into three sections based on Experiment I, II, and III.
Experiment I:
A 2 recirculation experiment with particle size dispersed in 1 µ filtered tap water, was executed with both microscopes and used as calibration data.
Firstly, the THVs were adjusted so that the known particle size can statistically be represented by the microscopes. A high ESV was selected only to determine the THVs on captured particles in focus. Determining the THV before the ESV is essential as the edge detection method depends on the THV.
Furthermore, for comparison of both microscopes’ performance, particle sizes were based on equivalent area diameter during post-processing of the data. Note that choosing another method of determining the particle diameter will most likely shift the distribution of measured sizes. The measurement in
Figure 4 presents the mean results of the first 100 images with particles captured at maximum ESV. The THV was then toggled in a range of 47–81, where at THV
and THV
, particles no longer fulfill the ESV’s criteria and are no longer counted. The Gaussian distribution in
Figure 4 represent the information analyzed by BS-Partikel:
µ
and
µ
.
Based on the results in
Figure 4, a THV
for C
2 was selected. For C
3, a THV
was selected. The analysis results with a THV
and maximum ESV
are shown in
Figure 5. A probability density function (PDF) is included on the graphs in
Figure 5, to visualize the performance of estimating the particle sizes.
Selecting the ESV is more subjective as it is a perception of what appears to be acceptably sharp.
Figure 6 presents eleven different particles captured by C
2 with a THV
, each with different acceptable ESV between 0–10. Associated information of each particle, presented in
Figure 6, is shown in
Table 3.
Furthermore,
Figure 7 and
Table 4 show an artifact of selecting a weak ESV. Even though an ESV of one outputs a measured size within
from the
of
, selecting an ESV
, the gradient of the particle is no longer accepted and is considered out-of-focus. Although a smaller, but stronger gradient within the particle, fulfills the requirement and is included as a particle.
Selecting a too high ESV will reduce the counted number of accepted particles in the analysis, or even worse, not represent other types of particles in the process that has a weaker gradient by nature. A further validation of selecting a proper ESV or FRV was executed by analyzing the corresponding results of the first three particles in
Table 3 that are within
from
of
(ESV: 5, 6, and 7). The results are shown in
Figure 8. The same evaluation procedure was executed on C
3.
Increasing the ESV from 5 to 7 counted 24% fewer particles in Experiment I. Another phenomenon occurs at ESV
, as the size histogram is skewed left due to the acceptable edge is found closer to the center of the particles as the peripheries are less in focus. Selecting between an ESV
or ESV
is a tradeoff between obtaining less countable particles and reducing the phenomena of underestimating the size of particles that are less in focus. ESV
was chosen for C
2, and FRV
was chosen for C
3. The calibration results of C
3 are shown in
Figure 9. The number of particles captured by C
3 was ∼5 times larger than the amount captured by C
2. This is primarily a result of higher image resolution and larger pixel size than C
2 as shown in
Table 1.
Experiment II:
A 2
validation test was executed by adding two additional particle sizes:
and
, along with
. The raw data from C
2 and C
3, without any classification, are shown in
Figure 10.
For each known particle size, 1
aqueous surfactant solution containing the polystyrene particles was added. Both microscopes observed a high number of
in the validation test, a result of containing more particles per volume than the others. By truncating the sample for each known particle size with a fixed range of
µ
from the known
of
,
, and
, the statistical information can be obtained as shown in
Table 5.
Experiment III:
A 9.5 experiment was carried out with the addition of oil after at least 30 between each OiW concentration of observation.
Figure 11 shows flow rates of the mainstream (Q
1) and sidestream (Q
2), and the opening degree of V
1. The sidestream was sufficiently kept constant at 3
during the entire experiment by manipulation of V
1.
Figure 12 shows the time series of measured OiW concentration from C
1. The vertical line marks a truncated time series of ∼3.5
that was left out. The DoF value in C
2 and C
3 was at that period adjusted to match a OiW concentration of 55
during the first 30
after the injection of oil into the process by use of (
4) and to classify oil droplets. The DoF value in C
2 and C
3 was estimated to be
µ
and
µ
, respectively.
During the last known OiW concentration of
, the mixed concentration can no longer be maintained in the setup according to the measurement of C
1.
Figure 13 presents an error bar of each OiW concentration measured in a duration of 30
for each online monitor. Each error bar represents the mean, minimum, and maximum OiW concentration measured, where C
1 was selected to have a sampling frequency of 0.1
, and C
2 and C
3 outputting an averaged concentration measurement every minute.
Table 6 shows the associated results related to
Figure 13.
Figure 14 and
Figure 15 present the PSD histogram of the measured OiW concentrations at 55
and 400
for each microscope. Furthermore, the
Figure 14 and
Figure 15 presents the mean of the log-normal distribution,
, the standard deviation of the log-normal distribution,
, and the number of counted particles,
.
Both OiW concentrations follow a log-normal size distribution measured by both microscopes, with only small changes in the distribution with respect to the concentration. Following the proposed calculation by Masuda and Iinoya [
21] in (
1)–(
3), for obtaining a sufficient statistical representation within 95% confidence level,
of the PSD can be calculated as shown in
Table 7.
A trailing moving average window of 1
was selected, for measuring the OiW concentration real-time by C
2 and C
3.
Figure 16 and
Figure 17 show the results of real-time measurements at the known OiW concentration of 55
from C
2 and C
3, respectively, together with real-time measurements from C
1. The calculated OiW concentration, based on post-processed equivalent volume measurements, are presented for each captured image.
4. Discussion
The process of calibrating the two online microscopes is highly based on the perception of what is considered in focus. To ensure consistent results when calibrating the three main calibration parameters: THV, ESV/FRV, and DoF, the calibrating procedure must be well documented to reduce the difference between each calibration execution.
Figure 4 and
Figure 8 presented one solution for selecting the THV and ESV/FRV that outputs an estimation that matches the Gaussian distribution of the known particle sizes based on the information from the manufacturer of the produced particles. The skewness observed in
Figure 8 when calibrating the ESV could result from selecting a slightly too low THV; thus, a lower ESV may have yielded better results. Although, based on the results shown in
Figure 10 the calibration procedure presented in Experiment I revealed only a small deviation from the known PSDs. To further improve online microscopy analyzers, auto-calibration or human-computer interaction (HCI) of the calibration procedure can lower the subjective perception of what is acceptably sharp.
The steady-state values of measuring OiW concentration in
Figure 13 showed promising results in relation to the known OiW concentrations. The measurement of C
1 confirmed that C
2 and C
3 are able to measure the mean OiW concentrations by estimating the DoF using (
4). The DoF for C
3 was close to the theoretical estimation of DoF from the manufacturer seen in
Table 1. The DoF for C
2 was
% lower than the theoretical estimation of DoF given by the manufacturer.
Moreover, C1 also confirmed that the setup was able to keep a well-mixed OiW concentration, apart from high OiW concentrations where a declining tendency of measured concentration was observed. The observed declining tendency may occur due to accumulation within the setup (dead volumes), natural separation in the supply tank, or pipelines and equipment becoming oil-wet.
Regardless of these sources of error, both microscopes were able to measure the mean OiW concentration of 30
in close relation to the known concentration. A high deviation was observed when outputting OiW concentration every minute in steady-state, especially considering the steady concentration measured by C
1 as shown in
Figure 12.
Figure 20 shows the results of extending the averaging time duration of both microscopes in steady-state. The results in
Figure 20, show that increasing the averaging time to 3
reduces the fluctuation of measured OiW concentration greatly. Meanwhile, the negative impact is that the time resolution of dynamics will be reduced accordingly.
Another solution is to measure the OiW concentration in real-time by incorporating a trailing moving average window. The size of the moving average window is highly based on the plant’s dynamics. Based on the results in
Figure 16,
Figure 17,
Figure 18 and
Figure 19, both microscopes are able to measure the the OiW concentration and thus they have the potential to track the transient behavior when it occurs in an IWT process. The fluctuating real-time results obtained with both microscopes, using a trailing moving average window of 1
, can be questioned to be sufficient to provide qualitative feedback to an operator. However, by extending the moving average window to 3
, the real-time measurements of OiW concentration from both microscopes are naturally more consistent, as shown in
Figure 21. Extending the moving average window is only valuable if the dynamic of the plant is not faster than the moving average window. The higher resolution of C
3 caused the microscope to capture twice as many particles compared to C
2, thus when using the moving average window, C
3 yields a more stable measurement.
One important procedure that has not been addressed in this paper is the classification within both microscopy analyzers. Oil droplets in Experiment III have been classified and withdrawn from the rest of the observed particles. The classification process has not been addressed as the procedure are very different in each microscope. C
2 entrust the operator in selecting the classification range of each measured variable. C
3 uses machine learning by letting an operator train the classification process by manually defining the captured particles. The different results of PSDs shown in
Figure 14 and
Figure 15, could be a result of classifying differently. The results obtained by C
3 are not representable around 13 µ
if the PSD truly follows a log-normal distribution. Another consequence to the bimodal PSD obtained by C
3 could be the installment right after the narrow view cell of C
2, which may change the PSD of oil droplets by breakups or coalescence. Another statical consideration was addressed when outputting a PSD, as the PSD should be trustworthy if an operator should draw any conclusion based on the presented PSD. Masuda and Iinoya [
21] presented an estimation for predicting the required number of particles to be counted for a certain accuracy. In this paper,
of the PSDs shown in
Figure 14 and
Figure 15 were calculated based on the number of particles obtained within 30
at two different OiW concentrations (
and
) for both microscopes. Based on the results in
Table 7, a
was observed within a 95% confidence interval of all the presentable PSDs. If the
is not enough more particles must be counted, i.e., at
roughly 100,000–130,000 for C
3 and 90,000–100,000 for C
2 must be counted based on the measured OiW concentrations, respectively. Seen from another point of view, if a PSD update every 30
is too slow, the allowed
in the presented PSD must be defined to determine the minimum number of particles required in order to achieve a defined level of accuracy.