非均匀介质地震AVO响应模拟及分析

地震AVO （amplitude versus offset）技术是一项利用振幅信息研究岩性、检测油气的地震勘探技术。常规方法基于Zoeppritz方程计算模型界面处的反射系数,而实际地下非均匀介质中地震反射特征不仅与入射角度、物性差异有关,还与入射波频率、地层厚度、薄互层结构等因素有关。为此,应用传播矩阵理论充分考虑与这些因素有关的调谐干涉等传播效应,针对实际数据计算非均匀地下介质的高精度合成地震记录,对比Zoeppritz方程、Shuey二项近似方程、Shuey三项近似方程以及传播矩阵算法的模拟效果。研究发现：在小角度入射时Zoeppritz方程、Shuey二项近似方程和Shuey三项近似方程的反射振幅和波形基本一致,大角度入射时Zoeppritz方程与Shuey三项近似方程接近;Zoeppritz算法的模拟结果在小角度入射和浅层情况下与传播矩阵算法差别较小,而在大角度入射和深层情况下与传播矩阵差别较大,说明层间多次波的调谐干涉等传播效应不可忽略。     

基于磁偶极子的频率域电磁系统几何误差分析

基于磁偶极子模型的频率域电磁法仪器在进行大地探测的过程中,由于收发线圈的安装误差、晃动以及材料的几何形变等因素影响,收发线圈的位置及姿态角发生变化,从而产生几何误差,而几何误差的定量分析对频率域电磁法仪器在制作、野外实验操作以及校正方面均具有很强的指导意义。本文根据误差来源的不同,将几何误差分为一次场误差和二次场误差两类;通过引入一次场张量矩阵和二次场张量矩阵,基于均匀大地模型正演得到9种线圈架构接收线圈接收到的一次场及二次场,并以微分的方式定量地分析了9种线圈架构收发线圈位置及姿态角变化引起几何误差的变化特征。对9种线圈架构的几何误差分析可知:在几何形变引入的一次场误差中,仅有PERyz、PERzy线圈架构不包含位置变化的一次项误差,仅有VCA、VCP及HCP线圈架构不包含姿态角变化的一次项误差;在几何形变引入的二次场误差中,PERzx、PERxz线圈架构最大,而VCA、VCP及HCP线圈架构受高度及收发距变化的影响很大,需要实时测量并进行校正,VCP线圈架构受收发距变化以及姿态角变化的影响最小。     

基于改进的散射积分算法的初至波走时层析

初至波走时层析是获取近地表速度结构的一种常用方法.随着采集技术的不断发展,可使用的数据量迅速增多,传统的基于射线追踪和解方程组的地震走时层析成像方法面临着内存占用大、方程求解不稳定等问题.为了解决这些问题,本文基于前人在波形反演研究中提出的一种改进的散射积分算法,提出了一种预条件最速下降法初至波走时层析.该方法无需存储核函数矩阵与Hessian矩阵即可方便地实现目标函数梯度的计算与预条件,且该方法计算效率高、求解稳定、易于并行.数值实验结果表明,该方法可以获得与传统方法精度相当的反演结果,但所占用的内存大幅减小.     

NWC通信台在电离层中激发电磁响应的时变特征

本文利用DEMETER卫星VLF频段电场和磁场频谱数据对DEMETER卫星运行期间2005年至2009年澳大利亚甚低频(Very Low Frequency)通信台NWC发射的通信信号造成的电离层电磁响应的日变化、季节变化及年变化特征进行了统计分析,统计结果表明电磁响应日变化显著,夜间电场强度明显增强可达40dB,磁场变化略小也可为15dB左右,而季节变化不显著,年变化主要受太阳活动的影响,太阳活动越强,电磁响应越小.为解释数据分析结果,对地-电离层电磁波传播过程采用传递矩阵方法进行了模拟计算,模拟结果与数据分析的结果一致.我们认为这种随时间变化的特点可能由250km以下电离层电子密度分布特征导致,因此研究250km以下的电离层电子密度变化可能对寻找地震电离层电磁异常有重要意义.     

北京市长安街沿线的扒窃案件高发区分析及防控对策

本文利用空间分析方法,对北京市长安街沿线的扒窃类案件的空间分布特征进行了研究。首先分析了长安街沿线的犯罪分布状况,发现该长安街沿线共存在3个异常明显的犯罪高发区域,分别对应于西单商业圈(A)、建国门路段(B)和大望路—四惠路段(C),通过将长安街沿线的常住人口与兴趣点(Point of Interest, POI)分布进行比较,发现案件与POI密度较高的区域存在较高的相关性。随后,利用核密度估计方法和时空热点矩阵方法分析了3个异常明显的案件高发区域内的犯罪活动时空分布模式,结果表明：高发区域A的犯罪高发时段为上午10时至下午18时,其中峰值出现在正午12时,犯罪活动主要集中在区域内的一些大型购物场所内;高发区域B和C的犯罪高发时间均为上午6时和下午18时前后,与早晚交通高峰时间比较一致,犯罪活动主要集中在区域内比较重要的公共交通站点附近。最后,根据长安街沿线案件高发区域内犯罪活动的时空规律性特征,提出了犯罪防控措施与建议。     

土地利用空间格局的置换变化与代换变化分析

为揭示土地利用空间格局的动态变化,本文基于土地利用转移矩阵,从类型间的空间转移过程出发,构建了置换变化与代换变化的计算模型,并对其2种变化进行了分析。置换变化是2种类型之间以相同面积发生空间位置的置换过程,且各种类型的数量结构保持不变;代换变化则是所研究类型的数量结构保持不变,而与之转换的其他类型数量上的增加或减少;从而将单一类型的总变化进一步细分为置换变化、代换变化与净变化。以晋江流域1985年与2006年2期土地利用格局的变化为例进行分析,结果显示：在1985-2006年间,晋江流域的土地利用空间格局发生了显著变化;园地、建设用地、未利用地呈现扩张态势,草地与旱地呈现萎缩变化,水域变化不大;这6种地类的交换变化面积相对较少,以置换变化为主。水田与林地则主要表现为空间位置的交换变化,其中,水田呈现出与林地、旱地2种类型之间的置换变化,而林地则主要是园地-林地与林地-草地这2对类型间的代换变化。结果表明,通过置换变化与代换变化的分析,更细化了各种类型之间的动态变化过程。     

基于选权拟合法的电离层电子密度层析重构

附加约束的电离层层析算法是解决电离层电子密度反演中不适定问题的主要方法,为避免此类方法中约束权阵的选取不当对电子密度分布重构产生的不良影响,本文将选权拟合法应用到电离层层析成像技术中.该方法特别设计了依据电子密度空间分布特性构造参数权矩阵的方案.新方法有明确的物理意义,挖掘了隐含的信息量,为解决电离层电子密度反演中由于观测数据的不足等因素引起的不适定问题提供了一种新途径,可以得到符合客观实际的结果.数值模拟实验和实测数据的反演结果证实了该算法有效性、可靠性和优越性.     

固体介质中孤立波的传播及演化特征

由于地球介质中广泛存在断裂、微裂缝等地质现象,实际地震资料中会出现形似孤立波的非线性地震现象。因此,对固体介质中孤立波的研究有利于解释这些非线性地震现象的形成机制。本文基于KdV方程,以雷克子波作为初始条件,采用伸缩子机理构建体力模型,利用有限差分的方法模拟了孤立波的演化过程。理论结果表明,非线性地震纵波可以从雷克子波逐渐演化成孤立波,而且地震波的初始振幅和频散系数对模拟结果也有重要影响。通过与实际资料的对比也能说明这种演化的可能性。同时根据方程系数矩阵中元素带状分布的特征,采用稀疏矩阵的存取方法,可以减小计算内存,提高计算效率。     

考虑吸力变化的非饱和高液限黏土-土工织物界面剪切特性研究

以高液限土为填料的加筋挡墙广泛在南方地区应用,加筋挡墙土体吸力经常随气候变化而发生变化,直接影响筋土界面强度特性,考虑吸力变化的非饱和高液限黏土-土工织物界面强度特性亟待深入研究。利用ZFY-1A型非饱和土应变控制式直剪仪对非饱和高液限黏土-土工织物界面进行控制吸力的剪切试验,探讨吸力循环变化（干湿循环）对筋-土界面剪切特性的影响。研究表明,高液限黏土-土工织物界面的抗剪强度随着净法向应力的增加而增大;随着循环次数的增加,界面强度呈现先升高而后衰减的变化规律。经历3次干湿循环的界面强度低于未经历干湿循环的强度。随着循环次数的增加,界面强度总体上呈降低趋势。吸力增加的速率 随吸力循环次数的增加先增大后减小,有效界面摩擦角 、界面有效黏聚力 、界面综合摩擦系数 、综合摩擦系数比K与此类似。 含水率变化量的最大值随着净法向应力、循环次数的增加而增大。最大剪胀量 随循环次数的增加先增大后减小,随净法向应力的增大而减小。最后通过对数拟合,得出经历若干次基质吸力循环后的界面强度计算公式。     

顶管推进引起施工场地竖向附加荷载分析

基于层状体系解析刚度矩阵理论解,结合5节点Gauss-Legendre求积公式,提出了层状地基中顶管施工正面附加推力、掘进机与土体之间摩擦力以及共同作用力引起的附加荷载计算方法,分析了顶管推进引起的土体竖向附加荷载分布规律,也研究了地基等效均质性、土层力学参数、计算点间距以及顶管埋深等因素对顶管施工诱发附加荷载的影响效应。研究结果表明,掘进摩擦力引起的附加荷载在掘进面前方迅速达到压应力峰值,其量值大小和影响范围均要大于正面附加推力,是顶管施工引起临近地层附加荷载的主要影响因素。此外,层状地基土体参数的改变会对顶管施工扰动地层的附加荷载产生一定影响,地基等效均质性、计算点间距以及顶管埋深等因素对附加荷载大小及分布均存在显著影响。成果可为合理制定顶管开挖对周围土工环境的保护措施提供一定理论依据,也可为其他盾构隧道工程提供一定的理论参考。