摘要
汶川地震后,紫坪铺水库蓄水是否触发了汶川地震在国内外学术界引起了广泛关注.除定性讨论外,许多学者也采用定量分析的方法进行了计算,但因计算结果不同而得出了不同的结论.本文从目前紫坪铺水库蓄水不同研究组定量计算中出现的争议为出发点,通过对水库蓄水定量计算基本原理和可能引起计算结果差异可能因素的分析,找出定量计算中的关键影响因素,了解目前水库蓄水定量计算中存在的不确定性问题所在.初步结果显示:计算方法、模型维数、扩散模型、震源参数和扩散系数等的取值不同是造成计算结果差异的主要因素,特别是裂隙岩体的扩散系数.在紫坪铺水库定量计算中模型维数的差别使得汶川地震震源处的库仑应力变化计算结果相差约3倍;仅考虑断层渗透率(把岩体渗透率视为无穷大)或仅考虑均匀各向同性的岩体渗透率(忽视断层渗透率),均具有片面性;震源机制解断层走向倾角的差异,会显著影响库仑应力大小计算结果,可到达2~7倍;不同扩散系数下,孔隙压力相差可达几百倍."紫坪铺水库蓄水是否能够触发汶川大地震的发生?",鉴于目前的研究成果,库仑应力变化在kPa量级,尚不能排除触发的可能性,但得出的蓄水震源处的库仑应力变化太低,在背景构造应力场不明确的情况下,也不能确定一定有联系.在未来的工作中需有针对性的进行野外考察和室内试验,改进模型,采用高性能模拟分析计算,并在此基础上对中国和世界多个水库地震触发机制进行对比研究,探讨不同机制下水库地震触发机制特点,进一步量化分析水库地震发生的力学机制及水库对构造活动的影响和作用机理.
Since the 12 may 2008 Wenchuan earthquake,many scholars have done a lot of work discussing the impact of Zipingpu reservoir on the MS8.0earthquake,which is less than 10 km from the epicenter.Different opinions occurred in the past seven years,not only in the conclusions of qualitative analysis but also the quantitative analysis.In order to comprehensively understand the impact of Zipingpu reservoir on the MS8.0earthquake,we analyze the main influence factors by systematically summarizing these researches.Following the impoundment of reservoir,there are three responses which might correspond to two fundamental mechanisms.Firstly,rapid increase of elastic stress may increase the stress promptly upon the fault to critical state,and failure occurs shortly.Secondly,pore pressure beneath the reservoir and its adjacent regions is redistributed by the porous elastic nature of lithosphere-fracture-fault system.Based on the basic principle of quantitative calculation of these responses,we analyze different models and main parameters which may cause the differences.Then we further study the key influence factors of quantitative calculation and know about the existing uncertainty elements during the numerical simulation.The preliminary results show that the calculation methods(analytical method or numerical method),dimension of models(2-D or 3-D),diffusion model,diffusion coefficient and focal mechanism are the main factors resulting in the differences,especially the diffusion coefficient of the fractured rock mass.Firstly,without considering the geological environment,the changes of elastic stress due to the impoundment of Zipingpu reservoir can be attained using half space Boussinesq's analysis method.In the same way,the pore pressure can be conveniently attained from Biot's consolidation equation using analytical method.Though,the geological structures of Longmenshan region are very complex.High performance computation provides a powerful tool to quantitatively evaluate stress,which can conveniently take the geological features into model.In addition,when using the numerical method,we should guarantee sufficient accuracy.The results attained from analytical method can be the calibration of numerical method.Secondly,the2-D model is a plane strain approximation,which takes infinite linear loads instead of finite loads in limited area.It may exaggerate the loading effect.The change of Coulomb failure stress at the epicenter of Wenchuan earthquake attained from 2-D model is about 3times greater than that of3-D model.Thirdly,it is not reasonable only considering the fault permeability(assuming the permeability of rock mass as infinity)or the homogeneous isotropic rock mass permeability(ignoring the fault permeability).Fourthly,as for the uncertainty of the diffusion coefficient of the deep rock and faults,the difference of Coulomb failure stresses can reach several hundreds times,when selecting different diffusion coefficients.Finally,the different focal mechanisms,especially the focal depth,also could dramatically affect the change of Coulomb failure stress.The differences can research 2~7times when adopting different focal mechanisms derived from different institutes.Although there are many influence factors affecting the quantitative analysis results,there are still more in common with the work of each group.The impoundment of reservoir makes the thrust fault more stable,whereas leads the normal fault to be more dangerous.And the instantaneous pore pressure and diffusion pressure accelerate the occurrence of the earthquake.At the hypocenter of the earthquake,the values of the Coulomb failure stress change of each group are close to several kPa,we could not rule out the possibility that the Zipingpu Reservoir may trigger the 2008 Wenchuan earthquake.However,for the background stress is not clear and Coulomb failure stress change is too little,we are also not sure there must be a connection between reservoir and earthquake.In future work,we should target on the basis of field survey and indoor experiment,in order to improve the model and develop high performance simulation.Based on this calculation method,we make a comparison with more typical reservoir-triggered earthquakes with different mechanisms both domestic and abroad to discuss the characteristics and further analyze the mechanics of reservoir-triggered earthquake.
引文
Bell M L,Nur A.1978.Strength changes due to reservoir-inducedpore pressure and stresses and application to lake Oroville.J.Geophysics.Res.,83(B9):4469-4483.
Chen H Q,Xu Z P,Lee M.2008.Wenchuan Earthquake andseismic safety of large dams.Shuili Xuebao(in Chinese),39(10):1158-1167.
Chen L,Talwani P.2001.Mechanism of initial seismicity followingimpoundment of the Monticello Reservoir,South Carolina.Bulletin of the Seismological Society of America,91(6):1582-1594.
Chen L,Talwani P.1998.Reservoir-induced seismicity in China.Pure and Applied Geophysics,153(1):133-149.
Chen P,Wang J D.2007.Wavelet estimation of the diffusioncoefficient in time dependent diffusion models.Science in ChinaSeries A:Mathematics,50(11):1597-1610.
Chen Y.2009.Did the reservoir impoundment trigger the Wenchuanearthquake?.Science in China Series D:Earth Sciences,52(4):431-433.
Chen Y,Huang T F,Liu N R.2009.Rock Physics(in Chinese).Hefei:Press of University of Science and Technology of China.
Cheng H H,Zhang H,Zhu B J,et al.2012.Finite elementinvestigation of the poroelastic effect on the XinfengjiangReservoir-triggered earthquake.Science in China Series D:Earth Sciences,55(12):1942-1952.
Cheng W Z,Zhang Z W,Ruan X.2010.Analysis of seismicity andits causes in Zipingpu reservior region at each water storagestage.Progress in Geophys.(in Chinese),25(3):759-767,doi:10.3969/j.issn.1004-2903.2010.03.003.
Deng K,Zhou S,Wang R,et al.2010.Evidence that the 2008Mw7.9Wenchuan Earthquake Could Not Have Been Induced bythe Zipingpu Reservoir.Bull.Seism.Soc.Am.,100(5B):2805-2814.
Ding Y Z.1989.The Reservoir Induced Earthquake(in Chinese).Beijing:Seismological Press.
Gahalaut K,Gahalaut V K.2010.Effect of the Zipingpu reservoirimpoundment on the occurrence of the 2008 Wenchuanearthquake and local seismicity.Geophysical JournalInternational,183(1):277-285.
Gahalaut K,Hassoup A.2012.Role of fluids in the earthquakeoccurrence around Aswan reservoir,Egypt.Journal ofGeophysical Research,117(B2):B02303,doi:10.1029/2011JB008796.
Ge S M,Liu M,Lu N,et al.2009.Did the Zipingpu Reservoirtrigger the 2008 Wenchuan earthquake?.Geophysical ResearchLetters,36(20):L20315,doi:10.1029/2009GL040349.
Gong G Y,Xie Y D.1989.A review of laboratory investigation ofrock permeability.Chinese Journal Rock Mechanics andEngineering(in Chinese),8(3):219-227.
Gong G Y,Xie Y D.1991.Research on the diffusion of pore fluidand in-situ hydraulic diffusivity in the epicentral region ofXinfengjiang reservoir earthquakes.Acta Seismologica Sinica(in Chinese),13(3):364-371.
Gough D I,Gough W L.1970.Stress and deflection in thelithosphere near Lake Kariba—I.Geophys.J.Roy.,21(1):65-78.
Gupta H K.2002.A review of recent studies of triggeredearthquakes by artificial water reservoirs with special emphasison earthquakes in Koyna,India.Earth-Science Reviews,58(3-4):279-310.
Hardebeck J L,Nazareth J J,Hauksson E.1998.The static stresschange triggering model:Constraints from two southernCalifornia aftershock sequences.J.Geophys.Res.,103(B10):24427-24437.
Harris R A.1998.Introduction to special section:Stress triggers,stress shadows,and implications for seismic hazard.J.Geophys.Res.,103(B10):24347-24358.
Heaton T H.1975.Tidal triggering of earthquakes.Geophys.J.Int.,43(2):307-326.
Jaeger J C,Cook N G W,Zimmerman R.2007.Fundamentals ofRock Mechanics.New York:John Wiley&Sons.
Jing H,Yu X,Zhang H,et al.2010.Numerical analysis onthermal conductivity of poly-mineral rock.Earthquake Science,23(3):223-232.
King G C P,Stein R S,Lin J.1994.Static stress changes and thetriggering of earthquakes.Bulletin of the Seismological Societyof America,84(3):935-953.
Klose C D.2012.Evidence for anthropogenic surface loading astrigger mechanism of the 2008 Wenchuan earthquake.Environmental Earth Sciences,66(5):1439-1447.
Lei X L.2010.Possible roles of the Zipingpu Reservoir in triggeringthe 2008 Wenchuan earthquake.Journal of Asian EarthSciences,40(4):844-854.
Lei X L,Ma S L,Wen X Z,et al.2008.Integrated analysis ofstress and regional seismicity by surface loading—A case studyof Zipingpu reservoir.Seismology and Geology(in Chinese),30(4):1046-1064.
Liu G P,Fu Z X.2001.A study on the great gulang earthquaketriggered probably by static stress change resulting from thegreat haiyuan earthquake.Chinese J.Geophys.(in Chinese),44(z1):107-115,doi:10.3321/j.issn:0001-5733.2001.z1.014.
LüJ,Su J R,Jin Y K,et al.2008.Discussion on relocation andseismo-tectonics of the Ms8.0Wenchuan earthquake sequences.Seismology and Geology,30(4):917-925.
Ma W T,Xu X W,Yu G H,et al.2012.Assess the reservoirinduced seismic hazard in the hubei section of the three gorgesreservoir using gray clustering method.Seismology andGeology(in Chinese),34(4):726-738.
Ma W T,Xu C P,Zhang X D,et al.2011.Study on therelationship between the reservoir-induced seismicity atzipingpu reservoir and the Ms8.0 wenchuan earthquake.Seismology and Geology(in Chinese),33(1):175-190.
Pauchet H,Rigo A,Rivera L,et al.1999.A detailed analysis ofthe February 1996aftershock sequence in the eastern Pyrenees,France.Geophys.J.Int.,137(1):107-127.
Perfettini H,Schmittbuhl J.2001.Periodic loading on a creepingfault:Implications for tides.Geophysical Research Letters,28(3):435-438.
Rice J R,Cleary M P.1976.Some basic stress diffusion solutionsfor fluid-saturated elastic porous media with compressibleconstituents.Reviews of Geophysics,14(2):227-241.
Roeloffs E A.1988.Fault stability changes induced beneath areservoir with cyclic variations in water level.J.Geophys.Res.,93(B3):2107-2124.
Shen L Y,Chang B Q,Wu M B.1992.Application of the stresspore pressure coupled theory of porous medium to the researchof inducing earthquakes at Xinfengjiang Reservoir.EarthquakeResearch China(in Chinese),8(3):36-43.
Simpson D W.1976.Seismicity changes associated with reservoirloading.Eng.Geol.,10(2-4):123-150.
Simpson D W,Leith W S,Scholz C H.1988.Two types ofreservoir-induced seismicity.Bull.Seism.Soc.Am.,78(6):2025-2040.
Stein R S,Lisowski M.1983.The 1979 Homestead Valleyearthquake sequence,California:Control of aftershocks andpostseismic deformation.J.Geophys.Res.,88(B8):6477-6490.
Stein R S,King G C P,Lin J.1992.Change in failure stress on thesouthern San Andreas fault system caused by the 1992magnitude=7.4 Landers earthquake.Science,258(5086):1328-1332.
Sun Y J,Kim I K,Zhang H,Shi Y L.2013.Simulation of tensiletest of the steel fiber reinforced concrete with the finite elementmethod.Earthquake(in Chinese),33(4):145-152.
Sun Y J,Zhang H,Dong S W,et al.2012.Study on effect of theZipingpu reservoir on the occurrence of the 2008 Wenchuanearthquake based on a 3D-poroelastic model.Chinese J.Geophys.(in Chinese),55(7):2353-2361,doi:10.6038/j.issn.0001-5733.2012.07.020.
Talwani P.1997.On the nature of reservoir-induced seismicity.Pure and Applied Geophysics,150(3-4):473-492.
Talwani P,Acree S.1984.Pore pressure diffusion and themechanism of reservoir-induced seismicity.Pure and AppliedGeophysics,122(6):947-965.
Talwani P,Chen L,Gahalaut K.2007.Seismogenic permeability,ks.J.Geophys.Res.,112(B7):B07309,doi:10.1029/2006JB004665.
Tao W,Masterlark T,Shen Z K,et al.2014.Triggering effect ofthe Zipingpu Reservoir on the 2008 Mw7.9 Wenchuan,China,Earthquake due to poroelastic coupling.Chinese J.Geophys.(in Chinese),57(10):3318-3331,doi:10.6038/cjg20141019.
Toda S,Stein R S,Reasenberg P A,et al.1998.Stress transferredby the1995 Mw=6.9 Kobe,Japan,shock:effect onaftershocks and future earthquake probabilities.J.Geophys.Res.,103(10):24543-24565.
Vidale J E,Agnew D C,Johnston M J S,et al.1998.Absence ofearthquake correlation with Earth tides:An indication of highpreseismic fault stress rate.J.Geophys.Res.,103(B10):24567-24572.
Wang C Y,Manga M.2010.Earthquakes and Water.New York:Springer.
Wang W M,Zhao L F,Li J,et al.2008.Rupture process of theMs8.0 Wenchuan earthquake of Sichuan,China.Chinese J.Geophys.(in Chinese),51(5):1403-1410.
Wu J L.2001.Mechanics of Elasticity(in Chinese).Beijing:ChinaHigher Education Press.
Zhang B,Shi Y L.2010.A discussion on the influences of Zipingpureservoir on the stability of faults in the neighborhood.Journalof the Graduate School of the Chinese Academy of Sciences,27(6):754-760.
Zhang Y,Feng W P,Xu L S,et al.2009.Spatio-temporal ruptureprocess of the 2008 great Wenchuan earthquake.Science inChina Series D:Earth Sciences,52(2):45-154.
Zhang Y,Wang L F,Wang J J,et al.2012.Electromagnetic surveyof the seismogenic structures beneath the Longtan reservoir inGuangxi Province.Chinese J.Geophys.(in Chinese),55(4):1400-1410,doi:10.6038/j.issn.0001-5733.2012.04.036.
Zheng L,Zhang H,Sun Y J,et al.2013.Comparison of 2Dand 3Dfinite element simulations in the research of reservoir inducedearthquakes.Earthquake(in Chinese),33(4):162-171.
Zhou B,Xue S F,Deng Z H,et al.2010.Relationship between theevolution of reservoir-induced seismicity in space-time and theprocess of reservoir water body load-unloading and waterinfiltration—A case study of Zipingpu reservoir.Chinese J.Geophys.(in Chinese),53(11):2651-2670,doi:10.3969/j.issn.0001-5733.2010.11.013.
Zhou L Q,Zhao C P,Chen Z L,et al.2012.Three-dimensional VPand VP/VS structure in the Longtan reservoir area by localearthquake tomography.Pure and Applied Geophysics,169(1-2):123-139.
Zhou S,Deng K.2011.Reply to'Comment on'Evidence that the2008 Mw7.9Wenchuan Earthquake could not have been inducedby the Zipingpu Reservoir'by Kai Deng,Shiyong Zhou,RuiWang,Russell Robinson,Cuiping Zhao,and Wanzheng Cheng"by Shemin Ge.Bull.Seism.Soc.Am.,101(6):3119-3120.
Zhu B,Liu C,Shi Y L,et al.2011.Apllication of flow driven porenetwork crack model to Zipingpi reservoir and Longmenshanslip.Science China(Physics,Mechanics&Astronomy),54(8):1532-1540.
Ziv A,Rubin A M.2000.Static stress transfer and earthquaketriggering:No lower threshold in sight?.J.Geophys.Res.,105(B6):13631-13642.
陈厚群,徐泽平,李敏.2008.汶川大地震和大坝抗震安全.水利学报,39(10):1158-1167.
陈颙.2009.汶川地震是由水库蓄水引起的吗?.中国科学D:地球科学,39(3):257-259.
陈颙,黄庭芳,刘恩儒.2009.岩石物理学.合肥:中国科学技术大学出版社.
程惠红,张怀,朱伯靖等.2012.新丰江水库触发地震的孔隙弹性耦合有限元模拟.中国科学D:地球科学,42(6):905-916.
程万正,张致伟,阮祥.2010.紫坪铺水库区不同蓄水阶段的地震活动及成因分析.地球物理学进展,25(3):759-767,doi:10.3969/j.issn.1004-2903.2010.03.003.
丁原章.1989.水库诱发地震.北京:地震出版社.
龚钢延,谢原定.1989.岩石渗透率变化的实验研究.岩石力学与工程学报,8(3):219-227
龚钢延,谢原定.1991.新丰江水库地震区内孔隙流体扩散那与原地水力扩散率的研究.地震学报,13(3):364-371.
雷兴林,马胜利,闻学泽等.2008.地表水体对断层应力与地震时空分布影响的综合分析———以紫坪铺水库为例.地震地质,30(4):1046-1064.
刘桂萍,傅征祥.2001.海原大地震对古浪大地震的静应力触发研究.地球物理学报,44(增刊):107-115,doi:10.3321/j.issn:0001-5733.2001.z1.014.
吕坚,苏金蓉,靳玉科等.2008.汶川8.0级地震序列重新定位及其发震构造初探.地震地质,30(4):917-925.
马文涛,徐锡伟,于贵华等.2012.使用灰色聚类方法评估长江三峡水库湖北不同库段水库诱发地震的震级上限.地震地质,34(4):726-738.
马文涛,徐长朋,张新东等.2011.紫坪铺水库与汶川地震关系的讨论.地震地质,33(1):175-190.
沈立英,常宝琦,吴名彬.1992.弹性多孔介质中应力-孔隙压耦联理论在新丰江水库地震研究中的应用.中国地震,8(3):36-43.
孙玉军,Kim I K,张怀等.2013.利用有限元方法模拟力学拉伸试验———以钢纤维混凝土为例.地震,33(4):145-152.
孙玉军,张怀,董树文等.2012.利用三维孔隙弹性模型探讨紫坪铺水库对汶川地震的影响.地球物理学报,55(7):2353-2361,doi:10.6038/j.issn.0001-5733.2012.07.020.
陶玮,Masterlark T,沈正康等.2014.紫坪铺水库造成孔隙弹性耦合变化及其对2008年汶川地震触发作用.地球物理学报,57(10):3318-3331,doi:10.6038/cjg20141019.
王卫民,李连锋,李娟等.2008.四川汶川8.0级地震震源过程.地球物理学报.51(5):1403-1410.
吴家龙.2001.弹性力学.北京:高等教育出版社.
张贝,石耀霖.2010.紫坪铺水库对附近断层稳定性影响的探讨.中国科学院研究生院学报,27(6):754-760.
张勇,冯万鹏,许力生等.2008.2008年汶川大地震的时空破裂过程.中国科学D:地球科学,38(10):1186-1194.
詹艳,王立凤,王继军等.2012.广西龙滩库区深部孕震结构大地电磁探测研究.地球物理学报,55(4):1400-1410,doi:10.6038/j.issn.0001-5733.2012.04.036.
郑亮,张怀,孙玉军等.2013.水库触发地震研究中二维与三维有限元模拟结果比较.地震,33(4):162-171.
周斌,薛世峰,邓志辉等.2010.水库诱发地震时空演化与库水加卸载及渗透过程的关系———以紫坪铺水库为例.地球物理学报,53(11):2651-2670,doi:10.3969/j.issn.0001-5733.2010.11.013.