岩性油气藏  2020, Vol. 32 Issue (3): 115-121       PDF    
×
地震伪井速度点宏观校正方法与应用——以珠江口盆地M气田为例
罗泽, 谢明英, 梁杰, 涂志勇, 侯凯     
中海石油(中国)有限公司 深圳分公司, 广东 深圳 518000
摘要: 地震勘探开发中速度研究方法的选取是时深转换的关键。鉴于珠江口盆地M气田断裂系统复杂,构造主体部位受断层影响导致横向速度变化明显,而以往时深转换方法诸如变速成图法或井时深关系拟合法得到的深度构造图误差偏大,很难满足开发阶段高精度构造研究的需求。由此在传统变速成图方法的基础上进行了改进,优化后的地震伪井速度点宏观校正方法是利用伪井速度时间函数与测井速度时间函数的差值对地震速度进行校正,将井纵向分辨率高和地震横向分辨率高的优势更好地结合在一起,从而提高了工作效率,同时也降低了构造的不确定性。新提出的地震伪井速度点宏观校正方法有效提高了该区时深转换的精度。
关键词: 时深转换    变速成图    伪井速度点    宏观校正    珠江口盆地    
Macro-correction method and application of seismic pseudo-well velocity point: a case study from M gas field in Pearl River Mouth Basin
LUO Ze, XIE Mingying, LIANG Jie, TU Zhiyong, HOU Kai     
Shenzhen Branch of CNOOC Ltd., Shenzhen 518000, Guangdong, China
Abstract: In seismic exploration and development, the selection of velocity research method is the key to timedepth conversion. In view of the complexity of fault system in M gas field of Pearl River Mouth Basin, the transverse velocity variation is obvious due to the influence of faults on the main part of structure, while the error of depth structure map obtained by previous time-depth conversion methods, such as variable velocity mapping or well-time-depth relationship fitting, is obvious. It is too large to meet the needs of high-precision structural research in the development stage. The traditional variable speed mapping method was improved. The optimized pseudo-well velocity point macro-correction method is to use the difference between pseudo-well velocity time function and logging velocity time function to correct the seismic velocity, which combines the advantages of high vertical resolution of wells and high lateral resolution of earthquakes better. It can improve work efficiency and reduce structural uncertainty at the same time. The actual drilling results show that the macro-correction method of pseudo-well velocity points effectively improves the accuracy of time-depth conversion in this area.
Key words: time-depth conversion    variable speed mapping    macro-correction of pseudo-well velocity point    Pearl River Mouth Basin    
0 引言

地震资料解释是油田勘探开发的一项重要工作,而时深转换是将时间域地震资料解释成果转化成构造图的桥梁。时深转换方法越精确,最后得到的构造图越能最真实地还原地下实际地质构造[1-2]。关于时深转换方法,国内学者[3-6]作了大量研究,探索出多种时深转换方法,如变速成图法、井时深关系拟合法、地质构造约束层速度模型法、地震速度场法和分层统一速度成图法等。这些常规时深转换方法在速度横向变化小的区域效果较好[7-10],但对于速度场非均质性较强的区域,这些方法却具有其局限性,且会导致最终构造图误差较大。为此,Slotnick[11]和Keydar等[12]提出了地震走时计算方法,由于底层速度随深度是线性变化的($v = {v_o} + az $),所以这些方法对构造图精度提高虽具有一定效果[9],但在开发阶段对构造的精细要求仍然存在一定差距。

时深转换的核心是速度的求取,速度一般来源于叠加速度和钻井速度。叠加速度具有全区均匀分布且横向分辨率高的特点[13-14],钻井速度具有纵向分辨率高且准确可靠的特点[15-17],若能将二者紧密结合起来,实现地震和测井速度协同研究,预计能得到三维空间高精度速度场[18]。鉴于此,凌云等[19-20]和梁卫等[21]提出了井震结合的速度研究方法,通过模型正演和实际研究区证实,井震结合时深转换方法能够较准确地反映实际地质构造的时深关系,最终得到精度较高的构造图,后续开发方案实钻井也证实了这种井震结合法得到的构造图误差较小且精度较高,大部分实钻井误差保持在2 m左右,对复杂地质环境的低幅度构造也有很好的应用效果[21-22]。在此基础上,针对井震结合时深转换方法展开进一步研究,对井震结合法进行优化和改进,并提出地震伪井速度点宏观校正时深转换方法,以期得到更高精度的构造图。

1 地震伪井速度点宏观校正时深转换方法

目前,比较常用且精度较高的时深转换方法是变速成图法,但是其对速度横向变化大的区域构造成图效果不佳,原因是速度不仅受埋藏深度的影响,还受到岩性、沉积相带等特殊地质因素的影响。对于地质条件复杂,速度横向非均质性强的地区,地震伪井速度点宏观校正是行之有效的时深转换方法,其主要分为层速度求取、平均速度求取和速度宏观趋势校正三大步骤。

(1) 层速度求取,利用高精度速度拾取和沿层速度分析。具体做法为先采用地震叠加剖面拾取宏观速度场,然后在参考层约束下沿层重新提取叠加速度和偏移速度等地震速度。

(2) 平均速度求取。先追踪出等时界面,确定层位倾角,再进行模型层析法处理,即在已知第n层层速度和反射界面的基础上,利用射线追踪原理迭代求取第n+1层的层速度并确定第n+1层反射界面,依次循环类推,逐层计算得到地震速度场。

(3) 速度宏观趋势校正。在对时间域地震资料和深度域测井资料进行标准参考层解释的基础上,假定时间域地震相对等时界面与相对应的深度域测井层位满足时深转换速度的概率统计分布特征[21],当ε趋于最小值时,可以求取时深转换中速度的拟合系数V0,即

$ \begin{array}{l} \varepsilon = \sum\limits_{j = 1}^M {\sum\limits_{i = 1}^N {} [(\frac{{2{z_i}}}{{{t_j}}}} ) - {v_0}(1 + {\beta _i}{z_i})]\\ i = 1,2, \ldots ,N;j = 1,2, \ldots ,M \end{array} $ (1)

式中:Zi为测井层位的深度,m;i为测井的点数,总点数为Ntj为井旁地震层位的双程旅行时,ms;j为参考层的数量,总层数为Mvo为顶面速度,m/s;βi为随深度变化而变化的系数。

利用式(1)可以求取最佳平均速度。地震伪井速度点宏观校正时深转换方法即是基于此理论求取的最佳时深转换速度,技术流程为:①由地震处理叠加速度体计算平均速度体;②利用解释的时间层位沿平均速度体提取层平均速度,再沿测井轨迹提取伪井速度,然后利用伪井速度时间函数与测井速度时间函数的差值对地震速度进行校正,再由校正的地震速度进行时深转换得到初始的深度域构造图;③结合地质认识,利用井点误差拟合一个横向上具有地质意义的误差趋势面;④将这个误差趋势面加到初始的深度域构造图上得到最终的高精度深度构造图(图 1)。

下载原图 图 1 地震伪井速度点宏观校正时深转换方法技术流程 Fig. 1 Technical flow chart of time-depth conversion method for macro-correction of pseudowell velocity point

新提出的地震伪井速度点宏观校正方法是在井震宏观速度校正方法的基础上进行改进而得到的。井震宏观速度校正方法是用井速度宏观校正地震速度网格,校正量为沿层地震速度时间函数与测井速度时间函数之间的差值。优化后的地震伪井速度点宏观校正方法是综合伪井速度与测井速度对地震速度进行校正,采用伪井速度时间函数与测井速度时间函数的差值对地震速度进行校正,将井纵向分辨率高和地震横向分辨率高的优势更好地结合在一起,在提高工作效率的同时也能降低构造的不确定性。

2 应用效果

珠江口盆地M气田是一个发育在基底隆起之上的翘倾半背斜气田,构造主体部位受雁行式排列正断层所控制,速度的横向变化受断层影响明显。区域研究认为,造成这种现象的机理主要为受断层挤压力学作用影响,沿垂直于断层方向存在明显的速度变化趋势,从而导致速度的横向变化受断层的影响明显。图 2是M气田时间域地震双程旅行时层位与深度域测井层位之间换算的各个井点的沿层平均速度散点图及井震拟合趋势线,黑色散点为伪井速度,蓝色散点为钻井实测VSP速度,浅绿色散点为进行地震伪井速度点宏观校正后得到的沿层平均速度,从图中可明显看出,地震速度与测井速度之间存在差异,通常这二者之间的速度差异是由地层各向异性引起的[23],这种差异在地震速度和测井速度之间是普遍存在的。鉴于此,进行地震伪井速度点宏观校正是十分必要的。

下载原图 图 2 珠江口盆地M气田速度场对比 Fig. 2 Velocity field contrast chart of M gas field in Pearl River Mouth Basin

针对M气田主力层H40层利用常规变速成图、井震宏观速度校正和地震伪井速度点宏观校正等3种方法得到沿层平均速度平面对比图(图 3),分析图 3发现:三者形态大体相似,但井震宏观速度校正和地震伪井速度点宏观校正有测井速度参与校正,在数值上多一个低频分量,反映的地质信息更多,相对于常规变速成图得到的平均速度更准确;而地震伪井速度点宏观校正方法能保证在井点处将地震速度精确校正到测井速度,校正后得到的速度更接近真实地层速度,是对井震宏观速度校正方法的优化和改进。利用3种方法得到层平均速度后,利用时间构造图求取得到初始的深度构造图,提取的勘探井的误差如表 1所列。从表 1可看出,勘探井的最大误差分别为93.30 m,44.17 m和29.55 m等,误差逐步减小。综合对比可知,地震伪井速度点宏观校正方法得到的构造图精度最高且效果最好。

下载原图 图 3 珠江口盆地M气田H40层不同方法求取平均速度平面对比图 Fig. 3 Comparisons of average velocity by different methods of H40 layer in M gas field of Pearl River Mouth Basin
下载CSV 表 1 珠江口盆地M气田H40层误差校正前探井误差统计 Table 1 Exploration well errors before correction of H40 layer in M gas field of Pearl River Mouth Basin

通过对比以上时深转换速度方法发现,地震伪井速度点宏观校正后,精度得到了明显改善,但误差趋势还存在东西分带现象,产生这种现象的原因可能是沉积体在横向分布上存在差异所导致的。为了进一步研究以消除这一现象,在研究区选择一条近东西向的联井地震剖面(图 4),发现研究区东部的T-3井区浅层存在一个厚度很大的低速带。

下载原图 图 4 珠江口盆地M气田联井三维地震剖面 Fig. 4 Three-dimensional seismic profile of combined wells in M gas field of Pearl River Mouth Basin

为了进一步落实误差趋势东西分带现象,研究了海底起伏和东部低速体空间展布情况。通过地震资料对海底、低速体顶和低速体底的精细解释进行追踪,最后得到了研究区的海底深度图和浅层低速体厚度图(图 56)。研究发现,研究区海底起伏和浅层低速体影响存在明显的北西西—南东东向趋势,鉴于这一特征,利用误差面板定量校正减小速度陷阱问题造成的影响,以此提高构造图的精度。具体做法是利用井资料的地质分层数据,结合此特征,最后拟合出一个具有地质意义的误差趋势面(图 7),以实现对构造图误差校正,最终得到深度构造图(图 8)。对比3种时深转换方法得到的构造图结果发现,变速成图方法得到的构造倾角最大,误差最大,构造形态与地质认识较符合;井震宏观速度校正后得到的构造图精度比变速成图高,构造倾角变缓,与地质认识更符合;地震伪井速度点宏观校正法对井震宏观校正进行了优化和改进,得到的构造图误差最小,精度最高,构造形态最平缓,与地质认识最吻合,效果最好。

下载原图 图 5 珠江口盆地M气田海底深度 Fig. 5 Seabed depth map of M gas field in Pearl River Mouth Basin
下载原图 图 6 珠江口盆地M气田浅层低速带厚度 Fig. 6 Thickness map of shallow low velocity zone in M gas field of Pearl River Mouth Basin
下载原图 图 7 地震伪井速度点宏观校正误差趋势面 Fig. 7 Macro-correction error trend chart of pseudo-well velocity point
下载原图 图 8 3种时深转换方法得到的珠江口盆地M气田H40层构造对比图 Fig. 8 Structural contrast before and after macro-correction of pseudo-well velocity point of H40 layer in M gas field of Pearl River Mouth Basin

基于珠江口盆地M气田的构造研究工作,后期在珠江口盆地M气田实施了4口生产开发井(k1— k4),其平面分布位置如图 9所示,实钻误差统计对比结果见表 2。从表 2可知:井震宏观速度趋势校正前最大误差高达14.28 m,误差普遍在10 m左右;井震宏观速度趋势校正后最大误差为1.37 m,误差控制在2 m以内,精度具有一定提高;地震伪井速度点宏观校正后最大误差仅有0.89 m,误差控制在1 m以内,相对比而言,误差最小,精度最高。综上可知,地震伪井速度点宏观校正后构造精度最高,总体上地震伪井速度点宏观校正方法得到的构造图误差最小,精度最高。

下载原图 图 9 珠江口盆地M气田H40层井点分布 Fig. 9 Distribution of well points of H40 Layer in M gas field of Pearl River Mouth Basin
下载CSV 表 2 珠江口盆地M气田H40层地震伪井速度点宏观校正前后开发井预测误差统计 Table 2 Development well prediction errors before and after macro-correction of pseudo-well velocity point of H40 layer in M gas field of Pearl River Mouth Basin
3 结论

(1) 地震伪井速度点宏观校正法是在井震精细标定和层位精细解释的基础上,以多口测井速度数据为基础,在宏观上校正地震速度,既利用了地震速度的横向高分辨率,又结合了测井速度的垂向高分辨率,得到的速度更接近地下真实地层速度。经过综合评价,证实了地震伪井速度点宏观校正时深转换方法对珠江口盆地M气田是适用的。

(2) 通过对珠江口盆地M气田构造的精细研究,并对传统变速成图时深转换方法进行了改进和完善,提出地震伪井速度点宏观校正方法,得到的构造图精度高,误差小。通过地震伪井速度点宏观校正方法可以得到精度较高的构造图,减小实钻深度与预测深度的误差,得到的高精度构造图很好地解决了开发阶段对构造的精细要求,为珠江口盆地M气田的高效开发提供了一定的指导作用。

(3) 实践证明地震伪井速度点宏观校正方法是一种行之有效的时深转换方法。此方法能为南海东部其他区域的油气田精细构造研究提供借鉴和帮助。

参考文献
[1]
孙希杰. 巨厚松散层煤田三维地震勘探时深转换方法探讨. 中国矿业, 2018, 27(1): 412-414.
SUN X J. Probe into coalfield 3D seismic prospecting timedepth conversion under thick soil layer. China Mining Magazine, 2018, 27(1): 412-414.
[2]
李军. 三维地震资料解释中时深转换方法对比分析及应用. 海油石油, 2014, 34(1): 36-40.
LI J. Analysis and applications of time-depth transform method in 3D seismic interpretation. Offshore Oil, 2014, 34(1): 36-40.
[3]
徐立恒, 鲜波, 薛玉英, 等. 高精度地震时深转换方法研究及应用. 吉林大学学报(地球科学版), 2014, 44(5): 1712-1718.
XU L H, XIAN B, XUE Y Y, et al. Study and application on seismic time-depth conversion with high-precision. Journal of Jilin University (Earth Science Edition), 2014, 44(5): 1712-1718.
[4]
陈可洋. 逆时成像技术在大庆探区复杂构造成像中的应用. 岩性油气藏, 2017, 29(6): 91-100.
CHEN K Y. Application of reverse-time migration technology to complex structural imaging in Daqing exploration area. Lithologic Reservoirs, 2017, 29(6): 91-100. DOI:10.3969/j.issn.1673-8926.2017.06.012
[5]
郭璃, 邓勇, 赵顺兰, 等. 特殊地质条件下的时深转换方法探讨:以涠西南凹陷南部斜坡带Z井为例. 地球物理学进展, 2017, 32(1): 152-158.
GUO L, DENG Y, ZHAO S L, et al. Research and discussion of time-depth conversion methods under special geological conditions. Progress in Geophysics, 2017, 32(1): 152-158.
[6]
李培培, 刘志国, 杨松岭, 等. 虚拟井技术在无井或少井条件下时深转换中的应用. 物探与化探, 2015, 39(5): 994-1000.
LI P P, LIU Z G, YANG S L, et al. The application of virtual well technique to time-depth conversion under the condition of no or few wells. Geophysical and Geochemical Exploration, 2015, 39(5): 994-1000.
[7]
刘文卿, 王孝, 胡书华, 等. 测井与全方位道集联合各向异性参数建模及成像. 岩性油气藏, 2018, 30(6): 83-88.
LIU W Q, WANG X, HU S H, et al. Well constrained anisotropic velocity model building based on full-azimuth angle gathers and imaging. Lithologic Reservoirs, 2018, 30(6): 83-88.
[8]
李文静, 王英民, 何敏, 等. 珠江口盆地中中新世陆架边缘三角洲的类型及控制因素. 岩性油气藏, 2018, 30(2): 58-66.
LI W J, WANG Y M, HE M, et al. Types and controlling factors of shelf margin delta of Middle Miocene in Pearl River Mouth Basin. Lithologic Reservoirs, 2018, 30(2): 58-66.
[9]
韩强, 杨子川, 赵渊. 塔里木盆地轮台地区低幅度构造圈闭落实技术及其应用. 石油与天然气地质, 2010, 31(1): 43-48.
HAN Q, YANG Z C, ZHAO Y. Identification of low-amplitude structural traps in Luntai area, the Tarim Basin and its application. Oil & Gas Geology, 2010, 31(1): 43-48.
[10]
蔡刚, 屈志毅. 构造复杂地区地震资料速度和成图方法研究与应用. 天然气地球科学, 2005, 16(2): 246-249.
CAI G, QU Z Y. Research and application of seismic data velocity and mapping method in complex structural areas. Natural Gas Geoscience, 2005, 16(2): 246-249. DOI:10.3969/j.issn.1672-1926.2005.02.025
[11]
SLOTNICK M M. On seismic computations, with applications Ⅱ. Geophysics, 1936, 1(3): 299-305. DOI:10.1190/1.1437111
[12]
KEYDAR S, KOREN Z, KOSLOFF D, et a1. Optimumtime-todepth conversion. Geophysics, 1989, 54(8): 1001-1005. DOI:10.1190/1.1442724
[13]
张军华, 王静, 郑旭刚, 等. 关于几种速度分析方法的讨论及效果分析. 石油物探, 2009, 48(4): 347-353.
ZHANG J H, WANG J, ZHENG X G, et al. Discussion and result analysis on several velocity analysis methods. Geophysical Prospecting for Petroleum, 2009, 48(4): 347-353. DOI:10.3969/j.issn.1000-1441.2009.04.005
[14]
蔡全升, 胡明毅, 陈孝红, 等. 小型断陷湖盆扇三角洲沉积特征与发育模式:以徐家围子断陷北部沙河子组为例. 岩性油气藏, 2018, 30(1): 86-96.
CAI Q S, HU M Y, CHEN X H, et al. Sedimentary characteristics and development model of fan delta in small faulted basin:a case of Shahezi Formation in northern Xujiaweizi Fault Depression, NE China. Lithologic Reservoirs, 2018, 30(1): 86-96. DOI:10.3969/j.issn.1673-8926.2018.01.009
[15]
刘浩杰, 王延光, 何惺华. 多尺度地球物理资料速度分析. 油气地球物理, 2006, 4(4): 19-22.
LIU H J, WANG Y G, HE X H. Velocity analysis for multi-scale geophysical data. Petroleum Geophysics, 2006, 4(4): 19-22.
[16]
武丽, 李剑峰, 施炜, 等. 巴楚夏河工区的三维速度建模方法. 石油地球物理勘探, 2006, 41(1): 87-92.
WU L, LI J F, SHI W, et al. 3-D velocity model-building method in Bachuxia river work zone. Oil Geophysical Prospecting, 2006, 41(1): 87-92. DOI:10.3321/j.issn:1000-7210.2006.01.017
[17]
杨占龙, 肖冬生, 周隶华, 等. 高分辨率层序格架下的陆相湖盆精细沉积体系研究:以吐哈盆地西缘侏罗系-古近系为例. 岩性油气藏, 2017, 29(5): 1-10.
YANG Z L, XIAO D S, ZHOU L H, et al. Depositional system of lacustrine basins within high-resolution sequence framework:a case of Jurassic to Paleogene in western Turpan-Kumul Basin. Lithologic Reservoirs, 2017, 29(5): 1-10. DOI:10.3969/j.issn.1673-8926.2017.05.001
[18]
王树华, 刘怀山, 张云银, 等. 变速成图方法及应用研究. 中国海洋大学学报, 2004, 34(1): 139-146.
WANG S H, LIU H S, ZHNAG Y Y, et al. Research on variable-velocity structure mapping and its application. Journal of Ocean University of China, 2004, 34(1): 139-146.
[19]
凌云, 郭建民, 郭向宇, 等. 油藏描述中的井震时深转换技术研究. 石油物探, 2011, 50(1): 1-13.
LING Y, GUO J M, GUO X Y, et al. Study on seismic timedepth conversion technique in reservoir description. Geophysical Prospecting for Petroleum, 2011, 50(1): 1-13.
[20]
郭向宇, 凌云, 高军, 等. 井地联合地震勘探技术研究. 石油物探, 2010, 49(5): 438-450.
GUO X Y, LING Y, GAO J, et al. Study on well-to-seismic combined seismic exploration technology. Geophysical Prospecting for Petroleum, 2010, 49(5): 438-450. DOI:10.3969/j.issn.1000-1441.2010.05.003
[21]
梁卫, 李熙盛, 罗东红, 等. 井震时深转换技术在低幅度构造评价中的应用. 中国海上油气, 2014, 26(3): 61-71.
LIANG W, LI X S, LUO D H, et al. Applying a technique of borehole seismic time-depth conversion to the evaluation of low relief structures. China Offshore Oil and Gas, 2014, 26(3): 61-71.
[22]
韩令贺, 胡自多, 冯会元, 等. 井震联合网格层析各向异性速度建模研究及应用. 岩性油气藏, 2018, 30(4): 91-97.
HAN L H, HU Z D, FENG H Y, et al. Grid tomography based on well-to-seismic integration in anisotropic velocity modeling and its application. Lithologic Reservoirs, 2018, 30(4): 91-97.
[23]
李伍志, 王璞珺, 张功成, 等. 珠江口盆地深部基底地层的地震时深转换研究. 地球物理学报, 2011, 54(2): 449-456.
LI W Z, WANG P J, ZHANG G C, et al. Researches on timedepth conversion of deep-seated basal strata of Pearl River Mouth Basin. Chinese Journal of Geophysics, 2011, 54(2): 449-456. DOI:10.3969/j.issn.0001-5733.2011.02.023