平面P1波斜入射下海底洞室地震响应解析分析

基于无黏性可压缩理想流体介质波动理论和Biot流体饱和多孔介质波动理论,考虑水下饱和土的流固耦合,借助Hankel函数积分变换法(HFITM)给出入射平面P1波在海底洞室周围散射问题的解析解。相比传统研究中的"大圆弧假定",Hankel函数积分变换法可以较好地处理半空间表面边界条件。利用该解析解,计算分析了洞室表面透水条件、入射角度、入射频率、海水水深和饱和土的孔隙率等因素对水-土交界面处水平位移、竖向位移和洞室表面动水压力、环周总应力的影响。结果表明:洞室表面透水条件对水-土交界面处水平位移和竖向位移影响较小;随着斜入射角的增加,水-土交界面处竖向位移减小;随着入射频率的增加,水-土交界面处水平位移随之增加;海水水深为2.5倍SV波的波长时,水-土交界面处水平位移及洞室表面动水压力最大值最大;随着孔隙率的增加,水-土交界面处水平位移、竖向位移和洞室表面环周总应力减小,而洞室表面动水压力随之增加。     

青岛市崂山区宅基地用地现状与对策

以青岛市崂山区3个街道的农村宅基地用地调研为基础,介绍崂山区农村宅基地用地现状,指出当前存在的问题,提出先行解决一部分住房困难村民的宅基地问题、加快编制报批村庄系列规划、强化执法监察、加大土地法律法规的宣传教育力度、积极推进农村建设用地整理等对策和建议。     

Data assimilation using the particle filter for identifying the elasto‐plastic material properties of geomaterials

A computational method, incorporating the finite element model (FEM) into data assimilation using the particle filter, is presented for identifying elasto‐plastic material properties based on sequential measurements under the known changing traction boundary conditions to overcome some difficulties in identifying the parameters for elasto‐plastic problems from which the existing inverse analysis strategies have suffered. A soil–water coupled problem, which uses the elasto‐plastic constitutive model, is dealt with as the geotechnical application. Measured data on the settlement and the pore pressure are obtained from a synthetic FEM computation as the forward problem under the known parameters to be identified for both the element tests and the ground behavior during the embankment construction sequence. Parameter identification for elasto‐plastic problems, such as soil behavior, should be made by considering the measurements of deformation and/or pore pressure step by step from the initial stage of construction and throughout the deformation history under the changing traction boundary conditions because of the embankment or the excavation because the ground behavior is highly dependent on the loading history. Thus, it appears that sequential data assimilation techniques, such as the particle filter, are the preferable tools that can provide estimates of the state variables, that is, deformation, pore pressure, and unknown parameters, for the constitutive model in geotechnical practice. The present paper discusses the priority of the particle filter in its application to initial/boundary value problems for elasto‐plastic materials and demonstrates a couple of numerical examples. Copyright © 2012 John Wiley & Sons, Ltd.     

非连续变形分析方法中的阻尼问题研究

定量研究了非连续变形分析(DDA)方法中的黏性阻尼与数值阻尼。首先,基于Newmark直接积分法,推导了块体系统的运动方程。其次,通过动力学中的黏性阻尼理论建立了DDA中动力系数、时间步长与黏性阻尼比的关系式,探讨了DDA中的常加速度积分方案的数值阻尼分区及阻尼比计算方法,进而得到两种阻尼共同作用时的阻尼比表达式,并分析了频域内阻尼比的分布情况。最后,以谐振激励下的块体振动为例,通过对比不同阻尼作用下块体位移的DDA计算值与理论解,验证了本文提出的阻尼比计算公式的正确性。研究表明：黏性阻尼对低频的衰减作用明显,数值阻尼则可以很快地消除高频干扰,而二者共同作用下可降低阻尼的频率相关性。该研究成果为DDA的振动、波动等动力计算的阻尼取值提供了理论依据。     

河道一维污染源控制反问题

在天然河道一维水质数值模拟中,对初始条件、边界条件及污染源条件分别构造基本解,以此为参数将任一断面的污染物浓度表示为初始条件、边界条件及各污染河段污染物排放量的线性函数.从模拟水域需满足的水环境功能提出约束条件,以污水处理费用最小为目标函数,构造污染源控制问题,并用遗传算法求解.     

东准噶尔巴塔玛依内山组火成岩锆石U-Pb年龄及其地质意义

以东准噶尔双井子地区巴塔玛依内山组和石钱滩组剖面为研究对象,通过岩矿鉴定、地球化学分析、锆石U-Pb定年和牙形刺挑选等方法,对巴塔玛依内山组和石钱滩组的年代学特征进行了研究。结果表明,巴塔玛依内山组火山岩的形成时间始于348. 7±2. 9Ma,318. 7±4. 2Ma之后结束,侵入岩的形成时间为255. 9±3. 5Ma,这些样品中普遍含有捕获锆石,这些锆石具有岩浆锆石的特征,年龄介于早石炭世杜内期和前寒武纪之间(354. 7～971. 2Ma),与盆地周边发育的岛弧火成岩和蛇绿岩年龄相仿,此外在石钱滩组灰岩内得到4颗牙形刺化石,分属于2个属,分别为Idiognathodus和Streptognathodus。以此为基础,通过搜集整理前人数据,得到如下结论:(1)巴塔玛依内山组火山岩地层是穿时的,从早石炭世到晚石炭世连续发育,形成时间介于348. 7±2. 9Ma～300. 4±1. 3Ma之间;(2)巴塔玛依内山组上部的石钱滩组形成于318. 7±4. 2Ma之后的晚石炭世;(3)卡拉麦里蛇绿岩所代表的准噶尔洋盆闭合时间在晚泥盆世法门期和348. 7Ma之间;(4)准噶尔盆地基底可能由岛弧和洋壳组成,且存在前寒武纪基底。本研究旨在确定巴塔玛依内山组和石钱滩组地层时代,并为卡拉麦里洋闭合时限以及盆地基底性质等问题的解决提供依据。     

青海柴达木盆地察尔汗盐湖首采区钾镁盐矿床开发及老卤排放对液体钾矿的影响现状

察尔汗盐湖首采区钾矿床开发利用所产生的水氯镁石——俗称“老卤”,没有合理的排放和开发,使矿床KCl品位降低,其造成的危害是一个长久持续的过程,解决矿床开采过程中造成的矿床老化问题的根本途径是如何将老卤合理开发利用。     

Integrating production data under uncertainty by parallel interacting Markov chains on a reduced dimensional space

In oil industry and subsurface hydrology, geostatistical models are often used to represent the porosity or the permeability field. In history matching of a geostatistical reservoir model, we attempt to find multiple realizations that are conditional to dynamic data and representative of the model uncertainty space. A relevant way to simulate the conditioned realizations is by generating Monte Carlo Markov chains (MCMC). The huge dimensions (number of parameters) of the model and the computational cost of each iteration are two important pitfalls for the use of MCMC. In practice, we have to stop the chain far before it has browsed the whole support of the posterior probability density function. Furthermore, as the relationship between the production data and the random field is highly nonlinear, the posterior can be strongly multimodal and the chain may stay stuck in one of the modes. In this work, we propose a methodology to enhance the sampling properties of classical single MCMC in history matching. We first show how to reduce the dimension of the problem by using a truncated Karhunen–Loève expansion of the random field of interest and assess the number of components to be kept. Then, we show how we can improve the mixing properties of MCMC, without increasing the global computational cost, by using parallel interacting Markov Chains. Finally, we show the encouraging results obtained when applying the method to a synthetic history matching case.     

基于线性互补的非连续变形分析

非连续变形分析(DDA)方法是一种新的用来分析块体系统运动和变形的非连续介质数值计算方法。研究的核心工作是致力于对现有DDA接触问题处理方法的改进。DDA主要采用罚函数法和Lagrange乘子法处理接触问题,合理设定罚参数很困难,此外,因开闭迭代而引起的刚度矩阵的不连续变化也会导致收敛方面的困难。为避免引入罚参数及传统意义上的开闭迭代,用混合线性互补模型(LCDDA)对DDA方法进行了重新描述。在此基础上,综合基于非光滑分析的Newton法的局部平方收敛和最速下降法的全局线性收敛的优势,提出求解LCDDA模型的有效算法。根据上述思想及理论研究成果编制了完整的计算程序,算例计算结果证明了方法的精度及可行性。     

岩溶红土地基中的软土类型及其工程特性

根据岩溶红土地基软土的形成条件划分出了岩沟型、层间型及洞中型三种不同的软土类型,在此基础上,通过取样测试并结合前人已有成果,对不同类型软土的工程地质特性进行了分析.研究结果表明:(1)岩溶红土地基软土是一种与岩溶作用关系密切的特殊软土;(2)岩溶红土地基软土的天然含水量和压缩性均明显小而重度又普遍高于泻湖相软土;(3)60%的岩沟型残积软土其天然含水量均小于液限含水量;(4)岩溶红土地基软土的工程地质问题主要表现为地基不均匀沉陷、边坡失稳和诱发土洞塌陷等.基于岩溶红土地基软土形成条件及其工程地质性质比较特殊,建议把它叫作"岩溶软土",以示与其它成因软土相区别.