Shenzhen,a major city in southern China,has experienced rapid advancements in Unmanned Aerial Vehicle(UAV)technology,resulting in extensive logistics networks with thousands of daily flights.However,frequent disruptio...Shenzhen,a major city in southern China,has experienced rapid advancements in Unmanned Aerial Vehicle(UAV)technology,resulting in extensive logistics networks with thousands of daily flights.However,frequent disruptions due to its subtropical monsoon climate,including typhoons and gusty winds,present ongoing challenges.Despite the growing focus on operational costs and third-party risks,research on low-altitude urban wind fields remains scarce.This study addresses this gap by integrating wind field analysis into UAV path planning,introducing key innovations to the classical model.First,UAV wind resistance and turbulence constraints are analyzed,mapping high-wind-speed and turbulence-prone zones in the airspace.Second,wind dynamics are incorporated into path planning by considering airspeed and groundspeed variation,optimizing waypoint selection and flight speed adjustments to improve overall energy efficiency.Additionally,a wind-aware Theta*algorithm is proposed,leveraging wind vectors to expedite search process,while Computational Fluid Dynamics(CFD)techniques are employed to calculate wind fields.A case study of Shenzhen,examining wind patterns over the past decade,demonstrates a 6.23%improvement in groundspeed and a 7.69%reduction in energy consumption compared to wind-agnostic models.This framework advances UAV logistics by enhancing route safety and energy efficiency,contributing to more cost-effective operations.展开更多
A downburst is a strong downdraft generated by intense thunderstorm clouds,producing radially divergent and highly destructive winds near the ground.Its characteristic scales are expressed through random variations in...A downburst is a strong downdraft generated by intense thunderstorm clouds,producing radially divergent and highly destructive winds near the ground.Its characteristic scales are expressed through random variations in jet height,velocity,and diameter during an event.In this study,a reduced-scale parked wind turbine is exposed to downburst wind fields to investigate the resulting extreme wind loads.The analysis emphasizes both the flow structure of downbursts and the variations of surface wind pressure on turbine blades under different jet parameters.Results show that increasing jet velocity markedly enhances the maximum horizontal wind speed,while greater jet height reduces the horizontal wind speed and shifts the peak velocity closer to the jet center.Increasing jet diameter primarily affects the radial position of the maximum horizontal wind speed.For the wind turbine,the maximum equivalent stress and blade displacement increase almost linearly with jet velocity,but exhibit the opposite trend with jet diameter.Specifically,as jet velocity rises from 10 m/s to 20 m/s,the surface pressure coefficient at the blade tip increases by approximately 4.5 times.Changes in jet diameter indirectly alter the turbine’s relative position within the wind field,leading to variations in wind load direction and exposure area.Conversely,increasing jet height extends the dissipation path of the downdraft,thereby reducing the intensity of the airflow acting on the blades.For example,when jet height increases from 0.3 m to 1.2 m,the surface pressure coefficient at the blade tip decreases by nearly 50%.展开更多
A three-dimensional wind field analysis sollware based on the Beigng-Gucheng dual-Doppler weather radar system has been built, and evaluated by using the numerical cloud model producing storm flow and hydrometeor fiel...A three-dimensional wind field analysis sollware based on the Beigng-Gucheng dual-Doppler weather radar system has been built, and evaluated by using the numerical cloud model producing storm flow and hydrometeor fields. The effects of observation noise and the spatial distribution of wind field analysis error are also investigated.展开更多
Stress raisers such as holes are inevitable in structures at which stress concentration occurs and the static as well as fatigue strength of the structures can be significantly weakened.Therefore,to accurately evaluat...Stress raisers such as holes are inevitable in structures at which stress concentration occurs and the static as well as fatigue strength of the structures can be significantly weakened.Therefore,to accurately evaluate the stress concentration factor and stress fields at holes is of essential importance for structure design and service life prediction.Although stress and strain concentration and fields at holes in finite thickness plates strongly change with and along the thickness,manuals of stress concentration for engineering design are mainly based on twodimensional theory and no explicit formula is available even for circular holes in finite thickness plates.Here we obtain for the first time a complete set of explicit formulae for stress and strain concentration factors and the out-of-plane constraint factor at circular as well as elliptical holes in finite thickness plates by integrating comprehensive three-dimensional finite element analyses and available theoretical solutions.The three-dimensional stress distributions ahead of holes can also be predicted by the obtained formulae.With their accuracy and the corresponding applicable range being analyzed and outlined in detail,the formulae can serve as an important fundamental solution for three-dimensional engineering structure design and guideline for developing threedimensional analytical methods.展开更多
In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through min...In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through minimizing a cost function defined as a radial observation term with the standard conjugate gradient method,avoiding a weighting parameter specification step.Compared with conventional dual-Doppler wind synthesis approaches,this variational method minimizes errors caused by interpolation from radar observation to analysis grid in the iterative solution process,which is one of the main sources of errors.Then,through the accelerated Liebmann method,the vertical velocity is further reestimated as an extra step by solving the Poisson equation with impermeable conditions imposed at the ground and near the tropopause.The Poisson equation defined by the second derivative of the vertical velocity is derived from the mass continuity equation.Compared with the method proposed by O’Brien,this method is less sensitive to the uncertainty of the boundary conditions and has better stability and reliability.Furthermore,the method proposed in this paper is applied to Doppler radar observation of a squall line process.It is shown that the retrieved vertical wind profile agrees well with the vertical profile obtained with the velocity–azimuth display(VAD)method,and the retrieved radial velocity as well as the analyzed positive and negative velocity centers and horizontal wind shear of the squall line are in accord with radar observations.There is a good correspondence between the divergence field of the derived wind field and the vertical velocity.And,the horizontal and vertical circulations within and around the squall line,as well as strong updrafts,the associated downdrafts,and associated rear inflow of the bow echo,are analyzed well.It is worth mentioning that the variational method in this paper can be applied to simultaneously synthesize the three-dimensional wind field from multiple-Doppler radar observations.展开更多
The theoretical solutions are obtained for the three-dimensional(3-D)stress field in an infinite isotropic elastic plate with a through-the-thickness circular hole subjected to shear load at far field by using Kane an...The theoretical solutions are obtained for the three-dimensional(3-D)stress field in an infinite isotropic elastic plate with a through-the-thickness circular hole subjected to shear load at far field by using Kane and Mindlin′s assumption based on the stress function method.Based on the present solutions,the characteristics of 3-D stress field are analyzed and the emphasis is placed on the effects of the plate thickness and Poisson′s ratio on the deviation of the present 3-D in-plane stress from the related plane stress solutions,the stress concentration and the out-of-plane constraint.The present solutions show that the stress concentration factor reaches its peak value of about 8.9% which is higher than that of the plane stress solutions.As expected,the out-of-plane stress constraint factor can reach 1on the surface of the hole when the plate is a very thick one.展开更多
This study employed a computational fluid dynamics model with an overset mesh technique to investigate the thrust and power of a floating offshore wind turbine(FOWT)under platform floating motion in the wind–rain fie...This study employed a computational fluid dynamics model with an overset mesh technique to investigate the thrust and power of a floating offshore wind turbine(FOWT)under platform floating motion in the wind–rain field.The impact of rainfall on aerodynamic performance was initially examined using a stationary turbine model in both wind and wind–rain conditions.Subsequently,the study compared the FOWT’s performance under various single degree-of-freedom(DOF)motions,including surge,pitch,heave,and yaw.Finally,the combined effects of wind–rain fields and platform motions involving two DOFs on the FOWT’s aerodynamics were analyzed and compared.The results demonstrate that rain negatively impacts the aerodynamic performance of both the stationary turbines and FOWTs.Pitch-dominated motions,whether involving single or multiple DOFs,caused significant fluctuations in the FOWT aerodynamics.The combination of surge and pitch motions created the most challenging operational environment for the FOWT in all tested scenarios.These findings highlighted the need for stronger construction materials and greater ultimate bearing capacity for FOWTs,as well as the importance of optimizing designs to mitigate excessive pitch and surge.展开更多
The accuracy of the Earth＇s gravitational field measured from the gravity field and steady-state ocean circulation explorer（GOCE）,up to 250 degrees,influenced by the radial gravity gradient V zz and three-dimension...The accuracy of the Earth＇s gravitational field measured from the gravity field and steady-state ocean circulation explorer（GOCE）,up to 250 degrees,influenced by the radial gravity gradient V zz and three-dimensional gravity gradient V ij from the satellite gravity gradiometry（SGG） are contrastively demonstrated based on the analytical error model and numerical simulation,respectively.Firstly,the new analytical error model of the cumulative geoid height,influenced by the radial gravity gradient V zz and three-dimensional gravity gradient V ij are established,respectively.In 250 degrees,the GOCE cumulative geoid height error measured by the radial gravity gradient V zz is about 2 1/2 times higher than that measured by the three-dimensional gravity gradient V ij.Secondly,the Earth＇s gravitational field from GOCE completely up to 250 degrees is recovered using the radial gravity gradient V zz and three-dimensional gravity gradient V ij by numerical simulation,respectively.The study results show that when the measurement error of the gravity gradient is 3×10 12 /s 2,the cumulative geoid height errors using the radial gravity gradient V zz and three-dimensional gravity gradient V ij are 12.319 cm and 9.295 cm at 250 degrees,respectively.The accuracy of the cumulative geoid height using the three-dimensional gravity gradient V ij is improved by 30%-40% on average compared with that using the radial gravity gradient V zz in 250 degrees.Finally,by mutual verification of the analytical error model and numerical simulation,the orders of magnitude from the accuracies of the Earth＇s gravitational field recovery make no substantial differences based on the radial and three-dimensional gravity gradients,respectively.Therefore,it is feasible to develop in advance a radial cold-atom interferometric gradiometer with a measurement accuracy of 10 13 /s 2-10 15 /s 2 for precisely producing the next-generation GOCE Follow-On Earth gravity field model with a high spatial resolution.展开更多
The wind and temperature fields at 20 to 55 km above the Antigua launch site(17°N,61°W)were analyzed by using sounding rocket data published by the research organization on Stratosphere-Troposphere Processes...The wind and temperature fields at 20 to 55 km above the Antigua launch site(17°N,61°W)were analyzed by using sounding rocket data published by the research organization on Stratosphere-Troposphere Processes and their Role in Climate(SPARC).The results showed distinct variations in the wind and temperature fields at different heights from the 1960s to the 1990s.The overall zonal wind speed showed a significant increasing trend with the year,and the overall change in meridional wind speed showed a falling trend from 1976 to 1990,whereas the temperature field showed a significant cooling trend from 1964 to 1990.The times the trends mutated varied at different levels.By taking the altitudes at 20,35,and 50 km as representatives,we found that the zonal wind speed trend had mutated in 1988,1986,and 1986,respectively;that the meridional wind speed trend had mutated in 1990,1986,and 1990,respectively;and that the temperature trend had mutated separately in 1977,1973,and 1967,respectively.Characteristics of the periodic wind and temperature field variations at different heights were also analyzed,and obvious differences were found in time scale variations across the different layers.The zonal and meridional wind fields were basically characterized as having a significant periodic variation of 5 years across the three layers,and each level was characterized as having a periodic variation of less than 5 years.Temperature field variation at the three levels was basically characterized as occurring in 10-year and 5-year cycles.展开更多
In the Jiaoshiba block of the Fuling shale gas field,the employed reserves and recovery factor by primary well pattern are low,no obvious barrier is found in the development layer series,and layered development is dif...In the Jiaoshiba block of the Fuling shale gas field,the employed reserves and recovery factor by primary well pattern are low,no obvious barrier is found in the development layer series,and layered development is difficult.Based on the understanding of the main factors controlling shale gas enrichment and high production,the theory and technology of shale gas three-dimensional development,such as fine description and modeling of shale gas reservoir,optimization of three-dimensional development strategy,highly efficient drilling with dense well pattern,precision fracturing and real-time control,are discussed.Three-dimensional development refers to the application of optimal and fast drilling and volume fracturing technologies,depending upon the sedimentary characteristics,reservoir characteristics and sweet spot distribution of shale gas,to form"artificial gas reservoir"in a multidimensional space,so as to maximize the employed reserves,recovery factor and yield rate of shale gas development.In the research on shale gas three-dimensional development,the geological+engineering sweet spot description is fundamental,the collaborative optimization of natural fractures and artificial fractures is critical,and the improvement of speed and efficiency in drilling and fracturing engineering is the guarantee.Through the implementation of three-dimensional development,the overall recovery factor in the Jiaoshiba block has increased from 12.6%to 23.3%,providing an important support for the continuous and stable production of the Fuling shale gas field.展开更多
Three-dimensional holographic vector of atomic interaction field（3D-HoVAIF） is used to describe the chemical structures of polychlorinated naphthalenes（PCNs）.After variable screening by stepwise multiple regressio...Three-dimensional holographic vector of atomic interaction field（3D-HoVAIF） is used to describe the chemical structures of polychlorinated naphthalenes（PCNs）.After variable screening by stepwise multiple regression（SMR） technique,the liner relationships between gas-chromatographic relative retention time（RRT）,298 K supercooled liquid pressures（logPL）,n-octanol/air partition coefficient（logKOA）,n-octanol/water partition coefficient（logKOW）,aqueous solubilities（logSW）,relative in vitro potency values（-logEROD） of PCNs and 3D-HoVAIF descriptors have been established by partial least-square（PLS） regression.The result shows that the 3D-HoVAIF descriptors can be well used to express the quantitative structure-property（activity） relationships of PCNs.Predictive capability of the models has also been demonstrated by leave-one-out cross-validation.Moreover,the predicted values have been presented for those PCNs which are lack of experimentally physico-chemical properties and biological activity by the optimum models.展开更多
The optimal path planning for fixed-wing unmanned aerial vehicles(UAVs) in multi-target surveillance tasks(MTST) in the presence of wind is concerned.To take into account the minimal turning radius of UAVs,the Dubins ...The optimal path planning for fixed-wing unmanned aerial vehicles(UAVs) in multi-target surveillance tasks(MTST) in the presence of wind is concerned.To take into account the minimal turning radius of UAVs,the Dubins model is used to approximate the dynamics of UAVs.Based on the assumption,the path planning problem of UAVs in MTST can be formulated as a Dubins traveling salesman problem(DTSP).By considering its prohibitively high computational cost,the Dubins paths under terminal heading relaxation are introduced,which leads to significant reduction of the optimization scale and difficulty of the whole problem.Meanwhile,in view of the impact of wind on UAVs' paths,the notion of virtual target is proposed.The application of the idea successfully converts the Dubins path planning problem from an initial configuration to a target in wind into a problem of finding the minimal root of a transcendental equation.Then,the Dubins tour is derived by using differential evolution(DE) algorithm which employs random-key encoding technique to optimize the visiting sequence of waypoints.Finally,the effectiveness and efficiency of the proposed algorithm are demonstrated through computational experiments.Numerical results exhibit that the proposed algorithm can produce high quality solutions to the problem.展开更多
An experimental investigation of three-dimensional flow field in a film-cooled turbine model is carried out by using particle image velocimeter （PIV） in a low-speed wind tunnel. The effects of different blowing rati...An experimental investigation of three-dimensional flow field in a film-cooled turbine model is carried out by using particle image velocimeter （PIV） in a low-speed wind tunnel. The effects of different blowing ratios （M=1.5, 2） on the flow field are studied. The experimental results reveal the classical phenomena of the formation of kidney vortex pair and secondary flow in wake region behind the jet hole. And the changes of the kidney vortex pair and the wake at different locations away from the hole on the suction and pressure sides are also studied. Compared with the flow field in stationary cascade, there are centrifugal force and Coriolis force existing in the flow field of rotating turbine, and these forces bring the radial velocity in the jet flow. The effect of rotatien on the flow field of the pressure side is more distinct than that on the suction side from the measured flow fields in Y-Z plane and radial velocity contours. The increase of blowing ratio makes the kidney vortex pair and the secondary flow in the wake region stronger and makes the range of the wake region enlarged.展开更多
Typhoon Chaba was the most intense typhoon to strike western Guangdong since Typhoon Mujigae in 2015.According to the National Disaster Reduction Center of China,in the morning of July 7,2022,over 1.5 million people i...Typhoon Chaba was the most intense typhoon to strike western Guangdong since Typhoon Mujigae in 2015.According to the National Disaster Reduction Center of China,in the morning of July 7,2022,over 1.5 million people in Guangdong,Guangxi,and Hainan were affected by Typhoon Chaba.The typhoon also caused the“Fukui 001”ship to be in distress in the waters near Yangjiang,Guangdong,on July 2,resulting in big casualties.Studies have indicated that wind field forecast for Typhoon Chaba was not accurate.To better simulate typhoon events and assess their impacts,we proposed the use of a model wind field(Fujita-Takahashi)integrated with the Copernicus Marine and Environmental Monitoring Service(CMEMS)data to reconstruct effectively the overall wind field of Typhoon Chaba.The simulation result aligns well with the observations,particularly at the Dashu Island Station,showing consistent trends in wind speed changes.However,certain limitations were noted.The model shows that the attenuation of wind speed is slower when typhoon neared land than that observed,indicating that the model has a high simulation accuracy for the ocean wind field,but may have deviations near coastal areas.The result is accurate for open sea but deviated for near land due to the land friction effect.Therefore,we recommend to adjust the model to improve the accuracy for near coasts.展开更多
Accurate short-term forecast of offshore wind fields is still challenging for numerical weather prediction models.Based on three years of 48-hour forecast data from the European Centre for Medium-Range Weather Forecas...Accurate short-term forecast of offshore wind fields is still challenging for numerical weather prediction models.Based on three years of 48-hour forecast data from the European Centre for Medium-Range Weather Forecasts Integrated Forecasting System global model(ECMWF-IFS)over 14 offshore weather stations along the coast of Shandong Province,this study introduces a multi-task learning(MTL)model(TabNet-MTL),which significantly improves the forecast bias of near-surface wind direction and speed simultaneously.TabNet-MTL adopts the feature engineering method,utilizes mean square error as the loss function,and employs the 5-fold cross validation method to ensure the generalization ability of the trained model.It demonstrates superior skills in wind field correction across different forecast lead times over all stations compared to its single-task version(TabNet-STL)and three other popular single-task learning models(Random Forest,LightGBM,and XGBoost).Results show that it significantly reduces root mean square error of the ECMWF-IFS wind speed forecast from 2.20 to 1.25 m s−1,and increases the forecast accuracy of wind direction from 50%to 65%.As an explainable deep learning model,the weather stations and long-term temporal statistics of near-surface wind speed are identified as the most influential variables for TabNet-MTL in constructing its feature engineering.展开更多
Proper knowledge of the nature of geomagnetic storms and their relationships with the conditions of the space environment at the outer part of the Earth's magnetosphere(bow shock nose) is essential to increase our...Proper knowledge of the nature of geomagnetic storms and their relationships with the conditions of the space environment at the outer part of the Earth's magnetosphere(bow shock nose) is essential to increase our resilience to space weather disturbances. In this article, we present an analysis of the interplanetary magnetic field(IMF) and solar wind parameters relevant to 100 geomagnetic storms in Solar Cycle 24. We revisit the relationship between the minimum disturbance storm time index(Dst_(min)), the minimum southward IMF(B_(S, min)), the maximum solar wind density(N_(SW, max)) and speed(V_(max)), and the lag time between the extrema(dT(B_(z), N),dT(B_(z), V)). We end with a regression formula that fits the data, with a coefficient of determination of 0.58, a root mean square error of 21.30 nT, and a mean absolute error of 15.87 nT. Even though more complex machine learning models can outperform this model, it serves as a theoretically sensible alternative for understanding and forecasting geomagnetic storms.展开更多
Stratospheric airships are lighter-than-air vehicles capable of continuous flying for months.The energy balance of the airship is the key to long-duration flights.The stratospheric airship is entirely powered by the s...Stratospheric airships are lighter-than-air vehicles capable of continuous flying for months.The energy balance of the airship is the key to long-duration flights.The stratospheric airship is entirely powered by the solar array.It is necessary to accurately predict the output power of the array for any flight state.Because of the uneven solar radiation received by the solar array,the traditional model based on components has a slow simulation speed.In this study,a data-driven surrogate modeling approach for prediction the output power of the solar array is proposed.The surrogate model is trained using the samples obtained from the high-accuracy simulation model.By using the input parameter preprocessor,the accuracy of the surrogate model in predicting the output power of the solar array is improved to 98.65%.In addition,the predictive speed of the surrogate model is ten million times faster than the traditional simulation model.Finally,the surrogate model is used to predict the energy balance of stratospheric airships flying throughout the year under actual global wind fields.展开更多
Stochastic simulation is an important means of acquiring fluctuating wind pressures for wind induced response analyses in structural engineering. The wind pressure acting on a large-span space structure can be charact...Stochastic simulation is an important means of acquiring fluctuating wind pressures for wind induced response analyses in structural engineering. The wind pressure acting on a large-span space structure can be characterized as a stationary non-Gaussian field. This paper reviews several simulation algorithms related to the Spectral Representation Method （SRM） and the Static Transformation Method （STM）. Polynomial and Exponential transformation functions （PSTM and ESTM） are discussed. Deficiencies in current algorithms, with respect to accuracy, stability and efficiency, are analyzed, and the algorithms are improved for better practical application. In order to verify the improved algorithm, wind pressure fields on a large-span roof are simulated and compared with wind tunnel data. The simulation results fit well with the wind tunnel data, and the algorithm accuracy, stability and efficiency are shown to be better than those of current algorithms.展开更多
The core sampling experiments were conducted after hydraulic fracturing in the three-dimensional development zone of Fuling shale gas.Six coring wells of different well types were systematically designed.Based on the ...The core sampling experiments were conducted after hydraulic fracturing in the three-dimensional development zone of Fuling shale gas.Six coring wells of different well types were systematically designed.Based on the integrated engineering technology of post-fracturing drilling,coring and monitoring of shale and the analysis of fracture source tracing,the evaluation of the fracture network after fracturing in the three-dimensional development of shale gas was conducted.The data of core fractures after fracturing indicate that three major types of fractures are formed after fracturing:natural fractures,hydraulic fractures,and fractures induced by external mechanical force,which are further classified into six subcategories:natural structural fractures,natural bedding fractures,hydraulic fractures,hydraulically activated fractures,drilling induced fractures,and fractures induced by core transportation.The forms of the artificial fracture network after fracturing are complex.Hydraulic fractures and hydraulically activated fractures interweave with each other,presenting eight forms of artificial fracture networks,among which the“一”-shaped fracture is the most common,accounting for approximately 70%of the total fractures.When the distance to the fractured wellbore is less than 35 m,the density of the artificial fracture network is relatively high;when it is 35–100 m,the density is lower;and when it is beyond 100 m,the density gradually increases.The results of the fracture tracing in the core sampling area confirm that the current fracturing technology can essentially achieve the differential transformation of the reservoir in the main area of Jiaoshiba block in Fuling.The three-layer three-dimensional development model can efficiently utilize shale gas reserves,although there is still room for improvement in the complexity and propagation uniformity of fractures.It is necessary to further optimize technologies such as close-cutting combined with temporary blocking and deflection within fractures or at fracture mouths,as well as limited flow perforation techniques,to promote the balanced initiation and extension of fractures.展开更多
We combined domestic ground-based and satellite magnetic measurements to create a regional three-dimensional surface Spline(3DSS)gradient model of the main geomagnetic field over the Chinese continent.To improve the p...We combined domestic ground-based and satellite magnetic measurements to create a regional three-dimensional surface Spline(3DSS)gradient model of the main geomagnetic field over the Chinese continent.To improve the precision of the model,we considered the data gap between the ground and satellite data.We compared and analyzed the results of the Taylor polynomial,surface Spline,and CHAOS-6(the CHAMP,?rsted and SAC-C model of Earth’s magnetic field)gradient models.Results showed that the gradients in the south-north and east-west directions of the four models were consistent.The 3DSS model was able to express not only gradients at different altitudes,but also average gradients inside the research area.The two Spline models were able to capture more information on gradient anomalies than were the fitted models.Strong local anomalies were observed in northern Xinjiang,Beijing,and the junction area between Jiangsu and Zhejiang,and the total intensity F decreased whereas the altitude increased.The gradient decreased by 21.69%in the south-north direction and increased by 11.78%in the east-west direction.In addition,the altitude gradient turned from negative to positive while the altitude increased.The Spline model and the two fitted models differed mainly in the field sources they expressed and the modeling theory.展开更多
基金supported by the National Natural Science Foundation of China(No.U2433214)。
文摘Shenzhen,a major city in southern China,has experienced rapid advancements in Unmanned Aerial Vehicle(UAV)technology,resulting in extensive logistics networks with thousands of daily flights.However,frequent disruptions due to its subtropical monsoon climate,including typhoons and gusty winds,present ongoing challenges.Despite the growing focus on operational costs and third-party risks,research on low-altitude urban wind fields remains scarce.This study addresses this gap by integrating wind field analysis into UAV path planning,introducing key innovations to the classical model.First,UAV wind resistance and turbulence constraints are analyzed,mapping high-wind-speed and turbulence-prone zones in the airspace.Second,wind dynamics are incorporated into path planning by considering airspeed and groundspeed variation,optimizing waypoint selection and flight speed adjustments to improve overall energy efficiency.Additionally,a wind-aware Theta*algorithm is proposed,leveraging wind vectors to expedite search process,while Computational Fluid Dynamics(CFD)techniques are employed to calculate wind fields.A case study of Shenzhen,examining wind patterns over the past decade,demonstrates a 6.23%improvement in groundspeed and a 7.69%reduction in energy consumption compared to wind-agnostic models.This framework advances UAV logistics by enhancing route safety and energy efficiency,contributing to more cost-effective operations.
基金the National Natural Science Foundation of China(Grant Nos.52276197, 52166014)Gansu Province Key Research and Development Program—Industrial Project(Grant No.23YFGA0069).
文摘A downburst is a strong downdraft generated by intense thunderstorm clouds,producing radially divergent and highly destructive winds near the ground.Its characteristic scales are expressed through random variations in jet height,velocity,and diameter during an event.In this study,a reduced-scale parked wind turbine is exposed to downburst wind fields to investigate the resulting extreme wind loads.The analysis emphasizes both the flow structure of downbursts and the variations of surface wind pressure on turbine blades under different jet parameters.Results show that increasing jet velocity markedly enhances the maximum horizontal wind speed,while greater jet height reduces the horizontal wind speed and shifts the peak velocity closer to the jet center.Increasing jet diameter primarily affects the radial position of the maximum horizontal wind speed.For the wind turbine,the maximum equivalent stress and blade displacement increase almost linearly with jet velocity,but exhibit the opposite trend with jet diameter.Specifically,as jet velocity rises from 10 m/s to 20 m/s,the surface pressure coefficient at the blade tip increases by approximately 4.5 times.Changes in jet diameter indirectly alter the turbine’s relative position within the wind field,leading to variations in wind load direction and exposure area.Conversely,increasing jet height extends the dissipation path of the downdraft,thereby reducing the intensity of the airflow acting on the blades.For example,when jet height increases from 0.3 m to 1.2 m,the surface pressure coefficient at the blade tip decreases by nearly 50%.
文摘A three-dimensional wind field analysis sollware based on the Beigng-Gucheng dual-Doppler weather radar system has been built, and evaluated by using the numerical cloud model producing storm flow and hydrometeor fields. The effects of observation noise and the spatial distribution of wind field analysis error are also investigated.
基金supported by the National Natural Science Foundation of China(51535005,51472117)the Research Fund of State Key Laboratory of Mechanics and Control of Mechanical Structures(MCMS-I-0418K01,MCMS-I-0418Y01,MCMS-0417G02,MCMS-0417G03)+1 种基金the Fundamental Research Funds for the Central Universities(NP2017101,NC2018001)a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions.The authors would like to thank Dr.Chongmin She for helpful discussions.
文摘Stress raisers such as holes are inevitable in structures at which stress concentration occurs and the static as well as fatigue strength of the structures can be significantly weakened.Therefore,to accurately evaluate the stress concentration factor and stress fields at holes is of essential importance for structure design and service life prediction.Although stress and strain concentration and fields at holes in finite thickness plates strongly change with and along the thickness,manuals of stress concentration for engineering design are mainly based on twodimensional theory and no explicit formula is available even for circular holes in finite thickness plates.Here we obtain for the first time a complete set of explicit formulae for stress and strain concentration factors and the out-of-plane constraint factor at circular as well as elliptical holes in finite thickness plates by integrating comprehensive three-dimensional finite element analyses and available theoretical solutions.The three-dimensional stress distributions ahead of holes can also be predicted by the obtained formulae.With their accuracy and the corresponding applicable range being analyzed and outlined in detail,the formulae can serve as an important fundamental solution for three-dimensional engineering structure design and guideline for developing threedimensional analytical methods.
基金supported by the National Key Research and Development Program of China(Grant No.2019YFC1510400)the National Natural Science Foundation of China(Grant Nos.41975054 and 41930967)the Special Fund for Forecasters of China Meteorological Administration(Grant No.CMAYBY2018-040)。
文摘In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through minimizing a cost function defined as a radial observation term with the standard conjugate gradient method,avoiding a weighting parameter specification step.Compared with conventional dual-Doppler wind synthesis approaches,this variational method minimizes errors caused by interpolation from radar observation to analysis grid in the iterative solution process,which is one of the main sources of errors.Then,through the accelerated Liebmann method,the vertical velocity is further reestimated as an extra step by solving the Poisson equation with impermeable conditions imposed at the ground and near the tropopause.The Poisson equation defined by the second derivative of the vertical velocity is derived from the mass continuity equation.Compared with the method proposed by O’Brien,this method is less sensitive to the uncertainty of the boundary conditions and has better stability and reliability.Furthermore,the method proposed in this paper is applied to Doppler radar observation of a squall line process.It is shown that the retrieved vertical wind profile agrees well with the vertical profile obtained with the velocity–azimuth display(VAD)method,and the retrieved radial velocity as well as the analyzed positive and negative velocity centers and horizontal wind shear of the squall line are in accord with radar observations.There is a good correspondence between the divergence field of the derived wind field and the vertical velocity.And,the horizontal and vertical circulations within and around the squall line,as well as strong updrafts,the associated downdrafts,and associated rear inflow of the bow echo,are analyzed well.It is worth mentioning that the variational method in this paper can be applied to simultaneously synthesize the three-dimensional wind field from multiple-Doppler radar observations.
基金Supported by the National Natural Science Foundation of China(11372269,10902057)
文摘The theoretical solutions are obtained for the three-dimensional(3-D)stress field in an infinite isotropic elastic plate with a through-the-thickness circular hole subjected to shear load at far field by using Kane and Mindlin′s assumption based on the stress function method.Based on the present solutions,the characteristics of 3-D stress field are analyzed and the emphasis is placed on the effects of the plate thickness and Poisson′s ratio on the deviation of the present 3-D in-plane stress from the related plane stress solutions,the stress concentration and the out-of-plane constraint.The present solutions show that the stress concentration factor reaches its peak value of about 8.9% which is higher than that of the plane stress solutions.As expected,the out-of-plane stress constraint factor can reach 1on the surface of the hole when the plate is a very thick one.
基金Supported by the National Natural Science Foundation of China(51679080 and 51379073)the Fundamental Research Funds for the Central Universities(B230205020).
文摘This study employed a computational fluid dynamics model with an overset mesh technique to investigate the thrust and power of a floating offshore wind turbine(FOWT)under platform floating motion in the wind–rain field.The impact of rainfall on aerodynamic performance was initially examined using a stationary turbine model in both wind and wind–rain conditions.Subsequently,the study compared the FOWT’s performance under various single degree-of-freedom(DOF)motions,including surge,pitch,heave,and yaw.Finally,the combined effects of wind–rain fields and platform motions involving two DOFs on the FOWT’s aerodynamics were analyzed and compared.The results demonstrate that rain negatively impacts the aerodynamic performance of both the stationary turbines and FOWTs.Pitch-dominated motions,whether involving single or multiple DOFs,caused significant fluctuations in the FOWT aerodynamics.The combination of surge and pitch motions created the most challenging operational environment for the FOWT in all tested scenarios.These findings highlighted the need for stronger construction materials and greater ultimate bearing capacity for FOWTs,as well as the importance of optimizing designs to mitigate excessive pitch and surge.
基金Project supported by the Main Direction Program of Knowledge Innovation of the Chinese Academy of Sciences for Distinguished Young Scholars (Grant No. KZCX2-EW-QN114)the National Natural Science Foundation of China for Young Scholars (GrantNos. 41004006,41131067,and 11173049)+3 种基金the Merit-Based Scientific Research Foundation of the State Ministry of Human Resources and Social Security of China for Returned Overseas Chinese Scholars (Grant No. 2011)the Open Research Fund Programof the Key Laboratory of Computational Geodynamics of the Chinese Academy of Sciences (Grant No. 2011-04)the Frontier Field Program of Knowledge Innovation of Institute of Geodesy and Geophysics of the Chinese Academy of Sciencesthe Open Fund of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Grant No. PLN1113)
文摘The accuracy of the Earth＇s gravitational field measured from the gravity field and steady-state ocean circulation explorer（GOCE）,up to 250 degrees,influenced by the radial gravity gradient V zz and three-dimensional gravity gradient V ij from the satellite gravity gradiometry（SGG） are contrastively demonstrated based on the analytical error model and numerical simulation,respectively.Firstly,the new analytical error model of the cumulative geoid height,influenced by the radial gravity gradient V zz and three-dimensional gravity gradient V ij are established,respectively.In 250 degrees,the GOCE cumulative geoid height error measured by the radial gravity gradient V zz is about 2 1/2 times higher than that measured by the three-dimensional gravity gradient V ij.Secondly,the Earth＇s gravitational field from GOCE completely up to 250 degrees is recovered using the radial gravity gradient V zz and three-dimensional gravity gradient V ij by numerical simulation,respectively.The study results show that when the measurement error of the gravity gradient is 3×10 12 /s 2,the cumulative geoid height errors using the radial gravity gradient V zz and three-dimensional gravity gradient V ij are 12.319 cm and 9.295 cm at 250 degrees,respectively.The accuracy of the cumulative geoid height using the three-dimensional gravity gradient V ij is improved by 30%-40% on average compared with that using the radial gravity gradient V zz in 250 degrees.Finally,by mutual verification of the analytical error model and numerical simulation,the orders of magnitude from the accuracies of the Earth＇s gravitational field recovery make no substantial differences based on the radial and three-dimensional gravity gradients,respectively.Therefore,it is feasible to develop in advance a radial cold-atom interferometric gradiometer with a measurement accuracy of 10 13 /s 2-10 15 /s 2 for precisely producing the next-generation GOCE Follow-On Earth gravity field model with a high spatial resolution.
基金supported by the National Natural Science Foundat ion of China(Grant no.41875045)
文摘The wind and temperature fields at 20 to 55 km above the Antigua launch site(17°N,61°W)were analyzed by using sounding rocket data published by the research organization on Stratosphere-Troposphere Processes and their Role in Climate(SPARC).The results showed distinct variations in the wind and temperature fields at different heights from the 1960s to the 1990s.The overall zonal wind speed showed a significant increasing trend with the year,and the overall change in meridional wind speed showed a falling trend from 1976 to 1990,whereas the temperature field showed a significant cooling trend from 1964 to 1990.The times the trends mutated varied at different levels.By taking the altitudes at 20,35,and 50 km as representatives,we found that the zonal wind speed trend had mutated in 1988,1986,and 1986,respectively;that the meridional wind speed trend had mutated in 1990,1986,and 1990,respectively;and that the temperature trend had mutated separately in 1977,1973,and 1967,respectively.Characteristics of the periodic wind and temperature field variations at different heights were also analyzed,and obvious differences were found in time scale variations across the different layers.The zonal and meridional wind fields were basically characterized as having a significant periodic variation of 5 years across the three layers,and each level was characterized as having a periodic variation of less than 5 years.Temperature field variation at the three levels was basically characterized as occurring in 10-year and 5-year cycles.
基金Supported by the Sinopec Science and Technology Project(P22183).
文摘In the Jiaoshiba block of the Fuling shale gas field,the employed reserves and recovery factor by primary well pattern are low,no obvious barrier is found in the development layer series,and layered development is difficult.Based on the understanding of the main factors controlling shale gas enrichment and high production,the theory and technology of shale gas three-dimensional development,such as fine description and modeling of shale gas reservoir,optimization of three-dimensional development strategy,highly efficient drilling with dense well pattern,precision fracturing and real-time control,are discussed.Three-dimensional development refers to the application of optimal and fast drilling and volume fracturing technologies,depending upon the sedimentary characteristics,reservoir characteristics and sweet spot distribution of shale gas,to form"artificial gas reservoir"in a multidimensional space,so as to maximize the employed reserves,recovery factor and yield rate of shale gas development.In the research on shale gas three-dimensional development,the geological+engineering sweet spot description is fundamental,the collaborative optimization of natural fractures and artificial fractures is critical,and the improvement of speed and efficiency in drilling and fracturing engineering is the guarantee.Through the implementation of three-dimensional development,the overall recovery factor in the Jiaoshiba block has increased from 12.6%to 23.3%,providing an important support for the continuous and stable production of the Fuling shale gas field.
基金supported by the Ministry of Science and Technology of China (2010DFA32680)the National Natural Science Foundation of China (21005062)the Fundamental Research Funds for the Central Universities (CDJRC10220010)
文摘Three-dimensional holographic vector of atomic interaction field（3D-HoVAIF） is used to describe the chemical structures of polychlorinated naphthalenes（PCNs）.After variable screening by stepwise multiple regression（SMR） technique,the liner relationships between gas-chromatographic relative retention time（RRT）,298 K supercooled liquid pressures（logPL）,n-octanol/air partition coefficient（logKOA）,n-octanol/water partition coefficient（logKOW）,aqueous solubilities（logSW）,relative in vitro potency values（-logEROD） of PCNs and 3D-HoVAIF descriptors have been established by partial least-square（PLS） regression.The result shows that the 3D-HoVAIF descriptors can be well used to express the quantitative structure-property（activity） relationships of PCNs.Predictive capability of the models has also been demonstrated by leave-one-out cross-validation.Moreover,the predicted values have been presented for those PCNs which are lack of experimentally physico-chemical properties and biological activity by the optimum models.
基金Project(61120106010)supported by the Projects of Major International(Regional)Joint Research Program Nature Science Foundation of ChinaProject(61304215,61203078)supported by National Natural Science Foundation of China+1 种基金Project(2013000704)supported by the Beijing Outstanding Ph.D.Program Mentor,ChinaProject(61321002)supported by the Foundation for Innovative Research Groups of the National Natural Science Foundation of China
文摘The optimal path planning for fixed-wing unmanned aerial vehicles(UAVs) in multi-target surveillance tasks(MTST) in the presence of wind is concerned.To take into account the minimal turning radius of UAVs,the Dubins model is used to approximate the dynamics of UAVs.Based on the assumption,the path planning problem of UAVs in MTST can be formulated as a Dubins traveling salesman problem(DTSP).By considering its prohibitively high computational cost,the Dubins paths under terminal heading relaxation are introduced,which leads to significant reduction of the optimization scale and difficulty of the whole problem.Meanwhile,in view of the impact of wind on UAVs' paths,the notion of virtual target is proposed.The application of the idea successfully converts the Dubins path planning problem from an initial configuration to a target in wind into a problem of finding the minimal root of a transcendental equation.Then,the Dubins tour is derived by using differential evolution(DE) algorithm which employs random-key encoding technique to optimize the visiting sequence of waypoints.Finally,the effectiveness and efficiency of the proposed algorithm are demonstrated through computational experiments.Numerical results exhibit that the proposed algorithm can produce high quality solutions to the problem.
基金This project is supported by National Natural Science Foundation ofChina(No. 50406017)
文摘An experimental investigation of three-dimensional flow field in a film-cooled turbine model is carried out by using particle image velocimeter （PIV） in a low-speed wind tunnel. The effects of different blowing ratios （M=1.5, 2） on the flow field are studied. The experimental results reveal the classical phenomena of the formation of kidney vortex pair and secondary flow in wake region behind the jet hole. And the changes of the kidney vortex pair and the wake at different locations away from the hole on the suction and pressure sides are also studied. Compared with the flow field in stationary cascade, there are centrifugal force and Coriolis force existing in the flow field of rotating turbine, and these forces bring the radial velocity in the jet flow. The effect of rotatien on the flow field of the pressure side is more distinct than that on the suction side from the measured flow fields in Y-Z plane and radial velocity contours. The increase of blowing ratio makes the kidney vortex pair and the secondary flow in the wake region stronger and makes the range of the wake region enlarged.
基金Supported by the National Key Research and Development Program of China(Nos.2021YFC3101801,2023YFC3008200)the National Natural Science Foundation of China(Nos.42476219,41976200)+6 种基金the National Foreign Experts Program(No.S20240134)the Innovative Team Plan of the Department of Education of Guangdong Province(No.2023KCXTD015)the Tropical Ocean Environment in Western Coastal Waters Observation and Research Station of Guangdong Province(No.2024B1212040008)the Independent Research Project of the Southern Ocean Laboratory(No.SML2022SP301)the Shandong Innovation and Development Research Institute Think Tank Projectthe Guangdong Ocean University Scientific Research Program(No.060302032106)the Start-up Fund for Ph D Researchers(No.060302032104)。
文摘Typhoon Chaba was the most intense typhoon to strike western Guangdong since Typhoon Mujigae in 2015.According to the National Disaster Reduction Center of China,in the morning of July 7,2022,over 1.5 million people in Guangdong,Guangxi,and Hainan were affected by Typhoon Chaba.The typhoon also caused the“Fukui 001”ship to be in distress in the waters near Yangjiang,Guangdong,on July 2,resulting in big casualties.Studies have indicated that wind field forecast for Typhoon Chaba was not accurate.To better simulate typhoon events and assess their impacts,we proposed the use of a model wind field(Fujita-Takahashi)integrated with the Copernicus Marine and Environmental Monitoring Service(CMEMS)data to reconstruct effectively the overall wind field of Typhoon Chaba.The simulation result aligns well with the observations,particularly at the Dashu Island Station,showing consistent trends in wind speed changes.However,certain limitations were noted.The model shows that the attenuation of wind speed is slower when typhoon neared land than that observed,indicating that the model has a high simulation accuracy for the ocean wind field,but may have deviations near coastal areas.The result is accurate for open sea but deviated for near land due to the land friction effect.Therefore,we recommend to adjust the model to improve the accuracy for near coasts.
基金the National Key Research and Development Plan of China[Grant No.2023YFB3002400]the Shanghai 2021 Natural Science Foundation[Grant Nos.21ZR1420400 and 21ZR1419800]+1 种基金the Shanghai 2023 Natural Science Foundation[Grant No.23ZR1463000]the Shandong Provincial Meteorological Bureau Scientific Research Project[Grant No.2023SDBD05].
文摘Accurate short-term forecast of offshore wind fields is still challenging for numerical weather prediction models.Based on three years of 48-hour forecast data from the European Centre for Medium-Range Weather Forecasts Integrated Forecasting System global model(ECMWF-IFS)over 14 offshore weather stations along the coast of Shandong Province,this study introduces a multi-task learning(MTL)model(TabNet-MTL),which significantly improves the forecast bias of near-surface wind direction and speed simultaneously.TabNet-MTL adopts the feature engineering method,utilizes mean square error as the loss function,and employs the 5-fold cross validation method to ensure the generalization ability of the trained model.It demonstrates superior skills in wind field correction across different forecast lead times over all stations compared to its single-task version(TabNet-STL)and three other popular single-task learning models(Random Forest,LightGBM,and XGBoost).Results show that it significantly reduces root mean square error of the ECMWF-IFS wind speed forecast from 2.20 to 1.25 m s−1,and increases the forecast accuracy of wind direction from 50%to 65%.As an explainable deep learning model,the weather stations and long-term temporal statistics of near-surface wind speed are identified as the most influential variables for TabNet-MTL in constructing its feature engineering.
文摘Proper knowledge of the nature of geomagnetic storms and their relationships with the conditions of the space environment at the outer part of the Earth's magnetosphere(bow shock nose) is essential to increase our resilience to space weather disturbances. In this article, we present an analysis of the interplanetary magnetic field(IMF) and solar wind parameters relevant to 100 geomagnetic storms in Solar Cycle 24. We revisit the relationship between the minimum disturbance storm time index(Dst_(min)), the minimum southward IMF(B_(S, min)), the maximum solar wind density(N_(SW, max)) and speed(V_(max)), and the lag time between the extrema(dT(B_(z), N),dT(B_(z), V)). We end with a regression formula that fits the data, with a coefficient of determination of 0.58, a root mean square error of 21.30 nT, and a mean absolute error of 15.87 nT. Even though more complex machine learning models can outperform this model, it serves as a theoretically sensible alternative for understanding and forecasting geomagnetic storms.
基金supported by the National Natural Science Foundation of China(Nos.51775021,52302511)the Fundamental Research Funds for the Central Universities,China(Nos.YWF-23-JC-01,YWF-23-JC-04,YWF-23-JC-09)。
文摘Stratospheric airships are lighter-than-air vehicles capable of continuous flying for months.The energy balance of the airship is the key to long-duration flights.The stratospheric airship is entirely powered by the solar array.It is necessary to accurately predict the output power of the array for any flight state.Because of the uneven solar radiation received by the solar array,the traditional model based on components has a slow simulation speed.In this study,a data-driven surrogate modeling approach for prediction the output power of the solar array is proposed.The surrogate model is trained using the samples obtained from the high-accuracy simulation model.By using the input parameter preprocessor,the accuracy of the surrogate model in predicting the output power of the solar array is improved to 98.65%.In addition,the predictive speed of the surrogate model is ten million times faster than the traditional simulation model.Finally,the surrogate model is used to predict the energy balance of stratospheric airships flying throughout the year under actual global wind fields.
基金National Natural Science Foundation of China under Grant Nos.51278160,51478155,51378147
文摘Stochastic simulation is an important means of acquiring fluctuating wind pressures for wind induced response analyses in structural engineering. The wind pressure acting on a large-span space structure can be characterized as a stationary non-Gaussian field. This paper reviews several simulation algorithms related to the Spectral Representation Method （SRM） and the Static Transformation Method （STM）. Polynomial and Exponential transformation functions （PSTM and ESTM） are discussed. Deficiencies in current algorithms, with respect to accuracy, stability and efficiency, are analyzed, and the algorithms are improved for better practical application. In order to verify the improved algorithm, wind pressure fields on a large-span roof are simulated and compared with wind tunnel data. The simulation results fit well with the wind tunnel data, and the algorithm accuracy, stability and efficiency are shown to be better than those of current algorithms.
文摘The core sampling experiments were conducted after hydraulic fracturing in the three-dimensional development zone of Fuling shale gas.Six coring wells of different well types were systematically designed.Based on the integrated engineering technology of post-fracturing drilling,coring and monitoring of shale and the analysis of fracture source tracing,the evaluation of the fracture network after fracturing in the three-dimensional development of shale gas was conducted.The data of core fractures after fracturing indicate that three major types of fractures are formed after fracturing:natural fractures,hydraulic fractures,and fractures induced by external mechanical force,which are further classified into six subcategories:natural structural fractures,natural bedding fractures,hydraulic fractures,hydraulically activated fractures,drilling induced fractures,and fractures induced by core transportation.The forms of the artificial fracture network after fracturing are complex.Hydraulic fractures and hydraulically activated fractures interweave with each other,presenting eight forms of artificial fracture networks,among which the“一”-shaped fracture is the most common,accounting for approximately 70%of the total fractures.When the distance to the fractured wellbore is less than 35 m,the density of the artificial fracture network is relatively high;when it is 35–100 m,the density is lower;and when it is beyond 100 m,the density gradually increases.The results of the fracture tracing in the core sampling area confirm that the current fracturing technology can essentially achieve the differential transformation of the reservoir in the main area of Jiaoshiba block in Fuling.The three-layer three-dimensional development model can efficiently utilize shale gas reserves,although there is still room for improvement in the complexity and propagation uniformity of fractures.It is necessary to further optimize technologies such as close-cutting combined with temporary blocking and deflection within fractures or at fracture mouths,as well as limited flow perforation techniques,to promote the balanced initiation and extension of fractures.
基金the support of the National Natural Science Foundation of China(Nos.41974073,41404053)the Macao Foundation and the pre-research project of Civil Aerospace Technologies(Nos.D020308 and D020303)+2 种基金funded by the National Space Administration of Chinathe opening fund of the State Key Laboratory of Lunar and Planetary Sciences(Macao University of Science and Technology,Macao Science and Technology Development Fund No.119/2017/A3)the Specialized Research Fund for State Key Laboratories,and the NUIST-UoR International Research Institute。
文摘We combined domestic ground-based and satellite magnetic measurements to create a regional three-dimensional surface Spline(3DSS)gradient model of the main geomagnetic field over the Chinese continent.To improve the precision of the model,we considered the data gap between the ground and satellite data.We compared and analyzed the results of the Taylor polynomial,surface Spline,and CHAOS-6(the CHAMP,?rsted and SAC-C model of Earth’s magnetic field)gradient models.Results showed that the gradients in the south-north and east-west directions of the four models were consistent.The 3DSS model was able to express not only gradients at different altitudes,but also average gradients inside the research area.The two Spline models were able to capture more information on gradient anomalies than were the fitted models.Strong local anomalies were observed in northern Xinjiang,Beijing,and the junction area between Jiangsu and Zhejiang,and the total intensity F decreased whereas the altitude increased.The gradient decreased by 21.69%in the south-north direction and increased by 11.78%in the east-west direction.In addition,the altitude gradient turned from negative to positive while the altitude increased.The Spline model and the two fitted models differed mainly in the field sources they expressed and the modeling theory.
