共查询到20条相似文献，搜索用时 140 毫秒
1.
2.
不同地形与取样数的Kriging插值精度对比研究——以舒城县土壤全氮空间分异为例 总被引：4，自引：0，他引：4
为分析地形与取样数对大尺度Kriging插值精度的关系,在舒城县的山区、丘陵、岗区、平原等4个不同地形单元采集土壤表层(0~20cm)样品611个,进行可变取样数条件下的土壤全氮空间分异插值精度比较。半变异函数分析表明,4个地形单元块金值与基台值比值的平均值<25%,具有强烈的空间相关性。Kriging插值精度标准均方根误差检验值(NRMSE)在各地形单元内都随着采样密度的增大而不断降低,插值精度提高;各地形单元之间,在相同的采样密度下,Kriging插值精度大小顺序为平原>岗区>丘陵>山区,随着地形复杂程度增加,Kriging插值精度降低。NRMSE检验值平均变化率分析表明,地形越复杂地区Kriging插值精度提高越快。地形条件和取样数是影响大尺度Kriging插值精度的重要因素。 相似文献
3.
对h型自适应自然单元法的自适应细化方案进行了初步研究。在ZZ误差分析的基础上实现了节点的自动加密,使得随着自适应细化的进行,求解误差减小,而且误差分布趋于均匀。在数值分析中,主要有两种误差来源--插值误差和积分误差。随着节点的加密,Delaunay三角形的尺寸随之减小,三角形内的应力场趋于线性分布,那么插值误差和积分误差也都会随之减小。因而,h型自适应分析可以同时减小上述两种误差而达到不断提高求解精度的目的。由于自然单元法求解依赖于求解域内离散节点的Voronoi结构,建议的细化方案中新节点的引入只需局部调整Delaunay结构,算法的实现极为容易,程序实现简单、高效。研究表明,建议的自适应方案是可行的,自然单元法特别适合进行h型自适应分析。 相似文献
4.
大地电磁(MT)数值模拟中通常使用有限单元法,通过伽辽金(Galerkin)法将微分方程转化为与其等价的泛函形式,对泛函求取极值并在单元上定义插值基函数,得到节点上电磁场值的线性方程组,最终形成大型复对称稀疏矩阵。要达到较高的有限元计算精度,一般采用密集的网格或高次插值的方法,这样做大大的减慢了正演的速度。结合两者的优点利用三次插值和h-型自适应相结合的有限元法来实现MT的正演算法。首先从一个粗网格出发并利用三次插值,通过后验误差估计方法局部加密网格,在计算量较小的情况获得较高的计算精度。这种方法可以针对目标区域和介质分界面发生突变处进行网格加密,不需要全局加密网格。最后通过对国际标准模型COMMEMI-2D1的模拟,分别比较二次插值与三次插值的自适应网格数量和数值模拟结果,证明了三次插值自适应有限元算法的可行性。 相似文献
5.
2.5D有限元方法在铁路路基动力响应研究领域中的应用渐趋广泛。针对其在求解随机不平顺条件下路基动力响应时计算效率显著下降的问题,构建了基于二维降阶Hermite插值的2.5D有限元路基动力响应快速计算框架。以路基在频率-波数域动力响应的基本特征为依据确定了插值原则,讨论了插值点分布和数量对插值精度的影响。研究表明：采用二维降阶Hermite插值方法可以实现随机不平顺条件下路基动力响应的快速计算。相比插值点非均匀分布,插值点均匀分布可以兼顾幅值和相位的插值精度,适应性更好。此外,该方法的计算效率仅与插值点数量相关,不受随机不平顺谐波数量的影响,在模拟随机不平顺条件下路基动力响应方面具备显著的优势。 相似文献
6.
为了研究采样和网格化方法对地球物理数据成图精度的影响,为野外数据采集布设提供一定的依据,采用数值模拟确定重力异常场场值,通过不同采样间距和不同插值方法计算重力异常绝对误差均方根值和节点处的绝对误差值,对比不同插值方法的误差,得到了如下认识：1）对于同一插值方法而言,存在小间距绝对误差均方根值小于大间距绝对误差均方根值的关系。2）对不同的插值方法而言:当采样间距小于最小异常地质体尺度时,绝对误差均方根值由小到大的顺序是径向基函数法、改进的谢别德法、克里金插值法、自然邻点法、反距离加权插值法、最近邻点法、最小曲率法,并且线性插值三角网法与自然邻点法具有几乎相同的数值;当采样间距大于最小异常地质体尺度时,绝对误差均方根值由小到大的顺序是径向基函数法、改进的谢别德法、克里金插值法、自然邻点法、最小曲率法、最近邻点法、反距离加权插值法,并且线性插值三角网法和自然邻点法具有几乎相同的数值。3）从绝对误差均方值看,径向基函数方法、改进的谢别德方法和克里金方法数值较小,其中径向基函数值绝对误差均方根值最小。4）从节点处绝对误差值来看,径向基函数方法、克里金方法、改进的谢别德方法相对其他插值方法具有更小的误差,不存在局部误差较小或较大的情况,是相对较好的插值方法,并且径向基函数方法是最好的。 相似文献
7.
位场延拓是重、磁位场数据处理的重要方法之一,高精度的位场延拓结果对后续的数据处理和解释尤为重要。笔者从平面位场延拓的基本公式出发,分析了空间域和频率域位场延拓结果精度的影响因素以及稳定性。通过理论模型测试比较了这些影响因素(场源体顶面埋深、剖面长度、扩边方法、窗口大小、点距和延拓高度)在空间域和频率域进行位场延拓时的异同性。经过测试表明,点距和延拓高度对延拓结果的影响最大,其次是剖面长度、扩边方法以及窗口大小,场源体的顶面埋深影响最小。随着顶面埋深的增大,会使延拓结果的精度降低;增大剖面长度和进行扩边会提高位场延拓结果的精度;选择合适滑动窗口可以提高计算效率,窗口越大延拓结果精度越高,窗口半径一般选择20倍延拓高度。空间域和频率域中位场延拓结果精度相当,但在空间域中,当延拓高度小于1倍点距时,延拓结果误差很大,此时需要用插值方法加密点距。 相似文献
8.
频率域激发极化法有限元数值模拟 总被引：2，自引：1，他引：1
在频率较低和忽略电磁效应的情况下,利用有限单元方法和Cole-Cole模型对频率域激发极化法进行数值模拟.首先在三维地电条件下,给出电场的边值问题和变分问题,将Cole-Cole模型的频率响应引入到地电模型中,运用有限单元法对模型进行单元剖分、插值、积分和总体合成等,通过解方程最后得到表征频率域极化强度的参数幅频率.通过改变地电模型的参数,得到了不同的幅频率响应曲线.模拟结果与实际情况符合,表明该方法是正确和适用的. 相似文献
9.
10.
岩土工程中经常会遇到无穷域问题,而采用无限单元可以实现对其有效地模拟。弱形式求积元法是一个有效的数值工具,它常通过提高积分阶次来提高计算精度。建立了无限弱形式求积单元并被应用于求解岩土工程中的无穷域问题,该单元基于坐标映射,将无穷域变换到标准域,在标准域上进行数值积分和数值微分,保留了传统弱形式求积元的积分点坐标和权系数。求解了瞬态渗流、固结和静力分析等数值算例,并与解析解或截断方法进行了对比。结果表明:基于坐标映射的无限弱形式求积单元使用简单,可以模拟各种类型的无穷域问题,仅需要将感兴趣的范围进行有限域划分并通过提高积分阶次来减小对极点位置的依赖,极大地节省了计算资源,提高了计算精度。 相似文献
11.
12.
Interpolation for geochemical surface reconstruction incorporating topographic catchment definitions
Geochemical surfaces are reconstructed by interpolating geochemical measurements obtained from stream-water and stream-sediment samples. The geographical region that influences (and therefore is represented by) the value of a geochemial sample is its topographic catchment area. However, standard convention is to treat and to record the stream sample in the database as a point location, and to reconstruct geochemical surfaces utilizing conventional point interpolation procedures. These interpolation procedures assume, generally, that a data point exerts geographical influence away from itself in all directions, and that influence declines with distance away from that data point. Conventional interpolation procedures are poorly suited for reconstructing geochemical surfaces from stream samples; they do not take into account the true geographic area that geochemical sample points represent (topographic catchments). In this paper we propose a method of interpolation which assumes that data points are representative of their topographic catchment areas. Experimental data indicates that a surface reconstruction procedure which preserves the areal character of geochemical samples provides results more meaningful than surfaces reconstructed using more conventional interpolation techniques. 相似文献
13.
提出了一种模拟裂纹扩展的水平集和无网格耦合方法。由于水平集和无网格方法都是建立在离散节点上,因而可以很自然地实现耦合。在该方法中,两个在裂尖处相互正交的水平集不仅用于描述裂纹的几何形态和裂尖位置,而且用于建立无网格伽辽金法（简称EFGM）不连续近似函数中的Heaviside跳跃项和裂尖处的Westergaard扩展项。当裂纹扩展时,则由水平集更新算法确定新裂纹的位置。水平集和无网格耦合法无需使用可视法、衍射法或透明法,克服了这些方法在裂尖处人为引入的不连续且能很好地再生 奇异场;而且节点影响域不受裂纹线切割的影响,在计算中往往使用较小的影响域,保持了整体刚度矩阵的带状、稀疏性;另外,水平集简化了扩展节点的选取和附加函数的建立,其更新过程无需求解演化方程,实现简单且易于编程。数值算例表明本文方法具有较高的计算精度,其模拟的裂纹扩展路径与试验结果吻合得很好,从而验证了本文方法的正确性和可行性。 相似文献
14.
Several kinds of data can provide information about a variable measured on a one- or two-dimensional space; at some points, the value is known to be equal to a certain number. At other points, the only information may be that the variable is greater or smaller than a given value. The theory of splines provides interpolating functions that can take into account both equality and inequality data. These interpolating functions are presented. The parallel between splines and kriging is reviewed, using the formalism of dual kriging. Coefficients of dual kriging can be obtained directly by minimizing a quadratic form. By adding some inequality constraints to this minimization, an interpolating function may be calculated which takes into account inequality data and is more general than a spline. The method is illustrated by some simple one-dimensional examples.Work performed at Sohio Petroleum Company 相似文献
15.
Shujun Peng Zhennan Zhang Jianye Mou Bing Zhao Zhiyuan Liu Jiading Wang 《国际地质力学数值与分析法杂志》2019,43(1):441-460
The improved element partition method (IEPM) is a newly developed fracture simulation approach. IEPM allows a fracture to run across an element without introducing extra degrees of freedom. It can also simulate any number of fractures in a prescribed mesh without remeshing. In this study, the IEPM is extended to hydraulic fracture simulation. First, the seepage and volumetric storage matrix of a cracked element are derived using virtual nodes (the intersection points of a crack with element edges). Subsequently, the fully coupled hydromechanical equation is derived for this cracked element. To eliminate the extra degrees of freedom (virtual nodal quantities), the water pressure and displacement of the virtual nodes are associated with their adjacent nodes through least squares interpolation. Finally, the fully coupled equation in terms of nodal quantities is obtained. The verification cases validate the method. By using this method, the field-scale hydraulic fracturing process is well simulated. The proposed approach is simple and efficient for field-scale hydraulic fracture simulation. 相似文献
16.
Some commonly used interpolation algorithms are analyzed briefly in this paper. Among all of the methods, biharmonic spline
interpolation, which is based on Green’s function and proposed by Sandwell, has become the mainstream method for its high
precision, simplicity and flexibility. However, the minimum curvature method has two flaws. First, it suffers from undesirable
oscillations between data points, which is solved by interpolation with splines in tension. Second, the computation time is
approximately proportional to the cube of the number of data constraints, making the method slow for situations with dense
data coverage. Focusing on the second problem, this paper introduces the moving surface spline interpolation method based
on Green’s function, and the interpolation error equations are deduced. Because the proposed method only chooses the nearest
data points by using the merge sort algorithm for interpolating, the computation time is greatly decreased. The optimal number
of the nearest points can be determined by using the interpolation error estimation equation. No matter how many data points
there are, this method can be implemented without difficulty. Examples show that the proposed method can obtain high interpolation
precision and high computation speed at the same time. 相似文献
17.
Because of the need for computational efficiency, bivariate interpolation methods applied to scattered observations often involve two stages. Initially the variable is estimated at regular grid nodes using a running subset of data (usually of fixed number). This, however, will produce discontinuities in the interpolated surface. Thus a second stage, curvilinear interpolation technique, is applied to estimated values to smooth out the effect of discontinuities. Such problems can be overcome efficiently in processing large data sets by interpolating over natural neighbor subsets. Interpolation procedures that generate discontinuities in the interpolated surface are inappropriate for geological applications, where dislocations due to structural complications may be present. 相似文献
18.
The development of stability problems related to classical mixed methods has recently been observed. In this study, a soil-water coupled boundary-value problem, one type of stability problem, is presented using the element-free Galerkin method (EFG method). In this soil-water coupled problem, anomalous behavior appears in the pressure field unless a stabilization technique is used. The remedy to such numerical instability has generally been to adopt a higher interpolation order for the displacements than for the pore pressure. As an alternative, however, an added stabilization term is incorporated into the equilibrium equation. The advantages of this stabilization procedure are as follows: (1) The interpolation order for the pore pressure is the same as that for the displacements. Therefore, the interpolation functions in the pore pressure field do not reduce the accuracy of the numerical results. (2) The stabilization term consists of first derivatives. The first derivatives of the interpolation functions for the EFG Method are smooth, and therefore, the solutions for pore pressure are accurate. In order to validate the above stabilization technique, some numerical results are given. It can be seen from the results that a good convergence is obtained with this stabilization term. 相似文献
19.
针对线弹性断裂力学问题,提出扩展径向点插值无网格法（X-RPIM）。该方法基于单位分解思想,在传统径向点插值无网格法的位移模式中加入扩展项来描述裂纹两侧的不连续位移场和裂尖奇异场。由于其形函数具有Kronecker ? 函数性质,易于施加本质边界条件。详细描述了X-RPIM不连续位移模式的建立,支配方程的离散形式以及J积分计算混合模式裂纹的应力强度因子的实现过程,讨论了不同积分区域对应力强度因子的影响。数值算例分析证明了该方法在求解断裂问题时的可行性和有效性,同时说明扩展径向点插值无网格法在模拟裂纹扩展问题时具有良好的前景。 相似文献
20.
岩土体的渗透破坏、地下工程的防渗设计等无不与渗流计算有关。针对渗流自由面问题,提出一种重心拉格朗日插值的配点型无网格方法。由于渗流自由面问题的求解区域是不规则区域,该方法通过将不规则求解区域嵌入一个正则矩形区域,在正则区域上采用重心拉格朗日插值近似未知函数,利用配点法离散渗流问题的控制方程,将重心拉格朗日插值的微分矩阵离散成代数方程表达的矩阵形式。将自由面上的边界条件通过重心拉格朗日插值离散,通过置换方程法和附加方程法施加边界条件,利用正则区域上的重心插值配点法,通过迭代确定最终自由面的位置。数值算例表明所提出的无网格方法对于求解渗流自由面问题的正确性和高精度。 相似文献