利用p~3软件进行GPS单点定位解算的若干问题

利用p3软件对大量实测数据进行了高精度单点定位有关问题的验证和分析;在IGS提供的精密星历和卫星钟差产品中,较为深入地比较分析了快速产品和几种最终产品的定位精度、收敛速度及p3软件的正反算结果精度。得出以下结论:快速产品与最终产品的定位精度和收敛速度相似;COD的收敛速度和定位精度最高;p3软件的反算精度明显优于正算,且精度也比较均匀。这些研究成果对PPP技术的实际应用具有较好的借鉴意义。     

Weak-contrast approximation of the elastic scattering matrix in anisotropic media

The plane-wave reflection and transmission coefficients at a plane interface between two anisotropic media constitute the elements of the elastic scattering matrix. For a 1-D anisotropic medium the eigenvector decomposition of the system matrix of the transformed elasto-dynamic equations is used to derive a general expression for the scattering matrix. Depending on the normalization of the eigenvectors, the expressions give scattering coefficients for amplitudes or for vertical energy flux.Computing the vertical slownesses and the corresponding polarizations, the eigenvector matrix and its inverse can be found. We give a simple formula for the inverse, regardless of the normalization of the eigenvectors. When the eigenvectors are normalized with respect to amplitudes of displacement (or velocity), the calculation of the scattering matrix for amplitudes is simplified.When the relative changes in all parameters are small, a weak-contrast approximation of the scattering matrix, based on the exactly determined polarization vectors in an average medium, is obtained. The same approximation is also derived directly from the transformed elasto-dynamic equations for a smooth vertically inhomogeneous medium, proving the consistency of the approximation.For monoclinic media, with the mirror symmetry plane parallel to the interface, the approximative scattering matrix is given in terms of analytic expressions for the non-normalized eigenvectors and vertical slownesses. For transversely isotropic media with a vertical axis of symmetry (VTI) and isotropic media, explicit solutions for the weak-contrast approximations of the scattering matrices have been obtained. The scattering matrix for amplitudes for isotropic media is well known. The scattering matrix for vertical energy flux may have applications in AVO analysis and inversion due to the reciprocity of the reflection coefficients for converted waves.Numerical examples for monoclinic and VTI media provide good agreement between the approximative and the exact reflection matrices. It is, however, expected that the approximations cannot be used when the symmetry properties of the two media are very different. This is because the approximation relies on a small relative contrast between the eigenvectors in the two media.Presented at the Workshop Meeting on Seismic Waves in Laterally Inhomogeneous Media, Castle of Trest, Czech Republic, May 22–27, 1995.     

鄂尔多斯盆地苏里格庙气田上古生界盒8段地质与地球物理响应特征

苏里格庙气田主要产气层为上古生界石盒子组盒8段，属于三角洲平原分流河道沉积。复合砂体近南北向展布；单砂体的厚度变化大，非均质性明显，储渗变化大。盒8段含气砂岩在AVO资料上有明显反映；在综合测井资料上，盒8段表现出自然伽玛测井值、密度测井值、补偿中子测井值都偏低，电阻率、声波时差测井值较高，出现了大幅度的自然电位异常。利用核磁共振测井技术测得苏里格庙区某井盒8段3555．9m～3559．6m为气层，3560．6m～3572m为气水层，与生产实践吻合得很好。     

Policy options to streamline the carbon market for agricultural nitrous oxide emissions

The majority of emissions of nitrous oxide – a potent greenhouse gas (GHG) – are from agricultural sources, particularly nitrogen fertilizer applications. A growing focus on these emission sources has led to the development in the United States of GHG offset protocols that could enable payment to farmers for reducing fertilizer use or implementing other nitrogen management strategies. Despite the development of several protocols, the current regional scope is narrow, adoption by farmers is low, and policy implementation of protocols has a significant time lag. Here we utilize existing research and policy structures to propose an ‘umbrella’ approach for nitrogen management GHG emissions protocols that has the potential to streamline the policy implementation and acceptance of such protocols. We suggest that the umbrella protocol could set forth standard definitions common across multiple protocol options, and then modules could be further developed as scientific evidence advances. Modules could be developed for specific crops, regions, and practices. We identify a policy process that could facilitate this development in concert with emerging scientific research and conclude by acknowledging potential benefits and limitations of the approach.

Key policy insights

• Agricultural greenhouse gas market options are growing, but are still underutilized

• Streamlining protocol development through an umbrella process could enable quicker development of protocols across new crops, regions, and practices

• Effective protocol development must not compromise best available science and should follow a rigorous pathway to ensure appropriate implementation

     

实时卫星钟差基准精化方法

在基于精密单点定位(PPP)的授时方法中,卫星钟差产品的高精度时间基准至关重要.针对实时卫星钟差产品时间基准不够稳定的问题,本文采用一组具有原子钟外部输入的国际全球卫星导航系统(GNSS)服务(IGS)跟踪站建立了顾及原子钟变化特性的基准精化方法.该方法首先采用阿伦方差对不同的IGS跟踪站外接原子钟进行稳定度分析,挑选出一组稳定度高的原子钟用以精化时间基准.在此基础上,利用阿伦方差分析各台原子钟的噪声参数特征,并确定不同原子钟之间的权比关系.最终,建立时间基准改正量的随机模型,并计算出精化后的时间基准.通过实例验证表明:与IGS事后精密钟差产品定义的时间基准比较,改正后的实时钟差基准单天内的标准差(STD)优于0.1 ns,相比于改正前最高提升了93%.同时,基准改正后的天内万秒稳达到10^-15量级,实现了一个量级的提高.此外,通过相对钟差精度的分析,表明钟差基准修正不影响PPP的定位精度.     

BDS卫星钟差半参数平差模型异常数据探测与处理

针对北斗卫星钟差数据中出现钟跳、粗差等异常数据频次较高的问题,将半参数平差模型引入BDS的钟差数据的处理过程中。首先,在考虑系统误差的同时,改进了常用的钟差模型;其次,综合考虑异常数据和系统误差,利用补偿最小二乘原理和极值求解方法,详细推导了分离异常数据前后参数和非参数估计值与相应观测值改正数的关系表达式,实现了异常数据的定值、参数求解和系统误差分离。在此基础上,引入Cook距离,给出了利用参数分量和非参数分量的Cook距离及混合Cook距离,去判断异常值的位置,并给出了一些参数的选取方法及相应的处理措施,实现了卫星钟差异常数据的定位和定值以及部分系统误差的分离。最后,采用武汉大学GNSS中心提供的采样间隔为5 min的北斗卫星精密钟差数据,将本文方法与常用方法进行了试验对比。试验结果表明,本文方法能够有效地识别并处理卫星钟差数据中的异常值,有效克服了基于经验阈值钟差异常数据探测方法的不足,且该方法对于量级较小的异常钟差数据也有很好的探测效果,一定程度丰富了现有的BDS钟差数据质量控制方法。     

利用GNSS测量射电望远镜参考点的仿真分析

高精度测量射电望远镜参考点和轴线偏差等参数,对建立天线指向模型和本地连接参数、提高测站坐标精度等具有重要意义。为完成新建射电望远镜参考点初始参考值的快速测定,根据望远镜的旋转模型,结合常规静态归心测量方法和随机动态测量方法,提出了一种利用GNSS天线代替测量靶标实现望远镜参考点测量的方法。通过仿真分析验证了该方法的合理性和有效性,并分析了数据点个数和数据点测量精度对天线参考点和轴线偏差解算精度的影响。     

等效偏移距方法及应用

当勘探目标地质构造复杂、地层倾角较陡或岩性横向突变时,在地表接收到的反射信号较弱、信噪比较低,甚至接收不到反射波,常规偏移成像方法对这类数据很难精确成像,因此可以考虑基于散射理论的等效偏移距方法来成像.该方法分两步进行：第一步是抽取共散射点道集;第二步是对共散射点道集进行克希霍夫积分偏移.文中数值模拟及实际资料的实例表明,该方法在复杂地质目标及低信噪比地震资料的成像处理中具有一定的实用价值.     

Offset dependent amplitudes: Fault curvature effects and dip move-out corrections

Equations have been presented previously which predict that reflector curvature can affect significantly seismic reflection amplitudes at both zero and nonzero source-receiver offsets. Here the fact that faults are generally concave-upward is used to examine the curvature effect for compaction-driven faults which the sediments have both exponential and logarithmic porosity decreases with increasing depth.The curvature effect is generally larger over upwardly concave reflectors where amplitudes can either increase with offset (exposed focus) or decrease with offset (buried focus). The magnitude of the effect depends on the ratio between the depth to the structure and the radius of curvature of the structure. A phase change of 90° also occurs at a critical offset in the case of an exposed focus concave-upward reflector, with decreasing amplitude at offsets larger than the critical value. Dip move-out (DMO) dominantly removes the amplitude variation with offset due to curvature. These results suggest that when looking for amplitude variations with offset in a fault prospect, DMO should be applied as a preprocessing step. Compaction-driven faults have an exposed focus and, for an exponential variation of porosity with depth, there is a maximum curvature effect at a depth roughly the same as the scaling depth for the porosity. Logarithmic porosity with depth variations suggest a continued increase in the curvature effect of faults with increasing depth, to a maximum amplitude increase of about 35% relative to a plane interface.     

温纳偏置电测深法压制表层电性不均匀体影响的研究

巷道底板浮煤等电性不均匀体及厚度变化是造成煤矿井下电测深曲线畸变的主要因素。通过模型实验研究了温纳偏置电测保法压制浮煤等表层电性不均匀体影响的能力。用感度分析法讨论了表层电性不均匀体产生地电干扰的机制。实测结果表明，该方法改善了电测深曲线质量，有效地提高了电测深法地质应用效果。