基于螺旋线上谱因式分解的地震波场隐式辛算法

均匀介质、复杂各向同性介质和各向异性介质中的地震波传播过程,可用统一形式的标量声波方程描述.考虑到在无损耗条件下,地震波方程描述了地震波场这一个无穷维的哈密顿体系随时间的演化过程,该过程为一个单参数连续辛变换,因而可以在其哈密顿形式表述下导出其辛格式.与显式辛算法相比,隐式辛格式对应的隐式辛几何算法具有无条件稳定的特点,可以允许较大的计算步长.但是由于隐式算法不可避免地面临高阶矩阵的求逆,其每一步的计算速度较慢.为实现矩阵快速求逆,文中采用了螺旋边界条件下谱因式分解的方法.在螺旋边界条件下,需要求逆的矩阵化为带状矩阵,而且其各列非零元素的位置和大小具有非常好的相似性,因而可以采用谱因式分解的方法实现快速LU分解.文中采用二阶精度的隐式蛙跳辛格式和谱因式分解方法,计算了常速度、层状介质和Marmousi模型中的波场.计算表明,隐式辛算法不失为波场计算的一种好方法.     

Experimental and computational validation and verification of the Stokes–Darcy and continuum pipe flow models for karst aquifers with dual porosity structure

In our previous study, we developed the Stokes–Darcy (SD) model was developed for flow in a karst aquifer with a conduit bedded in matrix, and the Beavers–Joseph (BJ) condition was used to describe the matrix–conduit interface. We also studied the mathematical well‐posedness of a coupled continuum pipe flow (CCPF) model as well as convergence rates of its finite element approximation. In this study, to compare the SD model with the CCPF model, we used numerical analyses to validate finite element discretisation methods for the two models. Using computational experiments, simulation codes implementing the finite element discretisations are then verified. Further model validation studies are based on the results of laboratory experiments. Comparing the results of computer simulations and experiments, we concluded that the SD model with the Beavers–Joseph interface condition is a valid model for conduit–matrix systems. On the other hand, the CCPF model with the value of the exchange parameter chosen within the range suggested in the literature perhaps does not result in good agreement with experimental observations. We then examined the sensitivity of the CCPF model with respect to the exchange parameter, concluding that, as has previously been noted, the model is highly sensitive for small values of the exchange parameter. However, for larger values, the model becomes less sensitive and, more important, also produces results that are in better agreement with experimental observations. This suggests that the CCPF model may also produce accurate simulation results, if one chooses larger values of the exchange parameter than those suggested in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

Microbialite of anoxic condition from Permian-Triassic transition in Guizhou,China

Spherical microfossils are present in the Hindeodus parvus zone of the Lower Triassic in Ziyun,Guizhou Province. They generally range from 0.15 to 0.30 mm across,with micritic wall and filled by micro-sparry calcites,and are evenly scattered in micritic matrix. Their abundance makes the rock as-signed to microbialites. The accompanied organisms include ostracods and algal mat,but no gastro-pods or bivalves. Presence of small (<7 μm) pyrite framboids indicates that this bed formed in anoxic conditions. In some sections,this bed is overlain and underlain by tidal-flat micritic limestone with microgastropods and small burrows. Occurrence only in deposits on reef top indicates that the micro-bial organism was benthic,and needed sunlight in life. The size of the microbial fossil exceeds that of any bacteria or cyanobacteria. Thus,it does not belong to these two phyla. It may belong to lower green algae,and is assigned to a new species in a new genus,Ziyunosphaeridium sinensis gen. et sp. nov. Up to now,the rocks described as microbialites from the Permian-Triassic transition include six types: (1) porous micritic limestone such as that descried from Japan by Sano and Nakashima (1997),(2) limestone with rich globular microfossils such as that described from the Nanpanjiang Basin in China by Lehrmann (1999),(3) dendroidal limestone such as that described from the vicinity of Chongqing by Kershaw et al. (1999),(4) spherical microbial limestone adapted to anoxic environments described from Ziyun,Guizhou reported in this paper,(5) limestone with cyanobacterial fossils such as that described by Wang et al. (2005),and (6) stromatolites. All these microbialites are not reefs in the proper sense,and the argument that Permian reef ecosystems extended into the Mesozoic is incorrect.     

A Methodology for Clustering Lakes in Alberta on the basis of Water Quality Parameters

In this study, a methodology for clustering 18 lakes in Alberta, Canada using the data of 19 water quality parameters for a period of 11 years (1988–2002) is presented. The methods consist of (i) principal component analysis (PCA) to determine the dominant water quality parameters, (ii) cluster analysis techniques to develop the characteristics of the clusters, and (iii) pattern‐match lakes to determine the appropriate cluster for each of the lakes. The PCA revealed that three principal components (PCs) were able to explain ～88% of the variability and the dominant water quality parameters were total dissolved solids, total phosphorus, and chlorophyll‐a. We obtained five clusters for the period 1994–1997 by using the dominant parameters with water quality deteriorating as the cluster number increased from 1 to 5. Upon matching cluster patterns with the entire dataset, it was observed that some of the lakes belonged to the same cluster all the time (e.g., cluster 1 for lakes Elkwater, Gregg, and Jarvis; cluster 3 for Sturgeon; cluster 4 for Moonshine; and cluster 5 for Saskatoon), while others changed with time. This methodology could be applied in other regions of the world to identify the most suitable source waters and prioritize their management. It could be helpful to analyze the natural controlling processes, pollution types, impact of seasonal changes and overall quality of source waters. This methodology could be used for monitoring water bodies in a cost effective and efficient way by sampling only less number of dominant parameters instead of using a large set of parameters.     

地下水机器学习方法研究

随着我国地下水监测工作的高速发展,高频率高密度水位监测数据的出现催生了对其进行深入信息挖掘的需求。在传统地下水模型研究中,地下水水位监测值常位于模型构建过程的下游,当水位监测的时空密度逐渐增大时,新增信息无法有效传导至模型的规划阶段并指导概念模型的修订。文章提出了一种地下水系统补排边界的识别方法,在不建立地下水数值模型的前提下,以监测井空间位置为节点,按照德劳内原则建立三角网格。在此网格系统中,首先定义一个水力梯度变换函数gradF,以求取网格中任意位置的水力梯度;借鉴机器学习领域的优化算法,使用水力梯度场驱动含水层中随机分布质点的运行轨迹,并以此推断和识别区域内地下水补给和排泄边界。在环境地学计算平台EnviFusion-CGS中实现,并构建了详细工作流程。以山东省青岛市大沽河中下游含水层为示范区,对含水系统的补给区和排泄区的空间分布及其动态变化进行了分析,取得了良好效果。本研究为构建和修订已有含水层概念模型提供了新思路。     

The Upper Devonian Kellwasser event recorded in a regressive sequence from inner shelf to lagoonal pond,Catalan Coastal Ranges,Spain

This study, conducted in the Catalan Coastal Ranges, north‐east Spain, describes the Upper Devonian Kellwasser event in a shallowing‐upward sequence of black shales, siltstones and quartz arenites. This sequence was deposited in a progradational and regressive coastal system where the sedimentary environment evolved from the inner shelf to a lagoonal pond located landward of the shoreline. Three anomalous succeeding steps have been identified by geochemical analysis. The first one, detected on the inner shelf, was characterized by oxygen depletion and high organic productivity. The second, detected in the nearshore, was caused by hydrothermal activity occurring under normal oxic conditions. The third and most intense step was identified in the muds of the lagoonal pond and has been linked to strong anoxic conditions, elevated clastic input derived from changes in the weathering regime at the source area and moderate hydrothermal activity. The Kellwasser event is thus defined in the study area as stepwise and multi‐causal. This is the first time that the Kellwasser event has been identified in a sedimentary environment behind the shoreline. It is also the first time that it has been reported in the Catalan Coastal Ranges.     

钻孔瞬变电磁方法探测越界开采采空区的应用

煤矿井下越界开采扰乱煤炭资源开发秩序,同时存在巨大的安全隐患,准确探测越界开采采空区及其边界对安全采掘极为重要。常规的矿井瞬变电磁探测方法容易受到巷道中的金属和电磁信号干扰,探测精度有限,难以满足采空区边界精准解释的需求。为了解决这一问题,提出一种钻孔中收发信号的钻孔瞬变电磁方法和施工工艺,利用其远离巷道和靠近异常体的特点,提高异常信号在实测数据中的占比;并研究可行域约束的OCCAM反演算法,进一步减小体积效应影响,实现了对越界开采采空区的电阻率精细成像。结合山西某煤矿的工程实践,对该方法精准探测邻矿越界开采采空区规模、性质的实用性和有效性进行了检验。结果表明,钻孔瞬变电磁方法是矿井地球物理技术与钻探工艺技术的高效结合和有益补充,基于可行域约束的OCCAM反演算法能够有效应用于越界开采采空区的边界精确成像解释。      

关于暴时电离层电流分布的南北半球不对称性

采用国际上广泛认可的高层大气和电离层经验模式提供的各种参数, 通过电离层电流连续方程, 计算出强磁暴条件下6月至日和12月至日内, 磁纬±72°和磁地方时00:00～24:00之间电离层电场、电流等的分布. 计算中考虑了地磁和地理坐标间的偏离; 除中性风场感生的发电机效应外, 还包含了磁层耦合(极盖区边界的晨昏电场和二区场向电流)的驱动外源. 结果表明, 6月至日时, 磁层扰动自极光区向中低纬的穿透情况在南、北半球内基本接近, 北半球内略强; 但12月至日时, 呈现明显的不对称性, 南半球的电流穿透远强于北半球, 而电场的穿透则是在北半球更强. 无论南北半球, 在中高纬地区, 午夜至黎明时段出现较强的东向电场分量, 其[WTHX]E×B[WTBZ]的向上漂移效应, 正是解释我们以往不少研究现象中所期盼的物理机制.     

THERMODYNAMIC PROPERTIES AND ROLE OF THERMALS IN THE STRATUS-TOPPED BOUNDARY LAYER

From results of two large-eddy simulation of stratus-topped boundary layer,the structure,thermodynamic prop-erties and role of thermals are investigated by using conditional sampling methods,which divided the thermals into thewarm/moist,cool/dry,warm/dry and cool/moist events.The results show that the main turbulent circulation in thestratus-topped boundary layer is composed of the warm/moist updraft and cool/dry downdraft.Below entrainmentregion,the warm/moist updrafts and cool/dry downdrafts are,respectively,positively and negatively buoyant andcontribute most to total fluxes and variances.Evaporative cooling has important effect on the structure and thermodynamic properties of thermals instratus-topped boundary layer.     

用变分法修正黑河试验数据的初步研究

为了对黑河试验观测资料进行订正,利用一个一维边界层垂直扩散模式和以拟牛顿极小化算法为基础的最优化方法建立了一个变分同化系统,其中,目标函数相对于控制变量的梯度用扰动法求取.并利用模式资料和黑河试验若干时刻的观测资料进行了一些数值试验.试验证明,将初始场变量与模式误差系数同时作为控制变量加以调整,能够有效调整初始场观测资料.并且,用变分同化方法能够较线性插值(时间、空间)方法更有效、合理的修正黑河试验观测资料场的种种缺陷.