抄袭与引用的细节 “天下文章一大抄，看你会抄不会抄。”虽然俗话这么说，但对于论文来说，这个怎么抄，就显得更为严肃。抄得好，是对研究领域熟悉，论证充分，研究扎实；抄不好，就是论文剽窃，是要掉学术脑袋的，从此滚出学术江湖。但什么是抄的好，什么是抄坏事了呢？合理引用就是抄得好，抄得妙，抄得文章呱呱叫；同样是抄，如果违反了组织纪律，就会沦为剽窃抄袭，成为臭老九、丧家之犬，被人们当面吐唾沫还不敢伸手去擦。写论文是必须要参考别人成果的，在参考别人的成果时要记得感谢人家，记上人家一笔。这就是，军功章有你的一半，也有我的一半啊。具体怎么操作呢？参考别人的观点，无论是否叙述了原文，都要以注释形式标注出作者和其具体出自哪本书或者刊物，详细到出版社、出版日期、页码。如果是转引，引用的观点和资料来自别人书中的引文，建议最好查到这条引文的最初出处，核对引文，标注它的最初出处。现在论文都进入数据库了，很多人写文章根本不看原始文献，拿着别人的文章拼拼凑凑，转引时不核对材料就直接标注最初出处。大家这样反复转引，很容易出现相声“马季下了个蛋”那样以讹传讹的情况。所以，转引时务必核对引文原文，核对后再标注原始出处。当然，有些材料比较珍稀，我们找不到原文，只能转引别人的引用。这种情况在引用时可以标注转引。比如：***，出处，转引自***引用只要标明了出处，就再也不怕沦为抄袭了。 （资料来源：中国论文网）
从事“凝聚态物理”和“材料物理与化学”方向的研究，围绕微纳结构材料研究开展科研。从事微纳加工工艺、薄膜技术和纳米材料合成工艺研究。 目前主要研究： 1．纳米多孔气凝胶类材料的制备与应用 2．微纳加工技术应用 奖励荣誉获奖 气凝胶材料的微结构控制机理与特性，2010年教育部自然科学奖二等奖，第三完成人 碳及掺碳气凝胶的结构形成机理与特性，2008年上海市自然科学奖二等奖，第三完成人 结构可控的纳米多孔气凝胶薄膜的研制及光学特性研究，2005年上海市科技进步奖二等奖，第三完成人 惯性约束聚变（ICF）分解实验用的平面调制靶和平面薄膜的研制，2004年上海市科技进步奖三等奖，第一完成人 荣誉 2008年教育部新世纪人才 2008年上海市青年科技启明星（跟踪） 2005年上海市青年科技启明星 2004年上海市高校优秀青年教师 学术兼职中国硅酸盐学会溶胶凝胶分会，副秘书长 中国核学会核靶专业委员会，副主任 主要论著论文：（1-8, 30更新） 杜艾，周斌（通信作者）等，空间探索用梯度密度气凝胶的合成与应用进展航空学报，2011,32:961- 归佳寅，周斌（通信作者）等，溶胶共凝法制备密度渐变SiO2气凝胶及界面研究航空学报，2011,32:941- 钟艳红，周斌（通信作者）等，密度渐变多层碳气凝胶靶型的制备研究强激光与粒子束，2011,23(3):657- 朱秀榕，周斌（通信作者）等，碳气凝胶/聚苯乙烯双介质柱状靶的研制，强激光与粒子束，2011,23 (7):1843- 陈珂，周斌（通信作者）等，多孔硅的模板限制还原法制备及性能 物理化学学报，2011,27(11):2719- 刘光武，周斌（通信作者）等，水玻璃为源的超疏水型SiO2 气凝胶块体制备与表征 硅酸盐学报，2012,40:160- 杜艾，周斌（通信作者）等 高效X光背光源靶材料的研究进展 原子能科学技术, 2012, 46(7): 842- 沈洋，周斌（通信作者）等 用于冲击波延时的密度梯度气凝胶飞片靶的制备及表征 原子能科学与技术, 2012, 46(6): 734- 李龙翔，周斌（通信作者）等 金套筒内SiO2 气凝胶原位成型技术研究 原子能科学与技术, 2012, 46(6): 731-
张玉君（地质矿产部航空物探遥感中心，北京100083）摘要：用图像处理技术对下述与油气田异常有关的问题进行了研究：快速提取高频弱磁异常，估算高频磁异常源的埋深，放射性异常性质的归一化处理，典型异常剖面的分析，以及油气远景评价，提出了4个有希望含油气的异常。经验证，这4个异常中有3个见到油气显示，其中1个为工业气流。关键词：油气田，烟囱效应，图像处理，航放Th归一化处理，油气远景评价。一、前言近10年，非地震物探方法寻找油气田的工作越来越受到人们的重视[1-8]，原因之一是地震勘探生产费用日益昂贵；另一方面，近10年航磁、航放测量技术发生了重大变化，航空光泵磁力仪的灵敏度从1nT提高到1～01nT，大体积方柱形NaI（Tl）晶体和数字电路被引入航空放射性测量仪，不仅使其灵敏度提高了10多倍，而且使强度测量被4道（K、Th、U及总道）或多道（256或512道）能谱测量所取代。这一切促使物探工作者在油气勘查中更加注意应用具有大面积低消耗特性的航空物探技术。柴达木中部地区的高灵敏度航磁、航空伽马能谱测量飞行线距1km，飞行高度约90m；安装在双水獭型飞机上的航放、航磁综合站包括GAD-6型4道能谱仪、MAP-4型质子旋进磁力仪、磁带收录机及多普勒雷达、航空摄影两种导航定位设备。能谱仪使用了3箱NaI（T1）晶体，总体积为48L，磁力仪的灵敏度为5nT。放射性伽马测量普查油气已有50多年的历史，实践表明，在油气田上方有偏低的放射性异常，并在油气构造边缘存在比本低高的异常峰值。1979年美国Donovan TJ等人在Cement油田发现了高频磁异常[13],10多年来人们公认了这种与油气田相关的高频磁异常[9]的存在。尽管对产生这些异常的机制存在不少争议，但最被接受的假说是烟囱效应，通常认为几乎所有的储油构造都从深部排出水份，烟囱效应将碳氢化合物和其他有机衍生物通过喷溢、溶液散逸、气体扩散等三种渗漏作用带到地表；地层中的液体和气体在从油、气层向上运移过程中，由于后成作用形成磁铁矿并沉淀放射性元素，从而形成可被航测发现的航磁航放异常，并作为判断地下赋存有碳氢化合物的标志。提出用图像处理技术来研究这一问题，基于两点考虑：航空物探的大面积宏观数据通过图像处理可以快速获得直观结果；图像系统具备既适合位场数据（航磁）又适合统计学性质数据（航放）的处理功能。本文所做基础图像网度为5×5m2。除航磁、航放外还利用了陆地卫星影像 MSS（7、5、4波段）数据。二、异常特征研究高频弱磁异常的快速提取柴达木盆地是我国重要的含油气盆地之一，在测区内已有多处井下油气显示，测区内人烟罕见，人文干扰较少，为研究油气田后成磁铁矿引起的高频弱磁异常特征的理想试验场地。测区磁场动态范围为0～533 nT，在0～127 nT的低端仅有12个像素，占全图面的5/100000，并集中在东北部的局部负异常内，故采取限幅拉伸的方法：张玉君地质勘查新方法研究论文集式中 z（x,y）为磁场异常，单位nT,l（x,y）为图像灰度级或称灰阶。每一灰阶代表598 nT，可见经限幅拉伸处理，使磁原始图像保留了更多的弱异常信息，再按（2）式[10]制作立体阴影图（彩版附图9（A））。张玉君地质勘查新方法研究论文集式中 θ为光源方位角；φ为光源高度角；λ为磁面法线矢量与光源矢量的夹角。测区磁场彩色立体阴影图显示出沿着从西北角到东南角对角线两侧成群地存在着高频弱磁异常，其宽度为1～2km，幅度为2～15 nT，在黑白图像上有3～20灰阶的变化。文献[2]认为，钻井套管和井盖可能引起幅度为2～5 nT的异常，但其波长更短（＜300m），也就是说在5×5m网度的图像上，套管和井盖所引起的异常仅表现为孤立点。为了突出有用弱磁异常，用图1中的模板通过褶积功能制作了水平二次导数图像，通过邻域滤波消除孤立异常点。彩版附图9（B）为航磁二次导数、经归一化处理U及油气异常分布叠加图像。青色点群为航磁水平二次导数异常；红色为经归一化处理 U；绿色和蓝色椭圆为已知异常位置；红色椭圆为新发现的有希望异常。10个已知油气显示全部都具有高频弱磁异常，这有力地表明这种仅用数分钟即可完成的磁场二次水平导数图对于提取与油气田有关的高频异常具有重要意义。图1 求导褶积模板高频弱磁异常的磁源埋深对磁场（经镜像对称扩边）进行快速富氏变换，然后对变换结果计算功率谱，从功率谱图像上取径向剖面值，最后做径向功率谱图（图2），为了做图方便直接用图像上的点数作为横坐标。该谱图表明全图由低频功率谱Pl和高频功率谱Ph组成，低频部分在该区主要由前泥盆系变质岩的磁性引起，高频部分主要由浅部后成作用局部矿化引起，用参数法求磁源埋深[14]张玉君地质勘查新方法研究论文集取对数张玉君地质勘查新方法研究论文集图2 径向功率谱如果用半对数坐标成图，式（4）为线性方程，－2H与－2h分别为低频段（1～8点）与高频段（12～50点）的直线斜率，lna与lnb分别为二直线的截距，△r为径向频率间隔。r=2π/（点数×点距）（5）本测区点数为512，点距为5 m，∆r=65×10-6m，用最小二乘法求得二直线的截距和斜率分别为a=1429,Kt=-5595, b=7212,Kh=-9182，张玉君地质勘查新方法研究论文集式中S被导入是因为图像处理所显示的对数功率谱值仅有相对意义，S的单位为m，其值可通过某已知异常源埋深求得。已知低频磁异常系由埋深为－5000至－15000m的前泥盆系变质岩引起，则张玉君地质勘查新方法研究论文集将H代入，S为0285～0855m。张玉君地质勘查新方法研究论文集即高频弱磁源埋深在3～1km之间。航放图像对油气藏的反映航放三元素合成图像[11-12]示于彩版附图9（C），叠加于其上的还有航磁二次水平导数异常。红为K，其含量变化为0～8%；绿为Th，其含量变化为0—39×10-6；蓝为U，其含量变化为0～9×10-6，彩版附图9（D）的上图为其色标；黄色点群为航磁水平二次导数异常。航放图像较之航磁水平二次导数对于油气显示的相关性要复杂得多，虽然并非所有已知油气显示都有清楚的航放异常，但还是有一些含油气构造反映较清楚，如彩版附图9（C）上的5号异常为已知鸭湖含油气穹隆。在双元素图像（彩版附图10（A），彩版附图9（D）的下图为其色标）上鸭湖构造也十分明显；在该图上部为R—Th,G—K,B—K合成图像，箭头所指为已知东台吉乃尔含油气构造；下部为R—Th,G—U,B—U合成图像，箭头所指为已知涩北1、2号含油气构造，它们都表现为周边高中间低的放射性异常。如何增强或突出航放数据中与油气有关的信息，关键在于减少或压制岩性影响，为此，我们采用了主分量分析（在图像处理中称之为KL变换）和Th归一化处理[15]。彩版附图10（B）为航放三元素主分量分析（即KL变换）图像，彩版附图9（D）的上色标亦为其图例；图中a、b、c分别为第1、第2、第3主分量，以KL（1）、KL（2）、KL（3）表示，d为其合成图像。各主分量的组成为：张玉君地质勘查新方法研究论文集式中，Xk、XTh、X。分别代表某像素K、Th、U之含量值。在第1主分量ξ1中以K为主，减少了Th的作用；由于在大部分沙岩、砾岩中K和Th有正相关性，故经此处理的KL（1）大大压制了岩性的影响，突出了K与Th不相关的地质过程的反映，如钾盐的沉积及烟囱效应等。ξ2综合了K、Th、U的变化，反映全区地质环境，对研究油气局部异常作用不明显。ξ3中以U为主，又减弱了Th和K的作用，大大压制了岩性影响，对研究油气局部异常很重要。为减少岩性影响，根据文献[15]，在整理地面或航空放射性测量数据时，常常除以Th含量值，称之为Th归一化处理。在做这种处理时，认为Th能较多地代表砂砾沉积物中岩性的变化，K、U与Th在这些沉积物中有一定正相关关系，除以Th值便减少了或压制了岩性的影响。我们的改进在于在归一化处理之前，先做相关分析，求出相关系数，再按相关系数做Th归一化处理。彩版附图10（C）即为K和U经归一化处理前后的对比图像，左边为原始图像，彩版附图9（D）的上图为其色标；右边为归一化处理结果。表1为已知及新发现异常统计表，共8项特征：K原始图像、K经KL变换、K经Th归一化处理、U原始图像、U经KL变换、U经Th归一化处理、U晕圈及MSS（7、5、4波段）图像上的晕圈反映。已知10处油气显示中有8个（占80%）有明显或较明显的航放异常；其典型特征是：在Th归一化图上中心部位K低、U低和U高值的外围晕圈。Th归一化处理的作用对下述含油气构造的异常十分明显：使V1、V2、V5的K高变为K低，使V1、V5的U高显著降低。故Th归一化处理对于统一异常性质有明显效果。典型异常的剖面图像以鸭湖（V2）和红三旱四号（V1）做为典型的已知异常进行剖面研究。鸭湖异常有以下典型特征（彩版附图10（D））:（1）U在归一化处理前后（彩版附图10（D）右图中之蓝色曲线）异常中心部位都是低值；（2）K在异常中心部位原为高值，经归一化处理后也获低值（彩版附图10（D）右图中之红色曲线）；（3）在异常周围U有明显的高峰和晕圈现象；（4）有明显的高频弱磁异常，且中心处更强（彩版附图10（D）右图中之绿色曲线）。表1 已知及新发现异常统计表红三旱四号的异常特征为：（1）U在归一化处理前异常中心部位为强高值，经归一化处理后显著降低；（2）K在归一化处理前为高值，经归一化处理后为低值；（3）经归一化处理后在异常周围U有明显的高峰，反映出晕圈现象；图3 红三旱四号地面能谱测量剖面（4）也有明显的高频弱磁异常。综合典型异常分析，与油气藏有关的航磁航放异常识别标志为：航磁局部高频弱变化；Th归一化处理后局部低钾、低铀（也有达不到低值）；以及铀在异常周边的高值晕圈现象。我们对红三旱四号异常进行了地面检查，用4道能谱仪进行了剖面测量，总长10km，点距100～200m。所获剖面（图3）与空中测量结果一致，在异常中心部位为K高、U高和Th高，从中心向两侧先降低后又升高。三、对含油气远景评价的讨论本测区内已知含油气构造或点共10处，如表1所列：红三旱四号（V1，气），鸭湖（V2，气），马新高点（V4，油），东台吉乃尔（V5，见工业油气流），南陵丘（V6，油气），驼峰山（V7，油气），涩北1号（V8，油气），涩北2号（V9，油气），涩聂北（V12，气），涩聂东（V3，气）。这10个已知油气显示无一例外都有明显的航磁水平二次导数异常。前8个有明显或较明显的航放异常显示，表现为在经Th归一化处理的K元素图像和U元素图像上为中心部位的低值及U的周边升高的晕圈现象。后两个已知区航放异常不够明显，主要原因是，这两个已知油气显示点处于钾盐类沉积丰富地段，钾引起的放射性变化强于与碳氢化合物有关的弱放射性异常，后者未能明显表现出来。故本测区已知油气显示构造或点与航磁水平二次导数异常的吻合率为100%，与航放K、U异常的吻合率为80%。李芦玲等李芦玲等柴达木盆地中部地区航空磁测详查结果报告地质矿产部航空物探遥感中心，1985。曾用传统的方法对本测区含油气远景进行了评价，圈出8个最有希望的含油气区，图像处理也完全证实了这一评价。经图像处理，我们认为还应补充提出以下4个有希望的含油气异常区：（1）依克雅乌汝背斜构造（V3），位于鸭湖含油气背斜构造东30km，台吉乃尔含油气背斜北15km。其航磁航放、卫片异常均很典型，且与鸭湖异常极为相似，异常特征依据充分，可做为新发现的最有希望的含油气异常点。（2）南陵丘西异常（V10），位于南陵丘已知含油气构造西约5km。有局部高频低幅磁异常、低钾、低铀及较明显可辨之周边闭合的高值U晕圈。可做为新发现有希望的含油气异常。（3）洪上异常（V11），位于台吉乃尔洪积扇上部。除有局部高频弱磁异常外，有低U及隐约可见的周边U高值晕圈，可做为新发现较有希望的含油气异常。（4）南八仙异常（V14），位于仙西南高点北缘。有局部高频弱磁、低钾和低铀异常。也可做为新发现之较有希望的含油气异常。如前所述彩版附图9（B）上展示了l0处已知（蓝、绿、黄色椭圆）和4处新发现油气异常（红色椭圆）。并与归一化处理后之U（红）以及航磁水平二次导数异常（青点）叠加在一起。上述远景评价于1990年完成，根据能源部1991年在柴达木中部地区进行的勘探成果，对比如下，上述4个含油气远景区中的3个见到了油气显示：即依克雅乌汝背斜构造（V3）于“伊中1”井见油气显示，洪上异常于“台中1”井见工业气流，南八仙异常于“仙3”井见油气显示。本工作所预测的4个含油气远景区中有3个见到了油气显示，说明航磁、航放等综合参数通过图像处理所做油气田异常特征的研究是有效的。常规方法没有发现的异常，却由本方法预测，并已得到验证，说明该方法的价值。参加野外验证的还有水恩海、史殿林、郭毅，屏幕图像均由杨星虹拍摄，能源部石油勘探开发科学研究院遥感地质所叶和飞给予了支持。在此表示诚挚的谢意。参考文献[1]Wold,R郭武林整理利用航磁和放射性测量直接寻找油气藏地质科技动态，1984,（24）:21～[2]管志宁利用航磁探测与化学剩磁有关的油气矿床国外地质勘探技术，1985,（7）:18～[3]王家林，万明浩，金国英磁测找油的进一步试验及有关问题的研究石油地球物理勘探，1985,（20）:397～[4]郭玉琨编译苏联地球物理直接找油方法概况石油物探译丛，1985,（1）:90～[5]王秀文摘译石油和天然气放射勘探的新远景地质科技动态，1984,（18）:17～[6]赵改善编译直接勘探油气田的放射性方法石油物探译丛，1987,（6）:95～[7]李淑仪核技术在石油天然气和地热勘探中的应用国外地质勘探技术，1986,（4）:22～[8]王锡田编油气勘查中用遥感方法检测烃类微渗漏地质科技动态，1987,（11）:13～16[9]蔡振京高精度大比例尺（1∶50000）航空物探在油气藏勘探中的应用物探与化探，1989（13）:401～[10]Teskey,D,Broome,JComputer programs for production of shaded relief and stereo shaded relief GSC 84—1 B,375—389,1984[11]Zhang YApplication of image processing techniques to airborne radiometric Abstract of the Second Symposium on Exploration Geophysics,Xian,China,525—526,1986[12]Zhang YDigital image processing of airborne radiometric and magnetic data from central chaidam An Overview of Exploration Geophysics in China,517—535,1989[13]Donovan,TJ,Forgey, RL, Roberts, AAAeromagnetic detection of diagenetic magnetite over oil Bulletin ofthe American Association ofPetroleum Geologists,63,245—248,1979[14]Donovan,TJ,Brien, DPBryan,JG, Cunningham, KINear surface magnetic indications of buried Aeromagnetic detection and separation of spurious signals, IGARSS' 87 Remote Sensing: understanding the earth as a system,219—232,1987[15]Alton VGRadiometrics—a practical exploration Weotern Oil Reporter'March,173,A STUDY OF AERO-MAGNETIC AND AERORADIOMETRIC ANOMALIES IN SOME OIL-GAS FIELDS BY IMAGE PROCESSING TECHNIQUESZhang Yu jun（Center of Aero-Geophysics and Remote Sensing, Ministry of Geology and Mineral Resources , Beijing 100083）AbstractThis work represents the results of study of the following problemsby image processing techniques: the quick extraction of weak magnetic anomalies with high frequency, the evaluation of depth for highfrequence magnetic anomalies,the normalization of radiometric anomalies, the analysis of typical anomalous profiles and the prediction of perspective anomalies to find oil and Four new perspective anomalies for oil gas were The study had been finished at the beginning of Comparing with the results of exploration in 1991, we have been encouraged by the fact that not only on 3 out of the 4 anomalieas the oil and gas werefound, but also on one of them the gas was Key words Oil-gas filed,The chimney effect,Image processing, The second holizontal gradient of magneticfield, The power spectrum, The Th normalization for aeroradiometrics, The perspective evaluation for oil-原载《地球物理学报》，1994,V37,N1。
陈小红1 潘懋1 徐旸2 汪艳梅3(北京大学地球与空间科学学院；北京市地质调查研究院；内蒙古自治区国土资源信息院)摘要 进行矿权范围内涉及国家出资地质勘查情况核查是矿政管理的一项重要工作。本文基于GIS技术，采用矿政空间数据与地质资料的结合，实现了矿权范围内涉及国家出资地质勘查情况筛查、重叠区域的核实等功能。该系统在内蒙古自治区国土资源厅实际工作中得到了应用，极大地提高了核查工作效率。关键词 国家出资地质勘查 GIS 空间分析 影像配准0 引言国家出资地质勘查，是指地质勘查资金由国家投入、地质勘查风险由出资者(国家)所承担的地质勘查。国家投入的地质勘查资金来源一般包括税收、矿产资源补偿款、矿业权使用费及价款、地质勘查基金(中央或者地方)[1～2]。国家出资地质勘查形成的成果以不同形成体现在矿业权中，因此在矿业权登记、转让等环节中，需要考虑国家出资地质勘查情况。矿政管理一般要承担矿业权的登记、变更、注销等工作以及矿产储量登记工作等。GIS技术已经广泛应用于矿政管理系统中，矿权管理系统、储量空间数据库管理系统、矿业权审批系统等[3～7]主要实现数据填报、矿权登记业务流程管理、网上审批、基于GIS功能的图形查询、统计、发布等功能；GIS应用于“一张图管矿”是业内研究的热点，主要将矿业权实地核查数据与矿产资源规划等多个数据库进行集成，实现矿权审批、监管、统计分析等[8～11]。本文基于GIS技术，结合地质档案资料，集成探矿权、采矿权、储量等空间信息，实现了矿权范围内涉及国家出资地质勘查情况筛查、重叠区域的核实等功能，为矿业权登记、转让等环节中价款核算提供支持。该系统在内蒙古自治区国土资源厅工作中得到应用，解决了传统核查中核查精度低、周期长的问题，提高了矿政管理效率。1 系统框架及功能1 设计思路矿业权(探矿权、采矿权、划定矿区范围批复等)审批、价款评估与缴费是矿政管理的一项重要工作，而进行矿业权范围内国家出资地质勘查情况核查工作是矿政管理及价款评估与收缴的重要依据之一。基于GIS的国家出资地质勘查核查系统建设应满足矿政管理的多层次需求，实现涉及矿政管理的基础数据管理、查询、图形核查、报告管理、统计分析等需求。系统建设应遵循总体设计、分步实施、技术先进、安全可靠的基本思路。总体设计，分步实施：系统设计按照矿政管理需求，全面统筹规划，做好顶层设计，进而分步骤实施。优先考虑系统架构的合理设计，在此基础上逐步进行各功能组件及系统接口的深化设计，使组件之间实现高内聚、低耦合，力求形成一个布局合理、功能完备、能力均衡、分工明确的平台。技术先进，安全可靠：系统建设充分借鉴信息技术发展前沿技术，综合运用数据库、地理信息系统(GIS)、分布式计算等多项信息技术，并遵循稳定性(运行稳定)、可扩展性(功能、数据、可视化、用户界面等)、可配置性(参数可配置)、标准性和规范性等设计原则。2 系统框架基于GIS的国家出资地质勘查核查系统框架包括数据库、技术支撑、应用系统、用户等，如图1所示。利用地质资料目录数据库、地质资料图文数据库、探矿权登记数据库、矿业权登记数据库及储量空间数据库，提取相关的数据项，构建数据库；采用空间数据库引擎进行数据访问；基于GIS平台二次开发，实现空间检索和图形核查等功能；基础组件实现显示、可视化、报表、报告模板、输出等等功能。图1 基于GIS的国家出资地质勘查核查系统框架3 主要功能基础数据管理：主要包括地理底图、地质资料目录数据、探矿权登记数据、矿业权登记数据、储量空间数据管理。除了具有数据导入、录入、修改、删除等功能外，还包括利用区域坐标(拐点坐标)生成图层功能，生成的图层为地质资料目录GIS图层、探矿权GIS图层、矿业权GIS图层、储量GIS图层；另外，区域坐标(拐点坐标)发生变化后，GIS图层将实现更新。条件查询：包括地质资料目录数据、探矿权登记数据、矿业权登记数据、储量空间数据查询。空间检索：除了具有点选检索、拉框检索、输入地理坐标检索功能外，还有两个图层交互检索功能。图形核查：功能包括两个方面：一是在地理底图基础上加载探矿权图层、矿业权图层、储量信息图层及地质资料目录图层，实现待核查矿权范围内地质资料检索，检索结果显示、追加、删除及资料查看等功能；二是实现待核查矿权与涉及国家出资地质勘查重叠区域比对功能，包括导入数据、配准、距离量算、面积量算、数据导出等功能。核查报告管理：核查信息保存到数据库后，系统具有提取核查信息并自动生成报告功能，并且还具有报告上传入库、查询、显示等功能。统计分析：按照矿权性质(招拍挂、探矿权、采矿权、划定矿区范围)统计，按矿权类别(探转采、延续、变更、转让、新立)统计，按行政区划查询(矿权所在的省、市、区行政区划)统计。2 系统主要关键技术及实现系统采用Oracle10 g数据库，把地质资料目录数据库、探矿权登记数据库、矿业权登记数据库及储量空间数据库提取的数据项按照空间-属性进行关联；系统采用C＃ WinForm，基于GIS平台二次开发；实现系统的关键技术主要为矿权涉及国家出资地质勘查项目筛查、矿权涉及国家出资地质勘查重叠区域核实、核查信息提取并生成核查报告等方面。1 空间数据与属性数据组织地质资料是一种国家投入巨额勘查资金获得的特有信息资源[3]，是将地质勘查形成的原始地质资料、实物地质资料等进行总结、编制而形成，地质资料需按照地质汇交管理办法进行汇交[12]。汇交的地质资料包含目录和电子文件。利用地质资料目录数据库、地质资料图文数据库、探矿权登记数据库、矿业权登记数据库及储量空间数据库，提取相关的数据项，再增加与核查特性相关的数据项，组成核查系统的数据内容。提取地质资料目录数据库中档号、报告名称、起始经度、终止经度、起始纬度、终止纬度等数据项，并增加是否涉及国家出资数据项；提取地质资料图文数据库的档案号及电子文件；提取探矿权登记数据库中矿权证号、项目名称、矿权性质、申请人、勘查单位、地理位置、东经起、东经止、北纬起、北纬止、区域坐标(拐点坐标)等数据项，并利用坐标信息组织空间数据图层；提取矿业权登记数据库中矿权证号、矿山名称、矿权性质、申请人、勘查单位、地理位置、东经起、东经止、北纬起、北纬止、区域坐标(拐点坐标)等数据项，并利用坐标信息组织空间数据图层；提取储量空间数据库中矿区名称、报告名称、行政区域、主矿种、矿区面积、地理位置、东经起、东经止、北纬起、北纬止、区域坐标(拐点坐标)等数据项，并利用坐标信息组织空间数据图层。另外，设计了保存核查结果的数据项，主要为：文号、项目名称、项目性质、矿权证号、矿权人、面积、行政区域、东经起、东经止、北纬起、北纬止、区域坐标(拐点坐标)、报告出具时间、是否国家出资等。2 矿权涉及国家出资地质勘查项目筛查在地理底图上，加载探矿权图层、采矿权图层、储量信息图层、地质资料目录图层，利用GIS的空间分析功能[13，14]，实现矿权范围内涉及国家出资地质勘查项目，具体步骤为：1)建立图形工作区，图形工作区为西安80 地理坐标，并加载探矿权图层、采矿权图层、储量信息图层、地质资料目录图层，其中，探矿权图层、采矿权图层、储量信息图层一般为高斯三度带或者高斯六度带投影，需进行坐标转换。2)地质资料目录图层中属性“是否涉及国家出资”数据项进行标识，分别标为“涉及国家出资”和“不涉及国家出资”。3)在探矿权图层、采矿权图层或储量信息图层中，选定待核查矿权，利用GIS空间分析功能，得出该矿权范围所有的地质勘查项目，如图2所示。4)利用“是否涉及国家出资”数据项进一步筛查，即可得出矿权范围内涉及国家出资地质勘查项目。5)筛查得出的涉及国家出资地质勘查项目，如有其他资料进行佐证，可以追加或者剔除。图2 待核查矿权涉及地质勘查项目示意图3 矿权与国家出资地质勘查重叠区域核实筛查得出的涉及国家出资地质勘查项目，与图文数据库链接，提取国家出资地质勘查项目数据，并利用GIS的矢量数据与影像数据配准功能[15～17]，实现矿权与国家出资地质勘查重叠区域核实。具体步骤为：1)建立图形工作区，图形工作采用高斯三度带或者高斯六度带投影，并加载待核查矿权图形。2)获取国家出资地质勘查工作区域影像数据，并以待核查矿权图形为参考进行配准，如图3所示。3)配准后，在图形工作区域利用距离或者面积量算工具，确定待核查矿权与涉及国家出资地质勘查工作区比例。图3 待核查矿权图形与涉及国家出资工作区域影像配准4 核查信息提取，自动形成核查报告核查中形成矿权人、矿权、位置(坐标)、面积、矿权范围涉及地质勘查情况等信息，提取后存入数据库中。利用Word模板，系统实现核查报告自动生成，该报告为矿权价款核算提供重要依据。3 系统应用该系统在开发完成后，依托内蒙古自治区矿政管理基础数据库及管理系统建设项目，在内蒙古自治区国土资源厅实际工作中得到应用，极大地提高了核查工作效率。图4为该系统主界面。图4 基于GIS的矿权涉及国家出资地质勘查核查系统工作界面以前核查一宗矿权涉及国家出资地质勘查需要调阅20 档左右的纸质资料，并需要采用手工的方式，将坐标转换后投影到矿权纸图上扣合后得出结论，这种方式费时、费力，效率低下。而该系统实现了核查流程自动化，利用GIS的空间分析功能，快速得出涉及国家出资地质勘查情况，而且可以直接调阅地质资料，并在计算机屏幕上进行核查，减少了大量反复调阅纸质地质资料的时间，核查效率提高30%以上。4 结论本文基于GIS技术，结合矿政信息管理实际需求，实现了矿权范围内涉及国家出资地质勘查核查功能。该系统在内蒙古自治区国土资源厅应用中取得了很好的效果，利用该系统，累计完成了超过近千宗的矿权核查，提高了矿政管理能力，取得了良好的社会和经济效益。参考文献[1]孙习稳，干飞国家出资地质勘查形成矿业权的经济关系研究[J]中国矿业，2007，16(5)：13～[2]徐婧，李慧洁，薛建明国家出资地质勘查形成的矿业权产权归属与收益分配[J]科学论坛，2012(5)：159～[3]赵善人，周桅，等矿业权管理信息系统建设[J]国土资源信息化，2003(2)：3～[4]周桅，赵善仁，等矿业权管理信息系统技术方案及其实现[J]国土资源信息化，2003(2)：7～[5]吴仲煜，查宗祥，等矿业权管理信息系统中图形子系统的设计及特点[J]国土资源信息化，2003(2)：11～[6]孙立双，马晓明，等矿业权登记管理系统图形辅助审查系统功能设计[J]国土资源信息化，2013，34(11)：95～[7]胡延晟矿业权项目审批管理信息系统的设计与实现[D]长春：吉林大学，[8]刘金星，李新通，等基于网络 GIS的福建矿业权图形辅助审查系统设计与建设[J]山西师范大学学报(自然科学版)，2011，25(1)：116～[9]杨文森，陆世东，等“一张图管矿”的数据组织与分类[J]地理空间信息，2012，10(1)：64～[10]杨文森，胡凯，等湖北省“一张图管矿”试点的研究与应用[J]国土资源情报，2010，(10)：70～[11]谭永杰，郭佳，等基于全国矿业权实地核查成果数据的矿政管理信息系统建设[J]中国矿业，2011，20(7)：50～[12]於顺然浅谈新的地质资料汇交内容及其制度[J]江苏地质，2002，26(3)：180～[13]张桂兰，杨文海浅谈《内藏古矿产资源储量空间数据库系统》建设与应用西北资源[J](2)：40～[14]吴立新，史文中编著地理信息系统原理与算法[M]北京：科学出版社，[15]郭仁忠著空间分析[M]武汉：武汉测绘科技大学出版社，[16]宋智礼图像配准技术及其应用的研究[D]上海：复旦大学，[17]魏晓敏图像配准算法研究与系统设计实现[D]南京：南京航空航天大学，
一般的文艺类的报刊杂志的稿件处理期为3个月，3个月没有得到通知，基本上可以自行处理，现在的杂志一般不退稿。如果是学术性论文，那就不好说了，因为现在潜规则在中国非常流行，可以找一下负责的编辑，让对方开出条件，否则，被毙的机率超高。一般正规的报刊，接到了好的稿件，第一时间会通知作者，告知可能采用意向。
《中国民用航空学报》属于一号多刊，不是正规期刊。但是有的院校检查不严，可以用。有的学校查的严格，就用不了啊。空中交通管理审稿是一个月。
ppt常用的话，最好做成模板首先准备好要做成模版的图片，打开powerpoint并新建一个空白的ppt文档。视图→母版→幻灯片母版，进入母版编辑状态。点击绘图工具栏上的“插入图片”按钮（或按“插入→图片→来自文件”），选中要做为模版的图片，确定。并调整图片大小，使之与母版大小一致。在图片上点击鼠标右键，叠放次序→置于底层，使图片不能影响对母版排版的编辑。现在我们已经制作好了我们需要的模版文档，退出母版视图后，接下来就是保存了。文件→保存，打开“另存为”对话框，在“保存类型”中选择“演示文稿设计模版”此时程序将打开默认的文件保存位置，不用更改它，在“文件名”中输入一个便于你自己记忆的名字，确定保存下来。现在关闭此ppt文档。再新建一个空白文档，看看你的“幻灯片设计”里，刚刚做好的模版文档是不是已经出现在那儿啦。
张玉君（地质矿产部航空物探遥感中心，北京100083）摘要：用图像处理技术对下述与油气田异常有关的问题进行了研究：快速提取高频弱磁异常，估算高频磁异常源的埋深，放射性异常性质的归一化处理，典型异常剖面的分析，以及油气远景评价，提出了4个有希望含油气的异常。经验证，这4个异常中有3个见到油气显示，其中1个为工业气流。关键词：油气田，烟囱效应，图像处理，航放Th归一化处理，油气远景评价。一、前言近10年，非地震物探方法寻找油气田的工作越来越受到人们的重视[1-8]，原因之一是地震勘探生产费用日益昂贵；另一方面，近10年航磁、航放测量技术发生了重大变化，航空光泵磁力仪的灵敏度从1nT提高到1～01nT，大体积方柱形NaI（Tl）晶体和数字电路被引入航空放射性测量仪，不仅使其灵敏度提高了10多倍，而且使强度测量被4道（K、Th、U及总道）或多道（256或512道）能谱测量所取代。这一切促使物探工作者在油气勘查中更加注意应用具有大面积低消耗特性的航空物探技术。柴达木中部地区的高灵敏度航磁、航空伽马能谱测量飞行线距1km，飞行高度约90m；安装在双水獭型飞机上的航放、航磁综合站包括GAD-6型4道能谱仪、MAP-4型质子旋进磁力仪、磁带收录机及多普勒雷达、航空摄影两种导航定位设备。能谱仪使用了3箱NaI（T1）晶体，总体积为48L，磁力仪的灵敏度为5nT。放射性伽马测量普查油气已有50多年的历史，实践表明，在油气田上方有偏低的放射性异常，并在油气构造边缘存在比本低高的异常峰值。1979年美国Donovan TJ等人在Cement油田发现了高频磁异常[13],10多年来人们公认了这种与油气田相关的高频磁异常[9]的存在。尽管对产生这些异常的机制存在不少争议，但最被接受的假说是烟囱效应，通常认为几乎所有的储油构造都从深部排出水份，烟囱效应将碳氢化合物和其他有机衍生物通过喷溢、溶液散逸、气体扩散等三种渗漏作用带到地表；地层中的液体和气体在从油、气层向上运移过程中，由于后成作用形成磁铁矿并沉淀放射性元素，从而形成可被航测发现的航磁航放异常，并作为判断地下赋存有碳氢化合物的标志。提出用图像处理技术来研究这一问题，基于两点考虑：航空物探的大面积宏观数据通过图像处理可以快速获得直观结果；图像系统具备既适合位场数据（航磁）又适合统计学性质数据（航放）的处理功能。本文所做基础图像网度为5×5m2。除航磁、航放外还利用了陆地卫星影像 MSS（7、5、4波段）数据。二、异常特征研究高频弱磁异常的快速提取柴达木盆地是我国重要的含油气盆地之一，在测区内已有多处井下油气显示，测区内人烟罕见，人文干扰较少，为研究油气田后成磁铁矿引起的高频弱磁异常特征的理想试验场地。测区磁场动态范围为0～533 nT，在0～127 nT的低端仅有12个像素，占全图面的5/100000，并集中在东北部的局部负异常内，故采取限幅拉伸的方法：张玉君地质勘查新方法研究论文集式中 z（x,y）为磁场异常，单位nT,l（x,y）为图像灰度级或称灰阶。每一灰阶代表598 nT，可见经限幅拉伸处理，使磁原始图像保留了更多的弱异常信息，再按（2）式[10]制作立体阴影图（彩版附图9（A））。张玉君地质勘查新方法研究论文集式中 θ为光源方位角；φ为光源高度角；λ为磁面法线矢量与光源矢量的夹角。测区磁场彩色立体阴影图显示出沿着从西北角到东南角对角线两侧成群地存在着高频弱磁异常，其宽度为1～2km，幅度为2～15 nT，在黑白图像上有3～20灰阶的变化。文献[2]认为，钻井套管和井盖可能引起幅度为2～5 nT的异常，但其波长更短（＜300m），也就是说在5×5m网度的图像上，套管和井盖所引起的异常仅表现为孤立点。为了突出有用弱磁异常，用图1中的模板通过褶积功能制作了水平二次导数图像，通过邻域滤波消除孤立异常点。彩版附图9（B）为航磁二次导数、经归一化处理U及油气异常分布叠加图像。青色点群为航磁水平二次导数异常；红色为经归一化处理 U；绿色和蓝色椭圆为已知异常位置；红色椭圆为新发现的有希望异常。10个已知油气显示全部都具有高频弱磁异常，这有力地表明这种仅用数分钟即可完成的磁场二次水平导数图对于提取与油气田有关的高频异常具有重要意义。图1 求导褶积模板高频弱磁异常的磁源埋深对磁场（经镜像对称扩边）进行快速富氏变换，然后对变换结果计算功率谱，从功率谱图像上取径向剖面值，最后做径向功率谱图（图2），为了做图方便直接用图像上的点数作为横坐标。该谱图表明全图由低频功率谱Pl和高频功率谱Ph组成，低频部分在该区主要由前泥盆系变质岩的磁性引起，高频部分主要由浅部后成作用局部矿化引起，用参数法求磁源埋深[14]张玉君地质勘查新方法研究论文集取对数张玉君地质勘查新方法研究论文集图2 径向功率谱如果用半对数坐标成图，式（4）为线性方程，－2H与－2h分别为低频段（1～8点）与高频段（12～50点）的直线斜率，lna与lnb分别为二直线的截距，△r为径向频率间隔。r=2π/（点数×点距）（5）本测区点数为512，点距为5 m，∆r=65×10-6m，用最小二乘法求得二直线的截距和斜率分别为a=1429,Kt=-5595, b=7212,Kh=-9182，张玉君地质勘查新方法研究论文集式中S被导入是因为图像处理所显示的对数功率谱值仅有相对意义，S的单位为m，其值可通过某已知异常源埋深求得。已知低频磁异常系由埋深为－5000至－15000m的前泥盆系变质岩引起，则张玉君地质勘查新方法研究论文集将H代入，S为0285～0855m。张玉君地质勘查新方法研究论文集即高频弱磁源埋深在3～1km之间。航放图像对油气藏的反映航放三元素合成图像[11-12]示于彩版附图9（C），叠加于其上的还有航磁二次水平导数异常。红为K，其含量变化为0～8%；绿为Th，其含量变化为0—39×10-6；蓝为U，其含量变化为0～9×10-6，彩版附图9（D）的上图为其色标；黄色点群为航磁水平二次导数异常。航放图像较之航磁水平二次导数对于油气显示的相关性要复杂得多，虽然并非所有已知油气显示都有清楚的航放异常，但还是有一些含油气构造反映较清楚，如彩版附图9（C）上的5号异常为已知鸭湖含油气穹隆。在双元素图像（彩版附图10（A），彩版附图9（D）的下图为其色标）上鸭湖构造也十分明显；在该图上部为R—Th,G—K,B—K合成图像，箭头所指为已知东台吉乃尔含油气构造；下部为R—Th,G—U,B—U合成图像，箭头所指为已知涩北1、2号含油气构造，它们都表现为周边高中间低的放射性异常。如何增强或突出航放数据中与油气有关的信息，关键在于减少或压制岩性影响，为此，我们采用了主分量分析（在图像处理中称之为KL变换）和Th归一化处理[15]。彩版附图10（B）为航放三元素主分量分析（即KL变换）图像，彩版附图9（D）的上色标亦为其图例；图中a、b、c分别为第1、第2、第3主分量，以KL（1）、KL（2）、KL（3）表示，d为其合成图像。各主分量的组成为：张玉君地质勘查新方法研究论文集式中，Xk、XTh、X。分别代表某像素K、Th、U之含量值。在第1主分量ξ1中以K为主，减少了Th的作用；由于在大部分沙岩、砾岩中K和Th有正相关性，故经此处理的KL（1）大大压制了岩性的影响，突出了K与Th不相关的地质过程的反映，如钾盐的沉积及烟囱效应等。ξ2综合了K、Th、U的变化，反映全区地质环境，对研究油气局部异常作用不明显。ξ3中以U为主，又减弱了Th和K的作用，大大压制了岩性影响，对研究油气局部异常很重要。为减少岩性影响，根据文献[15]，在整理地面或航空放射性测量数据时，常常除以Th含量值，称之为Th归一化处理。在做这种处理时，认为Th能较多地代表砂砾沉积物中岩性的变化，K、U与Th在这些沉积物中有一定正相关关系，除以Th值便减少了或压制了岩性的影响。我们的改进在于在归一化处理之前，先做相关分析，求出相关系数，再按相关系数做Th归一化处理。彩版附图10（C）即为K和U经归一化处理前后的对比图像，左边为原始图像，彩版附图9（D）的上图为其色标；右边为归一化处理结果。表1为已知及新发现异常统计表，共8项特征：K原始图像、K经KL变换、K经Th归一化处理、U原始图像、U经KL变换、U经Th归一化处理、U晕圈及MSS（7、5、4波段）图像上的晕圈反映。已知10处油气显示中有8个（占80%）有明显或较明显的航放异常；其典型特征是：在Th归一化图上中心部位K低、U低和U高值的外围晕圈。Th归一化处理的作用对下述含油气构造的异常十分明显：使V1、V2、V5的K高变为K低，使V1、V5的U高显著降低。故Th归一化处理对于统一异常性质有明显效果。典型异常的剖面图像以鸭湖（V2）和红三旱四号（V1）做为典型的已知异常进行剖面研究。鸭湖异常有以下典型特征（彩版附图10（D））:（1）U在归一化处理前后（彩版附图10（D）右图中之蓝色曲线）异常中心部位都是低值；（2）K在异常中心部位原为高值，经归一化处理后也获低值（彩版附图10（D）右图中之红色曲线）；（3）在异常周围U有明显的高峰和晕圈现象；（4）有明显的高频弱磁异常，且中心处更强（彩版附图10（D）右图中之绿色曲线）。表1 已知及新发现异常统计表红三旱四号的异常特征为：（1）U在归一化处理前异常中心部位为强高值，经归一化处理后显著降低；（2）K在归一化处理前为高值，经归一化处理后为低值；（3）经归一化处理后在异常周围U有明显的高峰，反映出晕圈现象；图3 红三旱四号地面能谱测量剖面（4）也有明显的高频弱磁异常。综合典型异常分析，与油气藏有关的航磁航放异常识别标志为：航磁局部高频弱变化；Th归一化处理后局部低钾、低铀（也有达不到低值）；以及铀在异常周边的高值晕圈现象。我们对红三旱四号异常进行了地面检查，用4道能谱仪进行了剖面测量，总长10km，点距100～200m。所获剖面（图3）与空中测量结果一致，在异常中心部位为K高、U高和Th高，从中心向两侧先降低后又升高。三、对含油气远景评价的讨论本测区内已知含油气构造或点共10处，如表1所列：红三旱四号（V1，气），鸭湖（V2，气），马新高点（V4，油），东台吉乃尔（V5，见工业油气流），南陵丘（V6，油气），驼峰山（V7，油气），涩北1号（V8，油气），涩北2号（V9，油气），涩聂北（V12，气），涩聂东（V3，气）。这10个已知油气显示无一例外都有明显的航磁水平二次导数异常。前8个有明显或较明显的航放异常显示，表现为在经Th归一化处理的K元素图像和U元素图像上为中心部位的低值及U的周边升高的晕圈现象。后两个已知区航放异常不够明显，主要原因是，这两个已知油气显示点处于钾盐类沉积丰富地段，钾引起的放射性变化强于与碳氢化合物有关的弱放射性异常，后者未能明显表现出来。故本测区已知油气显示构造或点与航磁水平二次导数异常的吻合率为100%，与航放K、U异常的吻合率为80%。李芦玲等李芦玲等柴达木盆地中部地区航空磁测详查结果报告地质矿产部航空物探遥感中心，1985。曾用传统的方法对本测区含油气远景进行了评价，圈出8个最有希望的含油气区，图像处理也完全证实了这一评价。经图像处理，我们认为还应补充提出以下4个有希望的含油气异常区：（1）依克雅乌汝背斜构造（V3），位于鸭湖含油气背斜构造东30km，台吉乃尔含油气背斜北15km。其航磁航放、卫片异常均很典型，且与鸭湖异常极为相似，异常特征依据充分，可做为新发现的最有希望的含油气异常点。（2）南陵丘西异常（V10），位于南陵丘已知含油气构造西约5km。有局部高频低幅磁异常、低钾、低铀及较明显可辨之周边闭合的高值U晕圈。可做为新发现有希望的含油气异常。（3）洪上异常（V11），位于台吉乃尔洪积扇上部。除有局部高频弱磁异常外，有低U及隐约可见的周边U高值晕圈，可做为新发现较有希望的含油气异常。（4）南八仙异常（V14），位于仙西南高点北缘。有局部高频弱磁、低钾和低铀异常。也可做为新发现之较有希望的含油气异常。如前所述彩版附图9（B）上展示了l0处已知（蓝、绿、黄色椭圆）和4处新发现油气异常（红色椭圆）。并与归一化处理后之U（红）以及航磁水平二次导数异常（青点）叠加在一起。上述远景评价于1990年完成，根据能源部1991年在柴达木中部地区进行的勘探成果，对比如下，上述4个含油气远景区中的3个见到了油气显示：即依克雅乌汝背斜构造（V3）于“伊中1”井见油气显示，洪上异常于“台中1”井见工业气流，南八仙异常于“仙3”井见油气显示。本工作所预测的4个含油气远景区中有3个见到了油气显示，说明航磁、航放等综合参数通过图像处理所做油气田异常特征的研究是有效的。常规方法没有发现的异常，却由本方法预测，并已得到验证，说明该方法的价值。参加野外验证的还有水恩海、史殿林、郭毅，屏幕图像均由杨星虹拍摄，能源部石油勘探开发科学研究院遥感地质所叶和飞给予了支持。在此表示诚挚的谢意。参考文献[1]Wold,R郭武林整理利用航磁和放射性测量直接寻找油气藏地质科技动态，1984,（24）:21～[2]管志宁利用航磁探测与化学剩磁有关的油气矿床国外地质勘探技术，1985,（7）:18～[3]王家林，万明浩，金国英磁测找油的进一步试验及有关问题的研究石油地球物理勘探，1985,（20）:397～[4]郭玉琨编译苏联地球物理直接找油方法概况石油物探译丛，1985,（1）:90～[5]王秀文摘译石油和天然气放射勘探的新远景地质科技动态，1984,（18）:17～[6]赵改善编译直接勘探油气田的放射性方法石油物探译丛，1987,（6）:95～[7]李淑仪核技术在石油天然气和地热勘探中的应用国外地质勘探技术，1986,（4）:22～[8]王锡田编油气勘查中用遥感方法检测烃类微渗漏地质科技动态，1987,（11）:13～16[9]蔡振京高精度大比例尺（1∶50000）航空物探在油气藏勘探中的应用物探与化探，1989（13）:401～[10]Teskey,D,Broome,JComputer programs for production of shaded relief and stereo shaded relief GSC 84—1 B,375—389,1984[11]Zhang YApplication of image processing techniques to airborne radiometric Abstract of the Second Symposium on Exploration Geophysics,Xian,China,525—526,1986[12]Zhang YDigital image processing of airborne radiometric and magnetic data from central chaidam An Overview of Exploration Geophysics in China,517—535,1989[13]Donovan,TJ,Forgey, RL, Roberts, AAAeromagnetic detection of diagenetic magnetite over oil Bulletin ofthe American Association ofPetroleum Geologists,63,245—248,1979[14]Donovan,TJ,Brien, DPBryan,JG, Cunningham, KINear surface magnetic indications of buried Aeromagnetic detection and separation of spurious signals, IGARSS' 87 Remote Sensing: understanding the earth as a system,219—232,1987[15]Alton VGRadiometrics—a practical exploration Weotern Oil Reporter'March,173,A STUDY OF AERO-MAGNETIC AND AERORADIOMETRIC ANOMALIES IN SOME OIL-GAS FIELDS BY IMAGE PROCESSING TECHNIQUESZhang Yu jun（Center of Aero-Geophysics and Remote Sensing, Ministry of Geology and Mineral Resources , Beijing 100083）AbstractThis work represents the results of study of the following problemsby image processing techniques: the quick extraction of weak magnetic anomalies with high frequency, the evaluation of depth for highfrequence magnetic anomalies,the normalization of radiometric anomalies, the analysis of typical anomalous profiles and the prediction of perspective anomalies to find oil and Four new perspective anomalies for oil gas were The study had been finished at the beginning of Comparing with the results of exploration in 1991, we have been encouraged by the fact that not only on 3 out of the 4 anomalieas the oil and gas werefound, but also on one of them the gas was Key words Oil-gas filed,The chimney effect,Image processing, The second holizontal gradient of magneticfield, The power spectrum, The Th normalization for aeroradiometrics, The perspective evaluation for oil-原载《地球物理学报》，1994,V37,N1。
张玉君（地质矿产部航空物探遥感中心研究所北京）摘要：本文介绍一种独特的航放数据图像复原方法。该方法的主要技术关键是：提出航放数据图像复原原理和理论基础；建立航放数据图像复原处理流程；制定重建数据网格文件的途径；进行航放数据图像复原效果及误差评价。关键词：航放数据，大气本底，图像处理，图像复原技术。一、引言自20世纪70年代初期方柱形NaI晶体进入机载综合航空站以来，航放测量的灵敏度和有效性显著提高，地质及地球物理勘探界对于航空放射性测量的需求产生了根本变化。在近20年的应用实际中，大气氡本底（简称大气本底）的改正很难准确，它始终是困扰该方法应用效果的主要难点；其后果是在图面上造成条带现象，从而严重地影响着图件的可用性及方法的效果。其原因[1]可概括为：空中所测放射性不仅仅来源于地下，而且受飞机硬件环境、宇宙射线、大气中氡及其子体的影响。后者称之为大气本底干扰，它又受气候、风力、风向、温度、季节及一天中何时测量等因素的影响。大气本底干扰的主要表现形式是架次与架次间本底水平不同。受干扰最大的是铀道，钾道次之，钍道和总道虽较小，但也不可忽视，（见彩版附图7中的图3、图4）。由于这种噪声的存在，来源于地质体的信息常常淹没于噪声之中。图3a（彩版附图7）为哈密土墩测区K（红）、Th（绿）、U（蓝）三元素复原图像，图3b为该测区航放原始数据合成图像，图4a（彩版附图7）表示各架次早、晚校准读数，图4b为总道原始数据图像。条带噪声的存在，可以形象地比喻为挂在有用信息图像前面的彩条窗帘，条带的严重性使得该工区原始航放数据无法绘制等值线图。航放图面条带问题是一个“世界性”的问题[2]。解决得较好的是加拿大，靠星罗棋布的湖泊河流等水域上空测量结果改正本底，取得较好的效果，而且不使用向上探头[1]。美国Geometrics公司及其他航空物探公司则是靠向上探头测量，作为本底改正的依据[3]。1986年Grasty[4]提出当测区内没有湖泊时，可用测线上无异常区的平均值代替本底。本文介绍的方法与国际上已采用的各种方法全然不同，该方法在数字图像处理学中可称之为航放图像复原技术。图像复原技术的主要目的，是要改善给定的图像。复原是一个过程，它试图利用蜕化现象的某种先验知识，把已经蜕化了的图像加以重建或恢复。因此，复原技术是把蜕化模型化，并运用相反的过程在某种程度上恢复原来的图像。Cannon博士[5]研究了一种图像复原技术或称图案去除技术，它适用于：从规则图案（如纺织品）上提取指纹图形，改善散焦图像，消除卫片图像探测器与探测器间的噪声，使在曝光过程中相机或物体平移造成模糊的图像清晰化等。Srinivasan也报道了这类研究[6]。张玉君等研究了深海锰结核照片光照不均匀等蜕化现象的图像复原问题[7]。航放数据图像复原处理是数字图像复原技术在地学界成功应用的又一实例，但航放数据图像所存在的蜕化问题与上述各例均不相同。该方法研究成功后，曾在6个测区得到验证。二、航放数据图像复原技术原理和理论基础航放所测到的是一幅蜕化了的图像G（x,y）它可视为由真实图像F（x,y）与干扰图像η（x,y）叠加而成，简化了的蜕化过程，见图1。航放图像蜕化现象的先验知识来自对航放测量过程及原始图像的分析。在测量过程中，来自地质体的有用信息是不随时间为转移的。而干扰在本质上是随时间变化的，但在图像上干扰已变为（x,y）的函数，因为：张玉君地质勘查新方法研究论文集图1 航放数据图像蜕化示意图η的变化可分为架次之间的跳变及架次之内的渐变，见图4（彩版附图7），在每一测线上此干扰大致为一个常数，如果将x（即图像上的列）表示垂直测线的方向，则η（x,y）简化为η（x），则有张玉君地质勘查新方法研究论文集航放图像复原的目的，就在于设法近似地求出η（x），从而近似地得到F（x,y）。为此，沿测线方向对原始图像进行多次单列多行窄长窗口褶积：张玉君地质勘查新方法研究论文集式中W为褶积模板，是由加权因子组成的矩阵。褶积过程是一种线性运算，其算子H不随空间变化。因为算子为线性的，则两个输入之和的响应等于两个响应之和。张玉君地质勘查新方法研究论文集由于假定了η仅与x有关，又由于褶积窗口为单列，则有：张玉君地质勘查新方法研究论文集现在分析HF（x,y）的性质，由于沿y方向的多次滑动平均，局部异常“淹没”于近区域特征之中，这种近区域特征表现为沿测线方向的低缓变化；如果用f（x,y）表示局部异常，用L（x,y）表示近区域场，于是：张玉君地质勘查新方法研究论文集再经如下处理张玉君地质勘查新方法研究论文集由（9）式可见，从原始图像中减去噪声图像后，所得复原图像f（x,y），它从局部异常角度是接近真实图像的，误差取决于所减去的“近区域背景值”在测线方向起伏变化的幅度。三、航放数据图像复原处理过程航放数据图像复原技术的研究是以多元统计为理论基础，以图像处理为工具来完成的，并体现了图像处理快速直观的特点，其流程如图2所示。图2 航放数据图像复原处理流程此方法假定航放噪声本底沿测线方向不变或呈线性变化。通过沿测线方向的多次滑动平均，使局部异常逐渐淹没在噪声本底之中，得出一幅与噪声本底线性相关的噪声图像。噪声图像尚需进行边缘影响补偿；对去除噪声后的图像，经中值滤波和空间变量反差增强，达到最终复原的效果。这一复原过程集中表示于图2的左半部。图2的右半部为数据网格文件的重建过程，它是实际应用所必不可少的。经过分类分区，求得各类别在复原前后的均值向量，经最小二乘拟合求出复原图像的元素含量或计数率值，重新建立为在主计算机上绘制等值线图用的网格文件。本研究曾试验通过沿测线方向取平均值做为噪声水平，结果不及上述方法理想。四、效果及误差评价航放数据图像复原的效果（1）图面直观效果的改善。可以形象地说，航放图像复原好比揭去一层条带窗帘，使原来透过此窗帘隐约可见的图像显示出了真面目，见图3a（彩版附图7）。图面直观效果的改善还表现在由于定位问题所引起的岩体边界上的锯齿状噪声得以消除，见图5（彩版附图7）。图5为总道对比图像，5a为原始数据，5b表示噪声图像，5c为去噪声后的图像，5d为复原图像。（2）用复原数据所做等值线图真实可信。以哈密土墩测区为例，原始资料由于条带干扰，在主计算机上，钾、钍、铀道都无法绘等值线图，仅提供了平剖图；只有总道提供了等值线图，但仍可看到条带的影响。经图像复原、重建网格文件，反馈回主计算机后绘制了TC、K、Th和U等值线图，现以K道经复原后数据等值线图为例示于图6（彩版附图7），与地质图对比，表明异常和地质体对应良好，各类岩性的放射性趋势也都吻合，证实这些等值线图的可靠性。利用复原图像所做分类图也证实了这一点，见图7（彩版附图7），图7中数字分别为：①超基性岩；②基性岩；③花岗岩；④闪长岩；⑤变质岩；⑥混合岩；⑦第四纪沉积；⑧第三、四纪沉积；⑨第三纪沉积。（3）有用信息增加。本研究利用多元统计的方法，对航放图像复原的效果给出定量评价。可用一幅图像有用信息构成的变异值的大小来对它做定量评价。为此应计算全图面总变异对于一个象素的平均值，即平均变异值。用C、C´和C＂分别表示原始图像中有用信息平均变异值、原始图像中干扰信息平均变异值和最终复原图像有用信息平均变异值。统计时以G´（x,y）近似代表η（x）；以[G（x,y）-G´（x,y）]近似代表F（x,y）；以P（x,y）表示最终复原图像，并假定它已无干扰存在。张玉君地质勘查新方法研究论文集式中，字母上加“—”表示平均值；M、N为图像的行、列数。表1为哈密土墩测区航放数据图像按上述各式所做定量评价的统计结果。表1从表1可见，K、Th、U、TC经图像复原后，有用信息都有十分显著的增长；就此工区而论，TC和K原始图像相对质量较好，Th和U较差。复原图像的准确度及误差评价复原图像的主要误差来源是“近区域背景值”L（x,y），它是在多次滑动平均时形成的。通过对干扰图像剖面数据的统计，得到以下准确度评价：K±16%（绝对含量）；Th+1 pp m；U±15 ppm；TC±6计数。五、结论（1）本文介绍的方法是在国内外首次提出的方法独特的航放数据图像的复原技术，并在多个工区验证了其可靠性和实用性。（2）本技术可以基本上消除由于大气本底及阈值变动所造成的图面条带现象，基本复原航放图像的真面目，为进一步图像处理（诸如：求导、增强、分类、逻辑运算等）做了准备，因此本技术也是一种快速预处理方法。（3）本方法改善了由于飞行往返定位位移所引起的某些地质体边缘呈锯齿状的图像噪声问题。（4）本研究建立了“有用信息平均变异值”做为定量评价航放数据图像复原效果的尺度。还讨论了图像复原做为一种预处理过程，对于元素含量值可能导入的绝对误差或称为方法的准确度。参考文献[1]Grasty, RL, Gamma ray spectrometric methods in uranium exploration—Thcory and operational procedures, Geophysics and Geochemistry in the Search for Metallic Ores,GSC,Ottawa,147-162,[2]Creen, AA,Leveling airborne gamma-radiation data using between-channel correlation information,Geophysics,52,1557-1562,[3]Foote, RS, Improvement in airborne gamma-radiation data analysis by removal of environmental and pedologic radiation changes, in S the Use of Nuclear Techniques in Prospecting and Developmcnt of Mineral Resources: IAtomic Energy Agency M,Buenos Aires,187-196,[4]Grasty,RL,Automated system for computing on-line atmospheric backgrounds,GSC paper,1-52,[5]Cannon,M,Lehar,Aand Preston,F,Background pattern removal by power spectral filtering,Applied Optics,22,777-779,[6]Srinivasan,R,Software image restoration techniques,Digital Design, 16,4,27-34, [7]张玉君，史鉴文深海多金属结核照片的图像复原和图像处理技术研究物探与化探，1989,（13）:435～致谢 林振民同志对本文提出了宝贵的意见，史鉴文同志参加了重复工区试验，张志民和谢欣同志分别编制了网格文件转换和最小二乘拟合程序，杨星虹同志拍摄了屏幕图片，水恩海同志搜集了试验工区校准资料，在此一并致谢。A STUDY ON IMAGE RESTORATION TECHNIQUES FOR AERORADIOMETRIC DATAZhang Yu jun（Research Institute, Center of Aero-Gcophysics and Remote Sensing,Ministry of Geology and Mineral Resources, Beijing）AbstractThis paper represents a specific methodfor restoration of images of airborne radiometric The main technical keys involved in this study are；the advancementof the principles and theory；the establishment of the flow-diagram for processing；the formulation of the means for reestablishment of the gridded data file；the evaluation of the restoration results and the errors, involved by the restoration Key words Aeroradiometric data, Atmospheric background, Image processing, Image restoration 原载《地球物理学报》，1990,V33,N4。
张玉君（地质矿产部航空物探遥感中心，北京100083）摘要：用图像处理技术对下述与油气田异常有关的问题进行了研究：快速提取高频弱磁异常，估算高频磁异常源的埋深，放射性异常性质的归一化处理，典型异常剖面的分析，以及油气远景评价，提出了4个有希望含油气的异常。经验证，这4个异常中有3个见到油气显示，其中1个为工业气流。关键词：油气田，烟囱效应，图像处理，航放Th归一化处理，油气远景评价。一、前言近10年，非地震物探方法寻找油气田的工作越来越受到人们的重视[1-8]，原因之一是地震勘探生产费用日益昂贵；另一方面，近10年航磁、航放测量技术发生了重大变化，航空光泵磁力仪的灵敏度从1nT提高到1～01nT，大体积方柱形NaI（Tl）晶体和数字电路被引入航空放射性测量仪，不仅使其灵敏度提高了10多倍，而且使强度测量被4道（K、Th、U及总道）或多道（256或512道）能谱测量所取代。这一切促使物探工作者在油气勘查中更加注意应用具有大面积低消耗特性的航空物探技术。柴达木中部地区的高灵敏度航磁、航空伽马能谱测量飞行线距1km，飞行高度约90m；安装在双水獭型飞机上的航放、航磁综合站包括GAD-6型4道能谱仪、MAP-4型质子旋进磁力仪、磁带收录机及多普勒雷达、航空摄影两种导航定位设备。能谱仪使用了3箱NaI（T1）晶体，总体积为48L，磁力仪的灵敏度为5nT。放射性伽马测量普查油气已有50多年的历史，实践表明，在油气田上方有偏低的放射性异常，并在油气构造边缘存在比本低高的异常峰值。1979年美国Donovan TJ等人在Cement油田发现了高频磁异常[13],10多年来人们公认了这种与油气田相关的高频磁异常[9]的存在。尽管对产生这些异常的机制存在不少争议，但最被接受的假说是烟囱效应，通常认为几乎所有的储油构造都从深部排出水份，烟囱效应将碳氢化合物和其他有机衍生物通过喷溢、溶液散逸、气体扩散等三种渗漏作用带到地表；地层中的液体和气体在从油、气层向上运移过程中，由于后成作用形成磁铁矿并沉淀放射性元素，从而形成可被航测发现的航磁航放异常，并作为判断地下赋存有碳氢化合物的标志。提出用图像处理技术来研究这一问题，基于两点考虑：航空物探的大面积宏观数据通过图像处理可以快速获得直观结果；图像系统具备既适合位场数据（航磁）又适合统计学性质数据（航放）的处理功能。本文所做基础图像网度为5×5m2。除航磁、航放外还利用了陆地卫星影像 MSS（7、5、4波段）数据。二、异常特征研究高频弱磁异常的快速提取柴达木盆地是我国重要的含油气盆地之一，在测区内已有多处井下油气显示，测区内人烟罕见，人文干扰较少，为研究油气田后成磁铁矿引起的高频弱磁异常特征的理想试验场地。测区磁场动态范围为0～533 nT，在0～127 nT的低端仅有12个像素，占全图面的5/100000，并集中在东北部的局部负异常内，故采取限幅拉伸的方法：张玉君地质勘查新方法研究论文集式中 z（x,y）为磁场异常，单位nT,l（x,y）为图像灰度级或称灰阶。每一灰阶代表598 nT，可见经限幅拉伸处理，使磁原始图像保留了更多的弱异常信息，再按（2）式[10]制作立体阴影图（彩版附图9（A））。张玉君地质勘查新方法研究论文集式中 θ为光源方位角；φ为光源高度角；λ为磁面法线矢量与光源矢量的夹角。测区磁场彩色立体阴影图显示出沿着从西北角到东南角对角线两侧成群地存在着高频弱磁异常，其宽度为1～2km，幅度为2～15 nT，在黑白图像上有3～20灰阶的变化。文献[2]认为，钻井套管和井盖可能引起幅度为2～5 nT的异常，但其波长更短（＜300m），也就是说在5×5m网度的图像上，套管和井盖所引起的异常仅表现为孤立点。为了突出有用弱磁异常，用图1中的模板通过褶积功能制作了水平二次导数图像，通过邻域滤波消除孤立异常点。彩版附图9（B）为航磁二次导数、经归一化处理U及油气异常分布叠加图像。青色点群为航磁水平二次导数异常；红色为经归一化处理 U；绿色和蓝色椭圆为已知异常位置；红色椭圆为新发现的有希望异常。10个已知油气显示全部都具有高频弱磁异常，这有力地表明这种仅用数分钟即可完成的磁场二次水平导数图对于提取与油气田有关的高频异常具有重要意义。图1 求导褶积模板高频弱磁异常的磁源埋深对磁场（经镜像对称扩边）进行快速富氏变换，然后对变换结果计算功率谱，从功率谱图像上取径向剖面值，最后做径向功率谱图（图2），为了做图方便直接用图像上的点数作为横坐标。该谱图表明全图由低频功率谱Pl和高频功率谱Ph组成，低频部分在该区主要由前泥盆系变质岩的磁性引起，高频部分主要由浅部后成作用局部矿化引起，用参数法求磁源埋深[14]张玉君地质勘查新方法研究论文集取对数张玉君地质勘查新方法研究论文集图2 径向功率谱如果用半对数坐标成图，式（4）为线性方程，－2H与－2h分别为低频段（1～8点）与高频段（12～50点）的直线斜率，lna与lnb分别为二直线的截距，△r为径向频率间隔。r=2π/（点数×点距）（5）本测区点数为512，点距为5 m，∆r=65×10-6m，用最小二乘法求得二直线的截距和斜率分别为a=1429,Kt=-5595, b=7212,Kh=-9182，张玉君地质勘查新方法研究论文集式中S被导入是因为图像处理所显示的对数功率谱值仅有相对意义，S的单位为m，其值可通过某已知异常源埋深求得。已知低频磁异常系由埋深为－5000至－15000m的前泥盆系变质岩引起，则张玉君地质勘查新方法研究论文集将H代入，S为0285～0855m。张玉君地质勘查新方法研究论文集即高频弱磁源埋深在3～1km之间。航放图像对油气藏的反映航放三元素合成图像[11-12]示于彩版附图9（C），叠加于其上的还有航磁二次水平导数异常。红为K，其含量变化为0～8%；绿为Th，其含量变化为0—39×10-6；蓝为U，其含量变化为0～9×10-6，彩版附图9（D）的上图为其色标；黄色点群为航磁水平二次导数异常。航放图像较之航磁水平二次导数对于油气显示的相关性要复杂得多，虽然并非所有已知油气显示都有清楚的航放异常，但还是有一些含油气构造反映较清楚，如彩版附图9（C）上的5号异常为已知鸭湖含油气穹隆。在双元素图像（彩版附图10（A），彩版附图9（D）的下图为其色标）上鸭湖构造也十分明显；在该图上部为R—Th,G—K,B—K合成图像，箭头所指为已知东台吉乃尔含油气构造；下部为R—Th,G—U,B—U合成图像，箭头所指为已知涩北1、2号含油气构造，它们都表现为周边高中间低的放射性异常。如何增强或突出航放数据中与油气有关的信息，关键在于减少或压制岩性影响，为此，我们采用了主分量分析（在图像处理中称之为KL变换）和Th归一化处理[15]。彩版附图10（B）为航放三元素主分量分析（即KL变换）图像，彩版附图9（D）的上色标亦为其图例；图中a、b、c分别为第1、第2、第3主分量，以KL（1）、KL（2）、KL（3）表示，d为其合成图像。各主分量的组成为：张玉君地质勘查新方法研究论文集式中，Xk、XTh、X。分别代表某像素K、Th、U之含量值。在第1主分量ξ1中以K为主，减少了Th的作用；由于在大部分沙岩、砾岩中K和Th有正相关性，故经此处理的KL（1）大大压制了岩性的影响，突出了K与Th不相关的地质过程的反映，如钾盐的沉积及烟囱效应等。ξ2综合了K、Th、U的变化，反映全区地质环境，对研究油气局部异常作用不明显。ξ3中以U为主，又减弱了Th和K的作用，大大压制了岩性影响，对研究油气局部异常很重要。为减少岩性影响，根据文献[15]，在整理地面或航空放射性测量数据时，常常除以Th含量值，称之为Th归一化处理。在做这种处理时，认为Th能较多地代表砂砾沉积物中岩性的变化，K、U与Th在这些沉积物中有一定正相关关系，除以Th值便减少了或压制了岩性的影响。我们的改进在于在归一化处理之前，先做相关分析，求出相关系数，再按相关系数做Th归一化处理。彩版附图10（C）即为K和U经归一化处理前后的对比图像，左边为原始图像，彩版附图9（D）的上图为其色标；右边为归一化处理结果。表1为已知及新发现异常统计表，共8项特征：K原始图像、K经KL变换、K经Th归一化处理、U原始图像、U经KL变换、U经Th归一化处理、U晕圈及MSS（7、5、4波段）图像上的晕圈反映。已知10处油气显示中有8个（占80%）有明显或较明显的航放异常；其典型特征是：在Th归一化图上中心部位K低、U低和U高值的外围晕圈。Th归一化处理的作用对下述含油气构造的异常十分明显：使V1、V2、V5的K高变为K低，使V1、V5的U高显著降低。故Th归一化处理对于统一异常性质有明显效果。典型异常的剖面图像以鸭湖（V2）和红三旱四号（V1）做为典型的已知异常进行剖面研究。鸭湖异常有以下典型特征（彩版附图10（D））:（1）U在归一化处理前后（彩版附图10（D）右图中之蓝色曲线）异常中心部位都是低值；（2）K在异常中心部位原为高值，经归一化处理后也获低值（彩版附图10（D）右图中之红色曲线）；（3）在异常周围U有明显的高峰和晕圈现象；（4）有明显的高频弱磁异常，且中心处更强（彩版附图10（D）右图中之绿色曲线）。表1 已知及新发现异常统计表红三旱四号的异常特征为：（1）U在归一化处理前异常中心部位为强高值，经归一化处理后显著降低；（2）K在归一化处理前为高值，经归一化处理后为低值；（3）经归一化处理后在异常周围U有明显的高峰，反映出晕圈现象；图3 红三旱四号地面能谱测量剖面（4）也有明显的高频弱磁异常。综合典型异常分析，与油气藏有关的航磁航放异常识别标志为：航磁局部高频弱变化；Th归一化处理后局部低钾、低铀（也有达不到低值）；以及铀在异常周边的高值晕圈现象。我们对红三旱四号异常进行了地面检查，用4道能谱仪进行了剖面测量，总长10km，点距100～200m。所获剖面（图3）与空中测量结果一致，在异常中心部位为K高、U高和Th高，从中心向两侧先降低后又升高。三、对含油气远景评价的讨论本测区内已知含油气构造或点共10处，如表1所列：红三旱四号（V1，气），鸭湖（V2，气），马新高点（V4，油），东台吉乃尔（V5，见工业油气流），南陵丘（V6，油气），驼峰山（V7，油气），涩北1号（V8，油气），涩北2号（V9，油气），涩聂北（V12，气），涩聂东（V3，气）。这10个已知油气显示无一例外都有明显的航磁水平二次导数异常。前8个有明显或较明显的航放异常显示，表现为在经Th归一化处理的K元素图像和U元素图像上为中心部位的低值及U的周边升高的晕圈现象。后两个已知区航放异常不够明显，主要原因是，这两个已知油气显示点处于钾盐类沉积丰富地段，钾引起的放射性变化强于与碳氢化合物有关的弱放射性异常，后者未能明显表现出来。故本测区已知油气显示构造或点与航磁水平二次导数异常的吻合率为100%，与航放K、U异常的吻合率为80%。李芦玲等李芦玲等柴达木盆地中部地区航空磁测详查结果报告地质矿产部航空物探遥感中心，1985。曾用传统的方法对本测区含油气远景进行了评价，圈出8个最有希望的含油气区，图像处理也完全证实了这一评价。经图像处理，我们认为还应补充提出以下4个有希望的含油气异常区：（1）依克雅乌汝背斜构造（V3），位于鸭湖含油气背斜构造东30km，台吉乃尔含油气背斜北15km。其航磁航放、卫片异常均很典型，且与鸭湖异常极为相似，异常特征依据充分，可做为新发现的最有希望的含油气异常点。（2）南陵丘西异常（V10），位于南陵丘已知含油气构造西约5km。有局部高频低幅磁异常、低钾、低铀及较明显可辨之周边闭合的高值U晕圈。可做为新发现有希望的含油气异常。（3）洪上异常（V11），位于台吉乃尔洪积扇上部。除有局部高频弱磁异常外，有低U及隐约可见的周边U高值晕圈，可做为新发现较有希望的含油气异常。（4）南八仙异常（V14），位于仙西南高点北缘。有局部高频弱磁、低钾和低铀异常。也可做为新发现之较有希望的含油气异常。如前所述彩版附图9（B）上展示了l0处已知（蓝、绿、黄色椭圆）和4处新发现油气异常（红色椭圆）。并与归一化处理后之U（红）以及航磁水平二次导数异常（青点）叠加在一起。上述远景评价于1990年完成，根据能源部1991年在柴达木中部地区进行的勘探成果，对比如下，上述4个含油气远景区中的3个见到了油气显示：即依克雅乌汝背斜构造（V3）于“伊中1”井见油气显示，洪上异常于“台中1”井见工业气流，南八仙异常于“仙3”井见油气显示。本工作所预测的4个含油气远景区中有3个见到了油气显示，说明航磁、航放等综合参数通过图像处理所做油气田异常特征的研究是有效的。常规方法没有发现的异常，却由本方法预测，并已得到验证，说明该方法的价值。参加野外验证的还有水恩海、史殿林、郭毅，屏幕图像均由杨星虹拍摄，能源部石油勘探开发科学研究院遥感地质所叶和飞给予了支持。在此表示诚挚的谢意。参考文献[1]Wold,R郭武林整理利用航磁和放射性测量直接寻找油气藏地质科技动态，1984,（24）:21～[2]管志宁利用航磁探测与化学剩磁有关的油气矿床国外地质勘探技术，1985,（7）:18～[3]王家林，万明浩，金国英磁测找油的进一步试验及有关问题的研究石油地球物理勘探，1985,（20）:397～[4]郭玉琨编译苏联地球物理直接找油方法概况石油物探译丛，1985,（1）:90～[5]王秀文摘译石油和天然气放射勘探的新远景地质科技动态，1984,（18）:17～[6]赵改善编译直接勘探油气田的放射性方法石油物探译丛，1987,（6）:95～[7]李淑仪核技术在石油天然气和地热勘探中的应用国外地质勘探技术，1986,（4）:22～[8]王锡田编油气勘查中用遥感方法检测烃类微渗漏地质科技动态，1987,（11）:13～16[9]蔡振京高精度大比例尺（1∶50000）航空物探在油气藏勘探中的应用物探与化探，1989（13）:401～[10]Teskey,D,Broome,JComputer programs for production of shaded relief and stereo shaded relief GSC 84—1 B,375—389,1984[11]Zhang YApplication of image processing techniques to airborne radiometric Abstract of the Second Symposium on Exploration Geophysics,Xian,China,525—526,1986[12]Zhang YDigital image processing of airborne radiometric and magnetic data from central chaidam An Overview of Exploration Geophysics in China,517—535,1989[13]Donovan,TJ,Forgey, RL, Roberts, AAAeromagnetic detection of diagenetic magnetite over oil Bulletin ofthe American Association ofPetroleum Geologists,63,245—248,1979[14]Donovan,TJ,Brien, DPBryan,JG, Cunningham, KINear surface magnetic indications of buried Aeromagnetic detection and separation of spurious signals, IGARSS' 87 Remote Sensing: understanding the earth as a system,219—232,1987[15]Alton VGRadiometrics—a practical exploration Weotern Oil Reporter'March,173,A STUDY OF AERO-MAGNETIC AND AERORADIOMETRIC ANOMALIES IN SOME OIL-GAS FIELDS BY IMAGE PROCESSING TECHNIQUESZhang Yu jun（Center of Aero-Geophysics and Remote Sensing, Ministry of Geology and Mineral Resources , Beijing 100083）AbstractThis work represents the results of study of the following problemsby image processing techniques: the quick extraction of weak magnetic anomalies with high frequency, the evaluation of depth for highfrequence magnetic anomalies,the normalization of radiometric anomalies, the analysis of typical anomalous profiles and the prediction of perspective anomalies to find oil and Four new perspective anomalies for oil gas were The study had been finished at the beginning of Comparing with the results of exploration in 1991, we have been encouraged by the fact that not only on 3 out of the 4 anomalieas the oil and gas werefound, but also on one of them the gas was Key words Oil-gas filed,The chimney effect,Image processing, The second holizontal gradient of magneticfield, The power spectrum, The Th normalization for aeroradiometrics, The perspective evaluation for oil-原载《地球物理学报》，1994,V37,N1。
淮 阴 工 学 院毕业设计（论文）报告题目 级 专业学 生 姓 名 指 导 教 师 起 讫 日 期 设 计 地 点 年 月 日（空2行） 目 次 （4号黑体，居中）1 引言（或绪论）（作为正文第1章，小4号宋体，行距18磅，下同）………… 12 ×××××× （正文第2章）…………………………………………………… Y1 ×××××× （正文第2章第1条） ………………………………………… Y2 ×××××× （正文第2章第2条） ………………………………………… YX ×××××× （正文第2章第X条） ………………………………………… Y3 ×××××（正文第3章） ……………………………………………… Y………………………………………（略）X ×××××（正文第X章）……………………………………………………… Y结论 …………………………………………………………………………………… Y致谢 …………………………………………………………………………………… Y参考文献………………………………………………………………………………Y附录A ××××（必要时） ………………………………………………………… Y附录B ××××（必要时） ………………………………………………………… Y图1 ×××××（必要时）………………………………………………………… Y图2 ×××××（必要时）………………………………………………………… Y表1 ×××××（必要时）………………………………………………………… Y表2 ×××××（必要时）………………………………………………………… Y注： 目录中的内容一般列出“章”、“条”二级标题即可；2．X、Y表示具体的阿拉伯数字；3．页眉中的页码用罗马数字（Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ……）表示。1 引言（或绪论）（可作为正文第1章标题，用小3号黑体，加粗，并留出上下间距为：段前5行，段后5行）×××××××××（小4号宋体，5倍行距）×××××××××××××××××××××………1．1 ××××××（作为正文2级标题，用4号黑体，加粗）×××××××××（小4号宋体）××××××…………1 ××××（作为正文3级标题，用小4号黑体，不加粗）×××××××××（小4号宋体）×××××××××××××××××××××××××××………2 ×××××××（作为正文第2章标题，用小3号黑体，加粗，并留出上下间距为：段前5行，段后5行）×××××××××（小4号宋体）×××××××××××××××××××××××××××××××××××………注：1．正文中表格与插图的字体一律用5号宋体；2．正文各页的格式请以此页为标准复制，页眉中的页码用阿拉伯数字表示（本文档的页码已设置成自动格式）；3．为保证打印效果，学生在打印前，请将全文字体的颜色统一设置成黑色。 （空2行） 结 论（小3号黑体，居中）×××××××××（小4号宋体，5倍行距）×××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××…………（空2行） 致 谢（小3号黑体，居中）×××××××××（小4号宋体，5倍行距）×××××××××××××××××××××…………（空2行） 参 考 文 献（小3号黑体，居中）1 ×××××××（小4号宋体，行距18磅）×××××2 ××××××××××××××××××××××××××××××××××××××××××3 ××××××××××××××××××××××…………例如：〔1〕刘谋佶，吕志咏，丘成昊，等边条翼与旋涡分离流[M)。北京：北京航空学院出版社，1988．24—27．〔2〕傅惠民．二项分布参数整体推断方法[J]．航空学报，2000，21(2)：155～158．〔3〕朱刚．新型流体有限元法及叶轮机械正反混合问题[D]．北京：清华大学，1996．〔4〕辛希孟．信息技术与信息服务国际研讨会论文集：A集[C]．北京:中国社会科学出版社，1994．〔5〕陈永康，李素循，李玉林．高超声速流绕双椭球的实验研究IA]．见:北京空气动力研究所编．第九届高超声速气动力会议论文集[C]．北京:北京空气动力研究所，1997：9—14．〔6〕孔祥福．FD-09风洞带地面板条件下的流场校测报告[R]．北京空气动力研究所技术报告B67-270，北京：北京空气动力研究所，1989．〔7〕黎志华，黎志军．反馈声抵消器[P]．中国专利：ZL85100748，1986—09—24．…………毕 业 论 文 中 文 摘 要（空1行） 摘要： ××××××××××××××××（小4号宋体，5倍行距）×××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××。（要求200~300字）（空2行）关键词 ××× ××× ××× ×××（小4号宋体）（小4号黑体）毕业设计说明书（论文）外文摘要Title ××××××××（4号宋体） Abstract×××××××××（小4号宋体，5倍行距，第一个字应顶格写）××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××（空2行）Keywords ××× ××× ××× ×××（小4号宋体）（小4号黑体）
1、论文格式的论文题目：（下附署名）要求准确、简练、醒目、新颖。 2、论文格式的目录 目录是论文中主要段落的简表。（短篇论文不必列目录） 3、论文格式的内容提要： 是文章主要内容的摘录，要求短、精、完整。字数少可几十字，多不超过三百字为宜。 4、论文格式的关键词或主题词 关键词是从论文的题名、提要和正文中选取出来的，是对表述论文的中心内容有实质意义的词汇。关键词是用作计算机系统标引论文内容特征的词语，便于信息系统汇集，以供读者检索。每篇论文一般选取3-8个词汇作为关键词，另起一行，排在“提要”的左下方。 主题词是经过规范化的词，在确定主题词时，要对论文进行主题分析，依照标引和组配规则转换成主题词表中的规范词语。（参见《汉语主题词表》和《世界汉语主题词表》）。 5、论文格式的论文正文： （1）引言：引言又称前言、序言和导言，用在论文的开头。引言一般要概括地写出作者意图，说明选题的目的和意义, 并指出论文写作的范围。引言要短小精悍、紧扣主题。 〈2）论文正文：正文是论文的主体，正文应包括论点、论据、论证过程和结论。主体部分包括以下内容： 提出问题-论点； 分析问题-论据和论证； 解决问题-论证方法与步骤； 结论。 6、论文格式的参考文献 一篇论文的参考文献是将论文在研究和写作中可参考或引证的主要文献资料，列于论文的末尾。参考文献应另起一页，标注方式按《GB7714-87文后参考文献著录规则》进行。 中文：标题--作者--出版物信息（版地、版者、版期） 英文：作者--标题--出版物信息 所列参考文献的要求是： （1）所列参考文献应是正式出版物，以便读者考证。 （2）所列举的参考文献要标明序号、著作或文章的标题、作者、出版物信息。
直接套用 著录格式： [序号] 编著者 书名[M] 版本，出版地：出版者，出版年 页码例：[1] 刘谋佶, 吕志咏, 丘成昊 边条翼与旋涡分离流[M] 北京: 北京航空学院出版社, 24~[2] Isidori A Nonlinear control systems[M] 2nd, New York: Springer Press, 32~注：初版书不标注版本，页码是可选项。（2） 期刊著录格式： [序号] 作者 题目[J] 刊名，年，卷（期）：页码例：[1] 傅惠民 二项分布参数整体推断方法[J] 航空学报，2000，21（2）： 155~[2] Moustafa G H Interaction of axisymmetric supersonic twin jets[J] AIAA J, 1995, 33(5): 871~注：外文期刊的刊名可用简称；请注意标注文章的年、卷、期、页，不要遗漏。（3） 学位论文 著录格式：[序号] 作者 题目[D] 地点：单位，年例：[1] 朱刚 新型流体有限元法及叶轮机械正反混合问题[D] 北京：清华大学，[2] Sun M A study of helicopter rotor aerodynamics in ground effect[D] Princeton: Princeton Univ, （4）论文集，会议录 著录格式：[序号] 主要责任者 题名[C] 出版地：出版者，出版年例：[1] 辛希孟 信息技术与信息服务国际研讨会论文集：A集[C] 北京：中国社会科学出版社，[2] 北京空气动力研究所 第九届高超声速气动力会议论文集[C] 北京：北京空气动力研究所，（5）论文集中析出的文献 著录格式：[序号] 作者 题目[A] 见：主编 论文集名[C] 论文集名 出版地：出版者，出版年：页码例：[1] 陈永康，李素循，李玉林 高超声速流绕双椭球的实验研究[A] 见：北京空气动力研究所编 第九届高超声速气动力会议论文集[C] 北京：北京空气动力研究所，1997：9～[2] Peng J, Luo X Z, Jin C J The study about the dynamics of the approach glide-down path control of the carrier aircraft[A] In: GONG Yao-nan Proceedings of the Second Asian-Pacific Conference on Aerospace Technology and Science[C] Beijing: Chinese Society of Aeronautics and Astronautics, 1997: 236～注：会议文集的出版者可能不是正式的出版社；出版地指出版者所在地，不一定是会议地点。（6）科技报告 著录格式：[序号] 作者 题名[R] 报告题名及编号，出版地：出版者，出版年例：[1] 孔祥福 FD-09风洞带地面板条件下的流场校测报告[R] 北京空气动力研究所技术报告 BG7-270，北京：北京空气动力研究所，[2] Carl E J Analysis of fatigue, fatigue-crack propagation and fracture data[R] NASA CR-132332, 注：对于NASA报告，AIAA Paper等航空航天领域知名报告，出版地和出版者可以省略。（7）国际、国家标准，行业规范 著录格式：[序号] 标准编号，标准名称[S] 出版地：出版者，出版年例：[1] MIL-E-5007 D, 航空涡轮喷气和涡轮风扇发动机通用规范[S] 美国空军，[2] GB 7713-87, 科学技术报告、学位论文和学术论文的编写格式[S]注：对于国标GB等，出版地、出版者和出版年可省略。（8）专利 著录格式：[序号] 设计人 专利题名[P] 专利国别：专利号，公告日例：[1] 黎志华，黎志军 反馈声抵消器[P] 中国专利：ZL85100748，1986－09－
朋友你好，根据我多年从事文字工作的经验，我认为：如果投稿更有针对性，命中率会更高一些。这就关系到，你是哪里的？干什么的？写的稿件是什么体裁？什么内容？如果说投稿的话，最好投当地的报刊、网络或者是你从事的职业报刊发表，要投哪个媒体首先要研究哪个媒体，看它需要什么内容、什么体裁、什么格式的稿件，“对症下药”，这样会更轻松一些、方便一些，命中率会更高一些。如果你能够告诉我你的具体情况（干什么工作，哪里的，写的小说的大致内容等），我可以给你一些建议。我1993年开始在部队时开始发表各类文章，包括：报告文学、新闻、诗歌、散文、小说、评论等体裁的，到目前，先后在《人民日报》《法制日报》《农民日报》《中国文化报》《法制文萃》《半月谈》《解放军报》《中国国防报》《中国绿色时报》《中国日报》《中国教育报》《人民公安报》《中国交通报》《中国安全生产报》《中国转业军官》《中国人事》《道路交通管理》等报刊发表的大约5000篇左右吧，有40多篇获奖。另外：投稿时，第一要有信心，第二要投对报刊媒体，这两点非常重要。祝你成功！
《中国民用航空学报》属于一号多刊，不是正规期刊。但是有的院校检查不严，可以用。有的学校查的严格，就用不了啊。空中交通管理审稿是一个月。