论文部分内容阅读
摘 要：综合运用沉积古地质分析法、印模法、残厚法等，对鄂尔多斯盆地苏里格东南区前石炭纪岩溶微地貌特征进行研究，并结合实际生产数据对其气水分布特征进行总结，发现研究区奥陶纪末岩溶古地貌形态表现为“西高东低、北高南低”的分布格局，发育2类二级古地貌单元和9类三级古地貌單元，其中岩溶斜坡为主体地貌单元，主要发育上斜坡、下斜坡、一级残丘和二级残丘，局部发育古洼地和古沟槽;研究区马家沟组上组合气水分布主要表现为西北部、东部及东南部地区富气，中西部地区富水的特征。结果表明：综合以本溪组顶部8#煤层为上基准面的石炭系印膜厚度、以K3标志层为下基准面的残余厚度及其各自与平均厚度间的相对差值，可对研究区岩溶微地貌进行有效刻画;古残丘及上斜坡地区是有利天然气富集高产区;水层、气水同产层因受燕山期地层翘倾影响，多分布于研究区中西部古岩溶上斜坡地区;古洼地、古沟槽溶蚀淋滤作用较弱，虽不利于天然气富集，但仍有少量低产气井分布。
关键词：苏里格东南部;古地貌恢复;印模法;残厚法;气水分布特征
中图分类号：TE 122.2 文献标志码：A
DOI：10.13800/j.cnki.xakjdxxb.2019.0316 文章编号：1672-9315（2019）03-0498-09
Abstract：By the comprehensive use of sedimentary paleogeological analysis，impression method，and residual thickness method，the karst micro-geomorphology of the Pre-Carboniferous area in southeast Sulige，Ordos Basin was studied.Combining with actual production data，the characteristics of gas-water distribution was summarized.The study found that the karst paleogeomorphology at the end of Ordovician in the study area was represented by the distribution pattern of “high in the west and south，low in the east and north”，and developed two types of second-order paleogeomorphic units and nine types of third-order paleogeomorphic units.The karst slope is the main physiognomy unit，which is mainly develop on the upper slope，lower slope，primary monadnock and secondary monadnock.Moreover the ancient depression and groove are developed in the local area.The distribution of gas and water in the Majiagou-formation of study area is mainly characterized by the rich gas in the northwest，east and southeast，and the rich water in the central and western regions.Results show that：The karst microtopography in the study area can effectively depicted，by combing the carboniferous film thickness of the 8# coal seam at the top of Benxi Formation and the residual thickness of the lower datum of K3 marker.The monadnock and upper slope area are favorable natural gas enrichment areas.Water is affected by the tilting of the Yanshan Strata，and is mostly distributed in the central and western regions of the upper karst slopes.The anticorrosive and leaching effects of ancient depressions and ancient trenches are weak，but there are still a few stripper well distributed.Key words：southeastern part of Sulige Area;impression method;residual thickness method;characteristics of gas and water distribution0 引 言据统计20%～30%的大气田发育于古风化壳岩溶储层之中[1]，20世纪70年代中后期至今，与古岩溶相关的古地貌研究逐渐受到国内外学者重视，并得以广泛应用。古地貌恢复是对历史时期盆地面貌变化进行分析的一种重要手段，是以构造变形作用为主，风化剥蚀、充填沉积、差异压实等综合作用的结果[2]。 目前常用的古地貌刻画方法有：地球物理法、印模法、残余厚度法、沉积学法、层序地层法及层拉平法、回剥法、构造趋势面转化法、计算机模拟法等[2-11]。苏中堂、王建民等运用地球物理法，通过提升测井、地震属性，建立古地貌的识别标志模型，最后运用井震等地球物理资料对模型进行验证、约束，直至接近最理想的程度[2-4]。该方法直观准确，但刻画精度较低。夏日元、魏新善等运用印模法及残厚法恢复古地貌[5-8]。印模法目前应用最为广泛，其原理是利用上覆地层与古地貌之间的“镜像”关系，通过上覆地层的厚度来半定量恢复古地貌形态。该方法简单易行，但基准面与剥蚀面之间地层的真实厚度恢复难度较大。残余厚度法与印模法原理相似，区别在于需在剥蚀开始之前所沉积的地层中，选取标志层作为基准面，并将其拉平，该基准面之上残余地层厚度即可代表当时的古地貌形态。该方法直观且刻画精细，但却未考虑构造及剥蚀差异对沉积厚度的影响。赵俊兴等运用沉积学法恢复古地貌，即借助沉积岩相古地理特征，从而定性地恢复当时的古地貌特征[9]。沉积学法恢复古地貌简单、易于操作，但影响因素复杂、工作量大，且精度较低。宋国奇、赵俊兴等运用层序地层法，将层序地层学引入古地貌恢复之中，在建立上覆地层等时层序格架的基础上，将其拉平，运用其与风化壳等势面之间的地层来反应古地貌形态[10-11]。该方法优点是上覆基准面相对等时性强，缺点是2个等时界面的选择较为困难。不同的学者虽然采用不同的方法对不同地区的古地貌进行恢复研究，但各种分析方法均各有利弊。此外，就鄂爾多斯盆地而言，盆地尺度、油气田尺度的岩溶古地貌研究已趋于成熟[12-15]，但井区尺度的沉积期古地貌研究仍处于探索之中[12]。因此在前人研究的基础上，结合研究区实际，文中首次采用综合法对苏里格东南部地区三级岩溶古地貌进行精细划分与解释，即：井震结合初步确定地层下古生界地层格架;统计奥陶系顶部地层出露情况，了解古构造格局及各地区剥蚀程度;合理选择基准面，采用印模法及残厚法初步恢复古地貌;最终依照每种方法的优弊进行取舍，多角度综合划定研究区三级古地貌划分标准，恢复岩溶古地貌形态，并分析不同地貌形态与气水分布间的关系。
1 地质背景鄂尔多斯盆地受加里东运动影响，于奥陶纪盆地中央呈现“L”形隆起构造雏形[16]。早奥陶世，盆地在经历短暂抬升之后，开始接受马家沟组碳酸盐岩沉积;随着加里东运动逐渐加剧，华北地台整体抬升，中奥陶及早石炭地层被广泛剥蚀，发育风化壳岩溶型储集体，与上覆的晚石炭统海陆交互源岩沉积层呈区域不整合接触[7]。后经晚三叠世至白垩纪内陆盆地发展，以及新生代周边断陷盆地发育形成了现今南北隆升、东抬西冲、盆内斜坡西倾与天环向斜相连的构造面貌[14-17]。苏里格东南部位于鄂尔多斯盆地陕北斜坡北部（图1），地区面积790 km2，完钻井468口，其中探井19口，在平缓西倾的单斜背景上发育一系列北东-北东东向复式鼻状构造[18-22]。奥陶纪末，研究区位于中央古隆起东侧，差异性抬升剥蚀致使该区风化壳储层发育，表现出较大的勘探开发潜力。
2 常规法岩溶古地貌恢复2.1 前石炭古地质特征在对研究区217口钻至奥陶系单井进行小层对比与划分的基础上，统计每口钻井所揭示的奥陶系顶面出露层位（表1），绘制研究区前石炭系古地质图（图2），宏观分析石炭系沉积前的古地貌特征。
2.2 “印模”地层分布特征及古地貌指示在选取标志层参考面时，按照全区稳定分布等时界面、与风化壳面距离较近、岩性界面易于识别3大原则[9，13]，选取本溪组顶部8#煤层为基准面（图3），采用填平补齐的沉积补偿原理[23]，对地层抬升剥蚀后古地貌轮廓进行刻画。本溪组由顶部煤层与底部大套泥岩夹砂层组成，且全区广泛分布。考虑压实减孔对砂泥岩及煤系地层厚度影响较大[24]，因此先对本溪组厚度进行压实校正，具体如下
2.2.1 求取压实系数[25]
2.2.2 求取去压实后目标层顶底深度[24]
2.2.3 求取去压实后目标层骨架高度及压实率
式中 Z为深度，m;为Z处对应孔隙度，%;0为地表孔隙度（表2），%;C为压实系数;D1，D2为分别表示目标地层当前的顶深和底深，m;D1′，D2′分别表示压实校正后目标层定顶深度，m;D′为校正后目标层厚度，m;D为当前目标层厚度，m;Φ为压实率，%.经计算，该区本溪组砂岩平均压实率为28.7%，泥岩及煤层压实率为63%左右。在对本溪组地层进行去压实校正之后，得到本溪组原始沉积厚度等值线图，进而表征古地貌形态（图4）。可知，其整体反映了奥陶末期苏里格东南部西高东低、北高南低的地貌形态，与古地质图分析结果较为一致。西北部及东北部地层厚度相当，多分布于10～17 m之间，向南地层厚度逐渐增大。南部地区表现出沟脊相间的特征，桃2-36-2井—苏93井、桃2-30-14井—桃47井、陕50井—桃2-30-25井均呈现为北西-南东向展布特征，地层最厚处可达26.2 m.西部地区平均地层厚度高于东部地区，其中桃2-16-16井区、桃2-16-20井区、桃2-1-8井区厚度最大，均大于30 m，对应在本溪组沉积时表现为低洼特征。
2.3 “残厚”地层分布特征及古地貌指示与印模法类似，残厚法恢复古地貌亦需根据已述3个条件合理选择基准面。参考前石炭古地质
在不同程度剥蚀现象，因此选取K3标志层（马五14底部凝灰岩层）作为基准面，将其底拉平，该标志层之上的残余地层厚度代表了古地貌展布形态特征[21-22]（图5，图6）。
可知，古地貌整体西北偏高，起伏强烈;东南、西南地区仍表现出近北西-南东向沟脊相间形态，与印模法分析结果较为一致。在加里东期强烈构造作用的影响下，残厚小于40 m的区域主要位于研究区西部、北部地区，厚度变化幅度大且连续性差，桃2-21-8井处残余地层厚度最小，为24.6 m.由此说明，西部及北部地区遭受风化及剥蚀作用强烈，属低地貌单元。大于60 m的地层分布于研究区内中部地区，面积小且起伏程度小，向西逐渐歼灭，其中桃2-12-13C1井、桃35井、桃2-16-19井、桃2-30-8井K3标志层之上残厚度均大于65 m，并且均有马六地层残余，属高地貌单元。40～60 m地层分布最为广泛，占全区面积3/4，其间厚度变化幅度较小，说明其地势起伏较为平缓，遭受剥蚀程度较小。3 综合法古地貌单元划分及其特征运用沉积学古地质分析，虽能较好确定古地貌形态，但其精度不够[2，4]。印模法容易操作，并能较为精细的恢复古地貌。残厚法对古地貌的刻画直观，也较为精细，但基准面的选择、压实校正的准确度等都会影响最终结果[4，9]。上述分析结果也体现出相同结论，各种方法之间虽结果相似，能相互印证，但局部结果仍有差异。因此，文中在古地貌恢复技术的优缺点分析的基础上，以鄂尔多斯盆地沉积背景为根本，通过去压实校正，将印模法作为重点，综合古地质分析以及残厚分析结果，将研究区奥陶系岩溶古地貌进一步精细划分为2个二级单元和9类三级单元（图7），具体划分指标见表3. 3.1 岩溶高地苏里格东南部地区古岩溶台地面积较小，主要分布于北部桃2-3-11井-桃2-7-14井区、桃2-3-25井-桃2-4-27井区，奥陶系地层剥蚀强烈，顶部主要出露马五14及马五12地层，去压实校区后石炭系厚度一般小于15 m.
3.2 岩溶斜坡岩溶斜坡是该地区主体古地貌单元，去压实后石炭系厚度为15～50 m，K3标志层之上残余厚度为30～55 m.以石炭系及残余地层相对厚度为标准，岩溶斜坡又可进一步划分为6个三级单元。
残丘处石炭系厚度多小于30 m，是古岩溶斜坡上呈孤立丘状分布的古高地，分布范围较小，主要位于桃2-5-8井区，桃2-16-18井区，桃2-8-24-陕166井区，在加里东期所经历的风化剥蚀时间较长。上斜坡、下斜坡为是残丘及洼地间的过渡地带，沿残丘及洼地周缘呈环带状分布，奥陶系顶部出露层系以马五21为主。上斜坡处坡降较大，垂向渗流与横向岩溶均发育，物质交换迅速。
下斜坡处坡降较小，降水滞留，随所经历的成岩作用较为充分，但介质流动性差，不利于溶蚀物质被迅速带走。
洼地是岩溶斜坡内的负相地形区，残余地层厚度多小于10 m，面积较小，多位于沟槽水系发育末端。由于该区域地势低、水体深，因此水体流动缓慢，化学沉淀作用强，不利于储集体的形成与保存。古沟槽位于上斜坡边缘，为条带状展布。研究区内共3条古沟槽经过，东南部2条沟槽为北西-南东向展布，西部地区沟槽为北东-南西走向长达15 km，对下伏地层下侵作用强烈，致使该处马五1+2地层几乎被剥蚀殆尽，奥陶系顶部出露马五41至马五23.
4 岩溶古地貌与气水分布关系
鄂尔多斯盆地奥陶系马家沟组顶部的岩溶古地貌控制着古岩溶作用类型及强度，进而改变了风化壳古岩溶储层的岩性及物性条件[5]。但受地层叠覆顺序及风化作用强度对储层影响的强度所限，前石炭纪古地貌对奥陶系上组合地层影响较大[8，28]，因此，文中重点讨论研究区前石炭纪末古地貌与奥陶系马家沟组上组合气水分布间的关系。苏里格东南部地区气水分布并不均匀，产气井多分布于西北、中东及南部地区，产水井较少且零星分布，气水同产井集中分布于中部地区（图7）。马五1+2亚段北部及东南部产气量高，其中桃2-9-24井试气无阻流量高于100×104 m3，日产气量高达25×104 m3;马五14小层西北及东部地区产气量相对较高，其中桃2-6-1井试气无阻流量为69.8×104 m3，日产气量高达15×104 m3.其产水量高于马五1+2亚段，例如桃38井、桃2-30-14井等日产水量均可达9 m3以上。
叠合古地貌与上组合气水分布（图7），并结合前人关于岩溶储层特征研究成果可知，研究区前石炭纪古地貌形态对奥陶系上组合气水差异性分布影响较大。古岩溶高地在研究区内面积较小，垂向渗流为主，多发育纵向溶洞且非均一性明显[28]，不利于形成有益于天然气聚集成藏的层状岩溶带，仅局部地势相对较高的残丘地区尚有天然气分布。古岩溶斜坡是古地貌主体单元，其中残丘风化壳岩石类型以白云岩为主，主要表现为淡水淋滤溶蚀作用[23-31]，具有良好的储集性能;一级残丘、二级残丘各自所包含的气井占比40.0%，26.7%;36.9%的差气层井位于二级残丘之上，其次为一级残丘含差气层井26.0%.气水同产井多位于二级残丘及上斜坡单元，分别为31.8%，27.3%.产水井则集中分布于上斜坡、二级残丘及下斜坡古地貌单元之内，占比20.8%，45.8%，12.5%.此外，沟槽两侧分布有3.1%的差气层井。古洼地为岩溶水的汇集区[28]，岩溶作用相对较弱，所含气层、水层均较少。综上，奥陶纪末期苏里格东南部地区主要呈现西高东低、北高南低的地貌形态，古岩溶高地分布局限，古岩溶斜坡广泛展布，其间三级地貌单元高低相间分布。至晚三叠世、中-晚侏罗世，上古煤系烃源岩进入生排烃高峰，天然气下穿奥陶系顶部铝土岩缺失部位发生垂向运移[22]，在奥陶系风化壳岩溶储层裂缝及孔隙中向下充注并发生侧向运移，在当时地势相对较高的一级残丘、二级残丘、上斜坡等部位聚集成藏。至燕山运动中期，整个盆地由东倾单斜转为西倾单斜。最初在西部古残丘及上斜坡圈闭中聚集的天然气，受侧向遮挡而被保留;其余低部位聚集的天然气在自身浮力及水动力作用下，向东部构造高部位发生大面积后期运移与调整，驱替高部位地层水向西部底势区调整。最终出现现今中东部及前石炭古地貌呈岩溶高地貌单元的西北、南部中富气，中西部地区富水的气水分布格局。
5 结 论1）采用沉积古地质分析，选取本溪组顶部8#煤层为基准面的印模法，以K3标志层（马五14底部凝灰岩层）作为基准面的残余厚度法，均可对研究区前石炭纪古地貌进行有效恢复。各方法分析结果类似，可相互印证，但单一的沉积学分析、印模法或残厚法仍存在局限性。沉积古地质分析随操作简单，但精度较差;印模法虽能较好的反应古地貌格局，但去压实校正工作繁琐，对微地貌刻画反映较差;残厚法对微地貌刻画较为精细，但分析结果与整体格局稍有出入。
2）综合法恢复古地貌形态，将印模法与残厚法二者有机结合，并参考印模地层厚度、残余地层厚度与平均厚度值间的相对变化情况，结合研究区实际将前石炭古地貌精细刻画为2类二级单元和9类三级地貌单元。苏里格东南部地区奥陶系岩溶古地貌形态为“西高东低、北高南低”的分布格局，以岩溶斜坡地貌单元为主，主要发育上斜坡、下斜坡、一级残丘和二级残丘，局部发育古洼地和古沟槽。3）苏里格东南部地区奥陶系上组合气水主要表现为西北、东部及东南部富气，中西部富水的特征。水層、气水同层因受燕山期地层翘转影响，多分布于研究区中西部古岩溶上斜坡地区。古洼地、古沟槽溶蚀淋滤作用较弱，不利于天然气富集，但仍有少量低产气井分布其中。岩溶斜坡古残丘及上斜坡地区则是天然气高产井的目标区。
参考文献（References）：
[1]赵永刚，王东旭，冯强汉，等.油气田古地貌恢复方法研究进展[J].地球科学与环境报，2017，39（4）：516-529.ZHAO Yong-gang，WANG Dong-xu，FENG Qiang-han，et al.Review on palaeomorphologic reconstruction methods in oil and gas fields[J].Journal of Earth Science and Environment，2013，39（4）：516-529. [2]苏中堂.鄂尔多斯盆地塔巴庙地区奥陶系马家沟组沉积相与古岩溶储层特征[D].成都：成都理工大学，2008.SU Zhong-tang.Sedimentary facies and character of palaeokarst cust reservoir in Majiagou Formation in Ordovician，Tabamiao Area，Ordos[D].Chengdu：Chengdu University of Technology，2008.
[3]王建民，王佳媛，沙建怀，等.鄂尔多斯盆地东部奥陶系风化壳岩溶古地貌特征及综合地质模型[J].吉林大学学报（地球科学版），2014，44（2）：409-418.WANG Jian-min，WANG Jia-yuan，SHA Jian-huai，et al.Karst paleogeomorphology and comprehensive geological model of the Ordovician weathering crust in the eastern Ordos Basin[J].Journal of Jilin University （Earth Science Edition），2014，44（2）：409-418.
[4]王建民.鄂尔多斯盆地东部奥陶系风化壳岩溶古地貌与储层特征[D].西安：西北大学，2011.WANG Jian-min.Karst paleogeomorphology and reservoir characteristics of Ordovician weathered crust in eastern Ordos Basin[D].Xi’an：Northwestern University，2011.
[5]夏日元，唐健生.碳酸盐岩油气田古岩溶研究及其在油气勘探开发中的应用[J].地球学报，2006，27（5）：503-509.XIA Ri-yuan，TANG Jian-sheng.Palaeo-karst research of the carbonate oil-gas field and its application to oil-gas exploration and development[J].Acta Geoscientia Sinica，2006，27（5）：503-509.
[6]魏新善，任军峰，赵俊兴，等.鄂尔多斯盆地东部奥陶系风化壳古地貌特征嬗变及地质意义[J].石油学报，2017，38（9）：999-1009.WEI Xin-shan，REN Jun-feng，ZHAO Jun-xing，et al.Paleo-geomorphologic characteristic evolution and geological significance of the Ordovician weathering crust in eastern Ordos Basin[J].Acta Petrolei Sinica，2013，38（9）：999-1009.
[7]朱 静，李文厚，韩永林，等.陇东地区前侏罗纪古地貌及油藏成藏规律研究[J].沉积学报，2010，28（6）：1229-1237.ZHU Jing，LI Wen-hou，HAN Yong-lin，et al.Research on palaeogeomorphic features of Pre-Jurassic and accumulation regularity of oil reservoir in the area of eastern Gansu，Ordos Basin[J].Acta Sedimentologica Sinica，2010，28（6）：1229-1237.
[8]杨 华，刘新社，张道锋.鄂尔多斯盆地奥陶系海相碳酸盐岩天然气成藏主控因素及勘探进展[J].天然气工业，2013，33（5）：1-12.YANG Hua，LIU Xin-she，ZHANG Dao-feng.Main controlling factors of gas pooling in Ordovician marine carbonate reservoirs in the Ordos Basin and advances in gas exploration[J].Natural Gas Industry，2013，33（5）：1-12.
[9]赵俊兴，陈洪德，时志强.古地貌恢复技术方法及其研究意义[J].成都理工大学学报，2001，28（3）：260-266.ZHAO Xing-jun，CHEN Hong-de，SHI Zhi-qiang.The way and implication of rebuilding palaeogeomorphology-taking the research of palaeogeomorphology of the Ordos Basin before Jurassic deposition as example[J].Journal of Chengdu University of Technology，2001，28（3）：260-266.
[10]宋國奇，徐春华，樊庆真.应用层序地层学方法恢复加里东期古地貌以济阳坳陷沽化地区为例[J].石油实验地质，2000，22（4）：350-353.SONG Guo-qi，XU Chun-hua，FAN Qing-zhen.Seduence stratigraphy applied in the rsetoration of the Caledonian Paleogeomorphology：a case study in Zhanhua Area，the Jiyang Depression[J].Petroleum Geology & Experiment，2000，22（4）：350-353. [11]赵俊兴，陈洪德，向 芳.高分辨层序地层学方法在沉积前古地貌恢复中的应用[J].成都理工大学学报（自然科學版），2003，30（1）：76-82.ZHAO Xing-jun，CHEN Hong-de，XIANG Fang.The possibility of rebuilding paleogeomorphology before basin deposition by high resolution sequence stratigraphy[J].Journal of Chengdu University of Technology （Science & Technology Edition），2003，30（1）：76-82.
[12]赵永刚.奥陶系碳酸盐岩古岩溶及其储层研究[D].成都：西南石油大学，2006.ZHAO Yong-gang.Study on palaeo-rock-solution of carbonate rock of Ordovician and its reservoir[D].Chengdu：Southwest Petroleum University，2013，39（4）：516-529.
[13]张云峰，谭 飞，屈海洲，等.塔里木盆地轮古地区奥陶系风化壳岩溶储集层精细刻画[J].石油勘探与开发，2017，44（5）：1-11.ZHANG Yun-feng，TAN Fei，QU Hai-zhou，et al.Fine characterization of Ordovician weathering crust karst reservoirs in Lungu Area，Tarim Basin，China[J].Petroleum Exploration and Development，2017，44（5）：1-11.
[14]郭 庆.靖边气田西侧下古生界风化壳岩溶储层及天然气富集规律研究[D].西安：西北大学，2013.GUO Qing.Research on natural gas enrichment rule of the weathering crust karst reservoir of Lower Palaeozoic erathem in western area of Jingbian Gasfield[D].Xi’an：Northwestern University，2013.
[15]申赵军.鄂尔多斯盆地靖边气田西侧前石炭纪古地貌及古岩溶研究[D].成都：成都理工大学，2012.SHEN Zhao-jun.The research Pre-Carboniferous palaeogeomorphology and paleokarst the west side of the gas field of Jingbian in Ordos Basin[D].Chengdu：Chengdu University of Technology，2012.
[16]侯科锋，李浮萍，罗川又，等.苏里格气田东南部碳酸盐岩储层特征及分类评价[J].西安科技大学学报，2019，39（2）：276-285.HOU Ke-feng，LI Fu-ping，LUO Chuan-you，et al.Characteristics and classification of carbonate reservoir in southeastern Sulige Gasfield[J].Journal of Xi’an University of Science and Technology，2019，39（2）：276-285.
[17]刘 森.苏里格桃2区块马家沟组主力气藏储层地质特征研究[D].西安：西安石油大学，2017.LIU Sen.The research on reservoir geological characteristics of the major gas reservoirs of Majiagou Formation in Tao 2 block，Sulige[D].Xi’an：Xi’an Shiyou University，2017.
[18]李 鑫.鄂尔多斯盆地东南部奥陶系烃源岩评价[J].西安科技大学学报，2015，35（1）：56-62.LI Xin.Evaluation to Ordovician hydrocarbon source rocks in the southeast of the Ordos Basin[J].Journal of Xi’an University of Scienceand Technology，2015，35（1）：56-62.
[19]费世祥，张保国，李 洋，等.苏里格东南召41井区盒8段沉积微相研究[J].西安科技大学学报，2016，36（2）：228-234.FEI Shi-xiang，ZHANG Bao-guo，LI Yang，et al.Sedimentary microfacies of He 8 Section in Zhao 41 Area of southeast Sulige[J].Journal of Xi’an University of Science and Technology，2016，36（2）：228-234.
[20]王若谷，廖友运，尚 婷，等.苏里格气田东二区山西组砂岩储层特征及主控因素分析[J].西安科技大学学报，2016，36（3）：385-392.WANG Ruo-gu，LIAO You-yun，SHANG Ting，et al.Controlling factors and characteristics of sandstone reservoir of Shanxi Formation in the East Ⅱ Block of Sulige Gasfield[J].Journal of Xi’an University of Science and Technology，2016，36（3）：385-392. [21]杨伟利，王起琮，刘佳玮，等.鄂尔多斯盆地奥陶系马家沟组标准化层序地层学研究[J].西安科技大學学报，2017，37（2）：234-241.YANG Wei-li，WANG Qi-cong，LIU Jia-wei，et al.Standardization of sequence stratigraphy in Ordovician Majiagou Formation，Ordos Basin[J].Journal of Xi’an University of Science and Technology，2017，37（2）：234-241.
[22]杜彦军，贺永红，董丽红，等.鄂尔多斯盆地中部长8成藏主控因素[J].西安科技大学学报，2017，37（3）：388-395.DU Yan-jun，HE Yong-hong，DONG Li-hong，et al.Main control factors of tight oil reservoir in Chang 8 oil-bearing formation，central part of Ordos Basin[J].Journal of Xi’an University of Science and Technology，2017，37（3）：388-395.
[23]姜正龙，邓宏文，林会喜，等.古地貌恢复方法及应用——以济阳坳陷桩西地区沙二段为例[J].现代地质，2009，23（5）：865-871.JIANG Zheng-long，DENG Hong-wen，LIN Hui-xi，et al.Methods and application of paleo-geomorphologies rebuilding： an example of the second member of Shahejie Formation，Zhuangxi Area，Jiyang Depression[J].Geoscience，2009，23（5）：865-871.
[24]王高平，王震亮，赵雪娇，等.鄂尔多斯盆地延安地区奥陶系风化壳古地貌恢复[J].沉积学报，2013，31（4）：563-570.WANG Gao-ping，WANG Zhen-liang，ZHAO Xue-jiao，et al.Palaeogeomorphology restoring of Ordovician wea-thering crust in Yan’an Area，Ordos Basin[J].Acta Sedimentologica Sinica，2013，31（4）：563-570.
[25]James N P.Paleokarst[C]//In Choquette P W：Science.New York：Springer-Verlag，1988：416.
[26]Al Shaieb Z，Lynch M.Paleokarstic features and thermaloverprints observed in some of the Arbuckle coresin Oklahoma[C]//Paleo karst Related Hydrocarbon Reservoirs SEPM Core Workshop，Seattle：Society for Sedimentary Geology，1993：11-59.
[27]Wang B Q，Aasm I S.Karst-conotrilled diagenesis and reservoir development： example from the reservoir carbonate rocks on the eastern margin of the Ordos Basin，China[J].AAPG Bulletin，2002，86（9）：1639-1658.
[28]汪 洋，李树同，牟炜卫，等.乌审旗-志丹地区奥陶系岩溶古地貌与马五14气水分布关系[J].岩性油气藏，2016，28（2）：64-71.WANG Yang，LI Shu-tong，MU Wei-wei，et al.Effect of karst paleogeomorphology of Ordovician on gas-water distribution of Ma514 in Wushenqi-Zhidan Area[J].Lithologic Reservoirs，2016，28（2）：64-71.
[29]董月霞，赵宗举，曹中宏，等.南堡凹陷奥陶系碳酸盐岩岩溶残丘圈闭勘探潜力及意义[J].石油学报，2015，36（6）：653-663.DONG Yue-xia，ZHAO Zong-ju，CAO Zhong-hong，et al.Exploration potential and significance of Ordovician carbonate karst monadnock traps in Nanpu Sag[J].Acta Petrolei Sinica，2015，36（6）：653-663.
关键词：苏里格东南部;古地貌恢复;印模法;残厚法;气水分布特征
中图分类号：TE 122.2 文献标志码：A
DOI：10.13800/j.cnki.xakjdxxb.2019.0316 文章编号：1672-9315（2019）03-0498-09
Abstract：By the comprehensive use of sedimentary paleogeological analysis，impression method，and residual thickness method，the karst micro-geomorphology of the Pre-Carboniferous area in southeast Sulige，Ordos Basin was studied.Combining with actual production data，the characteristics of gas-water distribution was summarized.The study found that the karst paleogeomorphology at the end of Ordovician in the study area was represented by the distribution pattern of “high in the west and south，low in the east and north”，and developed two types of second-order paleogeomorphic units and nine types of third-order paleogeomorphic units.The karst slope is the main physiognomy unit，which is mainly develop on the upper slope，lower slope，primary monadnock and secondary monadnock.Moreover the ancient depression and groove are developed in the local area.The distribution of gas and water in the Majiagou-formation of study area is mainly characterized by the rich gas in the northwest，east and southeast，and the rich water in the central and western regions.Results show that：The karst microtopography in the study area can effectively depicted，by combing the carboniferous film thickness of the 8# coal seam at the top of Benxi Formation and the residual thickness of the lower datum of K3 marker.The monadnock and upper slope area are favorable natural gas enrichment areas.Water is affected by the tilting of the Yanshan Strata，and is mostly distributed in the central and western regions of the upper karst slopes.The anticorrosive and leaching effects of ancient depressions and ancient trenches are weak，but there are still a few stripper well distributed.Key words：southeastern part of Sulige Area;impression method;residual thickness method;characteristics of gas and water distribution0 引 言据统计20%～30%的大气田发育于古风化壳岩溶储层之中[1]，20世纪70年代中后期至今，与古岩溶相关的古地貌研究逐渐受到国内外学者重视，并得以广泛应用。古地貌恢复是对历史时期盆地面貌变化进行分析的一种重要手段，是以构造变形作用为主，风化剥蚀、充填沉积、差异压实等综合作用的结果[2]。 目前常用的古地貌刻画方法有：地球物理法、印模法、残余厚度法、沉积学法、层序地层法及层拉平法、回剥法、构造趋势面转化法、计算机模拟法等[2-11]。苏中堂、王建民等运用地球物理法，通过提升测井、地震属性，建立古地貌的识别标志模型，最后运用井震等地球物理资料对模型进行验证、约束，直至接近最理想的程度[2-4]。该方法直观准确，但刻画精度较低。夏日元、魏新善等运用印模法及残厚法恢复古地貌[5-8]。印模法目前应用最为广泛，其原理是利用上覆地层与古地貌之间的“镜像”关系，通过上覆地层的厚度来半定量恢复古地貌形态。该方法简单易行，但基准面与剥蚀面之间地层的真实厚度恢复难度较大。残余厚度法与印模法原理相似，区别在于需在剥蚀开始之前所沉积的地层中，选取标志层作为基准面，并将其拉平，该基准面之上残余地层厚度即可代表当时的古地貌形态。该方法直观且刻画精细，但却未考虑构造及剥蚀差异对沉积厚度的影响。赵俊兴等运用沉积学法恢复古地貌，即借助沉积岩相古地理特征，从而定性地恢复当时的古地貌特征[9]。沉积学法恢复古地貌简单、易于操作，但影响因素复杂、工作量大，且精度较低。宋国奇、赵俊兴等运用层序地层法，将层序地层学引入古地貌恢复之中，在建立上覆地层等时层序格架的基础上，将其拉平，运用其与风化壳等势面之间的地层来反应古地貌形态[10-11]。该方法优点是上覆基准面相对等时性强，缺点是2个等时界面的选择较为困难。不同的学者虽然采用不同的方法对不同地区的古地貌进行恢复研究，但各种分析方法均各有利弊。此外，就鄂爾多斯盆地而言，盆地尺度、油气田尺度的岩溶古地貌研究已趋于成熟[12-15]，但井区尺度的沉积期古地貌研究仍处于探索之中[12]。因此在前人研究的基础上，结合研究区实际，文中首次采用综合法对苏里格东南部地区三级岩溶古地貌进行精细划分与解释，即：井震结合初步确定地层下古生界地层格架;统计奥陶系顶部地层出露情况，了解古构造格局及各地区剥蚀程度;合理选择基准面，采用印模法及残厚法初步恢复古地貌;最终依照每种方法的优弊进行取舍，多角度综合划定研究区三级古地貌划分标准，恢复岩溶古地貌形态，并分析不同地貌形态与气水分布间的关系。
1 地质背景鄂尔多斯盆地受加里东运动影响，于奥陶纪盆地中央呈现“L”形隆起构造雏形[16]。早奥陶世，盆地在经历短暂抬升之后，开始接受马家沟组碳酸盐岩沉积;随着加里东运动逐渐加剧，华北地台整体抬升，中奥陶及早石炭地层被广泛剥蚀，发育风化壳岩溶型储集体，与上覆的晚石炭统海陆交互源岩沉积层呈区域不整合接触[7]。后经晚三叠世至白垩纪内陆盆地发展，以及新生代周边断陷盆地发育形成了现今南北隆升、东抬西冲、盆内斜坡西倾与天环向斜相连的构造面貌[14-17]。苏里格东南部位于鄂尔多斯盆地陕北斜坡北部（图1），地区面积790 km2，完钻井468口，其中探井19口，在平缓西倾的单斜背景上发育一系列北东-北东东向复式鼻状构造[18-22]。奥陶纪末，研究区位于中央古隆起东侧，差异性抬升剥蚀致使该区风化壳储层发育，表现出较大的勘探开发潜力。
2 常规法岩溶古地貌恢复2.1 前石炭古地质特征在对研究区217口钻至奥陶系单井进行小层对比与划分的基础上，统计每口钻井所揭示的奥陶系顶面出露层位（表1），绘制研究区前石炭系古地质图（图2），宏观分析石炭系沉积前的古地貌特征。
2.2 “印模”地层分布特征及古地貌指示在选取标志层参考面时，按照全区稳定分布等时界面、与风化壳面距离较近、岩性界面易于识别3大原则[9，13]，选取本溪组顶部8#煤层为基准面（图3），采用填平补齐的沉积补偿原理[23]，对地层抬升剥蚀后古地貌轮廓进行刻画。本溪组由顶部煤层与底部大套泥岩夹砂层组成，且全区广泛分布。考虑压实减孔对砂泥岩及煤系地层厚度影响较大[24]，因此先对本溪组厚度进行压实校正，具体如下
2.2.1 求取压实系数[25]
2.2.2 求取去压实后目标层顶底深度[24]
2.2.3 求取去压实后目标层骨架高度及压实率
式中 Z为深度，m;为Z处对应孔隙度，%;0为地表孔隙度（表2），%;C为压实系数;D1，D2为分别表示目标地层当前的顶深和底深，m;D1′，D2′分别表示压实校正后目标层定顶深度，m;D′为校正后目标层厚度，m;D为当前目标层厚度，m;Φ为压实率，%.经计算，该区本溪组砂岩平均压实率为28.7%，泥岩及煤层压实率为63%左右。在对本溪组地层进行去压实校正之后，得到本溪组原始沉积厚度等值线图，进而表征古地貌形态（图4）。可知，其整体反映了奥陶末期苏里格东南部西高东低、北高南低的地貌形态，与古地质图分析结果较为一致。西北部及东北部地层厚度相当，多分布于10～17 m之间，向南地层厚度逐渐增大。南部地区表现出沟脊相间的特征，桃2-36-2井—苏93井、桃2-30-14井—桃47井、陕50井—桃2-30-25井均呈现为北西-南东向展布特征，地层最厚处可达26.2 m.西部地区平均地层厚度高于东部地区，其中桃2-16-16井区、桃2-16-20井区、桃2-1-8井区厚度最大，均大于30 m，对应在本溪组沉积时表现为低洼特征。
2.3 “残厚”地层分布特征及古地貌指示与印模法类似，残厚法恢复古地貌亦需根据已述3个条件合理选择基准面。参考前石炭古地质
在不同程度剥蚀现象，因此选取K3标志层（马五14底部凝灰岩层）作为基准面，将其底拉平，该标志层之上的残余地层厚度代表了古地貌展布形态特征[21-22]（图5，图6）。
可知，古地貌整体西北偏高，起伏强烈;东南、西南地区仍表现出近北西-南东向沟脊相间形态，与印模法分析结果较为一致。在加里东期强烈构造作用的影响下，残厚小于40 m的区域主要位于研究区西部、北部地区，厚度变化幅度大且连续性差，桃2-21-8井处残余地层厚度最小，为24.6 m.由此说明，西部及北部地区遭受风化及剥蚀作用强烈，属低地貌单元。大于60 m的地层分布于研究区内中部地区，面积小且起伏程度小，向西逐渐歼灭，其中桃2-12-13C1井、桃35井、桃2-16-19井、桃2-30-8井K3标志层之上残厚度均大于65 m，并且均有马六地层残余，属高地貌单元。40～60 m地层分布最为广泛，占全区面积3/4，其间厚度变化幅度较小，说明其地势起伏较为平缓，遭受剥蚀程度较小。3 综合法古地貌单元划分及其特征运用沉积学古地质分析，虽能较好确定古地貌形态，但其精度不够[2，4]。印模法容易操作，并能较为精细的恢复古地貌。残厚法对古地貌的刻画直观，也较为精细，但基准面的选择、压实校正的准确度等都会影响最终结果[4，9]。上述分析结果也体现出相同结论，各种方法之间虽结果相似，能相互印证，但局部结果仍有差异。因此，文中在古地貌恢复技术的优缺点分析的基础上，以鄂尔多斯盆地沉积背景为根本，通过去压实校正，将印模法作为重点，综合古地质分析以及残厚分析结果，将研究区奥陶系岩溶古地貌进一步精细划分为2个二级单元和9类三级单元（图7），具体划分指标见表3. 3.1 岩溶高地苏里格东南部地区古岩溶台地面积较小，主要分布于北部桃2-3-11井-桃2-7-14井区、桃2-3-25井-桃2-4-27井区，奥陶系地层剥蚀强烈，顶部主要出露马五14及马五12地层，去压实校区后石炭系厚度一般小于15 m.
3.2 岩溶斜坡岩溶斜坡是该地区主体古地貌单元，去压实后石炭系厚度为15～50 m，K3标志层之上残余厚度为30～55 m.以石炭系及残余地层相对厚度为标准，岩溶斜坡又可进一步划分为6个三级单元。
残丘处石炭系厚度多小于30 m，是古岩溶斜坡上呈孤立丘状分布的古高地，分布范围较小，主要位于桃2-5-8井区，桃2-16-18井区，桃2-8-24-陕166井区，在加里东期所经历的风化剥蚀时间较长。上斜坡、下斜坡为是残丘及洼地间的过渡地带，沿残丘及洼地周缘呈环带状分布，奥陶系顶部出露层系以马五21为主。上斜坡处坡降较大，垂向渗流与横向岩溶均发育，物质交换迅速。
下斜坡处坡降较小，降水滞留，随所经历的成岩作用较为充分，但介质流动性差，不利于溶蚀物质被迅速带走。
洼地是岩溶斜坡内的负相地形区，残余地层厚度多小于10 m，面积较小，多位于沟槽水系发育末端。由于该区域地势低、水体深，因此水体流动缓慢，化学沉淀作用强，不利于储集体的形成与保存。古沟槽位于上斜坡边缘，为条带状展布。研究区内共3条古沟槽经过，东南部2条沟槽为北西-南东向展布，西部地区沟槽为北东-南西走向长达15 km，对下伏地层下侵作用强烈，致使该处马五1+2地层几乎被剥蚀殆尽，奥陶系顶部出露马五41至马五23.
4 岩溶古地貌与气水分布关系
鄂尔多斯盆地奥陶系马家沟组顶部的岩溶古地貌控制着古岩溶作用类型及强度，进而改变了风化壳古岩溶储层的岩性及物性条件[5]。但受地层叠覆顺序及风化作用强度对储层影响的强度所限，前石炭纪古地貌对奥陶系上组合地层影响较大[8，28]，因此，文中重点讨论研究区前石炭纪末古地貌与奥陶系马家沟组上组合气水分布间的关系。苏里格东南部地区气水分布并不均匀，产气井多分布于西北、中东及南部地区，产水井较少且零星分布，气水同产井集中分布于中部地区（图7）。马五1+2亚段北部及东南部产气量高，其中桃2-9-24井试气无阻流量高于100×104 m3，日产气量高达25×104 m3;马五14小层西北及东部地区产气量相对较高，其中桃2-6-1井试气无阻流量为69.8×104 m3，日产气量高达15×104 m3.其产水量高于马五1+2亚段，例如桃38井、桃2-30-14井等日产水量均可达9 m3以上。
叠合古地貌与上组合气水分布（图7），并结合前人关于岩溶储层特征研究成果可知，研究区前石炭纪古地貌形态对奥陶系上组合气水差异性分布影响较大。古岩溶高地在研究区内面积较小，垂向渗流为主，多发育纵向溶洞且非均一性明显[28]，不利于形成有益于天然气聚集成藏的层状岩溶带，仅局部地势相对较高的残丘地区尚有天然气分布。古岩溶斜坡是古地貌主体单元，其中残丘风化壳岩石类型以白云岩为主，主要表现为淡水淋滤溶蚀作用[23-31]，具有良好的储集性能;一级残丘、二级残丘各自所包含的气井占比40.0%，26.7%;36.9%的差气层井位于二级残丘之上，其次为一级残丘含差气层井26.0%.气水同产井多位于二级残丘及上斜坡单元，分别为31.8%，27.3%.产水井则集中分布于上斜坡、二级残丘及下斜坡古地貌单元之内，占比20.8%，45.8%，12.5%.此外，沟槽两侧分布有3.1%的差气层井。古洼地为岩溶水的汇集区[28]，岩溶作用相对较弱，所含气层、水层均较少。综上，奥陶纪末期苏里格东南部地区主要呈现西高东低、北高南低的地貌形态，古岩溶高地分布局限，古岩溶斜坡广泛展布，其间三级地貌单元高低相间分布。至晚三叠世、中-晚侏罗世，上古煤系烃源岩进入生排烃高峰，天然气下穿奥陶系顶部铝土岩缺失部位发生垂向运移[22]，在奥陶系风化壳岩溶储层裂缝及孔隙中向下充注并发生侧向运移，在当时地势相对较高的一级残丘、二级残丘、上斜坡等部位聚集成藏。至燕山运动中期，整个盆地由东倾单斜转为西倾单斜。最初在西部古残丘及上斜坡圈闭中聚集的天然气，受侧向遮挡而被保留;其余低部位聚集的天然气在自身浮力及水动力作用下，向东部构造高部位发生大面积后期运移与调整，驱替高部位地层水向西部底势区调整。最终出现现今中东部及前石炭古地貌呈岩溶高地貌单元的西北、南部中富气，中西部地区富水的气水分布格局。
5 结 论1）采用沉积古地质分析，选取本溪组顶部8#煤层为基准面的印模法，以K3标志层（马五14底部凝灰岩层）作为基准面的残余厚度法，均可对研究区前石炭纪古地貌进行有效恢复。各方法分析结果类似，可相互印证，但单一的沉积学分析、印模法或残厚法仍存在局限性。沉积古地质分析随操作简单，但精度较差;印模法虽能较好的反应古地貌格局，但去压实校正工作繁琐，对微地貌刻画反映较差;残厚法对微地貌刻画较为精细，但分析结果与整体格局稍有出入。
2）综合法恢复古地貌形态，将印模法与残厚法二者有机结合，并参考印模地层厚度、残余地层厚度与平均厚度值间的相对变化情况，结合研究区实际将前石炭古地貌精细刻画为2类二级单元和9类三级地貌单元。苏里格东南部地区奥陶系岩溶古地貌形态为“西高东低、北高南低”的分布格局，以岩溶斜坡地貌单元为主，主要发育上斜坡、下斜坡、一级残丘和二级残丘，局部发育古洼地和古沟槽。3）苏里格东南部地区奥陶系上组合气水主要表现为西北、东部及东南部富气，中西部富水的特征。水層、气水同层因受燕山期地层翘转影响，多分布于研究区中西部古岩溶上斜坡地区。古洼地、古沟槽溶蚀淋滤作用较弱，不利于天然气富集，但仍有少量低产气井分布其中。岩溶斜坡古残丘及上斜坡地区则是天然气高产井的目标区。
参考文献（References）：
[1]赵永刚，王东旭，冯强汉，等.油气田古地貌恢复方法研究进展[J].地球科学与环境报，2017，39（4）：516-529.ZHAO Yong-gang，WANG Dong-xu，FENG Qiang-han，et al.Review on palaeomorphologic reconstruction methods in oil and gas fields[J].Journal of Earth Science and Environment，2013，39（4）：516-529. [2]苏中堂.鄂尔多斯盆地塔巴庙地区奥陶系马家沟组沉积相与古岩溶储层特征[D].成都：成都理工大学，2008.SU Zhong-tang.Sedimentary facies and character of palaeokarst cust reservoir in Majiagou Formation in Ordovician，Tabamiao Area，Ordos[D].Chengdu：Chengdu University of Technology，2008.
[3]王建民，王佳媛，沙建怀，等.鄂尔多斯盆地东部奥陶系风化壳岩溶古地貌特征及综合地质模型[J].吉林大学学报（地球科学版），2014，44（2）：409-418.WANG Jian-min，WANG Jia-yuan，SHA Jian-huai，et al.Karst paleogeomorphology and comprehensive geological model of the Ordovician weathering crust in the eastern Ordos Basin[J].Journal of Jilin University （Earth Science Edition），2014，44（2）：409-418.
[4]王建民.鄂尔多斯盆地东部奥陶系风化壳岩溶古地貌与储层特征[D].西安：西北大学，2011.WANG Jian-min.Karst paleogeomorphology and reservoir characteristics of Ordovician weathered crust in eastern Ordos Basin[D].Xi’an：Northwestern University，2011.
[5]夏日元，唐健生.碳酸盐岩油气田古岩溶研究及其在油气勘探开发中的应用[J].地球学报，2006，27（5）：503-509.XIA Ri-yuan，TANG Jian-sheng.Palaeo-karst research of the carbonate oil-gas field and its application to oil-gas exploration and development[J].Acta Geoscientia Sinica，2006，27（5）：503-509.
[6]魏新善，任军峰，赵俊兴，等.鄂尔多斯盆地东部奥陶系风化壳古地貌特征嬗变及地质意义[J].石油学报，2017，38（9）：999-1009.WEI Xin-shan，REN Jun-feng，ZHAO Jun-xing，et al.Paleo-geomorphologic characteristic evolution and geological significance of the Ordovician weathering crust in eastern Ordos Basin[J].Acta Petrolei Sinica，2013，38（9）：999-1009.
[7]朱 静，李文厚，韩永林，等.陇东地区前侏罗纪古地貌及油藏成藏规律研究[J].沉积学报，2010，28（6）：1229-1237.ZHU Jing，LI Wen-hou，HAN Yong-lin，et al.Research on palaeogeomorphic features of Pre-Jurassic and accumulation regularity of oil reservoir in the area of eastern Gansu，Ordos Basin[J].Acta Sedimentologica Sinica，2010，28（6）：1229-1237.
[8]杨 华，刘新社，张道锋.鄂尔多斯盆地奥陶系海相碳酸盐岩天然气成藏主控因素及勘探进展[J].天然气工业，2013，33（5）：1-12.YANG Hua，LIU Xin-she，ZHANG Dao-feng.Main controlling factors of gas pooling in Ordovician marine carbonate reservoirs in the Ordos Basin and advances in gas exploration[J].Natural Gas Industry，2013，33（5）：1-12.
[9]赵俊兴，陈洪德，时志强.古地貌恢复技术方法及其研究意义[J].成都理工大学学报，2001，28（3）：260-266.ZHAO Xing-jun，CHEN Hong-de，SHI Zhi-qiang.The way and implication of rebuilding palaeogeomorphology-taking the research of palaeogeomorphology of the Ordos Basin before Jurassic deposition as example[J].Journal of Chengdu University of Technology，2001，28（3）：260-266.
[10]宋國奇，徐春华，樊庆真.应用层序地层学方法恢复加里东期古地貌以济阳坳陷沽化地区为例[J].石油实验地质，2000，22（4）：350-353.SONG Guo-qi，XU Chun-hua，FAN Qing-zhen.Seduence stratigraphy applied in the rsetoration of the Caledonian Paleogeomorphology：a case study in Zhanhua Area，the Jiyang Depression[J].Petroleum Geology & Experiment，2000，22（4）：350-353. [11]赵俊兴，陈洪德，向 芳.高分辨层序地层学方法在沉积前古地貌恢复中的应用[J].成都理工大学学报（自然科學版），2003，30（1）：76-82.ZHAO Xing-jun，CHEN Hong-de，XIANG Fang.The possibility of rebuilding paleogeomorphology before basin deposition by high resolution sequence stratigraphy[J].Journal of Chengdu University of Technology （Science & Technology Edition），2003，30（1）：76-82.
[12]赵永刚.奥陶系碳酸盐岩古岩溶及其储层研究[D].成都：西南石油大学，2006.ZHAO Yong-gang.Study on palaeo-rock-solution of carbonate rock of Ordovician and its reservoir[D].Chengdu：Southwest Petroleum University，2013，39（4）：516-529.
[13]张云峰，谭 飞，屈海洲，等.塔里木盆地轮古地区奥陶系风化壳岩溶储集层精细刻画[J].石油勘探与开发，2017，44（5）：1-11.ZHANG Yun-feng，TAN Fei，QU Hai-zhou，et al.Fine characterization of Ordovician weathering crust karst reservoirs in Lungu Area，Tarim Basin，China[J].Petroleum Exploration and Development，2017，44（5）：1-11.
[14]郭 庆.靖边气田西侧下古生界风化壳岩溶储层及天然气富集规律研究[D].西安：西北大学，2013.GUO Qing.Research on natural gas enrichment rule of the weathering crust karst reservoir of Lower Palaeozoic erathem in western area of Jingbian Gasfield[D].Xi’an：Northwestern University，2013.
[15]申赵军.鄂尔多斯盆地靖边气田西侧前石炭纪古地貌及古岩溶研究[D].成都：成都理工大学，2012.SHEN Zhao-jun.The research Pre-Carboniferous palaeogeomorphology and paleokarst the west side of the gas field of Jingbian in Ordos Basin[D].Chengdu：Chengdu University of Technology，2012.
[16]侯科锋，李浮萍，罗川又，等.苏里格气田东南部碳酸盐岩储层特征及分类评价[J].西安科技大学学报，2019，39（2）：276-285.HOU Ke-feng，LI Fu-ping，LUO Chuan-you，et al.Characteristics and classification of carbonate reservoir in southeastern Sulige Gasfield[J].Journal of Xi’an University of Science and Technology，2019，39（2）：276-285.
[17]刘 森.苏里格桃2区块马家沟组主力气藏储层地质特征研究[D].西安：西安石油大学，2017.LIU Sen.The research on reservoir geological characteristics of the major gas reservoirs of Majiagou Formation in Tao 2 block，Sulige[D].Xi’an：Xi’an Shiyou University，2017.
[18]李 鑫.鄂尔多斯盆地东南部奥陶系烃源岩评价[J].西安科技大学学报，2015，35（1）：56-62.LI Xin.Evaluation to Ordovician hydrocarbon source rocks in the southeast of the Ordos Basin[J].Journal of Xi’an University of Scienceand Technology，2015，35（1）：56-62.
[19]费世祥，张保国，李 洋，等.苏里格东南召41井区盒8段沉积微相研究[J].西安科技大学学报，2016，36（2）：228-234.FEI Shi-xiang，ZHANG Bao-guo，LI Yang，et al.Sedimentary microfacies of He 8 Section in Zhao 41 Area of southeast Sulige[J].Journal of Xi’an University of Science and Technology，2016，36（2）：228-234.
[20]王若谷，廖友运，尚 婷，等.苏里格气田东二区山西组砂岩储层特征及主控因素分析[J].西安科技大学学报，2016，36（3）：385-392.WANG Ruo-gu，LIAO You-yun，SHANG Ting，et al.Controlling factors and characteristics of sandstone reservoir of Shanxi Formation in the East Ⅱ Block of Sulige Gasfield[J].Journal of Xi’an University of Science and Technology，2016，36（3）：385-392. [21]杨伟利，王起琮，刘佳玮，等.鄂尔多斯盆地奥陶系马家沟组标准化层序地层学研究[J].西安科技大學学报，2017，37（2）：234-241.YANG Wei-li，WANG Qi-cong，LIU Jia-wei，et al.Standardization of sequence stratigraphy in Ordovician Majiagou Formation，Ordos Basin[J].Journal of Xi’an University of Science and Technology，2017，37（2）：234-241.
[22]杜彦军，贺永红，董丽红，等.鄂尔多斯盆地中部长8成藏主控因素[J].西安科技大学学报，2017，37（3）：388-395.DU Yan-jun，HE Yong-hong，DONG Li-hong，et al.Main control factors of tight oil reservoir in Chang 8 oil-bearing formation，central part of Ordos Basin[J].Journal of Xi’an University of Science and Technology，2017，37（3）：388-395.
[23]姜正龙，邓宏文，林会喜，等.古地貌恢复方法及应用——以济阳坳陷桩西地区沙二段为例[J].现代地质，2009，23（5）：865-871.JIANG Zheng-long，DENG Hong-wen，LIN Hui-xi，et al.Methods and application of paleo-geomorphologies rebuilding： an example of the second member of Shahejie Formation，Zhuangxi Area，Jiyang Depression[J].Geoscience，2009，23（5）：865-871.
[24]王高平，王震亮，赵雪娇，等.鄂尔多斯盆地延安地区奥陶系风化壳古地貌恢复[J].沉积学报，2013，31（4）：563-570.WANG Gao-ping，WANG Zhen-liang，ZHAO Xue-jiao，et al.Palaeogeomorphology restoring of Ordovician wea-thering crust in Yan’an Area，Ordos Basin[J].Acta Sedimentologica Sinica，2013，31（4）：563-570.
[25]James N P.Paleokarst[C]//In Choquette P W：Science.New York：Springer-Verlag，1988：416.
[26]Al Shaieb Z，Lynch M.Paleokarstic features and thermaloverprints observed in some of the Arbuckle coresin Oklahoma[C]//Paleo karst Related Hydrocarbon Reservoirs SEPM Core Workshop，Seattle：Society for Sedimentary Geology，1993：11-59.
[27]Wang B Q，Aasm I S.Karst-conotrilled diagenesis and reservoir development： example from the reservoir carbonate rocks on the eastern margin of the Ordos Basin，China[J].AAPG Bulletin，2002，86（9）：1639-1658.
[28]汪 洋，李树同，牟炜卫，等.乌审旗-志丹地区奥陶系岩溶古地貌与马五14气水分布关系[J].岩性油气藏，2016，28（2）：64-71.WANG Yang，LI Shu-tong，MU Wei-wei，et al.Effect of karst paleogeomorphology of Ordovician on gas-water distribution of Ma514 in Wushenqi-Zhidan Area[J].Lithologic Reservoirs，2016，28（2）：64-71.
[29]董月霞，赵宗举，曹中宏，等.南堡凹陷奥陶系碳酸盐岩岩溶残丘圈闭勘探潜力及意义[J].石油学报，2015，36（6）：653-663.DONG Yue-xia，ZHAO Zong-ju，CAO Zhong-hong，et al.Exploration potential and significance of Ordovician carbonate karst monadnock traps in Nanpu Sag[J].Acta Petrolei Sinica，2015，36（6）：653-663.