[Objective] This study was to explore the difference of kriging interpolation and sequential Gaussian simulation on analyzing soil heavy metal pollution with a view to provide references for analyzing the heavy metal ...[Objective] This study was to explore the difference of kriging interpolation and sequential Gaussian simulation on analyzing soil heavy metal pollution with a view to provide references for analyzing the heavy metal pollution of soil. [Method] The sampling data of soil copper from a county of Liaocheng, Shandong Province was set as the study objective. Kriging interpolation and sequential Gaussian simu- lation were used to simulate the spatial distribution of soil copper. And 30 sampling points were selected as the cross-validation data set to compare the two interpola- tion methods. [Result] Kriging method and Gaussian sequential simulation have their own advantages on simulating mean segment and extreme segment, therefore, re- searchers should choose the proper method based on the characteristics of test data and application purposes. [Conclusion] Analysis of soil heavy metal pollution is the prerequisite for soil management and ecological restoration. The result of this study is of important significance for choosing different interpolating and simulating methods to analyze soil heavy metal pollution based on different purposes.展开更多
In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This st...In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This study combines the DIE method with the moving Kriging interpolation to present a boundary-type meshfree method, and the corresponding formulae of the MKIBNM are derived. In the present method, the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker＇s delta property, then the boundary conditions can be imposed directly and easily. To verify the accuracy and stability of the present formulation, three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.展开更多
Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the movi...Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.展开更多
The accuracy of spatial interpolation of precipitation data is determined by the actual spatial variability of the precipitation, the interpolation method, and the distribution of observatories whose selections are pa...The accuracy of spatial interpolation of precipitation data is determined by the actual spatial variability of the precipitation, the interpolation method, and the distribution of observatories whose selections are particularly important. In this paper, three spatial sampling programs, including spatial random sampling, spatial stratified sampling, and spatial sandwich sampling, are used to analyze the data from meteorological stations of northwestern China. We compared the accuracy of ordinary Kriging interpolation methods on the basis of the sampling results. The error values of the regional annual pre-cipitation interpolation based on spatial sandwich sampling, including ME （0.1513）, RMSE （95.91）, ASE （101.84）, MSE （？0.0036）, and RMSSE （1.0397）, were optimal under the premise of abundant prior knowledge. The result of spatial stratified sampling was poor, and spatial random sampling was even worse. Spatial sandwich sampling was the best sampling method, which minimized the error of regional precipitation estimation. It had a higher degree of accuracy compared with the other two methods and a wider scope of application.展开更多
Quality-controlled and serially complete daily air temperature data are essential to evaluating and modelling the influences of climate change on the permafrost in cold regions. Due to malfunctions and location chang...Quality-controlled and serially complete daily air temperature data are essential to evaluating and modelling the influences of climate change on the permafrost in cold regions. Due to malfunctions and location changes of observing stations, temporal gaps （i.e., missing data） are common in collected datasets. The objective of this study was to assess the efficacy of Kriging spatial interpolation for estimating missing data to fill the temporal gaps in daily air temperature data in northeast China. A cross-validation experiment was conducted. Daily air temperature series from 1960 to 2012 at each station were estimated by using the universal Kriging （UK） and Kriging with an external drift （KED）, as appropriate, as if all the ob-servations at a given station were completely missing. The temporal and spatial variation patterns of estimation uncertainties were also checked. Results showed that Kriging spatial interpolation was generally desirable for estimating missing data in daily air temperature, and in this study KED performed slightly better than UK. At most stations the correlation coefficients （R2） between the observed and estimated daily series were 〉0.98, and root mean square errors （RMSEs） of the estimated daily mean （Tmean）, maximum （Tmax）, and minimum （Tmin） of air temperature were 〈3 ℃. However, the estimation quality was strongly affected by seasonality and had spatial variation. In general, estimation uncertainties were small in summer and large in winter. On average, the RMSE in winter was approximately 1 ℃ higher than that in summer. In addition, estimation uncertainties in mountainous areas with complex terrain were significantly larger than those in plain areas.展开更多
The analysis of grey system, kriging interpolation, and integration selection index were employed to investigate the relationships between the flower yield/plant （FY） and 15 other quantitative traits of 20 rugosa ro...The analysis of grey system, kriging interpolation, and integration selection index were employed to investigate the relationships between the flower yield/plant （FY） and 15 other quantitative traits of 20 rugosa rose cultivars. The result showed that： The grey relational grade （GRG） of the number of flowers/plant （NF）, the number of branches/plant （NB）, the width of floral bud （WB）, and the weight/flower （WF） to the FY were larger （〉 0.5）; FY improved with the increase of NF and NB. Moreover, the indirect selection of either trait could not achieve improvement of FY. It is necessary to improve FY by multi-trait selection. The integration selection index （ISI） equation of FY was established with the characters NF, NB, WB, and WF： I= 0.3187x1 - 318.6x2 ＋ 670.1 x4 ＋ 6.3xa, index heritability = 0.8014, selective response of the integration breeding value = 245.8811. This will provide a theoretic basis for the genetic breeding of rugosa rose.展开更多
The shear behavior of large-scale weak intercalation shear zones(WISZs)often governs the stability of foundations,rock slopes,and underground structures.However,due to their wide distribution,undulating morphology,com...The shear behavior of large-scale weak intercalation shear zones(WISZs)often governs the stability of foundations,rock slopes,and underground structures.However,due to their wide distribution,undulating morphology,complex fabrics,and varying degrees of contact states,characterizing the shear behavior of natural and complex large-scale WISZs precisely is challenging.This study proposes an analytical method to address this issue,based on geological fieldwork and relevant experimental results.The analytical method utilizes the random field theory and Kriging interpolation technique to simplify the spatial uncertainties of the structural and fabric features for WISZs into the spatial correlation and variability of their mechanical parameters.The Kriging conditional random field of the friction angle of WISZs is embedded in the discrete element software 3DEC,enabling activation analysis of WISZ C2 in the underground caverns of the Baihetan hydropower station.The results indicate that the activation scope of WISZ C2 induced by the excavation of underground caverns is approximately 0.5e1 times the main powerhouse span,showing local activation.Furthermore,the overall safety factor of WISZ C2 follows a normal distribution with an average value of 3.697.展开更多
Formation pressure is the key parameter for the analysis of wellbore safety.With increasing drilling depth,how-ever,the behavior of this variable becomes increasingly complex.In this work,a 3D model of the formation p...Formation pressure is the key parameter for the analysis of wellbore safety.With increasing drilling depth,how-ever,the behavior of this variable becomes increasingly complex.In this work,a 3D model of the formation pres-sure under uncertainty is presented.Moreover a relevant algorithm is elaborated.First,the logging data of regional key drilling wells are collected and a one-dimensional formation pressure profile along the well depth is determined.Then,a 3D model of regional formation pressure of the hierarchical group layer is defined by using the Kriging interpolation algorithm relying on a support vector machine(SVM)and the formation pressure of the drilled wells.To validate the method,the formation pressure of one pre-drilled well is compared with the well logging results.The comparison reveals that the maximum relative error is less than 4.5%.The software based on this model is complemented by a computer visualization technology,which provides a relevant tool for under-standing and analyzing the 3D formation pressure.The outcomes of this study are intended to support the char-acterization of areas with missing or poor 3D seismic data and provide more accurate information for the analysis of wellbore integrity.展开更多
Located in Lu-Zong ore concentration area, middle-lower Yangtze metallogenic belt, ShaXi porphyry copper deposit is a typical hydrothermal deposit. To investigate the distribution of deep ore bodies and spatial charac...Located in Lu-Zong ore concentration area, middle-lower Yangtze metallogenic belt, ShaXi porphyry copper deposit is a typical hydrothermal deposit. To investigate the distribution of deep ore bodies and spatial characteristics of host structures, an AMT survey was conducted in mining area. Eighteen pseudo-2D resistivity sections were constructed through careful processing and inversion. These sections clearly show resistivity difference between the Silurian sandstones formation and quartz diorite porphyry and this porphyry copper formation was controlled by the highly resistive anticlines. Using 3D block Kxiging interpolation method and 3D visualization techniques, we constructed a detailed 3D resistivity model of quartz diorite porphyry which shows the shape and spatial distribution of deep ore bodies. This case study can serve as a good example for future ore prospecting in and around this mining area.展开更多
Geometric anisotropy is commonly assumed in the investigation of the spatial variations of geophysical parameters. However, this assumption is not always satisfied in practice. We propose a novel method to determine t...Geometric anisotropy is commonly assumed in the investigation of the spatial variations of geophysical parameters. However, this assumption is not always satisfied in practice. We propose a novel method to determine the anisotropy of geophysical parameters. In the proposed method, the variograms are first normalized in all directions. Then, the normalized samples are fitted by the unit range variation increment（URVI） function to estimate the intensities of the variograms in each direction, from which the anisotropy can be finally determined. The performance of the proposed method is validated using InSAR atmospheric delay measurements over the Shanghai region. The results show that the deviation of the method is 6.4%, and that of the geometric anisotropy-based method is 21.2%. In addition, the computational efficiency of the new method is much higher. Subsequently, the URVI- and the geometric anisotropy-based methods are cross-validated in the cross-validation experiments by using Kriging interpolation. The results demonstrate that the structure functions generated with the proposed method are more accurate and can better refl ect the spatial characteristics of the random fi eld. Therefore, the proposed method, which is more accurate and effi cient to determine the anisotropy than the conventional geometry anisotropy-based method, provides a better foundation to estimate the geophysical parameters of interest.展开更多
This research used both geostatistics and GIS approach to compare temporal change of soil salt between 1980 and 2003, to analyze the spatial distribution of surface soil salt, to developed methods for predicting soil ...This research used both geostatistics and GIS approach to compare temporal change of soil salt between 1980 and 2003, to analyze the spatial distribution of surface soil salt, to developed methods for predicting soil salinization potential based on recent improvements to the Dempster-Shafer theory, and to develop probability maps of potential salinization in Yucheng City, China. A semivariogram model of soil salt content was developed from the spherical model, and then employing kriging interpolation the spatial distribution of salt content in 2003 was obtained utilizing data from 100 soil sampling points. Potential salinization distribution was mapped using an approach that integrated soil data of the second general survey in 1980 in Yucheng City, which included groundwater salinity, groundwater depth, soil texture, soil organic matter content, and geomorphic maps. With the support of Dempster-Shafer theory and fuzzy set technique the factors that affected potential soil salinization were characterized and integrated; and then soil salinization was predicted. Finally a prognosis map of potential salinization distribution in the research area was obtained, with higher probability values indicating higher hazards to salinity processes. The distribution of the potential soil salinization probability was a successive surface.展开更多
Sand and dust storms （SDS） are common phenomena in arid and semi-arid areas. In recent years, SDS frequencies and intensities have increased significantly in Iran. A research on SDS sources is important for understa...Sand and dust storms （SDS） are common phenomena in arid and semi-arid areas. In recent years, SDS frequencies and intensities have increased significantly in Iran. A research on SDS sources is important for understanding the mechanisms of dust generation and assessing its socio-economic and environmental impacts. In this paper, we developed a new approach to identify SDS source areas in Iran using a combination of nine related datasets, namely drought events, temperature, precipitation, location of sandy soils, SDS frequency, hu- man-induced soil degradation （HISD）, human influence index （HII）, rain use efficiency （RUE） and net primary pro- ductivity （NPP） loss. To identify SDS source areas, we firstly normalized these datasets under uniform criteria in- cluding layer reprojection using Lambert conformal conic projection, data conversion from shapefile to raster, Min-Max Normalization with data range from 0 to 1, and data interpolation by Kriging and images resampling （resolution of 1 km）. After that, a score map for the possibility of SDS sources was generated through overlaying multiple datasets under average weight allocation criterion, in which each item obtained weight equally. In the score map, the higher the score, the more possible a specific area could be regarded as SDS source area. Exceptions mostly came from large cities, like Tehran and Isfahan. As a result, final SDS source areas were mapped out, and Al-Howizeh/Al-Azim marshes and Sistan Basin were identified as main SDS source areas in Iran. The SDS source area in Al-Howizeh/Al-Azim marshes still keeps expanding. In addition, Al-Howizeh/Al-Azim marshes are now suf- fering rapid land degradation due to natural and human-induced factors and might totally vanish in the near future. Sistan Basin also demonstrates the impacts of soil degradation and wind erosion. With appropriate intensity, dura- tion, wind speed and altitude of the dust storms, sand particles uplifting from this area might have developed into extreme dust storms, especially during the summer.展开更多
Rapid urbanization results in the conversion of natural soil to urban soil,and consequently,the storage and density of the soil carbon pools change.Taking Chongqing Municipality of China as a study case,this investiga...Rapid urbanization results in the conversion of natural soil to urban soil,and consequently,the storage and density of the soil carbon pools change.Taking Chongqing Municipality of China as a study case,this investigation attempts to better understand soil carbon pools in hilly cities.First,the vegetated areas in the study area were derived from QuickBird images.Then,topsoil data from 220 soil samples(0-20 cm) in the vegetated areas were collected and their soil organic carbon(SOC) densities were analyzed.Using the Kriging interpolation method,the spatial pattern of SOC was estimated.The results show that the SOC density exhibited high spatial variability in the urban topsoil of Chongqing.First,the SOC density in topsoil decreased according to slope in the order 2°-6° < 25°-90° < 0°-2° < 6°-15° < 15°-25°.Second,the newly developed areas during 2001-2010 had a lower SOC density than the areas built before 1988.Third,urban parks and gardens had a higher SOC density in topsoil,residential green land followed,and scattered street green land ranked last.For hilly cities,the variability of terrain affects the distribution of SOC.The Kriging results indicate that Kriging method combining slope with SOC density produced a high level of accuracy.The Kriging results show that the SOC density to the north of the Jialing River was higher than the south.The vegetated areas were estimated to amount to 73.5 km2 across the study area with an SOC storage of 0.192 Tg and an average density of 2.61 kg/m2.展开更多
The Thornthwaite moisture index, an index of the supply of water （precipitation） in an area relative to the climatic demand for water （potential evapotranspiration）, was used to examine the spatial and temporal va...The Thornthwaite moisture index, an index of the supply of water （precipitation） in an area relative to the climatic demand for water （potential evapotranspiration）, was used to examine the spatial and temporal variation of drought and to verify the influence of environmental factors on the drought in the Hengduan Mountains, China. Results indicate that the Thornthwaite moisture index in the Hengduan Mountains had been increasing since 1960 with a rate of 0.1938/yr. Annual Thomthwaite moisture index in Hengduan Mountains was between -97.47 and 67.43 and the spatial heterogeneity was obvious in different seasons. Thomthwaite moisture index was high in the north and low in the south, and the monsoon rainfall had a significant impact on its spatial distribution. The tendency rate of Thomthwaite moisture index variation varied in different seasons, and the increasing trends in spring were greater than that in summer and autumn. However, the Thomthwaite moisture index decreased in winter. Thomthwaite moisture index increased greatly in the north and there was a small growth in the south of Hengduan Mountains. The increase of precipitation and decrease of evaporation lead to the increase of Thomthwaite moisture index. Thornthwaite moisture index has strong correlation with vegetation coverage. It can be seen that the correlation between Normalized Difference Vegetation Index （NDVI） and Thomthwaite moisture index was positive in spring and summer, but negative in autumn and winter. Correlation between Thornthwaite moisture index and relative soil relative moisture content was positive in spring, summer and autumn, but negative in winter. The typical mountainous terrain affect the distribu- tion of temperature, precipitation, wind speed and other meteorological factors in this region, and then affect the spatial distribution of Thomthwaite moisture index. The unique ridge-gorge terrain caused the continuity of water-heat distribution from the north to south, and the water-heat was stronger than that from the east to west part, and thus determined the spatial distribution of Thornthwaite mois- ture index. The drought in the Hengduan Mountains area is mainly due to the unstable South Asian monsoon rainfall time.展开更多
The acquisition of precise soil data representative of the entire survey area, is a critical issue for many treatments such as irrigation or fertilization in precision agriculture. The aim of this study was to investi...The acquisition of precise soil data representative of the entire survey area, is a critical issue for many treatments such as irrigation or fertilization in precision agriculture. The aim of this study was to investigate the spatial variability of soil bulk electrical conductivity （ECb） in a coastal saline field and design an optimized spatial sampling scheme of ECb based on a sampling design algorithm, the variance quad-tree （VQT） method. Soil ECb data were collected from the field at 20 m interval in a regular grid scheme. The smooth contour map of the whole field was obtained by ordinary kriging interpolation, VQT algorithm was then used to split the smooth contour map into strata of different number desired, the sampling locations can be selected within each stratum in subsequent sampling. The result indicated that the probability of choosing representative sampling sites was increased significantly by using VQT method with the sampling number being greatly reduced compared to grid sampling design while retaining the same prediction accuracy. The advantage of the VQT method is that this scheme samples sparsely in fields where the spatial variability is relatively uniform and more intensive where the variability is large. Thus the sampling efficiency can be improved, hence facilitate an assessment methodology that can be applied in a rapid, practical and cost-effective manner.展开更多
A boundary-type meshless method called the scaled boundary node method （SBNM） is developed to directly evaluate mixed mode stress intensity factors （SIFs） without extra post-processing. The SBNM combines the scale...A boundary-type meshless method called the scaled boundary node method （SBNM） is developed to directly evaluate mixed mode stress intensity factors （SIFs） without extra post-processing. The SBNM combines the scaled boundary equations with the moving Kriging （MK） interpolation to retain the dimensionality advantage of the former and the meshless attribute of the latter. As a result, the SBNM requires only a set of scattered nodes on the boundary, and the displacement field is approximated by using the MK interpolation technique, which possesses the 5 function property. This makes the developed method efficient and straightforward in imposing the essential boundary conditions, and no special treatment techniques are required. Besides, the SBNM works by weakening the governing differential equations in the circumferential direction and then solving the weakened equations analytically in the radial direction. Therefore, the SBNM permits an accurate representation of the singularities in the radial direction when the scaling center is located at the crack tip. Numerical examples using the SBNM for computing the SIFs are presented. Good agreements with available results in the literature are obtained.展开更多
[Objective] The aim was to study the spatial information of temperature and precipitation data in Hengduan mountains. [Method] Considering GIS spatial interpolation and numerical statistics theory, spatial prediction ...[Objective] The aim was to study the spatial information of temperature and precipitation data in Hengduan mountains. [Method] Considering GIS spatial interpolation and numerical statistics theory, spatial prediction were carried out to the ten years average temperature and precipitation flux observation data in 109 sparse meteorological stations in Hengduan Mountains. Based on the spatial range of geographic position of Hengduan Mountains, and 1∶1 000 000 scale DEM as data sources, and using trend surface simulation and residual ordinary Kriging interpolation correction method, the spatial continuous surface for annual average temperature and precipitation were studied. [Result] It was scientific and reasonable to use certain unevenly distributed sparse climate observation station value, and by dint of trend simulation and residue interpolation method to get climate consecutive spatial result. This method can not only accurate the temperature and precipitation spatial distributions to grid point, but also can reflect macro and micro geography factors and topographic influence factor of variation. Furthermore, it can be predicted error term trend surface reasonable spatial distribution. Simulation results were basically in accordance with the objective law, and can be used for the region climate data spatial informatization simulation. [Conclusion] The study provided scientific spatial basic data for the further study of ecological and vegetation in Hengduan Mountains.展开更多
In abandoned mine sites, i.e., mine sites where mining operations have ended, wide spread contaminations are often evident, but the potential sources and pathways of contamination especially through the subsurface, ar...In abandoned mine sites, i.e., mine sites where mining operations have ended, wide spread contaminations are often evident, but the potential sources and pathways of contamination especially through the subsurface, are difficult to identify due to inadequate and sparse geochemical measurements available. Therefore, it is essential to design and implement a planned monitoring net-work to obtain essential information required for establishing the potential contamination source locations, i.e., waste dumps, tailing dams, pits and possible pathways through the subsurface, and to design a remediation strategy for rehabilitation. This study presents an illustrative application of modeling the flow and transport processes and monitoring network design in a study area hydrogeologically resembling an abandoned mine site in Queensland, Australia. In this preliminary study, the contaminant transport process modeled does not incorporate the reactive geochemistry of the contaminants. The transport process is modeled considering a generic conservative contaminant for the illustrative purpose of showing the potential application of an optimal monitoring design methodology. This study aims to design optimal monitoring network to: 1) minimize the contaminant solute mass estimation error;2) locate the plume boundary;3) select the monitoring locations with (potentially) high concentrations. A linked simulation optimization based methodology is utilized for optimal monitoring network design. The methodology is applied utilizing a recently developed software package CARE-GWMND, developed at James Cook University for optimal monitoring network design. Given the complexity of the groundwater systems and the sparsity of pollutant concentration observation data from the field, this software is capable of simulating the groundwater flow and solute transport with spatial interpolation of data from a sparse set of available data, and it utilizes the optimization algorithm to determine optimum locations for implementing monitoring wells.展开更多
Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes w...Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes with- out connectivity are divided into several subsets by cells with arbitrary shape. Local discrete equation is established over each cell by using moving Kriging interpolation, in which the nodes that located in the cell are used for approxima- tion. Then local discrete equations can be simplified by con- densation of degree of freedom, which transfers equations of inner nodes to equations of boundary nodes based on cells. The global dynamic system equations are obtained by as- sembling all local discrete equations and are solved by using the standard implicit Newmark＇s time integration scheme. In the scheme of present method, the calculation of each cell is carried out by meshfree method, and local search is imple- mented in interpolation. Numerical examples show that the present method has high computational efficiency and good accuracy in solving elastic dynamic problems.展开更多
Real-time performance and accuracy are two most challenging requirements in virtual surgery training.These difficulties limit the promotion of advanced models in virtual surgery,including many geometric and physical m...Real-time performance and accuracy are two most challenging requirements in virtual surgery training.These difficulties limit the promotion of advanced models in virtual surgery,including many geometric and physical models.This paper proposes a physical model of virtual soft tissue,which is a twist model based on the Kriging interpolation and membrane analogy.The proposed model can quickly locate spatial position through Kriging interpolation method and accurately compute the force change on the soft tissue through membrane analogy method.The virtual surgery simulation system is built with a PHANTOM OMNI haptic interaction device to simulate the torsion of virtual stomach and arm,and further verifies the real-time performance and simulation accuracy of the proposed model.The experimental results show that the proposed soft tissue model has high speed and accuracy,realistic deformation,and reliable haptic feedback.展开更多
基金Supported by the Science and Technology Development Program of Shandong Province (Soft Science) (2009RKB220),China~~
文摘[Objective] This study was to explore the difference of kriging interpolation and sequential Gaussian simulation on analyzing soil heavy metal pollution with a view to provide references for analyzing the heavy metal pollution of soil. [Method] The sampling data of soil copper from a county of Liaocheng, Shandong Province was set as the study objective. Kriging interpolation and sequential Gaussian simu- lation were used to simulate the spatial distribution of soil copper. And 30 sampling points were selected as the cross-validation data set to compare the two interpola- tion methods. [Result] Kriging method and Gaussian sequential simulation have their own advantages on simulating mean segment and extreme segment, therefore, re- searchers should choose the proper method based on the characteristics of test data and application purposes. [Conclusion] Analysis of soil heavy metal pollution is the prerequisite for soil management and ecological restoration. The result of this study is of important significance for choosing different interpolating and simulating methods to analyze soil heavy metal pollution based on different purposes.
基金Project supported by the Young Scientists Fund of the National Natural Science Foundation of China(Grant No.10902076)the Natural Science Foundation of Shanxi Province of China(Grant No.2007011009)+1 种基金the Scientific Research and Development Program of the Shanxi Higher Education Institutions(Grant No.20091131)the Doctoral Startup Foundation of Taiyuan University of Science and Technology(Grant No.200708)
文摘In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This study combines the DIE method with the moving Kriging interpolation to present a boundary-type meshfree method, and the corresponding formulae of the MKIBNM are derived. In the present method, the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker＇s delta property, then the boundary conditions can be imposed directly and easily. To verify the accuracy and stability of the present formulation, three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.
基金Project supported by the National Natural Science Foundation of China(Grant No.11072117)the Natural Science Foundation of Ningbo City,China(GrantNo.2013A610103)+2 种基金the Natural Science Foundation of Zhejiang Province,China(Grant No.Y6090131)the Disciplinary Project of Ningbo City,China(GrantNo.SZXL1067)the K.C.Wong Magna Fund in Ningbo University,China
文摘Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.
基金conducted within the National Major Scientific Research Project (No. 2013CBA01806)the National Natural Science Foundation of China (No. 41271085)the National Scientific and Technological Support Project (No. 2013BAB05B03)
文摘The accuracy of spatial interpolation of precipitation data is determined by the actual spatial variability of the precipitation, the interpolation method, and the distribution of observatories whose selections are particularly important. In this paper, three spatial sampling programs, including spatial random sampling, spatial stratified sampling, and spatial sandwich sampling, are used to analyze the data from meteorological stations of northwestern China. We compared the accuracy of ordinary Kriging interpolation methods on the basis of the sampling results. The error values of the regional annual pre-cipitation interpolation based on spatial sandwich sampling, including ME （0.1513）, RMSE （95.91）, ASE （101.84）, MSE （？0.0036）, and RMSSE （1.0397）, were optimal under the premise of abundant prior knowledge. The result of spatial stratified sampling was poor, and spatial random sampling was even worse. Spatial sandwich sampling was the best sampling method, which minimized the error of regional precipitation estimation. It had a higher degree of accuracy compared with the other two methods and a wider scope of application.
基金funded by the Chinese National Fund Projects (Nos. 41401028, 41201066)by the State Key Laboratory of Frozen Soils Engineering (Project No. SKLFSE201201)
文摘Quality-controlled and serially complete daily air temperature data are essential to evaluating and modelling the influences of climate change on the permafrost in cold regions. Due to malfunctions and location changes of observing stations, temporal gaps （i.e., missing data） are common in collected datasets. The objective of this study was to assess the efficacy of Kriging spatial interpolation for estimating missing data to fill the temporal gaps in daily air temperature data in northeast China. A cross-validation experiment was conducted. Daily air temperature series from 1960 to 2012 at each station were estimated by using the universal Kriging （UK） and Kriging with an external drift （KED）, as appropriate, as if all the ob-servations at a given station were completely missing. The temporal and spatial variation patterns of estimation uncertainties were also checked. Results showed that Kriging spatial interpolation was generally desirable for estimating missing data in daily air temperature, and in this study KED performed slightly better than UK. At most stations the correlation coefficients （R2） between the observed and estimated daily series were 〉0.98, and root mean square errors （RMSEs） of the estimated daily mean （Tmean）, maximum （Tmax）, and minimum （Tmin） of air temperature were 〈3 ℃. However, the estimation quality was strongly affected by seasonality and had spatial variation. In general, estimation uncertainties were small in summer and large in winter. On average, the RMSE in winter was approximately 1 ℃ higher than that in summer. In addition, estimation uncertainties in mountainous areas with complex terrain were significantly larger than those in plain areas.
文摘The analysis of grey system, kriging interpolation, and integration selection index were employed to investigate the relationships between the flower yield/plant （FY） and 15 other quantitative traits of 20 rugosa rose cultivars. The result showed that： The grey relational grade （GRG） of the number of flowers/plant （NF）, the number of branches/plant （NB）, the width of floral bud （WB）, and the weight/flower （WF） to the FY were larger （〉 0.5）; FY improved with the increase of NF and NB. Moreover, the indirect selection of either trait could not achieve improvement of FY. It is necessary to improve FY by multi-trait selection. The integration selection index （ISI） equation of FY was established with the characters NF, NB, WB, and WF： I= 0.3187x1 - 318.6x2 ＋ 670.1 x4 ＋ 6.3xa, index heritability = 0.8014, selective response of the integration breeding value = 245.8811. This will provide a theoretic basis for the genetic breeding of rugosa rose.
基金support from the Key Projects of the Yalong River Joint Fund of the National Natural Science Foundation of China(Grant No.U1865203)the Innovation Team of Changjiang River Scientific Research Institute(Grant Nos.CKSF2021715/YT and CKSF2023305/YT)。
文摘The shear behavior of large-scale weak intercalation shear zones(WISZs)often governs the stability of foundations,rock slopes,and underground structures.However,due to their wide distribution,undulating morphology,complex fabrics,and varying degrees of contact states,characterizing the shear behavior of natural and complex large-scale WISZs precisely is challenging.This study proposes an analytical method to address this issue,based on geological fieldwork and relevant experimental results.The analytical method utilizes the random field theory and Kriging interpolation technique to simplify the spatial uncertainties of the structural and fabric features for WISZs into the spatial correlation and variability of their mechanical parameters.The Kriging conditional random field of the friction angle of WISZs is embedded in the discrete element software 3DEC,enabling activation analysis of WISZ C2 in the underground caverns of the Baihetan hydropower station.The results indicate that the activation scope of WISZ C2 induced by the excavation of underground caverns is approximately 0.5e1 times the main powerhouse span,showing local activation.Furthermore,the overall safety factor of WISZ C2 follows a normal distribution with an average value of 3.697.
基金supported by Scientific Research and Technology Development Project of CNPC“Study on Exploration and Development Theory and Key Technology of Gulong Shale Oil in Daqing”(2021ZZ10-03)Scientific Research and Technology Development Project of CNPC“Development of Integrated Software(Smart Drilling)for Drilling and Completion Engineering Design and Optimization Decision”(2020B-4019)+1 种基金Scientific Research and Technology Development Project of CNPC“Integration and Experiment of Safe,Optimal and Fast Drilling and Completion Technology for Complex Ultra Deep Wells”(2020F-46)project funded by China Postdoctoral Science Foundation“Research on the Effect of Stress Distribution Difference on Acoustic Propagation Characteristics in Drill String”(2021M693508).
文摘Formation pressure is the key parameter for the analysis of wellbore safety.With increasing drilling depth,how-ever,the behavior of this variable becomes increasingly complex.In this work,a 3D model of the formation pres-sure under uncertainty is presented.Moreover a relevant algorithm is elaborated.First,the logging data of regional key drilling wells are collected and a one-dimensional formation pressure profile along the well depth is determined.Then,a 3D model of regional formation pressure of the hierarchical group layer is defined by using the Kriging interpolation algorithm relying on a support vector machine(SVM)and the formation pressure of the drilled wells.To validate the method,the formation pressure of one pre-drilled well is compared with the well logging results.The comparison reveals that the maximum relative error is less than 4.5%.The software based on this model is complemented by a computer visualization technology,which provides a relevant tool for under-standing and analyzing the 3D formation pressure.The outcomes of this study are intended to support the char-acterization of areas with missing or poor 3D seismic data and provide more accurate information for the analysis of wellbore integrity.
基金supported jointly by the National Natural Science Foundation Fund of China(Grant No.40930418)ChineseGovernment-funded Scientific Programmed of SinoProbe Deep Exploration in China(SinoProbe-03)the Basic Scientific Research-fund of Institute of Mineral Resources,Chinese Academy of Geological Sciences(Grant No.K1008)
文摘Located in Lu-Zong ore concentration area, middle-lower Yangtze metallogenic belt, ShaXi porphyry copper deposit is a typical hydrothermal deposit. To investigate the distribution of deep ore bodies and spatial characteristics of host structures, an AMT survey was conducted in mining area. Eighteen pseudo-2D resistivity sections were constructed through careful processing and inversion. These sections clearly show resistivity difference between the Silurian sandstones formation and quartz diorite porphyry and this porphyry copper formation was controlled by the highly resistive anticlines. Using 3D block Kxiging interpolation method and 3D visualization techniques, we constructed a detailed 3D resistivity model of quartz diorite porphyry which shows the shape and spatial distribution of deep ore bodies. This case study can serve as a good example for future ore prospecting in and around this mining area.
基金sponsored jointly by the National Hi-tech Research and Development Program of China(No.2012AA121301)National Basic Research Program of China(No.2012CB719903)+1 种基金the National Natural Science Foundation of China(Nos.41222027,41474007,and 41404013)Hunan Provincial Natural Science Foundation of China(No.13JJ1006)
文摘Geometric anisotropy is commonly assumed in the investigation of the spatial variations of geophysical parameters. However, this assumption is not always satisfied in practice. We propose a novel method to determine the anisotropy of geophysical parameters. In the proposed method, the variograms are first normalized in all directions. Then, the normalized samples are fitted by the unit range variation increment（URVI） function to estimate the intensities of the variograms in each direction, from which the anisotropy can be finally determined. The performance of the proposed method is validated using InSAR atmospheric delay measurements over the Shanghai region. The results show that the deviation of the method is 6.4%, and that of the geometric anisotropy-based method is 21.2%. In addition, the computational efficiency of the new method is much higher. Subsequently, the URVI- and the geometric anisotropy-based methods are cross-validated in the cross-validation experiments by using Kriging interpolation. The results demonstrate that the structure functions generated with the proposed method are more accurate and can better refl ect the spatial characteristics of the random fi eld. Therefore, the proposed method, which is more accurate and effi cient to determine the anisotropy than the conventional geometry anisotropy-based method, provides a better foundation to estimate the geophysical parameters of interest.
基金Project supported by the National Natural Science Foundation of China (No, 40371058), and the National Key BasicResearch Support Foundation of China (No. G1999011803)
文摘This research used both geostatistics and GIS approach to compare temporal change of soil salt between 1980 and 2003, to analyze the spatial distribution of surface soil salt, to developed methods for predicting soil salinization potential based on recent improvements to the Dempster-Shafer theory, and to develop probability maps of potential salinization in Yucheng City, China. A semivariogram model of soil salt content was developed from the spherical model, and then employing kriging interpolation the spatial distribution of salt content in 2003 was obtained utilizing data from 100 soil sampling points. Potential salinization distribution was mapped using an approach that integrated soil data of the second general survey in 1980 in Yucheng City, which included groundwater salinity, groundwater depth, soil texture, soil organic matter content, and geomorphic maps. With the support of Dempster-Shafer theory and fuzzy set technique the factors that affected potential soil salinization were characterized and integrated; and then soil salinization was predicted. Finally a prognosis map of potential salinization distribution in the research area was obtained, with higher probability values indicating higher hazards to salinity processes. The distribution of the potential soil salinization probability was a successive surface.
基金funded by the Small Scale Funding Agreement (UNEP/ROWA)
文摘Sand and dust storms （SDS） are common phenomena in arid and semi-arid areas. In recent years, SDS frequencies and intensities have increased significantly in Iran. A research on SDS sources is important for understanding the mechanisms of dust generation and assessing its socio-economic and environmental impacts. In this paper, we developed a new approach to identify SDS source areas in Iran using a combination of nine related datasets, namely drought events, temperature, precipitation, location of sandy soils, SDS frequency, hu- man-induced soil degradation （HISD）, human influence index （HII）, rain use efficiency （RUE） and net primary pro- ductivity （NPP） loss. To identify SDS source areas, we firstly normalized these datasets under uniform criteria in- cluding layer reprojection using Lambert conformal conic projection, data conversion from shapefile to raster, Min-Max Normalization with data range from 0 to 1, and data interpolation by Kriging and images resampling （resolution of 1 km）. After that, a score map for the possibility of SDS sources was generated through overlaying multiple datasets under average weight allocation criterion, in which each item obtained weight equally. In the score map, the higher the score, the more possible a specific area could be regarded as SDS source area. Exceptions mostly came from large cities, like Tehran and Isfahan. As a result, final SDS source areas were mapped out, and Al-Howizeh/Al-Azim marshes and Sistan Basin were identified as main SDS source areas in Iran. The SDS source area in Al-Howizeh/Al-Azim marshes still keeps expanding. In addition, Al-Howizeh/Al-Azim marshes are now suf- fering rapid land degradation due to natural and human-induced factors and might totally vanish in the near future. Sistan Basin also demonstrates the impacts of soil degradation and wind erosion. With appropriate intensity, dura- tion, wind speed and altitude of the dust storms, sand particles uplifting from this area might have developed into extreme dust storms, especially during the summer.
基金Under the auspices of the Specialized Research Fund for the Doctoral Program of Higher Education (No. 20090182120024)National Natural Science Foundation of China (No. 41101568)+1 种基金Natural Science Foundation Project of Chongqing Science & Technology Commission (No. cstcjjA00008)Fundamental Research Funds for the Central Universities (2012XZZX012)
文摘Rapid urbanization results in the conversion of natural soil to urban soil,and consequently,the storage and density of the soil carbon pools change.Taking Chongqing Municipality of China as a study case,this investigation attempts to better understand soil carbon pools in hilly cities.First,the vegetated areas in the study area were derived from QuickBird images.Then,topsoil data from 220 soil samples(0-20 cm) in the vegetated areas were collected and their soil organic carbon(SOC) densities were analyzed.Using the Kriging interpolation method,the spatial pattern of SOC was estimated.The results show that the SOC density exhibited high spatial variability in the urban topsoil of Chongqing.First,the SOC density in topsoil decreased according to slope in the order 2°-6° < 25°-90° < 0°-2° < 6°-15° < 15°-25°.Second,the newly developed areas during 2001-2010 had a lower SOC density than the areas built before 1988.Third,urban parks and gardens had a higher SOC density in topsoil,residential green land followed,and scattered street green land ranked last.For hilly cities,the variability of terrain affects the distribution of SOC.The Kriging results indicate that Kriging method combining slope with SOC density produced a high level of accuracy.The Kriging results show that the SOC density to the north of the Jialing River was higher than the south.The vegetated areas were estimated to amount to 73.5 km2 across the study area with an SOC storage of 0.192 Tg and an average density of 2.61 kg/m2.
基金Under the auspices of Chinese Postdoctoral Science Foundation(No.2015M570864)Open-ended Fund of State Key Laboratory of Cryosphere Sciences,Chinese Academy of Sciences(No.SKLCS-OP-2014-11)+2 种基金Northwest Normal University Young Teachers Scientific Research Ability Promotion Plan(No.NWNU-LKQN-13-10)National Natural Science Foundation of China(No.41273010,41271133)Major National Research Projects of China(No.2013CBA01808)
文摘The Thornthwaite moisture index, an index of the supply of water （precipitation） in an area relative to the climatic demand for water （potential evapotranspiration）, was used to examine the spatial and temporal variation of drought and to verify the influence of environmental factors on the drought in the Hengduan Mountains, China. Results indicate that the Thornthwaite moisture index in the Hengduan Mountains had been increasing since 1960 with a rate of 0.1938/yr. Annual Thomthwaite moisture index in Hengduan Mountains was between -97.47 and 67.43 and the spatial heterogeneity was obvious in different seasons. Thomthwaite moisture index was high in the north and low in the south, and the monsoon rainfall had a significant impact on its spatial distribution. The tendency rate of Thomthwaite moisture index variation varied in different seasons, and the increasing trends in spring were greater than that in summer and autumn. However, the Thomthwaite moisture index decreased in winter. Thomthwaite moisture index increased greatly in the north and there was a small growth in the south of Hengduan Mountains. The increase of precipitation and decrease of evaporation lead to the increase of Thomthwaite moisture index. Thornthwaite moisture index has strong correlation with vegetation coverage. It can be seen that the correlation between Normalized Difference Vegetation Index （NDVI） and Thomthwaite moisture index was positive in spring and summer, but negative in autumn and winter. Correlation between Thornthwaite moisture index and relative soil relative moisture content was positive in spring, summer and autumn, but negative in winter. The typical mountainous terrain affect the distribu- tion of temperature, precipitation, wind speed and other meteorological factors in this region, and then affect the spatial distribution of Thomthwaite moisture index. The unique ridge-gorge terrain caused the continuity of water-heat distribution from the north to south, and the water-heat was stronger than that from the east to west part, and thus determined the spatial distribution of Thornthwaite mois- ture index. The drought in the Hengduan Mountains area is mainly due to the unstable South Asian monsoon rainfall time.
基金We thank the financial support from the National Natural Science Foundation of China(40701007,40571066)the Postdoctoral Science Foundation of China(20060401048).
文摘The acquisition of precise soil data representative of the entire survey area, is a critical issue for many treatments such as irrigation or fertilization in precision agriculture. The aim of this study was to investigate the spatial variability of soil bulk electrical conductivity （ECb） in a coastal saline field and design an optimized spatial sampling scheme of ECb based on a sampling design algorithm, the variance quad-tree （VQT） method. Soil ECb data were collected from the field at 20 m interval in a regular grid scheme. The smooth contour map of the whole field was obtained by ordinary kriging interpolation, VQT algorithm was then used to split the smooth contour map into strata of different number desired, the sampling locations can be selected within each stratum in subsequent sampling. The result indicated that the probability of choosing representative sampling sites was increased significantly by using VQT method with the sampling number being greatly reduced compared to grid sampling design while retaining the same prediction accuracy. The advantage of the VQT method is that this scheme samples sparsely in fields where the spatial variability is relatively uniform and more intensive where the variability is large. Thus the sampling efficiency can be improved, hence facilitate an assessment methodology that can be applied in a rapid, practical and cost-effective manner.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11002054)
文摘A boundary-type meshless method called the scaled boundary node method （SBNM） is developed to directly evaluate mixed mode stress intensity factors （SIFs） without extra post-processing. The SBNM combines the scaled boundary equations with the moving Kriging （MK） interpolation to retain the dimensionality advantage of the former and the meshless attribute of the latter. As a result, the SBNM requires only a set of scattered nodes on the boundary, and the displacement field is approximated by using the MK interpolation technique, which possesses the 5 function property. This makes the developed method efficient and straightforward in imposing the essential boundary conditions, and no special treatment techniques are required. Besides, the SBNM works by weakening the governing differential equations in the circumferential direction and then solving the weakened equations analytically in the radial direction. Therefore, the SBNM permits an accurate representation of the singularities in the radial direction when the scaling center is located at the crack tip. Numerical examples using the SBNM for computing the SIFs are presented. Good agreements with available results in the literature are obtained.
基金Supported by Forest Management Key Subject Construction Project of Southwest Forestry University(XKZ200901)
文摘[Objective] The aim was to study the spatial information of temperature and precipitation data in Hengduan mountains. [Method] Considering GIS spatial interpolation and numerical statistics theory, spatial prediction were carried out to the ten years average temperature and precipitation flux observation data in 109 sparse meteorological stations in Hengduan Mountains. Based on the spatial range of geographic position of Hengduan Mountains, and 1∶1 000 000 scale DEM as data sources, and using trend surface simulation and residual ordinary Kriging interpolation correction method, the spatial continuous surface for annual average temperature and precipitation were studied. [Result] It was scientific and reasonable to use certain unevenly distributed sparse climate observation station value, and by dint of trend simulation and residue interpolation method to get climate consecutive spatial result. This method can not only accurate the temperature and precipitation spatial distributions to grid point, but also can reflect macro and micro geography factors and topographic influence factor of variation. Furthermore, it can be predicted error term trend surface reasonable spatial distribution. Simulation results were basically in accordance with the objective law, and can be used for the region climate data spatial informatization simulation. [Conclusion] The study provided scientific spatial basic data for the further study of ecological and vegetation in Hengduan Mountains.
文摘In abandoned mine sites, i.e., mine sites where mining operations have ended, wide spread contaminations are often evident, but the potential sources and pathways of contamination especially through the subsurface, are difficult to identify due to inadequate and sparse geochemical measurements available. Therefore, it is essential to design and implement a planned monitoring net-work to obtain essential information required for establishing the potential contamination source locations, i.e., waste dumps, tailing dams, pits and possible pathways through the subsurface, and to design a remediation strategy for rehabilitation. This study presents an illustrative application of modeling the flow and transport processes and monitoring network design in a study area hydrogeologically resembling an abandoned mine site in Queensland, Australia. In this preliminary study, the contaminant transport process modeled does not incorporate the reactive geochemistry of the contaminants. The transport process is modeled considering a generic conservative contaminant for the illustrative purpose of showing the potential application of an optimal monitoring design methodology. This study aims to design optimal monitoring network to: 1) minimize the contaminant solute mass estimation error;2) locate the plume boundary;3) select the monitoring locations with (potentially) high concentrations. A linked simulation optimization based methodology is utilized for optimal monitoring network design. The methodology is applied utilizing a recently developed software package CARE-GWMND, developed at James Cook University for optimal monitoring network design. Given the complexity of the groundwater systems and the sparsity of pollutant concentration observation data from the field, this software is capable of simulating the groundwater flow and solute transport with spatial interpolation of data from a sparse set of available data, and it utilizes the optimization algorithm to determine optimum locations for implementing monitoring wells.
基金supported by the National Natural Science Founda-tion of China(11272118)Open Found of State Key Laboratory of Explosion Science and Technology(KFJJ12-5M)
文摘Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes with- out connectivity are divided into several subsets by cells with arbitrary shape. Local discrete equation is established over each cell by using moving Kriging interpolation, in which the nodes that located in the cell are used for approxima- tion. Then local discrete equations can be simplified by con- densation of degree of freedom, which transfers equations of inner nodes to equations of boundary nodes based on cells. The global dynamic system equations are obtained by as- sembling all local discrete equations and are solved by using the standard implicit Newmark＇s time integration scheme. In the scheme of present method, the calculation of each cell is carried out by meshfree method, and local search is imple- mented in interpolation. Numerical examples show that the present method has high computational efficiency and good accuracy in solving elastic dynamic problems.
基金This work was supported in part by the National Nature Science Foundation of China(No.61502240,61502096,61304205,61773219)Natural Science Foundation of Jiangsu Province(BK20150634,BK20141002).
文摘Real-time performance and accuracy are two most challenging requirements in virtual surgery training.These difficulties limit the promotion of advanced models in virtual surgery,including many geometric and physical models.This paper proposes a physical model of virtual soft tissue,which is a twist model based on the Kriging interpolation and membrane analogy.The proposed model can quickly locate spatial position through Kriging interpolation method and accurately compute the force change on the soft tissue through membrane analogy method.The virtual surgery simulation system is built with a PHANTOM OMNI haptic interaction device to simulate the torsion of virtual stomach and arm,and further verifies the real-time performance and simulation accuracy of the proposed model.The experimental results show that the proposed soft tissue model has high speed and accuracy,realistic deformation,and reliable haptic feedback.
