李学杰万荣胜黄向青陈太浩
（广州海洋地质调查局 广州 510760）
第一作者简介：李学杰，男，1964年生，博士，教授级高工，主要从事海洋地质与第四纪地质研究工作，E-mai1：Xuejie1i＠yeah.net。
摘要 近岸极浅水区的水深复杂多变，因受调查条件所限，往往成为测量的盲区。本文利用Landsat ETM多光谱影像，结合实测数据，探讨其对近岸水深反演的可行性。结果表明，对北部湾海域，ETM2是适合水深反演的最佳波段，水深Z与光谱辐射值L2及深水辐射值LS2之间存在以下关系：Z=-17.19 1n（L2-LS2）+56.40。从反演效果来看，总体上能反映水深的变化趋势，同时一些小的地形变化，如潮沟等也能得到较好的体现。因此对近岸较清澈的海域，遥感影像的水深反演可以在一定程度上弥补实测的不足。
关键词 遥感 测深 北部湾
1 前言
水深是海洋环境的重要参数，长期以来，采用船载回声测深法进行测量，取得很好的效果。但对极浅水环境（如0～3m）的水深测量却是传统调查方法的盲区，因为调查船和测量人员均难以到达。而极浅水环境是真正的海陆共同作用的区域，环境脆弱且复杂多变。准确测量水深，不仅对了解地形地貌，而且对海岸保护和建设均十分重要。
可见光对水体有一定的穿透性，因此有可能利用其这种特性进行水深反演，尤其对海岸带，有可能形成对传统调查方法的补充。早在20世纪60年代末，美国密歇根环境研究所的一个小组就开始从事遥感测深的研究，利用MSS，TM和航片等多光谱数据和一些同步测量的海况数据进行测深模型研究，提出了基于地面反射的遥感测深理论，并长期对此方向进行探索［1～3］。此外，其他国家也进行类似的研究和应用［4，5］。近年来，相关的研究还在不断增加，应用领域也在拓宽［6，7］。
国内对于水深的定量研究起步较晚，20世纪90年代初，我国开始利用遥感进行水深方面的研究［8，9］。张鹰等［10］利用遥感研究近岸水深及潮滩的冲淤变化。地矿部航空物探遥感中心从1994年开始，利用南沙群岛海域18个景区26个时相的TM数据，开展遥感水深调查和制图，并取得良好的效果［11］。庞蕾等［12］介绍了水深遥感的不同方法。近年来国家海洋局对近岸海洋遥感做了大量工作，也包括对水深遥感的探讨［13，14］。
2 遥感测深的理论模型
利用星载多光谱数据进行浅海水深测量，其物理基础是可见光各个波段对于水体均具有一定的穿透力，如Landsat ETM1波段对水体的穿透深度最大，在清洁水的情形下可以穿透30m以上［12］；ETM2波段可达10～15m；但ETM3的水体穿透力则相对较弱。对于各种类型的水体，可见光的水体衰减系数最小值都出现在蓝绿波段之间，表明 Landsat ETM1，2，3波段是通常的遥感测深的最佳波段［4］。
遥感的水深模型有多种，主要包括解析法和统计法。前者是利用传感器所接收的辐射亮度建立其于底质反射率及水深的解析表达式从而计算水深，但许多参数计算，依赖于对大气影响的准确校正。后者是利用实测点回归得到辐射亮度与水深之间的关系，进而推求未知水深点的水深，根据利用波段数的不同，统计法又分为单波段法、双波段比值法、线性多波段法等三种。
基于海底反射的模型为
南海地质研究.2007
式中：Li是传感器接收到的第i波段的辐射值；Lsi是深水区辐射值，它反映的是水面反射、水体散射及大气散射等的总和，而不包含底质反射；Ci是与太阳辐射度、大气和水面透过率及水面折射有关的参数；Rbi是底反射率；ki是水体的衰减系数；f是水体路径长度（通常取2）；Z为水深。
将（1）式取自然对数得
Z=1n（CiRbi）/fki-1n（Li-Lsi）/fki （2）
假设底质反射率Rbi是常数，大气和海况是均一的，即衰减系数ki是常数，并设
Ii=1n（Li-Lsi），其中Li-Lsi代表海底发射值。
a=-1/fki
b=1n（CiRbi）/fki
那么：
Z=aIi+b （3）
其中系数a、b可用线性回归方法求得，这就是单波段线性回归模型。
将单波段和双波段模型推广到多波段，则有
Z=A0+A1I1+A2I2+⋯⋯+AnIn （4）
用多元回归方法求其系数，形成多波段模型。
从传感器第i波段的影像中，实测出一组Z～Xi的值，利用最小二乘法可以计算出a，b的值或A0，A1，A2⋯⋯An值，然后利用以上公式推算出其他未知水深点的水深值。
3 实验区的选择与采样分析
作为水深反演的海域，要求水体清澈，悬浮物、叶绿素及各种溶解有机质少，透光性好。广州海洋地质调查局于近年对大亚湾、大鹏湾、珠江口及北部湾等近岸海域进行调查，其中北部湾海域水体最清澈，较适合做水深反演，因此选择该区进行实验研究。
对北部湾钦州湾海域的调查是2006年进行的，由于沿测线方向测深点很密，测线之间的距离相对大得多，因此将测深点的数据抽稀（图1），以便对比。所测量的最小水深为3.3m。所采用的影像是2000年11月16日的Landsat7 ETM影像，尽管影像时间和实测数据采集时间有一定的差异，但考虑在这短期内水深的变化总体不大，可以适用。
图1 北部湾钦州湾海域位置及实测水深图
Fig.1 Location and bathymetry in Qingzhou Bay，Daya Bay
将水深测量数据的坐标与遥感影像坐标统一到UTM WGS84 坐标系，按一定的网格，选有实测水深的点在遥感影像上进行采样，读出对应该点影像的1，2，3波段的DN值，共采样88个点位。同时选择确定该影像最深水区1，2，3波段的DN值为73，44和31，分别代表这3个波段的深水辐射值，两者的差值（Li-Lsi）代表海底发射值。
图2 Landsat ETM1，2，3波段海底辐射值（Li-LSi）与水深关系
Fig.2 Relationship between the seabed reflectance of Landsat ETM band 1，2，3（Li-LSi）and water depth
从各波段底质发射值与实测水深关系（图2）来看，尽管波段1的水体穿透性最好，但本区与水深关系并不密切，这可能是由于该波段（蓝光）在本区受到的干扰较多之故，而波段2与水深关系最密切，因此采用该波段数据进行拟合（图2B）。
拟合结果：Z=-17.191n（L2-LS2）+56.40
其中Z为水深，L2是ETM2波段的DN值，LS2是ETM2深水辐射值，本文取44。
对该回归方程进行显著性检验，计算的剩余平方和Q=800.7，回归平方和U=1796.9，数据个数n=88，采用F检验：
南海地质研究.2007
在a=0.01，自由度为（1，86）条件下，查得其临界值F0.01（1，86）=6.94。F＞F0.01，表明在置信水平为99%，水深Z与Landsat ETM2的DN值L2与该波段的深水DN值LS2之差的对数，即1n（L2-LS2）之间是显著相关的，两者之间拟合的方程是有效的。
4 水深反演结果
利用上述实测数据与遥感影像的拟合结果，对遥感数据进行水深反演。从结果来看，总体效果不错，从近岸向外，水深呈增加趋势，尤其是较浅水海域，其效果更好，而且与实测结果基本吻合（图3）。
图3 钦州湾反演水深与实测等深线的对比
Fig.3 Comparison Water depth calculated from image to real one in Qingzhou area
图4 北部湾（大区）反演水深与实测水深的对比
Fig.4 Comparison water depth calculated from image to real one in Beibu Bay area
图5 北部湾西部反演水深
Fig.5 Water depth calculated from image in Western Beibu Bay
把反演的海域扩大，可以看出，其总体变化趋势依然与实际水深变化较吻合（图4）。北海东部银滩、北海港北部的浅水区域均得到体现，同样东南部海域因靠近其南部的涠州岛而水深变浅，也得到反映。
西部海域的反演效果似乎更好，永实岛南北的水深差异明显，南部水深明显大于北部（岛链内侧）应与实际吻合（图5）。永实岛之间的槽沟以及小岛周围的浅水区等均得到较好的体现，进一步表明该水深反演方法可能解决一些问题，成为实测的补充。
5 结论与讨论
通过结合实测数据，对北部湾的钦州湾海域进行遥感影像的水深反演表明，Landsat ETM2波段数据较适合于水深反演，且拟合的方程Z=-17.191n（L2-LS2）+56.40反演效果较好，与实际水深总体有较好的一致性，因此该方法可以在一定程度作为实测方法的补充。
同时也应该看到该方法的局限性，首先水体所含物质（包括悬浮物、叶绿素及溶解有色有机质等）对遥感辐射值有很大的影响，因此该方法只适合于清澈的水体。其次不同的底质，其反射率可能不同，对反演效果也将产生一定的影响。
此外，大气条件的空间差异以及影像几何校正的精度等均可能影响采样值，并因此影响拟合方程的效果。而且本文拟合方程时，缺乏小于3m的实测数据值，对拟合结果也产生一定影响。
参考文献
［1］Lyzenga D R.Passive remote sensing techniques for mapping water depth and bottom features.Applied Optics，1978，17（3）：379～383
［2］Lyzenga D R.Remote sensing of bottom reflectance and Water attenuation Parameters in shallow Water using aircraft and Landsat data.Int.J.Remote Sensing，1981，2（1）：71～82
［3］Nordman M E.Water Depth Extraction form Landsat-5 Imagery.Proc.Of 23th Int.Sym.On Remote Sensing of Env，1990
［4］Willian D P.Bathymetry mapping With Passive multispectral imagery.Applied Optics，1989，28（8）：1569～1578
［5］BierWirth P N，Lee T J，Burne R V.Shallow sea-floor reflectance and Water depth derived by UnmiXing multispectral imagery［J］.PE＆RS，1993，59（3）：331～338
［6］Isoun E，Fletcher C，Frazer N，Gradie J.Multi-spectral mapping of reef bathymetry and coral cover；Kailua Bay，Hawaii.Coral Reefs，2003，22：68～82
［7］Karpouzli E，Malthus T，Place C，et al.Underwater light characterization for correction of remotely sensed images.Int.J.Remote Sensing，2003，24（13）：2683～2702
［8］李铁芳等.卫星海洋遥感信息提取和应用.北京：海洋出版社，1990
［9］李铁芳等.浅海水下地形地貌遥感信息提取与应用.环境遥感，1991，6（1）：32～36
［10］张鹰，丁贤荣，王文.水深遥感与潮滩地形冲淤变化分析.港口工程，1998，（2）：26～30
［11］邸凯昌，丁谦，曹文玉.南沙群岛海域浅海水深提取及影像海图制作技术.国土资源遥感，1999，（3）：59～64
［12］庞蕾，聂志峰.星载多光谱浅海水深测量方法.山东理工大学学报（自然科学版），2003，17（6）：59～61
［13］恽才兴（主编）.海岸带及近海卫星遥感综合应用技术.北京：海洋出版社，2005
［14］杨晓梅，周成虎，杜云艳等.海岸带遥感综合技术与实力研究.北京：海洋出版社，2005
［15］韩震，恽才兴.伶仃洋大铲湾潮滩冲淤遥感反演研究，海洋学报，2003，（5）：58～64
Bathymetry in CoaStal area by LandSat ETM：Method and its Application in Beibu Bay
Li XuejieWan RongshengHuang XiangqingChen Taihao
（Guangzhou Marine Geological Survey，Guangzhou，510760）
Abstract：It is quite difficult to measure the water depth in very shallow area which is much variable，due to difficult arriving for survey ship.It is try to calculate Water depth by multi-spectral Landsat ETM image，combining the real measured data，in the Paper.The result suggested that the band 2 of Landsat ETM is better for calculating in the Beibu Bay and fitted formula is Z=-17.191n（L2-LS2）+56.40，Where Z is Water depth，L2 and LS2 is reflectance and deep Water reflectance of band 2 of band respectively.The calculating result can better fit for the real data，not only for the basic trend，but even for the tidal channel.Therefore it can be concluded that the method of calculating Water depth by remote sensing is suitable for costal clear Water area and be complementarity for real measure.
Key Words：Remote sensing Bathymetry Beibu Bay
山东理工大学学报(社会科学版)和山东理工大学学报(自然科学版)什么核心期刊都不是，就是最最普通的省级期刊。
影像科学与光化学 [1674-0475] 本刊收录在： 中国科学引文数据库(CSCD)来源期刊(2009-2010)提示: CSCD核心库(C)本刊收录在： 中国科学引文数据库(CSCD)来源期刊库(2013-2014)提示: CSCD核心库(C)本刊收录在： 中国科学引文数据库(CSCD)来源期刊核心库(2011-2012)本刊收录在： 中国科技期刊引证报告(2010年版)提示: 《引证报告》2010年版影响因子:0.413本刊收录在： 中国科技期刊引证报告(2011年版)提示: 《引证报告》2011年版影响因子:0.368本刊收录在： 中国科技期刊引证报告(2012年版)本刊收录在： 中国科技期刊引证报告(2013年版)提示: 《引证报告》2013年版影响因子:0.475本刊收录在： 中文核心期刊要目总览(2008年版)提示: 排序：化学、晶体学 - 第24位本刊收录在： 中文核心期刊要目总览(2011年版)提示: 排序：化学，晶体学类 - 第21位主题分类：O6,O7：化学，晶体学: O6,O7：化学，晶体学该刊物不属于EI检索，谢谢。
应该是1000-3231 感光科学与光化学，是EI的。EI《Ei Compendex》收录的中国期刊ISSN 期刊刊名0567-7718 Acta Mechanica Sinica/Lixue Xuebao1006-7191 Acta Metallurgica Sinica (English Letters)0253-4827 Applied Mathematics and Mechanics (English Edition)1004-5341 China Welding (English Edition)1004-9541 Chinese Journal of Chemical Engineering1022-4653 Chinese Journal of Electronics1000-9345 Chinese Journal of Mechanical Engineering (English Edition)1671-7694 Chinese Optics Letters1006-6748 High Technology Letters1004-0579 Journal of Beijing Institute of Technology (English Edition)1005-9784 Journal of Central South University of Technology (English Edition)1553-9105 Journal of Computational Information Systems1000-1484 Journal of Dong Hua University (English Edition)1005-9113 Journal of Harbin Institute of Technology (New Series)1001-6058 Journal of Hydrodynamics1548-7741 Journal of Information and Computational Science1005-0302 Journal of Materials Science and Technology1002-0721 Journal of Rare Earths1003-7985 Journal of Southeast University (English Edition)1004-4132 Journal of Systems Engineering and Electronics1003-2169 Journal of Thermal Science1005-8850 Journal of University of Science and Technology Beijing: MineralMetallurgy Materials (Eng Ed)1009-3095 Journal of Zhejiang University: Science1001-0521 Rare Metals1006-9291 Science in China, Series B: Chemistry1003-6326 Transactions of Nonferrous Metals Society of China (English Edition)1006-4982 Transactions of Tianjin University1007-0214 Tsinghua Science and Technology0253-4177 半导体学报1001-1455 爆炸与冲击1001-5965 北京航空航天大学学报1001-053X 北京科技大学学报1001-0645 北京理工大学学报1000-1522 北京林业大学学报1007-5321 北京邮电大学学报1005-0299 材料科学与工艺1009-6264 材料热处理学报1005-3093 材料研究学报1007-7294 船舶力学1000-8608 大连理工大学学报1000-2383 地球科学 — 中国地质大学学报1005-0388 电波科学学报1000-1026 电力系统自动化0372-2112 电子学报1005-3026 东北大学学报（自然科学版）1001-0505 东南大学学报 （自然科学版）1000-3851 复合材料学报1003-9015 高校化学工程学报1000-5773 高压物理学报1000-4750 工程力学1001-9731 功能材料1000-3819 固体电子学研究与进展1006-2793 固体火箭技术1005-0086 光电子.激光1004-924X 光学精密工程0253-2239 光学学报0454-5648 硅酸盐学报1006-7043 哈尔滨工程大学学报0367-6234 哈尔滨工业大学学报0253-360X 焊接学报1005-5053 航空材料学报1000-8055 航空动力学报1000-6893 航空学报0258-0926 核动力工程1001-9014 红外与毫米波学报1000-2472 湖南大学学报 （自然科学版）1000-565X 华南理工大学学报（自然科学版）1671-4512 华中科技大学学报（自然科学版）0438-1157 化工学报0577-6686 机械工程学报1671-5497 吉林大学学报（工学版）1006-5911 计算机集成制造系统0254-4164 计算机学报1007-4708 计算力学学报1001-246X 计算物理0258-1825 空气动力学学报1000-8152 控制理论与应用1001-0920 控制与决策1004-0595 摩擦学学报1000-0925 内燃机工程1000-0909 内燃机学报1005-2615 南京航空航天大学学报1005-9830 南京理工大学学报 （自然科学版）1001-4322 强激光与粒子束1000-0054 清华大学学报 （自然科学版）1006-8740 燃烧科学与技术1000-9825 软件学报0254-0150 润滑与密封1006-2467 上海交通大学学报1671-2021 沈阳建筑大学学报（自然科学版）0371-0025 声学学报1001-8719 石油学报:石油加工1009-3087 四川大学学报(工程科学版)0254-0096 太阳能学报0493-2137 天津大学学报0253-374X 同济大学学报 （自然科学版）1001-4055 推进技术1000-324X 无机材料学报1001-2400 西安电子科技大学学报0253-987X 西安交通大学学报1000-2758 西北工业大学学报1007-8827 新型 炭材料1000-6915 岩石力学与工程学报1000-7598 岩土力学1005-0930 应用基础与工程科学学报1000-6931 原子能科学技术1008-973X 浙江大学学报 （工学版）0253-9748 真空科学与技术学报1004-4523 振动工程学报1000-3835 振动与冲击0258-8013 中国电机工程学报1022-0666 中国航空太空学会 学刊0258-7025 中国激光1000-1964 中国矿业大学学报1001-4632 中国铁道科学0254-508X 中国造纸0529-6579 中山大学学报（自然科学版）0254-4156 自动化学报=========================================EI非核心版（《EiPageOne》）收录的中国期刊ISSN 期刊刊名1000-9116 Acta Seismologica Sinica English Edition1561-8625 Asian Journal of Control1005-9040 Chemical Research in Chinese Universities0890-5487 China Ocean Engineering1001-6279 International Journal of Sediment Research0254-9409 Journal of Computational Mathematics1007-6417 Journal of Shanghai University0257-9731 Journal of the Chinese Society of Mechanical Engineers, Transactionsof the Chinese Institute of Engineers1000-2413 Journal of Wuhan University of Technology -Materials Science Edition1001-8042 Nuclear Science and Techniques/Hewuli1002-0071 Progress in Natural Science1006-9283 Science in China (Series A: Mathematics)1007-1202 Wuhan University Journal of Natural Sciences1001-5868 半导体光电1000-1506 北方交通大学学报0479-8023 北京大学学报 (自然科学版)1007-2640 北京化工大学学报 （自然科学版）1000-1093 兵工学报1001-4381 材料工程1007-4112 长安大学学报1004-499X 弹道学报1000-6753 电工技术学报1001-0548 电子科技大学学报1005-9490 电子器件1009-5896 电子与信息学报1000-6761 动力工程1001-3784 粉末冶金技术1000-3436 辐射研究与辐射工艺学报1000-3231 感光科学与光化学0449-749X 钢铁0251-0790 高等学校化学学报1003-6520 高电压技术1000-7555 高分子材料科学与工程1002-0470 高技术通讯0254-3052 高能物理与核物理1001-1609 高压电器0253-231X 工程热物理学报1007-4252 功能材料与器件学报1003-501X 光电工程1000-0593 光谱学与光谱分析1002-1582 光学技术1004-4213 光子学报1001-2486 国防科技大学学报1009-606X 过程工程学报1006-9941 含能材料1002-0837 航天医学与医学工程1002-1396 合成树脂及塑料0258-0934 核电子学与探测技术0253-9950 核化学与放射化学0253-3219 核技术0254-6086 核聚变与等离子体物理1672-9102 湖南科技大学学报1006-5431 华北工学院学报1006-3080 华东理工大学学报 （自然科学版）1000-3932 化工自动化及仪表1001-7631 化学反应工程与工艺1005-9954 化学工程0250-3301 环境科学1002-0446 机器人1001-9669 机械强度1000-1158 计量学报1003-9775 计算机辅助设计与图形学学报1000-3428 计算机工程1000-1239 计算机研究与发展1007-9629 建筑材料学报1000-6869 建筑结构学报1006-852X 金刚石与磨料磨具工程0254-6051 金属热处理0412-1961 金属学报1001-5493 离子交换与吸附1008-0562 辽宁工程技术大学学报 （自然科学版）0253-2417 林产化学与工业1007-3124 流体力学实验与测量0253-9993 煤炭学报1003-6059 模式识别与人工智能1001-1935 耐火材料1000-1972 南京邮电学院学报1002-6819 农业工程学报1000-1298 农业机械学报0253-2409 燃料化学学报1001-2060 热能动力工程1000-985X 人工晶体学报1007-6735 上海理工大学学报1000-5870 石油大学学报 （自然科学版）1000-7210 石油地球物理勘探1000-8144 石油化工1006-396X 石油化工高等学校学报1000-0747 石油勘探与开发0253-2697 石油学报1004-9037 数据采集与处理1001-6791 水科学进展1003-1243 水力发电学报0559-9350 水利学报1007-2012 塑性工程学报1001-2249 特种铸造及有色合金1001-8360 铁道学报1003-8213 微细加工技术1000-050X 武汉大学学报 (信息科学版)1671-4431 武汉理工大学学报1006-2823 武汉理工大学学报 (交通科学与工程版)1000-3290 物理学报1006-7930 西安建筑科技大学学报1001-5361 西安石油 大学学报（自然科学版）0258-2724 西南交通大学学报1000-2634 西南石油学院学报 （自然科学版）1002-185X 稀有金属材料与工程1004-731X 系统仿真学报1000-6788 系统工程理论与实践1001-506X 系统工程与电子技术0253-4320 现代化工1004-2474 压电与声光1000-4548 岩土工程学报1000-7571 冶金分析0254-3087 仪器仪表学报1000-4939 应用力学学报1000-1328 宇航学报1008-9209 浙江大学学报 （农业与生命科学版）1004-6801 振动测试与诊断1001-7372 中国公路学报1000-6923 中国环境科学1004-132X 中国机械工程1004-0609 中国有色金属学报1000-6842 中国造纸学报1005-9792 中南工业大学学报 （自然科学版）1006-7329 重庆建筑大学学报1001-4977 铸造1000-8365 铸造技术Ei Compendex Web（EI网络版）是《Ei Compendex》和《Ei PageOne》合并而成的Internet版本。Ei来源期刊分三个档次 ：（1）核心期刊， Compendex 数据库收录重点是下列主要工程学科：化学工程，土木工程；电子/电气工程，机械工程，冶金、矿业、石油工程，计算机工程和软件。美国工程信息公司副总裁Katz认为：“这些是Compendex数据库的‘核心’领域。”目前，核心期刊约有1000种；每期所有论文均被录入Compendex。国内《金属学报》、《清华大学学报》等为核心期刊。（2）选择期刊，一些学科领域的期刊是有选择地收录，包括：农业工程、工业工程、纺织工程、应用化学、应用数学、应用力学、大气科学、造纸化学和技术，高等学校工程类学报等。Compendex 只选择与其主题范围有关的文章，并不是所以文章均被收录。目前，选择期刊约1600种，国内以上学科的期刊大多数为选择期刊。（3）扩充期刊，它只收录题录（Ei Page One）。在Ei的扩充版中，约2800种期刊。国内期刊比较多。
光化学的定义有不同的表述。C. H. Wells认为，光化学研究的是“吸收了紫外光或可见光的分子所经历的化学行为和物理过程”。N. J. Turro则认为“光化学研究的是电子激发态分子的化学行为和物理过程”。由于电子激发态通常由分子吸收紫外光或可见光形成，所以上述两种定义的实质是一样的。光化学是研究光与物质相互作用所引起的永久性化学效应的化学分支学科。由于历史的和实验技术方面的原因，光化学所涉及的光的波长范围为 100～1000纳米，即由紫外至近红外波段。比紫外波长更短的电磁辐射，如X或γ射线所引起的光电离和有关化学属于辐射化学的范畴。至于远红外或波长更长的电磁波，一般认为其光子能量不足以引起光化学过程，因此不属于光化学的研究范畴。观察到有些化学反应可以由高功率的红外激光所引发，但将其归属于红外激光化学的范畴。影像是人对视觉感知的物质再现。影像可以由光学设备获取，如照相机、镜子、望远镜及显微镜等；也可以人为创作，如手工绘画图像等。影像也是一种视觉符号。通过专业设计的影像，可以发展成人与人沟通的视觉语言，也可以了解世界美术中大量的平面绘画、立体雕塑与建筑。
学报自dyuoiu
刘林冯伟宋宪生宋哲李卫陈擎
（核工业二〇三研究所，陕西咸阳712000）
[摘要]介绍了查查香卡铀矿床勘查过程和矿区地质特征，系统总结了该矿床铀资源调查评价中取得的主要成果，提出了区域性北西西向断裂构造破碎带＋奥陶－志留系滩涧山群变质火山岩＋华力西－印支期中酸性岩浆－热液活动的“三位一体”成矿机理，认为查查香卡铀矿床具有形成铀多金属矿的良好前景。
[关键词]查查香卡铀矿床；成矿机理；铀多金属；潜力
查查香卡铀矿床属青海省乌兰县赛什克乡南柯柯村管辖，位于乌兰县城南的南坡，海拔4000～4200m。矿床南部有109国道，由398km里程碑向西北方向进入县级公路到达查查香卡农场，行程24km。由查查香卡农场北东有汽车便道行驶20km可直达矿床（图1）。矿床的西部及西北部从乌兰县或德令哈横穿牦牛山，到达查查香卡农场。
图1 查查香卡铀矿床交通位置图
1—城镇；2—铁路；3—国道；4—县级公路；5—河流；6—湖泊；7—查查看卡铀矿床
1发现和勘查过程
（1）第一阶段：矿床发现与初期勘查（1969～1971）
1969～1971年二机部西北一八二队对查查香卡地区进行了初步勘查。
1969年由二机部西北一八二队第三分队在1∶2.5万伽马普查中发现了17、18两条异常带，从而为查查香卡铀矿床（原324小型矿床）的发现提供了重要线索[1]。
1970年通过1∶10000伽马详查、1∶2000伽马详测、槽探系统控制长2000m，手掘硐探4个，确定324地区中段为主要矿化段[1]。
1971年机掘硐探1个，手掘2个，钻孔3个，控制深度200m，控制长度500m，了解了324地区中段深部矿化，发现1个工业矿孔。
1969～1971年完成的工作量见表1[2]。发现工业矿体13个，单个长35～70m，厚0.75～3.5m，品位0.054%～0.073%，平均长49.7m，平均品位0.062%，埋深0～50m；后备矿体26个，单个长35～210m，厚0.7～5.0m，品位0.035%～0.046%，平均长70.5m，平均厚1.83m，平均品位0.038%，埋深0～100m。
表1 1969～1971年完成的工作量统计
（2）第二阶段：停滞期（1972～2008）
根据核工业部三局南昌会议精神，遵照“富、近、浅、易”原则，结合324小型矿床的具体情况，经青海652大队研究，核工业部三局批准，1972年撤出324地区。
（3）第三阶段：评价期（2009～2012）
2009～2012年中国核工业地质局核工业二〇三研究所在本区开展了铀矿调查评价工作。
2009年在查查香卡324地区中段施工探槽10个，长度387.1m，发现4个地表工业矿体。
2010年在324地区中段施工了4条勘探线，5个钻孔，完成工作量1853.21m。发现2个工业铀矿孔，2个铀矿化孔。
2011年在查查香卡324地区完成实测地质剖面12条，约8.94km，发现4～5条铀异常带；完成1∶2000地面伽马能谱测量11条、面积12km2；施工钻孔10个，完成工作量3179.76m，发现1个工业铀矿孔，1个异常孔。
2012年在查查香卡地区完成1∶2000地质剖面4条，长度7.01km，发现Cu、Au等金属矿化；1∶1万激电测量60个点3条激电测深剖面，显示可能为两条金属硫化物异常带。表层矿化较弱，向深部500 m左右可能有较好多金属矿化；施工5个钻孔，完成1883.44m，发现1个工业铀矿孔、1个异常孔，在C76钻孔中见到两层黄铁矿化层（16.5～4.3m）。初步确定查查香卡矿化带为铀、多金属矿化带。
根据地表槽探圈定了4个工业铀矿体，钻孔在矿化带中部查证发现了4个工业铀矿孔、2个铀矿化孔，铀矿体5～7层，累计厚度10.82～8.12m，品位0.05%～0.25%。同时发现了该区铌、稀土等多金属矿化，且与铀呈正相关关系。在铀矿体周围发现铌矿体3个、矿化体11个，控制矿体长200～357m，厚2m左右，品位0.0883%～0.2498%。还发现有稀土矿化，圈定矿化体5个，主要为铈、镧、钇等矿化。
2矿床基本特征
矿床位于丁字口－埃姆尼克山－牦牛山新元古代－早古生代岩浆弧带与赛什腾－锡铁山－哇洪山新元古代－早古生代缝合带的结合部，夹持于泽日肯和陶力古隆起之间托莫尔日特蛇绿混杂岩带中。
2.1矿区地质
查查香卡铀矿床属赛什腾山－阿尔茨托山铀成矿带乌兰南山铀成矿亚带。位于柴达木北缘残山断褶带东段，华力西期下义山花岗闪长岩体外接触带中（图2），矿化赋存在奥陶－志留系滩涧山群（OST）变质火山岩中，受北西向挤压破碎带控制，长约9km，宽约2km。
矿床内主要矿化围岩为奥陶－志留系滩涧山群变质火山岩，主要以各类绿片岩为主，夹云母石英片岩及大理岩。地层产状北西西向。岩浆岩主要有华力西期中粗粒花岗闪长岩、条纹状混杂岩、中细粒闪长岩、条带状闪长岩等；脉岩有中细粒石英闪长岩脉、闪长斑岩脉、辉长岩脉、斑状角闪正长岩脉、石英脉等。
矿化赋存于北西向斜长角闪片岩挤压构造破碎带中，其中富矿地段复合于斜长角闪片岩之早期北西向构造带上，长2km，宽数米至百余米，倾向北东，倾角65°；矿后局部有北东向构造叠加，石英脉沿晚期北东向构造充填。该区伽马异常带断续长9km，呈北西向展布，在其中2km长的范围内，铀、钍矿化较好，矿化断续长几十米至几百米，宽几十厘米至十几米；该矿床同时为多种稀土元素的成矿区，主要有铌、铈、镧、铱等矿化[2,3]。
图2 查查香卡铀矿床矿体平面展布图
2.2含矿岩石岩性特征
查查香卡矿床围岩为碎裂斜长角闪片岩等（图3），近矿围岩蚀变强烈，主要有钠长石化、赤铁矿化、绿泥石化、碳酸盐化。铀含量较高，坑道口碎石一般为（200～300）×10-6，最高达878×10-。6 矿化岩石为灰黑色碎裂闪长岩，灰黑色，块状，主要成分为角闪石，镜下鉴定为斜长角闪片岩，岩石呈纤维状变晶结构、片状构造，说明矿石经过了一定的动力构造变质作用，片理化发育，岩石较破碎；矿石主要成分为角闪石（70%）、斜长石（20%），少量黑云母（3%）、石英（3%）、方解石（4%），此外还有微量的绢云母、阳起石和绿帘石等；岩石中大量角闪石呈不规则纤柱状晶体，彼此嵌连而平行排列，大小为（0.12～0.15）mm×0.11mm，组成较为密集的条带，角闪石不同程度遭受阳起石化蚀变；斜长石双晶不明显，局部有挤压破碎，具绿帘石化、绢云母化现象。
2.3铀矿化带特征
查查香卡铀矿床可以划分为西、中、东3段[2]。西段有3条异常带，长500～900m，宽2～10m，一般为（100～200）×0.1258nC/kg.h，最高大于1000×0.1258nC/kg.h；中段异常断续长200m，宽几十米至百余米（由多条带组成），一般（100～300）×0.1258nC/kg.h，最高大于1000×0.1258nC/kg.h（图4）；东段异常断续长3000m，一般（50～150）×0.1258nC/kg.h，最高大于1000×0.1258nC/kg.h。
铀矿化主要分布于含矿构造带闪长碎裂岩（或糜棱岩）及碎裂闪长岩中，富矿地段为碎裂岩或糜棱岩，岩石破碎强烈、弱硅化、片理发育；碳酸盐呈细脉或星点状分布多的地段矿化比较好，其他多为贫矿化地段。
镜下见沥青铀矿充填交代黄铁矿现象，铀矿化以浸染分散状或吸附形式存在于含矿岩石中。矿化基本有3种表现形式：碳酸盐成细脉状沿岩石片理充填成条状构造或呈星点状分布于含矿岩石中，一般矿化较富；白色长石－石英－碳酸盐脉成细脉状沿片理分布于含矿岩石中；岩石片理化发育，为黑绿色或紫色，肉眼看不到脉体，但矿化较好。对所采集的查查香卡矿床矿石、含矿岩石和围岩进行的铀钍、U4＋和U6＋分析表明，查查香卡矿床铀矿化蚀变碎裂岩的铀含量为0.1022%～0.11011%，而围岩斜长角闪片岩的铀含量较低，为（15～18）×10-6；微量钍分析表明，矿化岩石中钍含量比较高，为（2916～659）×10-6，而围岩的钍含量较低，只有116×10-6，矿床中铀以原始富集为主。另外通过矿化岩石U4＋和U6＋含量的对比也可以看出铀的富集特征，矿化岩石 U4＋的含量分别为0.10674%和124×10-6，而U6＋含量分别为0.10436%和9310×10-6, U4＋相对U6＋含量要高，U4+/U6＋比值分别为1∶1.55和1∶1.33，铀以U4＋为主，U6＋含量较低。对矿化岩石和非矿岩石的化学全分析表明，矿化岩石碱性元素中的Na2O含量为11 89%～2139%，非矿岩石Na2O 为0.174%, Na2O富集程度达2155～3123倍，说明查查香卡矿床有钠交代铀矿床的一些特征，可以与北祁连冷龙岭碱交代铀成矿带进行对比，显示了成矿特征的相似性[4]。
图3 查查香卡铀矿床纵剖面图
1—碎裂闪长岩；2—闪长岩；3—正长岩；4—斜长角闪片岩；5—碎裂斜长角闪片岩；6—斜长角闪岩；7—碎裂斜长角闪岩；8—角闪岩；9—碎裂角闪岩；10—辉绿岩；11—铀矿体；12—块段编号
图4 查查香卡铀矿床8号勘探线剖面图
3主要成果和创新点
查查香卡铀矿床经历了两个阶段的勘查，即1969～1971年的初步勘查和2008～2012年的调查评价。前一阶段主要为地表调查和槽探揭露，认为：①矿化主要分布于构造挤压破碎带闪长碎裂岩及碎裂闪长岩中，富矿地段为碎裂岩或糜棱岩，岩石破碎强烈，片理化发育，肉红色长石和碳酸盐呈细脉或星点状，其分布多的地段矿化最好；②岩石蚀变较强，主要有绿泥石化、绿帘石化、阳起石化、碳酸盐化、钠黝帘石化、钠长石化和弱硅化等；③区内成矿构造主要表现为片理化、碎裂岩化及糜棱岩化带，走向为NW W，倾向NNE；④由地表到深部钍矿化逐渐变贫，而铀矿化则相反，由地表到深部逐渐矿化增强；⑤根据伴生元素分析结果认为324地区属多种稀有金属成矿区，铀、钍、铌、铈、镧在矿化岩石中的含量明显高于非矿化岩石。
2008年以来，核工业二〇三研究所通过在柴达木地区实施中国核工业地质局和中国地质调查局铀矿调查评价项目，对查查香卡铀矿床进行了调查评价和研究，主要成果如下：
3.1通过对查查香卡异常带中段槽探揭露和钻探查证，圈定了铀、铌、稀土矿（化）体，初步估算3341铀、铌资源量
根据槽探圈定了4个工业铀矿体、6个矿化体，钻孔在矿化带中部查证发现了4个工业铀矿孔、2个铀矿化孔，矿床接近中型规模。
在铀矿体及周围发现铌矿体3个、矿化体11个。发现了该区铌、稀土等多金属矿化，且与铀呈正相关关系（图5），铀与铌、稀土元素伴生。按铀矿床伴生组分综合利用，圈定了13个矿体，圈定5个稀土矿化体（主要为铈、镧、钇等）[5]。
图5 查查香卡铀矿床铀、铌、稀土元素含量关系
3.2大致查明查查香卡铀矿床区域构造位置，属于古隆起与蛇绿混杂岩带的结合部，为断裂破碎带与多期次岩浆活动区
矿床处于丁字口－埃姆尼克山－牦牛山新元古代－早古生代岩浆弧带与赛什腾－锡铁山－哇洪山新元古代－早古生代缝合带的结合部，夹持于泽日肯和陶力古隆起之间托莫尔日特蛇绿混杂岩带中。
异常带位于泽日肯隆起与托莫尔日特蛇绿混杂岩带之间，受加里东—印支多期次中酸性岩浆活动和北西西向断裂构造破碎带控制[5]。长约9km，宽约2km。
3.3大致查明查查香卡铀矿床的含铀建造奥陶－志留系滩涧山群为一套蛇绿混杂岩，具有轻稀土元素富集特征
奥陶－志留系滩涧山群以各类绿片岩为主，部分为角闪片岩类，少量为变余中基—基性火山岩、火山碎屑岩夹云母石英片岩、大理岩，局部夹含放射虫硅质岩、海绿石砂岩及细碧岩等。绿片岩、角闪片岩恢复原岩为一套海相喷发沉积的中基—基性火山碎屑岩。总体上该层位为一套蛇绿混杂岩。
滩涧山群斜长角闪片岩稀土元素分析结果显示（图6）：在完整岩石中，稀土配分型式为右倾型，富轻稀土，∑Ce＞∑Y；稀土元素总量与铀矿化程度成正比，即稀土元素总量越大，铀矿化强度越大[4]。
图6 查查香卡地区斜长角闪片岩稀土元素配分型式
3.4大致查明查查香卡地区构造－岩浆活动具多旋回特征，铀矿化主要受印支—燕山期岩浆热液活动控制
区内可划分出吕梁、四堡－晋宁、加里东、华力西、印支－燕山和喜马拉雅等6个构造－岩浆活动旋回（表2）。断裂破碎带及华力西期、印支－燕山期中酸性岩浆活动为区内Cu、Au等多金属、U等矿化提供了空间和物质来源，其中印支－燕山期是铀矿化的重要阶段。
3.5发现了查查香卡地区有铜、金多金属矿找矿线索
地表调查中，在查查香卡地区含铀异常破碎带中发现了Cu、Au等金属矿化，见有孔雀石、铜蓝和明金[5],Cu含量0.6825%, Au含量560×10-9。矿化与石英脉、强褐铁矿化蚀变带关系密切。激电测深剖面发现了2～3条激电异常，与地表发现的控矿构造蚀变带相吻合，向深部500m左右可能有较好的金属矿化。钻孔中发现了2层黄铁矿化层，埋深在330.5～342.5m、396.6～400.6m之间，厚度分别为12.0m、4.0m。
表2 查查香卡地区构造、岩浆作用、变质作用与成矿作用一览表
3.6总结提出了查查香卡铀矿床“三位一体”成矿模式
综上所述，查查香卡铀矿床位于托莫尔日特蛇绿混杂岩带与由达肯达坂群和吕梁－印支期基性－中酸性岩浆岩组成的泽日肯隆起的结合部位，赋存于奥陶－志留系滩涧山群变质火山岩，受北西西向的区域断裂构造破碎带控制。印支－燕山期的酸、碱性岩浆－热液活动是区内重要的铀成矿作用。
控制查查香卡铀矿床的断裂破碎带是区域性长期活动的断裂带，发育于奥陶－志留系滩涧山群中，形成于加里东末期，于华力西期又有活动。断裂带内基性、中性和酸性脉体发育，铀矿化与正长岩脉、石英脉等关系密切，反映出该断裂带既是铀成矿溶液的导矿构造，亦是铀成矿的容矿构造。与成矿构造相邻的达肯达坂群（铀含量（30～40）×10-6）和滩涧山群（铀含量（15～18）×10-6），在华力西期—印支期中酸性岩浆活动中，地层中的铀不断活化向断裂带中迁移、富集，多期次的岩浆－热液活动的叠加形成了铀矿床。
综合分析认为，查查香卡铀矿床形成于区域北西西向断裂构造破碎带、奥陶－志留系滩涧山群变质火山岩与华力西－印支期岩浆－热液联合控制的“三位一体”[5,6]成矿模式，具备形成铀多金属矿床的条件。
4结束语
已有的地质调查表明，查查香卡伽马异常带断续长9km，经对中段2km长的范围内铀资源评价，铀资源量接近中型规模，铌资源量340.6t（按铀伴生铌的资源量1052.6t）。铀矿体呈扁豆体产出，往深部有钍变贫、铀变富的趋势；平衡系数1～117，地表偏镭，地下基本平衡。初显良好的找矿潜力。
查查香卡矿化带尚发现有较好铜、金等矿化线索，其北侧有与之相平行的已发现赛坝沟、托莫尔日特金矿床的托莫尔日特构造带，因此查查香卡铀矿床具有形成铀多金属综合矿床的可能。
因此，查查香卡铀多金属矿有待开展预查和综合研究等工作，可望达到大型铀多金属矿床规模。
参考文献
[1]核工业西北182队第三分队．324地区1971年地质工作总结[R],1971.
[2]核工业西北182队第三分队．324地区1972年工作设计专题报告[R],1972.
[3]孙圭，赵致和．中国北西部铀矿地质[M]．核工业西北地质局，1998.
[4]傅成铭，权志高，周伟．青海查查香卡矿床铀、稀土元素矿化特征及成矿潜力分析[J]．铀矿地质，2011,27（2）:103-107.
[5]刘林，宋宪生，冯伟，等．青海省柴达木盆地北缘鱼卡—查查香卡地区1∶25万铀资源区域评价[R],2010.
[6]刘林，冯伟，陈擎，等．柴达木东北缘铀成矿地质条件及找矿远景研究[J]．东华理工大学学报（自然科学版），2013.
我国铀矿勘查的重大进展和突破进-—入新世纪以来新发现和探明的铀矿床实例
[作者简介]刘林，男，1964年生。1987年7月毕业于华东地质学院，研究员级高级工程师，现任核工业二〇三研究所青海地质勘查院院长。1987年至今在核工业二〇三研究所工作，长期从事铀矿地质工作，近期主要承担柴达木盆地北缘预测评价及查查香卡铀矿床勘查工作。
核心期刊分二种，一种是北大核心，也叫中文核心；另一种是科技核心，也叫统计源期刊；我查询了目录在目录中没有东华理工大学学报(自然科学版)，所以目前不是核心期刊，如果在新的核心期刊目录中有，你需要向杂志社求证。
1、参与科研项目l 砂岩型铀矿…，参与，国防基础科研，2008-2010；l 铀矿资源三维电阻率…,参与，国际科技合作，2007-2009；l 九瑞矿集区重要靶区地球物理勘探和深部找矿预测，参与，科技部科技计划项目，2011-至今l 江西省德兴张家坂-先告山地区铜金矿远景调查V8频谱激电（SIP）测量，主持，江西省地质调查研究院，20112、发表的论文[1]探地雷达空气中介质回波问题讨论[A]，国际电磁会议[C]，2009:250-253[2]直流电阻率三维观测方式分析[A]，地震出版社[C]，地震出版社，2010:683[3]探地雷达空气中回波的f-k域处理技术[A]，吉林大学出版社[C]，吉林大学出版社，2005:647[4] 用神经网络法预测砂岩型铀矿体的深度[J]，铀矿冶，2009(1)[5] 探地雷达空气回波特点及识别方法[J]，物探化探计算技术，2009(5)[6] 探地雷达空气回波特点与处理方法[J]，世界核地质科学，2009(3)[7] 层间氧化带型砂岩铀矿自然电位二维正演计算[J]，铀矿冶，2009(4)[8]砂岩型铀矿地表物探测量的模糊聚类分析[J]，东华理工学院学报，2006(1)[9] 探地雷达复信号分析的几点讨论[J]，物探化探计算技术，2006(2)[10]探地雷达数据可视化方法[J]，东华理工学院学报，2006(4)[11]探地雷达谱白化法[J]，东华理工大学学报(自然科学版)，2008(2)3、软件著作权l GPRPro v1.0探地雷达数据处理解释，登记号：2008SR38835；l ResInv2D v1.0直流电阻率点源二维正反演，登记号：2010SR052047；l ResInv3D v1.0直流电阻率点源三维正反演，登记号：2011SR026474。
是的，是核心期刊
匝道就是从一般公路进出高速公路的缓冲专用路道，单向的，这个是个人的理解，我再帮你查了下百度的解释。如下： 匝道【读音】zādào【翻译】circuit【解释】立交桥和高架路上下两条道路相连接的路段，也指高速公路与邻近的辅路相连接的路段。高架路的匝道，进口路和出口路是分开的，只能顺行，不准掉头，车辆错过了下匝道，就不能从上匝道下路，只能从下一个下匝道下路。 立交桥的匝道，也是按照设定的标志行驶，谁也不能各行其是。 【定义】匝道，又称引道，是工程学上的术语，通常是指一小段提供车辆进出主干线（高速公路、高架道路、桥梁及行车隧道等）与邻近的辅路，或其他主干线的陆桥/斜道/引线连接道，以及集散道等之附属接驳路段。它是构成道路交流道的主要交通建设。1.在t型(y型)互通立交中,通常将相交的主要道路定义为主线,相交次要道路定义为引线,连接引线与主线互通的线路称为匝道源自: card/1系统在互通立交设计绘图上的... 《公路》 2003年 刘秋江,张晓元2.交叉口所谓“匝道”,是指在立交处连接立交上、下道而设置的单车道单方向的转弯道路.匝道的曲线元也是由直线段、圆曲线段和缓和曲线段组成的源自: 全站仪在公路施工放样中的应用 《山东农业大学学报(自然科学版)》 2001年 邱健壮,解明东,王继芳3.在线路立体交叉部位,线路的连接都是由不同种的曲线线形连接而成、称为匝道.由于匝道形式多样、其中桩的坐标计算就非常困难.笔者通过实际操作.摸索出将曲线分成各曲线元的方法来计算中桩坐标、以此来解决匝道各种线型的中桩坐标计算问题源自: 道路匝道的坐标计算方法 《森林工程》 1999年 王炳升,姜英杰,董永江4.t交图1普遍采用的互通式立体交叉形j℃出入高速公路的连接道路称为匝道,m道的ⅲ入口处与高速公路连接的平顺性影响着车辆的安全行驶,这里也是高速公路瓦通设计的难点和重点
汽车路口串道指的是一个路口中存在两条以上的车道，但是车辆行驶时需要按照一定的顺序通过这些车道，而不能直接随意变道。通常的情况是，道路标线会给出具体的指示，比如在转弯时必须保持左转车道，而不能从右转车道直接切入。汽车路口串道设计的主要目的是为了保证车辆行驶安全，防止因为违反规定变道或者争抢车道而发生交通事故。而且，路口串道还可以保持交通流畅，避免车辆之间出现拥堵或者排队等待的情况。如果不按照路口串道要求行驶，即违反交通规则，较大几率会出现交通事故。因此，在遇到路口串道时，一定要遵守规定行驶，保持车辆的安全和行驶效率。
1. 谢俊，张金亮. 《剩余油描述与预测》，石油工业出版社，2003，北京2. 张金亮，谢俊. 《储层沉积相》（教材），石油工业出版社，2008，北京3. 张金亮、谢俊，《油田开发地质学》[M]，石油工业出版社，2011，北京4. 2009, Xie Jun, Yu Jiangtao, Lin Peng,Meng Ningning. The Application of 3D Geological Modeling on below 7 oil group in Henan Shuanghe oilfield, 2009 Second ISECS International Collopuium on Computing, Communication, Control,and Manegement(CCCM2009), IEEE Technology Management Council,USA, 2009.8, pp474-477【EI收录】5. 2009, Xie Jun, Lin Peng, Yu Jiangtao, Meng Ningning. The Model Building Study of Henan Shuanghe Oilfield Updip Area, 2009 Second ISECS International Collopuium on Computing, Communication, Control,and Manegement (CCCM2009), IEEE Technology Management Council,USA, 2009.8, pp478-481【EI收录】6. 2009, Xie Jun, Meng Ningning, Yu Jiangtao, and Lin Peng. The Research of Fine Geological Modeling for a Single Well in 7 Lower Unit of Shuanghe Oilfield, 2009 Second International Conference on Future Information Technology and Management Engineering(FITME2009), IEEE Computatioanal Intelligence Society, 2009.12, pp555-558【EI收录】7. 2009, Hui-zhen Liang, Jun Xie, Jiang-tao Yu, Ning-ning Meng. The BP Network Study of the Time Series Overrolling Modek for Forecasting. 2009 IITA International Conference on Services Science, Management and Engineering (SSME2009), IEEE Computatioanal Society, 2009.7, pp307-310【EI收录】8. 2009，梁会珍，谢俊，张金亮. 下二门油田中上层系剩余油成因及可动油定量分布研究. 西安石油大学学报，24（6）9. 2009, 武英利，谢俊. 英吉苏凹陷盖层评价. 《石油地质与工程》，23（3）10. 2008，谢俊，张金亮，梁会珍. 濮城油田沙二上1砂组末端扇沉积体系研究. 西南石油大学学报， 30（2）11. 2008，谢俊，张金亮，梁会珍. 塔里木盆地志留系柯坪塔格组沉积相与沉积模式研究. 西安石油大学学报， 23（2）12. 2008，谢俊，张金亮，梁会珍等. 濮城油田末端扇储层隔夹层成因及分布特征. 中国海洋大学学报， 38（4）：653－65613. 2008，谢俊，张金亮，刘宝珺. 英吉苏凹陷天然气藏成岩演化模式及成藏机理. 西安石油大学学报， 23（3）14. KYPAM和MO-4000聚合物性能对比.特种油气藏，2008 15（1）15. 2008，李伟，谢俊，高晓辉. 吐鲁番坳陷侏罗系泥岩超压特征及其对尤其成藏的控制. 石油勘探与开发， 35（1）【EI收录】16. 2008，鲁明辉，谢俊，王金凯. 濮城油田西区沙二上1气顶油藏剩余油研究及方案调整. 山东科技大学学报（自然科学版），27（4）17. 2008，司学强，张金亮，谢俊. 成岩圈闭对气藏的影响，天然气工业，28（6）：27-3018. 2008，Zhang Jinliang，Xie Jun, Research on Reservoir Heterogeneity, the Interlay and Insulating layers and Controlling Factors on the 2+3 layer of the upper Second Member in the Shahejie Formation, the West of the Pucheng Oilfield. Petroleum Science ，Vol. 5, No.2，pp135-144【SCI收录】19. 2008, Jun Xie, Lin Peng, Huizhen Liang, and Minghui Lu. The improved rapid convergence algorithm of the connecting rights in the BP network, 2008 International Symposium on Knowledge Acquisition and Modeling (KAM2008), IEEE Computatioanal Intelligence Society, 2008.12, pp819-822【EI收录】20. 2008, Jun Xie, Minghui Lu, Huizhen Liang, and Peng Lin. The research and Application of bp Network Tracking Model for Forecasting Oil Well Yield, 2008 IEEE Pacific-Asia Workshop on Computational Intelligence and Industrial Application (PACIIA2008), IEEE Computatioanal Intelligence Society, 2008.12, pp18-22【EI收录】21. 2007, Xie Jun, Zhang Wei. Discussion on Particularities of Reservoir Heterogeneity in Later Development of Reservoirs. Journal of China University of Geosciences，2007.8，Vol.18【SCI收录】22. 2007，Si Xueqiang, Xie Jun, Zhang Jinliang. The Mechanism of Secondary Pore of the Beach-Bar Sandstones and Prediction of Favorable Region in the Upper Es4 of the Palaeogene,Boxing Sag. Journal of China University of Geosciences，2007.8，Vol.18【SCI收录】23. 2007，谢俊，张金亮，梁会珍. 适于高温高盐油藏的聚合物性能指标评价. 中国海洋大学学报， 37（2）24. 2007，谢俊，王金凯，刘法平. 胜坨油田聚合物驱转后续水驱技术对策研究. 山东科技大学学报（自然科学版），26（2）25. 2007，谢俊，刘法平. 用于提高油田采收率的聚合物降解与防护. 山东科技大学学报（自然科学版），26（5）26. 2007，谢俊，张金亮，梁会珍等. 胜坨油田聚合物驱油效果的影响因素分析. 复杂油气田，16（3）27. 2007，孙绪宝，谢俊. 基于频域方法混叠信噪比的研究. 电波科学学报，22（增刊）28. 2006，王金凯，谢俊. 油层韵律特性对驱油效果以及剩余油分布的影响. 山东科技大学学报（自然科学版），25（增刊）29. 2005，苏维，谢俊. 法哈牛油田法21断块剩余油分布研究. 山东科技大学学报（自然科学版），24（4）30. 2005，谢俊. 我国海上油田强化采油提高采收率配套技术研究. 我国近海油气勘探开发高技术发展研讨会文集，中国北京31. 2005，张金亮，谢俊. 我国高含水油田剩余油描述技术与预测方法. 我国近海油气勘探开发高技术发展研讨会文集，中国北京32. 2005，谢俊，张金亮，梁会珍. 渤海稠油油藏开采技术对策研究. 中国深水油气开发工程高技术论坛，中国上海33. 2005，张金亮，谢俊. 我国海上油田强化采油配套技术研究. 中国深水油气开发工程高技术论坛，中国上海34. 2005，张金亮，谢俊，梁会珍. 我国海洋石油剩余油描述技术研究. 中国深水油气开发工程高技术论坛，中国上海35. 2004，谢俊、张金亮. H31断块东二段水淹特征及剩余油分布研究. 山东科技大学学报(自然科学版)，23(1)：108-111.36. 2004，谢俊、张金亮、梁会珍. H26断块东二段水淹特征及剩余油分布研究. 河南石油，18(1)：25-28.37. 2004，谢俊，武英利，梁会珍等. 应用流动单元寻找油田有利挖潜区的新思路. 西北地质， 37（4）38. 2004，谢俊、张金亮、梁会珍. H31断块稠油油藏常规注水开发效果评价，全国矿井地质会议论文集，地质出版社，北京：357-36139. 2004，刘宝珺，谢俊，张金亮. 我国剩余油技术研究现状与进展. 西北地质， 37（4）40. 2003，谢俊，梁会珍. 黄原胶作为油田驱替剂的性能研究. 矿物岩石， 23（2）【EI收录】41. 2003，谢俊、张金亮. R11井山头油藏改建地下储气库方案设计. 中国海洋大学学报(自然科学版)，33(3)：413-419.42. 2003，谢俊、张金亮. PS剂驱油矿场先导性试验研究. 中国海洋大学学报(自然科学版)，33(1)：095-100.43. 2003，谢俊. 微积分法计算抽油井井底压力. 山东科技大学学报（自然科学版），22（1）44. 2003，谢俊、张金亮、梁会珍. 法哈牛油田法21断块剩余油分布规律研究. 特种油气藏，10(增刊)：53-55.45. 2003，谢俊、张金亮. 法21断块储层非均质性研究. 山东科技大学学报(自然科学版)，22(3)：10-13.46. 2003，谢俊，张金亮，梁会珍. 法哈牛油田法21断块低效开采技术对策研究. 复杂油气田，12（4）47. 2003，谢俊、张金亮、梁会珍. H31断块稠油油藏常规注水开发效果评价. 内蒙古石油化工，29(3)：136-138.48. 2002，谢俊，张曙光，梁会珍. PS剂驱油提高采收率室内评价研究. 矿物岩石，22（4）【EI收录】49. 2002，张曙光，谢俊. 基于广义Копытов模型及乙型水驱曲线的综合预测法. 矿物岩石， 22（4）【EI收录】50. 2002，谢俊，张金亮. 黑液木质素提高采收率现场试验研究. 《油藏地质与油藏表征》论文集，西安地图出版社，陕西西安233-27451. 2002，谢俊，张金亮. 京11断块提高采收率方法筛选及室内实验研究. 《油藏地质与油藏表征》论文集，西安地图出版社，陕西西安：275-30552. 2002，张金亮，谢俊等. 法哈牛油田法21断块沙三段储层地质. 《油藏地质与油藏表征》论文集，西安地图出版社，陕西西安：139-17653. 2002，张金亮，常象春，谢俊. 外海河油田东营组储层地质特征. 《油藏地质与油藏表征》论文集，西安地图出版社，陕西西安：198-23254. 2000，梁会珍，谢俊，李建召等. 微生物采油技术在华北油田的应用. 西安石油学院学报（自然科学版），15（4）【EI收录】55. 1998，谢俊. 剩余油饱和度平面分布方法研究及应用. 西安石油学院学报（自然科学版），13（4）56. 1998，谢俊、吴行才、宋富霞等. 华北油田提高采收率技术评价及潜力分布，华北油田勘探开发科技文选，石油工业出版社，北京：177-18657. 1997，李论，谢俊，薛中天. 应用人工神经网络预测油井产油量和产水量. 西安石油学院学报（自然科学版），12（1）
是的，详见封底截图
复合影响因子：0.454 综合影响因子：0.288
主办： 山东科技大学周期： 双月出版地：山东省青岛市语种： 中文;开本： 大16开ISSN： 1672-3767CN： 37-1357/N历史沿革：现用刊名：山东科技大学学报(自然科学版)曾用刊名：山东矿业学院学报创刊时间：1979该刊被以下数据库收录：CA 化学文摘(美)(2011)Pж(AJ) 文摘杂志(俄)(2011)核心期刊：中文核心期刊(1992)
山东科技大学学报自然科学版是比较好投的。山东科技大学学报出版规范，编委们对知识的理解全面、考虑周到，对待投稿的作者耐心讲解，提出改进的意见。审稿较快，版面费较低。
综合排名第7，依次是山东大学，中国海洋大学，中国石油大学，山东农业大学，青岛大学，山东师范大学，山东科技大学，青岛科技大学。理工科排名只有一个可参考（武书连的），1山东大学，2中国石油大学华东3青岛科技大学4山东科技大学5山东理工大学，这个5个学校进全国百强了。