THE CHARACTERISTICS OF SHALE RESERVOIR FROM QIONGZHUSI FORMATION IN LOWER CAMBRIAN IN QUJING AREA, EASTERN YUNNAN PROVINCE, CHINA:A CASE STUDY OF WELL QUYE NO.1
-
摘要: 下寒武统筇竹寺组页岩是南方页岩气勘查开发重点关注的层位之一,但滇东地区下寒武统筇竹寺组黑色页岩储层研究程度较低。样品采自2014年新钻的曲页1井,对下寒武统筇竹寺组黑色泥页岩进行全岩和粘土矿物X衍射分析、TOC、岩石物性、比表面积等测试和扫描电镜实验,综合分析了滇东曲靖地区筇竹寺组泥页岩储层特征及孔隙类型。结果表明:滇东曲靖地区下寒武统筇竹寺组泥页岩矿物组成主要是石英(44.75%)、粘土矿物(25.47%)、长石(16.44%);发育粒间孔、粒内孔和有机孔3种类型孔隙;有机质孔隙发育,以中孔为主,具有一定的页岩气勘探开发潜力;页岩储层自上而下孔隙体积有逐渐变小的趋势,但在井深400 m~500 m发育高孔隙体积;泥页岩储集性能的重要指标—孔隙度主要受有机质含量及发育程度的制约。Abstract: Shale from Qiongzhusi formation in Lower Cambrian is one of the most important focuses in shale gas exploration and development in South China; however, the black shale reservoir from Qiongzhusi formation in Lower Cambrian in eastern Yunnan province stays in a low level of study. Take Well Quye No.1 drilled in 2014 as an example, on the black shales from Qiongzhusi formation in Lower Cambrian, X-ray diffraction analysis of the whole rock and clay minerals are made; tests and SEM experiments like TOC, physical property and specific surface area are conducted; the characteristics and types of shale reservoir are analyzed comprehensively. The results show that shales are mainly composed of quartz (44.75%), clay minerals(25.47%)and feldspar (16.44%), consisting of three major pore types, which are intergranular pore, intragranular pore and organic-matter pore. It is worth mentioning that organic-matter pores are abundant and of great importance, in which mesopores are the main force with a certain potential of shale gas exploration and development. The volumes of pores decrease with the depth, and with increasing TOC, micropore volumes relatively increase. The porosity is mainly restricted by the organic-matter content and the development degree.
-
Key words:
- pore type /
- pore volume /
- characteristics of reservoir /
- shale /
- Qiongzhusi formation /
- Well Quye NO.1
-
0. 前言
隧道是在山岭地区修建高速公路常见的一种地下建筑物, 连拱和分离式隧道是目前高速公路隧道的两种形式。连拱隧道由于具有线形顺畅、使用占地面积小等特点, 而在城市道路用地受限、山区地形复杂、道路辗线困难等情况下获得了广泛的应用。从力学上讲, 连拱隧道施工是不可逆的非线性演化过程, 他的最终状态不是唯一的, 而是与应力路径或应力历史有关。其在施工过程中的受力演化机制非常复杂, 施工过程中围岩的稳定性、初衬和二衬的受力情况、中墙的受力特征和稳定性等已成为工程技术人员关注的焦点。有限元法在模拟隧道的分步开挖和适时支护方面具有独特的优点, 随着计算机技术的不断进步, 得到了越来越广泛的采用。当前, 有关单洞隧道施工力学行为的三维有限元分析的文献并不鲜见[1~3], 研究结果对隧道的设计和施工起了很积极的作用。对双连拱隧道的研究也已深入, 文献[4~6]用二维有限元对连拱隧道的力学行为进行了较为全面的研究, 文献[7~10]则采用三维数值模拟的方法进行研究, 但成果往往集中于某一方面, 如地表沉降、中墙受力状态的演变、塑性区的发展等, 如文献[11]那样用三维有限元对双连拱隧道浅埋偏压洞身段的施工力学响应行为进行系统研究和全面分析的文献尚不多见。
本文所做的工作是:在考虑洞身浅埋偏压与洞口仰坡耦合作用的条件下, 采用功能强大的MARC有限元软件, 对某高速公路浅埋偏压连拱隧道出口段施工全过程进行三维数值模拟, 得出塑性区分布和发展、拱顶下沉、正应力与剪应力的集中和转移、中隔墙竖向应力随施工过程的变化规律, 研究结果对同类隧道的设计和施工具有参考和借鉴意义。
1. 依托工程概况及工程地质条件
隧道全长215米, 采用曲中墙连体隧道结构形式, 纵坡为2. 499 %。隧道区处于区域上塔前-赋春逆冲推覆构造带之塔前~赋春主断裂构造下盘, 主断裂距隧道洞身不足百米。受区域地质构造环境的影响, 隧道区裂隙发育, 岩体破碎。洞口均由千枚状砂质板岩强风化带组成, 岩土结构以蠕动状松散结构为主, 碎块状松散结构次之, 稳定性较差。
2. 三维有限元模型的建立
在X方向左右边界取至自然山脊和山谷线, 至隧道左右边墙距离约为5倍的单洞跨度, 两边施加X向的水平位移约束。在Y方向自由上边界取至地面, 下边界取至距开挖洞底4倍洞高处并施加竖向位移约束。在Z方向选取长度40m, 整个Z=40m平面取Z向水平位移约束。考虑实际施工情况, 对Z=0平面(洞口)而言, 坡体前缘线5m以下和洞周5m开外的交集部分取Z向水平位移约束。见图 1。
二衬、初衬和岩石都采用8节点三维实体等参单元, 共35700个; 锚杆用三维杆单元(图 2), 共3120个; 锚杆单元节点与实体单元节点完全重合, 总节点36393个。
假设锚杆是线弹性材料, 隧洞围岩和混凝土是弹塑性材料, 采用Drucker—Plager屈服准则。初始荷载均为岩土体的自重荷载。
另外, 隧洞锚杆是按梅花形布置, 径向间距0. 8米, 但是建模时为减少单元数量, 单元的尺寸大于实际的锚杆分布间距, 所以锚杆截面积为折算面积。围岩及其它材料的参数如表 1所示。
表 1 模型材料力学参数一览表Table 1. Mechanical parameters of model materials工况的模拟应尽可能与实际施工过程一致。但由于实际过程比较复杂, 而且某些施工过程可能同时进行, 因此只能依照实际的施工进度, 对施工步骤进行简化, 以尽量减少数值模拟计算的工作量, 同时又能够反映施工过程中所关心的问题。由于实际施工中一次性开挖长度经常可达30m, 考虑到最危险情况, 将所选取段按通槽开挖进行模拟。开挖采取“杀死”单元, 初期支护和二衬施做采取改变单元属性来“激活”[12]。
采用中导坑台阶法施工, 施工工序为:Lcase 0, 自重作用下初始状态的模拟; Lcase1, 开挖中导洞; Lcase2, 中导洞支护; Lcase3, 浇注中墙; Lcase4, 右洞上台阶开挖; Lcase5, 右洞上台阶初期支护; Lcase6, 开挖左洞上台阶; Lcase7, 左洞上台阶初期支护; Lcase8, 右洞下台阶开挖; Lcase9, 右洞下台阶初期支护; Lcase10, 浇筑右洞二衬及仰拱; Lcase11, 左洞下台阶开挖; Lcase12, 左洞下台阶初期支护; Lcase13, 浇筑左洞二衬及仰拱。右洞靠山谷一侧, 浅埋; 左洞靠山脊, 相对深埋。相对位置可参见图 1。
3. 动态数值模拟结果分析[13~16]
3.1 拱顶下沉
图 3、图 4中:坐标0处为洞口边界, 埋深8m;坐标40处为洞内边界, 埋深40m。位移系MARC输出后, 用Excel消掉自重作用下的初始位移(Lcase 0)后生成。图 5、图 6中: Lcase 0的拱顶下沉为自重作用下的初始值, 与施工影响无关。左右洞拱顶节点等间距(8m)分布, 两洞节点号由洞口至洞内依次增大, 左洞为92 —4395, 右洞为40 —3019。拱顶下沉小结:
1.从拱顶下沉与工况步关系曲线图(图 5、图 6)看, 隧道位移释放主要发生在隧道中导洞开挖(Lcase1), 左、右洞各自上台阶(左洞Lcase6、右洞Lcase4)开挖之时, 其他施工步对拱顶下沉影响不大。
2.埋深相对较大的山脊侧左洞拱顶下沉比山谷侧右洞拱顶下沉大, 前者最大4. 72cm, 后者最大3. 81cm, 均发生在进洞40m处(图 3、图 4)。各工况步拱顶下沉由洞口至洞内有逐渐增大的趋势, 这很可能是由于埋深不断增大引起的。在浅埋条件下, 拱顶下沉有随埋深增大而增大的趋势, 洞内稳定不容忽视。
3.实测拱顶下沉位移一般2~3cm, 计算值相对偏大, 但这正好反应了计算模型和参数选取的合理性。开挖是引起拱顶下沉的主要因素, 而监测工作往往是滞后于开挖的, 开挖初期释放掉的部分位移, 未能在监测数据里反映出来。
3.2 应力集中分析
a剪应力(图 7):
剪应力集中小结:
1.中导洞支护(Lcase2)的剪应力集中在墙底和拱腰上角点; 在左右洞上台阶开挖完毕后(Lcase7), 中墙的倒梯形墙基剪应力也非常集中; 中墙墙底的剪应力集中区随着左右洞二衬和仰拱的施作而减弱, 而墙顶的则增强。
2.施工完毕, 仰拱成了剪应力集中的区域, 比二衬的拱墙和边墙发育, 说明仰拱在支护中发挥了很重要的作用, 应及时施作。
b正应力(图 8):
正应力集中小结:
1.按MARC规定, 应力“拉正压负”。中导洞开挖完毕(Lcase1)隧道顶部即出现拉应力区, 由于岩土体抗拉能力弱, 因此, 顶部是易发生塌方的部位, 应超前或及时支护。
2.左右洞上台阶开挖初支后(Lcase7)洞顶仍然出现明显的拉应力区, 山脊侧洞顶的拉应力区仍然较大。说明山脊侧更容易坍方, 初支还不能保证隧道的稳定, 应及时施作二衬。
3.随着左洞下台阶的开挖, 隧道两侧边墙都出现了明显的压应力区, 这显然是由于隧道的偏压推挤作用引起的。
4.随着二衬及仰拱的及时完成, 隧道周边应力转移, 压应力集中区基本消失, 二衬及仰拱成为正应力集中区。
3.3 塑性区分析
塑性区小结(图 9):
1.中导洞开挖初期支护后(Lcase2)两侧和墙脚皆有塑性区发育, 说明中导洞虽然截面小, 但其稳定性仍是一个不可忽视的问题。
2.左右洞上台阶开挖后(Lcase7), 中导洞的支护仍未拆除, 所以其下部仍然出现较大塑性区, 这是由于中导洞支护拱向两侧传递推力的缘故。说明:该阶段不拆除中导洞支护, 可以发挥其支护作用, 有利于中隔墙稳定。
3.施工支护完毕阶段(Lcase11、Lcase13)山脊侧中墙墙踵和隧道边墙塑性区比山谷侧发育, 说明施工完毕后山脊侧所受压力比山谷侧大, 失稳最易在山脊侧发生。
3.4 中墙竖向应力分析(图 10、图 11、图 12)
说明:图 11和图 12中应力拉正压负。由于施工中采用改变单元属性来“激活”支护, Lcase 0的节点σy值实质是岩体在该点的自重应力, Lcase1和Lcase2工况步由于中导洞开挖(杀死单元), 节点σy值为0。而中墙在Lcase3才施作, 所以, 分析中墙节点应力的变化应从Lcase3开始。为避免边界条件影响, 研究对象取距洞口20m的横断面。节点15362和14816分别为该横断面上部的山脊侧节点和山谷侧节点; 节点14753和15425分别为该横断面下部山脊侧节点和山谷侧节点。
中墙竖向应力小结:
1.中导洞开挖施作中墙后(Lcase3)σy基本为负值, 说明中隔墙主要受压, 但山脊侧受压更明显, 说明隧道偏压。
2.从σy分布云图看(图 10), 在开挖形态对称条件下(Lcase7、Lcase11-Lcase13), 中墙竖对称条件下(Lcase7、Lcase11-Lcase13), 中墙竖向应力处于相对对称状态。随着右洞下台阶开挖(Lcase8), 中墙山谷侧成了受压集中区, 说明开挖导致山体有向外移动的趋势, 其推挤作用使中墙受到明显的偏压作用。右洞初衬和二衬的施作(Lcase10)并没使中墙偏压应力作出很大的调整, 说明非对称开挖是引起中墙偏压的关键因素。
3.从中墙横断面(距洞口20m)上部两侧σy随工况步变化曲线看(图 11), 山脊侧上部节点15362在山谷侧右洞上台阶开挖后(Lcase3-Lcase4)出现了明显的拉应力, 山谷侧节点14816在山脊侧左洞开挖(Lcase6-Lcase7)后, 也出现了明显拉应力, 15362节点几乎同时由受拉变成受压状态, 而此时中墙下部节点(15425、14753)仍处于受压状态。
4.从中墙横断面(距洞口20m)下部两侧σy随工况步变化曲线看(图 12), 中墙下部山脊侧节点14753在山谷侧右洞下台阶开挖(Lcase7-Lcase8)后出现最大竖向拉应力0. 903MPa, 山谷侧节点15425受最大压应力3. 031MPa。隧道施作完成后(Lcase13)中墙两侧所受压应力相差很小, 说明中墙倾覆失稳、拉裂破坏或压致破坏最可能在非对称开挖的单侧隧道的下台阶开挖期间发生。
4. 主要结论
通过对连拱隧道的施工进行动态模拟, 深入理解了隧道的施工力学响应行为, 归纳起来, 可得到以下一些主要结论:
1) 实测拱顶下沉位移一般2~3cm, 计算最大值4~5cm, 相对偏大, 但这正好反应了计算模型和参数选取的合理性。开挖后、支护前的位移变化最大, 而监测工作相对滞后, 初期开挖已经释放掉部分位移, 未能在监测数据里反映出来。浅埋条件下, 拱顶下沉随埋深增大而增大, 洞内稳定不容忽视。
2) 施工完毕, 二衬和仰拱成了应力集中区, 是主要的受力结构, 施工中应及时施作以使其发挥作用。
3) 开挖完毕, 抗拉能力弱的隧道顶部岩土体出现拉应力区, 是易发生塌方的部位, 应及时支护或采用大管棚等措施超前支护。左右洞上台阶开挖初支后洞顶仍然存在明显的拉应力区, 初支还不能保证隧道的稳定, 应及时施作二衬和仰拱。随着施工完毕, 两侧边墙附近压应力集中区转移消失, 二衬及仰拱呈现拉应力集中的状态, 隧道受力状态得到了良好改善。
4) 中墙墙底的剪应力集中区和施工完毕后中墙墙踵出现的塑性区都说明中墙的基础稳定问题不可忽视。塑性区在左右洞上台阶开挖后最发育, 隧道在该施工段最易失稳, 施工时应作好相应的监测工作。由于偏压作用, 施工支护完毕, 中墙墙踵和隧道山脊侧边墙塑性区比山谷侧发育。可采用中墙基础锚固、在山脊侧施作偏压衬砌等方法, 来加强隧道的稳定性。
5) 中墙竖向应力状态在施工过程中经过了复杂的调整, 在开挖形态非对称条件下, 隧道出现最不利的偏心受压(受拉)状态, 初衬和二衬的施作并未让此状态得到有效改善。而开挖形态对称施工步的中墙应力则处于相对对称状态。因此, 对称开挖是防止中墙倾覆失稳、拉裂破坏或压致破坏的最有效措施。
-
表 1 曲页1井下寒武统筇竹寺组泥页岩孔隙度和渗透率
Table 1. Porosity vs. permeability of shales in Qiongzhusi formation in lower Cambrian, Well Quye NO. 1
样品
编号埋深
/m孔隙度
/%渗透率
/×10-3μm2QY1-Y57 276.03 1.78 0.0193 QY1-Y56 281.26 1.56 0.0184 QY1-J4 295.03 2.18 0.0136 QY1-J5 301.94 3.01 0.0275 QY1-Y52 302.95 1.94 0.0196 QY1-Y50 313.75 1.35 0.0476 QY1-Y48 324.29 1.73 0.0149 QY1-J9 326.42 2.45 0.0361 QY1-Y47 328.54 1.69 0.0190 QY1-Y46 333.39 1.14 0.0188 QY1-J13 362.48 2.57 0.0275 QY1-Y41 371.07 1.64 0.0170 QY1-Y40 378.8 2.24 0.0153 QY1-J16 388.19 2.00 0.0116 QY1-J18 402.37 1.18 0.0030 QY1-J19 415.83 0.77 0.0076 QY1-J21 430.6 1.82 0.0110 QY1-J23 442.31 0.86 0.0103 QY1-Y33 442.75 1.94 0.0429 QY1-J25 451.23 1.18 0.0121 QY1-J27 457.06 2.10 0.0088 QY1-Y32 457.66 2.34 0.0399 QY1-J29 463.02 1.59 0.0067 QY1-J31 469.14 1.61 0.0117 QY1-Y30 471.00 1.72 0.0308 QY1-J32 474.12 1.72 0.0078 QY1-Y28 484.47 2.17 0.0256 QY1-J34 484.64 1.81 0.0077 QY1-Y27 492.56 0.91 0.0211 QY1-J36 496.78 1.43 0.0113 QY1-Y26 498.44 0.38 0.0179 QY1-Y25 501.06 0.48 0.0279 QY1-J38 513.22 1.07 0.0103 QY1-Y22 523.81 1.88 0.0255 QY1-Y20 535.75 0.81 0.0276 QY1-J40 535.91 1.89 0.0073 QY1-Y18 547.95 0.26 0.0104 QY1-Y17 552.47 1.47 0.0150 QY1-Y16 556.93 0.37 0.0183 QY1-Y14 567.17 1.65 0.0162 QY1-Y13 572.19 1.15 0.0095 表 2 曲页1井下寒武统筇竹寺组泥页岩比表面积测试结果
Table 2. Test results of specific surface area of shales in Qiongzhusi formation in lower Cambrian, Well Quye NO. 1
样品编号 QY1-J1 QY1-J5 QY1-J9 QY1-J11 QY1-J13 QY1-J15 QY1-J19 QY1-J21 QY1-J23 QY1-J25 QY1-J27 QY1-J32 QY1-J36 QY1-J38 QY1-J40 井深/m 245.75 301.94 326.42 338.15 362.48 382.89 415.83 430.6 442.31 451.23 457.06 474.12 496.78 513.22 535.91 表面积/(m2/g) 5.3130 3.9960 4.0330 3.2980 2.0870 2.8000 10.5310 6.1420 1.8000 4.3920 4.2500 4.7990 8.2760 1.9500 2.6860 孔隙体积/(ml/g) 0.0103 0.0076 0.0074 0.0068 0.0050 0.0060 0.0147 0.0116 0.0030 0.0072 0.0096 0.0081 0.0118 0.0042 0.0054 平均孔径/nm 3.8690 3.7985 3.6923 4.1234 4.8119 4.2959 2.7823 3.7632 3.3439 3.2928 4.5177 3.3714 2.8579 4.2600 4.0343 -
[1] 李新民. 我国首个超千亿立方米页岩气田诞生[OL]. http://www.mlr.gov.cn/xwdt/mtsy/xinhuanet/201407/t20140722_1324552.htm.LI Xinmin. The birth of China's first billion cubic meters of shale gas field[OL]. http://www.mlr.gov.cn/xwdt/jrxw/201407/t20140722_1324457.htm.2014, 7.22. (in Chinese) [2] 郭旭升.南方海相页岩气"二元富集"规律-四川盆地及周缘龙马溪组页岩气勘探实践认识[J].地质学报, 2014, 88(7):1209~1218. https://mall.cnki.net/qikan-DZXE201407001.htmlGUO Xusheng. Rules of two~factor enrichiment for marine shale gas in Sourthern China-Understanding from the Longmaxi formation shale gas in Sichuan basin and its surroungding area[J]. Geological Journal, 2014, 88(7):1209~1218. (in Chinese with English abstract) https://mall.cnki.net/qikan-DZXE201407001.html [3] 郭彤楼, 张汉荣.四川盆地焦石坝页岩气田形成与富集高产模式[J].石油勘探与开发, 2014, 41(1):28~36. doi: 10.11698/PED.2014.01.03GUO Tonglou, ZHANG Hanrong. Formation and enrichment mode of Jiaoshiba shale gas field, Sichuan basin[J]. Petroleum Exploration and Development, 2014, 41(1):28~36. (in Chinese with English abstract) doi: 10.11698/PED.2014.01.03 [4] 王玉芳, 翟刚毅, 王劲铸, 等.四川盆地及周缘龙马溪组页岩产气效果影响因素[J].地质力学学报, 2017, 23(4):540~547. http://journal.geomech.ac.cn/ch/reader/view_abstract.aspx?file_no=20170405&flag=1WANG Yufang, ZHAI Gangyi, WANG Jinzhu, et al. Factors influencing gas production effectiveness of Longmaxi formation shale in Sichun basin and adjacent areas[J]. Journal of Geomechanics, 2017, 23(4):540~547. (in Chinese with English abstract) http://journal.geomech.ac.cn/ch/reader/view_abstract.aspx?file_no=20170405&flag=1 [5] 程克明, 王世谦, 董大忠, 等.上扬子区下寒武统筇竹寺组页岩气成藏条件[J]..天然气工业, 2009, 29(5):40~44. http://www.oalib.com/paper/4982341CHENG Keming, WANG Shiqian, DONG Dazhong, et al. Accumulation conditions of shale gas reservoirs in the lower Cambrian Qiongzhusi formation, the upper Yangtze region[J]. Natural Gas Industry, 2009, 29(5):40~44. (in Chinese with English abstract) http://www.oalib.com/paper/4982341 [6] 张聪, 包书景, 石砥石, 等.云南曲靖地区下寒武统筇竹寺组页岩气成藏地质条件及有利区预测[J].海相油气地质, 2017, 22(1):69~74. http://www.oalib.com/paper/4982341ZHANG Cong, BAO Shujing, SHI Dishi, et al. Geological conditions for shale gas accumulation and favorable area prediction in Qujing of Eastern Yunnan[J]. Marine Origin Petroleum Geology, 2017, 22(1):69~74. (in Chinese with English abstract) http://www.oalib.com/paper/4982341 [7] 王丽波, 久凯, 曾维特, 等.上扬子黔北地区下寒武统海相黑色泥页岩特征及页岩气远景区评价[J].岩石学报, 2013, 29(9):3263~3278. http://d.wanfangdata.com.cn/Periodical/ysxb98201309024WANG Libo, JIU kai, ZENG Weiye, et al. Characteristics of lower Cambrian marine black shales and evaluation of shale gas prospective area in Qianbei area, upper Yangtze region[J]. Acta Petrologica Sinica, 2013, 29(9):3263~3278. (in Chinese with English abstract) http://d.wanfangdata.com.cn/Periodical/ysxb98201309024 [8] 周新桂, 刘和甫, 王宗秀, 等.我国叠合盆地油气调查战略思考[J].地质力学学报, 2011, 17(2):111~120. http://journal.geomech.ac.cn/ch/reader/view_abstract.aspx?file_no=20110201&flag=1ZHOU Xingui, LIU Hepu, WANG Zongxiu, et al. Thoughts of strategic Hydrocabron survey on superimposed basins in China[J]. Journal of Geomechanics, 2011, 17(2):111~120. (in Chinese with English abstract) http://journal.geomech.ac.cn/ch/reader/view_abstract.aspx?file_no=20110201&flag=1 [9] 杨勤生.滇东下古生界页岩气成藏层位及远景[J].云南地质, 2010, 29(1):1~6. http://www.cqvip.com/qk/95791X/201001/33201747.htmlYANG Qinsheng. The horizon & potentiality of lower Palaeozoic shale gas deposit in East Yunnan[J]. Yunnan Geology, 2010, 29(1):1~6. (in Chinese with English abstract) http://www.cqvip.com/qk/95791X/201001/33201747.html [10] 梁兴, 叶熙, 张介辉, 等.滇黔北下古生界海相页岩气藏赋存条件评价[J].海相油气地质, 2011, 16(4):11~21. http://www.cnki.com.cn/Article/CJFDTotal-HXYQ201104004.htmLIANG Xing, YE Xi, ZHANG Jiehui, et al. Evaluation of preservation conditions of lower Paleozoic marine shale gas reservoirs in the Northern part of Dianqianbei depression[J]. Marine Origin Petroleum Geology, 2011, 16(4):11~21.(in Chinese with English abstract) http://www.cnki.com.cn/Article/CJFDTotal-HXYQ201104004.htm [11] 梁兴, 叶熙, 张介辉, 等.滇黔北坳陷威信凹陷页岩气成藏条件分析与有利区优选[J].石油勘探与开发, 2011, 38(6):641~6519. http://www.wenkuxiazai.com/doc/c2c0399ddd88d0d233d46a31-2.htmlLIANG Xing, YE Xi, ZHANG Jiehui, et al. Reservoir forming conditions and favorable exploration zones of shale gas in the Weixin sag, Dianqianbei depression[J]. Petroleum Exploration and Development, 2011, 38(6):693~699. (in Chinese with English abstract) http://www.wenkuxiazai.com/doc/c2c0399ddd88d0d233d46a31-2.html [12] 云南省地质矿产局.云南省岩石地层[M].北京:中国地质大学出版社, 1996. 21~45.Bureau of Geology and Mineral Resources of Yunnan Province. Rock Strata in Yunnan Province[M]. Beijing:China University of Geosciences Press, 1996, 21~45. (in Chinese) [13] 云南省地质矿产局.云南省区域地质志[M].北京:地质出版社, 1990. 1~728.Bureau of Geology and Mineral Resources of Yunnan Province. Regional Geology of Yunnan Province[M]. Beijing:Geological Publishing House, 1990. 1~728. (in Chinese) [14] 云南省地质矿产局.云南省岩相古地理图集[M].昆明:云南科技出版社, 1995. 53~60.Bureau of Geology and Mineral Resources of Yunnan Province. Atlas of the Palaeo Geography of Yunnan Province[M]. Kunming:Yunnan Science and Technology Press, 1995. 53~60. (in Chinese) [15] Rouquerol J, Avnir D, Fairbridge C W, et al. Recommendations for the characterization of porous solids[J]. Pure Applied Chemistry, 1994, 66(8):1739~1758. http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.460.3717 [16] Ross D J K, Bustin R M. The importance of shale composition and pore structure upon gas storage potential of shale gas reservoirs[J]. Marine and Petroleum Geology, 2009, 26(6):916~927. doi: 10.1016/j.marpetgeo.2008.06.004 [17] 蒲泊伶, 董大忠, 耳闯, 等.川南地区龙马溪组页岩有利储层发育特征及其影响因素[J].天然气工业, 2013, 33(12):41~47. http://www.doc88.com/p-4562199214457.htmlPU Boling, DONG Dazhong, ER Chuang, et al. Favorable reservoir characteristics of the Longmaxi shale in the southern Sichuan basin and their influencing factors[J]. Natural Gas Industry, 2013, 33(12):41~47. (in Chinese with English abstract) http://www.doc88.com/p-4562199214457.html [18] 杨峰, 宁正福, 胡昌蓬, 等.页岩储层微观孔隙结构特征[J].石油学报, 2013, 34(2):301~311. doi: 10.7623/syxb201302012YANG Feng, NING Zhengfu, HU Changpeng, et al. Characterization of microscopic pore structures in shale reservoirs[J]. Acta Petrolei Sinica, 2013, 34(2):301~311. (in Chinese with English abstract) doi: 10.7623/syxb201302012 [19] 于炳松.页岩气储层的特殊性及其评价思路和内容[J].地学前缘, 2012, 19(3):252~258. http://d.wanfangdata.com.cn/Periodical/dxqy201203027YU Bingsong. Particularity of shale gas reservoir and its evaluation[J]. Earth Science Frontiers, 2012, 19(3):252~258. (in Chinese with English abstract) http://d.wanfangdata.com.cn/Periodical/dxqy201203027 [20] 于炳松.页岩气储层孔隙分类与表征[J].地学前缘, 2013, 20(4):211~220. https://mall.cnki.net/qikan-CXYY201715138.htmlYU Bingsong. Classification and characterization of gas shale pore system[J]. Earth Science Frontiers, 2013, 20(4):211~220. (in Chinese with English abstract) https://mall.cnki.net/qikan-CXYY201715138.html -