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

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

BDS／GPS 组合 RTK 定位性能分析

针对单系统RTK存在可见卫星数少等问题,文中研究BDS/GPS站间单差的RTK算法模型,该模型采用二次型函数部分最小化及LAMBDA方法联合搜索模糊度。利用该模型分析BDS/GPS组合RTK的定位性能,通过短基线实测数据分析表明:站间单差RTK模型与双差模型是等价的;BDS/GPS组合系统相比于单一系统,明显提高定位的稳健性和精度,改善模糊度固定的成功率。     

从藏南陆-陆碰撞带深部结构构造演化探讨斑岩铜矿的成岩成矿问题

本文以INDEPTH项目对印度大陆与欧亚大陆碰撞带深部成像结果为基础,从构造演化角度探讨藏南陆-陆碰撞带冈底斯斑岩铜矿带的成矿作用问题。深部探测给出的碰撞带深部结构与侯增谦等地质学家提出的深部结构有较大的异同,如何协调起来以深化对藏南陆-陆碰撞条件下成矿作用的认识,这是本文讨论的中心。藏南碰撞带成矿实际上是在新特提斯大洋岩石圈俯冲形成的冈底斯岩浆弧成矿作用的基础上,再经过陆-陆碰撞挤压强烈改造后的再成矿。碰撞带的深部结构构造演化的特点是:(1)新特提斯大洋岩石圈板块向北连续俯冲了约120 Ma,形成的冈底斯陆缘火山岩浆弧带,这导致了陆缘带地壳增厚并含有大量的地幔岩浆流体物质(如南美安第斯成矿带那样);(2)在印度大陆与冈底斯陆缘弧接近碰撞时,在对挤中新特提斯大洋洋壳与大洋岩石圈地幔发生向上挤出与向下拆沉,并使部分洋壳残片和大洋岩石圈物质保存在中上地壳内;(3)两大陆岩石圈碰撞对接后,印度岩石圈地幔加深达70~80 km并沿地壳底部向北推进,并将加厚地壳内大量的成矿物质、钙碱性岩浆,洋壳及新生的下地壳,以及部分地幔物质从地壳底部将其围限起来,成为后期再成矿的物质基础;(4)查明了碰撞带深部壳/幔间产生了一层中间速度层(相当于MASH层),在中上地壳部位出现一层巨大的部分熔融层;(5)在碰撞挤压下冈底斯带内产生多组断裂构造,大型逆冲断裂系与背冲断裂,并引发了含矿岩浆的再活动,并在浮力(下地壳内)和挤压力作用下多次活动上升生成斑岩型铜矿床;(6)成矿后地表遭受过强烈的风化剥蚀作用,使矿床出露地表。     

涉及埃达克岩实验的若干基本概念与问题辨析

根据热力学原理,部分熔融是相变的一种表现形式。当岩石沿p-T轨迹穿过固相线时即可发生熔融。自然界岩石常见的不一致熔融行为可以导致残留相具有不同于源岩的整体成分与矿物组合。残留相与源岩可以属于不同的变质相,岩性也可以不同。"C型埃达克岩"含有含水矿物,它们不是干体系熔融的产物。实验岩石学研究表明,基性成分(SiO_250%)体系较中性、酸性体系在1.0~2.0 GPa压力条件下失水熔融时更易于形成"榴辉岩质"的残留相。源岩除SiO_2以外的其它主要氧化物会影响残留相中各矿物的比例,进而影响熔体的Sr、Y及HREE含量。因此,"C型埃达克岩"高Sr低HREE特征的形成取决于熔融温压条件以及源岩的主量元素、微量元素组成等多重因素。高钾含量(K_2O≈1.0%)的基性、中基性源岩形成的熔体成分与"C型埃达克岩"相比过于富Al或富Na。中等富钾的基性源岩的低程度熔融可以形成高硅的"C型埃达克岩",但无法形成中性的"C型埃达克岩"。     

重力方法在祁连山地区构造特征与岩体解释中的应用

依据地质、航磁等资料,结合实测岩石物性资料,以1∶250 000重力资料为基础对祁连山地区(张掖、刚察地区)进行综合研究。研究地质构造与地球物理特征,划分了大地构造单元,探讨了隐伏与半隐伏岩体分布。该区域经多年的重力工作,勘探面积已近3万多km²,根据重力资料研究布格重力场特征并进行分区; 对地质构造单元重新划分,对主要断裂进行探讨; 圈定隐伏与半隐伏岩体,推断其空间展布特征,并对其解释。结果表明: 祁连山地区地壳密度结构横向分区、纵向分层明显,布格重力异常场的分区特征与地质构造分区特征基本一致。结合2.5D重力剖面反演印证了新的推测结果,为综合地质、地球物理找矿提供参考。     

基于偏最小二乘回归的区域换式风速预报订正技术研究

利用偏最小二乘回归技术,将预报产品作为自变量,相应的自动气象站观测资料作为因变量,对2013-2016年冬季浙江省中尺度区域模式预报近地面风速进行订正和评估。所选956站中多数站点风速订正后有所改善,通过定量分析可知浙江西部地区整体改善效果最好,其中效果明显站点占91.7%;中部地区改善效果明显站点占86.5%;东部沿海地区改善效果略差,明显改善的站点占67%。各地级市整体表现均不错,除舟山地区为49.9%外,其他地区风速改善比例均超过50%。选取2017年1月20日浙江东北地区沿海大风过程分析发现订正后的风速与观测风速更为接近,在定海大岛站点(靠里)中表现尤为明显,订正后的结果具有显著参考价值。     

核心筒部分悬挂结构抗震分析的简化力学模型

提出了核心筒部分悬挂结构体系,在已有该体系振动台试验基础上,建立了该体系的抗震分析简化力学模型及其理论公式。按小震、中震和大震3种不同的地震作用情况,分析了该结构体系力学模型和动力方程的差异。利用有限元分析程序SAP2000,对这种结构体系进行了小震和大震作用下的结构地震反应分析,计算结果与试验结果吻合较好。     

提高机采水平的配套技术及应用

随着濮城油田的持续性高速开发,油田综合含水逐渐上升,油水井井况进一步恶化,作为油田生产主体的机械采油井躺井频繁、管杆泵使用寿命缩短。近几年中原油田采油二厂濮城油田年平均躺井数达560井次之多,导致机械采油井维护费用增加,管杆泵投入增大,严重制约着油田的生产和经营。通过加缓蚀剂在油管内壁形成保护油膜,起到润滑作用;应用特殊抽油杆防腐止磨;应用尼龙扶正器对抽油杆体扶正,避免了钢体与钢体的硬磨;应用旋转井口通过地面人力转动改变油管与抽油杆的偏磨面;配合合理生产参数;各种配套技术的合理应用,使濮城油田躺井减少,避免了频繁作业,增加长寿井,管杆泵投入节约了19.3%。     

pCO2 and carbon fluxes across sea-air interface in the Changjiang Estuary and Hangzhou Bay

Partial pressure of CO2(pCO2) was investigated in the Changjiang(Yangtze River) Estuary,Hangzhou Bay and their adjacent areas during a cruise in August 2004,China.The data show that pCO2 in surface waters of the studied area was higher than that in the atmosphere with only exception of a patch east of Zhoushan Archipelago.The pCO2 varied from 168 to 2 264 μatm,which fell in the low range compared with those of other estuaries in the world.The calculated sea-air CO2 fluxes decreased offshore and varied from -10.0 to 88.1 mmol m-2 d-1 in average of 24.4 ± 16.5 mmol m-2 d-1.Although the area studied was estimated only 2 × 104 km2,it emitted(5.9 ± 4.0) × 103 tons of carbon to the atmosphere every day.The estuaries and their plumes must be further studied for better understanding the role of coastal seas playing in the global oceanic carbon cycle.     

一种机载GNSS高精度定位算法

提出部分模糊度固定的加权电离层模型提高大范围全球卫星导航系统(GNSS)航空定位的精度、可靠性及连续性.该方法的主要思路包括:自适应调整大气扰动等误差影响以实现短基线与长基线两类解算模式之间的灵活切换;施加虚拟电离层观测约束信息,提高基线动态定位的浮点解精度;采用部分模糊度固定方法有效挖掘若干模糊度参数的整周约束.试验表明,提出的方法可提高模糊固定效率与定位精度,克服传统方法有效观测信息利用率不足、定位精度较差、可靠性不高以及连续性较差的问题.实验结果表明,部分模糊度固定算法可在2 min内固定95%以上宽巷模糊度解算与80%以上窄巷模糊度,约20 min后可固定所有模糊度.