地-空瞬变电磁法电阻率成像研究与应用

地-空瞬变电磁法（Semi-airborne transient electromagnetic,SATEM）凭借适应能力强、探测深度大、实时性强等特点,适用于湖泊、沼泽、山区等地形复杂地区观测.SATEM接收数据量大、精度要求高,传统成像方法基于地面视电阻率计算未考虑飞行高度,开展SATEM的高精度快速成像研究对实际运用有重要作用.针对SATEM中电性源发射装置感应电流分布,接收装置设置于空中以及两者不在相同平面的特点,本文考虑横向、纵向感应电流差异,飞行高度对视电阻率影响等因素,提出电性源SATEM电阻率成像方法.首先,基于电性源瞬变电磁场感应电流分布特征及各分量空间分布、扩散特性分析,推导出基于均匀半空间磁场与感应电压解析解,然后定义了电性源地-空瞬变电磁法早期、晚期、全期视电阻率;最后根据对感应电流分布分析定义了横向分量和纵向分量的成像深度,并研究了成像深度受飞行高度、发射磁矩、偏移距、偏移角度的影响;对煤炭采空区的实测结果表明,利用本文提出的电阻率成像方法可以更加精确定位目标体位置.

     

超磁致伸缩声波震源的时频特征与实测研究

宽频带、大功率、短余振声波震源是高分辨率探测大尺度复杂结构体的关键技术,为了解决探测距离与其分辨率的矛盾问题,采用超磁致伸缩材料制作的声波发射器作为震源,研制了一种新的声波探测仪器。该震源中心频率为5 kHz,激励电压300~600 V,具有辐射声功率大、余振短（3.5个周期）、频带宽（1~3 kHz）的特性;其最大的优点是每次激发的声波信号一致性好,可以采用单点多次激发,通过多次叠加的数据处理技术提高信噪比,较好地解决了声波探测距离与分辨率的矛盾,并成功地应用于煤矿采场顶煤厚度探测。现场探测结果表明：大功率声波探测系统能够有效提取4~7 m煤岩交界面的反射波信号,借助于小波多分辨分析信号处理方法,有效地提高了采场顶煤厚度探测的精度及其可靠性,为顶煤厚度精确探测提供了一条有效途径。     

CTBTO国际数据中心地震数据处理软件研究进展

全球台网地震震相有效关联与快速精确定位是地震监测和禁核试核查的核心问题之一。禁核试条约组织（CTBTO）筹备委员会临时技术秘书处（PTS）国际数据中心（IDC）自成立以来,一直致力于台网震相自动关联有效率及自动处理公报质量的提高,以减轻分析员的工作量。NET-VISA作为IDC下一代全球震相自动关联和事件检测定位软件,使用概率模型和贝叶斯方法找到解释检测观测结果的最佳事件集,进而从这些观测结果中生成自动处理地震公报。测试结果表明,NET-VISA产生的公报要比IDC目前使用的全球关联方法生成的自动处理公报SEL3更加完整和精准。本文介绍了NET-VISA方法的基本原理和测试进展,最后对CTBT北京国家数据中心参加NET-VISA测试中朝鲜6次核试验的NET-VISA关联结果作了分析对比。     

Characterisation of hydroclimatological trends and variability in the Lake Naivasha basin,Kenya

Recent hydro‐climatological trends and variability characteristics were investigated for the Lake Naivasha basin with the aim of understanding the changes in water balance components and their evolution over the past 50 years. Using a Bayesian change point analysis and modified Mann–Kendall tests, time series of annual mean, maximum, minimum, and seasonal precipitation and flow, as well as annual mean lake volumes, were analysed for the period 1960–2010 to uncover possible abrupt shifts and gradual trends. Double cumulative curve analysis was used to investigate the changes in hydrological response attributable to either human influence or climatic variability. The results indicate a significant decline in lake volumes at a mean rate of 9.35 × 10⁶ m³ year^?1. Most of the river gauging stations showed no evidence of trends in the annual mean and maximum flows as well as seasonal flows. Annual minimum flows, however, showed abrupt shifts and significant (upward/downward) trends at the main outlet stations. Precipitation in the basin showed no evidence of abrupt shifts, but a few stations showed gradual decline. The observed changes in precipitation could not explain the decline in both minimum flows and lake volumes. The findings show no evidence of any impact of climate change for the Lake Naivasha basin over the past 50 years. This implies that other factors, such as changes in land cover and infrastructure development, have been responsible for the observed changes in streamflow and lake volumes. Copyright © 2015 John Wiley & Sons, Ltd.     

Parallelized combined finite-discrete element (FDEM) procedure using multi-GPU with CUDA

This paper focuses on the efficiency of finite discrete element method (FDEM) algorithmic procedures in massive computers and analyzes the time-consuming part of contact detection and interaction computations in the numerical solution. A detailed operable GPU parallel procedure was designed for the element node force calculation, contact detection, and contact interaction with thread allocation and data access based on the CUDA computing. The emphasis is on the parallel optimization of time-consuming contact detection based on load balance and GPU architecture. A CUDA FDEM parallel program was developed with the overall speedup ratio over 53 times after the fracture from the efficiency and fidelity performance test of models of in situ stress, UCS, and BD simulations in Intel i7-7700K CPU and the NVIDIA TITAN Z GPU. The CUDA FDEM parallel computing improves the computational efficiency significantly compared with the CPU-based ones with the same reliability, providing conditions for achieving larger-scale simulations of fracture.     

北斗系统周跳探测与修复方法探究

北斗卫星导航系统是中国自行研制的继美国GPS和俄罗斯GLONASS之后第3个成熟的卫星导航系统。进行GNSS相位数据处理的过程中,由于某种原因,载波相位观测值突然出现整周数的跳变的现象称为周跳。周跳现象的发生将会严重影响后续的数据处理,因此,有必要在解算前完成周跳探测与修复工作。本文就北斗系统周跳探测与修复课题进行了较为全面的探究,先后介绍了高次差法、多普勒法、伪距载波相位组合法、电离层残差法、三频数据组合探测法等一系列周跳探测与修复方法,可分别应用于不同频段接收机、不同观测环境、不同数据处理模式。     

Retrieval of effective leaf area index (LAIe) and leaf area density (LAD) profile at individual tree level using high density multi-return airborne LiDAR

As an important canopy structure indicator, leaf area index (LAI) proved to be of considerable implications for forest ecosystem and ecological studies, and efficient techniques for accurate LAI acquisitions have long been highlighted. Airborne light detection and ranging (LiDAR), often termed as airborne laser scanning (ALS), once was extensively investigated for this task but showed limited performance due to its low sampling density. Now, ALS systems exhibit more competing capacities such as high density and multi-return sampling, and hence, people began to ask the questions like—“can ALS now work better on the task of LAI prediction?” As a re-examination, this study investigated the feasibility of LAI retrievals at the individual tree level based on high density and multi-return ALS, by directly considering the vertical distributions of laser points lying within each tree crown instead of by proposing feature variables such as quantiles involving laser point distribution modes at the plot level. The examination was operated in the case of four tree species (i.e. Picea abies, Pinus sylvestris, Populus tremula and Quercus robur) in a mixed forest, with their LAI-related reference data collected by using static terrestrial laser scanning (TLS). In light of the differences between ALS- and TLS-based LAI characterizations, the methods of voxelization of 3D scattered laser points, effective LAI (LAIe) that does not distinguish branches from canopies and unified cumulative LAI (ucLAI) that is often used to characterize the vertical profiles of crown leaf area densities (LADs) was used; then, the relationships between the ALS- and TLS-derived LAIes were determined, and so did ucLAIs. Tests indicated that the tree-level LAIes for the four tree species can be estimated based on the used airborne LiDAR (R² = 0.07, 0.26, 0.43 and 0.21, respectively) and their ucLAIs can also be derived. Overall, this study has validated the usage of the contemporary high density multi-return airborne LiDARs for LAIe and LAD profile retrievals at the individual tree level, and the contribution are of high potential for advancing forest ecosystem modeling and ecological understanding.     

多光谱遥感影像亮度空间相位一致性特征点检测

提出了一种基于亮度空间和相位一致性理论的多光谱遥感影像特征点检测算法。首先利用参数自适应的灰度变换函数建立影像亮度空间;然后结合相位一致性方法在影像亮度空间进行候选特征点检测,并将候选特征点映射到原始影像上进行非极大值抑制;最后在尺度空间计算特征点的特征尺度值。本文方法有效结合了亮度空间特征检测和相位一致性特征检测的优势,对多光谱遥感影像的辐射变化具有较强的稳健性。试验结果证明,与传统特征点检测算法相比,本文方法在特征重复率和重复特征数量方面都具有明显的优势。     

双频BDS实时精密单点定位算法和精度分析

为了对BDS实时精密单点定位性能进行评估,该文提出了一种适用于BDS系统的实时精密单点定位算法。采用无电离组合模型作为双频实时精密单点定位的数学模型,采用电离层残差法和Melbourne-Wübbena组合实时探测相位周跳,进而单历元实时估计坐标、模糊度等参数,实现了BDS双频实时精密单点定位算法。基于此算法,采用轨道钟差产品和采样间隔为1s的观测数据,模拟实时BDS双频精密单点定位算法,并评估其定位精度。实验结果表明:BDS双频实时定位的平面精度和三维精度均为0.2m左右。     

A synthetic aperture radar sea surface distribution estimation by n-order Bézier curve and its application in ship detection

To dates,most ship detection approaches for single-pol synthetic aperture radar(SAR) imagery try to ensure a constant false-alarm rate(CFAR).A high performance ship detector relies on two key components:an accurate estimation to a sea surface distribution and a fine designed CFAR algorithm.First,a novel nonparametric sea surface distribution estimation method is developed based on n-order Bézier curve.To estimate the sea surface distribution using n-order Bézier curve,an explicit analytical solution is derived based on a least square optimization,and the optimal selection also is presented to two essential parameters,the order n of Bézier curve and the number m of sample points.Next,to validate the ship detection performance of the estimated sea surface distribution,the estimated sea surface distribution by n-order Bézier curve is combined with a cell averaging CFAR(CA-CFAR).To eliminate the possible interfering ship targets in background window,an improved automatic censoring method is applied.Comprehensive experiments prove that in terms of sea surface estimation performance,the proposed method is as good as a traditional nonparametric Parzen window kernel method,and in most cases,outperforms two widely used parametric methods,K and G0 models.In terms of computation speed,a major advantage of the proposed estimation method is the time consuming only depended on the number m of sample points while independent of imagery size,which makes it can achieve a significant speed improvement to the Parzen window kernel method,and in some cases,it is even faster than two parametric methods.In terms of ship detection performance,the experiments show that the ship detector which constructed by the proposed sea surface distribution model and the given CA-CFAR algorithm has wide adaptability to different SAR sensors,resolutions and sea surface homogeneities and obtains a leading performance on the test dataset.