0 引言
中国西部塔里木盆地油气资源丰富,油气藏类型较多,其中碳酸盐岩缝洞型油气藏是主要的油气藏类型之一。关于塔里木盆地碳酸盐岩岩溶缝洞储层的研究一直是石油地质领域的热点和难点,特别是在缝洞储层的正演模拟方面前人做了大量工作,也得到了许多重要认识。胡中平[1]、吴俊峰等[2]、李凡异等[3]先后对缝洞体形成“串珠状”的机理及形成条件进行了分析,并进一步讨论了影响“串珠状”反射能量强弱的因素;闵小刚等[4]基于非均匀介质正演模拟技术,总结了溶洞储集体在叠加偏移剖面上的反射特征;姚姚等[5]利用随机模型模拟得出用常规频带地震数据可以探测一定规模的溶洞;孙东等[6]采用有限差分法求解波动方程,研究了塔中地区碳酸盐岩洞穴型储集层溶洞与串珠状反射的位置对应关系、多种串珠形态的成因机制;马灵伟等[7]定量分析了缝洞体储层纵横向发育规模、填充物类型及背景围岩对缝洞储集体反射特征的影响;前人主要从不同规模溶洞、储层各向异性、纵向多层洞穴、缝洞联合体等多个方面开展了碳酸盐岩储层反射特征的正演模拟研究,所建立的地质模型主要以球形特征的洞穴为主,对于呈管状结构的地下暗河的地震正演模拟研究很少。
塔里木盆地的轮古油田其主体为碳酸盐岩暗河型缝洞油藏[8]。该区奥陶系为典型的碳酸盐岩潜山,受长期强烈的风化、淋滤、剥蚀和溶蚀作用,发育大量与古暗河相关的缝洞储层,储层的非均质性非常强,该区的油气也主要富集在这些缝洞储层中[9]。由于前期对该区奥陶系古暗河系统的发育规律、暗河充填程度、暗河储层地震地质特征认识不清,造成该区在实际钻探中50%的井因未钻遇优质缝洞储层而失利。因此在轮古地区开展暗河系统特征的研究是十分的必要。本次研究在充分吸收前人认识的基础上,通过建立不同类型的暗河三维地质模型,利用三维正演数值模拟技术,研究不同长度、高度、充填物类型等因素对暗河地震响应特征的影响,为在地震上预测暗河的空间分布、寻找优质的与暗河相关的缝洞储层提供了一定的依据。
1 暗河地震特征
塔里木盆地轮古地区在地质年代史上经历了多期构造运动,特别是受到早海西期多次构造抬升作用[10],造成奥陶系碳酸盐岩上覆的泥盆系、志留系及中上奥陶统被剥蚀,同时受大气淡水淋滤作用的影响,形成了古潜山岩溶地貌,碳酸盐岩潜山内幕受大气淡水沿断裂、裂缝溶蚀作用,形成大量大小、长度不一的暗河和孤立洞穴,尤其暗河最为发育[11]。从现代岩溶发育特征来看,暗河是指碳酸盐岩中发育的地下河,是由于地表径流沿岩石裂隙渗入地下,岩石经过溶蚀、坍塌以及水的搬运,在地下形成了大小不同、长短不一、错综复杂的管道系统,而每一个洞穴也是一个大的暗河系统的一部分。洞穴不同部位宽度、高度变化较大,洞穴内低洼位置易被泥质充填[12]。
在塔里木盆地轮古地区有大量的钻井在奥陶系碳酸盐岩中钻遇了洞穴,洞穴里大多被砂泥质充填。轮古西地区的LG42钻井显示在井深5 810~5 831.5 m段钻遇一个砂泥岩充填的大洞穴,井眼处钻遇洞穴高度达21.5 m(图1a)。从测井曲线上看,钻遇暗河的位置测井曲线表现为高伽马、低电阻、低密度和高声波时差的特点[13](图1a),与围岩致密灰岩的测井曲线值有明显的差异。通过井震标定,发现井上钻遇洞穴位置在地震剖面上表现为强反射地震特征(图1b、d),且在地震上横向连续延伸较远,说明暗河区域与周围致密灰岩相比,存在较大的速度差,所以才在地震上表现为强反射的界面。在波阻抗剖面上,暗河表现为低阻抗特征(图1c),表明暗河内存在较多的空隙空间或低速岩性充填,造成暗河地层速度远小于周围致密灰岩的速度。
图1
图1
LG42井暗河井震标定
a—LG42井奥陶系柱状图;b—过LG42井地震剖面;c—过轮古42井反演剖面;d—LG42井区奥陶系潜山反演属性平面
Fig.1
Calibrating downhole seismic chart of underground river of well LG42
a—columnar chart of ordovician of well LG42;b—seismic sections through well LG42;c—inversion section through well LG42;d—the inversion attribute plane map of Ordovician buried hill in LG42 well area
根据目前已钻井钻遇洞穴的地震特征和地震属性预测,在轮古地区主要识别出3种暗河地震相特征: 连续强反射、杂乱点状强反射、杂乱点状弱反射(图2a)。
图2
图2
过暗河地震剖面(a)和LG15-2井区奥陶系潜山振幅属性平面(b)
Fig.2
Seismic profile along the underground rivers (a) and planar graph of amplitude attribute of Ordovician buried hill in LG15-2 well area(b)
2)杂乱点状强反射地震相特征。在地震剖面上表现为多个点状强反射纵向上杂乱分布,平面属性却在一条线上,该类地震相特征可能反映了沿暗河纵向上多层洞穴的特征(图2a、b)。这种特征暗河,钻井揭示暗河纵向规模较大,发育多层洞穴。
3)杂乱点状弱反射地震相特征。在地震剖面上表现为杂乱弱反射,产状杂乱,无明显特征,仅振幅比围岩略有增强,与围岩差异较小(图2a),平面属性上略增强的振幅也在一条线上,这类暗河钻井揭示暗河规模较小,以裂缝孔洞储层为主。
2 暗河系统地震正演模拟
为了进一步论证深埋于地下的暗河在地震资料上的反射特征,开展了暗河的三维正演模拟。在对研究区碳酸盐岩暗河地震相特征分析的基础上,采用基于井震标定的方法,建立符合轮古地区奥陶系古暗河系统典型地质模型,采用多规模、多方位暗河模拟条件,探讨不同类型地震相的成因,为该类储层预测由定性向定量发展提供理论指导。
2.1 暗河不同段宽度变化的三维正演模拟
图3为暗河三维正演模型和暗河三维正演速度模型。为了更能真实模拟深埋5 000 m以下的古暗河的地震响应特征,模型设计时完全参考实际钻井揭示的地层参数,暗河上覆地层和围岩模型完全按测井真实地层速度和地层厚度进行速度建模。设计的暗河长度为2.43 km,宽度为10~50 m,暗河洞穴高度都为25 m,暗河洞穴里设计为砂泥质充填,充填速度为4 000 m/s,暗河围岩的速度为6 100~6 300 m/s,暗河埋深设计跟研究区真实暗河埋深一致(5 500~6 000 m)。整个三维正演方法采用完全模拟真实地震采集的波动方程正演。图4为暗河三维正演模拟地震成果剖面,从图中可以看出,暗河发育的区域,地震剖面上表现为强波峰和强波谷特征(图4中①);在暗河河道宽的位置,表现为连续性强振幅特征,在横切暗河宽河道地震剖面上,地震上表现为强“串珠状”特征(图4中③);在洞穴狭窄部位振幅较弱,在横切暗河窄河道地震剖面上,地震上表现为弱“串珠”杂乱反射特征(图4中②)。为了进一步分析什么属性对地震剖面上暗河的反射识别最敏感,利用正演地震成果数据提取了均方根振幅、峰值频谱频率、平均瞬时相位等3种属性(图5),在正演地震成果均方根振幅属性平面图上(图5b),暗河位置表现为线性的强反射特征,振幅与围岩基质地震反射有明显差异,振幅属性平面上暗河轮廓与暗河地质模型平面特征十分吻合。这也进一步证明了深埋地下的古老暗河相对于围岩在地震上反映为连续性强振幅反射。在正演成果峰值频谱频率属性平面图和平均瞬时相位属性平面图上(图5c、d),这两种属性反映的暗河形态特征都十分清楚,但对暗河形态轮廓都有放大效应。暗河横向宽度大的地方,地震峰值频率高,地震平均瞬时相位低;暗河横向宽度变窄的地方,地震峰值频率低,地震平均瞬时相位变高。暗河与围岩交界处呈现一圈地震峰值频率最低、地震瞬时相位反相的特征,这种振幅、频率、相位变化特征为解释人员识别暗河空间分布特征提供了依据。
图3
图3
三维正演速度模型(a)和暗河三维地质模型(b)
Fig.3
Three dimensional forward velocity model(a)and three dimensional geological model of underground river(b)
图4
图4
正演成果地震剖面(a)、(b)和正演成果振幅属性平面(c)
Fig.4
Seismic profile of forward results(a)、(b)and planar map of amplitude attribute of forward results(c)
图5
图5
暗河地质模型平面图和正演地震成果振幅、频率、相位属性平面
a—暗河地质模型平面;b—正演成果振幅属性平面;c—正演成果峰值频谱频率属性平面;d—正演成果平均瞬时相位属性平面
Fig.5
Plan of underground river geological model and plan of amplitude, frequency and phase attributes of forward seismic results
a—geological model;b—planar map of amplitude attribute;c—planar map of peak spectral frequency;d—planar map of average instantaneous phase
2.2 不同充填程度、纵向多层暗河的正演模拟
根据现代岩溶发育特征可知,碳酸盐岩地区发育的地下暗河结构十分复杂,大小从不到1 m到上百米,空间上暗河多层叠置。并且深埋地下的古暗河内部大多被充填,有的地方被砂泥质充填,有的地方被垮塌的灰岩角砾岩充填,暗河内部非均质性强[14]。为了能更真实地模拟深埋地下的暗河特征,我们建立了一条形态变化大、各段充填程度不同、纵向多层的暗河模型(图6a),暗河最宽洞穴段高宽均为40 m,暗河最狭窄段洞穴高宽均为5 m。图中暗河不同的颜色代表洞穴充填速度不同,充填速度由3 000~6 000 m/s变化,代表未充填、砂泥质充填、灰岩角砾充填、灰岩破碎带等不同物性的充填。图6b为对应的正演成果地震剖面,从地震剖面上可以看出,暗河洞穴大、充填速度低的部位,地震响应最强,表现为连续性强振幅特征;暗河洞穴狭小部位,即使洞穴未充填,在地震剖面上地震响应特征也不明显,难以识别。当多层暗河洞穴纵向叠置时,如图6a地质模型中间位置,三层暗河洞穴纵向叠置,受地震纵向分辨率的限制,地震上仅表现为多了一个强波峰强波谷,从地震上难以分出这三层洞穴的界限。总的来说,暗河洞穴地震反射波实际上是暗河顶底面反射波及顶底面之间的多次波复合的结果[15],地震振幅越强,表明对应的暗河洞穴越大、充填物速度越低,在碳酸盐岩内幕属性平面上成线性分布的强振幅异常反射,可能是一条暗河,只是暗河不同段大小、充填情况不一样。
图6
图6
不同充填程度暗河地质模型和正演成果地震剖面
a—正演地质模型;b—正演地震剖面
Fig.6
Geological model of underground river with different filling degree and seismic profile of forward modeling results
a—forward geological model;b—forward seismic profile
2.3 暗河长度变化与地震反射特征关系
暗河空间纵、横向尺度变化会引起暗河的地震反射的振幅、频率、相位变化,研究这些参数与暗河空间大小之间的关系具有十分重要的意义,对实际地震资料中暗河型缝洞储层参数的估算及储层评价能够提供一定的依据。
图7为高度相同长度变化的暗河模型正演与对应地震属性曲线。在这次正演中,共设计了14个由短到长变化的圆形管道状暗河地质模型(图7a),固定暗河管道的直径(15 m),只变化暗河管道的长度。暗河的地质模型长度从左到右的变化分别是:15、25、35、60、75、100、125、140、160、200、240、315、500、630 m。暗河充填速度为4000 m/s,背景围岩速度为6 200~6 300 m/s,暗河埋深为6 000 m。图7b为管道状暗河正演地震偏移叠加剖面,从正演地震剖面中可以看出随着暗河长度的增加,暗河对应位置的地震特征有较大的差异。为了直观地描述暗河长度的变化对地震波场特征的影响,对理论模型中提取的地震属性参数与模型地震响应进行了相关对比分析(图7c、d、e)。从图7c可以看到:当暗河长度在100~125 m时,暗河对应的地震振幅最强(3 050);暗河长度小于125 m时,随着暗河长度的减小,暗河对应的地震振幅在减小;当暗河长度在15 m时,对应的地震振幅(900)减小到与地震背景振幅(850)相差无几,在地震上识别出暗河特征较为困难;暗河长度由125 m变长到200 m时,地震振幅在变小;暗河长度大于200 m后,地震振幅趋于稳定,振幅值为2700。在峰值频谱频率曲线图上(图7d),暗河长度从15 m增加到100 m时,暗河对应的地震峰值频谱频率在减小,暗河长度由100 m增加到200 m时,暗河对应的地震峰值频谱频率在增大,暗河长度大于200 m后,暗河对应的地震峰值频谱频率基本为一恒定值;在平均瞬时相位曲线图上(图7e),相位的变化与暗河地震峰值频谱频率具有相似的变化特征。暗河长度在地震记录上的分辨率取决于地震资料偏移处理的效果,由于一个绕射波偏移后的理想效果是能量集中在一个以1/4的波长为半径的圆内,也即地震数据实施偏移处理后的第一费涅耳带的半径为1/4的波长。在本次正演中,上覆地层平均速度为3 432 m/s,频率35 Hz时,波长可达98 m,在此情况下暗河的横向分辨率为49 m,也就是49 m长的暗河是可以准确识别的,图7c正演暗河地震振幅曲线上也可以看出,从暗河长度为60 m的模型开始振幅曲线的半幅点位置的宽度与暗河的长度是一致的。从图7c、图7d、图7e中可以看出振幅、频率、相位3种属性对应不同长度的暗河都有属性特征的变化,都可以用来识别暗河,但振幅属性的识别精度最高。
图7
图7
高度相同长度变化的暗河模型正演与对应地震属性曲线
a—地质模型;b—正演成果剖面;c—暗河地震反射振幅曲线;d—暗河地震反射峰值频谱频率曲线;e—暗河地震反射平均瞬时相位曲线
Fig.7
Forward modeling of underground river model with the same height and length variation and corresponding seismic attribute curve
a—geological model;b—forward modeling result;c—seismic reflection amplitude curve of underground rivers;d—peak spectrum frequency curve of underground rivers seismic reflection;e—average instantaneous phase curve of underground rivers seismic reflection
3 暗河系统地震识别与储层分布
3.1 暗河系统地震识别
前面利用井震标定和正演模拟,证实了碳酸盐岩中深埋地下的暗河在地震上表现为连续性强反射特征,地震上平均瞬时相位、峰值频谱频率、均方根振幅等属性都可以反映暗河空间的展布,但识别暗河精度最高、效果最好的地震属性还是均方根振幅属性。采用均方根振幅属性对轮古西—轮古7地区奥陶系暗河进行预测。图8为轮古西—轮古7地区奥陶系潜山以下0~100 ms振幅均方根属性平面,图上连续的呈管道状的强振幅和呈点状有规律线性展布的强串珠均是暗河,连续强振幅特征的暗河表明该段暗河管道大,暗河内可能仅为低速充填物(泥岩、砂岩)或未充填;振幅时强时弱,但平面上仍是呈一条线性分布,这样的暗河在振幅弱的地方可能是暗河管道小或暗河内部充填了垮塌的灰岩角砾,但仍是有暗河管道存在,在振幅强的位置,如果纵向发育多套强波峰波谷,该段暗河可能就是此区域暗河洞穴纵横向规模变大了,或者是多层洞穴叠置区域。在轮古西—轮古7地区,利用均方根振幅属性,识别出95条相暗河系统,为该区下步针对暗河型油藏的滚动开发奠定了基础。
图8
图8
轮古西—轮古7地区奥陶系潜山以下0~100 ms地震振幅均方根属性平面
Fig.8
Plan of root mean square attribute of 0~100 ms seismic amplitude under Ordovician buried hill in Lunguxi-Lungu 7 area
3.2 暗河系统储层分布
根据现代岩溶特征,通常暗河具有独立的补给、径流和排泄系统[16],大的暗河可形成地下河系,跟地表水系一样,通常有主河道和分支河道,具有很大的流量[17]。在轮古油田奥陶系潜山区,暗河大量发育,一些大的暗河十分复杂,形成暗河系统。图9为轮古15井区暗河系统平面、剖面图,从图中可以看出,该暗河系统是由主河道和分支河道组成,主暗河延伸较长,规模较大,如图9中①剖面,主暗河在整个暗河系统中所处位置最低,暗河管道大,地震剖面和振幅属性平面上都表现为连续的强振幅特征;暗河分支延伸短,位置较浅,规模小,部分暗河管道狭窄,地震剖面和属性平面图上表现为线性点状强反射,如图9中②剖面。每个暗河系统都有多个入水口和出水口,入水口主要在主暗河的高部位和暗河分支的高部位,出水口主要分布在主暗河低部位,暗河水最终汇聚到主暗河里流入明河。
图9
图9
轮古15井区暗河系统平面、剖面
a—LG15井区奥陶系潜山振幅属性平面;b—过主暗河地震剖面;c—过分支暗河地震剖面
Fig.9
Plan and section of underground river system in Lungu 15 well block
a—rms amplitude attribute map;b—seismic profile of the main underground river;c—seismic profile of the branches underground river
根据暗河系统的发育规律分析,暗河系统储层可能发育的位置为暗河主河道、暗河分支河道以及暗河洞顶破碎带或叫暗河洞顶缝[18]。暗河主河道是整个暗河系统水流的主要通道[19],由于其位于整个暗河系统的低部位,易被砂泥质充填。以轮古地区的LG15-2井所钻探的暗河主河道为例(图9的①)。暗河主河道上LG15-2井在井深5 810~5 843 m处钻遇一大型充填洞穴,洞穴高达33 m,洞内被砂泥质充填,储集空间有限,造成这口井产油较少,为低效井。统计显示轮古油田有11口井钻遇暗河主河道上,均被砂泥质充填,产油气较少,表明暗河主河道由于被砂泥质充填不是暗河储层发育的有利位置[20]。暗河的分支河道是主暗河的各个分支,也是主暗河水流的补给处。暗河分支河道位于整个暗河系统的较高位置,主要为暗河的水流补给通道,水动力强,不易被泥质充填,洞穴储层更为发育,也最容易成藏[21]。图9c为过分支暗河地震剖面,在该分支暗河上,钻探了2口开发井LG15-1、LG15-3,LG15-1位于分支暗河高部位,构造位置高且储层最发育,累产油23万t;LG15-3井构造位置相对于LG15-1井稍低,累产油6万t,从产能来看,说明LG15-1井所在的暗河储层十分发育,暗河没有被泥质充填,储集空间大。轮古15井区暗河系统上共钻探了5口井,开发过程中针对同一暗河系统上的井进行注水替油实验,发现LG15-2注水,分支暗河上LG15-6、LG15-3、LG15等井均出现了压力变化,说明分支暗河河道并没有完全充填,是可以通过暗河河道连通的。
暗河洞顶缝、洞顶破碎带也是暗河有利储层发育区[22]。轮古地区的LG100井在井深5 506~5 515 m处钻遇大型洞穴,但洞穴里岩性为泥质粉砂岩,已经被填充,为暗河河道位置。但在洞穴上方发育近40 m的Ⅱ、Ⅲ类洞顶缝储层,以孔洞型和裂缝孔洞型为主。洞顶缝储层与下部的洞穴至今均产油气,累计产油7.9 万t。这就表明:暗河洞顶破碎带、洞顶缝也是有利的储集空间。
4 认识与结论
通过正演模拟并结合地震、地质对轮古西地区奥陶系古岩溶暗河系统的综合分析研究,得出以下结论:
1)在致密灰岩背景下发育的地下暗河系统在地震剖面上表现为地震同相轴连续强反射特征,振幅属性平面图上呈线性展布;地震剖面上地震同相轴为点状强反射,振幅属性平面图上连续多个点状强反射成线性分布,这样的特征也是暗河的响应,只是暗河不同段大小、充填情况不一样。
2)模型正演显示暗河的长度、宽度的变化影响其地震振幅强弱的变化,暗河宽度一定的情况下,长度在一个波长以下时,地震响应振幅随暗河长度增加而增强,长度在一个波长时,振幅最强,随后随着长度增加振幅减弱,长度超过一个半波长后,地震响应振幅趋于稳定。暗河长度大于1/2地震波长时,地震资料可以较准确地刻画出暗河的长度。单纯根据地震振幅的强弱,不能判断暗河的孔隙度大小,但地震上暗河位置振幅越强,其所反映的储集空间也越大。
3)振幅属性可以较准确地刻画地震上暗河的横向范围,而频率、相位属性可以刻画地震上暗河的轮廓,但刻画的范围比实际暗河要大。
4)暗河系统与地表明河系统相似,主河道位置低,为汇水区,晚期易被泥质充填;分支暗河位置高、坡度陡,被充填的概率低,是暗河储层发育的主要部位,是油气富集的主要空间。
参考文献
溶洞地震波“串珠状”形成机理及识别方法
[J].
Mechanism and distinction method forthe seismic "string beads" characteristic
[J].
碳酸盐岩特殊孔洞型构造地震响应特征分析
[J].
Analysis on seismic response of special cavernous structure of carbonate
[J].
碳酸盐岩溶洞的“串珠”状地震反射特征形成机理研究
[J].
Research on the mechanism of the formation of the "beaded" seismic reflection characteristics of carbonate karst caves
[J].
塔河油田孔洞模型的波动方程正演模拟
[J].
Wave equation forward modeling of cavern models in Tahe oil-field
[J].
深层碳酸盐岩岩溶风化壳洞缝型油气藏可检测性的理论研究
[J].
Theoretical study on the detectability of deep carbonate karst weathered crust-cavity fractured oil and gas reservoirs
[J].
碳酸盐岩储层垂向长串珠形成机制
[J].
Formation mechanism of vertically long beads in carbonate reservoirs
[J].
应用随机介质正演模拟刻画深水区台缘礁碳酸盐岩储层
[J].
Using forward modeling of stochastic media to characterize carbonate reservoirs in platform margin reefs in deep water
[J].
轮古油田奥陶系潜山洞穴型储层发育特征及油气分布控制因素
[J].
Characteristics of Ordovician buried-hill cave reservoir and controlling factors of petroleum distribution of Lungu oilfield
[J].
塔北隆起西部缝洞型碳酸盐岩储层表征与评价
[J].
Characterization and evaluation of fracture-cavity type carbonate reservoir in the western part of Northern Tarim uplift
[J].
向源潜流侵蚀岩溶作用及其成因机理——以塔河油田早海西风化壳岩溶洞穴层为例
[J].
The mechanism of retrogressive erosion and karstification:A case study of cave formation in early hersinian weathered crust in tahe oilfield
[J].
塔里木盆地塔河油田岩溶峡谷区海西早期洞穴系统发育模式
[J].
Development model of Hercynian cave system in karst canyon area of Tahe Oilfield,Tarim Basin
[J].
塔河油田奥陶系古岩溶洞穴识别及预测
[J].
Identification and prediction of Ordovician ancient karst caves in Tahe Oilfield
[J].
测井资料在碳酸盐岩洞—裂缝型储层产能评价中的应用
[J].
Application of well logging data in productivity evaluation of carbonate cave-fracture reservoir
[J].
塔里木盆地轮古潜山暗河发育特征及其与油气富集的关系
[J].
Development characteristics of underground river and its relationship with oil and gas accumulation in Lungu buried-hill,Tarim Basin
[J].
AVO技术在识别充填流体溶洞中的应用
[J].
The application of avo to the predication of water-eroded caves filled with liquids for carbonate reservoirs
[J].
塔河油田西部古暗河缝洞结构特征及控制因素
[J].
Analysis on fracture-cave structure characteristics and its controlling factor of palaeo-subterranean rivers in the western Tahe Oilfield
[J].
桂林地区灰岩洞穴的溶蚀形态
[J].
Dissolution features of limestone caves in Guilin area
[J].
塔河地区奥陶系深埋岩溶洞穴特征及保存机制初探
[J].
Characteristics and preservation mechanism of the Ordovician deep burial karst caves in Tahe area
[J].
塔河油田奥陶系碳酸盐岩溶洞型储集体识别及定量表征
[J].
Ordovician carbonate fracture-cavity reservoirs identification and quantitative characterization in Tahe Oilfield
[J].
塔河油田奥陶系碳酸盐岩油藏溶洞发育规律
[J].
Development regularity of karstic caverns of Ordovician carbonate reservoirs in Tahe oilfield
[J].
塔河油田奥陶系缝洞型碳酸盐岩油藏的储层连通性及其油(气)水分布关系
[J].
Reservoir connectivity and oil-water relationship of rock dissolved Carbonate oil reservoir in Tahe Oilfield
[J].
