利用差分体像素CT由TWINS卫星ENA成像反演暴时环电流离子三维分布
摘要：本文将体像素CT方法应用于TWINS卫星高能中性原子立体成像数据,重建环电流离子通量强度三维分布.为了克服卫星ENA成像可能存在的仪器偏差造成反演误差增大甚至迭散无解的问题,发展了差分ENA体像素CT反演方法.采用环电流离子分布理论模型和实际卫星观测几何构形进行模拟检验,证明了该方法的可靠性和优越性.利用此方法分析2012年7月15日磁暴（min Dst=-133nT）主相期间TWINS两颗卫星的实测ENA数据（能量范围4keV～50keV）,反演得到环电流离子微分通量随L值/纬度/地方时的三维分布与能谱.所得离子通量呈现显著不对称的部分环电流特征,主要分布在低纬-赤道区磁午夜前后至黎明前数小时,L值约为3.5～6.5的区域;将反演得到的磁午夜后赤道区环电流离子通量能谱,与THEMIS卫星当地测量得到的此区域同时的能谱作比较,两者符合得很好;证明本文发展的差分体像素CT是由多卫星ENA二维图像重建暴时环电流离子分布的有效方法.基于CARISMA地磁台链观测的IPDP事件的统计分析
摘要：亚暴期间磁尾等离子体片离子注入内磁层能够激发电磁离子回旋（EMIC）波.对应于这种EMIC波,地面磁力仪可观测到周期逐渐减小的地磁脉动（IPDP）.利用GOES卫星数据,地磁指数和加拿大CARISMA地磁台站的数据,我们研究了IPDP事件的产生与亚暴磁尾注入的关系.同时利用CARISMA地磁台链中的MCMU和MSTK两个台站,从2005年4月到2014年5月期间的观测数据,统计分析了亚暴期间的IPDP事件,研究了IPDP事件的出现率关于季节和磁地方时的分布特征.我们总共获得128个两个台站同时观测的IPDP事件.该类事件关于季节分布的发生率,冬季最小,为13.28%,春季最大,为32.81%,结果表明IPDP事件关于季节分布的发生率受到电离层电导率及亚暴发生率的影响.两个台站同时观测到的IPDP事件最大出现率出现在15—18 MLT（磁地方时）,结果表明IPDP事件主要由亚暴期间产生的能量离子注入内磁层,西向漂移遇到等离子体层羽状结构（Plume）区的高密度等离子体所激发.宇宙线小波分析在大地磁暴预报中的应用
摘要：利用Morlet小波变换方法对北京宇宙线台站的地面宇宙线强度在地磁暴前后的变化特征进行分析,得到：1）在平静期,北京宇宙线数据存在准24h周期性的特征,且通过分析周期为12h的Morlet小波＂模＂,发现值稳定,且小于0.6;2）以90天为时间窗口,对2004年7月地磁暴前后的小波频谱变化进行详细分析,发现当发生大地磁暴时,宇宙线的静日准24h周期被打破,其他周期的波动开始增强.进一步研究发现,周期12h的波动在大地磁暴数小时到1天左右会出现显著增强,这一现象在2001、2002和2004年期间的大地磁暴前得到验证.3）Morlet小波＂模＂数据的急速增大是发生地磁暴的先兆特征,当小波模变化达到一定的阈值就可能发生大磁暴.本文分析了周期为12h时小波的模数据,对强地磁暴事件进行统计,选定阈值0.6,并通过2003年的6次大地磁暴进行预报验证,结果表明该方法不仅能够对大地磁暴事件进行预报,而且提前量满足预报需求,为基于宇宙线实测数据预报地磁暴方法提供了重要基础.中低纬地区电离层对CIR和CME响应的统计分析
摘要：本文利用中低纬日本地区（131°E,35°N）GPS-TEC格点化数据,分析了2001—2009年间109个共转相互作用区（CIR）事件、45个日冕物质抛射（CME）事件引起的地磁扰动期间电离层的响应.结果表明,电离层暴的类型随太阳活动的变化而有不同的变化,CIR事件引发的电离层正相暴、正负双相暴多发生在太阳活动下降年,负相暴多发生在高年,负正双相暴多发生在低年;CME事件引发的电离层正相暴和负相暴多发生在高年.CIR和CME引发的不同类型的电离层暴的季节性差异不大,在夏季多发生正负双相暴.电离层暴发生时间相对地磁暴的时延大部分在-6～6h之间,但CIR引发的电离层暴时延范围更广,在-12～24h之间,而CME引发的电离层暴时延主要在-6～6h之间.中低纬的电离层暴多发生在主相阶段,其中CIR引发的双相暴也会发生在初相阶段.电离层负暴多发生在AE最大值为800～1200nT之间.CIR引起的电离层扰动持续时间较长,一般在1～6天左右,而CME引起的电离层扰动持续时间一般在1～4天左右.NWC通信台在电离层中激发电磁响应的时变特征
摘要：本文利用DEMETER卫星VLF频段电场和磁场频谱数据对DEMETER卫星运行期间2005年至2009年澳大利亚甚低频（Very Low Frequency）通信台NWC发射的通信信号造成的电离层电磁响应的日变化、季节变化及年变化特征进行了统计分析,统计结果表明电磁响应日变化显著,夜间电场强度明显增强可达40dB,磁场变化略小也可为15dB左右,而季节变化不显著,年变化主要受太阳活动的影响,太阳活动越强,电磁响应越小.为解释数据分析结果,对地-电离层电磁波传播过程采用传递矩阵方法进行了模拟计算,模拟结果与数据分析的结果一致.我们认为这种随时间变化的特点可能由250km以下电离层电子密度分布特征导致,因此研究250km以下的电离层电子密度变化可能对寻找地震电离层电磁异常有重要意义.基于纯转动Raman激光雷达的中低空大气温度高精度探测
摘要：本文介绍了一套纯转动Raman测温激光雷达系统,通过高分辨光谱分光与滤光优化设计、收发精确匹配以及弱信号检测等技术,实现在武汉城市上空从10km至40km的中低空大气温度高精度探测.观测结果与同时段探空气球进行比对,在30km以下激光雷达探测温度与探空气球得到的温度数据吻合较好,最大偏差约为3.0K,表明了该激光雷达温度测量的可靠性.采用30min时间分辨率,在10～20km高度范围内温度统计误差约为0.3K（300m空间分辨）;20～30km统计误差约为0.8K（600m空间分辨）;30～40km统计误差约为3.0K（900m空间分辨）.通过整晚的温度廓线反演,为研究中低层大气中的波动现象提供依据.该转动Raman激光雷达实现了至40km高度的高精度大气温度探测,进一步可与Rayleigh测温激光雷达30～80km的高度衔接,为实现中低层大气连续观测研究提供了重要手段.一种能实现单频PPP-RTK的GNSS局域参考网数据处理算法
摘要：全球范围内大量布设的GNSS（Global Navigation Satellite System）参考网为精密定位、导航和授时等应用提供了丰富的数据资源.基于局域参考网,先后发展了若干侧重实现双频精密定位的技术,如NRTK（Network Real Time Kinematic）,PPP（Precise Point Positioning）和PPP-RTK等.其中,PPP-RTK融合了NRTK和PPP的技术优势,是目前相关研究的热点.本文改进了利用局域参考网提取各类改正信息的算法,以便于实现单频PPP-RTK,具体步骤包括：1）逐参考站实施非组合PPP,并固定已知站星距和卫星钟差,预估电离层延迟、浮点模糊度等参数;2）联合所有参考站的PPP模糊度预估值,通过重新参数化,形成一组双差整周模糊度和接收机、卫星相位偏差;3）固定双差整周模糊度,精化求解卫星相位偏差和各参考站PPP电离层延迟.基于网解中用到的卫星轨道和钟差,以及网解所提供的卫星相位偏差和（内插的）电离层延迟,参考网内的单频流动站即可实施PPP-RTK.基于澳大利亚某连续运行参考站网和流动站的实测数据,考察了：1）参考网数据处理中,双差模糊度的固定成功率（98.89%）和卫星相位偏差估值的时间稳定性（各连续弧段优于0.2周）;2）流动站处电离层延迟的内插精度（优于10cm）;3）单天内任一历元起算,固定静态（动态）单频PPP整周模糊度所需时长（均不超过10min）;4）模糊度固定前后,单频动态PPP的定位精度（模糊度固定后,平面和天顶RMS分别优于5cm和10cm;模糊度固定前,相应RMS仅为28～53cm）.一种InSAR大气相位建模与估计方法
摘要：为了削弱大气延迟对干涉结果的影响以提高InSAR的测量能力,本文在InSAR大气相位特征分析的基础上,研究了一种新的InSAR大气相位建模与估计方法.首先采用稳健估计确定大气垂直分层部分的模型参数,然后利用基于Matern模型的Kriging插值估计大气紊流部分,最后应用估计的大气垂直分层和紊流资料改正InSAR测量结果.利用覆盖河南义马地区的ASAR数据对本文提出的方法进行了验证,结果表明去除大气影响后,InSAR重建的DEM与参考DEM的高程差异的均方误差由19.5m降至5.3m,精度提高了约72%.同时,改正后的干涉图更合理地揭示了义马矿区的沉降漏斗情况,进一步验证了本文方法的有效性.基于GNSS时空数据的瞬态无震蠕滑信息检测
摘要：利用分布活动断层带的GNSS位移时空序列,基于卡尔曼滤波和主成分分析,研究了一套集GNSS时空噪声处理与断层瞬态无震蠕滑信息时空分布检测于一体的方法.基于一阶高斯马尔科夫理论表述了瞬时蠕滑信息与时空噪声的运动状态方程,通过卡尔曼滤波提高了时间域信噪比;根据断层形变高空间相关性的特点,通过主成分时空分析,提高了空间域信噪比.模拟试验结果表明,当断层滑移引起的地表位移大小至少与噪声水平相当时,均可检测得到断层瞬态无震蠕滑时空分布及其滑移特征.以Cascadia慢滑移事件为例,清晰地检测到了两次无震蠕滑信息,分析其蠕滑特征与相关研究结果吻合.通过对滇西2011—2013年陆态网GNSS连续站数据的处理,分析发现了期间有微弱的震后余滑信息,主要表现在澜沧江断裂、红河断裂和小江断裂的南段.其时空分布特征与2011年3月24日的缅甸MW7.2级地震相对应.因此,得出了2011—2013年期间云南地区的断层活动可能与缅甸地震带有着密切的关系的结论.2014年美国加州纳帕MW6.1地震断层参数的Sentinel-1A InSAR反演
摘要：2014年8月24日,在美国加州旧金山海湾北部的纳帕地区发生了MW6.1地震.发震断层是西纳帕断裂系统中的一部分,但是该断层之前并未被足够重视.本文利用欧洲空间局最近发射成功并刚刚投入使用的Sentinel-1A卫星获取的第一对同震干涉像对（20140807-20140831）,得到了该地震的地表同震形变场,结合震后24h内区域GPS同震形变资料作为约束条件,反演了纳帕地震的断层几何参数以及滑动分布.Sentinel-1A干涉结果表明,此次地震造成了明显的地面形变,视线向最大抬升和最大沉降量均达到了10cm.联合反演结果表明,该发震断层的走向为344°,倾角为80°.主要破裂以右旋走滑为主,平均倾滑角为-146.5°,最大倾滑量达到了1.1m,位于地表下约4km,存在明显的滑动亏损现象.此次地震,累计释放地震矩达1.5×10^18 N·m,约合矩震级MW6.1.该结果略小于InSAR单独约束结果,可能与Sentinel-1A像对中包含的快速震后形变分量有关.台湾碰撞带现今地壳变形场特征及其动力学成因的有限单元法模拟研究
摘要：台湾是菲律宾海板块与欧亚板块会聚、碰撞的产物,区内地质构造复杂,地震频发,变形强烈.为定量研究台湾地区变形特征及其动力学成因,文中运用二维不连续变形体弹性力学的有限单元计算方法,利用台湾地区1995—2005年GPS的观测结果作为边界约束,对台湾地区地壳变形场进行了模拟.计算结果显示,台湾中部地区主要呈压缩状态,但在台湾的东北部及南部地区出现了拉张的变形环境.变形程度最大的区域位于台湾东部海岸山脉及其附近海域,同时,计算给出台湾纵谷断裂滑移速率约为13.8～23.5mm/yr,由于纵谷断裂对变形的吸收,因此变形在纵谷断裂以西及西北地区迅速衰减.此外,计算结果还发现,计算给出的台湾岛上的速度值与GPS观测结果吻合得较好;计算给出的主应力方位与地应力观测及震源机制解结果也颇为一致.由此说明文中的有限元模型是合理的.此外,计算结果还表明,菲律宾海板块向北西方向的推挤,板块边界形状,观音、北康高地的抵阻,冲绳海槽扩张及琉球海沟的向南后撤以及断裂作用等,共同造成了台湾地区现今变形场的主要格局.查干凹陷中、新生代构造-热演化史
摘要：查干凹陷是中国内蒙古银根-额济纳旗盆地中最具有勘探潜力的凹陷.为了揭示其构造-热演化历史,本文利用35个磷灰石裂变径迹和119个镜质体反射率数据,采用耦合反演的方法恢复了查干凹陷白垩纪以来的热历史.结果显示查干凹陷白垩系具有高的古地温梯度,并且查干凹陷经历了地温梯度快速增加阶段（K1b-K1s）,地温梯度高峰阶段（K1y）,高地温延续阶段（K2w）和热沉降阶段（Cz）四个构造-热演化阶段.此外,基于热史恢复结果,正演获得查干凹陷高的古地温梯度有利于烃源岩成熟生烃,早白垩世的高古地温梯度控制着该地区的烃源岩有机质的热演化.本文的研究成果可以为下一步油气资源评价和勘探决策提供基础依据.川中古隆起超压分布与形成的地温场因素
摘要：温度和压力是沉积盆地两个重要的物理场,温度影响着超压的形成和分布.本文根据钻孔实测温度和压力数据分析了川中古隆起现今压力与温度的关系;在实验室对封闭流体进行了多组温-压关系实验;利用等效镜质体反射率和包裹体测温数据恢复了川中古隆起不同井区在白垩纪抬升之前的最大古地温,并在此基础上分析了温度降低对研究区超压的影响;最后探讨了生烃增压和欠压实超压形成过程中温度的作用.研究结果表明,川中古隆起现今超压层的压力系数与温度呈正相关关系;在绝对密封的条件下,当压力大于15 MPa时,温度每变化1℃,压力变化1.076 MPa.川中地区不同井区自晚白垩世以来的差异性降温是现今同一超压层系超压强度不同的主要因素,此外超压层还应发生了流体的横向压力传递和泄漏.下古生界原油裂解形成超压的时间是180～110 Ma;气态烃伴生的盐水包裹体均一温度暗示了在90 Ma超压发生调整.盆地模拟结果显示温度对上三叠统须家河组的欠压实增压影响微弱.两淮煤田大地热流分布及其构造控制
摘要：基于127块煤系地层岩石样品的热导率测定结果,并结合59个井田内可靠的系统测温数据,计算得出两淮煤田的大地热流值,并编制大地热流分布图,其结果表明：两淮煤田大地热流值变化范围为29.7～83.9mW·m^-2,平均值为58.3mW·m^-2,和其他沉积盆地存在一定的差异,且淮南煤田大地热流值（63.7mW·m^-2）远大于淮北煤田（55.2mW·m^-2）.综合分析得出,两淮煤田大地热流与其他盆地的差异以及淮南煤田热流值高于淮北煤田的现象为构造演化和区域地质背景的控制结果;而研究区内热流的分布不均主要是由于受地质构造对地温场的影响所致,推覆构造上下盘现今热流值的差异尤为突出.攀西地区重力场特征及地壳密度结构
摘要：攀西地区位于峨眉山大火成岩省中西部,构造和岩浆特征显著,地震活动强烈.通过对野外重力测量得到的云县—会东和普洱—七甸两条剖面的高精度重力观测数据进行处理和分析,构建了沿剖面的二维地壳密度结构,其中普洱—七甸剖面与孟连—马龙宽角地震剖面部分位置重合.同时结合区域重力异常特征及下地壳视密度填图结果,得到如下初步认识：红河断裂带是南北地震带南段地区重要的构造分界线,断裂带南北向密度结构和莫霍面分布形态存在较大差异,沿走向构造变化.云县—会东剖面上大姚—会东段下地壳底部存在密度较高的壳幔过渡层,结合研究区下地壳底部壳幔过渡层的密度分布特征,认为该过渡层不是攀西裂谷下的＂裂谷垫＂,而是由岩浆底侵作用造成的.四川盆地东西陆块中下地壳结构存在差异
摘要：四川盆地是中上扬子克拉通的主要组成部分.作为我国三大稳定克拉通之一,扬子克拉通经历了自太古代以来的长期演化,直到新元古代晚期与华夏板块发生碰撞拼合前,一直被认为是一个稳定的统一陆块.基底包括了新太古宙-新元古代岩层,其上广泛被新元古代晚期至显生宙地层覆盖,仅有～2.9—2.95Ga基底岩石零星出露于四川盆地的西缘、西南缘和三峡地区,使得对于沉积盖层之下的中下地壳的性质和分布规模的认识十分有限.重力异常则能够宏观揭示区域结构特征.本文通过刨除沉积盖层和莫霍面起伏引起的重力异常而获得了中下地壳的重力异常,反映了四川盆地东西陆块中下地壳存在结构差异,结合深地震反射资料、航磁异常和地球化学资料,证实了该分界线位于重庆—华蓥一线,故而推测中上扬子克拉通在太古宙-古元古代可能存在东西两个陆核.兰州盆地渐新世沉积物岩石磁学性质探究
摘要：我国西北地区出露大量富含哺乳动物化石的＂第三纪红层＂,探明这些地层的岩石磁学特征对进一步开展磁性地层、环境磁学和古气候学研究具有重要意义.本文对青藏高原东北缘兰州盆地渐新世地层进行了详细的岩石磁学研究,分别确定了咸水河组下部砂岩和泥岩中磁性矿物含量、种类及其在加热过程中的转变;并结合X射线衍射,漫反射光谱和粒度分析等非磁学手段,揭示了磁性矿物种类与沉积环境的关系.结果表明,剖面底部和顶部的砂岩为河道沉积,磁性矿物含量较低,主要为高矫顽力的赤铁矿和针铁矿;中部泥岩为湖相沉积,指示了一个相对稳定的静水沉积环境,其磁性矿物含量比砂岩要高,不仅含有高矫顽力的赤铁矿和针铁矿,也含有低矫顽力的磁铁矿.在岩石磁学和环境磁学结果的基础上,结合野外考察和古生物证据,我们推测兰州盆地在晚渐新世为半湿润半干旱的气候条件,这为赤铁矿的生成提供了有利条件,导致盆地中富含赤铁矿的红色泥岩广泛发育.西秦岭造山带（中段）及其两侧地块深部电性结构特征
摘要：本文对跨过西秦岭造山带（中段）的阿坝—若尔盖—临潭—兰州大地电磁剖面（WQL-L1）所采集到的数据进行了精细化处理分析和二维反演研究,结合跨过2013年岷县漳县地震区的WQL-L6剖面大地电磁探测结果和以往的地质与地球物理资料,对西秦岭造山带（中段）的深部电性结构、主要断裂带延伸状况以及与南北两侧地块的接触关系等进行了分析研究,结果表明：东昆仑断裂带塔藏段、迭部—白龙江断裂和光盖山—迭山断裂带共同组成了东昆仑断裂系统,分隔了松潘—甘孜地块和西秦岭造山带（中段）;西秦岭北缘断裂带为主要的高角度南倾大型电性边界带,延伸深度穿过莫霍面;临潭—宕昌断裂带具有电性边界带特征,其延伸情况具有东、西差异.西秦岭造山带（中段）自地表到深度约20km范围表现为东北和西南浅、中部深的倒＂梯形＂高阻层,在高阻层之下广泛发育低阻层,低阻层与高阻层相互契合,呈现相互挤压堆积的式样,其西南侧的松潘—甘孜地块中下地壳存在西南深、东北浅低阻层,其东北侧的陇西盆地具有稳定的成层性结构,显示出西秦岭造山带（中段）正处于松潘—甘孜地块向北挤压和陇西盆地向南的阻挡挤压作用中.松潘—甘孜地块从西南向东北推挤、东北侧陇西盆地相对阻挡的相互作用是2013年岷县漳县6.6级地震发生的外部动力学机制,同时地震震源区特殊介质属性是该次地震发生的内部因素.西秦岭造山带（中段）中上地壳倒＂梯形＂高阻体埋深西薄、东厚的分段差异与该段内部中强地震分布差异有关.东昆仑断裂玛沁段和塔藏段内部的深部电性结构差异和延伸状况与东昆仑断裂自西向东走滑速率减小有内在联系.
