油气勘探开发数据体平均值代表性分析

油气勘探开发涉及各类数据体,同一数据体可得多种平均值,目前尚无明确而有效的方法以判断何种平均值能客观反映数据体的典型水平。运用数据统计平均分析方法,对勘探开发实践中的孔隙度、渗透率、产量、成本等数据体进行系统分析,确立了加权中位数计算公式和平衡中位数法则。加权中位数计算公式适于分析不同领域正常有序数据体的基本特征,加权中位数为正常有序数据体的平衡点;平衡中位数法则适于确定有足够大数据容量、能满足统计分析基本要求、能选择合理权衡指标的正常有序数据体的典型水平;有明确物理意义的加权平均值也可确定数据体的典型水平。     

三维坐标转换的稳健加权总体最小二乘算法

针对当前加权对称相似变换中仅考虑观测值含有随机误差,而不考虑观测值含有粗差的情况,基于加权对称相似变换进一步验证其不具有抗差性,并在此基础上采用选权迭代法进行稳健加权对称相似变换。由中位数法求解具有抗差性的单位权中误差,构造标准化残差,求解权因子,进而求得可靠的参数解。实验结果表明,当观测值含有4~6个粗差时,采用该方法能够探测出更多可能存在粗差的数据,给出的权因子更合理,所得参数解精度更高、稳定性更强。     

地壳形变分析的抗差最小二乘配置迭代解法

利用最小二乘配置进行地壳形变分析,其结果的合理性关键在于经验协方差函数的拟合.考虑到观测数据存在粗差的情况,提出基于观测值中位数初值的抗差最小二乘配置方法和基于中位参数法的抗差最小二乘配置方法.两种方法首先分别利用观测值中位数给出观测值初始权阵以及利用中位参数法给出最小二乘配置初始解,然后均在给定协方差函数参数初始值的情况下,应用合适的等价权进行抗差估计并通过迭代计算,最终获得稳健的协方差函数参数估值及最小二乘配置解.利用本文提出的两种方法以及传统方法分别对庐山地震的GPS垂直位移数据和意大利L'Aquila地震的InSAR同震位移数据进行处理分析.结果表明:相对传统方法,基于观测值中位数初值的抗差最小二乘配置方法效果更好,更具稳健性.     

利用中位数进行光纤陀螺信号抗差估计

为了削弱光纤陀螺信号中异常噪声对寻北精度的影响,在滤除高频周期噪声的基础上,提出了一种处理陀螺信号的新算法。推导了陀螺信号的抗差解公式及其误差影响函数,取观测值的中位数作为抗差初值,利用一次抗差估计求得的观测残差再用中位数法求得均方差因子,并给出了高崩溃污染率的初值辅以IGGⅢ方案迭代解算的混合算法。计算结果表明,基于可靠初值和均方差因子估值的抗差估计可以有效削弱异常干扰对陀螺的影响。     

一种基于GNSS的可驯频标的设计与实现

重点介绍了GNSS可驯频标软件设计方法。在软件设计中,采用了滑动中位数方法以滤除测量数据中的奇异值;利用了Kalman滤波器对测量时差数据进行滤波以抑制测量噪声;提出了一种改进型的压控电压产生方法。采用GPS定时接收机1 PPS信号作为参考基准,对原来准确度为1×10-9恒温晶振进行驯服,驯服后频率准确度优于2×10-11。     

窄巷FCB估计方法改进及时变特性分析

窄巷波长相对较短,受观测数据质量和未模型化误差的影响较大,导致参考站计算的窄巷FCB估值的互差较大,窄巷FCB稳定性较差。针对窄巷FCB估值出现较大偏差的参考站,本文提出了基于抗差初值的窄巷FCB抗差估计方法,该方法取当前历元所有参考站窄巷FCB估值的中位数作为抗差初值,采用IGGⅢ方案降低了有偏差的窄巷FCB估值的权比。采用IGS全球跟踪站网的数据进行验证,改进的方法提高了窄巷FCB的精度和稳定性,采用每天估计的一组卫星端的窄巷FCB即可满足窄巷模糊度固定的精度要求。同时,采用本文改进方法估计的窄巷FCB,可以提高用户端仿动态单天解的窄巷模糊度固定成功率和定位精度。     

L1-范估计的无偏性

根据一元Laplace分布的概率密度函数及分布函数公式,推导了当观测值个数n为偶数时,中位数的概率密度函数;在此基础上证明了L1-范估计的无偏性.     

抗差加权整体最小二乘模型的牛顿-高斯算法

基于加权整体最小二乘的牛顿-高斯迭代算法,提出了一种抗差加权整体最小二乘模型。利用标准化残差构造权因子函数,并采用中位数法获得具有抗差性的单位权中误差估值,能同时实现观测空间和结构空间抗差。为获得标准化残差,利用线性近似的协因数传播律推导了加权整体最小二乘残差协因数阵的表达式,并给出模型的迭代计算方法。试验结果表明:对于加权整体最小二乘的粗差处理问题,本文提出的方法具有良好的抗差性能,参数估值与不含粗差时加权整体最小二乘的结果没有显著的差异,性能优于直接由残差构造的稳健加权整体最小二乘模型。     

小波变换辅助的陀螺信号抗差估计

在捷联寻北仪中,陀螺漂移是导致寻北精度下降的主要原因.在小波阚值消噪的基础上,提出了利用抗差估计提高光纤陀螺寻北仪的寻北精度,并给出了陀螺信号的抗差解式.首先对陀螺输出信号进行频谱分析,确定多分辨分析的尺度,对不同尺度的高频系数采取相应的处理措施;然后将抗差估计应用于陀螺数据的处理,采用高崩溃污染率的初值辅以IGGⅢ方案迭代解算的混合算法.计算结果表明,该方法能够有效地抵制异常扰动的影响,提高寻北精(密)度.     

一维参数期望估计

本文给出一维参数期望估计的简易算法。并探讨了在外业观测中的应用问题。