P波斜入射角对沉管隧道地震响应的影响

为了研究P波斜入射对沉管隧道地震响应的影响,以港珠澳大桥沉管隧道为工程背景,考虑上覆海水与海床、沉管隧道之间耦合作用,采用粘弹性边界和等效力的地震荷载输入方式,利用ADINA软件建立三维有限元模型进行地震响应分析。分析入射角为0°、20°、40°、50°、60°时P波对沉管隧道环向应力峰值（正应力峰值、剪应力峰值）和位移峰值的影响,结果表明:入射角为40°时,沉管隧道应力峰值最大;入射角为0°—40°时,隧道的应力峰值逐渐增大,入射角为40°—60°时,隧道的应力峰值逐渐减小;隧道截面4个转角处及隔墙与顶板、底板的连接处为隧道剪应力峰值最大处;隧道截面左侧剪应力峰值远大于右侧;隧道顶板正应力峰值最大,顶板的正应力峰值大约为底板的2倍;隧道截面左侧位移峰值远大于隧道截面右侧。     

多尺度裂缝储层地震预测方法研究

常用的叠前裂缝预测方法并不能明确区分不同尺度的裂缝,本文研究并提出基于裂缝密度的大、中尺度裂缝定量预测方法。该方法通过基于Curvelet变换的多方向相干分析技术进行大尺度裂缝(大于1/4波长)预测,通过地震方位各向异性方法和叠前衰减地震属性(如频率衰减梯度)进行中尺度裂缝(1/4～1/100波长)预测,然后将得到的两种不同尺度裂缝的发育强度以裂缝密度为基准进行融合,综合地评价不同尺度裂缝。将该方法运用到裂缝储层地震物理模型中,结果表明,该方法不但克服了叠前衰减属性的方位各向异性预测中尺度裂缝的同时,预测出断层和大尺度裂缝的密度分布不连续的问题,还能对裂缝的尺度进行区分,较准确的定量预测裂缝发育带。     

基于Google Earth的气象多参数综合显示系统

利用Google Earth平台,将天气过程中的雷电、红外云图、MICAPS高空图以及TMPA降水资料,通过编写KML/KMZ文件,建立了气象多参数综合显示系统。该系统平台的建立,可在三维空间中真实动态地再现天气的发生发展及演变过程。此平台将有利于灾害性天气成灾过程和物理机制的研究。     

1979-2016年祁连山地区大气水汽含量时空特征及其与降水的关系

采用1979-2016年ECMWF1.5°×1.5°逐月再分析资料及同期37个气象站点的降水资料,利用一元线性回归、累积距平、Kriging及IDW（反距离加权）等方法分析了祁连山地区大气水汽含量时空分布特征、降水转化率空间变化规律以及风场分布规律,并对比分析了中国西部不同地区降水转化率的变化趋势。结果表明：（1）1979-2016年祁连山地区大气水汽含量整体呈增加趋势,且季节变化明显。其中夏季是各层大气水汽含量最多的季节,高达329.24 mm,占多年平均大气水汽含量的48.1%。（2）近38 a来,祁连山地区的大气水汽含量呈东南多、西北少的空间分布,且随海拔的升高而逐渐减少,整层大气水汽主要集中在5 000 m以下。（3）祁连山地区的降水转化率从空间上表现出由东向西递减的趋势,说明该地区空中云水资源的开发潜力自东向西逐渐增强,空中云水资源的开发潜力区域差异明显;季风所携带的水汽对其影响区域的降水贡献率较高,西风所携带的水汽则对其影响区域的降水贡献率较低。（4）中国西部地区降水转化率呈向心式递减的趋势,且区域空间波动较大。     

遥感聚合图像的快速拓扑矢量化方法

鉴于当前遥感图像矢量化效率普遍不高,提出了一种遥感聚合图像的快速拓扑矢量化方法。该方法对聚合图像中的聚合对象进行唯一性标号,并将这个标号作为属性赋予生成的矢量多边形,通过矢量多边形的标号间映射关系便可以快速建立复杂岛屿多边形的拓扑包含关系。试验表明,该方法能够精确地完成斑块对象边界追踪及复杂岛屿多边形处理,并显著提高矢量化的速度。     

考虑二期荷载的高填明洞加筋减载土压力计算

随着城市建设的迅速发展,既有高填明洞填土顶面将不可避免地出现新增结构物等二期荷载的情况。为了明确明洞回填完成后新增二期荷载对其顶部垂直土压力的影响,在已有明洞减载结构土压力计算公式的基础上,利用土体沉降弹性理论计算公式建立了土工格栅变形与土体沉降变形的关系,并考虑二期荷载作用,进一步完善了高填明洞减载前后的洞顶垂直土压力计算公式。建立数值分析模型,将数值分析结果与文中公式计算结果进行了对比。结果表明：公式计算与数值分析结果较为接近,表现出相同的规律性;有无减载措施对明洞回填土压力影响较大,填土越高,减载越显著,实际作用在明洞顶土压力越小。新增二期荷载作用下的明洞顶垂直土压力计算,应考虑其对土体沉降变形的影响,而非简单叠加。     

直流电阻率2.5维有限元数值模拟

从稳定电流场满足的基本微分方程出发,推导了地面点源电场电位边值问题对应的等价变分问题.用网格离散积分区域时,以四边形为最小研究单元.针对电场在电源附近衰减快,变化梯度大,数值模拟结果在电源附近位置有较大误差问题,采用边界/奇异校正算法改善模拟结果在电源附近的奇异性问题,同时减少边界条件的影响,提高模拟结果精度.从波数域电位求取空间域电位,利用最优化离散波数进行反付氏变换.模型算例结果验证了该算法的可行性与有效性.     

矿柱稳定性影响因素敏感性分析及其应用研究

为研究矿柱稳定性影响因素的敏感程度,以招金大尹格庄金矿为研究背景,从矿柱载荷、强度、失稳形式及影响因素确定四方面推导出两种矿柱形式的安全系数计算公式,设计采用6因素5水平正交试验对矿柱稳定性影响因素敏感性进行分析,并研究主要影响因素与矿柱安全系数之间的函数关系。研究表明,矿柱宽度、矿体开采深度和矿房宽度3个主要因素对矿柱稳定性影响程度是最为剧烈的;矿柱安全系数随着矿柱宽度地增加呈指数形式递增且递增速率不断加大,随矿体开采深度地增加呈幂函数形式递减且递减速率逐渐变缓,随矿房宽度地增加呈负指数形式递减且降低速率逐渐变慢。利用DPS和Matlab建立的矿柱安全系数回归方程得出保证安全生产的合理矿柱尺寸,矿房宽度应不超过8 m,布置条形矿柱宽度应不小于3.6 m,布置方形矿柱宽度应不小于5.9 m。结合采场矿柱布置实际,应用效果较为理想,可为下中段矿体回采矿块布置提供较好的依据。     

基于有限单元法的二维／三维大地电磁正演模拟策略

对于二维和三维大地电磁正演问题,有限单元法最后形成了一个线性方程组KX=p。方程组中的K是大型稀疏的带状对称复系数矩阵,其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难。不完全LU分解处理的BICGSTAB算法,可用于该线性方程组的求解,并且具有速度快,精度高,稳定性好等优点。为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下,设计了非均匀网格剖分。在程序编制中,因只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间。通过对二维模型和三维模型电磁响应的计算,验证了该算法的正确性。     

一种新型CSAMT观测参数的计算方法

常规的可控源电磁法理论在计算视电阻率公式上,由于主要采用了其电磁场,所以难以直接反映全频域视电阻率的值,也不能直观地显现出地下介质的地质构造。这里采用水平偶极子激发的电磁场,提出了电场的全域精确表达式,并可直接计算出大地电阻率。通过与卡尼亚电阻率计算的对比结果表明：该方法的结果在远区等价于卡尼亚电阻率,在近区和过渡带则明显地改善了卡尼亚电阻率的非波场区场畸变,从而能更好地接近基底的真电阻率,更形象地反映了地下介质的垂向电性变化。