In the petroleum industry,detection of multi-phase fluid flow is very important in both surface and down-hole measurements.Accurate measurement of high rate of water or gas multi-phase flow has always been an academic...In the petroleum industry,detection of multi-phase fluid flow is very important in both surface and down-hole measurements.Accurate measurement of high rate of water or gas multi-phase flow has always been an academic and industrial focus.NMR is an efficient and accurate technique for the detection of fluids;it is widely used in the determination of fluid compositions and properties.This paper is aimed to quantitatively detect multi-phase flow in oil and gas wells and pipelines and to propose an innovative method for online nuclear magnetic resonance(NMR)detection.The online NMR data acquisition,processing and interpretation methods are proposed to fill the blank of traditional methods.A full-bore straight tube design without pressure drop,a Halbach magnet structure design with zero magnetic leakage outside the probe,a separate antenna structure design without flowing effects on NMR measurement and automatic control technology will achieve unattended operation.Through the innovation of this work,the application of NMR for the real-time and quantitative detection of multi-phase flow in oil and gas wells and pipelines can be implemented.展开更多
Most current lattice Boltzmann (LBM) models suffer from the deficiency that their parameters have to be obtained by fitting experimental results. In this paper, we propose a new method that integrates the molecular ...Most current lattice Boltzmann (LBM) models suffer from the deficiency that their parameters have to be obtained by fitting experimental results. In this paper, we propose a new method that integrates the molecular dynamics (MD) simulation and LBM to avoid such defect. The basic idea is to first construct a molecular model based on the actual components of the rock-fluid system, then to compute the interaction force between the rock and the fluid of different densities through the MD simulation. This calculated rock-fluid interaction force, combined with the fluid-fluid force determined from the equation of state, is then used in LBM modeling. Without parameter fitting, this study presents a new systematic approach for pore-scale modeling of multi-phase flow. We have validated this ap- proach by simulating a two-phase separation process and gas-liquid-solid three-phase contact angle. Based on an actual X-ray CT image of a reservoir core, we applied our workflow to calculate the absolute permeability of the core, vapor-liquid H20 relative permeability, and capillary pressure curves.展开更多
Hydraulic slotting in a gas drainage borehole is an effective method of enhancing gas drainage performance. However, it frequently occurs that a large amount of slotting products(mainly the coal slurry and gas) intens...Hydraulic slotting in a gas drainage borehole is an effective method of enhancing gas drainage performance. However, it frequently occurs that a large amount of slotting products(mainly the coal slurry and gas) intensely spurt out of the borehole during the slotting, which adversely affects the slotting efficiency. Despite extensive previous investigations on the mechanism and prevention-device design of the spurt during ordinary borehole drilling, a very few studies has focused on the spurt in the slotting process. The slotting spurt is mainly caused by two reasons: the coal and gas outburst in the borehole and the borehole deslagging blockage. This paper focuses on the second reason, and investigates the hydraulic deslagging flow patterns in the annular space between the drill pipe and borehole wall.Results show that there are six deslagging flow patterns when the drill pipe is still: pure slurry flow, pure gas flow, bubble flow, intermittent flow, layering flow and annular flow. When the drill pipe rotates, each of those six flow patterns changes due to the Taylor vortex effect. Outcomes of this study could help to better understand the slotting-spurt mechanism and provide guidance on the anti-spurt strategies through eliminating the borehole deslagging blockage.展开更多
Conventional conductivity methods for measuring the void fraction in gas-liquid multiphase systems are typically affected by accuracy problems due to the presence of fluid flow and salinity.This study presents a novel...Conventional conductivity methods for measuring the void fraction in gas-liquid multiphase systems are typically affected by accuracy problems due to the presence of fluid flow and salinity.This study presents a novel approach for determining the void fraction based on a reciprocating dynamic conductivity probe used to measure the liquid film thickness under forced annular-flow conditions.The measurement system comprises a cyclone,a conductivity probe,a probe reciprocating device,and a data acquisition and processing system.This method ensures that the flow pattern is adjusted to a forced annular flow,thereby minimizing the influence of complex and variable gas-liquid flow patterns on the measurement results;Moreover,it determines the liquid film thickness solely according to circuit connectivity rather than specific conductivity values,thereby mitigating the impact of salinity.The reliability of the measurement system is demonstrated through laboratory experiments.The experimental results indicate that,in a range of gas phase superficial velocities 5–20 m/s and liquid phase superficial velocities 0.079–0.48 m/s,the maximum measurement deviation for the void fraction is 4.23%.展开更多
In this paper, we present an adaptive moving mesh algorithm for meshesof unstructured polyhedra in three space dimensions. The algorithm automaticallyadjusts the size of the elements with time and position in the phys...In this paper, we present an adaptive moving mesh algorithm for meshesof unstructured polyhedra in three space dimensions. The algorithm automaticallyadjusts the size of the elements with time and position in the physical domain to resolvethe relevant scales in multiscale physical systems while minimizing computationalcosts. The algorithm is a generalization of the moving mesh methods basedon harmonic mappings developed by Li et al. [J. Comput. Phys., 170 (2001), pp. 562-588, and 177 (2002), pp. 365-393]. To make 3D moving mesh simulations possible,the key is to develop an efficient mesh redistribution procedure so that this part willcost as little as possible comparing with the solution evolution part. Since the meshredistribution procedure normally requires to solve large size matrix equations, wewill describe a procedure to decouple the matrix equation to a much simpler blocktridiagonaltype which can be efficiently solved by a particularly designed multi-gridmethod. To demonstrate the performance of the proposed 3D moving mesh strategy,the algorithm is implemented in finite element simulations of fluid-fluid interface interactionsin multiphase flows. To demonstrate the main ideas, we consider the formationof drops by using an energetic variational phase field model which describesthe motion of mixtures of two incompressible fluids. Numerical results on two- andthree-dimensional simulations will be presented.展开更多
A graphics-processing-unit(GPU)-parallel-based computational scheme is developed to realize the competitive growth process of converging bi-crystal in two-dimensional states in the presence of forced convection condit...A graphics-processing-unit(GPU)-parallel-based computational scheme is developed to realize the competitive growth process of converging bi-crystal in two-dimensional states in the presence of forced convection conditions by coupling a multi-phase field model and a lattice Boltzmann model.The elimination mechanism in the evolution process is analyzed for the three conformational schemes constituting converging bi-crystals under pure diffusion and forced convection conditions,respectively,expanding the research of the competitive growth of columnar dendrites under melt convection conditions.The results show that the elimination mechanism for the competitive growth of converging bi-crystals of all three configurations under pure diffusion conditions follows the conventional Walton-Chalmers model.When there is forced convection with lateral flow in the liquid phase,the anomalous elimination phenomenon of unfavorable dendrites eliminating favorable dendrites occurs in the grain boundaries.In particular,the anomalous elimination phenomenon is relatively strong in conformation 1 and conformation 2 when the orientation angle of unfavorable dendrites is small,and relatively weak in conformation 3.Moreover,the presence of convection increases the tip growth rate of both favorable and unfavorable dendrites in the grain boundary.In addition,the parallelization of the multi-phase-field-lattice Boltzmann model is achieved by designing the parallel computation of the model on the GPU platform concerning the computerunified-device-architecture parallel technique,and the results show that the parallel computation of this model based on the GPU has absolute advantages,and the parallel acceleration is more obvious as the computation area increases.展开更多
Multi-phase nitrides bonded silicon carbide lintel blocks were prepared using industrial SiC(SiC≥98 mass%,3-0.5,≤0.5 and≤0.044 mm),Si powder(Si≥98 mass%,≤0.044 mm),and SiO2 micropowder(SiO2≥96 mass%,d50=0.15 pm)...Multi-phase nitrides bonded silicon carbide lintel blocks were prepared using industrial SiC(SiC≥98 mass%,3-0.5,≤0.5 and≤0.044 mm),Si powder(Si≥98 mass%,≤0.044 mm),and SiO2 micropowder(SiO2≥96 mass%,d50=0.15 pm)as raw materials,and calcium lignosulfonate as the additive,batching,mixing,and molding on a vibration pressure molding machine,drying and then firing at 1420℃for 10 h in high-purity N2.The apparent porosity,the bulk density,the cold modulus of rupture,the hot modulus of rupture,and the linear expansion coefficient of the samples were tested.The phase composition and the microstructure of the samples at different nitriding depths(50,100,and 150 mm)were analyzed by XRD and SEM.The field application effects of the blocks were studied.The results show that:(1)the multi-phase nitrides bonded silicon carbide refractories can dynamically adjust their own phase composition and minimize structural and thermal stresses,improving the service life of key parts of dry quenching furnaces;(2)calcium lignosulfonate can improve the nitriding micro-environment of multi-phase nitrides bonded silicon carbide lintel blocks,successfully increasing the effective nitriding thickness of the blocks to 300 mm;(3)Sinosteel LI RR provides a unique concept in the design of materials and block types as well as the stable and scientific overall structure,promoting the industrialization process of dry quenching furnaces with long service life in China.展开更多
The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a...The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a preliminary proppant settling test,from which a solid volume fraction of 0.575 is calibrated for the proppant pack in the fracture.In the established workflow to investigate proppant flowback,a displacement is applied to the fracture surfaces to compact the generated proppant pack as well as further mimicking proppant embedment under closure stress.When a pressure gradient is applied to drive the fluid-particle flow,a critical aperture-to-diameter ratio of 4 is observed,above which the proppant pack would collapse.The results also show that the volumetric proppant flowback rate increases quadratically with the fracture aperture,while a linear variation between the particle flux and the pressure gradient is exhibited for a fixed fracture aperture.The research outcome contributes towards an improved understanding of proppant flowback in hydraulic fractures,which also supports an optimised proppant size selection for hydraulic fracturing operations.展开更多
Severe slugging flow is always challenging in oil & gas production, especially for the current offshore based production. The slugging flow can cause a lot of problems, such as those relevant to production safety~ fa...Severe slugging flow is always challenging in oil & gas production, especially for the current offshore based production. The slugging flow can cause a lot of problems, such as those relevant to production safety~ fatigue as well as capability. As one typical phenomenon in multi-phase flow dynamics, the slug can be avoided or eliminated by proper facility d^sign or control of operational conditions. Based on a testing facility which can emulate a pipeline-riser or a gas-lifted production well in a scaled-down manner, this paper experimentally studies the correlations of key operational parameters with severe Slugging flows. These correlations are reflected through an obtained stable surface in the parameter space, which is a natural extension of the bifurcation plot. The maximal production opportunity without compromising the stability is also studied. Relevant studies have already showed that the capability, performance and efficiency of anti-slug control can be dramatically improved if these stable surfaces can be experimentally determined beforehand. The paper concludes that obtaining the stable surface on the new developed map can significantly improve the production rate in a control scheme. Even though the production rate can be further improved by moving the stable surface using advanced control strategies, the constant inputs can in some cases be preferable due to the easier implementation.展开更多
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.展开更多
This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeabi...This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.51704327)
文摘In the petroleum industry,detection of multi-phase fluid flow is very important in both surface and down-hole measurements.Accurate measurement of high rate of water or gas multi-phase flow has always been an academic and industrial focus.NMR is an efficient and accurate technique for the detection of fluids;it is widely used in the determination of fluid compositions and properties.This paper is aimed to quantitatively detect multi-phase flow in oil and gas wells and pipelines and to propose an innovative method for online nuclear magnetic resonance(NMR)detection.The online NMR data acquisition,processing and interpretation methods are proposed to fill the blank of traditional methods.A full-bore straight tube design without pressure drop,a Halbach magnet structure design with zero magnetic leakage outside the probe,a separate antenna structure design without flowing effects on NMR measurement and automatic control technology will achieve unattended operation.Through the innovation of this work,the application of NMR for the real-time and quantitative detection of multi-phase flow in oil and gas wells and pipelines can be implemented.
文摘Most current lattice Boltzmann (LBM) models suffer from the deficiency that their parameters have to be obtained by fitting experimental results. In this paper, we propose a new method that integrates the molecular dynamics (MD) simulation and LBM to avoid such defect. The basic idea is to first construct a molecular model based on the actual components of the rock-fluid system, then to compute the interaction force between the rock and the fluid of different densities through the MD simulation. This calculated rock-fluid interaction force, combined with the fluid-fluid force determined from the equation of state, is then used in LBM modeling. Without parameter fitting, this study presents a new systematic approach for pore-scale modeling of multi-phase flow. We have validated this ap- proach by simulating a two-phase separation process and gas-liquid-solid three-phase contact angle. Based on an actual X-ray CT image of a reservoir core, we applied our workflow to calculate the absolute permeability of the core, vapor-liquid H20 relative permeability, and capillary pressure curves.
文摘Hydraulic slotting in a gas drainage borehole is an effective method of enhancing gas drainage performance. However, it frequently occurs that a large amount of slotting products(mainly the coal slurry and gas) intensely spurt out of the borehole during the slotting, which adversely affects the slotting efficiency. Despite extensive previous investigations on the mechanism and prevention-device design of the spurt during ordinary borehole drilling, a very few studies has focused on the spurt in the slotting process. The slotting spurt is mainly caused by two reasons: the coal and gas outburst in the borehole and the borehole deslagging blockage. This paper focuses on the second reason, and investigates the hydraulic deslagging flow patterns in the annular space between the drill pipe and borehole wall.Results show that there are six deslagging flow patterns when the drill pipe is still: pure slurry flow, pure gas flow, bubble flow, intermittent flow, layering flow and annular flow. When the drill pipe rotates, each of those six flow patterns changes due to the Taylor vortex effect. Outcomes of this study could help to better understand the slotting-spurt mechanism and provide guidance on the anti-spurt strategies through eliminating the borehole deslagging blockage.
基金the National Natural Science Foundation of China(No.62173049)the Open Fund of the Hubei Key Laboratory of Oil and Gas Drilling and Production Engineering(Yangtze University),YQZC202309.
文摘Conventional conductivity methods for measuring the void fraction in gas-liquid multiphase systems are typically affected by accuracy problems due to the presence of fluid flow and salinity.This study presents a novel approach for determining the void fraction based on a reciprocating dynamic conductivity probe used to measure the liquid film thickness under forced annular-flow conditions.The measurement system comprises a cyclone,a conductivity probe,a probe reciprocating device,and a data acquisition and processing system.This method ensures that the flow pattern is adjusted to a forced annular flow,thereby minimizing the influence of complex and variable gas-liquid flow patterns on the measurement results;Moreover,it determines the liquid film thickness solely according to circuit connectivity rather than specific conductivity values,thereby mitigating the impact of salinity.The reliability of the measurement system is demonstrated through laboratory experiments.The experimental results indicate that,in a range of gas phase superficial velocities 5–20 m/s and liquid phase superficial velocities 0.079–0.48 m/s,the maximum measurement deviation for the void fraction is 4.23%.
基金the Joint Applied Mathematics Research Institute of Peking University and Hong Kong Baptist University.Li was also partially supported by the National Basic Research Program of China under the grant 2005CB321701The research of Tang was supported by CERG Grants of Hong Kong Research Grant Council,FRG grants of Hong Kong Baptist University,and NSAF Grant#10476032 of National Science Foundation of China.He was supported in part by the Chinese Academy of Sciences while visiting its Institute of Computational Mathematics.
文摘In this paper, we present an adaptive moving mesh algorithm for meshesof unstructured polyhedra in three space dimensions. The algorithm automaticallyadjusts the size of the elements with time and position in the physical domain to resolvethe relevant scales in multiscale physical systems while minimizing computationalcosts. The algorithm is a generalization of the moving mesh methods basedon harmonic mappings developed by Li et al. [J. Comput. Phys., 170 (2001), pp. 562-588, and 177 (2002), pp. 365-393]. To make 3D moving mesh simulations possible,the key is to develop an efficient mesh redistribution procedure so that this part willcost as little as possible comparing with the solution evolution part. Since the meshredistribution procedure normally requires to solve large size matrix equations, wewill describe a procedure to decouple the matrix equation to a much simpler blocktridiagonaltype which can be efficiently solved by a particularly designed multi-gridmethod. To demonstrate the performance of the proposed 3D moving mesh strategy,the algorithm is implemented in finite element simulations of fluid-fluid interface interactionsin multiphase flows. To demonstrate the main ideas, we consider the formationof drops by using an energetic variational phase field model which describesthe motion of mixtures of two incompressible fluids. Numerical results on two- andthree-dimensional simulations will be presented.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.52161002,51661020,and 11364024)the Postdoctoral Science Foundation of China(Grant No.2014M560371)the Funds for Distinguished Young Scientists of Lanzhou University of Technology,China(Grant No.J201304).
文摘A graphics-processing-unit(GPU)-parallel-based computational scheme is developed to realize the competitive growth process of converging bi-crystal in two-dimensional states in the presence of forced convection conditions by coupling a multi-phase field model and a lattice Boltzmann model.The elimination mechanism in the evolution process is analyzed for the three conformational schemes constituting converging bi-crystals under pure diffusion and forced convection conditions,respectively,expanding the research of the competitive growth of columnar dendrites under melt convection conditions.The results show that the elimination mechanism for the competitive growth of converging bi-crystals of all three configurations under pure diffusion conditions follows the conventional Walton-Chalmers model.When there is forced convection with lateral flow in the liquid phase,the anomalous elimination phenomenon of unfavorable dendrites eliminating favorable dendrites occurs in the grain boundaries.In particular,the anomalous elimination phenomenon is relatively strong in conformation 1 and conformation 2 when the orientation angle of unfavorable dendrites is small,and relatively weak in conformation 3.Moreover,the presence of convection increases the tip growth rate of both favorable and unfavorable dendrites in the grain boundary.In addition,the parallelization of the multi-phase-field-lattice Boltzmann model is achieved by designing the parallel computation of the model on the GPU platform concerning the computerunified-device-architecture parallel technique,and the results show that the parallel computation of this model based on the GPU has absolute advantages,and the parallel acceleration is more obvious as the computation area increases.
文摘Multi-phase nitrides bonded silicon carbide lintel blocks were prepared using industrial SiC(SiC≥98 mass%,3-0.5,≤0.5 and≤0.044 mm),Si powder(Si≥98 mass%,≤0.044 mm),and SiO2 micropowder(SiO2≥96 mass%,d50=0.15 pm)as raw materials,and calcium lignosulfonate as the additive,batching,mixing,and molding on a vibration pressure molding machine,drying and then firing at 1420℃for 10 h in high-purity N2.The apparent porosity,the bulk density,the cold modulus of rupture,the hot modulus of rupture,and the linear expansion coefficient of the samples were tested.The phase composition and the microstructure of the samples at different nitriding depths(50,100,and 150 mm)were analyzed by XRD and SEM.The field application effects of the blocks were studied.The results show that:(1)the multi-phase nitrides bonded silicon carbide refractories can dynamically adjust their own phase composition and minimize structural and thermal stresses,improving the service life of key parts of dry quenching furnaces;(2)calcium lignosulfonate can improve the nitriding micro-environment of multi-phase nitrides bonded silicon carbide lintel blocks,successfully increasing the effective nitriding thickness of the blocks to 300 mm;(3)Sinosteel LI RR provides a unique concept in the design of materials and block types as well as the stable and scientific overall structure,promoting the industrialization process of dry quenching furnaces with long service life in China.
基金Funding support from Heilongjiang"Open Competition"project(Grant No.DQYT2022-JS-758)is greatly acknowledgedFinancial support from the National Natural Science Foundation of China(Grant Nos.52304025 and 52174025)is acknowledged+1 种基金supports from Northeast Petroleum University and Guangdong Basic and Applied Basic Research Foundationsupport from the Heilongjiang Touyan Innovation Team Program.
文摘The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a preliminary proppant settling test,from which a solid volume fraction of 0.575 is calibrated for the proppant pack in the fracture.In the established workflow to investigate proppant flowback,a displacement is applied to the fracture surfaces to compact the generated proppant pack as well as further mimicking proppant embedment under closure stress.When a pressure gradient is applied to drive the fluid-particle flow,a critical aperture-to-diameter ratio of 4 is observed,above which the proppant pack would collapse.The results also show that the volumetric proppant flowback rate increases quadratically with the fracture aperture,while a linear variation between the particle flux and the pressure gradient is exhibited for a fixed fracture aperture.The research outcome contributes towards an improved understanding of proppant flowback in hydraulic fractures,which also supports an optimised proppant size selection for hydraulic fracturing operations.
基金supported by Innovation Fund Denmark through the PDPWAC Project(No.95-2012-3)
文摘Severe slugging flow is always challenging in oil & gas production, especially for the current offshore based production. The slugging flow can cause a lot of problems, such as those relevant to production safety~ fatigue as well as capability. As one typical phenomenon in multi-phase flow dynamics, the slug can be avoided or eliminated by proper facility d^sign or control of operational conditions. Based on a testing facility which can emulate a pipeline-riser or a gas-lifted production well in a scaled-down manner, this paper experimentally studies the correlations of key operational parameters with severe Slugging flows. These correlations are reflected through an obtained stable surface in the parameter space, which is a natural extension of the bifurcation plot. The maximal production opportunity without compromising the stability is also studied. Relevant studies have already showed that the capability, performance and efficiency of anti-slug control can be dramatically improved if these stable surfaces can be experimentally determined beforehand. The paper concludes that obtaining the stable surface on the new developed map can significantly improve the production rate in a control scheme. Even though the production rate can be further improved by moving the stable surface using advanced control strategies, the constant inputs can in some cases be preferable due to the easier implementation.
基金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.
文摘This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.