起伏地表条件下各向异性地震波最短路径射线追踪

在地震波正反演研究中,考虑起伏地表和地震各向异性具有非常重要的理论意义和实际应用价值.本文在前人研究的基础上,将最短路径追踪算法引入到起伏地表各向异性介质模型的地震波走时计算中.模型剖分时,整体模型划分成正方形单元,起伏边界附近以不规则网格逼近,进而采用非规则节点布置实现非规则网格处的最短路径计算.追踪计算中采用Sena群速度近似公式,得到各向异性地震波的走时,实现了复杂地表情况下各向异性介质模型中地震波的射线追踪.理论模型计算结果显示,本文方法能够可靠地应用于复杂各向异性介质模型,具有较高的计算精度.     

SIRT法层析成像在某核电站岩溶区的应用

本文阐述了联合迭代法的基本原理,通过4种典型速度模型的数值模拟实验,验证该方法的有效性。以最短路径法做正演,联合迭代法做反演对某核电站岩溶区的井间声波CT探测资料进行成像,最后给出溶洞可能存在的位置。     

巩乃斯避暑旅游气象指数分析与评价

根据巩乃斯、和静站2009—2018年6—8月的逐日平均气温、最高气温、相对湿度、风速等气象观测资料,采用对比分析、通径分析及模糊综合评价等方法对巩乃斯地区夏季避暑旅游气象资源优势进行分析和评价.结果表明:巩乃斯夏季气候凉爽、湿度适宜、风速不大,人体感觉凉爽、舒适,适合避暑旅行.人体舒适度与空气温度呈显著正相关,与相对...     

青藏高原及周边2001年以来三次特大地震引起的形变场分布特征

基于地壳黏弹模型在GPS观测资料和地震位错数据为约束条件下,应用有限元数值模拟方法对2001年昆仑山8.1级、2008年汶川8.0级和2015年尼泊尔8.1级特大地震引起的地壳形变分布特征进行了模拟计算,获得了地震位移场和形变场.这3次大地震分别发生在青藏高原的北部、东部边界和南部边界,尽管震级大小基本相当,但发震区域和发震断层性质都各不相同,其各自产生的地壳形变场分布特征存在明显的差异,主要表现为形变场区域大小、幅值大小等的差异,以及在不同地壳深度也存在明显的差异,这些差异主要取决于地震断层性质和地下介质结构的不同.通过模拟计算,可以进一步了解大地震产生的应力加卸载区分布特征,对预测未来地震发生区域范围提供重要参考依据.     

汶川地震科学钻探3号井孔附近微震震源参数确定及其意义

基于微震监测仪器频带及波形记录的时频特征,我们筛选出汶川地震科学钻探3号井孔周围微震台阵2012年记录的218个t_s-t_p < 1 s的微震.通过盖戈法与和达法相结合确定微震震源的几何参数,发现这些微震分布大体呈NE-SW展布,与龙门山断裂带的走向基本一致.无论是采用《地震台站观测规范》中的量规函数（量规函数GF）还是李学政等（2003）的量规函数（量规函数LXZ）,近震震级均与矩震级呈现出较好的线性关系,向震级小的一端延伸时都表现为M_L < M_W,但采用李学政等的量规函数时该趋势更加明显.同时采用Brune和Boatwright震源谱衰减模型对观测震源谱的拟合表明拐角频率具有模型依赖性：基于Brune模型拟合得到的拐角频率大于基于Boatwright模型得到的拐角频率.无论基于哪一种模型,矩震级与拐角频率、破裂半径的对数线性关系均较弱,与应力降、视应力的对数则存在较好的线性关系.这些关系不支持应力降、折合能量、视应力为常数的观点,表明微震与大地震的震源物理过程存在差异.视应力与应力降成比例,比例系数小于0.5,表明破裂动力学模式符合Savage-Wood模式.近震震级与矩震级拟合关系M_w=a+bM_L中b的大小与应力降和地震矩的关系有关,Δσ∝M₀^γ,则b=1/（1+γ）,因此从b的大小可以粗略地判断应力降与地震矩的关系.本文对应于γ≈1的情况,与基于量规函数LXZ得到的近震震级与矩震级关系中系数b=0.53吻合,这说明仅从辐射能量的角度考虑量规函数LXZ较量规函数GF准确.

     

地震的应变张量观测与应用前景

地震发生时的动态应变场,在研究地震触发、地震破裂、地面破坏、水文和岩浆变化等方面都具有重要应用意义.地震的应变张量观测和现有的惯性地震仪观测的物理量不同.前者可以直接记录到地震发生时震源辐射的应变（应力）波,而后者记录到的是位移、速度或加速度.地震频率的应变测量在地震学中的应用前景主要表现在：①测量震源机制解理论预言的辐射4象限分布;②测量库仑应力变化;③换算成动态应力以评估地震烈度;④测量地震波的能量密度;⑤测量地震断层形变加速和形变局部化过程.用惯性地震仪的记录虽然在理论上也可以解算出动态应变值,然而种种原因导致计算结果的误差很大,往往不可接受.应变张量地震仪若能与现有的惯性地震仪配套起来,形成大规模台阵,则有可能推动应变地震学的诞生,在地震观测和地震学科领域引起重大革新.

     

钢管桩加固软土地基的数值模拟分析

针对广乐高速公路先行段,分别建立天然软土地基和钢管桩处理地基的两种平面应变有限元模型。通过施工阶段分析,分别计算得出在两种条件下地基在路堤施工过程中每阶段的位移场和应力场。通过对比分析揭示钢管桩加固软黏土地基的机理。结果表明软土地基经过桩处理后的整体性得到改善,桩通过承担很大的上部荷载减小了软弱土层的受力进而减小地基变形。此外,路堤中心处土层的应力和变形均较大,从经济角度考虑,可在路堤中心小范围内加密桩,两侧减少桩。     

华北克拉通辽北清原地体新太古代基性麻粒岩变质作用演化

对华北克拉通新太古代的构造演化模式有多种不同的认识, 需要进行更加深入的变质作用研究。通过对辽北清原地体基性麻粒岩进行系统的岩相学观察、矿物化学分析、相平衡模拟和锆石定年研究, 以阐明其变质演化过程和大地构造意义。研究选择的基性麻粒岩样品分为含石榴石域(19DJ07-GD)和不含石榴石域(19DJ07-NGD)两类, 含石榴石的区域呈条带状且分布不均匀。两种区域都发育两期麻粒岩相组合。在含石榴石域, 第一期变质矿物组合为石榴石+单斜辉石+斜方辉石+角闪石+黑云母+斜长石+石英。其中, 第一期斜长石(Pl1)发育复杂成分环带, 钙长石摩尔分数(xAn)从核部到幔部升高, 然后再向边部降低; 第一期角闪石(Amp1)的Ti成分环带同样为从核部到幔部升高再向边部降低。通过矿物组合和相应的成分环带推测第一期麻粒岩相变质作用具有逆时针型P-T轨迹, 包含峰期前升温升压阶段以及峰后降温降压阶段。通过相平衡模拟约束峰期温压条件为0.8~0.9 GPa/900~950 ℃, 达到高温—超高温(high-ultrahigh temperature)变质条件。锆石定年结果表明变质作用峰后冷却时间为2498±6.9 Ma(MSWD=0.39)。综合区域上的"穹隆-龙骨"构造、逆时针的变质轨迹以及和TTG岩浆活动晚期脉冲几乎一致的表壳岩变质时间, 表壳岩超高温麻粒岩相变质作用被认为受太古宙特有的垂向构造/沉落(sagduction)构造体制控制。第二期变质组合以局部生长的石榴石+石英±单斜辉石的后成合晶/冠状体为特征, 代表一期与古元古代造山事件有关的高压麻粒岩相变质作用。     

新时期榆林煤矿区生态保护修复与综合治理策略及路径探索

榆林作为陕北地区国家能源基地的重要组成部分,在过去30余年为中国经济社会发展作出了重要贡献。然而,煤炭资源开发不可避免地给当地生态环境带来了采空区塌陷、土地与水资源破坏、生态退化等诸多问题,且目前因煤炭规模开发的负面效应持续加剧着矿区水资源、土地资源的破坏以及生态系统韧性的下降,已然威胁到榆林经济社会的可持续发展及国家能源资源的长远安全。多年来,榆林煤矿区在地质环境治理方面形成了一批可复制、可推广的经验和模式,但目前仍然呈现出整体效果不显著、基金使用率低、实施阻力大、系统设计受限、顶层规划缺乏等“五个问题”的困境,多元化治理模式和相应的政策机制也尚未形成合力。为解决这一实际,本研究在剖析煤矿区生态问题现状及生态修复工作面临问题及原因的基础上,结合煤矿区治理现状及乡村高质量发展需求,创新形成了“1个目标、2个底线、3个空间、4种模式”的榆林市煤矿区综合治理系统性解决策略,建议了3项创新机制以及“三步走战略”的具体实践路径,总结了适用于榆林煤矿区生态修复+N的综合治理模式清单,最后从政策、技术、管理等方面提出了对策建议。研究结果有助于为下一步系统性、整体性指导榆林市各煤矿企业开展矿山生态保护修复工作提供借鉴,切实促进煤矿区生态保护修复走向多元化综合治理模式、走向社会–经济–生态协调的绿色发展新格局。     

断层地震孕育和发生的不稳定性模型

在以前的工作中,考虑直立走滑型断层地震,假设断层面微元破裂强度遵循Weibull概率分布,由细观力学方法推导出断层面的宏观本构关系是一个非线性函数,表现为弹性-软化塑性特征,在此基础上用稳定性理论研究了地震稳定性问题.而实际断层大多是倾斜的,为此,本文首先建立了由围岩和倾斜断层构成的平面地震力学模型,采用宏观的断层载荷-变形的全过程曲线,详细讨论了倾斜断层地震的不稳定性问题.结果表明,远场一旦施加位移,断层也同时错动,这可能与实际情况不符合.为了更好的模拟断层的初始能量累计过程,进一步对断层本构模型进行改进.考虑断层面破裂强度,采用Coulomb破裂准则,则断层表现为刚塑性本构关系,只有当断层面剪应力达到一个临界值时,断层才开始错动.研究表明,对于倾斜断层地震,与直立走滑型断层地震一样,系统刚度比β（围岩切线刚度与断层刚度最大值之比）是决定地震失稳的重要参数,只有当ββ≥1时,仅仅是断层无震滑动,不会发生地震.在远场应施以位移形式边界条件,以致地震失稳发生在平衡路径的位移转向点并伴有应力突跳.