地壳流变结构控制作用下的龙门山断裂带地震发生机理

青藏高原东缘低地形变速率的龙门山断裂带上相继发生了2008汶川M_w7.9级地震和2013芦山M_w6.6级地震.地震勘探与震源定位结果揭示了龙门山区域地震空间分布特征：纵向上,龙门山断裂带这两次地震主震均发生在龙门山断裂带上地壳的底部（14~19 km）,绝大部分余震均发生在上地壳范围（5~25 km）,而在其中、下地壳深度范围内鲜见余震发生;横向上,地震（M_w>3）在龙门山断裂带青藏高原一侧密集分布且曾有大震发生,而四川盆地地震稀少（M_w>3）.为探讨龙门山断裂带地震发生机理,并解释以上龙门山区域地震空间分布特征,本文建立了龙门山断裂带西南段跨芦山地震震中区域的四种不同流变结构的龙门山断裂带三维岩石圈模型,以地表GPS观测资料为约束边界条件,数值模拟龙门山断裂带岩石圈在数千年以上长期匀速构造挤压作用下的应力积累特征,探讨了地壳分层流变性质对地壳应力积累的影响,分析了该区域地震空间分布与构造应力积累速率的关系.计算结果表明：该区域在数千年的应力积累过程中,脆性上地壳中应力表现近于恒定值的线性增长趋势,龙门山断裂带上地壳底部出现应力集中积累现象,这一应力集中现象可以解释龙门山断裂带汶川地震与芦山地震主震的发生,及其大部分余震在脆性上地壳中的触发;青藏高原一侧上地壳应力积累速率远远高于四川盆地的应力积累速率,这一应力积累分布现象可以解释龙门山区域青藏高原一侧地震密集而四川盆地地震稀少的地震空间分布特征;通过比较不同流变结构模型中的应力积累状态,认为导致这一应力积累空间分布状态的重要控制因素在于青藏高原中、下地壳较低的黏滞系数与四川盆地中、下地壳较高的黏滞系数的差异.在柔性的中、下地壳内,应力增长近于指数形式,稳定状态之后其应力增长速率近于零,构造应力积累难以达到岩石破裂强度,因而鲜见地震发生.地壳各层位的应力增长率差异与地震成层分布的现象共同揭示了龙门山区域岩石圈分层流变结构：脆性上地壳、韧性中、下地壳（青藏高原一侧较弱,四川盆地一侧较强）、韧性岩石圈上地幔.     

四分量钻孔应变观测的实地相对标定

提出钻孔应变观测的实地标定问题有相对标定和绝对标定两层涵义;讨论了四分量钻孔应变观测的相对标定的归一化新方法：将前人标定应力观测的方法用于应变观测。一方面大大简化了计算公式,另一方面提出轮流以每个元件为标准计算灵敏度而后取平均值的方法。在统计上更为合理;并以昌平台RZB-1型钻孔应变观测为例进行实际应用,取得了令人满意的结果。     

冷泉形成的数值模拟研究

海底冷泉形成的一种可能机制是海平面下降引起天然气水合物的分解.本文基于对冷泉渗漏特征的分析,建立了二维轴对称模型,利用有限元方法定量分析了南海区域海平面下降对冷泉形成的影响.结果表明,末次冰盛期(26.5~19.0ka BP)海平面下降引起的冷泉活动可以持续到现在,但是从水合物停止分解至今,超孔隙压力的极值在持续减小,而流体向海底的渗漏达西速度先快速减小、然后缓慢减小.同时发现,流体向海底的渗漏达西速度与管状通道的渗透率、通道周围介质的渗透率以及通道的半径等有关,估计目前的冷泉活动还可以持续10000年以上.海平面下降引起的天然气水合物分解,可能是影响全球气候变化的一个重要因素.     

2015年4月25日尼泊尔M_S8.1大地震的同震效应

应用有限单元方法,计算了2015年尼泊尔MS8.1大地震发生产生的同震变形和应力变化.计算中考虑地球为球体以确保远场应力场变化得到可靠结果,采用PREM模型的地球分层模型,考虑了中国地震局(CEA)和美国地质调查局(USGS)各自提供的断层滑动模型.结果表明:尼泊尔MS8.1地震是一个比较典型的低角度逆冲地震,水平位移和应力降较大;地震造成南北方向上的水平位移最突出,且集中在首都加德满都附近区域.USGS断层滑动模型地表最大位移量达到3.5m,CEA滑动模型最大为1.2m;东西向和垂直方向上的同震位移相对较小;同震位移量级在0.1m的影响区域可达300km;地震造成尼泊尔地区最大库仑应力变化可达到MPa量级,地震危险性依然较大.此次MS8.1地震对我国西藏地区有一定影响,特别是雅鲁藏布江地区和拉萨块体南北走向的正断层,库仑应力变化为正,量级可达数千帕乃至十余千帕,应该注意该区被诱发中强震的可能性.     

用遗传有限单元法反演川滇下地壳流动对上地壳的拖曳作用

利用震源机制解和地质调查资料,运用伪三维遗传有限单元法反演了中国川滇部分地区（96°E～104°E, 22°N～30°N）受到的边界作用和该地区底部所受的剪切作用力. 对反演方法进行了讨论,获得了稳定的反演结果. 结果显示,该区川滇菱形块体上地壳底部受下地壳南南东向剪切力,与GPS反映的现今地块运动方向大体一致;而研究区域其他地区底部没有受到统一显著的剪切力作用. 结合该地区的研究资料,初步认为青藏高原物质受挤压向东和东南运动过程中,下地壳物质比上地壳更易于流动,从而对川滇菱形块体上地壳有拖曳作用. 从应力场反演的模型位移与GPS实测的现今位移资料也大体吻合,反映结果有合理的物理意义.     

各向同性膨胀点源模拟长白山火山区岩浆囊压力变形源

2002～2003年长白山火山区的水平运动缺乏中心对称性分布特征,本文采用同一深度沿三个互相垂直方向拉张的点源模式模拟了长白山火山区的岩浆囊压力变形源. 在介质参数为V_p＝6.700 km/s,V_s＝3.8700 km/s,ρ＝2900 kg/m³的均匀弹性半空间中, 联合利用2002～2003年长白山火山区的GPS观测数据和水准数据, 结合PSGRN/PSCMP代码和遗传算法反演长白山火山区的压力变形源特征. 反演结果表明,由三个垂直的不等量扩张方向确定的椭球状点源模型可以较好地模拟长白山火山区岩浆囊压力变形源. 岩浆囊深处9.2 km,在走向为34.4°,倾角为82.2°的断层面的法向方向体积扩张量最大,达到7000000 m³;在走向为3028°,倾角为78.5°的断层面的法向方向体积扩张量次之,达到6598071 m³;在走向为337.7°,倾角为14.0°的断层面的法向方向体积扩张量最小,达到5220160 m³.     

基于音频大地测深法研究乌拉嘎断裂

乌拉嘎断裂位于嘉荫凹陷的东缘,对乌拉嘎金矿的成矿的起着关键作用.地质推断的乌拉嘎断裂位于金矿的北部呈NNE展布,理论上可通过电阻率差异分辨出乌拉嘎断裂切实可行.本文采用音频大地电磁测深法(AMT)来观测乌拉嘎断裂的电性特征并获得了高分辨率的视电阻率剖面图,本次观测结果很好的揭示了乌拉嘎断裂的电阻率特征,我们试图通过研究乌拉嘎断裂进而对本区的深部成矿机制进行研究.本次工作初步厘定了乌拉嘎断裂的宽度、深度,实践证明在乌拉嘎地区采用音频大地电磁测深法探测大型断裂构造具有良好效果.     

2016年3月2日苏门答腊MS7.8地震同震位移和应力场数值模拟研究

苏门答腊地区发育多条左旋走滑性质断层,地震活动活跃.2006年3月2日该区西南海域发生了M_S7.8大地震.大地震的发生常常会引发区域位移场和应力场及周围断层应力状态发生变化.本文建立全球PREM有限元地球模型,据已有的断层滑动模型计算了此次苏门答腊地震引发的同震位移和应力及库仑应力变化,并进一步讨论了此次地震对周围断层的影响以及区域构造应力场对库仑应力变化计算的影响.初步结果表明此次苏门答腊M_S7.8地震造成较大的南北向水平位移且集中在探测者破裂区（Investigator Fracture Zone）,最大水平位移量约6.74 m,断层倾角接近垂直,下盘向北运动而上盘向南,进一步表明M_S7.8地震为典型的左旋走滑为主的地震,发生海啸的可能性较低;库仑应力变化达MPa量级的区域集中在震中,但近场大部分余震分布在库仑应力减小区域,有效摩擦系数变化和区域构造应力场的耦合作用可能是其原因;利用改进的库仑应力变化计算方法和最优破裂方向计算得出的结果显示库仑应力触发理论可较好地解释余震分布.     

外核黏性对地磁场发电机数值模型的影响

地磁场发电机过程由数学方程组和输入参数所控制.为研究发电机参数对系统的影响,本文使用MoSST模型模拟了在外核黏性ν变化下发电机模型的输出.通过使用跨越近3个量级的ν值,着重研究了各物理场及其典型尺度随黏性的变化.发现黏性变化显著影响流场;磁场随ν增加而近乎单调减小,但变化幅度不超过30%;温度扰动随ν增加而小幅（6%）单调增加.经过拟合,得到外核流速u和黏性ν的比例关系：u~ν^0.49.流速随黏性增加的现象本质上是由于黏性增加打破了Taylor-Proudman约束,使得临界Rayleigh数减小,从而在相同的驱动力下带来了流速的增加.此外,作用力平衡分析发现,随着ν的变化,系统在几种平衡模式间切换.通过与之前比例关系的研究对比,本研究支持在一定范围内,磁场与黏性关系不大的结论;但反对流速与外核黏性无关的假设.     

利用GPS观测资料划分现今地壳活动块体的方法

利用GPS观测资料研究地球动力学问题是目前地学中的热点。利用GPS观测的速率(或位移)计算网格点附近小区域的应变张量和转动的欧拉极．再运用聚类分析法对计算的欧拉极进行分类，由此划分出地壳活动的次级块体。最后以中国大陆GPS观测结果为例，将中国大陆划分为9个活动块体，通过比较发现，所得结果与目前其他方法给出的结果大体一致。