基于等价权抗差估计和线性规划的天线曲面拟合算法

现有天线曲面拟合算法均基于最小二乘法,难以消除多个粗差对拟合结果的影响。本文基于等价权抗差估计思想,利用线性规划计算残差初值,再进行选权迭代计算。通过算例,证明本文方法能较好地探查多个粗差,且计算结果精度更高。     

Bounding surface model for soil resistance to cyclic lateral pile displacements with arbitrary direction

The development of a two-surface elastic–plastic bounding surface P–Y model for cyclic lateral pile motions is described. The kinematic-hardening model is applicable to the analysis of pile foundations subjected to loading with arbitrary azimuths relative to the pile axis. The model realistically captures the hysteretic energy damping associated with dynamic loading of subsea foundations through physically correct plastic mechanisms and provides results consistent with those observed in physical tests including cyclic loading. Its performance is demonstrated in element states of stress and in pile foundation analyses. The development based on the incremental theory of plasticity results in more robust solutions than may be obtained using alternative elastic, variable moduli and deformation plasticity formulations.     

线图形简化与移位算子的协同方法

在地图综合过程中,线的图形简化和移位算子通常是分别执行的,图形简化和移位有时都会产生新的空间冲突。本文试图把这两种算子进行协同,避免在地图综合过程中进行空间冲突的多次探测,提高地图综合数据处理的效率。本文通过把线图形简化转换为线上的点移位,并构建邻近地图目标之间的移位传播路径,在考虑空间上下文关系和地图感受规则的前提下,使移位过程能考虑到线图形简化,并尽量保持有关地图目标的空间特征。最后,以道路和其周边的建筑物群为例,验证了该算法的有效性和可行性。     

应用两次展开法求解系泊柱体慢漂运动方程

在深海中系泊的海洋平台,如Spar平台,水下部分为带有系泊的圆柱结构,其水平方向运动响应往往具有较低的自振频率,容易在低频波浪力(源于非线性的差频效应)作用下发生共振响应,使结构发生大幅水平慢漂。当浮体的瞬时位置大幅偏离初始位置时,基于初始平衡位置的摄动展开法会存在较大误差。针对这一问题,采用两次展开方法,对大幅慢漂运动开展时域模拟研究。对双色波作用下自由漂浮圆柱的大幅运动响应问题进行数值分析,并与采用基于初始平衡位置的摄动展开法的计算结果进行了对比。结果表明,采用新的两次展开法可以计算出波浪遭遇频率的变化和波浪漂移阻尼,而这无法从基于初始平衡位置的常规摄动展开法中得到,体现了两次展开法在分析大幅慢漂问题上的优势。     

Application of Linear Mean-Square Estimation in Ocean Engineering

The attempt to obtain long-term observed data around some sea areas we concern is usually very hard or even impossible in practical offshore and ocean engineering situations. In this paper, by means of linear mean-square estimation method, a new way to extend short-term data to long-term ones is developed. The long-term data about concerning sea areas can be constructed via a series of long-term data obtained from neighbor oceanographic stations, through relevance analysis of different data series. It is effective to cover the insufficiency of time series prediction method’s overdependence upon the length of data series, as well as the limitation of variable numbers adopted in multiple linear regression model. The storm surge data collected from three oceanographic stations located in Shandong Peninsula are taken as examples to analyze the number-selection effect of reference oceanographic stations (adjacent to the concerning sea area) and the correlation coefficients between sea sites which are selected for reference and for engineering projects construction respectively. By comparing the N-year return-period values which are calculated from observed raw data and processed data which are extended from finite data series by means of the linear mean-square estimation method, one can draw a conclusion that this method can give considerably good estimation in practical ocean engineering, in spite of different extreme value distributions about raw and processed data.     

Application of Geochemistry and VNIR Spectroscopy in Mapping Heavy Metal Pollution of Stream Sediments in the Takab Mining Area, NW of Iran

This study considered the possibility of using visible and near infrared (VNIR) spectral absorption feature parameters (SAFPs) in predicting the concentration and mapping the distribution of heavy metals in sediments of the Takab area. In total, 60 sediment samples were collected along main streams draining from the mining districts and tailing sites, in order to measure the concentration of As, Co, V, Cu, Cr, Ni, Hg, Ti, Pb and Zn and the reflectance spectra (350–2500 nm). The quantitative relationship between SAFPs (Depth500nm, R610/500nm, R1344/778nm, Area500nm, Depth2200nm, Area2200nm, Asym2200nm) and geochemical data were assessed using stepwise multiple linear regression (SMLR) and enter multiple linear regression (EMLR) methods. The results showed a strong negative correlation between Ni and Cr with Area2200nm, a significant positive correlation between As and Asym2200nm, Ni and Co with Depth2200nm, as well as Co, V and total values with Depth500nm. The EMLR method eventuated in a significant prediction result for Ni, Cr, Co and As concentrations based on spectral parameters, whereas the prediction for Zn, V and total value was relatively weak. The spatial distribution pattern of geochemical data showed that mining activities, along with the natural weathering of base metal occurrences and rock units, has caused high concentrations of heavy metals in sediments of the Sarough River tributaries.     

矿石体重回归分析在新疆乌宗布拉克矿区的应用研究

在资源储量估算的过程中,矿石小体重的准确与否将直接影响到资源储量估算的客观程度。对矿石小体重的准确预测将需要通过大量的统计分析工作,传统的方法是利用矿石小体重的算术平均值进行估算矿床资源储量,并没有考虑到多矿种对矿石小体重的影响程度。本文利用Excel软件"数据分析"的"回归"功能模块,对实验室所测定的矿石小体重值与其对应品位进行二元线性回归分析,快速、准确地构建矿石体重与其品位的数学模型,从而为资源储量估算提供了更客观、更科学的矿石体重模型。     

确定砂土动剪切模量和阻尼比的方法对比

确定土体动剪切模量的常用方法有规范法、Kumar法和自相关函数法,确定相应阻尼比的方法有规范法、Das and Luo法、Kokusho法、Kumar法和互相关函数法,为了分析不同方法所产生差异,实现定量化对比分析,笔者以福建标准砂（粒径为0.5~1.0 mm）为研究对象,采用不排水的应力控制动三轴试验,探讨不同的确定土体动剪切模量和阻尼比方法的差异性,并给出了不同土体条件建议选用的方法。结果表明：1）3种方法确定动剪切模量的结果有一定的差异,随剪应变的增大结果的差异逐渐增大,有效围压对结果的差异性有所影响,当剪应变为4×10^-3,有效围压为100 kPa时,3种方法差异显著,相对误差最大接近20%;2）而5种方法确定阻尼比的结果差异显著,随着剪应变的增大,5种方法确定的阻尼比相对误差大体上均在迅速减小,只有规范法在有效围压为100 kPa时,其相对误差有较小的增大趋势;5种方法中,Kumar法确定的阻尼比最接近平均阻尼比,互相关函数法远高于平均阻尼比,Das and Luo法和Kokusho法确定的阻尼比基本一致但低于平均阻尼比。建议以后的工程应用中,加载方式为应力控制时,可采用自相关函数法确定动剪切模量,采用Kumar法确定阻尼比,二者确定的动剪切模量和阻尼比均最接近平均值。     

线性化物理过程对GRAPES 4DVAR同化的影响

线性化物理过程能够改善四维变分同化中极小化收敛的稳定性和增加极小化过程中对大气物理过程和动力更加精确的描述,它是四维变分同化中非常重要的一部分。通过在GRAPES全球模式中研究线性化物理过程,尤其是两个湿线性化物理过程,改善切线性模式预报精度,来提高GRAPES全球四维变分同化的分析和预报效果。线性化物理过程的开发首先需要简化原非线性化物理过程中的强非线性项,然后对线性化物理过程进行规约化,以抑制切线性扰动的异常增长。目前GRAEPS全球模式中的线性化物理过程主要包括次网格尺度地形参数化、垂直扩散、积云深对流和大尺度凝结。线性化物理过程预报精度的检验方法是通过选择合适大小的初始扰动（同化分析增量）,来比较非线性模式和切线性模式中的扰动演化的纬向平均误差。然后以绝热版本的切线性模式为基础,通过冬、夏两个个例试验来分别检验4个线性化物理过程的12 h预报效果。试验结果表明,通过添加次网格地形参数化和垂直扩散两个干线性化物理过程方案,可以有效抑制住绝热版本切线性模式低层扰动的异常增长,大幅度改善切线性模式预报效果。通过添加积云深对流和大尺度凝结两个湿线性化物理过程,可以在热带区域和中、高纬度地区提高切线性模式中湿变量和温度变量的近似精度,提高GRAPES全球四维变分同化的分析和预报效果。      

GRAPES全球奇异向量方法改进及试验分析

基于总能量模的奇异向量扰动常用于构造集合预报的初始条件。以建立GRAPES（Global and Regional Assimilation PrEdiction System）全球集合预报系统为目的,基于前期研发的GRAPES全球模式奇异向量方法,在GRAPES全球切线性模式和伴随模式2.0版的框架下,开展了引入线性化边界层方案来改善奇异向量结构,并提高奇异向量计算效率的研究。通过连续试验,从奇异向量的扰动能量结构、扰动能量谱及扰动空间分布等方面,综合分析改进GRAPES全球奇异向量的结构及演变特征。试验结果表明,改进后的GRAPES奇异向量方法有效抑制了之前扰动能量在近地面层不合理的快速增长,同时,奇异向量最优扰动的结构更客观地体现了中高纬度区域大气初始条件中的斜压不稳定扰动及其演变,如在初始时刻奇异向量扰动能量主要位于对流层中层,并呈现出随高度向西倾斜的大气斜压特征;经过线性化演变,扰动能量向较大水平尺度转移,并在垂直结构上表现出向对流层高层上传及向对流层低层下传的特征等。针对GRAPES奇异向量迭代求解中伴随模式计算耗时为主的情况,改进伴随模式中广义共轭余差方案的调用方式,并采用大内存存储法来提高其计算效率,进而将奇异向量总计算时间缩短了25%。总之,改进后的GRAPES奇异向量方法,可应用于构建面向业务应用的GRAPES全球集合预报系统。