吉黑东部山区构造应力场模拟.pdfVIP

第九届全国矿床会议论文集 吉黑东部山区构造应力场模拟术 邵兆刚 孟宪刚 朱大岗 韩建恩 余佳 吕荣平 (中国地质科学院地质力学研究所，北京100081) 自中生代开始，东亚滨太平洋地区即已经受着太平洋板块的向西挤压，由新华夏第一隆起带连结而成 的“花彩列岛”，以及新生代早期变故生成的沟一弧一盆系统，包括中国中部和东部的陆缘构造格架，莫不 与此相关。鉴于板块俯冲机制及其对研究区内生金铜成矿分析的认识，将研究区控矿构造和研究区内洋壳 深部俯冲作用相联系，探讨其深部动力成矿机理。 1模型的建立 地质模型根据区域地质调查和地震探测资料予以确定。在垂直方向共分三层：最上部为地壳层厚28 —43km，以板块俯冲带为界分为陆壳和洋壳，其下部为岩石圈地幔，深度为60kin，最下部为地幔层，由 49．65”～42。58 727．46”～131 712．68”与130047 壳幔俯冲带达到的深度，而在研究区南北两侧俯冲带相对滞后，总体形成在成矿区突出的犁形俯冲带。为 由于模型研究重点在板块深部俯冲带对模型中控矿构造的影响，在增大深部模拟范围的同时，于重点研究 区域对控矿构造体系进行细化。 数值模拟采用麦克斯韦体模型，参考已有地壳密度、弹性模量及黏性参数。对地壳表层弹性参数，重 点考虑洋陆地壳差异，增强大洋壳弹性，而在模型中部沉积盆地分布区弹性强度适当减小。有限元计算模 简化单元划分，提高计算效率，应力变化急剧的区域单元划分得小一些，反之则增大。单元的边长应尽量 接近，以提高计算精度。 边界条件的设定：根据区域地质资料在西部边界加载东西向固定约束，在底部结点加载垂向约束使模 型上下固定，在模型东部边界加载东西向挤压应力，应力大小随深度增加，并在28—55km岩石圈地幔处 于模型计算中重点考虑偏应力作用，因此忽略模型中体力产生的静岩压力。在模型南北两侧加载有应力约 束，保证计算结果收敛。模型加载时间尺度为年。 2模型位移场计算结果 构造应力场位移计算结果表明，模型存在整体向西位移，导致模型整体挤压收缩变形；北东向区域断 裂表现为右行扭动特征，而北西向断裂带表现为左行运动特征。在垂直方向+，模型整体为隆升增厚，但不 同区块隆升幅度存在明显差异。 ·本文得到中国地质调查局项目“吉黑东部山区内生金铜矿床找矿模型研究”(编号：1212010634003)资助。 第一作者简介：邵兆刚，男，1970年生，博士，副研究员，从事区域地质、构造地质与地球物理研究。 ·424· 邵兆刚等：吉黑东部山区构造应力场模拟 东西向位移计算结果表明，模型整体发育向西的位移分量(位移负值)，但在断裂带两侧存在明显幅 度变化，北东向断裂南侧和北西向断裂北侧位移明显增加，反映区域北东向断裂右行扭动、北西向断裂左 行扭动特征。在剖面上，位移分布在壳幔俯冲带发生偏折，反映俯冲带相对向下插入俯冲运动。 南北向位移计算结果表明，由于受到南北方向约束，模型整体南北向位移较小，但在断裂发育部位出 现强烈变化，总体表现为北东向断裂带西南部存在明显向北的正向位移，往东北逐渐向往南的负向位移过 渡，但幅度有所减弱，垂向上南北向位移变化不大。 模型垂向位移反映为整体挤压增厚变形，模型表面均表现为隆升的正向位移。但不同构造部位的隆升 幅度存在较大差异。邻近壳幔俯冲带西侧存在一个隆升幅度减弱的相对凹陷区，其部位与弧后盆地发育区 对应，与板块俯冲过程的拖拽作用有关。再往西至吉黑成矿区，隆升幅度加大，出现相对高隆起特征。俯 冲带东侧则为相对隆升高值区。在剖面图中可见在整体挤压增厚特征背景下，壳幔俯冲带上部分布着相对 隆升低值区。 3模型应力场计算结果 区域构造应力场计算结果反映应力分布明显受断裂控制，沿断裂发生应力集中和突变现象。在板块俯 冲的东西挤压应力作用下(取压应力为负值)，模型整体东西方向应力分量变化较为均匀，与东西向整体 变化的位移场对应，但在地壳表层构造发育部位出现显著变化，沿北东向断裂带，存在整体低应力显示， 有利于岩体和矿体沿断裂侵位、富集，北西向断裂同样存在低应力显示，但减低幅度小于北东向断裂，同 时在两组断裂的交汇处出现最低值。受两组构造断裂影响，断裂发育区域整体出现网状低值区。 在区域东西向挤压应力作用下，模型南北向应力变化较为均匀，其分布特征总体显示随深度逐渐