混合多位移模式的非连续变形分析法研究

传统的非连续变形分析法（DDA）法采用简单的线性位移模式计算效率高,描述大块体的高阶多项式位移模式在一定程度保留了该特点,并提高了计算精度。近年来流行的耦合有限元、自然单元的DDA法实质上是引入相应的插值形函数构成块体位移函数,计算相对低效,但具有计算更精细、更容易施加边界条件等优点。为结合传统DDA法与DDA耦合法各自的优点,建立了一种同时利用传统DDA法线性位移模式与耦合型DDA法非线性位移模式的混合法。该方法非线性模式主要针对大块体,采用了自然单元插值,缘于其具有一定无网格特征,且效率比有限元高。建立了混合模式下的整体矩阵并推导出接触等因素刚度子矩阵和荷载子向量的具体表达式。该方法建模更加方便合理,计算精度、效率介于线性模式的传统DDA法和非线性位移模式的耦合法之间。通过基本算例验证了混合法的有效性,并给出了节理围岩-隧道衬砌整体分析模型的计算结果,体现了新方法的优越性。     

灰线性加权非等距GM(1,1)形变预测模型

结合加权非等距GM(1,1)模型与线性回归理论,构建灰线性加权非等距GM(1,1)预测模型,并给出对模型预测精度起决定性作用的灰指数v和参数m的优化方法。与加权非等距GM(1,1)模型和线性回归预测模型相比,灰线性加权非等距GM(1,1)预测模型的精度更高,预测有效时间更长,模型的稳定性更好。优化v和m后,灰线性加权非等距GM(1,1)预测模型的实用性、稳定性进一步提高。     

基于改进线性光谱解混和 SCS-CN 模型的广州主城区降雨产流模拟

城市降雨产流特征的定量研究对于城市雨水资源开发利用、科学规划和管理具有十分重要的意义。以广州市主城区为研究对象,提出一种改进的线性光谱解混方法提取遥感影像的不透水面、植被和土壤盖度,并结合SCS-CN 产流模型进行城区降雨径流模拟。首先利用线性光谱解混方法提取 2015-10-18 Landsat 8 OLI 遥感影像的不透水面、植被和土壤盖度;然后利用归一化建筑物指数和归一化植被指数进一步优化解混结果;最后结合优化的不透水面、植被和土壤盖度计算合成的 CN 值,并应用 SCS-CN 产流模型分析研究区在不同降雨重现期的降雨产流分布特征。结果表明：改进的线性光谱解混方法能较好地提高不透水面、植被和土壤盖度的提取精度;不透水面盖度越高的区域,CN 值越高;CN 值对地表径流深模拟有显著的影响;在不同重现期降雨条件下,研究区地表径流深空间分布格局差异显著;总体上,降雨量越大,不透水面盖度越高,形成的地表径流深越深,内涝发生的可能性越大。     

基于断裂力学的岩石切削数值分析探讨

为研究岩石切削过程中刀具的力学特性以及刀具平均峰值作用力的预测,基于断裂力学和有限元理论建立了岩石切削的数值模型。模拟分析岩石对刀具的反作用过程,得到不同力学特性岩石和刀具作用条件下刀具所受的反作用力,结果表明,冲击速度在一定范围内对刀具平均峰值作用力影响较小;岩石体积呈阶跃性递减,刀具作用力大小与崩落岩石碎片体积变化一致。对模拟、实验以及已有理论计算得到的刀具平均峰值作用力进行线性回归分析,研究结果表明,数值方法计算结果与实验结果、Goktan[1]理论计算结果的线性相关系数分别为0.812、0.930,相关性较好,表明所建数值模型用于模拟岩石切削是正确和可靠的。     

电导率连续变化的线源FCSEM有限元正演模拟

为了更好地模拟地下介质连续变化及开展连续介质的反演,对二维电导率分块线性变化的线源频率域可控源电磁法进行了有限元正演模拟,在剖分单元内同时对电场及电导率参数线性插值,使电导率参数在剖分单元之间保持连续变化。首先,提出有限元正演模拟的边值问题及变分问题,并详细论述了有限元的剖分、插值、单元分析及总体合成的各个步骤;其次,采用稀疏存储及基于不完全LU分解的BICGSTAB算法求解复系数方程组,节省了内存并提高了计算速度;然后,对一个均匀半空间模型进行模拟,计算结果表明,低频及高频的有限元数值解都与解析解吻合,证明了算法的正确性;最后,对水平层状模型及垂直断层模型进行正演计算,视电阻率及相位的等值线图均较好地反映出了异常体,说明文中算法能够对电导率连续变化的线源可控源电磁法进行有效地模拟。     

1957～2002年南海—北印度洋海浪场波候特征分析

利用ERA-40海表10 m风场驱动第三代海浪数值模式WAVEWATCH-Ⅲ,得到南海—北印度洋1957年9月至2002年8月的海浪场,并分析其波候(风候)特征.研究发现如下主要特征:(1)该海域的波高波向、风速风向受季风影响显著;(2)北印度洋大部分海域的海表风速呈显著性逐年线性递增趋势,大约0.01～0.02 m/(s·a),南海线性递增的区域则较少,有效波高呈显著性逐年线性递增的区域主要集中在低纬度中东印度洋(约0.003～0.006 m/a)、索马里附近海域(大约0.002～0.005 m/a)、南海大部分海域(约0.002～0.004 m/a),线性递减的区域主要集中在孟加拉湾海域(约-0.002 m/a);(3)Nino3指数与南海—北印度洋的海表风场、浪场存在密切的关系;(4)南海—北印度洋的海表风速与有效波高存在5.2a左右的共同周期,南海的海表风速、有效波高还存在2.0a左右的共同周期,北印度洋的海表风速、有效波高还存在26.0a的长周期震荡.     

用混沌理论提高潮汐预报的准确度

为提高潮位预报的准确性,在具有较长潮汐观测数据的站点,基于混沌理论,对观测值与潮汐模型预测值之差所构成的余水位序列(即误差序列),采用局域线性模型的分析方法,给出可能误差预测,修正模型的预报结果,提高潮汐预报的准确性。所给例子,对预测跨度T=2 h,经局域法修正后,崇武站2007年12月份1个月预测水位统计的RMSE值减少74.7%,厦门站减少60.5%;对T=24 h,崇武、厦门两站RMSE值减小都在50%左右。     

西北太平洋海域海表风速长期变化趋势研究

利用来自ECMWF的ERA-40海表10m风场资料,对西北太平洋海域海表风速的长期变化趋势等进行研究,研究发现:(1)1958-2001年期间,西北太平洋第一岛链以内海域的海表风速不存在显著的逐年变化趋势,第一岛链以外地广阔洋面则基本表现出显著性逐年线性递增趋势,递增趋势约0.005-0.02 m·s1·a-1,呈显著性逐年线性递减的海域主要分布于一些零星海域;(2)近44年期间,西北太平洋海域的海表风速整体上以0.0072m·s-1·a-1的速度显著性逐年线性递增,在1958-1974年期间,海表风速的递增趋势较强,1975-2001年期间,西北太平洋海域的海表风速整体上变化趋势较为平缓,尤其是在1976-1983年期间,海表风速的走势甚为平缓;(3)西北太平洋海域的海表风速不存在显著的逐MAM、逐JJA变化趋势,逐SON和逐DJF则表现出显著的线性递增趋势,逐SON的递增趋势为0.0047 m· s-1· a-1,逐DJF的递增趋势为0.0079m·s-1·a-1;(4)西北太平洋海域的海表风速存在多种尺度的变化周期,具有明显的2.0-2.4年、4.3-5.2年以及26年以上的长周期震荡.     

基于两种数据集的全球SST变化趋势的对比研究

分别利用了来自NOAA和Hadley中心的SST(Sea Surface Temperature)资料,分析了全球海域每个网格点上SST的变化趋势,并将两种资料分析的结果进行比较。研究发现:近150a期间,全球大部分海域的SST以0.003℃.a-1左右的速度显著性逐年线性递增,仅在格陵兰南部海域和两极部分海域呈显著的逐年线性递减趋势;NOAA的SST变化趋势的区域性差异比Hadley的SST明显,但后者变化趋势的强度稍大于前者;不同海域SST的逐年变化趋势由不同的季节主导;在1910年为近150a来SST的最低点,之后缓慢上升,在1942年附近达到1个相对高点,而后略有下降,但很快又恢复上升趋势,上升趋势一直延续至今。     

利用三元组稀疏矩阵技术改进HASM算法——以全球平均气温模拟为例

 采用预处理共轭梯度法(PCG)进行求解的高精度曲面建模(HASM)算法的计算过程需要大量矩阵计算,采用稀疏矩阵方式可以压缩存储空间。三元组稀疏矩阵存储是较为传统的稀疏矩阵存储结构,这种存储结构的稀疏矩阵技术已被广泛使用。根据HASM-PCG的特点,本文通过改进三元组稀疏矩阵的部分计算方式,调整HASM-PCG算法中的部分计算顺序,从而舍弃部分不需要存储的非零元素,提高了计算效率。根据全球1998-2008年的近3000个气象观测台站的气温观测数据,以及全球DEM数据,对20世纪末至21世纪初的11年来5、6、7、8月份的全球平均气温进行数字模拟。数值结果表明,采用改进的三元组稀疏矩阵技术有效地提高了HASM方法的模拟效率。