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

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

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

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

基于ROMS和4DVAR的沿轨与网格化SSH数据同化效果评价

Remote sensing products are significant in the data assimilation of an ocean model. Considering the resolution and space coverage of different remote sensing data, two types of sea surface height(SSH) product are employed in the assimilation, including the gridded products from AVISO and the original along-track observations used in the generation. To explore their impact on the assimilation results, an experiment focus on the South China Sea(SCS) is conducted based on the Regional Ocean Modeling System(ROMS) and the four-dimensional variational data assimilation(4 DVAR) technology. The comparison with EN4 data set and Argo profile indicates that, the along-track SSH assimilation result presents to be more accurate than the gridded SSH assimilation, because some noises may have been introduced in the merging process. Moreover, the mesoscale eddy detection capability of the assimilation results is analyzed by a vector geometry–based algorithm. It is verified that, the assimilation of the gridded SSH shows superiority in describing the eddy's characteristics, since the complete structure of the ocean surface has been reconstructed by the original data merging.     

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

为揭示树木的不同空间配植方案对行人呼吸高度气流的影响,本文将树木视为均匀多孔介质,通过附加源项法从空气动力学角度用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地下水数值模拟研究结果表明采卤渠造成测区浅部卤水较大空缺,与研究区浅部含水量总体偏低的探测特征一致。研究结果表明地面核磁共振技术能准确探测卤水当前的赋存状态,地下水运移模拟技术可以深入认识采补活动引起卤水赋存状态的变化规律。探测和模拟结合的联合研究是确定卤水赋存规律的重要手段,研究成果可为盐湖资源合理开采高效利用提供基础资料和相关科学支撑。     

一种有发展前景的隧道及井巷超前预报技术——电磁导弹探测技术

对于隧道和煤矿井巷的安全掘进 ,超前预报显得特别重要。但由于掌子面前的特殊环境限制及多种干扰的严重影响 ,给隧道及井巷超前预报技术的研究带来许多困难。在瞬变电磁理论的基础上 ,引入新的电磁传播理论———电磁导弹技术 ,作为对该技术的初步探索 ,文章论证了该技术在隧道超前预报中应用的可能性 ,并提出了研究的主要内容和发展方向以及可能取得的新的技术突破和技术支撑     

FEASIBLE STUDY ON THE INTEGRATION SYSTEM FOR THE SPACE MONITORING OF MAJOR EARTHQUAKES AND VOLCANOES IN TERRESTRIAL LAND

With the rapid development of space technology, earth observation technology and sky observatory technolo-gy, they have played a more and more important part in monitoring and predicting of earthquakes and volcanoes in the terres-trial land. In recent years, the related agencies have done the experiments and researches on monitoring and predicting ofearthquakes and volcanoes in the forewarning period by means of many approaches, such as satellite thermal infrared re-mote sensing (TIRS), Global Positioning System (GPS), differential interferometric synthesis aperture radar (D-INSAR),astronomical time-latitude residual anomaly, and Geographic Information Systems (GIS), etc. A quite large number of re-search foundation has been built in the fundamental theories and application methods. The experiments and researcheshave shown that these technology is efficient methods for high frequency crust movement. If the existed separate scientificforces and results are possibly assembled together to form a more complete integration monitoring system with the combina-tion of space, sky observation, ground, deep geology and macro anomaly, it will come into a new stage of monitoring andpredicting of earthquakes and volcanic eruptions.     

Ab Initio Prediction of ^29Si-NMR Chemical Shifts

The ability of several ab initio models to predict experimental 29Si-NMR chemical shift is examined. The shielding values of trimethylsilyl chloride (A), t-butyldimethylsilyl chloride (B) and allyltrimethylsilane (C) are calculated by GIAO, CSGT and IGAIM methods, using HF/6-31G^*, B3LYP/6-31G^* , HF/6-311 G^**, B3LYP/6-311 G^** and MPWlPW91/6-311 G^** models respectively. The 29Si chemical shifts calculated by GIAO method using HF/6-311 G^** model are highly in agreement with those obtained experimentally. All of the models above reproduce the trends of chemical shifts in all cases studied, suggesting that the models are of practical value.     

Hamilton系统数值计算的新方法

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