那么有这五个技术环节，我们研究地震对象本身还可以分成三大过程，地震孕育的过程，地震破裂的过程，震后调整的过程，孕育特别临近破裂的时候，可能会有些非现象和现象，因此在一些前兆上可能有一些反应。那么这是一个几年、几十年孕育的过程，甚至百年、千年，地震破裂的话可能零点几秒、几秒，大地震在一二百秒也就完了，所以是一个很顺序的过程，那么震后调整的话在几个月或者几年之内的话有时候还有比较大的调整大，相当比这时候可能快，竖一个像级的差别，所以也是单独需求这个提出研究的一个阶段。 那我们就看看展望一下看看这三个方面我们都做的怎么样，孕震方面，那么你们特别有非地震专业的学生，那么你们可能到网上查过，为什么会有汶川地震发生，这是GPS测量出来的一个运动的方向和大小的这样一个图，那么到网上一查，这是我抄来的，印度板块与欧亚板块发生碰撞，那么碰撞的时候的话呢造成青藏高原快速隆升和物质向东缓慢流动，那么然后的话呢碰上了四川盆地的干性的板块的阻挡，造成各大因素其中，这就发生地震，那么这种说法对不对呢？是对的，但是够不够呢？这个应该说是不够的，因为你随便问就可以提出几个问题。 像在四川盆地这一次，我们看就没有大地震，在高原一侧有一些比较大的地震，但最大的地震反而在这里，为什么？地震发生在10-20公里的深度为主，为什么？为什么不是更深或者更浅，那么在南边我们主要是挤压逆掩断层，在东边的话除了逆掩断层之外还有水平错综的成份，为什么？随便提一些问题的话，你刚才那样简单的定性的结论的话就没法回答了，那么我们就根据刚才原则，根据我们对结构的了解，对岩性的了解，可以建立一些三维的模型来进行计算，计算的时候用的边界条件，可以用GPS测量的差值来给的边界条件，那么在这里我们要先说结果事先要说一个，因为有些人的话觉得你们搞模拟的话你们想算出什么，还不能辨别这辨别那就算出一个东西来，你们算的东西信得过吗？ 那么我们说这里的话我们是用一个弹性模型算的，那么这图中的红色表示应力高的，蓝色表示应力低的，那么这两张图的话表示垂直于断层走向的压应力，这两张图是表示右旋水平的剪应力，这张图是地震活动性的，我们看的出来地震活动性高原一侧的话在这里的话应力并不特别高，那么在剪应力在这种地方就有，这个地方根本没地震的，所以的话我们用这种弹性模型的话算出来的高应力其实跟地震活动性没有任何关系，也就是弹性模型没法解释地震活动性特征的，但是反过来，我们了解地下的结构，我们知道像这种地方，比较红的地方，就是地震波动比较慢的地方，就刚才这里是汶川，那么地震波速度慢的地方往往是温度比较高的地方，温度比较高的地方它的这个柔性也就比较高，岩制系数会比较低，所以同样的跟刚才一样的，别的都一样，只是说你认为理想弹性是不对的，我们明明知道高温高压下岩石是柔性的，你不能用这个理想弹性去算了，那么我们用马克斯岩性体去算，甚至跟刚才的图的话完全不一样了，那么这边的话应力很低，这是对应没有地震的，这边应力比较高，对应的是有地震的，这不叫应力，叫应力增长率，因为我们不知道应力的绝对值，我们它作用以后的应力增长率，应力增长较高地段是在上地壳的底部，大概10-20公里的地方，这是红的，所以对应的最大岩石活动性，那么挤压的话我们看出来在西南段是比较大的，所以这边是以断层为主，在东北段的话呢这个剪应力越来越增加，所以东北的话有一种减轻的分量。 那么在算这样的问题的时候，我们输出的是什么，是基于观测的结构模型，基于观测核试验事实验的物性，基于地表GPS观测的边界条件，不知道什么呢？深部的边界条件，我们只能做一些简化的假定，初始应力条件，不知道，我们之梦从0应力出发计算，边界位移的历史我们也不知道，所以假定从0起后是个恒定位移的数据，我们得到的是什么呢？稳态或者叫准稳态，应该说应力的增长率，而且发现这个应力增长率的高低与地震活动性是有相关性的，我们得不到什么呢？地震应力的绝对值，所以究竟哪里最危险我们不知道。 我们肯的出来整个龙门山上地壳的底部都是高应力增长区，现在我们只断了这三百公里，这有一百公里的话没有断，那么就是在我们做芦山地震之前的话我们也看的出来，这是个高应力增长率，那我们已经知道了，但是我们不知道它的应力演变历史，所以应该加强应力的绝对测量，就是应力状况，和应力变化叫做相对测量。 那么我们看从原来定义的高应力增长区，我们也可以算出一个地震以后会引起应力的变化，这样从而把危险区的话变的更危险的断层、地区勾画出来了，那么这里也是变的更危险的一个地区。所以的话芦山地震可以说呢我们一点也不觉得它意外，因为本来呢三百公里一撮，这一百公里没撮，而且这一百公里是，本来也就是一个应力增长率比较高的地区，何况汶川地震又加剧它的危险性，但到底三年震、五年震、二十震、一百年震我

那么我们做一些工作的话就是去定量计算一种所谓库仑应力，什么叫库仑应力呢？就是两个断层我们会给剪切应力，剪应力的话呢就会触动，但我要给正应力把它压的实一点呢？你就同样的剪应力可能就错不动了，所以这种压应力越大，正应力的话越大那么它就越安全，剪应力越大的话呢就会越危险，另外还有一个孔隙水压力的问题，那么如果固定存在全部应力，那么你压应力大的时候就困难，但往你孔隙水里充进了水，水承担了一部分这个应力，那么固体颗粒之间的压力，正应力变小的话它也变的更加容易错动，那么库仑应力是三种因素的综合，库仑应力变大了，意思就是更危险了，库仑应力变小了就更安全了。 那么我们呢做的话呢就是说我们把解析结，数字结，数字结的二维结和三维结，都做了比较完整的对照，而且我们的解析方法的话，数学方法的话呢要跟解析方法对应的上才时间我们数字计算是可靠的。那么在解析结里头我们看的出来，仅仅水库的载荷的话呢是让这个断层底下的断层变的更加困难了，更加安全了，所以单就蓄水来讲的话这个紫坪铺上方越压这个弹性载荷来讲它变的更安全，但是的话水可以渗进来，水渗进来的话总是让它变的更危险，压的地方是有些断层变的危险，有些断层变的安全。但不断什么样的断层，只要水渗进来的话总是变的越来越危险，那么在二维和三维我们也做了一些对比，因为我们是一个三维的问题，二维的时候你等于认为整个水库是一条线，那么就可以夸大危险的程度。所以二维的话大概夸大了三倍左右的结果，结果我们一般需要三维的计算的，那么三维的话因为我们进行百万单元的计算，所以我们可以算，那么算出来的话就是说加荷、载荷顺便会影响它的危险性，那么加荷反而变的更安全，加荷的时候反而往下更安全，卸荷的时候的话呢水位往下走马上变的更危险，但还另一个因素，水涨性因素单调的总是越来越有危险，所以这个线和这两个合起来的话有可能而产生这个，到接近1kpa应力，实际上变的危险，但这1个kpa是一个很小的应力，所以能不能触发的话就很难说，所以我们做了西洪江的，因为西洪江水库的话是一个公认的水库地震，但西洪江我们做出来的话也就是3个kpa，所以这个问题还不好说，3个kpa大概就像10个厘米高的水，加了10个厘米高的水或者加了几个厘米高的水触发的地震，这个东西，不是说加了这种水，加了相当压力，在断层那个地方的压力相当加了这么多水的压力或者减少这么多水的压力，它就地震了，这个好像自个儿也还不太敢信，但有一点是这样，我们算出来的像这个地方，西洪江应该是下沉18个毫米，实测的话是15毫米，所以算的东西跟观测是对得上的，那么我们算加了载荷以后，它的弹性应力增加的话呢大概是7×10 11次方焦耳，那个六级地震光地震波释放衡量的话大概就是10的14次焦耳，所以这东西连它的1%都不到，那么到汶川地震的话恐怕连百万分之都不到的。所以你说是人造这个地震的话是绝对不可能，还是大自然的构造应力的，构造应变能的积累和刚才讲的那个道理，但是不是有触发作用呢？对于汶川地震还很难说，因为它临近这个边界的临界值了，我们还需要做更多的研究，当然我们还算一些阿斯旺水库等等的，那么希望对有没有这样一个临界值，是多大我希望再做一些工作。 那么再一个是地震破裂的过程，那么像这是美国的工作，那么他们就是利用断层的一些本合关系，速度和摩擦应力的关系去算，那么就是没有地震的时候我可能是一年划一条线，到地震的时候可能就是一秒钟画一条线，那么看的出来，这一秒钟开始在这里，这里是地震发生，2秒钟、3秒钟、4秒钟、5秒钟，纵轴的话是位错量，就是这个位错在这里发生，很快的向南传播，那么然后变成这样一个，然后在若干年之内的话都没有变化，这些密的线，一年走一点，一年走一点，那么走到这里，走到这里以后又发生一次地震，白的地方就是地震的地方，就是他们可以模拟出这个断层上每二十多年发生一次6级地震来，那么大致可以模拟一下，可以模拟出实际上观测到的绝大多数地震都是在两端发生，向另一侧传播，有时候在南端，有时候在北端开始，但是的话完全实际模拟，真实的地震序列那么也还有一段距离，所以的话人们的话认识也还在不断做，我们国家也有一些人在做类似的工作。 那么还有震后的情况，震后的情况同样也讨论一个大地震造成的库仑应力变化是怎样的，那么这里也有一些解析的方法，我们要研究比较真实复杂的情况就要用这个数值的方法，参数的话就节理单元方法，那么怎么处理断层的话过去有两种方法，我们的学生都还做一些改进，我们有两个学生和博士后还提出了自己一些新方法，像过去传统的方法，那么我们也做一些这个计算和加深了解，那个我们就不多说了，那么还有方法是说断层的时候相当于它的地震的抗减模量减小了，发生位错，但是仅仅这么说的话呢，它的抗减模减少同时，抗压模也减少了，因此在压扭的时候这里会产生一些变化，而实际上呢断层的时候不会产生的，