A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes...A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.展开更多
Estimating the oil-water temperatures in flowlines is challenging especially in deepwater and ultra-deepwater offshore applications where issues of flow assurance and dramatic heat transfer are likely to occur due to ...Estimating the oil-water temperatures in flowlines is challenging especially in deepwater and ultra-deepwater offshore applications where issues of flow assurance and dramatic heat transfer are likely to occur due to the temperature difference between the fluids and the surroundings. Heat transfer analysis is very important for the prediction and prevention of deposits in oil and water flowlines, which could impede the flow and give rise to huge financial losses. Therefore, a 3D mathematical model of oil-water Newtonian flow under non-isothermal conditions is established to explore the complex mechanisms of the two-phase oil-water transportation and heat transfer in different flowline inclinations. In this work, a non-isothermal two-phase flow model is first modified and then implemented in the InterFoam solver by introducing the energy equation using OpenFOAM® code. The Low Reynolds Number (LRN) k-ε turbulence model is utilized to resolve the turbulence phenomena within the oil and water mixtures. The flow patterns and the local heat transfer coefficients (HTC) for two-phase oil-water flow at different flowlines inclinations (0°, +4°, +7°) are validated by the experimental literature results and the relative errors are also compared. Global sensitivity analysis is then conducted to determine the effect of the different parameters on the performance of the produced two-phase hydrocarbon systems for effective subsea fluid transportation. Thereafter, HTC and flow patterns for oil-water flows at downward inclinations of 4°, and 7° can be predicted by the models. The velocity distribution, pressure gradient, liquid holdup, and temperature variation at the flowline cross-sections are simulated and analyzed in detail. Consequently, the numerical model can be generally applied to compute the global properties of the fluid and other operating parameters that are beneficial in the management of two-phase oil-water transportation.展开更多
The pressure-sensitive effect on the pore structure of sandstone was investigated using X-ray computed micro-tomography and QEMSCAN quantitative mineral analysis. In a physical simulation study, we extracted the pore ...The pressure-sensitive effect on the pore structure of sandstone was investigated using X-ray computed micro-tomography and QEMSCAN quantitative mineral analysis. In a physical simulation study, we extracted the pore network model from digital cores at different confining pressures and evaluated the effect of pressure sensitivity on the multiphase displacement process. In both the pore network model and QEMSCAN scanning, the pore structure was observed to be damaged under a high confining pressure. Due to their different scales, the pores and throats exhibited inhomogeneous changes; further, the throats exhibited a significant variation compared to that exhibited by the pores. Meanwhile, the heterogeneity of the pore structure under the two aforementioned activities was aggravated by the elastic-plastic deformation of the pore structure.The pressure-sensitive effect increased the proportion of mineral particles, such as quartz(the main component of the core skeleton), and reduced the proportion of clay minerals. The clay minerals were originally attached to the pore walls or interspersed in the pores; however, as the pressure increased, the clay minerals accumulated in the pores resulting in blockage of the pores. While simulating the multiphase displacement process, increasing the confining pressure was observed to severely restrict the flowability of oil and water. This study promises to improve the efficiency of reservoir development in terms of oil and gas exploitation.展开更多
Food dye Brilliant Blue was introduced as the tracer in a dye-tracing experiment to obtain dye profile patterns of sandy loam soil, aeolian sandy soil, percolating paddy soil and permeable paddy soil. The dyed soil pr...Food dye Brilliant Blue was introduced as the tracer in a dye-tracing experiment to obtain dye profile patterns of sandy loam soil, aeolian sandy soil, percolating paddy soil and permeable paddy soil. The dyed soil profiles were then photographed and the photos were scanned into a computer. Edited with certain software, only the dyed areas were left on the profile photos, which indicted the preferential flow paths for water and solute transport. Fractal dimensions of the dye patterns were calculated according to Arnold's function. Soil particle size distribution was analyzed by pipette method. The regression analysis showed that there was significant relationship between soil clay content and fractal dimension D of the dye pattern of soil profile. Based on the experiment results, the possibility of introducing fractal dimension to estimation of soil sensitivity to preferential flow is discussed.展开更多
The heat transfer rate of the thermal Marangoni convective flow of a hybrid nanomaterial is optimized by using the response surface methodology(RSM).The thermal phenomenon is modeled in the presence of a variable incl...The heat transfer rate of the thermal Marangoni convective flow of a hybrid nanomaterial is optimized by using the response surface methodology(RSM).The thermal phenomenon is modeled in the presence of a variable inclined magnetic field,thermal radiation,and an exponential heat source.Experimentally estimated values of the thermal conductivity and viscosity of the hybrid nanomaterial are utilized in the calculation.The governing intricate nonlinear problem is treated numerically,and a parametric analysis is carried out by using graphical visualizations.A finite difference-based numerical scheme is utilized in conjunction with the 4-stage Lobatto IIIa formula to solve the nonlinear governing problem.The interactive effects of the pertinent parameters on the heat transfer rate are presented by plotting the response surfaces and the contours obtained from the RSM.The mono and hybrid nanomaterial flow fields are compared.The hybrid nanomaterial possesses enhanced thermal fields for nanoparticle volume fractions less than 2%.The irregular heat source and the thermal radiation enhance the temperature profiles.The high level of the thermal radiation and the low levels of the exponential heat source and the angle of inclination(of the magnetic field)lead to the optimized heat transfer rate(Nux=7.46275).展开更多
This letter reports traffic flow sensitivity to visco-elasticity, with the traffic flow modeling briefly described at first and then used to do traffic flow simulations whose results can reflect the properties of spat...This letter reports traffic flow sensitivity to visco-elasticity, with the traffic flow modeling briefly described at first and then used to do traffic flow simulations whose results can reflect the properties of spatial-temporal evolution of ring traffic flow. It reveals that visco-elasticity plays crucial role in formation of traffic flow patterns, implying that self-organization of traffic flow is crucial in determining traffic flow status.展开更多
Theoretical equations for computing sensitivity coefficients of wellbore pressures to estimate the reservoir parameters in low-permeability reservoirs conditioning to non-Darcy flow data at low velocity were obtained....Theoretical equations for computing sensitivity coefficients of wellbore pressures to estimate the reservoir parameters in low-permeability reservoirs conditioning to non-Darcy flow data at low velocity were obtained. It is shown by a lot of numerical calculations that the wellbore pressures are much more sensitive to permeability very near the well than to permeability a few gridblocks away from the well. When an initial pressure gradient existent sensitivity coefficients in the region are closer to the active well than to the observation well. Sensitivity coefficients of observation well at the line between the active well and the observation well are influenced greatly by the initial pressure gradient.展开更多
This paper presents a new analytical solution to investigate the mechanism of transient confinedunconfined flow in a confined aquifer induced by pumping with a large rate during mine drainage.The study focuses on unde...This paper presents a new analytical solution to investigate the mechanism of transient confinedunconfined flow in a confined aquifer induced by pumping with a large rate during mine drainage.The study focuses on understanding the impact of non-Darcian effect on flow towards a fully penetrated pumping well.The nonlinear relationship between specific discharge and the hydraulic gradient is described using Izbash's equation.A novel approximate method is developed to linearize the mathematical model,and the solution is derived using the Boltzmann transform.The proposed solution is validated by comparing it with previous works.The findings indicate that increased non-Darcian index,quasi-hydraulic conductivity,and specific storage have negatively affect the development of the unconfined region and aquifer drawdown,as greater turbulence flow accelerates recharge to the pumping well.Drawdown is found to be sensitive to the non-Darcian index,quasi-hydraulic conductivity,while it is unaffected by specific yield and specific storage.The conclusions provide valuable insights for mine drainage and the application of geological and hydrological conditions.展开更多
碎屑流是我国山区最危险的地质灾害之一,山区桥墩常受到碎屑流冲击而开裂、倾斜甚至倒塌,给山区桥梁建设、运营带来严重的安全隐患。采用离散元方法(discrete element method,DEM)和有限元方法(finite element method,FEM)耦合的三维数...碎屑流是我国山区最危险的地质灾害之一,山区桥墩常受到碎屑流冲击而开裂、倾斜甚至倒塌,给山区桥梁建设、运营带来严重的安全隐患。采用离散元方法(discrete element method,DEM)和有限元方法(finite element method,FEM)耦合的三维数值模拟方法模拟了碎屑流对双柱式桥墩的冲击效应,并结合斜槽试验,验证了耦合方法的准确性,进一步分析了碎屑流冲击坡度、距离和体积密度对桥墩冲击力的影响规律。结果表明,最大冲击力与碎屑流冲击坡度、距离和体积密度分别呈幂函数(指数大于1)、幂函数(指数小于1)和线性正相关。冲击坡度、距离和体积密度对最大冲击力的敏感度值分别为3.012、0.202、0.804,在桥梁碎屑流灾害防治时需重视冲击坡度和体积密度的影响。将冲击力的数值模拟值与流体动力学模型预测值对比分析表明,流体动力学模型理论公式能较好地预测桥墩所受的最大冲击力,最大预测误差低于23.6%。相关研究结果可为山区桥梁碎屑流灾害防治与设计提供一定的参考依据。展开更多
基金This study was supported by the National Natural Science Foundation of China(U22B2075,52274056,51974356).
文摘A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.
文摘Estimating the oil-water temperatures in flowlines is challenging especially in deepwater and ultra-deepwater offshore applications where issues of flow assurance and dramatic heat transfer are likely to occur due to the temperature difference between the fluids and the surroundings. Heat transfer analysis is very important for the prediction and prevention of deposits in oil and water flowlines, which could impede the flow and give rise to huge financial losses. Therefore, a 3D mathematical model of oil-water Newtonian flow under non-isothermal conditions is established to explore the complex mechanisms of the two-phase oil-water transportation and heat transfer in different flowline inclinations. In this work, a non-isothermal two-phase flow model is first modified and then implemented in the InterFoam solver by introducing the energy equation using OpenFOAM® code. The Low Reynolds Number (LRN) k-ε turbulence model is utilized to resolve the turbulence phenomena within the oil and water mixtures. The flow patterns and the local heat transfer coefficients (HTC) for two-phase oil-water flow at different flowlines inclinations (0°, +4°, +7°) are validated by the experimental literature results and the relative errors are also compared. Global sensitivity analysis is then conducted to determine the effect of the different parameters on the performance of the produced two-phase hydrocarbon systems for effective subsea fluid transportation. Thereafter, HTC and flow patterns for oil-water flows at downward inclinations of 4°, and 7° can be predicted by the models. The velocity distribution, pressure gradient, liquid holdup, and temperature variation at the flowline cross-sections are simulated and analyzed in detail. Consequently, the numerical model can be generally applied to compute the global properties of the fluid and other operating parameters that are beneficial in the management of two-phase oil-water transportation.
文摘The pressure-sensitive effect on the pore structure of sandstone was investigated using X-ray computed micro-tomography and QEMSCAN quantitative mineral analysis. In a physical simulation study, we extracted the pore network model from digital cores at different confining pressures and evaluated the effect of pressure sensitivity on the multiphase displacement process. In both the pore network model and QEMSCAN scanning, the pore structure was observed to be damaged under a high confining pressure. Due to their different scales, the pores and throats exhibited inhomogeneous changes; further, the throats exhibited a significant variation compared to that exhibited by the pores. Meanwhile, the heterogeneity of the pore structure under the two aforementioned activities was aggravated by the elastic-plastic deformation of the pore structure.The pressure-sensitive effect increased the proportion of mineral particles, such as quartz(the main component of the core skeleton), and reduced the proportion of clay minerals. The clay minerals were originally attached to the pore walls or interspersed in the pores; however, as the pressure increased, the clay minerals accumulated in the pores resulting in blockage of the pores. While simulating the multiphase displacement process, increasing the confining pressure was observed to severely restrict the flowability of oil and water. This study promises to improve the efficiency of reservoir development in terms of oil and gas exploitation.
基金Project supported by the National Key Basic Research Support FOundation(NKBRSF) of China(No.G19990ll708) and the Guangxi Uni,rsitv Science funds China(No.1701).
文摘Food dye Brilliant Blue was introduced as the tracer in a dye-tracing experiment to obtain dye profile patterns of sandy loam soil, aeolian sandy soil, percolating paddy soil and permeable paddy soil. The dyed soil profiles were then photographed and the photos were scanned into a computer. Edited with certain software, only the dyed areas were left on the profile photos, which indicted the preferential flow paths for water and solute transport. Fractal dimensions of the dye patterns were calculated according to Arnold's function. Soil particle size distribution was analyzed by pipette method. The regression analysis showed that there was significant relationship between soil clay content and fractal dimension D of the dye pattern of soil profile. Based on the experiment results, the possibility of introducing fractal dimension to estimation of soil sensitivity to preferential flow is discussed.
文摘The heat transfer rate of the thermal Marangoni convective flow of a hybrid nanomaterial is optimized by using the response surface methodology(RSM).The thermal phenomenon is modeled in the presence of a variable inclined magnetic field,thermal radiation,and an exponential heat source.Experimentally estimated values of the thermal conductivity and viscosity of the hybrid nanomaterial are utilized in the calculation.The governing intricate nonlinear problem is treated numerically,and a parametric analysis is carried out by using graphical visualizations.A finite difference-based numerical scheme is utilized in conjunction with the 4-stage Lobatto IIIa formula to solve the nonlinear governing problem.The interactive effects of the pertinent parameters on the heat transfer rate are presented by plotting the response surfaces and the contours obtained from the RSM.The mono and hybrid nanomaterial flow fields are compared.The hybrid nanomaterial possesses enhanced thermal fields for nanoparticle volume fractions less than 2%.The irregular heat source and the thermal radiation enhance the temperature profiles.The high level of the thermal radiation and the low levels of the exponential heat source and the angle of inclination(of the magnetic field)lead to the optimized heat transfer rate(Nux=7.46275).
基金support of Russian Foundation for Basic Research (RFBR 13-01-12056)National Natural Science Foundation of China (10972212)
文摘This letter reports traffic flow sensitivity to visco-elasticity, with the traffic flow modeling briefly described at first and then used to do traffic flow simulations whose results can reflect the properties of spatial-temporal evolution of ring traffic flow. It reveals that visco-elasticity plays crucial role in formation of traffic flow patterns, implying that self-organization of traffic flow is crucial in determining traffic flow status.
文摘Theoretical equations for computing sensitivity coefficients of wellbore pressures to estimate the reservoir parameters in low-permeability reservoirs conditioning to non-Darcy flow data at low velocity were obtained. It is shown by a lot of numerical calculations that the wellbore pressures are much more sensitive to permeability very near the well than to permeability a few gridblocks away from the well. When an initial pressure gradient existent sensitivity coefficients in the region are closer to the active well than to the observation well. Sensitivity coefficients of observation well at the line between the active well and the observation well are influenced greatly by the initial pressure gradient.
基金supported by the national natural science foundation of China(Grant Numbers 41807197,2017YFC0405900,and 51469002)the natural science foundation of Guangxi(Grant Numbers 2017GXNSFBA198087,2018GXNSFAA 138042,and GuiKeAB17195073)Hebei high level talent funding project(B2018003016).
文摘This paper presents a new analytical solution to investigate the mechanism of transient confinedunconfined flow in a confined aquifer induced by pumping with a large rate during mine drainage.The study focuses on understanding the impact of non-Darcian effect on flow towards a fully penetrated pumping well.The nonlinear relationship between specific discharge and the hydraulic gradient is described using Izbash's equation.A novel approximate method is developed to linearize the mathematical model,and the solution is derived using the Boltzmann transform.The proposed solution is validated by comparing it with previous works.The findings indicate that increased non-Darcian index,quasi-hydraulic conductivity,and specific storage have negatively affect the development of the unconfined region and aquifer drawdown,as greater turbulence flow accelerates recharge to the pumping well.Drawdown is found to be sensitive to the non-Darcian index,quasi-hydraulic conductivity,while it is unaffected by specific yield and specific storage.The conclusions provide valuable insights for mine drainage and the application of geological and hydrological conditions.
文摘碎屑流是我国山区最危险的地质灾害之一,山区桥墩常受到碎屑流冲击而开裂、倾斜甚至倒塌,给山区桥梁建设、运营带来严重的安全隐患。采用离散元方法(discrete element method,DEM)和有限元方法(finite element method,FEM)耦合的三维数值模拟方法模拟了碎屑流对双柱式桥墩的冲击效应,并结合斜槽试验,验证了耦合方法的准确性,进一步分析了碎屑流冲击坡度、距离和体积密度对桥墩冲击力的影响规律。结果表明,最大冲击力与碎屑流冲击坡度、距离和体积密度分别呈幂函数(指数大于1)、幂函数(指数小于1)和线性正相关。冲击坡度、距离和体积密度对最大冲击力的敏感度值分别为3.012、0.202、0.804,在桥梁碎屑流灾害防治时需重视冲击坡度和体积密度的影响。将冲击力的数值模拟值与流体动力学模型预测值对比分析表明,流体动力学模型理论公式能较好地预测桥墩所受的最大冲击力,最大预测误差低于23.6%。相关研究结果可为山区桥梁碎屑流灾害防治与设计提供一定的参考依据。