巢湖地区二叠系栖霞组和三叠系南陵湖组石灰岩生物标志物特征与生烃潜力
——以平顶山和马家山剖面为例
2015-11-23宋换新文志刚包建平
宋换新,文志刚,包建平
(1长江大学地球环境与水资源学院;2非常规油气湖北省协同创新中心)
巢湖地区二叠系栖霞组和三叠系南陵湖组石灰岩生物标志物特征与生烃潜力
——以平顶山和马家山剖面为例
宋换新1,2,文志刚1,包建平1
(1长江大学地球环境与水资源学院;2非常规油气湖北省协同创新中心)
对安徽巢湖地区平顶山和马家山剖面下二叠统栖霞组和下三叠统南陵湖组进行了地层沉积特征及暗色石灰岩生烃潜力的综合分析,并对这两个层位的样品进行了饱和烃生物标志物特征研究。结果表明,两个层位的石灰岩都属于有效烃源岩,均具有一定的生烃潜力。南陵湖组石灰岩饱和烃以高丰度的长链三环萜烷、重排藿烷、重排甾烷为特征,表明其形成于偏弱氧化的沉积环境。栖霞组沉积期水体能量相对较低,受陆源物质影响明显,栖霞组石灰岩萜烷系列为常规的分布模式,但其甾烷成熟度参数明显低于平衡值,认为高的热演化程度是造成栖霞组石灰岩20S/(20S+ 20R)-ααα-C29甾烷值和ββ/(ββ+αα)-C29甾烷值“倒转”的主要原因。
巢湖地区;南陵湖组;栖霞组;石灰岩;烃源岩;生物标志物;饱和烃;生烃潜力
安徽省中南部在二叠纪—三叠纪属于下扬子准地台,该时期沉积的海相碳酸盐岩地层在巢湖市北郊有两处发育典型、出露齐全的剖面——马家山剖面和平顶山剖面,其中的平顶山剖面是全球下三叠统印度阶—奥伦尼克阶界限层型首选剖面,即地学界俗称的“金钉子”(全球标准层型和点位)所在地。作为国际地层标准剖面,前人对于该区的沉积地层、古生物等方面的研究比较深入,取得了许多研究成果[1-3],但其他方面的研究涉及较少。本研究在野外剖面详细踏勘的基础上,选取巢湖地区马家山剖面和平顶山剖面进行系统的采样,对下三叠统南陵湖组和下二叠统栖霞组暗色海相碳酸盐岩(石灰岩)样品进行系统的地球化学分析,并在地层特征以及岩石生烃潜力综合分析的基础上,对这两个层位石灰岩的氯仿抽提物中饱和烃生物标志化合物特征进行对比研究,试图揭示其中包含的地质及地球化学信息,以期通过该研究能给相关地层沉积环境及油气勘探研究提供借鉴。
1 地层特征及生烃潜力评价
1.1 剖面地层特征
安徽省中南部地区的震旦系—三叠系为海相沉积,总厚度在14000m左右。该区地层沉积受江南断裂的控制,江南断裂以北属下扬子分区,以南属江南分区,巢湖地区属于下扬子地层分区[4]。本次研究,对平顶山和马家山两处野外剖面进行了详细观测,剖面上南陵湖组和栖霞组的具体沉积特征如下:
南陵湖组(T1n)马家山剖面,沉积厚度160m左右。岩性为灰色—深灰色石灰岩,夹薄层瘤状灰岩。上部地层白云质含量稍高,局部含白云质团块状灰岩,下部泥质含量稍高。生物化石以菊石为主,见少量有孔虫、牙形刺和鱼化石。南陵湖组在沉积早期为开阔海台地相,晚期水体变浅,盐度增加,为局限海台地相。
栖霞组(P1q)平顶山剖面,沉积厚度170m左右。岩性主要为灰色—灰黑色石灰岩、生物碎屑灰岩、燧石结核灰岩、硅质岩,夹有薄层的灰色—灰黑色页岩、碳质页岩及少量粉砂岩和煤线。生物以底栖的虫筳类、珊瑚、腕足类为主,植物化石丰富。研究区以局限海台地碳酸盐岩沉积为主,夹有滨海相、滨海沼泽相碎屑岩沉积。
采样剖面的岩性和古生物资料分析表明:研究区早二叠世栖霞期较早三叠世南陵湖期水体能量低,还原性强;栖霞组部分层段属于海陆过渡相沉积,受陆源物质的影响明显。
1.2 生烃潜力评价
通过野外剖面实测,在平顶山和马家山采石场新鲜剖面上进行了密集采样,其中马家山剖面上采集南陵湖组石灰岩样品46块,在平顶山剖面采集栖霞组石灰岩样品91块。对所采样品进行总有机碳(TOC)、岩石热解(Tmax)、干酪根碳同位素(δ13C)以及沥青质反射率(Rb)等相关测试(见表1)。结果表明,南陵湖组和栖霞组石灰岩的TOC平均值分别为0.32%和0.39%,考虑到较高的热演化程度,故认为这两个层位的石灰岩整体有机质丰度较高[5-6]。干酪根碳同位素和类型指数(TI)指标表明,两个层位的有机质类型都主要为Ⅰ型有机质。两个层位样品的H/C原子比平均值都相对较低,说明受到较高的热演化程度的影响。岩石热解烃最高峰温和沥青质反射率值表明,南陵湖组整体处于成熟热演化阶段,栖霞组整体处于高成熟热演化阶段。
表1 巢湖地区南陵湖组和栖霞组石灰岩样品生烃潜力评价参数表
考虑到研究区三叠系油气苗分布广泛,其中马家山剖面南陵湖组见多处以方解石脉晶洞原油和裂隙沥青形式产出的油苗[4,7],表明该层系石灰岩具有一定的生烃能力[8-9]。平顶山剖面栖霞组具有更高的残余TOC值,考虑到它具有相对更高的热演化程度,故推断其原始有机质丰度可能更高。整体认为:栖霞组和南陵湖组石灰岩都达到了有效烃源岩的标准,具有一定的生烃潜力[5-6]。
2 饱和烃生物标志物分布特征及环境意义
在烃源岩生烃潜力(表1)分析的基础上,选取在剖面上不同位置、TOC值较高的石灰岩样品(其中南陵湖组和栖霞组分别有5个和11个样品,其TOC均值为1.32%)进行了氯仿抽提、组分分离、饱和烃和芳烃色-质谱分析,以研究其生物标志物分布特征及其环境意义。
2.1 正构烷烃与类异戊二烯烷烃
研究区南陵湖组和栖霞组两个层位海相石灰岩的饱和烃中,正构烷烃整体分布面貌均为前峰型(图1),碳数分布范围一般为nC11~nC31,主峰碳为nC13~nC19,且以nC15~nC17为主(表2),无明显的奇偶优势;nC21-/nC22+值在南陵湖组5个样品的平均值为1.86,栖霞组11个样品的平均值为3.15(表2),显示了以低碳数饱和烃占优势的前峰型分布模式,符合研究区海相有机质的沉积以菌藻类为主要来源的特征[10]。在南陵湖组石灰岩样品饱和烃总离子流图(图1a)上,发现在nC25以后出现一个由环状甾、萜类等生物标志物组成的不可分辨的“鼓包”,这一现象在下扬子地区下古生界石灰岩中是很常见的,包建平等[11]认为这种现象可能与当时繁盛的浮游植物中富含环状化合物有关,即与特殊的有机质生源构成有关。
姥鲛烷(Pr)和植烷(Ph)是常用的表征古环境的生物标志物。一般认为,强还原、高含盐环境的沉积物,常常具有强烈的植烷优势,Pr/Ph值一般小于0.6。见表2中:南陵湖组上部的3个样品(样号MJ-78、MJ-67和MJ-65)的Pr/Ph均值为0.56,而下部的2个样品(样号MJ-b7和MJ-b5)的Pr/Ph值分别为1.18和1.39,这一显著变化反映了南陵湖组沉积晚期是一种水循环受阻、盐度偏高的局限台地相中的潮间—潮下、中—低能潟湖环境[12]。
图1 巢湖地区南陵湖组和栖霞组代表性石灰岩样品饱和烃总离子流图
表2 巢湖地区南陵湖组和栖霞组石灰岩正构烷烃与类异戊二烯烷烃参数统计表
Pr/nC17—Ph/nC18相关图(图2)表明,研究区南陵湖组和栖霞组石灰岩中均以海相来源有机质为主,同时也很好地体现了栖霞组石灰岩的热演化程度较高。Pr/nC17—Ph/nC18—Pr/Ph相关三角图(图3)同样体现了栖霞组石灰岩较高的成熟度,这3个指标的相似的变化规律在下扬子其他地区的海相地层中也有报道[11]。
图2 巢湖地区南陵湖组和栖霞组石灰岩Pr/nC17—Ph/nC18相关图
图3 巢湖地区南陵湖组和栖霞组石灰岩Pr/nC17—Ph/nC18—Pr/Ph相关三角图
2.2 三环萜烷
三环萜烷系列在地质体中非常常见。研究区南陵湖组和栖霞组石灰岩抽提物饱和烃中三环萜烷系列的分布样式不同。南陵湖组石灰岩中以长链三环萜烷系列分布完整(碳数分布范围一般为C19~C35,个别样品甚至到C37)、丰度较高为特征,以C24三环萜烷为主峰(图4a、4b),该分布样式比较少见。而栖霞组石灰岩抽提物中未见高丰度的长链三环萜烷分布,三萜烷碳数分布范围一般为C19~C29,以C23三环萜烷为主峰,属常见分布模式(图4c、4d)。两个组石灰岩中的四环萜烷均仅见C24四环萜烷。研究区南陵湖组和栖霞组石灰岩抽提物饱和烃中,三环萜烷和四环萜烷分布模式属于海相沉积物中的分布模式[10]。
图4 巢湖地区南陵湖组和栖霞组石灰岩饱和烃萜烷系列(m/z 191)分布图
研究发现,南陵湖组石灰岩中高丰度的长链三环萜烷系列、高丰度的重排藿烷系列与相对较低丰度的C24四环萜烷伴生(图4,图5),这一现象与其形成环境以及有机质来源有关。研究表明,三环萜烷系列,尤其是其中的长链三环萜烷,一般指示了低等生物藻类的贡献[13]。研究区南陵湖组石灰岩样品高丰度长链三环萜烷系列的存在表明,藻类等微生物可能是其主要有机质来源。而高丰度的C24四环萜烷一般认为与陆源有机质贡献及碳酸盐岩环境有关[10]。与南陵湖组石灰岩相比,栖霞组石灰岩中具有较高的C24四环萜烷/C26三环萜烷值以及较低的C24三环/C23三环萜烷值(图6),这表明南陵湖组有机质中低等藻类贡献较大,而栖霞组沉积有机质明显受到陆源有机质的影响。该认识与栖霞组地层中丰富的植物化石和数层煤线的特征相符,同时也表明栖霞组沉积期研究区频繁处于滨岸沼泽相等海陆过渡相的沉积环境,陆源供给相对充分。
图5 巢湖地区南陵湖组和栖霞组石灰岩Ts/(Ts+Tm)—C24四环/C26三环相关图
图6 巢湖地区南陵湖组和栖霞组石灰岩C24四环/C26三环—C24三环/C23三环相关图
2.3 藿烷系列
研究区南陵湖组和栖霞组石灰岩饱和烃中的藿烷系列也存在两种完全不同的分布模式:
南陵湖组石灰岩藿烷系列是以18α(H)-C27新藿烷(Ts)为基峰,包括18(H)-C29新藿烷(C29Ts)、17α(H)-21β(H)-C29降藿烷以及17α(H)-C29~32重排藿烷的特殊分布模式。正常的藿烷一般仅检测到C29降藿烷,且丰度较低,伽马蜡烷丰度极低。这种高丰度重排藿烷并伴随完整的长链三环萜烷系列的分布模式非常少见(图4a、4b)。
栖霞组石灰岩藿烷系列分布模式是以17α(H)-21β(H)-C30藿烷为基峰,包括成分为C27、C30~C35藿烷、C30莫烷、C29重排藿烷和重排莫烷以及Ts、C29Ts的常规分布模式。同时检测到一定丰度的伽马蜡烷(图4c、4d)。
对于南陵湖组石灰岩样品中这种高丰度重排藿烷的成因,现阶段还不十分清楚。Moldowan等[14]认为,Ts系列的前身物可能为里白烯和里白醇,而不是C35细菌藿烷多醇,即新藿烷和重排藿烷系列的前身物可能不同。朱杨明等[15]在研究川中龙场油田侏罗系湖相原油时,曾发现有类似的异常高丰度的以17α(H)重排藿烷为主峰的分布模式,并认为其与偏氧化性且富含黏土的沉积环境有密切关系。赵孟军等[16]对塔里木盆地原油的重排藿烷的研究表明,17α(H)重排藿烷可能与细菌类先母质有关,在氧化—亚氧化沉积环境和酸性介质、黏土矿物催化作用下易于形成。张文正等[17]对鄂尔多斯盆地延长组烃源岩中17α(H)重排藿烷的研究发现,亚氧化的环境、高含量的黏土矿物和有机质的分散状赋存形式更有利于17α(H)-15α甲基-C30重排藿烷的形成。刘海钰等[18]在研究松辽盆地梨树断陷下白垩统湖相烃源岩时,也发现了高丰度的17α(H)重排藿烷与长链的三环萜烷伴生的现象,认为可能与咸化环境下某些特殊藻类的贡献有关。本文结合研究区南陵湖组沉积环境分析认为:其高丰度重排藿烷与长链三环萜烷的分布模式可能与其偏弱氧化的沉积环境及该环境中咸水条件下某些特殊藻类的贡献有关。
2.4 甾烷系列
甾烷系列化合物主要来源于构成生物细胞壁的甾醇类物质[10]。南陵湖组和栖霞组石灰岩样品中,甾烷类生物标志物包括C21~C22孕甾烷系列、C27~C29规则甾烷系列和C27、C29的重排甾烷系列(图7)。就整体而言,南陵湖组样品的重排甾烷含量相对较高,5个样品C27重排/C27规则甾烷均值为0.81(栖霞组该参数的平均值为0.34)(表3),其中个别样品以C27重排甾烷为甾烷系列主峰(图7a)。研究表明,高含量的重排甾烷一般与黏土中酸性物质的催化作用有关,碳酸盐岩中高的重排甾烷含量较少见;也有研究者认为高的Eh值沉积物易形成重排甾烷,但其具体形成机理还不十分清楚[10]。南陵湖组石灰岩抽提物饱和烃中相对高丰度的重排藿烷系列,以及芳烃组分中相对高丰度的氧芴系列,均表明其沉积环境可能是一种高能量偏氧化性环境,由此推测南陵湖组石灰岩高含量的重排甾烷可能主要与偏氧化性的沉积环境有关。这种环境对于有机质的保存不利,其沉积有机质往往是部分被氧化改造的有机质[19]。
图7 巢湖地区南陵湖组和栖霞组石灰岩甾烷系列(m/z 217)分布图
表3 巢湖地区南陵湖组和栖霞组石灰岩甾烷、萜烷生物标志物参数统计表
南陵湖组和栖霞组石灰岩样品的C27~C29规则甾烷中,绝大部分以C27甾烷占优势,即整体上以水生藻类等低等生物来源为主。其中南陵湖组石灰岩中均以C27甾烷占优势,其C27/C29-20R-ααα-甾烷值均大于1;而栖霞组石灰岩的部分样品,该值小于1,即C29甾烷占优势(表3),表明陆源有机质对其有明显的影响。20R-ααα-C27—C28—C29甾烷三角图(图8)同样表明,这两个层位石灰岩中的有机质组成存在明显的差别。
图8 巢湖地区南陵湖组和栖霞组石灰岩20R-ααα-C27—C28—C29甾烷三角图
3 C29甾烷成熟度参数“倒转”及其成因
20S/(20S+20R)-ααα-C29甾烷和ββ/(ββ+αα)-C29甾烷是应用最广泛的生物标志物成熟度参数。研究区栖霞组高成熟石灰岩20S/(20S+20R)-ααα-C29甾烷和ββ/(ββ+αα)-C29甾烷的平均值(0.47和0.42)明显低于其平衡值(分别为0.52~0.55和0.67~0.71),也低于成熟阶段的南陵湖组石灰岩样品的值(图9,表3),栖霞组石灰岩样品这两个参数存在明显的“倒转”,这是一个值得研究的现象。
图9 巢湖地区南陵湖组和栖霞组石灰岩甾烷相关图
陈世加等[20]在研究四川盆地上三叠统香溪群、侏罗系沙溪庙组和大庆油田昌德地区的侏罗系含气砂岩储层沥青时发现,高成熟的储层沥青,其20S/(20S+20R)-ααα-C29甾烷和ββ/(ββ+αα)-C29甾烷的值明显低于平衡值,也低于该地区成熟度较低的源岩和储层沥青的值,存在“倒转”现象。Strachan等[21]在研究澳大利亚卡那封盆地侏罗系页岩时发现,20S/(20S+20R)-ααα-C29甾烷值在镜质体反射率(Ro)值小于1.12%时,是随成熟度增加而增加的,但在Ro值大于1.12%时,该值随成熟度增加反而减小。Peters等[22]通过页岩加水热解实验,也证实了C29甾烷成熟度参数在高成熟、过成熟阶段存在“倒转”的现象。
由上述分析表明,研究区栖霞组石灰岩整体处于高成熟阶段,而南陵湖组石灰岩整体处于成熟阶段。用这两个参数分别与Tmax参数作相关图(图10),发现它们均存在明显的负相关性,这表明高的热演化程度是造成研究区栖霞组石灰岩样品20S/(20S+20R)-ααα-C29甾烷和ββ/(ββ+αα)-C29甾烷参数“倒转”的主要原因。其成因可能是在高热演化阶段各异构体之间的热稳定性不同,例如高成熟阶段20R构型被优先破坏就可以导致20S/(20S+20R)-参数值的增大。同时,已有研究也发现,碳酸盐和(或)磷酸盐可能在高成熟阶段甾烷分子异构化方面扮演重要角色[22]。
图10 巢湖地区南陵湖组和栖霞组石灰岩甾烷与Tmax相关图
4 结论
(1)研究区下二叠统栖霞组沉积期较下三叠统南陵湖组沉积期水体能量低,栖霞组明显受陆源物质的影响。南陵湖组和栖霞组石灰岩TOC平均值分别为0.32%和0.39%;有机质类型都为Ⅰ型;南陵湖组整体处于成熟热演化阶段,栖霞组整体处于高成熟热演化阶段,两者都有一定的生烃潜力。
(2)南陵湖组石灰岩饱和烃以高丰度的长链三环萜烷和重排藿烷、重排甾烷及极低丰度的正常藿烷系列为特征,三环萜烷和四环萜烷系列烷分布模式比较特殊,而栖霞组石灰岩以C30藿烷为常规分布模式,表明南陵湖组可能形成于偏弱氧化的沉积环境,其沉积有机质可能与该环境中某些特殊藻类的贡献有关。
(3)栖霞组高成熟石灰岩的饱和烃20S/(20S+20R)-ααα-C29甾烷和ββ/(ββ+αα)-C29甾烷值存在明显的“倒转”现象。这两个指标与Tmax值呈明显的负相关关系,在高成熟阶段这两个参数随成熟度的增加而出现热稳定性的差异,可能是“倒转”发生的主要原因。
[1]童金南,赵来时,左景勋,等.安徽巢湖地区下三叠统综合层序[J].地球科学中国地质大学学报,2005,30(1):40-46.
[2]童金南,Zakharov Y D,吴顺宝.安徽巢湖地区早三叠世菊石序列[J].古生物学报,2004,43(2):192-204.
[3]殷鸿福,童金南.关于中国的海相三叠系建阶及下三叠统分阶界线[J].地球科学中国地质大学学报,2002,27(5):490-497.
[4]安徽省地质局区域地质调查队.安徽地层志:二叠系分册[M].合肥:安徽科学技术出版社,1983.
[5]张水昌,梁狄刚,张宝民,等.塔里木盆地海相油气的生成[M].北京:石油工业出版社,2004.
[6]陈安定.海相“有效烃源岩”定义及丰度下限问题讨论[J].石油勘探与开发,2005,32(2):23-25.
[7]宋换新,肖天君,文志刚,等.皖南巢县马家山地区油苗地球化学特征及油源[J].石油勘探与开发,2007,34(5):557-561.
[8]宋换新,文志刚,唐友军,等.南陵—无为地区下二叠统栖霞组烃源岩评价与标定[J].天然气地球科学,2005,16(3):334-337.
[9]曾艳涛,文志刚,宋换新.南陵—无为地区下三叠统海相烃源岩评价与标定[J].海相油气地质,2006,11(4):34-38.
[10]Peters K E,Walters C C,Moldowan J M.The Biomarker Guide[M].London:Cambridge University Press,2005.
[11]包建平,王铁冠,王金渝,等.下扬子地区海相中、古生界有机地球化学[M].重庆:重庆大学出版社,1996.
[12]侯读杰,冯子辉.油气地球化学[M].北京:石油工业出版社,2011.
[13]Mackenzie A S.Application of biomarkers in petroleum geochemistry[C]//Wehe D.Advances of petroleum geochemistry. London:Academic Press,1987.
[14]Moldowan J M,Fago F J,Carlson R M K,et al.Rearranged hopanes in sediments and petroleum[J].Geochimica et Cosmochimica Acta,1991,55:3333-3353.
[15]朱杨明,钟荣春,蔡勋育,等.川中侏罗系原油重排藿烷类化合物的组成及成因探讨[J].地球化学,2007,36(3):253-260.
[16]赵孟军,张水昌.17α(H)-重排藿烷在塔里木盆地中的指相意义[J].石油勘探与开发,2001,28(1):36-38.
[17]张文正,杨华,侯林慧,等.鄂尔多斯盆地延长组不同烃源岩17α(H)-重排藿烷的分布及其地质意义[J].中国科学:D辑,2009,39(10):1438-1445.
[18]刘海钰,张敏.梨树断陷下白垩统湖相烃源岩异常高丰度17α(H)-重排藿烷的成因探讨[J].石油天然气学报,2013,35(7):29-33.
[19]张水昌,张宝民,边立曾,等.中国海相烃源岩发育控制因素[J].地学前缘,2005,12(3):39-48.
[20]陈世加,王廷栋,黄清德,等.C29甾烷成熟度参数指标“倒转”及其地质意义[J].天然气地球科学,1997,8(1):28-30.
[21]Strachan M G,Alexander R,Subroto E A,etc.Constraints upon the use of 24-ethylcholestane diasteromer ratios as indicator of the maturity of petroleum[J].Org Geochem,1989,14(4):423-432.
[22]Peters K E,Moldowan J M,Sundararaman P.Effects of hydrous pyrolysis on biomarker thermal maturity parameters:Montereyphosphatic and siliceous members[J].Org Geochem,1990,15(3):249-265.
编辑:赵国宪
Characteristics of Biomarkers and Hydrocarbon Potential in Lower Permian Qixia and Lower Triassic Nanlinghu Limestone:Cases from Pingdingshan and Majiashan Outcrops in Chaohu,Anhui
Song Huanxin,Wen Zhigang,Bao Jianping
It is known that the classical and well-outcropping Lower Permian Qixia limestone and Lower Triassic Nanlinghu limestone develop at Pingdingshan and Majiashan outcrops in Chaohu,Anhui province,between which the Pingdingshan outcrop is the first-choice one of the Lower Triassic Induan-Olenekian boundary of GSSP in China. The outcrop samples were adopted from these outcrops for analyzing biomarkers in saturated hydrocarbon.The analysis result indicates that both the dark Qixia and Nanlinghu limestones are of effective source rock with a certain hydrocarbon generation potential.The Nanlinghu limestone is characterized by high abundance of long-chain tricyclic terpane,rearranged hopane and diasterane in saturation hydrocarbon,which implies an environment of weak-prone oxidization during sedimentation.Owing to the influence of low-energy terrigenous sediments,the Qixia limestone is characteristic of terpane series in a normal distribution patter,while the sterane shows an obviously lower maturity parameter than the balance value of it.It is suggested the reason why the Qixia limestone presents high thermal evolution just is the inversion of 20S/(20S+20R)-ααα-C29sterane and ββ/(ββ+αα)-C29sterane values.
Lower Permian;Lower Triassic;Limestone,Source rock;Biomarker;Saturated hydrocarbon;Southern Anhui
TE135+.2
A
10.3969/j.issn.1672-9854.2015.02.004
1672-9854(2015)-02-0021-08
2014-04-18;改回日期:2014-06-29
本文受国家科技重大专项课题“中上扬子海相碳酸盐岩成藏条件与有利勘探区带评价研究”(编号2011ZX05004-001-004)资助
宋换新:1980年生,讲师。主要从事油气地球化学及油气成藏方面的研究工作。通讯地址:430100湖北省武汉市蔡甸区大学路特1号;E-mail:shx@yangtzeu.edu.cn
Song Huanxin:Lecturer,Doctoral candidate.Add:Hubei Cooperative Innovation Center of Unconventional Oil and Gas,Yangtze University,Te-1 Daxue Rd.,Wuhan,Hubei,430100,China