安徽凤阳山区石英岩矿成矿规律分析

凤阳山区石英岩矿床成因类型主要有沉积变质型和沉积型两种,均属层控性矿床。其中沉积变质型石英岩矿是本区主要类型矿床,赋存构造体系为老青山单斜,控矿层位为凤阳群白云山组地层,矿石类型主要为石英岩,其次是石英砂岩、含白云母石英岩等;沉积型石英岩矿赋存构造体系为雷家户-宋集单斜,控矿层位为八公山群伍山组,矿石类型主要为石英砂岩。区内发育的褶皱构造为石英岩矿体提供了广阔的赋存空间,而断裂构造的发育却对矿体分布的连续性造成了不同程度的破坏,对石英岩矿石质量也有一定的影响。     

南极南设得兰群岛海域地质新进展:“海洋六号”南极科考初步成果

2016—2017年跨年度执行的中国第33次南极科考“海洋六号”航次, 首次获取了南极半岛南设得兰群岛北部海域、布兰斯菲尔德海峡以及欺骗岛三个区域的17 850 km2范围内多波束全覆盖高精度水深数据, 地形分辨率高于目前公开的国际数据。多道地震结果显示调查区多处存在天然气水合物地球物理标志BSR。地热测量热流值表明南设得兰群岛以北热流值与全球海洋热流值数据吻合, 不存在明显异常, 群岛以南热流值数据偏高。沉积物样品为典型的冰海沉积物, 在南设得兰群岛与南设得兰海沟之间, 沉积物以粉砂-砂为主, 夹杂砾石, 随水深增加, 沉积物粒径变细; 而在布兰斯菲尔德海区表层沉积物分布具有不对称性, 海槽南部陆架-上陆坡以及海槽东北和相邻陆架, 沉积物以砾石、砂质为主; 海槽和下陆坡为粉砂质硅质软泥。沉积物顶空气样品现场测试结果表明, 表层沉积物50 cm左右深度内的甲烷含量没有明显的差异性, 但在岛坡深水区两个测站甲烷含量在270 cm 左右深度时都存在局部高异常, 且两者总甲烷含量的变化趋势有一定的相似性, 可能与水合物分解迁移有关。     

基于小波分析与Mann-Kendall法的岩溶大泉动态研究

研究地下水动态是认识地下水资源的有效手段。根据1956－2013年济南岩溶泉域大气降水及地下水水位动态监测资料,采用小波分析法、Mann-Kendall趋势检验、突变检验法研究了58个水文年泉水位对大气降水的响应,可以看出:（1）大气降水和泉水位呈现出多尺度的变化特征,长时间尺度上两者的变化周期基本相同,变化周期为16年和12年,说明大气降水对泉水位有直接影响;（2）在1956—2013年,济南泉域地下水水位具有0.65 m·（10a）-1的年际显著下降趋势,但降水具有12.65 mm·（10a）-1的不显著上升趋势,说明在人为因素影响下泉水动态的影响因素的权重发生了变化;（3）大气降水在1999年发生突变,1999年之后年降水为增加趋势;而地下水水位突变年份为1967年,1967年以后水位持续降低,2004年以后水位快速上升,泉水位未来趋势应与降水保持一致,呈上升趋势,说明大气降水并非泉水动态的唯一影响因素;（4）通过建立不同时段的多元回归模型,表明近58年来地下水水位的主要影响因素由大气降水到人工开采之间的转换,同时验证了小波分析和Mann-Kendall法研究地下水动态的适宜性和可靠性,也为济南市的保泉提供了参考依据。      

新时代地热行业发展特点及几点思考

当前,中国特色社会主义进入了新时代,社会各项建设取得重大成就,地热产业也发展壮大。在大力倡导生态文明建设、推进绿色发展、提升环境保护、推进资源的节约和循环利用的社会发展形势下,一些重要规划和政策的实施、一批重点区域和重大项目的建设、大量地热能开发利用新技术的推广,又给地热产业带来了新的发展契机。不同形式、不同利用方式的地热能都将迎来更广阔的市场、更高的产业增长态势。为顺应新时代行业发展,地热产业应抓住机遇,迎接挑战,创新突破,成就未来。本文将从上位决策、项目建设、技术创新、管理等方面,分析当下地热行业发展的形势、特点和存在的问题,提出相关建议,以共同探讨,相互借鉴。     

A modified mix-differenced approach for estimating multi-GNSS real-time satellite clock offsets

A modified mixed-differenced approach for estimating multi-GNSS real-time clock offsets is presented. This approach, as compared to the earlier presented mixed-differenced approach which uses epoch-differenced and undifferenced observations, further adds a satellite-differenced process. The proposed approach, based on real-time orbit products and a mix of epoch-differenced and satellite-differenced observations to estimate only satellite clock offsets and tropospheric zenith wet delays, has fewer estimated parameters than other approaches, and thus its implementing procedure is efficient and can be performed and extended easily. To obtain high accuracy, the approach involves three steps. First, the high-accuracy tropospheric zenith wet delay of each station is estimated using mixed-differenced carrier phase observations. Second, satellite clock offset changes between adjacent epochs are estimated using also mixed-differenced carrier phase observations. Third, the satellite clock offsets at the initial epoch are estimated using satellite-differenced pseudorange observations. Finally, the initial epoch clock results and clock offset changes are concatenated to obtain the clock results of the current epoch. To validate the real-time satellite clock results, multi-GNSS post-processing clock products from IGS ACs were selected for comparison. From the comparison, the standard deviations of the GPS, GLONASS, BeiDou and Galileo systems clock results are approximately 0.1–0.4 ns, except for the BeiDou GEO satellites. The root mean squares are about 0.4–2.3 ns, which are similar to those of other international real-time products. When the clock estimates were assessed based on a pseudo-kinematic PPP procedure, the positioning accuracies in the East, North and Up components reach 5.6, 5.5 and 7.6 cm, respectively, which meet the centimeter level and are comparable to the application of other products.     

Estimating satellite phase fractional cycle biases based on Kalman filter

Phase fractional cycle biases (FCBs) originating from satellites and receivers destroy the integer nature of PPP carrier phase ambiguities. To achieve integer ambiguity resolution of PPP, FCBs of satellites are required. In former work, least squares methods are commonly adopted to isolate FCBs from a network of reference stations. However, it can be extremely time consuming concerning the large number of observations from hundreds of stations and thousands of epochs. In addition, iterations are required to deal with the one-cycle inconsistency among FCB measurements. We propose to estimate the FCB based on a Kalman filter. The large number of observations are handled epoch by epoch, which significantly reduces the dimension of the involved matrix and accelerates the computation. In addition, it is also suitable for real-time applications. As for the one-cycle inconsistency, a pre-elimination method is developed to avoid iterations and posterior adjustments. A globally distributed network consisting of about 200 IGS stations is selected to determine the GPS satellite FCBs. Observations recorded from DoY 52 to 61 in 2016 are processed to verify the proposed approach. The RMS of wide lane (WL) posterior residuals is 0.09 cycles while that of the narrow lane (NL) is about 0.05 cycles, which indicates a good internal accuracy. The estimated WL FCBs also have a good consistency with existing WL FCB products (e.g., CNES-GRG, WHU-SGG). The RMS of differences with respect to GRG and SGG products are 0.03 and 0.05 cycles. For satellite NL FCB estimates, 97.9% of the differences with respect to SGG products are within ±?0.1 cycles. The RMS of the difference is 0.05 cycles. These results prove the efficiency of the proposed approach.     

Mapping hourly dynamics of urban population using trajectories reconstructed from mobile phone records

Understanding the spatiotemporal dynamics of urban population is crucial for addressing a wide range of urban planning and management issues. Aggregated geospatial big data have been widely used to quantitatively estimate population distribution at fine spatial scales over a given time period. However, it is still a challenge to estimate population density at a fine temporal resolution over a large geographical space, mainly due to the temporal asynchrony of population movement and the challenges to acquiring a complete individual movement record. In this article, we propose a method to estimate hourly population density by examining the time‐series individual trajectories, which were reconstructed from call detail records using BP neural networks. We first used BP neural networks to predict the positions of mobile phone users at an hourly interval and then estimated the hourly population density using log‐linear regression at the cell tower level. The estimated population density is linearly correlated with population census data at the sub‐district level. Trajectory clustering results show five distinct diurnal dynamic patterns of population movement in the study area, revealing spatially explicit characteristics of the diurnal commuting flows, though the driving forces of the flows need further investigation.     

针对星敏感器非线性畸变的非参数化在轨检校方法

为了对星敏感器的非线性畸变进行在轨检校,本文提出了一种非参数化的检校方法,该方法以B样条格网的形式表达非线性畸变,不需要预先估计一个参数化的畸变模型,从而可以处理任何形式的畸变。针对星敏感器不同模式畸变的仿真实验表明,本文提出的方法有效、可靠并且具有鲁棒性。     

塔里木盆地周缘新元古代冰期事件与古老烃源岩发育

正四川盆地震旦系含油气系统已被证实为我国新元古界第一个具有商业价值的原生型气藏,而与扬子板块在同时期具有亲缘性的塔里木盆地新元古界油气资源一直未取得突破。塔里木盆地新元古界埋藏较深,钻遇南华系-震旦系的井位非常少,盆地深部层系的成像品质较差,同时新元古界本身缺乏可进行有效对比的古生物化石。诸多因素造成塔里木盆地新元古界划分与对比关系不清,地层格架模糊,直接导致后续深入研究无法进行。然而,