地下水开采引发的土体变形对污染物迁移影响研究

地下水开采引发的土体变形对污染物迁移具有重要影响。结合地下水渗流理论、太沙基一维固结理论和溶质运移对流-弥散理论,建立了变形土体中污染物迁移三维耦合数值模型,考虑了孔隙度、渗透系数和水动力弥散系数随有效应力的动态变化。通过对实际案例进行模拟,分析了地下水开采引发的土体变形对污染物迁移规律的影响。结果表明,地下水位下降,土体有效应力增加,土体骨架压缩变形,使得孔隙度、渗透系数减小,导致水动力弥散系数增幅减小,从而延缓了污染物迁移的对流扩散过程。     

饱和软土大变形非线性自重固结模型

基于分段线性方法,建立了饱和软土一维自重固结模型（简称SWC模型）。该模型能考虑自重固结过程中土体的大变形效应和材料参数的非线性变化。将该模型的计算结果与相关解析解、现场试验及室内试验结果进行了对比验证,证明了SWC模型能准确计算出大变形和非线性条件下饱和软土的自重固结过程,包括沉降量、平均固结度、孔隙比分布和超孔隙水压力分布等参数随时间的变化过程。随后,以现场试验为基础,采用SWC模型对饱和软土自重固结的4个主要影响因素（即土体初始高度、边界排水条件、初始孔隙比和土粒相对密度）进行了参数分析。结果表明,上述4个参数对软土自重固结过程均具有重要影响：土体初始高度越高,则自重固结沉降量和最终平均应变值越大;边界排水条件对土体自重固结的速度有重要影响,但不影响自重固结的最终沉降量;初始孔隙比越大,则自重固结沉降量越大,其完成自重固结所需时间越短;土粒相对密度越大,则土体的最终沉降量越大,完成自重固结所需时间越短。     

饱和软土大变形非线性自重固结模型

基于分段线性方法,本文建立了饱和软土一维自重固结模型(简称SWC模型)。该模型能考虑自重固结过程中土体的大变形效应和材料参数的非线性变化。将该模型的计算结果与相关解析解、现场试验及室内试验结果进行了对比验证,证明了SWC模型能准确计算出大变形和非线性条件下饱和软土的自重固结过程,包括沉降量、平均固结度、孔隙比分布和超孔隙水压力分布等参数随时间的变化过程。随后,以现场试验为基础,本文采用SWC模型对饱和软土自重固结的四个主要影响因素(即土体初始高度,边界排水条件,初始孔隙比和土粒比重)进行了参数分析。结果表明,上述四个参数对软土自重固结过程均具有重要影响:土体初始高度越高,则自重固结沉降量和最终平均应变值越大;边界排水条件对土体自重固结的速度有重要影响,但不影响自重固结的最终沉降量;初始孔隙比越大,则自重固结沉降量越大,其完成自重固结所需时间越短;土粒比重越大,则土体的最终沉降量越大,完成自重固结所需时间越短。     

饱和软土大变形非线性自重固结模型

基于分段线性方法,本文建立了饱和软土一维自重固结模型(简称SWC模型)。该模型能考虑自重固结过程中土体的大变形效应和材料参数的非线性变化。将该模型的计算结果与相关解析解、现场试验及室内试验结果进行了对比验证,证明了SWC模型能准确计算出大变形和非线性条件下饱和软土的自重固结过程,包括沉降量、平均固结度、孔隙比分布和超孔隙水压力分布等参数随时间的变化过程。随后,以现场试验为基础,本文采用SWC模型对饱和软土自重固结的四个主要影响因素(即土体初始高度,边界排水条件,初始孔隙比和土粒比重)进行了参数分析。结果表明,上述四个参数对软土自重固结过程均具有重要影响:土体初始高度越高,则自重固结沉降量和最终平均应变值越大;边界排水条件对土体自重固结的速度有重要影响,但不影响自重固结的最终沉降量;初始孔隙比越大,则自重固结沉降量越大,其完成自重固结所需时间越短;土粒比重越大,则土体的最终沉降量越大,完成自重固结所需时间越短。     

盾构施工引起地表固结沉降问题的研究

将隧道周围土体视为均质连续各向同性的饱和弹性介质,采用保角变换的方法将含有隧道的半无限平面映射为同心圆环计算域。根据Terzaghi-Rendulic二维固结理论,建立隧道在不透水的情况下周围土体超孔隙水压力分布的控制方程。然后,采用分离变量法计算得到土体超孔隙水压力分布的解析解,最后,根据弹性理论计算得出隧道中线上方地表固结沉降的计算公式。结合算例,分析盾构施工扰动程度、土体渗透系数、土体弹性模量及隧道埋深等因素对隧道中心上方地表处固结沉降的影响。研究结果表明,地表固结沉降的增加值与隧道外侧初始超孔隙水压力值C0的变化量成正比例关系,施工扰动程度越大所引起的固结沉降越大;土体的渗透系数越大固结沉降速度越快,但土体的渗透系数与最终的地表固结沉降量无关;土体的弹性模量越大,最终的地表固结沉降量越小;隧道埋深越深,地表固结沉降所需时间越长,最终的地表固结沉降量也越大。     

饱和软土成层地基一维非线性固结解析解

假定土体在固结过程中压缩性和渗透性的变化成正比,基于 - 及 - 关系,推导出饱和软土成层地基一维非线性固结解析解,分别给出了按沉降定义和按有效应力定义的每层土平均固结度及整个土层总固结度的计算公式。采用Fortran语言编制了相应的计算程序,将计算得到的结果与已有双层地基一维非线性固结解析解计算结果进行比较,验证该解析解的正确性。利用该程序分析成层地基一维非线性固结性状,分别讨论了初始竖向渗透系数、初始体积压缩系数、荷载值及土层厚度对地基固结性状的影响。分析结果表明：在成层地基一维非线性固结过程中,初始竖向渗透系数对超静孔压的影响较为复杂,对上下层地基固结速率影响不同;初始体积压缩系数增大,超静孔压增大,固结速率变小;所加荷载值越大,超静孔压消散越慢,固结发展越慢;超静孔压消散速率不仅取决于土层厚度,同时取决于各层土渗透性的相对大小。     

非等温分布条件下考虑半透水边界时饱和黏土的一维固结解析解

温度的变化会导致土体的物理力学性质改变,且在一些实际工程中,饱和黏土会处于非等温分布状态。为此,针对非等温分布条件下饱和黏土的一维固结问题,考虑了更具普遍性的半透水边界,通过某些假定推导了单级线性加荷形式下饱和黏土一维固结控制方程,并利用分离变量法求解得到了控制方程的解析解。通过将所提解析解分别与已有解析解和有限差分解展开对比分析,验证了所提解答的正确性。基于所提解析解,利用某一算例分析了温度梯度、半透水边界参数及加荷时间对固结性状的影响。结果表明：温度梯度 M 越大,土体的渗透性越大,土体的固结速率越快;半透水边界参数 R1和 R2越大,相同时间内土体的超孔隙水压力越小,土体的平均固结度越大;土体的平均固结度随加荷时间 tc 的增大而减小,这主要是由于加荷阶段所施加的外荷载小于最终荷载,但加荷时间tc的延长可一定程度减小土体中产生的最大超孔隙水压力。     

考虑涂抹区压缩及渗透系数变化的堆载预压固结解

为评估涂抹区土体压缩和渗透系数变化对含竖向排水体地基固结的影响,采用等体积应变假设,考虑涂抹区土体的压缩变形及其水平向渗透系数沿径向分别呈线性和抛物线分布,并考虑井阻作用以及地基附加球应力沿深度任意分布,推导了随时间线性堆载预压条件下固结微分方程的显式解析解答,分析了涂抹区半径、水平向渗透系数的分布模式、以及体积压缩系数对地基整体平均固结度的影响。结果表明,涂抹区土体采用均匀折减的水平向渗透系数明显低估了地基的固结速率,而当涂抹区半径较大时,不考虑涂抹区土体的压缩变形将会高估地基的固结速率。在含竖向排水体地基固结问题的分析中,这些影响不可忽视。     

饱和软土大应变自重固结非线性分析

饱和软土的自重固结涉及许多复杂的非线性问题,而现有分析方法一般基于线性假设,将固结控制方程简化,在沉降变形较大的情况下并不合理。基于Gibson大变形固结方程,推导在单面排水和双面排水条件下的边界条件,综合考虑饱和软土的自重特性与非线性压缩、渗透关系,采用稳定性较好的Crank-Nicolson型差分格式离散并求解非线性控制方程,并从孔隙比、沉降、孔隙水流速3个方面将数值解与CS2固结数值模型进行对比验证。结合工程实例,研究土层不同初始厚度和初始孔隙比对自重固结的影响。结果表明：初始厚度与初始孔隙比较大土层最终沉降也较大。与传统方法相比,该方法所得出的数值解精度较高,实现了高度非线性固结方程的快速求解,其成果更能广泛运用于实际工程的分析之中。     

考虑影响区真实形状的竖向排水井地基固结解

针对现有竖向排水井地基固结解析理论对影响区一律采用圆形等效假定的缺陷,研究了影响区真实形状为正六边形的按梅花形布置的竖向排水井地基的固结问题。通过建立新固结方程和引入新的边界条件,并考虑土体水平渗透系数变化,得到了相应的解析解。对于土体水平渗透系数的3种变化模式,分别给出了各种模式下的特殊解。在此基础上,分析了3种模式下3个主要的无量纲参数对地基固结性状的影响,并比较了计算结果和现有理论结果。分析结果表明,影响区和涂抹扰动区范围越大,固结越慢;土体的最大与最小水平渗透系数之比越小,固结越快;在相同条件下,考虑扰动区渗透系数线性变化的模式2固结最快,而假定扰动区渗透系数不变的模式1的解与现有理论解相当接近,验证了现有竖向排水井地基固结解析理论中对影响区采用圆形等效假定的合理性     

Distribution and enrichment of trace metals in marine sediments of Bay of Bengal,off Ennore,south-east coast of India

In order to avoid the pollution of trace metals in marine environment, it is necessary to establish the data and understand the mechanisms influencing the distribution of trace metals in marine environment. The concentration of heavy metals (Fe, Mn, Cr, Cu, Ni, Pb, Zn, Co and Cd) were studied in sediments of Ennore shelf, to understand the metal contamination due to heavily industrialized area of Ennore, south-east coast of India. Concentration of metals shows significant variability and range from 1.7 to 3.7% for Fe, 284–460 μg g⁻¹ for Mn, 148.6–243.2 μg g⁻¹ for Cr, 385–657 μg g⁻¹ for Cu, 19.8–53.4 μg g⁻¹ for Ni, 5.8–11.8 μg g⁻¹ for Co, 24.9–40 μg g⁻¹ for Pb, 71.3–201 μg g⁻¹ for Zn and 4.6–7.5 μg g⁻¹ for Cd. For various metals the contamination factor (CF) and geoaccumulation index (I _geo) has been calculated to assess the degree of pollution in sediments. The geoaccumulation index shows that Cd, Cr and Cu moderately to extremely pollute the sediments. This study shows that the major sources of metal contamination in the Ennore shelf are land-based anthropogenic ones, such as discharge of industrial wastewater, municipal sewage and run-off through the Ennore estuary. The intermetallic relationship revealed the identical behavior of metals during its transport in the marine environment.     

Kinetics of dehydration of Ca-montmorillonite

Isothermal thermogravimetric experiments have been carried out to determine the reaction kinetics of the dehydration processes in fuller's earth, a natural Ca-montmorillonite. Dehydration in swelling clays is a complex reaction, and analysis of the thermogravimetric data using empirical rate equations and time-transformation analysis reveals that the nature of the rate controlling mechanism is dependent upon both the temperature regime of the sample as well as the extent of reaction. For fuller's earth, we find that the dehydration kinetics are dominated by a nucleation and growth mechanism at low temperatures and fractions transformed (stage I), but above 90 °C the last stages of the reaction are diffusion controlled (stage II). The activation energy for dehydration during stage I is around 35 kJ · mol⁻¹, whereas the removal of water during stage II requires an activation energy of around 50 kJ · mol⁻¹. These two stages of dehydration are associated with primary collapse of the interlayer (stage I) and movement of water that is hydrated to cations within the interlayer (stage II). Received: 28 August 1998 / Revised, accepted: 27 January 1999     

A review of the stratigraphy of Marwar Supergroup of west-central Rajasthan

The lithostratigraphy, depositional environment and age of the Marwar Supergroup have been reviewed in the light of report of δ¹³C depletion recorded in the carbonates of the Bilara Group (middle part of Marwar Supergroup) and discovery of trilobite-like trace fossils from the ·Red bedsŽ of Nagaur Group (upper part of Marwar Supergroup). The δ¹³C depletion observed in Bilara carbonates is not a result of glaciation rather due to rapid burial and poor water circulation in the low energy water of the protected basin. Secondly, the trace fossils are, in fact, traces of notostracan crustaceans found in shallow fluvial and shallow lacustrine environment. The present paper also records a spiral, burrowing trace-fossil, possibly Gyrolithes, from a cross-bedded sandstone of the Jodhpur Group.     

P-V-T equation of state of lawsonite

A pressure-volume-temperature data set has been obtained for lawsonite [CaAl₂Si₂O₇(OH)₂.H₂O], using synchrotron X-ray diffraction and an externally heated diamond anvil cell. Unit-cell volumes were measured to 9.4 GPa and 767 K by angle dispersive X-ray diffraction using imaging plates. Phase changes were not observed within this pressure-temperature range, and lawsonite compressed almost isotropically at constant temperature. The P-V-T data have been analyzed using a Birch- Murnaghan equation of state and a linear equation of state expressed as β=–1/V₀ (∂V/∂P) _T . At room temperature, the derived equation of state parameters are: K ₀=124.1 (18) GPa K'₀ set to 4) and β^–1=142.0(24) GPa, respectively. Our results are intermediate between previously reported measurements. The high-temperature data show that the incompressibility of lawsonite decreases with increasing temperature to ∼500 K and then increases above. Hence, the second order temperature derivative of the bulk modulus is taken into account in the equation of state; a fit of the volume data yields K ₀=123.9(18) GPa, (∂K/∂T)_P=–0.111(3) GPa K^–1, (∂² K/∂T ²)_P=0.28(6) 10^–3 GPa K^–2, α₀=3.1(2) 10^–5 K^–1, assuming K'₀=4. Received: 2 June 1998 / Revised, accepted: 12 Ocotber 1998     

Thermomineral waters of Ca-chloride composition: review of diagnostics and of brine evolution

 Ca-chloride waters are defined as those in which Q=rCa/r(SO₄+HCO₃)>1, rNa/rCl<0.80, rMg/rCa<0.5 and wCl/wBr<286 (r=meq l^–1 and w=mg l^–1). Throughout the last 50 years, different models for the formation of such waters have been suggested. These models include: (1) filtration through semipermeable membranes under conditions of highly compacted argillaceous sediments, (2) deaquation of seawater by evaporation and/or by freezing followed by dolomitization, (3) hydrolysis of plagioclase and biotites in igneous metamorphic rock masses, (4) radiolytic modification of residual metamorphic fluids, and (5) dissolution of chalks followed by ion exchange on smectites. The better understanding of processes and of natural environments leading to the evolution and natural occurrence of such brines, is imperative for the prospection and further sustained exploitation of such waters. Received: 11 October 1996 · Accepted: 24 February 1997     

Effect of incorporation of iron and aluminum on the thermoelastic properties of magnesium silicate perovskite

In situ X-ray diffraction measurements of Fe- and Al-bearing MgSiO₃-rich perovskite (FeAl-Pv), which was synthesized from a natural orthopyroxene, were performed at pressures of 19–32 GPa and temperatures of 300–1,500 K using a combination of a Kawai-type apparatus with eight sintered-diamond anvils and synchrotron radiation. Two runs were performed using a high-pressure cell with two sample chambers, and both MgSiO₃ perovskite (Mg-Pv) and FeAl-Pv were synthesized simultaneously in the same cell. Thus we were able to measure specific volumes (V/V ₀) of Mg-Pv and FeAl-Pv at the same P−T conditions. At all the measurement conditions, values of the specific volume of FeAl-Pv are consistent with those of Mg-Pv within 2 Standard Deviation, strongly suggesting that effect of incorporation of iron and aluminum on the thermoelastic properties of magnesium silicate perovskite is undetectable in this composition, pressure, and temperature range. Two additional runs were performed using a high-pressure cell that has one sample chamber and unit-cell volumes of FeAl-Pv were measured at pressures and temperatures up to 32 GPa and 1,500 K, respectively. All the unit-cell volume data of FeAl-Pv perovskite were fitted to the high temperature Birch–Murnaghan equation of state and a complete set of thermoelastic parameters of this perovskite was determined with an assumption of K′ _300,0 = 4. The determined parameters are K _300,0 = 243(3) GPa, (∂K _T,0/∂T) _P = −0.030(8) GPa/K, a ₀ = 2.78(18) × 10⁻⁵ K⁻¹, and b ₀ = 0.88(28) × 10⁻⁸ K⁻², where a ₀ and b ₀ are the coefficients of the following expression describing the zero-pressure thermal expansion: α _T,0 = a ₀ + b ₀ T. The equation-of-state parameters of FeAl-Pv are in good agreement with those of MgSiO₃ perovskite at the conditions corresponding to the uppermost part of the lower mantle.     

Pressures of Crystallization of Icelandic Magmas

Iceland lies astride the Mid-Atlantic Ridge and was createdby seafloor spreading that began about 55 Ma. The crust is anomalouslythick (20–40 km), indicating higher melt productivityin the underlying mantle compared with normal ridge segmentsas a result of the presence of a mantle plume or upwelling centeredbeneath the northwestern edge of the Vatnajökull ice sheet.Seismic and volcanic activity is concentrated in 50 km wideneovolcanic or rift zones, which mark the subaerial Mid-AtlanticRidge, and in three flank zones. Geodetic and geophysical studiesprovide evidence for magma chambers located over a range ofdepths (1·5–21 km) in the crust, with shallow magmachambers beneath some volcanic centers (Katla, Grimsvötn,Eyjafjallajökull), and both shallow and deep chambers beneathothers (e.g. Krafla and Askja). We have compiled analyses ofbasalt glass with geochemical characteristics indicating crystallizationof ol–plag–cpx from 28 volcanic centers in the Western,Northern and Eastern rift zones as well as from the SouthernFlank Zone. Pressures of crystallization were calculated forthese glasses, and confirm that Icelandic magmas crystallizeover a wide range of pressures (0·001 to 1 GPa), equivalentto depths of 0–35 km. This range partly reflects crystallizationof melts en route to the surface, probably in dikes and conduits,after they leave intracrustal chambers. We find no evidencefor a shallow chamber beneath Katla, which probably indicatesthat the shallow chamber identified in other studies containssilica-rich magma rather than basalt. There is reasonably goodcorrelation between the depths of deep chambers (> 17 km)and geophysical estimates of Moho depth, indicating that magmaponds at the crust–mantle boundary. Shallow chambers (<7·1 km) are located in the upper crust, and probablyform at a level of neutral buoyancy. There are also discretechambers at intermediate depths (11 km beneath the rift zones),and there is strong evidence for cooling and crystallizing magmabodies or pockets throughout the middle and lower crust thatmight resemble a crystal mush. The results suggest that themiddle and lower crust is relatively hot and porous. It is suggestedthat crustal accretion occurs over a range of depths similarto those in recent models for accretionary processes at mid-oceanridges. The presence of multiple stacked chambers and hot, porouscrust suggests that magma evolution is complex and involvespolybaric crystallization, magma mixing, and assimilation. KEY WORDS: Iceland rift zones; cotectic crystallization; pressure; depth; magma chamber; volcanic glass     

Chemistry of rainwaters in the south Pacific area of Russia

On the south-eastern edge of Russia, the chemical composition of rainwater is controlled by sea salts, terrestrial material, as well as volcanic (Kuril islands volcanic area) and anthropogenic emissions, mostly in the southern part of the area. The predominant major ions of the Primorye, Sakhalin and the Kuril Islands rainwaters were respectively HCO₃⁻–SO₄²⁻, Ca–Na, and of Cl–Na. Concentration of trace elements changes within 1–2 orders of magnitude but some difference in the distribution of the elements between continental and island rainwater is found. The concentration of the chemical elements in the particulate fraction varies from < 10% to 90% of the total concentration (dissolved + particulate) with the following distribution: Tl, Na, Ca, Sr, Zn, Cd (< 10%)–Be, Th, Bi, Rb, U, K, Sc (10–20%)–Cu, Mn, Mg, Mo, Se, Ba, Ni, As, Ag, Cs, Co, Y, Ga, V (20–50%)–Sb, Pb, Ge, Cr, Fe, Al (50–90%).The concentration of elements of the particulate fraction of the rainwater usually is significantly different from concentrations in the crust, including both higher and lower concentrations. The terrestrial contribution to dissolved elements was evaluated and follows the decreasing order: Fe > K, Mg, Ca > Ba, Sr > Na (65–1%). Close order was found for total (dissolved and solid) concentrations. Sea salt contribution to dissolved element concentration in the rainwater decrease in the following order: Cl, Mg > K, SO₄ > Ca > HCO₃⁻, Ba, Fe (78–0.1%). Calculation of anthropogenic and volcanic inputs for two ions (Cl⁻ and SO₄²⁻) shows that anthropogenic inputs for the Vladivostok and Yuzno-Sakhalinsk cities can be evaluated as 15–20% of Cl⁻ and up to 80–90% of SO₄²⁻. Volcanic components in the Kuril Islands, where anthropogenic inputs are absent, can reach up to 76% of SO₄²⁻ and 36% of Cl⁻.     

Calculation of reflectance values for two models of texture of carbon materials

This paper presents two models of texture of carbon materials possessing porosity. For both models, R_MAX, R_INT and R_MIN values of reflectance have been calculated by modified Ting's 3A method [Ting, F.T.C., 1978. Petrographic techniques in coal analysis. In: C. Karr (Ed.), Analytical Methods for Coal and Coal Products, Vol. 1. Academic Press, New York, 1978, pp. 3–25]. It has been found that the reflectance indicating surface (RIS) for the studied models of texture changes from a biaxial negative up to a biaxial positive one (the type A texture) depending on the shape and size of pores. In particular case, for the type B texture and the isometry of pores, a uniaxial negative RIS is merely obtained.     

The Equation of State of Lakes

In recent years, a number of workers have studied the stability of deep lakes such as Lake Tanganyika, Lake Baikal and Lake Malawi. In this paper, the methods that can be used to determine the effect that the components of lakes have on the equation of state are examined. The PVT properties of Lakes have been determined by using apparent molal volume data for the major ionic components of the lake. The estimated PVT properties (densities, expansibility and compressibilities) of the lakes are found to be in good agreement with the PVT properties (P) of seawater diluted to the same salinity. This is similar to earlier work that showed that the PVT properties of rivers and estuarine waters could also be estimated from the properties of seawater.The measured densities of Lake Tanganyika were found to be in good agreement (± 2 × 10^-6 g cm^-3) with the values estimated from partial molal properties and the values of seawater at the same total salinity (S_T = 0.568). The increase in the densities of Lake Tanganyika waters increased due to changes in the composition of the waters. The measured increase in the measured density (45 × 10^-6 g cm^-3) is in good agreement (46 × 10^-6 g cm^-3) with the values calculated for the increase in Na⁺, HCO₃ ^-, Mg²⁺, Ca²⁺ and Si(OH)₄.Methods are described that can be used to determine the conductivity salinity of lakes using the equations developed for seawater. By combining these relationships with apparent molal volume data, one can relate the PVT properties of the lake to those of seawater.