多源数据协同定量遥感产品生产系统的领域模型

大数据技术中的一个难点是如何统筹多源异构的数据。在多源数据协同定量遥感产品生产系统中,不同传感器数据的协同使用为生产系统的各个环节都带来了许多的问题,例如异构数据文件的统一表示和调度,繁杂生产流程的统一抽象等。领域驱动设计是一种应对软件核心复杂性的设计方式。领域驱动是指在设计过程中经过不断的迭代,逐渐提炼出一套灵活优雅的领域模型。领域模型专注于分析问题领域本身,发掘重要的业务领域概念,并建立业务领域概念之间的关系。本文在不断实践的基础上提出了一组较成熟的领域模型,该模型用一种统一的方式解决了多源数据协同生产系统中各方面的问题,并显著地降低了系统集成的难度和工作量,且为新数据源的加入预留了灵活性。     

基于重力全张量数据的Parker-Oldenburg算法研究

相比重力数据,重力张量数据通常包含更多的异常信息。本文根据重力数据与重力张量数据的关系,利用位场转化技术,将重力张量数据应用于传统的Parker-Oldenburg密度界面反演算法中。通过模型试验,证明了在网格间距较大或者数据存在一定噪音时,使用本文算法进行反演能得到更好的效果。实验结果说明利用重力张量数据可以有效地提高密度界面反演的分辨率。     

Analysis of wave groups by wave envelope-phase and the Hilbert Huang transform methods

The wave groups are studied by both conventional wave analysis methods and by the non-stationary Hilbert Huang Transform (HHT) method. Full-scale wave records containing abnormal waves are used. Instantaneous quantities, such as envelope, phase and frequency, are adopted to study the wave grouping. A refined definition of wave group is proposed considering that the wave process is simultaneously amplitude and frequency modulated. The validation of the proposed definition is conducted by analysis of numerical simulation data. Group parameters are proposed based on the time-frequency distribution of energy. An attempt is made to find the relationship between the characteristics of abnormal waves and the group characteristics.     

A Method for Automatically Selecting Aids to Navigation Based on their Spatial Influence Domains

To address the limitations of manually selecting aids to navigation (AtNs) on charts, a method for automatically selecting AtNs based on their spatial influence domains (SIDs) is proposed. First, the associations between the spatial attributes of an AtN are analyzed. Second, an SID of the AtN is defined, and a model of the SID is constructed based on the associations between the spatial attributes. Third, the importance of the location of the AtN is weighted based on the SID model. Fourth, an algorithm to automatically select AtNs based on the maximum coverage of the SIDS of preselected AtNs is developed. Finally, several AtNs are selected automatically using the algorithm. The experimental results demonstrate that (1) the proposed method can automatically select AtNs and the results comply with the requirements; (2) the automatic selection can eliminate the human-induced errors or the inconsistent results of manual selections from different operators; and (3) the efficiency of the proposed method is higher than that of current manual methods.     

高阶边界元方法求解规则波在三维局部渗透海床上的传播

基于线性势流理论,利用高阶边界元法研究了规则波在三维局部渗透海床上的传播。根据Darcy渗透定律推导出渗透海床的控制方程,利用渗透海床顶部和海底处法向速度和压强连续条件得到渗透海床顶部满足的边界条件。根据绕射理论,利用满足自由水面条件的格林函数建立了求解渗透海床绕射势的边界积分方程,采用高阶边界元方法求解边界积分方程进而得到自由水面的绕射势和波浪在局部渗透海床上传播过程中幅值的变化情况。通过与已发表的波浪对圆柱形暗礁的时域全绕射结果对比,证明了本文建立的频域方法计算波幅的正确性和有效性。利用这一模型研究了三维矩形渗透海床区域上波浪的传播特性,并分析了入射波波长、海床渗透特性系数等参数对波浪传播的影响。     

Seismic attenuation estimation using a complete ensemble empirical mode decomposition-based method

The frequency attenuation gradient method can provide important information for hydrocarbon detection. In this paper, a method using Complete Ensemble Empirical Mode Decomposition (CEEMD), Hilbert transform and the least-squares curve-fitting is proposed for seismic attenuation estimation as an effective frequency attenuation gradient estimation approach. We first use CEEMD to obtain the different Intrinsic Mode Functions (IMFs), which have a narrow band and can enhance the physical meaning of instantaneous attributes trace by trace. The time-frequency spectrum, which is computed using a Hilbert transform of each IMF, is represented as a spectrum with a single-peak that has narrow side lobes, which is conducive to frequency attenuation gradient estimation. Second, for each time sample, the frequency-amplitude spectrum of each IMF trace is extracted from the time-frequency spectrum to conduct the attenuation gradient computation. Then, the logarithm operation is performed for each IMF trace. Due to the very narrow bands of some IMFs in some seismic traces, a variable frequency window is adopted along the IMF trace according to the local data characteristics. Finally, the attenuation gradient for each IMF in a seismic trace can be computed using least-squares fitting. A different IMF reflects a seismic trace with a different spatiotemporal scale and can highlight different geologic and stratigraphic information. The correlation weighted average operation is used to highlight some useful details in seismic trace and obtains the attenuation gradient for each seismic trace. Field data examples demonstrate our method and its effectiveness. The proposed method can stably estimate the frequency attenuation gradient.     

应用两次展开法求解系泊柱体慢漂运动方程

在深海中系泊的海洋平台,如Spar平台,水下部分为带有系泊的圆柱结构,其水平方向运动响应往往具有较低的自振频率,容易在低频波浪力(源于非线性的差频效应)作用下发生共振响应,使结构发生大幅水平慢漂。当浮体的瞬时位置大幅偏离初始位置时,基于初始平衡位置的摄动展开法会存在较大误差。针对这一问题,采用两次展开方法,对大幅慢漂运动开展时域模拟研究。对双色波作用下自由漂浮圆柱的大幅运动响应问题进行数值分析,并与采用基于初始平衡位置的摄动展开法的计算结果进行了对比。结果表明,采用新的两次展开法可以计算出波浪遭遇频率的变化和波浪漂移阻尼,而这无法从基于初始平衡位置的常规摄动展开法中得到,体现了两次展开法在分析大幅慢漂问题上的优势。     

Metallogeny of the Baiyangping Lead‐Zinc Polymetallic Ore Concentration Area,Northern Lanping Basin of Yunnan Province,China

The Lanping Basin in the Nujiang‐Lancangjiang‐Jinshajiang (the Sanjiang) area of northeastern margin of the Tibetan Plateau is an important part of eastern Tethyan metallogenic domain. This basin hosts a number of large unique sediment‐hosted Pb‐Zn polymetallic deposits or ore districts, such as the Baiyangping ore concentration area which is one of the representative ore district. The Baiyangping ore concentration area can be divided into the east and west ore belts, which were formed in a folded tectogene of the India‐Asia continental collisional setting and was controlled by a large reverse fault. Field observations reveal that the Mesozoic and Cenozoic sedimentary strata were outcropped in the mining area, and that the orebodies are obviously controlled by faults and hosted in sandstone and carbonate rocks. However, the ore‐forming elements in the east ore belt are mainly Pb‐Zn‐Sr‐Ag, while Pb‐Zn‐Ag‐Cu‐Co elements are dominant in the west ore belt. Comparative analysis of the C‐O‐Sr‐S‐Pb isotopic compositions suggest that both ore belts had a homogeneous carbon source, and the carbon in hydrothermal calcite is derived from the dissolution of carbonate rock strata; the ore‐forming fluids were originated from formation water and precipitate water, which belonged to basin brine fluid system; sulfur was from organic thermal chemical sulfate reduction and biological sulfate reduction; the metal mineralization material was from sedimentary strata and basement, but the difference of the material source of the basement and the strata and the superimposed mineralization of the west ore belt resulted in the difference of metallogenic elements between the eastern and western metallogenic belts. The Pb‐Zn mineralization age of both ore belts was contemporary and formed in the same metallogenetic event. Both thrust formed at the same time and occurred at the Early Oligocene, which is consistent with the age constrained by field geological relationship.     

基于小波分析与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法研究地下水动态的适宜性和可靠性,也为济南市的保泉提供了参考依据。      

Landslide susceptibility zonation method based on C5.0 decision tree and K-means cluster algorithms to improve the efficiency of risk management

Machine learning algorithms are an important measure with which to perform landslide susceptibility assessments,but most studies use GIS-based classification methods to conduct susceptibility zonation.This study presents a machine learning approach based on the C5.0 decision tree(DT)model and the K-means cluster algorithm to produce a regional landslide susceptibility map.Yanchang County,a typical landslide-prone area located in northwestern China,was taken as the area of interest to introduce the proposed application procedure.A landslide inventory containing 82 landslides was prepared and subse-quently randomly partitioned into two subsets:training data(70％landslide pixels)and validation data(30％landslide pixels).Fourteen landslide influencing factors were considered in the input dataset and were used to calculate the landslide occurrence probability based on the C5.0 decision tree model.Susceptibility zonation was implemented according to the cut-off values calculated by the K-means clus-ter algorithm.The validation results of the model performance analysis showed that the AUC(area under the receiver operating characteristic(ROC)curve)of the proposed model was the highest,reaching 0.88,compared with traditional models(support vector machine(SVM)=0.85,Bayesian network(BN)=0.81,frequency ratio(FR)=0.75,weight of evidence(WOE)=0.76).The landslide frequency ratio and fre-quency density of the high susceptibility zones were 6.76/km2 and 0.88/km2,respectively,which were much higher than those of the low susceptibility zones.The top 20％interval of landslide occurrence probability contained 89％of the historical landslides but only accounted for 10.3％of the total area.Our results indicate that the distribution of high susceptibility zones was more focused without contain-ing more"stable"pixels.Therefore,the obtained susceptibility map is suitable for application to landslide risk management practices.