Multiscale and Multiphysical Approaches to Fluids Flow in Unconventional Reservoirs
View this Special IssueResearch Article  Open Access
Jianhui Zeng, Yongchao Zhang, Jingsheng Ma, Changling Liu, Zhe Cao, "Investigating the Influences of PoreScale Characteristics on Tight Oil Migration by a TwoPhase Pore Network Model", Geofluids, vol. 2020, Article ID 4890758, 13 pages, 2020. https://doi.org/10.1155/2020/4890758
Investigating the Influences of PoreScale Characteristics on Tight Oil Migration by a TwoPhase Pore Network Model
Abstract
The migration of expelled hydrocarbon from source rock into unconventional tight reservoirs is subject to different porescale fluid transport mechanisms as opposed to the conventional counterparts and therefore plays a crucial role in controlling the hydrocarbon distribution and accumulation in the former. One of the different mechanisms is related to the formation of a more viscous boundary layer (BL) of brine, i.e., wetting phase fluid on pore surfaces, giving rise to the socalled BL effect. In this work, a twophase pore network model (PNM) that considers this BL effect is developed to study the influences of porescale characteristics on the oil migration process, manifested through the BL effect in tightsandstone media. Good agreements are reached between experimentally derived relative permeability curves and predicted ones, by applying this model to the porenetwork networks extracted from the same samples. Then, this validated model was used to evaluate the impacts of the following factors on the oil migration process: pore radius, coordination number, aspect ratio, brine viscosity, and wettability. The results show that all factors can influence the oil migration process but at different magnitudes. The applicability and significance of the developed tight oil migration PNM are discussed in this work.
1. Introduction
Following the shale gas breakthrough in the United States, tight oil has become the next focus of unconventional petroleum exploration worldwide [1, 2]. Tight oil from the Bakken Formation has accounted for more than 10% of the United States daily oil production in 2017 [3]. Unlike conventional reservoirs, tight oil reservoirs are usually characterized by low porosity, low permeability, complicated porethroat structures, and strong heterogeneity ([4–6]). The hydrocarbon migration is a process in which hydrocarbon expelled from lowpermeability source rocks finds tortuous porous paths through carrier beds into the traps where hydrocarbon can accumulate to form reservoirs in the geological timeframe [7]. Of many factors that control this process and therefore the resultant reservoirs in terms of their sizes, fluid saturations, and fluid properties within them, which are of crucial importance for hydrocarbon exploration and exploitation [8, 9], the tightness of pore space within carrier and reservoir rocks is of pivotal importance. Compared with conventional reservoirs, oil migration behaves distinctively in tight reservoirs in the following aspects [5, 10]: (1) there is no clear demarcation between the primary migration and secondary migration process, (2) the migration distance is typically short which induces hydrocarbon inner source accumulation or nearsource accumulation, and (3) the impact of capillary pressure is significantly important during the migration process while the effect of buoyancy is limited.
1.1. Relationships between PoreScale Characteristics and Oil Migration
The significance of porescale characteristics for the oil migration process has been assessed using experimental and geological methods in previous literature. Lai et al. [11] analyzed the porestructure characteristics of tight sandstones in the Sichuan Basin, China, by using experimental data and basinscale inferential methods. They concluded that the complex pore systems in tight sandstones are associated with the original depositional environments and subsequent diagenesis alteration; the coexistence of different types of pores results in the polymodal distribution of pore sizes. Relationships between the porestructure characteristics and diagenesis alteration effects are concluded in their published review paper [12]. Similar work was reported by Xiong et al. [13]. They classified the pore system of tight sandstone into four types: intergranular pores, dissolution pores, clayaggregated pores, and microfractures, the sizes of which are distributed in the different ranges varying from nano to microscale. The flow potential during the oil migration process and the oil storage capacity in terms of oil migration and accumulation are contributed by the pores with different sizes. Cao et al. [14] further identified the key geological controlling factors of microscale oil distribution during oil migration and accumulation by geochemical approaches combined with mesoscale petrophysical methods, which includes porosity, pore throat radius, and hydrocarbon generation capacity. Zhang et al. [15] conducted and analyzed hydrocarbon charging experiments using natural tight sandstones samples. The experimental results show that the starting pressure of oil migration and terminal oil saturation values both show good logarithmic relationships with the sample permeability. However, the growth curves of oil saturation () with charging pressure are irregular and distinctive for different samples, such characteristics cannot be explained without sufficient porescale information of rock samples. To get a better understanding of the oil migration critical conditions and controlling factors at the pore scale, further investigations were completed in our previous works [16–18] by using a series of microscopic observation methods including CT (Xray computed tomography), SEM (scanning electron microscopy), NMR (nuclear magnetic resonance), and casting thin section techniques. Several concluding remarks obtained from these works are as follows: (1) there exists a critical pore radius threshold for oil migration in tight cores, which ranges from 0.035 to 0.627 μm; (2) the oilcharging curves of tight reservoirs could be classified into four types, where the porestructure factors play crucial roles; and (3) the surface areas of effective pores/throats, pore quantity density, and coordination number are proved to be important factors controlling the oil migration process.
Despite some qualitative relationships and critical conditions obtained from these existing experiments, the influencing mechanisms of the porescale characteristics on the microscale oil migration behaviors are not clearly understood. Besides, due to the complex interactions of different factors and high costs of the forced fluid displacement experiments, the quantitative relationship between any single porescale parameter and oil migration cannot be achieved experimentally.
1.2. Pore Network Modeling of Physical Processes at the Pore Scale
PNM (pore network modeling) is referred to any analytical and/or numerical models of physiochemical processes prescribed on pore networks composed of pore elements that are connected through throat elements. Since pore networks can be constructed to resemble geometric and topologic characteristics of pore space of realistic porous sample, PNM has been proved to be able to capture the manifestation of porescale processes, through simulation, on sample scales and predict macroscopic transport properties in porous media [19, 20]. Fatt [21] pioneered PNM and proposed a 2D regular network model for predicting capillary pressure and relative permeability of the primary drainage process in which the radii of pores were randomly assigned. Dullien et al. [22] extended that method to the 3D network model with more realistic porescale representation. The subsequent development of PNM has allowed the packing of grains with different sizes and shapes to be modeled, simulating a wide range of flow processes in different porous sedimentary rocks [23–26]. The readers are referred to the reviews [13, 27] for a summary of the recent advances in PNM models applied to various processes. Recently, PNM has been applied to tight formations. Riepe et al. [28] presented a case history of combining CT imaging and PNM to establish a relationship between microscale structures and petrophysical properties of tight clastic rocks in Oman. That application shows the promise of an alternative approach for the evaluation of unconventional reservoirs. Ruspini et al. [29] introduced a multiscale pore network modeling workflow to compute the transport properties of clayrich tight sandstone samples with wide pore size distributions. The representative network is composed of macropores, defined as completely resolved porosity, and micropores which are defined as unresolved porosity and treated as the equivalent continuum. Wang and Sheng [30, 31] proposed a quasistatic PNM model to simulate the drainage process in shale and tight porous formations. The liquid nonDarcy flow mechanism was incorporated into the proposed PNM model using an empirical equation obtained from experiments [32]. Based on their simulations, Wang and Sheng concluded that the nonDarcy flow behaviors show significant influences on the absolute permeability predictions, but not on the relative permeability results.
PNM models have been developed for modeling the oil migration process in conventional formations, but not suitable for directly studying the oil migration process in tight formations. Compared with medium or highpermeability media, the fluid flow in tight media shows obvious nonDarcy flow (NDF) behaviors [33, 34], which is the flow velocity of fluid being not linearly dependent on the charging pressure gradient. Some researchers [35–37] proposed that these NDF behaviors were attributed to the existence of the boundary layer (BL) of water attached to the pore surface where the viscosity is much larger than the fluid in the pore center. The BL with different thicknesses in pores significantly reduces the effective seepage radius, which results in the NDF behaviors in tight media. In recent works, a parametric model of BL thickness of BL (see Equation (1)) was proposed by Cheng’s group [14, 38, 39] based on the experiments where coefficients are fitted using the particle swarm optimization algorithm [40]. The calculation equation is given as where is the thickness of BL, is the brine viscosity, and is the pressure difference between inlet and outlet per unit length. The three coefficients (, , and ) are determined to be 0.258, 0.261, and 0.419. The proposed method shows good predictions of BL thickness of numerous microtube experiment results conducted by Li et al. (2011) as shown in Figure 1. These experiments are conducted through simulating the porethroat system with microfused silica capillary tubes, the radius of which ranges from 2.5 μm to 10 μm. The influence of BL thickness on the fluid NDF is analyzed by measuring the flow rate, pressure gradient, and viscosity experimentally.
A PNM for simulating the oil migration behaviors is developed in this work, considering NDF behavior due to BL of water. Equation (1) is employed to calculate BL thicknesses in the network elements to account for the NDF behaviors during tight oil migration. We are aware that the developed model employed the BL theory may be neither universally applicable nor accurate for all the cases of tight formations. The principal objectives of this work are as follows: (1) to extend the application of PNM to analyze the oil migration in tight formations; (2) to study the influences of porescale characteristics on oil migration using PNM method, which includes the radius of pores, average coordination number, aspect ratio (the average diameter ratio of pores to throats), brine viscosity, and wettability conditions; and (3) to discuss the significances and limitations of the conducted PNM simulations.
2. A Quasistatic TwoPhase PNM with BL for Modeling the Oil Migration in Tight Formation
In this section, a quasistatic twophase porenetwork model with BL is developed for modeling oil migration in tight formation, by adapting a wellestablished quasistatic PNM model for the primary drainage process. This newly developed model is implemented to support a simulation workflow to be discussed in a later section.
2.1. The Model Development
The oil migration into initially waterfilled waterwet pore networks occurs in forms of pistonlike displacement [41]. When the driving force of oil phase exceeds the capillary entry pressure () of the neighboring waterfilled network elements, oil will displace water in pores/throats. We add a BL of immobile water of finite thickness, , to each network element whose radius and volume available to oil phase are reduced proportionally. So, as a result, we should expect higher capillary entry pressure at any oil pressure and higher oil relative permeability and lower water relative permeability at any water saturation.
The capillary entry pressure, for a cylindrical throat, is calculated as follows: where is the interfacial tension between the oil and water, is the contact angle that quantifies the wettability property, and is the pore radius. For an angular pore element, the capillary entry pressure can be adapted from that of Øren [42] as follows:
The parameter in Equations (3) and (4) is the shape factor for the network element, which can be calculated through Equation (6): where is the crosssectional area and is the perimeter of the element.
To calculate permeability of each connected phase from the inlet to outlet faces of a PNM model, the singlephase flow between two neighboring pore elements is given by a Poiseuilletype law [43]: where and are the pressure values of pore and pore and is the conductance between the two pores, which can be determined through the following Equation (8).
Here, , , and are the conductances of pores and and the throat connecting the pores.
In this work, the conductance of network elements is calculated following the functional form in the “3Rs” approach [44, 45] but takes into consideration BL, and it is expressed as where is the conductance constants, is the conductance exponent, is the length of the element. in Equation (9) is the average viscosity of the fluids in the pores and calculated by where and are the fractions of oil and water in the calculation element. Based on the calculations of Equation (7), the pressure field of the network model can be calculated by applying the mass conservation at each pore element:
When the pressure field for a phase in the network is obtained, we can calculate the relative permeability of each phase (oil or water phase) as follows:
The new model is implemented on an earlier version of the opensource PNM code (https://github.com/ahboujelben/numSCAL_basic) by extending its quasistatic twophase model. For technical details in that code, the reader is referred to Regaieg et al. [46] and Boujelben et al. [47].
2.2. Simulation Procedure
Simulations are carried according to the procedure illustrated by Figure 2. For a given pore network that may be extracted from Xray CT images of real rock samples or stochastically generated, (1) it is initialized 100% water saturation () at a reference pressure while wettability condition is applied to each network element; (2) explore the minimum and maximum radius (Min and Max ) of elements, then calculate the corresponding maximum and minimum entry capillary pressure (Max Pc and Min Pc) with considering the BL effect. The calculations of Max Pc and Min Pc were realized by the iterative method; (3) set an incremental charging pressure () which increases from Min Pc to Max Pc. The incremental step is set by (Max PcMin Pc)/simulation steps; (4) explore for and determine the elements to invade. An invasion event occurs when the charging pressure overcomes the entry pressure of element. The entry pressure of each element is calculated by Equations (3)–(5). (5) Update phase saturation and flow rate with the consideration of BL effect. (6) Check the terminal condition whether the charging pressure is larger than Max Pc. If no, increase the charging pressure by step and do loop for steps (4) to (6). If yes, output the phase saturation and calculate relative permeability for each phase.
3. Model Validations
Model validation has been conducted by making comparisons between the predicted relative permeability and experimental results for a Berea sandstone sample and two tight sandstone samples (Samples S5 and S6). The relative permeability data measured on the Berea sample and its matched CT image data can be downloaded from the Imperial College database (the website address is http://www.imperial.ac.uk/earthscience/research/researchgroups/perm/research/porescalemodelling/). Samples S5 and S6 are taken from the Chang7 members of the Yangchang Formation in the Ordos Basin, a typical tight oil basin in northwest China. Oil charging experiments are carried out by using Samples S5 and S6 to simulate the physical process of oil migration in tight media. Experimental apparatus, experimental procedures, and the properties of used water/oil have been described in our previous work [15, 48]. Both tight samples are scanned to obtain CT images using Zeiss Xradia500 Versa with a resolution of 0.99 μm, at the Stata Key Laboratory of Petroleum Resources and Respecting, China. Raw CT images are preprocessed and reconstructed to 16bit grey images using a multithreading software package “tomoRecon” [49]. The image filtering and phase segmentation are completed by using Avizo 9.7 software. An interactive thresholding module is used to segment the phases by selecting values of image intensity ranges for each phase [50]. The properties of the three experimental samples are shown in Table 1. It shows that the porosity, permeability, and pore connectivity of the tight sandstone samples (S5 and S6) are much smaller than those of the Berea sample. The pore size ranges of the tight samples are distributed in smaller and narrower scope compared to the Berea sample. Note that the lower limit of pore radius for tight samples is limited by the solution of CT images which means some pores with radius smaller than 1 μm may not be detected by CT scanning. More information concerning the porestructure features in the tight formations of the Ordos Basin can be found in our previous work [17, 18].

A modified maximal ball (MB) method developed by Dong and Blunt [51] was applied to extract pore networks from 3D CT images. The modified MB method explore and build an inscribed sphere at each void voxel of images first by a twostep searching algorithm, then remove those included in others. The rest spheres (MBs) are defined as pores and throats followed the clustering rule of family trees sorted by their size and rank. Geometrical properties are calculated for each pore and throat, including the size of pore and throat, volume, and shape factor. For more details of the pore network extraction, the readers are referred to Dong and Blunt’s original publication [51]. In our work, pore networks of tight sandstone samples were extracted using an opensource software package (website for download: https://github.com/aliraeini/pnextract) that implements those mentioned techniques. The flow modeling of these three samples is conducted followed the steps in Sections 2.1 and 2.2. For the three samples, the predicted relative permeability curves from the newly developed model with consideration of BL effect are compared with those without BL effect, while both the model predictions are compared to experimental data (Figures 3 and 4). Figure 3 indicates that both the PNM models with and without BL consideration show good agreements with the experimental data, which in turn partly proves that Equation (1) for BL calculation is basically reasonable for the pore size distribution in natural samples. Figure 4 shows that the predictions of the BL model are in better agreement with the experimental data than that without BL effect in tight sandstone samples. Compared to the model without BL effect, the BL model shows lower permeability predictions for each phase of which locates in a narrower and higher range. The differences between NBL and BL models are mainly attributed to the decrease of pore radius and connectivity in BL model caused by the existence of BL. Specific impacts of the single influencing factor on relative permeability predictions are further discussed in Sections 4 and 5.
4. Sensitivity Analyses of Five PoreScale Factors
Stochastically synthetic pore networks were constructed for the sensitivity analysis, with respect to five featured porescale factors (i.e., the radius of pores, average coordination number, aspect ratio, brine viscosity, and wettability), on their impacts on the oil migration. The changing scopes of these factors are determined by considering the porestructure characteristics of tight formations and referring to the measured results in several published works [4, 10, 53]. The resultant relative permeability curves are employed for showing the relative fluid flow capability during the oil migration, while oil charging curves for explaining the oil content accumulated. Technically, the water in the BL is simplified as dead water without flowing but count for in the calculation of relative permeability and . Note that each pore network is generated as a regular lattice network first and throats were then removed randomly to match the prescribed average coordination number. The size of pore elements follows a uniform distribution in the prescribed size range. Pore networks and their parameters set for sensitivity analysis are given in Table 2.

Runs 1~3 are used to assess the impact of radius on relative permeability and oil charge where the coordination number remains the same but the upper limits of pore size ranges increase from 1 to 5 μm with a constant lower limit of 0.1 μm. The calculated absolute permeabilities of the models are 0.001 mD, 0.03 mD, and 0.9 mD. Figure 5 shows that with the decrease of upper limit pore radius range the maximums of the oil/water phase relative permeability decrease. Meanwhile, the endpoints of at the axis increase from 3.5% to 20.8%. The predicted with small pores at the low (less than 35%) is higher than that with large pores but turns to be lower at a high condition. The predicted shows an opposite trend within the corresponding saturation range. The at the intersection points of and curves stayed constant for Runs 1~3. The oil charging curves in Figure 6 show that the calculated of models increases gradually with increasing charging pressure and terminates at the of 96.5%, 86.5%, and 79.2% for Runs 1~3, respectively. The required charging pressure for 0.1 ~ 1 μm model is much larger than that with 0.1 ~ 5 μm.
4.1. The Impact of Average Coordination Number
The coordination number () is defined as the average number of throats connected with each pore body, which quantifies the connectivity of the networks [54, 55]. A higher coordination number means a better connection condition of the pores in the network. Runs 4~6 set with different coordination numbers (i.e., , , and ) are used to evaluate the impact of coordination number on oil migration behaviors. The calculated absolute permeabilities of the models are 0.01 mD, 0.03 mD, and 0.06 mD. Figure 7 shows that the maximums of oil/water phase relative permeability decrease with the decreasing values. Meanwhile, the endpoints of curves at the axis are kept constant (%) for different simulations. The calculated in the poor connectivity model (Run 4) is slightly higher than that in the better connectivity models (Runs 5 and 6) at a low condition but turns to be lower at a high condition. The predicted curve shows an opposite trend. The corresponding axis values at the intersection points of and curves stayed constant while the corresponding axis values decrease with the decrease of values. The oil charging curves in Figure 8 show that the calculated of all three models increase with the charging pressure and terminates the same %. However, the required charging pressure reaching the maximum decreases with the increasing of the networks.
4.2. The Impact of Aspect Ratio
The impact of aspect ratio on oil migration is realized by increasing the pore radius while keeping the throat radius constant. Runs 7~9 with different aspect ratios are built and calculated in this section of whom the calculated absolute permeabilities are all 0.03 mD for the throat radii of models remain unchanged. Figure 9 shows that the influence of aspect ratio on relative permeability is relatively weak; the predicted of Run 9 is slightly lower than others (Runs 7~8) within a limited water saturation range of 0 ~ 35%. The oil charging curves in Figure 10 show that the maximums of in Runs 7~9 are 77.8%, 86.5%, and 89.2%, respectively. The required charging pressure reaching the maximum for three models are equal.
4.3. The Impact of Brine Viscosity
Runs 10~12 with different viscosity () settings are used to assess the impact of brine viscosity on oil migration behaviors. The calculated relative permeability results in Figure 11 show that the predicted and values both decrease with the increasing brine viscosity; under the condition of high water saturation, the impact of brine viscosity on relative permeability is more obvious which induces that the intersections point of and shift towards the bottomright area of the coordinates. The oil charging curves in Figure 12 show that the of all three models increase with the charging pressure and terminates the maximum of 63.2%, 74.6%, and 86.5%. The required charging pressure reaching the maximum for three models are equal.
4.4. The Impact of Wettability Condition
Sensitivity simulations of the wettability effects in conventional reservoir media have been performed in previous works [56, 57]. The wettability condition before the hydrocarbon migration process is assumed to be water wet. In our simulations, contact angle (CA) is used for qualifying the wettability of rock surface [58, 59]. Three models with different wettability settings (°, Run 13; °, Run 14; and °, Run 15) are modeled to assess the influence of wettability on oil migration, where CA is assigned to be uniformly distributed in the pores and throats. Figure 13 shows that the calculated becomes larger when CAs increase from 0° to 60°, while the wettability setting has little influence on the predictions. Figure 14 shows that of the three networks grow and terminate at the same maximum value of 86.7% with the increase of charging pressure. The required charging pressure reaching the maximum decreases with the increasing CA under the limited waterwet condition.
5. Discussions
5.1. The Applicability and Limitations of the Developed PNM Model
Numerous models have been developed in the past decades to simulate the fluid flow processes in porous media. The forces set in these models mainly involve capillary pressure, viscous force, and gravity force. These network models simulating the twophase flow system can be divided into two types: quasistatic model and dynamic model. The quasistatic model is suitable to simulate the fluid flow when the capillary pressure is dominating while the viscous force is negligible. When the nonnegligible viscous force is present; however, a dynamic model needs to be considered. When the capillary number is less than 10^{4}, the simulations of the quasistatic model are believed to be accurate enough [19]. In conventional reservoir formations, oil migration can also be divided into two types: steady oil migration and episodic oil migration [60]. The oil flow under the steady oil migration mode is assumed to be governed by capillary pressure, and the flow in the episodic migration is governed by viscous pressure while the natural buoyancy plays a role in the fluid flow. However, for the lowpermeability reservoir (e.g., tight reservoir and shale reservoir), whether both the steady oil migration and episodic oil migration can be identified to explain the flow phenomena in lowpermeability media is still a question. A common consensus is that the effect of buoyancy on fluid flow can be neglected because of the extremely small size of pores in tight formations [61, 62]. Due to the same reason, the effect of capillary pressure becomes significantly important, which means the fluid flow in tight media is more likely to follow the steady migration mode. This is the reason why the quasistatic method is employed in our present work. However, some other scholars [63] also proposed that present day reservoirs may not represent the conditions of the reservoir during the time of hydrocarbon expulsion and migration, indicating that the viscous force may have been the main controlling factor, and buoyancy may have been an important force forming the hydrocarbon distribution in tight oil/gas migrates and accumulates.
Some other limitations of the PNM methods should be noticed here. (1) The fluid inside the networks is considered incompressible and immiscible. (2) The real porethroat system of the rock is simplified as simple 3D geometric shapes following the porenetwork generation rules. (3) The chemical interactions between the fluids and solid matrix are neglected. (4) The wettability conditions are set as constantly distributed, which may not exactly match the real condition in nature [64]. (5) The gravity forces and influence of microfracture in the porethroat system are not considered in the developed model. For these limitations, the PNM model may not provide precise predictions for the actual geological conditions.
The permeability of the tested tight samples may be slightly larger than the permeability scope of tight sandstone by generally speaking. But the upper limit of tight media permeability is not standard. In some literature, tight oil is defined as oil resource reserved in very low permeability formation which cannot be produced at economic rates unless the formation is simulated by large hydraulic fracture treatments or produced by the use of horizontal wellbores [65]. The main reason we select these two samples with larger permeability is that the forced displacement experiments are hard to implement when the permeability of the sample is smaller than 0.5 mD. For the same reason, a complete series of experiments using different samples cannot be completed. Furthermore, the BL theory is employed in this work to explain the nonDarcy flow phenomenon in tight media. This theory is merely one of the numerous nonDarcy theories. Some other theory such as threshold pressure gradient theory can be also used for explaining the nonDarcy flow [66]. It may be acceptable to use the BL theory under most conditions, although the range of application, prediction accuracy, and limiting conditions of this theory remains to be further investigated.
5.2. The Simulation Results and Their Significances
Based on the reservoir exploration and production experiences, previous literature has concluded that the migration of tight oil behaves differently from the conventional reservoirs (as described in the Introduction). Some of these behaviors could be explained using the obtained simulation results in this work. For example, the initial water saturation (, the measured water saturation before the well drilling) in tight reservoirs is reported to be much higher than that in the conventional reservoirs [63]. Seen from the simulation results in Sections 4.1 and 4.2, the increases dramatically with the decrease of pore size and pore connectivity. The calculated with the pore size of 0.1 ~ 1 μm (absolute permeability mD) is 20.8% and the calculated with the pore size of 0.1 ~ 5 μm (absolute permeability mD) is only 3.5%. Note that the low here does not mean most of water could be replaced by the injected oil during the oil migration process in the natural system, which is only an idealized calculation value under the ideal hypothesis conditions.
The influential mechanisms of several factors on relative permeability and oil saturation during the oil migration are analyzed as follows. The decrease of pore radius leads to the decreases of both oil and water relative permeability and the increase of BL thickness. This increase of BL thickness further results in a higher irreducible water saturation and a lower effective seepage radius. Thus, the twophase seepage zone in the relative permeability curves decreases with the decreasing pore size; the maximum of oil saturation which corresponds to the irreducible water saturation in the oil saturation curves also shows a decreasing trend with the decrease of pore radius. Under the condition of the same pore radius distribution, the scope of twophase seepage remains unchanged with the change of pore connectivity. However, the decrease of coordination number will reduce the seepage area of the media which leads to the decrease of relative permeability for oil and water. Besides, under the condition of poor connectivity, water phase is easy to stay in the corner of the networks. This phenomenon may cause relative permeability to move towards the high water saturation direction. The changes of connectivity of the pore systems may change the pore space volume and water/oil content by absolute values but show no influence on the irreducible water saturation. Therefore, the maximum of oil saturation would keep constant at a high injection pressure. The differences between the relative permeability curves of Berea sample and tight samples (Figures 3 and 4) are mainly attributed to these aforementioned two factors. The impact of brine viscosity on oil migration is relatively simple. The thickness of the BL would increase with the increase of brine viscosity which causes the decrease of effective seepage radius for oil and water phases as well as the increase of irreducible water saturation. Meanwhile, a higher brine viscosity would directly reduce the permeability of water phase, which causes the decreasing trend of water relative permeability in Figure 11 to be much more obvious than that of oil with the increasing brine viscosity. Overall, the simulated five factors all show influence on the relative permeability and , but with different sensitivities. Compared with aspect ratio and hydrophilia of the rock, the pore radius, pore connectivity, and brine viscosity play more important roles in controlling the oil migration behaviors.
The improved oil migration PNM provides a new solution for analyzing the effects of complex porescale characteristics on tight oil migration and accumulation. Combining the quantitative results obtained from models and traditional geological analysis methods may be more helpful for studying the mechanisms of tight oil reservoir formation and evolution.
6. Conclusions
To better understand and quantitatively analyze the relationship between the oil migration characteristics and pore structures in tight formations, a twophase PNM with consideration of the BL effect is developed in this work. Several conclusions are drawn as follows: (1)This work proves that PNM method with consideration to the BL effect is applicable and effective for tight oil migration analyses, which shows better prediction accuracy than that without the BL effect(2)Five important factors including pore radius, coordination number, aspect ratio, brine viscosity, and wettability show influences on the tight oil migration, but sensitivities of these factors on oil migration are different(3)Based on the simulations, the influences of pore radius, coordination number, and brine viscosity on oil migration are more pronounced while the influences of aspect ratio and hydrophilia of the rock surface are relatively weak(4)With the decrease of pore radius, coordination number or the increase of brine viscosity in the oilphase relative permeability is weakened, which induces less oil accumulation in the pore system
Data Availability
The experimental and simulation data used to support the findings of this study have been deposited in the Mendeley repository (the download address is https://data.mendeley.com/datasets/j3467gsr3k/1). The data is also available on request.
Conflicts of Interest
The author(s) declare(s) that they have no conflicts of interest.
Acknowledgments
This work was financially supported by the National Natural Science Foundation of China (grants 41972147 and 41330319). Y. Zhang would show his sincere thanks to the HeriotWatt University, UK, for hosting his academic visit during which the work reported here was carried out in part. Ma thanks UK NERC (grant NE/R018022/1) for financial support.
References
 J. Boak and R. Kleinberg, “Shale gas, tight oil, shale oil and hydraulic fracturing,” Future Energy, vol. 6795, 2020. View at: Publisher Site  Google Scholar
 J. D. Hughes, “a reality check on the shale revolution,” Nature, vol. 494, no. 7437, pp. 307308, 2013. View at: Publisher Site  Google Scholar
 BP Energy, BP Energy Outlook, Energy, BP, 2018.
 J. Cai and X. Hu, Petrophysical Characterization and Fluids Transport in Unconventional Reservoirs, Elsevier, 2019.
 C. R. Clarkson and K. Pedersen, “Production analysis of Western Canadian unconventional light oil plays,” in Proceedings Canadian Unconventional Resources Conference, Calgary, Alberta, Canada, November, 2011, Society of Petroleum Engineers. View at: Publisher Site  Google Scholar
 P. H. Nelson, “Porethroat sizes in sandstones, tight sandstones, and shales,” AAPG Bulletin, vol. 93, no. 3, pp. 329–340, 2009. View at: Publisher Site  Google Scholar
 T. T. Schowalter, “Mechanics of secondary hydrocarbon migration and entrapment,” AAPG Bulletin, vol. 63, no. 5, pp. 723–760, 1979. View at: Google Scholar
 X. Luo, “Simulation and characterization of pathway heterogeneity of secondary hydrocarbon migration,” AAPG Bulletin, vol. 95, no. 6, pp. 881–898, 2011. View at: Publisher Site  Google Scholar
 P. Ungerer, J. Burrus, B. P. Y. C. Doligez, P. Y. Chenet, and F. Bessis, “Basin evaluation by integrated twodimensional modeling of heat transfer, fluid flow, hydrocarbon generation, and migration,” AAPG Bulletin, vol. 74, no. 3, pp. 309–335, 1990. View at: Google Scholar
 C. Zou, Unconventional petroleum geology, Elsevier, 2017. View at: Publisher Site
 J. Lai, G. Wang, Z. Fan et al., “Insight into the pore structure of tight sandstones using NMR and HPMI measurements,” Energy & Fuels, vol. 30, no. 12, pp. 10200–10214, 2016. View at: Publisher Site  Google Scholar
 J. Lai, G. Wang, Z. Wang et al., “A review on pore structure characterization in tight sandstones,” Earth Science Reviews, vol. 177, pp. 436–457, 2018. View at: Publisher Site  Google Scholar
 Q. Xiong, T. G. Baychev, and A. P. Jivkov, “Review of pore network modelling of porous media: experimental characterisations, network constructions and applications to reactive transport,” Journal of Contaminant Hydrology, vol. 192, pp. 101–117, 2016. View at: Publisher Site  Google Scholar
 Z. Cao, G. Liu, H. Zhan, Y. Kong, Z. Niu, and D. Zhao, “Geological control factors of micro oil distribution in tight reservoirs,” Marine and Petroleum Geology, vol. 77, pp. 1193–1205, 2016. View at: Publisher Site  Google Scholar
 Y. Zhang, J. Zeng, Z. Dai et al., “Experimental investigation on oil migration and accumulation in tight sandstones,” Journal of Petroleum Science Engineering, vol. 160, pp. 267–275, 2018. View at: Publisher Site  Google Scholar
 X. Feng, J. Zeng, S. Tao, J. Yang, S. Feng, and Z. Pang, “Oilcharging porethroat radius threshold of tight reservoirs: a comparison on multimethod calculation results,” Journal of Nanoscience and Nanotechnology, vol. 17, no. 9, pp. 6067–6076, 2017. View at: Publisher Site  Google Scholar
 J. Qiao, J. Zeng, X. Feng, Z. Yang, Y. Zhang, and S. Feng, “Characteristics of nanomicro pore networks and petroleum microscopic occurrence state in ultralow permeability (tight) sandstone reservoir,” Journal of Nanoscience and Nanotechnology, vol. 17, no. 9, pp. 6039–6050, 2017. View at: Publisher Site  Google Scholar
 J. Zeng, X. Feng, S. Feng, Y. Zhang, J. Qiao, and Z. Yang, “Influence of tight sandstone micronano porethroat structures on petroleum accumulation: evidence from experimental simulation combining Xray tomography,” Journal of Nanoscience and Nanotechnology, vol. 17, no. 9, pp. 6459–6469, 2017. View at: Publisher Site  Google Scholar
 M. J. Blunt, Multiphase Flow in Permeable Media: A PoreScale Perspective, Cambridge University Press, 2017. View at: Publisher Site
 J. Ma, X. Zhang, Z. Jiang, H. Ostadi, K. Jiang, and R. Chen, “Flow properties of an intact MPL from nanotomography and pore network modeling,” Fuel, vol. 136, pp. 307–315, 2014. View at: Publisher Site  Google Scholar
 I. Fatt, “The network model of porous media: capillary pressure characteristics,” Transactions of the AIME, vol. 207, no. 1, pp. 144–181, 2013. View at: Google Scholar
 F. A. L. Dullien, I. Chatzis, and M. S. El Sayed, “Modeling transport phenomena in porous media by networks consisting of nonuniform capillaries,” in Proceedings SPE Annual Fall Technical Conference and Exhibition, New Orleans, Louisiana, October, 1976, Society of Petroleum Engineers. View at: Publisher Site  Google Scholar
 M. J. Blunt, B. Bijeljic, H. Dong et al., “Porescale imaging and modelling,” Advances in Water Resources, vol. 51, pp. 197–216, 2013. View at: Publisher Site  Google Scholar
 J. Ma, J. P. Sanchez, K. Wu, G. D. Couples, and Z. Jiang, “A pore network model for simulating nonideal gas flow in micro and nanoporous materials,” Fuel, vol. 116, pp. 498–508, 2014. View at: Publisher Site  Google Scholar
 H. Sun, L. Duan, L. Liu et al., “The influence of microfractures on the flow in tight oil reservoirs based on porenetwork models,” Energies, vol. 12, no. 21, p. 4104, 2019. View at: Publisher Site  Google Scholar
 D. Wang, C. Wang, C. Li et al., “Effect of gas hydrate formation and decomposition on flow properties of finegrained quartz sand sediments using Xray CT based pore network model simulation,” Fuel, vol. 226, pp. 516–526, 2018. View at: Publisher Site  Google Scholar
 V. JoekarNiasar and S. M. Hassanizadeh, “Analysis of fundamentals of twophase flow in porous media using dynamic porenetwork models: a review,” Critical Reviews in Environmental Science and Technology, vol. 42, no. 18, pp. 1895–1976, 2012. View at: Publisher Site  Google Scholar
 L. Riepe, M. H. Suhaimi, M. Kumar, and M. A. Knackstedt, “Application of high resolution microCTimaging and pore network modeling (PNM) for the petrophysical characterization of tight gas reservoirsa case history from a deep clastic tight gas reservoir in Oman,” in Proceedings SPE Middle East Unconventional Gas Conference and Exhibition, Muscat, Oman, 2011, Society of Petroleum Engineers. View at: Publisher Site  Google Scholar
 L. Ruspini, G. Lindkvist, S. Bakke, L. Alberts, A. Carnerup, and P. E. Øren, “A multiscale imaging and modeling workflow for tight rocks,” in Proceedings SPE Low Perm Symposium, Denver, Colorado, USA, 2016, Society of Petroleum Engineers. View at: Publisher Site  Google Scholar
 X. Wang and J. J. Sheng, “Pore network modeling of the nonDarcy flows in shale and tight formations,” SPE Journal, vol. 163, pp. 511–518, 2018. View at: Publisher Site  Google Scholar
 X. Wang and J. J. Sheng, “Spontaneous imbibition analysis in shale reservoirs based on pore network modeling,” Journal of Petroleum Science and Engineering, vol. 169, pp. 663–672, 2018. View at: Publisher Site  Google Scholar
 X. Wang and J. J. Sheng, “Effect of lowvelocity nonDarcy flow on well production performance in shale and tight oil reservoirs,” Fuel, vol. 190, pp. 41–46, 2017. View at: Publisher Site  Google Scholar
 J. Cai, E. Perfect, C. L. Cheng, and X. Hu, “Generalized modeling of spontaneous imbibition based on HagenPoiseuille flow in tortuous capillaries with variably shaped apertures,” Langmuir, vol. 30, no. 18, pp. 5142–5151, 2014. View at: Publisher Site  Google Scholar
 Z. Zeng and R. Grigg, “A criterion for nonDarcy flow in porous media,” Transport in Porous Media, vol. 63, no. 1, pp. 57–69, 2006. View at: Publisher Site  Google Scholar
 Y. Huang, Z. Yang, Y. He, and X. Wang, “An overview on nonlinear porous flow in low permeability porous Media,” Theoretical and Applied Mechanics Letters, vol. 3, no. 2, article 022001, 2013. View at: Publisher Site  Google Scholar
 B. Zeng, L. Cheng, and C. Li, “Low velocity nonlinear flow in ultralow permeability reservoir,” Journal of Petroleum Science and Engineering, vol. 80, no. 1, pp. 1–6, 2011. View at: Publisher Site  Google Scholar
 Y. Yang, K. Wang, Q. Lv et al., “Flow simulation considering adsorption boundary layer based on digital rock and finite element method,” Petroleum Science, 2020. View at: Publisher Site  Google Scholar
 M. Chen, L. Cheng, R. Cao, C. Lv, J. Wu, and H. Liu, “A pore network model for studying boundary layer effect on fluid flow in tight formation (Russian),” in Proceedings SPE Annual Caspian Technical Conference and Exhibition, Baku, Azerbaijan, November, 2017, Society of Petroleum Engineers. View at: Publisher Site  Google Scholar
 X. Tian, L. Cheng, Y. Yan, H. Liu, W. Zhao, and Q. Guo, “An improved solution to estimate relative permeability in tight oil reservoirs,” Journal of Petroleum Exploration and Production Technology, vol. 5, no. 3, pp. 305–314, 2015. View at: Publisher Site  Google Scholar
 A. E. Bakyani, H. Sahebi, M. M. Ghiasi et al., “Prediction of CO_{2}–oil molecular diffusion using adaptive neurofuzzy inference system and particle swarm optimization technique,” Fuel, vol. 181, pp. 178–187, 2016. View at: Publisher Site  Google Scholar
 P. H. Valvatne and M. J. Blunt, “Predictive porescale modeling of twophase flow in mixed wet media,” Water Resources Research, vol. 40, no. 7, 2004. View at: Publisher Site  Google Scholar
 P. E. Oren, S. Bakke, and O. J. Arntzen, “Extending predictive capabilities to network models,” in Proceedings SPE Annual Technical Conference and Exhibition, San Antonio, Texas, October 1997, Society of Petroleum Engineers. View at: Publisher Site  Google Scholar
 J.C. Cai, “A fractal approach to low velocity nonDarcy flow in a low permeability porous medium,” Chinese Physics B, vol. 23, no. 4, article 044701, 2014. View at: Publisher Site  Google Scholar
 I. Bondino, S. R. McDougall, and G. Hamon, “Porescale Modelling of the effect of viscous pressure gradients during heavy oil depletion experiments,” Journal of Canadian Petroleum Technology, vol. 50, no. 2, pp. 45–55, 2013. View at: Publisher Site  Google Scholar
 J. Cruichshank, S. R. Mcdougall, and K. S. Sorbie, “Anchoring methodologies for porescale network models: application to relative permeability and capillary pressure prediction,” Petrophysics, vol. 43, p. 04, 2002. View at: Google Scholar
 M. Regaieg, S. R. McDougall, I. Bondino, and G. Hamon, “Finger thickening during extraheavy oil waterflooding: simulation and interpretation using porescale modeling,” PLoS One, vol. 12, no. 1, article e0169727, 2017. View at: Publisher Site  Google Scholar
 A. Boujelben, S. McDougall, M. Watson, I. Bondino, and N. Agenet, “Pore network modelling of low salinity water injection under unsteadystate flow conditions,” Journal of Petroleum Science and Engineering, vol. 165, pp. 462–476, 2018. View at: Publisher Site  Google Scholar
 J. Zeng, Y. Zhang, S. Zhang, J. Qiao, X. Feng, and S. Feng, “Experimental and theoretical characterization of the natural gas migration and accumulation mechanism in lowpermeability (tight) sandstone cores,” Journal of Natural Gas Science Engineering, vol. 33, pp. 1308–1315, 2016. View at: Publisher Site  Google Scholar
 M. L. Rivers, “tomoRecon: highspeed tomography reconstruction on workstations using multithreading,” in Developments in XRay Tomography VIII, vol. 8506, International Society for Optics and Photonics, 2012, 85060U. View at: Google Scholar
 C. Lore, User's GuideAvizo 9.7, 2018, https://assets.thermofisher.com/TFSAssets/MSD/ProductGuides/userguideavizosoftware.pdf.
 H. Dong and M. J. Blunt, “Porenetwork extraction from microcomputerizedtomography images,” Physical Review, vol. 80, no. 2, article 036307, 2009. View at: Google Scholar
 M. J. Oak, “Threephase relative permeability of waterwet Berea,” in Proceedings SPE/DOE Enhanced Oil Recovery Symposium, Tulsa, Oklahoma, April 1990, Society of Petroleum Engineers. View at: Publisher Site  Google Scholar
 H. Zhao, Z. Ning, Q. Wang et al., “Petrophysical characterization of tight oil reservoirs using pressurecontrolled porosimetry combined with ratecontrolled porosimetry,” Fuel, vol. 154, pp. 233–242, 2015. View at: Publisher Site  Google Scholar
 H. J. Vogel and K. Roth, “Quantitative morphology and network representation of soil pore structure,” Advances in Water Resources, vol. 24, no. 34, pp. 233–242, 2001. View at: Publisher Site  Google Scholar
 Y. Yang, Y. Li, J. Yao et al., “Formation damage evaluation of a sandstone reservoir via porescale Xray computed tomography analysis,” Journal of Petroleum Science and Engineering, vol. 183, p. 106356, 2019. View at: Publisher Site  Google Scholar
 M. J. Blunt, “Pore level modeling of the effects of wettability,” SPE Journal, vol. 2, no. 4, pp. 494–510, 1997. View at: Publisher Site  Google Scholar
 S. Zou, R. T. Armstrong, J. Y. Arns, C. H. Arns, and F. Hussain, “Experimental and theoretical evidence for increased ganglion dynamics during fractional flow in mixedwet porous media,” Water Resources Research, vol. 54, no. 5, pp. 3277–3289, 2018. View at: Publisher Site  Google Scholar
 J. S. Buckley, “Effective wettability of minerals exposed to crude oil,” Current Opinion in Colloid & Interface Science, vol. 6, no. 3, pp. 191–196, 2001. View at: Publisher Site  Google Scholar
 Y. Zhang, J. Zeng, J. Qiao, X. Feng, and Y. Dong, “Investigating the effect of the temperature and pressure on wettability in crude oil–brine–rock systems,” Energy & Fuels, vol. 32, no. 9, pp. 9010–9019, 2018. View at: Publisher Site  Google Scholar
 J. Zeng and Z. Jin, “Experimental investigation of episodic oil migration along fault systems,” Journal of Geochemical Exploration, vol. 7879, pp. 493–498, 2003. View at: Publisher Site  Google Scholar
 L. Chen, X. Zhu, L. Wang, H. Yang, D. Wang, and M. Fu, “Experimental study of effective amphiphilic graphene oxide flooding for an ultralowpermeability reservoir,” Energy & Fuels, vol. 32, no. 11, pp. 11269–11278, 2018. View at: Publisher Site  Google Scholar
 J. Ma, K. Wu, Z. Jiang, and G. D. Couples, “SHIFT: An implementation for lattice Boltzmann simulation in lowporosity porous media,” Physical Review E, vol. 81, no. 5, article 056702, 2010. View at: Publisher Site  Google Scholar
 W. K. Camp, “Porethroat sizes in sandstones, tight sandstones, and shales: discussion,” AAPG Bulletin, vol. 95, no. 8, pp. 1443–1447, 2011. View at: Publisher Site  Google Scholar
 Y. Liu, S. Iglauer, J. Cai, M. A. Amooie, and C. Qin, “Local instabilities during capillarydominated immiscible displacement in porous media,” Capillarity, vol. 2, no. 1, pp. 1–7, 2019. View at: Publisher Site  Google Scholar
 S. Holditch, “Tight gas sands,” Journal of Petroleum Technology, vol. 58, no. 6, pp. 86–93, 2006. View at: Publisher Site  Google Scholar
 G. Han, Y. Liu, K. Nawnit, and Y. Zhou, “Discussion on seepage governing equations for low permeability reservoirs with a threshold pressure gradient,” Advances in GeoEnergy Research, vol. 2, no. 3, pp. 245–259, 2018. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Jianhui Zeng et al. 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.