2. 大庆油田有限责任公司勘探开发研究院,大庆 163712
2. Exploration and Development Research Institute of Daqing Oilfield Company Ltd, Daqing 163712, China
徐家围子断陷是松辽盆地北部发育的30多个断陷盆地中规模较大的重要断陷盆地之一[1],2002年以来,徐家围子断陷火山岩勘探取得重大突破,在营城组发现了厚度深,面积广的火山岩气藏,具有良好的规模储量前景.该区营城组火山岩储层以裂缝-孔隙型储层[2]为主,火山岩内气孔、溶蚀孔洞、裂隙发育.火山岩中的裂缝除了本身具有一定的储集能力外,最主要的作用是能够改善储层孔、渗性能.孔、洞常常由裂缝贯通形成储集层,既可增加有效孔隙空间,也可连通孔隙作为油气运移和向井孔渗流的通道,对于气藏的勘探开发具有重要的意义[3].钻井岩芯及测井裂缝研究成果表明,火山岩裂缝形态具有多样性和多尺度性,不仅有构造缝、成岩缝而且有风化缝、溶蚀缝,在形态上表现为高角度缝、水平缝、斜交缝、网状缝[4, 5]等.
火山岩裂缝测井识别与评价研究成果较多[6~9],在气田开发中起到重要的作用,但还是无法评价未钻井区火山岩裂缝的发育情况.应用地震资料预测裂缝的井间分布规律具有实际应用意义.目前已经发展起来的裂缝性储层勘探技术有:S 波勘探、P-S转换波、多分量地震、多方位VSP(verticalseismicprofile)、纵波AVO 等[10~18],其中最有效的方法应属S波分裂技术.但S波受地表影响严重,信噪比和频率较低.多分量地震、多方位VSP、P-S 转换波技术有不错的效果,但要么勘探成本高,要么是非常规地震采集项目,在许多地区难以广泛应用.在各向异性介质中,P 波的AVO 梯度在平行于裂缝走向和垂直于裂缝走向上存在较大差异.这是用P 波检测裂缝的理论基础,与其他技术比较,P 波更低的采集成本和相对成熟的处理技术,使利用纵波来进行裂缝探测更具实用性.文献[19]中P 波方位AVO 对裂缝响应的研究,引发了利用P 波资料检测裂缝和各向异性的热潮[20~26].实践证明,纵波和横波的各向异性响应有很好的相关性[27].
由于火山岩具有岩性、岩相横纵向变化大,非均质性强,其岩石的矿物成分复杂,裂缝、溶蚀孔的类型、组合分布均表现出极强的各向异性等特点[28],使国内外针对火山岩复杂岩性的裂缝研究发展较晚,基本没有成功案例可供借鉴.所有针对火山岩裂缝裂缝的特性,利用叠前地震数据预测火山岩裂缝发育状况是一个值得探索的方法.
2 裂缝预测可行性分析 2.1 叠前地震预测裂缝理论基础地震叠前裂缝预测技术主要是基于地下介质各向异性理论.通常,可用5 个弹性参数(c11,c13,c33,c44 及c66)来描述横向各向同性介质,但是这些弹性参数没有明确的物理意义.Thomsen[29]在弱各向异性的前提下提出5个新的弹性参数来描述横向各向同性介质,这5 个弹性参数分别是:垂向P 波速度α、垂向S 波速度β、表示P 波各向异性程度的ε(V)、表示P波在横向和垂向之间各向异性变化快慢程度的δ(V) 及表示快、慢横波速度差异程度的γ.弹性参数和Thomsen 参数之间的关系为
(1) |
(2) |
(3) |
(4) |
(5) |
式中,c55 为弹性参数,且c55 =c66;ρ为介质的密度.
Rüger[30]在弱各向异性的概念下,推导出各向异性介质中反射系数与弹性参数之间的解析关系,当入射角较小时,近似表示为
其中,
(6) |
式中,Rpp(θ,φ)为与入射角和方位角有关的纵波反射系数;Z为波阻抗,且Z=ρα;G为切向模量,且G=ρβ2,Δ 为上、下介质岩性参数之差;上标"-"为上、下介质岩性参数的平均值;α 为纵波速度;β 为横波速度;γ,δ,ε 为3 个描述介质各向异性的参数;θ为地震波入射的角度;φsym 为裂隙走向的方位角;φ为激发点到接收点的方位角,P为纵波垂直入射时的反射振幅;Giso为振幅随偏移距的变化率(各向同性梯度),指示介质各向同性部分的影响;Gani为振幅随方位角的变化率(各向异性梯度),指示介质各向异性部分的影响.在固定入射角的情况下,式(6)中P波反射振幅响应可以进一步简化为
(7) |
利用式(7)可以方便地反演出A,B,φsym三个与裂隙有关的参数,且φsym指出了裂隙的走向方向,而B/A则可作为裂隙密度的相对度量.
2.2 裂缝地震响应正演模拟通过正演模型的研究,不仅可以帮助我们建立储层岩石物性的理论模型以及理论上该模型所产生的地震各向异性响应特征,而且有助于分析叠前地震资料预测裂缝的可行性以及方位角的叠加处理.在叠前正演模拟研究中,通过徐家围子断陷安达地区钻井的纵横波速度和岩石密度资料,在给定裂缝密度、所含流体、裂缝宽度/长度比等裂缝参数条件下,建立起与之对应的岩石弹性模量之间的关系,从而构建地震各向异性岩石物性模型,最终实现裂缝储层的地质模型及弹性介质模拟.图 1为DS8井进行裂缝正演模拟,通过试气资料获知其主要产气层深度范围为3064~3074m,位于火山岩上部.因此,我们对其给定裂缝参数,根据岩石物理模型计算井中的含有裂缝的火山岩段在含气情况下的岩石的弹性张量和各向异性等效指数,从而了解裂缝对该井产气层各向异性的岩石物理参数的影响.通过得到井中裂缝的各向异性参数,计算裂缝在井旁地震道的地震响应,如叠前地震反射在各个方位角的响应以及叠前AVO 特征,从而通过其各向异性实现地震反射振幅与裂缝的定向和定量描述.
通过对DS8井的裂缝地震响应正演模拟,表明充填流体的定向高角度开启裂缝发育带,会引起地震波出现异常,这些异常主要表现为地震波的振幅、频率、相位等特征发生变化:① 当入射角较小时,相当于地震波的垂直入射,各方位的反射振幅没有差异.图 1b 中入射角达到12°以上,各方位的反射振幅存在差异,入射角越大,差异越明显,这与理论公式一致.由此,叠前裂缝检测需要尽量利用大入射角即大偏移距的资料,以便更能观察到方位角的变化情况;②在固定方位角情况下,该井目的层的反射振幅随入射角的增大,反射振幅变小;③在固定入射角(入射角需大于0°)情况下,垂直于裂缝方向的振幅大于平行于裂缝方向的振幅,其他方位的振幅介于二者之间.方位角为90°(即正北方向)的振幅是椭圆的短轴,与裂缝方向一致.图 1中椭圆短轴代表裂缝方向;④入射角越大,各方位反射振幅差异越大,因此,在数据允许的情况下,适量抛弃一定的近偏移距数据,使得方位异性综合效应明显,有利于提高裂缝检测的精度.
3 火山岩裂缝叠前地震预测 3.1 叠前地震分方位角处理针对裂缝预测的地震资料处理,主要是在叠前常规处理的基础上,以正演模型的指导,进行偏移距部分截除和叠前分方位处理.
叠前地震裂缝检测主要应用地震数据的方位AVO 特性.为了凸显该特性,增强地震资料对裂缝的敏感度,在进行分方位处理前,需要对偏移距范围限制范围在200~2000m 之间,这样处理主要是基于以下两个因素的考虑:① 近偏移距数据来自于近垂直入射地震波,而该反射不能反映裂缝介质的方位异性,这样近偏移距数据参与叠加会削弱远偏移距数据携带的裂缝信息,需抛弃一定的近偏移距数据.但抛弃过多的近偏移距数据,会降低地震资料的覆盖次数和信噪比,所以在舍弃部分近偏移距数据时,需要平衡与覆盖次数的关系.②由于远偏移距数据方位角不全,在大于2000m 以后,几乎都集中在东西方向上,即方位角在90°附近,因此,为了与0°方位角数据匹配,方位角数据中也必须抛弃部分远偏移距数据.
本研究区观测系统的方位角-偏移距分布不完全均匀,在一定程度上限制了地震数据的方位处理,需要进行不等方位角划分(表 1).由于原始地震资料覆盖次数比较小,因此,在保证裂缝预测精度的基础上,适量减小方位角组数以保证每组地震资料拥有足够的覆盖次数,提高地震资料的信噪比和稳定性.
根据各向异性理论,纵波反射系数(绝对值)在平行于裂隙方向达到极大值,在垂直于裂隙方向出现极小值.裂隙的倾角会影响纵波反射系数的大小,倾角越小,纵波方位AVO 差别越小;而且裂隙密度的变化影响着纵波反射系数的大小.所以,定向裂缝带会导致纵波能量衰减和波场干涉等现象,火山岩中高角度裂缝发育带也势必会引起反射延时、强度降低、频率降低等现象,而且裂缝密度越大,现象越明显.
Rüger[31]通过不同的模型研究认为,对于较大的固定入射角,HTI介质的P 波反射系数随入射方位呈现为椭圆,P波反射系数总是在HTI介质的对称轴方向、或者排列片的延伸方向即裂缝方向最大或最小.椭圆的长、短轴大小之比即扁率始终与裂缝密度呈正比,裂缝密度越大,椭圆愈扁.椭圆的长轴还是短轴方位代表裂缝方位,这与实际地层的各种参数有关,可以通过井的模型正演实现.试验证明AVA 技术对裂缝发育的方向比较敏感,本文应用AVA 技术来预测裂缝的方向.
裂缝储层的地震散射理论研究表明[32, 33]:地震频率衰减和裂缝密度场的空间变化有关,沿裂缝走向传播时衰减慢,而垂直裂缝方向衰减快.因此可以采用多信号频率估算技术和基于小波变换的瞬时参数估算技术估算具有物理意义的地震属性参数,包括方位角的功率谱和方位角的瞬时频率随偏移距的变化,借此来确定裂缝的发育密度,裂缝发育密度亦是利用方向椭圆来决定,与裂缝发育方向的确定方法所不同的是始终采用椭圆的长轴与短轴之比来作为裂缝发育密度的度量,所以其最小值为1,表示无各向异性,即无裂缝发育.通过检测结果,本文综合应用AVA 技术和FVA 技术预测裂缝发育的方向和密度.
AVA 与FVA 技术就是在分方位角叠加数据体上检测波阻抗属性与频率衰减属性之间的各向异性,对提取的各向异性值拟合椭圆,求取裂缝发育的方向和密度.
图 2为综合利用各方位角的频率属性与振幅属性对研究区火山岩进行裂缝密度与方向的预测结果.结果表明,裂缝发育带主要集中在DS2-WS101连线一带以及DS8井周围,基本位于工区中部.DS2井与DS8 井裂缝发育方向为北北西-南南东向,WS101发育北北东-南南西向和北西-南东向两组裂缝.
根据裂缝密度预测连井剖面(图 3),裂缝发育带呈区域性串珠状、块状分布.工业气井WS101井附近裂缝较为发育,低产气井SS2附近裂缝发育程度较低,说明裂缝预测的发育状况与单井产能基本吻合.
FMI测井解释结果(表 2与图 4)显示DS2井与WS101井裂缝发育程度较高;DS2 发育一组北北东-南南西向裂缝,WS101 发育北北东-南南西向、北西-南东向和北东-南西向三组裂缝.对比预测结果,可以得到结论:裂缝预测结果与FMI测井解释结果基本一致.
根据图 2,DS2 井、DS8 井、WS101 井区裂缝较发育,SS2 井区裂缝发育程度较弱,结合表 3,DS2与WS101为工业气井,DS8井录井显示大部分为水层,SS2井录井显示为差气层和干层,说明预测结果和单井产能资料基本吻合.
图 5为地震叠前裂缝预测结果与FMI测井裂缝解释结果交汇图,从图上可以看出:
(1) 地震叠前裂缝预测显示,其强度在1.2(各向异性椭圆扁率,无量纲)以上时,裂缝较为发育,储层多为工业气层或是水层;位于该值以下时,裂缝不发育,储层多为干层或低产气层.裂缝预测结果与单井产能情况基本一致.
(2) 地震叠前裂缝密度预测与FMI测井裂缝具有一定关系,随着各向异性值的增大,成像测井显示其裂缝发育密度也呈增大趋势.地震叠前裂缝预测结果与钻井FMI裂缝发育基本一致.
4.2 裂缝发育地质因素分析在平面上,根据裂缝预测结果,工区中部裂缝比较发育,裂缝发育程度较高的地区是局部高点所在(如图 6),说明裂缝的产生与古火山活动和构造运动密切相关.
在图 7b剖面上(暖色代表裂缝相对发育区,冷色调代表相对不发育区),可以看出裂缝在剖面上并不是杂乱无章的排列,而是以区带性的方式呈珠状和块状分布,这是由于火山岩储层岩性变化较大引起的.
通过图 7 可以看出,在WS101 井与DS8 井区目的层裂缝发育程度较高,WS101 的单井产能为222538m3/d,DS8井差气层以下大段含水,分析其原因:WS101井正处在断层的正上方,存在构造运动形成的构造缝,而且本段火山岩大部分属于喷溢相上部亚相,原生孔隙比较发育.DS8井的构造位置位于火山机构顶部,存在流纹岩的原生孔隙和构造运动形成的构造缝[34, 35].综上所述,在剖面上预测结果符合地质规律,与钻井资料一致性很好.
综上所述,裂缝预测密度在平面和剖面上与钻井资料吻合度较高,裂缝预测的方向与钻井FMI测井解释结果一致性较好,符合地质规律,表明利用地震叠前数据预测火山岩裂缝效果较好,对裂缝主导型火山岩油气藏的勘探开发具有一定的指导作用.
5 结论与建议(1) 地震叠前裂缝预测技术实现过程无须井约束,人为干预少,结果相对客观.在安达地区裂缝预测中,预测结果与已钻井资料基本吻合,说明该方法在深层火山岩储层裂缝预测方面的可靠性,可为火山岩油气藏的勘探开发提供有效的技术支撑.
(2) 应用叠前地震数据开展深层火山岩储层裂缝预测,具有横向分辨率高、实现井间分布描述和宏观空间分布规律清楚等特点.火山岩裂缝发育带主要集中在火山机构主体构造高部位,裂缝的产生与火山活动以及构造运动密切相关.
(3) 应用地震叠前数据进行裂缝预测,最好在野外地震资料采集时,进行宽方位角及高覆盖次数采集;在处理过程中,尽最大能力做到相对保幅、保真处理,以及准确精细的地震构造解释,从而保证预测结果的可靠性.
[1] | 任延广, 朱德丰, 万传彪, 等. 松辽盆地北部深层地质特征与天然气勘探方向. 中国石油勘探 , 2004, 9(4): 12–19. Ren Y G, Zhu D F, Wan C B, et al. Geological characteristics of deep layers in northern part of Songliao Basin and orientation for natural gas exploration. China Petroleum Exploration (in Chinese) , 2004, 9(4): 12-19. |
[2] | 邵红梅, 毛庆云, 姜洪启, 等. 徐家围子断陷营城组火山岩气藏储层特征. 天然气工业 , 2006, 26(6): 29–32. Shao H M, Mao Q Y, Jiang H Q, et al. Reservoir characteristics of volcanic rock gas pool in Yingcheng formation in Xujiaweizi fault depression. Atural Gas Industry (in Chinese) , 2006, 26(6): 29-32. |
[3] | 杨辉, 张研, 邹才能, 等. 松辽盆地北部徐家围子断陷火山岩分布及天然气富集规律. 地球物理学报 , 2006, 4(49): 1136–1143. Yang H, Zhang Y, Zou C N, et al. Volcanic rock distribution and gas abundance regularity in Xujiaweizi faulted depression, Songliao basin. Chinese J. Geophys. (in Chinese) , 2006, 4(49): 1136-1143. |
[4] | 舒萍, 纪学雁, 丁日新, 等. 徐深气田火山岩储层的裂缝特征研究. 大庆石油地质与开发 , 2008, 2(27): 13–17. Shu P, Ji X Y, Ding R X, et al. Study on fracture features of volcanic reservoirs in Xushen gas field. Petroleum Geology & Oilfield Development in Daqing (in Chinese) , 2008, 2(27): 13-17. |
[5] | 姜传金, 陈树民, 初丽兰, 等. 徐家围子断陷营城组火山岩分布特征及火山喷发机制的新认识. 岩石学报 , 2010, 26(1): 63–72. Jiang C J, Chen S M, Chu L L, et al. A New Understanding about the Volcanic distribution characteristics and eruption mechanism of Yingchen formation in Xujiaweizi fault depression. Acta Petrologica Sinica (in Chinese) , 2010, 26(1): 63-72. |
[6] | 邓攀, 陈孟晋, 高哲荣, 等. 火山岩储层构造裂缝的测井识别及解释. 石油学报 , 2002, 23(6): 32–36. Deng J, Chen M J, Gao Zh R, et al. Log response and explanation of structural fractures in volcanic rock reservoir. Acta Petrolei Sinica (in Chinese) , 2002, 23(6): 32-36. |
[7] | 李佳阳, 夏宁, 秦启荣. 成像测井评价致密碎屑岩储层的裂缝与含气性. 测井技术 , 2007, 31(1): 17–20. Li J Y, Xia N, Qin Q R. Fracture and gas-bearing evaluation of tight clastic reservoir with image logging. Well Logging Technology (in Chinese) , 2007, 31(1): 17-20. |
[8] | 张树东, 司马立强, 刘萍英, 等. 基于测井新技术解释有效裂缝发育规律. 西南石油大学学报 , 2007, 29(1): 23–25. Zhang S D, Si Ma L Q, Liu P Y, et al. Interpret the law of effective fracture development based on new technology in well logging. Journal of Southwest Petroleum University (in Chinese) , 2007, 29(1): 23-25. |
[9] | 史謌, 何涛, 仵岳奇, 等. 用正演数值计算方法开展双侧向测井对裂缝的响应研究. 地球物理学报 , 2004, 47(2): 359–363. Shi G, He T, Wu Y Q, et al. A study on the dual laterolog response to fractures using the forward Numerical modeling. Chinese J. Geophys. (in Chinese) , 2004, 47(2): 359-363. |
[10] | 唐晓雪, 唐建候. 用P-SV 转换波资料研究储层各向异性. 石油地球物理勘探 , 1997, 32(1): 27–33. Tang X X, Tang J H. Study reservoiranisotropy by P-SV conversion data. Oil Geophysical Prospect (in Chinese) , 1997, 32(1): 27-33. |
[11] | Carcuz J R. A combined AVO analysis of P-P and P-S reflection data. In:Expanded Abstract of the 71st Annual International SEG Meeting , 2001. |
[12] | 朱培民, 王家映, 於文辉, 等. 用纵波AVO数据反演储层裂隙密度参数. 石油物探 , 2001, 40(2): 1–12. Zhu P M, Wang J Y, Yu W H, et al. Invert reservoir fracture-density parameters by P-wave AVO data. Geophysical Prospecting for Petroleum (in Chinese) , 2001, 40(2): 1-12. |
[13] | 曲寿利, 季玉新, 王鑫, 等. 全方位P波属性裂缝检测方法. 石油地球物理勘探 , 2001, 36(4): 390–397. Qu S L, Ji Y X, Wang X, et al. Seismic method for using full-azimuth P-wave attribution to detect fracture. Oil Geophysical Prospecting (in Chinese) , 2001, 36(4): 390-397. |
[14] | 李绪宣. 多波地震地层各向异性的初步研究和探讨. 天然气工业 , 2001, 21(5): 41–44. Li X X. An initiative study and discussion on rocks' anisotropy by using multiwave seismics. Natural Gas Industry (in Chinese) , 2001, 21(5): 41-44. |
[15] | 杜启振, 杨慧珠. 方位各向异性介质的裂缝预测方法研究. 石油大学学报(自然科学版) , 2003, 27(4): 32–36. Du Q Z, Yang H Z. Detection method for fractures in azimuthally anisotropic media. Journal of the University of Petroleum(Natural Science) (in Chinese) , 2003, 27(4): 32-36. |
[16] | 罗省贤, 李录明. 基于横波分裂的地层裂缝预测方法与应用. 成都理工大学学报(自然科学版) , 2003, 30(1): 52–59. Luo S X, Li L M. Method and application of predicting medium fracture based on shear-wave separating theory. Journal of Chengdu University of Technology(Natural Science) (in Chinese) , 2003, 30(1): 52-59. |
[17] | 王建民, 杨冬, 魏修成, 等. 多分量地震资料预测松辽盆地兴城地区深层火山岩与有利含气带. 地球物理学报 , 2007, 6(50): 1914–1923. Wang J M, Yang D, Wei X C, et al. Predicting deep volcanic rocks and favorable gas zone near Xingcheng area in Songliao Basin using multi-component seismic data. Chinese J.Geophys. (in Chinese) , 2007, 6(50): 1914-1923. |
[18] | 孙炜, 王彦春, 李梅, 等. 利用叠前地震数据预测火山岩储层裂缝. 物探与化探 , 2010, 34(2): 229–232. Sun W, Wang Y C, Li M, et al. The Detection of Fractures in the Volcanic Reservoir with Pre-Stack Seismic Data. Geophysical &Geochemical Exploration (in Chinese) , 2010, 34(2): 229-232. |
[19] | Mallick S, Frazer L N. Reflection transmission coefficients and azimuthal anisotropy in marine studies. Geophys.J. Int. , 1991, 105(1): 241-252. |
[20] | Johnson W E. Direct detection of gas in pre-Tertiary sediments. The Leading Edge , 1995, 14(2): 119-122. DOI:10.1190/1.1437100 |
[21] | Lynn H B, Simon K M. Azimuthal anisotropy in P-wave 3D multiazimuth data. The Leading Edge , 1996, 15(8): 923-928. DOI:10.1190/1.1437392 |
[22] | Mallick S, Craft K L, Meister L J, et al. Determination of the principal directions of azimuthal anisotropy from P-wave seismic data. Geophysics , 1998, 63(2): 692-706. DOI:10.1190/1.1444369 |
[23] | 阴可, 杨慧珠. 各向异性介质中的AVO. 地球物理学报 , 1998, 41(3): 382–391. Yin K, Yang H Z. AVO in anisotropic media. Chinese J. Geophys. (in Chinese) , 1998, 41(3): 382-391. |
[24] | 刘洋, 董敏煜. 各向异性介质中的方位AVO. 石油地球物理勘探 , 1999, 34(3): 260–267. Liu Y, Dong M Y. Azimuthal AVO in anisotropic medium. Oil Geophysical Prospecting (in Chinese) , 1999, 34(3): 260-267. |
[25] | Shen F, Sierra J, Burns D R, et al. Azimuthal offset-dependent attributes applied to fracture detection in a carbonate reservoir. Geophysics , 2002, 67(2): 355-364. DOI:10.1190/1.1468596 |
[26] | Hall S A, Kendall M J. Fracture characterization at Valhall:application of P-wave amplitude variation with offset and azimuth(AVOA) analysis to a 3D ocean-bottom data set. Geophysics , 2003, 68(4): 1150-1160. DOI:10.1190/1.1598107 |
[27] | Lynn H B, Simon K M, Bates C R. Correlation between P-wave AVOA and S-wave traveltime anisotropy in naturally fractured gas reservoir. The Leading Edge , 1996, 15(8): 931-935. DOI:10.1190/1.1437394 |
[28] | 张丽华, 潘保芝, 单刚义, 等. 长岭地区火山岩储层流体性质测井预测. 地球物理学进展 , 2009, 24(6): 2151–2155. Zhang L H, Pan B Z, Shan G Y, et al. Fluid property logging prediction of volcanic reservoirs in the Changling area. Progress in Geophysics (in Chinese) , 2009, 24(6): 2151-2155. |
[29] | Thomsen L. Weak elastic anisotropy. Geophysics , 1986, 51(10): 1954-1966. DOI:10.1190/1.1442051 |
[30] | Rüger A. P-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry. Geophysics , 1997, 62(3): 713-721. DOI:10.1190/1.1444181 |
[31] | Rüger A. Variation of P-wave reflectivity with offset and azimuth in anisotropic media. Geophysics , 1998, 63(3): 935-947. DOI:10.1190/1.1444405 |
[32] | Hudson J A. Wave speeds and attenuation of elastic waves in material containing cracks. Geophys.J.R.Astr.Soc. , 1981, 64: 133-150. |
[33] | Hudson J A. Overall elastic properties of isotropic materials with arbitrary distribution of circular cracks. Geophys.J.Int., , 1990, 102: 465-469. |
[34] | 王璞珺, 迟元林, 刘万洙, 等. 松辽盆地火山岩相:类型、特征和储层意义. 吉林大学学报(地球科学版) , 2003, 33(4): 449–456. Wang P J, Chi Y L, Liu W Z, et al. Volcanic facies of the Songliao basin:Classification,characteristics and reservoir significance. Journal of Jilin University (Earth Science Edition) (in Chinese) , 2003, 33(4): 449-456. |
[35] | 陈建文, 王德发, 张晓东, 等. 松辽盆地徐家围子断陷营城组火山岩相和火山机构分析. 地学前缘 , 2000, 3(4): 371–379. Chen J W, Wang D F, Zhang X S, et al. Analysis of volcanic facies and apparatus of Yingcheng formation in Xujiaweizi faulting depression,Songliao basin northeast China. Earth Science Frontiers (in Chinese) , 2000, 3(4): 371-379. |