Accurately characterizing the pore structure of Tamusu mudstone and simulating fluid flow within it are crucial for assessing underground disposal of high-level radioactive waste(HLW)in this formation.Modeling mudston...Accurately characterizing the pore structure of Tamusu mudstone and simulating fluid flow within it are crucial for assessing underground disposal of high-level radioactive waste(HLW)in this formation.Modeling mudstone presents challenges due to its multiscale pore structure and the necessity that accounts for the effects of high clay content.A method for constructing a dual-scale pore network model(PNM)for the Tamusu mudstone,which considers the hydrological expansion of clays,was proposed.This model integrates N2 adsorption data with focused ion beam/scanning electron microscopy(FIB/SEM)images and labels pores based on clay content.Simulations of single-phase flow were conducted to validate the proposed model.Additionally,the influences of cell number,connectivity,slip effects,and clay minerals on permeability were examined.The findings indicate that a configuration of 45×45×45 cells adequately represents the model.The permeability of the Tamusu mudstone,about 1020 m^(2),aligns with the experimental values.During the simulation,Knudsen diffusion is considered.Factors such as increased roughness,tortuosity,clay content,and water film thickness decrease the permeability,whereas increased connectivity enhances permeability.In the model,numerical coordination numbers 2 and 3 are deemed suitable for the Tamusu mudstone.The proposed model is effective as a tool for constructing and simulating fluid flow in the Tamusu mudstone.展开更多
In oceanic and atmospheric science,finer resolutions have become a prevailing trend in all aspects of development.For high-resolution fluid flow simulations,the computational costs of widely used numerical models incr...In oceanic and atmospheric science,finer resolutions have become a prevailing trend in all aspects of development.For high-resolution fluid flow simulations,the computational costs of widely used numerical models increase significantly with the resolution.Artificial intelligence methods have attracted increasing attention because of their high precision and fast computing speeds compared with traditional numerical model methods.The resolution-independent Fourier neural operator(FNO)presents a promising solution to the still challenging problem of high-resolution fluid flow simulations based on low-resolution data.Accordingly,we assess the potential of FNO for high-resolution fluid flow simulations using the vorticity equation as an example.We assess and compare the performance of FNO in multiple high-resolution tests varying the amounts of data and the evolution durations.When assessed with finer resolution data(even up to number of grid points with 1280×1280),the FNO model,trained at low resolution(number of grid points with 64×64)and with limited data,exhibits a stable overall error and good accuracy.Additionally,our work demonstrates that the FNO model takes less time than the traditional numerical method for high-resolution simulations.This suggests that FNO has the prospect of becoming a cost-effective and highly precise model for high-resolution simulations in the future.Moreover,FNO can make longer high-resolution predictions while training with less data by superimposing vorticity fields from previous time steps as input.A suitable initial learning rate can be set according to the frequency principle,and the time intervals of the dataset need to be adjusted according to the spatial resolution of the input when training the FNO model.Our findings can help optimize FNO for future fluid flow simulations.展开更多
With the widespread adoption of hydraulic fracturing technology in oil and gas resource development,improving the accuracy and efficiency of fracturing simulations has become a critical research focus.This paper propo...With the widespread adoption of hydraulic fracturing technology in oil and gas resource development,improving the accuracy and efficiency of fracturing simulations has become a critical research focus.This paper proposes an improved fluid flow algorithm,aiming to enhance the computational efficiency of hydraulic fracturing simulations while ensuring computational accuracy.The algorithm optimizes the aperture law and iteration criteria,focusing on improving the domain volume and crack pressure update strategy,thereby enabling precise capture of dynamic borehole pressure variations during injection tests.The effectiveness of the algorithm is verified through three flow-solid coupling cases.The study also analyzes the effects of borehole size,domain volume,and crack pressure update strategy on fracturing behavior.Furthermore,the performance of the improved algorithm in terms of crack propagation rate,micro-crack formation,and fluid pressure distribution was further evaluated.The results indicate that while large-size boreholes delay crack initiation,the cracks propagate more rapidly once formed.Additionally,the optimized domain volume calculation and crack pressure update strategy significantly shorten the pressure propagation stage,promote crack propagation,and improve computational efficiency.展开更多
The wave-induced fluid flow(WIFF) occurring in the ubiquitous layered porous media(e.g.,shales)usually causes the appreciable seismic energy dissipation,which further leads to the frequency dependence of wave velocity...The wave-induced fluid flow(WIFF) occurring in the ubiquitous layered porous media(e.g.,shales)usually causes the appreciable seismic energy dissipation,which further leads to the frequency dependence of wave velocity(i.e.,dispersion) and elastic anisotropy parameters.The relevant knowledge is of great importance for geofluid discrimination and hydrocarbon exploration in the porous shale reservoirs.We derive the wave equations for a periodic layered transversely isotropy medium with a vertical axis of symmetry(VTI) concurrently with the annular cracks(PLPC medium) based on the periodic-layered model and anisotropic Biot's theory,which simultaneously incorporate the effects of microscopic squirt fluid flow,mesoscopic interlayer fluid flow and macroscopic global fluid flow.Notably,the microscopic squirt shorten fluid flow emerges between the annular-shaped cracks and stiff pores,which generates one attenuation peak.Specifically,we first establish the stress-strain relationship and pore fluid pressure in a PLPC medium,and then use them to derive the wave equations by means of the Newton's second law.The plane analysis is implemented on the wave equations to yield the analytic solutions for phase velocities and attenuation factors of four waves,namely,fast P-wave,slow P-wave,SV-wave and SH-wave,and the anisotropy parameters can be therefore computed.Simulation results show that P-wave velocity have three attenuation peaks throughout the full frequency band,which respectively correspond to the influences of interlayer flow,the squirt flow and the Biot flow.Through the results of seismic velocity dispersion and attenuation at different incident angles,we find that the WIFF mechanism also has a significant impact on the dispersion characteristics of elastic anisotropy parameters within the low-mid frequency band.Moreover,it is shown that several poroelastic parameters,such as layer thickness ratio,crack aspect ratio and crack density have notable influence on seismic dispersion and attenuation.We compare the proposed modeled velocities with that given by the existing theory to confirm its validity.Our formulas and result can provide a better understanding of wave propagation in PLPC medium by considering the unified impacts of micro-,meso-and macro-scale WIFF mechanisms,which potentially lays a theoretical basis of rock physics for seismic interpretation.展开更多
A 1∶8 physical water model was constructed to investigate the fluid flow and mixing phenomena in the basic oxygen furnace(BOF)converter.The particle image velocimetry was employed to measure the velocity distribution...A 1∶8 physical water model was constructed to investigate the fluid flow and mixing phenomena in the basic oxygen furnace(BOF)converter.The particle image velocimetry was employed to measure the velocity distribution of the bath and the high-speed camera was applied to capture the cavity shape in the combined blowing BOF converter.The mixing time for varied operating conditions was measured by the stimulus-response approach.The cavity depth increased with the decrease in the lance height and the increase in the top gas flow rate while the bottom blowing gas had little influence on the cavity depth.The minimum cavity depth was obtained under the condition of a 69.8 m^(3)/h top gas flow rate,a 287.5 mm lance height and a 0.93 m^(3)/h bottom blowing gas flow rate,which was 161.2 mm.The mixing time decreased as the lance height decreased and the top blowing gas flow rate increased.The mixing time was first decreased and then increased with the increase in the bottom gas flow rate.With the condition of 69.8 m^(3)/h gas flow rate of top blowing,the 287.5 mm lance height and the 0.93 m^(3)/h gas flow rate of bottom blowing,the mixing time in the converter was 48.65 s.The empirical formula between the stirring power and the mixing time in the converter was calculated.展开更多
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.展开更多
The unstable fluid flow and severe free surface fluctuations in the wheel and belt caster can affect the quality of the cast bar.The lower level height tends to entrap inclusions in the molten metal.On the other hand,...The unstable fluid flow and severe free surface fluctuations in the wheel and belt caster can affect the quality of the cast bar.The lower level height tends to entrap inclusions in the molten metal.On the other hand,the higher level height makes the production process more dangerous due to the overflow of high temperature fluid from the mold.A computational model of the molten metal pouring process was established.The transient fluid flow and free surface fluctuations behavior were calculated using the three-dimensional large eddy simulation model and the volume of fluid model.The results show that the flow velocity of the main jet gradually decreases under the influence of the low kinetic energy fluid in the mold.There is an obvious oscillation in the tail of the jet,while the flow field is asymmetric in space.The jet is closer to the inside radius side due to the Coanda effect,and there is a recirculation zone on the inside radius and the outside radius respectively,according to the 10 s time-averaged results.Compared with the industrial observation and simulation results,the shape of the free surface is a wave that varies with time.In addition,the free surface height is lowest and the flow velocity is highest in the region near the jet.展开更多
Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on t...Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.展开更多
The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore st...The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore structure of the media, was constructed. With this model, three dimensional pore scale fluid flow among particles was simulated. Then the distributions of fluid flow velocity and pressure were analyzed and the hydraulic conductivity was calculated. The simulation results indicate the fluid flow behaviors are mainly dominated by the volume and topological structure of pore space. There exist obvious preferential flow and leaching blind zones simultaneously in the medium. The highest velocities generally occur in those narrow pores with high pressure drops. The hydraulic conductivity obtained by simulation is the same order of magnitude as the laboratory test result, which denotes the validity of the model. The pore-scale and macro-scale are combined and the established geometrical model can be used for the simulations of other phenomena during heap leaching process.展开更多
A two-equation K-ε turbulent fluid flow model is built to model the heat transfer and fluid flow in gas tungsten arc welding （GTAW） process of stainless steel S US310 and S US316. This model combines the buoyancy f...A two-equation K-ε turbulent fluid flow model is built to model the heat transfer and fluid flow in gas tungsten arc welding （GTAW） process of stainless steel S US310 and S US316. This model combines the buoyancy force, lorentz force and marangni force as the driving forces of thefluidflow in the weld pool. The material properties are functions of temperature in this model. The simulated results show that the molten metal flowing outward is mainly caused by the marangoni convection, which makes the weld pool become wider and shallower. The comparison of the weld pool shape of SUS310 and SUS316 shows that the slight differences of the value of thermal conductivity mainly attributes to the difference of the weld pool shape and the distinction of heat transport in laminar and turbulent model makes large diversity in the simulated results.展开更多
The bottom-simulating reflector（BSR） is a seismic indicator of the bottom of a gas hydrate stability zone. Its depth can be used to calculate the seafloor surface heat flow. The calculated BSR heat flow variations i...The bottom-simulating reflector（BSR） is a seismic indicator of the bottom of a gas hydrate stability zone. Its depth can be used to calculate the seafloor surface heat flow. The calculated BSR heat flow variations include disturbances from two important factors：（1） seafloor topography, which focuses the heat flow over regions of concave topography and defocuses it over regions of convex topography, and（2） the focused warm fluid flow within the accretionary prism coming from depths deeper than BSR. The focused fluid flow can be detected if the contribution of the topography to the BSR heat flow is removed. However, the analytical equation cannot solve the topographic effect at complex seafloor regions. We prove that 3D finite element method can model the topographic effect on the regional background heat flow with high accuracy, which can then be used to correct the topographic effect and obtain the BSR heat flow under the condition of perfectly flat topography. By comparing the corrected BSR heat flow with the regional background heat flow, focused fluid flow regions can be detected that are originally too small and cannot be detected using present-day equipment. This method was successfully applied to the midslope region of northern Cascadia subducting margin. The results suggest that the Cucumber Ridge and its neighboring area are positive heat flow anomalies, about 10%–20% higher than the background heat flow after 3D topographic correction. Moreover, the seismic imaging associated the positive heat flow anomaly areas with seabed fracture–cavity systems. This suggests flow of warm gas-carrying fluids along these high-permeability pathways, which could result in higher gas hydrate concentrations.展开更多
The continuum approach in fluid flow modeling is generally applied to porous geological media, but has limited applicability to fractured rocks. With the presence of a discrete fracture network relatively sparsely dis...The continuum approach in fluid flow modeling is generally applied to porous geological media, but has limited applicability to fractured rocks. With the presence of a discrete fracture network relatively sparsely distributed in the matrix, it may be difficult or erroneous to use a porous medium fluid flow model with continuum assumptions to describe the fluid flow in fractured rocks at small or even large field scales. A discrete fracture fluid flow approach incorporating a stochastic fracture network with numerical fluid flow simulations could have the capability of capturing fluid flow behaviors such as inhomogeneity and anisotropy while reflecting the changes of hydraulic features at different scales. Moreover, this approach can be implemented to estimate the size of the representative elementary volume (REV) in order to find out the scales at which a porous medium flow model could be applied, and then to determine the hydraulic conductivity tensor for fractured rocks. The following topics are focused on in this study: (a) conceptual discrete fracture fluid flow modeling incorporating a stochastic fracture network with numerical flow simulations; (b) estimation of REV and hydraulic conductivity tensor for fractured rocks utilizing a stochastic fracture network with numerical fluid flow simulations; (c) investigation of the effect of fracture orientation and density on the hydraulic conductivity and REV by implementing a stochastic fracture network with numerical fluid flow simulations, and (d) fluid flow conceptual models accounting for major and minor fractures in the 2 D or 3 D flow fields incorporating a stochastic fracture network with numerical fluid flow simulations.展开更多
Local fluid flow（LFF） at the mesoscopic scale is the main dissipation mechanism of seismic waves in heterogeneous porous media within the seismic frequency band.LFF is easily influenced by the structure and boundary...Local fluid flow（LFF） at the mesoscopic scale is the main dissipation mechanism of seismic waves in heterogeneous porous media within the seismic frequency band.LFF is easily influenced by the structure and boundary conditions of the porous media,which leads to different behaviors of the peak frequency of attenuation.The associated transition frequency can provide detailed information about the trend of LFF;therefore,research on the transition frequency of LFF and its relationship with the peak frequency of the corresponding attenuation（i.e.,inverse of quality factor） facilitates the detailed understanding of the effect of inner structures and boundary conditions in porous media.In this study,we firstly obtain the transition frequency of fluid flux based on Biot＇s theory of poroelasticity and the fast Fourier transform algorithm in a sample containing one repeating unit cell（RUC）.We then analyze changes of these two frequencies in porous media with different porous properties.Finally,we extend our analysis to the influence of the undrained boundary condition on the transition frequency and peak frequency in porous media with multiple RUCs.This setup can facilitate the understanding of the effect from the undrained boundary condition.Results demonstrate that these two frequencies have the same trend at low water saturation,but amplitude variations differ between the frequencies as the amount of saturation increases.However,for cases of high water saturation,both the trend and the amplitude variation of these two frequencies fit well with each other.展开更多
This paper presents the results of a set of numerical models focussing on structural controls on hydrothermal mineralization. We first give an overview of natural phenomena of structurally-controlled ore formation and...This paper presents the results of a set of numerical models focussing on structural controls on hydrothermal mineralization. We first give an overview of natural phenomena of structurally-controlled ore formation and the background theory and mechanisms for such controls. We then provide the results of a group of simple 2D numerical models validated through comparison with Cu-vein structure observed near the Shilu Copper deposit （Yangchun, Guangdong Province, China） and finally a case study of 3D numerical modelling applied to the Hodgkinson Province in North Queensland （Australia）. Two modelling approaches, discrete deformation modelling and continuum coupled deformation and fluid flow modelling, are involved. The 2D model-derived patterns are remarkably consistent with the Cu-vein structure from the Shilu Copper deposit, and show that both modelling approaches can realistically simulate the mechanical behaviours of shear and dilatant fractures. The continuum coupled deformation and fluid flow model indicates that pattern of the Cu- veins near the Shilu deposit is the result of shear strain localization, development of dilation and fluid focussing into the dilatant fracture segments. The 3D case-study models （with deformation and fluid flow coupling） on the Hodgkinson Province generated a number of potential gold mineralization展开更多
The Ordos Basin of North China is not only an important uranium mineralization province, but also a major producer of oil, gas and coal in China. The genetic relationship between uranium mineralization and hydrocarbon...The Ordos Basin of North China is not only an important uranium mineralization province, but also a major producer of oil, gas and coal in China. The genetic relationship between uranium mineralization and hydrocarbons has been recognized by a number of previous studies, but it has not been well understood in terms of the hydrodynamics of basin fluid flow. We have demonstrated in a previous study that the preferential localization of Cretaceous uranium mineralization in the upper part of the Ordos Jurassic section may have been related to the interface between an upward flowing, reducing fluid and a downward flowing, oxidizing fluid. This interface may have been controlled by the interplay between fluid overpressure related to disequilibrium sediment compaction and which drove the upward flow, and topographic relief, which drove the downward flow. In this study, we carried out numerical modeling for the contribution of oil and gas generation to the development of fluid overpressure, in addition to sedi- ment compaction and heating. Our results indicate that when hydrocarbon generation is taken into account, fluid overpressure during the Cretaceous was more than doubled in comparison with the simu- lation when hydrocarbon generation was not considered. Furthermore, fluid overpressure dissipation at the end of sedimentation slowed down relative to the no-hydrocarbon generation case. These results suggest that hydrocarbon generation may have played an important role in uranium mineralization, not only in providing reducing agents required for the mineralization, but also in contributing to the driving force to maintain the upward flow.展开更多
The Precambrian Dengying Formation is a set of large-scale, extensively dolomitized, carbonate reservoirs occurring within the Sichuan Basin. Petrographic and geochemical studies reveal dolomitization was a direct res...The Precambrian Dengying Formation is a set of large-scale, extensively dolomitized, carbonate reservoirs occurring within the Sichuan Basin. Petrographic and geochemical studies reveal dolomitization was a direct result of precipitation by chemically distinct fluids occurring at different times and at different intensities. Based on this evidence, dolomitization and multiple fluid flow events are analyzed, and three types of fluid evolution models are proposed. Results of analysis show that Precambrian Dengying Formation carbonates were deposited in a restricted peritidal environment （630- 542 Ma）. A high temperature and high Mg2＋ concentration seawater was a direct result of dolomitization for the micrite matrix, and for fibrous aragonite in primary pores. Geochemical evidence shows low δ18O values of micritic dolomite varying from -1.29%o to -4.52%o PDB, abundant light rare earth elements （REEs）, and low dolomite order degrees. Microbes and meteoric water significantly altered dolomite original chemical signatures, resulting in algal micritic dolomite and the fine-grained, granular, dolosparite dolomite having very negative δ18O values. Finely crystalline cement dolomite （536.3-280 Ma） and coarsely crystalline cement dolomite have a higher crystallization degree and higher order degree. The diagenetic sequence and fluid inclusion evidence imply a linear correlation between their burial depth and homogenization temperatures, which closely resemble the temperature of generated hydrocarbon. Compared with finely crystalline dolomite, precipitation of coarsely crystalline dolomite was more affected by restricted basinal fluids. In addition, there is a trend toward a more negative δ18O value, higher salinity, higher Fe and Mn concentrations, REE-rich. Two periods of hydrothermal fluids are identified, as the exceptionally high temperatures as opposed to the temperatures of burial history, in addition to the presence of high salinity fluid inclusions. The early hydrothermal fluid flow event was characterized by hot magnesium- and silicon-rich fluids, as demonstrated by the recrystallized matrix dolomite that is intimately associated with flint, opal, and microcrystalline quartz in intergranular or intercrystalline pores. This event was likely the result of a seafloor hydrothermal chimney eruption during Episode I of the Tongwan Movement （536.3-5.5 Ma）. In contrast, later hydrothermal fluids, which caused precipitation of saddle dolomite, were characterized by high salinity （15-16.05wt% NaCI equivalent） and homogenization temperatures （250 to 265℃）, δ18O values that were more enriched, and REE signatures. Geochemical data and the paragenetic sequence indicate that this hydrothermal fluid was related to extensive Permian large igneous province activity （360-280 Ma）. This study demonstrates the presence of complicated dolomitization processes occurring during various paleoclimates, tectonic cycles, and basinal fluids flow; results are a useful reference for these dolomitized Precambrian carbonates reservoirs.展开更多
The special gas wettability phenomenon of reservoir rocks has been recognized by more and more researchers.It has a significant effect on efficient development of unconventional reservoirs.First,based on the preferent...The special gas wettability phenomenon of reservoir rocks has been recognized by more and more researchers.It has a significant effect on efficient development of unconventional reservoirs.First,based on the preferentially gas-covered ability and surface free energy changes,definition and evaluation methods have been established.Second,a method for altering rock wettability and its mechanisms have been studied,surface oriented phenomena of functional groups with low surface energy are the fundamental reason for gas wettability alteration of rock.Third,the effect of gas wettability on the surface energy,electrical properties and dilatability are investigated.Last,the effects of gas wettability on capillary pressure,oil/gas/water distribution and flow are investigated with capillary tubes and etchedglass network models.The gas wettability theory of reservoir rocks has been initially established,which provides theoretical support for the efficient production of unconventional reservoirs and has great significance.展开更多
Following up the fluid flow simulation in a 60 t tundish, the trajectories of inclusions in the 60 t tundish without flow control are simulated by considering the force balance between the drag force and the inertial ...Following up the fluid flow simulation in a 60 t tundish, the trajectories of inclusions in the 60 t tundish without flow control are simulated by considering the force balance between the drag force and the inertial buoyancy force. The Stochastic model yields more accurate inclusion motion than the non-Stochastic model due to including the effect of the turbulent fluctuation. The average residence time of inclusions decreases with increasing size. The thermal buoyancy favors inclusions removal especially the small inclusions. Using solute transport like the dye injection in water model and copper addition in the real steel tundish cannot accurately study the motion of the inclusions. In the simulation, more than 68% inclusions bigger than 10μm are removed to the top, and less than 32% enters the mold. The thermal buoyancy has little effect on the fraction of inclusions moved to the top of the inlet zone, and it mainly favors the removal of inclusions smaller than 100μm to the top surface of the outlet zone. For inclusions bigger than 100μm, the effect of thermal buoyancy on their motion can be ignored compared to the inertial buoyancy effect.展开更多
The origin and migration of natural gas and the accumulation of gas hydrates within the Pearl River Mouth Basin of the northern South China Sea are poorly understood. Based on high-resolution 2D/3D seismic data, three...The origin and migration of natural gas and the accumulation of gas hydrates within the Pearl River Mouth Basin of the northern South China Sea are poorly understood. Based on high-resolution 2D/3D seismic data, three environments of focused fluid flow： gas chimneys, mud diapirs and active faults have been identified. Widespread gas chimneys that act as important conduits for fluid flow are located below bottom simulating reflections and above basal uplifts. The occurrence and evolution of gas chimneys can be divided into a violent eruptive stage and a quiet seepage stage. For most gas chimneys, the strong eruptions are deduced to have happened during the Dongsha Movement in the latest Miocene, which are observed below Pliocene strata and few active faults develop above the top of the Miocene. The formation pressures of the Baiyun Sag currently are considered to be normal, based on these terms： 1） Borehole pressure tests with pressure coefficients of 1.043-l.047; 2） The distribution of gas chimneys is limited to strata older than the Pliocene; 3） Disseminated methane hydrates, rather than fractured hydrates, are found in the hydrate samples; 4） The gas hydrate is mainly charged with biogenic gas rather than thermogenic gas based on the chemical tests from gas hydrates cores. However, periods of quiet focused fluid flow also enable the establishment of good conduits for the migration of abundant biogenic gas and lesser volumes ofthermogenic gas. A geological model goveming fluid flow has been proposed to interpret the release of overpressure, the migration of fluids and the formation of gas hydrates, in an integrated manner. This model suggests that gas chimneys positioned above basal uplifts were caused by the Dongsha Movement at about 5.5 Ma. Biogenic gas occupies the strata above the base of the middle Miocene and migrates slowly into the gas chimney columns. Some of the biogenic gas and small volumes ofthermogenic gas eventually contribute to the formation of the gas hydrates.展开更多
Understanding the behaviors of heat transfer and fluid flow in weld pool and their effects on the solidification microstructure are significant for performance improvement of laser welds.This paper develops a three-di...Understanding the behaviors of heat transfer and fluid flow in weld pool and their effects on the solidification microstructure are significant for performance improvement of laser welds.This paper develops a three-dimensional numerical model to understand the multi-physical processes such as heat transfer,melt convection and solidification behavior in full-penetration laser welding of thin 5083 aluminum sheet.Solidification parameters including temperature gradient G and solidification rate R,and their combined forms are evaluated to interpret solidification microstructure.The predicted weld dimensions and the microstructure morphology and scale agree well with experiments.Results indicate that heat conduction is the dominant mechanism of heat transfer in weld pool,and melt convection plays a critical role in microstructure scale.The mushy zone shape/size and solidification parameters can be modulated by changing process parameters.Dendritic structures form because of the low G/R value.The scale of dendritic structures can be reduced by increasing GR via decreasing heat input.The columnar to equiaxed transition is predicted quantitatively via the process related G^3/R.These findings illustrate how heat transfer and fluid flow affect the solidification parameters and hence the microstructure,and show how to improve microstructure by optimizing the process.展开更多
基金support of the National Natural Science Foundation of China(Grant Nos.42377179,U22A20595,12202463).
文摘Accurately characterizing the pore structure of Tamusu mudstone and simulating fluid flow within it are crucial for assessing underground disposal of high-level radioactive waste(HLW)in this formation.Modeling mudstone presents challenges due to its multiscale pore structure and the necessity that accounts for the effects of high clay content.A method for constructing a dual-scale pore network model(PNM)for the Tamusu mudstone,which considers the hydrological expansion of clays,was proposed.This model integrates N2 adsorption data with focused ion beam/scanning electron microscopy(FIB/SEM)images and labels pores based on clay content.Simulations of single-phase flow were conducted to validate the proposed model.Additionally,the influences of cell number,connectivity,slip effects,and clay minerals on permeability were examined.The findings indicate that a configuration of 45×45×45 cells adequately represents the model.The permeability of the Tamusu mudstone,about 1020 m^(2),aligns with the experimental values.During the simulation,Knudsen diffusion is considered.Factors such as increased roughness,tortuosity,clay content,and water film thickness decrease the permeability,whereas increased connectivity enhances permeability.In the model,numerical coordination numbers 2 and 3 are deemed suitable for the Tamusu mudstone.The proposed model is effective as a tool for constructing and simulating fluid flow in the Tamusu mudstone.
基金The National Natural Science Foundation of China under contract No.42425606the Basic Scientific Fund for the National Public Research Institute of China(Shu-Xingbei Young Talent Program)under contract No.2023S01+1 种基金the Ocean Decade International Cooperation Center Scientific and Technological Cooperation Project under contract No.GHKJ2024005China-Korea Joint Ocean Research Center Project under contract Nos PI-20240101(China)and 20220407(Korea).
文摘In oceanic and atmospheric science,finer resolutions have become a prevailing trend in all aspects of development.For high-resolution fluid flow simulations,the computational costs of widely used numerical models increase significantly with the resolution.Artificial intelligence methods have attracted increasing attention because of their high precision and fast computing speeds compared with traditional numerical model methods.The resolution-independent Fourier neural operator(FNO)presents a promising solution to the still challenging problem of high-resolution fluid flow simulations based on low-resolution data.Accordingly,we assess the potential of FNO for high-resolution fluid flow simulations using the vorticity equation as an example.We assess and compare the performance of FNO in multiple high-resolution tests varying the amounts of data and the evolution durations.When assessed with finer resolution data(even up to number of grid points with 1280×1280),the FNO model,trained at low resolution(number of grid points with 64×64)and with limited data,exhibits a stable overall error and good accuracy.Additionally,our work demonstrates that the FNO model takes less time than the traditional numerical method for high-resolution simulations.This suggests that FNO has the prospect of becoming a cost-effective and highly precise model for high-resolution simulations in the future.Moreover,FNO can make longer high-resolution predictions while training with less data by superimposing vorticity fields from previous time steps as input.A suitable initial learning rate can be set according to the frequency principle,and the time intervals of the dataset need to be adjusted according to the spatial resolution of the input when training the FNO model.Our findings can help optimize FNO for future fluid flow simulations.
基金supported by the National Natural Science Foundation of China(Nos.52164001,52064006,52004072 and 52364004)the Science and Technology Support Project of Guizhou(Nos.[2020]4Y044,[2021]N404 and[2021]N511)+1 种基金the Guizhou Provincial Science and Technology Foundation(No.GCC[2022]005-1),Talents of Guizhou University(No.201901)the Special Research Funds of Guizhou University(Nos.201903,202011,and 202012).
文摘With the widespread adoption of hydraulic fracturing technology in oil and gas resource development,improving the accuracy and efficiency of fracturing simulations has become a critical research focus.This paper proposes an improved fluid flow algorithm,aiming to enhance the computational efficiency of hydraulic fracturing simulations while ensuring computational accuracy.The algorithm optimizes the aperture law and iteration criteria,focusing on improving the domain volume and crack pressure update strategy,thereby enabling precise capture of dynamic borehole pressure variations during injection tests.The effectiveness of the algorithm is verified through three flow-solid coupling cases.The study also analyzes the effects of borehole size,domain volume,and crack pressure update strategy on fracturing behavior.Furthermore,the performance of the improved algorithm in terms of crack propagation rate,micro-crack formation,and fluid pressure distribution was further evaluated.The results indicate that while large-size boreholes delay crack initiation,the cracks propagate more rapidly once formed.Additionally,the optimized domain volume calculation and crack pressure update strategy significantly shorten the pressure propagation stage,promote crack propagation,and improve computational efficiency.
基金sponsorship of the National Natural Science Foundation of China (U24B2020,42174139)。
文摘The wave-induced fluid flow(WIFF) occurring in the ubiquitous layered porous media(e.g.,shales)usually causes the appreciable seismic energy dissipation,which further leads to the frequency dependence of wave velocity(i.e.,dispersion) and elastic anisotropy parameters.The relevant knowledge is of great importance for geofluid discrimination and hydrocarbon exploration in the porous shale reservoirs.We derive the wave equations for a periodic layered transversely isotropy medium with a vertical axis of symmetry(VTI) concurrently with the annular cracks(PLPC medium) based on the periodic-layered model and anisotropic Biot's theory,which simultaneously incorporate the effects of microscopic squirt fluid flow,mesoscopic interlayer fluid flow and macroscopic global fluid flow.Notably,the microscopic squirt shorten fluid flow emerges between the annular-shaped cracks and stiff pores,which generates one attenuation peak.Specifically,we first establish the stress-strain relationship and pore fluid pressure in a PLPC medium,and then use them to derive the wave equations by means of the Newton's second law.The plane analysis is implemented on the wave equations to yield the analytic solutions for phase velocities and attenuation factors of four waves,namely,fast P-wave,slow P-wave,SV-wave and SH-wave,and the anisotropy parameters can be therefore computed.Simulation results show that P-wave velocity have three attenuation peaks throughout the full frequency band,which respectively correspond to the influences of interlayer flow,the squirt flow and the Biot flow.Through the results of seismic velocity dispersion and attenuation at different incident angles,we find that the WIFF mechanism also has a significant impact on the dispersion characteristics of elastic anisotropy parameters within the low-mid frequency band.Moreover,it is shown that several poroelastic parameters,such as layer thickness ratio,crack aspect ratio and crack density have notable influence on seismic dispersion and attenuation.We compare the proposed modeled velocities with that given by the existing theory to confirm its validity.Our formulas and result can provide a better understanding of wave propagation in PLPC medium by considering the unified impacts of micro-,meso-and macro-scale WIFF mechanisms,which potentially lays a theoretical basis of rock physics for seismic interpretation.
基金support from the National Natural Science Foundation of China(U22A20171)the High Steel Center(HSC)at North China University of Technology and University of Science and Technology Beijing,China.
文摘A 1∶8 physical water model was constructed to investigate the fluid flow and mixing phenomena in the basic oxygen furnace(BOF)converter.The particle image velocimetry was employed to measure the velocity distribution of the bath and the high-speed camera was applied to capture the cavity shape in the combined blowing BOF converter.The mixing time for varied operating conditions was measured by the stimulus-response approach.The cavity depth increased with the decrease in the lance height and the increase in the top gas flow rate while the bottom blowing gas had little influence on the cavity depth.The minimum cavity depth was obtained under the condition of a 69.8 m^(3)/h top gas flow rate,a 287.5 mm lance height and a 0.93 m^(3)/h bottom blowing gas flow rate,which was 161.2 mm.The mixing time decreased as the lance height decreased and the top blowing gas flow rate increased.The mixing time was first decreased and then increased with the increase in the bottom gas flow rate.With the condition of 69.8 m^(3)/h gas flow rate of top blowing,the 287.5 mm lance height and the 0.93 m^(3)/h gas flow rate of bottom blowing,the mixing time in the converter was 48.65 s.The empirical formula between the stirring power and the mixing time in the converter was calculated.
文摘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 Innovative Research Groups Project of the Natural Science Foundation of Hebei Province(No.E2021203011)Central Government Guides Local Science and Technology Development Fund Projects(No.206Z1601G).
文摘The unstable fluid flow and severe free surface fluctuations in the wheel and belt caster can affect the quality of the cast bar.The lower level height tends to entrap inclusions in the molten metal.On the other hand,the higher level height makes the production process more dangerous due to the overflow of high temperature fluid from the mold.A computational model of the molten metal pouring process was established.The transient fluid flow and free surface fluctuations behavior were calculated using the three-dimensional large eddy simulation model and the volume of fluid model.The results show that the flow velocity of the main jet gradually decreases under the influence of the low kinetic energy fluid in the mold.There is an obvious oscillation in the tail of the jet,while the flow field is asymmetric in space.The jet is closer to the inside radius side due to the Coanda effect,and there is a recirculation zone on the inside radius and the outside radius respectively,according to the 10 s time-averaged results.Compared with the industrial observation and simulation results,the shape of the free surface is a wave that varies with time.In addition,the free surface height is lowest and the flow velocity is highest in the region near the jet.
文摘Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.
基金Projects (50934002, 51074013, 51104100) supported by the National Natural Science Foundation of China
文摘The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore structure of the media, was constructed. With this model, three dimensional pore scale fluid flow among particles was simulated. Then the distributions of fluid flow velocity and pressure were analyzed and the hydraulic conductivity was calculated. The simulation results indicate the fluid flow behaviors are mainly dominated by the volume and topological structure of pore space. There exist obvious preferential flow and leaching blind zones simultaneously in the medium. The highest velocities generally occur in those narrow pores with high pressure drops. The hydraulic conductivity obtained by simulation is the same order of magnitude as the laboratory test result, which denotes the validity of the model. The pore-scale and macro-scale are combined and the established geometrical model can be used for the simulations of other phenomena during heap leaching process.
基金The research is supported by China Postdoctoral Science Foundation （No. 20080430129 ） and National Key Technology R＆D Program （ No. 2007 BAE07 B07 ）.
文摘A two-equation K-ε turbulent fluid flow model is built to model the heat transfer and fluid flow in gas tungsten arc welding （GTAW） process of stainless steel S US310 and S US316. This model combines the buoyancy force, lorentz force and marangni force as the driving forces of thefluidflow in the weld pool. The material properties are functions of temperature in this model. The simulated results show that the molten metal flowing outward is mainly caused by the marangoni convection, which makes the weld pool become wider and shallower. The comparison of the weld pool shape of SUS310 and SUS316 shows that the slight differences of the value of thermal conductivity mainly attributes to the difference of the weld pool shape and the distinction of heat transport in laminar and turbulent model makes large diversity in the simulated results.
基金sponsored by the National Natural Science Foundation of China(Grant Nos.40904029 and 41274185)the Scientific Research Foundation for the Returned Overseas Chinese Scholars,State Education Ministry
文摘The bottom-simulating reflector（BSR） is a seismic indicator of the bottom of a gas hydrate stability zone. Its depth can be used to calculate the seafloor surface heat flow. The calculated BSR heat flow variations include disturbances from two important factors：（1） seafloor topography, which focuses the heat flow over regions of concave topography and defocuses it over regions of convex topography, and（2） the focused warm fluid flow within the accretionary prism coming from depths deeper than BSR. The focused fluid flow can be detected if the contribution of the topography to the BSR heat flow is removed. However, the analytical equation cannot solve the topographic effect at complex seafloor regions. We prove that 3D finite element method can model the topographic effect on the regional background heat flow with high accuracy, which can then be used to correct the topographic effect and obtain the BSR heat flow under the condition of perfectly flat topography. By comparing the corrected BSR heat flow with the regional background heat flow, focused fluid flow regions can be detected that are originally too small and cannot be detected using present-day equipment. This method was successfully applied to the midslope region of northern Cascadia subducting margin. The results suggest that the Cucumber Ridge and its neighboring area are positive heat flow anomalies, about 10%–20% higher than the background heat flow after 3D topographic correction. Moreover, the seismic imaging associated the positive heat flow anomaly areas with seabed fracture–cavity systems. This suggests flow of warm gas-carrying fluids along these high-permeability pathways, which could result in higher gas hydrate concentrations.
基金ChinaCommitteeofEducation theUniver sityofArizonaandtheMetropolitanWaterDistrictofSouthernCaliforni a.
文摘The continuum approach in fluid flow modeling is generally applied to porous geological media, but has limited applicability to fractured rocks. With the presence of a discrete fracture network relatively sparsely distributed in the matrix, it may be difficult or erroneous to use a porous medium fluid flow model with continuum assumptions to describe the fluid flow in fractured rocks at small or even large field scales. A discrete fracture fluid flow approach incorporating a stochastic fracture network with numerical fluid flow simulations could have the capability of capturing fluid flow behaviors such as inhomogeneity and anisotropy while reflecting the changes of hydraulic features at different scales. Moreover, this approach can be implemented to estimate the size of the representative elementary volume (REV) in order to find out the scales at which a porous medium flow model could be applied, and then to determine the hydraulic conductivity tensor for fractured rocks. The following topics are focused on in this study: (a) conceptual discrete fracture fluid flow modeling incorporating a stochastic fracture network with numerical flow simulations; (b) estimation of REV and hydraulic conductivity tensor for fractured rocks utilizing a stochastic fracture network with numerical fluid flow simulations; (c) investigation of the effect of fracture orientation and density on the hydraulic conductivity and REV by implementing a stochastic fracture network with numerical fluid flow simulations, and (d) fluid flow conceptual models accounting for major and minor fractures in the 2 D or 3 D flow fields incorporating a stochastic fracture network with numerical fluid flow simulations.
基金supported by National Natural Science Foundation of China(Grant No.41374116)the Fundamental Research Funds for Central Universities(Grant No.2014B39014)
文摘Local fluid flow（LFF） at the mesoscopic scale is the main dissipation mechanism of seismic waves in heterogeneous porous media within the seismic frequency band.LFF is easily influenced by the structure and boundary conditions of the porous media,which leads to different behaviors of the peak frequency of attenuation.The associated transition frequency can provide detailed information about the trend of LFF;therefore,research on the transition frequency of LFF and its relationship with the peak frequency of the corresponding attenuation（i.e.,inverse of quality factor） facilitates the detailed understanding of the effect of inner structures and boundary conditions in porous media.In this study,we firstly obtain the transition frequency of fluid flux based on Biot＇s theory of poroelasticity and the fast Fourier transform algorithm in a sample containing one repeating unit cell（RUC）.We then analyze changes of these two frequencies in porous media with different porous properties.Finally,we extend our analysis to the influence of the undrained boundary condition on the transition frequency and peak frequency in porous media with multiple RUCs.This setup can facilitate the understanding of the effect from the undrained boundary condition.Results demonstrate that these two frequencies have the same trend at low water saturation,but amplitude variations differ between the frequencies as the amount of saturation increases.However,for cases of high water saturation,both the trend and the amplitude variation of these two frequencies fit well with each other.
文摘This paper presents the results of a set of numerical models focussing on structural controls on hydrothermal mineralization. We first give an overview of natural phenomena of structurally-controlled ore formation and the background theory and mechanisms for such controls. We then provide the results of a group of simple 2D numerical models validated through comparison with Cu-vein structure observed near the Shilu Copper deposit （Yangchun, Guangdong Province, China） and finally a case study of 3D numerical modelling applied to the Hodgkinson Province in North Queensland （Australia）. Two modelling approaches, discrete deformation modelling and continuum coupled deformation and fluid flow modelling, are involved. The 2D model-derived patterns are remarkably consistent with the Cu-vein structure from the Shilu Copper deposit, and show that both modelling approaches can realistically simulate the mechanical behaviours of shear and dilatant fractures. The continuum coupled deformation and fluid flow model indicates that pattern of the Cu- veins near the Shilu deposit is the result of shear strain localization, development of dilation and fluid focussing into the dilatant fracture segments. The 3D case-study models （with deformation and fluid flow coupling） on the Hodgkinson Province generated a number of potential gold mineralization
基金supported by NSFC(41072069,40772061 and 40930423)State Basic Research Plan(2009CB421005)+1 种基金IRT (0755)111 Plan(B07011)
文摘The Ordos Basin of North China is not only an important uranium mineralization province, but also a major producer of oil, gas and coal in China. The genetic relationship between uranium mineralization and hydrocarbons has been recognized by a number of previous studies, but it has not been well understood in terms of the hydrodynamics of basin fluid flow. We have demonstrated in a previous study that the preferential localization of Cretaceous uranium mineralization in the upper part of the Ordos Jurassic section may have been related to the interface between an upward flowing, reducing fluid and a downward flowing, oxidizing fluid. This interface may have been controlled by the interplay between fluid overpressure related to disequilibrium sediment compaction and which drove the upward flow, and topographic relief, which drove the downward flow. In this study, we carried out numerical modeling for the contribution of oil and gas generation to the development of fluid overpressure, in addition to sedi- ment compaction and heating. Our results indicate that when hydrocarbon generation is taken into account, fluid overpressure during the Cretaceous was more than doubled in comparison with the simu- lation when hydrocarbon generation was not considered. Furthermore, fluid overpressure dissipation at the end of sedimentation slowed down relative to the no-hydrocarbon generation case. These results suggest that hydrocarbon generation may have played an important role in uranium mineralization, not only in providing reducing agents required for the mineralization, but also in contributing to the driving force to maintain the upward flow.
基金financially supported by the National Natural Science Foundation of China(grants No.41272159 and 41572099 and 41772272)China Geological Survey(grants No.1212011120964 and DD20180008)+2 种基金China Postdoctoral Science Foundation(Study on fluid diversity and genetic model of hydrothermal dolomitization,grant No.2016M601088)Liu Baojun Geoscience Foundation for Youthsthe Opening Foundation of Shandong Provincial Key Laboratory of Depositional Mineralization&Sedimentary Mineral(grant No.DMSM2017041)
文摘The Precambrian Dengying Formation is a set of large-scale, extensively dolomitized, carbonate reservoirs occurring within the Sichuan Basin. Petrographic and geochemical studies reveal dolomitization was a direct result of precipitation by chemically distinct fluids occurring at different times and at different intensities. Based on this evidence, dolomitization and multiple fluid flow events are analyzed, and three types of fluid evolution models are proposed. Results of analysis show that Precambrian Dengying Formation carbonates were deposited in a restricted peritidal environment （630- 542 Ma）. A high temperature and high Mg2＋ concentration seawater was a direct result of dolomitization for the micrite matrix, and for fibrous aragonite in primary pores. Geochemical evidence shows low δ18O values of micritic dolomite varying from -1.29%o to -4.52%o PDB, abundant light rare earth elements （REEs）, and low dolomite order degrees. Microbes and meteoric water significantly altered dolomite original chemical signatures, resulting in algal micritic dolomite and the fine-grained, granular, dolosparite dolomite having very negative δ18O values. Finely crystalline cement dolomite （536.3-280 Ma） and coarsely crystalline cement dolomite have a higher crystallization degree and higher order degree. The diagenetic sequence and fluid inclusion evidence imply a linear correlation between their burial depth and homogenization temperatures, which closely resemble the temperature of generated hydrocarbon. Compared with finely crystalline dolomite, precipitation of coarsely crystalline dolomite was more affected by restricted basinal fluids. In addition, there is a trend toward a more negative δ18O value, higher salinity, higher Fe and Mn concentrations, REE-rich. Two periods of hydrothermal fluids are identified, as the exceptionally high temperatures as opposed to the temperatures of burial history, in addition to the presence of high salinity fluid inclusions. The early hydrothermal fluid flow event was characterized by hot magnesium- and silicon-rich fluids, as demonstrated by the recrystallized matrix dolomite that is intimately associated with flint, opal, and microcrystalline quartz in intergranular or intercrystalline pores. This event was likely the result of a seafloor hydrothermal chimney eruption during Episode I of the Tongwan Movement （536.3-5.5 Ma）. In contrast, later hydrothermal fluids, which caused precipitation of saddle dolomite, were characterized by high salinity （15-16.05wt% NaCI equivalent） and homogenization temperatures （250 to 265℃）, δ18O values that were more enriched, and REE signatures. Geochemical data and the paragenetic sequence indicate that this hydrothermal fluid was related to extensive Permian large igneous province activity （360-280 Ma）. This study demonstrates the presence of complicated dolomitization processes occurring during various paleoclimates, tectonic cycles, and basinal fluids flow; results are a useful reference for these dolomitized Precambrian carbonates reservoirs.
基金supported by the Basic Research on Drilling & Completion of Critical Wells for Oil & Gas (Grant No. 51221003)National Science Fund for Petrochemical Industry (Project No. U1262201)+2 种基金"863" National Project (Project No. 2013AA064803)National Science Fund for Distinguished Young Scholars (Project No. 50925414)National Natural Science Foundation (Project No. 51074173)
文摘The special gas wettability phenomenon of reservoir rocks has been recognized by more and more researchers.It has a significant effect on efficient development of unconventional reservoirs.First,based on the preferentially gas-covered ability and surface free energy changes,definition and evaluation methods have been established.Second,a method for altering rock wettability and its mechanisms have been studied,surface oriented phenomena of functional groups with low surface energy are the fundamental reason for gas wettability alteration of rock.Third,the effect of gas wettability on the surface energy,electrical properties and dilatability are investigated.Last,the effects of gas wettability on capillary pressure,oil/gas/water distribution and flow are investigated with capillary tubes and etchedglass network models.The gas wettability theory of reservoir rocks has been initially established,which provides theoretical support for the efficient production of unconventional reservoirs and has great significance.
文摘Following up the fluid flow simulation in a 60 t tundish, the trajectories of inclusions in the 60 t tundish without flow control are simulated by considering the force balance between the drag force and the inertial buoyancy force. The Stochastic model yields more accurate inclusion motion than the non-Stochastic model due to including the effect of the turbulent fluctuation. The average residence time of inclusions decreases with increasing size. The thermal buoyancy favors inclusions removal especially the small inclusions. Using solute transport like the dye injection in water model and copper addition in the real steel tundish cannot accurately study the motion of the inclusions. In the simulation, more than 68% inclusions bigger than 10μm are removed to the top, and less than 32% enters the mold. The thermal buoyancy has little effect on the fraction of inclusions moved to the top of the inlet zone, and it mainly favors the removal of inclusions smaller than 100μm to the top surface of the outlet zone. For inclusions bigger than 100μm, the effect of thermal buoyancy on their motion can be ignored compared to the inertial buoyancy effect.
基金Supported by the National Natural Science Foundation of China (Nos.40930845 and 41006031)the International Science & Technology Cooperation Program of China (No. 2010DFA21740)the National Science and Technology Major Project (No. 2011ZX05026-004-06)
文摘The origin and migration of natural gas and the accumulation of gas hydrates within the Pearl River Mouth Basin of the northern South China Sea are poorly understood. Based on high-resolution 2D/3D seismic data, three environments of focused fluid flow： gas chimneys, mud diapirs and active faults have been identified. Widespread gas chimneys that act as important conduits for fluid flow are located below bottom simulating reflections and above basal uplifts. The occurrence and evolution of gas chimneys can be divided into a violent eruptive stage and a quiet seepage stage. For most gas chimneys, the strong eruptions are deduced to have happened during the Dongsha Movement in the latest Miocene, which are observed below Pliocene strata and few active faults develop above the top of the Miocene. The formation pressures of the Baiyun Sag currently are considered to be normal, based on these terms： 1） Borehole pressure tests with pressure coefficients of 1.043-l.047; 2） The distribution of gas chimneys is limited to strata older than the Pliocene; 3） Disseminated methane hydrates, rather than fractured hydrates, are found in the hydrate samples; 4） The gas hydrate is mainly charged with biogenic gas rather than thermogenic gas based on the chemical tests from gas hydrates cores. However, periods of quiet focused fluid flow also enable the establishment of good conduits for the migration of abundant biogenic gas and lesser volumes ofthermogenic gas. A geological model goveming fluid flow has been proposed to interpret the release of overpressure, the migration of fluids and the formation of gas hydrates, in an integrated manner. This model suggests that gas chimneys positioned above basal uplifts were caused by the Dongsha Movement at about 5.5 Ma. Biogenic gas occupies the strata above the base of the middle Miocene and migrates slowly into the gas chimney columns. Some of the biogenic gas and small volumes ofthermogenic gas eventually contribute to the formation of the gas hydrates.
基金the National Natural Science Foundation of China under Grant No.5181101756,51861165202 and No.51721092the Major Project of Science and Technology Innovation Special for Hubei Province under Grant No.2018AAA027+3 种基金the Fundamental Research Funds for the Central Universities,HUST:No.2018JYCXJJ034 and No.2019JYCXJJ025the Postdoctoral Science Foundation of China under Grant No.2018M632837the opening project of State Key Laboratory of Digital Manufacturing Equipment and Technology(HUST)under grant No.DMETKF2018001supported by the China Scholarship Council as a visiting scholar at the University of Virginia。
文摘Understanding the behaviors of heat transfer and fluid flow in weld pool and their effects on the solidification microstructure are significant for performance improvement of laser welds.This paper develops a three-dimensional numerical model to understand the multi-physical processes such as heat transfer,melt convection and solidification behavior in full-penetration laser welding of thin 5083 aluminum sheet.Solidification parameters including temperature gradient G and solidification rate R,and their combined forms are evaluated to interpret solidification microstructure.The predicted weld dimensions and the microstructure morphology and scale agree well with experiments.Results indicate that heat conduction is the dominant mechanism of heat transfer in weld pool,and melt convection plays a critical role in microstructure scale.The mushy zone shape/size and solidification parameters can be modulated by changing process parameters.Dendritic structures form because of the low G/R value.The scale of dendritic structures can be reduced by increasing GR via decreasing heat input.The columnar to equiaxed transition is predicted quantitatively via the process related G^3/R.These findings illustrate how heat transfer and fluid flow affect the solidification parameters and hence the microstructure,and show how to improve microstructure by optimizing the process.
