无网格局部Petrov-Galerkin法大地电磁二维正演

将无网格局部Petrov-Galerkin算法用于大地电磁二维正演。介绍了该方法的基本原理;从大地电磁二维边值问题出发,利用子域法详细推导了与之对应的局部Petrov-Galerkin弱式方程,并用高斯积分法将其离散化。论述了无网格局部Petrov-Galerkin法较无单元Galerkin法及有限元法的优缺点,最后通过二维模型的计算验证了算法的有效性。     

断层阶区对震源破裂传播过程的控制作用研究

地震破裂能否穿越断层阶区（stepover）引发更大震级的地震是震源动力学研究的重要内容.本文利用不连续变形体接触力学的动态有限单元方法,模拟断层阶区对地震破裂传播的控制作用.通过改变断层周边初始应力场、断层面上的摩擦本构关系以及断层阶区的间距大小来分析各个因素对破裂传播的影响,并定量分析产生这些影响的力学机制.模拟结果表明：断层面上的摩擦系数减小或断层周边区域内的初始剪应力增大,都将增加断层破裂跳跃阶区传播的可能性;此外,若断层阶区间距越小,断层破裂也越容易跳跃阶区传播.计算结果还显示：断层上的摩擦系数大、初始剪应力小、断层阶区间隔大,那么此阶区所在之处将可能是断层破裂的终止位置;相反,当断层面上的摩擦系数较小、初始剪应力较大、断层阶区间隔较小,破裂就容易穿越阶区而出现较大的地震.同时,从模拟结果可以看出,在发震断层破裂停止后,应力将继续向四周传播;当应力积累达到破裂极限时,触发断层阶区中的另一断层产生破裂,因此在破裂跳跃断层阶区的过程中存在一个时间延迟.最后,破裂能否跳跃断层阶区,可以利用库仑应力在空间的分布进行合理的解释.     

断层自发破裂动力过程的有限单元法模拟

断层自发破裂动力过程的研究对于认识地震过程及减轻地震灾害有着重要的科学意义.为合理地模拟断层的自发破裂过程,本文首先对经典的滑移弱化摩擦关系进行了改进,然后利用有限单元方法对破裂过程进行动态数值模拟.模拟结果表明,利用改进后的摩擦关系能够产生脉冲型（pulse-like）破裂模式,而经典的滑移弱化摩擦关系不能产生这种破裂形态.模拟结果还显示,断层自发破裂过程受初始应力场及摩擦关系影响,当初始应力场中剪应力水平较低时,容易产生脉冲型破裂;但当初始剪应力较高时,会产生裂纹型（crack-like）破裂.这个现象与在实验室里进行的岩石破裂实验结果是一致的.在相同的初始应力情况下,若滑移弱化摩擦本构关系中的动摩擦系数较大,断层将易于产生脉冲型破裂;若动摩擦系数较小,将倾向于产生裂纹型破裂.此外,本文也采用速率弱化摩擦关系对断层自发破裂过程进行了模拟,结果发现,在初始场及其他条件不变时,如果摩擦关系中的b-a值较小,容易产生脉冲型破裂;如果b-a值较大,会产生裂纹型破裂.     

华东片区Ⅳ级成矿单元划分及成矿地质特征

基于中国大陆东南部的构造—成矿背景及区域矿产分布地质特征,参考华东片区全国Ⅲ级成矿单元划分方案,进一步划分出28个Ⅳ级成矿单元,并论述了其中8个Ⅲ级成矿单元及19个Ⅳ级成矿单元的成矿地质特征.     

台湾碰撞带现今地壳变形场特征及其动力学成因的有限单元法模拟研究

台湾是菲律宾海板块与欧亚板块会聚、碰撞的产物,区内地质构造复杂,地震频发,变形强烈.为定量研究台湾地区变形特征及其动力学成因,文中运用二维不连续变形体弹性力学的有限单元计算方法,利用台湾地区1995—2005年GPS的观测结果作为边界约束,对台湾地区地壳变形场进行了模拟.计算结果显示,台湾中部地区主要呈压缩状态,但在台湾的东北部及南部地区出现了拉张的变形环境.变形程度最大的区域位于台湾东部海岸山脉及其附近海域,同时,计算给出台湾纵谷断裂滑移速率约为13.8~23.5 mm/yr,由于纵谷断裂对变形的吸收,因此变形在纵谷断裂以西及西北地区迅速衰减.此外,计算结果还发现,计算给出的台湾岛上的速度值与GPS观测结果吻合得较好;计算给出的主应力方位与地应力观测及震源机制解结果也颇为一致.由此说明文中的有限元模型是合理的.此外,计算结果还表明,菲律宾海板块向北西方向的推挤,板块边界形状,观音、北康高地的抵阻, 冲绳海槽扩张及琉球海沟的向南后撤以及断裂作用等,共同造成了台湾地区现今变形场的主要格局.     

电导率各向异性的海洋电磁三维有限单元法正演

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性, 文中选用二次场算法, 将背景电阻率设置为水平层状且各向异性, 场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组, 通过Galerkin加权余量法形成大型稀疏有限元方程, 采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位; 进而, 利用滑动平均方法得到二次矢量位和标量位在空间的导数, 由此得到二次电磁场; 通过一维模型对算法的可靠性进行验证, 与此同时, 针对实际复杂海洋电磁模型, 比较有限元模拟结果与积分方程模拟结果, 进一步验证算法精度.若干计算结果均表明, 文中算法具有良好的通用性, 适用于井中电磁、航空电磁, 环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.     

基于半监督学习的井震联合储层横向孔隙度预测方法

传统地震储层预测技术一般基于弹性参数反演和岩石物理建模的级联流程实现储层孔隙度预测, 其预测精度受到波动理论和岩石物理理论的近似假设、初始模型和二次反演累积误差等因素的影响.为缓解这些问题, 本文提出了一种基于双向门控递归单元神经网络的半监督学习井震联合孔隙度预测方法, 实现从地震数据直接预测储层横向孔隙度.通过少量的地震测井样本标签对和多目标函数约束建立智能化多尺度多信息融合孔隙度预测模型, 实现地震数据到孔隙度, 孔隙度再到生成地震数据的闭环映射.此外, 在网络模型每次迭代更新的过程中随机引入非井旁地震道参与网络训练, 非井旁地震道的波形匹配能在一定程度上保证井间孔隙度的预测精度.模型数据和实际数据测试结果表明, 本文提出的方法相比于有监督学习孔隙度预测方法能进一步提高储层孔隙度的预测准确性和横向连续性, 获得较为可靠的储层物性参数的空间分布.

     

基于OpenSEES的高强钢组合K形偏心支撑框架精细化建模分析

为研究高强钢组合K形偏心支撑框架在反复荷载作用下的受力性能,在已有拟静力试验的基础上,利用有限元软件OpenSEES对一个高强钢组合K形偏心支撑框架模型进行了精细化建模研究。首先,探讨了模型的单元选择,纤维截面划分和材料本构参数定义问题。同时,考虑到在地震荷载作用下,高强钢组合K形偏心支撑框架主要依靠消能梁段的剪切变形来耗能,因此,如何模拟结构中消能梁段的剪切效应成为研究的关键。本文采用OpenSEES提供的组合截面和零长度单元两种方法来对消能梁段进行建模,剪切材料分别考虑了Steel 02和滞回材料。将不同消能梁段建模方法的模拟滞回曲线和骨架曲线与试验结果进行比较。结果表明:精细化模型的数值模拟结果和试验结果具有较好的吻合度,验证了数值模型的正确性和可行性;采用零长度单元和Steel 02的组合方式来模拟消能梁段的剪切效应能够同时满足计算效率和计算精度的要求,为进一步利用OpenSEES进行高强钢组合K形偏心支撑框架的动态响应分析和复杂模型分析提供了建模参考依据。     

流动单元研究在乌石A油田扇三角洲沉积体系复杂储层的应用

乌石A油田为近年来在南海西部盆地新发现的油田之一,主要发育近源扇三角洲沉积体系,通过评价井产能测试和测压取资料表现出研究区储层渗流特征复杂,产能与测井解释储层物性相关性差,常规方法无法准确表征储层渗流特征。针对海上油田资料较匮乏的特点,应用流动单元研究思路,以微观孔隙结构特征分类为基础划分流动单元,总结归纳各类流动单元的微观与宏观特点,提出以分选性、电性、孔隙度3参数在全区范围内划分流动单元以及评价其产能,评价结果与实测结果吻合度较高,目前来看此方法能准确评价此类复杂储层,有效指导该油田后期开发部署。     

医用CT在三维椭球颗粒随机堆积结构研究中的应用

颗粒随机堆积问题是物理和工程领域的热门问题。椭球颗粒随机堆积体系作为球形体系的拓展,增加了旋转自由度,从而表现出更复杂的性质,由于椭球颗粒比球形颗粒更类似于自然界中普遍存在的非规则颗粒,因此其研究具有重要的实用价值。研究三维颗粒堆积结构问题一般需要CT成像技术,本文介绍了如何利用医用CT对三维椭球随机堆积体系进行成像研究的技术细节,包括如何通过图像处理的方法获取椭球颗粒的空间位置、取向和体积,如何利用接触数标度函数拟合的方法得到近邻颗粒间的接触信息等。本文的研究内容和方法为利用各种X射线CT成像方法研究三维颗粒体系提供了技术参考。