极坐标系有限差分中起伏地表边界条件处理

有限差分算法是地震学中重要的算法,在直角坐标系下同位网格有限差分中使用牵引力镜像方法,可以高效准确地处理起伏地表边界条件.当研究区域-全球尺度问题时需要考虑地球曲率影响,此时选择极坐标系更加直观方便,但已有方法无法在极坐标系下准确计算起伏地表影响.本文在极坐标系有限差分中引入贴体网格和牵引力镜像方法处理起伏地表边界条件,并在多个算例中验证算法正确性和适用范围,证明牵引力镜像法在极坐标系有限差分中有效.

     

基于交叉梯度结构约束的大地电磁主轴各向异性并行三维反演

各向异性介质对大地电磁观测数据的影响往往不可忽略,因此需要提高大地电磁各向异性三维反演的可靠性和有效性.为了满足大地电磁各向异性三维反演的需求,本文研究了一种基于交叉梯度结构约束的大地电磁主轴各向异性并行三维反演算法.根据大地电磁平面波理论假设,正演方程采用背景场与二次场分离的计算方式,二次场利用交错网格有限差分法求解.由于各向异性反演的多解性,本文将各向异性介质简化为主轴各向异性,并在此基础上进一步采用有限内存拟牛顿LBFGS法实现三维各向异性反演.为了提高各向异性反演的分辨率,反演目标函数中引入交叉梯度项,利用先验的结构信息,对三个方向的电阻率参数进行结构约束,最终的反演进一步利用MPI（Message Passing Interface,消息传递接口）技术实现分频并行计算,测试结果显示并行接近线性加速比.

     

地-空瞬变电磁法电阻率成像研究与应用

地-空瞬变电磁法（Semi-airborne transient electromagnetic,SATEM）凭借适应能力强、探测深度大、实时性强等特点,适用于湖泊、沼泽、山区等地形复杂地区观测.SATEM接收数据量大、精度要求高,传统成像方法基于地面视电阻率计算未考虑飞行高度,开展SATEM的高精度快速成像研究对实际运用有重要作用.针对SATEM中电性源发射装置感应电流分布,接收装置设置于空中以及两者不在相同平面的特点,本文考虑横向、纵向感应电流差异,飞行高度对视电阻率影响等因素,提出电性源SATEM电阻率成像方法.首先,基于电性源瞬变电磁场感应电流分布特征及各分量空间分布、扩散特性分析,推导出基于均匀半空间磁场与感应电压解析解,然后定义了电性源地-空瞬变电磁法早期、晚期、全期视电阻率;最后根据对感应电流分布分析定义了横向分量和纵向分量的成像深度,并研究了成像深度受飞行高度、发射磁矩、偏移距、偏移角度的影响;对煤炭采空区的实测结果表明,利用本文提出的电阻率成像方法可以更加精确定位目标体位置.

     

一种利用Nyström离散与FFT快速褶积的散射地震波并行计算方法

利用数值方法解Lippermann-Schwinger (L-S)方程的主要困难在于系数矩阵存储和线性方程组求解.这主要是因为L-S方程的积分部分是一个空间褶积,在离散后将导致一个满秩矩阵,进而形成一个大型或超大型代数方程组.因此,在利用L-S解决地震波散射问题时,一般是利用散射级数法而非数值方法.然而,散射级数法的计算精度和收敛性强烈地依赖于速度扰动的强度,而克服这种依赖性的一个可能的途径就是对现有的数值方法进行改进或是建立新的数值求解方案.在这种思想指导下,首先对L-S方程进行改写,得到一个与原L-S方程等价的积分方程(等价L-S方程).然后,对等价L-S方程进行逐点归一化处理,并利用Nyström法对经归一化处理的等价L-S方程(归一化等价L-S方程)进行离散,并用FFT计算空间褶积.之所以这样选择是由于归一化等价L-S方程经Nyström法离散生成的系数阵为一个Toeplitz阵,可利用其Toeplitz性质降低存储空间;而FFT可以将矩矢空间褶积转化为乘积,且积分核部分只要计算一次即可.进一步,为节约正演计算时间,设计了进程级和线程级相结合的MPI+OpenMP并行模式.数值试验表明,与传统的积分方程数值算法相比,利用等价L-S方程、Nyström离散和FFT快速褶积的计算方案可极大地降低存储需求,进而在保证精度的同时提高计算效率.

     

多尺度阵列嵌套组合反演宾川气枪源区横波速度结构

本文基于多尺度阵列嵌套组合的方式,利用频率-贝塞尔变换法（Frequency-Bessel,F-J方法）提取背景噪声面波频散信息,通过多个阵列融合的频散曲线反演得到宾川气枪发射台周边不同深度的横波速度结构.结果显示：浅层一阶面波频散信息的加入,使得基阶反演结果更加收敛,反演深度加深到8 km;深度在8 km以下的结构的研究利用多尺度阵列（密集台阵-宾川气枪台网-云南区域地震台网）嵌套组合的方式,面波基阶低频信息从0.55 Hz拓宽到0.008 Hz,使横波速度结构的反演深度显著增加,同时对反演过程提供约束,使得70 km深度以上的横波速度更收敛.由此本文所得的横波速度结构为该区地下结构的探测提供基础,多尺度阵列嵌套组合频散谱的研究方式也为以后区域结构的研究提供一种新的方法和思路.

     

任意各向异性介质三维非结构谱元法直流电阻率正演模拟研究

各向异性是地电异常解释中不可忽视的因素,广泛存在于裂隙或层理发育的地质环境中.本文针对任意各向异性条件下直流电阻率法三维正演问题进行研究,结合非结构谱元法建立模拟算法,充分利用谱方法的指数收敛性以及非结构有限元对地形和复杂异常体刻画能力,提高计算精度和效率.通过灵活的四面体网格剖分和高阶谱插值,实现了复杂介质任意各向异性模型电阻率响应的高精度数值模拟.我们首先通过层状各向异性模型验证本文非结构谱元法的计算精度,进而我们以半空间中立方体模型为例分析各向异性对电阻率响应的影响特征,并通过计算针对不同各向异性参数的视电阻率极性图,探究地下介质各向异性特征识别方法.最后,我们针对典型的山脊模型计算和分析存在地形效应条件下各向异性电流场分布及视电阻率特征.模型计算结果表明基于四面体网格的谱元法模拟带有复杂地形和异常体的任意各向异性模型具有很高的计算精度.本文的研究成果将在推进电阻率方法用于解决裂隙及层理等环境和工程地质问题中发挥积极作用.

     

在役公路梁式桥综合震害等级预测研究

为客观预测在役公路梁式桥综合震害状况,考虑在役桥梁在运营期存在的病害问题,从压力和承压两方面建立在役公路梁式桥综合震害预测评价指标体系.以桥梁作为承灾体,建立在役公路梁式桥综合震害物元可拓模型,运用熵权法进行赋权,确定桥梁的综合震害状况.以一座在役梁式桥为例,运用上述模型确定算例的综合震害状况.研究结果表明,该桥的综合...     

基于耐震时程法的进水塔结构抗震性能分析

进水塔结构作为水电站的进水建筑物,其抗震性能对水电站的安全运行有重要意义.以沙牌水电站的高耸独立进水塔结构为研究对象,应用3条耐震时程曲线对其进行抗震性能分析与评估,并将计算结果与天然地震动下的增量动力分析结果进行比较,对比分析进水塔塔顶位移、基底剪力及地震损伤历程.实例研究表明:耐震时程法可以较好地分析结构的抗震响应...     

瞬变电磁法在彬州市某煤矿采空区勘查中的应用

采用地面瞬变电磁法对拟建搬迁安置小区是否存在煤矿采空区进行勘查,运用未充水采空区、充水采空区分别在视电阻率断面图上呈现相对高阻、低阻异常的地球物理标志,结合地面地质灾害调查成果,较为准确的确定了彬州市某煤矿采空区具体位置,并推断塌陷过渡带范围。     

采石场大气污染物源强分析研究

采石场环境影响评价的关键在于生产过程中产生的大气污染TSP、PM10的影响,而TSP、PM10的影响又取决于其产生的源强的大小及扩散模式.以襄樊市襄阳城南岘山两处采石场为例,采用模拟实验、实测法对采石场大气污染物扩散模式、扩散参数及源强的确定进行分析研究.结果表明:在通常气象(平均风速)条件下,各类稳定度的TSP小时平均浓度在下风向距离800 m以上均较低.基本上与实测值的TSP值相吻合.在D类稳定度情况下,距离污染源800~900 mTSP最大小时平均浓度为0.13~0.11 mg/hm3.