钱塘江河口细颗粒泥沙絮凝沉降特性研究

钱塘江口河口上游河流和海域来沙多为细颗粒泥沙,粘性细颗粒泥沙由于其特殊的表面电化学性质遇到强电解质海水而产生絮凝沉降,是形成河口淤积的原因之一。影响絮凝的因素很多,除了电解质,还有泥沙粒径的大小、盐度、含沙量、PH值、温度、有机质含量、矿物成分、水流速度及紊动情况等。本文通过粒度分析、静水沉降、动水沉降等各种实验手段分析钱塘江口泥沙的基本特性,找出最佳絮凝盐度以及泥沙不淤流速等值,初步探讨了细颗粒泥沙的絮凝机理,为治理钱塘江口提供科学依据。     

Investigation on Settling Characteristics of Floes in the Changjiang Estuary, China

- Settling characteristics of floes, including relative settling velocity, relative flocculation coefficient and flocculation exponent, are obtained by the suspended load equations for different size fractions. Data of the Changjiang Estuary suggest that level of flocculation changes from river section, river mouth (turbidity maximum) to offshore area in sequence of low, very high and high. The settling characteristics of floes reflected by in situ estimation performs a similar feature as that obtained from still water experiment.     

海上溢油粒子追踪预测模型中的两种数值方法比较

在海上溢油粒子追踪预测模型中,关键的是对拉格朗日微分方程的求解。本文首先通过数值实验比较了欧拉法和龙格-库塔法求解拉格朗日溢油轨迹微分方程的优劣,然后将其应用到2005年4月3日发生在大连附近的“ARTEAGA”油轮溢油事故的油膜粒子追踪模型中。数值实验和应用结果表明,在近岸不均匀流场下,用龙格-库塔方法解拉格朗日油粒子微分方程比用欧拉法求解精度高,用龙格-库塔方法模拟“ARTEAGA”油轮轨迹及其扩散范围与实际观测更为接近,而用欧拉法模拟溢油扩散的面积偏大。     

咸淡水交汇区重金属粒径分配特征的模拟实验

利用从珠江莲花山河段采集的淡水和大洋海水经过滤后，按不同盐度混合，模拟河口咸淡水交汇区重金属在水中各微粒级的分配特征。实验结果表明，在各种盐度条件下，Ｆｅ，Ａｌ，Ｃｕ，Ｐｂ，Ｚｎ和Ｃｄ主要以粒径＞１．０μｍ和＜０．２２μｍ的微粒存在于水中。在加入Ｆｅ（３＋），Ａｌ（３＋）模拟实验中，重金属绝大部分含量形成粒径＞１．０μｍ的微粒，其中Ｆｅ达９９％以上，Ｃｕ，Ｐｂ约９５％，Ａｌ约有７０％－９０％，Ｃｄ约４０％，其余部分主要以＜０．２２μｍ微粒存在。加入与未加入Ｆｅ（３＋），Ａｌ（３＋）的实验比较，Ｆｅ（３＋），Ａｌ（３＋）的加入对水中Ｃｕ，Ｐｂ，Ｚｎ和Ｃｄ含量向大微粒聚集有很大影响，使其聚集在粒径＞１．０μｍ微粒的量增加３０％－４０％。另外从实验结果和理论分析可知，水的盐度与重金属在位径＜１．０μｍ微粒中的含量分布无明显关系，但咸水中大量离子将对河水中胶体絮凝和微粒增长有重大影响。     

Particulate n-alkanes and fatty acids in the Changjiang river system

Particulate samples were collected from the Changjiang river system during a flood period, in May 1997, and POC, stable isotope and lipids associated with particles were examined. Results showed the decrease (0.84% ～ 1.88%) of organic carbon content from the upper reaches to the estuary.δ13C values of particulate organic carbon was in the range of -24.9×10-3 to -26.6×10-3, which were close to the isotopic signature of continental C3 vegetation. Total particulate n-alkanes concentrations varied from 1.4 to 10.1μg/dm3,or from 23.7 to 107μg/g of total suspended matter. Fatty acids were present in all the samples, from 1.4 to 5.4μg/dm3, with saturated and unsaturated straight-chain and branched compounds in the carbon number range from C12 to C30. Both δ13C and the ratio of carbon content to nitrogen content indicate the predominance of terrestrial inputs (soil organic matter) among the particles. The biomarker approach has been used to identify the relative portion of terrigenous and autochthonous fraction in the particulate samples. The distribution of fatty acids suggests a striking phytoplanktonic and microbial signal in most particle samples. The terrestrial alkanes are used to estimate the contribution of terrestrial inputs along the mainstream.     

Numerical simulation of solitary wave propagation based on MPS method

The concept of candidate particle set is introduced in the MPS gridless numerical method to generate neighboring particle set matrix, which can reduce the CPU time to 1/11 of that before introduction. The Bi-CGSTAB (bi-conjugate gradient stabilized) algorithm is applied to solving the Poisson pressure equation, by which the solving speed is significantly accelerated. The process of solitary waves propagating over a numerical flume and interacting with a vertical wall is simulated. The simulated results of water surface elevation are in good agreement with the analytical solution as well as the measured data. The predicted maximum values of the run-up of solitary waves with various relative incident wave heights agree well with the measured results.     

Verification of linear theory for particle velocities in wind waves based on field measurements

Analyses are presented of field measurements of spectral transfer functions between surface elevation and subsurface three-dimensional particle velocity in wind-generated waves, in conditions ranging from young seas to old swells. The results are in agreement with the linear theory predictions to within the measurement error margin, which is estimated to be ± 5% for the gain functions and ± 4° in phase as far as the possible systemic errors are concerned. No correlation is found of the degree of agreement between measurements and linear theory with wave age or wave sleepness.     

Laboratory study of wave and turbulence velocities in a broad-banded irregular wave surf zone

The characteristics of wave and turbulence velocities created by a broad-banded irregular wave train breaking on a 1:35 slope were studied in a laboratory wave flume. Water particle velocities were measured simultaneously with wave elevations at three cross-shore locations inside the surf zone. The measured data were separated into low-frequency and high-frequency time series using a Fourier filter. The measured velocities were further separated into organized wave-induced velocities and turbulent velocity fluctuations by ensemble averaging. The broad-banded irregular waves created a wide surf zone that was dominated by spilling type breakers. A wave-by-wave analysis was carried out to obtain the probability distributions of individual wave heights, wave periods, peak wave velocities, and wave-averaged turbulent kinetic energies and Reynolds stresses. The results showed that there was a consistent increase in the kurtosis of the vertical velocity distribution from the surface to the bottom. The abnormally large downward velocities were produced by plunging breakers that occurred from time to time. It was found that the mean of the highest one-third wave-averaged turbulent kinetic energy values in the irregular waves was about the same as the time-averaged turbulent kinetic energy in a regular wave with similar deep-water wave height to wavelength ratio. It was also found that the correlation coefficient of the Reynolds stress varied strongly with turbulence intensity. Good correlation between u′ and w′ was obtained when the turbulence intensity was high; the correlation coefficient was about 0.3–0.5. The Reynolds stress correlation coefficient decreased over a wave cycle, and with distance from the water surface. Under the irregular breaking waves, turbulent kinetic energy was transported downward and landward by turbulent velocity fluctuations and wave velocities, and upward and seaward by the undertow. The undertow in the irregular waves was similar in vertical structure but lower in magnitude than in regular waves, and the horizontal velocity profiles under the low-frequency waves were approximately uniform.     

Generation of a realistic 3D sand assembly using X‐ray micro‐computed tomography and spherical harmonic‐based principal component analysis

Three‐dimensional particle morphology is a significant problem in the discrete element modeling of granular sand. The major technical challenge is generating a realistic 3D sand assembly that is composed of a large number of random‐shaped particles containing essential morphological features of natural sands. Based on X‐ray micro‐computed tomography data collected from a series of image processing techniques, we used the spherical harmonics (SH) analysis to represent and reconstruct the multi‐scale features of real 3D particle morphologies. The SH analysis was extended to some highly complex particles with sharp corners and surface cavities. We then proposed a statistical approach for the generation of realistic particle assembly of a given type of sand based on the principle component analysis (PCA). The PCA aims to identify the major pattern of the coefficient matrix, which is made up of the SH coefficients of all the particles involved in the analysis. This approach takes into account the particle size effect on the variation of particle morphology, which is observed from the available results of micro‐computed tomography and QICPIC analyses of sand particle morphology. Using the aforementioned approach, two virtual sand samples were generated, whose statistics of morphological parameters were compared with those measured from real sand particles. The comparison shows that the proposed approach is capable of generating a realistic sand assembly that retains the major morphological features of the mother sand. Copyright © 2016 John Wiley & Sons, Ltd.     

Impact of the spatial arrangement of land management practices on surface runoff for small catchments

Predicting the impact of land use changes on the hydrological response is crucial for water resource management. In the particular case of small catchments (1–10 km²), distributed models could provide useful answers regarding the effects of cultivation practices and man‐made works on water fluxes. However, the impacts of specific land use spatial arrangements are difficult to predict because of the prohibitive number of possible cases to consider. Focusing on surface runoff, this article describes a strategy based on a water particle tracking routine to be plugged‐in a distributed model that is designed to determine the spatial arrangements of land management practices that have the greatest impact on volume, peak discharge and lag time at the catchment outlet. A case study is described; the hydrological response of the Roujan catchment (Herault, France) is simulated with the MHYDAS model. The Roujan catchment contains a vineyard in a Mediterranean climate in a landscape in which weeding practices highly influence the partition between soil infiltration and runoff. Results showed that the proposed strategy is much more efficient than a random approach to design the spatial arrangements of the vineyard weeding practices with the greatest impact. Therefore, the proposed strategy may lead to innovative policies for the spatial planning of land management practices. Copyright © 2011 John Wiley & Sons, Ltd.