25～75℃酸性NaCl溶液中方铅矿的溶解动力学

在25～75℃、pH=0.43～2.45的1mol/LNaCl溶液中进行了方铅矿的溶解动力学实验。发现在远平衡条件下,方铅矿的溶解速率r与氢离子活度犤H+犦呈线性关系,溶解速率方程(速率定律)为:r=k犤H+犦,即对H+而言,溶解反应为一级。其中速率常数k为2.344×10-7mol/m2·s(25℃)、1.380×10-6mol/m2·s(50℃)、7.079×10-6mol/m2·s(75℃)。溶解反应的活化能为43.54kJ/mol,方铅矿的溶解机理为表面化学反应,速率决定步骤为表面配合物的离解。     

某盆地油气勘查成本预算和风险分析

某高原盆地的油气资源量较为丰富，可找到大型油气田。根据含油气盆地地质构造特征、生油岩的厚度与热演化程度、储集层和区域性盖层的发育情况等因素，确定了两个含油气远景区，拟作为重点找油气靶区。提出了两个可供选择的油气勘探布署方案。在此基础上，选用适当的成本预算基本参数，进行油气勘查成本估算，并对两个方案进行了优选对比与具体的风险分析，提出方案决策意见。     

Dissolved silica budget in the North basin of Lake Lugano

We studied the dissolved silica cycle in the water column of the North basin of Lake Lugano, Switzerland/Italy. Lake Lugano is a meromictic, eutrophic lake, permanently stratified below 100-m depth. A one-box model was used to calculate a silica mass-balance over 1993, based on various lake measurements, such as sediment traps, sediment cores, water analysis and biota countings. We found that the North basin of Lake Lugano is at steady state as far as dissolved silica is concerned. The primary source of dissolved silica in the lake is river input (about 80%), with diffusion from bottom sediments and groundwater input also playing a role. Atmospheric input is negligible. The main export of dissolved silica occurs via biogenic uptake by diatoms and final burial of their frustules in the bottom sediment. Loss of dissolved silica through the lake outflow only represents 15% of the total output. Of the total amount of Si exported to the lake bottom through diatom sinking, less than 20% is re-supplied to the surface water by diffusion. Thus, the lake acts as an important permanent sink for silica. The long residence time of dissolved silica in the lake (7 years) is related to the strong physical stratification of the lake. Only about 10% of the standing stock are available to phytoplankton uptake.     

有限弹性变形地壳中的地震波

地壳由半无限大的基岩上一层厚度为H^-的表土层组成，入射地震波为垂直的SH波，产生水平地面运动。当浅源大地震发生时，在极震区以外行波传播产生地面运动将使地壳介质有非线性的有限弹性变形。用小参数摄动法使非线性控制方程为线性化的小参数各阶控制方程，得出头两阶线性控制方程的解析解。     

Riparian alder fens — source or sink for nutrients and dissolved organic carbon? — 2. Major sources and sinks

Natural riparian forest wetlands are known to be effective in their ability to remove nitrate by denitrification and sediments with attached phosphorus via sedimentation. On the other hand, litter input and decomposition is a process of crucial importance in cycling of nitrogen and phosphorus in a forest ecosystem.In this study we investigated the amount of nitrogen and phosphorus entering the alder fen ecosystem through leaf litter and its decomposition and the removal capacity of nitrogen and phosphorus by measuring denitrification and sedimentation in the alder fen.We found an average input of leaf litter during fall 1998 of 226 g m⁻² yr⁻¹ DW with nutrient concentration of 0.17% P and 1.6% N. This means a yearly input of 0.4 g m⁻² yr⁻¹ P and 3.6 g m⁻² yr⁻¹ N. The decomposition of leaf litter using litter bags with small and large mesh size resulted in bags with macroinvertebrates (large mesh size) and without macroinvertebrates (small mesh size). After 57 days the litter bags with macroinvertebrates had a decomposition rate of 79%.Denitrification was measured in May and June of 1997 using the acetylene inhibition technique on intact soil cores and slurry-experiments. The average annual denitrification rate was 0.2 g m⁻² yr⁻¹ N using data from the core experiments. The denitrification rate was higher after addition of nitrate, indicating that denitrification in the riparian alder fen is mainly controlled by nitrate supply.The sedimentation rate in the investigated alder fen ranged from 0.47 kg m⁻² yr⁻¹ DW to 4.46 kg m⁻² yr⁻¹ DW in 1998 depending on the study site and method we used. Sedimentation rates were lower in newly designed plate traps than in cylinder traps. The alder fen also showed lower rates than the adjacent creek Briese. Average phosphorus removal rate was 0.33 g m⁻² yr⁻¹ P.Input sources for the surface water of the alder fen are sediment mineralization and decomposition of leaf litter; output sources are sedimentation and denitrification. This study showed that a nutrient input of 24.58 kg ha⁻¹ yr⁻¹ N, 8.8 kg ha⁻¹ yr⁻¹ P and 419 kg ha⁻¹ yr⁻¹ DOC into the surface water of the alder fen is possible. Alder fens cannot improve water quality of an adjacent river system. This is only true for a nearly pristine alder fen with the hydrology of 10 months flooded conditions and 2 months non-flooding conditions a year.     

剥蚀及地幔作用下青藏高原隆升过程的数值模拟

修改了England和Mckenzie的黏性薄层流变模型中控制大陆形变的连续性方程,将剥蚀作用对高原隆升演化的影响直接引入该方程,并考虑下伏地幔小尺度对流对增厚岩石层的搬离作用对高原隆升演化后期的影响,用有限差分法直接模拟青藏高原隆升过程. 数值模拟结果所显示的高原隆升演化过程与实际观测资料吻合较好,揭示了高原隆升演化过程的非平稳和多阶段的特性;同时还表明上地幔小尺度对流对岩石层底部的搬离作用可能是最近8Ma以来高原快速隆升的主导机制.     

基于波动方程的广义屏叠前深度偏移

地震波传播算子的计算效率和精度是制约三维叠前深度偏移的关键因素. 广义屏传播算子(GSP, Generalized Screen Propagator)是一种在双域中实现的广角单程波传播算子. 这一方法略去了在非均匀体之间发生的交混回响,但它可以正确处理包括聚焦、衍射、折射和干涉在内的各种多次前向散射现象. 通过背景速度下的相移和扰动速度下的陡倾角校正,广义屏算子能够适应地层速度的强烈横向变化. 这种算子可以直接应用于炮集叠前偏移,通过将广义屏算子作用于双平方根方程,还可以获得一种高效率、高精度的炮检距域叠前深度偏移方法,用于二维共炮检距道集和三维共方位角道集的深度域成像. 本文首先简述了炮检距域广义屏传播算子的理论,进而讨论了共照射角成像(CAI, Common Angle Imaging)条件,由此给出各个不同照射角(炮检距射线参数)下的成像结果,进而得到共照射角像集. 由于照射角和炮检距的对应关系,共照射角像集又为偏移速度分析和AVO(振幅随炮检距变化)分析等提供了有力工具.     

Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.     

Predicting average annual groundwater levels from climatic variables: an empirical model

On the basis of one-dimensional theoretical water flow model, we demonstrate that the groundwater level variation follows a pattern similar to recharge fluctuation, with a time delay that depends on the characteristics of aquifer, recharge pattern as well as the distance between the recharge and observation locations. On the basis of a water budget model and the groundwater flow model, we propose an empirical model that links climatic variables to groundwater level. The empirical model is tested using a partial data set from historical records of water levels from more than 80 wells in a monitoring network for the carbonate rock aquifer, southern Manitoba, Canada. The testing results show that the predicted groundwater levels are very close to the observed ones in most cases. The overall average correlation coefficient between the predicted and observed water levels is 0.92. This proposed empirical statistical model could be used to predict variations in groundwater level in response to different climate scenarios in a climate change impact assessment.     

贡嘎山东坡亚高山森林区蒸散力的估算

以海螺沟 30 0 0m气象站的观测资料为基础 ,运用Penman公式法　空气饱和差法和桑斯维特公式法 ,计算了贡嘎山东坡亚高山林林的年平均蒸散力 ,其结果为分别为 431 81mm、171 4mm和 44 6 4mm。通过分析蒸散力的影响因素 ,对这三种计算蒸散力的方法作了比较 ,对计算结果存在的差异作了比较合理的解释 ,认为用Penman公式法可以估计出本研究区的蒸散力。同时 ,对蒸散力及其影响因自进行了相关分析 ,指出湿度和风不是速制约研究区蒸散力的主导因子 ,并分析了蒸散力与水面蒸发的关系 ,由此推导出估算蒸散力的简便方程 :PE =6 6 77 0 6 91E60 1 0 75 1ITm 0 0 396Pm。