In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hy...In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.展开更多
In the process of using the original key stratum theory to predict the height of a water-flowing fractured zone(WFZ),the influence of rock strata outside the calculation range on the rock strata within the calculation...In the process of using the original key stratum theory to predict the height of a water-flowing fractured zone(WFZ),the influence of rock strata outside the calculation range on the rock strata within the calculation range as well as the fact that the shape of the overburden deformation area will change with the excavation length are ignored.In this paper,an improved key stratum theory(IKS theory)was proposed by fixing these two shortcomings.Then,a WFZ height prediction method based on IKS theory was established and applied.First,the range of overburden involved in the analysis was determined according to the tensile stress distribution range above the goaf.Second,the key stratum in the overburden involved in the analysis was identified through IKS theory.Finally,the tendency of the WFZ to develop upward was determined by judging whether or not the identified key stratum will break.The proposed method was applied and verified in a mining case study,and the reasons for the differences in the development patterns between the WFZs in coalfields in Northwest and East China were also fully explained by this method.展开更多
Ecoregion-based height-diameter models were developed in the present study for Scots pine(Pinus sylves-tris L.)stands in Turkiye and included several ecological factors derived from a pre-existing ecoregional classifi...Ecoregion-based height-diameter models were developed in the present study for Scots pine(Pinus sylves-tris L.)stands in Turkiye and included several ecological factors derived from a pre-existing ecoregional classification system.The data were obtained from 2831 sample trees in 292 sample plots.Ten generalized height–diameter models were developed,and the best model(HD10)was selected according to statistical criteria.Then,nonlinear mixed-effects modeling was applied to the best model.The R2 for the generalized height‒diameter model(Richards function)modified by Sharma and Parton is 0.951,and the final model included number of trees,dominant height,and diameter at breast height,with a random parameter associated with each ecoregion attached to the inverse of the mean basal area.The full model predictions using the nonlinear mixed-effects model and the reduced model(HD10)predictions were compared using the nonlinear sum of extra squares test,which revealed significant differences between ecore-gions;ecoregion-based height–diameter models were thus found to be suitable to use.In addition,using these models in appropriate ecoregions was very important for achieving reliable predictions with low prediction errors.展开更多
Non-pillar continuous mining(NPCM) is regarded as a high-efficient, high-level and one-step mining technology, which can be divided into two substopes. Back fill stability status in substope I, which directly influenc...Non-pillar continuous mining(NPCM) is regarded as a high-efficient, high-level and one-step mining technology, which can be divided into two substopes. Back fill stability status in substope I, which directly influence the loss rate and dilution rate, etc, will determine whether the experimental research is successful or not. By employing energy method of limit analysis and finite element numerical simulation method, the critical backfill height was determined under the prerequisite condition of its stability, which put forward theoretical basis for reasonable and correct selection of backfill’s parameters. The result showed that the first backfill could not keep stable for NPCM, while the other was able to.展开更多
To avoid impacts and vibrations during the processes of acceleration and deceleration while possessing flexible working ways for cable-suspended parallel robots(CSPRs),point-to-point trajectory planning demands an und...To avoid impacts and vibrations during the processes of acceleration and deceleration while possessing flexible working ways for cable-suspended parallel robots(CSPRs),point-to-point trajectory planning demands an under-constrained cable-suspended parallel robot(UCPR)with variable angle and height cable mast as described in this paper.The end-effector of the UCPR with three cables can achieve three translational degrees of freedom(DOFs).The inverse kinematic and dynamic modeling of the UCPR considering the angle and height of cable mast are completed.The motion trajectory of the end-effector comprising six segments is given.The connection points of the trajectory segments(except for point P3 in the X direction)are devised to have zero instantaneous velocities,which ensure that the acceleration has continuity and the planned acceleration curve achieves smooth transition.The trajectory is respectively planned using three algebraic methods,including fifth degree polynomial,cycloid trajectory,and double-S velocity curve.The results indicate that the trajectory planned by fifth degree polynomial method is much closer to the given trajectory of the end-effector.Numerical simulation and experiments are accomplished for the given trajectory based on fifth degree polynomial planning.At the points where the velocity suddenly changes,the length and tension variation curves of the planned and unplanned three cables are compared and analyzed.The OptiTrack motion capture system is adopted to track the end-effector of the UCPR during the experiment.The effectiveness and feasibility of fifth degree polynomial planning are validated.展开更多
This paper presents a refined method for estimating the annual extreme wave heights at a coastal or offshore project site on the basis of the data acquired at some nearby routine hydrographic stations. This method is ...This paper presents a refined method for estimating the annual extreme wave heights at a coastal or offshore project site on the basis of the data acquired at some nearby routine hydrographic stations. This method is based on the orthogonality principle in linear mean square estimation of stochastic processes. The error of the method is analyzed and compared with that of the conventional method. It is found that the method is able to effectively reduce the error so long as some feasible measures are adopted. A simulated test of the method has been conducted in a large scale wind wave flume. The test results are in good agreement with those given by theoretical error analysis. A scheme to implement the method is proposed on the basis of error analysis. The scheme is so designed as to reduce the estimation error as far as possible. This method is also suitable to utilizing satellite wave data for the estimation.展开更多
This paper presents a new method to estimate the height of the atmospheric boundary layer(ABL) by using COSMIC radio occultation bending angle(BA) data. Using the numerical differentiation method combined with the reg...This paper presents a new method to estimate the height of the atmospheric boundary layer(ABL) by using COSMIC radio occultation bending angle(BA) data. Using the numerical differentiation method combined with the regularization technique, the first derivative of BA profiles is retrieved, and the height at which the first derivative of BA has the global minimum is defined to be the ABL height. To reflect the reliability of estimated ABL heights, the sharpness parameter is introduced, according to the relative minimum of the BA derivative. Then, it is applied to four months of COSMIC BA data(January, April, July, and October in 2008), and the ABL heights estimated are compared with two kinds of ABL heights from COSMIC products and with the heights determined by the finite difference method upon the refractivity data. For sharp ABL tops(large sharpness parameters), there is little difference between the ABL heights determined by different methods, i.e.,the uncertainties are small; whereas, for non-sharp ABL tops(small sharpness parameters), big differences exist in the ABL heights obtained by different methods, which means large uncertainties for different methods. In addition, the new method can detect thin ABLs and provide a reference ABL height in the cases eliminated by other methods. Thus, the application of the numerical differentiation method combined with the regularization technique to COSMIC BA data is an appropriate choice and has further application value.展开更多
Significant wave height is an important criterion in designing coastal and offshore structures.Based on the orthogonality principle, the linear mean square estimation method is applied to calculate significant wave he...Significant wave height is an important criterion in designing coastal and offshore structures.Based on the orthogonality principle, the linear mean square estimation method is applied to calculate significant wave height in this paper.Twenty-eight-year time series of wave data collected from three ocean buoys near San Francisco along the California coast are analyzed.It is proved theoretically that the computation error will be reduced by using as many measured data as possible for the calculation of significant wave height.Measured significant wave height at one buoy location is compared with the calculated value based on the data from two other adjacent buoys.The results indicate that the linear mean square estimation method can be well applied to the calculation and prediction of significant wave height in coastal regions.展开更多
Tsunami ran-up height is a significant parameter for dimensions of coastal structures. In the present study, tsunami run-up heights are estimated by three different Artificial Neural Network （ANN） models, i.e. Feed ...Tsunami ran-up height is a significant parameter for dimensions of coastal structures. In the present study, tsunami run-up heights are estimated by three different Artificial Neural Network （ANN） models, i.e. Feed Forward Back Propagation （FFBP）, Radial Basis Functions （RBF） and Generalized Regression Neural Network （GRNN）. As the input for the ANN configuration, the wave height （H） values are employed. It is shown that the tsunami ran-up height values are closely approximated with all of the applied ANN methods. The ANN estimations are slightly superior to those of the empirical equation. It can be seen that the ANN applications are especially significant in the absence of adequate number of laboratory experiments. The results also prove that the available experiment data set can be extended with ANN simulations. This may be helpful to decrease the burden of the experimental studies and to supply results for comparisons.展开更多
Shallow layer method(SLM)based on the definition of the geoid can determine the gravity field inside the shallow layer.In this study,the orthometric height of Mount Everest(HME)is calculated based on SLM,in which the ...Shallow layer method(SLM)based on the definition of the geoid can determine the gravity field inside the shallow layer.In this study,the orthometric height of Mount Everest(HME)is calculated based on SLM,in which the key is to construct the shallow layer model.The top and bottom boundaries of the shallow layer model are the natural surface of the Earth and the surface at a certain depth below the reference geoid,respectively.The model-combined strategies to determine the geoid undulation(N)based on SLM are applied to calculate the HME by two approaches:(1)direct calculation by combining N and geodetic height(h);(2)calculation by the segment summation approach(SSA)using the gravity field inside the shallow layer.On December 8,2020,the Chinese and Nepalese governments announced an authoritative value of 8848.86 m,which is referred to a geoid determined by the International Height Reference System(IHRS)(i.e.,the geopotential is 62636853.4 m^(2) s^(-2)).Here,our results(combined strategies(1)EGM2008 and CRUST1.0,(2)EGM2008 and CRUST2.0,(3)EIGEN-6 C4 and CRUST1.0,and(4)EIGEN-6 C4 and CRUST2.0)are referred to the geoid defined by WGS84(i.e.,the geopotential is 62636851.7 m^(2) s^(-2)).The differences between our results and the authoritative value(8848.86 m)are 0.448 m,-0.009 m,-0.295 m,and -0.741 m by the first approach,and 0.539 m,0.083 m,-0.214 m,and -0.647 m by the second approach.When the reference surface WGS84 geoid is converted to the IHRS geoid,the differences are 0.620 m,0.163 m,-0.123 m,and -0.569 m by the first approach,and0.711 m,0.225 m,-0.042 m,and -0.475 m by the second approach.展开更多
A method to determine the design value of the extreme wave height is introduced in this paper. The method is, based on the theoretical distribution of wave heights or the experiential probability distribution of wave ...A method to determine the design value of the extreme wave height is introduced in this paper. The method is, based on the theoretical distribution of wave heights or the experiential probability distribution of wave heights, to calculate the mode of the extreme wave height and the design value with good assuredness. In fitting the experiential distribution of the observed wave heights, the paper adopts the simulating function of F(x) = 1 - exp[ ?axb exp( -cxd )] and the non-linear parameter in the function is determined by the Non-Linear Least Square Method. The practical application of the method has proved that it is not only simple and convenient, but also very effective in engineering design practice.展开更多
Tweek atmospherics are extremely low frequency and very low frequency pulse signals with frequency dispersion characteristics that originate from lightning discharges and that propagate in the Earth–ionosphere wavegu...Tweek atmospherics are extremely low frequency and very low frequency pulse signals with frequency dispersion characteristics that originate from lightning discharges and that propagate in the Earth–ionosphere waveguide over long distances.In this study,we developed an automatic method to recognize tweek atmospherics and diagnose the lower ionosphere based on the machine learning method.The differences(automatic−manual)in each ionosphere parameter between the automatic method and the manual method were−0.07±2.73 km,0.03±0.92 cm^(−3),and 91±1,068 km for the ionospheric reflection height(h),equivalent electron densities at reflection heights(Ne),and propagation distance(d),respectively.Moreover,the automatic method is capable of recognizing higher harmonic tweek sferics.The evaluation results of the model suggest that the automatic method is a powerful tool for investigating the long-term variations in the lower ionosphere.展开更多
文摘In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.
基金supported by the Key Projects of Natural Science Foundation of China(No.41931284)the Scientific Research Start-Up Fund for High-Level Introduced Talents of Anhui University of Science and Technology(No.2022yjrc21).
文摘In the process of using the original key stratum theory to predict the height of a water-flowing fractured zone(WFZ),the influence of rock strata outside the calculation range on the rock strata within the calculation range as well as the fact that the shape of the overburden deformation area will change with the excavation length are ignored.In this paper,an improved key stratum theory(IKS theory)was proposed by fixing these two shortcomings.Then,a WFZ height prediction method based on IKS theory was established and applied.First,the range of overburden involved in the analysis was determined according to the tensile stress distribution range above the goaf.Second,the key stratum in the overburden involved in the analysis was identified through IKS theory.Finally,the tendency of the WFZ to develop upward was determined by judging whether or not the identified key stratum will break.The proposed method was applied and verified in a mining case study,and the reasons for the differences in the development patterns between the WFZs in coalfields in Northwest and East China were also fully explained by this method.
基金supported by Scientific Research Projects Management Coordinator of Kastamonu University,under grant number KÜ-BAP01/2019-41.
文摘Ecoregion-based height-diameter models were developed in the present study for Scots pine(Pinus sylves-tris L.)stands in Turkiye and included several ecological factors derived from a pre-existing ecoregional classification system.The data were obtained from 2831 sample trees in 292 sample plots.Ten generalized height–diameter models were developed,and the best model(HD10)was selected according to statistical criteria.Then,nonlinear mixed-effects modeling was applied to the best model.The R2 for the generalized height‒diameter model(Richards function)modified by Sharma and Parton is 0.951,and the final model included number of trees,dominant height,and diameter at breast height,with a random parameter associated with each ecoregion attached to the inverse of the mean basal area.The full model predictions using the nonlinear mixed-effects model and the reduced model(HD10)predictions were compared using the nonlinear sum of extra squares test,which revealed significant differences between ecore-gions;ecoregion-based height–diameter models were thus found to be suitable to use.In addition,using these models in appropriate ecoregions was very important for achieving reliable predictions with low prediction errors.
文摘Non-pillar continuous mining(NPCM) is regarded as a high-efficient, high-level and one-step mining technology, which can be divided into two substopes. Back fill stability status in substope I, which directly influence the loss rate and dilution rate, etc, will determine whether the experimental research is successful or not. By employing energy method of limit analysis and finite element numerical simulation method, the critical backfill height was determined under the prerequisite condition of its stability, which put forward theoretical basis for reasonable and correct selection of backfill’s parameters. The result showed that the first backfill could not keep stable for NPCM, while the other was able to.
基金National Natural Science Foundation of China(Grant Nos.51925502,51575150).
文摘To avoid impacts and vibrations during the processes of acceleration and deceleration while possessing flexible working ways for cable-suspended parallel robots(CSPRs),point-to-point trajectory planning demands an under-constrained cable-suspended parallel robot(UCPR)with variable angle and height cable mast as described in this paper.The end-effector of the UCPR with three cables can achieve three translational degrees of freedom(DOFs).The inverse kinematic and dynamic modeling of the UCPR considering the angle and height of cable mast are completed.The motion trajectory of the end-effector comprising six segments is given.The connection points of the trajectory segments(except for point P3 in the X direction)are devised to have zero instantaneous velocities,which ensure that the acceleration has continuity and the planned acceleration curve achieves smooth transition.The trajectory is respectively planned using three algebraic methods,including fifth degree polynomial,cycloid trajectory,and double-S velocity curve.The results indicate that the trajectory planned by fifth degree polynomial method is much closer to the given trajectory of the end-effector.Numerical simulation and experiments are accomplished for the given trajectory based on fifth degree polynomial planning.At the points where the velocity suddenly changes,the length and tension variation curves of the planned and unplanned three cables are compared and analyzed.The OptiTrack motion capture system is adopted to track the end-effector of the UCPR during the experiment.The effectiveness and feasibility of fifth degree polynomial planning are validated.
文摘This paper presents a refined method for estimating the annual extreme wave heights at a coastal or offshore project site on the basis of the data acquired at some nearby routine hydrographic stations. This method is based on the orthogonality principle in linear mean square estimation of stochastic processes. The error of the method is analyzed and compared with that of the conventional method. It is found that the method is able to effectively reduce the error so long as some feasible measures are adopted. A simulated test of the method has been conducted in a large scale wind wave flume. The test results are in good agreement with those given by theoretical error analysis. A scheme to implement the method is proposed on the basis of error analysis. The scheme is so designed as to reduce the estimation error as far as possible. This method is also suitable to utilizing satellite wave data for the estimation.
基金supported by the National Natural Science Foundation of China (Grant No. 41475021)
文摘This paper presents a new method to estimate the height of the atmospheric boundary layer(ABL) by using COSMIC radio occultation bending angle(BA) data. Using the numerical differentiation method combined with the regularization technique, the first derivative of BA profiles is retrieved, and the height at which the first derivative of BA has the global minimum is defined to be the ABL height. To reflect the reliability of estimated ABL heights, the sharpness parameter is introduced, according to the relative minimum of the BA derivative. Then, it is applied to four months of COSMIC BA data(January, April, July, and October in 2008), and the ABL heights estimated are compared with two kinds of ABL heights from COSMIC products and with the heights determined by the finite difference method upon the refractivity data. For sharp ABL tops(large sharpness parameters), there is little difference between the ABL heights determined by different methods, i.e.,the uncertainties are small; whereas, for non-sharp ABL tops(small sharpness parameters), big differences exist in the ABL heights obtained by different methods, which means large uncertainties for different methods. In addition, the new method can detect thin ABLs and provide a reference ABL height in the cases eliminated by other methods. Thus, the application of the numerical differentiation method combined with the regularization technique to COSMIC BA data is an appropriate choice and has further application value.
基金support for this study was provided by the National Natural Science Foundation of China (No.40776006)Research Fund for the Doctoral Program of Higher Education of China (Grant No.20060423009)the Science and Technology Development Program of Shandong Province (Grant No.2008GGB01099)
文摘Significant wave height is an important criterion in designing coastal and offshore structures.Based on the orthogonality principle, the linear mean square estimation method is applied to calculate significant wave height in this paper.Twenty-eight-year time series of wave data collected from three ocean buoys near San Francisco along the California coast are analyzed.It is proved theoretically that the computation error will be reduced by using as many measured data as possible for the calculation of significant wave height.Measured significant wave height at one buoy location is compared with the calculated value based on the data from two other adjacent buoys.The results indicate that the linear mean square estimation method can be well applied to the calculation and prediction of significant wave height in coastal regions.
文摘Tsunami ran-up height is a significant parameter for dimensions of coastal structures. In the present study, tsunami run-up heights are estimated by three different Artificial Neural Network （ANN） models, i.e. Feed Forward Back Propagation （FFBP）, Radial Basis Functions （RBF） and Generalized Regression Neural Network （GRNN）. As the input for the ANN configuration, the wave height （H） values are employed. It is shown that the tsunami ran-up height values are closely approximated with all of the applied ANN methods. The ANN estimations are slightly superior to those of the empirical equation. It can be seen that the ANN applications are especially significant in the absence of adequate number of laboratory experiments. The results also prove that the available experiment data set can be extended with ANN simulations. This may be helpful to decrease the burden of the experimental studies and to supply results for comparisons.
基金supported in part by the National Natural Science Foundations of China under Grants Nos.41631072,42030105,41721003,41804012,and 41874023。
文摘Shallow layer method(SLM)based on the definition of the geoid can determine the gravity field inside the shallow layer.In this study,the orthometric height of Mount Everest(HME)is calculated based on SLM,in which the key is to construct the shallow layer model.The top and bottom boundaries of the shallow layer model are the natural surface of the Earth and the surface at a certain depth below the reference geoid,respectively.The model-combined strategies to determine the geoid undulation(N)based on SLM are applied to calculate the HME by two approaches:(1)direct calculation by combining N and geodetic height(h);(2)calculation by the segment summation approach(SSA)using the gravity field inside the shallow layer.On December 8,2020,the Chinese and Nepalese governments announced an authoritative value of 8848.86 m,which is referred to a geoid determined by the International Height Reference System(IHRS)(i.e.,the geopotential is 62636853.4 m^(2) s^(-2)).Here,our results(combined strategies(1)EGM2008 and CRUST1.0,(2)EGM2008 and CRUST2.0,(3)EIGEN-6 C4 and CRUST1.0,and(4)EIGEN-6 C4 and CRUST2.0)are referred to the geoid defined by WGS84(i.e.,the geopotential is 62636851.7 m^(2) s^(-2)).The differences between our results and the authoritative value(8848.86 m)are 0.448 m,-0.009 m,-0.295 m,and -0.741 m by the first approach,and 0.539 m,0.083 m,-0.214 m,and -0.647 m by the second approach.When the reference surface WGS84 geoid is converted to the IHRS geoid,the differences are 0.620 m,0.163 m,-0.123 m,and -0.569 m by the first approach,and0.711 m,0.225 m,-0.042 m,and -0.475 m by the second approach.
文摘A method to determine the design value of the extreme wave height is introduced in this paper. The method is, based on the theoretical distribution of wave heights or the experiential probability distribution of wave heights, to calculate the mode of the extreme wave height and the design value with good assuredness. In fitting the experiential distribution of the observed wave heights, the paper adopts the simulating function of F(x) = 1 - exp[ ?axb exp( -cxd )] and the non-linear parameter in the function is determined by the Non-Linear Least Square Method. The practical application of the method has proved that it is not only simple and convenient, but also very effective in engineering design practice.
基金supported by the Chinese Academy of Sciences(CAS)Project of Stable Support for Youth Team in Basic Research Field(Grant No.YSRR-018)the National Key R&D Program of China(Grant No.2019YFC1510103)+1 种基金the National Natural Science Foundation of China(Grant Nos.41875006 and U1938115)the Chinese Meridian Project,and the International Partnership Program of CAS(Grant No.183311KYSB20200003).
文摘Tweek atmospherics are extremely low frequency and very low frequency pulse signals with frequency dispersion characteristics that originate from lightning discharges and that propagate in the Earth–ionosphere waveguide over long distances.In this study,we developed an automatic method to recognize tweek atmospherics and diagnose the lower ionosphere based on the machine learning method.The differences(automatic−manual)in each ionosphere parameter between the automatic method and the manual method were−0.07±2.73 km,0.03±0.92 cm^(−3),and 91±1,068 km for the ionospheric reflection height(h),equivalent electron densities at reflection heights(Ne),and propagation distance(d),respectively.Moreover,the automatic method is capable of recognizing higher harmonic tweek sferics.The evaluation results of the model suggest that the automatic method is a powerful tool for investigating the long-term variations in the lower ionosphere.
