收藏本站
《中国地震局地质研究所》 2010年
收藏 | 手机打开
二维码
手机客户端打开本文

并行计算三维数值模拟在华北地区现今构造变形分析中的应用研究

胡勐乾  
【摘要】:华北地区历史悠久、人口稠密、经济文化发达,是中国的心脏地带。华北地区地质构造演化过程复杂、构造类型多样、构造变形独特、地震活动性强,历来是众多学者研究的焦点。数值模拟方法在地质学方面的应用时间并不长,但已成为研究地学问题不可或缺的手段。在华北地区现有的地质构造、地球物理研究资料基础上,利用高性能并行计算技术,建立逼近实际的华北地区三维有限元模型,以高精度的GPS观测研究结果为运动边界,通过数值模拟方法研究华北地区的构造变形特征、断层运动性质、地震活动规律以及大地震对华北地区的影响等具有十分重要的意义。 数值模拟方法对华北地区构造变形分析的数值模拟虽然取得了一系列非常重要的成果,但依然存在不少问题:1)由于计算机软硬件条件的限制,以往的模型网格较为粗糙,计算精度较低。2)断层在大陆岩石圈的新构造运动中起着重要的作用,华北地区的活动断裂分布密集,活动性强,而现有的数值模拟研究多数仅考虑华北地区的几条主要断裂。3)大地震对地表运动、形变及区域构造应力场影响很大,这些影响不仅仅局限于震中附近区域,有时会传递至很远区域。例如,2008年汶川8.0级地震对全国,特别是华北地区的构造变形和未来的地震活动有什么影响,是值得研究的一个问题。4)数值模拟实验中,模拟结果的对比是工作的重要环节之一,其用来检验数值模型的合理性及数据的可靠性,这在很大程度上决定了数值模拟结果的可信度。遗憾的是,这个环节常常被众多研究人员所忽视,或者仅进行某一方面的结果对比就认定数值模型的可靠性。 针对上述问题,本研究主要开展了如下几个方面的工作: (1)华北地区三维动力学模型建立 基于并行版ANSYS数值模拟软件平台,根据华北地区活动地块划分及活动断裂分布,结合GPS资料等,确定模型的几何边界范围是99.8°~121.4°E,27.9°~42.3°N,模型的边界范围包括了华北活动地块区的绝大部分以及周边青藏、西域、南华和东北亚地块的部分地区。 华北活动地块区由鄂尔多斯、华北平原和鲁东-黄海三个二级活动地块所构成。华北地区内部主要活动断裂带有鄂尔多斯西南弧形断裂束、银川-吉兰泰断陷带、河套断陷带、张家口-渤海断裂带、渭河断陷带、山西断陷带、郯庐断裂带、唐山-磁县断裂带及安阳-徐州断裂带。模型考虑了华北地区已知所有晚更新世以来的活动断裂,断裂宽度取5 km。 模型单元平均边长为25km,模型共划分为单元416 582个,节点582 392个。模型利用1999~2004、2004~2007两期GPS数据作为边界约束,计算得到华北地区构造变形的数值模拟结果。 (2)地壳运动速度模拟与分析对比 在GPS边界约束条件下,模拟了GPS站点运动速度,并与实际观测的GPS运动速度进行对比,对物理模型和数值模拟计算结果的可靠性进行评价。 模拟结果表明,华北地区地壳运动速率的整体分布呈现由东向西逐渐减小、自北向南逐渐增大的特点。华北地区北部、燕山地块的地壳运动速率总体较小,而在华北地区中部、南部以及华南地块,运动速率逐渐增大。此外,2004~2007年相比1999~2004年,华北地区各地块地壳运动速度加快了1mm/a左右。 分析对比结果显示,数值模拟结果与GPS观测运动速度大部分比较吻合,只有少数点在运动方向或速度大小上有较大差异。其中,差异较大的点主要分布在华北平原块体的东北部,其原因可能与2006年7月4日文安5.1级地震的孕育发生有关,也可能是数值模拟的边界效应和GPS观测误差所致。 (3)断层运动性质模拟与分析对比 华北地区有限元模型考虑了华北地区已知所有晚更新世以来的活动断裂,根据数值模拟的计算结果,将模拟结果所得的断层运动性质与活动断裂地质调查的运动性质进行比较,结果显示,在统计的48条具有走滑性质的活动断裂中,1999~2004模拟结果共有43条断裂活动性质符合,符合率为89.6%;2004~2007模拟结果共有45条断裂活动性质符合,符合率为93.8%。数值模拟结果能较好的反映断层的运动性质,说明华北地区数值模型具有一定的合理性。另外,利用1999~2004数据进行数值模拟得出的断层运动性质与杨国华等根据高精度GPS复测资料反演水平位移场得出的断层运动性质进行比较,结果显示数值模拟得到的断层运动性质与实测GPS资料得出的结果也吻合得较好。 (4)应力应变场模拟与分析对比 通过数值模拟,计算了1999~2004和2004~2007华北地区地表第一主应力、第三主应力(最大压应力)、第一主应变和第三主应变(最大压应变)的大小和方向。结果显示,华北地区的水平应变场总体上比较一致:第一主应变的方向为北北西-南南东,第三主应变的方向为北东东-南西西。这与前人利用震源机制结果、GPS测量结果等资料反演得到的华北地区主应变场比较一致。将数值模拟得出的华北及邻区主应变场与张静华等得出的华北及邻区主应变图进行比较,两者也具有比较好的一致性。 (5)汶川大地震对华北地区构造变形的影响分析 在模型的合理性得到验证的基础上,初步模拟了汶川大地震对华北地区应力应变场及断层活动方式的影响。通过比较汶川地震引起华北地区附加应力应变场和华北地区的基本应力应变场,发现在华北平原地块和鲁东-黄海地块,附加应力应变场与基本应力应变场的主应变方向差异较大(45°),说明汶川地震对这些地区应力应变场的影响起卸载作用;而在鄂尔多斯地块,两者的主应变方向差异较小,说明汶川地震对其区域应力应变场的影响起加载作用。模拟汶川地震引起的断层运动性质与断层原有运动性质相比,主要表现为同向运动的断裂带有银川-吉兰泰断陷带和山西断陷带,表现为反向运动的断裂带有河套断裂带、安阳-徐州断裂带和郯庐断裂带。 通过华北地区现今构造变形的三维数值模拟研究,得到如下主要结论: (1)华北地区地壳运动速率的整体分布呈现由西向东、自北向南逐渐增大的特点。华北地区数值模拟节点运动速度与观测GPS测点运动速度大部分比较吻合。少数差异较大的点主要分布在华北平原块体的东北部,其原因可能与2006年7月4日文安5.1级地震的孕育发生有关,也可能是数值模拟的边界效应和GPS观测误差所致。 (2)华北地区大部分地块模拟平均运动速度与GPS平均运动速度近似。2004~2007年相比1999~2004年,各地块无论是GPS平均运动速度还是平均模拟值都提高了1mm/a左右,反映了2004~2007年华北地区地壳运动速度加快。 (3)数值模拟结果能较好地反映断层的运动性质。在统计的48条具有走滑性质的活动断裂中,1999~2004模拟结果共有43条断裂活动性质与地质调查结果相符,符合率为89.6%;2004~2007模拟结果共有45条断裂活动性质符合,符合率达93.8%。 (4)华北地区的数值模拟水平应变场与前人利用震源机制结果、GPS测量结果等资料反演得到的华北地区主应变场一致:第一主应变的方向为北北西-南南东,第三主应变(最大压应变)的方向为北东东-南西西。模拟结果还显示,活动断裂对模拟结果的影响较大,活动断裂成为应力应变的较为集中的区域。 (5)模拟结果显示,汶川地震对华北平原地块和鲁东-黄海地块的区域应变场的影响起卸载作用;而对鄂尔多斯地块的区域应变场的影响起加载作用;汶川地震引起断层同向运动的断裂带有银川-吉兰泰断陷带和山西断陷带,引起断层反向运动的断裂带有河套断裂带、安阳-徐州断裂带和郯庐断裂带。
【学位授予单位】:中国地震局地质研究所
【学位级别】:硕士
【学位授予年份】:2010
【分类号】:P542

【相似文献】
中国期刊全文数据库 前10条
1 王春华,廖素琼,钱瑞华,谢卫江;应用全息干涉法定量分析华北地区现代构造应力场的分布特征[J];地震地质;1983年04期
2 杨卫平;;华北平原地区猛禽观察[J];大自然;1992年01期
3 王琳瑛,舒曦;1994年11、12月华北地区M_L≥3.0地震动态[J];国际地震动态;1995年05期
4 马京津;高晓清;曲迎乐;;华北地区春季和夏季降水特征及与气候相关的分析[J];气候与环境研究;2006年03期
5 胡海濤;戴广秀;;略談华北地区修建水库的几个工程地质问题[J];中国地质;1961年11期
6 杜恆勤;;麻蜥的生态观察[J];生物学通报;1964年01期
7 李绍柄;刘德林;许建恩;张振江;;华北几个地震带内强震迁移活动与断块掀斜运动[J];华北地震科学;1984年04期
8 金清波;王文林;李卫国;;刺猬的染色体组型分析[J];河南师范大学学报(自然科学版);1985年04期
9 孟宪东;龚平;金赤兵;;用烈度衰减的椭圆模型标度华北地区历史地震震级[J];东北地震研究;1987年03期
10 李月洪,刘煜;向外长波辐射与我国华北地区旱涝关系的初步研究[J];应用气象学报;1991年02期
中国重要会议论文全文数据库 前10条
1 王俊国;;华北地区目前的地震活动趋势[A];新方法和新技术在地震预报中的应用暨中国东部地区震情研讨会论文摘要集[C];2004年
2 潘滇;徐建芬;彭卫刚;;华北夏季降水与东亚夏季风[A];中国气象学会2007年年会气候学分会场论文集[C];2007年
3 姜明明;周仕勇;;多维应力释放统计模型(MSRM)在华北地区地震危险性区划中的应用[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
4 黄建华;姜盛松;管波;傅容珊;;华北地区地壳剖面粘弹性应力分析[A];1995年中国地球物理学会第十一届学术年会论文集[C];1995年
5 ;华北地区古生物图册[A];中国地质科学院文集(1981)[C];1983年
6 蒋溥;梁小华;;华北地区地震断错及其形变场和应力场[A];中国地震学会第三次全国地震科学学术讨论会论文摘要汇编[C];1986年
7 李振清;;回顾华北地区机械设计教学研究会[A];第五届全国机械设计及制造专业教学研讨会议论文集(卷2 教学史志)[C];1995年
8 赵秀兰;延晓冬;;近20年华北地区土壤水储量变化趋势及水分管理与调配对策[A];2008年北京气象学会科技优秀论文集[C];2008年
9 王华林;王健;梁更国;杨丽军;;中国华北地区大地震复发活动周期、迁移特征及趋势预测[A];1995年中国地球物理学会第十一届学术年会论文集[C];1995年
10 刘长海;;华北地区最近二次中强地震前Hurst指数H值的异常变化[A];1996年中国地球物理学会第十二届学术年会论文集[C];1996年
中国重要报纸全文数据库 前10条
1 记者 张建华 通讯员 陈宁;加强政风行风建设 切实维护群众利益[N];长春日报;2007年
2 记者 李贺林;全新定位 全新形象 规模空前[N];中国建材报;2006年
3 钟陆史;东北华北地区纪检监察工作座谈会召开[N];中国纪检监察报;2007年
4 张立岩;联动演练检验华北供气应急机制[N];中国石油报;2007年
5 本报记者 李卫中;“金虎”收购“早早”欲谋亮剑华北(下)[N];山西经济日报;2007年
6 孙德富;增强大局意识 确保思想不散[N];中国民航报;2008年
7 刘国安;华北地区15个科研单位“设坛论见”思想政治工作[N];中国石油报;2003年
8 记者 明慧;民航五项指标考评奥运预演[N];中国改革报;2006年
9 张湘;华北片区矿产资源潜力评价会议召开[N];中国矿业报;2007年
10 刘弘;东方海外投资保税物流[N];国际商报;2006年
中国博士学位论文全文数据库 前10条
1 时应敏;松辽盆地长岭断陷火山机构及天然气成藏特征研究[D];中国地质大学(北京);2012年
2 徐宏节;松辽盆地南部断陷层火山岩储层识别及评价研究[D];成都理工大学;2010年
3 刘峡;华北地区现今地壳运动及形变动力学数值模拟[D];中国科学技术大学;2007年
4 宫飞;华北地区结构型节水种植业模式及途径研究——以北京市顺义区为例[D];中国农业大学;2003年
5 王桂芝;中国华北地区东方蜜蜂多样性及其分类研究[D];山东农业大学;2007年
6 荣艳淑;大范围气候变化与华北干旱研究[D];南京信息工程大学;2004年
7 张利平;华北地区陆气水循环时空演化规律研究[D];武汉大学;2005年
8 温海清;金元之际的华北地方行政建置[D];复旦大学;2008年
9 张学民;华北地区壳幔剪切波速度结构及与晚新生代地质演化关系研究[D];中国地质大学(北京);2005年
10 郭锦超;近代天津和华北地区经济互动的系统研究(1880年代—1930年代)[D];南开大学;2004年
中国硕士学位论文全文数据库 前10条
1 胡勐乾;并行计算三维数值模拟在华北地区现今构造变形分析中的应用研究[D];中国地震局地质研究所;2010年
2 罗遵兰;华北地区冬小麦与夏玉米水分生产函数模型研究[D];中国农业大学;2005年
3 孙燕;华北地区降水异常及其与北太平洋海温、100hPa高度场异常的关系[D];南京气象学院;2002年
4 史晓娟;全面建设小康社会目标下的华北地区电力需求预测[D];华北电力大学(北京);2006年
5 郭士忠;清代华北地区粮食流通[D];山西大学;2008年
6 张艳;清代典当业税收问题[D];陕西师范大学;2007年
7 左艳艳;华北地区中老年针织毛衫感性消费研究[D];江南大学;2009年
8 范烨;华北地区气溶胶分布和层积云结构特征研究[D];南京信息工程大学;2006年
9 韩军彩;华北地区空中水汽含量的演变特征[D];南京信息工程大学;2011年
10 李坤;中国华北地区人群TAP1基因多态性与乙肝的关联性研究[D];山东大学;2005年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026