MODFLOW和FEFLOW在国内地下水数值模拟中的应用

数值模拟是研究地下水资源问题的主要手段。介绍MODFLOW和FEFLOW这两款地下水数值模型在国内地下水模拟中的应用,比较分析二者在基本原理、模型特点、以及应用领域等方面的差异,提出地下水数值模型的发展趋势。     

延边八家子镇上南沟“大理岩”地质地球化学特征及地质意义

传统观点认为,延边八家子镇上南沟"大理岩"为碳酸盐岩地层经过区域变质作用形成。通过详细的野外和镜下观察,发现该"大理岩"并不发育层理等典型沉积岩的特征,而主要呈岩株和岩墙产出,与围岩新东村岩组片岩、片麻岩呈港湾状侵入式接触;且其中发育大量围岩捕虏体,"大理岩"局部还发育有冷凝边;除上南沟"大理岩"外,在上南沟及周边地区并没有碳酸盐岩地层的发育。分析显示,上南沟"大理岩"具有与火成碳酸岩相似的微量元素和稀土元素组成特征,C、O同位素组成处于火成碳酸岩受低温蚀变的演化趋势之内。因此,初步认为延边八家子上南沟"大理岩"并不是由碳酸盐岩地层经区域变质作用形成,其成因可能存在以下2种解释:一种是上南沟"大理岩"为火成碳酸岩,很可能是华北克拉通破坏事件在北缘的响应;另一种可能是大规模富含CO2流体运移的产物,是大规模流体运移成矿的典型范例。     

东太平洋深海沉积物小于2μm组分的稀土元素地球化学特征

用等离子光谱测定了东太平洋９个深海沉积物样品的小于２μｍ组分及全岩的稀土元素,结果显示,小于２μｍ组分及全岩的稀土元素组成均受控于沉积物类型,而二者之间稀土元素组成的差异在钙质和硅质沉积物中出现相反的趋势。元素相分析和稀土元素赋存状态分析表明,这种差异主要由于沉积物物质组成的不同所致,表面在化学组成上,主要与ＭｎＯ,Ａｌ２Ｏ３Ｆｅ２Ｏ３,ＣａＯ含量的变化有关,大多数样品在小于２μｍ组分的稀土组成     

Late Jurassic-Early Cretaceous Northern Qaidam Basin, NW China: Implications for the earliest Cretaceous intracontinental tectonism

Formation of Mesozoic western China, which was dominated by tectonic amalgamation along its southern margin and associated intracontinental tectonisms, holds a key for interpreting the succedent Cenozoic evolution. This paper presents new data including lithology, sedimentary facies, stratigraphic contact, seismic interpretation and paleo-structures within the Upper Jurassic-Lower Cretaceous strata in the northern Qaidam Basin, NW China. These data all account for a contractional tectonic deformation in the earliest Cretaceous. The South Qilian Shan, according to the sedimentary features and provenance analysis, reactivated and exhumated during the deformation, controlling the deposition of the Lower Cretaceous sequences. A simplified model for the Late Jurassic-Early Cretaceous paleogeography and tectonics of the northern Qaidam Basin is accordingly proposed. The results also support a ∼25° clockwise rotation of the Qaidam Basin since the Early Cretaceous and a more accurate Mesozoic evolution process for the basin. This earliest Cretaceous deformation, associated with the reactivation of the South Qilian Shan at the time, are part of the intracontinental tectonisms in central Asia during the Mesozoic, and probably driven by both the closure of the Mongol-Okhostk Ocean to the north and the collision of the Lhasa and the Qiangtang blocks to the south.     

末次盛冰期东亚气候的成因检测

在国际古气候模拟比较计划设置的标准试验方案下,首先利用中国科学院大气物理研究所的全球大气环流模式（IAP-AGCM）模拟了末次盛冰期东亚气候状况,然后通过4组数值敏感性试验逐一模拟了大气CO2浓度、海洋表面温度（SST）和海冰、陆地冰盖和地形、东亚植被变化4项强迫因子的单独气候效应,进而对末次盛冰期东亚气候的成因进行了检测。结果表明,末次盛冰期除华南局部略有升温外,中国年均地表气温显著降低,降温幅度总体上向北增大,青藏高原处存在一个降温中心。其中,SST和海冰变化是华南局部略偏暖的主因,它同时导致了东亚其他区域地表气温的显著降低,特别是在东北亚地区;陆地冰盖和地形变化对于东亚地表气温的显著冷却作用主要体现在东亚的西北部;大气CO2浓度降低会引起东亚地区0.2～0.9℃的普遍降温;相对而言,东亚植被的降温作用（0.5～1.0℃）主要显现在中国40°N以南的区域。与此同时,SST和海冰变化能引起中国东部年均降水一定程度的减少,而大气CO2浓度、陆地冰盖和地形、东亚植被单独变化均不会显著影响东亚年均降水的分布状况,然而,上述四项因子的共同变化会通过协同作用引起中国东部年均降水的显著减少,西部地区降水则与现在差别不大。此外,末次盛冰期东亚夏季风的显著减弱源于SST和海冰变化,冬季风变化则可归因于SST和海冰、陆地冰盖和地形的变化。     

流体饱和两相多孔介质动力反应计算分析

基于流体饱和两相多孔介质的弹性波动方程组,运用显式逐步积分格式与局部透射人工边界相结合的时域显式有限元方法对该波动方程组进行求解,对两相介质在输入地震波作用下的弹性动力反应进行计算和分析;对在是否考虑孔隙流体渗流的两种情况下计算得到的两相介质弹性动力反应结果的差异进行对比研究,从而揭示孔隙流体渗流对两相介质动力反应性质的影响。计算结果表明：两相介质弹性动力反应时程的波形与入射地震波的波形相同,且弹性动力反应的峰值出现的时刻对应于入射地震波的峰值出现的时刻;孔隙流体的渗流将对两相介质的弹性动力反应性质产生显著的影响。数值计算同时表明,时域显式有限元方法是进行流体饱和两相多孔介质弹性动力反应计算分析的一种有效的方法。     

大跨度输水洞穿高速浅埋暗挖方案模拟比选研究

利用非线性有限元程序,对三跨拱涵开挖施工方案进行模拟。通过对两种较为可行的施工方案进行模拟,比较了开挖引起的地表沉降、衬砌的应力变化和洞周土体的塑性区等因素,综合考虑后利用模糊评判选出了优选方案。同时,由模拟的数据还发现了隧洞开挖时容易出现的各种危险情况,建设施工单位在危险区域采取适当的措施,得到很高的实际价值,对今后类似工程有借鉴意义。     

Pumping decisions for sustainable development of groundwater resources in areas of grassland degradation: a case study in Lanqi Banner, Inner Mongolia, China

The best planting alternatives for satisfying high water use demands of forage and fodder crops in a region of Inner Mongolia, China, were determined by a multiobjective distributed-parameter groundwater management model. These alternatives took account of different cropping patterns and pumping decisions associated with both temporal and spatial aspects of water allocation. The model was developed for phreatic, homogenous, and isotropic aquifers using the response matrix technique of quadratic programming theory and, in this case, using the alternative direction implicit (ADI) scheme. Model solutions using effective rainfall with a probability of 50%, show that average water table drawdown in the planning period (2006–2017) is 0.22 m and the groundwater fluctuation in each pumping well is very low. In order to evaluate the pumping decisions under an effective rainfall with a probability of 75%, a sensitive analysis was also conducted. Analysis shows that it is useful to apply the results from the proposed model to control the landscape degradation due to overgrazing and overpumping activities.     

土地开发整理项目规划设计变更的原因及对策

通过对临沂市部分县区土地开发整理项目规划设计变更实际情况的数据分析,对土地开发整理项目规划设计变更的类型进行了分类介绍,分析土地开发整理项目规划设计变更的主要原因是规划设计与项目区的实际情况脱节,并在此基础上提出了严格控制土地开发整理规划设计变更的对策。     

Simple models for dynamic melting in an upwelling heterogeneous mantle column: Analytical solutions

Several lines of evidence suggest that the melt generation and segregation regions of the mantle are heterogeneous, consisting of chemically and lithologically distinct domains of variable size and dimension. Partial melting of such heterogeneous mantle source regions gives rise to a diverse range of basaltic magmas. In order to better assess the role of source heterogeneity during mantle melting, we have undertaken a theoretical study of trace element distribution and fractionation during concurrent melting and melt migration in an upwelling, chemically heterogeneous, two-porosity double lithology melting column. Analytical solutions for the abundance of a trace element in the matrix and channel were obtained under the assumptions that the porosity, melt and solid velocities, and solid-melt partition coefficients are constant and uniform. For simplicity, we neglected diffusion and dispersion in the melt. Chemical source heterogeneities of arbitrary size and shape were integrated into the simple melting models by allowing trace element abundance in the source region to vary as a function of time and space. Concurrent melting and melt migration in an upwelling heterogeneous mantle may be approximated as a quasi-steady state problem in which time-dependent concentration patterns produced by melting of heterogeneous source regions are superimposed on a reference steady-state concentration distribution established by melting of the ambient or background mantle. Chromatographic fractionation is especially important for the matrix melt and solid when chemical heterogeneities are involved during melting and melt migration in the mantle, giving rise to significant phase-shift between two incompatible trace elements in the matrix melt and scattered correlations among incompatible trace elements in residual peridotites. Mixing is the chief mass transfer process in the dunite channel where the chromatographic effect is negligible for most of the incompatible trace elements. The lack of chromatographic fractionation among incompatible trace elements and isotopic ratios in MORB suggests either most MORB are channel melts or mixing in magma conduit and chamber is very efficient such that the phase-shift is averaged out during magma transport and storage processes. Advection brings melt produced by smaller-degree of melting in the deeper part of the melting column to the overlying melting region, increasing the incompatible trace element abundance in the matrix and the channel. This advection-induced self-enrichment is especially important when heterogeneous sources are involved and may account for some of the enriched incompatible trace element patterns observed in residual peridotite that were previously interpreted to be a result of mantle metasomatism. Systematic studies of high-resolution spatially correlated mantle samples may help to constrain the melting history and the size and nature of chemical heterogeneities in the mantle.