Geofluids

Volume 2017, Article ID 5395308, 8 pages

https://doi.org/10.1155/2017/5395308

## A New Dynamic Model for Sealed Coring Saturation Correction in Hydrocarbon Reservoir

^{1}China University of Petroleum, Beijing, China^{2}MOE Key Laboratory of Petroleum Engineering, China University of Petroleum, Beijing 102249, China

Correspondence should be addressed to Bin Tu; moc.361@tubin

Received 20 June 2017; Accepted 29 August 2017; Published 9 October 2017

Academic Editor: Zhenhua Rui

Copyright © 2017 Bin Tu and Jie Li. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

It is one of the most intuitive methods to analyze the formation of oil and water saturation by sealed coring. But this method is affected by a variety of factors, such as pore volume change, fluid compression, and separation of dissolved gas. In view of the impact of such factors on sealed coring test saturation, there are four kinds of correlation methods currently, namely, comprehensive coefficient correction method, experiment correlation method, mathematical statistics method, and mathematical model correction method, with their own advantages and disadvantages. Based on the analysis of physical change during core lifting process, this paper proposes the mathematical model for dynamic correction of core saturation based on dissolved gas drive filtration theory, multiphase flow fractional flow theory, and corresponding work flow. This method comprehensively considers the impact of relative permeability of three-phase fluid flow, elastic compression nature of fluid and rock, fluid viscosity, volume factor, solution GOR, and other high-pressure PVTs, so it features a wider adaptability, and the accuracy of model correction results satisfies the project requirements. This method provides a reliable basis for the true oil-water saturation of actual reservoir and has an important theoretical and practical significance.

#### 1. Introduction

The analysis of remaining oil saturation is the basis and foundation for water controlled reservoir potential evaluation and development plan adjustment.

There are numerous methods to calculate the remaining oil saturation of formation. The methods of indirect acquisition of the formation saturation, such as the logging interpretation method, have the problems of being affected by casing cement, casing collar, and other factors and having a lower resolution, and meanwhile they contain measurement errors and interpretation errors [1]. The method of direct acquisition of formation saturation is to obtain the formation rock by sealed coring and then obtain the remaining oil saturation [2] by test. Although this method requires a higher cost, the results obtained by such method are relatively accurate, intuitive, and reliable, so this method is utilized widely in China and provides reliable data of petrophysical properties [3, 4] for exploration and development and study and test of oilfield. But, the sealed coring method also has obvious systematic errors, and the actual formation fluid saturation [3] may be only obtained by saturation correction. Currently, there are four kinds of correction methods, namely, laboratory experiment correction method [5–9], correction coefficient method [10], mathematical statistics regression method [11, 12], and mathematical model correction method [13–17]. The indoor laboratory experiment method does not need to directly describe the complex processes of physical and chemical changes during the depressurization degassing process, with the results relatively accurate. But this method requires a heavy workload, the experiment results are only applicable to specific reservoirs, and the fitting results have no universality [5–9]. The mathematical statistics method is simple and easy to understand, but it considers neither the differences of core and fluid properties nor the state change in the coring process. The damage rates obtained from statistics are affected by different rock properties; the classification standards of sample points are not objective and contain larger errors; the correction coefficient method considers few factors and is quite simple, but it cannot obtain the accurate formation saturation results. In summary, the application of all the three methods is limited seriously. Only the mathematical model method has caused the extensive concern because it can describe the changes of core and fluids during the coring process by mathematical equations. The existing mathematical model correction methods [13–17] all consider the changes of rock pore volume and fluid volume and determine the changes of fluids in the coring process based on the fractional flow equation, but the established models have two problems: one is that the overall process coupling correction is not conducted for the changes of rock volume and changes of fluids; and the second is that the dissolution and separation of oil and gas are not considered in the overall coring process. In view of the above problems, this paper plans to establish a kind of method which considers both the rock deformation and separation of oil and gas in the overall process, so as to obtain a more accurate dynamic correction model for sealed coring saturation.

#### 2. Existing Correction Methods for Sealed Coring Saturation

From the published literature, there are currently four kinds of correction methods about pressure impacting saturation, such as depressurization degassing.

*(1) Comprehensive Coefficient Correction Method*. Kairong and Wanshun [4] and Rahtmell et al. [5] propose the following correction method: of which parameter is the comprehensive fluid contraction coefficient and refers to the *Parameter List* for the remaining parameters. The research of Kairong and Wanshun and Rahtmell et al. suggests that is within 1.10~1.15, and the research of Kazemi [6] suggests that is 1.062. It can be seen that this parameter has a wide value range, so this method is quite simple, but it cannot get accurate formation saturation results.

*(2) Experiment Correlation Method*. Shenglai et al. [7], Liang [8], and Liu [9] utilize the experiment method for correction of sealed coring saturation. The core of this method is to clean the core and inject a certain amount of simulated formation oil after ensuring that it is saturated with water and then simulate the drop of core barrel pressure to obtain the correlation formula of water saturation before core degassing and saturation after degassing: We can obtain coefficient and coefficient through regression; namely, we obtain the correction equation of water saturation of one reservoir. This method may not consider the complex processes of physical and chemical changes during the depressurization degassing process but has quite obvious disadvantages: ① experiment results are only applicable to specific reservoirs, and the fitting results have no universality; ② the simulated formation oil contains errors; ③ the experiment process contains errors.

*(3) Mathematical Statistics Method*. Egbogah and Amaefule [10] and Yijing et al. [11] consider that oil-water saturation losses of similar types of rocks should be basically the same, so the measured saturation of a large number of cores should satisfy the relation: of which stands for the loss rate of crude oil saturation; stands for the loss rate of formation water saturation; we can see that formula (3) is the equation of one line, so we can obtain the loss rate of crude oil saturation and that of formation water saturation as coefficients in this equation_{.}Xin et al. [13] and Huiming et al. [14] further discover that different correction coefficients can be obtained for different types of rocks. This method is simple and easy to understand, but the loss rate is impacted by different rock properties, and the classification standards of sample points are not objective and contain larger errors.

*(4) Static Zero-Dimensional Model Correction Method*. Mathematical models proposed by Xin et al. [13], Huiming et al. [14], and Yuhuan [15] are basically consistent and all establish the zero-dimensional mathematical model based on fractional flow equation and mass conservation in the basic theory of “the oil-gas-water is allocated according to fractional flow rate in the total volume of discharged total volume,” and consequently the derived correction equation for oil saturation is

There is an important assumed condition in these models; namely, it is considered that the fluid saturation in the core is completed instantaneously when the core is lifted to the ground, without considering that the precipitation and escape of dissolved gas in the lifting process are a changeable process and even without considering the fluid escapes from the core due to dissolved gas drive in the process.

#### 3. Establishment of Dynamic Correction Model for Sealed Coring Saturation

##### 3.1. Affecting Factors for Sealed Coring Saturation

A large number of study results show that the oil-water saturation measured by sealed coring are affected by the following 5 aspects [3–10]: ① rock pore volume change and fluid volume change; ② depressurization degassing of crude oil in core, and appearance of dissolved gas drive; ③ volatilization of oil and water in the transport process after the core is lifted to the ground; ④ mud invasion in drilling coring; ⑤ experiment error.

In the 5 factors, experiment errors are omnipresent, and the research of Xin et al. [13–15] suggests that the mud invasion and core transport process are not very important factors, so they are not key considerations in this paper. The combination of multiple researches indicates that the main affecting factors of core saturation of sealed coring are rock pore volume change, fluid volume change, and dissolved gas drive [1–10], so the correction of the three kinds of affecting factors is also the key problem focused by numerous researches.

##### 3.2. Dynamic Variation of Oil-Water Saturation in Coring Process

The well bottom is generally in the HTHP state in sealed coring. When the core is lifted to the ground, the change range of temperature in the core barrel is small due to a short time, while the pressure change is quite obvious, so the pressure change is the main reason to cause fluid saturation change in the core.

In the lifting process, the pressure continues to decline, and the fluid saturation in the core also changes constantly, which is a dynamic continuous process. Based on the discrete principle, we can utilize a finite number of stepped processes to approximately represent one continuous physical process; namely, we can replace the continuous pressure change in the core lifting process with stepped pressure, and the fluid in the core can have the following changes for reduction of each time:

① After pressure changes , if the pressure in the core barrel is higher than the saturation pressure, only the core pore volume and liquid volume in the core change. Under the combined action of the two, the liquid in the core is squeezed out, causing the changes of fluid saturation in the core. When the pressure drop is small, the oil-water volume squeezed out may be considered to be distributed according to the oil-water fractional flow [13–15].

② If the pressure in the core barrel is lower than the saturation pressure after the pressure changes , the dissolved gas separation phenomenon appears. So, besides the expansion of liquids and reduction of pore volume, the dissolved gas also expands and the volume change range is wider, and consequently, the fluid volume squeezed out is larger, and the change range of fluid volume in the core is also wider.

With the continuous lifting of core barrel, the pressure in the barrel reduces further. The saturation change in the core barrel may repeat the above ① or ② process, and when the pressure is much lower than the saturation pressure, the gas precipitation amount increases and becomes continuous phase. So most of the gas escapes from the core and carries liquid flowing out of the core, showing obvious characteristics of dissolved gas drive.

##### 3.3. Theoretical Basis and Computation Process for Dynamic Correction Model

###### 3.3.1. Theoretical Basis for Dynamic Model

The process of oil-gas two-phase fluid flow or oil-gas-water three-phase fluid flow is unstable, but every instant in the overall process can be approximately regarded as a stable state. In this way, the unstable state of the overall process can be treated as the superposition of countless stable states, and this method is called stable state successive displacement method [15]. When the value taken for time interval or pressure is very small, the results obtained by this method basically comply with the actual situation.

###### 3.3.2. Dynamic Correction Equation for Two-Phase Fluid

When the pressure declines , if the dissolved gas is not separated, the total fluids squeezed out due to expansion of liquids and rock pore volume compression cause is

According to the material balance principle,* residual crude oil amount after core depressurization = crude oil amount before depressurization − crude oil amount squeezed out*, and it is considered that the oil-water amount squeezed out is distributed according to the fractional flow equation, so we can get

Simplify and substitute (5) into the above equation: Of which, stands for water cut, and according to fractional flow equation: of which and are a function of saturation, and the corresponding saturation should be the average saturation before and after depressurization. Considering there are only oil-water two phases, .

Thus, (7) is an implicit function equation of saturation , and we can obtain with the iteration method.

###### 3.3.3. Correction Equation for Three-Phase Flow Saturation

When the pressure declines , if the dissolved gas is separated, the total fluids squeezed out due to expansion of liquids and rock pore volume compression cause is

When depressurization is very small, the seepage in the core can be approximately considered as steady seepage, so the oil-gas-water three phases squeezed out of the core may be approximately distributed according to fractional flows , , and .

According to the material balance principle,* residual oil amount after core depressurization = oil amount before core depressurization − discharged oil amount*.

In a similar way, we can get the correction equation of water saturation:

According to the material balance principle,* free gas quantity after depressurization = free gas quantity before depressurization + dissolved gas amount before depressurization − dissolved gas amount after depressurization − discharged gas amount*, so we can get the correction equation for gas saturation:

In the above three correction equations, the expressions for the fractional flows , , and for the oil, water, and gas are, respectively,

###### 3.3.4. Calculation Flow of Dynamic Model

According to the process analysis of oil-water saturation changes in the core lifting process, we can utilize the flow as shown in Figure 1 to fit the saturation of core test: