地球附近第23太阳活动周磁云和非磁云ICME的对比统计
摘要：行星际日冕物质抛射（ICME），作为影响地球空间天气的重要源头之一，根据其磁场结构特点可分为磁云（MC）和非磁云ICME两个子集.本文对第23周的磁云和非磁云ICME结构及其地磁效应进行对比统计研究.第23周ICME事件总数为317个，其中磁云占ICME比例为33.75%，非磁云ICME占66.25%.统计结果表明，非磁云ICME数与太阳黑子数呈现出非常好的正相关性，而磁云与太阳黑子数的这种相关性并不明显.相反，磁云占ICME的比率与太阳黑子数呈现出一定的反相关性.对磁云与非磁云ICME引起的地磁暴的比较研究表明：磁云及其鞘区引发的地磁暴平均水平要高于非磁云ICME及其鞘区.磁云和非磁云ICME的磁场强度、南向磁场强度和传播速度整体上都随地磁暴水平提升而增加.对磁云与非磁云ICME参数的进一步对比分析表明，磁云及其鞘区的平均磁场强度和南向磁场分量平均值都明显要比非磁云ICME的大；而二者的等离子体温度、密度和速度平均值相差并不明显.地形/均衡重力场模型构制及其拟稳分析
摘要：重力空白区数据填补的一个主要方法是基于地壳均衡理论进行的，该方法亦用于EGM系列模型的构建中.本文研究了地形数据在构制地形/均衡重力场模型中的应用，分析了补偿深度对Airy位模型和面凝聚位模型的影响，给出二者的最佳补偿深度分别为50 km和40 km.以纯卫星重力模型为参考，后者在前120阶的精度要高于前者，但在121~250阶的精度较低，组合模型精度高于单一模型精度.对地形/均衡地球重力场模型进行了EGM2008拟稳分析，研究了不同分辨率基准的拟稳效果，分析表明：30＇分辨率的拟稳基准所得拟稳模型对应的阶方差与参考阶方差曲线直到360阶都有较好的一致性，以EGM2008为基准，其相对累计大地水准面高误差在140阶时为6.83cm，相对累计重力异常误差在220阶时为1.10 mGal.秦岭造山带与沉积盆地和结晶基底地震波场及动力学响应
摘要：秦岭—大别造山带横贯中国大陆中部，并将我国东部分为南北两部；即华北克拉通和扬子克拉通.在南、北相向运动力系驱动下构成了一个极为复杂的复合、叠加构造带、成矿带和地震活动带.同时导致了该地域异常变化的沉积建造和强烈起伏的结晶基底.然而对它们形成的地球物理边界场响应，岩相和结构的异常变化尚不清晰，特别对盆山之间的耦合响应更缺乏深层动力过程的理解.为此本文通过该区榆林—铜川—涪陵长1000 km剖面的地震探测和研究结果提出：（1）沉积建造厚度变化为4~10 km，结晶基底起伏强烈，幅度可达4~6 km；（2）一系列基底断裂将该区切割为南鄂尔多斯盆地和秦岭北缘前陆盆地、秦岭—大巴造山带和南缘前陆盆地与东北四川盆地，其中前陆盆地为秦岭北渭河盆地和秦岭南通江—万源盆地；（3）秦岭造山带是北部华北克拉通向南推挤、南部扬子克拉通向北推挤下隆升的陆内山体，并构筑了其南、北前陆盆地；（4）秦岭造山带的南、北边界并非是一条边界断层，而应是包括前陆盆地在内的组合界带；（5）秦岭与大巴弧形山系源于同一深部结晶基底，即同根生.这一系列的新认识对深化理解秦岭—大巴造山带形成的深层动力过程和演化机理及厘定扬子克拉通的真实北界具有极为重要的意义.利用Rayleigh波相速度和群速度联合反演青藏高原东北缘S波速度结构
摘要：利用青藏高原东北缘地区固定和流动地震台网2007年8月到2012年1月期间记录的远震波形，运用小波变换频时分析方法分别测定了1216和653条周期从15到140 s的台站间基阶Rayleigh相速度和群速度频散曲线.通过对上述频散进行反演，重构了青藏高原东北缘分辨率高达0.5°×0.5°的2-D相速度和群速度分布图.然后通过对所提取到的每个格网点Rayleigh波相速度和群速度频散进行联合反演，得到了研究区下方一维S波速度结构.最后通过线性插值，得到了青藏高原东北缘下方地壳上地幔三维S波结构.结果表明，印度板块向北俯冲已经达到班公—怒江缝合带附近；在柴达木盆地北部祁连山下面我们发现了亚洲板块，且其没有表现出明显的向南俯冲的迹象；在两大板块中间，我们观测到延伸到250 km深度的低速异常，该低速异常可能是地幔物质底辟上涌现象造成的.中国东北地区噪声层析成像
摘要：中国东北地区是中国唯一的深震孕育区，获取该区的壳幔结构，对于研究板块俯冲、深震以及火山活动等有重要的科学意义.本文利用该区159个固定台站2011年1月至2012年6月和27个流动台站2011年1月至2011年6月间的垂向波形连续记录，计算了台站间的预估格林函数，并采用基于连续小波变换的时频分析方法，测量了双台路径上瑞雷波的群速度和相速度频散曲线.通过质量控制和筛选，最终得到了2204条路径上周期5~40 s的群/相速度频散曲线.检测板测试表明，反演结果可以达到2°×2°的分辨.利用Ditmar ＆ Yanovskaya反演方法，我们得到了研究区（105°E—135°E，39°N—52°N）周期8~30 s的瑞雷波的群速度和相速度分布图.不同周期的群/相速度分布图，反映了不同深度S波速度的横向变化情况.研究结果显示：中国东北地区的地壳上地幔S波速度结构存在横向非均匀性.短周期（如8 s）的群/相速度分布与地表地质构造具有明显的相关性，具体来说，山区显示为高速，沉积盆地显示为低速；随着周期的增大（如15 s，22 s），地形的控制作用相对减弱；较长周期（如30 s）的群/相速度分布与地壳厚度密切相关.基于背景噪声研究云南地区面波速度非均匀性和方位各向异性
摘要：基于中国地震科学探测台阵项目一期在南北地震带南段架设的300多个地震台站，利用2011年10月至2012年9月的连续观测记录，采用基于背景噪声互相关函数的面波层析成像技术，研究了青藏高原东南缘的云南地区面波群速度和方位各向异性分布.结果显示，地壳的面波快波方向呈现近南北向，整体表现出围绕东喜马拉雅构造结顺时针旋转的趋势，和地表GPS速度场以及S波分裂的快波方位较为一致.小江断裂东西两侧的快波方位有一定差异.对反映深度大概在下地壳和上地幔顶部的长周期面波，快波方向从近南北向逐渐向北西向过渡，在菱形块体附近，26°N以南，快波方向和红河断裂的走向趋于一致；其他区域相比上地壳的快波方向也有较大变化，这种结果较为支持青藏高原东南缘的云南地区壳幔变形的解耦.在滇西南，澜沧江向东弧形展布区域，中下地壳快波方位呈现局部的圆周旋转趋势，结合该区地震分布规律和应力主轴方向，推测这种现象和块体挤压及旋转具有一定相关性.云南高黎贡山典型板岩地震波各向异性研究
摘要：板岩作为一种浅变质岩在我国有着广泛的分布，对其地震波速度的研究将有助于对这类过渡性岩石的有效区分，对于浅层地壳的各向异性研究也具有重要意义.本文对采自云南丙中洛地区的板岩样品进行了地震波速度的室内实验研究，其中部分实验是在加拿大Dalhousie High Pressure Laboratory完成.实验获得了板岩在围压10~600 MPa条件下、不同构造主方向（X，Y和Z）上的地震波速度，在围压600 MPa时，X、Y、Z三个方向的P波速度分别为6.58、6.46、5.91 km/s，平均速度为6.30 km/s，S波平均速度约为3.62 km/s，VP/VS=1.74；并初步分析了板岩地震波速度、横波分裂及其波速各向异性随着围压的变化规律，发现所测量的板岩在较低围压（〈150 MPa）时波速的各向异性随围压升高而迅速减小，主要是由于其内部微裂隙的定向排列引起的，而随着围压的继续增加（〉150 MPa时）微裂隙基本闭合，黑云母、阳起石等片状矿物的定向排列成为其地震波各向异性的主导诱因，此时（围压为600 MPa）VP、VS的各向异性分别稳定在13%、16%左右.本研究所获取的基础实验数据及所探讨的板岩地震波性质将为确定地壳上部显微裂隙的优选定向、浅层地壳的各向异性分析、地球物理模型条件约束等提供基础.中国东北—华北地区地壳厚度与泊松比及其地质意义
摘要：本文通过收集和综合分析已有的接收函数H-k研究结果，给出了中国东北—华北地区的地壳厚度与波速比/泊松比分布图.本研究表明该区地壳最薄的地方出现在松辽盆地和华北平原地区（28~35 km）；大兴安岭、燕山—太行地区的地壳厚度介于36~45 km范围，其中燕山造山带地壳厚度由东向西逐渐增加；而最厚的地方则出现在鄂尔多斯盆地西南缘（~55 km）.研究区平均波速比为1.76±0.05，较全球大陆平均值明显偏高，这可能与中、新生代以来该区显著的岩石圈减薄与破坏过程相关.其中地壳波速比最高的地方出现在山西地堑、长白山、大同—张家口等新生代火山区，意味着这些地区可能具有较高的地壳温度或存在广泛的壳内部分熔融.本文研究显示，大兴安岭造山带地区地壳厚度与波速比/泊松比成负消长关系，推测大兴安岭在形成过程中，地壳的增厚以长英质上地壳增厚为主.与大兴安岭地区不同，松辽盆地及周边地区地壳厚度与泊松比没有明显的相关性，表明松辽盆地可能具有复杂的形成与演化过程.声发射矩张量反演
摘要：地震矩张量反演是获取震源过程的有效方法.岩石变形过程中的声发射与地震类似，均是弹性应变能快速释放.如假设条件得到满足，矩张量反演方法同样可用于了解声发射震源过程.声发射矩张量反演可使用P波位移进行计算.当样品尺寸远大于声发射波长，且声发射由微破裂产生，声发射源的尺度很小时，P波矩张量反演可采用远场近似.本文首先针对声发射的特点，实现根据远场P波反演声发射矩张量的算法，并通过人工声发射实验对算法的正确性和可靠性进行了检验.最后，用声发射矩张量反演方法对花岗岩单轴压缩实验的声发射源特征进行了分析.结果表明：对于纯剪切破裂模式，声发射矩张量反演可得断层面；对于非纯剪切破裂模式，如纵向挤压导致的横向张性劈裂，由于存在多解性不能得到断层面，但可通过矩张量的迹区分破裂模式.黏滑实验的震级评估和应力降分析
摘要：本文通过三种结构模型的黏滑地震模拟实验，利用高频速度连续观测系统获得了地震失稳过程的速度特征，讨论了最大位移量的选取方法，估算了实验室黏滑型地震的矩震级，探讨了黏滑类型、应力降大小与震级的关系.结果表明，黏滑型地震的应力降过程可能包含一次到多次高频振荡，对应若干次黏滑子事件.高频振荡的摆动幅度很大，包含有多种频率成分，峰值速度0.003~0.008 m·s-1.初步估计黏滑型地震的震级范围为-4.4~-3级，断层构造面的差异对各种黏滑模型的地震震级分布有明显影响.总体来看应力降与地震震级没有明显相关性，决定地震震级的主要因素应当是震源尺度.伊洛瓦底盆地热-沉降史模拟及构造-热演化特征
摘要：本文首先运用EASY% Ro反演法对伊洛瓦底盆地由北向南进行了热史的恢复，北部钦敦凹陷的平均古地温梯度为13.0~15.0 ℃/km，中部沙林凹陷的平均古地温梯度为18.0~22.0 ℃/km，南部三角洲凹陷的平均古地温梯度为33.0~37.0 ℃/km.从模拟结果可以看出，盆地由北向南地温梯度逐渐升高，生烃门限的深度由深变浅.然后模拟了盆地的构造沉降史.模拟结果表明，盆地具有幕式构造沉降特征，这反映了伊洛瓦底盆地可能处于弧间或弧后的构造背景.伊洛瓦底盆地北部和南部具有不同幕次的构造沉降史，北部在早始新世时期（53~51 Ma）经历了一幕拉伸过程，然后进入了热沉降期，并伴随局部的快速隆升；南部则经历了两幕拉张过程，分别是在早始新世时期（53~51 Ma）和中中新世时期（21~13 Ma）.盆地的这种南北构造沉降的差异很可能是造成盆地地温梯度北低南高的原因.间断有限元方法的数值频散分析及其波场模拟
摘要：数值求解波动方程是大尺度正演波场模拟、基于波动方程的地震偏移和反演成像的关键.本文针对求解二维声波方程的Runge-Kutta 间断有限元（RKDG）方法的数值频散问题，从理论推导和数值分析的角度进行了深入研究，并将其与近似解析离散化方法（Optimal Nearly Analytic Discrete Method，简称ONAD 方法）、Lax-Wendroff 修正方法、交错网格（Staggered-Grid，简称SG）方法的数值频散进行了比较研究.结果表明：RKDG方法以及近似解析离散化方法在压制数值频散方面要好于上述其他方法，特别是空间精度为3阶的RKDG方法，即使当空间步长取波长的一半，即一个波长内取2个网格点时，最大的频散误差也不超过1.67%.同时，我们也通过波场模拟对比研究了不同数值方法的数值频散问题，进一步直观地验证了数值频散的理论分析结果.多震源地震数据偏移成像方法
摘要：多震源地震技术是一种高效的地震数据采集方法技术，得到的地震记录是来自多个震源的混合地震数据.本文在多震源波场传播理论和地震波场满足线性叠加原理的基础上，提出了两种多震源地震数据的偏移成像方法.第一种方法是首先对多震源地震数据进行分离，得到各个单震源的地震数据，然后再利用常规的偏移成像方法进行处理；第二种方法是多震源地震数据的直接偏移成像.把本文提出的多震源偏移成像方法应用于数值模拟的多震源地震数据，验证了本文方法的正确性和有效性，直接偏移成像方法较分离后再偏移方法具有更高的计算效率.消除探地雷达数据的子波衰减和频散的反滤波方法
摘要：消除探地雷达数据的子波衰减和频散可以很好地提高探地雷达的勘探深度和勘探分辨率.常用的消除探地雷达数据的子波衰减和频散方法为反Q滤波方法.该方法需要利用地下介质的Q参数，但是正确求取地下介质的Q参数很困难.针对这一问题，本文提出了一种消除探地雷达数据的子波衰减和频散的反滤波方法.该方法以地下介质反射系数是随机数为前提，利用地下介质等效滤波器具有最小相位这个特性，通过求取等效滤波器的振幅谱来求取等效滤波器的反滤波器.最后，利用该反滤波器对探地雷达数据进行反滤波，实现消除探地雷达数据的子波衰减和频散.基于有限差分正演的带地形三维大地电磁反演方法
摘要：本研究实现了一套基于有限差分（FD）方法的大地电磁测深数据带地形三维反演算法及代码.其中，在大地电磁场正演数值模拟方面，开发了起伏地形条件下基于交错网格剖分、有限差分方法的大地电磁测深三维正演代码；在满足平面波场假设的前提下，使用长方体网格剖分模拟三维起伏地形，实现了带地形三维正演计算；并设计理论模型进行试算，经试算结果与前人的有限元法计算结果对比，验证了所研发的带地形三维正演计算的正确性与可靠性.在反演方面，本研究基于非线性共轭梯度方法编写了大地电磁测深带地形三维反演代码，试验了不同的共轭梯度搜索因子β，避免了目标函数对海森矩阵（参数二次导数矩阵）的显式计算和存储，初步实现了大地电磁资料的带地形三维反演.最后，对一系列理论模型进行正演计算，利用其生成的合成数据模拟实测数据进行反演，并与现有的不带地形大地电磁测深三维反演结果比较，检验了所研发的带地形三维反演计算的可靠性与稳定性.频率域航空电磁数据的加权横向约束反演
摘要：传统的一维反演技术已经被广泛应用于航空电磁数据解释中.然而，利用单点水平层状介质模型模拟地下复杂地电结构有时会遇到困难.突出表现在反演参数的横向不连续性，即使相邻测点的反演结果也会出现突变.本文针对航空电磁直升机吊舱系统可进行密集采样，相邻测点地下电性结构应具有某种程度连续性的特点，研究航空电磁数据横向约束反演理论，并提出参数加权约束方法.首先阐述频率域航空电磁正演和加权横向约束反演理论，着重介绍这种拟二维反演方法的基本原理和实施步骤，以及将该方法成功应用于频率域航空电磁数据反演处理的方法技术.最后，通过对理论和实测数据反演处理，并与传统的一维反演结果进行对比，验证加权横向约束反演方法的有效性.复杂地电结构条件下统计学建模法电阻率测深二维反演
摘要：为了提高复杂地层结构条件下，以统计学建模方法为基础的电阻率测深二维反演的准确性，本文改进了统计学建模的技术细节，建立统计学特征值与地质体空间位置、地层层参数的连接关系，并从建模环节开始加入先验信息.以常见的断裂模型为例，研究了用统计学判断地层结构建立初始模型、用加权统计的水平梯度确定断裂位置及用分段统计确定局部层参数，并在初始模型中加入先验信息.理论模型试算和对已知断裂上的实测高密度数据的反演表明，在复杂地层结构条件下，统计学建模方法可以建立较贴近实际的初始模型，通过反演获得清晰明确的结果.基于各向异性岩石物理的缝隙流体因子AVAZ反演
摘要：裂缝型储层表现出较强的各向异性特征.缝隙中充填不同流体时，裂缝储层的地震响应特征也不相同.本文从各向异性岩石物理模型出发，引入可有效识别缝隙流体的指示因子，并研究缝隙充填流体类型、饱和度以及缝隙纵横比与流体因子的相互关系，进而分析不同流体充填时介质的地震响应特征，并基于AVAZ反演方法估测缝隙流体指示因子.首先对缝隙流体因子的敏感性进行了分析，讨论当缝隙充填不同流体时，缝隙流体因子值的变化特征，同时研究了不同流体类型充填时裂缝储层反射系数随方位角和入射角的变化特征.某工区测井数据和复杂裂缝模型应用表明，基于各向异性岩石物理的缝隙流体因子AVAZ反演方法合理、可靠，且具有良好的抗噪性，即当对合成地震记录添加信噪比不小于1/2的随机噪声时，利用AVAZ反演方法估测所得流体因子值与真实值仍然吻合较好.
