水库蓄水应力场模拟孕震机理.doc

水库蓄水应力场模拟孕震机理 水库蓄水应力场模拟孕震机理预读: 摘要:最新杂志：武汉公安干部学院学天水行政学院学报天津市社会主义学院四川行政学院学报四川文化产业职业学四川统一战线四川省公安管理干部世纪风采时事资料手册生命与灾害上海支部生活陕西社会主义学院学三峡水库自2003年6月第一次蓄水以来,目前已成功实现了172m水位的试验性蓄水,标志着这个世界最大水利枢纽工程的防洪、发电、通航等功能均达到了设计要求,综合效益得到了充分发挥.然而自2003年6月蓄水以来,在三峡库区共发生近7000次地震,而且绝大多数M1.0以下地震属于非构造成因的水库诱发地震(车用太等,2009).然而,2008年9月第3期蓄水前后,2008年9月27日和11月22日在三峡库首区分别发生M3.2与MS4.1地震,MS4.1地震是三峡水库蓄水以来在库首区内发生的最大地震.陈蜀俊等(2006)和吴建超等(2009)曾对三峡水库蓄水后全位移场的变化进行了数值模拟,但仅局限于分析上地壳表层(Z=－1.5m)的位移场,没有讨论水库诱发地震更深层次的动力成因.本文从库首区上地壳应力场尤其是等效应力场分析的角度对水库诱发地震震例及地震活动性进行分析,对三峡水库诱发地震机理研究进行更进一步的探索.1胡家坪MS4.1地震概况2008年11月22日16时01分胡家坪MS4.1地震发生,震中位于宜昌市秭归县屈原镇胡家坪(31°N,110.8°E),震源深度约为8km,震中烈度为Ⅵ度.震源机制解揭示出震源破裂为NNW向断裂的逆冲左旋走滑错动,P轴压应力方位角为287°,仰角为33°,主要受压应力作用(图1).这次地震的有感范围较大,西自巴东官渡口镇,东至宜昌夷陵区,北自兴山县,南至秭归县杨林镇,总面积约400km2.车用太等(2009)研究认为,胡家坪MS4.1地震属构造型水库诱发地震,是第三次水库蓄水引起的水体荷载加剧了断层面上的剪应力作用,同时增大了断层面上的正应力作用,最终导致库仑破裂应力增大,从而引发了地震.田宗伟(2008)认为,本次地震属三峡地区正常背景条件下发生的构造地震.然而,他们都是基于有关地质资料和现场地震考察资料来进行分析讨论,并未从定量角度对蓄水导致的应力变化进行计算.因此,十分有必要对蓄水所导致的库首区孕震动力环境进行数值计算,从上地壳深处研究本次地震的发震背景(王博等,2008;廖武林等,2009).2库首区有限元模型的建立2.1几何参数模型本文所确定的库首区研究范围为110km×85km,厚度取至上地壳底部埋深12km处.因此,该几何模型可以看做是长轴沿NW—SE方向延伸的六面体(图2).以(31°56’35″N,110°00’00″E)为坐标原点,X轴方向为SE33°,Y轴方向为NE33°,Z轴垂直于水平面,向上为正.Z=0相当于海拔高程50m.2.2力学参数模型力学参数模型按地层岩性可划分为沉积盖层和基底构造层:第一层厚约7km,主要岩石类型为花岗岩、灰岩和砂岩;第二层厚约5km,主要岩石类型为花岗岩.考虑到其内部详细岩石类型及参数对模拟计算结果的影响,该模型共划分了7层,各层底面Z轴坐标自上而下分别为:－1.5、－3、－4、－5、－7和－12km,各层分别对应不同的岩石类型及参数(表1).本文采用连续介质力学方法对断层模型进行研究,即断层内为弱介质夹层,其厚度不为零,也不能忽略.为模拟断层应力场的变化途径和变化趋势,取断裂带的材料为弹塑性材料.依据前人研究及本区工程研究中得出的岩石力学参数,本文统计得出断裂带的参数分别为:弹性模量E=1.07×1010Pa,泊松比ν=0.24,密度ρ=2570kg/m3.2.3边界条件模型根据Wang等(2001)利用GPS数据获得的中国大陆水平方向速度场可知,包括研究区在内的整个华南块体总体向SE方向运动,地幔流总体向东运移、牵引.因此,首先对模型的NW向侧面作X向约束.三峡库首区的黄陵背斜为一古老的刚性地质体,根基深深扎入下地壳,为下元古界—上元古界基底,并在漫长的地质历史中表现为隆起和对区域应力作用产生强大的反作用力.因此,模型对黄陵背斜核部作X向和Y向约束.对模型的底面,即Z=－12km处作Z向的约束,模型顶面及其它侧面为自由面,不作任何约束.依据陈蜀俊等(2005,2006)和吴建超(2010)的研究,本文选用的加载方式为:NE50°方向压应力大小为790MPa,NW320°方向张应力大小为540MPa.采用三维线弹性模型,除作了约束的部位外,每个节点有三个自由度.以满足计算精度为要求作有限元网格划分,对断裂带附近和重点研究区域,例如巴东库段和香溪库段的网格再做细化.本次计算模型单元数共计17146个,节点数共计20264个.3蓄水达172m水位后等效应力场的数值模拟三峡水库蓄水后,提升了长江上游水位,并淹没了长江部分支流河段,库水影响范围内的水文地质条件和岩石物性参数必将发生变化.本文着重考虑了三峡水库蓄水达172m水位后,水体荷载作用、应力腐蚀作用和孔隙水压作用的影响.依据蓄水后GPS形变场、重力场实测结果对模型进行了调整,计算得到了蓄水前后三峡库首区等效应力场的三维数值模拟结果.3.1等效应力的含义弹塑性力学中,等效应力定义为在一般应力状态下的各应力分量经适当的组合而形成的与单向应力等效的应力,也称比较应力、统一应力和应力强度(陈明祥,2007;张林洪等,2002,2004),等效应力用主应力表示为式中,σ1,σ2,σ3分别为最大、中间、最小主应力.朱守彪和张培震(2009)对汶川地震的发生过程进行了有限元的数值模拟,讨论了等效应力与地震活动的关系,结果显示等效应力在空间分布上由地震前的分散状态逐渐向龙门山地区集聚并增大,空间上龙门山断裂带成为了等效应力积累单元.应力水平不断升高,直至超过摩擦极限,最终导致了汶川地震的发生.因此,本文选择三峡库首区上地壳深处的等效应力场作为研究对象,探讨蓄水后等效应力场的变化及其与胡家坪MS4.1地震的关系.3.2数值模拟结果分析本文的研究重点是上地壳深部等效应力场的变化,考虑到胡家坪MS4.1地震的震源深度为8km,因此,取Z=－8km作为输出层,得到蓄水达172m水位后等效应力场的结果(图3).三峡水库172m水位试验性蓄水之后,地壳深处8km附近沿新滩、归州和周坪三地出现了较大范围的等效应力高值区,等效应力值约为900MPa,长轴近似呈NE向.分析该等效应力高值区附近的地质构造,笔者认为该等效应力高值区与九畹溪断裂、仙女山断裂及新华—水田坝断裂的围限作用有关,同时还可能与该区域的库水渗透作用有密切关系.为了分析数值模拟结果的可靠性,笔者选择了重力实测资料作为验证(车用太等,2009).三峡水库蓄水达172m水位后,在香溪—茅坪之间水库北岸存在重力正异常,异常幅度达30×10－8m/s2(图4)(车用太等,2009).该重力场的异常区展布方向与水库走向一致.将数值模拟得到的蓄水后等效应力结果与重力实测结果对比,发现蓄水达172m水位后,库首区沿着新滩、归州和周坪三地出现了较大范围的等效应力高值区,不考虑边界效应和断裂端部的影响,这一等效应力高值区是研究范围内唯一比较明显的区域.等效应力高值区(图3)位于本文三维模型的中部,受底面边界条件的影响较小,模拟结果较为准确可靠.该区域与重力实测得到的异常区基本对应,且均分布于长江两岸约10km范围内.3.3等效应力的差值变化与胡家坪MS4.1地震关系等效应力的差值变化分析,集中反映了直接由蓄水引起的等效应力的变化及其分布特征,是蓄水引起库首区孕震环境变化的有力证据之一.为了得到三峡水库蓄水后等效应力场的差值变化,笔者将本文数值计算的172m水位与蓄水前的等效应力值相减,并采用局部多项式插值方法,得到了172m水位时相对蓄水前的等效应力差值图(图5).由图5可以发现,172m水位相对于蓄水前的等效应力差在巴东、归州和新滩镇出现了3个正高极值区,其中以新滩极值区的范围和量级最大,呈NE—SW向延伸,极值区的等效应力值由中心向两侧逐渐减小.香溪正高极值区的南东侧为一负低极值区,呈近SN向延伸,量值较小.三峡水库蓄水达172m水位后,胡家坪MS4.1地震就发生在正负异常的边界带上,也即香溪正高极值区和秭归负低极值区的过渡带上.此过渡带等效应力变化明显,易于应力和应变能的集中,加之此过渡带处于仙女山断裂北端段,容易促使断裂端部失稳错动而引发地震活动.4胡家坪MS4.1地震的孕震机理胡家坪MS4.1地震的震区,位于鄂西黄陵背斜与秭归向斜之间的沉积地层发育区.地震的宏观震中位于三叠系上统的灰岩地层,岩溶较为发育且渗透性较强.如图6所示,香溪库段ML≥2.0地震破裂面走向较为集中,主要为NW向和NNE向,与香溪库段的仙女山断裂及九畹溪断裂呈现了较好的对应关系.笔者研究认为,胡家坪MS4.1地震发生在蓄水后等效应力高值区与东部大范围低值区的过渡带上,等效应力差值约为400MPa.该过渡带属于等效应力变化最为剧烈、最易产生应变能积累和应力集中的区域.蓄水后引起等效应力场的变化的原因可能是库水荷载、下渗腐蚀背景下引起的水文地质条件变化及上地壳表层应力状态调整的结果.三峡水库172m水位试验性蓄水进一步加剧了断面上的压应力作用,同时库水下渗导致了断裂带上孔隙水压显著增大,最终导致库仑破裂应力增大,引发了胡家坪MS4.1构造型水库诱发地震.5结语(1)利用有限元方法计算得到了蓄水达172m三峡库首区8km深处的等效应力场,与流动重力场异常变化图像基本一致.沿新滩、归州和周坪三地出现了等效应力高值区,反映了172m水位蓄水造成的上地壳等效应力场变化的基本状况.(2)172m水位相对于蓄水前的等效应力差在巴东、归州和新滩镇出现了3个正高极值区,蓄水达172m水位后,胡家坪MS4.1地震发生在香溪正高极值区和秭归负低极值区的过渡带上.这样的地段等效应力变化明显,易于应力和应变能的集中,加之处于仙女山断裂北端段,容易促使断裂失稳错动而引发地震活动.(3)胡家坪MS4.1地震可能是库水荷载、下渗腐蚀背景下引起的水文地质条件变化的结果.三峡水库172m水位试验性蓄水进一步加剧了断面上的压应力作用,同时库水下渗导致了仙女山断裂带上孔隙水压显著增大,最终导致库仑破裂应力增大,引发了胡家坪MS4.1构造型水库诱发地震.。