Understanding the anisotropic creep behaviors of shale under direct shearing is a challenging issue.In this context,we conducted shear-creep and steady-creep tests on shale with five bedding orientations (i.e.0°,...Understanding the anisotropic creep behaviors of shale under direct shearing is a challenging issue.In this context,we conducted shear-creep and steady-creep tests on shale with five bedding orientations (i.e.0°,30°,45°,60°,and 90°),under multiple levels of direct shearing for the first time.The results show that the anisotropic creep of shale exhibits a significant stress-dependent behavior.Under a low shear stress,the creep compliance of shale increases linearly with the logarithm of time at all bedding orientations,and the increase depends on the bedding orientation and creep time.Under high shear stress conditions,the creep compliance of shale is minimal when the bedding orientation is 0°,and the steady-creep rate of shale increases significantly with increasing bedding orientations of 30°,45°,60°,and 90°.The stress-strain values corresponding to the inception of the accelerated creep stage show an increasing and then decreasing trend with the bedding orientation.A semilogarithmic model that could reflect the stress dependence of the steady-creep rate while considering the hardening and damage process is proposed.The model minimizes the deviation of the calculated steady-state creep rate from the observed value and reveals the behavior of the bedding orientation's influence on the steady-creep rate.The applicability of the five classical empirical creep models is quantitatively evaluated.It shows that the logarithmic model can well explain the experimental creep strain and creep rate,and it can accurately predict long-term shear creep deformation.Based on an improved logarithmic model,the variations in creep parameters with shear stress and bedding orientations are discussed.With abovementioned findings,a mathematical method for constructing an anisotropic shear creep model of shale is proposed,which can characterize the nonlinear dependence of the anisotropic shear creep behavior of shale on the bedding orientation.展开更多
The elliptic azimuthal anisotropy coefficient(v_(2))of the identified particles at midrapidity(|η|<0.8)was investigated in p-Pb collisions at√s_(NN)=5.02 TeV using a multi-phase transport model(AMPT).The calculat...The elliptic azimuthal anisotropy coefficient(v_(2))of the identified particles at midrapidity(|η|<0.8)was investigated in p-Pb collisions at√s_(NN)=5.02 TeV using a multi-phase transport model(AMPT).The calculations of differential v_(2)based on the advanced flow extraction method of light flavor hadrons(pions,kaons,protons,andΛ)in small collision systems were extended to a wider transverse momentum(p_(T))range of up to 8 GeV/c for the first time.The string-melting version of the AMPT model provides a good description of the measured p_(T)-differential v_(2)of the mesons but exhibits a slight deviation from the baryon v_(2).In addition,we observed the features of mass ordering at low p_(T)and the approximate number-of-constituentquark(NCQ)scaling at intermediate p_(T).Moreover,we demonstrate that hadronic rescattering does not have a significant impact on v_(2)in p-Pb collisions for different centrality selections,whereas partonic scattering dominates in generating the elliptic anisotropy of the final particles.This study provides further insight into the origin of collective-like behavior in small collision systems and has referential value for future measurements of azimuthal anisotropy.展开更多
Based on the tortuous capillary network model,the relationship between anisotropic permeability and rock normal strain,namely the anisotropic dynamic permeability model(ADPM),was derived and established.The model was ...Based on the tortuous capillary network model,the relationship between anisotropic permeability and rock normal strain,namely the anisotropic dynamic permeability model(ADPM),was derived and established.The model was verified using pore-scale flow simulation.The uniaxial strain process was calculated and the main factors affecting permeability changes in different directions in the deformation process were analyzed.In the process of uniaxial strain during the exploitation of layered oil and gas reservoirs,the effect of effective surface porosity on the permeability in all directions is consistent.With the decrease of effective surface porosity,the sensitivity of permeability to strain increases.The sensitivity of the permeability perpendicular to the direction of compression to the strain decreases with the increase of the tortuosity,while the sensitivity of the permeability in the direction of compression to the strain increases with the increase of the tortuosity.For layered reservoirs with the same initial tortuosity in all directions,the tortuosity plays a decisive role in the relative relationship between the variations of permeability in all directions during pressure drop.When the tortuosity is less than 1.6,the decrease rate of horizontal permeability is higher than that of vertical permeability,while the opposite is true when the tortuosity is greater than 1.6.This phenomenon cannot be represented by traditional dynamic permeability model.After the verification by experimental data of pore-scale simulation,the new model has high fitting accuracy and can effectively characterize the effects of deformation in different directions on the permeability in all directions.展开更多
Sand typically exhibits anisotropic internal structure which may significantly influence its mechanical behavior. The material point method (MPM) can eliminate mesh distortion and thus is suitable for investigating ge...Sand typically exhibits anisotropic internal structure which may significantly influence its mechanical behavior. The material point method (MPM) can eliminate mesh distortion and thus is suitable for investigating geotechnical problems with large deformation. In this study, an advanced anisotropic critical state theory (ACST)-based soil model is implemented in MPM to study the response of strip footing resting on anisotropic sand. The capability of the model is verified by simulating several element tests and strip footing tests with different soil densities and fabric bedding plane orientations. For the footing problem with a vertical load, as the fabric bedding plane orientation increases, the bearing capacity decreases and its corresponding settlement increases. The failure pattern becomes asymmetrical when the bedding plane orientation or the loading direction is inclined. A comparison between the simulation results predicted by the anisotropic and isotropic models is made, which demonstrates that neglecting the fabric anisotropy may lead to the overestimation of the bearing capacity.展开更多
Understanding the quantitative responses of anisotropic dynamic properties in organic-rich shale with different kerogen content(KC)is of great significance in hydrocarbon exploration and development.Conducting control...Understanding the quantitative responses of anisotropic dynamic properties in organic-rich shale with different kerogen content(KC)is of great significance in hydrocarbon exploration and development.Conducting controlled experiments with a single variable is challenging for natural shales due to their high variations in components,diagenesis conditions,or pore fluid.We employed the hot-pressing technique to construct 11 well-controlled artificial shale with varying KC.These artificial shale samples were successive machined into prismatic shape for ultrasonic measurements along different directions.Observations revealed bedding perpendicular P-wave velocities are more sensitive to the increasing KC than bedding paralleling velocities due to the preferential alignments of kerogen.All elastic stiffnesses except C_(13)are generally decreasing with the increasing KC,the variation of C_(1) and C_(33)on kerogen content are more sensitive than those of C_(44)and C_(66).Apparent dynamic mechanical parameters(v and E)were found to have linear correlation with the true ones from complete anisotropic equations independent of KC,which hold value towards the interpretation of well logs consistently across formations,Anisotropic mechanical parameters(ΔE and brittlenessΔB)tend to decrease with the reducing KC,withΔB showing great sensitivity to KC variations.In the range of low KC(<10%),the V_(P)/V_(S) ratio demonstrated a linearly negative correlation with KC,and the V_(P)/V_(S) ratio magnitude of less than 1.75may serve as a significant characterization for highly organic-rich(>10%)shale,compilation of data from natural organic rich-shales globally verified the similar systematic relationships that can be empirically used to predict the fraction of KC in shales.展开更多
The paper presents the results of comprehensive studies of filtration and capacitance properties of highly porous reservoir rocks of the aquifer of an underground gas storage facility.The geomechanical part of the res...The paper presents the results of comprehensive studies of filtration and capacitance properties of highly porous reservoir rocks of the aquifer of an underground gas storage facility.The geomechanical part of the research included studying the dependence of rock permeability on the stress-strain state in the vicinity of the wells,and physical modeling of the implementation of the method of increasing the permeability of the wellbore zone-the method of directional unloading of the reservoir.The digital part of the research included computed tomography(CT)-based computer analysis of the internal structure,pore space characteristics,and filtration properties before and after the tests.According to the results of physical modeling of deformation and filtration processes,it is found that the permeability of rocks before fracture depends on the stress-strain state insignificantly,and this influence is reversible.However,when downhole pressure reaches 7-8 MPa,macrocracks in the rock begin to grow,accompanied by irreversible permeability increase.Porosity,geodesic tortuosity and permeability values were obtained based on digital studies and numerical modeling.A weak degree of transversal anisotropy of the filtration properties of rocks was detected.Based on the analysis of pore size distribution,pressure field and flow velocities,high homogeneity and connectivity of the rock pore space is shown.The absence of pronounced changes in pore space characteristics and pore permeability after non-uniform triaxial loading rocks was shown.On the basis of geometrical analysis of pore space,the reasons for weak permeability anisotropy were identified.The filtration-capacitance properties obtained from the digital analysis showed very good agreement with the results of field and laboratory measurements.The physical modeling has confirmed the efficiency of application of the directional unloading method for the reservoir under study.The necessary parameters of its application were calculated:bottomhole geometry,stage of operation,stresses and pressure drawdown value.展开更多
Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell＇s equations for an anisotropic mo...Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell＇s equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic （EM） fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.展开更多
To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of ...To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.展开更多
This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in r...This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in real seismic data application. Our results reveal that there is an offset limitation for both PP- and PS-waves in elliptical anisotropy fitting, and that PS-waves show a wider applicable offset range and larger observable azimuthal anisotropy than PP-waves. The major axis of the elliptical fit to the amplitudes of the R-component is perpendicular to the fracture strike, which is opposite to that in PP-wave analysis. The azimuthal interval travel time of PS-waves shows a nearly elliptical distribution and the major axis of the fit ellipse is perpendicular to the fracture strike, which is same as that in PP-wave analysis. For data within the applicable offset range, the anisotropic magnitude obtained from amplitude and travel time attributes of PP- and PS-waves exhibits a dependence on fracture density, and the major to minor axis ratio of the fit ellipse may be used to infer the relative distribution of fracture densities.展开更多
Brittleness analysis becomes important when looking for sweet spots in tightoil sandstone reservoirs. Hence, appropriate indices are required as accurate brittleness evaluation criteria. We construct a seismic rock ph...Brittleness analysis becomes important when looking for sweet spots in tightoil sandstone reservoirs. Hence, appropriate indices are required as accurate brittleness evaluation criteria. We construct a seismic rock physics model for tight-oil sandstone reservoirs with vertical fractures. Because of the complexities in lithology and pore structure and the anisotropic characteristics of tight-oil sandstone reservoirs, the proposed model is based on the solid components, pore connectivity, pore type, and fractures to better describe the sandstone reservoir microstructure. Using the model, we analyze the brittleness sensitivity of the elastic parameters in an anisotropic medium and establish a new brittleness index. We show the applicability of the proposed brittleness index for tight-oil sandstone reservoirs by considering the brittleness sensitivity, the rock physics response characteristics, and cross-plots. Compared with conventional brittleness indexes, the new brittleness index has high brittleness sensitivity and it is the highest in oil-bearing brittle zones with relatively high porosity. The results also suggest that the new brittleness index is much more sensitive to elastic properties variations, and thus can presumably better predict the brittleness characteristics of sweet spots in tight-oil sandstone reservoirs.展开更多
The single crystal blade is one of the key technologies for improving the performance, durability and reliability of aero-engines and ground gas-turbine engines. However, the anisotropic mechanical properties of the s...The single crystal blade is one of the key technologies for improving the performance, durability and reliability of aero-engines and ground gas-turbine engines. However, the anisotropic mechanical properties of the single crystal material makes a great deal of difficulties on the development and the application of the single crystal blade, which is a challenge for the engineering application of the single crystal superalloy and the theoretic bases of the application. Some researches on the strength analysis and the life prediction of the anisotropic single crystal blade were carried out by the authors' research team. They are as follows. The crystallographic constitutive models for the plastic and the creep behaviors and the method of the rupture life prediction were established and verified. The tensile or the creep experiments for DD3 single crystal alloy with different orientations under different temperatures and different tensile rates or under different temperatures and different stress levels were carried out. The experimental data and the anisotropic properties at intermediate and high temperatures revealed by the experiments are significant for the application of the single crystal alloy. In addition, the experimental research for a kind of single crystal blade was also made. As the application of the researches the strength analysis and the life prediction were carried out for the single crystal blade of a certain aero-engine. In this part, the constitutive models and their applications are described, and the experimental research work will be described in part II.展开更多
This paper studies the critical behaviours and magnetic properties of three-dimensional bond and anisotropy dilution Blume-Capel model （BCM） in the presence of an applied field within the effective field theory. The...This paper studies the critical behaviours and magnetic properties of three-dimensional bond and anisotropy dilution Blume-Capel model （BCM） in the presence of an applied field within the effective field theory. The trajectory of tricritical point, reentrant transitions and degenerate patterns of anisotropy are obtained both for the bond and the anisotropy dilutions. The global phase diagrams demonstrate unusually reentrant phenomena. The temperature dependences of magnetization curves undergo remarkable spin glass behaviour at low temperatures, and transform from ferromagnetism to paramagnetism at high temperature in applied fields. Temperature dependence of magnetic susceptibility curve is in qualitative agreement with experimental result.展开更多
The seismic responses of the shale properties are critical for shale gas reservoir evaluation and production.It has been widely reported that the clay minerals have substantial influences on the seismic wave anisotrop...The seismic responses of the shale properties are critical for shale gas reservoir evaluation and production.It has been widely reported that the clay minerals have substantial influences on the seismic wave anisotropy and brittleness.Hence,knowing the seismic responses of the clay-rich shales and estimation of shale elastic properties are significant for the shale gas industry.A physical model containing two groups of shale blocks as the target formations is constructed in laboratory.The group S contains six shale blocks with different clay contents,and the group N contains six shale blocks with different porosity.The acquired 2D seismic data is used to analyze the seismic responses of two corresponding seismic lines.An anisotropic three-term inversion method is applied to one of the seismic inline to estimate the elastic properties the target shale blocks.The inversed attributes can be used to reveal the effects of shale clay contents.This study shows the substantial significance of using a physical model to observe the seismic responses of shale properties.The inversion results indicate that the anisotropic three-term inversion method could provide accurate results of elastic properties as well as the P-wave anisotropy parameter for shale formations.展开更多
Natural fracture data from one of the Carboniferous shale masses in the eastern Qaidam Basin were used to establish a stochastic model of a discrete fracture network and to perform discrete element simulation research...Natural fracture data from one of the Carboniferous shale masses in the eastern Qaidam Basin were used to establish a stochastic model of a discrete fracture network and to perform discrete element simulation research on the size efect and mechanical parameters of shale.Analytical solutions of fctitious joints in transversely isotropic media were derived,which made it possible for the proposed numerical model to simulate the bedding and natural fractures in shale masses.The results indicate that there are two main factors infuencing the representative elementary volume(REV)size of a shale mass.The frst and most decisive factor is the presence of natural fractures in the block itself.The second is the anisotropy ratio:the greater the anisotropy is,the larger the REV.The bedding angle has little infuence on the REV size,whereas it has a certain infuence on the mechanical parameters of the rock mass.When the bedding angle approaches the average orientation of the natural fractures,the mechanical parameters of the shale blocks decrease greatly.The REV representing the mechanical properties of the Carboniferous shale masses in the eastern Qaidam Basin were comprehensively identifed by considering the infuence of bedding and natural fractures.When the numerical model size is larger than the REV,the fractured rock mass discontinuities can be transformed into equivalent continuities,which provides a method for simulating shale with natural fractures and bedding to analyze the stability of a borehole wall in shale.展开更多
The construction of a shale rock physics model and the selection of an appropriate brittleness index （B/） are two significant steps that can influence the accuracy of brittleness prediction. On one hand, the existin...The construction of a shale rock physics model and the selection of an appropriate brittleness index （B/） are two significant steps that can influence the accuracy of brittleness prediction. On one hand, the existing models of kerogen-rich shale are controversial, so a reasonable rock physics model needs to be built. On the other hand, several types of equations already exist for predicting the BI whose feasibility needs to be carefully considered. This study constructed a kerogen-rich rock physics model by performing the self- consistent approximation and the differential effective medium theory to model intercoupled clay and kerogen mixtures. The feasibility of our model was confirmed by comparison with classical models, showing better accuracy. Templates were constructed based on our model to link physical properties and the BL Different equations for the BI had different sensitivities, making them suitable for different types of formations. Equations based on Young＇s Modulus were sensitive to variations in lithology, while those using Lame＇s Coefficients were sensitive to porosity and pore fluids. Physical information must be considered to improve brittleness prediction.展开更多
The Brazilian test is a widely used method for determining the tensile strength of rocks and for calibrating parameters in bondedparticle models（BPMs）. In previous studies, the Brazilian disc has typically been trim...The Brazilian test is a widely used method for determining the tensile strength of rocks and for calibrating parameters in bondedparticle models（BPMs）. In previous studies, the Brazilian disc has typically been trimmed from a compacted rectangular specimen. The present study shows that different tensile strength values are obtained depending on the compressive loading direction. Several measures are proposed to reduce the anisotropy of the disc. The results reveal that the anisotropy of the disc is significantly influenced by the compactibility of the specimen from which it is trimmed. A new method is proposed in which the Brazilian disc is directly generated with a particle boundary, effectively reducing the anisotropy. The stiffness（particle and bond） and strength（bond） of the boundary are set at less than and greater than those of the disc assembly, respectively,which significantly decreases the stress concentration at the boundary contacts and prevents breakage of the boundary particle bonds. This leads to a significant reduction in the anisotropy of the disc and the discreteness of the tensile strength. This method is more suitable for carrying out a realistic Brazilian test for homogeneous rock-like material in the BPM.展开更多
The velocity field in meandering compound channels with overhank flow is highly three dimensional. To date, its features have been investigated experimentally and little research has been undertaken to investigate the...The velocity field in meandering compound channels with overhank flow is highly three dimensional. To date, its features have been investigated experimentally and little research has been undertaken to investigate the feasibility of reproducing these velocity fields with computer models. If computer modeling were to prove successful in this context, it could become a useful prediction technique and research tool to enhance our understanding of natural river dynamics. A 3-D k-E turbulence hydrodynamic model in curvilinear coordinates is established to simulate the overhank flow. The bodyfitted coordinate is adopted in the horizontal plane, the part grid is adopted in the vertical direction, and the wall-function method is employed to simulate the bed resistance. The model is applied to the simulation of the meandering channel with straight flood plain banks, and the main velocities and secondary velocities for both the longitudinal and cross sections are presented. Comparison and analysis show that the results of simulation are fit to reflect the results of experiment. These results show the application value of the model to 3D overhank flow.展开更多
Series of testing on coarse grained soils were carried out with a true triaxial testing apparatus. The loads were applied from the major principal and minor principal directions, respectively, to simulate the construc...Series of testing on coarse grained soils were carried out with a true triaxial testing apparatus. The loads were applied from the major principal and minor principal directions, respectively, to simulate the construction and water impounding process of a rock fill dam. The stress and strain relationships induced by the different loading methods were investigated. A remarkable stress-induced anisotropy under complex stress state was observed. Contrary to popular assumptions in traditional numerical analysis and constitutive models, it was found that different elastic modulus and Poisson ratio exist in different principal directions in rock fill dams. From the testing results, an anisotropic constitutive model based on Duncan-Chang nonlinear model is presented to overcome the limitations of axi-symmetric assumptions in conventional triaxial experiments and constitutive models. Both models were then applied in FEM analysis of an under-construction earth core high rock soil filled dam with the focus on hydraulic fracturing. The study reveals the major biases that exist when numerical analysis and constitutive models do not give serious consideration to the intermediate principal stress and anisotropy effects in soil rock built structures.展开更多
When there exists anisotropy in underground media elastic parameters of the observed coordinate possibly do not coincide with that of the natural coordinate. According to the theory that the density of potential energ...When there exists anisotropy in underground media elastic parameters of the observed coordinate possibly do not coincide with that of the natural coordinate. According to the theory that the density of potential energy, dissipating energy is independent of the coordinate, the relationship of elastic parameters between two coordinates is derived for two-phase anisotropic media. Then, pseudospectral method to solve wave equations of two-phase anisotropic media is derived. At last, we use this method to simulate wave propagation in two-phase anisotropic media four types of waves are observed in the snapshots, i.e., fast P wave and slow P wave, fast S wave and slow S wave. Shear wave splitting, SV wave cusps and elastic wave reflection and transmission are also observed.展开更多
A three-dimensional(3D)multiple phase field model,which takes into account the grain boundary(GB)energy anisotropy caused by texture,is established based on real grain orientations and Read-Shockley model.The model is...A three-dimensional(3D)multiple phase field model,which takes into account the grain boundary(GB)energy anisotropy caused by texture,is established based on real grain orientations and Read-Shockley model.The model is applied to the grain growth process of polycrystalline Mg(ZK60)alloy to investigate the evolution characteristics in different systems with varying proportions of low-angle grain boundary(LAGB)caused by different texture levels.It is found that the GB energy anisotropy can cause the grain growth kinetics to change,namely,higher texture levels(also means higher LAGB proportion)result in lower kinetics,and vice versa.The simulation results also show that the topological characteristics,such as LAGB proportion and distribution of grain size,undergo different evolution characteristics in different systems,and a more serious grain size fluctuation can be caused by a higher texture level.The mechanism is mainly the slower evolution of textured grains in their accumulation area and the faster coarsening rate of non-textured grains.Therefore,weakening the texture level is an effective way for implementing a desired homogenized microstructure in ZK60 Mg alloy.The rules revealed by the simulation results should be of great significance for revealing how the GB anisotropy affects the evolution of polycrystalline during the grain growth after recrystallization and offer the ideas for processing the alloy and optimizing the microstructure.展开更多
基金funded by the National Natural Science Foundation of China(Grant Nos.U22A20166 and 12172230)the Guangdong Basic and Applied Basic Research Foundation(Grant No.2023A1515012654)+1 种基金funded by the National Natural Science Foundation of China(Grant Nos.U22A20166 and 12172230)the Guangdong Basic and Applied Basic Research Foundation(Grant No.2023A1515012654)。
文摘Understanding the anisotropic creep behaviors of shale under direct shearing is a challenging issue.In this context,we conducted shear-creep and steady-creep tests on shale with five bedding orientations (i.e.0°,30°,45°,60°,and 90°),under multiple levels of direct shearing for the first time.The results show that the anisotropic creep of shale exhibits a significant stress-dependent behavior.Under a low shear stress,the creep compliance of shale increases linearly with the logarithm of time at all bedding orientations,and the increase depends on the bedding orientation and creep time.Under high shear stress conditions,the creep compliance of shale is minimal when the bedding orientation is 0°,and the steady-creep rate of shale increases significantly with increasing bedding orientations of 30°,45°,60°,and 90°.The stress-strain values corresponding to the inception of the accelerated creep stage show an increasing and then decreasing trend with the bedding orientation.A semilogarithmic model that could reflect the stress dependence of the steady-creep rate while considering the hardening and damage process is proposed.The model minimizes the deviation of the calculated steady-state creep rate from the observed value and reveals the behavior of the bedding orientation's influence on the steady-creep rate.The applicability of the five classical empirical creep models is quantitatively evaluated.It shows that the logarithmic model can well explain the experimental creep strain and creep rate,and it can accurately predict long-term shear creep deformation.Based on an improved logarithmic model,the variations in creep parameters with shear stress and bedding orientations are discussed.With abovementioned findings,a mathematical method for constructing an anisotropic shear creep model of shale is proposed,which can characterize the nonlinear dependence of the anisotropic shear creep behavior of shale on the bedding orientation.
基金This work was supported by the Key Laboratory of Quark and Lepton Physics(MOE)in Central China Normal University(Nos.QLPL2022P01,QLPL202106)Natural Science Foundation of Hubei Provincial Education Department(No.Q20131603)+2 种基金National key research,development program of China(No.2018YFE0104700)National Natural Science Foundation of China(No.12175085)Fundamental research funds for the Central Universities(No.CCNU220N003).
文摘The elliptic azimuthal anisotropy coefficient(v_(2))of the identified particles at midrapidity(|η|<0.8)was investigated in p-Pb collisions at√s_(NN)=5.02 TeV using a multi-phase transport model(AMPT).The calculations of differential v_(2)based on the advanced flow extraction method of light flavor hadrons(pions,kaons,protons,andΛ)in small collision systems were extended to a wider transverse momentum(p_(T))range of up to 8 GeV/c for the first time.The string-melting version of the AMPT model provides a good description of the measured p_(T)-differential v_(2)of the mesons but exhibits a slight deviation from the baryon v_(2).In addition,we observed the features of mass ordering at low p_(T)and the approximate number-of-constituentquark(NCQ)scaling at intermediate p_(T).Moreover,we demonstrate that hadronic rescattering does not have a significant impact on v_(2)in p-Pb collisions for different centrality selections,whereas partonic scattering dominates in generating the elliptic anisotropy of the final particles.This study provides further insight into the origin of collective-like behavior in small collision systems and has referential value for future measurements of azimuthal anisotropy.
基金Supported by the National Natural Science Foundation of China(52274048)Beijing Natural Science Foundation Project of China(3222037)Shaanxi Provincial Technical Innovation Project of China(2023-YD-CGZH-02).
文摘Based on the tortuous capillary network model,the relationship between anisotropic permeability and rock normal strain,namely the anisotropic dynamic permeability model(ADPM),was derived and established.The model was verified using pore-scale flow simulation.The uniaxial strain process was calculated and the main factors affecting permeability changes in different directions in the deformation process were analyzed.In the process of uniaxial strain during the exploitation of layered oil and gas reservoirs,the effect of effective surface porosity on the permeability in all directions is consistent.With the decrease of effective surface porosity,the sensitivity of permeability to strain increases.The sensitivity of the permeability perpendicular to the direction of compression to the strain decreases with the increase of the tortuosity,while the sensitivity of the permeability in the direction of compression to the strain increases with the increase of the tortuosity.For layered reservoirs with the same initial tortuosity in all directions,the tortuosity plays a decisive role in the relative relationship between the variations of permeability in all directions during pressure drop.When the tortuosity is less than 1.6,the decrease rate of horizontal permeability is higher than that of vertical permeability,while the opposite is true when the tortuosity is greater than 1.6.This phenomenon cannot be represented by traditional dynamic permeability model.After the verification by experimental data of pore-scale simulation,the new model has high fitting accuracy and can effectively characterize the effects of deformation in different directions on the permeability in all directions.
基金supported by the National Natural Science Foundation of China(Grant No.52108359).
文摘Sand typically exhibits anisotropic internal structure which may significantly influence its mechanical behavior. The material point method (MPM) can eliminate mesh distortion and thus is suitable for investigating geotechnical problems with large deformation. In this study, an advanced anisotropic critical state theory (ACST)-based soil model is implemented in MPM to study the response of strip footing resting on anisotropic sand. The capability of the model is verified by simulating several element tests and strip footing tests with different soil densities and fabric bedding plane orientations. For the footing problem with a vertical load, as the fabric bedding plane orientation increases, the bearing capacity decreases and its corresponding settlement increases. The failure pattern becomes asymmetrical when the bedding plane orientation or the loading direction is inclined. A comparison between the simulation results predicted by the anisotropic and isotropic models is made, which demonstrates that neglecting the fabric anisotropy may lead to the overestimation of the bearing capacity.
基金supported by the National Natural Science Foundation of China(42004112,42274175,42030812,41974160)Natural Science Foundation of Sichuan Province(2023NSFSC0764)。
文摘Understanding the quantitative responses of anisotropic dynamic properties in organic-rich shale with different kerogen content(KC)is of great significance in hydrocarbon exploration and development.Conducting controlled experiments with a single variable is challenging for natural shales due to their high variations in components,diagenesis conditions,or pore fluid.We employed the hot-pressing technique to construct 11 well-controlled artificial shale with varying KC.These artificial shale samples were successive machined into prismatic shape for ultrasonic measurements along different directions.Observations revealed bedding perpendicular P-wave velocities are more sensitive to the increasing KC than bedding paralleling velocities due to the preferential alignments of kerogen.All elastic stiffnesses except C_(13)are generally decreasing with the increasing KC,the variation of C_(1) and C_(33)on kerogen content are more sensitive than those of C_(44)and C_(66).Apparent dynamic mechanical parameters(v and E)were found to have linear correlation with the true ones from complete anisotropic equations independent of KC,which hold value towards the interpretation of well logs consistently across formations,Anisotropic mechanical parameters(ΔE and brittlenessΔB)tend to decrease with the reducing KC,withΔB showing great sensitivity to KC variations.In the range of low KC(<10%),the V_(P)/V_(S) ratio demonstrated a linearly negative correlation with KC,and the V_(P)/V_(S) ratio magnitude of less than 1.75may serve as a significant characterization for highly organic-rich(>10%)shale,compilation of data from natural organic rich-shales globally verified the similar systematic relationships that can be empirically used to predict the fraction of KC in shales.
基金supported by the Russian Science Foundation(Grant No.22-11-00273).
文摘The paper presents the results of comprehensive studies of filtration and capacitance properties of highly porous reservoir rocks of the aquifer of an underground gas storage facility.The geomechanical part of the research included studying the dependence of rock permeability on the stress-strain state in the vicinity of the wells,and physical modeling of the implementation of the method of increasing the permeability of the wellbore zone-the method of directional unloading of the reservoir.The digital part of the research included computed tomography(CT)-based computer analysis of the internal structure,pore space characteristics,and filtration properties before and after the tests.According to the results of physical modeling of deformation and filtration processes,it is found that the permeability of rocks before fracture depends on the stress-strain state insignificantly,and this influence is reversible.However,when downhole pressure reaches 7-8 MPa,macrocracks in the rock begin to grow,accompanied by irreversible permeability increase.Porosity,geodesic tortuosity and permeability values were obtained based on digital studies and numerical modeling.A weak degree of transversal anisotropy of the filtration properties of rocks was detected.Based on the analysis of pore size distribution,pressure field and flow velocities,high homogeneity and connectivity of the rock pore space is shown.The absence of pronounced changes in pore space characteristics and pore permeability after non-uniform triaxial loading rocks was shown.On the basis of geometrical analysis of pore space,the reasons for weak permeability anisotropy were identified.The filtration-capacitance properties obtained from the digital analysis showed very good agreement with the results of field and laboratory measurements.The physical modeling has confirmed the efficiency of application of the directional unloading method for the reservoir under study.The necessary parameters of its application were calculated:bottomhole geometry,stage of operation,stresses and pressure drawdown value.
基金financially supported by the National Hi-tech Research and Development Program of China(863 Program)(No.2012AA09A20103)
文摘Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell＇s equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic （EM） fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.
文摘To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.
文摘This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in real seismic data application. Our results reveal that there is an offset limitation for both PP- and PS-waves in elliptical anisotropy fitting, and that PS-waves show a wider applicable offset range and larger observable azimuthal anisotropy than PP-waves. The major axis of the elliptical fit to the amplitudes of the R-component is perpendicular to the fracture strike, which is opposite to that in PP-wave analysis. The azimuthal interval travel time of PS-waves shows a nearly elliptical distribution and the major axis of the fit ellipse is perpendicular to the fracture strike, which is same as that in PP-wave analysis. For data within the applicable offset range, the anisotropic magnitude obtained from amplitude and travel time attributes of PP- and PS-waves exhibits a dependence on fracture density, and the major to minor axis ratio of the fit ellipse may be used to infer the relative distribution of fracture densities.
基金supported by the National 973 project(Nos.2014CB239006 and 2011CB202402)the National Natural Science Foundation of China(Nos.41104069 and 41274124)+1 种基金Sinopec project(No.KJWX2014-05)the Fundamental Research Funds for the Central Universities(No.R1401005A)
文摘Brittleness analysis becomes important when looking for sweet spots in tightoil sandstone reservoirs. Hence, appropriate indices are required as accurate brittleness evaluation criteria. We construct a seismic rock physics model for tight-oil sandstone reservoirs with vertical fractures. Because of the complexities in lithology and pore structure and the anisotropic characteristics of tight-oil sandstone reservoirs, the proposed model is based on the solid components, pore connectivity, pore type, and fractures to better describe the sandstone reservoir microstructure. Using the model, we analyze the brittleness sensitivity of the elastic parameters in an anisotropic medium and establish a new brittleness index. We show the applicability of the proposed brittleness index for tight-oil sandstone reservoirs by considering the brittleness sensitivity, the rock physics response characteristics, and cross-plots. Compared with conventional brittleness indexes, the new brittleness index has high brittleness sensitivity and it is the highest in oil-bearing brittle zones with relatively high porosity. The results also suggest that the new brittleness index is much more sensitive to elastic properties variations, and thus can presumably better predict the brittleness characteristics of sweet spots in tight-oil sandstone reservoirs.
文摘The single crystal blade is one of the key technologies for improving the performance, durability and reliability of aero-engines and ground gas-turbine engines. However, the anisotropic mechanical properties of the single crystal material makes a great deal of difficulties on the development and the application of the single crystal blade, which is a challenge for the engineering application of the single crystal superalloy and the theoretic bases of the application. Some researches on the strength analysis and the life prediction of the anisotropic single crystal blade were carried out by the authors' research team. They are as follows. The crystallographic constitutive models for the plastic and the creep behaviors and the method of the rupture life prediction were established and verified. The tensile or the creep experiments for DD3 single crystal alloy with different orientations under different temperatures and different tensile rates or under different temperatures and different stress levels were carried out. The experimental data and the anisotropic properties at intermediate and high temperatures revealed by the experiments are significant for the application of the single crystal alloy. In addition, the experimental research for a kind of single crystal blade was also made. As the application of the researches the strength analysis and the life prediction were carried out for the single crystal blade of a certain aero-engine. In this part, the constitutive models and their applications are described, and the experimental research work will be described in part II.
基金Project supported by the Education Bureau Key Natural Science Foundation of Jiangsu Province （Grant No 03KJA140117） and by the 0pen Foundation of Jiangsu Key Laboratory of Thin Film （Grant No K2022）.
文摘This paper studies the critical behaviours and magnetic properties of three-dimensional bond and anisotropy dilution Blume-Capel model （BCM） in the presence of an applied field within the effective field theory. The trajectory of tricritical point, reentrant transitions and degenerate patterns of anisotropy are obtained both for the bond and the anisotropy dilutions. The global phase diagrams demonstrate unusually reentrant phenomena. The temperature dependences of magnetization curves undergo remarkable spin glass behaviour at low temperatures, and transform from ferromagnetism to paramagnetism at high temperature in applied fields. Temperature dependence of magnetic susceptibility curve is in qualitative agreement with experimental result.
基金supported by the National Natural Science Fund Projects(U19B6003)Strategic Cooperation Technology Projects of CNPC and CUPB(ZLZX2020-03)the Science Foundation of China University of Petroleum(Beijing)(2462020YXZZ008)。
文摘The seismic responses of the shale properties are critical for shale gas reservoir evaluation and production.It has been widely reported that the clay minerals have substantial influences on the seismic wave anisotropy and brittleness.Hence,knowing the seismic responses of the clay-rich shales and estimation of shale elastic properties are significant for the shale gas industry.A physical model containing two groups of shale blocks as the target formations is constructed in laboratory.The group S contains six shale blocks with different clay contents,and the group N contains six shale blocks with different porosity.The acquired 2D seismic data is used to analyze the seismic responses of two corresponding seismic lines.An anisotropic three-term inversion method is applied to one of the seismic inline to estimate the elastic properties the target shale blocks.The inversed attributes can be used to reveal the effects of shale clay contents.This study shows the substantial significance of using a physical model to observe the seismic responses of shale properties.The inversion results indicate that the anisotropic three-term inversion method could provide accurate results of elastic properties as well as the P-wave anisotropy parameter for shale formations.
基金support of the National Natural Science Foundation of China(51604275)the Key Laboratory of Urban Under Ground Engineering of Ministry of Education(TUE2018-01)+1 种基金Yue Qi Young Scholar Project of China University of Mining&Technology,Beijingthe Fundamental Research Funds for the Central Universities(2016QL02).
文摘Natural fracture data from one of the Carboniferous shale masses in the eastern Qaidam Basin were used to establish a stochastic model of a discrete fracture network and to perform discrete element simulation research on the size efect and mechanical parameters of shale.Analytical solutions of fctitious joints in transversely isotropic media were derived,which made it possible for the proposed numerical model to simulate the bedding and natural fractures in shale masses.The results indicate that there are two main factors infuencing the representative elementary volume(REV)size of a shale mass.The frst and most decisive factor is the presence of natural fractures in the block itself.The second is the anisotropy ratio:the greater the anisotropy is,the larger the REV.The bedding angle has little infuence on the REV size,whereas it has a certain infuence on the mechanical parameters of the rock mass.When the bedding angle approaches the average orientation of the natural fractures,the mechanical parameters of the shale blocks decrease greatly.The REV representing the mechanical properties of the Carboniferous shale masses in the eastern Qaidam Basin were comprehensively identifed by considering the infuence of bedding and natural fractures.When the numerical model size is larger than the REV,the fractured rock mass discontinuities can be transformed into equivalent continuities,which provides a method for simulating shale with natural fractures and bedding to analyze the stability of a borehole wall in shale.
基金supported by the NSFC and Sinopec Joint Key Project(No.U1663207)National Science and Technology Major Project(No.2017ZX05049-002)National 973 Program(No.2014CB239104)
文摘The construction of a shale rock physics model and the selection of an appropriate brittleness index （B/） are two significant steps that can influence the accuracy of brittleness prediction. On one hand, the existing models of kerogen-rich shale are controversial, so a reasonable rock physics model needs to be built. On the other hand, several types of equations already exist for predicting the BI whose feasibility needs to be carefully considered. This study constructed a kerogen-rich rock physics model by performing the self- consistent approximation and the differential effective medium theory to model intercoupled clay and kerogen mixtures. The feasibility of our model was confirmed by comparison with classical models, showing better accuracy. Templates were constructed based on our model to link physical properties and the BL Different equations for the BI had different sensitivities, making them suitable for different types of formations. Equations based on Young＇s Modulus were sensitive to variations in lithology, while those using Lame＇s Coefficients were sensitive to porosity and pore fluids. Physical information must be considered to improve brittleness prediction.
基金Support provided by the National Basic Research Program of China (2015CB258500, 2015CB058102, 2014CB046904)
文摘The Brazilian test is a widely used method for determining the tensile strength of rocks and for calibrating parameters in bondedparticle models（BPMs）. In previous studies, the Brazilian disc has typically been trimmed from a compacted rectangular specimen. The present study shows that different tensile strength values are obtained depending on the compressive loading direction. Several measures are proposed to reduce the anisotropy of the disc. The results reveal that the anisotropy of the disc is significantly influenced by the compactibility of the specimen from which it is trimmed. A new method is proposed in which the Brazilian disc is directly generated with a particle boundary, effectively reducing the anisotropy. The stiffness（particle and bond） and strength（bond） of the boundary are set at less than and greater than those of the disc assembly, respectively,which significantly decreases the stress concentration at the boundary contacts and prevents breakage of the boundary particle bonds. This leads to a significant reduction in the anisotropy of the disc and the discreteness of the tensile strength. This method is more suitable for carrying out a realistic Brazilian test for homogeneous rock-like material in the BPM.
文摘The velocity field in meandering compound channels with overhank flow is highly three dimensional. To date, its features have been investigated experimentally and little research has been undertaken to investigate the feasibility of reproducing these velocity fields with computer models. If computer modeling were to prove successful in this context, it could become a useful prediction technique and research tool to enhance our understanding of natural river dynamics. A 3-D k-E turbulence hydrodynamic model in curvilinear coordinates is established to simulate the overhank flow. The bodyfitted coordinate is adopted in the horizontal plane, the part grid is adopted in the vertical direction, and the wall-function method is employed to simulate the bed resistance. The model is applied to the simulation of the meandering channel with straight flood plain banks, and the main velocities and secondary velocities for both the longitudinal and cross sections are presented. Comparison and analysis show that the results of simulation are fit to reflect the results of experiment. These results show the application value of the model to 3D overhank flow.
基金Project(50809023)supported by the National Natural Science Foundation of ChinaProject(2015B17714)supported by the Fundamental Research Funds for Central Universities,China
文摘Series of testing on coarse grained soils were carried out with a true triaxial testing apparatus. The loads were applied from the major principal and minor principal directions, respectively, to simulate the construction and water impounding process of a rock fill dam. The stress and strain relationships induced by the different loading methods were investigated. A remarkable stress-induced anisotropy under complex stress state was observed. Contrary to popular assumptions in traditional numerical analysis and constitutive models, it was found that different elastic modulus and Poisson ratio exist in different principal directions in rock fill dams. From the testing results, an anisotropic constitutive model based on Duncan-Chang nonlinear model is presented to overcome the limitations of axi-symmetric assumptions in conventional triaxial experiments and constitutive models. Both models were then applied in FEM analysis of an under-construction earth core high rock soil filled dam with the focus on hydraulic fracturing. The study reveals the major biases that exist when numerical analysis and constitutive models do not give serious consideration to the intermediate principal stress and anisotropy effects in soil rock built structures.
文摘When there exists anisotropy in underground media elastic parameters of the observed coordinate possibly do not coincide with that of the natural coordinate. According to the theory that the density of potential energy, dissipating energy is independent of the coordinate, the relationship of elastic parameters between two coordinates is derived for two-phase anisotropic media. Then, pseudospectral method to solve wave equations of two-phase anisotropic media is derived. At last, we use this method to simulate wave propagation in two-phase anisotropic media four types of waves are observed in the snapshots, i.e., fast P wave and slow P wave, fast S wave and slow S wave. Shear wave splitting, SV wave cusps and elastic wave reflection and transmission are also observed.
基金Project supported by the National Key Research and Development Program of China(Grant No.2016YFB0701204)the 111 Project,China(Grant No.B20029)+2 种基金the Fundamental Research Funds for the Central Universities,China(Grant Nos.N2002017 and N2007011)the National Natural Science Foundation of China(Grant No.51571055)the Science Fund from the Science and Technology Bureau of Jiangyin High-Tech Industrial Development Zone,China(Grant No.ZX20200062)。
文摘A three-dimensional(3D)multiple phase field model,which takes into account the grain boundary(GB)energy anisotropy caused by texture,is established based on real grain orientations and Read-Shockley model.The model is applied to the grain growth process of polycrystalline Mg(ZK60)alloy to investigate the evolution characteristics in different systems with varying proportions of low-angle grain boundary(LAGB)caused by different texture levels.It is found that the GB energy anisotropy can cause the grain growth kinetics to change,namely,higher texture levels(also means higher LAGB proportion)result in lower kinetics,and vice versa.The simulation results also show that the topological characteristics,such as LAGB proportion and distribution of grain size,undergo different evolution characteristics in different systems,and a more serious grain size fluctuation can be caused by a higher texture level.The mechanism is mainly the slower evolution of textured grains in their accumulation area and the faster coarsening rate of non-textured grains.Therefore,weakening the texture level is an effective way for implementing a desired homogenized microstructure in ZK60 Mg alloy.The rules revealed by the simulation results should be of great significance for revealing how the GB anisotropy affects the evolution of polycrystalline during the grain growth after recrystallization and offer the ideas for processing the alloy and optimizing the microstructure.
