页岩气地球物理测井评价综述

书书书 第２７卷 第４期 ２０１２年８月（页码：１６２４－１６３２） 地　球　物　理　学　进　展 ＰＲＯＧＲＥＳＳ　ＩＮ　ＧＥＯＰＨＹＳＩＣＳ Ｖｏｌ．２７，Ｎｏ．４ Ａｕｇ．，２０１２ 郝建飞，周灿灿，李　霞，等．页岩气地球物理测井评价综述．地球物理学进展，２０１２，２７（４）：１６２４－１６３２，ｄｏｉ：１０．６０３８／ｊ．ｉｓｓｎ． １００４－２９０３．２０１２．０４．０４０． ＨＡＯ　Ｊｉａｎ－ｆｅｉ，ＺＨＯＵ　Ｃａｎ－ｃａｎ，ＬＩ　Ｘｉａ，ｅｔ　ａｌ．Ｓｕｍｍａｒｙ　ｏｆ　ｓｈａｌｅ　ｇａｓ　ｅｖａｌｕａｔｉｏｎ　ａｐｐｌｙｉｎｇ　ｇｅｏｐｈｙｓｉｃａｌ　ｌｏｇｇｉｎｇ．Ｐｒｏｇｒｅｓｓ　ｉｎ Ｇｅｏｐｈｙｓ．（ｉｎ　Ｃｈｉｎｅｓｅ），２０１２，２７（４）：１６２４－１６３２，ｄｏｉ：１０．６０３８／ｊ．ｉｓｓｎ．１００４－２９０３．２０１２．０４．０４０． 页岩气地球物理测井评价综述 郝建飞，　周灿灿，　李　霞，　程相志，　李潮流，　宋连腾 （中国石油勘探开发研究院，北京１０００８３） 摘　要　页岩气是指生成、储集和封盖均发生于页岩体系中，或游离于基质孔隙和天然裂缝中，或吸附于有机质和 粘土矿物表面，或溶解于沥青和水中，在一定地质条件下聚集成藏并具有商业价值的生物成因和／或热解成因的天 然气．在页岩气勘探开发中，地球物理测井是识别、评价页岩气储层并为后期完井提供指导参数的重要手段．页岩 气属于极低孔极低渗的范畴，且具有很强的非均质性和各向异性，常规油气藏测井解释评价方法已不再适用，必须 建立新的页岩气测井解释评价体系才能够对页岩气藏做出准确评价．评价页岩气藏的潜力涉及对多种因素正反面 影响的权衡，包括页岩矿物组分和结构、粘土含量及类型、干酪根类型及成熟度、流体饱和度、吸附气和游离气存储 机制、埋藏深度、温度和孔隙压力等．其中，孔隙度、总有机碳含量和含气量等对于确定页岩储层是否具有进一步开 发价值非常重要．本文针对国外尤其是美国近期页岩气勘探开发的现状进行了广泛的文献调研，综述当前国外页 岩气地球物理测井技术的发展现状，针对勘探开发的不同阶段介绍常用的含气页岩的测井系列，然后 总结 初级经济法重点总结下载党员个人总结TXt高中句型全总结.doc高中句型全总结.doc理论力学知识点总结pdf 页岩气 测井响应特征，并详细论述了页岩气储层评价方法及储层评价的重要参数，包括有机碳含量、岩石矿物组分及含 量、孔隙度、含气量及岩石力学参数，最后提出我国页岩气地球物理测井研究存在的问题和发展趋势． 关键词　页岩气，测井识别，测井评价，测井技术 ｄｏｉ：１０．６０３８／ｊ．ｉｓｓｎ．１００４－２９０３．２０１２．０４．０４０　　　　中图分类号　Ｐ６３１　　　　文献标识码　Ａ 收稿日期　２０１１－０８－１１；　修回日期　２０１１－１１－１０．　　　　投稿网址　ｈｔｔｐ／／ｗｗｗ．ｐｒｏｇｅｏｐｈｙｓ．ｃｎ 基金项目　国家油气重大专项“复杂储层油气测井解释理论方法与处理技术（２０１１ＺＸ０５０２０－００８）”资助． 作者简介　郝建飞，男，１９８７年１月生，山东潍坊人，现为中国石油勘探开发研究院在读硕士研究生，研究方向为页岩气测井解释评价． （Ｅ－ｍａｉｌ：ｈａｏｊｉａｎｆｅｉ２００５＠ｙａｈｏｏ．ｃｏｍ．ｃｎ） Ｓｕｍｍａｒｙ　ｏｆ　ｓｈａｌｅ　ｇａｓ　ｅｖａｌｕａｔｉｏｎ　ａｐｐｌｙｉｎｇ　ｇｅｏｐｈｙｓｉｃａｌ　ｌｏｇｇｉｎｇ ＨＡＯ　Ｊｉａｎ－ｆｅｉ，　ＺＨＯＵ　Ｃａｎ－ｃａｎ，　ＬＩ　Ｘｉａ，　ＣＨＥＮＧ　Ｘｉａｎｇ－ｚｈｉ，　ＬＩ　Ｃｈａｏ－ｌｉｕ，　ＳＯＮＧ　Ｌｉａｎ－ｔｅｎｇ （Ｒｅｓｅａｒｃｈ　Ｉｎｓｔｉｔｕｔｅ　ｏｆ　Ｐｅｔｒｏｌｅｕｍ　Ｅｘｐｌｏｒａｔｉｏｎ　＆Ｄｅｖｅｌｏｐｍｅｎｔ，Ｂｅｉｊｉｎｇ１０００８３，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ　Ｉｎ　ｇａｓ　ｓｈａｌｅ，ｔｈｅ　ｇａｓ　ｉｓ　ｇｅｎｅｒａｔｅｄ　ｉｎ　ｐｌａｃｅ；ｔｈｅ　ｓｈａｌｅ　ａｃｔｓ　ａｓ　ｔｈｅ　ｓｏｕｒｃｅ　ｒｏｃｋ，ｔｈｅ　ｒｅｓｅｒｖｏｉｒ　ａｎｄ　ｔｈｅ　ｃａｐｐｉｎｇ ｆｏｒｍａｔｉｏｎ．Ｐａｒｔ　ｏｆ　ｔｈｅ　ｇａｓ　ｉｓ　ｓｔｏｒｅｄ　ｉｎ　ｍａｔｒｉｘ　ｐｏｒｏｓｉｔｙ　ａｎｄ　ｆｒａｃｔｕｒｅ　ａｓ　ｆｒｅｅ　ｇａｓ，ｔｈｅ　ｏｔｈｅｒ　ｐａｒｔ　ｉｓ　ａｄｓｏｒｂｅｄ　ｉｎ　ｔｈｅ ｓｕｒｆａｃｅ　ｏｆ　ｏｒｇａｎｉｃ　ａｎｄ　ｃｌａｙ　ａｓ　ａｄｓｏｒｂｅｄ　ｇａｓ，ａｎｏｔｈｅｒ　ｐａｒｔ　ｉｓ　ｄｉｓｓｏｌｖｅｄ　ｉｎ　ｂｉｔｕｍｅｎ　ａｎｄ　ｗａｔｅｒ　ａｓ　ｄｉｓｓｏｌｖｅｄ　ｇａｓ．Ｓｈａｌｅ ｇａｓ　ｓｙｓｔｅｍｓ　ａｒｅ　ｏｆ　ｔｗｏ　ｄｉｓｔｉｎｃｔ　ｔｙｐｅｓ：ｂｉｏｇｅｎｉｃ　ａｎｄ／ｏｒ　ｔｈｅｒｍｏｇｅｎｉｃ．Ｔｈｅ　ｇａｓ　ａｃｃｕｍｕｌａｔｅｄ　ｕｎｄｅｒ　ｓｏｍｅ　ｇｅｏｌｏｇｉｃ ｃｏｎｄｉｔｉｏｎ　ａｓ　ａ　ｃｏｍｍｅｒｃｉａｌ　ｕｎｃｏｎｖｅｎｔｉｏｎａｌ　ｒｅｓｅｒｖｏｉｒ．Ｉｎ　ｔｈｅ　ｅｘｐｌｏｒａｔｉｏｎ　ａｎｄ　ｅｘｐｌｏｉｔａｔｉｏｎ　ｏｆ　ｓｈａｌｅ　ｇａｓ，ｇｅｏｐｈｙｓｉｃａｌ ｌｏｇｇｉｎｇ　ｉｓ　ａｎ　ｉｍｐｏｒｔａｎｔ　ｍｅａｎｓ　ｆｏｒ　ｉｔｓ　ｉｄｅｎｔｉｆｉｃａｔｉｏｎ，ｅｖａｌｕａｔｉｏｎ　ａｎｄ　ｐｒｏｖｉｄｅｓ　ｓｏｍｅ　ｅｆｆｅｃｔｉｖｅ　ｐａｒａｍｅｔｅｒｓ　ｔｏ　ｄｉｒｅｃｔ　ｔｈｅ ｗｅｌｌ　ｃｏｍｐｌｅｔｉｏｎ．Ｇａｓ　ｓｈａｌｅ　ｉｓ　ｃｈａｒａｃｔｅｒｉｚｅｄ　ｂｙ　ｉｔｓ　ｅｘｔｒｅｍｅｌｙ　ｌｏｗ　ｐｏｒｏｓｉｔｙ，ｅｘｔｒｅｍｅｌｙ　ｌｏｗ　ｐｅｒｍｅａｂｉｌｉｔｙ，ｓｔｒｏｎｇ ｈｅｔｅｒｏｇｅｎｅｉｔｙ　ａｎｄ　ａｎｉｓｏｔｒｏｐｙ，ｔｈｕｓ　ｉｔ　ｉｓ　ｕｎｓｕｉｔａｂｌｅ　ｆｏｒ　ｗｅｌｌ　ｌｏｇｇｉｎｇ　ｉｎｔｅｒｐｒｅｔａｔｉｏｎ　ｏｆ　ｓｈａｌｅ　ｇａｓ　ａｐｐｌｙｉｎｇ　ｉｎｔｅｒｐｒｅｔａｔｉｏｎ ｍｅａｎｓ　ｏｆ　ｃｏｎｖｅｎｔｉｏｎａｌ　ｒｅｓｅｒｖｏｉｒ，ａｎｄ　ｉｔ　ｉｓ　ｎｅｃｅｓｓａｒｙ　ｔｏ　ｒｅｃｏｎｓｔｒｕｃｔ　ａ　ｎｅｗ　ｓｙｓｔｅｍ　ｆｏｒ　ｗｅｌｌ　ｌｏｇｇｉｎｇ　ｉｎｔｅｒｐｒｅｔａｔｉｏｎ　ｏｆ ｓｈａｌｅ　ｇａｓ．Ｅｖａｌｕａｔｉｎｇ　ｔｈｅ　ｒｅｓｅｒｖｏｉｒ　ｐｏｔｅｎｔｉａｌ　ｏｆ　ａ　ｇａｓ　ｓｈａｌｅ　ｉｎｖｏｌｖｅｓ　ｗｅｉｇｈｉｎｇ　ｐｏｓｉｔｉｖｅ　ｏｒ　ｎｅｇａｔｉｖｅ　ｃｏｎｔｒｉｂｕｔｉｏｎｓ　ｆｒｏｍ ａ　ｖａｒｉｅｔｙ　ｏｆ　ｆａｃｔｏｒｓ，ｉｎｃｌｕｄｉｎｇ　ｓｈａｌｅ　ｍｉｎｅｒａｌｏｇｙ　ａｎｄ　ｔｅｘｔｕｒｅ，ｔｈｅ　ｃｏｎｔｅｎｔ　ａｎｄ　ｔｙｐｅｓ　ｏｆ　ｃｌａｙ，ｋｅｒｏｇｅｎ　ｔｙｐｅ　ａｎｄ ｍａｔｕｒｉｔｙ，ｆｌｕｉｄ　ｓａｔｕｒａｔｉｏｎ，ａｄｓｏｒｂｅｄ　ａｎｄ　ｆｒｅｅ　ｇａｓ　ｓｔｏｒａｇｅ　ｍｅｃｈａｎｉｓｍｓ，ｂｕｒｉａｌ　ｄｅｐｔｈ，ｔｅｍｐｅｒａｔｕｒｅ　ａｎｄ　ｐｏｒｅ　ｐｒｅｓｓｕｒｅ， ｅｔｃ．Ｉｎ　ｐａｒｔｉｃｕｌａｒ，ｐｏｒｏｓｉｔｙ，ｏｒｇａｎｉｃ　ｃｏｎｔｅｎｔ　ａｎｄ　ｇａｓ　ｃｏｎｔｅｎｔ　ａｒｅ　ｉｍｐｏｒｔａｎｔ　ｆｏｒ　ｄｅｔｅｒｍｉｎｉｎｇ　ｗｈｅｔｈｅｒ　ａ　ｓｈａｌｅ　ｓｈｏｗｓ ｐｏｔｅｎｔｉａｌ　ｆｏｒ　ｆｕｒｔｈｅｒ　ｄｅｖｅｌｏｐｍｅｎｔ．Ｂａｓｅｄ　ｏｎ　ｔｈｅ　ｃｏｍｐｒｅｈｅｎｓｉｖｅ　ｒｅｓｅａｒｃｈ　ｏｆ　ｅｘｐｌｏｒａｔｉｏｎ　ａｎｄ　ｅｘｐｌｏｉｔａｔｉｏｎ　ｏｆ　ｓｈａｌｅ　ｇａｓ ａｂｒｏａｄ　ｒｅｃｅｎｔｌｙ，ｅｓｐｅｃｉａｌｌｙ　Ａｍｅｒｉｃａ，ｔｈｉｓ　ｐａｐｅｒ　ｓｕｍｍａｒｉｚｅｓ　ｔｈｅ　ｐｒｅｓｅｎｔ　ｓｉｔｕａｔｉｏｎ　ｏｆ　ｓｈａｌｅ　ｇａｓ　ｅｘｐｌｏｒａｔｉｏｎ　ａｎｄ ｅｘｐｌｏｉｔａｔｉｏｎ　ａｐｐｌｙｉｎｇ　ｇｅｏｐｈｙｓｉｃａｌ　ｌｏｇｇｉｎｇ　ａｎｄ　ｍａｄｅ　ａｎ　ａｎａｌｙｓｉｓ　ｏｆ　ｄａｔａ　ａｃｑｕｉｓｉｔｉｏｎ　ｓｅｒｉｅｓ　ｏｆ　ｗｅｌｌ　ｌｏｇｇｉｎｇ　ｉｎ　ｔｈｅ　ｖｉｅｗ　ｏｆ 　４期 郝建飞，等：页岩气地球物理测井评价综述 ｄｉｆｆｅｒｅｎｔ　ｓｔａｇｅ　ｏｆ　ｅｘｐｌｏｒａｔｉｏｎ．Ａｓ　ｔｏ　ｗｅｌｌ　ｌｏｇｇｉｎｇ　ｉｄｅｎｔｉｆｉｃａｔｉｏｎ　ｏｆ　ｓｈａｌｅ　ｇａｓ，ｔｈｉｓ　ｐａｐｅｒ　ｓｕｍｍａｒｉｚｅｓ　ｔｈｅ　ｌｏｇｇｉｎｇ ｒｅｓｐｏｎｓｅ　ｏｆ　ｇａｓ　ｓｈａｌｅ．Ｔｈｅｎ　ｔｈｅ　ｐａｐｅｒ　ｄｉｓｃｕｓｓｅｄ　ｓｏｍｅ　ｉｍｐｏｒｔａｎｔ　ｐａｒａｍｅｔｅｒｓ　ａｎｄ　ｉｔｓ　ｅｖａｌｕａｔｉｏｎ　ｍｅｔｈｏｄｓ　ｏｆ　ｓｈａｌｅ　ｇａｓ ｕｓｉｎｇ　ｇｅｏｐｈｙｓｉｃａｌ　ｌｏｇｇｉｎｇ　ｉｎ　ｄｅｔａｉｌ，ｉｎｃｌｕｄｉｎｇ　ｔｏｔａｌ　ｏｒｇａｎｉｃ　ｃａｒｂｏｎ，ｍｉｎｅｒａｌ　ｔｙｐｅｓ　ａｎｄ　ｃｏｎｔｅｎｔ　ｏｆ　ｆｏｒｍａｔｉｏｎ，ｐｏｒｏｓｉｔｙ， ｇａｓ　ｃｏｎｔｅｎｔ（ｆｒｅｅ　ｇａｓ　ａｎｄ　ａｄｓｏｒｂｅｄ　ｇａｓ）ａｎｄ　ｓｏｍｅ　ｒｏｃｋ　ｍｅｃｈａｎｉｃａｌ　ｐａｒａｍｅｔｅｒｓ．Ｉｎ　ｔｈｅ　ｌａｓｔ，ｔｈｉｓ　ｐａｐｅｒ　ｐｕｔｓ　ｆｏｒｗａｒｄ ｓｅｖｅｒａｌ　ｑｕｅｓｔｉｏｎｓ　ａｂｏｕｔ　ｔｈｅ　ｒｅｓｅａｒｃｈ　ｏｆ　ｇｅｏｐｈｙｓｉｃａｌ　ｌｏｇｇｉｎｇ　ｆｏｒ　ｓｈａｌｅ　ｇａｓ　ｉｎ　Ｃｈｉｎａ，ｔｈｅｎ　ｄｉｓｃｕｓｓｅｓ　ｔｈｅ　ｄｅｖｅｌｏｐｍｅｎｔ ｔｒｅｎｄ　ｏｆ　ｇｅｏｐｈｙｓｉｃａｌ　ｌｏｇｇｉｎｇ． Ｋｅｙｗｏｒｄｓ　ｓｈａｌｅ　ｇａｓ，Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ　ｏｆ　ｓｈａｌｅ　ｇａｓ　ｂｙ　ｗｅｌｌ　ｌｏｇｇｉｎｇ，Ｅｖａｌｕａｔｉｏｎ　ｏｆ　ｓｈａｌｅ　ｇａｓ　ｂｙ　ｗｅｌｌ　ｌｏｇｇｉｎｇ，Ｌｏｇｇｉｎｇ ｔｅｃｈｎｉｑｕｅｓ ０　引　言 页岩气是指生成、储集和封盖（有时）均发生于 页岩体系中，或游离于基质孔隙和天然裂缝中，或吸 附于有机质和粘土矿物表面，或溶解于沥青和水中 （很少），在一定地质条件下聚集成藏并具有商业价 值的生物成因和／或热解成因的天然气［１－２１］．这里说 的页岩不等同于地质概念上的页岩，它也包括富有 机质页岩所夹的粉砂岩与细砂岩．地球物理测井技 术是页岩气勘探开发中重要的技术之一，可以对含 气页岩进行识别，对其储层品质进行有效评价，并为 后期的完井作业提供重要的岩石力学等参数． 目前，我国页岩气地球物理测井技术仅仅处在 起步阶段，工作重点主要放在跟踪国外发展动态和 与国外测井服务公司合作对页岩气进行评价方面． ２００９年，中国石油第一口页岩气井 Ｗ２０１井的开 钻，标志着中国页岩气进入实质性勘探开发阶段，该 井于２０１０年２月７日由斯伦贝谢公司进行测井作 业，以斯伦贝谢为主，国内测井服务公司和研究机构 参与进行了测井解释评价． 目前，我国页岩气地球物理测井技术研究缺乏 系统性，还不能够满足页岩气勘探开发的需求．要进 一步加快我国页岩气的勘探开发步伐，需要测井工 作者加强页岩气岩石物理等基础研究工作，紧跟国 际测井发展趋势并有所创新，走出一条适合中国页 岩气勘探开发地球物理测井技术的新路． 页岩气藏无论在储层特征、气藏特征还是在成 藏机制方面都与常规油气藏有着很大的区别，完井 方式及其之后的开发方式亦有相当大的差别，由此 决定了页岩气藏对测井技术的需求与常规油气藏有 着很大不同［２２］． １　页岩气测井采集系列及其响应特征 １．１　测井采集系列 测井在页岩气藏勘探开发中有两大任务，一是 储层及含气量的评价，二是为完井服务提供指导参 数并在钻井中起地质导向作用，这其中包含了岩性、 孔隙、裂缝、有机碳、储层岩石力学等参数评价．勘探 和开发不同阶段达到上述目的采用的测井系列是不 同的，表１总结了国外针对不同井别采用的测井采 集系列．对于新区，一般而言，最经济的测井系列包 括自然伽马测井、自然电位测井、井径、岩性密度测 井、补偿中子测井、电阻率测井（双侧向或者阵列感 应测井）、元素俘获能谱测井和声波扫描测井． 从表１中可见，除了一些常规油气藏采用的测 井方法，在页岩气测井采集中还采用了一些测井新 技术，包括元素俘获能谱测井、核磁共振测井、微电 阻率成像测井和声波扫描测井，这些测井新技术的 应用在页岩气勘探开发的初期是非常有必要的，有 助于含气页岩储层特征的综合评价，也有助于指导 油公司后续的勘探开发，例如，运用微电阻率成像测 表１　页岩气不同井别采用的测井采集系列 Ｔａｂｌｅ　１　Ｌｏｇｇｉｎｇ　Ｓｅｒｉｅｓ　ｏｆ　Ｄｉｆｆｅｒｅｎｔ　Ｔｙｐｅｓ　ｏｆ　Ｓｈａｌｅ　Ｇａｓ　Ｗｅｌｌ 完井 方式 井别 必测井项目 可选测井项目 裸 　 眼 　 井 探井 自然伽马能谱测井 岩性密度测井 补偿中子测井 双侧向（阵列感应）测井 元素俘获能谱测井 微电阻率扫描成像测井 声波扫描测井 自然电位测井 井径、井斜、井温 核磁共振测井 评价井 自然伽马能谱测井 岩性密度测井 补偿中子测井 双侧向（阵列感应）测井 偶极声波测井 自然电位测井 井径、井斜、井温 核磁共振测井 元素俘获能谱测井 开发井 自然伽马能谱测井 岩性密度测井 补偿中子测井 阵列声波测井 井径、井斜、井温 微电阻率 扫描成像测井 套管井 自然伽马能谱测井 岩性密度测井 声波扫描测井 相对方位测井 热中子寿命测井 ５２６１ 地球物理学进展　ｈｔｔｐ／／ｗｗｗ．ｐｒｏｇｅｏｐｈｙｓ．ｃｎ　 ２７卷　 井、声波全波列测井和井下声波电视可以确定裂缝 性质［２３－２６］；用元素俘获能谱测井能够确定岩石矿物 含量并计算有机碳和无机碳含量［２７］． １．２　页岩气的测井响应特征 页岩气测井曲线呈现“三高两低”的特征，即高 自然伽马、高电阻率、高中子、低密度（相对普通页岩 而言）、低光电吸收截面（相对普通页岩而言）［２８－３０］． 图１为北美Ｂａｒｎｅｔｔ页岩某井的典型测井图． 图１　Ｂａｒｎｅｔｔ页岩某井典型测井图 Ｆｉｇ．１　Ｔｙｐｉｃａｌ　Ｓｈａｌｅ　Ｇａｓ　Ｌｏｇ　ｏｆ　Ｂａｒｎｅｔｔ　Ｓｈａｌｅ ２　页岩气测井评价 ２．１　总有机碳（ＴＯＣ）含量计算 总有机碳含量是页岩气评价一个很重要的参 数，它从侧面反映含气页岩中有机质含量的多少和 生烃潜力大小［３１］．对于可开发含气页岩，总有机碳 含量下限值是２％．Ｐａｓｓｅｙ等（１９９０）提出用电阻率 和孔隙度曲线叠加的方法来评价总有机碳含量，即 ΔｌｇＲ方法［３２－３４］；Ｄａｎｉｅｌ　Ｒｏｓｅ等（２００８）提出用密度 测井和自然伽马测井来对总有机碳含量进行评 价［３５］；Ｄａｎｉｅｌ　Ｃｏｏｐｅ等（２００９）提出了一种有限制的 矿物模型，在这种模型下运用元素俘获能谱测井和 核磁共振测井来计算总有机碳含量并得到有效应 用［３６］；Ｊａｃｏｂｉ等（２００９）等提出用密度测井和核磁共 振测井来计算有机碳含量［３７］；Ｒｉｃｈａｒｄ等（２００９）提 出用脉冲中子来计算有机碳含量［３８］．随着测井新技 术的推广应用，国外采用元素俘获能谱测井（ＥＣＳ） 对总有机碳含量进行计算，这种方法将干酪根作为 一种矿物，先计算干酪根的体积百分比，然后再通过 干酪根和总有机碳含量的关系来计算总有机碳含 量．对于以上几种方法而言，ΔｌｇＲ方法相对简单，元 素俘获能谱测井方法相对精确，在没有元素俘获能 谱测井数据的前提下，用ΔｌｇＲ 方法也能得到有效 的结果． ΔｌｇＲ方法的理论基础来源于阿尔奇公式、威利 时间公式及总有机碳含量与 ΔｌｇＲ 关系的经验公 式，具体方法如下： １）用自然伽马曲线和自然电位曲线剔除油层、 蒸发岩、火成岩、低空层段、欠压实的沉积物和井壁 垮塌严重的井段． ２）将刻度合适的孔隙度曲线（以声波时差曲线 为例）叠加在电阻率曲线上．在刻度时，声波时差曲 线刻度从左至右减小，右侧为０，左侧刻度为５０μｓ／ ｆｔ（１６４μｓ／ｍ）的倍数．电阻率从左至右增大．在叠加 时，应保持声波曲线不动，对电阻率曲线刻度进行变 化，但应保持电阻率曲线每两个数量级对应声波时 差为１００μｓ／ｆｔ（３２８μｓ／ｍ）间隔的原则．２条曲线在 一定的深度范围“一致”或完全重叠，重叠段即为基 线段． ３）计算ΔｌｇＲ的公式为 ΔｌｇＲ＝ｌｇ ＲＲ（ ）基线 ＋０．０２×（Δｔ－Δｔ基线）， （１） 式中，ΔｌｇＲ为电阻率曲线和孔隙度曲线幅度差在对 数坐标上数值；Ｒ为电阻率值，Ω·ｍ；Δｔ为声波时 差值，μｓ／ｆｔ；Ｒ基线 为基线处的电阻率值，Ω·ｍ； Δｔ基线为基线处的声波时差值，μｓ／ｆｔ． ４）利用ΔｌｇＲ和总有机碳含量的经验方程来计 算总有机碳含量的方法为 ＴＯＣ＝（ΔｌｇＲ）×１０（２．２９７－０．１６８８×Ｒｏ）， （２） 式中，Ｒｏ 为含气页岩的成熟度． 应该注意，在选择孔隙度曲线时，用声波时差与 电阻率曲线叠加比密度或者中子与电阻率曲线叠加 计算的总有机碳含量精度更高，因为井眼条件会影 响到密度和中子测井值．另外，ΔｌｇＲ方法对于过成 熟的含气页岩，计算得到的总有机碳含量偏低［３９］． ２．２　矿物含量计算 对于含气页岩，矿物含量计算显得尤为重要，因 ６２６１ 　４期 郝建飞，等：页岩气地球物理测井评价综述 为它直接影响着后期压裂／射孔井段的选择，其中脆 性矿物（石英、方解石等）和粘土矿物含量是很重要 的两个方面，在压裂时，应当选取粘土矿物含量相对 少、脆性矿物含量相对多的层段．应用常规测井方法 对矿物含量准确计算非常困难，近几年，随着元素俘 获测井（ＥＣＳ）的推广应用，计算矿物含量方法逐渐 成熟．特别是元素俘获能谱测井既可以应用于裸眼 井中，又可以在套管井中应用，而且不受井眼条件的 影响，这大大加大了其应用的空间． 元素俘获测井利用ＢＧＯ晶体探测器测量快中 子与地层中元素发生非弹性碰撞产生的伽马能谱， 经过解谱可以得到Ｃ、Ｏ、Ｓｉ、Ｃａ等元素的相对产额； 而对其中主要的俘获伽马能谱经过解谱处理可以得 到 Ｈ、Ｃｌ、Ｓｉ、Ｃａ、Ｓ、Ｋ、Ｆｅ、Ｔｉ和 Ｇｄ等元素相对产 额．有了元素的相对产额，通过特定的氧化物闭合模 型就能得到这些元素的百分含量．再通过元素丰度 和矿物含量（岩心实验）之间的统计关系得到方解 石、白云石、石英、长石、云母、硬石膏、黄铁矿等矿物 的含量．最后将元素俘获能谱测井与常规测井曲线 结合确定粘土矿物含量及类型［３７，４０，４１］． ２．３　孔隙度计算 孔隙度是计算页岩气中游离气含量需要的重要 参数之一，目前国外大部分可开发的页岩气藏孔隙 度下限值为２．对于含气页岩，孔隙存在于粒间孔 隙、粘土矿物、干酪根和微裂缝中及其周围［４２］．由于 含气页岩的孔隙度非常小，常规测井孔隙度计算方 法已不再适用，斯伦贝谢公司（２０１０）称核磁共振测 井在页岩气孔隙度计算中有着较好的应用；Ｍｉｃｈａｅｌ 等（２００２）年提出用元素俘获能谱测井对孔隙度进行 计算［４３］． 首先，通过４００多块样品的矿物分析和地球化 学实验，得到骨架密度值、中子值与元素含量的关系 式为 ρｍａ＝２．６２０＋０．０４９０ＷＳｉ－０．２２７４ＷＣａ＋ １．９９３ＷＦｅ＋１．１９３ＷＳ， （３） Ｎｍａ＝０．４０８－０．８８９ＷＳｉ－１．０１４ＷＣａ－ ０．２５７ＷＦｅ＋０．６７５ＷＳ， （４） 式中，ρｍａ为骨架密度值；Ｎｍａ为骨架中子值；ＷＳｉ为硅 元素含量百分比；ＷＣａ为钙元素含量百分比；ＷＦｅ为 铁元素含量百分比；ＷＳ为硫元素含量百分比． 然后应用计算得到的骨架密度值和中子值来计 算密度孔隙度、中子孔隙度 ｔｄ＝ρｍａ －ρｂ ρｍａ－ρｆ ， （５） ｔｎ＝ Ｎｍａ－ｎ Ｎｍａ－Ｎｆ ， （６） 式中，ｔｄ为密度算的总孔隙度；ｔｎ为中子算的总孔 隙度；ρｂ 测得密度值；ρｆ 为泥浆密度值；ｎ 为测得中 子值；Ｎｆ 为泥浆的中子值． 最后取页岩的总孔隙度为 ｔ＝ｍｉｎ（ｔｄ，０．６６６７ｔｄ＋０．３３３３ｔｎ）． （７） ２．４　含气量的计算 在含气页岩中，气分为三部分，游离气、吸附气 和溶解气，由于溶解气含量较少，而且与吸附气不宜 区分，因此通常将溶解气也归于吸附气［４４－４７］计算． Ｄａｖｉｄ等（２００４）通过实验得到总有机碳含量和总气 体含量之间的关系，通过总有机碳含量来计算总含 气量．一般说来，应该将吸附气和游离气含量分开 计算． ２．４．１　吸附气含量计算 Ｒｉｃｋ等（２００４）提出应用兰格缪尔（Ｌａｎｇｍｕｉｒ） 等温线（如图２）来计算含气页岩吸附气含量［４８，４９］； Ｄａｎｉｅｌ　Ｒｏｓｅ等（２００８）提出通过实验得到油藏温度 和吸附气含量之间的经验公式来计算吸附气含量． 兰格缪尔等温线是用来描述在特定温度下吸附在干 酪根表面上的甲烷与以气态形式的甲烷的一种平衡 状态．通过对含气页岩做等温吸附实验可以得到兰 格缪尔等温线［５０］，通常，对于一个盆地而言，只有一 条兰格缪尔曲线能够充分描述该盆地的含气页岩． 图２　兰格缪尔等温线 Ｆｉｇ．２　Ｌａｎｇｍｕｉｒ’ｓ　Ｉｓｏｔｈｅｒｍ 在已知岩心总有机碳含量下，可以通过等温吸 ７２６１ 地球物理学进展　ｈｔｔｐ／／ｗｗｗ．ｐｒｏｇｅｏｐｈｙｓ．ｃｎ　 ２７卷　 附实验得到兰格缪尔等温线．由于兰格缪尔曲线是 在特定的温度和有机碳含量的情况下得到的，因此 当将兰格缪尔等温线应用到测井评价上时，需要对 总有机碳含量、温度和压力做校正． 对于温度和压力校正，方法为 Ｖｌｔ＝１０（－ｃ１·（Ｔ＋ｃ２））， （８） Ｐｌｔ＝１０（ｃ３·（Ｔ＋ｃ４））， （９） ｃ４＝ｌｏｇＰｌ＋（－ｃ３·Ｔｉ）， （１０） ｃ２＝ｌｏｇＶｌ＋（ｃ１·Ｔｉ）， （１１） 式中，Ｖｌｔ为在油藏温度下的兰格缪尔体积（ｓｃｆ／ ｔｏｎ）；Ｐｌｔ为油藏温度下的兰格缪尔压力（ｐｓｉａ）；ｃ１为 常数，０．００２７；ｃ３为常数，０．００５；ｃ２和ｃ４为中间过 渡变量；Ｔ为油藏温度，℃；Ｔｉ为等温线的温度，℃． 然后对ＴＯＣ校正，方法为 Ｖｌｃ＝Ｖｌｔ· ＴＯＣｌｇ ＴＯＣｉｓｏ ， （１２） 式中，Ｖｌｃ为在油藏温度下经过ＴＯＣ校正后的兰格 缪尔体积；ＴＯＣｉｓｏ为得到兰格缪尔等温线时采用的 ＴＯＣ值；ＴＯＣｌｇ为测井得到的ＴＯＣ值． 最后，可以经过温度、压力和ＴＯＣ校正后的兰 格缪尔公式来计算吸附气含量 ｇｃ＝ Ｖｌｃｐ （ｐ＋Ｐｌｔ） ， （１３） 式中，ｇｃ 为吸附气含量（ｓｃｆ／ｔｏｎ）；ｐ 为油藏压力 （ｐｓｉａ）；Ｐｌｔ为油藏温度下的兰格缪尔压力（ｐｓｉａ）；Ｖｌｃ 为在油藏温度下经过 ＴＯＣ 校正后的兰格缪尔 体积． 值得注意的是，由兰格缪尔等温线计算得到的 气的含量是页岩能够容纳的吸附气含量的体积，而 不是页岩中所含有的吸附气的体积，即如果含气页 岩中气体出现逃逸现象，用兰格缪尔曲线得到的结 果会偏大，但页岩本身也是盖层，如果其裂缝不发 育，用这种方法计算吸附气含量效果较好． ２．４．２　游离气含量计算 游离气含量的大小取决于气藏的压力、有效孔 隙度和含气饱和度［４８］．游离气含量计算法为 Ｇｃｆｍ＝１Ｂｇ ·（ｅｆｆ（１－Ｓｗ））·ψ ρｂ ， （１４） Ｂｇ＝ ｐｓｃｚｓｃ（Ｔｓｃ＋４５９．６７［ ］）· ｚ （Ｔ＋４５９．６７）［ ］ｐ ， （１５） 式中，Ｇｃｆｍ 为游离气体积，ｓｃｆ／ｔｏｎ；Ｂｇ 为气体的地 层体积因素，ｃｆ（油藏）／ｓｃｆ；ｅｆｆ为有效孔隙度，ｖ／ｖ； Ｓｗ 为含水饱和度，ｖ／ｖ；ρｂ 为体积密度，ｇ／ｃｍ ３；ψ为 常数，３２．１０５２；ｐ为气藏压力，ｐｓｉａ；Ｔ为油藏温度， ℉；ｚ为气藏下的气体压缩因子，为一个无量纲参 数，可以通过实验来获得；ｐｓｃ为 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 状况下的压力， １４．６９６ｐｓｉａ；ｚｓｃ为标准状况下的温度，３２℉． ２．５　岩石力学参数计算 前已述及，除了对含气页岩进行储层参数评价 外，页岩气地球物理测井另外一个重要的任务就是 对水平井着陆点和后期的完井作业提供有效的指导 参数，来指导水平井着陆点和压裂／射孔井段的 选择． 国外测井服务公司，都采用阵列声波（斯伦贝谢 公司采用更先进的声波扫描测井）来评价岩石力学 特征，但评价方式及参数有所区别．贝克休斯和哈里 伯顿公司应用纵横波时差来计算岩石的杨氏模量和 泊松比，然后应用杨氏模量和泊松比来计算岩石脆 性程度［５１，５２］，通过岩石脆性程度来选择水平着陆点 和压裂／射孔层段；斯伦贝谢公司通过应用纵横波时 差和最小水平应力的各向异性模型来计算最小水平 应力，通过最小水平应力来优选水平着陆点和压裂／ 射孔层段． 页岩具有层状特征，造成其横向和纵向上的杨 氏模量和泊松比变化很大，具有很强的非均质性和 各向异性．与岩石脆性程度和应力的各向同性模型 相比，最小应力各向异性模型不仅考虑了岩石的纵 向各向异性而且也考虑了水平方向的各向异性，因 而其可靠性更强［４８，５３－５７］． 最小水平应力各向异性模型为 σＨｍｉｎ＝ ＥＨ Ｅｖ νｖ １－νＨ （σｖ－αｐｒ）＋αｐｒ， （１６） 式中，σＨｍｉｎ为最小水平应力，ｐｓｉ；α为比奥系数；ｐｒ 为孔隙压力，ｐｓｉａ；ＥＨ 为水平方向上的杨氏模量；Ｅｖ 为垂直方向上的杨氏模量；νｖ为垂直方向上的泊松 比；νＨ 为水平方向上的泊松比；σｖ 为垂直方向上的 应力，ｐｓｉａ．对于含气页岩，水平着陆点和后期压裂／ 射孔井段都应选择在应力低的地方． 最小水平应力模型的适用条件为： １）横向异性介质； ２）具有上覆岩层压力； ３）没有构造应力； ４）需要知道水平和垂直方向的杨氏模量和泊 松比． ３　我国页岩气地球物理测井研究存在的问 题及发展趋势 　　目前，对于地球物理测井而言，页岩气的岩石物 ８２６１ 　４期 郝建飞，等：页岩气地球物理测井评价综述 理实验研究和测井评价处于对国外页岩气岩石物理 实验和测井评价方法进行调研学习阶段，而且由于 我国地球物理测井技术与国外比较而言相对薄弱， 因此在立足我国国情下采用经济有效的方法对页岩 气进行评价非常困难． 面对这种困难，国内各大油公司也相继设立研 究课题对页岩气的测井评价方法进行研究，２０１１年 中国石油设立课题对页岩气的岩石物理特征和测井 评价方法进行研究，中国石化和中国海洋石油也相 继设立此类研究课题．在地球物理测井评价方面，首 先应该通过页岩的岩石物理及地球化学实验来对页 岩气产生和存储机理进行研究，然后采用经济有效 的测井系列和解释评价方法对页岩的总有机碳、孔 隙度和含气量进行评价，优选页岩的有利层段，最后 计算岩石的力学参数来指导后期压裂作业．针对目 前的研究现状，我认为我国的页岩气地球物理测井 研究存在如下问题和发展趋势． ３．１　我国页岩气地球物理测井研究存在的问题 我国页岩气的勘探开发和研究工作总体处于起 步状态，含气页岩的岩石物理实验分析技术和测井 评价关键技术属于空白，严重制约了我国页岩气地 球物理测井评价工作，目前我国页岩气地球物理测 井研究存在的主要问题有以下几方面： １）页岩气藏岩石物理基础研究和实验分析技术 尚属空白．岩石物理基础研究是实现页岩气合理准 确评价的基础，而其实验分析技术是基础中的基础． 国外针对含气页岩的岩石物理实验分析技术已经相 对成熟，而国内刚刚开始启动基础研究工作，对国外 先进的岩石物理实验分析技术尚不能掌握，也缺乏 与页岩气储层相适应的极低孔渗岩石的实验手段． 油公司首先要优先考虑页岩气藏测井技术人才的配 备和前期研究力度，同时应对岩石物理实验分析技 术加大投入，购买或者研发适用于含气页岩的岩石 物理实验设备，更好的服务于含气页岩的测井解释 评价． ２）有中国特色的页岩气测井系列还不明确．从 美国众多页岩气藏的分析可知，页岩气储层除了具 有共性特征之外，各个气藏的个性特征还是相差较 大．中国页岩气藏勘探刚刚起步，相信中国页岩气藏 和美国气藏以及中国不同地区气藏之间都会存在一 定差异性，如何根据气藏的岩石物理特征以及不同 的页岩气开采方式，采用不同的有针对性的测井系 列，甚至研究有针对性的测井方法和研发相应的测 量仪器，在借鉴国外测井成功经验的基础上走出自 己的技术之路，是一个必须探索的重要问题． ３）页岩气测井解释方法和软件缺乏．国外含气 页岩的测井解释评价日趋完善，各测井服务公司都 已经有了成型的软件系统，贝克休斯公司发展了一 套专门应用于页岩气评价的软件—页岩气专家系 统；斯伦贝谢也研发出了应用于解释评价的页岩气 模块．页岩气藏与常规油气藏在储层物性、电性等方 面有着很大的不同，因此不能将常规油气藏的测井 解释评价体系直接应用于含气页岩，对于含气页岩 的解释评价，需要在岩石物理研究的基础上，在借鉴 国外技术和经验的同时，建立一套适用于含气页岩 的测井处理解释评价体系．同时应该考虑，国外页岩 气的测井方法及测井解释评价技术是否适用于我国 含气页岩有待于进一步验证． ３．２　页岩气地球物理测井发展趋势 １）随钻测井的应用将大大推动页岩气的勘探开 发．页岩气一般采用水平井钻井，因为水平井井眼与 页岩接触面积更大，并且增加了相互交叉的多样的 裂缝的可能性［２５］，这一方面限制了电缆测井的应 用，另一方面推动了随钻测井在其中的发展．随钻测 井在水平井钻井中起着不可替代的地质导向作用， 而且随钻测井资料不受泥浆侵入的影响，资料真实 可靠并且更加经济． ２）井中和井间地球物理技术的结合将为页岩气 储层描述和完井质量提供更有效的评价．在页岩气 勘探开发中，阵列声波和声波扫描测井等井中地球 物理技术发挥了很重要的作用．若将井中和井间地 球物理技术（如井间微地震等）结合起来，必然会对 页岩气储层分布、品质及后期完井质量得到更为精 确的解释评价． ３）测井新技术在页岩气的定量评价中将有着更 为广阔的前景．应用常规方法对页岩气进行定量解 释评价非常困难，而低孔低渗、复杂岩性等储层的测 井实践证明现代测井技术是解决问题的主要手 段［５８］，因此也有理由相信测井新技术（如核磁共振 测井、微电阻率成像测井、元素俘获能谱测井等）在 页岩气定量解释评价方面应该有着更为广阔的 前景． ４　结论与建议 ４．１　加强页岩气测井基础理论研究，全面系统地开 展页岩气储层地球物理测井评价方法研究，来指导 地球物理测井在页岩气测井识别及解释评价中的 应用． ９２６１ 地球物理学进展　ｈｔｔｐ／／ｗｗｗ．ｐｒｏｇｅｏｐｈｙｓ．ｃｎ　 ２７卷　 ４．２　加快适用于页岩气的测井仪器的研发进程，并 引进先进的测井方法，服务于我国页岩气的勘探 开发． ４．３　页岩气储层具有非均质性和各项异性强的特 点，为满足页岩气勘探开发对地球物理测井技术提 出的要求，针对不同完井类型、井别及测井成本等情 况，优化不同的测井方法，为页岩气勘探开发建立强 有力的地球物理测井技术系列支持． ４．４　尽快建立起页岩气储层测井综合解释评价体 系、实现解释评价的定量化，并研发专门应用于页岩 气的测井解释评价软件． 参　考　文　献 （Ｒｅｆｅｒｅｎｃｅｓ）： ［１］　李新景，胡素云，程克明．北美裂缝性页岩气勘探开发的启 示．石油勘探与开发，２００７，３４（４）：３９２－４００． Ｌｉ　Ｘ　Ｊ，Ｈｕ　Ｓ　Ｙ，Ｃｈｅｎｇ　Ｋ　Ｍ．Ｓｕｇｇｅｓｔｉｏｎｓ　ｆｒｏｍ　ｔｈｅ ｄｅｖｅｌｏｐｍｅｎｔ　ｏｆ　ｆｒａｃｔｕｒｅｄ　ｓｈａｌｅ　ｇａｓ　ｉｎ　Ｎｏｒｔｈ　Ａｍｅｒｉｃａ． Ｐｅｔｒｏｌｅｕｍ　Ｅｘｐｌｏｒａｔｉｏｎ　ａｎｄ　Ｄｅｖｅｌｏｐｍｅｎｔ（ｉｎ　Ｃｈｉｎｅｓｅ），２００７， ３４（４）：３９２－４００． ［２］　叶军，曾华盛．川西须家河组泥页岩气成藏条件与勘探潜力． 天然气工业，２００８，２８（１２）：１８－２５． Ｙｅ　Ｊ，Ｚｅｎｇ　Ｈ　Ｓ．Ｐｏｏｌｉｎｇ　ｃｏｎｄｉｔｉｏｎｓ　ａｎｄ　ｅｘｐｌｏｒａｔｉｏｎ　ｐｒｏｓｐｅｃｔ ｏｆ　Ｓｈａｌｅ　Ｇａｓ　ｉｎ　Ｘｕｊｉａｈｅ　ｆｏｒｍａｔｉｏｎ　ｉｎ　ｗｅｓｔｅｒｎ　Ｓｉｃｈｕａｎ ｄｅｐｒｅｓｓｉｏｎ．Ｎａｔｕｒａｌ　Ｇａｓ　Ｉｎｄｕｓｔｒｙ （ｉｎ　Ｃｈｉｎｅｓｅ），２００８，２８ （１２）：１８－２５． ［３］　程克明，王世谦，董大忠等．上扬子区下寒武统筇竹寺组页岩 气成藏条件．天然气工业，２００９，２９（５）：４０－４４． Ｃｈｅｎｇ　Ｋ　Ｍ，Ｗａｎｇ　Ｓ　Ｑ，Ｄｏｎｇ　Ｄ　Ｚ，ｅｔ　ａｌ．Ａｃｃｕｍｕｌａｔｉｏｎ ｃｏｎｄｉｔｉｏｎｓ　ｏｆ　ｓｈａｌｅ　ｇａｓ　ｒｅｓｅｒｖｏｉｒｓ　ｉｎ　ｔｈｅ　Ｌｏｗｅｒ　Ｃａｍｂｒｉａｎ Ｑｉｏｎｇｚｈｕｓｉ　ｆｏｒｍａｔｉｏｎ，ｔｈｅ　Ｕｐｐｅｒ　Ｙａｎｇｔｚｅ　ｒｅｇｉｏｎ．Ｎａｔｕｒａｌ　Ｇａｓ Ｉｎｄｕｓｔｒｙ（ｉｎ　Ｃｈｉｎｅｓｅ），２００９，２９（５）：４０－４４． ［４］　董大忠，程克明，王世谦等．页岩气资源评价方法及其在四川 盆地的应用．天然气工业，２００９，２９（５）：３３－３９． Ｄｏｎｇ　Ｄ　Ｚ，Ｃｈｅｎｇ　Ｋ　Ｍ，Ｗａｎｇ　Ｓ　Ｑ，ｅｔ　ａｌ．Ａｎ　ｅｖａｌｕａｔｉｏｎ ｍｅｔｈｏｄ　ｏｆ　ｓｈａｌｅ　ｇａｓ　ｒｅｓｏｕｒｃｅ　ａｎｄ　ｉｔｓ　ａｐｐｌｉｃａｔｉｏｎ　ｉｎ　ｔｈｅ　Ｓｉｃｈｕａｎ ｂａｓｉｎ．Ｎａｔｕｒａｌ　Ｇａｓ　Ｉｎｄｕｓｔｒｙ（ｉｎ　Ｃｈｉｎｅｓｅ），２００９，２９（５）：３３－ ３９． ［５］　Ｌｏｕｃｋｓ　Ｒ　Ｇ，Ｒｕｐｐｅｌ　Ｓ　Ｃ．Ｍｉｓｓｉｓｓｉｐｐｉａｎ　Ｂａｒｎｅｔｔ　Ｓｈａｌｅ： Ｌｉｔｈｏｆａｃｉｅｓ　ａｎｄ　ｄｅｐｏｓｉｔｉｏｎａｌ　ｓｅｔｔｉｎｇ　ｏｆ　ａ　ｄｅｅｐ－ｗａｔｅｒ　ｓｈａｌｅ－ｇａｓ ｓｕｃｃｅｓｓｉｏｎ　ｉｎ　ｔｈｅ　Ｆｏｒｔ　Ｗｏｒｔｈ　Ｂａｓｉｎ，Ｔｅｘａｓ．ＡＡＰＧ　Ｂｕｌｌｅｔｉｎ， ２００７，９１（４）：５７９－６０１． ［６］　Ｊａｒｖｉｅ　Ｄ　Ｍ，Ｈｉｌｌ　Ｒ　Ｊ，Ｒｕｂｌｅ　Ｔ　Ｅ，ｅｔ　ａｌ．Ｕｎｃｏｎｖｅｎｔｉｏｎａｌ　ｓｈａｌｅ－ ｇａｓ　ｓｙｓｔｅｍｓ：Ｔｈｅ　Ｍｉｓｓｉｓｓｉｐｐｉａｎ　Ｂａｒｎｅｔｔ　Ｓｈａｌｅ　ｏｆ　ｎｏｒｔｈ－ｃｅｎｔｒａｌ Ｔｅｘａｓ　ａｓ　ｏｎｅ　ｍｏｄｅｌ　ｆｏｒ　ｔｈｅｒｍｏｇｅｎｉｃ　ｓｈａｌｅ－ｇａｓ　ａｓｓｅｓｓｍｅｎｔ． ＡＡＰＧ　Ｂｕｌｌｅｔｉｎ，２００７，９１（４）：４７５－４９９． ［７］　Ｃｕｒｔｉｓ　Ｊ　Ｂ．Ｆｒａｃｔｕｒｅｄ　ｓｈａｌｅ－ｇａｓ　ｓｙｓｔｅｍｓ．ＡＡＰＧ　Ｂｕｌｌｅｔｉｎ， ２００２，８６（１１）：１９２１－１９３８． ［８］　宋海斌，江为为，张文生等．天然气水合物的海洋地球物理研 究进展．地球物理学进展，２００２，１７（２）：２２４－２２９． Ｓｏｎｇ　Ｈ　Ｂ，Ｊｉａｎｇ　Ｗ　Ｗ，Ｚｈａｎｇ　Ｗ　Ｓ，ｅｔ　ａｌ．Ｐｒｏｇｒｅｓｓ　ｏｎ ｍａｒｉｎｅ　ｇｅｏｐｈｙｓｉｃａｌ　ｓｔｕｄｉｅｓ　ｏｆ　ｇａｓ　ｈｙｄｒａｔｅｓ．Ｐｒｏｇｒｅｓｓ　ｉｎ Ｇｅｏｐｈｙｓｉｃｓ（ｉｎ　Ｃｈｉｎｅｓｅ），２００２，１７（２）：２２４－２２９． ［９］　Ｂｏｗｋｅｒ　Ｋ　Ａ．Ｂａｒｎｅｔｔ　Ｓｈａｌｅ　ｇａｓ　ｐｒｏｄｕｃｔｉｏｎ，Ｆｏｒｔ　Ｗｏｒｔｈ　ｂａｓｉｎ： Ｉｓｓｕｅｓ　ａｎｄ　ｄｉｓｃｕｓｓｉｏｎ．ＡＡＰＧ　Ｂｕｌｌｅｔｉｎ，２００４，９１（４）：５２３－５３３． ［１０］　Ｋｉｎｌｅｙ　Ｔ　Ｊ，Ｃｏｏｋ　Ｌ　Ｗ，Ｂｒｅｙｅｒ　Ｊ　Ａ，ｅｔ　ａｌ．Ｈｙｄｒｏｃａｒｂｏｎ ｐｏｔｅｎｔｉａｌ　ｏｆ　ｔｈｅ　Ｂａｒｎｅｔｔ　Ｓｈａｌｅ （Ｍｉｓｓｉｓｓｉｐｐｉａｎ），Ｄｅｌａｗａｒｅ Ｂａｓｉｎ，ｗｅｓｔ　Ｔｅｘａｓ　ａｎｄ　ｓｏｕｔｈ　ｅａｓｔｅｒｎ　Ｎｅｗ　Ｍｅｘｉｃｏ．ＡＡＰＧ Ｂｕｌｌｅｔｉｎ，２００８，９２（８）：９６７－９９１． ［１１］　何雨丹，肖立志，毛志强等．测井评价“三低”油气藏面临的 挑战和发展方向．地球物理学进展，２００５，２０（２）：２８２－２８８． Ｈｅ　Ｙ　Ｄ，Ｘｉａｏ　Ｌ　Ｚ，Ｍａｏ　Ｚ　Ｑ，ｅｔ　ａｌ．Ｔｈｅ　ｐｒｅｓｅｎｔ　ｃｈａｌｌｅｎｇｅｓ ａｎｄ　ｒｅｓｅａｒｃｈ　ｄｉｒｅｃｔｉｏｎ　ｏｆ　ｗｅｌｌ　ｌｏｇｇｉｎｇ　ｅｖａｌｕａｔｉｏｎ　ｍｅｔｈｏｄｓ　ｉｎ ｌｏｗ　ｐｏｒｏｓｉｔｙ，ｌｏｗ　ｐｅｒｍｅａｂｉｌｉｔｙ　ａｎｄ　ｌｏｗ　ｓａｔｕｒａｔｉｏｎ　ｒｅｓｅｒｖｏｉｒ． Ｐｒｏｇｒｅｓｓ　ｉｎ　Ｇｅｏｐｈｙｓｉｃｓ（ｉｎ　Ｃｈｉｎｅｓｅ），２００５，２０（２）：２８２－２８８． ［１２］　Ｐａｓｓｅｙ　Ｑ　Ｒ，Ｂｏｈａｃｓ　Ｋ　Ｍ，Ｅｓｃｈ　Ｗ　Ｌ，ｅｔ　ａｌ．Ｆｒｏｍ　Ｏｉｌ－Ｐｒｏｎｅ Ｓｏｕｒｃｅ　Ｒｏｃｋ　ｔｏ　Ｇａｓ－Ｐｒｏｄｕｃｉｎｇ　Ｓｈａｌｅ　Ｒｅｓｅｒｖｏｉｒ－Ｇｅｏｌｏｇｉｃ　ａｎｄ Ｐｅｔｒｏｐｈｙｓｉｃａｌ　Ｃｈａｒａｃｔｅｒｉｚａｔｉｏｎ　ｏｆ　Ｕｎｃｏｎｖｅｎｔｉｏｎａｌ　Ｓｈａｌｅ　Ｇａｓ Ｒｅｓｅｒｖｏｉｒｓ．ＳＰＥ１３１３５０，２０１０． ［１３］　Ｓｏｎｄｅｒｇｅｌｄ　Ｃ　Ｈ，Ａｍｂｒｏｓｅ　Ｒ　Ｊ，Ｒａｉ　Ｃ　Ｓ，ｅｔ　ａｌ．Ｍｉｃｒｏ－ Ｓｔｒｕｃｔｕｒａｌ　Ｓｔｕｄｉｅｓ　ｏｆ　Ｇａｓ　Ｓｈａｌｅｓ．ＳＰＥ１３１７７１，２０１０． ［１４］　陈波，兰正凯．上扬子地区下寒武统页岩气资源潜力．中国 石油勘探，２００９，（３）：１０－１４． Ｃｈｅｎ　Ｂ，Ｌａｎ　Ｚ　Ｋ．Ｌｏｗｅｒ　ｃａｍｂｒｉａｎ　ｓｈａｌｅ　ｇａｓ　ｒｅｓｏｕｒｃｅ ｐｏｔｅｎｔｉａｌ　ｉｎ　ｕｐｐｅｒ　Ｙａｎｇｔｚｅ　Ｒｅｇｉｏｎ． Ｃｈｉｎａ　Ｐｅｔｒｏｌｅｕｍ Ｅｘｐｌｏｒａｔｉｏｎ（ｉｎ　Ｃｈｉｎｅｓｅ），２００９，（３）：１０－１４． ［１５］　王社教，王兰生，黄金亮等．上扬子区志留系页岩气成藏条 件．天然气工业，２００９，２９（５）：４５－５０． Ｗａｎｇ　Ｓ　Ｊ，Ｗａｎｇ　Ｌ　Ｓ，Ｈｕａｎｇ　Ｊ　Ｌ，ｅｔ　ａｌ．Ａｃｃｕｍｕｌａｔｉｏｎ ｃｏｎｄｉｔｉｏｎｓ　ｏｆ　ｓｈａｌｅ　ｇａｓ　ｒｅｓｅｒｖｏｉｒｓ　ｉｎ　Ｓｉｌｕｒｉａｎ　ｏｆ　ｔｈｅ　Ｕｐｐｅｒ Ｙａｎｇｔｚｅ　ｒｅｇｉｏｎ．Ｎａｔｕｒａｌ　Ｇａｓ　Ｉｎｄｕｓｔｒｙ（ｉｎ　Ｃｈｉｎｅｓｅ），２００９， ２９（５）：４５－５０． ［１６］　张抗，谭云冬．世界页岩气资源潜力和开采现状及中国页岩 气发展前景．当代石油化工，２００９，１７（３）：９－１８． Ｚｈａｎｇ　Ｋ，Ｔａｎ　Ｙ　Ｄ．Ｔｈｅ　ｓｔａｔｕｓ　ｏｆ　ｗｏｒｌｄ　ｓｈａｌｅ　ｇａｓ　ｒｅｓｏｕｒｃｅｓ ｐｏｔｅｎｔｉａｌ　ａｎｄ　ｐｒｏｄｕｃｔｉｏｎ　ｓｔａｔｕｓ　ａｓ　ｗｅｌｌ　ａｓ　ｄｅｖｅｌｏｐｍｅｎｔ ｐｒｏｓｐｅｃｔ　ｏｆ　Ｃｈｉｎａ＇ｓ　ｓｈａｌｅ　ｇａｓ．Ｐｅｔｒｏｌｅｕｍ　＆ Ｐｅｔｒｏｃｈｅｍｉｃａｌ Ｔｏｄａｙ（ｉｎ　Ｃｈｉｎｅｓｅ），２００９，１７（３）：９－１８． ［１７］　张金川，聂海宽，徐波等．四川盆地页岩气成藏地质条件． 天然气工业，２００８，２８（２）：１５１－１５６． Ｚｈａｎｇ　Ｊ　Ｃ，Ｎｉｅ　Ｈ　Ｋ，Ｘｕ　Ｂ，ｅｔ　ａｌ．Ｇｅｏｌｏｇｉｃａｌ　ｃｏｎｄｉｔｉｏｎ　ｏｆ ｓｈａｌｅ　ｇａｓ　ａｃｃｕｍｕｌａｔｉｏｎ　ｉｎ　Ｓｉｃｈｕａｎ　Ｂａｓｉｎ． Ｎａｔｕｒａｌ　Ｇａｓ Ｉｎｄｕｓｔｒｙ（ｉｎ　Ｃｈｉｎｅｓｅ），２００８，２８（２）：１５１－１５６． ［１８］　黄籍中．四川盆地页岩气与煤层气勘探前景分析．岩性油气 藏，２００９，２１（２）：１１６－１２０． Ｈｕａｎｇ　Ｊ　Ｚ．Ｅｘｐｌｏｒａｔｉｏｎ　ｐｒｏｓｐｅｃｔ　ｏｆ　ｓｈａｌｅ　ｇａｓ　ａｎｄ　ｃｏａｌ－ｂｅｄ ｍｅｔｈａｎｅ　ｉｎ　Ｓｉｃｈｕａｎ　Ｂａｓｉｎ． Ｌｉｔｈｏｌｏｇｉｃ　Ｒｅｓｅｒｖｏｉｒｓ （ｉｎ Ｃｈｉｎｅｓｅ），２００９，２１（２）：１１６－１２０． ［１９］　李登华，李建忠，王社教等．页岩气藏形成条件分析．天然 气工业，２００９，２９（５）：２２－２６． Ｌｉ　Ｄ　Ｈ，Ｌｉ　Ｊ　Ｚ，Ｗａｎｇ　Ｓ　Ｊ，ｅｔ　ａｌ．Ａｎａｌｙｓｉｓ　ｏｆ　ｃｏｎｔｒｏｌｓ　ｏｎ　ｇａｓ ｓｈａｌｅ　ｒｅｓｅｒｖｏｉｒｓ．Ｎａｔｕｒａｌ　Ｇａｓ　Ｉｎｄｕｓｔｒｙ（ｉｎ　Ｃｈｉｎｅｓｅ），２００９， ２９（５）：２２－２６． ０３６１ 　４期 郝建飞，等：页岩气地球物理测井评价综述 ［２０］　张金川，金之钧，袁明生．页岩气成藏机理和分布．天然气 工业，２００８，２４（７）：１５－１８． Ｚｈａｎｇ　Ｊ　Ｃ，Ｊｉｎ　Ｚ　Ｊ，Ｙｕａｎ　Ｍ　Ｓ．Ｒｅｓｅｒｖｏｉｒｉｎｇ　ｍｅｃｈａｎｉｓｍ　ｏｆ ｓｈａｌｅ　ｇａｓ　ａｎｄ　ｉｔｓ　ｄｉｓｔｒｉｂｕｔｉｏｎ．Ｎａｔｕｒａｌ　Ｇａｓ　Ｉｎｄｕｓｔｒｙ （ｉｎ Ｃｈｉｎｅｓｅ），２００８，２４（７）：１５－１８． ［２１］　张金川，薛会，卞昌蓉等．中国非常规天然气勘探雏议．天 然气工业，２００６，２６（１２）：５３－５６． Ｚｈａｎｇ　Ｊ　Ｃ，Ｘｕｅ　Ｈ，Ｂｉａｎ　Ｃ　Ｒ，ｅｔ　ａｌ．Ｒｅｍａｒｋｓ　ｏｎ ｕｎｃｏｎｖｅｎｔｉｏｎａｌ　ｇａｓ　ｅｘｐｌｏｒａｔｉｏｎ　ｉｎ　Ｃｈｉｎａ． Ｎａｔｕｒａｌ　Ｇａｓ Ｉｎｄｕｓｔｒｙ（ｉｎ　Ｃｈｉｎｅｓｅ），２００６，２６（１２）：５３－５６． ［２２］　王湘玉．煤层气藏与页岩气藏．试采技术，２００９，３０（３）：２６－３２． Ｗａｎｇ　Ｘ　Ｙ．Ｃｏａｌｂｅｄ　ｍｅｔｈａｎｅ　ｒｅｓｅｒｖｏｉｒ　ａｎｄ　ｓｈａｌｅ　ｇａｓ ｒｅｓｅｒｖｏｉｒ．Ｗｅｌｌ　Ｔｅｓｔｉｎｇ　ａｎｄ　Ｐｒｏｄｕｃｔｉｏｎ　Ｔｅｃｈｎｏｌｏｇｙ （ｉｎ Ｃｈｉｎｅｓｅ），２００９，３０（３）：２６－３２． ［２３］　Ｈｉｌｌ　Ｄ　Ｇ，Ｌｏｍｂａｒｄｉ　Ｔ　Ｅ，Ｍａｒｔｉｎ　Ｊ　Ｐ．Ｆｒａｃｔｕｒｅｄ　ｓｈａｌｅ　ｇａｓ ｐｏｔｅｎｔｉａｌ　ｉｎ　Ｎｅｗ　Ｙｏｒｋ． Ｎｏｒｔｈｅａｓｔｅｒｎ　Ｇｅｏｌｏｇｙ　ａｎｄ Ｅｎｖｉｒｏｎｍｅｎｔａｌ　Ｓｃｉｅｎｃｅｓ，２００４，２６（１－２）：５７－７８． ［２４］　王宣龙，李厚裕，冯红霞．利用声波和自然伽马能谱分析泥 岩裂缝储层．测井技术，１９９６，２０（６）：４３２－４３５． Ｗａｎｇ　Ｘ　Ｌ，Ｌｉ　Ｈ　Ｙ，Ｆｅｎｇ　Ｈ　Ｘ．Ａｎａｌｙｓｉｓ　ｏｆ　ｆｒａｃｔｕｒｅ－ｔｙｐｅ ｓｈａｌｅ　ｒｅｓｅｒｖｏｉｒ　ｕｓｉｎｇ　ａｃｏｕｓｔｉｃ　ａｎｄ　ｎａｔｕｒａｌ　Ｇａｍｍａ　ｒａｙ ｓｐｅｃｔｒｏｍｅｔｒｙ　ｌｏｇｓ．Ｗｅｌｌ　Ｌｏｇｇｉｎｇ　Ｔｅｃｈｎｏｌｏｇｙ（ｉｎ　Ｃｈｉｎｅｓｅ）， １９９６，２０（６）：４３２－４３５． ［２５］　蔡开平，张本健．用声波测井资料对有利储层分布的初探． 天然气勘探与开发，２００５，２８（２）：６－９． Ｃａｉ　Ｋ　Ｐ，Ｚｈａｎｇ　Ｂ　Ｊ．Ｓｕｒｖｅｙ　ｄｉｓｔｒｉｂｕｔｉｏｎ　ｏｆ　ｆａｖｏｒａｂｌｅ ｐｒｏｓｐｅｃｔｓ　ｂｙ　ｓｏｎｉｃ　ｌｏｇｇｉｎｇ　ｄａｔａ－ｔａｋｉｎｇ　Ｓｈａｘｉｍｉａｏ　ｇｒｏｕｐ　ｏｆ Ｂａｉｍａｍｉａｏ　ｇａｓ　ｆｉｅｌｄ　ａｓ　ｅｘａｍｐｌｅ．Ｎａｔｕｒａｌ　Ｃｔａｓ　Ｅｘｐｌｏｒａｉｔｏｎ　＆ ＤｅｖＥｌｏｐｍｅｎｔ（ｉｎ　Ｃｈｉｎｅｓｅ），２００５，２８（２）：６－９． ［２６］　张碧星，王克协．各向异性介质地层中多极源声波测井的现 状与发展．地球物理学进展，１９９８，１３（２）：１－１４． Ｚｈａｎｇ　Ｂ　Ｘ，Ｗａｎｇ　Ｋ　Ｘ．Ｓｔａｔｕｓ　ｑｕｏ　ａｎｄ　ｄｅｖｅｌｏｐｍｅｎｔ　ｏｆ　ｔｈｅ ａｃｏｕｓｔｉｃ　ｍｕｌｔｉｐｏｌｅ　ｌｏｇｇｉｎｇ　ｉｎ　ａｎｉｓｏｔｒｏｐｉｃ　ｆｏｒｍａｔｉｏｎ．Ｐｒｏｇｒｅｓｓ ｉｎ　Ｇｅｏｐｈｙｓｉｃｓ（ｉｎ　Ｃｈｉｎｅｓｅ），１９９８，１３（２）：１－１４． ［２７］　Ｑｕｉｎｎ　Ｔ　Ｈ，Ｄｗｙｅｒ　Ｊ，Ｗｏｌｆｅ　Ｃ，ｅｔ　ａｌ．Ｆｏｒｍａｔｉｏｎ　Ｅｖａｌｕａｔｉｏｎ Ｌｏｇｇｉｎｇ　Ｗｈｉｌｅ　Ｄｒｉｌｌｉｎｇ（ＬＷＤ）ｉｎ　Ｕｎｃｏｎｖｅｎｔｉｏｎａｌ　Ｒｅｓｅｒｖｏｉｒｓ ｆｏｒ　Ｐｒｏｄｕｃｔｉｏｎ　Ｏｐｔｉｍｉｚａｔｉｏｎ．ＳＰＥ１１９２２７，２００８． ［２８］　Ｂｏｙｅｒ　Ｃ，Ｋｉｓｃｈｎｉｃｋ　Ｊ，Ｌｅｗｉｓ　Ｒ　Ｅ，ｅｔ　ａｌ．Ｐｒｏｄｕｃｉｎｇ　Ｇａｓ　ｆｒｏｍ Ｉｔｓ　Ｓｏｕｒｃｅ．Ｏｉｌｆｉｅｌｄ　Ｒｅｖｉｅｗ，２００６：３６－４９． ［２９］　Ｚｈａｏ　Ｈ，Ｇｉｖｅｎｓ　Ｎ　Ｂ，Ｃｕｒｔｉｓ　Ｂ．Ｔｈｅｒｍａｌ　ｍａｔｕｒｉｔｙ　ｏｆ　ｔｈｅ Ｂａｒｎｅｔｔ　Ｓｈａｌｅ　ｄｅｔｅｒｍｉｎｅｄ　ｆｒｏｍ　ｗｅｌｌ－ｌｏｇ　ａｎａｌｙｓｉｓ．ＡＡＰＧ Ｂｕｌｌｅｔｉｎ，２００７，９１（４）：５３５－５４９． ［３０］　Ｄｕｂｅ　Ｈ　Ｇ，Ｃｈｒｉｓｔｉａｎｓｅｎ　Ｇ　Ｅ，Ｆｒａｎｔｚ　Ｊ　Ｈ，ｅｔ　ａｌ．Ｔｈｅ　Ｌｅｗｉｓ Ｓｈａｌｅ，Ｓａｎ　Ｊｕａｎ　Ｂａｓｉｎ：Ｗｈａｔ　Ｗｅ　Ｋｎｏｗ　Ｎｏｗ．ＳＰＥ６３０９１， ２０００． ［３１］　Ｋａｌｅ　Ｓ　Ｖ，Ｒａｉ　Ｃ　Ｓ，Ｓｏｎｄｅｒｇｅｌｄ　Ｃ　Ｈ，ｅｔ　ａｌ．Ｐｅｔｒｏｐｈｙｓｉｃａｌ Ｃｈａｒａｃｔｅｒｉｚａｔｉｏｎ　ｏｆ　Ｂａｒｎｅｔｔ　Ｓｈａｌｅ．ＳＰＥ１３１７７０，２０１０． ［３２］　Ｐａｓｓｅｙ　Ｑ　Ｒ，Ｃｒｅａｎｅｙ　Ｓ，Ｋｕｌｌａ　Ｊ　Ｂ，ｅｔ　ａｌ．Ａ　ｐｒａｃｔｉｃａｌ　ｍｏｄｅｌ ｆｏｒ　ｏｒｇａｎｉｃ　ｒｉｃｈｎｅｓｓ　ｆｒｏｍ　ｐｏｒｏｓｉｔｙ　ａｎｄ　ｒｅｓｉｓｔｉｖｉｔｙ　ｌｏｇｓ．ＡＡＰＧ Ｂｕｌｌｅｔｉｎ，１９９０，７４（１２）：１７７７－１７９４． ［３３］　谭茂金，张松扬．页岩气储层地球物理测井研究进展．地球 物理学进展，２０１０，２５（６）：２０２４－２０３０． Ｔａｎ　Ｍ　Ｊ，Ｚｈａｎｇ　Ｓ　Ｙ．Ｐｒｏｇｒｅｓｓ　ｉｎ　ｇｅｏｐｈｙｓｉｃａｌ　ｌｏｇｇｉｎｇ ｔｅｃｈｎｏｌｏｇｙ　ｆｏｒ　ｓｈａｌｅ　ｇａｓ　ｒｅｓｅａｒｖｏｉｒｓ．Ｐｒｏｇｒｅｓｓ　ｉｎ　Ｇｅｏｐｈｙｓｉｃｓ （ｉｎ　Ｃｈｉｎｅｓｅ），２０１０，２５（６）：２０２４－２０３０． ［３４］　潘仁芳，伍媛，宋争．页岩气勘探的地球化学指标及测井分 析方法初探．中国石油勘探，２００９，１４（３）：６－９，２８． Ｐａｎ　Ｒ　Ｆ，Ｗｕ　Ｙ，Ｓｏｎｇ　Ｚ．Ｇｅｏｃｈｅｍｉｃａｌ　ｐａｒａｍｅｔｅｒｓ　ｆｏｒ　ｓｈａｌｅ ｇａｓ　ｅｘｐｌｏｒａｔｉｏｎ　ａｎｄ　ｂａｓｉｃ　ｍｅｔｈｏｄｓ　ｆｏｒ　ｗｅｌｌ　ｌｏｇｇｉｎｇ　ａｎａｌｙｓｉｓ． Ｃｈｉｎａ　Ｐｅｔｒｏｌｅｕｍ　Ｅｘｐｌｏｒａｔｉｏｎ（ｉｎ　Ｃｈｉｎｅｓｅ），２００９，１４（３）：６－ ９，２８． ［３５］　Ｒｏｓｅ　Ｄ　Ｊ　Ｋ，Ｂｕｓｔｉｎ　Ｒ　Ｍ．Ｃｈａｒａｃｔｅｒｉｚｉｎｇ　ｔｈｅ　ｓｈａｌｅ　ｇａｓ　ｒｅｓｏｕｒｃｅ ｐｏｔｅｎｔｉａｌ　ｏｆ　Ｄｅｖｏｎｉａｎ－Ｍｉｓｓｉｓｓｉｐｐｉａｎ　ｓｔｒａｔａ　ｉｎ　ｔｈｅ　Ｗｅｓｔｅｒｎ Ｃａｎａｄａ　ｓｅｄｉｍｅｎｔａｒｙ　ｂａｓｉｎ：Ａｐｐｌｉｃａｔｉｏｎ　ｏｆ　ａｎ　ｉｎｔｅｇｒａｔｅｄ ｆｏｒｍａｔｉｏｎ　ｅｖａｌｕａｔｉｏｎ．ＡＡＰＧ　Ｂｕｌｌｅｔｉｎ，２００８，９２（１）：８７－１２５． ［３６］　Ｃｏｏｐｅ　Ｄ　Ｆ，Ｑｕｉｎｎ　Ｔ　Ｈ，Ｍａｎｎｉｎｇ　Ｍ　Ｊ，ｅｔ　ａｌ．Ａ　Ｒｏｃｋ　Ｍｏｄｅｌ ｆｏｒ　Ｓｈａｌｅ　Ｇａｓ　ａｎｄ　Ｉｔｓ　Ａｐｐｌｉｃａｔｉｏｎ　Ｕｓｉｎｇ　Ｍａｇｎｅｔｉｃ　Ｒｅｓｏｎａｎｃｅ ａｎｄ　Ｃｏｎｖｅｎｔｉｏｎａｌ　ＬＷＤ　Ｌｏｇｓ