Carbonate reservoirs are rich in oil and gas reserves; thus, they have great exploitation potential. Therefore research on the microscopic mechanisms of carbonate reservoirs is of great significance. Based on the thin section images of core castings of typical Well 555 and its pore and fracture features in actual reservoirs, this study designed three kinds of images representing the pore and throat structure of real rocks by applying image stitching and binarization processing methods. A microscopic pore model of carbonate rocks was then established using COMSOL numerical simulation software. The microscopic water flooding characteristics and residual oil distributions of different schemes were observed by designing different fracture development forms. The fractures that developed in parallel main lines showed a more obvious influence on water flooding characteristics compared to fractures that developed in vertical main lines. The cluster residual oil was the main residual oil type in the early stage of water flooding in the pure matrix model. With the progress of water flooding, the continuous cluster residual oil gradually turned into mainly discontinuous porous and columnar residual oils. Vertical mainline fractures reduced the amount of residual oil in clusters and replaced it with columnar residual oil. In contrast, parallel main line fractures expanded the unswept area, with the residual oil appearing in contiguous clusters. This study microscopically analyzed the law and characteristics of water flooding in carbonate reservoirs to provide key theoretical support for enhancing oil recovery.

The pore-throat radius of the shale oil reservoir is extremely small, and it is difficult to accurately obtain the absolute permeability and oil–water two-phase relative permeability of the actual oil reservoir through conventional core experiments. However, these parameters are very important for reservoir numerical simulation. In this paper, a method for characterizing flow parameters based on a pore network model that considers differential pressure flow and diffusion flow is proposed. Firstly, a digital core was reconstructed using focused ion beam scanning electron microscopy (FIB-SEM) from the Gulong shale reservoir in the Songliao Basin, China, and a pore network model was extracted. Secondly, quasi-static single-phase flow and two-phase flow equations considering diffusion were established in the pore network model. Finally, pore-throat parameters, absolute permeability, and oil–water two-phase permeability curves were calculated, respectively. The results show that the pore-throat distribution of the Gulong shale reservoir is mainly concentrated in the nanometer scale; the mean pore radius is 87 nm, the mean throat radius is 41 nm, and the mean coordination number is 3.97. The calculated permeability considering diffusion is 0.000124 mD, which is approximately twice the permeability calculated without considering diffusion. The irreducible water saturation of the Gulong shale reservoir is approximately 0.4, and the residual oil saturation is approximately 0.35. The method proposed in this paper can provide an important approach for characterizing the flow parameters of similar shale oil reservoirs.

In recent years, sand production has been frequently observed in offshore weakly consolidated sandstone reservoirs. Permeability changes due to sand migration seriously affect the confidence in well test interpretation, production forecasts, and oilfield development plan schedules. The purpose of this paper is to propose a comprehensive model of coupled sand migration, stress sensitivity, and high viscosity oil and to study the effect of sand production induced permeability zoning on transient pressure behavior by combining discrete boundary and discrete wellbore with the boundary element method. In this two-zone composite model, the reservoir can be divided into the inner zone with the improved permeability due to sand migration and the outer zone with initial reservoir permeability. The multifactor effects of stress-sensitive, highly viscous oil, sand migration, and horizontal well are included in this model. Thus, the seepage equation presents a highly nonlinear and difficult to obtain an accurate analytical solution. In this paper, the boundary element method (BEM) is introduced to separate the boundary and wellbore, and the semianalytical solution of the hybrid model is obtained. The comparative analysis of measured pressure curve fitting from a horizontal well, located in the eastern of the South China Sea, proves that this comprehensive model can be used for pressure transient analysis of the weakly consolidated sandstone reservoir. The flow regime analysis indicates that a two-zone composite system may develop seven flow regimes: the wellbore storage stage, early-time radial stage, first transition stage, inner linear stage, inner pseudoradial flow, transition flow from the inner area to the outer area, and outer pseudoradial flow. Sensitivity analysis indicates that the smaller the sand production radius, the shorter the duration of the transition flow from the inner to the outer zone, which suggests the well is mainly affected by the outer boundary in the later period. The larger the permeability ratio, the higher the pressure curves may move up.