1. Introduction
Pressure–volume–temperature (PVT) analysis has been living at the heart of oil and gas reservoir engineering since the early days of the 20th century when it was realized that reservoir fluid behavior against pressure and temperature would have an immense effect to the pressure state of the reservoir, hence its ability to keep producing fluids, as well as to the volume and quality of the produced fluids at surface. Thermodynamic effects including oil and gas shrinking and swelling with pressure depletion, gas dissolution in the oil phase and liquid dropout (condensation) are greatly responsible for the productivity and recovery factor of any specific field [
1]. Kinetic properties are directly related to the pressure drawdown required to establish the desired flow rate within the porous medium. When it comes to tertiary production, thermal properties such as heat capacity and the Joule–Thompson coefficient will also have a huge impact [
2,
3].
Strictly speaking, the fundamental phenomena applying in a reservoir are gravity, viscous flow and capillarity, with reservoir fluid properties being involved in all three of them. Fluid density, regarding oil, gas or water, is responsible for buoyancy effects and therefore for the vertical flow of fluids within the pore space [
4]. Isothermal compressibility is directly related to the pressure change once some volume has been withdrawn (or injected) to the reservoir. The more compressible the fluid, the more primary production is facilitated, thanks to the remaining fluid expansion [
5]. Dissolution effects are also involved in pressure maintenance as the more is the gas that is liberated in the reservoir, the more the total system compressibility and the more conservative the pressure reduction. Viscosity contributes to viscous flow as the expected flowrate for a fixed pressure drawdown is inversely linearly related to fluid viscosity. Finally, the rock wettability is determined by the rock–fluid interaction properties such as the relative permeability and the capillary pressure curves [
6].
When a reservoir is undersaturated, the contained fluid remains in a single phase (oil or gas); hence, there is no dissolution change along with pressure. However, when the pressure is locally or globally allowed to decline below the saturation one, gas is released from the reservoir oil, or condensate forms when a gas reservoir is considered [
7]. As a result, the remaining oil or gas composition and properties do vary over time as the pressure changes. Oil becomes heavier, more viscous and contains less dissolved gas, whereas gas becomes thinner. On the other hand, the liberated gas composition and properties keep changing along with pressure as different compounds get released at each pressure step due to the variance in their volatility [
8].
To model those effects, once an oil bottom hole sample (BHS) or a recombined surface oil sample (RSS) becomes available, it undergoes some standardized experimental PVT study, known as differential liberation (DL) or differential vaporization (DV) [
9]. The idea is quite simple: some amount of reservoir oil at saturation conditions is loaded in a PVT cell and it is gradually driven to atmospheric conditions by decreasing the pressure in steps under constant temperature equal to the reservoir one. At each step, the initially single phase, saturated oil turns into a diphasic one as gas forms a new phase at the top of the PVT cell. Subsequently, the gas is removed in an isobaric fashion so that a new saturated oil remains trapped in the cell. The process is repeated until pressure arrives at the atmospheric one. The gas collected at each step is driven at standard conditions and its volume, mass and composition are measured. The volume of the remaining saturated oil is also recorded. Finally, the residual oil is cooled down at standard temperature [
10]. The volumetric data collected at each step are utilized to calculate useful oil properties such as the formation volume factor (oil FVF,
), solution gas (
) and density (
), as well as released gas properties such as the formation factor (gas FVF,
), compressibility factor (
factor) and specific gravity (
) [
11]. The abovementioned volumetric properties are typically utilized to generate PVT tables [
12,
13] or to tune an equation of state (EoS) model so that a physically sound model can later be introduced into the compositional reservoir simulator [
14,
15,
16,
17].
Errors in the PVT results may arise mainly from three sources [
18]. Firstly, errors arise when the accuracy of raw data measured at the lab is compromised [
19]. Volume measurements, mass estimates and liberated gas composition analyses are often prone to minor or major lab errors which eventually lead to distorted PVT values found in the fluid report [
20]. Indeed, since PVT experiments are typically run using small volumes of reservoir fluids, a slight change in the volume at a specific pressure step—just a fraction of a cubic centimeter—can be comparable to the device measuring accuracy. Moreover, errors can arise during mass measurements of liberated gas at each pressure step, as well as during the measurement of the residual oil mass. These occur as both fluids are transferred by pressure difference through pipelines from the PVT cell to the measuring devices, where some mass may become trapped, thus introducing inaccuracies to the PVT values reported. Another source of error stems from differences between laboratory conditions and standard surface conditions (14.7 psi, 60 °F), necessitating volume conversions, which can further contribute to errors [
21,
22]. Additionally, the inadequate consideration of equipment calibration by lab staff is a potential source of error. This issue occurs when parts of the equipment, such as the gas chromatography (GC) column, the pump or the thermometer/manometer, are changed without subsequent re-calibration. Unfortunately, this practice is often carried out in an attempt to save time, overlooking the critical importance of ensuring that all equipment is properly calibrated after any modifications. This oversight can lead to significant inaccuracies and compromises in the reliability of the lab results. Finally, another cause of severe errors is the failure to achieve true thermodynamic equilibrium at each pressure step. This becomes particularly pronounced when dealing with fluids of high volatility, where the pressure step conditions closely approach critical points. In such cases, achieving thermodynamic equilibrium between liberated gas and saturated oil compositions requires extended time for the step flash to converge. Although this is a major concern when handling gas condensates, it still can introduce severe errors when analysing volatile oil samples.
Secondly, errors can occur during the handling of raw data when they are converted to the delivered PVT values. Although modern PVT data handling software takes care of various checks such as the material balance, various errors have been observed in real-life applications. For example, not thoroughly checking the curve fitting of the resulting PVT values may distort the original values and introduce severe inconsistency to the smoothed data [
23,
24,
25], as lab measurements taken at distinct pressure steps are typically interpolated using appropriate functions and delivered to the client in that form. Finally, an issue arises from the digitization of old paper-printed PVT studies which are often used either to tune an EoS model or to compare against measurements taken on newly collected samples. Poor optical character recognition (OCR) software performance [
26,
27,
28] or even mistakes while reporting the computed values (which had been performed manually until recently) may also introduce various value errors [
18].
The information which probably suffers the most is the GC-measured composition of the released gas at each pressure step which PVT labs also provide to their clients [
29]. Specifically, although the liberated fluid at each DL step is a gas, during the last few steps, when the pressure is low, heavier components’ concentration becomes important. It is then when the composition is often mismeasured because the appropriate handling of the heavy end is not part of the analysis workflow. A similar problem holds for the analysis of the residual oil composition where heavy components dominate the composition.
There has been a very long discussion on the usefulness of GC-measured composition data. For example, when it comes to compositional reservoir simulation, liberated gas composition demonstrates the preference of each component to get released from the oil with pressure; therefore, it might be used to further tune the EoS model (although, to the best of our knowledge, such an option is not directly available in most of the commercial PVT software). Alternatively, it can be used to evaluate the accuracy of the obtained lab data by incorporating it to a mass balance check per component [
30]. Clearly, if gas composition is to be utilized, it needs to be quality checked prior to any application. Nowadays, PVT labs either skip some measurement so that the quality control (QC) of the measurements by means of mass balance is not possible, or they may “smooth” the data to honor the mass balance [
23]. Legacy PVT reports typically violate those principles as neither the PVT lab nor the client pay proper attention to these issues.
Apart from the liberated gas compositions, people working on EoS tuning often utilize the residual oil composition to force the EoS model to reproduce accurately surface fluid densities. Note that almost all PVT data refer to operating conditions (reservoir temperature) and it is difficult for an EoS model to predict accurately surface oil densities and hence represent realistically PVT values such as
or
. Strictly speaking, only two compositions of surface fluids are available in a typical PVT study, namely the flash liquid one and the DL residual one, whereas the latter is obtained by mass balance calculations, as shown in
Section 3. Flash liquid composition is not considered as a good example to work with as it is obtained through an experimental procedure which is not thermodynamically consistent. Indeed, when the reservoir (or separator liquid) is flashed at standard conditions, this is conducted in a non-equilibrium manner since it is only composition recombination that is of interest rather than thermodynamics. As a result, the produced (and measured) flash oil composition cannot be reproduced by an EoS model and it is not recommended to impose its measured density value onto the tuning process. As a result, since it is the residual oil which could possibly be used, it is of the utmost importance that its calculated composition is as accurate as possible. If this composition is based on biased mass balance calculations, owing to inaccurate volumetrics or inaccurate DL gas GC analyses, it must not be incorporated into the EoS model.
Another issue arises from the fact that when incorrect or poorly quality-controlled PVT data are introduced to a reservoir simulator [
31,
32], the fluid model is distorted, thus forcing reservoir engineers to wrongly modify alternate properties (such as the pore volume or the relative permeability curves) to anticipate the deviation in the history-matching process [
33,
34]. Consequently, this leads to numerous incorrect calculations, including estimating recoverable reserves, analyzing fluid flow in reservoirs and wellbores, designing surface pipeline systems and selecting processing equipment. These errors ultimately result in significant errors in production forecasts, adversely affecting both short-term operational decisions and long-term field development strategies [
25,
35,
36,
37]. Finally, apart from validating a PVT report, a QC is also required [
38] when legacy data are collected and handled to form a PVT fluid database. Such databases are quite commonly used as the machine learning era is dominating the market [
39,
40,
41,
42,
43]. In all cases, it is of the utmost importance for the fluids engineer to obtain a deep understanding of the quality of the received data [
44,
45] and eventually map that to the sensitivity of the results obtained by their simulations, future development plans and operational decisions [
46].
In this work, a fully detailed workflow to validate the PVT data delivered in a DL experiment is analytically derived utilizing an exact experimental procedure, lab calculations and the mass balance principle. Subsequently, it is applied to a number of PVT studies and examples of “poor” performance are demonstrated. The proposed workflow can be incorporated into any standard industrial manual or software-based data validation procedure to ensure that only quality-controlled data will be further used for the field development.
The rest of the paper is organized as follows: The experimental work of the DL study is described in
Section 2, while the mass balance and Hoffman equations are developed in
Section 3.
Section 4 demonstrates the application of the proposed workflow over four sample studies to evaluate their accuracy and applicability for the tuning of an EoS model.
Section 5 presents a discussion on the findings and the paper concludes in
Section 6.
2. Differential Liberation Experiment and Data Handling
The differential liberation (DL) experiment, also referred to as the differential vaporization (DV) test, is a standard volumetric PVT experiment conducted in a laboratory [
47,
48]. This test is designed to simulate the depletion process of a black-oil or volatile-oil reservoir, closely approximating the gas–liquid separation process occurring within the hydrocarbon system below its saturation point pressure. Specifically, as the reservoir pressure declines below the saturation pressure, gas dissolved in the oil is liberated [
49,
50]. Once the saturation of the liberated gas reaches critical gas saturation, it begins to flow and separates from the oil that originally contained it. The DL experiment mimics the actual physical phenomenon by enabling the controlled removal of the liberated gas phase from contact with the liquid oil phase under various pressure conditions [
11].
Figure 1 presents a schematic overview of the DL experimental process. The experiment begins by loading a visual PVT cell, which is a high-pressure vessel, with the liquid sample, which is brought to a single, uniform phase at reservoir temperature [
10]. Next, pressure is decreased by moving a piston downward inside the cell until the fluid reaches its saturation point pressure and the oil volume,
, is recorded. Following this, the pressure within the cell is further reduced below the sample’s saturation pressure at reservoir temperature until equilibrium is established, causing the partial vaporization of the liquid [
11,
51]. All liberated gas is then removed from the cell using a needle valve to control the flow rate and directed to a metering device (e.g., gasometer). This process alters the overall composition of the oil sample, and the volume (
), moles (
) and specific gravity (
) of the expelled gas are measured at standard conditions. Note that the composition of the liberated gas is measured using a gas chromatograph (GC). Additionally, the volume of the remaining oil in the cell (
) is recorded. The whole process is repeated step-wisely a number of times at reservoir temperature until pressure arrives at the atmospheric one. Finally, the cell is cooled down to 60 °F, where the residual oil volume (
) and specific gravity (
) are determined [
52,
53]. Finally, the residual oil composition is analyzed using a GC, where it is vaporized into the gas phase before injection into the gas chromatograph.
In addition to the measured data, several properties, directly applicable to reservoir engineering applications, can be extracted at all stages of the DL test by applying mathematical operations on the raw data [
52,
54]. These include the following:
is calculated by dividing the recorded oil volume (
) of the oil at the
th pressure step by the residual oil volume (
):
is obtained by dividing the total amount of gas released from the current pressure step down to the atmospheric one (cumulative gas) at standard conditions by the volume of the residual oil (
):
Note that the subscript
indicates the final stage of the DL test at atmospheric pressure and reservoir temperature.
The oil specific gravity is the ratio of the equilibrium oil’s density at each pressure step to the density of water at 60 °F:
where
is determined from the measured density and volume of the residual oil. The moles of the gas,
, removed at each step, along with the gas molecular weight,
, typically in the form of its specific gravity,
, and the oil volume at each stage, are all measured and known.
The compressibility factor accounts for deviations of the released gas behavior from those of the ideal gas and can be determined by the real gas equation:
where
and
are the pressure and temperature of the hydrocarbon system,
is the volume of the liberated gas at
and
and
,
and
correspond to standard conditions.
From the compressibility factor, the gas formation volume factor (
) can be determined by definition as follows:
The gas specific gravity is an indication of its weight and is used to compute the amount of mass released at each pressure step:
where
is the air molar mass, i.e., 28.97 g/mole.
To relate reservoir oil volumes to produced oil and gas volumes, the differential oil formation volume factor (
) and solution gas–oil ratio (
) are converted to a stock-tank-oil basis according to Equations (7) and (8) [
55].
where
and
denote the oil formation volume factor and the solution gas–oil ratio at saturation point, respectively, as determined from a multi-stage separation test (MSST).
and
are the differential volume factor and solution gas–oil ratio at saturation point. Note that
reflects how much the oil volume expands at reservoir pressure compared to its volume under standard conditions, whereas
is a measure of how much dissolved gas comes out of the oil at a specific pressure.
3. Quality Control of the DL Study Results
To apply quality control (QC) on the differential liberation (DL) study results, one may proceed by considering two major aspects. The first one is related to the mass and mole balance of the products obtained at each pressure step by means of the delivered PVT data, i.e., oil and gas phase properties. This step can be thought of as the inverse calculation of what is supposed to have already been conducted at the lab. Indeed, mass balance is already incorporated in Equation (3), where the oil mass at step
is assumed to equal to the sum of the oil and gas masses obtained at the next step, i.e.,
. Clearly, if mass balance has been honored in the forward calculations, this should also guarantee its validity in the inverse calculations as well [
55].
The second aspect is related to the consistency and quality of the gas compositions obtained at each pressure step, their validity and consistency to the rules governing such compositions. In other words, in this aspect, we focus to the distribution of the gas to individual components rather than only relying on its total mass balance. Both aspects require that various variables, not delivered in the DL report, need to be calculated. Those values can be calculated using equations derived from the principle of conservation of mass and are further applied to a series of fluids with varying volatility. Subsequently, QC can be conducted on those results to determine whether they are consistent or not.
3.1. Total Mass Balance
The primary objective of applying the total mass balance equation is to estimate the density of the residual oil at the end of the DL experiment and compare it to the reported value. If no such value is available, the residual oil density can be evaluated for its physical soundness. In fact, the density of the residual oil is expected to be higher than that of the stock tank oil (STO) obtained after flashing the reservoir fluid to atmospheric conditions through the production separator train. This is due to the fact that the DL study is run in many steps at the reservoir temperature, i.e., at an elevated one, compared to the multi-stage separation test (MSST). On the other hand, during the oil production to atmospheric conditions, the procedure is run in a limited number of separation steps, allowing for some gas to stay dissolved in the produced crude. As a result, more gas is expected to be liberated in the DL test compared to the separator test, thus forcing an increasing residual oil density. When dealing with more volatile oil samples, the difference between these two density measurements is more pronounced, attributable to the higher quantity of dissolved gas present in these samples. Conversely, in the case of heavier oil samples characterized by a lower quantity of dissolved gas, the two density measures tend to converge, resulting in closer calculated values.
By writing down the mass conservation principle, Equation (9) is obtained, in which
stands for the pressure step. This equation asserts that the mass of the sample at the bubble point (
) equals the sum of the gas mass liberated at each pressure step (
) and the oil residual mass (
) received at the final step. The objective is to transform Equation (9) by implementing appropriate formulae based on volumetric data and eliminating the mass terms.
Equations (10)–(12) are used, where
and
are the oil density and volume at saturation pressure,
and
stand for the number of gas moles and their molar mass and
and
denote the density and volume of the residual oil:
By replacing Equations (9)–(11) with Equations (8) and (12), the following is obtained:
Subsequently,
,
and
are replaced by Equations (14)–(16), respectively, where
represents the oil formation volume factor at the bubble point pressure,
stands for the gas volume received at each pressure step at standard conditions and
is the gas molar volume at standard conditions.
is the solution gas-to-oil ratio and
denotes the liberated gas specific gravity.
By simplifying and solving the final expression for
, the residual oil density is obtained by Equation (17):
3.2. Mass Balance per Pressure Step
Setting up a mass balance equation per pressure step bears similarity to the total mass one discussed in the previous section. However, the key distinction lies in the scope of analysis: while the previous equations examined mass balance from the bubble point to the atmospheric pressure, the individual pressure step equations evaluate the conservation of mass at each step of the DL study. Replicating the previous approach, the same equation is initiated, resulting in a paired equation set, which is employed for every step of the study. The equality of the two components of the final equation should persist across each pressure step, as any significant deviation indicates a disruption in the mass conservation.
In this context, Equation (9) is reformulated in Equation (18) and applied to each specific pressure step. Equation (18) indicates that the oil mass at each step equals the combined mass of the liberated gas and the remaining oil:
By applying similar transformations to those of the previous section, the final equation (Equation (19)) is obtained:
By applying Equation (19) using the PVT values at each pressure step, both sides of the equation are anticipated to be equal, thereby ensuring a comprehensive mass balance assessment.
3.3. Component Mass Balance per Pressure Step
In this phase of the QC process, mass balance equations are used to calculate the composition of the remaining oil at each step of the DL study. The molecular composition of the fluid under reservoir conditions, which equals the fluid’s composition at the start of the DL study, along with the molecular composition of the gas released at each stage of the study are both available in the PVT report and are used in the calculations.
Once the oil compositions have been determined, they are plotted on a diagram for each pressure step and for each component. The changes in composition for each pressure value are observed, and any deviations from the expected behavior are analyzed. Under normal circumstances, with decreasing pressure, gas components should show decreasing concentrations, while the heavier oil components should exhibit increasing concentrations. When the pressure reaches the atmospheric level, the concentration of methane and other light components should be very close to zero. Moreover, the residual liquid composition can also be compared to the composition of the flashed liquid, readily accessible in the PVT report. As elucidated in
Section 3.1, it is anticipated that the residual oil extracted from the DL study will exhibit a heavier composition compared to its flashed counterpart.
Similarly to the previous approach, the deduction is initiated from Equation (20), which is also based on the principle of mass conservation. This equation asserts that for each component
, the total moles present in the sample equate to the sum of its moles in the gas phase and those in the remaining oil phase. Through a series of transformations, the ultimate expression is ultimately reached once again. Starting from
the mole expressions can be reformulated as the product of the total moles of each sample multiplied with the molar fraction of each component (
in the liquid phase and
in the gas phase).
The moles of the oil constituents can be subsequently substituted with Equation (22) and the result of these substitutions is shown in Equation (23).
The molecular weight of the oil sample at each pressure step can be computed by multiplying the sum of its oil component’s concentrations and of the sum of their individual molecular weights. The molecular weight of the oil sample in the next pressure step (
), however, cannot be computed in the same way; therefore, it should be computed with Equation (24).
In the following step, by applying similar transformations as those performed in
Section 3.1, Equation (25) is ultimately reached.
By solving Equation (25) for
and further on simplifying the resulting one, the final expression, Equation (27), is obtained:
where
A =
and
B =
.
By solving Equation (27) for at each pressure step and for each component j, the equilibrium oil’s composition is determined in terms of molecular fractions of the components.
3.4. k-Values Distribution and Consistency
The final step of the QC process involves the calculation of the equilibrium coefficients (k-values) for each component at every pressure step. k-values, which represent the equilibrium constants of individual components between the gas and oil phases at the equilibrium conditions, are calculated with the use of Equation (28) [
54,
56,
57]. In Equation (28),
represents the k-value of component
j at pressure step
i of the DL study. Similarly,
and
stand for the concentration of the liberated gas phase and of the remaining oil phase. The values of
were calculated using the method described in
Section 3.3.
Individual k-values are expected to tend towards unity at high pressures, specifically as the pressure approaches the convergence one.
In addition, plots contrasting the values of
for each pressure step against the
values are generated.
s, calculated from Equation (29), remain constant across all pressure steps, while
represents the
logarithm multiplied with the pressure at each step. Note that, in Equation (29),
and
represent the boiling point and the critical temperature of component
, whereas
represents the reservoir temperature.
According to Hoffman et al. [
56,
58], when it comes to hydrocarbon components ranging from C
1 to C
6, the plot of
exhibits a linear trend at all pressures. In addition, lighter components of the C
1 to C
6 range exhibit higher values due to their higher concentrations in the liberated gas phase. These plots need to be generated independently for each sample and for each pressure step of the DL study, comparing the values of the
of hydrocarbon components C
1–C
6 with the constant
values of the same components. As previously mentioned, the plot of
is expected to exhibit a linear trend at all pressure steps, with the
values decreasing for heavier components [
56,
58].
Finally, a combined plot of
can be created, integrating all pressure steps into a single diagram. Since the values of
remain constant in all pressure steps, the combined plot is expected to effectively illustrate the diverse trends of
on the
y-axis in relation to the values of
displayed on the
x-axis. According to Hoffman et al., the values of
are expected to converge at a common point with increasing values of
[
56,
58].
5. Discussion
The analysis of the DL test results for the selected oil samples highlights significant concerns regarding their reliability and accuracy. The study performed covered a wide range of oil volatility and revealed that while DL tests are a fundamental component of a PVT analysis, their results must be examined thoroughly before been utilized in reservoir engineering calculations.
The total mass balance study produced physically sound and reasonable results, regardless of the test fluid volatility, as far as the comparison between the density of the residual oil and that of the stock tank oil (STO) is concerned. Note that the residual oil density can be (indirectly) measured and should appear in every PVT study. However, PVT labs often skip this value to avoid quantitative check calculations and only allow for a rough qualitative test.
On the other hand, the pressure step mass balance results indicated various levels of divergence on the last steps of the process. The deviations were found to be more significant on high-volatile samples. The analysis confirms the previously mentioned findings and shows that the first steps of the test for both volatile and heavy oils are physically sound and can be trusted in general. However, caution is advised regarding the resulting oil composition at the final steps and especially that of the residual oil.
The findings indicate that volatile oils present greater difficulties in DL tests. The sample’s unrealistic results, presented in
Section 4.3.2, underline the challenges of accurately capturing gas liberation and compositional changes due to the rapid and dynamic nature of the volatile oil’s behavior under pressure depletion. This is especially evident in the last step of the test, where the residual oil often ends with a negative mole% methane concentration. Even the volatile sample in
Section 4.3.1, which demonstrated more acceptable results, exhibited inconsistencies that warrant caution, especially when comparing the composition of its residual oil to that of the stock tank oil. This variability could be attributed to the enhanced sensitivity of the volatile oils, especially near-critical ones, to minor errors in lab measurements and data handling, as well as the inherent complexities related to their volatility and the potential failure to maintain thermodynamic equilibrium.
Low-volatility oils exhibit more consistent and reliable test outcomes. The sample presented in
Section 4.3.4, though acceptable in general, still raises the need for careful interpretation to avoid potential fluid modeling risks. On the other hand, the sample presented in
Section 4.3.3, while exhibiting unrealistic behavior in the last steps, is found to be far less problematic compared to its volatile counterpart. The lower gas liberation rates and more stable properties of the low-volatility oils seem to contribute to fewer measurement errors and more trustworthy results. In addition to these, the saturation pressure of typical black oil is significantly lower than that of a volatile fluid. Consequently, fewer pressure steps are required to reach surface conditions, implying less measurements and a decreased risk of measurement errors. This is particularly important because measurement errors, in this context, are mostly systematic rather than random. These errors tend to accumulate rather than cancel out as pressure decreases towards atmospheric levels.
During this study, it was identified that the DL tests seem to be challenged by inaccuracies in GC gas composition measurements, with potential data processing errors and equipment calibration issues being evident as well. These problems highlight the need for extensive quality control (QC) procedures to take place across every step of the process. The necessity for rigorous QC is even more intense when dealing with volatile oils due to the potentially physically unsound results. Engineers should implement thorough QC checks and must validate the results though cross-checking with other PVT data, or even with historical data of the same well and field, if available, to ensure increased accuracy. In the case where discrepancies are detected, they could consider applying corrections to the affected data or even repeating the PVT experiment. Such manipulations can be applied mathematically to estimate the least possible modification which is needed to be applied to the liberated gas composition to honor all mass balance checks. By applying these measures, engineers can enhance the reliability of their fluid models and reservoir simulations, which in turn can lead to more informed and safer field management decisions.
The quality control methodology for DL PVT data developed in this work operates solely as a post-test analysis tool rather than a real-time QC tool. This is because the values for differential oil formation volume factor (Bo), differential solution gas–oil ratio (Rs), oil density (ρo) and gas specific gravity (γg) required for the material balance equation cannot be determined until the DL test is completed and the density of the stock tank oil is obtained, which is only defined at the final stage of the test. Consequently, the mass balance cannot be applied during the DL test to identify errors early; it can only be utilized once the test is complete. Once the raw measurements have been converted to PVT data through mathematical manipulations, the quality control procedures outlined in our study can be conducted. If quality control identifies unreliable PVT data at specific pressure steps, two solutions are possible: either repeating the test, which is both time consuming and expensive, or computationally correcting the volumes at these steps to ensure that the material balance is satisfied.
Future research could broaden the scope of this study by developing a data analysis technique that addresses measurement errors computationally. This technique would aim to identify the minimum adjustments to measured values needed to ensure all mass balances within the system are satisfied. Given the redundancy in the system (more unknowns than equations), the problem is inherently ill posed. However, incorporating specific constraints into the optimization problem could lead to a robust solution. One example may be implementing a weighting scheme that prioritizes adjustments to measurements from lower-pressure steps, where errors are known to be more prominent. Additionally, the framework could be expanded to incorporate other types of process constraints, leading to a more robust error mitigation strategy. This approach has the potential to significantly improve data analysis accuracy and reliability in various fields. By computationally accounting for measurement errors while maintaining mass balance integrity, researchers can extract more reliable information from their experiments.
6. Conclusions
The methodology for quality checking the results of a typical DL test analysis through mass balance verification was demonstrated on a selection of fluid samples with varying volatility. The comprehensive study that was conducted revealed some critical insights into the accuracy and reliability of the reported results. It further indicated that there was a general concern regarding the results of the DL test. While generally the results of the first steps can be trusted and considered rather safe, the later steps and especially the final step result in unreliable data. Volatile oils demonstrate more significant inconsistencies than their heavier counterparts. These inconsistencies, such as the unrealistic negative concentration of light components, highlight the challenges of properly and accurately capturing the composition of liberated gas and the compositional changes in each step.
Low-volatility oils exhibit more consistent and reliable results. The less intense gas liberation contributes to less error accumulation during subsequent pressure steps. However, issues may appear at the later steps, and although less severe than those in the volatile oils, the final step was still maybe plagued by inconsistency.
The study manages to underline the necessity for extensive and rigorous quality control (QC) procedures at every step of the DL process, especially for volatile oil, and engineers should not rely solely on DL results without performing checks and applying corrections when and if needed. Potential inaccuracy in gas composition measurements, data processing errors, equipment calibration or simply human errors are possible factors of error and could lead to highly inaccurate data, which could damage the performance and the accuracy of an EoS model, that of a reservoir simulation, as well as that of machine learning model developments based on these data. The proposed workflow for validating PVT data can be integrated into standard industrial procedures, ensuring that only quality-controlled data are used for field development and operational planning.
In future work, the developed workflow can be extended to the quality check of multi-stage separation tests (MSSTs) and constant volume depletion (CVD) PVT results. Step mass balance control can be applied to both sets of experimental PVT data, whereas the total mass balance cannot be applied to the CVD data as the residual fluid is discarded.