3月11日，第十八次李四光地质科学奖科研奖获得者、中国地质科学院矿产资源研究所（以下简称“资源所”）原副总地质师肖克炎研究员新申报的战略矿产项目“战略性矿产资源大数据综合信息预测与找矿勘查示范”通过评审。在过去的40年间，从数学地质到数字化矿产预测，从信息化大数据分析再到如今的智能化资源勘查，肖克炎始终围绕国家需求、紧盯科技前沿刻苦钻研，奋战在中国地质行业数字化、智能化研究的前沿，成为这一进程的缩影和地质科学工作者的榜样。

肖克炎研究员获李四光地质科学奖科研奖

理想：将“计算机+矿产预测”作为“一生的事业”

1980年，16岁的肖克炎考入湘潭矿业学院地质系煤田地质与勘探专业。如果不出意外，他的人生轨迹也会像他的许多大学同学那样，毕业进入煤炭行业，一生从事煤田地质方面的研究和实践。

不过，20世纪80年代初，改革开放的浪潮席卷中国大地，对“科学技术现代化”的渴望蓬勃生长。那时的湖南，一件大事正在发生。1983年12月22日，中国第一台每秒钟运算一亿次以上的“银河”巨型计算机，在位于长沙的国防科技大学计算机研究所问世。“银河”对青年肖克炎的影响是巨大的，从小头脑活跃、热爱数学的肖克炎看到了计算机在煤田矿产资源领域的巨大应用前景。

1983年年底，他为自己的实习论文圈定了两个方向——矿产资源趋势分析、计算机应用。学校没有计算机怎么办？在学校的支持下，肖克炎在核工业二三〇研究所（原长沙铀矿地质研究所）借了一台进口巨型计算机，它几乎占据了整个屋子。有了这个“宝贝”，肖克炎花费一个星期时间，对湖南湘中盆地煤田沉积谷地的矿产资源进行了5次趋势分析，圆满完成了实习论文。

这次成功尝试，点燃了他以数学研究地质的兴趣和志向。当然，那时的他，还不知道未来计算机技术将飞速发展，更不知道自己将以此为起点，与“计算机+矿产”打一辈子的交道。

1984年，肖克炎考取长春地质学院数学地质专业硕士研究生，师从著名数学地质专家王世称教授，并继续攻读该专业的博士学位。

数学地质是用数学分析及建模的手段，对地质问题进行规律性的研究。然而，把复杂的地质现象进行定量化研究实非易事，不仅需要把矿床成因、地球探测信息技术、计算机等多学科的知识进行交叉、融合，更需要具有宏观的科学视野和创造性的数理推算能力。在学习和实践中，肖克炎不断探索多元统计分析方法的应用、金铜异常区复合模型的建立以及综合信息矿产预测等领域，像海绵吸水一样，在课堂授课、实习项目和国内外资料中汲取知识和经验，并最终将自己的专业方向确定为矿产预测。这也是他坚定不移的选择，为他未来的研究奠定了坚实的基础。

在长春地质学院学习期间，肖克炎展现出了非凡的勤奋和才华。他发表了数十篇论文，远超学校规定的数量。期间，他还参与了由王世称教授主持的《综合信息矿产预测理论与方法体系》项目，并获得了国家科学技术进步奖二等奖。

1991年年底，肖克炎博士毕业后不久，中国地质大学开始招收博士后。在博士后流动站和长春地质学院之间，为了自己的梦想，他选择了前者，来到武汉，成为中国地质大学“矿产普查与勘探”博士后。师从我国著名数学地质学家、原地大总校校长赵鹏大教授，进一步深化对数学地质的研究和应用。1993年博士后出站，已晋升为副研究员的肖克炎，来到资源所成矿远景区划研究室（简称“区划室”）工作。

直到现在，肖克炎都很庆幸自己的“运气”，不仅能师从于数学地质领域有着“南赵北王”之称的两位泰斗，还能在中国矿产资源顶级科研院所找到专业最对口的工作岗位。当然，这份幸运并不是偶然的，而是他坚守初心、不懈努力的结果。 

肖克炎研究员野外工作

使命：“成矿区划”→“潜力评价”→“深部预测”

踏上工作岗位，正直20世纪90年代，全国地质工作陷入低谷期，而肖克炎的“好运”却并未离他而去。

肖克炎刚刚来到资源所，就赶上了一个大项目——全国第二轮成矿区划工作，并在其中从事成矿远景区划计算机化及GIS应用研究。巧合的是，肖克炎本科和硕士期间，也曾参与过全国第一轮成矿远景区划的部分工作。

这个项目由时任中国地质地科院院长和地质矿产部总工程师的陈毓川组织领导。此时，经过数十年的找矿，不少地区已转入隐伏矿床、盲矿床预测寻找和评价阶段，找矿难度大，必须筛选出具有较好成矿远景的重点区域，形成“以点带面、点面结合”的找矿突破。同时，当时的地质科技与成矿理论较之改革开放前有了明显的进步，急需使用新理论、新技术综合研究最新的资料，指导找矿。

陈毓川找到肖克炎，阐述了中国成矿体系和第二轮成矿区划工作的目标，以及应用计算机技术对区域成矿规律进行研究的设想。他强调，一定要注重开展全国和区域性物化资料综合研究，进行矿产预测和评价；要瞄准世界先进水平加快发展勘查技术，加快现代计算机技术的数据处理和成图技术研究，实现“地矿工作信息化”。

肖克炎始终认为，陈毓川院士是自己重要的伯乐和导师。通过陈毓川指导与项目实践，肖克炎的学术视野愈发开阔明朗，牢牢把成矿系列理论作为计算机技术应用的基础和出发点，成为他所坚守的、不同于其他矿产资源数字化、智能化研究团队的鲜明特点。

在第二轮成矿区划工作中，肖克炎作为主要骨干，带领团队逐步建立了以成矿系列理论为指导、以GIS计算机技术为支撑的综合信息矿产预测方法体系，开展了不同比例尺矿产资源方法学比较研究，提出了预测底图编图综合解释模型，完善了综合信息找矿评价模型的研究内容方法，编制了首张中国矿产成矿系列图和跨世纪找矿工作部署图，建立了全国成矿远景区划数据库。

时间紧、任务重。有时候，“活儿”要得急，肖克炎经常要加班到下半夜。好在，辛勤的汗水换来了扎实的成绩：第二轮成矿区划成果“中国成矿体系和区域成矿评价”，在2004年获得国家科学技术进步奖二等奖，肖克炎是主要完成人之一。

大项目催生大人才。1996年，肖克炎被破格晋升为研究员，1997年开始担任资源所成矿远景区划室总工，2000年，被任命为区划室主任。2006年，原国土资源部部署全国矿产资源潜力评价国情调查工作。这次，作为学科带头人的肖克炎，毫无悬念地成为其中“全国重要矿产总量预测”项目的总负责。

从2006年到2013年，肖克炎上满了“发条”，瞄准“总量预测”目标任务，创立了矿床模型综合地质信息矿产预测方法；研发出地质专家系统预测与综合信息预测的预测模型，总结了以计算机信息技术为手段，运用矿床模型地质参数法进行资源潜力估算的预测工作方法流程；首次提出区域预测综合信息编图的综合解释模型、矿产预测类型、矿产预测方法类型、最小预测区等概念，创新发展了矿床模型地质参数资源量估算方法；建立了全国矿产资源潜力预测评价成果数据库，为国家资源决策提供了权威国情数据。

肖克炎研究员在野外与国际知名专家共同野外查证

矿产资源潜力评价预测成果促进了地质工作顺应时代的发展：广泛应用于我国“十三五”各类矿产勘查规划部署，厘定了我国26个重要找矿部署区带，优选110余处国家重要整装勘查区和矿集区；支撑国家找矿突破战略行动和矿业企业的矿产勘查，推动新发现多处大中型矿产地，实现找矿工作重大突破。即使是现在，其成果仍在发挥着重要作用。

近年来，矿产勘查逐渐向深部转移，肖克炎智能勘查技术的研发目标也“由浅入深”，在科技部深地资源勘查开采专项“深部资源预测系统技术研究与示范”项目的支持下，围绕深部资源预测需求，系统开展了深部矿产资源预测理论研究、关键技术、平台系统研究，突破了成矿空间三维结构重建机制、深部矿化定位机制及深部预测途径等关键技术瓶颈，首次建立了找矿模型—三维建模—定量预测三元大数据深部矿产资源预测理论方法，自主研制了一套具有自主知识产权的深部综合信息预测评价平台系统，创新研发了深部成矿构造、地球物理、地球化学和地质异常的时空结构解析方法与深层次信息提取新技术，解决了矿集区深部成矿空间三维结构重建、成矿地质异常空间重构以及矿化空间定位等技术难题，推动了我国深部找矿、三维预测领域的技术进步。 

深部找矿勘查示范基地

一个个走在科技创新“风口浪尖”的大项目，有效地支撑服务了国家矿产资源安全保障和资源管理，引领了我国矿产资源潜力定量评价与预测科学进步。它们不仅赢得了近十个国家和省部级奖项，也为肖克炎带来了诸多荣誉：俄罗斯科学院与工程院外籍院士、享受国务院政府特殊津贴，以及部百名跨世纪科技人才、原国土资源部先进个人、自然资源部科技创新人才、湘江学者、楚天学者、珠峰人才等称号，而2023年获得的第十八次李四光地质科学奖科研奖，更是对他在地质科研方面突出贡献的充分肯定。

与此同时，肖克炎不仅在个人学术研究上取得了显著成就，更在团队建设方面展现出卓越的领导力。

在全国矿产资源潜力评价和国家重点研发计划等项目的实施过程中，他建立了一支集成矿规律、潜力评价、数学建模、软件工具研发等全链条、多层次的高水平科研团队，形成了全国领先的大数据智能找矿“理论—方法—系统平台”完全自主知识产权的开创性成果。

在他的指导下，团队突破学科壁垒，出色完成了全国25种重要矿产资源2千米以浅资源量预测研究，圈定了各类找矿靶区4万余处、成矿远景区2000余处、重要矿集区240余处，并科学估算了不同深度预测资源量。在山东莱州—招远金矿整装勘查区、四川会理—会东矿集区、湖南—贵州锰矿整装勘查区、贵州贞丰—普安金矿整装勘查区、甘肃玛曲—合作金矿整装勘查区、西藏冈底斯成矿带等地，开展了深部（3000 米以浅）三维预测评价及找矿预测示范，取得了多项找矿突破。由此，肖克炎科研团队在国际SCI期刊上发表了论文106篇，获得发明专利24项，软件著作权27项，6人获得省部级及以上荣誉称号；培养了博士后5人、博士50余人、硕士100余人，为相关行业、企业输送了大量高素质人才。同时，举办各类培训班40余次，培训了行业、企业技术骨干近千人，他本人还在商务部举办的“一带一路”非洲预测方法培训班中多次授课，近千名国外地质专家聆听相关课程。这些工作和成果在学术界产生了广泛的影响，为矿产预测定量化发展作出了重要贡献。

未来：发展人工智能预测技术是“创新的高地”

所谓矿产预测，就是在现代地质成矿理论指导下，对预测区的地质、地球物理、地球化学、重砂和遥感等各种成矿信息进行提取解译，并使用先进的多元统计方法和计算机技术进行综合，从而确定成矿有利地段和靶区，并最终估算其矿床数和资源量。

在如今的矿产预测领域，很少有人不知道MRAS（矿产资源综合信息评价系统）和MinExplorer（探矿者软件）。这套由肖克炎带领团队在全国矿产资源潜力评价项目实施期间自主开发并拥有全部核心技术与知识产权的矿产资源预测评价系统平台，打破了国外在此领域的技术壁垒，填补了我国矿产预测信息化领域的空白，已被全国32省（区、市）的上千家地质调查单位全面应用。 

探矿者软件(MRAS3.0)赠送仪式

科学的探索从来都是循序渐进、厚积薄发。

早在1996年，肖克炎申请承担的“基于GIS平台矿产资源评价辅助决策系统”，在国内首次提出将传统矿产资源评价与现代地学信息技术相结合的研究思路并研制了MRAS软件系统；2002年，肖克炎承担国家863项目“地学空间信息三维可视化系统”，通过基于GIS的资源评价分布式网格计算，开始了对“探矿者”软件核心技术的研发与积累；2006年~2013年，“探矿者”在全国危机矿山找矿专项中推广应用，形成了矿床勘探3DEM储量估算系统；之后，“探矿者”升级为v3.0版本，实现了对三维数字矿床模型及深部矿产资源的立体预测评价。

MRAS打通了矿产资源评价的GIS工作流程，实现了多源信息变量自动赋值及预测单元自动划分，彻底地改变了以往在预测单元选择、变量研究过程中进行矿产资源统计预测的费时费力的工作环境；而探矿者软件，则是一套适合地质矿产勘查、三维可视化智能分析、三维预测评价系统，其自主研发的三维矿产勘查评价智能地质制图、三维地质建模、三维可视化等底层核心技术，实现了从二维平面预测向三维空间立体预测转化，也使传统统计预测要几个月完成的工作，缩短到几天就可以完成。

这些，都颠覆性地提升了地质工作中矿产预测的工作模式、准确度以及经济性。

探矿者地质矿产勘查软件系统（MRAS v3.0）

30多年对矿产预测数字技术的研发与应用，让肖克炎深深地体会到了“大数据驱动”已成为新的科学范式，而对前沿科学极为敏感的他则清晰地意识到，发展新一代人工智能预测技术是当前资源评价领域创新高地，地质找矿工作智能化已是“奔腾的江河”，必然势不可挡、必然汹涌向前。

如今，在“人工智能+”找矿的大趋势中，肖克炎又针对战略性矿产资源开始了新的研究探索——通过战略矿产项目“战略性矿产资源大数据综合信息预测与找矿勘查示范”，针对锂、钴、镍等11个战略性关键矿产资源，创新大数据预测评价基本理论、方法流程、预测模型及成果表达等，研发智能地质预测模型和大数据精准预测技术，创新资源评价知识驱动和数据驱动双引擎，突破智能地质模型预测计算机关键技术，形成新一代智能战略性矿产评价软件系统，解决资源评价软件卡脖子难题。

又是一场硬仗。在今后的4年中，肖克炎将带领项目团队，完成战略性矿产超常富集机理及成矿模型研究、战略性矿产大数据综合信息提取与智能预测方法技术研究、战略性矿产数字矿床模型及预测评价系统平台研发、战略性矿产重点成矿区带资源潜力定量评价、战略性矿产重点矿集区三维预测与找矿勘查示范等课题，继续推进矿产资源大数据智能预测评价水平，为国家提供更多的找矿新靶区，尽快找大矿、找好矿、找急需的矿提供重要支撑。

熟悉肖克炎的人都知道，他有着地质人的豪迈大气，也有着数学家的精细严谨，既有身为学者的“勇攀高峰、敢为人先”，又有作为教师的“春风化雨、玉汝于成”。30多年来，被同事、学生称为“拼命三郎”的他，始终把李四光等地质先辈的科学家精神视为榜样和指引，以服务国家需求为己任，勇于创新、勤于实践，全身心地投入到了“推动矿产预测理论技术进步”的事业中，而未来，面对国家对能源矿产资源的全新需求和日新月异数字技术的不断召唤，他以科技创新驱动大数据智能预测方法研发的步伐也将更加坚定。

 

吉梨页油1井放喷点火现场

 

该成果由中国地质调查局油气资源调查中心牵头，中国石化东北油气分公司、中国石油渤海钻探工程公司共同参与完成。牵头完成人为徐兴友、张君峰。其主要进展及创新：

通过项目实施，创新形成3项陆相页岩气富集地质新认识：一是查明了松辽盆地南部火山活动期混积页岩形成机理；二是建立了松辽盆地南部纹层型混积页岩和夹层型硅质页岩两种页岩气富集模式；三是揭示了松辽盆地南部断陷湖盆陆相页岩气成藏规律，丰富和发展了中国陆相断陷湖盆页岩气地质理论。

集成创新了2项陆相页岩气关键工程技术：一是基于地质-地球物理建模，建立双向介质烃类检测陆相页岩气甜点预测模型，有效指导了陆相页岩气甜点优选；二是形成非均质性常压页岩气储层增能体积压裂和控压返排试气技术，实现了常压陆相页岩气高产。

项目圈定梨树断陷沙河子组陆相页岩核心区面积1286平方千米，应用体积法初步估算核心区页岩气地质资源量5558亿立方米，有望推动松辽盆地2.7万平方千米有利区陆相页岩气勘探开发，成为中国东北地区现实的清洁能源接替领域。

目前，深部矿产资源勘查已成为我国未来矿产勘查的重要方向。近年来，随着全国危机矿山接替资源找矿和老矿山接替资源找矿等项目的实施，深部找矿勘查工作取得了一系列重要进展。这些新发现和新成果，不仅颠覆了对已有矿床成矿规律（成矿系统的发育深度和不同类型矿床的分带、叠合规律等）和矿床成因的认识，推动了科学研究“向深部进军”，还为深部矿产资源勘查提供了重要依据。本文选择近年来深部找矿取得重要突破的矿山，介绍其深部勘查的新发现及意义，以期引起业界对矿床系统的理解和再认识，更好地服务于深部矿产资源勘查工作。同时，结合已有的深部找矿经验，对下一步深部矿产勘查工作提出建议。

深部找矿新发现促进矿床成因新认识

1. 西藏罗布莎铬铁矿

西藏罗布莎铬铁矿床是目前国内规模最大、研究和勘探程度相对最高的铬铁矿床，但经多年开采，矿山保有资源储量严重不足。

自2006年以来，通过开展矿区含矿构造岩相带与矿体空间分布规律研究，同时配合重磁电综合解释，在矿床深部发现厚大隐伏矿体，实现我国铬铁矿找矿的重大突破。钻探资料显示，钻孔ZKWT02在孔深286.3～341.1米处见累计视厚度达46.28米的4层矿体，Cr2O3含量50%以上；钻孔ZK西2004于孔深334.5～426.9米处见4层矿，累计视厚度46.1米，Cr2O3含量约50%，发现的Cr-80单矿体资源量达115万吨，成为目前国内发现单体规模最大的铬铁矿矿体。此外，通过“千米钻”深孔ZK东04的实施，在孔深679.5米处见蚀变破碎带中的浸染状铬铁矿，验证了深部有规模较大的中央含矿构造岩相带的存在。

通过深部找矿工作的实施和深部厚大铬铁矿体的发现，对该矿床的成矿规律和成矿模型取得了新的认识，打开了找矿空间。进一步验证了中央含矿构造岩相带对于深部铬铁矿体的控制意义，查明了豆荚状铬铁矿体主要赋存于罗布莎超基性岩体中北部的斜辉橄岩岩相带中，浸染状铬铁矿体赋存于岩体北部纯橄榄岩岩相带中；提出豆荚状铬铁矿形成的多阶段成因模式；矿体表现出成带分布、成群出现、分段集中的特点，与层状铬铁矿的形成机制完全不同，可能是含铬铁矿残浆受构造作用强力侵位的结果；建立了以西藏罗布莎铬铁矿为典型的“板块踫撞蛇绿混杂岩带超镁铁质岩岩相构造带控矿的找矿预测地质模型”。总之，深部找矿工作的开展为罗布莎铬铁矿以及国内同类型铬铁矿床下一步找矿勘查工作提供了新的思路。

2. 广东大宝山铁铜钼多金属矿

大宝山铁铜钼多金属矿床是一座以铁、铜、铅、锌、硫、钼为主的大型矿床。广东省地质局705地质大队曾于1958～1961年完成普查—勘探工作，并提交了相当于122b的铜多金属储量79.96万吨（品位0.86%），铅金属储量31.36万吨（品位1.77%），锌金属储量84.97万吨（品位4.44%），铁矿石储量8034万吨（品位49.69%），以及钼金属资源量（2M11，包括外围）2.26万吨。大宝山矿床经过半个多世纪开采，资源几尽枯竭，属严重危机矿山。

经全国危机矿山接替资源找矿（2006～2009年）和老矿山接替资源找矿（2014年）两轮找矿勘查工作，在大宝山矿区深部分别探获了一处大型规模的斑岩型钼矿体和一处厚大的斑岩型铜硫矿体，取得了重大找矿突破。其中，斑岩型钼矿体主要呈似层状、透镜状产于花岗闪长斑岩和次英安斑岩中，受岩体和构造联合控制。矿体围岩蚀变分带明显，从浅部至深部分别为粘土化带→青磐岩化带→绢英岩化带→钾化带（如图），显示出典型的斑岩型成矿特点。斑岩型铜硫矿体主要赋存于泥盆系棋梓桥组灰岩、粉砂岩和次英安斑岩中，矿石类型以黄铁矿型铜矿石、磁黄铁矿型铜矿石和黄铁矿型硫矿石为主。由此，大宝山铁铜钼多金属矿床的矿体空间分布可分为深部的花岗斑岩型钼（钨）、铜硫矿、矽卡岩型钼（钨）矿、浅部的脉状、似层状铜铅锌矿以及远端的菱铁矿。

关于大宝山铁铜钼多金属矿床的成因主要有3种观点，即与岩浆期后热液有关、与海相或陆相火山热液有关以及与沉积改造作用有关。大宝山深部斑岩型钼矿体和铜硫矿体的发现进一步改变了以往对该矿床成因的认识，为下一步找矿工作提供了重要线索和依据。祝新友等研究认为，大宝山矿床深部的斑岩型、矽卡岩钼钨铜矿化、浅部的脉状铜铅锌矿化以及远端的菱铁矿化构成了斑岩型或岩浆期后热液型成矿系统。应立娟等对层状铜硫矿体研究后认为，大宝山矿区及外围成矿作用与岩浆热液有关的，属同一成矿系统演化的产物。戴塔根等指出，大宝山多金属矿床斑岩型、矽卡岩型矿体与岩浆热液有关，而似层状铜铅锌矿体则与海底火山热液活动有关，并受后期岩浆热液的叠加改造。

大宝山矿床37线剖面图

3.江苏栖霞山铅锌矿

江苏栖霞山铅锌多金属矿床位于南京市东郊的栖霞镇境内，处于长江中下游成矿带宁镇矿集区。2012年～2014年，地质矿产调查评价专项实施“江苏省南京市栖霞山铅锌矿接替资源勘查”项目，对矿床开展了深部找矿工作。在研究石炭系高骊山组与黄龙组硅/钙面控矿和脉状矿体侧伏规律的基础上，项目利用坑内钻追索和控制了1号主矿体在走向和倾向上的延伸，并新发现多层厚大铅锌矿体。尤其值得注意的是，在钻孔KK4603和KK4201中均发现金、银、铜矿化。其中，钻孔KK4201中银、金、铜平均品位分别为金 2.19克/吨、银 319.73克/吨、铜0.49%。经估算，本次工作新增333+334资源量：铅锌58.13万吨，共（伴）生金7.64吨、银1113吨、铜1.53万吨。

栖霞山矿床是华东地区目前已发现规模最大的铅锌多金属矿床，前人开展了大量找矿研究工作。然而，关于该矿床的成因一直存在争议，即矿床的形成与岩浆热液成矿作用有关，还是与热水沉积成矿作用有关？此次通过深部找矿勘查工作，在栖霞山矿床深部发现银-金-铜矿化，一方面暗示矿床深部巨大的找矿前景，另一方面为矿床成因研究提供了新的线索和证据。根据这些深部矿化信息，张明超等对矿床开展了流体包裹体显微测温及氢-氧-硫-铅研究，指出栖霞山铅锌矿的成矿流体和成矿物质主要来源于岩浆，岩浆水和大气水的混合是导致矿质沉淀的重要因素。因此，栖霞山矿床深部勘查的新成果不仅促进了矿床成因方面的研究工作，取得了新认识，而且更为下一步找矿工作指明了方向。

4.云南澜沧铅矿

云南省澜沧老厂铅锌银多金属矿床一直被认为是以火山岩为容矿建造或与火山岩有关的块状硫化物矿床（VMS）。2006年实施矿山深部找矿工作以来，在矿区继ZK153101钻孔揭露了钼矿化带后，ZK14824、ZK14827钻孔均揭露了深部的钼矿(化)体。其中，ZK14827钻孔孔深1417.05米，圈定的工业矿体总长度696.25米，平均品位0.068%，矿体主要产于深部隐伏斑岩体（喜马拉雅期）和接触带内。深部隐伏斑岩体和岩浆型钼矿的发现引发了学者们对该矿床成因模式的争议。斑岩型或岩浆型成矿系统的存在已是事实，如果考虑前人的石炭纪的块状硫化物成矿系统，则应是2个成矿系统的叠加。

朱裕生等强调“成矿系列”的“缺位找矿”理论，在多年前建立了“同位成矿”理论。在前人的理论认识基础上，可以对上述新发现给予验证和深化。韩润生等认为，澜沧老厂矿床是2套主要成矿系统、多种有利成矿要素有机耦合的结果，具多成矿动力学环境、大时间间隔和多类型成矿作用同位叠加的显著特征，是多种构造环境演变过程中形成的火山喷流沉积+隐伏斑岩热液成因为主的矿床，并将其成矿模式归纳为“双成矿系统同位叠加模式”。

5.内蒙古毕力赫金矿

早期勘查和研究认为，毕力赫金矿床赋存于晚侏罗世玛尼图组火山、次火山、火山碎屑岩系和隐爆角砾岩中，从基性的玄武岩到中酸性的流纹岩构成一套含金火山建造。其中的蚀变火山岩、次火山岩、角砾岩是主要的含金地质体，又是赋矿的直接围岩。主要岩性有安山质角砾熔岩、安山玢岩质角砾熔岩、闪长玢岩、热液角砾岩、蚀变玄武质角砾岩等。矿体的形态特征也表明其与火山机构关系密切。卿敏等通过对矿区深部、外围勘查和典型矿床解剖研究后认为，毕力赫金矿的成矿作用与海西期岩浆活动有关，成矿流体的性质、矿体空间产出与浅成花岗闪长斑岩关系密切，并在此基础上建立了矿田“斑岩-构造蚀变岩-浅成石英脉型矿化类型一体化”成矿模式。

深部找矿新发现开拓区域找矿新思路 

1. 内蒙古维拉斯托锡多金属矿

维拉斯托锡多金属矿床位于大兴安岭南段，曾作为大型铜锌银多金属矿床进行开发和找矿工作。2013年，地质人员在勘查过程中发现矿区及外围的黑云斜长片麻岩中发育有含锡石英脉，进而推测深部可能存在隐伏矿体。2014年地质调查项目实施的老矿山找矿项目，在矿山深部发现石英脉型锡钨锌铜钼矿体、隐爆角砾岩型锡铜锌矿体和浸染状锡锌铷铌钽矿体。其中，隐爆角砾岩型矿体位于岩体顶部外接触带，呈筒状，角砾成分为主要为黑云斜长片麻岩，胶结物以石英、黑钨矿和锡石为主；石英脉型矿体赋存于下古生界锡林郭勒杂岩和石炭纪石英闪长岩中，总体倾向北东，矿物组合为石英、锡石、闪锌矿和黑云母等；浸染状锡锌铷铌钽矿体赋存于云英岩化和天河石化的岩体顶部，矿物组合主要为锡石、闪锌矿和黑钨矿等。项目已初步查明矿体30余条，其中，1号矿体为主矿体，长大于700米、宽大于300米，总体呈舒缓波状；矿体厚度（视厚）平均5.15米，Sn品位平均为0.89%；共提交锡金属量8.1万吨（333及以上），实现了大兴安岭南段锡矿找矿的重大突破。

 

矿床三维空间结构模型示意图

大兴安岭南段是我国重要的铅锌、银、铜、锡多金属成矿带。按照成矿物质来源和矿床组合，区内划分出3个成矿带，其中林西-甘珠尔庙一带以锡多金属成矿作用为主，而锡林浩特-霍林郭勒地区则以银铅锌成矿作用为主。因此，长期以来，大兴安岭南段的锡矿找矿工作一直以林西-甘珠尔庙地区为重点，但多年来始终未取得较大突破。维拉斯托矿床位于以银铅锌成矿作用为主的锡林浩特-霍林郭勒地区，其与拜仁达坝矿床一起曾作为大型银铅锌矿床开展找矿勘查和科学研究工作。然而，此次维拉斯托矿床深部锡多金属矿体的发现，显示出传统的银铅锌成矿区巨大的锡矿找矿潜力。最近，Liu et al、祝新友等和Wang et al研究认为，维拉斯托矿床的银铅锌矿体和锡矿体是同一成矿系统的产物，并指出浅部为热液型脉状铅锌银矿，深部则可能存在岩浆热液型锡（钨）多金属矿。因此，维拉斯托矿床深部锡多金属矿的发现不仅改变了人们以往对大兴安岭南段锡、银、铅锌成矿规律的认识，而且将影响该区域今后的锡矿找矿勘查思路，意义重大。

2. 四川拉拉铜矿

四川拉拉铜矿田位于扬子准地台康滇地轴中段，属东西走向的金沙江断裂褶皱带与川滇攀西大裂谷南北向构造带的交接复合部位。拉拉铜矿是我国西南重要的大型铜矿，也是四川省最大铜矿生产基地。矿田内金属矿产以铜为主，次为铁、镍，伴生金属有银、金、钴、钼、铂、钯等。2012年以前，经勘查的主要矿床（点）42个，其中大型矿床1个、中型矿床4个，其余为小型矿床和矿点，已累计探明铜资源量130万吨，铁资源量2600万吨，镍资源量5万吨。2012-2014年，老矿山接替资源勘查项目在拉拉落凼铜矿南部红泥坡矿区勘查取得了重大突破，新增333+334铜资源量63.22万吨，实现重大找矿突破。

拉拉地区铜矿矿体赋存于古元古界河口群火山沉积岩系中，具有明显的层控特征。但与以往拉拉地区铜矿均赋存在河口群中部火山沉积旋回落凼组地层中不同的是，本次发现的红泥坡铜矿体主要赋存在河口群上部火山沉积旋回长冲岩组上段，其主矿体东西最大延伸1950米，南北延伸2100米；矿体平均厚度为10.37米；Cu品位平均为1.36%。红泥坡铜矿矿石矿物有磁铁矿、黄铜矿、黄铁矿、斑铜矿、辉钼矿、磁黄铁矿等；脉石矿物有钠长石、白云母、黑云母、钾长石、石英、铁白云石等。与典型的拉拉铜矿条带状矿石构造不同，红泥坡铜矿矿石构造以脉状、网脉状、角砾状为主。

拉拉铜矿田是我国重要的铜矿集区，一直以来是研究的热点地区，众多学者对其进行了详细的研究并取得了丰富的成果。然而，由于矿床的复杂性，拉拉矿田的成因仍存在争议，目前主要的观点分为两类：一是一些学者根据其赋存于河口群落凼组海相火山岩中，矿体呈层状、似层状产出，矿石主要为块状、条带状构造，结合部分硫铅同位素数据显示成矿物质来自于河口群地层，据此认为矿床属于火山成因块状硫化物矿床（VMS），并认为成矿年龄与成岩年龄一致为16～17亿年，但拉拉铜矿目前较为可靠的辉钼矿Re-Os定年数据显示其可能的成矿年龄为10.86亿年，且该矿矿化元素仅有Cu和Fe，无Pb或Zn，显示其可能的VMS矿床成因证据仍显不足。

②也有学者认为其特征很类似IOCG矿床，其发育很强的钠长石化，出现大量的贫Ti磁铁矿，并富集Co、Au、P、F、REE等，而贫Pb、Zn；C-H-O-S同位素数据显示其岩浆热液来源，有部分成矿物质来自于地层；且根据辉钼矿的Re-Os同位素年龄得到其成矿年龄为10.86亿年，并与当时板内拉张环境密切相关。但与典型的IOCG相对比，落凼铜矿缺乏大范围分布的角砾状矿石，并且受断裂控制不明显，所以其可能的IOCG矿床成因也存在争议。红泥坡铜矿的发现及其揭示的矿床地质信息，为研究者研究拉拉铜矿田的成因机制提供了绝佳的契机，通过对其系统的解剖，有望最终确定拉拉铜矿田的矿床成因并建立找矿预测地质模型，进一步指明区域找矿方向。总之，红泥坡铜矿重大的找矿突破打开了拉拉地区广阔的找矿空间，新赋矿层位的发现明确了本地区下一步的找矿方向，并增强了矿山企业“拉拉外围找拉拉”的找矿信心，为川滇有色金属成矿带开展“拉拉式”铜矿找矿提供了可借鉴的成功经验。

3. 江西相山铀矿

江西相山铀矿田位于钦杭成矿带北东段，赣杭陆相火山岩铀成矿带西南端，德兴-遂川大断裂的南东侧，是我国重要的热液型铀矿田。以往对于相山铀矿田的研究工作主要侧重于铀的成矿作用机制和找矿预测等方面，且主要集中在矿田浅部。随着矿山深部找矿的实施，多个深孔均见到较好的铅锌铜等多金属矿化，且品位较高，已达综合利用价值。例如，钻孔CUSD3在孔深1095.4~1549.2米、1574.2~1577.5米分别见到铅锌铜矿化脉和铜多金属矿化脉。又如，钻孔ZK26-101在深部见约330米的铅、锌、银等多金属矿化，初步估算铅锌银矿资源量分别为Pb1.7万吨、Zn 0.88万吨、Ag76.2吨。地质特征显示，深部的铅锌银矿化主要赋存于流纹英安岩、碎斑熔岩及基底变质岩中，受断裂构造、岩层界面、火山塌陷构造复合控制；铜矿化主要分布于基底变质岩中。最新研究表明，相山铀矿田浅部的铀成矿作用与深部的铅锌铜成矿作用具有相同的成矿物质来源，属同一岩浆-热液成矿系统的产物。深部铅锌铜多金属矿化的发现构成了相山矿田上铀-中铅锌金-下铜的成矿空间模式，不仅加深了对矿田成矿规律的认识，拓展了深部找矿空间，而且开拓了区域找矿思路。

4. 南京市梅山铁矿

根据宁芜地区玢岩铁矿的成矿模式，梅山式铁矿为次火山岩体与火山岩接触带上下高温气液交代充填矿床和矿浆充填矿床，具透辉石-石榴子石-磷灰石-磁铁矿组合。近年来实施的梅山铁矿接替资源找矿项目中，在梅山铁矿层顶部硅化、黄铁矿化、碳酸盐化凝灰岩、沉积凝灰岩、凝灰角砾岩中发现了金矿体，同时这一层位也是梅山铁矿层顶部普遍存在的一个“近矿指示蚀变带”。据江苏省地矿局第一地质大队项目成果资料，ZK4101钻孔在392.2～407.8米处见到金处矿化15.6米，平均品位0.97克/吨；ZK4121钻孔在406.28～414.77米处见金矿层8.49米，平均品位1.88克/吨。初步估算新增333金矿资源量（金属量）2.85吨。

宁芜地区铜金矿主要有铜井地区的铜井式铜矿，矿床类型均属火山-次火山热液型矿床。此类矿床的矿化与偏碱性、碱性火山活动有关，受区域构造裂隙或火山构造裂隙的控制。矿体一般呈脉状沿构造破碎带或破火山口构造成群出现，陡倾斜产出，与围岩界线清楚。矿石组合多为铜、金、多金属矿物与石英（重晶石、碳酸盐矿物）组合，并以中-低温热液充填型为主。典型矿床有铜井金铜矿（中型）、观山铜铅矿（小型）、谷里铜矿（小型）、金驹山金矿（小型）等。梅山铁矿矿体顶部层位发现的含金蚀变带与铁矿的密切关系可能反映铁矿和金矿是同一岩浆热液不同演化阶段形成的产物，当然目前也不能排除后期热液叠加成矿的可能。尽管目前对铁矿化和金矿化成生关系的认识还有待于进一步深化，但这一发现对宁芜地区玢岩铁矿的找矿不乏指导意义。

5. 吉林夹皮沟金矿

夹皮沟金矿位于华北克拉通北缘，是一座资源几近枯竭的危机矿山。该矿床虽有近200年的开采历史，但其成因却一直存在争议。前人曾先后提出过绿岩型、层控型、韧性剪切带型、变质热液型、岩浆热液型等多种成因模式，尤其以韧性剪切带型成因模式最为主要。近年来，随着老矿山接替资源找矿项目的开展，在夹皮沟矿床深部取得了重要突破。经钻孔验证，在夹皮沟矿床深部742～754.9米处见10米厚的细脉浸染型金矿体，平均金品位为2.2克/吨。钻孔编录结果显示，细脉浸染型金矿体主要产于隐伏的石英闪长岩体顶部接触带，矿石矿物主要为黄铁矿、黄铜矿、方铅矿等，脉石矿物主要为石英、方解石等。此外，在八家子南西的头道溜河地区，还发现了与闪长玢岩体有关的爆破角砾岩型金矿床。多种金矿化类型的新发现，不仅加深了对夹皮沟金矿矿床系统的认识，而且为区域深部找矿提供了新的思路和方向。

对下一步深部找矿工作的建议 

1.持续推进矿集区找矿预测工作

老矿山是开展深部矿产资源勘查的理想场所，而矿集区是大中型矿山密集分布的地区。矿集区矿产地质调查程度和找矿预测直接关系到老矿山及其外围的深部找矿工作。然而，由于种种原因，许多矿集区尚缺乏系统的矿产地质调查及找矿预测，直接阻碍了深部找矿突破。矿集区找矿预测工作主要包括以下两方面内容——

一是在1∶5万矿产地质调查、地球物理测量、地球化学测量、矿产及异常检查、典型矿床研究、资源潜力评价等工作的基础上，确定重点工作区开展找矿预测。在重点工作区开展大比例尺专项地质填图（含修测）、专项物探、专项化探、专项样品采集及深部钻探探查等工作，系统研究工作区内主要矿床类型的成矿地质体、成矿构造和成矿结构面、成矿作用特征标志，构建找矿预测综合信息模型，预测矿体赋存位置，评价资源潜力，提交预测资源量，引导后续勘查。

二是围绕深部找矿预测需求，开展矿集区及老矿山“三位一体”（成矿地质体、成矿构造和成矿结构面、成矿作用特征标志）找矿预测理论与方法创新与示范，主要包括：1）成矿构造与成矿结构面研究方法与模式建立；2）开展矿床深部定量预测理论与方法研究；3）基于2DGIS/3D建模平台构建矿集区“成矿地质体-成矿构造与成矿结构面-成矿作用特征标志”找矿预测模型。

2.进一步加强勘查区找矿预测理论与方法研究

勘查区找矿预测理论方法体系的创建，首次提出成矿地质体、成矿结构面、成矿作用特征标志的概念，在实际使用中取得了良好效果，切实指导了深部找矿。但是，勘查区找矿预测理论目前仍不完善，需要在实际的应用过程中进行修正。

一是要结合矿山深部和外围找矿的新发现，重新认识各类矿床的成矿规律、成矿系统的发育深度和不同类型矿床的分带、叠合规律。以深部找矿为目标，通过建立典型矿集区脉、层、块、体矿化样式组成的上下、左右多元空间矿床矿化系统结构模型，突出反映找矿信息，进而指导矿山深部和外围找矿工作。

二是矿床模式的研究要从矿床的角度走向矿集区，从单个的矿床成矿模式发展为典型矿床成矿模式之间的组合模式，从构造体系控矿发展为构造成矿系列的阶段，这对于认识成矿系列控矿的规律，深入总结和认识矿床和矿集区成矿规律，提高对深部矿床成矿理论的认识和指导矿产资源勘查实现重大的突破具有重要意义。

3.加强科技引领，推进新方法和新技术的研发应用

从近年来矿山深部找矿的经验来看，重大成果的取得无一不是根据新现象，打破旧思维，结合实际情况进行理论、方法和技术创新的结果。科技创新无疑将对未来深部矿产资源勘查工作起到重要引领作用。因此，在开展深部找矿勘查工作的同时，还应重点研发和推广适合深部找矿的物化探等技术方法。

一是继续贯彻“三深一土”国土资源科技创新战略，针对矿集区3000米以浅的地下空间，重点研发覆盖区探测技术和深部地质结构与成矿系统探测技术，开展覆盖区物质组成识别标志研究、矿集区深大剖面探测、深部成矿系统蚀变标志研究等工作，查明各类重要成矿要素在深部空间的分布特征，尤其重点探明与成矿有关地质体、成矿构造与成矿结构面、成矿作用特征标志等关键成矿要素的空间展布规律和形态特征，建立符合矿集区深部找矿的地球物理和地球化学等技术的指标体系，构建矿集区深部地质三维结构模型。

二是对成熟的新技术新方法进行推广应用。在科学分析矿集区地质条件的基础上，根据深部新发现，充分利用KGR抗干扰电法仪激电测深、大比例尺低飞航磁测量、井-地磁测反演技术、构造地球化学测量、1∶5万抗干扰电法扫面等先进技术方法，进行推广应用，总结出一套适合本地区寻找同类型矿床的物化探方法组合，进一步指导矿集区及外围深部勘探工程布置，开展深部矿产勘查示范。

总之，深部找矿新发现进一步促进了深部成矿规律的认识，开拓了找矿思路，已成为持续推进矿产资源勘查“向深部进军”的强大动力。同时，深部找矿突破的实现，还需要系统的矿产地质调查、完善的找矿预测理论和先进的探测技术等作为强有力的支撑。

1 前言

近年由于常规天然气资源量和产量的下降，特别是在北美洲，非常规天然气得到了高度的重视。一些估计表明，全球非常规天然气资源量（不含水合物）超过30000万亿立方英尺，大约有50%的资源来自页岩气。Julander能源公司的首席执行官Fred Julander认为页岩气（SG）是“自发现石油以来最重要的能源进展”。

水平钻井技术的进步、水力压裂、相对高的天然气价格（相比2009年之前）和近来在巴内特页岩（Barnett Shale）和美国其他几个页岩气藏的商业成功都使页岩气在美国成为了热门能源，而且页岩气的勘探开发已开始蔓延到加拿大和世界其他几个地区。

由于页岩气远景的复杂性和广泛性，针对页岩气的应用不能采用普遍用于常规气和煤层气的应用技术，而需专门设计开发工具和方法。多名学者包括Gray等人（2007）和Harding（2008）认为基于确定性解决方案的决议不适用于页岩气开发，因其没有考虑与复杂成藏有关的风险和不确定性，且经常导致过于乐观的结果。

到目前为止，尽管在北美和欧洲的勘查活动活跃以及近期商品价格下降，页岩气远景分析工作也只完成了极少的部分。商品价格的下降使最高质量远景区的开发至关重要，这些区域的开发不仅最符合公司的利益，并且赋予公司与国外的低成本常规气田（即卡塔尔和沙特阿拉伯相关的天然气）竞争的最佳潜力。Williams-Kovacs和Clarkson（2011）提供了与非常规的远景分析有关的现有工作的回顾，并提供了一种专为页岩气应用而设计的综合的六阶段远景分析及开发评价方法（PADEM）。本文中，作者还展示了一个专门开发用以筛查页岩气远景区并且选择最适合详细分析远景的工具。本文以Williams-Kovacs和Clarkson的工作为基础，致力于远景评价并选择进行更深入分析的远景区的试点位置。

当前工作的目标是：①开发一种协助页岩气勘探开发阶段的方法和配套的分析工具；②演示已开发技术在加拿大西部致密砂岩/页岩远景区的应用。这项工作的主要贡献是开发与示范一种针对页岩气远景区的严格分析方法。当考虑共存关系时，基于先导试验井输入变量的不确定性，该方法能生成其预测的分布。以前所有的工作一直专注于全域开发方案，然而无法利用勘探开发早期阶段可获取的少量数据快速形成这种全域开发方案。

2 工具开发

在这项工作中，开发了一种用于分析页岩气远景的工具。该工具选择使用（以Williams-Kovacs和Clarkson提出的方法（2011）为例的）预筛选的方法。本文将重点放在该工具的开发和应用，分析某一远景区的不同区域，以确定它们是否是适合的试点项目，并描述了图1所示的PADEM工作流程的勘探阶段。勘探阶段的目的是对从更多的详细资料中筛选的远景进行调查，以增加对油藏流动性和碳氢化合物生成能力的了解。在这项工作中，我们对个别类型油井采用概率范围经济学（probabilistic scoping economics）作为勘探标准，以确定该远景区是否适合实行试点项目。表1中完整提供了Williams-Kovacs和Clarkson（2011）详细讨论整体勘探开发方法的总结。

表1  勘探开发方法概况

在这项应用中解析模型比数值模拟更适用，其原因在于应用程序自设置和初始化的时间很短，整合的蒙特卡罗模拟法简单易行，并且在勘探早期阶段不容易获得形成精准的数值模拟所需的详细数据。尽管数值模拟技术已得到改进，但解析方法在工业和文献中依然被大量使用。下文给出了开发工具的关键部分的概要。

2.1 属性图

勘查方法最关键的组成部分可能是关键储层、地质力学、岩石物理和地球化学特性的精确属性图的开发。从地质模型、产量不稳定分析（RTA）、压力不稳定分析（PTA）、岩石物理调查等组合中可以推导出这些属性图。这些属性图用于远景的可视化、区块选区以及单一区块的分析。天然气原始地质储量图（OGIP）、K_m-h图、压裂脆性图等有助于选择代表性区块以及具备更大开发潜力的区块，甚至高度非均质性区块。区块作为一种评价不同区块远景生产特性的方法，基于地质和岩石物理的观察，比较简单易于操作。采用区块方法不需要针对每个勘探网区块开发一种标准井进行分析，然而通过应用蒙特卡罗法依然解释了其变化性和不确定性。Clarkson和McGovern（2005）采用区块方法评价了煤层气（CBM）远景。通过输入X-Y坐标值以及Petrel^TM软件的储层属性Z值可以在Excel中创建储层属性图。随后，数据透视表程序被用于对数据排序，并利用二维绘图应用软件创建属性图。由于早期的岩石物理模型通常利用有限的数据集开发，单一区块在蒙特卡罗模拟中选择不确定的输入数据和参数范围可以解释模型参数的不确定性。这种解释不确定性的方法将在本文所示实例中进行演示。

2.2 水力压裂模型

该项工作中，水力压裂裂缝的半长采用Valko（2001）提出的在常规和致密气中应用的简单双翼压裂模型来预测。该模型采用基质渗透率、剪切模量（杨氏模量与泊松比的函数）以及其他储层参数作为输入数据，且如果建模的输入参数不确定，则都必须重新计算每次蒙特卡罗迭代。采用简单的关联（A_cm＝4x_fh）可将裂缝半长转换为与压裂有关的面积。这个压裂模型可能无法代表部分更复杂的页岩气裂缝。为了更好的表示引入到大部分页岩气储层的复杂压裂网，Xu（2009，2010）等人建立了一个更具有代表性的水力压裂模型，该模型将被结合到本次工作中所演示的更新版本的方法中。该区的微地震观测表明，在本文预测的远景区横向双翼压裂的假设是合理的。

作为所应用的速率预测模型中的关键组成部分必须估算裂缝半长，这一问题将在下面部分开展讨论。水力压裂裂缝半长在随机分析中作为不确定的输入量，其分布主要根据该地区的微地震事件或者其他方法来确定。

 

 

图1  非常规天然气勘探阶段的勘探/开发方法工作流程

2.3 速率预测

Clarkson（2013）提供了关于页岩气井生产分析和速率预测综合全面的概述。在该工作中，我们将页岩气井理想化为一个矩形双孔介质系统，气体从基质岩块流入到裂缝且储层不随着裂缝延展（如图2的概念模型）。该模型忽略了包括体积压裂（SRV）在内的影响，其他作者认为大部分低渗页岩气井在合理的时间内不会发生体积压裂。此外，图2所示的概念模型假设了一个均质的完井——Amborse等（2011）和Nobakht等（2011a）讨论了非均质储层完井的预测。

在本次工作中，该模型的解决方案首先由EI-Banbi（1998）提出来。人们普遍认为在页岩气藏中占主导地位的瞬时流动状态是从基质到裂缝的线性流。同时，也可能出现一个与水力压裂线性流动相关的线性流动周期，但是通常认为这个阶段持续时间很短，或者被水力压裂清理以及表皮效应所掩盖，而很少可用于分析。本项工作中，我们假设瞬时线性流（从基质到裂缝）之后是边界控制流，该流态与受表皮效应（见等式7）影响的线性流体模型存在早期偏差。压裂段之间的不渗透边界结构导致了边界控制流产生。由Wattenbarger等（1998）首先将早期线性到边界控制流体的假设引入到致密气的应用中，并且该假设被广泛应用于文献和页岩气行业的解析模型。

 

 

图2  从线性流到边界流的解的概念模型

2.3.1 瞬时线性流的速率预测

EI-Banbi（1998）提出通过恒定速率和恒定流体压力来描述瞬时线性流的公式。本项工作中采用恒定流体压力的条件，这也是本文其他部分的重点——该边界条件最接近大部分产生达到最大水位降低值的页岩气井的流动条件。Samandarli等人（2011）采用不同的流体压力迭代方法，对页岩气生产进行分析建模，但是他们表明在大部分情况下采用恒定流体压力的假设就可以了。

与常用于表征简单横向双翼压裂的裂缝半长（X_f）相比，相关储层面积（A_cm）能更好的表示完井措施和增产措施效果以及生成复杂裂缝的能力。因此，在这一分析中，采用相关的储层（气藏）面积（A_cm）取代裂缝半长（X_f）。许多业内专家相信由于页岩气藏超低的基质渗透率，复杂压裂对于页岩气的商业生产至关重要。

无因次时间，t_D,Acm，相关储层面积（A_cm）依据公式1在恒定压力条件下定义。

                           （1）

无因次速率，q_D,Acm，由无因次时间定义：

                                       （2）

基于储层特性的无因次速率表达式，如果可获得关于K_m和A_cm的估算值，通过公式（3）可确定气体流速。采用不稳定产量分析或者其他的模拟技术可估算K_m和A_cm。K_m也可以通过实验室技术单独确定。

                             （3）

Ibrahim和Wattenbarger（2006）认为线性流的性能受水位下降程度的影响，同时提出水位下降量修正因子（f_cp）。此次工作中采用的修正因子（f_cp）由公式4给出。

                            （4）

此处，

 

Nobakht等人2011a和Nobakht等人（2011b）通过分析中采用校正时间（本次工作未采用）提出一种更严格的校正水位下降量的方法。

将水位下降量修正因子应用到公式3得出公式5：

                       （5）

除了水位下降量的修正，这些公式经过进一步修改可直接应用于页岩气井。与致密气井相比，大部分页岩气井在时间曲线的平方根中表现出的较大截距（在致密气井中曲线通常穿过原点），而在流量和时间双对数曲线上页岩气井则呈现出的一半斜率的偏差。多名作者最初认为是裂缝的有限导流能力造成了这种偏差，但是Bello（2009）和Bello和Wattenbarger（2009，2010）认为这种偏差可以通过采用表面效应来更好的解释。Bello（2009）、Bello和Wattenbarger（2009）在恒定流量和恒定流体压力条件下完成了大量的受表皮效应（skin effect）影响的线性流分析，且推导出了恒定流体压力条件下的解析解。在他们的分析中，将表皮效应作为一个常量。Bello（2009）和Bello和Wattenbarger（2009）证明恒定流量情况下表皮是附加量，而恒定流体压力情况下表皮的作用是非线性的。由Bello和Wattenbarger（2009）提出的解析式可以使用下面的近似代数方程：

                （6）

从方程（6）可以看出，当t_D(t)值大时，包含表皮的项就会变小。

Nobakht等人（2012）研究了巴内特、马塞勒斯和蒙特利的大量页岩气井（这些气井在相对恒定的流压下产量不断降低），同时得出结论：通常这些页岩气井更多表现出恒定流量的情况而不是恒定流压的情况。作者假设这种意想不到的表现可能是由于Bello（2009）以及Bello和Wattenbarger（2009）提出的表皮模型太过理想化，因此无法代表野外条件。通过假设恒定的表皮效应，模型不能说明由压裂清理、压力敏感地层、变化的压裂导流能力、变化的井底流压、压力相关的流体性质、变化的井筒流体梯度、液体加载等导致的表皮改变。作为这项工作的结果，作者提出了一个可应用于公式（2）的替代表皮修正项：

                   （7）

包括水位最低量和表皮的影响，公式（1）、（5）、（7）能够利用预测的气体流量，作为时间的函数，在线性流区域可对K_m和A_cm给出独立的估测。

2.3.2 边界控制流的流量预测

上面描述的方法适用于有效的储层边界相互接触，边界控制流形成之前。基于图2所示的几何图形，边界控制流紧随着瞬时线性流的末期出现。当外部SRV的影响较为显著时，这一观点较为保守。Clarkson和Beierle（2011）认为如果遇到了其他的瞬时流区，则应采用多重分区的方法，此外，如果多级压裂井需要进行非均质性储层的完井（heterogeneous completion），早期线性流之后不会立刻发生真实边界控制流，且需要更复杂“混合”预测技术。如同下面叙述的，我们选择采用更为保守预测程序，假设线性流之后紧随边界控制流。

利用公式8计算达到线性流的拟稳态时间（或者是瞬时线性流的结束时间）：

                       （8）

正如图2中看到Y_e是压裂到储层边界的距离，计算公式如下：

                       （9）

多名作者已经提出了页岩气井拟稳态线性流的预测方法。包括Fraim和Wattenbarger（1987），Palacio和Blasingame（1993），Doublet等（1994），Agarwal等（1999）和Mattar和Anderson（2005）认为可采用物质平衡类模拟程序预测边界控制流。Clarkson和Pedersen（2010）将这种方法应用于致密油研究，同时本文也将采用这种方法。公式（10）给出采用物质平衡方法预测边界控制流的生产速度：

              （10）

此处q_pssi-Linear是边界控制流初始的页岩气流体速度，（Pri）pss是边界控制流初始的平均储层压力，且（Pwfi）pss边界控制流体初始时井筒流体压力。通过物质平衡计算平均储层实际气体拟压力。对于含有大量吸附气的页岩气开采（application），一般使用Clarkson和McGovern（2005）提出的MBE方法。而在以游离气为主的情况下，则使用定容气藏的常规MBE方法。物质平衡计算需要地质储量和气体特性（比如天然气压缩因子），这两者都是由关键PVT输入量和状态公式（EOS）确定的。

（a）	收入总额	（b）	收入总额
扣减	使用费	扣减	使用费
扣减	运营成本	扣减	运营成本
得出	税前运营现金收入（OCIBT）	扣减	资金成本补助（CCA）
扣减	收入税	扣减	加拿大开发费用（CDE）
得出	税后运营现金收入（OCIAT）	扣减	加拿大勘查费用(CEE)
扣减	资本支出	扣减	加拿大油气物业费（COGPE）
得出	税后现金流（CFAT）	得出	生产应税所得
贴现	税后贴现现金流（DCFAT）	乘	生产税率
		得出	应付税款
		扣减	免税额度
		得出	应付净所得税

图3  现金流分析：（a）现金流；（b）收入税（加拿大税制）

结合El-Banbi（1998）改进的瞬时线性流的无因次公式和边界控制流的物质平衡模拟方法，可以开发一种综合的预测方法：

1)        获取A_cm（或者X_f）和K_m（来源于微地震和/或RTA模拟/已有生产数据或者其他估计）的独立估算值。

2)        使用公式（1）和（7）作为时间函数计算t_D,Acm，q_D,Acm。

3)        线性流部分的数据利用公式（5）作为时间函数计算q_g。

4)        指定排放区（来源FMB模拟/已有的生产数据或者其他估算）。

5)        使用公式（8）和（9）计算t_PSS-Linear和Y_e。

6)        确定 。

7)        采用公式（10）通过废弃量（边界控制流）从t_PSS-Linear预测产量。

上面描述的解析模型是假设模型（最小变化）区块内的体积平均值参数是恒量，并从认为是不确定的参数的概率分布中选择一个值。每一次蒙特卡罗迭代将选择不同的值，导致不同的流量预测和不同的主要经济指标值。在许多参数高异质性水平的情况下，存在明显的不确定性，这种不确定性反映在关键输出参数的显著变化。

2.4 经济模块

将经济模块与速率预测集成来计算与生产相关的现金流。因为通常行业采用名义美元计算实际（通常的）现金流和名义（现行的）现金流，虽然采用实际的盈利指数计算项目的最低预期资本回收率，且通过不同的通货膨胀率来比较项目。采用图3中的业务流程计算现金流和收入税（加拿大税收制度）。

该模块中的天然气价格的确定实行了价格操纵，而非价格预测。采用价格操纵表明了项目十分稳定（不论是单独而言还是相较于其他项目），并且不再需要预测极不稳定的天然气价格，该模块中也设置了以价格预测为基础引导经济的选项。

方法中建立了多个实际盈利能力的指标，包括净现值（NPV）、内部收益率（IRR）和投资收益率（ROI），用来比较项目和公司设定的最低预期资本回收率，同时可给项目进行排序。

2.5 蒙特卡罗模拟的一体化

本次工作将蒙特卡罗模拟整合到方法开发中。采用@RISK^TM（Palisade Corporation，2010）对关键PVT和储层属性（原始参数）进行概率分布和模拟操作。概率分布的输入变量根据不同项目的数据数量和质量而变化。Clarkson和McGovern（2005），Haskett和Brown（2005）和Harding（2008）认为对数正态分布最能代表PVT、储层和经济特性，因此本文使用了这种分布类型。这些概率分布拟合按P10（低）、P50（中）和P90（高）不同的值输入各个不确定变量。这些输入值可能来自勘探/远景数据、个人经验、模拟数据等。缩减所有输入变量的分布保证每个实现只选择合理的数值（缩减分布将选择少量接近无穷大的数值，从而影响输出变量）。

上面讨论了@RISK^TM输出变量定义的关键经济参数，以及气体速率和累积天然气产量。由于每个输出变量允许量化与项目相关的不确定性，可对其生成一个概率分布，以便做出与远景选取和开发有关的明智决策。

通过在x轴上找到相应的最低预期资本回收率时的位置，向上垂直移动至曲线处，然后再水平投影到y轴，这样可以从累积概率分布计算出超过设定最低预期资本回收率的概率。用1减去y轴上求出的值，得出超过最低预期资本回收率的概率。这个方法在本文中将作为范例进行演示。

在这一应用中（如在孔隙度和渗透率之间），采用了拉丁超立方体抽样，如果有必要的话，还可合并相关性（如孔隙度与渗透率）。典型的多相（气+水）页岩气/致密气应用的主要参数如表2.3所示。在某些情况下，参数的依赖关系可使用行业普遍接受的经验模型进行解释，而在其他情况下会使用来自现场数据或者估算得到的基于方向的相关性（如较高的正相关关系）。例如，与压力有关的渗透率（绝对的渗透率比值）使用Yilmaz等人（1991）的方法可与储层压力和岩石力学特性关联。相反，束缚水饱和度与孔隙度密切正相关。可能的参数关系如表2所示。

蒙特卡罗模拟运用了一个类似于Clarkson和McGovern（2005）使用的煤层气气藏远景分析的方法。

表2  基本参数、可能的相关性和参数关系

3 该方法应用于远景勘探

本文中开发的方法广泛应用于SG远景将其分成区块进行分析的目的，以确定是否适合作为一个试点项目。由于SG试点和开发项目成本高，且其详细分析需要大量数据，页岩气远景勘探至关重要。

对于远景勘探应用而言，其方法的选择以当前远景数据和模拟数据相结合为基础。理想情况下，对于关键PVT和储层参数情况良好的估计，作为空间坐标的函数可用于远景勘探。如果事实并非如此，可以对模拟气藏或者其他数据源进行估算以获取数据，同时分析该方法带来的不确定性。

假定整个远景区PVT和其他储层特性不变，输入数据可用于生成主要储层特性图。关键生产指标图如OGIP和基质渗透率乘以可以开发的净投入（千米/小时），可用于区块的选择。区块的选择基于区域类似的关键生产指标的值。对页岩气储层而言，压裂的指标，如压裂指数或脆性也可能用于区块选择，同时许多作者表明建立复杂裂缝网的能力对于页岩气商业开采至关重要。

选择区块后，开始进行蒙特卡罗模拟，按照P10、P50、P90的概率预测和可以开发累积产气的区块，且结合使用关键经济指标的分析来确定区块能否适合一个试点项目。其他因素比如公司的经验，企业和商业策略，可用的资源和基础设施等都将纳入评估，以便为公司以及股东们确定哪些区域可以作为最佳试点选项作出明智的决策。

远景勘探方法工作流程见图4所示。

4 采用两段页岩开发模型的样本示例

为了进一步说明该方法的应用，对加拿大西部的某处致密砂岩/页岩（假定没有吸附气体）远景区的两段进行了分析。在之前的研究中，Petrel^TM开发的远景地质模型采用可用的岩石物性、储层和生产数据。图5所示研究区域内4口井的三维孔隙度模型和孔隙度相关的钻/录/测井记录。在该区域，存在两处可获益的产气水平井段（井段3和井段4）。

 

输入数据 关键储层属性的填图属性 PVT，其他储层和水力压裂属性 生产数据 经济投入

区块选择 根据OGIP或者其他关键属性确定区块

蒙特卡罗模拟

模拟输出 P10、P50、P90的概率预测和累积产气量 水力压裂运行情况 经济参数 可行的商业区块标志 其他

 

 

图4  远景勘探方法的工作流程

 

 

图5  三维孔隙度模型和孔隙度相关的测井

模型开发期间这个开发区拥有11口垂直井，2口倾斜井和4口水平井。最初钻完成垂直井，紧随其后的是开始于2008年的水平井。Clarkson和Beierle（2011）在该区选择一系列井进行不稳定产量试井（RTA）。模型开发中使用的水平井的总结显示在下面表3中，同时在图6中（在下面描述）该区域的天然气原始地质储量（OGIP）图上显示了井的近似轨迹。

表3  研究区水平井概况

井名	井向	进入层位	完井方式
1号井	水平	井段4	尾管注水泥
2号井	水平	井段3	自膨胀封隔器
3号井	水平	井段4	自膨胀封隔器
4号井	水平	井段4	自膨胀封隔器

所做的分析主要集中在大部分是水平井的井段4。为了简化分析，采用孔隙度下限为4%，通过Excel加权平均井段4层位，将Petrel^TM多层模型转换成一个单层模型。这一平均化过程是为了完成对基质的孔隙度、初始含水饱和度和渗透率的处理。利用孔隙度下限值还可以计算总有效收益和毛净收益（有效收益假设包括所有孔隙度下限值以上的层）。图7a和图8a显示了OGIP和K_m-h属性图。

模型采用的网格大小如表4所示。在整个开发过程中假设为常量的PVT、储层和生产参数如表5所示。

表4  网格属性

对于这种情况，人们认为井筒流动压力（pwf）为常量1750磅/平方英寸，接近开发区水平井最初的井筒流动压力。随着时间的推移井筒流动压力降低，后期模型中压力驱动力低于开发井，模拟气率并不乐观。这种情况下，在可获取日常生产和流动压力期间内，平均两个收益井的流动压力大约是1550磅/平方英寸，因此到开发后期之前，这种假设的影响并不很明显。在实际勘探中，该地区还没有投入生产，由于我们不需要将可用的生产数据与模型匹配，而是采用实际的流动压力估计值尝试得到一个准确的潜在生产能力估计值，所以这种假设的影响不是一个值得关注的问题。

表5  PVT常数、储层和生产投入参数

参数	值
PVT参数
气体比重	0.648
N2/%	0.46
CO2/%	0.2
H2S/%	0.0
温度/℉	166.5
Cw/磅/平方英寸-1	2.9×10-6
Cr/磅/平方英寸-1	5.6×10-6
VL/标准立方英尺/吨	N/A
PL/磅/平方英寸	N/A
储层参数
Pi/磅/平方英寸	3500
排放面积/Ac	80
生产参数
Pwf/磅/平方英寸	1750
rw/英尺	0.3

3个区块中假设关键属性的变化情况如表6所示。各属性的数值是每个区块的各个网格值的算术平均数。由于基质渗透率是蒙特卡罗输入量，且利用基质渗透率值可计算总压裂半径（虽然也可使用压裂分析模型在每次迭代时作为基质渗透率函数计算总压裂半径），故给出了一个基质渗透率值以显示区块之间总值的变化情况。

表6  储层变量和水力压裂输入参数

参数	区块1	区块2	区块3
储层参数
有效厚度/英尺	102	74	58
孔隙度/%	7.1	6.5	6.0
Sw/%	18	15	16
Km/毫达西	0.0084	0.0079	0.0077
水力压裂参数
剪切模量/磅/平方英寸	2×106	2×106	2×106
总压裂半径/英尺	1432	1477	1489

 

 

图6  研究区地质储量图呈现近似水平井轨迹

4.1 区块选择

利用从Petrel^TM多层模型开发的单层模型，其单层等量地质储量如图7a所示。根据类似颜色为代表的区域具有类似地质特征和岩石物理性质，通过视觉观察可选择区块。虽然已知气藏具有高度的横向非均质性，可以看到关键的地质和岩石物理性质明显凸出部分。该图形显示了更复杂的异质性模式的情况，需要更多的区块并且可能有必要用区块代表具有相似属性的不连续块段。图7b显示基于天然气原始地质储量选择的区块远景区。在计算天然气原始地质储量时，虽然该远景区吸附气体量很容易被包含其中，但还是假设其可以忽略不计。

 

 

图7  地质储量图：（a）地质储量；（b）选区

从图7b可以看出选取的三个区块中，区块1具有最高的天然气原始地质储量（红色和橙色），区块2具有的地质储量（光和暗绿色）次之，区块3具有的地质储量（紫色和蓝色）最低。从这幅图中可以推断出区块1将有最理想的属性，因此可能具有最高的产量，而区块3产气物性最不理想，因此可能具有最不理想产气量。如同气藏地质储量图（图7）一样，如果绘制K_m-h图我们也可以分辨出三个相似的区块。此次应用区块选区采用的天然气原始地质储量图和K_m-h图作为代表资源的程度/密度和储层特性的两个要素，这是工业上常用的评估致密砂岩和页岩远景好坏的关键因素。区块选区的属性根据不同项目而变化，取决于驱动特定资源类型远景的关键要素。

对于这种情况，假设简单的水平双翼压裂（如所使用的压裂模型所假定的）就足够了，因为微地震数据对同一区域的补充水平压裂井的解译说明复杂程度较低，如果不是水平情况，则进行压裂（图9）。采用水平和垂直观察井用以观察，同时采用双阵列处理会产生一个好的数据集。一般情况下，各个阶段仅出现一个水力压裂裂缝。水力压裂裂缝通常选择北东-南西方向，与加拿大西部沉积盆地（WCSB）部分最大水平应力方向一致。

 

 

图8  K_m-h图：（a）K_m-h；（b）选区

通过比较图6与图7b和8b可以看出在开发区所有水平井部分或全部在区块1范围内。因为这个原因，剩余的分析还将在区块1中开展。对区块1区域的水平井的预测情况而言稍微乐观，因为这些水平井水平延伸超出区块1区域进入地质储量和K_m-h更低的区域（该区水平井采用恒定的流体压力与（P_wf）i相比将获得相反的影响）。

 

 

图9  根据微地震数据解译的研究区内水平井水力压裂裂缝几何图形

4.2 经济分析

分析假设只有天然气价格是变量，而所有其他经济参数都保持常量。表7列出了其他主要经济参数的值（基于Magyar和Jordan的估算（2009））和表8介绍了主要的专利权使用费、税和贴现参数。

在本文的分析中，净现值（NPV）作为重要的收益经济指标且最低资本回报率为0。

分析远景的工作流程图如图4。

表7  资本和运营成本参数

表8  使用费、税收和折现率

4.3 蒙特卡罗模拟

在区块选择之后，本文进行了蒙特卡罗模拟研究。蒙特卡罗模拟中，基质渗透率（k_m）和页岩气价格不断变化，而所有其他的PVT、储层参数和经济参数保持不变。为了更好地进行说明，我们选择了将“不确定”的输入变量的数量显著限制在基本控制远景的油藏性能（储层渗透率）和经济情况（天然气价格）。基于P10、P50和P90值按照对数正态分布模拟参数。在大多数的勘探情况下，许多参数都是不确定的，可以通过这些参数的概率分布（见表2）来定义。对于需要使用概率分布进行定义的一些关键参数，可通过评估给定区块内重大变化的属性图来直接确定，或用更严格的统计技术，如采用区块内部数值计算变异系数（Cv）。由于基质渗透率是基质流动的主要控制要素，以及未来商品价格造成的天然气价格的高度不确定，针对这种情况，我们选择基于视觉观察的基质渗透率。

基质渗透率按照P10、P50和P90的值计算如下。通常情况下，可以通过岩石物理模型中的参数值拟合分布来生成概率分布，但是因为我们处理的是远景的早期评估，因此我们采用了替代的方法，即最大限度提高模型获取的不确定性来解释其他早期参数估算无法获取的变化性。如果需要，对其他不确定参数也可以使用相似的方法。

P10——区块1中比第十百分位值的基质渗透率低20%

P50——区块1中的基质渗透率值居中间数

P90——区块1中比基质渗透率的九十百分位值高20%

表9中定义了2个输入变量的分布。将模型内部不确定参数合并关联（见表2）也很重要。虽然孔隙度和渗透率之间的相关性被加入到原始岩石物理模型（幂律相关），并且压裂半径与剪切系数（正相关）、基质渗透率（负相关），净收益（负相关）和压裂模型井眼半径（负相关）相关，但是出于演示的目的，本文对这一方法进行了简化，使蒙特卡罗模拟中的主要变量之间没有相关性。由于压裂半径取决于基质渗透率，压裂模型必须在每次迭代时重新计算。气体流量，累积产气量和净现值被定义为@RISK^TM输出变量。

本文进行了5000次蒙特卡罗迭代，以确保蒙特卡罗输入变量充分覆盖样本空间。要求覆盖足够的样品空间，是为了确保每个模拟输入相同参数运行时，能得出同样的结论。出于演示的目的，用上述方法获得的迭代数并不是最优化。但是，通过将无限大（非常大）的样本输出分布与减少样本数量的输出分布比较，同时寻找要求充分重复“已知”输出分布的最小值，可以获得优化的迭代数。当进行多个模拟时，优化处理可用于减少处理时间和容量。

4.4 结果

图10显示了区块1中单口气井的确定产气量和累积产气量预测。这个“确定性”的基质渗透率的值来自于表9所示输入分布的斯旺森平均值（SM），假设这个值代表区块收益的平均水平（静态平均K_m=0.0095毫达西）。虽然Bickel等人（2011）指出了斯旺森平均值（SM）的缺点，但它仍然被广泛地用于工业，因此在这种情况下还将使用。此外，斯旺森平均值在输入分布的平均值的5%范围内（使用@RISK^TM计算），因此认为在这个例子中的平均值是准确的。另外，可以使用另一个估计的平均值（即分布平均值、区块值的算术平均值等）。图10a显示的产气速率与时间半对数图以及累积气体的产生与时间的笛卡尔曲线，而图10b显示了产气速率和时间的对数分布图。

 

 

图10  开发模型情况下的确定速率预测：(a)产气速率和时间、累积产气量和时间的半对数；（b）产气速率和时间的对数关系

图11显示了产气速率与时间的半对数图，图11b显示一个产气速率与时间的对数图和图11c显示预测（约14年）最初5000天累积产气与时间的笛卡尔曲线。

通过比较图10和图11，可以再次看到确定性预测与P50概率预测相比，具有更大的IP，持续的生产速度和累积产气量，表明确定性预测是比中位数情况稍微乐观，并且明显远超过P10的情况。这些结果再次支持使用概率分析取代非常规应用的确定性分析。

 

  

图11  开发模型情况中概率速率预测：（a）产气速度和时间的半对数关系；（b）产气速率和时间的对数关系；（c）累积产气量和时间

随后，P10、P50和P90产量预测与区块1内水平井可获取的生产数据进行对比，以测试开发方法的稳健性和准确性。在这个比较中，由于完井的复杂性，只有井3和井4可用，而井1表现不佳，且井2在此次分析区块外部。井3的产量被缩减了30天，以便使该井产量自然下降的初始时间与概率预测的一致（指修正井3）。生产的前430天的对比曲线如图12所示。

如图12所示，两口井的生产数据（修正井3和井4）普遍落在P10和P90之间（使用@RISK^TM生成的预测）。除了生产的前20天和第300天左右时的大约20天两个时间段（模型没有指出的操作问题导致的结果）外，约80%的数据点如预期处在P10和P90预测之间。初步预测产量可能更高，因为它不考虑压裂清理干扰、启动效应等，该模型增加了表皮效应来提高与IP的匹配程度。但是，在真正的勘探情况下表皮效应的大小无从得知，这是因为无法获取产气远景区域的数据且需要将其作为不确定的输入量以最大限度地提高模型的准确性。

 

 

图12  3号井和4号井生产数据和概率速率预测的对比：（a）产气速率和时间的半对数关系；（b）产气速率和时间的对数关系；（c）累积产气量和时间

虽然这不是一个令人满意的统计样本，只有一个关键属性（K_m）被认为是不确定的，但结果令人鼓舞。图13显示了净现值的增加的累积概率分布，直方图和回归系数托那多图。图13a再次显示超过最低预期资本回收率概率计算的累积概率分布图。

从图13a可以看出这个模拟平均净现值为53万美元，可能超过最低预期资本回收率的50%。然后，可将平均净现值和超过最低预期资本回收率的概率与相同远景的其他区块，以及与其他潜在远景的区块进行比较，从而确定哪些远景区域可提供最好的经济成功机会。这一分析显示了积极的NPV平均值和超过最低预期资本回收率的适度概率。基于这样的分析，可以得出结论：区块1的样品远景对于试点项目是极好的备选。这一分析支持了该地区的开发，但是这一测试中所采用的天然气价格网格假设对其结果影响极大。图13C中托那多图表明天然气价格对净现值带来的影响最大，基质渗透率给净现值带来的影响其次（区块1中最小的基质渗透率变化的结果）。这表明假设较高的气体价格（比如该区水平井钻探时期的气体价格）将提高远景的可取性。从图13b直方图可以看出模拟中大部分的净现值在300万美元和350万美元之间，众数等于-1.5万美元，相当于平均数53万美元左右。

 

 

图13  开发模型应用NPV法得出的经济结果：（a）累积概率分布；（b）柱状图；（c）回归系数的龙卷风图

此分析程序可在在开发区的其他2个区块内完成，以协助选择最适合公司的试点项目的位置。2号和3号区块的填图属性的直观观察（图7b和8b）表明，这些地区情况没有区块1理想，因此在本次分析所使用的气体价格假设中可能不适合作为试点项目。

5 结论

在本文中，开发了一种方法理论和基于excel的方法以协助页岩气和致密砂岩气藏的勘探。这个方法包含了来自不同来源的映射属性、一个用于估算水力压裂半径的简单的压裂模型、目前应用于页岩气井开采的速率预测技术、计算关键盈利能力指标的经济模块以及解释非常规资源中内在的风险和不确定性的蒙特卡罗模拟。本文所描述的方法和工具可被工业界用于评估远景区域内的各个区块和选择适合试点项目的地区。该方法较为严谨，以岩石物理、地质和现在产业应用的分析储层模型为基础，且通过重建现有实例的油藏动态来证明其准确性。由于不需要建立复杂的数值模型和详细的开发方案（所需数据是在开发早期通常无法获取），这种方法既简单又高效。

感谢代金友副教授对本文提出的宝贵意见。本文受中国地质调查“地学情报综合研究与产品研发”（121201015000150002）项目支持。

资料来源：Williams-Kovacs J. D., Clarkson C. R. A new tool for prospect evaluation in shale gas reservoirs. Journal of Natural Gas Science and Engineering，2014，18（5）：90-103.