L1范数与IQR统计量组合的GNSS坐标序列粗差探测算法

GNSS坐标时间序列中不可避免地含有粗差,未剔除的粗差将会导致参数估计有偏。因此,粗差探测与剔除是GNSS坐标序列分析中一项重要的数据预处理工作。针对GNSS坐标时间序列特点,提出了一种将L1范数(L1-norm)估计与四分位距统计量IQR(interquartile range)组合的移动开窗粗差探测算法,称之为L1_Mod IQR。该方法的主要思想是,首先利用L1范数估计得到较"真实"的残差,然后再对残差采用IQR统计量进行粗差探测。将L1_Mod IQR法与"3σ"法、基于最小二乘的τ检验法等粗差探测算法进行了模拟计算与对比,验证了该算法的有效性。进一步采用L1_Mod IQR算法对中国区域10个IGS站的高程时间序列进行了分析,结果表明中国区域IGS站高程序列的粗差剔除率最小为0.1%,最大为2.6%。并且以WUHN站为例与SOPAC提供的结果进行了对比,结果表明SOPAC提供的"Clean"数据仍含有大量的粗差,而L1_Mod IQR算法能够有效地剔除粗差。     

Provenance of Holocene beach sand in the Western Iberian margin: the use of the Kolmogorov–Smirnov test for the deciphering of sediment recycling in a modern coastal system

Detrital zircons from Holocene beach sand and igneous zircons from the Cretaceous syenite forming Cape Sines (Western Iberian margin) were dated using laser ablation – inductively coupled plasma – mass spectrometry. The U–Pb ages obtained were used for comparison with previous radiometric data from Carboniferous greywacke, Pliocene–Pleistocene sand and Cretaceous syenite forming the sea cliff at Cape Sines and the contiguous coast. New U–Pb dating of igneous morphologically simple and complex zircons from the syenite of the Sines pluton suggests that the history of zircon crystallization was more extensive (ca 87 to 74 Ma), in contrast to the findings of previous geochronology studies (ca 76 to 74 Ma). The U–Pb ages obtained in Holocene sand revealed a wide interval, ranging from the Cretaceous to the Archean, with predominance of Cretaceous (37%), Palaeozoic (35%) and Neoproterozoic (19%) detrital‐zircon ages. The paucity of round to sub‐rounded grains seems to indicate a short transportation history for most of the Cretaceous zircons (ca 95 to 73 Ma) which are more abundant in the beach sand that was sampled south of Cape Sines. Comparative analysis using the Kolmogorov–Smirnov statistical method, analysing sub‐populations separately, suggests that the zircon populations of the Carboniferous and Cretaceous rocks forming the sea cliff were reproduced faithfully in Quaternary sand, indicating sediment recycling. The similarity of the pre‐Cretaceous ages (>ca 280 Ma) of detrital zircons found in Holocene sand, as compared with Carboniferous greywacke and Pliocene–Pleistocene sand, provides support for the hypothesis that detritus was reworked into the beach from older sedimentary rocks exposed along the sea cliff. The largest percentage of Cretaceous zircons (<ca 95 Ma) found in Holocene sand, as compared with Pliocene–Pleistocene sand (secondary recycled source), suggests that the Sines pluton was the one of the primary sources that became progressively more exposed to erosion during Quaternary uplift. This work highlights the application of the Kolmogorov–Smirnov method in comparison of zircon age populations used to identify provenance and sediment recycling in modern and ancient detrital sedimentary sequences.     

朝不楞多金属矿外围1∶5万地电化学勘查应用

运用地电化学技术在内蒙古半干旱草原风成砂浅覆盖区朝不楞多金属矿外围开展了1∶50000面积性勘查工作,并与常规土壤测量开展了对比研究。结果表明:(1)从元素异常强度、衬度、规模及与已知矿体空间对应性等方面对比,地电化学测量比常规土壤测量具有显著优越性;(2)地电化学测量在区域上显现出Ag-Pb-Zn-Bi-Cd及Cu-Co-Cr-Fe-La-Mo-Ni两套元素异常组合,与已知见矿钻孔中探测到的矽卡岩型铁矿及热液型铅锌矿均有出现具有较好的对应性;(3)在工区东部及东南部未知区发现多处异常强度、规模较大的多元素组合异常,建议开展进一步的勘查工作,以实现该矿区找矿突破。     

松辽盆地大陆科钻二开段大井眼钻井液技术

CCSD-SK2井是一口钻穿松辽盆地白垩纪陆相沉积地层的中国大陆科学钻探井。二开段先以216mm口径钻至2806.20m,然后进行444.5mm口径扩孔。泥页岩地层大井眼钻进过程中粘土矿物含量高且易造浆,砂砾岩层渗透性强易粘卡,井壁稳定控制与钻屑携带难度大。选用了钾铵聚合物钻井液体系和聚磺钻井液体系,通过正交实验确定了钻井液配方。将钻井液控制在合适的密度、较低的失水量以及较强的封堵性,有效实现了井壁稳定;保持适当的环空返速、较高的动塑比以及定期打稠塞举砂等措施,有效实现了井眼净化。1086.45m~1147.56m与1182.74m~1256.01m井段累计收获岩心130.90m,取心率达到97.41%,平均机械钻速为0.84m/h;最终顺利钻至2806.20m。     

中国及周边海域对流云团的水平和垂直尺度

利用2007年1月-2010年12月CloudSat-CALIPSO二级云产品2B-CLDCLASS-LIDAR,统计中国及其周边海域对流云的发生频率,根据对流云发生频率的分布特征将中国及周边海域划分为青藏高原（TP）、东部陆地（EC）、南部海域（SO）和西北太平洋（WP）4个子区域,并研究了4个子区域积云团和深对流云团的水平尺度和垂直尺度。统计结果表明,海洋积云团的水平尺度约为2 km,陆地积云团的水平尺度约为1 km,海洋下垫面热力性质均匀,积云团尺度更大;陆地下垫面非均匀性强,积云团分布更为零散。深对流云团的水平尺度为10-50 km,东部陆地最大,约为45 km,西北太平洋最小,约为30 km。陆地深对流云团水平尺度较海洋上大,且多尺度特征显著,应该与深对流云发生的复杂天气背景有关。积云团的垂直尺度范围为0.24-2 km,4个区域无明显差异。垂直尺度海洋深对流云团大于陆地云团,其中在南部海域地区最大,约为15 km,青藏高原最小,约为10 km。与陆地云团相比,海洋深对流云团表现为水平尺度更小、垂直尺度更大的中尺度对流体特征。      

天山南麓雷暴日短期预报研究

利用阿克苏及邻近地区12个气象站1980—2013年雷暴资料,以及同期高空资料,统计了各站年均雷暴日数,对发生区域雷暴天气的环流形势进行分类,归纳出各型的入型指标。通过逐步回归法,建立阿克苏及邻近地区区域雷暴概率回归预报模型,并对2013年进行试预报。结果表明:(1)阿克苏及邻近地区区域雷暴的影响系统主要分为4类:巴湖低槽型、急流型、西北气流型和温度槽型。(2)对2002—2012年5—9月(共1683 d)历史资料进行判别,满足入型条件的样本数为876 d,消空率为48%;对2013年5—9月(共153 d)历史资料进行判别,入型样本数为80 d,消空率为48%。(3)对2002—2012年5—9月所有入型样本进行回代检验,平均准确率为72.0%(平均TS评分为30.1%);对2013年5—9月所有入型样本进行试预报,平均准确率为63.2%(平均TS评分为28.2%)。     

风云卫星数据服务评价指标的验证分析

简述了目前我国风云卫星数据服务以及气象服务效益评估业务的现状,利用层次分析法提出风云卫星数据服务评价的6个1级指标和19个2级指标。设计了风云卫星数据服务评测量表,并采用分半信度、因子分析等方法对量表作信度和效度分析,结果表明所设计的评测量表具有良好的信度和效度,确定了包括服务平台易用性、服务流程规范性、服务响应及时性、服务内容完整性、服务工作主动性、服务宣传广泛性等6个方面的风云卫星数据服务评价指标。     

天气尺度波列对长江中下游6月梅雨的影响

采用1979—2007年6月NCEP/NCAR2.5°×2.5°逐日再分析资料和中国743站逐日降水资料,利用相关分析、合成分析等方法,分析了天气尺度波列的特征及其对长江中下游6月梅雨的影响。结果表明:当长江中下游6月梅雨较少时,东亚及西太平洋区域存在一个天气尺度波列;该波列的延伸距离较短(从黄河河套地区经过长江中下游至南海、菲律宾海一带),维持时间也很短,且仅仅在500 h Pa以下较强。诊断及个例分析表明,当该波列异常显著时,长江中下游梅雨降水明显减少,而其南部区域降水则增多,说明该波列对预报长江中下游降水具有重要的指示意义。     

从平衡-封闭-静态转向非平衡-开放-动力学研究-热液矿床成因的若干新概念概述

热液金属矿床成因研究过程中,观测与实验始终是密切结合的。上世纪70年代,平衡热力学的实验数据的快速积累,使人们用热力学理论计算可以预测和反演矿石和岩石的成因。但是,没有矿物-水溶液的反应速率数据,又没有与流体力学的结合,搞清楚矿石成因是困难的。七、八十年代,开始研究矿物与水溶液的反应动力学实验。科学家们开始瞄准了从平衡-封闭-静态转向非平衡-开放-动力学研究的这个大方向。1992年我们建立地球化学动力学开放研究实验室。研究高温高压矿物与水反应速率,发现固液的开放体系的自组织现象。实验发现温度影响矿物的各个元素反应速率改变,发现在跨越水临界态时矿物与水反应速率涨落、在近临界的气-液两相不混溶区一些金属进入气相、超临界流体的氧化作用及特别的溶剂性能影响矿物溶解性质。实验证实:临界态区流体与矿石成因有关。水岩相互作用的反应动力学实验温度从低温到550度,揭示矿石的金属来源、迁移、金属与蚀变分带机制。一大批大于300度的矿物与水反应动力学实验在国际界是少有报道的。九十年代,超高压的科学发展,与同步辐射光源的技术进步的结合,使固体地球科学又迈向了地球深内部。我们发展了高温高压流体性质的原位直测(测量850℃水溶液)红外谱,发现深部流体的新性质:气液两相流体的新结构,在临界温度区(300~400℃),水分子氢键网络的破坏受压力影响不大(23MPa~3GPa),同时,出现水的高电导率。研发新仪器为开放-流动-非平衡的反应动力学实验与极端条件下物质性质的直接观测结合,在科学前沿领域开辟了创新道路。     

部分变量误差模型的整体抗差最小二乘估计

部分变量误差模型(partial EIV model)的加权整体最小二乘(weighted total least-squares,WTLS)估计不具备抵御粗差的能力。鉴于粗差可能同时出现在观测值和系数矩阵中,本文在提出部分变量误差模型WTLS估计的两步迭代解法的基础上,运用抗差M估计的等价权方法,发展了一种整体抗差最小二乘(TRLS)估计方法,并采用一致最大功效统计量确定降权因子。针对WTLS估计两步迭代解法的特点,设计了两个不同的降权方案:第1个方案是在估计系数矩阵元素时,不对观测值降权,仅对系数矩阵降权;第2个方案是在估计系数矩阵元素时,既对系数矩阵降权,同时也对观测值降权。通过对模拟2D仿射变换和线性拟合实例进行计算和分析,结果表明第1方案优于第2方案,并且优于基于残差和验后单位权方差的抗差估计和现有的变量误差模型抗差估计。