Generation and Development of Structures and Their Controls on Ore Mineralization in the Lehong Zinc-lead Deposit,Northeastern YunnanEI北大核心CSCD

The Lehong zinc-lead deposit is one of the new-found large Zn-Pb deposits in the Sichuan- Yunnan-Guizhou Zn-Pb poly-metallic mineralization area. The ore-bodies of the Lehong deposit are strictly controlled by structures, and the fault tectonites can be divided into four groups which include morbruk rock, cataclastic rock and tectonic breccia. The mechanical properties of the faults suggest that there are five tectonic systems formed in the Caledonian-early Hercynian, Indosinian-early Yanshanian, mid-Yanshanian, late Yanshanian and Himalayan periods, respectively. The Lehong fault and the Qiaojia-Lianfeng fault are two important ore-controlling structures which might have acted as channel ways of the ore fluids. The ores are hosted in the subordinary fracture zone, joint fracture zone and interlayer fault zones in trailing folds of the Lehong fault and the Baobaoshang syncline and Jinjiagou anticline. Hence, the Lehong deposit is a tectonic -controlled ore deposit. These research do not only lay a foundation for studying deposit genesis and ore-finding direction in the Lehong deposit, but also provides important enlightenments for the ore-finding prognosis of lead-zinc deposits in northwestern Guizhou ore concentration district. © 2018, Science Press. All right reserved.     

Geometric-constrained multi-view image matching method based on semi-global optimization

Targeting at a reliable image matching of multiple remote sensing images for the generation of digital surface models, this paper presents a geometric-constrained multi-view image matching method, based on an energy minimization framework. By employing a geometrical constraint, the cost value of the energy function was calculated from multiple images, and the cost value was aggregated in an image space using a semi-global optimization approach. A homography transform parameter calculation method is proposed for fast calculation of projection pixel on each image when calculating cost values. It is based on the known interior orientation parameters, exterior orientation parameters, and a given elevation value. For an efficient and reliable processing of multiple remote sensing images, the proposed matching method was performed via a coarse-to-fine strategy through image pyramid. Three sets of airborne remote sensing images were used to evaluate the performance of the proposed method. Results reveal that the multi-view image matching can improve matching reliability. Moreover, the experimental results show that the proposed method performs better than traditional methods.     

Position-domain integrity risk-based ambiguity validation for the integer bootstrap estimator

Integrity monitoring for ambiguity resolution is of significance for utilizing the high-precision carrier phase differential positioning for safety–critical navigational applications. The integer bootstrap estimator can provide an analytical probability density function, which enables the precise evaluation of the integrity risk for ambiguity validation. In order to monitor the effect of unknown ambiguity bias on the integer bootstrap estimator, the position-domain integrity risk of the integer bootstrapped baseline is evaluated under the complete failure modes by using the worst-case protection principle. Furthermore, a partial ambiguity resolution method is developed in order to satisfy the predefined integrity risk requirement. Static and kinematic experiments are carried out to test the proposed method by comparing with the traditional ratio test method and the protection level-based method. The static experimental result has shown that the proposed method can achieve a significant global availability improvement by 51% at most. The kinematic result reveals that the proposed method obtains the best balance between the positioning accuracy and the continuity performance.     

A modified three-dimensional ionospheric tomography algorithm with side rays

The three-dimensional ionospheric tomography (3DCIT) algorithm based on Global Navigation Satellite System (GNSS) observations have been developed into an effective tool for ionospheric monitoring in recent years. However, because the rays that come into or come out from the side of the inversion region cannot be used, the distribution of the rays in the edge and bottom part of the inversion region is scarce and the electron density cannot be effectively improved in the inversion process. We present a three-dimensional tomography algorithm with side rays (3DCIT-SR) applying the side rays to the inversion. The partial slant total electron content (STEC) of side rays in the inversion region is obtained based on the NeQuick2 model and GNSS-STEC. The simulation experiment results show that the algorithm can effectively improve the distribution of GNSS rays in the inversion region. Meanwhile, the iteration accuracy has also been significantly improved. After the same number of iterations, the iterative results of 3DCIT-SR are closer to the truth than 3DCIT, in particular, the inversion of the edge regions is improved noticeably. The GNSS data of the International GNSS Service (IGS) stations in Europe are used to perform real data experiments, and the inversion results show that the electron density profiles of 3DCIT-SR are closer to the ionosonde measurements. The accuracy improvement of 3DCIT-SR is up to 56.3% while the improvement is more obvious during the magnetic storm compared to the case of a calm ionospheric state .     

A new GPS/BDS tropospheric delay resolution approach for monitoring deformation in super high-rise buildings

A new approach for deformation monitoring of super high-rise building using GPS/BDS technology is proposed for the case when prior coordinates are known and the baseline is short but has a large height difference. The approach is based on the ambiguity function method (AFM). Considering that the double-differenced (DD) troposphere delay residual error cannot be ignored, the relative zenith tropospheric delay (RZTD) parameter is introduced into the original AFM equation. Thus, the RZTD and 3D coordinate parameters are together obtained through the modified AFM (MAFM). Due to the low computational efficiency of conventional AFM, an improved particle swarm optimization (IPSO) algorithm is used to search the four optimal parameters X/Y/Z/RZTD and replaces the grid search method. In this study, GPS/BDS deformation monitoring data for buildings with approximately 290 m height difference were used to verify the feasibility of the proposed MAFM. Numerical results show a single-epoch average computation time of approximately 0.3 s, which meets the requirements of near-real-time dynamic monitoring. The average accuracy of the GPS single-epoch RZTD solution is better than 1 cm, the combined GPS/BDS MAFM performance outperforms the GPS-only system, and using multi-epoch observations can further improve the accuracy of the RZTD solution. After RZTD correction, GPS/BDS monitoring precision can be improved, particularly the height dimension, whose precision is improved by approximately 6 cm.     

Analysis of estimated satellite clock biases and their effects on precise point positioning

This study provides a systemic analysis to identify the biases in estimated satellite clocks and illustrates their effects in precise point positioning (PPP). First, the precise satellite clock estimation method considering pseudorange and carrier phase hardware delays is derived. Two methods for satellite clock estimation are compared, and their equivalency is discussed. The results show that apart from the well-known constant code hardware biases, the time-variant phase hardware biases are also absorbed by the estimated clocks. Also, the satellite clocks contain biases caused by modeling errors. To analyze the effects of these biases, they are grouped into initial clock biases (ICBs) and time-dependent biases (TDBs). Then, a detailed analysis of the impact of the biases on PPP-based troposphere and coordinate estimates is conducted. The experimental analysis demonstrates that TDBs affect positioning and tropospheric estimates, and their impacts are more significant in the static mode. The ICBs affect coordinate accuracy, zenith total delay mean bias, and its standard deviations only at the millimeter level for kinematic and static PPP, which is negligible. However, the ICBs affect the convergence period for both static and real kinematic PPP, and the magnitude of their impact largely depends on data quality. Note that satellites clocks are generally estimated with the P1/P2 and L1/L2 ionospheric-free combinations, and that hardware-specific parts of ICBs and TDBs cancel if users employ the same type of observables as the clock providers. Otherwise, the effects of biases cannot be ignored, especially for triple-frequency applications. Also, modeling-specific parts of ICBs and TDBs are significant in real-time clocks, which also affect user applications. Our conclusion is applicable for understanding the effects of these biases.     

Integrated compression of vehicle spatio‐temporal trajectories under the road stroke network constraint

With fast growth of all kinds of trajectory datasets, how to effectively manage the trajectory data of moving objects has received a lot of attention. This study proposes a spatio‐temporal data integrated compression method of vehicle trajectories based on stroke paths coding compression under the road stroke network constraint. The road stroke network is first constructed according to the principle of continuous coherence in Gestalt psychology, and then two types of Huffman tree—a road strokes Huffman tree and a stroke paths Huffman tree—are built, based respectively on the importance function of road strokes and vehicle visiting frequency of stroke paths. After the vehicle trajectories are map matched to the spatial paths in the road network, the Huffman codes of the road strokes and stroke paths are used to compress the trajectory spatial paths. An opening window algorithm is used to simplify the trajectory temporal data depicted on a time–distance polyline by setting the maximum allowable speed difference as the threshold. Through analysis of the relative spatio‐temporal relationship between the preceding and latter feature tracking points, the spatio‐temporal data of the feature tracking points are all converted to binary codes together, accordingly achieving integrated compression of trajectory spatio‐temporal data. A series of comparative experiments between the proposed method and representative state‐of‐the‐art methods are carried out on a real massive taxi trajectory dataset from five aspects, and the experimental results indicate that our method has the highest compression ratio. Meanwhile, this method also has favorable performance in other aspects: compression and decompression time overhead, storage space overhead, and historical dataset training time overhead.     

华北克拉通东缘金岭杂岩体岩浆源区及构造背景:来自岩相学、岩石地球化学及年代学的证据

金岭杂岩体由细粒黑云角闪闪长岩、中细粒黑云角闪闪长岩、中细粒黑云角闪二长闪长岩和细粒角闪二长闪长岩组成,是鲁西地区典型矽卡岩型富铁矿(金岭铁矿)控矿岩体。本次研究对细粒黑云角闪闪长岩和细粒角闪二长闪长岩进行了锆石LA‐ICP‐MS U‐Pb定年,其结果分别为129.2±3.2 Ma和132.8±1.2 Ma,表明该杂岩体的侵位时代为早白垩世。样品SiO2、K2O和Na2O含量分别介于54.17%~63.73%、1.92%~4.76%和3.10%~5.41%之间,K2O/Na2O为0.58~0.94,A/CNK为0.60~0.93,具有轻稀土富集,重稀土亏损的右倾型稀土配分模式,轻、重稀土分馏程度中等((La/Yb)N=9.94~23.49),Eu异常不明显(δEu=0.84~1.10),具中—弱负Ce异常(δCe=0.56~0.92)。样品以富集大离子亲石元素(Ba、K、U、Pb等)、亏损高场强元素(Nb、Ta、Ti)以及高Sr/Y为特征。金岭杂岩体为燕山晚期岩浆活动产物,属准铝质高钾钙碱性系列,岩浆主要来源于富集岩石圈地幔的部分熔融,并在岩浆上升侵位的过程中有地壳物质的同化混染。燕山晚期华北克拉通在古太平洋板块俯冲后后撤引起的板内伸展环境下,增厚陆壳减薄阶段,岩浆上侵就位形成金岭杂岩体。     

漠河盆地西部二十二站组碎屑锆石U-Pb年龄、地球化学特征及其物源意义

大兴安岭北部漠河盆地内广泛发育二十二站组,其形成时代及物质来源一直以来备受争议.本文在野外地质调查及岩相学的基础上,利用碎屑锆石LA-ICP-MS U-Pb法测年,并结合岩石地球化学等方法,探讨了二十二站组形成时代、物源区大地构造环境及沉积物来源.测年结果显示132个测点年龄主要分布在4个群落,分别为508.3～492.3 Ma;486.9～435.5 Ma;278.3～254.2 Ma;209.5～178.1 Ma.本次获得二十二站组最小碎屑锆石年龄为178.1 Ma,并结合漠河盆地东西部碎屑锆石的年龄,综合认为二十二站组形成时代应为中晚侏罗世.地球化学分析结果显示,二十二站组砂岩的稀土元素与活动大陆边缘的稀土元素曲线吻合较好,且相对亏损高场强元素(Nb和Ta)和大离子亲石元素(Ba和Sr),相对富集高场强元素(La、Th、Zr、Hf),地球化学特征与活动大陆边缘相似.综合古流向、源区母岩分析、碎屑锆石测年分析、微量元素特征等诸多因素,本次研究认为本区二十二站组砂岩的母岩主要来源于盆地南缘广泛分布的前中生代沉积-变质基底(兴华渡口岩群、佳疙瘩组、额尔古纳河组等)及各期中酸性火成岩.     

Arching effect and displacement on theoretical estimation for lateral force acting on retaining wall

Natural Hazards - A new approach is proposed to evaluate the non-limit active earth pressure in c–φ soil based on the horizontal slices method and limit equilibrium method, which takes...