本文利用MAVEN卫星Langmuir Probe and Waves(LPW)仪器的在轨电子浓度探测数据,研究了火星电离层电子浓度随太阳天顶角(Solar Zenith Angle,SZA)的变化以及昼夜电子浓度变化的异同.基于2014年至2017年期间MAVEN的电子浓度数据,我们发现:在200km以下,白天电离层电子浓度主要受光化学平衡控制,由于白天光电离过程使得昼夜电子浓度差异较大,此...
我们采用COSMIC掩星系统提供的2009年1月1日至2010年12月31日的GPS信号信噪比(signal-tonoise ratio,SNR)的观测数据,通过计算SNR数据的归一化标准差对80~130km高度范围内的中国内陆偶发E层(Sporadic E,Es)的空间分布以及垂直厚度等特性进行分析研究,获得了Es发生率的空间分布以及发生率随着季节的变化特征.中国内陆30°N纬线以上区域Es分布...
为了提高空间磁场探测精度,磁测卫星一般将磁强计安装在长伸杆的顶端,同时通过整星磁洁净控制,来减小磁强计安装位置的磁场干扰量.针对低轨磁测卫星的特点,提出零磁场环境下卫星极弱干扰磁场的高精度测量方法,模拟在轨地磁场环境下卫星感应磁场的标定方法,以及磁力矩器剩磁不确定量的处理方法.通过以上方法对卫星本体产生的干扰磁场实施严格的标...
利用神经网络算法挖掘海量数据的规律已成为科技发展的一种趋势,本文针对卫星信号的天顶对流层延迟进行建模.对流层延迟是影响卫星定位精度的重要因素之一,建立精密区域对流层模型对高精度定位有着重要的意义.对区域测站对流层延迟数据的分析,考虑到实时建模中传统BP(Back Propagation)神经网络计算量大,易出现"过拟合"现象、不稳定等因素,...
地面摩擦和大尺度流场是影响锋面气旋结构演变的重要因子,本文使用WRF模式并采用湿物理方案,通过理想化试验,综合考虑陆面摩擦、气旋式扰动相对于急流位置和大尺度流场对锋面气旋结构变化的影响.结果表明:当仅考虑单一因子时,气旋式扰动位于急流南侧和辐合流流场有利于气旋形成Shapiro-Keyser(S-K)模型结构.当同时考虑地面摩擦和大尺度辐合时...
通过分析海空重力测量系统误差的形成机理,我们发现海空重力仪格值标定误差是引起系统性测量偏差的主要因素之一.本文简要介绍了重力仪格值的标定方法,分析论证了格值标定的精度要求,提出了利用东西正反向重复测线检测校正海空重力仪格值的计算模型和补偿方法,分析讨论了该方法的校正精度及其适用条件,利用航空重力实际观测网数据对该方法的合理...
大地震发生之后通常会诱发一系列的余震序列,对比1976年MS7.8唐山大地震和2001年MS8.1昆仑山大地震周边区域的地震事件可以看出,唐山大地震余震活动时间要明显长于昆仑山大地震余震活动时间.余震序列往往与震后形变密切相关,而影响震后形变的因素不仅与地震发震断层和震级有关,同时与岩石圈的结构有关.考虑到唐山大地震的发震区华北地块和昆仑山...
2017年11月18日,我国西藏林芝市米林县发生M6.9地震.其东南220km距离处,1950年发生了察隅8.6级大地震.二者同处喜马拉雅东构造结大拐弯处,周边分布多条活动断裂带.察隅地震作为20世纪我国内陆最大的地震,它的发生对周边断层活动性的影响目前研究尚少.本文基于分层半无限空间黏弹性地球模型,计算了察隅地震对周围活动断裂带和米林地震的影响,同时...
本文基于结合速率-状态摩擦定律(RSF)的二维准动力学数值模型,以半空间垂直走滑断层为研究对象,通过比较两种正应力随深度变化模型的模拟结果,研究了浅层正应力变化对断层演化参数、地震孕育过程、震后滑移传播等方面的影响.结果显示,我们的数值模型在给定模型参数和约束条件下,能够完整模拟出地震周期中震间、震前、同震以及震后多个特征阶段...
2012年发生的苏门答腊大地震激起了强烈的环形地球自由振荡,本文利用南极26个地震台站记录的此次大地震激发的0T10简正模,分析南极上地幔结构异常对长周期自由振荡的影响,我们将震后20h的自由振荡观测数据与利用PREM地球模型模拟的结果进行比较,结果表明:在震后较短时间内区域性(区域的范围远远小于简正模的波长)上地幔结构异常能显著影响长...
云南地区地处青藏高原东南缘,一直是地球科学研究的热点地区.目前,一些热点问题,如云南地区是否存在中下地壳低速流及其空间分布,仍有一定的争议.通过建立云南地区精细的地壳上地幔速度与各向异性结构,可为深入认识上述问题提供重要信息.本文利用天然地震波形记录,采用双台法提取了12~60s周期的Rayleigh和Love波相速度频散,并进一步反演获得了云...
利用中国数字测震台网(CDSN)记录到的台湾地区两个地震事件6°~30°震中距范围的三重震相波形资料,基于观测与理论波形拟合法,获到华南地区上地幔P波和S波的最佳波形拟合速度模型及其V_P/V_S比值.与AK135模型相比,华南地区410km深度上方存在明显低速层:S波低速区厚度约为70km,速度降为2%~5%;而P波低速区厚度为70~230km,速度降为5%~6%.另外,410k...
针对张渤地震带深部孕震构造环境和地幔岩浆对地壳底侵作用的问题,本文充分收集了华北地区区域数字地震台网176个固定台站记录到的观测数据,应用近震体波层析成像方法,获得了华北地区地壳的P波三维速度结构和VP/VS波速比.依据VP/VS波速比能反映岩石物理性质和流变学特征的特性,并结合人工地震测深剖面、大地电磁测深、地球化学等成果,讨论了张渤...
超慢速扩张的Mohns洋中脊共轭两侧的地球物理场与地壳结构具有显著的非对称性.利用我国第五次北极科学考察采集的水深、重力与磁力数据,结合历史资料,我们计算了14条垂直Mohns洋中脊剖面的扩张速率、剩余水深、剩余地幔布格重力异常(RMBA)、地壳厚度和非均衡地形.对洋中脊共轭两侧以上计算结果的进一步对比发现,Mohns洋中脊两侧整体(下文均指...
位场的解析延拓是实现不同高度海洋地磁场相互转换的主要途径,是构建海洋三维磁空间背景场模型的关键技术.针对位场向下延拓迭代法中最优正则化参数及最佳迭代次数难以确定问题,尝试引入微分进化法,以正则化参数及迭代次数为种群变量,以延拓结果的熵值为目标函数,以目标函数最小化为搜索准则,实现两种参数的并行全局寻优.采用实测数据对微分进化...
由于重力引起的岩石压实效应,一般来说,地震波传播速度由浅入深整体逐渐增大.梯形坐标系设计可耦合速度由浅入深逐渐增大的变化,该坐标系中均匀网格采样所对应的物理直角坐标系网格由浅入深逐渐增大,也即浅部低速区对应细网格,深部高速区对应粗网格.在梯形坐标系表征波动方程后利用有限差分求解,本文实现一种深度均匀采样、横向采样间隔随深度增...
对稀疏/非规则采样或者低信噪比数据,射线束提取困难并伴随有假频产生,对叠加剖面和道集造成严重干扰.为了提升射线束偏移在稀疏和低信噪比地震数据采集中的成像效果,本文提出基于三角滤波的局部倾斜叠加波束形成偏移假频压制方法.射线束偏移首先将地震数据划分为超道集,经过部分NMO后转化为以射线束中心定义的共偏移距数据,倾斜叠加和反假频操...
针对裂缝各向异性介质,本文提出一种非正交假设下的矢量波场分离方法.本文首先对多分量地震勘探中常见的波型泄漏现象进行了数学描述,提出在纵、横波波场分离的同时应该考虑恢复纵、横波的矢量振幅.为了对裂缝方位角与各向异性系数进行定量预测,本文将矢量波场分离拆分成三个步骤来实施:第一步,用Z、R两分量的仿射坐标系变换分离ZR平面内的P波...
时间域的波场延拓方法在本质上都可以归结为对一个空间-波数域算子的近似.本文基于一阶波数-空间混合域象征,提出一种新的方法求解解耦的二阶位移弹性波方程.该方法采用交错网格,连续使用两次一阶前向和后向拟微分算子,推导得到了解耦的二阶位移弹性波方程的波场延拓算子.由于该混合域象征在伪谱算子的基础上增加了一个依赖于速度模型的补偿项,...
叠前逆时偏移等基于波场互相关原理的地球物理方法存在极大的计算与存储需求,因此采用合适的波场重构方法显得尤为重要.常规的随机边界法容易产生成像噪声,而有效边界法在三维情况仍难以实现,检查点技术具有内存要求小的特点,但存在较高的重算率,因此本文提出了插值原理的检查点技术波场重构方法.在满足Nyquist采样定理的前提下对相邻检查点间的...
巨大的计算量是制约全波形反演(FWI)生产实用化的难题之一.为此,本文提出了一种高效的波场迭代解法,将其应用于频率域常密度声波方程FWI,并给出了详细的反演流程.通过建立用于波场迭代的目标函数,推导相应梯度、步长公式,新方法将反演中波场正传和残差波场反传过程转化为无约束优化问题,从理论上分析了新方法的计算效率显著高于常规FWI.在数值...
本文以基于改进BISQ模型的二维双相各向同性介质一阶速度-应力方程为基础,推导出了曲线坐标系下对应的方程,然后采用低频散、低耗散的同位网格MacCormack有限差分法来离散方程,并采用紧致的单边MacCormack差分格式结合牵引力镜像法来施加自由地表边界条件,实现了地震波场数值模拟.曲线网格有限差分法采用贴体网格来描述自由表面,地表的网格线紧...
槽波地震勘探利用槽波的频散特性反演煤层的结构特征,故理论频散曲线的计算是一个重要方面.使用水平层状模型假设下的面波频散曲线计算方法能够计算煤层厚度恒定模型地震槽波频散曲线;但当煤层厚度变化时该方法不再适用.基于前人水平层状均匀介质模型的面波理论频散曲线计算方法,对于含煤三层模型,本文发展了煤层厚度变化情况下的地震槽波理论频...
目前,对于大庆地区的地质构造研究成果仅局限在深部构造上,该地区从未开展过针对近地表隐伏断裂的探查工作.本文采用浅层反射地震勘探方法,查明了克山—大安断裂嫩江组以上地层的详细地层信息以及断裂的展布形态;同时,在主干断裂上覆的背斜构造中,发现了许多次级断裂,这些次级断裂在前人的成果中并未提出过,并且在本区的断裂-褶皱构造体系中,次...
基于不同测量高度重力场及其梯度数据可同时对应同一场源并用于反演场源位置的分析原理,拓展不同高度场数据在欧拉反褶积法中的应用范围.首先,立足于对欧拉反褶积方法的理论研究基础,提出不同高度数据融合联合欧拉反演公式.其次,在理论模型上对多种高度数据联合反演做了测试分析计算,验证了不同高度场数据融合联合欧拉反褶积法能够改善位场解释...
伊犁盆地中部区域地处多个构造单元的交汇区,构造复杂,缺乏高精度地质资料,制约了对盆地构造演化的进一步认识.本文联合使用多种先验信息作为约束条件,利用高精度重磁数据反演获得了研究区地下海拔-10km以浅的三维密度和磁性结构,增强了反演结果的可靠性.结果揭示以白石墩次凸为中心发育了一个"北断南超"的南西-北东向不对称型凹陷,北部沉积厚...
针对常规大地电磁(Magnetotelluric,MT)反演方法对电阻率异常体边界不太敏感的问题,本文尝试基于贝叶斯理论开展二维大地电磁电阻率尖锐边界反演研究.在反演中,模型参数由边界位置及内部电阻率组成,通过贝叶斯理论将模型参数与数据相联系,采用Markov Chain Monte Carlo(MCMC)的Metropolis-Hastings(MH)方法对后验概率密度函数(Posteriori...
针对基于互相关成像条件的探地雷达(GPR)逆时偏移计算效率低、存储量大及易产生低频假象的不足,本文将激发振幅成像条件应用于GPR逆时偏移成像中.通过在源点电磁波场正向传播过程计算每个网格点的能量密度,并保存最大能量密度的时刻和相应的电磁波场值;在接收点电磁波场逆向传播过程提取每个网格点最大能量密度时刻及对应的电磁波场值,并利用...
在多孔隙含水层中地面磁共振(surface nuclear magnetic resonance,SNMR)信号呈现多弛豫衰减特性,常规盲源分离方法和单指数拟合方法引起信号严重失真和信息缺失等问题.本文提出了基于稀疏表示的随机噪声背景下多弛豫SNMR信号的提取方法.根据SNMR信号的衰减特征,设计了精确刻画SNMR信号且与随机噪声不相关的离散衰减余弦冗余字典.其次,针对多...
为研究井周裂缝发育特征,本文提出一种新型方位侧向测井方法,利用三维有限元法,模拟裂缝的方位侧向测井响应.结果显示,深浅侧向电阻率幅度差异受裂缝倾角的控制,低角度缝为负差异,高角度缝为正差异;倾斜裂缝张开度的增大使测井响应值减小,方位电阻率差异增大;井周方位电阻率可反映裂缝方位产状,单一缝或裂缝密度较小时,沿裂缝走向的方位电阻率小...