Numerical modeling of reservoirs or pipes in groundwater seepage

In groundwater studies, the numerical modeling of complex boundary conditions can be made easier by considering reservoirs and pipes as “reservoir” finite elements, which can store or release large volumes of water at almost constant hydraulic head. The numerical code to be used must solve the complete conservation equation with unlinked functions for the water retention curve and the unsaturated hydraulic conductivity, and the ability to describe these as step functions. Four examples illustrate the performance of the “reservoir” element: reservoir pumping, laboratory variable-head permeability test, vertical flow in an open borehole, and pumping test with well storage.     

长偶极大功率可控源电磁波响应特征研究

地球物理学中关于电磁波勘探研究通常采用的是地球半空间模型.然而,对于几十公里的有限长电缆源（长偶极源）,远距离电磁波场探测必须要考虑电离层的影响,它是一个全空间问题.关于包含电离层、空气层和地球介质（我们称“地-电离层”模式）的电磁波场特征的研究在国外较少,国内几乎是空白.本文采用全空间积分方程法首先对小尺度的可控源电磁波场特征进行了研究,由于此时电离层的影响可忽略,它应该和半空间成熟的CSAMT模拟结果一致,对比结果表明,二者是一致的,验证了全空间模拟方法的可靠性和有效性.随后进行了50 km长电缆电离层和空气层高度都为100 km的“地-电离层”模式大尺度电磁波场模拟,以探讨大尺度可控源电磁波场的特征.给定频率的“地-电离层”模式电磁场的衰减曲线表明长电缆远距离电磁波场由于受电离层的作用存在衰减逐渐变小的过渡场和衰减变小的波导场.为了探讨复杂介质“地-电离层”模式电磁波特征,对“地-电离层”模式的典型地盾和地台多层介质模型进行了数值模拟,得到了偶极源长度50 km、电流200A、收发距离远达1600 km和2500 km的合理的电磁场结果.最后,对一简单含油储层结构模型进行了长偶极、大功率、远距离电磁波场响应计算.储层横向不均匀复杂结构模拟的结果表明,考虑电离层和大气层的“地-电离层”模式大尺度深层复杂介质模拟时,电磁场对深部目标体仍有很好的异常响应.     

Grampian migmatites in the Buchan Block,NE Scotland

Rocks exposed along the Scottish coast between Fraserburgh and Inzie Head contain information critical to understanding the evolution of the Buchan Block, the type locality for low‐P, high‐T regional metamorphism, and its relationship with the rest of the Grampian terrane, one of the major tectonostratigraphic components of the Scottish Caledonides. The ~8 km long section traverses a regional network of shear zones and, at the highest grades around Inzie Head, passes into the core of the Buchan Anticline, a large‐scale open fold that is commonly regarded as a late structure, post‐dating metamorphism. The metasedimentary rocks increase in grade from upper amphibolite to granulite facies and preserve unequivocal evidence for partial melting. The diatexite migmatites around Inzie Head, along with other gneissose units within the Buchan Block, have been regarded as allochthonous Precambrian basement rocks that were thrust into their current position during the Grampian orogenesis. However, field observations show that the onset of in situ partial melting in metapelitic rocks, which was associated with the formation of garnet‐bearing aplites and associated pegmatites, occurred around Fraserburgh, where shear fabrics are absent. Thus, the rocks preserve a continuous metamorphic field gradient that straddles the shear zone network. This observation supports an alternative interpretation that anatexis was the result of mid‐Ordovician (Grampian) metamorphism, rather than an older tectonothermal event, and that the Inzie Head gneisses are autochthonous. Using an average mid‐Dalradian pelite as a plausible representative protolith, phase equilibria modelling satisfactorily reproduces the observed appearance and disappearance of key minerals providing that peritectic garnet produced with the first formed melts (represented by the garnet‐bearing aplites) depleted the source rocks in Mn. The modelled metamorphic field gradient records a temperature increase of at least 150 °C (from ~650 °C near Fraserburgh to in excess of 800 °C at Inzie Head) but is isobaric at pressures of 2.7–2.8 kbar, suggesting the Buchan Anticline developed synchronous with partial melting. The Buchan Anticline is likely an expression of crustal thinning and asthenospheric upwelling, which produced voluminous gabbroic intrusions that supplied the heat for Buchan metamorphism.     

山东半岛典型冷涡暴雪个例对流云及风场特征的观测与模拟

利用卫星云图、 多普勒天气雷达资料、 常规观测资料和MM5模式模拟结果, 分析了2005年12月4日典型冷涡造成山东半岛暴雪的云带变化、 对流单体的生成与发展、 低空急流特征及雷达强回波带的变化特征。结果表明： 冬季冷涡携带强冷空气南下, 经过渤海暖湿下垫面后影响山东半岛时, 卫星云图上从渤海到山东半岛可以清楚地分析出造成暴雪的西南-东北向结构特征的横向云带; 雷达图上, 在渤海到山东半岛形成一强降雪回波带, PPI径向速度图上可清楚地分析出暴雪区的低空急流和中尺度对流的低空辐合, 雷达径向速度垂直剖面反映暴雪产生在零速度线附近和低空急流的存在; MM5模式模拟结果显示： 暴雪发生时山东半岛在1.5 km存在低空急流, 水平风的时间-高度垂直剖面反映暴雪发生时干冷空气经过暖湿的渤海海面在700 hPa以下出现暖平流, 垂直上升运动和散度的低层辐合、 高层辐散与雷达图像反映的强降雪回波是相对应的。     

对流层上传重力波的非线性演化

利用二维全隐欧拉格式对重力波在可压、非等温大气中的非线性传播过程进行了数值模拟和分析.分析结果表明,从对流层顶激发的重力波能稳定地经平流层传到中层顶,从而将能量和动量从一个区域带到另一个区域;在向上传播过程中,重力波经历了发展、位温翻转、对流直至最终破碎的演变;重力波的破碎是对流和小尺度波动的重要的源,对流不稳定和翻转是非线性现象的一个基本特征.计算还显示,扰动源的大小直接影响着重力波的非线性传播过程,当扰动源足够小时,重力波能稳定传播,而大振幅扰动可以加速重力波的破碎.     

Clarifying the distal to proximal tephrochronology of the Millennium (B–Tm) eruption,Changbaishan Volcano,northeast China

Tephra dispersed during the Millennium eruption (ME), Changbaishan Volcano, NE China provides one of the key stratigraphic links between Asia and Greenland for the synchronization of palaeoenvironmental records. However, controversy surrounds proximal-distal tephra correlations because (a) the proposed proximal correlatives of the distal ME tephra (i.e. B–Tm) lack an unequivocal chronostratigraphic context, and (b) the ME tephra deposits have not been chemically characterized for a full spectrum of element using grain-specific techniques. Here we present grain-specific glass chemistry, including for the first time, single grain trace element data, for a composite proximal sequence and a distal tephra from Lake Kushu, northern Japan (ca. 1100 km away from Changbaishan). We demonstrate a robust proximal-distal correlation and that the Kushu tephra is chemically associated with the ME/B–Tm. We propose that three of the proximal pyroclastic fall units were erupted as part of the ME. The radiocarbon chronology of the Kushu sedimentary record has been utilised to generate a Bayesian age-depth model, providing an age for the Kushu tephra which is consistent with high resolution ages determined for the eruption and therefore supports our geochemical correlation. Two further Bayesian age-depth models were independently constructed each incorporating one of two ice-core derived ages for the B–Tm tephra, providing Bayesian modelled ages of 933–949 and 944–947 cal AD (95.4%) for the Kushu tephra. The high resolution ice-core tephra ages imported into the deposition models help test and ultimately constrain the radiocarbon chronology in this interval of the Lake Kushu sedimentary record. The observed geochemical diversity between proximal and distal ME tephra deposits clearly evidences the interaction of two compositionally distinct magma batches during this caldera forming eruption.     

数值试验中时间步长的选取

建立一个二维的海陆风静力模式来讨论时间步长的选取，发现定时间步长不宜取得太大或者太小，必须是某个适当的值。变时间步长的模拟效果要比定步长理想。时间步长越大，海风中心强度模拟值变小。     

祁连山云系云微物理结构和人工增雨催化个例模拟研究

人工影响天气的学科基础是中小尺度天气动力学与云降水物理学,需要将天气-动力-云降水物理耦合为一体。考虑到目前将天气动力学性质的基础数值模式用于人工影响天气中的问题,从数值模式动力方程、模式分辨率、云物理过程、数值求解方案、初边值条件等方面系统地探索了发展人工影响天气数值模式中一些需要重点解决、且不可忽视的特色问题,并举例对相关问题提出了解决思路和方法。期望提出的问题有助于构思更适合于人工影响天气数值模式,使数值模式功能真正向满足人工影响天气的要求靠近一步。     

Numerical study of droplet size dependent chemistry in oceanic,wintertime stratus cloud at southern mid-latitudes

Cloud droplet chemistry is modelled for the first 150 m of rise in a wintertime, mid-latitude, marine stratus cloud using observations made at and near the Cape Grim Baseline Station as a source of input parameters. The emphasis in this work was to study the variation in droplet chemistry as a function of both droplet size and nucleus composition, with a particular focus on the way in which oxidation of dissolved sulfur dioxide varied.At 150 m above the condensation level, solute concentration as a function of droplet size was found to increase by as much as 2 to 3 orders of magnitude for only a factor of 2 increase in droplet radius, primarily as a consequence of the 1/r dependence in the droplet growth equation. This type of size dependence exists at all levels in the model cloud, and has a significant influence on oxidation rate of sulfur dioxide in droplets growing on sulfate nuclei, oxidation by ozone being favoured in the smallest droplets, but oxidation by hydrogen peroxide being favoured in the larger droplets. Oxidation by ozone is favoured at all sizes in droplets formed on sea-salt nuclei as a result of the initially high alkalinity of these droplets, and in the cloud overall is calculated to be the more important oxidation pathway. Although based on a simplified chemical scheme, these results suggest that both size-dependent and nucleus-dependent chemistry of cloud droplets may need to be considered explicitly in cloud modelling work.Volume-weighted mean pH values in the range 5 to 6 were predicted from sensitivity studies in which input variables were varied over reasonable ranges, in agreement with two sets of bulk cloud-water pH data obtained by aircraft near Cape Grim.     

水汽方程积分计算格式的数值试验

该文把Easter(1993)的水汽方程积分方案应用到球坐标的有限区模式（LAFS）中,并作了两组比较性试验。Easter的方案是把水汽的通量方程和连续方程联立求解。在计算中不会出现水汽负值。而且稳定性好。本文所做的两组试验表明。Easter的方案比LAFS的原方案所产生的计算性扩散小。因而降水预报更接近实况。但CPU时间增加了16%。