非线性参数估计的自适应松弛正则化算法

非线性方程参数估计存在的弊端在于非线性观测方程存在不适定问题时，以线性化平差估计和高斯牛顿为代表的经典数值算法会产生较强的不稳定特征。因此，针对传统非线性最小二乘求解不稳定且可靠性低的特点，基于稳定泛函极小准则最优化思想，提出了一种自适应松弛正则化数值算法。该算法采用正则化参数几何递增计算方法和残差最小步长准则，实现了正则参数和迭代步长计算的完全自适应，提高了非线性迭代收敛效率。以病态仿真数据和水下实测数据为例，验证了该方法的数值收敛解优于线性平差估计解，收敛效率优于迭代Tikhonov正则化方法。     

顾及地类转换差异的城市空间扩展元胞自动机模型及应用研究

元胞自动机模型已经成为城市空间扩展模拟研究的重要方法之一,并得到广泛应用。然而,现有的城市扩展元胞自动机模型仍存在不足。由于元胞状态设置较为简单,从而使模型转换规则中对不同用地类型向城市用地转换的差异与强度考虑不够。基于此本文在元胞自动机模型的框架下,设计了多元结构的元胞状态及转换规则,提出了顾及地类转换差异与强度的城市扩展元胞自动机模型。在计算非城市用地向城市用地转换的转换概率时,该模型考虑了3个方面的概率：① 地形地貌、经济发展等城市发展的驱动因素对城市用地扩展的影响概率,该概率采用logistics方法进行计算;② 邻域元胞的用地类型对中心元胞转换概率的影响,该概率采用扩展摩尔型方法进行计算;③ 不同类型的非城市用地（本研究中包括耕地、林地和裸地3种类型）向城市用地转换的强度,该概率由模拟基期土地利用数据与目标年份土地利用数据的叠加,得出不同类型的非城市用地在此时间段内向城市用地转换的规模,进而确定不同类型的非城市用地向城市用地转换的强度。最后,将以上3种概率的乘积作为元胞转换的概率。通过转换概率与转换阈值的对比判断中心元胞是否在下一个阶段转换为城市用地。经过迭代计算,不断增加城市用地元胞的数量。当模拟城市用地的结果与目标年份的城市用地规模差值在一定的范围内时停止模拟,得出最终结果。模型构建完成后,本文以长株潭城市群核心区为例进行了模拟实验。以2001年该地区的土地利用数据为基期数据,模拟2010年该地区的城市用地规模和空间分布。研究结果表明,根据本文提出的模型模拟的城市扩展结果与真实数据相比具有较高的一致性。模拟结果正确率达到68.66%,比基于传统logistics回归的元胞自动机模型的模拟精度提高了4.25%,Kappa系数为0.675。该模型较好地模拟了长株潭城市群核心区城市扩展,在城市空间扩展模拟中具有较好的适应性与有效性。     

弧形海岸裂流的数值模拟研究

弧形海岸波浪产生的裂流严重危害人类活动，但是目前对其特征缺乏充分认识。本文对Haller物理模型实验和三亚大东海的数值模拟表明FUNWAVE模式具有较好的裂流模拟能力。基于该模式进行了多种弧形海岸条件的裂流数值模拟，给出裂流的一些特征：（1）海岸弯曲度增大，裂流增强；（2）海岸坡度对裂流有比较大的影响，太陡或太平缓的海岸不利于形成裂流；（3）海岸尺寸减小，裂流减弱；（4）波高和波周期增大，裂流增强，但是对于某些海岸而言，0.4m波高可能就存在危害比较大的裂流。     

张掖盆地地下水硝酸盐污染与人体健康风险评价

地下水是张掖盆地的重要水资源，其硝酸盐污染尚未得到足够重视。对张掖盆地2004、2015年地下水硝酸盐浓度进行了系统分析，并采用美国环境保护署（USEPA）推荐的健康风险评价模型评估了地下水硝酸盐的健康风险。结果表明：自2004年以来张掖盆地地下水硝酸盐污染日趋严重。2015年硝酸盐浓度最高已达到283.32 mg·L^-1，17.61%的采样点硝酸盐氮浓度超过GB5749-2006《生活饮用水卫生标准》中饮用地下水限量值（20 mg·L^-1）。研究区人群经皮肤接触途径摄入硝酸盐的健康风险在可接受水平，而饮水摄入硝酸盐的健康风险较高，总风险中饮水途径引起健康风险的贡献率占99.40%，远大于皮肤接触途径。儿童经饮水摄入和皮肤接触两种途径的健康风险均显著高于成人，分别为成人的1.544倍和1.039倍。32.39%的采样点地下水硝酸盐对儿童的健康风险超出了可接受水平，14.79%的采样点地下水硝酸盐对成人的健康风险不可接受。甘州区城区、临泽县北部边缘及高台县城区周围硝酸盐浓度最高，这些区域内所有人群都面临硝酸盐引发的高健康风险，其余区域硝酸盐引发的健康风险相对较低。     

典型深街谷内树木空间配置对行人呼吸高度处气流的影响

为揭示树木的不同空间配植方案对行人呼吸高度气流的影响,本文将树木视为均匀多孔介质,通过附加源项法从空气动力学角度用CFD模拟了H/W=2的典型深街谷几何内4种树木配植情景,实验表明,不同空间配置下树木对街谷内行人呼吸高度处局地气流的影响强弱在空间分布模式上差异悬殊：① 均匀种植的树木对街谷内行人呼吸高度的气流起到阻碍作用,不均匀种植则有效提升街谷的整体流速。4种空间配植方案下树木对气流的影响程度不同,阻碍作用从大到小的顺序为均匀间距8 m（Spa8m）>均匀间距6 m（Spa6m）>均匀间距20 m（Spa20m）>不均匀配植（Non-uniform）;对应的平均气流增强指标顺序为$\bar{D}_{spa8}$（-19.31%）<$\bar{D}_{spa6}$（-16.14%）<$\bar{D}_{spa20}$（-10.73%）<$\bar{D}_{non-uniform}$（1.25%）。② 对比不均匀和均匀的种植方案,不均匀植树的街谷内部行人呼吸高度的气流流速比其对照案例（均匀植树Spa8m方案）整体增强了106.49%。街谷中部不种树,在街谷两端配置树木并预留足够的自由空间的不均匀植树方案,能够让角涡渗入街谷中部,促使街谷内部的垂直漩涡和两端的水平角涡运动,增强湍流和垂直交换,有效减少了街道两端“风口效应”和街道中部“风影效应”的区域,改善了整个街谷行人呼吸平面的风环境。④ 合理空间配置的树木能够改善街谷内部的行人风环境。街谷内行人呼吸高度处的气流对局地条件很敏感,树木的局部配置（空间簇集、密度）将引起强烈的空间变化。在既有城市建筑布局条件下,如何通过谨慎的景观设计,利用树木等城市绿化措施有效地改善城市的行人风环境,缓解污染扩散、疾病传播等问题,本文的方法可提供一定的参考。     

非线性波消波及波场分布研究

本文基于雷诺平均N-S方程,并结合RNG k-e方程建立了粘性数值波浪水槽,对不同波陡、不同相对水深、不同相对波高的非线性规则波的阻尼消波问题和波场分布进行研究。文中提出了两种描述消波区内部阻尼变化的阻尼函数,分别适用于小波陡情形和高波陡情形。研究结果表明,小波陡组消波区可设为一个波长,阻尼系数取10~4～10~5即可满足消波要求,计算结果与实验结果及造波理论吻合良好;高波陡组消波区可设为两个波长,阻尼系数取10~4～10~5亦可满足消波要求,计算结果与实验结果吻合良好。此外,当波陡较小时,波场内反射情况的小幅改变即可对整个波场造成影响,特别是当水深较浅时这种影响极为明显,需谨慎考虑。当波陡较大时,水波能量较高,整个波场沿水波传播方向可观测到明显的衰减现象,在具体试验中需进行考虑。     

察尔汗盐湖西北部典型采区地下卤水赋存状态的地面核磁探测与数值模拟研究

察尔汗盐湖地下晶间卤水蕴藏了宝贵的盐湖矿产资源,如何确定地下卤水的赋存规律以制定科学合理的卤水开采和补给方案,正成为盐湖资源大规模可持续利用亟需解决的重要问题。本文以察尔汗盐湖西北部的典型采卤区作为研究区域,1)利用新型的地面核磁共振找水仪,通过在研究区布设两条十字交叉的测线,探测并反演解译测区地下60米深度内的卤水静态赋存形式,2)结合OpenGeoSys(OGS)多物理场耦合地下水数值模拟软件,对测区地下卤水的动态运移形式进行数值模拟计算。地面核磁探测结果表明研究区含水量整体偏低,最大含水量约0.8%,平均含水量约0.4%。测点结果给出的含水层位分界面与邻近钻孔的地层层位分界面的一致性较高,浅部的层位误差在10%以下,达到了0.5米的精度。考虑到测区南北两边均有采卤渠正在汲取地下卤水,且在测区中心存在两条废弃的采卤渠,OGS地下水数值模拟研究结果表明采卤渠造成测区浅部卤水较大空缺,与研究区浅部含水量总体偏低的探测特征一致。研究结果表明地面核磁共振技术能准确探测卤水当前的赋存状态,地下水运移模拟技术可以深入认识采补活动引起卤水赋存状态的变化规律。探测和模拟结合的联合研究是确定卤水赋存规律的重要手段,研究成果可为盐湖资源合理开采高效利用提供基础资料和相关科学支撑。     

Hamilton系统数值计算的新方法

系统地介绍了近年来对Ｈａｍｉｌｔｏｎ系统数值计算新建立的辛算法和线性对称多步法，并对它们在动力天文中的应用作了一简要回顾。     

矩形虾池环流的流动特点与流动阻力

对矩形虾池在水车式增氧机作用下产生的环流特点进行理论分析，结果表明，在转角处会形成涡旋，边壁附近在一定条件下会产生回流；利用ANSYS软件的FLOTRAN CFD进行数值计算，获得与理论分析一致的结果，与现场实际情况吻合。并进一步对虾池阻力进行分析，提出了一些减阻措施。