土壤稳定入渗速率空间插值不确定性研究

以黄土丘陵沟壑第三副区的藉河流域为研究区,根据65个实测点数据,采用普通克里格法、反距离权重法、样条函数法等插值方法,分析了测点数量变化、栅格像元尺寸变化及插值方法的差异对土壤稳定入渗速率空间插值结果的影响,剖析了空间插值中的不确定性。结果表明:(1)参与插值站点越多,所得插值结果不确定性越小;(2)像元尺寸在25~800 m间变化对土壤稳定入渗速率的插值结果影响微弱;(3)不同插值方法对插值结果的精度影响较大,说明插值方法的差异对插值结果的不确定性有较大影响。

陕西省水蚀土壤因子指标插值方法比较研究

基于地统计学原理和半变异理论,采用地理信息系统技术定量研究了陕西省水蚀土壤因子指标的空间变异特征。分析比较了反距离权重法、样条函数法与普通克吕格法对陕西省土壤因子指标空间插值的精度。结果表明,研究区土壤抗冲系数具有中等强度的空间相关性,块金系数为32.29%,而稳渗速率、崩解速率、抗剪强度均表现为强烈的空间相关性,块金系数分别为13.19%,11.61%和12.98%。综合考虑平均相对误差、均方差及插值效果,认为普通克吕格法最好,更能反映土壤参数的空间特征并符合区域水土流失模型对数据的要求。对于普通克吕格法,稳渗速率的Lag步长为30 000 m,半方差理论模型为指数模型;抗冲系数、崩解速率、抗剪强度的Lag步长为55 000,半方差理论模型均为高斯模型。在空间分布上,各指标随土壤类型由北到南呈现明显的地带性规律。

基于等高线的黄土高原小流域地形虚拟

黄土高原地形三维虚拟是"数字黄土高原"的基础,可为区域水土保持生态建设提供科技支撑。针对直接在地理信息系统软件中观察三维场景存在的控制交互能力不足问题,提出综合利用地理信息系统软件的地形插值算法,基于MFC框架下的OpenGL程序设计的思路,实现地形的真实感三维虚拟。以黄土丘陵沟壑区康家沟小流域为例,等高线数据在AutoCAD和ArcView软件中处理,生成ASCII格式的规则网格DEM数据,依据它们绘制三角形带,采用加权平均法求得各点的法向量,设置光照与材质模式,添加动态天空背景,实现了该流域地形的真实感三维虚拟,并增加交互能力,完成自由漫游与多角度观察。

喷灌均匀系数的三次样条两次插值计算方法

为提高喷灌水量分布均匀性评价的准确性,当雨量筒径向布置时,为考虑所有测点数据对插值点降水深的影响,采用径向和周向两次的三次样条插值计算出未知点的降水深,从而计算喷灌均匀系数。以美国雨鸟30PSH型喷头雨量筒间隔为1m和2m的喷洒试验数据,计算网格点取1m和0.25m,分别采用三次样条两次插值法和邻近四点距离线性插值法计算了克里斯琴森均匀系数。结果表明,均匀系数由高至低的顺序依次为采样间隔为2m的线性插值、采样间隔为2m的三次样条两次插值、采样间隔为1m的线性插值和采样间隔为1m的三次样条两次插值。采样间隔2m比1m计算出的均匀系数总体高3~4个百分点,三次样条两次插值法比邻近点距离线性插值法略低1个百分点,2种计算网格点间距下的均匀系数差值小于1个百分点。结果证明,采样间距、插值方法、计算网格间距对均匀系数的影响依次降低,三次样条两次插值法可以用来评价喷灌组合均匀系数。

地形可视性分析中观察点设置问题研究

观察点设置问题是地形可视性分析中的一类重要问题，对该问题的研究可以在空间信息辅助决策、通信、旅游、野生动物保护等领域发挥重大作用。本文在对地形可视性分析中观察点设置问题现有研究成果总结和分析基础上对该问题展开深入研究。 首先，针对现有解决方法只从智能算法或地形数据表示方法单一角度进行分析和研究的局限性，提出了一种问题相关的智能算法和数据表示方法相结合的解决问题新框架。该框架考虑了解决观察点设置问题时智能算法的优点和数据表示方式的特点相互配合问题，目的是充分发挥二者各自的优势以提高观察点设置问题解决的准确度与效率。 其次，在深入分析观察点设置问题本身特点的基础上，结合隶属云理论的基本理论和方法，对经典模拟退火算法从退温函数设计、温度产生过程、状态生成过程三方面进行了问题相关的改进，提出了一种适于观察点设置问题的改进模拟退火算法(Improved Simulated Annealing algorithm, ISA)。该算法一方面保持了经典模拟退火算法的稳定倾向特性，保证了算法满足伴随退火温度的不断下降，对恶化的新状态越来越难于接受这一模拟退火算法的最基本特征；另一方面其退火温度的连续性随机变化特性和隐含的“回火升温”过程，则有利于算法有效拒绝恶化解，加速算法收敛，能够更好地满足观察点设置问题对于算法收敛速度的要求。 再次，在分析地形数据的精度、误差等因素对观察点设置问题的解决准确性和解决效率影响程度的基础上，提出了一种基于离散余弦变换的地形数据内插方法(Discrete Cosine Transformation Interpolation method, DCTI)。新方法将传统空域上的地形内插转换到变换域上进行，同时充分利用了离散余弦变换的熵保持特性和能量压缩特性，简化了变换域上的内插过程，提高了地形数据内插的效率和精度。DCTI方法与其他现有典型地形数据内插方法相比，对地形可视性信息获取的准确性和效率影响最小，为平衡观察点设置问题解决过程中时间效率和准确度之间的关系，最终有效地解决观察点设置问题提供了数据基础。 最后，从智能算法和地形数据相结合的角度出发，提出了一种基于ISA和DCTI相结合的观察点设置问题多分辨率处理方法(Multi-Resolution Processing method, MRP)。新方法将模拟退火算法的逐次退火特点和地形数据的多分辨率表示充分结合，达到了发挥算法数据相结合的综合优势的目的。与现有单纯基于模拟退火算法的解决方法相比，在问题解决准确度保持不变的前提下，基于MRP方法的观察点设置问题解决的平均耗时减少85%~95%，为实际工程应用问题的解决提供了一条重要途径。

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Application of the characteristic CIP method to a shallow water model on the sphere

Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the efficiency of the characteristic constrained interpolation profile (CIP) method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.

A multi-moment finite volume formulation for shallow water equations on unstructured mesh

A novel and accurate finite volume method has been presented to solve the shallow water equations on unstructured grid in plane geometry. In addition to the volume integrated average (VIA moment) for each mesh cell, the point values (PV moment) defined on cell boundary are also treated as the model variables. The volume integrated average is updated via a finite volume formulation, and thus is numerically conserved, while the point value is computed by a point-wise Riemann solver. The cell-wise local interpolation reconstruction is built based on both the VIA and the PV moments, which results in a scheme of almost third order accuracy. Efforts have also been made to formulate the source term of the bottom topography in a way to balance the numerical flux function to satisfy the so-called C-property. The proposed numerical model is validated by numerical tests in comparison with other methods reported in the literature. (C) 2010 Elsevier Inc. All rights reserved.

不同积分变分原理的统一

依据定量因果原理的数学表示,统一地导出了Lagrange量中含坐标关于时间一阶、二阶导数的积分型的Hamilton原理、Voss原理、Hlder原理和Maupertuis-Lagrange原理等,给出了这些原理的本质联系和统一描述.得出f0=0并不是通常的保持Euler-Lagrange方程不变的结果,而是满足定量因果原理的结果.还得出Lagrange量的所有的积分型变分原理等价地对应于两类满足定量因果原理的不变形式.同时发现所有积分型变分原理的运动方程都是Euler-Lagrange方程,但不同条件的变分原理所对应的不同群G作用下的守恒量是不同的.从而可对过去众多零散的积分型变分原理有一个系统和深入的理解,并使这些变分原理自然地成为定量因果原理的推论.

BRST invariant theory of a generalized 1+1 Dimensional nonlinear sigma model with topological term

We give a generalized Lagrangian density of 1 + 1 Dimensional O( 3) nonlinear sigma model with subsidiary constraints, different Lagrange multiplier fields and topological term, find a lost intrinsic constraint condition, convert the subsidiary constraints into inner constraints in the nonlinear sigma model, give the example of not introducing the lost constraint. N = 0, by comparing the example with the case of introducing the lost constraint, we obtain that when not introducing the lost constraint, one has to obtain a lot of various non-intrinsic constraints. We further deduce the gauge generator, give general BRST transformation of the model under the general conditions. It is discovered that there exists a gauge parameter beta originating from the freedom degree of BRST transformation in a general O( 3) nonlinear sigma model, and we gain the general commutation relations of ghost field.

Unified expressions of all integral variational principles

In terms of the quantitative causal principle, this paper obtains a general variational principle, gives unified expressions of the general, Hamilton, Voss, Holder, Maupertuis-Lagrange variational principles of integral style, the invariant quantities of the general, Voss, Holder, Maupertuis-Lagrange variational principles are given, finally the Noether conservation charges of the general, Voss, Holder, Maupertuis-Lagrange variational principles axe deduced, and the intrinsic relations among the invariant quantities and the Noether conservation charges of all the integral variational principles axe achieved.

分离扇回旋加速器的接受度及效率；Acceptance and efficiency of sector-separated cyclotron

以9.7MeV/u的238U36+,5.62MeV/u的70Zn10+为典型离子,分析并模拟了分离扇回旋加速器(SSC)的注入、加速和引出,得到了SSC在理论等时场下横向和纵向的接受度。为了研究SSC在实际情况下的接受度,在实测场的基础上采用Kr-Kb方法以及Lagrange插值方法建立了与实际比较符合的等时场,计算了该等时场下SSC横向和纵向的接受度,发现了导致SSC实际接受度和传输效率低的主要原因在于注入系统中的MSI3元件和引出系统中的MSE3元件设计存在缺陷。模拟结果显示,通过改变MSI3和MSE3的曲率或者垫铁改变元件内部的场分布可以改善SSC的实际接受度和传输效率。

Determination of Point of Zero Charge of Tunisian Kaolinites by Potentiometric and Mass Titration Methods

This paper deals with determining points of zero charge of natural and Na+-saturated mineral kaolinites using two methods: (1) acid-base potentiometric titration was employed to obtain the adsorption of H+ and OH- on amphoteric surfaces in solutions of varying ionic strengths in order to determinate graphically the point of zero net proton charge (PZNPC) defined equally as point of zero salt effect (PZSE); (2) mass titration curve at different electrolyte concentrations in order to estimate PZNPCs by interpolation and to compare with those determined by potentiometric titrations. The two methods involved points of zero charge approximately similar for the two kaolinites between 6.5-7.8, comparable to those reported previously and were in the range expected for these clay minerals. The comparison of potentiometric surface titration curves obtained at 25 °C and those published in the literature reveals significant discrepancies both in the shape and in the pH of PZNPCs values.

区域水土流失研究中降雨量空间插值方法选取

在区域水土流失模型研究中,空间插值可提供每个计算栅格的气象要素资料。考虑到研究区域降雨与高程相关性很弱,不宜采用梯度距离反比法(GIDS),故采用距离反比法(IDW)和普通克里格法(Kriging),对延安示范区及其周围共50个站点2000—2003年的5—10月逐月降雨量进行插值。交叉验证结果表明:对2种插值方法,二者经过对数变换后平均相对误差(MRE)为8.30%和7.67%,分别比原始数据插值后的MRE下降了23.17%和23.50%,说明插值精度得到了提升,对研究区域某一年逐月降水的插值Kriging方法比IDW方法更加精确。

基于直方图插值的均值移动小尺寸目标跟踪算法

小尺寸目标跟踪是视觉跟踪中的难题。该文首先指出了均值移动小尺寸目标跟踪算法中的两个主要问题:算法跟踪中断和丢失跟踪目标。然后,论文给出了相应的解决方法。对传统Parzen窗密度估计法加以改进,并用于对候选目标区域的直方图进行插值处理,较好地解决了算法跟踪中断问题。论文采用Kullback-Leibler距离作为目标模型和候选目标之间的新型相似性度量函数,并推导了其相应的权值和新位置计算公式,提高了算法的跟踪精度。多段视频序列的跟踪实验表明,该文提出的算法可以有效地跟踪小尺寸目标,能够成功跟踪只有6×12个像素的小目标,跟踪精度也有一定提高。