Research on Operation Optimization of Fluid Sampling in Wireline Formation Testing with Finite Volume Method

Preprint

Article

Research on Operation Optimization of Fluid Sampling in Wireline Formation Testing with Finite Volume Method

Altmetrics

Downloads

Views

Comments

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

18 June 2024

Posted:

18 June 2024

You are already at the latest version

Alerts

Abstract

Wireline Formation Testing is an important technique in the oil field of exploration and development. Not only it can make real fluid samples of the formation directly obtained to know exactly oil existed in formation or not; but it can also show flowing pressure change to see the production capacity of the formation. So it is important measurement for formation evaluation during the drilling process and supports for activities of the exploration and development in the oil field. A numerical simulation model in this article is researched and established based on the Finite Volume Method with considering the influence of sensitive parameters such as reservoir heterogeneity, probe suction area, and mud-filtrate invasion depth during the drilling. The model is capable to design and evaluate for formation fluid sampling operations with calculating hydrocarbon content and flowing pressure. Furthermore, through case application, the performance and effect in the process of Wireline Formation Testing was demonstrated. The results indicate that this technology can be served as an effective auxiliary tool for fluid sampling operations with the function of optimizing fluid sampling measures. It can improve the accuracy of predicting indicators such as hydrocarbon content and breakthrough time during the sampling process. This study provides important backups and technical guidance for professionals in geological exploration and oilfield development.

Keywords:

Subject: Engineering - Energy and Fuel Technology

1. Introduction

Wireline Formation Testing (WFT) plays an important role in exploration and development of oil and gas field, which is a logging method for fluid sampling, reservoir analysis, and fluid analysis. During drilling process, not only can it provide important basis for comprehensive geological analysis, but also it can provide necessary data for in-depth research on oil and gas reservoirs through fluid sampling [1]. WFT can complete formation fluid sampling, reservoir pressure analysis, and formation pressure gradient testing [2], that can combine with test data for calculating physical parameters such as formation permeability and fluid density [3]. Thus, it can directly, quickly, and accurately distinguish the fluid properties of the formation, assist in comprehensive reservoir evaluation, determine the oil-water interface, and interpret permeability and evaluate production capacity [4].

Therefore, a large number of scholars have also conducted related research on the optimization of fluid sampling operations in the process of Wireline formation testing. In 2007, a finite element method has been researched for the numerical simulation of WFT, which can describe the boundary condition of variable storage pipeline volume. Through the application of numerical simulation technology, it is important to providing theoretical guide for the parameter optimization design of WFT [5]. In order to solve the limitations of calculating pressure gradient from logging and sampling data with WFT due to the uncertainty of the method, an optimization method for the computational precision of pressure gradient was studied by numerical simulation, which can be applied in the analysis of the fluid property, the vertical connectivity of reservoirs, the vertical change of formation fluid and the assistant evaluation of reservoirs [6]. An inversion mothed was researched by Faruk O. et al. in 2008, which can analyze the in-situ petro physical properties of hydrocarbon bearing formations. Implemented the inversion process of rock physical parameters，a fully implicit finite-difference black-oil reservoir simulator is used to simulate fluid-flow phenomena with consideration of multi-physics data [7].

Mathieu N. et al. researched a new hypothesis during the WFT，investigated volumes can be related to the heterogeneity of the permeability field. In addition, a flow model was implemented in heterogeneous and anisotropic synthetic reservoirs, which not only explains the difference of drawdown and buildup mobility values inferred from buildup analysis, but also evaluates the impact of formation permeability near the WFT tool [8]. The sampling simulation method was studied by an extremely high precision grid model to optimize dynamic formation tester operation, which can confirm the main parameters affecting of fluid breakthrough time and predict fluid sampling purity, thereby forming fluid breakthrough volume prediction chart in the producing process. It is improved to operation efficiency [9]. Quan Z. et al. simulated the formation pressure response curve using numerical simulation technology, combined with different probe parameters of WFT instrument, pumping speed, geological features and fluid properties. From the results of the study, which can guide the selection of sampling probes and the setting of instrument parameters during the operation process [10]. On the whole, numerical simulation is an important part of fluid sampling analysis in WFT, and it is plays an important role in calculate hydrocarbon purity, formation permeability, pressure buildup curve and PVT parameters, the research results of a large number of scholars have provided important assistance for optimizing the design of WFT operations during drilling process.

However, many studies have tended to overlook the complex conditions of reservoir heterogeneity and probe suction area, specifically the impact of mud-filtrate invasion depth. This lack of consideration has resulted in a non-systematic understanding of fluid sampling during the WFT process, resulting in the optimization design of production measures and the prediction of key parameters will face many of problems.

In this paper, the technology of numerical simulation was used to establish a finite volume model for the operation optimization of fluid sampling in the process of WFT. Combining the characteristics information of WFT tools and matching probes, which can calculate the hydrocarbon content and flow pressure during the sampling process. Meanwhile, according to analyzing the pumping efficiency of different types of probes under the different geological conditions, it can provide a basis for optimizing fluid sampling operations. Finally, the practicality of the method is demonstrated by example application.

2. Materials and Methods

2.1. Introduction of Pumping Probe

WFT is one of the most direct and effective methods to confirm reservoir fluid properties and evaluate reservoirs in oil field exploration operations, playing an increasingly important role [11]. The Enhanced Formation Dynamic Tester (EFDT) is a modular, pump-extraction WFT instrument [12], designed for exploration and evaluation wells to accomplish formation fluid sampling tasks. It provides various probe options for sampling operations, allowing for the construction of different probe combinations to establish a connection between the internal pipeline of the formation tester and the formation fluid. This facilitates pump-extraction sampling and pressure measurement, hydrocarbon purity calculation, as well as analysis of reservoir fluid properties, fluid interfaces, and reservoir properties [13]. Therefore, under the same instrument configuration conditions, the efficiency of pump-extracted sampling varies depending on the shape and distribution of the probes [14]. This study mainly focuses on probes such as the Small Type Inlet Probe, Middle Type Inlet Probe, Ellipse Type Inlet Probe, and Large Type Inlet Probe of the EFDT. Relevant parameters are shown in Table 1.

The standard ratio refers to the ratio of the suction area of other probes to that of the middle-type inlet probe. In addition, the 3D probe is achieved by combining a fluid sampling pump with greater efficiency. Based on the Small Type Inlet Probe, Middle Type Inlet Probe, Ellipse Type Inlet Probe, and Large Type Inlet Probe, three probes of the same type are selected to achieve the application of a 3D probe with a uniformly distributed radial orientation of 120 degrees and a combination of radial push-and-connect.

2.2. Meshing Model

Grid model generation technology serves as a link between physical and computational models for numerical simulation. Depending on whether the generated grid units have regularity, it can be divided into structured grids and unstructured grids [15]. The quality of grid generation has a significant impact on factors such as accuracy and efficiency in subsequent numerical simulation calculations, playing a crucial role in the numerical simulation process [16]. This paper primarily focuses on the parameters such as the length and width of the probe suction mouth used in fluid pump-extraction sampling operations during WFT. Centered around the target operation well, a method for constructing a radial grid model with variable grid step length is developed. This method provides the spatial coordinates of each grid node and the connectivity of each grid node, supporting the implementation of the numerical simulation calculation model and providing a more detailed description of the operational process.

Therefore, in the process of dividing the three-dimensional spatial grid model, to simplify the calculation formula, it is assumed that the X and Y coordinates at the position of the wellbore center are both zero. The calculation formulas for the X and Y coordinates of the vertices of the grid model along the perimeter of the wellbore are shown in equations (1) and (2).

x (i, j, k) = x^{'} c o s (j α - 0.5 α)

(1)

y (i, j, k) = x^{'} s i n (j α - 0.5 α)

(2)

Where

x^{'} = (r e + \sum (n (i i) - 1) r_{s} (i i)) c o s (0.5 α)

(3)

α = 360 / n r

(4)

R = (r e + \sum n (i i) r_{s} (i i))

(5)

In the above formulas,

i

denotes the azimuthal node number,

i = 1,2 \dots, n c

，with

n c - 1

being the number of azimuthal grids;

j

denotes the radial node number, with

n r - 1

being the number of radial grids;

k

denotes the longitudinal node number,

k = 1,2, \dots, n k, \dots, 2 n k + 2

, with

2 n k + 1

being the number of longitudinal grids;

α

represents the azimuthal partition baseline angle;

r e

denotes the wellbore radius;

i i

represents the sequential number of the planar grid partition region;

n

denotes the number of partitions corresponding to each region on the plane;

r_{s}

represents the radial step length for each region;

R

represents the modeling control radius with the wellbore as the center on the plane.

Furthermore, based on the depth D where the probe is located, assuming the baseline node number of the longitudinal direction where the probe is located is used as the reference node,

k = n k

; The formula for calculating the longitudinal Z coordinate of the grid vertex with variable longitudinal step length is as follows:

z (i, j, k) = D - (2 (n k - k + 1) - 1) h - ∆ D (n k - k + 1)

(6)

In equations (5),

k = 1, \dots, n k + 1

D_{T} = z (i, j, 1)

(7)

To construct grid nodes downward and achieve the calculation of the longitudinal Z coordinate of grid vertices with variable longitudinal step length, the formula is as follows:

z (i, j, k) = D - (2 (n k - k + 1) - 1) h - ∆ D (n k - k + 1)

(8)

In equations (6),

k = n k + 2, \dots, 2 n k + 2

D_{B} = z (i, j, 2 n k + 2)

(9)

In the above formulas,

D

represents the depth of the probe operation,

h

denotes the reference grid thickness centered around the depth of the probe operation,

D_{T}

represents the depth of the top of the target reservoir, and

D_{B}

represents the depth of the bottom of the target reservoir. Based on the grid partitioning methods described above, both fixed-step and variable-step grid models are constructed, and the computed results are as follows:

In Figure 1,

R = 3

r e = 0.15

n c = 31

n r = 10

h = 0.05

∆ D = 0

n k = 18

. In Figure 2,

R = 3

r e = 0.15

n c = 31

n r = 10

h = 0.05

∆ D = 0.02

n k = 7

Through the analysis of the two cases above, it is evident that the total number of grids generated by the grid model shown in Figure 1 is 9*30*37, while the total number of grids generated by the grid model shown in Figure 2 is 9*30*15. Compared to the fixed-step grid partitioning method, by controlling the relevant parameters of the variable step, not only can the control of local grid volume be achieved, but also the total number of grid model nodes can be reduced more effectively for subsequent numerical simulation calculations, thereby improving computational efficiency. Through controlling the step length of the planar partition area, overlaying longitudinal step parameters, and controlling the width of the probe suction mouth, it is more convenient to establish a three-dimensional wellbore reservoir grid model considering information such as the probe suction area, wellbore radius, and mud invasion range.

2.3. Numerical Simulation Model

In the field of fluid mechanics, finite difference, finite element, and finite volume methods are commonly used numerical techniques for solving partial differential equations [17]. The basic principle of the finite volume method (FVM) is to discretize the computational domain of the physical phenomenon described by the differential equations into discrete finite volume elements. Then, the integral of properties such as energy and mass is performed over each volume element to construct a discrete system of algebraic equations. The numerical solution is obtained by solving this system of equations [18]. FVM directly discretizes integral forms of conservation equations in physical space, making it suitable for arbitrary grid forms. One significant advantage is its connection to the concept of conservation, automatically satisfying the discretization with conservatively. Therefore, in the context of numerical simulation of fluid sampling in wireline formation testing, considering the heterogeneity of reservoirs around the well, centimeter-scale grid resolution, and irregular grid models, the integral equation for mass conservation is formulated as follows:

\underset{s_{e} + s_{i}}{∯} (\sum V) d s + \sum Q = \underset{v}{∭} ∆ (\emptyset \sum (S ∙ ρ)) d v

(7)

In these equations,

\underset{s_{e} + s_{i}}{∯} (\sum V) d s

represents the surface flux integral over the control volume,

\sum Q

represents the volume integral of the source/sink terms, and

\underset{v}{∭} ∆ (\emptyset \sum (S ∙ ρ)) d v

represents the volume integral of the mass increment within the control volume.

To better understand the seepage during the pump extraction process, this study makes the following assumptions regarding the seepage of fluid into the probe during the wireline formation testing fluid sampling process, focusing on three main stages: pure water phase, oil-water mixed phase, and near-pure oil phase:

(a)Three-Dimensional Two-Phase Black Oil Model: In this model, the invading formation mud and formation water are collectively referred to as the water phase, while the formation oil is the hydrocarbon phase. The two phases are immiscible and do not undergo chemical reactions.

(b)Consideration of Various Factors: This includes the heterogeneity of the reservoir in both the horizontal and vertical planes, capillary pressure, gravity effects, rock compressibility, and fluid compressibility.

(c)Boundary Conditions: The invading mud around the wellbore is water-based, and the initial invasion is an editable known condition. The outer boundary condition of the reservoir is a constant pressure boundary, while the inner boundary condition of the probe fluid sampling model is a constant liquid operation.

Using the finite volume method for spatial discretization and the backward Euler method for temporal discretization, the mass conservation discretized control equation with pressure as the solving parameter is formulated as follows:

\sum_{j = 1}^{z (i)} T_{i j} (ρ_{w i j} λ_{o j}^{n} + ρ_{o i j} λ_{w j}^{n}) P_{h (i, j)}^{n + 1} - (\sum_{j = 1}^{z (i)} T_{i j} (ρ_{w i j} λ_{o j}^{n} + ρ_{o i j} λ_{w j}^{n}) + \frac{V ϕ_{i} ρ_{o i} ρ_{w i} (c_{o} S_{o}^{n} + c_{w} S_{w}^{n} + c_{f})}{Δ t}) P_{i}^{n + 1} = - g (\sum_{j = 1}^{z (i)} T_{i j} Δ D_{i j} ρ_{w i j} ρ_{o i j} (λ_{o j} + λ_{w j})) + \sum_{j = 1}^{z (i)} T_{i j} ρ_{o i j} λ_{w j}^{n} (P_{c o w j}^{n} - P_{c o w i}^{n}) - (\frac{V ϕ_{i} ρ_{o i} ρ_{w i} (c_{o} S_{o}^{n} + c_{w} S_{w}^{n} + c_{f}) P_{i}^{n}}{Δ t}) - Q

(8)

In the above equation:

T

is the transmissibility between the connected grid blocks,

T = \frac{1}{T_{i} + T_{j}}

T_{i} = \frac{L}{A ∙ N T G ∙ K ∙ d}

;

A

is the contact area of the current grid block in different connected directions;

L

is the distance from the center of the current grid block to different contact surfaces;

N T G

is the effective thickness ratio;

K

is the permeability;

d

is the probe area ratio coefficient, where

d = 1

for grid blocks not in contact with the probe inlet, and

d = \frac{A_{p}}{A_{c}}

for grid blocks in contact.

A_{p}

is the probe inlet area,

A_{c}

is the contact area of the grid block in operation;

n

is the time-step coefficient;

i, j

are the indices of grid blocks with a connection relationship;

z (i)

is the total number of grid connections for the current grid block

i

;

λ_{o} = \frac{k_{r o}}{μ_{o}} ρ_{o}

，

λ_{w} = \frac{k_{r w}}{μ_{w}} ρ_{w}

k_{r o}

is the relative permeability of the hydrocarbon phase;

k_{r w}

is the relative permeability of the water phase;

ρ_{o}

is the density of the hydrocarbon phase fluid;

ρ_{w}

is the density of the water phase fluid;

μ_{o}

is the viscosity of the hydrocarbon phase fluid;

μ_{w}

is the viscosity of the water phase fluid;

ϕ

is the porosity;

P

is the pressure;

V

is the grid block volume;

c_{o}

is the compressibility of the hydrocarbon phase fluid;

c_{w}

is the compressibility of the water phase fluid;

c_{f}

is the comprehensive compressibility;

Δ D

is the depth difference between grid blocks;

Δ t

is the iterative time step size;

g

is the gravitational acceleration;

Q

is the total sampling fluid rate at the probe, with

Q = 0

in the pressure equation for grid blocks not in contact with the probe.

The water phase saturation model for solving the above equation is given by:

S_{w i}^{n + 1} = S_{w i}^{n} + \frac{Δ t (\sum_{j = 1}^{z (i)} T_{i j} λ_{w i j} (P_{w j}^{n + 1} - P_{w i}^{n + 1} + ρ_{w i j} g Δ D_{i j}) + q_{w}) - V φ_{i} S_{w i}^{n} Δ P_{w i}^{n + 1} ρ_{w i} (c_{w} + c_{f})}{V φ_{i} ρ_{w i}}

(9)

S_{o i}^{n + 1} = 1 - S_{w i}^{n + 1}

(10)

Where,

S_{w}

is the water phase saturation and

S_{o}

is the hydrocarbon phase saturation.

Other major auxiliary equations include:

Relative permeability model:

K_{r w} = K_{r w} (S_{w})

(11)

K_{r o} = K_{r o} (S_{w})

(12)

Capillary pressure model:

P_{c o w} = P_{o} (S_{w}) - P_{w} (S_{w})

(13)

Mud invasion into the formation:

S_{w} = {S (x, y, z)}_{t = 0}

(14)

The effective thickness ratio of the formation:

N T G = {N T G (x, y, z)}_{t = 0}

(15)

The porosity of the formation:

ϕ = {ϕ (x, y, z)}_{t = 0}

(16)

The permeability of the formation:

K = {K (x, y, z)}_{t = 0}

(17)

Where,

P_{c o w}

is the capillary pressure,

P_{o}

is the hydrocarbon phase pressure,

P_{w}

is the water phase pressure,

N T G

is the effective thickness ratio, and

ϕ

is the porosity.

The main advantage of the finite volume method is its good adaptability to irregular and complex seepage simulation calculations, as well as its ability to handle complex physical phenomena such as phase change and multiphase flow occurring within the computational domain. Firstly, the computational domain can be divided into several non-overlapping control volumes or flow cells. Secondly, the mass conservation equations for each control volume or flow cell are discretized. Finally, using the discretized forms of the equations obtained, in conjunction with matrix solution methods, the numerical solution of the discretized equations for reservoir water saturation and formation pressure distribution during the fluid sampling operation process is computed.

2.4. Target Parameter Calculation Model

The calculation of hydrocarbon content during the wireline formation testing process while drilling plays a crucial role in reservoir analysis, productivity assessment, and implementation of operational measures. Therefore, to obtain the hydrocarbon content and flow pressure during the wireline formation testing process while drilling, considering factors such as skin factor, probe-to-grid contact area ratio coefficient, this paper integrates the pump rate. Based on spherical flow, the model for calculating flow pressure and hydrocarbon content at the probe inlet is constructed as follows:

The calculation model for hydrocarbon content in the fluid at the probe inlet:

f_{w} = \frac{1}{1 + \frac{u_{w}}{u_{o}} \cdot \frac{k_{r o}}{k_{r w}}}

(18)

q_{w} = Q \times f_{w}

(19)

q_{o} = Q - q_{w}

(20)

f_{o} = \frac{q_{o}}{Q}

(21)

The calculation model for the flow pressure at the probe inlet:

P_{w f} = \frac{q_{o}}{G_{o}} + P_{o}

(22)

The

G_{o}

solution model is as follows:

G_{o} = \frac{2 π {K k}_{r o} h d}{μ_{o} (l n (\frac{r_{e}}{r_{w}}) + s)}

(23)

Where,

f_{w}

represents the water phase fluid content ratio;

f_{o}

represents the hydrocarbon phase fluid content ratio;

q_{w}

is the flow rate of the water phase fluid;

q_{o}

is the flow rate of the hydrocarbon phase fluid;

P_{w f}

is the flowing pressure at the probe inlet;

G_{o}

is the comprehensive characterization coefficient;

r_{e}

is the wellbore radius;

r_{w}

is the effective fluid supply radius of the wellbore;

s

is the probe skin factor, which is dimensionless.

2.5. Solving Matrix Equation

The Stabilized BiConjugate Gradient (SBiCG) method is an iterative technique used to solve linear systems of equations with asymmetric matrices (Ax=b). In comparison to the conventional BiConjugate Gradient (BiCG) method, SBiCG introduces additional step parameters and a modified inverse preconditioner, leading to faster convergence rates and improved numerical stability [19]. Therefore, this paper primarily adopts the matrix pre-processing SBiCG solution method based on incomplete LU decomposition for solving pressure equations and saturation. The main steps of the solution process are outlined as shown in Figure 3.

The basic idea of the SBiCG iterative method is to maintain the orthogonality of vectors during the iteration process while using appropriate step parameters to enhance the convergence speed. This method first computes two conjugate directions and then uses them to generate two new search directions. These directions remain orthogonal to previous search directions, and new iterative solutions are generated using step parameters.

It also employs a modified inverse preconditioner. The modified inverse preconditioner utilizes an additional step parameter to balance the effects of preconditioning and incomplete LU decomposition. This helps prevent the accumulation of rounding errors in the preconditioner during the iteration process and aids in maintaining the orthogonality of search directions.

In the numerical simulation model solving process, SBiCG is an efficient method for solving asymmetric matrix linear systems of equations. It can maintain the orthogonality of search directions while achieving good numerical stability and rapid convergence speed. By characterizing the fluid seepage characteristics of the grid model and utilizing the computed pressure field values at each time step, combined with equations (9) and (10) mentioned above, the calculation of reservoir water and hydrocarbon saturation distributions can be further completed. Through calculation models such as phase permeability, hydrocarbon content, and flow pressure, the calculation of hydrocarbon phase fluid proportion and pressure is accomplished, providing important technical support for the design of WFT fluid sampling operations, process optimization analysis, and indicator prediction during drilling operations.

2.6. Model Analysis

To analyze the computational effectiveness of the WFT fluid sampling process numerical simulation model, numerical simulations were performed considering the pump-extraction sampling process of an elliptical probe with water-based drilling fluid invading the reservoir. A comprehensive analysis of the method's applicability was conducted by fitting and comparing the simulation results with measured data. The relevant basic parameters are shown in Table 2 and Table 3, and the experimental data for the relative permeability curves are presented in Table 4.

By considering the starting depth of 3286.8m and the ending depth of 3287.8m, with the probe located at a depth of 3286.9m, the grid longitudinal thickness is set to 0.026m and the longitudinal overlapping step is 0.005m. The planar modeling radius is 3m, with 37 radial grid divisions. The plane is divided into three regions with boundaries at 0.55m, 1.15m, and 3m, with division counts of 5, 3, and 5, respectively. The grid model construction is shown in Figure 4.

Based on the known water-based mud invasion depth of 0.85m, by fitting the actual measured hydrocarbon phase fluid proportion and pressure with the computed results, the initial distribution of water phase saturation within the invasion range is set as shown in the Figure 5 below.

The numerical simulation results compared with historical measured results are shown in Figure 6 and Figure 7.

In this calculation process, the water phase compressibility is 0.000015 1/bar, and the hydrocarbon phase compressibility is 0.000144 1/bar. The probe operating depth is 3286.9 meters, and the total operating duration is 7880 seconds. The actual fluid sampling speed for the first 200 seconds averages 0.1 cc/s, and from 200 seconds onward, it is 4.38 cc/s. The monitored pressure difference at the end of the actual operation is 1154.82 psi, with the hydrocarbon breakthrough time observed at 3400 seconds. Compared with the calculation results, the average fitting rate for the hydrocarbon fluid proportion is 98.5%, the average fitting rate for the pressure curve is 98.67%, and the final simulated operation production pressure difference is 1133.05 psi, with a fitting rate of 98.11%.

Therefore, it is evident from the results that a satisfactory computational performance has been achieved. This analysis holds significant guidance for activities such as analyzing the factors affecting WFT pump extraction, predicting hydrocarbon content, selecting probes, and optimizing pump extraction operations in drilling processes.

3. Results and Discussion

3.1. Parameters Sensitivity Analysis

First, based on the reservoir and fluid property parameters of the aforementioned numerical simulation model, a sensitivity study was conducted on factors such as reservoir permeability and fluid sampling operation speed. With a water-based mud invasion depth of 0.4m, a water phase saturation of 100% within the invasion range, and a longitudinal grid step length of 0.02m, two types of simulation scenarios were set up:

(a)For elliptical probe operation at a sampling speed of 4 cc/s with permeability of 10 MD, 30 MD, and 50 MD.

(b)For the same permeability of 30 MD with sampling speeds of 3 cc/s, 4 cc/s, and 5 cc/s.

Key indicators such as hydrocarbon proportion and pressure curves were calculated and analyzed. The aim was to explore the response characteristics of key parameters like hydrocarbon content and pressure during the fluid sampling process under different reservoir permeability and sampling speed conditions. The calculation results are shown below.

In summary, it can be seen that reservoir permeability significantly affects bottom hole pressure during the sampling operation. As shown in Figure 8 and Figure 10, higher reservoir permeability increases the Poisson ratio of hydrocarbon flow, resulting in a higher hydrocarbon sampling rate per unit time and better fluid permeability. Reasonably increasing the sampling speed can help accelerate the flow of hydrocarbon fluid from the reservoir to the sampling probe.

Furthermore, as shown in Figure 9 and Figure 11, higher reservoir permeability requires a smaller production pressure differential during the sampling process. This indicates that less energy is needed to create effective bottom-hole pressure, thereby reducing the difficulty of fluid sampling and improving the breakthrough time and operational efficiency of wireline fluid sampling for hydrocarbon fluids. Conversely, if the reservoir permeability is relatively low, the hydrocarbon fluid permeability is poor. In this case, the WFT operation pump requires higher power to provide a greater production pressure differential, forming higher pressure energy at the bottom-hole sampling point, thus facilitating the outflow of hydrocarbon fluid.

Therefore, for WFT fluid sampling operations in high-permeability reservoirs, increasing the pump speed under reasonable operational pressure differentials and maximum pump efficiency can help quickly extract hydrocarbon fluids from the reservoir. In low-permeability reservoirs, reducing the pump speed within a reasonable pump efficiency range can help maintain a stable operational pressure differential.

Evaluating the fluid sampling process for different permeability reservoirs provides direct and quantifiable guidance for assessing the effectiveness of operational measures during WFT fluid sampling. This approach offers practical insights for optimizing sampling operations in varying reservoir conditions.

3.2. Probe Operation Analysis

Secondly, this study examines the operational effectiveness of probes with different inlet sizes during fluid sampling processes to provide technical support for probe selection. Specifically, under a reservoir permeability of 10 MD, with a fluid sampling operation time of 8000 seconds and a final pressure differential of 841.15 psi, we analyzed the fluid sampling performance of four different inlet-sized probes.

Under the same operational pressure differential, the pump operation speeds were as follows: Small Type Inlet Probe at 0.545 cc/s, Middle Type Inlet Probe at 1.12 cc/s, Ellipse Type Inlet Probe at 2.2 cc/s, and Large Type Inlet Probe at 3 cc/s. The effectiveness of different pump operation measures is illustrated in the following Figure 12 and Figure 13.

This indicates that in low-permeability reservoirs, different types of probes, due to variations in the size of their inlet areas and the influence of pump extraction speed, exhibit differences in the time it takes to first encounter hydrocarbons. Under the known minimum control requirements for production pressure drop based on pump efficiency, probes with larger inlet areas tend to encounter hydrocarbons earlier during fluid sampling. Additionally, they extract a greater volume of hydrocarbon phase fluid within the same duration, indicating higher operational efficiency. Particularly in environments with lower permeability reservoirs, under feasible pressure drop conditions, probes with larger inlet areas demonstrate a more pronounced advantage in operational efficiency. This aligns with the requirements of field production operations.

3.3. Sampling Operations in Shallow Heavy Oil Formations

For a certain offshore heavy oil block characterized by high viscosity, high density, and high colloidal content, field experiments have shown that a lack of reasonable external force during fluid sampling can lead to difficulties such as the inability of oil samples to flow out of the sampling bucket and excessive pressure drop. Additionally, the reservoir in this area is loosely cemented, with weak bonding between sand grains, making sand detachment prone under large pressure gradients. Therefore, in an effective fluid sampling process, it is necessary not only to ensure efficient fluid sampling but also to avoid instrument blockage due to sanding.

Taking Well A in this block as an example, the fluid viscosity of the sampled crude oil in the formation is 43 mPa.s, density 9.15 g/cm³, permeability 500 MD, porosity 22%, NTG 64%, probe sampling depth 1362.2 m, reference pressure 1904.34 psi, maximum controlled pressure differential for effective operation is 900 psi, maximum pump efficiency extraction rate is 9 cc/s, and sampling duration is 12000 s. For three operation schemes involving elliptical probe operation (sampling speed 3.72 cc/s), large-suction probe operation (sampling speed 5 cc/s), and 3D elliptical probe operation (sampling speed 8.85 cc/s), optimization analysis and evaluation of probe selection and sampling efficiency in the WFT process of shallow heavy oil reservoirs were conducted. After numerical simulation calculations for different schemes, the hydrocarbon phase proportion and pressure curve calculation results for each operation scheme are as follows: Figure 14

In this analysis, based on the sampling operation target depth of the water-based mud invasion depth of 0.3m into the reservoir, with a saturation of 100% within the invasion range, longitudinal grid division with an equal step length of 0.03m was conducted. The depth range of the formation was from 1361.8m to 1362.6m, with a wellbore radius of 0.15m, and a radial grid division of 31. The modeling radius in the plane was 3m. Three regions were defined in the plane circumferentially, with boundary distances of 0.45m, 1.45m, and 3m, and respective divisions of 2, 4, and 4. While ensuring effective control of the maximum operating pressure differential, pump efficiency, and sampling speed, simulated analyses of operational measures were conducted. Within the same implementation period of 12000s for the operation, the simulated final pressure differentials for the three schemes shown in Figure 15 were 870psi, meeting the condition of a maximum operating pressure differential of 900psi.

Furthermore, by setting the hydrocarbon phase fluid proportion greater than 0.4% as the breakthrough threshold and 80% as the compliance analysis threshold, a timeliness analysis was conducted for elliptical probe operation, large-suction probe operation, and 3D elliptical probe operation. The statistical results are presented in the Table 5 below.

Under the operating pressure difference range of maximum pump efficiency, the 3D elliptical probe operation has advanced breakthrough time by 2300s and 3100s compared to elliptical probe operation and large suction nozzle probe operation, and the compliance time has been advanced by 3900s and 7000s. The efficiency of fluid sampling with the 3D probe is significantly better than the other two single-suction nozzle probe operation modes. Additionally, combined with Figure 14, it can be observed that after reaching the breakthrough time, the amount of hydrocarbon-phase fluid sampled per unit time with the 3D probe operation is far superior to the other two probes, and it is more conducive to extracting more hydrocarbon-phase fluid under stable pressure difference conditions.

Furthermore, for the two single-suction nozzle probe operations, if the sampling speed is increased, it will face the risk of insufficient pump efficiency, significant pressure drop, large operating pressure difference, and instrument blockage due to sand production from the reservoir. Therefore, the operation of the 3D probe in the process of sampling shallow heavy oil reservoirs can reduce pressure differences, increase flow area, ensure more crude oil flows out of the formation, improve fluid sampling efficiency, reduce the risk of instrument blockage due to sand production, and optimize pump-assisted fluid sampling efficiency, thereby enhancing timeliness with satisfactory results.

4. Conclusions

The present study unveils an advanced numerical simulation technology grounded in the Finite Volume Method, designed to optimize fluid sampling operations within WFT during drilling endeavors. It elucidates various key techniques, spanning from meshing methodologies to numerical simulation modeling, matrix equation resolution, and the analysis of influential factors. By integrating these methodologies with empirical case data, a thorough analysis of operational efficacy is conducted. The primary conclusions drawn from this investigation are as follows:

Combining the variable step-size radial grid division technique considering the wellbore radius and the probe suction mouth size, a reasonable representation of the matching between the probe suction mouth and the grid model contact area can be achieved. By adopting a variable step-size grid model, not only can the permeation characteristics of the formation fluid into the sampling probe be reasonably represented locally, but it can also reduce the total number of grid models. This approach has a positive impact on the stability and timeliness of subsequent numerical simulation calculations.
Based on the principles of the finite volume method and considering the influence of the ratio of the probe suction area to the contact area of the geological formation grid, a numerical simulation model for fluid sampling in WTF was constructed. This model can incorporate conditions such as probe type and mud invasion. Additionally, by using incomplete LU decomposition matrix preprocessing and SBiCG, the rapid solution of model target parameters was achieved. The validation and comparison of simulated calculation results with actual sampling data were completed, providing important support for the optimization design of WFT operations during drilling processes.
Combining the numerical simulation model, sensitivity analysis of key information such as reservoir permeability, fluid sampling speed, and probe type was conducted. This further verified the response characteristics of calculation result curves such as hydrocarbon phase ratio and pressure in different operational measures during fluid sampling processes under various permeability, fluid sampling speed, and probe type environments. Hence, reliable bases were provided for the analysis of factors affecting the efficiency of WFT fluid sampling, in conjunction with reservoir properties and fluid characteristics.
Addressing the optimization design requirements for fluid sampling operations in shallow heavy oil sand-bearing formations, the challenges faced in environments characterized by high viscosity, high density, and high colloidal content of heavy oil reservoirs were analyzed. These challenges include difficulties such as oil sample failure to flow out of the sampling bucket, excessive pressure drop during operations, and instrument clogging due to sanding. A comparative analysis was conducted between 3D probes and traditional single-suction probes, demonstrating advantages such as increased flow area and enhanced hydrocarbon phase fluid sampling under low-pressure differential operational requirements. This approach can effectively prevent instrument clogging caused by high-pressure gradients.

Author Contributions

Conceptualization, L.W. and H.L.; Methodology, H.L. and R.H.; Software, J.W. and C.Z.; Validation, H.X. and X.L. (Xiaodong Li); Formal analysis, L.W. and X.L. (Xuan Li); Investigation, J.L; Writing—original draft, L.W. and R.H.; Writing—review & editing, H.L. and R.H. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

Symbols	Meanings	Units
$P$	Pressure	psi
$S_{w}$	Water phase saturation	-
$S_{o}$	Hydrocarbon phase saturation	-
$r_{e}$	Wellbore radius	m
$r_{w}$	The effective fluid supply radius of the wellbore	m
$T$	Transmissibility of the grid blocks	-
$A$	Contact area of the current grid block	m²
$L$	The distance from the center of the grid block	m
$N T G$	The effective thickness ratio	m/m
$K$	Permeability	MD
$d$	Probe area ratio coefficient	m²/ m²
$A_{p}$	Probe inlet area	m²
$A_{c}$	Contact area of the grid block	m²
$k_{r o}$	Relative permeability of the hydrocarbon phase	-
$k_{r w}$	Relative permeability of the water phase	-
$ρ_{o}$	Density of the hydrocarbon phase fluid	g/cm³
$ρ_{w}$	Density of the water phase fluid	g/cm³
$μ_{o}$	Viscosity of the hydrocarbon phase fluid	mPa.s
$μ_{w}$	Viscosity of the water phase fluid	mPa.s
$ϕ$	Porosity	-
$V$	Grid volume	m³
$c_{o}$	The compressibility of the hydrocarbon phase fluid	1/bar
$c_{w}$	The compressibility of the water phase fluid	1/bar
$c_{f}$	The comprehensive compressibility	1/bar
$Δ D$	The depth difference between grid blocks	m
$Δ t$	Iterative time step size	s
$g$	Gravitational acceleration	m/s²
$Q$	Total sampling fluid rate at the probe	cc/s
$f_{w}$	Water phase fluid content ratio	-
$f_{o}$	Hydrocarbon phase fluid content ratio	-
$q_{w}$	The flow rate of the water phase fluid	cc/s
$q_{o}$	The flow rate of the hydrocarbon phase fluid	cc/s
$P_{w f}$	Flowing pressure at the probe inlet	psi

References

Jogn, A.; Agarwal, A.; Gaur, M. Challenges and opportunities of wireline formation testing in tight reservoirs: a case study from Barmer Basin, India. Journal of Petroleum Exploration and Production Technology 2017, 7, 33–42. [Google Scholar]
Haibo, L.; Lejun, W.; Meng, W.; Hongwei, Z.; Xiaofei, Q. Calculation of depth of mud filtrate invasion based on formation sampling. Journal of Southwest Petroleum University (Science & Technology Edition) 2023, 45, 97–106. [Google Scholar]
Mingying, L. ; Yantuo Song.; Yinghui Jiang.; Libo Jiao. Application of wireline formation tester in moderate and low permeability depleted zones. Well Logging Technology 2023, 47, 764–771. [Google Scholar]
Xiangr, W.; Cancan, Z.; Changxue, W.; Wenli, C. Review on application of the wireline formation tester. Progress in Geophysics 2008, 23, 1579–1585. [Google Scholar]
Bo, Z.; Xiuwen, M.; Guo, T. The numerical simulation of wireline formation tester with finite element method. Journal of Jilin University：Earth Science Edition 2007, 37, 629–632. [Google Scholar]
Changxue, W.; Wenli, C.; Xiangr, W. Pressure gradient computation and application of the wireline formation tester. Petroleum Exploration and Development 2008, 35, 476–481. [Google Scholar]
Faruk, O.; Carlos, T.; Tarek, M. Estimation of in-situ petrophysical properties from wireline formation tester and induction logging measurements: A joint inversion approach. Journal of Petroleum Science and Engineering 2008, 63, 1–17. [Google Scholar]
Mathieu, N.; Gerard, M.; Herve, J. On the use of Wireline Formation Testing (WFT) data: 2. Consequences of permeability anisotropy and heterogeneity on the WFT responses inferred flow modeling. Journal of Petroleum Science and Engineering 2015, 133, 776–784. [Google Scholar]
Yuanliang, Z.; Shengquan, G.; Chuang, H.; Kang, B.; Shichen, S.; Zhaoya, F.; Tao, Z. Oil-gas sampling purity prediction based on numerical simulation for modular dynamic formation tester. Well Logging Technology 2020, 44, 534–538+552. [Google Scholar]
Quan, Z.; Lejun, W.; Xiaofei, Q. Applicable conditions of wireline formation testing instrument. Petroleum Geology and Engineerinng 2022, 36, 115–118. [Google Scholar]
Jinxiu, X.; Hongzhi, L.; Yunjiang, C. An application analysis of wireline formation test in Bohai region. China Offshore Oil and Gas 2008, 20, 106–110. [Google Scholar]
MA Junguan. Application of EFDT in sand producing formation sampling. China Petroleum and Chemical Standards and Quality 2019, 39, 53–54. [Google Scholar]
Chunzhao, Z. Sampling probe optimization technology of the double hanging instrument for Wireline Formation Testing. Offshore Oil 2022, 42, 59–64. [Google Scholar]
Lejun, W.; Haibo, L.; Meng, W. Research on EFDT working system optimization and applicability. Petrochemical Technology 2019, 26, 342–343. [Google Scholar]
Jude, T.; David, M.; Andrew, C. Surface and hypersurface meshing techniques for space–time finite element method. Computer-Aided Design 2023, 163, 1–28. [Google Scholar]
Busto, S.; Dumbser, M.; Río-Martín, L. An Arbitrary-Lagrangian-Eulerian hybrid finite volume/finite element method on moving unstructured meshes for the Navier-Stokes equations. Applied Mathematics and Computation 2023, 437, 1–36. [Google Scholar] [CrossRef]
Van, J.; Van, H. Space–time discontinuous Galerkin finite element method with dynamic grid motion for inviscid compressible flows. Journal of Computational Physics 2002, 182, 546–585. [Google Scholar]
Riemslagh, K.; Vierendeels, J.; Dick, E. Arbitrary Lagrangian-Eulerian finite-volume method for the simulation of rotary displacement pump flow. Applied Numerical Mathematics 2000, 32, 419–433. [Google Scholar] [CrossRef]
Vorst, H. Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems. SIAM Journal on Scientific and Statistical Computing 1992, 13, 631–644. [Google Scholar] [CrossRef]

Figure 1. Fixed step length grid division in plane and longitudinal directions.

Figure 2. Variable step length grid division in plane and longitudinal directions.

Figure 3. Schematic diagram of the calculation process of the numerical simulation.

Figure 4. Numerical simulation grid model division.

Figure 5. Initial distribution of water phase saturation due to mud invasion.

Figure 6. Comparison of historical and calculated hydrocarbon proportion.

Figure 7. Comparison of historical and calculated pressures.

Figure 8. Results of hydrocarbon proportion at the same speed with different permeabilities. Figure 9.

Figure 9. Results of pressure at the same speed with different permeabilities.

Figure 10. Results of hydrocarbon proportion at different speeds with the same permeability.

Figure 11. Results of pressure at different speeds with the same permeability.

Figure 12. Results of hydrocarbon proportion during operations with different probes at the same differential pressure.

Figure 13. Results of pressure during operations with different probes at the same differential pressure.

Figure 14. Results of hydrocarbon proportion during operations with different schemes.

Figure 15. Results of pressure during operations with different schemes.

Table 1. Parameters of Fluid Sampling Probes.

Type	Name	Inlet Area(in²)	Standard Ratio
S	Small Type Inlet Probe	0.21	0.27
M	Middle Type Inlet Probe	0.79	1
E	Ellipse Type Inlet Probe	2.32	2.94
L	Large Type Inlet Probe	6.23	7.89
3D-E	3D-Ellipse Type Inlet Probe	2.32*3	2.94*3

Table 2. Basic parameters of model.

Porosity (%)	Permeability (MD)	Hole diameter (m)	Initial pressure (psi)	Invasion depth （m）	Ratio of effective thickness (m/m)
15.7	25	0.1556	4520.95	0.85	0.912

Table 3. Liquid parameters of model.

Drilling fluid viscosity (mPa.s)	Formation fluid viscosity (mPa.s)	Drilling fluid density (g/cm³)	Formation fluid density (g/cm³)
0.55	2.50	1.01	0.903

Table 4. Fluid relative permeability experimental data.

Water saturation (%)	Hydrocarbon relative permeability (-)	Water relative permeability (-)
0	1	0
32	0.88	0
35	0.82	0.01
45	0.50	0.12
52	0.30	0.25
60	0.10	0.48
65	0.02	0.61
70	0.01	0.71
80	0	0.88

Table 5. Statistical results of the aging time.

	Sampling speed (cc/s)	Operating time Hydrocarbon >0.4% (s)	Operating time Hydrocarbon >80% (s)
E-Probe	3.72	7800	11900
L-Probe	5.00	6000	8800
3D-E-Probe	8.85	3700	4900

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

Research on Operation Optimization of Fluid Sampling in Wireline Formation Testing with Finite Volume Method

Abstract

1. Introduction

2. Materials and Methods

2.1. Introduction of Pumping Probe

2.2. Meshing Model

2.3. Numerical Simulation Model

2.4. Target Parameter Calculation Model

2.5. Solving Matrix Equation

2.6. Model Analysis

3. Results and Discussion

3.1. Parameters Sensitivity Analysis

3.2. Probe Operation Analysis

3.3. Sampling Operations in Shallow Heavy Oil Formations

4. Conclusions

Author Contributions

Data Availability Statement

Conflicts of Interest

Nomenclature

References

MDPI Initiatives

Important Links

Subscribe