下载文档原格式
122为了满足复兴地区页岩气精细目标勘探的需求，充分发挥“两宽一高”地震资料处理的优势，进一步加强复兴侏罗系陆相页岩油气目标区复杂岩性目标的落实，针对复兴地区地震资料特点，作者对此展开宽方位矢量偏移距（OVT）处理技术[1]研究，首先介绍地震资料OVT处理技术的处理流程，对各项OVT处理技术及效果进行阐述，通过OVT处理技术的应用，在地震资料保真保幅的基础上有效提高了分辨率，提高了地震资料对复杂勘探目标的储层预测和地质体识别能力，满足了复兴地区岩性勘探、各向异性裂缝预测的需求。
1 OVT 处理技术根据复兴地区大连片地震资料特点，在常规地震资料处理流程的基础上，处理时强化各向异性分析，做好OVT域道集的方位各项异性处理，可选择进行OVT域道集识别裂缝方向再细分方位处理，抽取OVT域道集。
另外，OVT域规则化处理补齐OVT域内反射面元的空缺道，消除方位各向异性差异。
完成OVT域全方位偏移和分方位偏移，确保偏移成像准确，偏移后OVT域道集的AVO、AVAZ反射特征与实钻井相一致。
1.1 十字排列抽取及OVT 片划分合理选择划分方案对后续开展OVT域处理至关重要。
理论上划分后的OVT片数应全等于区域覆盖次数，由于本文研究的资料横纵比低，十字排列又具有3个一滚动的周期，对称法划分后产生12×10个面元，其中比较完整的数据量比较大的面元为7×10，横向两边上分别有个半片，在一个周期内滚动时分别生成6×10、7×10、6×10个OVT片；中心法划分后产生9×12个面元，其中比较完整的数据量比较大的面元为6×11，纵向两边上分别有个半片，在一个周期内滚动时分别生成6×11、7×11、6×11个OVT片。
按上述方法抽取的OVT道集的炮检距范围较大，是炮线距或检波线距的2倍。
如果采集时炮线距和检波线距在200m~300m，则OVT道集的炮检距在一个比较合理的范围内；如果炮线距和检波线距过大，OVT内的炮检距和方位角范围增大，失去OVT域处理的OVT处理技术在复兴地区大连片地震资料中的应用张颢 刘倩 乔小玲中国石油化工江汉油田分公司勘探开发研究院 湖北 武汉 430400摘要：针对复兴地区精细油气勘探的需求，基于“两宽一高”地震采集资料的特点，开展了宽方位矢量偏移距（OVT）域处理。
地震波形解释技术讲解地震波形解释是地球物理学中的一项重要工作，它通过分析地震记录中的波形信息，了解地下岩石结构、地震发生机制以及地震破裂过程等相关信息。
地震波形解释技术在地质勘探、矿产资源探测、地震监测等领域都有广泛应用。
本文将介绍地震波形解释技术的基本概念、方法和应用。
一、地震波形解释技术概述地震波形解释技术是根据地震波在地下介质中的传播和反射、折射等现象，通过分析波形记录来确定地下岩石的物理性质和结构。
地震波形记录中包含了地震波在地震源和地表接收点之间传播的信息，将这些信息进行处理和解释，就可以获取地下结构的相关信息。
二、地震波形解释技术的方法1. 震相分析法震相分析法是一种常用的地震波形解释方法。
它通过分析地震记录中的不同震相的到达时间和振幅，来推断不同岩石层界面的位置和性质。
震相分析法包括初动到时提取、振幅分析和速度分析等步骤。
通过田间实测和实验室分析，可以建立震相的速度表，利用速度表来解释地震记录中的波形信息。
2. 反射波形解释法反射波形解释法是根据地震记录中的反射波形特征，来推断地下界面的形态和属性。
在地震记录中，反射波是震源发射的地震波在地下岩石界面上发生反射后返回地表接收到的波形。
通过分析反射波的振幅、频率、相位等特征，可以判断反射面的位置、走向、倾角和反射系数等参数，从而得到地下构造的信息。
3. 折射波形解释法折射波形解释法主要应用于地下介质存在不均匀性的情况。
当地震波从一个介质传播到另一个介质时，波的传播方向会发生改变，这就是折射现象。
通过分析折射波的特征，可以计算出介质的折射系数、折射角度等参数，进而推断地下介质的物理特性。
三、地震波形解释技术的应用1. 地质勘探地震波形解释技术在地质勘探中有着广泛的应用。
通过分析地震记录中的波形信息，可以了解地下岩石的层序、岩性、构造等特征，为勘探活动提供重要信息。
地震波形解释技术广泛应用于石油、天然气、矿产资源等勘探项目中。
2. 地震监测地震波形解释技术在地震监测中也发挥着重要作用。
基于地震资料的薄互层储层精细地质建模石莉莉【摘要】某地区 C 油田钻井少、井距大，且具有储层砂体薄、横向连通差、砂体相变快的特点，仅依靠测井数据插值难以准确刻画砂体边界，建立的模型精度较低，直接影响到剩余油分布预测和油藏开发调整。
基于地震资料的精细地质建模，应用地震解释的断层和层位数据，建立储层精细构造格架；通过地震反演建立测井解释岩相与地震数据之间的联系，建立相模型时，利用协克里金方法，将地震反演的波阻抗作为第2变量约束井间砂体预测。
结果表明，井震结合研究断层更精确，地震解释面的约束使井间构造幅度保留了地震解释的相对趋势，更好地描述出井间微小构造；地震反演波阻抗的协同应用使模型在平面上融合了地震数据中观测到的大尺度结构和测井数据的地质特征，垂向上则在反演结果的基础上进一步趋于测井尺度，有效地降低了井间预测的不确定性，提高了岩相建模的精度。
【期刊名称】《长江大学学报（自然版）理工卷》【年(卷),期】2016(013)001【总页数】4页(P12-15)【关键词】井震结合;薄互层储层;精细地质建模;地震反演【作者】石莉莉【作者单位】中石化河南油田分公司采油一厂，河南桐柏 474780【正文语种】中文【中图分类】P631.44某地区C油田钻井少、井距大，整体上勘探程度较低，对研究区砂体空间展布规律认识不清；其储层为三角洲前缘沉积，单砂体厚度薄，砂泥叠置交互严重，为典型的薄互层储层，同时受储层非均质性强的影响，储层横向变化快，油气分布规律十分复杂[1,2]；亟需建立精细的储层地质模型来揭示砂体空间展布特征。
为此。
笔者运用基于地震资料的薄互层储层精细地质建模方法[3,4]对C油田进行研究。
基于地震资料的建模方法主要包括2个方面：井震结合构造框架模型建立和井震结合岩相模型、属性模型建立。
精确的构造框架模型为后期插值运算提供了边界限制，构造模型不准确会造成地震数据重采样到模型中时出现较大误差，从而使地震约束失去意义。
[收稿日期]20220107[基金项目]国家科技重大专项 川西源内成藏体系成藏机理研究 (2016Z X 05002-004-009)㊂ [第一作者]谢锐杰(1965),男,博士,教授,博士生导师,现主要从事层序地层学和储层地质方向的教学与科研工作,249539680@q q.c o m ㊂谢锐杰,曹永强,李瑶,等.丘东气田西山窑组波形指示反演砂体识别[J ].长江大学学报(自然科学版),2023,20(2):29-36.X I ERJ ,C A O Y Q ,L IY ,e t a l .S a n db o d y i d e n t i f i c a t i o nb y w a v e f o r mi n d i c a t i o n i n v e r s i o n i nX i s h a n y a oF o r m a t i o no fQ i u d o n g Ga s F i e l d [J ].J o u r n a l o fY a n g t z eU n i v e r s i t y (N a t u r a l S c i e n c eE d i t i o n ),2023,20(2):29-36.丘东气田西山窑组波形指示反演砂体识别谢锐杰1,2,曹永强1,2,李瑶1,2,唐欢欢1,21.长江大学地球物理与石油资源学院,湖北武汉4301002.油气资源与勘探技术教育部重点实验室(长江大学),湖北武汉430100[摘要]丘东气田西山窑组发育一套以扇三角洲前缘为主的灰色砂岩,其是主要的含气层段㊂为了查明砂体的分布规律,通过测井曲线敏感性分析发现:有效砂岩主要集中在高波阻抗区域,非储层主要集中在低波阻抗区域,可以明显区分砂泥岩,分界线在8000(g㊃m )/(s ㊃c m 3);深侧向电阻率R L L D 曲线相较于自然伽马㊁自然电位等测井曲线,对岩性的识别敏感度较大,因此选用R L L D 曲线重构进行波形指示反演㊂在平滑半径㊁有效样本数㊁频率参数等反演参数对比分析的基础上,对重构后的波阻抗曲线开展反演工作㊂反演结果表明,垂向上对较薄砂体的识别能力大大增强,反演精度较高,可以刻画出厚度5m 左右的单砂体;砂体分布从南往北逐渐减薄,东北部砂岩不发育,物源主要来自南部㊂经后期钻探验证,该方法有效增强了对砂岩的识别能力,缩小了岩性气藏的勘探范围,提高了钻探成功率㊂[关键词]波形指示反演;砂体识别;丘东气田;西山窑组[中图分类号]P 631.44[文献标志码]A [文章编号]16731409(2023)02002908S a n db o d y i d e n t i f i c a t i o nb y w a v e f o r mi n d i c a t i o n i n v e r s i o n i nX i s h a n ya o F o r m a t i o no f Q i u d o n g Ga s F i e l d X I ER u i j i e 1,2,C A O Y o n g q i a n g 1,2,L IY a o 1,2,T A N G H u a n h u a n 1,21.S c h o o l o fG e o p h y s i c s a n dP e t r o l e u m R e s o u r c e s ,Y a n g t z eU n i v e r s i t y,W u h a n430100,H u b e i 2.K e y L a b o r a t o r y o fE x p l o r a t i o n T e c h n o l o g i e s f o rO i l a n dG a sR e s o u r c e s ,M i n i s t r y o fE d u c a t i o n (Y a n g t z e U n i v e r s i t y ),W u h a n 430100,H u b e iA b s t r a c t :T h eX i s h a n y a oF o r m a t i o n o fQ i u d o n g G a sF i e l d d e v e l o p s a s e t o f g r e y s a n d s t o n e s d o m i n a t e db yt h e f r o n t o f f a nd e l t a ,w h i c h i s t h em a i n g a s -b e a r i n g i n t e r v a l .I n o r d e r t o f i n d o u t t h e d i s t r i b u t i o n l a wo f s a n db o d i e s ,t h r o u g h t h e s e n s i t i v i t y a n a l y s i s o f l o g g i n g c u r v e s ,i t i s f o u n d t h a t t h e e f f e c t i v e s a n d s t o n e i sm a i n l y c o n c e n t r a t e d i n t h eh i ghw a v e i m p e d a n c e a r e a ,a n dt h en o n -r e s e r v o i r i s m a i n l y c o n c e n t r a t e di nt h el o w w a v e i m p e d a n c ea r e a ,w h i c hc a nc l e a r l y d i s t i n g u i s hs a n da n dm u d s t o n ew i t h t h e b o u n d a r y a t 8000(g ㊃m )/(s ㊃c m 3).C o m p a r e dw i t h l o g g i n g cu r v e s s u c h a s n a t u r a l g a mm aa n ds p o n t a n e o u s p o t e n t i a l ,t h ed e e p l a t e r a l r e s i s t i v i t y c u r v eR L L D h a sa g r e a t e rs e n s i t i v i t y t ol i t h o l o g yi d e n t i f i c a t i o n .T h e r e f o r e ,R L L Dc u r v e r e c o n s t r u c t i o nw a ss e l e c t e df o rw a v e f o r mi n d i c a t i o n i n v e r s i o n .O nt h eb a s i so f t h e c o m p a r a t i v e a n a l y s i s o f i n v e r s i o n p a r a m e t e r s s u c h a s s m o o t h r a d i u s ,v a l i d s a m p l e n u m b e r a n d f r e q u e n c y pa r a m e t e r s ,t h e i n v e r s i o nw o r kw a s c a r r i e do u t o n t h e r e c o n s t r u c t e dw a v e i m p e d a n c e c u r v e s .T h e i n v e r s i o n r e s u l t s s h o wt h a t t h e ab i l i t y t o i d e n t i f y t h i n n e r s a n db o d i e s i n t h e v e r t ic a ld i re c t i o n i s g r e a t l y e n h a n c e d ,t h e i n v e r s i o n a c c u r a c y i s h i g h ,a n d t h e s i n g l e s a n db o d y w i t h a t h i c k n e s s of a b o u t 5mc a nb e d e s c r i b e d ;T h e s a n db o d y d i s t r i b u t i o ng r a d u a l l y d e c r e a s e s f r o ms o u th t on o r t h ,t h e s a n d s t o n ei nt h en o r t h e a s t i sn o td e v e l o p e d ,a n dt h e p r o v e n a n c e m a i n l y c o m e s f r o mt h es o u t h .A f t e r l a t e r d r i l l i n g v e r i f i c a t i o n ,t h i sm e t h o de f f e c t i v e l y e n h a n c e s t h e a b i l i t y t o i d e n t i f y s a n d s t o n e ,n a r r o w s t h ee x p l o r a t i o n s c o p e o f l i t h o l o g i c g a s r e s e r v o i r s ,a n d i m p r o v e s t h e d r i l l i n g su c c e s s r a t e .K e yw o r d s :w a v e f o r mi n d i c a t i o n i n v e r s i o n ;s a n db o d y i d e n t i f i c a t i o n ;Q i u d o n g G a sF i e l d ;X i s h a n y a oF o r m a t i o n ㊃92㊃长江大学学报(自然科学版) 2023年第20卷第2期J o u r n a l o fY a n g t z eU n i v e r s i t y (N a t u r a l S c i e n c eE d i t i o n ) 2023,V o l .20N o .2Copyright ©博看网. All Rights Reserved.㊃03㊃长江大学学报(自然科学版)2023年3月近年来,国内部分学者开展了大量的地震波阻抗反演研究㊂段中钰等[1]利用神经网络中的条件生成式对抗网络对密度曲线进行重构,生成了更符合实际情况的测井曲线㊂梁建刚等[2]通过对含矿最敏感的自然伽马与声波曲线重构波阻抗,发现储层对含铀信息的识别度大幅提升,弥补了传统反演对含铀砂体识别上的缺陷㊂张永健等[3]通过对复杂砂砾岩体数据的综合分析并对敏感曲线进行重构,建立了一套岩性识别图版,确立了各类岩性的定量解释标准㊂张远银等[4]利用地震约束建模方法基于线性反演去除地震资料的A V O效应获得纯纵波数据,然后结合测井数据计算波阻抗,以此来提高波阻抗反演质量㊂孔省吾等[5]开展了岩性敏感曲线重构和波形指示反演相结合的砂体预测研究,利用波形横向变化代替传统变差函数进行高频成分模拟来提高反演精度㊂李勇根等[6]提出了一种以多种地震属性加权为反演约束条件的薄储层及烃类高精度预测方法,通过改变模型的复杂度进行反演测试,能够有效判识薄层砂岩含气性及其饱和程度㊂陈彦虎等[7]通过地震波形动态聚类分析,建立地震波形结构与高频测井曲线结构的映射关系,提高了反演结果的分辨率;通过构建不同地震相类型的贝叶斯反演框架,实现了真正意义上的相控反演㊂丘东气田西山窑组的砂体横向变化快,砂体展布不清,且单砂体较薄,常规的波阻抗反演不能很好地识别出薄层砂岩㊂为此,笔者对丘东气田西山窑组的砂体进行曲线重构的波形指示反演[8-10],目的在于识别砂泥岩,查明砂体的分布规律,为寻找岩性气藏奠定基础㊂1波形指示反演技术1.1基本原理薄互层的 调谐效应 形成复波的反射界面,蕴含了一系列薄互层的综合地质信息,表现为反射系数结构和分布(反射系数的大小值域㊁垂向个数㊁横向分布和间距等),也就是说这些因素影响了㊁甚至决定了地震反射波形的调谐特征㊂在相似的沉积条件下,包括相似的沉积环境㊁相似的砂泥组合特征,反射系数结构和分布与地震波形可以形成内在匹配的关系㊂基于传统的地质统计学反演方法[11],通过地震反射波形的组合特征来优化反射系数的结构和分布,代替传统的变差函数分析,通过钻井和地震联合反演的手段,对储层的空间变化和分布进行模拟和表征,以达到提高地震反演分辨率的目的[12]㊂换句话说,由于地震资料横向分辨率高,地震反射波形变化定性地反映了沉积物的岩性变化特征,因此利用地震波形的横向变化进行储层横向表征,相较于变差函数分析,更能客观地反映岩性的变化[13]㊂在对无井区进行随机模拟时,变差函数是根据已钻井样本统计模拟分析,得到薄互层的概率分布㊁变程参数值,强调的是岩性空间变异程度和距离变化的因素;而波形指示反演强调的是波形组合的相似性和空间距离两个因素,突出了地震波形的横向变化因素,在空间上体现了岩性组合和沉积相带的约束[14](见图1)㊂1.2研究思路首先,在曲线重构的基础上,基于测井曲线和地震波形组合特征,优选钻井岩性变化和波形特征匹配度高的初始样本,建立初始的样本模型[15]㊂其次,根据砂泥岩岩石物理交会分析得到先验概率分布函数,通过地震波阻抗与初始样本模型进行匹配滤波处理,求得似然函数㊂例如,研究区若干口钻井的地震相特征类似,视其为沉积环境和砂泥组合类似㊂这种情况下,高频成分可能不同,但是低频有共性,其中测井曲线的共性频带范围超出了地震数据的有效频带㊂基于此,低频段反演结果的确定性增加,也约束了高频段范围,降低了反演的随机性㊂最后,在贝叶斯理论框架下,将先验概率分布函数和似然函数分布视为目标函数,通过模型参数的随机扰动,运用M-H抽样算法进行分布抽样,得到后验概率分布函数,其最大解得以有效实现,将若干次有效实现的随机解视为期望值输出:F(x0)= N i=1λi F(x i)式中:F(x0)为预测未知点的值;N为样本点的个数;λi为第i点处已知样本点对预测未知点的权重;Copyright©博看网. All Rights Reserved.F (x i)为样本井点处的地震波形优选的值㊂ 注:红色为参与反演井;紫色为验证井㊂图1 变差函数与波形指示反演的选择统计样本示意图F i g .1 S c h e m a t i c d i a g r a mo f s e l e c t e d s t a t i s t i c a l s a m pl e s f o r v a r i a t i o n f u n c t i o na n d w a v e f o r mi n d i c a t i o n i n v e r s i o n图2 波形指示反演流程图F i g.2 F l o wc h a r t o f t h ew a v e f o r mi n d i c a t i o n i n v e r s i o n 波形指示反演流程如图2所示㊂2 应用实例2.1 工区概况丘东气田位于新疆维吾尔自治区吐鲁番地区鄯善县境的七克台镇北东方向约40k m 处㊂南临温吉桑油气田的温八区块,西与米登油田呈斜鞍相接,北靠博格达山㊂地表为戈壁滩,地形总体平坦,呈北高南低趋势㊂地面海拔506~607m ㊂丘东背斜位于吐鲁番坳陷内的台北凹陷中部㊁温吉桑构造带东北端,北邻丘东次凹,南以丘东南断层与温八构造相接,西与米登油田呈斜鞍相连㊂丘东气田所处的台北凹陷面积大,巨厚的侏罗系是吐哈盆地主力生烃凹陷,也是盆内油气最富集的构造单元㊂丘东气田西山窑组发育一套以扇三角洲前缘为主的灰色砂岩,为区内主要的含气层段,砂体横向变化快且砂体展布不清㊂2.2 曲线敏感性分析及曲线重构敏感曲线重构的关键是对岩性变化(尤其是薄互层)敏感的测井曲线进行分频重构,以此提高薄层识别能力,这是高分辨率反演的重要突破口㊂不同类型的地球物理测井数据从不同角度反映了地层的岩石物理性质,不同类型的测井曲线相互关联是测井曲线重构反演的理论基础㊂测井曲线敏感性分析是重构反演的关键点㊂对测井曲线预处理后进行敏感性分析,分别对G R (自然伽马)拟声波波阻抗(Z G R )㊁S P (自然电位)拟声波波阻抗(Z S P )㊁R L L D (深侧向电阻率)拟声波波阻抗(Z R L L D )与A C (声波时差)进行交会㊂研究区共有71口井,其中40口井有岩心和测井数据,可作为关键井以测井曲线为基础开展曲线交会和储层敏感性分析㊂从Z G R 与A C 的交会图(见图3(a))上可知,有效砂岩主要集中在高波阻抗区域,非储层主要集中在低波阻抗区域,两者叠置区域较多,难以明显区分砂泥岩㊂从Z S P 与A C 的交会图(见图3(b ))上可知,有效砂岩主要集中在高波阻抗区域,非储层主要集中在低波阻抗区域,两者叠置区域较多,难以明显㊃13㊃第20卷第2期谢锐杰等:丘东气田西山窑组波形指示反演砂体识别Copyright ©博看网. All Rights Reserved.㊃23㊃长江大学学报(自然科学版)2023年3月区分砂泥岩㊂从Z R L L D与A C的交会图(见图3(c))上可知,有效砂岩主要集中在高波阻抗区域,非储层主要集中在低波阻抗区域,两者叠置区域少,可以明显区分砂泥岩,分界线在8000(g㊃m)/(s㊃c m3)㊂图3Z G R㊁Z S P㊁Z R L L D与A C交会图F i g.3C r o s sm a p o f ZG R,Z S P,Z R L L D a n dA C图4是R L L D曲线重构前后的波阻抗直方图,重构后的波阻抗不仅对岩性的敏感性增强,其差异性也有一定的提升,理论上满足了较薄砂岩预测的技术需求,更容易刻画出砂泥岩㊂综上所述,本次波形指示反演采用基于R L L D曲线重构拟A C和D E N(密度)测井数据㊂2.3合成地震记录制作在合成地震记录的制作过程中发现,地震剖面与子波的极性均为正极性,所以子波先后选取零相位的雷克子波(主频选取26H z)㊁井旁道子波㊂丘东气田利用上述方法共完成了40口井的地震合成记录(见图5),根据这些合成地震记录可以对标准反射层及油层㊁储层等进行精细标记㊂标定结果从波组关系到能量均具有较好的相似性,平均相关系数达到85%㊂同时为了更精确地判断砂体层位在地震剖面上的反射响应,使用经过放大的地震剖面与合成记录进行对比标定,从而能在单井上明确研究区内各界面的反射特征㊂2.4初始模型建立初始模型的建立过程就是井震相结合的过程,将井上的低频信息与地震频带宽度相结合,通过波阻抗曲线的内插外推来给反演结果提供低频分量㊂具体做法是先进行等时地震层位的追踪,在等时层位之间按等比例内插产生内插层位,选择适合于研究区地质情况的插值方法,在测井曲线和层位的共同约束下进行内插外推,建立初始模型㊂检查初始模型是否符合地质要求,如果不符合则需更改参数进行模型Copyright©博看网. All Rights Reserved.图4 R L L D 曲线重构前、后的波阻抗直方图对比F i g .4 C o m p a r i s o no f t h ew a v e i m p e d a n c e h i s t o gr a m s b e f o r e a n da f t e r t h eR L L Dc u r v e r e c o n s t r u c t i o n 注:X 11,X 21,X 31,X 41均为砂体的层位符号㊂图5 Q D 27井合成地震记录F i g .5 S y n t h e t i c s e i s m o gr a mo fw e l l Q D 27的重建,直到满意为止㊂初始模型精度受地震采样率㊁地震层位㊁标定精度和低频波阻抗模型的影响㊂以目的层顶底界面为约束层位,等比例内插层位,在测井曲线和层位的共同约束下进行内插外推,建立目的层的低频模型㊂2.5 反演关键参数设置2.5.1 平滑半径该参数表示储层的横向变化程度,其值越大,说明横向变化越快,反之越慢㊂通过实验对比分析,本次反演平滑半径定义为1㊂2.5.2 有效样本数该参数表征地震波形空间变化对储层的影响程度㊂由于通过井震统计的样本数逐渐增加,其地震相关性也会上升,当数据达到一定数量之后,其相关性趋于平衡不再上升,此时称为有效样本数量㊂结合工区资料,本次反演选取样本数为5㊂2.5.3 频率参数该参数控制了反演结果的有效频带,进而影响了反演的分辨率㊂结合工区实际情况,选择过Q D 6井的任意测线进行试验分析㊂结果表明,随着高频截止频率的上升,反演后得到的分辨率也逐渐上升;但当截止频率大于150H z 时,反演后产生的随机性变高,且钻井岩性匹配度变差,因此最终选取150H z作为最佳高频截止频率,兼顾了薄互层的识别精度,也保证了反演质量(见图6)㊂2.6 反演效果分析在上述反演参数对比分析的基础上,对重构后的波阻抗曲线开展反演工作,得到砂泥岩地震反演体㊂图7为过Q D 27井的波阻抗反演剖面,可以看出,砂岩的垂向叠置关系清楚,砂体与钻井剖面砂体符合程度较高,能识别5m 左右的单砂体,井间砂体变化快,识别特征明显,反演精度较高㊂反演结果发现,在井点位置电阻率反演结果与测井曲线符合程度较高,具有明显的优势,忠实于地震资料,垂向视分辨率明显提高,能够有效识别研究区较薄砂体㊂㊃33㊃第20卷第2期谢锐杰等:丘东气田西山窑组波形指示反演砂体识别Copyright ©博看网. All Rights Reserved.图6 井震相关性分析F i g .6 W e l l -s h o c k c o r r e l a t i o na n a l ys i s 在反演体中沿着目的层段提取电阻率均方根振幅属性,如图8所示,南部到中部地区为黄色红色,为中高波阻抗值,连片分布的范围相对较大,其余地区波阻抗为绿色蓝色,为低波阻抗值,可知研究区物源主要来自南部方向㊂将已钻井在目的砂体层位(X 11,X 21+X 22,X 41)的钻遇厚度与反演厚度进行对比,计算其相对误差,相对误差小于等于15%为符合储层预测精度要求,反之则不符合㊂统计17口井的X 11㊁X 21+X 22㊁X 41砂体共计51层的相对误差,小于等于15%的有45层,大于15%的有6层,符合度达88.2%(见表1),且在平面展布特征上符合研究区的砂体分布规律㊂综上所述,波形指示反演方法在研究区中取得了良好的应用效果㊂图7 过Q D 27井波阻抗反演剖面F i g .7 W a v e i m pe d a n c e i n v e r s i o n s e c t i o no fw e l l Q D 273 结论1)通过岩石物理分析,得出西山窑组泥岩波阻抗小于砂岩波阻抗,电阻率曲线R L L D 对岩性的识别敏感度较大,选用R L L D 曲线重构的拟声波阻抗曲线进行波阻抗反演,可以有效识别砂泥岩㊂㊃43㊃长江大学学报(自然科学版)2023年3月Copyright ©博看网. All Rights Reserved.图8 研究区目的层段电阻率均方根振幅平面图F i g .8 P l a no fR M S a m p l i t u d e o f t h e r e s i s t i v i t y o f t h e t a r g e t l a y e r s e g m e n t i n t h e s t u d y ar e a 表1 已钻井砂体厚度与预测砂体厚度对比T a b l e 1 C o m p a r i s o no f d r i l l e d s a n db o d y t h i c k n e s sw i t h p r e d i c t e d s a n db o d yt h i c k n e s s 井名X 11钻遇厚度/m 反演厚度/m 相对误差/%X 21+X 22钻遇厚度/m 反演厚度/m 相对误差/%X 41钻遇厚度/m 反演厚度/m 相对误差/%D S 2井14.0014.000.1615.6013.9010.9062.0057.806.71J 4H 井4.004.001.4511.3012.6011.5056.9055.202.87M 3井2.602.600.359.9011.0011.2855.4055.800.62Q D 22井20.2020.200.1015.6014.308.1756.4053.205.69Q D 23井15.2015.100.5818.9016.3013.8040.4040.000.95Q D 25井20.2020.501.4920.4017.6013.6544.9046.102.54Q D 27井4.504.704.6423.8020.3014.7038.8038.800.18Q D 3井8.508.302.008.8012.0036.7547.3053.3012.78Q D 31井40.1040.000.0720.9019.008.8152.8049.905.53Q D 32井23.6023.800.592.808.00185.7249.6051.103.00Q D 36井6.306.102.567.104.6034.7449.6051.804.35Q D 48井12.8012.105.7917.5017.400.6070.0059.7014.75Q D 52井21.8022.201.887.0011.9070.5750.2053.706.93Q D 58井26.6026.600.103.905.2032.3754.5051.505.34Q D 59H 井9.208.803.6028.5024.5014.1045.0044.900.29Q D 7井19.7019.500.997.0013.3089.8243.4046.406.81Q D 9井29.6029.500.3021.3021.902.5641.3037.409.332)利用基于敏感曲线重构的波形指示反演技术,开展了丘东气田西山窑组波阻抗反演和储层预测,反演结果在垂向上对薄砂体的识别能力大大增强,反演分辨率高,可以刻画出厚度5m 左右的单砂体,满足了储层预测的要求㊂3)研究区南部波阻抗值较大,砂岩发育,其余地区波阻抗值较小,从南往北砂体逐渐减薄,物源主要来自南部方向㊂㊃53㊃第20卷第2期谢锐杰等:丘东气田西山窑组波形指示反演砂体识别Copyright ©博看网. All Rights Reserved.㊃63㊃长江大学学报(自然科学版)2023年3月参考文献:[1]段中钰,吴俣,肖勇,等.基于C G A N与C N N-G R U组合模型的密度测井曲线重构方法[J].地球物理学进展,2022,37(5): 1941-1945.D U A NZY,WU Y,X I A O Y,e t a l.D e n s i t y l o g g i n g c u r v e r e c o n s t r u c t i o nm e t h o db a s e do nC G A Na n dC N N-G R Uc o m b i n e dm o d e l[J]. P r o g r e s s i nG e o p h y s i c s,2022,37(5):1941-1945.[2]梁建刚,秦喜林,匡海阳,等.自然伽马曲线重构波阻抗反演在勘探含铀有利成矿砂体中的尝试[J/O L].中国地质, 2020.h t t p s://k n s.c n k i.n e t/k c m s/d e t a i l/11.1167.p.20201124.1301.002.h t m1.L I A N GJG,Q I NXL,K U A N G H Y,e t a l.A n a t t e m p t t o f i n d f a v o r a b l e o r e-f o r m i n g s a n d b o d i e s b y r e c o n s t r u c t e dw a v e i m p e d a n c e i n v e r s i o n u s i n g n a t u r a l g a m m a c u r v e[J/O L].G e o l o g y i nC h i n a,2020.h t t p s://k n s.c n k i.n e t/k c m s/d e t a i l/11.1167.p.20201124.1301.002.h t m1.[3]张永健,汤子余,卢美月,等.基于曲线重构的复杂砂砾岩体岩性识别研究 以北16井区梧桐沟组为例[J].西安石油大学学报(自然科学版),2020,35(2):20-25,34.Z HA N G YJ,T A N GZY,L U M Y,e t a l.L o g g i n g i d e n t i f i c a t i o nm e t h o d o f c o m p l e x g l u t e n i t e b o d i e s b a s e d o n c u r v e r e c o n s t r u c t i o n:Ac a s e s t u d y o nW u t o n g g o u f o r m a t i o n i nB e i16w e l lb l o c k[J].J o u r n a lo fX i a nS h i y o u U n i v e r s i t y(N a t u r a lS c i e n c eE d i t i o n),2020,35(2): 20-25,34.[4]张远银,孙赞东,韩剑发,等.地震约束建模的强非均质碳酸盐岩储层波阻抗反演[J].石油地球物理勘探,2019,54(6): 1316-1323.Z HA N G Y Y,S U NZD,H A NJF,e ta l.I m p e d a n c e i n v e r s i o no f s t r o n g l y h e t e r o g e n e o u sc a r b o n a t er e s e r v o i rb a s e do na m p l i t u d e-c o n s t r a i n e dm ode l i n g[J].O i lG e o p h y s i c a l P r o s p e c t i n g,2019,54(6):1316-1323.[5]孔省吾,张云银,沈正春,等.波形指示反演在灰质发育区薄互层浊积岩预测中的应用 以牛庄洼陷沙三中亚段为例[J].物探与化探,2020,44(3):665-671.K O N GX W,Z HA N G Y Y,S H E N ZC,e ta l.T h ea p p l i c a t i o no fw a v e f o r mi n v e r s i o n p r e d i c t i o no f t h i nt u r b i d i t e r e s e r v o i r t oc a l c a r e o u sd e p o s i t i o n a l a r e a:Ac a s e s t u d y o f E3s23i nN i u z h u a n g S a g[J].G e o p h y s i c a l a n dG e o c h e m i c a l E x p l o r a t i o n,2020,44(3):665-671.[6]李勇根,张朝军,曾庆才,等.薄层砂岩含气性多维信息加权约束地震反演判识方法及应用 以鄂尔多斯盆地庆北地区山西组一段为例[J].地质学报,2021,95(8):2628-2640.L IY G,Z H A N GC J,Z E N GQC,e t a l.I d e n t i f i c a t i o nm e t h o d o f g a s b e a r i n g p r o p e r t y o f t h i n t i g h t s a n d s t o n e b a s e d o n c o m p r e h e n s i v e a n a l y s i s o f g e o l o g y a n d g e o p h y s i c s:Ac a s e s t u d y o n t h e1s tM e m b e r o f S h a n x i F o r m a t i o n i n n o r t hQ i n g c h e n g a r e a o fO r d o s B a s i n[J].A c t aG e o l o g i c a S i n i c a,2021,95(8):2628-2640.[7]陈彦虎,毕建军,邱小斌,等.地震波形指示反演方法及其应用[J].石油勘探与开发,2020,47(6):1149-1158.C H E N Y H,B I J J,Q I UXB,e t a l.Am e t h o d o f s e i s m i cm e m e i n v e r s i o n a n d i t s a p p l i c a t i o n[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, 2020,47(6):1149-1158.[8]曹孟起,王九拴,邵林海.叠前弹性波阻抗反演技术及应用[J].石油地球物理勘探,2006,41(3):323-326.C A O M Q,WA N GJ S,S HA OLH.P r e s t a c k e l a s t i cw a v e i m p e d a n c e i n v e r s i o n t e c h n i q u e a n d a p p l i c a t i o n[J].O i l G e o p h y s i c a l P r o s p e c t i n g, 2006,41(3):323-326.[9]栾颖,冯晅,刘财,等.波阻抗反演技术的研究现状及发展[J].吉林大学学报(地球科学版),2008,38(S1):94-98.L U A N Y,F E N G H,L I UC.T h e r e s e a r c h p r e s e n t a n d f u t u r e o fw a v e i m p e d a n c e i n v e r s i o n t e c h n i q u e[J].J o u r n a l o f J i l i nU n i v e r s i t y(E a r t h S c i e n c eE d i t i o n),2008,38(S1):94-98.[10]李国发,廖前进,王尚旭,等.合成地震记录层位标定若干问题的探讨[J].石油物探,2008,47(2):145-149.L IGF,L I A O QJ,WA N G SX,e t a l.D i s c u s s i o n s a b o u th o r i z o nc a l i b r a t i o nb a s e do nw e l l-l o g s y n t h e t i c s e i s m o g r a m[J].G e o p h y s i c a l P r o s p e c t i n g f o r P e t r o l e u m,2008,47(2):145-149.[11]张丰麒,刘俊州,刘兰锋,等.确定性反演协同约束的叠后随机地震反演方法[J].石油地球物理勘探,2021,56(5):1137-1149.Z H A N GFQ,L I UJ Z,L I ULF,e t a l.T h em e t h o d o l o g y o f a p o s t-s t a c k s t o c h a s t i c s e i s m i c i n v e r s i o nw i t h t h e c o-c o n s t r a i n t o f d e t e r m i n i s t i ci n v e r s i o n[J].O i l G e o p h y s i c a l P r o s p e c t i n g,2021,56(5):1137-1149.[12]王玉学,汤达祯,龚建华,等.拟声波反演技术在储层预测中的应用 以吐哈盆地台北凹陷胜北地区白垩系为例[J].石油实验地质,2004,26(3):303-307.WA N G Y X,T A N GDZ,G O N GJH,e t a l.A p p l i c a t i o n o f q u a s i-a c o u s t i c i n v e r s i o n t e c h n o l o g y t o r e s e r v o i r p r e d i c t i o n:Ac a s e s t u d y o f t h e c r e t a c e o u s i nS h e n g b e i a r e a o f t h eT a i b e iS a g,t h eT u h aB a s i n[J].P e t r o l e u m G e o l o g y a n dE x p e r i m e n t,2004,26(3):303-307.[13]尹正武,陈超,彭嫦姿.拟声波反演技术在优质泥页岩储层预测中的应用 以焦石坝页岩气田为例[J].天然气工业,2014,34(12):33-37.Y I NZ W,C H E N C,P E N G CZ.A p p l i c a t i o no f p s e u d o-a c o u s t i c i m p e d a n c e i n v e r s i o nt o q u a l i t y s h a l er e s e r v o i r p r e d i c t i o n:Ac a s e s t u d y o f t h e J i a o s h i b a s h a l e g a s f i e l d,S i c h u a nB a s i n[J].N a t u r a lG a s I n d u s t r y,2014,34(12):33-37.[14]廉桂辉,朱亚婷,王晓光,等.叠前反演技术在玛湖油田储层预测中的应用[J].特种油气藏,2022,29(1):80-84.L I A NG H,Z H U Y T,W A N G X G,e t a l.A p p l i c a t i o no f p r e-s t a c k i n v e r s i o nt e c h n o l o g y t or e s e r v o i r p r e d i c t i o no fM a h u O i l f i e l d[J].S p e c i a lO i l&G a sR e s e r v o i r s,2022,29(1):80-84.[15]张学芳,董月昌,慎国强,等.曲线重构技术在测井约束反演中的应用[J].石油勘探与开发,2005,32(3):70-72.Z H A N G XF,D O N G YC,S H E N G Q,e t a l.A p p l i c a t i o no f l o g r e b u i l d i n g t e c h n i q u e i n c o n s t r a i n i n v e r s i o n[J].P e t r o l e u m E x p l o r a t i o na n dD e v e l o p m e n t,2005,32(3):70-72.[编辑]龚丹Copyright©博看网. All Rights Reserved.。
(2)可控震源。
其激发信号作用时间较长、且为均衡振幅的连续扫描振动信号，例如双轴井下震动器Z-Trac。
特点：Z-Trac震源，可同时激发高频纵波、横波。
激发能量频带宽，且能稳定输出最低有效频率在30Hz，仪器稳定，保养周期长，单次保养，能保证100000万次激发正常，不对井壁及环空水泥产生影响。
1.2.2 检波器选择(1)基于流体耦合的多级水听器拖缆。
其特点为：布置简单、频率响应高，但缺乏推靠式传感器那种矢量波场的测量能力，受管波的影响严重。
(2)推靠式三分量多级检波器。
其特点：推靠臂工作，使检波器与地层耦合良好，记录井壁粒子运动形成的矢量波场，有抑制管波，但其记录效率比水听器拖缆差。
对井间地震来说，井下的干扰因素众多，波场复杂，在采集中首先要保证有效波场的高信噪比；另外，井间地震采集，不但有一般构造与地层的空间展布关系，而且还有各向异性与流体的问题，在施工时，应该综合考虑区块的地质情况、井况、施工要求、预算等多种因素，采用性价比最高、最适合该区块的井下检波器。
1.2.3 采集方法井间地震采集方法，最为常见的有两种：共检波器采集(CRG)及共炮点采集(CSG)。
以CRG为例，检波器固定深度不动，震源从底深度，定点上提至顶深度，该采集就构成了一个扇区。
通常情况下，检波器需要上移几个深度，每移动并固定到一个深度，则震源移动并激发，从而再次采集一个扇区。
以这样的模式，采集完所有扇区，则完成井间地震采集。
同理，CSG则相当于震源不动，检波器移动采集扇区，或者说相当于一个震源炮点的非零VSP采集。
1.2.4 采集参数设置(1)检波器、震源移动深度选择。
在施工之前，必须对施工井进行模拟实验，以确定震源、检波器的布放深度，必须保证目的层必须在覆盖范围内。
(2)采集深度间隔。
一般检波器级间距选用15m(或30m)，采集深度间隔为2.5m、5m、7.5m、15m(通常选择5m深度间隔)，震源采集深度间隔一般为5m。
