联合多代卫星测高和多源重力数据的局部大地水准面精化方法

本文研究了基于泊松小波径向基函数融合多代卫星测高及多源重力数据精化大地水准面模型的方法.分别以沿轨垂线偏差和大地水准面高高差作为卫星测高观测量,研究了使用不同类型测高数据对于大地水准面建模精度的影响.针对全球潮汐模型在浅水区域及部分开阔海域精度较低的问题,引入局部潮汐模型研究了不同潮汐模型对于大地水准面的影响.数值分析表明:相比于使用沿轨垂线偏差作为测高观测量,基于沿轨大地水准面高高差解算得到的大地水准面模型的精度更高,特别是在海域区域,其精度提高了2.3 cm.由于使用沿轨大地水准面高高差作为测高观测量削弱了潮汐模型长波误差的影响,采用不同潮汐模型对大地水准面解算的影响较小.总体而言,船载重力及测高观测数据在海洋重力场的确定中呈现互补性关系,联合两类重力场观测量可以提高局部重力场的建模精度.     

基于局部加密非结构网格的海洋可控源电磁法三维有限元正演

本文基于非结构网格实现了海洋可控源电磁法三维有限元正演模拟.该算法采用完全非结构网格剖分,可以模拟任意起伏地形和复杂地电模型.为了避免场源的奇异性,采用一次场/二次场分解算法,一次场由基于Schelkunoff势函数的一维解析公式得到.为了提高算法的精度和效率,采用对测点附近单元和异常体区域进行体积约束加密的方法,实现了非结构网格的局部加密.一、二维模型计算和分析表明,本文采用的局部加密方法能够明显地改善算法的精度,最大相对误差基本在1%以内.对三维模型计算及对比分析,说明了该算法对三维可控源电磁正演的实用性.复杂海底地形模型的正演模拟表明,海底地形对电磁场的影响很大,在进行海洋可控源电磁资料解释时,地形的影响有必要考虑在内.     

利用样本排序方法比较化探异常识别模型的效果

地球化学异常的有效识别是化探找矿成败的关键环节.利用样本排序方法对各种化探异常识别模型的处理效果和优劣性进行了比较.以云南个旧及其周边地区铜元素水系沉积物为例, 应用元素含量、异常衬度、局部奇异性分析三大类方法对数据进行了处理, 对化探异常指示变量的排序值进行了3个方面的对比: (1) 在高背景区和低背景区样品的排序特征; (2) 有矿出现位置处样品的样本排序特征; (3) 累积面积(等效于上侧样本排序值) 不同分组所圈定的异常范围和矿床(点) 之间的空间相关性.结果表明, 局部奇异性分析方法较按含量高低的全局性方法对弱缓异常识别更为有效, 也相对优于滑动衬值.按奇异性指数基于证据权法圈定的异常远景区更具预测意义, 显著优于用元素含量值所圈定的异常范围.局部奇异性方法原理清晰、方法简便、可操作性强, 在地球化学异常识别中可以用其替代滑动衬值方法.      

虚拟场地上杭州地铁隧洞开挖数值模拟

根据杭州市4 000多个综合钻孔资料,利用MapX组成的计算机虚拟场地,并在其上集成了ABAQUS等成熟的计算机软件,对正在建设的杭州市地铁1#线上任意处,在计算机上点击鼠标后,即进行隧洞、衬砌的应力、应变分析与地面沉降计算。由于分析、计算能随时随地地快速进行,故能同监理数据对照。对保证工程安全施工有较大的裨益。     

基于交替变量局部梯度法对土质边坡的三维稳定性分析

临界滑动面的获得是土质边坡稳定性分析中一个必要的过程,然而现有分析方法难以保证该临界滑面的准确性,而且部分分析方法在计算过程中会陷入局部最小值陷阱。基于单变量方法和最速梯度法,提出了交替变量局部梯度法,该方法能够突破最优化方法应用在边坡稳定性分析上的瓶颈。首先,通过斯宾塞极限平衡方法和网格搜索法（或者其他确定临界滑面的方法）,获得边坡初始椭球状临界滑面;其次,在该滑面上布置若干节点作为变量,目标函数为安全系数Fs关于空间节点坐标Zi的方程,为了使目标函数快速降低,沿负梯度方向循环优化每个节点,当满足一定精度要求时,即可获得三维土质边坡非圆弧状临界滑动面;最后,通过算例证明该方法可行,计算结果可靠。     

汶川地震地表破裂影响带调查与建筑场地避让宽度探讨

5.12汶川Ms8.0级地震沿龙门山中央断裂带形成长度约270km的地表破裂带,破裂带垂直和水平位移一般为1.0～5.0m,最大可达9.8m。在地表破裂带附近,建筑物破坏严重。笔者系统调查了龙门山中央断裂带地表破裂的垂直和水平位移以及地震地表破裂影响带宽度,并采用物探测试方法对部分影响带宽度进行了验证。统计结果表明,汶川地震地表破裂影响带宽度主要集中在16～60m,影响带宽度(D)与垂直位移(H)具有较好的线性关系,即D=10.11H+16.0。对于逆冲性质为主的地震地表破裂而言,上盘影响带宽度与下盘的比值一般为3:1～2:1,据此可确定地表破裂带上、下盘的避让带宽度。笔者根据现场调查并借鉴国外经验提出,上盘避让宽度下限一般不应低于15m,下盘避让宽度下限不应低于10m。本项研究对于汶川地震灾后重建场地选址具有重要的实际意义。     

西昌断陷盆地场地地震动差异分析

以汶川地震中的非发震断层安宁河断裂带周边西昌地区所获取的加速度记录为依据,研究了这些记录在峰值、频谱之间的差异,并通过这些地区的局部场地条件差异进行了初步解释。在此基础上通过数值模拟方法,利用基于显式有限元和局部透射人工边界的二维有限元模型对该地区的地面运动进行了模拟。数值试验结果表明：这些简单的理想化模型对这一地区地震动的本质特征进行重现,并对产生地面运动差异的机理进行合理解释。通过分析,初步揭示了由深部断层构造所控制的断陷盆地对地震动影响显著。     

局部地形对出平面运动谱特性的影响分析

本文利用显式有限元有限差分方法进行了脉冲和实际地震动输入下的粘弹性场地地形影响的分析计算,研究了坡地地形在高宽比不同、体波入射角度不同时出平面运动谱特性的差异。结果显示:地形的特征尺寸和入射角度不仅对地震动的峰值有影响,而且对地震动的谱特性有较大的影响,其中入射角度的影响尤为显著。另外,研究还表明,基于脉冲输入研究粘弹性地形谱特性的影响仍不失为一种有效的方法。     

单向传播子的局部平面波分解以及在定向照明分析和目标特征成像中的应用

本文基于波场的局部平面波分解,得到空间和方向都具有局部性的局部平面波小波束,并由此对Fourier有限差分传播算子进行局部平面波分解.这种用小波束分解的单向传播算子,可被用来对采集系统进行定向照明分析和采集倾角响应估算,以分析采集系统布局和上覆结构对面向目标的叠前偏移成像质量的影响.进而,利用局部平面波成像矩阵及其成像相册,根据目标结构的主倾角方向,抽取并叠加部分局部平面波像,可以进行目标结构定向化特征成像.最后,对具有速度强横向变化上覆结构的陡倾角目标断层的复杂二维SEG-EAGE盐丘模型数据,进行采集系统目标照明分析和目标结构定向化特征成像.     

硬场地实测记录下几种土层地震反应分析程序可靠性对比

以日本KiK-net强震观测台网中硬场地井下记录为样本,对传统等效线性化方法LSSRLI-1（频域）、SHAKE2000,时域非线性方法DEEPSOIL和频域一致等效线性化方法SOILQUAKE等几种计算程序在硬土场地地震反应分析中的可靠性进行对比检验。检验工况包括KiK-net井下台网中地表峰值加速度不小于0.05g的水平硬场地的总计344台次的加速度记录,涉及5个台站,土层厚度6~50m,地表峰值加速度范围0.050~0.805g。结果表明:认为硬土场地以往计算方法能够体现土层放大的认识是片面的,在中强地震动情况下,现有流行方法计算出的地表响应会偏小,强地震动下会严重偏小,会给工程抗震设计提供偏于危险的输入;硬场地中烈度8度以下（地表PGA在0.19g以下）,SHAKE2000和DEEPSOIL、LSSRLI-1（频域）计算结果与实际记录差距可以接受,但烈度8度以上,计算出地表响应均较实际记录明显偏小,且随地震动强度增加差距急剧增大;硬场地中,无论何种烈度,SOILQUAKE16计算的地表加速度响应与实测相当,可体现出土层放大作用。