二维浅水流动的一种普适的高性能格式——有限体积Osher格式

天然水体具有形状复杂的计算区域和水下地形。本文采用无结构的网格以适合这一情况,并能方便地根据工程应用的要求局部地和适应性地加密网格。相应地,对二维水流计算问题采用了有限体积法的数学表示。跨单元边界的法向数值通量通过用Osher格式求解黎曼问题得到。文中给出了适合于二维浅水方程组的有关单元界面和各种物理边界的法向数值通量公式。分析了这一格式所具有的优点,包括:普适性、守恒性、逆风性、单调保持性、高效性、对间断的高分辨率、边界处理和内部格式相容,以及不引入数值边界条件等。最后,通过长江口南支水流计算的实例,阐明其良好性能。     

利用宽频带远震体波波形研究青藏高原地壳上地幔速度结构的初步结果

利用根据中美合作研究青藏高原深部结构计划布设在青藏高原上的11个宽频带数字地震仪记录到的远震体波数据,采用接收函数(receiver function)反演的方法,对各台站下面地壳上地幔地震波速度结构进行了研究。台站接收函数是通过将三分向地震记录的两个水平分量旋转合成得到径向分量,然后在频率域除以垂直分量并变换回到时间域得到的,它仅与台站下面介质结构有关,而基本上与震源函数和传播路径无关。为压制噪声干扰,对来自同一方向上一定震中距范围内的远震记录得到的接收函数进行了叠加。采用分层弹性介质中弹性波传播矩阵理论,我们可以计算得到分层介质的理论接收函数以及它对各层弹性参数的偏导数,从而利用迭代线性反演可从观测接收函数得到台站下面的一维速度结构。本文给出了其中3个台,即温泉台、格尔木台和日喀则台的初步结果,它们分别位于高原的中部、北部和南部。从各台的接收函数中都可看到清晰的 Moho 面上的 P-S 转换波震相,其相对直达 P 波的走时延迟分别为:温泉台7.9s(东北方向结果),8.3s(东南方向结果);格尔木台8.2s;日喀则台9.0s,如此大的延迟表明高原地壳的巨厚.      

PREM-ZSCHAU滞弹地球模型对表面负荷的脉冲响应

利用PREM模型及ZSCHAU的地球内部粘滞性模型,解算了地球对表面负荷的脉冲响应问题,包括频率域和时间域的解.得到的复勒夫数及复格林函数表明,在ω>10~(-8)/s的频段内,地球的响应可视为基本上是弹性的,只有像冰后调整、地幔对流这样长期的运动,粘滞性才起重要作用;在ω<10~(-13)/s的频段,地球响应呈流体状态.     

南北地震带北段地震二维速度结构—ART方法的分辨率和解的精度

本文详细讨论了使用Pn波走时对南北地震带北段上地幔结构成象ART算法的分辨率和解的精度问题。分析了包含在实际资料中影响解精度的两大因素,给出了具体的解决办法。同时指出在该地区使用Pn波的可分辨尺度极限为25km。在分辨与方差合理折衷情况下,最佳网格化分尺寸为35—55km。计算了解的分辨矩阵和解的协方差矩阵,从而证实了文献[1]的结论。这些结果是初次在该区得到的,是令人满意的。     

概念性河口及长江口淤积机制分析模拟

为了探求河口淤积原因,从概念河口的角度出发,通过概化河口形态建立概念性河口模型,以二维水沙数学模型为工具,模拟分析了5种河口淤积机制;建立了长江口水沙二维数学模型,对长江口区域水沙特征及淤积状况进行了模拟分析。研究表明,建立概念性河口有助于研究不同形式河口的淤积机制,对于长江口水域,其不同区域淤积机制不一,其治理措施也应有所区别,找出河口不同区域起主导作用的泥沙淤积机制,有利于人们更合理地治理河口。     

三维地震勘探在滑动构造解释中的应用

滑动构造影响了煤层赋存特征,改变了煤层厚度和埋深,对煤矿开采不利。以河南某区滑动构造解释为实例,介绍了三维地震技术在精细刻画滑动构造中的解释应用,讨论了滑动面在时间剖面上和三维数据体上的识别特征。结果表明,三维地震勘探是一种有效的滑动构造探测手段。      

欢迎订阅2017年《矿产与地质》(双月刊)

《矿产与地质》创刊于1981年,是经国家新闻出版总署批准,由中国有色金属工业协会主管,中国有色桂林矿产地质研究院有限公司主办,全国公开发行的地质类学术刊物。是中国科技论文统计源期刊(中国科技核心期刊),《中国学术期刊综合评价数据库(CA-