发布时间：
安昌强，宋仲和，陈国英，陈立华，庄真，傅竹武，吕梓龄，胡家富.1993.中国西北地区剪切波三维速度结构.地球物理学报，36（3）：317～325
陈国英，宋仲和，安昌强，陈立华，庄真，傅竹武，吕梓龄，胡家富.1991.华北地区三维地壳上地幔结构.地球物理学报，34（2）：172～181
陈立华，宋仲和，安昌强，陈国英，庄真，傅竹武，胡家富，吕梓龄.1992.中国南北带地壳上地幔三维面波速度结构和各向异性.地球物理学报，35（5）：574～583
董颐珍.1997.秦巴地区QB-1测线地震测深成果介绍.见：袁学诚主编，阿尔泰-台湾地学断面论文集.武汉：中国地质大学出版社，165～175
高锐，吴功建.1996.地球物理综合解释模型与现今地球动力学过程。《青藏高原岩石圈结构构造和形成演化》.北京：地质出版社，138～148
阚荣举，韩源.1992.云南遮放至马龙地学断面说明书.北京：地震出版社
李立.1996.中国大陆地壳上地幔电性特征.地球物理学报，36（增刊）：136～139
李秋生，卢德源，高锐，李敬卫，范井义，熊贤明，张之英，刘文，李英康，闫全人，李德兴.2000.横跨西昆仑-塔里木接触带的爆炸地震探测.中国科学（D 辑），30（增刊）：16～21
林中洋，蔡文伯，陈学波，王椿镛.1992.青海门源至福建宁德地学断面说明书.北京：地震出版社
卢德源，李秋生，高锐，李英康，李德兴，刘文，张之英.2000.横跨天山的人工爆炸地震剖面.科学通报，45（9）：982～987
卢造勋，夏怀宽.1992.内蒙古东乌珠穆沁旗至辽宁东沟地学断面（说明书）.北京：地震出版社
马杏垣，刘昌栓，刘国栋.1991.江苏响水至内蒙古满都拉地学断面说明书.北京：地质出版社
潘裕生，孔祥儒主编.1998.青藏高原岩石圈结构演化和动力学.广州：广东科技出版社，1～20
彭聪，高锐.2000.中国大陆及邻近海域岩石圈/软流圈结构横向变化研究.北京：地震出版社
宋仲和，安昌强，陈国英，陈立华，庄真，傅竹武，吕梓龄，胡家富.1991.中国西部三维速度结构及其各向异性.地球物理学报，34（6）：694～707
宋仲和，陈国英，安昌强，陈立华，庄真，傅竹武，吕梓龄，胡家富.1992.中国东部及其相邻海域S波三维速度结构.地球物理学报，35（3）：316～330
孙武城，马宝林，宋松岩，胡鸿翔.1992.湖北随州至内蒙古喀喇沁旗地学断面说明书.北京：地震出版社
孙武城，徐杰，扬主恩，张先康.1992.上海奉贤至内蒙古阿拉善左旗地学断面说明书.北京：地震出版社
汪集旸.1997.深部地热岩石圈研究中的意义及应注意的几个问题.见：张炳喜，洪大卫，吴宣志编，岩石圈研究的现代方法.北京：原子能出版社，113～126
吴功建，高锐，余钦范，程庆云，孟令顺，董学斌，崔作舟，尹周勋，沈显杰，周姚秀.1991.青藏高原“亚东—格尔木地学断面”综合地球物理调查与研究.地球物理学报，34（5）：552～562
熊绍柏，赖明惠，刘若宾，于桂生.1993.屯溪—温州地带的岩石圈结构与速度分布.见：李继亮主编，东南大陆岩石圈结构与地质演化.北京：冶金工业出版社，250～256
胥颐，刘福田，刘建华，孙若昧.2000.天山地震带的地壳结构与强震构造环境.地球物理学报，43（2）：184～193
徐新忠，齐雄飞，扬长来，党淑娟.1993.青海花石峡—甘肃阿克塞人工爆破地震测深成果.中国地球物理学会年刊.北京：地震出版社，80p
徐新忠，王有学，蒋亚明，党淑娟，聂定平.1997.新疆可可托海-甘肃阿克塞人工爆破地震测深剖面的地壳速度结构研究及其大地构造单元的划分.见袁学诚主编，阿尔泰-台湾地学断面论文集.武汉：中国地质大学出版社，1～13
尹周勋，赖明惠，熊绍柏，刘宏兵，滕吉文，孔祥儒.1999.华南连县—博罗—港口地带地壳结构及速度分布的爆炸地震探测结果.地球物理学报，42（3）：383～392
袁学诚主编.1996.中国地球物理图集，国际岩石圈委员会201 号出版物，北京：地质出版社
张贻侠，孙运生，张兴洲，杨宝俊.1998.中国满洲里-绥芬河地学断面.北京：地质出版社
章惠芳.1988.广西柳州爆破观测和桂东地区地壳结构.中国大陆深部构造的研究与进展.北京：地质出版社
章惠芳.1988.中国东南沿海地区泉州爆破深地震测深资料的再解释.中国地震，4（2）：96～102
赵俊猛，李植纯，马宗晋.2003.天山分段性的地球物学分析.北京：地学前缘，第10特刊，125～131
郑晔，滕吉文.1989.随县—马鞍山地带地壳与上地幔结构及郯庐构造带南段的某些特征.地球物理学报，32（6）：48～659
朱介寿.1997.中国及其邻区地球三维结构初始模型的建立.地球物理学报，40（5）：623～638
庄真，傅竹武，吕梓龄，胡家富，宋仲和，陈国英，安昌强，陈立华.1992.青藏高原及邻近地区地壳与上地幔剪切波三维速度结构.地球物理学报，35（6）：694～709
可到它网站线上头稿，也看下投稿要求等等，图片放在WORD里就行，一般要求清晰
苏正1,2，曹运诚1，吴能友1,22，Lawrence ，陈多福1,2
苏正，（1980—），博士，助理研究员，主要从事天然气水合物及盆地流体活动的数值模拟研究，E-mail:。
注：本文曾发表于《地球物理学报》，2009,12:3124-3131，本次出版有修改。
1.中国科学院边缘海地质重点实验室/广州地球化学研究所，广州510640
2.中国科学院广州天然气水合物研究中心/可再生能源与天然气水合物重点实验室/广州能源研究所，广州510640
of Earth&Atmospheric Sciences,Cornell University,Ithaca,New York 14853-1504,USA
摘要：海洋环境中天然气水合物层是理想的毛细管封闭层，游离气被抑制在水合物层下，游离气层的气体压力随气体聚集和气层厚度的增加而升高，当气压超过封闭层的毛细管力时，游离气会克服毛细管进入压力、刺入上覆封闭层孔隙空间，毛细管封闭作用随之消失，从而形成水合物下伏游离气向海底的渗漏。通过对该过程进行的数值模拟计算表明：渗漏气体是以活塞式驱动上覆沉积层中的孔隙水向海底排出，水合物稳定带内流体渗漏速度随水流柱高度的减小而增加，当水流阻抗大于相应沉积层段的静岩压力时，沉积层将转变为流沙，流沙沉积被海流移除后便在海底留下凹陷麻坑。麻坑形成后流体运移通道演化为气体通道，气体快速排放。麻坑深度主要取决于游离气层的厚度和水合物封闭层（底界）的深度，而与沉积层的渗透率无关。麻坑深度一定程度上指示了渗漏前水合物层下伏游离气层的资源量。对布莱克海台海底麻坑的深度数值模拟计算表明，形成4 m深的海底麻坑需要至少22 m厚的游离气层。
关键词：天然气水合物；毛细管封闭；游离气渗漏；麻坑；布莱克海台
Numerical Computation and Case Analysis of the Venting Process of Free Gas Beneath Hydrate Layer
Su Zheng1,2,Cao Yuncheng1,2,Wu Nengyou1,2,Lawrence Duofu1,2
Key Laboratory of Marginal Sea Geology/Guangzhou Institute of Geochemistry,Chinese Academy of Sciences,Guangzhou 510640,China
Guangzhou Centerfor Gas Hydrate Research/CAS Key Laboratory of Renewable Energy and Gas Hydrate/Guangzhou Institute of Energy Conversion,Chinese Academy of Sciences,Guangzhou 510640,China
of Earth&Atmospheric Sciences,Cornell University,Ithaca,New York 14853-1504,USA
Abstract:A hydrate layer is an ideal capillary seal,beneath which free gas is overpressure increases as gas accumulates and gas column seals have the property that they fail completely when gas pressure reaches the point that they are invaded by gas,and thereafter they offer little resistance to gas the seepage is triggered,the venting gas will push the overlying water upward at increasingly higher velocities as the gas “piston”approaches the model shows that as the water velocity increases,the near surface sediments will become quick at a depth that the resistance of water flow exceeds the hydrostatic pressure of the sediment hosting the water quick sediments can then be removed by bottom ocean currents,leaving a hollow pockmark on the gas pathway isformed below the pockmarks and the reservoir gas drains pockmark depth is afunction of thickness of free gas column beneath the hydrate and depth of the hydrate seal (bottom of hydrate layer).Interestingly,pockmark depth does not depend on sediment depth implies the resource amount offree gas beneath hydrate model shows that a 22-m-thick free gas layer at least is needed toform a 4-m-deep pockmark on the rise of Blake ridge.
Key words:gas hydrate;capillary seal;gas seepage;pockrnarks; Blake ridge
0 引言
在海洋环境水合物稳定带内孔隙水溶解甲烷浓度超过甲烷水合物形成的溶解度时，溶解甲烷会结晶形成水合物，随着水合物含量的增加，形成水合物层圈闭，并在其之下发育游离气层[1-4]。在特定的条件水合物层之下的游离气沿通道向上渗漏进入海底，并在海底形成麻坑、自生碳酸盐岩、生物群落、气泡羽状体，如俄勒冈外海水合物脊[5]、布莱克海台等[6]、北刚果陆坡[7-8]、挪威外海[9]以及中国南海[10]。虽然水合物层下伏游离气向上渗漏活动在水合物发育区比较普遍，但是水合物层下伏游离气向上渗漏的机制和泄漏过程中的流体动力学特征，及流体渗漏对海底沉积地层的破坏（形成麻坑）过程并不清楚。
水合物层下伏游离气受到水合物层毛细管作用的封闭，随气体聚集和气层厚度增长，水合物下伏游离气的压力持续增加，当气体超压克服毛细管封闭作用后气体渗漏被激发，超压气体推动孔隙水向上排出，在海底形成麻坑，麻坑深度反映了流体的破坏强度和游离气层的超压幅度。因此，本文将应用水合物层毛细管封闭机理和沉积孔隙流体渗漏动力学，研究水合物稳定带之下游离气如何向上突破的动力学过程，建立游离气层压力状态与麻坑深度之间的数值模型，通过海底麻坑特征揭示水合物系统游离气层的演化规律。
1 毛细管封闭及游离气渗漏机理
海底沉积层中存在2种毛细管力封闭作用。第一类毛细管力封闭作用是存在于小型的气藏顶部的毛细管封闭作用，属于低渗透率的气体捕集封闭。封闭层的孔隙度和渗透率较低，而水更倾向存在于较小的孔隙空间，因此封闭层的孔隙空间完全被水占有，而封闭层之下含气层的孔隙度和渗透率相对较高[11]。碎屑沉积物孔隙介质一般为水润湿相，气液界面处的毛细管力阻止天然气进一步向上运移，使气体处于孔隙较大的沉积层段，但当气体压力超过相应孔隙的气体的毛细管进入压力时，超压气体将刺入封闭层的小孔隙，气藏开始排气，并在上覆沉积层中产生气体的渗漏通道。侵入毛细管压力由拉普拉斯方程给出[12]：
南海天然气水合物富集规律与开采基础研究专集
其中：γ为界面张力，取值 N/m[13],rf和rc分别代表小孔隙和大孔隙的有效孔隙半径。
第二类毛细管封闭作用存在于气－液二相共存的沉积孔隙中，气液二相均可流动[14-15]。由于整个沉积体是由沉积颗粒构成的孔隙介质，孔隙水优先占据并被吸附在孔隙的喉道位置，具有小孔径的孔喉部位产生的毛细管力抑制了孔隙腔中气体的流动。此类毛细管封闭条件是孔隙内2种流体共存，且二者均可流动。在渗漏活动初期这种情况出现在气流柱顶部和气柱周围的气－水混合的部位，沉积层中毛细管封闭线的位置随气柱的发育而变迁，这种毛细管封闭作用约束了气流柱的形状和发育，并使气流柱有一个相对平坦的顶部；同时也会形成一个相对稳定的通道直径，这意味着渗漏气柱顶部的气－水界面相对平坦，在理想均质介质中渗漏气体以“活塞”式向上推进。但是当渗漏气柱遇到渗透率在横向上不均匀或不连续（如断层）的沉积介质时会出现分支或扭曲的气体通道。
海洋环境扩散型水合物稳定带与下伏游离气之间属于第一类毛细管力的封闭，在水合物稳定带底部水合物含量最高[3,16]，水合物的形成降低了孔隙介质的有效孔隙度和渗透率，使水合物层的孔隙度低于下伏游离气层的孔隙度，水合物层的有效孔隙半径小于游离气层的有效孔隙半径。亲水性的水合物沉积层内除水合物外的其余孔隙空间被水占据，而下伏沉积体的孔隙空间完全被气体充填，水合物层与游离气层之间就存在一个上覆孔隙水与下伏游离气的界面。因此在水合物层与游离气层界面（大孔隙与小孔隙之间）上产生毛细管力，其方向指向孔隙半径较大的含气层，阻止下伏气体进入上覆含水层（水合物层），抑制气体向上运移。但是当下伏游离气层中的气体压力超过上覆水合物封闭层的毛细管力时，超压气体将刺入水合物封闭层，使水合物层的毛细管封闭作用完全失效或仅剩很小的封闭作用，气体泄漏开始。超压的气体渗漏进入水合物稳定带后，随着气柱的增长气体逐渐侵占原有孔隙水所占的孔隙空间，驱使孔隙水向上排出，并最终泄漏进入海底。水合物稳定带内气柱的增长过程受第二类毛细管封闭作用的控制，使气流柱以“活塞”式增长，而没有出现气流弯曲和分支，这与地球物理资料显示的近于垂直的流体渗漏通道（气囱）特征一致[8-9,17-19]。
图1给出了海洋水合物层下伏游离气渗漏过程。游离气在水合物层底界之下聚集，气层厚度和气体超压逐渐增加（A），当气体压力超过水合物封闭层的毛细管力时，高压气体会在封闭薄弱点或气层最顶端刺穿封闭，使水合物毛细管封闭失效（B）。气流柱在高压作用下向上推进，并驱使上覆沉积孔隙水向外排出。气流柱高度（hg）逐渐增长，而水流柱高度（hw）相应缩短（B到C过程）。如果气压驱动力保持相对恒定，由于岩层对水的黏滞力（或水流阻抗）远大于其对气的黏滞力（或气流阻抗），随水流柱高度hw减小，流体渗漏速度将越来越快，在单位长度水流柱上的压降（等于岩层对水流的黏滞力）随流体速度的增长而增加。在气流接近海底时流体速度明显增强，浅层水流阻抗（即水流对地层的作用力）超过相应沉积体的静岩压力，浅层含水沉积将被流沙化，当流沙化的沉积物被海底底流搬运后，便在海底形成“新鲜的”麻坑，此时麻坑下形成单一的气体运移通道（D）。由于气体黏度远小于水的黏度（约为1/60），气体排放异常迅速，游离气藏中气体会很快排干，流体渗漏通道中的气流逐渐退化（E），孔隙流体压力回归静水压力，孔隙水重新占据水合物封闭层和流体渗漏通道的孔隙空间，在气量通量减小体系温度降低的过程中伴随者水合物的生成（此文中不做详细论述），并因此减小了流体流动速度，少量气体仍可滞留在流体渗漏通道内，在地震记录上显示为气烟囱，水合物层底部的毛细管封闭作用恢复，水合物层之下游离气的聚集过程再次启动（F）。
图1 水合物下伏游离气渗漏概念模型示意图[11]
Z为海底以下深度，h为水合物稳定带厚度（或水合物封闭层深度）。黑色带表示毛细管封闭层，浅灰色表示气体所占据孔隙沉积层。A.气体被封闭在水合物层之下；B.气体刺穿封闭层开始泄漏C.气柱高度增加，推动水流向外排出，水流柱高度相应缩短，流体运移速度不断增加；D.含水流沉积中孔隙压力超过静岩压力，在海底出现麻坑，形成单一的气流通道；E.游离气藏中的天然气被逐渐排空，孔隙超压消失，流体通道中的气流柱逐渐退化；F.气流柱完全消失，在海底留下气烟囱，并有水合物生成，水合物封闭作用恢复，并开始新的气体聚集
2 游离气渗漏过程的数学模型
气体渗漏过程中（图1）气柱和水柱都是在游离气超压的驱动下流动，流体运移的总驱动力等于气体超压（ρw－ρg)gd。气流柱不断增大，并且以同一速度推动渗漏通道内的上覆孔隙水向上流动。假定水合物稳定带为一种均质孔隙介质，渗漏通道内流体（水和气）的渗漏速率相同，孔隙介质内流体渗漏模型可用达西定律描述为
南海天然气水合物富集规律与开采基础研究专集
其中：Δp为流体运移总推动力，是施加在气流柱和水流柱上的压降之和（Δpg＋Δpw），或者是气流阻抗与水流阻抗之和，等于气层底部的超压（ρw－ρg)gd；ρ为流体密度；d为游离气层的厚度；μ为流体黏度；V为流体速度；k为沉积体的渗透率；krg和krw分别为沉积体孔隙气和水的相对渗透率；hg和hw分别为气流柱和水流柱的高度。
假定气流柱中气的饱和度和水流柱中水的饱和度均为1，气和水的相对渗透率为1。由方程（1），流体（气体和水）的运移速度表示为
南海天然气水合物富集规律与开采基础研究专集
在方程（2）中，若 可知流体运移速度随气流柱高度（hg=h-hw）的增长而增加。对方程（2）进行积分得到气柱增长方程：
南海天然气水合物富集规律与开采基础研究专集
利用方程（3）既可以计算渗漏气流柱增长到某一高度所需要的时间，也可以计算某一时间点水合物稳定带内气流柱的高度。
由方程（1）和方程（2）可知，孔隙介质中单位长度流体柱所受阻抗随气流柱高度的增加（或水流柱高度的减小）而增加，也就是说沉积物格架所受流体的反作用力（流体阻抗）逐渐增加，当流体阻抗超过相应沉积体的静岩压力时，相应沉积层将被流体化而成为流沙[20]，渗漏流体速度须满足 。流沙沉积被海流移除后在海底形成麻坑，被流沙化沉积体的底界确定了麻坑深度。用 替换方程（2）中流体速度V，麻坑深度hpm替换水流柱高度hw，即可得到麻坑深度方程：
南海天然气水合物富集规律与开采基础研究专集
方程（4）中，若μw≌60μg、krw≌krg≌1（假定水流柱中水的饱和度和气流柱中气的饱和度近似为1），方程（4）可简化为
南海天然气水合物富集规律与开采基础研究专集
在一定的温压条件下流体密度和黏度为常数[12]。因此，方程（5）中麻坑深度可近似为水合物下伏的游离气层厚度（d）和水合物封闭层深度（h）的函数，与沉积体的渗透率无关。模型计算中所有参数取国际标准单位。
3 模型应用及讨论
美国卡罗莱纳外海的布莱克海台区是典型的水合物发育区，既有完美的BSR显示，又有游离气的渗漏活动及在海底形成的麻坑[6,21-22]。大洋钻探计划（ocean drilling program)1 64航次对布莱克海台进行了钻探取心研究，其中997站位钻至海底之下750 m，穿过了BSR（海底之下450 m），其中180～462 m 层段含水合物，水合物平均饱和度为6％，位于水合物稳定带底部（462 m）的水合物体积分数最高为24%[4]。996站位于布莱克底辟链的最南端，处于997站位西北98 km，最大钻孔深度为63 m，刚好位于麻坑之中，地震剖面显示该区BSR深度为440 m，深部底辟作用使上覆地层变形、形成小型断层，成为有利的流体渗漏通道，在海底发育有深4 m、直径50 m的麻坑，并且正在发生气体渗漏（图2），钻探获得的水合物体积分数高，最高达沉积孔隙的99%[6,21-23]。
驱动流体运移的气体超压取决于游离气层的厚度。如果下伏游离气层厚度达100 m（图1），其总的流体驱动力（等于气体超压）可达到 MPa；如果游离气层厚度为22 m，流体超压驱动力为 MPa（图3最左端A点）。渗漏开始时水流柱高度分数（等于hw/h）为1，总水流阻抗等于气体总超压，整个气流柱高度增加而降低。但是由于水流速度增加，施加在单位长度水流柱上的驱动力和相应的黏滞力增加，水流阻抗逐渐趋近海底相应深度沉积层静岩压力，且在水流柱高度分别小于40 m（对于游离气层厚度为100 m）和4 m（对于游离气层厚度为22 m）时水流阻抗超过沉积介质的质量（图3D点）。该位置以上的沉积物被流沙化[20]，转变成颗粒悬浮的液状混合体，这种流沙化沉积被海流搬运后在海底形成麻坑。利用方程（3）可以计算游离气从水合物稳定带底部渗漏到达海底所需的时间。假定渗漏率为10-12m2时， 100 m厚的游离气层泄漏到海底的时间大约为5 a。
图2 布莱克海台地震反射强度剖面揭示的BSR、底辟构造、海底麻坑及与ODP977站位揭示的BSR深度比较
a.地震反射强度显示布莱克海台水合物发育、气体聚集以及底辟构造顶端的流体渗漏[22];b.为ODP997站位BSR揭示的水合物封闭层深度[21]
图3 渗漏通道中的流体阻抗和含水沉积层的静岩压力曲线交点指示麻坑深度
水合物稳定带中气流柱高度增加（顶部坐标向右），水流柱高度减小（底部坐标向右），水流阻抗和静岩压力随之减小，水流阻抗大于静岩压力时发生流沙破坏，曲线交点位置指示麻坑深度（D点）。布莱克海台100 m的游离气层发生渗漏时在海底可形成40 m深的麻坑，而22 m厚的气层泄漏时可形成4 m深的海底麻坑（最右边灰色阴影）
方程（2）中流体渗漏速率与渗透率成正比，但方程（4）中麻坑深度不依赖于沉积体渗透率，只是水与气体相对渗透率比的函数，而相对渗透率决定于孔隙流体的饱和度[12]，因此沉积体渗透率控制流体渗漏速率，但不控制麻坑形成。实际上，渗透率越大，气体渗漏越快，麻坑形成越快；气体超压在水流柱和气流柱之间的分配不依赖于渗透率，而是决定于气体的超压幅度，以及流体黏度和气流柱高度（或水流柱高度）。
利用方程（5）可以简单计算海底麻坑深度，同时在已知水合物底界（封闭层）深度和麻坑深度，也可以通过方程（5）计算游离气层的厚度。图4显示麻坑深度与游离气层厚度和封闭层深度的关系。在给定封闭层深度，麻坑深度随游离气层厚度的增加而增大，相反较深的沉积层厚度削弱了渗漏流体对麻坑的挖掘作用，水合物封闭层越浅，形成一定深度的麻坑所需的游离气层厚度越小。
图4 水合物封闭层深度和麻坑深度与游离气层厚度的关系
麻坑深度主要决定于游离气层厚度和水合物封闭层埋深，与游离气层厚度呈正比，与水合物层埋深呈反比。如果水合物封闭层深700m，形成4m深的麻坑需要27m的游离气层，如果水合物封闭深度为440 m，则需要22 m的游离气层，如果水合物封闭层深100m，仅需要1l m厚的游离气层
地球物理显示布莱克海台ODP996站位周围的BSR深度为440 m，而在ODP996站位正下方游离气藏气体沿底辟构造上升至大约220 m（图2）处，在沿小断层渗漏至海底，由方程（5）可知麻坑深度与渗透率无关，取决于游离气藏的埋深和游离气层的厚度。对于海底4 m深的麻坑，计算表明在水合物层之下至少需要有22 m厚的游离气层。苏正和陈多福[4]计算了布莱克海台997站位的水合物和游离气体积分数分布，在水合物稳定带底界之下26 m处的气体饱和度为28％，底界之下74 m处气体饱和度为％，其中水合物体积分数分布与同一区域的ODP995站位是相近的[24]。28％的气体饱和度大于气体流动所需20％的饱和度，而底界之下74 m处％的气体饱和度不能流动，也不能传递孔隙气体压力。如果20％的饱和度指示可传递气层的底界，则气层的有效压力传递厚度约为30 m，这与笔者22 m厚的游离气层模型计算结果相近（图5）。实际上，该钻位水合物平均体积分数约为6%[4]，可封闭气层厚度为24 m（三角点所示），接近模型估计的22 m。此外，在水合物稳定带底部的水合物饱和度达24%[4]，其毛细管作用可封闭约33 m的游离气层（菱形点所示），与Flemings等[25]估计的极限破坏厚度29 m相似（虚线所示位置），接近但略小于30 m的参考厚度。然而，在996站位游离气发生泄漏后， 997站位扩散型水合物的体积分数仍在持续增加[26]，水合物层的封闭能力也相应增强，游离气层厚度不断增长，因此，997站位游离气厚度（30 m）大于996站位游离气发生泄漏时的22 m气层厚度是合理的。
图5 布莱克海台的水合物饱和度和所能封闭的游离气层厚度
气层厚度随水合物饱和度增加而增高，水平虚线与气层厚度曲线的交点（29 m）为Flemings等预测的997站位气层的临界水力压裂厚度[25]，圆形点标示约30 m的实际气层厚度，三角形点显示平均饱和度6％的水合物能封闭24m的气层，而饱和度24％的水合物可封闭33 m的游离气层（菱形点）
4 结语
本文构建了水合物层下伏游离气渗漏动力学过程的数学模型，游离气被水合物层的毛细管作用所圈闭，下伏游离气的超压随游离气层的增长而增加；当气体超压超过作用于水合物与游离气层界面的毛细管阻力时，游离气渗漏进入上覆水合物稳定带，并以“活塞式”驱动上覆孔隙水向外排出，渗漏速度随水流柱高度的减小而增加；当水流阻抗超过相应层段的静岩压力时沉积体变为流沙，流沙沉积被海流带走便在海底留下麻坑。模型显示麻坑深度为游离气层厚度和水合物封闭层埋深的函数，而与沉积介质的渗透率无关。游离气渗漏形成的海底麻坑对水合物下伏游离气层的厚度具有指示作用，在已知水合物封闭层深度和海底麻坑深度条件下，模型可以计算水合物层下伏游离气藏发生渗漏时的气层厚度，在布莱克海台海底发育有4 m深的麻坑，它的形成需要至少22 m厚的游离气层。
致谢：挪威国家石油公司Martin Hovland教授提供了全球麻坑基础资料和最新信息，表示感谢。
参考文献
[1]Xu W ,Ruppel the Occurrence,Distribution,and Evolution of Methane Gas Hydrate in Porous Marine Sediments[J].Journal of Geophysical Research,1999,104:5081-5095.
[2]Davie M K,Buffett B Steady State Model for Marine Hydrate Formation:Constraints on Methane Supply from Pore Water Sulfate Profiles[J].Journal of Geophysical Research,2003,108(B10): 2495,doi:.
[3]Chen Duo-Fu,Su Zheng,Cathles L of Gas Hydrates in Marine Environments and Their Thermodynamic Characteristics[J].Terrestrial Atmospheric and Oceanic Sciences,2006,17(4) :723-737.
[4]苏正，陈多福.海洋环境甲烷水合物溶解度及其对水合物发育特征的控制[J].地球物理学报，2007,50(5）: 1518-1526.
[5]Trehu A M,Bohrmann G,Rack F R,et ofthe Ocean Drilling Program,Initial Reports Volume 204[M].TX:Ocean Drilling Program,2003.
[6]Paull C K,Spiess F N,Ussler W Ⅲ，et Plumes on the Carolina Continental rise: Associations with Gas Hydrates[J].Geology,1995,23: 89-92.
[7]Sahling H,Bohrmann G,Spiess V,et in the Northern Congo Fan area,SW Africa: Complex seafloor features shaped by Fluid Flow[J].Marine Geology,2008,249 : 206-225.
[8]Gay A,Lopez M,Berndt C,et Controls on Focused Fluid Flow Associated with Seafloor Seeps in the Lower Congo Basin[J].Marine Geology,2007,244 (1/2/3/4):68-92.
[9]Hovland M,Svensen H,Forsberg C F,et Pockmarks with Carbonate-Ridges off Mid-Norway:Products of Sediment Degassing[J].Marine Geology,2005,218:191-206.
[10]陈多福，李绪宣，夏斌.南海琼东南盆地天然气水合物稳定域分布特征及资源预测[J].地球物理学报，2004,47(3）:483-489.
[11]Cathles L in Sub-Water Table Fluid Flow at the End of the Proterozoic and Its Implicationsfor Gas Pulsars and MVT Leadzinc Deposits[J].Geofluids,2007,7(2): 209-226.
[12]Bear of Fluids in Porous Media[M].New York:Elsevier,1972.
[13]Vigil G,Xu Z,Steinberg S,et of Silica Surfaces[J].J Colloid Interface Sci,1994,165:367.
[14]Cathles L Seals as a Cause of Pressure Com-partmentation in Sedimentary Basins:Presented at the Gulf Coast Section SEPM Foundation 21 st Annual Research Conference on Petroleum Systems of Deep-Water Basins,2001:561-572.
[15]Shosa J D,Cathles L Investigation of Capillary Blockage of Two-Phase Flow in Layered Porous Media,in Petroleum Systems of Deep-Water Basins: Global and Gulfof Mexico Experience: Proceedings ofthe GCSSEPM st Annual Bob Research Conference,2001:725-739.
[16]苏正，陈多福.海洋天然气水合物的类型及特征[J].大地构造与成矿学，2006,30(2).
[17]Hovland M,Judd A Pockmarks and on Geology,Biology and the Marine Ltd.,1988.
[18]Hovland M,Svensen Pingoes:Indicators of Shallow Gas Hydrates in a Pockmark at Nyegga,Norwegian Sea[J].Marine Geology,2006,228:15-23.
[19]Gay A,Lopez M,Cochonat P,et so1ated Seafloor Pockmarks Linked to BSRs,Fluid Chimneys,Polygonal Faults and Stacked Oligocene-Miocene Turbiditic Palaeochannels in the Lower Congo Basin[J].Marine Geology,2006,226(1/2):25-40.
[20]Nicholl M J,Karnowski Apparatus for the Demonstration of Quicksand[J].Journal of Geoscience Education,2006,54(5): 578-583.
[21]Matsumoto R,Paull C,Wallace Leg 164 Scientific party[C]//Gas hydrate sampling on the Blake Ridge and Carolina Rise: ODP,Leg 164 Preliminary Report,1996.
[22]Taylor M H,Dillon W P,Pecher I and Migration of Methane Associated with the Gas Hydrate Stability Zone at the Blake Ridge Diapir:New Insights from Seismic Data[J].Marine Geology,2000,164:79-89.
[23]Paull C K,Matsumoto R,Wallace P J,et [C].164:TX:Ocean Drilling Program,2000.
[24]王秀娟，吴时国，刘学伟.天然气水合物和游离气饱和度估算的影响因素[M].地球物理学报，2006,49(2）:504-511.
[25]Flemings P B,Liu X,Winters W Pressure and Multiphase Flow in Blake Ridge Gas Hydrates[J].Geology,2003,31: 1057-1060.
[26]Liu X,Flemings P Multiphase Flow Model of Hydrate Formation in Marine Sediments[J].Journal of Geophy sical Resea rch,2007,112,B03101,doi:.
罗伟东 李 刚
（广州海洋地质调查局 广州 510760）
作者简介：罗伟东，男，广东人，硕士，高级工程师，主要从事海洋地形地貌调查和研究工作。
E-mail：loky＿luo＠，电话：，。
摘要 本文介绍了地形地貌研究主要工具软件的基本功能，针对主要工具软件的特点和实际应用情况，从地形数据的浏览、编辑、提取和转换，数据的网格化，等深线和彩色阴影底图等的编辑，调色板的制作与获取，地形图和三维立体图的制作等方面探讨地形地貌研究工具软件应用的方法，为地形地貌研究工作的开展创造有利条件，大大提高了工作效率和质量。
关键词 地形地貌 工具软件 应用方法
1 前言
海洋地形地貌调查和研究是大多数海洋调查研究的基本内容和重要组成部分。我国自20世纪90年代初引进多波束测深系统后开展了我国近海和邻近海域多波束海底地形地貌调查和综合研究，积累了丰富的资料和研究成果。提高地形地貌研究水平和工作效率是研究人员不断探索的方向。地形地貌研究工作常用的工具软件发挥了重要作用，能进行原始数据处理、数据网格化、数据提取转换浏览、制图等，还能进行地形坡度、距离和面积计算，地形剖面的提取制作，方便对海底地形分区、分析分布特征和识别地貌单元等。熟练掌握地形地貌研究工作常用的工具软件是研究人员必须具备的基本技能，合理选择和恰当应用工具软件是提高工作效率的有效途径。目前，地形地貌研究可用的工具软件有许多，如MapGIS、ArcGIS、Golden Software Surfer、CorelDRAW、Global Mapper、Fledermaus、GMT、MB-System和Caris Hips and Sips等。我国的地形地貌研究人员在以往的使用过程中，对相关软件积累了一些经验，但是面对种类繁多的工具软件，要恰当合理综合应用以提高工作效率和质量却不是一件容易的事。在实践工作中，常常会因为没有掌握软件的使用方法或没有发挥各软件的特点，而遇到一些困难，降低工作效率和质量。作者通过充分应用软件已有的功能，开发其潜在功能，结合工作实际，扬长避短，综合应用，使相应软件更好地服务地形地貌调查和研究工作，发挥更大的作用。
2 地形地貌研究主要工具软件概述
MapGIS、ArcGIS和Golden Software Surfer功能概述
这三款软件，特别是MapGIS和Surfer在国内得到广泛应用，相关文章对用法和功能的介绍很多。地形地貌研究主要的应用是利用MapGIS、ArcGIS制作规范的地形图和地貌图等；Surfer主要的应用是制作坡度图和地形图、地貌图和三维立体图等。
GMT功能概述
GMT（Generic Mapping Tools）是一个通用地学制图工具软件，发展至今已在大气、海洋、地震等研究领域得到较为广泛的应用。GMT是全命令行软件，需要用户录入命令及其指定的输入输出数据和各种配置参数（如经度、纬度、颜色配置等）来进行操作。GMT虽然操作相对繁琐，但执行效率很高，而且图形文件输出的主要格式为PostScript文件，这种格式可提供高质量、跨平台的图形打印，这也是GMT得以广泛应用的原因。
实际应用中，主要是制作二维三维地形图。利用GMT强大的颜色渲染和标准图框输出功能，制作颜色逼真的彩色阴影图、报告和论文插图，大数据量的网格化功能也是其强项。
Global Mapper功能概述
Global Mapper（简称GM）是美国Global Mapper公司开发的一个简单、实用的图形管理与应用软件。它可以在指定投影和地理坐标的基础上，编辑、转换光栅和矢量地形图，绘制二维、三维地形图以及点、线自编图形，具有良好的绘图、编辑、显示以及数据输出界面。其主要功能特点和实际应用包括以下四个方面：
1）浏览、合成、输入、输出、显示大部分流行的栅格图形和高程及矢量数据
实际应用中，主要用来直接打开浏览各类型如Surfer和GMT网格文件；输入输出（特定范围）网格、XYZ文件及点、线文件；合成相同或不同格式、大小的数据（网格、XYZ）文件。
2）具有数据、图形的转换、编辑、拼接、打印功能和投影转换
主要进行数据、图形的格式转换，如将GMT网格文件转换成Surfer网格文件，重新网格化数据，将XYZ文件转换成网格文件，网格、XYZ文件转换（输出）成图片等。在网格文件中选取制作特定位置的地形剖面。进行坡度、距离和面积计算。结合不同地理位置的地图，设置和改变、转换图片的投影方式。
3）具有简单的地理信息功能和动态GPS接收功能
为没有位置信息的图片添加坐标信息，在同一个项目中打开的文件都会在相同投影下的大地坐标系中显示。
4）在实际应用中，GM也存在不足之处：等深线功能简单，不能调整标注间距和字体样式及大小；地形图图框、比例尺不规范，文字大小和标注不可调节。
Fledermaus功能概述
加拿大IVS 3D公司生产的Fledermaus是全球海测数据3D具体化的领先商业软件之一，是一套功能强大的交互式三维数值数据虚拟实境的系统。它可以帮助使用者完成包括海洋（海岸、海底）资源调查与制图、环境影响评估、采矿、地质调查以及各种研究等工作。Fledermaus可直接支持广泛的工业数据格式输入。可直接导入、显示数字地形图，点、线、多边形数据集合、卫星影像并进行分析。浏览器iView4D可以随时浏览处理过或分析后的数据结果。
实际应用中，主要是利用其三维可视化功能和复杂区域编辑模块，用于导入网格文件，对三维数据进行三维可视化显示和操作，可任意角度观察地形地貌特征，可方便进行海底地形分区、分析分布特征和识别地貌单元。
3 工具软件的综合应用
常用工具软件在地形地貌研究中发挥了重要作用，工具软件种类繁多，各有优点缺点，科研人员要全部掌握其功能和操作方法有一定难度也没有必要，应该根据实际工作需要，重点掌握几个主要软件，扬长避短，合理应用，其他软件的应用能满足工作需要即可。针对所从事的工作，主要的应用情况从以下几个方面进行阐述。
图1 GM软件主界面窗口
地形数据的浏览、编辑、提取和转换
海底地形分区、分析分布特征和识别地貌单元等是地形地貌的工作之一，这些工作都要对地形数据进行浏览、观察、描述，以往一般是通过打印的地形图和三维地形图等图件来进行识别和描述。随着计算机软硬件技术的发展，目前，更科学的方法是利用工具软件完成上述工作，也可结合图件打印的方法。利用工具软件观察更直观简便高效，可直接导入、显示数字地形图，调节着色、等值线间距和垂直比例，进行地形坡度、距离和面积计算，特定位置剖面显示提取和制作，多角度、随意缩放、显示任意位置水深和点线多边线标注等。GM软件在这方面具有强大的功能（图1），结合Fledermaus，Golden SoftwareVoxler软件能方便地实现上述常用功能，大大提高工作效率，非常适合地形地貌研究和海洋调查使用。
GM软件可方便地打开XYZ和多种格式的网格文件（包括大数据量文件），还能把数据转换成自己的网格格式Global mapper grid显示，文件小，占用内存少，可随意进行操作，如放大缩小编辑等，显示效果好，这是GM软件的主要特点之一，相比其他软件具有明显优势。显示二、三维地形图，具有标注等深线、测线、标注符号和简单的3D view功能，能从各个角度浏览地形；能实时显示当前或指定位置的水深和经纬度；具备生成地形剖面功能，可选取任意位置进行地形剖面提取和制作；能进行地形坡度、距离和面积计算，地形坡度的计算和显示是在应用地形剖面功能时，在地形剖面窗口中的Options/showpath detais查看。地形剖面可输出为CSV或XYZ数据和图片格式，CSV或XYZ数据可用Grapher等专业软件制作成剖面图；图片格式可直接在CorelDRAW中导入编辑，为保证输出后的图片质量和字体大小合适，字体编辑时要注意调节字体的大小，输出的图片长度应在2000像素左右，输出后在word中的字体大小约为小5号。剖面图如图2所示。
微地形地貌的浏览和描述方面，需要用到相关三维可视化系统，GM软件具有简单的3D view功能，并不能满足实际应用的要求，目前最好的三维可视化系统是Fledermaus，Golden Software Voxler也有不错的表现，实际应用时应结合使用。
图2 典型地形剖面
数据的网格化
随着多波束海底地形地貌调查和综合研究的技术发展，对研究人员的要求也越来越高，数据网格化和数据转换方法已经成为研究人员必须掌握的技能之一。过去认为数据的网格化是由专业处理人员完成的想法已经不符合实际。研究人员在进行研究时不是简单地利用网格文件，实际可能应用的情况：简单的数据编辑；数据（网格文件）的格式转换；重新网格化，生成不同网格间距的网格文件；提取特定区域网格数据等。
在进行XYZ数据的网格化时，不同的网格化（插值）方法的绘图效果是不一样的。一般地，各软件都会提供多种网格化方法供用户选择，如Surfer提供多达12种方法、GMT有4种方法和ArcGIS提供5种方法。合理选择网格化方法和工具软件能有效提高图件制作质量，地形地貌的制图具体使用的网格化方法，需要根据客观环境特征和数据本身的特点，进行相应的数据分析，才能绘制出准确、有意义的图件。一般地，中大数据量的数据用加权反距离法（Inverse distance to a power）或者最小曲率法效果好、效率高，小数据量数据可使用克里格法（Kriging）。单一的工具软件不能满足实际要求，Surfer软件只适应于网格数据量较小的数据（几百兆的数据文件），Surfer的工作效率很低；大数据量的数据一般要选用GMT或MB-System软件进行网格化；进行数据浏览转换时可使用GM软件。
GMT网格化采用gmtmbgrid和Surface命令；MB-System软件采用mbgrid命令；Sur-face是张性样条网格化法，是一种经过改良的标准最小曲率算法，允许用户把张量引入表面。执行surface程序前要对数据进行预处理，消除混淆现象。地形数据采用block-median对网格间隔框内的值返回中值。Gmtmbgrid算法跟mbgrid命令一样是采用加权反距离插值法。具体根据实际需要选择。MB-System和GMT软件一般配合使用。使用的具体命令如下：
blockmedian文件名.xyz -R113/120/11/13 -I200e -V ＞ ship＿
Surface文件名.xyz -R113/120/11/13 -I200e -G文件名.grd
mbgrid -I文件名 -0文件名.grd -R113/120/11/13 -A2 -C2 -E 200/200 -N -V
GM软件网格化采用三角剖分法（triangulation/liner interpolation），这种方法对大数据量的数据处理效率低，数据边缘效果稍差，具体操作流程：打开XYZ文件，出现Ge-neric ASCII Text File Import Options窗口，在Import type选项中，选择Elevation Grid from3D Point Data，其他选项可以采用默认。确定后，出现Elevation Grid Creation Options窗口，在Grid Spacing中可自定义网格大小；确认后即可完成网格化过程，然后在File/Ex-port Elevation grid format，选择合适的格式输出数据。图3是不同软件网格化及制图效果图。左图是用GM软件采用三角剖分法网格化绘制而成，中图是GMT软件用改良的标准最小曲率算法（surface命令）进行网格化，右图是Surfer采用加权平均算法；从图中可清楚地看出，三个软件绘制的地形图一致性很好；GMT软件在网格化效果和调色效果方面都是最出色的，能较好的反映微地形且色彩逼真。
等深线和彩色阴影底图等的编辑
多波束数据经过后处理后，数据质量一般都能达到要求，特殊情况时也会有局部小范围数据质量较差的数据，在制作成图件时，在彩色阴影图和等深线上都会有直接的反映。一般地，都要对其进行编辑，去伪存真。彩色阴影底图的制作一般用GMT或Surfer，彩色阴影底图的编辑一般选用Photoshop或者CorelDRAW，等深线的生成可用MapGIS、Surfer和GMT软件，MapGIS具有等深线编辑功能，也可以在CorelDRAW中进行编辑，效果都很不错。在CorelDRAW中进行编辑时，等深线一般是在Surfer生成*.emf格式文件，在Corel-DRAW中打开进行编辑，删除等深线时应选用虚拟段删除工具，CorelDRAW是把显示的整条线段定义为由多条虚拟段组成，用虚拟段删除工具选中整条线段或线段的一部分后，程序将会自动删除选中的部分，操作很方便。添加线段时选中手绘工具，调节合适粗细，进行曲线和直线绘制。图片输出时同样要注意像素的调节，确保图像质量和字体合适。图4为等深线和彩色阴影底图编辑前后对比图，底图是用Surfer制作的，编辑方法如上述。右图为编辑后的情况，效果非常明显。
图3 不同软件网格化及制图效果图
图4 等深线和彩色阴影底图编辑前后对比图
调色板的制作与获取
制图的美观程度主要取决于调色板的着色效果，制作美观实用的图件是科研工作者的不懈追求。根据实际需要选择合适的工具软件及制作或引用合适的调色板是制图工作的重中之重。不同的工具软件的调色效果是有区别的。实践应用中，GMT软件的调色效果最好，Surfer、MapGIS和ArcGIS调色效果稍差，GM的调色效果次之。这也是GMT软件被广泛应用的主要原因之一。因此，在实际应用中，正式图件及成果报告插图等可应用GMT软件制图或制作彩色阴影图。GMT软件自身提供了许多独有且全球流行的调色板，自身也提供相应的命令如makecpt和grd2cpt根据相应的调色板模板对特定数据进行调色板制作，但并不能满足实际应用的需求。因此，制作或者获取第三方调色板是可行的办法。
调色板模板制作与管理的有效可行办法：首先，根据工作实际，确定工作中常用的几种地形数据类型如海陆结合数据、浅水数据和深水数据等，每种数据类型制作几个常用的调色方案，所有的调色方案都做成各软件特定的格式，可直接或稍为更改即可引用。由于GMT是全命令行软件，调色板的制作并不像Surfer等友好界面软件那么直观方便。一般地，在利用GMT自带调色板模板稍加修改或者直接在Surfer软件中制作调色板后，导出后再转换为GMT格式的调色板。
获取调色板的方法有两种：一种为通过专业网站获取相关调色板，专业网站提供多个领域的世界上主流的调色板和大量用户自主设计的调色板。另一种方法是手头只有地形图图片而没有相应数据调色板的情况，这时最好的方法是利用相关软件如Photoshop或其他具有取色功能的软件进行调色板制作，着色效果跟原有图片基本一致。
地形图和三维立体图的制作
发挥各地学绘图软件的优势，结合应用才能做出规范且最佳效果的图件。地形图和三维立体图正式图件的制作，单一工具软件不能满足要求。Surfer软件功能较全面但特点不明显且没有投影功能，等深线制作和编辑功能较弱；GMT的网格化和调色功能强大，等深线制作也不错，但没有编辑功能；由于制作正式图件时，可采用由GMT软件制作地形图和三维立体图的彩色阴影底图，再导入MapGIS或ArcGIS制成图件的方法，制作时要注意投影方式的匹配问题，GMT制作彩色阴影图时应选择合适的投影，导入前应该把PS文件转换成BMP或JPG图片，且要根据打印图件的大小选择合适的输出像素。图5为西太平洋地形图，是由GMT软件制作彩色阴影图导入MapGIS进行中文标注等。正式图件的等深线可直接由MapGIS生成并编辑或由GMT生成导入等。
4 结 语
地形地貌研究常用的工具软件种类多、功能很强大。本文针对常用工具软件的特点和实际应用情况，总结了这套工具软件应用方法的主要特点：把各软件的优势应用到实际工作，扬长避短，最大限度服务地形地貌研究工作，提高工作效率和质量；根据实际应用制定相应操作流程，综合应用，简化操作，提高效率；挖掘新工具软件，把GM软件应用于地形地貌研究工作，起到非常好的应用效果。这套方法为地形地貌研究工作的开展创造有利条件，大大提高了工作效率和质量。
图5 西太平洋三维立体图
参考文献
［1］陈欢欢等.Surfer 等值线绘制中的十二种插值方法.工程地球物理学报，2007，2.
［2］刘方兰等.Global Mapper系统在海洋调查中的应用.海洋技术，2011，01.
The application method study of the general-used software for studying the topographic and geomorphic features
Luo Weidong Li Gang
（Guangzhou Marine Geological Survey，Guangzhou，510760）
Abstract：The paper introduces the basic function of the general-used software for studying thetopographic and geomorphic on the characteristic features and the comprehen-sive application in practical experience，talking about the usages of the general-used software forstudying the topographic and geomorphic features on scanning，editing，collecting and conversionof the geographic data，data gridding，edit for the bathymetric contour and the color shader basemaps，palette make and abtain，topographic map and Creates the advantages for studying the top-ographic and geomorphic features，and also increases the efficiency and the quality of the study.
Key words：Topographic and geomorphic features Software Application method
［1］华安增，孔圆波，李世平等.岩块降压破碎的能量分析.煤炭学报，1995，20（4）：389～392
［2］陈颙，姚孝新，耿乃光.应力途径、岩石的强度和体积膨胀.中国科学，1979，11：1093～1100
［3］姚孝新，耿乃光，陈颙.应力途径对岩石脆性和延性的影响.地球物理学报，1980，23（3）：312～319
［4］周火明，盛谦，李维树等.三峡船闸边坡卸荷扰动区范围及岩体力学性质弱化程度研究.岩石力学与工程学，2004，23（7）：1078～1081
［5］哈秋舲，李建林，张永兴等.节理岩体卸荷非线性岩体力学.北京：中国建筑工业出版社，1998
［6］吴刚，孙均.卸荷应力状态下裂隙岩体的变形和强度特性.岩石力学与工程学报，1998，17（6）：615～621
［7］张有天，周维垣.岩石高边坡的变形与稳定.北京：中国水利水电出版社，～247，303～333
［8］吴玉山.应力路径对凝灰岩力学特性的影响.岩土工程学报，1983，5（1）：112～120
［9］李天斌，王兰生.卸荷应力状态下玄武岩变形破坏特征的试验研究.岩石力学与工程学报，1993，12（4）：321～327
［10］Lau Joseph S O，Chandler N laboratory .，2004，41（8）：1427～1445
［11］汪浩译（Farmer I W 著）.岩石的工程性质.徐州：中国矿业大学出版社，～54
［12］尤明庆，华安增.卸围压法测量岩石材料的泊松比.实验力学，1997，12（2）：274～278
［13］朱之芳等.刚性试验机.北京：煤炭工业出版社，～53
［14］李通林，谭学术，刘伟.矿山岩石力学.重庆：重庆大学出版社，
［15］Byerlee J of ，116：615～626
［16］Rice J R，Ruina A of steady frictional .，1983，50：343～349
［17］Maksimovic shear strength components of a rough rock .，1996，33（8）：769～783
［18］Shimamoto pressure reduction .，1985，22（4）：227～236
［19］齐庆新，史元伟，刘天泉.冲击地压粘滑失稳机理的试验研究.煤炭学报，1997，22（2）：144～148
［20］尤明庆，华安增.岩石试样的强度准则及内摩擦系数.地质力学学报，2001，7（1）：53～60
［21］Swanson S R，Brown W observation of loading path independence of fracture in .，1971，8（3）：227～231
［22］尤明庆，华安增.应力路径对岩样强度的影响.岩土工程学报，1998，20（5）：101～104
［23］葛修润，卢应发.循环载荷作用下岩石疲劳破坏和不可逆变形问题的探讨.岩土工程学报，1992，14（3）：56～61
［24］耿乃光，许东俊.最小主应力减小引起岩石破坏时中间主应力的影响.地球物理学报，1985，28（2）：191～196
［25］许东俊，耿乃光.中等主应力变化引起的岩石破坏与地震.地震学报，1984，6（2）：159～165
［26］You M，Hua A，Li of loading path on supporting capacity and deformation property of ：Frontiersof Rock Mechanics and Sustainable Development in the 21st Century（2001 ISRM and 2nd ARMS）.Netherlands：A A Balkema，～59
［27］周小平，哈秋舲，张永兴等.峰前围压卸荷条件下岩石的应力-应变全过程分析和变形局部化研究.岩石力学与工程学报，2005，24（18）：3236～3245
［28］陈忠辉，林忠明，谢和平.三维应力状态下岩石损伤破坏的卸荷效应.煤炭学报，2004，29（1）：31～35
［29］任建喜，葛修润，蒲毅彬等.岩石卸荷损伤演化机理CT 实时分析初探.岩石力学与工程学报，2000，19（6）：697～701
［30］徐松林，吴文，王广印等.大理岩等围压三轴压缩全过程研究Ⅰ：三轴压缩全过程和峰前、峰后的卸围压全过程实验.岩石力学与工程学报，2001，20（6）：763～767
［31］尤明庆.复杂路径下岩样的强度、变形特性.岩石力学与工程学报，2002，21（1）：23～28
［32］尤明庆，勾攀峰.岩样三轴压缩的对角破坏和残余承载能力.矿山压力与顶板管理，2002，19（1）：98～100
［33］尤明庆，苏承东，徐涛.岩石试样的加载卸载过程及杨氏模量.岩土工程学报，2001，23（5）：588～592
［34］Tadokoro K，Ando for rapid fault healing derived from temporal changes in S wave .，2002，29（10）：1029
［35］Li Y G，Vidale J E，Day S M et fault healing on the rupture zone of the 1999 M Hector Mine，California earthquake .，2003，93：854～869
［36］Marone friction laws and their application to seismic .，1998，26：643～696
［37］Bos B，Spiers C healing processes in gouge-bearing faults：Insights from experiments on a rock analogue .，2002：2537～2566
［38］Yasuhara H，Marone C，Elsworth zone restrengthening and frictional healing：the role of ofGeophysical research，110，B06310，doi：
［39］Schenk O，Urai J L，Evans effect of water on recrystallization behavior and grain boundary morphology in calcite observations of nural marble of Structural Geology，2005，27（10）：1856～1872
［40］王仁.塑性力学基础.北京：科学出版社.～5
［41］杨显杰，罗艳，高庆，蔡力勋，刘宇杰.循环软化45 碳钢和循环硬化304 不锈钢的棘轮行为实验研究.固体力学学报，2005，26（2）：125～131
［42］尤明庆，苏承东.大理岩试样循环加载强化作用的试验研究.固体力学学报，2007，28（印刷中）
［43］苏承东，杨圣奇.岩样循环加载条件下强度和变形特征试验.河海大学学报，2006，34（6）：667～671
［44］尤明庆.两种晶粒大理岩力学性质的研究.岩土力学，2005，26（1）：91～96
［45］尤明庆.基于粘接、摩擦承载特性对岩石变形破坏的研究.地质力学学报，2005，11（3）：288～292
［46］陈颙，黄庭芳.岩石物理学.北京：北京大学出版社，～84
［47］葛修润，蒋宇，卢允德等.周期荷载作用下岩石疲劳变形特性试验研究.岩石力学与工程学报，2003，22（10）：1581～1585
［48］李海波，赵坚，李俊如等.基于裂纹扩展能量平衡的花岗岩动态本构模型研究.岩石力学与工程学报，2003，22（10）：1683～1688
［49］谢和平，鞠杨，黎立云.基于能量耗散与释放原理的岩石强度与整体破坏准则.岩石力学与工程学报，2005，24（17）：3003～3010
［50］蔡美峰，王金安，王双红.玲珑金矿深部开采岩体能量分析与岩爆综合预测.岩石力学与工程学报，2001，20（1）：38～42
［51］Hua A，You failure due to energy release during unloading and application to underground rock burst and Underground Space Technology，2001，16：242～246
［52］尤明庆，华安增.岩石试样破坏过程的能量分析.岩石力学与工程学报，2002，21（6）：778～781
［53］尤明庆，华安增，李玉寿.缺陷岩样的强度和变形特性.岩土工程学报，1998，20（2）：98～102
SCI论文是指被科学引文索引数据库（Science Citation Index, 简称 SCI）收录的论文。一般是将论文发表在被科学引文索引数据库收录的期刊，这些期刊被称为SCI期刊，SCI期刊发表您的论文后会将您的文章推送给SCI数据库，这样您的文章就被SCI数据库收录
属于社会调查性论文最好方法是设计调查问卷,根据调查结果,数据,进行分析.调查要接触到社会的各个层面,才有权威性.怎样写论文？一、前言部分前言部分也常用“引论”、“概论”、“问题背景”等做标题，在这部分中，主要介绍论文的选题。二、综述部分任何一个课题的研究或开发都是有学科基础或技术基础的。通过与中外研究成果的比较和评论，说明自己的选题是符合当前的研究方向并有所进展，或采用了当前的最新技术并有所改进，目的是使读者进一步了解选题的意义。三、方案论证在与文献资料中的方案进行比较时，首先要阐述自己的设计方案，说明为什么要选择或设计这样的方案，前面评述的优点在此方案中如何体现，不足之处又是如何得到了克服，最后完成的工作能达到什么性能水平，有什么创新之处（或有新意）。如果自己的题目是总方案的一部分，要明确说明自己承担的部分，及对整个任务的贡献。四、论文主体前面三个部分的篇幅大约占论文的1/3，主体部分要占2/3左右。在这部分中，要将整个研究开发工作的内容，包括理论分析、总体设计、模块划分、实现方法等进行详细的论述。论文主体部分的内容一般要分成几个章节来描述。在写作上，除了用文字描述外，还要善于利用各种原理图、流程图、表格、曲线等来说明问题，一篇条理清晰，图文并茂的论文才是一篇好的论文。六、结束语这一节篇幅不大，首先对整个论文工作做一个简单小结，然后将自己在研究开发工作中所做的贡献，或独立研究的成果列举出来，再对自己工作的进展、水平做一个实事求是的评论。但在用“首次提出”、“重大突破”、“重要价值”等自我评语时要慎重。七、后记在后记中，主要表达对导师和其他有关教师和同学的感谢之意。对此，仍要实事求是，过分的颂扬反而会带来消极影响。这一节也可用“致谢”做标题。八、参考文献中外文的参考文献应按照规范列举在论文最后。这一部分的编写反映作者的学术作风。编写参考文献要注意：（1）要严格按照规范编写，特别是外文文献，不要漏写、错写；（2）论文内容和参考文献要前后对应，正文中凡引用参考文献的地方应加注；（3）列出的文献资料应与论文课题相关，无关的文献只会使读者感到你的研究目标很分散；（4）选择的参考文献应主要是近期的。字数并不是最重要的，关键是论文的质量。摘 要1 摘要的概念和作用其基本要素包括研究目的、方法、结果和结论。具体地讲就是研究工作的主要对象和范围采用的手段和方法，得出的结果和重要的结论有时也包括具有情报价值的其它重要的信息。摘要应具有独立性和自明性，并且拥有与文献同等量的主要信息，即不阅读全文，就能获得必要的信息。主要功能：1) 让读者尽快了解论文的主要内容，以补充题名的不足。现代科技文献信息浩如烟海，读者检索到论文题名后是否会阅读全文，主要就是通过阅读摘要来判断；所以，摘要担负着吸引读者和将文章的主要内容介绍给读者的任务。2) 为科技情报文献检索数据库的建设和维护提供方便。33 摘要的写作注意事项1) 摘要中应排除本学科领域已成为常识的内容；切忌把应在引言中出现的内容写入摘要；一般也不要对论文内容作诠释和评论(尤其是自我评价)。2) 不得简单重复题名中已有的信息。比如一篇文章的题名是《几种中国兰种子试管培养根状茎发生的研究》，摘要的开头就不要再写：“为了……，对几种中国兰种子试管培养根状茎的发生进行了研究”3) 结构严谨，表达简明，语义确切。要按逻辑顺序来安排。句子之间要上下连贯，互相呼应。摘要慎用长句，句型应力求简单。每句话要表意明白，无空泛、笼统、含混之词，但摘要毕竟是一篇完整的短文，电报式的写法亦不足取。摘要不分段4) 用第三人称。建议采用“对……进行了研究”、“报告了……现状”、“进行了……调查”等记述方法标明一次文献的性质和文献主题，不必使用“本文”、“作者”等作为主语。5) 要使用规范化的名词术语，不用非公知公用的符号和术语。新术语或尚无合适汉文术语的，可用原文或译出后加括号注明原文。6) 除了实在无法变通以外，一般不用数学公式和化学结构式，不出现插图、表格。7) 不用引文，除非该文献证实或否定了他人已出版的著作。8) 缩略语、略称、代号，除了相邻专业的读者也能清楚理解的以外，在首次出现时必须加以说明。目前摘要编写中的主要问题有：要素不全，或缺目的，或缺方法；出现引文，无独立性与自明性；繁简失当。2 文献分类工作的程序和方法文献分类工作程序——按分类工作进行的先后次序安排的工作步骤。如果不按文献分类工作的程序进行，将造成重复劳动，并容易出现差错。文献分类工作的程序是：分析文献主题→归类→给号→校核。总体原则要求1、立论客观，具有独创性：所提出的问题在本专业学科领域内有一定的理论意义或实际意义，并通过独立研究，提出了自己一定的认知和看法。2、论据翔实，富有确证性：论文能够做到旁征博引，多方佐证，所用论据自己持何看法，有主证和旁证。论文中所用的材料应做到言必有据，准确可靠，精确无误。3、论证严密，富有逻辑性：作者提出问题、分析问题和解决问题， 要符合客观事物的发展规律，全篇论文形成一个有机的整体， 使判断与推理言之有序， 天衣无缝。4、体式明确，标注规范：论文必须以论点的形成构成全文的结构格局，以多方论证的内容组成文章丰满的整体，以较深的理论分析辉映全篇。此外，论文的整体结构和标注要求规范得体。5、语言准确、表达简明：论文最基本的要求是读者能看懂。 因此， 要求文章想的清，说的明，想的深，说的透，做到深入浅出，言简意赅。标准格式1、论文题目：（下附署名）要求准确、简练、醒目、新颖。2、目录目录是论文中主要段落的简表。（短篇论文不必列目录）3、内容提要：（摘要）是文章主要内容的摘录，要求短、精、完整。字数少可几十字，多不超过三百字为宜。4、关键词或主题词关键词是从论文的题名、提要和正文中选取出来的，是对表述论文的中心内容有实质意义的词汇。关键词是用作计算机系统标引论文内容特征的词语，便于信息系统汇集，以供读者检索。 每篇论文一般选取3-8个词汇作为关键词，另起一行，排在"提要"的左下方。主题词是经过规范化的词，在确定主题词时，要对论文进行主题分析，依照标引和组配规则转换成主题词表中的规范词语。（参见《汉语主题词表》和《世界汉语主题词表》）。5、论文正文：（1）引言：引言又称前言、序言和导言，用在论文的开头。 引言一般要概括地写出作者意图，说明选题的目的和意义, 并指出论文写作的范围。引言要短小精悍、紧扣主题。〈2）论文正文：正文是论文的主体，正文应包括论点、论据、 论证过程和结论。主体部分包括以下内容：a.提出问题-论点；b.分析问题-论据和论证；c.解决问题-论证方法与步骤；d.结论。6、参考文献一篇论文的参考文献是将论文在研究和写作中可参考或引证的主要文献资料，列于论文的末尾。参考文献应另起一页，标注方式按《GB7714-87文后参考文献著录规则》进行。中文：标题--作者--出版物信息（版地、版者、版期）英文：作者--标题--出版物信息所列参考文献的要求是：（1）所列参考文献应是正式出版物，以便读者考证。（2）所列举的参考文献要标明序号、著作或文章的标题、作者、出版物信息。以下例1—例5为通行的参考文献著录格式，各学员可根据不同情况分别选用。例1：张筑生．微分半动力系统的不变集：[学位论文]，北京：数学系统学研究所，1983例2：傅承义，陈运泰，祁贵中．地球物理学基础．北京：科学出版社，1985，447例3：华罗庚，王元．论一致分布与近似分析．中国科学，1973（4）：339～357例4：翟伦玉，陈亚光，彭军等．高等学校化学学报，1991，12（8）：1009例5：Borko H，Bernier C L．Indexing concepts and methads．New York： Academic Pr， 1978打印要求用计算机打印（16K纸；正文用四号楷体字，每页20行、每行28字；页码数字置于下方正中），并装订成册。
如果仅美欧日三个发达地区，地球物理学领域的代表性期刊应该是美国的《Journal of Geophysical Research》和《Geophysical Research Letters》，欧洲（德国）的《Annales Geophysicae》和日本的《Earth, Planets and Space》等。
地球物理学报、地球物理学进展、地震学报、中国地震、地震地质、地震、地震研究、石油地球物理勘探、物探与化探、石油物探、测井技术、石油勘探与开发、地震工程与工程振动、空间科学学报
地球物理方面：1.地球物理学报 2.地震学报 3.地震地质4.地震工程与工程振动 5.地震 6.中国地震 7.地震研究 8.地球物理学进展9.西北地震学报10.水文石油、天燃气方面：1.石油勘探与开发 2.石油学报 5.天然气工业 3.石油与天然气地质 6.石油化工 4.石油实验地质7.石油物探 8.中国石油大学学报.自然科学版 9.天然气地球科学 10.西南石油大学学报.自然科学版 11.石油钻采工艺 12.新疆石油地质 13.测井技术14. 油气地质与采收率 15.大庆石油地质与开发 16.钻采工艺 17.油田化学 18.石油钻探技术19.石油炼制与化工 20.石油地球物理勘探 21.特种油气藏 22.石油机械 23.西安石油大学学报.自然科学版 24.钻井液与完井液 25.石油学报.石油加工 26.大庆石油学院学报 27.油气田地面工程 28.海相油气地质 29.中国海上油气
中文核心，和CSSCI期刊都是重要期刊。
1908年（清光绪三十四年）6月25日，顾功叙出生于浙江省嘉善县。父亲是农村小学教员，幼年家境清寒，得到亲戚的资助，才有机会进中学和大学学习。1929年（民国十八年），在上海大同大学毕业后，到杭州浙江大学物理系任助教。1933年（民国二十二年），清华大学向全国招收公费留学生，顾功叙通过考试被录取，指定去美国学习地球物理勘探。从1933年底到1934年7月，顾功叙在北平清华大学为出国留学做准备，学习了与地球物理有关的基础课程，导师是翁文灏、袁复礼和叶企孙三位教授。1934年（民国二十三年）8月，顾功叙从上海乘海轮远涉重洋前去美国，进了科罗拉多州矿业学院，读地球物理勘探研究生。1936年（民国二十五年），通过论文答辩，取得硕士学位。他是中国第一位系统地掌握地球物理勘探学理论、方法和技术的学者。毕业后去加利福尼亚加州理工学院，在著名的地球物理和地震学家古登堡教授指导下做研究工作。1938年（民国二十七年），顾功叙中断在美国的研究工作回国，到抗战后方的云南省昆明，任当时已搬迁到昆明的北平研究院物理研究所研究员。那时科学研究和生活条件都很差，没有必要的仪器设备，他仅凭几台简陋的仪器在云南、贵州两省的10 余个地区，开展了中国较早的地球物理勘探工作。1947年（民国三十六年）2月，北平研究院物理研究所迁回北平，他与翁文波、傅承义、陈宗器、赵九章、王之卓等地球科学家们一起，发起和创建了中国地球物理学会，创刊编辑出版《地球物理学报》。顾功叙历任该学会常务理事、理事长、名誉理事长和《地球物理学报》主任委员、主编、名誉主编。中华人民共和国成立后，顾功叙任中国科学院地球物理研究所研究员兼副所长，主持地球物理勘探的研究工作。1949年下半年，他接受当时中央重工业部的委托，在北京举办了新中国最早的地球物理勘探训练班，从无到有培养了第一批地球物理勘探人员。1950年加入九三学社。1951年，国家地质工作计划指导委员会成立，顾功叙被任命为委员和地质矿产勘探局地球物理勘探处处长。1952年，中华人民共和国地质部成立后，顾功叙先后任该部地质矿产司副司长、地球物理勘探局副局长、总工程师。具体组织和指导地质部系统的地球物理勘探工作。1953年，顾功叙参加北京地质学院地球物理勘探系的建立，并在该系开办了全国各部门物探技术人员大规模的冬训班，不但进一步提高了中国物探人员的技术水平，而且为全国各部门物探队伍之间的沟通，为生产、科研和教学部门之间的横向联系打下了基础。地球物理勘探在中国是一门新兴的科学技术，还处于发展阶段，在进行实际勘探工作的同时，科学研究十分必要，顾功叙不仅在中国科学院地球物理研究所组织并指导这方面的科学研究，还把地质部地球物理勘探局的一个试验室发展组建成地球物理勘探研究所，1959年起任该所所长。1966年，河北省邢台大地震之后，遵从周恩来总理要求地质和地球物理工作者加强地震预测预防研究的指示，顾功叙就把主要精力转移到地震预报研究方面；主管中国科学院地球物理研究所的地震科学研究工作。1977年，被选为国际大地测量和地球物理联合会（IUGG）中国委员会主席。国家地震局成立后，1978年顾功叙任该局地球物理研究所副所长、研究员，继续指导深入开展地震预报及有关问题的研究工作。1979年顾功叙与傅承义、翁文波、马杏垣、张文佑、李善邦、秦馨菱、谢毓寿、曾融生等地球科学家们发起和创建了中国地震学会，创办编辑出版《地震学报》。顾功叙历任该学会理事长、名誉理事长和《地震学报》主编、名誉主编。1985年，顾功叙退出领导工作的第一线，任名誉所长，继续进行研究工作和培养指导研究生。 1992年1月14日，逝世于北京，享年84岁。
傅承义是一位爱国的、正直的科学家、教育家。不但具有卓越的才能，而且更有为人们所敬重的高尚品德。为人光明磊落，从不迎合潮流，随声附和，人云亦云。在是非问题上，直言不讳，刚正无私，傅承义认为是正确的，就敢于坚持。这曾使他吃了不少苦头，蒙受了不少冤屈。1958年大跃进时，有些人只从良好的愿望出发，不顾事物发展的客观规律，在地震预测工作中提出一些不切实际的口号，并形成一种声势，使得这一工作偏离了正确方向。傅承义以其科学家的责任感，反对这种做法，结果被扣上“反对搞地震预报”的帽子，直到“文化大革命”中还屡遭批判。1975年，海城地震之后，有人过高估计海城地震“预测成功”的经验，认为在中国解决地震预测问题已近在眼前。针对这种倾向，傅承义一针见血地指出，海城地震的预测是“歪打正着”，提醒人们不能头脑发热。20世纪60年代初，傅承义担任七室主任期间，亦曾因在工作中坚持正确的学术观点而被指责为“反对搞国防任务”，受到不公正的待遇。对于诸如此类的事情，傅承义都泰然处之。作为一个科学家，他把按照科学规律办事看成是自己的天职，坚持真理是他人生追求的最崇高目标，至于个人的荣辱得失，看得淡如清水，从不计较。傅承义在学术界有很高声望，但傅承义对别人的恭维却极为反感，也从不以自己的声望谋取私利。傅承义非常注重科学道德。在研究工作中，受过他的指导帮助的人很多，但是，不管傅承义对别人的研究工作出过多大力，傅承义从不在研究成果上署名，包括他指导完成的研究生论文。傅承义思路敏捷，在学术讨论中，直言快语，不讲情面，但从不以势压人。傅承义治学态度十分严谨。傅承义的论著，字字句句都经过仔细推敲，不仅内容深刻、丰富，而且文章结构严谨，条理分明，逻辑性强，文笔生动、流畅。傅承义的严谨学风也反映在他长期担任《地球物理学报》主编的工作中。凡是投到学报的稿件，一视同仁，傅承义都亲自审定，录用、退稿一定要有他的签字。对于有争议的稿件，处理更为慎重。英文版的每篇文章，他都要亲自把关。在傅承义的指导和带动下，《地球物理学报》在国内外赢得了广泛赞誉，成为中国被世界四大检索系统同时选用的九种刊物之一，连续被评为中国科学院优秀期刊。傅承义以其高尚的品德和卓越的科研、教学实践，为树立了一代人民科学家、教育家的风范。2000年1月8日，傅承仪与世长辞。刘光鼎院士挽之：“三篇文章开世界震波研究先河创新典范，一生耕耘育中华找矿精英大成风节长存”。
傅承义祖籍福建闽侯。1909年10月7日生于福建闽侯。祖父在清朝做过道台，伯父在北洋军阀时期的海军部供职。父亲傅仰贤长期在北洋军阀政府外交部及驻外使馆工作，曾任驻前苏联列宁格勒总领事。父亲虽然是旧官吏，但有强烈的爱国思想，亦比较开明。傅承义兄弟姊妹4人从小就受到良好的家庭教育，对其后来的长产生重要影响。哥哥傅鹰是著名化学家，担任过北京大学副校长。傅承义自幼记忆力极好，家里专门请了私塾先生教他和姐姐读《四书》、《五经》，聘请家庭教师教授数学和英语。但他用更多的时间博览群书，他童年虽然没有进过小学，但知识和能力远远超出了同龄小学生。1923年14岁时，在母亲提议下，傅承义跨入北京育英中学校门。一年之后，他感到功课太容易，便背着家里跳两级报考了汇文高级中学，结果考取了，但没有去读。初中三年，学习成绩年年名列全校第一，数学和英语成绩尤为突出。他多次参加学校组织的国语和英语讲演比赛，总是名列前茅。《福尔摩斯探案全集》是他最喜欢阅读的英文原著之一。他更喜欢逻辑推理，并勤于演算。傅承义不但学习成绩出类拔萃，还曾赢得全校三跳（跳高、跳远、三级跳远）及百米跑第一名。1926年，傅承义考入汇文高级中学。1929年，傅承义以获得理化、数学两项银杯奖的优异成绩结束中学时代的生活。傅承义原已考上燕京大学，但更喜欢清华大学的校风，选择了清华大学物理系。在这里，他的旺盛求知欲亦无法得到满足，更不满意当时教师队伍中的某些不正风气，痛切感受到教书育人者为人师表之重要性。傅承义和同窗好友王竹溪暗下决心，有朝一日教书育人，一定要立德、立言、立身。大学四年，傅承义以自学为主，从教师之讲解获益不多。虽考试成绩仍能保持在中、上之间，却因不重视教师的启发，多走了许多弯路，事倍功半，浪费了许多时间，事后检查，追悔不已。1933年，大学毕业后留校，傅承义先做一年研究生，后因教学工作需要，当上助教，从事核物理教学实验和研究工作。在此期间，他分别与黄子卿、赵忠尧合作，完成有关热力学研究和核物理实验方面的论文共4篇。1937年，抗日战争爆发，清华大学举校南迁，在昆明与北京大学、南开大学合并成立西南联合大学。1938年，傅承义应邀到西南联合大学继续任教。1939年，傅承义考取英“庚款”公费留学，在当时是全国少数优异学生才享有的殊荣，而地球物理专业仅此一个名额。由于第二次世界大战爆发，直到1940年，傅承义与林家翘、郭永怀、钱伟长等一行24人才转赴加拿大。傅承义进入麦吉尔大学物理系，师从当时最有声望的地球物理探矿学权威D·A·基斯（Keys）教授，进修地球物理勘探。1941年傅承义获得硕士学位。基斯教授对傅承义的成绩极为赞许，推荐傅承义到地球物理勘探领域里颇负盛名的美国科罗拉多矿冶学院继续攻读博士学位。可是傅承义对这里的专业方向不甚满意。此时，又赶上腰病发作，医生建议傅承义不要做野外工作，傅承义放弃了地球物理探矿专业。1942年，基斯教授又把傅承义推荐到加利福尼亚理工学院研究生院，师从近代地球物理学泰斗B·古登堡（Gutenberg）教授，攻读地球物理学及地震学。古登堡对傅承义在学习中和在学术论坛上表现出来的才能极为赞赏，将自己没能解决的一个理论问题——从理论上证明沿分界面传播的所谓“折射”地震波的存在，让傅承义去解决。傅承义凭借着深厚的物理学和数学基础，从数学上严密地论证了首波的存在，并从物理学上解释了首波与折射地震波之间的区别。此项研究成果得到古登堡高度评价，他也因此受到广大师生的推崇。1944年获该校地球物理学博士学位。随后受聘于几家石油、地球物理勘探公司，做技术咨询工作。傅承义以对地球物理学发展所作的贡献，赢得地球物理学界的普遍承认，1946年被聘为加利福尼亚理工学院地球物理学助理教授。在此期间，傅承义在地震波传播的研究领域里，发表了一系列具有创造性和开拓性的研究成果，成为地震波研究的先驱。傅承义发表在美国《地球物理》杂志上一组论文，系统地研究了地震体波、面波及首波的传播等问题。这些论著无论是在中国还是在美国、前苏联等国家都引起极大的重视。在1960年纪念该杂志创刊25周年之际，这组论文被评为地球物理学经典著作。1947年春天，傅承义收到大学时同窗好友、中央研究院气象研究所所长赵九章的来信，希望傅承义能回国主持气象研究所的地球物理研究工作。傅承义毫不犹豫，两周之后便启程回国，到气象研究所任高级研究员，并兼任中央大学物理系教授。1948年，国民党当局责令中央研究院气象研究所和历史语言研究所迁往台湾，傅承义与赵九章、陈宗器一起予以抵制，为新中国地球物理事业的发展保存了力量。1950年4月，中国科学院地球物理研究所成立，傅承义任研究员。1952年，国家决定从大学物理系抽调一批优秀毕业生，从事地球物理探矿工作，由傅承义主持对他们进行培训。1953年，中国科学院接受北京地质学院的请求，委托傅承义去该院任地球物理探矿教研室主任。当时物探教研室初建，傅承义面临的任务十分繁重，傅承义不仅要向大学生讲授“地球物理勘探”课，而且还要给教师（全部是物理系毕业生）系统讲课。为了使教师能尽早走上讲台，傅承义夜以继日地工作，为每一位教师修改、审定讲稿；为了使教师在讲台上能站得住、讲得好，傅承义还亲自去听课并作讲授示范。一次，一位实验员在准备磁法实验时，失手将刃口式磁秤掉到了地上，刃口出现了一个缺口，傅承义在得知这件事后，立即将教研室全体人员召集到实验室。傅承义首先指出，这不是件小事，而是个错误，特别是发生在教学实验室里。学校要培养学生爱护仪器，杜绝任何操作中的失误，否则就不能保证野外观测的质量。现在实验室里发生了这种事，对学生会有什么影响呢？接着，傅承义又作了自我检查，说自己对实验室工作抓得不严，没有像过去在清华大学工作时那样认真对待实验。最后傅承义表示，希望大家都以此为鉴，上下一心搞好教学工作，不再犯类似错误。傅承义处理这件事，实际上是对教研室全体人员的一次极其生动而又非常深刻的思想教育。傅承义在北京地质学院创建中国第一个地球物理教研室的3年（1953—1956）时间里，和教研室的同事们朝夕相处，言传身教，使每个人都深受其益。傅承义常告诫年轻的同事们，作为一个地球物理学家，既要有理论修养，又要能够动手实践，在实践中积累经验，再上升到理论高度上去把握这些源于实践的经验。而在实践中则必须学会根据地质条件去部署工作，正确地进行观测，并对所得资料作出符合客观地质情况的解释，等等。这些方法的传授同知识的传授一样给人以教益。然而比知识和方法的传授更重要的是，傅承义那严谨的治学态度、认真的学术作风和对己、对人的严格要求使大家懂得了作为一名科学家所应该具有的最根本的品质，这使得年轻的教师们在前进道路上少走了弯路。此后，傅承义教授于1956—1961年间在北京大学创建地球物理教研室，于1964—1966年间在中国科学技术大学创建地球物理教研室，分别主持领导这些教研室的工作，并担任第1任教研室主任。1973年兼任中国科学技术大学地球及空间科学系主任。在他的教学生涯中，他始终把高尚的科学道德、严谨的治学态度和献身科学事业的精神贯彻到教学工作中，深受广大师生的推崇和爱戴。傅承义在地球物理教育战线上辛勤耕耘30余载，在中国地球物理学界真可谓桃李满天下。傅承义在为发展中国地球物理教育事业的同时，潜心学习，使自己的学识水平处于地球物理学科的发展前沿。1956年，傅承义回国后的部分论著《地震面波的能量束》、《关于瑞雷波方程的无关根》、《平行介质中的弹性波之传播》、《地下薄地层自由振动》、《折射探矿法的研究》和《地表层的本质对于地震勘测的几种影响》等6篇文章，以“关于弹性波的传播理论和地震探矿的一些问题”项目，荣获国家自然科学三等奖。同年，傅承义参加中国12年科学技术发展远景规划制订工作，是第33项任务“中国地震活动性及其灾害防御研究”的两执笔人之一。傅承义率先提出在中国开展地震预报研究的长远规划，并指出解决这一问题的科学途径及实施方法。这项工作领先其他先进国家约5—10年时间。为了开展核爆炸地震侦察研究，并藉此全面提高地震学发展水平。1961年，地球物理研究所成立第七研究室，由傅承义担任室主任。该室在核爆炸地震观测和地震侦察工作中，为国家作出重要贡献，并对中国地震学与测震学的发展起了推动作用。1971年，傅承义提出地震成因的“红肿假说”。1972年，傅承义创建震源物理研究室，并领导震源物理研究工作。从此，中国的震源物理研究工作上升到有组织、有计划发展的新阶段。这一年傅承义发表专著《大陆漂移，海底扩张和板块构造》，把20世纪地球科学的最新理论成就——板块大地构造假说介绍到中国，为中国地球科学的发展指明了方向。1976年十年动乱结束之后，中国迎来了科学的春天，傅承义把主要精力投入到科学著述和研究生培养工作中。编著了《地球十讲》、《地球物理学基础》；主编了《中国大百科全书·固体地球物理学》，并亲自撰写其中的部分条目。傅承义是中国固体地球物理科学的主要奠基人和开拓者之一。1957年，傅承义当选为中国科学院学部委员（院士）。傅承义曾任中国科学院地学部常务委员，中国科学院地球物理研究所一级研究员、室主任、所负责人、所学术委员会主任、名誉所长。傅承义是中国地球物理学会（1947）的发起人之一，长期担任《地球物理学报》主编。曾任中国地球物理学会副理事长兼秘书长和中国地震学会副理事长，中国地质学会和中国声学学会理事，中国地球物理学会和中国石油物探学会名誉理事长，中国地震学会名誉理事，《中国大百科全书》总编辑委员会委员、固体地球物理编辑委员会主任，全国自然科学名词审定委员会委员、地球物理学名词审定委员会主任。曾被选为第三届全国人大代表，担任过第二、五、六届全国政协委员。1981年加入中国共产党。2000年1月逝世于北京。