发布时间：
你是文科还是理科的呢 我也想学地理学 但是不晓得以后能做什么 请赐教 374228520我的QQ
中国的水资源
上地理学报编辑部的网站，上面有详细的说明和要求，可以下载pdf来看。
（1）钟华平.城市化对水资源的影响[J].世界地质,1996（2）“十一五”期间区域经济发展的资源约束与环境承载力问题研究[J].领导决策信息，2004（3）左其亭.城市水资源承载能力——理论、方法、应用[M].北京：化学工业出版社，2005（4）《造福子孙后代》（5） 行政院环保署「安全饮用水手册」（6） 中国水资源公报 1999 中华人民共和国水利部（7） 中国城市排水建设与发展 （8） 全国人大常委会法制讲座第二十六讲：关于我国农业法制建设的几个问题（9） 水资源合理分配的基本路线 作者：关业祥 水利部规计司处长（10） 张岳，刘钰等．加快建立我国农业节水保障体系的对策和建议．中国水利，（11）刘 跃.中国水资源与可持续发展.《中国水资源与可持续发展》（12）楚泽涵，水资源问题应引起关注 《古地理学报》（13）李圭白，李星，水的良性循环与城市水资源 《中国工程科学》（14）王水燕 谈水资源的可持续发展, 科技资讯
苏正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:.
朱生旺1，2曲寿利1魏修成1刘春园3
（1.中国石化石油勘探开发研究院，北京100038；2.中国地质大学（北京），北京100083；3.中国石油大学（北京），北京102249）
摘要 研究复杂介质中地震波传播规律及地震响应特征需要在细小网格剖分下进行弹性波方程数值模拟计算，而小网格下的数值模拟计算将带来巨大的计算量问题，采用变网格计算是减少计算量的有效途径。本文给出一种变网格差分计算的实现方法，在局部复杂介质区域采用细网格计算，其余区域采用粗网格计算，在两种网格的过渡区通过改变差分算子和波场插值实现波传播的过渡衔接。该实现方法的特点是：①在所有计算节点均采用交错网格计算方式；②估算空间一阶导数时，保持差分计算节点的对称性；③可选择任意阶计算精度。该方法的计算结果表明，粗细网格的过渡不会对地震波传播模拟带来影响，从而达到了既减少计算量又保证计算精度的目的。
关键词 弹性波方程 数值模拟 变网格差分 频散
Elastic Modeling on A Staggered Grid with Varying Spacing
ZHU Sheng-wang1，2，QU Shou-li1，WEI Xiu-cheng1，LIU Chun-yuan3
（ ＆ Production Research lnstitute，SlNOPEC，Beijing100083； University of Geoscience，Beijing100083； University of Petroleum，Beijing102249）
Abstract It is necessary to use small grid in elastic wave equation modeling for researching the law of seismic wave propagation and feature of seismic response in complex is an effective scheme to calculate on a staggered grid with varying spacing which can reduce the enormous calculation amount caused by using the small grid calculation article proposes the finite difference method in which small grid is used for local complex medium，and big grid is used for transition region between two kinds of grid，transient conjugation of wave propagation is implemented by changing the difference operator and applying wave field trait of this way includes：（1）Keeping staggered-grid calculation mode in all operation node；（2）Maintaining The symmetry of operation node in estimating spatial first derivate；（3）Any precision order available in finite difference result of the way indicates that the simulation of seismic wave propagation will not be affected by using two kinds of grid，and as a result，not only the calculation amount decreases，but also the high calculation accuracy remains.
Key words elastic wave equation modeling staggered-grid dispersion
在复杂油气藏的勘探开发中，地震数值模拟（正演）发挥着越来越重要的作用。利用地震方法进行储层预测，首先要研究储层的地震响应特征，找出与储层最为相关的地震属性，进而优选储层预测的技术方法，地震数值模拟是研究储层地震响应特征不可缺少的手段。
对复杂储层或地质体进行地震数值模拟需在足够小的空间网格上进行计算，以保证地质模型的局部变化细节能够准确地反映到数值模拟结果中来，即确保数值模拟结果能够较为精确地反映介质的小尺度非均质性引起的波场变化。在正演计算中，为保证计算收敛，空间和时间方向的采样率必须满足收敛条件，由弹性波方程有限差分正演的收敛条件［1，2］知，减小空间计算网格尺度的同时，时间方向的递推计算步长也必须大体上做相应比例的减小，即空间计算网格点数增加 n倍，意味着计算量将增加 n2倍。因此，若对整个模型采用细小网格计算所面临的问题是计算量太大，即使是二维模型的正演计算，譬如使用1m×1m 或更小尺度的计算网格，当采用能够反映实际地震反射深度的地质模型时，模拟多次覆盖采集，其计算量亦十分惊人。实际中的数值模拟一般是针对某一特定的地质目标（或储层），且地质目标仅占据整个模型很小的一部分，对计算网格仅要求在该地质目标处细化，于是，采用变网格计算，即在目标地质体区域采用较小的计算网格，而在该区域以外采用较大的计算网格就可以大大地减少计算网格节点的数量，从而可在很大程度上降低计算量。
为减少计算量，提高数值模拟的计算效率，Jastram C和Tessmer E提出了单方向上变网格计算方法［3］，文献［2］亦强调采用变网格计算技术，但没有给出具体实现方法。采用变网格进行差分计算的关键是解决好网格尺度变化处的波场过渡衔接问题，保证在过渡区不会因网格尺度变化产生明显的计算噪声。本文基于二维弹性波方程数值模拟的交错网格差分算法，通过引入不等间距节点情况下一阶导数的高阶精度差分计算式，结合波场的高阶精度插值和粗细网格过渡区的计算节点选取技巧，给出弹性波方程变网格正演计算的实现方法，应用结果表明，该方法能够得到满意的计算精度，在过渡区不会产生明显的因网格尺度变化带来的附加计算噪声。
1 利用对称任意间距节点估算一阶导数的高阶精度差分近似式
在沿时间方向的递推计算过程中，提高波场函数空间一阶导数的估计精度是弹性波方程数值模拟研究的关键内容之一。对于固定差分网格，通过应用等距节点高阶精度的差分算子可实现空间导数的高精度估算［1，2］。对于局部采用细网格的非固定差分网格情况，在粗、细网格区域内部可分别应用上述等距节点的高阶精度的差分算子，而在粗、细网格过渡区，本文采用下面导出的具有任意节点间距的差分算子估算一阶导数值。
设在x方向上以点x0为中心对称分布2N个网格节点，其x坐标分别为x0-qNΔx/2，…，x0-q1Δx/2，x0+q1Δx/2，…，x0+qNΔx/2，如图1 所示，这里Δx为节点间的最小间距。由2N个节点处的函数值f（x0-qNΔx/2），…，f（x0-q1Δx/2），f（x0+q1Δx/2），…，f（x0+qNΔx/2）可给出函数f（x）在点x0处的一阶导数估值。由泰勒展开有：
图1 估算一阶导数的对称计算节点示意图
油气成藏理论与勘探开发技术
油气成藏理论与勘探开发技术
两式相减，得
油气成藏理论与勘探开发技术
略去（1）式中误差项O（（qiΔx/2）2N+1），将精确的f（i）（x0）换成f（i）（x0）的估计 （x0），并写成矩阵形式，则有
油气成藏理论与勘探开发技术
记
油气成藏理论与勘探开发技术
油气成藏理论与勘探开发技术
设A-1为A 的逆，即 AA-1=I，则（A-1）TAT=I，即（A-1）T为 AT的逆，从而有 AT（A-1）T=I，右乘向量e1得
AT（A-1）Te1=e1 （3）
（3）由方程（2）知
油气成藏理论与勘探开发技术
记向量（A-1）Te1=（c1，c2，…，cN）T，代入式（4）则有
油气成藏理论与勘探开发技术
且由式（3）及（A-1）Te1=（c1，c2，…，cN）T知cn（n=1，2，…，N）满足方程
油气成藏理论与勘探开发技术
式（5）即为用对称任意间距节点计算一阶导数的计算式，cn（n=1，2，…N）通过求解方程（6）得到。若取qn=2n-1（n=1，2，…，N），则cn（n=1，2，…，N）即为等距节点情况的一阶导数具有2N阶精度差分近似的差分系数。
将f（x）做傅氏变换变换到波数域，即f（x） F（k），对f（x）一阶导数进行傅氏变换则有f（1）（x） .i2πkF（k），即f（x）可看作是对f（x）进行的空间线性滤（1）波结果，其滤波响应为i2πk。用f（x）离散点的值估算∂f/∂x得 （x），亦看作是对f（x）进行的空间线性滤波结果，记其滤波响应为 （k），由式（5）知 （k）= cnsin（πkqnΔx）。记
油气成藏理论与勘探开发技术
显然，用 （x）作为f（1）（x）的估计，其精度可由 （k）逼近i2πk的程度，即e（k）来反映。
2 一阶导数的变网格有限差分估算
采用变网格计算，以x方向的一阶导数估算为例来说明两种网格过渡处的处理方法。设细网格节点间距为Δx，粗网格节点间距为mΔx，即粗网格节点间距是细网格节点间距的m倍。为表述方便，不妨取m=3，采用10阶精度，即N=5，如图2所示。图中 表示粗网格节点； 表示粗网格中的加密节点，其与粗网格节点一起组成细网格节点；“+”为求导点位置。①～⑥为粗细网格过渡区的6个一阶导数计算点，计算该6个点处的一阶导数所用的网格节点分别在图2中标示，其中①及其左边大网格区的各计算点按粗网格计算；⑥及其右边的细网格区各计算点按细网格计算；②～⑤4个点需按式（5）计算。
图2 两种网格过渡区计算示意图
一般取2N阶差分计算精度时，在粗、细网格分界的细网格一侧共有N-1个计算点需要采用如式（5）给出的对称不等间距节点差分算子进行一阶导数估算，每一个点的差分算子不同，加上粗网格和细网格内部计算点的差分算子，这样共有N+1个不同的差分算子。N+1个不同的差分算子可事先算好存于内存数组中供每次计算调用。
图3 10阶精度（N=5）不同计算点差分算子对应的滤波响应
①～⑥分别为图2中对应计算点的差分算子的滤波响应；⑦为理想的滤波响应
图3 给出的是图2中的①～⑥求导计算点的估算一阶导数差分算子的滤波响应。对于计算点⑥，由于参加计算的节点为细网格中的相邻2N个网格节点，显然其差分算子的滤波响应与理想的滤波响应最接近，精度最高；对于计算点①，参加计算的节点为粗网格中的相邻2N个网格节点，其差分算子的滤波响应仅在低波数与理想的滤波响应接近，精度最低；而对于计算点②～⑤，精度介于计算点①与⑥之间，且随着参加计算的近距离点增加，其精度越来越高。
在粗细网格过渡区一阶导数的估计精度不一致可能会引起反射，附加计算噪声，现对此做出分析。
用有限差分方法估算空间一阶导数的精度由差分频散决定，有限差分算子的滤波响应只能在低波数段逼近一阶导数算子的滤波响应，在高波数段将出现严重的频散。设kα为具有可接受差分频散水平的最高波数，即在k≤kα的低波数区， （k）与i2πk足够逼近，频散很小，对最终正演结果的影响可忽略不计。记kα=αkN，这里kN=1/（2Δx），0＜α＜1，α值由差分算子决定，如：对均匀网格时的10 阶精度的差分算子，可认为α约为。又设模拟地震波场的最大波数为kmax，它由子波最高频率fmax与最低速度vmin决定，即kmax=fmax/vmin。为了使有限差分数值模拟计算结果不受差分频散的显著影响，必须选取Δx满足α/（2Δx）=kα＞kmax=fmax/vmin，即
油气成藏理论与勘探开发技术
式（7）是保证数值模拟计算精度的空间离散步长Δx应满足的条件，只要满足条件（7），则在波数范围［0，kmax］内差分频散误差可忽略。
采用变网格计算，粗网格尺度也应该满足式（7），即
油气成藏理论与勘探开发技术
否则波场在粗网格区传播就会产生强的频散噪声。式（8）成立，则条件（7）必然满足，由前面对两种网格过渡区的差分算子精度的分析结论知，此时过渡区的差分算子均满足精度要求，即kα＞kmax。既然所有差分算子都能够满足不产生明显的差分频散噪声的要求，差分算子变化也将不会引起明显的计算误差。
因此，对于两种网格过渡区波场传播的衔接问题，上述分析的结论是只要粗网格计算满足精度要求，即不产生明显的差分频散，则按本文的计算方法，在两种网格的过渡区不会出现明显的反射噪声，波场传播的模拟精度基本不受计算网格变化的影响。
3 弹性波方程变网格有限差分正演模拟的实现
实现方法
油气成藏理论与勘探开发技术
以二维弹性波方程数值模拟计算来说明变网格有限差分法的实现。利用二维弹性波方程进行有限差分正演计算，在交错网格中的递推算式为
油气成藏理论与勘探开发技术
上式中， 和 分别表示 （x，z，t） 和 ∂z（x，z，t） 的有限差分估计，余类推。在式（10）中，ρ，c11，c13，c33及c55都是空变的，为书写简洁，略去了其空间下标。计算中对边界的处理是，在顶边界采用自由边界条件，在其他数值边界则采用PML吸收边界条件［4，6］。
图4 粗、细两种网格中计算节点的布排
在一定的网格尺度下，交错网格计算具有精度高的优点。当采用变网格，即在模型局部用细网格剖分进行计算时，关键是要合理布置计算节点，以保证在所有计算节点均能采用交错网格计算方式。图4 以粗细网格尺度比等于3 为例，给出本文采用的计算节点布排方式，不难看出，要保持交错网格的计算方式，x和z方向粗细网格尺度比 mx和 mz都必须为奇数。
对波场递推计算中空间一阶导数的有限差分估算，根据节点所处位置，分为 3种情况分别处理（为简化图示，不妨设差分算子长度N=2，即4阶精度）：
（1）对细网格区域之外（图4实线框之外）的计算节点采用粗网格进行计算；
（2）对细网格区域内部（图4虚线框内）的计算节点采用细网格进行计算；
（3）对细网格边界区域（图4虚线框与实线框之间）的计算节点按图2描述的方法进行计算，但在计算前对参与计算的某些点（图4 中空心点所在的位置）的波场值需要通过内插得到，本文采用拉格朗日插值方法，内插阶数与差分算子的阶数相同。
算例
首先，用一个均匀介质模型检验网格变化是否对模拟波场传播产生明显的影响。
如图 5 所示，模型大小为 1000m×1000m，vP=3000m/s，vS=1800m/s，ρ=，vz震源位于模型正中位置，采用两种网格划分进行正演计算。第一种网格划分如图5（a）所示，在225～275m深度段（图中阴影部分）用1m×1m的细网格，其余部分用5m×5m的粗网格，采用该变网格进行正演计算，图5（b）和（c）分别是160ms时刻的x分量和z分量波场切片；第二种网格划分如图5（d）所示，在725～775m深度段（图中阴影部分）用1m×1m的细网格，其余部分用5m×5m的粗网格，采用该变网格进行正演计算，图5（e）和（f）分别是160ms时刻的x分量和z分量波场切片。比较两组计算结果可以看到，图（b）与（e），（c）与（f）几乎没有差别，即局部采用细网格对波场传播没有带来明显的影响。为了看得更清楚，将沿过震源的垂直线上（图5（d）中虚线所示）的波传播情况显示于图6，从向上、向下两个方向传播的波的一致性可知，波通过细网格区域时没有产生反射噪声。
图5 变网格计算与波场切片
其次，通过一个含不同尺度圆形洞的介质模型正演结果观察变网格差分计算的效果。
如图7所示，模型大小为1500m×600m，vP=5000m/s，vS=3000m/s，ρ=，在模型中部500m深度上分布3个半径分别为5m，10m，20m的3个圆形洞，洞间相距200m，洞内vP=1800m/s，ρ=。模拟野外地震观测，炮间距10m，道间距5m，排列长度为1500m，中间激发，激发震源处于模型中部地表，激发子波主频为40Hz，采用变网格计算，粗网格尺度为5m×5m，细网格（图7上的白色虚线框内）尺度为1m×1m，共得150个炮记录，抽出零偏移距剖面，其结果如图8所示。图9（a）是激发点位于模型中部的一炮记录；图9（b）是与图9（a）相对应的对整个模型用1m×1m细网格计算的炮记录。将这两个炮记录进行对比可看到二者差别甚微，证实采用变网格计算基本没有降低计算精度。但采用变网格计算却极大地减少了计算量，就该模型来说，如果对整个模型采用1m×1m的网格进行计算，其计算量约为采用上述变网格计算的计算量的15倍。不难看出，如果对更大尺度的模型进行正演计算，采用变网格计算对减少计算量就显得更有必要。
图6 变网格对波场传播影响测试
图7 含圆洞的介质模型
4 结论
在对复杂介质进行数值模拟计算时，采用变网格计算是解决计算精度与计算量之间矛盾的有效途径。理论分析和数值模拟结果表明，本文给出的适应变网格计算的弹性波方程正演模拟实现方法保持了交错网格高阶差分的精度高特点，在粗细网格过渡区不产生人为反射噪声，是一种在保证计算精度前提下减少计算量的有效实用方法。
图8 圆洞模型的零偏移距剖面
图9 变网格与固定小网格计算的炮记录对比
参考文献
［1］董良国.一阶弹性波方程交错网格高阶差分解.地球物理学报，2000，43（3）：411～419.
［2］牟永光，裴正林.三维复杂介质地震数值模拟.北京：石油工业出版社，2005.
［3］Jastram C，Tessmer modeling on a grid with vertically varying .，1994，42：357～370.
［4］Francis Collino，Chrysoula of the perfectly matched absorbing layer model to the linear elastodynamic problem in an isotopic heterogeneous ，2001，66（1）：294～307
［5］Cunha C modeling in discontinuous ，1993，58（12）：1840～1851.
［6］Rune boundary conditions for elastic staggered-grid modeling ，2002，67（5）：1616～1623.
参考文献是文章或著作等写作过程中参考过的文献。
因参考文献的著录格式各刊不尽相同，投稿前作者应注意杂志稿约的有关规定，至少得先看看有关期刊发表的论文的参考文献是如何标注的，以了解有关期刊的参考文献的著录格式，以免出错。许多作者投递的稿件书写格式包括参考文献的著录格式与杂志所要求的不同。
坦率地讲，编辑和审稿专家也是人，工作中多少也有感情因素。如果拿到手中的是一篇书写格式不合要求的文章，别的暂且不论，就书写格式不规范这一条，就足以给编辑留下不好的印象，甚至让编辑做出退稿的决定。
就算最后没有被退稿，此类稿件较书写格式规范的稿件被录用的可能性大大降低。其实作者犯的是一个很低级的错误，让编辑很自然地联想到，该作者不太尊重期刊，还有期刊的编辑以及审稿专家。
因此，作者在投稿前一定要注意期刊参考文献的著录方式，以免产生不必要的负面影响。其实，并不复杂，只要稍稍留意即可。
地理教学论文参考文献
现如今，大家或多或少都会接触过论文吧，借助论文可以达到探讨问题进行学术研究的`目的。那要怎么写好论文呢？下面是我精心整理的地理教学论文参考文献，希望能够帮助到大家。
参考文献
[1]王建,张茂恒,白世彪.圈层相互作用与自然地理学[J].地理教育,2008(4):4-7
[2]张猛刚, 雷祥义.地球表层系统浅论[J].西北地质,2005,38(2):99-101
[3]马宗晋,杜品仁,卢苗安.地球的多圈层相互作用[J].地学前缘,2001,8(1):3-8
[4]狄丽颖,冯银珍,张爱国.地球系统科学研究进展[J].河南理工大学学报,2007,26(3):267-280
[5]袁道先.对地球系统科学的几点熟悉[J].高校地质学报,1999,5(1):1-6
[6]周俊.关于地球表层的讨论[J].地球科学进展,1993,8(1):72-77
[7]谢家泽.关于地球表层系统观的几个题目[J].地球科学进展,1995,10(5):432-435
[8]叶笃正,季劲钧,严中伟,等.简论人类圈(Anthroposphere)在地球系统中的作用[J].大气科学,2009,33(3):409-415
[9]孟凡刚,孙彦顺.地球的圈层结构[J].地理教育,2009 (4):11-11 [10]卡哈尔·尼牙孜,黄建华,布玛丽娅木·阿布拉.浅议地球系统科学与可持续发展[J].西部探矿工程,2010 (5):85-87
[11]陈志明, D. A.季莫菲耶夫,V. V.勃朗古里耶夫,等.试论地貌圈概念、属性及其与其它地表圈层关系[J].地球科学进展,
上地理学报编辑部的网站，上面有详细的说明和要求，可以下载pdf来看。
北京市科学技术委员会 .1997.可持续发展词语释义 ［M］ .北京: 学苑出版社
国家环境保护局自然保护司 .1995.中国生态环境补偿费的理论与实践 ［M］ .北京: 中国环境科学出版社
摆万奇 .2000.深圳市土地利用动态趋势分析 ［J］ .自然资源学报，15 (2): 112 ～116
摆万奇，赵士洞 .1997.土地利用和土地覆被变化研究模型综述 ［J］ .自然资源学报，12 (2): 169 ～175
毕晓丽，洪伟 .2001.生态环境综合评价方法的研究进展 ［J］ .农业系统科学与综合研究，17 (5): 122 ～125
卞建民，林年丰，汤洁 .2002.土地盐碱荒漠化预警理论及实证研究 ［J］ .农业环境科学学报 .22 (2): 207 ～209
蔡运龙 .2001.土地利用/土地覆被变化研究: 寻求新的综合途径 ［J］ .地理研究，20 (6): 645 ～652
常庆瑞，蒋平安，周勇等 .2004.遥感技术导论 ［M］ .北京: 科学出版社
畅军庆，张荣南 .1997.南水北调湖北库区水环境现状及问题分析 ［J］ .环境科学与技术，(2): 1 ～4
陈百明 .1997.试论中国土地利用和土地覆被变化及其人类驱动力研究 ［J］ .自然资源，(2): 31 ～36
陈国阶 .1996.对环境预警的探讨 ［J］ .重庆环境科学，18 (5): 1 ～4
陈国阶 .2002.论生态安全 ［J］ .重庆环境科学，24 (3): 1 ～4
陈国阶，何锦峰 .1999.生态环境预警的理论和方法探讨 ［J］ .重庆环境科学，21 (4): 8 ～11
陈美兰，吴延熊，周国模等 .1999.预测、监测和预警关系的初步探讨 ［J］ .浙江林学院学报，16 (1): 10 ～13
陈佑启，.中国土地利用-覆被的多尺度空间分布特征分析 ［J］ .地理科学，20 (3): 197 ～65
陈治谏，陈国阶 .1992.环境影响评价的预警系统研究 ［J］ .环境科学，13 (4): 20 ～26
丁一汇，王守荣主编 .2001.中国西北地区气候与生态环境概论 ［M］ .北京: 气象出版社
窦明，左其亭，胡彩虹 .2005.南水北调工程的生态环境影响评价研究 ［J］ .郑州大学学报 (工学版)，26 (2): 63～ 66
方创琳，张小雷 .2001.干旱区生态重建与经济可持续发展研究进展 ［J］ .生态学报，21 (7): 1163 ～1170
傅伯杰 .1991.区域生态环境预警的原理与方法 ［J］ .资源开发与保护，7 (3): 138 ～141
傅伯杰，陈利顶 .1996.景观多样性的类型及其生态意义 ［J］ .地理学报，51 (5): 454 ～461
傅伯杰，陈利顶，刘国华 .1999.中国生态区划的目的、任务及特点 ［J］ .生态学报，19 (5): 591 ～595
傅伯杰，陈利顶，马克明等 .2001.景观生态学原理及应用 ［M］ .北京: 科学出版社
高更和 .2002.南水北调中线工程源头区水土流失及其治理———南阳市汉水流域研究 ［J］ .水土保持通报，(3): 65 ～68
高志强，刘纪远，庄大方 .1999.基于遥感和 GIS 的中国土地利用/土地覆盖的现状研究 ［J］ .遥感学报，3 (2): 134～ 139
顾朝林 .1999.北京市土地利用/覆被变化机制研究 ［J］ .自然资源学报，14 (4): 307 ～312
关文彬，谢春华，李春平等 .2003.荒漠化危害预警原理与评价方法 ［J］ .北京林业大学学报 .25 (3): 79 ～83
关文彬，谢春华，孙保平等 .2001.荒漠化危害预警预景指标体系研究 ［J］ .北京林业大学学报 .23 (1): 44 ～47
郭中伟 .2001.建设国家生态安全预警系统与维护体系———面对严重的生态危机的对策 ［J］ .科技导报，(1): 54 ～56
胡安焱，郭生练，刘燕等 .2006.干旱区内陆河流域生态环境质量模糊综合评价 ［J］ .地质灾害与环境保护，17 (2):69 ～ 73
黄秉维 .1996.论地球系统科学与可持续发展战略科学基础 ［J］ .地理学，51 (4): 289 ～294
惠映河，蒋晓辉 .2001.水资源承载力评价指标体系研究 ［J］ .水土保持通报，21 (1): 30 ～34
蒋国富，白耕勤 .2004.南水北调中线工程水源区生态环境问题探析 ［J］ .南阳师范学院学报 (自然科学版)，3(9): 72 ～ 76
金岚主编 .1992.环境生态学 ［M］ .北京: 高等教育出版社
靳英华，赵东升，杨青山等 .2004.吉林省生态环境敏感性分区研究 ［J］ .东北师大学报自然科学版，36 (2): 68 ～74
康慕谊，江源，石瑞香 . 样带 1984 ～1996 年土地利用变化分析 ［J］ .地理科学，20 (2): 115 ～120
李红鹰 .2000.层次分析法在农业生态环境质量评价中的应用 ［J］ .监测与环境，(7): 30 ～31
李峻，董立克，慕金波等 .2000.区域环境质量预警研究 ［J］ .山东环境，(3): 19 ～20
李平，李秀彬，刘学军 .2001.我国现阶段土地利用变化驱动力的宏观分析 ［J］ .地理研究，20 (2): 129 ～138
李如忠 .2001.巢湖流域生态环境质量评价初步研究 ［J］ .合肥工业大学学报 (自然科学版)，24 (5): 987 ～990
李崧，邱微，赵庆良等 .2006.层次分析法应用于黑龙江省生态环境质量评价研究 ［J］ .环境科学，27 (5): 454 ～461
李晓兵 .1999.国际土地利用/土地覆盖变化的环境影响研究 ［J］ .地球科学进展，14 (4): 395 ～399
李晓秀 .1997.北京山区生态环境质量评价体系初探 ［J］ .自然资源，(5): 31 ～35
李秀彬 .1996.全球环境变化研究的核心领域 ［J］ .地理学报，51 (6): 553 ～556
李祚泳，邓兴民 .1995.综合环境质量的 B-P 网络二级评价 ［J］ .环境科学研究，8 (3): 32 ～35
李祚泳，彭荔红 .1999.基于人工神经网络的农业病虫害预测模型及其效果检验 ［J］ .生态学报，19 (5): 759 ～762
林年丰，汤洁 . 与环境模拟在环境地学研究中的作用和意义 ［J］ .土壤与环境，9 (4): 259 ～262
刘昌明 .2002.南水北调工程对生态环境的影响 ［J］ .海河水利，(1): 1 ～5
刘昌明，沈大军 .1997.南水北调工程的生态环境影响 ［J］ .大自然探索，16 (2): 1 ～6
刘彦随 .2001.山地土地类型的结构分析与优化利用 ［J］ .地理学报，56 (4): 426 ～436
陆登櫆 .1997.遥感技术在农业工程中应用 ［M］ .北京: 科技出版社
罗先香，邓伟 .2000.松嫩平原西部土壤盐渍化动态敏感性分析与预测 ［J］ .水土保持学报，14 (3): 36 ～40
马乃喜 .1998.区域生态环境评价中的几个理论问题 ［J］ .西北大学学报 (自然科学版)，28 (4): 330 ～334
毛汉英 .1996.山东省可持续发展指标体系初步研究 ［J］ .地理研究，15 (4): 16 ～23
梅宝玲，陈舜华 .2003.内蒙古生态环境预警指标体系研究 ［J］ .南京气象学院学报，26 (3): 384 ～393
苗鸿，王效科，欧阳志云 .2001.中国生态环境胁迫过程区划研究 ［J］ .生态学报，21 (1): 7 ～13
牛文元 .1994.持续发展导论 ［M］ .北京: 科学出版社，13 ～16
欧阳志云，王效科，苗鸿 .2000.中国生态环境敏感性及其区域差异规律初步研究 ［J］ .生态学报，20 (1): 9 ～12
冉圣宏，金建君，薛纪渝 .2002.生态脆弱区评价的理论与方法 ［J］ .自然资源学报，17 (1): 117 ～122
任志远，张艳芳 .2003.土地利用变化与生态安全评价 ［M］ .北京: 科学出版社
邵东国，李元红，王忠静等 .1996.基于神经网络的干旱内陆河流域生态环境预警方法研究 ［J］ .中国农村水利水电，(6): 10 ～ 12
盛连喜，冯江，王娓 .2005.环境生态学导论 ［M］ .北京: 高等教育出版社
史德明，梁音 .2002.我国脆弱生态环境的评估与保护 ［J］ .水土保持学报，16 (1): 6 ～10
史培军，陈晋，潘耀忠 .2000.深圳市土地利用变化机制分析 ［J］ .地理学报，55 (2): 151 ～160
史培军，宫鹏，李小兵等著 .2000.土地利用/覆盖变化研究方法与实践 ［M］ .科学出版社
史培军，江源，王静爱等 .2004.土地利用/覆盖变化与生态安全响应机制 ［M］ .北京: 科学出版社
苏维词，李久林 .1997.乌江流域生态环境预警评价初探 ［J］ .贵州科学，15 (3): 207 ～214
孙希华 .1999.山东农业生态环境质量评价与发展策略 ［J］ .山东师范大学学报 (自然科学版)，14 (2): 161 ～165
汤洁，林年丰，卞建民 .1999.农业生态地质环境质量综合评价———以吉林省乾安县为例 ［J］ .世界地质，18 (2):89 ～ 94
汤丽妮，张礼清，王卓 .2003.人工神经网络在生态环境质量评价中的应用 ［J］ .四川环境，22 (3): 69 ～72
唐剑武，叶文虎 .1998.环境承载力的本质及其定量化初步研究 ［J］ .中国环境科学，18 (3): 227 ～230
陶骏昌主编 .1994.农业预警概论 ［M］ .北京: 北京农业大学出版社，1 ～10
王根绪，钱鞠，程国栋 .2001.区域生态环境评价 (REA)的方法与应用 ［J］ .兰州大学学报 (自然科学版)，37(2): 131 ～ 140
王立群 .1991.关于建立我国森林资源预警系统的初步研究 ［J］ .林业经济，(5): 16 ～19
王良健，刘伟，包浩生 .1999.梧州市土地利用变化的驱动力研究 ［J］ .经济地理，19 (4): 74 ～79
王龙 .1995.山西煤炭开发与生态环境预警初探 (J).生态经济，(5):
王桥，杨一鹏，黄家柱等 .2005.环境遥感 ［M］ .北京: 科学出版社
王让会，樊自立 .1998.塔里木河流域生态脆弱性评价研究 ［J］ .干旱环境监测，12 (4): 218 ～323
王思敬，戴福初 .1997.环境工程地质评价、预测与对策分析 ［J］ .地质灾害与环境保护，8 (1): 27 ～34
王秀兰 .1999.基于遥感的呼伦贝尔盟农牧业土地利用变化及其对地区农业持续发展影响 ［J］ .地理科学进展，18(4): 322 ～ 330
文传甲 .1997.三峡库区大农业的自然环境现状与预警分析 ［J］ .长江流域资源与环境 6 (4): 340 ～345
文传甲 .1998.三峡库区农业生态经济系统的预警分析 ［J］ .山地研究，16 (1): 13 ～20
伍黎芝 .2000.生态脆弱区土地资源可持续利用问题 ［J］ .中国土地科学，14 (2): 13 ～16
武永峰，任志远 .2002.陕西省脆弱生态环境定量评价研究 ［J］ .干旱区资源与环境，16 (2): 10 ～14
夏军 .1999.区域水环境及生态质量评价 ［M］ .武汉: 武汉水利电力大学出版社
谢志仁，刘庄 .2001.江苏省区域生态环境综合评价研究 ［J］ .中国人口、资源与环境，11 (3): 85 ～88
徐福留 .2001.城市环境质量多级模糊综合评价 ［J］ .城市环境与城市生态，14 (2): 13 ～15
徐正春 .1993.森林资源预警系统初步研究 ［J］ .林业经济，(1): 58 ～62
许学工 .1996.黄河三角洲生态环境的评估和预警研究: (J).生态学报，16 (5): 461 ～468
阎积惠，康慧，陈怀亮 . 图像地质应用原理与方法 ［M］ .北京: 冶金工业出版社
颜卫忠 .2002.环境预警指标体系研究 ［J］ .长沙电力学院学报 (自然科学版)，17 (3): 87 ～90
杨桂山 .2001.长江三角洲近 50 年耕地数量变化的过程与驱动机制研究 ［J］ .自然资源学报，16 (2): 121 ～127
杨天柱 .2005.丹江口库区生态林现状与建设对策 ［J］ .河南林业科学，25 (2): 38 ～40
杨伟光，付怡 .1999.农业生态环境质量的指标体系与评价方法 ［J］ .监测与评价，(2): 42 ～43
杨志峰，徐俏等 .2002.城市生态敏感性分析 ［J］ .中国环境科学，22 (4): 360 ～364
姚建 . 法在县域生态环境质量评价中的应用 ［J］ .重庆环境科学，20 (2): 11 ～14
姚建，王燕，雷蕾等 .2006.岷江上游生态脆弱性的模糊评价 ［J］ .国土资源科技管理，23 (2): 90 ～92
叶文虎，栾胜基 .1994.环境质量评价学 ［M］ .北京: 高等教育出版社
叶亚平，刘鲁君 .2000.中国省域生态环境质量评价指标体系研究 ［J］ .环境科学研究，13 (3): 33 ～36
喻良，伊武军 .2002.层次分析法在城市生态环境质量评价中的应用 ［J］ .四川环境，21 (4): 38 ～40
张大任 .1991.洞庭湖生态环境预警 ［J］ .地理学与国土研究，7 (2): 42 ～44
张合平，刘云国 .2002.环境生态学 ［M］ .北京: 中国林业出版社
张明 .1999.区域土地利用结构及其驱动因子的统计分析 ［J］ .自然资源学报，14 (4): 381 ～384
张天义，张兴辽，张克伟等 .2003.遥感 .河南省国土资源综合调查与评价 ［M］ .北京: 地质出版社
张峥，张建文，李寅年等 .1999.湿地生态评价指标体系 ［J］ .农业环境保护，18 (6): 283 ～285
赵跃龙 .1999.中国脆弱生态环境类型分布及其综合治理 ［M］ .北京: 中国环境科学出版社
赵跃龙，刘燕华 .1994.中国脆弱生态环境类型划分及其范围确定 ［J］ .云南地理环境研究，6 (2): 34 ～44
赵跃龙，刘燕华 .1996.中国脆弱生态环境分布及其与贫困的关系 ［J］ .地球科学进展，11 (3): 245 ～251
赵跃龙，张玲娟 .1998.脆弱生态环境定量评价方法的研究 ［J］ .地理科学，18 (1): 73 ～79
郑新奇，王爱萍 .2000.基于 RS 与 GIS 的区域生态环境质量综合评价研究———以山东省为例 ［J］ .环境科学学报，20(4): 489 ～ 493
朱东红，上官铁梁，苏志珠等 .2000.区域生态环境质量评价方法 ［J］ .山西煤炭管理干部学院学报，(1): 64 ～67
朱发 .1999.层次分析法在石羊河流域生态环境现状评价中的应用 ［J］ .甘肃水利水电技术，(1): 66 ～68
朱晓华，杨秀春 .2002.层次分析法在区域生态环境质量评价中的应用研究 ［J］ .国土资源科技管理，18 (5): 43 ～46
庄锁法 .2000.基于层次分析法的综合评价模型 ［J］ .合肥工业大学学报 (自然科学版)，23 (4): 582 ～590
左伟，王桥，王文杰等 .2002.区域生态安全评价指标与标准研究 ［J］ .地理学与国土研究，18 (1): 67 ～71
Use and Land Cover Change: the LUCC Science Plan in Europe ［M］ .LUCC new sletter，(1): 7 ～ 10
Barnthouse role of models in ecological risk Toxicol Chem，(11): 1751 ～ 1760
Blaszczynski J Characterization with Geographic Information Engineering and Remote (2): 183 ～ 191
Boughton DA，Smith ER and O'Neill vulnerability: A conceptual Health，(5): 312 ～322
Edwards G， Uncertainty in Photointerpreted and Remote (4): 377 ～ 3
Farina and Methods in Landscape Ecology ［M］ .London: Chapman ＆ Hall，1 ～ 130
Lambin E F，Baulies X and Land-Cover Change Implementation Strategy ［A］ .IGBP Report ＆ IHDP Report ［C］ .Stochkholm: IGBP，11 ～ 16
TurnerⅡBL，Skole D，Sanderson S et and land-cover change science / research plan ［M］ .IGBP Report and HDP Report : IGBP
Turner I BL，David Skole， Steven， Use and Land Cover Change (LUCC ): Science / Research Plan［R］ .IGBP Reports，
Turner II types of global environmental changes: definitional and spatial scale issues in their human dimensions［J］ .Global Environmental Change: Humanand Policy Dimension，1 (1): 14 ～ 22
1] 冯健,周一星,程茂吉. 南京市流动人口研究[J]. 城市规划, 2001,(01) . [2] 俞德鹏. 现行城市外来人口管理方式的弊端[J]. 城市问题, 2001,(01) . [3] 顾朝林,蔡建明,张伟,马清裕,陈振光,李王鸣,沈道齐. 中国大中城市流动人口迁移规律研究[J]. 地理学报, 1999,(03) . [4] 程茂吉. 南京市流动人口的特征、影响与对策[J]. 城市问题, 1998,(02) . [5] 张启春. 谈谈进城务工人员的社会保障问题[J]. 江汉论坛, 2003,(04) . [6] 李建宁. 对强化外来流动人口综合管理的探索和设想[J]. 北京人民警察学院学报, 2001,(03) . [7] 石人炳. 浅谈对流动人口的计划生育管理[J]. 湖北大学学报(哲学社会科学版), 2000,(05) . [8] 张朝斌. 流动人口计划生育管理现状分析与对策[J]. 湖南环境生物职业技术学院学报, 2005,(04) . [9] 陈映芳. “农民工”:制度安排与身份认同[J]. 社会学研究, 2005,(03) . [10] 孟光辉. 浅析流动人口计划生育属地管理机制[J]. 市场与人口分析, 2001,(04) . [11] 辜胜阻,刘传江. 中国人口流动与城镇化的理论思考和政策选择[J]. 人口研究, 1996,(03) . [12] 辜胜阻,杨艳琳. 城市流动人口管理的思路与对策[J]. 探求, 1996,(03) . [13] 林晓红,孙跃明,朱梅芬. 流动人口现居住地管理的问题、思考和建议——江苏省昆山市个案分析[J]. 人口与计划生育, 2003,(06) . [14] 王旭. 流动人口犯罪控制论[J]. 云南公安高等专科学校学报, 2001,(01) . [15] 陈君武. 我国流动人口的新特点及其管理模式研究[J]. 云南公安高等专科学校学报, 2002,(03) . [16] 段成荣. 我国的“流动人口”[J]. 西北人口, 1999,(01) .
地理教学论文参考文献
现如今，大家或多或少都会接触过论文吧，借助论文可以达到探讨问题进行学术研究的`目的。那要怎么写好论文呢？下面是我精心整理的地理教学论文参考文献，希望能够帮助到大家。
参考文献
[1]王建,张茂恒,白世彪.圈层相互作用与自然地理学[J].地理教育,2008(4):4-7
[2]张猛刚, 雷祥义.地球表层系统浅论[J].西北地质,2005,38(2):99-101
[3]马宗晋,杜品仁,卢苗安.地球的多圈层相互作用[J].地学前缘,2001,8(1):3-8
[4]狄丽颖,冯银珍,张爱国.地球系统科学研究进展[J].河南理工大学学报,2007,26(3):267-280
[5]袁道先.对地球系统科学的几点熟悉[J].高校地质学报,1999,5(1):1-6
[6]周俊.关于地球表层的讨论[J].地球科学进展,1993,8(1):72-77
[7]谢家泽.关于地球表层系统观的几个题目[J].地球科学进展,1995,10(5):432-435
[8]叶笃正,季劲钧,严中伟,等.简论人类圈(Anthroposphere)在地球系统中的作用[J].大气科学,2009,33(3):409-415
[9]孟凡刚,孙彦顺.地球的圈层结构[J].地理教育,2009 (4):11-11 [10]卡哈尔·尼牙孜,黄建华,布玛丽娅木·阿布拉.浅议地球系统科学与可持续发展[J].西部探矿工程,2010 (5):85-87
[11]陈志明, D. A.季莫菲耶夫,V. V.勃朗古里耶夫,等.试论地貌圈概念、属性及其与其它地表圈层关系[J].地球科学进展,
其实你随便编个也没人知道当然不能太离谱我记得我当时写的时候还要5本外文参考文献,象我这种四级都过不了的人,只好充分发挥自己的想像能力编造了,结果居然通过了
参考文献著录项目：① 主要责任者(专著作者、论文集主编、学位申报人、专利申请人、报告撰写人、期刊文章作者、析出文章作者)。多个责任者之间以“，”分隔，注意在本项数据中不得出现缩写点“.”（英文作者请将作者名写全）。主要责任者只列姓名，其后不加“著”、“编”、“主编”、“合编”等责任说明。②. 文献题名及版本(初版省略)。③ 文献类型及载体类型标识。④ 出版项(出版地、出版者、出版年)。⑤ 文献出处或电子文献的可获得地址。⑥ 文献起止页码。⑦ 文献标准编号(标准号、专利号……)。如：[1] 李炳穆．理想的图书馆员和信息专家的素质与形象[J]．图书情报工作，2000(2):5-8.[2] 陶仁骥．密码学与数学[J]．自然杂志，1984，7(7):527．[3] 亚洲地质图编目组．亚洲地层与地质历史概述[J]．地质学报，1978,3:194-208．[4] DES MARAIS D J，STRAUSS H，SUMMONS R E，et a1．Carbon isotope evidence for the stepwise oxidation of the Proterozoic environment[J]．Nature，1992，359:605-609．[5] HEWITT J A．Technical services in 1 983[J]．Library Resource Services，1984，28(3):205-218．
大学本科毕业论文写作技巧
论文要求 Requirements of the Thesis
1) 论文应与本专业有关;
2) 用英文撰写;
3) 字数4000字左右(以word文档中“字数统计”栏第二行的字数为参照标准)
4) 如系翻译方面的研究，可以接受翻译评论，但不接受单纯的翻译文章
5) 学生务必在规定时间内提交正式论文定稿，过期不予受理。
3. 论文结构 Format of the Thesis
根据升达经贸管理学院外语系本专业特点，对毕业论文格式提出以下要求：
1) 所有毕业论文一律用计算机打印。具体格式见附件1;
2) 论文结构依次包括(装订顺序)：①封面(学校统一制作，由学生填写)，②英文标题，③学院、专业、学号、作者姓名、指导教师姓名，④英文论文摘要(约400字，主要概述论题的背景介绍、目的，以及研究所要解决的主要问题、可能产生的主要结果、结论和建议等)，⑤英文关键词(3 至6个)，⑥目录(可省)，⑦正文，⑧注释，(可省)，⑨参考文献，⑩中文论文题目、中文摘要、中文关键(此项单独一页);
3)所有论文要求有一定数量和分量的参考文献，文章中所引用的数据和观点应用恰当的方式表明具体出处(有关格式详见附件1的第(十)点及附件2);
4) 严禁任何抄袭或复制他人论文的不道德行为，违者按不及格处理;
5) 所有学生在提交打印定稿的同时必须提交一份word格式电子文档(文件名为：学号.DOC)，由指导教师统一交给系教学秘书。
4. 论文质量和评分标准 Quality of Thesis and Grading Standard
1) 论文评分依据：选题的理论价值和现实意义，内容的开拓性和深度、广度，分析方法，综合表达能力，文字质量，论文的格式，等。
2) 论文评分等级：成绩分为“优”、“良”、“及格”和“不及格”四个等级(每级不再分等)。所有学生优秀论文必须进行答辩，小组答辩成绩不是“优”的，小组成绩作为学生最终成绩;小组答辩成绩为“优”的，需参加系答辩委员会组织的.评优论文答辩，答辩成绩作为学生最终成绩(视情况而定)。毕业论文成绩为“不及格”的学生不能毕业。
附件1 毕业论文(设计)打印格式要求
一、毕业设计(论文)用纸、页面设置要求
毕业设计(论文)应按规定格式单面打印，纸张大小一律使用国际标准A4型复印纸。 页面设置：版心为297×210mm;页边距要求：每一面的上方(T) cm，下方(B)，左(L) cm，右(R) cm，装订线(T) cm，装订线位置(T)左，页眉，页脚，页码设置为：插入页码，居中;其余设置采取系统默认设置。
二、毕业设计(论文)内容打印要求
(一)论文题目：使用二号黑体字，加粗。居中放置
(二)学院、专业、学号、作者姓名、指导教师姓名(小四号宋体字，加粗)，依次排印在论文题下(上空二行，居中)。 ×××学院 ×××专业×××(学号) ×××(姓名) 指导教师 ×××
(三)摘要(上空二行，缩进2个汉字字符)
中文摘要采用楷体五号字，加粗，英文摘要采用五号“Time New Roman”字型，加粗;行距设置为：单倍行距。如：【摘 要】(黑体，五号字) ××××××××××××××××(楷体五号字)
(四)关键词
中文关键词采用楷体五号字，加粗，英文关键词采用五号“Time New Roman”字型，加粗;行距设置为：单倍行距。
【关键词】(黑体，五号字)××××;××××;×××;×××(楷体五号字)
(五)目录(可省项：根据实际情况确定是否需要此项)
采用五号字，其中每章题目用黑体字，每节题目用宋体字，并注明各章节起始页码，题目和页码用“„„”相连，如下所示：
目 录(黑体四号字)
(自然空二行)
1. ××××××× „„„„„„„„„„„„„„„„„„„(1)(黑体小四号字)
×××××× „„„„„„„„„„„„„„„„„„„(2) (宋体五号字)
×××××× „„„„„„„„„„„„„„„„„„„(6) (宋体五号字)
2. ×××××× „„„„„„„„„„„„„„„„„„„(40) (黑体小四号字)
(六) 正文字体要求、
每章题目居中、黑体四号字;每节题目左顶边、黑体小四字;每小节题目左顶边、黑体小四号。正文文字用宋体小四号汉字和小四号“Times New Roman”英文字体，每自然段首行缩进2个汉字字符。
(七) 行间距要求
正文行距设置：设倍行距。
每章题目与每节题目之间的行距设置：设单倍行距，段前1行、段后1 行。 每节题目与小节题目之间的行距设置：单倍行距，段后1 行。
(八) 正文章节序号编制
章，编写为：1. ，2.， 3.，„。 节，编写为：、„，、„。 小节，编写为：, „。
小节以下层次，先以括号及希腊数字为序，如(1)，(2)„;再以圈圈为序，如①, ②„。
〔注：英语写作中宜采用希腊数字为序，如(i)，(ii)„;上述的实为(1)，(2)„例证序号，而①, ②„则为脚注上标序号〕 (九) 毕业设计(论文)的公式、图与表 (略)
(十) 中文参考文献(以教务处的统一格式为准，英文的参考文献另见附件2) 参考文献(黑体，小四号字，上空二行)
[1] ×××××××××××××××××××××××××××××××××××××(宋体小四号字)
[2] ××××××××××××××××××××××××××××××××××××× 参考文献的编写方式为：
期刊文献的格式：“[编号](空一格)作者。文章题目名[类别]。期刊名,年份,卷号(期数):引文页码.”例如：
[1] 杨扬，王玉，周洲。油田储层非均质性研究[J]。高校地质学报，2000，4(3)：223-230。 图书文献的格式：“[编号] (空一格)作者。书名[类别]。出版社地址：出版社名，出版年份，引文页码。”例如：
[1] 戴功，纪律。油气储层地质学。天津：石油大学出版社，2000，110-119。
会议文献的格式：“[编号] (空一格)作者。文章题目名[类别]。论文集名[类别]。 出版社地址：出版社名，出版年份，引文页码。” 例如：
[1] 孙力，等。储层条件下水淹油层测井响应机理实验研究[A]。96国际学术讨论会论文集[C]。北京：石油工业出版社，2000,10-11。
注：1.作者姓名写到第三位，余者写“，等”;2.文献类别： A--论文集中的文章;J—期刊;C—论文集;M—书;D—学位论文;EB/OL—电子文档(这些类别代码同样适用于附件2的英文参考文献)
附件2 英文注释和参考文献格式要求(说明：本格式主要参阅了Applied Linguistics，《外语教学与研究》等杂志以及部分大学外语学院毕业论文格式要求。)
I.注释Endnotes and In-text Citations
1. 脚注(尾注)
在正文需详细注释处的右上方按顺序加注数码①②③„„，在全文之后写注文,每条加对应数码,
回行时与上一行注文对齐(注：脚注和尾注既可以用来标明、指示参考文献，也可以是对有关问题、论点做补充说明;而以下的夹注则主要针对参考文献)。
2. 夹注
某些引文或所依据的文献无须详细注释者，以夹注的形式随文在括号内注明。
1) 来自文章、专著的直接引语，作者姓名在文中已经出现： 格式：出版年份：页码
例： Rees said, “As key aspects of „in the process”(1986:241).
2 ) 来自文章、专著的直接引语，作者姓名在文中没有出现： 格式：作者姓名，出版年份：页码
例： The underlying assumption is that language “bound up with culture in multiple and complex ways”(Elli, 1968: 3).
3 ) 来自文章、专著的间接引语，作者姓名在文中已经出现： 格式：出版年份：引文页码
例：According to Alun Rees (1986:234)〔也可位于引语的最后〕, the writers focus on the unique contribution that each individual learner brings to the learning situation 〔(1986:234)〕.
4 ) 来自文章、专著的间接引语，作者姓名在文中没有提到：
格式：作者姓名 出版年份：引文页码
例：It may be true that in the appreciation of medieval art the attitude of the observer is of primary importance ( Robertson, 1987: 136). 5 ) 互联网资料： 格式：编号
例： “A deconstructive reading is a reading which analyses the specificity of a text’s critical difference from itself” (Net. 2 ). ( 注:：Net. 1, Net. 2, Net. 3,„„为了便于注明文内引语的出处，在参考文献中按编号顺序列出网址。)
II. 参考文献 Bibliography (or References)〔与夹注相对应〕
1.专著：
1) 作者的姓名(英文作者的姓，名)，书名，出版地：出版商，年份。如果同一作者两本以上同年出版的参考书，在年份后用a，b ，c 等标出。例：
Chomsky, N. Lectures on government and binding. Dordrecht:Foris. 1981a .
Chomsky, N. Theory of markedness in generative grammar. Pisa,
Italy: Scuola Normale Superiore. 1981b.
2) 书的主编：
例：Hall, David, ed. The Oxford book of American literary anecdotes[C]. New York: OUP, 1981.
3) 机构作者：
例：American library association. Intellectual freedom ed. Chicago: ALA, 1983.
4) 翻译著作：
例：Calvino, Ian. The uses of literature. (P. Creagh Trans.) San Diego: Harcourt, 1986.
2. 文章：
格式：作者姓名，篇名，刊名，刊物的卷号和期号，年份，文章的起止页码。
1) 期刊文章
例：Boling, D. The atomization of meaning[J]. Language 41, 1965:555-573
2) 论文集的文章：
例：Peters, M & T. B. Stephen. Interaction routines as cultural influences upon language acquisition[A]. In Schieffelin, B. B. & E. Ochs, eds. Language Socialization Across Cultures[C]. Cambridge: CUP, 1986: 80-96.
3. 文字资料：
1 ) 电脑软件
格式：制作人姓，名./姓名，软件名称。型号，年份。
2 ) 电影、电视、广播、录像、录音、艺术品等：
格式：片名(或节目名)，导演(或制作者)姓名。出品地(或播出台、地)，年份(或播出日期)。
3 ) 互联网资料：
格式：编号，作者名，篇名，网址
1、《工程地质学报》杂志社要求来稿请注明以下信息：作者姓名、职称/职务、工作单位、联系电话、电子邮箱、邮政编码。来稿文章字数以6000-8000为宜，摘要字数控制在300-500字左右，英文摘要控制在300 words左右。关键词3-5个。2、 图/表格：图最好为黑白色，表用三线表。图内、表内汉字为六号宋体，数字、字母为六号新罗马体。通栏图、表宽不超过15cm，半栏图、表不超过。图、表的中英文题目齐全，小五号字。坐标轴上的刻度线朝内。如：3、单位：有关记号的使用应符合国家标准，例如：sin-1应为arcsin, ctg应为cot, tg应为tan, ppm等表示法已要求停止使用；除特征数外，变量应使用单个字母表示（或带下标的单字母，否则由多个字母表示单个变量，易被误解为多个变量相乘）。4、《工程地质学报》编辑部提示注释主要用于对文章篇名、作者及文内某一特定内容作进一步的解释或说明，文中以注号（例①、②）表示。5、参考文献应按文中出现的先后顺序排列，并在文中相应处用上角标注明，未公开发表的文献请勿引用。
1.甘肃省地质局第六地质队.白家咀子硫化铜镍矿床地质.北京：地质出版社，1984,202～207.
2.汤中立，任端进.中国硫化镍矿床类型及成矿模式.地质学报，1987,61（4）:350～360.
Zhongli and Ren and metallogenic models of nickel sulfide deposits of China,Acta （2）:193～206.
Zhongli et al. Nickel deposits of Deposits of Publishing House,1992,～88.
model of the Jinchuan nickel-copper deposit,in Kirkham Duke Deposit Modeling:Geological Association of Canada,Special Paper 40,～401.
and ore deposites of the sudbury structure,Ontario,8th Iagod Sympostum[Field Trip 7].～22.
7.戈尔德施密特.地球化学.北京：科学出版社，1959,406～407.
8.戈德列夫斯基.含矿超基性岩浆的形成条件及演化.古方译.地质科技情报，1981,62～67.
9.毛景文等.桂北地区火成岩系列和锡多金属矿床成矿系列.北京：科学技术出版社，1988,8～18.