联合单站地基GNSS和空基掩星近实时反演电离层电子密度

采用非组合精密单点定位技术估算VTEC,并选取一系列与测站共址的COSMIC掩星电离层产品共同作为约束条件,利用函数模型近实时地探测GNSS测站上空电离层电子密度剖面特性。为了验证新方法的可行性,选取东亚11个IGS跟踪站2016年4个季节中3月、6月、9月和12月最后一周的观测数据及与测站共址的COSMIC掩星ionPrf产品进行实验。将每一个反演结果与同一时刻COSMIC掩星探测的电离层电子密度（作为真值）进行比较。结果表明,新方法实时反演的电子密度剖面与掩星电离层产品相比具有很好的一致性,所获得的f0F2与hmF2参数精度皆优于IRI2016模型。     

信号易遮挡地区GPS/BDS双频单历元短基线解算精度分析

信号易遮挡地区GPS、BDS单系统可见卫星个数少,无法达到相对定位的最低要求。讨论了GPS/BDS组合单历元基线解算模型,在多种模拟环境下将GPS/BDS联合解算结果与GPS和BDS单系统在可见卫星数、模糊度固定成功率、定位精度等方面进行对比分析。结果表明,相对于单系统单历元基线解算,GPS/BDS组合大大增加了可见卫星数,提高了模糊度固定成功率,并在极端的观测条件下有效地改善了定位精度,使信号易遮挡地区双频单历元基线解算精度可以达到cm级甚至mm级。     

采煤塌陷地复垦项目PPP融资模式应用研究

探索PPP融资的运作模式与双方合作机制,促进PPP融资模式在采煤塌陷地复垦项目中的应用。研究基于风险-收益均衡原则,从私人资本与政府的风险偏好特点出发,建立了项目总成本最小的最优风险分担模型;在风险分担系数已知的情况下,建立了Shapley值模型对双方收益进行合理分配。结果表明:私人投资主体的综合风险系数略大,双方基本上均担风险;私人投资主体承担绝大部分的项目资金投入,所获收益略大于政府。因此引入PPP模式进行采煤塌陷地复垦,对于政府而言,可以解决资金压力,分散不擅长的项目与金融、市场风险;对于私人投资主体而言,可以避免过多的政策与不确定性风险,获得较高的投资收益。该研究为采煤塌陷地复垦项目PPP融资谈判提供了理论依据,对加快塌陷区复垦进程具有非常重要的现实意义。     

基于间断自适应MRF的相位解缠算法

合成孔径干涉雷达（InSAR）技术已被广泛应用于地面高程信息的获取。相位解缠作为其关键步骤直接关系到地表高程的提取精度。然而,相位解缠受相位连续性假设的限制无法识别间断相位,影响在陡峭地形区域的高程反演精度。针对上述问题,本文提出了一种基于间断自适应马尔可夫随机场模型（DA-MRF）的相位解缠算法处理间断相位。首先构建DA-MRF模型的能量函数能够自适应地在连续相位处保持像素连接,同时在间断相位处阻断像素连接;然后再利用图割算法优化能量函数,完成相位解缠。仿真和实测数据试验证明,本文方法与常用的相位解缠算法相比,具有更强的间断相位保持能力,更适用于地形变化剧烈场景中的地形重建。     

陆架海岸台风沉积记录及信息提取

长时间尺度风暴强度–频率关系与气候变化相关联,而器测记录和历史记载难以提供充分的信息,因此从沉积记录中提取风暴信息成为一个前沿科学问题。在应用上,这项研究可为海岸带城市群应对未来气候和海面变化提供决策依据。本文回顾了台风沉积记录研究进展,显示陆架泥质沉积、海滩及海岸沙丘、潮滩、潟湖、巨砾是台风事件记录的良好载体,可通过层序形态和物质特性分析而识别。同时,还需进一步完善分析方法,以区分台风、冬季风暴、河流洪水和海啸等不同类型的极端事件沉积。在台风强度信息提取方面,陆架泥质沉积所含贝壳–粗颗粒沉积物可作为海底再悬浮强度的指标,但需更多实测数据的率定;海滩及海岸沙丘顶部的台风沉积分布高程指示了台风激浪流的上冲高度,而台风巨砾的重量可以与近岸波浪的波高建立联系。以上数据经过换算后可以得出台风强度的信息,虽然这些间接的沉积学信息还不足以建立风暴强度–频率关系,但有助于台风强度大数据的建立。潮滩、潟湖沉积连续性好,可构成台风事件的时间序列,然而关于台风强度却是多解的,台风最大风力、持续时间、移动路径、登陆地点的不同组合可能产生同样的事件沉积。我们建议,应发展台风信息提取的新方法,来解决这个问题。进行现代过程模拟,根据已知的台风事件资料构建沉积物输运堆积模型,使之能够复演事件沉积的特征;进行多个地点事件沉积的反演模拟,在此情形下,即便每个站位的结果是多解的,但针对多个站位上求取其解的交集之后,多解性将下降,这种模拟方法可称之为“解空间收缩法”;采用大数据融合方式,将其他来源的台风强度数据纳入模拟体系,可进一步降低风暴信息提取的不确定性。动力过程模拟与大数据融合方法的建立,有助于获得与沉积记录同样时间尺度的台风强度–频率关系曲线,进而分析台风动态与气候变化的关系。     

有界赤道大洋波包解及其年际年代际变率

Linearized shallow water perturbation equations with approximation in an equatorial β plane are used to obtain the analytical solution of wave packet anomalies in the upper bounded equatorial ocean. The main results are as follows. The wave packet is a superposition of eastward travelling Kelvin waves and westward travelling Rossby waves with the slowest speed, and satisfies the boundary conditions of eastern and western coasts, respectively.The decay coefficient of this solution to the north and south sides of the equator is inversely proportional only to the phase velocity of Kelvin waves in the upper water. The oscillation frequency of the wave packet, which is also the natural frequency of the ocean, is proportional to its mode number and the phase velocity of Kelvin waves and is inversely proportional to the length of the equatorial ocean in the east-west direction. The flow anomalies of the wave packet of Mode 1 most of the time appear as zonal flows with the same direction. They reach the maximum at the center of the equatorial ocean and decay rapidly away from the equator, manifested as equatorially trapped waves. The flow anomalies of the wave packet of Mode 2 appear as the zonal flows with the same direction most of the time in half of the ocean, and are always 0 at the center of the entire ocean which indicates stagnation, while decaying away from the equator with the same speed as that of Mode 1. The spatial structure and oscillation period of the wave packet solution of Mode 1 and Mode 2 are consistent with the changing periods of the surface spatial field and time coefficient of the first and second modes of complex empirical orthogonal function(EOF)analysis of flow anomalies in the actual equatorial ocean. This indicates that the solution does exist in the real ocean, and that El Ni?o-Southern Oscillation(ENSO) and Indian Ocean dipole(IOD) are both related to Mode 2.After considering the Indonesian throughflow, we can obtain the length of bounded equatorial ocean by taking the sum of that of the tropical Indian Ocean and the tropical Pacific Ocean, thus this wave packet can also explain the decadal variability(about 20 a) of the equatorial Pacific and Indian Oceans.     

三区复合凝析气藏多分支水平井压力动态特征

钻井技术的发展使多分支水平井应用于海上凝析气藏开发。根据凝析气藏开采过程中的相变特征,将气藏划分为3个区:Ⅰ区为油气同流区,Ⅱ区为气体流动、油不流动的过渡区,Ⅲ区为纯气区。建立三区复合凝析气藏多分支水平井的解析模型,采用点源函数、Fourier有限余弦积分变换、Green函数等方法求解模型,利用MATLAB软件绘制井底压力特征曲线,分析凝析气藏多分支水平井的渗流规律,对相关参数进行敏感性分析。结果表明:凝析气藏多分支水平井的渗流过程划分为9个阶段;应力敏感效应使第二径向流阶段后的压力特征曲线逐渐上翘;分支井长度增加,分支井线性流持续时间变长;分支井间距减小,两分支井间的压力波相互干扰,导致第二径向流水平段缩短;分支井数量减少使早期的压力特征曲线向上偏移;相邻两区流度比越接近1,曲线后期上移越明显;区域半径影响径向流阶段持续时间。该模型能更准确地模拟凝析气藏多分支水平井的生产动态,为凝析气藏开发提供指导。     

GNSS精密单点定位与惯性导航紧组合理论模型

针对精密单点定位(Precise Point Positioning,PPP)动态定位精度低、收敛速度慢等问题,本文采用PPP/INS紧组合系统来达到改善PPP动态定位性能的目的。本文对PPP/INS紧组合的观测方程、误差补偿模型、参数估计模型等进行详细推导。通过车载实验采集的卫星观测数据和不同等级的惯性数据,对动态PPP及PPP/INS紧组合的定位定姿性能进行分析,评估不同等级惯性传感器对PPP/INS紧组合定位精度和收敛速度的影响。实验结果表明:PPP/INS紧组合在北—东—高方向的位置RMS相对于PPP分别改善了70.2%、29.1%和16.8%,达到4.8 cm、12.3 cm和7.4 cm。在卫星跟踪条件良好时,惯性传感器性能对PPP/INS紧组合定位精度影响不大;而在卫星观测条件不足时,惯性传感器性能对PPP/INS紧组合定位精度影响明显。此外,仿真和恶劣条件下的数据结果表明,PPP/INS初始定位精度与收敛速度随惯性传感器性能提高而改善明显。     

熊俐治疗女性生殖道解脲支原体感染经验

介绍熊俐教授治疗女性生殖道解脲支原体（Uu）感染的临床经验。熊教授认为支原体感染为正气不固,湿毒之邪内侵,治当以除湿解毒,益气扶正,自拟除湿解毒扶正方联合阿奇霉素以治之,取得较好的临床疗效。附验案1则,以资佐证。