A comparative study on the illite crystallinity and the clay mineral reflectance spectral index for subdividing the very low-grade metamorphic belt along the Lizhou-Hekou geological section in the Youjiang sedimentary basin, Guangxi, China

To examine the application potential of hyperspectral remote sensing techniques in classifying very low-grade metamorphic belts, the composition of clay minerals and the cyrstallinity of illite from mudstones were measured using XRD and VIS-SWIR (400-2500 nm) reflectance spectroscopy. Based on the illite cyrstallinity, Kubler Index (KI), the Early Triassic LuoLou Group and the Middle Triassic lower Baifeng Formation were classified as the lower Epizone with KI△2θ° ranging from 0.22 to 0.25, the upper Baifeng Formation as upper anchizone with KI△2θ°ranging from 0.26 to 0.33, and the Hekou Formation as lower anchizone with KI△2θ° ranging from 0.38 to 0.40. According to a KI△2θ° value of 0.43, it is possible that there may exist a local diagenetic zone in the upper strata. The illite cyrstallinity Kubler index and the metamorphic grade increase from the bottom to the top of the stratigraphic sequence. The metamorphic grade boundaries nearly match the stratigraphic boundaries, indicating a burial metamorphism nature for the stratigraphic sequence. From the bottom to the top of the sequence, the spectral absorption band center of clay minerals from fresh rocks is around 2200 nm. The absorption band centers change towards shorter wavelengths: the Luolou Group being at 2220 nm, the Baifeng Formation at 2217-2213 nm, the lower member of the Hekou Formation at 2214-2206 nm, and the upper member of the Hekou Formation at 2205-2197 nm. The spectral absorption band center of illite shows the same change pattern. These results indicate that very low-grade metamorphic belts can be subdivided using spectral indices of clay minerals, which are measured by using field portable spectroradiometers. However, it may not work well with satellite and airborne sensors.     

Fitting the Linear Model of Coregionalization by Generalized Least Squares

In geostatistical studies, the fitting of the linear model of coregionalization (LMC) to direct and cross experimental semivariograms is usually performed with a weighted least-squares (WLS) procedure based on the number of pairs of observations at each lag. So far, no study has investigated the efficiency of other least-squares procedures, such as ordinary least squares (OLS), generalized least squares (GLS), and WLS with other weighing functions, in the context of the LMC. In this article, we compare the statistical properties of the sill estimators obtained with eight least-squares procedures for fitting the LMC: OLS, four WLS, and three GLS. The WLS procedures are based on approximations of the variance of semivariogram estimates at each distance lag. The GLS procedures use a variance–covariance matrix of semivariogram estimates that is (i) estimated using the fourth-order moments with sill estimates (GLS₁), (ii) calculated using the fourth-order moments with the theoretical sills (GLS₂), and (iii) based on an approximation using the correlation between semivariogram estimates in the case of spatial independence of the observations (GLS₃). The current algorithm for fitting the LMC by WLS while ensuring the positive semidefiniteness of sill matrix estimates is modified to include any least-squares procedure. A Monte Carlo study is performed for 16 scenarios corresponding to different combinations of the number of variables, number of spatial structures, values of ranges, and scale dependence of the correlations among variables. Simulation results show that the mean square error is accounted for mostly by the variance of the sill estimators instead of their squared bias. Overall, the estimated GLS₁ and theoretical GLS₂ are the most efficient, followed by the WLS procedure that is based on the number of pairs of observations and the average distance at each lag. On that basis, GLS₁ can be recommended for future studies using the LMC.     

中巴地球资源一号卫星CCD图像质量评价和交叉定标研究

使用定量分析方法评价了中巴地球资源卫星CCD相机遥感数据质量 ,在此基础上利用LANDSAT 7ETM 对CCD的对应波段进行交叉定标计算 ,得到了一组CCD的参考定标系数。在遥感数据质量评价研究中 ,使用统计分析方法计算得到了CCD相机图像数据的条纹强度 ;应用结构函数方法计算获得了CCD的噪声数据 ;通过功率谱计算 ,对比分析了CCD和LANDSAT 7ETM 对应波段图像空间纹理特征。该文通过交叉定标计算 ,得到了CCD相机前 4个波段的参考定标系数、动态范围和噪声等效反射率等参考指标 ,并对这组定标系数进行了验证 ,从验证结果中发现CCD1,2 ,3和LANDSAT 7ETM 对应波段反射率最大相差 1 98% ,最小0 0 3% ,平均差要低于 1 31% ,而波段 4最大也仅相差 4 4 1% ,平均差 3 0 2 % ,证明计算得到的交叉定标系数具有较高的参考价值     

元江-红河跨境水质环境问题和保护研究

在红河干流不同污染流域,设立5个环境监测断面,同时选定人烟稀少、无工业、自然条件良好的李仙江加禾断面为对照。利用1999～2002年对水质影响大的DO、COD、BOD5、NH 3-N、NO3-N、NO 2-N、Cu、Pb、Cd、T-P、石油类及SS 12项污染物,均在枯、丰、平3个不同水期采样监测,将监测结果做数理统计分析,研究污染物的迁移、转化规律,并得出:出境断面NH 3-N平均含量高达0.18～0.86m gL/,T-P高达0.045～0.315m g/L,SS高达69～3412m gL/,从而提出小流域治理、退耕还林、城市水网改造及控制含磷洗涤剂使用的建议。     

海峡两岸旅游管理比较研究

中国大陆和台湾的旅游业都有良好的发展前景,由于体制和区情的不同,两地在旅游发展和旅游管理上有着诸多差异;由于文化上的同根同源,两地在旅游发展与管理上也存在着互补共融的关系.如何处理好竞争与合作、经济利益和可持续发展之间的关系,成为两岸学术界和旅游业工作者共同关注的课题.通过重点比较两岸在旅游管理方面的特点,找出了两岸的差异及差异存在的部分原因,以期能为两岸旅游业取长补短、互相促进、共同发展提供有价值的参考.     

基于可持续发展的广东省城市供水价格改革探讨

针对目前城市供水价格计算与计量方式不合理和突出的水资源供需矛盾,提出了按全成本水价模式测算广东省城市供水价格,并采用阶梯水价、两部制水价和季节水价计量方式逐步达到全成本水价的措施与实施步骤,确保水资源的持续利用支撑广东省经济社会的可持续发展.     

Three‐dimensional shakedown solutions for anisotropic cohesive‐frictional materials under moving surface loads

Previous work on three‐dimensional shakedown analysis of cohesive‐frictional materials under moving surface loads has been entirely for isotropic materials. As a result, the effects of anisotropy, both elastic and plastic, of soil and pavement materials are ignored. This paper will, for the first time, develop three‐dimensional shakedown solutions to allow for the variation of elastic and plastic material properties with direction. Melan's lower‐bound shakedown theorem is used to derive shakedown solutions. In particular, a generalised, anisotropic Mohr–Coulomb yield criterion and cross‐anisotropic elastic stress fields are utilised to develop anisotropic shakedown solutions. It is found that shakedown solutions for anisotropic materials are dominated by Young's modulus ratio for the cases of subsurface failure and by shear modulus ratio for the cases of surface failure. Plastic anisotropy is mainly controlled by material cohesion ratio, the rise of which increases the shakedown limit until a maximum value is reached. The anisotropic shakedown limit varies with frictional coefficient, and the peak value may not occur for the case of normal loading only. Copyright © 2013 John Wiley & Sons, Ltd.     

东海PN断面透明胞外聚合颗粒物分布特征及来源研究

透明胞外聚合颗粒物（TEP）是凝聚网的重要组成部分甚至核心，在微尺度范围上形成了水环境结构的重要梯度，对于元素地球化学生物循环、碳沉降以及食物网有着举足轻重的作用。本文研究了东海典型断面PN透明胞外聚合颗粒物的分布特征及来源。结果表明:东海典型断面PN透明胞外聚合颗粒物含量介于28~376 μg Xeq./L之间，平均值为（115±67）μg Xeq./L；呈现出明显的夏季 > 冬季 > 秋季 > 春季的季节变化特征以及近岸 > 外海和底层 > 表层的分布趋势。通过对比透明胞外聚合颗粒物分布趋势和硅藻、甲藻分布趋势以及统计分析得出，东海典型断面PN透明胞外聚合颗粒物主要来源于硅藻，甲藻贡献不大；而外海TEP的来源可能主要由超微型浮游植物贡献。     

Nonplanar multi-modal vibrations of fluid-conveying risers under shear cross flows

Risers/pipes conveying fluid are a typical kind of slender structures commonly used in marine engineering. It is of great academic significance and application value for us to evaluate and understand the vibration characteristics and nonlinear responses of these risers under the combined action of internal and external fluid flows. In this paper, the nonplanar vibrations and multi-modal responses of pinned-pinned risers in shear cross flow are numerically studied. With this objective in mind, the van der Pol wake oscillators are used to simulate the dynamical behavior of the vortex shedding in the wake. Two nonlinear equations of motion of the riser are proposed to govern the lateral responses of the riser structure. The nonplanar nonlinear equations for the riser and wake are then discretized by employing Galerkin's method and solved by using a fourth-order Runge–Kutta integration algorithm. Theoretical results show that the coupled frequencies for cross-flow (CF) and in-line (IL) motions and the corresponding coupled damping ratio could be influenced by the external and/or internal fluid velocities. Based on extensive calculations, the dynamical behavior of the riser with various internal and external flow velocities are presented in the form of bifurcation diagrams, time traces, phase portraits, oscillation trajectories and response spectrum curves. It is shown that some interesting dynamical phenomena, such as ‘lock-in’ state, ‘figure-of-eight’ trajectory and quasi-periodic oscillation, could occur in such a fluid-structure interaction system. Our results also demonstrate that the shear parameter can significantly affect the dynamic responses of the riser. When the shear parameter of the cross flow is large, multi-modal quasi-periodic responses of the riser can be excited, showing some new features undetected in the system of fluid-conveying risers in uniform cross flow.     

Dynamic response of a pile considering the interaction of pile variable cross section with the surrounding layered soil

Assuming that the pile variable cross section interacts with the surrounding soil in the same way as the pile toe does with the bearing stratus, the interaction of pile variable cross section with the surrounding soil is represented by a Voigt model, which consists of a spring and a damper connected in parallel, and the spring constant and damper coefficient are derived. Thus, a more rigid pile–soil interaction model is proposed. The surrounding soil layers are modeled as axisymmetric continuum in which its vertical displacements are taken into account and the pile is assumed to be a Rayleigh–Love rod with material damping. Allowing for soil properties and pile defects, the pile–soil system is divided into several layers. By means of Laplace transform, the governing equations of soil layers are solved in frequency domain, and a new relationship linking the impedance functions at the variable‐section interface between the adjacent pile segments is derived using a Heaviside step function, which is called amended impedance function transfer method. On this basis, the impedance function at pile top is derived by amended impedance function transfer method proposed in this paper. Then, the velocity response at pile top can be obtained by means of inverse Fourier transform and convolution theorem. The effects of pile–soil system parameters are studied, and some conclusions are proposed. Then, an engineering example is given to confirm the rationality of the solution proposed in this paper. Copyright © 2017 John Wiley & Sons, Ltd.