Detecting biomarkers in body fluids by optical lateral flow immune assay(LFIA) technology provides rapid access to disease information for early diagnosis.LFIA is based on an antigen-antibody reaction and is rapidly b...Detecting biomarkers in body fluids by optical lateral flow immune assay(LFIA) technology provides rapid access to disease information for early diagnosis.LFIA is based on an antigen-antibody reaction and is rapidly becoming the preferred choice of physicians and patients for point-of-care testing due to its simplicity,cost-effectiveness,and rapid detection.Observing the optical signal change from the colloidal gold of the traditional LFIA strip has been widely applied for various biomarkers detection in body fluids.Despite the significant progress,rapid real-time detection of color changes in the colloidal gold by the naked eye still faces many limitations,such as large errors and the inability to quantify and accurately detect.New optical LFIA strip technology has emerged in recent years to extend its application scenarios for achieving quantitative detection such as fluorescence,afterglow,and chemiluminescence.Herein,we summarized the development of optical LFIA technology from single to hyphenated optical signals for biomarkers detection in body fluids from invasive and non-invasive sources.Moreover,the challenge and outlook of optical LFIA strip technology are highlighted to inspire the designing of next-generation diagnostic platforms.展开更多
The airflow mechanics in adult nasal airways,whether healthy or abnormal,are extensively studied and investigated,but the flow mechanics in child nasal airways remain underexplored.This study investigates the airflow ...The airflow mechanics in adult nasal airways,whether healthy or abnormal,are extensively studied and investigated,but the flow mechanics in child nasal airways remain underexplored.This study investigates the airflow mechanics in the child’s nasal upper airway with adenoid hypertrophy,with an adenoid nasopharyngeal ratio(AN of 0.9),under cyclic inhalation and exhalation.An inlet respiratory cycle with three different flow rates(3.2 L/min calm breathing,8.6 L/min normal breathing,and 19.3 L/min intensive breathing)was simulated by using the computational fluid dynamics approach.To better capture the interaction between airflow and the flexible airway tissue,fluid-structure interaction analysis was performed at the normal breathing rate.Comparing the airflow dynamics during inhalation and exhalation,the pressure drops,nasal resistance,and wall shear stress show significant differences in the nasopharyngeal region for all different flow rates.This observation suggests that the inertial effect associated with the transient flow is important during exhalation and inhalation.Furthermore,the considerable temporal variation in flow rate distribution across a specific cross-section of the nasal airway highlights the critical role of transient data in virtual surgery planning and data for clinical decisions.展开更多
Particle-fluid two-phase flows in rock fractures and fracture networks play a pivotal role in determining the efficiency and effectiveness of hydraulic fracturing operations,a vital component in unconventional oil and...Particle-fluid two-phase flows in rock fractures and fracture networks play a pivotal role in determining the efficiency and effectiveness of hydraulic fracturing operations,a vital component in unconventional oil and gas extraction.Central to this phenomenon is the transport of proppants,tiny solid particles injected into the fractures to prevent them from closing once the injection is stopped.However,effective transport and deposition of proppant is critical in keeping fracture pathways open,especially in lowpermeability reservoirs.This review explores,then quantifies,the important role of fluid inertia and turbulent flows in governing proppant transport.While traditional models predominantly assume and then characterise flow as laminar,this may not accurately capture the complexities inherent in realworld hydraulic fracturing and proppant emplacement.Recent investigations highlight the paramount importance of fluid inertia,especially at the high Reynolds numbers typically associated with fracturing operations.Fluid inertia,often overlooked,introduces crucial forces that influence particle settling velocities,particle-particle interactions,and the eventual deposition of proppants within fractures.With their inherent eddies and transient and chaotic nature,turbulent flows introduce additional complexities to proppant transport,crucially altering proppant settling velocities and dispersion patterns.The following comprehensive survey of experimental,numerical,and analytical studies elucidates controls on the intricate dynamics of proppant transport under fluid inertia and turbulence-towards providing a holistic understanding of the current state-of-the-art,guiding future research directions,and optimising hydraulic fracturing practices.展开更多
This study presents a theoretical and experimental analysis of laminar flow behavior of high-viscous non-Newtonian power-law and Bingham fluids in multiple structural bifurcation channel distributors.Theoretical model...This study presents a theoretical and experimental analysis of laminar flow behavior of high-viscous non-Newtonian power-law and Bingham fluids in multiple structural bifurcation channel distributors.Theoretical models are derived to describe velocity profiles,shear rate distribution,and pressure drop across specific channels.Modified pressure drop models are proposed based on experiments on 6 kinds of non-Newtonian fluids and 3 types of bifurcation channel distributors with different transition areas.Specifically,the deviations between the theoretical models and the experimental results are systematically analyzed,and models were modified with correction coefficients based on Reynolds number and dimensionless shape factor.The theoretical results,modified results and experimental data are compared and discussed to ensure the accuracy of the modified models,demonstrating a significant improvement in the prediction of pressure drops for high-viscous non-Newtonian fluids in certain types of bifurcation channels.Additional experiments and analyses were carried out to validate the modified models and the results suggest that the models offer a certain degree of universal applicability in bifurcation channel designs.It implies that for other types of bifurcation channel distributors,similar correction methods based on Reynolds number and shape factor may also be applicable even with different coefficients.展开更多
Chinaʼs continental shale exhibits favorable geological characteristics and substantial resource potential,yet oil recovery for natural energy extraction remains critically low.Investigating the mechanisms of hydrauli...Chinaʼs continental shale exhibits favorable geological characteristics and substantial resource potential,yet oil recovery for natural energy extraction remains critically low.Investigating the mechanisms of hydraulically induced bedding fracture to generate complex fracture networks in continental shale,and establishing effective flow systems,is of utmost importance.This study employs laboratory experiments and numerical simulations to investigate the flow capacity and percolation behavior of hydraulically induced bedding fractures by different fluids in full-diameter shale cores.Hydraulic stimulation using different fluids generates bedding plane fracture networks,establishing effective flow systems.Eroded and detached shale fragments support localized fractures,thereby increasing their opening and enhancing flow capacity.Cetyltrimethylammonium bromide(CTAB)solution and SiO2 solution reduce the hydration of the shale surface,preventing shale fragments from swelling and disintegrating,leading to more stable percolation behavior.Eroded and spalled shale fragments near the injection point are transported to farther locations,where they help support localized fractures.This process differs from conventional hydraulic fracturing.Under a constant injection rate,the velocity in the smaller flow paths near the closure is significantly higher than that in the main flow paths,leading to pronounced bypass flow behavior.This restricts the percolation of fluid during imbibition in shale cores.The results provide valuable insights into the mechanism of hydraulically induced bedding fracture in continental shale,offering guidance for the effective development of shale reservoirs.展开更多
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.展开更多
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.展开更多
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.展开更多
This study rigorously examines the interplay between viscous dissipation,magnetic effects,and thermal radiation on the flow behavior of a non-Newtonian Carreau squeezed fluid passing by a sensor surface within a micro...This study rigorously examines the interplay between viscous dissipation,magnetic effects,and thermal radiation on the flow behavior of a non-Newtonian Carreau squeezed fluid passing by a sensor surface within a micro cantilever channel,aiming to deepen our understanding of heat transport processes in complex fluid dynamics scenarios.The primary objective is to elucidate how physical operational parameters influence both the velocity of fluid flow and its temperature distribution,utilizing a comprehensive numerical approach.Employing a combination of mathematical modeling techniques,including similarity transformation,this investigation transforms complex partial differential equations into more manageable ordinary ones,subsequently solving them using the homotopy perturbation method.By analyzing the obtained solutions and presenting them graphically,alongside detailed analysis,the study sheds light on the pivotal role of significant parameters in shaping fluid movement and energy distribution.Noteworthy observations reveal a substantial increase in fluid velocity with escalating magnetic parameters,while conversely,a contrasting trend emerges in the temperature distribution,highlighting the intricate relationship between magnetic effects,flow dynamics,and thermal behavior in non-Newtonian fluids.Further,the suction velocity enhance both the local skin friction and Nusselt numbers,whereas theWeissenberg number reduces them,opposite to the effect of the power-law index.展开更多
In this study,we proposed a numerical technique for solving time-dependent partial differential equations that arise in the electro-osmotic flowofCarreau fluid across a stationary plate based on amodified exponential ...In this study,we proposed a numerical technique for solving time-dependent partial differential equations that arise in the electro-osmotic flowofCarreau fluid across a stationary plate based on amodified exponential integrator.The scheme is comprised of two explicit stages.One is the exponential integrator type stage,and the second is the Runge-Kutta type stage.The spatial-dependent terms are discretized using the compact technique.The compact scheme can achieve fourth or sixth-order spatial accuracy,while the proposed scheme attains second-order temporal accuracy.Also,a mathematical model for the electro-osmotic flow of Carreau fluid over the stationary sheet is presented with heat and mass transfer effects.The governing equations are transformed into dimensionless partial differential equations and solved by the proposed scheme.Simulation results reveal that increasing the Helmholtz-Smoluchowski velocityUHS by 400%leads to a 60%-75%rise in peak flowvelocity,while the electro-osmotic parameter me enhances near-wall acceleration.Conversely,velocity decreases significantly with higher Weissenberg numbers,indicating the Carreau fluid’s elastic resistance and increased magnetic field strength due to improved Lorentz forces.Temperature rises with the thermal conductivity parameter 2,while higher reaction ratesγdiminish concentration and local Sherwood number values.The simulation findings show the scheme’s correctness and efficacy in capturing the complicated interactions in non-Newtonian electro-osmotic transport by revealing the notable impact of electrokinetic factors on flowbehaviour.Theproposedmodel is particularly relevant for BiologicalMicro-Electro-Mechanical Systems(BioMEMS)applications,where precise control of electro-thermal transport in non-Newtonian fluids is critical for lab-on-a-chip diagnostics,drug delivery,and micro-scale thermal management.展开更多
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 reservoir landslide is typically characterized by high-speed movement of a particle-fluid mixture,and its flow and deposit mechanisms are complex.This paper presents the mechanism of submerged granular column coll...The reservoir landslide is typically characterized by high-speed movement of a particle-fluid mixture,and its flow and deposit mechanisms are complex.This paper presents the mechanism of submerged granular column collapse under different densities ambient fluids based on coupled computational fluid dynamics and discrete element method(CFD-DEM)analysis.Important fluid-particle interaction forces,such as the drag force and the buoyancy,are considered by exchanging interaction forces between the CFD and DEM computations.We focus on the flow and deposit characteristics of submerged granular column collapse,namely the runout distance,the tail end height,the particle velocity,the energy,and deposit morphology,which are analyzed qualitatively and quantitatively.The change in fluid field caused by submerged granular column collapse and the formation of eddies are also discussed.A relatively dense fluid can significantly hinder the motion of granular flow,but can improve the conversion efficiency of kinetic energy from the vertical to the horizontal direction.Moreover,the eddies caused by fluid turbulence erode the surface of the granular pile,which is especially marked in a high-density fluid.The findings can provide vital theoretical support for the flow and deposit characteristics of granular flow under fluid and offer insights for the study of reservoir landslides.展开更多
Understanding the complex interaction between heat and mass transfer in non-Newtonian microflows is essential for the development and optimization of efficient microfluidic and thermal management systems.This study in...Understanding the complex interaction between heat and mass transfer in non-Newtonian microflows is essential for the development and optimization of efficient microfluidic and thermal management systems.This study investigates the magnetohydrodynamic(MHD)thermosolutal convection of a Casson fluid within an inclined,porous microchannel subjected to convective boundary conditions.The nonlinear,coupled equations governing momentum,energy,and species transport are solved numerically using the MATLAB bvp4c solver,ensuring high numerical accuracy and stability.To identify the dominant parameters influencing flow behavior and to optimize transport performance,a comprehensive hybrid optimization framework—combining a modified Taguchi design,Grey Relational Analysis(GRA),and Principal Component Analysis(PCA)—is proposed.This integrated strategy enables the simultaneous assessment of skin friction,Nusselt number,and Sherwood number,providing a rigorous multi-objective evaluation of system performance.Comparative validation with benchmark results from the literature confirms the accuracy and reliability of the present formulation and its numerical implementation.The results highlight the intricate coupling among flow slip,buoyancy effects,and convective transport mechanisms.Increased slip flow enhances axial velocity,while a higher solutal Biot number intensifies concentration gradients near the channel walls.Conversely,a lower thermal Biot number diminishes the temperature field,indicating weaker heat transfer across the boundaries.PCA results reveal that the first principal component(PC1)accounts for most of the system variance,demonstrating the dominant influence of coupled flow and transport parameters on overall system performance.展开更多
The dynamics of inviscid multi-component relativistic fluids may be modeled by the relativistic Euler equations, augmented by one (or more) additional species equation(s). We use the high-resolution staggered central ...The dynamics of inviscid multi-component relativistic fluids may be modeled by the relativistic Euler equations, augmented by one (or more) additional species equation(s). We use the high-resolution staggered central schemes to solve these equations. The equilibrium states for each component are coupled in space and time to have a common temperature and velocity. The current schemes can handle strong shocks and the oscillations near the interfaces are negligible, which usually happens in the multi-component flows. The schemes also guarantee the exact mass conservation for each component, the exact conservation of total momentum, and energy in the whole particle system. The central schemes are robust, reliable, compact and easy to implement. Several one- and two-dimensional numerical test cases are included in this paper, which validate the application of these schemes to relativistic multi-component flows.展开更多
The newly formulated non-Newtonian rivulet flows streaming down an inclined planar surface,with additional periodic perturbations arising from the application of the 2nd Stokes problem to the investigation of rivulet ...The newly formulated non-Newtonian rivulet flows streaming down an inclined planar surface,with additional periodic perturbations arising from the application of the 2nd Stokes problem to the investigation of rivulet dynamics,are demonstrated in the current research.Hereby,the 2nd Stokes problem assumes that the surface,with a thin shared layer of the fluid on it,oscillates in a harmonic manner along the x-axis of the rivulet flow,which coincides with the main flow direction streaming down the underlying surface.We obtain the exact extension of the rivulet flow family,clarifying the structure of the pressure field,which fully absorbs the arising perturbation.The profile of the velocity field is assumed to be Gaussian-type with a non-zero level of plasticity.Hence,the absolutely non-Newtonian case of the viscoplastic flow solution,which satisfies the motion and continuity equations,is considered(with particular cases of exact solutions for pressure).The perturbed governing equations of motion for rivulet flows then result in the Riccati-type ordinary differential equation(ODE),describing the dynamics of the coordinate x(t).The approximated schematic dynamics are presented in graphical plots.展开更多
In this paper,we study the asymptotic behavior of the micropolar fluid flow through a thin domain,assuming zero Dirichlet boundary condition on the top boundary,which is rapidly oscillating,and non-standard boundary c...In this paper,we study the asymptotic behavior of the micropolar fluid flow through a thin domain,assuming zero Dirichlet boundary condition on the top boundary,which is rapidly oscillating,and non-standard boundary conditions on the flat bottom.Assuming“Reynolds roughness regime”,in which the thickness of the domain is very small compared to the wavelength of the roughness(i.e.a very slight roughness),we rigorously derive a generalized Reynolds equation for pressure,clearly showing the roughness-induced effects.Moreover,we give expressions for the average velocity and microrotation.展开更多
An unsteady numerical simulation is conducted to examine the dynamic runback characteristics of a water film flow driven by a boundary layer airflow over a solid surface pertinent to the dynamic glaze ice accretion pr...An unsteady numerical simulation is conducted to examine the dynamic runback characteristics of a water film flow driven by a boundary layer airflow over a solid surface pertinent to the dynamic glaze ice accretion process over aircraft wing surfaces.The multiphase flow simulation results of the wind-driven water runback(WDWR)flow are compared quantitatively with the experimental results in terms of the time-dependent variations of the water film thickness profiles and evolution of the front contact point of the runback water film flow.The underlying mechanism of the intermittent water runback behavior is elucidated by analyzing the time evolution of the airflow velocity and vorticity fields above the runback water film flow over the solid surface.To the best knowledge of the authors,the work presented here is the first successful attempt to numerically examine the transient runback characteristics of WDWR flows.It serves as an excellent benchmark case for the development of best practices to model the important micro-physical processes responsible for the transient water transport over aircraft wing surfaces.展开更多
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.展开更多
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.展开更多
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.展开更多
基金supported by the National Natural Science Foundation of China (Nos.22234005,22494632,22404081)the Natural Science Foundation of Jiangsu Province (Nos.BK20222015,BK20240534)。
文摘Detecting biomarkers in body fluids by optical lateral flow immune assay(LFIA) technology provides rapid access to disease information for early diagnosis.LFIA is based on an antigen-antibody reaction and is rapidly becoming the preferred choice of physicians and patients for point-of-care testing due to its simplicity,cost-effectiveness,and rapid detection.Observing the optical signal change from the colloidal gold of the traditional LFIA strip has been widely applied for various biomarkers detection in body fluids.Despite the significant progress,rapid real-time detection of color changes in the colloidal gold by the naked eye still faces many limitations,such as large errors and the inability to quantify and accurately detect.New optical LFIA strip technology has emerged in recent years to extend its application scenarios for achieving quantitative detection such as fluorescence,afterglow,and chemiluminescence.Herein,we summarized the development of optical LFIA technology from single to hyphenated optical signals for biomarkers detection in body fluids from invasive and non-invasive sources.Moreover,the challenge and outlook of optical LFIA strip technology are highlighted to inspire the designing of next-generation diagnostic platforms.
基金supported by the National Key Research and Development Program of China(Grant No.2022YFF0707601).
文摘The airflow mechanics in adult nasal airways,whether healthy or abnormal,are extensively studied and investigated,but the flow mechanics in child nasal airways remain underexplored.This study investigates the airflow mechanics in the child’s nasal upper airway with adenoid hypertrophy,with an adenoid nasopharyngeal ratio(AN of 0.9),under cyclic inhalation and exhalation.An inlet respiratory cycle with three different flow rates(3.2 L/min calm breathing,8.6 L/min normal breathing,and 19.3 L/min intensive breathing)was simulated by using the computational fluid dynamics approach.To better capture the interaction between airflow and the flexible airway tissue,fluid-structure interaction analysis was performed at the normal breathing rate.Comparing the airflow dynamics during inhalation and exhalation,the pressure drops,nasal resistance,and wall shear stress show significant differences in the nasopharyngeal region for all different flow rates.This observation suggests that the inertial effect associated with the transient flow is important during exhalation and inhalation.Furthermore,the considerable temporal variation in flow rate distribution across a specific cross-section of the nasal airway highlights the critical role of transient data in virtual surgery planning and data for clinical decisions.
基金the Australian Research Council Discovery Project(ARC DP 220100851)scheme and would acknowledge that.
文摘Particle-fluid two-phase flows in rock fractures and fracture networks play a pivotal role in determining the efficiency and effectiveness of hydraulic fracturing operations,a vital component in unconventional oil and gas extraction.Central to this phenomenon is the transport of proppants,tiny solid particles injected into the fractures to prevent them from closing once the injection is stopped.However,effective transport and deposition of proppant is critical in keeping fracture pathways open,especially in lowpermeability reservoirs.This review explores,then quantifies,the important role of fluid inertia and turbulent flows in governing proppant transport.While traditional models predominantly assume and then characterise flow as laminar,this may not accurately capture the complexities inherent in realworld hydraulic fracturing and proppant emplacement.Recent investigations highlight the paramount importance of fluid inertia,especially at the high Reynolds numbers typically associated with fracturing operations.Fluid inertia,often overlooked,introduces crucial forces that influence particle settling velocities,particle-particle interactions,and the eventual deposition of proppants within fractures.With their inherent eddies and transient and chaotic nature,turbulent flows introduce additional complexities to proppant transport,crucially altering proppant settling velocities and dispersion patterns.The following comprehensive survey of experimental,numerical,and analytical studies elucidates controls on the intricate dynamics of proppant transport under fluid inertia and turbulence-towards providing a holistic understanding of the current state-of-the-art,guiding future research directions,and optimising hydraulic fracturing practices.
基金the financial support from the National Science and Technology Major Project of China(2024XXXXX2700)the National Natural Science Foundation of China (22408099)+1 种基金the Key Research and Development Program of Xinjiang Uygur Autonomous Region (2022B01032)the National Ten Thousand Talents Program。
文摘This study presents a theoretical and experimental analysis of laminar flow behavior of high-viscous non-Newtonian power-law and Bingham fluids in multiple structural bifurcation channel distributors.Theoretical models are derived to describe velocity profiles,shear rate distribution,and pressure drop across specific channels.Modified pressure drop models are proposed based on experiments on 6 kinds of non-Newtonian fluids and 3 types of bifurcation channel distributors with different transition areas.Specifically,the deviations between the theoretical models and the experimental results are systematically analyzed,and models were modified with correction coefficients based on Reynolds number and dimensionless shape factor.The theoretical results,modified results and experimental data are compared and discussed to ensure the accuracy of the modified models,demonstrating a significant improvement in the prediction of pressure drops for high-viscous non-Newtonian fluids in certain types of bifurcation channels.Additional experiments and analyses were carried out to validate the modified models and the results suggest that the models offer a certain degree of universal applicability in bifurcation channel designs.It implies that for other types of bifurcation channel distributors,similar correction methods based on Reynolds number and shape factor may also be applicable even with different coefficients.
基金supported by the Frontier and Fundamental Research of Active Nanofluids Flooding for Enhanced Oil Recovery through Discontinuous and Variable-circle Modes in High Temperature and High Salinity Offshore Oilfields(U22B6005)National Natural Science Foundation of China(No.52274037)Study on Key Issues of Enhanced Oil Recovery of Gulong Shale Oil(DQYT-2022-JS-761).
文摘Chinaʼs continental shale exhibits favorable geological characteristics and substantial resource potential,yet oil recovery for natural energy extraction remains critically low.Investigating the mechanisms of hydraulically induced bedding fracture to generate complex fracture networks in continental shale,and establishing effective flow systems,is of utmost importance.This study employs laboratory experiments and numerical simulations to investigate the flow capacity and percolation behavior of hydraulically induced bedding fractures by different fluids in full-diameter shale cores.Hydraulic stimulation using different fluids generates bedding plane fracture networks,establishing effective flow systems.Eroded and detached shale fragments support localized fractures,thereby increasing their opening and enhancing flow capacity.Cetyltrimethylammonium bromide(CTAB)solution and SiO2 solution reduce the hydration of the shale surface,preventing shale fragments from swelling and disintegrating,leading to more stable percolation behavior.Eroded and spalled shale fragments near the injection point are transported to farther locations,where they help support localized fractures.This process differs from conventional hydraulic fracturing.Under a constant injection rate,the velocity in the smaller flow paths near the closure is significantly higher than that in the main flow paths,leading to pronounced bypass flow behavior.This restricts the percolation of fluid during imbibition in shale cores.The results provide valuable insights into the mechanism of hydraulically induced bedding fracture in continental shale,offering guidance for the effective development of shale reservoirs.
基金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.
基金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.
基金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.
文摘This study rigorously examines the interplay between viscous dissipation,magnetic effects,and thermal radiation on the flow behavior of a non-Newtonian Carreau squeezed fluid passing by a sensor surface within a micro cantilever channel,aiming to deepen our understanding of heat transport processes in complex fluid dynamics scenarios.The primary objective is to elucidate how physical operational parameters influence both the velocity of fluid flow and its temperature distribution,utilizing a comprehensive numerical approach.Employing a combination of mathematical modeling techniques,including similarity transformation,this investigation transforms complex partial differential equations into more manageable ordinary ones,subsequently solving them using the homotopy perturbation method.By analyzing the obtained solutions and presenting them graphically,alongside detailed analysis,the study sheds light on the pivotal role of significant parameters in shaping fluid movement and energy distribution.Noteworthy observations reveal a substantial increase in fluid velocity with escalating magnetic parameters,while conversely,a contrasting trend emerges in the temperature distribution,highlighting the intricate relationship between magnetic effects,flow dynamics,and thermal behavior in non-Newtonian fluids.Further,the suction velocity enhance both the local skin friction and Nusselt numbers,whereas theWeissenberg number reduces them,opposite to the effect of the power-law index.
基金supported and funded by the Deanship of Scientific Research at Imam Mohammad Ibn Saud Islamic University(IMSIU)(grant number IMSIU-DDRSP2503).
文摘In this study,we proposed a numerical technique for solving time-dependent partial differential equations that arise in the electro-osmotic flowofCarreau fluid across a stationary plate based on amodified exponential integrator.The scheme is comprised of two explicit stages.One is the exponential integrator type stage,and the second is the Runge-Kutta type stage.The spatial-dependent terms are discretized using the compact technique.The compact scheme can achieve fourth or sixth-order spatial accuracy,while the proposed scheme attains second-order temporal accuracy.Also,a mathematical model for the electro-osmotic flow of Carreau fluid over the stationary sheet is presented with heat and mass transfer effects.The governing equations are transformed into dimensionless partial differential equations and solved by the proposed scheme.Simulation results reveal that increasing the Helmholtz-Smoluchowski velocityUHS by 400%leads to a 60%-75%rise in peak flowvelocity,while the electro-osmotic parameter me enhances near-wall acceleration.Conversely,velocity decreases significantly with higher Weissenberg numbers,indicating the Carreau fluid’s elastic resistance and increased magnetic field strength due to improved Lorentz forces.Temperature rises with the thermal conductivity parameter 2,while higher reaction ratesγdiminish concentration and local Sherwood number values.The simulation findings show the scheme’s correctness and efficacy in capturing the complicated interactions in non-Newtonian electro-osmotic transport by revealing the notable impact of electrokinetic factors on flowbehaviour.Theproposedmodel is particularly relevant for BiologicalMicro-Electro-Mechanical Systems(BioMEMS)applications,where precise control of electro-thermal transport in non-Newtonian fluids is critical for lab-on-a-chip diagnostics,drug delivery,and micro-scale thermal management.
基金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.
基金supported by the National Natural Science Foundation of China(No.51825905).
文摘The reservoir landslide is typically characterized by high-speed movement of a particle-fluid mixture,and its flow and deposit mechanisms are complex.This paper presents the mechanism of submerged granular column collapse under different densities ambient fluids based on coupled computational fluid dynamics and discrete element method(CFD-DEM)analysis.Important fluid-particle interaction forces,such as the drag force and the buoyancy,are considered by exchanging interaction forces between the CFD and DEM computations.We focus on the flow and deposit characteristics of submerged granular column collapse,namely the runout distance,the tail end height,the particle velocity,the energy,and deposit morphology,which are analyzed qualitatively and quantitatively.The change in fluid field caused by submerged granular column collapse and the formation of eddies are also discussed.A relatively dense fluid can significantly hinder the motion of granular flow,but can improve the conversion efficiency of kinetic energy from the vertical to the horizontal direction.Moreover,the eddies caused by fluid turbulence erode the surface of the granular pile,which is especially marked in a high-density fluid.The findings can provide vital theoretical support for the flow and deposit characteristics of granular flow under fluid and offer insights for the study of reservoir landslides.
文摘Understanding the complex interaction between heat and mass transfer in non-Newtonian microflows is essential for the development and optimization of efficient microfluidic and thermal management systems.This study investigates the magnetohydrodynamic(MHD)thermosolutal convection of a Casson fluid within an inclined,porous microchannel subjected to convective boundary conditions.The nonlinear,coupled equations governing momentum,energy,and species transport are solved numerically using the MATLAB bvp4c solver,ensuring high numerical accuracy and stability.To identify the dominant parameters influencing flow behavior and to optimize transport performance,a comprehensive hybrid optimization framework—combining a modified Taguchi design,Grey Relational Analysis(GRA),and Principal Component Analysis(PCA)—is proposed.This integrated strategy enables the simultaneous assessment of skin friction,Nusselt number,and Sherwood number,providing a rigorous multi-objective evaluation of system performance.Comparative validation with benchmark results from the literature confirms the accuracy and reliability of the present formulation and its numerical implementation.The results highlight the intricate coupling among flow slip,buoyancy effects,and convective transport mechanisms.Increased slip flow enhances axial velocity,while a higher solutal Biot number intensifies concentration gradients near the channel walls.Conversely,a lower thermal Biot number diminishes the temperature field,indicating weaker heat transfer across the boundaries.PCA results reveal that the first principal component(PC1)accounts for most of the system variance,demonstrating the dominant influence of coupled flow and transport parameters on overall system performance.
文摘The dynamics of inviscid multi-component relativistic fluids may be modeled by the relativistic Euler equations, augmented by one (or more) additional species equation(s). We use the high-resolution staggered central schemes to solve these equations. The equilibrium states for each component are coupled in space and time to have a common temperature and velocity. The current schemes can handle strong shocks and the oscillations near the interfaces are negligible, which usually happens in the multi-component flows. The schemes also guarantee the exact mass conservation for each component, the exact conservation of total momentum, and energy in the whole particle system. The central schemes are robust, reliable, compact and easy to implement. Several one- and two-dimensional numerical test cases are included in this paper, which validate the application of these schemes to relativistic multi-component flows.
文摘The newly formulated non-Newtonian rivulet flows streaming down an inclined planar surface,with additional periodic perturbations arising from the application of the 2nd Stokes problem to the investigation of rivulet dynamics,are demonstrated in the current research.Hereby,the 2nd Stokes problem assumes that the surface,with a thin shared layer of the fluid on it,oscillates in a harmonic manner along the x-axis of the rivulet flow,which coincides with the main flow direction streaming down the underlying surface.We obtain the exact extension of the rivulet flow family,clarifying the structure of the pressure field,which fully absorbs the arising perturbation.The profile of the velocity field is assumed to be Gaussian-type with a non-zero level of plasticity.Hence,the absolutely non-Newtonian case of the viscoplastic flow solution,which satisfies the motion and continuity equations,is considered(with particular cases of exact solutions for pressure).The perturbed governing equations of motion for rivulet flows then result in the Riccati-type ordinary differential equation(ODE),describing the dynamics of the coordinate x(t).The approximated schematic dynamics are presented in graphical plots.
文摘In this paper,we study the asymptotic behavior of the micropolar fluid flow through a thin domain,assuming zero Dirichlet boundary condition on the top boundary,which is rapidly oscillating,and non-standard boundary conditions on the flat bottom.Assuming“Reynolds roughness regime”,in which the thickness of the domain is very small compared to the wavelength of the roughness(i.e.a very slight roughness),we rigorously derive a generalized Reynolds equation for pressure,clearly showing the roughness-induced effects.Moreover,we give expressions for the average velocity and microrotation.
基金supported by the National Science Foundation(NSF)of the USA(Grant Nos.TIP-2140489,CBET-2313310,and CBET-2415347).
文摘An unsteady numerical simulation is conducted to examine the dynamic runback characteristics of a water film flow driven by a boundary layer airflow over a solid surface pertinent to the dynamic glaze ice accretion process over aircraft wing surfaces.The multiphase flow simulation results of the wind-driven water runback(WDWR)flow are compared quantitatively with the experimental results in terms of the time-dependent variations of the water film thickness profiles and evolution of the front contact point of the runback water film flow.The underlying mechanism of the intermittent water runback behavior is elucidated by analyzing the time evolution of the airflow velocity and vorticity fields above the runback water film flow over the solid surface.To the best knowledge of the authors,the work presented here is the first successful attempt to numerically examine the transient runback characteristics of WDWR flows.It serves as an excellent benchmark case for the development of best practices to model the important micro-physical processes responsible for the transient water transport over aircraft wing surfaces.
基金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.
基金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.
基金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.
