散射问题中消除多次透射边界高频振荡失稳措施比较分析

多次透射边界是一种应用广泛的人工边界,类似于其它局部人工边界,它也存在计算稳定性问题。本文分析多次透射边界高频振荡失稳的原因,通过数值试验对已有的消除高频失稳措施进行比较分析,并给出在散射问题中已有消除高频失稳措施的适用性。研究表明,对于完全弹性问题,建议用滤波方法消除高频振荡失稳,但滤波系数的取值需要依据一定的经验;对于黏弹性问题,建议适当控制阻尼取值,利用差分格式本身的高能耗特性来抑制高频振荡失稳,该方法不存在参数选取的问题,更便于应用实施。     

场地地震波动模拟中透射边界稳定性问题研究

本文对场地地震波动模拟中透射边界的稳定性问题进行了一系列研究。对已有的透射边界稳定措施做了比较分析,并给出了不同措施的适用性。结合粘弹性边界的物理模型提出了一个新的消除透射边界飘移失稳的措施,参考高频滤波的思路探讨了一种新的消除透射边界葛频振荡失稳的措施,并用数值试验验证了这两种新的透射边界稳定方法的有效性。将显式行限元一有限差分方法结合透射边界的数值解法的计算结果与解析结果作对比分析,验证这弹一套数值解法的计算精度。利用二维模型模拟了自贡土层和山脊地形的场地效应,利用显戈有限元一有限差分方法结合透射边界这套数值解法进行数值模拟计算,应用汶川地震的实示强震观测记录来验证这套数值解法对于实际地震波动模拟的有效性。     

波动数值模拟中透射边界的稳定实现

从波动能量在计算区内累积增大的观点出发,通过简单的一维弹性波模型,系统地阐明了在近场波动数值模拟中透射边界两类数值失稳--"高频振荡"和"零频飘移"的机理:前者源于对波动数值模拟无意义的高频波动在人工边界上的放大和波在有限计算区内多次反射产生的反复放大;后者则源于透射边界允许零频和接近零频的分量不断进入计算区. 由此提出了稳定实现透射边界的完整方案包括两项简单措施:第一,在全部计算区内按文中建议的方法注入小阻尼,以消除高频振荡;第二, 给出一种具有明确物理意义的消除零频飘移的算子算法. 最后,提供了三维波源问题和散射问题的详细数值试验结果.     

透射边界高频失稳机理及稳定实现——P-SV波动

波动数值模拟的稳定性是获得可靠结果的前提.透射边界是一类具有高阶、高效等特点的人工边界,其引发的高频失稳是由内域格式和透射边界的不当耦合所致.本文针对P-SV波动有限元模拟中透射边界引发的失稳问题,基于GKS定理的群速度解释,通过对有限元和透射边界的频散分析揭示了数值失稳机理为透射边界和相邻内域格式支持了群速度指向内域的高频P波或SV波,波动能量将从边界进入内域引发数值失稳.同时,对比连续模型频散指出引发失稳的谐波是由有限元离散引入.本文采用修改的数值积分方法调整有限元刚度,以消除有限元中引发边界失稳的高频波动成分,从而稳定实现透射边界.理论分析和数值实验均表明本文稳定措施的有效性.     

一种改善多次透射边界稳定性的措施

本提出在边界区设置粘性阻尼的方法，在保持多次透射边界精度的同时可以消除波动数值模拟中的高频振荡失稳现象，中就出平面波动情形分析了采取和未采取这一措施时多次透射边界的稳定性和精度，给出相应的稳定性条件，所得结果表明了这一方法的有效性，中的算例验证了分析结果。     

消除多次透射公式高频振荡失稳的一种措施

多次透射公式（MTF）物理概念简单,便于在计算机上实现时空解藕的高精度波动数值模拟。然而,MTF与其它局部人工边界条件类似,存在数值模拟失稳问题,如高频振荡便是可能出现的失稳现象。本文在分析MTF高频振荡失稳机理的基础上,提出了在波动有限元数值模拟中消除MTF高频振荡失稳的一种措施,即在整个有限元数值模拟区内施加与应变速率成正比的较小粘性阻尼;同时,讨论了这一稳定措施的有效性及其对数值计算精度的影响,并通过数值试验检验了这一稳定措施的可行性。结果表明,消除高频振荡失稳的措施行之有效,且只对波动有限元数值模拟中无意义的高频分量具有抑制作用,而对有意义的较低频段内的波动有限元数值模拟精度影响较小。     

波动的有限元模拟——基本问题和基本研究方法

本文以紧凑形式给出了可以直接用于波动有限元计算的透射人工边界公式。为了阐明离散化误差和由人工边界引起的振荡失稳的基本概念,作者利用简单的一维模型详细分析了波动有限元模拟中遇到的这两个主要问题。基于在离散和连续模型中波动规律差别的分析,首先研究了离散化误差。这一研究导致对可能实现波动有限元模拟的频段的识别,以及采用一种兼备有限差和有限元优点的集中质量离散模型的建议。根据波在人工边界的放大效应和在有限离散模型中多次反射的概念在频域内阐明了这一常见的振荡失稳的机制。在此基础上提出了在保证数值计算精度的条件下通过对人工边界附近节点运动的修正来消除振荡失稳的措施。本文还以一维模型为例给出了人工边界条件的稳定性准则。     

人工边界高频振荡失稳机理的一点注记

进一步完善了人工边界高频振荡失稳机理的解释. 解析证明了对于一维波动模型的离散模型, 当人工波速大于1.5倍的离散空间步距与时间步距的比值时, 则在某一高频段内, 其稳态波动解在人工边界上反射系数的模大于1的命题. 在此基础上对失稳的高频段作了进一步的讨论, 揭示了高频振荡失稳发生于对波动数值模拟无意义的高频段.      

透射边界高频耦合失稳机理及稳定实现——SH波动

就大型近场波动的高效数值模拟而言,稳定实现高阶人工边界是一个尚未圆满解决的问题.本文针对使用多次透射公式的SH波动集中质量有限元模拟,依据GKS定理的群速度解释,进一步阐明了人工边界与内域离散格式耦合所导致高频失稳的机理,即两者支持群速度指向内域的外行高频平面谐波,波动能量自发地从人工边界进入內域,从而导致失稳,而这类谐波是由集中质量有限元离散引入的.本文提出了消除此种耦合失稳的一种方法:通过修改有限元刚度阵来改变内域离散格式,并保证修改格式的精度不低于原有格式的精度.理论分析和数值实验表明此法能稳定实现透射边界.本文研究结果具有推广应用前景.     

多次透射公式实用形式稳定性分析

本文以二维出平面波动方程为例,利用数值方法研究了多次透射公式的实用形式的稳定问题,结果表明：多次透射公式的实用形式与波动方程的集中质量显式有限元结合,不但存在零频漂移,而且同样存在高频振荡失稳,而与中心差分方法结合仅存在零频漂移失稳。     

兰州地震研究所发展的三段历史回顾

对兰州地震研究所的三段历史作回顾,分别为国家地震局、中国地震局兰州地震研究所的体制沿革变迁,中国科学院地球物理研究所;中国科学院兰州地球物理研究所观象台早期创建;兰州地震研究所硕士点早期阶段。     

山西大同数字遥测地震台网记录的震相特征

从地球及地壳的分层结构出发，使用较成熟的理论，从地震波传播的基本原理入手，结合多年的震相观察分析工作实践，选取了山西大同数字遥测地震台网近年来的一些典型震相进行了分析，得出本台网记录到震中距小于100km的地方震震相有Pg，P11，Sg，S11震中距在100km～500km的近震震相有Pn，Pg，P11，Sn，Sg,S11对比分析了在本台网内及周边地区发生的爆破与天然地震在震相上的差异性，并对其特征作了较详细的描述。     

淋巴结结核的影像学表现特点分析

目的:分析淋巴结结核的CT和MRI表现特点,并与病理学对照,探讨两种检查方法的应用价值。方法:回顾性分析2010年1月至2014年12月间34例经手术、穿刺病理或临床证实为淋巴结结核患者的CT或MRI影像学资料,分别从生长方式（融合/非融合）、边缘（光滑/欠光滑）、境界（清晰/欠清晰）、密度及信号、增强后强化方式（均匀/不均匀/环形规则/环形不规则）这5个方面对所受累淋巴结进行分析。结果:34例患者中共累及223枚淋巴结,其中颈部淋巴结67.7%,纵隔内淋巴结27.7%,肺门淋巴1.8%,肺内淋巴结0.5%,腹腔淋巴结0.5%,腹膜后淋巴结1.8%。所有受累淋巴结中,8.5%呈融合性生长,83.9%形态规则,79.4%边缘光滑。CT上发现的171枚病变淋巴结以等密度较多（53.8%）、低密度次之（43.3%）,MRI上发现的66枚病变淋巴结中98.5%为等T1长T2信号。所有CT及MRI观察到的病变淋巴结中,除57%的病灶密度或信号均匀外,仅伴钙化的为12.6%,仅伴部分坏死的26.9%,另外有3.6%同时伴钙化及坏死;增强扫描的204个病变淋巴结中,1%钙化显著、未见明显强化,38.7%均匀强化,伴坏死者内见不强化的低密度区,11.3%呈环形规则强化（均位于颈部）,8.3%环形不规则强化（纵隔内82.4%）,5.9%可见病变淋巴结内有分隔样强化（均位于颈部）。结论:CT与MRI对于淋巴结结核具有重要诊断价值,两者均能清晰显示病变部位、形态以及边缘情况,MRI对于病变边缘及周围情况显示较CT更明确,尤其是增强后,对于包膜及邻近组织受累情况的显示也更加敏感。      

Correction of twilight measurements of influence of night sky

Studia Geophysica et Geodaetica -     

山西夏县中心地震台温泉水化学特征

采集夏县中心地震台温泉水及周边水点样品进行水化学组分特征分析。利用矩形图、Na-K-Mg三角图等方法,初步分析温泉水的水质类型、水-岩平衡状态、热储温度以及循环深度等,并结合氢氧同位素组成特征,初步分析温泉水补给来源。研究表明,夏县中心地震台温泉水化类型属于Na-Cl·SO₄型,水-岩反应属于部分成熟水,热储温度为148.8℃,循环深度为3.91 km,补给源主要为大气降水,温泉为断裂型温泉。以上结果可为该台流体异常分析提供基础研究资料,为夏县地区水化研究奠定基础。     

山东数字地震台网子台背景噪声位移值的分布特征

利用山东数字地震台网的实时记录波形资料,采用在频率域积分方法,计算出各台站背景噪声的位移值。然后进行统计分析,得出了山东地区背景噪声位移值的空间分布和时间变化特征：台站背景噪声幅值与不同地区地质环境、岩性密切相关,位于第四系黄土覆盖层较厚的鲁西聊城、菏泽地区,鲁西北德州地区和鲁北东营、滨州地区各子台的背景噪声位移值较大,鲁中、鲁南山区和胶东半岛地区的基岩台站背景噪声位移值较小;同时,白天噪声值大于夜间噪声值,但深井观测台站的背景噪声值相对稳定。此项研究为绘制台网监测能力图提供基础数据,并为计算各台站的场地响应及震级修正值等测震学研究提供基础资料,促进台网的改建和发展。