孔隙介质的渗透特性初探

孔隙的存在是岩土类介质材料结构的本质特征,它不但改变了岩土体的力学特性,而且严重影响着岩土体的渗透特性。大多数经典的渗流理论中,多孔介质模型都假定孔隙率和渗透系数是与时间无关的材料常数。实际上由于淘涮、侵蚀、冲刷等原因,它们是随时间和坐标变化的,同时又与孔隙中的压力、流速等因素有关。基于孔隙率和损伤变量之间的定量关系,从连续损伤力学的角度对多孔介质岩土材料的渗流力学特性进行了研究。首先,对传统的达西定律形式进行修正,提出了孔隙介质完备有效的达西定律（模型）;然后,对该模型中渗透参数的特性进行了讨论和分析,得出了一些有益的结论。     

覆水场地地震反应分析

随着各种海洋结构物的兴建,覆水场地的地震反应逐渐成为研究热点。基于任意拉格朗日-欧拉描述,推导时变区域上的流体运动方程,给出流场、结构的接触条件和流场网格运动控制方法。对于平坦覆水场地,水平向地震动作用下,根据Couette流理论证明该类场地流体作用可以忽略;竖向地震动激励下,横向均匀场地可以通过动水压力公式准确考虑流体作用,横向非均匀场地则需要通过流固耦合方法考虑流体作用,以海底隧道为例加以说明。对于起伏场地,天然起伏场地在地震动激励下的动力反应具有明显的流固耦合特征,以三角形起伏场地为算例;结构物的兴建造成的人工起伏场地同样需要考虑流固耦合效应,以某沉管隧道在水平向地震动激励下的动力反应为算例,并根据结果初步提出该类结构物流固耦合分析的简化计算方法。     

基于ITRF2005的全球板块运动模型

利用最新的国际地球自转服务(IERS)发表的国际地球参考框架ITRF2005的速度场,建立了一个完全基于现代空间大地测量实测结果的现今全球板块运动模型ITRF2005VEL,并与地质模型NNR-NUVELlA及以往利用ITRF序列所建立的模型进行了比较.结果表明,从整体上看ITRF2005VEL与NNR-NUVELlA模型有较好的一致性,且较以往的ITRF序列的模型精度更高,但同时也存在着一定的差异性,这与观测台站、板块本身等相关因素都有关.本文还由ITRF2005VEL模型计算了构造板块的总角动量,结果不为零,说明ITRF2005框架有可能不满足无整体旋转的要求.     

欧拉反褶积中构造指数的计算与优化选取

构造指数是欧拉反演中最重要的参数,以往选取构造指数往往根据地质先验信息和人工经验试探选取,这样往往不能快速,准确地确定反演构造指数,从而影响反演效果,甚至反演失败.本文提出利用异常值自动估算构造指数,这样避免了经验选取构造指数产生的不确定因素的干扰,可以快速方便地选出较合理的构造指数,提高反演效率.模型实验显示了该方法良好的反演效果,证明了文中提出的方法的可行性.     

基于测绘期刊论文数据的可视分析

结合文献计量方法与可视化技术,以《测绘学报》《测绘通报》《测绘科学》《测绘工程》2003-2013年共收集的5863篇学术论文为数据源,进行了关键词、论文作者及论文作者所在地3方面的统计分析,并将统计结果以统计图表和地图的形式展现,直观地反映了我国测绘学科研究热点的演变,以及活跃作者群的人数、核心研究区的地理分布。研究显示,涉及大地测量、地理信息系统和地图制图3方面研究内容的论文占总论文数的59.9%,达一半以上;根据普赖斯定律确定的活跃作者群人数为131人;发文量排名前10的城市中,有90%为省会城市或直辖市。从研究结果可以得出,我国测绘学科的发展每年的研究热点都不相同,其中大地测量、地理信息系统和地图制图的研究较为广泛;测绘学科发展具有高质量的活跃作者群体,但学科的发展具有地域不均衡性。     

断层演化过程中的周期和非周期地震滑移及非地震滑移一维断层模型分析

本文首先根据Dieterich和Ruina提出的含速率和状态的摩擦定律(Dieterich-Ruina定律), 基于一维弹簧-滑块模型推导了地震复发周期的解析表达式, 然后将该近似解与数值模拟结果以及Barbot等的相关研究进行了对比分析． 此外, 本文还利用数值模拟与理论分析研究了断层周期和非周期演化的力学成因机制以及非地震滑移形成的另类力学机制, 并讨论了一维弹簧-滑块模型的优点及其局限性． 结果表明: ① 震后滑移和自加速/成核阶段的持续时间在整个演化过程中不能被忽略; ② 在修正后的复发周期模型中, 复发周期的长短除了与断层特征尺度、 作用于断层面上的有效正应力和远场加载速率相关外, 还受Dieterich-Ruina定律中摩擦参数的取值以及临界滑移距离的影响; ③ 当给定各个物理参数和几何参数时, 目前所得到的解析近似解可以很好地估计地震的复发周期, 其相对误差可小于5%; ④ 在断层演化过程中, 施加剪切应力加载会产生非周期的地震滑移, 而在自加速/成核阶段后期或震后滑移阶段早期, 施加较大的剪切应力加载, 则会出现非地震滑移．      

基于欧拉反褶积方法估算华北地区重力变化的等效源参数

利用华北地区2009—2014年绝对重力与相对重力多期重复观测资料, 得到不同时空尺度的华北区域重力场动态变化图像。 采用欧拉反褶积方法, 通过对理论模型试算, 获得最优反演参数, 对引起华北地区重力场变化的场源深度和空间分布规律进行了反演计算和分析。 结果表明: 当构造指数为1时, 适合对流动重力场变化数据进行反演。 实际资料的反演结果在2009—2014年间的场源位置集中于河套断裂带。 本文研究方法可以用于重力场反演和以场求源的定量研究。 本工作为流动重力变化信号的反演和解释提供一条新思路, 也为构建地震重力预报指标体系提供了定量依据。     

断层自发破裂传播的数值模拟:速率状态摩擦定律本构参数的影响

地震往往受控于滑动面的摩擦性质,这种摩擦性质可以由速率状态摩擦定律较好地描述.速率状态摩擦定律中的本构参数a和b与动态摩擦系数相关,从而影响着同震位移与剪切应力的时空演化.本文在前人工作的基础上,采用三维边界积分方程法模拟速率状态摩擦定律控制下均匀全空间中平面断层的自发破裂传播过程,并详细讨论了a和b对滑动速率、剪切应力和破裂传播速度的影响.数值结果表明a和b的不同取值将导致不同的破裂行为,b－a的值越大,断层越不稳定,这种不稳定性有利于裂纹的产生与扩展.但滑动速率的时空分布不只依赖b－a,而且还与a和b的具体取值有关,断层面上滑动速率峰值与剪切破裂强度均随着a的减小而增大,随着b的增大而增大.相关结果有助于加深对断层自发破裂传播的认识.     

基于垂向一阶导数与解析信号比值的欧拉反演方法

常规欧拉反褶积法中构造指数的选取以及分散解存在较多的问题,提出了基于联立垂向一阶导数与解析信号的欧拉齐次方程的RDAS-Euler反演方法。该方法可以更为精确的估计场源的范围及埋深,且不需考虑构造指数N的影响,避免了因构造指数不当而引起的反演误差。通过对单一地质体及组合地质体模型的实验证明本文方法能有效地完成目标体的反演工作,反演结果与理论值之间的误差小于10%,且相对于常规欧拉反褶积法更加稳定准确,能够更好的得到地质体边界及深度信息。将RDAS-Euler法应用于黑龙江省虎林盆地实测布格重力异常数据,获得了丰富的断裂信息,说明RDASEuler法增强了对断裂平面位置的识别能力。     

风沙运动的欧拉双流体模型参数研究

利用数值方法研究了双流体模型中动力黏度μ_s,kin、动量交换系数β和碰撞恢复系数e_ss对风沙运动的影响.结果显示:动力黏度μ_s,kin影响沙床面起沙率,动量交换系数β影响输沙率沿流程的变化,碰撞恢复系数e_ss影响输沙量的垂向分布.利用研究得到的参数,模拟风沙流运动,并与已有的诸多实验和数值模拟结果对比,结果吻合较好,说明该参数的选取可以有效地模拟风沙流运动.