煤岩动力灾害现象简介

1.1.1 煤岩动力灾害现象

煤岩动力灾害现象是煤岩体在外力作用下发生的具有动力效应和灾害后果的一种灾害现象。它具有范围广、形式多样的特点，如地震、火山爆发、隧道失稳、山体滑坡、煤与瓦斯突出、冲击矿压、顶板事故等均属于煤岩动力灾害现象。煤岩动力灾害现象危害巨大，对人们的生命和财产均会造成巨大的伤亡与损失。

面对飓风、旱灾、洪灾、地震、山体滑坡和雪崩等频发的自然灾害，印度建立了颇具特色的灾害管理体制。在国家、邦、县和区一级均有统一的灾害管理机构。印度还制定了一系列灾前备灾和部门发展计划，以及飓风地震等灾后重建计划［1］。

（1）地震［2～8］

众所周知，地球时时刻刻都处在运动变化之中。地球的运动变化会产生巨大的力，使地下的岩层发生变形。起初，变形很缓慢；但当受到的力太大，岩层不能承受时，就会突然破裂；岩层破裂所产生的振动传到地表就会引起地表的振动，这就是地震。一般情况下，变形区域越长、越宽，释放的能量就越多，构造地震的震级也将越大。如图1.1为一卫星图片，可以清晰地看到地震产生的断层，图1.2所示为1906年旧金山地震遗址照片，可以看到地震造成的破裂错断了一排原本连在一起的篱笆。

图1.1 地震断层卫星图片

图1.2 1906年旧金山地震遗址照片

地震灾害是群灾之首，它具有突发性和不可预测性，以及频度较高，并产生严重次生灾害，对社会也会产生很大影响等特点。强烈地震发生时会造成建筑物破坏，同时引起山崩、滑坡、泥石流、地裂、地裂缝、喷砂、冒水等地表的破坏和海啸。因地震的破坏而引起的一系列其他灾害，包括火灾、水灾和煤气、有毒气体泄漏，细菌、放射物扩散、瘟疫等对生命财产造成的灾害。

《竹书纪年》是迄今世界上对地震进行文字记载的最早书籍，书中将地震称为“泰山震”，记录的是公元前1830年发生的一次地震现象。全球每年发生地震500多万次！但只有不足2%的地震能被人们感觉到，其中每年100多个破坏性地震，给人类造成了巨大经济损失和人员伤亡。根据美国地震调查局国家地震信息中心（USGS NEIC）的互联网数据，2001年全世界因为地震死亡人数达到2.1万多人，其中死亡人数最多的一次地震是1月26日发生在印度古吉拉特邦的里氏7.7级大地震，导致2万人死亡。强烈地震往往以其猝不及防的突发性和巨大的破坏力给社会经济发展、人类生存安全和社会稳定、社会功能带来严重的危害。在各种自然灾害造成的损失中，地震损失占52%，是名副其实的“群害之首”。

中国位于世界两大地震带——环太平洋地震带与欧亚地震带之间，受太平洋板块、印度板块和菲律宾海板块的挤压，地震断裂带十分发育，主要地震带有23条。在刚刚过去的20世纪里，中国共发生6级以上地震近800次，遍布除贵州、浙江两省和香港特别行政区以外所有的省、自治区、直辖市，是世界上蒙受地震灾害最为严重的国家之一。在当前全球城市化进程急剧发展，尤其是我国城市数量剧增的情况下，更为明显。就我国国内的各种自然灾害而言，在地震灾害、气象灾害、海洋灾害、地质灾害、农作物生物灾害、林业灾害等各大灾种中，据新中国成立以来50多年来的资料统计，以各种灾害所造成的经济损失比例来看，气象灾害为57%，居灾害之首，但就各种自然灾害的人口死亡的统计来看，地震灾害占54%，为群害之首。

世界上地震的频繁发生以及造成的巨大损失给地震预测预报提出了严峻的研究课题，但从世界范围来看，目前地震预报仍然处于探索阶段，主要是根据多年积累的地震观测资料进行经验性的预测预报，在地震孕育机理发生发展规律上尚未完全清楚。因此，当今的地震预测预报具有很大的局限性。地震学家们对不同区域进行地震危险性分析，并对该地区未来一定时间内发生地震的可能性进行预测，划分危险区，对发生地震后可能造成的损失进行预测。在地震孕育过程和前兆机理方面的研究，自20世纪70年代以来有了很大进展，如美国和前苏联提出了膨胀-扩张模式和雪崩-不稳定模式，我国提出过震源组合模式和膨胀蠕动模式等。吴立新等［6～7］以构造地震孕震机制之二即断层粘滑发震为模拟对象，利用双轴加载实验系统和红外热像仪，对4类断层组合条件下双剪粘滑过程中的红外辐射温度场的时空演变特征进行了模拟实验研究。研究表明：断层表面温度场的时空演变不仅与应力水平有关，还与组成摩擦面的两侧岩性及其粗糙度有关。具体表现在：①应力集中和摩擦作用强的区域，其红外辐射较强；而应力松弛和摩擦作用弱的区域其红外辐射较弱。②当两断层岩性及摩擦面条件对称时，红外热像呈双蝶翼形；当摩擦面粗糙，红外热像呈非对称、非均衡时空演变特征，如串珠状、针状、倒针状、条带状、单蝶翼形或其依次变化。卫星热红外遥感与差分干涉雷达（D-INSAR）技术相结合，并以活动断层预滑及其后续粘滑为监测重点，可望成为构造地震遥感监测和地震短临时-空-强预报的新途径。这些研究虽然在一定程度上解释了地震孕育和前兆机理，但是由于地震的极其复杂性，各种理论仍不能圆满解释各种前兆现象的特征。

在防震减灾的实践中，人们逐渐认识到防震减灾涉及社会经济生活的各个方面和广大民众的生活，是一项复杂的社会系统工程。要有效地减轻地震灾害，一方面必须在政府领导下，普及地震科学知识、提高全民的防震减灾意识，依靠法制、依靠科技，动员、组织各方面的力量进行协调一致的努力；另一方面科研工作者要致力于地震预测预报的机理研究和综合预报的研究，从经验性预报向物理性预报方法转变，提高地震预报的准确性。

（2）火山喷发［8～9］

地壳之下100～150 km处，有一个“液态区”，区内存在着高温、高压下含气体挥发分的熔融状硅酸盐物质，即岩浆。岩浆从地壳薄弱的地段冲出地表，就形成了火山。火山爆发是地球最具威力的自然现象之一，释放出的能量有时比核爆炸还要大数千倍。

在火山喷出的物质中，既有固体物质如碎的岩石块、碎屑和火山灰粉，也有液体物质如熔融岩流、水以及各种泥流等，还夹杂有水蒸气和C、H、N、F、S的氧化物等气态物质，有时能喷射出可见或不可见的光、电、磁、声和放射性物质等，因而会致人死亡、使得电器仪表失灵、导致飞机轮船失事等事故。在过去的500年里，大约有几百万人死于火山之灾。人类历史上伤亡最大的一次火山发生于公元前1470年欧洲爱琴海的希腊桑托林岛，该火山喷发时喷出多达625亿立方米的物质，毁灭了一个名叫米诺斯的人类文明。

火山喷发虽然是突发的，但是也有其内在规律性，前兆比地震还要明显，人们可以预先得知而加以逃避。为了更加及时、准确地预测预报火山爆发，科学家们对此进行了不懈的努力，并成功地预报了一些火山爆发。如1979年美国在圣海伦斯山周围设置了观测站，从而成功地对圣海伦斯山火山爆发进行了预报。

（3）山崩、滑坡［10～11］

山崩通常是指土石或岩石在重力作用或因水的润滑而发生的快速下滑移动或坠落，是一种典型的块体快速运动形式。山崩会造成边坡失稳，产生的落石常会摧毁道路两旁的建筑，掩埋房屋、阻塞交通等事故，严重时还会击中过往车辆造成人员伤亡。

滑坡则是仅仅次于地震和洪水的一种严重地质灾害，不仅其出现的频率和广度远远大于地震事件，而且造成的损失也很巨大，并且会危及人们的生命财产安全。如表1.1为部分国家和地区的年平均滑坡灾害死亡人数。

表1.1 滑坡造成的平均死亡人数及年概率（20世纪末）

在我国，由于70%的地方属于山区，因此，滑坡的发生密度大，频率高，从而成为世界上受滑坡危害最严重的国家之一。迄今为止，有滑坡灾害报道的省市包括：南京、西安、宝鸡、延安、重庆、兰州、台湾以及香港等，随着经济的迅速发展，由于不合理的开发利用土地资源，危及人类安全的滑坡问题将更加突出和严重。因此，尽量减少对地质环境的破坏，治理边坡，预防边坡滑坡事故，对减少滑坡灾害造成的人们生命财产损失具有非常重要的实际意义。

（4）煤与瓦斯突出［12］

煤与瓦斯突出是矿井含瓦斯煤岩体在压力作用下从煤岩体中向采掘空间急剧运动并伴随着大量瓦斯气体喷出的一种煤岩动力灾害现象。煤岩动力灾害现象也是本书要研究的主要对象。

按动力现象力学特征可分为突出、压出和倾出，其主要作用力包括：地应力、瓦斯压力和重力。

按照动力现象的强度分类：

1）小型突出：强度小于50 t/次（突出后，经过几十分钟后瓦斯浓度可以恢复正常）；

2）中型突出：强度50～99 t/次（突出后，经过一个班后瓦斯浓度可以恢复正常）；

3）次大型突出：强度100～499 t/次（突出后，经过一天以上瓦斯浓度逐步恢复正常）；

4）大型突出：强度500～999 t/次（突出后，经过几天后瓦斯浓度可以恢复正常）；

5）特大型突出：强度大于1000 t/次（突出后，经过长时间瓦斯浓度才恢复正常）。

这种强大的动力灾害现象，给煤矿安全特别是井下工作人员的生命财产造成了极其严重的威胁。世界上最大的突出发生在1969年7月13日苏联顿巴斯加加林矿井-710 m水平石门揭煤过程中，突出煤量为14000 t，突出瓦斯量大于250万m3。

我国的煤与瓦斯突出主要是煤与甲烷突出，有4处矿井也曾发生过30多次煤、岩与二氧化碳突出，并且不同地区突出现象具有不同的特点。发生在我国最大的突出是1975年8月8日在四川天府矿务局三汇坝一矿井主平硐在用震动性放炮揭穿6号煤层时发生的，共突出煤岩12780 t，瓦斯140万m3。

（5）冲击矿压［13～17］

和煤与瓦斯突出相类似，冲击矿压（又称冲击地压、岩爆或矿震），也是一种强烈的动力灾害现象，只是没有瓦斯气体参与，其发生通常是由于煤、岩体内部应力达到一定程度后，聚集在其中的能量以急速、猛烈的形式释放，造成煤、岩体的破坏以及支架、设备、井巷的破坏，严重时会造成人员伤亡事故发生。

有记载的冲击矿压发生在1738年的英国南史塔福煤田。其后，德国、波兰、加拿大、前苏联、日本、法国以及我国相继发生了多起冲击矿压事故。我国冲击矿压事故产生的灾害十分严重。目前，我国有近50对矿井累计发生过冲击矿压4000多次，造成数以百计的人员伤亡，巷道破坏达30多公里。

冲击矿压具有突发性、瞬时震动性以及巨大的破坏性等三个明显的显现特征。根据煤岩体内部的应力状态不同，冲击矿压可以分为：重力型冲击矿压、构造应力型冲击矿压以及中间型或重力-构造型冲击矿压。根据冲击强度，冲击矿压可以分为：弹射、矿震、弱冲击和强冲击。根据震级的强度和煤岩体内部抛出的煤岩量，冲击矿压可以分为三级：轻微冲击（1级，指抛出煤岩量在10吨以下，震级也在1级以下）、中等冲击（2级，指抛出煤岩量在10～50 t，震级也在1～2级）、强烈冲击（3级，指抛出煤岩量在50 t以上，震级也在2级以上）冲击矿压。若根据冲击矿压发生的地点和位置又可以分为煤体冲击和围岩冲击两种。

冲击矿压的产生机理非常复杂，国内外学者基于现场实际调查和实验室研究结果，从不同角度对此进行了研究，提出了以下理论来解释冲击矿压产生的机理：强度理论、刚度理论、能量理论、冲击倾向理论、三准则理论和变形系统理论等。

1.1.2 煤岩动力灾害产生机理

煤岩动力灾害现象的机理非常复杂，众多学者对其产生机理与预测预报进行了诸多工作。因本书主要是对煤与瓦斯突出、冲击矿压等煤岩动力灾害过程力电耦合进行研究，因此，下面主要对煤与瓦斯突出、冲击矿压的产生机理进行阐述。

（1）煤与瓦斯突出机理［18～24］

煤与瓦斯突出是一种极其复杂的含瓦斯煤岩动力灾害现象，是矿井生产中重大自然灾害之一。对于煤与瓦斯突出等煤岩动力灾害现象，各国的研究者为认识突出机理付出了艰辛的努力，取得了很大的进步。然而突出的不断发生说明，人类对突出的认识远没有达到尽善尽美的地步。在前苏联学者斯柯钦斯基和霍多特提出典型的突出综合假说之前，已有100余种突出假说，但概括起来主要有以瓦斯为主的假说、以地压为主的假说和化学本质假说等单因素假说。到目前为止，化学本质假说在现场观察和实验室实验两个方面都没有得到支持，已为绝大多数研究者所抛弃。而以瓦斯为主和以地压为主的突出假说只是从一个侧面来说明突出的内在机制，还有许多不能解释的现象。

要在单因素假说阶段与综合假说之间划出一条时间界线是很困难的。苏联学者聂克拉索夫在20世纪50年代提出地压和瓦斯综合作用假说的前后，又有各式各样的单因素假说出现，而且聂克拉索夫的综合假说也不是经典的综合假说。实际上，综合假说应以20 世纪50年代中期苏联学者斯柯钦斯基根据突出煤层的经验和当时的科研成果提出的综合假说为标志。该假说认为，突出是地压、包含在煤体中的瓦斯、煤的物理力学性质、煤的微观结构、宏观结构、煤层构造及煤的自重力等因素综合作用的结果。此后，著名学者霍多特提出了综合假说的能量假说，使突出的综合假说更加完善。前苏联斯柯钦斯基矿业研究所的霍多特等人在实验室中对煤的孔隙结构、吸附性能、渗透性能和力学性质进行了大量的测定研究，并且在压力实验机上做了煤与瓦斯突出的模拟实验。在此基础上，阐述了煤与瓦斯突出的能量来源，从能量的观点出发，用数学方法计算了煤层的变形潜能、围岩的动能、瓦斯的膨胀功和造成突出所需的功，提出了能量假说。能量假说自问世以来对突出研究起了很大的促进作用，其中的大部分观点至今仍有指导意义。然而以霍多特为代表的综合作用假说的主要缺点是：没有说明含瓦斯煤的详细破坏过程及破坏条件，没有考虑时间因素在突出中所起的作用，因而也就无法解释石门的自行揭开和延期突出等现象。虽然在能量假说的前后提出了各式各样的综合假说，如前苏联马凯耶夫安全研究所的包布罗夫的应力分布不均匀假说，英国学者Pooley[19]和Farmer[20]等提出的动力效应假说，但都大同小异，且其共同点是承认突出是瓦斯、地应力和煤的物理力学性质三个因素综合作用的结果，其分歧点是哪种因素起主要作用。和霍多特的综合假说一样，他们都无一例外地忽略了时间因素对突出的影响，没有说明含瓦斯煤体的破坏过程和具体条件。他们也无法解释石门的自行揭开和延时突出等自然突出的现象。中国矿业大学的蒋承林等[21]对石门揭煤时的煤与瓦斯突出进行了实验研究，提出了煤与瓦斯突出机理的球壳失稳假说。张许良[22]研究了利用三维地震探测瓦斯突出危险带的技术，并在淮南潘三矿进行了应用和验证。研究结果表明：高分辨率三维地震可以探测落差大于3 m的小断层和相应幅度的小褶皱，断层和褶皱可以造成构造应力集中、煤体破坏、瓦斯聚集和突出，通过对断层破坏和褶皱变形程度的定量评定可以划分瓦斯突出危险带。应用三维地震技术探测瓦斯突出危险带具有广泛的应用前景。

在煤与瓦斯突出过程及其机理的研究过程中，采取了各种各样的方法，但归纳起来主要有观察统计法、逻辑演绎法、实验法、数值分析法和数值模拟法。早期对突出机理的认识主要是建立在观察统计的基础上，辅之以逻辑演绎，或这些方法的综合。对突出机理的认识程度取决于观察的深度和广度。突出的复杂性和不可重复性使观察统计十分困难，因此，基于个人观察统计基础上的突出机理不免带有某种片面性。随着观察统计数据的积累，实验室实验的开展，人们逐渐弄清了突出的大致轮廓。再加之相关学科（如岩石力学、岩石断裂力学、流变力学）的进展和高新技术（如SEM，计算机技术等）的应用，研究者在前人大量突出实例观察和实验基础上，综合数值模拟技术，加上新的实验室成果，提出了更为完善的突出机理，突出流变机理就是其中之一。

但是，上述假说大多数是以力学为基础从宏观方面建立的，很少考虑含瓦斯煤岩材料的详细变形及破裂过程，这是目前煤与瓦斯突出研究无法取得进一步发展的主要原因。含瓦斯煤岩是一种由空间分布不均匀的煤岩颗粒组成的多孔固体骨架-孔隙中的自由瓦斯气体-孔隙表面的准液态吸附层组成的典型的非均匀多相介质。这种介质的性态和行为不能用任何单一的理论来进行较全面和符合实际的描述。例如，固体力学或岩石力学一般是用纯力学理论来研究煤岩等固体孔隙介质的性态和行为，在考虑孔隙流体的作用时主要采用纯力学作用的有效应力原理；流体力学或渗流力学主要研究流体在管孔或裂隙中的运移规律，而对固体的作用没有加以考虑。

大量的实际现象和实验结果已经表明，含瓦斯煤岩这种特殊三相介质的性态变化和运动过程是力学（岩石力学和岩石断裂力学）、物理学、表面物理化学、电动力学、声学、热力学与传热学和渗流力学等的综合作用，它与单相或两相介质的性态变化和运动过程有本质上的不同。因此，必须采用宏观和微观相结合，将上述学科或理论相交叉的观点及方法，来描述和揭示含瓦斯煤岩这种特殊介质的性态及其动力灾害过程。

从普遍意义的角度来说，自然条件下的煤岩及各种混凝土建筑物都属于含孔隙流体的介质，只是在通常的条件下人们忽略了孔隙流体的作用，而将其简化为单相或两相介质。而且，这种简化是允许的，且偏差不大。但是，这种简化后的物质与实际情况是有偏差的，特别是在孔隙流体的性质比较活泼，压力较高的情况下，简化处理会导致错误的认识结果。因此，研究含瓦斯煤岩动力灾害过程不但对人们认识煤（岩）和瓦斯突出灾害机理、冲击矿压发生机理等有重要作用，而且对认识和揭示地震、滑坡机理、对大型混凝土坝基和高层建筑物的稳定性与寿命评价等具有重要的指导意义。

防治突出灾害的理论基础是突出机理，煤与瓦斯突出的综合假说得到了广泛的承认，该假说认为突出是地应力、瓦斯和煤的物理力学性质综合作用的结果。但是，由于突出的复杂性以及研究方法和研究手段的限制，突出的实验室模拟和现场实际观测都存在着无法克服的困难，我们尚未完全弄清其机理，还无法准确地预测和监测突出，也无法从根本上阐述突出所带来的危害。因此，迫切需要采用新的方法和手段研究煤与瓦斯突出等动力灾害现象。

（2）冲击矿压产生机理［25］

长期以来，冲击矿压作为岩石力学的重大难题之一，一直是国内外学术界和工程界所关注的问题。冲击矿压的发生机理也十分复杂，国内外学者在现场实际调查和实验室研究的基础上对此进行了一系列的研究。其产生机理可以归纳为以下几个理论。

强度理论

强度理论认为：当采场周围产生应力集中且应力达到煤（岩）强度的极限时，煤（岩）体突然破坏，从而形成冲击矿压。如20世纪30年代末被提出的拱顶理论和悬臂梁理论等近代强度理论均以“矿体-围岩”系统为研究对象，认为煤（岩）体的承载能力应是“煤体-围岩”系统的强度，导致煤（岩）体破坏的决定因素是：①应力值大小；②它与岩体强度的比值。70年代Brauner提出了煤（岩）体夹持理论，并导出煤体极限压应力的计算公式。

能量理论

能量理论是从能量转化的角度来对冲击矿压的产生机理进行解释。该理论认为“矿体-围岩”系统在其力学平衡状态遭破坏时释放能量大于消耗能量就发生冲击矿压。各国学者对此进行了不懈的努力，下面对一些主要成果进行介绍。

从能量的来源角度分析，Petukhov认为冲击矿压产生的能量由被破坏的煤（岩）积聚能与和邻接于煤柱或煤（岩）层边缘部分的弹性变形能组成。剩余能量理论于20世纪70年代有美国人提出，该理论认为剩余能量的存在是围岩动力失稳的力学原因。

从能量的判据角度分析，20世纪70年代Brauner提出了冲击矿压的能量判据，该判据考虑了冲击矿压发生时能量释放与时间的相关性，但没有考虑能量释放的空间效应。

尽管如此，对于能量理论的研究，还需要开展进一步的研究工作。

刚度理论

刚度理论最初由Cook等人根据刚性压力机理论得到。该理论认为：矿山结构的刚度大于矿山负载系统的刚度是发生冲击矿压的必要条件。近年来Petukhov则进一步将矿山结构的刚度明确为：达到峰值强度后其应力-应变曲线下降的刚度。

冲击倾向性理论

冲击倾向理论是由波兰和前苏联学者提出。冲击倾向性是指煤岩产生冲击破坏的固有属性，是产生冲击矿压的必要条件。我国学者对此通过大量研究，提出了用弹性能量指数、冲击能量指数以及煤岩动态破坏时间等三项指标来综合判断煤岩冲击倾向性的实验方法。反映冲击倾向性理论的还有顶板冲击倾向性的研究，包括顶板弯曲能指标和长壁开采方式下顶板断裂引起的煤层冲击等方面的研究。

此外，在试验方法、数据处理及综合评判等研究方面也取得了一定的进展。

稳定性理论

20世纪60年代中期NevilleCook首先应用稳定性理论来研究煤岩冲击矿压。其后，Lippmann提出了“初等理论”，章梦涛提出冲击矿压的失稳理论，并在实际中得到了初步的应用。

Bazant等分析了近壁裂纹扩展引起的能量耗散及尺度效应，使对冲击矿压的能量估算成为可能。张晓春等探讨了煤矿巷道附近围岩层裂区的形成和破坏机理，通过理论分析和试验模拟，建立了煤矿片帮型冲击矿压发生的层裂板结构失稳模型。

此外，材料破坏的分叉理论与近年来发展很快的突变理论均在冲击矿压发生机理研究中取得了一定的进展。

总的来说，冲击矿压是裂纹扩展及变形局部化导致的失稳现象，与具有裂纹的各向异性煤岩介质的力学性质和围岩在载荷作用下应力场和应变场的演化和失稳有密切的关系。

1988年，法国Nancy大学的J.L.Mallet推出了地质目标的计算机辅助设计（GOCAD）研究计划，目的是要开发一种新的地质目标计算机辅助设计方法和平台，以适应地质、地球物理和油藏工程的需要。国际勘探地球物理学家协会（SEG）和欧洲勘探地球物理学家协会（EAEG）对这项研究计划非常重视，在1992年末，成立了SEG/EAEG 3D建模委员会，开展了3D SEG/EAEG建模工程（SEM），在1994～1998年间，该委员会分别发表了三次SEM进展报告和修改报告以及两篇阶段报告。目前GOCAD研究计划和SEM工程的许多研究成果已经被国外的许多地球物理公司和石油公司使用，国外大型石油公司的地球科学专家预测，以模型为基础的数据采集、处理和解释已成为石油勘探技术各环节的纽带，对自然资源工业的新的变革将起到至关重要的作用。

SEG/EAEG工程目前努力的目标是帮助设计盐体和上冲断层的三维模型和在这些模型基础上的实际三维勘探的建模。考虑到这个工程的重要性和所涉及的范围，所以常常通过报告来通知全体成员有关它的进展，并结合通知征求资助。关于SEG/EAEG三维建模工程的第一份报告在SEG/EAEG三维建模委员会进展报告（1994）中和1994年第2期的《First Break》（P.57）中已给出。在SEG开罗会议、斯坦福勘查工程（SEP）和科罗里达矿业学院，波现象会议中心展示了他们的第二份报告成果。本节介绍第二份报告的核心内容。

1.6.1盐体模型工作进展

SEG/EAEG盐体模型的第一版本已经建立，并且可通过Internet获得它的GOCAD格式。这个版本主要描述岩体模型的结构组件，在X和Y方向上为9000ft，Z方向上为24000ft，图1.6显示岩体模型主要组成部分的三维可视化结果。

图1.6　盐丘模型（F.Aminzadeh等，1994）

盐体周围的速度是典型的墨西哥湾沉积岩速度，它通过基于k-V0（k值的空间变化）曲线上的致密梯度线和一个地压力面来描述。然而，由于所使用的有限差分软件和可用计算机资源的局限性，用一个常密度来约束盐体模型。这个常密度假设和简单的速度梯度导致在盐体周围的沉积岩中没有地震反射。图1.7表示的是目前正在通过二维有限差分建模来研究的能够产生反射的两种速度模型，第一种技术（“尖脉冲”技术）使每层边界最近处的有限差分网络单元的速度增加了百分之几，第二种技术（“块”技术）交替地先增加然后减小模型中相连结各层的k-V0函数中的速度，因为这两种技术都能改变模型的地震响应，使用价值二维有限差分建模测试就正演合成的地震响应和计算代价来说这两种方法哪一个更合适。使用“尖脉冲”技术，已经获得盐体模型的二维横截面的初步效果。从这些初步结果来看，似乎说明“尖脉冲”技术更合适，当然，还需要更多的测试。

图1.7　盐丘模型速度（F.Aminzadeh等，1994）

同时，一个使用“尖脉冲”技术的三维速度网已经产生并被提供给国家试验室（the National Labs）用来测试三维有限差分软件。

另一个广泛研究的领域是怎样利用美国能源部提供的资源在盐体模型上实现三维有限差分模拟。初步估计这所需要的计算资源远超过目前的预算，由于这个原因，决定使X、Y、Z方向上的模型尺寸减半。这样，总的尺寸减小8个因子。另外，记录时间减少2倍。炮点的数目减少4倍，总的计算量节省64倍。与模型边界有关的速度和结构保持不变。尺寸的减小使计算限制在速度网格范围内，速度网格的确定有一定的考虑，如果太小，计算量增加；如果太大，地震子波的中心频率减小，将导致盐体的具体细节特征不可见。目前使用的速度网格大小为80ft，期望得到大约15Hz的中心频率。另一个正在研究的问题是，使用推荐的二维有限差分模拟算子时，如何选择算子的次数。

除二维有限差分测试外，正在为产生三维零相位差分模拟作准备，许多现象诸如孔隙（光圈）、照明度、分辨率、阴影区域等问题将通过三维零相位模拟来确定。

分析了几个针对盐体模型的三维采集方案。当然关键问题是用减少炮点数来减少预算投资，将接收点的数目减少到与数据存取可接受的水平。所有这些必须实现，同时仍然要保证数据结果的质量。

1.6.2　逆掩断层模型工作报告

地质模型。逆掩断层模型描述了一个不整合于早期延伸裂缝序列之上的复杂上冲地层，这个三维逆冲结构体已经从二维对称平衡剖面上构造而成。它表现为两个相交的上冲断层加上一个附加的横向上逐渐消失的隐蔽冲断层，逆掩断层及底部的一些断层已经被利物浦大学的断层分析组确证。需要指出的是这个模型包括了各种复杂情况：中央一个冲断背斜，外部单斜褶皱和平台区域。逆掩断层（上冲断层）的顶部被风化，并被沉积物表层所覆盖，管道和展开在透镜体上的裂隙在一些层中出现，层的总数为17，模型的尺寸20km×20km×4km，三维表面模型已用GOCAD软件建立并以GOCAD格式存贮。层面模型用边长25m的立方体组成的规则网格表示，每个单元用该单元中心点的层位标记，以生成标记体网，图1.8表示了一个网格化模型的子域，模型中两相交的断层和它上部的通道清晰可见。

图1.8　超覆构造模型（F.Aminzadeh等，1994）

地震速度选取。速度场的定义是将模型用于实际的关键，使用逆掩断层模型的目标在于研究当前不同类型的复杂情况下速度估计问题，例如在地下深400m的浅部重要速度变化和复杂部位横向上的重要速度变化。由于事先是很难判断一个给定的复杂速度场是否合理，因此速度的选择实际上是一个迭代过程，每次循环都包括了由地质工作者定义的三维速度，并由此产生一个模型倾向上二维炮点的全排列及其处理。每层都给定一个与岩性和深度有关的速度，对大多数层位来说，存在范围在2500m·s-1到6000m·s-1的常速度的规律。此外，对表层和另外三层引入垂向和横向速度不均匀性：在表层中，速度横向上在1900m·s-1到3100m·s-1之间变化，相应的长度约为1km，结果造成信号的视周期序列的结构假象（在时间剖面中），在其他三个非均质层中，相应的长度是好几千米。对于这体网格化来说，两均匀同质层之间的界面的几何形态是通过一个由等于离散体网格跨距组成的不规则面来近似的，当这个网格间距约为波长的十分之一时，将导致离散界面的每一步产生绕射。理想情况下，需要对与两层之间的界面相交的那些单元进行计算，得到一个加权速度，实际所采用的方法是对模型的离散用偏移网格坐标轴来完成，用一个算子长度大约是50m的去假频滤波器来给出最终的速度网格。

数值参数的确定。逆掩断层模型的数值参数在二维中通过收敛测试来确定，在三维空间实际上是通过取自于逆掩断层模型的垂直速度剖面构筑而成的水平层状介质来确定的。利用于轴对称介质的比较算法，对结果进行比较，选择中心频率等于15Hz的雷克源，空间域十阶，时间域二阶，时间间隔1ms，空间间隔25m，这组参数得到的结果似乎可满足要求。如果我们考虑对10km×10km×4km大小的子域进行一次炮点模拟（中间放炮）并且记录延迟3秒的地震波曲线（这足够记录炮检距为4km的反射），在一个每秒10亿次浮点操作的超级向量计算机上应用这个参数集有以下性能：对330M内存来说，每个三维炮点需运算70min，这个估计已经在实际三维模拟中证实（见图1.9的地震曲线图）。

图1.9　超覆模型（F.Aminzadeh等，1994）

模拟中的主要阶段。鉴于数值计算的规模，必须把地震建模分成不同的阶段。在美国国家实验室的帮助下，前两个阶段在1994年完成，随着第一阶段结果的成功评估接着进行下阶段的工作。第一阶段的目的是检查模型与客观物体的相关性及研究倾向方向上合适的炮间距。研究的关键问题是三维效果和由于浅部速度变化而产生的失真问题。这阶段相当于在y=11300m处位于模型的复杂区域的倾向上（即x方向）获得一条线的炮点排列，建议记录240个炮点，这些点沿着这个区域面以炮点间距50m排列，波源的中心频率是15Hz（雷克子波），地震波场在一个子域内x和y方向上都是50m间隔的网格面上的每个网格点上以三维方式记录下来，这个子域相应的有一个最大4000m的炮检距。另外，一对深井VSP结果沿着这个剖面记录下来，15个炮点被定位在一个井中用来模拟一个相反的三维噪声监测（walk-away）。在1994年7月底完成综合结果。国家试验室和IFP各自计算一半的炮点，这个阶段的总的存贮需要估计是21GB。这是基于一个采样率为4ms的3s记录来估计的，最终也可能使用一个8ms的采样率。如使用4ms的采样率和单浮点存贮（32位），每道相应为3.1kb；总共将有670万道。第二阶段建议记录一个441个炮点的初始三维测量，这些炮点在一个12km的x和y方向间隔为600m规则网上击发，每个炮点由一个小的规则检波器网记录，记录区域是以炮点位置为中心的8km×8km的正方形，x方向和y方向的间隔是50m（每炮点160×160道）。这完整的叠加将等于49（7×7），总的存贮量估计是35GB（10.3兆道）。数据集的目的首先是为了获得三维叠后量和测试叠前偏移软件（炮点道集偏移），然而，由于没有边缘效应影响区域仅是8km×8km，这个数据集仍然是非常有限的，并且炮间距太大而不能有效测试大多数处理软件。

1.6.3　计算技术和算法工作进展

美国国家试验室的一个主要的研发努力是实现IFP三维反余弦波传播代码的并行版本。第一个通过IFP顺序代码获得来作为参考；第二个由SANDLA实验室使用Paragon产生，两图形很相似。随着时间的增加，差值减小，差的极值处与地震记录的极值相对应。例如，用时间和值表示在图形上的地震记录的最高幅值，最大差值也出现在同样时间并且相对误差小于0.0001。它们可通过在不同平台上的算法表示的差异来解释：一种是IBM的表示方法，另一种是IEEE的表示方法。结论是数据集的生成接近初始模型。

活动断层的定义随不同国家与地区及不同学者而有所不同。目前学界尚无统一标准。断层的活动具有时代性的消长，我国地质学界和工程地震学界普遍认为，活动断层是指晚第四纪以来有活动的断层。但由于各地区的地质环境不同，研究程度不同，各学科的研究目的和研究方法不同，使得国内外学者对活动断层的含义和时限认识也不尽相同（徐锡伟等，2006；景彦君等，2009）。

断裂构造可能会成为CO2泄漏通道，需要对断裂构造的特征进行调查。如存在活动断裂，可能会引起地层断裂、诱发地震的危险，对CO2地下储存库危害较大，因此必须开展断裂及其活动性调查。

（一）断裂调查

1.调查方法

采用大、中、小构造相结合，遥感解译与实地观察相结合的方法，首先确定断层是否存在，然后进一步收集有关资料。当已知或怀疑有断裂时，所需的调查应包括地层和地形分析、大地测量和地球物理调查、槽探、钻孔、沉积物或断层岩的年龄测定、当地的地震调查和任何其他用以弄清运动最近发生时间的适用技术，对在照片上由遥感成像显示的一切线性地形特征等，均应进行足够详细的调查，以解释它们的成因。

断层证据主要有：

1）地貌标志（断层崖、断层三角面、错断的山脊、水系、泉水的带状分布等）；

2）构造标志（线状或面状地质体的突然中断和错开、构造线不连续、岩层产状急变、节理化和劈理化窄带的突然出现、小褶皱剧增以及挤压破碎、擦痕等）；

3）地层标志（地层的缺失或不对称重复）；

4）岩浆活动和矿化作用（岩矿、矿化带或硅化等热液蚀变带沿一条线断续分布）；

5）岩相厚度标志（岩相和厚度的显著差异）。

2.调查内容

1）断层两盘的地层及其产状变化；

2）断层面产状（直接测量、根据断层“V”字形法判定，借助于伴生构造判定）；

3）断层两盘的相对运动方向（主要根据两盘地层的新老关系、牵引褶皱、擦痕、阶步、羽状节理、两侧小褶皱、断层角砾岩等）；

4）断层破碎带的宽度；

5）断层岩类型；

6）断层的组合形式（如正断层的地堑和地垒、阶梯状断层、箕状构造、逆断层的单冲型、背冲型、对冲型、楔冲型、双冲构造）。

（二）断裂活动性鉴定

1.断裂活动性鉴定对象

断裂活动性鉴定的对象是“主要断层”，一般是指：

1）区域地震构造图上有标示的区域性断层；

2）长度大于10km或大于15km的断层；

3）对其活动时代的认识有分歧，并且可能影响到场地地震危险分析结果的断层；

4）晚更新世以来有活动迹象的断层；

5）通过场址区并且与工程场址区安全性评价相关的断层；

6）与破坏性地震特别是M≥6级地震在空间位置上相关的断层；

7）与现代小震密集活动或条带状分布相关的断层；

8）可能延伸到近场区内的活动断层；

9）指向工程场地，并且可能对工程场址区安全性评价有所影响的断层。

2.活动性鉴定内容

1）断层的活动时代。断层活动时代的鉴定是判定该断层是否是发震构造，是否对场址区拟建工程产生重要影响，不能改变路由的管线工程是否采取相应的抗断措施的主要依据。

2）断层的活动性质。对于活动断层而言，其活动性质是划分相关潜在震源区并确定其震级上限的重要依据。潜在震源区范围与边界的确定，与活动断层的性质（包括产状）密切相关。在近场区发震构造评价工作中，应通过野外现场调查或采用成熟技术方法的探测，查明活动断层的活动性质，鉴别出正断层、逆断层、走滑断层、正－走滑断层、走滑－正断层、逆－走滑断层、走滑－逆断层等。

3）断层的运动特性。断层的运动包括“蠕滑”和“黏滑”两种特性。以地震的方式释放的能量往往只占活动断层应变积累的一部分。

4）断层分段性。断层的分段性是确定相应潜在震源区边界及其震级上限的主要依据。断层的分段性研究包括活动性分段和破裂分段两方面的内容。

活动性分段主要包括活动时代与活动性质分段。断层活动时代的差别是断层分段性活动最为显著的标志，在调查中，应当首先加以鉴别，判定“活动的段落”和“不活动的段落”。对于活动的段落，还应视工程的需要和可能性，进一步对其最新活动时代以及活动性质的差别加以细分。

破裂分段是一项难度很大，专业性更强的具有研究性的工作。由于它具有较大的不确定性，只有在工程必需的情况下，可进行专题性研究。

3.活动断裂调查鉴定技术

对目标区内的活动断裂进行详细探测和定期观测，调查其规模、性质、方向、活动强度、特征、地貌地质证据及其活动规律，并初步评价各活动断层的地震危险性。调查过程中应安排槽探、浅井工作，必要时施以地球物理勘探等手段，并采集样品进行地质年代测试。

我国活动断裂调查及研究方法研究较为成熟，调查研究技术手段有地球化学异常、地球物理异常等，并且尝试给出最佳的组合方法。邓起东等（2007）指出小间距钻探和槽探是研究断层新活动的有力手段，可以揭露活动断层最新活动和古地震错动历史的最好技术，并且中国地震局《活动断层探测（DB/T15—2009）》中给出了槽探、钻孔探测的精度适用范围及技术要求，《工程场地地震安全性评价》（GB17741—2005）也介绍了活动断裂调查鉴定技术。

1）进行主要断层活动性鉴定，应以地质地貌学的调查分析方法为主。在进行地质地貌调查与分析时，应注意：

①宏观人手。如断层所在地区的新构造活动背景、断层与第四纪新地层的关系、断层与地貌面的关系、断层的构造地貌特征等。

②微观取证。仅根据宏观现象说明断层的活动性是不足取信的，应选择典型地段和典型部位，通过现场调查，获得断层活动性确切的地质地貌证据。

③精细分析。对于活动断层，应采用断层地貌分析、断层活动性参数确定、古地震探槽、活动性分段、活断层填图、新年代学测定等多种技术方法进行现场调查取证，必要时进行活断层填图，详细鉴定其活动性；

④综合判定。应综合地震活动性、现代构造应力场等不同学科的资料，综合断层活动性的宏观及微观资料，进行断层活动性的综合判定。

2）断层最新活动时代的鉴定，在很大程度上要借助甚至依赖于新年代学测定技术。年代测定方法选择上应因地而异，有所侧重，同时又尽可能采用多种方法进行综合测年。一般来说，对有第四纪地层出露的地区，可采用放射性碳（14C）法、释光法、孢粉分析法；对基岩地区的断层泥的测年可采用释光法、电子自旋共振（ESR）法、钾－氩（K-Ar）法和电镜（SEM）扫描法等。

3）在覆盖区，已有资料不能确定已知主要断层的活动时代时，应选用地球物理、地球化学、地质钻探和测年等手段进行勘查。隐伏断层的活动性鉴定一般应遵循以下步骤：

①进行隐伏断层位置的初步探测。根据航、卫片判读和已有的地质、地貌、化探、物探、钻探资料进行综合分析，初步推测断层的位置、延伸和展布形态，然后选择适宜的探测手段，布置探测路线。

②进行隐伏断层的综合探测。在初步推测出断层的大体位置后，进一步按照先粗后细的原则，选择合适的物探或化探方法，初步确定断层位置。再进行浅层物探，如浅层地震勘探、地质雷达等，以查明隐伏断层的确切位置和断距。

③根据具体情况进行钻探和槽探，进一步帮助确定断距、断面、断错地层及上覆地层，并采集合适的样品，综合分析其活动性。

---