This paper has developed the first meshless numerical solver for the reservoir black-oil model. The meshless solver uses a point cloud that is virtually unrestricted by topological constraints during generation and the extended finite volume method to discretize the reservoir computational domain and the governing equations, respectively. It can directly use the mature nonlinear solvers in grid-based simulators to calculate the global nonlinear discrete equations and obtain the distribution of oil, gas, and water saturations and pressures. This paper has implemented four examples with regular and irregular computational domains to test the computational performance of the meshless simulator. The results of the examples demonstrate that, compared to grid-based methods, the meshless simulator can effectively manage heterogeneous reservoir models and significantly mitigate the challenges associated with discretizing complex reservoir computational domains. Moreover, it avoids the grid orientation effect of the finite volume method that utilizes Cartesian grids by reducing the computational error by more than six times and exhibits good convergence.

1.
Benito
,
J. J.
,
Urena
,
F.
, and
Gavete
,
L.
, “
Influence of several factors in the generalized finite difference method
,”
Appl. Math. Modell.
25
(
12
),
1039
1053
(
2001
).
2.
Benito
,
J. J.
,
Urena
,
F.
,
Gavete
,
L.
, and
Alvarez
,
R.
, “
An h-adaptive method in the generalized finite differences
,”
Comput. Methods Appl. Mech. Eng.
192
(
5–6
),
735
759
(
2003
).
3.
Chávez-Negrete
,
C.
,
Domínguez-Mota
,
F. J.
, and
Santana-Quinteros
,
D.
, “
Numerical solution of Richards' equation of water flow by generalized finite differences
,”
Comput. Geotech.
101
,
168
175
(
2018
).
4.
Coats
,
K. H.
,
Coats
,
K. H.
,
Thomas
,
L. K.
, and
Pierson
,
R. G.
, “
Compositional and black oil reservoir simulation
,” in
SPE Reservoir Simulation Conference?
(
SPE
,
1995
), Paper No.
SPE-29111
.
5.
Fu
,
Z. J.
,
Tang
,
Z. C.
,
Zhao
,
H. T.
et al, “
Numerical solutions of the coupled unsteady nonlinear convection-diffusion equations based on generalized finite difference method
,”
Eur. Phys. J. Plus
134
(
6
),
272
(
2019
).
6.
Gavete
,
L.
,
Gavete
,
M. L.
, and
Benito
,
J. J.
, “
Improvements of generalized finite difference method and comparison with other meshless method
,”
Appl. Math. Modell.
27
(
10
),
831
847
(
2003
).
7.
Gu
,
Y.
,
Wang
,
L.
,
Chen
,
W
et al, “
Application of the meshless generalized finite difference method to inverse heat source problems
,”
Int. J. Heat Mass Transfer
108
,
721
729
(
2017
).
8.
Heinemann
,
Z. E.
,
Brand
,
C.
,
Munka
,
M.
, and
Chen
,
Y. M.
, “
Modeling reservoir geometry with irregular grids
,” in
SPE Reservoir Simulation Conference?
(
SPE
,
1989
), Paper No. SPE-18412.
9.
Liszka
,
T.
, “
An interpolation method for an irregular net of nodes
,”
Int. J. Numer. Methods Eng.
20
(
9
),
1599
1612
(
1984
).
10.
Li
,
P. W.
,
Fan
,
C. M.
, and
Grabski
,
J. K.
, “
A meshless generalized finite difference method for solving shallow water equations with the flux limiter technique
,”
Eng. Anal. Boundary Elem.
131
,
159
173
(
2021
).
11.
Milewski
,
S.
, “
Meshless finite difference method with higher order approximation—Applications in mechanics
,”
Arch. Comput. Methods Eng.
19
(
1
),
1
49
(
2012
).
12.
Palagi
,
C. L.
and
Aziz
,
K.
, “
Use of Voronoi grid in reservoir simulation
,”
SPE Adv. Technol. Ser.
2
(
02
),
69
77
(
1994
).
13.
Qu
,
W.
and
He
,
H.
, “
A spatial–temporal GFDM with an additional condition for transient heat conduction analysis of FGMs
,”
Appl. Math. Lett.
110
,
106579
(
2020
).
14.
Rao
,
X.
, “
An upwind generalized finite difference method (GFDM) for meshless analysis of heat and mass transfer in porous media
,”
Comp. Part. Mech.
10
,
533
554
(
2023
).
15.
Rao
,
X.
,
He
,
X.
,
Kwak
,
H.
et al, “
A novel streamline simulation method for fractured reservoirs with full-tensor permeability
,”
Phys. Fluids
36
(
1
),
013107
(
2024a
).
16.
Rao
,
X.
,
Liu
,
Y.
, and
Zhao
,
H.
, “
An upwind generalized finite difference method for meshless solution of two-phase porous flow equations
,”
Eng. Anal. Boundary Elem.
137
,
105
118
(
2022a
).
17.
Rao
,
X.
,
Zhao
,
H.
, and
Liu
,
Y.
, “
A novel meshless method based on the virtual construction of node control domains for porous flow problems
,”
Eng. Comput.
40
(
1
),
171
211
(
2024b
).
18.
Rao
,
X.
,
Zhao
,
H.
, and
Liu
,
Y.
, “
A meshless numerical modeling method for fractured reservoirs based on extended finite volume method
,”
SPE J.
27
,
3525
3540
(
2022b
).
19.
Saucedo-Zendejo
,
F. R.
,
Reséndiz-Flores
,
E. O.
, and
Kuhnert
,
J.
, “
Three-dimensional flow prediction in mould filling processes using a GFDM
,”
Comp. Part. Mech.
6
(
3
),
411
425
(
2019
).
20.
Suchde
,
P.
,
Kuhnert
,
J.
, and
Tiwari
,
S.
, “
On meshfree GFDM solvers for the incompressible Navier–Stokes equations
,”
Comput. Fluids
165
,
1
12
(
2018
).
21.
Suchde
,
P.
, “
A meshfree Lagrangian method for flow on manifolds
,”
Int. J. Numer. Methods Fluids
93
(
6
),
1871
1894
(
2021
).
22.
Zheng
,
S.
,
Wang
,
X.
,
Chen
,
J.
et al, “
Multiphase flow simulation of fractured karst oil reservoirs applying three-dimensional network models
,”
Phys. Fluids
36
(
7
),
076610
(
2024
).
You do not currently have access to this content.