收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

华北地区现今地壳运动及形变动力学数值模拟

刘峡  
【摘要】: 华北地区位于我国北部,阴山、燕山以南,秦岭、大别山以北、鄂尔多斯高原以东延伸至沿海一带的广袤区域,地理坐标范围为东经112°~124°,北纬31°~42°,是我国政治、经济、文化中心。同时,该区域囿于太平洋板块、菲律宾板块和印度板块之间,新生代以来构造运动活跃,为我国大陆强地震频发区域之一。因此,研究该区域现今地壳运动及形变、构造应力场特征,继而探讨其动力学机理,具有重要理论和现实意义。 本文采用有限元方法在大型有限元软件平台ANSYS基础上,综合现有地质、地球物理资料构建华北地区岩石层(地壳和上地幔顶部)框架。以区域GPS观测、地形变观测以及震源机制解结果作为约束条件及检验标准,通过数值模拟的方式,研究该区域现今地壳运动和形变、构造应力场特征,并将表层运动与深部运动相联系,以区域动力学环境和驱动力为中心,探讨其动力学机制。 根据华北地区新生代构造发育特征和深部地球物理观测数据,本文选取了108°E~127°E,32°N~42°N矩形区域构建有限单元数值模拟的几何框架。与此同时,还将现今活动断层及断裂带分布、地形起伏以及地壳—上地幔速度结构等纳入模型中。本文构建和测试了三个基本模型:即三维模型3DCM_Ⅰ,3DCM_Ⅱ利二维模型2DDM_Ⅰ。模型3DCM_Ⅰ不考虑地壳—上地幔的实际分层,将模型沿纵向分为七个厚度均匀的水平层,分别代表上、中、下地壳和上地幔(四层)。模型3DCM_Ⅱ则依据该区域深反射地震研究结果进行水平分层。在模型3DCM_Ⅰ、3DCM_Ⅱ中,区域内主要活动断裂带处理为宽度5—6Km,深度不等(10~15km)的软弱带,模型含14736个节点,24829个三维实体单元。在模型2DDM_Ⅰ中断层及断裂带则处理为非连续接触边界,模型含1843个节点,3547个二维实体单元。 论文采用有限元弹性静力分析方法,以模型区周边GPS观测作为边界约束,计算华北岩石层形变的水平速度场和应力、应变增量场。并以模型区内部GPS观测作为标准,通过调整模型中断层的物理参数方式,寻求与GPS观测的最佳吻合的有限元解。计算显示,水平运动场的最佳数值解和GPS观测值的平均离散度为10171mm~2/year(三维模型)、1.176mm~2/year(二维模型)。与此同时,数值模拟预测的区内主要活动断层的错动方式、水平错动速率(三维模型结果为~0.1mm/year,二维模型结果为~0.3mm/year)与地质及跨断层形变观测结果基本一致。三维模型还显示,在10~15Km深度上,最大、最小主应力轴接近水平,主张应力为主压应力的2~8倍,主张应力方向为NNW,主压应力方向为NEE。这一结果与震源机制解、地应力观测结果比较吻合。三维模型得到的全区应变率较低约10~(-9)/year。并且,大致沿汾渭地堑断裂系、唐山—河间—磁县断裂带和郯庐断裂带存在三个北东走向的剪应变梯度带。不过二维模型预测的张家口—渤海断裂带水平运动速度为1.35~1.45mm/year,远高于其他断层以及三维模型结果。数值模拟结果可以推论:华北地区岩石层水平形变运动主要受控于周边大型构造块体的相对运动,其次才是岩石横向、纵向非均匀性的影响。同时,在周边块体运动的控制下,区内主要断裂带现今的错动方式基本是继承性的,且错动速率较低。只有张家口—渤海断裂带错动速度较高,值得关注。 论文考虑岩石应力—应变遵从幂指数本构关系,基于三维模型3DCM_Ⅰ,3DCM_Ⅱ的基本框架,采用弹性—蠕变静力学有限元方法对华北地区晚近时期(~4Ma)构造运动进行数值模拟,探讨华北地区现今地壳运动演化的动力学环境和驱动力因素。模拟时以模型周边GPS观测数据作为表层的边界约束,计算时顾及岩石的流变性、重力作用以及大变形导致的位移—应变非线性。本文通过选择侧面、底面边界约束方式,调整岩石的物性参以及分析板块运动、周边块体运动、断层活动、岩石层流变分层等因素的影响,寻求与区内GPS观测、构造应力场观测数据最佳符合的预测模型。论文计算了22个模型。结果表明,当侧面以深部稍快、浅部稍慢的方式运动,并考虑岩石层下部的拖曳运动时,则模拟预测果更接近实际观测。在该约束方式下,相应于模型3DCM_Ⅰ、3DCM_Ⅱ框架的预测与GPS观测的最低平均离散度分别为1.1529mm~2/year、1.1451mm~2/year,低于弹性模型,而且断层错动与地质研究结果相一致。在相同边界条件下,3DCM_Ⅱ结果优于3DCM_Ⅰ结果,考虑断层运动结果优于不考虑断层运动结果。模型预测应力场显示,深度10Km处的最小主压应力轴均接近水平,方向为NNW向,最大主压应力轴在南部区域接近水平,在北部区域则垂直于水平面,方向以NEE为主。这与其它研究显示的最大、最小应力主轴的方向基本吻合,但与现有的最大、最小主应力轴均接近水平的基本认识存在一定差别。预测应力场随深度的变化明显受控于边界加载方式(应力环境),而断层的影响仅局限于断层内部以及断层附近较小区域内。 为研究局部特别是隐伏构造对地表形变的影响,提取形变异常信息,论文发展了GPS差异形变分析方法。并用此方法,对华北地区和中国大陆GPS观测结果进行了GPS差异形变分析。华北地区GPS差异形变分析显示,张家口—渤海断裂带、唐山—河间—磁县断裂带为该地区现今构造活动较强烈区域,在未来若干年内需密切关注及防范强震发生。对中国大陆而言,论文分别使用三种有限元模型(均一、分块和多驱动力模型)而进行了GPS差异形变分析。结果表明,中国大陆现今构造运动是多种驱动力共同作用的结果。其中,印度板块向北的强烈推挤以及地幔对流对岩石层底部的拖曳作用占据着非常重要的地位。前者主导了以青藏高原为中心的大陆西部地区构造变形运动,而后者则对华北、华南地区的影响重大。同时,现今构造运动十分活跃大型断裂系(如阿尔金、张家口—渤海等)对大陆岩石层构造变形的影响也是不可忽视的。 综合数值模拟结果可以得到如下认识:华北地区岩石层形变的区域动力学环境和驱动力十分复杂。总体而言,一方面在太平洋板块的俯冲和印度—欧亚板块的碰撞挤压作用下,鄂尔多斯活动地块、华南活动地块以及东北亚活动地块的运动状态决定了华北地区表层运动的基本格局。另一方面,地幔对流对岩石层底部的拖曳作用将直接影响该区域岩石层形变运动。与此同时,以张家口—渤海断裂为代表的现今构造运动十分活跃大型断裂系以及岩石层内流变性非均匀分布(特别是中地壳软弱层的存在),对区域构造变形的影响也是不可忽视的。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 杨卫平;;华北平原地区猛禽观察[J];大自然;1992年01期
2 王琳瑛,舒曦;1994年11、12月华北地区M_L≥3.0地震动态[J];国际地震动态;1995年05期
3 马京津;高晓清;曲迎乐;;华北地区春季和夏季降水特征及与气候相关的分析[J];气候与环境研究;2006年03期
4 胡海濤;戴广秀;;略談华北地区修建水库的几个工程地质问题[J];中国地质;1961年11期
5 杜恆勤;;麻蜥的生态观察[J];生物学通报;1964年01期
6 李绍柄;刘德林;许建恩;张振江;;华北几个地震带内强震迁移活动与断块掀斜运动[J];华北地震科学;1984年04期
7 金清波;王文林;李卫国;;刺猬的染色体组型分析[J];河南师范大学学报(自然科学版);1985年04期
8 孟宪东;龚平;金赤兵;;用烈度衰减的椭圆模型标度华北地区历史地震震级[J];东北地震研究;1987年03期
9 李月洪,刘煜;向外长波辐射与我国华北地区旱涝关系的初步研究[J];应用气象学报;1991年02期
10 姚殿义;杜迎春;刘素英;;华北地区的震源机制[J];华北地震科学;1991年03期
11 张迺哲;;华北地区医学生化学术研讨会[J];生命的化学;1991年06期
12 黎彤,张西繁;华北花岗岩类的大地构造岩石化学特征[J];大地构造与成矿学;1992年04期
13 黄蓉良;1996年6~7月华北地区M_L≥3.0地震动态[J];国际地震动态;1996年11期
14 袁金荣,徐菊生,高士钧;利用GPS观测资料反演华北地区现今构造应力场[J];地球学报-中国地质科学院院报;1999年03期
15 范广洲,吕世华;陆面植被类型对华北地区夏季降水影响的数值模拟研究[J];高原气象;1999年04期
16 方剑;华北地区岩石圈三维密度分布[J];物探与化探;1999年03期
17 于静洁,任鸿遵;华北地区粮食生产与水供应情势分析[J];自然资源学报;2001年04期
18 范广洲,吕世华,程国栋;华北地区夏季水资源特征分析及其对气候变化的响应(Ⅱ):华北地区夏季水量丰、枯与气候变化的关系[J];高原气象;2002年01期
19 刘庆芳;2003年9月华北地区M_L≥3.0地震动态[J];国际地震动态;2003年12期
20 傅秉锋;华北地区社会经济发展对地质调查的需求分析[J];中国地质矿产经济;2003年01期
中国重要会议论文全文数据库 前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 刘峡;华北地区现今地壳运动及形变动力学数值模拟[D];中国科学技术大学;2007年
2 宫飞;华北地区结构型节水种植业模式及途径研究——以北京市顺义区为例[D];中国农业大学;2003年
3 王桂芝;中国华北地区东方蜜蜂多样性及其分类研究[D];山东农业大学;2007年
4 荣艳淑;大范围气候变化与华北干旱研究[D];南京信息工程大学;2004年
5 张利平;华北地区陆气水循环时空演化规律研究[D];武汉大学;2005年
6 温海清;金元之际的华北地方行政建置[D];复旦大学;2008年
7 张学民;华北地区壳幔剪切波速度结构及与晚新生代地质演化关系研究[D];中国地质大学(北京);2005年
8 郭锦超;近代天津和华北地区经济互动的系统研究(1880年代—1930年代)[D];南开大学;2004年
9 周连童;引起华北持续干旱的环流异常型及其与西北干旱区热力变化的关联[D];中国科学院研究生院(大气物理研究所);2007年
10 房立华;华北地区瑞利面波噪声层析成像研究[D];中国地震局地球物理研究所;2009年
中国硕士学位论文全文数据库 前10条
1 罗遵兰;华北地区冬小麦与夏玉米水分生产函数模型研究[D];中国农业大学;2005年
2 孙燕;华北地区降水异常及其与北太平洋海温、100hPa高度场异常的关系[D];南京气象学院;2002年
3 史晓娟;全面建设小康社会目标下的华北地区电力需求预测[D];华北电力大学(北京);2006年
4 郭士忠;清代华北地区粮食流通[D];山西大学;2008年
5 张艳;清代典当业税收问题[D];陕西师范大学;2007年
6 左艳艳;华北地区中老年针织毛衫感性消费研究[D];江南大学;2009年
7 范烨;华北地区气溶胶分布和层积云结构特征研究[D];南京信息工程大学;2006年
8 韩军彩;华北地区空中水汽含量的演变特征[D];南京信息工程大学;2011年
9 李坤;中国华北地区人群TAP1基因多态性与乙肝的关联性研究[D];山东大学;2005年
10 高明;北洋军阀统治时期华北地区兵变述论[D];扬州大学;2007年
中国重要报纸全文数据库 前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年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978