利用_logR技术计算柴达木盆地三湖地区第四系有机碳
- 格式:pdf
- 大小:962.45 KB
- 文档页数:8
柴达木盆地三湖地区第四系生物气藏水动力成因研究的开题报告摘要:柴达木盆地是我国重要的石油、天然气生产基地之一,其中三湖地区是盆地内石油、天然气勘探开发的重要区域。
本文选取三湖地区为研究对象,针对该地区第四系生物气藏水动力成因进行综合分析,包括地质背景、孔隙、渗透率、油藏水力压力、流体性质等方面,探讨了生物气藏的形成机制和水动力成因。
研究表明,三湖地区第四系生物气藏的水动力成因主要是由于沉积环境的影响,而孔隙和流体性质等因素也对水动力成因起到了一定的影响。
关键词:柴达木盆地;三湖地区;第四系生物气藏;水动力成因;地质背景一、研究背景柴达木盆地是我国重要的石油、天然气生产基地之一,其中三湖地区是盆地内石油、天然气勘探开发的重要区域。
该地区第四系生物气藏具有丰富的储量和开发价值,因此对于该地区的生物气藏的研究具有重要的意义。
二、研究目的本文旨在对三湖地区第四系生物气藏的水动力成因进行综合分析,探讨生物气藏的形成机制和水动力成因,为该地区的生物气藏开发提供科学的参考。
三、研究内容本文将从以下几个方面开展研究:(1)研究地质背景通过对三湖地区的地质情况进行调查和分析,了解该地区的地质背景,包括沉积环境、地质构造、岩性等方面的情况。
(2)研究孔隙、渗透率等方面对三湖地区生物气藏的孔隙、渗透率等方面的情况进行调查和分析,探讨不同孔隙、渗透率对水动力成因的影响。
(3)研究油藏水力压力通过现场调查和试验,了解三湖地区第四系生物气藏的油藏水力压力情况,分析水力压力对生物气藏的形成和演化的影响。
(4)研究流体性质对三湖地区第四系生物气藏的流体性质进行分析和研究,探讨不同流体性质对水动力成因的影响。
四、研究方法本文将采用现场调查、分析实验等方法进行研究,具体包括:(1)采集三湖地区第四系生物气藏的样本采集三湖地区第四系生物气藏的样本,进行物化性质分析和实验测试,以了解其孔隙、渗透率、流体性质等情况,为后续的水动力成因分析提供数据支撑。
山西农经/2019年5期DOI:10.16675/14-1065/f.2019.05.058柴达木盆地东部第四纪水文地质条件对当地农业的影响□徐汝汝摘要:柴达木盆地东部第四纪水文地质条件一直是国内学术界的研究热点。
该地有地层新、规模大、成因特殊等特点,从资源的角度出发,该地区的植物资源以及其他相关资源的勘探受到了地质学家的广泛关注。
天然气的运输和地层水的关系极为密切,而植物生存的基本条件就是有充足的地层水,且地层水间接关系到农业的发展,因此选择这一领域进行了研究。
论述了柴达木盆地东部地区第四纪水文地质条件对当地农业的影响,为有关部门提供参考。
关键词:柴达木盆地;东部;水文地质条件;当地农业文章编号:1004-7026(2019)05-0098-02中国图书分类号:S162.5文献标志码:A(山东师范大学地理与环境学院山东济南250358)柴达木盆地东部地区第四纪沉积以及沉降中心在构造上包括三级构造单元,分别是盐湖斜坡、三湖凹陷和格尔木斜坡。
在这一区域中,地面特点显著,盐湖分布广,地下水的水位高、盐度高、温度低,且出现了矿化现象[1],地理环境条件并不适合植物生存,但是却拥有丰富的植物气藏,存在一定的特殊性。
所以从水文地质条件入手,对于该地区的农业发展状况进行了相应的研究。
1柴达木盆地东部第四纪水文与地质发展状况1.1地质状况柴达木盆地东部地区第四纪沉降和沉积中心属凹陷区,凹陷中心位于西台吉乃尔湖南侧,整个凹陷带呈现出展布型特点[2],具体有鸭湖、南丘陵、驼峰山、盐湖等。
从断裂层面来看,呈现出逆断层的趋势,背斜顶部存在少数发育规模较小的正断层。
1.2水文状况从水系发育特点来讲,柴达木盆地气候干旱寒冷,盆地平均年降水量较低但蒸发量较高。
这一地区仍有湖泊存在,水源充足,这些水源主要为降水、降雪以及周边高山融化的雪水。
柴达木盆地东部周边山区有完善的供水系统[3],大致以无柴沟—大灶火为界限,西区主要依靠昆仑山水系供水,东区依靠于南北两侧山脉同时供水。
柴达木盆地第四系生物气的勘探历程与储量现状徐子远;谢丽;张道伟;姜桂凤;韩凤祥【摘要】以储量研究为主线,以地质认识为重点,通过对不同时期研究认识和勘探思路的阐述,展现了勘探对象由简到繁、从地面构造走向潜伏构造,地质认识由浅到深、从分析现象走向追求本质,储量研究由粗到精、从定性认识走向定量评价,储量规模由小到大、从小型气田走向大型气区的发展过程,对类似沉积盆地的油气勘探和储量研究工作具有积极的借鉴和指导意义.【期刊名称】《新疆石油地质》【年(卷),期】2005(026)004【总页数】4页(P437-440)【关键词】柴达木盆地;第四系;生物气;勘探历程;储量现状【作者】徐子远;谢丽;张道伟;姜桂凤;韩凤祥【作者单位】中国石油,青海油田分公司,勘探开发研究院,甘肃,敦煌,736202;中国石油,青海油田分公司,勘探开发研究院,甘肃,敦煌,736202;中国石油,青海油田分公司,勘探开发研究院,甘肃,敦煌,736202;中国石油,青海油田分公司,勘探开发研究院,甘肃,敦煌,736202;中国石油,青海油田分公司,勘探开发研究院,甘肃,敦煌,736202【正文语种】中文【中图分类】TE122.11柴达木盆地位于我国青藏高原北部,是一个为祁连山、阿尔金山和昆仑山所环绕的中、新生代大型内陆盆地,地理面积25×104km2,沉积岩分布面积12.1×104km2,最大沉积厚度1.72×104m.是我国七大内陆含油气盆地之一,油气总资源量超过40×108t当量,其中天然气资源24900×108m3,大致可分为西部古近-新近系成熟油型气、北缘侏罗系高成熟煤型气和东部第四系未成熟生物气(图1)。
目前盆地勘探发现的天然气储量,主要集中于盆地东部三湖地区的第四系[1]。
根据沉积相及成藏规律研究,三湖地区的第四系与新近系为连续沉积,已知气藏均为第四系和新近系自生自储的动态生物气藏。
第16卷 第1期2021年1月中国科技论文C H I N AS C I E N C E P A P E RV o l .16N o .1J a n .2021基于改进的Δl o g R 技术的鄂尔多斯盆地彭阳地区长7烃源岩有机碳含量评价王 龙1,伊 硕2,来又春3,陈国宁2,周立业3(1.中国海洋石油国际有限公司,北京100028;2.中海油研究总院有限责任公司,北京100028;3.中海石油有限公司天津分公司,天津300459)摘 要:鄂尔多斯盆地西南缘彭阳地区发育三叠系延长组长7烃源岩,受岩心测试样品有限的约束,仅通过样品实测有机碳含量评估烃源岩有机质丰度的方法已经无法满足致密油精细勘探的需求㊂针对研究区测井资料较为丰富的特点,结合实测数据与测井曲线,利用改进的Δl o g R 技术对研究区烃源岩进行评价㊂结果表明,改进的Δl o g R 技术相对于传统的Δl o g R 技术减少了人为误差与工作量,同时减少了系统误差,能较好地反应烃源岩总有机碳(t o t a l o r ga n i c c a rb o n ,T O C )含量㊂根据长7内烃源岩性质的差异,分别建立长7高阻泥岩和长7普通泥岩的T O C 含量评价模型对长7烃源岩进行评价,结果表明:改进的Δl o g R 技术与实测T O C 含量平均误差不超过11%,结果较为可靠;彭阳地区长7高阻泥岩T O C 含量整体较高,除西北部外,大部分地区都能达到优质烃源岩的级别;长7普通泥岩T O C 含量整体较低,大部分地区也能达到中等偏好的烃源岩级别㊂关键词:总有机碳含量;测井评价;改进的Δl o g R 技术;彭阳地区;鄂尔多斯盆地中图分类号:T E 122 文献标志码:A文章编号:20952783(2021)01004408开放科学(资源服务)标识码(O S I D ):E v a l u a t i o n o f T O C c o n t e n t o f C h a n g 7s o u r c e r o c k s i nP e n g y a n g ar e a o f O r d o s b a s i n b a s e d o n i m p r o v e dΔl o g R t e c h n o l o g yW A N GL o n g 1,Y I S h u o 2,L A IY o u c h u n 3,C H E NG u o n i n g 2,Z H O UL i ye 3(1.C N O O C I n t e r n a t i o n a l L t d .,B e i j i n g 100028,C h i n a ;2.C N O O CR e s e a r c h I n s t i t u d e ,B e i j i n g 100028,C h i n a ;3.T i a n j i n B r a n c h o f C N O O CC o .,L t d .,T i a n ji n 300459,C h i n a )A b s t r a c t :C h a n g 7s o u r c e r o c k s o f T r i a s s i cY a n c h a n g f o r m a t i o nw e r e d e v e l o p e d i nP e n g y a n g ar e a o f s o u t h -w e s tO r d o s b a s i n ,b e -c a u s e o f t h e c o n s t r a i n o f l i m i t e d c o r e t e s t s a m p l e s ,t h em e t h o d o f e v a l u a t i n g t h e o r g a n i cm a t t e r a b u n d a n c e o f s o u r c e r o c k o n l y b yt h em e a s u r e d t o t a l o r g a n i c c a r b o n (T O C )c o n t e n t o f t h e s a m p l e c a n n o l o n g e rm e e t t h e d e m a n d o f f i n e t i g h t o i l e x pl o r a t i o n .T h e l o g g i n g d a t a i n t h e s t u d y a r e a i s r e l a t i v e l y a b u n d a n t ,t h u s t h e i m p r o v e dΔl o g R t e c h n o l o g y wa s u s e d t o e v a l u a t e t h e s o u r c e r o c k sb yc o m b i n i n g w i t h t h em e a s u r ed d a t a a n dwe l l l o g g i n g c u r v e .C o m p a r e dw i t h t h e t r a d i t i o n a l Δl o g R t e c h n o l o g y,t h i sm e t h o d c a n r e d u c e t h ew o r k l o a d ,h u m a n e r r o r a n d t h e s y s t e me r r o r ,t h e T O C c o n t e n t o f s o u r c e r o c k c a n b ew e l l e v a l u a t e d t o o .A c c o r d i n gt o t h e d i f f e r e n c e o f s o u r c e r o c k s i nC h a n g 7,T O C c o n t e n t e v a l u a t i o nm o d e l o f C h a n g 7h i g h r e s i s t a n c em u d s t o n e a n dC h a n g 7or d i -n a r y m u d s t o n ew e r e e s t a b l i s h e d r e s p e c t i v e l y .R e s u l t s s h o w t h a t a v e r a g e e r r o r b e t w e e n i m p r o v e dΔl o g R t e c h n o l o g y a n d t h em e a s -u r e dT O C c o n t e n t i s l e s s t h a n 11%,t h e r e s u l t i sm o r e r e l i a b l e .T h e e v a l u a t i o n r e s u l t s h o w t h a t t h e C h a n g 7h i g h r e s i s t a n c em u d -s t o n e i nP e n g y a n g a r e a h a s h i g hT O Cc o n t e n t a n db e l o n g t oh i g h -q u a l i t y s o u r c e r o c k s i nm o s t a r e a e x c e pt t h en o r t h w e s t p a r t .C h a n g 7o r d i n a r y m u d s t o n e h a s r e l a t i v e l y l o wT O C c o n t e n t ,a n d b e l o n g t om i d d l e t o g o o d s o u r c e r o c k s i nm o s t p a r t o f P e n g y a n ga r e a .K e yw o r d s :t o t a l o r g a n i c c a r b o n (T O C )c o n t e n t ;l o g g i n g e v a l u a t i o n ;i m p r o v e dΔl o g R t e c h n i q u e ;P e n g y a n g a r e a ;O r d o s b a s i n 收稿日期:2020-03-16基金项目:国家科技重大专项(2017Z X 05032-04-01)第一作者:王龙(1986 ),男,高级工程师,主要研究方向为油气田勘探㊁开发地质,w a n g l o n g@c n o o c .c o m .c n 烃源岩评价是油气勘探的重要组成部分,常规烃源岩评价多是从岩心或岩屑入手,分析烃源岩的地球化学特征,但受限于岩心㊁岩屑样品的数量,难以实现烃源岩纵向或平面上连续变化特征的认识㊂由于测井资料在纵向上具有连续性,近些年一些学者开始利用测井资料识别与评价烃源岩,Δl o g R 技术就是其中一项发展较为成熟且应用比较广泛的技术㊂Δl o g R 技术是由Ex x o n /E s s o 公司于1979年开发的一种利用测井资料计算总有机碳(t o t a l o r ga n i c c a rb o n ,T O C )含量的技术[1-2]㊂该方法自2000年前后引入我国后,由于其在烃源岩T O C 含量评价上的便捷性等优势,很快得到了普遍的认可,并得到了广泛的应用㊂张志伟等[3]㊁朱光有等[4]㊁杜江民等[5]应用此方法开展了单井有机碳含量评价;郭泽清等[6]㊁王宗礼等[7]㊁金凤鸣等[8]㊁王振升等[9]利用此技术确定烃源岩的平面分布与规模;卢双舫等[10]用此方法评价了盆地内不同丰度烃源岩的体积㊂在应用该方第1期王龙,等:基于改进的Δl o g R技术的鄂尔多斯盆地彭阳地区长7烃源岩有机碳含量评价法的同时,学者们也根据研究区的实际地质特点对该方法进行了一定的改进与完善:曲彦胜等[11]探讨了应用Δl o g R技术计算T O C含量结果的影响因素;刘超等[12]提出了变系数Δl o g R技术,以解决Δl o g R技术中通过经验参数预测T O C含量所产生的较大误差和人工确定 基线 的不确定性与繁琐性㊂此外,改进的Δl o g R技术经过朱光有等[13]㊁刘超等[14]㊁胡慧婷等[15]的不断修正与应用而逐渐趋向于成熟㊂该方法能够实现自动优选叠合系数,在缺少成熟度参数㊁无需人为确定T O C含量背景值及基线值的情况下,通过计算机快速㊁准确地计算T O C 含量[14]㊂1烃源岩地球化学特征与测井响应彭阳地区位于鄂尔多斯盆地西南缘,构造上属于天环坳陷西翼,主要的勘探层位为三叠系延长组和侏罗系延安组㊂近些年来在延安组中发现了多个油藏,在延长组长3㊁长7㊁长8等多个油层组也发现了油气显示,显示了彭阳地区具有较好的勘探前景㊂鄂尔多斯盆地中生代晚三叠世延长期长7沉积期沉积的1套深湖-半深湖相富含有机质的暗色泥岩,是鄂尔多斯盆地中生界油气主要来源[16]㊂前人[17-20]对长7烃源岩的研究多集中在盆地内部,对位于盆地西南缘彭阳地区长7烃源岩的研究较少,同时由于实验经费㊁样品来源等因素的限制,该地区在长7段的烃源岩取样测试较少,只有30余个有效测试样品,且不论是在平面上还是纵向上的分布都比较离散,难以系统表征区域上长7烃源岩的T O C含量与平面分布;另一方面,该地区具有较为丰富的常规测井资料㊂因此,为了准确评价彭阳地区烃源岩的分布特征,本文采用地化与测井相结合的方法,利用烃源岩实测数据与改进的Δl o g R技术定量评价彭阳地区已钻井的T O C含量,然后依据单井评价结果开展烃源岩特征的平面分布研究,从而实现对鄂尔多斯盆地彭阳地烃源岩的系统评价,以期为该区中生界石油的进一步勘探提供参考依据㊂本研究区鄂尔多斯盆地彭阳地区区域构造位置,如图1[21]所示㊂1.1长7烃源岩地球化学特征岩心㊁地球化学数据分析表明:鄂尔多斯盆地彭阳地区长7高阻泥岩有机质丰度高,T O C含量分布范围为0.81%~29%,平均值为8.01%,氯仿沥青A 的平均含量为0.99%,生烃潜量的平均值为41.76m g/g;长7普通泥岩的T O C含量分布范围为0.21%~8.19%,平均值为0.93%,氯仿沥青 A 平均含量为0.071%,生烃潜量平均值为1.94m g/g,相对于长7高阻泥岩明显偏低㊂长7高阻泥岩的干酪根类型都属于Ⅰ型和Ⅱ1型,长7普通泥岩的干酪根类型多为Ⅱ2型和Ⅲ型;长7高阻泥岩镜质体反射率(R o)实测值为0.77~0.95,长7普通泥岩样品R o 值为0.73~0.91,均已经达到了成熟演化阶段㊂图1鄂尔多斯盆地彭阳地区区域构造位置[21]F i g.1 M a p s h o w i n g r e g i o n a l g e o l o g y l o c a t i o n o fP e n g y a n g a r e a,O r d o s b a s i n[21]总的来说,彭阳地区长7高阻泥岩有机质丰度高,类型好(Ⅰ型和Ⅱ1型),达到了成熟演化阶段,为优质烃源岩;延长组长7普通泥岩有机质丰度相对较低,有机质类型为Ⅱ2型和Ⅲ型,达到了成熟演化阶段,属于好到差的烃源岩㊂1.2长7烃源岩测井响应特征长7湖相暗色泥岩色括黑色泥岩㊁深灰色泥岩㊁灰黑色泥岩㊁灰黑色粉砂质泥岩㊁深灰色粉砂质泥岩等多种岩性,尤以长7黑色泥岩最为重要,其在测井曲线上具有高伽马(G R)㊁高电阻率(R t)㊁高声波时差(A C)和低密度(D E N)的 三高一低 特征,如某典型井长7烃源岩表现为高G R(150~500A P I)㊁高R t(60~580Ω㊃m)㊁高A C(230~370μs/m)和低D E N(1.9~2.5g/c m3)㊂彭阳地区典型井长7烃源岩测井响应特征如图2所示㊂2长7烃源岩测井评价2.1改进的Δl o g R技术原理改进的Δl o g R技术的理论依据是烃源岩层段在A C和R t测井曲线上的特殊响应特征,通过设定的叠合系数将对数坐标下的R t曲线和算术坐标下的A C叠合㊂通常叠合系数的取值为0.02f t/μs,即每个R t对数周期对应A C的刻度范围为50μs/f t,使2条曲线在细粒非烃源岩处重合,此时2条曲线之间的幅度差即为Δl o g R值㊂图3为Δl o g R技术识别地层示意图㊂由图3(a)54中国科技论文第16卷图2彭阳地区典型井长7烃源岩测井响应特征F i g.2L o g g i n g c h a r a c t e r o f C h a n g7s o u r c e r o c k o ft y p i c a l w e l l i nP e n g y a n g a r e a可见叠合后的2条测井曲线会出现以下几种情况: 1)A C和R t均呈现高异常的层段为成熟烃源岩段(图3(a)F),烃类的生成造成R t呈高值;2)A C和R t曲线相重合的层段为非烃源岩段(图3(a)A);3) A C表现为高异常和R t无异常的层段为未成熟的烃源岩段(图3(a)C),因为孔隙流体一般为地层水,所以R t呈低值;4)A C无异常和R t呈现高值的层段为含油储集层段(图3(a)D),由于储集层的孔隙空间中同时存在液态烃和地层水,因此液态烃的R t值高㊂在彭阳地区钻井实例中也证实了Δl o g R方法的有效性(图3(b)):1)长7上部R t曲线与A C曲线重合(图3(b)A),相当于图3(a)中的A段;2)长7中部R t低值㊁A C无异常(图3(b)B),与图3(a)中的水层相对应,解释结论证实该层为水层;3)长7底部高阻段,R t㊁A C均表现为高异常(图3(b)C),属于成熟烃源岩,与图3(a)中F段相对应;4)长8的R t呈高值㊁A C无异常(图3(b)D㊁E),是含油储层,与图3(a)中的D相对应,在解释结论中也证实该层含油㊂由A C㊁R t叠加得到Δl o g R技术的数学模型为Δl o g R=l o g(R/R基线)+0.02(Δt-Δt基线)㊂(1)式中:Δl o g R为2条曲线间的间距;R基线为基线的电阻率值,Ω㊃m;R为仪器实测的电阻率,Ω㊃m;Δt 为测井仪实测的声波时差,μs/f t;Δt基线为基线对应的声波时差,μs/f t;0.02为在对数坐标下的一个电阻率单位和算术坐标下一个声波时差周期50μs/f t的比值㊂Δl o g R与T O C含量具有线性相关关系,一般由Δl o g R计算T O C含量的经验公式为T O C=Δl o g Rˑ102.297-0.1688L O M㊂(2)图3 Δl o g R技术识别地层示意图F i g.3S t r a t i g r a p h i c f e a t u r e i d e n t i f i e d b yΔl o g R t e c h n o l o g y通过式(2)可计算出各井全井段的T O C值,但由以上公式可知,L O M选取不当时,对T O C的绝对含量会产生整体误差,同时,基线值的读取也较繁64第1期王龙,等:基于改进的Δl o g R技术的鄂尔多斯盆地彭阳地区长7烃源岩有机碳含量评价琐,易引入人为误差㊂因此,考虑对该方法模型进行改进㊂现将式(1)修改为Δl o g R=l o g(R/R基线)+K(Δt-Δt基线)㊂(3)其中,K=l o g(R m a x/R m i n)/(Δt m a x-Δt m i n)㊂(4)式中,K为单位声波时差所对应的在对数坐标下电阻率单位的个数㊂若设定对数坐标下的每个电阻率单位对应算术坐标下A C50μs/f t刻度,则K值为0.02[21-23]㊂确定基线后,有Δt基线=Δt m a x-l o g(R基线/R m i n)/K㊂(5)然后将式(4)和式(5)代入式(3),则有Δl o g R=l o g R+l o g(R m a x/R m i n)/(Δt m a x-Δt m i n)ˑ(Δt-Δt m a x)-l o g R m i n㊂(6)由式(6)可以看出:R值和Δt值已知,只需读取Δt m i n㊁Δt m a x(声波时差的刻度范围)和R m i n㊁R m a x(电阻率的刻度范围)便可求取Δl o g R值㊂若固定声波时差的刻度范围,并调整电阻率刻度范围使基线重合,使得K值不变,此时计算出的任意2点间Δl o g R的相对大小不变㊂故可以把式(2)修改为T O C=AˑΔl o g R+B㊂(7)式中:A㊁B分别为拟合公式的系数,其中A与公式(2)中的102.296-0.1688L O M相对应;B为常数,与Δl o g R有关㊂改进后的模型无需再读取A C及R t的基线值,减少了人为误差并大量减少了工作量;其次,无需再设定L O M值,减少引入系统误差,结果能较好地反应烃源岩T O C含量㊂此外,采用V i s u a l S t u d i o进行编程,使其应用性更加广泛,提高了工作效率㊂2.2实际应用延长组长7泥岩是鄂尔多斯盆地中生界最重要的烃源岩,由于长7发育一套高阻泥岩,与长7普通泥岩的测井响应特征差别较大,为了提高评价的精度,分别建立了长7高阻泥岩和长7普通泥岩2个模型(图4)㊂另外在实际处理过程中,首先对涉及取样的岩心段都进行了岩心归位,避免岩心取样点的深度与测井深度不一致;然后为了能够准确区分高阻泥岩和普通泥岩,应在各井的测井曲线上定性判别此井是否发育高阻烃源岩,确定高阻泥岩发育的深度范围;最后由于粉砂岩㊁砂岩等非烃源岩也会有Δl o g R的出现,因此在处理目的层段测井曲线之前,要根据G R㊁S P和C A L等曲线去除非烃源岩层段,确保计算结果的准确性㊂2.2.1长7高阻泥岩T O C预测模型长7高阻泥岩实测T O C样品共有31个㊂依据上述原理建立测井解释模型,首先计算出实测样品对应深度的Δl o g R值,然后将Δl o g R与实测T O C 绘制散点图,分析发现,Δl o g R与实测T O C呈线性相关,故建立二者的线性关系式(图5(a)),长7高阻泥岩T O C预测模型为T O C=26.807ˑΔl o g R-19.476㊂(8)Δl o g R与实测T O C的相关性为85%,吻合度较好㊂图4长7普通泥岩和长7高阻泥岩对比F i g.4C o m p a s i o n o f C h a n g7h i g h r e s i s t a n c em u da n dC h a n g7o r d i n a r y m ud图5长7高阻泥岩㊁长7普通泥岩T O C测井计算模型F i g.5L o g g i n g c a c u l a t e dm o d e l o f t o t a l o r g a n i c c a r b o n f o r C h a n g7h i g h r e s i s t a n c em u d a n dC h a n g7o r d i n a r y m u d2.2.2长7普通泥岩T O C预测模型与长7高阻泥岩相比,长7普通泥岩烃源岩的质量较差㊂长7普通泥岩实测样品共有39个㊂采用相同的方法建立长7普通泥岩T O C预测模型(图5(b)),可表示为T O C=8.4721ˑΔl o g R-5.0164㊂(9)Δl o g R与实测T O C的相关性为89%,吻合度也较好㊂长7高阻泥岩与长7普通泥岩不论是岩心的颜74中国科技论文第16卷色㊁密度,还是测井响应特征都有明显的差异,因此二者的T O C 预测模型也必然存在差异㊂同时,由于沉积时同一盆地不同地区的沉积条件存在一定的差异,导致沉积物必然也会存在一定的不同,因此有机碳预测模型往往具有较大的区域上的局限性㊂本文所建立的长7高阻泥岩T O C 预测模型与长7普通泥岩T O C 预测模型仅适用于彭阳地区延长组长7段,在鄂尔多斯盆地其他区域相同层位的适用性还需进一步验证㊂2.3 效果验证将长7实测样品的T O C 与样品对应深度的测井评价T O C 绘制散点图(图6),开展相关性分析㊂结果显示:长7高阻泥岩的计算T O C 与实测T O C的相关度较高,达到了89%以上(图6(a));长7普通泥岩的计算T O C 与实测T O C 的相关度也达到了91%以上,相关度较高(图6(b))㊂评价结果表明,2个模型的计算T O C 与实测T O C 的相关度均较高,误差不超过11%,计算出的各井目的层段的T O C 值与测试结果相比,也较为符合(图7),反映出所计算的T O C 含量比较准确㊂图6 长7高阻泥岩㊁长7普通泥岩计算T O C 与实测T O C 关系F i g .6 R e l a t i o n s h i p o f a c t u a l l y me a s u r e dT O Ca n d c a c u l a t e dT O Cf o r C h a ng 7hi gh r e s i s t a n c e m u d a n dC h a n g 7o r d i n a r y mud 图7 Z 120㊁Y 65长7烃源岩计算T O C 与实测T O C 关系F i g .7 A c t u a l l y m e a s u r e dT O Ca n d c a c u l a t e dT O Co f C h a n g 7so u r c e r o c k f o r Z 120a n dY 6584第1期王 龙,等:基于改进的Δl o g R 技术的鄂尔多斯盆地彭阳地区长7烃源岩有机碳含量评价由计算的T O C 结果来看,Z 120井长7高阻泥岩T O C 含量很高,最高达到了约30%(图7(a )),Y 65井长7高阻泥岩T O C 含量较高,最高达到了15%(图7(b )),说明彭阳地区长7烃源岩T O C 含量整体较高,烃源岩质量较好㊂此外,前人[21]T O C 研究表明,鄂尔多斯盆地长7优质烃源岩的有机质丰度很高,残余T O C 含量主要分布在6%~14%,最高可达30%~40%,与本文分析结果一致㊂综上所述,改进后的Δl o g R 方法可靠,能够适用于彭阳地区,进一步证实了该方法计算结果的可靠性㊂因此,改进后的Δl o g R 方法可以作为彭阳地区长7烃源岩评价的依据㊂3 长7烃源岩T O C 分布特征对研究区内100余口井长7泥岩进行T O C 的测井识别和评价,结果表明:长7高阻泥岩的T O C 含量大部分在2%~15%,其中优质烃源岩(T O C >2%)占90%(图8(a)),烃源岩丰度整体很高;长7普通泥岩的T O C 含量主要分布在0.5%~2.5%,中等级别的烃源岩占52%(0.5%<T O C <1%),好以上级别的烃源岩(T O C >1%)约占40%(图8(b)㊂图8 长7高阻泥岩㊁长7普通泥岩T O C 频率分布直方图F i g .8 F r e q u e n c y d i s t r i b u t i o n h i s t o g r a mo f T O C f o r C h a n g7h i g h r e s i s t a n c em u d a n dC h a n g 7o r d i n a r y mu d 依据测井评价烃源岩的成果,计算出单井长7普通泥岩和长7高阻泥岩T O C 含量的平均值,绘制T O C 分布图,如图9所示㊂图9 长7高阻泥岩㊁长7普通泥岩T O C 分布F i g .9 T O Cd i s t r i b u t i n g g r a p h o f C h a n g 73h i g h r e s i s t a n c em u d a n dC h a n g 7o r d i n a r y mu d 由图9(a )可见:长7高阻泥岩T O C 含量整体较高,除西北部外,大部分地区都能达到优质烃源岩的级别;东北部是高值区,T O C 含量基本大于5%,东北部的某井T O C 平均含量最高,可达21.5%;西南部T O C 平均含量明显较低,为2%~10%;西部特别是西北部T O C 平均含量最低,小于1%㊂由图9(b )可见:长7普通泥岩T O C 含量整体较低,同样表现出东北部及中部偏西地区的T O C 平均含量相对较高,其他地区较低,大部分地区也能达到中等以上的烃源岩级别;在平面上,中西部某井T O C 平均含量最高,可达3.0%㊂由此可见,长7高阻泥岩烃源岩的质量明显优于长7普通泥岩烃源岩㊂这说明在鄂尔多斯盆地西南部天环凹陷构造带的彭阳地区,其中生界主力烃源岩与盆地内部一致,依然是延长组长7高阻泥岩,同时彭阳地区的长7普通泥岩也具有一定的生烃潜力㊂94中国科技论文第16卷4 结 论1)改进的Δl o g R 技术无需读取A C 与R t 的基线值,同时无需设定L O M 值,不仅减少了人为误差与工作量,还减少引入系统误差,操作便捷㊂2)利用改进的Δl o g R 技术分别建立了长7高阻泥岩和长7普通泥岩T O C 评价模型,对研究区100余口井长7烃源岩T O C 的进行预测,平均误差不超过11%,较好地反映出彭阳地区长7烃源岩T O C 含量㊂3)评价结果显示,鄂尔多斯盆地彭阳地区长7高阻泥岩T O C 含量整体较高,平面上除西北部外,大部分地区都能达到优质烃源岩的级别,尤其以东北部T O C 含量最高;长7普通泥岩T O C 含量整体低于长7高阻泥岩,T O C 含量高值分布在东北部及中部偏西地区,平面上大部分地区也能达到中等以上的烃源岩级别㊂(由于印刷关系,查阅本文电子版请登录:h t t p:ʊw w w .p a p e r .e d u .c n /j o u r n a l /z g k j l w .s h t m l )[参考文献](R e f e r e n c e s)[1] P A S S E YQ R ,C R E A N E YS ,K U L L AJB ,e t a l .Ap r a c t i c a lm o d e l f o ro r g a n i c r i c h n e s s f r o m p o r o s i t y an d r e s i s t i v i t y l o g s [J ].A A P G B u l l e t i n ,1990,74(5):1777-1794.[2] 王贵文,朱振宇,朱广宇.烃源岩测井识别与评价方法研究[J ].石油勘探与开发,2002,29(4):50-53.W A N GG W ,Z H UZY ,Z H UGY .L o g g i n g i d e n t i f i -c a t i o n a n d e v a l u a t i o n o f C a n b r i a n -O r d o v i c i a n s o u r c er o c k s i n s y n e c l i s e o f T a r i mb a s i n [J ].P e t r o l e u mE x pl o -r a t i o n a n dD e v e l o pm e n t ,2002,29(4):50-53.(i nC h i -n e s e)[3] 张志伟,张龙海.测井评价烃源岩的方法及其应用效果[J ].石油勘探与开发,2000,27(3):84-87.Z H A N GZW ,Z H A N GLH.Am e t h o d o f s o u r c e r o c ke v a l u a t i o n b y w e l l -l o g g i i n g a n di t sa p pl i c a t i o nr e s u l t [J ].P e t r o l e u mE x p l o r a t i o n a n dD e v e l o p m e n t ,2000,27(3):84-87.(i nC h i n e s e )[4] 朱光有,金强.烃源岩的非均质性及其研究:以东营凹陷牛38井为例[J ].石油学报,2002,23(5):34-39.Z H UG Y ,J I N Q .S t u d y o ns o u r c e r o c kh e t e r o g e n e i -t y a c a s e o fN i u -38w e l l i nD o n g y i n g d e pr e s s i o n [J ].A c t a P e t r o l e i S i n i c a ,2002,23(5):34-39.(i nC h i n e s e)[5] 杜江民,张小莉,钟高润,等.致密油烃源岩有机碳含量测井评价方法优选及应用 以鄂尔多斯盆地延长组长7段烃源岩为例[J ].地球物理学进展,2016,31(6):2526-2533.D UJM ,Z H A N GXL ,Z H O N GGR ,e t a l .A n a l ys i s o n t h e o p t i m i z a t i o n a n d a p p l i c a t i o n o fw e l l l o g s i d e n t i f i -c a t i o n m e t h o d sf o ro r ga n i cc a rb o nc o n t e n ti ns o u r c e r o c k s o f t h e t i g h t o i l i l l u s t r a t ed b y t he e x a m pl e o f t h e s o u r c e r o c k s o f C h a n g 7m e m b e r o f Y a n c h a n gf o r m a t i o n i nO r d o s b a s i n [J ].P r og r e s s i nG e o ph ys i c s ,2016,31(6):2526-2533.(i nC h i n e s e)[6] 郭泽清,孙平,刘卫红.利用Δl o g R 技术计算柴达木盆地三湖地区第四系有机碳[J ].地球物理学进展,2012,27(2):626-633.G U OZQ ,S U NP ,L I U W H.T h e c a r b o n c a c u l a t i o n b y Δl o g R t e c h n o l o g y inS a n h ua r e ao fQ a i d a m b a s i n [J ].P r o g r e s s i nG e o p h y s i c s ,2012,27(2):626-633.(i nC h i n e s e )[7] 王宗礼,李君,罗强.利用测井资料精细评价冀中坳陷廊固凹陷烃源岩[J ].天然气地球科学,2012,23(3):430-437.W A N GZL ,L I J ,L U O Q .F i n e e v a l u a t i o no f s o u r c er o c kw i t h l o g g i n g d a t a i nL a n g g u s a g o f J i z h o n g d e pr e s -s i o n [J ].N a t u r a lG a sG e o s c i e n c e ,2012,23(3):430-437.(i nC h i n e s e)[8] 金凤鸣,师玉雷,罗强.廊固凹陷烃源岩精细评价研究及应用[J ].中国石油勘探,2012(6):23-29.J I NFM ,S H IYL ,L U OQ .R e s e a r c h a n d a p p l i c a t i o n o fd e t a i l e do nh y d r o c a r b o ns o u r c e r o c ko fL a n g gud e -p r e s s i o n [J ].C h i n aP e t r o l e u m E x pl o r a t i o n ,2012(6):23-29.(i nC h i n e s e)[9] 王振升,滑双军,于学敏,等.歧口凹陷沙河街组烃源岩分级评价及优质烃源岩分布[J ].天然气地球科学,2014,25(12):1896-1902.W A N GZ S ,H U AS J ,Y UXM ,e t a l .G r e d i n g e v a l u -a t i o n a n d h i g h q u a l i t y s o u r c e r o c kd i s t r i b u t i o n i nQ i k o u s a g [J ].N a t u r a l G a sG e o s c i e n c e ,2014,25(12):1896-1902.(i nC h i n e s e)[10]卢双舫,黄文彪,陈方文,等.页岩油气资源分级评价标准探讨[J ].石油勘探与开发,2012,39(2):249-256.L US F ,H U A N G WB ,C H E NFW ,e t a l .G r e d i n ge -v a l u a t i o na n dh i g h q u a l i t y so u r c er o c kd i s t r i b u t i o ni n Q i k o us a g [J ].P e t r o l e u m E x p l o r a t i o na n d D e v e l o p-m e n t ,2012,39(2):249-256.(i nC h i n e s e)[11]曲彦胜,钟宁宁,刘岩,等.烃源岩有机质丰度的测井计算方法及影响因素探讨[J ].岩性油气藏,2011,23(2):80-84.Q U YS ,Z H O N G N N ,L I U Y ,e t a l .U s i n g l o g g i n g m e t h o d s t o c a l c u l a t e o r ga n i cm a t t e r ab u n d a nc e o f s o u r c e r o c k s a nd i t s i n f l ue n c i n gf a c t o r s [J ].L i t h o l o gi cR e s e r -v o i r s ,2011,23(2):80-84.(i nC h i n e s e)[12]刘超,印长海,卢双舫.变系数Δl o g R 烃源岩测井评价技术关键参数厘定方法及应用[J ].天然气地球科学,2015,26(10):1925-1931.L I UC ,Y I NCH ,L USF .P r e d i c t i n g k e yp a r a m e t e r s f o r v a r i a b l e -c o e f f i c i e n t Δl o g R l o g g i n g t e c h n i qu e a n d i t s a p pl i c a t i o n i n s o u r c e r o c k s e v a l u a t i o n [J ].N a t u r a lG a s G e o s c i e n c e ,2015,26(10):1925-1931.(i nC h i n e s e)[13]朱光有,金强,张林晔.用测井信息获取烃源岩的地球化学参数研究[J ].测井技术,2003,27(2):104-109.Z H UGY ,J I N Q ,Z H A N GLY .U s i n g l o g i n f o r m a -t i o n t o a n a l ys e t h e g e o c h e m i c a l c h a r a c t e r i s t i c s o f s o u r c e r o c k s i n J i y a n g d e p r e s s i o n [J ].W e l l L o g g i n g Te c h n o l o -g y,2003,27(2):104-109.(i nC h i n e s e )[14]刘超,卢双舫,黄文彪.Δl o g R 技术改进及其在烃源岩评价中的应用[J ].大庆石油地质与开发,2011,305第1期王 龙,等:基于改进的Δl o g R 技术的鄂尔多斯盆地彭阳地区长7烃源岩有机碳含量评价(3):27-31.L I U C ,L U S F ,H U A N G W B .I m pr o v e m e n to f Δl o g R a n di t sa p p l i c a t i o ni ns o u r c er o c k se v a l u a t i o n [J ].P e t r o l e u m G e o l o g y &O i l f i e l d D e v e l o pm e n ti n D a q i n g,2011,30(3):27-31.(i nC h i n e s e )[15]胡慧婷,卢双舫,刘超,等.测井资料计算源岩有机碳含量模型对比及分析[J ].沉积学报,2011,29(6):1200-1205.H U H T ,L USF ,L I UC ,e t a l .M o d e l s f o r c a l c u l a t -i n g o r g a n i cc a r b o nc o n t e n t f r o ml o g g i n g in f o r m a t i o n :c o m p a r i s o n a n d a n a l y s i s [J ].A c t a S e d i m e n t o l o gi c a S i n -i c a ,2011,29(6):1200-1205.(i nC h i n e s e)[16]喻建,杨亚娟,杜金良.鄂尔多斯盆地晚三叠世延长组湖侵期沉积特征[J ].石油勘探与开发,2010,37(2):181-187.Y UJ ,Y A N GYJ ,D UJL .S e d i m e n t a t i o nd u r i n g t h e t r a n s g r e s s i o n p e r i o d i nL a t eT r i a s s i cY a n c h a n g f o r m a -t i o n ,O r d o s b a s i n [J ].P e t r o l e u m E x pl o r a t i o na n dD e -v e l o pe n t ,2010,37(2):181-187.(i nC h i n e s e )[17]成大伟,袁选俊,周川闽,等.测井岩性识别方法及应用 以鄂尔多斯盆地中西部长7油层组为例[J ].中国石油勘探,2016,21(5):117-126.C H E N GD W ,Y U A N XJ ,Z H O U C M ,e t a l .L o g -g i n g -l i t h o l o g y i d e n t i f i c a t i o nm e t h o d s a n d t h e i r a p pl i c a -t i o n :a c a s e s t u d y o nC h a n g 7Me m b e r i n c e n t r a l -w e s t -e r nO r d o sb a s i n ,N W C h i n a ,O r d o sb a s i n [J ].C h i n aP e t r o l e u m E x pl o r a t i o n ,2016,21(5):117-126.(i n C h i n e s e)[18]黄薇,张小莉,李浩,等.鄂尔多斯盆地中南部延长组7段页岩有机碳含量解释模型[J ].石油学报,2015,36(12):1508-1515.H U A N G W ,Z H A N GXL ,L IH ,e t a l .I n t e r pr e t a t i o n m o d e l o f o r ga n i c c a rb o nc o n t e n t o f s h a l e i nM e m b e r 7o f Y a n c h a n g fo r m a t i o n ,c e n t r a l -s o u t h e r n O r d o s b a s i n [J ].A c t aP e t r o l e i S i n i c a ,2015,36(12):1508-1515.(i nC h i n e s e)[19]杨华,张文正.论鄂尔多斯盆地长7段优质油源岩在低渗透油气成藏富集中的主导作用:地质地球化学特征[J ].地球化学,2005,34(2):147-154.Y A N G H ,Z H A N G WZ .L e a d i n ge f f e c t o f t h e s e v e n t h m e m b e r h i g h -q u a l i t y s o u r c e r o c k o f Y a n c h a n g f o r m a t i o n i nO r d o s b a s i nd u r i n g t h e e n r i c h m e n t o fL o w -pe n e t r a t -i n g l i n -g a s a c c u m u l a t i o n :g e o l o g y a n d g e o c h e m i s t r y[J ].G e o c h i m c a ,2005,34(2):147-154.(i nC h i n e s e )[20]于涛,于江涛,孟宁宁,等.鄂尔多斯盆地油房庄延长组长7油层组沉积环境的粒度特征研究[J ].中国科技论文,2017,24(2):256-262.Y UT ,Y UJT ,M E N G N N ,e t a l .G r a i n s i z e d i s t r i -b u t i o n s o f s e d i m e n t a r y e n v i r o n m e n to fC h a n g 7re s e r -v o i r o fY a n c h a n gf o r m a t i o n ,Y o u f a ng zh u a n g ar e a ,O e -d o s b a s i n [J ].C h i n aS c i e n c e p a pe r ,2017,24(2):256-262.(i nC h i n e s e)[21]王龙.鄂尔多斯盆地彭阳 崇信地区延长组石油勘探潜力研究[D ].青岛:中国石油大学(华东),2012:42.W A N GL .S t u d y o nP e t r o l e u mP o t e n t i a l i t y o f t h eY a n -c h a n g f o r m a t i o ni nP e n g y a n g -C h o n gx i na r e ao fO r d o s b a s i n [D ].Q i n g d a o :C h i n aU n i v e r s i t y ofP e t r o l e u m ,2012:42.(i nC h i n e s e)[22]郑兆惠.辽河滩海西部地区烃源岩测井识别[D ].北京:中国地质大学(北京),2009:23-24.Z H E N GZH.L o g r e c o gn i t i o n f o r s o u r c e r o c k s i nw e s t t i d a l z o n e ,L i a o h e o i l f i e l d [D ].B e i j i n g:C h i n aU n i v e r s i -t y of G e o s c i e n c e s ,2009:23-24.(i nC h i n e s e )[23]刘博.利用测井资料有效识别油页岩技术及其应用 以松辽盆地南部地区为例[D ].长春:吉林大学,2010:38-39.L I U B .T e c h n i q u ea n da p p l i c a t i o no fr e c o g n i z i n g oi l s h a l e e f f e c t i v e l y u s i n g l o g g i n g d a t a t h e e x a m pl e f r o m t h e a r e a i ns o u t h e r nS o n g l i a ob a s i n [D ].C h a n gc h u n :J i l i nU n i v e r s i t y,2010:38-39.(i nC h i n e s e )15。
柴达木盆地第四系疏松砂岩地层生物气含气检测方法探讨田继先;曾旭;王文卓;张绍胜;郭泽清;孔骅【摘要】柴达木盆地三湖地区第四系具有疏松及砂泥岩薄互层等特点,如何利用地震信息开展生物气含气检测一直是该地区研究的重点和难点问题.在梳理研究区含气检测现状的基础上,分析了影响生物气含气检测的两大主要因素,即地层特征和地震反射特征,指出疏松砂岩地层含气之后,地震频率变化最明显,利用频率信息可更好地预测生物气藏,地震资料需要做好保护低频处理.在此基础上提出利用瞬时频率包络加权法预测生物气藏的新方法,并取得了较好的应用效果.【期刊名称】《石油物探》【年(卷),期】2016(055)003【总页数】6页(P408-413)【关键词】第四系;生物气;含气检测;疏松砂岩;低频处理【作者】田继先;曾旭;王文卓;张绍胜;郭泽清;孔骅【作者单位】中国地质科学研究院地质力学研究所,北京100081;中国石油天然气股份有限公司勘探开发研究院廊坊分院,河北廊坊065007;中国石油天然气股份有限公司勘探开发研究院廊坊分院,河北廊坊065007;中国石油天然气股份有限公司青海油田分公司勘探开发研究院,甘肃敦煌736200;中国石油天然气股份有限公司勘探开发研究院廊坊分院,河北廊坊065007;中国石油天然气股份有限公司勘探开发研究院廊坊分院,河北廊坊065007;中国石油天然气股份有限公司勘探开发研究院廊坊分院,河北廊坊065007【正文语种】中文【中图分类】P631柴达木盆地是我国西部中新生代压扭性盆地,第四系沉积物是在新近纪末期新构造运动作用下,柴达木盆地沉积中心由西向东整体迁移的产物,盆地属于典型陆相山间闭塞的坳陷型盆地,地表多为风蚀残丘、沙丘、沙包、盐碱地、硬碱壳地以及盐沼等。
目前已发现台南,涩北一、二号等大型原生生物气田,并发现一批含气构造[1]。
柴达木盆地三湖地区已发现的生物气田在地震剖面上有明显的同相轴下拉异常现象。
柴达木盆地油气资源研究现状及问题刘顺宇;赵荣【摘要】柴达木盆地是我国西部重要的含油气盆地,地质条件非常复杂,发育了石炭系海陆过渡相、侏罗系湖相—湖沼相、古近系咸化湖相和第四系盐湖相4套不同类型的烃源岩层系.盆地的油气资源种类多样,柴西扎哈泉地区新近系上干柴沟组致密油、柴北缘鱼卡地区中下侏罗统页岩气、三湖坳陷第四系生物气等.在前人研究成果的基础上,总结了柴达木盆地油气资源的研究现状及问题,柴西地区一直是勘探的重点,常规的基岩油气藏和非常规的致密油共同开发,基岩油气藏存在埋藏深,地质构造复杂,勘探开发难度大,致密油具有有机质丰度、成熟度偏低,生烃转化率高的特点,因此难以应用常用的致密油评价方法对\"甜点\"进行预测;柴北缘地区页岩气发展迅速,但有机质丰度变化大、有机孔少,以次生孔隙为主,因此\"甜点\"预测难度较大;三湖地区生物气勘探时间长、分布广、资源量大,但鉴于特殊的地质背景,储层预测难度大.【期刊名称】《中州煤炭》【年(卷),期】2019(041)002【总页数】7页(P92-97,101)【关键词】柴达木盆地;致密油;页岩气;生物气【作者】刘顺宇;赵荣【作者单位】东北石油大学,黑龙江大庆 163318;东北石油大学,黑龙江大庆163318【正文语种】中文【中图分类】F426.220 引言柴达木盆地是我国西部重要的含油气盆地,地质条件非常复杂,存在许多勘探的不利因素,如在青藏高原隆升的作用下,受到喜山晚期构造运动的改造,油气保存条件差,柴西地区古今系—新近系所发育的咸化湖相烃源岩,具有有机质丰度、成熟度整体偏低(TOC多小于1.0%,Ro普遍低于1.0%)的特点[1],基于传统的认识,柴西地区生烃潜力小,难以形成大规模的油气田,柴北缘中下侏罗统在燕山晚期的构造作用下,遭受强烈的褶皱、抬升和剥蚀,优质烃源岩分布受限,给油气成藏的研究带来极大的困难,因此引起了许多学者的重视[2]。
第27卷第2期2012年4月(页码:0626-0633)地 球 物 理 学 进 展PROGRESS IN GEOPHYSICSVol.27,No.2Apr.,2012郭泽清,孙 平,刘卫红.利用ΔlogR技术计算柴达木盆地三湖地区第四系有机碳.地球物理学进展,2012,27(2):0626-0633,doi:10.6038/j.issn.1004-2903.2012.02.027.GUO Ze-qing,SUN Ping,LIU Wei-hong.The carbon calculation byΔlogRtechnology in Sanhu area of Qaidam basin.Progressin Geophys.(in Chinese),2012,27(2):0626-0633,doi:10.6038/j.issn.1004-2903.2012.02.027.利用ΔlogR技术计算柴达木盆地三湖地区第四系有机碳郭泽清1,2, 孙 平2, 刘卫红2(1.中国石油勘探开发研究院,北京100083; 2.中国石油勘探开发研究院廊坊分院,廊坊065007)摘 要 利用ΔlogR技术计算有机碳是一种不同于传统地球化学测试得到的有机碳方法.该方法利用烃源岩有机碳含量在测井曲线(声波时差、电阻率)上的响应特征,建立测井解释模型,计算TOC的连续分布值,从而得到准确的有机碳含量,补充取心资料的不足.本文应用ΔlogR技术,结合实测数据,建立起三湖地区测井资料与有机碳含量之间的定量关系模型,然后利用该模型对三湖地区25口重点探井进行有机碳计算.结果表明:有机碳总体上数值含量偏低,平均值为0.33%;在纵向上,具有2个峰值,其中k5-k6峰值段是第四系最大湖泛期形成的沉积;在平面上,有机碳高值区位于台南-涩北构造以北和台吉乃尔与南陵丘以南区域,而南斜坡和涩南1井以东区域有机碳含量较低.本次研究获取了较为准确的烃源岩评价参数,最终得到该区第四系天然气资源量得到油田认可.关键词 柴达木,三湖地区,第四系,ΔlogR方法,有机碳含量,资源量doi:10.6038/j.issn.1004-2903.2012.02.027 中图分类号 P631 文献标识码 A 收稿日期 2011-06-12; 修回日期 2011-11-12. 投稿网址 http//www.progeophys.cn基金项目 国家大型油气田及煤层气开发科技重大专项(2008ZX05)资助.作者简介 郭泽清,男,1976年生,山东沾化人,现为中国石油勘探开发研究院博士研究生,主要从事天然气地质研究.(E-mail:gzq1463@sina.com)The carbon calculation byΔlogRtechnologyin Sanhu area of Qaidam basinGUO Ze-qing1,2, SUN Ping2, LIU Wei-hong2(1.Research Institute of Petroleum Exploration and Development,CNPC,Beijing100083,China;2.Langfang Branch,Research Institute of Petroleum Exploration and Development,CNPC,Langfang065007,China)Abstract The carbon calculation byΔlogRtechnology is an organic carbon gained method different from theconventional geochemical testing.According to logging(interval transit time and resistivity)response characteristicsof the organic matter content in source rocks,the relation formula of organic carbon content is established by welllog,which is used to calculate the continuous TOC,thereby the accurate contents of organic carbon are obtained,and these data can make up the deficiency coring data.The quantitative interpretation model between well logs andcarbon content has been established by this technology and measured data in Sanhu area.Then,the organic carboncontents of 25wells were calculated and the characteristics of vertical and plane distribution were pointed out.Theaverage content of organic carbon is about 0.33%,and the carbon content in the vertical direction has 2peaks.Oneof the peak is the deposition formed during the maximum flooding period of Quaternary.The high value is in thearea which is located in the north of Tainan-Sebei and the south of the Taijinaier-Nanlingqiu,and the low value areais in the south slope and east of Senan1well.Therefore,more accurate evaluation parameters of source rock wereobtained,and according these data,the final result for the Quaternary natural gas resources extent is obtained.Keywords Qaidam,Sanhu area,Quaternary,ΔlogRtechnology,organic carbon content,resource extent 2期郭泽清,等:利用ΔlogR技术计算柴达木盆地三湖地区第四系有机碳0 引 言ΔlogR技术是埃克森(Exxon)和埃索(Esso)公司于1979年开发的一种烃源岩有机碳含量测井预测技术,经广泛试用证实,该技术既适用于碳酸盐岩烃源岩,也适用于碎屑岩烃源岩,而且能精确根据烃源岩成熟度计算有机碳含量[1].利用测井信息纵向分辨率高的特点以及有机质在测井曲线上的特殊响应特征,可以建立起测井资料与烃源岩有机质含量间的定量关系模型,进而计算出烃源岩段有机碳含量的连续分布值.这无疑可以有效地弥补实验室测样的不足,为烃源岩的评价研究提供更加准确、合理的参数值[2].ΔlogR技术已经应用于烃源岩评价[2-6]和层序地层划分[7-10]中.柴达木盆地三湖地区西起船形丘构造、东止南北霍布逊湖、北邻陵间断裂、南到昆仑山前,面积约30000km2,因其境内分布着东西台吉乃尔湖、涩聂湖和达布逊湖三个大型盐湖而得名[11].柴达木盆地第四系的形成是在第三纪末期新构造运动作用下,柴达木盆地沉积中心由西向东整体迁移的产物,最大沉积厚度超过2700m,是一套连续沉积的湖相沉积[12].在第四系中发现了世界上最大的生物气田[13,14].该区有机碳含量低,3000多块气源岩有机碳分析结果表明,有机碳含量普遍在0.1%~0.4%之间,平均仅为0.30%[15,16].而且有机质未成熟,由于本区第四系仍处于早期成岩阶段,镜质体反射率在0.2%~0.47%之间.由于有机碳含量是天然气资源评价的关键参数.因此建立起测井资料与机质含量间的定量关系模型,进而计算出烃源岩段有机碳含量在纵、横向上的分布规律,对于正确认识三湖地区第四系生物气源岩的有机质丰度,对于了解该区生物气生成条件、正确评价生物气资源及寻找有利的勘探领域具有重要意义.此外,该方法也首次应用于地层未成岩、源岩未成熟的地区.1 基本原理将刻度合适的孔隙度曲线(声波测井、补偿中子、密度)叠合在电阻率曲线上,对于贫含有机质层段,2条曲线相互重合或平行;而在富含有机质层段中两者的分离,这主要是由于低密度、低速率(高声波时差)的干酪根响应和地层流体在电阻率曲线上的反映[2].据声波时差-电阻率叠加计算ΔlogR的方程是ΔlogR=log(R/R基线)+(Δt-Δt基线)/164,(1)式(1)中,ΔlogR为实测曲线间距在对数电阻率坐标上的读数,单位为Ω·m;R为测井仪实测的电阻率,单位为Ω·m;R基线为基线对应的电阻率,单位为Ω·m;Δt为实测传播时间,单位为μs/m;Δt基线为基线对应的传播时间,单位为μs/m;系数1/164基于每164μs/m的声波时差Δt,相当于电阻率Rt的一个对数坐标单位.在应用时,电阻率曲线刻度为1个数量级的对数电阻率刻度,对应的声波时差为164μs/m(或100μs/ft)的间隔.两条曲线在一定深度范围内“一致”或完全重叠时即为基线.确定基线之后,用两条曲线之间的间距来识别富含有机质的层段.两条曲线的间距记为ΔlogR.在移动曲线刻度时,一般保持声波曲线不动,移动电阻率曲线,但应保持电阻率曲线刻度为每一个数量级对数电阻率刻度对应声波时差为164μs/m(或100μs/ft)间隔的原则.Passey等经过大量统计分析,提出了能够计算不同成熟度条件下有机碳含量的测井评价方法,并制订了有机碳计算图版[17].经后人改进,认为ΔlogR为成熟度Ro的函数,且与TOC(有机碳总量)呈线性相关[2]:TOC=ΔlogR*10-0.944Ro+1.5374,(2)Ro为镜质体反射率,单位为%.由于三湖地区第四系地层无论是泥岩还是砂岩均含一定量的有机碳,在电阻率测井曲线和自然伽玛曲线重叠段当作非烃源岩段,实际上含一定的有机碳背景值ΔTOC.因此,根据ΔlogR计算的TOC还需加一定的ΔTOC,其定量评价模型为TOC=ΔlogR*10-0.944Ro+1.5374+ΔTOC,(3)ΔTOC为有机碳背景值(%).实际上,对于样品较多的层位,有机碳(TOC)与分离量(ΔlogR)呈线性关系(一元一次方程):TOC=a×ΔlogR+ΔTOC.(4)2 三湖地区第四纪有机碳的计算涩3-2-4井是位于涩北二号气田的开发井,该井在520~580m井段测试有机碳数据247个,测试样品比较多,而且连续测量,可以用最小二乘法拟合出实测有机碳TOC和计算得到的ΔlogR之间的关系式:TOC=0.17×ΔlogR+0.23(a=0.17,ΔTOC=0.23),用该关系式得到的计算有机碳含量曲线与实测726地球物理学进展 http//www.progeophys.cn 27卷图1 涩3-2-4井520~580m井段有机碳计算值与实测值对比图Fig.1 The comparison between the calculated and measured value of organic carbon from 520mto 580min well Se3-2-4图2 涩3-2-4井1070m~1100m井段有机碳计算值与实测值对比图Fig.2 The comparison between the calculated and measured value of organiccarbon from 1070mto 1100min well Se3-2-4有机碳数值吻合很(图1),相似系数为R=0.91.涩3-2-4井在1070~1100m井段测试有机碳数据88个,用最小二乘法拟合出有机碳TOC和ΔlogR之间的关系式为TOC=0.33×ΔlogR+0.27(a=0.33,ΔTOC=0.27),用该关系式得到的计算有机碳含量曲线与实测有机碳数值吻合很好(图2),相似系数为R=0.90.826 2期郭泽清,等:利用ΔlogR技术计算柴达木盆地三湖地区第四系有机碳图3 台南5井800m~1500m井段有机碳计算值与实测值对比图Fig.3 The comparison between the calculated and measured value of organiccarbon from 800mto 1500min well Tainan5将上述公式应用于台南5(图3)、金达3、聂深1、涩东11、盐14、伊深1、和伊克2等有实测有机碳数值的井进行验证,表明计算值和实测值吻合.说明这两个公式普遍适用,但分别适用于不同的深度范围,分别是,当深度小于900m时,有TOC=0.17×ΔlogR+0.23,当深度大于900m时,有TOC=0.33×ΔlogR+0.27.因此该区有机碳TOC是ΔlogR的分段函数.原因是该地区以埋藏深度900m为界,在900m以上泥质岩压实程度小于砂质,而900m以下泥质岩压实程度超过砂质[18].因此,由于砂泥岩压实程度的差异,而导致计算公式的不同.2.1 纵向分布特征将上述公式应用于该区25口区域重点探井计算各层段有机碳含量,三湖地区七个泉组从上至下识别出k1、k2、k3、k4、k5、k6、k7、k9、k10、k11、k13(第四系底)共11层岩电标志层,将地层划分为10个层段[19].以位于涩北一号气田的涩23井(图4)和涩北二号气田的涩21井(图5)为例,可以看出有机碳分布具有2个峰值,分别是k5~k6层段和k1~k2层段.其中k5~k6层段是是第四系最大湖泛期形成的沉积,有效的证明了最大湖泛面的位置.k1~k2层段的峰值可能是由于地层欠压实导致声波时差增大,导致ΔlogR增大,使得计算结果偏大造成的;也可能是由于浅层有机质未达到产甲烷高峰,剩余有机碳比其它层段多而导致的.其它井的有机碳分布特征与此相似,这里不再赘述.从25口井的平均有机碳分布图(图6)来看,除了k1~k2层段,k5~k6层段仍然是有机碳最大值分布层段.总体上有机碳数值分布比较平均,平均值为0.33%,与实验测得的数值吻合.2.2 平面分布特征首先计算区域上25口重点探井纵向上有机碳数值,分层段统计有机碳平均含量,然后分层段平面926地球物理学进展 http//www.progeophys.cn 27卷图4 涩23井各层段有机碳含量Fig.4 Carbon content for each layer of well Se23图5 涩21井各层段有机碳含量Fig.5 Carbon content for each layer of well Se21图6 全区25口井各层段有机碳分布图Fig.6 Carbon content for each layer of 25wells036 2期郭泽清,等:利用ΔlogR技术计算柴达木盆地三湖地区第四系有机碳图7 三湖地区k5~k6层段有机碳平面分布特征Fig.7 Organic carbon plane distribution characteristics of k5~k6layer in Sanhu area图8 三湖地区第四系总生气强度图(单位:108/m3)Fig.8 Gas generation intensity of the Quaternary in Sanhu area136地球物理学进展 http//www.progeophys.cn 27卷 成图,完成10个层段有机碳平面分布图.这里以有机碳最大值分布层段k5~k6层段为例说明.有机碳高值区位于台南—涩北构造以北和台吉乃尔与南陵丘以南区域,南斜坡和涩南1井以东区域有机碳含量较低.其它层段有机碳分布规律与此相似,这里不再赘述.有机碳的这种分布规律与泥质含量分布规律一致[20].其原因是柴达木盆地三湖地区第四纪物源呈明显的不对称性,南部昆仑山的物源是主导物源,锡铁山和埃姆尼克山的物源是次要物源.尤其是在乌图美仁、格尔木分别发育了两个规模较大的水系进入盆地[20].受淡水供给的影响,靠近昆仑山的盆地南斜坡盐度偏低,有机质含量偏低,导致有机碳偏低,与此相反,远离物源的台南—涩北构造以北由于盐度偏高,有利于有机质保存,导致有机碳偏高.2.3 计算资源量本次研究采用盆地模拟方法计算资源量,该方法的关键参数为产气率图版、有机碳含量和运聚系数.产气率图版根据实验数据获得,有机碳含量由本次应用ΔlogR技术计算获得,运聚系数为生气量向资源量转化的关键参数,直接采用第三次资源评价时采用的参数,为0.92%[21].第四系按照10个层段,每一套地层形成一幅生烃强度图,然后将生烃强度累加,形成第四系累计生烃强度图(图8).模拟结果表明:生气强度在20~80×108 m3/km2,生气强度大于20×108 m3/km2的面积为27650km2.生气强度与面积的乘积就得到第四系总的生气量,再乘以运聚系数,就可以得到资源量.由于计算参数更加可靠,计算结果得到油田认可.3 结 论3.1 三湖地区可以用ΔlogR方法计算有机碳含量,吻合程度较高.有机碳TOC是ΔlogR的分段函数,当深度小于900m时,TOC=0.17×ΔlogR+0.23,当深度大于900m时,TOC=0.33×ΔlogR+0.27.3.2 总体上有机碳数值含量偏低,平均值为0.33%,在纵向上具有2个峰值,分别是k5~k6层段和k1~k2层段,其中k5~k6层段是是第四系最大湖泛期形成的沉积,有效的证明了最大湖泛面的位置.3.3 有机碳高值区位于台南—涩北构造以北和台吉乃尔与南陵丘以南区域,南斜坡和涩南1井以东区域有机碳含量较低,这种分布特征与湖水的盐度分布规律一致,有机碳含量受湖水盐度的控制.3.4 盆地模拟结果得到资源量.由于有机碳参数的获得更加可靠,计算一结果得到油田认可.参 考 文 献(References):[1] Orr F M Jr,Yu A D,Lien C L.Phase behavior of CO2andcrude oil in low-temperature reservoirs.SPE Journal,1981,21(4):480-492.[2] 于翔涛.测井技术在烃源岩评价中的应用.长江大学学报(自然科学版),2009,6(2):198-200.Yu X T.Application of well logging technology to sourcerocks evaluation.Journal of Yangtze University(NaturalScience Edition)(in Chinese),2009,6(2):198-200.[3] 徐思煌,朱义清.烃源岩有机碳含量的测井响应特征与定量预测模型—以珠江口盆地文昌组烃源岩为例.石油实验地质,2010,32(3):290-295.Xu X H,Zhu Y Q.Well logs response and prediction model oforganic carbon content in source rocks-a case study from thesource rock of Wenchang formation in the pearl mouth basin.Petroleum Geology &Experiment(in Chinese),2010,32(3):290-295.[4] 朱振宇,刘洪,李幼铭.ΔlogR技术在烃源岩识别中的应用与分析.地球物理学进展,2003,18(4):647-649.Zhu Z Y,Liu H,Li Y M.The analysis and application ofΔlogRmethod in the source rock’s identification.Progress inGeophysics(in Chinese),2003,18(4):647-649.[5] 张志伟,张龙海.测井评价烃源岩的方法及其应用效果.石油勘探与开发,2000,27(3):84-87.Zhang Z W,Zhang L H.A method of source rock evaluationby well-logging and its application result.PetroleumExploration and Development(in Chinese),2000,27(3):84-87.[6] 贺君玲,邓守伟,陈文龙等.利用测井技术评价松辽盆地南部油页岩.吉林大学学报(地球科学版),2006,36(6):909-914.He J L,Deng S W,Chen W L,et al.Evaluation of oil shale inthe southern Songliao basin using logging techniques.Journalof Jilin University(Earth Science Edition)(in Chinese),2006,36(6):909-914.[7] 闫建平,蔡进功,赵铭海.测井信息用于层序地层单元划分及对比研究综述.地层学杂志,2009,33(4):441-450.Yan J P,Cai J G,Zhao M H.Advances in the study ofsequence stratigraphic division and correlation using well loginformation.Journal of Stratigraphy(in Chinese),2009,33(4):441-450.[8] 操应长,姜在兴,夏斌等.利用测井资料识别层序地层界面的几种方法.石油大学学报(自然科学版),2003,27(2):23-26.Cao Y C,Jiang Z X,Xia B,et al.Some methods foridentifying sequence boundaries and condensation sectionsusing well logging.Journal of the University of Petroleum,China(Edition of Natural Science)(in Chinese),2003,27236 2期郭泽清,等:利用ΔlogR技术计算柴达木盆地三湖地区第四系有机碳(2):23-26.[9] 王卫红,姜在兴,操应长等.测井曲线识别层序边界的方法探讨.西南石油学院学报,2003,25(3):1-4.Wang W H,Jiang Z X,Cao Y C,et al.Discussion of themethods of using logging curves to recognize sequenceboundaries.Journal of Southwest Petroleum Institute(inChinese),2003,25(3):1-4.[10] 于均民,李红哲,刘震华等.应用测井资料识别层序地层界面的方法.天然气地球科学,2006,17(5):736-742.Yu J M,Li H Z,Liu Z H,et al.Identifying sequencestratigraphy boundaries using well logging.Natural GasGeoscience(in Chinese),2006,17(5):736-742.[11] 党玉琪,张道伟,徐子远等.柴达木盆地三湖地区第四系沉积相与生物气成藏.古地理学报,2004,6(1):110-118.Dang Y Q,Zhang D W,Xu Z Y,et al.Sedimentary faciesand biogenic gas pool of the Quaternary of Sanhu area inQaidam basin.Journal of Palaeogeography(in Chinese),2004,6(1):110-118.[12] 王永军,郭泽清,刘卫红等.柴达木盆地东部三湖地区四系米兰柯维奇旋回分析.地球物理学进展,2007,22(2):544-551.Wang Y J,Guo Z Q,Liu W H,et al.Analysis ofMilankovitch cycles of quaternary in Sanhu area,easternQaidam basin.Progress in Geophysics(in Chinese),2007,22(2):544-551.[13] 尚永生,杨长春,王真理等.柴达木盆地第四系气藏烃类检测研究.地球物理学进展,2007,22(3):831-835.Shang Y S,Yang C C,Wang Z L,et al.The research ofhydrocarbon indication on gas reservior of Qaidam basin.Progress in Geophysics(in Chinese),2007,22(3):831-835.[14] 苏明军,王西文,何亨华等.柴达木盆地三湖地区北斜坡岩性地层气藏形成的地质条件与预测技术研究.地球物理学进展,2008,23(4):1190-1198.Su M J,Wang X W,He H H,et al.The lithologic andstratigraphic gas reservoirs in the north slope of the Sanhuarea,Qaidam basin:a study on their genetic conditions andcharacterization techniques.Progress in Geophysics(inChinese),2008,23(4):1190-1198.[15] 戚厚发,关德师,钱贻伯等.中国生物气成藏条件.北京:石油工业出版社,1997.Qi H F,Guan D S,Qian Y B,et al.AccumulationConditions of Biological Gas in China(in Chinese).Beijing:Petroleum Industry Press,1997.[16] 顾树松.柴达木盆地东部第四系气田形成条件及勘探实践.北京:石油工业出版社,1993.Gu S S.The Forming Conditions and Exploration Practice ofQuaternary Gas Field in the Eastern Qaidam Basin(inChinese).Beijing:Petroleum Industry Press,1993.[17] Passey Q R,Creaney S,Kulla J B,et al.A practical modelfor organic richness from porosity and resistivity logs.AAPGBulletin,1990,74(12):1777-1794.[18] 张林,孙平,张绍胜等.柴达木盆地第四系岩性气藏成藏条件.天然气工业,2009,29(9):32-34.Zhang L,Sun P,Zhang S S,et al.An analysis ofhydrocarbon pooling conditions of the Quaternary lithologicgas reservoirs in Qaidam Basin.Natural Gas Industry(inChinese),2009,29(9):32-34.[19] 孙平,郭泽清,张林等.柴达木盆地三湖北斜坡岩性气藏勘探与发现.中国石油勘探,2011,16(1):25-31.Sun P,Guo Z Q,Zhang L,et al.Lithologic gas reservoirexploration and discovery in northern Sanhu slope,Qaidambasin.China Petroleum Exploration(in Chinese),2011,16(1):25-31.[20] 郭泽清,孙平,张林等.柴达木盆地东部三湖地区第四纪沉积物源分析.第四纪研究,2010,35(5):1044-1054.Guo Z Q,Sun P,Zhang L,et al.Source analysis ofquaternary sediments in the Sanhu area,Eastern QaidamBasin.Quaternary Sciences(in Chinese),2010,35(5):1044-1054.[21] 姜桂凤,孔红喜,侯泽生等.柴达木盆地生物气资源潜力评价.新疆石油地质,2005,26(4):363-366.Jiang G F,Kong H X,Hou Z S,et al.Assessment ofbiogenetic gas reservoir potential in Qaidam basin.XinjiangPetroleum Geology(in Chinese),2005,26(4):363-366.336。