|
|
An efficient method for drilling trajectory optimization based on seismological and geological guidance of VSP data during drilling and its application |
Jing-Ya YANG1, Xiang-Wen LI1, Yong-Lei LIU1, Bo XU1, Jiang-Tao GAO1, Mao WANG1, Jiang-Yong WU2, Quan ZHANG1 |
1. Bureau of Geophysical Prospecting ( BGP) Company of CNPC,Zhuozhou 072750,China 2. Exploration and Development Research Institute of Tarim Oilfield Company,Korla 841000,China |
|
|
Abstract There are several Ordovician carbonate fractured-vuggy type oil and gas fields in the Tarim Basin.With the continuous development and production of the reservoirs,the small and medium sized fractured-vuggy systems have become the main targets of capacity development.In the process of exploration and development,due to the influence of complex formation velocity,there is a certain error in seismic migration data,which often leads to the failure of drilling.In this paper,an efficient drilling track optimization method based on VSP data for seismic data processing is proposed.First,the VSP velocity and the new triangular mesh model are used to effectively correct the original seismic migration velocity field,and then the modified anisotropic parameter field is used for quick migration.Secondly,the structural characteristics of the whole strata are analyzed in the new data,the positions of the reservoirs are predicted,and the targets are determined.Finally,according to the engineering and geological conditions,the most economical optimal design of the drilling trajectory is used to drill into the target.The application shows that the relatively good effect is obtained in Tazhong area of Tarim Oilfield,and the imaging quality of seismic data has been improved obviously.After well track adjustment,the rates of drilled reservoir and the directly put on production are increased by 16% and 13.2% respectively.The low cost and effective development plan of the oilfield is effectively promoted and implemented.
|
Received: 22 May 2018
Published: 10 April 2019
|
|
|
|
|
|
The research route of VSP seismic geological guidance technology while drilling a—the general process of VSP seismic geological guidance technology while drilling,including;b—TTI anisotropic parameter field construction process and geological interpretation;c—trajectory adjustment process;the velocity in the process is the velocity of vertical propagation of P-wave;δ:the coefficient of variation (describing the degree of variation of vertical velocity);ε:the intensity of P-wave anisotropy (describing the difference between horizontal and vertical velocity);θ:the dip angle of formation;Φ:the azimuth angle
|
a—the results of elastic wave equation simulation and reverse time migration processing;b—the results obtained by the same processing method as figure a;c—the results obtained by the same processing method as figure a;d—the results obtained by the same processing method as figure a;e—the isotropic velocity field corresponds to figure b,the velocity in the black virtual frame is 98% of the value of the velocity field in figure d;f—the isotropic velocity field corresponds to figure c,the velocity in the black virtual frame is 102% of the value of the velocity field in figure d ">
|
Forward analysis of the effect of velocity errors on "beads" imaging in horizontal formation a—the results of elastic wave equation simulation and reverse time migration processing;b—the results obtained by the same processing method as figure a;c—the results obtained by the same processing method as figure a;d—the results obtained by the same processing method as figure a;e—the isotropic velocity field corresponds to figure b,the velocity in the black virtual frame is 98% of the value of the velocity field in figure d;f—the isotropic velocity field corresponds to figure c,the velocity in the black virtual frame is 102% of the value of the velocity field in figure d
|
|
A forward model for the effect of lateral velocity variation on migration and relocation of beads in shallow strata
|
|
A comparison of the results of forward model stacking depth migration showed in the figure 3
|
|
Comparison of modeling effect between triangular mesh algorithm and Petrel conventional construction a—based triangular mesh algorithm is used to construct a stratigraphic plane grid in the modeling results;b—the conventional structural modeling results using Petrel in the same as the stratigraphic plane grid shown in figure a;c—the constructing A-A’ section based on triangular mesh algorithm;d—conventional structure modeling A-A’ section using Petrel
|
|
Velocity contrast map (curve,profile) and well-seismic calibration profile based on VSP driving processing a—the stacking contrast of VSP velocity,pre-correction seismic velocity and post-correction seismic velocity;b—original seismic velocity field cross-well profile;c—the velocity field cross-well profile modified by VSP data;d—the VSP calibration results based on new migration imaging seismic data
|
|
Seismic migration profile and RMS amplitude attribute contrast map of WELL1 well before and after VSP drive processing a—seismic profile of migration processing data driven by VSP data;b—old seismic data profile;c—reservoir prediction plane pap of migration processing data modified by VSP data;d—reservoir prediction plane map based on old data
|
|
Contrast map of RMS amplitude attributes of data before and after multi-well VSP combined drive processing
|
序号 | 井号 | 年度 | 完钻层位 | 漏失 | 放空 | 酸压/m3 | 层段/m | 漏失量/m3 | 井段/m | 长度/m | 1 | W-511-7X | 2016 | O3l | 5405.94~5415 | 506 | | | | 2 | W-265 | 2016 | O2y | | 20 | | | 590 | 3 | W-43-3X | 2017 | O1-2y2 | 5452.6~5520 | 2994.3 | | | | 4 | W-293 | 2017 | O2y | 6184.4~6249 | 1306 | | | | 5 | W-112-2X | 2017 | O1-2y1 | 6472.6~6545 | 668.8 | 6534.95~6535.57 6528.68~6529.35 | 1.29 | | 6 | W-102-3X | 2017 | O3l | 6718.93~6751 | 625 | 6748.65~6749.9 6750.08~6751 | 2.17 | 633 | 7 | W-7-9X | 2017 | O3l | 5628.31~5813 | 1265.8 | | | | 8 | W-83-7X | 2017 | O3l | 5786~5791 | 73.1 | | | | 9 | W-291-5X | 2017 | O2y | 6270.12~6300.12 | 1035.5 | 6272.28~6273.56 6274.84~6275.38 | 1.82 | | 10 | W-1-1X | 2017 | O3l | 5902.75~5960 | 1123 | 5906.2~5906.8 | 0.6 | | 11 | W-801 | 2017 | O3l | 5850.11~5862 | 607.5 | | | 470 | 12 | W-8-2X | 2017 | O1-2y1 | 6155.63~6195 | 797.09 | | | 240 | 13 | W-296 | 2017 | O2y | | | | | 740 | 17 | W-295 | 2017 | O2y | 6078.32~6095.52 | 0.5 | 6192.47~6192.9 6192.9~6193.01 | 0.54 | | 18 | W-102-6X | 2017 | O3l | 6549.98~6570 | 79.2 | | | 740 | 19 | W-113 | 2017 | O1-2y2 | 6861.85~6880 | 1371 | | | | 20 | W-726-1X | 2018 | O3l | 5557.96~5558 | 1546.9 | | | | 21 | W-83-8X | 2018 | O1-2y1 | 5622~5635 | 766.9 | | | | 22 | W-6-6 | 2018 | O3l | 5782.3~5807 | 7.8 | | | | 23 | W-15-7X | 2018 | O2y | 6091.38~6124.36 | 1246.2 | | | | 24 | W-43-7 | 2018 | O1-2y2 | 5948.15~5976 | 456.7 | | | | 25 | W-113-2 | 2018 | O1-2y2 | 7204.56~7210 | 246.3 | | | | 26 | W-45-6 | 2018 | 设计O1-2y2 | —(正钻) | —(正钻) | | | |
|
Statistical table of lost circulation and emptying and reservoir modification during drilling with VSP-driven data processing in Tazhong area
|
[1] |
杜金虎, 周新源, 李启明 , 等. 塔里木盆地碳酸盐岩大油气区特征与主控因素[J]. 石油勘探与开发, 2011,38(6):652-661.
|
[1] |
Du J H, Zhou X Y, Li Q M , et al. Characteristics and controlling factors of the large carbonate petroleum province in the Tarim Basin,NW China[J]. Petroleum Geophysical Exploration, 2011,38(6):652-661.
|
[2] |
王绍先 . 随钻地震技术的研究及应用[D]. 北京:中国科学院地质与地球物理研究所, 2008.
|
[2] |
Wang S X . Research and application of seismic technology while drilling[D].Beijing:Institute of Geology and Geophysics, Chinese Academy of Sciences, 2008.
|
[3] |
Celniker G W . Performing geoscience interpretation with simulated data[P].US, US6256603B1, 2001.
|
[4] |
王成礼, 叶古清 . 利用非零井源距VSP资料计算层速度[J]. 石油地球物理勘探, 1988,23(5):616-619.
|
[4] |
Wang C L, Ye G Q . Calculating formation velocity using VSP data of non-zero well spacing[J]. Petroleum Geophysical Exploration, 1988,23(5):616-619.
|
[5] |
朱光明, 高静怀 . 用VSP资料估算Q值[C]//中国地球物理学会学术年会, 1991.
|
[5] |
Zhu G M, Gao J H . Estimating Q value with VSP data[C]//Annual Meeting of the Chinese Geophysical Society, 1991.
|
[6] |
张大伟, 孙赞东, 王学军 , 等. 利用零井源距VSP资料进行品质因子反演[J]. 石油地球物理勘探, 2011,46(s1):47-52.
|
[6] |
Zhang D W, Sun Z D, Wang X J , et al. Quality factor inversion using zero well spacing VSP data[J]. Petroleum Geophysical Exploration, 2011,46(s1):47-52.
|
[7] |
王正和, 崔永福, 向东 . 井控处理中的真振幅恢复与Q补偿方法及应用[J]. 物探与化探, 2008,32(4):434-437.
|
[7] |
Wang Z H, Cui Y F, Xiang D . The method and application of true amplitute recovery and Q-compensation in well-control processing[J]. Geophysical and Geochemical Exploration, 2008,32(4):434-437.
|
[8] |
Bakulin A, Slater C, Bunain H, 等. 在有横向非均匀性情况下由多方位变井源距VSP估算方位各向异性和裂缝参数[C]//SEG年会论文概要, 2000.
|
[8] |
Bakulin A, Slater C, Bunain H , et al. In the case of transverse heterogeneity,estimating azimuthal anisotropy and fracture parameters by multi-azimuthal variable well spacing VSP[C]//SEG Annual Conference Paper Summary, 2000.
|
[9] |
侯爱源, 耿伟峰, 张文波 . Walkaway VSP井旁各向异性速度模型反演方法[J]. 石油地球物理勘探, 2016,51(2):301-305.
|
[9] |
Hou A Y, Geng W F, Zhang W B . Inversion method of anisotropic velocity model near wellbore by Walkaway VSP[J]. Petroleum Geophysical Exploration, 2016,51(2):301-305.
|
[10] |
黄翼坚 . 多井源距VSP速度分析及逆时偏移[D]. 西安:长安大学, 2010.
|
[10] |
Huang Y J . Velocity analysis and reverse-time migration of multi-well source spacing VSP[D]. Xi’an:Chang’an University, 2010.
|
[11] |
罗小明, 凌云, 牛滨华 , 等. 井地联合提取VTI介质各向异性参数[J]. 地学前缘, 2005,12(4):576-580.
|
[11] |
Luo X M, Ling Y, Niu B H , et al. Joint extraction of anisotropic parameters of VTI media[J]. Geoscience Frontier, 2005,12(4):576-580.
|
[12] |
Liu Y M, Liang X H, Yin X Y , et al. Estimation of Thomsen’s anisotropic parameters from walkaway VSP and applications[J]. Applied Geophysics, 2014,11(1):23-30.
|
[13] |
谭仁春, 杜清运, 杨品福 , 等. 地形建模中不规则三角网构建的优化算法研究[J]. 武汉大学学报信息科学版, 2006,31(5):436-439.
|
[13] |
Tan R C, Du Q Y, Yang P F , et al. Study on the optimization algorithm of irregular triangular network construction in terrain modeling[J]. Journal of Wuhan University Information Science Edition, 2006,31(5):436-439.
|
[14] |
花卫华, 邓伟萍, 刘修国 , 等. 一种改进的不规则三角网格曲面切割算法[J]. 地球科学, 2006,31(5):619-623.
|
[14] |
Hua W H, Deng W P, Liu X G , et al. An improved algorithm for irregular triangular mesh surface cutting[J]. Geoscience, 2006,31(5):619-623.
|
[15] |
朱庆, 陈楚江 . 不规则三角网的快速建立及其动态更新[J]. 武汉测绘科技大学学报, 1998,23(3):204-207.
|
[15] |
Zhu Q, Chen C J . Rapid establishment and dynamic updating of irregular triangular networks[J]. Journal of Wuhan University of Surveying and Mapping Science and Technology, 1998,23(3):204-207.
|
[16] |
吴宇晓, 张登荣 . 生成三角网的快速合成算法[J]. 浙江大学学报理学版, 2004,31(3):343-348.
|
[16] |
Wu Y X, Zhang D R . Fast synjournal algorithm for generating triangular networks[J]. Journal of Zhejiang University, Science Edition, 2004,31(3):343-348.
|
[17] |
郭良辉, 孟小红, 郭志宏 , 等. 地球物理不规则分布数据的空间网格化法[J]. 物探与化探, 2005,29(5):438-442.
|
[17] |
Guo L H, Meng X H, Guo Z H , et al. Gridding methods of geophysical irregular data in space domain[J]. Geophysical and Geochemical Exploration, 2005,29(5):438-442.
|
[18] |
彭更新, 但光箭, 郑多明 , 等. 塔里木盆地哈拉哈塘地区三维叠前深度偏移与储层定量雕刻[J]. 中国石油勘探, 2011,16(5):49-56.
|
[18] |
Peng G X, Dan G J, Zheng D M , et al. 3D Prestack depth migration and quantitative characterization of carbonate reservoirs in Halahatang Area,Tarim Basin[J]. China Petroleum Exploration, 2011,16(5):49-56.
|
[19] |
黄元溢, 罗仁泽, 王进海 , 等. 几种叠前深度偏移技术效果的对比[J]. 物探与化探, 2011,35(6):798-803.
|
[19] |
Huang Y Y, Luo R Z, Wang J H , et al. A comparison of the application effects of several prestack depth migration techniques[J]. Geophysical and Geochemical Exploration, 2011,35(6):798-803.
|
[20] |
常建华, 张玉娥, 林卫 , 等. 三维叠前深度偏移处理技术应用[J]. 物探与化探, 1999,23(4):271-276.
|
[20] |
Chang J H, Zhang Y E, Lin W , et al. Application of 3-D prestack depth migration processing technology[J]. Geophysical and Geochemical Exploration, 1999,23(4):271-276.
|
[21] |
胡英, 姚逢昌 . 用于叠前深度偏移速度建模的一种新方法[J]. 石油勘探与开发, 2000,27(2):62-64.
|
[21] |
Hu Y, Yao F C . A new method for prestack depth migration velocity modeling[J]. Petroleum Exploration and Development, 2000,27(2):62-64.
|
[22] |
陈生昌 . 波动方程叠前深度偏移成像与速度建模方法研究[D]. 上海:同济大学, 2002.
|
[22] |
Chen S C . Wave equation prestack depth migration imaging and velocity modeling method[D]. Shanghai:Tongji University, 2002.
|
[23] |
马彦彦, 李国发, 张星宇 , 等. 叠前深度偏移速度建模方法分析[J]. 石油地球物理勘探, 2014,49(4):687-693.
|
[23] |
Ma Y Y, Li G F, Zhang X Y , et al. Prestack depth migration velocity modeling method analysis[J]. Petroleum Geophysical Exploration, 2014,49(4):687-693.
|
[24] |
高书琴, 冯许魁, 郭向宇 , 等. 波动方程叠前深度偏移处理技术开发与应用[J]. 天然气工业, 2007,27(s1):205-207.
|
[24] |
Gao S Q, Feng X K, Guo X Y , et al. Development and application of wave equation prestack depth migration processing technology[J]. Natural Gas Industry, 2007,27(s1):205-207.
|
[25] |
谢会文, 许永忠, 郑多明 , 等. 火成岩速度反演及三维速度场建立的关键技术及效果分析[J]. 物探与化探, 2013,37(6):1071-1079.
|
[25] |
Xie H W, Xu Y Z, Zheng D M , et al. An analysis of the key techniques and effects of the igneous rock velocity inversion and the establishment of 3D velocity field[J]. Geophysical and Geochemical Exploration, 2013,37(6):1071-1079.
|
[1] |
CUI Zhi-Qiang, XU Zhi-Li, LI Fei. The application of high-precision aeromagnetic survey to the investigation of oil-gas geological structure in Southwest Tarim Basin[J]. Geophysical and Geochemical Exploration, 2021, 45(4): 846-858. |
[2] |
Li-Jun GAO, Zong-Jie LI, Hai-Ying LI, Hong WANG, Chao HUANG. The deep karst fissure and cavern reservoir in Tarim basin carved in three dimensions:Research and application of "five-step method" quantitative description technology[J]. Geophysical and Geochemical Exploration, 2020, 44(3): 691-697. |
|
|
|
|