三维模式约束的电离层掩星反演方法

目前电离层掩星数据反演是基于电离层电子密度分布局部球对称近似的Abel反演方法,实际电离层的非球对称性会给电子密度反演结果带来误差.本文研究利用三维电离层模式来提供电子密度水平变化的先验信息约束电离层掩星反演的方法,即三维模式约束法;并将该方法应用于模拟掩星观测数据和实测掩星数据的反演.模拟观测数据的反演结果表明,与Abel反演方法相比,三维模式约束法能够减小反演误差.采用IRI2001模式作为约束,对COSMIC电离层掩星实测数据反演,将反演结果与全球的垂测仪数据进行比较,结果表明,三维模式约束法和Abel反演方法都能很好地反演电离层掩星.     

电离层GPS掩星反演技术研究

GPS无线电掩星技术是崭新的、高效的地球大气层和电离层探测技术,但仍在发展和完善之中.本文详细推导了Abel积分和绝对TEC电离层反演方法,研究了如何解决Abel积分产生的上下限异常问题;用COSMIC发布的GPS原始数据进行了反演计算,将结果与地面电离层测高仪数据进行了比较,最后讨论了周跳对反演结果的影响问题.结果表明：（1）在较高轨道高度（约800 km）,Abel积分与绝对TEC方法的反演结果基本一致,都与电离层测高仪反演结果符合良好;在较低轨道高度（约500 km）,绝对TEC反演精度优于Abel积分反演精度;（2）绝对TEC反演的最大电子密度N_m较Abel积分法获得的结果更接近于电离层测高获得的峰值电子密度N_mF₂,绝对TEC反演法更加严密和有效;（3）周跳对绝对TEC反演结果的影响较Abel积分反演结果的影响更为敏感,但无论哪种方法,周跳对反演精度都造成严重损失.综合而言,绝对TEC反演法是更优的方法.     

地基GPS VTEC约束的电离层掩星反演方法

本文在分析GPS电离层掩星Abel积分变换反演的基础上,介绍了一种顾及电离层水平梯度信息的反演新方法,即地基GPS VTEC约束的电离层掩星反演法,并将该方法应用于COSMIC低纬掩星观测资料的反演.初步的统计结果表明,与经典Abel变换反演相比,地基GPS VTEC约束反演的临界频率f_oF₂更接近于电离层测高仪观测的结果.     

GNSS掩星中大气水汽的非线性反演

研究利用无线电掩星探测数据反演大气水汽参数,文中使用一维非线性方法,反演对流层水汽压廓线.首先,定义了一维目标函数,并以几何光学假设为前提,反演出电波弯曲角;将电波弯曲角作为掩星探测量,大气模式输出的温、湿参量作为初始场,代入目标函数,对目标函数求最优反演出水汽压廓线.文中给出了部分反演结果,经讨论分析认为:非线性方法反演水汽压对初始场的精度依赖较小,能够同时反演出大气温度\水汽压廓线及相应的误差分析,对数值天气预报的应用及数据同化研究具有重要意义.     

地基GNSS VTEC约束的大气掩星电离层误差修正方法

介绍了一种考虑电离层沿GNSS掩星射线路径分布不对称且兼顾一阶和二阶项的大气掩星电离层误差修正新方法,它综合地基GNSS VTEC的水平变化信息和电离层模式的垂直变化信息,在沿"入射线"与"出射线"双边局部球对称假设下,估算GNSS掩星射线路径上的电子密度;进而计算一阶和二阶项弯曲角电离层误差廓线.采用太阳活动低年的2008年7月15日和太阳活动较高年的2013年7月15日两天的MetOp-A和GRACE掩星观测资料和IGS的GNSS VTEC数据,计算了GNSS大气掩星弯曲角一阶和二阶项电离层误差廓线.对比分析表明：新方法经验模型和理论模型的二阶项弯曲角电离层残差,以及经验模型与实测数据的一阶和二阶项弯曲角电离层误差均具有良好的一致性,因此该方法可用于L2信号数据质量较差或L2信号中断的掩星事件,同时修正大气掩星一阶和二阶电离层误差,从而提高弯曲角精度和掩星观测资料的利用率.     

电离层GPS掩星观测反演技术

在电离层局部地区球对称假设下,推导了利用双频和单频无线电掩星观测数据,反演电离层电子密度剖面的两种方法. 双频反演的误差来自于载波相位的观测误差,单频反演误差则主要由伪距的观测精度决定. 由于载波相位测量精度比伪距测量精度高两个量级,因此双频反演的精度一般比单频反演的高些. 不过,两载波信号L1和L2之间的传播路径差异会给双频方法带来误差. 利用三维射线追踪的程序模拟的无线电掩星数据来评估这些方法,结果表明,反演出的电离层剖面与给定的模式电离层非常吻合,验证了两种方法的可靠性和准确性. 将这两种反演方法应用于处理实测的GPS/MET掩星观测数据,均能获取合理的电离层剖面信息. 且单频方法得到的反演剖面与双频方法相当一致, 这为利用LEO星载单频GPS接收机进行电离层掩星观测提供了理论基础.     

Abel电离层掩星反演方法及误差分析

GNSS-LEO电离层无线电掩星技术是近年来发展的电离层探测新技术.为消除LEO轨道以上的电离层影响,改正TEC反演方法采用非掩星侧的观测数据进行电离层掩星反演.本文首次提出了一种新方法--基于历元差分的电离层反演方法;并将改正TEC与历元差分两种反演方法应用于模拟掩星观测数据反演,随后基于反演结果及误差分析得到一些有益的结论:历元差分反演精度较改正TEC反演精度均有所提高;不管是哪种方法,高轨(约800 km)反演结果要优于低轨(约500 km)的反演结果;随着剖面高度的降低,反演精度随之下降;反演误差主要集中在8至18时(当地时),主要分布在磁纬-30°至30°之间.     

大气掩星反演误差特性初步分析

GPS大气掩星探测技术可以获得全球大气折射率、气压、密度、温度和水汽压等气象参数，该技术基本原理是基于几何光学近似的Abel积分反演.地球扁率、电离层传播时间延迟、大气大尺度水平梯度、多路径传播现象等因素在某些高度范围影响大气反演的精度.本文采用模拟的方法，分析其中地球扁率及电离层对反演结果的影响，并讨论局部圆弧修正及电离层修正的效果.利用CHAMP掩星实测轨道数据和有关电离层和大气经验模式、采用三维射线追踪方法模拟计算几种情形下的GPS掩星观测附加相位数据，对模拟数据进行反演，将反演气象参量剖面与模拟时给定模式剖面进行比较，得到了0～60 km高度范围内的反演误差.误差统计分析结果表明，局部圆弧中心的修正以及电离层修正，对于高精度的GPS掩星反演是非常重要的；电离层修正残差仍是制约30～60 km高度范围内反演精度的重要因素，进一步完善和优化大气掩星反演需要发展新的电离层修正算法.     

电离层GPS掩星观测改正TEC反演方法

电离层掩星观测中，当低轨卫星（LEO）轨道高度较低时，轨道以上的电离层电子总含量（TEC）对掩星反演的影响不能忽略.目前，一般采用指数函数等外推方法来处理该问题，对反演结果可能引起较大误差.为提高电离层掩星反演精度，本文研究利用LEO处于非掩星一侧GPS观测数据的改正TEC新反演方法.用三维射线追踪程序计算出电离层掩星观测模拟数据，分别应用改正TEC方法和外推方法进行反演，将反演结果与所用模式值进行比较.结果表明：对于轨道高度约800km的GPS/MET掩星模拟数据，外推方法和改正TEC方法反演结果都与模式值基本一致；对于轨道高度约400km的CHAMP掩星模拟数据，外推方法误差较大，改正TEC方法反演结果与模式值相符得较好.将改正TEC方法应用于GPS/MET实测数据的反演，取得了合理的结果.这些说明，改正TEC算法是一种有效的电离层掩星反演方法，尤其是对于轨道较低的LEO的电离层掩星观测反演特别有用.     

GPS地面台网和掩星观测结合的时变三维电离层层析

本文给出GPS地面台网和掩星观测结合的时变三维电离层层析的原理、算法和基于实测数据的反演结果.反演结果的比较表明,联合地基GPS与掩星观测数据进行重建,电子密度整体图像的重建质量特别是其垂直结构的重建质量得到了明显改善.在平静日和磁暴期间两种条件下利用实测数据的重建结果表明,GPS地面台网和掩星观测结合的电离层层析可以获得电离层电子密度在高度-纬度-经度-时间四维空间中的变化.重建结果清晰地显示了磁暴期间电离层负相暴效应,表明结合GPS地面台网和掩星观测的时变三维电离层层析可以有效地监测扰动条件下的大尺度电离层结构.     

模型约束三维折射静校正方法研究

本文详细的介绍了模型约束三维折射静校正技术的基本原理和计算步骤,具体给出了三维抽线,目的层质量监控,三维速度内插,三维底界内插等技术的实现方法,并结合了西部某工区的具体实例加以说明.结果表明,这种由微测井、小折射求取近地表深度-时间速度模型,由大炮记录折射初至旅行时求取各点延迟时,两者联合反演低降速层底界的厚度和高程,最终求取低降速层静校正量的方法,可明显提高低降速带的反演精度,较好地解决表层结构复杂地区的静校正问题.     

大气GPS掩星观测反演方法

大气GPS掩星观测可获得全球的大气气象参量剖面信息.本文阐述了地球大气GPS掩星观测反演原理,详细介绍了其几何光学近似反演方法和全谱反演方法,提出了将几何光学反演方法和全谱反演方法结合以形成可以处理多路径掩星数据的反演新方案.该方案和几何光学反演方法应用于GPS/MET和CHAMP大气掩星数据反演,成功地获得了大气参量剖面.结果表明,新反演方案是可行的、有效的GPS大气掩星反演方案.     

地震危险性估计的一种尝试

本文根据三角网统一平差和跨断层测量资料，利用位错理论的反演方法，计算了唐山和北京，怀来地区几条主要断层滑动速率和滑动方向，其结果与我们用有限元方法计算的断层应力场配套。并获得三点结果：（１）断层滑动速率大，则断层上庆力高，反之亦然；（２）断层滑动方向与断层应力场方向配套；（３）当跨断层短线大度年变化持续出现线笥关系时，断层浅部与深部滑动一致，而且断层浅部滑速率大，此却表示有发生地震的危险性。     

Comments on “Steady- and transient-state inversion in hydrogeology by successive flux estimation” by P. Pasquier and D. Marcotte

Pasquier and Marcotte [Pasquier P, Marcotte D. Steady- and transient-state inversion in hydrogeology by successive flux estimation. Adv Wat Res 2006;29:1934–52] propose some modifications to the Comparison Model Method (CMM), in order to apply it to transient 3D ground water flow data for conductivity identification. We present some remarks on that paper to improve the comprehension of the basic features of the CMM and of the real value of the novelties introduced by Pasquier and Marcotte.     

大地电磁资料精细处理和二维反演解释技术研究(三)——构建二维反演初始模型的印模法

在主流的线性最优化大地电磁二维反演中,如何合理构建初始模型是一个亟待解决的问题.常用的是采用均匀半空间或一维反演结果构建初始模型,不易获得稳定可靠的反演效果.实践表明,尽管基于不同初始模型的大地电磁二维反演结果差别较大,但均较初始模型更为接近真实模型.基于这样一种认识,经过反复的理论和实践探索,我们提出构建大地电磁二维反演初始模型的印模法.印模法的基本思想是依据已有反演结果和均匀半空间模型之间的加权来确定下一步二维反演的初始模型,它一方面保留了已有反演结果中关于真实模型的宏观轮廓信息,另一方面,保证了深部电性结构的均匀性,从而满足大地电磁二维正演所要求的底边界条件.基于印模法,本文进一步提出了迭代重构的反演思想.通过多个理论模型和实测数据的反演计算,验证了上述方法可在很大程度上压制初始模型对反演结果的影响.     

基于广义模型约束的时间域航空电磁反演研究

由于航空电磁具有海量数据,因此快速有效的成像和反演手段至关重要.本文针对层状介质模型推导与实现了广义模型约束条件下时间域航空电磁一维反演.从正则化反演的目标函数出发,通过改变模型约束项构造Lp范数反演和聚焦反演,进而通过改变模型求解域构造出基于小波变换的稀疏约束反演.针对不同反演方法目标函数的构建方式,本文进一步从数学原理上分析不同反演方法的预期效果,并通过理论模型和实测数据进行验证.结果表明L0.8范数反演、聚焦反演和基于小波变换的稀疏约束反演可以得到更符合地下层状介质陡变界面的反演结果.     

可控源音频大地电磁数据的反演方法

从反演方程、构造目标函数和求解三方面对用于可控源音频大地电磁法(CSAMT)的实用反演方法中的四种进行了描述．水平层状地层CSAMT法资料的直接反演法首次尝试了一维空间的全资料CSAMT反演，效果较好，但该方法尚难应用于2D、3D复杂介质中；奥克姆反演方法既考虑了横向的光滑函数，又考虑了纵向的光滑函数，得到比较光滑的横向、纵向变化的背景电性结果，但有可能把一些小构造光滑掉．快速松驰反演算法和共轭梯度算法由于计算速度快，占内存少而被用于三维反演中，二者相比，快速松驰算法在求解雅可比矩阵时只做一次正演计算，在更新模型时解小型方程组，所以在速度上更胜一筹．在后三种算法中，由于复杂电性结构无解析解，正演计算都采用数值计算．数值计算的可靠性、速度影响着反演算法的有效性，这方面的研究也将是2D、3D复杂电性结构反演的研究方向之一.     

多道瞬变电磁法共中心点道集数据联合反演

多道瞬变电磁法是目前地面资源勘查领域的研究热点,开展该方法观测数据的一维反演对完成装备研发及技术推广具有积极的推动作用.针对多道瞬变电磁法特殊的源波形和数据采集方式,本文对多道瞬变电磁数据的响应提取和一维反演进行了研究.首先,采用反卷积方法从伪随机二进制序列全波形响应数据中提取大地脉冲响应.为了避免脉冲响应中空气波的干扰,提出采用反向积分算法,把脉冲响应转换成下降沿阶跃响应,并将其整理成共中点道集.最后,为充分利用各个偏移距下的响应数据对不同深度范围目标体的分辨能力,采用基于OCCAM算法的联合反演方法对共中心点道集数据进行联合反演.结果表明在不同目标体深度以及有噪声干扰的情况下,共中心点道集数据联合反演计算均可获得较好的结果.     

反射地震数据的逐层波形反演

本文针对层状介质并结合梯度法波形反演，提出逐层波形反演的方法. 首先给出介质扰动响应的概念，并在此基础上分析了梯度法波形反演方法. 波形反演实质上是将实测地震记录和预测地震记录的波形残差信息转化为实际地质模型和预测地质模型的模型残差信息. 波形反演的优点是利用大量振幅相位信息得到高分辨率的反演结果， 其缺点是运行耗时大；当初始模型和实际模型相差较大时，迭代算法容易陷入局部极小点，这是因为目标函数和初始模型同实际模型间的差异是非线性的关系. 逐层波形反演方法是使自上而下每一层的目标函数最小，这样总的目标函数也是最小的. 利用二分法速度扫描确定每一层速度不仅提高了运算速度也避免了迭代算法陷入局部极小点的问题. 结合介质扰动响应和目标函数值变化可以更为准确迅速地确定每一层速度和该层界面位置.     

地球物理学中的电磁场正演与反演

本文在近年来众多的地球物理研究者的研究基础上，总结了当前地电磁模型正反演已有成果，定量分析了各种主要正反演的性能测试，指出不同正反演法的优点、缺点以及应用范围局限，提出了各种方法的发展趋势以及当前计算地球物理领域的核心内容，指出了计算地球物理领域的数值模拟发展方向.