Cyclic supercritical multi-thermal fluid stimulation (CSMTFS) is a novel technology that can efficiently recover heavy oil, while the heating effect, production and heat loss characteristics of CSMTFS have not been discussed. In this study, a physical simulation experiment of CSMTFS is conducted with a three-dimensional experimental system. The results of the study indicate that the whole process of CSMTFS can be divided into four stages, namely, the preheating stage, production increase stage, production stable stage and production decline stage, of which the production stable stage is the main oil production stage, and the production decline stage is the secondary oil production stage. In the first two stages of the CSMTFS process, there is no supercritical multi-thermal fluid chamber, and only a relatively small supercritical multi-thermal fluid chamber is formed in the last stage of the CSMTFS process. Out of the supercritical multi-thermal fluid chamber, supercritical water in the thermal fluids condensates to hot water and flows downward to heat the subjacent oil layer. At the same time, the non-condensate gas in the thermal fluids accumulates to the upper part of the oil layer and reduces heat loss. The analysis of heat loss shows that the heat loss rate gradually increases at first and then tends to be stable. Compared with conventional thermal fluid, the CSMTFS can more effectively reduce heat loss. The enthalpy value of supercritical multi-thermal fluid is significantly increased compared with that of multi-thermal fluid, which effectively solves the problem of insufficient heat carrying capacity of multi-thermal fluid. Overall, cyclic supercritical multi-thermal fluid stimulation can effectively solve the problems of conventional heavy oil thermal recovery technology in offshore heavy oil recovery; it is indeed a new improved-oil-recovery technique for offshore heavy oil. The findings of this study can help in better understanding the cyclic supercritical multi-thermal fluid stimulation process. This study is significant and helpful for application of CSMTFS technology in heavy oil recovery.

Complex fracture systems including natural fractures and hydraulic fractures exist in shale gas reservoir with fractured horizontal well development. The flow of shale gas is a multi-scale flow process from microscopic nanometer pores to macroscopic large fractures. Due to the complexity of seepage mechanism and fracture parameters, it is difficult to realize fine numerical simulation for fractured horizontal wells in shale gas reservoirs. Mechanisms of adsorption–desorption on the surface of shale pores, slippage and Knudsen diffusion in the nanometer pores, Darcy and non-Darcy seepage in the matrix block and fractures are considered comprehensively in this paper. Through fine description of the complex fracture systems after horizontal well fracturing in shale gas reservoir, the problems of conventional corner point grids which are inflexible, directional, difficult to geometrically discretize arbitrarily oriented fractures are overcome. Discrete fracture network model based on unstructured perpendicular bisection grids is built in the numerical simulation. The results indicate that the discrete fracture network model can accurately describe fracture parameters including length, azimuth and density, and that the influences of fracture parameters on development effect of fractured horizontal well can be finely simulated. Cumulative production rate of shale gas is positively related to fracture half-length, fracture segments and fracture conductivity. When total fracture length is constant, fracturing effect is better if single fracture half-length or penetration ratio is relatively large and fracturing segments are moderate. Research results provide theoretical support for optimal design of fractured horizontal well in shale gas reservoir.

Air–foam flooding is an important mean to reduce water cut, improve oil production rate, enhance oil displacement efficiency and oil recovery in old oilfield. Because of high cost of core displacement experiment, numerical simulation is an important method to get optimal injection-production parameters of air–foam flooding, which can improve displacement effect and direct field test in Wuliwan District 1 of Jing’an oilfield. Based on laboratory experiment of air–foam flooding, parameters including model component, low temperature oxidation reaction, physical/chemical phenomena, foam interpolation function were set, and then numerical simulation concept model of air–foam flooding was established, which realized accurate fitting for 1-D core flooding experiment. According to numerical simulation concept model of air–foam flooding, the injection-production parameters including injection volume, gas–liquid ratio, and injection time are optimized. Results indicated that optimum surfactant concentration was 0.5%, optimum injection volume of air–foam was 0.25PV, optimum gas–liquid ratio was 1.5:1, optimum injection time was 65% water cut when air–foam flooding began. The researches were applied in Wellbock ZJ53, Wuliwan District 1 of Jing’an Oilfield. There were 15 well groups with air–foam injection and 63 producers, with annual oil production rate scale of 50000 tons. The oil increase effect was very obvious and the application in Jing’an Oilfield was successful.

The use of polymer gel for water control and oil addition is a common technical method in oilfield development. The polymer and hydrated bentonite react under the action of an organic crosslinking agent to form a composite gel. The particle-size change and microstructure of the composite gel were analyzed via shear thinning, thixotropic, viscoelastic, and start-up stress rheology experiments. The experimental results show that the polyacrylamide/bentonite organic crosslinked composite gel was a gel system with bentonite as the core aggregate structure, and the large particle-size distribution was mostly increased with increasing crosslinker content. The composite gel presented shear thinning characteristics, the content of bentonite or crosslinking agent was increased, and the shear resistance was stronger at a high shear rate. The composite gel exhibited positive thixotropic properties, and the thixotropy increased with increasing bentonite content. The composite gel had good viscoelastic characteristics, the elastic characteristics of the composite gel showed more significantly with bentonite increases, and the viscosity of the composite gel showed its characteristics more significantly with the crosslinking agent increased. After loading at a rate on the composite gel, the shear stress increased significantly with time and reached its maximum value, and then the shear stress decreased and gradually stabilized.

The adsorption of polymer microspheres in a stratum can directly affect its action mode and performance in the actual application process. Understanding the adsorption pattern of polymer microspheres and their adsorption mechanism can facilitate optimization of the application mode and enhance the use efficiency. Ultraviolet spectrophotometry was employed to measure the static adsorption characteristics of polymer microspheres (PMS) on the surface of quartz sand. The PMS adsorption capacity on the surface of quartz sand increased with increasing concentration. When the concentration was 1000 mg/L, the static equilibrium adsorption capacity was 402 μg/g, and monolayer adsorption was dominant. The effect of the contact time on the adsorption was investigated, and the fitting was performed using the isothermal adsorption thermodynamic equilibrium model and the adsorption kinetic model. The adsorption of 800 mg/L PMS tended to equilibrate after 0.8 h of adsorption on the surface of quartz sand, and the adsorption of 1400 mg/L PMS tended to equilibrate after 1 h of adsorption on the surface of quartz sand. Good fitting results of the kinetic adsorption process were obtained using the pseudo-first-order (PFO) model, pseudo-second-order (PSO) model, Elovich model, and mixed-order (MO) model. The effects of the temperature, particle size of the quartz sand, solid–liquid ratio, and salinity on the adsorption of PMS on the surface of quartz sand were examined. The PMS adsorption capacity on the surface of quartz sand decreased with increasing environmental temperature. The adsorption of PMS at the solid–liquid interface was an exothermic process, and the enthalpy of adsorption was negative. As the mass of the quartz sand in the solid–liquid ratio increased, the adsorption capacity decreased; a low salinity and neutral pH were conducive to the adsorption of PMS on the surface of quartz sand.

To evaluate the production performance of shale gas wells accurately, a quadruple-porosity medium model based on kerogen, the inorganic matrix, and natural-fracture–hydraulic-fracture network is proposed to simulate shale gas production. This model includes an apparent-kerogen permeability model to characterize the combined effects of gas slippage, Knudsen diffusion, surface diffusion, stress sensitivity, and matrix shrinkage. The fractal theory is used to characterize the non-equidistant distribution of secondary fractures in the hydraulic fracture network quantitatively. In addition, the effect of the stress sensitivity of the hydraulic fracture network on the production is considered, and a numerical method for the quadruple-porosity model is presented. After the validation of the model with field production data, the sensitivity parameters were analyzed to study the effects of the apparent kerogen permeability, Langmuir volume, fractal dimension, primary and secondary fracture conductivities, geomechanical effect, stimulated reservoir volume (SRV), and inorganic-matrix permeability on the gas production. This study provides a theoretical model for a more accurate evaluation of the shale gas production and improves the shale gas seepage theory.