非均匀各向异性储层地震AVAZ响应模拟及分析

基于波动理论的各向异性反射率方法能够充分考虑地震波动力学特征。对单界面理论模型、实际非均匀各向异性模型分别进行正演模拟及分析,比较Ruger理论与反射率方法模拟结果的波动力学特征。结果表明,对于界面模型两种方法均能得到相同的计算结果,但对于实际非均匀各向异性介质,基于Ruger理论的算法受时间采样率的影响,而基于波动理论的反射率方法输入实际空间域模型,计算结果不受时间采样率的限制。     

北淮阳东段西汤池岩体地球化学特征、锆石U-Pb定年及成因

西汤池岩体位于安徽北淮阳地区舒城县西汤池一带。西汤池似斑状二长花岗岩具高硅、富碱质特征,为铝饱和碱性系列岩石。稀土元素含量中等,(La/Yb)_N与HREE/LREE值较高,重稀土相对轻稀土亏损明显,Eu呈弱负异常,属轻稀土富集型。大离子亲石元素(Rb、K、La、Nd等)明显富集,而高场强元素(Nb、Ta、P、Ti等)亏损。通过LA-ICP-MS锆石U-Pb定年,获得西汤池二长花岗岩~(206)Pb/~(238)U年龄为(125.5±1.6)Ma,为早白垩世岩浆活动的产物。基于前人关于北淮阳地区区域地质背景方面的研究成果和本次地球化学构造环境判别方面的研究,认为西汤池二长花岗岩属A_1型花岗岩,形成于后碰撞的伸展构造环境。     

Sinian-Cambrian Tectonic Framework in the Tarim Basin and its Influences on the Paleogeography of the Early CambrianEI北大核心CSCD

Due to strong remodification and sparse trace remnants, tectonic restoration is critical for the reconstructing of the paleography and basin-prototype in the multiple-cycles superimposed basin, which is also fundamental for the hydrocarbon exploration of the deep marine carbonate in the Tarim Basin. This study presents the tectonic framework of the transition phase from the Sinian to the Cambrian Period and its constrains on the paleography of the Cambrian in the Tarim Basin based on comprehensive analysis of new geochronology data, seismic data and regional tectonic data. The results showed that (1) there is a large regional unconformity between the Cambrian and Precambrian, suggesting a discontinuous deposition from the Sinian to the Cambrian that is contrary to the major view; (2) a broad flattened paleotopography formed before Cambrian Period, which is favorable for the wide epicontinental sea environment and a gentle homoclinal ramp platform development in the early Cambrian, and the residual basement low relief uplifts that influence the microfacies differentiation in the carbonate platform; (3) the Cambrian carbonate platform has been involved much into the marginal orogenic belts, and the proto-platform is probably of 200 km out of the present basin in S-N direction; (4) there is a weak extensional setting in the Cambrian-Early Ordovician rather than a strong rifting setting with a large aulacogen into the platform, in which a inter-platform shallow depression and stable ramp platform developed in the Cambrian; (5) there is not troughs, but rather low relief uplifts developed in the Early Cambrian in the southwestern Tarim Basin, and lack of large uplifts along the carbonate platform margins; (6) it is showed a broad architecture of platform-wide slope-ocean basin from the inner Tarim plate to its margin in the Cambrian Period. Considering the basin prototype, the restoration of the tectono-paleography in the early Cambrian in the Tarim Basin is distinct from those of the previously proposed. We propose that a gentle pattern with “two platforms and one inter-depression” in E-W striking other than multiple small platform constrained in the basin, and further subdivision of the western platform with inner platform low uplift and sag in N-S striking, which is possibly inherited from the basement architecture. On the outer of the united platform, there is probably broad gentle slope for transition to ocean basin. We therefore argue that tectonic restoration is crucial for reconstruction of the paleogeography and basin prototype, especially for basins experienced multiphase of tectonic cycles. © 2018, Science Press. All right reserved.     

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.     

Research Progress and Prospect on Reactivation of Ancient Landslides

The ancient landslide has endured long-term slope evolution which results in its complicated material and special rock-soil properties. The risk of ancient landslide reactivation is substantially increasing due to the increase of intensified human engineering activities and the frequency of extreme weather events. Many ancient landslides have been reactivated all over the world and led to serious fatalities and severe damage to many important engineering facilities such as transportation and hydropower engineering projects. On the basis of the analysis of the research situation about the ancient landslides at home and abroad, the main research advances were summarized including the regional developing laws and recognizing of the ancient landslides, the mechanics properties of ancient landslide body and related sliding zone, reactivation mechanism of ancient landslides, reactivating process and modeling analysis of ancient landslides, early recognization of ancient landslide reactivation, etc. To meet the demands of disaster prevention and reduction, three key scientific issues were put forward to be solved: ①automaticaly establishing the methodology and identification criterions for recognition of ancient landslide; ②revealing the reactivation mechanism of ancient landslide based on a new strength theory; ③establishing the early rapid recognition method and predictive model for ancient landslide reactivation. Solving the above mentioned scientific theory and methodology will facilitate the planning and site selection of major projects as well as the disaster prevention and reduction in ancient landslide developing areas.     

一种盾构施工引起的地面沉降预测方法

针对盾构施工引起的地面沉降的影响因素众多、关系错综复杂,常规数学模型难以准确预测沉降量的问题,该文提出了采用主成分分析法和粒子群优化的支持向量机方法来建立预测模型。并结合工程实例将预测结果与常用的多元线性回归模型和基于Levenberg-Marquardt算法改进的BP神经网络模型(LM-BP)的预测结果做了比较。结果表明:PCA-PSO-SVM的预测结果精度较多元线性回归模型和LM-BP神经网络有很大的提高,证明了该研究方法具有一定的理论意义和参考价值。     

Beidou satellite maneuver thrust force estimation for precise orbit determination

Beidou satellites, especially geostationary earth orbit (GEO) and inclined geosynchronous orbit (IGSO) satellites, need to be frequently maneuvered to keep them in position due to various perturbations. The satellite ephemerides are not available during such maneuver periods. Precise estimation of thrust forces acting on satellites would provide continuous ephemerides during maneuver periods and could significantly improve orbit accuracy immediately after the maneuver. This would increase satellite usability for both real-time and post-processing applications. Using 1 year of observations from the Multi-GNSS Experiment network (MGEX), we estimate the precise maneuver periods for all Beidou satellites and the thrust forces. On average, GEO and IGSO satellites in the Beidou constellation are maneuvered 12 and 2 times, respectively, each year. For GEO satellites, the maneuvers are mainly in-plane, while out-of-plane maneuvers are observed for IGSO satellites and a small number of GEO satellites. In most cases, the Beidou satellite maneuver periods last 15–25 min, but can be as much as 2 h for the few out-of-plane maneuvers of GEO satellites. The thrust forces acting on Beidou satellites are normally in the order of 0.1–0.7 mm/s². This can cause changes in velocity of GEO/IGSO satellites in the order of several decimeters per second. In the extreme cases of GEO out-of-plane maneuvers, very large cross-track velocity changes are observed, namely 28 m/s, induced by 5.4 mm/s² thrust forces. Also, we demonstrate that by applying the estimated thrust forces in orbit integration, the orbit errors can be estimated at decimeter level in along- and cross-track directions during normal maneuver periods, and 1–2 m in all the orbital directions for the enormous GEO out-of-plane maneuver.     

A new multipath mitigation method based on adaptive thresholding wavelet denoising and double reference shift strategy

Multipath disturbance is one of the major error sources in high-accuracy positioning for global navigation satellite system (GNSS). Although various methods based on software and hardware strategies have been developed to mitigate this error, they are still limited by different kinds of factors and the effect is not ideal. After analyzing the existing methods, a new single-difference sidereal filtering method, based on adaptive thresholding wavelet denoising and double reference shift strategy (ATDR), is proposed to mitigate multipath effects for static short-baseline GNSS applications. The key idea of the proposed method is the use of both the adaptive thresholding wavelet denoising to extract an accurate multipath correction model from the reference Day and the double reference shift strategy to mitigate multipath for subsequent Day 2 more accurately and efficiently. By applying the introduced adaptive thresholding method, the average improvement rate of the RMS values of the single-difference residuals can reach about 15.82% compared with the constant thresholding method. Moreover, after applying the proposed ATDR method, the 3D positioning precision is improved by about 37.73% for the single epoch mode with 30 s data sampling rate and about 31.22% for the continuous mode with 1 s high sampling rate compared with the original results. Even compared with the constant thresholding single orbital reference (CTSR) method, the improvement percentage is about 33.94% in single epoch mode and about 25.40% in continuous mode for 3D positioning precision, respectively. In conclusion, the results of the two experiments indicate that the proposed ATDR method performs much better than the CTSR method in mitigating multipath for different sampling rates and different processing modes in the measurement domain for GNSS static short-baseline postprocessing applications.     

Revisit the calibration errors on experimental slant total electron content (TEC) determined with GPS

The calibration errors on experimental slant total electron content (TEC) determined with global positioning system (GPS) observations is revisited. Instead of the analysis of the calibration errors on the carrier phase leveled to code ionospheric observable, we focus on the accuracy analysis of the undifferenced ambiguity-fixed carrier phase ionospheric observable determined from a global distribution of permanent receivers. The results achieved are: (1) using data from an entire month within the last solar cycle maximum, the undifferenced ambiguity-fixed carrier phase ionospheric observable is found to be over one order of magnitude more accurate than the carrier phase leveled to code ionospheric observable and the raw code ionospheric observable. The observation error of the undifferenced ambiguity-fixed carrier phase ionospheric observable ranges from 0.05 to 0.11 total electron content unit (TECU) while that of the carrier phase leveled to code and the raw code ionospheric observable is from 0.65 to 1.65 and 3.14 to 7.48 TECU, respectively. (2) The time-varying receiver differential code bias (DCB), which presents clear day boundary discontinuity and intra-day variability pattern, contributes the most part of the observation error. This contribution is assessed by the short-term stability of the between-receiver DCB, which ranges from 0.06 to 0.17 TECU in a single day. (3) The remaining part of the observation errors presents a sidereal time cycle pattern, indicating the effects of the multipath. Further, the magnitude of the remaining part implies that the code multipath effects are much reduced. (4) The intra-day variation of the between-receiver DCB of the collocated stations suggests that estimating DCBs as a daily constant can have a mis-modeling error of at least several tenths of 1 TECU.     

An enhanced approach for surface flow routing over drainage‐constrained triangulated irregular networks

The accuracy and efficiency of the simulations in distributed hydrological models must depend on the proper estimation of flow directions and paths. Numerous studies have been carried out to delineate the drainage patterns based on gridded digital elevation models (DEMs). The triangulated irregular network (TIN) has been increasingly applied in hydrological applications due to the advantages of high storage efficiency and multi‐scale adaptive performance. Much of the previous literature focuses mainly on filling the depressions on gridded DEMs rather than treating the special cases in TIN structures, which has hampered its applications to hydrological models. This study proposes a triangulation‐based solution for the removal of flat areas and pits to enhance the simulation of flow routing on triangulated facet networks. Based on the drainage‐constrained TIN generated from only a gridded DEM by the compound point extraction (CPE) method, the inconsistent situations including flat triangles, V‐shape flat edges and sink nodes are respectively identified and rectified. The optimization algorithm is an iterative process of TIN reconstruction, in which the flat areas are generalized into their center points and the pits are rectified by embedding break lines. To verify the proposed algorithm and investigate the potential for flow routing, flow paths of steepest descent are derived by the vector‐based tracking algorithm based on the optimized TIN. A case study of TIN optimization and flow path tracking was performed on a real‐world DEM. The outcomes indicate that the proposed approach can effectively solve the problem of inconsistencies without a significant loss in accuracy of the terrain model.