中国地震科学实验场断层系统的三维计算网格模型

doi:10.12196/j.issn.1000-3274.2023.03.002

地震 ›› 2023, Vol. 43 ›› Issue (3): 18-33.doi: 10.12196/j.issn.1000-3274.2023.03.002

中国地震科学实验场断层系统的三维计算网格模型

张熔鑫^1,2,3, 邢会林^1,2,3, 舒涛^1,2,3, 刘骏标^1,2,3, 郭志伟^1,2,3, 王建超^1,2,3, 谭玉阳^1,2,3

1.深海圈层与地球系统前沿科学中心, 海底科学与探测技术教育部重点实验室, 中国海洋大学海洋地球科学学院, 山东青岛 266100;
2.青岛海洋科学与技术试点国家实验室, 山东青岛 266237;
3.中国海洋大学海底科学与工程计算国际中心, 山东青岛 266100

收稿日期:2023-03-02 修回日期:2023-05-12 出版日期:2023-07-31 发布日期:2023-08-28
通讯作者: 邢会林, 教授。 E-mail: h.xing@ouc.edu.cn
作者简介:张熔鑫(1998-), 男, 山东淄博人, 在读硕士研究生, 主要从事断层系统有限元模型构建与计算研究。
基金资助:
地震数值预测联合实验室开放基金(2020NEF02)

3D Mesh for the Fault System in China Seismic Experimental Site

ZHANG Rong-xin^1,2,3, XING Hui-lin^1,2,3, SHU Tao^1,2,3, LIU Jun-biao^1,2,3, GUO Zhi-wei^1,2,3, WANG Jian-chao^1,2,3, TAN Yu-yang^1,2,3

1. Frontiers Science Center for Deep Ocean Multispheres and Earth System, Key Lab of Submarine Geosciences and Prospecting Techniques, MOE and College of Marine Geosciences, Ocean University of China, Qingdao 266100, China;
2. Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266237, China;
3. International Center for Submarine Geosciences and Geoengineering Computing (iGeoComp), Ocean University of China, Qingdao 266100, China

Received:2023-03-02 Revised:2023-05-12 Online:2023-07-31 Published:2023-08-28

摘要/Abstract

摘要： 川滇地区受印度板块与欧亚板块碰撞运动影响地震频发, 是开展地震研究的理想实验场。为实现基于几何模型的三维断层系统有限元模型快速构建, 本文以川滇地区断层系统的主要断层几何模型为基础进行了以下研究。 ① 改进了推进波前法, 使其能够基于几何模型生成三维曲面的三角形网格, 并自动识别需要局部加密的区域以及进行网格加密。 ② 提出了交叉断层的交叉线识别算法, 根据非参数曲面间的空间几何关系识别不同面(包括断层面与研究区域边界)之间的相交线。 ③ 完善了三维断层系统有限元计算网格模型快速构建方法: 对原始断层几何模型进行延长、连接等操作, 使其能在几何上代表断层真实形状; 将研究区域外边界与内部断层面整合, 并识别它们之间的交叉线; 以交叉线为约束, 使用改进的推进波前法对断层面进行网格重划分, 修复网格拓扑关系和提高网格质量; 最后, 以断层面新生成的三角形面网格为约束, 自动生成含断层的四面体计算网格模型。 ④ 将上述方法应用于中国地震科学实验场, 构建了川滇地区断层系统的三维计算网格模型, 为利用有限元数值模拟研究该区域地震动力学过程奠定了计算网格模型基础。

关键词: 川滇地区, 断层系统, 有限元网格生成, 推进波前法

Abstract: The Sichuan-Yunnan region is subjected to frequent earthquakes due to the collision of the Indian and Eurasian plates, rendering it a ideal experimental site for seismic research. This paper presents a rapid method for constructing a finite element mesh model with complex faults based on a geometric mesh model of the region, using the example of the geometry model of major faults in fault system of Sichuan-Yunnan region. The following research has been conducted: ① The Advancing Front Technique (AFT) has been improved to generate a three-dimensional triangle mesh based on the geometric model, and automatically identify areas that require local refinement and conduct mesh refinement. ② An algorithm for identifying the intersection lines of intersecting faults has been proposed, and it can recognize the intersection lines between different surfaces (including fault surfaces and study area's boundaries) based on the spatial geometric relationship between non-parametric surfaces. ③ The rapid method for constructing a finite element mesh model with complex faults based on a geometric mesh model of the region has been improved. First, the original fault model is manipulated using expansion and connection operations to enable it to represent the true shape of the fault geometrically. Subsequently, the boundary of the study area and the fault mesh are combine. The intersection lines between different fault planes or between fault planes and the study area's boundary are then identified. Then, using the improved advancing front technique, the fault surface is remesh with the intersection lines as constraints to repair the mesh topology and improve the mesh quality. Finally, a tetrahedral calculation mesh model containing faults is automatically generated using the newly generated triangular mesh of the fault surface as constraints. ④ The above method has been applied to the China seismic experimental site to construct a three-dimensional calculation mesh model containing the main fault system which lays the foundation for using finite element numerical simulations to study the seismic dynamics in the study area.

Key words: Sichuan-Yunnan region, Fault system, Mesh generation, Advancing front technique

中图分类号:

P315.7

张熔鑫, 邢会林, 舒涛, 刘骏标, 郭志伟, 王建超, 谭玉阳. 中国地震科学实验场断层系统的三维计算网格模型[J]. 地震, 2023, 43(3): 18-33.

ZHANG Rong-xin, XING Hui-lin, SHU Tao, LIU Jun-biao, GUO Zhi-wei, WANG Jian-chao, TAN Yu-yang. 3D Mesh for the Fault System in China Seismic Experimental Site[J]. EARTHQUAKE, 2023, 43(3): 18-33.

参考文献

[1] 徐锡伟, 闻学泽, 郑荣章, 等. 川滇地区活动块体最新构造变动样式及其动力来源[J]. 中国科学(D辑), 2003, 33(S1): 151-162.
XU Xi-wei, WEN Xue-ze, ZHENG Rong-zhang, et al. The latest tectonic change patterns and dynamic sources of active blocks in Sichuan-Yunnan region[J]. Science in China (Series D), 2003, 33(S1): 151-162 (in Chinese).
[2] 吴中海, 龙长兴, 范桃园, 等. 青藏高原东南缘弧形旋扭活动构造体系及其动力学特征与机制[J]. 地质通报, 2015, 34(1): 1-31.
WU Zhong-hai, LONG Chang-xing, FAN Tao-yuan, et al. The arc rotational-shear active tectonic system on the southeastern margin of Tibetan Plateau and its dynamic characteristics and mechanism[J]. Geological Bulletin of China, 2015, 34(1): 1-31 (in Chinese).
[3] 张培震, 邓起东, 张国民, 等. 中国大陆的强震活动与活动地块[J]. 中国科学(D辑), 2003, 33(S1): 12-20.
ZHANG Pei-zhen, DENG Qi-dong, ZHANG Guo-min, et al. Strong earthquake activity and active blocks in Chinese Mainland[J]. Science in China (Series D), 2003, 33(S1): 12-20 (in Chinese).
[4] Allen C R, Luo Z L, Qian H, et al. Field study of a highly active fault zone: The Xianshuihe fault of southwestern China[J]. Geological Society of America Bulletin, 2016, 103(9): 1178-1199.
[5] Royden L H, Burchfiel B C, King R W, et al. Surface deformation and lower crustal flow in eastern Tibet[J]. Science, 1997, 276(5313): 788-790.
[6] Papadimitriou E, Wen X Z, Karakostas V, et al. Earthquake triggering along the Xianshuihe fault zone of western Sichuan, China[J]. Pure and Applied Geophysics, 2004, 161(8): 1683-1707.
[7] Wen X Z, Ma S L, Xu X W, et al. Historical pattern and behavior of earthquake ruptures along the eastern boundary of the Sichuan-Yunnan faulted-block, southwestern China[J]. Physics of the Earth and Planetary Interiors, 2008, 168(1-2): 16-36.
[8] Luo G, Liu M. Stress evolution and fault interactions before and after the 2008 Great Wenchuan earthquake[J]. Tectonophysics, 2010, 491(1-4): 127-140.
[9] Xing H L, Makinouchi A. Three dimensional finite element modeling of thermomechanical frictional contact between finite deformation bodies using R-minimum strategy[J]. Computer Methods in Applied Mechanics and Engineering, 2002, 191(37-38): 4193-4214.
[10] Xing H L, Makinouchi A, Mora P. Finite element modeling of interacting fault systems[J]. Physics of the Earth and Planetary Interiors, 2007, 163(1-4): 106-121.
[11] Xing H L, Mora P, Makinouchi A. A unified friction description and its application to the simulation of frictional instability using the finite element method[J]. Philosophical Magazine, 2006, 86(21-22): 3453-3475.
[12] Lu R Q, Wang M M, Li Z G. The fault model in China Seismic Experimental Site[Z]. CSES Scientific Products, 2019.
[13] Lu R Q, Liu Y D, Xu X W, et al. Three-dimensional model of the lithospheric structure under the eastern Tibetan Plateau: Implications for the active tectonics and seismic hazards[J]. Tectonics, 2019, 38(4): 1292-1307.
[14] He D F, Lu R Q, Huang H Y, et al. Tectonic and geological setting of the earthquake hazards in the Changning shale gas development zone, Sichuan Basin, SW China[J]. Petroleum Exploration and Development, 2019, 46(5): 1051-1064.
[15] Tan X B, Liu Y D, Lee Y H, et al. Parallelism between the maximum exhumation belt and the Moho ramp along the eastern Tibetan Plateau margin: Coincidence or consequence?[J]. Earth and Planetary Science Letters, 2019, 507: 73-84.
[16] Mavriplis D J. An advancing front Delaunay triangulation algorithm designed for robustness[J]. Journal of Computational Physics, 1995, 117(1): 90-101.
[17] Cuillière J C. An adaptive method for the automatic triangulation of 3D parametric surfaces[J]. Computer-Aided Design, 1998, 30(2): 139-149.
[18] Lan T, Lo S H. Finite element mesh generation over analytical curved surfaces[J]. Computers and Structures, 1996, 59(2): 301-309.
[19] Vollmer J, Mencl R, Müeller H. Improved laplacian smoothing of noisy surface meshes[J]. Computer Graphics Forum, 1999, 18(3): 131-138.
[20] Liu Y, Xing H L. An effective 3D meshing approach for fractured rocks[J]. International Journal for Numerical Methods in Engineering, 2016, 107(5): 363-376.
[21] 尹迪, 董培育, 曹建玲, 等. 川滇地区地震危险性数值分析[J]. 地球物理学报, 2022, 65(5): 1612-1627.
YIN Di, DONG Pei-yu, CAO Jian-ling, et al. Numerical analysis of the seismic hazard in Sichuan-Yunnan region[J]. Chinese Journal of Geophysics, 2022, 65(5): 1612-1627 (in Chinese).

中国地震科学实验场断层系统的三维计算网格模型

3D Mesh for the Fault System in China Seismic Experimental Site

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	郑秋月, 王青华, 刘东, 黄江培, 陈石. 基于时变重力数据的川滇地区场源模型反演和解释[J]. 地震, 2021, 41(1): 205-218.
[2]	徐伟民, 陈石, 阮明明, 杨锦玲, 韩建成, 李红蕾, 卢红艳. 川滇地区陆地流动重力测网场源分辨能力评估[J]. 地震, 2021, 41(1): 191-204.
[3]	李红蕾, 韩建成, 陈石, 侍文, 张贝. 基于等效源的重力异常空间导数计算方法及其在川滇地区的应用[J]. 地震, 2021, 41(1): 129-140.
[4]	陈学忠, 李艳娥. 川滇地区地震活动与地球自转速率变化之间的关系[J]. 地震, 2019, 39(1): 126-135.
[5]	巩浩波,李光科,廖欣,陈敏. 川滇地区井水位潮汐响应特性分析[J]. 地震, 2017, 37(1): 20-30.
[6]	王俊, 邵志刚, 孙小龙, 王熠熙, 向阳, 王喜龙. 川滇地区强震前流体异常特征与预测指标体系初探[J]. 地震, 2016, 36(4): 109-119.
[7]	刘月, 吕晓健, 田勤俭. 基于“区域-时间-长度算法”分析川滇地区强震前地震活动性变化^*[J]. 地震, 2016, 36(2): 94-104.
[8]	魏文薪, 江在森, 刘晓霞, 邹镇宇. 川滇地区应变率场分布及其变化特征研究[J]. 地震, 2015, 35(4): 11-20.
[9]	董曼, 杨天青, 郑通彦, 陈通. 川滇地区走滑型地震烈度衰减特征研究[J]. 地震, 2015, 35(3): 147-154.
[10]	何康, 晏锐, 王春, 郑海刚, 李军辉, 方震. 基于Wigner-Ville分布的地磁前兆异常提取方法及其在川滇地区的应用[J]. 地震, 2014, 34(4): 88-99.
[11]	杨文, 刘杰. 2012年4月11日苏门答腊8.6级地震前后川滇地区波速变化研究[J]. 地震, 2013, 33(4): 248-256.
[12]	康春丽. 川滇地区长波辐射场变化与地震活动关系研究[J]. 地震, 2008, 28(3): 43-48.
[13]	马宏生, 张国民, 周龙泉, 刘杰, 邵志刚, 夏红. 川滇地区中小震重新定位与速度结构的联合反演研究[J]. 地震, 2008, 28(2): 29-38.
[14]	韩渭宾, 蒋国芳. 川滇地区较长时间尺度的地震活动盛衰交替性[J]. 地震, 2005, 25(1): 51-57.
[15]	付虹, 刘丽芳, 张晓东. 川滇地区成组强震活动判别指标研究[J]. 地震, 2005, 25(1): 8-14.