The mesoscopic failure mechanism and the macro-mechanical characteristics of soil-rock mixture(S-RM) under external load are largely controlled by S-RM's meso-structural features. The objective of this work is to ...The mesoscopic failure mechanism and the macro-mechanical characteristics of soil-rock mixture(S-RM) under external load are largely controlled by S-RM's meso-structural features. The objective of this work is to improve the three-dimensional technology for the generation of the random meso-structural models of S-RM, for randomly generating irregular rock blocks in S-RM with different shapes, sizes, and distributions according to the characteristics of the rock blocks' size distribution. Based on the new improved technology, a software system named as R-SRM3 D for generation and visualization of S-RM is developed. Using R-SRM3 D, a three-dimensional meso-structural model of S-RM is generated and used to study the meso-mechanical behavior through a series of true-triaxial numerical tests. From the numerical tests, the following conclusions are obtained. The meso-stress field of S-RM is influenced by the distribution of the internal rock blocks, and the macro-mechanical characteristics of S-RM are anisotropic in 3D; the intermediate principal stress and the soil-rock interface properties have significant influence on the macro strength of S-RM.展开更多
A large number of breakwaters have been constructed along coasts to protect humans and infrastructures from tsunamis.There is a risk that foundation soils of these structures may liquefy,or partially liquefy during th...A large number of breakwaters have been constructed along coasts to protect humans and infrastructures from tsunamis.There is a risk that foundation soils of these structures may liquefy,or partially liquefy during the earthquake preceding a tsunami,which would greatly reduce the structures’capacity to resist the tsunami.It is necessary to consider not only the soil’s liquefaction behavior due to earthquake motions but also its post-liquefaction behavior because this behavior will affect the breakwater’s capacity to resist an incoming tsunami.In this study,numerical tests based on a sophisticated constitutive model and a soil-water coupled finite element method are used to predict the mechanical behavior of breakwaters and the surrounding soils.Two real breakwaters subjected to two different seismic excitations are examined through numerical simulation.The simulation results show that,earthquakes affect not only the immediate behavior of breakwaters and the surrounding soils but also their long-term settlements due to post-earthquake consolidation.A soil profile with thick clayey layers beneath liquefied soil is more vulnerable to tsunami than a soil profile with only sandy layers.Therefore,quantitatively evaluating the seismic behavior of breakwaters and surrounding soils is important for the design of breakwater structures to resist tsunamis.展开更多
Depending on the numerical test approach on a computer, the relationships among relevant parameters, eg branch number, node number, mesh number, computation accuracy, preliminary value of airflow rate, iteration numbe...Depending on the numerical test approach on a computer, the relationships among relevant parameters, eg branch number, node number, mesh number, computation accuracy, preliminary value of airflow rate, iteration number, computation time and convergence in a mine ventilation network analysis, were investigated based on 5 mine ventilation systems. The results show that a higher computation accuracy greatly influences the iteration number. When the accuracy reaches 10-6m3·s-1 for solving a complicated mine ventilation network, the running time is too long though a high-speed computer is used. The preliminary value of airflow rate in the range of 1100m3·s-1 has little effects the iteration number. The structure of network also has some effect on the iteration number.展开更多
The impact trend of reck-coal system was studied by the method of accumulating and releasing of deformation energy and interaction of rock-coal system. The system model of roof-coal-floor was established. Based on the...The impact trend of reck-coal system was studied by the method of accumulating and releasing of deformation energy and interaction of rock-coal system. The system model of roof-coal-floor was established. Based on the RFPA software, rock fracture process analysis system, the numerical test of deformation, fracture and energy transmission of nonlinear and nonhomogeneous rock-coal system, and the numerical test and evaluating method of impact trend of reck-coal system were achieved. When the same coal seam was in different roof and floor conditions, the fracture process of reck-coal system can be classified as gradual, sudden and delayed fracture three kinds, and their impact trend can be classified as void, intense and medium correspondingly. The rock-coal system＇s impact trend is evaluated by the system impact index p and burst expanding forms. The criteria μ are μ〈1.0, 1.0≤μ〈l .5 and μ≥1.5 when the impact trend is void, intense or medium, which are tested and verified by the No.2 and No.4 coal seams in Sun- cun mine.展开更多
A two-dimensional mesoscale model is used to investigate the effect of background wind field andcharacter of the underlying surface on the coastal frontogenesis along the east continental coast when cold air outbreak....A two-dimensional mesoscale model is used to investigate the effect of background wind field andcharacter of the underlying surface on the coastal frontogenesis along the east continental coast when cold air outbreak. The results of numerical experiments indicate that the coastal frontogenesis mainly appears near ground atnight. When the background wind is easterly, it is favorable to coastal frontogenesis in the lower layer of the atmosphere. When the background wind is westerly, it is unfavorable to coastal frontogenesis in the lower layer of theatmosphere. The vertical shear of the background wind has important influence on the coastal frontogenesis. Thecoastal frontogenesis at lower layer at night is strengthened with increasing temperature difference between sea andair. The effect of coastal sloping terrain on the coastal frontogenesis is little when the topography sloping angle islesser.The analysis of the physical mechanism of frontogenesis shows that the main factor influencing the coastalfrontogenesis is the deformation field of horizontal wind speed.展开更多
Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer...Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.展开更多
In order to better understand the mechanical properties of graded crushed rocks （GCRs） and to optimize the relevant design, a numerical test method based on the particle flow modeling technique PFC2D is developed fo...In order to better understand the mechanical properties of graded crushed rocks （GCRs） and to optimize the relevant design, a numerical test method based on the particle flow modeling technique PFC2D is developed for the California bearing ratio （CBR） test on GGRs. The effects of different testing conditions and micro-mechanical parameters used in the model on the CBR numerical results have been systematically studied. The reliability of the numerical technique is verified. The numerical results suggest that the influences of the loading rate and Poisson＇s ratio on the CBR numerical test results are not significant. As such, a loading rate of 1.0-3.0 mm/min, a piston diameter of 5 cm, a specimen height of 15 cm and a specimen diameter of 15 cm are adopted for the CBR numerical test. The numerical results reveal that the GBR values increase with the friction coefficient at the contact and shear modulus of the rocks, while the influence of Poisson＇s ratio on the GBR values is insignificant. The close agreement between the CBR numerical results and experimental results suggests that the numerical simulation of the CBR values is promising to help assess the mechanical properties of GGRs and to optimize the grading design. Be- sides, the numerical study can provide useful insights on the mesoscopic mechanism.展开更多
For the improved two-sided projected quasi-Newton algorithms, which were presented in PartI, we prove in this paper that they are locally one-step or two-step superlinearly convergent. Numerical tests are reported the...For the improved two-sided projected quasi-Newton algorithms, which were presented in PartI, we prove in this paper that they are locally one-step or two-step superlinearly convergent. Numerical tests are reported thereafter. Results by solving a set of typical problems selectedfrom literature have demonstrated the extreme importance of these modifications in making Nocedal& Overton's original methon practical. Furthermore, these results show that the improved algoritnmsare very competitive in comparison with some highly praised sequential quadratic programmingmethods.展开更多
In order to understand the dynamics of granular flow on an erodible base soil,in this paper,a series of material point method-based granular column collapse tests were conducted to investigate numerically the mobility...In order to understand the dynamics of granular flow on an erodible base soil,in this paper,a series of material point method-based granular column collapse tests were conducted to investigate numerically the mobility and dynamic erosion process of granular flow subjected to the complex settings,i.e.,the aspect ratio,granular mass,friction and dilatancy resistance,gravity and presence of water.A set of power scaling laws were proposed to describe the final deposit characteristics of granular flow by the relations of the normalized run-out distance and the normalized final height of granular flow against the aspect ratio,being greatly affected by the complex geological settings,e.g.,granular mass,the friction and dilatancy resistance of granular soil,and presence of water in granular flow.An index of the coefficient of friction of granular soil was defined as a ratio of the target coefficient of friction over the initial coefficient of friction to quantify the scaling extent of friction change(i.e.,friction strengthening or weakening).There is a characteristic aspect ratio of granular column corresponding to the maximum mobility of granular flow with the minimum index of the apparent coefficient of friction.The index of the repose coefficient of friction of granular flow decreased gradually with the increase in aspect ratio because higher potential energy of granular column at a larger aspect ratio causes a larger kinetic energy of granular soil to weaken the friction of granular soil as a kind of velocity-related friction weakening.An increase in granular mass reduces gradually the indexes of the apparent and repose coefficients of friction of granular soil to enhance the mobility of granular flow.The mobility of granular flow increases gradually with the decrease in friction angle or increase in dilatancy angle of granular soil.However,the increase of gravity accelerates granular flow but showing the same final deposit profile without any dependence on gravity.The mobility of granular flow increases gradually by lowering the indexes of the apparent and repose coefficients of friction of granular flow while changing the surroundings,in turn,the dry soil,submerged soil and saturated soil,implying a gradually increased excessive mobility of granular flow with the friction weakening of granular soil.Presence of water in granular flow may be a potential catalyzer to yield a long run-out granular flow,as revealed in comparison of water-absent and water-present granular flows.In addition,the dynamic erosion and entrainment of based soil induced by granular flow subjected to the complex geological settings,i.e.,the aspect ratio,granular mass,gravity,friction and dilatancy resistance,and presence of water,were comprehensively investigated as well.展开更多
In order to advance the understanding of the impact dynamics of granular flow in complex geological settings,this paper studied the impact dynamics of granular flow on rigid barriers with a number of Material Point Me...In order to advance the understanding of the impact dynamics of granular flow in complex geological settings,this paper studied the impact dynamics of granular flow on rigid barriers with a number of Material Point Method(MPM)numerical tests.The impact behavior of granular flow on a rigid barrier was characterized by the initial dynamic impact stage,dynamic surge impact stage,compression impact stage and static stage of granular flow,where the impact force of granular flow was comprised of the dynamic and static forces of granular flow.The impact behavior of granular flow on a rigid barrier was characterized by the states of the fast or slow impact behavior of granular flow.The angle of slope and aspect ratio of granular soil greatly affected the impact behavior of granular flow on a column rigid barrier,where a power model was proposed to quantify the residual(Fnr)-over-maximum(Fnmax)normal impact force ratio of granular flow Fnr⁄Fnmax incorporating the effects of the angle of slope and aspect ratio of granular soil.With the increase of the column rigid barrier up to the semi-infinite column rigid barrier,the impact dynamics of granular flow gradually increased up to a maximum by progressively transforming the overflow into the dynamic surge impact of the incoming flow on the rigid barrier to capture more granular soil of granular flow against the rigid barrier.Presence of water in granular flow,i.e.,a mixture of solid and liquid in granular flow,yielded a dynamic coupling contribution of the solid and liquid,being accompanied by the whole dynamic process of granular flow,on the impact behavior of granular flow on a rigid barrier,where the liquid-phase material of granular flow,i.e.,the water,was predominant to contribute on the normal impact force of granular flow in comparison with the solid-phase material of granular flow.In addition,other factors,e.g.,the shape of rigid barrier(i.e.,the column barrier,arch barrier and circle barrier),and the gravity(i.e.,in the gravitational environment of the Moon,Earth and Mars),greatly affected the impact behavior of granular flow on a rigid barrier as well.展开更多
In this work,we study the error estimates of the fully discrete Fourier pseudospectral numerical scheme for solving the nonlocal volume-conserved Allen-Cahn(AC)equation.The time marching method of the numerical scheme...In this work,we study the error estimates of the fully discrete Fourier pseudospectral numerical scheme for solving the nonlocal volume-conserved Allen-Cahn(AC)equation.The time marching method of the numerical scheme is based on the well-known Invariant Energy Quadratization(IEQ)method.We demonstrate that the proposed fully discrete numerical method is uniquely solvable,unconditionally energy stable,and obtain the optimal error estimate of the scheme for both space and time.Additionally,we conduct several numerical tests to verify the theoretical results.展开更多
Almost all sandstone reservoirs contain interlayers. The identification and characterization of these interlayers iscritical for minimizing the uncertainty associated with oilfield development and improving oil and ga...Almost all sandstone reservoirs contain interlayers. The identification and characterization of these interlayers iscritical for minimizing the uncertainty associated with oilfield development and improving oil and gas recovery.Identifying interlayers outside wells using identification methods based on logging data and machine learning isdifficult and seismic-based identification techniques are expensive. Herein, a numerical model based on seepageand well-testing theories is introduced to identify interlayers using transient pressure data. The proposed modelrelies on the open-source MATLAB Reservoir Simulation Toolbox. The effects of the interlayer thickness, position,and width on the pressure response are thoroughly investigated. A procedure for inverting interlayer parametersin the reservoir using the bottom-hole pressure is also proposed. This method uses only transient pressuredata during well testing and can effectively identify the interlayer distribution near the wellbore at an extremelylow cost. The reliability of the model is verified using effective oilfield examples.展开更多
The cretaceous gas reservoir in Kelasu Gas Field of the Tarim Basin is a rare ultra-deep and ultra-high pressure fractured tight sandstone gas reservoir where multi-scale discrete fractures of matrix,fracture and faul...The cretaceous gas reservoir in Kelasu Gas Field of the Tarim Basin is a rare ultra-deep and ultra-high pressure fractured tight sandstone gas reservoir where multi-scale discrete fractures of matrix,fracture and fault are developed,so its development cannot be conducted just based on static and dynamic reservoir description.In order to solve this problem,this paper establishes a numerical well test model of vertical wells based on matrix,fractures and faults(large fractures and small faults)by combining the random generation of natural fracture networks with the unstructured discrete fracture modeling method to break through the traditional continuous medium well test model.In addition,the model is solved by using the finite element method with mixed element,and the typical well test type curves under different random fracture networks are obtained.And the following research results are obtained.First,based on the observed data,the fracture network distribution modes of fractured tight sandstone gas reservoirs are classified into three categories.The influence of random generation of fracture networks on typical well test type curves is discussed.The results of discrete fracture well test model are compared with those of the traditional continuous medium well test model,and the applicable conditions of the traditional continuous medium well test model is determined.Second,there are great differences between the results of discrete fracture model and those of dual porosity medium model.1 The dual porosity medium model is a special case of the discrete fracture model,in which the fractures are evenly distributed within infinitely small spacing.Third,the characteristics of well test type curves under three fracture network distribution modes are discussed.The well test type curves that cannot be interpreted by the conventional dual/triple porosity continuous medium model are successfully interpreted by using the established well test interpretation model of random discrete fracture.The curve matching effect is ideal and the interpreted parameters are reasonable.In conclusion,the new model and the new method reveal the development mechanism of step-by-step production and coordinated gas supply between media of different scales,explain the development characteristics of large inter-well productivity difference and abnormal rapid inter-well pressure response,and provide a reference for the development of similar gas reservoirs.展开更多
In order to explore the thermal conductivity of the natural poly-mineral rock,numerical tests of rock models with randomly-distributed components were conducted and compared with each other.Elaborately designed Monte ...In order to explore the thermal conductivity of the natural poly-mineral rock,numerical tests of rock models with randomly-distributed components were conducted and compared with each other.Elaborately designed Monte Carlo method was adopted to ingratiate the requirement of the random characteristics of grain size and the grains＇spatial distribution.This requirement was fulfilled by clustering the randomly generated unstructured tetrahedral elements in full three dimensions.Natural rocks are consisted of randomly distributed crystal particles or intergranular minerals.Our primary results verify that the thermal conductivity of the rock is strongly sensitive to the ingredients＇ volume fraction and their spatial distribution.Furthermore,we proved that,in order to reduce the measurement error to a reasonable range,the numerical specimen must be large enough or include sufficient number of mineral particles.Our numerical test results are in accordance with a variety of empirical formulas which are currently employed in petrology.展开更多
From the mathematical principles, the generalized potential theory can be employed to create constitutive model of geomaterial directly. The similar Cam-clay model, which is created based on the generalized potential ...From the mathematical principles, the generalized potential theory can be employed to create constitutive model of geomaterial directly. The similar Cam-clay model, which is created based on the generalized potential theory, has less assumptions,clearer mathematical basis, and better computational accuracy. Theoretically, it is more scientific than the traditional Cam-clay models. The particle flow code PFC3 D was used to make numerical tests to verify the rationality and practicality of the similar Cam-clay model. The verification process was as follows: 1) creating the soil sample for numerical test in PFC3 D, and then simulating the conventional triaxial compression test, isotropic compression test, and isotropic unloading test by PFC3D; 2)determining the parameters of the similar Cam-clay model from the results of above tests; 3) predicting the sample's behavior in triaxial tests under different stress paths by the similar Cam-clay model, and comparing the predicting results with predictions by the Cam-clay model and the modified Cam-clay model. The analysis results show that the similar Cam-clay model has relatively high prediction accuracy, as well as good practical value.展开更多
If a traditional explicit numerical integration algorithm is used to solve motion equation in the finite element simulation of wave motion, the time-step used by numerical integration is the smallest time-step restric...If a traditional explicit numerical integration algorithm is used to solve motion equation in the finite element simulation of wave motion, the time-step used by numerical integration is the smallest time-step restricted by the stability criterion in computational region. However, the excessively small time-step is usually unnecessary for a large portion of computational region. In this paper, a varying time-step explicit numerical integration algorithm is introduced, and its basic idea is to use different time-step restricted by the stability criterion in different computational region. Finally, the feasibility of the algorithm and its effect on calculating precision are verified by numerical test.展开更多
The hot or cold processing would induce the change and the inhomogeneous of the material mechanical properties in the local processing region of the structure,and it is difficult to obtain the specific mechanical prop...The hot or cold processing would induce the change and the inhomogeneous of the material mechanical properties in the local processing region of the structure,and it is difficult to obtain the specific mechanical properties in these regions by using the traditional material tensile test.To accurately get actual material mechanical properties in the local region of structure,a micro-indentation test system incorporated by an electronic universal material test device has been established.An indenter displacement sensor and a group of special micro-indenter assemblies are estab-lished.A numerical indentation inversion analysis method by using ABAQUS software is also proposed in this study.Based on the above test system and analysis platform,an approach to obtaining material mechanical properties in the local region of structures is proposed and established.The ball indentation test is performed and combined with the energy method by using various changed mechanical properties of 316L austenitic stainless steel under differ-ent elongations.The investigated results indicate that the material mechanical properties and the micro-indentation morphological changes have evidently relevance.Compared with the tensile test results,the deviations of material mechanical parameters,such as hardness H,the hardening exponent n,the yield strength σy and others are within 5%obtained through the indentation test and the finite element analysis.It provides an effective and convenient method for obtaining the actual material mechanical properties in the local processing region of the structure.展开更多
Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed. Bottom hole pressure response curves and the pressur...Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed. Bottom hole pressure response curves and the pressure field are obtained by solving the model equations with the finite-element method. By analyzing bottom hole pressure curves and the fluid flow in the pressure field, seven flow stages can be recognized on the curves. An upscaling method is developed to compare with the dual-porosity model （DPM）. The comparisons results show that the DPM overestimates the inter-porosity coefficient ）, and the storage factor w. The analysis results show that fracture conductivity plays a leading role in the fluid flow. Matrix permeability influences the beginning time of flow from the matrix to fractures. Fractures density is another important parameter controlling the flow. The fracture linear flow is hidden under the large fracture density. The pressure propagation is slower in the direction of larger fracture density.展开更多
Mitchell's solution is commonly used to determine the required strength of vertically exposed cemented backfill in mines. Developed for drained backfill, Mitchell model assumed a zero friction angle for the backfi...Mitchell's solution is commonly used to determine the required strength of vertically exposed cemented backfill in mines. Developed for drained backfill, Mitchell model assumed a zero friction angle for the backfill. Physical model tests were performed. Good agreements were obtained between the required strengths predicted by the analytical solution and experimental results. However, it is well-known that zero friction angle can only be possible in terms of total stresses when geomaterials are submitted to unconsolidated and undrained conditions. A revisit to Mitchell's physical model tests reveals that both the laboratory tests performed for obtaining the shear strength parameters of the cemented backfill and the box stability tests were conducted under a condition close to undrained condition. This explains well the good agreement between Mitchell's solution and experimental results. Good agreements are equally obtained between Mitchell's experimental results and FLAC3 D numerical modeling of shortterm stability analyses of exposed cemented backfill.展开更多
According to the preliminary concept of the high-level radioactive waste （HLW） repository in China, a large-scale mock-up facility, named China-Mock-Up was constructed in the laboratory of Beijing Research Institute...According to the preliminary concept of the high-level radioactive waste （HLW） repository in China, a large-scale mock-up facility, named China-Mock-Up was constructed in the laboratory of Beijing Research Institute of Uranium Geology （BRIUG）. A heater, which simulates a container of radioactive waste, is placed inside the compacted Gaomiaozi （GMZ）-Na-bentonite blocks and pellets. Water inflow through the barrier from its outer surface is used to simulate the intake of groundwater. The numbers of water injection pipes, injection pressure and the insulation layer were determined based on the nu- merical modeling simulations. The current experimental data of the facility are herein analyzed. The experiment is intended to evaluate the thermo-hydro-mechano-chemical （THMC） processes occurring in the compacted bentonite-buffer during the early stage of HLW disposal and to provide a reliable database for numerical modeling and further investigation of engineered barrier system （EBS）, and the design of HLW repository.展开更多
基金Project(51109117)supported by the National Natural Science Foundation of ChinaProject(20111081125)supported by the Independent Research Plan of Tsinghua University,ChinaProject(2013-KY-4)supported by the State Key Laboratory of Hydroscience and Engineering Project,China
文摘The mesoscopic failure mechanism and the macro-mechanical characteristics of soil-rock mixture(S-RM) under external load are largely controlled by S-RM's meso-structural features. The objective of this work is to improve the three-dimensional technology for the generation of the random meso-structural models of S-RM, for randomly generating irregular rock blocks in S-RM with different shapes, sizes, and distributions according to the characteristics of the rock blocks' size distribution. Based on the new improved technology, a software system named as R-SRM3 D for generation and visualization of S-RM is developed. Using R-SRM3 D, a three-dimensional meso-structural model of S-RM is generated and used to study the meso-mechanical behavior through a series of true-triaxial numerical tests. From the numerical tests, the following conclusions are obtained. The meso-stress field of S-RM is influenced by the distribution of the internal rock blocks, and the macro-mechanical characteristics of S-RM are anisotropic in 3D; the intermediate principal stress and the soil-rock interface properties have significant influence on the macro strength of S-RM.
基金National Natural Science Foundation of China under Grant Nos.51678369,41627801 and 41372284Technical Innovation Foundation of Shenzhen under Grant No.JCYJ20170302143610976+2 种基金Doctoral Fund of Shandong Province under Grant No.ZR2017BEE071China Postdoctoral Science Foundation under Grant No.2017M612227the Special Project Fund of Taishan Scholars of Shandong Province under Grant No.2015-212
文摘A large number of breakwaters have been constructed along coasts to protect humans and infrastructures from tsunamis.There is a risk that foundation soils of these structures may liquefy,or partially liquefy during the earthquake preceding a tsunami,which would greatly reduce the structures’capacity to resist the tsunami.It is necessary to consider not only the soil’s liquefaction behavior due to earthquake motions but also its post-liquefaction behavior because this behavior will affect the breakwater’s capacity to resist an incoming tsunami.In this study,numerical tests based on a sophisticated constitutive model and a soil-water coupled finite element method are used to predict the mechanical behavior of breakwaters and the surrounding soils.Two real breakwaters subjected to two different seismic excitations are examined through numerical simulation.The simulation results show that,earthquakes affect not only the immediate behavior of breakwaters and the surrounding soils but also their long-term settlements due to post-earthquake consolidation.A soil profile with thick clayey layers beneath liquefied soil is more vulnerable to tsunami than a soil profile with only sandy layers.Therefore,quantitatively evaluating the seismic behavior of breakwaters and surrounding soils is important for the design of breakwater structures to resist tsunamis.
基金Project (50474050) supported by the National Natural Science Foundation of China
文摘Depending on the numerical test approach on a computer, the relationships among relevant parameters, eg branch number, node number, mesh number, computation accuracy, preliminary value of airflow rate, iteration number, computation time and convergence in a mine ventilation network analysis, were investigated based on 5 mine ventilation systems. The results show that a higher computation accuracy greatly influences the iteration number. When the accuracy reaches 10-6m3·s-1 for solving a complicated mine ventilation network, the running time is too long though a high-speed computer is used. The preliminary value of airflow rate in the range of 1100m3·s-1 has little effects the iteration number. The structure of network also has some effect on the iteration number.
基金the Ministry of Education Backbone Teachers Funded Projects
文摘The impact trend of reck-coal system was studied by the method of accumulating and releasing of deformation energy and interaction of rock-coal system. The system model of roof-coal-floor was established. Based on the RFPA software, rock fracture process analysis system, the numerical test of deformation, fracture and energy transmission of nonlinear and nonhomogeneous rock-coal system, and the numerical test and evaluating method of impact trend of reck-coal system were achieved. When the same coal seam was in different roof and floor conditions, the fracture process of reck-coal system can be classified as gradual, sudden and delayed fracture three kinds, and their impact trend can be classified as void, intense and medium correspondingly. The rock-coal system＇s impact trend is evaluated by the system impact index p and burst expanding forms. The criteria μ are μ〈1.0, 1.0≤μ〈l .5 and μ≥1.5 when the impact trend is void, intense or medium, which are tested and verified by the No.2 and No.4 coal seams in Sun- cun mine.
文摘A two-dimensional mesoscale model is used to investigate the effect of background wind field andcharacter of the underlying surface on the coastal frontogenesis along the east continental coast when cold air outbreak. The results of numerical experiments indicate that the coastal frontogenesis mainly appears near ground atnight. When the background wind is easterly, it is favorable to coastal frontogenesis in the lower layer of the atmosphere. When the background wind is westerly, it is unfavorable to coastal frontogenesis in the lower layer of theatmosphere. The vertical shear of the background wind has important influence on the coastal frontogenesis. Thecoastal frontogenesis at lower layer at night is strengthened with increasing temperature difference between sea andair. The effect of coastal sloping terrain on the coastal frontogenesis is little when the topography sloping angle islesser.The analysis of the physical mechanism of frontogenesis shows that the main factor influencing the coastalfrontogenesis is the deformation field of horizontal wind speed.
基金supported by the National Natural Science Foundation of China(52104049)the Young Elite Scientist Sponsorship Program by Beijing Association for Science and Technology(BYESS2023262)。
文摘Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.
基金supported by the Program for New Century Excellent Talents in University (NCET-08-0749)Fundamental Research Funds for the Central Universities (CHD2012JC054)
文摘In order to better understand the mechanical properties of graded crushed rocks （GCRs） and to optimize the relevant design, a numerical test method based on the particle flow modeling technique PFC2D is developed for the California bearing ratio （CBR） test on GGRs. The effects of different testing conditions and micro-mechanical parameters used in the model on the CBR numerical results have been systematically studied. The reliability of the numerical technique is verified. The numerical results suggest that the influences of the loading rate and Poisson＇s ratio on the CBR numerical test results are not significant. As such, a loading rate of 1.0-3.0 mm/min, a piston diameter of 5 cm, a specimen height of 15 cm and a specimen diameter of 15 cm are adopted for the CBR numerical test. The numerical results reveal that the GBR values increase with the friction coefficient at the contact and shear modulus of the rocks, while the influence of Poisson＇s ratio on the GBR values is insignificant. The close agreement between the CBR numerical results and experimental results suggests that the numerical simulation of the CBR values is promising to help assess the mechanical properties of GGRs and to optimize the grading design. Be- sides, the numerical study can provide useful insights on the mesoscopic mechanism.
基金This research was supported in part by the National Natural Science Foundation of china
文摘For the improved two-sided projected quasi-Newton algorithms, which were presented in PartI, we prove in this paper that they are locally one-step or two-step superlinearly convergent. Numerical tests are reported thereafter. Results by solving a set of typical problems selectedfrom literature have demonstrated the extreme importance of these modifications in making Nocedal& Overton's original methon practical. Furthermore, these results show that the improved algoritnmsare very competitive in comparison with some highly praised sequential quadratic programmingmethods.
基金This work was supported by the National Natural Science Foundation of China(Grant no.U22A20603)Sichuan Science and Technology Program-China(Grant No.2023ZYD0149)CAS"Light of West China"Program-China(Grant No.Fangwei Yu).In addition,a special acknowledgement should be expressed to a famous Chinese television drama:My Chief and My Regiment that accompanied me(Dr.Fangwei Yu)through the loneliness time of completing this study.
文摘In order to understand the dynamics of granular flow on an erodible base soil,in this paper,a series of material point method-based granular column collapse tests were conducted to investigate numerically the mobility and dynamic erosion process of granular flow subjected to the complex settings,i.e.,the aspect ratio,granular mass,friction and dilatancy resistance,gravity and presence of water.A set of power scaling laws were proposed to describe the final deposit characteristics of granular flow by the relations of the normalized run-out distance and the normalized final height of granular flow against the aspect ratio,being greatly affected by the complex geological settings,e.g.,granular mass,the friction and dilatancy resistance of granular soil,and presence of water in granular flow.An index of the coefficient of friction of granular soil was defined as a ratio of the target coefficient of friction over the initial coefficient of friction to quantify the scaling extent of friction change(i.e.,friction strengthening or weakening).There is a characteristic aspect ratio of granular column corresponding to the maximum mobility of granular flow with the minimum index of the apparent coefficient of friction.The index of the repose coefficient of friction of granular flow decreased gradually with the increase in aspect ratio because higher potential energy of granular column at a larger aspect ratio causes a larger kinetic energy of granular soil to weaken the friction of granular soil as a kind of velocity-related friction weakening.An increase in granular mass reduces gradually the indexes of the apparent and repose coefficients of friction of granular soil to enhance the mobility of granular flow.The mobility of granular flow increases gradually with the decrease in friction angle or increase in dilatancy angle of granular soil.However,the increase of gravity accelerates granular flow but showing the same final deposit profile without any dependence on gravity.The mobility of granular flow increases gradually by lowering the indexes of the apparent and repose coefficients of friction of granular flow while changing the surroundings,in turn,the dry soil,submerged soil and saturated soil,implying a gradually increased excessive mobility of granular flow with the friction weakening of granular soil.Presence of water in granular flow may be a potential catalyzer to yield a long run-out granular flow,as revealed in comparison of water-absent and water-present granular flows.In addition,the dynamic erosion and entrainment of based soil induced by granular flow subjected to the complex geological settings,i.e.,the aspect ratio,granular mass,gravity,friction and dilatancy resistance,and presence of water,were comprehensively investigated as well.
基金supported by the Sichuan Science and Technology Program - China (Grant no. 2023ZYD0149)National Natural Science Foundation of China (Grant no. U22A20603)CAS "Light of West China" Program - China (Grant No. Fangwei Yu)
文摘In order to advance the understanding of the impact dynamics of granular flow in complex geological settings,this paper studied the impact dynamics of granular flow on rigid barriers with a number of Material Point Method(MPM)numerical tests.The impact behavior of granular flow on a rigid barrier was characterized by the initial dynamic impact stage,dynamic surge impact stage,compression impact stage and static stage of granular flow,where the impact force of granular flow was comprised of the dynamic and static forces of granular flow.The impact behavior of granular flow on a rigid barrier was characterized by the states of the fast or slow impact behavior of granular flow.The angle of slope and aspect ratio of granular soil greatly affected the impact behavior of granular flow on a column rigid barrier,where a power model was proposed to quantify the residual(Fnr)-over-maximum(Fnmax)normal impact force ratio of granular flow Fnr⁄Fnmax incorporating the effects of the angle of slope and aspect ratio of granular soil.With the increase of the column rigid barrier up to the semi-infinite column rigid barrier,the impact dynamics of granular flow gradually increased up to a maximum by progressively transforming the overflow into the dynamic surge impact of the incoming flow on the rigid barrier to capture more granular soil of granular flow against the rigid barrier.Presence of water in granular flow,i.e.,a mixture of solid and liquid in granular flow,yielded a dynamic coupling contribution of the solid and liquid,being accompanied by the whole dynamic process of granular flow,on the impact behavior of granular flow on a rigid barrier,where the liquid-phase material of granular flow,i.e.,the water,was predominant to contribute on the normal impact force of granular flow in comparison with the solid-phase material of granular flow.In addition,other factors,e.g.,the shape of rigid barrier(i.e.,the column barrier,arch barrier and circle barrier),and the gravity(i.e.,in the gravitational environment of the Moon,Earth and Mars),greatly affected the impact behavior of granular flow on a rigid barrier as well.
基金Supported by the National Natural Science Foundation of China(Grant Nos.12261017,62062018)the Foundation of Science and Technology of Guizhou Province(Grant No.ZK[2022]031)the Scientific Research Foundation of Guizhou University of Finance and Economics(Grant Nos.2022KYYB08,2022ZCZX077)。
文摘In this work,we study the error estimates of the fully discrete Fourier pseudospectral numerical scheme for solving the nonlocal volume-conserved Allen-Cahn(AC)equation.The time marching method of the numerical scheme is based on the well-known Invariant Energy Quadratization(IEQ)method.We demonstrate that the proposed fully discrete numerical method is uniquely solvable,unconditionally energy stable,and obtain the optimal error estimate of the scheme for both space and time.Additionally,we conduct several numerical tests to verify the theoretical results.
文摘Almost all sandstone reservoirs contain interlayers. The identification and characterization of these interlayers iscritical for minimizing the uncertainty associated with oilfield development and improving oil and gas recovery.Identifying interlayers outside wells using identification methods based on logging data and machine learning isdifficult and seismic-based identification techniques are expensive. Herein, a numerical model based on seepageand well-testing theories is introduced to identify interlayers using transient pressure data. The proposed modelrelies on the open-source MATLAB Reservoir Simulation Toolbox. The effects of the interlayer thickness, position,and width on the pressure response are thoroughly investigated. A procedure for inverting interlayer parametersin the reservoir using the bottom-hole pressure is also proposed. This method uses only transient pressuredata during well testing and can effectively identify the interlayer distribution near the wellbore at an extremelylow cost. The reliability of the model is verified using effective oilfield examples.
基金supported by the Major Science and Technology Project of PetroChina Company Limited“Research and application of key technologies for development of deep and ultra-deep gas reservoirs in Kuqa Depression(No.2018E-1803).
文摘The cretaceous gas reservoir in Kelasu Gas Field of the Tarim Basin is a rare ultra-deep and ultra-high pressure fractured tight sandstone gas reservoir where multi-scale discrete fractures of matrix,fracture and fault are developed,so its development cannot be conducted just based on static and dynamic reservoir description.In order to solve this problem,this paper establishes a numerical well test model of vertical wells based on matrix,fractures and faults(large fractures and small faults)by combining the random generation of natural fracture networks with the unstructured discrete fracture modeling method to break through the traditional continuous medium well test model.In addition,the model is solved by using the finite element method with mixed element,and the typical well test type curves under different random fracture networks are obtained.And the following research results are obtained.First,based on the observed data,the fracture network distribution modes of fractured tight sandstone gas reservoirs are classified into three categories.The influence of random generation of fracture networks on typical well test type curves is discussed.The results of discrete fracture well test model are compared with those of the traditional continuous medium well test model,and the applicable conditions of the traditional continuous medium well test model is determined.Second,there are great differences between the results of discrete fracture model and those of dual porosity medium model.1 The dual porosity medium model is a special case of the discrete fracture model,in which the fractures are evenly distributed within infinitely small spacing.Third,the characteristics of well test type curves under three fracture network distribution modes are discussed.The well test type curves that cannot be interpreted by the conventional dual/triple porosity continuous medium model are successfully interpreted by using the established well test interpretation model of random discrete fracture.The curve matching effect is ideal and the interpreted parameters are reasonable.In conclusion,the new model and the new method reveal the development mechanism of step-by-step production and coordinated gas supply between media of different scales,explain the development characteristics of large inter-well productivity difference and abnormal rapid inter-well pressure response,and provide a reference for the development of similar gas reservoirs.
基金supported by the National Science and Technology Project (No. Si-noProbe-07)National Basic Research Program of China (No. 2008cb425701)+1 种基金National Natural Science Foundation of China (No. 40774049)Institute of Earthquake Science,China Earthquake Administration,who provided us with the support from the collaborative project named as "Pre-study of several issues of testing ground for numerical forecast earthquakes in North China"
文摘In order to explore the thermal conductivity of the natural poly-mineral rock,numerical tests of rock models with randomly-distributed components were conducted and compared with each other.Elaborately designed Monte Carlo method was adopted to ingratiate the requirement of the random characteristics of grain size and the grains＇spatial distribution.This requirement was fulfilled by clustering the randomly generated unstructured tetrahedral elements in full three dimensions.Natural rocks are consisted of randomly distributed crystal particles or intergranular minerals.Our primary results verify that the thermal conductivity of the rock is strongly sensitive to the ingredients＇ volume fraction and their spatial distribution.Furthermore,we proved that,in order to reduce the measurement error to a reasonable range,the numerical specimen must be large enough or include sufficient number of mineral particles.Our numerical test results are in accordance with a variety of empirical formulas which are currently employed in petrology.
基金Projects(51378131,51378403)supported by the National Natural Science Foundation of ChinaProject(2012210020203)supported by the Fundamental Research Funds for the Central Universities,China
文摘From the mathematical principles, the generalized potential theory can be employed to create constitutive model of geomaterial directly. The similar Cam-clay model, which is created based on the generalized potential theory, has less assumptions,clearer mathematical basis, and better computational accuracy. Theoretically, it is more scientific than the traditional Cam-clay models. The particle flow code PFC3 D was used to make numerical tests to verify the rationality and practicality of the similar Cam-clay model. The verification process was as follows: 1) creating the soil sample for numerical test in PFC3 D, and then simulating the conventional triaxial compression test, isotropic compression test, and isotropic unloading test by PFC3D; 2)determining the parameters of the similar Cam-clay model from the results of above tests; 3) predicting the sample's behavior in triaxial tests under different stress paths by the similar Cam-clay model, and comparing the predicting results with predictions by the Cam-clay model and the modified Cam-clay model. The analysis results show that the similar Cam-clay model has relatively high prediction accuracy, as well as good practical value.
基金National Natural Science Foundation of China (50178065), 973 Program (2002CB412706), National Social Com-monweal Research Foundation (2002DIB30076) and Joint Seismological Science Foundation (101066).
文摘If a traditional explicit numerical integration algorithm is used to solve motion equation in the finite element simulation of wave motion, the time-step used by numerical integration is the smallest time-step restricted by the stability criterion in computational region. However, the excessively small time-step is usually unnecessary for a large portion of computational region. In this paper, a varying time-step explicit numerical integration algorithm is introduced, and its basic idea is to use different time-step restricted by the stability criterion in different computational region. Finally, the feasibility of the algorithm and its effect on calculating precision are verified by numerical test.
基金Supported by National Natural Science Foundation of China(Grant No.52075434)Key R&D Projects in Shaanxi Province(Grant No.2021KW-36).
文摘The hot or cold processing would induce the change and the inhomogeneous of the material mechanical properties in the local processing region of the structure,and it is difficult to obtain the specific mechanical properties in these regions by using the traditional material tensile test.To accurately get actual material mechanical properties in the local region of structure,a micro-indentation test system incorporated by an electronic universal material test device has been established.An indenter displacement sensor and a group of special micro-indenter assemblies are estab-lished.A numerical indentation inversion analysis method by using ABAQUS software is also proposed in this study.Based on the above test system and analysis platform,an approach to obtaining material mechanical properties in the local region of structures is proposed and established.The ball indentation test is performed and combined with the energy method by using various changed mechanical properties of 316L austenitic stainless steel under differ-ent elongations.The investigated results indicate that the material mechanical properties and the micro-indentation morphological changes have evidently relevance.Compared with the tensile test results,the deviations of material mechanical parameters,such as hardness H,the hardening exponent n,the yield strength σy and others are within 5%obtained through the indentation test and the finite element analysis.It provides an effective and convenient method for obtaining the actual material mechanical properties in the local processing region of the structure.
基金Project supported by the National Natural Science Foundation of China(No.5140232)the National Science and Technology Major Project(No.2011ZX05038003)the China Postdoctoral Science Foundation(No.2014M561074)
文摘Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed. Bottom hole pressure response curves and the pressure field are obtained by solving the model equations with the finite-element method. By analyzing bottom hole pressure curves and the fluid flow in the pressure field, seven flow stages can be recognized on the curves. An upscaling method is developed to compare with the dual-porosity model （DPM）. The comparisons results show that the DPM overestimates the inter-porosity coefficient ）, and the storage factor w. The analysis results show that fracture conductivity plays a leading role in the fluid flow. Matrix permeability influences the beginning time of flow from the matrix to fractures. Fractures density is another important parameter controlling the flow. The fracture linear flow is hidden under the large fracture density. The pressure propagation is slower in the direction of larger fracture density.
基金financial support from China Scholarship Council(CSC)under the Grant CSC No.201406460041financial support from the Natural Sciences and Engineering Research Council of Canada(NSERC 402318)+4 种基金the Institut de Recherche Robert-Sauvéen Santéet en Sécuritédu Travail(IRSST 2013-0029)Fonds de Recherche du Québec-Nature et Technologies(FRQNT 2015-MI-191676)the industrial partners of Research Institute on Mines and Environment(RIME UQAT-Polytechnique)The financial support from the National Science and Technology Support Program of China(No.2013BAB02B02)the Scientific Research Fund of Beijing General Research Institute of Mining and Metallurgy of China(No.YJ201507)
文摘Mitchell's solution is commonly used to determine the required strength of vertically exposed cemented backfill in mines. Developed for drained backfill, Mitchell model assumed a zero friction angle for the backfill. Physical model tests were performed. Good agreements were obtained between the required strengths predicted by the analytical solution and experimental results. However, it is well-known that zero friction angle can only be possible in terms of total stresses when geomaterials are submitted to unconsolidated and undrained conditions. A revisit to Mitchell's physical model tests reveals that both the laboratory tests performed for obtaining the shear strength parameters of the cemented backfill and the box stability tests were conducted under a condition close to undrained condition. This explains well the good agreement between Mitchell's solution and experimental results. Good agreements are equally obtained between Mitchell's experimental results and FLAC3 D numerical modeling of shortterm stability analyses of exposed cemented backfill.
文摘According to the preliminary concept of the high-level radioactive waste （HLW） repository in China, a large-scale mock-up facility, named China-Mock-Up was constructed in the laboratory of Beijing Research Institute of Uranium Geology （BRIUG）. A heater, which simulates a container of radioactive waste, is placed inside the compacted Gaomiaozi （GMZ）-Na-bentonite blocks and pellets. Water inflow through the barrier from its outer surface is used to simulate the intake of groundwater. The numbers of water injection pipes, injection pressure and the insulation layer were determined based on the nu- merical modeling simulations. The current experimental data of the facility are herein analyzed. The experiment is intended to evaluate the thermo-hydro-mechano-chemical （THMC） processes occurring in the compacted bentonite-buffer during the early stage of HLW disposal and to provide a reliable database for numerical modeling and further investigation of engineered barrier system （EBS）, and the design of HLW repository.
