藏北北羌塘盆地晚三叠世古风化壳地质地球化学特征及其意义

对最近在藏北羌塘盆地菊花山、肖茶卡、石水河等地区新发现的肖茶卡组顶部的粘土质角砾岩的地质、地球化学特征研究表明:其岩石矿物成分、结构构造、角砾层产状及剖面演化序列等与古风化壳的形成及其特征较为一致;相对于基岩,古风化壳的主量元素表现为CaO、CO2的淋失和SiO2、Al2O3、Fe2O3、K2O、Na2O、MnO、P2O5等不同程度的相对富集;古风化壳中的粘土质风化残积层的稀土元素具有明显的富集,表现为LREE含量远大于HREE的特征,存在明显的Ce的正异常以及Eu的负异常;微量元素除Mo、Sn和U外,粘土质风化残积层的其余元素含量均有明显的增加。这些特征从另一角度证实了北羌塘盆地晚三叠世存在明显的沉积间断。晚三叠世肖茶卡组顶部古风化壳的发现,为了解区域地壳发展史、古气候、古地理演化、层序地层的划分与对比、羌塘含油气盆地分析等提供了重要的依据。     

第四纪以来天津北部燕山隆升与山前盆地的形成*

文中报道了天津市北部燕山山前地带共4个高岩心采取率钻孔的沉积学、地层学、年代学和微体古生物综合研究。420个热退磁样品的磁倾角曲线显示,具147 m厚松散地层的标准孔Z3孔下部地层属于2.58 Ma的松山(Matuyama)极性带。随着地壳沉降,只是从早更新世,冲积扇侧缘泥质沉积物才开始充填了底砾岩层之上的沉积空间,形成广义的山前盆地中的泛滥平原;与传统上视为上新统特征沉积的红棕色黏土类似的亮黄橙-黄橙色-亮红棕色黏土,延续出现到早更新世中期,当时区域沉降中心位于距离山前断裂不足4 km处。0.78 Ma以来的布容(Brunhes)极性带厚57.8 m。谨慎地依据120 ka的布莱克(Blake)亚时确定上更新统后,可见中更新统大部缺失,与天津市区沧县隆起和北京永定河冲积平原南部所见此时的地层间断十分相近。3个钻孔的7个¹⁴C测年数据证实,距离山前断裂不足2 km的钻孔45 ka以来的地层最厚达41 m;显示自MIS3初始期,沉降中心向北迁移更靠近山前断裂,且在构造断裂作用下形成了新的狭义的山前盆地。研究区可见15 ka富有机质黏土及上覆泥炭,表明末次冰消期开始的气候转暖首先影响到华北平原山前地区。此项工作为研究晚新生代以来“源到汇”过程中渤海湾西岸的沉降提供了一个基准点。     

没水倾斜板式防波堤消波性能分析

提出了一种将波浪中倾斜板问题等效化简为波浪中水平板单元组的方法,该方法建立在使用分离变量法求解水工结构边值问题的基础上,并使用伽辽金法精确求解连续边界条件,确定考虑衰减波态的速度势函数,从而求解没水倾斜板结构的消波性能。等效化简法计算精度于边界元法相当,且计算单元数量少、开销低。基于二维线性势波理论,对没水倾斜板式防波堤消波性能分析显示,没水板的倾斜角度、没水深度与板长是结构消波性能的控制因素:没水倾斜板防波堤的消波性能优于没水水平板防波堤的消波性能,随着没水板结构的倾斜角度增大,没水倾斜板结构的波浪透射系数显著减小,且长板优于短板,浅板优于深板;与前人的水槽实验对比显示,相对没水深度与波陡影响结构的消波性能,且波陡造成的波浪破碎贡献了显著的波能消耗。该结论对板式防波堤的结构配置、优化设计有重要意义。     

起伏地表弹性波传播的间断Galerkin有限元数值模拟方法

间断Galerkin有限元法（DG-FEM）作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳定性保持龙格库塔（SSP Runge-Kutta）算法,将DG-FEM推广至时间任意高阶精度.另外,借鉴近最佳匹配层（NPML）的思想,基于复频移（CFS）拉伸坐标变换推导了一种新的PML吸收边界条件（简称为CFS-NPML）,该CFS-NPML能够与DG-FEM算法很好地结合,形成有效的起伏地表地震波传播数值模拟技术.数值试验结果表明,DG-FEM具有高阶精度,可以适应任意复杂起伏地表和复杂构造情况下的弹性波传播数值模拟.同时,CFS-NPML对包括面波等震相的人为边界反射都具有良好的吸收效果.     

间断有限元方法的数值频散分析及其波场模拟

数值求解波动方程是大尺度正演波场模拟、基于波动方程的地震偏移和反演成像的关键.本文针对求解二维声波方程的Runge-Kutta 间断有限元（RKDG）方法的数值频散问题,从理论推导和数值分析的角度进行了深入研究,并将其与近似解析离散化方法（Optimal Nearly Analytic Discrete Method,简称ONAD 方法）、Lax-Wendroff 修正方法、交错网格（Staggered-Grid,简称SG）方法的数值频散进行了比较研究.结果表明：RKDG方法以及近似解析离散化方法在压制数值频散方面要好于上述其他方法,特别是空间精度为3阶的RKDG方法,即使当空间步长取波长的一半,即一个波长内取2个网格点时,最大的频散误差也不超过1.67%.同时,我们也通过波场模拟对比研究了不同数值方法的数值频散问题,进一步直观地验证了数值频散的理论分析结果.     

东北地区古生代地层间断

东北地区古生代在南部有一个地层间断面,北部有三个地层间断面：奥陶纪、晚志留世末-早泥盆世和早石炭世晚期-晚石炭世,相应于古亚洲洋构造域发展的不同阶段,结束于早三叠世末.最后的间断面以后,东北南部和北部古生代地层进入同步发展.     

非均匀弹性介质中旋转交错网格有限差分与任意高阶间断有限元地震波场模拟方法研究

复杂介质中的地震波数值模拟对于地震勘探非常重要.实际应用中经常遇到介质参数剧烈变化的情况,必须选择合适的数值方法进行正演,使波场模拟精度和计算效率满足要求.本文研究了旋转交错网格有限差分法与任意高阶间断有限元法在非均匀弹性介质中波场模拟的精度与计算效率,分析了界面两侧介质参数相对变化量以及界面倾角对上述两种方法数值模拟结果的影响.对于水平界面,界面两侧介质参数一定范围内的改变对旋转交错网格有限差分法的振幅精度和相位精度没有影响;在介质参数存在强反差的情况下,任意高阶间断有限元法需要使用高阶多项式基函数来达到较高的相位精度.有限元法的相位精度优于有限差分法,但需要更多的计算量.对于倾斜界面,当单位波长内含有14个网格点时,界面倾角的变化对旋转交错网格有限差分法的振幅精度及相位精度没有影响,且其精度与任意高阶间断有限元法的精度相接近.     

基于时间可达性与伽萨法则的大卖场区位探讨——以上海市中心城区为例

目前多数基于GIS的商业区位模型与技术研究,未全面考虑交通网络、道路等级、网点吸引力与竞争因素。针对此局限,借助GIS软件,依托可达性方法与伽萨法则,尝试综合考虑上述因素探讨大卖场可达性与商圈特征。以上海中心城区大卖场为例,采用同心圆法、扇形法与最近邻域法探讨卖场空间特征;利用行进成本分析法计算卖场可达性,并根据伽萨法则叠加了卖场引力因素,得到伽萨法则商圈。结果表明:大卖场集中分布于距市中心4.5~10km范围,各方向扩张不均衡。城区整体可达性较好,浦西优于浦东。外环附近商圈面积大,外环附近与市中心区商圈市场规模大,内环附近商圈面积与市场规模小。伽萨法则商圈市场规模"两极分化"既表明卖场市场竞争激烈地区与网点短缺区并存,也说明网点市场规模与职能存在等级差异。     

基于间断位移场的隧道围岩承载能力研究

岩体是岩石块体和分割这些块体的边界的组合体。沿边界面(结构面)的局部剪切变形使得位移场产生间断,对岩体的力学性能产生巨大的影响。采用光滑函数组合来描述滑移面上的间断,得到反映岩体沿着边界局部剪切的数学模型。利用此局部剪切模型得出了应力–滑移关系具有强化阶段时圆形坑道围岩自身承载能力的解析表达式,研究了岩体滑移面对坑道围岩自身承载能力的影响。通过分析得到,自承能力在深部坑道围岩的稳定中发挥了重要作用,应力–滑移关系的强化段对围岩承载力有着很大的影响;围岩压力与坑道的埋深、坑道的尺寸以及围岩的物理力学性质等具有直接关系;围岩出现不稳定的动力现象(岩爆)与此岩体参数m关系极大,而岩体参数与接触面的软化程度、能量转移、耗散结构的几何尺寸以及弹性模量有关。     

砂土液化变形的有限元-无网格耦合方法

通过构造Biot固结理论u-p方程的无网格伽辽金-有限元耦合方法,对砂土液化变形问题进行了数值模拟。对于饱和砂土,采用Oka等提出的弹塑性本构模型,同时采用更新的Lagrange计算格式推导了控制方程。耦合方法能够发挥有限元和无网格各自的优点,既避免了由于单元变形扭曲而引起的计算中断,也可节约计算时间,算例验证了该方法在地震液化问题中的有效性。