厄瓜多尔俯冲区震间闭锁与粘弹性变形研究北大核心CSCD

以厄瓜多尔俯冲区为例,建立区域三维纯弹性和粘弹性模型,并基于有限元数值模拟方法分析其震间变形过程,同时以十余年的震间GPS速度场为约束,估算不同模型下的最佳震间闭锁深度。结果表明,在厄瓜多尔俯冲区,基于弹性模型的最佳震间闭锁深度约为30~35 km,而基于粘弹性模型的最佳闭锁深度为25 km。同时,相比于弹性模型,粘弹性模型能够更好地解释震间地表变形,说明该区域深部粘弹性圈层在俯冲板块的推挤作用下会引起较为明显的上地壳变形。     

基于正交贴体网格的VTI介质地震波数值模拟

为了解决规则网格有限差分法在处理起伏地表自由边界条件时需进行复杂坐标转换和插值运算的问题,以及对起伏地表进行阶梯状离散近似所产生的虚假散射波问题,将计算流体力学中的正交贴体网格生成技术引入到起伏地表下VTI介质的网格剖分中,采用基于同位网格的DRP/opt MacCormack有限差分法模拟起伏地表下VTI介质中的地震波场。网格的正交性使实施自由边界条件时无需做复杂的坐标转换和插值运算。数值算例表明,本文提出的网格剖分法能够有效消除阶梯状网格产生的虚假散射波,从而提高起伏地表下VTI介质地震波数值模拟精度。     

基于降雨强度-历时评价边坡稳定性——以安徽省3个边坡为例

为了探究降雨条件对边坡稳定性的影响,确定降雨阈值曲线,利用FLAC^2D 7.0数值模拟软件对安徽省3个小型边坡进行饱和-非饱和渗流及稳定性分析,获取了边坡安全系数随降雨强度和历时的变化趋势。边坡监测位置处孔隙水压力与饱和度的变化结果显示：降雨对边坡的影响是有限的,降雨历时越长,极限降雨强度越小并趋于稳定,边坡趋于饱和,此时安全系数达到最小值,后续降雨对边坡安全系数影响很小;为确保安全,假设临界安全系数为1.15,对应降雨强度为临界降雨强度,降雨历时越长,临界降雨强度越小并趋于稳定,3个边坡临界降雨强度趋于稳定时分别为13、12和16 mm·d^-1,对应的降雨历时分别为5、9和12 d,表明前期降雨对边坡稳定性影响很大。利用这3个边坡的临界降雨强度数据进行降雨阈值曲线拟合,并以435组安徽省2008～2022年历史灾害降雨数据验证该降雨阈值曲线。结果表明：95%以上的灾害数据点可以被该降雨阈值曲线预测。     

变井径环空分形表征及其对固井顶替效率的影响

为分析井径非线性变化对油气井注水泥顶替效率的影响,利用W-M分形函数、体积填充等理论,研究变井径环空几何形态表征及注水泥顶替效率计算,建立实测井径曲线的分形函数特征参数计算方法及变井径环空顶替流动分析模型,开展注水泥顶替流动数值模拟。结果表明：变井径环空的顶替效率低于均匀井径环空的顶替效率;标度参数、分形维数增大,井壁壁面凸起、凹陷高度增大,不规则程度增高,顶替效率降低;高标度参数、高分形维数范围内,顶替效率降低幅度高,顶替效率相对变化率大;顶替注入流速增大,顶替效率相对变化率减小;井斜角增大,顶替效率相对变化率增大;计算顶替效率时,当标度参数、分形维数低于阈值时,可将井眼环空近似考虑为均匀井径环空,而当标度参数、分形维数较高时,采用均匀井径环空简化处理会造成很大偏差;工程中,应注意处理计算模型选择及有效解决顶替效率问题。该结果为解决注水泥顶替效率计算的准确性、推进科学化固井设计与施工提供基础。     

星系形成的超大型计算机模拟

扼要地就数值模拟在星系形成研究中的作用，数值模拟与解析模型之间的关系，数值模拟的局限性等几方面的问题作了综述，然后介绍了马普小组与日本东京大学最近的合作情况，包括合作得到的模拟样本，研究成晨和正在开展的工作。     

环渤海区域及邻区现今地壳构造运动形变特征分析

将环渤海区域划分为5个次级块体,利用整体旋转与线性应变模型分析2009~2014年中国大陆构造环境监测网络的GPS对地观测速度场数据,得到相对于欧亚板块0.5°×0.5°的水平运动速度场和应变场。分析形变速度场和应变场的空间变化发现,环渤海区域整体上呈现0.25×10^-9/a的NW-SE 111.3°的双向趋势性扩张运动,太行块体、山西块体分别呈NW-SE 116.3°、NW-SE 130°的微弱扩张运动,胶辽块体、冀鲁块体和阴山-燕山块体分别存在NW-SE 144.3°、NE-SW 39.5°和NW-SE 155.6°的微弱压缩运动,其中太行块体的扩张运动量级相对于其他次级块体较大。尽管环渤海区域内部各次级块体的形变在空间上存在以上差异,但整体上环渤海区域主压应变轴方向基本为NEE-SWW,主张应变轴方向为NNW-SSE,且主压应变轴走向为NE 47.70°～89.74°,与利用地球物理方法得到的主压应变轴优势方向大体一致,表明环渤海区域现今地壳运动相对稳定。     

多级次地下水流系统的最小能耗率原理初探

20世纪60年代初期,Tóth基于定水头上边界条件推导出解析解,得出多级次地下水流系统,是水文地质学里程碑式的突破,成功地解决了一系列理论和实际问题。但Tóth解析解存在的缺陷也长期沿袭：单纯重视数学模拟而忽视物理机制;将地形控制地下水位看成是普适性规律;忽视给定水头上边界数学模拟的失真。这些缺陷,尤其是忽视物理机制探究,不仅妨碍Tóth理论自身发展,而且导致地下水流系统理论尚未被国际水文地质界普遍接受。参照河流动力学中应用的最小能耗率原理,类比提出地下水流最小能耗率的表达式。基于已有的通量上边界地下水流模式数值模拟结果,进一步探究物理机制,归纳得出地下水流系统遵循最小能耗率原理的结论。     

弹性波正演模拟中PML吸收边界条件的改进

在弹性波有限差分正演模拟中, 完全匹配层(PML) 吸收边界条件是使用广泛、吸收效果最好的吸收边界条件.在目前的两种PML实现方法中, 分裂形式的完全匹配层(SPML) 方法计算存储量大、编程实现复杂; 非分裂形式的完全匹配层(NPML) 方法计算效率低、计算过程复杂.针对传统PML吸收边界条件在实现过程中存在的问题, 推导出了一种简洁有效的非卷积实现的NPML吸收边界条件, 既不需要对场分量进行分裂, 也不需要做复杂的卷积运算.分析结果表明, 本文实现的NPML吸收边界条件不仅具有良好的吸收衰减性能, 而且计算方程简单, 编程实现容易, 占有内存更小.      

利用量子力学势能数据预测温度、压力、盐度和沉积物孔径对甲烷水合物形成和分解的影响

介绍一个预测不同温度、压力、盐度和沉积物毛细管孔径条件下甲烷水合物 溶液 气体多相平衡模型。该模型以Van der Waals和 Platteeuw热力学模型、量子力学从头算粒子相互作用势能、DMW 92状态方程和Pitzer电解质理论为基础,能在很宽广温压范围内预测温度、压力、盐度和毛细管力对甲烷水合物形成和分解的影响。通过对比本模型的预测结果与实验数据,可知本模型能够准确地预测海水和多孔介质中甲烷水合物的相平衡条件。对于一定盐度下多孔介质中甲烷水合物的形成温压条件的在线计算可浏览: www.geochem model.org/models.htm。     

辽河西部凹陷南段油气运移史研究与有利聚集区预测

通过流体包裹体均一温度测试、含油气包裹体指数GOI统计和盆地模拟,研究了渤海湾盆地辽河西部凹陷南段的油气运移史,预测了油气聚集的有利地区.包裹体均一温度资料表明,本区有三期成藏:东二期、东一期和明化镇期;以前两期为主.盆地模拟研究表明,尽管沙河街组烃源岩在沙一期末开始生烃,但在东营期才开始大规模排烃和成藏.GOI的统计分析表明,GOI>5%的区域为油气聚集的有利地区,其中的圈闭应成为下步首选的钻探目标,已发现的欢喜岭油田和双南油田就分布在GOI>5%的区域内.